串列双圆球流场涡结构.pdfVIP

串列双圆球的流场涡结构‘ 邹建锋‘任安禄。邓见 (浙江大学力学系流体工程研究所，浙江杭州310027) 球的绕流场。研究发现，在小间距比(L／D≈I．5)的情况下，由于上下游圆球尾流区的相互抑制消除了压力不 稳定性，整个流场呈现稳态轴对称特征。间距比为2．0时，周向压力梯度诱发出流体的周向输运，流场呈 现稳态非对称性，但流场中存在特定的对称面。间距比增大到2．5后，绕流场开始周期振荡，原有的对称 面依I同存在。文中详细研究了间距比为3．0时的流场三维特征，给出了一个流动周期内的涡结构脱落情况． 发现下游圆球上表面的涡脱频率为下表面的2倍频。进而对上下游尾流场频谱特征的分析明确的给出了不 |司删距比下占优频率的竞争机制。间距比增至7．0时．整个流场恢复稳态特征，两圊球尾部同时出现双线 涡．这时流场对称面的位置发生丁变动。本文的绝大多数数值结果在浙江大学工程与科学计算中心的SGI Origin3900并行机Jj完成，在此致谢。 关键词：串列双网球绕流；涡结构；流场频谱；虚拟边界法 引言 三维圆球粘性绕流是流体力学特别是多相流动力学中的基础性问题，该问题的研究直接 涉及剑复杂流动特性的揭示。相关文献中这方面的研究集中于实验测量和数值模拟两方面。 随着大规模高性能计算机的发展，直接数值模拟在复杂流动研究中所占的比重越来越大。 对于单圆球的尾流场特性，实验测量取得了比较一致的结果，这方面的研究可参考 流动分离，产生轴对称的涡环，这里胎数由pUD／／z定义(U为进口处的均匀流速度，D 为圆球直径)；Rez210时流场经历～次常规分叉，轴对称特征消失，流场呈现平面对称性， 但依IH保持稳态，试验中尾流场显示为双线涡结构；Re在270附近时，Hopf分又出现，所 谓的发夹状涡结构开始形成并规则脱落，流场中的对称面继续存在；胎数介于420与800 之问，尾涡开始不规则脱落并向湍流转捩，涡结构在向下游输运过程中发生低频旋转，对称 面不复存在。需要注意，Taneda的早期实验在Re=130左右就观察到尾流场中存在震荡现象， 商这与其他广为认可的实验和数值计算结果有出入。 进入90年代，有关单圆球绕流数值模拟的文献出了很多。然而由于计算速度和数值方 法的限制，文献中大都采用了轴对称的假定。由实验研究结果可知，这些文献中尺P数超过 Lee～等通过求解原始的三维流动方 2lO的数值结果都是没有意义的。而Johnson一和Sungsu 程，模拟了非稳态的圆球绕流，得到了与试验结论吻合良好的结果。 至于若干圆球尾流场相互作用的研究，在文献中很少发现。80年代初期，Yutaka‘4通过 实验手段研究双圆球的相互影响，取得了一些定性结果。最近Schouveile?。对并列双圆球的 尾涡相互作用进行了细致的实验研究，提供了不同间距下双涡的耦合模式。 由于多连通区域内结构化网格生成的困难，很少发现用差分方法模拟双圆球绕流的文 ’罐金项H：国家自然科学基助 。联系人．Email地址：lucIc,tuke@zju．edu．cn 6联系人．Email地址：renantu@zju．edu．cn 137 献。KimI”利用对称性条件，通过模拟四分之一流场，数值研究了低胎数(150)时不同间距 l?并列坝圆球的尾流场。串列双圆球的数值模拟未见文献报道。 本文对Goldstein01提出的虚拟边界法作了细微改动并将之推广到三维流场的数值模拟， 称为虚拟体法。对Re=250时不同问距比L／D下串列双圆球的绕流场作了研究，给出了尾涡 结构和流场频谱随间距的变化情况。这里L为圆球中心之间的距离。 1数值方法描述 虚拟体积法的实现 Goldstein提出的虚拟边界法让流体充满绕流物所在的空间，物体的边界效应通过在流 体运动方程的右端跗加一个反馈力场加以实现。反馈力场只在表示物体边界的虚拟点上有 值，其表达式可以写成： 一一 ■一一 _一 一一 一～ F(x一。0t“I，【“rx一，砂一Vrx一，砂】卉+∥【“rx一，0-Vrx一，纠(1) ，℃(1)中，a与B是两个绝对值足够大的负常数，x，标识物面位置，u