基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.      

UDMH/NTO火箭发动机尾焰流场特性数值仿真

采用AUSM(advection upstream splitting method)空间离散格式,探究了常用的3种k-ε湍流模型和2阶、3阶迎风格式对激波捕捉效果的影响。采用标准k-ε双方程模型和2阶迎风格式对不同高度和不同来流马赫数下的偏二甲肼/四氧化二氮(UDMH/NTO)火箭发动机尾焰流场进行喷管-尾焰流场一体化仿真。复燃反应采用12组分18步化学反应模型,喷管入口参数由热力计算给出。结果表明:随着来流马赫数的增加,波节数逐渐降低;随着飞行高度的上升,尾焰影响区域逐渐扩大;复燃反应造成混合区中的O2、N2质量分数大幅下降,而混合区中O、OH、NO质量分数则有所上升。      

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

空-空导弹用二元混压超声速进气道数值研究

对空空导弹固体火箭冲压发动机二元混压超声速进气道进行了数值研究。应用ＮＮＤ格式耦合对流迎风矢通量分裂技术（ＡＵＳＭ）,按照ＭａｃＣｏｒｍａｃｋ时间分裂方法对贴体坐标系下二维雷诺平均Ｎ－Ｓ方程进行有限差分离散。数值模拟了此进气道多种工况下的流动状态。结果表明,对于捕获激波和模拟附面层内流动均有很好的效果。在此基础上分析了来流马赫数、进气道工作高度及出口反压对进气道性能参数的影响。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

基于高阶逆风通量差分裂格式的直升机旋翼前飞流场模拟

描述了三维雷诺平均NS方程数值模拟前飞状态旋翼非定常流场的方法和模型.在该模型中,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,无需添加人工粘性,因而可有效地减少Jameson二阶中心差分有限体积法导致的旋翼尾迹数值耗散.为充分考虑旋翼尾迹对流场的影响和计入桨叶的旋转、挥舞、变距运动,发展了一种运动嵌套网格.针对嵌套网格中关键的贡献单元搜寻问题,提出了一种新的"伪贡献单元搜寻法".应用上述方法,对Caradonna模型旋翼在周期变距和挥舞运动条件下的前飞流场进行了计算,与采用Jameson二阶中心差分方法相比,计算结果更贴近实验值.并针对几种不同桨尖形状的旋翼前飞流场进行计算,表明后掠桨尖可有效减弱超临界流动.     

强制偏流喷管的数值仿真

采用薄层Ｎ－Ｓ方程成功地求解了强制流喷管内的三维复杂流场，对通量项和压力项用Ｖａｎ Ｌｅｅｒ发展的矢通量分裂法进行了离散，对离散方程采用线超松弛方法迭代求解，在亚音速与超音速气流共同存在的出口边界分情况进行了处理。计算结果验证了强制偏喷管具有自动补偿外界压力变化的能力，与实验结果基本吻合。     

基于计算电磁学方法的旋翼雷达散射截面特性

以Maxwell方程组为主控方程,时间离散采用4步Runge-Kutta法,空间离散使用MUSCL(monotonic upstream-centered scheme for conversation laws)格式与Steger-Warming通量分裂结合的方法,建立了一套基于时域有限体积法的旋翼RCS(雷达散射截面)特性数值计算方法.采用代数方法生成围绕旋翼的O型贴体、正交网格.在验证方法有效性基础上,着重分析了旋翼平面和翼型RCS极化、入射角、电尺寸等影响特性,并给出了翼型几何特点对RCS的影响规律;然后运用线性加权和法进行了旋翼翼型隐身性能的综合筛选.研究表明:旋翼RCS动态包络线是一个强散射水平的连续振荡区域;选取3片桨叶旋翼的最大雷达探测距离为2片桨叶的82.2%,且有利于控制旋翼的散射最大峰值.同时,HH02旋翼翼型在4个重要的散射角域,最大探测距离分别为NACA0012翼型的105.1%,99.4%,86.6%和83.5%,表明综合雷达隐身性能最好.     

耦合高效配平策略的旋翼气动特性分析方法

建立了一个耦合高效配平策略和计算流体力学（CFD）技术的旋翼气动特性分析方法.采用CFD求解器进行雅克比矩阵和旋翼流场的计算，建立了一个包含直接配平法的旋翼气动特性分析平台；引入差量法思想对配平策略进行简化，发展了一种差量配平法，解决了直接配平法中雅克比矩阵计算效率低下的问题，显著缩短了计算时间.旋翼流场求解中控制方程采用Navier-Stokes方程，空间离散格式选取2阶迎风Roe格式，时间推进采用隐式lower-upper symmetric Gauss-Seidel (LU-SGS)格式，湍流模型为Baldwin-Lomax (B-L)模型.分别采用建立的基于直接配平方法和差量配平方法的气动特性分析平台对AH-1G和Helishape 7A直升机的旋翼进行了配平计算和研究，对比分析了两种配平方法的计算效率和精度及背景网格的网格量对配平结果的影响.结果表明:与直接配平法相比，耦合差量配平法的气动特性分析方法计算圈数更少，可节省至少60%的计算时间；使用疏化的背景网格虽然不利于捕捉流场细节，但是可以在兼顾配平结果准确性的同时，降低配平计算代价，提高旋翼配平和流场计算的效率.      

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

倒置喷管的湍流流场计算

数值模拟采用基于 Van Leer的矢通量分裂格式 ;应用线 Gauss- Seidel迭代法求解 ;并结合局部时间步长技术对收敛过程加速。对 Re=1 .66× 1 0 6 的倒置喷管流场进行数值计算的结果与实验符合     

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

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

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

旋翼气动力数值模拟准确度影响因素计算

针对旋翼气动力CFD计算时的一些常见影响因素的选择问题,进行了数值计算讨论.计算求解器选择自研代码,该代码采用基于共享内存(OMP)的并行化方法,在嵌套网格系统下求解Navier-Stokes方程,空间方向上使用2阶迎风格式,时间方向采用双时间方法.讨论的问题包括背景网格数目对悬停桨尖涡捕捉、对前飞状态激波捕捉和气动力的影响,远场动量理论修正作用对计算结果影响等.对于讨论的不同网格数目情况下,计算的悬停状态旋翼拉力差异在6%左右,扭矩差异2.6%左右,旋翼前飞阻力和扭矩最大差异为2.58%和3.28%.是否采用动量修正旋翼拉力相差3%~5%.通过计算,给出了针对以上因素计算结果差异的范围,可以为旋翼气动力计算工作提供一定的参考.