The gas dynamics equations, coupled with a static gravitational field, admit the hydrostatic balance where the flux produced by the pressure is exactly canceled by the gravitational source term. Many astrophysical problems involve the hydrodynamical evolution in a gravitational field, therefore it is essential to correctly capture the effect of gravitational force in the simulations. Improper treatment of the gravitational force can lead to a solution which either oscillates around the equilibrium, or deviates from the equilibrium after a long time run. In this paper we design high order well-balanced finite difference WENO schemes to this system, which can preserve the hydrostatic balance state exactly and at the same time can maintain genuine high order accuracy. Numerical tests are performed to verify high order accuracy, well-balanced property, and good resolution for smooth and discontinuous solutions. The main purpose of the well-balanced schemes designed in this paper is to well resolve small perturbations of the hydrostatic balance state on coarse meshes. The more difficult issue of convergence towards such hydrostatic balance state from an arbitrary initial condition is not addressed in this paper.

The convergence to steady state solutions of the Euler equations for weighted compact nonlinear schemes (WCNS) [Deng X. and Zhang H. (2000), J. Comput. Phys. 165, 22-44 and Zhang S., Jiang S and Shu C.-W. (2008), J. Comput. Phys. 227, 7294-7321] is studied through numerical tests. Like most other shock capturing schemes, WCNS also suffers from the problem that the residue can not settle down to machine zero for the computation of the steady state solution which contains shock waves but hangs at the truncation error level. In this paper, the techniques studied in [Zhang S. and Shu. C.-W. (2007), J. Sci. Comput. 31, 273-305 and Zhang S., Jiang S and Shu. C.-W. (2011), J. Sci. Comput. 47, 216-238], to improve the convergence to steady state solutions for WENO schemes, are generalized to the WCNS. Detailed numerical studies in one and two dimensional cases are performed. Numerical tests demonstrate the effectiveness of these techniques when applied to WCNS. The residue of various order WCNS can settle down to machine zero for typical cases while the small post-shock oscillations can be removed.

We propose second order accurate discontinuous Galerkin (DG) schemes which satisfy a strict maximum principle for general nonlinear convection-diffusion equations on unstructured triangular meshes. Motivated by genuinely high order maximum-principle-satisfying DG schemes for hyperbolic conservation laws, we prove that under suitable time step restriction for forward Euler time stepping, for general nonlinear convection-diffusion equations, the same scaling limiter coupled with second order DG methods preserves the physical bounds indicated by the initial condition while maintaining uniform second order accuracy. Similar to the purely convection cases, the limiters are mass conservative and easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. Following the idea in earlier work, we extend the schemes to two-dimensional convection-diffusion equations on triangular meshes. There are no geometric constraints on the mesh such as angle acuteness. Numerical results including incompressible Navier-Stokes equations are presented to validate and demonstrate the effectiveness of the numerical methods.

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 gaseousdetonations. 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.

In this paper, we provide a priori $L^2$ error estimates for the semi-discrete discontinuous Galerkin method and the local discontinuous Galerkin method for one- and two-dimensional nonlinear Hamilton-Jacobi equations with smooth solutions. With a special Gauss-Radau projection, the optimal error estimates on rectangular meshes are obtained.

In this paper, we develop and analyze discontinuous Galerkin (DG) methods to solve hyperbolic equations involving $\delta$-singularities. Negative-order norm error estimates for the accuracy of DG approximations to $\delta$-singularities are investigated. We first consider linear hyperbolic conservation laws in one space dimension with singular initial data. We prove that, by using piecewise $k$-th degree polynomials, at time $t$, the error in the $H^{-(k+2)}$ norm over the whole domain is $(k+1/2)$-th order, and the error in the $H^{-(k+1)}(\mathbb{R}\backslash\mathcal{R}_t)$ norm is $(2k+1)$-th order, where $\mathcal{R}_t$ is the pollution region due to the initial singularity with the width of order $\mathcal{O}(h^{1/2} \log (1/h))$ and $h$ is the maximum cell length. As an application of the negative-order norm error estimates, we convolve the numerical solution with a suitable kernel which is a linear combination of B-splines, to obtain $L^2$ error estimate of $(2k+1)$-th order for the post-processed solution. Moreover, we also obtain high order superconvergence error estimates for linear hyperbolic conservation laws with singular source terms by applying Duhamel's principle. Numerical examples including an acoustic equation and the nonlinear rendez-vous algorithms are given to demonstrate the good performance of DG methods for solving hyperbolic equations involving $\delta$-singularities.

We develop an alternative formulation of conservative finite difference weighted essentially non-oscillatory (WENO) schemes to solve conservation laws. In this formulation, the WENO interpolation of the solution and its derivatives are used to directly construct the numerical flux, instead of the usual practice of reconstructing the flux functions. Even though this formulation is more expensive than the standard formulation, it does have several advantages. The first advantage is that arbitrary monotone fluxes can be used in this framework, while the traditional practice of reconstructing flux functions can be applied only to smooth flux splitting. The second advantage, which is fully explored in this paper, is that a narrower effective stencil is used compared with previous high order finite difference WENO schemes based on the reconstruction of flux functions, with a Lax-Wendroff time discretization. We will describe the scheme formulation and present numerical tests for one- and two-dimensional scalar and system conservation laws demonstrating the designed high order accuracy and non-oscillatory performance of the schemes constructed in this paper.

No abstract uploaded!

In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method detects critical numerical fluxes which may lead to negative density and pressure, and then imposes a simple flux limiter combining the high-order numerical flux with the first-order Lax-Friedrichs flux to satisfy a sufficient condition for preserving positivity. Though an extra time-step size condition is required to maintain the formal order of accuracy, it is less restrictive than those in previous works. A number of numerical examples suggest that this method, when applied on an essentially non-oscillatory scheme, can be used to prevent positivity failure when the flow involves vacuum or near vacuum and very strong discontinuities.

Wave propagation problems arise in a wide range of applications. The energy conserving property is one of the guiding principles for numerical algorithms, in order to minimize the phase or shape errors after long time integration. In this paper, we develop and analyze a local discontinuous Galerkin (LDG) method for solving the wave equation. We prove optimal error estimates, superconvergence toward a particular projection of the exact solution, and the energy conserving property for the semi-discrete formulation. The analysis is extended to the fully discrete LDG scheme, with the centered second-order time discretization (the leap-frog scheme). Our numerical experiments demonstrate optimal rates of convergence and superconvergence. We also show that the shape of the solution, after long time integration, is well preserved due to the energy conserving property.