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.

Difference approximations of hyperbolic partial differential equations with highly oscillatory coefficients and initial values are studied. Analysis of strong and weak convergence is carried out in the practically interesting case when the discretization step sizes are essentially independent of the oscillatory wave lengths. © 1993 John Wiley & Sons, Inc.

The nonlinear stability in the {L^p} -norm, p ≥ 1 , of stationary weak discrete shocks for the Lax-Friedrichs scheme approximating general m × m systems of nonlinear hyperbolic conservation laws is proved, provided that the summations of the initial perturbations equal zero. The result is proved by using both a weighted estimate and characteristic energy method based on the internal structures of the discrete shocks and the essential monotonicity of the Lax-Friedrichs scheme.

In this paper we study the asymptotic nonlinear stability of discrete shocks for the Lax-Friedrichs scheme for approximating general m×m systems of nonlinear hyperbolic conservation laws. It is shown that weak single discrete shocks for such a scheme are nonlinearly stable in the Lp-norm for all p ≧ 1, provided that the sums of the initial perturbations equal zero. These results should shed light on the convergence of the numerical solution constructed by the Lax-Friedrichs scheme for the single-shock solution of system of hyperbolic conservation laws. If the Riemann solution corresponding to the given far-field states is a superposition of m single shocks from each characteristic family, we show that the corresponding multiple discrete shocks are nonlinearly stable in Lp (P ≧ 2). These results are proved by using both a weighted estimate and a characteristic energy method based on the internal structures of the discrete shocks and the essential monotonicity of the Lax-Friedrichs scheme.

This is the first of a series of papers on the subject of projection methods for viscous incompressible flow calculations. The purpose of these papers is to provide a thorough understanding of the numerical phenomena involved in the projection methods, particularly when boundaries are present, and point to ways of designing more efficient, robust, and accurate numerical methods based on the primitive variable formulation. This paper contains the following topics: 1. convergence and optimal error estimates for both velocity and pressure up to the boundary; 2. explicit characterization of the numerical boundary layers in the pressure approximations and the intermediate velocity fields; 3. the effect of choosing different numerical boundary conditions at the projection step. We will show that a different choice of boundary conditions gives rise to different boundary layer structures. In particular, the straightforward Dirichlet boundary condition for the pressure leads to O(1) numerical boundary layers in the pressure and deteriorates the accuracy in the interior; and 4. postprocessing the numerical solutions to get more accurate approximations for the pressure.