Search

We prove global well-posedness of the ocean primitive equations coupled to advection-diffusion equations of the oceanic tracers temperature and salinity that are supplemented by the eddy parametrization model due to Gent-McWilliams and Redi. This parametrization forms a milestone in global ocean modelling and constitutes a central part of any general ocean circulation model computation. The eddy parametrization adds a secondary transport velocity to the tracer equation and renders the original Laplacian operators in the advection-diffusion equations nonlinear, with a diffusion matrix that depends via the equation of state in a nonlinear fashion on both tracers simultaneously. The eddy parametrization of Gent-McWilliams-Redi augments the complexity of the mathematical analysis of the whole system which we present here. We show first that weak solutions exist globally in time, provided the parametrization uses a regularized density. Then we prove by a detailed analysis of the eddy operators the global well-posedness. Our results apply also to the ``small-slope approximation'' that is commonly used in global ocean simulations., Comment: 31 pages

We consider the three-dimensional incompressible Euler equations on a bounded domain $Ω$ with $C^3$ boundary. We prove that if the velocity field $u \in C^{0,α} (Ω)$ with $α> 0$ (where we are omitting the time dependence), it follows that the pressure $p \in C^{0, α} (Ω)$. In order to prove this result we use a local parametrisation of the boundary and a very weak formulation of the boundary condition for the pressure, as was introduced in [C. Bardos and E.S. Titi, Philos. Trans. Royal Soc. A, 380 (2022), 20210073]. Moreover, we provide an example illustrating the necessity of this new very weak formulation of the boundary condition for the pressure. This result is of importance for the proof of the first half of the Onsager Conjecture, the sufficient conditions for energy conservation of weak solutions to the three-dimensional incompressible Euler equations in bounded domains., 35 pages

This paper is devoted to investigating the rotating Boussinesq equations of inviscid, incompressible flows with both fast Rossby waves and fast internal gravity waves. The main objective is to establish a rigorous derivation and justification of a new generalized quasi-geostrophic approximation in a channel domain with no normal flow at the upper and lower solid boundaries, taking into account the resonance terms due to the fast and slow waves interactions. Under these circumstances, We are able to obtain uniform estimates and compactness without the requirement of either well-prepared initial data (as in [10]) or domain with no boundary (as in [17]). In particular, the nonlinear resonances and the new limit system, which takes into account the fast waves correction to the slow waves dynamics, are also identified without introducing Fourier series expansion. The key ingredient includes the introduction of (full) generalized potential vorticity.

A compressible, viscous and heat conducting fluid is confined between two parallel plates maintained at a constant temperature and subject to a strong stratification due to the gravitational force. We consider the asymptotic limit, where the Mach number and the Froude number are of the same order proportional to a small parameter. We show the limit problem can be identified with Majda's model of layered ``stack-of-pancake'' flow., Comment: arXiv admin note: text overlap with arXiv:2206.14041

The purpose of this article is to give a complete proof of a [Formula: see text] regularity result for the pressure for weak solutions of the two-dimensional 'incompressible Euler equations' when the fluid velocity enjoys the same type of regularity in a compact simply connected domain with [Formula: see text] boundary. To accomplish our result, we realize that it is compulsory to introduce a new weak formulation for the boundary condition of the pressure, which is consistent with, and equivalent to, that of classical solutions. This article is part of the theme issue 'Scaling the turbulence edifice (part 1)'.

In this paper, we provide rigorous justification of the hydrostatic approximation and the derivation of primitive equations as the small aspect ratio limit of the incompressible three-dimensional Navier-Stokes equations in the anisotropic horizontal viscosity regime. Setting $\varepsilon >0$ to be the small aspect ratio of the vertical to the horizontal scales of the domain, we investigate the case when the horizontal and vertical viscosities in the incompressible three-dimensional Navier-Stokes equations are of orders $O(1)$ and $O(\varepsilon^\alpha)$, respectively, with $\alpha>2$, for which the limiting system is the primitive equations with only horizontal viscosity as $\varepsilon$ tends to zero. In particular we show that for "well prepared" initial data the solutions of the scaled incompressible three-dimensional Navier-Stokes equations converge strongly, in any finite interval of time, to the corresponding solutions of the anisotropic primitive equations with only horizontal viscosities, as $\varepsilon$ tends to zero, and that the convergence rate is of order $O\left(\varepsilon^\frac\beta2\right)$, where $\beta=\min\{\alpha-2,2\}$. Note that this result is different from the case $\alpha=2$ studied in [Li, J.; Titi, E.S.: \emph{The primitive equations as the small aspect ratio limit of the Navier-Stokes equations: Rigorous justification of the hydrostatic approximation}, J. Math. Pures Appl., \textbf{124} \rm(2019), 30--58], where the limiting system is the primitive equations with full viscosities and the convergence is globally in time and its rate of order $O\left(\varepsilon\right)$.

Differential buoyancy surface sources in the ocean may induce a density-driven flow that joins faster flow components to create a multi-scale, 3D flow. Potential temperature and salinity are active tracers that determine the ocean's potential density: their distribution strongly affects the density-driven component while the overall flow affects their distribution. We present a robust framework that allows one to study the effects of a general 3D flow on a density-driven velocity component, by constructing a modular observation-based 3D model of intermediate complexity. The model contains an incompressible velocity that couples two advection-diffusion equations, for temperature and salinity. Instead of solving the Navier-Stokes equations for the velocity, we consider a flow composed of several temporally separated, spatially predetermined modes. One of these modes models the density-driven flow: its spatial form describes the density-driven flow structure and its strength is determined dynamically by average density differences. The other modes are completely predetermined, consisting of any incompressible, possibly unsteady, 3D flow, e.g. as determined by kinematic models, observations, or simulations. The model is a non-linear, weakly coupled system of two non-local PDEs. We prove its well-posedness in the sense of Hadamard, and obtain rigorous bounds regarding analytical solutions. The model's relevance to oceanic systems is demonstrated by tuning the model to mimic the North Atlantic ocean's dynamics. In one limit the model recovers a simplified oceanic box model and in another limit a kinematic model of oceanic chaotic advection, suggesting it can be utilized to study spatially dependent feedback processes in the ocean., Comment: 15 pages text, 20 pages appendix, 5 figures

We construct a determining form for the 2D Rayleigh-B��nard (RB) system in a strip with solid horizontal boundaries, in the cases of no-slip and stress-free boundary conditions. The determining form is an ODE in a Banach space of trajectories whose steady states comprise the long-time dynamics of the RB system. In fact, solutions on the global attractor of the RB system can be further identified through the zeros of a scalar equation to which the ODE reduces for each initial trajectory. The twist in this work is that the trajectories are for the velocity field only, which in turn determines the corresponding trajectories of the temperature.

The Rayleigh-B��nard system with stress-free boundary conditions is shown to have a global attractor in each affine space where velocity has fixed spatial average. The physical problem is shown to be equivalent to one with periodic boundary conditions and certain symmetries. This enables a Gronwall estimate on enstrophy. That estimate is then used to bound the $L^2$ norm of the temperature gradient on the global attractor, which, in turn, is used to find a bounding region for the attractor in the enstrophy, palinstrophy-plane. All final bounds are algebraic in the viscosity and thermal diffusivity, a significant improvement over previously established estimates. The sharpness of the bounds are tested with numerical simulations.

We consider fully discrete numerical schemes for a downscaling data assimilation algorithm aimed at approximating the velocity field of the 2D Navier-Stokes equations corresponding to given coarse mesh observational measurements. The time discretization is done by considering semi- and fully-implicit Euler schemes, and the spatial discretization is based on a spectral Galerkin method. The two fully discrete algorithms are shown to be unconditionally stable, with respect to the size of the time step, number of time steps and the number of Galerkin modes. Moreover, explicit, uniform in time error estimates between the fully discrete solution and the reference solution corresponding to the observational coarse mesh measurements are obtained, in both the $L^2$ and $H^1$ norms. Notably, the two-dimensional Navier-Stokes equations, subject to the no-slip Dirichlet or periodic boundary conditions, are used in this work as a paradigm. The complete analysis that is presented here can be extended to other two- and three-dimensional dissipative systems under the assumption of global existence and uniqueness.

Based on a previously introduced downscaling data assimilation algorithm, which employs a nudging term to synchronize the coarse mesh spatial scales, we construct a determining map for recovering the full trajectories from their corresponding coarse mesh spatial trajectories, and investigate its properties. This map is then used to develop a downscaling data assimilation scheme for statistical solutions of the two-dimensional Navier-Stokes equations, where the coarse mesh spatial statistics of the system is obtained from discrete spatial measurements. As a corollary, we deduce that statistical solutions for the Navier-Stokes equations are determined by their coarse mesh spatial distributions. Notably, we present our results in the context of the Navier-Stokes equations; however, the tools are general enough to be implemented for other dissipative evolution equations.

In this paper, we consider the 3D primitive equations of oceanic and atmospheric dynamics with only horizontal eddy viscosities in the horizontal momentum equations and only vertical diffusivity in the temperature equation. Global well-posedness of strong solutions is established for any initial data such that the initial horizontal velocity $v_0\in H^2(\Omega)$ and the initial temperature $T_0\in H^1(\Omega)\cap L^\infty(\Omega)$ with $\nabla_HT_0\in L^q(\Omega)$, for some $q\in(2,\infty)$. Moreover, the strong solutions enjoy correspondingly more regularities if the initial temperature belongs to $H^2(\Omega)$. The main difficulties are the absence of the vertical viscosity and the lack of the horizontal diffusivity, which, interact with each other, thus causing the "\,mismatching\," of regularities between the horizontal momentum and temperature equations. To handle this "mismatching" of regularities, we introduce several auxiliary functions, i.e., $\eta, \theta, \varphi,$ and $\psi$ in the paper, which are the horizontal curls or some appropriate combinations of the temperature with the horizontal divergences of the horizontal velocity $v$ or its vertical derivative $\partial_zv$. To overcome the difficulties caused by the absence of the horizontal diffusivity, which leads to the requirement of some $L^1_t(W^{1,\infty}_\textbf{x})$-type a priori estimates on $v$, we decompose the velocity into the "temperature-independent" and temperature-dependent parts and deal with them in different ways, by using the logarithmic Sobolev inequalities of the Br\'ezis-Gallouet-Wainger and Beale-Kato-Majda types, respectively. Specifically, a logarithmic Sobolev inequality of the limiting type, introduced in our previous work [12], is used, and a new logarithmic type Gronwall inequality is exploited.

We construct a determining form for the surface quasi-geostrophic (SQG) equation with subcritical dissipation. In particular, we show that the global attractor for this equation can be embedded in the long-time dynamics of an ordinary differential equation (ODE) called a determining form. Indeed, there is a one-to-one correspondence between the trajectories in the global attractor of the SQG equation and the steady state solutions of the determining form. The determining form is a true ODE in the sense that its vector field is Lipschitz. This is shown by combining De Giorgi techniques and elementary harmonic analysis. Finally, we provide elementary proofs of the existence of time-periodic solutions, steady state solutions, as well as the existence of finitely many determining parameters for the SQG equation., Comment: 34 pages

The paper studies the issue of stability of solutions to the Navier-Stokes and damped Euler systems in periodic boxes. We show that under action of fast oscillating-in- time external forces all two dimensional regular solutions converge to a time periodic flow. Unexpectedly, effects of stabilization can be also obtained for systems with stationary forces with large total momentum (average of the velocity). Thanks to the Galilean transformation and space boundary conditions, the stationary force changes into one with time oscillations. In the three dimensional case we show an analogical result for weak solutions to the Navier- Stokes equations.

We apply the Postprocessing Galerkin method to a recently introduced continuous data assimilation (downscaling) algorithm for obtaining a numerical approximation of the solution of the two-dimensional Navier-Stokes equations corresponding to given measurements from a coarse spatial mesh. Under suitable conditions on the relaxation (nudging) parameter, the resolution of the coarse spatial mesh and the resolution of the numerical scheme, we obtain uniform in time estimates for the error between the numerical approximation given by the Postprocessing Galerkin method and the reference solution corresponding to the measurements. Our results are valid for a large class of interpolant operators, including low Fourier modes and local averages over finite volume elements. Notably, we use here the 2D Navier-Stokes equations as a paradigm, but our results apply equally to other evolution equations, such as the Boussinesq system of Benard convection and other oceanic and atmospheric circulation models.

An algorithm for continuous data assimilation for the two- dimensional B��nard convection problem is introduced and analyzed. It is inspired by the data assimilation algorithm developed for the Navier-Stokes equations, which allows for the implementation of variety of observables: low Fourier modes, nodal values, finite volume averages, and finite elements. The novelty here is that the observed data is obtained for the velocity field alone; i.e. no temperature measurements are needed for this algorithm. We provide conditions on the spatial resolution of the observed data, under the assumption that the observed data is free of noise, which are sufficient to show that the solution of the algorithm approaches, at an exponential rate, the unique exact unknown solution of the B��nard convection problem associated with the observed (finite dimensional projection of) velocity.

Motivated by the presence of a finite number of determining parameters (degrees of freedom) such as modes, nodes and local spatial averages for dissipative dynamical systems, we present a continuous data assimilation algorithm for the three-dimensional Navier-Stokes-$��$ model. This algorithm consists of introducing a nudging process through general type of approximation interpolation operator (that is constructed from observational measurements) that synchronizes the large spatial scales of the approximate solutions with those of the unknown solutions the Navier-Stokes-$��$ equations that corresponds to these measurements. Our main result provides conditions on the finite-dimensional spatial resolution of the collected data, sufficient to guarantee that the approximating solution, that is obtained from these collected data, converges to the unkown reference solution (physical reality) over time. These conditions are given in terms of some physical parameters, such as kinematic viscosity, the size of the domain and the forcing term.

In this article we consider the Euler-$��$ system as a regularization of the incompressible Euler equations in a smooth, two-dimensional, bounded domain. For the limiting Euler system we consider the usual non-penetration boundary condition, while, for the Euler-$��$ regularization, we use velocity vanishing at the boundary. We also assume that the initial velocities for the Euler-$��$ system approximate, in a suitable sense, as the regularization parameter $��\to 0$, the initial velocity for the limiting Euler system. For small values of $��$, this situation leads to a boundary layer, which is the main concern of this work. Our main result is that, under appropriate regularity assumptions, and despite the presence of this boundary layer, the solutions of the Euler-$��$ system converge, as $��\to 0$, to the corresponding solution of the Euler equations, in $L^2$ in space, uniformly in time. We also present an example involving parallel flows, in order to illustrate the indifference to the boundary layer of the $��\to 0$ limit, which underlies our work., 22pages

This paper is concerned with the cubic Szeg�� equation $$ i\partial_t u=��(|u|^2 u), $$ defined on the $L^2$ Hardy space on the one-dimensional torus $\mathbb T$, where $��: L^2(\mathbb T)\rightarrow L^2_+(\mathbb T)$ is the Szeg�� projector onto the non-negative frequencies. For analytic initial data, it is shown that the solution remains spatial analytic for all time $t\in (-\infty,\infty)$. In addition, we find a lower bound for the radius of analyticity of the solution. Our method involves energy-like estimates of the special Gevrey class of analytic functions based on the $\ell^1$ norm of Fourier transforms (the Wiener algebra).

We establish global existence and uniqueness theorems for the two-dimensional non-diffusive Boussinesq system with viscosity only in the horizontal direction, which arises in Ocean dynamics. This work improves the global well-posedness results established recently by R. Danchin and M. Paicu for the Boussinesq system with anisotropic viscosity and zero diffusion. Although we follow some of their ideas, in proving the uniqueness result, we have used an alternative approach by writing the transported temperature (density) as $��= ����$ and adapting the techniques of V. Yudovich for the 2D incompressible Euler equations. This new idea allows us to establish uniqueness results with fewer assumptions on the initial data for the transported quantity $��$. Furthermore, this new technique allows us to establish uniqueness results without having to resort to the paraproduct calculus of J. Bony. We also propose an inviscid $��$-regularization for the two-dimensional inviscid, non-diffusive Boussinesq system of equations, which we call the Boussinesq-Voigt equations. Global regularity of this system is established. Moreover, we establish the convergence of solutions of the Boussinesq-Voigt model to the corresponding solutions of the two-dimensional Boussinesq system of equations for inviscid flow without heat (density) diffusion on the interval of existence of the latter. Furthermore, we derive a criterion for finite-time blow-up of the solutions to the inviscid, non-diffusive 2D Boussinesq system based on this inviscid Voigt regularization. Finally, we propose a Voigt-$��$ regularization for the inviscid 3D Boussinesq equations with diffusion, and prove its global well-posedness. It is worth mentioning that our results are also valid in the presence of the $��$-plane approximation of the Coriolis force.

In this paper we develop and use successive averaging methods for explaining the regularization mechanism in the the periodic Korteweg--de Vries (KdV) equation in the homogeneous Sobolev spaces $\dot{H}^s$, for $s\ge0$. Specifically, we prove the global existence, uniqueness, and Lipschitz continuous dependence on the initial data of the solutions of the periodic KdV. For the case where the initial data is in $L_2$ we also show the Lipschitz continuous dependence of these solutions with respect to the initial data as maps from $\dot{H}^s$ to $\dot{H}^s$, for $s\in(-1,0]$.

The global existence of weak solutions for the three-dimensional axisymmetric Euler-$��$ (also known as Lagrangian-averaged Euler-$��$) equations, without swirl, is established, whenever the initial unfiltered velocity $v_0$ satisfies $\frac{\nabla \times v_0}{r}$ is a finite Randon measure with compact support. Furthermore, the global existence and uniqueness, is also established in this case provided $\frac{\nabla \times v_0}{r} \in L^p_c(\mathbb{R}^3)$ with $p>{3/2}$. It is worth mention that no such results are known to be available, so far, for the three-dimensional Euler equations of ideal incompressible flows.

We present an alpha-regularization of the Birkhoff-Rott equation, induced by the two-dimensional Euler-alpha equations, for the vortex sheet dynamics. We show the convergence of the solutions of Euler-alpha equations to a weak solution of the Euler equations for initial vorticity being a finite Radon measure of fixed sign, which includes the vortex sheets case. We also show that, provided the initial density of vorticity is an integrable function over the curve with respect to the arc-length measure, (i) an initially Lipschitz chord arc vortex sheet (curve), evolving under the BR-alpha equation, remains Lipschitz for all times, (ii) an initially Holder C^{1,beta}, 0, 36 pages, no figures, will appear in CPAM, corrected typos, minor clarifications added

The nonlinear Schr��dinger-Helmholtz (SH) equation in $N$ space dimensions with $2��$ nonlinear power was proposed as a regularization of the classical nonlinear Schr��dinger (NLS) equation. It was shown that the SH equation has a larger regime ($1\le��, new article, 20 pages, 10 figures, to appear in Numerical Functional Analysis and Optimization

We display a continous equation for the deconvolution process that generalizes the Van Cittert algorithm in the case of oceanic boundary conditions for a given fixed wind. We deduce a LES model for which we have existence and uniqueness of a strong solution. Finally, we display several numerical simulations showing the practical interest of the model.

We present data from high-resolution numerical simulations of the Navier-Stokes-$��$ and the Leray-$��$ models for two-dimensional turbulence. It was shown previously (Lunasin et al., J. Turbulence, 8, (2007), 751-778), that for wavenumbers $k$ such that $k��\gg 1$, the energy spectrum of the smoothed velocity field for the two-dimensional Navier-Stokes-$��$ (NS-$��$) model scales as $k^{-7}$. This result is in agreement with the scaling deduced by dimensional analysis of the flux of the conserved enstrophy using its characteristic time scale. We therefore hypothesize that the spectral scaling of any $��$-model in the sub-$��$ spatial scales must depend only on the characteristic time scale and dynamics of the dominant cascading quantity in that regime of scales. The data presented here, from simulations of the two-dimensional Leray-$��$ model, confirm our hypothesis. We show that for $k��\gg 1$, the energy spectrum for the two-dimensional Leray-$��$ scales as $k^{-5}$, as expected by the characteristic time scale for the flux of the conserved enstrophy of the Leray-$��$ model. These results lead to our conclusion that the dominant directly cascading quantity of the model equations must determine the scaling of the energy spectrum., 11 pages, 4 figures

Recently, the Navier-Stokes-Voight (NSV) model of viscoelastic incompressible fluid has been proposed as a regularization of the 3D Navier-Stokes equations for the purpose of direct numerical simulations. In this work we prove that the global attractor of the $D NSV equations, driven by an analytic forcing, consists of analytic functions. A consequence of this result is that the spectrum of the solutions of the 3D NSV system, lying on the global attractor, have exponentially decaying tail, despite the fact that the equations behave like a damped hyperbolic system, rather than the parabolic one. This result provides an additional evidence that the 3D NSV with the small regularization parameter enjoys similar statistical properties as the 3D Navier-Stokes equations. Finally, we calculate a lower bound for the exponential decaying scale -- the scale at which the spectrum of the solution start to decay exponentially, and establish a similar bound for the steady state solutions of the 3D NSV and 3D Navier-Stokes equations. Our estimate coincides with similar available lower bound for the smallest dissipation length scale of solutions of the 3D Navier-Stokes equations.

Inspired by recent developments in Berdina-like models for turbulence, we propose an inviscid regularization for the surface quasi-geostrophic (SQG) equations. We are particularly interested in the celebrated question of blowup in finite time of the solution gradient of the SQG equations. The new regularization yields a necessary and sufficient condition, satisfied by the regularized solution, when a regularization parameter $\alpha$ tends to zero, for the solution of the original SQG equations to develop a singularity in finite time. As opposed to the commonly used viscous regularization, the inviscid equations derived here conserve a modified energy. Therefore, the new regularization provides an attractive numerical procedure for finite time blow up testing. In particular, we prove that, if the initial condition is smooth, then the regularized solution remains as smooth as the initial data for all times. Moreover, much like the original problem, the inviscid regularization admits a maximum principle.

In this paper we consider the three-dimensional Navier-Stokes equations in an infinite channel. We provide a sufficient condition, in terms of $\partial_z p$, where $p$ is the pressure, for the global existence of the strong solutions to the three-dimensional Navier-Stokes equations.

In this paper we study a three dimensional thermocline planetary geostrophic ``horizontal" hyper--diffusion model of the gyre-scale midlatitude ocean. We show the global existence and uniqueness of the weak and strong solutions to this model. Moreover, we establish the existence of a finite dimensional global attractor to this dissipative evolution system. Preliminary computational tests indicate that our hyper--diffusion model does not exhibit any of the nonphysical instabilities near the literal boundary which are observed numerically in other models., Comment: 32 pages