Skip to main content

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

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.

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 [14].

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 [57]. 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 [1315]. 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 [1820] 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.1 Suppose 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)
  1. (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)
  1. (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)dy.

Notations Throughout this paper, C>0 denotes the generic positive constant independent of time. L p (R) (1p<) 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) (k0) stands for the space of L 2 (R)-functions f whose derivatives (in the sense of distribution) D x l f (lk) 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,l0) 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 λ and C such 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.3 If the initial data satisfies (1.8) and ω φ 10 , ω ψ 10 , ω φ 20 , ω ψ 20 , ω E 0 L 2 ( R + ), there exists a positive constant T such 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.1 Suppose that the same assumptions as in Theorem  1.1 hold.

  1. (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 α and T. Then there exist positive constants ε 0 and C such 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.

  1. (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 ζ and T. 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.2 Suppose 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 dxdτ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γ(s1)|C | s 1 | 2 for |s1|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 dxdτ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 dxdτ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.3 There 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 dxCχ(t) ( φ 1 x ( t , 0 ) 2 + φ 2 x ( t , 0 ) 2 )
(3.21)

and

0 χ [ ( ln ρ 1 ln ρ 2 ) E ] x dxCχ ( φ 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 [2124], in the resultant equality, and take ε and δ suitably small. These computations together with (3.3) give the desired estimate (3.13). □

Lemma 3.4 There 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.1 Suppose 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 α and T. Then there exist positive constants ε 0 and C such 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 constant C such that

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

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

Lemma 4.3 There 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

2E= 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 ) dxdτ
(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.4 There 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 φ i t ψ i ρ i ρ ˜ x ψ i 2 ρ i ) ] x + W P ( ρ i ) φ i x 2 ρ i 2 ] = i = 1 2 [ W t ( φ i x 2 2 ρ i 3 + φ i x ψ i ρ i ) + W ( ψ i x 2 + ( 1 ) i + 1 E φ i x ρ i + R 21 + R 22 ) W x ( u i φ i x 2 2 ρ i 3 φ i t ψ i ρ i ρ ˜ x ψ i 2 ρ i ) ] .
(4.15)

Integrating (4.15) over R + ×(0,t) and substituting (4.2) gives the estimate for φ 1 x and φ 2 x as

(4.16)

Here, we have used the inequalities

( ( 1 + B x ) β ) x | 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 | C δ ( 1 + B x ) β ( φ 1 x 2 + φ 2 x 2 ) + C ( 1 + B x ) β ( ψ 1 x 2 + ψ 2 x 2 ) + C β ( 1 + B x ) β 2 ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 )

and

( 1 + B x ) β | R 21 + R 22 | ( ε + C δ ) ( 1 + B x ) β ( φ 1 x 2 + φ 2 x 2 ) + C ε ( 1 + B x ) β ( ψ 1 x 2 + ψ 2 x 2 ) + C δ ( 1 + B x ) β ( φ 1 2 + ψ 1 2 + φ 2 2 + ψ 2 2 ) ,

where ε is an arbitrary positive constant. We note that the third term on the right-hand side of the above inequality is estimated by applying the Poincaré-type inequality (3.23) for the case of β=0.

Next, we prove the estimate for ( ψ 1 x , ψ 2 x ). Multiply (3.25) by W= ( 1 + t ) ξ ( 1 + B x ) β to get

( W ψ 1 x 2 2 + W ψ 2 x 2 2 ) t i = 1 2 ( W ψ i t ψ i x + W u i ψ i x 2 2 ) x + W i = 1 2 ψ i x x 2 ρ i = 1 2 W t i = 1 2 ψ i x 2 2 W x i = 1 2 ( ψ i t ψ i x + u i ψ i x 2 2 ) E W ψ 1 x x + E W ψ 2 x x + W R 3 .
(4.17)

Integrate (4.17) in R + ×(0,t) and substitute (4.2) and (4.16) in the resultant equality with the inequalities

( ( 1 + B x ) β ) x | ψ 1 t ψ 1 x + u 1 ψ 1 x 2 2 + ψ 2 t ψ 2 x + u 2 ψ 2 x 2 2 | + ( 1 + B x ) β | R 3 | ε ( 1 + B x ) β i = 1 2 ψ i x x 2 + C ε ( 1 + B x ) β i = 1 2 ( φ i x 2 + ψ i x 2 + ψ i x 4 ) + C ε δ ( 1 + B x ) β 4 i = 1 2 ( φ i 2 + ψ i 2 )

and

0 ( 1 + B x ) β ( ψ 1 x 4 + ψ 2 x 4 ) dxC N 1 (t) ( ψ 1 x , ψ 1 x x , ψ 2 x , ψ 2 x x ) a , β 2 .

This procedure yields

( 1 + t ) ξ ( ψ 1 x , ψ 2 x ) a , β 2 + 0 t ( 1 + τ ) ξ ( ψ 1 x x , ψ 2 x x ) a , β 2 d τ C ( ( φ 10 , ψ 10 , φ 20 , ψ 20 , E 0 ) a , β 2 + ( φ 10 x , ψ 10 x , φ 20 x , ψ 20 x ) a , β 2 ) + C ξ 0 t ( 1 + τ ) ξ 1 i = 1 2 ( φ i , ψ i , φ i x , ψ i x ) a , β 2 d τ + C ( N 1 ( t ) + δ ) M 1 ( t ) 2 .
(4.18)

Finally, adding (4.16) to (4.18) and taking N 1 (t)+δ suitably small give the desired estimate (4.14). □

By the same inductive argument as in deriving (4.1), we can prove Proposition 4.1, which immediately yields the decay estimate (1.12).

References

  1. Besse C, Claudel J, Degond P, et al.: A model hierarchy for ionospheric plasma modeling. Math. Models Methods Appl. Sci. 2004, 14: 393–415. 10.1142/S0218202504003283

    Article  MathSciNet  Google Scholar 

  2. Degond P: Mathematical modelling of microelectronics semiconductor devices. AMS/IP Stud. Adv. Math. 15. In Some Current Topics on Nonlinear Conservation Laws. Am. Math. Soc., Providence; 2000:77–110.

    Google Scholar 

  3. Jüngel A Progress in Nonlinear Differential Equations and Their Applications. In Quasi-hydrodynamic Semiconductor Equations. Birkhäuser, Basel; 2001.

    Chapter  Google Scholar 

  4. Sitnko A, Malnev V: Plasma Physics Theory. Chapman & Hall, London; 1995.

    Google Scholar 

  5. Ducomet B: A remark about global existence for the Navier-Stokes-Poisson system. Appl. Math. Lett. 1999, 12: 31–37.

    Article  MathSciNet  Google Scholar 

  6. Ducomet B: Local and global existence for the coupled Navier-Stokes-Poisson problem. Q. Appl. Math. 2003, 61: 345–361.

    Google Scholar 

  7. Zhang Y-H, Tan Z: On the existence of solutions to the Navier-Stokes-Poisson equations of a two-dimensional compressible flow. Math. Methods Appl. Sci. 2007, 30: 305–329. 10.1002/mma.786

    Article  MathSciNet  Google Scholar 

  8. Chan D:On the nonexistence of global weak solutions to the Navier-Stokes-Poisson equations in R N . Commun. Partial Differ. Equ. 2010, 35: 535–557. 10.1080/03605300903473418

    Article  Google Scholar 

  9. Hao C-C, Li H-L: Global existence for compressible Navier-Stokes-Poisson equations in three and higher dimensions. J. Differ. Equ. 2009, 246: 4791–4812. 10.1016/j.jde.2008.11.019

    Article  MathSciNet  Google Scholar 

  10. Li H-L, Matsumura A, Zhang G-J:Optimal decay rate of the compressible Navier-Stokes-Poisson system in R 3 . Arch. Ration. Mech. Anal. 2010, 196: 681–713. 10.1007/s00205-009-0255-4

    Article  MathSciNet  Google Scholar 

  11. Zhang G-J, Li H-L, Zhu C-J:Optimal decay rate of the non-isentropic compressible Navier-Stokes-Poisson system in R 3 . J. Differ. Equ. 2011, 250: 866–891. 10.1016/j.jde.2010.07.035

    Article  MathSciNet  Google Scholar 

  12. Wang W-K, Wu Z-G: Pointwise estimates of solution for the Navier-Stokes-Poisson equations in multi-dimensions. J. Differ. Equ. 2010, 248: 1617–1636. 10.1016/j.jde.2010.01.003

    Article  Google Scholar 

  13. Donatelli D, Marcati P: A quasineutral type limit for the Navier-Stokes-Poisson system with large data. Nonlinearity 2008, 21: 135–148. 10.1088/0951-7715/21/1/008

    Article  MathSciNet  Google Scholar 

  14. Ju Q-C, Li F-C, Li Y, Wang S: Rate of convergence from the Navier-Stokes-Poisson system to the incompressible Euler equations. J. Math. Phys. 2009., 50: Article ID 013533

    Google Scholar 

  15. Wang S, Jiang S: The convergence of the Navier-Stokes-Poisson system to the incompressible Euler equations. Commun. Partial Differ. Equ. 2006, 31: 571–591. 10.1080/03605300500361487

    Article  Google Scholar 

  16. Li H-L, Yang T, Zou C: Time asymptotic behavior of the bipolar Navier-Stokes-Poisson system. Acta Math. Sci., Ser. B 2009, 29: 1721–1736. 10.1016/S0252-9602(10)60013-6

    Article  MathSciNet  Google Scholar 

  17. Duan R-J, Yang X-F: Stability of rarefaction wave and boundary layer for outflow problem on the two-fluid Navier-Stokes-Poisson equations. Commun. Pure Appl. Anal. 2013, 12: 985–1014.

    Article  MathSciNet  Google Scholar 

  18. Gasser I, Hsiao L, Li H-L: Large time behavior of solutions of the bipolar hydrodynamical model for semiconductors. J. Differ. Equ. 2003, 192: 326–359. 10.1016/S0022-0396(03)00122-0

    Article  MathSciNet  Google Scholar 

  19. Huang F-M, Li Y-P: Large time behavior and quasineutral limit of solutions to a bipolar hydrodynamic model with large data and vacuum. Discrete Contin. Dyn. Syst. 2009, 24: 455–470.

    Article  MathSciNet  Google Scholar 

  20. Huang F-M, Mei M, Wang Y: Large time behavior of solution to n -dimensional bipolar hydrodynamic model for semiconductors. SIAM J. Math. Anal. 2011, 43: 1595–1630. 10.1137/100810228

    Article  MathSciNet  Google Scholar 

  21. Kawashima S, Nishibata S, Zhu P: Asymptotic stability of the stationary solution to the compressible Navier-Stokes-Poisson equations in half space. Commun. Math. Phys. 2003, 240: 483–500.

    Article  MathSciNet  Google Scholar 

  22. Kawashima S, Zhu P: Asymptotic stability of nonlinear wave for the compressible Navier-Stokes-Poisson equations in half space. J. Differ. Equ. 2008, 224: 3151–3179.

    Article  MathSciNet  Google Scholar 

  23. Nakamura T, Nishibata S, Yuge T: Convergence rate of solutions toward stationary solutions to the compressible Navier-Stokes equation in a half line. J. Differ. Equ. 2007, 241: 94–111. 10.1016/j.jde.2007.06.016

    Article  MathSciNet  Google Scholar 

  24. Nikkuni Y, Kawashima S: Stability of stationary solutions to the half-space problem for the discrete Boltzmann equation with multiple collisions. Kyushu J. Math. 2000, 54: 233–255. 10.2206/kyushujm.54.233

    Article  MathSciNet  Google Scholar 

  25. Matsumura A, Nishihara K: Asymptotic stability of traveling waves for scalar viscous conservation laws with non-convex nonlinearity. Commun. Math. Phys. 1994, 165: 83–96. 10.1007/BF02099739

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The research of Li is partially supported by the National Science Foundation of China (Grant No. 11171223) and the Innovation Program of Shanghai Municipal Education Commission (Grant No. 13ZZ109).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yeping Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors typed, read and approved the final manuscript.

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Zhou, F., Li, Y. Convergence rate of solutions toward stationary solutions to the bipolar Navier-Stokes-Poisson equations in a half line. Bound Value Probl 2013, 124 (2013). https://doi.org/10.1186/1687-2770-2013-124

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-124

Keywords