粘性流体力学第三章.

流动的解为：
2 b dp y u 1 2 dx b
（3－7）
可 以 看 出 ： 有 压 梯 度 的 Couette 流 动 是 简 单 Couette流动和Poiseuille流动的叠加。
10
(2)充分发展的管流－圆管中的Poiseuille流动
h2 dp 令 B 2U dx
（3－6）
式中:
y y* h
u u* U
图3.2 两平行直壁之间的库埃特流动
7
B值不同，流动曲线不同 （1）B 0,
dp 0 顺流压力梯度为零时： dx u y U h
dp
流速为线性分布称为简单的Couette流动。
（2）当B>0， 0 ，压力顺流递减称为顺压梯度，在 dx 整个断面上流速为正值，当 B 值很大时，流动接近 Poiseuille流动的抛物线分布。
u y y y （3）当B＝－1时： U h h (1 h ) u * y 令 u , y* 则 * 2 * u u * *2 * U h u y 2y 2 * 2 * y y
u 从连续方程可以得出 0 ，因此对于平行流动(二 x
2
u u ( y , z , t ) v0 w0
（3－1）
p p 0, 0 ，压强p为P（x） 利用N－S方程可以得到 y z
2u 2u u p 2 2 t x y z
精确解本质上是层流解。从方程上看精确解尽管 在高雷诺数下其数学关系是正确的，但是在高雷诺数 时流体运动不稳定，在物理上数学解不存在。 精确解虽然简单，数量少，但却有重要的理论和 实践意义： 揭示粘性流动的一些本质特征； 应用于发展新的数值计算方法； 作为研究复杂问题初步估算和求解的基础； 探求新理论。
管道很长时，除了进口段，可以认为管流为二维 流动，采用圆柱坐标 (r , , z ) 系，连续方程为： u r r u u x 0 rr r x
其中， ur u 均 为0。只有 u x 不 为零，令 u = u x u 可以看出 0， x 即流速分布沿管的 轴线x是相同的。
p 2u 2 y x p 0 y
（3－3）
4
y
dp / dx 0
U
h
x
图3－1 平行平板间的流动
5
p p 只是y的函数，那么 x x
由于，p只是x的函数；又由于u只是y的函数，故
= 常数。
（3－4）
2u 2 C y
边界条件为：
（
8
Leabharlann Baidu
可见 u* f ( y*)曲线为凹曲线，在 y* 0 时，曲线与 dp 2U y* 轴相切。 2 时为流动要产生回流的临界状态。
dx h
(4) 在 B 1 ， 流动在靠近下壁为负值 有回流出现。这就是说明由于流体的带动上壁的运动 速度传到下壁附近时，不足以抵抗逆压梯度的作用， 而产生反向回流。
图3.3 圆管中泊肃叶流动
11
N－S方程
1 p 0 r 1 1 p 0 r u 1 dp 2u 1 u u 0 2 dx r r r x
（3－8）
dp 由于 只能是常数 dx
式（3－8）为：
第三章 层流流动的精确解
第一节 第二节 第三节 第四节 第五节 平行流动 驻点附近的平面运动 旋转盘引起的流动 缓慢流动的N－S方程的近似解 滑动轴承内的流动
由于N－S方程的非线性，一般情况下在数学上寻 求其精确解有巨大的困难。大多数实际问题要引入不 同程度的物理或数学上的近似求近似解。随着计算机 的发展，数值求解越来越重要。
1
第一节 平行流动
粘性流动的动量方程应包括粘性项，是二阶偏微 分方程，应采用物体表面上流速为零的边界条件。 平行流动是流动中最简单的一种。平行流动中， 所有的质点均沿同一方向流动，即只有一个速度分量 不等于零，令其为 x 方向，即 u≠0 ，而另外两个 y ， z 方向上速度分量v，w 均为零。 阶线性偏微分方程 )
（3－2）
式（3－2）为二阶线性偏微分方程。
3
1、库埃特(Couette)流动
两个平行壁面间的平行流动，一个壁面静止不动， 另一个壁面以速度U沿x轴运动（图3－1）。由于粘性， 运动壁面将带动流体运动。通过流体的内摩擦，这个运 动的影响传播到整个流动区域。设上下两个壁面的宽度 为无穷远，流动为二维定常平行流动，因而 u f ( y) ， 方程（3－2）将有以下形式
dp 2U 2 dx h
2、泊肃叶(Poiseuille)流动
(1) 平面Poiseuille流动
9
在两个平行平板之间充满粘性流体，上下两板均 dp 静止不动，而顺压梯度 dx const ，坐标系仍如图3－1 所示。方程仍如（3－3）式，边界条件为：
u 0 u 0
2
y b y b
u max
断面平均速度
r02 dp 4 dx
（3－11） （3－12）
积分时， 代入边界条件：
d 2 u du 1 dp r 2 r dr dx dr
（3－9）
du r0 dr u 0 r r0
12
圆管中Poiseuille流动的速度分布：
1 dp 2 u r r02 4 dx
圆管中心处最大流速


（3－10）
u 0 u U
y0 yh
（3－4）式积分并代入边界条件则得：
u y h2 dp y y 1 U h 2U dx h h
（3－5）
6
为量纲为1的压力梯度称为Brinkman 数。解(3-5)的量纲为1的形式为：
u y y y u* B 1 y * By * 1 y * U h h h