纳维－斯托克斯方程

纳维－斯托克斯方程（Navier-Stokes equations），以克劳德-路易·纳维(Claude-Louis Navier)和乔治·加布里埃尔·斯托克斯命名，是一组描述像液体和空气这样的流体物质的方程。这些方程建立了流体的粒子动量的改变率(加速度)和作用在液体内部的压力的变化和耗散粘滞力(类似于摩擦力)以及重力之间的关系。这些粘滞力产生于分子的相互作用，能告诉我们液体有多粘。这样，纳维-斯托克斯方程描述作用于液体任意给定区域的力的动态平衡。

他们是最有用的一组方程之一，因为它们描述了大量对学术和经济有用的现象的物理过程。它们可以用于模拟天气，洋流，管道中的水流，星系中恒星的运动，翼型周围的气流。它们也可以用于飞行器和车辆的设计，血液循环的研究，电站的设计，污染效应的分析，等等。

纳维-斯托克斯方程依赖微分方程来描述流体的运动。这些方程，和代数方程不同，不寻求建立所研究的变量（譬如速度和压力）的关系，而是建立这些量的变化率或通量之间的关系。用数学术语来讲，这些变化率对应于变量的导数。这样，最简单情况的0粘滞度的理想流体的纳维-斯托克斯方程表明加速度（速度的导数，或者说变化率）是和内部压力的导数成正比的。

这表示对于给定的物理问题的纳维-斯托克斯方程的解必须用微积分的帮助才能取得。实用上，只有最简单的情况才能用这种方法解答，而它们的确切答案是已知的。这些情况通常涉及稳定态（流场不随时间变化）的非湍流，其中流体的粘滞系数很大或者其速度很小（小的雷诺数）。

对于更复杂的情形，例如厄尔尼诺这样的全球性气象系统或机翼的升力，纳维-斯托克斯方程的解必须借助计算机。这本身是一个科学领域，称为计算流体力学。

虽然湍流是日常经验中就可以遇到的，但这类问题极难求解。一个$1,000,000的大奖由克雷数学学院于2000年5月设立，奖给对于能够帮助理解这一现象的数学理论作出实质性进展的任何人。

基本假设

在解释纳维-斯托克斯方程的细节之前，我们必须首先对流体的性质作几个假设。第一个假设是流体是连续的。这强调它不包含形成内部的空隙，例如，溶解的气体的气泡，而且它不包含雾状粒子的聚合。另一个必要的假设是所有涉及到的场，全部是可微的，例如压强，速度，密度，温度，等等。

该方程从质量，动量，和能量的守恒的基本原理导出。对此，有时必须考虑一个有限的任意体积，称为控制体积，在其上这些原理很容易应用。该有限体积记为$\Omega$，而其表面记为$\partial\Omega$。该控制体积可以在空间中固定，也可能随着流体运动。这会导致一些特殊的结果，我们将在下节看到。

随体导数

运动流体的属性的变化，譬如大气中的风速的变化，可以有两种不同的方法来测量。可以用气象站或者气象气球上的风速仪来测量。显然，第一种情况下风速仪测量的速度是所有运动的粒子经过一个固定点的速度，而第二种情况下，仪器在测量它随着流体运动时速度的变化。同样的论证对于密度、温度、等等的测量也是成立的。因此，当作微分时必须区分两种情况。第一种情况称为空间导数或者欧拉导数。第二种情况称为实质或拉格朗日导数。例子请参看随体导数条目。

随体导数定义为算子（operator）：

$\frac{\mathrm{D}}{\mathrm{D}t}(\star) = \frac{\partial(\star)}{\partial t} + (\mathbf{v}\cdot\nabla)(\star)$

其中$\mathbf{v}$是流体的速度。方程右边的第一项是普通的欧拉导数（也就是在静止参照系中的导数）而第二项表示由于流体的运动带来的变化。这个效应称为移流（advection）。

L的守恒定律在一个控制体积上的积分形式是：

$\frac{\mathrm{d}}{\mathrm{d}t}\int_\Omega \mathbf{L}\,\mathrm{d}\Omega = 0$

因为Ω是共动的，它随着时间而改变，所以我们不能将时间导数和积分简单的交换。

$\frac{\mathrm{d}}{\mathrm{d}t}\int_\Omega \mathbf{L}\,\mathrm{d}\Omega=\int_{\Omega}\frac{\partial}{\partial t}\mathbf{L}\,\mathrm{d}\Omega+\int_{\partial \Omega} \mathbf{L} \left(\mathbf{v}\cdot \mathbf{n}\,\mathrm{d}\partial\Omega\right)=\int_\Omega \left[\frac{\partial}{\partial t}\mathbf{L}+\nabla\cdot\left(\mathbf{L}\mathbf{v}\right)\right]\mathrm{d}\Omega=0$

因为这个表达式对于所有$\Omega$成立，它可以简化为：

$\frac{\mathrm{D}}{\mathrm{D}t}\mathbf{L} + \left(\nabla\cdot \mathbf{v}\right) \mathbf{L} = \frac{\partial}{\partial t} \mathbf{L}+ \nabla\cdot\left(\mathbf{v} \mathbf{L}\right) = 0$

对于不是密度的量（因而它不必在空间中积分），$\frac{\mathrm{D}}{\mathrm{D}t}$给出了正确的共动时间导数。

守恒定律

NS方程可以从守恒定律通过上述变换导出，并且需要用状态定律来闭合。

在控制体积上，使用上述变换，下列的量视为守恒：

• 质量
• 能量
• 动量
• 角动量

连续性方程

质量的守恒写作：

$\frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf{v})= 0$

其中

$\rho$ 是流体的密度。

在不可压缩流体的情况 $\rho$不是时间或空间的函数。方程简化为：

$\nabla\cdot\mathbf{v} = 0$

动量守恒

动量守恒写作：

$\frac{\partial}{\partial t}\left(\rho\mathbf{v}\right) + \nabla(\rho\mathbf{v}\otimes\mathbf{v}) = \sum\rho\mathbf{f}$

注意$\mathbf{v}\otimes\mathbf{v}$是一个张量，$\otimes$代表张量积。

我们可以进一步简化，利用连续性方程，这成为：

$\rho\frac{D\mathbf{v}}{D t}=\sum\rho\mathbf{f}$

我们可以认出这就是通常的F=ma。

方程组

一般形式

方程组的形式

纳维-斯托克斯方程的一般形式是：

$\rho\frac{\mathrm{D}\mathbf{v}}{\mathrm{D} t} = \nabla \cdot\mathbb{P} + \rho\mathbf{f}$

关于动量守恒。张量$\mathbb{P}$代表施加在一个流体粒子上的表面力（应力张量）。除非流体是由象旋涡这样的旋转自由度组成，$\mathbb{P}$是一个对称张量。一般来讲，我们有如下形式：

$\mathbb{P} = \begin{pmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{pmatrix} = - \begin{pmatrix} p&0&0\\ 0&p&0\\ 0&0&p \end{pmatrix} + \begin{pmatrix} \sigma_{xx}+p & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy}+p & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz}+p \end{pmatrix}$

其中$\sigma$是法向约束，而$\tau$是切向约束。

迹 $\sigma_{xx}+\sigma_{yy}+\sigma_{zz}$在流体处于平衡态时为0。这等价于流体粒子上的法向力的积分为0。

我们再加上连续性方程：

$\frac{\mathrm{D}\rho}{\mathrm{D}t} + \rho\nabla\cdot\mathbf{v} = 0$

对于处于平衡的液体，$\mathbb{P}$的迹是3p。

其中

p是压强

最后，我们得到：

$\rho\frac{\mathrm{D}\mathbf{v}}{\mathrm{D} t} = -\nabla p + \nabla \cdot\mathbb{T} + \rho\mathbf{f}$

其中$\mathbb{T}$是$\mathbb{P}$的非对角线部分。

闭合问题

这些方程是不完整的。要对它们进行完备化，必须对$\mathbb{P}$的形式作一些假设。例如在理想流体的情况$\tau$分量为0。用于完备方程组的方程是状态方程。

再如，压强可以主要是密度和温度的函数。

要求解的变量是速度的各个分量，流体密度，静压力，和温度。流场假定为可微并连续，使得这些平衡得以用偏微分方程表达。这些方程可以转化为涡度和流函数这些次变量的威尔金森方程组。解依赖于流体的性质（例如粘滞度、比热、和热导率），并且依赖于所研究的区域的边界条件。

$\mathbb{P}$的分量是流体的一个无穷小元上面的约束。它们代表垂直和剪切约束。$\mathbb{P}$是对称的，除非存在非零的自旋密度。

所谓非牛顿流体是就是其中该张量没有特殊性质使得方程的特殊解出现的流体

特殊形式

这些是问题的特定的常见简化，有时解是已知的。

牛顿流体

在牛顿流体中，如下假设成立:

$\tau_{ij}=\mu\left(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}\delta_{ij}\nabla\cdot\mathbf{v}\right)$

其中

$\mu$是液体的粘滞度。

$\rho \left(\frac{\partial \mathbf{v}}{\partial t}+\nabla_{\mathbf{v}}\mathbf{v}\right)=\rho \mathbf{f}-\nabla p +\mu\left(\Delta \mathbf{v}+\frac{1}{3}\nabla\left(\nabla\cdot \mathbf{v}\right)\right)$

$\rho \left(\frac{\partial v_i}{\partial t}+v_j\frac{\partial v_i}{\partial x_j}\right)=\rho f_i-\frac{\partial p}{\partial x_i}+\mu\left(\frac{\partial ^2 v_i}{\partial x_j \partial x_j}+\frac{1}{3}\frac{\partial ^2 v_j}{\partial x_i \partial x_j}\right)$

其中为简化书写，对脚标使用了爱因斯坦求和约定。

不采用简化书写的完整形式非常繁琐，分别为：

动量守恒：

$\rho \cdot \left({\partial u \over \partial t}+ u {\partial u \over \partial x}+ v {\partial u \over \partial y}+ w {\partial u \over \partial z}\right) = k_x \partial p \over \partial x} + {\partial \over \partial x} \left[ \mu \cdot \left(2 \cdot {\partial u \over \partial x\frac{2}{3}\cdot (\nabla \cdot \vec{v} \right) \right] + {\partial \over \partial y}\left[\mu \cdot \left({\partial u \over \partial y} + {\partial v \over \partial x} \right) \right] + {\partial \over \partial z}\left[\mu \cdot \left({\partial w \over \partial x} + {\partial u \over \partial z} \right) \right]$

$\rho \cdot \left({\partial v \over \partial t}+ u {\partial v \over \partial x}+ v {\partial v \over \partial y}+ w {\partial v \over \partial z}\right) = k_y \partial p \over \partial y} + {\partial \over \partial y} \left[ \mu \cdot \left(2 \cdot {\partial v \over \partial y\frac{2}{3}\cdot (\nabla \cdot \vec{v} \right) \right] + {\partial \over \partial z}\left[\mu \cdot \left({\partial v \over \partial z} + {\partial w \over \partial y} \right) \right] + {\partial \over \partial x}\left[\mu \cdot \left({\partial u \over \partial y} + {\partial v \over \partial x} \right) \right]$

$\rho \cdot \left({\partial w \over \partial t}+ u {\partial w \over \partial x}+ v {\partial w \over \partial y}+ w {\partial w \over \partial z}\right) = k_z \partial p \over \partial z} + {\partial \over \partial z} \left[ \mu \cdot \left(2 \cdot {\partial w \over \partial z\frac{2}{3}\cdot (\nabla \cdot \vec{v} \right) \right] + {\partial \over \partial x}\left[\mu \cdot \left({\partial w \over \partial x} + {\partial u \over \partial z} \right) \right] + {\partial \over \partial y}\left[\mu \cdot \left({\partial v \over \partial z} + {\partial w \over \partial y} \right) \right]$

质量守恒：

${\partial \rho \over \partial t} + {\partial (\rho \cdot u) \over \partial x}+{\partial (\rho \cdot v) \over \partial y}+{\partial (\rho \cdot w) \over \partial z}=0$

因为密度是一个未知数，我们需要另一个方程。

能量守恒：

$\rho \left({\partial e \over \partial t}+ u {\partial e \over \partial x}+ v {\partial e \over \partial y}+ w {\partial e \over \partial z}\right) = \left( {\partial \over \partial x} \left(\lambda \cdot {\partial T \over \partial x} \right) + {\partial \over \partial y} \left(\lambda \cdot {\partial T \over \partial y} \right) + {\partial \over \partial z} \left(\lambda \cdot {\partial T \over \partial z} \right) \right) - p \cdot \left( \nabla \cdot \vec{v} \right) + \vec{k} \cdot \vec{v} + \rho \cdot \dot{q}_s + \mu \cdot \Phi$

其中：

$\Phi = 2 \cdot \left[ \left({\partial u \over \partial x} \right)^2+\left({\partial v \over \partial y}\right)^2+\left({\partial w \over \partial z}\right)^2 \right] + \left({\partial v \over \partial x}+{\partial u \over \partial y} \right)^2 + \left({\partial w \over \partial y}+{\partial v \over \partial z} \right)^2 + \left({\partial u \over \partial z}+{\partial w \over \partial x} \right)^2 -\frac{2}{3} \cdot \left({\partial u \over \partial x}+{\partial v \over \partial y}+{\partial w \over \partial z} \right)^2$

假设一个理想气体：

$e = c_p \cdot T - \frac{p}{\rho}$

上面是一个６个方程６个未知数的系统。（u， v， w， T， e 以及　$\rho$）。

宾汉（Bingham）流体

在宾汉流体中，我们有稍微不同的假设：

$\tau_{ij}=\tau_0 + \mu\frac{\partial v_i}{\partial x_j},\;\frac{\partial v_i}{\partial x_j}>0$

那些流体在开始流动之前能够承受一定的剪切。牙膏是一个例子。

幂律流体

这是一种理想化的流体，其剪切应力，$\tau$，由下式给出

$\tau = K \left( \frac {\partial u} {\partial y} \right)^n$

该形式对于模拟各种一般流体有用。

不可壓缩流體

其纳维－斯托克斯方程(Navier-Stoke equation)为

$\overbrace{\rho \Big( \underbrace{\frac{\partial \mathbf{v}}{\partial t}}_{ \begin{smallmatrix} \text{Unsteady}\\ \text{acceleration} \end{smallmatrix}} + \underbrace{\mathbf{v} \cdot \nabla \mathbf{v}}_{ \begin{smallmatrix} \text{Convective} \\ \text{acceleration} \end{smallmatrix}}\Big)}^{\text{Inertia}} = \underbrace{-\nabla p}_{ \begin{smallmatrix} \text{Pressure} \\ \text{gradient} \end{smallmatrix}} + \underbrace{\mu \nabla^2 \mathbf{v}}_{\text{Viscosity}} + \underbrace{\mathbf{f}}_{ \begin{smallmatrix} \text{Other} \\ \text{forces} \end{smallmatrix}}$

動量守恆和

$\nabla\cdot\mathbf{v}=0$

質量守恆。

其中，對不可壓縮牛頓流體來說，只有對流項(convective terms)為非線性形式。對流加速度(convective acceleration)來自於流體流動隨空間之變化所產生的速度改變，例如：當流體通過一個漸縮噴嘴(convergent nozzle)時，流體產生加速之情況。由於此項的存在，對於暫態運動中的流體來說，其流場速度變化不再單是時間的函數，亦與空間有關。

另外一個重要的觀察重點，在於黏滯力(viscosity)在流場中的以流體速度作拉普拉斯運算來表現。這暗示了在牛頓流體中，黏滯力為動量擴散(diffusion of momentum)，與熱擴散方程式非常類似。

$e_{ij}=\frac{1}{2}\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right)$;

$\Delta=e_{ii}$ 是散度，

$\delta_{ij}$是克罗内克记号。

若$\mu$在整个流体上均匀，动量方程简化为

$\rho\frac{Du_i}{Dt}=\rho f_i-\frac{\partial p}{\partial x_i} +\mu \left( \frac{\partial^2u_i}{\partial x_j\partial x_j}+ \frac{1}{3}\frac{\partial\Delta}{\partial x_i}\right)$

（若 $\mu=0$ 这个方程称为欧拉方程；那里的重点是可压缩流和冲击波）。

如果现在再有$\rho$为常数，我们得到如下系统：

$\rho \left({\partial v_x \over \partial t}+ v_x {\partial v_x \over \partial x}+ v_y {\partial v_x \over \partial y}+ v_z {\partial v_x \over \partial z}\right)= \mu \left[{\partial^2 v_x \over \partial x^2}+{\partial^2 v_x \over \partial y^2}+{\partial^2 v_x \over \partial z^2}\right]-{\partial p \over \partial x} +\rho g_x$

$\rho \left({\partial v_y \over \partial t}+ v_x {\partial v_y \over \partial x}+ v_y {\partial v_y \over \partial y}+ v_z {\partial v_y \over \partial z}\right)= \mu \left[{\partial^2 v_y \over \partial x^2}+{\partial^2 v_y \over \partial y^2}+{\partial^2 v_y \over \partial z^2}\right]-{\partial p \over \partial y} +\rho g_y$

$\rho \left({\partial v_z \over \partial t}+ v_x {\partial v_z \over \partial x}+ v_y {\partial v_z \over \partial y}+ v_z {\partial v_z \over \partial z}\right)= \mu \left[{\partial^2 v_z \over \partial x^2}+{\partial^2 v_z \over \partial y^2}+{\partial^2 v_z \over \partial z^2}\right]-{\partial p \over \partial z} +\rho g_z$

连续性方程（假设不可压缩性）：

${\partial v_x \over \partial x}+{\partial v_y \over \partial y}+{\partial v_z \over \partial z}=0$

N-S方程的简化版本。采用《不可压缩流》，Ronald Panton所著第二版

注意纳维－斯托克斯方程仅可近似描述液体流，而且在非常小的尺度或极端条件下，由离散的分子和其他物质（例如悬浮粒子和溶解的气体）的混合体组成的真实流体，会产生和纳维－斯托克斯方程所描述的连续并且齐性的液体不同的结果。依赖于问题的纳森数，统计力学可能是一个更合适的方法。但是，纳维－斯托克斯方程对于很大范围的实际问题是有效的，只要记住他们的缺陷是天生的就可以了。

参看

• 雷诺数
• 马赫数
• 雷诺平均纳维－斯托克斯方程

参考文献

• Inge L. Rhyming Dynamique des fluides, 1991 PPUR.
• Polyanin A.D., Kutepov A.M., Vyazmin A.V., Kazenin D.A., Hydrodynamics, Mass and Heat Transfer in Chemical Engineering, Taylor & Francis, London, 2002. ISBN 0-415-27237-8.

外部链接

• 克雷数学研究院纳维-斯托克斯方程大奖
  • 该问题的正式命题
• 纳维-斯托克斯方程的一个推导
• 纳维-斯托克斯方程的推导
• NASA关于纳维-斯托克斯方程的网页
• 纳维-斯托克斯方程（一些精确解），位于EqWorld:数学方程的世界