仅测角定轨问题迭代格式的推导与应用

针对由仅测角信息确定卫星轨道的问题，建模为求解相应适应度函数零点的一类方法，分析了求解该问题的高斯-牛顿迭代法中涉及的观测时间步长、动力学方程求解步长和迭代格式步长三者之间的区别与联系。针对一般精定轨方法由于适应度函数高度非线性而难以构造有效迭代格式的问题，提出了初定轨与精定轨方法结合的定轨方法，并在理论上分析了该方法的可行性。通过数值实验验证了所提求解方法的准确性、有效性和高效性。     

Roe格式中不同类型熵修正性能分析

针对Roe格式中较为流行的3类熵修正:Muller型、Harten-Yee型和Harten-Hyman型熵修正,从理论分析入手, 辅以Euler方程为控制方程的激波管问题,前台阶流动和运动激波的双马赫反射3个数值实验,对不同熵修正的性能做了深入研究,得出如下结论: Muller型和Harten-Yee型熵修正方法作用于激波和非物理的膨胀激波2种情况:激波情况下引入的数值耗散,能根本改善"Carbuncle"现象,膨胀激波时其数值粘性也足以使膨胀激波得以耗散;而类Harten-Hyman型熵修正只对膨胀过程起修正作用,对激波情况无效,不能改善"Carbuncle"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.     

背景太阳风数值模拟的熵守恒格式

背景太阳风研究是根据行星际扰动传播情况预测空间天气状况的基础,磁流体（MHD）模拟是背景太阳风研究的重要手段.采用一种新的数值计算方式,利用Ideal GLM-MHD将计算过程中产生的磁场散度以c_h的速度向计算区域外传播,从而消去磁场散度;重构部分使用受约束的最小二乘法,将磁场散度作为约束条件添加到重构中,进一步对重构后的磁场梯度进行修正;通量计算采用满足热力学第二定律的熵守恒格式,该格式能够确保在计算过程中熵不增,保证数值稳定.研究结果表明,该方法应用于太阳风数值模拟的求解得到了更加稳定的结果.      

Euler方程的分裂型通量分裂双时间步隐式方法

传统隐式方法有格式复杂、计算量大等缺点,在Euler方程的差分离散过程中,利用算子分裂思想,结合通量分裂法、双时间步法等隐式离散方法,构造了一种更简单的分裂型隐式计算方法.通过对典型空气动力学问题的计算,检验了该方法的有效性和可靠性,并对其性能做了具体讨论.该方法具有稳定性好、时间步长约束小等隐式格式的普遍优点,同时具有格式简单、程序易实现等优点;避免了传统隐式方法单步推进时的方程组常规求解及矩阵求逆过程,计算量小;比LU-SGS方法收敛速度快.      

基于改进多目标遗传算法的舰尾紊流模拟方法

为提高舰尾紊流自由大气紊流分量仿真的可信度,提出了一种紊流数值模拟的新方法.首先,使用带有修正系数的Euler前向差分格式表示各个方向的紊流序列,同时结合智能算法的思想,把紊流相关性检验中的均方差误差和相关函数误差作为待优化目标函数,将修正系数的选择看成一个多目标优化问题,并采用改进的多目标遗传算法进行求解.最后,通过仿真算例验证了本文方法的正确性与合理性,计算结果表明该方法可以根据不同的采样步长灵活地生成所需紊流.尤其在小步长情况下,亦可得到很好符合理论值的紊流序列,可以满足虚拟飞行实验的要求.      

稳定因子对FDTD数值计算的影响

讨论了FDTD(Finite-Difference Time-Domain)迭代计算中稳定因子的作用.推导了含有稳定因子的二阶Mur条件修正式,并以线源为例研究了稳定因子对吸收边界性能以及数值稳定性的影响.在三维计算空间中以软激励源形式加入平面波,计算了稳定因子不同时散射场区的RCS以及收敛性能的变化,结果表明在满足数值稳定条件的前提下,增加时间步长对计算更为有利.      

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

基于数据链的空战对抗建模及增援决策分析

研究了空战双方力量的变化问题.将飞机的空战能力指数转化为平均战斗力水平,进而以蓝彻斯特方程为基础,考虑数据链对空战的影响和作战中有增援力量的加入,建立了对抗双方战机数量变化的微分方程模型,给出了其数值计算方法,指出了离散计算的时间步长所应满足的条件,并根据增援模式的不同,将模型细化为不同形式.仿真实验表明,存在有效增援时间范围,一次性增援和匀减速增援具有时间优势,并且通过构造评价函数有效衡量了这两者具有较高的增援效率,为空中作战应用和决策提供了参考依据.      

气动谐振管加热机理数值研究

放置在高速喷流中的半开口气动谐振管在一定的条件下会出现加热现象.针对这一现象,利用数值计算方法计算了一维非定常非线性Euler方程组,采用MacCormark格式,对不同管口条件下谐振管内的振荡过程进行了数值研究.研究发现,当管口压力变化频率处在谐振管的线性谐振频率附近时,管内形成的激波以及激波的反射增强导致了谐振管端壁区域明显的加热效果.对不同的驱动气体和不同的谐振管型式(圆柱形和圆锥形)的数值研究也给出了有价值的结论.研究结果与实验事实或现象是相吻合的.      

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

Characterisation of residual ionospheric errors in bending angles using GNSS RO end-to-end simulations

Global Navigation Satellite System (GNSS) radio occultation (RO) is an innovative meteorological remote sensing technique for measuring atmospheric parameters such as refractivity, temperature, water vapour and pressure for the improvement of numerical weather prediction (NWP) and global climate monitoring (GCM). GNSS RO has many unique characteristics including global coverage, long-term stability of observations, as well as high accuracy and high vertical resolution of the derived atmospheric profiles. One of the main error sources in GNSS RO observations that significantly affect the accuracy of the derived atmospheric parameters in the stratosphere is the ionospheric error. In order to mitigate the effect of this error, the linear ionospheric correction approach for dual-frequency GNSS RO observations is commonly used. However, the residual ionospheric errors (RIEs) can be still significant, especially when large ionospheric disturbances occur and prevail such as during the periods of active space weather. In this study, the RIEs were investigated under different local time, propagation direction and solar activity conditions and their effects on RO bending angles are characterised using end-to-end simulations. A three-step simulation study was designed to investigate the characteristics of the RIEs through comparing the bending angles with and without the effects of the RIEs. This research forms an important step forward in improving the accuracy of the atmospheric profiles derived from the GNSS RO technique.     

正态云模型雾化性质统计分析

对于超熵较大情况下的正态云模型,说明了云模型雾化过程.通过统计分析云滴离散的整体趋势,说明超熵增大过程中,云滴整体趋于离散.通过分析各论域区间内云滴离散趋势,说明靠近概念核心的云滴的离散速度相对缓慢.归纳云模型雾化性质:在超熵取值持续增大的过程中(HeEn/3),正态云表示的概念的论域范围持续增大,呈雾化状态,但靠近概念核心的论域区间内的云滴不失数量优势.雾化性质适用于建模偏离正态分布、缺乏共识的定量数据,期望表示概念语义值核心,熵描述概念语义的离散程度,超熵表示各种语义的共识程度,扩展了云模型知识表示的应用范围.     

守恒律方程的导数人工压缩法

将高分辨率差分格式用于守恒律方程的导数方程,可以克服传统高分辨率格式在极值点精度退化的缺点,类似于Harten的人工压缩法,新方法称为导数人工压缩法.本方法既能提高间断的分辨率,又能提高极值点处的分辨率,是一种高分辨率低耗散低扩散格式.用单个守恒律方程带间断和多极值的初值问题和一维激波管问题进行了验证,比较了Harten-TVD,人工压缩,导数人工压缩方法的在间断和极值点的分辨率问题.     

高超音速双锥绕流数值模拟的格式效应分析

通过模拟双锥外形高超音速绕流,对各种计算流体力学空间离散格式的性能进行分析.研究了Roe格式、van Leer格式、AUSM(Advection Up-stream Splitting Method)格式、AUSMDV格式、AUSM+格式、上风TVD(Total Variation Diminishing)格式等对双锥高超音速流场特征的模拟能力.结果表明:各空间格式激波分辨率都较高,但粘性分辨率相差很大.其特性与格式耗散性有关,耗散性越小,粘性分辨率越高;格式耗散性对热流极值影响显著,耗散性越小,热流极值越大,其位置越靠近前缘.      

两种流体控制方案矢量喷管内流场计算及分析

在激波诱导和喉道倾斜2种流体控制方案下,对推力矢量喷管的二维内流场进行了数值模拟.用局部加密的四边形结构网格对流场进行网格划分之后,采用二阶迎风离散格式和二方程湍流模型求解强守恒形式的Navier-Stokes方程,通过对内流场的模拟获得流体控制推力矢量喷管方案的推力矢量特性.计算结果表明:对于激波诱导矢量控制方案,在小落压比大注气量且后缝注气的情况下所产生的矢量角最大;对于喉道倾斜矢量控制方案,在扩张片上以与主流成某一角度的方向注入气流产生的矢量角最大且没有带来较大的推力损失;喉道倾斜矢量控制方案因其推力损失较小且能够降低喷管的重量与造价而更具发展前景.      

绳驱空间机械臂动力学建模与仿真效率分析

针对一种适用于在轨服务精细操作的绳驱空间机械臂进行讨论。建立该机械臂的动力学模型，将机械臂视为链式多刚体系统，认为绳索满足线弹性假设并忽略绳索的质量，利用空间算子代数方法推导了机械臂的动力学方程。分析动力学模型的数值仿真效率，求解出机械臂微振动的振动频率，将最高阶振动频率作为特征系数代入到微分方程数值算法的稳定区域，计算得到仿真步长的临界值；对动力学模型中的绳索模型进行改进，通过增大仿真时间步长临界值来提高动力学模型的数值仿真效率。仿真算例表明，改进后动力学模型的仿真效率最多可提高200多倍，并可以实现实时仿真。