非结构网格上p型多重网格法流场数值模拟

在二维、三维非结构网榕上,针对间断Galerkin方法计算量大、收敛慢的缺点将p型多重网格方法应用于该方法求解跨音速Euler方程,提高计算效率。p型多重网格方法是通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛。文中对高阶近似(p＞0)使用显式格式,最低阶近似(p=0)采用隐式格式。NACA0012翼型和O...     

用混合网格及各向异性多重网格法求解三维可压紊流流动

采用混合网格求解紊流Navier Stokes方程。在物面附近采用柱状网格 ,其他区域则采用完全非结构网格。方程的求解采用Jamson的有限体积法 ,紊流模型采用两层Baldwin Lomax代数紊流模型。用各向异性多重网格法来加速解的收敛。数值算例表明 ,用混合网格及各向异性多重网格求解紊流流动是非常有效的     

三维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解三维非定常变系数对流扩散方程的一种高精度全隐紧致差分格式,该格式在空间上具有四阶精度,时间具有二阶精度。为了克服传统迭代法在每一个时间步上迭代求解隐格式时收敛速度慢的缺点,采用多重网格加速技术,建立了适用于本文高精度全隐紧致格式的多重网格算法,从而大大加快了迭代收敛速度。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

基于聚合多重网格方法的跨音速颤振的数值模拟

采用流固耦合方法对跨音速颤振进行了数值模拟。流体方面在非结构网格上用有限体积方法求解了Euler方程;结构方面则求解了后掠机翼典型剖面的结构模态方程。时间推进采用双时间步长：对每一真实时间步,都通过基于聚合多重网格方法的伪时间步推进,对流体和结构方程交替迭代．得到一个稳态的流固耦合的解。文章最后给出了NACA64A010翼型剖面的跨音速颤振边界．与相关文献的计算结果符合良好。     

一种有效的广义特征值分析方法

提出了一种适合于自适应有限元分析中求解广义特征值问题的多重网格方法．这种方法充分利用了初始网格下的结果，通过插值或最小二乘拟合技术来得到网格变化后的新的近似特征向量，然后由多重网格迭代过程实现对结构广义特征值问题的求解．在多重网格迭代的光滑步中，选择了收敛梯度法以提高其收敛率；在粗网格校正步中，则导出了一种近似求解特征向量误差的方程．这种方法将网格离散过程和数值求解过程很好地相结合，建立了一个网格细分后广义特征值问题的快速重分析方法，与传统有限元方法相比较，具有计算简便、计算量少等特点，可以作为结构动力问题自适应有限元分析的一种十分有效的工具．     

基于高阶和高效格式的悬停旋翼可压缩无粘绕流的计算

发展了一种基于高精度和高效格式计算悬停旋翼复杂绕流的隐式有限体积方法。控制方程为Euler方程,其中对流项通量的左右状态量采用五阶加权基本无振荡(WENO)格式重构,时间推进应用隐式LU-SGS算法,为进一步加速收敛,采用三层V循环多重网格松弛方法。考虑到多重网格方法的思想以及五阶WENO格式涉及6个网格单元,建议仅在最细网格上使用WENO格式。计算结果表明本文方法能有效捕捉激波,对尾迹也有较高分辨率,基于隐式LU-SGS算法的多重网格方法能有效提高计算效率。     

二维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解二维非定常变系数对流扩散方程的一种时间二阶、空间四阶精度的三层全隐紧致差分格式。为了加快迭代求解隐格式时在每一个时间步上的收敛速度，采用多重网格加速技术，建立了适用于本文高精度金隐紧致格式的多重网格算法。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

Walsh函数有限体积法的多重网格特征研究

Walsh函数有限体积法(FVM-WBF)是一种能够在网格内部捕捉间断的新型数值方法. 持续增加Walsh基函数数目能够稳步提高FVM-WBF方法的求解分辨率, 但计算量暴发式增长和收敛速度下降的问题也会同步出现. 针对Walsh基函数数目增加而引起的计算效率问题, 本文分析了Walsh基函数及其系数所能影响的网格单元局部均值区域尺度, 发现其中隐含类似多重网格的尺度特征, 据此提出一种结合多重网格策略的FVM-WBF方法. 在定常流场计算中根据各级Walsh基函数影响尺度的不同, 对每级Walsh基函数设置满足其稳定性约束的时间步长, 在时间推进求解的过程中快速消除不同波长的数值误差, 实现多重网格的加速收敛效果. 选取NACA0012翼型和二维圆柱的定常无黏绕流问题作为算例, 对引入多重网格策略的FVM-WBF方法和不考虑多重网格策略的FVM-WBF方法进行对比测试. 数值结果证实: 新发展的FVM-WBF方法具备多重网格的关键特征, 在不增加任何特殊处理和计算量的情况下, 只需通过时间步长的调整, 就能够达到多重网格的加速效果, 显著提升计算效率.      

求解三维流固耦合问题的一种全隐全耦合区域分解并行算法

三维流固耦合问题的非结构网格数值算法在很多工程领域都有重要应用, 目前现有的数值方法主要基于分区算法, 即流体和固体区域分别进行求解, 因此存在收敛速度较慢以及附加质量导致的稳定性问题, 此外, 该类算法的并行可扩展性不高, 在大规模应用计算方面也受到一定限制．本文针对三维非定常流固耦合问题, 提出一种基于区域分解的全隐全耦合可扩展并行算法．首先基于任意拉格朗日?欧拉框架建立流固耦合控制方程, 然后时间方向采用二阶向后差分隐式格式、空间方向采用非结构稳定化有限元方法进行离散．对于大规模非线性离散系统, 构造一种结合非精确Newton法、Krylov子空间迭代法与区域分解Schwarz预条件子的Newton-Krylov-Schwarz (NKS) 并行求解算法, 实现流体、固体和动网格方程的一次性整体求解．采用弹性障碍物绕流的标准测试算例对数值方法的准确性进行了验证, 数值性能测试结果显示本文构造的全隐全耦合算法具有良好的稳定性, 在不同的物理参数下具有良好的鲁棒性, 在“天河二号”超级计算机上, 当并行规模从192增加到3072个处理器核时获得了91%的并行效率．性能测试结果表明本文构造的NKS算法有望应用于复杂区域流固耦合问题的大规模数值模拟研究中．      

多重网格法与Newton-Raphson法耦合求线弹流数值解

利用多重网格算法,在最粗网格上,采用Newton-Raphson迭代法模式精解线接触弹流非线性方程组,充分利用了多重网格法与Newton-Raphson法各自的优点,计算实践表明,求解弹流问题的数值解过程在收敛与稳定性方面均有较大改善,且有相当宽的载荷参数适用范围。     

Symmetric Gauss-Seidel multigrid solution of the Euler equations on structured and unstructured grids

The efficient symmetric Gauss-Seidel (SGS) algorithm for solving the Euler equations of inviscid, compressible flow on structured grids, developed in collaboration with Jameson of Stanford University, is extended to unstructured grids. The algorithm uses a nonlinear formulation of an SGS solver, implemented within the framework of multigrid. The earlier form of the algorithm used the natural (lexicographic) ordering of the mesh cells available on structured grids for the SGS sweeps, but a number of features of the method that are believed to contribute to its success can also be implemented for computations on unstructured grids. The present paper reviews, the features of the SGS multigrid solver for structured gr0ids, including its nonlinear implementation, its use of “absolute” Jacobian matrix preconditioning, and its incorporation of multigrid, and then describes the incorporation of these features into an algorithm suitable for computations on unstructured grids. The implementation on unstructured grids is based on the agglomerated multigrid method developed by Sørensen, which uses an explicit Runge-Kutta smoothing algorithm. Results of computations for steady, transonic flows past two-dimensional airfoils are presented, and the efficiency of the method is evaluated for computations on both structured and unstructured meshes.     

A Hybrid Domain Decomposition and Multigrid Method for the Acceleration of Compressible Viscous Flow Calculations on Unstructured Triangular Meshes

This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.     

Turbulent flow calculations using unstructured and adaptive meshes

A method of efficiently computing turbulent compressible flow over complex two-dimensional configurations is presented. The method makes use of fully unstructured meshes throughout the entire flow field, thus enabling the treatment of arbitrarily complex geometries and the use of adaptive meshing techniques throughout both viscous and inviscid regions of the flow field. Mesh generation is based on a locally mapped Delaunay technique in order to generate unstructured meshes with highly stretched elements in the viscous regions. The flow equations are discretized using a finite element Navier-Stokes solver, and rapid convergence to steady state is achieved using an unstructured multigrid algorithm. Turbulence modelling is performed using an inexpensive algebraic model, implemented for use on unstructured and adaptive meshes. Compressible turbulent flow solutions about multiple-element aerofoil geometries are computed and compared with experimental data.     

Unstructured volume-agglomeration MG: Solution of the Poisson equation

We are interested in solving second-order PDEs with multigrid and unstructured meshes. The multigrid strategy we present here is adapted from the generalized finite volume agglomeration multigrid algorithm we have developed recently for the solution of the Euler equations. We now focus on Poisson's equation. A strategy is defined by introducing a correction factor for the diffusive terms, and some illustrating results are given.     

Multigrid flow solutions in complex two-dimensional geometries

A finite difference solution algorithm is described for use on two-dimensional curvilinear meshes generated by the solution of the transformed Laplace equation. The efficiency of the algorithm is improved through the use of a full approximation scheme (FAS) multigrid algorithm using an extended pressure correction scheme as smoother. The multigrid algorithm is implemented as a fixed V-cycle through the grid levels with a constant number of sweeps being performed at each grid level. The accuracy and efficiency of the numerical code are validated using comparisons of the flow over two backward step configurations. Results show close agreement with previous numerical predictions and experimental data. Using a standard Cartesian co-ordinate flow solver, the multigrid efficiency obtainable in a rectangular system is shown to be reproducible in two-dimensional body-fitted curvilinear co-ordinates. Comparisons with a standard one-grid method show the multigrid method, on curvilinear meshes, to give reductions in CPU time of up to 93%.     

A Coupled Level Set and Volume of Fluid method for automotive exterior water management applications

Motivated by the need for practical, high fidelity, simulation of water over surface features of road vehicles a Coupled Level Set Volume of Fluid (CLSVOF) method has been implemented into a general purpose CFD code. It has been implemented such that it can be used with unstructured and non-orthogonal meshes. The interface reconstruction step needed for CLSVOF has been implemented using an iterative ‘clipping and capping’ algorithm for arbitrary cell shapes and a re-initialisation algorithm suitable for unstructured meshes is also presented. Successful verification tests of interface capturing on orthogonal and tetrahedral meshes are presented. Two macroscopic contact angle models have been implemented and the method is seen to give very good agreement with experimental data for a droplet impinging on a flat plate for both orthogonal and non-orthogonal meshes. Finally the flow of a droplet over a round edged channel is simulated in order to demonstrate the ability of the method developed to simulate surface flows over the sort of curved geometry that makes the use of a non-orthogonal grid desirable.     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

An analytical interface reconstruction algorithm in the PLIC‐VOF method for 2D polygonal unstructured meshes

Piecewise linear interface calculation (PLIC) schemes have been extensively employed in the volume‐of‐fluid (VOF) method for interface capturing in numerical simulations of multiphase flows. Polygonal unstructured meshes are often adopted because of their geometric flexibility and superiority in gradient calculation. An analytical interface reconstruction algorithm in the PLIC‐VOF method for arbitrary convex polygonal cells has been proposed in this study. The line interface at a given orientation within a polygonal cell is located by an analytical technique. It has been tested successfully for four different geometric shapes that are common in polygonal meshes. The computational efficiency of the present algorithm has been compared with several published schemes in the literature. The proposed algorithm has been shown to yield higher accuracy with reduction in computational complexity. A numerical simulation of a dam‐breaking problem has been performed using the proposed analytical PLIC technique on polygonal meshes. The results are in good agreement with experimental data available in the literature, which serves as a demonstration of its performance in a real multiphase flow.     

An algebraic multigrid solver for Navier-Stokes problems

An efficient numerical method is presented for solving the equations of motion for viscous fluids. The equations are discretized on the basis of unstructured finite element meshes and then solved by direct iteration. Advective fluxes are temporarily fixed at each iteration to provide a linearized set of coupled equations which are then also solved by iteration using a fully implicit algebraic multigrid (AMG) scheme. A rapid convergence to machine accuracy is achieved that is almost mesh-independent. The scaling of computing time with mesh size is therefore close to the optimum.     

非定常激波三维双楔面反射的数值研究

采用基于MUSCL-Hancock插值的有限体积方法, 在非结构自适应网格上求解三维Euler方程, 研究了非定常激波在三维垂直双楔面上的反射现象. 研究结果表明, 由于三维效应的影响, 通过二维非定常激波反射理论预测的三维马赫反射区域范围存在着一定的局限, 在该区域范围内出现了一种非典型的三维突起结构, 即第二类三维马赫干. 另外, 对于不同的激波马赫数和楔面倾角组合, 非定常激波在三维双楔面上会形成具有四波结构的三维马赫反射或具有三波结构的三维规则反射.