Dispersive and diffusive-dispersive shock waves for nonconvex conservation laws

We consider two physically and mathematically distinct regularization mechanisms of scalar hyperbolic conservation laws. When the flux is convex, the combination of diffusion and dispersion are known to give rise to monotonic and oscillatory traveling waves that approximate shock waves. The zero-diffusion limits of these traveling waves are dynamically expanding dispersive shock waves (DSWs). A richer set of wave solutions can be found when the flux is non-convex. This review compares the structure of solutions of Riemann problems for a conservation law with non-convex, cubic flux regularized by two different mechanisms: 1) dispersion in the modified Korteweg--de Vries (mKdV) equation; and 2) a combination of diffusion and dispersion in the mKdV-Burgers equation. In the first case, the possible dynamics involve two qualitatively different types of DSWs, rarefaction waves (RWs) and kinks (monotonic fronts). In the second case, in addition to RWs, there are traveling wave solutions approximating both classical (Lax) and non-classical (undercompressive) shock waves. Despite the singular nature of the zero-diffusion limit and rather differing analytical approaches employed in the descriptions of dispersive and diffusive-dispersive regularization, the resulting comparison of the two cases reveals a number of striking parallels. In contrast to the case of convex flux, the mKdVB to mKdV mapping is not one-to-one. The mKdV kink solution is identified as an undercompressive DSW. Other prominent features, such as shock-rarefactions, also find their purely dispersive counterparts involving special contact DSWs, which exhibit features analogous to contact discontinuities. This review describes an important link between two major areas of applied mathematics, hyperbolic conservation laws and nonlinear dispersive waves.Comment: Revision from v2; 57 pages, 19 figure

Multidimensional Conservation Laws: Overview, Problems, and Perspective

Some of recent important developments are overviewed, several longstanding open problems are discussed, and a perspective is presented for the mathematical theory of multidimensional conservation laws. Some basic features and phenomena of multidimensional hyperbolic conservation laws are revealed, and some samples of multidimensional systems/models and related important problems are presented and analyzed with emphasis on the prototypes that have been solved or may be expected to be solved rigorously at least for some cases. In particular, multidimensional steady supersonic problems and transonic problems, shock reflection-diffraction problems, and related effective nonlinear approaches are analyzed. A theory of divergence-measure vector fields and related analytical frameworks for the analysis of entropy solutions are discussed.Comment: 43 pages, 3 figure

Unique continuation property with partial information for two-dimensional anisotropic elasticity systems

In this paper, we establish a novel unique continuation property for two-dimensional anisotropic elasticity systems with partial information. More precisely, given a homogeneous elasticity system in a domain, we investigate the unique continuation by assuming only the vanishing of one component of the solution in a subdomain. Using the corresponding Riemann function, we prove that the solution vanishes in the whole domain provided that the other component vanishes at one point up to its second derivatives. Further, we construct several examples showing the possibility of further reducing the additional information of the other component. This result possesses remarkable significance in both theoretical and practical aspects because the required data is almost halved for the unique determination of the whole solution.Comment: 14 pages, 1 figur

Coupling techniques for nonlinear hyperbolic equations. IV. Multi-component coupling and multidimensional well-balanced schemes

This series of papers is devoted to the formulation and the approximation of coupling problems for nonlinear hyperbolic equations. The coupling across an interface in the physical space is formulated in term of an augmented system of partial differential equations. In an earlier work, this strategy allowed us to develop a regularization method based on a thick interface model in one space variable. In the present paper, we significantly extend this framework and, in addition, encompass equations in several space variables. This new formulation includes the coupling of several distinct conservation laws and allows for a possible covering in space. Our main contributions are, on one hand, the design and analysis of a well-balanced finite volume method on general triangulations and, on the other hand, a proof of convergence of this method toward entropy solutions, extending Coquel, Cockburn, and LeFloch's theory (restricted to a single conservation law without coupling). The core of our analysis is, first, the derivation of entropy inequalities as well as a discrete entropy dissipation estimate and, second, a proof of convergence toward the entropy solution of the coupling problem.Comment: 37 page

A Positive and Entropy-Satisfying Finite Volume Scheme for the Baer-Nunziato Model

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer-Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in [16] for the isentropic Baer-Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer-Nunziato model, namely Schwendeman-Wahle-Kapila's Godunov-type scheme [39] and Toro-Tokareva's HLLC scheme [42]. The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions

Bi-partite entanglement entropy in integrable models with backscattering

In this paper we generalise the main result of a recent work by J. L. Cardy and the present authors concerning the bi-partite entanglement entropy between a connected region and its complement. There the expression of the leading order correction to saturation in the large distance regime was obtained for integrable quantum field theories possessing diagonal scattering matrices. It was observed to depend only on the mass spectrum of the model and not on the specific structure of the diagonal scattering matrix. Here we extend that result to integrable models with backscattering (i.e. with non-diagonal scattering matrices). We use again the replica method, which connects the entanglement entropy to partition functions on Riemann surfaces with two branch points. Our main conclusion is that the mentioned infrared correction takes exactly the same form for theories with and without backscattering. In order to give further support to this result, we provide a detailed analysis in the sine-Gordon model in the coupling regime in which no bound states (breathers) occur. As a consequence, we obtain the leading correction to the sine-Gordon partition function on a Riemann surface in the large distance regime. Observations are made concerning the limit of large number of sheets.Comment: 22 pages, 2 figure

Riemann solvers and undercompressive shocks of convex FPU chains

We consider FPU-type atomic chains with general convex potentials. The naive continuum limit in the hyperbolic space-time scaling is the p-system of mass and momentum conservation. We systematically compare Riemann solutions to the p-system with numerical solutions to discrete Riemann problems in FPU chains, and argue that the latter can be described by modified p-system Riemann solvers. We allow the flux to have a turning point, and observe a third type of elementary wave (conservative shocks) in the atomistic simulations. These waves are heteroclinic travelling waves and correspond to non-classical, undercompressive shocks of the p-system. We analyse such shocks for fluxes with one or more turning points. Depending on the convexity properties of the flux we propose FPU-Riemann solvers. Our numerical simulations confirm that Lax-shocks are replaced by so called dispersive shocks. For convex-concave flux we provide numerical evidence that convex FPU chains follow the p-system in generating conservative shocks that are supersonic. For concave-convex flux, however, the conservative shocks of the p-system are subsonic and do not appear in FPU-Riemann solutions

A note on regularity and failure of regularity for systems of conservation laws via Lagrangian formulation

The paper recalls two of the regularity results for Burgers' equation, and discusses what happens in the case of genuinely nonlinear, strictly hyperbolic systems of conservation laws. The first regularity result which is considered is Oleinik-Ambroso-De Lellis SBV estimate: it provides bounds on the x-derivative of u when u is an entropy solution of the Cauchy problem for Burgers' equation with bounded initial data. Its extensions to the case of systems is then mentioned. The second regularity result of debate is Schaeffer's theorem: entropy solutions to Burgers' equation with smooth and generic, in a Baire category sense, initial data are piecewise smooth. The failure of the same regularity for general genuinely nonlinear systems is next described. The main focus of this paper is indeed including heuristically an original counterexample where a kind of stability of a shock pattern made by infinitely many shocks shows up, referring to [Caravenna-Spinolo] for the rigorous result.Comment: 10 pages, 1 figur