流体绕多个钝体不稳定分离流动数值仿真.docVIP

吴文权( 华东工业大学动力工程学院)摘 要 采用离散涡方法对绕多个钝体不稳定 吴文权 ( 华东工业大学动力工程学院) 摘 要 采用离散涡方法对绕多个钝体不稳定 、分离流动进行了数值仿真 。数值 试验结果表明 ,在来流对称的条件下 ,绕双圆柱流动是双稳态的 ,这和实验是一 致的 。绕多个钝体的流动更为复杂 ,前排对后排的流动存在很大的影响 ,尤其在 前排不仅是一个物体的情况下 。计算中还包括了物体的不同截面形状 。这些计算 结果将有助于各种换热器 、海洋平台 、建筑群布置等正确设计 ,改变以往仅按平 均 、稳定流动的设计 。 关键词 绕多个钝体流动 ;不稳定流 ;分离流动 ;离散涡方法 中图法分类号 O 357 . 52 绕多个钝体流动是动力工程中广泛存在的现象 ,例如叶栅流动就是一个典型的绕多个 物体的流动 。再如各种换热器中流体绕圆柱群流动 ,海水绕海洋平台的圆柱群绕流 ,空气绕 建筑群流动等等 。因而对绕多个钝体流动的研究具有很实用的意义 。同时 ,无论内部流动或 外部流动 ,不稳定 、分离流动是当前的一个重要的研究方向 。而绕多体流动往往存在着不稳 定 、分离现象 。例如绕双圆柱的流动 ,即使在来流均匀 、完全对称的条件下 ,也存在着双稳态 流动 。这是在实验中发现的非常典型的不稳定流动 ,也是流体力学中一个经典难题 。再如我 们熟知的 Kar man 涡街 ,也是一种不稳定流动 。除此以外 ,这种流动中往往还存在着旋涡 、分 离等 ,所以流动特别复杂 。除了理论分析 、实验以外 ,对其进行数值仿真是一个重要方面 。近 年来发展的离散涡方法对此提供了一个有效的手段 。 本文中结合工程实际 ,对双圆柱 、三圆柱 、多方块等多种物体形状 ,在不同的来流方向采 用了离散涡方法进行数值实验 ,得到了一系列有意义的结果 。 1 离散涡方法 基于二维 、不可压的假定 ,以旋度ω 为框架的 N2S 方程为 Dω/ D t =υΔω 式中 υ, 为粘性扩散系数进一步将连续的旋涡场离散为一系列涡泡的集合 ,每一个涡泡的直径充分小 。在离散过 程中保持其形成的环量不变 ,有N 式中 υ, 为粘性扩散系数 进一步将连续的旋涡场离散为一系列涡泡的集合 ,每一个涡泡的直径充分小 。在离散过 程中保持其形成的环量不变 ,有 N Γ = ∫ω( x , y) d x d y = ∫ω( a , b) d ad b = ∑ω \ - i ( a , b) dσ i = 1 随后在拉格朗日坐标下 ,计算各个涡泡的运动 。当各个涡泡位置确定后 ,就可以得到由其诱 导的流体速度及流场全部信息 。这种方法的特点是在计算中不需要计算网格 ,而且自动将计 算集中在旋涡存在的区域 ,那里正是我们感兴趣的区域 ,可以得到随时间发展的非定常流 场 。有关这种方法的详细介绍见文献1 ] 。 2 绕双圆柱流动 绕圆柱流动一直是流体力学中的一个经典研究课题 。即使绕单圆柱流动 ,由于其丰富多 采 ,至今仍是一个研究热点 。我们最近的一些研究2 ] 获得一系列有趣的结果 ,例如旋涡的产 生 ,脱落 ,振荡等等 ,这些都是最基本的流动现象 。绕双圆柱的流动自然较此远为复杂 。我们 计算均匀 、零角度来流的情况 ,如图 1a 所示 。其计算 Re 数为 1 ×107 。按流动条件下看 ,这种 流动似乎应该是对称的 。其实不然 ,我们对其进行了长时间序列的计算 ,其结果如图 1a 、2a (a) 旋涡分布 图 1 ( b) 流动流谱 绕双圆柱流动 tt = 11 . 05 (无量纲时间) (a) 旋涡分布 图 2 ( b) 流动流谱 绕双圆柱流动 tt = 26 . 00 (无量纲时间) ? 1994-2013 China Academic Journal Electronic Publishing House. All rights reserved. ki. 其代表的物理时间就长 ,反之亦然 。从图 其代表的物理时间就长 ,反之亦然 。 从图 1 、2 中我们可以看到 ,流动是非定常的 ,存在着剧烈的分离 ,而且是双稳态的 。对比 图 1a 和图 2a ,可以看到在不同时刻 ,圆柱后面的旋涡处于不同的位置 。可见它是不稳定的平 衡态 ,这个结果与实验是一致的 。在每一个时刻 ,除了流谱以外 ,还给出了旋涡的分布 ,如图 1b 、2b 所示 。我们可以由流谱得出涡谱 ,反之 ,也可以由涡谱得出流谱 ,这是一个微分与积分 的关系 ,两者是统一的 。事实上 ,在离散涡方法中是先得出涡的分布 ,然后直接得出由旋涡所 诱导的流场 。对比流谱与涡谱可见 ,旋涡区域也就是我们离散的旋涡元集中的地方 。我们可 以将计算集中于最值得注意的地方 ,这正是离散涡方法