ADER Schemes for Nonlinear Systems of Stiff Advection�CDiffusion�CReaction Equations

In this article we extend the high order ADER finite volume schemes introduced for stiff hyperbolic balance laws by Dumbser, Enaux and Toro (J. Comput. Phys. 227:3971?C4001, 2008) to nonlinear systems of advection?Cdiffusion?Creaction equations with stiff algebraic source terms. We derive a new efficient formulation of the local space-time discontinuous Galerkin predictor using a nodal approach whose interpolation points are tensor-products of Gauss?CLegendre quadrature points. Furthermore, we propose a new simple and efficient strategy to compute the initial guess of the locally implicit space-time DG scheme: the Gauss?CLegendre points are initialized sequentially in time by a second order accurate MUSCL-type approach for the flux term combined with a Crank?CNicholson method for the stiff source terms. We provide numerical evidence that when starting with this initial guess, the final iterative scheme for the solution of the nonlinear algebraic equations of the local space-time DG predictor method becomes more efficient. We apply our new numerical method to some systems of advection?Cdiffusion?Creaction equations with particular emphasis on the asymptotic preserving property for linear model systems and the compressible Navier?CStokes equations with chemical reactions.     

Hermite WENO Schemes and Their Application as Limiters for Runge-Kutta Discontinuous Galerkin Method,III: Unstructured Meshes

In [J. Comput. Phys. 193:115–135, 2004] and [Comput. Fluids 34:642–663, 2005], Qiu and Shu developed a class of high order weighted essentially non-oscillatory (WENO) schemes based on Hermite polynomials, termed HWENO (Hermite WENO) schemes, for solving nonlinear hyperbolic conservation law systems, and applied them as limiters for the Runge-Kutta discontinuous Galerkin (RKDG) methods on structured meshes. In this continuation paper, we extend the method to solve two dimensional problems on unstructured meshes. The emphasis is again on the application of such HWENO finite volume methodology as limiters for RKDG methods to maintain compactness of RKDG methods. Numerical experiments for two dimensional Burgers’ equation and Euler equations of compressible gas dynamics are presented to show the effectiveness of these methods. The research was partially supported by the European project ADIGMA on the development of innovative solution algorithms for aerodynamic simulations, NSFC grant 10671091 and JSNSF BK2006511.     

A Simple Extension of the Osher Riemann Solver to?Non-conservative Hyperbolic Systems

We propose a simple extension of the well-known Riemann solver of Osher and Solomon (Math. Comput. 38:339?C374, 1982) to a certain class of hyperbolic systems in non-conservative form, in particular to shallow-water-type and multi-phase flow models. To this end we apply the formalism of path-conservative schemes introduced by Parés (SIAM J. Numer. Anal. 44:300?C321, 2006) and Castro et al. (Math. Comput. 75:1103?C1134, 2006). For the sake of generality and simplicity, we suggest to compute the inherent path integral numerically using a Gaussian quadrature rule of sufficient accuracy. Published path-conservative schemes to date are based on either the Roe upwind method or on centered approaches. In comparison to these, the proposed new path-conservative Osher-type scheme has several advantages. First, it does not need an entropy fix, in contrast to Roe-type path-conservative schemes. Second, our proposed non-conservative Osher scheme is very simple to implement and nonetheless constitutes a complete Riemann solver in the sense that it attributes a different numerical viscosity to each characteristic field present in the relevant Riemann problem; this is in contrast to centered methods or incomplete Riemann solvers that usually neglect intermediate characteristic fields, hence leading to excessive numerical diffusion. Finally, the interface jump term is differentiable with respect to its arguments, which is useful for steady-state computations in implicit schemes. We also indicate how to extend the method to general unstructured meshes in multiple space dimensions. We show applications of the first order version of the proposed path-conservative Osher-type scheme to the shallow water equations with variable bottom topography and to the two-fluid debris flow model of Pitman & Le. Then, we apply the higher-order multi-dimensional version of the method to the Baer?CNunziato model of compressible multi-phase flow. We also clearly emphasize the limitations of our approach in a special chapter at the end of this article.     

A simple construction of very high order non-oscillatory compact schemes on unstructured meshes

In [Abgrall R, Roe PL. High order fluctuation schemes on triangular meshes. J Sci Comput 2003;19(1-3):3-36] have been constructed very high order residual distribution schemes for scalar problems. They were using triangle unstructured meshes. However, the construction was quite involved and was not very flexible. Here, following [Abgrall R. Essentially non-oscillatory residual distribution schemes for hyperbolic problems. J Comput Phys 2006;214(2):773-808], we develop a systematic way of constructing very high order non-oscillatory schemes for such meshes. Applications to scalar and systems problems are given.     

Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier-Stokes equations

In this paper, we propose a new unified family of arbitrary high order accurate explicit one-step finite volume and discontinuous Galerkin schemes on unstructured triangular and tetrahedral meshes for the solution of the compressible Navier-Stokes equations. This new family of numerical methods has first been proposed in [16] for purely hyperbolic systems and has been called P_NP_M schemes, where N indicates the polynomial degree of the test functions and M is the degree of the polynomials used for flux and source computation. A particular feature of the general P_NP_M schemes is that they contain classical high order accurate finite volume schemes (N=0) as well as standard discontinuous Galerkin methods (M=N) just as special cases, which therefore allows for a direct efficiency comparison.In the application section of this paper we first show numerical convergence results on unstructured meshes obtained for the compressible Navier-Stokes equations with Sutherland’s viscosity law, comparing all third to sixth order accurate P_NP_M schemes with each other. In order to validate the method also in practice we show several classical steady and unsteady CFD applications, such as the laminar boundary layer flow over a flat plate at high Reynolds numbers, flow past a NACA0012 airfoil, the unsteady flows past a circular cylinder and a sphere, the unsteady flows of a compressible mixing layer in two space dimensions and finally we also show applications to supersonic flows with shock Mach numbers up to M_s=10.     

Discontinuous Galerkin error estimation for hyperbolic problems on unstructured triangular meshes

We extend the error analysis of Adjerid and Baccouch [1], [2] for the discontinuous Galerkin discretization error to variable-coefficient linear hyperbolic problems as well as nonlinear hyperbolic problems on unstructured meshes. We further extend this analysis to transient hyperbolic problems and prove that the local superconvergence results presented in [1] still hold for both steady and transient variable-coefficient linear and nonlinear problems. This local error analysis allows us to construct asymptotically correct a posteriori error estimates by solving local hyperbolic problems with no boundary conditions on each element of general unstructured meshes. We illustrate the superconvergence and the efficiency of our a posteriori error estimates by showing computational results for several linear and nonlinear numerical examples.     

Maximum-Principle-Satisfying and Positivity-Preserving High Order Discontinuous Galerkin Schemes for?Conservation Laws on Triangular Meshes

In Zhang and Shu (J. Comput. Phys. 229:3091–3120, 2010), two of the authors constructed uniformly high order accurate finite volume and discontinuous Galerkin (DG) schemes satisfying a strict maximum principle for scalar conservation laws on rectangular meshes. The technique is generalized to positivity preserving (of density and pressure) high order DG or finite volume schemes for compressible Euler equations in Zhang and Shu (J. Comput. Phys. 229:8918–8934, 2010). The extension of these schemes to triangular meshes is conceptually plausible but highly nontrivial. In this paper, we first introduce a special quadrature rule which is exact for two-variable polynomials over a triangle of a given degree and satisfy a few other conditions, by which we can construct high order maximum principle satisfying finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or DG method solving two dimensional scalar conservation laws on triangular meshes. The same method can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. We also obtain positivity preserving (for density and pressure) high order DG or finite volume schemes solving compressible Euler equations on triangular meshes. Numerical tests for the third order Runge-Kutta DG (RKDG) method on unstructured meshes are reported.     

Selecting the Numerical Flux in Discontinuous Galerkin Methods for Diffusion Problems

In this paper we present numerical investigations of four different formulations of the discontinuous Galerkin method for diffusion problems. Our focus is to determine, through numerical experimentation, practical guidelines as to which numerical flux choice should be used when applying discontinuous Galerkin methods to such problems. We examine first an inconsistent and weakly unstable scheme analyzed in Zhang and Shu, Math. Models Meth. Appl. Sci. (M ³ AS) 13, 395–413 (2003), and then proceed to examine three consistent and stable schemes: the Bassi–Rebay scheme (J. Comput. Phys. 131, 267 (1997)), the local discontinuous Galerkin scheme (SIAM J. Numer. Anal. 35, 2440–2463 (1998)) and the Baumann–Oden scheme (Comput. Math. Appl. Mech. Eng. 175, 311–341 (1999)). For an one-dimensional model problem, we examine the stencil width, h-convergence properties, p-convergence properties, eigenspectra and system conditioning when different flux choices are applied. We also examine the ramifications of adding stabilization to these schemes. We conclude by providing the pros and cons of the different flux choices based upon our numerical experiments.This revised version was published online in July 2005 with corrected volume and issue numbers.     

WENO Scheme with Subcell Resolution for Computing Nonconservative Euler Equations with Applications to One-Dimensional Compressible Two-Medium Flows

High order path-conservative schemes have been developed for solving nonconservative hyperbolic systems in (Parés, SIAM J.?Numer. Anal. 44:300?C321, 2006; Castro et al., Math. Comput. 75:1103?C1134, 2006; J.?Sci. Comput. 39:67?C114, 2009). Recently, it has been observed in (Abgrall and Karni, J.?Comput. Phys. 229:2759?C2763, 2010) that this approach may have some computational issues and shortcomings. In this paper, a modification to the high order path-conservative scheme in (Castro et al., Math. Comput. 75:1103?C1134, 2006) is proposed to improve its computational performance and to overcome some of the shortcomings. This modification is based on the high order finite volume WENO scheme with subcell resolution and it uses an exact Riemann solver to catch the right paths at the discontinuities. An application to one-dimensional compressible two-medium flows of nonconservative or primitive Euler equations is carried out to show the effectiveness of this new approach.     

Simulation of a viscous compressible flow past a circular cylinder with high-order discontinuous Galerkin methods

This paper is devoted to the simulation of viscous compressible flows with high-order accurate discontinuous Galerkin methods on bidimensional unstructured meshes. The formulation for the solution of the Navier–Stokes equations is due to Oden et al. [An hp-adaptive discontinuous finite element method for computational fluid dynamics. PhD thesis, The University of Texas at Austin, 1997; J Comput Phys 1998;146:491–519]. It involves a weak imposition of continuity conditions on the state variables and on fluxes across interelement boundaries. It does not make use of any auxiliary variables and then the cost for the implementation is reasonable. The method is coupled with a limiting procedure recently developed by the authors to suppress oscillations near large gradients. The limiter is totally free of problem dependence and maintains the convergence order for errors measured in the L¹-norm. This paper presents detailed numerical results of a viscous compressible flow past a circular cylinder at a Reynolds number of 100 for the cases of subsonic and supersonic regimes. The proposed simulations suggest that the method is very robust and is able to produce very accurate results on unstructured meshes.     

A discontinuous Galerkin solver for Boltzmann–Poisson systems in nano devices

In this paper, we present results of a discontinuous Galerkin (DG) scheme applied to deterministic computations of the transients for the Boltzmann–Poisson system describing electron transport in semiconductor devices. The collisional term models optical-phonon interactions which become dominant under strong energetic conditions corresponding to nano-scale active regions under applied bias. The proposed numerical technique is a finite element method using discontinuous piecewise polynomials as basis functions on unstructured meshes. It is applied to simulate hot electron transport in bulk silicon, in a silicon n⁺–n–n⁺ diode and in a double gated 12 nm MOSFET. Additionally, the obtained results are compared to those of a high order WENO scheme simulation and DSMC (Discrete Simulation Monte Carlo) solvers.     

Direct Discontinuous Galerkin Method and Its Variations for Second Order Elliptic Equations

In this paper, we study direct discontinuous Galerkin method (Liu and Yan in SIAM J Numer Anal 47(1):475–698, 2009) and its variations (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010; Vidden and Yan in J Comput Math 31(6):638–662, 2013; Yan in J Sci Comput 54(2–3):663–683, 2013) for 2nd order elliptic problems. A priori error estimate under energy norm is established for all four methods. Optimal error estimate under \(L^2\) norm is obtained for DDG method with interface correction (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010) and symmetric DDG method (Vidden and Yan in J Comput Math 31(6):638–662, 2013). A series of numerical examples are carried out to illustrate the accuracy and capability of the schemes. Numerically we obtain optimal \((k+1)\)th order convergence for DDG method with interface correction and symmetric DDG method on nonuniform and unstructured triangular meshes. An interface problem with discontinuous diffusion coefficients is investigated and optimal \((k+1)\)th order accuracy is obtained. Peak solutions with sharp transitions are captured well. Highly oscillatory wave solutions of Helmholz equation are well resolved.     

An estimate on the convergence of MKZ–Bézier operators

The pointwise approximation properties of the MKZ–Bézier operators M_n,α(f,x) for α≥1 have been studied in [X.M. Zeng, Rates of approximation of bounded variation functions by two generalized Meyer–König–Zeller type operators, Comput. Math. Appl. 39 (2000) 1–13]. The aim of this paper is to study the pointwise approximation of the operators M_n,α(f,x) for the other case 0<α<1. By means of some new estimate techniques and a result of Guo and Qi [S. Guo, Q. Qi, The moments for the Meyer–König and Zeller operators, Appl. Math. Lett. 20 (2007) 719–722], we establish an estimate formula on the rate of convergence of the operators M_n,α(f,x) for the case 0<α<1.     

Local-Structure-Preserving Discontinuous Galerkin Methods with Lax-Wendroff Type Time Discretizations for Hamilton-Jacobi Equations

In this paper, a family of high order numerical methods are designed to solve the Hamilton-Jacobi equation for the viscosity solution. In particular, the methods start with a hyperbolic conservation law system closely related to the Hamilton-Jacobi equation. The compact one-step one-stage Lax-Wendroff type time discretization is then applied together with the local-structure-preserving discontinuous Galerkin spatial discretization. The resulting methods have lower computational complexity and memory usage on both structured and unstructured meshes compared with some standard numerical methods, while they are capable of capturing the viscosity solutions of Hamilton-Jacobi equations accurately and reliably. A collection of numerical experiments is presented to illustrate the performance of the methods.     

A Conservative Method for the Simulation of?the?Isothermal Euler System with the van-der-Waals Equation of State

In this article, we are interested in the simulation of phase transition in compressible flows, with the isothermal Euler system, closed by the van-der-Waals model. We formulate the problem as an hyperbolic system, with a source term located at the interface between liquid and vapour. The numerical scheme is based on (Abgrall and Saurel, J. Comput. Phys. 186(2):361?C396, 2003; Le Métayer et al., J. Comput. Phys. 205(2):567?C610, 2005). Compared with previous discretizations of the van-der-Waals system, the novelty of this algorithm is that it is fully conservative. Its Godunov-type formulation allows an easy implementation on multi-dimensional unstructured meshes.     

Smoothness-Increasing Accuracy-Conserving (SIAC) Filters for Discontinuous Galerkin Solutions: Application to Structured Tetrahedral Meshes

In this paper, we attempt to address the potential usefulness of smoothness-increasing accuracy-conserving (SIAC) filters when applied to real-world simulations. SIAC filters as a class of post-processors were initially developed in Bramble and Schatz (Math Comput 31:94, 1977) and later applied to discontinuous Galerkin (DG) solutions of linear hyperbolic partial differential equations by Cockburn et al. (Math Comput 72:577, 2003), and are successful in raising the order of accuracy from $k+1$ to $2k+1$ in the $L^2$ —norm when applied to a locally translation-invariant mesh. While there have been several attempts to demonstrate the usefulness of this filtering technique to nontrivial mesh structures (Curtis et al. in SIAM J Sci Comput 30(1):272, 2007; Mirzaee et al. in SIAM J Numer Anal 49:1899, 2011; King et al. in J Sci Comput, 2012), the application of the SIAC filter never exceeded beyond two-space dimensions. As tetrahedral meshes are often the type considered in more realistic simulations, we contribute to the class of SIAC post-processors by demonstrating the effectiveness of SIAC filtering when applied to structured tetrahedral meshes. These types of meshes are generated by tetrahedralizing uniform hexahedra and therefore, while maintaining the structured nature of a hexahedral mesh, they exhibit an unstructured tessellation within each hexahedral element. Moreover, we address the computationally intensive task of performing numerical integrations when one considers tetrahedral elements for SIAC filtering and provide guidelines on how to ameliorate these challenges through the use of more general cubature rules. We consider two examples of a hyperbolic equation and confirm the usefulness of SIAC filters in obtaining the superconvergence accuracy of $2k+1$ when applied to structured tetrahedral meshes. Additionally, the DG methodology merely requires weak constraints on the fluxes between elements. As SIAC filters improve this weak continuity to $\mathcal{C }^{k-1}$ —continuity at the element interfaces, we provide results that show how post-processing is useful in extracting smooth isosurfaces of DG fields.     

A discontinuous finite difference streamline diffusion method for time-dependent hyperbolic problems

In this article, a new finite element method, discontinuous finite difference streamline diffusion method (DFDSD), is constructed and studied for first-order linear hyperbolic problems. This method combines the benefit of the discontinuous Galerkin method and the streamline diffusion finite element method. Two fully discrete DFDSD schemes (Euler DFDSD and Crank–Nicolson (CN) DFDSD) are constructed by making use of the difference discrete method for time variables and the discontinuous streamline diffusion method for space variables. The stability and optimal L² norm error estimates are established for the constructed schemes. This method makes contributions to the discontinuous methods. Finally, a numerical example is provided to show the benefit of high efficiency and simple implementation of the schemes.     

Solutions of Multi-dimensional Hyperbolic Systems of Conservation Laws by Square Entropy Condition Satisfying Discontinuous Galerkin Method

In this paper, we study formally high-order accurate discontinuous Galerkin methods on general arbitrary grid for multi-dimensional hyperbolic systems of conservation laws [Cockburn, B., and Shu, C.-W. (1989, Math. Comput. 52, 411–435, 1998, J. Comput. Phys. 141, 199–224); Cockburn et al. (1989, J. Comput. Phys. 84, 90–113; 1990, Math. Comput. 54, 545–581). We extend the notion of E-flux [Osher (1985) SIAM J. Numer. Anal. 22, 947–961] from scalar to system, and found that after flux splitting upwind flux [Cockburn et al. (1989) J. Comput. Phys. 84, 90–113] is a Riemann solver free E-flux for systems. Therefore, we are able to show that the discontinuous Galerkin methods satisfy a cell entropy inequality for square entropy (in semidiscrete sense) if the multi-dimensional systems are symmetric. Similar result [Jiang and Shu (1994) Math. Comput. 62, 531–538] was obtained for scalar equations in multi-dimensions. We also developed a second-order finite difference version of the discontinuous Galerkin methods. Numerical experiments have been obtained with excellent results.      

ADER: Arbitrary High Order Godunov Approach

This paper concerns the construction of non-oscillatory schemes of very high order of accuracy in space and time, to solve non-linear hyperbolic conservation laws. The schemes result from extending the ADER approach, which is related to the ENO/WENO methodology. Our schemes are conservative, one-step, explicit and fully discrete, requiring only the computation of the inter-cell fluxes to advance the solution by a full time step; the schemes have optimal stability condition. To compute the intercell flux in one space dimension we solve a generalised Riemann problem by reducing it to the solution a sequence of m conventional Riemann problems for the kth spatial derivatives of the solution, with k=0, 1,..., m–1, where m is arbitrary and is the order of the accuracy of the resulting scheme. We provide numerical examples using schemes of up to fifth order of accuracy in both time and space.