基于Navier Stokes数值计算的风洞收缩段设计

运用Navier—Stokes数值模拟对不同收缩曲线和设计方案的风洞收缩段和试验段的流动进行模拟，并对收缩段的分离特性、出口速度均匀性以及试验段的速度均匀性和附面层厚度作出评价。计算格式在空间上采用中心有限体积离散，在时间上采用多步Runge-Kutta时间步长格式进行积分。     

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

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

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

MPI＋OpenMP算法在三维可压缩流场计算中的应用

在多核CPU集群并行体系结构下,采用MPI＋OpenMP的混合并行算法,对可压缩流场进行数值模拟,并在计算时间上与MPI算法进行比较。流场计算的控制方程为Euler方程,空间离散采用Jameson中心差分格式,时间离散采用R-K法,并行模式采用Master/Slave模式。通过对M6机翼和某弹丸采用多种并行方法进行流场的数值模拟,得出MPI＋OpenMP混合算法在一定条件下具有高效性的结论。     

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

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

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

预处理方法在低速流场数值模拟中的应用

自20世纪70年代以来，随着高分辨率差分格式和高计算效率时间推进格式的建立和发展，可压缩流动数值模拟取得了很大进展。并广泛采用时间推进法。但是。对于低速不可压流动，采用时间推进法效果并不太好。目前对低速不可压流动的数值模拟主要有三种途径：（1）基于压力的算法，如著名的SIMPLE算法（压力偶合方程组的半隐式方法）；     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

三维侧压高超声速进气道试验与计算研究

针对三维侧压高超声速进气道开展了试验与计算研究。试验在中国航天空气动力技术研究院的高超声速脉冲风洞中进行,在试验中,来流马赫数为6,文中采用三维有限体积并行计算程序对风洞试验条件下的流动进行了数值模拟,采用k-g湍流模型,时间迭代采用LU隐式格式,对流通量采用了AUSM+、AUSMPW两种不同的空间离散格式进行了对比计算研究,结果表明两个不同的格式都能较好地模拟进气道流场。顶板中心线和离开中心线0.02m位置处的压力分布计算值和试验结果基本一致,表明AUSM+、AUSMPW格式的激波分辨率较高,粘性耗散比较低。计算结果还显示在侧压板壁面出现激波诱导边界层分离现象。     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

高超声速马赫反射

通过实验和数值计算模拟了激波绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项，粘 英和化学源项分别采用迎风ＴＶＤ格式，中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

二元矢量喷管流场及性能数值模拟

用有限差分法对二元矢量喷管的内流场及推力性能进行数值模拟。流场计算采用ＭａｃＣｏｒ－ｍａｃｋ两步格式，运用两层代数紊流模型，对初场、网格进行处理，采用加快收敛的当地时间步长，使计算精度和时间都能较好地满足工程应用的要求，得到了与实验数据较一致的喷管内特性。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

跨声速风洞机翼半模型试验洞壁边界层影响的数值模拟研究

采用数值模拟方法，模拟垫块法风洞试验二维垫块厚度对半模型试验结果的影响，计算模型采用M6机翼半模型。控制方程采用全场三维可压、雷诺平均的Navier—Stokes方程；采用Jameson中心有限体积、多步Runge—Kutta时间步长格式进行求解；湍流模型采用Baldwin—Lomax模型。通过本文的研究，可以为跨声速风洞半模型试验垫块厚度的选取提供理论计算方法和设计依据。     

跨音速翼型试验侧壁干扰数值模拟研究

采用分区网格技术，数值模拟NACA0012翼型在二维风洞中的侧壁干扰特性。控制方程采用全场三维可压、雷诺平均的Navier—Stokes方程；计算格式采用Jameson中心有限体积、多步Runge—Kutta时间步长格式进行求解；湍流模型采用Baldwin—Lomax模型。通过分析不同迎角、马赫数和展弦比情况下的流动状态，反映了由于侧壁干扰而产生的二维风洞流场的三维性。