Search

This study addresses the existence and uniqueness of solutions to the Hoogendoorn–Bovy (HB) pedestrian flow model, which describes the dynamic user-optimal pedestrian flow assignment problem in continuous space and time. The HB model consists of a forward conservation law (CL) equation that governs density and a backward Hamilton–Jacobi–Bellman (HJB) equation that contains a maximum admissible speed constraint (MASC), in which the flow direction is determined by the path-choice strategy. The existence and uniqueness of solutions are significantly more difficult to determine when the HJB equation contains an MASC; however, we prove that the HB model can be formulated as a forward CL equation and backward Hamilton–Jacobi (HJ) equation in which the MASC is non-binding if suitable model parameters are selected. This model is formulated as a fixed-point problem upon the simultaneous satisfaction of both equations. To verify the existence and uniqueness results, we first demonstrate the existence and uniqueness of the solutions to the CL and HJ equations, and then show that the coupled HB model is well-posed and has a unique solution. A numerical example is presented to illustrate the properties of the HB model. [ABSTRACT FROM AUTHOR]

In this paper we propose and analyze finite element discontinuous Galerkin methods for the one- and two-dimensional stochastic Maxwell equations with multiplicative noise. The discrete energy law of the semi-discrete DG methods were studied. Optimal error estimate of the semi-discrete method is obtained for the one-dimensional case, and the two-dimensional case on both rectangular meshes and triangular meshes under certain mesh assumptions. Strong Taylor 2.0 scheme is used as the temporal discretization. Both one- and two-dimensional numerical results are presented to validate the theoretical analysis results. [ABSTRACT FROM AUTHOR]

A second order method and a limited variation WENO numerical scheme are proposed for the Scanning Hysteresis Model (SHM) which describes continuum flows with hysteretic looping fluxes with local memory. The SHM is a nonconservative hyperbolic system with discontinuous coefficients. The limited variation WENO method is fifth order accurate in spatial coordinate and third order accurate in time steps in smooth and monotone regions. The proposed high resolution methods produce the desired solutions more efficiently than the first order upwinding method. These methods are total variation diminishing for the flux, and monotonicity-preserving for the flux and the hysteresis state variable. The flux and the hysteresis state variable produced by these two methods satisfy maximum principle. These two methods also have invariant regions. [ABSTRACT FROM AUTHOR]

In this paper, we are concerned about a time-domain carpet cloak model, which was originally derived in our previous work Li et al. (SIAM J. Appl. Math., 74(4), pp. 1136–1151, 2014). Some finite element schemes have been developed for this model and used to simulate the cloaking phenomenon in Li et al. (SIAM J. Appl. Math., 74(4), pp. 1136–1151, 2014) and Li et al. (Methods Appl. Math., 19(2), pp. 359–378, 2019). However, numerical stabilities for those proposed explicit schemes are only proved under the time step constrain τ = O(h2), which is impractical and too restricted. To overcome this disadvantage, we propose two new finite element schemes for solving this carpet cloak model: one is the implicit Crank-Nicolson (CN) scheme, and another one is the explicit leap-frog (LF) scheme. Inspired by a totally new energy developed for the continuous model, we prove the unconditional stability for the CN scheme and conditional stability for the LF scheme under the usual CFL constraint τ = O(h). Both numerical stabilities inherit the exact form as the continuous stability. Optimal error estimate is also established for the LF scheme. Finally, numerical results using the LF scheme are presented to support our analysis and demonstrate the cloaking phenomenon. [ABSTRACT FROM AUTHOR]

The three-temperature (3-T) radiation hydrodynamics (RH) equations play an important role in the high-energy-density-physics fields, such as astrophysics and inertial confinement fusion (ICF). It describes the interaction between radiation and high-energy-density plasmas including electron and ion in the assumption that radiation, electron and ion are in their own equilibrium state, which means they can be characterized by their own temperatures. The 3-T RH system consists of the density, momentum and three internal energy (electron, ion and radiation) equations. In this paper, we propose a high order conservative finite difference weighted essentially non-oscillatory (WENO) scheme solving one-dimensional (1D) and two-dimensional (2D) 3-T RH equations respectively. Following our previous paper [7] , we introduce the three new energy variables, and then design a finite difference scheme with both the conservative property and arbitrary high order accuracy. Based on the WENO interpolation and the strong stability preserving (SSP) high order time discretizations, taken as an example, we design a class of fifth order conservative finite difference schemes in space and third order in time. Compared with the Lagrangian method we proposed in [7] , which can only reach second order accuracy for 2D 3-T RH equations if straight-line edged meshes are used, the finite difference scheme can be easily designed to arbitrary high order accuracy for multi-dimensional 3-T RH equations. The finite difference formulation is also much less expensive in multi-dimensions than finite volume schemes used in [7]. Furthermore, our method can handle fluids with large deformation easily. Numerous 1D and 2D numerical examples are presented to verify the desired properties of the high order finite difference WENO schemes such as high order accuracy, non-oscillation, conservation and adaptation to severely distorted single-material radiation hydrodynamics problems. [ABSTRACT FROM AUTHOR]

Since the classical WENO schemes [27] might suffer from slight post-shock oscillations (which are responsible for the numerical residual to hang at a truncation error level) and the new high-order multi-resolution WENO schemes [59] are successful to solve for steady-state problems, we apply these high-order finite volume multi-resolution WENO techniques to serve as limiters for high-order Runge-Kutta discontinuous Galerkin (RKDG) methods in simulating steady-state problems. Firstly, a new troubled cell indicator is designed to precisely detect the cells which would need further limiting procedures. Then the high-order multi-resolution WENO limiting procedures are adopted on a sequence of hierarchical L 2 projection polynomials of the DG solution within the troubled cell itself. By doing so, these RKDG methods with multi-resolution WENO limiters could gradually degrade from the optimal high-order accuracy to the first-order accuracy near strong discontinuities, suppress the slight post-shock oscillations, and push the numerical residual to settle down to machine zero in steady-state simulations. These new multi-resolution WENO limiters are very simple to construct and can be easily implemented to arbitrary high-order accuracy for solving steady-state problems in multi-dimensions. [ABSTRACT FROM AUTHOR]

In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L2-stability of the numerical scheme for fully nonlinear equations. Optimal error estimates (O(h(k+1))) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method. [ABSTRACT FROM AUTHOR]

In this article, we study the spectral volume (SV) methods for scalar hyperbolic conservation laws with a class of subdivision points under the Petrov–Galerkin framework. Due to the strong connection between the DG method and the SV method with the appropriate choice of the subdivision points, it is natural to analyze the SV method in the Galerkin form and derive the analogous theoretical results as in the DG method. This article considers a class of SV methods, whose subdivision points are the zeros of a specific polynomial with a parameter in it. Properties of the piecewise constant functions under this subdivision, including the orthogonality between the trial solution space and test function space, are provided. With the aid of these properties, we are able to derive the energy stability, optimal a priori error estimates of SV methods with arbitrary high order accuracy. We also study the superconvergence of the numerical solution with the correction function technique, and show the order of superconvergence would be different with different choices of the subdivision points. In the numerical experiments, by choosing different parameters in the SV method, the theoretical findings are confirmed by the numerical results. [ABSTRACT FROM AUTHOR]

The weighted essentially non-oscillatory technique using a stencil of 2 r points (WENO-2 r) is an interpolatory method that consists in obtaining a higher approximation order from the non-linear combination of interpolants of r + 1 nodes. The result is an interpolant of order 2 r at the smooth parts and order r + 1 when an isolated discontinuity falls at any grid interval of the large stencil except at the central one. Recently, a new WENO method based on Aitken-Neville's algorithm has been designed for interpolation of equally spaced data at the mid-points and presents progressive order of accuracy close to discontinuities. This paper is devoted to constructing a general progressive WENO method for non-necessarily uniformly spaced data and several variables interpolating in any point of the central interval. Also, we provide explicit formulas for linear and non-linear weights and prove the order obtained. Finally, some numerical experiments are presented to check the theoretical results. [ABSTRACT FROM AUTHOR]

In this paper, high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (WENO) limiters are designed for solving hyperbolic conservation laws on triangular meshes. These multi-resolution WENO limiters are new extensions of the associated multi-resolution WENO finite volume schemes [49,50] which serve as limiters for RKDG methods from structured meshes [47] to triangular meshes. Such new WENO limiters use information of the DG solution essentially only within the troubled cell itself which is identified by a new modified version of the original KXRCF indicator [24] , to build a sequence of hierarchical L 2 projection polynomials from zeroth degree to the highest degree of the RKDG method. The second-order, third-order, and fourth-order RKDG methods with associated multi-resolution WENO limiters are developed as examples, which could maintain the original order of accuracy in smooth regions and keep essentially non-oscillatory property near strong shocks or contact discontinuities by gradually degrading from the highest order to the first order. The linear weights inside the procedure of the new multi-resolution WENO limiters can be any positive numbers on the condition that their sum equals one. This is the first time that a series of polynomials of different degrees within the troubled cell itself are applied in a WENO fashion to modify the DG solutions in the troubled cell on triangular meshes. These new WENO limiters are very simple to construct, and can be easily implemented to arbitrary high-order accuracy and in higher dimensions on unstructured meshes. Such spatial reconstruction methodology improves the robustness in the simulation on the same compact spatial stencil of the original DG methods on triangular meshes. Extensive one-dimensional (run as two-dimensional problems on triangular meshes) and two-dimensional tests are performed to demonstrate the effectiveness of these RKDG methods with the new multi-resolution WENO limiters. [ABSTRACT FROM AUTHOR]

In this paper, we study the optimal error estimates of the classical discontinuous Galerkin method for time-dependent 2-D hyperbolic equations using Pk elements on uniform Cartesian meshes, and prove that the error in the L2 norm achieves optimal (k + 1)th order convergence when upwind fluxes are used. For the linear constant coefficient case, the results hold true for arbitrary piecewise polynomials of degree k ≥ 0. For variable coefficient and nonlinear cases, we give the proof for piecewise polynomials of degree k = 0, 1, 2, 3 and k = 2, 3, respectively, under the condition that the wind direction does not change. The theoretical results are verified by numerical examples. [ABSTRACT FROM AUTHOR]

Two new classes of compact weighted essentially nonoscillatory (WENO) polynomial limiters are presented for second-, third-, fourth-, and fifth-order discontinuous Galerkin (DG) schemes on irregular simplex elements. The presented WENO-DG procedures are extensions of the high-order WENO finite-volume and finite-difference schemes of Zhu and Shu (2017) [25] , (2019) [26] to high-order unstructured DG schemes. A compact positivity preserving limiter is applied to the solutions to ensure pressure and density remain within physical ranges at all time. It is then verified that the bounded WENO-DG maintains the formal order of accuracy of the underlying DG schemes in the smooth regions. The performance of the proposed WENO-DG is also demonstrated with inviscid test cases including the classical Riemann problems, shock-turbulence interaction, scramjet, blunt body flows, and the double Mach Reflection problems. • Compact high-order WENO for DG methods on triangular elements. • Bounded polynomial limiter with a positivity preserving. • (k + 1) th order of accuracy on irregular triangular elements. • Nonintrusive and simple implementation procedures. [ABSTRACT FROM AUTHOR]

• High order energy conserving DG methods which converge in optimal order are designed for linear, symmetric hyperbolic systems. • These DG methods are compared with upwind DG methods and energy conserving DG methods with central fluxes to demonstrate their superior performance especially for long time simulation. • Arbitrary high-order Lax-Wendroff type, energy-conserving time stepping methods are developed. We propose energy-conserving discontinuous Galerkin (DG) methods for symmetric linear hyperbolic systems on general unstructured meshes. Optimal a priori error estimates of order k + 1 are obtained for the semi-discrete scheme in one dimension, and in multi-dimensions on Cartesian meshes when tensor-product polynomials of degree k are used. A high-order energy-conserving Lax-Wendroff time discretization is also presented. Extensive numerical results in one dimension, and two dimensions on both rectangular and triangular meshes are presented to support the theoretical findings and to assess the new methods. One particular method (with the doubling of unknowns) is found to be optimally convergent on triangular meshes for all the examples considered in this paper. The method is also compared with the classical (dissipative) upwinding DG method and (conservative) DG method with a central flux. It is numerically observed for the new method to have a superior performance for long time simulations. [ABSTRACT FROM AUTHOR]

Recently, there has been a series of works on the positivity-preserving discontinuous Galerkin methods for stationary hyperbolic equations, where the notion of mass conservation follows from a straightforward analogy of that of time-dependent problems, i.e. conserving the mass = preserving cell averages during limiting. Based on such a notion, the implementations and theoretical proofs of positivity-preserving limited methods for stationary equations are unnecessarily complicated and constrained. As will be shown in this paper, in some extreme cases, their convergence could even be problematic. In this work, we clarify a more appropriate definition of mass conservation for limiters applied to stationary hyperbolic equations and establish the genuinely conservative high-order positivity-preserving limited discontinuous Galerkin methods based on this definition. The new methods are able to preserve the positivity of solutions of scalar linear equations and scalar nonlinear equations with invariant wind direction, with much simpler implementations and easier proofs for accuracy and the Lax-Wendroff theorem, compared with the existing methods. Two types of positivity-preserving limiters preserving the local mass of stationary equations are developed to accommodate for the new definition of conservation and their accuracy are investigated. We would like to emphasize that a major advantage of the original DG scheme presented in [24] is a sweeping procedure, which allows for the computation of conservative steady-state solutions explicitly, cell by cell, without iterations, even for nonlinear equations as long as the wind direction is fixed. The main contribution of this paper is to introduce a limiting procedure to enforce positivity without changing the conservative property of this original DG scheme. The good performance of the algorithms for stationary hyperbolic equations and their applications in time-dependent problems are demonstrated by ample numerical tests. • A new definition of local conservation is given for stationary hyperbolic systems. • This allows the design of positivity-preserving discontinuous Galerkin (DG) schemes in more general cases than before. • Such high order positivity-preserving DG schemes are more general than before. [ABSTRACT FROM AUTHOR]

In this paper, we continue our work in [46] and propose a new type of high-order finite volume multi-resolution weighted essentially non-oscillatory (WENO) schemes to solve hyperbolic conservation laws on triangular meshes. Although termed "multi-resolution WENO schemes", we only use the information defined on a hierarchy of nested central spatial stencils and do not introduce any equivalent multi-resolution representation. We construct new third-order, fourth-order, and fifth-order WENO schemes using three or four unequal-sized central spatial stencils, different from the classical WENO procedure using equal-sized biased/central spatial stencils for the spatial reconstruction. The new WENO schemes could obtain the optimal order of accuracy in smooth regions, and could degrade gradually to first-order of accuracy so as to suppress spurious oscillations near strong discontinuities. This is the first time that only a series of unequal-sized hierarchical central spatial stencils are used in designing arbitrary high-order finite volume WENO schemes on triangular meshes. The main advantages of these schemes are their compactness, robustness, and their ability to maintain good convergence property for steady-state computation. The linear weights of such WENO schemes can be any positive numbers on the condition that they sum to one. Extensive numerical results are provided to illustrate the good performance of these new finite volume WENO schemes. • A new type of multi-resolution WENO schemes on triangular meshes has been developed. • High order WENO schemes using unequal-sized central spatial stencils have been designed. • The main advantages are compactness, robustness, and good steady-state convergence. [ABSTRACT FROM AUTHOR]

Abstract We propose an energy-conserving ultra-weak discontinuous Galerkin (DG) method for the generalized Korteweg– de Vries (KdV) equation in one dimension. Optimal a priori error estimate of order k + 1 is obtained for the semi-discrete scheme for the KdV equation without the convection term on general nonuniform meshes when polynomials of degree k ≥ 2 are used. We also numerically observe optimal convergence of the method for the KdV equation with linear or nonlinear convection terms. It is numerically observed for the new method to have a superior performance for long-time simulations over existing DG methods. [ABSTRACT FROM AUTHOR]

Abstract In this paper, a new type of high-order finite difference and finite volume multi-resolution weighted essentially non-oscillatory (WENO) schemes is presented for solving hyperbolic conservation laws. We only use the information defined on a hierarchy of nested central spatial stencils and do not introduce any equivalent multi-resolution representation. These new WENO schemes use the same large stencils as the classical WENO schemes in [25,45] , could obtain the optimal order of accuracy in smooth regions, and could simultaneously suppress spurious oscillations near discontinuities. The linear weights of such WENO schemes can be any positive numbers on the condition that their sum equals one. This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finite difference and finite volume WENO schemes. These new WENO schemes are simple to construct and can be easily implemented to arbitrary high order of accuracy and in higher dimensions. Benchmark examples are given to demonstrate the robustness and good performance of these new WENO schemes. Highlights • A new class of high order finite difference and finite volume WENO schemes are constructed. • These schemes are based on the multi-resolution idea, and a series of unequal-sized hierarchical central spatial stencils. • These schemes can use arbitrary positive linear weights, and are easy to implement for one and multi-dimensions. • These schemes have a gradual degrading of accuracy near discontinuities. [ABSTRACT FROM AUTHOR]

This paper deals with the numerical resolution of kinetic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. We focus on the kinetic model considered in [P. Degond and S. Motsch, Continuum limit of self-driven particles with orientation interaction, Math. Models Methods Appl. Sci. 18 (2008) 1193-1215; P. Degond, J.-G. Liu, S. Motsch and V. Panferov, Hydrodynamic models of self-organized dynamics: Derivation and existence theory, Methods Appl. Anal. 20 (2013) 89-114] where alignment is taken into account in addition to an attraction-repulsion interaction potential. We apply a discontinuous Galerkin method for the free transport and non-local drift velocity together with a spectral method for the velocity variable. Then, we analyze consistency and stability of the semi-discrete scheme. We propose several numerical experiments which provide a solid validation of the method and illustrate its underlying concepts. [ABSTRACT FROM AUTHOR]

In this paper, we develop bound-preserving modified exponential Runge–Kutta (RK) discontinuous Galerkin (DG) schemes to solve scalar hyperbolic equations with stiff source terms by extending the idea in Zhang and Shu [43] . Exponential strong stability preserving (SSP) high order time discretizations are constructed and then modified to overcome the stiffness and preserve the bound of the numerical solutions. It is also straightforward to extend the method to two dimensions on rectangular and triangular meshes. Even though we only discuss the bound-preserving limiter for DG schemes, it can also be applied to high order finite volume schemes, such as weighted essentially non-oscillatory (WENO) finite volume schemes as well. [ABSTRACT FROM AUTHOR]

We present a discontinuous Galerkin (DG) scheme with suitable quadrature rules [15] for ideal compressible magnetohydrodynamic (MHD) equations on structural meshes. The semi-discrete scheme is analyzed to be entropy stable by using the symmetrizable version of the equations as introduced by Godunov [32] , the entropy stable DG framework with suitable quadrature rules [15] , the entropy conservative flux in [14] inside each cell and the entropy dissipative approximate Godunov type numerical flux at cell interfaces to make the scheme entropy stable. The main difficulty in the generalization of the results in [15] is the appearance of the non-conservative “source terms” added in the modified MHD model introduced by Godunov [32] , which do not exist in the general hyperbolic system studied in [15] . Special care must be taken to discretize these “source terms” adequately so that the resulting DG scheme satisfies entropy stability. Total variation diminishing / bounded (TVD/TVB) limiters and bound-preserving limiters are applied to control spurious oscillations. We demonstrate the accuracy and robustness of this new scheme on standard MHD examples. [ABSTRACT FROM AUTHOR]

A new class of high order weighted essentially non-oscillatory (WENO) schemes (Zhu and Qiu, 2016, [50] ) is applied to solve Euler equations with steady state solutions. It is known that the classical WENO schemes (Jiang and Shu, 1996, [23] ) might suffer from slight post-shock oscillations. Even though such post-shock oscillations are small enough in magnitude and do not visually affect the essentially non-oscillatory property, they are truly responsible for the residue to hang at a truncation error level instead of converging to machine zero. With the application of this new class of WENO schemes, such slight post-shock oscillations are essentially removed and the residue can settle down to machine zero in steady state simulations. This new class of WENO schemes uses a convex combination of a quartic polynomial with two linear polynomials on unequal size spatial stencils in one dimension and is extended to two dimensions in a dimension-by-dimension fashion. By doing so, such WENO schemes use the same information as the classical WENO schemes in Jiang and Shu (1996) [23] and yield the same formal order of accuracy in smooth regions, yet they could converge to steady state solutions with very tiny residue close to machine zero for our extensive list of test problems including shocks, contact discontinuities, rarefaction waves or their interactions, and with these complex waves passing through the boundaries of the computational domain. [ABSTRACT FROM AUTHOR]

We introduce a new troubled-cell indicator for the discontinuous Galerkin (DG) methods for solving hyperbolic conservation laws. This indicator can be defined on unstructured meshes for high order DG methods and depends only on data from the target cell and its immediate neighbors. It is able to identify shocks without PDE sensitive parameters to tune. Extensive one- and two-dimensional simulations on the hyperbolic systems of Euler equations indicate the good performance of this new troubled-cell indicator coupled with a simple minmod-type TVD limiter for the Runge–Kutta DG (RKDG) methods. [ABSTRACT FROM AUTHOR]

In this paper, we consider the linear convection-diffusion equation in one dimension with periodic boundary conditions, and analyze the stability of fully discrete methods that are defined with local discontinuous Galerkin (LDG) methods in space and several implicit-explicit (IMEX) Runge-Kutta methods in time. By using the forward temporal differences and backward temporal differences, respectively, we establish two general frameworks of the energy-method based stability analysis. From here, the fully discrete schemes being considered are shown to have monotonicity stability, i.e. the L^2 norm of the numerical solution does not increase in time, under the time step condition \tau \le \mathcal {F}(h/c, d/c^2), with the convection coefficient c, the diffusion coefficient d, and the mesh size h. The function \mathcal {F} depends on the specific IMEX temporal method, the polynomial degree k of the discrete space, and the mesh regularity parameter. Moreover, the time step condition becomes \tau \lesssim h/c in the convection-dominated regime and it becomes \tau \lesssim d/c^2 in the diffusion-dominated regime. The result is improved for a first order IMEX-LDG method. To complement the theoretical analysis, numerical experiments are further carried out, leading to slightly stricter time step conditions that can be used by practitioners. Uniform stability with respect to the strength of the convection and diffusion effects can especially be relevant to guide the choice of time step sizes in practice, e.g. when the convection-diffusion equations are convection-dominated in some sub-regions. [ABSTRACT FROM AUTHOR]

It is well known that semi-discrete high order discontinuous Galerkin (DG) methods satisfy cell entropy inequalities for the square entropy for both scalar conservation laws (Jiang and Shu (1994) [39] ) and symmetric hyperbolic systems (Hou and Liu (2007) [36] ), in any space dimension and for any triangulations. However, this property holds only for the square entropy and the integrations in the DG methods must be exact. It is significantly more difficult to design DG methods to satisfy entropy inequalities for a non-square convex entropy, and/or when the integration is approximated by a numerical quadrature. In this paper, we develop a unified framework for designing high order DG methods which will satisfy entropy inequalities for any given single convex entropy, through suitable numerical quadrature which is specific to this given entropy. Our framework applies from one-dimensional scalar cases all the way to multi-dimensional systems of conservation laws. For the one-dimensional case, our numerical quadrature is based on the methodology established in Carpenter et al. (2014) [5] and Gassner (2013) [19] . The main ingredients are summation-by-parts (SBP) operators derived from Legendre Gauss–Lobatto quadrature, the entropy conservative flux within elements, and the entropy stable flux at element interfaces. We then generalize the scheme to two-dimensional triangular meshes by constructing SBP operators on triangles based on a special quadrature rule. A local discontinuous Galerkin (LDG) type treatment is also incorporated to achieve the generalization to convection–diffusion equations. Extensive numerical experiments are performed to validate the accuracy and shock capturing efficacy of these entropy stable DG methods. [ABSTRACT FROM AUTHOR]

In this paper, a new type of finite difference Hermite weighted essentially non-oscillatory (HWENO) schemes are constructed for solving Hamilton–Jacobi (HJ) equations. Point values of both the solution and its first derivatives are used in the HWENO reconstruction and evolved via time advancing. While the evolution of the solution is still through the classical numerical fluxes to ensure convergence to weak solutions, the evolution of the first derivatives of the solution is through a simple dimension-by-dimension non-conservative procedure to gain efficiency. The main advantages of this new scheme include its compactness in the spatial field and its simplicity in the reconstructions. Extensive numerical experiments in one and two dimensional cases are performed to verify the accuracy, high resolution and efficiency of this new scheme. [ABSTRACT FROM AUTHOR]

This paper is concerned with an interior penalty discontinuous Galerkin (IPDG) method based on a flexible type of non-polynomial local approximation space for the Helmholtz equation with varying wavenumber. The local approximation space consists of multiple polynomial-modulated phase functions which can be chosen according to the phase information of the solution. We obtain some approximation properties for this space and a priori L 2 error estimates for the h -convergence of the IPDG method using duality argument. We also provide ample numerical examples to show that, building phase information into the local spaces often gives more accurate results comparing to using the standard polynomial spaces. [ABSTRACT FROM AUTHOR]

In this paper, we analyze the Lax–Wendroff discontinuous Galerkin (LWDG) method for solving linear conservation laws. The method was originally proposed by Guo et al. in [W. Guo, J.-M. Qiu and J. Qiu, J. Sci. Comput. 65 (2015) 299–326], where they applied local discontinuous Galerkin (LDG) techniques to approximate high order spatial derivatives in the Lax–Wendroff time discretization. We show that, under the standard CFL condition τ ≤ λh (where τ and h are the time step and the maximum element length respectively and λ > 0 is a constant) and uniform or non-increasing time steps, the second order schemes with piecewise linear elements and the third order schemes with arbitrary piecewise polynomial elements are stable in the L2 norm. The specific type of stability may differ with different choices of numerical fluxes. Our stability analysis includes multidimensional problems with divergence-free coefficients. Besides solving the equation itself, the LWDG method also gives approximations to its time derivative simultaneously. We obtain optimal error estimates for both the solution u and its first order time derivative ut in one dimension, and numerical examples are given to validate our analysis. [ABSTRACT FROM AUTHOR]

In this paper, we develop a second-order asymptotic-preserving and positivity-preserving discontinuous Galerkin (DG) scheme for the Kerr-Debye model. By using the approach first introduced by Zhang and Shu in [Q. Zhang and C.-W. Shu, Error estimates to smooth solutions of Runge-Kutta discontinuous Galerkin methods for scalar conservation laws, SIAM J. Numer. Anal. 42 (2004) 641-666.] with an energy estimate and Taylor expansion, the asymptotic-preserving property of the semi-discrete DG methods is proved rigorously. In addition, we propose a class of unconditional positivity-preserving implicit-explicit (IMEX) Runge-Kutta methods for the system of ordinary differential equations arising from the semi-discretization of the Kerr-Debye model. The new IMEX Runge-Kutta methods are based on the modification of the strong-stability-preserving (SSP) implicit Runge-Kutta method and have second-order accuracy. The numerical results validate our analysis. [ABSTRACT FROM AUTHOR]

In this article, we focus on error estimates to smooth solutions of semi-discrete discontinuous Galerkin (DG) methods with quadrature rules for scalar conservation laws. The main techniques we use are energy estimate and Taylor expansion first introduced by Zhang and Shu in (Zhang and Shu, SIAM J Num Anal 42 (2004), 641-666). We show that, with [ABSTRACT FROM AUTHOR]

Numerical simulation of seismic wave equations has attracted much attention and plays a significant role in exploration seismology. As one of seismic wave models, the diffusive–viscous wave theory usually describes the attenuation of seismic wave propagating in fluid-saturated medium. In this paper, we focus on the design of numerical methods for the diffusive–viscous wave equations with variable coefficients. We develop a local discontinuous Galerkin (LDG) method, in which numerical fluxes are chosen carefully to maintain stability and accuracy. Moreover, we also prove the optimal error estimates for both the energy norm and the L 2 norm. Numerical experiments are provided to demonstrate the optimal convergence rate and effectiveness of the proposed LDG method. • Diffusive-viscous wave equations with variable coefficient is a key model in geophysics. • For this model, we design a high order LDG method, and analyzed the energy stability. • We provide the optimal error estimates both for the energy norm and L 2 norm. • Numerical examples demonstrate accuracy and effectiveness of the proposed method. [ABSTRACT FROM AUTHOR]

There have been intensive studies on maximum-principle-satisfying and positivity-preserving methods for hyperbolic conservation laws. Most of them are based on the method of lines type time marching approaches, e.g. the Runge-Kutta methods, multi-step methods and backward Euler method. As an alternative, the Lax-Wendroff time marching approach utilizes the information of PDEs in the Taylor expansion of the solution in time, hence it is a high order and single-stage method. In this work, we propose third order maximum-principle-satisfying and positivity-preserving schemes for scalar conservation laws and the Euler equations based on the Lax-Wendroff time discretization and discontinuous Galerkin spatial discretization. The accuracy and effectiveness of the maximum-principle-satisfying and positivity-preserving techniques are demonstrated by ample numerical tests. • Third order positivity-preserving Lax-Wendroff DG schemes for conservation laws are designed. • Lax-Wendroff procedure without certain mixed derivatives and suitable usage of DDG techniques ensure cell-average positivity. • The Zhang-Shu framework is then adopted with a scaling limiter which maintains high order accuracy when enforcing positivity. [ABSTRACT FROM AUTHOR]

In this paper, we improve the simple WENO limiter for the Runge-Kutta discontinuous Galerkin (RKDG) method in [34] for solving two-dimensional hyperbolic systems on unstructured meshes. The major improvement is reducing the number of polynomials transformed to the characteristic fields for each direction. For the triangular mesh, this new simple WENO limiter transforms only two polynomials to the characteristic fields for each direction, while the original simple WENO limiter [34] uses four polynomials. Thus the improved simple WENO limiter reduces the computational cost and improves the efficiency. It provides a simpler and more practical and efficient way to the characteristic-wise limiting procedure, while still simultaneously maintaining uniform high-order accuracy in smooth regions and controlling spurious nonphysical oscillations near discontinuities. We also apply this improved simple WENO limiter to a high order method constructed for solving hyperbolic conservation laws on arbitrarily distributed point clouds [10] , where polygonal meshes are constructed based on random points and the traditional DG method was adopted on the constructed polygonal mesh. For such a complex polygonal mesh, this limiter still transforms only two polynomials to the characteristic fields. Thus the simplicity and efficiency of this new limiter is more evident in this case. Numerical results are provided to illustrate the accuracy and effectiveness of this procedure. For some examples, the improved limiter has less smearing and higher resolution than the original one. • Improve the simple WENO limiter for hyperbolic systems on structured meshes. • Reduces the number of polynomials transformed in each direction. • Limiter for DG methods on complex polygonal mesh. • Compact stencil and simplicity of linear weights. • High-order accuracy and non-oscillatory property. [ABSTRACT FROM AUTHOR]

This is a follow-up work of Yuan et al. (2016) [19] and Ling et al. (2018) [13] that further investigates the positivity-preserving discontinuous Galerkin (DG) methods for stationary hyperbolic equations. In 2016, Yuan et al. proposed a high order positivity-preserving DG method for stationary hyperbolic equations with constant coefficients, but the scheme has to be used in combination with a non-conservative rotational limiter introduced in case of negative cell averages. Ling et al. (2018) improved the results in one dimensional space by rigorously proving the positivity of cell averages of the unmodulated DG scheme, which allows the conservative scaling limiter in Zhang and Shu (2010) [22] to be used to maintain positivity without affecting accuracy, but extension to two space dimensions requires an augmentation of the DG space and works only in the second order case. Considering that the aforementioned works only address stationary hyperbolic equations with constant coefficients and higher than second order conservative methods are still unavailable in two and three space dimensions, we propose high order conservative positivity-preserving DG methods for variable coefficient and nonlinear stationary hyperbolic equations in one dimension and constant coefficient stationary hyperbolic equations in two and three dimensions, via a suitable quadrature in the DG framework. We show the good performance of the algorithms by ample numerical experiments, including their applications in time-dependent problems. • High order conservative positivity-preserving DG for stationary hyperbolic equations. • In 1D, schemes are designed for variable coefficient and nonlinear equations. • In 2D and 3D, high order schemes are designed for constant coefficient equations. [ABSTRACT FROM AUTHOR]

For solving time-dependent convection-dominated partial differential equations (PDEs), which arise frequently in computational physics, high order numerical methods, including finite difference, finite volume, finite element and spectral methods, have been undergoing rapid developments over the past decades. In this article we give a brief survey of two selected classes of high order methods, namely the weighted essentially non-oscillatory (WENO) finite difference and finite volume schemes and discontinuous Galerkin (DG) finite element methods, emphasizing several of their recent developments: bound-preserving limiters for DG, finite volume and finite difference schemes, which address issues in robustness and accuracy; WENO limiters for DG methods, which address issues in non-oscillatory performance when there are strong shocks, and inverse Lax–Wendroff type boundary treatments for finite difference schemes, which address issues in solving complex geometry problems using Cartesian meshes. [ABSTRACT FROM AUTHOR]

In this paper, we develop a discontinuous Galerkin (DG) method to solve the ideal special relativistic hydrodynamics (RHD) and design a bound-preserving (BP) limiter for this scheme by extending the idea in X. Zhang and C.-W. Shu, (2010) [56] . For RHD, the density and pressure are positive and the velocity is bounded by the speed of light. One difficulty in numerically solving the RHD in its conservative form is that the failure of preserving these physical bounds will result in ill-posedness of the problem and blowup of the code, especially in extreme relativistic cases. The standard way in dealing with this difficulty is to add extra numerical dissipation, while in doing so there is no guarantee of maintaining the high order of accuracy. Our BP limiter has the following features. It can theoretically guarantee to preserve the physical bounds for the numerical solution and maintain its designed high order accuracy. The limiter is local to the cell and hence is very easy to implement. Moreover, it renders L 1 -stability to the numerical scheme. Numerical experiments are performed to demonstrate the good performance of this bound-preserving DG scheme. Even though we only discuss the BP limiter for DG schemes, it can be applied to high order finite volume schemes, such as weighted essentially non-oscillatory (WENO) finite volume schemes as well. [ABSTRACT FROM AUTHOR]

In this paper, we consider linear stability issues for one-dimensional hyperbolic conservation laws using a class of conservative high order upwind-biased finite difference schemes, which is a prototype for the weighted essentially non-oscillatory (WENO) schemes, for initial–boundary value problems (IBVP). The inflow boundary is treated by the so-called inverse Lax–Wendroff (ILW) or simplified inverse Lax–Wendroff (SILW) procedure, and the outflow boundary is treated by the classical high order extrapolation. A third order total variation diminishing (TVD) Runge–Kutta time discretization is used in the fully discrete case. Both GKS (Gustafsson, Kreiss and Sundström) and eigenvalue analyses are performed for both semi-discrete and fully discrete schemes. The two different analysis techniques yield consistent results. Numerical tests are performed to demonstrate the stability results predicted by the analysis. [ABSTRACT FROM AUTHOR]

This paper is the second part of a series of two. It follows [44] , in which the positivity-preservation property of methods solving one-dimensional Lagrangian gas dynamics equations, from first-order to high-orders of accuracy, was addressed. This article aims at extending this analysis to the two-dimensional case. This study is performed on a general first-order cell-centered finite volume formulation based on polygonal meshes defined either by straight line edges, conical edges, or any high-order curvilinear edges. Such formulation covers the numerical methods introduced in [6,32,5,41,43] . This positivity study is then extended to high-orders of accuracy. Through this new procedure, scheme robustness is highly improved and hence new problems can be tackled. Numerical results are provided to demonstrate the effectiveness of these methods. It is important to point out that even if this paper is concerned with purely Lagrangian schemes, the theory developed is of fundamental importance for any methods relying on a purely Lagrangian step, as ALE methods or non-direct Euler schemes. [ABSTRACT FROM AUTHOR]

One of the main issues in the field of numerical schemes is to ally robustness with accuracy. Considering gas dynamics, numerical approximations may generate negative density or pressure, which may lead to nonlinear instability and crash of the code. This phenomenon is even more critical using a Lagrangian formalism, the grid moving and being deformed during the calculation. Furthermore, most of the problems studied in this framework contain very intense rarefaction and shock waves. In this paper, the admissibility of numerical solutions obtained by high-order finite-volume-scheme-based methods, such as the discontinuous Galerkin (DG) method, the essentially non-oscillatory (ENO) and the weighted ENO (WENO) finite volume schemes, is addressed in the one-dimensional Lagrangian gas dynamics framework. After briefly recalling how to derive Lagrangian forms of the 1D gas dynamics system of equations, a discussion on positivity-preserving approximate Riemann solvers, ensuring first-order finite volume schemes to be positive, is then given. This study is conducted for both ideal gas and non-ideal gas equations of state (EOS), such as the Jones–Wilkins–Lee (JWL) EOS or the Mie–Grüneisen (MG) EOS, and relies on two different techniques: either a particular definition of the local approximation of the acoustic impedances arising from the approximate Riemann solver, or an additional time step constraint relative to the cell volume variation. Then, making use of the work presented in [89,90,22] , this positivity study is extended to high-orders of accuracy, where new time step constraints are obtained, and proper limitation is required. Through this new procedure, scheme robustness is highly improved and hence new problems can be tackled. Numerical results are provided to demonstrate the effectiveness of these methods. This paper is the first part of a series of two. The whole analysis presented here is extended to the two-dimensional case in [85] , and proves to fit a wide range of numerical schemes in the literature, such as those presented in [19,64,15,82,84] . [ABSTRACT FROM AUTHOR]

The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with implicit–explicit (IMEX) time discretization schemes, for solving one-dimensional convection–diffusion equations with a nonlinear convection. Both Runge–Kutta and multi-step IMEX methods are considered. By the aid of the energy method, we show that the IMEX LDG schemes are unconditionally stable for the nonlinear problems, in the sense that the time-step τ is only required to be upper-bounded by a positive constant which depends on the flow velocity and the diffusion coefficient, but is independent of the mesh size h . We also give optimal error estimates for the IMEX LDG schemes, under the same temporal condition, if a monotone numerical flux is adopted for the convection. Numerical experiments are given to verify our main results. [ABSTRACT FROM AUTHOR]

One- and multi-dimensional stochastic Maxwell equations with additive noise are considered in this paper. It is known that such system can be written in the multi-symplectic structure, and the stochastic energy increases linearly in time. High order discontinuous Galerkin methods are designed for the stochastic Maxwell equations with additive noise, and we show that the proposed methods satisfy the discrete form of the stochastic energy linear growth property and preserve the multi-symplectic structure on the discrete level. Optimal error estimate of the semi-discrete DG method is also analyzed. The fully discrete methods are obtained by coupling with symplectic temporal discretizations. One- and two-dimensional numerical results are provided to demonstrate the performance of the proposed methods, and optimal error estimates and linear growth of the discrete energy can be observed for all cases. [ABSTRACT FROM AUTHOR]

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes [45] , to the correction procedure via reconstruction (CPR) framework for solving conservation laws. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper is particularly simple to implement and will not harm the conservativeness of the CPR framework. Also, it uses information only from the target cell and its immediate neighbors, and thus can maintain the compactness of the CPR framework. Since the CPR framework with the WENO limiter does not in general preserve positivity of the solution, we also extend the positivity-preserving limiters in [43,44,42] to the CPR framework. Numerical results in one and two dimensions are provided to illustrate the good behavior of this procedure. [ABSTRACT FROM AUTHOR]

In this paper, we discuss high-order finite difference weighted essentially non-oscillatory schemes, coupled with total variation diminishing (TVD) Runge-Kutta (RK) temporal integration, for solving the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology. Since the solutions to this system are non-negative, we discuss a positivity-preserving limiter without compromising accuracy. Analysis is performed to justify the maintenance of third-order spatial/temporal accuracy when the limiters are applied to a third-order finite difference scheme and third-order TVD-RK time discretization for solving this model. Numerical results are also provided to demonstrate these methods up to fifth-order accuracy. [ABSTRACT FROM AUTHOR]

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes on two-dimensional unstructured triangular meshes [39] , to the correction procedure via reconstruction (CPR) framework for solving nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes with straight or curved edges. This is an extension of our earlier work [4] in which the WENO limiter was designed for the CPR framework on regular meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper uses information only from the target cell and its immediate neighbors. Hence, it is particularly simple to implement and will not harm the conservativeness and compactness of the CPR framework. Since the CPR framework with this WENO limiter does not in general satisfy the positivity preserving property, we also extend the positivity-preserving limiters [36,33] to the CPR framework. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of this procedure. [ABSTRACT FROM AUTHOR]

Abstract: In applications such as astrophysics and inertial confinement fusion, there are many three-dimensional cylindrical-symmetric multi-material problems which are usually simulated by Lagrangian schemes in the two-dimensional cylindrical coordinates. For this type of simulation, a critical issue for the schemes is to keep spherical symmetry in the cylindrical coordinate system if the original physical problem has this symmetry. In the past decades, several Lagrangian schemes with such symmetry property have been developed, but all of them are only first order accurate. In this paper, we develop a second order cell-centered Lagrangian scheme for solving compressible Euler equations in cylindrical coordinates, based on the control volume discretizations, which is designed to have uniformly second order accuracy and capability to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. The scheme maintains several good properties such as conservation for mass, momentum and total energy, and the geometric conservation law. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of accuracy, symmetry, non-oscillation and robustness. The advantage of higher order accuracy is demonstrated in these examples. [Copyright &y& Elsevier]

Abstract: Solving wave propagation problems within heterogeneous media has been of great interest and has a wide range of applications in physics and engineering. The design of numerical methods for such general wave propagation problems is challenging because the energy conserving property has to be incorporated in the numerical algorithms in order to minimize the phase or shape errors after long time integration. In this paper, we focus on multi-dimensional wave problems and consider linear second-order wave equation in heterogeneous media. We develop and analyze an LDG method, in which numerical fluxes are carefully designed to maintain the energy conserving property and accuracy. Compatible high order energy conserving time integrators are also proposed. The optimal error estimates and the energy conserving property are proved for the semi-discrete methods. Our numerical experiments demonstrate optimal rates of convergence, and show that the errors of the numerical solutions do not grow significantly in time due to the energy conserving property. [Copyright &y& Elsevier]

Abstract: Gibbs phenomenon is the particular manner how a global spectral approximation of a piecewise analytic function behaves at the jump discontinuity. The truncated spectral series has large oscillations near the jump, and the overshoot does not decay as the number of terms in the truncated series increases. There is therefore no convergence in the maximum norm, and convergence in smooth regions away from the discontinuity is also slow. In Gottlieb and Shu (1995) [5], a methodology is proposed to completely overcome this difficulty in the context of spectral collocation methods, resulting in the recovery of exponential accuracy from collocation point values of a piecewise analytic function. In this paper, we extend this methodology to handle spectral collocation methods for functions which are analytic in the open interval but have singularities at end-points. With this extension, we are able to obtain exponential accuracy from collocation point values of such functions. Similar to Gottlieb and Shu (1995) [5], the proof is constructive and uses the Gegenbauer polynomials . The result implies that the Gibbs phenomenon can be overcome for smooth functions with endpoint singularities. [Copyright &y& Elsevier]

Abstract: Robustness of numerical methods has attracted an increasing interest in the community of computational fluid dynamics. One mathematical aspect of robustness for numerical methods is the positivity-preserving property. At high Mach numbers or for flows near vacuum, solving the conservative Euler equations may generate negative density or internal energy numerically, which may lead to nonlinear instability and crash of the code. This difficulty is particularly profound for high order methods, for multi-material flows and for problems with moving meshes, such as the Lagrangian methods. In this paper, we construct both first order and uniformly high order accurate conservative Lagrangian schemes which preserve positivity of physically positive variables such as density and internal energy in the simulation of compressible multi-material flows with general equations of state (EOS). We first develop a positivity-preserving approximate Riemann solver for the Lagrangian scheme solving compressible Euler equations with both ideal and non-ideal EOS. Then we design a class of high order positivity-preserving and conservative Lagrangian schemes by using the essentially non-oscillatory (ENO) reconstruction, the strong stability preserving (SSP) high order time discretizations and the positivity-preserving scaling limiter which can be proven to maintain conservation and uniformly high order accuracy and is easy to implement. One-dimensional and two-dimensional numerical tests for the positivity-preserving Lagrangian schemes are provided to demonstrate the effectiveness of these methods. [Copyright &y& Elsevier]

Abstract: In this paper, we develop a three-dimensional parallel solver using the fifth order high-resolution weighted essentially non-oscillatory (WENO) finite difference scheme to perform extensive simulation for three-dimensional gaseous detonations. A careful study is conducted for the propagation modes of three-dimensional gaseous detonation wave-front structures in a long square duct. The numerical results indicate that, the instability of detonation, overdrive factor, transverse dimension and initial perturbation are the main factors that influence the appearance of spinning detonation. For unstable detonation with an overdrive factor f =1.0, with the duct width decreasing, the transverse modes tend to be the lower ones. When the duct width is sufficiently small, detonation propagation is coupled to the lowest transverse mode in order to sustain itself. There is a maximum duct width for the appearance of spinning detonation. When the duct width is smaller than this value, the unstable detonation propagates in a spinning mode. When the duct width is larger than this value, spinning detonation disappears. For highly unstable detonation with an overdrive factor f =1.2, the formation of the spinning detonation does not depend on the initial perturbation and it only depends on the high instability of detonation and channel dimension in the transverse directions. Both initial rectangular modes and initial diagonal modes can eventually evolve into spinning modes. For mildly unstable detonation with an overdrive factor f =1.6, it is easier for the initial diagonal mode to trigger the spinning mode than for the initial rectangular mode. The formation of spinning detonation thus not only depends on the high instability of detonation and channel dimension in the transverse directions, but also depends on the initial perturbation. For the stable detonation with an overdrive factor f =1.1, when the initial perturbation is in the transverse cosine mode, the front of the stable detonation has a rectangular structure and slapping waves appear on the walls. When the initial perturbation is in the symmetrical mode along the diagonals of the detonation front, the diagonal detonation is formed and the slapping waves on the walls disappear. [Copyright &y& Elsevier]