基于Newton/Gauss-Seidel迭代的DGM隐式方法

在Newton迭代方法的基础上,对高阶精度间断Galerkin有限元方法(DGM)的时间隐式格式进行了研究. Newton迭代 法的优势在于收敛效率高效,并且定常和非定常问题能够统一处理,对于非定常问题无需引入双时间步策略. 为了避免大型矩阵的求逆,采用一步Gauss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵,从而只需计算和存储对角(块)矩阵. 对角(块)矩阵采用数值方法计算. 空间离散采用Taylor基,其优势在于对于任意形状的网格,基函数的形式是一致的,有利于在混合网格上推广. 利用该方法,数值模拟了Bump绕流和NACA0012翼型绕流. 计算结果表明,与显式的Runge-Kutta时间格式相比,隐式格式所需的迭代步数和CPU时间均在很大程度上得到减少,计算效率能够提高1～ 2个量级.     

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

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

高雷诺数下求解NS方程的无网格算法

提出了一种适合高雷诺数NS方程求解的隐式无网格算法。针对高雷诺数粘性流动的特点,在附面层内的粘性影响区域采用法向层次推进布点的方法形成离散点云,在附面层外的计算区域内实行填充式布点的方法形成离散点云。根据附面层内外点云的不同构造特点,推导出运用格林公式和最小二乘曲面拟合方法求取空间导数的统一形式,在此基础上运用AUSM _up格式求得数值通量,并引入BL湍流模型对雷诺平均NS方程的湍流应力项进行封闭。时间推进格式方面,采用了计算效率较高的隐式高斯-赛德尔迭代算法。为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型低速流动、RAE2822翼型跨音速绕流和二维圆柱的分离流动进行了数值模拟。     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

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

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

基于全隐式无分裂算法求解三维N-S方程

基于多块结构网格,本文研究和发展了三维N-S方程的全隐式无分裂算法.对流项的离散运用Roe格式,粘性项的离散利用中心型格式.在每一次隐式时间迭代中,运用GMRES方法直接求解隐式离散引起的大型稀疏线性方组.为了降低内存需求以及矩阵与向量之间的运算操作数,Jacobian矩阵的一种逼近方法被应用在本文的算法之中.计算结果与实验结果基本吻合,表明本文的全隐式无分裂方法是有效和可行的.     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

基于非结构/混合网格模拟黏性流的高阶精度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 混合格式提升了计算效率和稳定性, 具有良好的应用前景.      

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

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

An implicit high-order preconditioned flux reconstruction method for low-Mach-number flow simulation with dynamic meshes

A fully implicit high-order preconditioned flux reconstruction/correction procedure via reconstruction (FR/CPR) method is developed to solve the compressible Navier-Stokes equations at low Mach numbers. A dual-time stepping approach with the second-order backward differentiation formula (BDF2) is employed to ensure temporal accuracy for unsteady flow simulation. When dynamic meshes are used to handle moving/deforming domains, the geometric conservation law is implicitly enforced to eliminate errors due to the resolution discrepancy between BDF2 and the spatial FR/CPR discretization. The large linear system resulted from the spatial and temporal discretizations is tackled with the restarted generalized minimal residual solver in the PETSc (portable, extensible toolkit for scientific computation) library. Through several benchmark steady and unsteady numerical tests, the preconditioned FR/CPR methods have demonstrated good convergence and accuracy for simulating flows at low Mach numbers. The new flow solver is then used to study the effects of Mach number on unsteady force generation over a plunging airfoil when operating in low-Mach-number flows. It is observed that weak compressibility has a significant impact on thrust generation but has a negligible effect on lift generation of an oscillating airfoil.     

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.     

A domain decomposition approach to finite volume solutions of the Euler equations on unstructured triangular meshes

We report on our recent efforts on the formulation and the evaluation of a domain decomposition algorithm for the parallel solution of two‐dimensional compressible inviscid flows. The starting point is a flow solver for the Euler equations, which is based on a mixed finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi‐discrete equations is obtained using a linearized backward Euler implicit scheme. As a result, each pseudo‐time step requires the solution of a sparse linear system for the flow variables. In this study, a non‐overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. First, we formulate an additive Schwarz algorithm using appropriate matching conditions at the subdomain interfaces. In accordance with the hyperbolic nature of the Euler equations, these transmission conditions are Dirichlet conditions for the characteristic variables corresponding to incoming waves. Then, we introduce interface operators that allow us to express the domain decomposition algorithm as a Richardson‐type iteration on the interface unknowns. Algebraically speaking, the Schwarz algorithm is equivalent to a Jacobi iteration applied to 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 our approach, the interface unknowns are numerical (normal) fluxes. Copyright © 2001 John Wiley & Sons, Ltd.     

An implicit meshless method for application in computational fluid dynamics

An implicit meshless scheme is developed for solving the Euler equations, as well as the laminar and Reynolds‐averaged Navier–Stokes equations. Spatial derivatives are approximated using a least squares method on clouds of points. The system of equations is linearised, and solved implicitly using approximate, analytical Jacobian matrices and a preconditioned Krylov subspace iterative method. The details of the spatial discretisation, linear solver and construction of the Jacobian matrix are discussed; and results that demonstrate the performance of the scheme are presented for steady and unsteady two dimensional fluid flows. Copyright © 2012 John Wiley & Sons, Ltd.     

基于WENO重构的真正多维黎曼求解器

具有良好守恒性与网格适应性的有限体积格式在流体力学的数值计算中占有重要地位。其中，求解数值流通量是实施有限体积法的关键步骤。一维情形下，通过求解局部黎曼问题来获得数值流通量的相关理论已经比较成熟。但是在计算多维问题时，传统的维度分裂方法仅考虑沿界面法向传播的信息，这不仅影响格式的精度，还可能会造成数值不稳定性从而诱发非物理现象。本文基于对流-压力通量分裂方法来构造真正多维的黎曼求解器，通过求解网格顶点处的多维黎曼问题来实现格式的多维特性。采用五阶WENO重构方法来获得空间的高阶精度，时间离散采用三阶TVD龙格-库塔格式。一系列数值实验的结果表明，真正多维的黎曼求解器不仅具有更高的分辨率还能有效克服多维强激波模拟中的数值不稳定性。     

An efficient algorithm for solving the incompressible fluid flow equations

The present paper reports on a modified pressure implicit predictor corrector type scheme for solving the flow governing equations, in which a consistent formulation is combined with a multi-grid solver for the pressure correction. In addition a parabolic sublayer (PSL) approach for the treatment of the flow in the vicinity of solid walls is critically evaluated in terms of accuracy and computational efficiency. The lid-driven cavity flow is chosen as the test case and results are presented for Reynolds numbers ranging from 100 to 1000. Predictions with the proposed scheme indicate substantial computational savings and fairly good agreement when compared with previous work. The PSL approach reduces the computing time, but with increasing Reynolds numbers the accuracy of the solutions tends to deteriorate.     

Positivity-preserving scheme for the flux reconstruction framework with a novel implementation

The positivity-preserving scheme for the flux reconstruction (FR) framework is studied in this paper. A simple and direct implementation of the scheme is also proposed. Owing to the scheme and the implementation, FR or other high-order nodal methods can provide completely physically meaningful solutions at each time step. Their robustness is greatly enhanced. In the numerical tests, the solving of Euler equations coupled with artificial viscosity and Navier–Stokes equations is focused. Satisfying results are obtained in problems that include strong shocks or vacuum states.     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

非稳态Hamilton-Jacobi方程的7阶加权紧致非线性格式

Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.