Skip to main content

On the free boundary value problem for one-dimensional compressible Navier-Stokes equations with constant exterior pressure

Abstract

In this paper, we consider the free boundary value problem (FBVP) for one-dimensional isentropic compressible Navier-Stokes equations (CNS) with density-dependent viscosity coefficient and constant exterior pressure. Under certain assumptions imposed on the initial data, the global existence and uniqueness of a strong solution to FBVP for CNS are established, in particular, the strong solution tends pointwise to a non-vacuum equilibrium state at an exponential time-rate as the time tends to infinity.

1 Introduction

In the present paper, we consider the free boundary value problem to one-dimensional isentropic compressible Navier-Stokes equations with density-dependent viscosity coefficient for regular initial data in the case that across the free surface the stress tensor is balanced by constant exterior pressure. In general, one-dimensional isentropic compressible Navier-Stokes equations with density-dependent viscosity coefficient can be written as

{ ρ t + ( ρ u ) x = 0 , ( ρ u ) t + ( ρ u 2 + P ( ρ ) ) x = ( μ ( ρ ) u x ) x , ( x , t ) R × [ 0 , T ] ,
(1.1)

where ρ>0, u and P(ρ)= ρ γ (γ>1) stand for the flow density, velocity and pressure, respectively, and the viscosity coefficient is μ(ρ)= ρ α with α>0. Note here that the case γ=2 and α=1 in (1.1) corresponds to the viscous Saint-Venant system for shallow water.

Recently, there have been much significant progress achieved on the compressible Navier-Stokes equations with density-dependent viscosity coefficients. For instance, the mathematical derivations are obtained in the simulation of flow surface in shallow region [1, 2]. The existence of solutions for the 2D shallow water equations is investigated by Bresch and Desjardins [3, 4]. The well-posedness of solutions to the free boundary value problem with initial finite mass and the flow density being connected with the infinite vacuum either continuously or via jump discontinuity is considered by many authors; refer to [515] and the references therein. The global existence of classical solutions is shown by Mellet and Vasseur [16]. The qualitative behaviors of global solutions and dynamical asymptotics of vacuum states are also made, such as the finite time vanishing of finite vacuum or the asymptotical formation of vacuum in large time, the dynamical behaviors of vacuum boundary, the large time convergence to rarefaction wave with vacuum, and the stability of shock profile with large shock strength; refer to [1722] and the references therein.

In addition, some important progress has been made about free boundary value problems for multi-dimensional compressible viscous Navier-Stokes equations with constant viscosity coefficients for either barotropic or heat-conducive fluids by many authors; for example, in the case that across the free surface stress tensor is balanced by a constant exterior pressure and/or the surface tension, classical solutions with strictly positive densities in the fluid regions to FBVP for CNS (1.1) with constant viscosity coefficients are proved locally in time for either barotropic flows [2325] or heat-conductive flows [2628]. In the case that across the free surface the stress tensor is balanced by exterior pressure [25], surface tension [29], or both exterior pressure and surface tension [30], respectively, as the initial data is assumed to be near to a non-vacuum equilibrium state, the global existence of classical solutions with small amplitude and positive densities in fluid region to the FBVP for CNS (1.1) with constant viscosity coefficients is obtained. The global existence of classical solutions to FBVP for compressible viscous and heat-conductive fluids is also established with the stress tensor balanced by the exterior pressure and/or surface tension across the free surface; refer to [31, 32] and the references therein.

In this paper, we consider the free boundary value problem (FBVP) for one-dimensional isentropic compressible Navier-Stokes equations with density-dependent viscosity coefficient and constant exterior pressure, and we focus on the existence, regularities and dynamical behaviors of a global strong solution, etc. As γ>1 we show that the free boundary value problem with regular initial data admits a unique global strong solution which tends pointwise to a non-vacuum equilibrium state at an exponential time-rate as the time tends to infinity (refer to Theorem 2.1 for details).

The rest of the paper is arranged as follows. In Section 2, the main results about the existence and dynamical behaviors of a global strong solution to FBVP with two different initial data for compressible Navier-Stokes equations are stated. Then, some important a priori estimates are given in Section 3 and the theorem is proven in Section 4.

2 Main results

We are interested in the global existence and dynamics of the free boundary value problem for (1.1) with the following initial data and boundary conditions:

{ ( ρ , u ) ( x , 0 ) = ( ρ 0 , u 0 ) , x [ a 0 , b 0 ] , ( ρ γ ρ α u x ) ( a ( t ) , t ) = p e , ( ρ γ ρ α u x ) ( b ( t ) , t ) = p e , t > 0 ,
(2.1)

where x=a(t) and x=b(t) are the free boundaries defined by

{ d d t a ( t ) = u ( a ( t ) , t ) , a ( 0 ) = a 0 , d d t b ( t ) = u ( b ( t ) , t ) , b ( 0 ) = b 0 , t > 0 ,
(2.2)

and the positive constant p e is the exterior pressure.

Without loss of generality, the total initial mass is renormalized to be one, i.e.,

a ( t ) b ( t ) ρ(x,t)dx= a 0 b 0 ρ 0 (x)dx:=1.
(2.3)

Define

E 0 := 1 2 a 0 b 0 ρ 0 u 0 2 dx+ a 0 b 0 ρ 0 ( 1 γ 1 ( ρ 0 γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 0 1 ρ ¯ 1 ) ) dx
(2.4)

and

E 1 : = 1 2 a 0 b 0 ρ 0 ( u 0 + 1 α ( ρ 0 α ) x ) 2 d x + a 0 b 0 ρ 0 ( 1 γ 1 ( ρ 0 γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 0 1 ρ ¯ 1 ) ) d x ,
(2.5)

where ρ ¯ := p e 1 γ >0.

Firstly, we consider the initial data satisfying

{ inf [ a 0 , b 0 ] ρ 0 ρ ̲ > 0 , ( ρ 0 , u 0 ) W 1 , ( [ a 0 , b 0 ] ) , ( ρ 0 γ ρ 0 α u 0 x ) ( a 0 ) = ρ ¯ γ , ( ρ 0 γ ρ 0 α u 0 x ) ( b 0 ) = ρ ¯ γ , | ρ 0 α ( a 0 ) ρ ¯ α | δ , | ρ 0 α ( b 0 ) ρ ¯ α | δ ,
(2.6)

where ρ ̲ and δ>0 are positive constants, δ is bounded, and we also consider the other initial data satisfying

{ inf [ a 0 , b 0 ] ρ 0 ρ ̲ > 0 , ( ρ 0 , u 0 ) W 1 , ( [ a 0 , b 0 ] ) , ( ρ 0 γ ρ 0 α u 0 x ) ( a 0 ) = ρ ¯ γ , ( ρ 0 γ ρ 0 α u 0 x ) ( b 0 ) = ρ ¯ γ , ρ 0 ( b 0 ) ρ 0 ( a 0 ) ρ ¯ ,
(2.7)

note that the compatibility conditions between the initial data and boundary conditions hold. Then we have the global existence and time-asymptotical behavior of a strong solution as follows.

Theorem 2.1 (FBVP)

Let γ>1. Assume that the initial data satisfies (2.6) for 0<α 1 2 , and satisfies (2.6) together with

E 0 1 2 ( E 0 + 2 E 1 ) 1 2 ν 2 ( 2 α 1 ) ρ ¯ γ + 2 α 1 2
(2.8)

for α> 1 2 or the initial data satisfies (2.7) for 0<α 1 2 , and satisfies (2.7) together with

E 0 1 2 ( E 0 + E 1 + ( ρ 0 γ ( b 0 ) ρ ¯ γ ) b 0 ( ρ 0 γ ( a 0 ) ρ ¯ γ ) a 0 ) 1 2 ν 2 ( 2 α 1 ) ρ ¯ γ + 2 α 1 2
(2.9)

for 1 2 <αγ, where ν is a positive constant. Then there exists a unique global strong solution (ρ,u,a,b) to FBVP (1.1) and (2.1) satisfying, for T>0,

{ c ρ L ( 0 , T ; H 1 ( [ a ( t ) , b ( t ) ] ) ) C 0 ( [ a ( t ) , b ( t ) ] × [ 0 , T ] ) , u L ( 0 , T ; H 1 ( [ a ( t ) , b ( t ) ] ) ) L 2 ( 0 , T ; H 2 ( [ a ( t ) , b ( t ) ] ) ) , a ( t ) , b ( t ) H 1 ( [ 0 , T ] ) , ( ρ γ ρ α u x ) L ( 0 , T ; L 2 ( [ a ( t ) , b ( t ) ] ) ) ,
(2.10)

with c>0 being a constant independent of time.

If it further holds that u 0 H 2 ([ a 0 , b 0 ]), then (ρ,u,a,b) satisfies

{ ( ρ , u ) C 0 ( [ a ( t ) , b ( t ) ] × [ 0 , T ] ) , ρ L ( 0 , T ; H 1 ( [ 0 , a ( t ) ] ) ) , ρ t L ( 0 , T ; L 2 ( [ a ( t ) , b ( t ) ] ) ) , u L ( 0 , T ; H 2 ( [ a ( t ) , b ( t ) ] ) ) L 2 ( 0 , T ; H 3 ( [ a ( t ) , b ( t ) ] ) ) , u t L ( 0 , T ; L 2 ( [ a ( t ) , b ( t ) ] ) ) L 2 ( 0 , T ; H 1 ( [ a ( t ) , b ( t ) ] ) ) , a ( t ) , b ( t ) H 2 ( [ 0 , T ] ) , ( ρ γ ρ α u x ) C 0 ( [ 0 , T ] × ( [ a ( t ) , b ( t ) ] ) ) .
(2.11)

The solution (ρ,u) tends to the non-vacuum equilibrium state exponentially

( ρ ρ ¯ , u ) ( , t ) L ( [ a ( t ) , b ( t ) ] ) C 1 e C 2 t ,t>0,
(2.12)

where C 1 and C 2 are positive constants independent of time.

Remark 2.1 Theorem 2.1 holds for the one-dimensional Saint-Venant model for shallow water, i.e., γ=2, α=1.

Remark 2.2 The initial constraints (2.8) and (2.9) do not always require that the perturbation of the initial data around the equilibrium state ( ρ ¯ ,0) is small. Indeed, it can be large provided that the state ρ ¯ >0 is large enough.

3 The a priori estimates

It is convenient to make use of the Lagrange coordinates in order to establish the a priori estimates. Take the Lagrange coordinates transform

ξ= a ( t ) x ρ(y,t)dy,τ=t.
(3.1)

Since the conservation of total mass holds, the boundaries x=a(t) and x=b(t) are transformed into ξ=0 and ξ=1, respectively, and the domain [a(t),b(t)] is transformed into [0,1]. FBVP (1.1) and (2.1) is reformulated into

{ ρ τ + ρ 2 u ξ = 0 , u τ + ( ρ γ ) ξ = ( ρ 1 + α u ξ ) ξ , ( ρ 0 , u 0 ) = ( ρ 0 , u 0 ) ( ξ ) , ξ [ 0 , 1 ] , ( ρ γ ρ 1 + α u ξ ) ( 0 , τ ) = ρ ¯ γ , ( ρ γ ρ 1 + α u ξ ) ( 1 , τ ) = ρ ¯ γ , τ [ 0 , T ] ,
(3.2)

where the initial data satisfies

{ inf [ 0 , 1 ] ρ 0 ρ ̲ > 0 , ( ρ 0 , u 0 ) W 1 , ( [ 0 , 1 ] ) , ( ρ 0 γ ρ 0 1 + α u 0 x ) ( 0 ) = ρ ¯ γ , ( ρ 0 γ ρ 0 1 + α u 0 x ) ( 1 ) = ρ ¯ γ , | ρ 0 α ( 0 ) ρ ¯ α | δ , | ρ 0 α ( 1 ) ρ ¯ α | δ
(3.3)

or

{ inf [ 0 , 1 ] ρ 0 ρ ̲ > 0 , ( ρ 0 , u 0 ) W 1 , ( [ 0 , 1 ] ) , ( ρ 0 γ ρ 0 1 + α u 0 x ) ( 0 ) = ρ ¯ γ , ( ρ 0 γ ρ 0 1 + α u 0 x ) ( 1 ) = ρ ¯ γ , ρ 0 ( 1 ) ρ 0 ( 0 ) ρ ¯ ,
(3.4)

and the consistencies between the initial data and boundary conditions hold.

Next, we deduce the a priori estimates for the solution (ρ,u) to FBVP (3.2).

Lemma 3.1 Let T>0. Under the assumptions of Theorem  2.1, it holds for any strong solution (ρ,u) to FBVP (3.2) that

0 1 ( 1 2 u 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ + 0 τ 0 1 ρ 1 + α u ξ 2 d ξ d s = E 0 , τ [ 0 , T ] .
(3.5)

Proof Taking the product of (3.2)2 with u, integrating on [0,1] and using boundary conditions, we have

1 2 d d τ 0 1 u 2 d ξ + 0 1 ρ 1 + α u ξ 2 d ξ = 0 1 u ξ ( ρ γ ρ ¯ γ ) d ξ = 0 1 ( ρ γ 2 ρ ¯ γ ρ 2 ) ( ρ τ ) d ξ = d d τ 0 1 ( 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ ,
(3.6)

which leads to (3.5) after the integration with respect to τ[0,T]. □

Lemma 3.2 Let T>0. Under the assumptions of Theorem  2.1, it holds for any strong solution (ρ,u) to FBVP (3.2) with the initial data satisfying (3.3) that

0 1 ( 1 2 ( u + 1 α ( ρ α ) ξ ) 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ + γ 0 τ 0 1 ρ γ + 2 α 1 ρ ξ 2 d ξ d s 2 E 1 , τ [ 0 , T ] .
(3.7)

As the initial data satisfies (3.4), it holds

0 1 ( 1 2 ( u + 1 α ( ρ α ) ξ ) 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ + γ 0 τ 0 1 ρ γ + 2 α 1 ρ ξ 2 d ξ d s + ( ρ γ ( 1 , τ ) ρ ¯ γ ) b ( τ ) ( ρ γ ( 0 , τ ) ρ ¯ γ ) a ( τ ) + γ α 0 τ ( ρ γ α ( 1 , s ) ( ρ γ ( 1 , s ) ρ ¯ γ ) b ( s ) ρ γ α ( 0 , s ) ( ρ γ ( 0 , s ) ρ ¯ γ ) a ( s ) ) d s = E 1 + ( ρ 0 γ ( 1 ) ρ ¯ γ ) b 0 ( ρ 0 γ ( 0 ) ρ ¯ γ ) a 0 , τ [ 0 , T ] ,
(3.8)

where a(τ) satisfies a (τ)=u(0,τ) and a(0)= a 0 , b(τ) satisfies b (τ)=u(1,τ) and b(0)= b 0 .

Proof Multiplying (3.2)1 by ρ α 1 leads to

( ρ α ) τ /α+ ρ 1 + α u ξ =0,
(3.9)

which gives

( ρ α ) τ ξ /α+ ( ρ 1 + α u ξ ) ξ =0.
(3.10)

Summing (3.2)1 and (3.10), we deduce

( u + ( ρ α ) ξ / α ) τ + ( ρ γ ) ξ =0.
(3.11)

Multiplying (3.11) by (u+ ( ρ α ) ξ /α) and integrating the result over [0,1]×[0,τ], we have

d d τ 0 1 ( 1 2 ( u + 1 α ( ρ α ) ξ ) 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ + γ 0 1 ρ γ + 2 α 1 ρ ξ 2 d ξ + ( ( ρ γ ( 1 , τ ) ρ ¯ γ ) u ( 1 , τ ) ( ρ γ ( 0 , τ ) ρ ¯ γ ) u ( 0 , τ ) ) = 0 .
(3.12)

Next, we prove (3.7) in the case that the initial data satisfies (3.3) firstly, to obtain (3.7). We assume a priori that there are constants ρ ± >0 so that

0< ρ ρ(ξ,τ), ρ ¯ ρ + ,(ξ,τ)[0,1]×[0,T],T>0.
(3.13)

By means of (3.2)1 and the boundary condition, we have

( ρ α ( 1 , τ ) ρ ¯ α ) τ + ( ρ γ ( 1 , τ ) ρ ¯ γ ) =0,
(3.14)

which yields

ρ α (1,τ) ρ ¯ α = ( ρ 0 α ( 1 ) ρ ¯ α ) e 0 τ ρ γ ( 1 , s ) ρ ¯ γ ρ α ( 1 , s ) ρ ¯ α d s .
(3.15)

From (3.13), we can obtain

γ α ρ γ α ρ γ ( 1 , τ ) ρ ¯ γ ρ α ( 1 , τ ) ρ ¯ α γ α ρ + γ α .
(3.16)

It holds from (3.5), (3.13), (3.15) and (3.16) that

| 0 τ ( ρ γ ( 1 , s ) ρ ¯ γ ) u ( 1 , s ) d s | C ˆ 0 τ ρ + γ α | ρ α ( 1 , s ) ρ ¯ α | ( 0 1 u 2 d x + 0 1 u x 2 d x ) 1 2 d s C ˆ 0 τ | ρ 0 α ( 1 ) ρ ¯ α | e γ α ρ γ α s ( ρ + 2 ( γ α ) + 0 1 u 2 d ξ + ρ ( 1 + α ) 0 1 ρ 1 + α u x 2 d ξ ) d s C ˆ | ρ 0 α ( 1 ) ρ ¯ α | ( ρ α γ ρ + 2 ( γ α ) + ρ α γ + ρ ( 1 + α ) ) 1 2 E 1 ,
(3.17)

where C ˆ is a positive constant independent of time, and we assume that

| ρ 0 α (1) ρ ¯ α |min { ρ γ α ρ + 2 ( α γ ) 2 C ˆ E 1 , ρ γ α 2 C ˆ E 1 , ρ 1 + α 2 C ˆ E 1 } := δ 1 .
(3.18)

Using the same method we can obtain

| 0 τ ( ρ γ ( 0 , s ) ρ ¯ γ ) u(0,s)ds| 1 2 E 1 .
(3.19)

From (3.12), (3.17) and (3.19), we have (3.7).

Then, we prove (3.8) as the initial data satisfies (3.4), where we have the fact

0 τ ( ( ρ γ ( 1 , s ) ρ ¯ γ ) u ( 1 , s ) ( ρ γ ( 0 , s ) ρ ¯ γ ) u ( 0 , s ) ) d s = 0 τ ( ρ γ ( 1 , s ) ρ ¯ γ ) b ( s ) d s 0 τ ( ρ γ ( 0 , s ) ρ ¯ γ ) a ( s ) d s = ( ρ γ ( 1 , τ ) ρ ¯ γ ) b ( τ ) ( ρ 0 γ ( 1 ) ρ ¯ γ ) b 0 ( ρ γ ( 0 , τ ) ρ ¯ γ ) a ( τ ) + ( ρ 0 γ ( 0 ) ρ ¯ γ ) a 0 γ α 0 τ b ( s ) ρ γ α ( 1 , s ) ( ρ α ( 1 , s ) ) τ d s + γ α 0 τ a ( s ) ρ γ α ( 0 , s ) ( ρ α ( 0 , s ) ) τ d s = ( ρ γ ( 1 , τ ) ρ ¯ γ ) b ( τ ) ( ρ γ ( 0 , τ ) ρ ¯ γ ) a ( τ ) ( ρ 0 γ ( 1 ) ρ ¯ γ ) b 0 + ( ρ 0 γ ( 0 ) ρ ¯ γ ) a 0 + γ α 0 τ ρ γ α ( 1 , s ) ( ρ γ ( 1 , s ) ρ ¯ γ ) b ( s ) d s γ α 0 τ ρ γ α ( 0 , s ) ( ρ γ ( 0 , s ) ρ ¯ γ ) a ( s ) d s .
(3.20)

Applying equation (3.2)1 and the boundary condition, we have that

( ρ α ( 1 , τ ) ρ α ( 0 , τ ) ) τ + ( ρ γ ( 1 , τ ) ρ γ ( 0 , τ ) ) =0
(3.21)

and

ρ α (1,τ) ρ α (0,τ)= ( ρ 0 α ( 1 ) ρ 0 α ( 0 ) ) e 0 τ ρ γ ( 1 , s ) ρ γ ( 0 , s ) ρ α ( 1 , s ) ρ α ( 0 , s ) d s 0,
(3.22)

which together with

ρ α (1,τ) ρ ¯ α = ( ρ 0 α ( 1 ) ρ ¯ α ) e 0 τ ρ γ ( 1 , s ) ρ ¯ γ ρ α ( 1 , s ) ρ ¯ α d s 0,
(3.23)
ρ α (0,τ) ρ ¯ α = ( ρ 0 α ( 0 ) ρ ¯ α ) e 0 τ ρ γ ( 0 , s ) ρ ¯ γ ρ α ( 0 , s ) ρ ¯ α d s 0
(3.24)

and

b(τ)a(τ)= 0 1 1 ρ ( ζ , τ ) dζ>0
(3.25)

gives rise to (3.8). □

Lemma 3.3 Let T>0. Under the assumptions of Theorem  2.1, it holds

0< ρ ρ(ξ,τ) ρ ,(ξ,τ)[0,1]×[0,T],T>0,
(3.26)

where ρ and ρ are positive constants independent of time.

Proof Define

φ(ρ):= 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 )
(3.27)

and

ψ(ρ):= ρ ¯ ρ φ ( η ) 1 2 η α 1 dη.
(3.28)

It is easy to verify that φ(ρ)0 and ψ (ρ)0. In addition, it holds as ρ+ that

lim ρ + ψ ( ρ ) ( γ 1 ) 1 2 lim ρ + ρ ¯ ρ η γ + 2 α 3 2 d η = lim ρ + 2 ( γ + 2 α 1 ) γ 1 ( ρ γ + 2 α 1 2 ρ ¯ γ + 2 α 1 2 ) + ,
(3.29)

and as ρ0 that

lim ρ 0 ψ ( ρ ) lim ρ 0 ρ ¯ ρ ρ ¯ γ 2 η α 3 2 d η = { lim ρ 0 2 2 α 1 ρ ¯ γ 2 ( ρ α 1 2 ρ ¯ α 1 2 ) , α 1 2 , lim ρ 0 ρ ¯ γ 2 ( ln ρ ln ρ ¯ ) , α = 1 2 { , 0 < α 1 2 , 2 2 α 1 ρ ¯ γ + 2 α 1 2 , α > 1 2 .
(3.30)

As the initial data satisfies (3.3), it follows from (3.5) and (3.7) that

| ψ ( ρ ( ξ , τ ) ) | | ψ ( ρ ( 1 , τ ) ) | + | 0 1 ξ ψ ( ρ ) d ξ | C ˜ | ρ α ( 1 , τ ) ρ ¯ α | + | 0 1 φ ( ρ ) 1 2 ρ α 1 ρ ξ d ξ | C ˜ | ρ 0 α ( 1 ) ρ ¯ α | + 1 ν ( 0 1 φ ( ρ ) d ξ ) 1 2 ( 0 1 1 α 2 ( ρ α ) ξ 2 d ξ ) 1 2 4 ν E 0 1 2 ( E 0 + 2 E 1 ) 1 2 ,
(3.31)

where C ˜ and ν are positive constants independent of time, and we assume that

| ρ 0 α (1) ρ ¯ α |min { 2 C ˜ ν E 0 1 2 ( E 0 + 2 E 1 ) 1 2 , δ 1 } :=δ.
(3.32)

As 0<α 1 2 and as α> 1 2 , from the condition

E 0 1 2 ( E 0 + 2 E 1 ) 1 2 ν 2 ( 2 α 1 ) ρ ¯ γ + 2 α 1 2 ,
(3.33)

we can find that there are two positive constants ρ and ρ independent of time and choose

ρ = ρ 2 , ρ + =2 ρ ,
(3.34)

such that

0< ρ < ρ ρ(ξ,τ) ρ < ρ + ,(ξ,τ)[0,1]×[0,T],T>0.
(3.35)

As the initial data satisfies (3.4), it follows from (3.5) and (3.8) that

| ψ ( ρ ( ξ , τ ) ) | | 0 1 ψ ( ρ ) d ξ | + | 0 1 ξ ψ ( ρ ) d ξ | | 0 1 φ ( ρ ) 1 2 ρ α d ξ | + | 0 1 φ ( ρ ) 1 2 ρ α 1 ρ ξ d ξ | 1 ν ( 0 1 φ ( ρ ) d ξ ) 1 2 ( 0 1 ρ 2 α d ξ ) 1 2 + 1 ν ( 0 1 φ ( ρ ) d ξ ) 1 2 ( 0 1 1 α 2 ( ρ α ) ξ 2 d ξ ) 1 2 2 ν ( 0 1 φ ( ρ ) d ξ ) 1 2 ( 0 1 1 α 2 ( ρ α ) ξ 2 d ξ ) 1 2 4 ν E 0 1 2 ( E 0 + E 1 + ( ρ 0 γ ( b 0 ) ρ ¯ γ ) b 0 ( ρ 0 γ ( a 0 ) ρ ¯ γ ) a 0 ) 1 2 ,
(3.36)

where ν is a positive constant independent of time, and we use the fact that

ρ α ( ξ , τ ) L 2 2 ρ α ( ξ , τ ) L 2 C ρ α ( ξ , τ ) L 2 ( ρ α ) ξ ( ξ , τ ) L 2 ,
(3.37)

which together with Young’s inequality gives

ρ α ( ξ , τ ) L 2 2 C ( ρ α ) ξ ( ξ , τ ) L 2 2 ,
(3.38)

where C is a positive constant independent of time. As 0<α 1 2 and as 1 2 <αγ, by means of the condition

E 0 1 2 ( E 0 + E 1 + ( ρ 0 γ ( b 0 ) ρ ¯ γ ) b 0 ( ρ 0 γ ( a 0 ) ρ ¯ γ ) a 0 ) 1 2 ν 2 ( 2 α 1 ) ρ ¯ γ + 2 α 1 2 ,
(3.39)

we can obtain two positive constants ρ and ρ independent of time such that

0< ρ ρ(ξ,τ) ρ ,(ξ,τ)[0,1]×[0,T],T>0.
(3.40)

The proof of this lemma is completed. □

We also have the regularity estimates for the solution (ρ,u) to FBVP (3.2) as follows.

Lemma 3.4 Let T>0. Under the assumptions of Theorem  2.1, it holds for any strong solution (ρ,u) to FBVP (3.2) that

{ ρ L ( 0 , T ; H 1 ( [ 0 , 1 ] ) ) C 0 ( [ 0 , 1 ] × [ 0 , T ] ) , u L ( 0 , T ; H 1 ( [ 0 , 1 ] ) ) L 2 ( 0 , T ; H 2 ( [ 0 , 1 ] ) ) , a ( τ ) , b ( τ ) H 1 ( [ 0 , T ] ) , ( ρ γ ρ 1 + α u ξ ) L ( 0 , T ; L 2 ( [ 0 , 1 ] ) ) .
(3.41)

If it is also satisfied that

u 0 H 2 ( [ 0 , 1 ] ) ,
(3.42)

then the strong solution (ρ,u) has the regularities

{ ( ρ , u ) C 0 ( [ 0 , 1 ] × [ 0 , T ] ) , ρ L ( 0 , T ; H 1 ( [ 0 , 1 ] ) ) , ρ τ L ( 0 , T ; L 2 ( [ 0 , 1 ] ) ) , u L ( 0 , T ; H 2 ( [ 0 , 1 ] ) ) L 2 ( 0 , T ; H 3 ( [ 0 , 1 ] ) ) , u τ L ( 0 , T ; L 2 ( [ 0 , 1 ] ) ) L 2 ( 0 , T ; H 1 ( [ 0 , 1 ] ) ) , a ( τ ) , b ( τ ) H 2 ( [ 0 , T ] ) , ( ρ γ ρ 1 + α u ξ ) C 0 ( [ 0 , T ] × ( [ 0 , 1 ] ) ) .
(3.43)

Proof Multiplying (3.2)2 by ρ ( 1 + α ) u τ , integrating the result over [0,1] and making use of the boundary conditions, after a direct computation and recombination, we have

d d τ 0 1 ( 1 2 u ξ 2 ( ρ γ ρ ¯ γ ) ρ ( 1 + α ) u ξ ) d ξ + 0 1 ρ ( 1 + α ) u τ 2 d ξ = [ γ ( 1 + α ) ] 0 1 ρ γ α u ξ 2 d ξ + ( 1 + α ) ρ ¯ γ 0 1 ρ α u ξ 2 d ξ ( 1 + α ) 0 1 ( ρ γ ρ ¯ γ ) ρ ( 2 + α ) ρ ξ u τ d ξ + ( 1 + α ) 0 1 ρ 1 ρ ξ u ξ u τ d ξ .
(3.44)

Integrating equation (3.44) over [0,τ], from (3.5), (3.7), (3.8) and (3.26), it is easily verified that

0 1 u ξ 2 d ξ + 0 τ 0 1 u s 2 d ξ d s C + C 0 τ 0 1 u ξ 2 d ξ d s + C 0 τ 0 1 ρ ξ 2 d ξ d s + C 0 τ 0 1 ρ ξ 2 u ξ 2 d ξ d s ,
(3.45)

where C denotes a constant independent of time. From (3.2)2, (3.5), (3.7), (3.8) and (3.26), we deduce

C 0 τ 0 1 ρ ξ 2 u ξ 2 d ξ d s C 3 0 τ 0 1 ρ ξ 2 u ξ 2 d ξ d s + 1 6 0 τ 0 1 u s 2 d ξ d s + C 0 τ 0 1 ρ ξ 2 d ξ d s + C 0 τ 0 1 u ξ 2 d ξ d s .
(3.46)

Using (3.46), we can obtain that

0 1 u ξ 2 dξ+ 0 τ 0 1 u s 2 dξdsC,
(3.47)

which together with (3.2)2 implies

0 1 u ξ 2 dξ+ 0 τ 0 1 u s 2 dξ+ 0 τ 0 1 u ξ ξ 2 dξdτC.
(3.48)

Differentiating (3.2)2 with respect to τ, we get

u τ τ + ( ρ γ ) ξ τ = ( ρ 1 + α u ξ ) ξ τ .
(3.49)

Taking product between (3.49) and u τ , integrating the results over [0,1] and using the boundary conditions (3.2)4, we have

1 2 d d τ 0 1 u τ 2 dξ= 0 1 ( ρ γ ρ ¯ γ ) τ u ξ τ dξ 0 1 ( ρ 1 + α u ξ ) τ u ξ τ dξ.
(3.50)

The terms on the right-hand side of (3.50) can be bounded respectively as follows:

0 1 ( ρ γ ρ ¯ γ ) τ u ξ τ d ξ = 0 1 γ ρ γ + 1 u ξ u ξ τ d ξ γ 2 d d τ 0 1 ρ γ + 1 u ξ 2 d ξ + C 0 1 ( ρ 1 + α u ξ 2 + ρ 2 γ α + 3 u ξ 4 ) d ξ
(3.51)

and

0 1 ( ρ 1 + α u ξ ) τ u ξ τ d ξ = 0 1 ( ( 1 + α ) ρ α ρ τ u ξ + ρ 1 + α u ξ τ ) u ξ τ d ξ 1 2 0 1 ρ 1 + α u ξ τ 2 d ξ + C 0 1 ρ 3 + α u ξ 4 d ξ .
(3.52)

Summing (3.50)-(3.52) together and making use of (3.26) and (3.48), we obtain

1 2 d d τ 0 1 u τ 2 d ξ + γ 2 d d τ 0 1 ρ γ + 1 u ξ 2 d ξ + 1 2 0 1 ρ 1 + α u ξ τ 2 d ξ C 0 1 u ξ 2 d ξ + C u ξ L ( [ 0 , 1 ] ) 2 0 1 u ξ 2 d ξ C + C u ξ L ( [ 0 , 1 ] ) 2 0 1 u ξ 2 d ξ .
(3.53)

Integrating equation (3.53) over [0,τ], we get from (3.2)2, (3.5) and (3.26) that

0 1 u τ 2 d ξ + 0 1 u ξ 2 d ξ + 0 τ 0 1 u ξ s 2 d ξ d s C + C sup τ [ 0 , T ] u ξ L 2 C + C sup τ [ 0 , T ] ( 0 1 u ξ 2 d ξ ) 1 2 ( 0 1 u ξ ξ 2 d ξ ) 1 2 C + 1 2 0 1 u τ 2 d ξ ,
(3.54)

which gives

0 1 u τ 2 dξ+ 0 1 u ξ 2 dξ+ 0 T 0 1 u ξ τ 2 dξdτC,
(3.55)

which implies ( ρ γ ρ 1 + α u ξ ) L (0,T; H 1 ([0,1])), and it follows from the definition of a (τ)=u(0,τ) and b (τ)=u(1,τ) that a(τ),b(τ) H 2 ([0,T]). The proof of this lemma is completed. □

Finally, we give the large time behaviors of the strong solution as follows.

Lemma 3.5 Let T>0. Under the assumptions of Theorem  2.1, it holds for any strong solution (ρ,u) to FBVP (3.2) that

( ρ ρ ¯ , u ) ( , τ ) L ( [ 0 , 1 ] ) C 1 e C 2 τ ,τ>0,
(3.56)

where C 1 and C 2 denote two positive constants independent of time.

Proof Applying (3.6) and (3.12) with modification, we can obtain

d d τ 0 1 ( 1 2 u 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) dξ+ 0 1 u ξ 2 dξ0
(3.57)

and

d d τ 0 1 ( 1 2 ( u + 1 α ( ρ α ) ξ ) 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ + 0 1 ρ ξ 2 d ξ ( ( ρ γ ( 1 , τ ) ρ ¯ γ ) u ( 1 , τ ) ( ρ γ ( 0 , τ ) ρ ¯ γ ) u ( 0 , τ ) ) .
(3.58)

We have from (3.5), (3.7), (3.26), the Gagliardo-Nirenberg-Sobolev inequality and Young’s inequality that

0 1 ( 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) d ξ C 0 1 ( ρ ρ ¯ ) 2 d ξ C 0 1 ρ ξ 2 d ξ ,
(3.59)
0 1 u 2 dξC 0 1 u ξ 2 dξ,
(3.60)

and define

E(τ):= 0 1 ( 1 4 u 2 + 1 4 ( u + r ( ρ α ) ξ ) 2 + 1 γ 1 ( ρ γ 1 ρ ¯ γ 1 ) + ρ ¯ γ ( ρ 1 ρ ¯ 1 ) ) dξ.
(3.61)

As the initial data satisfies (3.3), we have

| ( ρ γ ( 1 , τ ) ρ ¯ γ ) u(1,τ)|+| ( ρ γ ( 0 , τ ) ρ ¯ γ ) u(0,τ)|C e C ¯ τ ,
(3.62)

where C and C ¯ are positive constants independent of time. By (3.57)-(3.62), a complicated computation gives rise to

d d τ E(τ)+ C 0 E(τ)C e C ¯ τ ,
(3.63)

where C 0 < C ¯ is a positive constant independent of time. From (3.63), we have

E(τ)E(0) e C 0 τ +C e C 0 τ 0 τ e ( C ¯ C 0 ) s dsC e C 0 τ .
(3.64)

As the initial data satisfies (3.4), from (3.22)-(3.24), we obtain

ρ(1,τ)ρ(0,τ) ρ ¯ >0,
(3.65)

which implies that the domain (b(τ)a(τ)) expands as t grows up, so that we have

b (τ) a (τ)=u(1,τ)u(0,τ)0.
(3.66)

Then it holds from (3.65) and (3.66) that

( ρ γ ( 1 , τ ) ρ ¯ γ ) u(1,τ) ( ρ γ ( 0 , τ ) ρ ¯ γ ) u(0,τ)0.
(3.67)

Using (3.58), after a complicated computation, we have

d d τ E(τ)+ C 0 E(τ)0,
(3.68)

where C 0 is a positive constant independent of time, and we have

E(τ)E(0) e C 0 τ C e C 0 τ ,
(3.69)

where C is a positive constant independent of time.

By the fact

E(τ)c ( ( ρ ρ ¯ , u ) L 2 ( [ 0 , 1 ] ) 2 + ( ρ α ) ξ L 2 ( [ 0 , 1 ] ) 2 ) ,
(3.70)

where c>0 is a constant independent of time, and the Gagliardo-Nirenberg-Sobolev inequality

( ρ ρ ¯ , u ) L ( [ 0 , 1 ] ) ( ρ ρ ¯ , u ) L 2 ( [ 0 , 1 ] ) 1 2 ( ρ ρ ¯ , u ) ξ L 2 ( [ 0 , 1 ] ) 1 2 ,
(3.71)

we can deduce (3.56). □

4 Proof of the main results

Proof The global existence of a unique strong solution to FBVP (1.1) and (2.1) can be established in terms of the short time existence carried out as in [7], the uniform a priori estimates and the analysis of regularities, which indeed follow from Lemmas 3.1-3.4. We omit the details. The large time behaviors follow from Lemma 3.5 directly. The proof of Theorem 2.1 is completed. □

References

  1. Gerbeau JF, Perthame B: Derivation of viscous Saint-Venant system for laminar shallow water, numerical validation. Discrete Contin. Dyn. Syst., Ser. B 2001, 1: 89-102.

    Article  MathSciNet  MATH  Google Scholar 

  2. Marche F: Derivation of a new two-dimensional viscous shallow water model with varying topography, bottom friction and capillary effects. Eur. J. Mech. B, Fluids 2007, 26: 49-63. 10.1016/j.euromechflu.2006.04.007

    Article  MathSciNet  MATH  Google Scholar 

  3. Bresch D, Desjardins B: Existence of global weak solutions for a 2D viscous shallow water equations and convergence to the quasi-geostrophic model. Commun. Math. Phys. 2003, 238: 211-223.

    Article  MathSciNet  MATH  Google Scholar 

  4. Bresch D, Desjardins B: On the construction of approximate solutions for the 2D viscous shallow water model and for compressible Navier-Stokes models. J. Math. Pures Appl. 2006, 86: 362-368. 10.1016/j.matpur.2006.06.005

    Article  MathSciNet  MATH  Google Scholar 

  5. Fang D, Zhang T: Global solutions of the Navier-Stokes equations for compressible flow with density-dependent viscosity and discontinuous initial data. J. Differ. Equ. 2006, 222: 63-94. 10.1016/j.jde.2005.07.011

    Article  MathSciNet  MATH  Google Scholar 

  6. Guo Z, Li H, Xin Z: Lagrange structure and dynamics for solutions to the spherically symmetric compressible Navier-Stokes equations. Commun. Math. Phys. 2012, 309: 371-412. 10.1007/s00220-011-1334-6

    Article  MathSciNet  MATH  Google Scholar 

  7. Jiang S, Xin Z, Zhang P: Global weak solutions to 1D compressible isentropy Navier-Stokes with density-dependent viscosity. Methods Appl. Anal. 2005, 12: 239-252.

    MathSciNet  MATH  Google Scholar 

  8. Lian R, Guo Z, Li H: Dynamical behaviors of vacuum states for 1D compressible Navier-Stokes equations. J. Differ. Equ. 2010, 248: 1926-1954. 10.1016/j.jde.2009.11.029

    Article  MathSciNet  MATH  Google Scholar 

  9. Liu J: Local existence of solution to free boundary value problem for compressible Navier-Stokes equations. Acta Math. Sci., Ser. B 2012, 32: 1298-1320. 10.1016/S0252-9602(12)60100-3

    Article  MathSciNet  MATH  Google Scholar 

  10. Liu T, Xin Z, Yang T: Vacuum states for compressible flow. Discrete Contin. Dyn. Syst. 1998, 4: 1-32.

    MathSciNet  MATH  Google Scholar 

  11. Okada M, Nečasoá ŠM, Makino T: Free boundary problem for the equation of one-dimensional motion of compressible gas with density-dependent viscosity. Ann. Univ. Ferrara, Sez. 7: Sci. Mat. 2002, 48: 1-20.

    MathSciNet  MATH  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  MATH  Google Scholar 

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

    Article  MathSciNet  MATH  Google Scholar 

  14. 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  MathSciNet  MATH  Google Scholar 

  15. 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  MATH  Google Scholar 

  16. Mellet A, Vasseur A: Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations. SIAM J. Math. Anal. 2008, 39: 1344-1365. 10.1137/060658199

    Article  MathSciNet  MATH  Google Scholar 

  17. Guo Z, Jiu Q, Xin Z: Spherically symmetric isentropic compressible flows with density-dependent viscosity coefficients. SIAM J. Math. Anal. 2008, 39: 1402-1427. 10.1137/070680333

    Article  MathSciNet  MATH  Google Scholar 

  18. Jiu Q, Wang Y, Xin Z: Stability of rarefaction waves to the 1D compressible Navier-Stokes equations with density-dependent viscosity. Commun. Partial Differ. Equ. 2011, 36: 602-634. 10.1080/03605302.2010.516785

    Article  MathSciNet  MATH  Google Scholar 

  19. Jiu Q, Xin Z: The Cauchy problem for 1D compressible flows with density-dependent viscosity coefficients. Kinet. Relat. Models 2008, 1: 313-330.

    Article  MathSciNet  MATH  Google Scholar 

  20. Li H, Li J, Xin Z: Vanishing of vacuum states and blow-up phenomena of the compressible Navier-Stokes equations. Commun. Math. Phys. 2008, 281: 401-444. 10.1007/s00220-008-0495-4

    Article  MathSciNet  MATH  Google Scholar 

  21. Lian R, Liu J, Li H, Xiao L: Cauchy problem for the one-dimensional compressible Navier-Stokes equations. Acta Math. Sci., Ser. B 2012, 32: 315-324. 10.1016/S0252-9602(12)60019-8

    Article  MathSciNet  MATH  Google Scholar 

  22. Xin Z: Blowup of smooth solutions to the compressible Navier-Stokes equation with compact density. Commun. Pure Appl. Math. 1998, 51: 229-240. 10.1002/(SICI)1097-0312(199803)51:3<229::AID-CPA1>3.0.CO;2-C

    Article  MathSciNet  MATH  Google Scholar 

  23. Solonnikov VA, Tani A: Free boundary problem for a viscous compressible flow with a surface tension. II. In Constantin Carathéodory: An International Tribute. World Scientific, Teaneck; 1991:1270-1303.

    Chapter  Google Scholar 

  24. Zajaczkowski WM: Existence of local solutions for free boundary problems for viscous compressible barotropic fluids. Ann. Pol. Math. 1995, 60: 255-287.

    MathSciNet  MATH  Google Scholar 

  25. Zajaczkowski WM: On nonstationary motion of a compressible barotropic viscous fluid bounded by a free surface. Diss. Math. 1993, 324: 1-101.

    MathSciNet  MATH  Google Scholar 

  26. Secchi P, Valli A: A free boundary problem for compressible viscous fluids. J. Reine Angew. Math. 1983, 341: 1-31.

    MathSciNet  MATH  Google Scholar 

  27. Tani A: On the free boundary value problem for compressible viscous fluid motion. J. Math. Kyoto Univ. 1981, 21: 839-859.

    MathSciNet  MATH  Google Scholar 

  28. Zadrzyńska E, Zajaczkowski WM: On local motion of a general compressible viscous heat conducting fluid bounded by a free surface. Ann. Pol. Math. 1994, 59: 133-170.

    MathSciNet  MATH  Google Scholar 

  29. Solonnikov VA, Tani A: Evolution free boundary problem for equations of motion of viscous compressible barotropic liquid. Lecture Notes in Math. 1530. In The Navier-Stokes Equations II - Theory and Numerical Methods. Springer, Berlin; 1992:30-55. Oberwolfach, 1991

    Chapter  Google Scholar 

  30. Zajaczkowski WM: On nonstationary motion of a compressible barotropic viscous capillary fluid bounded by a free surface. SIAM J. Math. Anal. 1994, 25: 1-84. 10.1137/S0036141089173207

    Article  MathSciNet  MATH  Google Scholar 

  31. Zadrzyńska E: Evolution free boundary problem for equations of viscous compressible heat-conducting capillary fluids. Math. Methods Appl. Sci. 2001, 24: 713-743. 10.1002/mma.238

    Article  MathSciNet  MATH  Google Scholar 

  32. Zadrzyńska E, Zajaczkowski WM: On the global existence theorem for a free boundary problem for equations of a viscous compressible heat conducting capillary fluid. J. Appl. Anal. 1996, 2: 125-169.

    MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors thank the referee for the helpful comments and suggestions on the paper. The research of Ruxu Lian is supported by NNSFC No. 11101145, China Postdoctoral Science Foundation No. 2012M520360, Doctoral Foundation of North China University of Water Sources and Electric Power No. 201032, Innovation Scientists and Technicians Troop Construction Projects of Henan Province. The research of Jian Liu is supported by NNSFC No. 11326140, the Doctoral Starting up Foundation of QuZhou University No. BSYJ201314.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ruxu Lian.

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. All authors 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lian, R., Liu, J. On the free boundary value problem for one-dimensional compressible Navier-Stokes equations with constant exterior pressure. Bound Value Probl 2014, 93 (2014). https://doi.org/10.1186/1687-2770-2014-93

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-93

Keywords