Lagrangian methods are widely used in many fields for multi-material compressible flow simulations such as in astrophysics and inertial confinement fusion (ICF), due to their distinguished advantage in capturing material interfaces automatically. In some of these applications, multiple internal energy equations such as those for electron, ion and radiation are involved. In the past decades, several staggered-grid based Lagrangian schemes have been developed which are designed to solve the internal energy equation directly. These schemes can be easily extended to solve problems with multiple internal energy equations. However such schemes are typically not conservative for the total energy. Recently, significant progress has been made in developing cell-centered Lagrangian schemes which have several good properties such as conservation for all the conserved variables and easiness for remapping. However, these schemes are commonly designed to solve the Euler equations in the form of the total energy, therefore they cannot be directly applied to the solution of either the single internal energy equation or the multiple internal energy equations without significant modifications. Such modifications, if not designed carefully, may lead to the loss of some of the nice properties of the original schemes such as conservation of the total energy. In this paper, we establish an equivalency relationship between the cell-centered discretizations of the Euler equations in the forms of the total energy and of the internal energy. By a carefully designed modification in the implementation, the cell-centered Lagrangian scheme can be used to solve the compressible fluid flow with one or multiple internal energy equations and meanwhile it does not lose its total energy conservation property. An advantage of this approach is that it can be easily applied to many existing large application codes which are based on the framework of solving multiple internal energy equations. Several two dimensional numerical examples for both Euler equations and three-temperature hydrodynamic equations in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry preserving, accuracy and non-oscillatory performance.

High order path-conservative schemes have been developed for solving nonconservative hyperbolic systems. Recently, it has been observed that this approach may have some computational issues and shortcomings. In this paper, a modification to the high order path-conservative scheme in 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.

We introduce and analyze two new semi-discrete numerical methods for the multi-dimensional Vlasov-Poisson system. The schemes are constructed by combing a discontinuous Galerkin approximation to the Vlasov equation together with a mixed finite element method for the Poisson problem. We show optimal error estimates in the case of smooth compactly supported initial data. We propose a scheme that preserves the total energy of the system.

In this paper, an analysis of the superconvergence property of the semidiscrete discontinuous Galerkin (DG) method applied to one-dimensional time-dependent nonlinear scalar conservation laws is carried out. We prove that the error between the DG solution and a particular projection of the exact solution achieves $\left(k + \frac32\right)$-th order superconvergence when upwind fluxes are used. The results hold true for arbitrary nonuniform regular meshes and for piecewise polynomials of degree $k$ ($k \ge 1$), under the condition that $|f'(u)|$ possesses a uniform positive lower bound. Numerical experiments are provided to show that the superconvergence property actually holds true for nonlinear conservation laws with general flux functions, indicating that the restriction on $f(u)$ is artificial.

In this paper, we study the superconvergence of the error for the discontinuous Galerkin (DG) finite element method for linear conservation laws when upwind fluxes are used. We prove that if we apply piecewise $k$-th degree polynomials, the error between the DG solution and the exact solution is ($k+2$)-th order superconvergent at the downwind-biased Radau points with suitable initial discretization. Moreover, we also prove the DG solution is ($k+2$)-th order superconvergent both for the cell averages and for the error to a particular projection of the exact solution. The superconvergence result in this paper leads to a new {\em a posteriori} error estimate. Our analysis is valid for arbitrary regular meshes and for $\mathcal{P}^k$ polynomials with arbitrary $k\geq1$, and for both periodic boundary conditions and for initial-boundary value problems. We perform numerical experiments to demonstrate that the superconvergence rate proved in this paper is optimal.