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

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

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。     

平面不可压后台阶湍流流动的数值模拟——壁面温度的影响

对平面不可压后台阶湍流流动进行了数值模拟,研究了台阶后下壁面温度对流场的影响.数值模拟采用SIMPLE算法,湍流采用k-ω湍流模型模拟.研究结果表明台阶后通道温度剖面分为大梯度区、过渡区和稳定区,对于温度剖面的过渡区在后台阶流动的回流区具有凸的特性,而在恢复区具有凹的特性.另外,数值解表明湍流模型在模拟具有分离再附着流动特性的后台阶湍流流动时具有很高的精度.     

分离流动问题的三层模型解法

在三层模型的理论范围内,用底层方程数值求解后向台阶和压缩拐角的分离流动问题,得到了与N-S方程数值方法及风洞实验值较一致的结果。文中采用三阶迎风格式和加松驰因子的Gauss-Seidel迭代,对提高数值算法的精度和收敛速度是可行的.     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

涡轮高压导叶流场结构及损失分析

采用具有TVD性质有三阶精度Godumov格式，数值模拟了某型涡轮高压导叶流场，并对流场损失及涡系结构进行分析，结果表明由于强烈的径向流，使低能流体大量地迁移到下端壁，使得下端部的损失显著增大，从而引起了总损失的增加。因此在高压导叶设计中，当采用低展弦比叶片时，应注意削弱叶片出口附近流场的严重径向窜流。     

压气机叶片排内三维紊流流动的隐式高分辨率数值分析

对跨音速压气机叶片排内具有复杂三维激波的紊流流场所进行的数值分析，采用了空间三阶精度无波动格式以保证对紊流特征的正确模拟，利用ＬＵ－ＳＧＳ算法避免了一般隐式算法的繁琐运算，使之与显式方法推进一步的运算工作量基本相当，从而构成了一种既准确又高效的跨音速压气机紊流流动数值求解系统。计算与实验进行了对比，结果是令人满意的。     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

球心非对称的钝圆头体近区绕流场气动特性分析研究

本文采用数值模拟对球心非对称的钝圆头体附近流场气动特性进行分析研究。数值求解三维可压缩NS方程，采用NND格式，以修正的Baldwin-Lomax代数模型对不同飞行状态进行数值模拟，分析了球心非对称钝圆头体对流场的影响，与部分试验比较结果满意。     

叶片扩压器内三维流场数值计算与实验研究

本文采用加罚有限元方法对离心压缩机叶片扩压器内流场进行了数值分析。在数值求解过程中,为了避免流动对流占优所引起的压力场波动,控制方程离散时借鉴有限差分法的技巧,进行权函数修正,通过迎风参数来控制沿流动方向上所施加的流线迎风量,保证数值求解的稳定性。叶片扩压器内紊流流动的数值计算结果表明:在其他项的假扩散被降至最小时,压力项的震荡亦被消弱。紊流场的计算采用k-ε紊流模型。同时本文对叶片扩压器内流场进行了实验测量,通过数值计算结果与实验测量值的比较分析,结果是令人满意的。最后就提高叶片扩压器的静压效率和压缩机整机性能提出了有益的建议。      

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

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

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

台阶后横喷氢气超音速燃烧流场数值模拟研究

采用ＭｃＣｏｒｍａｃｋ格式、代数涡粘性湍流模型及有限速率化学动力学模型，用数值模拟方法研究了台阶后横喷氢气二元燃烧流场。数值计算结果表明：台阶的作用不仅能扩大火焰稳定性，而且增加氢射流对主流的穿透深度，提高燃烧效率；当进口气流Ｍ数越高，进口温度越高，油气当量比越接近于恰当比，壁面温度对流场的影响越大。还提出了在超音速气流中，横喷氢自动着火时滞的火焰稳定机理新观点，由此可更准确地预估自动着火点的位置及火焰稳定性。     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

带扰流片的翼型绕流流场的计算

文章用基于伪压缩方法的有效不可压Navier-Stokes解算器，（该解算器可推广应用于非定常问题），以时间精度的方式对带扰流片的翼型绕流流场进行了数值模拟。对大多数计算而言，空间离散多采用三阶通量分裂格式，对得到的矩阵方程用点SGS松弛格式来求解。用单区和多区网格生成方法分别对流场进行了数值计算，并对得到的数值计算结果进行了比较。该程序由流过NACA 0012翼型的定常湍流进行了进一步验证，并且计算得到的表面压力系数与物理试验和其它数值计算的结果符合得很好。而且该程序还显示出很好的矢量化和并行化特点得到进一步证明。计算用基本翼型是先进的波音运输机预研翼型。其几何尺寸在网格生成时稍微有些改动。文章还研究了扰流片偏转角固定为60。。并研究了攻角对流场的影响。扰流片端部和翼型后缘非定常涡脱落过程也得到了很好的模拟，同时还显示了在不同时间步时的流型。在湍流模型的局限内，计算得到的平均量和物理试验结果符合得很好。     

微型飞行器小展弦比机翼的低雷诺数气动力特性分析

以固定机翼式微型飞行器的气动设计为研究背景,对小展弦比薄机翼的低雷诺数气动力特性进行了数值分析.不可压粘性绕流的求解采用了人工压缩性方法,其中对流项的离散应用了三阶迎风格式.湍流模型为Baldwin-Barth一方程模型.一无弯度薄翼的数值计算结果表明本文的数值方法是成功的.通过数值分析,揭示了弯度和展弦比对气动力特性的影响,也给出了不同弯度机翼流动结构的差别.     

高速飞行弹丸诱导的三维爆轰波数值模拟

采用三维非定常化学反应控制方程组对驻定在带攻角高速飞行弹丸上的斜爆轰波流场进行了数值模拟。对流项和化学反应项采用时间分裂法,分别用TVD格式和基元反应模型进行处理。计算结果表明三维数值模拟可以更准确的描述弹丸飞行攻角对爆轰波的形成及结构的影响,是进一步研究超驱爆轰推进技术的重要手段。     

乘方和 QUICK 格式在非对称扩压器内流动计算的比较

采用Ｓｐｅｚｉａｌｅ的非线性ｋ－ε模型计算了非对称扩压器内的流动,着重用数值方法考察了两种典型的对流差分格式（乘方格式和ＱＵＩＣＫ格式）对模型预测性能的影响,得出两种紊流模型对对流差分格式并不敏感的结论。     

适用于直升机飞行力学分析的三维空间大气紊流模型

发展了一种新的适用于直升机飞行力学分析的三维空间大气紊流模型。采用Dryden紊流模型生成给定空间点的时域离散紊流速度,在速度坐标系框架下,横向和垂向采用统计特性守恒的高斯插值法进行紊流场空间扩展,而沿速度的纵向采用时间序列延迟法扩展紊流场,最终形成覆盖直升机各个气动面的三维空间大气紊流速度场。将本文计算结果和国外经过飞行仿真验证的二维空间紊流模型结果进行对比分析,验证本文模型的有效性和正确性,研究二维和三维空间紊流模型的区别。结果表明:本文模型改善了二维空间大气紊流模型对俯仰角速度、横向速度和航向角速度响应计算幅值偏小的缺陷;紊流引起的直升机响应幅值随飞行速度和紊流强度的增加而增大,随飞行高度的增加而减小。     

用蒙特卡罗法仿真生成三维空间大气紊流场

基于三维紊流场相关函数矩阵用蒙特卡罗法仿真生成三维空间大气紊流场,并用理论证明和紊流场数值检验的方法,证明这样的空间大气紊流场较好地符合所选用的大气紊流模型 (如 Dryden模型),且具有良好的均匀性和各向同性。三维空间大气紊流场可用于飞机飞行仿真,尤其是多架飞机在大气紊流中飞行的仿真 (如空中加油、编队飞行)。