Research

# Convergence rate of solutions toward stationary solutions to the bipolar Navier-Stokes-Poisson equations in a half line

Fang Zhou1 and Yeping Li2*

Author Affiliations

1 Department of Mathematics, Hubei University of Science and Technology, Xianning, 437100, P.R. China

2 Department of Mathematics, Shanghai Normal University, Shanghai, 200234, P.R. China

For all author emails, please log on.

Boundary Value Problems 2013, 2013:124  doi:10.1186/1687-2770-2013-124

 Received: 2 November 2012 Accepted: 26 April 2013 Published: 14 May 2013

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

In this paper, we show the convergence rate of a solution toward the stationary solution to the initial boundary value problem for the one-dimensional bipolar compressible Navier-Stokes-Poisson equations. For the supersonic flow at spatial infinity, if an initial perturbation decays with the algebraic or the exponential rate in the spatial asymptotic point, the solution converges to the corresponding stationary solution with the same rate in time as time tends to infinity. For the transonic flow at spatial infinity, the solution converges to the stationary solution in time with the lower rate than that of the initial perturbation in the spatial. These results are proved by the weighted energy method.

MSC: 35M31, 35Q35.

##### Keywords:
convergence rate; Navier-Stokes-Poisson equation; stationary wave; weighted energy method

### 1 Introduction

In this paper, we are concerned with the following bipolar Navier-Stokes-Poisson equations:

{ t ρ 1 + x ( ρ 1 u 1 ) = 0 , t ( ρ 1 u 1 ) + x ( ρ 1 u 1 2 + P 1 ( ρ 1 ) ) = μ 1 u 1 x x + ρ 1 E , t ρ 2 + x ( ρ 2 u 2 ) = 0 , t ( ρ 2 u 2 ) + x ( ρ 2 u 2 2 + P 2 ( ρ 2 ) ) = μ 2 u 2 x x ρ 2 E , E x = ρ 1 ρ 2 , (1.1)

in a one-dimensional half space R + : = ( 0 , ) . Here the unknown functions are the densities ρ i ( i = 1 , 2 ) > 0 , the velocities u i ( i = 1 , 2 ), and the electron field E. P i ( ρ i ) ( i = 1 , 2 ) is the pressure depending only on the density. μ i ( i = 1 , 2 ) is viscosity coefficient. Throughout this paper, we assume that two fluids of electrons and ions have the same equation of state P 1 ( ) = P 2 ( ) = P ( ) with P ( ρ ) = K ρ γ for K > 0 and γ 1 , and also they have the same viscosity coefficients μ 1 = μ 2 = 1 . The bipolar Navier-Stokes-Poisson system is used to simulate the transport of charged particles (e.g., electrons and ions). It consists of the compressible Navier-Stokes equation of two-fluid under the influence of the electro-static potential force governed by the self-consisted Poisson equation. Note that when we only consider one particle in the fluids, we also have the unipolar Navier-Stokes-Poisson equations. For more details, we can refer to [1-4].

Recently, some important progress was made for the compressible unipolar Navier-Stokes-Poisson system. The local and/or global existence of a renormalized weak solution for the Cauchy problem of the multi-dimensional compressible Navier-Stokes-Poisson system were proved in [5-7]. Chan [8] also considered the nonexistence of global weak solutions to the Navier-Stokes-Poisson equations in R N . Hao and Li [9] established the global strong solutions of the initial value problem for the multi-dimensional compressible Navier-Stokes-Poisson system in a Besov space. The global existence and L 2 -decay rate of the smooth solution of the initial value problem for the compressible Navier-Stokes-Poisson system in R 3 were achieved by Li and his collaborators in [10,11]. The pointwise estimates of the smooth solutions for the three-dimensional isentropic compressible Navier-Stokes-Poisson equation were obtained in [12]. The quasineutral limit of the compressible Navier-Stokes-Poisson system was studied in [13-15]. However, the results about the bipolar Navier-Stokes-Poisson equations are very few. Lastly, Li et al.[16] showed the global existence and asymptotic behavior of smooth solutions for the initial value problem of the bipolar Navier-Stokes-Poisson equations. Duan and Yang [17] studied the unique existence and asymptotic stability of a stationary solution for the initial boundary value problem, and they showed that the large-time behavior of solutions for the bipolar Navier-Stokes-Poisson equations coincided with the one for the single Navier-Stokes system in the absence of the electric field. The consistency is also observed and proved between the bipolar Euler-Poisson system and the single damped Euler equation; for example, see [18-20] and the references therein.

In this paper, we are going to discuss the initial-boundary value problem for the one-dimensional bipolar Navier-Stokes-Poisson equations. Now we give the initial condition

( ρ 1 , u 1 , ρ 2 , u 2 ) ( x , 0 ) = ( ρ 10 , u 10 , ρ 20 , u 20 ) ( x ) ( ρ + , u + , ρ + , u + ) as  x , (1.2)

and the boundary date

u 1 ( 0 , t ) = u 2 ( 0 , t ) = u b < 0 . (1.3)

Here, we suppose inf x R + ρ i 0 ( i = 1 , 2 ) > 0 and further the compatibility condition u b = u 10 ( 0 ) = u 20 ( 0 ) . Moreover, for the unique existence, we also assume

E ( + , t ) = 0 . (1.4)

In [17], the authors showed that the solution to (1.1)-(1.4) converges to the corresponding stationary solution of the single Navier-Stokes system in the absence of the electric field

{ t ρ + x ( ρ u ) = 0 , t ( ρ u ) + x ( ρ u + P ( ρ ) ) = u x x , (1.5)

as time tends to infinity. Then, let ( ρ ˜ , u ˜ ) ( x ) be the stationary solution to the system (1.5). We know that the stationary solution ( ρ ˜ , u ˜ ) satisfies

{ ( ρ ˜ u ˜ ) x = 0 , ( ρ ˜ u ˜ 2 + P ( ρ ˜ ) ) x = μ ρ ˜ x x , (1.6)

and the boundary and spatial asymptotic conditions

u ˜ ( 0 ) = u b , lim x ( ρ ˜ , u ˜ ) = ( ρ + , u + ) , inf x R + ρ ˜ ( x ) > 0 . (1.7)

In this paper, we are mainly concerned with the decay rate of solutions to (1.1)-(1.4) toward the stationary solution ( ρ ˜ , u ˜ , ρ ˜ , u ˜ , 0 ) . Now we state the main result in the following theorem.

Theorem 1.1Suppose that M + 1 and u b < u hold. The initial data ( ρ 10 , u 10 , ρ 20 , u 20 , E 0 ) ( x ) is supposed to satisfy

( ρ i 0 , u i 0 ) ( i = 1 , 2 ) H 1 ( R + ) , E 0 ( x ) L 2 ( R + ) , inf x R + ( ρ 10 , ρ 20 ) > 0 , (1.8)

and there exists a positive constant ε 0 such that

( ρ 10 ρ ˜ , u 10 u ˜ , ρ 20 ρ ˜ , u 20 u ˜ ) 1 + E 0 + δ < ε 0 . (1.9)

(i) When M + > 1 , in addition, the initial data also satisfies ( 1 + x ) α 2 ( ρ 10 ρ ˜ ) , ( 1 + x ) α 2 ( u 10 u ˜ ) , ( 1 + x ) α 2 ( ρ 20 ρ ˜ ) , ( 1 + x ) α 2 ( u 20 u ˜ ) , ( 1 + x ) α 2 E 0 L 2 ( R + ) for a certain positive constantα, then the solution ( ρ 1 , u 1 , ρ 2 , u 2 , E ) to (1.1)-(1.3) satisfies the decay estimate

( ρ 1 ρ ˜ , u 1 u ˜ , ρ 2 ρ ˜ , u 2 u ˜ , E ) L C ( 1 + x ) α 2 . (1.10)

On the other hand, if the initial data satisfies e ζ 2 x ( ρ 10 ρ ˜ ) , e ζ 2 x ( u 10 u ˜ ) , e ζ 2 x ( ρ 20 ρ ˜ ) , e ζ 2 x ( u 20 u ˜ ) , e ζ 2 x E 0 L 2 ( R + ) for a certain positive constantζ, then there exists a positive constantαsuch that the solution ( ρ 1 , u 1 , ρ 2 , u 2 , E ) to (1.1)-(1.3) satisfies

( ρ 1 ρ ˜ , u 1 u ˜ , ρ 2 ρ ˜ , u 2 u ˜ , E ) L C e α t . (1.11)

(ii) When M + = 1 , and there exists a positive constant ε 0 such that if the initial data also satisfies ( 1 + x ) α 2 ( ρ 10 ρ ˜ , u 10 u ˜ , ρ 20 ρ ˜ , u 20 u ˜ ) 1 + ( 1 + x ) α 2 E 0 < ε 0 for a certain constantαsatisfying α [ 2 , α ) , where α is a constant defined by

α ( α 2 ) = 4 γ + 1 and α > 0 ,

then the solution ( ρ 1 , u 1 , ρ 2 , u 2 , E ) to (1.1)-(1.3) satisfies

( ρ 1 ρ ˜ , u 1 u ˜ , ρ 2 ρ ˜ , u 2 u ˜ , E ) L C ( 1 + t ) α 4 , (1.12)

where M + , u andδare defined in Section 2, and E 0 ( x ) = x ( ρ 10 ρ 20 ) ( y ) d y .

Notations Throughout this paper, C > 0 denotes the generic positive constant independent of time. L p ( R ) ( 1 p < ) denotes the space of measurable functions with the finite norm L p = ( R | | p d x ) 1 p , and L is the space of bounded measurable functions on ℝ with the norm L = ess sup x | | . We use to denote the L 2 -norm. H k ( R ) ( k 0 ) stands for the space of L 2 ( R ) -functions f whose derivatives (in the sense of distribution) D x l f ( l k ) are also L 2 ( R ) -functions with the norm k = ( l = 0 k D l 2 ) 1 2 . Moreover, C k ( [ 0 , T ] ; H l ( R ) ) ( k , l 0 ) denotes the space of the k-times continuously differentiable functions on the interval [ 0 , T ] with values in H l ( R ) .

The rest of the paper is organized as follows. In Section 2, we review the results of the stationary solution and the non-stationary solutions, then we reformulate our problem. Finally, we give the a priori estimates for the cases M + > 1 and M + = 1 in Section 3 and 4, respectively.

### 2 Stationary solution and global existence of non-stationary solution

In this section we mainly review the property of a stationary solution, and the unique existence and asymptotic behavior of non-stationary solutions for (1.1)-(1.3). To begin with, we recall the stationary equation

{ ( ρ ˜ u ˜ ) x = 0 , ( ρ ˜ u ˜ 2 + P ( ρ ˜ ) ) x = μ ρ ˜ x x (2.1)

with

u ˜ ( 0 ) = u b < 0 , lim x ( ρ ˜ , u ˜ ) = ( ρ + , u + ) , inf x R + ρ ˜ ( x ) > 0 . (2.2)

Integrating (2.1)1 over ( x , ) yields ρ ˜ ( x ) = ρ + u + ( u ˜ ( x ) ) 1 , which implies by letting x 0 + , ρ b = ρ ˜ ( 0 ) = ρ + u + ( u b ) 1 . Namely, u ˜ = u + v + v ˜ ( v ˜ = 1 ρ ˜ , v + = 1 ρ + ), which together with (2.1) implies

u + = v + v ( 0 ) u b < 0 . (2.3)

Thus, the condition u + < 0 has to be assumed whenever the outflow problem, i.e., the case u b < 0 , is consider. Moreover, let the strength of the boundary layer ( ρ ˜ , u ˜ ) ( x ) be measured by δ : = | u + u b | . Finally, we also define ( v , u ) as follows:

u + = u + v + v ( v + v ) [ P ( 1 v ) P ( 1 v + ) ] , u = u + v + v , (2.4)

and denote the Mach number at infinity M + = : | u + | P ( ρ + ) . Then one has the following lemma.

Lemma 2.1 (see [21,22])

Assume that the condition (2.3) holds. The boundary problem (2.1)-(2.2) has a smooth solution ( ρ ˜ , u ˜ ) ( x ) , if and only if M + 1 and u b < u . Moreover, if M + > 1 , there exist two positive constantsλandCsuch that the stationary solution ( ρ ˜ , u ˜ ) satisfies the estimate

| x k ( ρ ˜ ( x ) ρ + , u ˜ ( x ) u + ) | C δ e λ x for   k = 0 , 1 , 2 , . (2.5)

If M + = 1 , the stationary solution ( ρ ˜ , u ˜ ) satisfies

| x k ( ρ ˜ ( x ) ρ + , u ˜ u + ) | C δ k + 1 ( 1 + δ x ) k + 1 for   k = 0 , 1 , 2 , . (2.6)

As to the stability of the stationary solution of (1.1)-(1.4), Duan and Yang showed the following results in [17].

Lemma 2.2 (see [17])

Suppose that M + 1 and u b < u hold. In addition, the initial data ( ρ 10 , u 10 , ρ 20 , u 20 , E 0 ) is supposed to satisfy

( ρ 10 ρ ˜ , u 10 u ˜ , ρ 20 ρ ˜ , u 20 u ˜ ) H 1 ( R + ) , E 0 L 2 ( R + ) , inf x R + ( ρ 10 , ρ 20 ) > 0 .

Then there exists a positive constant ε 0 such that if

( ρ 10 ρ ˜ , u 10 u ˜ , ρ 20 ρ ˜ , u 20 u ˜ ) 1 + E 0 + δ < ε 0 ,

the initial boundary value problem (1.1)-(1.3) has a unique solution ( ρ 1 , u 1 , ρ 2 , u 2 , E ) X ( 0 , T ) for arbitrary T > 0 . Moreover, the solution ( ρ 1 , u 1 , ρ 2 , u 2 , E ) converges to the stationary solution ( ρ ˜ , u ˜ , ρ ˜ , u ˜ , 0 ) as time tends to infinity:

lim t sup x R + | ρ 1 ρ ˜ , u 1 u ˜ , ρ 2 ρ ˜ , u 2 u ˜ , E | = 0 .

Here the solution space X ( 0 , T ) is defined by

X ( 0 , T ) = { ( ρ 1 , u 1 , ρ 2 , u 2 , E ) : ρ 1 ρ ˜ , u 1 u ˜ , ρ 2 ρ ˜ , u 2 u ˜ C ( 0 , T ; H 1 ) , ( ρ 1 ρ ˜ ) x , ( ρ 2 ρ ˜ ) x L 2 ( 0 , T ; L 2 ) , ( u 1 u ˜ ) x , ( u 2 u ˜ ) x L 2 ( 0 , T ; H 1 ) , E C ( 0 , T ; L 2 ) , ( u 1 u ˜ ) ( t , 0 ) = ( u 2 u ˜ ) ( t , 0 ) = 0 ( 0 t T ) } .

Finally, to enclose this section, we reformulate the original problem in terms of the perturbed variables. Set ( φ 1 , ψ 1 , φ 2 , ψ 2 ) from the stationary solution as

φ i = ρ i ρ ˜ , ψ i = u i u ˜ , i = 1 , 2 .

Due to (1.1) and (1.6), we have the system of equations for ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) as

{ t φ 1 + u 1 x φ 1 + ρ 1 x ψ 1 = f 1 , ρ 1 t ψ 1 + P ( ρ 1 ) x φ 1 + ρ 1 u 1 x ψ 1 = ψ 1 x x g 1 + ρ 1 E , t φ 2 + u 2 x φ 2 + ρ 2 x ψ 2 = f 2 , ρ 2 t ψ 2 + P ( ρ 2 ) x φ 2 + ρ 2 u 2 x ψ 2 = ψ 2 x x g 2 ρ 2 E , E = x ( φ 1 φ 2 ) ( y , t ) d y , (2.7)

where the nonlinear terms f i ( i = 1 , 2 ) and g i ( i = 1 , 2 ) are given by

f i = u ˜ x φ i + ρ ˜ x ψ i , g i = u ˜ x ( u ˜ φ i + ρ i ψ i ) + ρ ˜ x ( P ( ρ i ) P ( ρ ˜ ) ) .

The initial and boundary condition to (2.7) are derived from (1.2), (1.3) and (1.4) as follows:

(2.8)

(2.9)

The uniform bound of the solutions in the weighted Sobolev space is derived later in Sections 3 and 4. For this purpose, we introduce the function spaces X ω ( 0 , T ) and X ω 1 ( 0 , T ) defined by

X ω ( 0 , T ) = { ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) : ( ω φ 1 , ω ψ 1 , ω φ 2 , ω ψ 2 , ω E ) C ( 0 , T ; L 2 ( R + ) ) }

and

X ω 1 ( 0 , T ) = { ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) : ( ω φ 1 , ω ψ 1 , ω φ 2 , ω ψ 2 ) C ( 0 , T ; H 1 ( R + ) ) ; ω E C ( 0 , T ; L 2 ( R + ) ) } .

Here the two types of weight functions are considered: ω ( x ) : = ( 1 + x ) α , or ω ( x ) = e α x . Also, we use the norms | | 2 , ω , | | a , α , and | | e , α defined by

| f | 2 , ω : = ( 0 ω ( x ) f ( x ) 2 d x ) 1 2 , | f | a , α : = | f | 2 , ( 1 + x ) α , | f | e , α : = | f | 2 , e α x .

The following lemma, concerning the existence of the solution locally in time, is proved by the standard iteration method. Hence we omit the proof.

Lemma 2.3If the initial data satisfies (1.8) and ω φ 10 , ω ψ 10 , ω φ 20 , ω ψ 20 , ω E 0 L 2 ( R + ) , there exists a positive constantTsuch that the initial boundary value problem (2.7)-(2.9) has a unique solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) X ω ( 0 , T ) . Moreover, if the initial data satisfies (1.8), (1.9) and ω φ 10 , ω ψ 10 , ω φ 20 , ω ψ 20 H 1 ( R + ) and ω E 0 L 2 ( R + ) , there exists a unique solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) in X ω 1 ( 0 , T ) .

### 3 A priori estimates for M + > 1

In this section, we derive the a priori estimates of the solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) for the case that M + > 1 holds in some Sobolev space. To summarize the a priori estimate, we use the following notation (see [23]) for a weight function W ( x , t ) = χ ( t ) ω ( x ) :

N ( t ) = sup 0 τ t ( φ 1 , ψ 1 , φ 2 , ψ 2 ) ( τ ) 1 , M ( t ) 2 = 0 t χ ( τ ) ( ( φ 1 x , φ 2 x , E x ) ( τ ) 2 + ( ψ 1 x , ψ 2 x ) ( τ ) 1 2 ) d τ + 0 t χ ( τ ) ( φ 1 ( τ , 0 ) 2 + φ 2 ( τ , 0 ) 2 + E ( τ , 0 ) 2 ) d τ , L ( t ) 2 = 0 t χ t ( τ ) ( | ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) ( τ ) | 2 , ω 2 + ( φ 1 x , ψ 1 x , φ 2 x , ψ 2 x ) ( τ ) 2 ) + χ ( τ ) ( | ( ψ 1 , ψ 2 ) ( τ ) | 2 , ω x x 2 + | ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) ( τ ) | 2 , | u ˜ x | ω 2 ) d τ .

Proposition 3.1Suppose that the same assumptions as in Theorem 1.1 hold.

(i) (Algebraic decay) Suppose that ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) X ( 1 + x ) α ( 0 , T ) is a solution to (2.7)-(2.9) for certain positive constantsαandT. Then there exist positive constants ε 0 andCsuch that if N ( T ) + δ ε 0 , then the solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) satisfies the estimate

( 1 + t ) α + ε ( ( φ 1 , ψ 1 , φ 2 , ψ 2 ) ( t ) 1 2 + E ( t ) 2 ) + 0 t ( 1 + τ ) α + ε ( ( φ 1 x , φ 2 x , E x ) ( τ ) 2 + ( ψ 1 x , ψ 2 x ) ( τ ) 1 2 ) d τ + 0 t ( 1 + τ ) α + ε | ( φ 1 , φ 2 , E , φ 1 x , φ 2 x ) ( τ , 0 ) | 2 d τ C ( ( φ 10 , ψ 10 , φ 20 , ψ 20 ) 1 2 + ( φ 10 , ψ 10 , φ 20 , ψ 20 , E 0 ) a , α 2 ) ( 1 + t ) ε (3.1)

for arbitrary t [ 0 , T ] and ε > 0 .

(ii) (Exponential decay) Suppose that ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) X e ζ x ( 0 , T ) is a solution to (2.7)-(2.9) for certain positive constantsζandT. Then there exist positive constants ε 0 , C, β (<ζ) andαsatisfying α β such that if N ( T ) + δ ε 0 , then the solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) satisfies the estimate

(3.2)

For the sake of clarity, we divide the proof of Proposition 4.1 into the following lemmas. We first derive the basic energy estimate.

Lemma 3.2Suppose that the same assumptions as in Theorem 1.1 hold. Then there exists a positive constant ε 0 such that N ( T ) + δ < ε 0 , it holds that

(3.3)

Proof From (2.7), a direct computation yields

( ρ 1 E 1 + ρ 2 E 2 + 1 2 E 2 ) t G 1 x + ψ 1 x 2 + ψ 2 x 2 = ( ψ 1 ψ 1 x + ψ 2 ψ 2 x ) x + R 1 , (3.4)

here

E i : = E ( ρ i , u i ) ( i = 1 , 2 ) , E ( ρ , u ) : = Φ ( ρ , ρ ˜ ) + 1 2 | u u ˜ | 2 , Φ ( ρ , ρ ˜ ) : = ρ ˜ ρ P ( s ) P ( ρ ˜ ) s 2 d s , G 1 = ρ 1 u 1 E 1 ρ 2 u 2 E 2 ( P ( ρ 1 ) P ( ρ ˜ ) ) ψ 1 ( P ( ρ 2 ) P ( ρ ˜ ) ) ψ 2 u ˜ 2 E 2 , R 1 = u ˜ x [ P ( ρ 1 ) P ( ρ ˜ ) P ( ρ ˜ ) φ 1 + ( ρ 1 u 1 ρ ˜ u ˜ ) ψ 1 + P ( ρ 2 ) P ( ρ ˜ ) P ( ρ ˜ ) φ 2 + ( ρ 2 u 2 ρ ˜ u ˜ ) ψ 2 1 2 E 2 ] 1 ρ ˜ P ( ρ ˜ ) x φ 1 ψ 1 1 ρ ˜ P ( ρ ˜ ) x φ 2 ψ 2 .

Owing to Lemmas 2.1 and 2.2, we see that the energy form E ( ρ , u ) is equivalent to | ( ρ ρ ˜ , u u ˜ ) | 2 . That is, there exist positive constants c and C such that

c ( φ i 2 + ψ i 2 ) E i C ( φ i 2 + ψ i 2 ) , i = 1 , 2 . (3.5)

We also have positive bounds of ρ i ( i = 1 , 2 ) as

0 < c ρ i ( i = 1 , 2 ) C , ( t , x ) [ 0 , T ] × R + . (3.6)

Further, multiplying (3.4) by a weight function W ( t , x ) = χ ( t ) ω ( x ) , we have

( W ρ 1 E 1 + W ρ 2 E 2 + 1 2 W E 2 ) t ( W G 1 ) x + W x G 1 + W ψ 1 x 2 + W ψ 2 x 2 = W t ( ρ 1 E 1 + ρ 2 E 2 + 1 2 E 2 ) + ( i = 1 2 W ψ i ψ i x 1 2 W x ψ i 2 ) x + 1 2 W x x ( ψ 1 2 + ψ 2 2 ) + W R 1 . (3.7)

Due to the boundary conditions (1.3) and (2.9), the integration of the second term on the left-hand side of (3.7) over R +

R + ( W [ ρ 1 u 1 E 1 + ρ 2 u 2 E 2 + ( P ( ρ 1 ) P ( ρ ˜ ) ) ψ 1 + ( P ( ρ 1 ) P ( ρ ˜ ) ) ψ 2 + u ˜ 2 E 2 ] ) x d x = χ ( t ) ρ 1 ( t , 0 ) u b E 1 ( t , 0 ) χ ( t ) ρ 2 ( t , 0 ) u b E 2 ( t , 0 ) 1 2 u b χ ( t ) E ( t , 0 ) 2 C χ ( t ) ( φ 1 ( t , 0 ) 2 + φ 2 ( t , 0 ) 2 + E ( t , 0 ) 2 ) , (3.8)

where we have used the estimates (3.5) and (3.6). Next, G 1 can be computed as

G 1 = G 11 + G 12 (3.9)

with

G 11 = K γ ρ + γ 2 | u + | 2 ( φ 1 2 + φ 2 2 ) + ρ + | u + | 2 ( ψ 1 2 + ψ 2 2 ) K γ ρ + γ 1 ( φ 1 ψ 1 + φ 2 ψ 2 ) u + 2 E 2

and

G 12 = i = 1 2 ( K γ ρ + u + 2 ρ i 2 ( ρ ˜ γ 1 ρ + γ 3 ρ i 2 ) φ i 2 + K ρ + u + ρ ˜ γ 1 [ Φ ( ρ ˜ ρ i ) γ 2 ( ρ ˜ ρ i 1 ) 2 ] + K γ ( ρ ˜ γ 1 ρ + γ 1 ) φ i ψ i + K ρ ˜ γ [ ( ρ i ρ ˜ ) γ 1 γ ( ρ i ρ ˜ 1 ) ] ψ i + ( ρ i u i ρ + u + ) E i ) 1 2 ( u ˜ u + ) E 2 .

The conditions M + > 1 and u + < 0 yield that the quadratic form G 11 is positive definite since

G 11 = i = 1 2 [ ( P ( ρ + ) 3 2 2 ρ + φ i 2 + ρ + ρ + 2 ψ i 2 ) ( M + 1 ) + P ( ρ + ) 2 ρ + ( P ( ρ + ) φ i ρ + ψ i ) 2 ] 1 2 u + E 2 C ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 + E 2 ) ,

which yields

0 t R + W x G 11 d x d τ C 0 t χ ( τ ) | ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) | 2 , ω x 2 d τ . (3.10)

Using (2.5), (3.5) and the inequalities | Φ ( s ) r 2 ( s 1 ) 2 | C | s 1 | 2 , | s γ 1 γ ( s 1 ) | C | s 1 | 2 for | s 1 | 1 , we have the estimate for G 12 as

| G 12 | C ( N ( t ) + δ ) ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 + E 2 ) ,

which implies

0 t R + W x G 12 d x d τ C ( N ( t ) + δ ) 0 t χ ( τ ) | ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) | 2 , ω x 2 d τ . (3.11)

Moreover, the positive bound of ρ i ( i = 1 , 2 ), (3.6) and the Schwarz inequality yield the estimate for R 1 as

| R 1 | C | u ˜ x | ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 + E 2 ) ,

then we have

0 t R + W R 1 d x d τ C 0 t χ ( t ) | ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) | 2 , | u ˜ x | ω 2 d τ . (3.12)

Therefore, integrating (3.7) over R + × ( 0 , t ) , substituting the above inequalities (3.8)-(3.12) into the resultant equality and then taking N ( T ) + δ suitably small, we obtain the desired estimate (3.3). □

Next, we obtain the estimate for the first-order derivatives of the solution for (2.7)-(2.9). As the existence of the higher-order derivatives of the solution is not supposed, we need to use the difference quotient for the rigorous derivation of the higher-order estimates. Since the argument using the difference quotient is similar to that in the paper [21,22], we omit the details and proceed with the proof as if it verifies

( φ 1 , ψ 1 , φ 2 , ψ 2 ) C ( [ 0 , T ] ; H 2 ( R + ) ) , φ 1 x , φ 2 x L 2 ( [ 0 , T ] ; H 1 ( R + ) ) , E x L 2 ( [ 0 , T ] ; L 2 ( R + ) ) , ψ 1 x , ψ 2 x L 2 ( [ 0 , T ] ; H 2 ( R + ) ) .

Lemma 3.3There exists a positive constant ε 0 such that if N ( T ) + δ < ε 0 , then

χ ( t ) ( φ 1 x , φ 2 x ) 2 + 0 t χ ( τ ) ( ( φ 1 x , φ 2 x , E x ) 2 + φ 1 x ( τ , 0 ) 2 + φ 2 x ( τ , 0 ) 2 ) d τ C ( ( φ 10 , φ 20 ) x 2 + | ( φ 10 , φ 20 , ψ 10 , ψ 20 , E 0 ) | 2 , ω 2 ) + C L ( t ) 2 + C δ M ( t ) 2 . (3.13)

Proof By differentiating the first and third equations of (2.7) in x, and then multiplying them by φ 1 x ρ 1 3 and φ 2 x ρ 2 3 , respectively, one has for i = 1 , 2 ,

( φ i x 2 2 ρ i 3 ) t + ( u i φ i x 2 2 ρ i 3 ) x u ˜ x φ i x 2 ρ i 3 + ρ ˜ x ρ i 3 φ i x ψ i x + 1 ρ i 2 φ i x ψ i x x = f i x φ i x ρ i 3 ,

which yields

( φ 1 x 2 2 ρ 1 3 + φ 2 x 2 2 ρ 2 3 ) t + ( u 1 φ 1 x 2 2 ρ 1 3 + u 2 φ 2 x 2 2 ρ 2 3 ) x + 1 ρ 1 2 φ 1 x ψ 1 x x + 1 ρ 2 2 φ 2 x ψ 2 x x = R 21 (3.14)

with R 21 = u ˜ x φ 1 x 2 ρ 1 3 ρ ˜ x ρ 1 3 φ 1 x ψ 1 x + u ˜ x φ 2 x 2 ρ 2 3 ρ ˜ x ρ 2 3 φ 2 x ψ 2 x f 1 x φ 1 x ρ 1 3 f 2 x φ 2 x ρ 2 3 . On the other hand, multiplying the second and fourth equations of (2.7) by φ 1 x ρ 1 2 and φ 2 x ρ 2 2 , respectively, gives

( ψ i ρ i φ i x ) t ( ψ i ρ i φ i t + ρ ˜ x ψ i 2 ρ i ) x ψ i x 2 + P ( ρ i ) ρ i 2 φ i x 2 1 ρ i u ˜ x φ i ψ i x + ρ ˜ x x ρ i ψ i 2 + 2 ρ ˜ x ρ i ψ i ψ i x + u ˜ x ( ρ ˜ x φ i ρ ˜ φ i x ) ψ i ρ i 2 = g i φ i x ρ i 2 + ( 1 ) i 1 E ρ i φ i x + φ i x ψ i x x ρ i 2 , i = 1 , 2 .

Further, we have

( ψ 1 ρ 1 φ 1 x + ψ 2 ρ 2 φ 2 x ) t ( ψ 1 ρ 1 φ 1 t + ρ ˜ x ψ 1 2 ρ 1 + ψ 2 ρ 2 φ 2 t + ρ ˜ x ψ 2 2 ρ 2 ) x ψ 1 x 2 ψ 2 x 2 + P ( ρ 1 ) ρ 1 2 φ 1 x 2 + P ( ρ 2 ) ρ 2 2 φ 2 x 2 = E ρ 1 φ 1 x E ρ 2 φ 2 x + R 22 + φ 1 x ψ 1 x x ρ 1 2 + φ 2 x ψ 2 x x ρ 2 2 , (3.15)

here

R 22 = 1 ρ 1 u ˜ x φ 1 ψ 1 x ρ ˜ x x ρ 1 ψ 1 2 2 ρ ˜ x ρ 1 ψ 1 ψ 1 x u ˜ x ( ρ ˜ x φ 1 ρ ˜ φ 1 x ) ψ 1 ρ 1 2 g 1 φ 1 x ρ 1 2 + 1 ρ 2 u ˜ x φ 2 ψ 2 x ρ ˜ x x ρ 2 ψ 2 2 u ˜ x ( ρ ˜ x φ 2 ρ ˜ φ 2 x ) ψ 2 ρ 2 2 g 2 φ 2 x ρ 2 2 .

Combining (3.14) and (3.15), we have

( φ 1 x 2 2 ρ 1 3 + φ 1 x ψ 1 ρ 1 + φ 2 x 2 2 ρ 2 3 + φ 2 x ψ 2 ρ 2 ) t + ( u 1 φ 1 x 2 2 ρ 1 3 φ 1 t ψ 1 ρ 1 ρ ˜ x ψ 1 2 ρ 1 + u 2 φ 2 x 2 2 ρ 2 3 φ 2 t ψ 2 ρ 2 ρ ˜ x ψ 2 2 ρ 2 ) x + P ( ρ 1 ) φ 1 x 2 ρ 1 2 + P ( ρ 2 ) φ 2 x 2 ρ 2 2 = ψ 1 x 2 + ψ 2 x 2 + E ( φ 1 x ρ 1 φ 2 x ρ 2 ) + R 21 + R 22 . (3.16)

The second term on the right-hand side of (3.16) can be rewritten as

E ( φ 1 x ρ 1 φ 2 x ρ 2 ) = [ ( ln ρ 1 ln ρ 2 ) x ] E ( ρ 1 ρ 2 ) E ρ ˜ x ρ 1 ρ 2 = [ ( ln ρ 1 ln ρ 2 ) E ] x ( ln ρ 1 ln ρ 2 ) E x + E x E ρ ˜ x ρ 1 ρ 2 .

Under the assumption (3.6) on the densities, it holds that

( ln ρ i ln ρ e ) E x C E x 2 , | E x E ρ ˜ x ρ 1 ρ 2 | C | ρ ˜ x | ( E 2 + E x 2 ) . (3.17)

Moreover, owing to the Schwarz inequality with the aid of (2.5), R 21 is estimated as

R 21 | u ˜ x | ( φ 1 x 2 + ψ 1 x 2 + φ 2 x 2 + ψ 2 x 2 ) + | ρ ˜ x | ( φ 1 x 2 + ψ 1 x 2 + φ 2 x 2 + ψ 2 x 2 ) + | ρ ˜ x x | ( ψ 1 2 + φ 1 x 2 + ψ 2 2 + φ 2 x 2 ) + | u ˜ x x | ( φ 1 2 + φ 1 x 2 + φ 2 2 + φ 2 x 2 ) C δ ( φ 1 2 + φ 2 2 + ψ 1 2 + ψ 2 2 + φ 1 x 2 + φ 2 x 2 + ψ 1 x 2 + ψ 2 x 2 ) . (3.18)

Similarly, we have

R 22 ε ( φ 1 x 2 + φ 2 x 2 ) + C ε ( ψ 1 x 2 ψ 1 2 + ψ 1 4 + ψ 2 x 2 ψ 2 2 + ψ 2 4 ) + C δ ( φ 1 2 + φ 2 2 + ψ 1 2 + ψ 2 2 + φ 1 x 2 + φ 2 x 2 + ψ 1 x 2 + ψ 2 x 2 ) . (3.19)

Multiplying (3.16) by a weight function χ ( t ) , we get

(3.20)

The boundary condition (2.9) gives

0 [ χ ( t ) i = 1 2 ( u i φ i x 2 2 ρ i 3 φ i t ψ i ρ i ρ ˜ ψ i 2 ρ i ) ] x d x C χ ( t ) ( φ 1 x ( t , 0 ) 2 + φ 2 x ( t , 0 ) 2 ) (3.21)

and

0 χ [ ( ln ρ 1 ln ρ 2 ) E ] x d x C χ ( φ 1 ( t , 0 ) 2 + φ 2 ( t , 0 ) 2 + E ( t , 0 ) 2 ) . (3.22)

Integrating (3.20) over [ 0 , t ] × R + , substituting (3.17), (3.18), (3.19), (3.21), (3.22) and the estimate

| φ 1 ( t , x ) | + | φ 2 ( t , x ) | + | ψ 1 ( t , x ) | + | ψ 2 ( t , x ) | + | E ( t , x ) | | φ 1 ( t , 0 ) | + | φ 2 ( t , 0 ) | + | E ( t , 0 ) | + x ( φ 1 x , φ 2 x , ψ 1 x , ψ 2 x , E x ) ( t ) , (3.23)

which is proved by the similar computation as in [21-24], in the resultant equality, and take ε and δ suitably small. These computations together with (3.3) give the desired estimate (3.13). □

Lemma 3.4There exists a positive constant ε 0 such that if N ( T ) + δ < ε 0 , then

χ ( t ) ( ψ 1 x , ψ 2 x ) 2 + 0 t χ ( τ ) ( ψ 1 x x , ψ 2 x x ) 2 d τ C ( ( ψ 10 , ψ 20 ) x 2 + | ( φ 10 , φ 20 , ψ 10 , ψ 20 , E 0 ) | 2 , ω 2 ) + C L ( t ) 2 + C ( N ( t ) + δ ) M ( t ) 2 . (3.24)

Proof Multiplying (2.7)2 by ψ 1 x x ρ 1 , and (2.7)4 by ψ 2 x x ρ 2 , respectively, we have for i = 1 , 2 ,

( ψ i x 2 2 ) t ( ψ i t ψ i x + u i ψ i x 2 2 ) x + ψ i x x 2 ρ i = ψ i x 3 2 u ˜ x ψ i x 2 2 + P ( ρ i ) ρ i φ i x ψ i x x + g i ψ i x x ρ i + ( 1 ) i E ψ i x x ,

which yields

( ψ 1 x 2 2 + ψ 2 x 2 2 ) t i = 1 2 ( ψ i t ψ i x + u i ψ i x 2 2 ) x + ψ 1 x x 2 ρ 1 + ψ 2 x x 2 ρ 2 = E ( ψ 1 x x ψ 2 x x ) + R 3 (3.25)

with

R 3 = ψ 1 x 3 2 u ˜ x ψ 1 x 2 2 + P ( ρ 1 ) ρ 1 φ 1 x ψ 1 x x ψ 2 x 3 2 u ˜ x ψ 2 x 2 2 + P ( ρ 2 ) ρ 2 φ 2 x ψ 2 x x + g 1 ψ 1 x x ρ 1 + g 2 ψ 2 x x ρ 2 .

Note that E ψ 1 x x + E ψ 2 x x = ( E ( ψ 1 ψ 2 ) x ) x + E x ( ψ 1 ψ 2 ) x , and the function R 3 is estimated by using (2.5) and Schwarz inequality as

R 3 ε ( ψ 1 x x 2 + ψ 2 x x 2 ) + C ε ( ψ 1 x 2 + ψ 2 x 2 + φ 1 x 2 + φ 2 x 2 + ψ 1 x 4 + ψ 2 x 4 ) + C ε | u ˜ x | ( φ 1 2 + φ 2 2 + ψ 1 2 + ψ 2 2 + ψ 1 x 2 + ψ 2 x 2 ) + C ε | ρ ˜ x | ( φ 1 2 + φ 2 2 ) , (3.26)

where ε is an arbitrary positive constant and C ε is a positive constant depending on ε. Then, multiplying (3.25) by a weight function χ ( t ) , we get

i = 1 2 [ ( χ ψ i x 2 2 ) t ( χ ( ψ i t ψ i x + u i ψ i x 2 2 ) ) x + χ ψ i x x 2 ρ i ] = χ t i = 1 2 ψ i x 2 4 χ [ E ( ψ 1 x x ψ 2 x x ) R 3 ] . (3.27)

Integrate (3.27) over [ 0 , t ] × R + , substitute (3.26) as well as the estimate

0 ( ψ 1 x 4 + ψ 2 x 4 ) d x C ( ψ 1 x 1 2 ψ 1 x 2 + ψ 2 x 1 2 ψ 2 x 2 ) C N ( t ) ( ψ 1 x , ψ 1 x x , ψ 2 x , ψ 2 x x ) 2

and

0 t χ | E ( ψ 1 ψ 2 ) x ( 0 , τ ) | d τ 0 t χ | E ( 0 , τ ) | 2 d τ + 0 t χ | ( ψ 1 ψ 2 ) x ( 0 , τ ) | 2 d τ 0 t χ | E ( 0 , τ ) | 2 d τ + C ε 0 t χ ( ψ 1 , ψ 2 ) x 2 d τ + ε 0 t χ ( ψ 1 , ψ 2 ) x x 2 d τ

in the resultant equality, and take ε suitably small. These computations together with (3.3), (3.13) and (3.23) give the desired estimate (3.24). □

Proof of Proposition 3.1 Summing up the estimates (3.3), (3.13) and (3.24) and taking N ( T ) + δ suitably small, we have

(3.28)

First, we prove the estimate (3.1). Noting the Poincaré-type inequality (3.23), and substituting ω ( x ) = ( 1 + x ) β and χ ( t ) = ( 1 + t ) ξ in (3.28) for β [ 0 , α ] and ξ 0 gives

(3.29)

Therefore, applying an induction to (3.29) gives the desired estimate (3.1). Since this computation is similar those in [23,25], we omit the details.

Next, we prove the estimate (3.2). Substitute ω ( x ) = e β x and χ ( t ) = e α t in (3.28) for β < λ to obtain

(3.30)

Here, we have used the Poincaré-type inequality (3.23) again. Thus, taking δ, β and α suitably small, we obtain the desired a priori estimate (3.2). □

### 4 A priori estimate for M + = 1

In the section we proceed to consider the transonic case M + = 1 . To state the a priori estimate of the solution precisely, here we use the notations:

N 1 ( t ) = sup 0 τ t ( 1 + x ) α / 2 ( φ 1 , ψ 1 , φ 2 , ψ 2 ) ( τ ) 1 , M 1 ( t ) 2 = 0 t ( 1 + t ) ξ ( φ 1 x , φ 2 x , E x , ψ 1 x , ψ 2 x , ψ 1 x x , ψ 2 x x ) ( τ ) a , β 2 d τ .

Proposition 4.1Suppose that the same assumption as in Theorem 1.1 holds. Suppose that ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) X ( 1 + x ) α 1 ( 0 , T ) is a solution to (2.7)-(2.9) for certain positive constantsαandT. Then there exist positive constants ε 0 andCsuch that if N 1 ( T ) + δ ε 0 , then the solution ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) satisfies the estimate

( 1 + t ) α / 2 + ε ( ( φ 1 , ψ 1 , φ 2 , ψ 2 ) ( t ) 1 2 + E ( t ) 2 ) + 0 t ( 1 + τ ) α / 2 + ε | ( φ 1 , φ 2 , φ 1 x , φ 2 x , E ) ( τ , 0 ) | 2 d τ + 0 t ( 1 + τ ) α / 2 + ε ( ( φ 1 x , φ 2 x , E x ) ( τ ) 2 + ( ψ 1 x , ψ 2 x ) ( τ ) 1 2 ) d τ C ( ( φ 10 , ψ 10 , φ 20 , ψ 20 ) x 2 + ( φ 10 , ψ 10 , φ 20 , ψ 20 , E 0 ) a , α 2 ) ( 1 + t ) ε . (4.1)

In order to prove Proposition 4.1, we need to get a lower estimate for u ˜ x and the Mach number M ˜ on the stationary solution ( ρ ˜ , u ˜ ) defined by M ˜ ( x ) : = | u ˜ ( x ) | P ( ρ ˜ ( x ) ) .

Lemma 4.2 (see [23])

The stationary solution u ˜ ( x ) satisfies

u ˜ x ( x ) A ( u + u b ) γ + 2 δ 2 ( 1 + B x ) 2 , A : = ( γ + 1 ) ρ + 2 , B : = δ A

for x ( 0 , ) . Moreover, there exists a positive constantCsuch that

γ + 1 2 | u + | δ 1 + B x C δ 2 ( 1 + B x ) 2 M ˜ ( x ) 1 C δ 1 + B x .

Based on Lemma 4.2, we obtain the weighted L 2 estimate of ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) .

Lemma 4.3There exists a positive constant ε 0 such that if N 1 ( T ) + δ < ε 0 , then

(4.2)

for β [ 0 , α ] and ξ 0 .

Proof First, from (2.7), similar as (3.4), we also have

(4.3)

Notice that from the second and fourth equations of (1.1), one has

2 E = t ( u 1 u 2 ) + 1 2 x ( u 1 2 u 2 2 ) + [ P ( ρ 1 ) x ρ 1 ρ 1 P ( ρ 2 ) x ρ 2 ρ 2 ] [ u 1 x x ρ 1 u 2 x x ρ 2 ] ,

which implies

u ˜ x 2 E 2 = ( u ˜ x 4 ( u 1 u 2 ) E ) t u ˜ x 4 ( u 1 u 2 ) E t + u ˜ x 8 E x ( u 1 2 u 2 2 ) + A γ u ¯ x 4 ( γ 1 ) ( ρ 1 γ 1 ρ 2 γ 1 ) x E + u ˜ x 4 E [ u 1 x x ρ 1 u 2 x x ρ 2 ] = ( u ˜ x 4 ( u 1 u 2 ) E ) t + ( u ˜ u ˜ x 4 E ( ψ 1 ψ 2 ) ) x + u ˜ 8 ( ρ 1 + ρ 2 ) ( ψ 1 ψ 2 ) 2 u ˜ u ˜ x x 4 E ( ψ 1 ψ 2 ) + u ˜ x 8 E x ( ψ 1 2 ψ 2 2 ) + A γ u ¯ x 4 ( γ 1 ) ( ρ 1 γ 1 ρ 2 γ 1 ) x E + u ˜ x 4 E [ ( ψ 1 2 2 ψ 2 2 2 ) x + ( u 1 x x ρ 1 u 2 x x ρ 2 ) ] .

Plugging the above equality into (4.3), we arrive at

(4.4)

where

R 4 = u ˜ u ˜ x x 4 E ( ψ 1 ψ 2 ) + u ˜ x 8 E x ( ψ 1 2 ψ 2 2 ) + u ˜ x 4 E ( ψ 1 ψ 1 x ψ 2 ψ 2 x ) + A γ u ˜ x 4 ( γ 1 ) ( ρ 1 γ 1 ρ 2 γ 1 ) x E + u ˜ x 4 [ u 1 x x ρ 1 u 2 x x ρ 2 ] E = : R 41 + R 42 + R 43 + R 44 + R 45 .

Further, multiplying (4.4) by a weight function W ( t , x ) : = ( 1 + B x ) β ( 1 + t ) ξ , we have

(4.5)

Where G 2 and G 3 are defined

G 2 = ρ 1 u 1 E 1 ρ 2 u 2 E 2 ( P ( ρ 1 ) P ( ρ ˜ ) ) ψ 1 ( P ( ρ 2 ) P ( ρ ˜ ) ) ψ 2 u ˜ u ˜ x 4 ( ψ 1 ψ 2 ) E u ˜ 2 E 2

and

G 3 : = W u ˜ x [ P ( ρ 1 ) P ( ρ ¯ ) P ( ρ ˜ ) φ 1 + ρ 1 ψ 1 2 ρ 1 8 ( ψ 1 ψ 2 ) 2 + P ( ρ 2 ) P ( ρ ¯ ) P ( ρ ˜ ) φ 2 + ρ 2 ψ 2 2 ρ 2 8 ( ψ 1 ψ 2 ) 2 ] 1 2 W x x ψ 1 2 1 2 W x x ψ 2 2 .

By the same computation as in deriving (3.9), we rewrite the terms G 2 and G 3 to G 2 = G 21 + G 22 , G 3 = G 31 + G 32 with

G 21 = i = 1 2 [ ( P ( ρ + ) 3 / 2 2 ρ + φ i 2 + ρ + P ( ρ + ) 2 ψ i 2 ) ( M ˜ 1 ) + P ( ρ ˜ ) 2 ρ ˜ ( P ( ρ ˜ ) φ i ρ ˜ ψ i ) 2 ] u ˜ u ˜ x 4 ( ψ 1 ψ 2 ) E 1 2 u + E 2 , G 22 = i = 1 2 [ ρ ˜ P ( ρ ˜ ) u ˜ 2 ( 1 ρ i 2 1 ρ ˜ 2 ) φ i 2 + P ( ρ ˜ ) u ˜ [ Φ ( ρ ˜ ρ i ) γ 2 ( ρ ˜ ρ i 1 ) 2 ] + ( ρ i u i ρ ˜ u ˜ ) E i + ( P ( ρ i ) P ( ρ ˜ ) P ( ρ ˜ ) φ i ) ψ i [ ( P ( ρ ˜ ) 3 / 2 2 ρ ˜ P ( ρ + ) 3 / 2 2 ρ + ) φ i 2 + ( ρ ˜ P ( ρ ˜ ) 2 ρ + P ( ρ + ) 2 ) ψ i 2 ] ( M ˜ 1 ) ] 1 2 ( u ˜ u + ) E 2 , G 31 = W u ˜ x ( ρ + ψ 1 2 + 1 2 P ( ρ + ) φ 1 2 ) 1 2 W x x ψ 1 2 + W u ˜ x ( ρ + ψ 2 2 + 1 2 P ( ρ + ) φ 2 2 ) 1 2 W x x ψ 2 2

and

G 32 = W u ˜ x i = 1 2 [ ( ρ i ρ + ) ψ i 2 + 1 2 ( P ( ρ ˜ ) P ( ρ + ) ) φ i 2 + P ( ρ i ) P ( ρ ˜ ) P ( ρ ˜ ) φ i 1 2 P ( ρ ˜ ) φ i 2 ] .

By utilizing Lemma 4.2 with the aid of the fact that β < α and u + < 0 , we obtain the lower estimate of W x G 21 + G 31 as

W x G 21 + G 31 i = 1 2 [ K γ ρ + γ 2 A 4 [ ( γ + 1 ) β + 2 ( u + u b ) γ + 2 ( γ 1 ) ] δ 2 ( 1 + t ) ξ ( 1 + B x ) β 2 φ i 2 + ρ + A 4 [ 4 ( u + u b ) γ + 2 ( γ + 1 ) β ( β 2 ) ] δ 2 ( 1 + t ) ξ ( 1 + B x ) β 2 ψ i 2 ] C β δ 3 ( 1 + t ) ξ ( 1 + B x ) β 3 ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 ) C δ 2 ( 1 C δ ) ( 1 + t ) ξ ( 1 + B x ) β 2 ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 ) + C β ( 1 + t ) ξ ( 1 + B x ) β 1 E 2 (4.6)

for β ( 0 , α ] . On the other hand, the estimates (2.6), (3.5) and (3.6) yield

| W x G 22 + G 32 | C ( N 1 ( t ) + δ 2 ) δ ( 1 + t ) ξ ( 1 + B x ) β 2 ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 ) + C δ β ( 1 + t ) ξ ( 1 + B x ) β 1 E 2 . (4.7)

For the first term on the right-hand side of (4.5), we estimate it as

(4.8)

Similarly, we get

(4.9)

In the same way, we estimate R 42 and R 43 as follows:

| 0 t R + W R 42 d x d τ | C δ 0 t ( 1 + τ ) ξ R + ( 1 + B x ) β ( ψ 1 x 2 + ψ 2 x 2 + E x 2 ) d x d τ (4.10)

and

(4.11)

Since ( ρ 1 γ 1 ρ 2 γ 1 ) x = ( γ 1 ) [ ( ρ 1 γ 2 ρ 2 γ 2 ) ρ ˜ x + ρ 1 γ 2 φ 1 x ρ 2 γ 2 φ 2 x ] , we have

(4.12)

Moreover, it is easy to compute R 45 = u ˜ x ( E ψ 1 x x ρ 1 E ψ 2 x x ρ 2 ) + u ˜ x x u ˜ x ρ 1 ρ 2 E ( φ 1 φ 2 ) , which implies

(4.13)

Finally, integrate (4.5) over R + × ( 0 , t ) , substitute (4.6)-(4.13) in the resultant equality, and take N 1 ( t ) and δ suitably small. This procedure yields the desired estimate (4.2) for β ( 0 , α ] .

Next, we prove (4.2) for β = 0 . Substituting W = ( 1 + t ) ξ in (4.5) and integrating the resultant equality over R + × ( 0 , t ) , we get

( 1 + t ) ξ ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) ( t ) 2 + 0 t ( 1 + τ ) ξ ( i = 1 2 φ i ( τ , 0 ) 2 + E ( τ , 0 ) 2 + ( ψ 1 x , ψ 2 x ) 2 ) d τ C ( φ 10 , ψ 10 , φ 10 , ψ 10 , E 0 ) 2 + C ξ 0 t ( 1 + τ ) ξ 1 ( φ 1 , ψ 1 , φ 2 , ψ 2 , E ) 2 d τ + C δ 0 t ( 1 + τ ) ξ ( φ 1 x , φ 2 x , E x , ψ 1 x x , ψ 2 x x ) 2 d x d τ .

Here, we have used the fact that G 2 0 holds. Therefore, we obtain the estimate (4.2) for the case of β = 0 . □

In order to complete the proof of Proposition 4.1, we need to obtain the weighted estimate of ( φ 1 x , ψ 1 x , φ 2 x , ψ 2 x ) .

Lemma 4.4There exists a positive constant ε 0 such that if N 1 ( T ) + δ < ε 0 , then for β [ 0 , α ] and ξ 0 ,

(4.14)

Proof Since the derivation of the estimate (4.14) is similar to that of (3.13) and (3.24), we only give the outline of the proof. Multiplying (3.16) by W = ( 1 + B x ) β ( 1 + t ) ξ , we have

i = 1 2 [ [ W ( φ i x 2 2 ρ i 3 + φ i x ψ i ρ i ) ] t + [ W ( u i φ i x 2 2 ρ i 3