混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

三维复杂外形横截面拓扑结构随迎角变化研究

通过数值求解抛物化的椭圆型方程生成复杂外形三维网格;采用交替方向隐式分解的隐式NND格式求解全N S方程模拟"类升力体"外形在高超声速下的大迎角流动;采用张涵信发展的流动拓扑结构理论分析给出了"类升力体"外形垂直于体轴的横截面流线随迎角变化的拓扑结构;此外,通过计算发现:大于20°迎角后,在部分横截面背风对称线上出现结构不稳定的鞍点相连现象。     

分区对接网格算法的应用研究

本文应用分区对接网格技术，由于采用了对接面上网格点完全对应的方法，能自动满足通量守恒和计算的稳定性。采用有限体积算法下的原始变量型NND格式，求解NS方程，数值模拟了复杂外形的高超声速层流流场，计算结果与实验测量符合较好。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

两种空气/煤油燃气发生器富油燃烧组织方案对比

根据空气/煤油富油燃烧的特点,提出了两种空气/煤油燃气发生器富油燃烧组织方案,设计了采用钝体稳定火焰和旋流空气、二次喷注空气稳定火焰的两种燃气发生器.为了对比两种方案的点火和燃烧特性,对两种燃气发生器进行了一系列热试,结果表明余氧系数是燃气发生器最重要的工况参数.随着余氧系数的增加,燃气发生器的状态逐渐从启动失败变为中途熄火,最终呈正常启动状态.采用钝体稳定火焰的燃气发生器稳定工作的余氧系数边界为0.518,采用旋流空气和二次喷注空气稳定火焰可将该边界延伸到0.237,极大地扩大了燃气发生器的工作范围.与钝体稳定火焰的燃气发生器相比,旋流空气和二次喷注空气稳定火焰的燃气发生器的富油燃烧的燃烧效率提高了20%.两种方案结构复杂性相当,旋流空气和二次喷注空气稳定火焰的燃气发生器不需要冷却火焰稳定器,可提高燃气发生器的工作时间.      

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

Adaptive Gaussian sum squared-root cubature Kalman filter with split-merge scheme for state estimation

The paper deals with state estimation problem of nonlinear non-Gaussian discrete dynamic systems for improvement of accuracy and consistency. An efficient new algorithm called the adaptive Gaussian-sum square-root cubature Kalman filter（AGSSCKF） with a split-merge scheme is proposed. It is developed based on the squared-root extension of newly introduced cubature Kalman filter（SCKF） and is built within a Gaussian-sum framework. Based on the condition that the probability density functions of process noises and initial state are denoted by a Gaussian sum using optimization method, a bank of SCKF are used as the sub-filters to estimate state of system with the corresponding weights respectively, which is adaptively updated. The new algorithm consists of an adaptive splitting and merging procedure according to a proposed split-decision model based on the nonlinearity degree of measurement. The results of two simulation scenarios（one-dimensional state estimation and bearings-only tracking） show that the proposed filter demonstrates comparable performance to the particle filter with significantly reduced computational cost.     

基于地月L1点的载人登月飞行方案分析

为解决传统载人登月方案中重型运载火箭难以研制、交会对接窗口选择困难等问题,提出了符合我国现有技术条件的基于地月L1点的载人登月方案设想。基于对圆形限制性三体模型和地月三,点特性的分析,将载人登月任务划分为不同的飞行阶段,研究了载人登月任务的规模、系统组成以及各舱段的质量分配等。通过与基于环月轨道方案各参数的分析和对比,提出通过付出较小代价将基于环月轨道的载人登月方案转化为基于地月,L1点的载人登月方案的方法。结果分析表明,经方案转化,不但能充分发挥地月L1点的优势,而且在一定条件下可以使节省燃料成为可能。