The three limits of the hydrostatic approximation
Ken Furukawa, Yoshikazu Giga, Matthias Hieber, Amru Hussein, Takahito Kashiwabara, Marc Wrona
TL;DR
This work analyzes the hydrostatic approximation by studying the three formal limits arising from anisotropic scaling ν_z=ε^2δ in a vertically thin domain for the 3D Navier–Stokes equations. Using maximal L^2-regularity and a quadratic-inequality framework, it proves norm-convergence to (i) the primitive equations with full viscosity (δ finite, ε→0), (ii) the primitive equations with only horizontal viscosity (δ→0), and (iii) the 2D-Navier–Stokes equations (δ→∞) as ε→0, including convergence for ε→0 with δ fixed and, in some regimes, explicit rates like O(ε+δ) or δ^{-1/4}. The analysis relies on uniform linear estimates for the scaled Stokes problem, preservation of vertical regularity, and detailed nonlinear estimates in anisotropic function spaces, enabling a flexible treatment of multiple coupled limits. The results substantiate the stability of reduced models under parameter variations and provide a method applicable to related geophysical models with anisotropic diffusion and transport. The work also emphasizes the versatility of the maximal-regularity approach over purely energy-based methods in capturing sharp convergence behavior across regimes.
Abstract
The primitive equations are derived from the $3D$-Navier-Stokes equations by the hydrostatic approximation. Formally, assuming an $\varepsilon$-thin domain and anisotropic viscosities with vertical viscosity $ν_z=\mathcal{O}(\varepsilon^γ)$ where $γ=2$, one obtains the primitive equations with full viscosity as $\varepsilon\to 0$. Here, we take two more limit equations into consideration: For $γ<2$ the $2D$-Navier-Stokes equations are obtained. For $γ>2$ the primitive equations with only horizontal viscosity $-Δ_H$ as $\varepsilon\to 0$. Thus, there are three possible limits of the hydrostatic approximation depending on the assumption on the vertical viscosity. The latter convergence has been proven recently by Li, Titi, and Yuan using energy estimates. Here, we consider more generally $ν_z=\varepsilon^2 δ$ and show how maximal regularity methods and quadratic inequalities can be an efficient approach to the same end for $\varepsilon,δ\to 0$. The flexibility of our methods is also illustrated by the convergence for $δ\to \infty$ and $\varepsilon\to 0$ to the $2D$-Navier-Stokes equations.
