用双流体模型模拟湍流两相流场

首先针对颗粒浓度在两相流动中的重要性，建立了颗粒体积浓度方程，阐述了颗粒浓度在求解过程中的特点及颗粒体积浓度方程在计算方法中的贯穿。其次计算了沙层上的稀相可压气固两相流、直管内的不同压液固两相流及叶片扩压器中的密相可压力固两相流的流动，发现了颗粒浓度两种截然不同的分布，并验证了用不可压或稀相的模型来计算可压或密相的模型会带来误差，说明了在计算可压密相两相流动时，可压及密相的因素必须考虑到。     

二维凹槽超声速湍流流动数值模拟

本文采用时间、空间均为二阶精度的NND差分格式,以及Menterk ωSST二方程湍流模型,数值模拟了二维开式凹槽超声速粘性流动。用组分扩散方程计算了凹槽内粒子驻留时间。本文给出了M∞=3,长深比L/D=3,后壁倾斜角为90°和30°两种几何形状凹槽的计算结果,与相应的实验符合甚好。计算亦表明,该模型较高地估计了涡粘性,从而使流动更稳定。     

细长体大迎角湍流流场的数值模拟

通过数值方法对大迎角细长体超声速流场的模拟，探求正确模拟导弹大迎角绕流湍流流场的简单而有效的方法。数值模拟求解一般曲线坐标系下的三维可压缩Navier-Stokes方程，时间离散采用Euler向后差分，无粘项的空间离散采用二阶TVD格式，分别研究了Bald-win-Lomax代数模型及其修正形式(BLDS)对大迎角分离流动的模拟能力。数值试验表明修正的Baldwin-Lomax模型更精确地预测了流场的旋涡与分离情况，给出合理的表面压力分布，因而更准确地模拟了存在横向分离的导弹流场。     

铝合金反重力铸造充型过程数值模拟

介绍了反重力铸造(CGC)技术,分析了CGC条件下合金液充填铸型过程数值模拟的数学模型,利用建立的数学模型进行了反重力铸件充填过程模拟,并进行叶轮铸件研制和模拟结果的对比分析.     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

矢量喷管壁温分布的数值计算研究

计算了矢量喷管的对流和辐射换热及壁温的二维分布,考虑了对各热部件采用相应的冷却方式,比较了不同的气膜冷却模型对喷管扩张调节片的冷却效果,计算喷管的辐射换热时分光谱进行,考虑喷管内燃气参数的三维分布,编制了对轴对称矢量喷管、圆转方二元矢量喷管、球面过渡段矢量喷管在加力与非加力状态下均适用的计算程序,并给出了几个算例的计算结果。结果表明:偏转矢量喷管的壁温沿周向有明显变化,沿偏转方向的壁温较低,相反方向的壁温较高。      

双侧进气突扩燃烧室流场数值模拟

采用有限体积法和SIMPLE算法对冲压发动机突扩燃烧室流场进行了数值模拟,湍流模型采用标准二方程湍流模型,燃烧模型采用EDM模型。给出了流场中速度、速度流函数、湍流强度以及温度的分布。分析了流场特征,并与无燃烧状态下的流场进行了比较。     

直升机旋翼/机体耦合非线性系统的动稳定性分析

为了准确反映直升机旋翼／机体耦合系统的动稳定性，建立了旋翼／机体耦合非线性动力学微分方程，在时域内求解微分方程得到各片桨叶的挥舞、摆振及机体的响应用以对系统进行数值模拟；为了获得系统稳定性的定量值，用快速傅立叶变换(FFT)确定模态频率，用基于傅立叶级数的移动矩形窗方法得到模态阻尼。地面共振分析表明，时域分析与特征值分析结果具有良好的相关性，并与试验值吻合，从而验证了该方法的有效性。大总距时，用时域分析得到的模态阻尼与试验值吻合得更好，该方法可用于具有非线性减摆器的直升机旋翼／机体耦合系统的动稳定性分析。     

Gravity Field Recovery from GRACE: Unique Aspects of the High Precision Inter-Satellite Data and Analysis Methods

The very high accuracy of the Doppler and range measurements between the two low-flying and co-orbiting spacecraft of the GRACE mission, which will be at the μm/sec and ≈10 μm levels respectively, requires that special procedures be applied in the processing of these data. Parts of the existing orbit determination and gravity field parameters retrieval methods and software must be modified in order to fully benefit from the capabilities of this mission. This is being done in the following areas: (i) numerical integration of the equations of motion (summed form, accuracy of the predictor-corrector loop, Encke's formulation): (ii) special inter-satellite dynamical parameterization for very short arcs; (iii) accurate solution of large least-squares problems (normal equations vs. orthogonal decomposition of observation equations); (iv) handling the observation equations with high accuracy. Theoretical concepts and first tests of some of the newly implemented algorithms are presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

卫星推进系统发动机启动过程数值仿真

以一维管道瞬变流理论和特征线方法为基础,建立了仿真系统部件流体动力学模型。计算方法采用带内插的特征线方法和四阶龙格库塔法。对燃烧室模型,考虑燃烧时滞的影响,计算了卫星推进系统发动机在启动过程中各参数的变化。仿真结果与卫星推进系统热试车试验结果基本吻合。结果表明计算模型较好地描述了发动机启动过程中的水击和管流振荡现象,采用的仿真方法符合精度要求,在工程上可为卫星推进系统的设计与试验提供指导。