超音速流在正规反射激波系后的参数优化特性

研究了超音速气流中正规反射波加一正激波结构下的参数优化特性，给出了求解这种优化结构的方法，通过计算得到了优化条件存在的区域和优化解的特性，它用可于指导超音速飞行器的进气道及相关装置的设计。     

斜切反喷管性能分析

固体火箭发动机前端斜切反喷管，其结构简单、作用时间短、气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。本文从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，结合采用Ｂｏｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，利用时间相关法及ＭａｃＯｃｏｒｍａｒｋ两步显格式求解，模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分布相一致。该方法可应用于斜切反喷管的     

导弹构形的初步设计与筛选程序

介绍非圆形弹身导弹的气动系数的计算程序。该程序速率快，通用性好，便于应用，因而特别适合导导构架的初步设计。为便于识别有开发前景的候选构形，该程序的精度符合分析多种构形和飞行条件的要求。     

统一分区加强解算法解物体穿越激波非定常流场

使用统一分区加强解算法求解非定常可压缩Navier-Stokes方程，模拟了物体分别以超音速和亚声速穿越与之同向运动激波的全过程，并且通过和单区计算结果的详细对比，说明统一分区加强解算法能够简便、有效的处理复杂几何外形和提高网格分辨率。     

有烧蚀气体产物影响的化学非平衡粘性流数值计算研究

本文是针对高超声速绕流钝锥体再入飞行器，研究有烧蚀气体产物影响的化学非平衡反应粘性流ＮＳ方程数值计算方法。通过带有化学源项的ＮＳ方程对流动作数值模拟，总的连续方程由单个组元守恒方程组代替。对１３个化学组元，２６个气相化学反应，以及２个气相－固相化学反应，研究烧烛碳－空气的化学反应系统，并与纯空气化学反应系统的计算结果进行比较。利用二阶迎风ＴＶＤ差分格式和激波捕捉法，计算球锥体头部和身部化学非平衡粘性流场     

冷气推进器流场仿真和分析

采用Fluent软件对火箭姿态控制系统小推力冷气推进器的内流场进行了数值仿真,得到了流场压力与马赫数的曲线图。气流在推进器控制口达到音速,通过控制口后加速并相互撞击产生激波,流束相对集中,超音速流经过激波变成亚音速流;亚音速流在直管段扩张达到超音速,在收缩段受压缩减速,流束集中后再次扩张加速,喷嘴喉部气流略大于音速,但喷嘴扩张段内没有激波,推进器工作状态正常。仿真结果与高模试验结果相吻合。     

冲压发动机进气道流场数值模拟

采用有限体积法对积分形式的Navier-Stokes方程组进行空间离散，所获得的常微分方程用多步Runge-Kutta方法求解。数值模拟了某导弹冲压发动机进气道的三维流场，给出了不同反压条件下该进气道的内外流场结构和表面压强分布，并分析了进气道内激波与附面层的相互作用现象。     

轴对称体外部流场的有旋特征线法

轴对称体的无粘，可压缩，定常，超音速流动的外部流场，在计算时，常采用流动为无旋的假设，当头部激波弯曲较大时，会产生一定的误差。本文对绕轴对称体的这类流动，用有旋特征线理论进行了数值计算。结果与实验数据的比较表明：用有旋特征线理论对绕轴对称体的流场进行计算是完全可行的。该方法的应用对轴对称体的气动力计算具有实际应用价值。     

轴对称横向喷流强干扰流场的巨型机向量并行计算

本文在国产巨型计算机ＹＨ－１和ＹＨ－２上对ＭａｃＣｏｒｍａｃｋ显格式解轴对称横向喷流强干扰流场的计算程序进行了重构和向量优化使其能并行执行。得到该程序在ＹＨ－１上的向量加速比约为６．０，ＹＨ－２单机向量运行速度约为ＹＨ－１上的向量运行速度的３倍。与传统的串行计算相比较，并行计算大大缩短了程序运行的ＣＰＵ时间。     

导弹气动特性工程计算通用程序设计与研究

以三级以下弹箭常见的气动外形为计算模型，编制了气动特性工程计算通用程序（ADM系统），该软件适用于计算马赫数0～10、飞行攻角0～30°范围内多种火箭和导弹的气动参数。利用风洞实验数据详细检验了计算误差，除跨音速段外，升力、阻力系数计算误差分别小于3％和3．5％，压心系数误差小于1％。经过防雹火箭和探空火箭使用验证，该程序用于火箭设计是可靠的。     

采用敏度分析和计算流体力学的气动力设计最优化

介绍一种以计算流体力学（ＣＦＤ）－敏度分析算法为基础的气动设计优化方法。这种方法新颖而有效，可用于设计轴向推力最佳的超音速冲压式喷气发动机后体外形。用无粘流分析方法求解欧拉方程，依次产生目标函数和各种约束。因此，ＣＦＤ分析是与最优化过程采用的约束极小化方法相结合的。对于最优化所需的敏度系数，即目标函数和约束的梯度采用半解析法获得，而不采用传统的有限差分平滑方法。在一维寻优过程中，为了降低计算成本，     

固体火箭发动机喷管瞬态流场特性分析

为准确预测某型固体火箭发动机喷管的流场特性,建立了发动机燃烧室-喷管一体化三维流场模型,考虑了上游流场-燃烧室对喷管流场的影响,应用有限体积法,仿真计算出了发动机点火启动过程中喷管内激波的存在及变化趋势,仿真结果与一维等熵函数理论分析结果相一致。结果表明,在发动机点燃初期,喷管内燃气呈现亚音速流动,随着时间推移,在喷管扩张段出现了一道激波,燃气流动出现壅塞,随着燃烧室内燃气压力升高,激波逐渐移出喷管,喷管内呈现超音速顺畅流动。     

基于流固耦合的部分进气涡轮数值模拟研究

发展了三维线性插值算法用于CSD/CFD耦合计算数据交换,对某型液体火箭发动机部分进气涡轮进行了气/热/固多学科耦合数值仿真。结果表明,发展的三维线性插值程序对网格类型限制性小,计算简单,计算量小,插值结果能够满足耦合计算要求。仿真结果表明,某型火箭发动机涡轮由于其部分进气结构设计和叶轮高速旋转,设计工况下在涡轮转子入口处产生了较强的激波,激波与边界层干涉不仅使涡轮转子叶片的载荷分布出现了强烈的不均匀性,同时在叶轮的高速旋转下,该涡轮转子受到强烈的气动、热交变力冲击,其结构强度问题变得尤为突出。耦合计算分析认为设计工况下,该型涡轮结构设计,转子强度能够满足要求。     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用，重点研究了波涡相互作用引起的激波结构的变化过程。应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Elder方程，对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究。漩涡与激波间相对强度的大小对激波结构的影响明显，弱相互作用诱导激波变形，强相互作用诱导激波变形并产生分叉。计算结果表明，应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息。     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用,重点研究了波涡相互作用引起的激波结构的变化过程.应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Euler方程,对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究.漩涡与激波间相对强度的大小对激波结构的影响明显,弱相互作用诱导激波变形,强相互作用诱导激波变形并产生分叉.计算结果表明,应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息.     

基于高阶WENO格式的喷管动态特性仿真分析

喷管的动态特性分析是进行高品质控制系统设计、控制参数优化、控制规律验证的基础。为了分析喷管的动态特性,针对轴对称收扩喷管,基于一维变截面非稳态无粘流动控制方程,利用有限体积Roe类型的高精度WENO格式对控制方程进行离散,采用三阶强稳定型龙格库塔方法进行时间方向的积分,利用激波探测函数实时捕获激波位置,建立了喷管的一维动态特性仿真程序。通过与理论解析解的对比验证了仿真程序和激波探测函数的正确性。在欠膨胀状态下,获得了喷管入口总压、总温扰动下,喷管进出口压力和推力的阶跃响应特性,对比分析了各扰动因素对推力响应特性的影响。进一步分析了扩张段存在激波的情况下,在不同的扰动下激波运动和推力等参数的变化规律,能够为推力控制系统分析提供参考。     

应用TTM网格研究冲压增程弹丸进气道内外流场

基于TTM贴体网格生成法,研究了某冲压增程弹丸超音速进气道内外流场数值模拟中单区域网格生成问题,流场求解采用了二阶隐式TVD格式。NACA0012翼型跨音速流场算例结果表明,所采用的算法和编制的程序可靠,能够适用于具有亚、跨、超音速三种流态共存的复杂流场数值计算。对超音速进气道数值模拟结果显示,得到的流场复杂波系结构是合理的,来流马赫数和攻角变化显著影响着进气道性能。     

导弹翼面结构一次二阶矩法可靠性设计研究

给出了导弹翼面结构可靠性预测与可靠性设计的随机有限元一次二阶矩法的运算公式；提出了翼面壁板厚度与加强肋宽度作可靠性设计的思路；编制了用于微机的计算程序。试验证明：本文的方法适用于常用导弹翼面结构的可靠性预测与可靠性设计     

卫星回收舱再入过渡流区气动热数值计算

返回式卫星进入地球大气层95km高度左右时,速度达到Ma20,回收舱被高超声速稀薄来流形成的弓形激波环绕,气动加热问题非常明显,准确预测过渡流区的气动热成为一个十分突出的问题。过渡区由于气体空气分子仍然比较密集,基于分子动力学的直接蒙特卡洛模拟方法(DSMC)极其耗费计算资源,而求解Navier-Stokes(N-S)方程的方法误差较大。采用添加二阶滑移条件的N-S方程求解过渡流区气动热,并与开源DSMC2V程序计算结果对比,研究了舱体母线热流、压强系数变化;通过对比两种方法下舱体前缘弓形激波及流场参数变化,分析了滑移条件影响壁面和流场参数的机理。结果表明:回收舱再入过渡区时,钝头驻点区为高压、高热流区,锥身区气动热和壁面压强保持在较低水平。带滑移条件的计算机流体力学方法(ComputationalFluidDynamics,CFD)计算得到的热流、压强系数与DSMC结果吻合良好,具有计算效率高、精度较高的优势。对流场压力、温度、速度等参数分析显示,滑移条件中壁面速度滑移和温度跳跃的加入,改变了壁面流动参数,进而改善了壁面热流和压强的准确模拟能力。滑移条件对外流场参数的影响极小,没有改善稀薄流区流场的模拟能力和激波捕获能力,模拟得到的温度、速度等流场参数及激波位置、激波层厚度等与DSMC结果仍有差别。可以认为,文章采用方法能够满足工程上快速高效预测卫星回收舱再入过渡流区气动热需要。     

惯性组合-计算机数学模型辨识方法研究

根据频率响应测试的结果，本文讨论了离散传递函数最小二乘解辨识方法，并对某战术导弹捷联惯导控制系统的惯性组合－计算机的数学模型进行了辨识。结果表明：该方法辨识精度高，算法简单，易于计算机编程。