基于FLUENT的方块绕流问题探究.docxVIP

基于FLUENT的方块绕流问题探究 研究模型 图1 模型简图 如图1所示，一长方形流场中放置一边长为a的方块，左侧为速度入口，流场宽度b与a的关系为b/a=10，分别研究在Re=1，10，100时，流场的不同。 本模型中取a=0.05m，b=0.5m，流场长度为1m，方块位于靠近入口约为全长四分之一的地方。 网格划分 对本模型划分网格用ICEM软件进行结构网格： 图2 网格划分 如图2，对方块周围流场进行加密，以保证计算精度，并能顺利捕捉到涡街，网格数量约为14000。 参数设置 温度为20oC水的运动粘度υ＝1.007e-6m2/s，由雷诺数定义Re＝ud／υ，Re=1，10，100对应的流速为0.00002m/s，0.0002 m/s，0.002 m/s.具体设置步骤如图。 将画好的网格文件导入FLUENT，首先对网格文件check，确认没有负体积出现。 图3 网格检查图 设置求解器： 图4 求解器设置 此模型应选择非定常，即Transient。 模型选择： 图5 模型选择 此模型雷诺数明显应该选择层流。 流体性质： 图6 调出水的性质 图7 工质选择液态水 工作介质选择液态水，其参数直接从材料库中调出。 边界条件： 图8 入口设置 入口设置为速度入口，速度根据雷诺数不同分别设为0.00002m/s，0.0002 m/s，0.002 m/s。 图9 出口设置 出口边界设置为outflow。 求解方法： 图10 求解方法设置 求解算法选用最为普遍的SIMPLE算法，离散方法对压力和动量均采用二阶精度，瞬态方程采用二阶隐式，以提高计算精度。 残差设置： 图11 残差设置 需要收敛的残差设置为1e-05。 迭代设置： 图12 迭代设置 时间步长设置为0.5s，每个时间步长内的迭代数为40步，保证能在时间步内收敛。 计算结果 图13给出不同雷诺数下的速度场 Re=1 （b）Re=10 （c）Re=100 图13 速度分布云图 由不同雷诺数的速度分布云图可以看出： 在Re=1时，方块上下两侧速度场有所增大，方块四周速度由于阻力降低，总的来说，方块对速度场影响不大； 在Re=10时，方块上下两侧流速增大的影响区域增大，直至出口才趋于平缓，方块后流速降低的区域也有所增加； 在Re=100时，在方块后，可以看到明显的涡从方块后的流场尾迹处脱落并呈上下交替。 图14 给出不同雷诺数的涡??云图 （a）Re=1 （b）Re=10 （c）Re=100 图14 涡量分布云图 从涡量分布云图中可以看出： 在Re=1时，涡量变化主要集中在方块周围，且影响区域很小，随着雷诺数增大，Re=10时，涡量变化区域向下游延伸，但依然对称，当Re=100时，可以明显看到涡的脱落，上下交替。 由速度和涡量在不同雷诺数下的分布云图，可以看出在雷诺数大到一定数值时，会在方块后的区域出现涡的生成和脱落现象，出现“卡门涡街”。 为了看出不同雷诺数下涡街对流场的影响，图15为流线图： Re=1 Re=10 Re=100 图15 流线图 图15中可以看出，小雷诺数时，流线光滑平顺，方块后只有很小的回流区域，而随着雷诺数增大，雷诺数为100时，出现涡街后，流线受涡街影响，波浪形上下摆动，回流区域增大。 图16为Re=100时速度矢量图： （a） （b） 图16 Re=100时，速度矢量的分布 由图16可以看出在方块的后方会出现回流，方块前流线密集，速度较大，而方块后流线稀疏，流速较小。流动的方向随着涡街上下摆动，其他地方流线光滑平顺。 总结 本作业对Re=1，10，100这3种不同情况下，方块绕流问题进行探究，对比不同雷诺数下的速度和压力云图，得出以下结论： 流体流动经过方块形障碍物时，方块两侧流???增加，并随着雷诺数的增大，高速区越大，方块后流速很低，并出现回流，雷诺数越大，回流区域越大； 随着雷诺数的增大，当Re大于某个值时，方块后会出现上下交替的脱落的涡，出现“卡门涡街”； 小雷诺数时，流线光滑平顺，涡量变化集中在方块附近，而当雷诺数为100时，流线受涡量脱落扰动，呈波浪形上下摆动。