Discrete shock profiles for systems of conservation laws

The existence of discrete shock profiles for difference schemes approximating a system of conservation laws is the major topic studied in this paper. The basic theorem established here applies to first-order accurate difference schemes; for weak shocks, this theorem provides necessary and sufficient conditions involving the truncation error of the linearized scheme which guarantee entropy satisfying or entropy violating discrete shock profiles. Several explicit difference schemes are used as examples illustrating the interplay between the entropy condition, monotonicity, and linearized stability. Entropy violating stationary shocks for second-order accurate Lax-Wendroff schemes approximating systems are also constructed. The only tools used in the proofs are local analysis and the center manifold theorem.     

Geometrical optics and post shock behavior for nonlinear conservation laws

In this paper we present an approach to the study of nonlinear waves and shocks associated with signaling problems for hyperbolic systems of conservation laws. Our approach employs a nonlinear phase variable, treats problems for which the flux function (and hence the matrix coefficients in the quasilinear system of partial differential equations) has explicit spatial dependence, and provides the post shock representation for the solution.     

Pointwise convergence to shock waves for viscous conservation laws

We are interested in the pointwise behavior of the perturbations of shock waves for viscous conservation laws. It is shown that, besides a translation of the shock waves and of linear and nonlinear diffusion waves of heat and Burgers equations, a perturbation also gives rise to algebraically decaying terms, which measure the coupling of waves of different characteristic families. Our technique is a combination of time-asymptotic expansion, construction of approximate Green functions, and analysis of nonlinear wave interactions. The pointwise estimates yield optimal L_p convergence of the perturbation to the shock and diffusion waves, 1 ≤ p ≤ ∞. The new approach of obtaining pointwise estimates based on the Green functions for the linearized system and the analysis of nonlinear wave interactions is also useful for studying the stability of waves of distinct types and nonclassical shocks. These are being explored elsewhere. © 1997 John Wiley & Sons, Inc.     

Asymptotic stability of relaxation shock profiles for hyperbolic conservation laws

This paper studies the asymptotic stability of traveling relaxation shock profiles for hyperbolic systems of conservation laws. Under a stability condition of subcharacteristic type the large time relaxation dynamics on the level of shocks is shown to be determined by the equilibrium conservation laws. The proof is due to the energy principle, using the weighted norms, the interaction of waves from various modes is treated by imposing suitable weight matrix.     

Continuum shock profiles for discrete conservation laws II: Stability

We show that the continuum shock profiles for dissipative difference schemes constructed in Part I are nonlinearly stable. It is shown first that the profiles have the conservation property, obtained as the limit of the discrete version for profiles with nearby rational, quasi‐Diophantine speeds. This allows us to formulate antidifferencing of the schemes and to apply a generalization of the pointwise approach for viscous conservation laws for the stability analysis. © 1999 John Wiley & Sons, Inc.     

Continuum shock profiles for discrete conservation laws I: Construction

We construct continuum shock profiles for finite difference schemes for hyperbolic conservation laws. Our analysis is based on the time‐asymptotic estimates for solutions of the difference schemes. Our result applies to dissipative schemes with the nonresonance property, such as the Lax‐Friedrichs and Godunov schemes. Discrete profiles have been constructed for shocks with rational speed using a fixed‐point and central‐manifold approach. For such an approach the strength of the shock is required to be small as compared to the denominator of its rational speed. Thus it does not apply to shocks with irrational speed. Our new approach yields continuum profiles for shocks with speed satisfying the Diophantine property. © 1999 John Wiley & Sons, Inc.     

Asymptotic stability of viscous shock wave profiles for viscous conservation laws

This paper is concerned with the asymptotic stability of travelling wave solution to the two-dimensional steady isentropic irrotational flow with artificial viscosity. We prove that there exists a unique travelling wave solution up to a shift to the system if the end states satisfy both the Rankine–Hugoniot condition and Lax's shock condition, and that the travelling wave solution is stable if the initial disturbance is small.     

Multiple shock fronts for hyperbolic conservation laws in higher dimensional space

The local existence of multiple shock fronts for hyperbolic conservation laws in higher dimensional space is established under the assumption that its frozen problem produces multiple uniformly stable planar shock fronts.     

Interaction of elementary waves on a boundary for a hyperbolic system of conservation laws

In this paper, we study the interaction of elementary waves including delta‐shock waves on a boundary for a hyperbolic system of conservation laws. A boundary entropy condition is derived, thanks to the results of Dubois and Le Floch (J. Differ. Equations 1988; 71 :93–122) by taking a suitable entropy–flux pair. We obtain the solutions of the initial‐boundary value problem for the system constructively, in which initial‐boundary data are piecewise constant states. Copyright © 2008 John Wiley & Sons, Ltd.     

Multiresolution schemes for conservation laws

Summary. In recent years a variety of high–order schemes for the numerical solution of conservation laws has been developed. In general, these numerical methods involve expensive flux evaluations in order to resolve discontinuities accurately. But in large parts of the flow domain the solution is smooth. Hence in these regions an unexpensive finite difference scheme suffices. In order to reduce the number of expensive flux evaluations we employ a multiresolution strategy which is similar in spirit to an approach that has been proposed by A. Harten several years ago. Concrete ingredients of this methodology have been described so far essentially for problems in a single space dimension. In order to realize such concepts for problems with several spatial dimensions and boundary fitted meshes essential deviations from previous investigations appear to be necessary though. This concerns handling the more complex interrelations of fluxes across cell interfaces, the derivation of appropriate evolution equations for multiscale representations of cell averages, stability and convergence, quantifying the compression effects by suitable adapted multiscale transformations and last but not least laying grounds for ultimately avoiding the storage of data corresponding to a full global mesh for the highest level of resolution. The objective of this paper is to develop such ingredients for any spatial dimension and block structured meshes obtained as parametric images of Cartesian grids. We conclude with some numerical results for the two–dimensional Euler equations modeling hypersonic flow around a blunt body. Received June 24, 1998 / Revised version received February 21, 2000 / Published online November 8, 2000     

Interactions of delta shock waves for a class of nonstrictly hyperbolic system of conservation laws

In this paper, we study the perturbed Riemann problem for a class of nonstrictly hyperbolic system of conservation laws, and focuse on the interactions of delta shock waves with the shock waves and the rarefaction waves. The global solutions are constructed completely with the method of splitting delta function. In solutions, we find a new kind of nonclassical wave, which is called delta contact discontinuity with Dirac delta function in both components. It is quite different from the previous ones on which only one state variable contains the Dirac delta function. Moreover, by letting perturbed parameter $\varepsilon$ tend to zero, we analyze the stability of Riemann solutions.     

Existence of discrete shock profiles of a class of monotonicity preserving schemes for conservation laws

When shock speed times is rational, the existence of solutions of shock profile equations on bounded intervals for monotonicity preserving schemes with continuous numerical flux is proved. A sufficient condition under which the above solutions can be extended to , implying the existence of discrete shock profiles of numerical schemes, is provided. A class of monotonicity preserving schemes, including all monotonicity preserving schemes with numerical flux functions, the second order upwinding flux based MUSCL scheme, the second order flux based MUSCL scheme with Lax-Friedrichs' splitting, and the Godunov scheme for scalar conservation laws are found to satisfy this condition. Thus, the existence of discrete shock profiles for these schemes is established when is rational.

     

On the convergence of a local third order shock capturing method for hyperbolic conservation laws

The Piecewise Polynomial Harmonic Method (PPHM) is a local third order accurate shock capturing method for hyperbolic conservation laws. In this paper, theoretical stability properties are presented. Using these properties, the convergence of the scheme for scalar conservation laws is obtained. A direct adaptation of these results should be used to derive the convergence of the PHM (Marquina, SIAM J. Sci. Comput. 15(4), 892–915, 1994). Finally, the method is tested in several classical problems in order to explore its numerical behavior.     

Delta shock waves with Dirac delta function in both components for systems of conservation laws

We study a class of non-strictly and weakly hyperbolic systems of conservation laws which contain the equations of geometrical optics as a prototype. The Riemann problems are constructively solved. The Riemann solutions include two kinds of interesting structures. One involves a cavitation where both state variables tend to zero forming a singularity, the other is a delta shock wave in which both state variables contain Dirac delta function simultaneously. The generalized Rankine–Hugoniot relation and entropy condition are proposed to solve the delta shock wave. Moreover, with the limiting viscosity approach, we show all of the existence, uniqueness and stability of solution involving the delta shock wave. The generalized Rankine–Hugoniot relation is also confirmed. Then our theory is successfully applied to two typical systems including the geometric optics equations. Finally, we present the numerical results coinciding with the theoretical analysis.