Isothermal Limit of Entropy Solutions for the Euler Equations with Source Term

We investigate the isothermal limit of entropy solutions for isentropic Euler equations with a special source term. It is shown that the

L

∞

entropy solutions of the isentropic Euler equations converge strongly to the corresponding entropy solutions of the isothermal Euler equations as the adiabatic exponent

γ→1

. The proof is based on the homogeneous case which is achieved by combining the entropy analysis and refined kinetic formulation with compensated compactness framework. We apply a maximum principle to obtain the uniform estimates with respect to

γ

for both density and momentum.

• Entropy Solutions
• Isothermal Limit
• Isentropic Euler Equations
• Isothermal Euler Equations
• Strong Convergence

1. Introduction

In this paper, we consider the one-dimensional compressible Euler equations for isentropic gas dynamics with a special source term as follows:

(1.1) { ∂ t ρ + ∂ x m = 0 , ∂ t m + ∂ x ( m 2 ρ + p ( ρ ) ) = F ( t , x ) ρ + H ( t , x ) m ,

where ρ , m = ρ u represent the density and momentum respectively, u is the velocity and p ( ρ ) = ρ γ γ is the pressure of the gas with γ ≥ 1 denoting the adiabatic exponent. The source terms F ( t , x ) and H ( t , x ) are C 1 functions on ℝ + 2 = ℝ + × ℝ , satisfying | F ( t , x ) | , | H ( t , x ) | ≤ M 0 for some positive constant M 0 for all ( t , x ) .

System (1.1) is equipped with the following initial data:

(1.2) ( ρ , m ) t = 0 = ( ρ 0 ( x ) , m 0 ( x ) ) ∈ L ∞ .

In [1], Chen, Huang and Wang first obtained that the weak entropy solution of isentropic Euler equations converges strongly to the corresponding isothermal Euler equations for homogeneous case. And there have been many results for the existence of entropy solutions of Euler equations with source terms. Cao, Huang, Li and Yu in [2] acquired the global existence of entropy solutions independent of time to isentropic compressible Euler equations with the source term same as in this context. Cao, Huang and Yuan [3] obtained the uniform bound (independent of time) of approximate solutions both for isentropic ( 1 < γ < 3 ) and isothermal Euler equations in a general nozzle. Chen and Luo in [4] obtained a convergence theorem of the fractional step Lax-Friedrichs scheme and Godunov scheme for a general source term ( U ( ρ , u , x , t ) , V ( ρ , u , x , t ) ) and 1 < γ ≤ 5 3 by using compensated compactness framework. When H ( t , x ) = 0 , Naoki Tsuge [5] proved the existence of a global solution without any boundary condition on the external force by employing an invariant region through a modified difference scheme. When F ( t , x ) = E ( t , x ) , where E ( t , x ) stands for electric field, Li, Cheng and Yu in [6] got the global existence and large time behavior of entropy solutions. Fang and Yu in [7] proved that the L ∞ weak solutions derived by Lax-Friedrichs scheme are uniformly bounded in time when F ( t , x ) = E ( t , x ) and H ( t , x ) = − 1 .

For isentropic Euler equations with source terms, existing results have primarily focused on establishing the existence of solutions, while investigations concerning the isothermal limit remain limited. This article is devoted to extending the findings from [1] to a class of Euler equations with special source terms. We show that the L ∞ entropy solutions of the Euler equations with such a source term converge strongly to the corresponding entropy solutions of the isothermal Euler equations as the adiabatic exponent γ → 1 . This is achieved by combining entropy analysis, a refined kinetic formulation and the compensated compactness argument to obtain the necessary uniform estimates for the limit. For the case considered here, the point is to obtain a uniform bound (independent of γ ) for the weak entropy solutions when γ > 1 . Inspired by the approach in [2], we set a new control function, using the maximum principle, and finally get the uniform bound independent of γ for the solution.

2. Preliminary Knowledge and Main Results

In this section, we provide some known results about hyperbolic conservation laws and present two main theorems of this article.

Firstly, system (1.1)-(1.2) can be written as a hyperbolic system of balance laws of the form:

(2.1) { ∂ t U + ∂ x F ( θ ) ( U ) = G ( U ) , U 0 ( x ) = ( ρ 0 ( x ) , m 0 ( x ) ) ,

where U = ( ρ , m ) ⊤ , F ( θ ) ( U ) = ( m , m 2 ρ + p ( θ ) ( ρ ) ) ⊤ , G ( U ) = ( 0 , F ( t , x ) ρ + H ( t , x ) m ) ⊤ and θ = γ − 1 2 , p ( θ ) ( ρ ) = ρ 2 θ + 1 2 θ + 1 . System (2.1) represents the isentropic gas when θ > 0 and the isothermal case when θ = 0 .

For system (2.1), a general entropy pair ( η , q ) satisfies the following hyperbolic system:

∇ q ( ρ , m ) = ∇ η ( ρ , m ) ∇ F ( θ ) ( ρ , m ) .

A weak entropy is an entropy η ( ρ , m ) that vanishes at ρ = 0 (vacuum), the corresponding pair ( η , q ) is then called a weak entropy pair.

Next, we introduce the concept of weak entropy solutions for system (1.1)-(1.2).

Definition2.1.Let θ ≥ 0 . Ifforanytestfunction φ ∈ C 0 ∞ ( ℝ + 2 ) , function U ( t , x ) = ( ρ , m ) ( t , x ) ∈ L ∞ ( ℝ + 2 ) satisfies

(2.2) ∫ ℝ + 2 ( U ∂ t φ + F ( θ ) ( U ) ∂ x φ + G ( U ) φ ) d x d t + ∫ ℝ U 0 ( x ) φ ( 0 , x ) d x = 0 ,

andifwefurtherhavefor φ ≥ 0 , foranyweakentropy-entropyflux η and q :

(2.3) ∫ ℝ + 2 ( η ( U ) ∂ t φ + q ( U ) ∂ x φ + ∇ η ( U ) G ( U ) φ ) d x d t +   ∫ ℝ η ( U 0 ) φ ( 0 , x ) d x ≥ 0 ,

thenthefunction U ( t , x ) = ( ρ , m ) ( t , x ) iscalledanentropysolutionoftheCauchyproblem (1.1)-(1.2).

We now discuss the isentropic and isothermal cases respectively.

2.1. Isentropic Case ( θ > 0 )

The Riemann invariants of system (2.1) are

w j ( θ ) = m ρ + ( − 1 ) j + 1 ρ θ θ     for   j = 1 , 2.

Weak entropy-entropy flux pairs can be expressed as follows:

Lemma2.1.For θ > 0 , theweakentropypairsof (2.1) canberepresentedasthefollowingform:

η ( θ ) ( ρ , m ; ψ ) = ∫ ℝ χ ( θ ) ( ρ ; s − m ρ ) ψ ( s ) d s , q ( θ ) ( ρ , m ; ψ ) = ∫ ℝ ( θ s + ( 1 − θ ) m ρ ) χ ( θ ) ( ρ ; s − m ρ ) ψ ( s ) d s ,

forany ψ ∈ C 2 ( ℝ ) , where

(2.4) χ ( θ ) ( ρ ; s − u ) = a θ [ ( ρ θ θ ) 2 − ( s − u ) 2 ] + 1 − θ 2 θ = a θ [ ( w 1 − s ) ( s − w 2 ) ] + 1 − θ 2 θ ,

istheweakentropykernelof (2.1), a θ = θ 1 θ ( ∫ − 1 1 [ 1 − τ 2 ] + 1 − θ 2 θ d τ ) − 1 and [ x ] + = max ( 0 , x ) (cf. [8]).

Using the standard change of variable s = u + ρ θ θ τ , we have

(2.5) η ( θ ) ( ρ , m ; ψ ) = a θ ∫ ℝ [ ( ρ θ θ ) 2 − ( u − s ) 2 ] + 1 − θ 2 θ ψ ( s ) d s = ρ ∫ ℝ ψ ( u + ρ θ θ τ ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ ∫ ℝ [ 1 − τ 2 ] + 1 − θ 2 θ d τ , (2.6) q ( θ ) ( ρ , m ; ψ ) = a θ ∫ ℝ ( θ s + ( 1 − θ ) u ) [ ( ρ θ θ ) 2 − ( u − s ) 2 ] + 1 − θ 2 θ ψ ( s ) d s = ρ ∫ ℝ ( u + ρ θ τ ) ψ ( u + ρ θ θ τ ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ ∫ ℝ [ 1 − τ 2 ] + 1 − θ 2 θ d τ .

The weak entropy pairs ( η ξ ( θ ) , q ξ ( θ ) ) of system (2.1) with θ > 0 under consideration are

(2.7) η ξ ( θ ) ( ρ , m ) = ρ ∫ ℝ e ξ 1 − ξ 2 ( m ρ + ρ θ θ τ ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ ∫ ℝ e ξ τ θ ( 1 − ξ 2 ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ ≥ 0 , (2.8) q ξ ( θ ) ( ρ , m ) = ρ ∫ ℝ ( m ρ + ρ θ τ ) e ξ 1 − ξ 2 ( m ρ + ρ θ θ τ ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ ∫ ℝ e ξ τ θ ( 1 − ξ 2 ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ .

These entropy pairs are obtained by choosing ψ ( s ) in (2.5)-(2.6) as

ψ ξ ( θ ) ( s ) = ∫ ℝ [ 1 − τ 2 ] + 1 − θ 2 θ d τ ∫ ℝ e ξ τ θ ( 1 − ξ 2 ) [ 1 − τ 2 ] + 1 − θ 2 θ d τ e ξ 1 − ξ 2 s .

This choice is the same as in [1].

Lemma2.2. (Existencetheoremfor θ > 0 ). Forany θ > 0 , thereexistsaboundedentropysolution ( ρ , m ) oftheCauchyproblem (1.1)-(1.2) on ℝ + 2 satisfying

(2.9) | m ( t , x ) | ρ ( t , x ) + ρ ( t , x ) θ − 1 θ ≤ C ( θ ) ,

where C ( θ ) isaconstantdependingon θ (cf. [2][3][8]-[12]).

2.2. Isothermal Case ( θ = 0 )

The Riemann invariants are

w j ( 0 ) = ρ e ( − 1 ) j + 1 m ρ     for   j = 1 , 2.

The family of weak entropy pairs can be shown as:

Lemma2.3.When θ = 0 , thetypicalentropypairs ( η ξ , q ξ ) are:

(2.10) η ξ ( ρ , m ) = ρ 1 1 − ξ 2 e ξ 1 − ξ 2 m ρ ,     q ξ ( ρ , m ) = ( m ρ + ξ ) η ξ     for   ξ ∈ ( − 1 , 1 ) ,

Lemma2.4. (CompactnessFrameworkfor θ = 0 ). Let ( ρ ε , m ε ) beasequenceofapproximatesolutionsof (2.1) with θ = 0 satisfying

(2.11) 0 ≤ ρ ε ( t , x ) ≤ C ,   | m ε ( t , x ) | ≤ ρ ε ( t , x ) ( | ln ρ ε ( t , x ) | + C )

where C > 0 isaconstantindependentof θ .Assumethereexistsasmallconstant δ > 0 suchthat, forany ξ ∈ ( − δ , δ ) ,

(2.12) ∂ t η ξ ( ρ ε , m ε ) + ∂ x q ξ ( ρ ε , m ε )

where ( η ξ , q ξ ) aretheweakentropypairsdefinedin (2.10).Thenthereexistsbothasequence (stilldenoted) ( ρ ε , m ε ) andavectorfunction ( ρ , m ) ( t , x ) suchthat

( ρ ε , m ε ) ( t , x ) → ( ρ , m ) ( t , x ) in L l o c p ( ℝ + 2 ) for all p ∈ [ 1 , ∞ ) as ε → 0 .

See Theorem 2.1 and Theorem 2.2 in [13] for details.

2.3. Main Theorem

Since we focus on the limit θ → 0 in the isentropic case, we take θ ∈ ( 0 , 1 2 ] throughout this paper. Two main results of this paper are as follows:

Theorem2.1. (uniformboundwithrespectto θ ). Let 0 < θ ≤ 1 2 . AssumethatthereexistsapositiveconstantCsuchthattheinitialdata ( ρ 0 ( θ ) ( t , x ) , m 0 ( θ ) ( t , x ) ) satisfies

(2.13) | m 0 ( θ ) ( t , x ) | ρ 0 ( θ ) ( t , x ) + ( ρ 0 ( θ ) ( t , x ) ) θ − 1 θ ≤ C ,

thenthereexistsaconstant C ′ dependingonlyontheinitialdata (independentof θ ), suchthattheweakentropysolutions ( ρ ( θ ) ( t , x ) , m ( θ ) ( t , x ) ) satisfies

(2.14) 0 ≤ ρ ( θ ) ( t , x ) ≤ e C ′ ,   | m ( θ ) ( t , x ) | ≤ ρ ( θ ) ( t , x ) ( | ln ρ ( θ ) ( t , x ) | + C ′ ) .

Theorem2.2 (Thestrongconvergence). Forany θ > 0 , ( ρ ( θ ) , m ( θ ) ) isaweakentropysolutionof (2.1) satisfyingcondition (2.9), andinitialdata ( ρ 0 ( θ ) , m 0 ( θ ) ) satisfies (2.13), thenthereexistsasequence (stilldenoted) ( ρ ( θ ) , m ( θ ) ) andavectorfunction ( ρ , m ) suchthat

( ρ ( θ ) , m ( θ ) ) ( t , x ) → ( ρ , m ) ( t , x ) in L l o c p ( ℝ + 2 ) for all p ∈ [ 1 , ∞ ) as θ → 0 .

Moreover, ( ρ , m ) isanentropysolutionof (2.1) when θ = 0 , satisfying

0 ≤ ρ ( t , x ) ≤ C ,   | m ( t , x ) | ≤ ρ ( t , x ) ( | ln ρ ( t , x ) | + C ) a.e. ( t , x ) ∈ ℝ + 2 ,

where C > 0 isaconstantindependentof θ .Furthermore, foranyentropypair ( η ξ , q ξ ) , ξ ∈ ( − 2 + 1 , 2 − 1 ) , U = ( ρ , m ) satisfiestheentropyinequality

∂ t η ξ ( U ) + ∂ x q ξ ( U ) + ∇ η ξ ( U ) G ( U ) ≤ 0 ,

inthesenseofdistributions.

3. The Uniform Estimate Independent of θ

In this section, we focus on proving Theorem 2.1.

Before the proof, we first recall some basic knowledge of system (2.1) when θ > 0 . The eigenvalues are

(3.1) λ 1 = m ρ − ρ θ ,     λ 2 = m ρ + ρ θ ,

and the corresponding right eigenvectors are as follows:

(3.2) r 1 = [ 1 λ 1 ] ,     r 2 = [ 1 λ 2 ] .

The Riemann invariants ( w , z ) are given by

(3.3) w = m ρ + ρ θ − 1 θ ,     z = m ρ − ρ θ − 1 θ .

The Riemann invariants chosen here are different from [2] because in this paper we focus on the solution with θ → 0 , one of the advantages of this choice is that it can ensure the boundary of ρ is independent of θ .

Next we introduce a maximum principle which will be used in the proof and the corresponding proof process is outlined in [2].

Lemma3.1 (Maximumprinciple) Let ( p , q ) ( x , t ) , ( x , t ) ∈ ℝ × [ 0 , T ] beanyboundedclassicalsolutionofthefollowingquasilinearparabolicsystem

{ p t + μ 1 p x = ε p x x + a 11 p + a 12 q + R 1 , q t + μ 2 q x = ε q x x + a 21 p + a 22 q + R 2 ,

withinitialdata p ( x , 0 ) ≤ 0 , q ( x , 0 ) ≥ 0 , wherethecoefficients μ i and a i j areboundedwithrespectto ( x , t ) andmaydependonpandq.Thesourceterms R i mayalsodependon p , q and x , t .Assumethat a 12 , a 21 ≤ 0 , R 1 ≤ 0 , R 2 ≥ 0 .Thenforany ( x , t ) , wehave p ( x , t ) ≤ 0 , q ( x , t ) ≥ 0 .

Remark. Lemma 3.1 holds true for ε = 0 . When ε > 0 , we can first construct an approximate solution and then prove this solution is bounded when ε → 0 . Detailed steps can be referred to [2].

ProofofTheorem2.1

By the formulas of Riemann invariants (3.3), we can decouple system (1.1) as

(3.4) { w t + λ 2 w x = F + H 2 ( w + z ) , z t + λ 1 z x = F + H 2 ( w + z ) .

Set a control function ϕ as follows:

ϕ = M + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t ,

where M ≥ 0 is a constant satisfying M ≥ max { | w ( x , 0 ) | , | z ( x , 0 ) | } . A direct calculation tells us that

ϕ x = 0 ,     ϕ t = | H | L ∞ e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ )                             = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] .

Define the new Riemann invariants ( w ¯ , z ¯ ) as

(3.5) w ¯ = w − ϕ ,     z ¯ = z + ϕ ,

so inserting (3.5) to (3.4), we can get a decoupled system of ( w ¯ , z ¯ )

(3.6) { w ¯ t + λ 2 w ¯ x = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 ( w ¯ + z ¯ ) , z ¯ t + λ 1 z ¯ x = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 ( w ¯ + z ¯ ) .

When z ¯ ≤ 2 ϕ , we have

(3.7) sup { − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 z ¯ } = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + | F | L ∞         + | H | L ∞ [ M + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t ] = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) + | F | L ∞ + M | H | L ∞ = ( | F | L ∞ + M | H | L ∞ ) ( 1 − e | H | L ∞ t ) ≤ 0 , w ¯ ≥ − 2 ϕ (3.8) inf { e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 w ¯ } = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] − | F | ∞         − | H | ∞ [ M + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t ] = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) − ( | F | L ∞ + M | H | L ∞ ) = ( | F | L ∞ + M | H | L ∞ ) ( e | H | L ∞ t − 1 ) ≥ 0.

Having obtained an upper bound for w and a lower bound for z , we next need to establish w is bounded below and z is bounded above. So we define another new invariants as follows:

(3.9) w ˜ = w + ϕ ,     z ˜ = z − ϕ .

Again the steps above, we can get another decoupled system of ( w ˜ , z ˜ )

(3.10) { w ˜ t + λ 2 w ˜ x = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 ( w ˜ + z ˜ ) , z ˜ t + λ 1 z ˜ x = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 ( w ˜ + z ˜ ) .

When z ¯ ≥ − 2 ϕ , we have

(3.11) inf { e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 w ˜ } = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] − | F | L ∞         − | H | L ∞ [ M + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t ] = e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) − ( | F | L ∞ + M | H | L ∞ ) = ( | F | L ∞ + M | H | L ∞ ) ( e | H | L ∞ t − 1 ) ≥ 0 , w ¯ ≤ 2 ϕ (3.12) sup { − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + F + H 2 z ˜ } = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) [ | H | L ∞ t + 1 ] + | F | L ∞         + | H | L ∞ [ M + e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) t ] = − e | H | L ∞ t ( | F | L ∞ + M | H | L ∞ ) + | F | L ∞ + M | H | L ∞ = ( | F | L ∞ + M | H | L ∞ ) ( 1 − e | H | L ∞ t ) ≤ 0.

To ensure the condition of initial data, we need

M ≥ max { | w ( x , 0 ) | , | z ( x , 0 ) | } .

Applying the maximum principle Lemma 3.1, we have the estimate

(3.13) | w ( x , t ) | ≤ C ,     | z ( x , t ) | ≤ C ,

where C is a constant independent of θ .

Using Lemma 3.1 of reference [1], we extend these uniform bounds from the Riemann invariants to the solution ( ρ ( θ ) , m ( θ ) ) , hence (2.14) holds.

Figure 1

4. Proof of Theorem 2.2

With the uniform bounds on the solution established, we now prove Theorem 2.2. We divide the proof into two steps.

Step1. H l o c − 1 compactnessoftheentropypair. For any θ > 0 , let ( ρ ( θ ) , m ( θ ) ) be the corresponding entropy solutions of (2.1) constructed in Lemma 2.2. It follows from (2.14) that the solution sequence ( ρ ( θ ) , m ( θ ) ) is uniformly bounded with respect to θ , which satisfies the condition (2.11) in Lemma 2.4.

We now show that the solution sequence ( ρ ( θ ) , m ( θ ) ) satisfies the condition (2.12) in Lemma 2.4 for any δ ∈ ( 0 , 2 − 1 ) . We will apply Murat lemma to achieve the goal.

Lemma4.1. (MuratLemma) Let Ω ∈ ℝ n beanopenset, then

( c o m p a c t   s e t   o f   W l o c − 1 , q ( Ω ) ) ∩ ( b o u n d e d   s e t   o f   W l o c − 1 , r ( Ω ) ) ⊂ ( c o m p a c t   s e t   o f   H l o c − 1 ( Ω ) ) ,

where 1 < q ≤ 2 < r .

For any weak entropy pairs given in Lemma 2.1, multiplying (1.1) by ∇ η ξ ( θ ) with η ξ ( θ ) defined in (2.7), and a calculation tells us that

η ξ m ( θ ) ( ρ , m ) = ξ 1 − ξ 2 1 ρ η ξ ( θ ) ,

so we have

η t + q x = η ξ m ( θ ) ( F ( x , t ) ρ + H ( x , t ) m ) ≤ C ( ρ M ′ + m M ′ ) ,

where C , M ′ > 0 are two constants. This implies that η t + q x is bounded in L l o c − 1 ( Ω ) . Therefore

η t + q x is compact in W l o c − 1 , α ( Ω ) with some 1 < α < 2 .

On the other hand, since ρ and m are uniformly bounded, we have

η t + q x is bounded in W l o c − 1 , ∞ ( Ω ) .

We conclude that

(4.1) η t + q x

for all weak entropy pairs with the help of the Murat Lemma.

By (4.1) and the compactness framework Lemma 2.4, we can prove that there exists a subsequence of ( ρ ( θ ) , m ( θ ) ) (still denoted by ( ρ ( θ ) , m ( θ ) ) ) such that

( ρ ( θ ) , m ( θ ) ) ( t , x ) → ( ρ , m ) ( t , x ) in L l o c p ( ℝ + 2 ) for all p ∈ [ 1 , ∞ ) as θ → 0 .

It’s easy to see that ( ρ , m ) is a weak solution to the Cauchy problem (1.1)-(1.2) so we omit here.

Step2.Entropyinequality. We next prove that ( ρ , m ) ( t , x ) satisfies the entropy inequality: for any ( η ξ , q ξ ) , there is

(4.2) ∂ t η ξ ( ρ , m ) + ∂ x q ξ ( ρ , m ) + ∇ η ξ ( ρ , m ) ⋅ G ( ρ , m ) ≤ 0 ,

in the sense of distributions.

As we know that ( ρ ( θ ) , m ( θ ) ) is the entropy solution of system (2.1) when θ > 0 , so for any ( ρ ξ ( θ ) , m ξ ( θ ) ) and ( η ξ ( θ ) , q ξ ( θ ) ) , there is

(4.3) ∂ t η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∂ x q ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∇ η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) ≤ 0 ,

in the sense of distributions. Next we prove that as θ → 0 , (4.3) converges to (4.2), that is

(4.4) | ∂ t η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∂ x q ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∇ η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) − [ ∂ t η ξ ( ρ , m ) + ∂ x q ξ ( ρ , m ) + ∇ η ξ ( ρ , m ) ⋅ G ( ρ , m ) ] | → 0 ,     as   θ → 0.

Before this part of proof, we state a lemma which has been proved in [1].

Lemma4.2 (RelationBetween ( η ( θ ) , q ( θ ) ) )

Givenany θ > 0 andanyfunction ( ρ ( θ ) , m ( θ ) ) satisfying (2.14), thenforanyinterval [ a , b ] ⋐ ( − 2 + 1 , 2 − 1 ) , thereexists C > 0 independentof θ suchthat

(4.5) | η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − η ξ ( ρ ( θ ) , m ( θ ) ) | + | q ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − q ξ ( ρ ( θ ) , m ( θ ) ) | ≤ C θ ,

forany ξ ∈ [ a , b ] .

First we have

(4.6) | ∂ t η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∂ x q ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) + ∇ η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) − [ ∂ t η ξ ( ρ , m ) + ∂ x q ξ ( ρ , m ) + ∇ η ξ ( ρ , m ) ⋅ G ( ρ , m ) ] | ≤ | ∂ t η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − ∂ t η ξ ( ρ , m ) | + | ∂ x q ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − ∂ x q ξ ( ρ , m ) |         + | ∇ η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) − ∇ η ξ ( ρ , m ) ⋅ G ( ρ , m ) | = G 1 + G 2 + G 3 ,

where

(4.7) G 1 ≤ | ∂ t η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − ∂ t η ξ ( ρ ( θ ) , m ( θ ) ) |               + | ∂ t η ξ ( ρ ( θ ) , m ( θ ) ) − ∂ t η ξ ( ρ , m ) | → 0 ,     as   θ → 0

The inequality holds since Lemma 4.2 and ( ρ ( θ ) , m ( θ ) ) → ( ρ , m ) (as θ → 0 ). Similarly, we have

(4.8) G 2 → 0 ,     as   θ → 0 G 3 G 3 ≤ | ∇ η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) − ∇ η ξ ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) |       + | ∇ η ξ ( ρ ( θ ) , m ( θ ) ) ⋅ G ( ρ ( θ ) , m ( θ ) ) − ∇ η ξ ( ρ , m ) ⋅ G ( ρ , m ) | = G 31 + G 32 .

A calculation tells us that

η ξ m ( ρ , m ) = ξ 1 − ξ 2 1 ρ η ξ .

So

G 31 = | η ξ m ( θ ) ( ρ ( θ ) , m ( θ ) ) ( F ( x , t ) ρ ( θ ) + H ( x , t ) m ( θ ) )       − η ξ m ( ρ ( θ ) , m ( θ ) ) ( F ( x , t ) ρ ( θ ) + H ( x , t ) m ( θ ) ) | = | F ( x , t ) ρ ( θ ) + H ( x , t ) m ( θ ) | | η ξ m ( θ ) ( ρ ( θ ) , m ( θ ) ) − η ξ m ( ρ ( θ ) , m ( θ ) ) | = | F ( x , t ) ρ ( θ ) + H ( x , t ) m ( θ ) | ξ 1 − ξ 2 1 ρ ( θ ) | η ξ ( θ ) ( ρ ( θ ) , m ( θ ) ) − η ξ ( ρ ( θ ) , m ( θ ) ) | ≤ | M 0 + M 0 C | ξ 1 − ξ 2 C θ → 0     as   θ → 0 ,

where C and M 0 are constants independent of θ .

As for G 32 , a similar estimate together with the strong convergence ( ρ ( θ ) , m ( θ ) ) → ( ρ , m ) , as θ → 0 yields G 32 → 0 , as θ → 0 .

Combining the estimates for G 1 , G 2 , G 3 , we conclude that the right-hand side of (4.6) converges to 0 as θ → 0 . Hence the entropy inequality holds for the limit function ( ρ , m ) .

This completes the proof of Theorem 2.2.

Figure 2

5. Conclusion

In this paper, we extend the result of [1] to a class of isentropic Euler equations with a special source term. To apply the proving framework introduced in [1], we need to show that solutions of isentropic Euler equations are independent of the adiabatic exponent γ . Therefore, we employ a maximum principle to achieve this. Next, we can easily prove the strong convergence and get the entropy inequality.

References

1.. Chen, G.G., Huang, F. and Wang, T. (2024) Isothermal Limit of Entropy Solutions of the Euler Equations for Isentropic Gas Dynamics. SIAM Journal on Mathematical Analysis, 56, 1300-1320. https://doi.org/10.1137/23m1549948 10.1137/23m1549948https://doi.org/10.1137/23m1549948 Chen, G.G. Huang, F. Wang, T. 2024 Isothermal Limit of Entropy Solutions of the Euler Equations for Isentropic Gas Dynamics SIAM Journal on Mathematical Analysis 56 10.1137/23m1549948

2.. Cao, W.T., Huang, F.M., Li, T.H. and Yu, H.M. (2016) Global Entropy Solutions to an Inhomogeneous Isentropic Compressible Euler System. Acta Mathematica Scientia, 36, 1215-1224. https://doi.org/10.1016/s0252-9602(16)30063-7 10.1016/s0252-9602(16)30063-7https://doi.org/10.1016/s0252-9602(16)30063-7 Cao, W.T. Huang, F.M. Li, T.H. Yu, H.M. 2016 Global Entropy Solutions to an Inhomogeneous Isentropic Compressible Euler System Acta Mathematica Scientia 9602 16 10.1016/s0252-9602(16)30063-7

3.. Cao, W.T., Huang, F.M. and Yuan, D.F. (2019) Global Entropy Solutions to the Gas Flow in General Nozzle. SIAM Journal on Mathematical Analysis, 51, 3276-3297. https://doi.org/10.1137/19m1249436 10.1137/19m1249436https://doi.org/10.1137/19m1249436 Cao, W.T. Huang, F.M. Yuan, D.F. 2019 Global Entropy Solutions to the Gas Flow in General Nozzle SIAM Journal on Mathematical Analysis 51 10.1137/19m1249436

4.. Ding, X.X., Chen, G.Q. and Luo, P.Z. (1989) Convergence of the Fractional Step Lax-Friedrichs Scheme and Godunov Scheme for the Isentropic System of Gas Dynamics. CommunicationsinMathematicalPhysics, 121, 63-84. https://doi.org/10.1007/bf01218624 10.1007/bf01218624https://doi.org/10.1007/bf01218624 Ding, X.X. Chen, G.Q. Luo, P.Z. 1989 Convergence of the Fractional Step Lax-Friedrichs Scheme and Godunov Scheme for the Isentropic System of Gas Dynamics Communications in Mathematical Physics 121 10.1007/bf01218624

5.. Tsuge, N. (2016) Existence and Stability of Solutions to the Compressible Euler Equations with an Outer Force. Nonlinear Analysis: Real World Applications, 27, 203-220. https://doi.org/10.1016/j.nonrwa.2015.07.017 10.1016/j.nonrwa.2015.07.017https://doi.org/10.1016/j.nonrwa.2015.07.017 Tsuge, N. 2016 Existence and Stability of Solutions to the Compressible Euler Equations with an Outer Force Nonlinear Analysis: Real World Applications 27 10.1016/j.nonrwa.2015.07.017

6.. Li, Y., Cheng, Y. and Yu, H. (2020) Existence and Large Time Behavior of Entropy Solutions to One-Dimensional Unipolar Hydrodynamic Model for Semiconductor Devices with Variable Coefficient Damping. Advances in Mathematical Physics, 2020, 1-7. https://doi.org/10.1155/2020/8881435 10.1155/2020/8881435https://doi.org/10.1155/2020/8881435 Li, Y. Cheng, Y. Yu, H. 2020 Existence and Large Time Behavior of Entropy Solutions to One-Dimensional Unipolar Hydrodynamic Model for Semiconductor Devices with Variable Coefficient Damping Advances in Mathematical Physics 2020 10.1155/2020/8881435

7.. Fang, X.X. and Yu, H.M. (2018) Uniform Boundedness in Weak Solutions to a Specific Dissipative System. Journal of Mathematical Analysis and Applications, 461, 1153-1164. https://doi.org/10.1016/j.jmaa.2018.01.047 10.1016/j.jmaa.2018.01.047https://doi.org/10.1016/j.jmaa.2018.01.047 Fang, X.X. Yu, H.M. 2018 Uniform Boundedness in Weak Solutions to a Specific Dissipative System Journal of Mathematical Analysis and Applications 461 10.1016/j.jmaa.2018.01.047

8.. Lions, P.L., Perthame, B. and Tadmor, E. (1994) Kinetic Formulation of the Isentropic Gas Dynamics and p-Systems. Communications in Mathematical Physics, 163, 415-431. https://doi.org/10.1007/bf02102014 10.1007/bf02102014https://doi.org/10.1007/bf02102014 Lions, P.L. Perthame, B. Tadmor, E. 1994 Kinetic Formulation of the Isentropic Gas Dynamics and p-Systems Communications in Mathematical Physics 163 10.1007/bf02102014

9.. DiPerna, R.J. (1983) Convergence of the Viscosity Method for Isentropic Gas Dynamics. Communications in Mathematical Physics, 91, 1-30. https://doi.org/10.1007/bf01206047 10.1007/bf01206047https://doi.org/10.1007/bf01206047 DiPerna, R.J. 1983 Convergence of the Viscosity Method for Isentropic Gas Dynamics Communications in Mathematical Physics 91 10.1007/bf01206047

10.. Ding, X.X., Chen, G.Q. and Luo, P.Z. (1985) Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (I). ActaMathematicaScientia, 5, 415-432. https://doi.org/10.1016/s0252-9602(18)30542-3 10.1016/s0252-9602(18)30542-3https://doi.org/10.1016/s0252-9602(18)30542-3 Ding, X.X. Chen, G.Q. Luo, P.Z. 1985 Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (I) Acta Mathematica Scientia 9602 18 10.1016/s0252-9602(18)30542-3

11.. Ding, X.X., Chen, G.Q. and Luo, P.Z. (1985) Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (II). ActaMathematicaScientia, 5, 433-472. https://doi.org/10.1016/s0252-9602(18)30543-5 10.1016/s0252-9602(18)30543-5https://doi.org/10.1016/s0252-9602(18)30543-5 Ding, X.X. Chen, G.Q. Luo, P.Z. 1985 Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (II) Acta Mathematica Scientia 9602 18 10.1016/s0252-9602(18)30543-5

12.. Chen, G.Q. (1986) Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (III). ActaMathematicaScientia, 6, 75-120. https://doi.org/10.1016/s0252-9602(18)30535-6 10.1016/s0252-9602(18)30535-6https://doi.org/10.1016/s0252-9602(18)30535-6 Chen, G.Q. 1986 Convergence of the Lax-Friedrichs Scheme for Isentropic Gas Dynamics (III) Acta Mathematica Scientia 9602 18 10.1016/s0252-9602(18)30535-6

13.. Huang, F.M. and Wang, Z. (2003) Convergence of Viscosity Solutions for Isothermal Gas Dynamics. SIAM Journal on Mathematical Analysis, 34, 595-610. https://doi.org/10.1137/s0036141002405819 10.1137/s0036141002405819https://doi.org/10.1137/s0036141002405819 Huang, F.M. Wang, Z. 2003 Convergence of Viscosity Solutions for Isothermal Gas Dynamics SIAM Journal on Mathematical Analysis 34 10.1137/s0036141002405819