SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Regularity of global solution to atmospheric circulation equations with humidity effect

Hong Luo

Author Affiliations

College of Mathematics and Software Science, Sichuan Normal University, Chengdu, Sichuan, 610066, China

Boundary Value Problems 2012, 2012:143  doi:10.1186/1687-2770-2012-143


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2012/1/143


Received:1 June 2012
Accepted:14 November 2012
Published:5 December 2012

© 2012 Luo; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this article, the regularity of the global solutions to atmospheric circulation equations with humidity effect is considered. Firstly, the formula of the global solutions is obtained by using the theory of linear operator semigroups. Secondly, the regularity of the global solutions to atmospheric circulation equations is presented by using mathematical induction and regularity estimates for the linear semigroups.

MSC: 35D35, 35K20, 35Q35.

Keywords:
global solution; regularity; atmospheric circulation equations; humidity effect

1 Introduction

This paper is concerned with the regularity of solutions to the following initial-boundary problem of atmospheric circulation equations involving unknown functions ( u , T , q , p ) at ( x , t ) = ( x 1 , x 2 , t ) Ω × ( 0 , ) ( Ω = ( 0 , 2 π ) × ( 0 , 1 ) is a period of C field ( , + ) × ( 0 , 1 ) ):

(1.1)

(1.2)

(1.3)

(1.4)

where P r > 0 , R > 0 , R ˜ , L e > 0 are constants, u = ( u 1 , u 2 ) , T, q, p denote velocity field, temperature, humidity and pressure respectively, Q, G are known functions, and σ is a constant matrix

σ = ( σ 0 ω ω σ 1 ) .

The problems (1.1)-(1.4) are supplemented with the following Dirichlet boundary condition at x 2 = 0 , 1 and the periodic condition for x 1 :

(1.5)

(1.6)

and initial value conditions

( u , T , q ) = ( u 0 , T 0 , q 0 ) , t = 0 . (1.7)

Partial differential equations (1.1)-(1.7) are presented in atmospheric circulation with humidity effect [1]. Atmospheric circulation is one of the main factors affecting the global climate so it is very necessary to understand and master its mysteries and laws. Atmospheric circulation is an important mechanism to complete the transports and balance of atmospheric heat and moisture and the conversion between various energies. Moreover, it is also the important result of these physical transports, balance and conversion. Thus, it is of necessity to study the characteristics, formation, preservation, change and effects of the atmospheric circulation and master its evolution law, which is not only the essential part of human’s understanding of nature, but also a helpful method of changing and improving the accuracy of weather forecasts, exploring global climate change, and making effective use of climate resources.

The atmosphere and ocean around the earth are rotating geophysical fluids, which are also two important components of the climate system. The phenomena of the atmosphere and ocean are extremely rich in their organization and complexity, and a lot of them cannot be produced by laboratory experiments. The atmosphere or the ocean or the couple atmosphere and ocean can be viewed as initial and boundary value problems [2-5], or an infinite dimensional dynamical system [6-8]. We deduce atmospheric circulation models which are able to show the features of atmospheric circulation and are easy to be studied from the very complex atmospheric circulation model based on the actual background and meteorological data, and we present global solutions of atmospheric circulation equations with the use of the T-weakly continuous operator [1]. In [9], the steady state solutions to atmospheric circulation equations with humidity effect are studied. A sufficient condition of the existence of the steady state solutions to atmospheric circulation equations is obtained, and the regularity of the steady state solutions is verified. In this article, we investigate the regularity of the solutions to atmospheric circulation equations (1.1)-(1.7).

The paper is organized as follows. In Section 2, we present preliminary results. In Section 3, we present the formula of the solution to the atmospheric circulation equations. In Section 4, we obtain the regularity of the solutions to equations (1.1)-(1.7).

X denotes the norm of the space X, and C, C i are variable constants.

2 Preliminaries

We consider the divergence form of the linear elliptic equation

L u = D j ( a i j D i u ) + b i D i u + c u = f , (2.1)

where a i j , b i , c L ( Ω ) , f L 2 ( Ω ) , a i j = a j i , ( a i j ) is uniformly elliptic, i.e., there exist constants 0 < λ 1 λ 2 such that

λ 1 | ξ | 2 a i j ( x ) ξ i ξ j λ 2 | ξ | 2 , ξ R n , x Ω .

The problem (2.1) is supplemented with the following Dirichlet boundary condition:

u | Ω = φ . (2.2)

Lemma 2.1[10] (Theory of linear elliptic equations)

Let Ω R n be a C 2 field, a i j C 0 ( Ω ) , b i , c L ( Ω ) , f L p ( Ω ) , φ W 2 , p ( Ω ) . If u W 2 , p is a solution of Eqs. (2.1), (2.2), then

u W 2 , p C ( u L p + f L p + φ W 2 , p ) ,

where C > 0 depends onn, p, λ, Ω and C 0 , α -norm or L -norm of the coefficient functions.

We consider the Stokes equation

{ μ u + p = f ( x ) , div u = 0 , u | Ω = φ . (2.3)

Lemma 2.2[11,12] (ADN theory of the Stokes equation)

Let f W k , p ( Ω , R n ) , φ W k + 2 , p ( Ω , R n ) , k 0 . If ( u , p ) W 2 , p ( Ω , R n ) × W 1 , p ( Ω ) ( 1 < p < ) is a solution of Eq. (2.3), then the solution ( u , p ) W k + 2 , p ( Ω , R n ) × W k + 1 , p ( Ω ) , and

u W k + 2 , p + p W k + 1 , p C ( f W k , p + ( u , p ) L p + φ W k + 2 , p ) ,

where C > 0 depends onμ, n, k, α, Ω.

Let X be a linear space, X 1 , X 2 be two separable reflexive Banach spaces, and H be a Hilbert space. X 1 , X 2 , and H are completion spaces of X under the respective norm. X 1 , X 2 H are dense embedding. F : X 2 × ( 0 , ) X 1 is a continuous mapping. We consider the abstract equation

{ d u d t = F u , 0 < t < , u ( 0 ) = φ , (2.4)

where φ H , u : [ 0 , + ) H is unknown.

Definition 2.3 Let φ H be a given initial value. u L p ( ( 0 , T ) , X 2 ) L ( ( 0 , T ) , H ) ( 0 < T < ) is called a global solution of Eq. (2.4) if u satisfies

( u ( t ) , v ) H = 0 t F u , v d t + ( φ , v ) H , v X 1 H .

Definition 2.4 Let u n , u 0 L p ( ( 0 , T ) , X 2 ) . u n u 0 is called uniformly weak convergence in L p ( ( 0 , T ) , X 2 ) if { u n } L ( ( 0 , T ) , H ) is bounded, and

{ u n u 0 , in  L p ( ( 0 , T ) , X 2 ) , lim n 0 T | u n u 0 , v H | 2 d t = 0 , v H . (2.5)

Definition 2.5 A mapping F : X 2 × ( 0 , ) X 1 is called T-weakly continuous if for p = ( p 1 , , p m ) , 0 < T < and u n uniformly weakly converges to u 0 , we have

lim n 0 T F u n , v d t = 0 T F u 0 , v d t , v X 1 .

Lemma 2.6[3]

Assume F : X 2 × ( 0 , ) X 1 isT-weakly continuous and satisfies:

(A1) there exists p = ( p 1 , , p m ) , p i > 1 ( 1 i m ), such that

F u , u C 1 u X 2 p + C 2 u H 2 + f ( t ) , u X ,

where C 1 , C 2 are constants, f L 1 ( 0 , T ) ( 0 < T < ), X 2 p = i = 1 m | | i p i , | | i is a seminorm of X 2 , X 2 = i = 1 m | | i ,

(A2) there exists 0 < α < 1 for any 0 < h < 1 and u C 1 ( [ 0 , ) , X ) ,

| t t + h F u , v d t | C h α , v X   and   0 t < T ,

where C > 0 depends only onT, v X 1 , 0 t u X 2 p d t , and sup 0 t T u H .

Then for any φ H , Eq. (2.4) has a global weak solution

u L ( ( 0 , T ) , H ) L p ( ( 0 , T ) , X 2 ) , 0 < T < , p   in  (A1) .

If F : X 2 × ( 0 , ) X 1 is Frechét differentiable, then the regular solution can be presented under some condition.

We introduce a space sequence

X H ˜ X 3 X 1 H 1 H ,

where X, X 1 , H are such as in Lemma 2.6, X 3 is a Banach space, H 1 is a Hilbert space, and H 1 H are compact including. There exist a constant C 1 0 and a nonnegative function α L 1 ( 0 , T ) ( 0 < T < ) such that

(2.6)

(2.7)

Lemma 2.7In addition to the assumptions about the existence of a global solution in Lemma 2.6, if F : X 2 × ( 0 , ) X 1 is Frechét differentiable and satisfies (2.6), (2.7), then Eq. (2.4) has a unique global solution

u W 1 , ( ( 0 , T ) , H ) W 1 , 2 ( ( 0 , T ) , H 1 ) , 0 < T < ,

for all φ X 3 .

Lemma 2.8[13]

LetLbe a generator of a strongly continuous semigroup Φ ( t ) . If u ( , φ ) L 1 ( ( 0 , T ) , Y ) is a weak solution to the equation

{ d u d t = L u + F ( u ) , 0 < t < , u ( 0 ) = φ , (2.8)

and F ( u ( , φ ) ) L 1 ( ( 0 , T ) , X ) , then the solution u ( t , φ ) can be read as

u = Φ ( t ) φ + 0 t Φ ( t τ ) F ( u ) d τ , 0 t T .

Note that we used to assume that the linear operator L in (2.8) is a sectorial operator which generates an analytic semigroup Φ ( t ) . It is known that there exists a constant λ 0 such that L λ I generates the fractional power operators L α and fractional order spaces H α for α R 1 , where L = ( L λ I ) . Without loss of generality, we assume that ℒ generates the fractional power operators L α and fractional order spaces H α as follows:

L α = ( L ) α : H α H , α R 1 ,

where H α = D ( L α ) is the domain of L α . By the semigroup theory of linear operators (Pazy [13]), we know that H β H α is a compact inclusion for any β > α .

Lemma 2.9[13-15] (Imbedding theorem of factional order spaces)

Let Ω R n be a Lipschitz field, L : W m , p L p ( Ω ) be a sectorial operator, m 2 , and 1 p < . Then for 0 α 1 , the fractional order spaces H α = D ( L α ) satisfy the following relations:

and the inequalities

For sectorial operators, we also have the following properties which can be found in [13].

Lemma 2.10Let L : H 1 H be a sectorial operator which generates an analytic semigroup T ( t ) = e t L . If all eigenvaluesλofLsatisfy Re λ < λ 0 for some real number λ 0 > 0 , then for L α ( L = L ) , we have

(1) T ( t ) : H H α is bounded for all α R 1 and t > 0 ,

(2) T ( t ) L α x = L α T ( t ) x , x H α ,

(3) for each t > 0 , L α T ( t ) : H H is bounded and

L α T ( t ) C α t α e δ t ,

where some δ > 0 , C α > 0 is a constant only depending onα,

(4) the H α -norm can be defined by

x H α = L α x H , (2.9)

(5) ifis symmetric, for any α , β R 1 , we have

L α u , v H = L α β u , L β v H .

3 Formula of global solutions

We introduce the spaces

Let

Then Eqs. (1.1)-(1.7) can be rewritten as an abstract equation

d ϕ d t = L ϕ + F ( ϕ ) .

Theorem 3.1If ϕ 0 = ( u 0 , T 0 , q 0 ) H , Q , G L 2 ( Ω ) , then the global solutionϕof Eqs. (1.1)-(1.7) can be read as

ϕ ( x , t ) = Φ ( t ) ϕ 0 + 0 t Φ ( t τ ) P ( F ( ϕ ) ) d τ , (3.1)

where Φ ( t ) is an analytic semigroup generated byL, and P : L 2 ( Ω , R 4 ) H is a Leray projection.

Proof As ϕ = ( u , T , q ) L ( ( 0 , T ) , H ) L 2 ( ( 0 , T ) , H 1 ) is a weak solution to Eqs. (1.1)-(1.7) [1], from the Hölder inequality and the Sobolev imbedding theorem, it follows that

0 T [ Ω | ( u ) u | 3 2 d x ] 2 3 d t 0 T [ Ω | u | 3 2 | u | 3 2 d x ] 2 3 d t 0 T [ ( Ω | u | 6 d x ) 1 6 ( Ω | D u | 2 d x ) 1 2 ] d t C 0 T Ω | D u | 2 d x d t C u L 2 ( 0 , T ; H 1 ) 2 .

Then ( u ) u L 1 ( ( 0 , T ) , L 3 2 ( Ω ) ) . Hence,

P r σ u + P r ( R T R ˜ q ) κ ( u ) u L 1 ( ( 0 , T ) , L 3 2 ( Ω ) ) . (3.2)

From the Hölder inequality and the Sobolev imbedding theorem, we see

0 T [ Ω | ( u ) T | 3 2 d x ] 2 3 d t 0 T [ Ω | u | 3 2 | D T | 3 2 d x ] 2 3 d t 0 T [ ( Ω | u | 6 d x ) 1 6 ( Ω | D T | 2 d x ) 1 2 ] d t C 0 T [ ( Ω | D u | 2 d x ) 1 2 ( Ω | D T | 2 d x ) 1 2 ] d t C ( 0 T Ω | D u | 2 d x d t ) 1 2 ( 0 T Ω | D T | 2 d x d t ) 1 2 C T L 2 ( 0 , T ; H 1 ) u L 2 ( 0 , T ; H 1 ) .

Then ( u ) T L 1 ( ( 0 , T ) , L 3 2 ( Ω ) ) . Thus,

u 2 ( u ) T + Q L 1 ( ( 0 , T ) , L 3 2 ( Ω ) ) . (3.3)

Similarly, we have

u 2 ( u ) q + G L 1 ( ( 0 , T ) , L 3 2 ( Ω ) ) . (3.4)

According to the ADN theory and the theory of linear elliptic equations, we have that

L ( ϕ , p ) = ( P r ( Δ u p ) Δ T L e Δ q )

is a sectorial operator and

L : W 2 , q ( Ω , R 4 ) × W 1 , q ( Ω ) L q ( Ω , R 4 ) , 1 < q < .

Therefore, L generates the analytic semigroup Φ ( t ) .

It follows from (3.2), (3.3), and (3.4) that

F ( ϕ ) L 1 ( ( 0 , T ) , L 3 2 ( Ω , R 4 ) ) .

Applying Lemma 2.8 yields

ϕ ( x , t ) = Φ ( t ) ϕ 0 + 0 t Φ ( t τ ) P ( F ( ϕ ) ) d τ .

 □

Remark 3.2 The analytic semigroup Φ ( t ) generated by L can be read as

Φ ( t ) = ( Φ 1 ( t ) Φ 2 ( t ) Φ 3 ( t ) ) .

Remark 3.3 The semigroup generated by Eqs. (1.1)-(1.7) can be rewritten as

( u T q ) = ( Φ 1 ( t ) ( ϕ 0 ) + 0 t Φ 1 ( t τ ) P ( F 1 ) d τ Φ 2 ( t ) ( ϕ 0 ) + 0 t Φ 2 ( t τ ) P ( F 2 ) d τ Φ 3 ( t ) ( ϕ 0 ) + 0 t Φ 3 ( t τ ) P ( F 3 ) d τ ) .

4 Regularity of global solution

Theorem 4.1If ϕ 0 = ( u 0 , T 0 , q 0 ) H 2 ( Ω , R 4 ) H 1 , Q , G L 2 ( Ω ) , then Eqs. (1.1)-(1.7) have a unique solution ( ϕ , p ) = ( u , T , q , p ) , and

for all 0 < T < .

Proof Let X 2 = X 1 = H 1 and X 3 = H 2 ( Ω , R 4 ) H 1 . Define F : H 1 H 1 as

F ϕ , ψ = Ω [ P r u v P r σ u u + P r ( R T R ˜ q ) v 2 ( u ) u v T S + u 2 S ( u ) T S + Q S L e q z + u 2 z ( u ) q z + G z ] d x , for   ψ = ( v , S , z ) H 1 .

Then

F ϕ , ψ = Ω [ P r Δ u P r σ u + P r ( R T R ˜ q ) κ ˜ ( u ) u ] v + [ Δ T + u 2 ( u ) T + Q ] S + [ L e Δ q + u 2 ( u ) q + G ] z d x 1 2 Ω ( | v | 2 + | S | 2 + | z | 2 ) d x + C Ω ( | Δ u | 2 + | T | 2 + | q | 2 + | u D u | 2 + | Δ T | 2 + | u 2 | 2 + | u D T | 2 + | Q | 2 + | Δ q | 2 + | u D q | 2 + | G | 2 ) d x 1 2 ψ H 2 + g ( ϕ )

for any ψ H 1 and ϕ X 3 . Then (2.7) holds.

We prove (2.6).

By the interpolation inequality [16], we see

v X 1 4 = v H 1 2 C v L 2 1 2 v H 1 1 2 . (4.1)

By the imbedding theorem of factional order spaces, we have

( Ω | v | 4 d x ) 1 4 C v X 1 4 . (4.2)

Then it follows from (4.1) and 4.2) that

Since ϕ = ( u , T , q ) L 2 ( ( 0 , T ) , H 1 ) is a weak solution to Eqs. (1.1)-(1.7), α ( t ) = C ( u H 1 2 + T H 1 2 + q H 1 2 + 1 ) L 1 ( 0 , T ) , 0 < T < . Then (2.6) and (2.7) hold. From Lemma 2.7, we deduce that the solution ϕ is unique and

ϕ t = ( u t , T t , q t ) L ( ( 0 , T ) , H ) L 2 ( ( 0 , T ) , H 1 ) . (4.3)

Multiplying (1.1) by u and integrating over Ω, we get

P r Ω | u | 2 d x + P r Ω σ u u d x = Ω [ P r ( R T R ˜ q ) κ u u t u ] d x .

Using the Young inequality, we obtain

P r Ω | u | 2 d x + P r σ ˜ Ω u 2 d x ε Ω | u | 2 d x + C ε 1 Ω [ | T | 2 + | q | 2 + | u t | 2 ] d x ,

where ε > 0 is a real constant satisfying P r σ ˜ ε . Then there exists a constant C > 0 such that

Ω | u | 2 d x C ( T L 2 2 + q L 2 2 + u t L 2 2 ) .

Thanks to (4.3), we have

u L ( ( 0 , T ) , H 1 ( Ω , R 2 ) ) . (4.4)

We consider the Stokes equation

{ Δ u + p = g ( x , t ) , g = σ u + ( R T R ˜ q ) κ 1 P r ( u ) u u t , div u = 0 , u = 0 , x 2 = 0 , 1 , u ( 0 , x 2 ) = u ( 2 π , x 2 ) . (4.5)

From (4.3), (4.4), and the Sobolev imbedding theorem, we find that g ( x , t ) L q ( Ω , R 2 ) , 1 < q < 2 . By the ADN theorem, Eq. (4.5) has a solution

u L ( ( 0 , T ) , W 2 , q ( Ω , R 2 ) ) .

Then ( u ) u L 2 ( Ω , R 2 ) and g ( x , t ) L 2 ( Ω , R 2 ) . Using the ADN theorem, we obtain

{ u L ( ( 0 , T ) , H 2 ( Ω , R 2 ) ) , p L ( ( 0 , T ) , H 1 ( Ω ) ) . (4.6)

Multiplying (1.2) by T and integrating over Ω, we get

Ω | T | 2 d x = Ω ( u 2 T + Q T T t T ) d x ε Ω | T | 2 d x + C ε 1 Ω [ | u | 2 + | Q | 2 + | T t | 2 ] d x ,

where ε > 0 is a constant. Then there exists a constant C > 0 such that

Ω | T | 2 d x C ( u L 2 2 + Q L 2 2 + T t L 2 2 ) .

Using (4.3), we have

T L ( ( 0 , T ) , H 1 ) . (4.7)

We consider the elliptic equation

{ Δ T = f 1 , f 1 = u 2 ( u ) T + Q T t , T = 0 , x 2 = 0 , 1 , T ( 0 , x 2 ) = T ( 2 π , x 2 ) . (4.8)

It follows from (4.3), (4.7), and the Sobolev imbedding theorem that f 1 ( x , t ) L q ( Ω ) , 1 < q < 2 . Using the theory of linear elliptic equations, Eq. (4.8) has a solution

T L ( ( 0 , T ) , W 2 , q ( Ω ) ) .

Then ( u ) T L 2 ( Ω , R 2 ) and f 1 ( x , t ) L 2 ( Ω , R 2 ) . Using the theory of linear elliptic equations, we have that

T L ( ( 0 , T ) , H 2 ( Ω ) ) . (4.9)

Multiplying (1.3) by q and integrating over Ω, we get

Ω | q | 2 d x = Ω ( u 2 q + G q q t ) d x ε Ω | q | 2 d x + C ε 1 Ω [ | u | 2 + | G | 2 + | q t | 2 ] d x ,

where ε > 0 is a constant. Then there exists a constant C > 0 such that

Ω | q | 2 d x C ( u L 2 2 + Q L 2 2 + q t L 2 2 ) .

Using (4.3), we have

q L ( ( 0 , T ) , H 1 ) . (4.10)

We consider the elliptic equation

{ Δ q = f 2 , f 2 = u 2 ( u ) q + G q t , q = 0 , x 2 = 0 , 1 , q ( 2 π , x 2 ) = q ( 0 , x 2 ) . (4.11)

Using the arguments similar to those for (4.8), we get

q L ( ( 0 , T ) , H 2 ( Ω ) ) . (4.12)

It follows from (4.6), (4.9), and (4.12) that

 □

Theorem 4.2If ϕ 0 = ( u 0 , T 0 , q 0 ) H k + 2 ( Ω , R 4 ) H 1 , Q , G H k ( Ω ) , then Eqs. (1.1)-(1.7) have a higher regular solution ( ϕ , p ) = ( u , T , q , p ) and

for 0 < T < .

Proof We prove the theorem using mathematical induction.

If k = 1 , ϕ 0 = ( u 0 , T 0 , q 0 ) H 3 ( Ω , R 4 ) H 1 , Q , G H 1 ( Ω ) , then ϕ 0 = ( u 0 , T 0 , q 0 ) H 2 ( Ω , R 4 ) H 1 , Q , G L 2 ( Ω ) . Using Theorem 4.1, we find that ϕ = ( u , T , q ) H 2 ( Ω , R 4 ) .

Thanks to the Sobolev imbedding theorem, H 2 ( Ω ) W 1 , 4 ( Ω ) C 0 ( Ω ) if n = 2 . We obtain

Then ( u ) u H 1 ( Ω ) and F 1 ( ϕ ) H 1 ( Ω ) .

We have from the formula (3.1)

u t = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

u t H 1 α = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) H 1 α L 1 Φ 1 ( t ) u 0 H 1 α + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ H 1 α + F 1 ( Φ ) H 1 α C + 0 t L 1 3 α 2 Φ 1 ( t τ ) F 1 ( Φ ) L 2 d τ C + 0 t L 1 1 α 2 Φ 1 ( t τ ) L 1 1 2 F 1 ( Φ ) L 2 d τ C .

Thus, u t H 1 α ( Ω ) . Then g H 1 α ( Ω ) in Eq. (4.5). By the ADN theory, u H 3 α ( Ω ) . Thus, ( u ) u H 2 α ( Ω ) and F 1 ( ϕ ) H 2 α ( Ω ) .

We have

u t H 1 = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) H 1 L 1 Φ 1 ( t ) u 0 H 1 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ H 1 + F 1 ( Φ ) H 1 C + 0 t L 1 3 2 Φ 1 ( t τ ) F 1 ( Φ ) L 2 d τ C + 0 t L 1 1 + α 2 Φ 1 ( t τ ) L 1 1 α 2 F 1 ( Φ ) L 2 d τ C ,

which implies u t H 1 ( Ω ) . Then g H 1 ( Ω ) in Eq. (4.5). Using the ADN theory, u H 3 ( Ω ) and p H 2 ( Ω ) . Thus, u L ( ( 0 , T ) , H 3 ( Ω , R 2 ) ) and p L ( ( 0 , T ) , H 2 ( Ω ) ) . Then u C ( [ 0 , T ] , H 3 ( Ω , R 2 ) H 1 ) and u t C ( ( 0 , T ] , H 1 ( Ω , R 2 ) ) from the formula (3.1).

Similarly,

Ω | ( u ) T | 2 d x C Ω | u | 4 + | T | 2 | u | 2 d x C ( u W 1 , 4 4 + sup Ω | u | 2 T H 2 2 ) C ( u H 2 4 + u H 2 2 T H 2 2 ) .

Then ( u ) T H 1 ( Ω ) and F 2 ( ϕ , p ) H 1 ( Ω ) .

We have from the formula (3.1)

T t = L 2 Φ 2 ( t ) u 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

T t H 1 α = L 2 Φ 2 ( t ) T 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) H 1 α L 2 Φ 2 ( t ) T 0 H 1 α + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ H 1 α + F 2 ( Φ ) H 1 α C + 0 t L 2 3 α 2 Φ 2 ( t τ ) F 2 ( Φ ) L 2 d τ C + 0 t L 2 1 α 2 Φ 2 ( t τ ) L 2 1 2 F 2 ( Φ ) L 2 d τ C ,

which implies T t H 1 α ( Ω ) . Then f 1 H 1 α ( Ω ) in Eq. (4.8). It follows from the linear elliptic equation T H 3 α ( Ω ) . Thus, ( u ) T H 2 α ( Ω ) and F 2 ( ϕ ) H 2 α ( Ω ) .

Then

T t H 1 = L 2 Φ 2 ( t ) T 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) H 1 L 2 Φ 2 ( t ) T 0 H 1 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ H 1 + F 2 ( Φ ) H 1 C + 0 t L 2 3 2 Φ 2 ( t τ ) F 2 ( Φ ) L 2 d τ C + 0 t L 2 1 + α 2 Φ 2 ( t τ ) L 2 1 α 2 F 2 ( Φ ) L 2 d τ C ,

which implies T t H 1 ( Ω ) . We obtain that f 1 H 1 ( Ω ) in Eq. (4.8). Then T H 3 ( Ω ) from the theory of linear elliptic equations. Thus, T L ( ( 0 , T ) , H 3 ( Ω , R 2 ) ) . From the formula (3.1), T C ( [ 0 , T ] , H 3 ( Ω ) ) and T t C ( ( 0 , T ] , H 1 ( Ω ) ) .

Similarly,

Ω | ( u ) q | 2 d x C Ω | u | 4 + | q | 2 | u | 2 d x C ( u W 1 , 4 4 + sup Ω | u | 2 q H 2 2 ) C ( u H 2 4 + u H 2 2 q H 2 2 ) .

Then ( u ) q H 1 ( Ω ) and F 3 ( ϕ , p ) H 1 ( Ω ) .

We have from the formula (3.1)

q t = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

q t H 1 α = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) H 1 α L 3 Φ 3 ( t ) u 0 H 1 α + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ H 1 α + F 3 ( Φ ) H 1 α C + 0 t L 3 3 α 2 Φ 3 ( t τ ) F 3 ( Φ ) L 2 d τ C + 0 t L 3 1 α 2 Φ 3 ( t τ ) L 3 1 2 F 3 ( Φ ) L 2 d τ C ,

which implies q t H 1 α ( Ω ) . Then f 2 H 1 α ( Ω ) in Eq. (4.11). Thus, q H 3 α ( Ω ) from the theory of linear elliptic equations. Then ( u ) q H 2 α ( Ω ) and F 3 ( ϕ ) H 2 α ( Ω ) .

Thus,

q t H 1 = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) H 1 L 3 Φ 3 ( t ) q 0 H 1 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ H 1 + F 3 ( Φ ) H 1 C + 0 t L 3 3 2 Φ 3 ( t τ ) F 3 ( Φ ) L 2 d τ C + 0 t L 3 1 + α 2 Φ 3 ( t τ ) L 3 1 α 2 F 3 ( Φ ) L 2 d τ C ,

which implies q t H 1 ( Ω ) . We see f 2 H 1 ( Ω ) in Eq. (4.11). Then q H 3 ( Ω ) from the theory of linear elliptic equations. Thus, q L ( ( 0 , T ) , H 3 ( Ω ) ) . We have q C ( [ 0 , T ] , H 3 ( Ω ) ) and q t C ( ( 0 , T ] , H 1 ( Ω ) ) from the formula (3.1).

It follows from Eq. (4.5) that

p = Δ u σ u + ( R T R ˜ q ) κ 1 P r ( u ) u u t .

Clearly, the right-hand side of the above equality is continuous in H 1 ( Ω ) . Thus,

p C ( ( 0 , T ) , H 2 ( Ω ) ) .

If k = m , ϕ 0 = ( u 0 , T 0 , q 0 ) H m + 2 ( Ω , R 4 ) H 1 , and Q , G H m ( Ω ) , then ϕ 0 = ( u 0 , T 0 , q 0 ) H m + 1 ( Ω , R 4 ) H 1 and Q , G H m 1 ( Ω ) . From the hypothesis of mathematical induction, we see ϕ = ( u , T , q ) H m + 1 ( Ω , R 4 ) .

By the Sobolev imbedding theorem, we have H m + 1 ( Ω ) W m , 4 ( Ω ) C m 1 ( Ω ) if n = 2 . Then it follows from the Sobolev imbedding theorem and the interpolation inequality that

Ω | m ( u ) u | 2 d x C Ω [ i = 0 m | i u | 2 | m + 1 i u | 2 ] d x C Ω [ i = 1 m | i u | 2 | m + 1 i u | 2 + | u | 2 | m + 1 u | 2 ] d x C [ i = 1 m ( Ω | i u | 4 d x ) 1 2 ( Ω | m + 1 i u | 4 d x ) 1 2 + sup Ω | u | 2 Ω | m + 1 u | 2 d x ] C ( u W 1 , 4 2 u W m , 4 2 + sup Ω | u | 2 u H m + 1 2 ) , C u H m + 1 4 .

Then ( u ) u H m ( Ω ) and F 1 ( ϕ ) H m ( Ω ) .

We have from the formula (3.1)

u t = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

u t H m α = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) H m α L 1 Φ 1 ( t ) u 0 H m α + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ H m α + F 1 ( Φ ) H m α C + 0 t L 1 2 + m α 2 Φ 1 ( t τ ) F 1 ( Φ ) L 2 d τ C + 0 t L 1 1 α 2 Φ 1 ( t τ ) L 1 m 2 F 1 ( Φ ) L 2 d τ C .

Then u t H m α ( Ω ) . We see that g H m α ( Ω ) in Eq. (4.5). Thus, u H m + 2 α ( Ω ) from the ADN theory. Hence, ( u ) u H m + 1 α ( Ω ) and F 1 ( ϕ ) H m + 1 α ( Ω ) . Then

u t H m = L 1 Φ 1 ( t ) u 0 + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ + F 1 ( Φ ) H m L 1 Φ 1 ( t ) u 0 H m + 0 t L 1 Φ 1 ( t τ ) F 1 ( Φ ) d τ H m + F 1 ( Φ ) H m C + 0 t L 1 1 + m 2 Φ 1 ( t τ ) F 1 ( Φ ) L 2 d τ C + 0 t L 1 1 + α 2 Φ 1 ( t τ ) L 1 m + 1 α 2 F 1 ( Φ ) L 2 d τ C ,

which implies u t H m ( Ω ) . Then g H m ( Ω ) in Eq. (4.5). Using the ADN theory, u H m + 2 ( Ω ) , and p H m + 1 ( Ω ) , we get

u L ( ( 0 , T ) , H m + 2 ( Ω , R 2 ) ) , p L ( ( 0 , T ) , H m + 1 ( Ω ) ) .

From the formula (3.1), we have

u C ( [ 0 , T ] , H m + 2 ( Ω , R 2 ) ) , u t C ( ( 0 , T ] , H m ( Ω , R 2 ) ) .

Similarly,

Then ( u ) T H m ( Ω ) and F 2 ( ϕ ) H m ( Ω ) .

We have from the formula (3.1)

T t = L 2 Φ 2 ( t ) T 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

T t H m α = L 2 Φ 2 ( t ) T 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) H m α L 2 Φ 2 ( t ) T 0 H m α + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ H m α + F 2 ( Φ ) H m α C + 0 t L 2 2 + m α 2 Φ 2 ( t τ ) F 2 ( Φ ) L 2 d τ C + 0 t L 2 1 α 2 Φ 2 ( t τ ) L 2 m 2 F 2 ( Φ ) L 2 d τ C ,

which implies T t H m α ( Ω ) . Then f 1 H m α ( Ω ) in Eq. (4.8). It follows from the linear elliptic equation T H m + 2 α ( Ω ) that ( u ) T H m + 1 α ( Ω ) and F 2 ( ϕ ) H m + 1 α ( Ω ) . We obtain

T t H m = L 2 Φ 2 ( t ) T 0 + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ + F 2 ( Φ ) H m L 2 Φ 2 ( t ) u 0 H m + 0 t L 2 Φ 2 ( t τ ) F 2 ( Φ ) d τ H m + F 2 ( Φ ) H m C + 0 t L 2 1 + m 2 Φ 2 ( t τ ) F 2 ( Φ ) L 2 d τ C + 0 t L 2 1 + α 2 Φ 2 ( t τ ) L 2 m + 1 α 2 F 2 ( Φ ) L 2 d τ C .

Then T t H m ( Ω ) . We have f 1 H m ( Ω ) in Eq. (4.8). Then T H m + 2 ( Ω ) from the theory of linear elliptic equations. Thus,

T L ( ( 0 , T ) , H m + 2 ( Ω ) ) .

From the formula (3.1), we induce

T C ( [ 0 , T ] , H m + 2 ( Ω ) ) , u t C ( ( 0 , T ] , H m ( Ω ) ) .

Similarly,

Then ( u ) q H m ( Ω ) and F 3 ( ϕ ) H m ( Ω ) .

We have from the formula (3.1)

q t = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) .

Then there exists α satisfying 0 < α < 1 such that

q t H m α = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) H m α L 3 Φ 3 ( t ) u 0 H m α + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ H m α + F 3 ( Φ ) H m α C + 0 t L 3 2 + m α 2 Φ 3 ( t τ ) F 3 ( Φ ) L 2 d τ C + 0 t L 3 1 α 2 Φ 3 ( t τ ) L 3 m 2 F 3 ( Φ ) L 2 d τ C ,

which implies q t H m α ( Ω ) . Then f 2 H m α ( Ω ) in Eq. (4.11). Thus, q H m + 2 α ( Ω ) from the theory of linear elliptic equations. Then ( u ) q H m + 1 α ( Ω ) and F 3 ( ϕ ) H m + 1 α ( Ω ) . Thus,

q t H m = L 3 Φ 3 ( t ) q 0 + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ + F 3 ( Φ ) H m L 3 Φ 3 ( t ) q 0 H m + 0 t L 3 Φ 3 ( t τ ) F 3 ( Φ ) d τ H m + F 3 ( Φ ) H m C + 0 t L 3 1 + m 2 Φ 3 ( t τ ) F 3 ( Φ ) L 2 d τ C + 0 t L 3 1 + α 2 Φ 3 ( t τ ) L 3 m + 1 α 2 F 3 ( Φ ) L 2 d τ C ,

which implies q t H m ( Ω ) . We find f 2 H m ( Ω ) in Eq. (4.11). Then q H m + 2 ( Ω ) from the theory of linear elliptic equations. We have

q L ( ( 0 , T ) , H m + 2 ( Ω ) ) .

From the formula (3.1), we see

q C ( [ 0 , T ] , H m + 2 ( Ω ) ) , q t C ( ( 0 , T ] , H m ( Ω ) ) .

It follows from Eq. (4.5) that

p = Δ u σ u + ( R T R ˜ q ) κ 1 P r ( u ) u u t .

Clearly, the right-hand side of the above equality is continuous in H m ( Ω ) . Then

p C ( ( 0 , T ) , H m + 1 ( Ω ) ) .

The proof is completed. □

Since the differentiability of time and of space can be transformed into each other, we obtain

Remark 4.3 If ϕ 0 = ( u 0 , T 0 , q 0 ) H k + 2 ( Ω , R 4 ) H 1 , Q , G H k ( Ω ) , then Eqs. (1.1)-(1.7) have a higher regular solution ( ϕ , p ) = ( u , T , q , p ) , and

for 0 < T < , where l, r, α, β are positive integers satisfying 2 l + r = k + 2 and 2 α + β = k + 1 .

Competing interests

The author declares that he has no competing interests.

Acknowledgements

The author is very grateful to the anonymous referees whose careful reading of the manuscript and valuable comments enhanced presentation of the manuscript. The project is supported by the National Natural Science Foundation of China (11271271), the NSF of Sichuan Science and Technology Department of China (2010JY0057) and the NSF of Sichuan Education Department of China (11ZA102).

References

  1. Luo, H: Global solution of atmospheric circulation equations with humidity effect (submitted)

  2. Ma, T, Wang, SH: Phase Transition Dynamics in Nonlinear Sciences, Springer, New York (2012)

  3. Ma, T: Theories and Methods in Partial Differential Equations, Science Press, China (2011) in Chinese

  4. Phillips, NA: The general circulation of the atmosphere: a numerical experiment. Q. J. R. Meteorol. Soc.. 82, 123–164 (1956). Publisher Full Text OpenURL

  5. Rossby, CG: On the solution of problems of atmospheric motion by means of model experiment. Mon. Weather Rev.. 54, 237–240 (1926). Publisher Full Text OpenURL

  6. Lions, JL, Temam, R, Wang, SH: New formulations of the primitive equations of atmosphere and applications. Nonlinearity. 5(2), 237–288 (1992). Publisher Full Text OpenURL

  7. Lions, JL, Temam, R, Wang, SH: On the equations of the large-scale ocean. Nonlinearity. 5(5), 1007–1053 (1992). Publisher Full Text OpenURL

  8. Lions, JL, Temam, R, Wang, SH: Models for the coupled atmosphere and ocean. (CAO I,II). Comput. Mech. Adv.. 1(1), 5–54 (1993)

  9. Luo, H: Steady state solution to atmospheric circulation equations with humidity effect. J. Appl. Math. doi:10.1155/2012/867310 (2012)

  10. Evens, LC: Partial Differential Equations, Am. Math. Soc., Providence (1998)

  11. Temam, R: Navier-Stokes Equation and Nonlinear Functional Analysis, SIAM, Philadelphia (1983)

  12. Temam, R: Navier-Stokes Equations: Theory and Numerical Analysis, North-Holland, Amsterdam (1979)

  13. Pazy, A: Semigroups of Linear Operators and Applications to Partial Differential Equations, Springer, Berlin (2006)

  14. Ma, T, Wang, SH: Stability and Bifurcation of Nonlinear Evolution Equations, Science Press, China (2007) in Chinese

  15. Ma, T, Wang, SH: Bifurcation Theory and Applications, World Scientific, Singapore (2005)

  16. Temam, R: Infinite-Dimensional Dynamical Systems in Mechanics and Physics, Springer, New York (1997)