非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

燃烧室室壁冲击／气膜复合冷却壁温分布研究（二）理论部分

在考虑轴向导热的基础上，结合壁温分布实验数据，用纯气膜及冲击换热的研究成果，发展了计算纯气膜，夹层气膜及／气膜复合冷却火焰筒一维壁温分布计算程序，计算结果与实验结果二者基本吻合，研究表明，在冲击／气膜复合冷却计算中，应充分考虑Ｃｍ的值随Ｇｃ增大而增大的情况，且冲击／气膜冷却中的Ｃｍ取值高于纯气膜及夹层气膜冷却计算的Ｃｍ取值，本程序可用于纯气膜、夹层气膜及冲击／气膜复合冷却火焰筒一维壁温分布的初步分     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

一体化的重量计算程序设计浅析

提出编制重量计算一体化程序的必要性，通过介绍重量计算程序的设计思想，设计要点，特点及应用示例，对一体化的重量计算程序设计做了实阐述，可供工程技术人员进一步开发应用重量计算程序参考。     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

矢阵及其在航天器测控参数计算中的应用

本文将矢阵的概念及方法引入以航天器测控参数计算中，包括各种参数的坐标系统转换计算方法，航天器姿态矩阵的计算，姿态确定方法等，使其概念更清晰，表达更明了，推导更简捷，算法更适合于数值计算。     

机敏性指标RSD的计算与分析

对飞机敏性指标ＲＳＤ（后方分离距离尺度）进行了计算，通过计算得出，ＲＳＤ值受到飞机的滚转能力，俯仰能力及滚转时飞机阻力特性的综合影响。阻力特性中，由副翼偏度及滚转角速度引起的阻力对ＲＳＤ值起了重大作用，为此，从气动实验及计算中确定滚转中的飞机阻力特性的重要性是不可忽视的。     

固体发动机及火箭固有振动特性预估

用有限元方法计算了选用的固体火箭各级发动机壳体，药柱及整机固有频率和振型，并计算了火箭固有振动特性，为发动机冲击试验，今后的模态和振动试验提供了参考数据。这些数据是选择测试设备及量程，激励方式及位置，测点布置及试验夹具设计等测试技术所必需的。     

三维层流分离的N—S方程隐式上风格式数值模拟

数值模拟了绕半圆球柱M∞＝1．2，a＝19°的粘性分离流和涡流．用隐式LU－SGS及Roe的通量差分分裂格式数值求解了薄层N－S方程，分析计算了所得的流场结构，包括物体背风面的涡流及三维分离方式，计算结果同实验结果定性及定量吻合很好，另外，还用流动拓扑理论分析了计算所得的流场形态。     

波瓣喷管双层壁扩压器流场的数值分析

采用三维贴体曲线坐标网格，波瓣边界网格加密且正交，在整个计算区域进行全场计算。文中采用了Chen-Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法，对带有双层壁扩压器的波瓣喷管流动进行了冷热态数值计算和分析。对波瓣及双层壁，采用大粘性的方法解决流固耦合，计算结果表明：在双层壁间有外界冷气流被引射进入，形成了壁面的冷却气流，相对单层壁扩压器，双层壁扩压器的壁面温度明显降低；自波瓣出口截面沿流向产生的环流速度场，强化了主次流的掺混，速度分布渐趋均匀。通过试验测试，计算结果与实验数据符合良好，二者在离开波瓣420mm的混合管内相对主流速度的最大误差为15．35％。     

反导战术C4I系统中的拦截决策建模与仿真

在多枚来袭导弹情况下，形成反导拦截决策很容易引起计算内存混乱及CPU计算速率慢的问题。研究了应用面向对象及多线程编程技术，合理封装各个战术C^4I系统组成的对象，并利用多CPU的计算能力进行并行拦截决策。仿真表明，采用多线程偏程技术可使整个拦截决策过程更加可靠、迅速。     

某大展弦比机翼的刚度特性计算

通过对某新型民用飞机机翼结构进行分析、简化，建立计算模型。采用薄壁结构力学剖面刚度计算的理论进行编程计算，通过计算获得了翼面的刚心位置和弯曲刚度及扭转刚度，为某型飞机机翼结构设计和气动弹性分析提供了重要依据。     

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

舰载倾转旋翼机着舰流场及气动干扰特性数值计算研究

针对舰载倾转旋翼机甲板着舰时的气动干扰问题，开展了着舰流场及气动干扰特性数值计算研究。首先，采用CFD技术建立了基于滑移网格方法的舰船与倾转旋翼计算模型，并通过试飞及风洞试验结果对比，验证了计算方法的合理性。然后，进行了倾转旋翼机流场及舰机复合流场数值计算，研究了风向及气流速度对倾转旋翼机着舰气动干扰特性的影响规律。结果表明：来流风速越大，气流经机库陡壁下冲速度越高，导致旋翼有效迎角减小，旋翼拉力减小；全机滚转力矩受侧向气流干扰左侧升力增大时，引起向右的滚转力矩；全机俯仰力矩受气流偏角的影响，气流从前方流动产生抬头力矩，气流从后方流动产生低头力矩。     

战斗机燃油系统流体网络的数值计算

论述了战斗机燃油系统自供能源式流体网络的组成及工作原理，并对流体网络中的关键部件燃油增压泵，液动涡轮泵及射流泵等数学建模，而构成自供能源式流体网络的计算模型。在此基础上，通过分析几种状态下流体网络的工作情况，提出了计算仿真流体网络中流量与压力分布的方法，并对计算结果进行了分析和总结。     

直升机主减速装置传动效率飞行试验测定

介绍了一种直升机主减速装置传动效率的测量及计算方法。根据能量守恒定律，利用试验测量的主减速装置进、出口滑油温度和流量及输入轴功率等数据建立数学模型，计算出主减速装置的传动效率。该方法计算得出的结果与实际较符合，精度比较高，且测量方法简单可靠，可为主减速装置传动效率的鉴定提供参考。     

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

一种计算气膜冷却效率的新模型

提出了等效边界层的概念及“质量喷出率”的参数，并在此基础上提出了一种新的离散孔气膜冷却的掺混模型，解决了冷气喷出后边界层内速度分布波动大及高吹风比时冷气穿透边界层引起的计算困难，计算表明，与以往气膜冷却的掺混模型相比，新的掺混模型能适用于较高吹风比情况下的冷却效率的计算问题。     

航空发动机转子热弯曲稳态响应计算方法研究

研究航空发动机多转子系统热弯曲稳态响应计算方法，导出能计及诸多因素的热弯曲复数传递矩阵，给出热弯曲振动位能计算方法，并将系统位能作为热弯曲振动水平的一种衡量指标。建立了多转子系统热弯曲稳态响应和不平衡响应计算通用程序，计算了Jeffcoff转子、单转子试验器及多转子发动机系统等多种算例的热弯曲稳态响应及系统位能。并与不平衡响应规律进行比较分析。本计算方法有较大的工程应用价值。     

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

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