It is always dangerous to venture an opinion as to why a problem is hard (cf. Clarke’s first law), but I’m going to stick my neck out on this one, because (a) it seems that there has been a lot of effort expended on this problem recently, sometimes perhaps without full awareness of the main difficulties, and (b) I would love to be proved wrong on this opinion :-) .
The global regularity problem for Navier-Stokes is of course a Clay Millennium Prize problem and it would be redundant to describe it again here. I will note, however, that it asks for existence of global smooth solutions to a Cauchy problem for a nonlinear PDE. There are countless other global regularity results of this type for many (but certainly not all) other nonlinear PDE; for instance, global regularity is known for Navier-Stokes in two spatial dimensions rather than three (this result essentially dates all the way back to Leray’s thesis in 1933!). Why is the three-dimensional Navier-Stokes global regularity problem considered so hard, when global regularity for so many other equations is easy, or at least achievable?
(For this post, I am only considering the global regularity problem for Navier-Stokes, from a purely mathematical viewpoint, and in the precise formulation given by the Clay Institute; I will not discuss at all the question as to what implications a rigorous solution (either positive or negative) to this problem would have for physics, computational fluid dynamics, or other disciplines, as these are beyond my area of expertise. But if anyone qualified in these fields wants to make a comment along these lines, by all means do so.)
The standard response to this question is turbulence – the behaviour of three-dimensional Navier-Stokes equations at fine scales is much more nonlinear (and hence unstable) than at coarse scales. I would phrase the obstruction slightly differently, as supercriticality. Or more precisely, all of the globally controlled quantities for Navier-Stokes evolution which we are aware of (and we are not aware of very many) are either supercritical with respect to scaling, which means that they are much weaker at controlling fine-scale behaviour than controlling coarse-scale behaviour, or they are non-coercive, which means that they do not really control the solution at all, either at coarse scales or at fine. (I’ll define these terms more precisely later.) At present, all known methods for obtaining global smooth solutions to a (deterministic) nonlinear PDE Cauchy problem require either
- Exact and explicit solutions (or at least an exact, explicit transformation to a significantly simpler PDE or ODE);
- Perturbative hypotheses (e.g. small data, data close to a special solution, or more generally a hypothesis which involves an
somewhere); or
- One or more globally controlled quantities (such as the total energy) which are both coercive and either critical or subcritical.
(Note that the presence of (1), (2), or (3) are currently necessary conditions for a global regularity result, but far from sufficient; otherwise, papers on the global regularity problem for various nonlinear PDEs would be substantially shorter :-) . In particular, there have been many good, deep, and highly non-trivial papers recently on global regularity for Navier-Stokes, but they all assume either (1), (2) or (3) via additional hypotheses on the data or solution. For instance, in recent years we have seen good results on global regularity assuming (2), as well as good results on global regularity assuming (3); a complete bibilography of recent results is unfortunately too lengthy to be given here.)
The Navier-Stokes global regularity problem for arbitrary large smooth data lacks all of these three ingredients. Reinstating (2) is impossible without changing the statement of the problem, or adding some additional hypotheses; also, in perturbative situations the Navier-Stokes equation evolves almost linearly, while in the non-perturbative setting it behaves very nonlinearly, so there is basically no chance of a reduction of the non-perturbative case to the perturbative one unless one comes up with a highly nonlinear transform to achieve this (e.g. a naive scaling argument cannot possibly work). Thus, one is left with only three possible strategies if one wants to solve the full problem:
- Solve the Navier-Stokes equation exactly and explicitly (or at least transform this equation exactly and explicitly to a simpler equation);
- Discover a new globally controlled quantity which is both coercive and either critical or subcritical; or
- Discover a new method which yields global smooth solutions even in the absence of the ingredients (1), (2), and (3) above.
For the rest of this post I refer to these strategies as “Strategy 1″, “Strategy 2″, and “Strategy 3″.
Much effort has been expended here, especially on Strategy 3, but the supercriticality of the equation presents a truly significant obstacle which already defeats all known methods. Strategy 1 is probably hopeless; the last century of experience has shown that (with the very notable exception of completely integrable systems, of which the Navier-Stokes equations is not an example) most nonlinear PDE, even those arising from physics, do not enjoy explicit formulae for solutions from arbitrary data (although it may well be the case that there are interesting exact solutions from special (e.g. symmetric) data). Strategy 2 may have a little more hope; after all, the Poincaré conjecture became solvable (though still very far from trivial) after Perelman introduced a new globally controlled quantity for Ricci flow (the Perelman entropy) which turned out to be both coercive and critical. (See also my exposition of this topic.) But we are still not very good at discovering new globally controlled quantities; to quote Klainerman, “the discovery of any new bound, stronger than that provided by the energy, for general solutions of any of our basic physical equations would have the significance of a major event” (emphasis mine).
I will return to Strategy 2 later, but let us now discuss Strategy 3. The first basic observation is that the Navier-Stokes equation, like many other of our basic model equations, obeys a scale invariance: specifically, given any scaling parameter , and any smooth velocity field
solving the Navier-Stokes equations for some time T, one can form a new velocity field
to the Navier-Stokes equation up to time
, by the formula
(Strictly speaking, this scaling invariance is only present as stated in the absence of an external force, and with the non-periodic domain rather than the periodic domain
. One can adapt the arguments here to these other settings with some minor effort, the key point being that an approximate scale invariance can play the role of a perfect scale invariance in the considerations below. The pressure field
gets rescaled too, to
, but we will not need to study the pressure here. The viscosity
remains unchanged.)
We shall think of the rescaling parameter as being large (e.g.
). One should then think of the transformation from u to
as a kind of “magnifying glass”, taking fine-scale behaviour of u and matching it with an identical (but rescaled, and slowed down) coarse-scale behaviour of
. The point of this magnifying glass is that it allows us to treat both fine-scale and coarse-scale behaviour on an equal footing, by identifying both types of behaviour with something that goes on at a fixed scale (e.g. the unit scale). Observe that the scaling suggests that fine-scale behaviour should play out on much smaller time scales than coarse-scale behaviour (T versus
). Thus, for instance, if a unit-scale solution does something funny at time 1, then the rescaled fine-scale solution will exhibit something similarly funny at spatial scales
and at time
. Blowup can occur when the solution shifts its energy into increasingly finer and finer scales, thus evolving more and more rapidly and eventually reaching a singularity in which the scale in both space and time on which the bulk of the evolution is occuring has shrunk to zero. In order to prevent blowup, therefore, we must arrest this motion of energy from coarse scales (or low frequencies) to fine scales (or high frequencies). (There are many ways in which to make these statements rigorous, for instance using Littlewood-Paley theory, which we will not discuss here, preferring instead to leave terms such as “coarse-scale” and “fine-scale” undefined.)
Now, let us take an arbitrary large-data smooth solution to Navier-Stokes, and let it evolve over a very long period of time [0,T), assuming that it stays smooth except possibly at time T. At very late times of the evolution, such as those near to the final time T, there is no reason to expect the solution to resemble the initial data any more (except in perturbative regimes, but these are not available in the arbitrary large-data case). Indeed, the only control we are likely to have on the late-time stages of the solution are those provided by globally controlled quantities of the evolution. Barring a breakthrough in Strategy 2, we only have two really useful globally controlled (i.e. bounded even for very large T) quantities:
- The maximum kinetic energy
; and
- The cumulative energy dissipation
.
Indeed, the energy conservation law implies that these quantities are both bounded by the initial kinetic energy E, which could be large (we are assuming our data could be large) but is at least finite by hypothesis.
The above two quantities are coercive, in the sense that control of these quantities imply that the solution, even at very late times, stays in a bounded region of some function space. However, this is basically the only thing we know about the solution at late times (other than that it is smooth until time T, but this is a qualitative assumption and gives no bounds). So, unless there is a breakthrough in Strategy 2, we cannot rule out the worst-case scenario that the solution near time T is essentially an arbitrary smooth divergence-free vector field which is bounded both in kinetic energy and in cumulative energy dissipation by E. In particular, near time T the solution could be concentrating the bulk of its energy into fine-scale behaviour, say at some spatial scale . (Of course, cumulative energy dissipation is not a function of a single time, but is an integral over all time; let me suppress this fact for the sake of the current discussion.)
Now, let us take our magnifying glass and blow up this fine-scale behaviour by to create a coarse-scale solution to Navier-Stokes. Given that the fine-scale solution could (in the worst-case scenario) be as bad as an arbitrary smooth vector field with kinetic energy and cumulative energy dissipation at most E, the rescaled unit-scale solution can be as bad as an arbitrary smooth vector field with kinetic energy and cumulative energy dissipation at most
, as a simple change-of-variables shows. Note that the control given by our two key quantities has worsened by a factor of
; because of this worsening, we say that these quantities are supercritical – they become increasingly useless for controlling the solution as one moves to finer and finer scales. This should be contrasted with critical quantities (such as the energy for two-dimensional Navier-Stokes), which are invariant under scaling and thus control all scales equally well (or equally poorly), and subcritical quantities, control of which becomes increasingly powerful at fine scales (and increasingly useless at very coarse scales).
Now, suppose we know of examples of unit-scale solutions whose kinetic energy and cumulative energy dissipation are as large as , but which can shift their energy to the next finer scale, e.g. a half-unit scale, in a bounded amount O(1) of time. Given the previous discussion, we cannot rule out the possibility that our rescaled solution behaves like this example. Undoing the scaling, this means that we cannot rule out the possibility that the original solution will shift its energy from spatial scale
to spatial scale
in time
. If this bad scenario repeates over and over again, then convergence of geometric series shows that the solution may in fact blow up in finite time. Note that the bad scenarios do not have to happen immediately after each other (the self-similar blowup scenario); the solution could shift from scale
to
, wait for a little bit (in rescaled time) to “mix up” the system and return to an “arbitrary” (and thus potentially “worst-case”) state, and then shift to
, and so forth. While the cumulative energy dissipation bound can provide a little bit of a bound on how long the system can “wait” in such a “holding pattern”, it is far too weak to prevent blowup in finite time. To put it another way, we have no rigorous, deterministic way of preventing Maxwell’s demon from plaguing the solution at increasingly frequent (in absolute time) intervals, invoking various rescalings of the above scenario to nudge the energy of the solution into increasingly finer scales, until blowup is attained.
Thus, in order for Strategy 3 to be successful, we basically need to rule out the scenario in which unit-scale solutions with arbitrarily large kinetic energy and cumulative energy dissipation shift their energy to the next highest scale. But every single analytic technique we are aware of (except for those involving exact solutions, i.e. Strategy 1) requires at least one bound on the size of solution in order to have any chance at all. Basically, one needs at least one bound in order to control all nonlinear errors – and any strategy we know of which does not proceed via exact solutions will have at least one nonlinear error that needs to be controlled. The only thing we have here is a bound on the scale of the solution, which is not a bound in the sense that a norm of the solution is bounded; and so we are stuck.
To summarise, any argument which claims to yield global regularity for Navier-Stokes via Strategy 3 must inevitably (via the scale invariance) provide a radically new method for providing non-trivial control of nonlinear unit-scale solutions of arbitrary large size for unit time, which looks impossible without new breakthroughs on Strategy 1 or Strategy 2. (There are a couple of loopholes that one might try to exploit: one can instead try to refine the control on the “waiting time” or “amount of mixing” between each shift to the next finer scale, or try to exploit the fact that each such shift requires a certain amount of energy dissipation, but one can use similar scaling arguments to the preceding to show that these types of loopholes cannot be exploited without a new bound along the lines of Strategy 2, or some sort of argument which works for arbitrarily large data at unit scales.)
To rephrase in an even more jargon-heavy manner: the “energy surface” on which the dynamics is known to live in, can be quotiented by the scale invariance. After this quotienting, the solution can stray arbitrarily far from the origin even at unit scales, and so we lose all control of the solution unless we have exact control (Strategy 1) or can significantly shrink the energy surface (Strategy 2).
The above was a general critique of Strategy 3. Now I’ll turn to some known specific attempts to implement Strategy 3, and discuss where the difficulty lies with these:
- Using weaker or approximate notions of solution (e.g. viscosity solutions, penalised solutions, super- or sub- solutions, etc.). This type of approach dates all the way back to Leray. It has long been known that by weakening the nonlinear portion of Navier-Stokes (e.g. taming the nonlinearity), or strengthening the linear portion (e.g. introducing hyperdissipation), or by performing a discretisation or regularisation of spatial scales, or by relaxing the notion of a “solution”, one can get global solutions to approximate Navier-Stokes equations. The hope is then to take limits and recover a smooth solution, as opposed to a mere global weak solution, which was already constructed by Leray for Navier-Stokes all the way back in 1933. But in order to ensure the limit is smooth, we need convergence in a strong topology. In fact, the same type of scaling arguments used before basically require that we obtain convergence in either a critical or subcritical topology. Absent a breakthrough in Strategy 2, the only type of convergences we have are in very rough – in particular, in supercritical – topologies. Attempting to upgrade such convergence to critical or subcritical topologies is the qualitative analogue of the quantitative problems discussed earlier, and ultimately faces the same problem (albeit in very different language) of trying to control unit-scale solutions of arbitrarily large size. Working in a purely qualitative setting (using limits, etc.) instead of a quantitative one (using estimates, etc.) can disguise these problems (and, unfortunately, can lead to errors if limits are manipulated carelessly), but the qualitative formalism does not magically make these problems disappear. Note that weak solutions are already known to be badly behaved for the closely related Euler equation. More generally, by recasting the problem in a sufficiently abstract formalism (e.g. formal limits of near-solutions), there are a number of ways to create an abstract object which could be considered as a kind of generalised solution, but the moment one tries to establish actual control on the regularity of this generalised solution one will encounter all the supercriticality difficulties mentioned earlier.
- Iterative methods (e.g. contraction mapping principle, Nash-Moser iteration, power series, etc.) in a function space. These methods are perturbative, and require something to be small: either the data has to be small, the nonlinearity has to be small, or the time of existence desired has to be small. These methods are excellent for constructing local solutions for large data, or global solutions for small data, but cannot handle global solutions for large data (running into the same problems as any other Strategy 3 approach). These approaches are also typically rather insensitive to the specific structure of the equation, which is already a major warning sign since one can easily construct (rather artificial) systems similar to Navier-Stokes for which blowup is known to occur. The optimal perturbative result is probably very close to that established by Koch-Tataru, for reasons discussed in that paper.
- Exploiting blowup criteria. Perturbative theory can yield some highly non-trivial blowup criteria – that certain norms of the solution must diverge if the solution is to blow up. For instance, a celebrated result of Beale-Kato-Majda shows that the maximal vorticity must have a divergent time integral at the blowup point. However, all such blowup criteria are subcritical or critical in nature, and thus, barring a breakthrough in Strategy 2, the known globally controlled quantities cannot be used to reach a contradiction. Scaling arguments similar to those given above show that perturbative methods cannot achieve a supercritical blowup criterion.
- Asymptotic analysis of the blowup point(s). Another proposal is to rescale the solution near a blowup point and take some sort of limit, and then continue the analysis until a contradiction ensues. This type of approach is useful in many other contexts (for instance, in understanding Ricci flow). However, in order to actually extract a useful limit (in particular, one which still solves Navier-Stokes in a strong sense, and does collapse to the trivial solution), one needs to uniformly control all rescalings of the solution – or in other words, one needs a breakthrough in Strategy 2. Another major difficulty with this approach is that blowup can occur not just at one point, but can conceivably blow up on a one-dimensional set; this is another manifestation of supercriticality.
- Analysis of a minimal blowup solution. This is a strategy, initiated by Bourgain, which has recently been very successful in establishing large data global regularity for a variety of equations with a critical conserved quantity, namely to assume for contradiction that a blowup solution exists, and then extract a minimal blowup solution which minimises the conserved quantity. This strategy (which basically pushes the perturbative theory to its natural limit) seems set to become the standard method for dealing with large data critical equations. It has the appealing feature that there is enough compactness (or almost periodicity) in the minimal blowup solution (once one quotients out by the scaling symmetry) that one can begin to use subcritical and supercritical conservation laws and monotonicity formulae as well (see my survey on this topic). Unfortunately, as the strategy is currently understood, it does not seem to be directly applicable to a supercritical situation (unless one simply assumes that some critical norm is globally bounded) because it is impossible, in view of the scale invariance, to minimise a non-scale-invariant quantity.
- Abstract approaches (avoiding the use of properties specific to the Navier-Stokes equation). At its best, abstraction can efficiently organise and capture the key difficulties of a problem, placing the problem in a framework which allows for a direct and natural resolution of these difficulties without being distracted by irrelevant concrete details. (Kato’s semigroup method is a good example of this in nonlinear PDE; regrettably for this discussion, it is limited to subcritical situations.) At its worst, abstraction conceals the difficulty within some subtle notation or concept (e.g. in various types of convergence to a limit), thus incurring the risk that the difficulty is “magically” avoided by an inconspicuous error in the abstract manipulations. An abstract approach which manages to breezily ignore the supercritical nature of the problem thus looks very suspicious. More substantively, there are many equations which enjoy a coercive conservation law yet still can exhibit finite time blowup (e.g. the mass-critical focusing NLS equation); an abstract approach thus would have to exploit some subtle feature of Navier-Stokes which is not present in all the examples in which blowup is known to be possible. Such a feature is unlikely to be discovered abstractly before it is first discovered concretely; the field of PDE has proven to be the type of mathematics where progress generally starts in the concrete and then flows to the abstract, rather than vice versa.
If we abandon Strategy 1 and Strategy 3, we are thus left with Strategy 2 – discovering new bounds, stronger than those provided by the (supercritical) energy. This is not a priori impossible, but there is a huge gap between simply wishing for a new bound and actually discovering and then rigorously establishing one. Simply sticking in the existing energy bounds into the Navier-Stokes equation and seeing what comes out will provide a few more bounds, but they will all be supercritical, as a scaling argument quickly reveals. The only other way we know of to create global non-perturbative deterministic bounds is to discover a new conserved or monotone quantity. In the past, when such quantities have been discovered, they have always been connected either to geometry (symplectic, Riemmanian, complex, etc.), to physics, or to some consistently favourable (defocusing) sign in the nonlinearity (or in various “curvatures” in the system). There appears to be very little usable geometry in the equation; on the one hand, the Euclidean structure enters the equation via the diffusive term and by the divergence-free nature of the vector field, but the nonlinearity is instead describing transport by the velocity vector field, which is basically just an arbitrary volume-preserving infinitesimal diffeomorphism (and in particular does not respect the Euclidean structure at all). One can try to quotient out by this diffeomorphism (i.e. work in material coordinates) but there are very few geometric invariants left to play with when one does so. (In the case of the Euler equations, the vorticity vector field is preserved modulo this diffeomorphism, as observed for instance by Li, but this invariant is very far from coercive, being almost purely topological in nature.) The Navier-Stokes equation, being a system rather than a scalar equation, also appears to have almost no favourable sign properties, in particular ruling out the type of bounds which the maximum principle or similar comparison principles can give. This leaves physics, but apart from the energy, it is not clear if there are any physical quantities of fluids which are deterministically monotone. (Things look better on the stochastic level, in which the laws of thermodynamics might play a role, but the Navier-Stokes problem, as defined by the Clay institute, is deterministic, and so we have Maxwell’s demon to contend with.) It would of course be fantastic to obtain a fourth source of non-perturbative controlled quantities, not arising from geometry, physics, or favourable signs, but this looks somewhat of a long shot at present. Indeed given the turbulent, unstable, and chaotic nature of Navier-Stokes, it is quite conceivable that in fact no reasonable globally controlled quantities exist beyond that which arise from the energy.
Of course, given how hard it is to show global regularity, one might try instead to establish finite time blowup instead (this also is acceptable for the Millennium prize). Unfortunately, even though the Navier-Stokes equation is known to be very unstable, it is not clear at all how to pass from this to a rigorous demonstration of a blowup solution. All the rigorous finite time blowup results (as opposed to mere instability results) that I am aware of rely on one or more of the following ingredients:
- Exact blowup solutions (or at least an exact transformation to a significantly simpler PDE or ODE, for which blowup can be established);
- An ansatz for a blowup solution (or approximate solution), combined with some nonlinear stability theory for that ansatz;
- A comparison principle argument, dominating the solution by another object which blows up in finite time, taking the solution with it; or
- An indirect argument, constructing a functional of the solution which must attain an impossible value in finite time (e.g. a quantity which is manifestly non-negative for smooth solutions, but must become negative in finite time).
It may well be that there is some exotic symmetry reduction which gives (1), but no-one has located any good exactly solvable special case of Navier-Stokes (in fact, those which have been found, are known to have global smooth solutions). (2) is problematic for two reasons: firstly, we do not have a good ansatz for a blowup solution, but perhaps more importantly it seems hopeless to establish a stability theory for any such ansatz thus created, as this problem is essentially a more difficult version of the global regularity problem, and in particular subject to the main difficulty, namely controlling the highly nonlinear behaviour at fine scales. (One of the ironies in pursuing method (2) is that in order to establish rigorous blowup in some sense, one must first establish rigorous stability in some other (renormalised) sense.) Method (3) would require a comparison principle, which as noted before appears to be absent for the non-scalar Navier-Stokes equations. Method (4) suffers from the same problem, ultimately coming back to the “Strategy 2″ problem that we have virtually no globally monotone quantities in this system to play with (other than energy monotonicity, which clearly looks insufficient by itself). Obtaining a new type of mechanism to force blowup other than (1)-(4) above would be quite revolutionary, not just for Navier-Stokes; but I am unaware of even any proposals in these directions, though perhaps topological methods might have some effectiveness.
So, after all this negativity, do I have any positive suggestions for how to solve this problem? My opinion is that Strategy 1 is impossible, and Strategy 2 would require either some exceptionally good intuition from physics, or else an incredible stroke of luck. Which leaves Strategy 3 (and indeed, I think one of the main reasons why the Navier-Stokes problem is interesting is that it forces us to create a Strategy 3 technique). Given how difficult this strategy seems to be, as discussed above, I only have some extremely tentative and speculative thoughts in this direction, all of which I would classify as “blue-sky” long shots:
- Work with ensembles of data, rather than a single initial datum. All of our current theory for deterministic evolution equations deals only with a single solution from a single initial datum. It may be more effective to work with parameterised familes of data and solutions, or perhaps probability measures (e.g. Gibbs measures or other invariant measures). One obvious partial result to shoot for is to try to establish global regularity for generic large data rather than all large data; in other words, acknowledge that Maxwell’s demon might exist, but show that the probability of it actually intervening is very small. The problem is that we have virtually no tools for dealing with generic (average-case) data other than by treating all (worst-case) data; the enemy is that the Navier-Stokes flow itself might have some perverse entropy-reducing property which somehow makes the average case drift towards (or at least recur near) the worst case over long periods of time. This is incredibly unlikely to be the truth, but we have no tools to prevent it from happening at present.
- Work with a much simpler (but still supercritical) toy model. The Navier-Stokes model is parabolic, which is nice, but is complicated in many other ways, being relatively high-dimensional and also non-scalar in nature. It may make sense to work with other, simplified models which still contain the key difficulty that the only globally controlled quantities are supercritical. Examples include the Katz-Pavlovic dyadic model for the Euler equations (for which blowup can be demonstrated by a monotonicity argument; see this survey for more details), or the spherically symmetric defocusing supercritical nonlinear wave equation.
- Develop non-perturbative tools to control deterministic non-integrable dynamical systems. Throughout this post we have been discussing PDEs, but actually there are similar issues arising in the nominally simpler context of finite-dimensional dynamical systems (ODEs). Except in perturbative contexts (such as the neighbourhood of a fixed point or invariant torus), the long-time evolution of a dynamical system for deterministic data is still largely only controllable by the classical tools of exact solutions, conservation laws and monotonicity formulae; a discovery of a new and effective tool for this purpose would be a major breakthrough. One natural place to start is to better understand the long-time, non-perturbative dynamics of the classical three-body problem, for which there are still fundamental unsolved questions.
- Establish really good bounds for critical or nearly-critical problems. Recently, I showed that having a very good bound for a critical equation essentially implies that one also has a global regularity result for a slightly supercritical equation. The idea is to use a monotonicity formula which does weaken very slightly as one passes to finer and finer scales, but such that each such passage to a finer scale costs a significant amount of monotonicity; since there is only a bounded amount of monotonicity to go around, it turns out that the latter effect just barely manages to overcome the former in my equation to recover global regularity (though by doing so, the bounds worsen from polynomial in the critical case to double exponential in my logarithmically supercritical case). I severely doubt that my method can push to non-logarithmically supercritical equations, but it does illustrate that having very strong bounds at the critical level may lead to some modest progress on the problem.
- Try a topological method. This is a special case of (1). It may well be that a primarily topological argument may be used either to construct solutions, or to establish blowup; there are some precedents for this type of construction in elliptic theory. Such methods are very global by nature, and thus not restricted to perturbative or nearly-linear regimes. However, there is no obvious topology here (except possibly for that generated by the vortex filaments) and as far as I know, there is not even a “proof-of-concept” version of this idea for any evolution equation. So this is really more of a wish than any sort of concrete strategy.
- Understand pseudorandomness. This is an incredibly vague statement; but part of the difficulty with this problem, which also exists in one form or another in many other famous problems (e.g. Riemann hypothesis,
,
, twin prime and Goldbach conjectures, normality of digits of
, Collatz conjecture, etc.) is that we expect any sufficiently complex (but deterministic) dynamical system to behave “chaotically” or “pseudorandomly”, but we still have very few tools for actually making this intuition precise, especially if one is considering deterministic initial data rather than generic data. Understanding pseudorandomness in other contexts, even dramatically different ones, may indirectly shed some insight on the turbulent behaviour of Navier-Stokes.
In conclusion, while it is good to occasionally have a crack at impossible problems, just to try one’s luck, I would personally spend much more of my time on other, more tractable PDE problems than the Clay prize problem, though one should certainly keep that problem in mind if, in the course on working on other problems, one indeed does stumble upon something that smells like a breakthrough in Strategy 1, 2, or 3 above. (In particular, there are many other serious and interesting questions in fluid equations that are not anywhere near as difficult as global regularity for Navier-Stokes, but still highly worthwhile to resolve.)

127 comments
Comments feed for this article
18 March, 2007 at 7:14 pm
Greg Kuperberg
This is an interesting-looking review of a topic that I don’t know beans about. But I notice that you do not mention models of fluid flow with quantized vortices. These models approximate continuous 3+1-dimensional fluid equations by non-local 1+1-dimensional equations. It seems conceivable, although I know nothing about the analytic issues, that you could hope for some kind of convergence as the quantization parameter goes to zero.
My wife, Rena Zieve, studies fluids with quantized vortices. In her case, the motivation is not that the model approximates classical fluids, but that it is a true description of superfluid helium.
18 March, 2007 at 9:07 pm
Terence Tao
Dear Greg,
It is indeed the analytic issues (and specifically, the establishing of estimates which control the solution, or convergence of approximate solutions to the true solution) which are the heart of the matter. Suppose for instance that we do manage to construct a global solution
to a quantised Navier-Stokes for every
(this is certainly conceivable, as one can do similar things for other relaxations, regularisations, or discretisations of Navier-Stokes). Now, as you say, one now needs to establish convergence in some topology of these quantised solutions in the limit
. In order to get any useful sort of limit, as a bare minimum one needs the sequence
to enjoy some sort of bound in some function space norm, uniformly in
, otherwise all sorts of bad things could happen (e.g. breakdown of the energy conservation identity). The key phrase here is uniformity in
. If one has some argument that bounds solutions to the quantised Navier-Stokes uniformly in
, then by taking limits one expects the same argument to work directly in the limit
. In other words, one can dispense with the quantised Navier-Stokes and work with the original Navier-Stokes equation directly.
Thus, for the point of view of establishing estimates (which is the key problem), perturbing the equation by varying a small parameter does not really help. Such tricks are however useful for more qualitative or formal aspects of the theory. For instance, to justify things like the energy conservation identity for very rough solutions (for which one cannot directly justify things like differentiation under the integral sign), a typical trick is to first relax the solution to a smoother solution (and perhaps a smoother equation), establish the relevant identity in that smooth setting, and then take limits to recover the identity in the original solution. These sorts of perturbations are also useful for a number of simple topological arguments, such as an application of the continuity method (showing that all solutions obey some property P by showing first that the set of solutions obeying P is open, closed, and non-empty, and then using the connectedness of the solution space).
From a physical viewpoint, it may well be that one of these modified equations is in fact a more realistic model for fluids than Navier-Stokes. But for the narrow purposes of solving the Clay Prize Problem, we’re stuck with the original Navier-Stokes equation :-) .
18 March, 2007 at 9:13 pm
Greg Kuperberg
You make it sound like it could be worthwhile to try to disprove the conjecture.
19 March, 2007 at 5:01 am
Nets Katz
Terry,
I think you are a little pessimistic about strategy 1 for proving blow-up – working with ensembles. For that strategy to work it is not absolutely essential that the typical solution blow up. In order for Navier Stokes to blow up, you just to have each scale “activated”, that is have energy flow sufficiently into each scale.
Suppose you could prove with positive probabiity that energy cascades into a generic situation in which it is one scale higher. Generic has to mean that we get positive probability of flowing to the next scale. You might get a measure zero set of solutions with blow-up but nevertheless use a probabilistic argument
to prove that it occurs.
19 March, 2007 at 8:50 am
Terence Tao
Greg, I certainly think one should pursue the blowup direction as well as the regularity problem. As noted in my main post, though, our technologies for establishing blowup are rather limited at present. Which brings me to Nets’ interesting idea. This idea may at least be able to establish a “norm inflation” scenario, in which some high regularity Sobolev norm is shown to increase quite rapidly in a bounded amount of time. That would be enough to disprove a fairly strong version of global regularity, namely that one can bound some critical global spacetime norm uniformly by, say, something depending only on the
norm of the inital data. (This would imply as a corollary by standard persistence-of-regularity arguments that any Sobolev norm of the solution at a late time is controlled by the Sobolev norm at time zero, and the
norm of the initial data.) I suspect that this type of global regularity result, while very common for critical problems, is probably false for supercritical problems, and might be disprovable by some sort of contradiction argument (e.g. by Bourgain’s induction on energy method). If one was extremely optimistic one might then hope to run a Baire category type argument to create an actual blowup, but this looks somewhat unlikely to me.
19 March, 2007 at 7:19 pm
Terence Tao
Dear Nets,
I thought about it a bit more, and this strategy to establish blowup may run up against the same supercriticality issues which plague the global regularity problem.
Thanks to the recent work of Escauriaza-Seregin-Sverak and others, we know that blowup solutions to Navier-Stokes must in fact blow up in the critical norm
. This rules out self-similar type solutions in which the solution shifts all of its energy from one scale to the next while keeping critical norms under control. It is plausible that these results can be localised, leading one to also rule out “self-similar + radiation” solutions (such as those appearing in the recent works of Merle and Raphael for NLS, as well as even more recent work on wave maps by Krieger-Schlag-Tataru and Rodnianski-Sterbenz) in which a lot of
mass is radiated away to coarse scales but the concentrating portion of the solution stays bounded in
. If these scenarios are ruled out, then the blowup solution must become increasingly large (and hence increasingly nonlinear) in critical norms at fine scales as one approaches the blowup time. This creates a lot of scope for an anti-Maxwell’s demon (Maxwell’s angel?) to cause trouble – one may start with a promisingly randomly distributed ensemble at coarse scales, but as one pushes the ensemble into fine scales with large critical norm, and the evolution becomes very nonlinear and unpredictable, Maxwell’s angel could conceivably sneak in and drain a lot of entropy out of the ensemble, and eventually nudge the entire ensemble into states in which the energy all dissipates harmlessly and one has global regularity. (Now, if we could obtain some sort of rigorous version of the second law of thermodynamics here, one might be able to prevent this from happening, but this looks very remote at present.)
19 March, 2007 at 11:20 pm
Nets Katz
Terry,
That’s really unfortunate. For this to work, nonlinearity needs to be a
friend and not an enemy. In fact the test problem should first to prove blow up
for 3D Euler with finite energy.
Some advantages there: Energy never dissipates away, cascading always
occurs exactly by vortex stretching so that some structure seems to be preserved (for instance the locations of zeroes just convect)
Of course the passage to high scales should happen at time scales so startlingly
fast that the linear term is having very little effect at all so that the blow up for
Euler should indicate the blow up for Navier Stokes.
20 March, 2007 at 12:21 am
Nets Katz
In fact, so that I can understand your objection better, let’s restrict
our attention to 3D Euler. The critical norm which you say we should be
desperately clinging to is L^{infty} of vorticity. We don’t need a recent result
here. Beale-Kato-Majda already prove that if we have blow up of Euler
then the L^{infty} of vorticity must blow up. Does that mean the equation
becomes “increasingly unstable”?
By contrast let me mention a 2D problem for Euler which is critical. There
global solvability is guaranteed since vorticity is convected so that L^{\infty}
of vorticity is conserved. Beale-Kato-Majda show that the growth of
Sobolev norms in time is doubly exponential. The question, perhaps somewhat less glamorous than global solvability, is whether this estimate
is sharp and double exponential growth actually occurs. As we’ve discussed in the past, if one looks at the highest order part of the growth
exclusively, in a certain position, one is forced to consider a system
of ODE for SL(2) valued variables, one variable for each active scale. The
criticality enforces that different scales could be of equal strength and
that many different scales have to work together to achieve double exponential growth. More over if we look at this system of N equations
which governs the growth of the frequencies that are like 2^{-N}, this
system is highly nonlinear, in that, in the relevant time scale which
is 1/N, you could have all of the N/2 highest scales change completely.
You can make very precise that this critical setting becomes increasingly
unstable as the N gets larger. In particular the system has increasing
numbers of variables more remote from each other in scale, all of which
interact. That is tough going.
On the other hand, in the supercritical case, if you posit any reasonable
worst cascading of energy, what happens is that the highest order scale
dominates exponentially the lower order scales. The cascading process
should indeed be nonlinear, but it seems self-similar as frequencies increase. In the time we can activate the N+1st scale, the Nth scale might get all messed up – that’s the nonlinearity – but we don’t care so
much what happens at slightly lower scales. The problem is essentially local in scale.
Now you might argue that this local process could have a derandomizing effect on the ensemble. (And the more times you do it, the more derandomizing it gets.) But if you could show that were true, it would
probably be the most amazing theorem in fluid mechanics ever. In any
case, in this setting, it seems to me the supercritical problem might be
easier to understand than the critical one. Am I wrong?
20 March, 2007 at 8:26 pm
Terence Tao
Dear Nets,
It’s a bit trickier to use scaling analysis for Euler, since it has a two-parameter scaling symmetry rather than a one-parameter one; the spatial scale, temporal scale, and velocity magnitude are related only by a single constraint rather than two. Nevertheless I agree that the problem of improving the double-exponential Beale-Kato-Majda bound for 2D Euler is a great problem – it is perhaps the one place where the wall separating the impossibly supercritical problems from the feasible critical problems is thinnest. It is clear that some improvement should be possible, and whatever technique is used to accomplish this will undoubtedly be interesting.
I also agree that it is against all physical intuition that fluid equations could derandomise ensembles of solutions – but we don’t seem to have any way at present to make this intuition rigorous. If we have to admit the existence of derandomisers to move the solution to either the best-case scenario or worst-case scenario at whim, then any supercritical equation could conceivably do just about anything from generic blowup to global regularity, depending on the mood of the derandomiser. This is why I feel that “understanding pseudorandomness” needs to be part of the solution (unless we find a new monotonicity formula or something, which lets us gain more control on either the best case or worst case and allows us to establish results even if the equation tries its best to derandomise the flow).
20 March, 2007 at 8:37 pm
Ars Mathematica » Blog Archive » Why Navier-Stokes is Hard
[...] Tao has a thoughtful post that explains why proving existence results for Navier-Stokes equations is so [...]
23 March, 2007 at 8:18 am
Not Even Wrong » Blog Archive » All Sorts of Links
[...] each have fascinating blog entries on Millenium problems. Terry Tao writes a long explanation of Why Global Regularity for Navier-Stokes is Hard. He also comments about the recent New York Times piece about him and about math education issues. [...]
29 March, 2007 at 5:52 am
Gil Kalai
Dear Terry
This is a wonderful post. I have a (rather generic) question: is looking at the problem at high dimensions (even asymptotically somehow) can make things easier regarding the negative direction or even regarding some aspects of the positive direction?
29 March, 2007 at 7:19 am
Terence Tao
Dear Gil,
Intuitively, the equation should behave worse in higher dimensions, because the relationship between the energy and the scaling becomes increasingly unfavourable. Certainly the equation is more unstable; but it is hard to convert instability to actual blowup (it’s the “Maxwell’s angel” problem discussed earlier – the instability might have the improbable effect of always rescuing the solution just before it blows up). Conversely, in two dimensions, where the energy becomes critical, the result of Beale-Kato-Majda shows that one has global smooth solutions for Navier Stokes.
Another popular model to study is the hyperdissipation model, in which the dispersive factor
in the Navier-Stokes equation is replaced by a power
, where
is a parameter which serves a similar purpose to dimension, in that it determines the relationship between energy and scaling. The level
is the threshold beyond which the energy becomes subcritical and one has global regularity. (I’m not sure what happens exactly at 5/4, but perhaps Nets does.) In this paper it is shown that as
moves from 5/4 continuously down to 1, the (upper bound on the) dimension of the singular set at blowup time moves continuously from 0 to 1.
31 March, 2007 at 2:10 pm
Hongjie
Dear Terry,
3D Navier-Stokes equation is globally well-posed when \alpha=5/4. More generally, n D NSE is globally well-posed when \alpha=(n+2)/4. One proof of this can be found, for example, in the paper http://arxiv.org/abs/math/0104199 you mentioned. But it is not explicitly stated there…
4 April, 2007 at 10:21 am
Y. Charles Li
Hi, Terry
This is a very nice article, will be helpful especially to fresh Ph.D.’s.
About the scaling argument, simple as it is, it does dictate important
estimates, e.g. for 3D NS, Leray obtained: For any initial condition u(0),
when
t > T = C \nu^-5 || u(0) ||^4_L^2
no more singularity. Using the scaling,
|| u(0) ||^4_L^2 —> \lambda^2 || u(0) ||^4_L^2 , T —> \lambda^2 T
as dictated by the scaling t —> \lambda^2 t. Pretend this is 2D
(of course, 2D global regularity has easier argument),
|| u(0) ||^4_L^2 —> || u(0) ||^4_L^2 , T —> T
while t —> \lambda^2 t by scaling. So any local solution can be rescaled to global.
I’m also glad you start to appreciate “turbulence”. I wrote a thing on
Mathematical Intelligencer:
http://www.math.missouri.edu/~cli/Nature-T.pdf
It should have appeared in hard copy. It is a follow-up of the famous article by
Ruelle and Takens:
http://www.springerlink.com/content/h1760361517×10h2/
I hope it brings chaos “closer” to turbulence. How close? it is hard to judge.
From a different perspective, the case is like the “global warming” problem,
the debate of whether turbulence is chaos is more or less over, the question is:
What can one do about chaos (turbulence) in terms of physical “reachable” description?
I wrote something along this line which is too immature to make public. I asked
Joel Lebowitz to read it for me and he is reading it.
Best Regards
Charles
4 April, 2007 at 10:26 am
Y. Charles Li
Ruelle and Takens paper link is not working properly:
D. Ruelle, F. Takens, On the nature of turbulence, Comm. Math. Phys., 20 (167-192),
23 (243-244), 1971.
Y. Charles Li
7 April, 2007 at 10:22 am
Dave Purvance
In the time evolution of a viscous, incompressible flow’s Fourier modes convection’s quadratic nonlinearity is expressed as a convolution of an unknown flow unshifted and shifted in wavenumber. When the unshifted portion is left to be some unknown function and its shifted counterpart is assumed to be a time series of finite order, and, when a discrete span of wavenumbers are considered simultaneously, then the three-space Navier-Stokes evolution equations become a large matrix differential equation. As the order of the time series and the number of wavenumbers considered get large, this Navier-Stokes matrix differential equation becomes ever closer to the continuous three-space Navier-Stokes evolution equations. Using conjugate symmetry and similarity transforms, it is argued at “http://arxiv.org/ftp/math/papers/0610/0610086.pdf” that the exact solution to the Navier-Stokes matrix differential equation is a stable matrix exponential and that the assumed time series is its truncated Taylor expansion in time.
Care to comment?
7 April, 2007 at 10:45 am
Terence Tao
Dear Dave,
This is an instance of the type 1 of Strategy 3 listed in my post above: “Using weaker or approximate notions of solution”. A localisation of the frequency space is essentially equivalent via the uncertainty principle to a discretisation of physical space. As discussed above, there is indeed no difficulty creating a global solution for the approximated equation. When however one tries to take a limit to create a global strong solution for the original equation, it is necessary to obtain a uniform subcritical or critical bound on all approximate solutions (i.e. to achieve Strategy 2) in order to extract a smooth limit. Otherwise, the best one can achieve are the global weak solutions of Leray.
8 April, 2007 at 6:24 am
Dave Purvance
Dear Terence,
Thanks for your reply.
Unless something strange happens to a matrix eigendecomposition as the matrix size becomes infinitely large, convection remains oscillatory (purely imaginary eigenvalues) for all time. Convection is therefore absolutely bounded by the magnitude of the initial condition. Viscous shear (negative real eigenvalues) diffuses these convective oscillations to zero as time goes to infinity. The only problem I see in taking a limit here is that convective oscillations may become extremely (infinitely?) fast.
Is this the limit problem you are talking about ?
8 April, 2007 at 10:12 am
Terence Tao
Dear Dave, this is part of the problem, but more importantly is the fact that as one transitions from finite-dimensional state spaces to infinite-dimensional ones, being bounded in one norm no longer implies being bounded in another. Thus for instance it is conceivable that limit remains bounded in energy, but is no longer smooth, which is the whole point of the Navier-Stokes global regularity property. (Bounded-energy global weak solutions were constructed all the way back in 1933 by Leray, essentially by a variant of the method you describe, but it is known that weak solutions need not be smooth.) In particular, the energy may concentrate in finer and finer length scales until a singularity occurs in finite time. This singularity is not entirely visible from the finite dimensional approximations, as they only have finitely many scales in the first place; it is a phenomenon which emerges in the limit.
In order to obtain regularity control on the limit, as opposed to merely a weak solution, it is not enough to bound the energy (which is a supercritical quantity); one also needs to bound a critical or subcritical quantity, and no such global bounds for arbitrarily large data are currently known at present (this is the “Strategy 2″ I discuss above).
8 April, 2007 at 12:02 pm
Simons Lecture III: Structure and randomness in PDE « What’s new
[...] For , one has global smooth solutions for small data with either sign. For large data in the focusing case, finite time blowup is possible. For large data in the defocusing case, the existence of global smooth solutions are unknown even for spherically symmetric data, indeed this problem, being supercritical, is of comparable difficulty to the Navier-Stokes global regularity problem. [...]
9 April, 2007 at 5:45 am
One side of "our" difficult problem (i.e. turbulence, did you think about something else?) « Alex’s Blog
[...] Why global regularity for Navier-Stokes is hard « What’s new [...]
10 April, 2007 at 3:24 pm
Dave Purvance
With the best of respect, I really believe the new bound provided by the Navier-Stokes “matrix” differential equation is the [-1,1] bound on convection (purely imaginary eigenvalues). Convection “chirps up” in frequency in all powers of increasing time, but stays bounded in magnitude by [-1,1]. Viscous shear (real negative eigenvalues) smoothly diffuses these highly nonlinear oscillatory chirps to zero as time goes to infinity. There needs to be no mention of energy.
And with this I’ll quit babbling…
15 April, 2007 at 3:11 am
sgrajeev
What a thoughtful and inspiring discussion!
Some questions
Are there other examples of supercritical PDEs that are better understood, or are the difficulties you mention generic to all such cases? e.g. \waveoperator\phi+\lambda\phi^p=0 for p big enough.
The scaling arguments you start with also appear in
quantum field theory where supercritical means something like `non-renormalizable’. Is there more to this analogy? If there is, it is bad news for proving regularity.
There is some geometry in the Euler equations: the fluid follows a geodesic w.r.t. the L^2 metric on the group of volume preserving Diffeomoprhisms (Arnold). Does this suggest anything for Navier-Stokes?
15 April, 2007 at 6:52 am
Terence Tao
Dear sgrajeev,
Most supercritical evolution equations are just as poorly understood as the Navier-Stokes equations, unfortunately, with the notable exception of the supercritical elliptic and parabolic equations for which there is a favourable sign which allows for some sort of maximum principle or comparison principle to take hold; basically, in such cases, the nonlinearity is extremely strong, but it always acting in a favourable direction, using its strength to reduce the size of the solution rather than increase it. However, in oscillatory evolution equations such as wave or Schrodinger equations it appears that there is no similar phenomenon taking place; the nonlinearity can “try” to reduce the size of its solution, but it is so strong that it can “overshoot” and end up changing the sign or direction of the solution and making it much larger at the same time.
There do seem to be analogies between classical PDE and quantum field theory (which can be viewed as a kind of quantum PDE) but this is definitely an underexplored area of study. For instance, the Cauchy problem for quantum field theory has not been studied much, even in linear models (perhaps it is a bad question to ask).
It is true that the Euler equation has some nice geometrical structure, and I would indeed think this structure will be key in the future understanding of this equation, and thus indirectly to Navier-Stokes. There is however a major obstruction with using the diffeomorphism group structure to understand Navier-Stokes, namely that this equation contains the Laplacian (via the viscosity term), which relies on the Euclidean (or Riemannian) geometry of space. This type of geometry is not preserved at all by diffeomorphisms. Because of the presence of two very different types of geometry in Navier-Stokes it seems significantly less likely that we get the same type of “geometric miracles” (e.g. unexpected monotonicity formulae) that we do in, say, Ricci flow, which only involves one type of geometry. But there could still be many surprises in this equation.
15 April, 2007 at 1:02 pm
Rajeev
Thanks.
). Euler equations are the geodesic equations for this metric, according to Arnold.
. Often dissipative equations can be thought of as a Hamiltonian system with a complex valued hamilton
. (Sorry to plug my work here, but it is relatively recent, so may not be known: arXiv:quant-ph/0701141).
The Euclidean metric on the domain does go into the Euler equations as well; e.g., in determining the metric on the group of volume preserving diffeomorphisms. The length^2 of a tangent vector to the group is just the length^2 of the velocity field w.r.t. the Euclidean metric integrated on the domain of the fluid ( which is also the kinetic energy
Viscosity introduces the Laplacian in a more direct way. The viscous force is the gradient of the function
Not claiming that this gives a better control on the scaling properties of Navier-Stokes. But, it looks like Navier-Stokes has a geometrical meaning too, using complex geometry. Might be of independent interest.
13 May, 2007 at 4:35 am
Rajeev’s Journal » Blog Archive » Fuzzy Fluids
[...] Terrence Tao has made some deep observations on why the regularity of three dimensional Navier-Stokes is such a hard problem. He has gone on to many other equally fascinating topics, I remain fixated on his main point there: that Navier-Stokes is `supercritical’. The nonlinearities become stronger at small distance scales, making it impossible to know (using present techniques) whether solutions remain smooth for all time. Thus it is crucial to understand the scale dependence of non-linearities in fluid mechanics. [...]
22 May, 2007 at 8:00 am
Smooth Solution to the 3-Space Navier-Stokes Equations by David Purvance « Smooth Chaos in Fluid Dynamics
[...] post argues that viscous shear is both the coercive and critical quantity that assures smooth solutions to the 3-space-periodic Navier-Stokes [...]
22 May, 2007 at 8:19 am
DavePurvance
This post argues that viscous shear is both the “coercive” and “critical” quantity that assures smooth solutions to the 3-space Navier-Stokes equations. Isn’t this what physicists have long been expected?
23 May, 2007 at 7:44 am
math student
Hi Prof. Tao,
I am curious that if we assume the solution is periodic , any partial result ?
23 May, 2007 at 8:33 pm
Terence Tao
Dear math student,
The difficulty in global regularity for Navier-Stokes lies at the fine spatial scales – in particular, scales much smaller than the period L. Intuitively, we expect the behaviour at such scales to not “see” the periodicity, and so whatever difficulties which are present in the non-periodic case will also be expected to be present in the periodic case.
Dear Dave,
There is an error in your post when you attempt to derive (42) from (40), trying to diagonalise different matrices A_n simultaneously. The rotation matrices T_n used to diagonalise A_n depend on n, so you cannot deduce (42) from (7) by conjugating by a single matrix.
More generally, there is no quantitative advantage gained in discretising the problem by introducing some truncation parameters (such as the parameters N and M in your post). The only bounds in the discretised model which would be of use to the original continuous model would be those bounds which are uniform in the truncation parameters, since these are the only bounds which will survive the passage back to the limit. But if there was a bound in the discretised model that was independent of those parameters, one could also have proven it directly in the continuous model simply by taking the limit of the _proof_ rather than of the result. Truncated models are useful for justifying some qualitative statements (e.g. ensuring that all sums and integrals converge, etc.) but do not progress towards the heart of the matter, which is to establish global bounds on solutions to the continuous equation.
25 May, 2007 at 7:08 am
DavePurvance
Let’s try this again…
becomes known only after
are known for
, so a simultaneous diagonalization of all convective matrices is not possible. I don’t see how diagonalizing convective matrices one at a time in increasing time-order
invalidates the Navier-Stokes matrix differential equation (7) or its matrix exponential solution (45). Could you explain? Also, the coercive, critical (Gaussian low-pass) nature of viscous shear holds independent of
and
and therefore should remain in the limit as
.
Dear Terence,
Even though it only becomes clear in the paragraphs below (47) of my post, convective matrix
I hope I am understanding your welcomed criticisms.
25 May, 2007 at 8:41 am
Terence Tao
Dear Dave,
Each time you apply one of the transformations T_n to diagonalise one of the matrices A_n in (7), you undo the diagonalisation of the preceding matrices. Thus, as I said above, the formulation (42) of your truncated Navier-Stokes equation is not valid (indeed, as you have just commented, the matrices A_n are not simultaneously diagonalisable in general, so there is no way that (42) can be deduced from (7)).
28 May, 2007 at 10:10 am
Frederick
Dear Terence,
Even one puts N-S equation in a cubic box, the essential difficulty remains, because the high frequency (ultra-violet) part remains. Can we truncate the ultra-violet part? Numerical (lattice) simulation tries to smear it, so lattice simulation is still an ostrich policy.
I guess it is the meaning of some of your meaning. I think this is right and the main difficulty of nonlinear analysis. Actually I think the whole subject of nonlinear equations is try to find some ways to deal with the ultra-violet part. Because in a nonlinear equation, one could not omit even the 1 billionth Fourier mode. Such as x^2, even x is a periodic function, whose property is known, nearly all the information about x^2 is lost. This is mystic to me, why even one nonlinear term could destroy all the information and brings lots of randomness? There seems no efficient method to calculate Fourier coefficients of x^2 and to understand property of x^2.
The above is my understanding of nonlinear equations from the point of Fourier analysis.
My question are:
1. Is there some other Fourier series other than sine, cosine that is friendly to nonlinear terms? I mean, when one expands x in some new Fourier series, then x^2 could be easily calculated. If one expands x in terms of exp(i*n*omega*t), then one has to take into consideration of all the terms such as exp(i*(m-n)*omega*t), even m=1 billion+1 and n=1 billion. Could we expand x in terms of NewFourier(n, omega, t), when we calculate x^2, maybe only several m and n that are needed to calculate accurately x^2? Have we find such new Fourier series? Or there is just no such series at all? Just as quintic polynomial equations have no general solution in radicals according to Abel–Ruffini theorem.
2. Someone saw a analogy between ultra-violet part of N-S equation and ultra-violet divergences of quantum field theory. Is the analogy proper and accurate? Surely if it is accurate, renormalization group is quite a great tool to extract information.
Frederick
29 May, 2007 at 10:34 am
DavePurvance
Dear Terence,
I have added two extra steps to equation (42) to make clear that
is used as a similarity transform on
to change, if necessary, its wavenumber flow basis from
to
.
does not effect any lower order matrix
for
. For one to argue that there is something wrong with (42), I would think one would have to argue that the 3-space Naiver-Stokes equations and its solution cannot be described in terms of a single wavenumber flow basis. I haven’t seen anyone argue this yet. Have you?
You are right, though. These similarity transforms cannot be “deduced” from (7) . They don’t invalidate (42) either.
30 May, 2007 at 4:15 am
DavePurvance
Dear Terence,
I might also add that the similarity transform
is necessary in my post’s (42) to keep the flow from “diffusing” into the one forbidden flow direction
as a passive scalar can do in the advection-diffusion equation.
(Excuse my less-than-stellar grammer. I believe in my last comment “effect” should read “affect”. Also, WordPress promises me that sooner or later they will get all of the bugs out of rendering latex in blog comments.)
3 August, 2007 at 9:31 am
2006 ICM: Étienne Ghys, “Knots and dynamics” « What’s new
[...] have no discernible structure other than that of a general diffeomorphism. As I discussed in my own post on Navier-Stokes, the discovery of a new conserved quantity for fluid equations could potentially be extremely [...]
16 August, 2007 at 3:51 pm
Stephen Montgomery-Smith
My feeling is that there aren’t any undiscovered globally controlled quantities of the type described in strategy (2). I admit that my reasoning is full of holes. But the energy estimate for the Navier-Stokes equation gives rise (via a heuristic argument) to the Kolmogorov 5/3 power law for the spectrum. However in 2D we have enstrophy estimates giving rise to the Kraichnan power of 3 law. My understanding is that both experiment and numerics confirm this, at least to some extent.
Similarly, if there were some strange, new, undiscovered globally controlled quantity, perhaps one might expect this to also give rise to a law different to the 5/3 law. But then experiment would have picked this up.
I do admit that one of the huge holes in my argument is that energy and enstrophy in 3D and 2D come equipped with well defined dissipation rates. So in a way, it is just a thought. But, for example, if one could make helicity somehow “monotone”, then this would definitely qualify.
17 August, 2007 at 8:44 am
Terence Tao
Dear Stephen,
That’s a fairly plausible argument. (One could perhaps hypothesise that there are quantities which are globally controlled, but only start becoming strong enough to affect ensemble distributions such as the power law at very small scales, beyond what experiment and numerics can detect. Of course, if that was the case, then it would be a mathematical artefact rather than a physically or computationally useful quantity, but I presume it would still be technically legal to use it in order to claim the Millennium prize. :-) )
There is also the possibility that there is some sort of “adaptive” controlled quantity, which relies on the prior history of the solution as well as on the current state, and so is not cutting down the energy surface in a way which would distort the power law. This seems to go against thermodynamical intuition though; fluids don’t seem to have much “memory”. On balance I would agree with you that Strategy (2) is unlikely to work, though it can’t be ruled out completely, and it’s not as if the other strategies have a significantly higher probability of success right now anyway. :-)
17 August, 2007 at 11:26 am
Stephen Montgomery-Smith
I know that the Navier-Stokes equation is “memoryless” but whenever one actually looks at numerics, the solution really does look like it has a lot of memory. I have observed this quite dramatically with the 2D Euler equation. The vorticity scalar is pushed around by the flow, and after a certain amount of time seems to arrange itself in “puff pastry” like thin layers of positive and negative vorticity, to the point that one could look at an instance of a flow, and make a rather good guess as to how long the flow has evolved.
From watching the Science Channel, I get the impression that a similar thing happens with the solar magnetic field. Presumably this satisfies some kind of memoryless hydromagnetic equation. But according to the T.V. shows I watch, the magnetic fields get all twisted up around each other. At some point they “snap” (presumably get stretched so much that some kind of dissipation kicks in), and this causes release of huge amounts of energy, giving rise to solar flares that cause inconvenience to man made communications satellites.
And I know that you have spent some time with the “magnetization variable” formulation of the Navier-Stokes equation. While everyone who seems to have worked on it seems to have got no-where, nevertheless the magnetization variable differs from the velocity field by the gradient of a scalar field, and perhaps that scalar field codes in a kind of abstract memory.
(For those of you unfamiliar with the magnetization variable formulation, I did a short write up at http://www.qeden.com/wiki/Navier-Stokes_Existence_and_Smoothness).
And I guess that this paper:
http://www.iumj.indiana.edu/IUMJ/fulltext.php?artid=42034&year=1993&volume=42
is an example of what an adaptive controlled quantity would be.
18 August, 2007 at 11:33 am
Stephen Montgomery-Smith
I was looking back at the old comments (I only found this page a few days ago). I saw Hongjie’s comment from March 31st that the hyperviscous NS is globally well posed for alpha=5/4. Actually this has a rather short proof, because in this situation we have the boundedness of
int_0^infty || (-Lap)^{5/8} u ||_2^2 dt.
Then the proof proceeds following the same lines as, say, the proof that the Prodi-Serrin conditions are sufficient for global well-posedness for regular NS.
20 August, 2007 at 3:31 pm
“Math Doesn’t Suck”, and the Chayes-McKellar-Winn theorem « What’s new
[...] In that case, Markov chain theory lets one conclude that if the solution started out at a fixed total energy E, and the system S was isolated, then the limiting distribution of microstates would just be the uniform distribution on the energy surface ; every state on this surface is equally likely to occur at any given instant of time (this is known as the fundamental postulate of statistical mechanics, though in this simple Markov chain model we can actually derive this postulate rigorously). This distribution is known as the microcanonical ensemble of S at energy E. It is remarkable that this ensemble is largely independent of the actual values of the transition probabilities; it is only the energy E and the function H which are relevant. (This analysis is perfectly rigorous in the Markov chain model, but in more realistic models such as Hamiltonian mechanics or quantum mechanics, it is much more difficult to rigorously justify convergence to the microcanonical ensemble. The trouble is that while these models appear to have a chaotic dynamics, which should thus exhibit very pseudorandom behaviour (similar to the genuinely random behaviour of a Markov chain model), it is very difficult to demonstrate this pseudorandomness rigorously; the same difficulty, incidentally, is present in the Navier-Stokes regularity problem.) [...]
7 September, 2007 at 6:23 am
Stephen Montgomery-Smith
I am having doubts about my earlier remark that the Kolmogorov 5/3 law suggests the non-existence of undiscovered global controlling quantities. The reason I say this is because in the 2D case we have the Kraichnan power of 3 law that comes from the enstrophy. But in the 2D case, there are also better controlling quantities than the enstrophy, namely the L_p norms of the vorticity for all p greater than 2.
2 October, 2007 at 4:15 am
Notes : Wiki page on NS equations « jtstnsp
[...] Now ,the next job is to take a look at Terry Tao’s article on NS equations. [...]
9 October, 2007 at 11:23 am
A quantitative formulation of the global regularity problem for the periodic Navier-Stokes equation « What’s new
[...] Tuesday, October 9th, 2007 in math.AP, paper Tags: compactness, Navier-Stokes equations I have just uploaded to the arXiv my paper “A quantitative formulation of the global regularity problem for the periodic Navier-Stokes equation”, submitted to Dynamics of PDE. This is a short note on one formulation of the Clay Millennium prize problem, namely that there exists a global smooth solution to the Navier-Stokes equation on the torus given any smooth divergence-free data. (I should emphasise right off the bat that I am not claiming any major breakthrough on this problem, which remains extremely challenging in my opinion.) [...]
29 December, 2007 at 12:36 am
nOnoscience » Blog Archive » Turbulence
[...] If you can handle comfortably some mathematical language, read the excellent post titled “Why Global Regularity for Navier Stokes is Hard” by Terrence [...]
16 March, 2008 at 6:21 am
cgjoh@csc.kth.se
For computational evidence of blowup of incompressible Euler solutions, see the recent article (to appear in BIT Numerical Mathematics)
Blowup of incompressible Euler solutions
available on http://www.csc.kth.se/~cgjoh
20 March, 2008 at 8:47 am
DavePurvance
Dear Terence,
In your comments above I believe you have confirmed my argument that for any spatially periodic flow
of initial finite value
the incompressible Navier-Stokes equations can be posed as the nonlinear matrix differential equation
. (1)
is a function of the unknown flow
. When
in
is expanded as a time series
(2),
. (3)
have a common flow basis, i.e., when
commute, then the solution to (1) is the stable matrix exponential
. (4)
must be assumed a priori and, therefore, it is a mistake to assert (4) is a general solution to the Navier-Stokes (1).
. (5)
solves for the unknown flow coefficients
. For
they are
. (6)
. (7)
and
converge even when
do not commute. Obviously, when
commute,
. (8)
do not commute, the illusive bound to the periodic Navier-Stokes flow is
(9)
using a white-noise initial flow. These numerical results have given me the courage to ask you again for your opinion of my recent findings (even though it was a humbling experience last time). Thank you.
This equation is nonlinear because matrix
then (1) becomes
When
However, as you correctly pointed out in your comments above, commutativity of
To remedy this commutativity problem, I have recently argued in this arXiv paper that using the time expanded flow (2) everywhere in (3) gives the equation
Matching coefficients in (5) when
Without assuming commutativity, let the Taylor expansion of the matrix exponential function (4) be
This is a stable function of
What I have also argued in this same paper is that when
The paper numerically confirms this bound (9) for
Dave Purvance
20 March, 2008 at 9:41 am
A bound for periodic Navier-Stokes flows « Smooth Chaotic Fluids
[...] as Terry Tao correctly pointed out, commutativity of must be assumed a priori and, therefore, it is a mistake [...]
20 March, 2008 at 11:10 am
DavePurvance
Dear Anonymous,
It’s only a mistake to hold that my matrix exponential function is a general solution to the periodic Navier-Stokes equation. The bound that I have identified is this same matrix exponential function, but expanded without the commutativity assumption. It is NOT claimed to be a solution, just a bound.
Dave Purvance
22 March, 2008 at 2:14 pm
Terence Tao
Dear Dave,
Unfortunately there is an error in deriving (26) in your paper. You seem to be treating various monomials involving the
as if they were non-negative numbers, when in fact they are matrices of indefinite signature. As such it is not necessary the case that a difference is bounded in magnitude by a sum; for instance,
might be larger than
(e.g. consider the case when
), and similarly
can be larger than
. Also, it is not true that matrix product order has no influence on the “bound” of the product: for instance, it is entirely possible for
to be different from
. (Both are bounded above by
, but this is not the type of expression you have on the right-hand side of (26). Similarly,
is bounded above by
, but this is again not the type of expression you have on the right-hand side of (26), in which the norm is outside the sum and product, rather than inside.)
More generally, one cannot control a non-commutative flow by its commutative counterpart, because cancellations in the latter do not necessarily force cancellations in the former. It’s easiest to illustrate this phenomenon with a discrete equation
rather than a continuous equation
. It is an easy matter to find non-commuting matrices A, B whose commutator
has an eigenvalue larger than 1. If we then let U(t) vary periodically between A, B, -A, and -B in turn and let u(0) be an eigenvector of
corresponding to that eigenvalue, then
grows exponentially in time, even though the abelianisation
stays bounded. It is easy to convert this discrete-time example to a continuous-time one but I will leave this as an exercise.
[The one time in which one can safely control a non-commutative flow by a commutative one is if the latter has no cancellation in it, in which case Gronwall's inequality applies.]
18 April, 2008 at 9:16 am
DavePurvance
Dear Terence,
While it may turn out that “commutative flows”, as you term them, may be the only flows physically possible, I did not mean my mathematical bounds to be derived from a commutative flow. My bounds, which I now denote as
with corresponding matrices
, use the commutative solution’s matrix exponential function, but this function is evaluated with a noncommutative solution’s
.
are therefore not solution coefficients to the Navier-Stokes matrix differential equation or any other differential equation that I know of, unless, of course,
become commutative.
Relating to your criticism that the inequality
(1)
(2)
with
this latter inequality asserts
. (3)
and
and the contested inequality (1) must accommodate both commutative and noncommutative cases.
is comprised of 512 matrix product terms and each matrix is 1029×1029. Using numerous “white-noise” initial flows, the sum of these matrix product terms always agrees with the bounds derived from (1). If (1) is not true, then how does one explain these (near miraculous) numerical results?
is not established, isn’t its opposite
contradictory? For instance, in the commutative case when
Finally, the numerical tests. The 10th order
18 April, 2008 at 11:35 pm
Jonas
Dave,
As far as I can tell there is no miracle going on here. That the inequality holds for white noise initial data says very little. There are numerous examples of inequalities in mathematics which hold for a typical random sample of the objects in question, however it is also often the case that the examples where the inequality is false are highly structured objects which a random sample have essentially no chance of ever turning up. A proper test of the bound would e.g. start out wit ha single initial data and then try to modify it locally to make the inequality come closer to equality, then modify this data and repeat this until, possibly, data is found which violates the inequality.
In you discussion above about inequalities (1)+(2) you miss the alternative that in general none of the inequalities need to be true. There could very well be matrices for which (1) holds and others for which (2) holds.
For matrices in general let A be the 3X3-matrix which has a zero in the upper left corner and the rest is 0, and let B be the matrix which has -1 at the last two position in the top row the the rest is 0. If we let A be B_n and B be S_n then inequality (1) fails, since the norm of A-B is 3^0.5 > 2, and the norm of A is 1.
19 April, 2008 at 7:00 am
DavePurvance
Dear Jonas,
Your comments on my numerical tests may be true, and I would love to have the computing power to try what you describe as a valid numerical test. In fact, I hope someone is interested enough to try, even though if nothing negative is found, someone else might criticize that you just didn’t probe hard enough.
One point that I believe adds credence to my numerical tests is that both solution matrices
and bound matrices
are made from the same system matrices
, and,
are made from the same flow coefficients
. So , if
is random, so is the bound
, and if
is structured, so is
.
Also, I’m sure one can always produce simple 3×3 matrices that support any argument. The challenge is to make your argument with matrices that represent the periodic 3D incompressible Navier-Stokes equation.
21 April, 2008 at 10:56 am
285G, Lecture 7: Rescaling of Ricci flows and kappa-noncollapsing « What’s new
[...] to begin the surgery program discussed above. (In contrast, the main reason why questions such as Navier-Stokes global regularity are so difficult is that no controlled quantity which is both coercive and critical or subcritical is known.) The [...]
24 April, 2008 at 12:15 pm
Lower and Upper Bounds for Periodic Navier-Stokes Flows « Smooth Chaotic Fluids
[...] differential equation also has a general noncommutative time series solution. These all have been recognized by the world-renown expert Terry Tao (read in his blog’s comments). Terry also happens to be a very decent and fair person in my [...]
9 May, 2008 at 3:41 am
DavePurvance
Terence et al,
I seemed to have messed up the links.
My proposed solution to the periodic 3D Navier-Stokes problem will always be arXived here. My WordPress blog discussing this solution will always be here.
Your criticisms and comments are always welcome.
30 May, 2008 at 4:29 am
DavePurvance
Dear Terence,
Thanks for acknowledging in your 22 March 2008 comment above my commutative
and noncommutative
solutions to my periodic Navier-Stokes matrix differential equation
having initial condition
.
In the same comment, though, you misunderstood the pair of matrix polynomial inequalities
I used to put bounds
on the general noncommutative solution coefficients
. I take full responsibility for this misunderstanding because I didn’t really prove (2). My arXiv paper now does.
At the time of your comment I denoted in (2)
as
and in (3)
as
. Perhaps this notation was the source of the misunderstanding. Whatever the case, I am sure you understand that this is an important issue and I ask you for the opportunity to counter your criticism of (2).
The inequalities in (2) do not mean that a commutative flow
of one initial condition
somehow can control another noncommutative flow
of another initial condition
.
Matrix polynomials
are defined by
Bound coefficients
are the Taylor time expansion coefficients of the matrix exponential
assuming the same noncommutative
as in (4) and resulting in matrix polynomials
with
Again, as noted in (4) and (6),
and
are functions of the same
and therefore both are a function of the same initial condition
. The difference between
and
is that
is a solution to (1) because
assumes commutative
.
is not as solution to (1) because
assumes noncommutative
. They both, however, are stable functions of
independent of their commutative properties and
and
simplify to the same polynomial when
commute.
Inequalities (2) are generalizations of a pair of equivalent commutator inequalities, which for illustration, I will derive here for the two simple matrix monomials making up the matrix commutator
By the the definition of any matrix norm
Adding the two halves of (8) gives
Again by the definition of any matrix norm
Adding (9) and (10) and simplifying gives
Doubling both halves of (8) and subtracting the result of each from (11) gives
which upon rearrangement yield the desired inequalities
So actually, inequalities (13) are just two operations away from the well known matrix commutator inequalities (11). And again, the reasoning for deriving the Navier-Stokes (2) parallels the reasoning provided in (7)-(13).
Inequalities (13) can be easily checked numerically for arbitrary square matrices
. I used random complex matrices and MATLAB’s “gallery” function to mix and match different types of
. The experiments ranged in dimension from 2 to 2000 and the inequalities (13) never failed. And, of course, as detailed in my paper, I have also numerically checked the Navier-Stokes (2) through time order 10 using both random and structured initial flows. These checks have never failed either.
Sorry for the confusion and thanks for letting me clear it up.
5 June, 2008 at 7:39 pm
lib.mexmat.ru/forum
Тао http://terrytao.wordpress.com/2007/03/18/why-globa…..ment-26446 писал(а):
At present, all known methods for obtaining global smooth solutions to a (deterministic) nonlinear PDE Cauchy problem require either {В настоящее время, все известные методы получения глобальных гладких решений (детерминированных) нелинейных PDE (Cauchy проблема) требуют то или иное (также?)}
1. Exact and explicit solutions (or at least an exact, explicit transformation to a significantly simpler PDE or ODE); -{Точные и явные решения (или по крайней мере точное, явное преобразование к значительно более простому PDE или ОDЕ)}
2. ;Perturbative hypotheses (e.g. small data, data close to a special solution, or more generally a hypothesis which involves an somewhere); or
3. One or more globally controlled quantities (such as the total energy) which are both coercive and either critical or subcritica.- {Одно или более глобально управляемые количества, (типа полной энергии), которые являются принудительными и также критическими или подкритическими.}
Или, может быть, самого автора стоит об этом спросить, воспользовавшись окном для комментариев его блога по указанному в цитате адресу. Возможно, он, действительно, как и Вы, почему-то молчаливо подразумевает п. 1 в обязательном сочетании с п. 3, хотя из приведенной цитаты это, кажется, никак не вытекает. При этом заодно можно предложить Тао и его коллегам по обсуждению сформулировать свое отношение к высказываниям Л.Д. Ландау и М. Клайна, о чем мы здесь столько спорим.
12 June, 2008 at 9:23 pm
Shock waves « Hydrobates
[...] If viscosity is included in the description of fluids then the Euler equations are replaced by the Navier-Stokes equations. There is reason to suspect that in this case shock waves are smoothed out and a smooth initial configuration remains smooth in the course of the evolution, for all time. There are simple examples where this can be seen but there is still no global regularity result for Navier-Stokes (and no counterexample). The Clay Foundation has offered a prize of one million dollars for the solution of this problem in either direction. The fact that the prize has not yet been collected is a sign of the difficulty of the problem. For a discussion of this question and its broader mathematical significance I recommend the excellent account of Tao. [...]
16 July, 2008 at 7:28 am
Anonymous
I am looking for more information about the use of Strichartz estimates in solving wave equations. Can anyone provide me some reference? I am not looking for research papers yet. I want to grasp the main ideas. The info on DispersiveWiki is not particularly useful. And I don’t seem anything revelant except a lot of research papers which I want to avoid at the first stage. Thanks.
17 July, 2008 at 2:42 pm
Terence Tao
Dear Anonymous,
You might try Chris Sogge’s book on nonlinear wave equations, or my own on nonlinear dispersive equations.
17 July, 2008 at 4:50 pm
Anonymous
Thanks. I’ll have a look at both.
19 July, 2008 at 3:54 pm
Global existence and uniqueness results for weak solutions of the focusing mass-critical non-linear Schrödinger equation « What’s new
[...] smooth solutions to the Navier-Stokes equation is one of the Clay Millennium problems that I have blogged about before, but global existence of weak solutions is quite easy with today’s technology and was first [...]
22 July, 2008 at 11:19 am
Claes Johnson
A resolution of the Clay Navier-Stokes problem is proposed in the article
Blowup of Incompressible Euler Solutions, published online July 19 2008
in BIT Numerical Mathematics.
22 July, 2008 at 7:20 pm
Anonymous
To Claes Johnson
I don’t think this paper is called “resolution”… it is just some numerical evidence. For incompressible Euler, finite time blow-up or not, has been a long time controversial question. Both sides have numerical “evidence”…
For incompressible Navier-Stokes, physicists tend to believe no finite time blow-up.
23 July, 2008 at 12:21 am
Claes Johnson
To Anonymous:
This discussion is important. I hope e.g. Terence will take part.
You express a common misconception, addressed in the article, that
computational evidence is not mathematical evidence.
The central concepts are wellposedness and turbulence, and in this
context a computational solution is as much a solution as anything.
Please read the article with an open mind. Looking forward to
your comments.
Claes
29 July, 2008 at 6:34 pm
Anonymous
Claes,
The issue is that the NS equation is highly numerically unstable, so simulations showing blow up tell us very little; the blow up may be an artifact of the simulation rather than an indication that NS actually blows up.
[different anonymous from 7:20pm above]
29 July, 2008 at 11:21 pm
Claes Johnson
To Anonymous 2:
You are right that an unstable numerical scheme for a stable problem
can give artifacts. But our computational Euler solution is validated
by a posteriori output error estimation by Euler residuals multiplied
by stability factors obtained by solving dual linearized problems.
A computed Euler solution is thus a representative solution and
not an artifact. Since the computed solution shows blowup and
is representative, there is blowup.
Correspondence (with e.g. Terence) on the blowup problem is published on
http://www.nada.kth.se/~jhoffman/pmwiki/pmwiki.php?n=Forum.Clay
Claes
30 July, 2008 at 5:20 am
Anonymous
Claes, take a look at [[Goodstein's theorem]] in Wikipedia for how fast a function can grow and still be considered bounded in a mathematical sense. Double exponential is not bad at all. Triple or quadruple or 57-times-iterated exponential is only a little bit worse. Then you get functions like Ackermann’s, which are worse than n-times-iterated exponential for arbitrary large n, but for which you can still write down an formula in terms of induction on more than one variable. And then there are functions that grow so fast that no expression can even be written down–one can only prove indirectly that they are still finite. Goodstein’s theorem says that starting with any n, if you iterate a certain calculation for enough steps you’ll eventually reach zero. But it turns out (proved by Kirby and Paris) that as a function of n, the number of steps before you hit zero is so large, that it’s impossible to write down a formula for it. Yet that it’s always finite (the process never goes on without bound) is a remarkable fact of pure mathematics, even though even for values as small as n=4, f(n) could not be written down as a decimal number if you could engrave a billion digits on each electron in the physical universe.
A more practical example: suppose you have a computer program with array references like a[x*y + z]. It is provably impossible in general for a compiler to tell without running the program whether there will ever be a subscript-out-of-range error. Suppose instead the references are like a[37*y++19*z+3*w]–that is, you never multiply two variables together in a subscript, you can only add them or multiply by constants. Then compile time checking is possible (this is called Presberger arithmetic). The running time formula is a tower of iterated exponentials in the size of the expressions, which is absolutely intractable in the worst case, but it’s a remarkable theoretical result that it’s decidable at all. And it turns out that for most real practical cases, the worst case exponential tower is avoided. And even without that, the decidability is of interest to pure mathematicians. Maybe they’re not of interest to anyone else, but pure mathematicians went into that field precisely because they’re into that sort of thing. It’s not up to anyone else to talk them out of it.
Who knows, maybe global stability will involve growth functions like that. If so, the Clay math people want to know about it. In a physical sense it’s as bad as a singularity but from a pure math point of view it’s not the same at all. And as Terence and others have said, reformulations of the problem may be very interesting from a physics perspective, but the pure math question posed as a Clay problem is very specific.
30 July, 2008 at 5:27 am
Anonymous
Re the above: an amusing question occurs to me for the math folks here (I’m just a computer guy). What happens if someone turns in a proof of NS global stability, that is rigorous but depends on something like a large cardinal axiom? Is that even imaginable, mathematically speaking?
30 July, 2008 at 7:48 am
Claes Johnson
To Anonymous:
You touch an essential point: There is a big difference between
100 and 10^100 = googol. If you do not make this distinction
in a quantitative analysis, it is not an analysis. I am surprised
to see that some (pure) mathematicans do not make this distinction,
and I wonder what mathematical tradition it can reflect.
Claes
8 August, 2008 at 5:39 am
westy31
Hi turbulence friends,
I happened to be thinking about the Navier Stokes equation and turbulence, when I found this discussion. I have some questions and comments.
I recently put a free n-dimensional Navier Stokes simulator on the web at:
http://www.xs4all.nl/~westy31/CellFlow/CellFlow.html
One question I am working on is how turbulence scales in n dimensions. It was remarked by Terrence in this discussion that 3D turbulence is supercritical, while 2D turbulence is critical. This I would like to understand.
First, I will comment on scaling. I think the scaling

can be considered a special case of scaling behaviour that is most easily understood in terms of dimensional analyses. If you are given data about a fluid flow, without mention of the units used, you would not be able to deduce the correct units. This is because classical physics has a symmetry that corresponds to the fact that there are no natural length, time and mass scales. The scaling law for a quantity is easily written down by just writing out the dimension of the quantity:
etc,
The
scaling is a special case:
.
Deriving an energy spectrum from dimensional analyses only, was just what Kolmogorov did in his 1941 paper. The outcome, which is easy to check yourself, is

See also the Wikipedia page:
http://en.wikipedia.org/wiki/Turbulence
This outcome is independent of the dimension of the flow, (Note that k stands for the norm |k|, not the vector) which seems to contradict the claim that 3D and 2D have a different criticality. Personally, I have a hard time believing that the energy spectrum does not depend on space dimension, so I would welcome any arguments to that end. I heard (see Stephen’s comment) that in 2D, you have ‘Kraichnan power of [minus] 3 law’, but it is unclear to me how this relates to Kolmogorov’s theory. Clarification would be welcome here.
While thinking about the -5/3 law, I decided they need to be modified anyhow, to meet some consistency requirements. If you try to integrate the dissipation (
) over k, to check that you retrieve
, you run into a diverging integral. To solve this, I propose:

, dissipation takes over. With the ultraviolet exponential cut-off, we can integrate the dissipation and get a consistent answer.
Note that this formula is scale invariant.
The double exponential cut-off arises because for wave number higher than the Kolmogorov scale
Next, try to integrate E(k) to get the total turbulent kinetic energy. Again, we get infinity, but now because of the infrared limit. This is solved by an abrupt low wave number cut-off k_0, which is justified because we have a fluid in a finite container.
.
Interestingly, we can now for a dimensionless number out off the ratio between the high and low wave number cut-offs: Re=
This dimensionless number can be interpreted as the Reynolds number, if we insert the Kolmogorov velocity for velocity, and 1/k_0 as the scale in the regular expression for the Reynolds number.
The idea of this Reynolds number makes sense. After all, a laminar flow is quite smooth, and this should somehow be reflected in the energy spectrum formula. The formula I propose automatically makes a smooth (=steeply declining in wave number) energy spectrum if the Reynolds number is low.
Because of the viscous cut-off, I personally believe that the Navier Stokes equation does not blow-up, ie has zero energy in infinite wave numbers.
Gerard
23 August, 2008 at 6:48 am
PDEbeginner
Dear Prof. Tao,
This article is very nice and very helpful for me to understand the last term’s lecture on NSE I attended. As for the supercritical and critical things (I didn’t understand it in the lecture), I was wondering if my present understanding is right:
When the energy transfers to the finer and finer scales, the velocity will concentrate on smaller and smaller space. This will possibly make the critical and subcritical norm of velocity blow up.
On the other hand, some people works on partial regularity of the solutions, suitable solutions, and considers Hausdorff dimension of singular points. I was wondering if this is helpful for solving the problem.
Thanks a lot!
All the best,
PDEbeginner
24 August, 2008 at 9:15 am
Alberto G. P.
Hi,
The String 2008 conferences were held at CERN this week. The talk by Minwalla was about ‘Nonlinear Fluid Dynamics from Gravity’ and seem me very interesting:
(pdf here: http://indico.cern.ch/getFile.py/access?contribId=24&resId=0&materialId=slides&confId=21917)
I have read the Terence’s post and I think that he didn’t include this technique in the Strategy 3. Unfortunately I am not an expert on this subject, more exactly, I know very little about physics and mathematics. Nevertheless I am very interested in the following questions:
1) Could we know something about Navier Stokes equations (NSE) by means of the proposed duality between NSE and EFE (Einstein Field Equations)?.
2) We could discovery analytic EFE solutions in AdS space and then obtain his dual NSE solutions ( ?).
3) Is the global regularity problem for EFE resolved? If so it can be use the proposed duality to resolve the global regularity problem for NSE.
Can somebody help me? I asked about Minwalla’s talk in some physics’s blog but nobody reply me. Perhaps I am doing senseless questions (:-)sorry.
25 August, 2008 at 2:25 am
Pedro Lauridsen Ribeiro
Hi Alberto,
The aforementioned transparencies are a bit hard to read,
as some of the math symbols there are truncated. However,
I can make some comments on your questions proper.
First, let’s set the stage with a bit more of detail. Minwalla’s
proposed duality uses a particular instance of the
Maldacena-Witten AdS/CFT correspondence which relates a
particular class of solutions of Einstein’s equations with negative cosmological constant in d+1 dimensions which possess “locally” a conformal infinity of anti-de Sitter (AdS) type with solutions of
“generalized” Navier-Stokes equations in d dimensions. For that,
he uses the Fefferman-Graham asymptotic expansion of the
conformally rescaled bulk metric in a collar neighbourhood of
conformal infinity, whose coefficients satisfy, due to the Einstein
equations, a family of recursion relations which constitute a non-linear
Fuchsian system (i.e., a system of differential equations whose
coefficients present isolated algebraic singularities for one or
more coordinates, here the distance to the boundary in geodesic
normal coordinates for the conformally rescaled metric). Such
system determines uniquely the asymptotic expansion of the
bulk metric near infinity (hence, its _long-distance_ or
_long-wavelength_ behaviour) from the boundary metric and
a symmetric rank-2 tensor on the boundary called the
boundary stress-energy tensor. When conserved and traceless,
this tensor reduces to the rescaled electric part of the bulk
Weyl tensor, but this is not always the case – in general, the
divergence and the trace of the boundary stress-energy tensor
give rise to conformal invariants at the boundary which can often
be computed exactly, due to the Fefferman-Graham equations.
It’s the inhomogeneous conservation law of the boundary
stress-energy tensor that gives rise to the “generalized”
Navier-Stokes equation of Minwalla, for a particular class of
bulk metrics.
Now, about your questions (in some places without pretenses to rigour):
1.) The AdS-CFT correspondence exchanges large scales in the bulk
with small scales in the boundary, but knowing the behaviour of the
boundary conservation law in the large involves geometric information
deep in the bulk. The class of solutions of the bulk Einstein’s equations
which give rise to the “generalized” Navier-Stokes equations has
a small-codimension singular locus (a “black brane”) deep inside.
In this case, one may expect to obtain blowup for these “generalized” Navier-Stokes equations, but Minwalla presents only circunstantial
evidence to the conjecture that the “generalized” Navier-Stokes
equation always arises from the procedure above – this involves proving
global existence and uniqueness for the Fefferman-Graham system for
the bulk metric.
2.) Using the Fefferman-Graham formulae for the boundary stress-energy
tensor in terms of the bulk metric, one can obtain the actual form of the coefficients of the “generalized” Navier-Stokes equations, but not its actual solution, unless the form of the resulting equations is, of course,
sufficiently simplified. Otherwise, the equation is as hard to solve as
ever.
3.) When the boundary metric and the boundary stress-energy tensor
are real-analytic, the Fefferman-Graham system admits a unique
solution and, moreover, the asymptotic expansion actually converges
in a sufficiently small collar neighbourhood of the conformal infinity
(this expansion is a power series in the geodesic distance to the
boundary if and only if the boundary stress-energy tensor is conserved
and traceless, otherwise the expansion involves logarithmic terms as
well). This was proved by Kichenassamy. As far as I know, almost nothing
is known in general about the maximum domain of convergence.
For non-analytic boundary data, there seems to be so far no sufficiently
strong estimate of the solution in terms of the boundary data. This
question was studied by Michael T. Anderson in several papers,
almost all of them in Euclidean signature, because then one can use
techniques from nonlinear boundary value problems (Leray-Schauder
degree, etc.), but the problem is quite different from a PDE perspective,
and Wick rotation (to move back and forth from Lorentzian to Euclidean signature) is a very tricky procedure in curved spacetimes, which globally
in general cannot be done, specially if the metric is not stationary.
Only very few of Anderson’s papers deal with Lorentzian signature.
In one of these, he shows that a sequence of Cauchy data for
geodesically complete, asymptotically AdS spacetimes (i.e., subject
to the boundary condition that the conformal infinity should be that of
AdS, or, more precisely, of its universal cover) which become asymptotically
stationary at infinite times actually give rise to a _globally_ stationary
spacetime. This is rather different from the asymptotically flat case,
since nonstationary (albeit sufficiently small) perturbations of the
Cauchy data for Minkowski spacetime give rise to nontrivial, smooth
and geodesically complete spacetimes, as proven by Christodoulou
and Klainerman. The difference is that the boundary conditions needed
for the Einstein equations in the case of a negative cosmological constant
(as here the “unperturbed” solution – i.e. AdS – is not globally hyperbolic)
do not let the nonlinear effects disperse at large times – they are focused
back in an “almost-perodic” manner. The result of Anderson gives evidence
to the possibility that AdS spacetime is globally nonlinearly unstable – according to his result, nonstationary perturbations will eventually lead
to geodesic incompleteness. But here there’s a caveat: to prove his result,
Anderson needs a “unique continuation” property for the rescaled linearized Einstein equations across the boundary. Here he borrows his intuition from the Euclidean case, where such property holds even in the nonanalytic case. However, the linearized Einstein equations are of wave type (modulo
some gauge fixing), and in the nonanalytic case (which is even beyond the
reach of Kichanassamy’s approach), there are counterexamples to
unique continuation by employing nonanalytic perturbations of
the coefficients of lower order (Cohen, Hörmander, Alinhac-Baouendi)
if, in particular, the boundary is totally geodesic (i.e., gliding – null geodesics tangent to the boundary always have infinite order contact), which is
_always_ the case for conformal infinities. A proof of global nonlinear
instability of AdS would probably involve proving that such perturbations always lead to geodesically incomplete bulk metrics, which might be the
case for all we know, but it remains an open (albeit perhaps solvable in
this form) problem. Now, to move from this back to Minwalla’s “generalized”
Navier-Stokes equations, one needs estimates relating boundary quantities to bulk quantities (which is also needed if one wants to extend
Kichenassamy’s existence result to the nonanalytic case and prove
on the other hand the _stablility_ of AdS, which does hold locally in time in
4 dimensions, as proven by Friedrich). Such estimates exist in the
case of Euclidean signature, but it’s hardly likely they carry through
to Lorentzian signature, specially in nonanalytic cases. Since the boundary stress-energy tensor involves the Weyl tensor, and so do the estimates obtained by Christodoulou and Klainerman in their proof of the global nonlinear stability of Minkowski spacetime, there may be such estimates
in spite of that.
You see, it’s not that your questions are senseless, it’s just that they are
too hard, but nevertheless always worth asking…;-)
Cheers,
Pedro
25 August, 2008 at 8:31 am
Pedro Lauridsen Ribeiro
Ah, yet on Alberto’s third question: for large asymptotically
flat data, Christodoulou recently proved that, under some
reasonable technical assumptions, solutions of Einstein’s
equations develop trapped surfaces in finite time, which
implies gravitational collapse by Penrose’s singularity
theorem. For those brave souls who wish to endeavor
going through the proof (594 pages!!!), the preprint
can be found at http://arxiv.org/abs/0805.3880.
25 August, 2008 at 9:47 am
Alberto G. P.
Let me make some comments about I have understood
We can start with a bulk metric that define an AdS space. The bulk metric is an exact solution of EFE with negative constant. There are an unique asintotic expansion of the bulk metric valid only in a certain region belongs to AdS (near infinity). The coefficients of this expansion are used to buid a new tensor defined on the boundary of AdS (stress-energy tensor). If it is imposed some conservation law to this tensor then the generalized NSE arise. But although we know the stress-energy tensor we don´t know dual generalized NSE solutions (i.e.it hasn´t got the solutions yet).
The global regularity for EFE is an open problem. Moreover It can not be used the AdS/CFT correspondence to solve global regularity for generalized NSE because nobody know how to transform the bulk quantities to boundary quantities.
As thanks let me conclude with a quotation from Don Quixote:
“nunca fuera caballero
de damas tan bien servido
como fuera Don Quijote
cuando de su aldea vino
doncellas curaban dél;
princesas de su rocino,”
;-).
25 August, 2008 at 5:08 pm
Terence Tao
Dear PDEBeginner,
You are correct that the difficulty in the Navier-Stokes regularity problem is in preventing passage of energy to higher scales, which will cause subcritical norms to blow up (and critical norms to either blow up, concentrate, or otherwise develop a singularity). Thus far, the partial regularity results that are known do not prevent this blowup from occurring. To oversimplify, the above blowup scenario is consistent with concentration on any set of (parabolic) dimension 1 or less, but not on dimension higher than 1. And indeed, the known partial regularity theory, in particular the result of Kohn, Caffarelli, and Nirenberg excludes singularities with dimension higher than 1. (I am not sure if singularities of dimension exactly 1 are known to be excluded. Given that in two dimensions no singularities occur, this suggests that 1-dimensional singularities do not occur in three dimensions, though this is nowhere near a rigorous argument.)
26 August, 2008 at 2:39 pm
westy31
Hi again,
Question:
Would the Energy spectrum contain enough information to decide if the equation blows up?
For example, say the energy spectrum is k^-5/3. This would mean, I presume, that the velocity does not fluctuate infinitely fast in space, because high frequancy components of the field get progressively smaller. Not so for higher derivatives of the field though, since these get multiplied by k each time you differentiate.
As a second example, take k^-5/3*exp(-k/k_c).
This one stays small for large k, for all derivatives. It would seem to me that this guarantees smoothness.
Gerard
26 August, 2008 at 3:50 pm
Terence Tao
Dear Gerard,
Yes, smoothness of the velocity field is equivalent to the energy spectrum decaying faster than any power of k as
. So if the Kolmogorov
power law managed to extend itself to arbitrarily high frequencies then regularity of Navier-Stokes would break down; but I believe the expectation is that the viscosity term in Navier-Stokes should cause one to leave the Kolmogorov regime at sufficiently fine scales. (The derivation of the Kolmogorov law implicitly assumes that energy is spread evenly throughout space at each frequency scale, and that the rate of energy flows remain constant, whereas it is known that blowup only occurs if the energy concentrates on a low dimensional set, and if the rate of energy flow becomes infinite. So blowup is not expected to be caused by the same mechanism that gives the Kolmogorov law, or at least the two phenomena should occur at different scales.)
28 August, 2008 at 1:26 pm
On the Convergence of Periodic Navier-Stokes Flows « Smooth Chaotic Navier-Stokes Flows
[...] On the Convergence of Periodic Navier-Stokes Flows Filed under: Mathematics — DavePurvance @ 2:26 pm It’s not that the 3D spatially periodic, incompressible Navier-Stokes equation can be posed as a nonlinear matrix differential equation. It’s not that this matrix differential equation has a smooth and bounded matrix exponential solution, found when the flow in this differential equation is Taylor expanded in time and the resulting system matrices are all commutative. And, it’s not that this time expanded matrix differential equation also has a more general time series solution. It is called the noncommutative solution because its coefficients involve noncommutative matrix polynomials. This matrix differential equation and its commutative and noncommutative solutions have been recognized by the world-renown expert Terence Tao. Read Terence’s 22 March 2008 comments under his blog Why global regularity for the Navier-Stokes is so hard. [...]
28 August, 2008 at 10:26 pm
Claes Johnson
Dear Gerard:
It is necessary to give a quantitative meaning “smooth”, which means that the size of Sobolev norms must be taken into account. If the C^1 norm of a function is a number of moderate size, say 10^2, you can say that the function is C^1-smooth, but if the C^1-norm is of size 10^100 = googol, then it is not C^1-smooth. The reason to introduce Sobolev spaces with norms is to measure norms, and to measure norms is to make a distinction between 10^2
and 10^100. It is very strange that this aspect seems to be completely lost in
the discussion of the Clay problem. In order to answer if a smooth solution exists, you first have to define what “smooth” is. This is not done now and no progress is being made. So Terence: what is your definition o “smooth”? Is the size of Sobolev norms not relevant? If not, how do you make a distinction between a smooth and a non-smooth function? What is then the meaning of Sobolev norms?
15 September, 2008 at 10:23 am
Bill Layton
Hi All,
This is a very interesting discussion that I’ve enjoyed reading!
Might helicity play a key role in NSE uniqueness / non-uniqueness and deserves much more study?
It is the other key invariant of the 3d Euler equations and a control quantity for the 3d NSE. Of course it is much harder than enstrophy in 2d because it has 2 signs and dissipation reduces helicity mode by mode rather than by global magnitude.
The plan is one that surely many have thought about. If the NSE nonlinearity is zero (or sufficiently small) uniqueness is trivial. If the helicity is zero or sufficiently small, then the flow is also unique and smooth (a recent result of Berselli). The “only” problem is then to interpolate between these two cases locally.
This is surely a technically intricate problem. But, to me it seems like one with more hope of success for the first few steps than squeezing a “new” a priori bound on ||grad u|| from what are in many ways essentially the same estimates that Leray derived in his seminar papers.
From the computational point of view (and computational insight is very useful since it’s generally easier to prove something that is true and understood at some level), helicity is a scalar so much easier to visualize.
So, …
Any ideas on this path among the many possibilities??
Bill
26 September, 2008 at 5:09 am
David Collins
Dear Terence,
I am a physics layman and math undergraduate from germany on a ODE level. I have a bulk of wild ideas, some probably persued before and many false. Some at least seem to fit to some problems outlined above.
Smooth … ? … Discrete
a) Instead of only using discrete models in order to find smooth limit tactics, maybe Navier Stokes equations lead in the other direction to a ‘self-quantization’. A superposition of data ensembles might ‘wash out’ the many probable topological knots of vortex filaments (’Shaking Strings in a box leads to Knot formation’ http://www.pnas.org/cgi/content/abstract/104/42/16432 ) (long shots 1 & 5)
b) Velocity field modeled as a rapid growth phenomena.
A toy model which would be less rigid and would allow a bit more complexity than triangular or cubic lattices, and that is appears frequently in nature could be a recently further generalized model of primordia growth, the Snow map http://maven.smith.edu/~phyllo/Assets/pdf/snow.pdf which allows continuous diffeomorphism via irregular rhombic tilings (which are periodic weak attractors) between disjunct discrete lattices (strong asymptotic attractors). In it’s current form it is of course only an analog to laminar flow of molecule ensembles. But there is a striking interplay between random divergence angles in local topology & pseudorandomness in global topology. What happens when taking the limit in this model? (long shots 2, 5, 6)
c) With the unproved assumption that nature follows a ‘least crypticity princible’ and thus the simplest mathematics possible one could heuristicly ask ‘If a smooth continuous solution to fluid dynamics does exist, why does nature bother with quanta (x > 0)?’
d) Maybe the mathematical superposition of all possible solutions lead to temporal and spacial aliasing effects, thus leading to quasilattice geometry. Thus blowups or quasi-blowups in Ackermann/Goodstein scenarios. (Very wild guess)
Dimensional analysis shows that fluid entities have analogs in Plancks Laws:
Dynamic Viscosity n (eta)[kg/ms] ~ Spectral Energy Density u
Surface Tension s (sigma) [kg/s²] ~ Spectral Radiance per Frequency I(f)
Pressure Change dP/dt [kg/ms²] ~ Spectral Radiance per Wavelength I(r)
Dimensional analysis also shows possible uncertainties or complements:
Dyn. Viscosity & Volume: n * V > h
Kinematic Visc. & Mass v * m > h
(measuring Volume of superfluid H3 with exact viscosity n = 0 ?)
Again a very wild guess:
A ‘bubble’ or more precisely a ‘Spherical Hill Vortex’ might be the missing link between ‘Hollow black body’ and ‘Fluid’, having n, s, dP/dt aka u, I(f), I(r).
h) A vector that keeps appearing in this dimensional analysis and in Navier Stokes is w := 1/rt or [1/ms]. For instance h=mw and Du = w in the viscose term.
i) First attempt to interpretate w: This could be interpreted either as rate of curvature change (dK/dt with K=1/r) or as change of spacial frequency. A 3-d Snow model might show which and when is the appropriate interpretation via the parameters D (spacial frequency, periodic attractors), the ‘crookedness’ of the primordial front (local discrete curvature) i.e. the divergence angles and/or the mean parastechy angle (global curvature). Actually a common formulation of NS has the (‘fluid kinematic energy’) dimensions
F/V = mw² (Navier Stokes)
n = mw (Dyn. Viscosity as ‘fluid impuls’)
U ~ mw (Spect. Radiance per frequency)
P= mw/dt (Pressure as ‘fluid force’) etc.
j) Second attempt: After asking my calculus professor about the nature of this ‘dual’ vector, she gave me the hint to look up the Poisson resummation formula which relates a summation over the numbers x_n = n, to a summation over the numbers y_n = 1/n .
This in turn can be generalized via Pontryagin duality to dual lattices, and more general via Selberg trace leads directly into the heart of pseudo/randomness, Riemanns Zeta hypothesis.
Thanks for your time and patience
Dave
28 October, 2008 at 9:53 am
Turbulence « Unruled Notebook
[...] If you can handle comfortably some mathematical language, read the excellent post titled “Why Global Regularity for Navier Stokes is Hard” by Terrence [...]
15 November, 2008 at 4:29 am
isaac
Dear Terence,
What do you do during your free time? Any hobbies?
2 January, 2009 at 12:53 pm
Liangyu
Hi, terry, if the solution of the Navier Stokes equation can develop sigularity in finite time then the Navier Stokes equation surely cannot be an accurate model of 3D fluid dynamics. Which would implie that some physical assumption in the formulation of the Navier-Stokes equation is wrong?
3 January, 2009 at 10:32 pm
Terence Tao
Dear Liangyu,
The derivation of the Navier-Stokes equations from more fundamental laws of physics (e.g. Newton’s laws of motion) involves a number of simplifying assumptions (most notably, treating a fluid as a continuum rather than as consisting of a large number of atoms). Presumably, if these equations lead to singularity, then these assumptions would break down before the actual singularity is reached. (Note also that it is known that singularity can only occur if the velocity goes to infinity; in practice, of course, relativistic effects (among other things) would kick in long before then.)
26 January, 2009 at 7:34 am
Michael Nielsen » Doing science online
[...] information and insight. To understand how valuable Tao’s blog is, let’s look at a example post, about the Navier-Stokes equations. As many of you know, these are the standard equations used by [...]
27 January, 2009 at 7:14 pm
km
Hi,
In the following link:
http://knol.google.com/k/claes-johnson/the-clay-navier-stokes-millennium/yvfu3xg7d7wt/14#
a comment says that a blowup solution in Navier-Stokes means turbulence and “A turbulent solution is a non-smooth solution, and blowup from smooth data initial data is the same as transition from laminar to turbulent flow”.
So, if the above is true, then Navier-Stokes can be used to describe turbulence only if it does give non-smooth solution. Is that right?
27 January, 2009 at 10:11 pm
tmr
Is it known for blow up to occur in 4 dimensional Navier-Stokes?
28 January, 2009 at 1:10 am
Mathematics, Science, and Blogs « Combinatorics and more
[...] and blogs. Michael’s primary example is a post over Terry Tao’s blog about the Navier-Stokes equation and he suggests blogs as a way of scaling up scientific conversation. Michael is writing a book [...]
28 January, 2009 at 9:29 am
Terence Tao
Dear km: Unlike blowup (which means that the solution ceases to be smooth after a finite amount of time), which is a precise mathematical concept, turbulence (which means that energy shifts to higher and higher frequencies, or equivalently to finer and finer scales) does not have a canonical mathematical definition. But, at an intuitive level at least (and oversimplifying somewhat), one can think of blowup as an infinite amount of turbulence: a significant portion of the energy moves into infinitely fine scales within a finite amount of time, causing singularity. So if there is no blowup for the Navier-Stokes equation, it means that the solution can become somewhat turbulent for a finite amount of time (as is of course seen in real-life fluids), but eventually the viscosity effects dominate and damp out the turbulence. If there was blowup, then the mathematical Navier-Stokes equation would become singular after some time, though in the physical setting what would happen instead is that the solution becomes so turbulent at so fine scales that the simplifying assumptions used to derive the Navier-Stokes equation (e.g. assuming that a fluid is a continuum, rather than being made of molecues) break down, and some other dynamics take over. This does not seem to occur in practice (fluids do become turbulent, but not all the way to the molecular scale; at some point, viscosity effects assert themselves to remove the turbulent energy from the system). Whether this is always the case is, literally, a million-dollar question.
Dear tmr: No blowup result is known for any higher-dimensional Navier-Stokes equation as far as I am aware. Heuristically, the higher-dimensional equations should be more unstable and so would have a better chance of exhibiting blowup, but by the same token the instability may make it even harder to rigorously demonstrate blowup in higher dimensions than in lower ones. But perhaps some low-dimensional reduction of a high-dimensional Navier-Stokes equation might be amenable to study. (In the periodic setting, one can view a low-dimensional Navier-Stokes solution as a solution of a high-dimensional Navier-Stokes equation also, simply by extending the solution trivially in some additional spatial dimensions, so in that case, at least, global regularity of high-dimensional NS implies the same for low-dmiensional NS.)
28 January, 2009 at 4:05 pm
km
Hi Terence,
Thanks for the explanation. So, reality suggests that Navier-Stokes equation should not give blow-up solutions if we believe Navier-Stokes equation describes it accurately.
29 January, 2009 at 3:51 am
big bangs spectator
The Laplacian is physically problematic, however, that we haven’t observed a blowup doesn’t mean it can’t occur. E.g., there have been reports of strange new experiments, e.g., a hexagonal pattern in a rotating bucket of fluid. Ps a quick off-topic perspective statistic: the Physical Review Series database gives 4844, 3271, 1997, 299, 804, 564, 44 citations with the words Schrodinger, Josephson, Ginzburg-Landau, Landau-Ginzburg, sine-Gordon, Navier-Stokes, KdV in the Title/Abstract.
29 January, 2009 at 4:32 am
km
First, I have to say that I am an novice in turbulence or Navier-Stokes equation (NS equation).
Second, I think it is impossible to observe blowup or infinite amount of turbulence in real life. An object in real life with singular properties is probably the black-hole. I can think of no other physical phenomena in reality that possess singularities like infinite velocity or energy.
Third, your comment suggests a possibility that not observing a blowup or infinite amount of turbulence in real life doesn’t mean NS equation can’t give blowup solution.
29 January, 2009 at 11:17 am
tmr
Thanks a lot prof Tao for your reply!
11 February, 2009 at 5:22 am
naresh
heloo frnds
i think everyone over here are sticking themselves to naviers stoke equation , now lets think different, is navier stokes equation true????????????? i dont think so specially for a 3d case, reason is
as far as my knowledge, in an flow over an aircraft wing huge turbulance exist, first of all we need to understand how viscosity varies in 3 dimension with speed, below mach 1 its a very simple case we can say viscosity is proportional to speed but after mach 1 shock waves begin to form so here comes the case , hence if we understand the variation of viscosity with speed keeping the surface in mind then turbulance can be solved,
naviers stokes equation is wrong because!!!!!!!!!!!!!!!!!!!!
it relates velocity with pressure only what about temperature?????
for god sake viscosity of air changes with temperature hence we need to consider it/// please think about this hence we shldnt find the grad of kinematic viscosity instead we should use the variation of temperature with space coordinates in terms of viscosity, next thing is co effcient of viscosity must be considered , but it is nt done in navier stoke equation
so my opinion is navier stokes equation is absolutely wrong and hence it is unprooved
frnds and professors ive got the solution for this problem but i dnt know how to move about pls help me
thanks a lot for sparing your precious time sin reading my opinion
thankyou and please reply me ur questions and comments
im eager to hear from you
thanks again
7 March, 2009 at 7:51 pm
Tarun
Thanks for the great run down, Professor Tao.
I was reading about the function space BMO (Bounded Mean Oscillation) when I came across a paper which said that the regularity of BMO is ideal for Navier-Stokes solutions. I was wondering why this was and what role does (or maybe, should) Fefferman duality (
) play in the regularity of Navier-Stokes solutions?
Thanks a lot,
Tarun
7 March, 2009 at 8:52 pm
Terence Tao
Dear Tarun,
The space BMO (or more precisely the space
) is generally believed to be the weakest function space in which one can plausibly hope to (locally) solve the Navier-Stokes equation by perturbation theory methods, at least without first applying some further renormalisation of the equation. Basically, in order for the equation to make sense even as a distribution, the velocity field u should be locally square-integrable in spacetime; and in particular the effect of the nonlinearity
(viewed as a forcing term) should have this square-integrability. If one follows the perturbation theory philosophy of assuming that the nonlinear solution u behaves like its linear counterpart
, and then one inspects the characterisation of BMO functions in terms of Carleson measure properties of their heat extensions, we see that what one is asking for is essentially that
lies in BMO.
This seems to not be too directly related to the Fefferman-Stein duality between BMO and H^1, although of course the Carleson measure characterisation of BMO dualises to give a version of the atomic decomposition for H^1.
20 March, 2009 at 12:28 am
Student
Dear Professor Tao:
I am a student studying stochastic differential equation and am a novice with physics or applied math. Some elaboration on your sentence below would be very much appreciated:
“Things look better on the stochastic level, in which the laws of thermodynamics might play a role, but the Navier-Stokes problem, as defined by the Clay institute, is deterministic, and so we have Maxwell’s demon to contend with.)”
Dear Tarun:
I would be interested to know the paper which argued that the space of BMO is ideal for NSE solution.
Thank you.
21 March, 2009 at 10:31 pm
tmr
Probably that is a Koch-Tataru paper.
22 March, 2009 at 3:07 pm
Student
Dear tmr:
Thanks. That’s kind of what I thought, but I was wondering if there are any others that can be recommended.
26 March, 2009 at 9:27 am
Anonymous
Dear Terry,
Did you used to have notes on the Nash-Moser iteration scheme somewhere? I looked on your website and could not find it. Thanks!
29 March, 2009 at 1:37 pm
Terence Tao
Dear anonymous,
The Nash-Moser notes are at http://www.math.ucla.edu/%7Etao/preprints/Expository/nashmoser.dvi
29 March, 2009 at 8:52 pm
Anonymous
Dear Prof. Tao,
Excuse me my question is not related to the NS equations but since I do not know under which post we should write our general questions, I am writing it here.
My question is that is it possible to have two dependent random variables such that one has binomial distribution and the other has normal distribution?
thanks
17 April, 2009 at 5:51 pm
Student
I would really appreciate if you could explain why
“Riesz transforms are bounded in BMO, but not in L^oo”
a claim made by “Limiting Case of the Sobolev Inequality in BMO, with Application to the Euler Equations” by Kozono and Taniuchi.
Best regards,
17 April, 2009 at 6:01 pm
Terence Tao
Dear Student,
If you test the Riesz transform R_j against the indicator function
, you will obtain a logarithmic divergence. (This is easiest to see in one dimension, i.e. applying the Hilbert transform to the Heaviside function.) So
is not bounded on
. The bound on BMO is more complicated to show, and basically requires the use of Calderon-Zygmund theory; I would suggest Stein’s “Harmonic analysis” as a reference.
18 April, 2009 at 12:05 pm
anonymous
I have a question about strategy 2, discovering globally controlled quantities. Controlled quantities seem to be connected to invariance of the equations. For example, with the (linear) Schrodinger equation a phase rotation corresponds to the conserved quantity of mass, a space translation corresponds to momentum, etc. However, these invariances are precisely the defects to compactness of the linear Schrodinger operator, in the sense of the profile decomposition. In fact, one can think of the profile decomposition as the statement that the known invariance of the Schrodinger operator are the only (at least, noncompact) ones. It seems reasonable that we might be able to reverse this logic, and conclude that the known conserved quantities are the only ones. I don’t know very much about Navier-Stokes, but would it be plausible to attempt to classify the conserved quantities in a similar manner?
18 April, 2009 at 1:47 pm
Terence Tao
Dear anonymous,
For Hamiltonian equations (such as the Schrodinger equation), Noether’s theorem provides a very close connection between conserved quantities and invariances of the equation. But conservation laws aren’t necessarily the only ways to get control, and in any event Navier-Stokes is not a Hamiltonian system (it is not time-reversible).
Another obvious source of controlled quantities are monotonicity formulas – a formula that demonstrates that some coercive quantity is non-increasing in time. For instance, energy is a monotone quantity in Navier-Stokes. But unlike the situation with conservation laws of Hamiltonian systems, where we have a nice characterisation, there is no known general way to systematically classify all monotone quantities associated with an equation such as Navier-Stokes. But if one could find more such monotone quantities (and if they had a favourable scaling), then this would be major progress on the Navier-Stokes problem. (The solution to the Poincare conjecture, incidentally, was only made possible by the discovery by Perelman of certain scale-invariant monotone quantities for Ricci flow; one could conceivably imagine a similar miracle leading to a breakthrough for Navier-Stokes, though given that this equation is much less geometric in nature than Ricci flow, such miracles are perhaps less likely.)
24 April, 2009 at 11:29 pm
Student
Dear Professor Tao:
Thank you very much for your reply above. I have the follow-up questions.
I understand that this improvement below was made (Sorry for not using LaTeX; I’ll retry if it’s not acceptable):
From Beale-Kato-Majda (1984):
||grad u||_{L inf}
<=
C(1 + ||vort u||_{L inf}(1 + log^{+}||u||_{W^{s+1,p}) + ||vort u||_{L^2})
To Kozono, Taniuchi (2000):
For 1 < p n/p, with C(n, s, p)
||f||_{L inf}
<=
C(1 + ||f||_{BMO}(1 + log^{+}||f||_{W^{s,p}}))
In words of Ohkitani Koji (2008), that is,
“At the level of the BKM theorem which deals wit ha sup norm of the vorticity, it was not possible to rule out a possibility that vorticity blows up mildly as a logarithmic function of space variables. According to the recent updates using a BMO norm, we can safely exclude such possibilities”
So what is the next step? Unfortunately, Kozono and Taniuchi did not really leave any thought on the subsequent possible advancement. Looking back, it seems, at least to my Professor, BMO was the obvious choice. But he thinks it’s hard to see any more improvement possible beyond.
Elaboration would be very much appreciated. I will reread Tataru’s paper you mentioned; it seems to be connected somewhat.
Finally, please share the source of the paper you are talking about, in reply to a comment by Katz below:
“Thanks to the recent work of Escauriaza-Seregin-Sverak and others, we know that blowup solutions to Navier-Stokes must in fact blow up in the critical norm L^3({\Bbb R}^3). “
6 May, 2009 at 9:07 pm
Student
“Try a topological method. This is a special case of (1). It may well be that a primarily topological argument may be used either to construct solutions, or to establish blowup; there are some precedents for this type of construction in elliptic theory.”
Please give us your specific recommended references (books, articles) to examples of these “precedents.”
Best regards,
7 May, 2009 at 3:51 am
Terence Tao
Dear Student,
One example of such methods that comes to mind is Gromov’s topological method to construct J-holomorphic curves in his celebrated paper on non-squeezing. There are also a number of methods for elliptic variational problems with a somewhat topological flavour (e.g. the mountain pass lemma).
18 May, 2009 at 9:56 pm
NSE Beginner
Dear Professor Tao:
1. Aside from using Leray-Projection operator, one way to get rid of the pressure is to work on the vorticity equation. Solving this vorticity equation should be equivalent to solving the Navier-Stokes Equation. Is that correct? And if so, any pros and cons for this idea?
2. Beale-Kato-Majda showed in 1984 that if the vorticity does not blow up, then the solution can be extended beyond time T. It seems obvious to me that if the gradient, first derivative, of the vorticity does not blow up, then the same result applies (I am asking this question because showing that the gradient of the vorticity does not blow up seems easier from looking at the vorticity equation). Is that correct?
Thank you.
19 May, 2009 at 4:16 pm
Terence Tao
Dear NSE beginner,
Well, the vorticity equation still involves the velocity u, which requires the Biot-Savart law to recover from the vorticity, so one is pretty much back to where one started. But it does seem clear that the dynamics of the vorticity (and in particular, the behaviour of vortex stream lines, vorticity shearing, etc.) is going to be an important aspect of the problem.
At fine scales (which is where all the difficulties to the problem lie), the derivative of vorticity is going to be much larger than the vorticity itself, so it is extremely likely that the BKM result extends to the derivative. By the same token, though, this derivative is going to be harder to control than the raw vorticity; in particular it will be sensitive to shearing effects caused by irregularity in the velocity field u, which as mentioned earlier cannot be so easily eliminated from the vorticity equation.
23 May, 2009 at 9:15 am
Student
Dear Professor Tao:
There was nowhere else for me to write this; at least there was somebody else who asked you about the Nash-Moser Notes on this blog.
Below I am writing about your Nash-Moser Notes
http://www.math.ucla.edu/%7Etao/preprints/Expository/nashmoser.dvi
On page 3, in the middle you state
“we conclude from (5) (with x = u_n …”
but I don’t see any x in (5). You mean u = u_n?
On page 4 bottom, (11), I think the power of N should be 40, not 20.
In the middle page 6, you define the Littlewood-Paley operators, the case for N. I think there is one extra f^ that you do not need, for the case of > N that is.
Finally, my question: on the bottom of page 6 you claim a certain identity “Suppose for now that u_1, . . . , u_{n+1} \in 5B. Then from the identity…”
I do not know how you obtain the third equality; i.e. going from second line to third. Please explain.
23 May, 2009 at 10:32 am
Terence Tao
Thanks for the corrections! In the application of (5), x and h should be u and v (similarly the h in (5) should also be a v).
On page 6, on the third line of the equality,
(you can see this also by looking at the next displayed inequality). The point is that
is equal to
.
26 May, 2009 at 2:31 pm
PDEbeginner
Dear Prof. Tao,
I read your nice note on the Nash-Moser iteration scheme http://www.math.ucla.edu/%7Etao/preprints/Expository/nashmoser.dvi as being shown above. But it is a little abstract, could you give some nice literatures having some examples? (I tried to find some examples in your lecture notes on Nonlinear Dispersive and Wave Equations at New Mexico State University, but failed).
Thanks in advance!
26 May, 2009 at 3:55 pm
timur
Dear PDEbeginner,
To me, the notes you mentioned was the most readable account on Nash-Moser, the others being very complicated, abstract, long, or the main idea was buried in details. If you think you lack intuition on what is “frequency” of a function you might want to look at first couple of the following notes (which are again by prof. Tao):
http://www.math.ucla.edu/~tao/254a.1.01w/
26 May, 2009 at 11:33 pm
PDEbeginner
Dear Timur,
Thanks a lot for your kindness and your reference! I have downloaded part of the lecture notes provided by you and plan to read them (especially on Littlewood-Paley projection part).
27 May, 2009 at 9:22 am
Anonymous
Hamilton’s paper, referenced in Tao’s notes, has a number of examples of how the Nash-Moser inverse function theorem can be used.
3 June, 2009 at 8:09 am
hegel triad
Hi Prof. Tao,
I would like your advice on the following. I am at a conference in France, and the famous Prof. Temam has just shared with his audience a model used in oceanography – the Lagrangian-averaged Navier-Stokes-alpha. He said that it regularizes the NS by modifying the nonlinearity so that scales smaller than alpha are swept by larger scales, while preserving their fundamental inviscid transport properties such as convection of vorticity. I have misgivings as the audience was shell-shocked. Please express your opinion on the model.
4 June, 2009 at 12:53 pm
The Vlasov-Poisson system « Hydrobates
[...] of the equations. For a discussion of the significance of scaling properties in general see Terry Tao’s post on the Navier-Stokes equations. In the case the potential and kinetic energies satisfy an inequality of the form and this plays [...]
13 June, 2009 at 3:05 am
Student
Dear Professor Tao:
You stated:
“There are a couple of loopholes that one might try to exploit: one can instead try to refine the control on the “waiting time” or “amount of mixing” between each shift to the next finer scale, or try to exploit the fact that each such shift requires a certain amount of energy dissipation.”
I’m sorry but I do not have a good grasp of what you mean here. It seems as if there is some periodic finer rescaling from time to time and shift between each. I can see why if there is such thing, it would be a problem but how do you know that there IS such periodic rescaling and how should it be described mathematically? Or are you suggesting to suppose that there is such periodic rescaling and try to overcome it?
And why does this shift require any energy dissipation? I think of the Laplacian in the NSE when you say dissipation, but I’m afraid that’s not what you mean.
Any elaboration would be appreciated.
13 June, 2009 at 7:35 am
Terence Tao
Well, if I if I could prove that a solution could periodically shift its scales to higher scales indefinitely, then I would have disproven the regularity conjecture; my point is that this is an “enemy” that one will have to address at some point if one is to solve the conjecture. I also don’t know whether the dissipation effect of the Laplacian will be strong enough to prevent an infinite number of such shifts, but this is one potential way to stop such a cascade from occurring. I used some ideas distantly related to this when establishing global regularity for the logarithmically supercritical wave equation, but this equation was only barely supercritical, and is unlikely to be adapted as is for the Navier-Stokes equation (though perhaps some sort of “logarithmically supercritical Navier-Stokes equation” is within reach by these sorts of methods).
19 June, 2009 at 9:08 pm
Global regularity for a logarithmically supercritical hyperdissipative Navier-Stokes equation « What’s new
[...] supercritical and thus establishing global regularity beyond the reach of most known methods (see my earlier blog post for more [...]