Search

We consider a modification of the three-dimensional Navier--Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as $\ue ^{|k|/\kd}$ at high wavenumbers $|k|$. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than $\ue ^{-C(k/\kd) \ln (|k|/\kd)}$ for any $C<1/(2\ln 2)$. The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with $C= C_\star =1/\ln2$. The same behavior with a universal constant $C_\star$ is conjectured for the Navier--Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier--Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed., Comment: 29 pages, 3 figures, Comm. Math. Phys., in press

The nonlinear Schr\"{o}dinger-Helmholtz (SH) equation in $N$ space dimensions with $2\sigma$ nonlinear power was proposed as a regularization of the classical nonlinear Schr\"{o}dinger (NLS) equation. It was shown that the SH equation has a larger regime ($1\le\sigma<\frac{4}{N}$) of global existence and uniqueness of solutions compared to that of the classical NLS ($0<\sigma<\frac{2}{N}$). In the limiting case where the Schr\"{o}dinger-Helmholtz equation is viewed as a perturbed system of the classical NLS equation, we apply modulation theory to the classical critical case ($\sigma=1,\:N=2$) and show that the regularization prevents the formation of singularities of the NLS equation. Our theoretical results are supported by numerical simulations, Comment: new article, 20 pages, 10 figures, to appear in Numerical Functional Analysis and Optimization

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

We prove uniqueness and existence of the weak solutions of Euler equations with helical symmetry, with initial vorticity in $L^{\infty}$ under "no vorticity stretching" geometric constraint. Our article follows the argument of the seminal work of Yudovich. We adjust the argument to resolve the difficulties which are specific to the helical symmetry.

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

We study a class of quadratic, infinite-dimensional dynamical systems, inspired by models for viscoelastic fluids. We prove that these equations define a semi-flow on the cone of positive, essentially bounded functions. As time tends to infinity, the solutions tend to an equilibrium manifold in the $L^2$-norm. Convergence to a particular function on the equilibrium manifold is only proved under additional assumptions. We discuss several possible generalizations.

We present an alpha-regularization of the Birkhoff-Rott equation, induced by the two-dimensional Euler-alpha equations, for the vortex sheet dynamics. We show that initially smooth self-avoiding vortex sheet remains smooth for all times under the alpha-regularized dynamics, provided the initial density of vorticity is an integrable function over the curve with respect to the arc-length measure., Comment: 7 pages, no figures, submitted to the proceedings of the Conference EE250, corrected typos, references added

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

We derive a criterion for the breakdown of solutions to the Oldroyd-B model in $\R^3$ in the limit of zero Reynolds number (creeping flow). If the initial stress field is in the Sobolev space $H^m$, $m> 5/2$, then either a unique solution exists within this space indefinitely, or, at the time where the solution breaks down, the time integral of the $L^\infty$-norm of the stress tensor must diverge. This result is analogous to the celebrated Beale-Kato-Majda breakdown criterion for the inviscid Eluer equations of incompressible fluids.

A regularized $\alpha-$system of the Nonlinear Schr\"{o}dinger Equation (NLS) with $2\sigma$ nonlinear power in dimension $N$ is studied. We prove existence and uniqueness of local solution in the case $1 \le \sigma <\frac{4}{N-2}$ and existence and uniqueness of global solution in the case $1 \le \sigma < \frac{4}{N}$. When $\alpha \to 0^+$, this regularized system will converge to the classical NLS in the appropriate range. In particular, the purpose of this numerical regularization is to shed light on the profile of the blow up solutions of the original Nonlinear Schr\"{o}dinger Equation in the range $\frac{2}{N}\le \sigma <\frac{4}{N}$, and in particular for the critical case $\sigma = \frac{2}{N}$.

We investigate the long-term dynamics of the three-dimensional Navier-Stokes-Voight model of viscoelastic incompressible fluid. Specifically, we derive upper bounds for the number of determining modes for the 3D Navier-Stokes-Voight equations and for the dimension of a global attractor of a semigroup generated by these equations. Viewed from the numerical analysis point of view we consider the Navier-Stokes-Voight model as a non-viscous (inviscid) regularization of the three-dimensional Navier-Stokes equations. Furthermore, we also show that the weak solutions of the Navier-Stokes-Voight equations converge, in the appropriate norm, to the weak solutions of the inviscid simplified Bardina model, as the viscosity coefficient $\nu \to 0$.

This article is a survey concerning the state-of-the-art mathematical theory of the Euler equations of incompressible homogenous ideal fluid. Emphasis is put on the different types of emerging instability, and how they may be related to the description of turbulence.

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

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

We obtain a logarithmically sharp estimate for the space-analyticity radius of the solutions of the damped-driven 2D Navier-Stokes equations with periodic boundary conditions and relate this to the small scales in this system. This system is inspired by the Stommel--Charney barotropic ocean circulation model.

In a recent paper it was proposed that for some nonlinear shell models of turbulence one can construct a linear advection model for an auxiliary field such that the scaling exponents of all the structure functions of the linear and nonlinear fields coincide. The argument depended on an assumption of continuity of the solutions as a function of a parameter. The aim of this paper is to provide a rigorous proof for the validity of the assumption. In addition we clarify here when the swap of a nonlinear model by a linear one will not work., Comment: 7 pages, 4 figures, submitted to Nonlinearity

In this work we derive a lower bounds for the Hausdorff and fractal dimensions of the global attractor of the Sabra shell model of turbulence in different regimes of parameters. We show that for a particular choice of the forcing and for sufficiently small viscosity term $\nu$, the Sabra shell model has a global attractor of large Hausdorff and fractal dimensions proportional to $\log_\lambda \nu^{-1}$ for all values of the governing parameter $\epsilon$, except for $\epsilon=1$. The obtained lower bounds are sharp, matching the upper bounds for the dimension of the global attractor obtained in our previous work. Moreover, we show different scenarios of the transition to chaos for different parameters regime and for specific forcing. In the ``three-dimensional'' regime of parameters this scenario changes when the parameter $\epsilon$ becomes sufficiently close to 0 or to 1. We also show that in the ``two-dimensional'' regime of parameters for a certain non-zero forcing term the long-time dynamics of the model becomes trivial for any value of the viscosity.

In this paper we consider the three-dimensional Navier-Stokes equations in infinite channel. We provide a regularity criterion for solutions of the three-dimensional Navier-Stokes equations in terms of the vertical component of the velocity field.

In this paper we consider the r\^ole that numerical computations -- in particular Galerkin approximations -- can play in problems modelled by the 3d Navier-Stokes equations, for which no rigorous proof of the existence of unique solutions is currently available. We prove a robustness theorem for strong solutions, from which we derive an {\it a posteriori} check that can be applied to a numerical solution to guarantee the existence of a strong solution of the corresponding exact problem. We then consider Galerkin approximations, and show that {\it if} a strong solution exists the Galerkin approximations will converge to it; thus if one is prepared to assume that the Navier-Stokes equations are regular one can justify this particular numerical method rigorously. Combining these two results we show that if a strong solution of the exact problem exists then this can be verified numerically using an algorithm that can be guaranteed to terminate in a finite time. We thus introduce the possibility of rigorous computations of the solutions of the 3d Navier-Stokes equations (despite the lack of rigorous existence and uniqueness results), and demonstrate that numerical investigation can be used to rule out the occurrence of possible singularities in particular examples.

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.

In this paper we present an analytical study of a subgrid scale turbulence model of the three-dimensional magnetohydrodynamic (MHD) equations, inspired by the Navier-Stokes-alpha (also known as the viscous Camassa-Holm equations or the Lagrangian-averaged Navier-Stokes-alpha model). Specifically, we show the global well-posedness and regularity of solutions of a certain MHD-alpha model (which is a particular case of the Lagrangian averaged magnetohydrodynamic-alpha model without enhancing the dissipation for the magnetic field). We also introduce other subgrid scale turbulence models, inspired by the Leray-alpha and the modified Leray-alpha models of turbulence. Finally, we discuss the relation of the MHD-alpha model to the MHD equations by proving a convergence theorem, that is, as the length scale alpha tends to zero, a subsequence of solutions of the MHD-alpha equations converges to a certain solution (a Leray-Hopf solution) of the three-dimensional MHD equations., Comment: 26 pages, no figures, will appear in Journal of Math Physics; corrected typos, updated references

We show that the only locally integrable stationary solutions to the integrated Kuramoto-Sivashinsky equation in $R$ and $R^2$ are the trivial constant solutions. We extend our technique and prove similar results to other nonlinear elliptic problems in $\R^N$., Comment: 10 pages, Journal Article

In this paper we study analytically the viscous `sabra' shell model of energy turbulent cascade. We prove the global regularity of solutions and show that the shell model has finitely many asymptotic degrees of freedom, specifically: a finite dimensional global attractor and globally invariant inertial manifolds. Moreover, we establish the existence of exponentially decaying energy dissipation range for the sufficiently smooth forcing.

We derive upper bounds for the number of asymptotic degrees (determining modes and nodes) of freedom for the two-dimensional Navier--Stokes system and Navier-Stokes system with damping. In the first case we obtain the previously known estimates in an explicit form, which are larger than the fractal dimension of the global attractor. However, for the Navier--Stokes system with damping our estimates for the number of the determining modes and nodes are comparable to the sharp estimates for the fractal dimension of the global attractor. Our investigation of the damped-driven 2D Navier--Stokes system is inspired by the Stommel--Charney barotropic model of ocean circulation where the damping represents the Rayleigh friction. We remark that our results equally apply to the Stommel--Charney model.

In this paper we prove the global existence and uniqueness (regularity) of strong solutions to the three-dimensional viscous primitive equations, which model large scale ocean and atmosphere dynamics.

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

We propose the viscous Camassa-Holm equations as a closure approximation for the Reynolds-averaged equations of the incompressible Navier-Stokes fluid. This approximation is tested on turbulent channel flows with steady mean. Analytical solutions for the mean velocity and the Reynolds shear stress across the entire channel are obtained, showing good agreement with experimental measurements and direct numerical simulations. As Reynolds number varies, these analytical mean velocity profiles form a family of curves whose envelopes are shown to have either power law, or logarithmic behavior, depending on the choice of drag law., Comment: 4 pages, 4 figs, revtex