Abstract
Incompressible Navier–Stokes equations on a thin spherical domain \(Q_\varepsilon \) along with free boundary conditions under a random forcing are considered. The convergence of the martingale solution of these equations to the martingale solution of the stochastic Navier–Stokes equations on a sphere \(\mathbb {S}^2\) as the thickness converges to zero is established.
Similar content being viewed by others
1 Introduction
For various motivations, partial differential equations in thin domains have been studied extensively in the last few decades; e.g. Babin and Vishik [4], Ciarlet [16], Ghidaglia and Temam [18], Marsden et al. [37] and references there in. The study of the Navier–Stokes equations (NSE) on thin domains originates in a series of papers by Hale and Raugel [20,21,22] concerning the reaction-diffusion and damped wave equations on thin domains. Raugel and Sell [44, 45] proved the global existence of strong solutions to NSE on thin domains for large initial data and forcing terms, in the case of purely periodic and periodic-Dirichlet boundary conditions. Later, by applying a contraction principle argument and carefully analysing the dependence of the solution on the first eigenvalue of the corresponding Laplace operator, Arvin [2] showed global existence of strong solutions of the Navier–Stokes equations on thin three-dimensional domains for large data. Temam and Ziane [52] generalised the results of [44, 45] to other boundary conditions. Moise et al. [41] proved global existence of strong solutions for initial data larger than in [45]. Iftimie [26] showed the existence and uniqueness of solutions for less regular initial data which was further improved by Iftimie and Raugel [27] by reducing the regularity and increasing the size of initial data and forcing.
In the context of thin spherical shells, large-scale atmospheric dynamics that play an important role in global climate models and weather prediction can be described by the 3-dimensional Navier–Stokes equations in a thin rotating spherical shell [34, 35]. Temam and Ziane in [53] gave the mathematical justification for the primitive equations of the atmosphere and the oceans which are known to be the fundamental equations of meteorology and oceanography [36, 43]. The atmosphere is a compressible fluid occupying a thin layer around the Earth and whose dynamics can be described by the 3D compressible Navier–Stokes equations in thin layers. In [53] it was assumed that the atmosphere is incompressible and hence a 3D incompressible NSE on thin spherical shells could be used as a mathematical model. They proved that the averages in the radial direction of the strong solutions (whose existence for physically relevant initial data was established in the same article) to the NSE on the thin spherical shells converge to the solution of the NSE on the sphere as the thickness converges to zero. In a recent paper Saito [46] studied the 3D Boussinesq equations in thin spherical domains and proved the convergence of the average of weak solutions of the 3D Boussinesq equations to a 2D problem. More recent work on incompressible viscous fluid flows in a thin spherical shell was carried out in [23,24,25].
For the deterministic NSE on the sphere, Il’in and Filatov [28,29,30] considered the existence and uniqueness of solutions while Temam and Wang [51] considered inertial forms of NSE on spheres. Brzeźniak et al. proved the existence and uniqueness of the solutions to the stochastic NSE on the rotating two dimensional sphere and also proved the existence of an asymptotically compact random dynamical system [9]. Recently, Brzeźniak et al. established [10] the existence of random attractors for the NSE on two dimensional sphere under random forcing irregular in space and time deducing the existence of an invariant measure.
The main objective of this article is to establish the convergence of the martingale solution of the stochastic Navier–Stokes equations (SNSE) on a thin spherical domain \(Q_\varepsilon \), whose existence can be established as in the forthcoming paper [7] to the martingale solution of the stochastic Navier–Stokes equations on a two dimensional sphere \({\mathbb {S}^2}\) [9] as thickness \(\varepsilon \) of the spherical domain converges to zero. In this way we also give another proof for the existence of a martingale solution for stochastic NSE on the unit sphere \({\mathbb {S}^2}\).
We study the stochastic Navier–Stokes equations (SNSE) for incompressible fluid
in thin spherical shells
along with free boundary conditions
In the above, \(\widetilde{u}_\varepsilon =(\widetilde{u}_\varepsilon ^r, \widetilde{u}_\varepsilon ^\lambda , \widetilde{u}_\varepsilon ^\varphi )\) is the fluid velocity field, p is the pressure, \(\nu >0\) is a (fixed) kinematic viscosity, \(\widetilde{u}_0^\varepsilon \) is a divergence free vector field on \(Q_\varepsilon \) and \(\mathbf {n}\) is the unit outer normal vector to the boundary \(\partial Q_\varepsilon \) and \(\widetilde{W}_\varepsilon (t)\), \(t \ge 0\) is an \(\mathbb {R}^N\)-valued Wiener process in some probability space \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}\right) \) to be defined precisely later.
The main result of this article is Theorem 3, which establishes the convergence of the radial averages of the martingale solution (see Definition 1) of the 3D stochastic equations (1)–(5), as the thickness of the shell \(\varepsilon \rightarrow 0\), to a martingale solution u (see Definition 2) of the following stochastic Navier–Stokes equations on the unit sphere \({\mathbb {S}^2}\subset {\mathbb R}^3\):
where \(u=(u_\lambda , u_\varphi )\) and \(\varvec{{\Delta }}^\prime \, \), \(\nabla ^\prime \) are the Laplace–de Rham operator and the surface gradient on \({\mathbb {S}^2}\) respectively. Assumptions on initial data and external forcing will be specified later.
The paper is organised as follows. We introduce necessary functional spaces in Sect. 2. In Sect. 3, we define some averaging operators and give their properties. Stochastic Navier–Stokes equations on thin spherical domains are introduced in Sect. 4 and a priori estimates for the radially averaged velocity are obtained which are later used to prove the convergence of the radial average of a martingale solution of stochastic NSE on thin spherical shell (see (1)–(5)) to a martingale solution of the stochastic NSE on the sphere (see (6)–(8)) with vanishing thickness.
2 Preliminaries
A point \(\mathbf{y}\in Q_\varepsilon \) could be represented by the Cartesian coordinates \(\mathbf{y}= \left( x, y, z\right) \) or \(\mathbf{y}= \left( r, \lambda , \varphi \right) \) in spherical coordinates, where
for \(r \in (1, 1+\varepsilon )\), \(\lambda \in [0,\pi ]\) and \(\varphi \in [0, 2 \pi )\).
For \(p \in [1, \infty )\), by \(L^p(Q_\varepsilon )\), we denote the Banach space of (equivalence-classes of) Lebesgue measurable \(\mathbb {R}\)-valued pth power integrable functions on \(Q_\varepsilon \). The \(\mathbb {R}^3\)-valued pth power integrable vector fields will be denoted by \(\mathbb {L}^p(Q_\varepsilon )\). The norm in \(\mathbb {L}^p(Q_\varepsilon )\) is given by
If \(p=2\), then \(\mathbb {L}^2(Q_\varepsilon )\) is a Hilbert space with the inner product given by
By \(\mathbb {H}^1(Q_\varepsilon ) = \mathbb {W}^{1,2}(Q_\varepsilon )\), we will denote the Sobolev space consisting of all \(u \in {\mathbb {L}^2(Q_\varepsilon )}\) for which there exist weak derivatives \(D_i u \in {\mathbb {L}^2(Q_\varepsilon )}\), \(i = 1,2,3\). It is a Hilbert space with the inner product given by
where
The Lebesgue and Sobolev spaces on the sphere \({\mathbb {S}^2}\) will be denoted by \(\mathbb {L}^p({\mathbb {S}^2})\) and \(\mathbb {W}^{s,q}({\mathbb {S}^2})\) respectively for \(p, q \ge 1\) and \(s \ge 0\). In particular, we will write \(\mathbb {H}^1({\mathbb {S}^2})\) for \(\mathbb {W}^{1,2}({\mathbb {S}^2})\).
2.1 Functional Setting on the Shell \(Q_\varepsilon \)
We will use the following classical spaces on \(Q_\varepsilon :\)
On \(\mathrm { H}_\varepsilon \), we consider the inner product and the norm inherited from \({\mathbb {L}^2(Q_\varepsilon )}\) and denote them by \(\left( \cdot , \cdot \right) _{\mathrm { H}_\varepsilon }\) and \(\Vert \cdot \Vert _{\mathrm { H}_\varepsilon }\) respectively, that is
Let us define a bilinear map \(\mathrm {a}_\varepsilon :\mathrm { V}_\varepsilon \times \mathrm { V}_\varepsilon \rightarrow \mathbb {R}\) by
where
and for \(u \in \mathrm { V}_\varepsilon \), we define
Note that for \(u \in \mathrm { V}_\varepsilon \), \(\Vert u\Vert _{\mathrm { V}_\varepsilon } = 0\) implies that u is a constant vector and \(u\cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \) i.e., u is tangent to \(Q_\varepsilon \) for every \(\mathbf{y}\in \partial Q_\varepsilon \), and thus must be 0. Hence \(\Vert \cdot \Vert _{\mathrm { V}_\varepsilon }\) is a norm on \(\mathrm { V}_\varepsilon \) (other properties can be verified easily). Under this norm \(\mathrm { V}_\varepsilon \) is a Hilbert space with the inner product given by
We denote the dual pairing between \(\mathrm { V}_\varepsilon \) and \(\mathrm { V}^\prime _\varepsilon \) by \(\langle \cdot , \cdot \rangle _\varepsilon \), that is \(\langle \cdot , \cdot \rangle _\varepsilon := {}_{\mathrm { V}_\varepsilon ^\prime }\langle \cdot , \cdot \rangle _{\mathrm { V}_\varepsilon }\). By the Lax–Milgram theorem, there exists a unique bounded linear operator \(\mathcal {A}_\varepsilon :\mathrm { V}_\varepsilon \rightarrow \mathrm { V}_\varepsilon ^\prime \) such that we have the following equality\(:\)
The operator \(\mathcal {A}_\varepsilon \) is closely related to the Stokes operator \(\mathrm { A}_\varepsilon \) defined by
The Stokes operator \(\mathrm { A}_\varepsilon \) is a non-negative self-adjoint operator in \(\mathrm { H}_\varepsilon \) (see Appendix B). Also note that
We recall the Leray–Helmholtz projection operator \(P_\varepsilon \), which is the orthogonal projector of \({\mathbb {L}^2(Q_\varepsilon )}\) onto \(\mathrm { H}_\varepsilon \). Using this, the Stokes operator \(\mathrm { A}_\varepsilon \) can be characterised as follows\(:\)
We also have the following characterisation of the Stokes operator \(\mathrm { A}_\varepsilon \) [53, Lemma 1.1]\(:\)
For \(u \in \mathrm { V}_\varepsilon \), \({\mathrm {v}}\in \mathrm {D}(\mathrm { A}_\varepsilon )\), we have the following identity (see Lemma 28)
Let \(b_\varepsilon \) be the continuous trilinear from on \(\mathrm { V}_\varepsilon \) defined by\(:\)
We denote by \(B_\varepsilon \) the bilinear mapping from \(\mathrm { V}_\varepsilon \times \mathrm { V}_\varepsilon \) to \(\mathrm { V}_\varepsilon ^\prime \) by
and we set
Let us also recall the following properties of the form \(b_\varepsilon \), which directly follows from the definition of \(b_\varepsilon \) \(:\)
In particular,
2.2 Functional Setting on the Sphere \({\mathbb {S}^2}\)
Let \(s\ge 0\). The Sobolev space \(H^s({\mathbb {S}^2})\) is the space of all scalar functions \(\psi \in L^2({\mathbb {S}^2})\) such that \((-{\varDelta ^\prime })^{s/2} \psi \in L^2({\mathbb {S}^2})\), where \({\varDelta ^\prime }\) is the Laplace–Beltrami operator on the sphere (see (177)). We similarly define \(\mathbb {H}^s({\mathbb {S}^2})\) as the space of all vector fields \(u \in \mathbb {L}^2({\mathbb {S}^2})\) such that \((-\varvec{{\Delta }}^\prime \, )^{s/2} u \in {\mathbb {L}^2({\mathbb {S}^2})}\), where \(\varvec{{\Delta }}^\prime \, \) is the Laplace–de Rham operator on the sphere (see (180)).
For \(s \ge 0\), \(\left( H^s({\mathbb {S}^2}), \Vert \cdot \Vert _{H^{s}({\mathbb {S}^2})}\right) \) and \(\left( \mathbb {H}^s({\mathbb {S}^2}), \Vert \cdot \Vert _{\mathbb {H}^{s}({\mathbb {S}^2})}\right) \) are Hilbert spaces under the respective norms, where
and
By the Hodge decomposition theorem [3, Theorem 1.72] the space of \(C^\infty \) smooth vector fields on \({\mathbb {S}^2}\) can be decomposed into three components:
where
and \(\mathcal {H}\) is the finite-dimensional space of harmonic vector fields. Since the sphere is simply connected, \(\mathcal {H}= \{0\}\). We introduce the following spaces
Note that it is known (see [50])
Given a tangential vector field u on \({\mathbb {S}^2}\), we can find vector field \(\tilde{u}\) defined on some neighbourhood of \({\mathbb {S}^2}\) such that their restriction to \({\mathbb {S}^2}\) is equal to u, that is \(\tilde{u}\vert _{{\mathbb {S}^2}} = u \in T{\mathbb {S}^2}\). Then we define
Since \(\mathbf{x}\) is orthogonal to the tangent plane \(T_{\mathbf{x}} {\mathbb {S}^2}\), \(\mathrm{curl}^\prime \; u\) is the normal component of \(\nabla \times \tilde{u}\). It could be identified with a normal vector field when needed.
We define the bilinear form \(a: \mathrm { V}\times \mathrm { V}\rightarrow \mathbb {R}\) by
The bilinear from a satisfies \(a(u,{\mathrm {v}}) \le \Vert u\Vert _{\mathbb {H}^1({\mathbb {S}^2})} \Vert {\mathrm {v}}\Vert _{\mathbb {H}^1({\mathbb {S}^2})}\) and hence is continuous on \(\mathrm { V}\). So by the Riesz representation theorem, there exists a unique operator \(\mathcal {A}: \mathrm { V}\rightarrow \mathrm { V}^\prime \) such that \(a(u, {\mathrm {v}}) = {}_{\mathrm { V}^\prime }\langle \mathcal {A}u, {\mathrm {v}}\rangle _{\mathrm { V}}\) for \(u, {\mathrm {v}}\in \mathrm { V}\). Using the Poincaré inequality, we also have \(a(u,u) \ge \alpha \Vert u\Vert ^2_{\mathrm { V}}\), for some positive constant \(\alpha \), which means a is coercive in \(\mathrm { V}\). Hence, by the Lax–Milgram theorem, the operator \(\mathcal {A}: \mathrm { V}\rightarrow \mathrm { V}^\prime \) is an isomorphism.
Next we define an operator \(\mathrm { A}\) in \(\mathrm { H}\) as follows:
By Cattabriga [15], see also Temam [49, p. 56], one can show that \(\mathrm { A}\) is a non-negative self-adjoint operator in \(\mathrm { H}\). Moreover, \(\mathrm { V}= \mathrm {D}(\mathrm { A}^{1/2})\), see [49, p. 57].
Let P be the orthogonal projection from \({\mathbb {L}^2({\mathbb {S}^2})}\) to \(\mathrm { H}\), called the Leray–Helmholtz projection. It can be shown, see [19, p. 104], that
\(\mathrm {D}(\mathrm { A})\) along with the graph norm
forms a Hilbert space with the inner product
Note that \(\mathrm {D}(\mathrm { A})\)-norm is equivalent to \(\mathbb {H}^2({\mathbb {S}^2})\)-norm. For more details about the Stokes operator on the sphere and fractional power \(\mathrm { A}^{s}\) for \(s \ge 0\), see [9, Sec. 2.2].
Given two tangential vector fields u and \({\mathrm {v}}\) on \({\mathbb {S}^2}\), we can find vector fields \(\tilde{u}\) and \(\tilde{{\mathrm {v}}}\) defined on some neighbourhood of \({\mathbb {S}^2}\) such that their restrictions to \({\mathbb {S}^2}\) are equal to, respectively, u and \({\mathrm {v}}\). Then we define the covariant derivative
where \(\pi _\mathbf{x}\) is the orthogonal projection from \(\mathbb {R}^3\) onto the tangent space \(T_\mathbf{x}{\mathbb {S}^2}\) to \({\mathbb {S}^2}\) at \(\mathbf{x}\). By decomposing \(\tilde{u}\) and \(\tilde{{\mathrm {v}}}\) into tangential and normal components and using orthogonality, one can show that
where in the last equality, we use the fact that \(\mathbf{x}\cdot {\mathrm {v}}= 0\) for any tangential vector \({\mathrm {v}}\).
We set \({\mathrm {v}}= u\) and use the formula
to obtain
Using (26) for the vector fields \(\tilde{u}\) and \(\tilde{{\mathrm {v}}} = \nabla \times \tilde{u} =\mathrm{curl}\;\tilde{u}\), we have
Thus
We consider the trilinear form b on \(\mathrm { V}\times \mathrm { V}\times \mathrm { V}\), defined by
where \(d \sigma (\mathbf {x})\) is the surface measure on \(\mathbb {S}^2\).
3 Averaging Operators and Their Properties
In this section we recall the averaging operators which were first introduced by Raugel and Sell [44, 45] for thin domains. Later, Temam and Ziane [53] adapted those averaging operators to thin spherical domains, introduced some additional operators and proved their properties using the spherical coordinate system. Recently, Saito [46] used these averaging operators to study Boussinesq equations in thin spherical domains. We closely follow [46, 53] to describe our averaging operators and provide proofs for some of the properties mentioned below.
Let \(M_\varepsilon :\mathcal {C}(Q_\varepsilon , \mathbb {R}) \rightarrow \mathcal {C}({\mathbb {S}^2}, \mathbb {R})\) be a map that projects functions defined on \(Q_\varepsilon \) to functions defined on \({\mathbb {S}^2}\) and is defined by
Remark 1
We will use the Cartesian and spherical coordinates interchangeably in this paper. For example, if \(\mathbf{x}\in {\mathbb {S}^2}\) then we will identify it by \(\mathbf{x}= (\lambda , \varphi )\) where \(\lambda \in [0, \pi ]\) and \(\varphi \in [0, 2 \pi )\).
Lemma 1
The map \(M_\varepsilon \) as defined in (28) is continuous (and linear) w.r.t norms \(L^2(Q_\varepsilon )\) and \(L^2({\mathbb {S}^2})\). Moreover,
Proof
Take \(\psi \in \mathcal {C}({Q}_\varepsilon )\) then by the definition of \(M_\varepsilon \) we have
Thus, using the Cauchy–Schwarz inequality we have
where the last equality follows from the fact that
is the volume integral over the spherical shell \(Q_\varepsilon \) in spherical coordinates, with
being the Lebesgue measure over a unit sphere. Therefore, we obtain
and hence the map is bounded and we can infer (29). \(\square \)
Corollary 1
The map \(M_\varepsilon \) as defined in (28) has a unique extension, which without the abuse of notation will be denoted by the same symbol \(M_\varepsilon :L^2(Q_\varepsilon ) \rightarrow L^2({\mathbb {S}^2})\).
Proof
Since \(\mathcal {C}(Q_\varepsilon )\) is dense in \(L^2(Q_\varepsilon )\) and \(M_\varepsilon :L^2(Q_\varepsilon ) \rightarrow L^2({\mathbb {S}^2})\) is a bounded map thus by the Riesz representation theorem there exists a unique extension. \(\square \)
Lemma 2
The following map
is bounded and
Proof
It is sufficient to consider \(\psi \in \mathcal {C}({\mathbb {S}^2})\). For \(\psi \in \mathcal {C}({\mathbb {S}^2})\), we have
where \({\mathrm {v}}(\mathbf{y}) = \frac{1}{|\mathbf{y}|}\psi \left( \frac{\mathbf{y}}{|\mathbf{y}|}\right) \). But, for \(\mathbf{x}\in {\mathbb {S}^2}\)
So
thus, showing that the map \({R}_\varepsilon \) is bounded w.r.t. \(L^2({\mathbb {S}^2})\) and \(L^2(Q_\varepsilon )\) norms. \(\square \)
Lemma 3
Let \(\psi \in W^{1,p}({\mathbb {S}^2})\) for \(p \ge 2\). Then for \(\varepsilon \in (0,1)\) there exists a constant \(C > 0\) independent of \(\varepsilon \) such that
Proof
By the definition of the map \({R}_\varepsilon \) (see (31)), and identities (172), (176) for the scalar function \(\psi \in W^{1,p}({\mathbb {S}^2})\), we have for \(Q_\varepsilon \ni \mathbf{y}= r\mathbf{x}\), \(r \in (1,1+\varepsilon )\) and \(\mathbf{x}\in {\mathbb {S}^2}\),
Hence,
\(\square \)
Lemma 4
Let \(\psi \in H^2({\mathbb {S}^2})\). Then for \(\varepsilon \in (0,1)\)
Proof
Let \(\psi \in H^2({\mathbb {S}^2})\), then
Therefore, for every \(Q_\varepsilon \ni \mathbf{y}= r\mathbf{x}\), \(r \in (1, 1+\varepsilon )\) and \(x \in {\mathbb {S}^2}\), we have (see (171) and (177) for the definition of Laplace–Beltrami operator)
Hence
Since \(\varepsilon \in (0,1)\), the inequality (32) holds. \(\square \)
Remark 2
It is easy to check that the dual operator \(M_\varepsilon ^*:L^2({\mathbb {S}^2}) \rightarrow L^2(Q_\varepsilon )\) is given by
Next we define another map
Courtesy of Corollary 1 and Lemma 2, \(\widehat{{M}}_\varepsilon \) is well-defined and bounded. Using definitions of maps \(R_\varepsilon \) and \(M_\varepsilon \) , we have
Lemma 5
Let \(\psi \in L^2(Q_\varepsilon )\), then we have the following scaling property
Proof
Let \(\psi \in L^2(Q_\varepsilon )\). Then by the defintion of the map \(\widehat{{M}}_\varepsilon \), we have
\(\square \)
The normal component of a function \(\psi \) defined on \(Q_\varepsilon \) when projected to \({\mathbb {S}^2}\) is given by the map \(\widehat{N}_\varepsilon \) which is defined by
i.e.
The following result establishes an important property of the map \(\widehat{N}_\varepsilon \).
Lemma 6
Let \(\psi \in L^2(Q_\varepsilon )\), then
Proof
Let us choose and fix \(\psi \in L^2(Q_\varepsilon )\). Then by the definitions of the operators involved we have the following equality in \(L^2({\mathbb {S}^2})\):
Therefore ,we deduce that in order to prove equality (37), it is sufficient to show that
Hence, by taking into account definitions (36) of \(\widehat{N}_\varepsilon \) and (34) of \(\widehat{{M}}_\varepsilon \) , we infer that it is sufficient to prove that
Let us choose \(\psi \in \mathcal {C}(Q_\varepsilon )\) and put \(\phi = M_\varepsilon \psi \), i.e.
Note that
Thus, we infer that
Thus, we proved \(M_\varepsilon \circ {R}_\varepsilon \circ M_\varepsilon \psi = M_\varepsilon \psi \) for every \(\psi \in \mathcal {C}({Q}_\varepsilon )\). Since \(\mathcal {C}({Q}_\varepsilon )\) is dense in \(L^2(Q_\varepsilon )\) and the maps \(M_\varepsilon \) and \(M_\varepsilon \circ {R}_\varepsilon \circ M_\varepsilon \) are bounded in \(L^2(Q_\varepsilon )\), we conclude that we have proved (37). \(\square \)
Lemma 7
For all \(\psi , \xi \in L^2(Q_\varepsilon )\), we have
Proof
Let \(\psi , \xi \in L^2(Q_\varepsilon )\), then
By the definition (34) of the map \(\widehat{{M}}_\varepsilon \) and by Lemma 6, we infer that
\(\square \)
Next we define projection operators for \(\mathbb {R}^3\)-valued vector fields using the above maps (for scalar functions), as follows \(:\)
Lemma 8
Let \(u \in {\mathbb {L}^2(Q_\varepsilon )}\). Then
Moreover, if u satisfies the boundary condition \(u \cdot \mathbf {n} = 0\), then
Proof
The normal vector \(\mathbf {n}\) to \(\partial Q_\varepsilon \) is given by \(\mathbf {n} = \left( 1, 0, 0\right) \). Thus by the definition of \(\widetilde{M}_\varepsilon \) we have
Now for the second part, from the definition of \(\widetilde{N}_\varepsilon \) we have
\(\square \)
We also have the following generalisation of Lemma 6.
Lemma 9
Let \(u \in \mathbb {L}^2(Q_\varepsilon )\), then
The following Lemma makes sense only for vector fields.
Lemma 10
Let \(u \in \mathrm { H}_\varepsilon \), then
Proof
Let \(u \in \mathrm{{U}}_\varepsilon := \left\{ {\mathrm {v}}\in \mathcal {C}({Q_\varepsilon }) :\mathrm{div}\;{\mathrm {v}}= 0\, \text{ in } \, Q_\varepsilon \, \text{ and } \, {\mathrm {v}}\cdot \mathbf {n} = 0\, \text{ on } \partial Q_\varepsilon \right\} \), then using the definition of divergence for a vector field \({\mathrm {v}}= \left( {\mathrm {v}}_r, {\mathrm {v}}_\lambda , {\mathrm {v}}_\varphi \right) \) in spherical coordinates (see (174)), we get for \(Q_\varepsilon \ni \mathbf{y}= r\mathbf{x}\), \(\mathbf{x}\in {\mathbb {S}^2}\), \(r \in (1, 1+\varepsilon )\),
Now considering each of the terms individually, we have
Using (43) and (44) in the equality (42), we obtain
Since \(u \in \mathrm { U}_\varepsilon \), \(\mathrm{div}\;u = 0\) in \(Q_\varepsilon \), which implies
Using this in (45), we get
Thus, we have proved that \(\mathrm{div}\;\widetilde{M}_\varepsilon u = 0\), for every \(u \in \mathrm { U}_\varepsilon \). Since, \(\mathrm { U}_\varepsilon \) is dense in \(\mathrm { H}_\varepsilon \), it holds true for every \(u \in \mathrm { H}_\varepsilon \) too. The second part follows from the definition of \(\widetilde{N}_\varepsilon \) and \(\mathrm { H}_\varepsilon \). \(\square \)
From Lemmas 8 and 10, we infer the following corollary\(:\)
Corollary 2
If \(u \in \mathrm { H}_\varepsilon \) then \(\widetilde{M}_\varepsilon u\) and \(\widetilde{N}_\varepsilon u\) belong to \(\mathrm { H}_\varepsilon \).
Using the definition of maps \(\widetilde{M}_\varepsilon \) and \(\widetilde{N}_\varepsilon \) and Lemma 7, we conclude:
Proposition 1
For all \(u, {\mathrm {v}}\in {\mathbb {L}^2(Q_\varepsilon )}\), we have
Moreover,
Finally we define a projection operator that projects \(\mathbb {R}^3\)-valued vector fields defined on \(Q_\varepsilon \) to the “tangent” vector fields on sphere \({\mathbb {S}^2}\).
Lemma 11
Let \(u \in {\mathbb {L}^2(Q_\varepsilon )}\), then
Proof
Let \(u \in {\mathbb {L}^2(Q_\varepsilon )}\), then
\(\square \)
Remark 3
Similar to the scalar case, one can prove that the dual operator 
is given by
Indeed, for \(u \in {\mathbb {L}^2(Q_\varepsilon )}, {\mathrm {v}}\in {\mathbb {L}^2({\mathbb {S}^2})}\)
Using the identities (168)–(170), we can show that for a divergence free smooth vector field \({\mathrm {u}}\)
We define a weighted \(L^2\)-product on \(\mathrm { H}_\varepsilon \) by
and the corresponding norm will be denoted by \(\Vert \cdot \Vert _r\) which is equivalent to \(\Vert \cdot \Vert _{{\mathbb {L}^2(Q_\varepsilon )}}\), uniformly for \(\varepsilon \in (0, \frac{1}{2})\) \(:\)
We end this section by recalling a lemma and some Poincaré type inequalities from [53].
Lemma 12
[53, Lemma 1.2] For \(u, {\mathrm {v}}\in \mathrm { V}_\varepsilon \), we have
Moreover,
Corollary 3
Let \(\varepsilon \in (0, \frac{1}{2})\) and \(u \in \mathrm { V}_\varepsilon \). Then
Proof
Let \(\varepsilon \in (0, \frac{1}{2})\) and \(u \in \mathrm { V}_\varepsilon \). Then, by relation (50), equivalence of norms (52) and Eq. (53), we have
The second inequality can be proved similarly. \(\square \)
The following two lemmas are taken from [53]. For the sake of completeness and convenience of the reader we have provided the proof in Appendix C.
Lemma 13
(Poincaré inequality in thin spherical shells) [53, Lemma 2.1] For \(0 < \varepsilon \le \tfrac{1}{2}\), we have
Lemma 14
(Ladyzhenskaya’s inequality) [53, Lemma 2.3] There exists a constant \(c_1\), independent of \(\varepsilon \), such that
Corollary 4
For \(\varepsilon \in (0, \frac{1}{2})\), there exists a constant \(c_2 > 0\) such that
Proof
Let \(u \in \mathrm { V}_\varepsilon \), then by the Hölder inequality, we have
Thus, by Lemmas 13 and 14, we get
\(\square \)
In the following lemma we enlist some properties of operators \(\widehat{{M}}_\varepsilon \), \(\widehat{N}_\varepsilon \), \(\widetilde{M}_\varepsilon \) and \(\widetilde{N}_\varepsilon \).
Lemma 15
Let \(\varepsilon > 0\). Then
-
(i)
for \(\psi \in L^2(Q_\varepsilon )\)
$$\begin{aligned}&\widehat{{M}}_\varepsilon \left( \widehat{{M}}_\varepsilon \psi \right) = \widehat{{M}}_\varepsilon \psi , \end{aligned}$$(57)$$\begin{aligned}&\widehat{N}_\varepsilon \left( \widehat{N}_\varepsilon \psi \right) = \widehat{N}_\varepsilon \psi , \end{aligned}$$(58)$$\begin{aligned}&\widehat{{M}}_\varepsilon \left( \widehat{N}_\varepsilon \psi \right) = 0, \qquad \text{ and } \qquad \widehat{N}_\varepsilon \left( \widehat{{M}}_\varepsilon \psi \right) = 0, \end{aligned}$$(59) -
(ii)
and for \(u \in {\mathbb {L}^2(Q_\varepsilon )}\)
$$\begin{aligned}&\widetilde{M}_\varepsilon \left( \widetilde{M}_\varepsilon u\right) = \widetilde{M}_\varepsilon u, \end{aligned}$$(60)$$\begin{aligned}&\widetilde{N}_\varepsilon \left( \widetilde{N}_\varepsilon u \right) = \widetilde{N}_\varepsilon u, \end{aligned}$$(61)$$\begin{aligned}&\widetilde{M}_\varepsilon \left( \widetilde{N}_\varepsilon u\right) = 0, \qquad \text{ and } \qquad \widetilde{N}_\varepsilon \left( \widetilde{M}_\varepsilon u\right) = 0. \end{aligned}$$(62)
Proof
Let \(\psi \in \mathcal {C}({Q}_\varepsilon )\). Put
i.e. for \(\mathbf{y}\in Q_\varepsilon \)
Next for \(\mathbf{z}\in Q_\varepsilon \)
Hence, we proved (57) for every \(\psi \in \mathcal {C}({Q}_\varepsilon )\). Since \(\mathcal {C}({Q}_\varepsilon )\) is dense in \(L^2(Q_\varepsilon )\), it holds true for every \(\psi \in L^2(Q_\varepsilon )\).
Proof of first part of (59). Let \(\psi \in \mathcal {C}({Q}_\varepsilon )\). Put \(\phi = \widehat{N}_\varepsilon \psi \in \mathcal {C}({Q}_\varepsilon )\). By Lemma 6
Therefore, for \(\mathbf{y}\in Q_\varepsilon \),
Therefore, we infer that
for all \(\mathbf{y}\in Q_\varepsilon \). Thus, we have established first part of (59) for all \(\psi \in \mathcal {C}({Q}_\varepsilon )\). Using the density argument, we can prove it for all \(\psi \in L^2(Q_\varepsilon )\).
Now for (58), by the definition of \(\widehat{N}_\varepsilon \) and (59), we obtain
Again using the definition of \(\widehat{N}_\varepsilon \) and Eq.(57), we have
concluding the proof of second part of (59).
Proof of (60). Let \(u \in \mathcal {C}({Q}_\varepsilon , \mathbb {R}^3)\). Write \(u = \left( u_r, u_\lambda , u_\varphi \right) \). Put \({\mathrm {v}}= \widetilde{M}_\varepsilon u\), i.e.
where
Thus, by the definition of \(\widetilde{M}_\varepsilon \) and identity (57)
We can extend this to \(u \in {\mathbb {L}^2(Q_\varepsilon )}\) by the density argument. The remaining identities can be also established similarly as in the case of scalar functions. \(\square \)
Later in the proof of Theorem 3, in order to pass to the limit we will use an operator 
defined by
where
Using the definition of map \({R}_\varepsilon \) from Lemma 2, we can rewrite 
as
Note that 
is a bounded linear map from \({\mathbb {L}^2({\mathbb {S}^2})}\) to \({\mathbb {L}^2(Q_\varepsilon )}\).
This operator 
is retract of 
, i.e. a map 
such that
One can easily show that if \(u \in \mathrm {D}(\mathrm { A})\) then 
. In particular, for \(u \in \mathrm { H}\), 
. Next we establish certain scaling properties for the map 
.
Lemma 16
Let \(\varepsilon > 0\), then
Proof
Let \(\varepsilon > 0\) and consider \({\mathbb {L}^2({\mathbb {S}^2})}\ni u = \left( 0, u_\lambda , u_\varphi \right) \). Then, by the definition of the retract operator 
and \({\mathbb {L}^2(Q_\varepsilon )}\)-norm we have
\(\square \)
Using the definition of the map 
and Lemmas 3, 4 we can deduce the following two lemmas (we provide the detailed proof of the latter in Appendix C):
Lemma 17
Let \(u \in \mathbb {W}^{1,p}({\mathbb {S}^2})\) for \(p \ge 2\). Then for \(\varepsilon \in (0,1)\) there exists a constant \(C > 0\) independent of \(\varepsilon \) such that
Lemma 18
Let \(u \in \mathbb {H}^2({\mathbb {S}^2}) \cap \mathrm { H}\) and \(\varepsilon \in (0,1)\). Then
where \(\varvec{{\Delta }}^\prime \, \) is defined in (180).
4 Stochastic NSE on Thin Spherical Domains
This section deals with the proof of our main result, Theorem 3. First we introduce our two systems; stochastic NSE in thin spherical domain and stochastic NSE on the sphere, then we present the definition of martingale solutions for both systems. We also state the assumptions under which we prove our result. In Sect. 4.1, we obtain a priori estimates (formally) which we further use to establish some tightness criterion (see Sect. 4.2) which along with Jakubowski’s generalisation of Skorokhod Theorem gives us a converging (in \(\varepsilon \)) subsequence. At the end of this section we show that the limiting object of the previously obtained converging subsequence is a martingale solution of stochastic NSE on the sphere (see Sect. 4.3).
In thin spherical domain \(Q_\varepsilon \), which was introduced in (3), we consider the following stochastic Navier–Stokes equations (SNSE)
Recall that, \(\widetilde{u}_\varepsilon =(\widetilde{u}_\varepsilon ^r, \widetilde{u}_\varepsilon ^\lambda , \widetilde{u}_\varepsilon ^\varphi )\) is the fluid velocity field, p is the pressure, \(\nu >0\) is a (fixed) kinematic viscosity, \(\widetilde{u}_0^\varepsilon \) is a divergence free vector field on \(Q_\varepsilon \) and \(\mathbf {n}\) is the unit normal outer vector to the boundary \(\partial Q_\varepsilon \). We assume thatFootnote 1\(N \in \mathbb {N}\). We consider a family of maps
such that
for some \(\widetilde{g}_\varepsilon ^j: \mathbb {R}_+\rightarrow \mathrm { H}_\varepsilon \), \(j=1, \ldots ,N\). The Hilbert–Schmidt norm of \(\widetilde{G}_\varepsilon \) is given by
Finally we assume that \(\widetilde{W}_\varepsilon (t)\), \(t \ge 0\) is an \(\mathbb {R}^N\)-valued Wiener process defined on the probability space \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}\right) \). We assume that \(\left( \beta _j\right) _{j=1}^N\) are i.i.d real valued Brownian motions such that \(W(t)=\bigl (\beta _j(t)\bigr )_{j=1}^N\), \(t\ge 0\).
In this section, we shall establish convergence of the radial averages of the martingale solution of the 3D stochastic equations (69)–(72), as the thickness of the shell \(\varepsilon \rightarrow 0\), to a martingale solution u of the following stochastic Navier–Stokes equations on the sphere \({\mathbb {S}^2}\):
where \(u=(0, u_\lambda , u_\varphi )\) and \(\varvec{{\Delta }}^\prime \, \), \(\nabla ^\prime \) are as defined in (176)–(180). Assumptions on initial data and external forcing will be specified later (see Assumptions 1, 2). Here, \(G : \mathbb {R}_+\rightarrow \mathcal {T}_2(\mathbb {R}^N; \mathrm { H})\) and W(t), \(t \ge 0\) is an \(\mathbb {R}^N\)-valued Wiener process such that
where \(N \in \mathbb {N}\), \(\left( \beta _j\right) _{j=1}^N\) are i.i.d real valued Brownian motions as before and \(\left\{ g^j\right\} _{j=1}^N\) are elements of \(\mathrm { H}\), with certain relation to \(\widetilde{g}_\varepsilon ^j\), which is specified later in Assumption 2.
Remark 4
We are aware of other formulations of the Laplacian in (75) such as the one with an additional Ricci tensor term [47, 48]. However, as it was written in [47, p. 144], “Deriving appropriate equations of motion involves dynamical considerations which do not seem adapted to Riemannian space; in particular it is not evident how to formulate the principle of conservation of momentum.” Therefore, in this paper, we follow the approach presented in [53], that the Navier–Stokes equations on the sphere is the thin shell limit of the 3-dimensional Navier–Stokes equations defined on a thin spherical shell.
Now, we specify assumptions on the initial data \(\widetilde{u}_0^\varepsilon \) and external forcing \(\widetilde{f}_\varepsilon \), \(\widetilde{g}_\varepsilon ^j\).
Assumption 1
Let \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}\right) \) be the given filtered probability space. Let us assume that \(p \ge 2\) and that \(\widetilde{u}_0^\varepsilon \in \mathrm { H}_\varepsilon \), for \(\varepsilon \in (0,1]\), such that for some \(C_1 > 0\)
We also assume that \(\widetilde{f}_\varepsilon \in L^p([0,T]; \mathrm { V}_\varepsilon ^\prime )\), for \(\varepsilon \in (0,1]\), such that for some \(C_2 > 0\),
Let \(\widetilde{W}^\varepsilon \) be an \(\mathbb {R}^N\)-valued Wiener process as before and assume that
such that, using convention (73), for each \(j=1,\ldots ,N\),
Projecting the stochastic NSE (on thin spherical shell) (69)–(72) onto \(\mathrm { H}_\varepsilon \) using the Leray–Helmholtz projection operator and using the definitions of operators from Sect. 2.1, we obtain the following abstract Itô equation in \(\mathrm { H}_\varepsilon \), \(t \in [0,T]\)
Definition 1
Let \(\varepsilon \in (0,1]\). A martingale solution to (82) is a system
where \(\left( \varOmega , \mathcal {F}, \mathbb {P}\right) \) is a probability space and \(\mathbb {F}= \left( \mathcal {F}_t\right) _{t \ge 0}\) is a filtration on it, such that
-
\(\widetilde{W}_\varepsilon \) is a \(\mathbb {R}^N\)-valued Wiener process on \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}\right) \),
-
\(\widetilde{u}_\varepsilon \) is \(\mathrm { V}_\varepsilon \)-valued progressively measurable process, \(\mathrm { H}_\varepsilon \)-valued weakly continuous \(\mathbb {F}\)-adapted process such thatFootnote 2\(\mathbb {P}\)-a.s.
$$\begin{aligned}&\widetilde{u}_\varepsilon (\cdot , \omega ) \in \mathcal {C}([0,T],\mathrm{H}_{\varepsilon }^w) \cap L^2(0,T; \mathrm { V}_\varepsilon ),\\&\mathbb {E}\left[ \frac{1}{2} \sup _{0\le s \le T} \Vert \widetilde{u}_\varepsilon (s) \Vert ^2_{{\mathbb {L}^2(Q_\varepsilon )}} + \nu \int _0^T \Vert \mathrm{curl}\;\widetilde{u}_\varepsilon (s) \Vert ^2_{{\mathbb {L}^2(Q_\varepsilon )}}\,ds \right] < \infty \end{aligned}$$and, for all \(t \in [0,T]\) and \({\mathrm {v}}\in \mathrm { V}_\varepsilon \), \(\mathbb {P}\)-a.s.,
$$\begin{aligned}&(\widetilde{u}_\varepsilon (t),{\mathrm {v}})_{{\mathbb {L}^2(Q_\varepsilon )}} + \nu \int _0^t (\mathrm{curl}\;\widetilde{u}_\varepsilon (s), \mathrm{curl}\;{\mathrm {v}})_{{\mathbb {L}^2(Q_\varepsilon )}} \,ds + \int _0^t \left\langle B_\varepsilon ( \widetilde{u}_\varepsilon (s), \widetilde{u}_\varepsilon (s)), {\mathrm {v}}\right\rangle _\varepsilon \,ds \nonumber \\&\quad = (\widetilde{u}_0^\varepsilon , {\mathrm {v}})_{{\mathbb {L}^2(Q_\varepsilon )}} + \int _0^t \left\langle \widetilde{f}_\varepsilon (s),{\mathrm {v}}\right\rangle _{\varepsilon }\,ds + \left( \int _0^t \widetilde{G}_\varepsilon (s)\,d\widetilde{W}_\varepsilon (s),{\mathrm {v}}\right) _{\mathbb {L}^2(Q_\varepsilon )}. \end{aligned}$$(83)
In the following remark we show that a martingale solution \(\widetilde{u}_\varepsilon \) of (82), as defined above, satisfies an equivalent equation in the weak form.
Remark 5
Let \(\widetilde{u}_\varepsilon = \widetilde{u}_\varepsilon (t),\) \(t \ge 0\) be a martingale solution of (82). We will use the following notations
and also from Lemma 15 we have
Then, for \(\phi \in \mathrm {D}(\mathrm { A})\), we have
and using Lemma 6, Proposition 1 and Lemma 12, we can rewrite the weak formulation identity (83) as follows.
where \(\langle \cdot , \cdot \rangle \) denotes the duality between \(\mathrm { V}^\prime \) and \(\mathrm { V}\).
Next, we present the definition of a martingale solution for stochastic NSE on \({\mathbb {S}^2}\).
Definition 2
A martingale solution to equation (75)–(77) is a system
where \(\left( \widehat{\varOmega }, \widehat{\mathcal {F}},\widehat{\mathbb {P}} \right) \) is a probability space and \(\widehat{\mathbb {F}} = \left( \widehat{\mathcal {F}}_t\right) _{t \ge 0}\) is a filtration on it, such that
-
\(\widehat{W}\) is an \(\mathbb {R}^N\)-valued Wiener process on \(\left( \widehat{\varOmega }, \widehat{\mathcal {F}}, \widehat{\mathbb {F}}, \widehat{\mathbb {P}}\right) \),
-
\(\widehat{u}\) is \(\mathrm { V}\)-valued progressively measurable process, \(\mathrm { H}\)-valued continuous \(\widehat{\mathbb {F}}\)-adapted process such that
$$\begin{aligned}&\widehat{u}(\cdot , \omega ) \in \mathcal {C}([0,T],\mathrm{H}) \cap L^2(0,T; \mathrm { V}),\\&\quad \widehat{\mathbb {E}} \left[ \sup _{0\le s \le T} \Vert \widehat{u}(s) \Vert ^2_{{\mathbb {L}^2({\mathbb {S}^2})}} + \nu \int _0^T \Vert \mathrm{{curl}^\prime }\widehat{u}(s) \Vert ^2_{{\mathbb {L}^2({\mathbb {S}^2})}}\,ds \right] < \infty \end{aligned}$$and
$$\begin{aligned}&\left( \widehat{u}(t), \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}+ \nu \int _0^t \left( \mathrm{{curl}^\prime }\widehat{u}(s), \mathrm{{curl}^\prime }\phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}\,ds + \int _0^{t} \left( \left[ \widehat{u}(s) \cdot \nabla ^\prime \right] \widehat{u}(s), \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}\,ds \nonumber \\&\quad = \left( u_0, \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}+ \int _0^{t} \left( f(s), \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}\,ds + \left( \int _0^t G(s) d \widehat{W}(s), \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}, \end{aligned}$$(86)for all \(t \in [0,T]\) and \(\phi \in \mathrm { V}\).
Assumption 2
Let \(p \ge 2\). Let \(\left( \widehat{\varOmega }, \widehat{\mathcal {F}}, \widehat{\mathbb {F}}, \widehat{\mathbb {P}}\right) \) be the given probability space, \({\mathrm {u}}_0 \in \mathrm { H}\) such that
Let \(f \in L^p([0,T]; \mathrm { V}^\prime )\), such that for every \(s \in [0,T]\),
And finally, we assume that \(G \in L^p(0,T; \mathcal {T}_2(\mathbb {R}^N; \mathrm { H}))\), such that for each \(j=1,\ldots ,N\) and \(s \in [0,T]\), 
converges weakly to \(g^j(s)\) in \(\mathbb {L}^2({\mathbb {S}^2})\) as \(\varepsilon \rightarrow 0\) and
for some \(M > 0\).
Remark 6
(Existence of martingale solutions) In a companion paper [7] we will address an easier question about the existence of a martingale solution for (1)–(5) in a more general setting with multiplicative noise. The key idea of the proof is taken from [11], where authors prove existence of a martingale solution for stochastic NSE in unbounded 3D domains.
The existence of a pathwise unique strong solution (hence a martingale solution) for the stochastic NSE on a sphere \({\mathbb {S}^2}\) is already established by two of the authors and Goldys in [9]. Through this article we give another proof of the existence of a martingale solution for such a system.
We end this subsection by stating the main theorem of this article.
Theorem 3
Let the given data \(\widetilde{u}_0^\varepsilon \), \(u_0\), \(\widetilde{f}_\varepsilon \), f, \(\widetilde{g}_\varepsilon ^j\), \(g^j\), \(j \in \{1, \ldots , N\}\) satisfy Assumptions 1 and 2. Let \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}, \widetilde{W}_\varepsilon , \widetilde{u}_\varepsilon \right) \) be a martingale solution of (69)–(72) as defined in Definition 1. Then, the averages in the radial direction of this martingale solution i.e. 
converge to a martingale solution, \(\left( \widehat{\varOmega }, \widehat{\mathcal {F}}, \widehat{\mathbb {F}}, \widehat{\mathbb {P}}, \widehat{W}, \widehat{u}\right) \), of (75)–(77) in \(L^2(\widehat{\varOmega }\times [0,T]\times {\mathbb {S}^2})\).
Remark 7
According to Remark 6, for every \(\varepsilon \in [0,1]\) there exists a martingale solution of (69)–(72) as defined in Definition 1, i.e. we will obtain a tuple \(\left( \varOmega _\varepsilon , \mathcal {F}_\varepsilon , \mathbb {F}_\varepsilon , \mathbb {P}_\varepsilon , \widetilde{W}_\varepsilon , \widetilde{u}_\varepsilon \right) \) as a martingale solution. It was shown in [31] that is enough to consider only one probability space, namely,
where \(\mathcal {L}\) denotes the Lebesgue measure on [0, 1]. Thus, it is justified to consider the probability space \(\left( \varOmega , \mathcal {F}, \mathbb {P}\right) \) independent of \(\varepsilon \) in Theorem 3.
4.1 Estimates
From this point onward we will assume that for every \(\varepsilon \in (0,1]\) there exists a martingale solution \(\left( \varOmega , \mathcal {F}, \mathbb {F}, \mathbb {P}, \widetilde{W}_\varepsilon , \widetilde{u}_\varepsilon \right) \) of (82). Please note that we do not claim neither we use the uniqueness of this solution.
The main aim of this subsection is to obtain estimates for \(\alpha _\varepsilon \) and \(\widetilde{\beta }_\varepsilon \) uniform in \(\varepsilon \) using the estimates for the process \(\widetilde{u}_\varepsilon \).
The energy inequality (90) and the higher-order estimates (105)–(106), satisfied by the process \(\widetilde{u}_\varepsilon \), as obtained in Lemmas 19 and 22 is actually a consequence (essential by-product) of the existence proof. In principle, one obtains these estimates (uniform in the approximation parameter N) for the finite-dimensional process \(\widetilde{u}_\varepsilon ^{(N)}\) (using Galerkin approximation) with the help of the Itô lemma. Then, using the lower semi-continuity of norms, convergence result (\(\widetilde{u}_\varepsilon ^{(N)} \rightarrow \widetilde{u}_\varepsilon \) in some sense), one can establish the estimates for the limiting process. Such a methodology was employed in a proof of Theorem 4.8 in the recent paper [13] by the first named author, Motyl and Ondreját.
In Lemmas 19 and 22 we present a formal proof where we assume that one can apply (ignoring the existence of Lebesgue and stochastic integrals) the Itô lemma to the infinite dimensional process \(\widetilde{u}_\varepsilon \). The idea is to showcase (though standard) the techniques involved in establishing such estimates.
Lemma 19
Let \(\widetilde{u}_0^\varepsilon \in \mathrm { H}_\varepsilon \), \(\widetilde{f}_\varepsilon \in \mathbb {L}^2(\varOmega \times [0,T]; \mathrm { V}^\prime _\varepsilon )\) and \(\widetilde{G}_\varepsilon \in \mathbb {L}^2(\varOmega \times [0,T]; \mathcal {T}_2(\mathbb {R}^N; \mathrm { H}_\varepsilon ))\). Then, the martingale solution \(\widetilde{u}_\varepsilon \) of (82) satisfies the following energy inequality
where K is some positive constant independent of \(\varepsilon \).
Proof
Using the Itô formula for the function \(\Vert \xi \Vert ^2_{{\mathbb {L}^2(Q_\varepsilon )}}\) with the process \(\widetilde{u}_\varepsilon \), for a fixed \(t \in [0,T]\) we have
Using the Cauchy–Schwarz inequality and the Young inequality, we get the following estimate
which we use in (91), to obtain
Using the Burkholder–Davis–Gundy inequality (see [32, Prop. 2.12]), we have
Taking the supremum of (92) over the interval [0, T], then taking expectation and using inequality (93) we infer the energy inequality (90). \(\square \)
Let us recall the following notations, which we introduced earlier, for \(t \in [0,T]\)
Lemma 20
Let \(\widetilde{u}_\varepsilon \) be a martingale solution of (82) and Assumption 1 hold, in particular, for \(p =2\). Then
where \(C_1,C_2\) are positive constants from (79) and (80) and \(C_3 > 0\) (determined within the proof) is another constant independent of \(\varepsilon \).
Proof
Let \(\widetilde{u}_\varepsilon \) be a martingale solution of (82), then it satisfies the energy inequality (90). From Eq. (47), we have
Moreover, by Corollary 3
Therefore, using (96) and (97) in the energy inequality (90), we get
and hence from the scaling property, Lemma 11, we have
By the assumptions on \(\widetilde{g}_\varepsilon ^j\) (81), there exists a positive constant c such that for every \(j \in \{1, \ldots , N\}\)
Therefore, using Assumption 1 and (99) in (98), cancelling \(\varepsilon \) on both sides and defining \(C_3 = NKc\), we infer inequality (95). \(\square \)
From the results of Lemma 20, we deduce that
Since \(\mathrm { V}\) can be embedded into \(\mathbb {L}^6({\mathbb {S}^2})\), by using interpolation between \(L^\infty (0,T; \mathrm { H})\) and \(L^2(0,T; \mathbb {L}^6({\mathbb {S}^2}))\) we obtain
Lemma 21
Let \(\widetilde{u}_\varepsilon \) be a martingale solution of (82) and Assumption 1 hold, in particular, for \(p =2\). Then
Proof
Let \(\widetilde{u}_\varepsilon \) be a martingale solution of (82), then it satisfies the energy inequality (90). From (47), we have
Thus, by Corollary 3
Therefore, using Assumption 1, (99), inequalities (103)–(104), in the energy inequality (90), we infer (102). \(\square \)
In the following lemma we obtain some higher order estimates (on a formal level) for the martingale solution \(\widetilde{u}_\varepsilon \), which will be used to obtain the higher order estimates for the processes \({\alpha }_\varepsilon \) and \(\widetilde{\beta }_\varepsilon \).
Lemma 22
Let Assumption 1 hold true and \(\widetilde{u}_\varepsilon \) be a martingale solution of (82). Then, for \(p > 2\) we have following estimates
and
where
and \(K_1\) is a constant from the Burkholder–Davis–Gundy inequality.
Proof
Let \(F(x) = \Vert x\Vert ^p_{{\mathbb {L}^2(Q_\varepsilon )}}\) then
and
Applying the Itô lemma with F(x) and process \(\widetilde{u}_\varepsilon \) for \(t \in [0,T]\), we have
Using the fact that \(\langle B_\varepsilon (\widetilde{u}_\varepsilon ,\widetilde{u}_\varepsilon ), \widetilde{u}_\varepsilon \rangle _\varepsilon = 0\) and \( \langle A_\varepsilon \widetilde{u}_\varepsilon ,\widetilde{u}_\varepsilon \rangle _\varepsilon = \Vert \widetilde{u}_\varepsilon \Vert _{\mathrm { V}_\varepsilon }^2 \) we arrive at
Using (107) and the Cauchy–Schwarz inequality, we get
where we recall
Using the generalised Young inequality \(abc \le a^q/q + b^r/r + c^s/s\) (where \(1/q+1/r+1/s=1\)) with \(a = \sqrt{\nu } \Vert \widetilde{u}_\varepsilon \Vert _{{\mathbb {L}^2(Q_\varepsilon )}}^{p/2-1} \Vert \widetilde{u}_\varepsilon \Vert _{\mathrm { V}_\varepsilon }\), \(b=\Vert \widetilde{u}_\varepsilon \Vert _{{\mathbb {L}^2(Q_\varepsilon )}}^{p/2-1}\), \(c = \frac{1}{\sqrt{\nu }}\Vert f_\varepsilon \Vert _{\mathrm { V}'_\varepsilon }\) and exponents \(q=2, r=p, s=2p/(p-2)\) we get
Again using the Young inequality with exponents \(p/(p-2)\), p/2 we get
Using (108) and (109) we obtain
Since \(\widetilde{u}_\varepsilon \) is a martingale solution of (82) it satisfies the energy inequality (90), hence the real-valued random variable
is a \(\mathcal {F}_t\)-martingale. Taking expectation both sides of (110) we obtain
Therefore, by the Gronwall lemma we obtain
where
By Burkholder–Davis–Gundy inequality, we have
where in the last step we have used the Young inequality with exponents \(p/(p-2)\) and p/2.
Taking supremum over \(0\le s\le t\) in (110) and using (112) we get
Thus using the Gronwall lemma, we obtain
where \(C_2\left( p, \widetilde{u}_0^\varepsilon , \widetilde{f}_\varepsilon ,\widetilde{G}_\varepsilon \right) \) and \(K_p\) are the constants as defined in the statement of lemma. We deduce (106) from (113) and (105). \(\square \)
In the following lemma we will use the estimates from previous lemma to obtain higher order estimates for \(\alpha _\varepsilon \) and \(\widetilde{\beta }_\varepsilon \).
Lemma 23
Let \(p > 2\). Let \(\widetilde{u}_\varepsilon \) be a martingale solution of (82) and Assumption 1 hold with the chosen p. Then, the processes \({\alpha }_\varepsilon \) and \(\widetilde{\beta }_\varepsilon \) (as defined in (94)) satisfy the following estimates
and
where \(K(\nu ,p)\) is a positive constant independent of \(\varepsilon \) and \(K_p\) is defined in Lemma 22.
Proof
The lemma can be proved following the steps of Lemmas 20 and 21 with the use of Proposition 1, scaling property from Lemma 11, the Cauchy–Schwarz inequality, Assumptions 1, 2 and the estimates obtained in Lemma 22. \(\square \)
4.2 Tightness
In this subsection we will prove that the family of laws induced by the processes \(\alpha _\varepsilon \) is tight on an appropriately chosen topological space \(\mathcal {Z}_T\). In order to do so we will consider the following functional spaces for fixed \(T > 0\):
\(\mathcal {C}([0,T]; {\mathrm {D}(\mathrm { A}^{-1})}) :=\) the space of continuous functions \(u: [0,T] \rightarrow {\mathrm {D}(\mathrm { A}^{-1})}\) with the topology \(\mathbf {T}_1\) induced by the norm \(\Vert u\Vert _{\mathcal {C}([0,T]; {\mathrm {D}(\mathrm { A}^{-1})})} := \sup _{t \in [0,T]}\Vert u(t)\Vert _{{\mathrm {D}(\mathrm { A}^{-1})}}\),
\(L^2_{\mathrm {w}}(0,T; \mathrm { V}) :=\) the space \(L^2(0,T; \mathrm { V})\) with the weak topology \(\mathbf {T}_2\),
\(L^2(0,T; \mathrm { H}) :=\) the space of measurable functions \(u : [0,T] \rightarrow \mathrm {H}\) such that
with the topology \(\mathbf {T}_3\) induced by the norm \(\Vert u\Vert _{L^2(0,T; \mathrm {H})}\).
Let \(\mathrm { H}_\mathrm {w}\) denote the Hilbert space \(\mathrm { H}\) endowed with the weak topology.
\(\mathcal {C}([0,T]; \mathrm {H_w}) :=\) the space of weakly continuous functions \(u: [0,T] \rightarrow \mathrm {H}\) endowed with the weakest topology \(\mathbf {T}_4\) such that for all \(h \in \mathrm {H}\) the mappings
are continuous. In particular \(u_n \rightarrow u\) in \(\mathcal {C}([0,T]; \mathrm { H}_\mathrm {w})\) iff for all \(h \in \mathrm { H}:\)
Let
and let \(\mathcal {T}\) be the supremumFootnote 3 of the corresponding topologies.
Lemma 24
The set of measures \(\left\{ \mathcal {L}(\alpha _\varepsilon ),\, \varepsilon \in (0,1]\right\} \) is tight on \(\left( \mathcal {Z}_T, \mathcal {T}\right) \).
Proof
Let \(\widetilde{u}_\varepsilon \), for some fixed \(\varepsilon >0\), be a martingale solution of problem (82). Let us choose and fix \(\phi \in \mathrm {D}(\mathrm { A})\). Then, recalling the definition (40) of the operator \(\widetilde{N}_\varepsilon \), by Lemma 9, we infer that for \(t\in [0,T]\) we have
Similarly we have, also for \(t\in [0,T]\),
Thus, by Proposition 1, identity (83), equalities (117), (118), and the notations from (94), we infer that martingale solution \(\widetilde{u}_\varepsilon \) satisfies the following equality, for \(t \in [0,T]\), \(\mathbb {P}\)-a.s.
The proof of lemma turns out to be a direct application of Corollary 6. Indeed, by Lemma 20, assumptions (a) and (b) of Corollary 6 are satisfied and therefore, it is sufficient to show that the sequence \(\left( \alpha _\varepsilon \right) _{\varepsilon > 0}\) satisfies the Aldous condition \([\mathbf{A} ]\), see Definition 6, in space \(\mathrm {D}(\mathrm { A}^{-1})\).
Let \(\theta \in (0,T)\) and \(\left( \tau _\varepsilon \right) _{\varepsilon > 0}\) be a sequence of stopping times such that \(0 \le \tau _\varepsilon \le \tau _\varepsilon +\theta \le T\). We start by estimating each term in the R.H.S. of (119). We will use the Hölder inequality, the scaling property from Lemma 11, the Poincaré type inequality (54), the Ladyzhenskaya inequality (55), inequality (56), the a priori estimates from Lemmas 20, 21, results from Lemmas 17 and 18.
In what follows, we will prove that each of the eight process from equality (119) satisfies the Aldous condition \([\mathbf{A} ]\). In order to help the reader, we will divide the following part of the proof into eight parts.
-
For the first term, we obtain
(120) -
Similarly for the second term we have
(121) -
Now we consider the first non-linear term.
(122) -
Similarly for the second non-linear term, we have
(123) -
Now as in the previous case, for the next mixed non-linear term, we obtain
(124) -
Finally, for the last non-linear term, we get
(125) -
Now for the term corresponding to the external forcing \(\widetilde{f}_\varepsilon \), we have using the radial invariance of \(M_\varepsilon \widetilde{f}_\varepsilon \) and assumption (80)
(126) -
At the very end we are left to deal with the last term corresponding to the stochastic forcing. Using the radial invariance of \(M_\varepsilon \widetilde{g}_\varepsilon ^j\), Itô isometry, scaling (see Lemma 11) and assumption (81), we get
(127)$$\begin{aligned}&\le \mathbb {E}\left( \int _{\tau _\varepsilon }^{{\tau _\varepsilon }+\theta } \sum _{j=1}^N \Vert \widetilde{g}_\varepsilon ^j(s)\Vert ^2_{\mathbb {L}^2(Q_\varepsilon )}\,ds\right) \Vert \phi \Vert _{\mathbb {L}^2({\mathbb {S}^2})}^2 \nonumber \\&\le \varepsilon Nc \Vert \phi \Vert _{\mathrm {D}(\mathrm { A})}^2 \theta := c_8 \varepsilon \cdot \theta \Vert \phi \Vert _{\mathrm {D}(\mathrm { A})}^2 . \end{aligned}$$(128)
After having proved what we had promised, we are ready to conclude the proof of Lemma 24. Since for every \(t > 0\)
one has for \(\phi \in \mathrm {D}(\mathrm { A})\),
Let us fix \(\kappa > 0\) and \(\gamma > 0\). By equality (119), the sigma additivity property of probability measure and (129), we have
Using the Chebyshev’s inequality, we get
Thus, using estimates (120)–(128) in (130), we get
Let \(\delta _i = \left( \dfrac{\kappa }{8 c_i} \gamma \right) ^2\), for \(i = 1, \ldots , 7\) and \(\delta _8 = \dfrac{\kappa ^2}{8 c_8} \gamma \). Choose \( \delta = \max _{i \in \{1, \ldots , 8\}}\delta _i\). Hence,
Since \({\alpha }_\varepsilon \) satisfies the Aldous condition \([\mathbf{A} ]\) in \({\mathrm {D}(\mathrm { A}^{-1})}\), we conclude the proof of Lemma 24 by invoking Corollary 6. \(\square \)
4.3 Proof of Theorem 3
For every \(\varepsilon > 0\), let us define the following intersection of spaces
Now, choose a countable subsequence \(\left\{ \varepsilon _k\right\} _{k \in \mathbb {N}}\) converging to 0. For this subsequence define a product space \(\mathcal {Y}_T\) given by
and \(\eta :\varOmega \rightarrow \mathcal {Y}_T\) by
Now with this \(\mathcal {Y}_T\)-valued function we define a constant \(\mathcal {Y}_T\)-sequence
Then by Lemma 24 and the definition of sequence \(\eta _k\), the set of measures \(\left\{ \mathcal {L}\left( \alpha _{\varepsilon _k}, \eta _k\right) , k \in \mathbb {N}\right\} \) is tight on \(\mathcal {Z}_T \times \mathcal {Y}_T\).
Thus, by the Jakubowski–Skorohod theoremFootnote 4 there exists a subsequence \(\left( k_n\right) _{n \in \mathbb {N}}\), a probability space \((\widehat{\varOmega }, \widehat{\mathcal {F}}, \widehat{\mathbb {P}})\) and, on this probability space, \(\mathcal {Z}_T \times \mathcal {Y}_T \times \mathcal {C}([0,T]; \mathbb {R}^N)\)-valued random variables \((\widehat{u}, \widehat{\eta }, \widehat{W})\), \(\left( \widehat{\alpha }_{\varepsilon _{k_n}}, \widehat{\eta }_{k_n}, \widehat{W}_{\varepsilon _{k_n}}\right) , n \in \mathbb {N}\) such that
and
In particular, using marginal laws, and definition of the process \(\eta _k\), we have
where \(\widehat{\beta }_{\varepsilon _{k_n}}\) is the \(k_n\)th component of \(\mathcal {Y}_T\)-valued random variable \(\widehat{\eta }_{k_n}\). We are not interested in the limiting process \(\widehat{\eta }\) and hence will not discuss it further.
Using the equivalence of law of \(\widehat{W}_{\varepsilon _{k_n}}\) and \(\widetilde{W}\) on \(\mathcal {C}([0,T];\mathbb {R}^N)\) for \(n \in \mathbb {N}\) one can show that \(\widehat{W}\) and \(\widehat{W}_{\varepsilon _{k_n}}\) are \(\mathbb {R}^N\)-valued Wiener processes (see [8, Lemma 5.2 and Proof] for details).
\(\widehat{\alpha }_{\varepsilon _{k_n}} \rightarrow \widehat{u}\) in \(\mathcal {Z}_T \), \(\widehat{\mathbb {P}}\text{-a.s. }\) precisely means that
and
Let us denote the subsequence \((\widehat{\alpha }_{\varepsilon _{k_n}}, \widehat{\beta }_{\varepsilon _{k_n}}, \widehat{W}_{\varepsilon _{k_n}})\) again by \((\widehat{\alpha }_\varepsilon , \widehat{\beta }_\varepsilon , \widehat{W}_\varepsilon )_{\varepsilon \in (0,1]}\).
Note that since \(\mathcal {B}\left( \mathcal {Z}_T \times \mathcal {Y}_T \times \mathcal {C}([0,T]; \mathbb {R}^N) \right) \subset \mathcal {B}(\mathcal {Z}_T) \times \mathcal {B}(\mathcal {Y}_T)\times \mathcal {B}\left( \mathcal {C}([0,T]; \mathbb {R}^N)\right) \), the functions \(\widehat{u}\), \(\widehat{\eta }\) are \(\mathcal {Z}_T\), \(\mathcal {Y}_T\) Borel random variables respectively.
Using the retract operator 
as defined in (63)–(65), we define new processes \(\underline{\widehat{\alpha }_\varepsilon }\) corresponding to old processes \(\widetilde{\alpha }_\varepsilon \) on the new probability space as follows
Moreover, by Lemma 16 we have the following scaling property for these new processes, i.e.
The following auxiliary result which is needed in the proof of Theorem 3, cannot be deduced directly from the Kuratowski Theorem (see Theorem 7).
Lemma 25
Let \(T > 0\) and \(\mathcal {Z}_T\) be as defined in (116). Then the following sets \(\mathcal {C}([0,T];\mathrm { H}) \cap \mathcal {Z}_T\), \(L^2(0,T; \mathrm { V}) \cap \mathcal {Z}_T\) are Borel subsets of \(\mathcal {Z}_T\).
Proof
See Appendix E.2. \(\square \)
By Lemma 25, \(\mathcal {C}([0,T]; \mathrm { H})\) is a Borel subset of \(\mathcal {Z}_T\). Since \({\alpha }_\varepsilon \in \mathcal {C}([0,T]; \mathrm { H})\), \(\widehat{\mathbb {P}}\)-a.s. and \(\widehat{\alpha }_\varepsilon \), \({\alpha }_\varepsilon \) have the same laws on \(\mathcal {Z}_T\), thus
and from estimates (95) and (114), for \(p \in [2, \infty )\)
Since \(L^2(0,T; \mathrm {V}) \cap \mathcal {Z}_T\) is a Borel subset of \(\mathcal {Z}_T\) (Lemma 25), \({\alpha }_\varepsilon \) and \(\widehat{\alpha }_\varepsilon \) have same laws on \(\mathcal {Z}_T\); from (95), we have
Since the laws of \(\eta _{k_n}\) and \(\widehat{\eta }_{k_n}\) are equal on \(\mathcal {Y}_T\), we infer that the corresponding marginal laws are also equal. In other words, the laws on \(\mathcal {B}\left( \mathcal {Y}_T^{\varepsilon _{k_n}}\right) \) of \(\mathcal {L}(\widehat{\beta }_{\varepsilon _{k_n}})\) and \(\mathcal {L}(\widetilde{\beta }_{\varepsilon _{k_n}})\) are equal for every \({k_n}\).
Therefore, from the estimates (102) and (115) we infer for \(p \in [2, \infty )\)
and
By inequality (140) we infer that the sequence \((\widehat{\alpha }_\varepsilon )_{\varepsilon > 0}\) contains a subsequence, still denoted by \((\widehat{\alpha }_\varepsilon )_{\varepsilon > 0}\) convergent weakly (along the sequence \(\varepsilon _{k_n}\)) in the space \(L^2([0,T] \times \widehat{\varOmega }; \mathrm { V})\). Since \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) in \(\mathcal {Z}_T\) \(\widehat{\mathbb {P}}\)-a.s., we conclude that \(\widehat{u}\in L^2([0,T] \times \widehat{\varOmega }; \mathrm { V})\), i.e.
Similarly by inequality (139), for every \(p \in [2, \infty )\) we can choose a subsequence of \((\widehat{\alpha }_\varepsilon )_{\varepsilon >0}\) convergent weak star (along the sequence \(\varepsilon _{k_n}\)) in the space \(L^p(\widehat{\varOmega }; L^\infty (0,T; \mathrm { H}))\) and, using (134), we infer that
Using the convergence from (134) and estimates (139)–(144) we will prove certain term-by-term convergences which will be used later to prove Theorem 3. In order to simplify the notation, in the result below we write \(\lim _{\varepsilon \rightarrow 0}\) but we mean \(\lim _{k_n \rightarrow \infty }\).
Before stating the next lemma, we introduce a new functional space \(\mathbb {U}\) as the space of compactly supported, smooth divergence free vector fields on \({\mathbb {S}^2}\):
Lemma 26
For all \(t \in [0,T]\), and \(\phi \in \mathbb {U},\) we have (along the sequence \(\varepsilon _{k_n}\))
-
(a)
\(\lim _{\varepsilon \rightarrow 0} \widehat{\mathbb {E}}\left[ \int _0^T \left| \left( \widehat{\alpha }_\varepsilon (t) - \widehat{u}(t),\phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}\right| \,dt\right] = 0\),
-
(b)
\(\lim _{\varepsilon \rightarrow 0} \;\left( \widehat{\alpha }_\varepsilon (0) - u_0,\phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}= 0\),
-
(c)
\(\lim _{\varepsilon \rightarrow 0} \widehat{\mathbb {E}}\left[ \int _0^T \left| \int _0^t \left( \dfrac{\nu }{1+\varepsilon } \mathrm{{curl}^\prime }\widehat{\alpha }_\varepsilon (s) - \mathrm{{curl}^\prime }\widehat{u}(s), \mathrm{curl}\;\phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}ds\right| \,dt\right] = 0\),
-
(d)
\(\lim _{\varepsilon \rightarrow 0} \widehat{\mathbb {E}}\left[ \int _0^T \left| \int _0^t \left( \dfrac{1}{1{+}\varepsilon } \left[ \widehat{\alpha }_\varepsilon (s) \cdot \nabla ^\prime \right] \widehat{\alpha }_\varepsilon (s) {-} \left[ \widehat{u}(s) \cdot \nabla ^\prime \right] \widehat{u}(s), \phi \right) _{\mathbb {L}^2({\mathbb {S}^2})}ds\right| \,dt\right] = 0\),
-
(e)
-
(f)
Proof
Let us fix \(\phi \in \mathbb {U}\).
(a) We know that \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) in \(\mathcal {Z}_T\). In particular,
Hence, for \(t \in [0,T]\)
Since by (139), for every \(\varepsilon > 0\), \(\sup _{t\in [0,T]}\Vert \widehat{\alpha }_\varepsilon (s)\Vert ^2_{\mathbb {L}^2({\mathbb {S}^2})}\le K_1(2)\), \(\widehat{\mathbb {P}}\)-a.s., using the dominated convergence theorem we infer that
By the Hölder inequality, (139) and (144) for every \(\varepsilon > 0\) and every \(r > 1\)
where c is some positive constant. To conclude the proof of assertion (a) it is sufficient to use (147), (148) and the Vitali’s convergence theorem.
(b) Since \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) in \(\mathcal {C}([0,T]; \mathrm { H}_w)\) \(\widehat{\mathbb {P}}\)-a.s. we infer that
Also, note that by condition (87) in Assumption 2, 
converges weakly to \(u_0\) in \(\mathbb {L}^2({\mathbb {S}^2})\).
On the other hand, by (133) we infer that the laws of \(\widehat{\alpha }_\varepsilon (0)\) and \(\alpha _\varepsilon (0)\) on \(\mathrm { H}\) are equal. Since \(\alpha _\varepsilon (0)\) is a constant random variable on the old probability space, we infer that \(\widehat{\alpha }_\varepsilon (0)\) is also a constant random variable (on the new probability space) and hence, by (72) and (94), we infer that 
almost surely (on the new probability space). Therefore we infer that
concluding the proof of assertion (b).
(c) Since \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) in \(\mathcal {C}([0,T]; \mathrm { H}_w)\) \(\widehat{\mathbb {P}}\)-a.s.,
The Cauchy–Schwarz inequality and estimate (140) infer that for all \(t \in [0,T]\) and \(\varepsilon \in (0,1]\)
for some constant \(c > 0\). By (150), (151) and the Vitali’s convergence theorem we conclude that for all \(t \in [0,T]\)
Assertion (c) follows now from (140), (143) and the dominated convergence theorem.
(d) For the non-linear term, using the Sobolev embedding \(\mathbb {H}^2({\mathbb {S}^2}) \hookrightarrow \mathbb {L}^\infty ({\mathbb {S}^2})\), we have
The first term converges to zero as \(\varepsilon \rightarrow 0\), since \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) strongly in \(L^2(0,T; \mathrm { H})\) \(\widehat{\mathbb {P}}\)-a.s., \(\widehat{u} \in L^2(0,T; \mathrm { V})\) and the second term converges to zero too as \(\varepsilon \rightarrow 0\) because \(\widehat{\alpha }_\varepsilon \rightarrow \widehat{u}\) weakly in \(L^2(0,T; \mathrm { V})\). Using the Hölder inequality, estimates (139) and the embedding \(\mathbb {H}^2({\mathbb {S}^2}) \hookrightarrow \mathbb {L}^\infty ({\mathbb {S}^2})\) we infer that for all \(t \in [0,T]\), \(\varepsilon \in (0,1]\), the following inequalities hold
By (152), (153) and the Vitali’s convergence theorem we obtain for all \(t \in [0,T]\),
Using the Hölder inequality and estimates (139), (144), we obtain for all \(t \in [0,T]\), \(\varepsilon \in (0,1]\)
where \(c, \widetilde{c} > 0\) are constants. Hence by (154) and the dominated convergence theorem, we infer assertion (d).
(e) Assertion (e) follows because by Assumption 2 the sequence 
converges weakly in \(L^2(0,T; {\mathbb {L}^2({\mathbb {S}^2})})\) to f.
(f) By the definition of maps \(\widetilde{G}_\varepsilon \) and G, we have
Since, by Assumption 2, for every \(j \in \{1, \ldots , N\}\), and \(s \in [0,t]\), 
converges weakly to \(g^j(s)\) in \({\mathbb {L}^2({\mathbb {S}^2})}\) as \(\varepsilon \rightarrow 0\), we get
By assumptions on \(\widetilde{g}_\varepsilon ^j\), we obtain the following inequalities for every \(t \in [0,T]\) and \(\varepsilon \in (0,1]\)
where \(c, \widetilde{c} > 0\) are some constants. Using the Vitali’s convergence theorem, by (155) and (156) we infer
Hence, by the properties of the Itô integral we deduce that for all \(t \in [0,T]\),
By the Itô isometry and assumptions on \(\widetilde{g}_\varepsilon ^j\) and \(g^j\) we have for all \(t \in [0,T]\) and \(\varepsilon \in (0,1]\)
where \(c > 0\) is a constant. Thus, by (158), (159) and the dominated convergence theorem assertion (f) holds. \(\square \)
Lemma 27
For all \(t \in [0,T]\) and \(\phi \in \mathbb {U}\), we have (along the sequence \(\varepsilon _{k_n}\))
-
(a)
-
(b)
-
(c)
-
(d)
where the process \(\underline{\widehat{\alpha }_\varepsilon }\) is defined in (136).
Proof
Let us fix \(\phi \in \mathbb {U}\).
(a) Let \(t \in [0,T]\), then by the Hölder inequality, Lemma 18, Poincaré inequality and estimate (142), we have the following inequalities
Thus
We infer assertion (a) by dominated convergence theorem, estimate (142) and convergence (160).
(b) Using the Hölder inequality, scaling property (Lemma 11), Corollary 4, relations (67), (68), and estimates (139), (142), we get for \(t \in [0,T]\)
Thus
We infer assertion (b) by dominated convergence theorem, estimates (139), (142) and convergence (161). Assertion (c) can be proved similarly.
(d) Now for the last one, using the Hölder inequality, Corollary 4, relations (67), (68), and estimates (141), (142), we get for \(t \in [0,T]\)
Thus
We infer assertion (d) by dominated convergence theorem, estimates (141), (142) and convergence (162). \(\square \)
Finally, to finish the proof of Theorem 3, we will follow the methodology as in [42] and introduce some auxiliary notations (along sequence \(\varepsilon _{k_n}\))
Corollary 5
Let \(\phi \in \mathbb {U}\). Then (along the sequence \(\varepsilon _{k_n}\))
and
Proof
Assertion (165) follows from the equality
and Lemma 26 (a). To prove assertion (166), note that by the Fubini Theorem, we have
To conclude the proof of the corollary, it is sufficient to note that by Lemma 26\((b)-(f)\) and Lemma 27, each term on the right hand side of (163) tends at least in \(L^1(\widehat{\varOmega } \times [0,T])\) to the corresponding term (to zero in certain cases) in (164). \(\square \)
Conclusion of proof of Theorem 3
Let us fix \(\phi \in \mathbb {U}\). Since \({\alpha }_\varepsilon \) is a solution of (85), for all \(t \in [0,T]\),
In particular,
Since \(\mathcal {L}({\alpha }_\varepsilon , \eta _{k_n}, \widetilde{W}_\varepsilon ) = \mathcal {L}(\widehat{\alpha }_\varepsilon , \widehat{\eta }_{k_n}, \widehat{W}_\varepsilon )\), on \(\mathcal {B}\left( \mathcal {Z}_T \times \mathcal {Y}_T \times \mathcal {C}([0,T]; \mathbb {R}^N)\right) \) (along the sequence \(\varepsilon _{k_n}\)),
Therefore by Corollary 5 and the definition of \(\varLambda \), for almost all \(t \in [0,T]\) and \(\widehat{\mathbb {P}}\)-almost all \(\omega \in \widehat{\varOmega }\)
i.e. for almost all \(t \in [0,T]\) and \(\widehat{\mathbb {P}}\)-almost all \(\omega \in \widehat{\varOmega }\)
Hence (167) holds for every \(\phi \in \mathbb {U}\). Since \(\widehat{u}\) is a.s. \(\mathrm { H}\)-valued continuous process, by a standard density argument, we infer that (167) holds for every \(\phi \in \mathrm { V}\) (\(\mathbb {U}\) is dense in \(\mathrm { V}\)).
Putting \(\widehat{\mathcal {U}} := \left( \widehat{\varOmega }, \widehat{\mathcal {F}}, \widehat{\mathbb {F}}, \widehat{\mathbb {P}}\right) \), we infer that the system \(\left( \widehat{\mathcal {U}}, \widehat{W}, \widehat{u}\right) \) is a martingale solution to (75)–(77). \(\square \)
Notes
We could have considered the case \(N = \infty \). This case will be considered in the companion paper [7].
The space \(X^w\) denotes a topological space X with weak topology. In particular, \(\mathcal {C}([0,T]; X^w)\) is the space of weakly continuous functions \({\mathrm {v}}:[0,T] \rightarrow X\).
\(\mathcal {T}\) is the supremum of topologies \(\mathbf {T}_1\), \(\mathbf {T}_2\), \(\mathbf {T}_3\) and \(\mathbf {T}_4\), i.e. it is the coarsest topology on \(\mathcal {Z}_T\) that is finer than each of \(\mathbf {T}_1\), \(\mathbf {T}_2\), \(\mathbf {T}_3\) and \(\mathbf {T}_4\).
The space \(\mathcal {Z}_T \times \mathcal {Y}_T \times \mathcal {C}([0,T]; \mathbb {R}^N)\) satisfies the assumption of Theorem 6. Indeed, since \(Z_T\) and \(\mathcal {Y}_T^\varepsilon \), \(\varepsilon > 0\) satisfies the assumptions (see [5, Lemma 4.10]) and \(\mathcal {C}([0,T]; \mathbb {R}^N)\) is a Polish space and thus automatically satisfying the required assumptions.
References
Aldous, D.: Stopping times and tightness. Ann. Probab. 6(2), 335–340 (1978)
Avrin, J.D.: Large-eigenvalue global existence and regularity results for the Navier-Stokes equation. J. Differ. Equ. 127, 365–390 (1996)
Aubin, T.: Some Nonlinear Problems in Riemannian Geometry. Springer, New York (1998)
Babin, A.V., Vishik, M.I.: Attractors of partial differential equations and estimate of their dimension. Russ. Math. Surv. 38, 151–213 (1983)
Brzeźniak, Z., Dhariwal, G.: Stochastic constrained Navier-Stokes equations on \(\mathbb{T}^{2}\). Submitted (2019). arXiv:1701.01385
Brzeźniak, Z., Dhariwal, G.: Stochastic tamed Navier-Stokes equations on \(\mathbb{R}^{3}\): the existence and the uniqueness of solutions and the existence of an invariant measure. To appear in J. Math. Fluid Mech. (2020). arXiv:1904.13295
Brzeźniak, Z., Dhariwal, G., Le Gia, Q.T.: Stochastic Navier-Stokes equations on a thin spherical domain: Existence of a martingale solution (In preparation)
Brzeźniak, Z., Goldys, B., Jegaraj, T.: Weak solutions of a stochastic Landau-Lifshitz-Gilbert equation. Appl. Math. Res. eXpress 2013(1), 1–33 (2013)
Brzeźniak, Z., Goldys, B., Le Gia, Q.T.: Random dynamical systems generated by stochastic Navier-Stokes equations on a rotating sphere. J. Math. Anal. Appl. 426, 505–545 (2015)
Brzeźniak, Z., Goldys, B., Le Gia, Q.T.: Random attractors for the stochastic Navier-Stokes equations on the 2D unit sphere. J. Math. Fluid. Mech. 20, 227–253 (2018)
Brzeźniak, Z., Motyl, E.: Existence of a martingale solution of the stochastic Navier-Stokes equations in unbounded 2D and 3D domains. J. Differ. Equ. 254(4), 1627–1685 (2013)
Brzeźniak, Z., Motyl, E.: The existence of martingale solutions to the stochastic Boussinesq equations. Glob. Stoch. Anal. 1(2), 175–216 (2014)
Brzeźniak, Z., Motyl, E., Ondreját, M.: Invariant measure for the stochastic Navier-Stokes equations in unbounded 2D domains. Ann. Probab. 45(5), 3145–3201 (2017)
Brzeźniak, Z., Ondreját, M.: Stochastic wave equations with values in Riemanninan manifolds. Stochastic partial differential equations and applications. Quaderni di Matematica 25, 65–97 (2011)
Cattabriga, L.: Su un problema al contorno relativo al sistema di equazioni di Stokes. Rend. Semin. Mat. Univ. Padova 31, 308–340 (1961)
Ciarlet, P.G.: Plates and Junctions in Elastic Multi-structures. An Asymptotic Analysis. Masson, Paris (1990)
Duvaut, G., Lions, J.L.: Inequalities in Mechanics and Physics. Springer, Berlin (1976)
Ghidaglia, J.M., Temam, R.: Lower bound on the dimension of the attractor for the Navier-Stokes equations in space dimension 3. In: Mechanics, Analysis and Geometry: 200 Years After Lagrange, pp. 33–60, North-Hollan Delta Ser., North-Holland, Amsterdam (1991)
Grigoryan, A.: Heat Kernel and Analysis on Manifolds. Amer. Math. Soc, Providence, RI (2009)
Hale, J.K., Raugel, G.: A damped hyperbolic equation on thin domains. Trans. Am. Math. Soc. 329, 185–219 (1992)
Hale, J.K., Raugel, G.: Partial differential equations on thin domains. In: Differential Equations and Mathematical Physics (Birmingham, AL, 1990), pp. 63–97, Math. Sci. Engrg., vol. 186. Academic Press, Boston (1992)
Hale, J.K., Raugel, G.: Reaction-diffusion equation on thin domains. J. Math. Pures Appl. 71, 33–95 (1992)
Ibragimov, R.N., Pelinovsky, D.E.: Incompressible viscous fluid flows in a thin spherical shell. J. Math. Fluid Mech. 11, 60–90 (2009)
Ibragimov, R.N.: Nonlinear viscous fluid patterns in a thin rotating spherical domain and applications. Phys. Fluids 23, 123102 (2011)
Ibragimov, N.H., Ibragimov, R.N.: Integration by quadratures of the nonlinear Euler equations modeling atmospheric flows in a thin rotating spherical shell. Phys. Lett. A 375, 3858 (2011)
Iftimie, D.: The 3D Navier-Stokes equations seen as a perturbation of the 2D Navier-Stokes equations. Bull. Soc. Math. France 127, 473–517 (1999)
Iftimie, D., Raugel, G.: Some results on the Navier-Stokes equations in thin 3D domains. J. Differ. Equ. 169, 281–331 (2001)
Il’in, A.A.: The Navier-Stokes and Euler equations on two dimensional manifolds. Math. USSR Sb. 69, 559–579 (1991)
Il’in, A.A.: Partially dissipative semigroups generated by the Navier-Stokes system on two dimensional manifolds and their attractors. Russ. Acad. Sci. Sb. Math. 78, 47–76 (1994)
Il’in, A.A., Filatov, A.N.: On unique solvability of the Navier-Stokes equations on the two dimensional sphere. Sov. Math. Dokl. 38, 9–13 (1989)
Jakubowski, A.: The almost sure Skorokhod representation for subsequences in nonmetric spaces. Teor. Veroyatn. Primen. 42(1), 209–216 (1998); translation in Theory Probab. Appl. 42(1), 167–174 (1998)
Kruse, R.: Strong and Weak Approximation of Semilinear Stochastic Evolution Equations. Lecture Notes Math., vol. 2093. Springer, Cham (2014)
Kuratowski, K.: Topologie, Vol. I (French)3’eme Ed. Monografie Matematyczne, Tom XX, Polskie Towarzystwo Matematyczne, Warszawa (1952)
Lions, J.L., Temam, R., Wang, S.: New formulations of the primitive equations of atmosphere and applications. Nonlinearity 5, 237–288 (1992)
Lions, J.L., Temam, R., Wang, S.: On the equations of the large-scale ocean. Nonlinearity 5, 1007–1053 (1992)
Lions, J.L., Temam, R., Wang, S.H.: Mathematical theory for the coupled atmosphere-ocean models. J. Math. Pures Appl. (9) 74(2), 105–163 (1995)
Marsden, J.E., Raitu, T.S., Raugel, G.: Les équation d’Euler dans des coques sphériques minces. C. R. Acad. Sci. Paris 321, 1201–1206 (1995)
Métivier, M.: Semimartingales: A Course on Stochastic Processes. Walter de Gruyter & and Co., Berlin (1982)
Métivier, M.: Stochastic Partial Differential Equations in Infinite Dimensions, p. 142. Scuola Normale Superiore, Pisa (1988)
Mikulevicius, R., Rozovskii, B.L.: Global \(L^2\)-solutions of stochastic Navier-Stokes equations. Ann. Prob. 33(1), 137–176 (2005)
Moise, I., Temam, R., Ziane, M.: Asymptotic analysis of the Navier-Stokes equations in thin domains. Topol. Methods Nonlinear Anal. 10, 249–282 (1997)
Motyl, E.: Stochastic hydrodynamic-type evolution equations driven by Lévy noise in 3D unbounded domains—abstract framework and applications. Stoch. Process. Appl. 124, 2052–2097 (2014)
Pedlosky, J.: Geophysical Fluid Dynamics, 2nd edn, p. xiv+710. Springer, New York (1987)
Raugel, G., Sell, G.R.: Navier-Stokes equations on thin 3D domains I. Global attractors and global regularity of solutions. J. Am. Math. Soc. 6, 503–568 (1993)
Raugel, G., Sell, G.R.: Navier-Stokes equations on thin 3D domains II. Global regularity of spatially periodic solutions. In: Nonlinear Partial Differential Equations and Their Applications, Collège de France Seminar, vol. XI, pp. 205–247, Longman, Harlow (1994)
Saito, J.: Boussinesq equations in thin spherical domains. Kyushu J. Math. 59, 443–465 (2005)
Serrin, J.: Mathematical Principles of Classical Fluid Mechanics, Encly. of Physics, vol. 8, pp. 125–263. Springer, New York (1959)
Taylor, M.E.: Analysis on Morrey spaces and applications to Navier-Stokes and other evolution equations. Commun. Partial Differ. Equ. 17, 1407–1456 (1992)
Temam, R.: Infinite-Dimensional Dynamical Systems in Mechanics and Physics, 2nd edn, p. xxi+650. Springer, New York (1997)
Temam, R.: Navier-Stokes Equations: Theory and Numerical Analysis. American Mathematical Society, UK (2000)
Temam, R., Wang, S.: Inertial forms of Navier-Stokes equations on the sphere. J. Funct. Anal. 117, 215–242 (1993)
Temam, R., Ziane, M.: Navier-Stokes equations in three-dimensional thin domains with various boundary conditions. Adv. Differ. Equ. 1, 499–546 (1996)
Temam, R., Ziane, M.: Navier-Stokes equations in thin spherical domains. Contemp. Math. 209, 281–314 (1997)
Acknowledgements
Open access funding provided by Austrian Science Fund (FWF).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The research of all three authors is partially supported by Australian Research Council Discover Project Grant DP180100506, “Uncertainty on Spheres and Shells: Mathematics and Methods for Applications”. Zdzisław Brzeźniak has been supported by the Leverhulme Project Grant Ref No RPG-2012-514 and by Australian Research Council Discover Project Grant DP160101755. The research of Gaurav Dhariwal was supported by Department of Mathematics, University of York and is partially supported by the Austrian Science Fund (FWF) Grants P30000, W1245, and F65.
Appendices
Vector Analysis in Spherical Coordinates
In this appendix we collect some basic results from vector algebra and formulas for Laplace and gradient of scalar function and vector fields in spherical coordinates.
The following identities are very well known [53, Appendix] in vector algebra\(:\) Let \({\mathrm {u}}, {\mathrm {v}}\) and \({\mathrm {w}}\) be \(\mathbb {R}^3\)-valued smooth vector fields then
The Laplace–Beltrami operator of a scalar function \(\psi \), in spherical coordinates \((r, \lambda , \varphi )\) is given by
and its gradient is given by
For a vector field u written in the spherical coordinates, \(u = u_r \widehat{e_r} + u_\lambda \widehat{e_\lambda } + u_\varphi \widehat{e_\varphi }\), the curl and the divergence are given as follows
The Laplacian of a vector field in spherical coordinates is
where \(\nabla ^2 u_r\), \(\nabla ^2 u_\lambda \) and \(\nabla ^2 u_\varphi \) are as in (171).
We recall some standard differential operators on the unit sphere \({\mathbb {S}^2}\). For \(\psi \) a scalar function defined on \({\mathbb {S}^2}\), the tangential gradient is given by
The Laplace–Beltrami of a scalar function \(\psi \) is
For a tangential vector field \({\mathrm {v}}\) defined on \({\mathbb {S}^2}\), \({\mathrm {v}}= {\mathrm {v}}_\lambda \hat{e}_\lambda + {\mathrm {v}}_\varphi \hat{e}_\varphi \), the tangential divergence is expressed by
and the \(\mathrm{{curl}}^\prime {\mathrm {v}}\) is the scalar function defined by
The Laplace–de Rham operator applied to a vector field \({\mathrm {v}}\) is given by
where \({\varDelta ^\prime }{\mathrm {v}}_\varphi \) and \({\varDelta ^\prime }{\mathrm {v}}_\lambda \) are as in (177).
The curl and the Stokes operator
In this section we present a integration by parts formula corresponding to \(\mathrm {curl}\) operator and later we use it to give a relation between the Stokes operator \(\mathrm { A}_\varepsilon \) and \(\mathrm{curl}\;\).
Let \(\mathcal {O} \subset \mathbb {R}^3\) be a bounded domain with a regular boundary \(\partial \mathcal {O}\). Define
\(\mathrm { H}(\mathrm {curl})\) with the graph norm
is a Hilbert space.
The following theorem is a reformulation of Lemma 4.2 from [17, p. 341].
Theorem 4
Assume that \(\mathcal {O} \subset \mathbb {R}^3\) be a bounded domain with a regular boundary and \(\mathbf {n}\) be unit normal vector field on \(\varGamma = \partial \mathcal {O}\) (directed towards exterior of \(\mathcal {O}\)). Then there exists a unique bounded linear map
such that
if \(u \in C^1_0(\overline{\mathcal {O}})\) and
for every \(u \in \mathrm { H}(\mathrm {curl})\) and \({\mathrm {v}}\in \mathbb {H}^1(\mathcal {O})\).
Remark 8
We will call formula (184) the generalised Stokes formula. From now on we will write \(n \times u \big |_\varGamma \) instead of \((n \times \cdot )(u)\big |_\varGamma \) for \(u \in \mathrm { H}(\mathrm {curl})\).
Recall that
and the Stokes operator is given by
Remark 9
To define \(n \times \mathrm{curl}\;u\) as an element of \(\mathbb {H}^{-1/2}(\varGamma )\), we need to know that \(\mathrm{curl}\;(\mathrm{curl}\;u) \in \mathbb {L}^2(\mathcal {O})\). But if \(u \in \mathbb {H}^2(\mathcal {O})\), then obviously this condition is satisfied.
Theorem 5
A is self-adjoint and non-negative on \(\mathrm { H}\).
Proof
Here we will only show that A is symmetric and non-negative. Let \(u, {\mathrm {v}}\in D(A)\), then
Recall that (from (168)) for smooth \(\mathbb {R}^3\)-valued vector fields,
Using the above identity in (188) along with the fact that \(u \in D(A)\), in particular, \(\mathrm{div}\;u = 0\) and generalised Stokes formula (184), we have
Similarly
This establishes that A is symmetric on \(\mathrm { H}\). The non-negativity follows from the above identity by taking \({\mathrm {v}}= u \in D(A)\). \(\square \)
Using the definition of D(A), we can characterise \(D(A^{1/2})\) as
By Theorem 6.1 [17, Pg 358] we have
We use the following relation repeatedly in our calculations.
Lemma 28
Let \(u \in D(A^{1/2})\) and \({\mathrm {v}}\in D(A)\). Then
Proof
Note that for \(u \in D(A^{1/2})\) and \({\mathrm {v}}\in D(A)\), the LHS makes sense. Using the generalised Stokes formula (184), we get
To finish the proof we need the following lemma: \(\square \)
Lemma 29
Let \(u \in \mathbb {H}^1(\mathcal {O})\) such that \(u \cdot \mathbf {n} = 0\) on \(\varGamma \) and \(\varPhi \in \mathbb {H}^1(\mathcal {O})\) with \(\mathbf {n}\times \varPhi = 0\) on \(\varGamma \). Then
Proof
It is sufficient to prove (189) for \(u \in C^1_0(\overline{\mathcal {O}})\) with \(u \cdot \mathbf {n} =0\) on \(\varGamma \). In this case for all \(x \in \varGamma \)
\(\square \)
The proof of Lemma 28 is finished by observing that
from the proof of Theorem 5. \(\square \)
Let us consider an abstract framework. Let H be a Hilbert space and A be a non-negative self-adjoint operator on H.
Lemma 30
Let \(u \in D(A^{1/2})\) and \({\mathrm {v}}\in D(A)\). Then
Proof
Take \(u \in D(A^{1/2})\), \({\mathrm {v}}\in D(A)\). Then \(A^{1/2} {\mathrm {v}}\in D(A^{1/2})\). So by self-adjointness of \(A^{1/2}\) and \(A^{1/2}A^{1/2} = A\),
\(\square \)
Lemma 31
Let \(u, {\mathrm {v}}\in D(A^{1/2})\), then by Lemmas 30 and 28,
Proof
Let \(u, {\mathrm {v}}\in D(A)\), then
By density argument, this is true for all \(u, {\mathrm {v}}\in D(A^{1/2})\). \(\square \)
As a consequence of the above lemma we have
Proof of the Poincaré and the Ladyzhenskaya Inequalities
Proof of Lemma 13
We will establish the Poincaré inequality (54) following the footsteps of Lemma 2.1 [53] with all the details. By density argument, it is enough to prove (54) for smooth functions. Let \(\psi \in \mathcal {C}(\overline{Q_\varepsilon })\) be a real continuous function. We write for any \(\xi , \eta \in [1, 1 + \varepsilon ]\):
with \(\mathbf{x}= \frac{\mathbf{y}}{|\mathbf{y}|} \in {\mathbb {S}^2}\). We fix \(\xi \) and integrate w.r.t. \(\eta \in [1, 1+\varepsilon ]\) to obtain
With \(\psi = u_r\) and \(\xi = 1\), observing that \(u_r(1, \mathbf{x}) = 0\) (because of the boundary condition \(u \cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \)) from (191) we obtain
Applying (191) with \(\psi = \widehat{N}_\varepsilon u_\lambda \), we get
Observing from Lemma 6 for every \(\psi \in L^2(Q_\varepsilon )\)
and since the first term on the LHS of (193) is positive we can simplify (193) as follows
Similarly for \(\psi = \widehat{N}_\varepsilon u_\varphi \), we have
Thus, using (192), (194) and (195) for each of the cases \(\psi = u_r\), \(\psi = \widehat{N}_\varepsilon u_\lambda \) and \(\psi = u_\varphi \), we obtain
Using the Cauchy–Schwarz inequality, we find
the last inequality follows since \(r \ge 1\). On rearranging, we obtain
which implies for \(0 \le \varepsilon < \frac{1}{2}\)
We then integrate w.r.t. \(\lambda \) and \(\varphi \) to obtain
Adding (200) for \(\psi = u_r\), \(\psi = \widehat{N}_\varepsilon u_\lambda \) and \(\psi = \widehat{N}_\varepsilon u_\varphi \); using finally
we conclude the proof of the inequality (54). \(\square \)
Proof of Lemma 14
We will prove the lemma for smooth vector fields \(u \in \mathcal {C}^\infty (Q_\varepsilon )\). By Lemma 2.3 [53] there exists a constant \(c_0 > 0\) s.t.
By Lemma 6.1 [17, Eq. 6.11, Pg 359] for vector fields \({\mathrm {v}}\in \mathcal {C}^1(Q_\varepsilon )\) with \({\mathrm {v}}\cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \), we have a constant \(c_2 > 0\) s.t.
Also by Poincaré inequality (54) for all \(u \in \mathrm { V}_\varepsilon \), we have
Using (202) with \({\mathrm {v}}= \widetilde{N}_\varepsilon u\) along with the fact that \(\mathrm{div}\;{\mathrm {v}}= 0\) if \(\mathrm{div}\;u = 0\) and \({\mathrm {v}}\cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \) if \(u \cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \) and combining it with the Poincaré inequality (203), we obtain
Therefore, choosing \(c_1 = \sqrt{2} c_0 \left( 1 + 4 c_2 \right) ^{1/2}\) we establish (55) for every \(u \in \mathcal {C}^1(Q_\varepsilon )\) with \(\mathrm{div}\;u = 0\) on \(Q_\varepsilon \) and \(u \cdot \mathbf {n} = 0\) on \(\partial Q_\varepsilon \). We finish the proof using density argument. \(\square \)
Proof of Lemma 18
Let u be a tangential vector field defined on \({\mathbb {S}^2}\), \(u = (0, u_\lambda , u_\varphi )\). Then using the definition of the map 
(see (63)), for \(Q_\varepsilon \ni \mathbf{y}= r \mathbf{x}\), \(r \in (1,1+\varepsilon )\) and \(\mathbf{x}\in {\mathbb {S}^2}\), we have
Using the definitions of Laplace (\(\varDelta \)) for vector fields in spherical coordinates, Laplace–Beltrami (\(\nabla ^2\)) for scalars, tangential Laplace (\(\varvec{{\Delta }}^\prime \, \)) for tangential vector fields and Laplace–Beltrami (\(\varDelta ^\prime \)) for the scalar defined on \({\mathbb {S}^2}\), we have following relations:
where
If \(\mathrm{{div}^\prime }u = 0\), then by the definition of \(\mathrm{{div}^\prime }\),
which is equivalent to
Hence using all the above relations, we obtain
Since \(\varepsilon \in (0,1)\), the inequality (68) holds. \(\square \)
Compactness
1.1 Skorohod Theorem and Aldous Condition
Let E be a separable Banach space with the norm \(\Vert \cdot \Vert _E\) and let \(\mathcal {B}(E)\) be its Borel \(\sigma \)-field. The family of probability measures on \((E, \mathcal {B}(E))\) will be denoted by \(\mathcal {P}\). The set of all bounded and continuous E-valued functions is denoted by \(\mathcal {C}_b(E)\).
Definition 3
The family \(\mathcal {P}\) of probability measures on \(\left( E, \mathcal {B}(E)\right) \) is said to be tight if for arbitrary \(\varepsilon > 0\) there exists a compact set \(K_\varepsilon \subset E\) such that
We used the following Jakubowski’s generalisation of the Skorokhod Theorem, in the form given by Brzeźniak and Ondreját [14, Theorem C.1], see also [31], as our topological space \(\mathcal {Z}_T\) is not a metric space.
Theorem 6
Let \(\mathcal {X}\) be a topological space such that there exists a sequence \(\{f_m\}_{m \in \mathbb {N}}\) of continuous functions \(f_m : \mathcal {X} \rightarrow \mathbb {R}\) that separates points of \(\mathcal {X}\). Let us denote by \(\mathcal {S}\) the \(\sigma \)-algebra generated by the maps \(\{f_m\}\). Then
-
(a)
every compact subset of \(\mathcal {X}\) is metrizable,
-
(b)
if \((\mu _m)_{m \in \mathbb {N}}\) is a tight sequence of probability measures on \((\mathcal {X}, \mathcal {S})\), then there exists a subsequence \((m_k)_{k \in \mathbb {N}}\), a probability space \((\varOmega , \mathcal {F}, \mathbb {P})\) with \(\mathcal {X}\)-valued Borel measurable variables \(\xi _k, \xi \) such that \(\mu _{m_k}\) is the law of \(\xi _k\) and \(\xi _k\) converges to \(\xi \) almost surely on \(\varOmega \).
Let \((\mathbb {S}, \varrho )\) be a separable and complete metric space.
Definition 4
Let \(u \in \mathcal {C}([0,T]; \mathbb {S})\). The modulus of continuity of u on [0, T] is defined by
Let \((\varOmega , \mathcal {F}, \mathbb {P})\) be a probability space with filtration \(\mathbb {F}:= (\mathcal {F}_t)_{t \in [0,T]}\) satisfying the usual conditions, see [38], and let \((X_n)_{n \in \mathbb {N}}\) be a sequence of continuous \(\mathbb {F}\)-adapted \(\mathbb {S}\)-valued processes.
Definition 5
We say that the sequence \((X_n)_{n \in \mathbb {N}}\) of \(\mathbb {S}\)-valued random variables satisfies condition \([\mathbf {T}]\) iff \(\forall \, \varepsilon>0, \forall \, \eta> 0,\, \exists \, \delta > 0\):
Lemma 32
[12, Lemma 2.4] Assume that \((X_n)_{n \in \mathbb {N}}\) satisfies condition \([\mathbf {T}]\). Let \(\mathbb {P}_n\) be the law of \(X_n\) on \(\mathcal {C}([0,T]; \mathbb {S})\), \(n \in \mathbb {N}\). Then for every \(\varepsilon > 0\) there exists a subset \(A_\varepsilon \subset \mathcal {C}([0,T]; \mathbb {S})\) such that
and
Now we recall the Aldous condition \([\mathbf {A}]\), which is connected with condition \([\mathbf {T}]\) (see [39] and [1]). This condition allows to investigate the modulus of continuity for the sequence of stochastic processes by means of stopped processes.
Definition 6
(Aldous condition) A sequence \((X_n)_{n \in \mathbb {N}}\) satisfies condition \([\mathbf {A}]\) iff \(\,\forall \, \varepsilon > 0\), \(\forall \, \eta > 0\), \(\exists \, \delta > 0\) such that for every sequence \((\tau _n)_{n \in \mathbb {N}}\) of \(\mathbb {F}\)-stopping times with \(\tau _n \le T\) one has
Lemma 33
[39, Theorem 3.2] Conditions \([\mathbf {A}]\) and \([\mathbf {T}]\) are equivalent.
1.2 Tightness Criterion
Now we formulate the compactness criterion analogous to the result due to Mikulevicus and Rozowskii [40], Brzeźniak and Motyl [12] for the space \(\mathcal {Z}_T\), see also [5, Lemma 4.2].
Lemma 34
Let \(\mathcal {Z}_T\), \(\mathcal {T}\) be as defined in (116). Then a set \(\mathcal {K}\subset \mathcal {Z}_T\) is \(\mathcal {T}\)-relatively compact if the following three conditions hold
-
(a)
\(\sup _{u \in \mathcal {K}} \sup _{s \in [0,T]} \Vert u(s)\Vert _{{\mathbb {L}^2({\mathbb {S}^2})}} < \infty ,\)
-
(b)
\(\sup _{u \in \mathcal {K}} \int _0^T \Vert u(s)\Vert ^2_\mathrm { V}\,ds < \infty \,\), i.e. \(\mathcal {K}\) is bounded in \(L^2(0,T; \mathrm { V})\),
-
(c)
\(\lim _{\delta \rightarrow 0} \sup _{u \in \mathcal {K}} \sup _{\underset{|t-s| \le \delta }{s,t \in [0,T]}}\Vert u(t) - u(s)\Vert _{\mathrm {D}(\mathrm { A}^{-1})} = 0.\)
Using Sect. D.1 and the compactness criterion from Lemma 34 we obtain the following corollary.
Corollary 6
(Tightness criterion) Let \((\alpha _\varepsilon )_{\varepsilon > 0}\) be a sequence of continuous \(\mathbb {F}\)-adapted \(\mathrm {H}\)-valued processes such that
-
(a)
there exists a constant \(C_1 > 0\) such that
$$\begin{aligned}\sup _{\varepsilon > 0}\,\mathbb {E}\left[ \sup _{s \in [0,T]} \Vert {\alpha }_\varepsilon (s)\Vert ^2_{\mathrm {H}} \right] \le C_1,\end{aligned}$$ -
(b)
there exists a constant \(C_2 > 0\) such that
$$\begin{aligned}\sup _{\varepsilon > 0}\,\mathbb {E}\left[ \int _0^T \Vert \mathrm{{curl}^\prime }{\alpha }_\varepsilon (s)\Vert ^2_{{\mathbb {L}^2({\mathbb {S}^2})}}\,ds \right] \le C_2,\end{aligned}$$ -
(c)
\(({\alpha }_\varepsilon )_{\varepsilon > 0}\) satisfies the Aldous condition \([\mathbf {A}]\) in \(\mathrm {D}(\mathrm { A}^{-1})\).
Let \(\widetilde{\mathbb {P}}_\varepsilon \) be the law of \({\alpha }_\varepsilon \) on \(\mathcal {Z}_T\). Then for every \(\delta > 0\) there exists a compact subset \(K_\delta \) of \(\mathcal {Z}_T\) such that
Kuratowski Theorem and Proof of Lemma 25
This appendix is dedicated to the proof of Lemma 25. We will first recall the Kuratowski Theorem [33] in the next subsection and prove some related results which will be used later to prove Lemma 25 in Sect. E.2.
1.1 Kuratowski Theorem and Related Results
Theorem 7
Assume that \(X_1, X_2\) are the Polish spaces with their Borel \(\sigma \)-fields denoted respectively by \(\mathcal {B}(X_1), \mathcal {B}(X_2)\). If \(\varphi :X_1 \rightarrow X_2\) is an injective Borel measurable map then for any \(E_1 \in \mathcal {B}(X_1)\), \(E_2 := \varphi (E_1) \in \mathcal {B}(X_2)\).
Next two lemmas are the main results of this appendix. For the proof of Lemma 35 please see [6, Appendix B].
Lemma 35
Let \(X_1, X_2\) and Z be topological spaces such that \(X_1\) is a Borel subset of \(X_2\). Then \(X_1 \cap Z\) is a Borel subset of \(X_2 \cap Z\), where \(X_2 \cap Z\) is a topological space too, with the topology given by
1.2 Proof of Lemma 25
In this subsection we recall Lemma 25 and prove it using the results from previous subsection.
Lemma 36
Let \(T > 0\) and \(\mathcal {Z}_T\) be as defined in (116). Then, the following sets \(\mathcal {C}([0,T];\mathrm { H}) \cap \mathcal {Z}_T\), \(L^2(0,T; \mathrm { V}) \cap \mathcal {Z}_T\) are Borel subsets of \(\mathcal {Z}_T\).
Proof
First of all \(\mathcal {C}([0,T]; \mathrm { H}) \subset \mathcal {C}([0,T]; \mathrm {D}(\mathrm { A}^{-1})) \cap L^2(0,T; \mathrm { H})\). Secondly, \(\mathcal {C}([0,T]; \mathrm { H})\) and \(\mathcal {C}([0,T]; \mathrm {D}(\mathrm { A}^{-1})) \cap L^2(0,T; \mathrm { H})\) are Polish spaces. And finally, since \(\mathrm { H}\) is continuously embedded in \(\mathrm {D}(\mathrm { A}^{-1})\), the map
is continuous and hence Borel. Thus, by application of the Kuratowski Theorem (see Theorem 7), \(\mathcal {C}([0,T]; \mathrm { H})\) is a Borel subset of \(\mathcal {C}([0,T]; \mathrm {D}(\mathrm { A}^{-1})) \cap L^2(0,T; \mathrm { H})\). Therefore, by Lemma 35, \(\mathcal {C}([0,T]; \mathrm { H}) \cap {\mathcal {Z}}_T\) is a Borel subset of \(\mathcal {C}([0,T]; \mathrm {D}(\mathrm { A}^{-1})) \cap L^2(0,T; \mathrm { H}) \cap {\mathcal {Z}}_T\) which is equal to \({\mathcal {Z}}_T\).
Similarly we can show that \(L^2(0,T; \mathrm { V}) \cap {\mathcal {Z}}_T\) is a Borel subset of \({\mathcal {Z}}_T\). \(L^2(0,T; \mathrm { V}) \hookrightarrow L^2(0,T; \mathrm { H})\) and both are Polish spaces thus by application of the Kuratowski Theorem, \(L^2(0,T; \mathrm { V})\) is a Borel subset of \(L^2(0,T; \mathrm { H})\). Finally, we can conclude the proof of lemma by Lemma 35. \(\square \)
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
Brzeźniak, Z., Dhariwal, G. & Le Gia, Q.T. Stochastic Navier–Stokes Equations on a Thin Spherical Domain. Appl Math Optim 84, 1971–2035 (2021). https://doi.org/10.1007/s00245-020-09702-2
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00245-020-09702-2