第四节管路系统的计算-欢迎访问河南城建学院化学与化学工程

11.4流体流动阻力本节重点：直管阻力与局部阻力的计算，摩擦系数的影响因素。难点：用因次分析法解决工程实际问题。流动阻力的大小与流体本身的物理性质、流动状况及壁面的形状等因素有关。化工管路系统主要由两部分组成，一部分是直管，另一部分是管件、阀门等。相应流体流动阻力也分为两种：直管阻力：流体流经一定直径的直管时由于内摩擦而产生的阻力；局部阻力：流体流经管件、阀门等局部地方由于流速大小及方向的改变而引起的阻力。1.4.1流体在直管中的流动阻力1.阻力的表现形式如图1-24所示，流体在水平等径直管中作定态流动。在1-1′和2-2′截面间列柏努利方程，fWpugzpugz222212112121因是直径相同的水平管，21uu21zz21ppWf（1-34）若管道为倾斜管，则)()(2211gzpgzpWf（1-34a）由此可见，无论是水平安装，还是倾斜安装，流体的流动阻力均表现为静压能的减少，仅当水平安装时，流动阻力恰好等于两截面的静压能之差。2.直管阻力的通式在图1-24中，对1-1′和2-2′截面间流体进行受力分析：由压力差而产生的推动力为4221dpp与流体流动方向相同流体的摩擦力为dlAF与流体流动方向相反。2流体在管内作定态流动，在流动方向上所受合力必定为零。dldpp4)(221整理得dlpp421（1-35）将式(1-35)代入式（1-34）中，得dlWf4（1-36）将式（1-36）变形，把能量损失fW表示为动能22u的某一倍数。2822udluWf令28u则22udlWf（1-37）式（1-37）为流体在直管内流动阻力的通式，称为范宁（Fanning）公式。式中为无因次系数，称为摩擦系数或摩擦因数，与流体流动的Re及管壁状况有关。根据柏努利方程的其它形式，也可写出相应的范宁公式表示式：压头损失gudlhf22（1-37a）压力损失22udlpf（1-37b）值得注意的是，压力损失fp是流体流动能量损失的一种表示形式，与两截面间的压力差)(21ppp意义不同，只有当管路为水平时，二者才相等。应当指出，范宁公式对层流与湍流均适用，只是两种情况下摩擦系数不同。以下对层流与湍流时摩擦系数分别讨论。3.层流时的摩擦系数3流体在直管中作层流流动时，管中心最大速度如式（1-35）所示。将平均速度max21uu及2dR代入上式中，可得22132)(dlupp232dlupf（1-38）式（1-38）称为哈根-泊谡叶（Hagen-Poiseuille）方程，是流体在直管内作层流流动时压力损失的计算式。结合式（1-34），流体在直管内层流流动时能量损失或阻力的计算式为232dluWf（1-39）表明层流时阻力与速度的一次方成正比。式（1-39）也可改写为2Re6426432222udludluddluWf（1-39a）将式（1-39a）与式（1-37）比较，可得层流时摩擦系数的计算式Re64（1-40）即层流时摩擦系数λ是雷诺数Re的函数。4．湍流时的摩擦系数（1）因次分析法层流时阻力的计算式是根据理论推导所得，湍流时由于情况要复杂得多，目前尚不能得到理论计算式，但通过实验研究，可获得经验关系式，这种实验研究方法是化工中常用的方法。在实验时，每次只能改变一个变量，而将其它变量固定，如过程涉及的变量很多，工作量必然很大，而且将实验结果关联成形式简单便于应用的公式也很困难。若采用化工中常用的工程研究方法——因次分析法，可将几个变量组合成一个无因次数群（如雷诺数Re即是由d、ρ、u、μ四个变量组成的无因次数群），用无因次数群代替个别的变量进行实验，由于数群的数目总是比变量的数目少，就可以大大减少实验的次数，关联数据的工作也会有所简化，而且可将在实验室规模的小设备中用某种物料实验所得的结果应用到其它物料及实际的化工4设备中去。因次分析法的基础是因次一致性原则，即每一个物理方程式的两边不仅数值相等，而且每一项都应具有相同的因次。因次分析法的基本定理是白金汉（Buckinghan）的π定理：设影响某一物理现象的独立变量数为n个，这些变量的基本因次数为m个，则该物理现象可用N＝（n－m）个独立的无因次数群表示。根据对摩擦阻力性质的理解和实验研究的综合分析，认为流体在湍流流动时，由于内摩擦力而产生的压力损失fp与流体的密度ρ、粘度μ、平均速度u、管径d、管长l及管壁的粗糙度ε有关，即,,,,,ldufpf（1-41）7个变量的因次分别为：[p]=M-2L-1[]=ML-3[u]=M-1[d]=L[l]=L[]=L[]=M-1L-1基本因次有3个。根据π定理，无因次数群的数目N=n-m=7-3=4个将式（1-41）写成幂函数的形式：fedcbafulkdp因次关系式：fedcbaLMLMLLLLLM)1()()(13112根据因次一致性原则：对于M：ed＋1对于L：fedcba31＋对于：ec2设b，e，f已知，解得：5edecfeba12feeebfebfulkdp12febfduddlkup2即ddludupf,,2（1-42）式中ud——雷诺数Re，2upf——欧拉（Euler）准数，也是无因次数群。dl、d均为简单的无因次比值，前者反映了管子的几何尺寸对流动阻力的影响，后者称为相对粗糙度，反映了管壁粗糙度对流动阻力的影响。式（1-42）具体的函数关系通常由实验确定。根据实验可知，流体流动阻力与管长l成正比，该式可改写为：ddlupfRe,2（1-43）或2Re,uddlpWff（1-43a）与范宁公式（1-37）相对照，可得)(Re,d（1-44）即湍流时摩擦系数λ是Re和相对粗糙度d的函数，如图1-25所示，称为莫狄（Moody）摩擦系数图。6根据Re不同，图1-25可分为四个区域；（1）层流区（Re≤2000），λ与d无关，与Re为直线关系，即Re64，此时uWf，即fW与u的一次方成正比。（2）过渡区（2000Re4000），在此区域内层流或湍流的λ～Re曲线均可应用，对于阻力计算，宁可估计大一些，一般将湍流时的曲线延伸，以查取λ值。（3）湍流区（Re≥4000以及虚线以下的区域），此时λ与Re、d都有关，当d一定时，λ随Re的增大而减小，Re增大至某一数值后，λ下降缓慢；当Re一定时，λ随d的增加而增大。（4）完全湍流区（虚线以上的区域），此区域内各曲线都趋近于水平线，即λ与Re无关，只与d有关。对于特定管路d一定，λ为常数，根据直管阻力通式可知，2uWf，所以此区域又称为阻力平方区。从图中也可以看出，相对粗糙度d愈大，达到阻力平方区的Re值愈低。对于湍流时的摩擦系数λ，除了用Moody图查取外，还可以利用一些经验公式计算。这里介绍适用于光滑管的柏拉修斯（Blasius）式：25.0Re3164.0（1-45）其适用范围为Re＝5×103～105。此时能量损失fW约与速度u的1.75次方成正比。图1-25摩擦系数λ与雷诺数Re及相对粗糙度d的关系7考莱布鲁克（Colebrook）式此式适用于湍流区的光滑管与粗糙管直至完全湍流区。（2）管壁粗糙度对摩擦系数的影响光滑管：玻璃管、铜管、铅管及塑料管等称为光滑管；粗糙管：钢管、铸铁管等。管道壁面凸出部分的平均高度，称为绝对粗糙度，以表示。绝对粗糙度与管径的比值即d，称为相对粗糙度。工业管道的绝对粗糙度数值见教材。管壁粗糙度对流动阻力或摩擦系数的影响，主要是由于流体在管道中流动时，流体质点与管壁凸出部分相碰撞而增加了流体的能量损失，其影响程度与管径的大小有关，因此在摩擦系数图中用相对粗糙度d，而不是绝对粗糙度。流体作层流流动时，流体层平行于管轴流动，层流层掩盖了管壁的粗糙面，同时流体的流动速度也比较缓慢，对管壁凸出部分没有什么碰撞作用，所以层流时的流动阻力或摩擦系数与管壁粗糙度无关，只与Re有关。流体作湍流流动时，靠近壁面处总是存在着层流内层。如果层流内层的厚度δL大于管壁的绝对粗糙度，即δL时，如图1-26（a）所示，此时管壁粗糙度对流动阻力的影响与层流时相近，此为水力光滑管。随Re的增加，层流内层的厚度逐渐减薄，当δL时，如图1-26（b）所示，壁面凸出部分伸入湍流主体区，与流体质点发生碰撞，使流动阻力增加。当Re大到一定程度时，层流内层可薄得足以使壁面凸出部分都伸到湍流主体中，质点碰撞加剧，致使粘性力不再起作用，而包括粘度μ在内的Re不再影响摩擦系数的大小，流动进入了完全湍流区，此为完全湍流粗糙管。图1-26流体流过管壁面的情况Re7.182log274.11d8例分别计算下列情况下，流体流过φ76×3mm、长10m的水平钢管的能量损失、压头损失及压力损失。（1）密度为910kg/、粘度为72cP的油品，流速为1.1m/s；（2）20℃的水，流速为2.2m/s。解：（1）油品：200097310721.191007.0Re3ud流动为层流。摩擦系数可从图1-25上查取，也可用式（1-40）计算：0658.097364Re64所以能量损失J/kg69.521.107.0100658.0222udlWf压头损失m58.081.969.5gWhff压力损失Pa517869.5910ffWp（2）20℃水的物性：3kg/m2.998，310005.1Pa·s531053.110005.12.22.99807.0udRe流动为湍流。求摩擦系数尚需知道相对粗糙度d，查表1-2，取钢管的绝对粗糙度为0.2mm，则00286.0702.0d根据Re=1.53×105及d＝0.00286查图1-25，得λ＝0.027所以能量损失J/kg33.922.207.010027.0222udlWf压头损失m95.081.933.9gWhff压力损失Pa931333.92.998ffWp5.非圆形管道的流动阻力9对于非圆形管内的湍流流动，仍可用在圆形管内流动阻力的计算式，但需用非圆形管道的当量直径代替圆管直径。当量直径定义为Ade44＝润湿周边流通截面积(1-46)对于套管环隙，当内管的外径为d1，外管的内径为d2时，其当量直径为1212212244ddddddde对于边长分别为a、b的矩形管，其当量直径为baabbaabde2)(24在层流情况下，当采用当量直径计算阻力时，还应对式（1-40）进行修正，改写为ReC（1-47）式中C为无因次常数。一些非圆形管的C值见教材。注意，当量直径只用于非圆形管道流动阻力的计算，而不能用于流通面积及流速的计算。1.4.2局部阻力局部阻力有两种计算方法：阻力系数法和当量长度法。1.阻力系数法克服局部阻力所消耗的机械能，可以表示为动能的某一倍数，即22'uWf（1-48）或guhf22'（1-48a）式中ζ称为局部阻力系数，一般由实验测定。常用管件及阀门的局部阻力系数见教材。注意表中当管截面突然扩大和突然缩小时，式（1-48）及（1-48a）中的速度u均以小管中的速度计。当流体自容器进入管内，5.0进口，称为进口阻力系数；当流体自管子进入容器或从10管子排放到管外空间，1出口，称为出口阻力系数。当流体从管子直接排放到管外空间时，管出口内侧截面上的压强可取为与管外空间相同，但出口截面上的动能及出口阻力应与截面选取相匹配。若截面取管出口内侧，则表示流体并未离开管路，此时截面上仍有动能，