PRIMITIVE EQUATIONS WITH NON-ISOTHERMAL TURBULENT PRESSURE 5
1.2. On the physical derivations. Besides the symmetry of the relations (1.1c)-(1.1d), to mo-
tivate the presence of the non-isothermal balance (1.1d), in Section 2we provide two physical
derivations of (1.1). In both derivations the condition (1.1d) appears naturally. Following the
strategy used in the deterministic framework, we derive (1.1) by employing suitable stochastic
variants of the Boussinesq and the hydrostatic approximations. In both cases, the main ideas are
in the Boussinesq approximation. In fluid dynamics, the Boussinesq approximation is employed
in the study of buoyancy-driven flows (also referred to as natural convection), and it is typically
a good approximation in the context of oceanic flows. Roughly speaking, the idea behind the
Boussinesq approximation is that, in a natural convection regime, the role of the compressibility
is negligible in the inertial and the convection terms, but not in the gravity term. More precisely,
in the compressible Navier-Stokes equations one assumes
(1.5) pρ´ρrq`BtU´ pU¨∇qU˘«0
for some reference density ρrą0. Here, Uand ρdenote the velocity and density of the fluid,
respectively. In our first approach to derive (1.1), borrowing some ideas from stochastic climate
modeling (see e.g. [MTVE01]), we replace the right hand side in (1.5) by a noisy term:
(1.6) pρ´ρrq`BtU´ pU¨∇qU˘«ÿ
ně1“pρ´ρrqkn´∇r
Qn‰9
βn
t.
Here knPR3is given and r
Qn’s are turbulent pressures that make the modelling assumption on the
right-hand side in (1.6) compatible with the divergence-free condition which follows from assuming
ρ«ρrin the density balance, cf. (1.1e) and (2.3b).
At least formally, the right-hand side in (1.6) has zero expectation (if we interpret the noise in
the Itˆo formulation). Hence, the approximation in (1.6) is consistent with (1.5) when considering
expected values, and it can be seen as a refinement of the usual Boussinesq approximation. Em-
ploying the approximation (1.6) and the hydrostatic approximation used in the deterministic case
(see e.g. [AG01,FGH`20,LT19]) one obtains (1.1) where σn“ ´λk3
nfor some λPR, where k3
nis
the third component of knPR3. The reader is referred to Subsection 2.2 for more details.
Our second derivation of (1.1) is based on a two-scale interpretation of the primitive equations.
Indeed, as the small aspect ratio limit suggests, in the context of the primitive equations the
horizontal and the vertical directions can be thought of as small and large scales, respectively.
Hence, as usual in the literature (see e.g. [BE12,DP24,FP20,MTVE01]), it is physically reasonable
to consider an additive noise (per unit of mass) on the small-scale dynamics. Eventually, such
choice and a further variant of the Boussinesq and hydrostatic approximations lead to the system
(1.1). Details on this approach can be found in Subsection 2.3.
1.3. Comments on the literature. Here we collect further references to the literature on prim-
itive equations. Since the literature is extensive, we restrict to literature particularly relevant to
this work, referring to the references in the cited works for a more extensive and complete overview.
In the deterministic setting, the primitive equations were first studied by J. L. Lions, R. Teman,
and S. Wang in a series of articles [LTW92a,LTW92b,LTW93]. There, the authors proved the
existence of global Leray-Hopf type solutions for initial data v0PL2. As for the Navier-Stokes
equations, the uniqueness of such solutions is still open. Under additional regularity assumptions
uniqueness holds, see [Ju17]. In the deterministic setting, a breakthrough result has been proven
independently by C. Cao and E.S. Titi [CT07] and R.M. Kobelkov [Kob07] where they proved the
global well-posedness of the primitive equations via L8
tpH1
xq X L2
tpH2
xqa-priori estimates provided
v0PH1. See also [KZ07] for other boundary conditions. The results of [CT07,Kob07] have been
extended to the Lp-setting by the second author and T. Kashiwabara in [HK16]. Further results
can be found in [GGH`20a,GGH`20b,GGH`21]. See also [HH20] for an overview.
Stochastic versions of the primitive equations have been studied by several authors. Global
well-posedness for pathwise strong solutions has been established for multiplicative white noise
in time by A. Debussche, N. Glatt-Holtz and R. Temam in [DGHT11] and the same authors
with M. Ziane in [DGHTZ12]. There, the authors used a Galerkin approach to first show the
existence of martingale solutions, and then strong existence is deduced via pathwise uniqueness