Abstract
This work is devoted to establishing the local-in-time well-posedness of strong solutions to the three-dimensional compressible primitive equations of atmospheric dynamics. It is shown that strong solutions exist, are unique, and depend continuously on the initial data, for a short time in two cases: with gravity but without vacuum, and with vacuum but without gravity.
Similar content being viewed by others
1 Introduction
1.1 The Compressible Primitive Equations
The general hydrodynamic and thermodynamic equations (see, e.g., [36]) with Coriolis force and gravity are used to model the motion and state of the atmosphere, which is a specific compressible fluid. However, such equations are extremely complicated and prohibitively expensive computationally. However, since the vertical scale of the atmosphere is significantly smaller than the planetary horizontal scale, the authors in [17] take advantage, as it is commonly done in planetary scale geophysical models, of the smallness of this aspect ratio between these two orthogonal directions to formally derive the compressible primitive equations (CPE) from the compressible Navier–Stokes equations. Specifically, in the CPE the vertical component of the momentum in the compressible Navier–Stokes equations is replaced by the hydrostatic balance equation (1.1)\(_{3}\), below, which is also known as the quasi-static equilibrium equation. It turns out that the hydrostatic approximation equation is accurate enough for practical applications and has become a fundamental equation in atmospheric science. It is the starting point of many large scale models in the theoretical investigations and practical weather predictions (see, e.g., [35]). This has also been observed by meteorologists (see, e.g., [39, 44]). In fact, such an approximation is reliable and useful in the sense that the balance of gravity and pressure dominates the dynamic in the vertical direction and that the vertical velocity is usually hard to observe in reality. In many simplified models, it is assumed that the atmosphere is under adiabatic process and therefore the entropy remains unchanged along the particle path. In particular, if the entropy is constant in the spatial variables initially, it remains so in later time. On the other hand, instead of the molecular viscosity, eddy viscosity is used to model the statistical effect of turbulent motion in the atmosphere. The observations above and more perceptions from the meteorological point of view can be found in [39, Chapter 4]. Therefore, under the above assumptions, one can write down the isentropic compressible primitive equations as in (1.1), below. Moreover, we also study the problem by further neglecting the gravity in (1.2), below. We remark here that, although it does not cause any additional difficulty, we have omitted the Coriolis force in this work for the convenience of presentation. That is, the local well-posedness theorems still work for systems (1.1) and (1.2) with the Coriolis force.
The first mathematical treatment of the compressible primitive equations (CPE) can be tracked back to Lions, Temam and Wang [35]. Actually, the authors formulated the compressible primitive equations in the pressure coordinates (p-coordinates) and show that in the new coordinate system, the equations are in the form of classical primitive equations (called primitive equations, or PE hereafter) with the incompressibility condition. In yet another work [34], the authors modeled the nearly incompressible ocean by the PE. It is formulated as the hydrostatic approximation of the Boussinesq equations. The authors show the existence of global weak solutions and therefore indirectly study the CPE (see, e.g., [32, 33] for additional work by the authors). Notably, the PE have been the subject of intensive mathematical research. For instance, Guillén-González, Masmoudi and Rodríguez-Bellido in [21] study the local existence of strong solutions and global existence of strong solutions for small initial data to the PE. In [47] the authors address the global existence of strong solutions to PE in a domain with small depth for restricted large initial data depending on the depth. In [38], the authors study the Sobolev and Gevrey regularity of the solutions to PE. The first breakthrough concerning the global well-posedness of PE is obtained by Cao and Titi in [8], in which the authors show the existence of unique global strong solutions (see, also, [9, 22,23,24,25, 27,28,29, 46] and the references therein for related study). On the other hand, with partial anisotropic diffusion and viscosity, Cao, Li and Titi in [3,4,5,6,7, 10] establish the global well-posedness of strong solutions to PE. For the inviscid primitive equations, or hydrostatic incompressible Euler equations, in [1, 26, 37], the authors show the short time existence of solutions in the analytic function space and in \(H^s \) space. More recently, the authors in [2, 45] construct finite-time blowup for the inviscid PE in the absence of rotation. Also, in [20], the authors establish the Gevrey regularity of hydrostatic Navier–Stokes equations with only vertical viscosity.
Despite the fruitful study of the primitive equations, it still remains interesting to study the compressible equations. On the one hand, it is a more direct model to study the atmosphere and perform practical weather predictions. On the other hand, the former deviation of the PE from the CPE in the p-coordinates did not treat the corresponding derivation of the boundary conditions. In fact, due to the change of pressure on the boundary, the appropriate studying domain for the PE should be evolving together with the flows in order to recover the solutions to the CPE. Thus, even though the formulation of the PE significantly simplifies the equations of the CPE, the boundary conditions are more complicated than before in order to study the motion of the atmosphere. We believe that this might be one of the reasons that is responsible for the not-completely successful prediction of the weather by using the PE.
Recently, Gatapov, Kazhikhov, Ersoy, Ngom construct a global weak solution to some variant of two-dimentional compressible primitive equations in [16, 19]. Meanwhile, Ersoy, Ngom, Sy, Tang, Gao study the stability of weak solutions to the CPE in [17, 41] in the sense that a sequence of weak solutions satisfying some entropy conditions contains a subsequence converging to another weak solution. In recent work, we show the existence of such weak solutions in [31]. See also [43].
In this and subsequent works, we aim to address several problems concerning the compressible primitive equations. In this work, we start by studying the local well-posedness of strong solutions to the CPE. That is, we will establish the local strong solutions to (1.1) and (1.2), below, in the domain \( \Omega = \Omega _h \times (0,1) \), with \( \Omega _h = {\mathbb {T}}^2 = [0,1]^2 \subset {\mathbb {R}}^2 \) being the fundamental periodic domain. In comparison with the compressible Navier–Stokes equations [18], the absence of evolutionary equations for the vertical velocity (vertical momentum) causes the main difficulty. This is the same difficulty as in the case of the PE. In fact, the procedure of recovering the vertical velocity is a classical one in the modeling of the atmosphere [39, Chapter 5]. This is done with the help of the hydrostatic equation, which causes the stratification of density profiles in the CPE. On the one hand, in (1.1), as one will see later, the hydrostatic equation implies that if there is vacuum in the physical domain \( \Omega \), the sound speed will be at most 1/2-Hölder continuous. Thus the \( H^2 \) estimate of the density is not available in the presence of vacuum. However in (1.2), such an obstacle no longer exists. For this reason, the local well-posedness established in this work doesn’t allow vacuum in the presence of gravity, but vacuum is allowed in the case without gravity. On the other hand, the hydrostatic equation does have some benefits. Indeed, such a relation yields that the density admits a stratified profile along the vertical direction. This fact will help us recover the vertical velocity from the continuity equation (see (1.8) and (1.13), below).
In this work, we will first reformulate the compressible primitive equations (1.1), (1.2) by making use of the stratified density profile. Then we will study the local well-posedness of the reformulated systems under the assumption that there is no vacuum initially. This is done via a fixed point argument. Next, in order to obtain the existence of strong solutions to (1.2) with non-negative density, we establish some uniform estimates independent of the lower bound of the density. We point out that in comparison to the compressible Navier–Stokes equations (see, e.g., [12,13,14,15]), we will require \( H^2 \) estimate of \( \rho ^{1/2} \) in order to derive the above mentioned uniform estimates. Such estimates are not available in the case with gravity (1.1). To this end, continuity arguments are used to establish the solutions with vacuum. We also study the continuous dependence on the initial data and the uniqueness of the strong solutions.
Through out this work, we will use \( {x} := (x,y,z)^\top , {x}_h := (x,y)^\top \) to represent the coordinates in \( \Omega \) and \( \Omega _h \), respectively. In addition, we will use the following notations to denote the differential operators in the horizontal direction:
The isentropic compressible primitive equations with gravity are governed by the system
with \( P := \rho ^\gamma \). We will study in this work only the case when \( \gamma = 2 \) in (1.1) for the sake of simplifying our presentation. For general \( \gamma > 1 \), we refer to Remark 1, below.
On the other hand, the isentropic compressible primitive equations without gravity are governed by the system
with \( P := \rho ^\gamma \) and \( \gamma > 1 \).
In the above systems, (1.1) and (1.2), the viscosity coefficients \( \mu , \lambda \) are assumed to be strictly positive. Also, (1.1) and (1.2) are supplemented with the following boundary conditions:
The rest of this paper will be organized as follows: in section 1.2, we present a reformulation of (1.1) and (1.2) by making use of the stratified density profiles. Also, we present the formula for recovering the vertical velocity and the main theorems of this work. After listing some useful inequalities and notations, we study in section 2 the existence theory. Next, in section 3 we show the continuous dependence on the initial data and the uniqueness of strong solutions.
1.2 Reformulation, Analysis and Main Theorems
In this section, we will reformulate (1.1) and (1.2) and point out how to recover the vertical velocity in terms of the density and the horizontal velocity.
1.2.1 The Case with Gravity and \( \gamma = 2 \)
We first consider (1.1). From (1.1)\(_{3}\), one has
Denote by \( \xi = \xi ({x}_h,t) : = \rho ^{\gamma -1}({x}_h,0,t) \). The continuity equation (1.1)\(_{1}\) implies
In particular, since \( \gamma = 2 \), we have \( \rho ({x},t) = \xi ({x}_h,t) + \dfrac{1}{2} gz \) and
(1.1) can be written as
Hereafter, we denote, for any \( f :\Omega \mapsto {\mathbb {R}} \),
Then averaging over the vertical variable in (1.5)\(_{1}\) yields, thanks to (1.3),
Then comparing (1.7) with (1.5)\(_{1}\) implies
Therefore, the vertical velocity w is determined, thanks to the boundary condition (1.3), by the relation
System (1.5) is complemented with the initial data
with \( \xi _0, v_0 \in H^2(\Omega ) \). Also the following compatible conditions are imposed:
Also, we will denote the bounds
Theorem 1
Suppose the initial data \( (\rho _0, v_0) = (\xi _0 + \frac{1}{2}gz, v_0) \) satisfy (1.11) and the compatible conditions (1.10). Then there is a unique strong solution \( (\rho ,v) \) to system (1.1), with the boundary condition (1.3), in \( \Omega \times (0,T) \), for some positive constant \( T = T(B_{g,1},B_{g,2},{{\underline{\rho }}}) > 0 \). Also, the solution satisfies
Furthermore, for some positive constant \( {\mathcal {C}}(B_{g,1},B_{g,2},{{\underline{\rho }}}) \),
Moreover, for any two solutions \( (\rho _i, v_i), i = 1,2 \) with initial data \( (\rho _{i,0},v_{i,0}) , i = 1,2 \) satisfying the conditions mentioned above, we have the following inequality
for some positive constant \( C_{\mu , \lambda , B_{g,1}, B_{g,2},{{\underline{\rho }}}, T} \).
Remark 1
For general \( \gamma > 1 \), after multiplying (1.4) with
and averaging the resultant in the z-variable, one obtains
Consequently, by eliminating \( \partial _t\xi \) from the above equation and (1.4), it follows that
Therefore, from above, and as in the case when \( \gamma =2 \), the vertical velocity w can be represented in the form
similarly to (1.8), for an explicit function \(H(\cdot ) \). Notably, the arguments and proofs, below, apply equally, and similar conclusion of Theorem 1 also holds for \( \gamma > 1 \).
1.2.2 The Case Without Gravity and \( \gamma > 1 \)
Concerning system (1.2), since (1.2)\(_{3}\) already yields the independence of the density of the vertical variable, after taking the vertical average of (1.2)\(_{1}\), as before, one has
Comparing (1.12) with (1.2)\(_{1}\) yields, thanks to the boundary condition (1.3), that the vertical velocity w is determined by the relation
In particular, by denoting \( \sigma := \rho ^{1/2} \), from (1.12) and (1.13), one has either \( \sigma = 0 \) or
In fact, for \( (\sigma , v) \) regular enough, (1.14), (1.15) hold regardless of whether \( \sigma = 0 \) or not. See also the justification in the beginning of section 3.2.
System (1.2) is complemented with the initial data
with \( \sigma _0 = \rho _0^{1/2}, v_0 \in H^2(\Omega ) \), and the initial total mass and physical energy satisfy
Also the following compatible conditions are imposed:
Also, we will denote the bounds
Moreover, if \( \rho = \sigma ^2 > 0 \), (1.2) can be written as
Theorem 2
Suppose the initial data \( (\rho _0,v_0) = (\sigma _0^2, v_0) \) satisfy (1.17), (1.19) and the compatible conditions (1.18). Then there is a unique strong solution \( (\rho , v) \) to system (1.2), with the boundary condition (1.3), in \( \Omega \times (0,T^*) \), for some positive constant \( T^* = T^*(B_1,B_2) > 0 \). Also, the solution satisfies
Furthermore, for some positive constant \( {\mathcal {C}}(B_1,B_2) \),
Moreover, for any two strong solutions \( (\rho _i, v_i), i = 1,2 \), with initial data \( (\rho _{i,0},v_{i,0}) , i = 1,2 \), satisfying the conditions mentioned above, we have the inequality
for some positive constant \( C_{\mu , \lambda , B_{1}, B_{2}, T^*} \).
1.3 Preliminaries
We will use \( \bigl | \cdot \bigr |_{}, \bigl \Vert \cdot \bigr \Vert _{} \) to denote norms in \( \Omega _h \subset {\mathbb {R}}^2 \) and \( \Omega \subset {\mathbb {R}}^3 \), respectively. After applying Ladyzhenskaya’s and Agmon’s inequalities in \( \Omega _h \) and \( \Omega \), directly we have
for any function f with bounded right-hand sides. Also, \( \bigl | {{\overline{f}}} \bigr |_{L^{p}}, \bigl \Vert {{\widetilde{f}}} \bigr \Vert _{L^{p}} \leqq C \bigl \Vert f \bigr \Vert _{L^{p}} \), for every \( p \geqq 1 \). Considering any quantities A, B, we use the notation \( A \lesssim B \) to denote \( A \leqq CB \) for some generic positive constant C, which may be different from line to line. In what follows \( \delta , \omega > 0 \) are arbitrary constants which will be chosen later in the relevant paragraphs to be adequately small. \( C_q \) represents a positive constant depending on the quantity q. We will also need the following classical inequality:
Lemma 1
Let \( 2\leqq p \leqq 6 \), and \( \rho \geqq 0 \) such that \( 0< \int _{\Omega }\rho \,\hbox {d}{x}= M <\infty \), and \( \int _{\Omega }\rho ^\gamma \,\hbox {d}{x}\leqq E_0 \), for some \( \gamma \in (1, \infty ) \). Then one has
for some constant \( C = C(M, E_0) \), provided the right-hand side is finite.
Proof
This is standard. See, e.g., [18, Lemma 3.2]. \(\quad \square \)
2 Associated Linear Systems and Existence Theory
In this section, we will establish the local existence theory of (1.1) and (1.2). To do this, we will first study the local existence of solutions to (1.5) and (1.20) via the Schauder–Tchonoff fixed point theorem capitalizing on some a priori estimates. In fact, under the assumption that
we will first introduce linear systems and the function spaces \( {\mathfrak {Y}} \) associated with (1.5) and (1.20) with some given input states \( (\xi ^o,v^o) \) and \((\sigma ^o,v^o)\), respectively, in section 2.1 and 2.2. Here, \( {\mathfrak {Y}} \) are compactly embedded in some corresponding spaces \( {\mathfrak {V}} \). Also, we will show that the maps \( {\mathcal {T}}: {\mathfrak {X}} \mapsto {\mathfrak {X}} \), for some convex bounded subsets \( {\mathfrak {X}} \) of \( {\mathfrak {Y}} \), given by
are well-defined; observing that \( {\mathfrak {X}} \) are convex subsets of \({\mathfrak {Y}} \) and hence compact in \( {\mathfrak {V}} \). We will use the same notations \( {\mathfrak {X}}, {\mathfrak {Y}}, {\mathfrak {V}}, {\mathcal {T}} \) to denote the convex bounded sets, the compact function spaces, the embedded function spaces and the constructed maps in both cases. We summarize the relevant regularity estimates in section 2.3 and show that the Schauder-Tchonoff fixed point theorem will yield the existence of solutions to (1.5) and (1.20) in the corresponding set. Recall that the Schauder-Tchonoff fixed point theorem states that for a Banach space V with a convex compact subset \( X \subset V \), if \(F : X \mapsto X\) is continuous, then F has at least one fixed point in X. In our case, we will take \( X = {\mathfrak {X}} \) and \( V = {\mathfrak {V}} := \lbrace (\xi ,v)| \xi , v \in L^\infty (0,T;L^2(\Omega )), \nabla v\in L^2(0,T;L^2(\Omega )) \rbrace \) in the case with gravity, or \( V= {\mathfrak {V}} := \lbrace (\sigma ,v)| \sigma , v \in L^\infty (0,T;L^2(\Omega )), \nabla v\in L^2(0,T;L^2(\Omega )) \rbrace \) in the case without gravity, with the corresponding norms.
We will only sketch the key steps in this paper. For more detailed calculation, we refer to our preprint [30].
2.1 The Case with Gravity and \( \gamma = 2 \)
2.1.1 Associated Linear Inhomogeneous System
Consider a finite positive time T, which will be determined later. Let \( {\mathfrak {Y}} = {\mathfrak {Y}}_T \) be the function space defined by
with the norm
Notice that, thanks to the Aubin compactness theorem (see, e.g., [42, Theorem 2.1] and [11, 40]), every bounded subset of \( {\mathfrak {Y}} \) is a compact subset of the space
with the norm
Let \( {\mathfrak {X}} = {\mathfrak {X}}_{T} \) be a bounded subset of \( {\mathfrak {Y}} \) defined by
where \( M_0, M_1 \) are the bounds of initial data in (2.9) and \(C_1 = C_1(M_0,\mu ,\lambda ,{{\underline{\rho }}})\), \( C_2 = C_2(M_0,C_1M_1) \) are given below in (2.34), (2.17), respectively. Notice, for \( (\xi , v) \in {\mathfrak {X}} \),
Let \( (\xi ^o,v^o) \in {\mathfrak {X}} \). The following inhomogeneous linear system is inferred from (1.5) using \((\xi ^o,v^o) \) as an input:
Here \( w^o \) is given by (1.8) with \( (\xi ^o, v^o) \) instead of \( (\xi , v) \), i.e.,
Notice that (2.6)\(_{1}\) is inferred from (1.7). For details, see the deviation from (1.5) to (1.8). Hereafter, denote by \( \rho ^o := \xi ^o + \dfrac{1}{2}gz \) . The initial and boundary conditions for the linear system (2.6) are given by
The compatible conditions in (1.10) are still imposed and we require
Recall that \( V_1 \) is given in (1.10), essentially, \( V_1 = v_t|_{t=0} \).
Then the map \( {\mathcal {T}} \), in this case, is defined as
where \( (\xi , v) \) is the unique solution to the linear system (2.6) with \( (\xi ^o, v^o) \in {\mathfrak {X}} \). We claim that \( {\mathcal {T}} \) is a well defined map from \( {\mathfrak {X}} \) to \( {\mathfrak {X}} \), which is the consequence of the following two propositions:
Proposition 1
For given \( (\xi ^o,v^o) \in C^\infty ({\overline{\Omega }}\times [0,T]) \cap {\mathfrak {X}}_T \), there is a unique strong solution \( (\xi , v) \in {\mathfrak {Y}}_T \) of system (2.6) with the initial and boundary conditions (2.8).
Suppose, in addition, that \( \xi _0, v_0 \in H^3(\Omega ) \). One will have the following regularity of the unique solution \( (\xi , v) \) of system (2.6):
Proposition 2
Consider the initial data with the bounds \( M_0, M_1 \) in (2.9) and \( (\xi ^o, v^o) \in {\mathfrak {X}} = {\mathfrak {X}}_T \). There is a \( T_g = T_g(M_0,M_1,\mu ,\lambda ,{{\underline{\rho }}}) > 0 \) sufficiently small such that for any \( T \in (0, T_g) \), there exists a unique solution to (2.6). Moreover, the solution belongs to \( {\mathfrak {X}} = {\mathfrak {X}}_T \). Therefore, for any such T, the map \( {\mathcal {T}} \) in (2.10) is a well defined map from \( {\mathfrak {X}} \) into \( {\mathfrak {X}} \).
We omit the proof of Proposition 1, and refer to [30] for the details. The existence of solutions in Proposition 2 follows from Proposition 1 and a standard approximating argument. We only show the required a priori estimates in the rest of this subsection, which are sufficient to establish Proposition 2. See Propositions 3 and 4, below.
2.1.2 A Priori Estimates for the Inhomogeneous Linear System
Hereafter, we assume that the solution \( (\xi , v) \) to the linear system (2.6) is smooth enough so that the following estimates are rigorous.
We start by establishing some estimates for the solutions of (2.6)\(_{1}\). In particular, we will establish the following:
Proposition 3
There exists a \( T'= T'(M_0,C_1 M_1,{{\underline{\rho }}}) > 0 \) sufficiently small such that for any \( T \in (0,T'] \), the solution \( \xi \) to (2.6)\(_{1}\) satisfies that
where \( M_0 \) is as in (2.9).
The lower bound for \( \xi \)
In order to derive the lower bound of \( \xi \), we employ the following Stampaccia-like argument. Let \( M = M(t) > 0 \) be a nonnegative integrable function to be determined later. Consider \( \eta = \eta (x,y,t) := \xi - {{\underline{\rho }}} + \int _0^t M(s) \,\hbox {d}s \). Then according to (2.6)\(_{1}\), \( \eta \) satisfies the equation
Let
and denote by \( \eta _- := - \eta \mathbb {1}_{\lbrace \eta < 0 \rbrace } \geqq 0 \). Observe that since \( \xi \in H^1(\Omega \times [0,T]) \), so is \( \eta _- \). Thus, one has
Now, let \( 0< M(t) := C \max \lbrace \bigl | \overline{\mathrm {div}_h\,v^o} \bigr |_{L^{\infty }}, \bigl | \overline{z \mathrm {div}_h\,v^o} \bigr |_{L^{\infty }} \rbrace \leqq C \bigl \Vert v^o \bigr \Vert _{H^{3}} < \infty , a.e.,\), for some constant \( C>0 \). Then the integrand on the right-hand side of the above equation satisfies
provided C is large enough and T is small enough such that
Therefore, we have
which, after integrating over \( [0,t_0] \) for any \( t_0 \in [0,T] \), thanks to the fact \( \eta _-(0) \equiv 0 \), yields
Hence \( \eta _- = 0 ~ \text {in} ~ \Omega _h \times [0,T] \). That is, \( \eta (t) = \xi (t) - {{\underline{\rho }}} + \int _0^t M(s) \,\hbox {d}s \geqq 0 \) and
for \( t \in [0,T], ~ T \leqq T_1 \), with \( T_1= T_1(C_1M_1, {{\underline{\rho }}}) \) sufficiently small.
The \( H^2(\Omega ) \) norm for \(\xi \)
Applying the standard \( H^2 \) estimate of linear transport equations to (2.6)\(_{1}\) yields, since \( \bigl | \xi \bigr |_{H^{2}} = \bigl \Vert \xi \bigr \Vert _{H^{2}} \),
Thanks to the Grönwall and the Hölder inequalities, one has
Then for \( T \in (0, T_2] \) with \( T_2(M_0,C_1M_1) \) sufficiently small, (2.14) yields
The \( H^1(\Omega ) \) norm for \( \partial _t\xi \) Using equation (2.6)\(_{1}\), \( \partial _t\xi \) can be represented in terms of the spatial derivatives of \( \xi , v^o \). Then applying the Hölder and the Sobolev embedding inequalities implies
where \( C_2 = C_2(M_0,C_1M_1) \) is given by
Proof of Proposition 3
By choosing \( T' = \min \lbrace T_1,T_2 \rbrace \), the proof of Proposition 3 follows from (2.13), (2.15) and (2.16). \(\quad \square \)
Next, we show
Proposition 4
There exists a \( T''=T''(M_0,M_1,C_1,C_2,{{\underline{\rho }}}) \in (0,\infty ) \), sufficiently small, such that for every \( T \in (0, T''] \), the solution v to (2.6)\(_{2}\) satisfies
Horizontal spatial derivative estimates for v
Applying \( \partial _{hh} \) to (2.6)\(_{2}\) will yield
where \( \rho ^o = \xi ^o + \frac{1}{2} gz \). After taking the inner product of (2.18) with \( \partial _{hh} v \) and integrating by parts, one has
While we only will omit the detailed estimates, which are standard, we list below the estimates for the \( I_i\) terms (see [30] for details). We will use the fact \( \bigl \Vert \rho ^o \bigr \Vert _{H^{2}}^2 \leqq C \bigl \Vert \xi ^o \bigr \Vert _{H^{2}}^2 + C g^2 \leqq C M_0 + C \), \( \bigl \Vert \partial _t\rho ^o \bigr \Vert _{L^{2}}^2 = \bigl \Vert \partial _t\xi ^o \bigr \Vert _{L^{2}}^2 \leqq C_2\). Also, hereafter the estimates hold for every \( \delta , \omega >0 \) which will be chosen later to be adequately small. Correspondingly, \( C_{\delta }, C_{\omega }, C_{\delta ,\omega } \) are some positive constants depending on \( \delta , \omega \):
In order to estimate \( I_4 \), we shall plug in (2.7). One has
Then applying the Minkowski and the Sobolev embedding inequalities yields
and, similarly,
where we have employed (1.21). Summing up the above inequalities, with \( \delta \) small enough, yields the following estimate
Hereafter, \( {\mathcal {H}} \) will be used to denote a polynomial quantity of its arguments (i.e., the norms of the initial data and \(\xi ^o, v^o\)) which may be different from line to line. Also \( c_{\mu ,\lambda }, C_\omega \) denote positive constants depending on \( \mu , \lambda \) and \( \omega \), respectively. Similar arguments also hold for the lower order derivatives. Then after suitable choice of \( \omega \), one has
Time derivative estimates for v
Applying \( \partial _t\) to (2.6)\(_{2}\) yields
Consequently, one has
Then applying similar estimates as before to the right-hand side of (2.22) (see [30] for details) and making use of the fact that \( w^o \) is given by (2.7), one has
Vertical derivative estimates for v
Taking the inner product of (2.6)\(_{2}\) with \( v_t \) and applying similar estimates as before to the resultant implies
On the other hand, (2.6)\(_{2}\) can be written as
Taking the inner product of (2.25) with \( \partial _t\partial _{zz}v \) and applying similar estimates as before will yield,
Next, applying \( \partial \in \lbrace \partial _x, \partial _y, \partial _z \rbrace \) to (2.25) yields
This implies
Also, by employing the Minkowski’s inequality, one obtains
and, similarly,
Therefore, we have
Proof of Proposition 4
From (2.20), (2.23), (2.24) and (2.26), there is a constant \( c_{\mu , \lambda , {{\underline{\rho }}}} \) such that
where
Notice that for some positive constants \( C_{i,\mu , \lambda , {{\underline{\rho }}}, M_0}, i = 1,2 \), depending on \( \mu , \lambda , {{\underline{\rho }}}, M_0 \), we have
For \( 0 < \omega \leqq \frac{c_{\mu ,\lambda ,{{\underline{\rho }}}}}{2}\), one infers from (2.29), that
Therefore, applying the Grönwall’s inequality yields
Now, we integrate with respect to the time variable inequality (2.29). It follows, since \( 0< \omega < \frac{c_{\mu ,\lambda ,{{\underline{\rho }}}}}{2} \), that
Additionally, from (2.28), we have
Therefore, we conclude that
for some positive constant \( C_{3,\mu ,\lambda ,{{\underline{\rho }}},M_0} \), depending on \( \mu , \lambda , {{\underline{\rho }}},M_0 \). Now fix \( \omega = \frac{1}{2} \min \lbrace \frac{c_{\mu ,\lambda ,{{\underline{\rho }}}}}{2}, \frac{1}{C_1} \rbrace \) and let \( T \in (0,T''] \), where \( T'' = T''(M_0,M_1,C_1,C_2,{{\underline{\rho }}}) \) is small enough and satisfying
Then (2.32) yields
where \( C_1 \) is given by
This concludes the proof. \(\quad \square \)
2.2 The Case Without Gravity and \(\gamma > 1\)
Consider a finite positive time T, which will be determined later. Let \( {\mathfrak {Y}} = {\mathfrak {Y}}_T \) be the function space defined by
Notice that, thanks to Aubin compactness theorem (see [42, Theorem 2.1] and [11, 40]), every bounded subset of \( {\mathfrak {Y}} \) is a compact subset of the space
Let \( {\mathfrak {X}} = {\mathfrak {X}}_T \) be a bounded subset of \( {\mathfrak {Y}} \) defined by
where \({{\underline{\rho }}}\) is the positive lower bound of initial density profile as in (2.1), for some positive constants \(C_1 = C_1(M_0,\mu ,\lambda ,{{\underline{\rho }}})\), \( C_2 = C_2(M_0,C_1M_1)\). Notice, for \( (\sigma , v) \in {\mathfrak {X}} \), \( \int _0^1 \mathrm {div}_h\,(\sigma ^2 {\widetilde{v}}) \,\hbox {d}z =0. \) Let \( (\sigma ^o, v^o) \in {\mathfrak {X}} \). The linear system inspired by (1.20) with \( (\sigma ^o, v^o ) \) as an input is given by
where \( \rho ^o = (\sigma ^o)^2 \) and \( w^o \) is determined, as in (1.15), by
The initial and boundary conditions for system (2.38) are given by
Here, in addition to the compatible conditions in (1.18), we require \( \rho _0 \geqq {{\underline{\rho }}} > 0 \) , for some positive constant \( {{\underline{\rho }}} \) as in (2.1). Also, we denote by \( V_1 : = h_1 / \rho _0^{1/2} \). Recall that \( h_1 \) is given in (1.18). Then \( V_1 \in L^2 (\Omega ) \) and we require \( \bigl \Vert \sigma _0 \bigr \Vert _{H^{2}}^2 \leqq M_0, \bigl \Vert v_0 \bigr \Vert _{H^{2}}^2 + \bigl \Vert V_1 \bigr \Vert _{L^{2}}^2 \leqq M_1 \). Essentially \( V_1 = v_t|_{t=0} \).
Then the map \( {\mathcal {T}} \), in the case without gravity, is defined as
where \( (\sigma , v) \) is the unique solution to (2.38) for given \( (\sigma ^o, v^o) \in {\mathfrak {X}} \).
Proposition 5
There is a \( T_v = T_v(M_0,M_1,\mu ,\lambda ,{{\underline{\rho }}}) >0 \), sufficiently small, such that for every \( T \in (0, T_v] \), there is a unique solution \( ( \sigma , v) \) to (2.38) in the set \( {\mathfrak {X}} = {\mathfrak {X}}_T \). Therefore, for such T, the map \( {\mathcal {T}} \) defined in (2.41) is a well defined map from \( {\mathfrak {X}} \) into \( {\mathfrak {X}} \).
The proof is similar as Proposition 1 and Proposition 2 in section 2.1 and therefore is omitted.
2.3 Existence Theory
In this subsection, we will establish the existence of solutions for (1.1) and (1.2) for given corresponding initial data and boundary conditions.
2.3.1 The Case With Gravity and \( \gamma = 2 \), but Without Vacuum
We will apply the Schauder-Tchonoff fixed point theorem to establish the existence of strong solutions to (1.1). With Proposition 2, it is sufficient to verify that \( {\mathcal {T}} \), defined by (2.10), is continuous in \( {\mathfrak {V}} = {\mathfrak {V}}_T \) given in (2.3), where the norm is given by
In order to show this, let \( M_0 = B_{g,1} , M_1 = B_{g,2} \) in \( {\mathfrak {X}}_T \) and \( T \in (0, T_g] \), with \( T_g \) given in Proposition 2. Here \( B_{g,1}, B_{g,2} \) are given in (1.11). We denote \( (\xi _1^o,v_1^o), (\xi _2^o,v_2^o) \in {\mathfrak {X}}_T \) and
Denote by \( \xi _{12} : = \xi _1 - \xi _2, v_{12} := v_1 - v_2, \xi _{12}^o : = \xi _1^o - \xi _2^o, v^o_{12} := v_1^o - v_2^o \). Then \( (\xi _{12}, v_{12})|_{t=0} = 0 \). By taking the differences of the equations satisfied by \( (\xi _i, v_i), i = 1,2 \), we have
Now we perform standard \( L^2 \) estimates for (2.42). Take the \( L^2 \)-inner product of (2.42)\(_{1}\) with \( 2 \xi _{12} \) and take the \( L^2 \)-inner product of (2.42)\(_{2}\) with \( 2 v_{12} \). Then applying similar estimates as before yields, together with the Grönwall inequality, that
This implies the continuity of \( {\mathcal {T}} \) in \({\mathfrak {V}}\). Therefore, after applying the fixed point theorem mentioned before, we have the following:
Proposition 6
Consider
given in (1.9) satisfying (1.10) and (1.11). There is a positive constant T depending on the initial data such that there is a strong solution \( (\rho , v) = (\xi + \frac{1}{2}gz, v) \) to (1.1) (or equivalently (1.5)) with the boundary conditions (1.3) and with \( (\xi , v) \in {\mathfrak {X}}_T \).
2.3.2 The Case With Vacuum and \( \gamma > 1 \), but Without Gravity
When \( \rho _0 \geqq {\underline{\rho }} > 0 \), the existence of strong solutions to (1.2) follows from the estimates in section 2.2 and similar arguments to those in section 2.3.1. In fact, taking \( M_0 = B_1, M_1 = B_2 + {{\underline{\rho }}}^{-1} B_2 \), we have the following:
Proposition 7
Suppose that (1.17), (1.18), (1.19) hold for the given initial data (1.16) with \( \rho _0 \geqq {{\underline{\rho }}} > 0 \). Then there is a positive constant T, depending on the initial data and \( {{\underline{\rho }}} \), such that there exists a strong solution \((\rho , v) = (\sigma ^2, v) \) to (1.2) (or equivalently (1.20)) satisfying the boundary conditions (1.3) and that \( (\sigma ,v) \in {\mathfrak {X}}_T \).
In the following, we shall present some estimates independent of \({\underline{\rho }} \) and show that for a given non-negative initial density \( \rho _0 \geqq 0 \), there are strong solutions to equations (1.2). We will use here the notation \( \sigma ^2 = \rho \) and the alternative form of equations (1.20), as well as (1.2). Meanwhile, let us assume that
for given \( K_1, K_2 > 0 \). Recall essentially \( h_1 = (\sigma v_t)|_{t=0} \) from (1.18). Also, taking inner product of (1.2)\(_{2}\) with v yields, the conservation of physical energy
where (1.2)\(_{1}\) is also applied. Also, integrating (1.2)\(_{1}\) over \( \Omega \times (0,T) \) yields the conservation of total mass
These facts are important when applying (1.22) in what follows.
A priori assumptions
Let \( (\sigma , v) \) be the solution to (1.20) given in Proposition 7. We assume first, for some constants \( C_d \geqq K_2^2 \), \( T_d \) (may depend on \( {\underline{\rho }} \)),
In what follows, we will derive some a priori estimates independent of \( {\underline{\rho }} \). Also, we set \( T \in (0, T_d] \) to be determined later. We emphasize that the smallness of T in what follows is independent of \( {{\underline{\rho }}} \).
\( {\underline{\rho }} \)-independent lower bound: non-negativity of \( \rho \)
We will use the same Stampaccia-like argument as before to derive the lower bound of \( \rho \). Consider
Due to (1.12), \( \eta \) satisfies
for every \( t \in [0,T] \) with \( T \in (0,T_1] \), and \( T_1 \) sufficiently small such that
Denote by \( \eta _-: = - \eta \mathbb {1}_{\lbrace \eta < 0 \rbrace } \geqq 0 \). Then multiplying the above equation with \( - \mathbb {1}_{\lbrace \eta < 0 \rbrace } \) and integrating the resultant in the spatial variable yield
Hence, \( \eta _- = 0 \) in \( \Omega _h \times (0,T] \), since \( \eta _-(0) \equiv 0 \). Therefore, \( \eta \geqq 0 \) and
\({\underline{\rho }} \)-independent estimate: \( H^2(\Omega ) \) for \( \sigma = \rho ^{1/2} \)
After performing standard \( H^2 \) estimate of (1.14) and applying the Grönwall inequality to the result, one has
for all \( T \in (0,T_2] \), provided \( T_2 \) is sufficiently small.
\( {\underline{\rho }} \)-independent estimate: \( H^1(\Omega ) \) for \( \sigma _t \)
It directly follows from (1.14) that
for all \( t \in [0,T] \).
\( {\underline{\rho }} \)-independent estimate: \( L^2(\Omega ) \) for \( v_t \)
Taking the time derivative of (1.20)\(_{2}\) yields
Taking the \(L^2\)-inner product of (2.51) with \( \partial _tv \) gives
Then one has the following estimates to the terms in the right-hand side of (2.52) (see [30] for details):
We have applied above the Hölder inequality and (1.22). Notice that \( \sigma = \rho ^{1/2} \) and that the conservations of energy and mass (2.45), (2.46) hold. In order to estimate \( L_3 \) term, we substitute (1.13) and integrate by parts. Then
Now we use (1.21), the Minkowski and Hölder inequalities,
and, similarly,
After summing the above inequalities, (2.52) implies
where, as before, \( {\mathcal {H}} \) denotes a polynomial quantity of its arguments.
\( {\underline{\rho }} \)-independent estimate: spatial derivatives of v
Now we are able to derive the estimates on the spatial derivatives of v. Standard \( L^2 \) estimate of (1.20)\(_{2}\) yields that
Furthermore, taking the \( L^2 \)-inner product of (1.20)\(_{2}\) with \( v_t \) yields, after applying similar estimates as before,
Next we estimate the second order spatial derivatives. Taking the \( L^2 \)-inner product of (1.20)\(_{2}\) with \( \partial _{zz} v_t \) yields
At the same time, taking the \( L^2 \)-inner product of (1.20)\(_{2}\) with \( \Delta _hv_t \) yields
Applying similar estimates as before to the right-hand sides of (2.56) and (2.57) yields, after summing up the results,
Finally, we provide estimates for the third spatial derivative of v. Applying \( \partial \in \lbrace \partial _x, \partial _y,\partial _z \rbrace \) to (1.20)\(_{2}\) yields
Taking the \( L^2 \)-inner product of (2.59) with \( \Delta _h\partial v \), for \( \partial \in \lbrace \partial _h, \partial _z \rbrace \), and integrating by parts imply, after applying the Cauchy-Schwarz inequality,
Then, similarly to (2.28), noticing the fact that w are given by (1.15), (2.60) yields
where \( C_d \) is as in (2.47)
We summarize the estimates obtained, so far, in this section in the following:
Proposition 8
Consider the solution \((\sigma ,v) = (\rho ^{1/2},v) \) to (1.2) with the bound (2.47) and initial data satisfying (1.17), (2.44). There is a positive constant \( T^*=T^*(C_d,K_1,K_2) \), sufficiently small, such that \( (\sigma ,v) \) admits the following bounds, for \( T = \min \lbrace T^*, T_d \rbrace \),
where \( K_4 = \sqrt{2C_{\mu ,\lambda }} K_2, K_3 = CK_1K_4 \) are given in (2.63) and (2.64). Notably, the bounds in these estimates depend only on the initial bounds \( K_1, K_2 \) and do not depend on the lower bound of density. Also, the smallness of \( T^* \) does not depend on \( {{\underline{\rho }}} \), even though \( T_d \) may depend on \( {{\underline{\rho }}} \).
Proof
Denote by
Then from (2.53), (2.54), (2.55) and (2.58), we have
Then integrating the above inequality yields, for \( T \in (0, T_d] \), where \( T_d \) is as in (2.47),
Then together with (2.61), we have, after choosing \( \omega \) small enough and then T sufficiently small,
where we have employed inequality (1.22) and the fact for some positive constant \( C_{\mu ,\lambda } > 0 \) we have
Then plugging in (2.63) back into (2.50) implies
Thus the conclusion is drawn from (2.48), (2.49), (2.63) and (2.64). \(\quad \square \)
Existence of strong solutions with vacuum but no gravity and \( \gamma > 1 \)
Now we are in the place to remove the strict positivity (of the initial density profile) assumption in Proposition 7. In order to do so, we introduce a sequence of approximating initial data \( (\rho _{0,n}, v_{0,n} ) \) satisfying, in addition to (1.17), (1.18), (1.19),
such that
in \( H^2(\Omega ) \), as \( n \rightarrow \infty \), where \( (\rho _0, v_0 ) \) (or equivalently \( (\sigma _0, v_0) \)) is given in (1.16) satisfying (1.18) and (1.19).
We require that the initial physical energy and total mass given in (1.17) with \( \rho _0, v_0 \) replaced by \( \rho _{0,n}, v_{0,n} \) satisfy
uniformly in n, so that when we apply inequality (1.22), the constant in the inequality is independent of n.
Now we apply Proposition 7 with the initial data \( (\rho _{0,n},v_{0,n}) \). Indeed, consider \( M_0 = B_1 \), \( M_1 = B_2 + n B_2 \). Then Proposition 7 guarantees that there is a \( T_1 = T_1(n,B_1,B_2) \) such that (1.2) admits a strong solution \( (\rho _n, v_n) = (\sigma _n^2 , v_n) \) satisfying
Sequently, we apply Proposition 8 with \( K_1 = B_1^{1/2}, K_2 = 2B_2^{1/2} \), \( C_d = (1 + 2 n) {\mathcal {C}}_2(B_1,B_2,n) \) and \( T_d = T_1 \). It yields that there is a \( T_2 = T_2(B_1,B_2,n) \leqq T_1 \) such that the following bounds are satisfied
Next, let \( (\sigma _n,v_n)|_{t=T_2} \) as a new initial data for (1.2). The same arguments as above yield the bound (2.65) with lower bound of \( \rho _n \) replaced by \( \frac{1}{4n} \), \( B_1 \) replaced by \( 4 B_1 \) and \( B_2 \) replaced by \( {\mathcal {C}}_4(B_1,B_2) \). That is, for some \( \delta T =\delta T(B_1,B_2,n) > 0 \),
Now we apply Proposition 8 with \( T_d = T_2 + \delta T \) and \( C_d = {\mathcal {C}}_{4}(4B_1,{\mathcal {C}}_{4}(B_1,B_2)) \) in the time interval \( (0,T_d) \). This will yield that there is a \( T^* = T^*(B_1,B_2) \) and \( T_3 := \min \lbrace T^*, T_2 + \delta T \rbrace \), the bounds in (2.65) hold with \( T_2 \) replaced by \( T_3 \).
If \( T_3 = T^* \), we have got an existence time independent of n and this finishes the job. Otherwise, let \( (\sigma _n,v_n) |_{t= T_3} \) as a new initial data and repeat the arguments above to get the bounds in (2.65) with \( T_2 \) replaced by \( T_4:= \min \lbrace T^*, T_3 + \delta T\rbrace = \min \lbrace T^*, T_2 + 2 \delta T\rbrace \). Keep repeating this process, one will eventually get that there is a \( m \in {\mathbb {Z}}^+ \) sufficiently large that \( T_m := \min \lbrace T^*, T_2 + (m-2) \delta T \rbrace = T^* \).
Therefore, we have got a sequence of approximating solutions \( (\rho _n, v_n) = (\sigma _n^2, v_n) \) with a uniform existence time \( T^* \) independent of n for the approximating initial data \( (\rho _{0,n}, v_{0,n}) \) constructed above. In particular, \( (\sigma _n, v_n) \) satisfies the bounds in (2.65) with \( T_2 \) replaced by \( T^* \). Thus by taking \( n\rightarrow \infty \), it is straightforward to check that we have got a strong solution \( (\rho , v) = (\sigma ^2, v) \) to (1.2). In fact, we have the following:
Proposition 9
Consider the initial data \( (\rho _0, v_0) \) (or equivalently \( (\sigma _0, v_0) \)) given in (1.16) satisfying (1.17), (1.18) and (1.19). There is a constant \( T^* > 0 \) such that there exists a solution \( (\rho , v) = (\sigma ^{2}, v) \) to equation (1.2) satisfying
and
for some constant \( {\mathcal {C}}_3 = {\mathcal {C}}_3(B_1,B_2), {\mathcal {C}}_4 = {\mathcal {C}}_4(B_1,B_2) \).
3 Continuous Dependence on Initial Data and Uniqueness
In this section, we will show the continuous dependence of the solutions of (1.1) and (1.2) on the initial data. This will also imply the uniqueness of strong solutions constructed in Proposition 6 and Proposition 9.
3.1 The Case With Gravity and \( \gamma = 2 \), but Without Vacuum
Consider two sets of initial data \( (\rho _{i,0},v_{i,0}) = (\xi _{i,0}+\frac{1}{2}gz ,v_{i,0}), i = 1,2, \) in (1.9) for (1.1) satisfying (1.10), (1.11). Denote \( (\rho _i,v_i) = (\xi _i+\frac{1}{2} gz, v_i), i = 1,2, \) as the corresponding strong solutions constructed in Proposition 6 in the interval [0, T] for some \( T > 0 \). Then we have \((\xi _i, v_i) \in {\mathfrak {X}}_T, i=1,2 \). Throughout this section we will denote the constant \( C > 0 \) which may be different from line to line and depends on \( \mu , \lambda , B_{g,1}, B_{g,2},{{\underline{\rho }}}, T \). Also, we will use the notations
Taking the difference of the equations satisfied by \( (\xi _i, v_i), i=1,2 \), as in (2.42), then \( (\xi _{12}, v_{12}) \) satisfies
Then after applying standard \( L^2 \) estimates to the above system and applying the Grönwall’s inequality to the resultant, one can show the following:
Proposition 10
Given two sets of initial data \( (\rho _{i,0},v_{i,0}) = (\xi _{i,0}+\frac{1}{2}gz ,v_{i,0}), i = 1,2 \), satisfying (1.10) and (1.11), the corresponding strong solutions \( (\rho _i,v_i) = (\xi _i+\frac{1}{2} gz, v_i), i = 1,2 \), of (1.1) constructed in Proposition 6 in the interval [0, T] , for some \( T > 0 \), satisfy
In particular, if \( \rho _{1,0} = \rho _{2,0}, v_{1,0} = v_{2,0} \), we have \( \rho _1 =\rho _2, v_1 = v_2 \) in [0, T] .
3.2 The Case With Vacuum and \( \gamma > 1 \), but Without Gravity
First, we claim that any solution \( (\rho , v) = (\sigma ^2, v) \) to (1.2) satisfying (2.66) with the bounds in (2.67) will also satisfy the following equations:
To show this claim, we first consider the non-degenerate variable \( \rho + \varepsilon = \sigma ^2 + \varepsilon \), for some constant \( \varepsilon > 0 \). From (1.12), one has
Then after dividing \( (\rho + \varepsilon )^{1/2} \), one has
Now it is easy to verify that (3.1) will converge to 1.2’\(_{1}\) in the sense of distribution, as \( \varepsilon \rightarrow 0 \). On the other hand, from (1.13), one has
We define
Then \( \sigma \sigma w_\sigma = \rho w \) and we will use hereafter the notation \( \sigma w = \sigma w_\sigma \). As before it is easy to verify that 1.2’\(_{3}\) is equivalent to 1.2’\(_{2}\) in the sense of distribution. Summing up the facts above, we have shown that the solutions to (1.2) satisfying the (2.66) regularity with the bounds in (2.67) are also solutions to 1.2’.
Consider two sets of initial data \( ( \rho _{i,0},v_{i,0}) = (\sigma _{i,0}^2, v_{i,0}), i = 1,2 \), in (1.16) for (1.2) satisfying (1.18) and (1.19). Denote \( (\rho _i,v_i) = (\sigma _i^2, v_i), i = 1,2 \), as the corresponding strong solutions constructed in Proposition 9 in the interval \( [0,T^*] \), for some \( T^* > 0 \). Then we have \( (\sigma _i, v_i), i=1,2 \), satisfying the bounds in (2.67). Also \( (\sigma _i, v_i), i=1,2 \), are solutions to 1.2’. Throughout this section, we will denote the constant \( C > 0 \) which may be different from line to line and depends on \( \mu , \lambda , B_1,B_2, T^* \). Also, we will use the notations
Taking the difference of the equations satisfied by \( (\sigma _i, v_i), i=1,2 \), we have
Then as before, one can derive
Therefore, after employing (1.22) and noticing the fact that we can interchange \( (\sigma _1, v_1), (\sigma _2,v_2) \) in the previous arguments, we will have the following:
Proposition 11
Given two sets of initial data \( (\rho _{i,0},v_{i,0}) = (\sigma _{i,0}^2 ,v_{i,0}), i = 1,2 \), for (1.2) satisfying (1.17), (1.18) and (1.19), the corresponding strong solutions \( (\rho _i,v_i) = (\sigma _i^2, v_i), i = 1,2 \), constructed in Proposition 9 in the interval \([0,T^*] \), for some \( T^* > 0 \), satisfy
In particular, if \( \rho _{1,0} = \rho _{2,0}, v_{1,0} = v_{2,0} \), we have \( \rho _1 =\rho _2, v_1 = v_2 \) in \( [0,T^*] \).
3.3 Proofs of the Main Theorems
Theorem 1 follows from Proposition 6 and Proposition 10. Theorem 2 follows from Proposition 9 and Proposition 11.
References
Brenier, Y.: Homogeneous hydrostatic flows with convex velocity profiles. Nonlinearity 12(3), 495–512, 1999
Cao, C., Ibrahim, S., Nakanishi, K., Titi, E.S.: Finite-time blowup for the inviscid primitive equations of oceanic and atmospheric dynamics. Commun. Math. Phys. 337(2), 473–482, 2015
Cao, C., Li, J., Titi, E.S.: Global well-posedness of strong solutions to the 3D primitive equations with horizontal eddy diffusivity. J. Differ. Equ. 257(11), 4108–4132, 2014
Cao, C., Titi, E.S.: Local and global well-posedness of strong solutions to the 3D primitive equations with vertical eddy diffusivity. Arch. Ration. Mech. Anal. 214(1), 35–76, 2014
Cao, C., Li, J., Titi, E.S.: Global well-posedness of the three-dimensional primitive equations with only horizontal viscosity and diffusion. Commun. Pure Appl. Math. 69(8), 1492–1531, 2016
Cao, C., Li, J., Titi, E.S.: Strong solutions to the 3D primitive equations with only horizontal dissipation: near \(H^1\) initial data. J. Funct. Anal. 272(11), 4606–4641, 2017
Cao, C., Li, J., Titi, E.S.: Global well-posedness of the 3D primitive equations with horizontal viscosity and vertical diffusivity. Phys. D: Nonlinear Phenomena 412, 132606, 2020
Cao, C., Titi, E.S.: Global well-posedness of the three-dimensional viscous primitive equations of large scale ocean and atmosphere dynamics. Ann. Math. 166(1), 245–267, 2007
Cao, C., Titi, E.S.: Global well-posedness and finite-dimensional global attractor for a 3-D planetary geostrophic viscous model. Commun. Pure Appl. Math. 56(2), 198–233, 2003
Cao, C., Titi, E.S.: Global well-posedness of the 3D primitive equations with partial vertical turbulence mixing heat diffusion. Commun. Math. Phys. 310(2), 537–568, 2012
Chen, X., Liu, J.-G.: Two nonlinear compactness theorems in \( L^p (0, T; B) \). Appl. Math. Lett. 25(12), 2252–2257, 2012
Cho, Y., Choe, H.J., Kim, H.: Unique solvability of the initial boundary value problems for compressible viscous fluids. J. Math. Pures Appl. 83(2), 243–275, 2004
Cho, Y., Kim, H.: Existence results for viscous polytropic fluids with vacuum. J. Differ. Equ. 228(2), 377–411, 2006
Cho, Y., Kim, H.: On classical solutions of the compressible Navier–Stokes equations with nonnegative initial densities. Manuscr. Math. 120(1), 91–129, 2006
Choe, H.J., Kim, H.: Strong solutions of the Navier–Stokes equations for isentropic compressible fluids. J. Differ. Equ. 190(2), 504–523, 2003
Ersoy, M., Ngom, T.: Existence of a global weak solution to compressible primitive equations. C. R. Acad. Sci. Paris Ser. I 350(7–8), 379–382, 2012
Ersoy, M., Ngom, T., Sy, M.: Compressible primitive equations: formal derivation and stability of weak solutions. Nonlinearity 24(1), 79–96, 2011
Feireisl, E.: Dynamics of Viscous Compressible Fluids. Oxford Lecture Series in Mathematics and Its Applications, vol. 26. Oxford University Press (2004).
Gatapov, B.V., Kazhikhov, A.V.: Existence of a global solution to one model problem of atmosphere dynamics. Sib. Math. J. 46(5), 805–812, 2005
Gerard-Varet, D., Masmoudi, N., Vicol, V.: Well-posedness of the hydrostatic Navier–Stokes equations. Anal. PDE 13(5), 1417–1455, 2020
Guillén-González, F.M., Masmoudi, N., Bellido, M.Á.R.: Anisotropic estimates and strong solutions of the primitive equations. Differ. Integr. Equ. 14(11), 1381–1408, 2001
Hieber, M., Kashiwabara, T.: Global strong well-posedness of the three dimensional primitive equations in \({L^p}\)-Spaces. Arch. Ration. Mech. Anal. 221(3), 1077–1115, 2016
Ignatova, M., Kukavica, I., Ziane, M.: Local existence of solutions to the free boundary value problem for the primitive equations of the ocean. J. Math. Phys. 53(10), 103101, 2012
Kobelkov, G.M.: Existence of a solution ‘in the large’ for the 3D large-scale ocean dynamics equations. C. R. Acad. Sci. Paris Ser. I 343(4), 283–286, 2006
Kukavica, I., Pei, Y., Rusin, W., Ziane, M.: Primitive equations with continuous initial data. Nonlinearity 27(6), 1135–1155, 2014
Kukavica, I., Temam, R., Vicol, V.C., Ziane, M.: Local existence and uniqueness for the hydrostatic Euler equations on a bounded domain. J. Differ. Equ. 250(3), 1719–1746, 2011
Kukavica, I., Ziane, M.: On the regularity of the primitive equations of the ocean. Nonlinearity 20(12), 2739–2753, 2007
Kukavica, I., Ziane, M.: The regularity of solutions of the primitive equations of the ocean in space dimension three. C. R. Acad. Sci. Paris Ser. I 345(5), 257–260, 2007
Li, J., Titi , E.S.: Existence and uniqueness of weak solutions to viscous primitive equations for a certain class of discontinuous initial data. SIAM J. Math. Anal. 49(1), 1–28, 2017
Liu, X., Titi, E.S.: Local well-posedness of strong solutions to the three-dimensional compressible primitive equations. arXiv:1806.09868 (2018)
Liu, X., Titi, E.S.: Global existence of weak solutions to the compressible primitive equations of atmospheric dynamics with degenerate viscosities. SIAM J. Math. Anal. 51(3), 1913–1964, 2019
Lions, J.-L., Temam, R., Wang, S.: On mathematical problems for the primitive equations of the ocean: the mesoscale midlatitude case. Nonlinear Anal. Theory Methods Appl. 40(1–8), 439–482, 2000
Lions, J.-L., Temam, R., Wang, S.: Geostrophic asymptotics of the primitive equations of the atmosphere. Topol. Methods Nonlinear Anal. 4(2), 253–287, 1994
Lions, J.-L., Temam, R., Wang, S.: On the equations of the large-scale ocean. Nonlinearity 5(5), 1007–1053, 1992
Lions , J.-L., Temam , R., Wang , S.: New formulations of the primitive equations of atmosphere and applications. Nonlinearity 5(2), 237–288, 1992
Lions, P.-L.: Mathematical Topics in Fluid Mechanics, Incompressible Models. Oxford Lecture Series in Mathematics and Its Applications, vol. 1. Oxford University Press, 1996
Masmoudi, N., Wong, T.K.: On the \(H^s\) theory of hydrostatic Euler equations. Arch. Ration. Mech. Anal. 204(1), 231–271, 2012
Petcu , M., Wirosoetisno , D.: Sobolev and Gevrey regularity results for the primitive equations in three space dimensions. Appl. Anal. 84(8), 769–788, 2005
Richardson, L.F.: Weather Prediction by Numerical Process. Cambridge University Press, 2007
Simon, J.: Compact sets in the space \( L^p(0, T;B) \). J. Annali di Matematica pura ed applicata 146(1), 65–96, 1986
Tang , T., Gao , H.: On the stability of weak solution for compressible primitive equations. Acta Appl. Math. 140(1), 133–145, 2015
Temam, R.: Navier–Stokes Equations: Theory and Numerical Analysis. Studies in Mathematics and Its Applications. Elsevier Science, 2016
Wang , F., Dou , C., Jiu , Q.: Global weak solutions to 3D compressible Primitive equations with density-dependent viscosity. J. Math. Phys. 61, 021507, 2020
Washington, W.M., Parkinson, C.L.: University Science Books. An Introduction to Three-Dimensional Climate Modeling, 2005
Wong , T.K.: Blowup of solutions of the hydrostatic Euler equations. Proc. Am. Math. Soc. 143(3), 1119–1125, 2014
Zelati , M.C., Huang , A., Kukavica , I., Temam , R., Ziane , M.: The primitive equations of the atmosphere in presence of vapour saturation. Nonlinearity 28(3), 625–668, 2015
Ziane , M., Temam , R., Changbing , H.: The primitive equations on the large scale ocean under the small depth hypothesis. Discrete Contin. Dyn. Syst. 9(1), 97–131, 2002
Acknowledgements
The authors would like to thank the referee for the useful remarks that led to the improvement of the paper. They are also thankful to the École Polytechnique for its kind hospitality, where this work was completed, and the École Polytechnique Foundation for its partial financial support through the 2017–2018 “Gaspard Monge Visiting Professor" Program. This work is supported in part by the NSF grant number DMS-1516866 and by the ONR grant N00014-15-1-2333. The work of E.S.T. was also supported in part by the Einstein Stiftung/Foundation—Berlin, through the Einstein Visiting Fellow Program.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by C. Dafermos
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Liu, X., Titi, E.S. Local Well-Posedness of Strong Solutions to the Three-Dimensional Compressible Primitive Equations. Arch Rational Mech Anal 241, 729–764 (2021). https://doi.org/10.1007/s00205-021-01662-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00205-021-01662-3