两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

高过载条件下绝热层烧蚀实验方法研究 (Ⅰ)方案论证及数值模拟

分析了高过载对固体火箭发动机流场和绝热层烧蚀的影响规律，提出一种新的研究思路：采用数值模拟方法来预示发动机三维两相流场，建立高过载流场模拟实验装置，开展绝热层烧蚀实验，建立高过载条件下的绝热层烧蚀模型，在此基础上发展高过载发动机绝热层设计和烧蚀预示方法。其中关键技术是高过载流场的模拟，对粒子加入法和弯管分离法两种方案进行了论证，认为弯管分离法原理上是可行的。为了验证这种方案的可行性，开展了弯管通道两相流的数值模拟研究，计算结果表明弯管装置具有使凝相粒子聚集形成高浓度粒子流的功能。     

光电瞄准设备通用计量平台

提出了光电瞄准设备通用计量平台的设计方案，论述了平台的硬件、软件构成，并着重分析了双分划双光源照准与两维同步动态测角、读数图像自动识别、多点位多通道实时检测的方法，最后对系统的不确定度进行了分析和评估。     

航空发动机燃烧室主燃区特性的研究

采用FLUENT6．1软件对某型机燃烧室主燃区的各个典型工作状态进行三维两相湍流燃烧的速度场、温度场、油滴轨迹进行了数值模拟计算。通过对计算结果的分析，得出了主燃区的特性：主燃区的冷热区决定了出口温场冷热区，其效率决定了燃烧室的效率；随着负荷增加，油滴轨迹明显呈缩短趋势，模拟计算结果与试验结果及燃烧室工作特点比较吻合。     

非壅塞火箭冲压发动机补燃室两相流数值模拟

采用颗粒轨道模型进行了非壅塞火箭冲压发动机补燃室两相流的数值模拟,其中铝颗粒的燃烧模型采用的是Davis扩散控制燃速的燃烧模型,建立了发动机补然室简单反应流模型,并在该模型下对某实验发动机进行了模拟,得出颗粒在补燃室内的分布,结果表明:与颗粒确定轨道模型相比,颗粒随机轨道模型更加适合模拟冲压发动机中的两相流动,并且颗粒的初始直径对颗粒燃烧效率有很大的影响。     

以目标机场为初始目的地备降场的二次放行

对二次放行飞行计划时能否以最终目标机场做初始目标机场的备降场等问题进行了讨论,通过理论分析和实际计算证明:以最终目标机场做初始目标机场的备降场的二次放行,一般不能增加业载,而且会使业载减少.当以最终目标机场做初始目标机场的备降场时,如果在二次放行点的剩余油量少于正常情况的再次放行所需油量,但满足做无备降场的二次放行的油量和气象条件的要求,可以由二次放行点再次放行到最终目标机场.     

影响Al2O3凝相尺寸分布的因素

用激光全息装置测量含铝复合推进剂燃烧场凝聚相的分布，分辨率为5μm。分析了滞留时间、铝粉含量（2％－16％）和燃烧室压强（0.1MPa-4MPa）对凝相粒子尺寸分布的影响，实验表明，粒子总数随滞留时间减少，粒子越大，减少的相对比例越大；粒子尺寸分布不随燃烧室压强和推进剂铝粉含量变化。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

航天飞机自动着陆轨迹优化设计

采用沿约束轨迹的推演计算方法及优化原理对航天飞机自动着陆飞行轨迹进行优化设计。首先根据过载、垂直下降率、连续性和平稳性要求确定飞行轨迹的基本形状，然后采用两点边值优化方法对轨迹参数进行优化设计以获得理想的接地动压。与传统的和制导系统联系在一起的轨迹设计方法相比，该方法具有简便、快速、可重复使用、鲁棒性强等特点。     

TC11合金高温变形行为及其机理

采用Gleeble-1500热模拟机研究了TC11合金在800～1 050℃、应变速率0.005～5/s条件下的高温变形行为.根据动力学分析,确定了不同温度区间的热激活能和热变形方程.结合变形微观组织观察确定了TC11合金的高温变形机制.结果显示:TC11合金在(α β)两相区和β相区的热变形激活能分别为285.38和141.98 kJ/mol,表明不同温度区间的热变形机理不同;在两相区变形主要发生片状组织的球化,在β相区变形时低应变速率下(0.005～0.05/s)主要发生β相的动态再结晶,高应变速率下(0.05～5/s)主要发生动态回复.研究结果为确定该合金的最佳变形工艺参数提供了理论依据.