基于浸入边界算法的振动二维矩形柱绕流模拟研究  

Numerical simulation on the flow around an oscillating two-dimensinoal rectangular cylinder based on the immersed boundary method

作　　者：杨青 曹曙阳[2] 周军文[3] YANG Qing;CAO Shuyang;ZHOU Junwen(College of Civil Engineering,Anhui Jianzhu University,Hefei 230601,China;State Key Laboratory of Disaster Reduction in Civil Engineering,Tongji University,Shanghai 20092,China;Changzhou Key Lab of Construction Engineering Structure and Material Properties,Changzhou Institute of Technology,Changzhou 213032,China)

机构地区：[1]安徽建筑大学土木工程学院,合肥230601 [2]同济大学土木工程防灾国家重点实验室,上海200092 [3]常州工学院常州市建设工程结构与材料性能研究重点实验室,江苏常州213032

出　　处：《振动与冲击》2018年第14期254-261,共8页Journal of Vibration and Shock

基　　金：国家自然科学基金(51478358);建设部科技项目(2014-K2-014)

摘　　要：利用浸入边界算法,编写出可在静止网格中计算动边界绕流问题的数值模拟程序,对竖向正弦强迫振动矩形柱(Re=UD/v=10~3、振幅恒定、边长比B/D=1,3,5)展开绕流模拟,以探究边长比效应对矩形柱气动特性的影响。研究发现:柱体锁定区随边长比增大而发生变化,在B/D=3时,出现多级锁定现象;不同边长比矩形柱运动均会抑制同方向涡的产生;由于流动再附,B/D=3,5时,矩形柱侧面压力分布出现曲线特征;振动频率越高,柱体侧面前缘负压特征越明显,其侧面压力曲线分布特征亦同时增强;不同边长比矩形柱阻力系数均在接近自然涡脱频率处发生峰值效应。The immersed boundary method（ IBM） was adopted to compile a program to compute the moving boundary on a stable grid. The numerical simulation on the flow around a forced transversal sinusoidal oscillating rectangular cylinder（ Re = UD/v = 10~3,constant oscillating amplitude,side ratio B/D = 1,3,5） was carried out to investigate the influence of side ratio on the aerodynamic characteristics of the cylinder. The results indicate that the length of lock-in region changes with the varying of side ratio. When the side ratio B/D = 3,the synchronization zone is classified into two stages. The fluid flow near the surface in the direction consisting with the direction of cylinder displacement is suppressed by the motion of rectangular cylinder,which is a common flow feature near border,in all the three side ratio cases. When the side ratio B/D = 3 ～ 5,the curve characteristic of mean pressure distribution appears on the lateral surface of rectangular cylinder because of the flow reattachment. The negative surface pressure on the leading corner will become larger at high oscillating frequency. The mean drag coefficient in the three cases（ B/D = 1,3,5） all reaches its peak value near the natural vortex frequency under stable condition.

关 键 词：浸入边界算法 振动矩形柱 边长比 锁定区 表面压力 

分 类 号：O35[理学—流体力学]