方柱绕流和圆柱绕流旋涡脱落数值模拟.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：086020 方柱绕流和圆柱绕流旋涡脱落的数值模拟 周云龙 邓冬 刁成东 洪文鹏 （东北电力大学能源与机械学院，吉林，吉林 132012 ） (Tel:0432-4806655, Email: zyl@mail.nedu.edu.cn) 摘要：本文采用了有限容积法，数值模拟了二维圆柱、方柱绕流的流动分离、旋涡生成、发展和脱落的瞬 态过程。通过模拟发现：圆柱绕流和方柱绕流的流体流动特征存在很多共同点和不同点。两者最大共同点 是都存在旋涡周期性脱落，但是方柱绕流的旋涡脱落频率比圆柱绕流的要小；最大的不同点是流体绕流分 离点的不同，圆柱绕流没有其固定分离点而方柱绕流固定分离点在其前后锐边上。本文数值结果与文献结 果吻合较好。 关键词：圆柱绕流 方柱绕流 旋涡脱落 数值模拟 0 引言 流体绕流物体的流动问题在工程实际中常可遇到，如风对各种建筑物的绕流，河水流过 桥墩，各种飞行器的设计，海洋石油工程中的开采平台、钻杆、水下输油管道等，在工业设 备中绕流现象更是经常发生，如各类管壳式换热器。因此掌握流体绕物体流动的特性对工程 实际和工业设备的设计非常重要，长期以来一直是学者们研究的热点问题。 对圆柱绕流尾迹结构和特性的研究已有半个世纪的历史。在实验研究方面，Roshoko[1] 最早在实验中发现圆柱绕流存在与流场 Re 数相关的三个不同的尾流状态，Taneda[2][5]应用流 场显示技术观测到圆柱尾迹随着下游距离的增加而不断演化成比初生涡频率要低的二次涡 等，人们完成了大量的实验研究。与此同时，随着计算机技术的发展，数值研究工作也取得 了很大的进展。Breuer[3]用大涡模拟（LES ）方法对大雷诺数（Re ＝140,000）的圆柱绕流问 题进行了数值模拟，Phuocloc Ta 和 Bouard R[4]数值模拟了绕流圆柱初期流动的二次涡结构, 得到与实验数据十分接近的结果。周云龙[6-7] 对气液两相流绕单个圆柱和多个串列，错列圆 柱旋涡脱落频率、圆柱表面升阻力系数等进行了大量的数值模拟，模拟结果与实验结果十分 相近。 与对圆柱绕流进行的大量研究相比，矩形柱体（包括方柱）绕流的研究相对少。Okajima[8] 调查了各个类型矩形柱体的漩涡脱落频率，并得出 Strouhal 数是雷诺数的函数。Saha,Biswas Muralidhar Senda et al.[9]用实验的方法研究了截面边长比分别为 0.5 和 1.0 矩形柱体，得到 了柱体尾涡脱落结构。Senda et al [10]建立三维模型，对雷诺数 50-500 的方柱绕流进行计算旨 在调查尾涡演变过程。 本文采用有限容积法，对不同雷诺数下的圆柱体柱体和方柱体绕流产生的旋涡脱落特性 分别进行了数值模拟，计算得到两种类型柱体绕流的升阻力系数、漩涡脱落频率、Strouhal 数、流函数、涡量等结果，并与文献结果进行了比较。 1 基本方程和数值方法 1.1 基本方程 二维不可压缩粘性流体非定常流动基本方程组为 ∂ρ ∂ ( ) + ρuj 0 （1） t ∂x ∂ j ∂τ p ∂( ) ∂ ( ) ∂ ij ρ ρu + u u − +