求解不可压缩Navier-Stokes方程的GLS平面八节点单元
The Planar Element with Eight Nodes for Solving Incompressible Navier-Stokes Equations
DOI: 10.12677/IJFD.2013.12003, PDF, HTML, 下载: 2,621  浏览: 10,115 
作者: 魏晓娟:四川电力职业技术学院动力工程系,成都;刘力菱:西南交通大学材料科学与工程学院,成都;谢凌志, 易丽清, 魏泳涛*:四川大学建筑与环境学院应用力学系,成都
关键词: Galerkin/最小二乘有限元稳定因子单元长度逆估计常数Galerkin/Least Squares FEM; Stabilization Factor; Element Length; Inverse Estimate Constant
摘要: 给出了求解不可压缩Navier-Stokes方程的GLS平面八节点单元的“单元长度”的定义和相应的逆估计常数,由此可计算出该单元的GLS稳定因子。给出的“单元长度”只依赖单元形状且易于计算;各种形状的直边四边形单元的逆估计常数的最大值由遗传算法得出。对Reynolds数为20,000的方腔上盖板流的数值模拟表明,本文所提出的单元长度和逆估计常数的正确性。
Abstract: The definition of the element length and the inverse estimate constant of straight-edged quadrilateral element with eight nodes for solving incompressible Navier-Stokes equations was proposed, based on this, the stabilization factor was readily calculated. The element length was only dependent on the quadrilateral geometry and easy to compute, and the largest inverse estimate constant among various quadrilateral was obtained by means of the optimization method of genetic algorithm (GA). The numerical solutions of the classical lid-driven cavity flow problem were ob- tained for Reynolds number of 20,000 which demonstrated the validity of the element length and the inverse estimate constant proposed.
文章引用:魏晓娟, 刘力菱, 谢凌志, 易丽清, 魏泳涛. 求解不可压缩Navier-Stokes方程的GLS平面八节点单元[J]. 流体动力学, 2013, 1(2): 15-20. http://dx.doi.org/10.12677/IJFD.2013.12003

参考文献

[1] 1. Babuska I. Error bounds for finite element method. Numerische Mathematik 1971; 16:322–333.
[2] 2. Brezzu F. On the existence, uniqueness, and approximation of saddle–point problems arising from Lagrange multiplier. RAIRO B-R2, 1974; 129–151.
[3] 3. Brooks AN, Hughes TJR. Streamline upwind/Petrov–Galerkin formulation for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations. Computer Methods in Applied Mechanics and Engineering 1982; 32:199–249.
[4] 4. Franca LP, Frey SL. Stabilized finite element methods: II. The incompressible Navier–Stokes equation. Computer Methods in Applied Mechanics and Engineering 1992; 99:209–233.
[5] 5. Burda P, Novotny J, Sistek J. On a modification of GLS stabilized FEM for solving incompressible viscous flows. International Journal for Numerical Methods in Fluids 2006; 51:1001–1016
[6] 6. Harari I, Hughes TJR. What are C and h: inequalities for the analysis and design of finite element methods.Computer Methods in Applied Mechanics and Engineering 1992; 97:157–192.
[7] 7. Franca LP, Madureria AL. Element diameter free stability parameters for stabilized methods applied to fluids. Computer Methods in Applied Mechanics and Engineering 1993; 105:395–403.
[8] 8.Tezduyar TE, Osawa Y. Finite element stabilization parameters computed from element matrices and vectors. Computer Methods in Applied Mechanics and Engineering 2000; 190:411–430.
[9] 9. Erturk E, Corke TC, Gokco1 C. Numerical solutions of 2-D steady incompressible driven cavity flow at high Reynolds numbers. International Journal for Numerical Methods in Fluids 2005; 48:747–774.
[10] 10. Y.T. Wei, P.H.Geubelle, A comparative study of GLS finite element with velocity and pressure equally interpolated for solving incompressible viscous flows. International Journal for Numerical Methods in Fluids. 2009; 61:514–535
[11] 11.Wei YT, Yu JH, Chen JK. Object-oriented approach to the finite element programming: basic data classes.Journal of Sichuan University (Engineering Science Edition) 2001; 33(2):17–21.[魏泳涛、于建华、陈君楷,面向对象有限元程序设计—基本数据类 [J],四川大学学报:工程科学版,2001,33(2):1721]
[12] 12.Wei YT, Yu JH, Chen JK. Object-oriented approach to the finite element programming: design of element procedure. Journal of Sichuan University (Engineering Science Edition) 2001; 33(3): 9–12. [魏泳涛、于建华、陈君楷,面向对象有限元程序设计—单元过程设计 [J],四川大学学报:工程科学版,2001,33(3):912]
[13] 13. Wei YT, Yu JH, Chen JK. Object-oriented approach to the finite element programming:the applicationarchitecture. Journal of Sichuan University (Engineering Science Edition) 2001; 33(4): 21–25. [魏泳涛、于建华、陈君楷,面向对象有限元程序设计—程序构架 [J],四川大学学报:工程科学版,2001,33(4):2125]