摘要

谐波平衡法是一种高效周期性非定常流动频域计算方法.本文研究可变周期谐波平衡法,通过求解Navier-Stokes方程模拟低速不可压条件下的二维钝体绕流周期性非定常涡脱落问题.对于这类流动变化周期未知的非定常问题,将涡脱落周期T作为变量,采用基于残差导数的可变周期计算方法推进求解.以圆柱绕流和方柱绕流为例,研究考察了谐波平衡法的计算精度和效率,并分析研究了不同参数对计算结果的影响.针对圆柱绕流问题,采用三种不同优化方法进行周期T的迭代计算,对比研究了它们的计算精度和效率.计算结果表明:谐波平衡法采用3个谐波数就可以准确模拟周期性非定常涡脱落问题,辨识的涡脱落频率和阻力系数与实验值及其他数值结果一致,与时域方法相比该方法具有较高的计算效率.不同优化方法的计算结果相同,共轭梯度法和牛顿法的收敛速度与最速下降法采用较大搜索步长时的收敛速度一致.由于牛顿法没有参数问题,因此该方法在工程计算中更有优势.