模拟真空羽流场的特征线法

介绍了特征线法（MOC法）对真空羽流场进行数值模拟的过程，用MOC法进行了喷管出口流动不均匀情况真空羽流场算例的求解，所得结果与理论分析相符合，与DSMC法计算结果吻合，提高了计算效率，节省了计算时间。结果表明：用MOC法可代替DSMC法对连续介质区进行模拟，并可在用N-S和DSMC耦合数值模拟中确定耦合边界。     

尾喷管内外超声速流场数值模拟

采用有限元法数值模拟带二次流可调收－扩喷管内外跨声速和超声速流场，在计算中采用四边形等参元确定位移函数，利用加权余量法中的Galerkin法建立有限元方程，数值研究三种型式喷管和四种飞行工况下尾喷管内外流场。计算结果与实验数据符合较好，说明可用此数值模拟技术帮助开展实验研究，进行尾喷管优化设计。     

基于高阶间断伽辽金数值方法的浸没边界法

为了降低复杂边界条件下贴体网格生成难度，近年来基于笛卡尔网格的浸没边界法逐渐成为研究此类问题的主要数值方法之一，然而计算精度和计算效率仍然是此类方法目前面临的挑战。与传统的空间2阶精度的有限体积格式相比，空间精度为3阶或3阶以上的高精度方法具有空间精度高、数值分辨率高、数值耗散小的优点，而作为高精度数值方法之一的间断伽辽金方法在浸没边界方面的应用仍较少。本文将高阶间断伽辽金方法的高精度优势与浸没边界法无需贴体网格的优势结合起来，提出了适用于可压缩流动的高精度浸没边界法。其中边界条件采用体积惩罚方式实现，同时采用牛顿法迭代以及MPI并行提高计算效率，物面的数据重建采用插值点处的逆距离权重（IDW-IP）方法代替高阶格式下多项式插值方法。本文基于笛卡尔网格测试验证了二维定常和非定常情况下的数值模拟效果，并与传统贴体网格的计算结果进行了对比。     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

台阶后紊流场的大涡模拟

用大涡模拟法对台阶后复杂的紊流场进行了数值研究。用有限体积法求解滤波N－S方程，亚格子紊流模型采用结构模型。数值计算中，对流项，扩散项和时间分别采用三阶SHARP中式，二阶中心差分格式和三阶Runge-Kutta格式，给出了台阶后紊流场的时间统计特性，瞬时特性和紊流强度分布，并与试验结果进行了比较。     

低温液体火箭燃烧室内射流雾化研究进展

从理论、实验和数值计算方面，阐述和总结了低温推进剂射流雾化的研究现状和进展。对射流雾化机理及形态进行了分类，重点从对比常温流体与低温流体射流雾化特征参数不同的角度，梳理了射流雾化方面的数值研究方法，并介绍了低温射流雾化的相关实验研究。结果表明:空气扰动破碎机理仍是最广为接受的理论，而超临界工况下射流雾化形态不符合现有的射流雾化分类；数值计算方法逐渐由流体体积法和水平集法向直接数值模拟转变，对闪蒸现象的建模是研究重难点之一；低温流体射流喷雾实验数据匮乏，需通过丰富测量手段，提高测量精度的途径获得完整、准确的实验数据。      

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

近地点变轨发动机高空喷管性能预示研究

对长征－2E火箭近地点变轨发动机EPKM高空喷管的地面性能和高空性能进行了预示。在喷管跨声速区，气相采用显式MacCormack差分格式、颗粒相采用特征线法，而在喷管超声速区采用特征线法，数值求解轴对称二维无粘两相流动模型。结合喷管内的气流分离准则预测发动机地面工作时的性能，同时根据地面试验数据外推发动机的高空性能，与实测性能数据比较，平均推力相对误差约为5.6%和1.5%左右；而直接对发动机高空工作时的满流状态喷管进行数值模拟所得的发动机平均推力与实测性能数据比较相对误差约为1.7％左右。研究表明，所采用的流动模型、气流分离准则和数值方法对高空喷管不同工况下的性能进行工程预示是有效的。     

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

燃烧室温度计算方法分析及应用

介绍了航空发动机燃烧室出口温度常用的两种计算方法，结合批产某型发动机试车数据进行计算，对两种计算结果进行对比分析。用涡轮流量法对某燃气发生器的燃烧室出口温度进行计算，其结果与外方给的数值基本一致。提出用涡轮流量法计算燃烧室出口温度较切合实际，为发动机寿命评估及其工作可靠性提供依据。     

基于能量法和特征值法的颤振预测数值方法研究

为研究适用于工程设计的颤振数值预测方法,采用能量法和特征值法对出现颤振现象的真实叶片进行预测。通过将2种数值方法的预测结果与试验结果进行对比,验证数值预测方法的精度和工程适用性。结果表明:2种方法都较准确的判断了某典型风扇转子叶片的气弹稳定性,与试验结果基本吻合,具有较高的工程实用价值。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

基于单层数值特征曲线法的复合材料接头强度分析

提出了基于单层的特征曲线法，利用Patran／Nastran有限元软件对复合材料机械连接接头强度进行了分析计算。根据提出的基于单层的数值特征曲线法，各个角度铺层的特征长度采用数值方法分别计算，根据各层的特征曲线计算各层的破坏载荷，最后确定接头的破坏载荷。通过算例分析，计算结果跟试验值吻合很好，较基于层合板特征曲线法所得结果更精确，说明了该方法的有效性。     

基于阻抗法的密封动力特性系数实验识别

搭建并优化了密封动力特性系数实验识别装置，应用阻抗法对梳齿密封动力特性系数开展实验识别，研究转速、进口压力和涡动频率对密封动力特性系数的影响。应用基于微元理论的密封动力特性系数识别方法对密封动力学特性进行了数值计算分析，与实验结果进行对比分析。结果表明：实验识别梳齿密封动力特性系数整体上略高于数值计算结果，但趋势上吻合良好，特别是表征系统稳定性的密封有效阻尼系数，数值与实验具有良好的一致性。梳齿密封动力特性系数在转子低频(<100 Hz)涡动时存在一定的频率依赖性，高频(>100 Hz)时阻尼系数的频率依赖性较弱。实验和数值模拟均表明，相较于转速，压比的变化对密封有效阻尼系数的影响更大。     

粘性流场三维NS方程有限体积法数值解

提出了有限体积法中的隐式ＮＮＤ格式数值求解三维全ＮＳ方程的薄层三维ＮＳ方程的方法，对绕带座舱的再入飞行器高超声速粘性流动进行数值模拟。     

三维带操纵面机翼欧拉方程数值模拟

本文给出了用欧拉方程求解带操纵面机翼的数值结果，采用有限体积法，四步Ｒｏｎｇｅ－Ｋｕｔｔａ时间推进，结合不连续面通量守恒传递技术，对有操纵面偏转的三维机翼流场进行求解，计算结果与实验符合较好。     

三维离散矢量场中的涡核识别

基于临界点理论，使用特征矢量方法在三维离散矢量场中进行了涡核寻找的研究。该方法能够快速准确地找到三维矢量场中的涡核，为CFD工作提供了一个有效的数值分析工具，搜寻结果与其它方法（如流线法等）进行了比较。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。