Navier-Stokes方程是什么:流体方程的含义与数学形式
从简单直觉到数学形式,清晰介绍流体运动的偏微分方程
纳维-斯托克斯方程是什么
纳维-斯托克斯方程是一个偏微分方程组,描述粘性流体(如水和空气)的运动。
它们用于描述管道中的水、飞机机翼周围的空气、动脉中的血液以及无数其他流动现象。
从高层次来看,它们表达的是:流体因压力、粘性和作用于其上的力而改变其运动。压力推动流体,粘性平滑运动中的急剧差异,而重力等外力可以驱动流动。
这些方程不仅仅是物理学口号。它们是流体动力学、工程学和计算模拟中大部分工作的语言。
纳维-斯托克斯方程是粘性牛顿流体的动量平衡方程。在不可压缩情形下,它们通过非线性偏微分方程组将速度场 $u(x,t)$ 和压力 $p(x,t)$ 耦合起来。
严格来说,纳维-斯托克斯是一个方程组而非单个方程。它们将动量守恒与本构定律(粘性应力与对称速度梯度成正比)一起建模。不可压缩性增加了局部体积守恒的约束。
对于克雷千禧年问题以及本网站的大部分内容,相关情形是三维不可压缩系统。
简化形式的纳维-斯托克斯方程
在最简单的常见形式中,方程组如下所示:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$
$$\nabla \cdot u = 0$$
其中:
- $u$ 是流体的速度
- $p$ 是压力
- $\nu$ 是运动粘度
- $f$ 是任何外力,例如重力
左侧描述速度如何随时间变化以及流体如何输运其自身的运动。右侧包含推动和平滑流动的力。
定义在 $\mathbb{R}^3$ 上的不可压缩Navier-Stokes方程组为
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f,$$
$$\nabla \cdot u = 0, \qquad u(x,0)=u_0(x).$$
各项具有标准解释:
- $\partial_t u$:局部时间变化
- $(u \cdot \nabla)u$:非线性对流,即流体输运其自身的速度
- $-\nabla p$:压力力
- $\nu \Delta u$:粘性扩散
- $f$:外部强迫
- $\nabla \cdot u = 0$:不可压缩性约束
这是本网站关于千禧年问题、困难性和证明策略页面中使用的形式。
纳维-斯托克斯方程的来源
这些方程源于一个简单的思想:将牛顿第二定律应用于流体的微小质点。该质点的质量乘以其加速度必须等于作用于其上的总力。
对于粘性流体,这些力主要来自压力和内摩擦。当你在流体的每一点仔细写出这种平衡时,就得到了纳维-斯托克斯方程。
因此这些方程并非任意的。它们是力等于质量乘以加速度的连续介质力学版本。关于完整的逐步推导,请参阅纳维-斯托克斯方程的推导。
推导从连续介质的动量守恒开始。在物质体积上写出线动量平衡,然后将该恒等式局部化以获得偏微分方程。
对于牛顿不可压缩流体,Cauchy应力张量具有如下形式
$$T = -pI + 2\mu D(u), \qquad D(u)=\frac{1}{2}(\nabla u + \nabla u^T),$$
其中 $\mu$ 是动力粘度。将此本构定律代入动量方程并除以密度,得到运动粘度为 $\nu = \mu/\rho$ 的常见不可压缩方程组。
在标准的常密度不可压缩设定中,连续性方程简化为 $\nabla \cdot u = 0$。关于从动量平衡和牛顿本构定律出发的完整推导,参见推导。
为什么纳维-斯托克斯方程很困难
困难的部分是非线性项 $(u \cdot \nabla)u$。流体不仅响应外力;它还推动自己。这种反馈正是使湍流和看似混乱的运动成为可能的原因。
在二维空间中,方程的表现要好得多。在三维中,我们仍然不知道每个光滑的初始流动是否永远保持光滑。
这就是为什么这些方程远超工程学范围而闻名:它们直接导致了纳维-斯托克斯千禧年问题。
它们的用途
纳维-斯托克斯方程每天都在科学和工程中使用。典型应用包括:
- 机翼和车辆周围的气流
- 天气和气候模型
- 海洋环流
- 工业流体输运
- 血流和其他生物输运问题
实际上,人们通过数值方法求解这些方程的近似,通常还带有额外的建模假设。这种实际成功正是剩余数学问题如此引人注目的原因之一。
应用工作通常在特定情形下使用纳维-斯托克斯或相关模型的数值近似:不可压缩流、可压缩流、湍流闭合、边界层近似和简化模型。
直接数值模拟、大涡模拟和雷诺平均闭合都可以追溯到相同的连续介质偏微分方程框架,但它们并未消除三维中的基础正则性问题。
实用有效性与不完整理论之间的这种分裂,正是使这一学科如此引人入胜的部分原因。
接下来读什么
如果你的主要问题是该问题是否已解决,请从纳维-斯托克斯问题解决了吗?开始。
如果你想了解广泛的数学意义,请继续阅读千禧年问题。
如果你想理解纳维-斯托克斯方程与无粘欧拉方程的区别以及粘性为何重要,请阅读欧拉方程与纳维-斯托克斯方程。
如果你想了解湍流和小尺度背后的物理直觉,请阅读雷诺数、湍流以及为什么小尺度重要。
如果你想了解主要障碍,请转到为什么它很困难。
本网站的自然后续步骤:
- 纳维-斯托克斯问题解决了吗? 了解现状问题以及弱解存在性与全局光滑理论的区别
- 千禧年问题 了解克雷的精确表述
- 欧拉方程与纳维-斯托克斯方程 了解粘性的作用以及无粘和粘性系统之间的关键差异
- 不可压缩与可压缩纳维-斯托克斯方程 了解常密度和变密度表述之间的差异
- 雷诺数、湍流以及为什么小尺度重要 了解多尺度流动背后的量级水平直觉
- 为什么它很困难 了解标度、超临界性和爆破情景