沿程损失阻力系数的FLUENT数值模拟计算报告 李济然.docx
- 文档编号:27900777
- 上传时间:2023-07-06
- 格式:DOCX
- 页数:10
- 大小:87.73KB
沿程损失阻力系数的FLUENT数值模拟计算报告 李济然.docx
《沿程损失阻力系数的FLUENT数值模拟计算报告 李济然.docx》由会员分享,可在线阅读,更多相关《沿程损失阻力系数的FLUENT数值模拟计算报告 李济然.docx(10页珍藏版)》请在冰豆网上搜索。
沿程损失阻力系数的FLUENT数值模拟计算报告李济然
沿程损失阻力系数
工程力学2007级李济然20071210114
一概述:
沿程损失水流流动过程中,由于固体壁面的阻滞作用而引起的摩擦阻力所造成的水头损失。
流体流动中为克服摩擦阻力而损耗的能量称为沿程损失。
沿程阻力损失与长度、粗糙度及流速的平方成正比,而与管径成反比,沿程能量损失的计算公式是:
hr=λv2/(2dg)
其中:
为管长,
为沿程损失系数,
为管道内径,v2/(2g)为单位重力流体的动压头(速度水头),
为流体的运动粘度系数。
沿程损失能量损失的计算公式由带粘性的伯努利方程:
v12/(2g)+p1/(ρg)+z1=v22/(2g)+p2/(ρg)+z2+hf
推出:
hf=(p1-p2)/(ρg)
其中:
v22/(2g)——单位质量流体的动能(速度水头)。
流体静止时为0。
Z——单位质量流体的势能(位置水头)。
p/(ρg)——单位质量流体的压力能(压强水头)。
又由量纲分析的
定理,得出Δp/(ρV2/2)=λL/d,计算出达西摩擦因子λ=2Δpd/(LρV2),则hf=λLV2/(2gD)
由于Re=Vd/v和v=μ/ρ,则λ=f(Red)
湍流光滑管的沿程损失系数按卡门一普朗特(Karmn-Prandtl)公式:
1/λ1/2=2lg(Reλ1/2)-0.8
当105<Re<3×106时,尼古拉兹的计算公式为:
λ=0.0032+0.221Re-0.237
1.湍流粗糙管过渡区:
26.98(d/ε)8/7<Re<2308(d/ε)0.85为湍流粗糙管过渡区。
该区域的沿程损失系数与按洛巴耶夫(Б.H.Лo6aeв)的公式进行计算,即
λ=1.42[lg(dRe/ε)]-2=1.42[lg(1.273qv/vε)]-2
2.湍流粗糙管平方阻力区:
2308(d/ε)0.85 沿程损失系数与雷诺数无关,只与相对粗糙度有关。 平方阻力区的沿程能量损失可按尼占拉兹公式: 1/λ1/2=2lg(d/2ε)+1.74 关于沿程损失最著名的是尼古拉茨在1932~1933年问所做的实验(右图为实验装置图)。 其测得曲线,从此得出了几个重要结论: 1.层流区: Re<2320为层流区。 在该区域内,管壁的相对粗糙度对沿程损失系数没有影响。 2.过渡区: 2320<Re<4000为由层流向湍流的转换区,可能是层流,也可能是湍流,实验数据分散,无一定规律,如图中的区域所示。 3.湍流光滑管区: 4000<Re<26,98(d/ε)8/7,为湍流光滑管区。 勃拉修斯(p.Blasius)1911年用解析方法证明了该区沿程损失系数与相对粗糙度无关,只与雷诺数有关,并借助量纲分析得出了4×10e3<Re<10e5范围内的勃拉休斯的计算公式为 λ=0.3164/Re0.25 二实验目的: 1、测量管道两端压差ΔP 2、应用粘性的伯努利方程,求出管道的沿程能量损失hf 3、根椐hf=λlv2/(2gd),求出λ值,并绘出λ-f(Red)曲线。 三实验原理: 由连通器原理,在管道上安装一根U形管,U形管两端与测压孔连通,U形管中间有一阀门3与大气相通(见装置图),将要测量之溶液注入U形管内,首先打开阀门3,用泵将管道中空气排出后,然后关闭阀门3。 U形管内溶液的密度 为已知,当管道两端测压孔处存在压差时,U形管中的二液面将有高差 ,由此,根据水静力学基本方程 ,由测出的高差 ,可求出管道两端测压孔处的压差 。 四实验装置图 五实验过程 1、打开水伐,将水箱放满水,打开管路中的阀门一和阀门二,然后再打开泵的电源开关,检查整个实验装置水路是否畅通。 2、当实验装置水路稳定后,再利用阀门三将U形管中的气体排出,随后观察U形管中之液面的变化。 3、计录U形管中两端的水压差高度 。 并记入记录表格。 测读时应尽量减少误差。 4、记录表格 λ1为布拉休斯解 Q(m3/h) Re v △h ΔP λ λ1 误差 4.14 29023 0.8301 0.045 440.56 0.0179 0.0242 26.07% 4.26 29864 0.8541 0.041 401.4 0.0154 0.0241 35.93% 4.53 31757 0.9083 0.045 440.56 0.015 0.0237 36.85% 4.76 33369 0.9544 0.049 479.72 0.0148 0.0234 36.94% 5.08 35613 1.0185 0.060 587.41 0.0159 0.0230 31.09% 5.40 37856 1.0827 0.066 646.15 0.0154 0.0227 31.89% 5.56 38978 1.1148 0.069 675.52 0.0152 0.0225 32.34% 5.85 41011 1.1729 0.175 1713.3 0.0349 0.0222 -56.99% 6.12 42904 1.2270 0.177 1732.9 0.0323 0.022 -46.73% 管径d=42mm管长L=3000mmμ=0.0012kg/m*sρ=999kg/m3 六计算模型: 研究湍流水力光滑区的达西摩擦因子与Re的关系。 在FLUENT中通过改变流速或者粘度系数控制Re,并进行数值模拟,从而计算出管中试验段两端的压力的差值,即可得到沿程损失阻力系数,再将所得的值与水力光滑区曲线对比,判断其是否正确。 模型: 建立一个半径r=21mm,长l=5m的圆截面直管,其中前2m是前置段,用来让湍流充分发展,后3m为实验段。 假设其材料是光滑的,没有摩擦,内部流体为水。 设水的ρ为1000 ,粘度系数μ为0.001 。 左图就是试算的velocityinlet后端 云图,说明在试验段之前设置前置段还是十分有必要的。 由于液体的粘性力作用,在壁面附近有比较大的速度梯度,而且在入口端是湍流发展段,所以需要端面使用边界层网格加密,轴向在入口处加密。 具体步骤是: 1.做半径为0.22的圆。 2.做出x=0.022,y=0,z=5的点,并连接圆上与其对应的两点。 3.为该线linemesh,选择ratio1.04,让网格在圆附近加密。 在此同时将将端面的圆分成60等分线网格。 4.使用sweep命令,选上withmesh选项,让直线绕圆周旋转成圆柱面,并且将网格自动画好。 5.端面上创建边界层网格,firstpercentage在这里取了14,rows取5层。 6.为端面直接画面网格,由于之前端面的圆已经分好了网格和边界层网格,不用设定参数gambit自动画网格,完成后如下图。 7.在生成体的选项中选择sweep,勾选withmesh选项,让圆端面沿管轴线方向扫过,即可完成体网格的绘制。 8.最后选择求解器(solver)Fluent5/6,设置入口、壁面和出口。 注意: 在画边界层网格时有个方向选择问题,打开edge的list里面,每个edge其实可以点多次,具体多少次看该edge属于多少个face,通过试验,就可以看到边界层具体会向哪个方向生成。 七数值模拟及数据处理 由于光滑圆管,则达西摩擦因子λ只是Re的函数。 Re<2000时,圆管中的流动属于层流,在泊肃叶的实验中,经验公式为: λ=64/Red Red定义为dv/v,在这里V为距入口10m(即试验段的起始端)的截面平均流速,湍流时的Re也如此定义。 在圆管流动中雷诺数Re>2000时才进入湍流状态,并且在2000 为了更好的与尼古拉茨试验的比对,选择Red=3*103~105内的10个值3500、4000、4500、5000、6000、7000、9000、12000、15000、20000作为入口的Re,具体的Red需模拟后才能得出,再将这几个数值作出曲线和误差分析。 使用ANSYS12.0中的Fluent作为流场模拟的软件,在这里圆管属于细长结构中的流动用双精度模式模拟精确一些。 准备使用k-epsilon和S-A湍流模式分别计算。 Rein为入口雷诺数,v为入口速度,p1是试验段起始端的压力,p2是试验段结束端的压力。 Red为实验段起始处雷诺数。 λ1和λ2分别为模拟算出的达西摩擦因子和用布拉休斯公式算出的达西摩擦因子。 标准k-epsilon在壁面区使用了不够精确的近壁函数的半经验公式,以及工况中流场为层流向湍流的过度区都有可能导致误差过大。 将FLUENT中对K-epsilon做修改(见附录二)。 再通过此表数据作出拟合曲线与布拉休斯公式的解对比,分析误差。 增强壁面函数的K-epsilon湍流模式计算结果 ν p1 p2 误差 3500 0.0833 -10.265813 -14.505771 0.051313 0.041174 24.76% 4000 0.0952 -12.627162 -17.870213 0.048533 0.039734 22.14% 4500 0.1071 -15.308632 -21.644656 0.046434 0.038687 20.10% 5000 0.119 -18.074271 -25.591423 0.044594 0.037623 18.50% 6000 0.1428 -24.361341 -34.508725 0.041824 0.035345 16.26% 7000 0.1666 -31.505652 -44.656134 0.039753 0.034524 15.05% 9000 0.2142 -48.221476 -68.535343 0.037123 0.032456 14.48% 12000 0.2856 -79.798474 -113.648412 0.034532 0.030566 15.30% 15000 0.357 -119.07143 -170.06142 0.033654 0.028354 17.54% 20000 0.476 -190.21973 -271.87977 0.030234 0.026676 13.78% 对比标准k-epsilon的精度高多了,但是仍然不够精确。 如果将网格划分得更精细些,将更好的控制误差。 S-A湍流模式计算 ν p1 p2 误差 3500 0.0833 -8.285567 -18.98961 0.043193 0.041136 5.00% 4000 0.0952 -10.27976 -23.63747 0.041268 0.039785 3.73% 4500 0.1071 -12.46471 -28.73887 0.039726 0.038631 2.84% 5000 0.119 -14.83306 -34.28072 0.038453 0.037627 2.20% 6000 0.1428 -20.10822 -46.64936 0.036444 0.03595 1.37% 7000 0.1666 -26.07841 -60.67434 0.034901 0.034591 0.90% 9000 0.2142 -40.20676 -93.95073 0.032798 0.032484 0.97% 12000 0.2856 -65.64599 -153.9034 0.030297 0.03023 0.22% 15000 0.357 -96.69733 -227.0509 0.028638 0.02859 0.17% 20000 0.476 -162.0646 -381.9362 0.027171 0.026606 2.13% 以上是两种湍流模式模拟的曲线图,很明显S-A模式模拟出的结果优于k-epslion模式的。 附录一 [参考文献] [1]周光垌,严宗毅,许世雄,章克本.流体力学[M].北京;高等教育出版社,2000. [2]周志军,林震,周俊虎,刘建忠,岑可法.不同湍流模型在管道流动阻力计算中的应用和比较[J].热力发电。 2007。 [3]宁方飞,徐力平.Spalart—Allmaras湍流模型在内流流场数值模拟中的应用[J].工程热物理学报。 2001。 附录二
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 沿程损失阻力系数的FLUENT数值模拟计算报告 李济然 损失 阻力 系数 FLUENT 数值 模拟 计算 报告