第二章 连续性方程与运动方程
合集下载
流体力学连续性方程和恒定总流动量方程
qv qv
( (
2 v2 x 2V2 y
1v1x ) 1v1y )
Fx Fy
qv
(
2v2
z
1v1z
)
Fz
外力项
不包括惯 性力
输入项 19
恒定总流的动量方程
6 未知力的方向可以假定，若计算为正值，则说明假定正确； 反之，则说明实际力的方向和假定相反。
7 动量方程只能求解一个未知数，如果未知数的数目多于一， 必须联合其他方程（连续方程、或能量程）方可求解。
恒定总流的动量方程
上述积分问题 的解决
用断面平均流速v 代替点流速。定
义V的大小为v ，
方向为u的方向 。
uudA 2 A A
造成的误差用动 量修正系数 来 修正。
uudA
A
2A
2021/4/10
14
恒定总流的动量方程
引入动量修正系数后：
K1-1' dt A1 u1u1dA1 dt1v12 A1 qv11v1dt K 2-2' dt A2 u2u2dA2 dt2 v22 A2 qv22 v2dt
整理得： FR 0qvV0
2021/4/10
26
arctg FRz
FRx
23
恒定总流的动量方程的应用
管轴水平放置 1 1
重力与水流方 向垂直，可忽 略。
FP1=p1A1
V1 FRy FR FRx
V2
2
2
沿x方向列动量方程：
y
x
FP2=p2A·2
沿y方向列动量方程为：
p1A1 FRx qv (0 1V1)
FRx p1A1 1qvV1
K1'2'
流体的连续性方程和动量方程
流体的连续性方程和动量方程流体力学是研究流体运动和流体力学性质的学科。
在流体力学中，连续性方程和动量方程是两个重要的基本方程。
本文将详细介绍流体的连续性方程和动量方程的定义和应用。
一、流体的连续性方程连续性方程描述了流体的质量守恒原理，表达了流体在空间和时间上的连续性。
连续性方程的数学表达形式为：∂ρ/∂t + ∇·(ρv) = 0其中，ρ是流体的密度，t是时间，v是流体的速度矢量，∇·(ρv)表示速度矢量的散度。
该方程表示，流体的密度在一个闭合曲面上的变化率等于通过该曲面的质量流量。
连续性方程是基于质量守恒原理推导得出的。
它表明，在稳定流动条件下，流体在通道中的截面积变化时，速度会发生相应的变化，以保持质量的守恒。
根据连续性方程，我们可以推导出管道中的速度分布。
在管道的收缩段，速度增加，截面积减小，密度保持不变，从而保证质量守恒。
这也是为什么水管收缩后出水流速增加的原因。
二、流体的动量方程动量方程描述了流体运动的力学性质，表达了流体在空间和时间上的动量守恒。
动量方程的数学表达形式为：ρ(dv/dt) = -∇p + μ∇^2v + F其中，ρ是流体的密度，t是时间，v是流体的速度矢量，p是压强，μ是流体的粘度，∇p表示压强的梯度，∇^2v表示速度的拉普拉斯算子，F是外力的合力。
动量方程由牛顿第二定律推导而来。
它表示，在流体中，流体质点的动量变化等于合外力对质点的作用力。
动量方程用于描述流体在受力作用下的运动状态，通过求解动量方程，可以得到流体的速度分布。
根据动量方程，我们可以推导出流体中的压力分布。
在水管中，如果水流速度增大，则根据动量方程中的负梯度项，压力会降低。
这是因为速度增大会导致动能的增加，压力会减少以保持动量守恒。
综上所述，流体的连续性方程和动量方程是流体力学中的两个基本方程。
连续性方程描述了质量守恒原理，动量方程描述了动量守恒原理。
通过求解这两个方程，我们可以获得流体在空间和时间上的运动状态和力学性质。
化工传递(第二章)2014
对流传递 ——由于流体质点的宏观流动引起， 是动量的主体流动过程。
2.1 动量传递概述
1.分子动量传递
分子动量传递的通量由牛顿黏性定律描述： dux d（ ux） yx = dy dy
2.1 动量传递概述
2.涡流动量传递
涡流中大小不等的微团在各流层之间交换。 1877年，波希尼斯克（boussinesq）提出了涡流 通量表达式 d( ux ) r
第二章 动量传递的微分方程
本章先讨论动量传递的基本概念，动量传
递的基本方式：扩散传递和对流动量传递。然
后推导不考虑组分浓度变化的连续性方程和动
量传递的微分方程——运动方程。
2.1 动量传递概述
分子传递 —— 因流场中存在速度 梯度，分子随机运动引 扩散传递 起的动量传递过程。 动量 传递
涡流传递 ——湍流中质点的随机 脉动引起的动量传递。
ux ux
ux dA
第二章 动量传递的变化方程
2.2 连续性方程
一、连续性方程的推导 二、连续性方程的分析和简化 三、柱坐标与球坐标系方程
一、 连续性方程的推导
对单组分流体系统（如水）或组成均匀的多组 分混合物系统（如空气）中，运用质量守恒原理进 行微分质量衡算，所得方程称为连续性方程。 质量守恒定律 流入质量速率—流出质量速率＝积累质量速率 采用欧拉观点 在流场中选一微分控制体。
一、 连续性方程的推导
空间M（x,y,z）点处取微元控制体 dV=dxdydz u x, y, z, 该点流速：
ux
( ux ) x
流体密度： x, y,z, 设流体在M点的质量 通量为 u u 在坐标x,y,z方向分量： 微分质量衡算 ux，uy，uz uz 。 u 沿坐标x,y,z方向分量： ux、 u y、
2.1 动量传递概述
1.分子动量传递
分子动量传递的通量由牛顿黏性定律描述： dux d（ ux） yx = dy dy
2.1 动量传递概述
2.涡流动量传递
涡流中大小不等的微团在各流层之间交换。 1877年，波希尼斯克（boussinesq）提出了涡流 通量表达式 d( ux ) r
第二章 动量传递的微分方程
本章先讨论动量传递的基本概念，动量传
递的基本方式：扩散传递和对流动量传递。然
后推导不考虑组分浓度变化的连续性方程和动
量传递的微分方程——运动方程。
2.1 动量传递概述
分子传递 —— 因流场中存在速度 梯度，分子随机运动引 扩散传递 起的动量传递过程。 动量 传递
涡流传递 ——湍流中质点的随机 脉动引起的动量传递。
ux ux
ux dA
第二章 动量传递的变化方程
2.2 连续性方程
一、连续性方程的推导 二、连续性方程的分析和简化 三、柱坐标与球坐标系方程
一、 连续性方程的推导
对单组分流体系统（如水）或组成均匀的多组 分混合物系统（如空气）中，运用质量守恒原理进 行微分质量衡算，所得方程称为连续性方程。 质量守恒定律 流入质量速率—流出质量速率＝积累质量速率 采用欧拉观点 在流场中选一微分控制体。
一、 连续性方程的推导
空间M（x,y,z）点处取微元控制体 dV=dxdydz u x, y, z, 该点流速：
ux
( ux ) x
流体密度： x, y,z, 设流体在M点的质量 通量为 u u 在坐标x,y,z方向分量： 微分质量衡算 ux，uy，uz uz 。 u 沿坐标x,y,z方向分量： ux、 u y、
《流体力学》流体力学基本方程
2.2 描述流体运动的一些基本概念
2.2.1定常流与非定常流
流场中所有的运动 要素不随时间变化
u u(x, y, z)
(x, y, z)
p p(x, y, z)
u 0 t p 0 t
0
t
流场中有运动 要素随时间变化
u u(x, y, z,t)
(x, y, z,t)
p p(x, y, z,t)
p p(x, y, z,t) (x, y, z,t)
x, y, z ,t--欧拉变量，其中x，y，z与时间t有关。
欧拉法是常用的方法。
5
16 October 2021
欧拉法中的加速度 -- 质点速度矢量对时间的变化率。
a
u t
ux
u x
uy
u y
uz
u z
三个分量：
ax
ux t
ux
ux x
拉格朗日法 从流体质点的运动着手，描述每一个流体质点自始至 终的运动过程。如果知道了所有流体质点的运动规律，那么整个流 体的运动规律也就清楚了。是质点--时间描述法。
质点运动的轨迹
x x(a,b,c,t)
y y(a,b,c,t)
z z(a,b,c,t)
a, b, c --- t = t0 时刻质点所在的空间位置坐标， 称为拉格朗日变量，用来指定质点。
ln x t ln y t ln c
(x t)(y t) c
将 t = 0，x = -1，y = -1 代入，得瞬时流线 xy = 1, 流线是双曲线。
y x
12
16 October 2021
2. 求迹线
将已知速度分布代入式(2.2.1)可得
dx x t, dy ( y t), dz 0
流体力学-第二章 基本方程
h
0
xy
z
经流体柱后侧流入的流体质量应为：
流入质量=
h
0
uy
z
同时，经流体柱前侧流出的质量为：
z
流出质量=
h
0
uy
z
x
h
0
uy
z
x
O
x u u x
x
y
u
h y
x
Chen Haishan NIM NUIST
流出质量减去流入质量 ＝柱体内质量的减少。
柱体内的净流出量
（流入质量减去流出质量 ＝柱体内质量的增加）
pnx nx pxx ny pyx nz pzx
pny nx pxy ny pyy nz pzy
pnz
nx pxz
ny pyz
nz pzz
Chen Haishan
NIM NUIST
z
pzz
z
pzx
pz pzy
pxz
px
pxx
pxy
pyy
pyx
py
P Pnz n
Pny
y Pnx o
Chen Haishan NIM NUIST
通过体积分，作用于体积为 的流体块上的质量力：
Fd =作用于流体的质量力
Chen Haishan NIM NUIST
② 表面力
表面力：是指流体内部之间或者流体与其他物体之 间的接触面上所受到的相互作用力。
如流体内部的粘性应力和压力、流体与固体接触面 上的摩擦力等。
x y
n n
cosn, cosn,
x y
nxn n y n
z n cosn, z nzn
Chen Haishan NIM NUIST
地下水数值模拟02_地下水运动的数学模型
2
H 0
n 2
——隔水边界
第三类边界条件 H aH b n
例：弱透水边界
K H Hn H 0 n m1 / K1
溶质运移问题的边界条件
第一类边界条件
c(x,
y, z,t) 1
c1(x,
y, z,t)
——给定浓度边界
第二类边界条件 c
Di, j x j ni 2 f2 (xi , t)
u(x, y, z,t) t0 0(x, y, z)
• 2、边界条件
第一类边界条件 u(x, y, z,t) 1 1(x, y, z,t)
第二类边界条件
u n
2
1(x, y, z,t)
第三类边界条件
u
u n
3
3x,
y, z,t
水流问题的边界条件
Reynolds数小于1~10
• 有些情况下，用液体压强表示更为方便
– 例如：油水两相流动
vx
K
H x
vy
K
H y
vz
K
H z
K g k
H z p
g
k p
vx
x
v y
k
p y
vz
k
K ( d
)
dhc
C
t
x
K( )
x
y
K
(
)
y
z
K (
动力气象学第2章描写大气运动的基本方程组
Fi
i
→ ：单位质量空气质点受到的真实力
→ 广义牛顿粘性假设，有
→
左边：加速度项； 右边：引起大气运动变化的原因
用
近似表示，a是地球半径 －－万有引力
－－惯性离心力项 万有引力＋惯性离心力＝重力
垂直地面向下
－－压力梯度力
－－科氏力 （地转偏向力）
：分子粘性力
重力： 保守力 科氏力：不做功，只改变运动方向
sin
u)
k
(2
cos
u)
其中: ~f 2 cos , f 104 s1; 7.292 105 s1
u
t
u
u x
Hale Waihona Puke vu yw u z
1
p x
~fw
fv
Fx
v
t
u
v x
v
v y
w
v z
1
p y
fu
Fy
w
t
u
w x
v
w y
w w z
1
p z
~fu
g
Fz
u v w (u v w) 0
（运动形式） 分子粘性力：耗散 驱动大气运动的主要动力：压力梯度力
从以上讨论可见： 物理上－－压力梯度力是驱动大气运动 的主要因子，而压力的变化与热力与动 力过程相关联，因此描写大气过程必须 考虑热力过程。
数学上：运动方程：1个（矢量） 3个（分量）
未知量：速度、气压、密度 必须寻找描写气压、密度变化的方程 －－方程才能闭合
第二章 描写大气运动的基本方程组
一、运动方程： 牛顿第二定律：（单位质量的气团）
成立条件：绝对（惯性）坐标系
【流变力学 精】1.1流变学基础理论
连续介质模型
流体的连续介质模型，应包含两个内容：
• 其一，流体是由排列的流体质点所组成，即空间每 一点都被确定的流体质点所占据，其中并无间隙。 于是流体的任一物理参数B可以表达成空间坐标(x,y,z) 及时间t的连续函数B=B(x,y,z,t)； • 其二，在充满连续介质的空间中，物理参数B不单是 x,y,z的连续函数，而且是连续可微函数。 严格来说，这一模型是一种假定，但这是被实践证明 了的ቤተ መጻሕፍቲ ባይዱ确的假定。
• 对于另一流体质点(a2,b2,c2)在t时刻的坐标为：
•
x2=x2(a2,b2,c2,t)
•
y2=y2(a2,b2,c2,t)
•
z2=z2(a2,b2,c2,t)
拉格朗日法(lagrange)－质点法（3）
• 但标当值研(a，究b任，意c)流不体一质样点，的因位此置各时质，点由在于任各意个时质刻点的在空t间=t位0时置刻，的将坐 是a，b，c，t这四个量的函数。 x=x (a,b,c,t) y=y (a,b,c,t) z=z (a,b,c,t)
目录
• 第三章 1. 两平行平板之间的流变过程 2. 在缝模中的流变过程的分析
第一章
1. 流动运动的描述
流场：流体所占据的空间 连续介质模型：流体所占据的空间
连续介质模型
•连续介质模型 不考虑微观分子结构，把流体视为由 无数多个充满所在空间、相互间无任 何间隙的质点所组成，相邻质点宏观 物理量的变化是连续的。
特征体积ΔV0
图1-1(b)告诉我们，在包含P点的微小体积△V向某一个 △V0逐渐缩小的过程中，其平均密度ρ1,ρ2，……，ρi， ρi+1……，的变化逐渐缓慢，而且当体积△V连续收缩达到 △V0时其平均密度ρ0不再变化，而且到一个恒定的极限值， 这个△V0尺度就是连续介质模型中的质点的尺度。这时， △V0内包括的分子数与随机走出△V0的分子数趋于常数， 即随机进入△V0的分子数与随机走出△V0的分子数趋于平 衡。若将△V体积再进一步缩小，此时△V内的分子数已 减小到不能随机平衡的地步。由于分子数的随机波动，从 而引起△V体积内流体平均密度的随机波动，分子的个性 即随机性与不连续性将显示出来，导致密度发生急剧的变 化。
传递现象基本方程
3、偏导数
偏导数是指在固定位置处某物理量对时间的导数，记为 偏导数又称为局部导数或当地导数。
10:50
t
8
三种导数的物理意义：
偏导数的物理意义：空间中固定位置处观察到的某物理量随时 间的变化率。 全导数的物理意义： 当观察者以任意速度运动时，某物理量随 时间的变化率。 随体导数的物理意义：当观察者随流体一起运动时，某物理量随 时间的变化率。
A u
10:50
式中，rA表示单位体积由于化学反应引起的A的质量变化速率。 这里规定A生成，rA取正值；A消耗，rA取负值
21
将菲克扩散定律的表达式 j A DAB 带入上式，有
D A A u DAB 2 A rA Dt
2、多组分体系连续性方程的简化
r —— 为径向距离 θ——为方位角 z —— 为轴向距离
10:50
18
球坐标下的连续性方程
1 1 1 2 2 ( r ur ) ( u sin ) ( u ) 0 t r r r sin r sin
r —— 为径向距离
θ —— 为余纬度
u u( x0 , y0 , z0 , t )
注意：这里的x0，y0，z0是质点标号，是 t = 0时刻质点的位置， 不是场点坐标。对于不同的流体质点x0，y0，z0有不同的数值。
10:50 3
拉格朗日观点常用于微分动量方程和能量传递微分方程的推 导，推导时研究对象选择的是空间中固定质量的流体微元， 其位置和体积可以变化的。此外，拉格朗日观点也常用于理 论分析当中。
( ux ) d xdy d z x
10:50 12
同理，在 y 方向和 z 方向上的净质量流率分别为
化工传递过程基础(第三版)
※ 湍流流体在流向上的动量，分子传递+涡流传递。
第二节 动量、热量与质量传递的类似性
1. 分子间动量传递
※ 牛顿粘性定律
dux
dy
2. 分子间热量传递 —— 热传导
※ 傅立叶定律
q k dt
A
dy
高温
低温
3. 分子间质量传递 ——分子扩散
※ 费克定律
jA
DAB
d A
dy
一、分子传递的基本定律
图0-1 McCabe-Thiele图
二、本课程的学习内容？
✓ 物理过程的速率和传递机理的探讨
• 动量传递
• 热量传递
• 质量传递
推动力：速度差 推动力：温度差 推动力：浓度差
第一章 传递过程概论
第一节 流体流动导论
※ 流体：气体和液体的统称
一、静止流体的特性 （一）流体的密度（ρ）
均质流体：
※ 非均质流体： f x, y, z
dx
ux d
dy
uy d
dz
uz d
流率：单位时间内流体通过流动截面的量
[m/s]
※ 以流体的体积计量称为体积流率（流量，Vs）m3/s ※ 以质量计量称为质量流率（w），kg/s
计算：在流动截面上任取一微分面积dA，其点流速为ux，则通过该微元面积 的体积流率dVs？通过整个流动截面积A的体积流率Vs？
p0
0
流体静力学方程
p p0 gh
h p p0
g
※对于一定密度的液体，压力差与深度h成正比，故 液柱高度h可用来表示压力差的大小（mmHg,mH2O)
二、流体流动的基本概念
（一）流速与流率
流速：流体流动的速度，表示为 u u f (x, y, z, )
第二节 动量、热量与质量传递的类似性
1. 分子间动量传递
※ 牛顿粘性定律
dux
dy
2. 分子间热量传递 —— 热传导
※ 傅立叶定律
q k dt
A
dy
高温
低温
3. 分子间质量传递 ——分子扩散
※ 费克定律
jA
DAB
d A
dy
一、分子传递的基本定律
图0-1 McCabe-Thiele图
二、本课程的学习内容？
✓ 物理过程的速率和传递机理的探讨
• 动量传递
• 热量传递
• 质量传递
推动力：速度差 推动力：温度差 推动力：浓度差
第一章 传递过程概论
第一节 流体流动导论
※ 流体：气体和液体的统称
一、静止流体的特性 （一）流体的密度（ρ）
均质流体：
※ 非均质流体： f x, y, z
dx
ux d
dy
uy d
dz
uz d
流率：单位时间内流体通过流动截面的量
[m/s]
※ 以流体的体积计量称为体积流率（流量，Vs）m3/s ※ 以质量计量称为质量流率（w），kg/s
计算：在流动截面上任取一微分面积dA，其点流速为ux，则通过该微元面积 的体积流率dVs？通过整个流动截面积A的体积流率Vs？
p0
0
流体静力学方程
p p0 gh
h p p0
g
※对于一定密度的液体，压力差与深度h成正比，故 液柱高度h可用来表示压力差的大小（mmHg,mH2O)
二、流体流动的基本概念
（一）流速与流率
流速：流体流动的速度，表示为 u u f (x, y, z, )
第2章 流体运动方程组
第二章 流体运动方程组§1.连续方程§2.作用于流体上的力§3.流体运动方程及其简化形式 §4.能量方程§5.Naver-Stokes 方程的简单解本章重点：流体的三大守恒定律，作用在流体上的力，Naver-Stokes 方程。
作为物体的形态之一，流体也遵循基本的物理规律：质量守恒定律、动量守恒定律和能量守恒定律，其分别对应本章的连续方程、运动方程和能量方程。
§1.连续方程1. 连续方程 设流体块体积δτδδδ=x y z ，则质量δρδτm =。
由于质量守恒，有：()dd δm =0t （2.1）()d d ρδτ=0t（2.1´）’展开：()d d d d δτρδτρ+=0t t（2.2） 同除δτ，得：()d d d d δτρρδτ+=0t t（2.3）∵()d d δτδτ∇⋅ =1V t（体胀速度）∴（2.3）式可变为:——连续方程（速度散度形式） （2.4）∵d d ρρρ∂⋅∇∂+=V t t，而()ρρ⋅∇∇⋅∇⋅ +ρ V V =V ，则（2.4）式可改写为:——连续方程（质量散度形式） （2.5） 其中称为（速度）散度，表示单位体积的流体通量。
而∇⋅V ()ρ∇⋅V 称为质量散度，表示单位体积的流体质量通量。
质量有净流入,→()ρ∇⋅V <0ρ∂∂>0t（密度增大）；质量有净流出,()ρ∇⋅ V >0→ρ∂∂<0t（密度减小）。
2. 有关流体密度的几种近似 1） 若d d ρ=0t，即const ρ()x,y,z,t =，称为不可压(缩）流体。
∵d d ρ=0t →∇⋅ V =0，∴不可压流体＝（三维）无辐散流体。
2）若d d ρ=0t，且各处的ρ（常数）也一样， const ρ=（不随而异），称为均质（均匀）不可压（缩）流体。
x,y,z,t 3）若ρ∂∂=0t ，则ρ与无关，称为（密度）定常流体（不同于定常流场）。
第二章 第1节连续方程
可见，当
）或不定常性（
0 t
）。
（流体质点在运动中保持密度不变）时，流体场内
的流体密度可以不定常和不均匀，只要这两者的变化对流体块 的贡献正好抵消就可以。那么就有了以下的特例：
* 0 ：均质流体。 流体密度空间处处相等，但可能随时间在变 化，只不过，不管怎么变化，密度处处相同，或说密度处处发生 着相同的变化。 * ，或 不可压缩流体。某一流体质点在运动中
第二章 基本方程
首先回忆一下描写流体运动的两种观点： 拉格朗日观点和欧拉观点
1
从这一章开始介绍流体动力学。
首先将根据力学中的普遍规律——质量守恒、动量守恒、 能量守恒定律推导出流体力学基本方程------连续方程、运 动方程、能量方程。
从现在开始，都使用欧拉观点描述流体运动。
2
§1 连续方程 （由质量守恒定律推出的） 一、连续方程的推导： （1）拉氏观点的连续方程 设在流体中取一块由若干流点组成的小流体团，在这个流体团 运动中，不管流体团发生何种形变，总是由那些流点组成， 则这个流体团的质量始终不变。 这个流体团的体积是 密度是 质量是 （2.1）
该点密度必然减少 ，反之。可见，流动改变了流体质量的分布。
6二Βιβλιοθήκη 讨论的流体称为无辐散流体，也称为不可压缩流体。 此时根据（2.2） 即： 拉氏观点 表示流体质点在运动中保持密度不变。但这时 不一定为零 可知： （2.8） 欧拉观点 和
，表示了密度的不均匀性。（若
表示密度在空间处处都一样）。
7
表示了流体的定常（
4
（2）欧拉观点的连续方程
上一章我们已推出：
则：
带入（2.2）,并合并就得到： （2.6） 欧拉观点的连续方程
这是欧拉观点的连续方程，也是我们常用的形式。
）或不定常性（
0 t
）。
（流体质点在运动中保持密度不变）时，流体场内
的流体密度可以不定常和不均匀，只要这两者的变化对流体块 的贡献正好抵消就可以。那么就有了以下的特例：
* 0 ：均质流体。 流体密度空间处处相等，但可能随时间在变 化，只不过，不管怎么变化，密度处处相同，或说密度处处发生 着相同的变化。 * ，或 不可压缩流体。某一流体质点在运动中
第二章 基本方程
首先回忆一下描写流体运动的两种观点： 拉格朗日观点和欧拉观点
1
从这一章开始介绍流体动力学。
首先将根据力学中的普遍规律——质量守恒、动量守恒、 能量守恒定律推导出流体力学基本方程------连续方程、运 动方程、能量方程。
从现在开始，都使用欧拉观点描述流体运动。
2
§1 连续方程 （由质量守恒定律推出的） 一、连续方程的推导： （1）拉氏观点的连续方程 设在流体中取一块由若干流点组成的小流体团，在这个流体团 运动中，不管流体团发生何种形变，总是由那些流点组成， 则这个流体团的质量始终不变。 这个流体团的体积是 密度是 质量是 （2.1）
该点密度必然减少 ，反之。可见，流动改变了流体质量的分布。
6二Βιβλιοθήκη 讨论的流体称为无辐散流体，也称为不可压缩流体。 此时根据（2.2） 即： 拉氏观点 表示流体质点在运动中保持密度不变。但这时 不一定为零 可知： （2.8） 欧拉观点 和
，表示了密度的不均匀性。（若
表示密度在空间处处都一样）。
7
表示了流体的定常（
4
（2）欧拉观点的连续方程
上一章我们已推出：
则：
带入（2.2）,并合并就得到： （2.6） 欧拉观点的连续方程
这是欧拉观点的连续方程，也是我们常用的形式。
传递过程基础总结
cp k
Pr 同时存在动量、热量传递 。
DAB
DAB
k c p DAB
Sc 同时存在动量、质量传递 。 Le 同时存在热量、质量传递 。
DAB
若三个数均等于 1，则表示同时进行的两种传递过程可以类比。 3、传递过程、分子传递和涡流传递概念。 传递过程——质量、能量、动量等具有强度性质的物理量可由高强度向低强
化工传递过程基础总结
化研 1205 班
宁鹏
4、势函数的定义式、势函数存在的判据。 ①定义：对于不可压缩流体的平面二维流动，若存在速度势 ( x, y ) ，且满足
u x u y x y
，则 ( x, y ) 称为势函数。
②存在的判据：理想流体做无旋运动，或有势运动时，势函数存在判断旋度 u u x y 。 为 0 的方法：二维 y x
因为 y 0时，u x umax ，所以 umax
y 2 从而得出： u x umax 1 y 0
1 p 2 y0 2 x
第 4 页 共 28 页
化工传递过程基础总结
化研 1205 班
宁鹏
若在 x 方向取单位宽度的流通截面 A 2 y0 1 ，则通过该界面的体积流率 Vs 为： Vs u x dy 2 u x dy
1、什么是欧拉研究方法？ 在流场内某一固定位置， 找一固定体积的流体微元，但该微元的质量可随时 间改变， 观察者分析该流体微元的流动状态，并由此获得整个流场流体运动的规 律。 特点：流体微元的位置和体积不随时间变化，而质量随时间变化。 2、什么是拉格朗日研究方法？ 在流场内选择一固定质量的流体微元，观察者追随流体微元一起运动，并研 究其运动规律，据此获得整个流场内流体的运动规律。 特点：流体微元的质量不随时间变化，而而位置和体积随时间改变。 3、随体导数、全导数、偏导数的定义式和物理意义。 以流体密度ρ为例： 定义式： 偏导数： 全导数：
第二章 连续性方程与运动方程
二、连续性方程的简化
1. 稳态流动
( ρux ) ( ρu y ) ( ρuz ) 0 x y z
2. 不可压缩流体
ux u y uz 0 x y z
u 0
三、柱坐标与球坐标系方程
1. 柱坐标系
ρ 1 1 ( ρrur ) ( ρuθ ) ( ρuz ) 0 ' θ r r r θ z
第二节 连续性方程
一、连续性方程的推导 二、连续性方程的简化 三、柱坐标与球坐标系方程
一、 连续性方程的推导
于单组分流体系统（如水）或组成均匀的多组分 混合物系统（如空气）中，运用质量守恒原理进行 微分质量衡算，所得方程称为连续性方程。 质量守恒定律 流出质量速率－流入质量速率+ 积累质量速率＝0 采用欧拉观点 在流场中选一微分控制体。
D d D d f ( , x, y, z )
d dx dy dz d x d y d z d
风洞实验
车辆行驶的风洞(Wind tunnel)试验
随体导数
物理意义： 流场中流体质点上的物理量（如温度）随时 间和空间的变化率。
第二章 连续性方程与运动方程
θ －时间； r －径向坐标； －方位角； θ－余纬度； ur , uφ , uθ －各方向的 速度分量。
球 坐 标 系
例 题
例.对于在 r 平面内的不可压缩流体的流动，r 方向的速度分量为
ur A cos / r
试确定
2
方向的速度分量 u 的表达式。
第二章 动量传递的变化方程
第一篇 动量传递
第二章 连续性方程与运动方程 第一节 描述流动问题的两种观点 第二节 连续性方程 第三节 运动方程 第三章 运动方程的应用
第二章 连续方程和运动方程
第二节 连续性方程（微分质量方程） 连续性方程（微分质量方程）
一、连续性方程的推导： 连续性方程的推导： 研究方法： 研究方法：欧拉观点 理论依据： 理论依据：质量守恒定律 计算依据：输出-输入+累积=0 计算依据：输出-输入+累积=0 (*) ∂ ( ρu x ) x方向：输出 - 输入 = 方向： d xd yd z (1) ∂x ∂ ( ρu y ) y方向：输出 - 输入 = 方向： d xd yd z (2) ∂y
拉格朗日法
着眼于流体质点， 着眼于流体质点，跟踪 质点描述其运动历程
欧拉法
是描述液体运动 常用的一种方法。 常用的一种方法。
着眼于空间点，研究 质点流经空间各固定 点的运动特性
二、物理量的时间倒数：偏导数、全导数和随体导数(真实导数) 物理量的时间倒数：偏导数、全导数和随体导数(真实导数) 1. 三者定义(以流体密度ρ为例) 三者定义(以流体密度ρ为例) ∂ρ ∂ρ ∂ρ ∂ρ dθ + dx + dy + dz y z θ 若 ρ = f (x， ， ， ) ⇒ d ρ = ∂θ ∂x ∂y ∂z
∂(ρux ) ∂(ρuy ) ∂(ρuz ) ∂ρ (2-1) + + + =0 ∂x ∂y ∂z ∂θ ∂u x ∂u y ∂uz ∂ρ ∂ρ ∂ρ ∂ρ ) + ux + + uy + uz + ⇒ ρ( + =0 ∂y ∂z ∂x ∂x ∂y ∂z ∂θ r Dρ ∂u x ∂ u y ∂uz Dρ )+ (2-2) ⇒ ρ( + + = 0 或 ρ∇ ⋅ u + =0 Dθ Dθ ∂x ∂y ∂z
∂ ( ρu z ) z方向：输出 - 输入 = 方向： d x d yd z ∂z (3) z
连续性方程和运动方程的公式运用
– U—流体的速度向量；θ—时间。
• 惯性力＝外力＝（质量）＊加速度
42
• 拉格朗日法：在流体运动的空间内，选择 某一固定质量的流体微元（M为常数），观 察者追随此流体微元且一起运动（在相对 坐标系下，可以用随体导数的概念来描 述），并根据此运动流体微元的变化状况 来研究整个流场中流体运动规律。
• 固定质量的流体微元： 体积 dv dxdydz
测者在流体中与流体流速完全相同的速度运动。
此时：
ux
dx
d
;uy
dy
d
;uz
dz
d
24
随体导数
d dx dy dz d x d y d z d
d d
ux
x
uy y uz z
D D
ux
x
uy
y
uz
z
对温度t、浓度c等也有类似表达式
Dt
D
t
ux
t x
uy
t y
uz
• 特例：在速度为零或无限大的空间点上例外，速 度为零的点称为驻点，速度无限大的点称为奇点。
• ③流线的形状和位置，在稳态流动中不随时间变 化，在非稳态流动中，一般要随时间变化。
• ④对于稳态流动，流场中任何参数均不随时间变 化，故流线方程与轨线方程重合。
16
例题; 已知流体运动速度为
ux x uy y uz 0
z (x,y,z)或
θ （r,Φ,θ）
y θ
y Φ
x
x
38
柱坐标和球坐标连续性方程式
• 柱坐标： • 球坐标：
1 r
(
ru
r
r
)
1 r
(u
)
流变学的基础方程-北京化工大学
或者说：净流出与质量变化之和为零 方程性质：这是一个标量方程
三、特殊形式
1、稳定流场（ ∂ρ = 0 ） ∂t
( ) ∇ • ( ρ v ) = 0 或 ∂ ( ρ v x ) + ∂ ρ v y + ∂ ( ρ v z ) = 0 净流出之和为零
∂x
∂y
∂z
2、不可压缩流体（ ∂ρ ≡ 0 ， ρ 为常数） ∂t
⎛ ⎜ ⎝
∂vz ∂t
+
vx
∂vz ∂x
+
vy
∂vz ∂y
+
vz
∂vz ∂z
⎞ ⎟ ⎠
=
−
∂P ∂z
+
∂τ xz ∂x
+
∂τ yz ∂y
+
∂τ zz ∂z
+
ρ gz
矢量微分表示式：
ρ Dv = −∇ P + ∇ • [τ ] + υ g
Dt
二、物理意义
流体体系在 x、y、z 方向上单位体积动量的变化等于流体体系在 x、y、z 方向上压 力变化、应力变化和重力之和。或者：
沿 x－x’面净流出流体质量： ∂ ρ v x dxdydz ∂x
沿 y－y’面净流出流体质量： ∂ ρ v y dxdydz ∂y
沿 z－z’面净流出流体质量： ∂ ρ v z dxdydz ∂z
于是：
( ) ∂ρ dxdydz ∂t
=
−
⎡ ⎢ ⎣
∂
(ρ v
∂x
x
)
+
∂
ρv y ∂y
+
∂ (ρv z
1、掌握的内容
⑴ 连续性方程的数学表示式（矢量微分式子及其在直角坐标系下的展开式）、 物理意义；不可压缩流体连续性方程的数学表示式；
三、特殊形式
1、稳定流场（ ∂ρ = 0 ） ∂t
( ) ∇ • ( ρ v ) = 0 或 ∂ ( ρ v x ) + ∂ ρ v y + ∂ ( ρ v z ) = 0 净流出之和为零
∂x
∂y
∂z
2、不可压缩流体（ ∂ρ ≡ 0 ， ρ 为常数） ∂t
⎛ ⎜ ⎝
∂vz ∂t
+
vx
∂vz ∂x
+
vy
∂vz ∂y
+
vz
∂vz ∂z
⎞ ⎟ ⎠
=
−
∂P ∂z
+
∂τ xz ∂x
+
∂τ yz ∂y
+
∂τ zz ∂z
+
ρ gz
矢量微分表示式：
ρ Dv = −∇ P + ∇ • [τ ] + υ g
Dt
二、物理意义
流体体系在 x、y、z 方向上单位体积动量的变化等于流体体系在 x、y、z 方向上压 力变化、应力变化和重力之和。或者：
沿 x－x’面净流出流体质量： ∂ ρ v x dxdydz ∂x
沿 y－y’面净流出流体质量： ∂ ρ v y dxdydz ∂y
沿 z－z’面净流出流体质量： ∂ ρ v z dxdydz ∂z
于是：
( ) ∂ρ dxdydz ∂t
=
−
⎡ ⎢ ⎣
∂
(ρ v
∂x
x
)
+
∂
ρv y ∂y
+
∂ (ρv z
1、掌握的内容
⑴ 连续性方程的数学表示式（矢量微分式子及其在直角坐标系下的展开式）、 物理意义；不可压缩流体连续性方程的数学表示式；
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1 r
rur
r
1 r
u
uz
z
'
0
同理，可得球坐标连续性方程
1 r2
(rur )
r
1
r sin
(u sin )
1
r sin
(
u )
0
运动方程
通过微分动量衡算，可以导出流体的运 动方程。运动方程与连续性方程结合起 来，可以处理许多流体流动问题。同时 运动方程在动量、热量与质量传递过程 中也是求解大量有实际意义问题的基础 方程。本节在推导运动方程时采用拉格 朗日观点。
连续性方程的推导
单组份系统：
y
(x,y,z)
（输出的质量流率）—（输入的质量流率）
＋累积的质量速率＝0
dy
在x左侧面：
输入微元体积的质量流率 uxdydz
z
输出微元体积的质量流率
dy
ρux
(ux
(ux
x
)
dx)dydz
dz
dz dx
x
ux
(ux
x
)
dx
dy
dz
连续性方程的推导
θ （r,Φ,θ）
y θ
y Φ
x
x
公式回顾：
sV s V s V
s
1 h1
s q1
,
1 h2
s q2
,
1 h3
s q3
V
1 h1h2h3
V1h2h3
q1
V2h3h1
q2
V3h2h1
) ux
x
uy
y
uz
z
0
由于密度ρ是空间（x，y，z）和时间θ的连续函数，
即ρ＝f(x,y,z,θ），那么ρ的随体导数
D D
x
ux
y
uy
z
uz
于是 D ( (ux ) (uy ) (uz ) ) 即
面上的法向应力与切线 应力都是相应地大小相
yx(下）
dx
dz
zx
zx z
dz(后）
等、方向相反的。故只
x
需采用9个机械应力就可 以完全表达：3个法向分
z
量，6个切线分量。
现将上图中的流体微元在x—y
平面的一个相应的平面分离出
来加以考察。环绕该平面四周 所作用的4个剪应力，可表示 在右图中。由图可见，假如有
q3
连续性方程式为：
u 0
u u u
r
1 1
rቤተ መጻሕፍቲ ባይዱ
1 r
z
1 1
z
r
r
1 r
z
z
u
rur
u
x
y
z
(ux ) (uy ) (uz ) 0
x
y
z
写成向量形式 (u) 0
连续性方程的进一步分析
将连续性方程 (u) 0 展开，得到
( ux
x
u y y
uz z
用应力表示的运动方程
任何物体的运动，都遵循动量守恒定律即牛顿第二定律，流体
的运动也不例外。将牛顿第二定律应用于运动着的流体时，可
理解为：流体的动量随时间的变化率应等于作用在该流体上的
诸外力向量之和，即
F
d (Mu)
d
式中 F——诸外力向量之和
M——流体的质量；
u——流体的速度向量
θ——时间。
v D D
D D
u
1 Dv u
v D
上式的左侧
表示流体微元的体积膨胀速率或形变速率，右
侧是速度向量的散度
几种特殊情况下连续方程简化
稳态流动，密度不随时间变化，即
0
( ux ) ( u y ) ( uz )
(
) 0
体积力
令fB表示单位质量流体所受的质量力，其在直角坐
标x，y，z方向上的分量分别为X，Y和Z，则
根据上述定义，可知所考察的流体微元上所受的质
量力为
dFB dFBxi dFBy j dFBzk fBdxdydz
写成坐标分量形式为
表面力
y
单位面积上的表面力定义为表面应力
或机械应力，表面应力亦可分解为法
u r
uz
z
1
r
ur r r
ur
r
r
1 r
u
u
uz z
uz
z
1 rur 1 u uz r r r z
于是柱坐标连续性方程为
u z r z
1 urr 1 u uz
r r r z
u urr u uz r r r z
u u u
r
ur r r
r
u
uz z
ur
r
到另一点时所发生的变化。； 的物理意义为：
当流体质点在dθ时间内由空间的一点（x, y, z)
移动到另一点（x+dx, y+dy, z+dz)时，流体密度ρ
随时间的变化率。
由于 v 1
将上式对时间求随体导数，亦即
Dv 0 D
Dv v D 0 D D
1 Dv 1 D 0
zuz
r
r
1 r
z
z
ur
r
u r
uz
z
u
1
1 r
1
ur
r r
1
u 11
uz
1 z
r
1 urr r r
u
y
yx
yx y
dy 2
一根平行于z轴的轴线或z轴本
身穿过该流体微元的形心O点 时，显然，出于上述这四个剪
xy
xy x
D
x y z
为连续性方程的另一表达形式。
D u D
随体导数 的意义
D D
ux
x
uy
y
uz
z
局部导数： 间的变化；
表示ρ在空间的一个固定点处随时
对流导数：( x
ux
y
uy
z
uz )
表示密度由一点移动
即
x
y
z
上式可简化为：(ux ) (u y ) (uz ) 0
x
y
z
u 0
对于不可压缩流体，ρ＝常数，则无论稳态还是非稳态：
ux u y uz 0 x y z
• u 0
上式为不可压缩流体的连续性方程。
全导数的表达式可由对t进行全微分得到
随体导数：观测者随流体随波逐流运动，即观测者在流体
中与流体流速完全相同的速度运动。第三种测量大气温
度的方法是将测温计装于探空气球上。此时探空气球随
空气一起漂动，其速度与周围大气的速度相同。观察者 记录下不同时刻的大气温度。如此获得的温度t随时间θ 的变化称为随体导数(Substantial derivatives)亦称拉格 朗日导数(Largrangian derivatives)，以Dt/Dθ表示。
小微元流体在运动时，由于法向应力
和剪应力的存在，使其发生形变。
τxx
x
六个表面，每一表面的
机械应力均可分解成三
个平行于x、y、z三个坐
y
标轴的应力分量3×6=18
yx
yx y
dy(上）
个在x、y、z方向上各有
六个。当小微元体体积
zx (前）
xx
xx x
dx(右）
缩小为一点时，相对表 xx(左）dy
向应力和剪应力，一般记为τ。
τxy
τxy 第一个下标表示应力分量的作用面
与x轴垂直。第二个下标x、y、z表示
应力方向为 x轴、y轴和z轴方向。 τxx 表示法向 应力分量。拉伸方向（向外）
τxz
为正，压缩方向（向内）为负。这三
个表面应力分量中．一个是法向应力
z
分量τxx另外两个是剪应力分量τxy和τxz。
由于采用拉格朗日观点，故在推导微分动且衡算方程时，可在流 场中选一固定质量的流体微元即微元系统，如图2－3所示，考察 该微元系统随环境流体一起流动过程中的动量变化。
设在某一时刻θ，此微元系统的体积为dv＝dxdydz(注意其体积
和位置是随时间改变的)，将牛顿第二定律应用于此微元系统得