k-ε-g紊流燃烧模型数值研究

本文对k-ε-g紊流燃烧模型中模型系数C_R,C_(g1),C_(g2)的选定和计算浓度脉动均方值g的不同方法进行数值研究。计算结果表明模型系数值的选定直接与化学反应速率有关,而其中C_R影响最大。g的计算虽然可有不同方法,而其中微分方程求解最佳。并用实验数据来验证计算结果。     

三维侧压高超声速进气道不启动流场试验与数值模拟研究

对某构型三维侧压高超声速进气道开展了Ma4的自由射流试验和数值仿真,研究 了低马赫数下不启动流场的流动机理。观测到了具有“对涡”结构的底板油流图案,并得到 了分离区的范围和内部流动特征。分析得出,“对涡”结构油流图案和流场中部分离区的形 成是唇口激波和内收缩形成的逆压梯度作用于侧板激波形成的流场中部低能流区的结果。根 据试验和数值模拟结果给出了流场结构示意图,并为下一步进气道构形设计和性能改善工作 提出了建议。     

液氮贮箱自增压实验与数值仿真

为研究低温推进剂在常温下的自增压过程,设计了以液氮为模拟介质可视化低温玻璃贮箱自增压实验系统,研究了自增压过程压力和温度的变化规律及体积充填率对压力和温度变化的影响。实验结果表明：气枕区和液体区存在显著的轴向温度分层,液体区温度的上升速率低于压力引起饱和温度的上升速率。压力上升分为有典型意义的三段：初始段、过渡段和稳定段,稳定段的压力上升速率随体积充填率增加而增加。液体区的对流运动在自增压过程受到抑制,气液界面逐渐进入准静止状态。并以实验测得温度作为边界条件,采用流体体积（VOF）模型对整个自增压过程进行了175 s的数值仿真。仿真得到的压力曲线变化规律与实验结果基本一致,稳定段的压力上升速率是实验值的1.58倍。本文得到的自增压物理参数变化规律,为低温推进剂的贮存和贮箱的热防护设计提供参考。     

载人航天器生活舱内湿度场的稳态数值模拟

研究载人航天器生活舱内部的湿度情况具有重要意义。由于生活舱设备外形和结构相当复杂,难以建模,因此按照速度梯度、湿度梯度较大的优先原则,合理选择简化舱体结构及舱内设备。采用基于有限元网格的控制容积法对计算区域离散,对离散方程组求解得出水汽质量浓度分布后,结合舱内温度场通过一定的转换关系最终得到了生活舱内部的湿度分布。数值分析结果表明生活舱内处于相对湿度较高的状态下,提示还需进一步采取措施以改善舱内湿度水平。     

喷雾运动初始条件问题的研究

根据所完成的流动环境中颗粒在初始区域内速度分布的测量结果，对颗粒运动初始条件的问题进行了全面的研究，提出了统计性，随机性和确定性的三种初始条件及其数学处理方法，在此基础上，应用喷雾两相流的数值方法，计算了喷雾的运动，比较了不同初始条件对计算结果的影响，并将所得计算结果与实验结果做了比较，两者符合较好。     

飞机数字化装配系统中机身交点孔加工方法研究

针对传统交点孔加工方法加工精度低且存在回弹变形的问题,设计了一种满足飞机大部件数字化装配要求的机身交点孔精加工方法,根据交点孔的加工特性完成了机身定位与固持、数控加工中心定位以及加工工艺的总体设计,确定了交点孔加工程序原点的计算方法和编程原则,最后通过交点孔试切加工试验验证了加工方法的可行性。     

A CFD-based numerical virtual flight simulator and its application in control law design of a maneuverable missile model

A CFD-based Numerical Virtual Flight (NVF) simulator is presented, which integrates an unsteady flow solver on moving hybrid grids, a Rigid-Body Dynamics (RBD) solver and a module of the Flight Control System (FCS). A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing, with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces (e.g. the afterbody rudders of a missile). For the flow/kinematic coupled problems, the 6 Degree-Of-Freedom (DOF) equations are solved by an explicit or implicit method coupled with the URANS CFD solver. The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model. A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model. Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law. Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.     

A new global version of M(3000)F2 prediction model based on artificial neural networks

A new version of global empirical model for the ionospheric propagation factor, M(3000)F2 prediction is presented. Artificial neural network (ANN) technique was employed by considering the relevant geophysical input parameters which are known to influence the M(3000)F2 parameter. This new version is an update to the previous neural network based M(3000)F2 global model developed by Oyeyemi et al. (2007), and aims to address the inadequacy of the International Reference Ionosphere (IRI) M(3000)F2 model (the International Radio Consultative Committee (CCIR) M(3000)F2 model). The M(3000)F2 has been found to be relatively inaccurate in representing the diurnal structure of the low latitude region and the equatorial ionosphere. In particular, the existing hmF2 IRI model is unable to reproduce the sharp post-sunset drop in M(3000)F2 values, which correspond to a sharp post-sunset peak in the peak height of the F2 layer, hmF2. Data from 80 ionospheric stations globally, including a good number of stations in the low latitude region were considered for this work. M(3000)F2 hourly values from 1987 to 2008, spanning all periods of low and high solar activity were used for model development and verification process. The ability of the new model to predict the M(3000)F2 parameter especially in the low latitude and equatorial regions, which is known to be problematic for the existing IRI model is demonstrated.     

预旋进气小尺寸涡轮叶片冷却的流场研究

为了了解和掌握一种具有直通式冷气预旋进气系统的小型燃气轮机涡轮转子叶片的流场,在旋转雷诺数Reθ=4.66×106和冷却空气的无量纲质量流量Cw=1750时改变预旋角θ的大小,使其在15°～90°变化,通过数值研究得到了预旋角对涡轮盘腔、连管和涡轮叶片内冷却空气的流动以及叶栅通道中燃气的流动的影响。结果表明:(1)预旋角的变化会改变涡轮盘腔、连管和涡轮叶片冷气进口附近局部区域的流场,但是对涡轮叶片内其它区域和叶栅通道中的流动基本没有影响。(2)随着预旋角的增大,涡轮盘腔内预旋进气冷气射流的轴向穿透深度先增大后减小;当θ<45°时冷却空气沿外围屏流向转盘接收孔,而当θ>45°时冷却空气沿内围屏流向转盘接收孔;气流的周向速度随着预旋角的增大而减小。(3)垂直进气时连管内存在多个回流区和很大的涡流,流动损失较大,而采用预旋进气能够减弱或消除这些流动结构,存在最优预旋角θopt,θopt≈45°,此时连管的有效流通面积最大。     

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验