柱体绕流阻力理论研究和数值模拟.docx
- 文档编号:27920277
- 上传时间:2023-07-06
- 格式:DOCX
- 页数:9
- 大小:663.86KB
柱体绕流阻力理论研究和数值模拟.docx
《柱体绕流阻力理论研究和数值模拟.docx》由会员分享,可在线阅读,更多相关《柱体绕流阻力理论研究和数值模拟.docx(9页珍藏版)》请在冰豆网上搜索。
柱体绕流阻力理论研究和数值模拟
柱体绕流阻力理论研究和数值模拟
李金钊
2011-11-27
柱体绕流阻力理论研究和数值模拟
摘要:
柱体绕流是流体力学中的经典问题之一,而对于绕流阻力的研究又是该问题的尖键
之处。
本文对柱体绕流阻力产生的原因进行了理论分析,并对国内外尖于柱体绕流阻力的研究成果进行了归纳总结,指出了研究的方向和对前景的展望。
同时本文借助于
Fluent软件,针对二维圆柱和方柱绕流进行了数值模拟,得出了绕流阻力系数并与相矢
试验结果进行了比较分析。
尖键词:
绕流阻力、研究成果、数值模拟
1刖吞1.1柱体绕流阻力研究意义及应用背景
流体绕结构物流动时,过流断面收缩,流速沿程增加,压强沿程减小。
由于黏性力的存在,就会在物体周围形成边界层的分离,形成绕流。
而由于结构物的存在,会在物体迎水面产生雍水现象,同时也增加了结构物的受力,使得绕流问题变的十分复杂。
目前相尖理论研究成果较贫乏,因此对绕流现象进行研究具有重要的理论基础意义。
研究结构物绕流问题在工程实际中也具有重要意义O如在水流对桥下部结构的作用中,风对桥塔、索缆的作用中,都有重要的工程应用背景。
因此对结构物绕流进行深入研究,掌握其流动机理和水动力学规律,不仅具有理论意义,还有明显的社会经济效益。
1-2绕流阻力理论分析
水流流经柱体时,作用于物体上的力可分为两类:
摩擦阻力和压差阻力。
其中流体作用于物体表面的摩擦力在水流方向上的投影就是摩擦阻力。
压差阻力的产生是由于物体表面边界层产生分离。
边界层的分离常常伴随着涡旋的产生和能量的损失,从而物体前后面压强发生变化,产生了压强差了,增加了流动的阻力。
压差阻力主要取决于物体的形状,因此也称为形状阻力。
对于细长物体,例如顺水流放置的平板或翼型,则摩擦阻力占主导地位;而钝,性物体的绕流,例如圆球、桥墩等,则主要是压差阻力。
液体对物体的绕流阻力可用下列公式计算[1]
(1)
fd二Cd
式中,Cd为绕流阻力系数;A为物体在流速垂直方向的迎流投影面积
2研究进展
鉴于柱体绕流阻力问题在理论和工程实际中的重要意义,国内外许多学者对此问题进行了许多研究工作,大多数是通过对柱体受到的作用力和流速的测量,来确定其绕流阻力系数。
研究的方法包括试验研究、原型观测和数值模拟等。
2.1试验研究与原型观测
1985年,N.E•伊杰里奇克通过试验得出了在均匀来流条件下,光滑圆柱体
阻力系数Cd与雷诺数Re有很大尖系[2],即Cd二f(Re,并得出了相应的矢系曲
线,如图1所示
图1阻力系数与雷诺数矢系曲线图
由上图曲线可看出:
当Re很小时(例如Rev0・5),惯性力与粘滞力相比可以忽
略,阻力与Uo成正比。
阻力系数Cd则与Re成反比,如图中的直线部分,这时流动称为蠕动。
当Re增加,在圆柱表面产生了层流边界层,而一旦发生了
边界层的分离(Re壯5),则压差阻力大大增加。
当Re肚200,发生卡门涡街。
当Re增加到104时,压差阻力占主要,相对比摩擦阻力很小,故绕流阻力与Re几乎无尖。
当Re
增加至310时,由于圆柱表面的层流边界层开始转变为紊流边界层,动能增大,Cd突然
下降。
该结果对单柱体绕流阻力的研究起到了重要意义。
为便于对比,本文将国内外其他研究成果列于下表1。
表1国内外各阻力系数汇总表[3]
080
073
i.O
073
选用费料名竹園柱但力攝4tCp址
3柱体绕流数值模拟
本文采用有限容积法(FVM),借助于Fluent软件,对二维圆柱体和方柱体绕流进
行数值模拟,计算得到了两种类型绕流流线图和阻力系数等结果,并与相
尖文献结果进行了比较。
3.1数值计算
式中i,j=1,2,3;「为流体密度「••为运动粘度系数
2)计算域和边界条件
3.5DI
图2圆柱绕流计算区域
■
1
J
1
U=lt.v=()
3.5D
ED
图3方柱绕流计算区域
显然要指定来流面、出流面、结构物和区域表面的边界条件。
来流面:
为保证流速得到充分发展,使用UDF定义其入口流速边界条件;
(4)
u=1.5%[1一岸一1)2][1]
出流面:
指定为压力出口;
柱体表面和区域上下界面:
指定为壁面(wall),u=°v°,也称无滑移边界条件。
3)离散格式及求解
控制方程使用有限体积法离散,对流项采用二阶迎风格式,扩散项采用中心差分格式,速度和压力耦合采用SIMPLE(算法,最后进行迭代求解。
3.2计算结果与分析
本文采用2D层流模型,流动为定常稳态流,分别计算了雷诺数Re=40、200、800时的圆柱、方柱绕流问题。
图5、6分别显示了不同雷诺数下的圆柱、方柱绕流流线图。
从图中可看出,雷诺数Re=40时,柱体后面产生了一组排列对称的涡旋,这是由于边界层的分离致使;雷诺数Re二200时,柱体后面出现了涡脱落现象,涡街为层流涡街;雷诺数Re=800时,边界层仍为层流分离,而尾流已转为紊流涡街。
表2列出了圆柱绕流阻力系数Cd的计算结果,通过与上述文献试验资料[2]比较分析可看出,数值模拟的计算结果与试验结果存在一定的偏差。
这是由于本文数值模拟是二维情况下的柱体绕流,而试验结果是在三维的条件下得出的,所以简化的二维模拟与实际情况会出现偏差。
表2阻力系数Cd
阻力系数Cd
Re=40
Re=200
Re=800
数值模拟
3.2
2.5
1.1
试验结果
2.4
1.4
0.9
偏差
33.3%
78.6%
22.2%
Re=200
图5圆柱绕流流线图
Re=40
图6方柱绕流流线图
Re=800
4柱体绕流研究方向展望
传统的试验验研究由于受其本身局限性(如模型尺寸、试验周期、资金投入等)的影响,很多试验研究很难实现;由于流体力学方程为非线性方程,无法得到其理论解析解,使得理论流体力学的发展受到限制。
随着CFD技术的发
展,不仅克服了传统试验研究的限制,又由于其强大的计算能力能够解决解析法无法解决的微分方程,因此CFD技术的发展大大促进了流体力学研究的进步o同时,我们也应该看到在数值计算中如何有效的划分网格,确定离散格式和算法是保证数值模拟结果合理正确性的尖键。
由于实际工程实践的复杂性,绕流问题更多的是三维流态的情况,采用2D进行简化模拟有时欠妥,因此如何使用3D模型对柱体绕流进行正确合理模拟还需要深入研究。
随着测试技术的发展,粒子图像测试(简称PIV)[7]与CFD技术的结合已被实践证明是一种高效、可靠的流体流动研究方法。
工程界已开始利用PIV与
CFD结合的方法设计反应器,为目前仅以经验为主进行反应器放大设计的方法开拓了一条新路°
5结语
本文对已有国内外研究资料进行了总结并对二维柱体绕流进行了数值模拟,可看出简化后的二维模拟与试验结果有所偏差,而尖于三维绕流阻力问题的理论研究不多,数值模拟也大多是二维绕流的情况。
因此,对该问题进行深入研究和完善具有重要意义。
参考文献
[1]李玉柱,贺五洲•工程流体力学[M].北京:
清华大学出版社,2006.
[2]N.E伊杰里奇克•水利摩阻手册[M]•北京:
航空发动机编辑部,1985.
[3]杨纪伟,付晓丽・圆柱绕流研究进展・中国水运[J],2005.5,第08卷第5期
1992.12‘A辑第七增刊.
[5]Breuer,M.Achallengingtestcaseforlargeeddysimulation:
highReynolds
numbercircularcylinderflow.IntJ.HeatandFluidFlow.2000,21:
648-654
[6]姚熊亮等•基于LES方法圆柱绕流三维数值模拟•水动力学与进展[J],A辑第22
卷第5期
[7]张孝棣等•圆柱体绕流尾迹的PIV测量•试验流体力学[J],第19卷第2期
[8]Okajima,A.Strouhalnumberofrectangularcylinders.JournalofFluidMechanics.1982,123:
379-398
[9]白桦,李家武,夏勇•低雷诺数圆柱绕流数值模拟及控制措施•建筑科学与工程
学报[J],2010.12,第27卷第四期
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 柱体 阻力 理论研究 数值 模拟