A high-order discontinuous Galerkin method for the incompressible Navier-Stokes equations on arbitrary grids

A high-order discontinuous Galerkin (DG) method is proposed in this work for solving the two-dimensional steady and unsteady incompressible Navier-Stokes (INS) equations written in conservative form on arbitrary grids. In order to construct the interface inviscid fluxes both in the continuity and in the momentum equations, an artificial compressibility term has been added to the continuity equation for relaxing the incompressibility constraint. Then, as the hyperbolic nature of the INS equations has been recovered, the local Lax-Friedrichs (LLF) flux, which was previously developed in the context of hyperbolic conservation laws, is applied to discretize the inviscid term. Unlike the traditional artificial compressibility method, in this work, the artificial compressibility is introduced only for the construction of the inviscid numerical fluxes; therefore, a consistent discretization of the INS equations is obtained, irrespective of the amount of artificial compressibility used. What is more, as the LLF flux can be obtained directly and straightforward, no numerical iteration for solving an exact Riemann problem is entailed in our method. The viscous term is discretized by the direct DG method, which was developed based on the weak formulation of the scalar diffusion problems on structured grids. The performance and the accuracy of the method are demonstrated by computing a number of benchmark test cases, including both steady and unsteady incompressible flow problems. Due to its simplicity in implementation, our method provides an attractive alternative for solving the INS equations on arbitrary grids.     

An immersed discontinuous Galerkin method for compressible Navier-Stokes equations on unstructured meshes

We introduce an immersed high-order discontinuous Galerkin method for solving the compressible Navier-Stokes equations on non–boundary-fitted meshes. The flow equations are discretised with a mixed discontinuous Galerkin formulation and are advanced in time with an explicit time marching scheme. The discretisation meshes may contain simplicial (triangular or tetrahedral) elements of different sizes and need not be structured. On the discretisation mesh, the fluid domain boundary is represented with an implicit signed distance function. The cut-elements partially covered by the solid domain are integrated after tessellation with the marching triangle or tetrahedra algorithms. Two alternative techniques are introduced to overcome the excessive stable time step restrictions imposed by cut-elements. In the first approach, the cut-basis functions are replaced with the extrapolated basis functions from the nearest largest element. In the second approach, the cut-basis functions are simply scaled proportionally to the fraction of the cut-element covered by the solid. To achieve high-order accuracy, additional nodes are introduced on the element faces abutting the solid boundary. Subsequently, the faces are curved by projecting the introduced nodes to the boundary. The proposed approach is verified and validated with several two- and three-dimensional subsonic and hypersonic low Reynolds number flow applications, including the flow over a cylinder, a space capsule, and an aerospace vehicle.     

Large eddy simulation of incompressible free round jet with discontinuous Galerkin method

In the paper, discontinuous Galerkin method is applied to simulation of incompressible free round turbulent jet using large eddy simulation with eddy viscosity approach. The solution algorithm is based on the classical projection method, but instead of the solution of the Poisson equation, a parabolic equation is advanced in pseudo‐time, which provides the pressure field ensuring the proper pressure–velocity coupling. For time and pseudo‐time integration, explicit Runge–Kutta method is employed. The computational meshes consist of hexahedral elements with flat faces. Within a given finite element, all flow variables are expressed with modal expansions of the same order (including velocity and pressure). Discretisation of the viscous terms in the Navier–Stokes equations and Laplacian in the Poisson equation is stabilised with mixed finite element approach. The correctness of the solution algorithm is verified in a commonly used test case of laminar flow in 3D lid‐driven cavity. The results of computations of the free jet are compared with experimental and numerical reference data, the latter obtained from the high‐order pseudospectral code. The statistics of centerline flow velocity – mean velocity and its fluctuations – show satisfactory agreement with the reference data. Copyright © 2015 John Wiley & Sons, Ltd.     

A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows

Both compressible and incompressible Navier-Stokes solvers can be used and are used to solve incompressible turbulent flow problems. In the compressible case, the Mach number is then considered as a solver parameter that is set to a small value, M ≈0.1, in order to mimic incompressible flows. This strategy is widely used for high-order discontinuous Galerkin (DG) discretizations of the compressible Navier-Stokes equations. The present work raises the question regarding the computational efficiency of compressible DG solvers as compared to an incompressible formulation. Our contributions to the state of the art are twofold: Firstly, we present a high-performance DG solver for the compressible Navier-Stokes equations based on a highly efficient matrix-free implementation that targets modern cache-based multicore architectures with Flop/Byte ratios significantly larger than 1. The performance results presented in this work focus on the node-level performance, and our results suggest that there is great potential for further performance improvements for current state-of-the-art DG implementations of the compressible Navier-Stokes equations. Secondly, this compressible Navier-Stokes solver is put into perspective by comparing it to an incompressible DG solver that uses the same matrix-free implementation. We discuss algorithmic differences between both solution strategies and present an in-depth numerical investigation of the performance. The considered benchmark test cases are the three-dimensional Taylor-Green vortex problem as a representative of transitional flows and the turbulent channel flow problem as a representative of wall-bounded turbulent flows. The results indicate a clear performance advantage of the incompressible formulation over the compressible one.     

A reconstructed discontinuous Galerkin method for multi-material hydrodynamics with sharp interfaces

Discontinuous Galerkin (DG) methods have been well established for single-material hydrodynamics. However, consistent DG discretizations for non-equilibrium multi-material (more than two materials) hydrodynamics have not been extensively studied. In this work, a novel reconstructed DG (rDG) method for the single-velocity multi-material system is presented. The multi-material system being considered assumes stiff velocity relaxation, but does not assume pressure and temperature equilibrium between the multiple materials. A second-order DG(P₁) method and a third-order least-squares based rDG(P₁P₂) are used to discretize this system in space, and a third-order total variation diminishing (TVD) Runge-Kutta method is used to integrate in time. A well-balanced DG discretization of the non-conservative system is presented and is verified by numerical test problems. Furthermore, a consistent interface treatment is implemented, which ensures strict conservation of material masses and total energy. Numerical tests indicate that the DG and rDG methods are, indeed, the second- and third-order accurate. Comparisons with the second-order finite volume method show that the DG and rDG methods are able to capture the interfaces more sharply. The DG and rDG methods are also more accurate in the single-material regions of the flow. This work focuses on the general multidimensional rDG formulation of the non-equilibrium multi-material system and a study of properties of the method via one-dimensional numerical experiments. The results from this research will be the foundation for a multidimensional high-order rDG method for multi-material hydrodynamics.     

A discontinuous Galerkin method for the Navier–Stokes equations

The foundations of a new discontinuous Galerkin method for simulating compressible viscous flows with shocks on standard unstructured grids are presented in this paper. The new method is based on a discontinuous Galerkin formulation both for the advective and the diffusive contributions. High‐order accuracy is achieved by using a recently developed hierarchical spectral basis. This basis is formed by combining Jacobi polynomials of high‐order weights written in a new co‐ordinate system. It retains a tensor‐product property, and provides accurate numerical quadrature. The formulation is conservative, and monotonicity is enforced by appropriately lowering the basis order and performing h‐refinement around discontinuities. Convergence results are shown for analytical two‐ and three‐dimensional solutions of diffusion and Navier–Stokes equations that demonstrate exponential convergence of the new method, even for highly distorted elements. Flow simulations for subsonic, transonic and supersonic flows are also presented that demonstrate discretization flexibility using hp‐type refinement. Unlike other high‐order methods, the new method uses standard finite volume grids consisting of arbitrary triangulizations and tetrahedrizations. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical comparison of hybridized discontinuous Galerkin and finite volume methods for incompressible flow

In this study, a hybridizable discontinuous Galerkin method is presented for solving the incompressible Navier–Stokes equation. In our formulation, the convective part is linearized using a Picard iteration, for which there exists a necessary criterion for convergence. We show that our novel hybridized implementation can be used as an alternative method for solving a range of problems in the field of incompressible fluid dynamics. We demonstrate this by comparing the performance of our method with standard finite volume solvers, specifically the well‐established finite volume method of second order in space, such as the icoFoam and simpleFoam of the OpenFOAM package for three typical fluid problems. These are the Taylor–Green vortex, the 180‐degree fence case and the DFG benchmark. Our careful comparison yields convincing evidence for the use of hybridizable discontinuous Galerkin method as a competitive alternative because of their high accuracy and better stability properties. Copyright © 2014 John Wiley & Sons, Ltd.     

FOURIER NONLINEAR GALERKIN APPROXIMATION FOR THE TWO DIMENSIONAL NAVIER－STOKES EQUATIONS

ＦＯＵＲＩＥＲＮＯＮＬＩＮＥＡＲＧＡＬＥＲＫＩＮＡＰＰＲＯＸＩＭＡＴＩＯＮＦＯＲＴＨＥＴＷＯＤＩＭＥＮＳＩＯＮＡＬＮＡＶＩＥＲ－ＳＴＯＫＥＳＥＱＵＡＴＩＯＮＳＨｏｕＹａｎｒｅｎ（侯延仁）（ＲｅｃｅｉｖｅｄＪｕｎｅ１．１９９５．Ｃｏｍｍｕｎｉｃａｔｅｄ...     

Revisit of dilation-based shock capturing for discontinuous Galerkin methods

The idea of using velocity dilation for shock capturing is revisited in this paper, combined with the discontinuous Galerkin method. The value of artificial viscosity is determined using direct dilation instead of its higher order derivatives to reduce cost and degree of difficulty in computing derivatives. Alternative methods for estimating the element size of large aspect ratio and smooth artificial viscosity are proposed to further improve robustness and accuracy of the model. Several benchmark tests are conducted, ranging from subsonic to hypersonic flows involving strong shocks. Instead of adjusting empirical parameters to achieve optimum results for each case, all tests use a constant parameter for the model with reasonable success, indicating excellent robustness of the method. The model is only limited to third-order accuracy for smooth flows. This limitation may be relaxed by using a switch or a wall function. Overall, the model is a good candidate for compressible flows with potentials of further improvement.     

A locally conservative and energy-stable finite-element method for the Navier-Stokes problem on time-dependent domains

We present a finite-element method for the incompressible Navier-Stokes problem that is locally conservative, energy-stable, and pressure-robust on time-dependent domains. To achieve this, the space-time formulation of the Navier-Stokes problem is considered. The space-time domain is partitioned into space-time slabs, which in turn are partitioned into space-time simplices. A combined discontinuous Galerkin method across space-time slabs and space-time hybridized discontinuous Galerkin method within a space-time slab results in an approximate velocity field that is H(div)-conforming and exactly divergence-free, even on time-dependent domains. Numerical examples demonstrate the convergence properties and performance of the method.     

A high-order nodal discontinuous Galerkin method to solve preconditioned multiphase Euler/Navier-Stokes equations for inviscid/viscous cavitating flows

In this study, a high-order accurate numerical method is applied and examined for the simulation of the inviscid/viscous cavitating flows by solving the preconditioned multiphase Euler/Navier-Stokes equations on triangle elements. The formulation used here is based on the homogeneous equilibrium model considering the continuity and momentum equations together with the transport equation for the vapor phase with applying appropriate mass transfer terms for calculating the evaporation/condensation of the liquid/vapor phase. The spatial derivative terms in the resulting system of equations are discretized by the nodal discontinuous Galerkin method (NDGM) and an implicit dual-time stepping method is used for the time integration. An artificial viscosity approach is implemented and assessed for capturing the steep discontinuities in the interface between the two phases. The accuracy and robustness of the proposed method in solving the preconditioned multiphase Euler/Navier-Stokes equations are examined by the simulation of different two-dimensional and axisymmetric cavitating flows. A sensitivity study is also performed to examine the effects of different numerical parameters on the accuracy and performance of the solution of the NDGM. Indications are that the solution methodology proposed and applied here is based on the NDGM with the implicit dual-time stepping method and the artificial viscosity approach is accurate and robust for the simulation of the inviscid and viscous cavitating flows.     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

Divergence‐free discontinuous Galerkin schemes for the Stokes equations and the MAC scheme

We show that recently studied discontinuous Galerkin discretizations in their lowest order version are very similar to the marker and cell (MAC) finite difference scheme. Indeed, applying a slight modification, the exact MAC scheme can be recovered. Therefore, the analysis applied to the DG methods applies to the MAC scheme as well and the DG methods provide a natural generalization of the MAC scheme to higher order and irregular meshes. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonlinear Galerkin method for the exterior nonstationary Navier-Stokes equations

ＩｎｔｒｏｄｕｃｔｉｏｎＬｅｔΩｃｏｎｔａｉｎｉｎｇｚｅｒｏｐｏｉｎｔｂｅａｓｉｍｐｌｙ＿ｃｏｎｎｅｃｔｅｄｂｏｕｎｄｅｄｏｐｅｎｓｅｔｏｆＲ2 ｗｉｔｈｓｍｏｏｔｈｂｏｕｎｄａｒｙΓａｎｄｌｅｔΩ′ｄｅｎｏｔｅｔｈｅｃｏｍｐｌｅｍｅｎｔｏｆΩ ∪Γ .ＴｈｅｅｘｔｅｒｉｏｒｎｏｎｓｔａｔｉｏｎａｒｙＮａｖｉｅｒ＿ＳｔｏｋｅｓｐｒｏｂｌｅｍｆｏｒａｆｌｕｉｄｏｃｃｕｐｙｉｎｇΩ′ｃｏｎｓｉｓｔｓｉｎｆｉｎｄｉｎｇｔｈｅｖｅｌｏｃｉｔｙ ｕ(ｘ,ｔ)ｏｆｔｈｅｆｌｕｉｄａｎｄｉｔｓｐｒｅｓｓｕｒｅ ｐ(ｘ ,…     

A high-order characteristics/finite element method for the incompressible Navier-Stokes equations

In this paper we consider a discretization of the incompressible Navier-Stokes equations involving a second-order time scheme based on the characteristics method and a spatial discretization of finite element type. Theoretical and numerical analyses are detailed and we obtain stability results abnd optimal eror estimates on the velocity and pressure under a time step restriction less stringent than the standard Courant-Freidrichs-Levy condition. Finally, some numerical results obtained wiht the code N3S are shown which justify the interest of this scheme and its advantages with respect to an analogous first-order time scheme. © 1997 John Wiley & Sons, Ltd.     

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q ₁ ^rot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H ¹-norm and the pressure in the L ²-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.     

An extension of the SIMPLE based discontinuous Galerkin solver to unsteady incompressible flows

In this paper, we present a SIMPLE based algorithm in the context of the discontinuous Galerkin method for unsteady incompressible flows. Time discretization is done fully implicit using backward differentiation formulae (BDF) of varying order from 1 to 4. We show that the original equation for the pressure correction can be modified by using an equivalent operator stemming from the symmetric interior penalty (SIP) method leading to a reduced stencil size. To assess the accuracy as well as the stability and the performance of the scheme, three different test cases are carried out: the Taylor vortex flow, the Orr‐Sommerfeld stability problem for plane Poiseuille flow and the flow past a square cylinder. (1) Simulating the Taylor vortex flow, we verify the temporal accuracy for the different BDF schemes. Using the mixed‐order formulation, a spatial convergence study yields convergence rates of k + 1 and k in the L²‐norm for velocity and pressure, respectively. For the equal‐order formulation, we obtain approximately the same convergence rates, while the absolute error is smaller. (2) The stability of our method is examined by simulating the Orr–Sommerfeld stability problem. Using the mixed‐order formulation and adjusting the penalty parameter of the symmetric interior penalty method for the discretization of the viscous part, we can demonstrate the long‐term stability of the algorithm. Using pressure stabilization the equal‐order formulation is stable without changing the penalty parameter. (3) Finally, the results for the flow past a square cylinder show excellent agreement with numerical reference solutions as well as experiments. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessment of shock capturing schemes for discontinuous Galerkin method

This paper carries out systematical investigations on the performance of several typical shock-capturing schemes for the discontinuous Galerkin （DG） method, including the total variation bounded （TVB） limiter and three artificial diffusivity schemes （the basis function-based （BF） scheme, the face residual-based （FR） scheme, and the element residual-based （ER） scheme）. Shock-dominated flows （the Sod problem, the Shu- Osher problem, the double Mach reflection problem, and the transonic NACA0012 flow） are considered, addressing the issues of accuracy, non-oscillatory property, dependence on user-specified constants, resolution of discontinuities, and capability for steady solutions. Numerical results indicate that the TVB limiter is more efficient and robust, while the artificial diffusivity schemes are able to preserve small-scale flow structures better. In high order cases, the artificial diffusivity schemes have demonstrated superior performance over the TVB limiter.     

A nonlinear Galerkin/Petrov-least squares mixed element method for the stationary Navier-Stokes equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｏｎｌｉｎｅａｒＧａｌｅｒｋｉｎｍｅｔｈｏｄｉｓａｍｕｌｔｉ＿ｌｅｖｅｌｓｃｈｅｍｅｔｏｆｉｎｄｔｈｅａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｆｏｒｔｈｅｄｉｓｓｉｐａｔｉｖｅＰＤＥ .ＴｈｉｓｍｅｔｈｏｄｈａｓｆｉｒｓｔｍａｉｎｌｙｂｅｅｎａｄｄｒｅｓｓｅｄｂｙＦｏｉａｓ＿Ｍａｎｌｅｙ＿Ｔｅｍａｍ[1],Ｍａｒｉｏｎ＿Ｔｅｍａｍ[2 ],Ｆｏｉａｓ＿Ｊｏｌｌｙ＿Ｋｅｖｒｅｋｉｄｉｓ＿Ｔｉｔｉ[3]ａｎｄＤｅｖｕｌｄｅｒ＿Ｍａｒｉｏｎ＿Ｔｉｔｉ[4 ]ｉｎｔｈｅｃａｓｅｏｆｓｐｅｃｔ…