弯管方腔圆柱绕流报告.docx

1、弯管方腔圆柱绕流报告混合弯管、方腔、圆柱绕流上机实习报告一、混合弯管1.建模使用Gambit软件建立2D混合弯管模型1）建立大管模型，管径为16，弯曲处角度为90；2）在大管弯曲处建立小管与大管相接，管径为4；3）划分网格：直管段处网格大小为1，弯曲处则加密使计算更精确，网格数设置为52；图14）选择求解器Fluent5/6；5）选择大管上边界，设置为outflow；选择大管左边界，命名为inflow-cold，选择小管下边界，命名为inflow-hot，由于该求解器中无inflow的设置，可先设置为wall，再在fluent中修改；6）导出msh文件。2.计算使用Fluent软件计算 1）导

2、入msh文件，检查确定无负体积； 2）设置材料为water-liquid，密度等参数使用默认值；3）设置边界条件，将inflow-cold边界设为velocity-inlet，速度为5m/s，温度为280K，inflow-hot边界也设为velocity-inlet，速度为8m/s，温度为350K；4）计算精度默认为10e-3；5）初始化条件；6）进行迭代计算。3.计算结果从fluent中导出图像：速度矢量分布图：图2速度大小分布图：图3流线图：图4温度分布图：图5可以看到在大管和小管流体混合区域，温度的分布是沿着流线的。二、方腔环流1.建模 使用Gambit软件建立方腔模型 1）建一个边长为

3、1的正方形；2）将每边均分为200份，即网格数200200，网格大小为0.005；图63）选择求解器Fluent5/6；4）选择正方形的上边界，设置为wall，命名movewall；5）导出msh文件。2.计算使用Fluent软件计算 1）导入msh文件，检查确定无负体积； 2）设置材料为water-liquid，密度为1000kg/m3，粘性系数为0.001Ns/m2；3）设置movewall边界条件，设置为moving wall，绝对速度为0.1m/s，则雷诺数为Re=10e5m2/s；4）计算精度默认为10e-3；5）初始化条件；6）进行迭代计算。3.计算结果调出stream funct

4、ion即流线图，如下图所示：图7可以看到方腔中间形成环流，方腔的右下角、左下角及左上角分别有一个涡。4.对比分析分别改变上边界速度u、方腔边长D、粘性系数，得到不同的流线图，对比如下：1）保持u及D不变，改变粘性系数，则雷诺数Re将随之改变：增大粘性系数，雷诺数减小D=1，=0.01，u=0.1，Re=10e4图8缩小粘性系数，雷诺数增大D=1，=0.0001，u=0.1，Re=10e6图9 可以看到在D=1，u=0.1的条件下，及雷诺数的改变对流线改变不大，流线图几乎没有变化，最大流函数值保持不变。2）保持u及不变，改变方腔尺寸D，雷诺数Re随之改变：缩小方腔尺寸D，雷诺数减小D=0.1，=

5、0.001，u=0.1，Re=10e4图10注：由于流线图level有限，故只显示部分流线，使边角的涡显示出来，下同。增大方腔尺寸D，雷诺数增大D=10，=0.001，u=0.1，Re=10e6图11可以看到在=0.001，u=0.1的条件下，Re较小时，只能看到右下角的涡，且涡的范围较大，由于level有限，左下角涡难以显现，中心区域的环流并不很圆，而是沿着上板速度方向有变形；Re较大时，三个涡更为明显，右下角涡范围变小，左上角涡范围变大，中心区域的环流变圆，区域变大。3）保持D及不变，改变上板速度u，雷诺数Re随之改变：减小上板速度u，雷诺数减小D=0.1，=0.001，u=0.01，Re

6、=10e4图12增大上板速度u，雷诺数增大D=0.1，=0.001，u=1，Re=10e6图13 进一步增大上板速度D=0.1，=0.001，u=10，Re=10e6图14在=0.001，D=0.1的条件下，环流及涡随雷诺数的变化趋势与=0.001，u=0.1时一致，当u=10m/s，即雷诺数达到10e7 m2/s时，右下角涡消失，左下角涡及左上角涡也变小，其中左上角涡形状变得较为规整，中间环流几乎占满整个方腔，中心已成为圆形流。4）保持D不变，通过改变上板速度u和粘性系数，使雷诺数Re同样保持不变： 减小上板速度u和粘性系数D=1，=0.0001，u=0.01，Re=10e5图15增大上板速

7、度u和粘性系数D=1，=0.01，u=1，Re=10e5图16雷诺数和方腔尺寸不变而速度有所改变时的流函数变化趋势与前两项讨论一致。5）保持粘性系数不变，通过改变上板速度u和方腔尺寸D，使雷诺数Re同样保持不变： 减小上板速度u，增大方腔尺寸DD=10，=0.001，u=0.01，Re=10e5图17增大上板速度u，减小方腔尺寸DD=0.1，=0.001，u=1，Re=10e5图18雷诺数和粘性系数不变而速度有所改变时，流线改变不大，流线图几乎没有变化，最大流函数值变化较小。6）保持上板速度u不变，通过改变粘性系数和方腔尺寸D，使雷诺数Re同样保持不变： 减小粘性系数和方腔尺寸DD=0.1，=

8、0.0001，u=0.1，Re=10e5图19增大粘性系数和方腔尺寸DD=10，=0.01，u=0.1，Re=10e5图20雷诺数和上板速度不变而方腔尺寸和粘性系数有所改变时的流函数变化趋势与前相同。7）改变流体材料，将流体材料设置为air： 保持方腔尺寸D及上板速度u不变D=1，=1.225，=1.7894e-5，u=0.1，Re=6846图21保持方腔尺寸D及雷诺数Re不变D=1，=1.225，=1.7894e-5，u=0.146，Re=10e5图22 由图可知改变材料后流线图变化亦不大。各工况具体数据见下表：表一编号12345678网格数200200200200200200200200h

9、0.0050.0050.0050.00050.050.0050.0050.005D1110.110111100010001000100010001000100010000.0010.010.00010.0010.0010.0010.0010.001u0.10.10.10.10.10.01110Re10e510e410e610e410e610e410e610e7不变量u、Du、Du、u、D、D、D、迭代次数19271927192715902323137722322022max9.98e-39.98e-39.98e-31.24e-35.35e-21.2e-35.32e-20.21表一续编号9101

10、11213141516网格数200200200200200200200200h0.0050.0050.00050.050.00050.050.0050.005D110.1100.110111000100010001000100010001.2251.2250.00010.010.0010.0010.00010.011.79e-51.79e-5u0.01110.010.10.10.10.146Re10e510e510e510e510e510e5684610e5不变量D、ReD、Re、Re、Reu、Reu、Re迭代次数13772232193718841590232319272016max1.2e-

11、35.32e-29.96e-39.98e-31.24e-35.35e-29.98e-31.35e-2由以上各项对比可以看出初步的规律：当上板速度u和方腔尺寸D均不变或雷诺数Re和粘性系数均不变时，流线图几乎不变；当四项主要参数中有两项变小，另两项保持不变时，流函数值减小，流线图中涡变得不明显；当四项主要参数中有两项变大，另两项保持不变时，流函数值增大，流线图中三个涡变得明显，中间环流趋于圆形。三、圆柱绕流1.建模使用Gambit软件建立圆柱绕流模型1）流场区域为x（-11.5,20），y（-12.5,12.5）；2）圆柱直径为1，圆心位于原点；3）划分网格：将整个流场区域划分为9个部分；其中圆

12、柱周围再细分为4部分，每小部分网格数为2525；其余矩形部分网格大小为0.2；图23图244）选择求解器Fluent5/6；5）选择流场区域左边界，设置为inlet，先设置为wall，再在fluent中修改；选择右边界，命名为outlet，设置为outflow，中间圆柱边界设置为wall；6）导出msh文件。2.计算使用Fluent软件计算 1）导入msh文件，检查确定无负体积； 2）设置材料密度为150kg/m3，粘性系数为1Ns/m2，进口流体速度为1m/s，则雷诺数为Re=150m2/s；3）采用定常算法，计算精度默认为10e-3；也可采用非定常算法，time step size可设为0.2；4）初始化条件；5）进行迭代计算。3.计算结果定常算法涡线图：图25流线图：图26非定常算法：涡线图图27流线图图28可以看到在圆柱后方出现了交替的两排涡，非定常算法所得的结果更精确。四、小结1.可使用Gambit建模，得到mesh文件，再导入Fluent中计算，用Gambit建模更方便简易；2.在Fluent中，可通过改变材料密度、粘性系数、流体速度等改变雷诺数，以改变流体性质；可调出速度大小分布图、矢量图、流线图、涡线图等，通过改变level数调整数据分级，以改变图像的显示，如方腔环流中只显示一部分流线以使右下角的涡得以显示。