基于制造解的非结构二阶有限体积离散格式的精度测试与验证

随着计算机技术的飞速进步,计算流体力学得到迅猛发展,数值计算虽能够快速得到离散结果,但是数值结果的正确性与精度则需要通过严谨的方法来进行验证和确认.制造解方法和网格收敛性研究作为验证与确认的重要手段已经广泛应用于计算流体力学代码验证、精度分析、边界条件验证等方面.本文在实现标量制造解和分量制造解方法的基础上,通过将制造解方法精度测试结果与经典精确解(二维无黏等熵涡)精度测试结果进行对比,进一步证实了制造解精度测试方法的有效性,并将两种制造解方法应用于非结构网格二阶精度有限体积离散格式的精度测试与验证,对各种常用的梯度重构方法、对流通量格式、扩散通量格式进行了网格收敛性精度测试.结果显示,基于Green-Gauss公式的梯度重构方法在不规则网格上会出现精度降阶的情况,导致流动模拟精度严重下降,而基于最小二乘(least squares)的梯度重构方法对网格是否规则并不敏感.对流通量格式的精度测试显示,所测试的各种对流通量格式均能达到二阶精度,且各方法精度几乎相同;而扩散通量离散中界面梯度求解方法的选择对流动模拟精度有显著影响.     

非结构网格二阶有限体积法中黏性通量离散格式精度分析与改进

非结构网格二阶有限体积离散方法广泛应用于计算流体力学工程实践中,研究非结构网格二阶精度有限体积离散方法的计算精度具有现实意义. 计算精度主要受到网格和计算方法的影响,本文从单元梯度重构方法、黏性通量中的界面梯度计算方法两个方面考察黏性流动模拟精度的影响因素. 首先从理论上分析了黏性通量离散中的“奇偶失联”问题,并通过基于标量扩散方程的制造解方法验证了“奇偶失联”导致的精度下降现象,进一步通过引入差分修正项消除了“奇偶失联”并提高了扩散方程计算精度;其次,在不同类型、不同质量的网格上进行基于扩散方程的制造解精度测试,考察单元梯度重构方法、界面梯度计算方法对扩散方程计算精度的影响,结果显示,单元梯度重构精度和界面梯度计算方法均对扩散方程计算精度起重要作用;最后对三个黏性流动算例(二维层流平板、二维湍流平板和二维翼型近尾迹流动)进行网格收敛性研究,初步验证了本文的结论,得到了计算精度和网格收敛性均较好的黏性通量计算格式.      

非结构网格二阶有限体积法中黏性通量离散格式精度分析与改进

非结构网格二阶有限体积离散方法广泛应用于计算流体力学工程实践中,研究非结构网格二阶精度有限体积离散方法的计算精度具有现实意义.计算精度主要受到网格和计算方法的影响,本文从单元梯度重构方法、黏性通量中的界面梯度计算方法两个方面考察黏性流动模拟精度的影响因素.首先从理论上分析了黏性通量离散中的"奇偶失联"问题,并通过基于标量扩散方程的制造解方法验证了"奇偶失联"导致的精度下降现象,进一步通过引入差分修正项消除了"奇偶失联"并提高了扩散方程计算精度;其次,在不同类型、不同质量的网格上进行基于扩散方程的制造解精度测试,考察单元梯度重构方法、界面梯度计算方法对扩散方程计算精度的影响,结果显示,单元梯度重构精度和界面梯度计算方法均对扩散方程计算精度起重要作用;最后对三个黏性流动算例(二维层流平板、二维湍流平板和二维翼型近尾迹流动)进行网格收敛性研究,初步验证了本文的结论,得到了计算精度和网格收敛性均较好的黏性通量计算格式.     

基于混合网格多维梯度重构的热流预测方法研究

高超声速气动热环境的数值计算对算法和网格的敏感度极高. 随着高超声速飞行器外形日益复杂, 生成高质量的结构网格时间成本呈指数增加, 难以满足工程应用的需求. 非结构/混合网格因具有很强的复杂外形适应能力, 为了缩短任务周期, 有必要在非结构/混合网格上开展高精度的气动热环境数值计算方法研究. 梯度重构方法是影响非结构/混合网格热流计算精度的重要因素之一. 本文通过引入多维梯度重构方法, 发展了基于常规的非结构/混合网格的高精度热流计算方法, 对典型的高超声速Benchmark算例(二维圆柱)进行了模拟, 并与气动力计算广泛采用的Green-Gauss类方法和最小二乘类方法进行了对比. 计算结果表明, 多维梯度重构方法能有效提高非结构/混合网格热流预测精度, 其鲁棒性和收敛性更好. 最后将多维梯度重构方法应用于常规混合网格的三维圆柱和三维双椭球绕流问题, 得到了与实验值吻合较好的热流计算结果, 展现了良好的应用前景.      

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD) 方法和软件已经在航空航天飞行器设计中发挥了重要的作用, 但是由于二阶精度格式的耗散和色散较大, 对于湍流、分离等多尺度流动现象的模拟, 现有成熟的CFD 软件仍难以给出满意的结果, 为此CFD 工作者发展了众多的高阶精度计算格式. 如果以适应的计算网格来分类, 一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法. 由于非结构/混合网格具有良好的几何适应性, 基于非结构/混合网格的高阶精度格式近年来备受关注. 本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展, 重点介绍了空间离散方法, 主要包括k-Exact 和ENO/WENO 等有限体积方法, 间断伽辽金(DG) 有限元方法, 有限谱体积(SV) 和有限谱差分(SD) 方法, 以及近来发展的各种DG/FV 混合算法和将各种方法统一在一个框架内的CPR (correctionprocedure via reconstruction) 方法等. 随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题, 包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等. 在综述过程中, 介绍了各种方法的优势与不足, 其间介绍了作者发展的基于"静动态混合重构" 的DG/FV 混合算法. 最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.     

动网格生成技术及非定常计算方法进展综述

对应用于飞行器非定常运动的数值计算方法(包括动态网格技术和相应的数值离散格式)进行了综述.根据网格拓扑结构的不同,重点论述了基于结构网格的非定常计算方法和基于非结构/混合网格的非定常计算方法,比较了各种方法的优缺点.在基于结构网格的非定常计算方法中,重点介绍了刚性运动网格技术、超限插值动态网格技术、重叠动网格技术、滑移动网格技术等动态结构网格生成方法,同时介绍了惯性系和非惯性系下的控制方程,讨论了非定常时间离散方法、动网格计算的几何守恒律等问题.在基于非结构/混合网格的非定常计算方法中,重点介绍了重叠非结构动网格技术、重构非结构动网格技术、变形非结构动网格技术以及变形/重构耦合动态混合网格技术等方法,以及相应的计算格式,包括非定常时间离散、几何守恒律计算方法、可压缩和不可压缩非定常流动的计算方法、各种加速收敛技术等.在介绍国内外进展的同时,介绍了作者在动态混合网格生成技术和相应的非定常方法方面的研究与应用工作.     

基于非结构/混合网格模拟黏性流的高阶精度DDG/FV混合方法

间断Galerkin有限元方法(discontinuous Galerkin method, DGM) 因具有计算精度高、模板紧致、易于并行等优点, 近年来已成为非结构/混合网格上广泛研究的高阶精度数值方法. 但其计算量和内存需求量巨大, 特别是对于网格规模达到百万甚至数千万的大型三维实际复杂外形问题, 其计算量和存储量对计算资源的消耗是难以承受的. 基于“混合重构”的DG/FV 格式可以有效降低DGM 的计算量和存储量. 本文将DDG 黏性项离散方法推广应用于DG/FV 混合算法, 得到新的DDG/FV混合格式, 以进一步提高DG/FV混合算法对于黏性流动模拟的计算效率. 通过Couette流动、层流平板边界层、定常圆柱绕流, 非定常圆柱绕流和NACA0012 翼型绕流等二维黏性流算例, 优化了DDG 通量公式中的参数选择, 验证了DDG/FV 混合格式对定常和非定常黏性流模拟的精度和计算效率, 并与广泛使用的BR2-DG 格式的计算结果和效率进行对比研究. 一系列数值实验结果表明, 本文构造的DDG/FV混合格式在二维非结构/混合网格的Navier-Stokes 方程求解中, 在达到相同的数值精度阶的前提下, 相比BR2-DG格式, 对于隐式时间离散的定常问题计算效率提高了2 倍以上, 对于显式时间离散的非定常问题计算效率提高1.6 倍, 并且在一些算例中, 混合格式具有更优良的计算稳定性. DDG/FV 混合格式提升了计算效率和稳定性, 具有良好的应用前景.      

AUFS格式在无网格方法中的应用

将计算量小,激波分辨率高的AUFS(artificially upstream flux vector splitting)格式应用于无网格方法.所发展算法基于多项式基函数最小二乘无网格方法,采用线性基函数曲面拟合及AUFS格式计算各离散点的空间导数,应用四阶Runge--Kutta法进行时间显式推进.为验证算法健壮性、精度以及计算效率,对Riemann问题、超音速平面流动,以及不同攻角NACA0012翼型跨音速流场进行了数值模拟,其结果同采用HLLC(Harten-Lax-van Leer-contact)格式的无网格方法以及文献报道结果吻合较好,并且计算量较形式简单HLLC格式减少约15%.     

A high-order Runge-Kutta discontinuous Galerkin method with a subcell limiter on adaptive unstructured grids for two-dimensional compressible inviscid flows

A robust, adaptive unstructured mesh refinement strategy for high-order Runge-Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two-dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock-affected cells to control numerical spurious oscillations. An efficient cell-by-cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks.     

非结构混合网格在气动热标模模拟中的应用

非结构混合网格消除了结构网格节点的结构性限制,可以较好地处理边界,同时兼顾了粘性边界层模拟的需求,具有灵活性大、对复杂外形适应能力强和生成耗时短等优点,在飞行器气动特性模拟中得到广泛应用.本文针对非结构混合网格的特点,把前期针对非结构混合网格气动力高精度模拟发展改进的梯度计算方法和Roe格式熵修正方法推广应用到气动热流的数值模拟.以典型钝锥标模外形的高超声速绕流为研究对象,开展了不同网格形式和第一层网格不同间距的影响研究.结果 表明,热流计算时,头部物面网格最好采用四边形或四边形交叉剖分得到的三角形网格,物面法向的网格雷诺数取20左右,为热流计算时非结构混合网格的生成提供了指导,同时验证了计算方法的有效性和可靠性.     

A second-order upwind finite-volume method for the Euler solution on unstructured triangular meshes

A scheme for the numerical solution of the two-dimensional (2D) Euler equations on unstructured triangular meshes has been developed. The basic first-order scheme is a cell-centred upwind finite-volume scheme utilizing Roe's approximate Riemann solver. To obtain second-order accuracy, a new gradient based on the weighted average of Barth and Jespersen's three-point support gradient model is used to reconstruct the cell interface values. Characteristic variables in the direction of local pressure gradient are used in the limiter to minimize the numerical oscillation around solution discontinuities. An Approximate LU (ALU) factorization scheme originally developed for structured grid methods is adopted for implicit time integration and shows good convergence characterisitics in the test. To eliminate the data dependency which prohibits vectorization in the inversion process, a black-gray-white colouring and numbering technique on unstructured triangular meshes is developed for the ALU factorization scheme. This results in a high degree of vectorization of the final code. Numerical experiments on transonic Ringleb flow, transonic channel flow with circular bump, supersonic shock reflection flow and subsonic flow over multielement aerofoils are calculated to validate the methodology.     

An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

An h-adaptive solution of the spherical blast wave problem

Shock waves and contact discontinuities usually appear in compressible flows, requiring a fine mesh in order to achieve an acceptable accuracy of the numerical solution. The usage of a mesh adaptation strategy is convenient as uniform refinement of the whole mesh becomes prohibitive in three-dimensional (3D) problems. An unsteady h-adaptive strategy for unstructured finite element meshes is introduced. Non-conformity of the refined mesh and a bounded decrease in the geometrical quality of the elements are some features of the refinement algorithm. A 3D extension of the well-known refinement constraint for 2D meshes is used to enforce a smooth size transition among neighbour elements with different levels of refinement. A density-based gradient indicator is used to track discontinuities. The solution procedure is partially parallelised, i.e. the inviscid flow equations are solved in parallel with a finite element SUPG formulation with shock capturing terms while the adaptation of the mesh is sequentially performed. Results are presented for a spherical blast wave driven by a point-like explosion with an initial pressure jump of 10⁵ atmospheres. The adapted solution is compared to that computed on a fixed mesh. Also, the results provided by the theory of self-similar solutions are considered for the analysis. In this particular problem, adapting the mesh to the solution accounts for approximately 4% of the total simulation time and the refinement algorithm scales almost linearly with the size of the problem.     

A novel multidimensional solution reconstruction and edge‐based limiting procedure for unstructured cell‐centered finite volumes with application to shallow water dynamics

On unstructured meshes, the cell‐centered finite volume (CCFV) formulation, where the finite control volumes are the mesh elements themselves, is probably the most used formulation for numerically solving the two‐dimensional nonlinear shallow water equations and hyperbolic conservation laws in general. Within this CCFV framework, second‐order spatial accuracy is achieved with a Monotone Upstream‐centered Schemes for Conservation Laws‐type (MUSCL) linear reconstruction technique, where a novel edge‐based multidimensional limiting procedure is derived for the control of the total variation of the reconstructed field. To this end, a relatively simple, but very effective modification to a reconstruction procedure for CCFV schemes, is introduced, which takes into account geometrical characteristics of computational triangular meshes. The proposed strategy is shown not to suffer from loss of accuracy on grids with poor connectivity. We apply this reconstruction in the development of a second‐order well‐balanced Godunov‐type scheme for the simulation of unsteady two‐dimensional flows over arbitrary topography with wetting and drying on triangular meshes. Although the proposed limited reconstruction is independent from the Riemann solver used, the well‐known approximate Riemann solver of Roe is utilized to compute the numerical fluxes, whereas the Green–Gauss divergence formulation for gradient computations is implemented. Two different stencils for the Green–Gauss gradient computations are implemented and critically tested, in conjunction with the proposed limiting strategy, on various grid types, for smooth and nonsmooth flow conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Multi‐material incompressible flow simulation using the moment‐of‐fluid method

This paper compares the numerical performance of the moment‐of‐fluid (MOF) interface reconstruction technique with Youngs, LVIRA, power diagram (PD), and Swartz interface reconstruction techniques in the context of a volume‐of‐fluid (VOF) based finite element projection method for the numerical simulation of variable‐density incompressible viscous flows. In pure advection tests with multiple materials MOF shows dramatic improvements in accuracy compared with the other methods. In incompressible flows where density differences determine the flow evolution, all the methods perform similarly for two material flows on structured grids. On unstructured grids, the second‐order MOF, LVIRA, and Swartz methods perform similarly and show improvement over the first‐order Youngs' and PD methods. For flow simulations with more than two materials, MOF shows increased accuracy in interface positions on coarse meshes. In most cases, the convergence and accuracy of the computed flow solution was not strongly affected by interface reconstruction method. Published in 2009 by John Wiley & Sons, Ltd.