What’s new

Updates on my research and expository papers, discussion of open problems, and other maths-related topics. By Terence Tao

# Localisation and compactness properties of the Navier-Stokes global regularity problem

4 August, 2011 in math.AP, math.MP, paper | Tags: concentration compactness, energy estimates, global well-posedness, homogenisation, localisation, Navier-Stokes equations

I’ve just uploaded to the arXiv my paper “Localisation and compactness properties of the Navier-Stokes global regularity problem“, submitted to Analysis and PDE. This paper concerns the global regularity problem for the Navier-Stokes system of equations

$\displaystyle \partial_t u + (u \cdot \nabla) u = \Delta u - \nabla p + f \ \ \ \ \ (1)$

$\displaystyle \nabla \cdot u = 0 \ \ \ \ \ (2)$

$\displaystyle u(0,\cdot) = u_0 \ \ \ \ \ (3)$

in three dimensions. Thus, we specify initial data ${(u_0,f,T)}$, where ${0 < T < \infty}$ is a time, ${u_0: {\bf R}^3 \rightarrow {\bf R}^3}$ is the initial velocity field (which, in order to be compatible with (2), (3), is required to be divergence-free), ${f: [0,T] \times {\bf R}^3 \rightarrow {\bf R}^3}$ is the forcing term, and then seek to extend this initial data to a solution ${(u,p,u_0,f,T)}$ with this data, where the velocity field ${u: [0,T] \times {\bf R}^3 \rightarrow {\bf R}^3}$ and pressure term ${p: [0,T] \times {\bf R}^3 \rightarrow {\bf R}}$ are the unknown fields.

Roughly speaking, the global regularity problem asserts that given every smooth set of initial data ${(u_0,f,T)}$, there exists a smooth solution ${(u,p,u_0,f,T)}$ to the Navier-Stokes equation with this data. However, this is not a good formulation of the problem because it does not exclude the possibility that one or more of the fields ${u_0, f, u, p}$ grows too fast at spatial infinity. This problem is evident even for the much simpler heat equation

$\displaystyle \partial_t u = \Delta u$

$\displaystyle u(0,\cdot) = u_0.$

As long as one has some mild conditions at infinity on the smooth initial data ${u_0: {\bf R}^3 \rightarrow {\bf R}}$ (e.g. polynomial growth at spatial infinity), then one can solve this equation using the fundamental solution of the heat equation:

$\displaystyle u(t,x) = \frac{1}{(4\pi t)^{3/2}} \int_{{\bf R}^3} u_0(y) e^{-|x-y|^2/4t}\ dy.$

If furthermore ${u}$ is a tempered distribution, one can use Fourier-analytic methods to show that this is the unique solution to the heat equation with this data. But once one allows sufficiently rapid growth at spatial infinity, existence and uniqueness can break down. Consider for instance the backwards heat kernel

$\displaystyle u(t,x) = \frac{1}{(4\pi(T-t))^{3/2}} e^{|x|^2/(T-t)}$

for some ${T>0}$, which is smooth (albeit exponentially growing) at time zero, and is a smooth solution to the heat equation for ${0 \leq t < T}$, but develops a dramatic singularity at time ${t=T}$. A famous example of Tychonoff from 1935, based on a power series construction, also shows that uniqueness for the heat equation can also fail once growth conditions are removed. An explicit example of non-uniqueness for the heat equation is given by the contour integral

$\displaystyle u(t,x_1,x_2,x_3) = \int_\gamma \exp(e^{\pi i/4} x_1 z + e^{5\pi i/8} z^{3/2} - itz^2)\ dz$

where ${\gamma}$ is the ${L}$-shaped contour consisting of the positive real axis and the upper imaginary axis, with ${z^{3/2}}$ being interpreted with the standard branch (with cut on the negative axis). One can show by contour integration that this function solves the heat equation and is smooth (but rapidly growing at infinity), and vanishes for ${t<0}$, but is not identically zero for ${t>0}$.

Thus, in order to obtain a meaningful (and physically realistic) problem, one needs to impose some decay (or at least limited growth) hypotheses on the data ${u_0,f}$ and solution ${u,p}$ in addition to smoothness. For the data, one can impose a variety of such hypotheses, including the following:

• (Finite energy data) One has ${\|u_0\|_{L^2_x({\bf R}^3)} < \infty}$ and ${\| f \|_{L^\infty_t L^2_x([0,T] \times {\bf R}^3)} < \infty}$.
• (${H^1}$ data) One has ${\|u_0\|_{H^1_x({\bf R}^3)} < \infty}$ and ${\| f \|_{L^\infty_t H^1_x([0,T] \times {\bf R}^3)} < \infty}$.
• (Schwartz data) One has ${\sup_{x \in {\bf R}^3} ||x|^m \nabla_x^k u_0(x)| < \infty}$ and ${\sup_{(t,x) \in [0,T] \times {\bf R}^3} ||x|^m \nabla_x^k \partial_t^l f(t,x)| < \infty}$ for all ${m,k,l \geq 0}$.
• (Periodic data) There is some ${0 < L < \infty}$ such that ${u_0(x+Lk) = u_0(x)}$ and ${f(t,x+Lk) = f(t,x)}$ for all ${(t,x) \in [0,T] \times {\bf R}^3}$ and ${k \in {\bf Z}^3}$.
• (Homogeneous data) ${f=0}$.

Note that smoothness alone does not necessarily imply finite energy, ${H^1}$, or the Schwartz property. For instance, the (scalar) function ${u(x) = \exp( i |x|^{10} ) (1+|x|)^{-2}}$ is smooth and finite energy, but not in ${H^1}$ or Schwartz. Periodicity is of course incompatible with finite energy, ${H^1}$, or the Schwartz property, except in the trivial case when the data is identically zero.

Similarly, one can impose conditions at spatial infinity on the solution, such as the following:

• (Finite energy solution) One has ${\| u \|_{L^\infty_t L^2_x([0,T] \times {\bf R}^3)} < \infty}$.
• (${H^1}$ solution) One has ${\| u \|_{L^\infty_t H^1_x([0,T] \times {\bf R}^3)} < \infty}$ and ${\| u \|_{L^2_t H^2_x([0,T] \times {\bf R}^3)} < \infty}$.
• (Partially periodic solution) There is some ${0 < L < \infty}$ such that ${u(t,x+Lk) = u(t,x)}$ for all ${(t,x) \in [0,T] \times {\bf R}^3}$ and ${k \in {\bf Z}^3}$.
• (Fully periodic solution) There is some ${0 < L < \infty}$ such that ${u(t,x+Lk) = u(t,x)}$ and ${p(t,x+Lk) = p(t,x)}$ for all ${(t,x) \in [0,T] \times {\bf R}^3}$ and ${k \in {\bf Z}^3}$.

(The ${L^2_t H^2_x}$ component of the ${H^1}$ solution is for technical reasons, and should not be paid too much attention for this discussion.) Note that we do not consider the notion of a Schwartz solution; as we shall see shortly, this is too restrictive a concept of solution to the Navier-Stokes equation.

Finally, one can downgrade the regularity of the solution down from smoothness. There are many ways to do so; two such examples include

• (${H^1}$ mild solutions) The solution is not smooth, but is ${H^1}$ (in the preceding sense) and solves the equation (1) in the sense that the Duhamel formula

$\displaystyle u(t) = e^{t\Delta} u_0 + \int_0^t e^{(t-t')\Delta} (-(u\cdot\nabla) u-\nabla p+f)(t')\ dt'$

holds.

• (Leray-Hopf weak solution) The solution ${u}$ is not smooth, but lies in ${L^\infty_t L^2_x \cap L^2_t H^1_x}$, solves (1) in the sense of distributions (after rewriting the system in divergence form), and obeys an energy inequality.

Finally, one can ask for two types of global regularity results on the Navier-Stokes problem: a qualitative regularity result, in which one merely provides existence of a smooth solution without any explicit bounds on that solution, and a quantitative regularity result, which provides bounds on the solution in terms of the initial data, e.g. a bound of the form

$\displaystyle \| u \|_{L^\infty_t H^1_x([0,T] \times {\bf R}^3)} \leq F( \|u_0\|_{H^1_x({\bf R}^3)} + \|f\|_{L^\infty_t H^1_x([0,T] \times {\bf R}^3)}, T )$

for some function ${F: {\bf R}^+ \times {\bf R}^+ \rightarrow {\bf R}^+}$. One can make a further distinction between local quantitative results, in which ${F}$ is allowed to depend on ${T}$, and global quantitative results, in which there is no dependence on ${T}$ (the latter is only reasonable though in the homogeneous case, or if ${f}$ has some decay in time).

By combining these various hypotheses and conclusions, we see that one can write down quite a large number of slightly different variants of the global regularity problem. In the official formulation of the regularity problem for the Clay Millennium prize, a positive correct solution to either of the following two problems would be accepted for the prize:

• Conjecture 1.4 (Qualitative regularity for homogeneous periodic data) If ${(u_0,0,T)}$ is periodic, smooth, and homogeneous, then there exists a smooth partially periodic solution ${(u,p,u_0,0,T)}$ with this data.
• Conjecture 1.3 (Qualitative regularity for homogeneous Schwartz data) If ${(u_0,0,T)}$ is Schwartz and homogeneous, then there exists a smooth finite energy solution ${(u,p,u_0,0,T)}$ with this data.

(The numbering here corresponds to the numbering in the paper.)

Furthermore, a negative correct solution to either of the following two problems would also be accepted for the prize:

• Conjecture 1.6 (Qualitative regularity for periodic data) If ${(u_0,f,T)}$ is periodic and smooth, then there exists a smooth partially periodic solution ${(u,p,u_0,f,T)}$ with this data.
• Conjecture 1.5 (Qualitative regularity for Schwartz data) If ${(u_0,f,T)}$ is Schwartz, then there exists a smooth finite energy solution ${(u,p,u_0,f,T)}$ with this data.

I am not announcing any major progress on these conjectures here. What my paper does study, though, is the question of whether the answer to these conjectures is somehow sensitive to the choice of formulation. For instance:

1. Note in the periodic formulations of the Clay prize problem that the solution is only required to be partially periodic, rather than fully periodic; thus the pressure has no periodicity hypothesis. One can ask the extent to which the above problems change if one also requires pressure periodicity.
2. In another direction, one can ask the extent to which quantitative formulations of the Navier-Stokes problem are stronger than their qualitative counterparts; in particular, whether it is possible that each choice of initial data in a certain class leads to a smooth solution, but with no uniform bound on that solution in terms of various natural norms of the data.
3. Finally, one can ask the extent to which the conjecture depends on the category of data. For instance, could it be that global regularity is true for smooth periodic data but false for Schwartz data? True for Schwartz data but false for smooth ${H^1}$ data? And so forth.

One motivation for the final question (which was posed to me by my colleague, Andrea Bertozzi) is that the Schwartz property on the initial data ${u_0}$ tends to be instantly destroyed by the Navier-Stokes flow. This can be seen by introducing the vorticity ${\omega := \nabla \times u}$. If ${u(t)}$ is Schwartz, then from Stokes’ theorem we necessarily have vanishing of certain moments of the vorticity, for instance:

$\displaystyle \int_{{\bf R}^3} \omega_1 (x_2^2-x_3^2)\ dx = 0.$

On the other hand, some integration by parts using (1) reveals that such moments are usually not preserved by the flow; for instance, one has the law

$\displaystyle \partial_t \int_{{\bf R}^3} \omega_1(t,x) (x_2^2-x_3^2)\ dx = 4\int_{{\bf R}^3} u_2(t,x) u_3(t,x)\ dx,$

and one can easily concoct examples for which the right-hand side is non-zero at time zero. This suggests that the Schwartz class may be unnecessarily restrictive for Conjecture 1.3 or Conjecture 1.5.

My paper arose out of an attempt to address these three questions, and ended up obtaining partial results in all three directions. Roughly speaking, the results that address these three questions are as follows:

1. (Homogenisation) If one only assumes partial periodicity instead of full periodicity, then the forcing term ${f}$ becomes irrelevant. In particular, Conjecture 1.4 and Conjecture 1.6 are equivalent.
2. (Concentration compactness) In the ${H^1}$ category (both periodic and nonperiodic, homogeneous or nonhomogeneous), the qualitative and quantitative formulations of the Navier-Stokes global regularity problem are essentially equivalent.
3. (Localisation) The (inhomogeneous) Navier-Stokes problems in the Schwartz, smooth ${H^1}$, and finite energy categories are essentially equivalent to each other, and are also implied by the (fully) periodic version of these problems.

The first two of these families of results are relatively routine, drawing on existing methods in the literature; the localisation results though are somewhat more novel, and introduce some new local energy and local enstrophy estimates which may be of independent interest.

Broadly speaking, the moral to draw from these results is that the precise formulation of the Navier-Stokes equation global regularity problem is only of secondary importance; modulo a number of caveats and technicalities, the various formulations are close to being equivalent, and a breakthrough on any one of the formulations is likely to lead (either directly or indirectly) to a comparable breakthrough on any of the others.

This is only a caricature of the actual implications, though. Below is the diagram from the paper indicating the various formulations of the Navier-Stokes equations, and the known implications between them:

The above three streams of results are discussed in more detail below the fold.

— 1. Homogenisation —

We first discuss the homogenisation results. Let us first give the fully periodic version of the periodic regularity conjectures:

• Conjecture 1.13 (Qualitative fully periodic regularity for periodic data) If ${(u_0,f,T)}$ is periodic and smooth, then there exists a smooth fully periodic solution ${(u,p,u_0,f,T)}$ with this data.
• Conjecture 7.5 (Qualitative fully periodic regularity for homogeneous periodic data) If ${(u_0,0,T)}$ is periodic, homogeneous, and smooth, then there exists a smooth fully periodic solution ${(u,p,u_0,0,T)}$ with this data.

The implications are then that Conjecture 1.4, Conjecture 1.6, and Conjecture 7.5 are equivalent, with Conjecture 1.13 implying any of the previous three conjectures.

The equivalence of Conjecture 1.4 and Conjecture 7.5 is fairly well known (it is explicitly made for instance in my previous paper). The only remaining nontrivial implication is the deduction of Conjecture 1.6 from Conjecture 1.4, i.e. the deduction of global regularity for the inhomogeneous periodic problem from global regularity for the homogeneous periodic problem.

The deduction relies on the trick of using an asymptotic version of the Galilean symmetry ${(u,p,u_0,f,T) \mapsto (\tilde u,\tilde p,\tilde u_0, \tilde f, T)}$, where

$\displaystyle \tilde u(t,x) := u(t,x-\int_0^t v(s)\ ds) + v(t)$

$\displaystyle \tilde p(t,x) := p(t,x-\int_0^t v(s)\ ds) - x \cdot v'(t)$

$\displaystyle \tilde u_0(x) := u_0(x) + v(0)$

$\displaystyle \tilde f(t,x) := f(t,x-\int_0^t v(s)\ ds),$

and ${v: [0,T]\rightarrow {\bf R}^3}$ is an arbitrary smooth velocity function. One can verify that this symmetry preserves the Navier-Stokes system of equations. One can apply this symmetry to some rapidly growing velocity, say ${v = 2wt}$ for some large vector ${w \in {\bf R}^3}$. Then the transformation does not change the initial data ${u_0}$, but introduces a rapid time oscillation to the forcing term ${f}$, replacing it by ${\tilde f(t,x) := f(t,x-wt^2)}$. As one sends ${w}$ to infinity in a suitable “irrational” direction, ${\tilde f}$ converges weakly to a function that is constant in space (this is an application of the Riemann-Lebesgue lemma). It turns out that the Navier-Stokes system of equations has enough smoothing (or compactness) properties that we may then (for ${w}$ large and irrational enough) approximate ${f}$ effectively as if it was a constant. Applying a further Galilean-type symmetry

$\displaystyle \tilde u(t,x) := u(t,x-\int_0^t v(s)\ ds) + v(t)$

$\displaystyle \tilde p(t,x) := p(t,x-\int_0^t v(s)\ ds)$

$\displaystyle \tilde u_0(x) := u_0(x) + v(0)$

$\displaystyle \tilde f(t,x) := f(t,x-\int_0^t v(s)\ ds) + v'(t),$

for a suitable ${v}$ (not the same as the previous ${v}$), we can then set ${\tilde f}$ to be identically zero, without adjusting the initial data ${u_0}$. Thus, we have effectively transformed the inhomogeneous problem to a homogeneous one, which gives the deduction of Conjecture 1.6 from Conjecture 1.4.

Note that at one stage we subtracted a linear term ${x \cdot v'(t)}$ from the pressure. This would destroy any periodicity properties the pressure had, and so this trick does not work in the fully periodic setting.

One interpretation of this result is that the partially periodic setting is not the right setting to discuss the inhomogeneous problem. (For the homogeneous problem, it is not difficult to use Galilean transformations to show that the global regularity problem for the partially periodic and fully periodic problem are equivalent.) Indeed, the Galilean invariance reveals a breakdown of uniqueness for the Navier-Stokes flow in the partially periodic setting. Once one fixes the pressure to be periodic also, it is possible (by a variety of means, e.g. energy methods) to recover uniqueness for the flow (modulo the rather trivial caveat that the Navier-Stokes equations only determine the pressure ${p}$ up to a time-dependent constant, which does not have any net impact on (1) since ${p}$ only appears through its gradient).

— 2. Concentration compactness —

The compactness results in the paper relate qualitative and quantitative versions of the Navier-Stokes regularity problem. This portion of the argument is based on my previous paper in this topic, which dealt with the periodic homogeneous case; in the non-periodic setting it also uses the concentration-compactness method as developed by Bahouri-Gerard, by Gerard, and by Gallagher.

In the periodic homogeneous setting, the results in my previous paper show the equivalence of Conjecture 1.4 (and hence Conjectures 1.6 and 7.5) with the following conjectures:

• Conjecture 7.4 (Global existence of ${H^1}$ mild solutions for homogeneous periodic data) If ${(u_0,0,T)}$ is periodic, homogeneous, and ${H^1}$, then there exists a periodic ${H^1}$ mild solution ${(u,p,u_0,0,T)}$ with this data.
• Conjecture 7.2 (Local quantitative regularity for homogeneous periodic ${H^1}$ data) If ${(u,p,u_0,0,T)}$ is a periodic homogeneous smooth solution, then one can bound ${\|u\|_{L^\infty_t H^1_x}}$ by a function of ${\|u_0\|_{H^1_x}}$ and ${T}$.
• Conjecture 7.3 (Global quantitative regularity for homogeneous periodic ${H^1}$ data) If ${(u,p,u_0,0,T)}$ is a periodic homogeneous smooth solution, then one can bound ${\|u\|_{L^\infty_t H^1_x}}$ by a function of ${\|u_0\|_{H^1_x}}$.

The equivalences were largely based on a compactness property of the periodic Navier-Stokes flow (which was inherited by a corresponding compactness property of the periodic heat equation), in that for short times, the flow map mapped weakly convergent data in ${H^1}$ to strongly convergent solutions, and in particular maps closed balls in ${H^1}$ (which are weakly compact) to compact subsets; because continuous functions on compact sets are bounded, we can then convert qualitative estimates to quantitative ones. The global form of quantitative regularity was obtained by using the fact (from energy estimates) that the energy dissipation (which is essentially the ${H^1}$ norm of the solution ${u}$) must eventually become small.

It turns out that the same method partially applies in the inhomogeneous setting (after stating the conjectures properly, see Conjectures 1.14 and 1.15 in the text); we will not detail this here. Note that the inhomogeneous equation can add energy to the system, as well as dissipate it out, and so we do not have the equivalence between local and global quantitative estimates in this setting.

The method also works in the homogeneous non-periodic setting (with the right formulations of the conjectures, which are Conjectures 1.9, 1.17, 1.18, 1.19 in the text), but with the key additional difficulty that the Navier-Stokes flow is no longer compact, but merely concentration-compact, due to the non-compact translation symmetry that is available for the space ${H^1_x({\bf R}^3)}$. (Concentration compactness is discussed in these previous blog posts.) One then has to deal with sequences of data that are not strongly convergent, but are essentially the superposition of a number of profiles that are being translated off to infinity in different directions. However, the theory for dealing with this is well-developed; in particular, the paper of Gallagher already treats this sort of profile decomposition for Navier-Stokes in the critical regularity regime, which is more difficult than the subcritical regularity regime considered here. As such, we were able to use standard techniques to obtain the equivalences here. The main point here is that if one superimposes two profiles that are translated to be sufficiently far from each other, then the nonlinear interactions between the two profiles are weak enough that one can use perturbative techniques to obtain an approximate principle of superposition; the evolution of the superimposed profiles is essentially just the superposition of the individual profiles.

— 3. Localisation —

The most novel arguments in the paper concern the the localisation results, which connect the periodic, Schwartz, smooth finite energy, and smooth ${H^1}$ categories to each other. Before stating the results precisely, let us consider a simpler (but still unsolved) situation, namely that of the supercritical nonlinear wave equation

$\displaystyle - \partial_{tt} u + c^2 \Delta u = u^7, \ \ \ \ \ (4)$

for some scalar field ${u: [0,T] \times {\bf R}^3 \rightarrow {\bf R}}$ with a given smooth initial position ${u(0,x) = u_0(x)}$ and (for simplicity) zero initial velocity ${\partial_t u(0,x) = 0}$, and where ${c>0}$ is a fixed constant (the speed of light). One can phrase the global regularity problem in this setting for periodic smooth ${u_0}$, Schwartz ${u_0}$, smooth finite energy ${u_0}$ (where the energy is now the ${H^1_x({\bf R}^3)}$ norm rather than the ${L^2_x({\bf R}^3)}$ norm), and so forth. But for this equation, all formulations of the problem are logically equivalent, thanks to the finite speed of propagation property. Among other things, this property asserts that if one has two solutions ${u, \tilde u}$ to the equation (4) that initially agree on a ball ${B(x_0,R)}$, then for subsequent times ${t>0}$, the solutions will continue to agree on a slightly smaller ball ${B(x_0,R-ct)}$ (of course, this statement becomes vacuous once ${t}$ is large enough). This property allows us, for instance, to deduce global regularity for arbitrary smooth data (of unlimited growth) from, say, the Schwartz data case. Indeed, to obtain a solution up to time ${T}$ for smooth data ${u_0}$, one could cover the domain ${{\bf R}^3}$ by balls ${B(x_i,T)}$ of radius ${T}$, then for each such ball ${B(x_i,T)}$, smoothly truncate the data ${u_0}$ to the ball ${B(x_i,3T)}$ in such a way that it still agreed with ${u_0}$ on ${B(x_i,2T)}$. This truncated data is Schwartz, and so by hypothesis can be extended to a smooth solution up to time ${[0,T] \times B(x_i,T)}$; from finite speed of propagation we see that these partially defined solutions agree with each other on their common domain of definition, and can thus be glued together to form a global solution for the original data. The same sort of argument (combined with the trick of embedding a ball such as ${B(x_i,T)}$ into a sufficiently large torus) lets one deduce the non-periodic global regularity conjecture for (4) from the periodic one.

The finite speed of propagation property for (4) is proven by energy estimates; one basically computes the energy at time ${t}$ on the ball ${B(x_0,R-ct)}$ and shows via integration by parts that this local energy is non-increasing in time.

The Navier-Stokes equation (1) unfortunately does not enjoy finite speed of propagation. However, due to the transport term ${(u \cdot \nabla)u}$ in (1), it is reasonable to expect that the solution propagates at velocity ${u}$ (this is of course consistent with the physical interpretation of ${u}$ as the velocity field, though with the caveat that ${u}$ is actually the particle velocity rather than the group velocity). As such, if one had an a priori bound on the ${L^1_t L^\infty_x}$ norm of ${u}$, this would suggest that solutions to the Navier-Stokes equation only propagate themselves by a bounded distance. This would, heuristically at least, allow one to repeat the above types of arguments to equate the various forms of the Navier-Stokes conjecture.

It turns out, somewhat remarkably, that an a priori bound on the norm ${L^1_t L^\infty_x}$ of ${u}$ is indeed available. This does not resolve the main difficulty in the global regularity problem – the lack of a controlled coercive quantity that is either subcritical or critical, as discussed in this post – because this norm is supercritical. Actually, one can guess at the existence of such a bound by using the amplitude-frequency heuristics, discussed in this post. If the Navier-Stokes solution ${u}$ has an amplitude ${A}$ and a frequency ${N}$ for a period of time ${\tau}$ over a volume ${V}$ of space, then the energy dissipation estimate

$\displaystyle \int_0^T \int_{{\bf R}^3} |\nabla u(t,x)|^2\ dx \leq \int_{{\bf R}^3} \frac{1}{2} |u_0(t,x)|^2\ dx$

(setting ${f=0}$ for simplicity) suggests that

$\displaystyle \tau N^2 A^2 V \lesssim 1,$

while the uncertainty principle (discussed in this post) suggests that ${V \gtrsim N^{-3}}$. Finally, the linear term ${\Delta u}$ and nonlinear term ${(u \cdot \nabla) u}$ have heuristic magnitudes about ${N^2 A}$ and ${A^2 N}$ respectively, so nonlinear-dominant behaviour should only occur when ${A \gtrsim N}$. Putting all this together, one soon calculates heuristically that

$\displaystyle \tau A \lesssim 1,$

which then predicts the ${L^1_t L^\infty_x}$ bound. It turns out that one can make this argument rigorous by a routine application of Littlewood-Paley theory.

Even with the rigorous bound on ${\|u\|_{L^1_t L^\infty_x}}$, though, there is still work to do to make the localisation results rigorous. The main technical tool used in this paper is a localisation result for the enstrophy

$\displaystyle \frac{1}{2} \int_{{\bf R}^3} |\omega(t,x)|^2\ dx.$

The precise local enstrophy estimate we derive is technical, but roughly speaking, this result asserts that if the enstrophy is initially small on some ball ${B(x_0,R)}$, then it will remain small on the ball ${B(x_0,R-O(1))}$. It is proven somewhat similarly to finite speed of propagation for wave equations, in that one computes the local enstrophy adapted to a ball ${B(x_0,R(t))}$ of shrinking radius and tries to keep this enstrophy from blowing up. Whereas ${R(t)}$ was previously shrinking at a constant speed ${c}$, though, now one needs to shrink the radius at a speed proportional to ${\|u(t)\|_{L^\infty_x}}$.

A key difficulty though arises from non-local effects. If one takes the curl of (1) one obtains the vorticity equation

$\displaystyle \partial_t \omega + (u \cdot \nabla) \omega = (\omega \cdot \nabla) u + \nabla \times f.$

The problem is with the nonlinear term ${(\omega \cdot \nabla) u}$. Using the divergence-free nature of ${u}$, one can solve for ${u}$ in terms of ${\omega}$ via the Biot-Savart law

$\displaystyle u = \Delta^{-1} \nabla \times \omega.$

Unfortunately this law is non-local; the value of ${u}$ at a given point ${x}$ can be influenced by the vorticity at quite distant parts of space. In particular, the velocity field ${u}$ inside the ball ${B(x_0,R(t))}$ is influenced by the vorticity outside of ${B(x_0,R(t))}$, which will not be controlled by local enstrophy. One can still control this influence using other expressions, such as the total energy, but these are supercritical quantities and if one relies on them too heavily, then it will be impossible to keep the local enstrophy under control.

To resolve this, one has to perform some delicate harmonic analysis (in particular, a Whitney decomposition of the ball ${B(x_0,R(t))}$, and local elliptic regularity), and to carefully choose the right cutoffs to define local enstrophy properly; it turns out that one can (barely) reduce the non-local effects mentioned earlier to the point where they can actually be controlled satisfactorily by the quantities one has to play with (namely, local enstrophy and the total energy).

Using the local enstrophy inequality, one can deduce global regularity for Schwartz or smooth ${H^1}$ data from the periodic global well-posedness result in ${H^1}$. The idea is to locate a large ball ${B(0,R)}$ outside of which the data has small enstrophy; using the local enstrophy inequality, combined with standard partial regularity methods, this keeps the solution uniformly smooth (spatially, at least) outside of a slightly larger ball, say ${B(0,2R)}$. One can then truncate this solution to be compactly supported in, say, ${B(0,3R)}$, at the cost of introducing a forcing term to compensate for the error terms introduced by the truncation. One can then embed the resulting object into a sufficiently large torus to deduce the smooth ${H^1}$ regularity problem from the periodic theory. A related argument (using weak compactness, as in the construction of Leray-Hopf weak solutions) almost allows one to construct global smooth solutions from finite energy data once one has smooth solutions from finite enstrophy data.

In both of these results, though, there is an annoying technical quirk that prevents the results from being stated as cleanly as one might initially hope, which is that the partial regularity methods alluded to earlier give plenty of regularity in space in regions of small enstrophy, but very little regularity in time. This seems to be an inherent feature of the Navier-Stokes problem: if a singularity or near-singularity occurs at one point ${x_0}$ in time, then this instantaneously propagates (via the incompressibility of the fluid) to create discontinuities in time (but not in space) for the pressure at points ${x}$ very far from ${x_0}$. (Think of how, for instance, a heart beat can cause instantaneous change in blood pressure throughout the body, or how a collision in Newton’s cradle instantaneously creates an effect at the other end of the cradle.) This effect prevents one from completely localising the effect of a singularity to a bounded region of space, and causes some unpleasant technicalities to be introduced into the various implications one can prove by this method. For instance, the solutions ${u}$ we construct from smooth finite energy data are not completely smooth; they are smooth for all positive times (thanks to parabolic regularising effects) and are ${C^1_t C^\infty_x}$ in a neighbourhood of the initial time ${t=0}$, but I was unable to get more regularity in time beyond ${C^1}$.

### Recent Comments

Sandeep Murthy on An elementary non-commutative…
Luqing Ye on 245A, Notes 2: The Lebesgue…
Frank on Soft analysis, hard analysis,…
andrescaicedo on Soft analysis, hard analysis,…
Richard Palais on Pythagoras’ theorem
The Coffee Stains in… on Does one have to be a genius t…
Benoît Régent-Kloeck… on (Ingrid Daubechies) Planning f…
Luqing Ye on 245B, Notes 7: Well-ordered se…
Luqing Ye on 245B, Notes 7: Well-ordered se…
Arjun Jain on 245B, Notes 7: Well-ordered se…
%anchor_text% on Books
Luqing Ye on 245B, Notes 7: Well-ordered se…
Arjun Jain on 245B, Notes 7: Well-ordered se…
Luqing Ye on 245A, Notes 2: The Lebesgue…

## 9 comments

missed some word in the first sentence.:) [Corrected, thanks - T.]

5 August, 2011 at 2:54 am

mfrasca

Hi Terry,

In the first sentence of point 3.localisation there is a double “the”.

Best,

Marco

[Corrected, thanks - T.]

We have seen attempts at creating a “pseudo” pressure grad term that is actually the sum of the pressure grad and a “hidden” external force. Perhaps it would have been better to have put limits on the external supply of energy rather than the force which has turned out to be awkward to bind. For example, any point cannot exceed a predefined external energy extraction rate per unit volume (power draw density). This would be coupled with a total external energy supply limit.

6 August, 2011 at 11:51 pm

vydonhu

Thank Terry!

[...] and compactness properties of the Navier-Stokes global regularity problem“, discussed in this previous blog post.  As it turns out, I was somewhat impatient to finalise the paper and move on to other things, and [...]

Hi Terry,

I have very much enjoyed your posts on the N-S equations. I have been thinking about these ‘annoying technical quirks’ of the NS equations. As a fluid physicist, it is at times deeply unsettling that I tacitly allow for effects to propagate instantaneously in the formulations as written and analyzed. I wanted to see if you (or others) worry about the non-physicality of these mathematical representations of fluids, and whether non-regularity could simply be an inevitable result of this physical contradiction.

This problem seems to arise when using the kinematic approximation of the Biot-Savart (valid only for a steady flow) in the dynamic NS equations for an unsteady process. I see more complexity being needed in the definition of vorticity in these flows. e.g. replacing Biot with something like: http://en.wikipedia.org/wiki/Jefimenko%27s_equations

This would be akin to how E-M fields were proven to change each other in time and space (Maxwell’s original 1861 proof was in terms of ‘molecular vorticies’!) I believe accurately describing the mechanism of local flow organization (in space and time) via the vorticity could be key to unlocking this process.

Well, just about any useful mathematical model makes nonphysical assumptions – for instance, fluids are almost always modeled by a continuum, when in reality they are composed of a huge number of interacting particles. But if the model is robust enough, one can still expect it to give an accurate prediction of reality, even if at an ontological level it is quite distinct.

With regard to the instantaneous propagation of the pressure, one could use the compressible Navier-Stokes equation instead of the incompressible one, which would presumably be a more faithful model, and one without the instantaneous pressure issue. But in the incompressible limit I believe that it is expected that solutions of the former equation should converge to that of the latter (assuming, of course, that solutions to the latter are regular). There is quite a bit of literature on these topics; see for instance this paper of Desjardins, Grenier, Lions, and Masmoudi.

My feeling though is that this non-locality of pressure, while technically annoying, is not the main factor in determining the regularity, stability, or validity of the Navier-Stokes flow. The nonlocal aspects of the Biot-Savart law primarily affect the low frequency (i.e. coarse scale) components of the flow, whereas all the interesting nonlinear and turbulent behaviour is occuring at fine scales. Of course, in this highly nonlinear equation the coarse scales and fine scales are not completely decoupled, and so the non-locality of the Biot-Savart law does have some impact on the fine-scale dynamics, but it is indirect and I would be surprised if it ends up playing a decisive role in the regularity theory (except at the initial time t=0, where as noted in the post, slowly decaying fine-scale behaviour in the velocity can cause instantaneous loss of time smoothness for the pressure).

Think EthanEagle may have a point, though I feel a bit like an acolyte commenting to the masters.

While I am certainly nowhere near of your level in mathematics, could we not posit that the mathematical model need not be treated as a per se continuum or at the level of abstraction that it is. Based on the 2-d NS solution it would seem that the pressure terms are the problem.

Just as a thought model suppose we remove the concept of pressure and treat the fluid mass as a bounded variable step size n dimensional matrix populated with point sources (containing all the variables with which a molecule would be endowed), with n degrees of freedom and thermally defined mean free path step sizes.

Introduce a planer heat source and planer cold sink on either side of the fluid mass and one would simultaneously have an expansion wave propagating through the medium from the planar heat source, an increase in the mean free path of the area surrounding the heat source, and a decrease in the mean free path of the area surrounding the cold sink. This would introduce a density gradient from one side of the mass to the other. The random motion of the molecules within the mass would engender a momentum flux in the sparse direction. If this were occurring with gravity the expansion wave would propagate across the upper portion of the mass and the momentum flux from dense to sparse would propagate along the lower portion, i.e. a natural convective loop. The steeper the gradient the stronger the flow except in the region where the expansion and momentum force interact and counterbalance each other in the middle of the counter rotating flow which would also be the region of minimum bulk velocity and maximum turbulence provided the walls of the “container” did not have physical properties.

With this thought model it is a natural extension that turbulent regions are a function of the presence of introduced local density imbalances combined with energy conditions under which the directionally chaotic energy of the molecules expresses itself in structured turbulent motion simply because the energy of the flow near the source of turbulence is insufficient to apparently suppress the random motion (except at mach the random motion would never be suppressed) but instead “interacts” with it imparting a degree of local structure to the essentially chaotic motion. This thought model would satisfy the no-slip condition and the formation of turbulent boundary layers near physical objects.

Simplified the model would parse the flow into uni-directional and omni-directional momentum flux components based on local and global density and energy conditions. From the continuum perspective this would not be dissimilar to the pressure-velocity relationship, i.e. when pressure is at a maximum omni-directional flux is at a maximum and when velocity is at a maximum unidirectional flux is at a maximum.

I have been working on ambient flows (mainly wind) for the last 6 years and the pressure conceptual model is particularly problematic in accurately explaining or predicting the dynamics of these types of flows under certain conditions in the absence of a large body of empirical data.

It may be worth noting that the conceptual framework of incompressible vs. compressible flows comes from late 19th and early 20th century aerodynamics and forced flow environments such as wind tunnels, planes, or pipes, i.e. either driving a fluid through or at a physical object or driving a physical object through a fluid. It was convenient when analyzing those types of flows (which are most of the flows we need to analyze) to deal with subsonic flows as incompressible. It makes sense as any “compressible” effects are negligible compared to the energy contained in the forced flow except in the boundary layer where, notably, turbulence develops.

Perhaps the solution to NS may be to reformulate the state equations without the pressure term. That would of course mean stepping way back and having to reformulate all the associated equations sans pressure.

Such an approach would also inherently solve the coarse vs fine issue based on the resolution of the applied matrix. Only real problem would be introducing the right basis functions to generate the chaotic turbulent behavior in regions of high turbulence. It would never be fully accurate at describing the “reality” of the turbulence as the turbulence in this model is inherently chaotic but it should do a decent job of coming quite close at the fine scale and be very accurate at the coarse scale.

12 September, 2011 at 5:13 am

Tim Nguyen

Hi Terry, it seems that your blog posts are not updating on your buzz (at least for me).

Cancel
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 227, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "github"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "english_confidence": 0.9227301478385925, "perplexity_flag": "head"}