沿程损失阻力系数的FLUENT数值模拟(计算流体力学作业)

1、计算流体力学课程作业作业题目：沿程损失阻力系数的FLUENT数值模拟学生姓名：易鹏学生学号：专业年级：动力工程及工程热物理12级学院名称：机械与运载工程学院 2012年5月2日沿程损失阻力系数的FLUENT数值模拟一、 引言沿程损失（pipeline friction loss）是指管道内径不变的情况下，管内流体流过一段距离后的水头损失。其中边界对水流的阻力是产生水头损失的外因，液体的粘滞性是产生水头损失的内因，也是根本原因。沿程能量损失的计算公式是：。其中：为管长，为沿程损失系数，为管道内径，为单位重力流体的动压头（速度水头），为流体的运动粘度系数。粘性流体在管道中流动时，呈现出两种流动状态

2、，管道中的流速(为层流向湍流转变的临界流速)为层流，此时整个流场呈一簇互相平行的流线。则时为湍流，流场中的流体质点作复杂的无规则的运动。沿程损失与流动状态有关，故计算各种流体通道的沿程损失，必须首先判别流体的流动状态。沿程损失能量损失的计算公式由带粘性的伯努利方程 推出，可知， 其中： 单位质量流体的动能（速度水头）。流体静止时为0。 单位质量流体的势能（位置水头）。单位质量流体的压力能（压强水头）。又由量纲分析的定理，得出 ，计算出达西摩擦因子， 则,由于，则。关于沿程损失最著名的是尼古拉茨在1932 1933年问所做的实验(右图为实验装置图)。其测得曲线如图1，从此得出了几个重要结论：1层

3、流区Re2320为层流区。在该区域内，管壁的相对粗糙度对沿程损失系数没有影响。2过渡区2320Re4000为由层流向湍流的转换区，可能是层流，也可能是湍流，实验数据分散，无一定规律。3湍流光滑管区4000Re2698（d/）8/7，为湍流光滑管区。勃拉修斯（p.Blasius）1911年用解析方法证明了该区沿程损失系数与相对粗糙度无关，只与雷诺数有关，并借助量纲分析得出了410e3Re10e5范围内的勃拉休斯的计算公式为湍流光滑管的沿程损失系数也可按卡门一普朗特（Karmn-Prandtl）公式进行计算。当105Re3106时，尼古拉兹的计算公式为4湍流粗糙管过渡区2698（d/）8/7Re2

4、308（d/）0.85为湍流粗糙管过渡区。该区域的沿程损失系数与按洛巴耶夫（.H.o6ae）的公式进行计算，即5湍流粗糙管平方阻力区2308(d/)0.85Re为湍流粗糙管平方阻力区。沿程损失系数与雷诺数无关，只与相对粗糙度有关。平方阻力区的沿程能量损失可按尼占拉兹公式 进行计算。图1尼古拉茨曲线三、数值模拟 1、前处理因为层流有精确解所以在此不做讨论，而湍流状态下如果考虑圆管的粗糙度也是十分复杂，而且在粗糙表面的流动很难模拟，所以在此我们重点研究湍流水力光滑区的达西摩擦因子与Re的关系。FLUENT中通过改变流速或者粘度系数来控制Re，并进行数值模拟，计算出管中试验段两端的压力的差值，即可得

5、到沿程损失阻力系数（达西摩擦因子），再将所得的值与上图水力光滑区曲线或布拉休斯公式对比，判断其是否正确。 1模型建立一个半径r=21mm，长l=3m的圆截面直管，其中前2m是前置段，用来让湍流充分发展，后1m为实验段。假设其材料是光滑的，没有摩擦，内部流体为水。设水的为1000，粘度系数为0.001。下图就是使用K-epsilon湍流模式试算的velocity inlet后端云图，说明在试验段之前设置前置段还是十分有必要的。使用gambit可以很容易的建模，直接使用cylinder命令建立方案中的模型。但是在此未使用这种方法，由于液体的粘性力作用，在壁面附近有比较大的速度梯度，而且在入口端是湍

6、流发展段，所以需要端面使用边界层网格加密，轴向在入口处加密。具体步骤是：1做半径为0.021的圆。2做出x=0.021，y=0，z=5的点，并连接圆上与其对应的两点。3为该线mesh，选择ratio 1.05，让线网格在入口处加密。在此同时将将入口端面的圆分成50等分线网格（数目自定，但是这样已经足够）4使用sweep命令，选上with mesh选项，让直线绕圆周旋转成圆柱面，并且将网格自动画好。如右图。5.端面上创建边界层网格，first percentage(第一层边界层网格的高度关于宽度的百分比)在这里取了15,rows取5层，Growth factor取1.1。(注意一个问题，就是在画

7、边界层网格时有个方向选择问题，打开edge的list里面，每个edge其实可以点多次，具体多少次看该edge属于多少个face，通过试验，就可以看到边界层具体会向哪个方向生成)。具体设置如右图。6.为端面直接画面网格，由于之前 端面的圆已经分好了网格和边界层网格，不用设定参数gambit自动画网格，完成后如下图。7在生成体的选项中选择sweep，勾选with mesh选项，让圆端面沿管轴线方向扫过，即可完成体网格的绘制。8最后选择求解器（solver）Fluent 5/6，设置z=0处端面为Velocity in、圆柱面为wall和z=3处端面为outflow。9Export mesh。注意：

8、不要选择2D模型输出的选项。下图是网格完成后的模型。一共生成了50800个体网格。四、数值模拟及数据处理由于是光滑圆管（或水力光滑），则达西摩擦因子只是Re的函数。而在Re2000才进入湍流状态，并且在2000Re4000时，为层流向湍流过渡区。为了更好的与尼古拉茨试验的比对，选择内的10个值3500、4000、4500、5000、6000、7000、9000、12000、15000、20000作为入口的Re，具体的需模拟后才能得出，再将这几个数值作出曲线和误差分析。使用ANSYS12.0中的Fluent作为流场模拟的软件，在这里圆管属于细长结构中的流动用双精度（Double Precisio

9、n）模式模拟较精确。准备使用k-epsilon，增强壁面函数的k-epsilon和S-A湍流模式分别计算。而且由于流动是湍流，并且网格在内部并不是和流速垂直的，所以使用二阶迎风格式，来提高精度，并且设置残差到10e-5，以提高精度。通过Fluent的Report菜单中的surface Integrals命令可以获得入口和出口的压力和速度的平均值。下图为Surface Integrals的窗口，其中Inexp是实验段的入口截面，out就是出截面：为入口雷诺数，v为入口速度，p1是试验段起始端的压力，p2是试验段结束端的压力。为实验段起始处雷诺数。和分别为模拟算出的达西摩擦因子和用布拉休斯公式算出

10、的达西摩擦因子。再通过此表数据作出拟合曲线与布拉休斯公式的解对比，分析误差。K-epsilon湍流模式计算结果p1p2误差35000.0833-22.77689-48.1338420.10230.0411148.74%40000.0952-27.525547-58.7379340.096430.0398142.38%45000.1071-32.735207-70.200340.09150.0386136.74%50000.1190-38.257824-82.3440090.08720.0376131.67%60000.1428-50.02042-108.35730.08010.0359122.

11、82%70000.1666-62.615845-136.41620.07450.0346115.23%90000.2142-90.401306-198.274140.06580.0325102.65%120000.2856-136.95033-302.464230.05680.030287.95%150000.3570-188.45871-417.953130.05040.028676.35%200000.4760-282.94147-630.133730.04290.026661.26%可见误差相当之大，究其原因，应该是标准k-epsilon在壁面区使用了不够精确的近壁函数的半经验公式，以及

12、工况中流场为层流向湍流的过度区。在FLUENT中对K-epsilon做如下修改：增强壁面函数的K-epsilon湍流模式计算结果p1p2误差35000.0833-10.265804-14.5057140.0513270.04113624.77%40000.0952-12.627173-17.8702810.0485950.03978522.15%45000.1071-15.308632-21.6446250.04640.03863120.11%50000.119-18.074289-25.5914860.044590.03762718.51%60000.1428-24.361353-34.50

13、87090.04180.0359516.27%70000.1666-31.505663-44.6561430.0397990.03459115.06%90000.2142-48.221485-68.5353010.037190.03248414.49%120000.2856-79.798485-113.648510.034860.0302315.31%150000.357-119.07154-170.061510.0336070.0285917.55%200000.476-190.21982-271.879880.0302740.02660613.79%对比标准k-epsilon的精度高多了，

14、但是仍然不够精确。如果将网格划分得更精细些，将更好的控制误差。 S-A湍流模式计算p1p2误差35000.0833-8.285567-18.989610.0431930.0411365.00%40000.0952-10.27976-23.637470.0412680.0397853.73%45000.1071-12.46471-28.738870.0397260.0386312.84%50000.119-14.83306-34.280720.0384530.0376272.20%60000.1428-20.10822-46.649360.0364440.035951.37%70000.1666

15、-26.07841-60.674340.0349010.0345910.90%90000.2142-40.20676-93.950730.0327980.0324840.97%120000.2856-65.64599-153.90340.0302970.030230.22%150000.357-96.69733-227.05090.0286380.028590.17%200000.476-162.0646-381.93620.0271710.0266062.13%以上是两种湍流模式模拟的曲线图，很明显S-A模式模拟出的结果优于k-epslion模式的。 五、总结 虽然模拟的结果和经验公式还是有误差，但是经验公式本身也是不精确的，而且由于湍流模式和数值计算中都有不可避免的误差，然而能做到5%以内已经满足工程上的需要了。而为什么S-A这种适用于低雷诺数的湍流模式，反而在高雷诺数情况下的计算结果更精于k-，通过查找文献，有以下的解释，一个湍流模型要想精确地求解出流动阻力，必须考虑近壁区低雷诺数的