圆柱绕流的数值模拟讲解.doc

圆柱绕流的数值模拟 张玉静 20070360204 过控（2）班 化工与能源学院 摘 要:使用计算流体力学软件FLUENT,模拟均匀来流绕固定圆柱的流动,模拟雷诺数为5，20,40,100时的绕流流动, 得到流场的流函数等值线图和速度矢量图。计算结果表明:当雷诺数增加时,流动表现出一系列不同的构造。当Re=5时，流动不发生分离,其后未形成旋涡，当Re=20,40,100时，流体发生分离，其后形成旋涡，且旋涡随着Re的增大而增大。利用计算流体力学软件FLUENT可以成功地模拟圆柱绕流问题,反映出流动特性。 关键词:圆柱绕流；FLUENT；雷诺数 Abstract:Uniform flow around a mounting cylinder is simulated with the application of FLUENT software while Reynolds number is 5,20,40,100. Stream function and velocity vector distributions are indicated. The results show that a series of construction appears as Reynolds number increases. When Re is 5, Flow separation does not occur, and it does not form vortex . When Re is 20,40,100, Flow separation occurs, and it forms vortex. Vortex increases with the increase of Re. Using computational fluid dynamics software FLUENT can successfully simulate flow around cylindrical, reflect the flow characteristic. Key words：Flow around a circular cylinder；FLUENT；Reynolds number （1） 式中，u、v、w表示速度在x、y、z方向上的分量。 （2）动量方程 （4） 式中，xi（i=1,2,3）为坐标系坐标，xj(j=1,2,3)为坐标系坐标，ui(i=1,2,3)为沿i方向的速度分量，fi为沿i方向的质量力，P是静压，ρ为空气密度。Gambit中网格划分结果如图2所示。 图1 计算区域 图2 网格划分 4边界条件与求解设置 进口为Velocity-inlet，出口为Outflow，其他边界默认为Wall，计算区域设为fluid，在gambit中设置界面如图3所示。 图3边界条件置界面 图4指定流体区域 求解设置如下： （1）控制方程采用连续方程和动量方程，而不考虑能量方程。 （）求解器选择COUPLED IMPLICIT SOLVER，选择Steady状态。 （）为了提高计算精度，差分格式采用二阶迎风格式。 （）速度场和压力场的求解基于Simple算法，松弛因子默认。 （）模型选择。 （）残差设为e-6。fluent中的操作具体如下： （1）读入从gambit中导出的msh文件，检查网格质量，并确定模型单位，如图5所示。 （2）设置求解器，如图6所示。 （3）选择层流模型，如图7所示。 图5 确定单位 图6 选择求解器 图7 选择层流模型 （4）改变流体介质，如图8所示。 （5）设置入口速度，如图9所示。 （6）设置残差，如图10所示。 （7）初如化流场，进行迭代求解。 图8 改变流体介质 图9 设置入口速度 图10 设置残差 5数值计算与结果讨论 雷诺数由圆柱体直径和自由来流速度确定。通过改变粘度,得到不同的雷诺值(5，20,40,100),Re =ρUD/μ。 边界条件:进口和侧面边界始终为自由来流条件U=1.0m/s,流动出口为零法向梯度出口边界。 流体特性假定为常数,结果见表1。 表1 流体特性 Re 5 20 40 100 密度ρ/(kg/m3) 1 1 1 1 动力粘度μ/(Pa.s) 0．4 0.1 0.05 0.02 Re=20时的残差曲线如图11所示。 图11 Re=20时计算的残差曲线 (a) Re=5时的速度分布 （b）Re=20时的速度分布 (c) Re=40时的速度分布 (d) Re=100