多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。

跨声速条件下轴对称收扩喷管内外流场的数值研究

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用不同湍流模型、壁面函数对轴对称喷管内外流场进行三维数值模拟.计算结果与实验数据对比表明,采用涡粘修正的RNG模型,内流场使用增强壁面函数,外流场使用标准壁面函数可以得到准确的计算结果.研究表明:在亚声速条件下,喷管外壁面压力分布随落压比的变化不大;超声速条件下,喷管外壁面激波随着落压比的增大前移.

二维定常湍流计算中的GMRES算法

在以前工作的基础上将广义极小残差（Ｇｅｎｅｒａｌｉｚｅｄ Ｍｉｎｉｍｕｍ ＲＥＳｉｄｕａｌ（ＧＭＲＥＳ））算法发展到用于求解二维可压 Ｆａｖｉｅｒ平均 Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组．控制方程经 Ｎｅｗｔｏｎ线化处理后构成近似的线性系统，然后采用分别耦合了ＬＵＳＧＳ和ＩＬＵ两种预处理矩阵的ＧＭＲＥＳ算法求解．Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ湍流模型被用来封闭流体控制方程组，采用与流体控制方程非耦合的方式，使用ＬＵＳＧＳ方法求解．对ＧＭＲＥＳ算法中矩阵－向量的乘积采用了有限插分方法，从而避免了精确的左端系数矩阵的计算和存储．对预处理矩阵的两种使用方法（左预处理和右预处理）进行了分析和讨论．用两个算例对ＬＵＳＧＳ和ＩＬＵ两种预处理矩阵进行了比较，同时探讨了左预处理和右预处理各自的优缺点．通过对Ｓａｊｂｅｎ扩压器和ＮＡＣＡ００１２有攻角流动的计算，表明带有预处理的ＧＭＲＥＳ算法在二维定常跨音黏性流动计算中相比于得到广泛应用的ＤＤＡＤＩ方法具有很大优势，左预处理要优于右预处理．