基于浸入式边界方法的串联双矩形柱绕流数值模拟.pdfVIP

物 理 学 报 ActaPhys．Sin．Vo1．63，No．21(2014) 214702 基于浸入式边界方法的串联双矩形柱绕流 数值模拟术 杨青 曹曙阳十 刘十一 (同济大学土木工程防灾国家重点实验室，上海 200092) (2014年4月25日收到；2014年6月5日收到修改稿) 基于浸入式边界算法 (VirtualBoundaryMethod)中力源反馈边界的思想，改进其原有内部流体处理方 法 以减少计算耗费，并结合非等间距网格 以便工程应用计算，模拟雷诺数范围内(Re=20 100)串联双矩 形柱绕流，研究表明：Re：2o0—300时，前柱尾流涡脱处于双剪切层控制阶段；柱间涡街为Karman类涡街， 在小间距条件下被抑制，形成涡环；前柱对后柱屏蔽效应体现为后柱阻力系数远小于前柱；临界间距时柱间涡 街充分发展，后柱阻力系数等气动参数亦在此发生跃升，但仍小于前柱值；随雷诺数升高，尾流涡街尺寸缩小， 临界间距及跃升幅度变小．Re 400时，前柱尾流涡脱进入冲击剪切层控制阶段，阻力系数不再呈现规律性 振荡；此后随雷诺数升高，冲击剪切层逐步完善，前柱流动分离使其表面产生更多附着涡，导致尾流旋涡尺寸 进一步减小，屏蔽效应消失，涡脱更为剧烈，进而对后柱产生脉动冲击效应；适当间距比条件下此类脉动冲击 效应使得后柱阻力系数发生跃升，并略高于前柱． 关键词：浸入式边界算法，串联双矩形柱，屏蔽效应，临界间距 PAcS：47．11．_j，47．27．ek，47．15．n，47．32．cd DOI：10．7498／aps．63．214702 布置临界间距与并联偏流均能较好地符合实验结 1 引 言 果．Liu等 7【】在雷诺数2×100Re1．6×104范 围内进行串联方柱绕流实验，得出串联临界间距随 物体绕流是流体力学基础研究的重中之重．同 雷诺数变化的规律．然而如上所述，目前多钝体绕 单钝体绕流相 比[1I2】，多钝体绕流流场更是 由于其 流研究仍集 中于方柱和圆柱，对于桥梁、船舶等工 相邻影响会变得更加复杂．对其研究可以直接涉及 程 中具有重要研究价值的矩形柱 串联流场研究却 到工程实际流场的本质，例如建筑结构群以及复杂 鲜有涉及． 结构的流场特性． 本文利用刚性边界浸入式数值模拟方法fVir— 较早对此进行研究的学者是Igarashi[3J’他根 tualBoundaryMethod1[rl中力源反馈物体边界 据流场特征将不同间隔距离下低雷诺数串联圆柱 的思想，用 以避开传统CFD方法 中边界网格构建 绕流划分为六个阶段．随后为了更好地理解 串联 所带来的计算耗费问题 【，驯，并相应改动Virtual 钝体流场气动干扰机理，Slaouti【]和Meneghini[5】 BoundaryMethod原有 内部流体处理方式，进一步 等分别采用离散涡方法和有限元法模拟 了二维 将算法拓展至有利于提高工程计算效率的非等间 串联双圆柱绕流．陈素琴等 6【]则基于标准网格法 距网格，从而展开雷诺数范围Re=200—100内不 (MAC)数值模拟出高雷诺数下 (Re=104)二维串 同串联间距的矩形柱绕流数值模拟，通过细化分析 联和并联两种布置方式的方柱绕流，其得出的串联 矩形柱阻力系数、升力系数、表面压力分布等气动 $国家重点基础研究发展计划 (973计划)项 目(批准号：2013CB036301)、国家自然科学基金 (批准号和土木工程防灾国 家重点实验室自主研究课题基金f批准号：SLDR