Navier-Stokes方程的弱解、强解与光滑解：正则性的差距

情况是这样的。千禧年大奖问题悬赏一百万美元,用于解决三维Navier-Stokes方程是否总是存在永久持续的光滑解,还是可能发生爆破。光滑意味着速度场完全良态:没有突然的跳跃,没有无限的速度,没有数学失效的点。但近一个世纪以来,任何人能够证明的最佳存在性结果,只能保证一些更弱的东西。这些被称为弱解。

那么什么是弱解?它不是近似解。它不是"差不多对"。它是方程的精确解,但遵循的是放松的规则。正常的("经典")解要求速度足够光滑,以便你能在每一个点计算其变化率。弱解跳过了这个要求。它不是逐点检验方程,而是在空间区域上"平均地"检验它们。

这里有个类比。经典解是一个学生,逐步展示所有计算过程来解决每道考题。弱解是一个学生,无法向你展示中间步骤,但其最终答案对于你可能提出的每个问题都是可证正确的。你看不到他们的解题过程,但答案总是能验证通过。

为什么要接受这种做法?因为有时方程太过剧烈,经典解无法应对。流体可能发展出速度变化极其剧烈的区域,以至于你根本无法在那里计算变化率。数学崩溃了。弱解让你能够在经典解放弃的地方继续前进。它们是在情况变糟时维持方程存活的安全网。

代价是:弱解可能不唯一。你可能从完全相同的流动出发得到多个弱解,而没人能告诉你哪一个是"真正的答案"。这是个问题,因为物理学说流体应该做一件特定的事,而非多件。而且弱解可能不光滑。光滑性正是千禧年大奖问题所要求的,也是没人能够证明的。

Navier-Stokes方程的弱解将逐点偏微分方程替换为分布意义下的形式。考虑$\mathbb{R}^3 \times (0,T)$上的不可压缩系统:

$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$

一个无散度向量场$u \in L^2_{\mathrm{loc}}(\mathbb{R}^3 \times [0,T))$是弱解,如果对每个光滑、紧支撑、无散度的测试函数$\varphi$:

$\int_0^T \!\int_{\mathbb{R}^3} \bigl[-u \cdot \partial_t \varphi - (u \otimes u) : \nabla \varphi - \nu \, u \cdot \Delta \varphi \bigr] \, dx \, dt = \int_{\mathbb{R}^3} u_0 \cdot \varphi(x,0) \, dx.$

压力项从此形式中消失,因为$\varphi$是无散度的。所有导数都通过分部积分从$u$转移到了$\varphi$上。要点在于:$u$不需要在经典意义下可微。它只需要足够可积,使得这些积分收敛。

弱形式不是近似。逐点满足偏微分方程的经典解也对每个测试函数满足弱形式(反向分部积分)。反之则不成立:弱解可能不够光滑,无法逐点满足偏微分方程。

这很重要,因为:(1) 对于$L^2$初始数据,弱解全局存在(Leray 1934),而经典解在三维情况下仅已知局部时间存在;(2) 弱解一般不已知是唯一的;(3) 弱解是否总是光滑的问题,对于适当的初始数据,等价于克雷千禧年问题。

1934年,Jean Leray做了一件至今仍定义该领域的事情。在一篇73页的论文中,他证明了三维Navier-Stokes方程的弱解对于任何合理的初始流动都全时存在。任何。只要起始速度不是无限能量或物理上荒谬的,Leray保证你会得到一个永久持续的解。这是第一次有人证明三维方程的全局存在性结果,九十多年后,它仍然是我们拥有的最强无条件存在性定理。

他的策略很巧妙。由于流体速度如何反馈到自身(这是非线性),实际方程太难直接求解。所以Leray稍微模糊了方程,就像给图像添加微小的高斯滤波器。模糊后的方程足够温和,可以求解。然后他将模糊度逐渐降至零,并证明解不会飞散。它们稳定到某个在弱意义下满足原始、未模糊方程的东西。

但这是Leray没有证明的。唯一性。他的方法产生至少一个弱解,但可能有其他从相同流动出发的解。他无法排除这种可能。他也没有证明光滑性。他的解具有有限能量并满足能量不等式:摩擦可以耗散能量,但能量不能凭空出现。就是这样。没有更多。

Leray本人怀疑可能形成奇点。他勾勒了一个可能的样子:流体向一个点坍缩,越来越快,将所有能量集中到越来越小的区域,就像一个旋涡以无限速度收缩到一个点。1996年,Nečas、Růžička和Šverák证明了这种确切的自相似坍缩不可能发生。Leray关于潜在爆破形状的猜测是错误的。爆破是否以任何形式发生?没人知道。

1951年,Eberhard Hopf将Leray的构造扩展到有界容器中的流体(不仅仅是整个无限空间),由此产生的类别被称为Leray-Hopf弱解:满足能量不等式的弱解。这是标准概念。当研究人员说"弱解"而不作进一步限定时,他们几乎总是指这个。

还有一件事。即使在Leray-Hopf弱解中,也有一个更挑剔的子类叫做适当弱解。这些不仅在全局上满足能量不等式(总能量不增长)。它们在局部上也满足:能量不能在流体的一个角落里秘密堆积而从另一个角落流失。Caffarelli、Kohn和Nirenberg(CKN)在1982年专门针对这个更小的类证明了他们著名的部分正则性结果。不要混淆这两者:CKN适用于适当弱解,而不是所有Leray-Hopf解。

Leray的1934年论文确立了以下内容:对于任何无散度的$u_0 \in L^2(\mathbb{R}^3)$,存在至少一个弱解$u$到$\mathbb{R}^3 \times (0,\infty)$上的Navier-Stokes方程,满足:

• $u \in L^\infty(0,\infty; L^2(\mathbb{R}^3)) \cap L^2(0,\infty; \dot{H}^1(\mathbb{R}^3))$
• 能量不等式:对于几乎处处的$t > 0$,有$\|u(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u(s)\|_{L^2}^2 \, ds \leq \|u_0\|_{L^2}^2$

构造通过磨光进行。将非线性项$(u \cdot \nabla)u$替换为$(u_\varepsilon \cdot \nabla)u$,其中$u_\varepsilon = J_\varepsilon * u$是空间磨光。正则化系统有全局光滑解(磨光消除了最坏的非线性相互作用)。Leray获得了正则化解的一致能量界,然后提取了弱收敛子序列。极限满足弱形式和能量不等式。

Leray没有建立唯一性。紧性论证给出了至少一个聚点的存在性;不同的子序列可能收敛到不同的极限。三维Leray-Hopf弱解的唯一性至今仍是开放问题。

Hopf(1951)将构造适配到带Dirichlet边界条件的有界区域$\Omega \subset \mathbb{R}^3$,使用Galerkin近似(投影到有限维子空间)而非磨光。由此产生的类别,即满足能量不等式的弱解,以两人的名字命名:Leray-Hopf弱解。

Caffarelli-Kohn-Nirenberg定理(1982)涉及一个更严格的类:适当弱解,它额外满足形式为

$\partial_t \left(\frac{|u|^2}{2}\right) + \nabla \cdot \left(\left(\frac{|u|^2}{2} + p\right)u\right) + \nu |\nabla u|^2 \leq \nu \Delta \left(\frac{|u|^2}{2}\right)$

的局部能量不等式,在分布意义下成立。CKN证明了对于任何适当弱解,时空奇异集的一维抛物Hausdorff测度为零。这意味着奇点,如果存在,是极其稀疏的(它们的一维抛物Hausdorff测度为零)。但该定理没有说明奇点是否实际发生,并且它仅适用于适当弱解,而不是所有Leray-Hopf解。

Leray本人考虑了形式为$u(x,t) = (T-t)^{-1/2} U(x / (T-t)^{1/2})$的自相似爆破的可能性。Nečas、Růžička和Šverák(1996)证明了对于$L^3(\mathbb{R}^3)$中的解,不存在这样的自相似爆破,Tsai(1998)在相应假设下排除了某些渐近自相似爆破情景。潜在奇点的形状,如果有的话,仍然未知。

弱解全局存在。但它们可能不唯一,也可能不光滑。我们能做得更好吗?

可以,但只是暂时的。强解是升级版:它们有足够的正则性,使方程几乎处处成立,而不仅仅是"平均地"成立。再往上一层,光滑解或经典解才是方程逐点成立的情形。对于三维的光滑初始数据,强解在短时间内存在。多短?这取决于起始流动有多剧烈。平静、温和的流动获得更长的保证。剧烈、湍流的起始条件?微秒级。

而且没人能证明这些强解最终不会爆破。

1962年,James Serrin证明了类似于晋升规则的东西。是这样的:如果一个弱解恰好保持足够良态(不太大,不将能量集中到越来越小的区域),那么它实际上一直都是光滑的。你可以晋升它。而且根据弱强唯一性原理,它也是具有这些起始条件的唯一弱解。一个解,光滑且唯一,案件结束。但如果你无法验证解保持温和?什么也没有。你被困住了。

这是一个条件性结果。如果解不太剧烈,那么它是完全良态的。全部困难在于证明"如果"。

在二维中,能量估计足够强,以至于每个弱解自动通过Serrin的测试。完成。这就是为什么二维已解决。在三维中,估计恰好略低于所需,而弥合这个差距就是整个博弈。

研究人员也找到了其他条件性测试,每一个都是不同的攻击角度:"证明关于解的这一件特定事情,我就免费给你光滑性。"无条件地证明其中任何一个都将解决千禧年问题。没人成功过。关于人们尝试过的不同证明策略的综述,有专门的页面。

Navier-Stokes方程的强解是具有足够正则性的解,使得偏微分方程几乎处处逐点成立,并且非线性项$(u \cdot \nabla)u$作为函数而非仅作为分布是良定义的。通常这意味着$u \in L^\infty(0,T; H^1(\mathbb{R}^3)) \cap L^2(0,T; H^2(\mathbb{R}^3))$。一个相关但不同的框架是Fujita-Kato(1964),它在临界空间中构造局部温和解。

对于足够正则的Sobolev数据,强解的局部存在性已充分确立。在临界空间中,Fujita-Kato(1964)框架通过不动点论证构造局部温和解:

$u(t) = e^{\nu t \Delta} u_0 - \int_0^t e^{\nu(t-s)\Delta} \mathbb{P} \nabla \cdot (u \otimes u)(s) \, ds,$

其中$\mathbb{P}$是到无散度场的Leray投影。该积分方程通过适当函数空间中的压缩映射原理具有唯一局部解。对于临界空间中的小数据($L^3$,$\dot{H}^{1/2}$,$BMO^{-1}$),解是全局的。

问题在于大数据强解是否永久持续。关键的条件性结果是Serrin(1962)的:如果Leray-Hopf弱解满足$u \in L^p(0,T; L^q(\mathbb{R}^3))$,其中

$\frac{2}{p} + \frac{3}{q} \leq 1, \qquad q > 3,$

那么$u$在$(0,T] \times \mathbb{R}^3$上光滑,并且是具有给定初始数据的唯一Leray-Hopf弱解。这些被称为Prodi-Serrin条件(Prodi 1959建立了相关结果)。

端点$q = 3$($p = \infty$)由Escauriaza、Seregin和Šverák(2003)解决:如果$u \in L^\infty(0,T; L^3(\mathbb{R}^3))$,则$u$在时间$T$不爆破。这是Prodi-Serrin尺度的端点,也是已知最锐利的延拓准则之一。

在二维中,能量界$u \in L^\infty_t L^2_x \cap L^2_t H^1_x$结合Ladyzhenskaya不等式$\|f\|_{L^4}^2 \leq C\|f\|_{L^2}\|\nabla f\|_{L^2}$(二维特有)给出$u \in L^4_t L^4_x$,满足Serrin条件$2/4 + 2/4 = 1$(在二维版本中带$2/p + 2/q \leq 1$)。在三维中,能量界通过Sobolev嵌入给出$u \in L^{10/3}_t L^{10/3}_x$,满足$2/(10/3) + 3/(10/3) = 3/2 > 1$。Serrin条件恰好以对应于超临界差距的幅度失效。关于这一结构性障碍的更多信息,见为什么Navier-Stokes困难。

光滑解是黄金标准。速度场在任何地方、任何时间都是完全良态的。没有突然的跳跃。没有无限的速度。无论放大多少,解都保持良好。

克雷千禧年大奖问题,由Charles Fefferman在2000年提出,提出了一个能写在索引卡上的问题。从任何填充三维空间的光滑、物理上合理的速度场开始。Navier-Stokes方程是否总是产生一个永久持续的光滑解,还是你能找到一个起始流动,其解最终爆破?

任一答案都值一百万美元。

这是我们的现状。没人证明过光滑解在三维中总是全局存在,也没人构造出爆破。自Leray 1934年论文以来,我们一直困在中间,超过九十年来数学中最困难的开放问题之一,我们仍然不知道答案落在哪一边。

短期?没问题。对于光滑起始数据,方程确实在一段时间内产生光滑解。流体开始运动,数学有效,一切都很干净。但后来会发生什么?解是否永远保持光滑,还是会达到速度冲向无穷的点?

如果它确实保持光滑,就会发生好事。那个光滑解也自动满足弱解的放松规则,所以它是弱解。而根据弱强唯一性,具有这些起始条件的其他弱解不可能存在。所以如果有人证明了全局光滑性,整个层次结构就会坍缩:弱、强和光滑都将被证明是同一件事,一个永久完全良态的单一唯一解。这就是使这个问题如此吸引人和如此困难的原因。我们能证明存在的(弱解)和我们想要的(光滑解)之间的差距,正是这百万美元问题的内容。

克雷千禧年问题(Fefferman 2000)问:对于任何满足对所有$\alpha, K$有$|\partial^\alpha u_0(x)| \leq C_{\alpha K}(1 + |x|)^{-K}$的无散度$u_0 \in C^\infty(\mathbb{R}^3)$,是否存在$u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$和$p \in C^\infty(\mathbb{R}^3 \times [0,\infty))$满足Navier-Stokes方程,对所有$t \geq 0$有$\int_{\mathbb{R}^3} |u(x,t)|^2 \, dx < C$?

备选方案(也值得奖金):找到上述类中的$u_0$,使得不存在这样的光滑解。

已知内容:

• 局部存在性:对于$u_0 \in H^s(\mathbb{R}^3)$,$s \geq 1/2$,唯一的局部温和解/强解存在于$[0,T^*)$上,对于某个依赖于$\|u_0\|_{H^s}$的$T^* > 0$,并且对于每个正时间$t > 0$它是光滑的。如果$T^* < \infty$,则当$t \to T^*$时$\|u(t)\|_{H^s} \to \infty$。
• 小数据全局存在性:如果$\|u_0\|_{\dot{H}^{1/2}}$(或$\|u_0\|_{L^3}$,或$\|u_0\|_{BMO^{-1}}$)小于一个通用常数,解是全局和光滑的。关键参考文献:$\dot{H}^{1/2}$(Fujita-Kato 1964),$BMO^{-1}$(Koch-Tataru 2001),在$L^3$中有类似结果(Kato 1984)。
• 弱强唯一性:如果强解存在于$[0,T]$上,则每个具有相同初始数据的Leray-Hopf弱解在$[0,T]$上与之重合。这由Serrin(1962)建立,并由后续工作改进。这意味着证明正则性也解决了Leray-Hopf类中的唯一性。

层次结构向上坍缩:光滑$\Rightarrow$强$\Rightarrow$弱,而弱强唯一性意味着光滑解,如果存在,就是唯一的Leray-Hopf弱解。因此千禧年问题等价于问:所有具有光滑初始数据的Leray-Hopf弱解本身是否光滑?弱解存在"(Leray 1934)和"光滑解存在"(开放)之间的差距正是大奖问题。

如果弱解存在并描述流体,为什么还要关心光滑性?

三个原因。

第一,唯一性。物理学要求一个答案。给我流体的初始状态,我应该能准确告诉你它接下来会做什么。不是"这里有几种可能性,随便选一个"。但弱解不保证这一点。多个弱解可能从相同的起始流动出现,无法判断哪一个是真实流体遵循的。方程会变成菜单而非配方。那不是物理学。

第二,数值可靠性。许多重要的流体模拟基于Navier-Stokes或密切相关的模型:天气预报、空气动力学、通过动脉的血流等。使网格更精细,模拟应该收敛到真实答案。没有光滑性和唯一性保证?没有定理说这在每个三维场景中实际发生。模拟有效。我们无法完全解释为什么。

第三,极端物理。如果奇点可以形成,那是自然在向我们发送消息。Navier-Stokes方程,我们最好的流体运动模型,将有一个内置的失效日期:在某个极端尺度上,模型本身停止工作,方程在告诉我们,"你需要新物理学。"

这不是技术细节。这是贯穿一切的断层线。一边是存在性(Leray,1934,完成)。另一边是光滑性(开放,一百万美元)。为什么跨越如此困难?三维中的能量估计恰好略低于所需,九十年的努力,数百篇论文,耗尽整个职业生涯的尝试,没人能弥合这个差距。

目前正在进行的每个证明策略都是试图弥合这一鸿沟的尝试。证明弱解是光滑的。或者证明它们不是。两个字还是三个字。

三维不可压缩Navier-Stokes方程的解层次结构为:

$\text{smooth} \subset \text{strong} \subset \text{Leray-Hopf weak} \subset \text{distributional weak}$

每个层次的已知内容:

千禧年问题位于全局Leray-Hopf弱存在性和全局光滑正则性之间的差距。核心困难:能量不等式提供$u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$,这比临界缩放$\dot{H}^{1/2}$低半阶导数。弥合这一超临界差距等价于证明全局正则性。

弱强唯一性意味着如果光滑解全局存在,层次结构就会坍缩。更精确地:如果对于给定的$u_0 \in C^\infty$,光滑解$u$存在于$[0,T]$上,则每个具有相同数据的Leray-Hopf弱解在$[0,T]$上等于$u$。因此全局光滑性$\Rightarrow$Leray-Hopf类中的全局唯一性。

反过来,对于光滑可容许初始数据的Leray-Hopf弱解的非唯一性将意味着全局光滑解不能对该数据类持续(因为光滑解会强制唯一性)。关于凸积分的最近工作(建立在De Lellis-Székelyhidi对Euler的工作之上,由Buckmaster-Vicol 2019扩展到在Leray-Hopf正则性以下构造Navier-Stokes的非唯一弱解)表明分布弱解可能高度非唯一。这种非唯一性是否延伸到Leray-Hopf类是一个重大开放问题,对千禧年问题有直接影响。

关于攻击这一差距的证明策略的当前状态,以及它如此顽固的结构性原因,见为什么Navier-Stokes困难。