Skip to main content

3-D flow of a compressible viscous micropolar fluid with spherical symmetry: a local existence theorem

Abstract

We consider nonstationary 3-D flow of a compressible viscous heat-conducting micropolar fluid in the domain to be the subset of R 3 bounded with two concentric spheres that present solid thermoinsulated walls. In thermodynamical sense fluid is perfect and polytropic. Assuming that the initial density and temperature are strictly positive we will prove that for smooth enough spherically symmetric initial data there exists a spherically symmetric generalized solution locally in time.

1 Introduction

The theory of micropolar fluids is introduced in [1] by Eringen. Various problems with different initial and boundary conditions for incompressible micropolar fluid are presented in [2], but the theory of compressible micropolar fluid is still in the beginning. N. Mujakovic in [3] developed the model for one-dimensional isotropic, viscous, compressible micropolar fluid which is in thermodynamical sense perfect and polytropic. In the same work, the local existence of the solution for homogeneous boundary conditions is proved. N. Mujakovic in [4] and in the references cited therein proved the local and global existence of inhomogeneous boundary conditions for velocity and microrotation as well as stabilization and regularity. In [5] the Cauchy problem for the described problem was also considered. In the last years we find some interesting works with different kind of problems concerning micropolar fluid, e.g., [6, 7], but till now the described model of compressible micropolar fluid in three-dimensional case has not been considered.

In this work we consider the three-dimensional model with spherical symmetry. The first article in which the problem of spherical symmetry was described is [8], but for classical fluid. The spherical symmetry for classical fluid is also considered in articles [912].

In the setting of the field equations we use the Eulerian description.

In what follows we use the notation:

ρ - mass density

v - velocity

p - pressure

T - stress tensor

T x - an axial vector with the Cartesian components ( T x ) i = ε i j k T j k , where ε i j k is Levi-Civita alternating tensora

symT= 1 2 (T+ T T ), skwT= 1 2 (T T T )

ω - microrotation velocity

ω skw - skew tensor with Cartesian components ( ω skw ) i j = ε k i j ω k

j I - microinertia density (a positive scalar field)

C - couple stress tensor

θ - absolute temperature

E - internal energy density

q - heat flux density vector

f - body force density

g - body couple density

δ - body heat density

The problem we consider here is based on local forms of the conservation laws for mass, momentum, momentum moment and energy, which are stated respectively as follows:

ρ ˙ +ρdivv=0,
(1)
ρ v ˙ =divT+ρf,
(2)
ρ j I ω ˙ =divC+ T x +ρg,
(3)
ρ E ˙ =T:v+C:ω T x ωdivq+ρδ,
(4)

where a ˙ denotes material derivative of a field a:

a ˙ = a t +(a)v.

The scalar product of tensors A and B is defined by

A:B=tr ( A T B ) .

The linear constitutive equations for stress tensor, couple stress tensor and heat flux density vector are respectively in the forms:

T=(p+λdivv)I+2μsymv2 μ r skwv2 μ r ω skw ,
(5)
C= c 0 (divω)I+2 c d symω2 c a skwω,
(6)
q=kθ,
(7)

where

λ, μ - coefficients of viscosity,

μ r , c 0 , c d , c a - coefficients of microviscosity,

k - heat conduction coefficient

are constants with the properties

μ0,3λ+2μ0, μ r 0,
(8)
c d 0,3 c 0 +2 c d 0,| c d c a | c d + c a ,k0.
(9)

Assuming that the fluid is perfect and polytropic, for pressure and internal energy we have the equations

p=Rρθ,
(10)
E= c v θ,
(11)

where R and c v are positive constants.

Let Ω={x R 3 ,a<|x|<b}, a>0, denote the domain bounded by two concentric spheres with radii a and b. The boundary of the described domain is Ω={x R 3 ,|x|=a or |x|=b}. We shall consider the problem (1)-(11) in the region Q T =Ω×]0,T[ (where T>0 is arbitrary) with the following initial conditions:

ρ(x,0)= ρ 0 (x),v(x,0)= v 0 (x),ω(x,0)= ω 0 (x),θ(x,0)= θ 0 (x),
(12)

for xΩ and boundary conditions

v | Ω =0,ω | Ω =0, θ ν | Ω =0,
(13)

for 0<t<T; the vector ν is an exterior unit normal vector.

For simplicity we also assume that f=g=0 and δ=0.

The initial boundary problems for the system (1)-(13) so far have not been considered in three-dimensional case. The same and similar models in one-dimensional case were considered in [3, 5, 13] and [4]. In [2] the three-dimensional model was considered but for an incompressible micropolar fluid.

In this paper we prove the local existence of generalized spherically symmetric solution to the problem (1)-(13) in the domain Ω, assuming that the initial functions are also spherically symmetric. In the proof we use the Faedo-Galerkin method. We follow some ideas of [14] where this method was applied to a classical fluid (where microrotation is equal to zero) in one-dimensional case as well as the ideas from [3] and [13] where the same result as here was provided for one-dimensional case.

The paper is organized as follows. In the second section, we derive a spherically symmetric form of (1)-(4), introduce Lagrangian description, and present the main result. In the third section, we consider an approximate problem and get an approximate solution for each nN. In the forth section, we prove uniform a priori estimates for the approximate solutions. The proof of the main result is given in the fifth section.

2 Spherically symmetric form and the main result

We first derive the spherically symmetric form of (1)-(7) and (10)-(11). A spherically symmetric solution of (1)-(7) has the form:

v i ( x , t ) = x i r v ( r , t ) , ω i ( x , t ) = x i r ω ( r , t ) , i = 1 , 2 , 3 , ρ ( x , t ) = ρ ( r , t ) , θ ( x , t ) = θ ( r , t ) ,
(14)

where x=( x 1 , x 2 , x 3 ) R 3 , r=|x|, v=( v 1 , v 2 , v 3 ) and ω=( ω 1 , ω 2 , ω 3 ). We assume that

ρ 0 (x)= ρ 0 (r), v 0 (x)= x r v 0 (r), ω 0 (x)= x r ω 0 (r), θ 0 (x)= θ 0 (r),
(15)

where ρ 0 , v 0 , ω 0 and θ 0 are known real functions defined on ]a,b[, and thus we reduce system (1)-(7) and conditions (10)-(13) to the following equations for ρ(r,t), v(r,t), ω(r,t) and θ(r,t) of the form:

ρ t + r (vρ)+ 2 ρ r v=0,
(16)
ρ ( v t + v v r ) =R r (ρθ)+(λ+2μ) r ( v r + 2 v r ) ,
(17)
ρ j I ( ω t + v ω r ) =4 μ r ω+( c 0 +2 c d ) r ( ω r + 2 ω r ) ,
(18)
(19)

with the following initial and boundary conditions

ρ ( r , 0 ) = ρ 0 ( r ) , v ( r , 0 ) = v 0 ( r ) , ω ( r , 0 ) = ω 0 ( r ) , θ ( r , 0 ) = θ 0 ( r ) , r ] a , b [ ,
(20)
v ( a , t ) = v ( b , t ) = 0 , ω ( a , t ) = ω ( b , t ) = 0 , θ r ( a , t ) = θ r ( b , t ) = 0 , 0 < t < T .
(21)

To investigate the local existence, it is convenient to transform the system (16)-(19) to that in Lagrangian coordinates. The Eulerian coordinates (r,t) are connected to the Lagrangian coordinates (ξ,t) by the relation

r(ξ,t)= r 0 (ξ)+ 0 t v ˜ (ξ,t)dτ, r 0 (ξ)=r(ξ,0)=ξ,
(22)

where v ˜ (ξ,t) is defined by

v ˜ (ξ,t)=v ( r ( ξ , t ) , t ) .
(23)

We introduce the new function η by

η(ξ)= a ξ s 2 ρ 0 (s)ds.
(24)

Note that if ρ 0 (s)>0 for s]a,b[ (which is assumed in Theorem 2.1 later), then there exists an inverse function η 1 . Let the constant L be defined as

η(b)= a b s 2 ρ 0 (s)ds=L.
(25)

From (16) we can easily get the equation

ρ r ξ r 2 = ρ 0 (ξ) ξ 2 ,
(26)

i.e.,

a r ( ξ , t ) ρ(s,t) s 2 ds= a ξ ρ 0 (s) s 2 ds.
(27)

It is useful to introduce the next coordinate

x = L 1 η(ξ)
(28)

and the following functions

ρ ( x , t ) = ρ ˜ ( η 1 ( x L ) , t ) , v ( x , t ) = v ˜ ( η 1 ( x L ) , t ) , ω ( x , t ) = ω ˜ ( η 1 ( x L ) , t ) , θ ( x , t ) = θ ˜ ( η 1 ( x L ) , t ) , r ( x , t ) = r ( η 1 ( x L ) , t ) , ρ 0 ( x ) = ρ 0 ( η 1 ( x L ) ) , v 0 ( x ) = v 0 ( η 1 ( x L ) ) , ω 0 ( x ) = ω 0 ( η 1 ( x L ) ) , θ 0 ( x ) = θ 0 ( η 1 ( x L ) ) , r 0 ( x ) = r 0 ( η 1 ( x L ) ) = η 1 ( x L ) .

Similarly as in [15], for a new coordinate x we get

0= L 1 η(a) x L 1 η(b)=1.
(29)

Taking into account (26) and (24), we obtain that the functions ρ , v , ω , θ and r satisfy the system that we write omitting the primes for simplicity:

ρ t = 1 L ρ 2 x ( r 2 v ) ,
(30)
v t = R L r 2 x (ρθ)+ λ + 2 μ L 2 r 2 x ( ρ x ( r 2 v ) ) ,
(31)
ρ ω t = 4 μ r j I ω+ c 0 + 2 c d j I L 2 r 2 ρ x ( ρ x ( r 2 ω ) ) ,
(32)
(33)

in ]0,1[×]0,T[, where T>0 is arbitrary. Now we have the following boundary and initial conditions

v(0,t)=v(1,t)=0,
(34)
ω(0,t)=ω(1,t)=0,
(35)
θ x (0,t)= θ x (1,t)=0,
(36)

for t]0,T[,

ρ(x,0)= ρ 0 (x),v(x,0)= v 0 (x),
(37)
ω(x,0)= ω 0 (x),θ(x,0)= θ 0 (x),
(38)

for x]0,1[. We also have

r(x,t)= r 0 (x)+ 0 t v(x,τ)dτ,(x,t)]0,1[×]0,T[.
(39)

From

r x (x,t)= L ρ ( x , t ) r 2 ( x , t )

putting t=0 and integrating over ]0,x[, we get

r 0 (x)= ( a 3 + 3 L 0 x 1 ρ 0 ( y ) d y ) 1 3 ,x]0,1[,
(40)

where a>0 is a radius of the smaller boundary sphere.

We assume the inequalities

ρ 0 (x)m, θ 0 (x)mfor x]0,1[,
(41)

where m R + .

Before stating the main result, we introduce the following definition.

Definition 2.1 A generalized solution of the problem (30)-(38) in the domain Q T =]0,1[×]0,T[ is a function

(x,t)(ρ,v,ω,θ)(x,t),(x,t) Q T ,
(42)

where

ρ L ( 0 , T ; H 1 ( ] 0 , 1 [ ) ) H 1 ( Q T ), inf Q T ρ>0,
(43)
v,ω,θ L ( 0 , T ; H 1 ( ] 0 , 1 [ ) ) H 1 ( Q T ) L 2 ( 0 , T ; H 2 ( ] 0 , 1 [ ) ) ,
(44)

that satisfies Equations (30)-(33) a.e. in Q T and conditions (34)-(38) in the sense of traces.

Remark 2.1 From the embedding and interpolation theorems (e.g., [16] and [17]) one can conclude that from (43) and (44) it follows:

ρ L ( 0 , T ; C ( [ 0 , 1 ] ) ) C ( [ 0 , T ] , L 2 ( ] 0 , 1 [ ) ) ,
(45)
v,ω,θ L 2 ( 0 , T ; C ( 1 ) ( [ 0 , 1 ] ) ) C ( [ 0 , T ] , H 1 ( ] 0 , 1 [ ) ) ,
(46)
v,ω,θC( Q ¯ T ).
(47)

It is easy to check that the solution (42) with properties (43)-(44) satisfies the condition for a strong solution of the described problem.

The aim of this paper is to prove the following statements.

Theorem 2.1 Let the functions

ρ 0 , θ 0 H 1 ( ] 0 , 1 [ ) ,
(48)
v 0 , ω 0 H 0 1 ( ] 0 , 1 [ )
(49)

satisfy conditions (41). Then there exists T 0 , 0< T 0 T, such that the problem (30)-(38) has a generalized solution in Q 0 = Q T 0 , having the property

θ>0in  Q ¯ 0 .
(50)

For the function r, it holds

r L ( 0 , T ; H 2 ( ] 0 , 1 [ ) ) H 2 ( Q 0 )C( Q ¯ 0 ),
(51)
a 2 r2Min  Q ¯ 0 .
(52)

Remark 2.2 Notice that the function r 0 introduced by (40) belongs to H 2 (]0,1[). Because of the embedding H 1 (]0,1[)C([0,1]) we can conclude that there exists M R + such that

ρ 0 (x), v 0 (x), ω 0 (x), θ 0 (x)M,x[0,1].
(53)

From (40) and (41) we get

r 0 C ( 1 ) ( [ 0 , 1 ] ) ,
(54)
0<a r 0 (x)M,0< a 1 r 0 (x) M 1 ,x[0,1],
(55)

where a 1 =L M 3 and M 1 =L ( m a 2 ) 1 .

The proof of Theorem 2.1 is essentially based on a careful examination of a priori estimates and a limit procedure. We first study, for each nN, an approximate problem and derive the a priori estimates for approximate solutions independent of n by utilizing a technique of Kazhikov [14, 18] and Mujakovic [3, 13] for one-dimensional case. Using the obtained a priori estimates and results of weak compactness, we extract the subsequence of approximate solutions, which, when n tends to infinity, has limit in the same weak sense on ]0,1[×]0, T 0 [ for sufficiently small T 0 , 0< T 0 T. Finally, we show that this limit is the solution to our problem.

3 Approximate solutions

We shall find a local generalized solution to the problem (30)-(38) as a limit of approximate solutions

( ρ n , v n , ω n , θ n ) ,nN,
(56)

obtained in what follows. First, we introduce the approximations v n and r n of the functions v and r by

v n (x,t)= i = 1 n v i n (t)sin(πix),
(57)
r n (x,t)= r 0 (x)+ 0 t v n (x,τ)dτ,
(58)

where r 0 (x) is defined by (40) and v i n , i=1,2,,n are unknown smooth functions defined on an interval [0, T n ], T n T.

Then, we can write the solution ρ n to the problem

ρ n t + L 1 ( ρ n ) 2 x ( ( r n ) 2 v n ) =0, ρ n (x,0)= ρ 0 (x),
(59)

in the similar way as in [3] and [13] in the form

ρ n (x,t)= L ρ 0 ( x ) L + ρ 0 ( x ) x 0 t ( r n ) 2 v n d τ .
(60)

Since r n and v n are sufficiently smooth functions, we can conclude that the function ρ n is continuous on the rectangle [0,1]×[0, T n ] with the property ρ n (x,0)= ρ 0 (x)m>0. Because of aforementioned, we can conclude that there exists such T n , 0< T n T that

ρ n (x,t)>0,for (x,t)[0,1]×[0, T n ].
(61)

We also introduce the approximations ω n and θ n of the functions ω and θ respectively by

ω n (x,t)= j = 1 n ω j n (t)sin(πjx),
(62)
θ n (x,t)= k = 0 n θ k n (t)cos(πkx);
(63)

where ω j n and θ k n are again unknown smooth functions defined on an interval [0, T n ], T n T.

Evidently, the boundary conditions

v n (0,t)= v n (1,t)= ω n (0,t)= ω n (1,t)= θ n x (0,t)= θ n x (1,t)=0
(64)

for t]0, T n [ are satisfied.

According to the Faedo-Galerkin method, we take the following approximate conditions:

0 1 [ v n t + R L ( r n ) 2 x ( ρ n θ n ) λ + 2 μ L 2 ( r n ) 2 x ( ρ n x ( ( r n ) 2 v n ) ) ] sin ( π i x ) d x = 0 ,
(65)
0 1 [ ω n t + 4 μ r j I ω n ρ n c 0 + 2 c d j I L 2 ( r n ) 2 x ( ρ n x ( ( r n ) 2 ω n ) ) ] sin ( π j x ) d x = 0 ,
(66)
0 1 [ θ n t k c v L 2 x ( ( r n ) 4 ρ n θ n x ) + R c v L ρ n θ n x ( ( r n ) 2 v n ) λ + 2 μ c v L 2 ρ n [ x ( ( r n ) 2 v n ) ] 2 + 4 μ c v L x ( r n ( v n ) 2 ) c 0 + 2 c d c v L 2 ρ n [ x ( ( r n ) 2 ω n ) ] 2 + 4 c d c v L x ( r n ( ω n ) 2 ) 4 μ r c v ( ω n ) 2 ρ n ] cos ( π k x ) d x = 0
(67)

for i,j=1,,n, k=0,1,,n.

Let v 0 i , ω 0 j , and θ 0 k be the Fourier coefficients of the functions v 0 , ω 0 , and θ 0 respectively:

v 0 i = 2 0 1 v 0 ( x ) sin ( π i x ) d x , i = 1 , , n , ω 0 j = 2 0 1 ω 0 ( x ) sin ( π j x ) d x , j = 1 , , n , θ 00 = 0 1 θ 0 ( x ) d x , θ 0 k = 2 0 1 θ 0 ( x ) cos ( π k x ) d x , k = 1 , , n .

Let v 0 n , ω 0 n and θ 0 n be

v 0 n (x)= i = 1 n v 0 i sin(πix),
(68)
ω 0 n (x)= j = 1 n ω 0 j sin(πjx),
(69)
θ 0 n (x)= k = 0 n θ 0 k cos(πkx).
(70)

We take the initial conditions for v n , ω n and θ n in the form

v n (x,0)= v 0 n (x),
(71)
ω n (x,0)= ω 0 n (x),
(72)
θ n (x,0)= θ 0 n (x).
(73)

Let z m n , λ p q n and μ s l g n be

z m n (t)= 0 t v m n (τ)dτ,m=1,,n,
(74)
λ p q n (t)= 0 t z p n (τ) v q n (τ)dτ,p,q=1,,n,
(75)
μ s l g n (t)= 0 t z l n (τ) z s n (τ) v g n (τ)dτ,s,l,g=1,,n
(76)

then we have

r n (x,t)= r 0 (x)+ m = 1 n z m n (t)sin(πmx),
(77)
(78)

where r 0 (x) and ρ 0 (x) are known functions. Taking into account (57), (62), (63), (74)-(78), from (65)-(67) we obtain for {( v i n , ω j n , θ k n , z m n , λ p q n , μ s l g n ):i,j,m,p,q,s,l,g=1,,n,k=0,1,,n}, the following Cauchy problem:

(79)
(80)
(81)
z ˙ m n (t)= v m n ,
(82)
λ ˙ p q n (t)= z p n v q n ,
(83)
μ ˙ s l g n (t)= z s n z l n v g n ,
(84)
v i n (0)= v 0 i , ω j n (0)= ω 0 j , θ k n (0)= θ 0 k ,
(85)
z m n (0)=0, λ p q n (0)=0, μ s l g n (0)=0.
(86)

Here we have λ 0 =1, λ k =2 for k=1,2,,n and

Φ i n = 2 0 1 [ λ + 2 μ L 2 ( r n ) 2 x ( ρ n x ( ( r n ) 2 v n ) ) R L ( r n ) 2 x ( ρ n θ n ) ] sin ( π i x ) d x ,
(87)
Ψ j n = 2 0 1 [ c 0 + 2 c d j I L 2 ( r n ) 2 x ( ρ n x ( ( r n ) 2 ω n ) ) 4 μ r j I ω n ρ n ] sin ( π j x ) d x ,
(88)
Π k n = 0 1 [ k c v L 2 x ( ( r n ) 4 ρ n θ n x ) R c v L ρ n θ n x ( ( r n ) 2 v n ) + λ + 2 μ c v L 2 ρ n [ x ( ( r n ) 2 v n ) ] 2 4 μ c v L x ( r n ( v n ) 2 ) + c 0 + 2 c d c v L 2 ρ n [ x ( ( r n ) 2 ω n ) ] 2 4 c d c v L x ( r n ( ω n ) 2 ) + 4 μ r c v ( ω n ) 2 ρ n ] cos ( π k x ) d x .
(89)

Notice that the functions on the right-hand side of (79)-(84) satisfy the conditions of the Cauchy-Picard theorem [19, 20] and we can easily conclude that the following statements are valid.

Lemma 3.1 For each nN there exists such T n , 0< T n T that the Cauchy problem (79)-(86) has a unique solution defined on [0, T n ]. The functions v n , ω n and θ n defined by the formulas (57), (62) and (63) belong to the class C ( Q ¯ n ), Q n =]0,1[×]0, T n [ and satisfy conditions (71)-(73).

From (77) and (78) we can also easily conclude that

ρ n C( Q ¯ n ),
(90)
r n C ( 1 ) ( Q ¯ n ).
(91)

Lemma 3.2 There exists T n , 0< T n T, such that the functions ρ n , r n and r n x satisfy the conditions

m 2 ρ n (x,t)2M,
(92)
a 2 r n (x,t)2M,
(93)
a 1 2 r n x (x,t)2 M 1 ,
(94)

on Q ¯ n . The constants m, a, a 1 , M and M 1 are introduced by (40), (41), (53) and (55).

Proof The statements follow from (90)-(91), (41), (53) and (55). □

4 A priori estimates

Our purpose is to find T 0 , 0< T 0 T such that for each nN there exists a solution to the problem (79)-(86), defined on [0, T 0 ]. It will be sufficient to find uniform (in nN) a priori estimates for the solution ( ρ n , v n , ω n , θ n ) defined through Lemmas 3.1 and 3.2.

In what follows we denote by C>0 or C i >0 (i=1,2,) a generic constant, not depending on nN, which may have different values at different places.

We also use the notation

f= f L 2 ( ] 0 , 1 [ ) .

Some of our considerations are very similar or identical to that of [3] or [13]. In these cases we omit proofs or details of proofs making references to corresponding pages of the articles [3] or [13].

Lemma 4.1 For t[0, T n ] it holds

2 r n x 2 ( t ) 2 C ( 1 + 0 t 2 v n x 2 ( τ ) 2 d τ ) .
(95)

Proof From (58) follows

2 r n x 2 = r 0 + 0 t 2 v n x 2 dτ

and using Remark 2.2 we get (95) immediately. □

Lemma 4.2 For t[0, T n ], the following inequality holds:

ω n ( t ) 2 + 0 t ( ω n ( τ ) 2 + x ( r n 2 ω n ) ( τ ) 2 d τ ) C.
(96)

Proof Multiplying (66) by ω j n and summing over j=1,,n, after integration by parts, we obtain

1 2 d d t ω n ( t ) 2 + 4 μ r j I 0 1 ( ω n ) 2 ρ n dx+ c 0 + 2 c d j I L 2 0 1 ρ n [ x ( ( r n ) 2 ω n ) ] 2 dx=0.

Integrating over [0,t], 0<t T n , and taking into account (72), we obtain

1 2 ω n ( t ) 2 + 4 μ r j I 0 t 0 1 ( ω n ) 2 ρ n d x d τ + c 0 + 2 c d j I L 2 0 t 0 1 ρ n [ x ( ( r n ) 2 ω n ) ] 2 d x d τ = 1 2 ω 0 n 2 1 2 ω 0 2 .

Using (92), we get (96). □

In what follows, we use the inequalities

| f | 2 2f f x ,| f x | 2 2 f x 2 f x 2
(97)

(for a function f vanishing at x=0 and x=1 and with the first derivative vanishing at some point x[0,1]) that satisfy the functions v n , ω n and θ n x .

Lemma 4.3 For t[0, T n ], the following inequality holds:

| 0 1 θ n (x,t)dx|C ( 1 + v n x ( t ) 2 ) .
(98)

Proof Multiplying (65) by c v 1 v i n and summing over i=1,,n, after integration by parts and adding to (67) for k=0, we have

d d t ( 1 2 c v v n ( t ) 2 + 0 1 θ n ( x , t ) d x ) = c 0 + 2 c d c v L 2 0 1 ρ n [ x ( ( r n ) 2 ω n ) ] 2 d x + 4 μ r c v 0 1 ( ω n ) 2 ρ n d x .

Integrating over [0,t], 0<t T n and using (92) we get

( 1 2 c v v n ( t ) 2 + 0 1 θ n ( x , t ) d x ) C 0 t ( ω n ( t ) 2 + x ( ( r n ) 2 ω n ) ( τ ) 2 ) d τ + 1 2 c v v 0 n ( t ) 2 + θ 0 n ( t ) .

Taking into account (96), (71), (73), and (97) we obtain (98). □

Lemma 4.4 ([3], Lemma 5.3)

For (x,t) Q ¯ n , the following inequality holds:

| θ n (x,t)|C ( 1 + θ n x ( t ) + v n x ( t ) 2 ) .
(99)

Lemma 4.5 For t[0, T n ], the following inequality holds:

ρ n x ( t ) 2 C ( 1 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 2 ) .
(100)

Proof Taking the derivative of the function ρ n with respect to x and using the estimates (92)-(94), we obtain

| ρ n x |C ( 1 + 0 t ( | v n | + | v n | | 2 r n x 2 | + | v n x | + | 2 v n x 2 | ) d τ ) .

With the help of (97) applied to the function v n , the Hoelder and Young inequalities as well as (95), we get (100). □

Lemma 4.6 For t[0, T n ] it holds

d d t ( v n x ( t ) 2 + ω n x ( t ) 2 + θ n x ( t ) 2 ) + B ( 2 v n x 2 ( t ) 2 + 2 ω n x 2 ( t ) 2 + 2 θ n x 2 ( t ) 2 ) C ( 1 + v n x ( t ) 16 + ω n x ( t ) 16 + θ n x ( t ) 16 + ( B 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) ,
(101)

where

B= m a 4 32 min { λ + 2 μ L 2 , c 0 + 2 c d j I L 2 , k c v L 2 } .

Proof As in [3] Lemma 5.5, [14] pp.63-66 and in [13] Lemma 5.6, multiplying (65), (66) and (67) respectively by ( π i ) 2 v i n , ( π j ) 2 ω j n and ( π k ) 2 θ k n and taking into account (57), (62) and (63), after summation over i,j,k=1,2,,n and addition of the obtained equations, we get

1 2 d d t ( v n x 2 + ω n x 2 + θ n x 2 ) + λ + 2 μ L 2 0 1 ρ n ( r n ) 4 ( 2 v n x 2 ) 2 d x + c 0 + 2 c d j I L 2 0 1 ρ n ( r n ) 4 ( 2 ω n x 2 ) 2 d x + k c v L 2 0 1 ρ n ( r n ) 4 ( 2 θ n x 2 ) 2 d x = p = 1 28 I p ( t ) ,
(102)

where

I 1 = 2 ( λ + 2 μ ) L 2 0 1 ( r n ) 3 ρ n x r n x v n 2 v n x 2 d x , I 2 = 2 ( λ + 2 μ ) L 2 0 1 ρ n ( r n ) 2 ( r n x ) 2 v n 2 v n x 2 d x , I 3 = 2 ( λ + 2 μ ) L 2 0 1 ρ n ( r n ) 3 2 r n x 2 v n 2 v n x 2 d x , I 4 = 4 ( λ + 2 μ ) L 2 0 1 ρ n ( r n ) 3 r n x v n x 2 v n x 2 d x , I 5 = λ + 2 μ L 2 0 1 ( r n ) 4 ρ n x v n x 2 v n x 2 d x , I 6 = R L 0 1 ( r n ) 2 θ n ρ n x 2 v n x 2 d x , I 7 = R L 0 1 ( r n ) 2 ρ n θ n x 2 v n x 2 d x , I 8 = 4 μ r j I 0 1 ω n ρ n 2 ω n x 2 d x , I 9 = 2 ( c 0 + 2 c d ) j I L 2 0 1 ( r n ) 3 ρ n x r n x ω n 2 ω n x 2 d x , I 10 = 2 ( c 0 + 2 c d ) j I L 2 0 1 ρ n ( r n ) 2 ( r n x ) 2 ω n 2 ω n x 2 d x , I 11 = 2 ( c 0 + 2 c d ) j I L 2 0 1 ρ n ( r n ) 3 2 r n x 2 ω n 2 ω n x 2 d x , I 12 = 4 ( c 0 + 2 c d ) j I L 2 0 1 ρ n ( r n ) 3 r n x ω n x 2 ω n x 2 d x , I 13 = c 0 + 2 c d j I L 2 0 1 ( r n ) 4 ρ n x ω n x 2 ω n x 2 d x , I 14 = 4 k c v L 2 0 1 ρ n ( r n ) 3 r n x θ n x 2 θ n x 2 d x , I 15 = k c v L 2 0 1 ( r n ) 4 ρ n x θ n x 2 θ n x 2 d x , I 16 = 2 R c v L 0 1 ρ n r n r n x v n θ n 2 θ n x 2 d x , I 17 = R c v L 0 1 ρ n ( r n ) 2 θ n v n x 2 θ n x 2 d x , I 18 = 4 ( λ + 2 μ ) c v L 2 0 1 ρ n ( r n ) 2 ( r n x ) 2 ( v n ) 2 2 θ n x 2 d x , I 19 = 4 ( λ + 2 μ ) c v L 2 0 1 ρ n ( r n ) 3 r n x v n v n x 2 θ n x 2 d x , I 20 = λ + 2 μ c v L 2 0 1 ρ n ( r n ) 4 ( v n x ) 2 2 θ n x 2 d x , I 21 = 4 μ c v L 0 1 ( v n ) 2 r n x 2 θ n x 2 d x , I 22 = 8 μ c v L 0 1 r n v n v n x 2 θ n x 2 d x , I 23 = 4 ( c 0 + 2 c d ) c v L 2 0 1 ρ n ( r n ) 2 ( r n x ) 2 ( ω n ) 2 2 θ n x 2 d x , I 24 = 4 ( c 0 + 2 c d ) c v L 2 0 1 ρ n ( r n ) 3 r n x ω n x ω n 2 θ n x 2 d x , I 25 = c 0 + 2 c d c v L 2 0 1 ρ n ( r n ) 4 ( ω n x ) 2 2 θ n x 2 d x , I 26 = 4 c d c v L 0 1 ( ω n ) 2 r n x 2 θ n x 2 d x , I 27 = 8 c d c v L 0 1 r n ω n ω n x 2 θ n x 2 d x , I 28 = 4 μ r c v 0 1 ( ω n ) 2 ρ n 2 θ n x 2 d x .

Taking into account (92)-(94) and (95)-(100), we estimate the terms on the right-hand side of (102). For instance,

| I 1 | = 2 ( λ + 2 μ ) L 2 | 0 1 ( r n ) 3 ρ n x r n x v n 2 v n x 2 d x | C max x [ 0 , 1 ] | v n ( x , t ) | ρ n x ( t ) 2 v n x 2 ( t ) C v n x ( t ) ρ n x ( t ) 2 v n x 2 ( t ) .

Applying the Young inequality, we get

| I 1 | ε 2 v n x 2 ( t ) 2 + C ( v n x ( t ) 4 + ρ n x ( t ) 4 ) ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) ,

where ε>0 is arbitrary. In the analogous way, we obtain the following inequalities:

| I 2 | ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 3 | ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 4 | ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 5 | ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) d τ ) 8 ) , | I 6 | ε 2 v n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 + v n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 7 | ε 2 v n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 ) , | I 8 | ε 2 ω n x 2 ( t ) 2 + C , | I 9 | ε 2 ω n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 10 | ε 2 ω n x 2 ( t ) 2 + C , | I 11 | ε 2 ω n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 12 | ε 2 ω n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 13 | ε 2 ω n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 14 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 ) , | I 15 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 + ( 0 t 2 v n x 2 ( τ ) 2 d τ ) 8 ) , | I 16 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 + v n x ( t ) 16 ) , | I 17 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + θ n x ( t ) 16 + v n x ( t ) 16 ) , | I 18 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 19 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 20 | ε 2 θ n x 2 ( t ) 2 + ε 2 v n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 21 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 22 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + v n x ( t ) 16 ) , | I 23 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 24 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 25 | ε 2 θ n x 2 ( t ) 2 + ε 2 ω n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 26 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 27 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) , | I 28 | ε 2 θ n x 2 ( t ) 2 + C ( 1 + ω n x ( t ) 16 ) .

Using these inequalities with sufficiently small ε and estimates (92)-(94), from (102) we get (101). □

Lemma 4.7 There exists T 0 , (0< T 0 T) such that for each nN the Cauchy problem (79)-(86) has a unique solution defined on [0, T 0 ]. Moreover, the functions v n , ω n , θ n , ρ n and r n satisfy the inequalities

(103)
a 2 r n (x,t)2M, a 1 2 r n x (x,t)2 M 1 ,
(104)
m 2 ρ n (x,t)2M,(x,t) Q ¯ 0 , Q 0 = Q T 0 ,
(105)
max t [ 0 , T 0 ] ρ n x ( t ) C,
(106)
max t [ 0 , T 0 ] 2 r n x 2 ( t ) C
(107)

(a, a 1 , m and M are defined by (41) and (53)-(55)).

Proof To get the estimate (103) we use an approach similar to that in [3] (Lemma 5.6) and [14] (pp.64-67). First, we introduce the function

y n (t)= v n x ( t ) 2 + ω n x ( t ) 2 + θ n x ( t ) 2 +B 0 t 2 v n x 2 ( τ ) 2 dτ.
(108)

Using Lemma 4.6, we find that the function y n satisfies the differential inequality

y ˙ n (t)C ( 1 + y n 8 ( t ) ) .
(109)

There exists a constant C ¯ such that

C ¯ = d v 0 d x 2 + d ω 0 d x 2 + d θ 0 d x 2 ,

and we can conclude that

y n (0)= d v 0 n d x 2 + d ω 0 n d x 2 + d θ 0 n d x 2 C ¯ .
(110)

We compare the solution of the problem (109)-(110) with the solution of the Cauchy problem

y ˙ (t)=C ( 1 + y 8 ( t ) ) ,
(111)
y(0)= C ¯ .
(112)

Let [0, T [, 0< T T be an existence interval of the solution of the problem (111)-(112). From (109)-(112) we conclude that

y n (t)y(t),t[0, T [.
(113)

Let T 0 be such that 0< T 0 T . From (113) and (108) we obtain

max t [ 0 , T 0 ] ( v n x ( t ) 2 + ω n x ( t ) 2 + θ n x ( t ) 2 ) + 0 T 0 2 v n x 2 ( τ ) 2 dτC
(114)

and from (101) it follows

d d t ( v n x ( t ) 2 + ω n x ( t ) 2 + θ n x ( t ) 2 ) + B ( 2 v n x 2 ( t ) 2 + 2 ω n x 2 ( t ) 2 + 2 θ n x 2 ( t ) 2 ) C .
(115)

Integrating (101) over [0,t], 0<t T 0 and using estimates (110) and (115), we immediately get (103).

Now, using the inequalities (97) for the function v, we easily get

| v n (x,t)|2 v n x ( t ) , v n x ( t ) 2 2 v n x 2 ( t ) ,| v n x (t)|2 2 v n x 2 ( t ) .
(116)

Using (116), we derive the following estimates:

0 T 0 | v n (x,t)|dτ4 0 T 0 2 v n x 2 ( t ) dτ4 ( C B 1 ) 1 2 T 0 1 2 ,
(117)
0 T 0 | v n x (t)|dτ2 0 T 0 2 v n x 2 ( t ) dτ2 ( C B 1 ) 1 2 T 0 1 2 ,
(118)

where C and B are from (103). Assuming that

T 0 <min { T , a 2 B 64 C , a 1 2 B 16 C , ( L B 1 2 16 M 2 ( 4 M 1 + M ) C 1 2 ) 2 }

and using (117) and (118) from (58) and (60), we get (104)-(105).

Because of (57), (62) and (63), from (103) and (98), we easily get that for t[0, T 0 ] hold

i = 1 n ( | v i n ( t ) | + | ω i n ( t ) | + | θ i n ( t ) | ) C,
(119)
| θ 0 n (t)|C,
(120)

and we can conclude that the solution of the problem (79)-(86) is defined on [0, T 0 ]. □

From (119) and (120), we can easily conclude that

max t [ 0 , T 0 ] ( v n ( t ) 2 + ω n ( t ) 2 + θ n ( t ) 2 ) C,
(121)

and from (95), (100) and (99) it follows

max t [ 0 , T 0 ] 2 r n x 2 ( t ) C, max t [ 0 , T 0 ] 2 ρ n x 2 ( t ) C, max ( x , t ) Q ¯ 0 | θ n (t)|C.
(122)

Lemma 4.8 Let T 0 be defined by Lemma  4.7. Then for each nN the inequalities

0 T 0 ( v n t ( τ ) 2 + ω n t ( τ ) 2 + θ n t ( τ ) 2 ) dτC,
(123)
max t [ 0 , T 0 ] ρ n t ( t ) C,
(124)
max t [ 0 , T 0 ] r n t ( t ) C, max t [ 0 , T 0 ] 2 r n x t ( t ) C, 0 T 0 2 r n t 2 ( τ ) 2 dτC
(125)

hold true.

Proof Multiplying (65) by d v i n d t (t), summing over i=1,2,,n and using (104)-(105), we obtain

v n t ( t ) 2 = R L 0 1 ( r n ) 2 x ( ρ n θ n ) v n t d x + λ + 2 μ L 2 0 1 ( r n ) 2 x ( ρ n x ( ( r n ) 2 v n ) ) v n t d x C ( max ( x , t ) Q ¯ 0 | θ n ( x , t ) | ρ n x ( t ) v n t ( t ) + θ n x ( t ) v n t ( t ) + max ( x , t ) Q ¯ 0 | v n ( x , t ) | ρ n x ( t ) v n t ( t ) + v n ( t ) v n t ( t ) + max ( x , t ) Q ¯ 0 | v n ( x , t ) | 2 r n x 2 ( t ) v n t ( t ) + v n t ( t ) v n x ( t ) + max ( x , t ) Q ¯ 0 | v n x ( x , t ) | ρ n x ( t ) v n t ( t ) + 2 v n x 2 ( t ) v n x ( t ) ) .
(126)

Using (121), (122), (103), (116) and applying the Young inequality, we get

v n t ( t ) 2 C ( 1 + 2 v n x 2 2 ) +ε v n t ( t ) 2 .
(127)

Taking into account (103) for sufficiently small ε>0 from (127), we obtain

0 T 0 v n t ( τ ) 2 dτC.
(128)

In the same way, from (66) and (67) we obtain the estimates for ω n t and θ n t respectively. The estimates (124) and (125) follow directly from (59) and (58). □

From Lemmas 4.7 and 4.8 we easily derive the following statements.

Proposition 4.1 Let T 0 be defined by Lemma  4.7. Then for the sequence {( r n , ρ n , v n , ω n , θ n ):nN} the following properties are satisfied:

  1. (i)

    { r n } is bounded in L ( Q 0 ), L (0, T 0 ; H 2 (]0,1[)) and H 2 ( Q 0 );

  2. (ii)

    { r n x } is bounded in L ( Q 0 );

  3. (iii)

    { ρ n } is bounded in L ( Q 0 ), L (0, T 0 ; H 1 (]0,1[)) and H 1 ( Q 0 );

  4. (iv)

    { v n }, { ω n }, { θ n } are bounded in L (0, T 0 ; H 1 (]0,1[)), H 1 ( Q 0 ) and L 2 (0, T 0 ; H 2 (]0,1[)).

5 The proof of Theorem 2.1

In the proofs we use some well-known facts of functional analysis (e.g., [21]). Let T 0 R + be defined by Lemma 4.7. Theorem 2.1 is a consequence of the following lemmas.

Lemma 5.1 There exists a function

r L ( 0 , T 0 ; H 2 ( ] 0 , 1 [ ) ) H 2 ( Q 0 )C( Q ¯ 0 )
(129)

and the subsequence of { r n } (for simplicity reasons denoted again as { r n }) such that

r n rweakly- in  L ( 0 , T 0 ; H 2 ( ] 0 , 1 [ ) ) ,
(130)
r n rweakly in  H 2 ( Q 0 ),
(131)
r n rstrongly in C( Q ¯ 0 ),
(132)
r n x r x strongly in C( Q ¯ 0 ).
(133)

The function r satisfies the conditions

a 2 r2Min  Q ¯ 0 ,
(134)
r(x,0)= r 0 (x),x[0,1],
(135)

where r 0 is defined by (40).

Proof The conclusions (130) and (131) follow immediately from Proposition 4.1. Let (x,t), ( x , t ) belong to Q ¯ 0 . Then we have

| r n ( x , t ) r n ( x , t ) | | r n ( x , t ) r n ( x , t ) | + | r n ( x , t ) r n ( x , t ) | , | r n x ( x , t ) r n x ( x , t ) | | r n x ( x , t ) r n x ( x , t ) | + | r n x ( x , t ) r n x ( x , t ) | .

Using (104), (58), (116), (103) and (107), we obtain

| r n ( x , t ) r n ( x , t ) | x x | r n x ( ξ , t ) | dξC | x x | ,
(136)
(137)
| r n x ( x , t ) r n x ( x , t ) | 2 r n x 2 ( t ) | x x | 1 / 2 C | x x | 1 / 2 ,
(138)
(139)

and we can conclude that the sequences { r n } and { r n x } satisfy the conditions of Arzelà-Ascoli theorem. Applying that theorem, we get the strong convergence in (132) and (133). Because of (132) and (104) we have

a 2 ε< r n (x,t)ε<r(x,t)< r n (x,t)+ε<2M+ε,(x,t) Q ¯ 0
(140)

for each ε>0 and sufficiently big nN. From (140) we can easily conclude that (134) is satisfied. From (132) it follows

lim n max x [ 0 , 1 ] | r n (x,0)r(x,0)|= lim n max x [ 0 , 1 ] | r 0 (x)r(x,0)|=0
(141)

and because of that we have (135). □

Lemma 5.2 There exists a function

ρ L ( 0 , T 0 ; H 1 ( ] 0 , 1 [ ) ) H 1 ( Q 0 )C( Q ¯ 0 )
(142)

and the subsequence of { ρ n } (denoted again as { ρ n }) such that

ρ n ρ weakly- in  L ( 0 , T 0 ; H 1 ( ] 0 , 1 [ ) ) ,
(143)
ρ n ρweakly in  H 1 ( Q 0 ),
(144)
ρ n ρstrongly in C( Q ¯ 0 ).
(145)

The function ρ satisfies the conditions

m 2 ρ(x,t)2Min  Q ¯ 0 ,
(146)
ρ(x,0)= ρ 0 (x),x[0,1].
(147)

Proof Taking into account Proposition 4.1, estimates (103)-(106) and the Arzelà-Ascoli theorem, we prove in the same way as in the previous lemma the properties (143)-(147). □

Lemma 5.3 There exist functions

v,ω,θ L ( 0 , T 0 ; H 1 ( ] 0 , 1 [ ) ) H 1 ( Q 0 ) L 2 ( 0 , T 0 ; H 2 ( ] 0 , 1 [ ) )

and the subsequence of { v n , ω n , θ n } (denoted again as { v n , ω n , θ n }) such that

( v n , ω n , θ n ) (v,ω,θ)weakly- in  ( L ( 0 , T 0 ; H 1 ( ] 0 , 1 [ ) ) ) 3 ,
(148)
( v n , ω n , θ n ) (v,ω,θ)weakly in  ( H 1 ( Q 0 ) ) 3 ,
(149)
( v n , ω n , θ n ) (v,ω,θ)weakly in  ( L 2 ( 0 , T 0 ; H 2 ( ] 0 , 1 [ ) ) ) 3 ,
(150)
( v n , ω n , θ n ) (v,ω,θ)strongly in  ( L 2 ( Q 0 ) ) 3 .
(151)

The functions v, ω and θ satisfy the conditions

v(0,t)=v(1,t)=ω(0,t)=ω(1,t)=0,t[0, T 0 ],
(152)
θ x (0,t)= θ x (1,t)=0,a.e. in ]0, T 0 [,
(153)
v(x,0)= v 0 (x),ω(x,0)= ω 0 (x)=0,θ(x,0)= θ 0 (x),x[0,1].
(154)

Proof The conclusions (148)-(151) follow from Proposition 4.1 and embedding properties (see Remark 2.1). For verification of the boundary and initial conditions (152), (153) and (154), we use the Green formula as follows.

Let φ be a function of C ([0, T 0 ]) equal to zero in a neighborhood of T 0 , with φ(0)0 and u H 1 (]0,1[). Using the integration by parts we have for v n and v (e.g.) the following equalities:

0 T 0 0 1 v n t ( x , t ) u ( x ) φ ( t ) d x d t + 0 T 0 0 1 v n ( x , t ) u ( x ) d φ d t ( t ) d x d t = φ ( 0 ) 0 1 v 0 n ( x ) u ( x ) d x ,
(155)
0 T 0 0 1 v t ( x , t ) u ( x ) φ ( t ) d x d t + 0 T 0 0 1 v ( x , t ) u ( x ) d φ d t ( t ) d x d t = φ ( 0 ) 0 1 v ( x , 0 ) u ( x ) d x .
(156)

Passing n in (155) and comparing (155) and (156), we obtain

0 1 v(x,0)u(x)dx= 0 1 v 0 (x)u(x)dx,u H 1 ( ] 0 , 1 [ )
(157)

and conclude

v(x,0)= v 0 (x),x[0,1].
(158)

In the similar way, we get all the remaining properties in (152)-(154). □

Lemma 5.4 The functions r, ρ, v, ω, θ, defined by Lemmas  5.1, 5.2 and  5.3 satisfy the Equations (30)-(33) a.e. in Q 0 .

Proof Let {( r n , ρ n , v n , ω n , θ n ):nN} be subsequence defined by Lemmas 5.1, 5.2 and 5.3. Taking into account (144), (149) and strong convergencies (132), (133), (145) and (151) we get that (30) follows immediately from (59). We can write Equation (65) in the form

0 T 0 0 1 v n t ( x , t ) sin ( π i x ) φ ( t ) d x d t R L 0 T 0 0 1 [ ( r n ) 2 ρ n x θ n + ( r n ) 2 ρ n θ n x ] sin ( π i x ) φ ( t ) d x d t λ + 2 μ L 2 0 T 0 0 1 [ 2 ( r n ) 3 ρ n x r n x v n + 2 ( r n ) 2 ρ n ( r n x ) 2 v n + 2 ( r n ) 3 ρ n 2 r n x 2 v n + 4 ( r n ) 3 ρ n r n x v n x + ( r n ) 4 ρ n x v n x + ( r n ) 4 ρ n 2 v n x 2 ] sin ( π i x ) φ ( t ) d x d t = 0
(159)

for each φD(]0, T 0 [), where D denotes the space of test functions. Now we consider the convergence of one integrand. For example, we will prove the convergence

0 T 0 0 1 ( r n ) 3 ρ n 2 r n x 2 v n sin(πix)φ(t)dxdt 0 T 0 r 3 ρ 2 r x 2 vsin(πix)φ(t)dxdt
(160)

when n. Taking into account

| 0 T 0 0 1 [ ( r n ) 3 ρ n 2 r n x 2 v n r 3 ρ 2 r x 2 v ] sin ( π i x ) φ ( t ) d x d t | C max ( x , t ) Q ¯ 0 | ( r n ) 3 r 3 | 2 r n x 2 ( t ) v n ( t ) + C max ( x , t ) Q ¯ 0 | ρ n ρ | 2 r n x 2 ( t ) v n ( t ) + C v n ( t ) | 0 T 0 0 1 ( 2 r n x 2 2 r x 2 ) sin ( π i x ) φ ( t ) d x d t | + C 2 r ( t ) x 2 ( v n v ) ( t )
(161)

and using already mentioned convergences, we can easily conclude that (160) is satisfied. In the same way, we can derive the convergences of other integrals in (159). Analogously, we get that (32) and (33) follow from (66) and (67). □

Remark 5.1 Taking into account (150) and (58), we can easily prove that the function r defined by Lemma 5.1 has the form

r(x,t)= r 0 (x)+ 0 t v(x,t)dτ,(x,t) Q ¯ 0 ,
(162)

where v is from Lemma 5.3.

Lemma 5.5 There exists T 0 , 0< T 0 T such that the function θ defined by Lemma  5.3 satisfies the condition

θ>0in  Q ¯ 0 .
(163)

Proof Because of the inclusion θC( Q ¯ 0 ) (see Remark 2.1), in the same way as in [3], we conclude that for each ε>0 there exists T 0 , T 0 T, such that for (x,t) Q ¯ 0 holds

| θ ( x , t ) θ ( x , 0 ) | = | θ ( x , t ) θ 0 ( x ) | < ε , θ ( x , t ) > θ 0 ( x ) ε m ε .

The conclusions of Theorem 2.1 are an immediate consequence of the above lemmas. □

endnote

We assume Einstein notation for summation.

References

  1. Eringen CA: Simple microfluids. Int. J. Eng. Sci. 1964, 2(2):205-217. 10.1016/0020-7225(64)90005-9

    Article  MATH  MathSciNet  Google Scholar 

  2. Lukaszewicz G Modeling and Simulation in Science, Engineering and Technology. In Micropolar Fluids: Theory and Applications. Birkhäuser, Boston; 1999.

    Chapter  Google Scholar 

  3. Mujakovic N: One-dimensional flow of a compressible viscous micropolar fluid: a local existence theorem. Glas. Mat. 1998, 33: 71-91.

    MATH  MathSciNet  Google Scholar 

  4. Mujakovic N: Nonhomogeneous boundary value problem for one-dimensional compressible viscous micropolar fluid model: regularity of the solution. Bound. Value Probl. 2008., 2008: Article ID 189748

    Google Scholar 

  5. Mujakovic N: One-dimensional compressible viscous micropolar fluid model: stabilization of the solution for the Cauchy problem. Bound. Value Probl. 2010., 2010: Article ID 796065

    Google Scholar 

  6. Easwaran C, Majumdar S: A uniqueness theorem for compressible micropolar flows. Acta Mech. 1987, 68: 185-191. 10.1007/BF01190882

    Article  MATH  MathSciNet  Google Scholar 

  7. Chen M: Blowup criterion for viscous, compressible micropolar fluids with vacuum. Nonlinear Anal., Real World Appl. 2012, 13(2):850-859. 10.1016/j.nonrwa.2011.08.021

    Article  MATH  MathSciNet  Google Scholar 

  8. Hoff D: Spherically symmetric solutions of the Navier-Stokes equations for compressible, isothermal flow with large, discontinuous initial data. Indiana Univ. Math. J. 1992, 41: 1225-1302. 10.1512/iumj.1992.41.41060

    Article  MATH  MathSciNet  Google Scholar 

  9. Jiang S: Global spherically symmetric solutions to the equations of a viscous polytropic ideal gas in an exterior domain. Commun. Math. Phys. 1996, 178: 339-374. 10.1007/BF02099452

    Article  MATH  Google Scholar 

  10. Jiang S, Zhang P: On spherically symmetric solutions of the compressible isentropic Navier-Stokes equations. Commun. Math. Phys. 2001, 215: 559-581. 10.1007/PL00005543

    Article  MATH  Google Scholar 

  11. Fujita-Yashima H, Benabidallah R: Equation à symétrie sphérique d’un gaz visqueux et calorifère avec la surface libre. Ann. Mat. Pura Appl. 1995, 168: 75-117. 10.1007/BF01759255

    Article  MATH  MathSciNet  Google Scholar 

  12. Yanagi S: Asymptotic stability of the spherically symmetric solutions for a viscous polytropic gas in a field of external forces. Transp. Theory Stat. Phys. 2000, 29(3-5):333-353. 10.1080/00411450008205878

    Article  MathSciNet  Google Scholar 

  13. Mujakovic N: Non-homogeneous boundary value problem for one-dimensional compressible viscous micropolar fluid model: a local existence theorem. Ann. Univ. Ferrara 2007, 53: 361-379. 10.1007/s11565-007-0023-z

    Article  MATH  MathSciNet  Google Scholar 

  14. Antontsev SN, Kazhikhov AV, Monakhov VN Studies in Mathematics and its Applications 22. In Boundary Value Problems in Mechanics of Nonhomogeneous Fluids. North-Holland, Amsterdam; 1990.

    Google Scholar 

  15. Weigant V: Example of non-existence in the large for the problem of the existence of solutions of Navier-Stokes equations for compressible viscous barotropic fluids. Dokl. Akad. Nauk SSSR 1994, 339: 155-156.

    Google Scholar 

  16. Lions JL, Dautray R Mathematical Analysis and Numerical Methods for Science and Technology 2. In Functional and Variational Methods. Springer, Berlin; 2000.

    Google Scholar 

  17. Lions JL, Dautray R Mathematical Analysis and Numerical Methods for Science and Technology 5. In Evolution Problems I. Springer, Berlin; 2000.

    Google Scholar 

  18. Kazhikhov AV: To a theory of boundary value problems for equation of one-dimensional non-stationary motion of viscous heat-conductive gases. Bound. Value Probl. Hydrodyn. Equ. 1981, 50: 327-362.

    Google Scholar 

  19. Arnold VI: Ordinary Differential Equations. MIT Press, Cambridge; 1978.

    Google Scholar 

  20. Petrowski IG: Vorlesungen über die Theorie der gewöhnlichen Differentialgleichungen. Teubner, Leipzig; 1954.

    MATH  Google Scholar 

  21. Brezis H: Analyse fonctionnelle: théorie et applications. Masson, Paris; 1983.

    MATH  Google Scholar 

Download references

Acknowledgement

The paper was made with financial support of the scientific project Mathematical analysis of composite and thin structures (037-0693014-2765), Ministry of Science, Education and Sports, Republic of Croatia.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ivan Dražić.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The paper is the result of joint work of both authors who contributed equally to the final version of the paper. 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

Cite this article

Dražić, I., Mujaković, N. 3-D flow of a compressible viscous micropolar fluid with spherical symmetry: a local existence theorem. Bound Value Probl 2012, 69 (2012). https://doi.org/10.1186/1687-2770-2012-69

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-69

Keywords