A new variation for the relativistic Euler equations

The Glimm scheme is one of the so famous techniques for getting solutions of the general initial value problem by building a convergent sequence of approximate solutions. The approximation scheme is based on the solution of the Riemann problem. In this paper, we use a new strength function in order to present a new kind of total variation of a solution. Based on this new variation, we use the Glimm scheme to prove the global existence of weak solutions for the nonlinear ultra-relativistic Euler equations for a class of large initial data that involve the interaction of nonlinear waves.

Global bounded variation solutions describing Fanno–Rayleigh fluid flows in nozzles

In this paper, we investigate the initial-boundary value problem of compressible Euler equations including friction and heating that model the transonic Fanno–Rayleigh flows through symmetric variable area nozzles. In particular, the case of contracting nozzles is considered. A new version of a generalized Glimm scheme (GGS) is presented for establishing the global existence of entropy solutions with bounded variation. Modified Riemann and boundary Riemann solutions are applied to design this GGS, which is constructed using the contraction matrices acting on the homogeneous Riemann (or boundary-Riemann) solutions. The extended Glimm–Goodman’s type of wave interaction estimates are investigated to determine the stability of the scheme and the positivity of gas velocity that results in the existence of the weak solution. The limit of approximation solutions serves as an entropy solution. Moreover, a quantitative relation between the shape of the nozzle, friction, and heat is proposed for the global existence result in the contracting nozzle. Numerical simulations of the contraction-expansion and expansion-contraction nozzles are presented to validate the scheme.

A remark on the Glimm scheme for inhomogeneous hyperbolic systems of balance laws

General hyperbolic systems of balance laws with inhomogeneous flux and source are studied. Global existence of entropy weak solutions to the Cauchy problem is established for small BV data under appropriate assumptions on the decay of the flux and the source with respect to space and time. There is neither a hypothesis about equilibrium solution nor about the dependence of the source on the state vector as previous results have assumed.

On a quadratic functional for scalar conservation laws

We prove a quadratic interaction estimate for approximate solutions to scalar conservation laws obtained by the wavefront tracking approximation or the Glimm scheme. This quadratic estimate has been used in the literature to prove the convergence rate of the Glimm scheme. The proof is based on the introduction of a quadratic functional 𝔔(t), decreasing at every interaction, and such that its total variation in time is bounded. Differently from other interaction potentials present in the literature, the form of this functional is the natural extension of the original Glimm functional, and coincides with it in the genuinely nonlinear case.

Numerical Approximation of a Compressible Multiphase System

The numerical simulation of non conservative system is a difficult challenge for two reasons at least. The first one is that it is not possible to derive jump relations directly from conservation principles, so that in general, if the model description is non ambiguous for smooth solutions, this is no longer the case for discontinuous solutions. From the numerical view point, this leads to the following situation: if a scheme is stable, its limit for mesh convergence will depend on its dissipative structure. This is well known since at least [1]. In this paper we are interested in the “dual” problem: given a system in non conservative form and consistent jump relations, how can we construct a numerical scheme that will, for mesh convergence, provide limit solutions that are the exact solution of the problem. In order to investigate this problem, we consider a multiphase flow model for which jump relations are known. Our scheme is an hybridation of Glimm scheme and Roe scheme.