Weakly Nonlinear Geometric Optics for Hyperbolic Systems of Conservation Laws

We present a new approach to analyze the validation of weakly nonlinear geometric optics for entropy solutions of nonlinear hyperbolic systems of conservation laws whose eigenvalues are allowed to have constant multiplicity and corresponding characteristic fields to be linearly degenerate. The approach is based on our careful construction of more accurate auxiliary approximation to weakly nonlinear geometric optics, the properties of wave front-tracking approximate solutions, the behavior of solutions to the approximate asymptotic equations, and the standard semigroup estimates. To illustrate this approach more clearly, we focus first on the Cauchy problem for the hyperbolic systems with compact support initial data of small bounded variation and establish that the L ¹-estimate between the entropy solution and the geometric optics expansion function is bounded by O(?²), independent of the time variable. This implies that the simpler geometric optics expansion functions can be employed to study the behavior of general entropy solutions to hyperbolic systems of conservation laws. Finally, we extend the results to the case with non-compact support initial data of bounded variation.     

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.     

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.     

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.     

Existence of global L∞ solutions to a generalized n × n hyperbolic system of LeRoux type

In this article, we give the existence of global L^∞ bounded entropy solutions to the Cauchy problem of a generalized n × n hyperbolic system of LeRoux type. The main difficulty lies in establishing some compactness estimates of the viscosity solutions because the system has been generalized from 2 × 2 to n × n and more linearly degenerate characteristic fields emerged, and the emergence of singularity in the region {v₁=0} is another difficulty. We obtain the existence of the global weak solutions using the compensated compactness method coupled with the construction of entropy-entropy flux and BV estimates on viscous solutions.     

Solution of a system of nonstrictly hyperbolic conservation laws

In this paper we study a special case of the initial value problem for a 2×2 system of nonstrictly hyperbolic conservation laws studied by Lefloch, whose solution does not belong to the class ofL ^∞ functions always but may contain δ-measures as well: Lefloch's theory leaves open the possibility of nonuniqueness for some initial data. We give here a uniqueness criteria to select the entropy solution for the Riemann problem. We write the system in a matrix form and use a finite difference scheme of Lax to the initial value problem and obtain an explicit formula for the approximate solution. Then the solution of initial value problem is obtained as the limit of this approximate solution.     

Entropy solutions for linearly degenerate hyperbolic systems of rich type

Consider a linearly degenerate hyperbolic system of rich type. Assuming that each eigenvalue of the system has a constant multiplicity, we construct a representation formula of entropy solutions in L^∞ to the Cauchy problem. This formula depends on the solution of an autonomous system of ordinary differential equations taking x as parameter. We prove that for smooth initial data, the Cauchy problem for such an autonomous system admits a unique global solution. By using this formula together with classical compactness arguments, we give a very simple proof on the global existence of entropy solutions. Moreover, in a particular case of the system, we obtain an another explicit expression and the uniqueness of the entropy solution. Applications include the one-dimensional Born–Infeld system and linear Lagrangian systems.     

Numerical approximation of entropy solutions for hyperbolic integro-differential equations

Summary. Scalar hyperbolic integro-differential equations arise as models for e.g. radiating or self-gravitating fluid flow. We present finite volume schemes on unstructured grids applied to the Cauchy problem for such equations. For a rather general class of integral operators we show convergence of the approximate solutions to a possibly discontinuous entropy solution of the problem. For a specific model problem in radiative hydrodynamics we introduce a convergent fully discrete finite volume scheme. Under the assumption of sufficiently fast spatial decay of the entropy solution we can even establish the convergence rate h^1/4|ln(h)| where h denotes the grid parameter. The convergence proofs rely on appropriate variants of the classical Kruzhkov method for local balance laws together with a truncation technique to cope with the nonlocal character of the integral operator.Mathematics Subject Classification (2000): 35L65, 35Q35, 65M15     

A splitting positive definite mixed finite element method for two classes of integro-differential equations

In this article, we establish a new mixed finite element procedure to solve the second-order hyperbolic and pseudo-hyperbolic integro-differential equations, in which the mixed element system is symmetric positive definite without requiring the LBB consistency condition. Convergence analysis shows that the method yields the approximate solutions with optimal accuracy in L ²(??) norm for u and in H(?div;??) norm for the flux???. Numerical experiments are given to verify the theoretical results.     

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.     

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.     

Large time behavior of entropy solutions to some hyperbolic system with dissipative structure

In this paper the large time behavior of the global L∞ entropy solutions to the hyperbolic system with dissipative structure is investigated. It is proved that as t →∞ the entropy solutions tend to a constant equilibrium state in L2 norm with exponential decay even when the initial values are arbitrarily large. As an illustration, a class of 2 × 2 system is studied.     

Global existence and asymptotic behavior of classical solutions of quasilinear hyperbolic systems with linearly degenerate characteristic fields

In this paper, we study the global existence and the asymptotic behavior of classical solution of the Cauchy problem for quasilinear hyperbolic system with constant multiple and linearly degenerate characteristic fields. We prove that the global C¹ solution exists uniquely if the BV norm of the initial data is sufficiently small. Based on the existence result on the global classical solution, we show that, when the time t tends to the infinity, the solution approaches a combination of C¹ traveling wave solutions. Finally, we give an application to the equation for time-like extremal surfaces in the Minkowski space-time R1+n.     

Approximation of entropy on hyperbolic sets for one-dimensional maps and their multidimensional perturbations

We consider piecewise monotone (not necessarily, strictly) piecewise C ² maps on the interval with positive topological entropy. For such a map f we prove that its topological entropy h _top(f) can be approximated (with any required accuracy) by restriction on a compact strictly f-invariant hyperbolic set disjoint from some neighborhood of prescribed set consisting of periodic attractors, nonhyperbolic intervals and endpoints of monotonicity intervals. By using this result we are able to generalize main theorem from [1] on chaotic behavior of multidimensional perturbations of solutions for difference equations which depend on two variables at nonperturbed value of parameter.     

Asymptotic behavior of classical solutions of reducible quasilinear hyperbolic systems with characteristic boundaries

In this paper, we investigate the asymptotic behavior of classical solutions of reducible quasilinear hyperbolic systems with characteristic boundaries. Under some suitable assumptions, we prove that the solution approaches a combination of Lipschitz continuous and piecewise C¹ traveling wave solution. As an application, we apply the result to the equation for time-like extremal surfaces in the Minkowski space-time R1+(1+n).     

Hyperbolic conservation laws on manifolds. An error estimate for finite volume schemes

Following Ben-Artzi and LeFloch, we consider nonlinear hyperbolic conservation laws posed on a Riemannian manifold, and we establish an L1-error estimate for a class of finite volume schemes allowing for the approximation of entropy solutions to the initial value problem. The error in the L1 norm is of order h1/4 at most, where h represents the maximal diameter of elements in the family of geodesic triangulations. The proof relies on a suitable generalization of Cockburn, Coquel, and LeFloch＇s theory which was originally developed in the Euclidian setting. We extend the arguments to curved manifolds, by taking into account the effects to the geometry and overcoming several new technical difficulties.     

On a time-splitting method for a scalar conservation law with a multiplicative stochastic perturbation and numerical experiments

In this paper, we present a numerical scheme for a first-order hyperbolic equation of nonlinear type perturbed by a multiplicative noise. The problem is set in a bounded domain D of ${\mathbb{R}^{d}}$ and with homogeneous Dirichlet boundary condition. Using a time-splitting method, we are able to show the existence of an approximate solution. The result of convergence of such a sequence is based on the work of Bauzet–Vallet–Wittbold (J Funct Anal, 2013), where the authors used the concept of measure-valued solution and Kruzhkov’s entropy formulation to show the existence and uniqueness of the stochastic weak entropy solution. Then, we propose numerical experiments by applying this scheme to the stochastic Burgers’ equation in the one-dimensional case.     

Asymptotic behavior of the solutions to the damped compressible Euler equations with vacuum

The asymptotic behavior of solutions of the damped compressible Euler equations is conjectured to obey to the famous porous media equations (PMES). The previous works on this topic concern the case away from vacuum where the system is strictly hyperbolic. In present paper, we prove that the L^∞ entropy weak solution with vacuum, obtained by the compensated compactness theory, converges strongly in space to the unique similarity solution of the related PME, as time goes to infinity.     

Second‐order finite‐volume schemes for a non‐linear hyperbolic equation: error estimate

We study second‐order finite‐volume schemes for the non‐linear hyperbolic equation u_t(x, t) + div F(x, t, u(x, t)) = 0 with initial condition u₀. The main result is the error estimate between the approximate solution given by the scheme and the entropy solution. It is based on some stability properties verified by the scheme and on a discrete entropy inequality. If u₀ ∈ L^∞ ∩ BV_loc(ℝ^N), we get an error estimate of order h^1/4, where h defines the size of the mesh. Copyright © 2000 John Wiley & Sons, Ltd.     

Shallow water equations: viscous solutions and inviscid limit

We establish the inviscid limit of the viscous shallow water equations to the Saint-Venant system. For the viscous equations, the viscosity terms are more degenerate when the shallow water is close to the bottom, in comparison with the classical Navier-Stokes equations for barotropic gases; thus, the analysis in our earlier work for the classical Navier-Stokes equations does not apply directly, which require new estimates to deal with the additional degeneracy. We first introduce a notion of entropy solutions to the viscous shallow water equations and develop an approach to establish the global existence of such solutions and their uniform energy-type estimates with respect to the viscosity coefficient. These uniform estimates yield the existence of measure-valued solutions to the Saint-Venant system generated by the viscous solutions. Based on the uniform energy-type estimates and the features of the Saint-Venant system, we further establish that the entropy dissipation measures of the viscous solutions for weak entropy-entropy flux pairs, generated by compactly supported C ² test-functions, are confined in a compact set in H ^?1, which yields that the measure-valued solutions are confined by the Tartar-Murat commutator relation. Then, the reduction theorem established in Chen and Perepelitsa [5] for the measure-valued solutions with unbounded support leads to the convergence of the viscous solutions to a finite-energy entropy solution of the Saint-Venant system with finite-energy initial data, which is relative with respect to the different end-states of the bottom topography of the shallow water at infinity. The analysis also applies to the inviscid limit problem for the Saint-Venant system in the presence of friction.