1 Introduction
We consider the primitive equations
|
|
|
(1.5) |
where is a vector field, is a scalar function, is the horizontal gradient, and is the flat torus.
We impose the periodic boundary conditions.
The second equation in (1.5) is called the hydrostatic approximation.
The vertical component is give by the formula
|
|
|
(1.6) |
where is the horizontal divergence.
This formula is from the divergence-free condition.
We also impose .
We always assume this assumption to the horizontal component of three-dimensional divergence-free vector fields in this paper.
In the case of the Neumann boundary conditions in the domain are reduced to the periodic case in by the even and odd extension for and , respectively.
We invoke that satisfies the nonlinear parabolic equation
|
|
|
(1.7) |
for , where and .
Note that (1.7) differs from the equation for derived in Proposition 4.6 of [3] in that no vertical derivatives appear in the right-hand-side terms.
The derivation is described in Appendix A.
Existence of the global weak solution to the primitive equations on spherical shells was proved by Lions, Temam and Wang [13].
Local well-posedness was proved by GuillΓ©n-GonzΓ‘lez, Masmoudi and RodrΓguez-Bellido [10] for initial data, where is the Sobolev space with .
Cao and Titi [2] proved a energy bound to establish the global well-posedness.
Hieber and Kashiwabara [11] extended this result and proved the global well-posedness in Lebesgue spaces -settings for .
Recently, Giga, Gries, Hieber, Hussein, and Kashiwabara [5] showed the global well-posedness in - settings under the periodic, Neumann, Dirichlet, Dirichlet-Neumann mixed boundary conditions.
Giga, Gries, Hieber, Hussein, and Kashiwabara [8] showed the global well-posedness in , where for denotes an anisotropic space equipped with the norm
|
|
|
(1.8) |
Giga, Gries, Hieber, Hussein, and Kashiwabara [7] also proved the global well-posedness in for and under the Dirichlet-Neumann mixed boundary conditions.
An advantage of -approach is that one need not assume smoothness for initial data.
The aim of this paper is to give a mathematically rigorous justification of the hydrostatic approximation for the primitive equations under less smoothness assumptions than the previous works.
We first introduce a brief derivation of the primitive equations.
The primitive equation is derived by the Navier-Stokes equations with anisotropic viscosity, which are horizontally and vertically .
Applying a scaling to equations, we obtain the scaled Navier-Stokes equations
|
|
|
(1.13) |
see [3], [4], and [12] for the details.
If we multiply to the seconde equation of (1.13) and take formal limit , then we obtain (1.5).
To justify this formal derivation we have to show the difference between the solutions to (1.5) and (1.13) converges to zero in some topologies.
We put
|
|
|
Then we see that satisfies
|
|
|
(1.18) |
where
|
|
|
(1.19) |
Note that
|
|
|
where for a vector field .
The justification of the hydrostatic approximation is reduced to showing that converges to zero.
AzΓ©rad and GuillΓ©n [1] proved the weak convergence in the energy space.
Li and Titi [12] showed the strong convergence in the energy space.
They also proved global well-posedness to (1.13) for small compared to the initial data.
The authors together with Giga, Hieber, Hussein, and Wrona [3] extended Li and Titiβs result to - settings under the Neumann boundary conditions.
The authors together with Giga [4] showed the strong convergence in - settings under the Dirichlet boundary conditions.
We consider the solution to (1.18) in the sense of a mild solution, namely
|
|
|
(1.24) |
where is the anisotropic Helmholtz projection which maps from -vector fields to -free -vector fields.
The first main result of this paper is the global well-posedness to (1.18) in setting for small compared to the initial data of the primitive equations.
We write equipped with the norm (1.8) and equipped with the norm
|
|
|
(1.25) |
for and an interval .
Theorem 1.1.
Let , , satisfy and , and .
Let be a solution to (1.5) with initial data .
Then there exists such that, if the equation of the differences (1.18) admits a unique solution
|
|
|
|
|
|
|
|
|
|
|
|
(1.26) |
where is independent of .
The global well-posedness to the primitive equations in , Theorem 1.1, and the definition of imply
Corollary 1.3.
Under the same assumptions for , if , then there exists a unique solution to (1.13) such that
|
|
|
|
|
|
|
|
|
(1.27) |
where are constants independent of .
The proof of Theorem 1.1 based on the contraction mapping principle.
Note that the initial data of (1.24) is zero and the external force can be small for small .
To estimate the right-hand-side of (1.24), we need bound for the composite operator for .
Since the Riesz operator is not bounded in , we estimate by direct calculations for their kernel.
The formula (2.10) is a key observation.
To estimate we need a control of since has a term such as for and .
In the semi-group approach, such kind of term cannot be estimated without additional regularity assumptions.
Note that if then for by proof of [8] and the integral formulation for the primitive equations.
For this reason we assume .
To prove our main theorem, we first show short time well-posedness to (1.18) and obtain the estimate (1.1) in a short interval for .
We can extended this solution to the interval since is also small.
Since is finite, we can establish global well-posedness to (1.18) for small .
This argument is used in [3] and [4].
We introduce some notation in this paper.
For , a domain , and , we write to denote the Lebesgue space equipped with the norm
|
|
|
We use the usual modification when .
For , domains and , and , we write to denote anisotropic Lebesgue spaces equipped with the norm
|
|
|
If and are the horizontal and vertical variable, respectively, then we write to denote .
For vector fields and , we denote their tensor product by .
For an integrable function on , we denote its vertical average by .
We write and to denote the Fourier and Fourier inverse transform for a integrable function , respectively.
We denote by the horizontal Laplace operator.
This paper is organized as follows.
In Section 2, we show linear estimate for the heat semi-group and the composite operator .
We also show some estimates for composite operators having fractional derivatives.
In Section 3 and Section 4, we show non-linear estimates and the external force , respectively, from the linear estimates of Section 2.
In Section 5, we prove our main theorem.
2 Linear Estimates
We give --estimate for the hear semi-group.
The reader refers to Section 4 of Grafakosβs book [9] for properties of the heat semi-group on .
Let be the heat kernel on for such that
|
|
|
where is the Gaussian of the form
|
|
|
For a integrable function , we denote by the heat semi-group on .
It is known that
|
|
|
|
|
|
for any multi-index and some constant , see [9] for the proof.
We find from interpolation inequalities that
|
|
|
(2.1) |
for all .
We next consider the composite operator with fraction derivative for and a integrable function on .
We write to denote the kernels of the corresponding composite operator on , respectively, namely
|
|
|
We know from Proposition 4.2 of Giga et al. [8] that
|
|
|
for some constant .
Thus we obtain
|
|
|
for all .
The above observations and the Poisson summation formula yield
|
|
|
(2.2) |
The same way as above we define the horizontal and vertical composite operators and .
It worth mentioning that, in [8], they define the vertical operator using the Caputo derivative since they consider the well-posedness in the anisotropic domain .
However, we do not use the Caputo derivative since we consider in .
Proposition 2.1.
Let .
Let () be a multi-index and .
Then there exists a constant such that
|
|
|
|
(2.3) |
|
|
|
|
(2.4) |
|
|
|
(2.5) |
for all and .
Proof.
The Young inequality and (2.1) imply
|
|
|
|
|
|
|
|
|
|
|
|
In view of (2.1) and (2.2), the inequalities (2.3) and (2.5) can be proved by the same way as above.
β
Let be the anisotropic Helmholtz projection on with the matrix-valued symbol
|
|
|
where and denotes the symbol of a multiplier operator .
The projection is an unbounded operator on since the Riesz operator is unbounded on for all .
However, the composite operators for and are bounded on .
We can rewrite the anisotropic Helmholtz projection as
|
|
|
(2.10) |
This is a key formula to show the boundedness for .
We denote by the anisotropic Riesz operator with symbol for .
We write .
We prove elementally estimates, which is used to show -independent bounds for composite operators.
It may be somewhat prolix, but we show calculation to clarify dependence of for the estimates.
Proposition 2.2.
Let and .
Then there exists a constant such that
|
|
|
(2.11) |
|
|
|
(2.12) |
|
|
|
(2.13) |
|
|
|
(2.14) |
for all and .
Proof.
We first prove (2.11).
The change of variable and the inequality yield
|
|
|
|
|
|
|
|
For (2.12) we divide the integral interval to see that
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.15) |
Similar to the first inequality, we use the change of variable to estimate
|
|
|
|
|
|
|
|
|
|
|
|
We proved (2.13).
For the last inequality, we apply the change of variables and the estimate to get
|
|
|
|
|
|
|
|
|
|
|
|
We obtain (2.14).
β
Proposition 2.3.
Let , and .
Then there exists a constant such that
|
|
|
|
(2.16) |
|
|
|
|
(2.17) |
|
|
|
|
(2.18) |
for all , , , and .
Proof.
We prove the first inequality.
The second and third inequalities can be proved by the completely same way combining with Propositions 2.1 and 2.2.
The formula (2.10) and
|
|
|
from Section 4 of [8], where and is the gamma function, lead to
|
|
|
|
|
|
|
|
|
|
|
|
The operator norm of the first term from to is bounded by Proposition 2.1.
We use Proposition 2.3 to estimate
|
|
|
|
|
|
We use Proposition 2.2 to see that
|
|
|
(2.19) |
Proposition 2.2 leads to
|
|
|
|
|
|
Thus, Proposition 2.2 implies
|
|
|
(2.20) |
The conclusion follows from (2.19) and (2.20).
We proved Proposition 2.3.
β
Corollary 2.5.
Let and .
Then there exists a constant such that
|
|
|
|
for all , , and .
Appendix B Decomposition of the solution around initial time
In this appendix we briefly show that the solution to (1.5) can be decomposed such that (4.1) and (4.2).
The proof is quite similar to the proof Theorem 1.1, we do not repeat the things for simplicity.
We decomposed the initial data satisfying the assumptions of Theorem 1.1 such that
|
|
|
satisfying
|
|
|
(B.1) |
for some constant , and small .
We know that there exist and a unique solution to (1.5) satisfying
|
|
|
(B.2) |
for some constant and small .
The reader refers to [5].
Let be the hydrostatic Helmholtz projection on .
Proposition B.1.
Let , , and .
Then there exits a constant such that
|
|
|
(B.3) |
for all two-dimensional vector fields and divergence-free satisfying .
Proof.
The proof is essentially same as Lemma 6.1 of [8].
We know from the Lemma that
|
|
|
(B.4) |
We apply the interpolation inequality to find
|
|
|
|
|
|
|
|
|
|
|
|
The seconde term in (B.4) is bounded as
|
|
|
|
|
|
|
|
|
|
|
|
and
|
|
|
|
|
|
|
|
|
Thus we have (B.3).
Note that the Proposition B.1 also holds if we change the role of and .
β
We now show the decomposition (4.1).
We only consider the case in (4.1) for simplicity.
Since and also has more regularity for the horizontal direction, it not difficult to improve the regularity to the case .
Put
|
|
|
where and is give by (1.6).
To show the decomposition, it is enough to show that there exists a solution to the equation satisfying the second estimate of (4.2).
We apply Proposition 6.2 in [8], see also the proof of Theorem 2.1, to get
|
|
|
|
|
|
and
|
|
|
|
|
|
Proposition B.1 and (B.2) imply
|
|
|
|
|
|
and
|
|
|
|
|
|
|
|
|
|
|
|
where we took for the first and second estimates, respectively.
If we set
|
|
|
the above estimates lead the quadratic estimate
|
|
|
for some constants and small .
If we take and sufficiently small beforehand, can be bounded small for small .
Since this argument is same as the proof of Theorem 1.1, we omit details here.
By the similar way, we see that becomes a contraction mapping for small .
Thus we can obtain the desired solution by the contraction mapping principle.