Search

We establish the local and global well-posedness of strong solutions to the two- and three-dimensional anelastic equations of stratified viscous flows. In this model, the interaction of the density profile with the velocity field is taken into account, and the density background profile is permitted to have physical vacuum singularity. The existing time of the solutions is infinite in two dimensions, with general initial data, and in three dimensions with small initial data., Comment: arXiv admin note: substantial text overlap with arXiv:1903.01523

This work concerns the zero Mach number limit of the compressible primitive equations. The primitive equations with the incompressibility condition are identified as the limiting equations. The convergence with well-prepared initial data (i.e., initial data without acoustic oscillations) is rigorously justified, and the convergence rate is shown to be of order $ \mathcal O(\varepsilon) $, as $ \varepsilon \rightarrow 0^+ $, where $ \varepsilon $ represents the Mach number. As a byproduct, we construct a class of global solutions to the compressible primitive equations, which are close to the incompressible flows.

The Rayleigh-B\'enard 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.

In this paper we prove the local well-posedness and global well-posedness with small initial data of the strong solution to the reduced $3D$ primitive geostrophic adjustment model with weak dissipation. The term reduced model stems from the fact that the relevant physical quantities depends only on two spatial variables. The additional weak dissipation helps us overcome the ill-posedness of original model. We also prove the global well-posedness of the strong solution to the Voigt $\alpha$-regularization of this model, and establish the convergence of the strong solution of the Voigt $\alpha$-regularized model to the corresponding solution of original model. Furthermore, we derive a criterion for finite-time blow-up of reduced $3D$ primitive geostrophic adjustment model with weak dissipation based on Voigt $\alpha$-regularization., Comment: 35 pages

We show that weak solutions of general conservation laws in bounded domains conserve their generalized entropy, and other respective companion laws, if they possess a certain fractional differentiability of order 1/3 in the interior of the domain, and if the normal component of the corresponding fluxes tend to zero as one approaches the boundary. This extends various recent results of the authors.

We study a three-dimensional fluid model describing rapidly rotating convection that takes place in tall columnar structures. The purpose of this model is to investigate the cyclonic and anticyclonic coherent structures. Global existence, uniqueness, continuous dependence on initial data, and large-time behavior of strong solutions are shown provided the model is regularized by a weak dissipation term.

We describe a spectrally-filtered discrete-in-time downscaling data assimilation algorithm and prove, in the context of the two-dimensional Navier--Stokes equations, that this algorithm works for a general class of interpolants, such as those based on local spatial averages as well as point measurements of the velocity. Our algorithm is based on the classical technique of inserting new observational data directly into the dynamical model as it is being evolved over time, rather than nudging, and extends previous results in which the observations were defined directly in terms of an orthogonal projection onto the large-scale (lower) Fourier modes. In particular, our analysis does not require the interpolant to be represented by an orthogonal projection, but requires only the interpolant to satisfy a natural approximation of the identity., Comment: 27 pages, 1 appendix

An intrinsic property of almost any physical measuring device is that it makes observations which are slightly blurred in time. We consider a nudging-based approach for data assimilation that constructs an approximate solution based on a feedback control mechanism that is designed to account for observations that have been blurred by a moving time average. Analysis of this nudging model in the context of the subcritical surface quasi-geostrophic equation shows, provided the time-averaging window is sufficiently small and the resolution of the observations sufficiently fine, that the approximating solution converges exponentially fast to the observed solution over time. In particular, we demonstrate that observational data with a small blur in time possess no significant obstructions to data assimilation provided that the nudging properly takes the time averaging into account. Two key ingredients in our analysis are additional boundedness properties for the relevant interpolant observation operators and a non-local Gronwall inequality., Comment: 44 pages

Using the convex integration technique for the three-dimensional Navier-Stokes equations introduced by T. Buckmaster and V. Vicol, it is shown the existence of non-unique weak solutions for the 3D Navier-Stokes equations with fractional hyperviscosity $(-\Delta)^{\theta}$, whenever the exponent $\theta$ is less than J.-L. Lions' exponent $5/4$, i.e., when $\theta < 5/4$.

We show the existence of global weak solutions to the three-dimensional compressible primitive equations of atmospheric dynamics with degenerate viscosities. In analogy with the case of the compressible Navier-Stokes equations, the weak solutions satisfy the basic energy inequality, the Bresh-Desjardins entropy inequality and the Mellet-Vasseur estimate. These estimates play an important role in establishing the compactness of the vertical velocity of the approximating solutions, and therefore are essential to recover the vertical velocity in the weak solutions.

In this paper we analyze a finite element method applied to a continuous downscaling data assimilation algorithm for the numerical approximation of the two and three dimensional Navier-Stokes equations corresponding to given measurements on a coarse spatial scale. For representing the coarse mesh measurements we consider different types of interpolation operators including a Lagrange interpolant. We obtain uniform-in-time estimates for the error between a finite element approximation and the reference solution corresponding to the coarse mesh measurements. We consider both the case of a plain Galerkin method and a Galerkin method with grad-div stabilization. For the stabilized method we prove error bounds in which the constants do not depend on inverse powers of the viscosity. Some numerical experiments illustrate the theoretical results.

This work is devoted to establishing the local-in-time well-posedness of strong solutions to the three-dimensional compressible primitive equations of atmospheric dynamics. It is shown that strong solutions exist, unique, and depend continuously on the initial data, for a short time in two cases: with gravity but without vacuum, and with vacuum but without gravity. We also introduce the free boundary problem for the compressible primitive equations.

The aim of this work is to extend and prove the Onsager conjecture for a class of conservation laws that possess generalized entropy. One of the main findings of this work is the "universality" of the Onsager exponent, $\alpha > 1/3$, concerning the regularity of the solutions, say in $C^{0,\alpha}$, that guarantees the conservation of the generalized entropy; regardless of the structure of the genuine nonlinearity in the underlying system.

In this paper we show the global well-posedness of solutions to a three-dimensional magnetohydrodynamical (MHD) model in porous media. Compared to the classical MHD equations, our system involves a nonlinear damping term in the momentum equations due to the "Brinkman-Forcheimer-extended-Darcy" law of flow in porous media.

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.

We study the three-dimensional Hasegawa-Mima model of turbulent magnetized plasma with horizontal viscous terms and a weak vertical dissipative term. In particular, we establish the global existence and uniqueness of strong solutions for this model.

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.

We introduce a continuous (downscaling) data assimilation algorithm for the 2D B\'enard convection problem using vorticity or local circulation measurements only. In this algorithm, a nudging term is added to the vorticity equation to constrain the model. Our numerical results indicate that the approximate solution of the algorithm is converging to the unknown reference solution (vorticity and temperature) corresponding to the measurements of the 2D B\'enard convection problem when only spatial coarse-grain measurements of vorticity are assimilated. Moreover, this convergence is realized using data which is much more coarse than the resolution needed to satisfy rigorous analytical estimates.

The goal of this note is to show that, also in a bounded domain $\Omega \subset \mathbb{R}^n$, with $\partial \Omega\in C^2$, any weak solution, $(u(x,t),p(x,t))$, of the Euler equations of ideal incompressible fluid in $\Omega\times (0,T) \subset \mathbb{R}^n\times\mathbb{R}_t$, with the impermeability boundary condition: $u\cdot \vec n =0$ on $\partial\Omega\times(0,T)$, is of constant energy on the interval $(0,T)$ provided the velocity field $u \in L^3((0,T); C^{0,\alpha}(\overline{\Omega}))$, with $\alpha>\frac13\,.$

An important feature of the planetary oceanic dynamics is that the aspect ratio (the ratio of the depth to horizontal width) is very small. As a result, the hydrostatic approximation (balance), derived by performing the formal small aspect ratio limit to the Navier-Stokes equations, is considered as a fundamental component in the primitive equations of the large-scale ocean. In this paper, we justify rigorously the small aspect ratio limit of the Navier-Stokes equations to the primitive equations. Specifically, we prove that the Navier-Stokes equations, after being scaled appropriately by the small aspect ratio parameter of the physical domain, converge strongly to the primitive equations, globally and uniformly in time, and the convergence rate is of the same order as the aspect ratio parameter. This result validates the hydrostatic approximation for the large-scale oceanic dynamics. Notably, only the weak convergence of this small aspect ratio limit was rigorously justified before., Comment: arXiv admin note: text overlap with arXiv:1604.01695

In this paper we introduce a finite-parameters feedback control algorithm for stabilizing solutions of the Navier-Stokes-Voigt equations, the strongly damped nonlinear wave equations and the nonlinear wave equation with nonlinear damping term, the Benjamin-Bona-Mahony-Burgers equation and the KdV-Burgers equation. This algorithm capitalizes on the fact that such infinite-dimensional dissipative dynamical systems posses finite-dimensional long-time behavior which is represented by, for instance, the finitely many determining parameters of their long-time dynamics, such as determining Fourier modes, determining volume elements, determining nodes , etc..The algorithm utilizes these finite parameters in the form of feedback control to stabilize the relevant solutions. For the sake of clarity, and in order to fix ideas, we focus in this work on the case of low Fourier modes feedback controller, however, our results and tools are equally valid for using other feedback controllers employing other spatial coarse mesh interpolants.

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

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.

In this paper we survey the various implementations of a new data assimilation (downscaling) algorithm based on spatial coarse mesh measurements. As a paradigm, we demonstrate the application of this algorithm to the 3D Leray-$\alpha$ subgrid scale turbulence model. Most importantly, we use this paradigm to show that it is not always necessary that one has to collect coarse mesh measurements of all the state variables, that are involved in the underlying evolutionary system, in order to recover the corresponding exact reference solution. Specifically, we show that in the case of the 3D Leray$-\alpha$ model of turbulence the solutions of the algorithm, constructed using only coarse mesh observations of any two components of the three-dimensional velocity field, and without any information of the third component, converge, at an exponential rate in time, to the corresponding exact reference solution of the 3D Leray$-\alpha$ model. This study serves as an addendum to our recent work on abridged continuous data assimilation for the 2D Navier-Stokes equations. Notably, similar results have also been recently established for the 3D viscous Planetary Geostrophic circulation model in which we show that coarse mesh measurements of the temperature alone are sufficient for recovering, through our data assimilation algorithm, the full solution; viz. the three components of velocity vector field and the temperature. Consequently, this proves the Charney conjecture for the 3D Planetary Geostrophic model; namely, that the history of the large spatial scales of temperature is sufficient for determining all the other quantities (state variables) of the model.