横向排列双圆柱绕流的临界Re数和St数跳跃研究.pdfVIP

横向排列双圆柱绕流的临界Re数和St数跳跃· 刘坤，马东军，孙德军，尹协远 中围科学技术人学，力学和机械T程系，合肥230027 小扰动方程，求得系统的优势特征值及其特征向量。这种流动稳定性的数值分析方法，町以处理复杂流动，具自． 商的精度和实用性。本文将此方法应用于研究横向排列双厕柱绕流的整体稳定性．分析J，流动的分甜情7兕，得到 了临界Reynolds数和频率(St数)随无量纲圆心距T／D变化曲线。 关键词：横向排列双圆柱，谱元法，整体稳定性分析，临界Reynolds数 识，同类研究目前尚很少见剑报道J 0EJI言 1 流动稳定性的数值计算方法 对丁流向排列和横向排列舣圆柱绕流的研究历 史可以追溯到很早的年代，且研究势头始终不衰， 对于较为简单的流动，如平行流，用LE则模法 ofFluids和J．Fluids 仅近五年以来在Physics 进行线性稳定性分析是通行的做法。但对于复杂流 Structures等杂忠上发表的论文，不完生统计就有J0动，正则模法就显现出它的局限性。目前，利用数 篇之多【1．7】，有兴趣的读者可以从中找到更多的实值方法研究流动稳定性日渐成为流行的做法。它可 验和数值研究的文献。大多数已有的研究是在中等 以．}}j于求出非平行流的整体稳定性特征。例如，孙 Re数(100~1000)下进行的，研究者主要关心的是：(1) 流动图像及其分类，(2)升阻力系数特性．(3)烈稳匀来流中静止和旋转圆柱的整体稳定性，马乃i军 态和滞迥现象等。其计算结果与实验符合很好，并 【10，ll】利用潜元法计算了旋转圆柱绕流和方腔流的 表明Re=1000以F的数值计算就足以重现实验中更 ‘ 高Re得剑的流动现象。 丁精度较低，只能得剑定性止确、定链合理的结果。 最近，Kang(2003)[3]求解了低Re数(40≤Rc≤l60)最近发展的基了：谱元法的稳定性数值分析法可以得 横向排列情况下双杠的流动。他把尾迹图缘分成6 剑很高的精度。而基本DNS的烛非线性的整体稳定 类：反相同步开．J．(antiphase．synchronized)，同相位同步性的数值分析方法和周期基本流的Floquet稳,定性分 析也j!：企迅速发展中。 型tin。phase．synchronized)，暇稳触发型(flip． 理论上讲，任恿时间推进格式剪法都可以进行 flopping)t斜偏犁(deflected)，单钝体型tsingle wake 流动分义分析而不筒进彳j：额外修改：首宄通过计舜 bluff-body),tq!定常迹蚕-．(steadypauems)，这些类 得钊稳定定常解，再逐渐改变控制参数，n剑得到 !i：．J．LjRe数和间隙I：E(L／D)宵关，并认为对后者的依赖 性更人。 流动模式发生变化，此时的控制参数对麻’r流动分 MizushimaSuehiro(2005)[2]对流向排列双圆岔点。但这种处理方法收敛速度低，会耗费人城时 托J【{j数值力‘法进仃了粘体线陛稳定性分折，给出了间和硬件资源．而H时问推进格式算法无法得剑不 流动分近情况，解释了中等削甄旧寸流动的删稳现象。 稳定的定常解及其扰动增K率，所得临界参数不够 准确。 其后，TasakaKon铬人(20·16)11I实验证实了12i的 结粟。 州数值方法{17f究流动稳定性问题，其土要思想 本文的目的楚}lj数值方法研究横向排列舣测枉是：直接数值求解出稳定或不稳定的定常基本流，