纳维-斯托克斯方程

来源：互联网发布：linux开机启动命令行编辑：程序博客网时间：2024/04/30 14:41

纳维-斯托克斯方程[编辑]

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

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

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

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

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

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

1 性质
- 1.1 非线性
- 1.2 湍流
- 1.3 适用性
2 基本假设
- 2.1 随体导数
- 2.2 守恒定律
  - 2.2.1 连续性方程
  - 2.2.2 动量守恒
3 方程组
- 3.1 一般形式
  - 3.1.1 方程组的形式
  - 3.1.2 闭合问题
4 特殊形式
- 4.1 牛顿流体
- 4.2 宾汉（Bingham）流体
- 4.3 幂律流体
- 4.4 不可压缩流体
5 参看
6 参考文献
7 外部链接

性质[编辑]

非线性[编辑]

纳维－斯托克斯方程在大多数实际情况下是非线性的偏微分方程。^[1]^[2]在某些情况下，一维流和斯托克斯流（或蠕动流）方程可以简化为线性方程组。非线性使得很多的问题很难或者没法解决，但是它却是动荡方程组模型的重要影响因素。非线性是由于对流加速（与点速度变化相关联的加速度），因此，任何对流无论湍流与否都会涉及非线性。

湍流[编辑]

湍流是时变的无序行为，这种行为常见于许多流体流动中。人们普遍认为，这是由于惯性流体被看作为一个整体而产生的：时间依赖性和对流加速度;因此惯性的影响很小的流体往往是层流（雷诺数量化了流所受惯性的大小）。据不精确了解，纳维－斯托克斯方程能够合理的描述湍流。^[3]纳维-斯托克斯方程关于湍流的数值解是非常难得到的，而且由于湍流的显著不同的混合长度，需要精细分辨率的稳定解的计算时间变得不可行（见直接数值模拟）。尝试使用层流来解决湍流问题的方法通常会导致时间不稳定的解，它不能够适当收敛。为了解决这个问题，时间平均方程，如雷诺平均纳维-斯托克斯方程（RANS），再辅以湍流模型，被用作实际的计算流体动力学模拟湍流时（CFD）的应用程序。例如Spalart-Allmaras湍流，K-ω湍流，k-ε湍流和SST模型。另一种方法解决数值的纳维－斯托克斯方程是大涡模拟（LES）。这种方法在计算上是比RANS方法更昂贵（在时间和计算机存储器上），但因为较大的湍流尺度被明确地解决所以会产生更好的结果。

适用性[编辑]

补充方程（例如，质量守恒定律)和良好的边界条件，使纳维－斯托克斯方程似乎可以进行流体运动的精确建模；甚至湍流（平均上）也符合实际观察结果。纳维－斯托克斯方程假定所研究的流体是连续（它是无限可分的，而不是由粒子组成），并且是相对静止的。在非常小的尺度或极端条件下，由离散分子组成的真实流体将产生与由纳维－斯托克斯方程建模的连续流体不同的结果。根据Knudsen数的问题，统计力学甚至分子动力学可能是一个更合适的方法。另一个限制是方程的复杂性。对于普遍流体家族的公式是存在的，但纳维－斯托克斯方程对不常见流系的应用，往往会导致非常复杂的构成，而且该配方是目前研究的一个领域。出于这个原因，这些方程通常用于描述牛顿流体。研究这种液体是“简单”的，因为粘度模型最终被线性化;真正描述其他类型流体（如血液）的普遍的模型截至2012年还不存在。

基本假设[编辑]

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

该方程从质量，动量，和能量的守恒的基本原理导出。对此，有时必须考虑一个有限的任意体积，称为控制体积，在其上这些原理很容易应用。该有限体积记为 $\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}

上面是一个6个方程6个未知数的系统。（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所著第二版

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

参看[编辑]

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

参考文献[编辑]

^ Fluid Mechanics（Schaum's Series）, M. Potter, D.C. Wiggert, Schaum's Outlines, McGraw-Hill (USA), 2008, ISBN 978-0-07-148781-8
^ Vectors, Tensors, and the basic Equations of Fluid Mechanics, R. Aris, Dover Publications, 1989, ISBN (10) 0-486-66110-5
^ Encyclopaedia of Physics (2nd Edition), R.G. Lerner, G.L. Trigg, VHC publishers, 1991, ISBN (Verlagsgesellschaft) 3-527-26954-1, ISBN（VHC Inc.）0-89573-752-3

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:数学方程的世界

0 0