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

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

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

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

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

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

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

低马赫数流动中的预处理Euler方程的收敛特性

对预处理后Euler方程在低马赫数流动数值模拟中的收敛特性进行了细致研究。首先采用量级分析法研究了Weiss-Smith预处理方法对低速情形下Euler方程收敛性的影响;然后针对低速情形下预处理后的Euler方程组中各方程的收敛性,提出了一个改进的预处理矩阵。数值模拟结果表明：Weiss-Smith预处理方法可以很好地改善Euler方程在低速流动时的收敛特性,并且改进后的预处理矩阵在连续方程的收敛特性方面效果明显。     

运动边界非定常全速流场的数值模拟

将预处理引入到ALE格式的非定常可压缩N-S方程中,发展了一种能有效模拟含运动边界的三维非定常全速流场的数值方法.控制方程在非结构动态网格上用双时间步推进并用隐式、LU-SGS迭代求解;提出了一种特征变量的边界条件以适合预处理后的特征系统.几个典型算例的结果表明,预处理方法将可压缩流方法的计算能力拓展到了低马赫数范围,有效地提高了收敛速度和计算精度.     

低耗散TVD格式及叶轮机内低马赫数流动模拟

推导了条件化预处理矩阵在守恒变量下的形式,以及三维任意曲线坐标系下的相应的系统特征值与特征矩阵。根据预处理技术修正了任意曲线坐标下的Harten-TVD格式,得到低耗散TVD格式,并讨论了抑制了压力速度失耦的方法。分别使用修正后低耗散格式与未修正的格式,对高负荷低压涡轮叶栅T106内低速流动:进口马赫数分别为0.1,0.01和0.001,进行数值模拟。计算表明,未修正格式不能得到有意义的结果,而修正后的低耗散格式,即使在极低的马赫数流动下,也工作良好。同时可以注意到:对于无粘流动,三种进口马赫数条件下,流场内等压线分布几乎是一致的。      

新型高低速流统一算法及其叶轮机应用

 对于可压缩/不可压缩混合流动或高低速混合流动问题，目前没有较好的解决方法。基于Roe格式及预处理方法，提出了全速度Roe格式。全速度Roe格式是将Roe格式与低速Roe格式通过一个与本地马赫数相关的函数联合而得到的。与标准Roe格式相比，对亚音速流动，全速度Roe格式中声速对数值耗散的影响随着马赫数变化而变化。与预处理Roe格式相比，全速度Roe格式简单易于实现，而且没有计算稳定性的限制，不需要采用全局截断处理，从而在理论上对高低速流动具有一致精度。事实上，全速度Roe格式意味着可以通过直接改变Roe格式的非线性特征值，解决可压不可压流动统一计算精度的问题。这是对统一算法的新观点。数值模拟表明,全速度Roe格式具有与预处理Roe格式一样的激波捕获能力、极低速流动准确模拟能力。     

计算全马赫数下粘性流动的隐式算法

把预处理方法引入到求解三维粘性流动的Navier-Stokes方程隐式格式算法中,并应用这种方法对平板算例在层流和湍流两种不同流态下及四个不同Ma数下进行了数值模拟,计算结果与理论值及相关公式吻合,同时对比分析了有无预处理情况下的收敛性,证明了本文所采用的预处理隐式算法对于任意马赫数都是适合的,而且可以提高残差精度量级,加快收敛。     

基于控制理论和NS方程的气动设计方法研究

研究了基于控制理论和NS方程的气动设计方法,针对给定的目标函数表达形式,应用该设计理论在计算坐标下详细推导了相应的共轭方程及边界条件具体表达形式,以及梯度方程求解表达式,通过合理的数学变换,导出了共轭方程在笛卡尔坐标系下的直观表达形式,发展了有效的共轭方程数值求解方法,通过流动控制方程数值求解、共轭方程数值求解、目标函数对设计变量的梯度求解和优化算法等方面的有效结合,研究与发展了一种新的气动设计方法,以二维机翼气动设计为例,成功进行了亚、跨音速情形下的相关设计算例研究,研究结果表明应用控制理论和NS方程的气动设计方法在设计理论、适用性以及时间花费等方面都有着很好的特色和优点,且设计结果也更为可靠.     

低速预处理方法在悬停旋翼模拟中的应用

 在可压缩方程的求解当中引入时间导数低速预处理技术,构造了能够同时处理低速和高速流场的统一数值方法,采用该方法对Caradonna-Tung旋翼的亚声速悬停状态进行了模拟,其中方程采用在旋转坐标系下求解,网格采用重叠网格进行模拟。数值结果表明,在加入低速预处理方法以后,流场的收敛速度得到了明显提高,更为适合对下洗速度场等细节的模拟;另外,在加入低速预处理以后,方程的各个特征值都处在同一量级,避开了因特征值差异带来的数值耗散过大的问题,提高了数值计算的精度。     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

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

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

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

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

A low-diffusion robust flux splitting scheme towards wide-ranging Mach number flows

This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows. Based on Toro-Vázquez splitting approach, the new scheme splits inviscid flux into convective system and pressure system. This method introduces Mach number splitting function and numerical sound speed to evaluate advection system. Meanwhile, pressure diffusion term, pressure momentum flux, interface pressure and interface velocity are constructed to measure pressure system. Then, typical test problems are utilized to systematically assess the robustness and accuracy of the resulting scheme. Matrix stability analysis and a series of numerical cases, such as double shear-layer problem and hypersonic viscous flow over blunt cone, demonstrate that TVAP scheme achieves excellent low diffusion, shock stability, contact discontinuity and low-speed resolution, and is potentially a good candidate for wide-ranging Mach number flows.     

粘性极低速流动的预处理方法研究

讨论了采用预处理方法计算低速流动时遇到的参数选择、计算格式、舍入误差等问题,提出了一种适用于粘性计算的预处理参数计算方法.计算结果表明,使用高速流的计算格式计算低速问题时,即使残值收敛,也有可能得不到正确的计算结果;AUSM+-up格式适用于低速计算;马赫数小于等于10-6时,舍入误差会造成数值解无法收敛.