摘要:为克服传统有限元方法求解流场时由于对流项占优而引起的求解振荡和基函数选择困难的问题,推导了基于特征线方程的粘性不可压N-S方程的增量型和非增量型分离算法的公式及求解步骤,并讨论了两者的Babuska-Brezzi条件.沿着特征线方向,N-S方程的对流项消失,方程矩阵是自伴随矩阵,可以自动满足有限单元法中能量泛函最小的要求,并可以给出合理的粘性耗散项.动量方程的求解采用非增量型分离算法,压力和速度可以采用任意阶次的插值函数,离散后的方程自动满足Babuska-Brezzi条件.为验证算法的可靠性,采用T3P3空间等阶次U-p单元计算了方柱绕流.结果表明:基于特征线方程的非增量型分离算法可以很好地应用于粘性不可压流场的计算,与目前流行的其它方法相比,该算法有明显的优势.
关键词:分离算法 计算流体力学 有限元方法 特征线方程 方柱绕流
单位:清华大学热能工程系水沙科学与水利水电工程国家重点实验室 北京100084
注:因版权方要求,不能公开全文,如需全文,请咨询杂志社