利用FLUENT软件模拟圆柱绕流.docx
《利用FLUENT软件模拟圆柱绕流.docx》由会员分享,可在线阅读,更多相关《利用FLUENT软件模拟圆柱绕流.docx(9页珍藏版)》请在冰豆网上搜索。
![利用FLUENT软件模拟圆柱绕流.docx](https://file1.bdocx.com/fileroot1/2023-2/7/e1e858b5-8bf8-4bd0-a3ea-4648a1086cf7/e1e858b5-8bf8-4bd0-a3ea-4648a1086cf71.gif)
利用FLUENT软件模拟圆柱绕流
利用FLUENT软件模拟圆柱绕流
摘 要:
使用计算流体力学软件FLUENT,模拟均匀来流绕固定圆柱的流动,模拟雷诺数为20,40,100时的绕流流动,得到流场的流函数等值线图和速度矢量图。
计算结果表明:
当雷诺数增加时,流动表现出一系列不同的构造。
在雷诺数约为40前后流场有明显变化。
小于这个数时,存在一对位置固定的旋涡。
大于40时,流场开始变得不稳定,旋涡扩大、脱落、又生成,逐渐发展成两排周期性摆动和交错的旋涡。
并与实验及数值模拟结果比较,确认FLUENT能够很好地预测流动结构。
关键词:
圆柱绕流 FLUENT软件 雷诺数
1圆柱绕流理论分析研究的状况
一个世纪以来,圆柱绕流一直是众多理论分析、实验研究及数值模拟对象。
但迄今对该流动现象物理本质的理解仍是不完整的。
圆柱绕流中,起决定作用的是雷诺数,但还受到许多因素,如阻塞比,来流湍流度,下游边界条件等的影响。
随着雷诺数的增加,粘性不可压缩流体绕圆柱的流动会呈现各种不同的流动状态,在小雷诺数时,流动是定常的,随着雷诺数的增加,圆柱后会出现一对尾涡。
当雷诺数较大时,尾流首先失稳,出现周期性的振荡。
而后附着涡交替脱落,泻入尾流形成Karman涡街,随着雷诺数的增加,流动变得越来越复杂,最后发展为湍流。
White[1]认为“圆柱涡流具有经典性的重要意义”。
一般认为圆柱绕流有2种定常的流动图案:
雷诺数为较小时,圆柱后无尾涡;当雷诺数为较大时,圆柱后有一对对称的尾涡。
关于定常流失稳以及出现湍流的临界雷诺数主要是通过应用流场显示技术观察流动形态得到的,所以不是准确值。
对于分界点雷诺数就有不同的见解,Kovasznay、Roshko等认为定常流动失稳的临界雷诺数大约为40。
而从周期性尾流到湍流的详细的转变过程的实验研究似乎还是空白。
对均匀来流绕固定圆柱的二维平面流动,国内外许多学者进行过大量的研究。
决定圆柱绕流流态的是雷诺数
的值,
时,流动不发生分离,
在圆柱体后面出现一对位置固定的旋涡;
旋涡扩大,然后有一个旋涡开始脱落,接着另一个也脱落,在圆柱体后面又生成新的旋涡,这样逐渐发展成两排周期性摆动和交错的旋涡,即Karman涡街。
涡街是层流,
旋涡由层流向湍流转变。
称为亚临界区。
此时,柱体表面上的边界层为层流,而柱体后面的涡街已完全转变为湍流,并按一定的频率发放旋涡。
称为过渡区。
此时,柱体表面上的边界层也转变为湍流,分离点后移,阻力显著下降,旋涡的发放不规则,发放频率是宽频率随机的。
称为超临界区。
此时,重新建立起比较规则的准周期性发放的涡街。
目前,对圆柱绕流的三维数值研究则不多,然而试验表明,即使均匀来流垂直流过等截面的圆柱体,当Re足够大时,也会呈现三维的流动状态。
圆柱绕流属于非定常分离流动问题,在工业工程中的应用非常广泛,数值模拟是研究这类问题的有力工具。
但是,控制方程的高度非线性以及边界条件的多样性严重阻碍了这类问题的解决。
到目前为止,有限差分、有限元、大涡模拟(LES)、直接数值模拟(DNS)等方法都先后被应用。
但是,这些方法都需要增加网格生成的附加工作,尤其是在遇到复杂几何区域时更需要巨大的网格数量,大大影响计算效率。
本文利用N-S方程,对固定圆柱绕流进行了三维数值模拟,利用计算流体力学软件FLUENT,选取Re=20,40,100,模拟了圆柱周围的流场,得到流场的涡量等值线和速度矢量图,并与实验结果进行了比较。
2数学模型与数值方法
对不可压缩粘性流体,在直角坐标系下,其运动规律可用N-S方程来描述,连续性方程和动量方程分别为
圆柱周围的流场利用FLUENT求得,采用有限体积法和SMPLEC算法。
3数值计算与结果讨论
一个无穷长直径为
2.0m的圆面积柱体,放置在无穷远来流速度为1.0m/s不受干扰的均匀横流中,其他参数如图1、图2所示。
图1 计算域
雷诺数由圆柱体直径和自由来流速度确定。
通过改变粘度,得到不同的雷诺值(20,40,100),
当雷诺数为100时,试验的Strouhal数是0.165,
。
边界条件:
进口和侧面边界始终为自由来流条件
流动出口为零法向梯度出口边界。
流体特性假定为常数,结果见表1。
表1 流体特性
Re
20
40
100
密度ρ/(kg/m3)
1
1
1
动力粘度μ/(Pa2s)
0.1
0.05
0.02
网格示意图如图2所示。
图2 网格剖分
Re=20和Re=40时的流函数和速度矢量如图3、图4、图5、图6所示,流动是对称的。
图3 流函数等值线图(Re=20)
图4 速度向量图(Re=20)
图5 流函数等值线图(Re=40)
图6 速度向量图(Re=40)
图7 流函数等值线图(Re=100)
图8 速度向量图(Re=100)
4结论
(1)计算结果表明:
当雷诺数增加时,流动表现出一系列不同的构造。
在Re<40时,存在一对位置固定的旋涡。
Re>40时,流场开始变得不稳定,旋涡扩大、脱落、又生成,逐渐发展成两排周期性摆动和交错的旋涡。
与实验及数值模拟结果比较,确认FLUENT能够很好地预测流动结构。
(2)利用计算流体力学软件FLUENT可以成功地模拟圆柱绕流问题,反映出流动特性。
(3)可以对大于100的雷诺数做进一步的研究,探讨其机理。
FLUENT应用发展现状
FLUENT软件于1998年进人中国市场。
FLUENT软件在我国应用的也很广泛。
FLUENT软件在空气流场中的应用已比较成熟。
田铖等[2]模拟了地铁专用轴流风机的内部流场;李磊等[3]模拟了两种风向条件下街区风速分布和CO质量浓度的分布情况;杨伟等[4]计算了大气边界层中单栋高层建筑的定常风流场;赵琴[5]对空调室内的温度场进行了计算机模拟。
FLUENT软件在水流流场中的应用不同于单纯的空气流场。
水利工程问题边界复杂,水气混掺,存在多个自由水面,需要采用二相流模型,计算机模拟的成果还不多。
Salaheldin等[9]对垂直圆墩周围的三维流场进行了计算机模拟;刘加海等[6]在水槽中对二维规则波进行了计算机模拟;韩龙喜等[6]在环形水槽中对三维流场进行了计算机模拟。
Chen等[8]采用双方程紊流模型,对阶梯溢流坝流场进行了三维紊流计算机模拟,得出溢流坝沿程的自由水面位置、速度大小和分以及阶梯面上的压力等重要流场特性。
刁明军等[9]采用VOF法对挑流消能从库区到下游水垫塘进行了水气二相流计算机模拟,得到射流的空中轨迹线、水垫塘的自由水面和计算域的压力、流速、紊动能及紊动能耗散率分布。
杨忠超等[10]结合向家坝水电工程,采用VOF法和RNG紊流模型对多股多层水平淹没射流进行三维计算机模拟,得出中孔和表孔联合泄洪时消力池的流速和紊动能分布。
陈云良等[11]根据某电站进水口1:
50模型试验结果,对进水口流场进行了计算,研究来流、复杂地形边界等对进水口流态的影响。
张挺等[12-13]对索风营电站X形宽尾墩和阶梯溢流坝联合消能的三维流场进行数值计算,给出压力特性、流速分布、阶梯及墩后水气两相流的部分特性。
周勤等[14]对小南海水库双弯曲型溢洪道原型流场进行三维计算机模拟,得出溢洪道沿程水面线和速度等重要水力参数。
叶茂等[15]
结合典型工程算例,就复杂水气交界面的处理、模型建体和网格划分等问题作详细介绍。
罗挺等[16]以windows2000操作系统和通用CFD软件工具FLUENT为基础,研究构建了不同的并行计算硬件平台的力法,计算了典型的山谷污染物扩散过程的湍流流场。
可见,FLUENT软件在水力学中的应用已逐步展开,计算机模拟同试验研究相比,具有成本低、方案变化快、无测量仪器干扰、无比尺效应和数据信息完整等优势,可作为模型试验的有力补充,如进水口模型试验在高水位时,无法用常规仪器测量进口周围的流速,而数值计算则可得到全域流场特性。
另外,由于其减少了研究者在计算方法、编程、前后处理等方而投入的重复、低效的劳动,将更多的精力和时间投入到考虑问题的物理本质,优化算法选用,参数的设定,因而提高了工作效率,其必然会获得越来越多的应用,应用效果也必将越来越好。
参考文献:
[1]FrankM.White.ViscousFluidFlow[M].McGraw2HillInc,1974.
[2]田铖,张欢,由世俊,等.利用FLUENT软件模拟地铁专用轴流风机的内部流场
(一):
对称翼叶片轴流风机[J].流体机械,2003,31(11):
13-15.
[3]李磊,胡非,程雪玲,等.FLUENT在城市街区大气环境中的一个应用[J].中国科学院研究生院学报,2004,21(4):
476-480.
[4]杨伟,顾明.高层建筑三维定常风场数值模拟[J].同济大学学报,2003,31(6):
647-651.
[5]赵琴.FLUENT软件的技术特点及其在暖通空调领域的应[J].计算机应用,2003,23(12):
424-425.
[5]SALAHELDINTM,IMRANJ.CHAUDHRYMH.Numericalmodelingofthree-dimensionalflowfieldaroundcircularpiers[J].JournalofHydraulicEngineering,2003,130
(2):
91-100.
[6]刘加海,杨永全,张洪雨,等.二维数值水槽波浪生成过程及波浪形态分析[J].四川大学学报:
工程科学版,2004,36(6):
28-31.
[7]韩龙喜,WAJWing-hong.加长型环形水槽水流特性的数值模拟[J].河海大学学报:
自然科学版,2003,31(6):
639-643.
[8]CHENQun,DAIGuang-qing,LIUHao-wu.Volumeoffluidmodelforturbulencenumericalsimulationofsteppedspillwayoverflow[J].JournalofHydraulicEngineering,2002,128(7):
683-688.
[9]刁明军,杨水全,王玉蓉,等.挑流消能水气二相流数值模拟[J].水利学报,2003,9:
77-82.
[10]杨忠超,邓军,杨永全,等.多股多层水平淹没射流数值模拟研究[J].水利学2004,5:
31-38.
[11]陈云良,伍超,叶茂.水电站进水口水流[12]陈云良,伍超,叶茂.水电站进水口水流流态的研究[J].水动力学研展,2005,20(3):
340-345.
[13]张挺,伍超,卢红,等.X型宽尾墩+阶梯滋流坝流场三维数值模拟[J].水利学报,2004,8:
15-20.
[14]ZHANGTing,WUChao,LIAOHue-sheng,etal.3Dnumericalsimulationonwaterandairtwo-phaseflowsofthestepsendflaringgatepier[J].JournalofHydrodynamics:
Ser.B,2005,17(3):
338-343.
[15]周勤,伍超,赵元弘,等."S"型溢洪道水流特性试验与数值模拟研究[J].水力发电学报,2005,24(3):
88-92.