流体力学公式
第二章 流体的主要物理性质
流体的可压缩性计算、牛顿内摩擦定律的计算、粘度的三种表示方法。
1．密度 ρ = m /V 2．重度 γ = G /V
3．流体的密度和重度有以下的关系：γ = ρ g 或 ρ = γ/ g 4．密度的倒数称为比体积，以υ表示υ = 1/ ρ = V/m 5．流体的相对密度：d = γ流 /γ水 = ρ流 /ρ水 6．热膨胀性
7．压缩性. 体积压缩率κ
8．体积模量
9．流体层接触面上的内摩擦力 T
V V ??=
1αp
V V ??-
=1κV
P V K ??-
=κ
1
n
A F d d υ
μ=
10．单位面积上的内摩擦力（切应力）（牛顿内摩擦定律）
11．.动力粘度μ：
12．运动粘度ν ：ν = μ/ρ
13．恩氏粘度°E ：°E = t 1 / t 2
第三章 流体静力学
重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。
1．常见的质量力： 重力ΔW = Δmg 、
直线运动惯性力ΔFI = Δm·a 离心惯性力ΔFR = Δm·r ω2 . 2．质量力为F 。：F = m ·am = m (f xi+f yj+f zk)
dn
d v
μ
τ±=n
v d /d τ
μ=
am = F /m = f xi+f yj+f zk 为单位质量力，在数值上就等于加速度 实例：重力场中的流体只受到地球引力的作用，取z 轴铅垂向上，xoy 为水平面，则单位质量力在x 、y 、 z 轴上的分量为 fx = 0 , fy = 0 , fz = -mg /m = -g 式中负号表示重力加速度g 与坐标轴z 方向相反
3流体静压强不是矢量，而是标量，仅是坐标的连续函数。即：p = p (x ,y ,z )，由此得静压强的全微分为:
4．欧拉平衡微分方程式
单位质量流体的力平衡方程为：
z
z p
y y p x x
p
p d d d d ??????++=d d d d d d 0
x p f x y z x y z x
??-=ρd d d d d d 0y p
f x y z x y z y
??-=ρd d d d d d 0z p
f x y z x y z z
??-
=ρ0
1=??-x
p f x ρ
5．压强差公式（欧拉平衡微分方程式综合形式）
6．质量力的势函数
7．重力场中平衡流体的质量力势函数
积分得：U = -gz + c
8．等压 .面微分方程式 .fx d x + fy d y + fz d z = 0
10y p
f y
??-
=ρ01=??-
z
p
f z ρz
z
p y y
p x x
p z f y f x f z y x d d d )d d d (??+??+??=++ρ)
d d d (d z f y f x f p z y x ++=ρd (d d d )x y z p f x f y f z dU ρ=++=ρd d d d x y z U U U U x y z =f dx f dy f dz
x
y
z
gdz
??????=++++=-
9．流体静力学基本方程
对于不可压缩流体，ρ = 常数。积分得： p + ρgz = c
形式一
形式二
形式三
10．压强基本公式p = p 0+ρ g h 11．.静压强的计量单位
应力单位：Pa 、N/m2、bar 液柱高单位：mH2O 、mmHg
标准大气压：1 atm = 760 mmHg =10.33 mH2O = 101325 Pa ≈ 1bar
第四章 流体运动学基础
1拉格朗日法：流体质点的运动速度的拉格朗日描述为
1
2
12p p c
+=
+=gz gz ρ
ρ
121
2p p
c g g
+
=+=z z ρρ??
?
??===),,,(),,,(),,,(t c b a w w t c b a t c b a u u υυ
压强 p 的拉格朗日描述是：p =p (a ,b ,c ,t ) 2．欧拉法
流速场
压强场：p =p (x,y,z
,t)
加速度场
简写为
??
?
??===),,,(),,,()
,,,(t z y x w w t z y x t z y x u u υυv ui v j wk
=++ (,,,)x y z a a x y z t a i a j a k ==++ d d (,,)d d d d (,,)d d d d (,,)d d x y
z u u x y z,t u u u u
a u w t t t x y z x y z,t a u w t t t x y z w w x y z,t w w w w
a u w t t t x y z ????????υυ?υ?υ?υ?υυ????????υ?????===+++??
?===+++???===+++??
υυ
υυ
)(??+??=t
a
时变加速度：
位变加速度
3．流线微分方程：.在流线任意一点处取微小线段d l = d xi + d yj + d zk ，
该点速度为：v = ui + vj + wk ，由于v 与d l 方向一致，所以有： d l × v = 0 4．流量计算：
单位时间内通过d A 的微小流量为 d q v=u d A 通过整个过流断面流量
相应的质量流量为
5．平均流速
6．连续性方程的基本形式
t
??υυυ)(????==A
A
u q q d d v v ?==A m A
u q q d v ρρv d A
v u A q A A
q vA
υ==
=?
2
1
2211d d d A A V u A u A V
t
?ρρ?-=?
??ρ
对于定常流动 有 即
ρ1A 1υ1= ρ2A 2υ2
对于不可压缩流体，ρ1 = ρ2 =c ，有
即A 1υ1=A 2υ2= q v
7．三元流动连续性方程式
定常流动
不可压缩流体定常或非定常流：ρ = c
8．雷诺数
0t
??=ρA u A u A A d d 2
1
2211??ρ=ρA
u A u A A d d 2
1
21??=()()()0u w t x y z
ρρρυρ????+++=????()()()
0u w x
y z
ρρυρ???+
+=???0u w
x y z
υ???++=???ud
Re ρμ
=
对于圆管内的流动：
Re <2000 时，流动总是层流型态，称为层流区； Re >4000时，一般出现湍流型态，称为湍流区；
2000 10．剪切应力，或称内摩擦力， N/m2 11．动力黏性系数 12．运动黏度 m2/s F U A y μ=x d d u y τμ =-x d d u y τμ=- ,μνρ = 13．.临界雷诺数 14．进口段长度 第五章 流体动力学基础 1.欧拉运动微分方程式 2.欧拉平衡微分方程式 c 0 xc x u Re ρμ =e l d 1y p dv f y dt ρ?- =?1x p du f x dt ρ?- =?1z p dw f z dt ρ?-=?10y p f y ρ?-=?10x p f x ρ?-=?10z p f z ρ?-=? 3.理想流体的运动微分方程式 4. 理想不可压缩流体重力作用下沿流线的伯努利方程式：三个式子 5．理想流体总流的伯努利方程式 6．总流的伯努利方程 1x p u u u u f u w x t x y z υρ?????-=+++?????1y p f u w y t x y z υυυυ υρ?????- =+++?????1z p w w w w f u w z t x y z υρ?????-=+++?????2 2 p v gz c ρ + + =2 2p v z c g g ρ+ +=22 11221222p v p v z z c g g g g ρρ+ +=++=22 1112221222p v p v z z g g g g ααρρ++=++ g V g p z g V g p z 222 2 2 2221111αρ αρ ++=++ 7．实际流体总流的伯努利方程式 8．粘性流体的伯努利方程 9．.总流的动量方程 10．总流的动量矩方程 11．叶轮机械的欧拉方程 22 1112221222f p v p v z z h g g g g ααρρ++=+++22 1 12212L 22p v p v z z h g g γγ++=+++∑=-F V Q V Q 111222ρβρβ∑?=?-?F r V r Q V r Q 11112222ρβρβ) cos cos (111222ααρr V r V Q M -=0dW d P=dt dt W Md M M M θ θθθω ====?功 功率 12．洒水器 第七章 流体在管路中的流动 1．临界雷诺数 临界雷诺数=2000，小于2000，流动为层流 大于2000，流动为湍流 2．沿程水头损失 当流动为层流时沿程水头损失hf 为， V(1.0) ； 当流动为湍流时沿程水头损失hf 为， V(1.75～2.0) 3．水力半径 相当直径 θωωθρcos 0 )cos (22R V R VR Q =∴=-V d Vd Re ρμ υ = = 12 f p p p h γ γ -?= = h A r P = h h 4d r = 4．圆管断面上的流量 5．平均流速 6．局部阻力因数为 7．管道沿程摩阻因数 8．沿程水头损失的计算 第九章 1．.薄壁孔口特征：L /d ≤2 厚壁孔口特征：2＜L /d ≤4 2．流速系数 4π 8Q GR μ = 2 max 2max 21π12 π82 R v Q G V R v A R μ====0 f 21 2 c V τρ= f 644c Re λ== f 2 8p Gl l h V R μγγγ?===226422l V l V Vd d g d g λρμ =??=??c v 11ζ+= C .3。流量系数Cd = CcCv D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ 1流体中稳定流动和均匀流动的区别 （1）①根据当地加速度是否为0，即流体运动要素是否随时间变化，流体分为 稳定流动和不稳定流动。 ②根据迁移加速度是否为0，即流体运动要素是否随空间参数变化，流体 分为均匀流和非均匀流。（非均匀流又分为缓变流和急变流） （2）稳定流动是流场中流体质点通过空间点时所有的运动要素都不随时间改变 的流动。 （3）均匀流动是指流场中同一直线上的各流体质点的运动要素沿程不变（不随 空间参数变化）的流动。 （4）稳定流的流线可以为曲线。均匀流的流线不能为曲线，只能是一元流动。 2迹线方程最后是写成多个还是整合成一个？ 答：如果迹线方程可以合并为一个，尽量合并为一个，并且尽量消掉参数t 。如果不能合并，就不用合并。理论上说都是可以的，但是从考试的答案来说，基本上都是合并的。 流体力学基本公式 1.牛顿内摩擦定律 （1）表达式： dy du μτ±=。 （2）内摩擦定律与三个因素相关，粘性切应力与流体粘度和速度梯度有关，与 压力的大小关系不大。 （3）适用条件：牛顿流体的层流运动。 2.欧拉平衡微分方程 （1）01=??-x p X ρ，01=??-y p Y ρ，01=??-z p Z ρ （2）适用于绝对静止状态和相对静止状态，可压缩流体和不可压缩流体。 3.静力学基本方程式 （1） g p z g p z ρρ2 211+=+ （2）适用条件：重力作用下、静止的、连通的、均质流体。 （3）几何意义：静止流体中，各点的测压管水头为常数。 （4）物理意义：静止流体中，各点的总比能为常数。 4.连续性方程 （1）适用于系统的质量守恒定律在控制体上的应用。 （2）三种形式：一般形式，恒定流，不可压缩流。 ①一般形式：0)()()(=??+??+??+??z u y u x u t z y x ρρρρ ②恒定流：0)()()(=??+??+??z u y u x u z y x ρρρ ③不可压缩流体：0=??+??+??z u y u x u z y x 5.欧拉运动方程 （1） dt du z p Z dt du y p Y dt du x p X z y x =??-=??-=??-ρρρ1,1,1 （2）适用条件：所有理想流体。 6．理想流体的伯努利方程 （1）2211221222p u p u z z g g g g ρρ++=++ （2）适用条件：理想流体；不可压缩流体；质量力只有重力；沿稳定流的流线 或微小流束。 （3）几何意义：沿流线总水头为常数。 （4）物理意义：沿流线总比能为常数。 7.实际流体总流的伯努利方程 （1）221112221222w p v p v z z h g g g g ααρρ++=+++ （2）适用条件：实际流体稳定流；不可压缩流体；质量力只有重力；所取断面 为缓变流断面。 （3）动能修正系数α：总流有效断面上的实际动能与按平均流速算出的假想动 能的比值。1α>，由断面上的速度分布不均匀引起，不均匀性越大，α越大。 8.动量方程 （1）() 21=Q F v v ρ-∑ 第二章 流体的主要物理性质 1．密度 ρ = m /V 7．压缩系数 T p V V ???? ? ?-=δδκ 体积模量 6．体胀系数 P V T V V ??? ??=δδα 9．牛顿内摩擦定律 h Av F /μ= dy dv x μ τ= 动力黏度：μ 运动黏度 ρμν= 第三章 流体静力学 重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。 1. 01=??-x p f x ρ 01=?-p ρf 2. 压强差公式 )(dz f dy f dx f dp z y x ++=ρ 等压面：dp =0 3.重力场中流体的平衡 4.帕斯卡定理 ()gh p z z g p p ρρ+=-+=000 5. 真空度 p p p a v -= 6. 等加速直线运动容器内液体的相对平衡 7.等角速度旋转容器中液体的相对平衡 C z g r g p +??? ? ??-=222ωρ 外加边界条件确定C 如：0,0,0p p z r === V P V K ??-=κ1 自由液面上某点的铅直坐标：g r Zs 22 2ω= 8.静止液体作用在平面上的总压力 9.静止液体作用在曲面上的总压力 水平方向的作用力：z x ghdA ghdA dF dF ρθρθ===cos cos 垂直方向的作用力 x z ghdA ghdA dF dF ρθρθ===sin sin 总压力 22y x F F F += z x F F tg = θ 第四章 流体运动学基础 1.．欧拉法 加速度场 简写为 当地加速度： 迁移加速度 2. 拉格朗日法：流体质点的运动速度的拉格朗日描述为 3.流线微分方程： 4．流量计算： 单位时间内通过d A 的微小流量为 d q v=u d A 通过整个过流断面流量 平均流速 5. 水力半径 ：总流的有效截面积与湿周之比 χ A R h = 6． ???' =V dV N ηρ 连续性方程 对于定常流动 r 1A 1u 1= r 2A 2u 2 对于不可压缩流体，r1 = r 2 =c A 1u 1=A 2u 2= q v υυ)(????==A A u q q d d v v 第1章 流体力学的基本概念 流体力学是研究流体的运动规律及其与物体相互作用的机理的一门专门学科。本章叙述在以后章节中经常用到的一些基础知识，对于其它基础内容在本科的流体力学或水力学中已作介绍，这里不再叙述。 连续介质与流体物理量 连续介质 流体和任何物质一样，都是由分子组成的，分子与分子之间是不连续而有空隙的。例如，常温下每立方厘米水中约含有3×1022 个水分子，相邻分子间距离约为3×10－8 厘米。因而，从微观结构上说，流体是有空隙的、不连续的介质。 但是，详细研究分子的微观运动不是流体力学的任务，我们所关心的不是个别分子的微观运动，而是大量分子“集体”所显示的特性，也就是所谓的宏观特性或宏观量，这是因为分子间的孔隙与实际所研究的流体尺度相比是极其微小的。因此，可以设想把所讨论的流体分割成为无数无限小的基元个体，相当于微小的分子集团，称之为流体的“质点”。从而认为，流体就是由这样的一个紧挨着一个的连续的质点所组成的，没有任何空隙的连续体，即所谓的“连续介质”。同时认为，流体的物理力学性质，例如密度、速度、压强和能量等，具有随同位置而连续变化的特性，即视为空间坐标和时间的连续函数。因此，不再从那些永远运动的分子出发，而是在宏观上从质点出发来研究流体的运动规律，从而可以利用连续函数的分析方法。长期的实践和科学实验证明，利用连续介质假定所得出的有关流体运动规律的基本理论与客观实际是符合的。 所谓流体质点，是指微小体积内所有流体分子的总体，而该微小体积是几何尺寸很小（但远大于分子平均自由行程）但包含足够多分子的特征体积，其宏观特性就是大量分子的统计平均特性，且具有确定性。 流体物理量 根据流体连续介质模型，任一时刻流体所在空间的每一点都为相应的流体质点所占据。流体的物理量是指反映流体宏观特性的物理量，如密度、速度、压强、温度和能量等。对于流体物理量，如流体质点的密度，可以地定义为微小特征体积内大量数目分子的统计质量除以该特征体积所得的平均值，即 V M V V ??=?→?'lim ρ （1-1） 式中，M ?表示体积V ?中所含流体的质量。 按数学的定义，空间一点的流体密度为 V M V ??=→?0 lim ρ （1-2） C3.6.2 达西摩擦因子 为了确定λ与Re 的关系,人们作了大量实验和理论研究，下面介绍有代表性的结果。 1．尼古拉兹实验 尼古拉兹（J.Nikuradse，1932）分析了达西的圆管沿程阻力实验数据后，发现壁面粗糙度对λ的影响很大，决定用人工粗糙度方法实现对粗糙度的控制。他用当地黄砂砂粒经筛选后分类均匀粘贴在管内壁上，相对粗糙度ε/d 从1/30—1/1014分6种，测得λ与Re 的关系，得到尼古拉兹图（图C3.6.1）。 2． 常用计算公式 从尼古拉兹图中看到在不同Re 数和ε/d 值的区域，λ有不同的变化规律。 图C3.6.1 （1）层流区 由泊肃叶定律推导的沿程水头损失（C3.4.10）式可得 代入达西公式（C3.6.3）式，可得层流区λ的解析式 上式表明层流区λ与管壁粗糙度无关，写成常用对数形式为 上式在双对数坐标系中是一条直线，与尼古拉兹图吻合。 （2）过渡区 该区是层流向湍流的转捩区（2000 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ，3 m kg 29.1=空气 ρ 牛顿内摩擦定律：剪切力：dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度：ρυμ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱m m 73610/9800012===m m N at 2/1013251m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D += 1 ) () 2(32121h h h h L e ++= 3 2L e y D = = （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形 12 3bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力 z x P P P += 与水平面的夹角 x z P P arctan =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ?? ?????+??+??+??=??+??+??+??= ??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ 工程流体力学公式总结 第二章流体得主要物理性质 ?流体得可压缩性计算、牛顿内摩擦定律得计算、粘度得三种表示方法。１.密度ρ= m／V 2.重度γ= G /V 3．流体得密度与重度有以下得关系:γ= ρg或ρ= γ/ g ４.密度得倒数称为比体积,以υ表示υ= １/ ρ= V/m 5.流体得相对密度:d = γ流/γ水= ρ流/ρ水 6．热膨胀性 ７.压缩性、体积压缩率κ 8．体积模量 ９.流体层接触面上得内摩擦力 10．单位面积上得内摩擦力(切应力)(牛顿内摩擦定律) 1１.、动力粘度μ： 12.运动粘度ν:ν=μ/ρ 1３．恩氏粘度°E:°E = t 1 /t 2 第三章流体静力学 ?重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体得压强计算、流体静压力得计算（压力体)。 1.常见得质量力: 重力ΔW = Δmg、 直线运动惯性力ΔFI ＝Δm·ａ 离心惯性力ΔＦR ＝Δm·rω２、 ２.质量力为F。:F= ｍ·aｍ= m(ｆｘi+f yj＋ｆzk) am ＝F/m = f xi＋f yj＋fｚk为单位质量力,在数值上就等于加速度 实例:重力场中得流体只受到地球引力得作用，取z轴铅垂向上,xoy为水平面，则单位质量力在x、y、z轴上得分量为 fｘ= ０,ｆy＝0 , fz=-mg/ｍ= －ｇ式中负号表示重力加速度g与坐标轴z方向相反 ３流体静压强不就是矢量,而就是标量,仅就是坐标得连续函数。即:p=p(x,y,z),由此得静压强得全微分为： 4.欧拉平衡微分方程式 单位质量流体得力平衡方程为: 1、单位质量力：m F f B B = 2、流体的运动粘度：ρ μ=v （μ[动力]粘度，ρ密度） 3、压缩系数：dp d dp dV V ρρκ?=?-=11（κ的单位是N m 2）体积模量为压缩系数的倒数 4、体积膨胀系数：dT d dT dV V v ρρα?-=?=11（v α的单位是C K ?1,1） 5、牛顿内摩擦定律：为液体厚）为运动速度，以应力表示为y u dy du dy du A T (,μτμ== 6、静止液体某点压强：为该点到液面的距离）h gh p z z g p p ()(000ρρ+=-+= 7、静水总压力： )h (为受压面积，为受压面形心淹没深度为静水总压力，A p ghA A p p c ρ== 8、元流伯努利方程;'2221112w h g p z g u g p z ++=++ρρ('w h 为粘性流体元流单位重量流体由过流断面1-1运动至过流断面2-2的机械能损失，z 为某点的位置高度或位置水头，g p ρ为测压管高度或压强水头，g u ρ2是单位流体具有的动能，u gh g p p g u 22'=-=ρ，u gh C g p p g C u 22'=-=ρC 是修正系数，数值接近于1） 9、总流伯努利方程:w h g v g p z g v g p z +++=++222 221221111αραρ(α为修正系数通常取1） 10、文丘里流量计测管道流量：)21)(41()()(42 122211g d d d k h k g p z g p z k Q -=?=+-+=πμρρμ 11、沿程水头损失一般表达式：g v d l h f 22 λ=（l 为管长，d 为管径，v 为断面平均流速，g 为重力加速度，λ为沿程阻力系数） 第一章流体力学基本知识 学习本章的目的和意义：流体力学基础知识是讲授建筑给排水的专业基础知识，只有掌握了该部分知识才能更好的理解建筑给排水课程中的相关内容。 §1-1 流体的主要物理性质 1．本节教学内容和要求： 1.1本节教学内容： 流体的4个主要物理性质。 1.2教学要求： （1）掌握并理解流体的几个主要物理性质 （2）应用流体的几个物理性质解决工程实践中的一些问题。 1.3教学难点和重点： 难点：流体的粘滞性和粘滞力 重点：牛顿运动定律的理解。 2．教学内容和知识要点： 2.1 易流动性 （1）基本概念：易流动性——流体在静止时不能承受切力抵抗剪切变形的性质称易流动性。 流体也被认为是只能抵抗压力而不能抵抗拉力。 易流动性为流体区别与固体的特性 2.2密度和重度 （1）基本概念：密度——单位体积的质量，称为流体的密度即： M ρ= V M——流体的质量，kg ; V——流体的体积，m3。 常温，一个标准大气压下Ρ水=1×103kg/ m3 Ρ水银=13.6×103kg/ m3 基本概念：重度：单位体积的重量，称为流体的重度。重度也称为容重。 G γ= V G——流体的重量，N ; V——流体的体积，m3。 ∵G=mg ∴γ=ρg 常温，一个标准大气压下γ水=9.8×103kg/ m3 γ水银=133.28×103kg/ m3密度和重度随外界压强和温度的变化而变化 液体的密度随压强和温度变化很小，可视为常数，而气体的密度随温度压强变化较大。 2..3 粘滞性 （1）粘滞性的表象 基本概念：流体在运动时抵抗剪切变形的性质称为粘滞性。当某一流层对相邻流层发生位移而引起体积变形时，在流体中产生的切力就是这一性质的表 现。 为了说明粘滞性由流体在管道中的运动速度实验加以分析说明。用流速仪测出管道中某一断面的流速分布如图一所示 设某一流层的速度为u,则与其相邻的流层为u+du,du为相邻流层的速度增值，设相邻流层的厚度为dy，则du/dy叫速度梯度。 由于各流层之间的速度不同，相邻流层间有相对运动，便在接触面上产生一种相互作用的剪切力，这个力叫做流体的内摩擦力，或粘滞力。 平板实验 （2）牛顿内摩擦定律 基本概念：牛顿在平板实验的基础上于1867年在所著的《自然哲学的数学原理》中提出了流体内摩擦力的假说——牛顿内摩擦定律： 当切应力一定时，粘性越大，剪切变形的速度越小，所以粘性又可定义为流体 工程流体力学公式总结 第二章 流体的主要物理性质 流体的可压缩性计算、牛顿内摩擦定律的计算、粘度的三种表示方法。 1．密度 ρ = m /V 2．重度 γ = G /V 3．流体的密度和重度有以下的关系：γ = ρ g 或 ρ = γ/ g 4．密度的倒数称为比体积，以υ表示υ = 1/ ρ = V/m 5．流体的相对密度：d = γ流 /γ水 = ρ流 /ρ水 6．热膨胀性 7．压缩性. 体积压缩率κ 8．体积模量 9．流体层接触面上的内摩擦力 10．单位面积上的内摩擦力（切应力）（牛顿内摩擦定律） 11．.动力粘度μ： 12．运动粘度ν ：ν = μ/ρ 13．恩氏粘度°E ：°E = t 1 / t 2 第三章 流体静力学 重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。 1．常见的质量力： 重力ΔW = Δmg 、 直线运动惯性力ΔFI = Δm·a 离心惯性力ΔFR = Δm·r ω2 . T V V ??=1αp V V ??-=1κV P V K ??-=κ1n A F d d υμ=dn d v μτ±=n v d /d τμ= 2．质量力为F 。：F = m ·am = m (f xi+f yj+f zk) am = F /m = f xi+f yj+f zk 为单位质量力，在数值上就等于加速度 实例：重力场中的流体只受到地球引力的作用，取z 轴铅垂向上，xoy 为水平面，则单位质量力在x 、y 、 z 轴上的分量为 fx = 0 , fy = 0 , fz = -mg /m = -g 式中负号表示重力加速度g 与坐标轴z 方向相反 3流体静压强不是矢量，而是标量，仅是坐标的连续函数。即：p = p (x ,y ,z )，由此得静压强的全微分为: 4．欧拉平衡微分方程式 单位质量流体的力平衡方程为： 5．压强差公式（欧拉平衡微分方程式综合形式） 6．质量力的势函数 7．重力场中平衡流体的质量力势函数 z z p y y p x x p p d d d d ??????++=d d d d d d 0x p f x y z x y z x ??-=ρd d d d d d 0y p f x y z x y z y ??-=ρd d d d d d 0z p f x y z x y z z ??-=ρ0 1=??-x p f x ρ10y p f y ??-=ρ01=??-z p f z ρz z p y y p x x p z f y f x f z y x d d d )d d d (??+??+??=++ρ) d d d (d z f y f x f p z y x ++=ρd (d d d )x y z p f x f y f z dU ρ=++=ρd d d d x y z U U U U x y z =f dx f dy f dz x y z gdz ??????=++++=- 计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区 《流体力学与流体机械》(上)主要公式及方程式 1．流体的体积压缩系数计算式：β1dρ p=-1dV Vdp=ρdp 流体的体积弹性系数计算式：E=-Vdpdp dV=ρdρ 流体的体积膨胀系数计算式：βdV T=1 VdT=-1dρ ρdT 2．等压条件下气体密度与温度的关系式：ρ0 t=ρ 1+βt，其中β=1 273。 3T=±μAdu dy 或τ=Tdu A=±μdy 恩氏粘度与运动粘度的转换式：ν=(0.0731E-0.0631 E)?10-4 f1?p? x-ρ?x=0?fr-1?p=0? ?ρ?r?? 4．欧拉平衡微分方程式： f? y-1?p ρ?y=0??和fθ-1?p ρ=0? f1?p?r?θ ρ?z=0?? ??f1?p? z-z-ρ?z=0?? 欧拉平衡微分方程的全微分式：dp=ρ(fxdx+fydy+fzdz) dp=ρ(frdr+fθrdθ+fzdz) 5 fxdx+fydy+fzdz=0 frdr+fθrdθ+fzdz=0 6p γ+z=C 或 p1 γ+zp21=γ+z2 或p1+ρgz1=p2+ρgz2 相对于大气时：pm+(ρ-ρa)gz=C 或pm1+(ρ-ρa)gz1=pm2+(ρ-ρa)gz2 7p=p0+γh，其中p0为自由液面上的压力。 8．水平等加速运动液体静压力分布式：p=p0-ρ(ax+gz)；等压面方程式： ax+gz=C；自由液面方程式：ax+gz=0。注意：p0为自由液面上的压力。 1 9．等角速度旋转液体静压力分布式：p=p0+γ(ω2r2 2g-z)；等压面方程式：ω2r2 2-gz=C；自由液面方程式：ω2r2 2-gz=0。注意：p0为自由液面上的压力。 10．静止液体作用在平面上的总压力计算式：P=(p0+γhc)A=pcA，其中p0为自由液面上的相对压力。压力中心计算式：yD=yc+γsinαIxc (p0+γycsinα)A Ixc ycA或yD-yc=Ixc ycA。当自由液面上的压力为大气压时：yD=yc+ 矩形截面的惯性矩Ixc计算式：Ixc= 圆形截面的惯性矩Ixc计算式：Ixc11bh3；三角形截面的惯性矩Ixc计算式：Ixc=bh3 1236π4=d 64 11．静止液体作用在曲面上的总压力的垂直分力计算式：Pz=p0Az+γVP，注意：式中p0应为自由液面上的相对压力。 12 ?ux?ux?ux?ux?+ux+uy+uz?τ?x?y?z???uy?uy?uy?uy?+ux+uy+uz直角坐标系：ay=? ?τ?x?y?z??u?uz?uz?uz?az=z+ux+uy+uz?τ?x?y?z??ax= ?ur?ur?ur?uruθ2ar=+ur+uθ+uz-?τ?rr?θ?zr ?u?u?u?uuu圆柱坐标系：aθ=θ+urθ+uθθ+uzθ+rθ ?τ?rr?θ?zr ?u?uz?uz?uzaz=z+ur+uθ+uz?τ?rr?θ?z????????? 流体质点的压力、密度等流动参量对时间的变化率计算式： dp?p?p?p?p=+ux+uy+uzdτ?τ?x?y?z dρ?ρ?ρ?ρ?ρ=+ux+uy+uz?τ?x?y?z dτ 13 drrdθdzdxdydz==== 及uxuyuzuruθuz2 ?ρ?(ρux)?(ρuy)?(ρuz)14．三维连续性方程式的一般式：+++=0 ?τ?x?y?z ?ρρur?(ρur)?(ρuθ)?(ρuz)++++=0 ?τr?rr?θ?z ?ux?uy?uz15．不可压缩流体的三维连续性方程式：++=0 ?x?y?z ur?ur?uθ?uz+++=0?rr?θ?z r 16M=ρ11A1=ρ22A2 对于不可压缩流体： Q=1A1=2A2 随体倒数 ()D u D t t ααα?= +??? ()() u u i v j w k i j k u v w x y z x y z ??????????=++?++=++ ????????? 雷诺输运定理：对系统的随体倒数求法 [()][ )]V V k V V k D dv u dv D t t D dv u dv D t t x φφφφφφ?=+????? = +?????? （ ij i j e e δ=? ()i j k i jkl l jkl il jki ijk e e e e e εεδεε??=?=== i j ijk k e e e ε?= ()()()() i j i j i j i j i i e e e e x x x x x x φ φφφ?????????=?=?=?????? ()i i i i e e x x φφφ???==?? ()i i j j i i a a e a e x x ??????=?= ????? ()()j j k i j j i j ijk k ijk i i i i j a a a a e a e e e e e x x x x εε??????=?=?==???? 1、i j u x ?? ?? ?????? ：速度梯度张量 应变率张量：表示微团的变形运动 1122112211 22ij u u v u w x y x z x v u v v w s x y y z y w u w v w x z y z z ?? ?? ???????++ ? ? ?????????? ? ? ?? ?? ????? ?=++ ? ?????? ? ???? ? ?? ??????? ? ++ ? ? ???????? ?? ? ? 旋转张量：表示旋转 3231 210 0 0ij a ωωωωωω-?? ?= ? ?-?? － 质量守恒： ()0k k u t x ρρ??+=?? 0k k u D D t x ρρ ?+=? 第二那诺雷诺输运定律： V V D D dv dv D t D t αραρ= ? ? 动量守恒定律：() u u u f t ρ ρρ?+??=??+?σ ij i i j D u f D t x σρ ρ?= +? ij i i j i j j u u u f t x x σρ ρρ???+= +??? D u f D t ρ ρ=??+σ 能量守恒定律：()1 2i i i j ij i i i i q D e u u u u f D t x x ρ σρ???? +=+- ????? 231a ω=-312 a ω=-123a ω=-ij ijk k a εω=- 2、流体的运动粘度: ［动力］粘度， 密 度） 5、牛顿内摩擦定律： T A ，以应力表示为 （u 为运动速度，y 为液体厚） dy dy 6、静止液体某点压强： p P o g （z o z ） p o gh （h 为该点到液面的距离） 7、静水总压力: 10、文丘里流量计测管道流量: -、 2g) 1 11、沿程水头损失一般表达式： h f 1 V （ l 为管长，d 为管径，v 为断面平均流速，g d 2g 1单位质量力: F B 3、压缩系数: 1?dV V dp 丄？d dp 的单位是m %）体积模量为压缩系数的倒数 4、体积膨胀系数: v 1?dV V dT （V 的单位是 1K ,1 C ） p P c A ghA （p 为静水总压力, h 为受压面形心淹没深度 ，A 为受压面积） 8、元流伯努利方程；乙旦 g 2 U 1 2g Z 2 虫 h w' （h w'为粘性流体元流单位重量流体由过流 g 断面1-1运动至过流断面2-2 的机械能损失，z 为某点的位置高度或位置水头, 管高度或压强水头, 2 —是单位流体具有的动能, g u fg 晋丽， u C 2g p g p C 2gh u C 是修正系数，数值接近于 9、总流伯努利方程 2 1 V 1 z Z 2 2g R L g 2 2V 2 h w （为修正系数通常取1） (Z 2 为重力加速度， 为沿程阻力系数） 12、局部水头损失一般表达式: 2 h j —（为局部水头损失系数， v 为 对应的断面平均流速） J 2g .-pl 13、圆管流雷诺数：R e 一（u 为流速，V 为运动粘度，d 为圆管直径） V uR 14、非圆管道流雷诺数： R e （R 为水力半径，水力半径R V 渠宽度，h 为明渠水深） 力坡度，J 牛） 半径，J '为所取流束的水力坡度，与总水流坡度相等） 17、过流断面上的流速分布的解析式： u J （r ； r 2） 4 18、平均流速：v Q A Q 2 r 。 8 r0 ，断面平均流速与最大流速的关系： 1 v U max 2 19、沿程水头损失： h f 64 l v 2 l 2 爲g ，其中为沿程摩阻系数 ，沿程摩阻系数 Re d 2g 64 Re 20、谢才公式：V 8g . RJ C ? RJ （v 为断面平均流速，R 为水力半径，J 为水力坡 度，C 为谢才系数） A A 为过流断面面积，x 为过流断面上流体与固体接触的周界, 矩形断面明渠流的水力半径: R 一 ，b 为明 b 2h 15、均匀流动方程式： h f l gA gR? gRJ （R 为水力半径，J 为水 16、流束的均匀流动方程: gRJ （为所取流束表面的剪应力, R'为所取流束的水力 21、曼宁公式: 1 -R n 1 0.5 6（吹） （n 为综合反映壁面对水流阻滞作用的系数，称为粗糙 第二章 流体的主要物理性质 流体的可压缩性计算、牛顿内摩擦定律的计算、粘度的三种表示方法。 1．密度 ρ = m /V 2．重度 γ = G /V 3．流体的密度和重度有以下的关系：γ = ρ g 或 ρ = γ/ g 4．密度的倒数称为比体积，以υ表示υ = 1/ ρ = V/m 5．流体的相对密度：d = γ流 /γ水 = ρ流 /ρ水 6．热膨胀性 7．压缩性. 体积压缩率κ 8．体积模量 9．流体层接触面上的内摩擦力 10．单位面积上的内摩擦力（切应力）（牛顿内摩擦定律） 11．.动力粘度μ： 12．运动粘度ν ：ν = μ/ρ 13．恩氏粘度°E ：°E = t 1 / t 2 第三章 流体静力学 重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、 压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。 1．常见的质量力： 重力ΔW = Δmg 、 直线运动惯性力ΔFI = Δm·a 离心惯性力ΔFR = Δm·r ω2 . 2．质量力为F 。：F = m ·am = m (f xi+f yj+f zk) am = F /m = f xi+f yj+f zk 为单位质量力，在数值上就等于加速度 实例：重力场中的流体只受到地球引力的作用，取z 轴铅垂向上，xoy 为水平面，则单位质量力在x 、y 、 z 轴上的分量为： fx = 0 , fy = 0 , fz = -mg /m = -g 式中负号表示重力加速度g 与坐标轴z 方向相反 3流体静压强不是矢量，而是标量，仅是坐标的连续函数。即：p = p (x ,y ,z )，由此得静压强的全微分为: 4．欧拉平衡微分方程式 T V V ??=1αp V V ??-=1κV P V K ??- = κ1 n A F d d υ μ=dn d v μτ±=n v d /d τ μ=z z p y y p x x p p d d d d ??????+ + = d d d d d d 0 x p f x y z x y z x ??- =ρd d d d d d 0 y p f x y z x y z y ??- =ρd d d d d d 0 z p f x y z x y z z ??- =ρ 流体力学 一、流体的主要物性与流体静力学 1、静止状态下的流体不能承受剪应力，不能抵抗剪切变形。 2、粘性：内摩擦力的特性就是粘性，也是运动流体抵抗剪切变形的能力，是运动流体产生机械能损失的根源；主要与流体的种类和温度有关，温度上升粘性减小，与压强没关系。 3、牛顿内摩擦定律： du F A dy μ= F d u A d y τμ== 相关因素：粘性系数、面积、速度、距离；与接触面的压力没有关系。 例1：如图6-1所示，平板与固体壁面间间距为1mm,流体的动力黏滞系数为0.1Pa.S, 以50N 的力拖动，速度为1m/s,平板的面积是（ ）m 2。 解：F F A du dy δ μνμ= ==0.5 例2：如图6-2所示，已知活塞直径d=100mm,长l=100mm 气缸直径D=100.4mm,其间充满黏滞系数为0.1Pa·s 的油，活塞以2m/s 的速度运动时，需要的拉力F 为（ ）N 。 解：332 0.1[(10010)0.1]31.40.210 du F A N dy μπ--==?????=? 4、记忆个参数，常温下空气的密度3 1.205/m kg ρ=。 5、表面力作用在流体隔离体表面上，起大小和作用面积成正比，如正压力、剪切力；质量 力作用在流体隔离体内每个流体微团上，其大小与流体质量成正比，如重力、惯性力，单位质量力的单位与加速度相同，是2 /m s 。 6、流体静压强的特征： A 、垂直指向作用面，即静压强的方向与作用面的内法线方向相同； B 、任一点的静压强与作用面的方位无关，与该点为位置、流体的种类、当地重力加速度等因素有关。 7、流体静力学基本方程 0p p gh ρ=+ 2198/98at kN m kPa == 流体微团运动分析 加速度 : 欧拉法的加速度三个分量 z u u y u u x u u t u Dt Du a y z y y y x y y y ??+??+??+??= = z u u y u u x u u t u Dt Du a z z z y z x z z z ??+??+??+??= = z u u y u u x u u t u Dt Du a x z x y x x x x x ??+??+??+??= = u u t u Dt u D a (??+??= = 哈密顿算子 t k t j t i ??+?? +??=? 1. 线变形 (1线应变率(线变形速度 : (2面积扩张率 : 流体面元的面积在平面内的局部瞬时相对扩张速率(3体积膨胀率 :流体体元的体积在空间的局部瞬时相对膨胀速率 x u x xx ??= εy u y yy ??= ε z u z zz ??= ε y u x u u y x ??+ ??=?? z u y u x u u z y x ??+ ??+??=??????????+??=y u x u x y xy 21ε??? ? ????+??=y u z u z y yz 21ε?? ? ????+??=z u x u x z zx 21ε2. 角变形速度:单位时间直角边的偏转角度之半为流体微团的的角变形速度。 3 流体的旋转(旋转运动 ? 旋转角速度 : 两正交线元在 xy 面内 绕一点的旋转角速度平均值??? ? ????- ??=y u x u x y z 21ω(规定逆时针方向为正 流体力学的计算公式 众所周知，锅炉之类的热力设备可以采用水封管作为限压装置。设备工作时，水封管内的水将设备与大气隔绝，使设备保持一定的压力。水封管内水的多少，即水封管内的实际水位的高低与设备汽压大小有关。汽压大，水位就高；汽压小，水位就低。当设备压力升至额定蒸汽压力时，水封管内的水位应该达到最高水位；而在设备压力稍微大于额定蒸汽压力时，水封管内的水应能立即被冲除掉，使设备内的蒸汽能够迅速得到排放。但是，水封管的内径必须足够大；以保证它的排汽能力大于设备的最大蒸发量，从而防止设备发生超压事故。《蒸汽锅炉安全技术监察规程》第132条规定：“对于额定蒸汽压力小于等于0.1MPa的锅炉可以采用静重式安全阀或水封式安全装置。水封装置的水封管内径不应小于25mm”。《小型和常压热水锅炉安全监察规定》第二十八条要求：“水封管的内径应当根据锅炉的额定容量和压力确定，且内径不得小于25毫米”。 但是，水封管的内径究竟应该根据什么样的公式去计算？这里没有具体加以规定。水封管的内径为何不能小于25mm？这里没有加以说明。 1 确定水封管内径的一种错误观点 有人说，可以按安全阀的排放量计算公式去间接确定水封管的内径。笔者认为这是不对的。理由如下： 其一，按流体力学，安全阀的局部阻碍因为其进出口呈直角型式，边壁是突变的，主流与边壁之间形成大尺度旋涡，蒸汽排放时能量损失很大；而水封管的边壁是渐变的，又不出现减速增压现象的部位，故蒸汽排放时的能量损失很小。其二，安全阀的出口之外肯定还有排汽管，它的排放量计算公式应该考虑这一额外的因素；而水封管却只相当于它的排汽管。其三，安全阀在排汽时，汽流还自始至终受到弹簧或者重锤施加的反方向作用力，能量损失非常之大；而水封管在排汽时却不会受到这样的反作用力，无此项能量损失。显然，安全阀排汽与水封管排汽区别太大，不能按安全阀排放量计算公式去确定水封管的内径。 2 用流体力学理论建立水封管内径计算方法 为了正确确定水封管的应有内径，我们有必要按流体力学的理论来分析一下水封管的流动阻力和能量损失。蒸汽流过水封管时的阻力损失ΔP等于水封管进口压力P1与出口压力P2之差。管道的阻力损失 ΔP＝Pa（1） 式中λ—管道摩擦系数； d —管道内径，m； ρ—流体的密度，k g／m3； u —流体在管内的流速，m／s； L—平直管段的总长度，m； Σξ—管道中各种局部阻力之和。 管道的摩擦系数λ值取决于流体流动的雷诺数Re 和管壁的相对粗糙度。管壁的相对粗糙度等于管壁的绝对粗糙度（即管壁内凸起高度）K与管内径d之比。雷诺数Re =(2) 式中ν－流体运动粘滞系数，m2／s； 其他符号说明同上。 由于水封管进出口压差通常大于40000Pa，排汽时的流速都大于100m／s ，雷诺数Re 都大于80000000，所以蒸汽的流动状态为紊流流动，而且位于阻力平方区。在此区域内，摩擦阻力系数λ值仅取决于管壁的相对粗糙度，可以用尼古拉兹公式进行计算：λ＝（3）流体力学复习要点(计算公式)
流体力学基本公式
工程流体力学公式
第1章流体力学的基本概念
流体力学计算公式
流体力学复习要点(计算公式)
流体力学公式总结
流体力学计算公式
流体力学基础学习知识知识
流体力学公式总结
计算流体力学常用数值方法简介[1]
《流体力学》Ⅰ主要公式及方程式讲解
流体力学-公式
流体力学计算公式
流体力学公式
流体力学总结笔记
流体力学公式总结.
流体力学的计算公式