流体流动阻力的测定

  • 格式:docx
  • 大小:24.92 KB
  • 文档页数:7

实验名称:流体流动阻力的测定

一、实验目的及任务:

1. 掌握测定流体流动阻力实验的一般方法。

2. 测定直管的摩擦阻力系数及突然扩大管的局部阻力系数。

3. 验证湍流区内摩擦阻力系数为雷诺数和相对粗糙度的函数。

4. 将所得光滑管的方程与Blasius方程相比较。

二、实验原理:

流体输送的管路由直管和阀门、弯头、流量计等部件组成。由于粘性和涡流作用,流体在输送过程中会有机械能损失。这些能量损失包括流体流经直管时的直管阻力和流经管道部件时的局部阻力,统称为流体流动阻力。

1. 根据机械能衡算方程,测量不可压缩流体直管或局部的阻力

𝐻𝑓=(𝑔𝑧1+𝑝1𝜌+𝑢122)−(𝑔𝑧2+𝑝2𝜌+𝑢222)+𝐻𝑒

如果管道无变径,没有外加能量,无论水平或倾斜放置,上式可简化为:

𝐻𝑓=𝑝1′−𝑝2′𝜌=?𝑝𝜌

Δp为截面1到2之间直管段的虚拟压强差,即单位体积流体的总势能差,通过压差传感器直接测量得到。

2. 流体流动阻力与流体性质、流道的几何尺寸以及流动状态有关,可表示为:

?𝑝=𝑓(𝑑,𝑙,𝑢,𝜌,𝜇,𝜀)

由量纲分析可以得到四个无量纲数群:

欧拉数Eu=?𝑝/𝜌𝑢2,雷诺数𝑅𝑒=𝑑𝑢𝜌/𝜇,相对粗糙度ε/d和长径比𝑙/𝑑

从而有

?𝑝𝜌𝑢2=𝛹(𝑑𝑢𝜌𝜇,𝜀𝑑,𝑙𝑑)

取𝜆=𝛷(𝑅𝑒,𝜀/𝑑),可得摩擦系数与阻力损失之间的关系: 𝐻𝑓=?𝑝𝜌=𝜆𝑙𝑑×𝑢22

从而得到实验中摩擦系数的计算式

𝜆=2?𝑝𝑑𝜌𝑢2𝑙

当流体在管径为d的圆形管中流动时,选取两个截面,用压差传感器测出两个截面的静压差,即可求出流体的流动阻力。根据伯努利方程摩擦系数与静压差的关系,可以求出摩擦系数。改变流速可测得不同Re下的λ,可以求出某一相对粗糙度下的λ-Re关系。

在湍流区内摩擦系数𝜆=𝛷(𝑅𝑒,𝜀/𝑑),对于光滑管(水力学光滑),大量实验证明,Re在103~105氛围内,λ与Re的关系遵循Blasius关系式,即

𝜆=0.3163/𝑅𝑒0.25

对于粗糙管,λ与Re的关系以图来表示。

3. 对局部阻力,可用局部阻力系数法表示:

4. 𝐻𝑓= 𝜁𝑢22

对于扩大和缩小的直管,式中的流速按照细管的流速来计算。

对一段突然扩大的圆直管,局部阻力远大于其直管阻力。由忽略直管阻力时的伯努利方程

𝐻𝑓= 𝜁𝑢122=(𝑝1𝜌+𝑢122)−(𝑝2𝜌+𝑢222)

可以得到局部阻力系数的计算式:

𝜁=1−𝑢22+2?𝑝/𝜌𝑢12

式中,𝑢1、𝑢2分别为细管和粗管中的平均流速,?𝑝为2,1截面的压差。

突然扩大管的理论计算式为:ζ=(1−𝐴1/𝐴2)2 ,𝐴1、𝐴2分别为细管和粗管的流通截面积。

三、实验流程: 本实验装置如图1所示,管道水平安装,水循环使用,其中管5为不锈钢管,测压点之间距1.50m,内径21.0mm;管6为镀锌钢管,测压点间距离1.50m,内径22..5mm;管7为突然扩大管,由16.0mm扩大至42.0mm。各测量元件由测压口与压差传感器相连,通过管口的球阀切换被测管路,系统流量由涡轮流量计3调节,离心泵的功率由变频器通过改变输入频率控制转速来实现控制。

四、实验操作要点:

1.开泵:在关闭所有阀门的情况下,打开电源,启动变频器至50Hz,固定转速,观察泵出口压力稳定后,即可进行排气。

2.排气:在对某一管路进行实验之前,排尽设备主管和该管路及对应测压管路内的空气,每切换管路都要排一次气。关闭其他控制阀,打开对应管路的控制阀、测压阀和排气阀,在50Hz下,调节流量至1-2𝑚3/ℎ,待2min以上,压差传感器示数稳定后,关闭排气阀和流量调节阀,在流量为0下观察压差传感器示数是否为0,若有较大偏差则气未排尽,若偏差较小且稳定则记录初始偏差值。

3.实验数据测取:确定排气完毕且其余管路切换阀和测压阀关闭后,调节变频器至25Hz左右。对于直管阻力,按照流量由大到小的顺序,测取10组数据,控制压差在3.7~0.7kpa之间。对于突然扩大管的阻力,可测取3组数据。测取数据时,每个数据点取值应等待2min以上且压差和流量稳定为某值或在很小范围内波动。波动时可取其中点。

五、原始数据及处理:

1.原始数据记录

水的物理性质:

测定光滑管时,25℃下,ρ=996.95kg/m3,μ=0.8973mPa/s

测定粗糙管及突然扩大管时,26.5℃下,ρ=996.58g/m3,μ=0.8641mPa/s

(1)光滑管和粗糙管实验数据

光滑管数据:不锈钢管,l=1.500m,d=21.0mm,ε≈0.02mm,零点误差p0=-0.04kPa。 粗糙管数据:镀锌钢管,l=1.500m,d=22.5mm,ε≈0.2mm,零点误差p0=-0.06kPa。

表1 光滑管和粗糙管原始数据记录表

序号 光滑管 粗糙管

流量/(m3*h-1) 压差/(kPa) 流量/(m3*h-1) 压差/(kPa)

1 3.69 7.14 3.69 6.61

2 3.34 5.99 3.36 5.56

3 3.1 5.2 3.06 4.65

4 2.7 4.02 2.75 3.75

5 2.45 3.39 2.44 3.02

6 2.08 2.5 2.02 2.13

7 1.71 1.75 1.73 1.6

8 1.36 1.14 1.37 1.05

9 1.03 0.7 1.01 0.59

10 0.69 0.32 0.69 0.3

(2)突然扩大局部阻力系数测定数据

突扩管: d1=16.0mm,d2=42.0mm,初始误差p0=-0.04kPa。

表2 突然扩大局部阻力系数数据记录表

序号 流量/(m3*h-1) 压差/(kPa)

1 3.57

3.16

2 2.51 1.54

3 1.46 0.49

2.数据处理

表3 光滑管数据处理表

序号 流量/(m3*h-1) 流速/(m*s-1) 实际压差/(kPa) Re λ λb

1 3.69 2.96 7.18 69034 0.02303 0.01951

2 3.34 2.68 6.03 62486 0.02360 0.02001

3 3.1 2.49 5.24 57996 0.02381 0.02038

4 2.7 2.17 4.06 50513 0.02432 0.02110

5 2.45 1.96 3.43 45836 0.02495 0.02162

6 2.08 1.67 2.54 38913 0.02564 0.02252

7 1.71 1.37 1.79 31991 0.02673 0.02365

8 1.36 1.09 1.18 25443 0.02786 0.02504

9 1.03 0.83 0.74 19270 0.03046 0.02685

10 0.69 0.55 0.36 12909 0.03302 0.02967

其中,λb项为根据Blasius公式计算的理论摩擦系数值。

直管阻力系数的计算示例:

由表3中第1组数据为例, u=𝑞𝑣𝐴=4𝑞𝑣𝜋𝑑2=4×3.693.14×(21.0×10−3)2×13600𝑚/𝑠=2.96𝑚/𝑠

𝑅𝑒=𝑑𝑢𝜌𝜇=(21.0×10−3)×2.96×996.950.8973×10−3=69034

𝜆=2?𝑝𝑑𝜌𝑢2𝑙=2×(7.18×103)×(21.0×10−3)996.95×2.962×1.5=0.02303

𝜆𝑏=0.3163𝑅𝑒0.25=0.3163690340.25=0.01951

表4 粗糙管数据处理表

序号 流量/(m3*h-1) 流速/(m*s-1) 实际压差/(kPa) Re λ

1 3.69 2.58 6.67 66896 0.03021

2 3.36 2.35 5.62 60913 0.03070

3 3.06 2.14 4.71 55474 0.03102

4 2.75 1.92 3.81 49854 0.03107

5 2.44 1.70 3.08 44234 0.03191

6 2.02 1.41 2.19 36620 0.03310

7 1.73 1.21 1.66 31363 0.03421

8 1.37 0.96 1.11 24837 0.03648

9 1.01 0.71 0.65 18310 0.03930

10 0.69 0.48 0.36 12509 0.04664

图2 光滑管和粗糙管的λ-Re关系曲线

曲线分析:

(a)光滑管和粗糙管的摩擦系数均随Re的增大而减小,且随着Re的增大,摩擦系数减小的趋势趋缓。

(b)在同一Re下,相对粗糙度更高的粗糙管比光滑管的摩擦系数更大,说明ε/d越大,摩擦系数越大。

(c)在同一Re下,光滑管的摩擦系数大于水力学光滑摩擦系数的理论值,说明实验用的光滑管和理论光滑有一定差距。

表5 突然扩大管数据处理表

序号 流量/(m3*h-1) 压差/(kPa) 细管流速/(m*s-1) 粗管流速/(m*s-1) ζ

1 3.57 3.20 4.93 0.72 0.7149

2 2.51 1.58 3.47 0.50 0.7152

3 1.46 0.53 2.02 0.29 0.7175

局部阻力的计算示例: 以表5中第1组数据为例,

𝑢1=𝑞𝑣𝐴1=4𝑞𝑣𝜋𝑑12=4×3.573.14×(16.0×10−3)2×13600𝑚/𝑠=4.93𝑚/𝑠

𝑢2=𝑞𝑣𝐴2=4𝑞𝑣𝜋𝑑22=4×3.573.14×(42.0×10−3)2×13600𝑚/𝑠=0.72𝑚/𝑠

𝜁=1−𝑢22+2?𝑝𝜌𝑢12=1−0.722+2×3.20×103996.584.932=0.7149

𝜁=∑𝜁𝑖3=0.7159

理论值ζ𝑡=(1−𝐴1/𝐴2)2=(1−𝑑12/𝑑22)2=(1−162/422)2=0.7308

相对偏差

δ=|𝜁−ζ𝑡ζ𝑡|×100%=|0.7159−0.73080.7308|×100%=2.04%

测量值与理论值基本符合,但存在一定误差。

五、结果讨论分析

1.本次曲线拟合的相对大小比较准确,但是其中表现的趋势不明显,并未得到随着雷诺数增大,摩擦系数趋近于某一值的结论。可能是测定的摩擦系数和雷诺数范围较小,如果增大测定的雷诺数上限,即在更高的流速下做实验,可以看到更好的趋势。

2.测定的局部阻力系数和理论值接近,说明实验结果较好。实验值低于理论值,可能是实验设备本身存在损耗,细管在高流量下腐蚀变粗的结果。可以看到随着流量增大有上升趋势,而的三次结果的差值应该是被忽略的直管阻力的影响,因而随着流量增大,表观的局部阻力系数应该增大而不是减小,可能是实验记录和计算舍入的影响。