一般状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型和高精度多介质Piecewise Parabolic Method(PPM)算法模拟一般形式Mie-Grüneisen状态方程的复杂多流体可压缩流动.根据多流体接触界面无振荡原则设计高精度计算格式,采用双波近似和两层迭代算法求解一般状态方程的Riemann问题,利用Strang分裂方法将一维算法直接推广到高维情况,不需要动能校正.二维多流体接触间断平移问题的模拟结果表明本算法在接触间断附近压力和速度没有振荡,界面数值耗散很小,被控制在两到三个网格之内.本文还模拟了二维Riemann问题和激波同多流体交界面的相互作用,算例结果表明本文方法可以有效地处理接触间断、激波等物理问题,且具有耗散小精度高的特点.     

SPH方法对气液两相流自由界面运动的追踪模拟

气液两相流动是自然界中常见的流动现象,对其进行数值模拟时要求必须能够准确跟踪界面运动变形。本文利用光滑粒子流体动力学方法,结合微可压缩模型(SCM),引入界面控制方法XSPH速度修正以及Van Derwaals状态方程修正,对典型的气液两相流动如二维溃坝、气泡上浮等问题进行了数值模拟,分析了空气和水相互作用机理以及界面运动规律,并同实验结果和其他数值计算结果进行了比较。结果表明,该方法在模拟多相介质界面运动问题准确有效,可用于处理更为复杂的多相流动工程问题。     

间断伽辽金方法在可压缩流数值模拟中的应用研究综述

本文对近三十年来,国内外对于高精度数值方法研究中的热点——间断伽辽金方法在可压缩流数值模拟方面的应用研究进行了综述。首先对间断伽辽金方法的基本概念和特点作了简单介绍,然后对应用该方法解决双曲型及椭圆型问题的发展历程进行了回顾,并重点梳理了其在计算流体力学领域可压缩流数值模拟方面的应用发展以及研究现状,之后对该方法在对应的网格技术、激波捕捉方法、湍流流动模拟以及计算量需求方面目前仍然存在的研究难点和可能的发展趋势做出了总结和分析。最后给出了间断伽辽金方法在可压缩流数值模拟中的若干应用实例。     

气液两相流混合模型代数重构方法

为提高国家数值风洞工程气液两相流软件混合模型求解器的自由界面模拟分辨率，增强多介质流体界面问题的模拟能力，基于代数重构方法，应用双曲正切函数界面捕捉格式对体积分数进行了插值重构，实现两相流混合模型与自由界面代数重构方法的耦合，抑制原混合模型下自由界面的非物理耗散现象，保持界面清晰尖锐；引入预处理技术，解决低速流动中由于特征值量级差距导致的计算刚性问题，提高收敛速度；应用双时间进行非定常时间推进。最后，采用方块对流、Zalesak盘刚性旋转和溃坝算例进行考核验证，数值结果显示，耦合了双曲正切函数界面捕捉格式的两相流混合模型在保持质量守恒的同时，获得了更高的自由界面分辨率，以及更丰富的多相流动细节。     

国家数值风洞(NNW)工程中的CFD基础科学问题研究进展

计算流体力学（CFD）基础理论方面的持续创新对数值模拟软件的功能拓展及推广应用具有十分重要的意义,国家数值风洞（NNW）工程在CFD相关的若干关键基础科学问题上开展了研究。简要综述了国内各参研团队近3年的研究进展,NNW工程在转捩与湍流模型及计算方法、多相多介质计算模型与方法、多物理场耦合计算模型与方法、高精度数值计算方法等方面形成了一系列阶段性的研究成果,建立了多种原创算法与模型。这些算法与模型中具有较高成熟度的将被集成至NNW工程的相关软件,并面向全国发布。     

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

罗茨鼓风机内部气流脉动的非定常数值分析

运用计算流体动力学方法对渐开线三叶罗茨鼓风机排气过程中的脉动进行了可压缩非定常流动的数值模拟, 在动静耦合区域采用了非一致网格界面, 用动网格方法实现了排气区域流场的容积变化, 通过将数值模拟结果与半经验理论公式计算结果的对比分析, 探讨了运转周期内气流脉动的分布特性和对排气流量的影响.      

自适应的高精度中心﹣WE NO混合格式

为对可压缩流动进行高精度高分辨率的数值模拟,基于模板的光滑度量函数,构造了新型的加权因子,在此基础上将7阶精度的WENO格式与8阶中心格式进行加权组合,构造了一种自适应低耗散的中心-WENO混合格式(H-WENO7-CD8),并采用Fourier方法对离散格式数学特性进行了理论分析。该自适应的高精度中心-WENO混合格式相比于七阶的WENO格式具有更小的耗散误差。通过对激波-密度干扰问题、Rayleigh-Taylor不稳定性问题和双Mach反射问题的数值模拟,并与WENO-JS格式的计算结果对比,结果表明:该格式结合了WENO格式和中心格式的特点,能更好地捕捉激波,对物理脉动也具有较高的分辨率,适合于复杂流体流动的数值模拟。     

气体动理学格式及其在再入问题中的应用

介绍了气体动理学格式GKS在飞行器高超声速再入过程涉及的跨流域稀薄及湍流等复杂流动问题中的拓展和应用。GKS利用BGK方程的通解来计算单元界面上的通量,耦合了分子的自由运动和相互碰撞,内涵多尺度特性。本文发展了耦合多种湍流模式的拓展GKS,能对典型高超声速湍流进行有效模拟。基于通量重构CPR框架发展了非结构网格上的高精度格式CPR-GKS,在可压缩黏性流动中可以同时保持高精度和良好的激波捕捉能力。为模拟航天器再入过程中的跨流域稀薄流动,发展了适合大规模并行计算的高效UGKS,并在再入问题中得到了很好的应用。研究揭示了GKS在再入问题模拟中的优异性能及其广阔的应用前景。     

高精度有限差分格式的色散优化及耗散控制

本文综述了我们在高精度有限差分格式的色散优化和耗散控制方面的研究进展。首先,我们提出了半离散有限差分格式色散和耗散相互独立的充分条件,实现了差分格式色散和耗散特性的独立调节。在此基础上,提出了色散最小、耗散可控的高精度差分格式,称为MDCD格式。MDCD格式已经得到了广泛应用,取得了很好的计算效果,但其主要缺点是耗散的调节依赖于经验。为了解决这一问题,我们进一步提出耗散的自适应调节方法。具有自适应耗散特性的高精度有限差分格式的基本特征是,差分格式的耗散能够随解的局部尺度自适应调节。为了构造这类格式,我们提出了一种新型的尺度识别器,它能够以等效无量纲波数的形式来定量衡量数值解的局部长度尺度。在此基础上,设计差分格式耗散参数与尺度识别器得到的等效无量纲波数之间的关系,从而构造了一类色散最小、耗散自适应的差分格式,称为MDAD格式。为了计算含有间断的问题,同时保持在光滑区的良好耗散特性,我们利用尺度识别器对一种经典的激波探测器进行改进,提出了一种新的激波探测器,并将自适应耗散格式与对应的WENO格式相混合,得到自适应耗散混合格式。近似色散关系显示该混合格式兼具高分辨率和鲁棒性。多个含间断流场的标准算例测试结果显示,自适应耗散混合格式具有良好的分辨率和激波捕捉能力。     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm,GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。     

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...     

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

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

High-order methods for the Euler and Navier–Stokes equations on unstructured grids

This article reviews several unstructured grid-based high-order methods for the compressible Euler and Navier–Stokes equations. We treat the spatial and temporal discretizations separately, hoping that it is easier to spot the similarities and differences of various types of methods. Our main focus is to present the basic design principles of each method, and highlight its pros and cons when appropriate. Sample computational results are shown to illustrate the capability of selected methods. These high-order methods are expected to be more efficient than low-order methods for problems requiring high accuracy, such as wave propagation problems, vortex-dominated flows including high-lift configuration, helicopter blade vortex interaction, as well as large eddy simulation and direct numerical simulation of turbulence. We conclude the paper with several current challenges in the proliferation of high-order methods in the aerospace community.     

High-order accurate,low numerical diffusion methods for aerodynamics

In recent years numerical methods have been widely used to effectively resolve complex flow features of aerodynamics flows with meshes that are reasonable for today's computers. High-order numerical methods were used mainly in direct numerical simulations and aeroacoustics. For many aeronautical applications, accurate computation of vortex-dominated flows is important because the vorticity in the flow field and the wake of swept wings at an incidence and rotor blades largely determines the distribution of loading. The main deficiency of widely available, second-order accurate methods for the accurate computation of these flows is the numerical diffusion of vorticity to unacceptable levels. Application of high-order accurate, low-diffusion numerical methods can significantly alleviate this deficiency of traditional second order methods. Furthermore, higher-order space discretizations have the potential to improve detached eddy simulation predictions of separated flows with significant unsteadiness. Recently developed high-order accurate finite-difference, finite-volume, and finite-element methods are reviewed. These methods can be used as an attractive alternative of traditional low-order central and upwind computational fluid dynamics methods for improved predictions of vortical and other complex, separated, unsteady flows. The main features of these methods are summarized, from a practical user's point of view, their applicability and relative strength is indicated, and examples from recent applications are presented to illustrate their performance on selected problems.     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

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

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