No abstract uploaded!

In this addendum note we fill in the gap left in the description of 2D homogeneous solutions to the stationary Euler system initiated in the previous publication. This completes the classification of all homogeneous stationary solutions. The note includes updated classification tables.

In this paper, we propose new spectral viscosity methods based on the generalized Hermite functions for the solution of nonlinear scalar conservation laws in the whole line. It is shown rigorously that these schemes converge to the unique entropy solution by using compensated compactness arguments, under some conditions. The numerical experiments of the inviscid Burger's equation support our result, and it verifies the reasonableness of the conditions.

It is well-known that sparse grid algorithm has been widely accepted as an efficient tool to overcome the “curse of dimensionality" in some degree. In this note, we give the error estimate of hyperbolic cross (HC) approximations with all sorts of Askey polynomials. These polynomials are useful in generalized polynomial chaos (gPC) in the field of uncertainty quantification. The exponential convergences in both regular and optimized HC approximations have been shown under the condition that the random variable depends on the random inputs smoothly in some degree. Moreover, we apply gPC to numerically solve the ordinary differential equations with slightly higher dimensional random inputs. Both regular and optimized HC have been investigated with Laguerre-chaos, Charlier-chaos and Hermite-chaos in the numerical experiment. The discussion of the connection between the standard ANOVA approximation and Galerkin approximation is in the appendix.

In this follow-up work, we extend the discontinuous Galerkin (DG) methods previously developed on rectangular meshes \cite{previous} to triangular meshes. The DG schemes in \cite{previous} are both optimally convergent and energy conserving. However, as we shall see in the numerical results section, the DG schemes on triangular meshes only have suboptimal convergence rate. We prove the energy conservation and an error estimate for the semi-discrete schemes. The stability of the fully discrete scheme is proved and its error estimate is stated. We present extensive numerical results with convergence consistent of our error estimate, and simulations of wave propagation in Drude metamaterials to demonstrate the flexibility of triangular meshes.

To avoid the order reduction when third order implicit-explicit Runge-Kutta time discretization is used together with the local discontinuous Galerkin (LDG) spatial discretization, for solving convection-diffusion problems with time-dependent Dirichlet boundary conditions, we propose a strategy of boundary treatment at each intermediate stage in this paper. The proposed strategy can achieve optimal order of accuracy by numerical verification. Also by suitably setting numerical flux on the boundary in the LDG methods, and by establishing an important relationship between the gradient and interface jump of the numerical solution with the independent numerical solution of the gradient and the given boundary conditions, we build up the unconditional stability of the corresponding scheme, in the sense that the time step is only required to be upper bounded by a suitable positive constant, which is independent of the mesh size.

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in \cite{yuan}. The DG methods in \cite{yuan} can maintain positivity and high order accuracy, but they rely both on the scaling limiter in \cite{zhang1} and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax-Wendroff type theorem is proved in \cite{yuan}, guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $P^k$ of piecewise polynomials of degree at most $k$. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in \cite{zhang1} can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $P^k$ or $Q^k$ spaces (piecewise $k$-th degree polynomials or piecewise tensor-product $k$-th degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in \cite{zhang1}. Computational results are provided to demonstrate the good performance of our DG schemes.

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 \cite{ZhSh:10max}. 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.

No abstract uploaded!

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 \cite{DM, DLMP} where alignment is taken into account in addition of 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 analyse consistency and stability of the semi-discrete scheme. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

In this paper, we study the superconvergence behavior of the semi-discrete discontinuous Galerkin (DG) method for {\color{red} scalar} nonlinear hyperbolic equations in one spatial dimension. Superconvergence results for problems with fixed and alternating wind directions are established. On the one hand, we prove that, if the wind direction is fixed (i.e., the derivative of the flux function is bounded away from zero), both the cell average error and numerical flux error at cell interfaces converge at a rate of $2k+1$ when upwind fluxes and piecewise polynomials of degree $k$ are used. Moreover, we also prove that the function value approximation of the DG solution is superconvergent at interior right Radau points, and the derivative value approximation is superconvergent at interior left Radau points, with an order of $k+2$ and $k+1$, respectively. As a byproduct, we show a $k+2$-th order superconvergence of the DG solution towards the Gauss-Radau projection of the exact solution. On the other hand, superconvergence results for problems with alternating wind directions (i.e., the derivative of the flux function either changes sign or otherwise achieves the value zero in the domain) are also established. To be more precise, we first prove that the DG flux function is superconvergent towards a particular flux function of the exact solution, with an order of $k+2$, when Godunov fluxes are used. We then prove that the highest superconvergence rate of the DG solution itself is $k+\frac32$ when sonic points (i.e., the derivative of the flux function achieves zero) appear in the computational domain. As byproducts, we obtain superconvergence properties for the DG solution and the DG flux function at special points and for cell average. Numerical experiments demonstrate that most of our results are optimal, i.e., the superconvergence rates are sharp.

We analyze the central discontinuous Galerkin (DG) method for time-dependent linear conservation laws. In one dimension, optimal a priori $L^2$ error estimates of order $k+1$ are obtained for the semidiscrete scheme when piecewise polynomials of degree at most $k$ ($k\geq0$) are used on overlapping uniform meshes. %Our analysis is valid for both periodic boundary conditions and %for inflow-outflow boundary conditions. We then extend the analysis to multidimensions on uniform Cartesian meshes when piecewise tensor product polynomials are used on overlapping meshes. Numerical experiments are given to demonstrate the theoretical results.

We present a discontinuous Galerkin (DG) scheme with suitable quadrature rules \cite{Tchen2017entropy} 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 \cite{godunov1972symmetric}, the entropy stable DG framework with suitable quadrature rules \cite{Tchen2017entropy}, the entropy conservative flux in \cite{chandrashekar2016entropy} inside each cell and the entropy dissipative approximate Godunov type numerical flux at cell interfaces to make the scheme entropy stable. {\color{blue}{ The main difficulty in the generalization of the results in \cite{Tchen2017entropy} is the appearance of the non-conservative ``source terms'' added in the modified MHD model introduced by Godunov \cite{godunov1972symmetric}, which do not exist in the general hyperbolic system studied in \cite{Tchen2017entropy}. 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.

Positivity-preserving discontinuous Galerkin (DG) methods for solving hyperbolic conservation laws have been extensively studied in the last several years. But nearly all the developed schemes are coupled with explicit time discretizations. Explicit discretizations suffer from the constraint for the Courant-Friedrichs-Levis (CFL) number. This makes explicit methods impractical for problems involving unstructured and extremely varying meshes or long-time simulations. Instead, implicit DG schemes are often popular in practice, especially in the computational fluid dynamics (CFD) community. In this paper we develop a high-order positivity-preserving DG method with the backward Euler time discretization for conservation laws. We focus on one spatial dimension, however the result easily generalizes to multidimensional tensor product meshes and polynomial spaces. This work is based on a generalization of the positivity-preserving limiters in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), pp.~3091--3120) and (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), pp.~8918--8934) to implicit time discretizations. Both the analysis and numerical experiments indicate that a lower bound for the CFL number is required to obtain the positivity-preserving property. The proposed scheme not only preserves the positivity of the numerical approximation without compromising the designed high-order accuracy, but also helps accelerate the convergence towards the steady-state solution and add robustness to the nonlinear solver. Numerical experiments are provided to support these conclusions.

We consider a class of time dependent second order partial differential equations governed by a decaying entropy. The solution usually corresponds to a density distribution, hence positivity (non-negativity) is expected. This class of problems covers important cases such as Fokker-Planck type equations and aggregation models, which have been studied intensively in the past decades. In this paper, we design a high order discontinuous Galerkin method for such problems. If the interaction potential is not involved, or the interaction is defined by a smooth kernel, our semi-discrete scheme admits an entropy inequality on the discrete level. Furthermore, by applying the positivity-preserving limiter, our fully discretized scheme produces non-negative solutions for all cases under a time step constraint. Our method also applies to two dimensional problems on Cartesian meshes. Numerical examples are given to confirm the high order accuracy for smooth test cases and to demonstrate the effectiveness for preserving long time asymptotics.

The main purpose of this paper is to study the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with multi-step implicit-explicit (IMEX) time discretization schemes, for solving time-dependent incompressible fluid flows. We will give theoretical analysis for the Oseen equation, and assess the performance of the schemes for incompressible Navier-Stokes equations numerically. For the Oseen equation, using first order IMEX time discretization as an example, we show that the IMEX-LDG scheme is unconditionally stable for $\Qcal_k$ elements on cartesian meshes, in the sense that the time-step $\dt$ is only required to be {bounded from above} by a positive constant independent of the spatial mesh size $h$. Furthermore, by the aid of the \emph{Stokes projection} and an elaborate energy analysis, we obtain the $L^{\infty}(L^2)$ optimal error estimates for both the velocity and the stress (gradient of velocity), in both space and time. By the \emph{inf-sup} argument, we also obtain the $L^{\infty}(L^2)$ optimal error estimates for the pressure. Numerical experiments are given to validate our main results.

In this paper we introduce a definition of the local conservation property for numerical methods solving time dependent conservation laws, which generalizes the classical local conservation definition. The motivation of our definition is the Lax-Wendroff theorem, and thus we prove it for locally conservative numerical schemes per our definition in one and two space dimensions. Several numerical methods, including continuous Galerkin methods and compact schemes, which do not fit the classical local conservation definition, are given as examples of locally conservative methods under our generalized definition.

This is an extension of our earlier work \cite{DS} in which a high order stable method was constructed for solving hyperbolic conservation laws on arbitrarily distributed point clouds. %The initial condition is given on such a point cloud, %and the algorithm solves for point values of the solution at %later time also on this point cloud. An algorithm of building a suitable polygonal mesh based on the random points was given and the traditional discontinuous Galerkin (DG) method was adopted on the constructed polygonal mesh. Numerical results in \cite{DS} show that the current scheme will generate spurious numerical oscillations when dealing with solutions containing strong shocks. In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for DG schemes on two-dimensional unstructured triangular meshes \cite{ZZSQ}, to our high order method on polygonal meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the original method in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt is particularly simple to implement and will not harm the conservativeness and compactness of the original method. Moreover, we also extend the maximum-principle-satisfying limiter for the scalar case and the positivity-preserving limiter for the Euler system to our method. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of these limiters.

For a Lagrangian scheme solving the compressible Euler equations in cylindrical coordinates, two important issues are whether the scheme can maintain spherical symmetry (symmetry-preserving) and whether the scheme can maintain positivity of density and internal energy (positivity-preserving). While there were previous results in the literature either for symmetry-preserving in the cylindrical coordinates or for positivity-preserving in cartesian coordinates, the design of a Lagrangian scheme in cylindrical coordinates, which is high order in one-dimension and second order in two-dimensions, and can maintain both spherical symmetry-preservation and positivity-preservation simultaneously, is challenging. In this paper we design such a Lagrangian scheme and provide numerical results to demonstrate its good behavior.

A new class of high order weighted essentially non-oscillatory (WENO) schemes [J. Comput. Phys., 318 (2016), 110-121] is applied to solve Euler equations with steady state solutions. It is known that the classical WENO schemes [J. Comput. Phys., 126 (1996), 202-228] 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 [J. Comput. Phys., 126 (1996), 202-228] 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.

In this article, we present a second-order in time implicit-explicit (IMEX) local discontinuous Galerkin (LDG) method for computing the Cahn-Hilliard equation, which describes the phase separation phenomenon. It is well-known that the Cahn-Hilliard equation has a nonlinear stability property, i.e., the free-energy functional decreases with respect to time. The discretized Cahn-Hilliard system modeled by the IMEX LDG method can inherit the nonlinear stability of the continuous model. We apply a stabilization technique and prove unconditional energy stability of our scheme. Numerical experiments are performed to validate the analysis. Computational efficiency can be significantly enhanced by using this IMEX LDG method with a large time step.

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.

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 (\RefOne{Jiang and Shu (1994)} \cite{jiang1994cell}) and symmetric hyperbolic systems (\RefOne{Hou and Liu (2007)} \cite{hou2007solutions}), 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 \RefOne{Carpenter et al (2014)} \cite{carpenter2014entropy} and \RefOne{Gassner (2013)} \cite{gassner2013skew}. The main ingredients are summation-by-parts (SBP) operators derived from Legendre Gauss-Lobatto quadrature, the \RefThree{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.

In this paper, we study the superconvergence behavior of discontinuous Galerkin methods using upwind numerical fluxes for one-dimensional linear hyperbolic equations with degenerate variable coefficients. The study establishes superconvergence results for the flux function approximation as well as for the DG solution itself. To be more precise, we first prove that the DG flux function is superconvergent towards a particular flux function of the exact solution, with an order of $O(h^{k+2})$, when piecewise polynomials of degree $k$ are used. We then prove that the highest superconvergence rate of the DG solution itself is $O(h^{k+\frac 32})$ as the variable coefficient degenerates or achieves the value zero in the domain. As byproducts, we obtain superconvergence properties for the DG solution and the DG flux function at special points and for cell averages. All theoretical findings are confirmed by numerical experiments.

The purpose of this work is to incorporate numerically, in a discontinuous Galerkin (DG) solver of a Boltzmann-Poisson model for hot electron transport, an electronic conduction band whose values are obtained by the spherical averaging of the full band structure given by a local empirical pseudopotential method (EPM) around a local minimum of the conduction band for silicon, as a midpoint between a radial band model and an anisotropic full band, in order to provide a more accurate physical description of the electron group velocity and conduction energy band structure in a semiconductor. This gives a better quantitative description of the transport and collision phenomena that fundamentally define the behaviour of the Boltzmann - Poisson model for electron transport used in this work. The numerical values of the derivatives of this conduction energy band, needed for the description of the electron group velocity, are obtained by means of a cubic spline interpolation. The EPM-Boltzmann-Poisson transport with this spherically averaged EPM calculated energy surface is numerically simulated and compared to the output of traditional analytic band models such as the parabolic and Kane bands, numerically implemented too, for the case of 1D $n^+-n-n^+$ silicon diodes with $400nm$ and $50nm$ channels. Quantitative differences are observed in the kinetic moments related to the conduction energy band used, such as mean velocity, average energy, and electric current (momentum).