Skip to main content

Unique solvability for the non-Newtonian magneto-micropolar fluid

Abstract

In this paper, a motion of an incompressible non-Newtonian magneto-micropolar fluid is considered. We assume that the stress tensor has a p-structure, and we establish the global in time existence and uniqueness of the weak solutions with p 5 2 in three dimensions.

1 Introduction and main results

This paper is concerned about the existence and uniqueness of the weak solutions to the non-Newtonian magneto-micropolar fluid equations in Ω×(0,T), which are described by

{ u t + ( u ) u div ( | e ( u ) | p 2 e ( u ) ) + ( π + 1 2 | b | 2 ) = χ rot ω + ( b ) b , ω t + ( u ) ω μ Δ ω + 2 χ ω = χ rot u , b t λ Δ b + ( u ) b ( b ) u = 0 , div u = div b = 0 ,
(1.1)

here Ω R 3 is an open-bounded domain with Lipschitz boundaries, and the unknowns u, ω, b, π denote the velocity of the fluid, the micro-rotational velocity, magnetic field and hydrostatic pressure, respectively. χ, μ, λ are positive numbers associated with properties of the material: χ is the vortex viscosity, μ is spin viscosity and 1 λ is the magnetic Reynold. In (1.1), e(u) is the symmetric part of the velocity gradient, i.e.,

2e(u)=u+ ( u ) T .

To (1.1) we append the following initial and boundary conditions

u(x,0)= u 0 (x),ω(x,0)= ω 0 (x),b(x,0)= b 0 (x),xΩ,
(1.2)
u(x,t)=ω(x,t)=b(x,t)=0,(x,t) Σ T ,
(1.3)

where Σ T =Ω×(0,T).

The theory of micropolar fluid was first proposed by Eringen [1] in 1966, which enabled us to consider some physical phenomena that cannot be treated by the classical Navier-Stokes equations for the viscous incompressible fluid, for example, the motion of animal blood, liquid crystals and dilute aqueous polymer solutions, etc.

If p=2, system (1.1) reduces to the classical magneto-micropolar fluid equations, and there are many earlier results concerning the weak and strong solvability in a bounded domain Ω R 3 . For strong solutions, Galdi and Rionero [2] stated, without proof, the results of existence and uniqueness of strong solutions. Rojas-Medar [3] studied it and established the local in time existence and uniqueness of strong solutions by the spectral Galerkin method. In 1999, Ortega-Torres and Rojas-Medar [4] proved global existence of strong solutions with the small initial values. For weak solutions, Rojas-Medar and Boldrini [5] proved the existence of weak solutions, and in the 2D case, also proved the uniqueness of the weak solutions.

On the other hand, there are few existence results about the non-Newtonian magneto-micropolar fluid, i.e., the p2 case. In a recent work, Gunzburger et al. [6] studied the reduced problem (with both ω=0 and χ=0), and gave the global unique solvability of the first initial-boundary value problem in a bounded two or three-dimensional domain. Improved results are proved for the periodic boundary condition case.

In this paper, we will prove the global existence and uniqueness of the weak solutions for the full system (1.1)-(1.3) under the condition that p 5 2 . These results are based on the Galerkin method and a series of uniform estimates, which do not depend on the parameters.

Throughout this work, we use a standard notation L p (Ω) (normed p ) for Lebesgue L p -spaces, as well as W k , p (Ω) (normed k , p ) for the usual Sobolev spaces. As usual, C 0 (Ω) denotes the set of all C -functions with the compact support in Ω. Given T>0 and a Banach space X, we denote by L q (0,T;X) Bochner spaces, which are equipped with the norm

L q ( 0 , T ; X ) := ( 0 T X q d s ) 1 q .

We also introduce the following functional vector spaces:

V { u C 0 ( Ω ) , div u = 0 } , H the closure of  V  in  L 2 ( Ω ) , V p the closure of  V  in  W 1 , p ( Ω ) .

We next introduce the definition of a weak solution for problems (1.1)-(1.3).

Definition 1.1 We say that (u,ω,b) is a weak solution to problems (1.1)-(1.3) if

u L ( 0 , T ; V p ) , u t L 2 ( 0 , T ; H ) , ω L ( 0 , T ; W 1 , 2 ( Ω ) ) L 2 ( 0 , T ; W 2 , 2 ( Ω ) ) , b L ( 0 , T ; V 2 ) L 2 ( 0 , T ; W 2 , 2 ( Ω ) ) , ω t L 2 ( 0 , T ; W 1 , 2 ( Ω ) ) , b t L 2 ( 0 , T ; H )
(1.4)

satisfy

u t , ϕ + | e ( u ) | p 2 e ( u ) , e ( ϕ ) = χ rot ω , ϕ + ( b ) b , ϕ ( u ) u , ϕ for all  ϕ V p ,
(1.5)
ω t , ψ + μ ω , ψ = χ rot u , ψ ( u ) ω , ψ 2 χ ω , ψ for all  ψ W 0 1 , 2 ( Ω ) ,
(1.6)
b t , η + λ b , η = ( b ) u , η ( u ) b , η for all  η V 2 ,
(1.7)

where the symbol , denotes a generic duality pairing.

The following theorem gives the main results of this paper.

Theorem 1.1 Let Ω R 3 be an open-bounded domain with a Lipschitz boundary Ω. Assume that p 5 2 , u 0 W 1 , p (Ω), ω 0 , b 0 W 1 , 2 (Ω). Then, for T(0,+), there exists a unique weak solution to problem (1.1)-(1.3) in the sense of Definition  1.1.

Remark 1.1 If (1.4)-(1.5) hold, it could be easy to introduce the pressure π L (0,T; L p (Ω)), p =p/(p1). This will be done at the end of Section 3.

For latter use, let us state some useful inequalities.

Lemma 1.1 (See [7]) (Korn’s inequality)

Let 1<p<. Then there exists a constant C p = C p (Ω) such that

C p v 1 , p e ( v ) p for all v W 0 1 , p (Ω),
(1.8)

where Ω R n is open and bounded with a Lipschitz boundary.

Lemma 1.2 (See [8]) (On negative norm)

Let 1<p<, and let v W 0 1 , p (Ω). Then there exists a constant C such that

C v p v 1 , p + v 1 , p .

Lemma 1.3 (See [9])

Let v,w W 1 , p (Ω). For each 2<p<, there exists a constant C = C (p)>0 such that

e ( v ) e ( w ) p p C Ω ( | e ( v ) | p 2 e ( v ) | e ( w ) | p 2 e ( w ) ) ( e ( v ) e ( w ) ) .

By using Hölder’s inequality and the imbedding inequality, we could arrive at

u m C(m,r) u r + C ˜ (m,r) u 2
(1.9)

with

{ m 3 r / ( 3 r ) for  r [ 1 , 3 ) , for any  m < for  r = 3 , m = for  r > 3 .

Here, C ˜ (m,r)=0 if u | Ω =0 or Ω udx=0. We will also apply the so-called multiplicative inequalities

u q C(q) u 2 α u 2 1 α + C ˆ (q) u 2
(1.10)

with

α=3 ( 1 2 1 q ) [0,1],q[2,6].

If u | Ω =0 or Ω udx=0, then C ˆ (q)=0.

Finally, the paper is organized as follows. In Section 2, we focus on the derivation of the priori estimates for the smooth solutions. On the bases of these estimates, in Section 3, we get the existence result with the help of the Galerkin method. The aim of Section 4 is to give the uniqueness criterion.

2 The priori estimates

Let (u,ω,b) be a smooth solution to system (1.1)-(1.3). The goal of this section is to derive some priori estimates about it. In all the following sections, we always assume that p 5 2 holds.

Setting ϕ=u in (1.5), ψ=ω in (1.6), η=b in (1.7), and observing that (u)u,u=(u)ω,ω=(u)b,b=0, we obtain

1 2 d d t u 2 2 + e ( u ) p p = χ rot ω , u + ( b ) b , u , 1 2 d d t ω 2 2 + μ ω 2 2 + 2 χ ω 2 2 = χ rot u , ω , 1 2 d d t b 2 2 + λ b 2 2 = ( b ) u , b .

Adding the identities above, noting that (b)b,u+(b)u,b=0, and Korn’s inequality (1.8), we get

1 2 d d t ( u 2 2 + ω 2 2 + b 2 2 ) + ( C p u p p + μ ω 2 2 + 2 χ ω 2 2 + λ b 2 2 ) = χ rot ω , u + χ rot u , ω C χ ( ω 2 u 2 + u 2 ω 2 ) ε ( ω 2 2 + u p p ) + C ε ( ω 2 2 + u 2 2 ) ( for  p 2 ) .

After choosing ε properly small, integrating over (0,t), t(0,T], the Gronwall’s inequality yields that

sup t ( 0 , T ) ( u 2 2 + ω 2 2 + b 2 2 ) + 0 T ( u p p + ω 2 2 + b 2 2 ) C 1 ,
(2.1)

where C 1 is a constant depending on the time T and u 0 2 , w 0 2 , b 0 2 .

Next, we derive the higher order estimates for ω and b. Setting ψ=Δω in (1.6), we find

1 2 d d t ω 2 2 + μ Δ ω 2 2 = ( u ) ω , Δ ω + 2 χ ω , Δ ω χ rot u , Δ ω | ( u ) ω , Δ ω | + ε Δ ω 2 2 + C ε ( ω 2 2 + u 2 2 )
(2.2)

for the first term on the right hand side, we compute by the divergence free conditions

| ( u ) ω , Δ ω | = | i j k k u i i ω j , k ω j | C Ω | u | | ω | 2 C u p ω 2 p p 1 2 C u p ( ω 2 2 p 3 p Δ ω 2 3 p + ω 2 2 ) ε Δ ω 2 2 + C ε u p 2 p 2 p 3 ω 2 2 + C u p ω 2 2 ,
(2.3)

where Hölder’s, Young’s inequality and (1.9), (1.10) have been used.

Inserting (2.3) into (2.2), choosing ε= μ 4 and integrating over (0,t), t(0,T], we have

1 2 ω 2 2 + μ 2 0 t Δ ω 2 2 C 0 t ( u p 2 p 2 p 3 + u p ) ω 2 2 + C 0 t ( ω 2 2 + u 2 2 ) + 1 2 ω 0 2 2 (2.1) C 0 t ( u p 2 p 2 p 3 + u p ) ω 2 2 + 1 2 ω 0 2 2 + C 1 ,
(2.4)

since p 5 2 , we have

2 p 2 p 3 p.
(2.5)

Gronwall’s inequality and estimate (2.1) now provide the bound

sup t ( 0 , T ) ω 2 2 + 0 T Δ ω 2 2 C 2 ,
(2.6)

where C 2 is a constant depending on the time T, C 1 and ω 0 2 .

Next, set η=Δb in (1.7) to discover

1 2 d d t b 2 2 + λ Δ b 2 2 = ( u ) b , Δ b ( b ) u , Δ b | ( u ) b , Δ b | + | ( b ) u , Δ b | .
(2.7)

Reasoning similar to (2.3), we could find

| ( u ) b , Δ b | ε Δ b 2 2 + C ε u p 2 p 2 p 3 b 2 2 +C u p b 2 2 .
(2.8)

For the second term on the right hand side of (2.7), we compute

| ( b ) u , Δ b | ε Δ b 2 2 + C ε ( b ) u 2 2 ε Δ b 2 2 + C ε u p 2 b 2 p p 2 2 ε Δ b 2 2 + C ε u p 2 b 6 p 5 p 6 2 ε Δ b 2 2 + C ε u p 2 ( b 2 2 ( 2 p 3 ) p Δ b 2 2 ( 3 p ) p + b 2 2 ) ε Δ b 2 2 + δ Δ b 2 2 + C ε δ u p 2 p 2 p 3 b 2 2 + C ε u p 2 b 2 2 ,
(2.9)

where we have used Hölder’s, Young’s inequality and (1.9), (1.10). Choosing ε and δ properly small, inserting (2.8)-(2.9) into (2.7) and integrating it over (0,t), t(0,T], we find

1 2 b 2 2 + λ 2 0 t Δ b 2 2 C 0 t ( u p 2 p 2 p 3 + u p 2 ) b 2 2 + 1 2 b 0 2 2 .
(2.10)

Observing (2.5) and estimate (2.1), then Gronwall’s inequality yields

sup t ( 0 , T ) b 2 2 + 0 T Δ b 2 2 C 3 ,
(2.11)

where C 3 is a constant depending on the time T, C 1 and b 0 2 .

Reasoning analogously to (2.6) and (2.11), it is easy to see that identity (1.6) with ψ= ω t , (1.7) with η= b t , with the help of (2.6) and (2.11), guarantee the estimate

0 T ω t 2 2 C 4 , 0 T b t 2 2 C 5 ,
(2.12)

where C 4 , C 5 are both constants depending only on the time T and some norm of the initial values.

In fact (we here only take b t as an example), set η= b t in (1.7), we deduce that

b t 2 2 + λ 2 d d t b 2 2 | ( b ) u , b t | + | ( u ) b , b t | ε b t 2 2 + C ε ( ( b ) u 2 2 + ( u ) b 2 2 ) .
(2.13)

Now, we compute, by using Hölder’s, Young’s inequality and (1.9), (1.10)

( b ) u 2 2 u p 2 b 2 p p 2 2 C ( u p p + b 2 p p 2 2 p p 2 ) C ( u p p + b 6 p 5 p 6 2 p p 2 ) C u p p + C ( b 2 2 ( 2 p 3 ) p 2 Δ b 2 2 ( 3 p ) p 2 + b 2 2 p p 2 )
(2.14)

and

( u ) b 2 2 u 3 p 3 p 2 b 6 p 5 p 6 2 C u p 2 b 6 p 5 p 6 2 C ( u p p + b 6 p 5 p 6 2 p p 2 ) C u p p + C ( b 2 2 ( 2 p 3 ) p 2 Δ b 2 2 ( 3 p ) p 2 + b 2 2 p p 2 ) .
(2.15)

Combining (2.13)-(2.15), by choosing ε= 1 2 , we arrive at

1 2 b t 2 2 + λ 2 d d t b 2 2 C ( u p p + b 2 2 ( 2 p 3 ) p 2 Δ b 2 2 ( 3 p ) p 2 + b 2 2 p p 2 ) ,

noting that p 5 2 , so 2(3p)/(p2)2, and now estimate (2.1), (2.11) and Gronwall’s inequality imply the estimate of b t in (2.12).

In the following, we will derive the bound for u t . Setting η= u t in (1.5), we deduce that

u t 2 2 + 1 p d d t e ( u ) p p | ( b ) b , u t | + | ( u ) u , u t | + χ | rot ω , u t | ε u t 2 2 + C ε ( ( b ) b 2 2 + ( u ) u 2 2 + ω 2 2 ) .

Integrating it over (0,t), t(0,T], by choosing ε= 1 2 and Korn’s inequality, we have

1 2 0 t u t 2 2 + C p u p p C 0 t ( ( b ) b 2 2 + ( u ) u 2 2 ) + C 1 .
(2.16)

Now, we compute, by (2.11)

0 t ( b ) b 2 2 0 t b 6 2 b 3 2 C 0 t b 2 2 ( b 2 Δ b 2 + b 2 2 ) C sup τ ( 0 , t ) b 2 3 0 t Δ b 2 + C sup τ ( 0 , t ) b 2 4 t C 3 4 ( 1 + T ) ,
(2.17)
0 t ( u ) u 2 2 0 t u 2 p p 2 2 u p 2 C 0 t u p 4 .
(2.18)

Inserting (2.17)-(2.18) into (2.16), by appealing to Korn’s inequality, it follows that

u p p + 0 t u t 2 2 C 0 t u p p u p 4 p + C 6 +C u 0 p p ,
(2.19)

where C 6 depends on T and C 3 . Now, Gronwall’s inequality and (2.1) yield that

sup t ( 0 , T ) u p + 0 T u t 2 2 C 7 ,
(2.20)

where C 7 depends on the time T, C 1 , C 3 and u 0 p .

3 Approximate solutions and existence result

In this section, we show the existence of a weak solution to the system (1.1)-(1.3) via the Galerkin approximations. For this purpose, we take the set { ϕ i } i = 1 formed by the eigenvectors ϕ i , i=1,2, , of the Stokes operator and the set { ψ i } i = 1 formed by the eigenvectors ψ i , i=1,2, , of the Laplace operator. According to the Appendix of [7], the functions { ϕ i } i = 1 form a basis in the space V p , and V 2 W 2 , 2 (Ω). Setting R k =span{ ϕ 1 , ϕ 2 ,, ϕ k } and S k =span{ ψ 1 , ψ 2 ,, ψ k }, we construct the Galerkin approximations { u k , ω k , b k } being of the form

u k (x,t)= i = 1 k a i k (t) ϕ i (x); ω k (x,t)= i = 1 k c i k (t) ψ i (x); b k (x,t)= i = 1 k d i k (t) ϕ i (x),

where a k :=( a 1 k ,, a k k ), c k :=( c 1 k ,, c k k ), d k :=( d 1 k ,, d k k ) solve the system of ordinary equations

u t k , ϕ + | e ( u k ) | p 2 e ( u k ) , e ( ϕ ) = χ rot ω k , ϕ + ( b k ) b k , ϕ ( u k ) u k , ϕ for all  ϕ R k ,
(3.1)
ω t k , ψ + μ ω k , ψ = χ rot u k , ψ ( u k ) ω k , ψ 2 χ ω k , ψ for all  ψ S k ,
(3.2)
b t k , ϕ + λ b k , ϕ = ( b k ) u k , ϕ ( u k ) b k , ϕ for all  ϕ R k .
(3.3)

Moreover, we require that u k , ω k , b k satisfy the following initial conditions

u k | t = 0 = i = 1 k ( u 0 , ϕ i ) ϕ i , ω k | t = 0 = i = 1 k ( ω 0 , ψ i ) ψ i , b k | t = 0 = i = 1 k ( b 0 , ϕ i ) ϕ i .
(3.4)

The local solvability is guaranteed by the Carathéodory theorem, and the global unique solvability follows from the fact that

sup t ( 0 , T ) u k ( t ) 2 2 = sup t ( 0 , T ) i = 1 k ( a i k ( t ) ) 2 C , sup t ( 0 , T ) ω k ( t ) 2 2 = sup t ( 0 , T ) i = 1 k ( c i k ( t ) ) 2 C , sup t ( 0 , T ) b k ( t ) 2 2 = sup t ( 0 , T ) i = 1 k ( d i k ( t ) ) 2 C

with upper bounds C that do not depend on k. Moreover, we have for u k , ω k and b k the same estimates for all norms we have obtained in Section 2. More precisely, we have

sup t ( 0 , T ) u k 2 , sup t ( 0 , T ) u k p , 0 T u t 2 2 C , sup t ( 0 , T ) ω k 1 , 2 , 0 t ω k 2 , 2 2 , 0 T ω t 2 2 C , sup t ( 0 , T ) b k 1 , 2 , 0 t b k 2 , 2 2 , 0 T b t 2 2 C
(3.5)

with a constant C that does not depend on k.

Uniform estimates (3.5) imply that there exists a subsequence of { u k }, { ω k } and { b k } (not relabeled) such that

u k u , weak-* in  L ( 0 , T ; H ) L ( 0 , T ; V p ) , ω k ω , weakly in  L 2 ( 0 , T ; W 2 , 2 ( Ω ) )  and  weak-* in  L ( 0 , T ; W 1 , 2 ( Ω ) ) , b k b , weakly in  L 2 ( 0 , T ; V 2 W 2 , 2 ( Ω ) )  and  weak-* in  L ( 0 , T ; V 2 ) , u t k u t , b t k b t , weakly in  L 2 ( 0 , T ; H ) , ω t k ω t , weakly in  L 2 ( 0 , T ; L 2 ( Ω ) ) , | e ( u k ) | p 2 e ( u k ) Λ , weakly in  L p ( 0 , T ; L p ( Ω ) ) ,

where p =p/(p1). Therefore, by making use of the Aubin-Lions lemma (see Lions [10], Theorem 1.5.1), we have

u k u , strongly in  L 2 ( 0 , T ; H ) , ω k ω , strongly in  L 2 ( 0 , T ; W 1 , 2 ( Ω ) ) , b k b , strongly in  L 2 ( 0 , T ; V 2 ) .

With the convergence above, it is easy to pass to the limit as k in (3.1)-(3.3) to find

u t , ϕ + Λ , e ( ϕ ) = χ rot ω , ϕ + ( b ) b , ϕ ( u ) u , ϕ for all  ϕ V p ,
(3.6)
ω t , ψ + μ ω , ψ = χ rot u , ψ ( u ) ω , ψ 2 χ ω , ψ for all  ψ W 2 , 2 ( Ω ) ,
(3.7)
b t , ϕ + λ b , ϕ = ( b ) u , ϕ ( u ) b , ϕ for all  ϕ V 2 .
(3.8)

Next, to complete the existence proof, we need to verify that

Λ= | e ( u ) | p 2 e(u).
(3.9)

By Lemma 1.3, we have

e ( u k ) e ( u ) p p C Ω ( | e ( u k ) | p 2 e ( u k ) | e ( u ) | p 2 e ( u ) ) ( e ( u k ) e ( u ) ) = C Ω | e ( u k ) | p 2 e ( u k ) e ( u k ) Ω | e ( u ) | p 2 e ( u ) ( e ( u k ) e ( u ) ) Ω | e ( u k ) | p 2 e ( u k ) e ( u ) = (3.1) χ rot ω k , u k + ( b k ) b k , u k u t k , u k Ω | e ( u ) | p 2 e ( u ) ( e ( u k ) e ( u ) ) Ω | e ( u k ) | p 2 e ( u k ) e ( u ) .

Considering lim k of this identity together with (3.6) implies that

lim k e ( u k ) e ( u ) p p 0,

and thus (3.9) follows.

Having the estimates

u L ( 0 , T ; W 0 1 , p ( Ω ) ) , ω L ( 0 , T ; W 0 1 , 2 ( Ω ) ) , b L ( 0 , T ; W 0 1 , 2 ( Ω ) ) C,
(3.10)

we can now introduce the pressure from (1.5). For t(0,T], define the functional F W 1 , p (Ω) as

F , ξ : = div ( | e ( u ) | p 2 e ( u ) ) , ξ + χ rot ω , ξ b × rot b , ξ ( u ) u , ξ u t , ξ .

We have

F,ξ=0,ξ V p , a.e. t(0,T].

By using De Rahm’s theorem (see [11], Lemma 2.7), we obtain a function π L p (Ω) such that

F=π,a.e. t(0,T].

Moreover, due to estimates (3.10),

π W 1 , p ( Ω ) C,a.e. t(0,T].

Then, by Lemma 1.2, there is a generic constant C, depending only on the data such that

π p C,a.e. t(0,T].

Now, we complete the proof of the existence part of Theorem 1.1.

4 Uniqueness criterion

Let ( u 1 , ω 1 , b 1 ) and ( u 2 , ω 2 , b 2 ) be both solutions of the problem. Then, for their difference u ¯ = u 1 u 2 , ω ¯ = ω 1 ω 2 , b ¯ = b 1 b 2 , we have

u ¯ t , ϕ + | e ( u 1 ) | p 2 e ( u 1 ) | e ( u 2 ) | p 2 e ( u 2 ) , e ( ϕ ) = χ rot ω ¯ , ϕ + ( b 1 ) b ¯ , ϕ + ( b ¯ ) b 2 , ϕ ( u 1 ) u ¯ , ϕ ( u ¯ ) u 2 , ϕ for all  ϕ V p ,
(4.1)
ω ¯ t , ψ + μ ω ¯ , ψ + 2 χ ω ¯ , ψ = χ rot u ¯ , ψ ( u 1 ) ω ¯ , ψ ( u ¯ ) ω 2 , ψ for all  ψ W 0 1 , 2 ( Ω ) ,
(4.2)
b ¯ t , η + λ b ¯ , η = ( b 1 ) u ¯ , η + ( b ¯ ) u 2 , η ( u 1 ) b ¯ , η ( u ¯ ) b 2 , η for all  η V 2 .
(4.3)

Taking ϕ= u ¯ in (4.1), by Lemma 1.3 and the fact that ( u 1 ) u ¯ , u ¯ =0 and p>2, we obtain

1 2 d d t u ¯ 2 2 + C p u ¯ p p χ rot ω ¯ , u ¯ + ( b ¯ ) b 2 , u ¯ ( u ¯ ) u 2 , u ¯ + ( b 1 ) b ¯ , u ¯ i = 1 3 J i + ( b 1 ) b ¯ , u ¯

for each J i , i=1,2,3, it follows from Hölder’s, Young’s inequality and (1.9), (1.10) that

J 1 χ | rot ω ¯ , u ¯ | χ ω ¯ 2 u ¯ 2 ε ω ¯ 2 2 + C ε u ¯ 2 2 , J 2 | ( b ¯ ) b 2 , u ¯ | b 2 6 b ¯ 3 u ¯ 2 C ( Δ b 2 2 + b 2 2 ) b ¯ 2 u ¯ 2 ε b ¯ 2 2 + C ε ( Δ b 2 2 2 + b 2 2 2 ) u ¯ 2 2 , J 3 | ( u ¯ ) u 2 , u ¯ | u 2 p u ¯ 2 p p 1 2 u 2 p u ¯ 2 2 p 3 p u ¯ 2 3 p ε u ¯ 2 2 + C ε u 2 p 2 p 2 p 3 u ¯ 2 2 ,

so, we have

1 2 d d t u ¯ 2 2 + C p u ¯ p p ε ( u ¯ 2 2 + ω ¯ 2 2 + b ¯ 2 2 ) + C ε ( 1 + b 2 2 2 + Δ b 2 2 2 + u 2 p 2 p 2 p 3 ) u ¯ 2 2 + ( b 1 ) b ¯ , u ¯ .
(4.4)

Taking ψ= ω ¯ in (4.2) and noting that ( u 1 ) ω ¯ , ω ¯ =0, it follows

1 2 d d t ω ¯ 2 2 + μ ω ¯ 2 2 + 2 χ ω ¯ 2 2 = χ rot u ¯ , ω ¯ ( u ¯ ) ω 2 , ω ¯ i = 1 2 I i ,

and for I i , i=1,2,

I 1 χ u ¯ 2 ω ¯ 2 ε u ¯ 2 2 + C ε ω ¯ 2 2 , I 2 ω 2 6 ω ¯ 3 u ¯ 2 C ( Δ ω 2 2 + ω 2 2 ) ω ¯ 2 u ¯ 2 ε ω ¯ 2 2 + C ε ( Δ ω 2 2 2 + ω 2 2 2 ) u ¯ 2 2 ,

thus, we obtain

1 2 d d t ω ¯ 2 2 + μ ω ¯ 2 2 ε ( u ¯ 2 2 + ω ¯ 2 2 ) + C ε ω ¯ 2 2 1 2 d d t ω ¯ 2 2 + μ ω ¯ 2 2 + C ε ( Δ ω 2 2 2 + ω 2 2 2 ) u ¯ 2 2 .
(4.5)

Similarly, by taking η= b ¯ in (4.3), reasoning analogous as above, we could get

1 2 d d t b ¯ 2 2 + λ b ¯ 2 2 ε b ¯ 2 2 + C ε ( Δ b 2 2 2 + b 2 2 2 ) u ¯ 2 2 + C ε u 2 p 2 p 2 p 3 b ¯ 2 2 + ( b 1 ) u ¯ , b ¯ .
(4.6)

Adding (4.4)-(4.6) and observing that ( b 1 ) b ¯ , u ¯ +( b 1 ) u ¯ , b ¯ =0, after choosing ε properly small, we finally get

1 2 d d t ( u ¯ 2 2 + ω ¯ 2 2 + b ¯ 2 2 ) + C ( u ¯ 2 2 + ω ¯ 2 2 + b ¯ 2 2 ) C F ( t ) ( u ¯ 2 2 + ω ¯ 2 2 + b ¯ 2 2 )

with

F(t)=1+ b 2 2 2 + Δ b 2 2 2 + u 2 p 2 p 2 p 3 + ω 2 2 2 + Δ ω 2 2 2 .

Since 2p/(2p3)p for p 5 2 , then Gronwall’s inequality and the estimates obtained in Section 2 yield that

u ¯ = ω ¯ = b ¯ =0for t[0,T].

This completes the proof of the theorem.

References

  1. Eringen AC: Theory of micropolar fluids. J. Math. Mech. 1966, 16: 1-18.

    MathSciNet  Google Scholar 

  2. Galdi GP, Rionero S: A note on the existence and uniqueness of solutions of the micropolar fluid equations. Int. J. Eng. Sci. 1977, 15: 105-108. 10.1016/0020-7225(77)90025-8

    Article  MathSciNet  Google Scholar 

  3. Rojas-Medar MA: Magneto-micropolar fluid motion: existence and uniqueness of strong solution. Math. Nachr. 1997, 188: 301-319. 10.1002/mana.19971880116

    Article  MathSciNet  Google Scholar 

  4. Ortega-Torres EE, Rojas-Medar MA: Magneto-micropolar fluid motion: global existence of strong solutions. Abstr. Appl. Anal. 1999, 4: 109-125. 10.1155/S1085337599000287

    Article  MathSciNet  Google Scholar 

  5. Rojas-Medar MA, Boldrini JL: Magneto-micropolar fluid motion: existence of weak solutions. Rev. Mat. Complut. 1998, 11: 443-460.

    Article  MathSciNet  Google Scholar 

  6. Gunzburger MD, Ladyzhenskaya OA, Peterson JS: On the global unique solvability of initial-boundary value problems for the coupled modified Navier-Stokes and Maxwell equations. J. Math. Fluid Mech. 2004, 6: 462-482.

    MathSciNet  Google Scholar 

  7. Málek J, Něcas J, Rokyta M, Ružička M: Weak and Measure-Valued Solutions to Evolutionary PDEs. Chapman & Hall, London; 1996.

    Book  Google Scholar 

  8. Bellout H, Bloom F, Nečas J: Solutions for incompressible non-Newtonian fluids. C. R. Math. Acad. Sci. Paris, Sér. I 1993, 317: 795-800.

    Google Scholar 

  9. Ebmeyer C, Urbano JM: Quasi-steady Stokes flow of multiphase fluids with shear-dependent viscosity. Eur. J. Appl. Math. 2007, 18: 417-434. 10.1017/S0956792507006948

    Article  MathSciNet  Google Scholar 

  10. Lions JL: Quelques méthodes de Résolution de Problèmes aux limites non Linéaires. Dunod, Paris; 1969.

    Google Scholar 

  11. Amrouche C, Girault V: Decomposition of vector spaces and application to the Stokes problem in arbitrary dimension. Czechoslov. Math. J. 1994, 44: 109-140.

    MathSciNet  Google Scholar 

Download references

Acknowledgements

The author would like to thank the referees for valuable comments and suggestions for improving this paper. This work was partially supported by the National NSF (Grant No. 10971080) of China.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Changjia Wang.

Additional information

Competing interests

The author declares that they have no competing interests.

Authors’ contributions

The author completed the paper. The author 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

Wang, C. Unique solvability for the non-Newtonian magneto-micropolar fluid. Bound Value Probl 2013, 182 (2013). https://doi.org/10.1186/1687-2770-2013-182

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords