Source author record

Edriss S. Titi

Edriss S. Titi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

77works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

77 published item(s)

preprint2026arXiv

Plastic limit of a viscoplastic Burgers equation -- A toy model for sea-ice dynamics

We study the plastic Burgers equation in one space dimension, i.e., the Burgers equation featuring an additional term formally given by the p-Laplacian with p=1, or rather, by the multivalued subdifferential of the total variation functional. Our study highlights that the interplay of the advection term with the stresses given by the multivalued 1-Laplacian is a crucial feature of this model. Eventhough it is an interesting model in itsef, it can also be regarded as a one-dimensional version of the momentum balance of Hibler's model for sea-ice dynamics. Therein, the stress tensor is given by a term with similar properties as the 1-Laplacian in order to account for plastic effects of the ice. For our analysis we start out from a viscoplastic Burgers equation, i.e., a suitably regularized version of the plastic Burgers equation with a small regularization parameter $\varepsilon>0$. For the viscoplastic Burgers equation, we construct a global BV-solution. In the singular limit $\varepsilon\to0$ we deduce the existence of a BV-solution for the plastic Burgers equation. In addition we show that the term arising as the limit of the regularized stresses is indeed related to an element of the subdifferential of the total variation functional.

preprint2024arXiv

Data Assimilation in Chaotic Systems Using Deep Reinforcement Learning

Data assimilation (DA) plays a pivotal role in diverse applications, ranging from climate predictions and weather forecasts to trajectory planning for autonomous vehicles. A prime example is the widely used ensemble Kalman filter (EnKF), which relies on linear updates to minimize variance among the ensemble of forecast states. Recent advancements have seen the emergence of deep learning approaches in this domain, primarily within a supervised learning framework. However, the adaptability of such models to untrained scenarios remains a challenge. In this study, we introduce a novel DA strategy that utilizes reinforcement learning (RL) to apply state corrections using full or partial observations of the state variables. Our investigation focuses on demonstrating this approach to the chaotic Lorenz '63 system, where the agent's objective is to minimize the root-mean-squared error between the observations and corresponding forecast states. Consequently, the agent develops a correction strategy, enhancing model forecasts based on available system state observations. Our strategy employs a stochastic action policy, enabling a Monte Carlo-based DA framework that relies on randomly sampling the policy to generate an ensemble of assimilated realizations. Results demonstrate that the developed RL algorithm performs favorably when compared to the EnKF. Additionally, we illustrate the agent's capability to assimilate non-Gaussian data, addressing a significant limitation of the EnKF.

preprint2022arXiv

A dynamic-kinematic 3D model for density-driven ocean circulation flows: Construction, global well-posedness and dynamics

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.

preprint2022arXiv

Enhanced Simulation of the Indian Summer Monsoon Rainfall Using Regional Climate Modeling and Continuous Data Assimilation

This study assesses a Continuous Data Assimilation (CDA) dynamical-downscaling algorithm for enhancing the simulation of the Indian summer monsoon (ISM) system. CDA is a mathematically rigorous technique that has been recently introduced to constrain the large-scale features of high-resolution atmospheric models with coarse spatial scale data. It is similar to spectral nudging but does not require any spectral decomposition for scales separation. This is expected to be particularly relevant for ISM, which involves various interactions between large-scale circulations and regional physical processes. Along with a control simulation, several downscaling simulations were conducted with the Weather Research and Forecasting (WRF) model using CDA, spectral (retaining different wavenumbers) and grid nudging for three ISM seasons: normal (2016), excess (2013), and drought (2009). The simulations are nested within the NCEP Final Analysis and the model outputs are evaluated against the observations. Compared to grid and spectral nudging, the simulations using CDA produce enhanced ISM features over the Indian subcontinent including the low-level jet, tropical easterly jet, easterly wind shear, and rainfall distributions for all investigated ISM seasons. The major ISM processes, in particular the monsoon inversion over the Arabian Sea, tropospheric temperature gradients and moist static energy over central India, and zonal wind shear over the monsoon region, are all better simulated with CDA. Spectral nudging outputs are found to be sensitive to the choice of the wavenumber, requiring careful tuning to provide robust simulations of the ISM system. In contrast, control and grid nudging generally fail to well reproduce some of the main ISM features.

preprint2022arXiv

Global well-posedness of a three-dimensional Brinkman-Forchheimer-Bénard convection model in porous media

We consider three-dimensional (3D) Boussinesq convection system of an incompressible fluid in a closed sample of a porous medium. Specifically, we introduce and analyze a 3D Brinkman-Forchheimer-Bénard convection problem describing the behavior of an incompressible fluid in a porous medium between two plates heated from the bottom and cooled from the top. We show the existence and uniqueness of global in-time solutions, and the existence of absorbing balls in $L^2$ and $H^1$. Eventually, we comment on the applicability of a data assimilation algorithm to our system.

preprint2022arXiv

On the effect of fast rotation and vertical viscosity on the lifespan of the $3D$ primitive equations

We study the effect of the fast rotation and vertical viscosity on the lifespan of solutions to the three-dimensional primitive equations (also known as the hydrostatic Navier-Stokes equations) with impermeable and stress-free boundary conditions. Firstly, for a short time interval, independent of the rate of rotation $|Ω|$, we establish the local well-posedness of solutions with initial data that is analytic in the horizontal variables and only $L^2$ in the vertical variable. Moreover, it is shown that the solutions immediately become analytic in all the variables with increasing-in-time (at least linearly) radius of analyticity in the vertical variable for as long as the solutions exist. On the other hand, the radius of analyticity in the horizontal variables might decrease with time, but as long as it remains positive the solution exists. Secondly, with fast rotation, i.e., large $|Ω|$, we show that the existence time of the solution can be prolonged, with "well-prepared" initial data. Finally, in the case of two spatial dimensions with $Ω=0$, we establish the global well-posedness provided that the initial data is small enough. The smallness condition on the initial data depends on the vertical viscosity and the initial radius of analyticity in the horizontal variables.

preprint2022arXiv

Well-posedness of Hibler's dynamical sea-ice model

This paper establishes the local-in-time well-posedness of solutions to an approximating system constructed by mildly regularizing the dynamical sea-ice model of {\it W.D. Hibler, Journal of Physical Oceanography, 1979}. Our choice of regularization has been carefully designed, prompted by physical considerations, to retain the original coupled hyperbolic-parabolic character of Hibler's model. Various regularized versions of this model have been used widely for the numerical simulation of the circulation and thickness of the Arctic ice cover. However, due to the singularity in the ice rheology, the notion of solutions to the original model is unclear. Instead, an approximating system, which captures current numerical study, is proposed. The well-posedness theory of such a system provides a first-step groundwork in both numerical study and future analytical study.

preprint2022arXiv

Zero Mach Number Limit of the Compressible Primitive Equations: Ill-prepared Initial Data

In the work, we consider the zero Mach number limit of compressible primitive equations in the domain $\mathbb{R}^2 \times 2\mathbb{T}$ or $\mathbb{T}^2 \times 2\mathbb{T}$. We identify the limit equations to be the primitive equations with the incompressible condition. The convergence behaviors are studied in both $\mathbb{R}^2 \times 2\mathbb{T}$ and $\mathbb{T}^2 \times 2\mathbb{T}$, respectively. This paper takes into account the high oscillating acoustic waves and is an extension of our previous work by X. Liu and E.S. Titi, Arch. Rational Mech. Anal., 238, 705-747, 2020.

preprint2021arXiv

$C^{0,α}$ boundary regularity for the pressure in weak solutions of the $2d$ Euler equations

The purpose of this note is to give a complete proof of a $C^{0,α}$ 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 $C^2$ 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.

preprint2020arXiv

Algebraic Bounds on the Rayleigh-Bénard attractor

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.

preprint2020arXiv

Asymptotic Expansions in Time for Rotating Incompressible Viscous Fluids

We study the three-dimensional Navier--Stokes equations of rotating incompressible viscous fluids with periodic boundary conditions. The asymptotic expansions, as time goes to infinity, are derived in all Gevrey spaces for any Leray-Hopf weak solutions in terms of oscillating, exponentially decaying functions. The results are established for all non-zero rotation speeds, and for both cases with and without the zero spatial average of the solutions. Our method makes use of the Poincaré waves to rewrite the equations, and then implements the Gevrey norm techniques to deal with the resulting time-dependent bi-linear form. Special solutions are also found which form infinite dimensional invariant linear manifolds.

preprint2020arXiv

Finite-time Blowup and Ill-posedness in Sobolev Spaces of the Inviscid Primitive Equations with Rotation

Large scale dynamics of the oceans and the atmosphere are governed by the primitive equations (PEs). It is well-known that the three-dimensional viscous PEs is globally well-posed in Sobolev spaces. On the other hand, the inviscid PEs without rotation is known to be ill-posed in Sobolev spaces, and its smooth solutions can form singularity in finite time. In this paper, we extend the above results in the presence of rotation. First, we construct finite-time blowup solutions to the inviscid PEs with rotation, and establish that the inviscid PEs with rotation is ill-posed in Sobolev spaces in the sense that its perturbation around a certain steady state background flow is both linearly and nonlinearly ill-posed in Sobolev spaces. Its linear instability is of the Kelvin-Helmholtz type similar to the one appears in the context of vortex sheets problem. This implies that the inviscid PEs is also linearly ill-posed in Gevrey class of order $s > 1$, and suggests that a suitable space for the well-posedness is Gevrey class of order $s = 1$, which is exactly the space of analytic functions.

preprint2020arXiv

Global well-posedness for a rapidly rotating convection model of tall columnar structure in the limit of infinite Prandtl number

We analyze a three-dimensional rapidly rotating convection model of tall columnar structure in the limit of infinite Prandtl number, i.e., when the momentum diffusivity is much more dominant than the thermal diffusivity. Consequently, the dynamics of the velocity field takes place at a much faster time scale than the temperature fluctuation, and at the limit the velocity field formally adjusts instantaneously to the thermal fluctuation. We prove the global well-posedness of weak solutions and strong solutions to this model.

preprint2020arXiv

Global Well-posedness for the Primitive Equations Coupled to Nonlinear Moisture Dynamics with Phase Changes

In this work we study the global solvability of the primitive equations for the atmosphere coupled to moisture dynamics with phase changes for warm clouds, where water is present in the form of water vapor and in the liquid state as cloud water and rain water. This moisture model contains closures for the phase changes condensation and evaporation, as well as the processes of autoconversion of cloud water into rainwater and the collection of cloud water by the falling rain droplets. It has been used by Klein and Majda in \cite{KM} and corresponds to a basic form of the bulk microphysics closure in the spirit of Kessler \cite{Ke} and Grabowski and Smolarkiewicz \cite{GS}. The moisture balances are strongly coupled to the thermodynamic equation via the latent heat associated to the phase changes. In \cite{HKLT} we assumed the velocity field to be given and proved rigorously the global existence and uniqueness of uniformly bounded solutions of the moisture balances coupled to the thermodynamic equation. In this paper we present the solvability of a full moist atmospheric flow model, where the moisture model is coupled to the primitive equations of atmospherical dynamics governing the velocity field. For the derivation of a priori estimates for the velocity field we thereby use the ideas of Cao and Titi \cite{CT}, who succeeded in proving the global solvability of the primitive equations.

preprint2020arXiv

Non-uniqueness of Weak Solutions to Hyperviscous Navier-Stokes Equations -- On Sharpness of J.-L. Lions Exponent

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 $(-Δ)^θ$, whenever the exponent $θ$ is less than J.-L. Lions' exponent $5/4$, i.e., when $θ< 5/4$.

preprint2019arXiv

On the Well-posedness of Reduced $3D$ Primitive Geostrophic Adjustment Model with Weak Dissipation

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 $α$-regularization of this model, and establish the convergence of the strong solution of the Voigt $α$-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 $α$-regularization.

preprint2019arXiv

Well-posedness of strong solutions to the anelastic equations of stratified viscous flows

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.

preprint2019arXiv

Zero Mach Number Limit of the Compressible Primitive Equations Part I: Well-prepared Initial Data

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.

preprint2016arXiv

A discrete data assimilation scheme for the solutions of the 2D Navier-Stokes equations and their statistics

We adapt a previously introduced continuous in time data assimilation (downscaling) algorithm for the 2D Navier-Stokes equations to the more realistic case when the measurements are obtained discretely in time and may be contaminated by systematic errors. Our algorithm is designed to work with a general class of observables, such as low Fourier modes and local spatial averages over finite volume elements. Under suitable conditions on the relaxation (nudging) parameter, the spatial mesh resolution and the time step between successive measurements, we obtain an asymptotic in time estimate of the difference between the approximating solution and the unknown reference solution corresponding to the measurements, in an appropriate norm, which shows exponential convergence up to a term which depends on the size of the errors. A stationary statistical analysis of our discrete data assimilation algorithm is also provided.

preprint2016arXiv

Non-viscous Regularization of the Davey-Stewartson Equations: Analysis and Modulation Theory

In the present study we are interested in the Davey-Stewartson equations (DSE) that model packets of surface and capillary-gravity waves. We focus on the elliptic-elliptic case, for which it is known that DSE may develop a finite-time singularity. We propose three systems of non-viscous regularization to the DSE in variety of parameter regimes under which the finite blow-up of solutions to the DSE occurs. We establish the global well-posedness of the regularized systems for all initial data. The regularized systems, which are inspired by the $α$-models of turbulence and therefore are called the $α$-regularized DSE, are also viewed as unbounded, singularly perturbed DSE. Therefore, we also derive reduced systems of ordinary differential equations for the $α$-regularized DSE by using the modulation theory to investigate the mechanism with which the proposed non-viscous regularization prevents the formation of the singularities in the regularized DSE. This is a follow-up of the work of Cao, Musslimani and Titi on the non-viscous $α$-regularization of the nonlinear Schrödinger equation.

preprint2016arXiv

On the Charney Conjecture of Data Assimilation Employing Temperature Measurements Alone: The Paradigm of 3D Planetary Geostrophic Model

Analyzing the validity and success of a data assimilation algorithm when some state variable observations are not available is an important problem in meteorology and engineering. We present an improved data assimilation algorithm for recovering the exact full reference solution (i.e. the velocity and temperature) of the 3D Planetary Geostrophic model, at an exponential rate in time, by employing coarse spatial mesh observations of the temperature alone. This provides, in the case of this paradigm, a rigorous justification to an earlier conjecture of Charney which states that temperature history of the atmosphere, for certain simple atmospheric models, determines all other state variables.

preprint2016arXiv

Postprocessing Galerkin method applied to a data assimilation algorithm: a uniform in time error estimate

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.

preprint2016arXiv

Recent Advances Concerning Certain Class of Geophysical Flows

This paper is devoted to reviewing several recent developments concerning certain class of geophysical models, including the primitive equations (PEs) of atmospheric and oceanic dynamics and a tropical atmosphere model. The PEs for large-scale oceanic and atmospheric dynamics are derived from the Navier-Stokes equations coupled to the heat convection by adopting the Boussinesq and hydrostatic approximations, while the tropical atmosphere model considered here is a nonlinear interaction system between the barotropic mode and the first baroclinic mode of the tropical atmosphere with moisture. We are mainly concerned with the global well-posedness of strong solutions to these systems, with full or partial viscosity, as well as certain singular perturbation small parameter limits related to these systems, including the small aspect ratio limit from the Navier-Stokes equations to the PEs, and a small relaxation-parameter in the tropical atmosphere model. These limits provide a rigorous justification to the hydrostatic balance in the PEs, and to the relaxation limit of the tropical atmosphere model, respectively. Some conditional uniqueness of weak solutions, and the global well-posedness of weak solutions with certain class of discontinuous initial data, to the PEs are also presented.

preprint2016arXiv

Strong solutions to the 3D primitive equations with only horizontal dissipation: near $H^1$ initial data

In this paper, we consider the initial-boundary value problem of the three-dimensional primitive equations for oceanic and atmospheric dynamics with only horizontal viscosity and horizontal diffusivity. We establish the local, in time, well-posedness of strong solutions, for any initial data $(v_0, T_0)\in H^1$, by using the local, in space, type energy estimate. We also establish the global well-posedness of strong solutions for this system, with any initial data $(v_0, T_0)\in H^1\cap L^\infty$, such that $\partial_zv_0\in L^m$, for some $m\in(2,\infty)$, by using the logarithmic type anisotropic Sobolev inequality and a logarithmic type Gronwall inequality. This paper improves the previous results obtained in [Cao, C.; Li, J.; Titi, E.S.: Global well-posedness of the 3D primitive equations with only horizontal viscosity and diffusivity, Comm. Pure Appl.Math., Vol. 69 (2016), 1492-1531.], where the initial data $(v_0, T_0)$ was assumed to have $H^2$ regularity.

preprint2015arXiv

A Blow-Up Criterion for the 3D Euler Equations Via the Euler-Voigt Inviscid Regularization

We propose a new blow-up criterion for the 3D Euler equations of incompressible fluid flows, based on the 3D Euler-Voigt inviscid regularization. This criterion is similar in character to a criterion proposed in a previous work by the authors, but it is stronger, and better adapted for computational tests. The 3D Euler-Voigt equations enjoy global well-posedness, and moreover are more tractable to simulate than the 3D Euler equations. A major advantage of these new criteria is that one only needs to simulate the 3D Euler-Voigt, and not the 3D Euler equations, to test the blow-up criteria, for the 3D Euler equations, computationally.

preprint2015arXiv

A Computational Study of a Data Assimilation Algorithm for the Two-dimensional Navier--Stokes Equations

We study the numerical performance of a continuous data assimilation (downscaling) algorithm, based on ideas from feedback control theory, in the context of the two-dimensional incompressible Navier--Stokes equations. Our model problem is to recover an unknown reference solution, asymptotically in time, by using continuous-in-time coarse-mesh nodal-point observational measurements of the velocity field of this reference solution (subsampling), as might be measured by an array of weather vane anemometers. Our calculations show that the required nodal observation density is remarkably less that what is suggested by the analytical study; and is in fact comparable to the {\it number of numerically determining Fourier modes}, which was reported in an earlier computational study by the authors. Thus, this method is computationally efficient and performs far better than the analytical estimates suggest.

preprint2015arXiv

A Tropical Atmosphere Model with Moisture: Global Well-posedness and Relaxation Limit

In this paper, we consider a nonlinear interaction system between the barotropic mode and the first baroclinic mode of the tropical atmosphere with moisture; that was derived in [Frierson, D.M.W.; Majda, A.J.; Pauluis, O.M.: Dynamics of precipitation fronts in the tropical atmosphere: a novel relaxation limit, Commum. Math. Sci., 2 (2004), 591-626.] We establish the global existence and uniqueness of strong solutions to this system, with initial data in $H^1$, for each fixed convective adjustment relaxation time parameter $\varepsilon>0$. Moreover, if the initial data enjoy slightly more regularity than $H^1$, then the unique strong solution depends continuously on the initial data. Furthermore, by establishing several appropriate $\varepsilon$-independent estimates, we prove that the system converges to a limiting system, as the relaxation time parameter $\varepsilon$ tends to zero, with convergence rate of the order $O(\sqrt\varepsilon)$. Moreover, the limiting system has a unique global strong solution, for any initial data in $H^1$, and such unique strong solution depends continuously on the initial data if the the initial data posses slightly more regularity than $H^1$. Notably, this solves the VISCOUS VERSION of an open problem proposed in the above mentioned paper of Frierson, Majda and Pauluis.

preprint2015arXiv

Abridged continuous data assimilation for the 2D Navier-Stokes equations utilizing measurements of only one component of the velocity field

We introduce a continuous data assimilation (downscaling) algorithm for the two-dimensional Navier-Stokes equations employing coarse mesh measurements of only one component of the velocity field. This algorithm can be implemented with a variety of finitely many observables: low Fourier modes, nodal values, finite volume averages, or finite elements. 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 asymptotically in time, to the unique exact unknown reference solution, of the 2D Navier-Stokes equations, associated with the observed (finite dimensional projection of) velocity.

preprint2015arXiv

Continuous Data Assimilation for the 2D Bénard Convection through Velocity Measurements Alone

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.

preprint2015arXiv

Continuous data assimilation for the three-dimensional Brinkman-Forchheimer-extended Darcy model

In this paper we introduce and analyze an algorithm for continuous data assimilation for a three-dimensional Brinkman-Forchheimer-extended Darcy (3D BFeD) model of porous media. This model is believed to be accurate when the flow velocity is too large for Darcy's law to be valid, and additionally the porosity is not too small. The algorithm is inspired by ideas developed for designing finite-parameters feedback control for dissipative systems. It aims to obtaining improved estimates of the state of the physical system by incorporating deterministic or noisy measurements and observations. Specifically, the algorithm involves a feedback control that nudges the large scales of the approximate solution toward those of the reference solution associated with the spatial measurements. In the first part of the paper, we present few results of existence and uniqueness of weak and strong solutions of the 3D BFeD system. The second part is devoted to the setting and convergence analysis of the data assimilation algorithm.

preprint2015arXiv

Data Assimilation algorithm for 3D Bénard convection in porous media employing only temperature measurements

In this paper we propose a continuous data assimilation (downscaling) algorithm for the Bénard convection in porous media using only coarse mesh measurements of the temperature. In this algorithm, we incorporate the observables as a feedback (nudging) term in the evolution equation of the temperature. We show that under an appropriate choice of the nudging parameter and the size of the mesh, and under the assumption that the observed data is error free, the solution of the proposed algorithm converges at an exponential rate, asymptotically in time, to the unique exact unknown reference solution of the original system, associated with the observed (finite dimensional projection of) temperature data. Moreover, we note that in the case where the observational measurements are not error free, one can estimate the error between the solution of the algorithm and the exact reference solution of the system in terms of the error in the measurements.

preprint2015arXiv

Determining form and data assimilation algorithm for weakly damped and driven Korteweg-de Vries equaton- Fourier modes case

We show that the global attractor of a weakly damped and driven Korteweg-de Vries equation (KdV) is embedded in the long-time dynamics of an ordinary differential equation called a determining form. In particular, there is a one-to-one identification of the trajectories in the global attractor of the damped and driven KdV and the steady state solutions of the determining form. Moreover, we analyze a data assimilation algorithm (down-scaling) for the weakly damped and driven KdV. We show that given a certain number of low Fourier modes of a reference solution of the KdV equation, the algorithm recovers the full reference solution at an exponential rate in time.

preprint2015arXiv

Existence and uniqueness of weak solutions to viscous primitive equations for certain class of discontinuous initial data

We establish some conditional uniqueness of weak solutions to the viscous primitive equations, and as an application, we prove the global existence and uniqueness of weak solutions, with the initial data taken as small $L^\infty$ perturbations of functions in the space $X=\left\{v\in (L^6(Ω))^2|\partial_zv\in (L^2(Ω))^2\right\}$; in particular, the initial data are allowed to be discontinuous. Our result generalizes in a uniform way the result on the uniqueness of weak solutions with continuous initial data and that of the so-called $z$-weak solutions.

preprint2015arXiv

Finite-Parameters Feedback Control for Stabilizing Damped Nonlinear Wave Equations

In this paper we introduce a finite-parameters feedback control algorithm for stabilizing solutions of various classes of damped nonlinear wave equations. Specifically, stabilization the zero steady state solution of initial boundary value problems for nonlinear weakly and strongly damped wave equations, nonlinear wave equation with nonlinear damping term and some related nonlinear wave equations, introducing a feedback control terms that employ parameters, such as, finitely many Fourier modes, finitely many volume elements and finitely many nodal observables and controllers. In addition, we also establish the stabilization of the zero steady state solution to initial boundary value problem for the damped nonlinear wave equation with a controller acting in a proper subdomain. Notably, the feedback controllers proposed here can be equally applied for stabilizing other solutions of the underlying equations.

preprint2015arXiv

Global well-posedness of strong solutions to a tropical climate model

In this paper, we consider the Cauchy problem to the TROPIC CLIMATE MODEL derived by Frierson-Majda-Pauluis in [Comm. Math. Sci, Vol. 2 (2004)] which is a coupled system of the barotropic and the first baroclinic modes of the velocity and the typical midtropospheric temperature. The system considered in this paper has viscosities in the momentum equations, but no diffusivity in the temperature equation. We establish here the global well-posedness of strong solutions to this model. In proving the global existence of strong solutions, to overcome the difficulty caused by the absence of the diffusivity in the temperature equation, we introduce a new velocity $w$ (called the pseudo baroclinic velocity), which has more regularities than the original baroclinic mode of the velocity. An auxiliary function $ϕ$, which looks like the effective viscous flux for the compressible Navier-Stokes equations, is also introduced to obtain the $L^\infty$ bound of the temperature. Regarding the uniqueness, we use the idea of performing suitable energy estimates at level one order lower than the natural basic energy estimates for the system.

preprint2015arXiv

Global well-posedness of the 2D Boussinesq equations with vertical dissipation

We prove the global well-posedness of the two-dimensional Boussinesq equations with only vertical dissipation. The initial data $(u_0,θ_0)$ are required to be only in the space $X=\{f\in L^2(\mathbb R^2)\,|\,\partial_xf\in L^2(\mathbb R^2)\}$, and thus our result generalizes that in [C. Cao, J. Wu, Global regularity for the two-dimensional anisotropic Boussinesq equations with vertical dissipation, Arch. Rational Mech. Anal., Vol. 208 (2013), 985-1004], where the initial data are assumed to be in $H^2(\mathbb R^2)$. The assumption on the initial data is at the minimal level that is required to guarantee the uniqueness of the solutions. A logarithmic type limiting Sobolev embedding inequality for the $L^\infty(\mathbb R^2)$ norm, in terms of anisotropic Sobolev norms, and a logarithmic type Gronwall inequality are established to obtain the global in time a priori estimates, which guarantee the local solution to be a global one.

preprint2015arXiv

On the backward behavior of some dissipative evolution equations

We prove that every solution of a KdV-Burgers-Sivashinsky type equation blows up in the energy space, backward in time, provided the solution does not belong to the global attractor. This is a phenomenon contrast to the backward behavior of the periodic 2D Navier-Stokes equations studied by Constantin-Foias-Kukavica-Majda [18], but analogous to the backward behavior of the Kuramoto-Sivashinsky equation discovered by Kukavica-Malcok [50]. Also we study the backward behavior of solutions to the damped driven nonlinear Schrodinger equation, the complex Ginzburg-Landau equation, and the hyperviscous Navier-Stokes equations. In addition, we provide some physical interpretation of various backward behaviors of several perturbations of the KdV equation by studying explicit cnoidal wave solutions. Furthermore, we discuss the connection between the backward behavior and the energy spectra of the solutions. The study of backward behavior of dissipative evolution equations is motivated by the investigation of the Bardos-Tartar conjecture stated in [5].

preprint2014arXiv

A determining form for the damped driven Nonlinear Schrödinger Equation- Fourier modes case

In this paper we show that the global attractor of the 1D damped, driven, nonlinear Schrödinger equation (NLS) is embedded in the long-time dynamics of a determining form. The determining form is an ordinary differential equation in a space of trajectories $X=C_b^1(\mathbb{R}, P_mH^2)$ where $P_m$ is the $L^2$-projector onto the span of the first $m$ Fourier modes. There is a one-to-one identification with the trajectories in the global attractor of the NLS and the steady states of the determining form. We also give an improved estimate for the number of the determining modes.

preprint2014arXiv

Approximation of 2D Euler Equations by the Second-Grade Fluid Equations with Dirichlet Boundary Conditions

The second-grade fluid equations are a model for viscoelastic fluids, with two parameters: $α> 0$, corresponding to the elastic response, and $ν> 0$, corresponding to viscosity. Formally setting these parameters to $0$ reduces the equations to the incompressible Euler equations of ideal fluid flow. In this article we study the limits $α, ν\to 0$ of solutions of the second-grade fluid system, in a smooth, bounded, two-dimensional domain with no-slip boundary conditions. This class of problems interpolates between the Euler-$α$ model ($ν= 0$), for which the authors recently proved convergence to the solution of the incompressible Euler equations, and the Navier-Stokes case ($α= 0$), for which the vanishing viscosity limit is an important open problem. We prove three results. First, we establish convergence of the solutions of the second-grade model to those of the Euler equations provided $ν= \mathcal{O}(α^2)$, as $α\to 0$, extending the main result in [19]. Second, we prove equivalence between convergence (of the second-grade fluid equations to the Euler equations) and vanishing of the energy dissipation in a suitably thin region near the boundary, in the asymptotic regime $ν= \mathcal{O}(α^{6/5})$, $ν/α^2 \to \infty$ as $α\to 0$. This amounts to a convergence criterion similar to the well-known Kato criterion for the vanishing viscosity limit of the Navier-Stokes equations to the Euler equations. Finally, we obtain an extension of Kato's classical criterion to the second-grade fluid model, valid if $α= \mathcal{O}(ν^{3/2})$, as $ν\to 0$. The proof of all these results relies on energy estimates and boundary correctors, following the original idea by Kato.

preprint2014arXiv

Convergence of the 2D Euler-$α$ to Euler equations in the Dirichlet case: indifference to boundary layers

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.

preprint2014arXiv

Feedback Control of Nonlinear Dissipative Systems by Finite Determining Parameters - A Reaction-diffusion Paradigm

We introduce here a simple finite-dimensional feedback control scheme for stabilizing solutions of infinite-dimensional dissipative evolution equations, such as reaction-diffusion systems, the Navier-Stokes equations and the Kuramoto-Sivashinsky equation. The designed feedback control scheme takes advantage of the fact that such systems possess finite number of determining parameters (degrees of freedom), namely, finite number of determining Fourier modes, determining nodes, and determining interpolants and projections. In particular, the feedback control scheme uses finitely many of such observables and controllers. This observation is of a particular interest since it implies that our approach has far more reaching applications, in particular, in data assimilation. Moreover, we emphasize that our scheme treats all kinds of the determining projections, as well as, the various dissipative equations with one unified approach. However, for the sake of simplicity we demonstrate our approach in this paper to a one-dimensional reaction-diffusion equation paradigm.

preprint2014arXiv

Global Regularity vs. Finite-Time Singularities: Some Paradigms on the Effect of Boundary Conditions and Certain Perturbations

In light of the question of finite-time blow-up vs. global well-posedness of solutions to problems involving nonlinear partial differential equations, we provide several cautionary examples which indicate that modifications to the boundary conditions or to the nonlinearity of the equations can effect whether the equations develop finite-time singularities. In particular, we aim to underscore the idea that in analytical and computational investigations of the blow-up of three-dimensional Euler and Navier-Stokes equations, the boundary conditions may need to be taken into greater account. We also examine a perturbation of the nonlinearity by dropping the advection term in the evolution of the derivative of the solutions to the viscous Burgers equation, which leads to the development of singularities not present in the original equation, and indicates that there is a regularizing mechanism in part of the nonlinearity. This simple analytical example corroborates recent computational observations in the singularity formation of fluid equations.

preprint2014arXiv

Global Well-posedness of the 3D Primitive Equations with Only Horizontal Viscosity and Diffusion

In this paper, we consider the initial-boundary value problem of the 3D primitive equations for planetary oceanic and atmospheric dynamics with only horizontal eddy viscosity in the horizontal momentum equations and only horizontal diffusion in the temperature equation. Global well-posedness of strong solution is established for any $H^2$ initial data. An $N$-dimensional logarithmic Sobolev embedding inequality, which bounds the $L^\infty$ norm in terms of the $L^q$ norms up to a logarithm of the $L^p$-norm, for $p>N$, of the first order derivatives, and a system version of the classic Gronwall inequality are exploited to establish the required a priori $H^2$ estimates for the global regularity.

preprint2014arXiv

Hadamard well-posedness for a hyperbolic equation of viscoelasticity with supercritical sources and damping

Presented here is a study of a viscoelastic wave equation with supercritical source and damping terms. We employ the theory of monotone operators and nonlinear semigroups, combined with energy methods to establish the existence of a unique local weak solution. In addition, it is shown that the solution depends continuously on the initial data and is global provided the damping dominates the source in an appropriate sense.

preprint2014arXiv

Inertial Manifolds for Certain Sub-Grid Scale $α$-Models of Turbulence

In this note we prove the existence of an inertial manifold, i.e., a global invariant, exponentially attracting, finite-dimensional smooth manifold, for two different sub-grid scale $α$-models of turbulence: the simplified Bardina model and the modified Leray-$α$ model, in two-dimensional space. That is, we show the existence of an exact rule that parameterizes the dynamics of small spatial scales in terms of the dynamics of the large ones. In particular, this implies that the long-time dynamics of these turbulence models is equivalent to that of a finite-dimensional system of ordinary differential equations.

preprint2014arXiv

On the uniqueness of weak solutions to the Ericksen-Leslie liquid crystal model in $\mathbb R^2$

This paper concerns the uniqueness of weak solutions to the Cauchy problem to the Ericksen-Leslie system of liquid crystal models in $\mathbb R^2$, with both general Leslie stress tensors and general Oseen-Frank density. It is shown here that such a system admits a unique weak solution provided that the Frank coefficients are close to some positive constant, which solves an interesting open problem. One of the main ideas of our proof is to perform suitable energy estimates at level one order lower than the natural basic energy estimates for the Ericksen-Leslie system.

preprint2013arXiv

A unified approach to determining forms for the 2D Navier-Stokes equations - the general interpolants case

In this paper we show that the long time dynamics (the global attractor) of the 2D Navier-Stokes equation is embedded in the long time dynamics of an ordinary differential equation, named {\it determining form}, in a space of trajectories which is isomorphic to $C^1_b(\bR; \bR^N)$, for $N$ large enough depending on the physical parameters of the Navier-Stokes equations. We present a unified approach based on interpolant operators that are induced by any of the determining parameters for the Navier-Stokes equations, namely, determining nodal values, Fourier modes, finite volume elements, finite elements, etc... There are two immediate and interesting consequences of this unified approach. The first is that the constructed determining form has a Lyapunov function, thus its solutions converge, as time goes to infinity, to the set of steady states of the determining form. The second is that these steady states of the determining form are identified, one-to-one, with the trajectories on the global attractor of the Navier-Stokes equations. It is worth adding that this unified approach is general enough that it applies, in an almost straightforward manner, to a whole class of dissipative dynamical systems.

preprint2013arXiv

Continuous Data Assimilation Using General Interpolant Observables

We present a new continuous data assimilation algorithm based on ideas that have been developed for designing finite-dimensional feedback controls for dissipative dynamical systems, in particular, in the context of the incompressible two-dimensional Navier--Stokes equations. These ideas are motivated by the fact that dissipative dynamical systems possess finite numbers of determining parameters (degrees of freedom) such as modes, nodes and local spatial averages which govern their long-term behavior. Therefore, our algorithm allows the use of any type of measurement data for which a general type of approximation interpolation operator exists. Our main result provides conditions, on the finite-dimensional spatial resolution of the collected data, sufficient to guarantee that the approximating solution, obtained by our algorithm from the measurement data, converges to the unknown reference solution over time. Our algorithm is also applicable in the context of signal synchronization in which one can recover, asymptotically in time, the solution (signal) of the underlying dissipative system that is corresponding to a continuously transmitted partial data.

preprint2013arXiv

Dissipative length scale estimates for turbulent flows - a Wiener algebra approach

In this paper, a lower bound estimate on the uniform radius of spatial analyticity is established for solutions to the incompressible, forced Navier-Stokes system on an n-torus. This estimate improves or matches previously known estimates provided that certain bounds on the initial data are satisfied. It is argued that for 2D or 3D turbulent flows, the initial data is guaranteed to satisfy these hypothesized bounds on a significant portion of the 2D global attractor or the 3D weak attractor. In these scenarios, the estimate obtained for 3D generalizes and improves upon that of [Doering-Titi], while in 2D, the estimate matches the best known one found in [Kukavica]. A key feature in the approach taken here is the choice of the Wiener algebra as the phase space, i.e., the Banach algebra of functions with absolutely convergent Fourier series, whose structure is suitable for the use of the so-called Gevrey norms.

preprint2013arXiv

Global Regularity for an Inviscid Three-dimensional Slow Limiting Ocean Dynamics Model

We establish, for smooth enough initial data, the global well-posedness (existence, uniqueness and continuous dependence on initial data) of solutions, for an inviscid three-dimensional {\it slow limiting ocean dynamics} model. This model was derived as a strong rotation limit of the rotating and stratified Boussinesg equations with periodic boundary conditions. To establish our results we utilize the tools developed for investigating the two-dimensional incompressible Euler equations and linear transport equations. Using a weaker formulation of the model we also show the global existence and uniqueness of solutions, for less regular initial data.

preprint2013arXiv

Global Well-posedness of a System of Nonlinearly Coupled KdV equations of Majda and Biello

This paper addresses the problem of global well-posedness of a coupled system of Korteweg-de Vries equations, derived by Majda and Biello in the context of nonlinear resonant interaction of Rossby waves, in a periodic setting in homogeneous Sobolev spaces $\dot H^s$, for $s\geq 0$. Our approach is based on a successive time-averaging method developed by Babin, Ilyin and Titi [1].

preprint2013arXiv

Mathematics and Turbulence: where do we stand?

This contribution covers the topics presented by the authors at the {\it ``Fundamental Problems of Turbulence, 50 Years after the Marseille Conference 1961"} meeting that took place in Marseille in 2011. It focuses on some of the mathematical approaches to fluid dynamics and turbulence. This contribution does not pretend to cover or answer, as the reader may discover, the fundamental questions in turbulence, however, it aims toward presenting some of the most recent advances in attacking these questions using rigorous mathematical tools. Moreover, we consider that the proofs of the mathematical statements (concerning, for instance, finite time regularity, weak solutions and vanishing viscosity) may contain information as relevant, to the understanding of the underlying problem, as the statements themselves.

preprint2013arXiv

On the Global Regularity of a Helical-decimated Version of the 3D Navier-Stokes Equations

We study the global regularity, for all time and all initial data in $H^{1/2}$, of a recently introduced decimated version of the incompressible 3D Navier-Stokes (dNS) equations. The model is based on a projection of the dynamical evolution of Navier-Stokes (NS) equations into the subspace where helicity (the $L^2-$scalar product of velocity and vorticity) is sign-definite. The presence of a second (beside energy) sign-definite inviscid conserved quadratic quantity, which is equivalent to the $H^{1/2}-$Sobolev norm, allows us to demonstrate global existence and uniqueness, of space-periodic solutions, together with continuity with respect to the initial conditions, for this decimated 3D model. This is achieved thanks to the establishment of two new estimates, for this 3D model, which show that the $H^{1/2}$ and the time average of the square of the $H^{3/2}$ norms of the velocity field remain finite. Such two additional bounds are known, in the spirit of the work of H. Fujita and T. Kato \cite{kato1,kato2}, to be sufficient for showing well-posedness for the 3D NS equations. Furthermore, they are directly linked to the helicity evolution for the dNS model, and therefore with a clear physical meaning and consequences.

preprint2013arXiv

On the Radius of Analyticity of Solutions to the Cubic Szegö Equation

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).

preprint2013arXiv

Persistency of Analyticity for Nonlinear Wave Equations: An Energy-like Approach

We study the persistence of the Gevrey class regularity of solutions to nonlinear wave equations with real analytic nonlinearity. Specifically, it is proven that the solution remains in a Gevrey class, with respect to some of its spatial variables, during its whole life-span, provided the initial data is from the same Gevrey class with respect to these spatial variables. In addition, for the special Gevrey class of analytic functions, we find a lower bound for the radius of the spatial analyticity of the solution that might shrink either algebraically or exponentially, in time, depending on the structure of the nonlinearity. The standard $L^2$ theory for the Gevrey class regularity is employed; we also employ energy-like methods for a generalized version of Gevrey classes based on the $\ell^1$ norm of Fourier transforms (Wiener algebra). After careful comparisons, we observe an indication that the $\ell^1$ approach provides a better lower bound for the radius of analyticity of the solutions than the $L^2$ approach. We present our results in the case of period boundary conditions, however, by employing exactly the same tools and proofs one can obtain similar results for the nonlinear wave equations and the nonlinear Schrödinger equation, with real analytic nonlinearity, in certain domains and manifolds without physical boundaries, such as the whole space $\mathbb{R}^n$, or on the sphere $\mathbb{S}^{n-1}$.

preprint2013arXiv

Planar limits of three-dimensional incompressible flows with helical symmetry

Helical symmetry is invariance under a one-dimensional group of rigid motions generated by a simultaneous rotation around a fixed axis and translation along the same axis. The key parameter in helical symmetry is the step or pitch, the magnitude of the translation after rotating one full turn around the symmetry axis. In this article we study the limits of three-dimensional helical viscous and inviscid incompressible flows in an infinite circular pipe, with respectively no-slip and no-penetration boundary conditions, as the step approaches infinity. We show that, as the step becomes large, the three-dimensional helical flow approaches a planar flow, which is governed by the so-called two-and-half Navier-Stokes and Euler equations, respectively.

preprint2013arXiv

The 3D incompressible Euler equations with a passive scalar: a road to blow-up?

The 3D incompressible Euler equations with a passive scalar $θ$ are considered in a smooth domain $Ω\subset \mathbb{R}^{3}$ with no-normal-flow boundary conditions $\bu\cdot\bhn|_{\partialΩ} = 0$. It is shown that smooth solutions blow up in a finite time if a null (zero) point develops in the vector $\bB = \nabla q\times\nablaθ$, provided $\bB$ has no null points initially\,: $\bom = \mbox{curl}\,\bu$ is the vorticity and $q = \bom\cdot\nablaθ$ is a potential vorticity. The presence of the passive scalar concentration $θ$ is an essential component of this criterion in detecting the formation of a singularity. The problem is discussed in the light of a kinematic result by Graham and Henyey (2000) on the non-existence of Clebsch potentials in the neighbourhood of null points.

preprint2012arXiv

Analysis of a Mixture Model of Tumor Growth

We study an initial-boundary value problem (IBVP) for a coupled Cahn-Hilliard-Hele-Shaw system that models tumor growth. For large initial data with finite energy, we prove global (local resp.) existence, uniqueness, higher order spatial regularity and Gevrey spatial regularity of strong solutions to the IBVP in 2D (3D resp.). Asymptotically in time, we show that the solution converges to a constant state exponentially fast as time tends to infinity under certain assumptions.

preprint2012arXiv

Finite-time Blowup for the Inviscid Primitive Equations of Oceanic and Atmospheric Dynamics

In an earlier work we have shown the global (for all initial data and all time) well-posedness of strong solutions to the three-dimensional viscous primitive equations of large scale oceanic and atmospheric dynamics. In this paper we show that for certain class of initial data the corresponding smooth solutions of the inviscid (non-viscous) primitive equations blow up in finite time. Specifically, we consider the three-dimensional inviscid primitive equations in a three-dimensional infinite horizontal channel, subject to periodic boundary conditions in the horizontal directions, and with no-normal flow boundary conditions on the solid, top and bottom, boundaries. For certain class of initial data we reduce this system into the two-dimensional system of primitive equations in an infinite horizontal strip with the same type of boundary conditions; and then show that for specific sub-class of initial data the corresponding smooth solutions of the reduced inviscid two-dimensional system develop singularities in finite time.

preprint2012arXiv

Long-time Behavior of a Two-layer Model of Baroclinic Quasi-geostrophic Turbulence

We study a viscous two-layer quasi-geostrophic beta-plane model that is forced by imposition of a spatially uniform vertical shear in the eastward (zonal) component of the layer flows, or equivalently a spatially uniform north-south temperature gradient. We prove that the model is linearly unstable, but that non-linear solutions are bounded in time by a bound which is independent of the initial data and is determined only by the physical parameters of the model. We further prove, using arguments first presented in the study of the Kuramoto-Sivashinsky equation, the existence of an absorbing ball in appropriate function spaces, and in fact the existence of a compact finite-dimensional attractor, and provide upper bounds for the fractal and Hausdorff dimensions of the attractor. Finally, we show the existence of an inertial manifold for the dynamical system generated by the model's solution operator. Our results provide rigorous justification for observations made by Panetta based on long-time numerical integrations of the model equations.

preprint2012arXiv

Navier-Stokes equations, determining forms, determining modes, inertial manifolds, dissipative dynamical systems

The determining modes for the two-dimensional incompressible Navier-Stokes equations (NSE) are shown to satisfy an ordinary differential equation of the form $dv/dt=F(v)$, in the Banach space, $X$, of all bounded continuous functions of the variable $s\in\mathbb{R}$ with values in certain finite-dimensional linear space. This new evolution ODE, named {\it determining form}, induces an infinite-dimensional dynamical system in the space $X$ which is noteworthy for two reasons. One is that $F$ is globally Lipschitz from $X$ into itself. The other is that the long-term dynamics of the determining form contains that of the NSE; the traveling wave solutions of the determining form, i.e., those of the form $v(t,s)=v_0(t+s)$, correspond exactly to initial data $v_0$ that are projections of solutions of the global attractor of the NSE onto the determining modes. The determining form is also shown to be dissipative; an estimate for the radius of an absorbing ball is derived in terms of the number of determining modes and the Grashof number (a dimensionless physical parameter). Finally, a unified approach is outlined for an ODE satisfied by a variety of other determining parameters such as nodal values, finite volumes, and finite elements.

preprint2012arXiv

Stability of Two-dimensional Viscous Incompressible Flows Under Three-dimensional Perturbations and Inviscid Symmetry Breaking

In this article we consider weak solutions of the three-dimensional incompressible fluid flow equations with initial data admitting a one-dimensional symmetry group. We examine both the viscous and inviscid cases. For the case of viscous flows, we prove that Leray-Hopf weak solutions of the three-dimensional Navier-Stokes equations preserve initially imposed symmetry and that such symmetric flows are stable under general three-dimensional perturbations, globally in time. We work in three different contexts: two-and-a-half-dimensional, helical and axi-symmetric flows. In the inviscid case, we observe that, as a consequence of recent work by De Lellis and Székelyhidi, there are genuinely three-dimensional weak solutions of the Euler equations with two-dimensional initial data. We also present two partial results where restrictions on the set of initial data, and on the set of admissible solutions rule out spontaneous symmetry breaking; one is due to P.-L. Lions and the other is a consequence of our viscous stability result.

preprint2012arXiv

Vanishing viscosity as a selection principle for the Euler equations: The case of 3D shear flow

We show that for a certain family of initial data, there exist non-unique weak solutions to the 3D incompressible Euler equations satisfying the weak energy inequality, whereas the weak limit of every sequence of Leray-Hopf weak solutions for the Navier-Stokes equations, with the same initial data, and as the viscosity tends to zero, is uniquely determined and equals the shear flow solution of the Euler equations. This simple example suggests that, also in more general situations, the vanishing viscosity limit of the Navier-Stokes equations could serve as a uniqueness criterion for weak solutions of the Euler equations.

preprint2011arXiv

Global Well-posedness of an Inviscid Three-dimensional Pseudo-Hasegawa-Mima Model

The three-dimensional inviscid Hasegawa-Mima model is one of the fundamental models that describe plasma turbulence. The model also appears as a simplified reduced Rayleigh-Bénard convection model. The mathematical analysis the Hasegawa-Mima equation is challenging due to the absence of any smoothing viscous terms, as well as to the presence of an analogue of the vortex stretching terms. In this paper, we introduce and study a model which is inspired by the inviscid Hasegawa-Mima model, which we call a pseudo-Hasegawa-Mima model. The introduced model is easier to investigate analytically than the original inviscid Hasegawa-Mima model, as it has a nicer mathematical structure. The resemblance between this model and the Euler equations of inviscid incompressible fluids inspired us to adapt the techniques and ideas introduced for the two-dimensional and the three-dimensional Euler equations to prove the global existence and uniqueness of solutions for our model. Moreover, we prove the continuous dependence on initial data of solutions for the pseudo-Hasegawa-Mima model. These are the first results on existence and uniqueness of solutions for a model that is related to the three-dimensional inviscid Hasegawa-Mima equations.

preprint2011arXiv

Higher-Order Global Regularity of an Inviscid Voigt-Regularization of the Three-Dimensional Inviscid Resistive Magnetohydrodynamic Equations

We prove existence, uniqueness, and higher-order global regularity of strong solutions to a particular Voigt-regularization of the three-dimensional inviscid resistive Magnetohydrodynamic (MHD) equations. Specifically, the coupling of a resistive magnetic field to the Euler-Voigt model is introduced to form an inviscid regularization of the inviscid resistive MHD system. The results hold in both the whole space $\nR^3$ and in the context of periodic boundary conditions. Weak solutions for this regularized model are also considered, and proven to exist globally in time, but the question of uniqueness for weak solutions is still open. Since the main purpose of this line of research is to introduce a reliable and stable inviscid numerical regularization of the underlying model we, in particular, show that the solutions of the Voigt regularized system converge, as the regularization parameter $α\maps0$, to strong solutions of the original inviscid resistive MHD, on the corresponding time interval of existence of the latter. Moreover, we also establish a new criterion for blow-up of solutions to the original MHD system inspired by this Voigt regularization. This type of regularization, and the corresponding results, are valid for, and can also be applied to, a wide class of hydrodynamic models.

preprint2010arXiv

Discrete Data Assimilation in the Lorenz and 2D Navier--Stokes Equations

Consider a continuous dynamical system for which partial information about its current state is observed at a sequence of discrete times. Discrete data assimilation inserts these observational measurements of the reference dynamical system into an approximate solution by means of an impulsive forcing. In this way the approximating solution is coupled to the reference solution at a discrete sequence of points in time. This paper studies discrete data assimilation for the Lorenz equations and the incompressible two-dimensional Navier--Stokes equations. In both cases we obtain bounds on the time interval h between subsequent observations which guarantee the convergence of the approximating solution obtained by discrete data assimilation to the reference solution.

preprint2010arXiv

Global regularity criterion for the 3D Navier-Stokes equations involving one entry of the velocity gradient tensor

In this paper we provide a sufficient condition, in terms of only one of the nine entries of the gradient tensor, i.e., the Jacobian matrix of the velocity vector field, for the global regularity of strong solutions to the three-dimensional Navier-Stokes equations in the whole space, as well as for the case of periodic boundary conditions.

preprint2010arXiv

Global Well-posedness for The 2D Boussinesq System Without Heat Diffusion and With Either Anisotropic Viscosity or Inviscid Voigt-$α$ Regularization

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.

preprint2010arXiv

Global Well-posedness of the 3D Primitive Equations With Partial Vertical Turbulence Mixing Heat Diffusion

The three--dimensional incompressible viscous Boussinesq equations, under the assumption of hydrostatic balance, govern the large scale dynamics of atmospheric and oceanic motion, and are commonly called the primitive equations. To overcome the turbulence mixing a partial vertical diffusion is usually added to the temperature advection (or density stratification) equation. In this paper we prove the global regularity of strong solutions to this model in a three-dimensional infinite horizontal channel, subject to periodic boundary conditions in the horizontal directions, and with no-penetration and stress-free boundary conditions on the solid, top and bottom, boundaries. Specifically, we show that short time strong solutions to the above problem exist globally in time, and that they depend continuously on the initial data.

preprint2010arXiv

On the Higher-Order Global Regularity of the Inviscid Voigt-Regularization of Three-Dimensional Hydrodynamic Models

We prove higher-order and a Gevrey class (spatial analytic) regularity of solutions to the Euler-Voigt inviscid $α$-regularization of the three-dimensional Euler equations of ideal incompressible fluids. Moreover, we establish the convergence of strong solutions of the Euler-Voigt model to the corresponding solution of the three-dimensional Euler equations for inviscid flow on the interval of existence of the latter. Furthermore, we derive a criterion for finite-time blow-up of the Euler equations based on this inviscid regularization. The coupling of a magnetic field to the Euler-Voigt model is introduced to form an inviscid regularization of the inviscid irresistive magneto-hydrodynamic (MHD) system. Global regularity of the regularized MHD system is also established.

preprint2010arXiv

On the regularization mechanism for the periodic Korteweg-de Vries equation

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]$.

preprint2010arXiv

Regularity "in Large" for the 3D Salmon's Planetary Geostrophic Model of Ocean Dynamics

It is well known, by now, that the three-dimensional non-viscous planetary geostrophic model, with vertical hydrostatic balance and horizontal Rayleigh friction, coupled to the heat diffusion and transport, is mathematically ill-posed. This is because the no-normal flow physical boundary condition implicitly produces an additional boundary condition for the temperature at the literal boundary. This additional boundary condition is different, because of the Coriolis forcing term, than the no heat flux physical boundary condition. Consequently, the second order parabolic heat equation is over determined with two different boundary conditions. In a previous work we proposed one remedy to this problem by introducing a fourth-order artificial hyper-diffusion to the heat transport equation and proved global regularity for the proposed model. Another remedy for this problem was suggested by R. Salmon by introducing an additional Rayleigh-like friction term for the vertical component of the velocity in the hydrostatic balance equation. In this paper we prove the global, for all time and all initial data, well-posedness of strong solutions to the three-dimensional Salmon's planetary geostrophic model of ocean dynamics. That is, we show global existence, uniqueness and continuous dependence of the strong solutions on initial data for this model.

preprint2009arXiv

On the convergence rate of the Euler-$α$, an inviscid second-grade complex fluid, model to the Euler equations

We study the convergence rate of the solutions of the incompressible Euler-$α$, an inviscid second-grade complex fluid, equations to the corresponding solutions of the Euler equations, as the regularization parameter $α$ approaches zero. First we show the convergence in $H^{s}$, $s>n/2+1$, in the whole space, and that the smooth Euler-$α$ solutions exist at least as long as the corresponding solution of the Euler equations. Next we estimate the convergence rate for two-dimensional vortex patch with smooth boundaries.

preprint2009arXiv

Stochastic attractors for shell phenomenological models of turbulence

Recently, it has been proposed that the Navier-Stokes equations and a relevant linear advection model have the same long-time statistical properties, in particular, they have the same scaling exponents of their structure functions. This assertion has been investigate rigorously in the context of certain nonlinear deterministic phenomenological shell model, the Sabra shell model, of turbulence and its corresponding linear advection counterpart model. This relationship has been established through a "homotopy-like" coefficient $λ$ which bridges continuously between the two systems. That is, for $λ=1$ one obtains the full nonlinear model, and the corresponding linear advection model is achieved for $λ=0$. In this paper, we investigate the validity of this assertion for certain stochastic phenomenological shell models of turbulence driven by an additive noise. We prove the continuous dependence of the solutions with respect to the parameter $λ$. Moreover, we show the existence of a finite-dimensional random attractor for each value of $λ$ and establish the upper semicontinuity property of this random attractors, with respect to the parameter $λ$. This property is proved by a pathwise argument. Our study aims toward the development of basic results and techniques that may contribute to the understanding of the relation between the long-time statistical properties of the nonlinear and linear models.

preprint2006arXiv

A Note on the Regularity of Inviscid Shell Model of Turbulence

In this paper we continue the analytical study of the sabra shell model of energy turbulent cascade initiated in \cite{CLT05}. We prove the global existence of weak solutions of the inviscid sabra shell model, and show that these solutions are unique for some short interval of time. In addition, we prove that the solutions conserve the energy, provided that the components of the solution satisfy $|{u_n}| \le C k_n^{-1/3} (\sqrt{n} \log(n+1))^{-1}$, for some positive absolute constant $C$, which is the analogue of the Onsager's conjecture for the Euler's equations. Moreover, we give a Beal-Kato-Majda type criterion for the blow-up of solutions of the inviscid sabra shell model and show the global regularity of the solutions in the ``two-dimensional'' parameters regime.