Convergence of a second-order shock-capturing scheme for the system of isentropic gas dynamics with $ {L^\infty }$ initial data is established by analyzing the entropy dissipation measures. This scheme is modified from the classical MUSCL scheme to treat the vacuum problem in gas fluids and to capture local entropy near shock waves. Convergence of this scheme for the piston problem is also discussed.

Simulation of electromagnetic wave propagation in metamaterials leads to more complicated time domain Maxwell's equations than the standard Maxwell's equations in free space. In this paper, we develop and analyze a non-dissipative discontinuous Galerkin (DG) method for solving the Maxwell's equations in Drude metamaterials. Previous discontinuous Galerkin methods in the literature for electromagnetic wave propagation in metamaterials were either non-dissipative but sub-optimal, or dissipative and optimal. Our method uses a different and simple choice of numerical fluxes, achieving provable non-dissipative stability and optimal error estimates simultaneously. We prove the stability and optimal error estimates for both semi- and fully discrete DG schemes, with the leap-frog time discretization for the fully discrete case. Numerical results are given to demonstrate that the DG method can solve metamaterial Maxwell's equations effectively.

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high

Wormhole propagation, arising in petroleum engineering, is used to describe the distribution of acid and the increase of porosity in carbonate reservoir under dissolution of injected acid. The important physical features of porosity and acid concentration include their boundedness between $0$ and $1$, as well as the monotone increasing for porosity. How to keep these properties in the simulation is crucial to the robustness of the numerical algorithm. In this paper, we propose high-order bound-preserving discontinuous Galerkin methods to keep these important physical properties. The main technique is to introduce a new variable $r$ to replace the original acid concentration and use a consistent flux pair to deduce a ghost equation such that the positive-preserving technique can be applied on both original and deduced equations. A high-order slope limiter is used to keep a polynomial upper bound which changes over time for $r$. Moreover, the high-order accuracy is attained by the flux limiter. Numerical examples are given to demonstrate the high-order accuracy and bound-preserving property of the numerical technique.

Local discontinuous Galerkin (LDG) methods are popular for convection-diffusion equations. In LDG methods, we introduce an auxiliary variable $p$ to represent the derivative of the primary variable $u$, and solve them on the same mesh. In this paper, we will introduce a new LDG method, and solve $u$ and $p$ on different meshes. The stability and error estimates will be investigated. The new algorithm is more flexible and flux-free for pure diffusion equations without introducing additional computational cost compared with the original LDG methods, since it is not necessary to solve each equation twice. Moreover, it is possible to construct third-order maximum-principle-preserving schemes based on the new algorithm. However, one cannot anticipate optimal accuracy in some special cases. In this paper, we will find out the reason for accuracy degeneration which further leads to several alternatives to obtain optimal convergence rates. Finally, several numerical experiments will be given to demonstrate the stability and optimal accuracy of the new algorithm.