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.

We prove the convergence of vortex blob methods to classical weak solutions for the twodimensional incompressible Euler equations with initial data satisfying the conditions that the vorticity is a finite Radon measure of distinguished sign and the kinetic energy is locally bounded. This includes the important example of vortex sheets. The result is valid as long as the computational grid size h does not exceed the smoothing blob size E, i.e., h/~ 5 C. 0 1995 John Wiley & Sons, Inc.

The development of shocks in nonlinear hyperbolic conservation laws may be regularized through either diffusion or relaxation. However, we have observed surprisingly that for some physical problems, when both of the smoothing factors - diffusion and relaxation - coexist, under appropriate asymptotic assumptions, the dispersive waves are enhanced. This phenomenon is studied asymptotically in the sense of the Chapman-Enskog expansion and demonstrated numerically.