N-waves in elasticity

We study the large time behavior of solutions of the one-dimensional equations of elasticity, typically a nonconvex system of conservation laws. We show that solutions with compact initial support converge in L₁ to a superposition of N-waves, at algebraic rate. Likewise, we show that the perturbation expansion of weakly nonlinear geometric optics is uniformly valid to second order. Our analysis uses approximate scalar laws, together with L₁ stability of scalar solutions and decay in total variation of the wave speed.     

Validity of nonlinear geometric optics for entropy solutions of multidimensional scalar conservation laws

Nonlinear geometric optics with various frequencies for entropy solutions only in L^∞ of multidimensional scalar conservation laws is analyzed. A new approach to validate nonlinear geometric optics is developed via entropy dissipation through scaling, compactness, homogenization, and L¹-stability. New multidimensional features are recognized, especially including nonlinear propagations of oscillations with high frequencies. The validity of nonlinear geometric optics for entropy solutions in L^∞ of multidimensional scalar conservation laws is justified.     

Regularized Conservation Laws and Nonlinear Geometric Optics

 This paper is devoted to the study of Cauchy problems for regularized conservation laws in Colombeau algebras of generalized functions. The existence and uniqueness of generalized solutions to these Cauchy problems are obtained. Further, we develop a generalized variant of nonlinear geometric optics for the regularized problems. Consistency with the classical results is shown to hold for scalar conservation laws with bounded variation initial data in one space variable. Received 6 November 1996; in revised form 5 August 1997     

Wave fans are special

It is shown that self-similar BV solutions of genuinely nonlinear strictly hyperbolic systems of conservation laws are special functions of bounded variation, with vanishing Cantor part.     

Well-posedness theory for geometry-compatible hyperbolic conservation laws on manifolds

Motivated by many applications (geophysical flows, general relativity), we attempt to set the foundations for a study of entropy solutions to non-linear hyperbolic conservation laws posed on a (Riemannian or Lorentzian) manifold. The flux of the conservation laws is viewed as a vector-field on the manifold and depends on the unknown function as a parameter. We introduce notions of entropy solutions in the class of bounded measurable functions and in the class of measure-valued mappings. We establish the well-posedness theory for conservation laws on a manifold, by generalizing both Kruzkov's and DiPerna's theories originally developed in the Euclidian setting. The class of geometry-compatible (as we call it) conservation laws is singled out as an important case of interest, which leads to robust L^p estimates independent of the geometry of the manifold. On the other hand, general conservation laws solely enjoy the L¹ contraction property and leads to a unique contractive semi-group of entropy solutions. Our framework allows us to construct entropy solutions on a manifold via the vanishing diffusion method or the finite volume method.     

A new nonlinear functional for general scalar conservation laws

The approach based on the construction of some nonlinear functionals was proved to be robust in the study of the well-posedness theories of hyperbolic conservation laws, especially in one space dimensional case. In particular, a generalized entropy functional was constructed in [T.-P. Liu, T. Yang, A new entropy functional for scalar conservation laws, Comm. Pure Appl. Math. 52 (1999) 1427-1442] for the L¹ stability of weak solutions. However, this generalized functional is so far only defined for scalar equations with convex flux function. In this paper, we introduce a new nonlinear functional which is motivated by the new Glimm functional introduced in [J.-L. Hua, Z.-H. Jiang, T. Yang, A new Glimm functional and convergence rate of Glimm scheme for general systems of hyperbolic conservation laws, preprint] for general scalar conservation laws. This functional improves the one given in [H.-X. Liu, T. Yang, A nonlinear functional for general scalar hyperbolic conservation laws, J. Differential Equations 235 (2) (2007) 658-667] and it can be viewed as a better attempt for the generalized entropy functional for general equations.     

A nonlinear functional for general scalar hyperbolic conservation laws

A generalized entropy functional was introduced in [T.-P. Liu, T. Yang, A new entropy functional for scalar conservation laws, Comm. Pure Appl. Math. 52 (1999) 1427-1442] for the scalar hyperbolic conservation laws with convex flux function. This functional was crucially used in the functional approach to the L¹ stability study on the system of hyperbolic conservation laws when each characteristic field is either genuinely nonlinear or linearly degenerate. However, how to construct the generalized entropy functional for scalar conservation laws with general flux, and then how to apply the functional approach to the L¹ study on general systems are still open. In this paper, we construct a new nonlinear functional which gives some partial answer to this question and we expect the analysis will shed some light on the future investigation in this direction.     

A two-dimensional hyperbolic system of nonlinear conservation laws with functional solutions

A two-dimensional hyperbolic system of nonlinear conservation laws is considered for any piecewise constant initial data having two discontinuity rays with the origin as vertex. One kind of new waves, which is labeled the Dirac-contact wave, appears in the solution. The entropy conditions for the Dirac-contact waves are given. The solutions on the Dirac-contact waves can be viewed as the bounded linear functionals onC ₀ ^∞ (R ² ×R ₊). Supported by CNSF and a grant from Academia Sinica Author’s current address: CMAP, Ecole Polytechnique, 91128 Palaiseau Cedex, France     

A non-parameterized entropy correction for Roe's approximate Riemann solver

Summary. The main drawback with Roe's approximate Riemann solver is that non-physical expansion shocks can occur in the vicinity of sonic points. Previous work aimed at enforcing the entropy condition is based on the representation of sonic rarefaction waves. We propose a new non-parameterized approach which is based on a nonlinear Hermite interpolation of an approximate flux function and the exact resolution of non convex scalar Riemann problems. Convergence and consistency with the entropy condition are proved for scalar convex conservation laws with arbitrarily large initial data. When considering strictly hyperbolic systems of conservation laws, consistency of the resulting scheme with the entropy condition is also proved for initial data sufficiently close to a constant. Numerical results on a one-dimensional shock-tube and a two-dimensional supersonic forward facing step confirm our theoretical results. Received March 1, 1993 / Revised version received August 26, 1994     

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.     

Hyperbolic conservation laws with stiff relaxation terms and entropy

We study the limiting behavior of systems of hyperbolic conservation laws with stiff relaxation terms. Reduced systems, inviscid and viscous local conservation laws, and weakly nonlinear limits are derived through asymptotic expansions. An entropy condition is introduced for N × N systems that ensures the hyperbolicity of the reduced inviscid system. The resulting characteristic speeds are shown to be interlaced with those of the original system. Moreover, the first correction to the reduced system is shown to be dissipative. A partial converse is proved for 2 × 2 systems. This structure is then applied to study the convergence to the reduced dynamics for the 2 × 2 case. © 1994 John Wiley & Sons, Inc.     

Asymptotic behavior of smooth solutions for partially dissipative hyperbolic systems with a convex entropy

We study the asymptotic time behavior of global smooth solutions to general entropy, dissipative, hyperbolic systems of balance laws in m space dimensions, under the Shizuta‐Kawashima condition. We show that these solutions approach a constant equilibrium state in the L^p‐norm at a rate O(t^{? (m/2)(1 ? 1/p)}) as t → ∞ for p ∈ [min{m, 2}, ∞]. Moreover, we can show that we can approximate, with a faster order of convergence, the conservative part of the solution in terms of the linearized hyperbolic operator for m ≥ 2, and by a parabolic equation, in the spirit of Chapman‐Enskog expansion in every space dimension. The main tool is given by a detailed analysis of the Green function for the linearized problem. © 2007 Wiley Periodicals, Inc.     

Upwind finite volume schemes for weakly coupled hyperbolic systems of conservation laws in 2D

Summary. Systems of nonlinear hyperbolic conservation laws in two space dimensions are considered which are characterized by the fact that the coupling of the equations is only due to source terms. To solve these weakly coupled systems numerically a class of explicit and implicit upwind finite volume methods on unstructured grids is presented. Provided an unique entropy solution of the system of conservation laws exists we prove that the approximations obtained by these schemes converge for vanishing discretization parameter to this entropy solution. These results are applied to examples from combustion theory and hydrology where the existence of entropy solutions can be shown. The proofs rely on an extension of a result due to DiPerna concerning measure valued solutions to the case of weakly coupled hyperbolic systems. Received April 29, 1997     

Numerical viscosity and the entropy condition

Weak solutions of hyperbolic conservation laws are not uniquely determined by their initial values; an entropy condition is needed to pick out the physically relevant solutions. The question arises whether finite difference approximations converge to this particular solution. It is known that this is not always the case with the standard Lax-Wendroff (L-W) difference scheme. In this paper a simple variant of the L-W scheme is devised which retains its desirable computational features—conservation form, three point scheme, second-order accuracy on smooth solutions, but which has the additional property that limit solutions satisfy the entropy condition. This variant is constructed by adding a simple nonlinear artificial viscosity to the usual L-W operator. The nature of the viscosity is deduced by first analyzing a model differential equation derived from the truncation error for the L-W operator, keeping only terms of order (Δx)². Furthermore, this viscosity is “switched on” only when sufficiently steep discrete gradients develop in the approximate solution: The full L-W scheme is then shown to have the desired property provided that the Courant-Friedrichs-Lewy restriction |λf′(u)|≤0.14 is satisfied.     

Irrotational approximation to the quasi-1-d gas flow

Under the assumptions that both initial data and the cross-section have sufficiently small total variation and that the initial data are supersonic (or are subsonic respectively), we prove that in any bounded domain the L¹ norm of the difference between solutions of the hyperbolic system of balance laws and the potential flow system of balance laws with the same initial data can be bounded by the cube of the sum over total variations of the initial data and the cross-section.     

Entropy flux-splittings for hyperbolic conservation laws part I: General framework

A general framework is proposed for the derivation and analysis of flux-splittings and the corresponding flux-splitting schemes for systems of conservation laws endowed with a strictly convex entropy. The approach leads to several new properties of the existing flux-splittings and to a method for the construction of entropy flux-splittings for general situations. A large family of genuine entropy flux-splittings is derived for several significant examples: the scalar conservation laws, the p-system, and the Euler system of isentropic gas dynamics. In particular, for the isentropic Euler system, we obtain a family of splittings that satisfy the entropy inequality associated with the mechanical energy. For this system, it is proved that there exists a unique genuine entropy flux-splitting that satisfies all of the entropy inequalities, which is also the unique diagonalizable splitting. This splitting can be also derived by the so-called kinetic formulation. Simple and useful difference schemes are derived from the flux-splittings for hyperbolic systems. Such entropy flux-splitting schemes are shown to satisfy a discrete cell entropy inequality. For the diagonalizable splitting schemes, an a priori L^∞ estimate is provided by applying the principle of bounded invariant regions. The convergence of entropy flux-splitting schemes is proved for the 2 × 2 systems of conservation laws and the isentropic Euler system. ©1995 John Wiley & Sons, Inc.     

Discrete approximations of BV solutions to doubly nonlinear degenerate parabolic equations

Summary. In this paper we present and analyse certain discrete approximations of solutions to scalar, doubly nonlinear degenerate, parabolic problems of the form under the very general structural condition . To mention only a few examples: the heat equation, the porous medium equation, the two-phase flow equation, hyperbolic conservation laws and equations arising from the theory of non-Newtonian fluids are all special cases of (P). Since the diffusion terms a(s) and b(s) are allowed to degenerate on intervals, shock waves will in general appear in the solutions of (P). Furthermore, weak solutions are not uniquely determined by their data. For these reasons we work within the framework of weak solutions that are of bounded variation (in space and time) and, in addition, satisfy an entropy condition. The well-posedness of the Cauchy problem (P) in this class of so-called BV entropy weak solutions follows from a work of Yin [18]. The discrete approximations are shown to converge to the unique BV entropy weak solution of (P). Received November 10, 1998 / Revised version received June 10, 1999 / Published online June 8, 2000     

Entropy stable shock capturing space–time discontinuous Galerkin schemes for systems of conservation laws

We present a streamline diffusion shock capturing spacetime discontinuous Galerkin (DG) method to approximate nonlinear systems of conservation laws in several space dimensions. The degrees of freedom are in terms of the entropy variables and the numerical flux functions are the entropy stable finite volume fluxes. We show entropy stability of the (formally) arbitrarily high order accurate method for a general system of conservation laws. Furthermore, we prove that the approximate solutions converge to the entropy measure valued solutions for nonlinear systems of conservation laws. Convergence to entropy solutions for scalar conservation laws and for linear symmetrizable systems is also shown. Numerical experiments are presented to illustrate the robustness of the proposed schemes.     

Gauss‐Green theorem for weakly differentiable vector fields,sets of finite perimeter,and balance laws

We analyze a class of weakly differentiable vector fields F : ?ⁿ → ?ⁿ with the property that F ∈ L^∞ and div F is a (signed) Radon measure. These fields are called bounded divergence‐measure fields. The primary focus of our investigation is to introduce a suitable notion of the normal trace of any divergence‐measure field F over the boundary of an arbitrary set of finite perimeter that ensures the validity of the Gauss‐Green theorem. To achieve this, we first establish a fundamental approximation theorem which states that, given a (signed) Radon measure μ that is absolutely continuous with respect to ??^{N ? 1} on ?^N, any set of finite perimeter can be approximated by a family of sets with smooth boundary essentially from the measure‐theoretic interior of the set with respect to the measure ||μ||, the total variation measure. We employ this approximation theorem to derive the normal trace of F on the boundary of any set of finite perimeter E as the limit of the normal traces of F on the boundaries of the approximate sets with smooth boundary so that the Gauss‐Green theorem for F holds on E. With these results, we analyze the Cauchy flux that is bounded by a nonnegative Radon measure over any oriented surface (i.e., an (N ? 1)‐dimensional surface that is a part of the boundary of a set of finite perimeter) and thereby develop a general mathematical formulation of the physical principle of the balance law through the Cauchy flux. Finally, we apply this framework to the derivation of systems of balance laws with measure‐valued source terms from the formulation of the balance law. This framework also allows the recovery of Cauchy entropy flux through the Lax entropy inequality for entropy solutions of hyperbolic conservation laws. © 2008 Wiley Periodicals, Inc.