Skip to main content

Global behavior of 1D compressible isentropic Navier-Stokes equations with a non-autonomous external force

Abstract

In this paper, we study a free boundary problem for compressible Navier-Stokes equations with density-dependent viscosity and a non-autonomous external force. The viscosity coefficient μ is proportional to ρθwith 0 < θ < 1, where ρ is the density. Under certain assumptions imposed on the initial data and external force f, we obtain the global existence and regularity. Some ideas and more delicate estimates are introduced to prove these results.

1 Introduction

We study a free boundary problem for compressible Navier-Stokes equations with density-dependent viscosity and a non-autonomous external force, which can be written in Eulerian coordinates as:

ρ τ + ( ρ u ) ξ = 0 , τ > 0
(1.1)
( ρ u ) τ + ( ρ u 2 + P ( ρ ) ) ξ = ( μ u ξ ) ξ + ρ f , a ( τ ) < ξ < b ( τ )
(1.2)

with initial data

( ρ , u ) ( ξ , 0 ) = ( ρ 0 , u 0 ) ( ξ ) , a = a ( 0 ) ξ b ( 0 ) = b ,
(1.3)

where ρ = ρ (ξ,τ), u = u(ξ,τ), P = P(ρ) and f = f(ξ,t) denote the density, velocity, pressure and a given external force, respectively, μ = μ(ρ) is the viscosity coefficient. a(τ) and b(τ) are the free boundaries with the following property:

d d τ a ( τ ) = u ( a ( τ ) , τ ) , d d τ b ( τ ) = u ( b ( τ ) , τ ) ,
(1.4)
( - P ( ρ ) + μ ( ρ ) u ξ ) ( ξ , τ ) = 0 , ξ = a ( τ ) , b ( τ ) .
(1.5)

The investigation in [1] showed that the continuous dependence on the initial data of the solutions to the compressible Navier-Stokes equations with vacuum failed. The main reason for the failure at the vacuum is because of kinematic viscosity coefficient being independent of the density. On the other hand, we know that the Navier-Stokes equations can be derived from the Boltzmann equation through Chapman-Enskog expansion to the second order, and the viscosity coefficient is a function of temperature. For the hard sphere model, it is proportional to the square-root of the temperature. If we consider the isentropic gas flow, this dependence is reduced to the dependence on the density function by using the second law of thermal dynamics.

For simplicity of presentation, we consider only the polytropic gas, i.e. P(ρ) = γwith A > 0 being constants. Our main assumption is that the viscosity coefficient μ is assumed to be a functional of the density ρ, i.e. μ = θ, where c and θ are positive constants. Without loss of generality, we assume A = 1 and c = 1.

Since the boundaries x = a(τ) and x = b(τ) are unknown in Euler coordinates, we will convert them to fixed boundaries by using Lagrangian coordinates. We introduce the following coordinate transformation

x = a ( τ ) ξ ρ ( y , τ ) d y , t = τ ,
(1.6)

then the free boundaries ξ = a(τ) and ξ = b(τ) become

x = 0 , x = a ( τ ) b ( τ ) ρ ( z , τ ) d z = a b ρ 0 ( z ) d z
(1.7)

where a b ρ 0 ( z ) d z is the total initial mass, and without loss of generality, we can normalize it to 1. So in terms of Lagrangian coordinates, the free boundaries become fixed. Under the coordinate transformation, Eqs. (1.1)-(1.2) are now transformed into

ρ t + ρ 2 u x = 0 , t > 0 ,
(1.8)
u t + P ( ρ ) x = ( ρ μ ( ρ ) u x ) x + f ( r , t ) , 0 < x < 1
(1.9)

where r= 0 x ρ - 1 ( y , t ) d y. The boundary conditions (1.4)-(1.5) become

( - ρ γ + ρ 1 + θ u x ) ( d , t ) = 0 , d = 0 , 1 ,
(1.10)

and the initial data (1.3) become

( ρ , u ) ( x , 0 ) = ( ρ 0 , u 0 ) ( x ) , x [ 0 , 1 ] .
(1.11)

Now let us first recall some previous works in this direction. When the external force f ≡ 0, there have been many works (see, e.g., [29]) on the existence and uniqueness of global weak solutions, based on the assumption that the gas connects to vacuum with jump discontinuities, and the density of the gas has compact support. Among them, Liu et al. [4] established the local well-posedness of weak solutions to the Navier-Stokes equations; Okada et al. [5] obtained the global existence of weak solutions when 0 < θ < 1/3 with the same property. This result was later generalized to the case when 0 < θ < 1/2 and 0 < θ < 1 by Yang et al. [7] and Jiang et al. [3], respectively. Later on, Qin et al. [8, 9] proved the regularity of weak solutions and existence of classical solution. Fang and Zhang [2] proved the global existence of weak solutions to the compressible Navier-Stokes equations when the initial density is a piece-wise smooth function, having only a finite number of jump discontinuities.

For the related degenerated density function and viscosity coefficient at free boundaries, see Yang and Zhao [10], Yang and Zhu [11], Vong et al. [12], Fang and Zhang [13, 14], Qin et al. [15], authors studied the global existence and uniqueness under some assumptions on initial data.

When f ≠ 0, Qin and Zhao [16] proved the global existence and asymptotic behavior for γ = 1 and μ = const with boundary conditions u(0,t) = u(1,t) = 0; Zhang and Fang [17] established the global behavior of the Equations (1.1)-(1.2) with boundary conditions u(0,t) = ρ(1,t) = 0. In this paper, we obtain the global existence of the weak solutions and regularity with boundary conditions (1.4)-(1.5). In order to obtain existence and higher regularity of global solutions, there are many complicated estimates on external force and higher derivations of solution to be involved, this is our difficulty. To overcome this difficulty, we should use some proper embedding theorems, the interpolation techniques as well as many delicate estimates. This is the novelty of the paper.

The notation in this paper will be as follows:

L p ,1p+, W m , p ,mN, H 1 = W 1 , 2 , H 0 1 = W 0 1 , 2 denote the usual (Sobolev) spaces on [0,1]. In addition, || · || B denotes the norm in the space B; we also put ||||=||| | L 2 ( [ 0 , 1 ] ) .

The rest of this paper is organized as follows. In Section 2, we shall prove the global existence in H1. In Section 3, we shall establish the global existence in H2. In Section 4, we give the detailed proof of Theorem 4.1.

2 Global existence of solutions in H1

In this section, we shall establish the global existence of solutions in H1.

Theorem 2.1 Let 0 < θ < 1, γ > 1, and assume that the initial data (ρ0,u0) satisfies inf [ 0 , 1 ] ρ 0 >0, ρ 0 W 1 , 2 n , u 0 H 1 and external force f satisfies f(r(x,·),·) L2n([0,T], L2n[0,1]) for some n N satisfying n(2n - 1)/(2n2 + 2n - 1) > θ, then there exists a unique global solution (ρ (x,t),u(x,t)) to problem (1.8)-(1.11), such that for any T > 0,

0 < C 1 - 1 ( T ) ρ ( x , t ) C 1 ( T ) , ρ L ( [ 0 , T ] , H 1 [ 0 , 1 ] ) , u L ( [ 0 , T ] , H 1 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 2 [ 0 , 1 ] ) , u t L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) .

The proof of Theorem 2.1 can be done by a series of lemmas as follows.

Lemma 2.1 Under conditions of Theorem 2.1, the following estimates hold

0 1 1 2 u 2 + 1 γ - 1 ρ γ - 1 d x + 0 t 0 1 ρ 1 + θ u x 2 ( x , s ) d x d s C 1 ( T ) ,
(2.1)
ρ ( x , t ) C 1 ( T ) , ( x , t ) [ 0 , 1 ] × [ 0 , T ] ,
(2.2)
0 1 u 2 n d x + n ( 2 n - 1 ) 0 t 0 1 ρ 1 + θ u 2 n - 2 u x 2 ( x , s ) d x d s C 1 ( T )
(2.3)

where C1(T) denotes generic positive constant depending only on|| ρ 0 | | W 1 , 2 n [ 0 , 1 ] ,|| u 0 | | H 1 [ 0 , 1 ] , time T and||f| | L 2 n ( [ 0 , T ] , L 2 n [ 0 , 1 ] ) .

Proof Multiplying (1.8) and (1.9) by ργ-2 and u, respectively, using integration by parts, and considering the boundary conditions (1.10), we have

d d t 0 1 1 2 u 2 + 1 γ - 1 ρ γ - 1 d x + 0 1 ρ 1 + θ u x 2 d x = 0 1 f u d x
(2.4)

Integrating (2.4) with respect to t over [0,t], using Young's inequality, we have

0 1 1 2 u 2 + 1 γ - 1 ρ γ - 1 d x + 0 t 0 1 ρ 1 + θ u x 2 d x d s C 1 ( T ) + 1 2 0 t 0 1 u 2 d x + C 1 0 t 0 1 f 2 d x d s 1 2 0 t 0 1 u 2 d x + C 1 ( T )

which, by virtue of Gronwall's inequality and assumption f(r(x,·),·) L2n([0,T], L2n[0,1]), gives (2.1).

We derive from (1.8) that

( ρ θ ) t = - θ ρ 1 + θ u x
(2.5)

Integrating (2.5) with respect to t over [0,t] yields

ρ θ ( x , t ) = ρ 0 θ - θ 0 t ρ 1 + θ u x ( x , s ) d s .
(2.6)

Integrating (1.9) with respect to x, applying the boundary conditions (1.10), we obtain

ρ 1 + θ u x = 0 x u t d y + ρ γ - 0 x f ( r ( y , t ) , t ) d y
(2.7)

Inserting (2.7) into (2.6) gives

ρ θ + θ 0 t ρ γ d s = ρ 0 θ + θ 0 t 0 x f ( r , ( y , s ) , s ) d y d s - θ 0 x ( u - u 0 ) d y
(2.8)

Thus, the Hölder inequality and (2.1) imply

0 x u ( y , t ) d y C 1
(2.9)

and (2.2) follows from (2.8) and (2.9).

Multiplying (1.9) by 2nu2n-1and integrating over x and t, applying the boundary conditions (1.10), we have

0 1 u 2 n d x + 2 n ( 2 n - 1 ) 0 t 0 1 u 2 n - 2 ρ 1 + θ u x 2 d x d s = 0 1 u 0 2 n d x + 2 n ( 2 n - 1 ) 0 t 0 1 u 2 n - 2 ρ γ u x d x d s + 2 n 0 t 0 1 f u 2 n - 1 d x d s .
(2.10)

Applying the Young inequality and condition f(r(x, ·), ·) L2n([0,T],L2n[0,1]) to the last two terms in (2.10) yields

0 1 u 2 n d x + n ( 2 n - 1 ) 0 t 0 1 u 2 n - 2 ρ 1 + θ u x 2 d x d s C 1 + 0 t 0 1 f 2 n d x d s + ( 2 n - 1 ) 0 t 0 1 u 2 n d x d s + n ( 2 n - 1 ) 0 t 0 1 u 2 n - 2 ρ 2 γ - 1 - θ d x d s C 1 ( T ) + n ( 2 n - 1 ) 0 t 0 1 1 n ρ ( 2 γ - 1 - θ ) n + n - 1 n u 2 n d x d s + ( 2 n - 1 ) 0 t 0 1 u 2 n d x d s C 1 ( T ) + n ( 2 n - 1 ) 0 t 0 1 u 2 n d x d s .
(2.11)

Applying Gronwall's inequality, we conclude

0 1 u 2 n d x C 1 ( T )
(2.12)

, which, along with (2.11), yields (2.3). The proof of Lemma 2.1 is complete.

Lemma 2.2 Under conditions of Theorem 2.1, the following estimates hold

0 1 ( ρ θ ) x 2 n d x C 1 ( T ) ,
(2.13)
ρ ( x , t ) C 1 - 1 ( T ) > 0 .
(2.14)

Proof We derive from (2.5) and (1.9) that

( ρ θ ) x t = - θ ( u t + ( ρ γ ) x - f ) .
(2.15)

Integrating it with respect to t over [0,t], we obtain

( ρ θ ) x = ( ρ 0 θ ) x - θ ( u - u 0 ) - θ 0 t ( ρ γ ) x d s + θ 0 t f d s .
(2.16)

Multiplying (2.16) by [(ρθ) x ]2n-1, and integrating the resultant with respect to x to get

0 1 ( ρ θ ) x 2 n d x = 0 1 ( ρ θ ) x 2 n - 1 ( ρ 0 θ ) x d x - θ 0 1 ( u - u 0 ) + 0 t ( ρ γ ) x d s - 0 t f d s ( ρ θ ) x 2 n - 1 d x C 0 1 ( ρ θ ) x 2 n d x 2 n - 1 2 n 0 1 ( ρ 0 θ ) x 2 n d x 1 2 n + u - u 0 L 2 n + 0 1 0 t ( ρ γ ) x d s 2 n d x 1 2 n + 0 1 0 t f d s 2 n d x 1 2 n C 0 1 ( ρ θ ) x 2 n d x 2 n - 1 2 n 0 1 ( ρ 0 θ ) x 2 n d x 1 2 n + u - u 0 L 2 n + 0 t 0 1 ( ρ γ ) x 2 n d s 1 2 n d x + 0 t 0 1 f 2 n d x 1 2 n d s
(2.17)

here, we use the inequality g ( , s ) L p g ( , s ) L p d s . Using Young's inequality and assumptions of external of f, we get from (2.17) that

0 1 ( ρ θ ) x 2 n d x 1 2 0 1 ( ρ θ ) 2 n x d x + C 0 t 0 1 ( ρ γ ) x 2 n d x d s + C 0 t 0 1 f 2 n d x d s + C 1 ( T ) 1 2 0 1 ( ρ θ ) x 2 n d x + C 1 ( T ) 0 t 0 1 ( ρ γ ) x 2 n d x d s + C 1 ( T ) .
(30)

Hence,

0 1 ( ρ θ ) x 2 n d x C 1 ( T ) + C 1 ( T ) 0 t 0 1 ( ρ γ ) x 2 n d x d s
(2.18)

Using the Gronwall inequality to (2.18), we obtain (2.13).

The proof of (2.14) can be found in [3], please refer to Lemma 2.3 in [3] for detail.

Lemma 2.3 Under the assumptions in Theorem 2.1, for any 0 ≤ tT, we have the following estimate

u x ( t ) 2 + 0 t u t ( s ) 2 d s C 1 ( T ) .
(2.19)

Proof Multiplying (1.9) by u t , then integrating over [0,1] × [0,t], we obtain

0 t 0 1 u t 2 d x d s = 0 t 0 1 u t ( ρ 1 + θ u x - ρ γ ) x d x d s + 0 t 0 1 u t f d x d s .
(2.20)

Using integration by parts, (1.8) and the boundary conditions (1.10), we have

0 t 0 1 u 1 ρ 1 + θ u x - ρ γ x d x d s = 0 t 0 1 u t x ( ρ γ - ρ 1 + θ u x ) d x d s = 0 1 u x ρ γ - 1 2 ρ 1 + θ u x - u 0 x ρ 0 γ - 1 2 ρ 0 1 + θ u 0 x d x + 0 t 0 1 γ u x 2 ρ γ + 1 - 1 + θ 2 u x 3 ρ 2 + θ d x d s .

Thus,

0 t 0 1 u t 2 d x d s + 1 2 0 1 ρ 1 + θ u x 2 d x = 0 1 u x ρ γ - u 0 x ρ 0 γ - 1 2 ρ 0 1 + θ u 0 x d x + 0 t 0 1 γ u x 2 ρ γ + 1 - 1 + θ 2 u x 3 ρ 2 + θ d x d s + 0 t 0 1 u t f d x d s C 1 ( T ) + 0 1 1 4 ρ 1 + θ u x 2 + ρ 2 γ - 1 - θ d x + C 1 ( T ) 0 t sup [ 0 , 1 ] ρ γ - θ 0 1 ρ 1 + θ u x 2 d x d s + C 1 ( T ) 0 t 0 1 ρ 1 + θ | u x | 3 d x d s + 1 4 0 t 0 1 u t 2 d x d s + C 1 ( T ) 0 t 0 1 f 2 d x d s .

Using Lemmas 2.1-2.2, we derive

0 1 u x 2 d x + 0 t 0 1 u t 2 d x d s C 1 ( T ) + C 1 ( T ) 0 t 0 1 ρ 1 + θ | u x | 3 d x d s
(2.21)

The last term on the right-hand side of (2.21) can be estimated as follows, using (1.8), conditions (1.10) and Lemmas 2.1-2.2,

C 1 ( T ) 0 t 0 1 ρ 1 + θ | u x | 3 d x d s C 1 ( T ) 0 t max [ 0 , 1 ] | ρ 1 + θ u x | u x 2 d x d s C 1 ( T ) 0 t max [ 0 , 1 ] | ρ 1 + θ u x - ρ γ | 0 1 u x 2 d x d s + C 1 ( T ) 0 t 0 1 u x 2 d x d s C 1 ( T ) + C 1 ( T ) 0 t 0 1 | ( ρ 1 + θ u x - ρ γ ) x d s 0 1 u x 2 d x d s C 1 ( T ) + C 1 ( T ) 0 t 0 1 | u t —d s 0 1 u x 2 d x d s + C 1 ( T ) 0 t 0 1 | f | d s 0 1 u x 2 d x d s C 1 ( T ) + 1 4 0 t 0 1 u t 2 d x d s + C 1 ( T ) 0 t 0 1 f 2 d x d s + C 1 ( T ) 0 t 0 1 u x 2 d x 2 d s C 1 ( T ) + 1 4 0 t 0 1 u t 2 d x d s + C 1 ( T ) 0 t 0 1 u x 2 d x 2 d s .
(2.22)

Inserting the above estimate into (2.21),

0 1 u x 2 d x + 0 t 0 1 u t 2 d x d s C 1 ( T ) + C 0 t u x 2 0 1 u x 2 d x d s .

which, by virtue of Gronwall's inequality, (2.1) and (2.14), gives (2.19).

Proof of Theorem 2.1 By Lemmas 2.1-2.3, we complete the proof of Theorem 2.1.

3 Global existence of solutions in H2

For external force f(r, t), we suppose

f ( r , t ) L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , f r ( r , t ) L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , f t ( r , t ) L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] )
(3.1)

Constant C2(T) denotes generic positive constant depending only on the H2-norm of initial data ( ρ 0 , u 0 ) , f L ( [ 0 , T ] ) , L 2 [ 0 , 1 ] ) , f r L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , f t L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , time T and constant C1(T).

Remark 3.1 By (3.1), we easily know that assumptions (3.1) is equivalent to the following conditions

f ( r ( x , t ) , t ) L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) ,
(3.2)
f r ( r ( x , t ) , t ) L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , f t ( r ( x , t ) , t ) L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) .
(3.3)

Therefore, the generic constant C2(T) depending only on the norm of initial data (ρ0,u0) in H2, the norms of f in the class of functions in (3.2)-(3.3) and time T.

Theorem 3.1 Let 0 < θ < 1, γ > 1, and assume that the initial data satisfies (ρ0,u0) H2and external force f satisfies conditions (3.1), then there exists a unique global solution (ρ (x,t),u(x,t)) to problem (1.8)-(1.11), such that for any T > 0,

ρ L ( [ 0 , T ] , H 2 [ 0 , 1 ] ) , u L ( [ 0 , T ] , H 2 [ 0 , 1 ] L 2 ( [ 0 , T ] , H 3 [ 0 , 1 ] ) ,
(3.4)
u t L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 1 [ 0 , 1 ] ) .
(3.5)

The proof of Theorem 3.1 can be divided into the following several lemmas.

Lemma 3.2 Under the assumptions in Theorem 3.1, for any 0 ≤ tT, we have the following estimates

u t ( t ) 2 + 0 t 0 1 u t x 2 ( x , s ) d x d s C 2 ( T ) ,
(3.6)
u x ( t ) L 2 + u x x ( t ) 2 d x C 2 ( T ) .
(3.7)

Proof Differentiating (1.9) with respect to t, multiplying the resulting equation by u t in L2[0,1], performing an integration by parts, and using Lemma 2.1, we have

1 2 d d t u t 2 + 0 1 ρ 1 + θ u t x 2 d x = 0 1 ( θ + 1 ) ρ θ + 2 u x 2 - γ ρ γ + 1 u x + f t u t x d x 1 2 0 1 ρ 1 + θ u t x 2 d x + C 1 ( T ) 0 1 ρ 2 θ + 3 u x 4 + ρ 2 γ + 1 - θ u x 2 d x + C 1 ( T ) 0 1 ( f r r t ) 2 + f t 2 d x .
(3.8)

Integrating (3.8) with respect to t, applying the interpolation inequality, we conclude

u t ( t ) 2 + 0 t 0 1 ρ 1 + θ u t x 2 d x d s u t ( x , 0 ) + C 1 ( T ) 0 t 0 1 u x 4 + u x 2 + f r 2 u 2 + f t 2 d x d s u t ( x , 0 ) + C 1 ( T ) 0 t u x 2 + u x x 1 4 u x 3 4 + u x 4 ( s ) d s + 0 t u L 2 0 1 f r 2 d x d s + C 1 ( t ) 0 t 0 1 f t 2 d x d s .
(3.9)

On the other hand, by (1.9), we get

u 0 t = - γ ρ 0 γ - 1 ρ 0 x + ρ 0 θ + 1 u 0 x x + ( θ + 1 ) ρ 0 θ ρ 0 x u 0 x + f ( r 0 , 0 ) .
(3.10)

We derive from assumption (3.1) and (3.10) that

0 1 u 0 t 2 ( x ) d x C 2 ( T ) .
(3.11)

Inserting (3.11) into (3.9), by virtue of Lemmas 2.1-2.3 and assumption (3.1), we obtain (3.6). We infer from (1.9),

u t = - γ ρ γ - 1 ρ x + ρ θ + 1 u x x + ( θ + 1 ) ρ θ ρ x u x + f ( r , t ) .
(3.12)

Multiplying (3.12) by u xx in L2[0,1], we deduce

0 1 ρ θ + 1 u x x 2 d x = 0 1 u x x u t + γ ρ γ - 1 ρ x - ( θ + 1 ) ρ θ ρ x u x - f ( r , t ) d x .
(3.13)

Using Young's inequality and Sobolev's embedding theorem W1,1 W, Lemma 2.1 and (3.6), we deduce from (3.13) that

0 1 u x x 2 d x C 1 ( T ) 0 1 u t 2 + ρ x 2 + ρ x 2 u x 2 + f 2 d x + 1 4 0 1 u x x 2 d x C 2 ( T ) + C 1 ( T ) u x L 2 0 1 ρ x 2 d x + 1 4 0 1 u x x 2 d x C 2 ( T ) + 1 2 0 1 u x x 2 d x

whence

0 1 u x x 2 d x C 2 ( T ) .
(3.14)

Applying embedding theorem, we derive from (3.14) that

u x L 2 C 1 ( T ) u x 2 + u x x 2 C 2 ( T )

which, along with (3.14), gives (3.7). The proof is complete.

Lemma 3.3 Under the assumptions in Theorem 3.1, for any 0 ≤ tT, we have the following estimates

ρ x x ( t ) 2 + 0 t ρ x x ( s ) 2 d s C 2 ( T ) ,
(3.15)
0 t u x x x ( s ) 2 d x C 2 ( T ) .
(3.16)

Proof Differentiating (1.9) with respect to x, exploiting (1.8), we have

u t x = - ρ γ + ρ 1 + θ u x x x + d f d x = - γ ( γ - 1 ) ρ γ - 2 ρ x 2 - γ ρ γ - 1 ρ x x + ( θ + 1 ) θ ρ θ - 1 ρ x 2 u x + ( θ + 1 ) ρ θ ρ x x u x + 2 ( θ + 1 ) ρ θ ρ x u x x + ρ θ + 1 u x x x + 1 ρ f r
(3.17)

which gives

ρ θ - 1 ρ x x t + P ρ ρ x x = E ( x , t ) ,
(3.18)

with

E ( x , t ) = - P ρ ρ ρ x 2 - 2 ( 1 - θ ) ρ θ ρ x u x x + ( 1 + θ ) θ ρ θ - 1 ρ x 2 u x - 2 ρ θ - 1 ρ x 2 u x - u t x + 1 ρ f r .

Multiplying (3.18) by ρθ-1ρ xx , integrating the resultant over [0,1], using condition (3.1), Young's inequality, Lemma 3.2 and Theorem 2.1, we deduce

d d t ρ θ - 1 ρ x x 2 + 0 1 γ ρ γ + θ - 2 ρ x x 2 d x C 1 ( T ) 0 1 ρ x 4 + u t x 2 + ρ x 4 u x 2 + ρ x 2 u x x 2 + f x 2 d x .
(3.19)

Integrating (3.19) with respect to t over [0,t], using Theorem 2.1, Lemma 3.2 and the interpolation inequality, we derive

ρ x x ( t ) 2 + 0 t ρ x x ( s ) 2 d s C 2 ( T ) + C 1 ( T ) 0 t u x L 2 0 1 ρ x 2 d x d s + C 1 ( T ) 0 t 0 1 ρ x 4 + u t x 2 d x d s + C 1 ( T ) 0 t ρ x L 2 0 1 u x x 2 d x d s + C 1 ( T ) 0 t 0 1 f r 2 d x d s C 2 ( T ) + C 1 ( T ) 0 t 0 1 ρ x 2 d x d s + 1 2 0 t ρ x x ( s ) 2 d s
(3.20)

which, along with Lemma 2.1, gives estimate (3.15).

Differentiating (1.9) with respect to x, we can obtain

u x x x = ρ - 1 - θ u t x + γ ( γ - 1 ) ρ γ - 2 ρ x 2 + γ ρ γ - 1 ρ x x - ( θ + 1 ) ρ θ ρ x x u x + 2 ( θ + 1 ) ρ θ ρ x u x x + ( θ + 1 ) θ ρ θ - 1 ρ x 2 u x - f x .
(3.21)

Integrating (3.21) with respect to x and t over [0,1] × [0,t], applying the embedding theorem, Lemmas 2.1-2.3 and Lemma 3.1, and the estimate (3.15), we conclude

0 t 0 1 u x x x 2 d x d s C 1 ( T ) 0 t 0 1 u t x 2 + ρ x 4 + ρ x x 2 + ρ x 2 u x x 2 + ρ x 4 u x 2 + ρ x x 2 u x 2 + f r 2 d x d s C 1 ( T ) 0 t u x L 2 0 1 ρ x x 2 + ρ x 4 d x d s + C 1 ( T ) 0 t ρ x L 2 u x x 2 d s + C 1 ( T ) 0 t 0 1 ρ x 4 + u t x 2 + ρ x x 2 + f r 2 d x d s C 2 ( T ) .
(3.22)

The proof is complete.

Proof of Theorem 3.1 By Lemmas 3.2-3.3, Theorem 2.1 and Sobolev's embedding theorem, we complete the proof of Theorem 3.1.

4 Global existence of solutions in H4

For external force f(r,t), besides (3.1), we assume that

f r , f t , f r t L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) , f r r , f r t , f t t , f r r r L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) .
(4.1)

Remark 4.1 By (4.1), we easily know that assumptions (4.1) is equivalent to the following conditions

f r ( r ( x , t ) , t ) , f t ( r ( x , t ) , t ) , f r r ( r ( x , t ) , t ) L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) ,
(4.2)
f r r ( r ( x , t ) , t ) , f r t ( r ( x , t ) , t ) , f t t ( r ( x , t ) , t ) , f r r r ( r ( x , t ) , t ) L 2 ( [ 0 , T ] , L 2 [ 0 , 1 ] ) .
(4.3)

Therefore, the generic constant C4(T) depending only on the norm of initial data (ρ0,u0) in H4, the norms of f in the class of functions in (4.2)-(4.3) and time T.

Theorem 4.1 Let 0 < θ < 1, γ > 1, and assume that the initial data satisfies (ρ0,u0) H4and external force f satisfies conditions (4.1), then there exists a unique global solution (ρ (x,t),u(x,t)) to problem (1.8)-(1.11), such that for any T > 0,

ρ L ( [ 0 , T ] , H 4 [ 0 , 1 ] ) , ρ t L ( [ 0 , T ] , H 3 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 4 [ 0 , 1 ] ) ,
(4.4)
ρ t t L ( [ 0 , T ] , H 1 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 2 [ 0 , 1 ] ) ,
(4.5)
u L ( [ 0 , T ] , H 4 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 5 [ 0 , 1 ] ) ,
(4.6)
u t L ( [ 0 , T ] , H 2 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 3 [ 0 , 1 ] ) ,
(4.7)
u t t L ( [ 0 , T ] , L 2 [ 0 , 1 ] ) L 2 ( [ 0 , T ] , H 1 [ 0 , 1 ] ) .
(4.8)

The proof of Theorem 4.1 can be divided into the following several lemmas.

Lemma 4.2 Under the assumptions of Theorem 4.1, the following estimates hold for any t [0,T],

u t x ( x , 0 ) + u t x x ( x , 0 ) + u t t ( x , 0 ) C 4 ( T ) ,
(4.9)
u t t ( t ) 2 + 0 t u t t x ( s ) 2 d s C 4 ( T ) .
(4.10)

Proof We easily infer from (1.9) and Theorem 2.1, Theorem 3.1 that

u t ( t ) C 2 ( T ) ( u x ( t ) H 1 + ρ x ( t ) + f ( t ) ) .
(4.11)

Differentiating (1.9) with respect to x and exploiting Lemmas 2.1-2.3, we have

u t x ( t ) C 2 ( T ) ( u x ( t ) H 2 + ρ x ( t ) H 1 + f r ( t ) ) ,
(4.12)

or

u x x x ( t ) C 2 ( T ) ( u x ( t ) H 1 + ρ x ( t ) H 1 + u t x ( t ) + f r ( t ) ) .
(4.13)

Differentiating (1.9) with respect to x twice, using Lemmas 2.1-2.3, 3.2-3.3 and the embedding theorem, we have

u t x x ( t ) C 2 ( T ) ( u x ( t ) H 3 + ρ x ( t ) H 2 + f r ( t ) + f r r ( t ) ) ,
(4.14)

or

u x x x x ( t ) C 2 ( T ) ( u x ( t ) H 2 + ρ x ( t ) H 2 + u t x x ( t ) + f r ( t ) + f r r ( t ) ) .
(4.15)

Differentiating (1.9) with respect to t, and using Lemmas 2.1-2.3 and (1.8), we deduce that

u t t ( t ) C 2 ( T ) ( u t x ( t ) + u x ( t ) H 1 + ρ x ( t ) + u t x x ( t ) + f r ( t ) + f t ( t ) )
(4.16)

which together with (4.12) and (4.14) implies

u t t ( t ) C 2 ( T ) ( u x ( t ) H 3 + ρ x ( t ) H 2 + f r ( t ) + f t ( t ) + f r r ( t ) ) .
(4.17)

Thus, estimate (4.9) follows from (4.12), (4.14), (4.17) and condition (4.1).

Now differentiating (1.9) with respect to t twice, multiplying the resulting equation by u tt in L2([0,1]), and using integration by parts, (1.8) and the boundary condition (1.10), we deduce

0 1 u t t t u t t d x = 0 1 ( - ρ γ + ρ 1 + θ u x ) t t x + d 2 f d t 2 u t t d x = - 0 1 ( - ρ γ + ρ 1 + θ u x ) t t u t t x d x + 0 1 f r r r t 2 + f r r t t + f r t + f r t r t + f t t u t t d x - 0 1 ρ 1 + θ u t t x 2 d x + 1 2 0 1 ρ 1 + θ u t t x 2 d x + 1 2 0 1 u t t 2 d x + C 1 ( T ) 0 1 u x 4 + u t x 2 + u x 2 u t x 2 + u x 6 + f r 2 u t 2 + f r r 2 + f t 2 + f t t 2 d x
(4.18)

here, we use d 2 f d t 2 = f r r r t 2 + f r t r t + f r r 1 2 + f t t . Integrating (4.18) with respect to t, applying assumption (4.1) and (4.9), we have

u t t ( t ) 2 + 0 t 0 1 ρ 1 + θ u t t x 2 d x d s C 4 ( T ) + 1 2 0 t u t t ( s ) 2 d s + C 1 ( T ) 0 t u x 2 + u t x 2 + u x L 6 6 + u x L 2 u t x 2 + f r 2 u t L 2 ( s ) d s C 4 ( T ) + 1 2 0 t u t t ( s ) 2 d s + C 2 ( T ) 0 t u x H 1 2 + u t H 1 2 ( s ) d s

which, with Lemmas 2.1-2.3 and Theorem 3.1, implies

u t t ( t ) 2 + 0 t 0 1 ρ 1 + θ u t t x 2 d x d s C 4 ( T ) + 1 2 0 t u t t ( s ) 2 d s , t [ 0 , T ] .
(4.19)

If we apply Gronwall's inequality to (4.19), we conclude (4.11). The proof is complete.

Lemma 4.3 Under the assumptions of Theorem 4.1, the following estimate holds for any t [0,T],

u t x ( t ) 2 + 0 t u t x x ( s ) 2 d s C 4 ( T ) .
(4.20)

Proof Differentiating (1.9) with respect to x and t, multiplying the resulting equation by u tx in L2[0,1], and integrating by parts, we deduce that

0 1 u t t x u t x d x = 0 1 ( - ρ γ + ρ 1 + θ u x ) t x x + 2 f t x u t x d x = ( - ρ γ + ρ 1 + θ u x ) t x u t x | 0 1 - 0 1 ( - ρ γ + ρ 1 + θ u x ) t x u t x x d x + 0 1 ( f r r r t r x + f r r t x + f r t r x ) u t x d x = B 1 + B 2 + B 3
(4.21)

where

B 1 = ( - ρ γ + ρ 1 + θ u x ) t x u t x 0 , 1 B 2 = - 0 1 ( - ρ γ + ρ 1 + θ u x ) t x u t x x d x , B 3 = 0 1 ( f r r r t r x + f r r t x + f r t r x ) u t x d x .

Employing Theorem 2.1, Theorem 3.1 Lemma 4.2 and the interpolation inequality, we conclude

B 1 C 2 ( T ) ( u x x L + ρ x L u x L + ρ x L u t x L + u t x x L + u x L u x x L + ρ x L u x L 2 ) u t x L C 2 ( T ) ( B 01 + B 02 ) u t x 1 2 u t x x 1 2
(4.22)

with

B 0 1 = | | u x | | H 2 + | | ρ x | | H 1 , B 0 2 = | | u t x | | 1 2 | | u t x x | | 1 2 + | | u t x x | | 1 2 | | u t x x x | | 1 2 .

Applying Young's inequality several times, we have that for any ε (0,1),

C 2 ( T ) B 0 1 | | u t x | | 1 2 | | u t x x | | 1 2 ε 2 2 | | u t x x | | 2 + C 2 ( T ) ε - 3 ( | | u t x | | 2 + | | u x | | H 2 2 = + | | ρ x | | H 1 2 ) ,
(4.23)

and

C 2 ( T ) B 0 2 | | u t x | | 1 2 | | u t x x | | 1 2 ε 2 2 | | u t x x | | 2 + ε 2 | | u t x x x | | 2 + | | C 2 ( T ) ε - 6 | | u t x | | 2 .
(4.24)

Thus we infer from (4.22)-(4.24) that

B 1 ε 2 | | u t x x | | 2 + | | u t x x x | | 2 + C 2 ( T ) ε - 6 | | u t x | | 2 + | | u x | | H 2 2 + | | ρ x | | H 1 2
(4.25)

which, together with Theorem 2.1, Theorem 3.1 and Lemma 4.2, implies

0 t B 1 ( s ) d s C 2 ( T ) + ε 2 0 t | | u t x x | | 2 + | | u t x x x | | 2 ( s ) d s .
(4.26)

On the other hand, differentiating (1.9) with respect to x and t, and using Theorem 3.1 and Lemma 4.2, we derive

| | u t x x x ( t ) | | 2 C 2 ( T ) | | u x | | H 2 2 + | | ρ x | | H 1 2 + | | u t x | | H 1 2 + | | u t t x | | 2 + | | 2 f x t | | 2 C 2 ( T ) | | u x | | H 2 2 + | | ρ x | | H 1 2 + | | u t x | | H 1 2 + | | u t t x | | 2 + | | f r r | | 2 + | | f r | | 2 | | u x | | L 2 + | | f r t | | 2
(4.27)

Inserting (4.27) into (4.26), employing Theorem 2.1, Theorem 3.1 and Lemma 4.2, we conclude

0 t B 1 ( s ) d s C 4 ( T ) + ε 2 0 t | | u t x x ( s ) | | 2 d s .
(4.28)

Similarly, by Theorem 2.1, Theorem 3.1, Lemma 4.2 and the embedding theorem, we get that for any ε (0,1),

B 2 - 0 1 ρ 1 + θ u t x x 2 d x + ε 0 1 ρ 1 + θ u t x x 2 d x + C 2 ( T ) 0 1 ρ x 2 u x 4 + ρ x 2 u t x 2 + ρ x 2 u x 2 + u x x 2 d x .
(4.29)

By virtue of assumption (4.1), Theorem 2.1 and Theorem 3.1, we derive that

B 3 = C 1 ( T ) ( u L f r r + u x L f r + f r t ) C 2 ( T ) ( f r r + f r + f r t )

which, combined with (4.21) and (4.27)-(4.29), gives

d dt | | u t x ( t ) | | 2 + 0 1 ρ 1 + θ u t x x 2 d x ε 2 ( | | u t x x | | 2 + | | u t x x x | | 2 ) + C 2 ( T ) | | u t x | | 2 + | | u x | | H 2 2 + | | ρ x | | H 1 2 + | | f r r | | 2 + | | f r | | 2 + | | f r t | | 2 .
(4.30)

Integrating (4.30) with respect to t, picking ε small enough, using Theorem 2.1 and Theorem 3.1, Lemma 4.2 and assumption (4.1), we complete the proof of estimate (4.20).

Lemma 4.4 Under the assumptions of Theorem 4.1, the following estimates hold for any t [0,T],

| | ρ x x x ( t ) | | 2 + | | ρ x x x x ( t ) | | C 4 ( T ) ,
(4.31)
| | u x x x ( t ) | | H 1 2 + | | u t x x ( t ) | | 2 + 0 t | | u t t | | H 1 2 + | | u t x x | | H 1 2 ( s ) d s C 4 ( T ) ,
(4.32)
0 t | | u x x x x ( s ) | | H 1 2 d s C 4 ( T ) .
(4.33)

Proof Differentiating (3.18) with respect to x, we have

( ρ θ - 1 ρ x x x ) t + P ρ ρ x x x = E 1 ( x , t )
(4.34)

where

E 1 ( x , t ) = E x ( x , t ) - P ρ ρ ρ x ρ x x - ( θ - 1 ) ( ρ θ - 2 ρ x ρ x x ) t .
(4.35)

An easy calculation with the interpolation inequality, Theorem 2.1 and Theorem 3.1, gives

| | E x ( t ) | | C 2 ( T ) ( | | ρ x ( t ) | | L 6 3 + | | ρ x ρ x x | | + | | ρ x u x x x | | + | | ρ x x u x x | | + | | ρ x | | L 3 | | u x | | + | | ρ x | | L 2 | | u x x | | + | | u t x x | | + + | | ρ x | | L | | f r | + | | f r r | | ) C 2 ( T ) ( | | ρ x ( t ) | | H 1 + | | u x ( t ) | | H 2 + | | u t x x | | + | | f r | | + | | f r r | | ) ,
(4.36)

and

| | E 1 | | C 2 ( T ) ( | | ρ x ( t ) | | H 1 + | | u x ( t ) | | H 2 + | | u t x x ( t ) | | + | | f r | | + | | f r r | | ) .
(4.37)

By virtue of Theorem 2.1 and Theorem 3.1, we infer from (4.36)-(4.37), (4.20) and assumption (4.1) that

0 t | | E 1 ( s ) | | 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.38)

Now multiplying (4.34) by ρθ-1ρ xxx in L2[0,1], we obtain

d d t | | ρ θ - 1 ρ x x x | | 2 + | | ρ x x x ( t ) | | 2 C 1 ( T ) | | E 1 ( T ) | | 2 .
(4.39)

Integrating (4.39) with respect to t, using Theorem 2.1 and Theorem 3.1, assumption (4.1) and (4.38), we can get

| | ρ x x x ( t ) | | 2 + 0 t | | ρ x x x ( s ) | | 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.40)

By virtue of Theorem 2.1 and Theorem 3.1, we infer from (4.10), (4.15) and (4.40) that

| | u x x x ( t ) | | 2 + 0 t | | u x x x ( s ) | | H 1 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.41)

Differentiating (1.9) with respect to t, using Theorem 2.1 and Theorem 3.1 and Lemmas 4.2-4.3, we infer that for any t [0, T],

| | u t x x ( t ) | | C 2 ( T ) | | u t t ( t ) | | + C 2 ( T ) ( | | u x ( t ) | | H 1 + | | u t x ( t ) | | + | | ρ x ( t ) | | ) C 4 ( T )
(4.42)

which, combined with (4.15), (4.40) and (4.42), gives

| | u x x x x ( t ) | | 2 + 0 t | | u t x x ( s ) | | 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.43)

Differentiating (4.34) with respect to x, we see that

ρ θ - 1 ρ x x x x t + P ρ ρ x x x x = E 2 ( x , t ) ,
(4.44)

with

E 2 ( x , t ) = E 1 x ( x , t ) - P ρ ρ ρ x ρ x x x - ( θ - 1 ) ρ θ - 2 ρ x ρ x x x t

and

E 1 x ( x , t ) = E x x ( x , t ) P ρ ρ ρ x ρ x x ) x ( θ 1 ) ( ρ θ 2 ρ x ρ x x ) t x .

Using the embedding theorem, (1.8), Theorem 2.1, Theorem 3.1 and Lemmas 4.1-4.2, we can deduce that

| | E x x ( t ) | | C 4 ( T ) ( | | u x ( t ) | | H 3 + | | ρ x ( t ) | | H 2 + | | f r | | + | | f r r | | + | | f r r r | | ) ,
(4.45)
| | E 1 x ( t ) | | C 4 ( T ) ( | | u x ( t ) | | H 3 + | | ρ x ( t ) | | H 2 + | | u t x ( t ) | | H 2 + | | f r | | + | | f r r | | + | | f r r r | | ) ,
(4.46)
| | E 2 ( t ) | | C 4 ( T ) ( | | u x ( t ) | | H 3 + | | ρ x ( t ) | | H 2 + | | u t x ( t ) | | H 2 + | | f r | | + | | f r r | | + | | f r r r | | ) .
(4.47)

Inserting (4.46) into (4.47), we have

| | E 2 ( t ) | | C 4 ( T ) ( | | u x ( t ) | | H 3 + | | ρ x ( t ) | | H 2 + | | u t x ( t ) | | H 1 ) + | | u t t x ( t ) | | + | | f r ( t ) | | + | | f r r ( t ) | | + | | f r r r ( t ) | | ) .
(4.48)

By virtue of Theorems 2.1, 3.1, Lemmas 4.2-4.3, we derive from (4.40)-(4.43) and assumption (4.1) that

0 t | | E 2 ( s ) | | 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.49)

Multiplying (4.44) by ρθ-1ρ xxxx in L2[0,1], we get

d dt | | ρ θ - 1 ρ x x x x | | 2 + | | ρ x x x x ( t ) | | 2 C 1 ( T ) | | E 2 ( t ) | | 2 .
(4.50)

Integrating (4.50) with respect to t, using condition (4.1) and (4.49), we conclude

| | ρ x x x x ( t ) | | 2 + 0 t | | ρ x x x x ( s ) | | 2 d s C 4 ( T ) , t [ 0 , T ] .
(4.51)

Differentiating (1.9) with respect to x three times, using Theorems 2.1, 3.1, Lemmas 4.2-4.3 and the interpolation inequality, we infer

| | u x x x x x ( t ) | | C 4 ( T ) ( | | u t x x x ( t ) | | + | | u x ( t ) | | H 3 + | | ρ x ( t ) | | H 3 + | | f r ( t ) | | + | | f r r ( t ) | | + | | f r r r ( t ) | | ) .
(4.52)

Thus we conclude from (1.8), (4.27), (4.41), (4.43), (4.51) and assumption (4.1) that

0 t ( | | u x x x x x | | 2 + | | u t x x x | | 2 ) ( s ) d s C 4 ( T ) , t [ 0 , T ] .
(4.53)

Thus (4.31) follows from (4.40) and (4.51), we can derive estimate (4.32)-(4.33) from Theorem 2.1, Theorem 3.1, Lemmas 4.2-4.3, (4.41), (4,43) and (4.53). The proof is complete.

Proof of Theorem 4.1 Using (1.8),Theorem 2.1, 3.1 and Lemmas 4.2-4.4 and the proper interpolation inequality, we readily get estimate (4.4)-(4.8) and complete the proof from Theorem 4.1.

Corollary 4.5 Under assumptions of Theorem 4.1 and some suitable compatibility conditions, the global solution (ρ (x,t),u(x,t)) to problem (1.8)-(1.11) is the classical solution verifying

| | ρ ( t ) | | C 3 + 1 2 + | | u ( t ) | | C 3 + 1 2 C 4 ( T ) .
(116)

Proof By the embedding theorem, we easily prove the corollary from Theorem 4.1.

References

  1. Hoff D, Serre D: The failure of continuous dependence on initial data for the Navier- Stokes equations of compressible flow. SIMA J Appl Math 1991, 51: 887-898. 10.1137/0151043

    Article  MathSciNet  Google Scholar 

  2. Fang D, Zhang T: Discontinuous solutions of the compressible Navier-Stokes equations with degenerate viscosity coefficient and vacuum. J Math Anal Appl 2006, 318: 224-245. 10.1016/j.jmaa.2005.05.050

    Article  MathSciNet  Google Scholar 

  3. Jiang S, Xin Z, Zhang P: Global weak solutions to 1D compressible isentropic Navier-stokes equations with density-dependent viscosity. Meth Appl Anal 2005, 12: 239-252.

    MathSciNet  Google Scholar 

  4. Liu T, Xin Z, Yang T: Vacuum states of compressible flow. Discret Cont Dyn Syst 1998, 4: 1-32.

    Article  MathSciNet  Google Scholar 

  5. Okada M, Matušů-Nečasová Š, Makino T: Free boundary problem for the equation of one-dimensional motion of compressible gas with density-dependent viscosity. Ann Univ Ferrara Sez VII (N.S.) 2002, 48: 1-20.

    MathSciNet  Google Scholar 

  6. Xin Z, Yao Z: The existence, uniqueness and regularity for one-dimensional compressible Navier-Stokes equations. (Preprint)

  7. Yang T, Yao Z, Zhu C: Compressible Navier-Stokes equations with degenerate viscosity coefficient and vacuum. Commun Partial Differ Equ 2001, 26: 965-981. 10.1081/PDE-100002385

    Article  MathSciNet  Google Scholar 

  8. Qin Y, Huang L, Yao Z: Regularity of 1D compressible isentropic Navier-Stokes equations with density-dependent viscosity. J Differ Equ 2008, 245: 3956-3973. 10.1016/j.jde.2008.03.014

    Article  MathSciNet  Google Scholar 

  9. Qin Y, Huang L, Yao Z: A remark on regularity of 1D compressible isentropic Navier- Stokes equations with density-dependent viscosity. J Math Anal Appl 2009, 351: 497-508. 10.1016/j.jmaa.2008.10.044

    Article  MathSciNet  Google Scholar 

  10. Yang T, Zhao H: A vacuum problem for the one-dimensional compressible Navier-Stokes equations with density-dependent viscosity. J Differ Equ 2002, 184: 163-184. 10.1006/jdeq.2001.4140

    Article  Google Scholar 

  11. Yang T, Zhu C: Compressible Navier-Stokes equations with degenerate viscosity coefficient and vacuum. Commun Math Phys 2002, 230: 329-363. 10.1007/s00220-002-0703-6

    Article  MathSciNet  Google Scholar 

  12. Vong S, Yang T, Zhu C: Compressible Navier-Stokes equations with degenerate viscosity coefficient and vacuum(II). J Differ Equ 2003, 192: 475-501. 10.1016/S0022-0396(03)00060-3

    Article  MathSciNet  Google Scholar 

  13. Fang D, Zhang T: Compressible Navier-Stokes equations with vacuum state in one dimension. Commun Pure Appl Anal 2004, 3: 675-694.

    Article  MathSciNet  Google Scholar 

  14. Fang D, Zhang T: A note on compressible Navier-Stokes equations with vacuum state in one dimension. Nonlinear Anal 2004, 58: 719-731. 10.1016/j.na.2004.05.016

    Article  MathSciNet  Google Scholar 

  15. Qin Y, Huang L, Deng S, Ma Z, Su X, Yang X: Interior regularity of the compressible Navier-Stokes equations with degenerate viscosity coefficient and vacuum. Discret Cont Dyn Syst Ser S 2009, 2: 163-192.

    Article  MathSciNet  Google Scholar 

  16. Qin Y, Zhao Y: Global existence and asymptotic behavior of the compressible Navier- Stokes equations for a 1D isothermal viscous gas. Math Models Methods Appl Sci 2008, 18: 1383-1408. 10.1142/S0218202508003078

    Article  MathSciNet  Google Scholar 

  17. Zhang T, Fang D: Global behavior of compressible Navier-Stokes equations with a degenerate viscosity coefficient. Arch Ration Mech Anal 2006, 182: 223-253. 10.1007/s00205-006-0425-6

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The work is in part supported by Doctoral Foundation of North China University of Water Sources and Electric Power (No. 201087), the Natural Science Foundation of Henan Province of China (No. 112300410040) and the NNSF of China (No. 11101145).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lan Huang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors contributed to each part of this work equally.

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

Huang, L., Lian, R. Global behavior of 1D compressible isentropic Navier-Stokes equations with a non-autonomous external force. Bound Value Probl 2011, 43 (2011). https://doi.org/10.1186/1687-2770-2011-43

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2011-43

Keywords