为什么Navier-Stokes问题困难：非线性、压力与湍流

人们在物理学中最初遇到的许多方程都是线性的:输入加倍,响应也加倍。纳维-斯托克斯方程并非如此。

纳维-斯托克斯?非线性。流体的速度影响其自身的变化率,这意味着流体在推动自己。想象一下,试图预测人群将去往何处,当每个人的移动都取决于周围所有人在做什么,而那些人在做什么又取决于他们周围的所有人,如此螺旋式向外无限延伸。这就是你正面对的情况。

罪魁祸首是自相互作用项$(u \cdot \nabla)u$。它创造了反馈循环,小扰动被放大为大扰动,这就是为什么流体湍流如此复杂混乱(更多内容见子问题)。

对流非线性项$(u \cdot \nabla)u$ 是根本障碍。在涡量形式$\omega = \nabla \times u$中,方程变为

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$

涡拉伸项$(\omega \cdot \nabla)u$ 没有符号。它可以无界地放大涡量。在二维情况下该项消失(因为$\omega$ 是垂直于流动的标量),这就是为什么二维全局正则性已知 (Ladyzhenskaya, 1969)。在三维情况下,涡拉伸是有限时间爆破的主要候选机制。

问题的关键在于:非线性项关于二次$u$$H^1$ 能量估计给出$\|\nabla u\|_{L^2}$,但控制$(u \cdot \nabla)u$$L^2$ 通常需要更强的信息,例如$u \in L^\infty$ 加上$\nabla u \in L^2$,或等价的临界尺度控制;仅能量类本身无法提供这一点。

纳维-斯托克斯方程具有尺度对称性。放大一个解,按正确的比例使所有东西变小变快,你会得到另一个完全有效的解。这种对称性在数学上很自然,但在分析上却很危险。

为什么?我们唯一能可靠控制的量是流体的总能量,而它处于完全错误的尺度上,告诉我们大局情况,但对爆破实际形成的微观尺度上发生的事情却毫无说明。

想象一下,监控一个城市的总用电量来检测一个火花。有用吗?当然。足够精细吗?差得远呢。这个缺口就是整个问题所在。

在自然尺度变换$u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$下,临界Sobolev空间是$\dot{H}^{1/2}(\mathbb{R}^3)$ (等价于$L^3$)。如果一个受控范数在放大到更小尺度时反而变小,那么相对于这类控制,方程就是超临界的:这个界低于看见小尺度集中的临界尺度。

能量不等式给出$u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$的控制。在能量层面,三维问题是超临界的:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

当$\lambda \to \infty$(即放大到更小尺度)时,这些界相对于不变的量$\|u\|_{L^3}$和$\|u\|_{\dot{H}^{1/2}}$会变弱。非线性项和粘性项在标度下是同步变化的;问题在于已知的先验估计落在临界尺度之下。从能量类走到临界范数,这正是核心困难。

观察一条河流。流体运动变得混乱。湍流。大涡旋破碎成小涡旋,小涡旋又破碎成更小的涡旋,一路级联到微观尺度,最终粘性将其平滑。

Kolmogorov在1941年描述了这种能量级联,纳维-斯托克斯方程完美地捕捉了它。但这里有一个让人夜不能寐的问题:如果能量集中到越来越小的区域的速度快于粘性耗散它的速度,会怎样?那就是爆破。

这真的会发生吗?还是粘性总是会获胜?这就是悬而未决的问题,句号。关于从雷诺数到这一小尺度图景的物理桥梁,见雷诺数、湍流与小尺度的重要性。

Kolmogorov的K41理论预测惯性区间$E(k) \sim \varepsilon^{2/3} k^{-5/3}$中的能量谱$k_f \ll k \ll k_\eta$,其中$k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ 是Kolmogorov耗散波数。跨尺度的恒定能量通量。这是清晰的图景。

正则性问题询问更深层的东西:这种级联能否退化?耗散尺度$k_\eta^{-1}$ 能否在有限时间内收缩到零,使得$\|\nabla u\|_{L^2} \to \infty$ 而总能量保持有限?

Onsager的耗散反常 猜想(1949)在某种意义上说是的:在消失粘性极限$\nu \to 0$中,能量耗散持续存在,而Euler方程的弱解 可以耗散能量。对于Hölder指数低于$1/3$ 的情况已得到确认(Isett, 2018; Buckmaster et al., 2018)。但这对Navier-Stokes正则性意味着什么仍然真正不清楚。关于区域层面的直观理解,见Reynolds数、湍流以及为什么小尺度重要

纳维-斯托克斯方程中的压力很奇怪。它根本不是一个独立变量;速度通过一个约束完全决定它:流体是不可压缩的,所以它不能被挤压。

这使压力具有非局部性。在不可压缩模型中,改变一个区域的速度会全局地影响压力场,因为压力是由该时刻的整个速度场决定的。

对于分析来说,这是毁灭性的。你不能研究单个点发生的事情而不考虑整个流体。局部推理?忘掉它吧。

不可压缩性约束$\nabla \cdot u = 0$ 通过Poisson方程确定压力

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j)$$

因此$p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$,涉及Riesz变换(奇异积分算子)。压力是速度的非局部函数,这种非局部性是逐点或空间局部估计的核心障碍。

标准最大值原理论证在这里失效:尽管粘性项$\nu \Delta u$ 是耗散的,但压力梯度$-\nabla p$ 可以从远处区域集中能量。Caffarelli-Kohn-Nirenberg理论通过抛物柱面上的局部能量不等式处理这个问题,但从中提取逐点正则性仍然是困难的步骤。

对于标准设置下的二维不可压缩纳维-斯托克斯方程,全局光滑解是已知的;这在包括Ladyzhenskaya(1969)的经典工作中已经建立。见为什么二维更简单。

三维?一切都崩溃了,原因归结为一个机制:涡拉伸。在二维中,涡旋可以旋转和合并,但它们不能拉伸。在三维中,流体可以抓住涡管并将它们拉得越来越细,潜在地将每一点能量集中到一根无限细的丝中。

这种拉伸能否在有限时间内失控到无穷大,还是粘性总是会介入?这就是百万美元问题。确实如此。

二分法是尖锐的。

二维: 涡量$\omega$ 是满足$\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$ 的标量。最大值原理给出$\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$,BKM蕴含全局正则性,而涡拉伸项$(\omega \cdot \nabla)u$ 恒为零。

三维: 涡量$\omega \in \mathbb{R}^3$ 满足$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$,其中拉伸项$(\omega \cdot \nabla)u$ 可以超线性地放大$|\omega|$ (通过Biot-Savart形式上有$\sim |\omega|^2$)。没有可用的最大值原理。

涡能$\|\omega\|_{L^2}^2$ 满足

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

控制$\|\nabla u\|_{L^\infty}$ 需要$\omega \in L^\infty$,这需要控制$\|\nabla u\|_{L^\infty}$。循环论证。现有技术都未能打破这个循环。

本文是问题的一部分。

这些障碍促使数学家们将问题分解为子问题,并为每个问题开发专门方法。

关于粘性项为何有帮助但还不够的背景,见欧拉方程与纳维-斯托克斯方程。

本文是问题的一部分。

正则性问题可分解为可处理的部分;见子问题。为攻克这些障碍而构建的分析工具在方法中涵盖。粘性项为何有帮助但仍未能弥合缺口?欧拉方程与纳维-斯托克斯方程。