纳维-斯托克斯方程(N-S方程)详细推导
纳维-斯托克斯方程(N-S方程)详细推导
针对非牛解方法 ，以揭示其复杂的流动行为和机理。
THANKS FOR WATCHING
感谢您的观看
N-S方程的改进和发展
数值方法
为了解决N-S方程的求解问题， 研究者们发展出了许多数值方法，
如有限差分法、有限元法、谱方 法等。
近似模型
针对某些特定流动，研究者们提出 了许多近似模型，如雷诺平均N-S 方程、湍流模型等，以简化求解过 程。
多物理场耦合
随着计算技术的发展，多物理场耦 合成为研究流体流动的重要方向， 如流固耦合、流热耦合等。
应力张量
01
应力张量是描述流体内部应力的二阶张量，包括正应力和剪切 应力。
02
正应力表示流体在单位面积上受到的压力，而剪切应力表示流
体在单位面积上受到的切向力。
应力张量是流体的状态函数，其值取决于流体的状态和所处的
03
边界条件。
03 纳维-斯托克斯方程的推 导
纳维方程的推导
01
02
03
从质量守恒、动量守恒 和牛顿第二定律出发， 推导出描述流体运动的
考虑流体的粘性和惯性
02
N-S方程中包含了流体的粘性和惯性力，能够描述粘性流体在运
动过程中的受力情况和运动规律。
涉及到复杂的数学处理
03
N-S方程的推导涉及到复杂的数学处理，包括微积分、线性代数
和偏微分方程等。
02 流体的基本性质
流体的定义和分类
流体是能够流动的物质，具有连续性和 不可压缩性。根据其流动特性，流体可 分为牛顿流体和非牛顿流体两大类。
04 N-S方程的应用和限制
N-S方程的应用领域
流体力学
N-S方程是描述流体运动的基本方程，广泛应用于航空、航海、 气象、环境等领域。
纳维-斯托克斯方程(Navier-Stokes equations)
纳维-斯托克斯方程(Navier-Stokes equations)名称由来Navier-Stokes equations描述粘性不可压缩流体动量守恒的运动方程。
简称N-S方程。
因1821年由C.-L.-M.-H.纳维和1845年由G.G.斯托克斯分别导出而得名。
该方程是可压缩流体的N-S方程。
其中，Δ是拉普拉斯算子；ρ是流体密度；pN-S方程意义后人在此基础上又导出适用于可压缩流体的N-S方程。
N-S方程反映了粘性流体（又称真实流体）流动的基本力学规律，在流体力学中有十分重要的意义。
它是一个非线性偏微分方程，求解非常困难和复杂，目前只有在某些十分简单的流动问题上能求得精确解；但在有些情况下，可以简化方程而得到近似解。
例如当雷诺数Re1时，绕流物体边界层外，粘性力远小于惯性力，方程中粘性项可以忽略，N-S方程简化为理想流动中的欧拉方程（=－Ñp+ρF）；而在边界层内，N-S方程又可简化为边界层方程，等等。
在计算机问世和迅速发展以后，N-S方程的数值求解才有了很大的发展。
基本假设在解释纳维-斯托克斯方程的细节之前，首先，必须对流体作几个假设。
第一个是流体是连续的。
这强调它不包含形成内部的空隙，例如，溶解的气体的气泡，而且它不包含雾状粒子的聚合。
另一个必要的假设是所有涉及到的场，全部是可微的，例如压强P，速度v，密度，温度Q，等等。
该方程从质量，动量，和能量的守恒的基本原理导出。
对此，有时必须考虑一个有限的任意体积，称为控制体积，在其上这些原理很容易应用。
该有限体积记为\Omega，而其表面记为\partial\Omega。
该控制体积可以在空间中固定，也可能随着流体运动。
纳维-斯托克斯方程（Navier-Stokes equations），以克劳德-路易·纳维(Claude-Louis Navier)和乔治·盖伯利尔·斯托克斯命名，是一组描述象液体和空气这样的流体物质的方程。
流体力学-N-S方程
实际流体的运动微分 方程
——纳维-斯托克斯方程式 （N-S方程式）
以应力表示的黏性流体运动微分方程式
• 一、作用在流体微元上的应力 在粘性不起作用的平衡流体 中，或者在没有粘性的理想运动 流体中，作用在流体微元表面上 的表面力只有与表面相垂直的压 应力，而且压应力又具有一点上 各向同性的性质。
图一
v x x v y
（6）
由式（6）可以看出，由于各个方向的直线应变速 度不见得相等，因而这种由于粘性阻碍作用所产生的 法向应力也是各向不等的，p'xxp'yyp'zz统称为一点上的 各项异性压强。 • 于是在实际流体运动时，一点上的法向应力除了由 于分子运动统计平均的各向同性压强p之外，还需加上 由于粘性影响而与直线变形有关的各向异性压强，最 后可以得到法向应力与直线应变速度之间的关系为
（9）
此式说明一点上的各向同性压强也就是不可 压缩实际流体中不同方向压强的算术平均值。这 给具体计算实际流体中的压强带来很大的方便， 我们无需进一步研究各向异性压强，只要找出各 向同性压强与其他流动参数之间的关系，则据此 算出的各向同性压强事实上也就是不可压缩实际 运动流体一点上的流体动压强。
p的含义
但是在运动着的实际流体中取出边长dx、dy、 dz的六面体微元，如右图1多示，由于粘性影响，当 微元有剪切变形时，作用在微元体ABCDEFGH上的表 面力就不仅有压应力p，而且也有切应力τ 。当微元 有直线变形时，一点上的压应力也不再具有各项同 性的性质了。
纳维斯托克方程
纳维-斯托克斯方程（Navier-Stokes equations），简称N-S 方程，是一组描述流体运动的方程。
它们由克劳德-路易·纳维和乔治·斯托克斯提出，可以看作是流体运动的牛顿第二定律。
这些方程建立了流体的粒子动量的改变率（力）和作用在液体内部的压力的变化和耗散粘滞力（类似于摩擦力）以及引力之间的关系。
这些粘滞力产生于分子的相互作用，能告诉我们液体有多粘。
这样，纳维-斯托克斯方程描述作用于液体任意给定区域的力的动态平衡。
纳维-斯托克斯方程总是要与连续性方程同时进行求解：纳维-斯托克斯方程表示动量守恒，而连续性方程则表示质量守恒。
由于这些方程本身的复杂性，我们只能得到非常有限的解析解。
不过，在特定的边界条件（如入口、出口和壁）下，可以求解这些方程来预测给定几何体中的流体速度和压力。
它们被广泛应用于模拟天气、洋流、管道中的水流、星系中恒星的运动、翼型周围的气流等领域，也被用于飞行器和车辆的设计、血液循环的研究、电站的设计、污染效应的分析等。
纳维-斯托克斯方程(N-S方程)详细推导
• 线变形运动 微团左、右两侧的 A 点和 C 点沿 x 方向的速
度差为
，当这速度差值为正时，微团沿 x 方向发生
伸长变形；当它为负时，微团沿 x 方向发生缩短变形。
• 线变形速度 单位时间，单位长度的线变形称为线变形速 度。流体微团沿 x 方向的线变形速度：
Hale Waihona Puke 本构方程和NS方程粘性流体动力学基础
旋转角速度 把对角线的旋转角速度定义为整个流体微团在平 面上的旋转角速度。
dxdydzdt
t
dxdydzdt
或：
(vx ) (vy ) (vz ) 0
t x
y
z
连续性方程
矢量形式：
r
g()
0
t
（适用于层流、湍流、 牛顿、非牛顿流体）
连续方程物理意义：流体在单位时间内流经单位体积空间输 出与输入的质量差与其内部质量变化的代数和为零。
本构方程和NS方程
粘性流体动力学基础
本构方程和NS方程
方程的物理意义：
粘性流体动力学基础
方程左边是：任意时刻t通过考察点A的流体质点
加速度的三个分量；
Dvx / Dt ax
方程右边是：作用在单位体积流体上的表面力和体 积力在各坐标上的分量。
方程可简略表示成：
r ur
a F
这就是以单位体积的流体质量为基准的牛顿第二运动定律
本构方程和NS方程
dy
zx
zx z
dz
xx yx
z y xz
yz
xy
xy x
dx
xx
xx x
dx
x
zx
zz
yx
yx y
dy
应力状态：
N-S方程推导
作用在隔离流体（也就是所取的研究流体单元）的表面，和作用 的面积成正比的力。分为垂直于作用面的压力和沿作用面方向的切 力。表面力可以使作用于流体界面的压力、切力，也可以是一部分流 体质点作用于相邻另一部分流体质点的压力、切力。单位作用面的压 （第一个下标表示作用面的法线方向， 应力、 切应力即为图 1 中的 σ 、τ 第二个下标表示力的方向）。 以 x 方向为例，流体单元受到的力：
2 əτzy dz — τzy+ —— əz 2 əτyz dy — τyz+ —— əy 2 əσyy dy — σyy+ —— əy 2
y
əτxy dx — τxy+ —— əx 2
əτyx dy — τyx+ —— əy 2
x
图 1 作用在单元体上的力 作用力有两类，即质量力和表面力。
1.1 质量力
( )
( )
( )
( )
( )
（22） 同理： r du ∂p 1 ∂ ρ x = − + µ∇ 2u x + µ ∇ u + ρ X （23） dt ∂x 3 ∂x
( )
ρ
ρ
du y dt
=−
r ∂p 1 ∂ + µ∇ 2u y + µ ∇ u + ρY （24） ∂y 3 ∂y
( )
r du z ∂p 1 ∂ = − + µ∇ 2u z + µ ∇ u + ρ Z （25） dt ∂z 3 ∂z 矢量形式： v v 1 v uv du ρ = −∇p + µ∇ 2 u + µ∇ ∇ u + ρ F （26） dt 3
引言
N-S方程讲解
黏性流体动量平衡方程−纳维-斯托克斯方程（Navier-Stokes equations ） 1.动量平衡的定义流体在流动过程中遵守能量守恒定律，称为能量平衡根据牛顿第二定律：⎩⎨⎧≠∑=∑，运动，动力平衡，静止，静力平衡0F 0F 作用力的合力 = 单位时间内动量的变化量作用力形式 动量形式[动量传入量] - [动量传出量] +[系统作用力的总和] = [动量蓄积量][动量传入量] - [动量传出量] + [系统作用力的总和] = 0稳定流动系统：不稳定流动系统：动量收支差量动量收支差量⒉ 动量传递方式1 黏性动量传输dydv x yx μτ-= 2 对流动量传输对流动量传输vvρ⒊ 作用力的形式体积力表面力压力重力作用力⒋ 动量平衡方程的推导元体分析法牛顿第二定律分析法建立方法建立依据在直角坐标系中由于有三个方向的分速度，所以共有九个动量通量。
⎪⎩⎪⎨⎧⋅⋅⋅⋅⋅⋅⋅⋅⋅zz yz xz z y y y x y z x y x xx v v v v v v v v v v v v v v v v v ρρρρρρρρρv 以v x动量通量收支差量⑴ 对流动量收支差量x 方向的速度、x 方向的动量通量对流动量收支差量为同理，以v x 为准，y 方向、z 方向的对流动量收支差量：以v x 为准，元体对流动量收支差量为同理，以v y 、v z 为准，元体对流动量收支差量为 v x → v y 、v z⑵黏性动量收支差量黏性动量通量同样由九个分量组成以v x为准，C、D黏性动量通量收支差量黏性动量收支差量同理，v x在y、x以v x为准，元体黏性动量收支差量为同理，以⑶作用力的总和zxgxddydρzxgyddydρzxgzddydρx方向：P Ax方向合压力为x方向的总压力为同理，y、z方向的总压力为x →y、z重力⑷ 动量蓄积量z 方向x 方向y 方向 单位时间内元体动量的变化量[动量传入量] - [动量传出量] +[系统作用力的总和] = [动量蓄积量]⒌ 动量平衡方程式将以上式子代入下式，整理得：N-S 方程简化：const=ρ，连续性方程⑵const=μ，牛顿黏性定律⑴动量收支差xx x x x x x x g x p zv y v x v z v v y v v x v v v ρμτρ+∂∂-∂∂+∂∂+∂∂=∂∂+∂∂+∂∂+∂∂)()(222222z y x yy y y y y y y g xp zv yv xv zv v yv v xv v v ρμτρ+∂∂-∂∂+∂∂+∂∂=∂∂+∂∂+∂∂+∂∂)()(222222zyxzz z z z z z z g xp zv yv xv zv v yv v xv v v ρμτρ+∂∂-∂∂+∂∂+∂∂=∂∂+∂∂+∂∂+∂∂)()(222222zyx黏性力引起压力 体积力积累动量收支差量⒍ 动量平衡方程的讨论x2x 22x 22x 2x z x y x x x g x P z v y v x v z v y v v x v v ρμτρ+∂∂-⎥⎦⎤⎢⎣⎡∂∂+∂∂+∂∂=⎥⎦⎤⎢⎣⎡∂∂+∂∂+∂∂+∂∂v v 对流动量动量蓄积量黏性动量压重（1）方程的物理意义：运动的流体能量守恒的表现⎩⎨⎧作用力形式动量形式z zv y y v x x v v v d d d d d ⋅∂∂+⋅∂∂+⋅∂∂+⋅∂∂=ττz y x v zv v y v v x v v ⋅∂∂+⋅∂∂+⋅∂∂+∂∂=v d d ττz y x v zvv y v v x v v a ⋅∂∂+⋅∂∂+⋅∂∂+∂∂=τz x y x x x x x v zvv y v v x v v a ⋅∂∂+⋅∂∂+⋅∂∂+∂∂=τ全微分)z ,y ,x ,(v v τ=x2x 22x 22x 2g x Pz v y v xv ρμρ+∂∂-⎥⎦⎤⎢⎣⎡∂∂+∂∂+∂∂=xa y 2y 22y 22y2g Pz v y v xv ρμρ+∂∂-⎥⎥⎦⎤⎢⎢⎣⎡∂∂+∂∂+∂∂=ya y z 2z 22z 22z 2g Pz v y v xv ρμρ+∂∂-⎥⎦⎤⎢⎣⎡∂∂+∂∂+∂∂=z az惯性力黏性力压力重力流体在运动中以作用力及动量形式表现能量平衡 关系是统一的⑵ 适用条件黏性流体、不稳定流动、不可压缩流体(元体范围内)、层流流动理想流体：=μ没有黏性的流体简化： 0v =∂∂τ② 稳定流动， ③ 单位质量流体 0=μ①时，N-S 方程简化为欧拉方程理想流体、稳定流动、不可压缩流体(元体范围内)流动微分方程的应用求解步骤(1)根据问题特点对一般形式的运动方程进行简化，获得针对具体问题的微分方程或方程组。
N-S(纳维斯托克斯)方程推导过程
很多人一听到N-S 方程就有点头皮发麻，因为涉及到流体力学的知识比较多，如果没有一个完整有逻辑的思路，理解N-S 方程是有点困难。
其中涉及到欧拉法，场论，随体导数，流体力学连续性方程（即质量守恒方程），流体力学N-S 方程（即动量方程），动量方程在流体力学中有两种，一种是理想流体动量方程，一种是粘性流体动量方程，粘性流体的动量方程也叫纳维-斯托克斯方程，也简称N-S 方程。
我试图想把N-S 方程弄清楚点，所以写了一点东西，分享一下。
首先要讲一下流体力学的欧拉法，在课本中还讲了拉格朗斯法，因为连续性方程和N-S 方程是用欧拉法得出的，和拉格朗日法没什么关系。
我就不讲拉格朗日法，以免产生混乱。
欧拉方法的着眼点不是流体质点而是空间点。
设法在空间中的每一点上描述出流体运动随时间的变化状况。
如果每一点的流体运动都已知道，则整个流体的运动状况也就清楚了。
欧拉方法中流体质点的运动规律数学上可表示为下列矢量形式：假设空间一点的坐标(x,y,z,t)，其中x,y,z 是该空间的坐标，t 是此刻时间。
u,v,w 是这一空间点的三个方向速度。
p,ρ,T 是这一空间点的压力，密度和温度。
这样就有了每一个点的速度，压力，密度，温度，就可以描述运动流体的状态。
这里需要强调一点的是下面这六个式子，可以换一个角度把他们看成方程，对后面理解连续性方程和N-S 方程有帮助，比如u=x+2y+3z),,,();,,,();,,,();,,,();,,,();,,,(t z y x T T t z y x t z y x p p t z y x w w t z y x v v t z y x u u ======ρρ因为后面需要随体导数的概念，还需要把速度函数表示成矢量的形式。
前面u,v,w 是标量，是ν在(x,y,z,t)直角坐标系三个方向的速度。
),(t rνν=M 点（x,y,z,t ），速度为),(t M ν ，过了t ∆之后，在M '点，速度为),(t t M ∆+'ν。
NS方程推导
代入上面加速度公式，得到
d dt
(M ,t) t
t0,M和M 靠近，
MMMlMi的m变0化会(引M起,三tM个) M方向速(M度的,变t)化
用M点速度
du dt
u(x, y, t
z,t)
u(x,
y,
z,t)
u(x, y, x
z,t)
v(x,
y, z,t)
u(x, y, z,t) y
w(
x,
y,
z,
t
)
u(
x, y, z
z,
t
)
u u u v u w u t x y z
du dt
u t
u
u x
v
u y
w
u z
;
dv dt
v t
u
v x
v
v y
w
v z
;
dw w u w v w w w dt t x y z
至此已经用欧拉法推到出了流体速度和加速度(即随体导数)的公式。随体导数也可以用复合
以上就已经得到了连续性方程。 对不可压缩流体，连续性方程可以简化，可以得到以下简化的连续性方程：
u x
v y
பைடு நூலகம்
纳维-斯托克斯方程(N-S方程)详细推导
yx xy
yz zy
zx xz
16
本构方程和NS方程
粘性流体动力学基础
微元体表面力的总力分量
X方向的表面力：
x dxdydz yx dydxdz zx dzdxdy
x
y
z
x yx zx dxdydz
x y z
Y方向的表面力：
xy
x
yy
y
zy
z
dxdydz
Z方向的表面力：
xz
x
yz
y
zz
z
dxdydz
17
本构方程和NS方程
动量流量及动量变化率
粘性流体动力学基础
z
vz vx
vz vx z
dz
dy
vyvx
vy vx y
dy
dx
动量流量
动量通量 x 流通面积
vx vx
dz
vxvx
vxvx x
dx
＝ 动量流量
y
vyvx vzvx
vz x
vx z
24
本构方程和NS方程
本构方程的讨论：
正应力与线变形速率：
线变形率与流体流动：
正应力中的粘性应力：
粘性流体动力学基础
流体正应力与三个速度偏导数有关 (即:线变形率)，同固体力学中的虎 克定律。
从流体流动角度看，线变形率的正负 反映了流体的流动是加速还是减速； 体变形率的正负反映了流动过程中流 体体积是增加还是减少。
粘性流体动力学基础
牛顿流体的本构方程：
xx
p 2
x
x
2 3
vx x
vy y
vz z
yy
p 2
【精编】纳维-斯托克斯方程(N-S方程)详细推导
一维流动的连续方程 1 A1 2 A2
连续性方程是流体流动微分方程最基本的方程 之一。任何流体的连续运动均必须满足。
本构方程和NS方程
粘性流体动力学基础
( vy ) ( vx ) ( vz ) dxdydzdt dxdydzdt dxdydzdt x y z ( vx ) ( v y ) ( vz ) dxdydzdt y z x
本构方程和NS方程
与微团内各点速度的变化有关。 设方形流体微团中心 M 的流速
分量为 ux 和 uy ，则微团各侧边
的中点 A 、 B 、 C 、 D 的流速 分量分别为：
微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。
本构方程和NS方程
粘性流体动力学基础
• 平移运动速度 微团上各点公有的分速度 ux 和uy ，使它们 在 dt 时间内均沿 x 方向移动一距离 uxdt , 沿 y 方向移动一 距离 uydt 。因而，把中心点 M 的速度 ux和 uy ，定义为流 体微团的平移运动速度。
1 u z u y u y u x ( ) y ( x z ) 2 y z 2 z x
1 u y u x z ( ) 2 x y
本构方程和NS方程
粘性流体动力学基础
亥姆霍兹速度分解定理
整理推 广得
本构方程和NS方程
粘性流体动力学基础
不可压缩流体连续性微分方程
理想流体的运动微分方程
理想流体运动微分方程式是研究流体运动学的重要理论基 础。可以用牛顿第二定律加以推导。 a F 受力分析：
纳维-斯托克斯方程(N-S方程)详细推导ppt课件
dy
zx
zx z
dz
xx yx
z y xz
Z方向：
( v y
y
) dxdydzdt
2、dt时间内，整个六面体内输入与输出的质量差：
(vx ) dxdydzdt (vy ) dxdydzdt (vz ) dxdydzdt
x
y
z
( vx
x
)
(vy
y
)
( vz
z
)
dxdydzdt
10
本构方程和NS方程
粘性流体动力学基础
3、微元体内的质量变化： dxdydzdt
12
本构方程和NS方程
粘性流体动力学基础
理想流体的运动微分方程
理想流体运动微分方程式是研究流体运动学的重要理论基
础。可以用牛顿第二定律加以推导。
受力分析：
r ur
a F
1、质量力： fxρdxdydz x轴正方向
2、表面力：
切向应力＝0（理想流体） 法向应力＝压强
p p dx p p dx
微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。5
本构方程和NS方程
粘性流体动力学基础
• 平移运动速度 微团上各点公有的分速度 ux 和uy ，使它们 在 dt 时间内均沿 x 方向移动一距离 uxdt , 沿 y 方向移动一 距离 uydt 。因而，把中心点 M 的速度 ux和 uy ，定义为流 体微团的平移运动速度。
3
本构方程和NS方程
粘性流体动力学基础
平移运动、旋转运动、线变形运动和角变形运动
右图为任意t时刻在平面流场中所取的一个正方形流体微团。由 于流体微团上各点的运动速度不一致，经过微小的时间间隔后， 该流体微团的形状和大小会发生变化，变成了斜四边形。
纳威—斯托克斯(N—S)方程推导的专题教学改革探索
纳威—斯托克斯（N—S）方程推导的专题教学改革探索作者：吴磊马孝义来源：《课程教育研究·下》2014年第11期【摘要】水力学是一门技术基础科学，它是力学的一个分支。
N-S（Navier-Stokes）方程是水力学课程的重要内容，是联系本科生水力学与研究生高等流体力学课程的关键环节。
但现实授课过程中经常听到学生反映很难理解N-S方程，并且觉得没有应用价值，原因是本科教学中没有用其做过习题。
因此，为了更好地让学生理解N-S方程的原理，掌握方程的推导过程，我们对此开展专题性教学改革探索，进而推进水力学重、难点的教学改革与创新，使学生更愿意学，更容易懂，更能真正理解和应用。
【关键词】水力学 N-S方程研究型本科专题教学教学探索【中图分类号】G642.0 【文献标识码】A 【文章编号】2095-3089（2014）11 -0060-02水力学是高等工科院校水利、土木、环境等诸多专业都需要掌握的一门重要技术基础课，是学生专业能力形成与未来职业发展必不可少的依托[1]。
水力学的任务主要是研究液体（主要是水）的平衡和机械运动规律及其工程应用[2]，是一门既有较强理论性又有较强工程实际意义的课程，具有理论推导复杂，概念多且易混淆特点[3]。
另外，水力学具有对学生的高等数学、大学物理以及综合分析问题的能力要求较高的特点[4]。
作为专业技术基础，它不仅为后续专业课学习提供理论基础，而且注重理论联系实际[5]。
比如用于描述粘性不可压缩流体动量守恒的纳维-斯托克斯方程（N-S方程），是不可压缩流体动力学最有用的一组方程之一，也是学习高等流体力学、计算流体力学等课程的重要基础[6]。
N-S方程可用于描述与分析大量对学术和经济有用现象的物理过程，包括模拟管道中的水流、翼型周围的气流、污染效应的分析、电站的设计、天气、洋流、星系中恒星的运动、飞行器和车辆的设计、血液循环的研究等[7]。
因此，考虑到N-S方程在上述专业领域本科教学过程中的重要地位，本文从专题教学的角度积极探索本科生容易接受的N-S方程推导方法，把握水力学研究型本科教学内涵，着力提升学生的理论基础和实践能力，为学生创新精神和研究能力的培养提供可靠支撑。
纳维-斯托克斯方程(N-S方程)详细推导
粘性流体动力学基础
本构方程及N-S方程
李连侠
水力学与山区河流开发保护国家重点实验室 2009年4月
本构方程和NS方程
粘性流体动力学基础
内容提要
• 流体运动分析及理想流体基本方程 • 真实流体受力分析 • 利用张量理论推导本构方程和粘性流体力学基本方程
本构方程和NS方程
粘性流体动力学基础
vz vx z dz
v y vx y
z
dy
v y vx
dy
动量流量
动量通量
dx dz
vx vx
vx vx
vx vx x
x
流通面积
dx
＝ 动量流量
图中标注的是动量的输入或 输出方向，而动量或其通量 本身的方向均指向 x方向，即 分速度vx的方向。
粘性流体动力学基础
3、微元体内的质量变化： dxdydzdt
t
从而有：
( vx ) ( v y ) ( vz ) dxdydzdt dxdydzdt y z t x
或：
( vx ) ( vy ) ( vz ) 0 连续性方程 t x y z
与微团内各点速度的变化有关。 设方形流体微团中心 M 的流速
分量为 ux 和 uy ，则微团各侧边
的中点 A 、 B 、 C 、 D 的流速 分量分别为：
微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。
本构方程和NS方程
粘性流体动力学基础
• 平移运动速度 微团上各点公有的分速度 ux 和uy ，使它们 在 dt 时间内均沿 x 方向移动一距离 uxdt , 沿 y 方向移动一 距离 uydt 。因而，把中心点 M 的速度 ux和 uy ，定义为流 体微团的平移运动速度。
纳维-斯托克斯方程详细推导
本构方程和NS方程
粘性流体动力学基础
平移运动、旋转运动、线变形运动和角变形运动 右图为任意t时刻在平面流场中所取的一个正方形流体微团。由 于流体微团上各点的运动速度不一致，经过微小的时间间隔后， 该流体微团的形状和大小会发生变化，变成了斜四边形。
本构方程和NS方程
粘性流体动力学基础
流体微团的运动形式
粘性流体动力学基础
1、x方向：dt时间内沿从六面体 x 处与 x+dx 处输入与输出的 质量差：
( v ) ( v ) x x v d y d z d t v d x d y d z d t d x d y d z d t x x x x
流体质点运动的分析
•分析流场中任意流体微团运动是研究整个流场运动的基础。 •流体运动要比刚体运动复杂得多，流体微团基本运动形式 有平移运动、旋转运动、线变形和角变形运动等。实际运 动也可能遇到只有其中的某几种形式所组成。 •当流体微团无限小而变成质点时，其运动也是由平动、线 变形、角变形及旋转四种基本形式所组成。
粘性流体动力学基础
理想流体的运动微分方程
根据牛顿第二定律得x轴方向的运动微分方程
du p d x p d x x f d x d y d z p d y d z p d y d z d x d y d z x x 2 x 2 d
与微团内各点速度的变化有关。 设方形流体微团中心 M 的流速
分量为 ux 和 uy ，则微团各侧边
的中点 A 、 B 、 C 、 D 的流速 分量分别为：
微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
xy
xy x
dx
粘性流场中任意一点的应
xx
xx x
dx
力有9个分量，包括3个正
应力分量和6个切应力分量：
x
zx
zz
yx
yx y
dy
切应力互等定律
微元体上X和Z方向的表 在6个切应力分量中，互换
面力
下标的每一对切应力是相
yx xy
yz zy
zx 等 x的z 。
微元体表面力的总力分量
线变形运动 微团左、右两侧的 A 点和 C 点沿 x
方向的速度差为
，当这速度差值为正时
，微团沿 x 方向发生伸长变形；当它为负时，微
团沿 x 方向发生缩短变形。
线变形速度 单位时间，单位长度的线变形称为线 变形速度。流体微团沿 x 方向的线变形速度：
旋转角速度 把对角线的旋转角速度定义为整个流体微团在平 面上的旋转角速度。
vx
dydzdt
vx
(vx
x
)
dx
dydzdt
(vx
x
)
dxdydzdt
Y方向： (
v
z
z
)
dxdydzdt
；
Z方向
( v y
y
) dxdydzdt
：
2、dt时间内，整个六面体内输入与输出的质量差：
(vx ) dxdydzdt (vy ) dxdydzdt (vz ) dxdydzdt
x
y
当流体微团无限小而变成质点时，其运动也是由 平动、线变形、角变形及旋转四种基本形式所组成 。
平移运动、旋转运动、线变形运动和角变形运动
右图为任意t时刻在平面流场中所取的一个正方形流体 微团。由于流体微团上各点的运动速度不一致，经过微 小的时间间隔后，该流体微团的形状和大小会发生变化 ，变成了斜四边形。
ux t
ux x
ux
ux y
uy
ux z
uzt
u y x
ux
u y y
uy
u y z
uz
Y
1
p y
uz t
uz x
ux
uz y
uy
uz z
uz
Z
1
p z
理想流体的运动微分方程
即欧拉运动微分方程
粘性流体的运动微分方程
以流体微元为分析对象，流体的运动方程可 写为如下的矢量形式：
DV F P
z
( vx
x
)
(vy
y
)
( vz
z
)
dxdydzdt
3、微元体内的质量变化： dxdydzdt
t
从而有：
( vx
x
)
(vy
y
)
( vz
z
)
dxdydzdt
t
dxdydzdt
或： (vx ) (vy ) (vz ) 0
t x
y
z
连续性方程
矢量形式：
()
t
0
（适用于层流、湍流、 牛顿、非牛顿流体）
X方向的表面力：
x
x
dxdydz
yx
y
dydxdz
zx
z
dzdxdy
x yx zx dxdydz
x y z
Y方向的表面力： xy
x
yy
y
zy
z
dxdydz
Z方向的表面力：
xz
x
yz
y
zz
z
dxdydz
动量流量及动量变化率
z
vz vx
vz vx z
流体微团的运动形式 与微团内各点速度的变化有关。 设方形流体微团中心 M 的流速 分量为 ux 和 uy ，则微团各侧边 的中点 A 、 B 、 C 、 D 的流速 分量分别为：
微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。
平移运动速度 微团上各点公有的分速度 ux 和 uy ，使它们在 dt 时间内均沿 x 方向移动一距离 uxdt , 沿 y 方向移动一距离 uydt 。因而，把中 心点 M 的速度 ux和 uy ，定义为流体微团的平移 运动速度。
本构方程及N-S方程
李连侠
水力学与山区河流开发保护国家重点实验室 2009年4月
内容提要
流体运动分析及理想流体基本方程 真实流体受力分析 利用张量理论推导本构方程和粘性流体力学基本
方程
流体质点运动的分析
分析流场中任意流体微团运动是研究整个流场运 动的基础。
流体运动要比刚体运动复杂得多，流体微团基本 运动形式有平移运动、旋转运动、线变形和角变形 运动等。实际运动也可能遇到只有其中的某几种形 式所组成。
连续方程物理意义：流体在单位时间内流经单位体 积空间输出与输入的质量差与其内部质量变化的代
数和为零。
若流体不可压缩：vx vy vz 0
x y z
上式表明，对于不可压缩液体，单位时间单位体积 空间内流入与流出的液体体积之差等于零，即液体
体积守恒。 适用范围： 恒定流或非恒定流；理想液体或实际液
整理 推广
得
不可压缩流体连续性微分方程
直角坐标系中的连续性方程质量守恒
z
dy
vx dydz
dz
dx
输入微元体－ 输出微元体
的质量流量 的质量流量
vx
vx
x
dx
dydz
＝ 微元体内的
质量变化率
x
y
微元体及其表面的质量
通量
1、x方向：dt时间内沿从六面体 x 处与 x+dx 处输入
与输出的质量差：
体。 一维流动的连续方程1 A1 2 A2
连续性方程是流体流动微分方程最基本的方程 之一。任何流体的连续运动均必须满足。
理想流体的运动微分方程
理想流体运动微分方程式是研究流体运动学的重
要理论基础。可以用牛顿第二定律加以推导。
受力分析：
a F
1、质量力f：xρdxdydz x轴正方向
2、表面力： 切向应力＝0（理想流体） 法向应力＝压强
p p dx p p dx
x 2
x 2
x轴正方向
x轴负方向
理想流体的运动微分方程
根据牛顿第二定律得x轴方向的运动微分方程
f x dxdydz
p
p x
dx 2
dydz
p
p x
dx 2
dydz
dxdydz
dux dt
fx
1
p x
dux dt
fy
1
p y
duy dt
fz
1
p z
duz dt
；
；
x
1 ( uz 2 y
u y z
)
y
1 ( ux 2 z
uz ) x
z
1 ( u y 2 x
u x y
)
•角变形速度：直角边 AMC （或BMD）与对角线 EMF 的夹角的变形速度
x
1 ( uz 2 y
u y ) z
y
y ( ux 2 z
uz ) x
z
1
u (
y
2 x
u x ) y
亥姆霍兹速度分解定理
dz
dy dx
vyvx
(1)
Dt
这里 ：
DV V V V
Dt t
(2)
是流体微团D的Dt 加t 速V 度 ，t 微 Vi 分xi 符号：
(3)
称为物质导数或随体导数，它表示流体微团的 某性质
应力状态及切应力互等定律
zz
zz z
dz
应力状态：
yz
yz y
dy
zx
zx z
dz
xx yx
z y xz
yz