低速流数值模拟的低耗散预处理格式

在预处理Roe格式中加入可调参数,给出了一种低耗散预处理Roe格式,并分析了可调参数的取值规律.在极低速区,低耗散预处理Roe格式反应了流动的实际耗散;在与来流速度相近的区域,低耗散预处理R0e格式可以保证流场的稳定性.对低马赫数、低雷诺数定常圆柱绕流和低马赫数、高雷诺数NACA0012翼型绕流所做的数值实验,得到了比...     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

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

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

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

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

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

A novel high-order scheme for numerical simulation of wake flow over helicopter rotors in hover

Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor. A new high-order scheme(we call it WENO-K) proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors. WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction, and its advantage is to reduce dissipation in smoo...     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。     

烧蚀外形转捩流动的数值模拟

本文采用雷诺平均Navier-Stokes方程和修正Baldwin-Lomax代数湍流模型,以Roe二阶流通量差分分裂格式进行离散,数值模拟了典型烧蚀外形的高超声速流场,重点讨论了流动转捩的判定及其与复杂流态的关联.结果表明:数值模拟出的球钝锥外形转捩点位置与工程估算结果符合较好;对工程估算中存在困难的奶头状凹陷烧蚀外形,数值模拟也可给出符合流场结构的转捩点位置.     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.