双圆柱绕流诱发振动的数值模拟(PartⅠ横向振动)

采用任意拉格朗日-欧拉(ALE)方法数值模拟圆柱在尾流中的流体诱发振动特性.重点分析了圆柱的动力学响应特性,包括升阻力、位移振幅、拍和锁定等现象;另外也详细分析了圆柱的尾涡结构.研究结果表明,在串列下,大质量比时,圆柱的振动会受到抑制,小质量比时,圆柱的振动则会被放大;较之孤立柱的情形,圆柱共振发生在低于且接近于1.0的频率比带内且不易随质量比的变化发生偏移;小质量比下圆柱振幅及共振带都要比大质量比下的大得多;不同频率比和间距比下,圆柱的动力学特性存在明显的差异,相应的涡脱落模态呈现出2P、P S和2S模态甚至是2P S模态,各种模态之间互相竞争促进流固耦合的不断变化发展,导致涡间距和涡街宽度的变化.     

用高精度紧致差分格式分块耦合求解二维粘性不可压缩复杂流场

本文提出了张量形式的粘性不可压缩流场方程的高精度分块耦合求解方法。在分块算法中用重叠一层的方法来实现子块间的数据交换，同时克服了块剖分所造成的数值奇性。利用本文所提出的算法成功地模拟了椭圆柱绕流问题的时间发展过程，如尾流，旋涡的产生，脱落过程以及Karman涡街的形成等。     

ALE方法求解圆柱的涡致振动

数值求解基于拉格朗日-欧拉(ALE)描述的不可压缩流体的N-S方程,计算了较小雷诺数下圆柱的涡致振动现象.N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散,计算网格采用H-O非交错网格系统,并结合分块耦合方法.柱体运动简化为弹簧-阻尼-质量系统,柱体运动方程采用经典龙格-库塔方法求解.通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"锁定"、"拍"和"相位开关"等现象,并与试验数据相吻合.另外,本文详细分析了Re＝200时圆柱尾涡形态、升力、阻力以及圆柱位移等随圆柱自振频率变化的过程,捕捉到尾涡结构变化的频率转折点.     

复杂回旋流场的数值模拟和实验研究

本文通过四角喷射炉膛内回旋流场的数值计算和实验测试的对比研究，比较了几种计算方法的预估结果与实验结果的符合程序。结果表明：为克服当流动方向与网格线交角过大时所产和珠过分数值扩散，引入流向平衡粘性方法与κ－ε湍流模型相结合可获得与实验结果相一致的计算结果，而采用混合分格式加κ－ε模型所得到的计算结果比流向平衡粘性方法加等效常湍流粘性系数这样的简化处理的计算结果还要差。