纳维-斯托克斯方程
性质
非线性
纳维-斯托克斯方程在大多数实际情况下是非线性的偏微分方程。 在某些情况下,一维流和斯托克斯流(或蠕动流)方程可以简化为线性方程组。非线性使得很多的问题很难或者没法解决,但是它却是动荡方程组模型的重要影响因素。非线性是由于对流加速(与点速度变化相关联的加速度),因此,任何对流无论湍流与否都会涉及非线性。
湍流
湍流是时变的无序行为,这种行为常见于许多流体流动中。人们普遍认为,这是由于惯性流体被看作为一个整体而产生的:时间依赖性和对流加速度;因此惯性的影响很小的流体往往是层流(雷诺数量化了流所受惯性的大小)。据不精确了解,纳维-斯托克斯方程能够合理的描述湍流。 纳维-斯托克斯方程关于湍流的数值解是非常难得到的,而且由于湍流的显著不同的混合长度,需要精细分辨率的稳定解的计算时间变得不可行(见直接数值模拟)。尝试使用层流来解决湍流问题的方法通常会导致时间不稳定的解,它不能够适当收敛。为了解决这个问题,时间平均方程,如雷诺平均纳维-斯托克斯方程(RANS),再辅以湍流模型,被用作实际的计算流体动力学模拟湍流时(CFD)的应用程序。例如Spalart-Allmaras湍流,K-ω湍流,k-ε湍流和SST模型。另一种方法解决数值的纳维-斯托克斯方程是大涡模拟(LES)。这种方法在计算上是比RANS方法更昂贵(在时间和计算机存储器上),但因为较大的湍流尺度被明确地解决所以会产生更好的结果。
适用性
补充方程(例如,质量守恒定律)和良好的边界条件,使纳维-斯托克斯方程似乎可以进行流体运动的精确建模;甚至湍流(平均上)也符合实际观察结果。纳维-斯托克斯方程假定所研究的流体是连续(它是无限可分的,而不是由粒子组成),并且是相对静止的。在非常小的尺度或极端条件下,由离散分子组成的真实流体将产生与由纳维-斯托克斯方程建模的连续流体不同的结果。根据Knudsen数的问题,统计力学甚至分子动力学可能是一个更合适的方法。另一个限制是方程的复杂性。对于普遍流体家族的公式是存在的,但纳维-斯托克斯方程对不常见流系的应用,往往会导致非常复杂的构成,而且该配方是目前研究的一个领域。出于这个原因,这些方程通常用于描述牛顿流体。研究这种液体是“简单”的,因为粘度模型最终被线性化;真正描述其他类型流体(如血液)的普遍的模型截至2012年还不存在。
基本假设
在解释纳维-斯托克斯方程的细节之前,我们必须首先对流体的性质作几个假设。第一个假设是流体是连续的。这强调它不包含形成内部的空隙,例如,溶解的气体的气泡,而且它不包含雾状粒子的聚合。另一个必要的假设是所有涉及到的场,全部是可微的,例如压强,速度,密度,温度,等等。
该方程从质量,动量,和能量的守恒的基本原理导出。对此,有时必须考虑一个有限的任意体积,称为控制体积,在其上这些原理很容易应用。该有限体积记为 Ω Ω --> {\displaystyle \Omega } ,而其表面记为 ∂ ∂ --> Ω Ω --> {\displaystyle \partial \Omega } 。该控制体积可以在空间中固定,也可能随着流体运动。这会导致一些特殊的结果,我们将在下节看到。
随质导数
运动流体的属性的变化,譬如大气中的风速的变化,可以有两种不同的方法来测量。可以用气象站或者气象气球上的风速仪来测量。显然,第一种情况下风速仪测量的速度是所有运动的粒子经过一个固定点的速度,而第二种情况下,仪器在测量它随着流体运动时速度的变化。同样的论证对于密度、温度、等等的测量也是成立的。因此,当作微分时必须区分两种情况。第一种情况称为 空间导数 或者 欧拉导数 。第二种情况称为 实质 或 拉格朗日 导数。
随质导数定义为算子(operator):
其中 v {\displaystyle \mathbf {v} } 是流体的速度。方程右边的第一项是普通的欧拉导数(也就是在静止参照系中的导数)而第二项表示由于流体的运动带来的变化。这个效应称为移流(advection)。
L 的守恒定律在一个控制体积上的积分形式是:
因为Ω是共动的,它随着时间而改变,所以我们不能将时间导数和积分简单的交换。
因为这个表达式对于所有 Ω Ω --> {\displaystyle \Omega } 成立,它可以简化为:
对于不是密度的量(因而它不必在空间中积分), D D t {\displaystyle {\frac {\mathrm {D} }{\mathrm {D} t}}} 给出了正确的共动时间导数。
守恒定律
NS方程可以从守恒定律通过上述变换导出,并且需要用状态定律来闭合。
在控制体积上,使用上述变换,下列的量视为守恒:
质量
能量
动量
角动量
连续性方程
质量的守恒写作:
其中
在不可压缩流体的情况 ρ ρ --> {\displaystyle \rho } 不是时间或空间的函数。方程简化为:
动量守恒
动量守恒写作:
注意 v ⊗ ⊗ --> v {\displaystyle \mathbf {v} \otimes \mathbf {v} } 是一个张量, ⊗ ⊗ --> {\displaystyle \otimes } 代表张量积。
我们可以进一步简化,利用连续性方程,这成为:
我们可以认出这就是通常的 F =m a 。
方程组
一般形式
方程组的形式
纳维-斯托克斯方程的一般形式是:
关于动量守恒。张量 P {\displaystyle \mathbb {P} } 代表施加在一个流体粒子上的表面力(应力张量)。除非流体是由象旋涡这样的旋转自由度组成, P {\displaystyle \mathbb {P} } 是一个对称张量。一般来讲,我们有如下形式:
其中 σ σ --> {\displaystyle \sigma } 是法向约束,而 τ τ --> {\displaystyle \tau } 是切向约束。
迹 σ σ --> x x + σ σ --> y y + σ σ --> z z {\displaystyle \sigma _{xx}+\sigma _{yy}+\sigma _{zz}} 在流体处于平衡态时为0。这等价于流体粒子上的法向力的积分为0。
我们再加上连续性方程:
对于 处于平衡 的液体, P {\displaystyle \mathbb {P} } 的迹是3 p 。
其中
最后,我们得到:
其中 T {\displaystyle \mathbb {T} } 是 P {\displaystyle \mathbb {P} } 的非对角线部分。
闭合问题
这些方程是不完整的。要对它们进行完备化,必须对 P {\displaystyle \mathbb {P} } 的形式作一些假设。例如在理想流体的情况 τ τ --> {\displaystyle \tau } 分量为0。用于完备方程组的方程是状态方程。
再如,压强可以主要是密度和温度的函数。
要求解的变量是速度的各个分量,流体密度,静压力,和温度。流场假定为可微并连续,使得这些平衡得以用偏微分方程表达。这些方程可以转化为涡度和流函数这些次变量的威尔金森方程组。解依赖于流体的性质(例如粘滞度、比热、和热导率),并且依赖于所研究的区域的边界条件。
P {\displaystyle \mathbb {P} } 的分量是流体的一个无穷小元上面的约束。它们代表垂直和剪切约束。 P {\displaystyle \mathbb {P} } 是对称的,除非存在非零的自旋密度。
所谓非牛顿流体是就是其中该张量没有特殊性质使得方程的特殊解出现的流体
特殊形式
这些是问题的特定的常见简化,有时解是已知的。
牛顿流体
在牛顿流体中,如下假设成立:
其中
其中为简化书写,对脚标使用了爱因斯坦求和约定。
不采用简化书写的完整形式非常繁琐,分别为:
动量守恒:
质量守恒:
因为 密度 是一个未知数,我们需要另一个方程。
能量守恒:
其中:
假设一个理想气体:
上面是一共6个方程6个未知数的系统。(u,v,w,T,e以及 ρ ρ --> {\displaystyle \rho } )。
宾汉(Bingham)流体
在宾汉流体中,我们有稍微不同的假设:
那些流体在开始流动之前能够承受一定的剪切。牙膏是一个例子。
幂律流体
这是一种理想化的流体,其剪切应力, τ τ --> {\displaystyle \tau } ,由下式给出
该形式对于模拟各种一般流体有用。
不可压缩流体
其纳维-斯托克斯方程(Navier-Stoke equation)分为动量守恒公式
和质量守恒公式
其中,对不可压缩牛顿流体来说,只有对流项(convective terms)为非线性形式。对流加速度(convective acceleration)来自于流体流动随空间之变化所产生的速度改变,例如:当流体通过一个渐缩喷嘴(convergent nozzle)时,流体产生加速之情况。由于此项的存在,对于暂态运动中的流体来说,其流场速度变化不再单是时间的函数,亦与空间有关。
另外一个重要的观察重点,在于黏滞力(viscosity)在流场中的以流体速度作拉普拉斯运算来表现。这暗示了在牛顿流体中,黏滞力为动量扩散(diffusion of momentum),与热扩散方程式非常类似。
若 μ μ --> {\displaystyle \mu } 在整个流体上均匀,动量方程简化为
(若 μ μ --> = 0 {\displaystyle \mu =0} 这个方程称为欧拉方程;那里的重点是可压冲击波冲击波)。
如果现在再有 ρ ρ --> {\displaystyle \rho } 为常数,我们得到如下系统:
连续性方程(假设不可压缩性):
注意纳维-斯托克斯方程仅可近似描述液体流,而且在非常小的尺度或极端条件下,由离散的分子和其他物质(例如悬浮粒子和溶解的气体)的混合体组成的真实流体,会产生和纳维-斯托克斯方程所描述的连续并且齐性的液体不同的结果。依赖于问题的纳森数,统计力学可能是一个更合适的方法。但是,纳维-斯托克斯方程对于很大范围的实际问题是有效的,只要记住他们的缺陷是天生的就可以了。
参看
雷诺数
马赫数
雷诺平均纳维-斯托克斯方程
参考文献
^ 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
免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。
- 有价值
- 一般般
- 没价值