Based on the double shockwave approximation procedure, a local Riemann solver for strongly nonlinear equations of state (EOS) such as the Jones-Wilkins-Lee (JWL) EOS is presented, which has suppressed successfully numerical oscillation caused by high-density ratio and high-pressure ratio across the interface between explosion products and air. The real ghost fluid method (RGFM) and the level set method have been used for converting multi-medium flows into pure flows and for implicitly tracking the interface, respectively. A fifth order finite difference weighted essentially non-oscillatory (WENO) scheme and a third order TVD Runge-Kutta method are utilized for the spatial discretization and the time advance, respectively. An enclosed-type MPI-based parallel methodology for the RGFM procedure on a uniform structured mesh is presented to realize the parallelization of three-dimensional {\color{red}(3D)} air explosion. The overall process of 3D air explosion of both TNT and aluminized explosives has been successfully simulated.

In this paper, 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. We show that, with ${P}^k$ (piecewise polynomials of degree $k$) finite elements in 1D problems, if the quadrature over elements is exact for polynomials of degree $(2k)$, error estimates of $O(h^{k+1/2})$ are obtained for general monotone fluxes, and optimal estimates of $O(h^{k+1})$ are obtained for upwind fluxes. For multidimensional problems, if in addition quadrature over edges is exact for polynomials of degree $(2k+1)$, error estimates of $O(h^k)$ are obtained for general monotone fluxes, and $O(h^{k+1/2})$ are obtained for monotone and sufficiently smooth numerical fluxes. Numerical results validate our analysis.

In this paper we generalize a new type of compact Hermite weighted essentially non-oscillatory (HWENO) limiter for the Runge-Kutta discontinuous Galerkin (RKDG) method, which was recently developed for structured meshes, to two dimensional unstructured meshes. The main idea of this HWENO limiter is to reconstruct the new polynomial by the usage of the entire polynomials of the DG solution from the target cell and its neighboring cells in a least squares fashion while maintaining the conservative property, then use the classical WENO methodology to form a convex combination of these reconstructed polynomials based on the smoothness indicators and associated nonlinear weights. The main advantage of this new HWENO limiter is the robustness for very strong shocks and simplicity in implementation especially for the unstructured meshes considered in this paper, since only information from the target cell and its immediate neighbors is needed. Numerical results for both scalar and system equations are provided to test and verify the good performance of this new limiter.

Finite difference WENO schemes have established themselves as very worthy performers for entire classes of applications that involve hyperbolic conservation laws. In this paper we report on two major advances that make finite difference WENO schemes more efficient.

H(div) conforming and discontinuous Galerkin (DG) methods are designed for incompressible Euler's equation in two and three dimension. Error estimates are proved for both the semi-discrete method and fully-discrete method using backward Euler time stepping. Numerical examples exhibiting the performance of the methods are given.