三维圆管流动状况的数值模拟分析详解

1、三维圆管流动状况的数值模拟分析在工程和生活中，圆管内的流动是最常见也是最简单的一种流动，圆管流动有层流和紊流两种流动状况。层流，即液体质点作有序的线状运动，彼此互不混掺的流动；紊流，即液体质点流动的轨迹极为紊乱，质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD软件来模拟研究三维圆管的层流和紊流流动状况，主要对流速分布和压强分布作出分析。1物理模型三维圆管长l2000mm，直径d100mm。流体介质：水，其运动粘度系数1106m2/sInlet流速入，0.005m/s，2Olm/sOutlet压强出口Wall:光滑壁面，无滑移2在ICEMCFD中建立模型2.1首先建立三

2、维圆管的几何模型Geometry2.2做Blocking因为截面为圆形，故需做“O”型网格。2.3划分网格mesh注意检查网格质量。在未加密的情况下，网格质量不是很好，如下图YDynamicsmode-hitF9toreturntoselectionmode18126O0.20.30.40.50.60.80.9entitiesentitiesentitiesentities因管流存在边界层，故需对边界进行加密，网格质量有所提升，如下图Tl打u2.4生成非结构化网格，输出fluent.ms等相关文件3数值模拟原理3.1层流流动.u当水流以流速0.005m/s，从Inlet方向流入圆管，可计算出雷

3、诺数Re500,故圆管内1v流动为层流。假设水的粘性为常数（运动粘度系数1106m2/s）、不可压流体，圆管光滑，则流动的控制方程如下：t(u)X动量守恒方程(u)(uu)(uv)(uw)tXyZ(v)(vu)(vv)(vw)tXyZ(w)(wu)(wv)(ww)tXyZ质量守恒方程：(v)(w)0(1-1)yZ(u)-(-u)(-u)p(1-2)XXyyZZX(v)-(-v)(-v)-p(1-3)XXyyZZy(w)(w)(w)p(1-4)XXyyZZZ动量守恒方程：（u）t亠(uu)(uv)y(uv)(v)(vu)tX(uv)(vv)y(v2)(uw)Z(uw)Z(vw)Z(vw)_（上）

4、（上）（XXyyZU)Z-（丄）一（丄）一（丄）XXyyZZ(1-6)(1-7)式中，为密度，u、W是流速矢量在X、y和Z方向的分量，p为流体微元体上的压强。方程求解：对于细长管流，FLUENT建议选用双精度求解器，流场计算采用SIMPLE算法，属于压强修正法的一种。3.2紊流流动ud当以水流以流速0-1m/s，从Inle方向流入圆管，可计算出雷诺数Re10000，故圆管内2v流动为紊流。假设水的粘性为常数（运动粘度系数1106m2/s）、不可压流体，圆管光滑，则流动的控制方程如下：质量守恒方程：(1-5)(u)(v)(w)(wu)(wv)(ww)(xW)(xy-(丄)zz(uw)(vw)(1

5、-8)湍动能方程：(k)(ku)(kv)(kw)(x一)k)xk(y一)k)yktxyz-(zk一)k)Gzkk湍能耗散率方程()(u)(v)(w)()xk()yktxyzxLyxyzz(w2)(-JJkC2Tk(1-9)(1-10)式中，为密度，u、W是流速矢量在X、y和z方向的分量，p为流体微元体上的压强。方程求解：采用双精度求解器，定常流动，标准k模型，SIMPLEC算法。4在FLUENT中求解计算层流流动4.1导入并检查网格注意调整Scale大小。因在ICEM中作网格时，已采用的是以“米”为单位的长度，故不需更换单位。网格显示流动沿X方向，共存在283575hexahedralcel范

6、围DomainExtents:X-coordinate:min(m)=0.000000e+000,maX(m)=2.000000e+000y-coordinate:min(m)=-4.995393e-002,max(m)=4.995393e-002z-coordinate:min(m)=-4.995393e-002,max(m)=4.995393e-002GridAug20,2010FLUENT6.3(3d,pbns,lam)4.2设置求解器本模型基于压强计算，可采取绝对流速计算，Solve俅解器可采取默认设置。u雷诺数Re500，故圆管内流动为层流，Viscous!设置为Laminar。v4

7、.3定义材料因本研究采用水流动，故需使Materialty定义为Fluent设置成水。4.4设置边界条件4.4.1将Solid定义为Fluid并设置成水4.4.2定义进口InletInle定义为流速入口Velocity-in，e并设置入流速为0.005m/s4.4.3定义出Outlet出为压强出PressureOut,默认设置。4.4.4定义壁面Wall。设置为默认。4.5设置操作条件因为圆管截面较小，故可不考虑重力选项。压强选项默认为一个大气压。4.6求解方法的设置与控制4.6.1求解参数的设置、在Solutioncontr中s将Momentum设置为Secondorderupwind其他保

8、持默认。4.6.2设置监视残差注意点选Plot4.6.3流场初始化Computefrom设置为Inle。4.7监视切面4.7.1首先切取所需面以网格Grid为单位，在X方向，在0到2m之间，每隔0.2m切一平面，以来监视流速和压强的变化在Y方向，取Y=0的位置切面，相当于横剖圆柱截面；在Z方向，取Z=0的位置切面，相当于沿X轴方向竖剖圆柱截面。注意标清切面名称，以供查找。4.7.2设置监视窗口因不需监视太多所需切面,故建立4个监视窗即可，需将Plot和Write选取,设为时间步长，再Define内容。sSBBgPlotPrintWriteEveryrr厂厂monitor-3monitor-4N

9、ame|monitor-1lmonitor-2两两1717弓勻自厂厂厂厂WhenDefine.!Define.Define.Define.例监视1,监视Inle切面的流速，可设置为:厂DefineSurfaceIonitorNameVelocityJVelocityMagnitudeJSurfacese|=inletaint_solidoutletwallx-coordinate-0.2x-coordinate-0.4VjFileName4.8开始迭代设置迭代次数为200，实际比这个更少，迭代收敛时会自动停止。电Iterate|XIterationIterateApplyCloseHelpSu

10、rfaceMonitorspj5层流计算结果及分析计算120步后，已收敛，自动停止运算。残差监视窗口为0204060801e+01r1e+00-1e-01-1e-02-1e-03-1e-04-1e-05-1e-06Iterations100120Residualscontinuity-x-velocityy-velocityz-velocitvScaledResidualsAug20,2010FLUENT6.3(3d,dp,pbns,lam)5.1显示流速等值线图打开Display*Contours选择Velocit和Velocitymagnitude5.1.1入口和出口截面的流速分布图分布在

11、Surface里选择inle及outlet(1)Velocityofinlet8.48e-038.06e-037.63e-037.21e-036.79e-036.36e-035.94e-035.51e-035.09e-034.67e-034.24e-033.82e-033.39e-032.97e-032.54e-032.12e-031.70e-031.27e-038.48e-044.24e-040.00e+00ContoursofVelocityMagnitude(m/s)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)可见，入处流速分布不明显，基乎都等于入流速10.00

12、5m/s，只是外层靠近壁面处流速几乎为零，符合圆管层流流动规律，也符合边界层理论。(2)Velocityofoutlet出口截面流速分布较为明显，显同心圆分布，内层流速偏大，外层靠近壁面处流速几乎为零，边界层很薄。分层更为严重，层流显现的更为明显，且趋于稳定状态。5.1.2圆管内不同截面的流速分布图下述截面均为距inlet从0.2m到1.8m的截面（1）Velocityofinlet-0.23j36e-D37J63e-O3721e-H36.79&-036J&-035S4i-O35J51e-O35ee-(J3iJ5?e-035-0339e-032S?e-032544-032.126-031.7i

13、Dfe-01276-03S.43&-Dii2fe-04CommofVeloclVMagiltHdenA)Aig20,2010FLUENT0d,中,pb浮.lan)om&400Velocityofinlet-0.4&咏4133D6e-i3?J63e-E0721e-l36.796-035544-035j51e-(J351396-03lJ5?e-(J33J82e-E03施切2S7e-O32544-03031.7De-O31276-03&486-Qli2fe-liom&400CommofVeloclVMagiltHdenA)Aig20,2010FLUENT0d,中,pb浮.lan)Velocityof

14、inlet-0.6 )Velocityofinlet-131356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33倆阿1244-033J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400Co他册BVeloc哼Magiltide如的Aig2Q.201DFLUENT6-3lam)Velocityofinlet-0.8S.4Se-D31356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33倆阿1244-0

15、33J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400CoitalrrofVeloclVMa9kn：|如的Aig2Q.201DFLUENT6-3lam)9)Velocityofinlet-1. )Velocityofinlet-1.431356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33倆阿1244-033J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400Co他册BVeloc哼Magi

16、ltide如的Aig2Q.201DFLUENT6-3lam)6)Velocityofinlet-1.2S.4Se-D31356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33倆阿1244-033J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400CoitalrrofVeloclVMa9kn：|如的Aig2Q.201DFLUENT6-3lam)31356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33

17、倆阿1244-033J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400Co他册BVeloc哼Magiltide如的Aig2Q.201DFLUENT6-3lam)8)Velocityofinlet-1.6S.4Se-D31356-03?J63e-0?216-036.T9e-O363fe-(J353fe-E05j51e-l35j39ed33倆阿1244-033J826-E03396-032J97e-O32Sfe-O3031.7Ee-l312?e-i3S.4Se-!4OID&400CoitalrrofVeloclVMa9kn：|

18、如的Aig2Q.201DFLUENT6-3lam)8.43e-03SJCBe-037J63&-O3721e-036.79&-D363&-035J9te-03551edO35jQ9e-n34j67e4Bi2te-033J52e-03339&-032S?e-D32J5fe-03031.7D&-O3127e4B8.436-04i24e-0iOIDe-HXlCoitalnoTVeloctvM用hHide如朗Aig20,2010FLUENT6L30d.中,ptilS,lam)上述图像为圆管内部X轴方向不同截面的流速分布，可看出流速在截面上从入到出的变化。水流在圆管内部的流速分层很明显，靠近壁面处流速接近

19、于零，有一很薄的边界层，流速在边界层内很快上升，到最大流速；在圆管中央的一大片圆形区域内，流速基本一致，达到最大，且中心流速最大，为max。流速在截面的变化规律可以看出，在0到1.2m之间，每个截面的流速分布都不同，当离Inlet1.2远之后，流速在截面的分布基本一致，说明层流达到了稳定状态，这符合圆管流动进段及流中层流分布规律。以上图像因只能看到沿X轴截面的流速分布，故下面讨论从Y轴和Z轴方向看圆管的整体流速分布。5.1.3Y轴和Z轴方向流速截面截面若均沿圆管长度X方向截取，可看到对称的效果。1)Velocityofy-0整根圆管：8.48e-038.06e-037.63e-037.21e-

20、036.79e-036.36e-035.94e-035.51e-035.09e-034.67e-034.24e-033.39e-032.97e-032.12e-031.70e-038.48e-04ContoursofVelocityMagnitude(m/s)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)4.24e-040.00e+00入口段：8.48e-038.06e-037.63e-037.21e-036.79e-036.36e-035.94e-035.51e-035.09e-034.67e-034.24e-033.82e-033.39e-032.97e-032.5

21、4e-032.12e-031.70e-031.27e-038.48e-044.24e-040.00e+00ContoursofVelocityMagnitude(m/s)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)出口段：8.48e-038.06e-037.63e-037.21e-036.79e-036.36e-035.94e-035.51e-035.09e-034.67e-034.24e-033.82e-033.39e-032.97e-032.54e-032.12e-031.70e-031.27e-038.48e-044.24e-040.00e+00Contours

22、ofVelocityMagnitude(m/s)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)（2）Velocityofz-0整根圆管：8.48e-038.06e-037.63e-037.21e-036.79e-036.36e-035.94e-035.51e-035.09e-034.67e-034.24e-033.82e-033.39e-032.97e-032.54e-032.12e-031.70e-031.27e-038.48e-044.24e-040.00e+00ContoursofVelocityMagnitude(m/s)Aug20,2010FLUENT6.3(

23、3d,dp,pbns,lam)YzX以上两个截面流速分布图的效果是一样的，可以看出圆管水流入口段及之后的流速发展趋势，而且显示流速变化的规律更为明显。由数值模拟实验设置了入口均匀流速，可以认为在进口处的流速分布是均匀的，进入管内后，靠近壁面的流动受到阻滞，流速降低，形成边界层，且边界层的厚度逐渐加大，以致尚未受管壁影响的中心部分的流速加快。进口段的流动是流速分布不断变化的非均匀流动，且边界层的厚度在进口段逐渐增加，之后的流动是各个截面流速分布均相同的均匀流动，由于为层流流动，故流速分层现象很明显。但平均流速为多少？最大流速u为多少？进段长度L为多少?等等问题需要再进行讨论。max5.2轴向流速

24、的变化沿X轴截取轴线执行PloAXYPlot,选择YAxisFunctic里的Velocit和VelocityMagnitu,e选择Surface里圆管的0.20.40.60.82682对称轴line-x可得到轴向流速分布散点图。line-x8.50000000e-038.00000000e-03-7.50000000e-03-7.00000000e-03-Position(m)5.50000000e-03-4_5.00000000e-03丿4.50000000e-03Velocity50000000e-3Magnitude(m/S.00000000e-03VelocityMagnitudeA

25、ug20,2010FLUENT6.3(3d,dp,pbns,lam)由上图可以看出，在圆管的轴上，进段流速分布变化较大，从进流速0.005m/s急剧上升到最大流速u0.00848m/s。max层流入段长度有经验公式可以算的，即(1-11)L0.058dRe可算得入段长度约为1.18m，由上图显示效果可以看出，流速在离入口1.5到1.2m之间，即入段长度约为1厂1.2m，符合理论计算结果。5.3截面流速分布散点图取流动充分发展后，离Inlet1.6m远的截面x-coordinate-，其流速分布如下图（注意Plotdirection的选取），可以看处流速沿半径Y方向成抛物线分布，与理论公式抛物面

26、公式相符，即Jr24or2)(1-12)9.LILILILILILILILIe-D3Velocity1Maqnitude(rn/s)8.DDDDDDDDe-037.LlLlLlLlLlLlULle-Ll3ti.LILILILILILILILIe-LI：35.LlLlLlLlLlLlULle-Ll34.DDDDDDDDe-D33.LILILILILILILILIe-032.DDDDDDDDe-D31.LILILILILILILILIe-03LI.DDDDDDDLIe-HZlLIAig33.2010FLUENTE3fM,中,pt比ism)取沿Y方向中心轴线的流速分布，即9.DDDDDDDDe-D3

27、8.D00DDDD0e-D37.LlLlLlLlLlLlULle-D3ti.LILILILILILIULIe-035.000DD000e-D34.LlLlLlLlLlLlULle-D33.LILILILILILIULIe-D32.DDD00DDDe-031.LiLiLiLiLiLiULie-D3LI.LILILILILILlDDe-HlILIVeincityMagnitude(rn/s)Aig23,2|10FLUENTlam)5.4显示压强分布图在Contours里选取Pressur和Staticpressure在Surfaces里选择int-soli（即管道内部流体整体，以两种方式显示:Pr

28、essureofintsolidtop5.44e-025.16e-024.89e-024.62e-024.35e-024.08e-023.81e-023.53e-023.26e-022.99e-022.72e-022.45e-022.17e-021.90e-021.63e-021.36e-021.09e-028.15e-035.44e-032.72e-030.00e+00ContoursofStaticPressure(pascal)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)Pressureofintsolidisometrie5.44e-025.16e-024.8

29、9e-024.62e-024.35e-024.08e-023.81e-023.53e-023.26e-022.99e-022.72e-022.45e-022.17e-021.90e-021.63e-021.36e-021.09e-028.15e-035.44e-032.72e-030.00e+00ContoursofStaticPressure(pascal)Aug20,2010FLUENT6.3(3d,dp,pbns,lam)由以上两图可以看出圆管内部压强分布从管口处向延伸方向逐渐减小，可知流速相应增大，符合流速大，压强小的流动定律，也符合圆管流动压降的原理。另外从入口处的压强分布可以看出，

30、在圆管任何截面上，其压强分布也不是均匀的，也有分层现象。5.5轴向压强的变化执行PloAXYPlot,选择YAxisFunctic里的Pressur和PressureMagnitude选择Surface里圆管的对称轴line-x可得到轴向压强分布散点图。圆管层流中的压降，理论上存在下述公式128LqpV(1-13)d4即压降与流体的粘度、管道长度、流体的流量成正比，在本模拟实验中，由于流体的粘度、流体的流量不变，可认为压降与长度成正比，即P与L成正比。由上图可以看出，除了入段压强分布因流速急剧上升而下降过快外，其余部分均可看做是一条直线，即P随L的增加而降低，是正比关系。5.6总结报告5.6.

31、1系统总流量MassFlowRate(kg/s)inlet0.039138829int_solid-7.8335983outlet-0.039138853wall0Net-2.4280432e-085.6.2入出口流速积分IntegralVelocityMagnitude(m/s)(m2)inlet3.8362443e-05outlet3.9154264e-05Net7.7516706e-055.6.3入口出口压强积分IntegralStaticPressure(pascal)(m2)inletoutlet0.000388095040Net0.000388095046在FLUENT中求解计算紊

32、流流动6.1FLUENT设置除以下设置为紊流所必须设置的外，其余选项和层流相同，不再详述。Viscous设置U雷诺数Re10000，故圆管内流动为紊流，Viscous设置为RealizableK-epsi模型，其余默认。vBoundary设置Inle设置为速度入口为20.1m/s，Turbulenc锻置为IntensityandHydTicDia方法r即Outlet殳置为自由出DOutflow如设置成压力出，则之后计算会存在问题（已验证）。Solutioii殳置采用双精度求解器，定常流动，Realizablek模型，SIMPLEC算法。6.2开始迭代设置迭代次数为300，实际比这个少，迭代收敛

33、时会自动停止。IterationIterateApplyCloseHelp7紊流计算结果及分析计算293步后，已收敛，自动停止运算。残差监视窗口为Residualscontinuity-x-velocityy-velocity-z-velocitykepsilonScaledResidualsAug22,2010FLUENT6.3(3d,dp,pbns,rke)7.1显示流速等值线图7.1.1入口和出口截面的流速分布图分布在Surface里选择inle及outletVelocityofinlet1336-011236-011.16B4311.106-011J33e-019j99e-O23j9t

34、te-O2322e-027536-026S5e-026.1fe-O25.43e4Ei.79e-O2i.116-023.426-022.7fe-O22J35e-ffi1J?e-02cm4006S5e-O3Co他册BVeloc哼Magiltide如的Aig22,2010FLUENT阴閃中,帥庞）可见，入处流速分布不明显，基乎都等于入流速2OJm/s，只是外层靠近壁面处流速几乎为零。（2）Velocityofoutlet1376-D113-4311J3e-i11.16B4311.1QS4319J396-4323306-4323226-432T536-4326J85e4E6.1&-4325.4Se-l

35、2i.He-ffi3.426-4322.7-432snae-ffi13?e-4326J85e-E0om&400CommofVeloclVMagiltHdenA)Aig22,2010FLUENT6i3Qd.中pO河.rte)可见，出截面流速分布较为明显，和层流一样，显同心圆分布，内层流速偏大，外层靠近壁面处流速几乎为零，分层更为严重，边界层很薄。7.1.2Y轴和Z轴方向流速截面圆管内各个截面的流速分布均不相同，可以认为紊流还没达到稳定状态，在此不再分析各个截面的流速分布，仅对整个圆管的流速作出分析。截面沿圆管长度X方向截取，可看到对称的效果。（1）Velocityofy-0整根圆管：1.37e-

36、011.30e-011.23e-011.16e-011.10e-011.03e-019.59e-028.90e-028.22e-027.53e-026.85e-026.16e-025.48e-024.79e-024.11e-023.42e-022.74e-022.05e-021.37e-026.85e-030.00e+00ContoursofVelocityMagnitude(m/s)Aug22,2010FLUENT6.3(3d,dp,pbns,rke)（2）Velocityofz-0整根圆管：I3?e43l1336-4311236-111.16B431lJ33e-ilasae-reSjgce

37、-ce6S5&-4326.1&-4321.796-4321.116-4323.426-432Co他册BVeloclVMagiltide如的Aig22,2010FLUENT阴閃中,帥庞）2IEe-l2価沁6S5e-O3OID&400以上两个截面流速分布图的效果是一样的，可以看出圆管水流紊流入口段及之后的流速发展趋势，而且显示流速变化的规律更为明显（3）入口段1.37e-011.30e-011.23e-011.16e-011.10e-011.03e-019.59e-028.90e-028.22e-027.53e-026.85e-026.16e-025.48e-024.79e-024.11e-023

38、.42e-022.74e-022.05e-021.37e-026.85e-030.00e+00ContoursofVelocityMagnitude(m/s)Aug22,2010FLUENT6.3(3d,dp,pbns,rke)与层流入口段的流速分布相比，可以明显的看出紊流入口段的流速分布不太明显，且基本没有分层，符合紊流流动的基本规律。流速分布也不像层流流速那样显明显抛物线分布，而是更加平滑，越超后发展发展越平滑，到底是什么曲面，之后再加分析。紊流过流断面的流速对数分布比层流的抛物面分布均匀的u1多，符合InyC的规律，即uK4）出口段1.37e-011.30e-012.74e-022.05

39、e-021.37e-02ContoursofVelocityMagnitude(m/s)Aug22,2010FLUENT6.3(3d,dp,pbns,rke)6.85e-030.00e+00出口段的流层分布很明显，切趋于均匀，但仔细观察圆管轴心的速度，其实速度分布并未达到均匀，可见紊流并未达到充分发展的状况。7.1.3轴向流速的变化执行PloAXYPlot,选择YAxisFunctic里的Velocit和VelocityMagnitu，e选择Surface里圆管的对称轴line-x可得到轴向流速分布散点图。由上图可以看出，在圆管的轴上，进段流速分布变化较大，从进流速2lm/s急剧上升到最大流速

40、u0.1369m/s。之后又下降。但实际经验表明，紊流应该在进段后达到稳定状态，轴向流速max应该趋于恒定，可见此模拟实验设置长度不够，使流动并未达到充分紊流。紊流入口段长度有经验公式可以算的，即L2540d(1-14)由此可见，紊流的边界层厚度的增长比层流边界层要快，因此紊流的进口段要短些，而且长度主要受来流扰动的程度有关，与雷诺数无关，扰动越大，进段越短。可算得入段长度约为3m,由上图显示效果可以看出，轴向流速一直在变化，并未达到最大且稳定的速度，故紊流未发展充分。改进实验应加大圆管长度。7.1.4出截面的流速分布散点图因紊流并未充分，故选取出截面来进行分析(注意Plotdirect的选取)可见截面流速分布已很平滑，与层流出截面的流速分布截然不同。若紊流充分发展，则截面流速散点图最高处几乎为一条直线，说明圆管内大多数流体流速趋于稳定，几乎没有分层。取沿Y方向中