Stationary flow fields prediction of variable physical domain based on proper orthogonal decomposition and kriging surrogate model

In this paper a new flow field prediction method which is independent of the governing equations, is developed to predict stationary flow fields of variable physical domain. Predicted flow fields come from linear superposition of selected basis modes generated by proper orthogonal decomposition(POD). Instead of traditional projection methods, kriging surrogate model is used to calculate the superposition coefficients through building approximate function relationships between profile geometry parameters of physical domain and these coefficients. In this context,the problem which troubles the traditional POD-projection method due to viscosity and compressibility has been avoided in the whole process. Moreover, there are no constraints for the inner product form, so two forms of simple ones are applied to improving computational efficiency and cope with variable physical domain problem. An iterative algorithm is developed to determine how many basis modes ranking front should be used in the prediction. Testing results prove the feasibility of this new method for subsonic flow field, but also prove that it is not proper for transonic flow field because of the poor predicted shock waves.     

Assessment study of ionosphere correction model using single- and multi-shell algorithms approach over sub-Saharan African region

It is a known fact that ionosphere is the largest and the least predictable among the sources of error limiting the reliability and accuracy of Global Navigation Satellite Systems (GNSS) and its regional augmentation systems like Satellite Based Augmentation System (SBAS) in a safety-of-life application. The situation becomes worse in the Equatorial Ionization Anomaly (EIA) region, where the daytime ionization distribution is modified by the fountain effect that develops a crest of electron density at around ±15° to ±20° of the magnetic equator and a trough at the magnetic equator during the local noon hours. Related to this phenomenon is the appearance of ionosphere irregularities and plasma bubbles after local sunset. These may degrade further the quality of service obtained from the GNSS/SBAS system of the said periods. Considering the present operational augmentation systems, the accuracy and integrity of the ionosphere corrections estimate decreases as the level of disturbances increases. In order to provide a correct ionosphere correction to the user of GNSS operating in African EIA region and meet the integrity requirements, a certified ionosphere correction model that accurately characterizes EIA gradient with the full capacity to over-bound the residual error will be needed. An irregularities detector and a decorrelation adaptor are essential in an algorithm usable for African sub-Saharan SBAS operation. The algorithm should be able to cater to the equatorial plasma vertical drifts, diurnal and seasonal variability of the ionosphere electron density and also should take into account the large spatial and temporal gradients in the region. This study presents the assessment of the ionosphere threat model with single and multi-layer algorithm, using modified planar fit and Kriging approaches.     

基于 Kriging与多点取样算法的压气机特性建模

针对航空燃气轮机流量系数插值建模过程中,获取样本数据时的实验成本较高以及插值模型精度偏低这 2个问题,根据样本点在设计空间中的分布与模型精度的关系,在 Kriging插值模型的基础上,采用把垂距作为设计变量取值标准的取样算法。首先,将样本点与相邻 2点之间连线的垂距与垂距阈值进行比较,筛选出符合条件的数据点,形成候选点集;其次,选用高斯函数判定基础点集和候选点集的相关性;最后,使用该算法筛选后的样本点构建 Kriging插值模型。实验结果表明,该方法使用的样本点为原来数量的 70%左右,并且预测精度得到了提高,该方法是有效可行的。     

基于代理模型的航天器振动夹具优化方法

航天器振动夹具的一般设计原则是在避免安装干涉的前提下,尽量提高其刚度,同时减小自身重量对试验条件制定的影响。文章从提高振动试验夹具设计的效率出发,在优化过程中引入Kriging代理模型来代替原有的动力学输入输出关系,最大程度地减少优化迭代过程中有限元模型重构与求解次数。为进一步提高代理模型更新的效率,提出了一种混沌邻域搜索的多点加点准则,在发掘当前最优点的同时,有选择性地加入了更多的预测点,用以快速改善代理模型的局部代理精度。数值仿真针对典型航天器力学环境试验夹具的几何参量进行优化,优化结果验证了该方法的有效性与高效性。     

2009年1月平流层爆发性增温期间全球电离层响应的研究

2009年1月平流层爆发性增温(Stratospheric Sudden Warming, SSW)事件是有记录以来最强、持续时间最长的一次主增温事件(Major Warming Event, MWE), 期间太阳活动和地磁活动均处于较低的水平, 因此非常有利于研究电离层对平流层增温事件的响应情况. 本文利用COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate)系统提供的掩星数据, 使用Kriging方法分别构建了此次SSW期间及平静期的电离层N_mF₂, h_mF₂和110～750km高度范围的垂直积分TEC (简称VTEC)地图, 同时从全球定位导航卫星系统服务组织(International GNSS Service, IGS)发布的全球电离层TEC地图(Global Ionospheric Maps, GIMs)中提取了日固坐标系(Sun-fixed, 地磁纬度和地方时)下全球TEC地图. 通过对比发现, SSW期间与平静期相比, 地磁纬度中低纬电离层参数存在早晨上升, 下午和夜间下降的现象. 利用OSTM/JASON-2卫星高度计观测值进行验证后的结果显示, 此前研究均未有提及的夜间时段电离层参数N_mF₂, h_mF₂和TEC (VTEC和IGS TEC)的下降现象的确存在.      

基于统计学方法的卫星在轨热变形快速分析

为使用在轨实测温度数据辅助完成成像载荷视轴指向的在轨标定,首次将应用于地统计学中的普通克里格法用于卫星在轨温度场的无偏估计,提出了一种基于统计学方法的在轨热变形快速分析方法。应用热试验过程中星上少量温度实测数据进行卫星温度场的无偏估计,并将温度计算结果赋值到有限元模型,进而完成了热变形分析。应用此种方法,相比通过热物理方法进行温度场分析,可将热变形分析的工期缩短数日,且载荷变形的分析结果与试验实测结果的误差在1′左右,可较准确地得到载荷视轴在热变形作用下的指向变化。文章可为卫星在轨热变形的快速分析提供参考。     

Kriging模型在机翼气动外形优化中的应用

针对粒子群等随机优化算法计算量大的缺点,发展了基于Kriging模型的优化方法.采用改进的量子粒子群算法对Kriging模型的相关模型参数进行优化,以提高代理模型预测精度,并与具有双层结构的粒子群算法相结合.采用雷诺平均N-S方程流场求解器与多目标非线性适应值加权方法,对高维度多目标多约束的跨声速机翼进行了优化,设计的机翼具有理想的压力分布,降低了机翼阻力系数,并且有效控制了低头力矩和翼根弯矩,表明该方法具有较强的工程实用性.     

基于单循环方法的涡轮叶片可靠性及多学科设计优化

为了得到适用于涡轮叶片复杂结构并同时考虑可靠性的多学科设计优化方法,将基于单循环方法的可靠性分析（SLBRA）与并行子空间设计优化方法 （CSSO）相结合,提出了一种基于可靠性的多学科设计优化（RBMDO）方法。在优化过程中使用Kriging近似模型并不断提高模型精度。该方法在计算最可能失效点（MPP）的过程中避免了优化迭代,提高了计算效率。以涡轮叶片的设计优化为例,对该方法进行了验证并与传统双循环方法进行了对比。结果表明,优化结果满足可靠性的要求,与双循环方法相比优化效率明显提高,证明了该方法在工程应用中的可行性和有效性。     

民用飞机航材支援网络选址因素权衡方法研究

民用飞机航材支援网络作为主制造商为客户提供航材的必要媒介,合理地搭建该航材支援网络有利于提高主制造商的服务品质和客户运营能力。基于Kriging方法,提出一种定性与定量相结合的民用飞机航材支援网络选址因素权衡方法,对选址因素进行梳理和分类,获取选址定量分析数据需求,并运用所提方法进行民用飞机航材支援网络选址因素权衡,进而对选址因素进行重要度评价;将该方法与层次分析法、响应面法、蒙特卡洛法进行对比分析,以验证其有效性和适用性。结果表明：本文所提方法可有效进行选址因素重要度评价,为民用飞机航材支援网络的搭建奠定基础。     

Temporal extrapolation of TEC using WNN during 2007–2018 and comparison against GIM,IRI2016 and Kriging

In this paper, a new method of temporal extrapolation of the ionosphere total electron content (TEC) is proposed. Using 3-layer wavelet neural networks (WNNs) and particle swarm optimization (PSO) training algorithm, TEC time series are modeled. The TEC temporal variations for next times are extrapolated with the help of training model. To evaluate the proposed model, observations of Tehran GNSS station (35.69°N, 51.33°E) from 2007 to 2018 are used. The efficiency of the proposed model has been evaluated in both low and high solar activity periods. All observations of the 2015 and 2018 have been removed from the training step to test the proposed model. On the other hand, observations of these 2 years are not used in network training. According to the F10.7, the 2015 has high solar activity and the 2018 has quiet conditions. The results of the proposed model are compared with the global ionosphere maps (GIMs) as a traditional ionosphere model, international reference ionosphere 2016 (IRI2016), Kriging and artificial neural network (ANN) models. The root mean square error (RMSE), bias, dVTEC = |VTEC_GPS ? VTEC_Model| and correlation coefficient are used to assess the accuracy of the proposed method. Also, for more accurate evaluation, a single-frequency precise point positioning (PPP) approach is used. According to the results of 2015, the maximum values of the RMSE for the WNN, ANN, Kriging, GIM and IRI2016 models are 5.49, 6.02, 6.34, 6.19 and 13.60 TECU, respectively. Also, the maximum values of the RMSE at 2018 for the WNN, ANN, Kriging, GIM and IRI2016 models are 2.47, 2.49, 2.50, 4.36 and 6.01 TECU, respectively. Comparing the results of the bias and correlation coefficient shows the higher accuracy of the proposed model in quiet and severe solar activity periods. The PPP analysis with the WNN model also shows an improvement of 1 to 12 mm in coordinate components. The results of the analyzes of this paper show that the WNN is a reliable, accurate and fast model for predicting the behavior of the ionosphere in different solar conditions.