Discontinuous Galerkin spectral element lattice Boltzmann method on triangular element

Discontinuous Galerkin spectral element method is used to solve the lattice Boltzmann equation (LBE) in the discrete velocity space. The triangular elements are adopted because of their flexibility to deal with complex geometries. The flow past a circular cylinder is simulated by the proposed scheme. The results are consistent with those obtained from the previous numerical methods and experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

Connections between the discontinuous Galerkin method and high‐order flux reconstruction schemes

With high‐order methods becoming more widely adopted throughout the field of computational fluid dynamics, the development of new computationally efficient algorithms has increased tremendously in recent years. One of the most recent methods to be developed is the flux reconstruction approach, which allows various well‐known high‐order schemes to be cast within a single unifying framework. Whilst a connection between flux reconstruction and the more widely adopted discontinuous Galerkin method has been established elsewhere, it still remains to fully investigate the explicit connections between the many popular variants of the discontinuous Galerkin method and the flux reconstruction approach. In this work, we closely examine the connections between three nodal versions of tensor‐product discontinuous Galerkin spectral element approximations and two types of flux reconstruction schemes for solving systems of conservation laws on quadrilateral meshes. The different types of discontinuous Galerkin approximations arise from the choice of the solution nodes of the Lagrange basis representing the solution and from the quadrature approximation used to integrate the mass matrix and the other terms of the discretization. By considering both linear and nonlinear advection equations on a regular grid, we examine the mathematical properties that connect these discretizations. These arguments are further confirmed by the results of an empirical numerical study. Copyright © 2014 John Wiley & Sons, Ltd.     

A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications

We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 10⁵?10⁶. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.     

A triangular spectral/hp discontinuous Galerkin method for modelling 2D shallow water equations

We present a spectral/hp element discontinuous Galerkin model for simulating shallow water flows on unstructured triangular meshes. The model uses an orthogonal modal expansion basis of arbitrary order for the spatial discretization and a third‐order Runge–Kutta scheme to advance in time. The local elements are coupled together by numerical fluxes, evaluated using the HLLC Riemann solver. We apply the model to test cases involving smooth flows and demonstrate the exponentially fast convergence with regard to polynomial order. We also illustrate that even for results of ‘engineering accuracy’ the computational efficiency increases with increasing order of the model and time of integration. The model is found to be robust in the presence of shocks where Gibbs oscillations can be suppressed by slope limiting. Copyright 2004 John Wiley & Sons, Ltd.     

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.     

High‐order discontinuous Galerkin spectral element methods for transitional and turbulent flow simulations

In this paper, we investigate the accuracy and efficiency of discontinuous Galerkin spectral method simulations of under‐resolved transitional and turbulent flows at moderate Reynolds numbers, where the accurate prediction of closely coupled laminar regions, transition and developed turbulence presents a great challenge to large eddy simulation modelling. We take full advantage of the low numerical errors and associated superior scale resolving capabilities of high‐order spectral methods by using high‐order ansatz functions up to 12th order. We employ polynomial de‐aliasing techniques to prevent instabilities arising from inexact quadrature of nonlinearities. Without the need for any additional filtering, explicit or implicit modelling, or artificial dissipation, our high‐order schemes capture the turbulent flow at the considered Reynolds number range very well. Three classical large eddy simulation benchmark problems are considered: a circular cylinder flow at Re_D=3900, a confined periodic hill flow at Re_h=2800 and the transitional flow over a SD7003 airfoil at Re_c=60,000. For all computations, the total number of degrees of freedom used for the discontinuous Galerkin spectral method simulations is chosen to be equal or considerably less than the reported data in literature. In all three cases, we achieve an equal or better match to direct numerical simulation results, compared with other schemes of lower order with explicitly or implicitly added subgrid scale models. Copyright © 2014 John Wiley & Sons, Ltd.     

固体非傅立叶温度场的时域间断Galerkin有限元法

运用时域间断Galerkin有限元法[1],对高频非傅立叶热波动问题[2-3]进行分析。其主要特点是:取温度及温度的时间导数为基本未知量,对其分别进行3次Hermite插值和线性插值。在保证节点温度自动保持连续的基础上,温度的时间导数在离散时域存在间断。数值结果表明所提出的方法能够滤掉虚假的数值震荡,能够良好地模拟固体中的非傅立叶热波动行为。     

Coupling of finite element method and discontinuous Galerkin method to simulate viscoelastic flows

In this paper, a numerical method, which is about the coupling of continuous and discontinuous Galerkin method based on the splitting scheme, is presented for the calculation of viscoelastic flows of the Oldroyd‐B fluid. The momentum equation is discretized in time by using the Adams‐Bashforth second‐order algorithm, and then decoupled via the splitting approach. Considering the Oldroyd‐B constitutive equation, the second‐order Runge‐Kutta approach is selected to complete the temporal discretization. As for the spatial discretizations, the fundamental purpose is to make the best of finite element method (FEM) and discontinuous Galerkin (DG) method to handle different types of equations. Specifically speaking, for the subequations, FEM is chosen to treat the Poisson and Helmholtz equations, and DG is employed to deal with the nonlinear convective term. In addition, because of the hyperbolic nature, DG is also utilized to discretize the Oldroyd‐B constitutive equation spatially. This coupled method avoids resorting to extra stabilization technique occurred in standard FEM framework even for moderately high values of Weissenberg number and also reduces the complexity compared with unified DG scheme. The Oldroyd‐B model is applied to investigate several typical and challenging benchmarks, such as the 4:1 planar contraction flow and the lid‐driven cavity flow, with a wide range of Weissenberg number to illustrate the feasibility, robustness, and validity of our coupled method.     

A moving discontinuous Galerkin finite element method for flows with interfaces

A moving discontinuous Galerkin finite element method with interface condition enforcement is formulated for flows with discontinuous interfaces. The underlying weak formulation enforces the interface condition separately from the conservation law, so that the residual only vanishes upon satisfaction of both. In this formulation, the discrete grid geometry is treated as a variable, so that, in contrast to the standard discontinuous Galerkin method, this method has both the means to detect interfaces, via interface condition enforcement, and to satisfy, via grid movement, the conservation law and its associated interface condition. The method therefore directly fits interfaces, including shocks, preserving a high-order representation up to the interface without requiring shock capturing or an upwind numerical flux to achieve stability. It can be generalized to flows with a priori unknown interfaces with nontrivial topology and curved interface geometry as well as to an arbitrary number of spatial dimensions. Unsteady flows are represented in a manner similar to steady flows using a space-time formulation. In addition to computing flows with interfaces, the method can represent point singularities in a flow field by degenerating cuboid elements. In general, the method works in conjunction with standard local grid operations, including edge collapse, to ensure that degenerate cells are removed. Test cases are presented for up to three-dimensional flows that provide an initial assessment of the stability and accuracy of the method.     

A high‐order accurate discontinuous Galerkin finite element method for laminar low Mach number flows

In this paper we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of laminar flow simulations at low Mach numbers using an implicit scheme. The algorithm is based on the flux preconditioning approach, which modifies only the dissipative terms of the numerical flux. This formulation is quite simple to implement in existing implicit DG codes, it overcomes the time‐stepping restrictions of explicit multistage algorithms, is consistent in time and thus applicable to unsteady flows. The performance of the method is demonstrated by solving the flow around a NACA0012 airfoil and on a flat plate, at different low Mach numbers using various degrees of polynomial approximations. Computations with and without flux preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers. The time accurate solution of unsteady flow is also demonstrated by solving the vortex shedding behind a circular cylinder at the Reynolds number of 100. Copyright © 2012 John Wiley & Sons, Ltd.     

A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

基于时域间断有限元方法的生物组织非傅立叶热行为数值分析

在热传导分析中,当热流与温度梯度存在时间延迟时,需采用非傅立叶热传导模型进行分析。生物组织具有较强的热松弛时间系数,承受激光、微波及烧烫等作用时,其呈现出较强的非傅立叶行为。本文对脉冲热源作用下生物组织的非傅立叶热传导进行研究,针对强脉冲引起的温度场在空间域的高梯度变化、波阵面的间断行为以及通用传统时域数值方法会带来虚假数值振荡的特点,提出采用所发展的时域间断Galerkin有限元法（DG-FEM ）进行求解计算。对多种脉冲热源作用下的非傅立叶热传导过程进行数值模拟,通过考量强脉冲作用下温度场分布和热致生物组织损伤行为的影响,表明了本文所发展的DGFEM 能够有效、准确地描述温度场空间分布和热传导过程以及非傅立叶行为下的生物热损伤更为明显,在生物组织热行为分析中应该受到重视。     

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.     

From h to p efficiently: optimal implementation strategies for explicit time‐dependent problems using the spectral/hp element method

We investigate the relative performance of a second‐order Adams–Bashforth scheme and second‐order and fourth‐order Runge–Kutta schemes when time stepping a 2D linear advection problem discretised using a spectral/hp element technique for a range of different mesh sizes and polynomial orders. Numerical experiments explore the effects of short (two wavelengths) and long (32 wavelengths) time integration for sets of uniform and non‐uniform meshes. The choice of time‐integration scheme and discretisation together fixes a CFL limit that imposes a restriction on the maximum time step, which can be taken to ensure numerical stability. The number of steps, together with the order of the scheme, affects not only the runtime but also the accuracy of the solution. Through numerical experiments, we systematically highlight the relative effects of spatial resolution and choice of time integration on performance and provide general guidelines on how best to achieve the minimal execution time in order to obtain a prescribed solution accuracy. The significant role played by higher polynomial orders in reducing CPU time while preserving accuracy becomes more evident, especially for uniform meshes, compared with what has been typically considered when studying this type of problem.© 2014. The Authors. International Journal for Numerical Methods in Fluids published by John Wiley & Sons, Ltd.     

On the discontinuous Galerkin method for the numerical solution of the Navier–Stokes equations

The paper deals with the use of the discontinuous Galerkin finite element method (DGFEM) for the numerical solution of viscous compressible flows. We start with a scalar convection–diffusion equation and present a discretization with the aid of the non‐symmetric variant of DGFEM with interior and boundary penalty terms. We also mention some theoretical results. Then we extend the scheme to the system of the Navier–Stokes equations and discuss the treatment of stabilization terms. Several numerical examples are presented. Copyright © 2004 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.     

A numerical investigation of a spectral‐type nodal collocation discontinuous Galerkin approximation of the Euler and Navier–Stokes equations

In this paper, we focus on the applicability of spectral‐type collocation discontinuous Galerkin methods to the steady state numerical solution of the inviscid and viscous Navier–Stokes equations on meshes consisting of curved quadrilateral elements. The solution is approximated with piecewise Lagrange polynomials based on both Legendre–Gauss and Legendre–Gauss–Lobatto interpolation nodes. For the sake of computational efficiency, the interpolation nodes can be used also as quadrature points. In this case, however, the effect of the nonlinearities in the equations and/or curved elements leads to aliasing and/or commutation errors that may result in inaccurate or unstable computations. By a thorough numerical testing on a set of well known test cases available in the literature, it is here shown that the two sets of nodes behave very differently, with a clear advantage of the Legendre–Gauss nodes, which always displayed an accurate and robust behaviour in all the test cases considered.Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

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.