Search

We prove finite-time vorticity blowup for smooth solutions of the 2D compressible Euler equations with smooth, localized, and non-vacuous initial data. The vorticity blowup occurs at the time of the first singularity, and is accompanied by an axisymmetric implosion in which the swirl velocity enjoys full stability, as opposed to finite co-dimension stability.

We construct a fundamental piece of the boundary of the maximal globally hyperbolic development (MGHD) of Cauchy data for the multi-dimensional compressible Euler equations, which is necessary for the local shock development problem. For an open set of compressive and generic $H^7$ initial data, we construct unique $H^7$ solutions to the Euler equations in the maximal spacetime region below a given time-slice, beyond the time of the first singularity; at any point in this spacetime, the solution can be smoothly and uniquely computed by tracing both the fast and slow acoustic characteristic surfaces backward-in-time, until reaching the Cauchy data prescribed along the initial time-slice. The future temporal boundary of this spacetime region is a singular hypersurface, containing the union of three sets: first, a co-dimension-$2$ surface of ``first singularities'' called the pre-shock; second, a downstream hypersurface called the singular set emanating from the pre-shock, on which the Euler solution experiences a continuum of gradient catastrophes; third, an upstream hypersurface consisting of a Cauchy horizon emanating from the pre-shock, which the Euler solution cannot reach. We develop a new geometric framework for the description of the acoustic characteristic surfaces which is based on the Arbitrary Lagrangian Eulerian (ALE) framework, and combine this with a new type of differentiated Riemann variables which are linear combinations of gradients of velocity, sound speed, and the curvature of the fast acoustic characteristic surfaces. With these new variables, we establish uniform $H^7$ Sobolev bounds for solutions to the Euler equations without derivative loss and with optimal regularity., Comment: 218 pages, 20 figures, minor edits

For every $\alpha < 1/3$, we construct an explicit divergence-free vector field $\mathbf{b}(t,x)$ which is periodic in space and time and belongs to $C^0_t C^{\alpha}_x \cap C^{\alpha}_t C^0_x$ such that the corresponding scalar advection-diffusion equation $$\partial_t \theta^\kappa + \mathbf{b} \cdot \nabla \theta^\kappa - \kappa \Delta \theta^\kappa = 0$$ exhibits anomalous dissipation of scalar variance for arbitrary $H^1$ initial data: $$\limsup_{\kappa \to 0} \int_0^{1} \int_{\mathbb{T}^d} \kappa | \nabla \theta^\kappa (t,x) |^2 \,dx\,dt >0.$$ The vector field is deterministic and has a fractal structure, with periodic shear flows alternating in time between different directions serving as the base fractal. These shear flows are repeatedly inserted at infinitely many scales in suitable Lagrangian coordinates. Using an argument based on ideas from quantitative homogenization, the corresponding advection-diffusion equation with small $\kappa$ is progressively renormalized, one scale at a time, starting from the (very small) length scale determined by the molecular diffusivity up to the macroscopic (unit) scale. At each renormalization step, the effective diffusivity is enhanced by the influence of advection on that scale. By iterating this procedure across many scales, the effective diffusivity on the macroscopic scale is shown to be of order one., Comment: 111 pages. Minor remark added. Announcement at: https://www.scottnarmstrong.com/2023/05/anomalous-diffusion-fractal-homogenization/

From an open set of initial data, we construct a family of classical solutions to the 1D nonisentropic compressible Euler equations which form $C^{0,\nu}$ cusps as a first singularity, for any $\nu \in [1/2,1)$. For this range of $\nu$, this is the first result demonstrating the stable formation of such $C^{0,\nu}$ cusp-type singularities, also known as pre-shocks. The proof uses a new formulation of the differentiated Euler equations along the fast acoustic characteristic, and relies on a novel set of $L^p$ energy estimates for all $1 < p < \infty$, which may be of independent interest.

We provide a detailed analysis of the shock formation process for the non-isentropic 2d Euler equations in azimuthal symmetry. We prove that from an open set of smooth and generic initial data, solutions of Euler form a first singularity or gradient blow-up or shock. This first singularity is termed a H\"{o}lder $C^{\frac{1}{3}}$ pre-shock, and our analysis provides the first detailed description of this cusp solution. The novelty of this work relative to [Buckmaster-Drivas-Shkoller-Vicol, 2022] is that we herein consider a much larger class of initial data, allow for a non-constant initial entropy, allow for a non-trivial sub-dominant Riemann variable, and introduce a host of new identities to avoid apparent derivative loss due to entropy gradients. The method of proof is also new and robust, exploring the transversality of the three different characteristic families to transform space derivatives into time derivatives. Our main result provides a fractional series expansion of the Euler solution about the pre-shock, whose coefficients are computed from the initial data.

For any regularity exponent $\beta<\frac 12$, we construct non-conservative weak solutions to the 3D incompressible Euler equations in the class $C^0_t (H^{\beta} \cap L^{\frac{1}{(1-2\beta)}})$. By interpolation, such solutions belong to $C^0_tB^{s}_{3,\infty}$ for $s$ approaching $\frac 13$ as $\beta$ approaches $\frac 12$. Hence this result provides a new proof of the flexible side of the $L^3$-based Onsager conjecture. Of equal importance is that the intermittent nature of our solutions matches that of turbulent flows, which are observed to possess an $L^2$-based regularity index exceeding $\frac 13$. Thus our result does not imply, and is not implied by, the work of Isett [A proof of Onsager's conjecture, Annals of Mathematics, 188(3):871, 2018], who gave a proof of the H\"older-based Onsager conjecture. Our proof builds on the authors' previous joint work with Buckmaster et al. (Intermittent convex integration for the 3D Euler equations: (AMS-217), Princeton University Press, 2023), in which an intermittent convex integration scheme is developed for the 3D incompressible Euler equations. We employ a scheme with higher-order Reynolds stresses, which are corrected via a combinatorial placement of intermittent pipe flows of optimal relative intermittency., Comment: 54 pages, no figures, published version available at https://link.springer.com/article/10.1007/s00222-023-01185-6. arXiv admin note: text overlap with arXiv:2101.09278

We address the problem of controllability of the MHD system in a rectangular domain with a control prescribed on the side boundary. We identify a necessary and sufficient condition on the data to be null controllable, i.e., can be driven to the zero state. We also show that the validity of this condition allows the states to be stirred to each other. If the condition is not satisfied, one can move from one state to another with the help of a simple shear external magnetic force., Comment: 14 pages

We establish the validity of the Euler$+$Prandtl approximation for solutions of the Navier-Stokes equations in the half plane with Dirichlet boundary conditions, in the vanishing viscosity limit, for initial data which are analytic only near the boundary, and Sobolev smooth away from the boundary. Our proof does not require higher order correctors, and works directly by estimating an $L^{1}$-type norm for the vorticity of the error term in the expansion Navier-Stokes$-($Euler$+$Prandtl$)$. An important ingredient in the proof is the propagation of local analyticity for the Euler equation, a result of independent interest.

A fundamental question in fluid dynamics concerns the formation of discontinuous shock waves from smooth initial data. We prove that from smooth initial data, smooth solutions to the 2d Euler equations in azimuthal symmetry form a first singularity, the so-called $C^{\frac{1}{3}} $ pre-shock. The solution in the vicinity of this pre-shock is shown to have a fractional series expansion with coefficients computed from the data. Using this precise description of the pre-shock, we prove that a discontinuous shock instantaneously develops after the pre-shock. This regular shock solution is shown to be unique in a class of entropy solutions with azimuthal symmetry and regularity determined by the pre-shock expansion. Simultaneous to the development of the shock front, two other characteristic surfaces of cusp-type singularities emerge from the pre-shock. These surfaces have been termed weak discontinuities by Landau & Lifschitz [Chapter IX, \S 96], who conjectured some type of singular behavior of derivatives along such surfaces. We prove that along the slowest surface, all fluid variables except the entropy have $C^{1, {\frac{1}{2}} }$ one-sided cusps from the shock side, and that the normal velocity is decreasing in the direction of its motion; we thus term this surface a weak rarefaction wave. Along the surface moving with the fluid velocity, density and entropy form $C^{1, {\frac{1}{2}} }$ one-sided cusps while the pressure and normal velocity remain $C^2$; as such, we term this surface a weak contact discontinuity., Comment: 150 pages, 15 figures, typos corrected

We investigate the stability properties for a family of equations introduced by Moffatt to model magnetic relaxation. These models preserve the topology of magnetic streamlines, contain a cubic nonlinearity, and yet have a favorable $L^2$ energy structure. We consider the local and global in time well-posedness of these models and establish a difference between the behavior as $t\to \infty$ with respect to weak and strong norms., Comment: 24 pages

For any positive regularity parameter $\beta < \frac 12$, we construct non-conservative weak solutions of the 3D incompressible Euler equations which lie in $H^{\beta}$ uniformly in time. In particular, we construct solutions which have an $L^2$-based regularity index \emph{strictly larger} than $\frac 13$, thus deviating from the $H^{\frac{1}{3}}$-regularity corresponding to the Kolmogorov-Obhukov $\frac 53$ power spectrum in the inertial range., Comment: 155 pages, 11 figures, typos corrected

We analyze the shock formation process for the 3d non-isentropic Euler equations with the ideal gas law, in which sounds waves interact with entropy waves to produce vorticity. Building on our theory for isentropic flows in [3,4], we give a constructive proof of shock formation from smooth initial data. Specifically, we prove that there exist smooth solutions to the non-isentropic Euler equations which form a generic stable shock with explicitly computable blowup time, location, and direction. This is achieved by establishing the asymptotic stability of a generic shock profile in modulated self-similar variables, controlling the interaction of wave families via: (i) pointwise bounds along Lagrangian trajectories, (ii) geometric vorticity structure, and (iii) high-order energy estimates in Sobolev spaces., Comment: 87 Pages, this paper builds on methods developed in arXiv:1912.04429

We consider the 3D isentropic compressible Euler equations with the ideal gas law. We provide a constructive proof of shock formation from smooth initial datum of finite energy, with no vacuum regions, with nontrivial vorticity present at the shock, and under no symmetry assumptions. We prove that for an open set of Sobolev-class initial data which are a small $L^ \infty $ perturbation of a constant state, there exist smooth solutions to the Euler equations which form a generic stable shock in finite time. The blow up time and location can be explicitly computed, and solutions at the blow up time are smooth except for a single point, where they are of cusp-type with H\"{o}lder $C^ {\frac{1}{3}}$ regularity. Our proof is based on the use of modulated self-similar variables that are used to enforce a number of constraints on the blow up profile, necessary to establish the stability in self-similar variables of the generic shock profile., Comment: 94 pages, 1 figure, minor typos corrected, description of asymptotic profile added

We construct weak solutions to the ideal magneto-hydrodynamic (MHD) equations which have finite total energy, and whose magnetic helicity is not a constant function of time. In view of Taylor's conjecture, this proves that there exist finite energy weak solutions to ideal MHD which cannot be attained in the infinite conductivity and zero viscosity limit. Our proof is based on a Nash-type convex integration scheme with intermittent building blocks adapted to the geometry of the MHD system., Comment: 31 pages

We consider the 2D isentropic compressible Euler equations, with pressure law $p(\rho) = (\sfrac{1}{\gamma}) \rho^\gamma$, with $\gamma >1$. We provide an elementary constructive proof of shock formation from smooth initial datum of finite energy, with no vacuum regions, and with {nontrivial vorticity}. We prove that for initial data which has minimum slope $- {\sfrac{1}{ \eps}}$, for $ \eps>0$ taken sufficiently small relative to the $\OO(1)$ amplitude, there exist smooth solutions to the Euler equations which form a shock in time $\OO(\eps)$. The blowup time and location can be explicitly computed and solutions at the blowup time are of cusp-type, with H\"{o}lder $C^ {\sfrac{1}{3}}$ regularity. Our objective is the construction of solutions with inherent $\OO(1)$ vorticity at the shock. As such, rather than perturbing from an irrotational regime, we instead construct solutions with dynamics dominated by purely azimuthal wave motion. We consider homogenous solutions to the Euler equations and use Riemann-type variables to obtain a system of forced transport equations. Using a transformation to modulated self-similar variables and pointwise estimates for the ensuing system of transport equations, we show the global stability, in self-similar time, of a smooth blowup profile., Comment: 39 Pages

The Triple Deck model is a classical high order boundary layer model that has been proposed to describe flow regimes where the Prandtl theory is expected to fail. At first sight the model appears to lose two derivatives through the pressure-displacement relation which links pressure to the tangential slip. In order to overcome this, we split the Triple Deck system into two coupled equations: a Prandtl type system on $\mathbb{H}$ and a Benjamin-Ono type equation on $\mathbb{R}$. This splitting enables us to extract a crucial leading order cancellation at the top of the lower deck. We develop a functional framework to subsequently extend this cancellation into the interior of the lower deck, which enables us to prove the local well-posedness of the model in tangentially real analytic spaces.

We address the inviscid limit for the Navier-Stokes equations in a half space, with initial datum that is analytic only close to the boundary of the domain, and has finite Sobolev regularity in the complement. We prove that for such data the solution of the Navier-Stokes equations converges in the vanishing viscosity limit to the solution of the Euler equation, on a constant time interval., Comment: 36 pages, 1 figure

We describe a method to construct smooth and compactly supported solutions of 3D incompressible Euler equations and related models. The method is based on localizable Grad-Shafranov equations and is inspired by the recent result \cite{gav}.

In this review article we discuss a number of recent results concerning wild weak solutions of the incompressible Euler and Navier-Stokes equations. These results build on the groundbreaking works of De Lellis and Sz\'ekelyhidi Jr., who extended Nash's fundamental ideas on $C^1$ flexible isometric embeddings, into the realm of fluid dynamics. These techniques, which go under the umbrella name convex integration, have fundamental analogies the phenomenological theories of hydrodynamic turbulence. Mathematical problems arising in turbulence (such as the Onsager conjecture) have not only sparked new interest in convex integration, but certain experimentally observed features of turbulent flows (such as intermittency) have also informed new convex integration constructions. First, we give an elementary construction of nonconservative $C^{0+}_{x,t}$ weak solutions of the Euler equations, first proven by De Lellis-Sz\'ekelyhidi Jr.. Second, we present Isett's recent resolution of the flexible side of the Onsager conjecture. Here, we in fact follow the joint work of De Lellis-Sz\'ekelyhidi Jr. and the authors of this paper, in which weak solutions of the Euler equations in the regularity class $C^{\frac 13-}_{x,t}$ are constructed, attaining any energy profile. Third, we give a concise proof of the authors' recent result, which proves the existence of infinitely many weak solutions of the Navier-Stokes in the regularity class $C^0_t L^{2+}_x \cap C^0_t W^{1,1+}_x$. We conclude the article by mentioning a number of open problems at the intersection of convex integration and hydrodynamic turbulence., Comment: 76 pages, 10 figures, minor corrections

We consider a sequence of Leray-Hopf weak solutions of the 2D Navier-Stokes equations on a bounded domain, in the vanishing viscosity limit. We provide sufficient conditions on the associated vorticity measures, away from the boundary, which ensure that as the viscosity vanishes the sequence converges to a weak solution of the Euler equations. These assumptions are consistent with vortex sheet solutions of the Euler equations., Comment: 14 pages

We prove non-uniqueness for a class of weak solutions to the Navier-Stokes equations which have bounded kinetic energy, integrable vorticity, and are smooth outside a fractal set of singular times with Hausdorff dimension strictly less than 1., Comment: 31 pages, minor corrections, to appear in JEMS

We address the local well-posedness of the hydrostatic Navier-Stokes equations. These equations, sometimes called reduced Navier-Stokes/Prandtl, appear as a formal limit of the Navier-Stokes system in thin domains, under certain constraints on the aspect ratio and the Reynolds number. It is known that without any structural assumption on the initial data, real-analyticity is both necessary and sufficient for the local well-posedness of the system. In this paper we prove that for convex initial data, local well-posedness holds under simple Gevrey regularity.

Coherent vortices are often observed to persist for long times in turbulent 2D flows even at very high Reynolds numbers and are observed in experiments and computer simulations to potentially be asymptotically stable in a weak sense for the 2D Euler equations. We consider the incompressible 2D Euler equations linearized around a radially symmetric, strictly monotone decreasing vorticity distribution. For sufficiently regular data, we prove the inviscid damping of the $\theta$-dependent radial and angular velocity fields with the optimal rates $\|u^r(t)\| \lesssim \langle t \rangle^{-1}$ and $\|u^\theta(t)\| \lesssim \langle t \rangle^{-2}$ in the appropriate radially weighted $L^2$ spaces. We moreover prove that the vorticity weakly converges back to radial symmetry as $t \rightarrow \infty$, a phenomenon known as vortex axisymmetrization in the physics literature, and characterize the dynamics in higher Sobolev spaces. Furthermore, we prove that the $\theta$-dependent angular Fourier modes in the vorticity are ejected from the origin as $t \to \infty$, resulting in faster inviscid damping rates than those possible with passive scalar evolution. This non-local effect is called vorticity depletion. Our work appears to be the first to find vorticity depletion relevant for the dynamics of vortices., Comment: 124 pages, 1 figure

For initial datum of finite kinetic energy, Leray has proven in 1934 that there exists at least one global in time finite energy weak solution of the 3D Navier-Stokes equations. In this paper we prove that weak solutions of the 3D Navier-Stokes equations are not unique in the class of weak solutions with finite kinetic energy. Moreover, we prove that Holder continuous dissipative weak solutions of the 3D Euler equations may be obtained as a strong vanishing viscosity limit of a sequence of finite energy weak solutions of the 3D Navier-Stokes equations., Comment: 36 pages; minor corrections per referee comments; version accepted by Annals of Mathematics

We discuss a lubrication approximation model of the interface between two immiscible fluids in a Hele-Shaw cell, derived in \cite{CDGKSZ93} and widely studied since. The model consists of a single one dimensional evolution equation for the thickness $2h = 2h(x,t)$ of a thin neck of fluid, \[ \partial_t h + \partial_x( h \, \partial_x^3 h) = 0\, , \] for $x\in (-1,1)$ and $t\ge 0$. The boundary conditions fix the neck height and the pressure jump: \[ h(\pm 1,t) = 1, \qquad \partial_{x}^2 h(\pm 1,t) = P>0. \] We prove that starting from smooth and positive $h$, as long as $h(x,t) >0$, for $x\in [-1,1], \; t\in [0,T]$, no singularity can arise in the solution up to time $T$. As a consequence, we prove for any $P>2$ and any smooth and positive initial datum that the solution pinches off in either finite or infinite time, i.e., $\inf_{[-1,1]\times[0,T_*)} h = 0$, for some $T_* \in (0,\infty]$. These facts have been long anticipated on the basis of numerical and theoretical studies.

We prove that any weak space-time $L^2$ vanishing viscosity limit of a sequence of strong solutions of Navier-Stokes equations in a bounded domain of ${\mathbb{R}}^2$ satisfies the Euler equation if the solutions' local enstrophies are uniformly bounded. We also prove that $t-a.e.$ weak $L^2$ inviscid limits of solutions of 3D Navier-Stokes equations in bounded domains are weak solutions of the Euler equation if they locally satisfy a scaling property of their second order structure function. The conditions imposed are far away from boundaries, and wild solutions of Euler equations are not a priori excluded in the limit.

We prove that given any $\beta<1/3$, a time interval $[0,T]$, and given any smooth energy profile $e \colon [0,T] \to (0,\infty)$, there exists a weak solution $v$ of the three-dimensional Euler equations such that $v \in C^{\beta}([0,T]\times \mathbb{T}^3)$, with $e(t) = \int_{\mathbb{T}^3} |v(x,t)|^2 dx$ for all $t\in [0,T]$. Moreover, we show that a suitable $h$-principle holds in the regularity class $C^\beta_{t,x}$, for any $\beta<1/3$. The implication of this is that the dissipative solutions we construct are in a sense typical in the appropriate space of subsolutions as opposed to just isolated examples., Comment: 36 pages, 1 figure

We prove that weak solutions of the inviscid SQG equations are not unique, thereby answering Open Problem 11 in the survey arXiv:1111.2700 by De Lellis and Sz\'ekelyhidi Jr. Moreover, we also show that weak solutions of the dissipative SQG equation are not unique, even if the fractional dissipation is stronger than the square root of the Laplacian., Comment: 49 pages

We consider the 2D Navier-Stokes equation on $\mathbb T \times \mathbb R$, with initial datum that is $\varepsilon$-close in $H^N$ to a shear flow $(U(y),0)$, where $\| U(y) - y\|_{H^{N+4}} \ll 1$ and $N>1$. We prove that if $\varepsilon \ll \nu^{1/2}$, where $\nu$ denotes the inverse Reynolds number, then the solution of the Navier-Stokes equation remains $\varepsilon$-close in $H^1$ to $(e^{t \nu \partial_{yy}}U(y),0)$ for all $t>0$. Moreover, the solution converges to a decaying shear flow for times $t \gg \nu^{-1/3}$ by a mixing-enhanced dissipation effect, and experiences a transient growth of gradients. In particular, this shows that the stability threshold in finite regularity scales no worse than $\nu^{1/2}$ for 2D shear flows close to the Couette flow.

In 1980, van Dommelen and Shen provided a numerical simulation that predicts the spontaneous generation of a singularity in the Prandtl boundary layer equations from a smooth initial datum, for a nontrivial Euler background. In this paper we provide a proof of this numerical conjecture by rigorously establishing the finite time blowup of the boundary layer thickness.

We consider the vanishing viscosity limit of the Navier-Stokes equations in a half space, with Dirichlet boundary conditions. We prove that the inviscid limit holds in the energy norm if the product of the components of the Navier-Stokes solutions are equicontinuous at $x_2=0$. A sufficient condition for this to hold is that the tangential Navier-Stokes velocity remains uniformly bounded and has a uniformly integrable tangential gradient near the boundary., Comment: We thank Gung-Min Gie and Jim Kelliher for many useful comments about the paper, which are incorporated in this new version

We consider a model of electroconvection motivated by studies of the motion of a two dimensional annular suspended smectic film under the influence of an electric potential maintained at the boundary by two cylindrical electrodes. We prove that this electroconvection model has global in time unique smooth solutions.