Skip to main content

Global existence and convergence rates for the smooth solutions to the compressible magnetohydrodynamic equations in the half space

Abstract

Background

With the characteristics of low pollution and low energy consumption, the magnetohydrodynamics has made widely attention. This paper provides the standard energy method to solve the magnetohydrodynamic equations (MHD) in the half space + 3 . It proves the global existence for the compressible (MHD) by combining the careful a priori estimates and the local existence result. This study also considers the large time behaviors of the solutions.

Results

The interactions between the viscous, compressible fluid motion and the magnetic field are modeled by the magnetohydrodynamic system which describes the coupling between the compressible Navier-Stokes equations and the magnetic equations. This study has applied the analytical method to obtain the solutions to (MHD) in + 3 . It proves that under the assumption that the initial data are close to the constant state, the global existence of smooth solutions can be established. Moreover, the various decay rates of such solutions in Lp-norm with 2≤p≤+ and their derivatives in L2-norm can also be derived from combining the decay estimates of the linearized system and the energy method.

Conclusions

This study demonstrates that the global existence and the decay rates for the compressible (MHD) can be established under the similar initial assumptions as for the compressible Navier-Stokes equations. Especially, the results suggest that if the initial velocity is small, the velocity decays at a certain rate. This implies that only under the initial assumption that the data are large, it may reach the requirements of (MHD) power generation, which can be used to achieve the value of industrial application and environmental protection.

Background

Magnetohydrodynamics, which combines the environmental fluid mechanics and electrodynamics theories to study the interaction discipline between the conduction fluid and electromagnetic, is the theory of the macroscopic, and it has spanned a very large range of applications (Gerebeau et al. [2006]). Due to the lower environmental pollution, especially in energy industry, magnetohydrodynamic power generation is used to conserve energy and mitigate pollution in order to protect the environment. In virtue of the industrial importance and theoretical challenges, the study on (MHD) has attracted many scientists. In the present paper, we are interested in the well-posedness theory of (MHD). Many results concerning the existence and uniqueness of (weak, strong or smooth) solutions in one dimension can be found in (Chen and Wang [2002], [2003]; Kawashima and Okada [1982]) and the references cited therein. In multi-dimensional case the global existence of weak solutions for the bounded domains has been established recently in (Ducomet and Feireisl [2006]; Tan and Wang [2009]). The local unique strong solution has been obtained in (Fan and Yu [2009]). In (Chen and Tan [2010], [2012]) we established the global existence and decay rates of the smooth solutions for the Cauchy problem. However, many fundamental problems for the compressible (MHD) in the half space are still open. In this paper, we will extend our results (Chen and Tan [2010]) to the initial boundary problem in the half space.

Results

In this paper, we will consider the initial boundary value problem for the compressible magnetohydrodynamic equations (MHD) in the half space + 3 = x = ( x , x 3 ) : x 2 , x 3 > 0 (cf. (Gerebeau et al. [2006])):

ρ t + div ( ρ u ) = 0 , ( ρ u ) t + div ( ρ u u P ) = μ 0 curl H × H , H t curl ( u × H ) + 1 σ μ 0 curlcurl H = 0 , div H = 0 .
(1)

Here ρ,u=(u1,u2,u3),H=(H1,H2,H3) represent the density, velocity of the fluid and the magnetic field respectively. μ0>0 stands for permeability of free space, and σ>0 is the electric conductivity. The stress tensor P is given by

P = p I + μ u + u T + λ div u I ,
(2)

where p=p(ρ) is the pressure and the viscosity coefficients λ, μ satisfy

2 μ + 3 λ > 0 and μ > 0 .
(3)

For convenience, we reformulate the system (1) as

ρ t + div ( ρ u ) = 0 , ρ u t + ρ u · u + p = Δ u + div u + curl H × H , H t + curl ( u × H ) Δ H = 0 , div H = 0 ,
(4)

in (0,)× + 3 . Notice that we have normalized some physical constants to be unit but without reducing any essential difficulties for our analysis. We complement (2) the initial condition

(ρ,u,H)(0,x)=( ρ 0 (x), u 0 (x), H 0 (x)),x + 3 ,
(5)

and the following boundary conditions

u | { x 3 = 0 } =0,H | { x 3 = 0 } =0,
(6)

or

u | { x 3 = 0 } =0,H·n | { x 3 = 0 } =0,curlH×n | { x 3 = 0 } =0,
(7)

where n=(0,0,−1) is the normal vector of + 3 . We assume that throughout the paper the initial boundary data satisfy certain compatibility conditions as usual in (Matsumura and Nishida [1983]).

Before stating out our results, we shall introduce some standard notations.

Notations. We denote by Lp, Wm, p the usual Lebesgue and Sobolev spaces on + 3 and Hm=Wm,2, with norms |· | L p ,|· | W m , p ,|· | H m respectively. For the sake of conciseness, we do not precise in functional space names when they are concerned with scalar-valued or vector-valued functions. We denote = x =(1,2,3)t, where i = x i , and put x l f= l f= l 1 f for l=1,2,3,. We assume that C be a positive generic constant throughout this paper that may vary at different places and the integration domain + 3 will be always omitted without any ambiguity. Now our main results can be formulated as the following theorems. Firstly we state the results on the global existence and uniqueness of smooth solutions as:

Theorem 1.

Assume that the initial data are close enough to the constant state ( ρ ̄ ,0,0) with ρ ̄ >0, i.e., there exists a constant δ0 such that

| ρ 0 ρ ̄ , u 0 , H 0 | H 3 δ 0 .
(8)

Then there exists a unique globally smooth solution (ρ,u,H) of the initial boundary problem (2)–(4) or (2), (3) and (5) such that for any t[ 0,), it holds

| ( ρ ρ ̄ , u , H ) ( · , t ) | H 3 2 + 0 t | x ρ ( · , s ) | H 2 2 + | ( x u , x H ) ( · , s ) | H 3 2 ds C | ( ρ 0 ρ ̄ , u 0 , H 0 ) | H 3 2 .
(9)

Remark 1.

We will only prove Theorem 1 under the boundary condition (4). Due to the special geometry of the boundary of the half space, a simple calculation shows that the boundary conditions on H in (5) are equivalent to the following Dirichlet-Neumann boundary conditions:

3 H 1 | { x 3 = 0 } = 3 H 2 | { x 3 = 0 } = 0 , H 3 | { x 3 = 0 } = 0 .
(10)

Hence to treat H1,H2 as in (Matsumura and Nishida [1983]), we can also prove Theorem 1 under the boundary conditions (5) in the similar way as we will proceed.

By imposing some additional conditions on the initial data we will establish the following various decay rates of the solutions obtained in Theorem 1:

Theorem 2.

Let (ρ,u, H) be the solution obtained in Theorem 1 and assume in addition that the initial data ( ρ 0 ρ ̄ , u 0 , H 0 ) L 1 + 3 and there exists δ1>0 such that

|( ρ 0 ρ ̄ , u 0 , H 0 ) | L 1 +|( ρ 0 ρ ̄ , u 0 , H 0 ) | H 3 < δ 1 ,
(11)

then for all t≥0, it holds that

|(ρ ρ ̄ ,u,H)(t) | L p =O ( 1 + t ) 3 2 1 1 p ,p[2,],
(12)
| x (ρ ρ ̄ ,u,H)(t) | L 2 =O ( 1 + t ) 5 4 ,
(13)

and

|(ρ ρ ̄ ,u,H)(t) | H 3 C( δ 1 ) ( 1 + t ) 3 4 ε 1 ,0< ε 1 < ε ~ ,
(14)

where ε ~ is some positive number.

Moreover, if the data satisfy ( ρ 0 ρ ̄ , u 0 , H 0 ) H 4 + 3 and there exists δ2>0 such that

|( ρ 0 ρ ̄ , u 0 , H 0 ) | L 1 +|( ρ 0 ρ ̄ , u 0 , H 0 ) | H 4 < δ 2 ,
(15)

then

| x (ρ ρ ̄ ,u,H)(t) | H 3 C( δ 2 ) ( 1 + t ) 5 4 ε 2 ,0< ε 2 < ε ̂ ,
(16)

for all t≥0, where ε ̂ < ε ~ is a positive number. In fact, it holds

| x 2 ρ , x 2 u , x 2 H , x 3 H (t) | L 2 C( δ 2 ) ( 1 + t ) 5 4 ,
(17)

for all t≥0.

We will prove the global existence of smooth solutions by the standard energy method in spirit of (Matsumura and Nishida [1983], [1979], [1980]). And we remark that we can also obtain the global existence of strong solutions for the small initial data in the H2-framework, which can be proved in the similar way. On the other hand, the L2- L decay rates for the smooth solutions and the L2 decay rates for the derivatives of first order are optimal since (9)–(10) concerning (ρ ρ ̄ ,u) and H are the same as the optimal decay rates for the compressible Navier-Stokes equations (Kagei and Kobayashi [2005]) and the heat equation respectively. Related convergence rates of solutions for the Navier-Stokes equations on the unbounded domain can be found in (Kobayashi [2002]; Kobayashi and Shibata [1999]; Matsumura and Nishida [1979]) and the references cited therein. Although our proofs are in spirit of those for the Navier-Stokes equations, (Kagei and Kobayashi [2005]; Kobayashi [2002]; Kobayashi and Shibata [1999]; Matsumura and Nishida [1979], [1980]), we should derive the new estimates arising from the presence of the magnetic field and overcome the strong coupling between the mass, momentum equations and magnetic equation. However, it is easy to obtain the optimal decay rates of H and its first derivatives by the properties of heat kernel. Indeed, we can rewrite the equation (19) 3 analogously to the form of (19) 1–(19) 2, i.e.,

0 t + div H = 0 , H t + 0 Δ H div H = S 3 .
(18)

Thus, we can get the estimates of (0, H) which are similar to (ϱ,v) in the system (19) 1–(19) 2. Moreover, we will get the better decay rate of the magnetic field by the elliptic system.

Discussion

As well known, the heavy emissions of Greenhouse gases, such as C O2, C H4, N2O, S F6 cause global warming, and also result in a great deal of harm to the environment. It has been a hot topic and widespread concern to study on how to strictly control the greenhouse gases emissions. In order to profoundly reduce the environment pollution, we must focus on energy structure adjustment. Without any course of mechanical motion, Magnetohydrodynamics (MHD) power generation technology, also called plasma power generation technology, transforms thermal energy and kinetic energy directly into electricity. Thus by applying (MHD) power generation technology, we can realize the desulphuriz and reduce the production of N O x effectively, so as to achieve the effect of high efficiency and low pollution.

To complete the (MHD) generation process, which is of high industrial application value, a conductive gas (plasma) will be directed through a magnetic field with a large velocity, under a high temperature condition. In this situation, how to control the initial velocity of the conductive gas has to be considered. From the results in Section Results, we can conclude that if we assume that the initial data are close enough to the constant state, then there exists a unique globally solution to the (MHD) system and the solution decays at some rates. This indicates that if the initial velocity is sufficiently small, although the solution to the (MHD) system exists globally, then the velocity will decays and never be large, which implies that it may never reach the requirements of (MHD) power generation. However, the problem of the global existence of the solutions with the large initial data is still open.

Conclusions

In this paper, we demonstrate that the global existence and the decay rates for the compressible (MHD) in + 3 can be established under the similar initial assumptions as for the compressible Navier-Stokes equations which can be seen in (Matsumura and Nishida [1983]). It implies that the magnetic field does not affect the decay rates of the velocity. Indeed, the results (9)–(13) in Theorem 2 suggest that the decay rates for the derivatives of the magnetic field are the same as the velocity’s. And in (14), we cannot get the estimates for x 3 u but x 3 H. Furthermore, the results suggest that if the initial velocity is small, the velocity decays at the optimal rate. This implies that it may never reach the requirements of (MHD) power generation unless giving the gas an large initial velocity.

Methods

Proof of theorem 1

In this section, we will prove the existence part of Theorem 1 and the uniqueness is standard so it will be omitted.

Some elementary inequalities

The first bright idea to reduce many complicated computations lies in that we just need to do the lowest-order and highest-order energy estimates for the solutions. This is motivated by the following observation:

|f | H k 2 C| f , x k f | L 2 2 ,f H k + 3 .
(19)

The inequality (15) can be easily proved by combing Young’s inequality and Gagliardo-Nirenberg’s inequality

| x i f | L p C(p)|f | L q α | x k f | L r ( 1 α ) ,f H k + 3
(20)

where 1 p i 3 = 1 q α+ 1 r k 3 (1α) with ik. Indeed (16) can be proved by the extension technique together with Gagliardo-Nirenberg’s inequality in the whole space. We can obtain the following useful inequality by Hölder inequality and (16):

| x k ( fg ) | L 2 C | f | L | x k g | L 2 + | x k f | L 2 | g | L , f , g H k + 3 , k 2
(21)

and the general form of (3) can be deduced directly in the following:

| f 1 f 2 f s | k C j = 1 s | f 1 | L | f j 1 | L | f j | k | f j + 1 | L | f s | L , f j H k + 3 , k 2 , for j = 1 , 2 , , s.
(22)

The linearized system

We will linearize the problem (2)–(4) as follows. Setting γ= P ( ρ ̄ ) ,μ=1/ ρ ̄ and introducing new variables by

ϱ = ρ ρ ̄ , v = 1 μγ u , H = H .
(23)

Hence the initial boundary value problem (2)–(4) can be reformulated as

ϱ t + γ div v = S 1 , v t + γ ϱ μΔ v μ div v = S 2 , H t Δ H = S 3 , div H = 0 , ( ϱ , v , H ) ( x , 0 ) = ( ϱ 0 , v 0 , H 0 ) ( x ) = ( ρ 0 ρ ̄ , v 0 , H 0 ) , v | { x 3 = 0 } = H | { x 3 = 0 } = 0 ,
(24)

where

S 1 = μγ div ( ϱ v ) ,
(25)
S 2 = 1 μγρ curl H × H + 1 ρ 1 ρ ̄ Δ v + 1 ρ 1 ρ ̄ div v μγ v · v 1 μγ P ( ρ ) ρ P ( ρ ̄ ) ρ ̄ ϱ ,
(26)

and

S 3 = μγ curl ( v × H ) .
(27)

In order to state our results more concisely, we define an energy functional as:

N ( t 1 , t 2 ) = sup t 1 t t 2 | ( ϱ , v , H ) ( t ) | 3 2 + t 1 t 2 | x ϱ ( s ) | 2 2 + | ( x v , x H ) ( s ) | 3 2 dx 1 2 ,
(28)

and change the condition of the initial data (6) as

|( ϱ 0 , v 0 , H 0 ) | H 3 δ 0 =max 1 2 μγ , 1 δ 0 .
(29)

Local and global existence

We will finish the proof of Theorem 1 in this subsection. First we state out the local existence without proof, since it can be proved in a standard way (Matsumura and Nishida [1980]) or can be found in (Ströhmer [1990], Vol’pert and Hudjaev [1972]):

Theorem 3.

(local existence) Under the assumption (23), there exists a positive constant T such that the initial boundary value problem (19) has a unique solution (ϱ,v,H) which is continuous in [0,T]× 3 together with its derivatives of first order in t and of second order in x. Moreover, there exists a constant C1>1 such that it holds N(0,t)≤C1N(0,0), for any t[ 0,T].

We will prove in this subsection the following a priori estimate:

Theorem 4.

(a priori estimate) There exists a constant δ1 such that if N(0,T)≤δ, then there exists a constant C2>1 such that N(0,T)≤C2N(0,0).

The global existence of smooth solutions will be proved via a continued argument by combining the local existence theorem and the a priori estimate theorem. We shall state the global existence of smooth solutions to the linearized problem (19) as follows.

Proposition 3.1.

(global existence) Under the assumptions of Theorem 1, the initial boundary value problem (19) has a unique global solution such that for t[ 0,), it holds N(0,t)≤C N(0,0). Thus (ρ,u,H) which satisfies (4) uniquely solves the initial boundary value problem (2)–(4) for all time.

Proof.

See in (Chen and Tan [2010]). □

A priori estimates

We observe that the a priori assumption in Theorem 4 and the embedding inequality together with the continuity equation (19) 1 imply

sup 0 t T | ( ϱ , ϱ t , x ϱ , v , x v , H , x H ) ( t ) | C sup 0 t T | ( ρ ρ ̄ , v , H ) ( · , t ) | H 3 Cδ.
(30)

In particular

ρ ̄ 2 ρ = ϱ + ρ ̄ 2 ρ ̄ .
(31)

In the sequel, we will always use the smallness assumption of δ and (24)–(25).

Next we shall do some preparatory work from Lemma 3.1 to Lemma 3.7. Firstly, we regard the equations (19) 2–(19) 3 as the elliptic system with respect to x variables, i.e.,

Δ v + div v = 1 μ v t + γ μ ρ 1 μ S 2 , Δ H = H t S 3 , v | { x 3 = 0 } = 0 , H | { x 3 = 0 } = 0 .
(32)

Thus we have the following estimates which we can found in (Cho et al. [2004]):

Lemma 3.1.

Under the assumptions of Theorem 4, we have for k=2,3,4 that

| x k v | L 2 C | v t | k 2 + | x ϱ | k 2 + | S 2 | k 2 + | v | L 2 ,
(33)

and

| x k H | L 2 C | H t | k 2 + | S 3 | k 2 + | x H | L 2 .
(34)

Next we derive the following stokes equation from the equations (19) 1–(19) 2:

γ div v = g = dt μγϱ div v , μΔ v + γ ρ = h = v t μ γ g + S 2 , v | { x 3 = 0 } = 0 .
(35)

We have the following estimates which can be found in (Galdi et al. [1994]):

Lemma 3.2.

Under the assumptions of Theorem 4, we have for k=2,3,4 that

| x k v | L 2 2 + | x k 1 ρ | L 2 2 C | g | k 1 2 + | h | k 2 2 .
(36)

Next we shall do the estimates for the terms contained in N(0,t).

Lemma 3.3.

Under the assumptions of Theorem 4, we have that

| ( ϱ , v , H ) | L 2 2 + 0 t | x ( v , H ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 ,
(37)

and

| ( ϱ t , v t , H t ) | L 2 2 + 0 t | x ( v t , H t ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(38)
Proof.

By multiplying the equations (19) 1–(19) 3 by ϱ, v and H respectively, integrating over + 3 and adding the resulting equations, we have

1 2 d dt | ( ϱ , v , H ) | L 2 2 + μ | ( div , x v ) | L 2 2 + | x H | L 2 = S 1 , ϱ + S 2 , v + S 3 , H .
(39)

We shall estimate the terms in the right-hand side term by term as:

S 1 , ϱ = γμ ( ϱ v ) · x ϱdx C | ϱ | L 3 | v | L 6 | x ϱ | L 2 C | ϱ | H 1 | ( x v , x ϱ ) | L 2 2 | ( x v , x ϱ ) | L 2 2 , S 2 , v C | ( ϱ , v , H ) | L 3 | ( x ϱ , x v , x 2 v , x H ) | L 2 | v | L 6 C | ( ϱ , v , H ) | H 1 | ( x ϱ , x v , x 2 v , x H ) | L 2 2 | ( x ϱ , x v , x 2 v , x H ) | L 2 2 , S 3 , H C | ( v , H ) | L 3 | ( x v , x H ) | L 6 | H | L 2 C | ( v , H ) | H 1 | ( x v , x H ) | L 6 | H | L 2 | ( x v , x H ) | L 2 2 ,
(40)

where Hölder’s inequality and Sobolev’s inequality are used. Thus by integrating the above inequality in time and the definition of N(0,t), we have got (30).

Similarly, by taking t to (19) 1–(19) 3, multiplying by ϱ t , v t and H t respectively and integrating the resulting inequalities, we obtain

| ( ϱ t , v t , H t ) | L 2 2 + 0 t | x ( v t , H t ) | L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t t S 1 , ϱ t + t S 2 , v t + t S 3 , H t ds .
(41)

Thus we have to estimate the terms in the right-hand side of the above inequality.

x S 1 , ϱ t = γμ t ( ϱ div v + ϱ · v ) t ϱdx C | ( ϱ , x ϱ , v , x v ) | L 3 | ( ϱ t , x ϱ t , v t , x v t ) | L 2 | ϱ t | L 6 C | ( ϱ , x ϱ , v , x v ) | H 1 | ( ϱ t , x ϱ t , v t , x v t ) | L 2 2 | ( x ϱ , x 2 ϱ , x v , x 2 v , x 3 v , x H , x 2 H ) | L 2 2 CδN ( 0 , t ) 2 ,
(42)

where we can get the L2-norm estimates of ϱ t , x ϱ t , v t , and x v t with the aid of the equation (19) 1–(19) 2 and (24). Similarly we have that

t S 2 , v t C | t ( ϱ curl H × H , ϱΔ v , ϱ div v , v · v , ϱ · ϱ ) , v t | C | ( ϱ , ϱ t , v , x v , H , x H ) | L 3 × | ( x ϱ , x ϱ t , v t , x 2 v , x v t , x 2 v t , x H , H t , x H t ) | L 2 | v t | L 6 C | ( ϱ , x ϱ , v , x v , H , x H ) | H 1 × | ( x ϱ , x 2 ϱ , x 3 ϱ , x v , x 2 v , x 3 v , x 4 v , x H , x 2 H , x 3 H ) | L 2 2 | ( x ϱ , x 2 ϱ , x 3 ϱ , x v , x 2 v , x 3 v , x 4 v , x H , x 2 H , x 3 H ) | L 2 2 ,
(43)

and

t S 3 , H t C | t curl ( v × H ) , H t | C | ( v , x v , H , x H ) | L 3 | ( v t , x v t , H t , x H t ) | L 2 | H t | L 6 C | ( v , H ) | H 2 | ( x ϱ , x 2 ϱ , x v , x 2 v , x 3 v , x H , x 2 H , x 3 H ) | L 2 2 | ( x ϱ , x 2 ϱ , x v , x 2 v , x 3 v , x H , x 2 H , x 3 H ) | L 2 2 .
(44)

Together with these inequalities, we can deduce the inequality (31). Hence the proof of Lemma 3.3 is complete. □

Next we estimate the L2-norm of the first derivatives of v and H.

Lemma 3.4.

Under the assumptions of Theorem 4, we have

| x ( v , H ) | L 2 2 + 0 t | ( ϱ t , v t , H t ) ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 ,
(45)

and

| x ( v t , H t ) | L 2 2 + 0 t | ( ϱ tt , v tt , H tt ) ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(46)
Proof.

By multiplying ϱ t , v t and H t to the equations (19) 1–(19) 3 respectively, integrating over + 3 and adding the resulting equations, we have

d dt μ 2 | x v | 2 + μ 2 | div v | 2 + 1 2 | x H | 2 γϱ div v dx + ϱ t 2 + | v t | 2 + | H t | 2 + 2 γ ϱ t div v dx = S 1 , ϱ t + S 2 , v t + S 3 , H t dx ε | ϱ t | L 2 2 + | v t | L 2 2 + | H t | L 2 2 + C ( ε ) | S 1 | L 2 2 + | S 2 | L 2 2 + | S 3 | L 2 2 .
(47)

Since

| S 1 | L 2 2 C | ( ϱ , v ) | L | x ( ϱ , v ) | L 2 2 | x ( ϱ , v ) | L 2 2 ,
(48)
| S 2 | L 2 2 C | ( ϱ , v , H ) | L | ( x ϱ , x v , x 2 v , x H ) | L 2 2 | ( x ϱ , x v , x 2 v , x H ) | L 2 2 ,
(49)

and

| S 3 | L 2 2 C | ( v , H ) | L | x ( v , H ) | L 2 2 | x ( v , H ) | L 2 2 ,
(50)

thus by integrating (34) in time and with the aid of the smallness of ε and Cauchy’s inequality, we obtain

| x ( v , H ) | L 2 2 + 0 t | ( ϱ t , v t , H t ) ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C | ϱ | L 2 2 + 0 t | x v ( s ) | L 2 2 ds ,
(51)

which together with (30) yields (32). And we can also get (33) in the similar way. The proof of Lemma 3.4 is complete. □

Next we estimate the L2-Norms of the first derivatives of ϱ. We shall divide the estimates into two parts. Firstly we denote the tangential derivatives by =(1,2). And it is easy to see that the tangential derivatives of the solution of (19) satisfy the same boundary conditions in (19).

By taking to the equations (19) 1–(19) 2 and multiplying them with ϱ and v respectively, we can deduce that

1 2 d dt | ∂ϱ | 2 + | v | 2 dx + μ | x v | 2 + μ | div v | 2 dx = S 1 , ∂ϱ + S 2 , v .
(52)

Define the material derivative dt = ρ t +ρ·u. Then by the continuity equation (19) 1 and the formula of variable substitution (18), we have

Dt = γ div v μγϱ div v ,
(53)

where Dt = dt = ϱ t +μγϱ·v. By taking to the equation (36), we can obtain the following inequality:

Dt L 2 2 C x v | L 2 2 + | ( ϱ div v ) | L 2 2 ,
(54)

then by multiplying a small enough constant α to it, together with (35), and integrating in time, we have

| ( ϱ , v ) | L 2 2 + 0 t Dt , x v L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t S 1 , ∂ϱ + S 2 , v ds + C 0 t | ( ϱ div v ) | L 2 2 ds CN ( 0 , 0 ) 2 + CN ( 0 , t ) 2 + C 0 t S 1 , ∂ϱ + S 2 , v ds .
(55)

Similarly, we can obtain the following lemma and we omit the proof of it.

Lemma 3.5.

Under the assumptions of Theorem 4, we have for k = 1, 2, 3 that

| k ( ϱ , v ) | L 2 2 + 0 t k Dt , x v L 2 2 ds CN ( 0 , 0 ) 2 + CN ( 0 , t ) 2 + C 0 t k S 1 , k ϱ + k S 2 , k v ds .
(56)

Next we have to obtain the estimates for the normal derivatives of solution. We can derive the following equations from (19) 1–(19) 2.

3 ϱ t + γ 3 div v = μγ 3 ( ϱ div v ) μγ 3 ( ϱ · v ) , v t 3 + γ 3 ϱ μΔ v 3 μ 3 div v = S 2 3 ,
(57)

where S 2 = S 2 1 , S 2 2 , S 2 3 T . To eliminate the term 33v3, we have the following equality from the above equalities

2 μ γ 3 ϱ t + γ 3 ϱ = v t 3 + μ 11 v 3 + 22 v 3 + 13 v 1 + 23 v 2 2 μ 2 3 ( ϱ div v ) 2 μ 2 3 ( ϱ · v ) + S 2 3 .
(58)

And we can also derive the following equality from (36)

2 μ γ 3 Dt + γ 3 ϱ = v t 3 + μ 11 v 3 + 22 v 3 + 13 v 1 + 23 v 2 2 μ 2 3 ( ϱ div v ) + S 2 3 .
(59)

By multiplying (38) with 3ϱ and integrating on + 3 , we have

μ γ d dt | 3 ϱ | L 2 2 + γ | 3 ϱ | L 2 2 = v t 3 + μ 11 v 3 + 22 v 3 + 13 v 1 + 23 v 2 2 μ 2 3 ( ϱ div v ) 2 μ 2 3 ( ϱ · v ) + S 2 3 3 ϱdx γ 2 | 3 ϱ | L 2 2 + C ( γ ) | v t , x v , 3 ( ϱ div v ) , S 2 | L 2 2 + C 3 ( ϱ · v ) 3 ϱdx ,
(60)

then integrating this in time, it implies

| 3 ϱ | L 2 2 + 0 t | 3 ϱ ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t | ( v t , x v , 3 ( ϱ div v ) , S 2 ) | L 2 2 + | 3 ( ϱ · v ) 3 ϱdx | ds CN ( 0 , 0 ) 2 + C 0 t | ( v t , x v ) | L 2 2 ds + CδN ( 0 , t ) 2 .
(61)

Similarly, by multiplying 3 Dt to (39) and integrating it, we have

| 3 ϱ | L 2 2 + 0 t | 3 Dt ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t | ( v t , x v , 3 ( ϱ div v ) , S 2 ) | L 2 2 CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t | ( v t , x v ) | L 2 2 .
(62)

This together with (41) yields

| 3 ϱ | L 2 2 + | 3 ϱ ( s ) | L 2 2 + | 3 Dt ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t | ( v t , x v ) | L 2 2 ds.
(63)

Similarly, by taking k 3 l to (38) and (39), multiplying with k 3 1 + l ϱ and k 3 1 + l Dt respectively, we will get the general form of (42) as follows:

Lemma 3.6.

Under the assumptions of Theorem 4, we have for k+l=0,1,2 that

| k 3 1 + l ϱ | L 2 2 + 0 t | k 3 1 + l ϱ | L 2 2 + | k 3 1 + l Dt ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t | k 3 l ( v t , x v ) | L 2 2 ds.
(64)
Proof.

Here we only to prove the case when k+l=2. As in (42), we can obtain

| k 3 1 + l ϱ | L 2 2 + 0 t | k 3 1 + l ϱ | L 2 2 + | k 3 1 + l Dt ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t | k 3 l ( v t , x v , 3 ( ϱ div v ) , S 2 ) | L 2 2 + | k 3 1 + l ( ϱ · v ) · k 3 1 + l ϱdx | ds.
(65)

By (17), we have

0 t | x 3 ( ϱ div v ) | L 2 2 ds C 0 t | ϱ , div v | L 2 | x 3 ϱ , x 3 div v | L 2 2 ds CδN ( 0 , t ) 2 .
(66)

Similarly,

0 t | x 2 S 2 | L 2 2 ds C 0 t | ϱ , v, H | L 2 | x 3 ϱ , x 3 v, H | L 2 2 ds + C 0 t | ϱ , x ϱ | L 2 | x 2 v | 1 2 + | x 2 ϱ | L 4 2 | x 2 v | L 4 2 ds CδN ( 0 , t ) 2 + C 0 t | x 2 ϱ | 1 2 | x 2 v | 1 2 ds CδN ( 0 , t ) 2 .
(67)

Moreover, by using the identity

( k 3 1 + l ϱ · v ) · k 3 1 + l ϱ = | k 3 1 + l ϱ | 2 2 · v
(68)

and integration by parts, we shall obtain

0 t k 3 1 + l ( ϱ · v ) · k 3 1 + l ϱdxds CδN ( 0 , t ) 2 .
(69)

Substituting these inequalities into (44), we get (43). □

Last by taking x l k to (28), and by Lemma 3.2, we have

Lemma 3.7.

Under the assumptions of Theorem 4, we have for k+l=0,1,2 that

0 t | x 2 + l k v | L 2 + | x 1 + l k ϱ | L 2 ds C 0 t | k v t | l + | k Dt | 1 + l ds + CδN ( 0 , t ) 2 .
(70)

Now we will finish the proof of Theorem 4 by doing the estimates for the lowest-order and highest-order derivatives. In the sequel, we divide the a priori estimates into three parts.

Part 1: estimates for the lowest derivatives of ϱ,v,H

Proposition 3.2.

Under the assumptions of Theorem 4, we have

| ( ϱ , v , H ) | L 2 2 + 0 t | x ( ϱ , v , H ) ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(71)
Proof.

Due to (30), we have only to estimate 0 t | x ϱ(s) | L 2 2 ds. First by Lemma 3.5, k=0,1, we have

0 t | Dt ( s ) | L 2 2 + | Dt ( s ) | L 2 2 + | x v | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + 0 t | ( S 1 , S 2 , S 3 ) ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(72)

And by (42), we have

0 t | 3 Dt | L 2 2 ds CN ( 0 , 0 ) 2 + C 0 t | ( v t , x v ) | L 2 2 ds + CδN ( 0 , t ) 2 .
(73)

Thus (47)–(48) together with (32) implies

0 t | Dt ( s ) | 1 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(74)

By Lemma 3.4, k=l=0, and integrating in time, we obtain

0 t | x ϱ | L 2 2 ds C 0 t | v t | L 2 2 + | Dt | 1 2 ds + CδN ( 0 , t ) 2 .
(75)

This combining with (49) and (32), yields the estimate for 0 t | x ϱ(s) | L 2 2 ds. □

Part 2: estimates for the highest derivatives of ϱ,v,H

Proposition 3.3.

Under the assumptions of Theorem 4, we have

| x 3 ( ϱ , v , H ) | L 2 2 + 0 t | x 3 ϱ , x 4 v , x 4 H ( s ) | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(76)
Proof.

We divide the proof into three steps as follows.

i) By Lemma 2.8, k=3 and Lemma 2.9, k=2,l=0, we have

2 x ϱ L 2 2 + 0 t 2 3 ϱ , 2 x Dt , 3 x v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t k S 1 , k ϱ + k S 2 , k v ds + C 0 t | 2 v t | L 2 2 CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t | 2 v t | L 2 2 ds ,
(77)

where we get the last inequality as in the proof of Lemma 3.6. By Lemma 3.7, k=2,l=0, we have

0 t x 2 2 v , x 2 ϱ L 2 2 ds C 0 t 2 v t , 2 x Dt L 2 2 ds + CδN ( 0 , t ) 2 .
(78)

This together with (51) yields

2 x ϱ L 2 2 + 0 t 2 x ϱ , 2 x Dt , 2 x 2 v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t 2 v t L 2 2 ds.
(79)

By Lemma 3.6, k=l=1, and (52), we have

x 2 ϱ L 2 2 + 0 t x 2 ϱ , x 2 Dt , 2 x 2 v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t L 2 2 ds.
(80)

Then by Lemma 3.7, k=l=1, and with the help of (15), (31) and (47), we have

0 t x 3 v , x 2 ∂ϱ L 2 2 ds C 0 t x v t , x 2 v t , dt , x dt , x 2 Dt L 2 2 ds + CδN ( 0 , t ) 2 CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t , x 2 Dt L 2 2 ds ,
(81)

which together with (53) yields

x 2 ϱ L 2 2 + 0 t x 2 ϱ , x 2 Dt , x 3 v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t L 2 2 ds.
(82)

By Lemma 3.6, k=0,l=2, and by (54), we have

x 3 ϱ L 2 2 + 0 t x 3 ϱ , x 3 Dt , x 3 v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t L 2 2 ds.
(83)

Thus by Lemma 3.7, k=0,l=2, we have

0 t x 4 v , x 3 ϱ L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t , x 3 Dt L 2 2 ds ,
(84)

together it with (55), then we have got

x 3 ϱ L 2 2 + 0 t x 3 ϱ , x 3 Dt , x 4 v L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t L 2 2 ds.
(85)

i i) By Lemma 3.1, k=3, we have

x 3 ( v , H ) L 2 2 C x ϱ , v t , H t 1 2 + ( S 2 , S 3 ) 1 2 + v , x H L 2 2 CδN ( 0 , t ) 2 + C ϱ , x 3 ϱ , v , v t , x v t , x H , H t , x H t L 2 2 ,
(86)

which together with (30)–(33), (46) and (56) implies

x 3 v , H L 2 2 CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t L 2 2 ds.
(87)

i i i) By Lemma 3.1, k=4, we have

0 t x 4 H L 2 2 ds C 0 t H t 2 2 + | S 3 | 2 2 + | x H | L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + 0 t x 2 H t L 2 2 ds ,
(88)

this together with (56) and (58), then it implies

x 3 ( ϱ , v , H ) L 2 2 + 0 t x 3 ϱ , x 4 v , x 4 H ( s ) L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 + C 0 t x 2 v t , H t L 2 2 ds.
(89)

We can derive from Lemma 3.1 that

x 2 ( v t , H t ) L 2 2 C ( x ϱ t , v t , v tt , x H t , H tt , t S 2 , t S 3 L 2 2 .
(90)

Here with the help of (19) 1, we have

0 t x ϱ t L 2 2 ds C 0 t | x 2 ϱ | L 2 2 + x 2 v L 2 2 ds C ( ε ) 0 t x ϱ , x v L 2 2 ds + ε 0 t | x 3 ϱ , x 4 v | L 2 2 ds ,
(91)

where the terms in the right-hand side can be absorbed by (46) and (59). Thus by (30)–(33), we obtain

0 t x 2 ( v t , H t ) L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 ,
(92)

which together with (59) implies (50). □

Part 3: ConclusionCombining the inequalities in Proposition 3.2 and Proposition 3.3, it yields

ϱ , x 3 ϱ , v , x 3 v , H , x 3 H L 2 2 + 0 t x ϱ , x 3 ϱ , x v , x 4 v , x H , x 4 H L 2 2 ds CN ( 0 , 0 ) 2 + CδN ( 0 , t ) 2 .
(93)

Thanks to (15), we can obtain that the left-hand side of (60) is equivalent to N(0,t)2. Thus by the smallness of δ, it finished the proof of Theorem 4. □

Proof of Theorem 2

In this section we shall prove the decay rates of the solution obtained in Theorem 1 to finish the proof of Theorem 2.

Some elementary decay-in-time estimates

We shall consider the convergence rates of the solution (ϱ,v,H) for the linearized problem (19). To use the Lp- Lq estimates of the linear problem for the nonlinear problem (19) 1–(19) 2, we rewrite the solution of (19) 1–(19) 2 as

U ( t ) = E ( t ) U 0 + 0 t E ( t s ) F ( U ( s ) , H ( s ) ) ds ,
(94)

where we have used the notations

U = [ ϱ , v ] T , U 0 = [ ϱ 0 , v 0 ] T , F = [ S 1 , S 2 ] T ,
(95)

and the fact that E(t) is the solution semigroup defined by E(t)=etA, t≥0, with A being a matrix-valued differential operator given by

A = 0 γ T γ μΔ μ div .
(96)

The semigroup E(t) has the following decay-in-time properties which can be found in (Kagei and Kobayashi [2005]; Kobayashi [2002]).

Lemma 3.8.

There exist positive constants C and C0 such that for any t≥1 and l=1,2, we havei)

E ( t ) U 0 L 2 C t 3 4 U 0 L 1 + e C 0 t ϱ 0 | L 2 + | v 0 L 2 ,
(97)
x l E ( t ) U 0 L 2 C t 3 4 l 2 | U 0 | L 1 + e C 0 t ( | ϱ 0 | l + | v 0 | l 1 ) ,
(98)
x 3 E ( t ) U 0 L 2 C t 15 8 | U 0 | L 1 + e C 0 t ( | ϱ 0 | 3 + | v 0 | 2 ) ,
(99)

and

E ( t ) U 0 L C t 3 2 | U 0 | L 1 + e C 0 t ( | ϱ 0 | 2 + | v 0 | 1 ) .
(100)

ii)

x l t E ( t ) U 0 L 2 C t 5 4 l 2 | U 0 | L 1 + e C 0 t ( | ϱ 0 | l + | v 0 | l 1 ) .
(101)
Lemma 3.9.

For any k, k≥0, the following inequalities hold uniformly in 0<t≤1,

E ( t ) U 0 H 1 L 2 C | U 0 | H 1 L 2 ,
(102)
x k x E ( t ) U 0 L 2 C t 1 2 | U 0 | H k + 1 H k ,
(103)

and

E ( t ) U 0 L C t ( 1 ε ̄ ) | U 0 | H 2 H 1 ,
(104)

for some ε ̄ >0. Here |·|XY=|·| X +|·| Y .

To treat the magnetic field, we notice that the solution to the heat equation (19) 3 has the following convergence estimates which one can refer to (Kagei and Kobayashi [2005]; Kobayashi [2002]).

Lemma 3.10.

For the solution H to the heat equation (19) 3 with the initial data H(x,0)=H0 and the boundary condition H | { x 3 = 0 } =0, there exists a constant C such that

| H | L q C ( 1 + t ) 3 2 1 1 q | H 0 | L 1 L q + C 0 t ( 1 + t s ) 3 2 ( 1 1 q ) | S 3 ( · , s ) | L 1 L q ds ,
(105)

and

| x H | L 2 C ( 1 + t ) 5 4 | H 0 | L 1 H 1 + C t 1 t ( t s ) 1 2 | S 3 | L 2 ds + C 0 t 1 ( 1 + t s ) 5 4 | S 3 ( · , s ) | L 1 ds ,
(106)

for any t≥1, 2≤q≤+.

Convergence rates of the lower-order derivatives

Proposition 3.4.

Let s≥3. Under the assumptions of Proposition 3.1, if there exists δ 1 >0 such that the initial data (ϱ0,v0,H0)HsL1 and

| ( ϱ 0 , v 0 , H 0 ) | H s L 1 δ 1 ,
(107)

then the solution (ϱ,v,H) of (19) satisfies

| ( ϱ , v , H ) ( t ) | L p = O t 3 2 ( 1 1 p ) , p [ 2 , ] .
(108)

and

| x ( ϱ , v , H ) ( t ) | L 2 = O t 5 4
(109)

as t.

Proof.

We firstly estimate the L2-estimate of the solution and its derivatives in Proposition 3.4. For simplicity, we shall introduce some notation in the sequel. Set

M σ s ( t ) = sup 0 τ t ( 1 + τ ) σ | x s ( ϱ , v , H ) ( τ ) | L 2 ,
(110)

We will show that

M 3 4 0 ( t ) + M 5 4 1 ( t ) C | ( ϱ 0 , v 0 , H 0 ) | L 1 + N ( 0 , t ) + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(111)

where t≥0. Thus we can derive the L2 estimates in Proposition 3.4 from (73) by using a standard method under assumption of the smallness of initial data.

Now we shall consider the L2 estimates of U(t) and H(t) respectively. Thanks to Proposition 3.1, we only have to show the decay rate part of Proposition 3.4 for the case t≥1. Thus we assume that t≥1 and decompose U(t) as

U ( t ) = E ( t ) U 0 + t 1 t E ( t s ) F ( U ( s ) , H ( s ) ) ds + 0 t 1 E ( t s ) F ( U ( s ) , H ( s ) ) ds = I 0 ( t ) + I 1 ( t ) + I 2 ( t ) .
(112)

By Lemma 3.8, we have

x l I 0 ( t ) L 2 C t 3 4 l 2 U 0 L 1 H 1 C ( 1 + t ) 3 4 l 2 × | U 0 | L 1 + N ( 0 , t ) , l = 0 , 1 .
(113)

Lemma 3.11.

Under the assumptions of Proposition 3.4, the following inequalities hold for all t≥0,

| S 1 , x S 1 , S 2 , S 3 | L 2 C ( 1 + t ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 ,
(114)

and

| S 1 , S 2 , S 3 | L 1 C ( 1 + t ) 11 8 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 .
(115)
Proof.

Since by (16) and (17), we have

| S 1 , S 3 | L 2 C | ( ϱ , v , H ) | L | x ( ϱ , v , H ) | L 2 C ( 1 + t ) 5 4 M 5 4 1 ( t ) N ( 0 , t ) C ( 1 + t ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 ,
(116)
| x S 1 | L 2 C | x ( ϱ , v ) | L | x ( ϱ , v ) | L 2 + ϱ , v L 6 x 2 ( ϱ , v ) L 6 1 2 x 2 ( ϱ , v ) L 2 1 2 C ( 1 + t ) 5 4 M 5 4 1 ( t ) N ( 0 , t ) + | x ( ϱ , v ) | L 2 x 2 ϱ , x 3 ϱ , x 2 v , x 3 v L 2 C ( 1 + t ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 ,
(117)

and

| S 2 | L 2 C | ( ϱ , v , H ) | L | x ( ϱ , v , H ) | L 2 + C x 2 v L 2 1 2 x 2 v L 6 1 2 ϱ L 6 C ( 1 + t ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 ,
(118)

together with these inequalities, we can easily get (76). Moreover, (77) follows from

| S 1 , S 3 | L 1 C | ( ϱ , v , H ) | L 2 | x ( ϱ , v , H ) | L 2 C ( 1 + t ) 2 M 3 4 0 ( t ) M 5 4 1 ( t ) C ( 1 + t ) 2 × M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(119)

and

| S 2 | L 1 C | ( ϱ , v , H ) | L 2 | ( x ϱ , x v , x H ) | L 2 + C | ϱ | L 2 | x 2 v L 2 C ( 1 + t ) 2 M 3 4 0 ( t ) M 5 4 1 ( t ) + C ϱ L 2 x v L 2 1 2 x 3 v | L 2 1 2 C ( 1 + t ) 2 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + C ( 1 + t ) 11 8 M 3 4 0 ( t ) M 5 4 1 ( t ) 1 2 N 1 2 ( 0 , t ) C ( 1 + t ) 11 8 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 .
(120)

Now we have to estimate I1(t) and I2(t). For l=0,1, and by Lemma 3.9 and Lemma 3.11, we have

| x l I 1 ( t ) | L 2 C t 1 t ( t τ ) l 2 ( | S 1 | 1 + | S 2 | L 2 ) C t 1 t ( t τ ) l 2 ( 1 + τ ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 C ( 1 + t ) 5 4 M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 ,
(121)

where the last inequality is obtained by taking l=1 and l=2 respectively.

For I2(t), we can derive from Lemma 3.8 and Lemma 3.9 that

| x l I 2 ( t ) | L 2 C 0 t 1 ( t τ ) 3 2 l 2 | F ( U ( τ ) , H ( τ ) ) | L 1 + e C 0 ( t τ ) | U ( τ ) | 1 C 0 t 1 ( t τ ) 3 2 l 2 ( 1 + τ ) 11 8 × M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 + C 0 t 1 e C 0 ( t τ ) ( 1 + τ ) 5 4 × M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 C ( 1 + t ) 3 4 l 2 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 .
(122)

Now we turn to do the estimates for H(t). By (71) in Lemma 3.10, we can deduce that

| H | L 2 C ( 1 + t ) 3 4 | H 0 | L 1 L 2 + 0 t ( 1 + t s ) 3 4 | S 3 | L 1 L 2 ds C ( 1 + t ) 3 4 | H 0 | L 1 + N ( 0 , t ) + C M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + N ( 0 , t ) 2 × 0 t ( 1 + t s ) 3 4 ( 1 + s ) 11 8 ds C ( 1 + t ) 3 4 | H 0 | L 1 + N ( 0 , t ) + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(123)

and similarly by (72), we have

| x H | L 2 C ( 1 + t ) 5 4 | H 0 | L 1 H 1 + t 1 t ( t s ) 1 2 | S 3 | L 2 ds + 0 t 1 ( 1 + t s ) 5 4 | S 3 | L 1 ds C ( 1 + t ) 5 4 | H 0 | L 1 + N ( 0 , t ) + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 .
(124)

Combining these with (75), (78) and (79), it implies (73).

Now we shall do the L estimate. Define

M ( t ) = sup 0 τ t ( 1 + τ ) 3 2 | ( ϱ , v, H ) ( τ ) | L ,
(125)

and

σ s ( t ) = k = 0 s M σ k ( t ) = k = 0 s sup 0 τ t ( 1 + τ ) σ | x k ( ϱ , v, H ) ( τ ) | L 2 .
(126)

By Lemma 3.8, we have

| I 0 | L C t 3 2 ( | U 0 | L 1 + | ϱ 0 | 2 + | v 0 | 1 ) C ( 1 + t ) 3 2 ( | U 0 | L 1 + N ( 0 , t ) ) .
(127)

To estimate I1(t) and I2(t), we state the following estimates for S1 and S2.

Lemma 3.12.

Under the assumptions of Proposition 3.4, there exists a sufficiently small constant ε ~ >0 such that for any 0<ε< ε ~ , the following inequalities hold uniformly in t≥0.i)

| S 1 | 2 , | S 2 | 1 C ( 1 + t ) 3 2 M ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 ,
(128)

and

| S 1 , S 2 | L 1 C t ( 7 4 ε ) M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(129)

ii)

| S 1 | 3 , | S 2 | 2 C ( 1 + t ) 3 2 M ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 4 ( t ) 2 ,
(130)

and

| t S 1 | 1 , | t S 2 | L 2 C ( 1 + t ) 5 4 K 3 4 ε 4 ( t ) 2 + M ( t ) 2 + M 5 4 1 ( t ) 2 ,
(131)

here we define

K σ s ( t ) = σ s ( t ) + j = 1 [ s 2 ] sup 0 τ t ( 1 + τ ) σ | τ j ( v, H ) ( τ ) | s 2 j + j = 1 [ s + 1 2 ] sup 0 τ t ( 1 + τ ) σ | τ j ϱ ( τ ) | s + 1 2 j .
(132)
Proof.

For l=0,1,2, we have

| x l S 1 | L 2 C | x l + 1 ( ϱ v ) | L 2 C | ϱ , v | L | x ( ϱ , v ) | 2 + | x ( ϱ , v ) | L 2 1 2 | x 2 ( ϱ , v ) | 1 3 2 C ( 1 + t ) 9 4 ε M ( t ) 3 4 ε 3 ( t ) + ( 1 + t ) 7 4 3 2 ε M 5 4 1 ( t ) 1 2 3 4 ε 3 ( t ) 3 2 C ( 1 + t ) 3 2 M ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 ,
(133)

and for l=0,1, we have

| x l S 2 | L 2 C | ϱ , v , H | L | x ( ϱ , v , H ) | 2 + | x ( ϱ , v , H ) | L 2 1 2 | x ( ϱ , v , H ) | 1 3 2 C ( 1 + t ) 3 2 M 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(134)

Then we get (81). Similarly,

| S 1 | L 1 C | ϱ , v | L 2 | x ( ϱ , v ) | L 2 C ( 1 + t ) 2 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(135)

and

| S 2 | L 1 C | ϱ , v , H | L 2 | x ( ϱ , v , H ) | L 2 + | ϱ | L 2 | x 2 v | L 2 C | ϱ , v , H | L 2 | x ( ϱ , v , H ) | L 2 + | ϱ | L 2 | x v | L 2 1 2 | x 3 v | L 2 1 2 C ( 1 + t ) 2 M 3 4 0 ( t ) M 5 4 1 ( t ) + ( 1 + t ) 7 4 ε M 3 4 0 ( t ) M 5 4 1 ( t ) 3 4 ε 3 ( t ) C ( 1 + t ) ( 7 4 ε ) M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(136)

Thus we obtain (82). We shall get the inequalities in i i) in the similar way. Here we only estimate | x 3 S 1 | L 2 and | t (ϱΔv) | L 2 as follows.

| x 3 S 1 | L 2 C | ϱ , v | L | x 4 ( ϱ , v ) | L 2 C ( 1 + t ) 9 4 ε M ( t ) K 3 4 ε 4 ( t ) C ( 1 + t ) 3 2 M 2 + K 3 4 ε 4 ( t ) 2 ,
(137)

and

| t ( ϱΔ v ) | L 2 | ϱ | L | Δ v t | L 2 + | ϱ t | 1 | x 3 v | L 2 C ( 1 + t ) 9 4 ε M ( t ) K 3 4 ε 4 ( t ) + ( 1 + t ) 3 2 2 ε K 3 4 ε 4 ( t ) 2 C ( 1 + t ) 5 4 M 2 + K 3 4 ε 4 ( t ) 2 .
(138)

By Lemma 3.9 and (81), we have

| I 1 | L C t 1 t ( t τ ) ( 1 ε ̄ ) ( | S 1 | 2 + | S 2 | 1 ) C M 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 × t 1 t ( t τ ) ( 1 ε ̄ ) ( 1 + τ ) 3 2 C ( 1 + t ) 3 2 M 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(139)

As for I2, we can deduce from Lemma 3.8 that

| I 2 ( t ) | L C 0 t 1 ( t τ ) 3 2 | ( S 1 , S 2 ) ( τ ) | L 1 + 0 t 1 e C 0 ( t τ ) ( | S 1 ( τ ) | 2 + | S 2 ( τ ) | 1 ) C M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 × 0 t 1 ( t τ ) 3 2 ( 1 + τ ) 7 4 ε 1 + C M 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 × 0 t 1 e C 0 ( t τ ) ( 1 + τ ) 3 2 C ( 1 + t ) 3 2 M 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(140)

At last, we have to estimate |H | L to finish the L estimate. By Lemma 3.3, For q=, we have

| H | L C ( 1 + t ) 3 2 | H 0 | L 1 L + C 0 t ( 1 + t τ ) 3 2 | S 3 ( τ ) | L 1 L C ( 1 + t ) 3 2 | H 0 | L 1 + N ( 0 , t ) + C M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + M ( t ) 2 + N ( 0 , t ) 2 × 0 t ( 1 + t τ ) 3 2 ( 1 + τ ) 3 2 C ( 1 + t ) 3 2 | H 0 | L 1 + N ( 0 , t ) + M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(141)

where we have used

| S 3 | L 1 C ( 1 + t ) 2 M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 ,
(142)

which is in the proof of Lemma 3.12, and

| S 3 | L | v, H | L | v, H | 3 C ( 1 + t ) 3 2 M ( t ) 2 + N ( 0 , t ) 2 .
(143)

Thus, (80) and (85)–(87) imply

M C | ϱ 0 , v 0 , H 0 | L 1 + N ( 0 , t ) + M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + 3 4 ε 3 ( t ) 2 .
(144)

This together with (73), the estimate of K 3 4 ε 4 (t) in Proposition 3.5 in the next subsection, it yields the L-estimate in Proposition 3.4. Hence by interpolation, we get the Lp-estimate for 2≤p. The proof of Proposition 3.4 is complete. □

Convergence rates of the higher-order derivatives

Now we shall do the estimates for the higher-order derivatives to finish the proof of Theorem 2.

Proposition 3.5.

There exists a positive number ε ~ such that under the assumptions of Proposition 3.4, the solution (ϱ,v,H) of (19) satisfies

| ( ϱ , v , H ) ( t ) | s C ( | ϱ 0 , v 0 , H 0 | L 1 + | ϱ 0 , v 0 , H 0 | s ) × ( 1 + t ) 3 4 ε 1 , 0 < ε 1 < ε ~
(145)

for all t≥0.

Set

σ s = 0 t ( 1 + τ ) 2 σ ( | x ( v , H ) | s + 1 + | x ϱ | s ) + j = 1 [ s + 1 2 ] 0 t ( 1 + τ ) 2 σ | τ j ( ϱ , v , H ) | s + 1 2 j 1 2 ,
(146)

and

N σ s = j = 0 [ s 1 2 ] sup 0 τ t ( 1 + τ ) 2 σ | τ j S 1 ( τ ) | s 1 2 j 2 + j = 0 [ s 2 ] 0 t ( 1 + τ ) 2 σ | τ j ( ϱ div v ) ( τ ) | s 2 j 2 + j = 0 [ s 1 2 ] 0 t ( 1 + τ ) 2 σ | τ j ( v · ϱ ) | s 1 2 j 2 + j = 0 [ s 2 2 ] sup 0 τ t ( 1 + τ ) 2 σ | τ j ( S 2 , S 3 ) ( τ ) | s 2 2 j 2 + j = 0 [ s 1 2 ] 0 t ( 1 + τ ) 2 σ | τ ( S 2 , S 3 ) | s 1 2 j 2 + 0 k + 2 j s 0 t ( 1 + τ ) 2 σ | k τ j S , k τ j V | + 0 k + 2 j s 1 0 t ( 1 + τ ) 2 σ | k τ j S , k τ j + 1 V | + 0 k + l + 2 j s 1 0 t ( 1 + τ ) 2 σ γμ | k 3 l + 1 τ j × ( v · ϱ ) , k 3 l + 1 τ j ϱ | 1 2 ,
(147)

where S=(S1,S2,S3)T and V=(ϱ,v,H)T. As in the proof of (11.6) in (Kagei and Kobayashi [2005]), we have the following inequality in the similar way:

Lemma 3.13.

Under the assumptions of Proposition 3.4, the following inequality holds uniformly in t≥0 and for σ 1 2 ,

K σ s ( t ) 2 + σ s ( t ) 2 C | ϱ 0 , v 0 , H 0 | s 2 + C 0 s ( t ) 2 + C N σ s ( t ) 2 + 0 t ( 1 + τ ) 2 σ 1 | ( ϱ , v , H ) ( τ ) | L 2 2 .
(148)
Proof.

Here we only estimate the magnetic field. We shall show that,

( 1 + t ) 2 σ | H ( t ) | s + j = 1 [ s 2 ] ( 1 + t ) 2 σ | t j H ( t ) | s 2 j + j = 1 [ s + 1 2 ] 0 t ( 1 + τ ) 2 σ | τ j H ( t ) | s + 1 2 j C | ϱ 0 , v 0 , H 0 | s 2 + 0 s ( t ) 2 + N σ s ( t ) 2 + 0 t ( 1 + τ ) 2 σ 1 | H | L 2 2 .
(149)

Then we estimate these as follows:

i) Taking t j to (19) 3, multiplying ( 1 + t ) 2 σ t j H, and integrating on + 3 , it implies

1 2 d dt ( 1 + t ) 2 σ | t j H ( t ) | 2 dx + ( 1 + t ) 2 σ | x τ j H | 2 dx = σ ( 1 + t ) 2 σ 1 | t j H ( t ) | 2 dx + ( 1 + t ) 2 σ t j S 3 , t j H .
(150)

Integrating this from 0 to t yields

( 1 + t ) 2 σ | t j H ( t ) | L 2 2 + 0 t ( 1 + τ ) 2 σ | x τ j H | L 2 2 C | ϱ 0 , v 0 , H 0 | s 2 + C 0 t σ ( 1 + τ ) 2 σ 1 | τ j H | L 2 2 + C 0 t ( 1 + τ ) 2 σ τ j S 3 , τ j H ,
(151)

for 2js.

i i) Taking t j to (19) 3, multiplying ( 1 + t ) 2 σ t j + 1 H, and integrating on + 3 , it holds

( 1 + t ) 2 σ | t j + 1 H ( t ) | 2 dx + 1 2 d dt ( 1 + t ) 2 σ | x t j H ( t ) | 2 dx = σ ( 1 + t ) 2 σ 1 | x t j H ( t ) | 2 dx + ( 1 + t ) 2 σ t j S 3 , t j + 1 H .
(152)

By integrating this from 0 to t, we get

( 1 + t ) 2 σ | x t j H ( t ) | L 2 2 + 0 t ( 1 + τ ) 2 σ | τ j + 1 H | L 2 2 C | ϱ 0 , v 0 , H 0 | s 2 + C 0 t σ ( 1 + τ ) 2 σ 1 | x τ j H | L 2 2 + C 0 t ( 1 + τ ) 2 σ | τ j S 3 , τ j + 1 H | ,
(153)

for 2j+1≤s.

i i i) Due to Lemma 3.7, we have

( 1 + t ) 2 σ x k + 2 t j H ( t ) L 2 2 C ( 1 + t ) 2 σ × t j H k 2 + t j S 3 k 2 + x t j H L 2 2 , for k + 2 j + 2 s
(154)

and

0 t ( 1 + τ ) 2 σ | x k + 2 τ j H | L 2 2 C 0 t ( 1 + τ ) 2 σ × | τ j H | k 2 + | τ j S 3 | k 2 + | x τ j H | L 2 2 ,
(155)

for k+2j+1≤s.

By Cauchy inequality and σ 1 2 , we have (1+t)2σ−1C(α)+α(1+t)2σ,α>0. This combining with (90)–(91), we get

k = 0 s 2 ( 1 + t ) 2 σ | t k H ( t ) | L 2 2 + 0 t ( 1 + τ ) 2 σ | x τ k H | L 2 2 + j = 0 s 1 2 ( 1 + t ) 2 σ | x t j H ( t ) | L 2 2 + 0 t ( 1 + τ ) 2 σ | τ j + 1 H | L 2 2 C | ϱ 0 , v 0 , H 0 | s 2 + 0 t ( 1 + τ ) 2 σ 1 | H | L 2 2 + k = 1 s 2 0 t | τ k H | L 2 2 + j = 0 s 1 2 0 t | x τ j H | L 2 2 + k = 0 s 2 0 t ( 1 + τ ) 2 σ | τ k S 3 , τ k H | + j = 0 s 1 2 0 t ( 1 + τ ) 2 σ | τ j S 3 , τ j + 1 H | .
(156)

Multiplying (92) and (93) for k=0,1 with a sufficiently positive number β, and combining with (94), it yields

2 j + k s , k 3 ( 1 + t ) 2 σ | t j H ( t ) | k 2 + 2 j + k s + 1 , k 3 0 t ( 1 + τ ) 2 σ | τ j H | k 2 C | ϱ 0 , v 0 , H 0 | s 2 + 0 t ( 1 + τ ) 2 σ 1 | H L 2 2 + k = 1 [ s 2 ] 0 t | τ k H | L 2 2 + j = 0 [ s 1 2 ] 0 t | x τ j H L 2 2 + k = 0 [ s 2 ] 0 t ( 1 + τ ) 2 σ | τ k S 3 , τ k H | + j = 0 [ s 1 2 ] 0 t ( 1 + τ ) 2 σ | τ j S 3 , τ j + 1 H | + 2 j + k s , 2 k 3 ( 1 + t ) 2 σ | τ j S 3 | k 2 2 + 2 j + k + 1 s , 2 k 3 0 t ( 1 + τ ) 2 σ | τ j S 3 | k 2 2 .
(157)

For k≥2, the first term in the right-hand side of (92) and (93) can be absorbed by the left-hand side of (92) and (93) for k=k−2. Thus together with (95), we get (89). □

In the proof of Subsubsection A priori estimates, we can obtain that K 0 s ( t ) 2 + 0 s ( t ) 2 C| ϱ 0 , v 0 , H 0 | s 2 . Combining this with Proposition 3.4, we can derive from Lemma 3.13, for σ= 3 4 ε 1 that

K 3 4 ε 1 s ( t ) 2 + 3 4 ε 1 s ( t ) 2 C ( δ 1 ) + C N 3 4 ε 1 s ( t ) 2 ,
(158)

here C( δ 1 )0, as δ 1 0. Hence it remains to estimate N 3 4 ε 1 s (t). We shall show that

N 3 4 ε 1 s ( t ) 2 ε 3 4 ε 1 s ( t ) 2 + C ( ε ) K 3 4 ε 1 s ( t ) 4 + K 3 4 ε 1 s ( t ) 2 0 s ( t ) 2 .
(159)

Thus Proposition 3.5 follows. Since the other terms can be estimated in the similar way or we can see in (Kagei and Kobayashi [2005], here we only estimate 0 t k τ j 1 ρ curl H × H , k τ j v with k=s−2j and s≥3 as

l k τ j 1 ρ curl H × H , k τ j v = s 1 2 j τ j ( 1 ρ curl H × H ) , s + 1 2 j τ j v | s 1 2 j τ j 1 ρ curl H × H | L 2 | | s + 1 2 j τ j v | L 2 C ( ε ) | s 1 2 j τ j ( 1 ρ curl H × H ) | L 2 2 + ε 2 | s + 1 2 j τ j v | L 2 2 .
(160)

In virtue of (17) and its general form, and with the help of the smallness assumption in Proposition 3.4, we have

i) While j=0,

| s 1 1 ρ curl H × H | L 2 2 1 ρ L 2 | H | L 2 | x H | s 1 2 + 1 ρ L 2 | H | s 1 2 | x H | L 2 + x 1 ρ L 2 | H | L 2 | x H | s 2 2 + x 1 ρ L 2 | H | s 2 2 | x H | L 2 + x 1 ρ s 2 2 | H | L 2 | x H | L 2 C | ϱ , H | s 1 2 | x H | s 1 2 .
(161)

i i) While 1j s 1 2 ,

s 1 2 j t j 1 ρ curl H × H L 2 2 1 j 1 + j 2 + + j k j f ρ ; j 1 , j 2 , , j k 2 sgn ( j 1 ) t j 1 ϱ × sgn ( j k 2 ) t j k 2 ϱ × curl t j k 1 H × t j k H s 1 2 j 2 C j = 1 [ s 1 2 ] | t j ( ϱ , H ) | s 1 2 j 2 l = 1 [ s 1 2 ] | x t l H | s 1 2 l 2 ,
(162)

where f= c ρ m , c is a nonzero integer and m=m(j1,j2,,jk−2) is some positive integer. And the last inequality in i i) can be derived from the continuous equation (19) 1 and the smallness assumption of initial data. This together with the definition of 0 s (t), 3 4 ε 1 s (t) and K 3 4 ε 1 s (t), it implies

0 t ( 1 + τ ) 2 3 4 ε 1 t k τ j 1 ρ curl H × H , k τ j v ε 3 4 ε 1 s ( t ) 2 + C ( ε ) K 3 4 ε 1 s ( t ) 2 0 s ( t ) 2 .
(163)

Now we shall finish the proof of the last part of Theorem 2.

Proposition 3.6.

Let s≥4. Under the assumptions of Proposition 3.5, if there exists a positive constants δ2 such that | ϱ 0 , v 0 , H 0 | s +| ϱ 0 , v 0 , H 0 | L 1 δ 2 , then we have

| x ( ϱ , v , H ) | 2 C δ 2 ( 1 + t ) 5 4 ε 2 , 0 < ε 2 < ε ̂
(164)

for all t≥0, where ε ̂ ε ~ is some positive number.

In order to prove Proposition 3.6, we first set

M ~ σ 3 ( t ) = sup 0 τ t ( 1 + τ ) σ x 2 ( ϱ , v , H ) ( τ ) | L 2 + | τ ( ϱ , v , H ) ( τ ) | 1 ,
(165)
K ~ σ s = k = 3 s M σ k + j = 1 [ s 2 ] sup 0 τ t ( 1 + τ ) 3 4 ε | τ j v ( τ ) | s 2 j + j = 1 s + 1 2 sup 0 τ t ( 1 + τ ) 3 4 ε | τ j ϱ ( τ ) | s + 1 2 j 1 2 ,
(166)
~ σ s = 0 t ( 1 + τ ) 2 σ ( | x v | s + | 3 v | s 2 2 + | x 2 ϱ | s 2 + | x τ ( ϱ , v ) | s 2 2 + j = 2 s + 1 2 0 t ( 1 + τ ) 2 σ | τ j ( v , ϱ ) | s + 1 2 j ) 1 2 ,
(167)

and

N ~ σ s = sup 0 τ t ( 1 + τ ) 2 σ | x 2 S 1 ( τ ) | s 3 2 + j = 1 s 1 2 sup 0 τ t ( 1 + τ ) 2 σ | τ j S 1 ( τ ) | s 1 2 j 2 + 0 t ( 1 + τ ) 2 σ | x ( ϱ div v ) ( τ ) | s 1 2 + j = 1 s 2 | τ j ϱ div v ( τ ) | s 2 j 2 + 0 t ( 1 + τ ) 2 σ | x 2 ( v · ϱ ) ( τ ) | s 3 2 + j = 1 s 1 2 | τ j ( v · ϱ ) ( τ ) | s 1 2 j 2 + ( 1 + τ ) 2 σ | x ( S 2 , S 3 ) ( τ ) | s 3 2 + j = 1 s 2 2 sup 0 τ t ( 1 + τ ) 2 σ | τ j ( S 2 , S 3 ) ( τ ) | s 2 2 j 2 + 0 t ( 1 + τ ) 2 σ | x S 2 | s 2 2 + j = 1 s 1 2 | τ j ( S 2 , S 3 ) | s 1 2 j 2 + 0 k + 2 j s , ( k , j ) ( 0 , 0 ) 0 t ( 1 + τ ) 2 σ k τ j S , k τ j V + 0 k + 2 j s 1 , ( k , j ) ( 0 , 0 ) 0 t ( 1 + τ ) 2 σ k τ j S , k τ j + 1 V + 0 k + l + 2 j s 1 , ( k , l , j ) ( 0 , 0 , 0 ) 0 t ( 1 + τ ) 2 σ γμ k 3 l + 1 τ j × ( v · ϱ ) , k 3 l + 1 τ j ϱ 1 2 .
(168)

As in the proof of Lemma 3.13 or we can see in (Kagei and Kobayashi [2005]), the following inequality can be easily deduced:

Lemma 3.14.

Let s≥2. Then under the assumptions of Proposition 3.5, the following inequalities hold uniformly in t≥0,

K ~ 5 4 ε 2 s ( t ) 2 + ~ 5 4 ε 2 s ( t ) 2 C K 0 s ( t ) 2 + 0 s ( t ) 2 + N ~ 5 4 ε 2 s ( t ) 2 + 0 t ( 1 + τ ) 2 5 4 ε 2 1 | x ( ϱ , v , H ) | L 2 2 + | τ ( ϱ , v , H ) | L 2 2 .
(169)

And we estimate N ~ 5 4 ε 2 s (t) as:

Lemma 3.15.

Let s≥2. Then under the assumptions of Proposition 3.5, the following inequalities hold uniformly in t≥0,

N ~ 5 4 ε 2 s ( t ) 2 C K 3 4 ε 1 s ( t ) 4 + K 3 4 ε 1 s ( t ) 2 3 4 ε 1 s ( t ) 2 + C M ( t ) + M 5 4 1 ( t ) + M ~ 5 4 3 ( t ) + K ~ 5 4 ε 2 s ( t ) × 3 4 ε 1 s ( t ) 2 .
(170)
Proof.

Set N ~ 5 4 ε 2 s (t)= J 1 + J 2 ++ J 10 . Here we only consider s≥4 since the case while 2≤s<4 can be deduced more easily. Thus we will estimate J1- J10 term by term.

J 1 , J 2 C sup 0 τ t ( 1 + τ ) 2 5 4 ε 2 | x 3 ( ϱ v ) | s 3 2 + j = 1 s 1 2 | τ j ( ϱ v ) | s 2 j 2 C sup 0 τ t ( 1 + τ ) 2 5 4 ε 2 | ϱ , v | L 2 | ϱ , v | s 2 + j = 1 s 1 2 | ϱ , v | L 2 | τ j ( ϱ , v ) | s 2 j 2 + j = 1 s 1 2 k = 1 j 1 τ k ϱ , τ j k v | L 2 | τ k ϱ , τ j k v s 2 j 2 C sup 0 τ t ( 1 + τ ) 2 5 4 ε 2 | ϱ , v | s 4 + | ϱ , v | s 2 j = 1 s 1 2 | τ j ( ϱ , v ) | s 2 j 2 + k = 1 s 1 2 1 | τ k ϱ | s 2 k 1 2 l = 1 s 1 2 1 | τ l v | s 2 l 1 2 C K 3 4 ε 1 s ( t ) 4 .
(171)

And this holds similarly for J5,J6. For the terms J3, J4 and J7, we only estimate 0 t ( 1 + τ ) 2 ( 5 4 ε 2 ) | x (ϱdivv) | s 1 2 contained in J3. The estimates of the other terms can arrive in the similar way.

0 t ( 1 + τ ) 2 5 4 ε 2 | x ( ϱ div v ) | s 1 2 C 0 t ( 1 + τ ) 2 5 4 ε 2 | ϱ div v | s 2 C 0 t ( 1 + τ ) 2 5 4 ε 2 | ϱ | L 2 | x v | s 2 + | ϱ | s 2 | x v | L 2 C 0 t ( 1 + τ ) 2 5 4 ε 2 | ϱ | s 2 | x v | s 2 C K 3 4 ε 1 s ( t ) 2 3 4 ε 1 s ( t ) 2 .
(172)

For J8- J10, we only estimate the term 1 k s 1 0 t ( 1 + τ ) 2 5 4 ε 2 | x k S 3 || x k τ H| as

1 k s 1 0 t ( 1 + τ ) 2 5 4 ε 2 | x k S 3 | | x k τ H | C 1 k s 1 0 t ( 1 + τ ) 2 5 4 ε 2 | x k + 1 ( v × H ) | | x k τ H | C 1 k s 1 0 t ( 1 + τ ) 2 5 4 ε 2 | v, H | L | x k + 1 ( v,H ) | L 2 + | x ( v , H ) | L 3 | x k ( v,H ) | L 6 + 2 j k 1 | x j ( v,H ) | L 6 | x k + 1 j ( v,H ) | L 3 | x k τ H | L 2 C 1 k s 1 0 t ( 1 + τ ) 2 5 4 ε 2 | v, H | L + | x ( v , H ) | L 2 + | x 2 ( v , H ) | L 2 + 3 j k | x j ( v,H ) | L 2 | x k + 1 ( v,H ) | L 2 | x k τ H | L 2 C M ( t ) + M 5 4 1 ( t ) + M ~ 5 4 3 ( t ) + K ~ 5 4 ε 2 s ( t ) 3 4 ε 1 s ( t ) 2 ,
(173)

which finished the proof of Lemma 3.15. □

By virtue of Lemma 3.14 and Lemma 3.15, it remains to estimate M ~ 5 4 ε 2 3 (t) to finish the proof of Proposition 3.6.

Lemma 3.16.

Let s≥4. Then under the assumptions of Lemma 3.13, the following inequalities hold uniformly in t≥0,

M ~ 5 4 3 ( t ) 2 C K 3 4 ε 1 s ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + M ( t ) 2 .
(174)
Proof.

By Lemma 3.8 i i), we have

| x l t I 0 ( t ) | L 2 C t 5 4 | U 0 | L 1 + | ϱ 0 | 2 + | v 0 | 1 ,
(175)

for t≥1 and l=0,1. Now we shall estimate x l t I 1 (t). Since

t I 1 ( t ) = E ( t τ ) F ( U , H ) ( τ ) | t 1 t + t 1 t E t ( t τ ) F ( U , H ) ( τ ) = E ( t τ ) F ( U , H ) ( τ ) | t 1 t t 1 t E τ ( t τ ) F ( U , H ) ( τ ) = t 1 t E ( t τ ) τ F ( U , H ) ( τ ) ,
(176)

then by Lemma 3.9 and Lemma 3.12, we have

| x l t I 1 ( t ) | L 2 C t 1 t ( t τ ) l 2 | τ F ( U , H ) ( τ ) | H 1 × L 2 C t 1 t ( t τ ) l 2 { | t S 1 | 1 + | t S 2 | L 2 } C ( 1 + t ) 5 4 K 3 4 ε 4 ( t ) 2 + M ( t ) 2 + M 5 4 1 ( t ) 2 .
(177)

Similarly, we can get

t I 2 ( t ) = E ( 1 ) F ( U , H ) ( t 1 ) + 0 t 1 t E ( t τ ) τ F ( U , H ) ( τ ) dτ.
(178)

Thus by Lemma 3.8, Lemma 3.11 and Lemma 3.12,

| x l t I 2 ( t ) | L 2 C | F ( U , H ) ( t 1 ) | 1 + C 0 t 1 ( t τ ) 5 4 | F ( U , H ) ( τ ) | L 1 + e C 0 ( t τ ) | F ( U , H ) ( τ ) | H 1 × L 2 C t 1 t τ 1 2 | t S 1 | 1 + | t S 2 | L 2 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 3 ( t ) 2 + N ( 0 , t ) 2 .
(179)

It then follows that

| t U ( t ) | 1 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + N ( 0 , t ) 2 .
(180)

This together with the magnetic term which can be treated in the similar way, one has

| t ( ϱ , v, H ) | 1 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + N ( 0 , t ) 2 .
(181)

Next we estimate | x 2 U(t) | L 2 and | x 2 v | L 2 . For k=1,2, we see from Lemma 3.8 that

| k x I 0 ( t ) | L 2 C t 7 4 | U 0 | L 1 + | ϱ 0 | k + 1 + | v 0 | k .
(182)

And by Lemma 3.9 and Lemma 3.12,

| k x I 1 ( t ) | L 2 C t 1 t ( t τ ) 1 2 ( | S 1 | k + 1 + | S 2 | k ) C M ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2 × t 1 t ( t τ ) 1 2 ( 1 + τ ) 3 2 C ( 1 + t ) 3 2 M ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2 .
(183)

By Lemma 3.8 and Lemma 3.12,

| k x I 2 ( t ) | L 2 C 0 t 1 ( t τ ) 5 4 | ( S 1 , S 2 ) ( τ ) | L 1 + C 0 t 1 e C 0 ( t τ ) ( | S 1 ( τ ) | k + 1 + | S 2 ( τ ) | k ) C M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2 × 0 t 1 ( t τ ) 5 4 ( 1 + τ ) 7 4 ε 1 + C M ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2 × 0 t 1 e C 0 ( t τ ) ( 1 + τ ) 3 2 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2 .
(184)

Combining these with the magnetic term which has the similar estimates, it yields

| k x ( ϱ , v , H ) | L 2 C ( 1 + t ) 5 4 δ 1 + M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 k + 2 ( t ) 2
(185)

for k=1,2.

Next, we consider 3 2 (ϱ,v). From the equation (19) 2, we have

2 μ 3 2 v = v t + γ ϱ μ ( · ) v μ∂ div v μ 3 ( 1 v 1 + 2 v 2 ) S 2 .
(186)

Thus by (100) and (101), it follows that

| l 3 2 v | L 2 C | l v t | L 2 + | l x ϱ | L 2 + | l + 1 x v | L 2 + | l S 2 | L 2 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + δ 1 .
(187)

Similarly by (19) 1–(19) 2, we have

3 ϱ t + γ 2 2 μ 3 ϱ = γ 2 μ v t 3 γ 2 3 ( 1 v 1 + 2 v 2 ) + γ 2 ( · ) v 3 + 3 S 1 + γ 2 μ S 2 .
(188)

Differentiating this in x3, multiplying 3 2 ϱ to the result equation, and integrating on + 3 , we have

d dt | 3 ϱ | L 2 + γ μ | 3 ϱ | L 2 C ( | x v t , x 2 v , x 2 S 1 , x S 2 | L 2 ) ,
(189)

Thus by Gronwall’s inequality, and with the help of Lemma 3.12, and (87)–(89), we obtain

| 3 2 ϱ | L 2 e γ 4 μ ( t 1 ) | 3 2 ϱ ( 1 ) | L 2 + C 1 t e γ 2 μ ( t τ ) × | x v t , x 2 v , x 2 S 1 , x S 2 | L 2 CN ( 0 , t ) e γ 4 μ ( t 1 ) + C M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + δ 1 × 1 t e γ 2 μ ( t τ ) ( 1 + τ ) 5 4 C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + δ 1 ,
(190)

this together with (101) and (102), we get the estimate of | x 2 (ϱ,v) | L 2 . Hence it remains to estimate | x 2 H | L 2 . By Lemma 2.4 and (100), we have

| x 2 H | L 2 C ( | H t | L 2 + | S 3 | L 2 + | x H | L 2 ) C ( 1 + t ) 5 4 M ( t ) 2 + M 3 4 0 ( t ) 2 + M 5 4 1 ( t ) 2 + K 3 4 ε 1 4 ( t ) 2 + η 1 2 .
(191)

Combining this with the above estimates, Lemma 3.16 follows. □

As in (103), we can deduce the decay rates of higher order derivatives for the magnetic field by Lemma 3.8, Lemma 3.16 and (100). Thus we state the result which is better than one in Proposition 3.6 in the following:

Proposition 3.7.

Under the assumptions of Proposition 3.6, the solution H of (19) 3 satisfies that for all t≥0,

| x 3 H ( t ) | L 2 C ( 1 + t ) 5 4 .
(192)

In view of the above established, the proof of Theorem 2 is complete.

References

  • Chen Q, Tan Z: Global existence and convergence rates of smooth solutions for the compressible magnetohydrodynamic equations. Nonlinear Anal Theory, Methods Appl 2010, 72: 4438–4451. 10.1016/j.na.2010.02.019

    Article  Google Scholar 

  • Chen Q, Tan Z: Global existence in critical spaces for the compressible magnetohydrodynamic equations. Kinet Relat Models 2012, 5: 743–787. 10.3934/krm.2012.5.743

    Article  Google Scholar 

  • Chen GQ, Wang DH: Global solutions of nonlinear magneto-hydrodynamics with large initial data. J Differential Equations 2002, 182: 344–376. 10.1006/jdeq.2001.4111

    Article  Google Scholar 

  • Chen GQ, Wang DH: Existence and continuous dependence of large solutions for the magneto-hydrodynamics equations. Z Angew Math Phys 2003, 54: 608–632. 10.1007/s00033-003-1017-z

    Article  Google Scholar 

  • Cho Y, Choe HJ, Kim H: Unique solvability of the initial boundary value problems for compressible viscous fluids. J Math Pures Appl 2004, 83: 243–275. 10.1016/j.matpur.2003.11.004

    Article  Google Scholar 

  • Ducomet B, Feireisl E: The equations of magneto-hydrodynamics: on the interaction between matter and radiation in the evolution of gaseous stars. Commun Math Phys 2006, 266: 595–629. 10.1007/s00220-006-0052-y

    Article  Google Scholar 

  • Fan JS, Yu W: Strong solution to the compressible magneto-hydrodynamic equations with vacuum. Nonlinear Anal Real World Appl 2009, 10: 392–409. 10.1016/j.nonrwa.2007.10.001

    Article  Google Scholar 

  • Galdi GP, Bris CL, Lelievre T: An introduction to the mathematical theory of the Navier-Stokes equations. Springer-Verlag, New York; 1994.

    Google Scholar 

  • Gerebeau JF, Bris CL, Lelievre T: Mathematical methods for the magneto-hydrodynamics of liquid metals. Oxford University Press, Oxford; 2006.

    Book  Google Scholar 

  • Kagei Y, Kobayashi T: Asymptotic behavior of solutions of the compressible Navier-Stokes equations on the half space. Arch Ration Mech Anal 2005, 177: 231–330. 10.1007/s00205-005-0365-6

    Article  Google Scholar 

  • Kawashima S, Okada M: Smooth global solutions for the one-dimensional equations in magneto-hydrodynamics. Proc Japan Acad Ser A Math Sci 1982, 58: 384–387. 10.3792/pjaa.58.384

    Article  Google Scholar 

  • Kobayashi T:Some estimates of solutions for the equations of motion of compressible viscous fluid in an exterior domain in. J Differential Equations 2002, 184: 587–619. 10.1006/jdeq.2002.4158

    Article  Google Scholar 

  • Kobayashi T, Shibata Y:Decay estimates of solutions for the equations of motion of compressible viscous and heat-conductive gases in an exterior domain inComm Math Phys 1999, 200: 621–659. 10.1007/s002200050543

    Article  Google Scholar 

  • Matsumura A, Nishida T: The initial value problem for the equations of motion of compressible viscous and heatconductive fluids. Proc Japan Acad Ser A 1979, 55: 337–342. 10.3792/pjaa.55.337

    Article  Google Scholar 

  • Matsumura A, Nishida T: The initial value problem for the equations of motion of viscous and heat-conductive gases. J Math Kyoto Univ 1980, 20: 67–104.

    Google Scholar 

  • Matsumura A, Nishida T: Initial boundary value problems for the equations of motion of compressible viscous and heat-conductive fluids. Commun Math Phys 1983, 89: 445–464. 10.1007/BF01214738

    Article  Google Scholar 

  • Ströhmer S: About compressible viscous fluid flow in a bounded regione. Pacific J Math 1990, 143: 359–375. 10.2140/pjm.1990.143.359

    Article  Google Scholar 

  • Tan Z, Wang YJ: Global existence and large-time behavior of weak solutions to the compressible magneto-hydrodynamic equations with coulomb force. Nonlinear Anal Theory, Methods Appl 2009, 71: 5866–5884. 10.1016/j.na.2009.05.012

    Article  Google Scholar 

  • Vol’pert AI, Hudjaev SI: On the Cauchy problem for composite systems of nonlinear differential equations. Math USSR-Sb 1972, 16: 517–544. 10.1070/SM1972v016n04ABEH001438

    Article  Google Scholar 

Download references

Acknowledgements

Q. Chen’s research is supported in part by National Natural Science Foundation of China-NSAF (Nos. 11226174, 11301439). Z. Tan’s research is supported in part by National Natural Science Foundation of China-NSAF (No. 11271305).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Qing Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Author QC composed this paper, and others have revised it many times for publication. Both authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, Q., Tan, Z. Global existence and convergence rates for the smooth solutions to the compressible magnetohydrodynamic equations in the half space. Environ Syst Res 3, 20 (2014). https://doi.org/10.1186/s40068-014-0020-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40068-014-0020-z

Keywords