束解非定常方柱绕流的可变周期时间谱方法

摘要

针对周期性非定常流动特点,采用时间谱方法求解了无限大区域内低雷诺数下二维不可压缩层流绕方柱流动问题.由于方柱绕流的非定常运动周期事先未知,在计算过程中加入了基于残差梯度的可变周期计算方法,使变量和周期可以同时迭代计算.研究了时间谱方法的计算精度和求解速度,分析了不同求解参数对最终结果的影响.计算结果表明,时间谱方法可以准确地模拟方柱后漩涡的非定常运动现象,而计算所需时间只是传统的基于有限差分的双时间推进方法的10%～25%,流体非定常运动的Strouhal数和方柱的阻力特性证实了这种方法的正确性.可变周期计算方法对求解参数的变化不敏感,可以准确地预测运动周期.