arxiv: v1 [math.fa] 30 Jun 2011

Transcription

1 Existence of strong solutions for the compressible arxiv: v1 [math.fa] 3 Jun 211 Ericksen-Leslie model Xiangao Liu, Lanming Liu, Yihang Hao School of Mathematic Sciences, Fudan University, Shanghai, 2433, P.R.China Abstract: In this paper, we prove the existence and uniqueness of local strong solutions of the hydrodynamics of nematic liquid crystals system under the initial data satisfying a natural compatibility condition. Also the global strong solutions of the system with small initial data are obtained. Keywords: Compressible liquid crystals; Strong solutions; Vacuum. Mathematics Subject Classification 35M1, 76N1, 82D3. 1 Introduction It is well known liquid crystals are states of matter which are capable of flow and in which the molecular arrangement gives rise to a preferred direction. By Ericksen-Leslie theory, the compressible nematic liquid crystals reads the following system: ρ t +div(ρu) =, (1.1) (ρu) t +div(ρu u)+ p = µ u λdiv( d d 1 2 ( d 2 +F(d))I), (1.2) d t +u d = ν( d f(d)), (1.3) xgliu@fudan.edu.cn lanmingliu@gmail.com @fudan.edu.cn 1

2 2 in (,T), for a bounded smooth domain R 3. ρ(t,x) is density, u(t,x) the velocity field, d(t, x) orientation parameter of the liquid crystal and p(ρ) pressure with p = p( ) C 1 [, ), p() =. (1.4) The viscosity coefficients µ,λ,ν are positive constants. The unusual term d d denotes the 3 3 matrix whose (i,j)-th element is given by 3 k=1 x i d k xj d k. I is the unite matrix. f(d) = 1 σ 2( d 2 1)d and F(d) = 1 4σ 2( d 2 1) 2. We are interested in the initial data ρ(,x) = ρ, u(,x) = u, d(,x) = d (x), x, (1.5) and the boundary condition u(t,x) =, d(t,x) = d (x), d (x) = 1, (t,x) (,T). (1.6) The first author[1] firstly considers the model(1.1)-(1.4) and gives a global existence of finite energy weak solutions by using Lions s technique(see[2,3]). The incompressible model of a similar simplified Ericksen-Leslie model has been studied in the papers [4 8]. Recently, for the incompressible Ericksen-Leslie model, Wen-Ding [9] have gave the proof of local strong solutions in two and three dimensions, and Lin-Lin-Wang [1] have established the existence of global (in time) weak solutions on a bounded smooth domain in two dimensions. Throughout this paper, we use the following simplified notations L q = L q (), W k,q = W k,q (), H k = W k,2 (), H 1 = W1,2 (). (1.7) In order to obtain strong solutions, we need (ρ,u,d ) satisfies the regularity ρ W 1,6, u H 1 H2, d H 3, (1.8)

3 3 and a compatibility condition µ u λdiv( d d 1 2 ( d 2 +F(d ))I) p = ρ 1 2 g, x, (1.9) for some g L 2. The compatibility condition is a natural request to insure ρu t () L 2 bounded, and a compensation to the lack of a positive lower bound of ρ. It is firstly introduced by Salvi, Străskraba [11] and the authors of [12] independently. Our main result is the following theorem: Theorem 1. Part I: (Local existence) Under the assumptions of the regularity condition (1.8) and the compatibility condition (1.9), there exists a small T (,T) such that the system (1.1)-(1.6) has a unique local strong solution (ρ, u, d) satisfying ρ C([,T ];W 1,6 ), ρ t C([,T ];L 6 ), u C([,T ];H 1 H2 ) L 2 (,T ;W 2,6 ), d C([,T ];H 3 ), u t L 2 (,T ;H 1), d t C([,T ];H 1) L2 (,T;H 2 ), ρut C([,T ];L 2 ), d tt L 2 (,T ;L 2 ). (1.1) Part II: (Continuity of initial data) Suppose that ( ρ,ũ, d) is another solution of (1.1)-(1.4) with the following initial boundary conditions ρ(,x) = ρ, ũ(,x) = ũ, d(,x) = d (x), x, ũ(t,x) =, d(t,x) = d (x), (t,x) (,T), then for any t (,T ], the following quantities tend to zero when ( ρ,ũ, d ) (ρ,u,d ) in W 1,6 H 2 H 3 : d d H 2(t), d d L 2 (,T ;H 3 ), d t d t L 2(t), d t d t L 2 (,T ;H 1 ), u ũ H 1(t), u ũ L 2 (,T ;H 2 ), ρ(u t ũ t ) L 2 (,T ;L 2 ), ρ ρ L 6(t). In general, as Navier-Stokes equations we have not the global existence with large initial data. However, a blow-up criterion is obtained in our article [13] and for small initial data, we have the following global existence.

4 4 Theorem 2. (Global existence with small initial data) Let α be a nonnegative constant and m a constant unite vector in R 3. Then there exists a positive constant θ small such that if the initial data satisfies max{ ρ α W 1,6, u H 2, d m H 3, g 2 L2} θ, (1.11) for all θ (,θ ], then the system(1.1)-(1.6) has a unique global strong solution. Remark 1. In this paper, we only consider the case α =, because the other case (α > ) can be induced to the problem with a positive initial data of density. The methods to prove the Theorem 1 and Theorem 2 are the successive approximation in Sobolev spaces (see [14, 15]). We generalize this method to two variables, which is based on the following careful observation on coupling terms of d and u: (1) d 2 H 2 can be deduced by u 2 H 2 dτ and u 2 H 1 ; (2) u 2 H 2 can be derived by d 2 H 2 dτ. Since the initial density has vacuum, the movement equation(1.2) becomes a degenerate parabolic-elliptic couples system. To overcome this difficulty, as usual, the technique is to approximate the nonnegative initial data of density by a positive initial data. For a special linear equations of the system (1.1)-(1.3), we prove not only the local existence of strong solutions with large initial data but also the global existence of strong solutions with small initial data. Employing energy law and higher energy inequalities, we can prove both the uniqueness and the continuous dependence on the initial data. This paper is written as follows. In section 2, after establishing a linear problem of the nonlinear problem (1.1)-(1.3), we prove local existence of a strong solution to the special linear problem with a positive initial data of density. And we also establish some uniform a prior estimates, which imply the existence of a local strong solution to the linear problem when the initial data of density allows vacuum in section 2.3. In section 3, after constructing a sequence of approximate solutions, a strong solution of the nonlinear problem (1.1)-(1.3) is obtained. In section 4, the uniqueness and continuous

5 5 dependence on initial data are proved. 2 A linear problem 2.1 Linearization At the beginning, we linearize the equations (1.1)-(1.3) as following: ρ t +div(ρv) =, (2.1) d t +v d = ν( d 1 σ2(n+m) (d m)n), (2.2) ρu t +ρv u+ p(ρ) = µ u λ( d) T ( d f(d)), (2.3) where m is a constant unite vector and (v, n) satisfies the following regularity v C([,T];H 1 H2 ) L 2 (,T;W 2,6 ), v t L 2 (,T;H 1 ), n C([,T];H 2 ) L 2 (,T;H 3 ), n t C([,T];H 1 ) L2 (,T;H 2 ). Assume further that v(,x) = u (x), n(,x) = d (x), x, v(t,x) =, n(t,x) = d (x), (t,x) (,T). (2.4) (2.5) 2.2 The existence of approximate solutions For each δ (,1), let u δ be the solution to the boundary value problem: µ u δ λdiv( d d 1 2 ( d 2 +F(d ))I) p(ρ δ ) = (ρδ )1 2g, u δ =, x, (2.6) where ρ δ = ρ +δ. Then u δ u in H 1 H2 as δ. For the linear equations (2.1)-(2.3), we have the following theorem. Theorem 2.1. Assume ρ, u and d satisfy the regularity (1.8) and (v,n) satisfies the above conditions (2.4) and (2.5). Then there exists a unique strong solution (ρ, u, d)

6 6 of the linear equations (2.1)-(2.3) with the initial data (ρ δ,uδ,d ) and the boundary condition (1.6) such that ρ C([,T];W 1,6 ), ρ t C([,T];L 6 ), u C([,T];H 1 H2 ) L 2 (,T;W 2,6 ), u t C([,T];L 2 ) L 2 (,T;H 1 ), d C([,T];H 2 ) L 2 (,T;H 3 ), d t C([,T];H 1 ) L 2 (,T;H 2 ), (2.7) u tt L 2 (,T;H 1 ), d tt L 2 (,T;L 2 ). We will prove the theorem by the following three lemmas. Firstly, suppose the constants c, c 1 and c 2 satisfying c > 1+ ρ W 1,6 + d H 3 + u H 2 + g L 2, (2.8) c 1 > sup ( v H 1 + n H 1 + n t H 1 )+ ( v t 2 L 2 + v 2 W 2,6 t T + 2 n t 2 L 2 + n 2 H3)dt, (2.9) c 2 > sup ( 2 v L n L 2), (2.1) t T c 2 > c 1 > c > 1. (2.11) T Lemma 2.1. Assume ρ and v satisfy the regularities (1.8) and (2.4) respectively. Then the problem (2.1) and (1.5) has a global unique strong solution such that ρ C([,T];W 1,6 ), ρ t C([,T];L 6 ). (2.12) Moreover, if v,n satisfy (2.8)-(2.11), then there is a T 1 = min{c 1 1,T} such that t [,T 1 ], ρ W 1,6(t) Cc, ρ t L 6(t) Cc c 2. (2.13) In particular, p L 6(t) CM(c ), p L 6(t) CM(c )c, p t L 6(t) CM(c )c c 2, (2.14) where the constant M(c ) is defined by (2.18).

7 7 Remark 2. The above lemma is also true to the problem (2.1) with the initial data ρ δ and the estimates also hold for all small δ. Proof. Now let s start to prove Lemma 2.1. The existence is obvious, based on the classical method of characteristics (see [15]). So ρ(t,x) = ρ (y(,t,x))exp( divv(τ, y(τ, t, x))dτ), (2.15) where y(τ,t,x)( C([,T] [,T] )) is a solution of the initial value problem: Moreover, we have y(τ,t,x) = v(τ,y(τ,t,x)), τ T, τ y(t,t,x) = x, t T, x. ρ W 1,6(t) ρ W 1,6exp( v W 1,6(τ)dτ), t [,T]. (2.16) So that for all t [,T 1 ], ρ W 1,6(t) ρ W 1,6exp( v W 1,6(τ)dτ) c exp(t 1 2 ( v 2 W 1,6(τ)dτ)1 2 ) ec. (2.17) Because of ρ L C ρ W 1,6, set We can obtain for all t [,T 1 ], M(c ) = sup Cec (1+p( )+p ( )). (2.18) ρ t L 6(t) = v ρ ρdivv L 6 v L ρ L 6 + ρ L v L 6 v H 2 ρ W 1,6 + ρ W 1,6 v H 2 Cc c 2.

8 8 Similarly, for all t [,T 1 ], p L 6(t) CM(c ), p L 6(t) CM(c )c, p t L 6(t) CM(c )c 2 c. Lemma 2.2. Under the hypotheses of Theorem 2.1, then the equation (2.2) with the initial boundary conditions (1.5)-(1.6) has a global unique strong solution d satisfying (2.7). Moreover, if v,n satisfy (2.8)-(2.11), then there is a T 3 = min{c 22 2,T} such that sup ( d H 1 +c d L 2 +c 3 2 d t H 1 +c 2 2 c 1 1 d H 2) t T 3 + T3 d t 2 H 2 +c3 d 2 H 3dt Cc3. (2.19) Proof. By the classic Galerkin method to the linear parabolic equation (2.2) with (1.5)- (1.6), the existence and regularity of d described in (2.7) can be obtained. Differentiating (2.2) with respect to time, multiplying by d t and then integrating over, we can deduce that 1 d d t 2 dx+ν d t 2 dx 2dt C( v t L 6 d L 2 d t L 3 + v L d t L 2 d t L 2 + n L 3 n t L 6 d m L 3 d t L 6 + n+m L 3 n t L 6 d m L 3 d t L 6 + n L 6 n+m L 6 d t 2 L 3) 5 = I i. Here i=1 (2.2) I 1 Cη v t 2 H 1 d 2 L 2 +Cη 1 d t 2 L + ν 2 5 d t 2 L 2, I 2 C v 2 H 2 d t 2 L + ν 2 5 d t 2 L 2, I 3 +I 4 C(( n 2 L 2 +1)( n 2 H +1) n t 4 1 L 6 d m 2 L 2 + d 2 L 2 + d m 2 L 2) + ν 5 d t 2 L 2, I 5 C n 2 H 1 n+m 2 H 1 d t 2 L + ν 2 5 d t 2 L 2,

9 9 where the small positive constant η will be fixed later. Since d dt d dt d 2 dx = 2 d m 2 dx combining (2.2), (2.21) and (2.22), we get d d t 2 + d 2 + d m 2 dx+ dt d : d t dx d 2 dx+ d t 2 dx, (2.21) d m 2 dx+ d t 2 dx, (2.22) d t 2 dx C( d 2 L 2 + d m 2 L 2 + d t 2 L 2) (η 1 + v 2 H 2 + n 2 H 1 n+m 2 H 1 +η v t 2 H 1 +( n 2 L 2 +1)( n 2 H +1) n t 4 1 L2 +1). (2.23) From the equation (2.2), we can deduce d t L 2() C( d L 2 + u H 2 d L 2 + d +m L 2 d m L 2 d L 2) Cc 3. (2.24) Hence, by Gronwall s inequality, we can deduce from (2.23) d t 2 + d 2 + d m 2 dx+ d t 2 dxdτ Cc 6 exp(c η 1 + v 2 H 2 + n 2 H 1 + n+m 2 H +η v t 2 1 H 1 +( n 2 L 2 +1)( n 2 H +1) n t 4 1 L2dτ). (2.25) Taking η = c 1 1 and using the assumption (2.8)-(2.11), we obtain T2 d t 2 + d 2 + d m 2 dx+ d t 2 dxdτ Cc 6, (2.26) sup t T 2 where T 2 = min{c 8 2,T 1 }. From the equation (2.2) and using the elliptic estimates, we get d m H 2 C( d m H 2 + d t L 2 + v d L 2 + (n+m) (d m)n L 2) C(c 3 + v H 1 d 1 2 L 2 d m 1 2 H 2 + n+m L 6 d m L 6 n L 6) C(c 3 + v 2 H 1 d L 2 + n+m H 1 d m H 1 n H 1)+ 1 2 d H 2.

10 1 So 2 d L 2 C(c 3 +c 2 1c 3 +c 2 1c 3 ) Cc 2 1c 3. (2.27) Applying the operator to the linear equation (2.2), we get ν ( d) = d t + (v d)+ ν ((n+m) (d m)n). (2.28) σ2 By the elliptic estimates, we can estimate the term (d m) H 2 as follows (d m) H 2 C( d t L 2 + (v d) L 2 + ν σ 2 [(n+m) (d m)n] L 2 + d m H 3) C( d t L 2 + v H 1 d 1 2 L 2 d 1 2 H 1 + v H 2 2 d L 2 + n H 1 n H 1 d m H 1 + n 2 H 2 d L 2 + d m H 3) C( d t L 2 +c 2 2 c 1c 3 ), (2.29) where we use the assumption (2.8)-(2.11) and (2.26)-(2.27). Differentiating (2.2) with respect to time and taking inner product with d t, then we can derive 1 d d t 2 dx+ ν 2dt C d t 2 H 2 1 d d t 2 dx+ν d t 2 dx 2dt = (v t d) d t dx+ (v d t ) d t dx+ ν (n σ 2 t (d m))n d t dx + ν (n d σ 2 t )n d t dx+ ν ((n+m) (d m))n σ 2 t d t dx = 5 J j, j=1 where we use the elliptic estimate d t 2 H 2 C d t 2 L 2. (2.3)

11 11 Here J 1 = (v t d) d t dx v t L 2 d L 6 d t L 3 + v t L 3 2 d L 6 d t L 2 η v t 2 L 2 +η 1 d 2 H 1 d t L 2 d t H 1 +η v t 2 L 3 d t 2 L 2 +η 1 2 d 2 L 6 η v t 2 L 2 +Cη 2 ε 1 c 8 1 c12 d t 2 L 2 +ε d t 2 H 2 +η v t 2 H 1 d t 2 L 2 +Cη 1 ( d t 2 L 2 +c4 2 c2 1 c6 ), J 2 = (v d t ) d t dx Cε 1 c 2 2 d t 2 L +ε d t 2 2 L 2, J 4 Cε 1 n 4 L d t 2 L +ε d t 2 2 L 2 Cε 1 c 4 2c 6 +ε d t 2 L 2, J 3 + J 5 Cε 1 ( n 2 L +1) n t 2 L 2 d m 2 L +ε d t 2 L 2 Cε 1 c 2 2c 6 1c 6 +ε d t 2 L 2, where the small positive constants ε and η will be fixed later. Substituting the above J 1 J 5 into (2.3) and taking ε small enough, we get d d t 2 dx+ d t 2 H dt A η(t) d 2 t 2 L +B η(t), (2.31) 2 where A η (t) = C(η 2 c 8 1 c12 +η v t 2 H 1 +η 1 +c 2 2 ), B η (t) = C(η v t 2 H 1 +η 1 c 4 2 c2 1 c8 +η 1 c 2 2 c6 1 c3 +c4 2 c6 ). In view of the inequalities (2.26)-(2.27) and taking η = c 1 2, we have for all t [,T 2 ], A η(s)ds C +Cc 2 2c 8 1c 12 t, B η(s)ds C +Cc 5 2c 4 1c 6 t. (2.32) From the equation (2.28), we can estimate the term d t (τ) L 2 at time d t () L 2 C( d m H 3 + u L 3 2 d L 6 + u L 3 2 d L 6 + d L 6 d m L 6 d L 6 + d 2 L 6 d L 6) Cc 3. (2.33)

12 12 Applying Gronwall s inequality to (2.31) and using (2.32)-(2.33), we can obtain T3 sup d t 2 dx+ d t 2 H 2dt C(c3 +c 5 2c 4 1c 6 T 3 )exp(c +Cc 2 2c 8 1c 12 T 3 ) t T 3 Cc 3. (2.34) From (2.29) and (2.34), we have d H 2(t) C(c 3 2 +c 2 2c 1 c 3 ), t [,T 3 ]. (2.35) Hence d 2 H 3dτ C, t [,T 3]. (2.36) Since (2.26), (2.27), (2.34) and (2.36) can deduce the estimate (2.19), we complete the proof. Lemma 2.3. Under the hypotheses of Theorem 2.1, then there exists a global unique solution u of the equation (2.3) with the initial data u δ and the boundary condition (1.6), and the solution u has the regularity in (2.7). Moreover, if v, n satisfy (2.8)-(2.11), then, for all small δ, sup (M(c )c 13 u H 1 +M(c )c 8 c u L 2 +M(c )c 13 ρu t L 2) t T 3 + T3 c 8 u t 2 L 2 + u 2 W2,6dt Cc18 M2 (c ). (2.37) Proof. Because ρ δ δ >, it follows from the representation (2.15) that ρ(t,x) δexp( where δ is a positive constant. v(τ) W 1,6dτ) δ, (t,x) [,T], (2.38) Thanks to (2.38), we change (2.3) into the following form: u t +v u+ 1 ρ p = µ ρ u λ ρ ( d)t ( d f(d)).

13 13 Applying the Galerkin method again to the above equation with the initial data u δ and the boundary condition (1.6), we can deduce the existence and regularity of u described in (2.7). Differentiating (2.3) with respect to time t, multiplying by u t and then integrating over, we can derive Here = = 1 d ρ u t 2 dx+µ u t 2 dx 2dt ( p t ρ t v u 2ρv u t ρv t u))u t dx λ ( d t ) T ( d f(d))u t dx λ ( d) T ( d f(d)) t u t dx 6 K k. (2.39) k=1 K 1 C p t 2 L 2 + µ 8 u t 2 L 2 CM2 (c )c 2 c µ 8 u t 2 L 2, K 2 C ρ t 2 L 6 v 2 L 6 u 2 L 2 + µ 8 u t 2 L 2 Cc2 2 c2 1 c2 u 2 L 2 + µ 8 u t 2 L 2, K 3 C ρ L v 2 L ρu t 2 L 2 + µ 8 u t 2 L 2 Cc2 2c ρu t 2 L 2 + µ 8 u t 2 L 2, K 4 η v t 2 L 2 ρu t 2 L 2 +η 1 ρ L u L 2 u H 1 η v t 2 L 2 ρu t 2 L 2 +η 2 c 2 u 2 L 2 + u 2 H 1, K 5 C( 2 d L 3 d t L 2 u t L 6 + d t L 2 d L 3 u t L 6 d 2 L + d t L 2 d L 3 u t L 6) C 2 d L 2 2 d L 6 d t 2 L 2 +C d t 2 L 2 d L 2 d H 1 d 4 H 2 +C d t 2 L 2 d L 2 d H 1 + µ 8λ u t 2 L 2 Cc 2 2c 6 1c 21 + µ 8λ u t 2 L 2,

14 14 K 6 C( 2 d L 3 d t L 2 u t L 6 + d L d t L 2 u t L 2 + d L 2 d t L 3 u t L 6 + d L 3 d t L 6 d 2 L 6 u t L 6) C( 2 d L 2 2 d L 6 d t 2 L 2 + d 2 W 1,6 d t 2 L 2 + d t H 1 d t L 2 d 2 L 2 + d L 2 d H 1 d t 2 L 2 d 4 H 1) On the other hand + µ 8λ u t 2 L 2 Cc 4 2 c2 1 c17 + µ 8λ u t 2 L 2. d dt u 2 L 2 = u : u t dx C u 2 L + µ 2 8 u t 2 L2. (2.4) Substituting K 1 K 6 into (2.39) and combing with (2.4), we obtain d (ρ u t 2 + u 2 )dx+ u t 2 dx dt where A η (t)( ρu t 2 L 2 + u 2 L 2)+B(t)+C u 2 H1, (2.41) A η (t) = C(c 2 2c 2 1c 2 +η v t 2 L 2 +η 2 c 2 ), B(t) = Cc 4 2c 4 1c 21 M 2 (c ). Taking η = c 1 2 and using the estimates (2.13)-(2.14) and (2.19), we have t [,T 3 ], A η(s)ds C +Cc 4 2t, B(s)ds Cc4 2c 4 1c 21 M 2 (c )t+c. (2.42) Multiplying (2.3) by u t, then integrating it over and using the Young s inequality, we can obtain ρ u t 2 dx(τ) C ρ v 2 u 2 + ρ 1 2 (µ u λdiv( d d 1 2 ( d 2 +F(d))I) p) 2 dx(τ). Hence lim sup τ + ρ u t 2 dx(τ) C(c 5 + g 2 L 2) Cc5. (2.43)

15 15 Integrating (2.41) with respect to time (τ,t) and letting τ +, thanks to (2.43) and Gronwall s inequality and the inequality (2.42), we have for all t [,T 3 ], (ρ u t 2 + u 2 )dx(t)+ u t 2 dxdτ Cc 7 M2 (c )+C u 2 H1dτ. (2.44) Using the elliptic regularity result to the linear movement equation (2.3), we can estimate the term 2 u L 2 as follows u H 1 C( ρu t L 2 + ρv u L 2 + p L 2 + u L 2 + ( d) T ( d f(d)) L 2). (2.45) From the assumption (2.8)-(2.11) and the estimates (2.13), (2.14) and (2.19), we can derive ( d) T d L 2 C ( d) T (d t +v d) L 2 +C ( d) T [(n+m) (d m)]n L 2 C( d 1 2 L 2 d 1 2 H 1 d t L 2 + d H 1 v H 1 d H 1 + d H 2 n+m L 6 d m H 2 n L 6) Cc 6 1c 6, (2.46) ( d) T f(d) L 2 C d H 1 d+m L 6 d m L 6 d H 2 Cc 4 1c 12 (2.47) and ρu t L 2 + ρv u L 2 + p L 2 + u L 2 ρ L ρu t L 2 + ρ L v L 6 u L 3 +CM(c )c + u L 2 Cc 1 2 ρu t L 2 +Cc c 1 u 1 2 L 2 u 1 2 H 1 +CM(c )c + u L 2 CM(c )c +Cc 4 1( u L 2 + ρu t L 2)+ 1 2 u H1. (2.48) Taking (2.46)-(2.48) into (2.45), we can deduce u H 1 Cc 6 1 c1 +CM(c )c +Cc 4 1 ( u L 2 + ρu t L 2). (2.49)

16 16 Taking (2.49)into (2.44) andusing Gronwall s inequality, we canderive forall t [,T 3 ] (ρ u t 2 + u 2 )dx+ u t 2 dxdτ Cc 1 M 2 (c ). (2.5) So substituting (2.5) into (2.49), we get u H 1 Cc 6 1 c1 M(c ). (2.51) Using the elliptic regularity result to the linear movement equation(2.3), the term 2 u L 6 can be controlled as follows 2 u L 6 C( ρu t L 6 + ρv u L 6 + p L 6 + ( d) T ( d f(d)) L 6 + u L 6). (2.52) It follows from the assumption (2.8)-(2.11) and the estimates (2.13), (2.14) and (2.51) that ρu t L 6 + ρv u L 6 + p L 6 + u L 6 ρ L u t L 6 + ρ v L u L 6 +CM(c )c +Cc 6 1 c1 M(c ) Cc u t L 2 +Cc 2 c 6 1c 11 M(c ). (2.53) It follows from the estimates (2.27) and (2.35) that ( d) T ( d f(d)) L 6 C d 2 H 3 +C d L 6 d m H 2 d+m H 2 d H 2 Cc 4 2c 2 1c 6 +Cc 8 1c 12 Cc 4 2 c4 1 c12. (2.54) Taking (2.53) and (2.54) into (2.52), integrating it over time and using the estimate (2.5), we can derive for all t [,T 3 ], 2 u 2 L 6dτ Cc2 u t 2 L2dτ +Cc18 M2 (c ) Cc 18 M2 (c ). It is obvious that Lemma 2.1-Lemma 2.3 imply Theorem 2.1.

17 Local existence of a solution to the linear problem (2.1)- (2.3) Since the estimates obtained in the Lemma 2.1-Lemma 2.3 are uniform for all small δ, we have the following theorem: Theorem 2.2. If the initial condition (ρ,u,d ) satisfies the regularity (1.8) and the compatibility condition (1.9), then there exists a unique strong solution (ρ, u, d) to the linear equations (2.1)-(2.3) with initial boundary value (1.5)-(1.6) such that ρ C([,T 3 ];W 1,6 ), ρ t C([,T 3 ];L 6 ), u C([,T 3 ];H 1 H2 ) L 2 (,T 3 ;W 2,6 ), u t L 2 (,T 3 ;H 1 ), d C([,T 3 ];H 2 ) L 2 (,T 3 ;H 3 ), d t C([,T 3 ];H 1 ) L2 (,T 3 ;H 2 ), ρut C([,T 3 ];L 2 ). (2.55) Moreover, (ρ, u, d) also satisfies the inequalities (2.13)-(2.14), (2.19) and (2.37). Before proof, we give the following two classical lemmas which are proved in the book [16]. Lemma 2.4. Let Y = {v v L α (,T;X ),v t L α 1 (,T;X 1 )} with norm v Y = v L α (,T;X )+ v t L α 1(,T;X1 ) where 1 < α,α 1 <,X X X 1 are Banachspaces and X,X 1 are reflexive. Suppose that the injections X X X 1 are continuous, and the injection from X into X is compact. Then the injection from Y into L α (,T;X) is compact. Lemma 2.5. (An Interpolation Theorem) Let V, H, V be three Hilbert spaces, each space included in the following one as V H H V, V being the dual of V. If a function u belongs to L 2 (,T;V) and its derivative u t belongs to L 2 (,T;V ), then u is almost everywhere equal to a function continuous form [,T] into H. Proof. We now start to prove Theorem 2.2. From the Theorem 2.1 and the estimates (2.13), (2.14), (2.19) and (2.37), by the compactness Lemma 2.4, there exists (ρ, d, u)

18 18 such that (ρ δ,d δ,u δ ) (ρ,d,u) in L 2 (,T 3 ;L r H 2 W 1,q ), (2.56) p δ p in L (,T 3 ;W 1,6 ), as δ, (2.57) where r (1,+ ) and q [2,+ ). Hence p = p, a.e.. Because of the lower semi-continuity of various norms, the estimates (2.13), (2.19) and (2.37) also hold for (ρ,u,d). So for almost every (t,x)( [,T 3 ] ), (ρ,u,d) satisfies the system (2.1)-(2.3) which means (ρ,u,d) is a strong solution to the linear equation (2.1)-(2.3) with initial-boundary conditions (1.5)-(1.6). The solution (ρ,u,d) is unique: From the Lemma 2.1, ρ is the unique solution of the linear equation (2.1). Using the same method as section 4, we can prove d and u are the unique solution to the linear equations (2.2) and (2.3) respectively. Finally, we will prove the time continuity of the solution (ρ,u,d). The solution from the Lemma 2.1 is the same as from the approximation (2.56) due to the uniqueness of solution. So we get From the linear equation (2.1), we easily show By the interpolation Lemma 2.5, (2.37) can deduce that ρ C([,T];W 1,6 ). (2.58) ρ t C([,T 3 ];L 6 ). (2.59) d t L 2 (,T 3 ;H 2 ), d L 2 (,T 3 ;H 3 ) d C([,T 3 ];H 2 ). (2.6) Differentiating the linear equation (2.2) with respect to time and space, we get d tt + (v d) t = ν( d t 1 σ 2 [(n+m) (d m)n] t ). By the estimate (2.19), we deduce d tt L 2 (,T 3 ;H 1 ). Because d t L 2 (,T 3 ;H 1 ), by the interpolation Lemma 2.5 again, we have d t C([,T 3 ];H 1 ). (2.61)

19 19 From the linear equation (2.2), we get d C([,T 3 ];H 1 ). (2.62) By the interpolation Lemma 2.5, (2.37) can deduce that u t L 2 (,T 3 ;H 1 ), u L 2 (,T 3 ;W 2,6 ) u C([,T 3 ];H 1 ). (2.63) From the linear equation (2.3) and the estimates (2.13), (2.19) and (2.37), we obtain (ρu t,(ρu t ) t ) L 2 (,T 3 ;H 1 ) L 2 (,T 3 ;H 1 ) ρu t C([,T 3 ];L 2 ). (2.64) From the linear equation (2.3) and the elliptic regularity estimate 2 u L 2 C u L 2, we have u C([,T 3 ];H 2 ). (2.65) So we get time-continuity of (ρ, u, d) from (2.58)-(2.65). 3 Iteration and existence in Theorem 1 Set c 1 = Cc 18 M 2 (c ), c 2 = c 7 1, T 3 = min{c 22 2,T}. (3.1) At the beginning, let s choose a initial data of iteration (u (t,x),d (t,x)). u (t,x) satisfies the following heat equation φ t φ =, with φ t= = u, φ =, and d (t,x) = d (x). Because (c 1,c 2,T 3 ) depends only on c, we can choose a T ( [,T 3 ]) so small that (2.9)-(2.11) hold for u and d with T instead of T. Replacing (v,n) by (u,d ) and using the Theorem 2.2, we can obtain the solution (ρ 1,u 1,d 1 ) of (2.1)-(2.3) with (1.5)-(1.6), and it satisfies the estimates (2.13), (2.19)

20 2 and (2.37). Inductively, for all k N +, replacing (v,n) by (u k 1,d k 1 ) and using the Theorem 2.2, we can obtain a sequence (ρ k,u k,d k ) of solution of (2.1)-(2.3) with (1.5)- (1.6), and they satisfy the following estimates with the same (c,c 1,c 2,T ) independent of k N + : and sup ( u k H 1 + d k H 1 + d k t H 1 +c 6 1 ( 2 u k L d k L 2)) t T + T u k t 2 L 2 + uk 2 W 2,6 + 2 d k t 2 L 2 + dk 2 H 3dt Cc 22 M2 (c ) (3.2) sup ( ρ k W 1,6 + ρ k t L 6) Cc 2 c, sup ρ k u k t L 2 CM(c )c 5, (3.3) t T t T sup ( p k W 1,6 + p k t L 6) CM(c )c 2 c, sup d k H 2 Cc 2 2 c 1c 3. (3.4) t T t T We will show (ρ k,u k,d k ) converges to a strong solution to the original nonlinear problem (1.1)-(1.3). Define have ρ k+1 = ρ k+1 ρ k, d k+1 = d k+1 d k, u k+1 = u k+1 u k. (3.5) Since (ρ k,u k,d k ) and (ρ k+1,u k+1,d k+1 ) satisfy the linear equations (2.1)-(2.3), we ρ k+1 t +div(ρ k+1 u k )+div(ρ k u k ) =, (3.6) d k+1 t ν d k+1 = u k d k+1 u k 1 d k+1 ν σ 2((dk +m)(d k+1 m))d k ν σ 2((dk +m) d k+1 )d k 1 ν σ 2(dk (d k m))d k 1, (3.7) ρ k+1 u k+1 t +ρ k+1 u k u k+1 µ u k+1 + (p k+1 p k ) = λ( d k+1 ) T ( d k+1 f(d k+1 )) λ( d k ) T ( d k+1 (f(d k+1 ) f(d k ))) ρ k+1 (u k 1 u k +u k t) ρ k+1 u k u k. (3.8)

21 21 Define where Ψ k+1 = ρ k+1 2 L 2 + dk+1 2 L 2 + dk+1 2 L + ρ k+1 u k L2. (3.9) Before estimates, we introduce two small positive unfixed constants η and ǫ. From the first equation (3.6), we can derive d dt ρk+1 2 L 2 Ak η(t) ρ k+1 2 L 2 +η uk 2 L2, (3.1) A k η(t) = C u k (t) W 1,6 +η 1 C( ρ k (t) 2 L 3 + ρk (t) 2 L ). Using the uniform estimates (3.2)-(3.4), we obtain A k η (s)ds C + C η t, t [,T ]. (3.11) Multiplying (3.7) by d k+1 and integrating it over, we have d d k+1 2 dx+ d k+1 2 dx dt where for all t [,T ], B k η (t) dk+1 2 L 2 +η( uk 2 L 2 + dk 2 L2), (3.12) B k η (t) = Cη 1 d k+1 2 L 2 +C uk 1 2 L +Cη 1 d k +m 2 L 6 dk m 2 L 6 +C d k +m L d k 1 L +Cη 1 d k 1 2 L 6 dk m 2 L 6. The uniform estimates (3.2)-(3.4) implies B k η (s)ds C(1+ 1 )t. (3.13) η Multiplying (3.7) by d k+1, integrating it over and using the elliptic estimate 2 d k+1 L 2 C d k+1 L 2, we can deduce by integration by parts d d k+1 2 dx+ 2 d k+1 2 dx C dt 13 i=1 L i, (3.14)

22 22 where L 1 = u k L 2 d k+1 L d k+1 L 2 Cη 1 d k+1 2 H 2 dk+1 2 L 2 +η uk 2 L 2, L 2 = u k L 6 2 d k+1 L 3 d k+1 L 2 Cη 1 2 d k+1 2 H 1 dk+1 2 L 2 +η uk 2 L 2, L 3 = u k 1 W 1,6 d k+1 2 L 2, L 4 = u k 1 L 2 d k+1 L 2 d k+1 L 2 Cǫ 1 u k 1 2 H 2 dk+1 2 L 2 +ǫ 2 d k+1 2 L 2, L 5 = d k L 3 d k+1 m L d k L 6 d k+1 L 2 Cη 1 d k 2 H 1 dk+1 m 2 H 2 dk+1 2 L 2 +η dk 2 L 2, L 6 = d k L (d k+1 m) L 3 d k L 6 d k+1 L 2 Cη 1 d k 2 H 2 (dk+1 m) 2 H 1 dk+1 2 L 2 +η dk 2 L 2, L 7 = d k L d k+1 m L d k L 2 d k+1 L 2 Cη 1 d k 2 H 2 dk+1 m 2 H 2 dk+1 2 L 2 +η dk 2 L 2, L 8 = (d k +m) L 3 d k+1 L 6 d k 1 L d k+1 L 2 C (d k +m) 2 H 1 dk 1 H 2 d k+1 2 L 2, L 9 = d k +m H 2 d k 1 H 2 d k+1 2 L 2, L 1 = d k +m L d k+1 L 6 d k 1 L 3 d k+1 L 2 C d k +m 2 H 2 dk 1 H 1 d k+1 2 L 2,

23 23 L 11 = d k L 2 d k m L d k 1 L d k+1 L 2 Cη 1 d k m 2 H 2 dk 1 2 H 2 dk+1 2 L 2 +η dk 2 L 2, L 12 = d k L 6 (d k m) L 3 d k 1 L d k+1 L 2 Cη 1 (d k m) 2 H 1 dk 1 2 H 2 dk+1 2 L 2 +η dk 2 L 2, L 13 = d k L 6 d k+1 m L d k 1 L 3 d k+1 L 2 Cη 1 d k m 2 H 2 dk 1 2 H 1 dk+1 2 L 2 +η dk 2 L 2. Let s ǫ small enough, the inequality (3.14) becomes d d k+1 2 dx+ 2 d k+1 2 dx dt CC k η(t) d k+1 2 L 2 +Cη( dk 2 L 2 + uk 2 L2), (3.15) where C k η (t) = η 1 d k+1 2 H 2 + uk 1 W 1,6 + u k 1 2 H 2 +η 1 d k 2 H 1 dk+1 m 2 H 2 +η 1 d k 2 H 2 (dk+1 m) 2 H 1 +η 1 d k 2 H 2 dk+1 m 2 H 2 + (d k +m) 2 H 1 dk 1 H 2 + d k +m H 2 d k 1 H 2 + d k +m 2 H 2 dk 1 H 1 +η 1 (d k m) 2 H 1 dk 1 2 H 2 +η 1 (d k m) 2 H 1 dk 1 2 H 2 +η 1 d k m 2 H 2 dk 1 2 H 1. The uniform estimates (3.2)-(3.4) implies C k η(s)ds C +(C + C η )t, t [,T ]. (3.16) Multiplying the movement equation (3.8) by u k+1 and integrating over, we can

24 24 deduce 1 d ρ k+1 u k+1 2 dx+µ u k+1 2 dx 2dt ρ k+1 ( u k 1 u k + u k t ) uk+1 dx+ ρ k+1 u k u k u k+1 dx +λ d k+1 d k+1 f(d k+1 ) u k+1 dx+λ 2 d k d k+1 u k+1 dx +λ d k d k+1 u k+1 dx+ 1 d k d k+1 +d k d k+1 d k+1 u k+1 dx σ d k d k +m d k m d k+1 u k+1 dx+ p k+1 p k u k+1 dx σ 2 = 8 M i, i=1 where we have used (3.17) f(d k+1 ) f(d k ) = 1 σ 2(dk+1 (d k+1 +d k )d k+1 +( d k 2 1)d k+1 ). Here M 1 ρ k+1 L 2 u k+1 L 6( u k 1 L 6 u k L 6 + u k t L 3) C ρ k+1 2 L 2( uk 1 2 H 1 uk 2 W 2,6 + uk t L 2 u k t L 2)+ µ 9 uk+1 2 L 2, M 2 ρ k+1 L 6 u k L 6 u k L 6 ρ k+1 u k+1 L 2 η 1 ρ k+1 2 L 6 uk 2 H 2 ρ k+1 u k+1 2 L 2 +η uk 2 L 2, M 3 C d k+1 L 2 d k+1 f(d k+1 ) L 3 u k+1 L 6 µ 9 uk+1 2 L 2 +C dk+1 2 L 2( 2 d k+1 L 2 2 d k+1 L 6 + d k+1 +m 2 H 1 dk+1 m 2 H 1 dk+1 2 H 2),

25 25 M 4 C 2 d k L 3 d k+1 L 2 u k+1 L 6, C 2 d k L 6 2 d k L 2 d k+1 2 L 2 + µ 9 uk+1 2 L 2 M 5 C d k L d k+1 L 2 u k+1 L 2 C d k 2 W 1,6 dk+1 2 L 2 + µ 9 uk+1 2 L 2, M 6 C d k L 3 d k+1 +d k L 6 d k+1 L 6 d k+1 L 6 u k+1 L 6 C d k L 2 d k H 1 d k+1 +d k 2 H 1 dk+1 2 H 1 dk+1 2 L 2 + µ 9 uk+1 2 L 2, M 7 C d k L 3 d k +m L 6 d k m L 6 d k+1 L 6 u k+1 L 6 C d k L 2 d k H 1 d k +m 2 H 1 dk m 2 H 1 dk+1 2 L 2 + µ 9 uk+1 2 L 2, M 8 p k+1 p k L 2 u k+1 L 2 CM 2 (c ) ρ k+1 2 L 2 + µ 9 uk+1 2 L 2. Taking M 1 M 8 into (3.17), we obtain d ρ k+1 u k+1 2 dx+ u k+1 2 dx dt where CD k η (t)( ρk+1 2 L + ρ k+1 u k L 2 + dk+1 2 L 2)+Cη uk 2 L2, (3.18) D k η(t) = u k 1 2 H 1 uk 2 W 2,6 + uk t L 2 u k t L 2 +η 1 ρ k+1 2 L 6 uk 2 H d k+1 L 2 2 d k+1 L 6 + d k+1 +m 2 H 1 dk+1 m 2 H 1 dk+1 2 H d k L 6 2 d k L 2 + d k 2 W 1,6 + dk L 2 d k H 1 ( d k+1 +d k 2 H 1 dk+1 2 H 1 + dk +m 2 H 1 dk m 2 H 1)+M2 (c ). From the uniform estimates (3.2)-(3.4), we have D k η (s)ds C +(C + C η )t, t [,T ]. (3.19)

26 26 where Summing (3.1), (3.12), (3.15) and (3.18), we obtain d dt Ψk+1 +( d k+1 2 L d k+1 2 L 2 + uk+1 2 L 2) CE k η(t)ψ k+1 +Cη( u k 2 L 2 + dk 2 L2), (3.2) E k η (t) = Ak η (t)+bk η (t)+ck η (t)+dk η (t). Using the uniform estimates (3.2)-(3.4), we obtain again E k η(s)ds C +C(1+η 1 )t, t [,T ]. (3.21) Applying Gronwall s inequality to (3.2), we can deduce Cη Ψ k+1 (t)+ ( d k+1 2 L d k+1 2 L 2 + uk+1 2 L 2)ds ( u k 2 L 2 + dk 2 L 2)dτ exp(c +C(1+η 1 )t). (3.22) Hence choose small constants η, T (< T ), so that Cηexp(C +C(1+η 1 )t) 1 2, t [,T ]. (3.23) We easily deduce that sup Ψ k+1 (t)+ t T k=1 C <. k=1 T ( d k+1 2 L d k+1 2 L 2 + uk+1 2 L 2)ds (3.24) (3.24) implies that the full sequence (ρ k,d k,u k ) converges to a limit (ρ,d,u) in the following strong sense ρ k ρ in L (,T ;L 2 ), u k u in L 2 (,T ;H 1 ), d k d in L (,T ;H 1 ) L 2 (,T ;H 2 ). (3.25) Hence the problem (1.1)-(1.3) with initial boundary data (1.5)-(1.6) has a weak solution (ρ, u, d). Furthermore using the estimates (3.2)-(3.4), we obtain that a subsequence of

27 27 (ρ k,u k,d k ) converges to (ρ,u,d) in an obvious weak or weak* sense. Due to the lower semi-continuity of various norms, from (3.2)-(3.4), (ρ, u, d) also satisfies the following regularity estimates sup ( ρ W 1,6 + ρ t L 6 + u H 1 t T + p W 1,6 + p t L 6 + d H 1 + d t H u L d L 2 + d H 2) T + ρu t 2 + u t 2 L 2 + u 2 W 2,6 + 2 d t 2 L 2 + d 2 H 3dt C. (3.26) Hence (ρ,u,d) is also a strong solution to the problem (1.1)-(1.3). 4 Uniqueness and continuity in Theorem 1 In this section, we will use energy method to prove the uniqueness and continuity in Theorem 1. For simplicity, we introduce some notations ρ = ρ ρ, u = u ũ, d = d d. Define Ψ(t) = ρ 2 L 2 + d 2 L 2 + d 2 L 2 + ρu 2 L 2. Using the similar process in the section 3, we can obtain the following estimate (see (3.22)) where F(t) L 1 (,T ). d dt Ψ+( d 2 L d 2 L 2 + u 2 L2) CF(t)Ψ, (4.1) Applying the Gronwall s inequality to (4.1), we get for all t [,T ], Ψ(t)+ ( 2 d 2 L 2 + u 2 L2)dτ CΨ(), (4.2)

28 28 which implies the uniqueness, and for all t [,T ], we have ( ρ 2 L 2 + d 2 H 1 + ρu 2 L 2)(t)+ as ( ρ,ũ, d ) (ρ,u,d ) in W 1,6 H 2 H 3. Because d and d satisfy (1.3), we obtain, similar to (3.7), d d 2 dx+ d t 2 dx dt ( 2 d 2 L 2 + u 2 L2)dτ (4.3) C( u 2 L 2 + d 2 L 2 ũ 2 L 6 + d 2 L 2)+ ũ 2 L d 2 L 2. Applying Gronwall s inequality to the above inequality, and using the inequality (4.3) and the elliptic estimate d H 2 C d L 2, we have, for all t [,T ], d H 2(t)+ as ( ρ,ũ, d ) (ρ,u,d ) in W 1,6 H 2 H 3. d t 2 L2dτ (4.4) Similarly from (3.8), using Gronwall s inequality, (3.26) and the convergence (4.3)- (4.4), we obtain, for all t [,T ], u H 1(t)+ T as ( ρ,ũ, d ) (ρ,u,d ) in W 1,6 H 2 H 3. ρu t 2 L2dt (4.5) Multiplying the difference between the continuity equations by 6ρ 5, integrating over (,t) and then using Gronwall s inequality, it follows from the estimate (3.26) and the convergence (4.5) that for all t [,T ], ρ ρ L 6(t). (4.6) From the equations (1.2) and (1.3), by a simple discussion, we can obtain for all t [,T ], d L 2(t), d L 2 (,T ;H 3 ), u L 2 (,T ;H 2 ) (4.7) as ( ρ,ũ, d ) (ρ,u,d ) in W 1,6 H 2 H 3. In conclusion, (4.3)-(4.7) complete the proof of the continuity in Theorem 1.

29 29 5 Proof of Theorem 2 Suppose that there are two positive constants θ(< 1) and C such that max{ ρ W 1,6, u H 2, d m H 3, g 2 L2} < θ, (5.1) sup ( v H 2 + n m H 2 + n t H 1 )+ t T T v t 2 L 2 + v 2 W 2,6 + 2 n t 2 L 2 + n 2 H 3dt < C. (5.2) In this section, we assume the genuine constant C, maybe depending on the constant M(1) which is defined by (2.18). By Lemma 2.1, there exists a small θ 1 (< 1) so that t [,T], θ (,θ 1 ], ρ W 1,6 Cθ 1 2, ρ t L 6 Cθ 1 3, p W 1,6 Cθ 1 2, p t L 6 Cθ 1 3, (5.3) From Lemma 2.2, we can find a small θ 2 (< 1) so that θ (,θ 2 ], d t 2 H 1(t), d m 2 H 2(t), d m 2 H 3dτ Cθ 1 2, t [,T]. (5.4) By Lemma 2.3, a small θ 3 ( min{θ 1,θ 2 }) also can be found so that θ (,θ 3 ], u 2 H 2(t), ρu t 2 L 2, u t 2 H 1dτ, u 2 W 2,6dτ Cθ 1 6, t [,T]. (5.5) Thanks to the estimates (5.3)-(5.5), using the Theorem 2.2, we can obtain the global strong solution of the linear system (2.1)-(2.3) with initial boundary value (1.5) and (1.6) provided max{ ρ W 1,6, u H 2, d m H 3, g 2 L2} θ, where θ (,θ 3 ]. Now let s talk about the iteration. First, we notice that if θ 3 is taken so small that Cθ 3 C, then the process of iteration can be continued for the same θ 3. Next we will pay attention to the convergence of the iteration.

30 3 As the same process in the section 3, using the estimates (5.3)-(5.5), we can obtain, like (3.22), Ψ k+1 (t)+ ( d k+1 2 L d k+1 2 L 2 + uk+1 2 L 2)ds Cηexp(C(t+θ t+θ 6 +η 1 θ 1 2 t+η 1 θ 1 2 )) ( u k 2 L 2 + dk 2 L2)dτ, (5.6) where Ψ k+1 = ρ k+1 2 L 2 + dk+1 2 L 2 + dk+1 2 L 2 + ρ k+1 u k+1 2 L 2. Hence choose small constants η, θ, so that θ (,θ ] and t [,T], We easily deduce that t [,T], k=1 sup t T C. Cηexp(C(t+θ 1 6 t+θ 1 6 +η 1 θ 1 2 t+η 1 θ 1 2 )) 1 2, (5.7) Ψ k+1 (t)+ k=1 T ( d k+1 2 L d k+1 2 L 2 + uk+1 2 L 2)ds (5.8) So we complete the proof of Theorem 2. References [1] C. Liu, X.G. Liu, J. Qing, Existence of globally weak solutions to the flow of compressible liquid crystals system, preprint. [2] P.L. Lions, Mathematical topics in fluid dynamics, Vol.2. Compressible models. Oxford Science Publication, Oxford, [3] E. Feireisl, Dynamics of viscous compressible fluids, Oxford University Press, Oxford, 24.

31 31 [4] F.H. Lin, C. Liu, Nonparabolic dissipative systems modeling the flow of liquid crystals, Comm. Pure Appl. Math. 48(1995) [5] M.C. Calderer, C. Liu, Liquid crystal flow: dynamic and static configurations, SIAM J. Appl. Math. 6(2) [6] F.H. Lin, C. Liu, Existence of solutions for the Ericksen-Leslie system, Arch. Rational Mech. Annal. 154(2) [7] F.H. Lin, C. Liu, Static and dynamic theories of liquid crystals, J. Partial Differential Equations 14(21) [8] D. Coutand, S. Shkoller, Well posedness of the full Ericksen-Leslie model of nematic liquid crystals, Note C.R.A.S, t. Serie I 333 (21) [9] H.Y. Wen, S.J. Ding, Solutions of incompressible hydrodynamic flow of liquid crystals, Nonlinear Analysis: Real World Applications, article in press(21). [1] F.H. Lin, J. Lin, C.Y. Wang, Liquid crystals flows in two dimensions, Arch. Rational Mech. Anal. 197(21) [11] R. Salvi, I. Străskraba, Global existence for viscous compressible fluids and their behavior as t, J. Fac. Sci. Univ. Tokyo Sect. IA Math. 4(1993) [12] H.J. Choe, H. Kim, Strong solutions of the Navier-Stokes equations for isentropic compressible fluids, J. Differential Equations 19(23) [13] X.G. Liu, L.M. Liu, A blow-up criterion for the compressible liquid crystals system, arxiv: v2(21). [14] Y. Cho, H.J. Choe, H. Kim, Unique solvability of the initial boundary value problems for compressible viscous fluids, J. Math. Pures Appl. 83(24)

32 32 [15] Y. Cho, H. Kim, Existence results for viscous polytropic fluids with vacuum, J. Differential Equations 228(26) [16] R. Temam, Navier-Stokes equations, rev.ed., Studies in Mathematics and its Applications 2, Amsterdam, 1977.