用预处理方法求解全马赫数范围内的无粘/粘性流动

摘要

用时间相关法、LU-SGS隐式求解经过预处理后的二维可压Euler/Navier-Stokes方程，数值模拟了全赫数范围内直管道下壁面带凸包(Bump)的无粘流动和低速情况下不同参考马赫数、不同雷诺数的粘性空腔流动。为了使时间相关法能应用于计算低速不可压流动及加快收敛，对控制方程进行了时间导数预处理。数值模拟结果表明该文选取的预处理矩阵在低马赫数时大大加快了计算的收敛性，计算有较高的效率。所得结果与可能的计算结果进行了比较。