Euler equations for isentropic gas dynamics with general pressure law

In this work, we explore the limiting behavior of Riemann solutions to the Euler equations in isentropic gas dynamics with general pressure law. We demonstrate that in the distributional sense the delta wave of zero-pressure gas dynamics is formed by a limit solution. Finally, to present the concentration phenomena, we also offer some numerical outcomes.


Introduction
The primary theme of the paper is the Euler equations of isentropic gas dynamics in Eulerian coordinates, with Riemann initial values (γ , v)| τ =0 = ⎧ ⎨ ⎩ (γ -, v -), η < 0, (γ + , v + ), η > 0, (2) where v, p(γ ), and γ are the velocity, pressure, and density, respectively, We assume that v -> v + and the pressure function is with p satisfying ⎧ ⎨ ⎩ p (γ ) > 0, Equations (1) can be expressed by the hyperbolic system for conservation laws: As → 0, (1) becomes the zero-pressure limiting scheme Many scientists have researched the alternatives to the zero-pressure gas dynamics [1][2][3][4][5][6][7][8]. Brenier and Grenier [9] provided uncountable viscous particles possessing momentum and mass preserved under collision to enhance the dynamics of one-dimensional pressureless gas (4). Bouchut [1] proved the presence of the system entropy solution (6) for the Riemann case, and the associated numerical findings affirmed that the alternatives have the concentration of the density, that is, the delta wave. Grenier and Brenier [9] and Rykov and Sinai [2] initially obtained work delta wave entropy alternatives for the Cauchy problem irrespective of the pressureless gas; see Huang, Ding, and Wang [7] with distinct original information classes. Wang and Ding [8] first demonstrated the uniqueness when the original information were features. With a distinct technique, Bouchut [10] had the same outcome. Finally, Wang and Huang [5,11,12] demonstrated the uniqueness when measuring radon in the original information under circumstances of energy and Oleinik entropy. See Huang [13][14][15] for overall pressure-free system. Regarding the two-dimensional Riemann case of pressureless gas dynamics, see [6,[16][17][18], and the references therein.
Let us now turn to the Euler scheme isentropic (1). The Euler system officially refers as a dynamics of zero gas pressure (6) in a situation whereby pressure gets close to zero or nonzero. Chen and Lin [19] initially demonstrated the development of vacuum states and δ-shocks of the Riemann solutions to the Euler equations (1) as the stress comes to zero, which describes the concentration and cavitation phenomenon exactly in mathematics. We demonstrated the same concentration phenomenon in [20][21][22] as the stress tends to a nonzero constant. All operations of [19] are carried out only for the gas dynamics of γ -law. In current research, we demonstrate the same concentration phenomenon in a general pressure mode for P = p(γ ) as → 0, where the pressure p satisfies (4). To the best of our knowledge, the scheme Riemann problem (1) was not previously considered in the literature as → 0.
The rest of the work is structured as follows. In Sect. 2, we report some preliminaries on the delta wave solution for the pressure-free gas dynamical system (6). In Sect. 3, we establish the classical Riemann solutions to the Euler equations. In Sect. 4, we depict the limiting behavior of the Riemann solution to the Euler equations (1) in the forming of delta wave. Section 5 reports the Lax-Friedrichs numerical scheme [23,24] of first order.

Preliminaries
Alternatives to the pressure-free method have been explored in (6) involving delta wave, backed on a sleek curve, [6], the solution is a delta function.
Suppose that S = {(η(s), τ (s)) : c < s < d} is a smooth curve, we define for φ ∈ C ∞ 0 (R 2 + ). The delta function weight is denoted by w(·). For the Riemann case with v -> v + , the Dirac-measured solution with parameter σ can be expressed by For across, the discontinuity [h] := h +hrepresents the function jump. The Diracmeasured alternative (γ , γ v) established above is said to be a "delta wave" to the pressureless gas dynamics (6) if for all ψ ∈ C ∞ 0 (R 2 + ). The insight into the distribution is taken as Thus σ is determined in a unique way as under the entropy axiom v -> σ > v + [6].
Remark 2.1 In a geometrical sense the entropy axiom v -> σ > v + implies that any characteristic lines generated in the x-axis enter into the Delta wave η = σ τ ; see Fig. 1. The results give us the plethora of particles concentrated on η = σ τ with the increase in time τ . Thus by (10) we get [6,13].

Riemann solutions
Next, we can take the Riemann case (1) and (2) in the scenario v -> v + . The eigenvalues of system (1) are There should be a general pressure under assumption (4). We can take a particular smooth solution along the line S = {(σ τ , τ ) : 0 < τ < ∞} with bounded jump. Then for the discontinuity of system (1), by the Rankine-Hugoniot relation we have so that we obtain with σ , (γ -, v -), and (γ + , v + ) showing the shock speed and the left and right states, respectively.
The Lax entropy inequality [25] yields So by applying the Rankine-Hugoniot condition, we arrive at a contradiction to (15). Therefore we obtain the shock curve S 1 (γ -, v -) in the phase plane for the first family (see Fig. 2 In this way, we can obtain the shock curve S 2 (γ -, v -) for the second family: We are especially interested in the case with the shock speeda respectively. The RIE solution in this regard is

Formation of delta wave
Assume that v -> v + for the limiting manner of Riemann solutions for the Euler system of the isentropic fluids with initial data as tends to 0. So we have to show that the limit is the delta wave solution to the pressureless gas dynamics (6).

Lemma 4.2 lim
In the right-hand side of (23), taking the limit, we obtain which leads to a = ( Moreover, we have

Proposition 4.1 We have
and Proof By simple calculation we get Remark 4. 1 The results show that the delta wave velocity σ is the limiting value of the speed particle v * and coupled shocks velocity σ 1 , σ 2 . We can illustrate it as the catching up of 1-shock with the coupled shock, whereas overlap of the coupled shocks forms a delta wave. In addition, multiplying (26) with τ , we may readily notice a mass concentration situation; see also Remark 2.1.

Numerical results
Here we show a set of representative numerical outcomes for the Euler scheme (1) with original Riemann information (2). A number of iterative numerical tests have been conducted to verify that the proof does not coincide with numerical objects. The initial values in (1)-(2) are solved by p(γ ) = γ 2 + γ 3 . This model of pressure satisfies general assumption (4). We take the following initial data: For (5) to be discretized, we apply the following explicit conservative approach: where η is the spatial size, τ is time step, and f n j+ 1 2 is the numerical flux. For the simplicity, we use the local first-order Lax-Friedrichs scheme [23]: with c = p (γ ) denoting the velocity sound. Using the stability axiom, we got the time where CFL is the short form of the Courant-Friedrichs-Lewyor number [21]. As demonstrated in Figs. 3-4, the intermediate state of the density approaches infinity, and the velocity becomes constant as is close to zero. As the limit is applied, the overlap of the coupled shocks generates the delta wave, which solves the zeropressure gas dynamics (6), (2). Thus we conclude that the infinite amplitude is retained by the solutions if the strict hyperbolicity degenerates. In this situation the system of equations (1) is a model. For the numerous selection of (i.e., = 0.1, 0.007, 0.0015, 0.00009, and the time τ = 0.2), the numerical simulation is done by taking CFL = 1 with 200 points and is shown in Figs. 4-5. The density producing a Dirac measure is depicted on the left sides of Figs. 4-5. Similarly, the right sides of Figs. 5-6 present the changes in velocity as decreasing using the space variable η as horizontal and velocity v as vertical, which guarantee a step function in the following. It is illustrated that the delta wave for the zero gas dynamics (6) appears in a way of an efficient numerical approach for the Euler system by considering the limit of pressureless in the solutions. In addition to this technique, we can explore a lot of attractive physical phenomena by employing the proposed approach to evaluate solutions for different types of problems.

Conclusion
Finally, we conclude that the limiting values of Riemann solutions to the Euler equations are derived in isentropic gas dynamics fulfilling the general pressure law. We have also provided the zero-pressure gas dynamics of delta wave in distributional form made by the limiting values of the required solution. We have given some numerical approximations for the phenomena of concentration representation. As a future suggestion, the stability of the vacuum state and the case of nonisentropic Euler equations can be analyzed and investigated.