方柱绕流的数值模拟.docVIP

方柱绕流的数值模拟 动力工程 易玉微 2013G0502007 摘要：为了研究方柱绕流中的流动状况，本文采用CFD软件，对当Re=100时，与来流方向呈0°和45°夹角两种情况进行了数值模拟。从模拟结果可以得出夹角为45°的方柱对尾部的涡流影响更大。在不同时间步时的速度矢量图，可以看出旋涡的生产、成长和脱落的全过程。结果验证了特征长度越大对柱体尾涡影响越大和同一雷诺数下速度越小涡的周期越长。 关键字：方柱绕流；特征长度；数值模拟 1、引言 流体撞击障碍物之后，其流场情况研究具有很大的工程价值，如，可通过模拟风载荷或水荷载对建筑物的作用力、剪切力等分布，可以对建筑物的形状优化、结构选择等产生很大的理论价值。国内外许多学者对这个问题进行了研究。如Antoine Placzek等[1]对震荡中的主体在低雷诺数的流场进行了数值模拟，Yoshida等[2]计算在流场中央的方柱的流场速度、压力等分布与实验结果进行对比结果很好的符合，李玲等[3]对处于浅水层直径为1.27m圆柱进行数字图像测速实验，对浅水层柱体绕流进行很好地解释。本文采用FLUENT软件，对与来流方向呈0°和45°夹角的方柱绕流进行数值模拟，得出其流场情况并作出分析，与现有的文献进行比较。 2、数值模拟 2.1、计算区域及边界条件 如图1和图2，二维计算区域为20D×50D，其中D为方柱变长。在GAMBIT中进行建模与网格划分，其中D=2mm。定义入口为velocity inlet，出口为outflow，其余边界为无滑移固壁边界条件。 图1 与来流方向呈0°的计算区域 图2 与来流方向呈45°的计算区域 2.2控制方程 本文方柱绕流场内的相关流体动力学控制方程，如下： 连续性方程： (1) x方向动量方程： (2) y方向动量方程： (3) 其中，u为x方向速度，v为y方向速度，P为压力，ρ为流体密度，μ为动力粘度。方柱绕流的Re=100，其公式如下： (4) 其中，l为特征长度，u为特征速度。对于与来流呈0°的方柱其l=D=2mm，与来流呈45°的方柱其。通过计算模型一的来流速度，=；模型二的来流速度，。 对于方柱绕流的相关动力学参数，包括阻力系数，升力系数和斯特劳哈尔数St。 (5) (6) (7) 式中，为阻力，为升力，D为方柱边长，为来流速度，为旋涡脱落频率。 2.3数值方法 本文的求解器采用FLUENT软件，流动采用层流模型求解，压力速度耦合采用SIMPLCE算法，空间梯度差分采用Least Squares Cell Based格式，压力采用Second Order格式，动量方程采用二阶迎风格式格式求解。时间步长等于最小网格长度除以流速或者旋转流动速度，保证每次迭代都在一个网格的范围内，因为计算一段时间之后可能加快计算速度，会将计算步长设置变大，所以要将时间步长设置小些，本文案例计算时间步长均设置为0.001s。计算区域内的流体为水流，密度为998.2 kg/m3，动力粘度为0.001003 Pa·s。 2.4网格无关性验证 通过保证连续性方程及x、y方向的动量方程收敛来判断数值模拟是否收敛，残差曲线如图3。 图3 残差曲线图 为了保证计算精度，证明计算结果与网格数量没有关系，对网格进行无关性验证。本文对与来流方向呈0°的模型一进行网格划分，其网格数分别为23600、95600、599600进行计算。在计算区域建立一条过方柱中心的竖直线，计算等到改线上的速度大小分布，如图4。从图4可以看出，在网格数为23600、95600、599600时，竖直线上的速度分布几乎重合。由此可知，计算结果与网格粗密无关。 图4 竖直线上速度分布 3、结果和分析 图5为模型一的升力系数和阻力系数随时间变化曲线。图6为模型二的升力系数和阻力系数随时间变化曲线。每个旋涡的生产、成长和脱落定为一个周期，其周期时间即升力的周期时间。观察图5和图6可知，升力系数在计算到一定时间后围绕着0开始有规律的上下脉动，模型一的脉动周期大概为0.