沿程损失阻力系数的FLUENT数值模拟计算流体力学作业.docx
- 文档编号:7872548
- 上传时间:2023-01-26
- 格式:DOCX
- 页数:14
- 大小:332.86KB
沿程损失阻力系数的FLUENT数值模拟计算流体力学作业.docx
《沿程损失阻力系数的FLUENT数值模拟计算流体力学作业.docx》由会员分享,可在线阅读,更多相关《沿程损失阻力系数的FLUENT数值模拟计算流体力学作业.docx(14页珍藏版)》请在冰豆网上搜索。
沿程损失阻力系数的FLUENT数值模拟计算流体力学作业
计算流体力学课程作业
沿程损失阻力系数的
FLUENT数值模拟作业题目:
易鹏:
名生姓学
:
生学号学
12级动力工程及工程热物理:
级业年专
机械与运载工程学院:
名称学院
日2月5年2012.
沿程损失阻力系数的
FLUENT数值模拟
一、引言
沿程损失(pipelinefrictionloss)是指管道内径不变的情况下,管内流体流过一段距离后的水头损失。
其中边界对水流的阻力是产生水头损失的外因,液体的粘滞性是产生水头损失的内因,也是根2vl本原因。
沿程能量损失的计算公式是:
。
其中:
为管长,=λh?
l
fd2g2v为沿程损失系数,为管道内径,为单位重力流体的动压头(速度d
2g水头),为流体的运动粘度系数。
粘性流体在管道中流动时,呈现出v两种流动状态,管道中的流速(为层流向湍流转变的临界流速)vvv?
crcr为层流,此时整个流场呈一簇互相平行的流线。
则时为湍流,vv?
cr流场中的流体质点作复杂的无规则的运动。
沿程损失与流动状态有关,故计算各种流体通道的沿程损失,必须首先判别流体的流动状态。
沿程损失能量损失的计算公式由带粘性的伯努利方程
22ppvvP-P2211hz++++z=+21=h推出,可知,其中:
f212gρg2gρg
fρg2v——单位质量流体的动能(速度水头)。
流体静止时为0。
2g。
——单位质量流体的势能(位置水头)zp。
——单位质量流体的压力能(压强水头)
gρ1
ΔpL=λ,计算出达西摩擦因子又由量纲分析的定理,得出?
1d2Vρ
22LVμVd2Δpdh=λν==Re=λ,,由于。
,则,则)=f(Reλ
f
22gDρLρVνd关于沿程损失最著名的
1933~是尼古拉茨在1932右图为实(年问所做的实验其测得曲线如图。
验装置图),从此得出了几个重要结1论:
为层流区。
在该区域内,管壁的相对粗糙<2320.层流区Re1
度对沿程损失系数没有影响。
为由层流向湍流的转换区,可能是4000Re2320<<2.过渡区
层流,也可能是湍流,实验数据分散,无一定规律。
,为湍流光滑ε)8/79826.(d/Re.湍流光滑管区34000<<年用解析方法证明了该区沿程损1911)管区。
勃拉修斯(p.Blasius并借助量纲分析得出了只与雷诺数有关,失系数与相对粗糙度无关,范围内的勃拉休斯的计算公式为10e510e3×<Re<40.3164?
?
0.25Re2
湍流光滑管的沿程损失系数也可按卡门一普朗特(Karmn-Prandtl)公式
11/2?
)?
2lg(Re?
0.8
1/2?
进行计算。
当105<Re<3×106时,尼古拉兹的计算公式为
?
?
0.0032?
0.221Re?
0.237
4.湍流粗糙管过渡区26.98(d/ε)8/7<Re<2308(d/ε)0.85为湍流粗糙管过渡区。
该区域的沿程损失系数与按洛巴耶夫(Б.H.Лo6aeв)的公式进行计算,即
?
2?
?
?
?
qd?
?
?
?
?
1.42?
V?
?
1.273Re?
1.42lglg?
?
?
?
?
?
?
?
v?
?
?
?
?
?
?
?
?
?
?
?
5.湍流粗糙管平方阻力区2308(d/ε)0.85 沿程损失系数与雷诺数无关,只与相对粗糙度有关。 平方阻力区的沿程能量损失可按尼占拉兹公式 1d? 1.74? 2lg ? 1/22? 进行计算。 3 尼古拉茨曲线图1三、数值模拟、前处理1 因为层流有精确解所以在此不做讨论,而湍流状态下如果考虑圆所以在管的粗糙度也是十分复杂,而且在粗糙表面的流动很难模拟,的关系。 此我们重点研究湍流水力光滑区的达西摩擦因子与Re,并进行数值模拟,中通过改变流速或者粘度系数来控制ReFLUENT计算出管中试验段两端的压力的差值,即可得到沿程损失阻力系数,再将所得的值与上图水力光滑区曲线或布拉休斯(达西摩擦因子)公式对比,判断其是否正确。 的圆截面直管,l=3m建立一个半径r=21mm,长1.模型为实验段。 假设1m2m其中前是前置段,用来让湍流充分发展,后4 其材料是光滑的,没有摩擦,内部流体为水。 设水的ρ为3kg/(m*s)kg/m。 下图就是使用1000为,粘度系数μ0.001v云图,说明在velocityinlet后端K-epsilon湍流模式试算的y试验段之前设置前置段还是十分有必要的。 命令建立可以很容易的建模,直接使用cylinder使用gambit由于液体的粘性力作用,但是在此未使用这种方法,方案中的模型。 在壁面附近有比较大的速度梯度,而且在入口端是湍流发展段,所具体步骤是: 以需要端面使用边界层网格加密,轴向在入口处加密。 的圆。 .做半径为0.0211的点,并连接圆上与其对应的两,z=5x=0.0212.做出,y=0点。 ratiomesh,选择3.为该线,让线网格在入口处加密。 1.05在此同时将将入口端面的圆分等分线网格(数目自定,50成但是这样已经足够)选项,让直线绕圆周旋命令,选上.使用4sweepwithmesh5 如转成圆柱面,并且将网格自动画好。 右图。 first5.端面上创建边界层网格,第一层边界层网格的高度percentage(了取在这里度的百分比)关于宽。 1.1factor取15,rows取5层,Growth 注意一个问题,就是在画边界层网格(的时有个方向选择问题,打开edge其实可以点多次,edge里面,list每个个多少该edge属于具体多少次看就可以看到边界层具通过试验,face,。 具体设置如右体会向哪个方向生成)图。 为端面直接画面网格,由于之6.端面的圆已经分好了网格和边界前 自动画层网格,不用设定参数gambit网格,完成后如下图。 ,sweep7.在生成体的选项中选择选项,让圆端面沿管轴线方向扫过,即可完成体withmesh勾选网格的绘制。 处端面z=0,设置solver8.最后选择求解器()Fluent5/6outflowz=3wallVelocityin为、圆柱面为和处端面为。 6 9.Exportmesh。 注意: 不要选择2D模型输出的选项。 下图是网格完成后的模型。 一共生成了50800个体网格。 四、数值模拟及数据处理的Re由于是光滑圆管(或水力光滑),则达西摩擦因子λ只是时,圆管中的流动属于层流,泊肃叶也做过此范Re<2000函数。 而在64dV在这里定义为围内流动的大量实验,得出经验公式,,=λRed νRedRe10m(即试验段的起始端)的截面平均流速,湍流时的V为距入口才进入湍流状态,也如此定义。 又因为在圆管流动中雷诺数Re>2000时,为层流向湍流过渡区。 为了更好的与尼古2000 ? Re、、4500个值拉茨试验的比对,选择3500、4000内的103~10d,具Re、1500020000作为入口的、、、、、500060007000900012000体的需模拟后才能得出,再将这几个数值作出曲线和误差分析。 Red7 使用ANSYS12.0中的Fluent作为流场模拟的软件,在这里圆管属于细长结构中的流动用双精度(DoublePrecision)模式模拟较精确。 准备使用k-epsilon,增强壁面函数的k-epsilon和S-A湍流模式分别计算。 而且由于流动是湍流,并且网格在内部并不是和流速垂直的,所以使用二阶迎风格式,来提高精度,并且设置残差到10e-5,以提高精度。 通过Fluent的Report菜单中的surfaceIntegrals命令可以获得入口和出口的压力和速度的平均值。 下图为Surface Integrals的窗口,其中Inexp是实验段的入口截面,out就是出截面: 是试验段起始端的压力,p1v为入口速度,为入口雷诺数,Rein? ? 分别p2是试验段结束端的压力。 为实验段起始处雷诺数。 和Red21为模拟算出的达西摩擦因子和用布拉休斯公式算出的达西摩擦因子。 8 再通过此表数据作出拟合曲线与布拉休斯公式的解对比,分析误差。 K-epsilon湍流模式计算结果 Rein ν p1 p2 ? 1 ? 2 误差 3500 0.0833 -22.77689 -48.133842 0.1023 0.0411 148.74% 4000 0.0952 -27.525547 -58.737934 0.09643 0.0398 142.38% 4500 0.1071 -32.735207 -70.20034 0.0915 0.0386 136.74% 5000 0.1190 -38.257824 -82.344009 0.0872 0.0376 131.67% 6000 0.1428 -50.02042 -108.3573 0.0801 0.0359 122.82% 7000 0.1666 -62.615845 -136.4162 0.0745 0.0346 115.23% 9000 0.2142 -90.401306 -198.27414 0.0658 0.0325 102.65% 12000 0.2856 -136.95033 -302.46423 0.0568 0.0302 87.95% 15000 0.3570 -188.45871 -417.95313 0.0504 0.0286 76.35% 20000 0.4760 -282.94147 -630.13373 0.0429 0.0266 61.26% 可见误差相当之大,究其原因,应该是标准k-epsilon在壁面区使用了不够精确的近壁函数的半经验公式,以及工况中流场为层流向湍流的过度区。 在FLUENT中对K-epsilon做如下修改: 增9 强壁面函数的K-epsilon湍流模式计算结果 Rein ν p1 p2 ? 1 ? 2 误差 3500 0.0833 -10.265804 -14.505714 0.051327 0.041136 24.77% 4000 0.0952 -12.627173 -17.870281 0.048595 0.039785 22.15% 4500 0.1071 -15.308632 -21.644625 0.0464 0.038631 20.11% 5000 0.119 -18.074289 -25.591486 0.04459 0.037627 18.51% 6000 0.1428 -24.361353 -34.508709 0.0418 0.03595 16.27% 7000 0.1666 -31.505663 -44.656143 0.039799 0.034591 15.06% 9000 0.2142 -48.221485 -68.535301 0.03719 0.032484 14.49% 12000 0.2856 -79.798485 -113.64851 0.03486 0.03023 15.31% 15000 0.357 -119.07154 -170.06151 0.033607 0.02859 17.55% 20000 0.476 -190.21982 -271.87988 0.030274 0.026606 13.79% 对比标准k-epsilon的精度高多了,但是仍然不够精确。 如果将网格划分得更精细些,将更好的控制误差。 S-A湍流模式计算 Rein ν p1 p2 ? 1 ? 2 误差 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% 10 模式模拟出的结果优很明显S-A以上是两种湍流模式模拟的曲线图,模式的。 于k-epslion 五、总结虽然模拟的结果和经验公式还是有误差,但是经验公式本身也 是不精确的,而且由于湍流模式和数值计算中都有不可避免的误差,这种适用S-A5%以内已经满足工程上的需要了。 而为什么然而能做到k-于低雷诺数的湍流模式,反而在高雷诺数情况下的计算结果更精于一个湍流模型要想精确地求解出ε,通过查找文献,有以下的解释,特别是必须能很好地模必须考虑近壁区低雷诺数的影响,流动阻力,湍流模型可很好地Allmaras拟出近壁区的时均速度轮廓。 Spalart—湍流模式最早是用于解决—Allmaras满足上述要求,要知道Spalart湍流模型在模拟Spalart-Allmaras飞行器阻力问题,而且结果也证明εk-S-A流动阻力方面的优势,但是模式的计算成本也较高。 标准的11 模型在近壁区采用壁面函数的半经验果公式,其误差较大。 如使用修正的k-ε模式,即在近壁区的壁面函数做了修正,通过对比修正k-ε的确可以提高精度,如果使用了更精细网格会有更不错的结果。 12
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 损失 阻力 系数 FLUENT 数值 模拟 计算 流体力学 作业