NS方法数值预测跨音速气动颤振边界(英文)

本文通过耦合求解三维薄层Navier-Stokes方程与结构运动方程数值模拟了跨音速气动颤振现象。其中,流动控制方程的求解采用具有三阶精度的HLLEW(Harten-Lax-vanLeer-EinfeldtWada)迎风TVD空间离散格式和LU-SGS内迭代时间推进方法,同样的内迭代方法用于结构运动方程的求解。在每一步实时推进计算中,通过内迭代,使整个耦合计算的时间精度达到二阶。针对每一时间步的结构变形,发展了一种自适应网格变形方法,在中等结构变形的情况下,该方法能保证变形后的网格具有原网格的质量。为检验发展的跨音速气动颤振计算程序,对一标准气动弹性机翼的跨音速气动弹性边界进行了计算,获得了与实验一致的结果。另外,还详细研究了网格数、时间步长及内迭代步数对气动颤振计算的影响。