Well‐balanced discontinuous Galerkin method and finite volume WENO scheme based on hydrostatic reconstruction for blood flow model in arteries

The blood flow model in arteries admits the steady state solutions, for which the flux gradient is nonzero, and is exactly balanced by the source term. In this paper, by means of hydrostatic reconstruction, we construct a high order discontinuous Galerkin method, which exactly preserves the dead‐man steady state, which is characterized by a discharge equal to zero (analogue to hydrostatic equilibrium). Moreover, the method maintains genuine high order of accuracy. Subsequently, we apply the key idea to finite volume weighted essentially non‐oscillatory schemes and obtain a well‐balanced finite volume weighted essentially non‐oscillatory scheme. Extensive numerical experiments are performed to verify the well‐balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.     

High‐order continuous and discontinuous Galerkin methods for wave problems

Three Galerkin methods—continuous Galerkin, Compact Discontinuous Galerkin, and hybridizable discontinuous Galerkin—are compared in terms of performance and computational efficiency in 2‐D scattering problems for low and high‐order polynomial approximations. The total number of DOFs and the total runtime are used for this correlation as well as the corresponding precision. The comparison is carried out through various numerical examples. The superior performance of high‐order elements is shown. At the same time, similar capabilities are shown for continuous Galerkin and hybridizable discontinuous Galerkin, when high‐order elements are adopted, both of them clearly outperforming compact discontinuous Galerkin. Copyright © 2013 John Wiley & Sons, Ltd.     

Non‐linear flux‐splitting schemes with imposed discrete maximum principle for elliptic equations with highly anisotropic coefficients

Families of flux‐continuous, locally conservative, finite‐volume schemes have been developed for solving the general tensor pressure equation of petroleum reservoir simulation on structured and unstructured grids. The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization. Improved convergence using the quadrature parameterization has been established for the family of flux‐continuous schemes. When applied to strongly anisotropic full‐tensor permeability fields the schemes can fail to satisfy a maximum principle (as with other FEM and finite‐volume methods) and result in spurious oscillations in the numerical pressure solution. This paper presents new non‐linear flux‐splitting techniques that are designed to compute solutions that are free of spurious oscillations. Results are presented for a series of test‐cases with strong full‐tensor anisotropy ratios. In all cases the non‐linear flux‐splitting methods yield pressure solutions that are free of spurious oscillations. Copyright © 2010 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

Efficiency of high‐performance discontinuous Galerkin spectral element methods for under‐resolved turbulent incompressible flows

The present paper addresses the numerical solution of turbulent flows with high‐order discontinuous Galerkin methods for discretizing the incompressible Navier‐Stokes equations. The efficiency of high‐order methods when applied to under‐resolved problems is an open issue in the literature. This topic is carefully investigated in the present work by the example of the three‐dimensional Taylor‐Green vortex problem. Our implementation is based on a generic high‐performance framework for matrix‐free evaluation of finite element operators with one of the best realizations currently known. We present a methodology to systematically analyze the efficiency of the incompressible Navier‐Stokes solver for high polynomial degrees. Due to the absence of optimal rates of convergence in the under‐resolved regime, our results reveal that demonstrating improved efficiency of high‐order methods is a challenging task and that optimal computational complexity of solvers and preconditioners as well as matrix‐free implementations are necessary ingredients in achieving the goal of better solution quality at the same computational costs already for a geometrically simple problem such as the Taylor‐Green vortex. Although the analysis is performed for a Cartesian geometry, our approach is generic and can be applied to arbitrary geometries. We present excellent performance numbers on modern cache‐based computer architectures achieving a throughput for operator evaluation of 3·10⁸ up to 1·10⁹ DoFs/s (degrees of freedom per second) on one Intel Haswell node with 28 cores. Compared to performance results published within the last five years for high‐order discontinuous Galerkin discretizations of the compressible Navier‐Stokes equations, our approach reduces computational costs by more than one order of magnitude for the same setup.     

A family of multi‐point flux approximation schemes for general element types in two and three dimensions with convergence performance

A family of flux‐continuous, locally conservative, control‐volume‐distributed multi‐point flux approximation (CVD‐MPFA) schemes has been developed for solving the general geometry‐permeability tensor pressure equation on structured and unstructured grids. These schemes are applicable to the full‐tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full‐tensor flow approximation. The family of flux‐continuous schemes is characterized by a quadrature parameterization. Improved numerical convergence for the family of CVD‐MPFA schemes using the quadrature parameterization has been observed for structured and unstructured grids in two dimensions. The CVD‐MPFA family cell‐vertex formulation is extended to classical general element types in 3‐D including prisms, pyramids, hexahedra and tetrahedra. A numerical convergence study of the CVD‐MPFA schemes on general unstructured grids comprising of triangular elements in 2‐D and prismatic, pyramidal, hexahedral and tetrahedral shape elements in 3‐D is presented. Copyright © 2011 John Wiley & Sons, Ltd.     

High-order Finite Difference and Finite Volume WENO Schemes and Discontinuous Galerkin Methods for CFD

In recent years, high order numerical methods have been widely used in computational fluid dynamics (CFD), to effectively resolve complex flow features using meshes which are reasonable for today's computers. In this paper, we review and compare three types of high order methods being used in CFD, namely the weighted essentially non-oscillatory (WENO) finite difference methods, the WENO finite volume methods, and the discontinuous Galerkin (DG) finite element methods. We summarize the main features of these methods, from a practical user's point of view, indicate their applicability and relative strength, and show a few selected numerical examples to demonstrate their performance on illustrative model CFD problems.     

The third‐order polynomial method for two‐dimensional convection and diffusion

Using the upstream polynomial approximation a series of accurate two‐dimensional explicit numerical schemes is developed for the solution of the convection–diffusion equation. A third‐order polynomial approximation (TOP) of the convection term and a consistent second‐order approximation of the diffusion term are combined in a single‐step flux‐difference algorithm. Stability analysis confirms that the TOP‐12 scheme satisfies the CFL condition for two dimensions. Using smaller and narrower flux stencils compared to algorithms of similar accuracy, the TOP‐12 scheme is more efficient in terms of computations per single node. Numerical tests and comparison with other well‐known algorithms show a high performance of the developed schemes. Copyright © 2003 John Wiley & Sons, Ltd.     

High‐order accurate p‐multigrid discontinuous Galerkin solution of the Euler equations

Discontinuous Galerkin (DG) methods have proven to be perfectly suited for the construction of very high‐order accurate numerical schemes on arbitrary unstructured and possibly nonconforming grids for a wide variety of applications, but are rather demanding in terms of computational resources. In order to improve the computational efficiency of this class of methods a p‐multigrid solution strategy has been developed, which is based on a semi‐implicit Runge–Kutta smoother for high‐order polynomial approximations and the implicit Backward Euler smoother for piecewise constant approximations. The effectiveness of the proposed approach is demonstrated by comparison with p‐multigrid schemes employing purely explicit smoothing operators for several 2D inviscid test cases. Copyright © 2008 John Wiley & Sons, Ltd.     

One‐dimensional shock‐capturing for high‐order discontinuous Galerkin methods

Discontinuous Galerkin methods have emerged in recent years as an alternative for nonlinear conservation equations. In particular, their inherent structure (a numerical flux based on a suitable approximate Riemann solver introduces some stabilization) suggests that they are specially adapted to capture shocks. However, numerical fluxes are not sufficient to stabilize the solution in the presence of shocks. Thus, slope limiter methods, which are extensions of finite volume methods, have been proposed. These techniques require, in practice, mesh adaption to localize the shock structure. This is is more obvious for large elements typical of high‐order approximations. Here, a new approach based on the introduction of artificial diffusion into the original equations is presented. The order is not systematically decreased to one in the presence of the shock, large high‐order elements can be used, and several linear and nonlinear tests demonstrate the efficiency of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids

A novel class of weighted essentially nonoscillatory (WENO) schemes based on Hermite polynomials,termed as HWENO schemes,is developed and applied as limiters for high order discontinuous Galerkin (DG) method on triangular grids.The developed HWENO methodology utilizes high-order derivative information to keep WENO reconstruction stencils in the von Neumann neighborhood.A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils,making higher-order scheme stable and simplifying the reconstruction process at the same time.The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement.Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy,the designed HWENO limiters can simultaneously obtain uniform high order accuracy and sharp,essentially non-oscillatory shock transition.     

Spectral p‐multigrid discontinuous Galerkin solution of the Navier–Stokes equations

Discontinuous Galerkin (DG) methods are very well suited for the construction of very high‐order approximations of the Euler and Navier–Stokes equations on unstructured and possibly nonconforming grids, but are rather demanding in terms of computational resources. In order to improve the computational efficiency of this class of methods, a high‐order spectral element DG approximation of the Navier–Stokes equations coupled with a p‐multigrid solution strategy based on a semi‐implicit Runge–Kutta smoother is considered here. The effectiveness of the proposed approach in the solution of compressible shockless flow problems is demonstrated on 2D inviscid and viscous test cases by comparison with both a p‐multigrid scheme with non‐spectral elements and a spectral element DG approach with an implicit time integration scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

A γ‐DGBGK scheme for compressible multi‐fluids

The paper presents a Discontinuous Galerkin γ‐BGK (γ‐DGBGK) method for compressible multicomponent flow simulations by coupling the discontinuous Galerkin method with a γ‐BGK scheme based on WENO limiters. In this γ‐DGBGK method, the construction of the flux in the DG method is based on the kinetic scheme which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous terms in the flux formulation at cell interfaces. WENO limiters are used to obtain uniform high‐order accuracy and sharp non‐oscillatory shock transition, and time accuracy obtained by integration for the flux function at the cell interface. Numerical examples in one and two space dimensions are presented to illustrate the robust and accuracy of the present scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

discontinuous galerkin method for discontinuous temperature field problems

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允 许插值函数在单元边界上存在跳变. 在建立有限元方程时,通过在单元边界上引入数值通量 项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出. 数值算例 表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局 部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高. 该方法也为研 究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

适用于间断Galerkin方法的限制器研究

开发了一种适用于高精度间断Galerkin方法的斜率（多项式系数）限制器。与现有的斜率限制器不同,该限制器实施过程不考虑网格单元类型（三角形或四边形）,通过全微分方法构造新的多项式系数,因此,该限制器能够适用于各种类型网格——结构化网格、具有单一单元的非结构化网格和具有混合单元的非结构化网格。由于该限制器能够方便地应用于具有混合单元的非结构化网格,因此,本文使用的程序能够方便地求解具有复杂几何结构的流动问题。本文利用一些典型算例对其性能进行了验证,表明该限制器适用于不同类型的网格单元,能够在光滑解区保证高的精度,并能够在阊断区抑帛3非物理振荡。     

Implicit high-order discontinuous Galerkin method with HWENO type limiters for steady viscous flow simulations

Two types of implicit algorithms have been improved for high order discontinuous Galerkin (DG) method to solve compressible Navier-Stokes (NS) equations on triangular grids. A block lower-upper symmetric Gauss-Seidel (BLU-SGS) approach is implemented as a nonlinear iterative scheme. And a modified LU-SGS (LLU-SGS) approach is suggested to reduce the memory requirements while retain the good convergence performance of the original LU-SGS approach. Both implicit schemes have the significant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory (HWENO) limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock transition and the designed high-order accuracy simultaneously.     

High‐order implicit time integration for unsteady incompressible flows

The spatial discretization of unsteady incompressible Navier–Stokes equations is stated as a system of differential algebraic equations, corresponding to the conservation of momentum equation plus the constraint due to the incompressibility condition. Asymptotic stability of Runge–Kutta and Rosenbrock methods applied to the solution of the resulting index‐2 differential algebraic equations system is analyzed. A critical comparison of Rosenbrock, semi‐implicit, and fully implicit Runge–Kutta methods is performed in terms of order of convergence and stability. Numerical examples, considering a discontinuous Galerkin formulation with piecewise solenoidal approximation, demonstrate the applicability of the approaches and compare their performance with classical methods for incompressible flows. Copyright © 2011 John Wiley & Sons, Ltd.     

A third‐order implicit discontinuous Galerkin method based on a Hermite WENO reconstruction for time‐accurate solution of the compressible Navier–Stokes equations

A space and time third‐order discontinuous Galerkin method based on a Hermite weighted essentially non‐oscillatory reconstruction is presented for the unsteady compressible Euler and Navier–Stokes equations. At each time step, a lower‐upper symmetric Gauss–Seidel preconditioned generalized minimal residual solver is used to solve the systems of linear equations arising from an explicit first stage, single diagonal coefficient, diagonally implicit Runge–Kutta time integration scheme. The performance of the developed method is assessed through a variety of unsteady flow problems. Numerical results indicate that this method is able to deliver the designed third‐order accuracy of convergence in both space and time, while requiring remarkably less storage than the standard third‐order discontinous Galerkin methods, and less computing time than the lower‐order discontinous Galerkin methods to achieve the same level of temporal accuracy for computing unsteady flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Anisotropy favoring triangulation CVD(MPFA) finite‐volume approximations

A family of flux‐continuous, control‐volume distributed multi‐point flux approximation schemes CVD (MPFA) have been developed for solving the general geometry‐permeability tensor pressure equation on structured and unstructured grids (Comput. Geo. 1998; 2 : 259–290, Comput. Geo. 2002; 6 : 433–452). The locally conservative schemes are applicable to the diagonal and full‐tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full‐tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization. Improved numerical convergence for the family of CVD(MPFA) schemes for specified quadrature points has been observed for lower anisotropy ratios for both structured and unstructured grids in two dimensions. However, for strong full‐tensor anisotropy fields the quadrilateral schemes can induce strong spurious oscillations in the numerical solution. This paper motivates and demonstrates the benefit of using anisotropy favoring triangulation for treating such cases. Test examples involving strong full‐tensor anisotropy fields are presented in 2‐D and 3‐D, which show that the family of schemes on anisotropy favoring triangulation (prisms in 3‐D) yield well‐resolved pressure fields with little or no spurious oscillations. Copyright © 2010 John Wiley & Sons, Ltd.     

Implicit FEM‐FCT algorithms and discrete Newton methods for transient convection problems

A new generalization of the flux‐corrected transport (FCT) methodology to implicit finite element discretizations is proposed. The underlying high‐order scheme is supposed to be unconditionally stable and produce time‐accurate solutions to evolutionary convection problems. Its nonoscillatory low‐order counterpart is constructed by means of mass lumping followed by elimination of negative off‐diagonal entries from the discrete transport operator. The raw antidiffusive fluxes, which represent the difference between the high‐ and low‐order schemes, are updated and limited within an outer fixed‐point iteration. The upper bound for the magnitude of each antidiffusive flux is evaluated using a single sweep of the multidimensional FCT limiter at the first outer iteration. This semi‐implicit limiting strategy makes it possible to enforce the positivity constraint in a very robust and efficient manner. Moreover, the computation of an intermediate low‐order solution can be avoided. The nonlinear algebraic systems are solved either by a standard defect correction scheme or by means of a discrete Newton approach, whereby the approximate Jacobian matrix is assembled edge by edge. Numerical examples are presented for two‐dimensional benchmark problems discretized by the standard Galerkin finite element method combined with the Crank–Nicolson time stepping. Copyright © 2007 John Wiley & Sons, Ltd.