最小二乘有限元法和有限体积法在CFD中应用比较.docVIP

最小二乘有限元法和有限体积法在CFD中应用比较 　　摘要： 为比较最小二乘有限元法（Least Square Finite Element Method，LSFEM）和有限体积法在CFD应用中的优劣，采用最小二乘法离散不可压N-S方程的有限元模型，得到正定对称线性系统，采用高效的预处理共轭梯度法求解方程组；利用LSFEM和基于有限体积法的FLUENT分别计算Kovasznay流动、定常二维和三维后台阶流动以及非定常圆柱绕流等4个实例并比较计算结果. 结果表明，LSFEM比有限体积法的收敛性和精确性更好，在CFD领域的应用价值很高. 　　关键词： 　　最小二乘有限元法； 有限体积法； CFD； 预处理共轭梯度法； Kovasznay流动； 后台阶?流动； 圆柱绕流 　　中图分类号： O35；TB115.1 　　 文献标志码： ?A 　　 　　Application comparison of least square finite element method and finite volume method in CFD 　　TAO Sha?1, YANG Zhigang?1, JIANG Bonan?2, GU Wenjun?1 　　(1. Shanghai Automotive Wind Tunnel Center, Tongji University, Shanghai 201804, China; 　　 2. Department of Mathematics and Statistics, Oakland University, Miami 48309, USA) 　　Abstract： To compare the advantages and disadvantages of Least Square Finite Element Method(LSFEM) and finite volume method in actural applications, a least square method is used to discrete the finite element model of incompressible N-S equation, the positive definite symmetric linear system is obtained, and the equations are solved by efficient preconditioned conjugate gradient method; the least square finite element method and FLUENT which is based on finite volume method are separately used to calculate the examples such as Kovasznay flow, steady 2D and 3D backward-facing step flows, and unsteady flow around cylinder. The comparison results indicate that the least square finite element method performs better convergence rate and accuracy than finite volume method and shows good application value in CFD. 　　Key words： 　　least square finite element method; finite volume method; CFD; preconditioned conjugate gradient method; Kovasznay flow; backward-facing step flow; flow around cylinder 　　 　　0 引 言 　　近30年来，最小二乘有限元法（Least Square Finite Element Method，LSFEM）在流体动力学中的应用得到越来越多的关注，其引入残量最小化概念求解N-S方程，具有很强的稳定性，不需要验证LBB条件［1］，??更灵活地选择有限元计算域.此外，LSFEM求解各种条件下N-S方程的离散形式都正定对称，有利于选择高效、统一的求解方法而不需要任何迎风格式和调节参数，使得基于LSFEM的计算代码具有很强的通用性.［2］ 　　目前，大多数针对LSFEM的研究主要侧重于可行性探索和论证［2-8］，尚未见对其实际