An efficient, high-order finite element method using the nodal averaging technique for incompressible fluid flows

Abstract A new finite element method is presented for use of quadrilateral nine-node elements in the solution of the incompressible Navier–Stokes equations. In a conventional predictor–corrector scheme, the method applies the nodal averaging technique to discretize the Poisson equation used for the simultaneous relaxation of velocity and pressure. Additionally, efficient approximation procedures are devised to increase the speed of computation without deteriorating solution accuracy. The proposed numerical schemes are evaluated on two-dimensional test problems including a classical lid-driven cavity flow and a flow over a backward-facing step in a flow channel. The results show good accuracy even when distorted elements are used for calculation.

[1]  László Könözsy,et al.  A unified fractional-step, artificial compressibility and pressure-projection formulation for solving the incompressible Navier-Stokes equations , 2014 .

[2]  K. Bathe.  Finite Element Procedures , 1995 .

[3]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[4]  C. W. Hirt,et al.  Calculating three-dimensional flows around structures and over rough terrain☆ , 1972 .

[5]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[6]  Klaus-Jürgen Bathe,et al.  A flow‐condition‐based interpolation finite element procedure for triangular grids , 2006 .

[7]  W. Rider,et al.  High-Resolution Methods for Incompressible and Low-Speed Flows , 2004 .

[8]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[9]  Evgeniy Shapiro,et al.  Artificial compressibility, characteristics-based schemes for variable density, incompressible, multi-species flows. Part I. Derivation of different formulations and constant density limit , 2005 .

[10]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[11]  T. Tanahashi,et al.  GSMAC FINITE ELEMENT METHOD FOR UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS AT HIGH REYNOLDS NUMBERS , 1990 .

[12]  Evgeniy Shapiro,et al.  Artificial compressibility, characteristics-based schemes for variable-density, incompressible, multispecies flows: Part II. Multigrid implementation and numerical tests , 2005 .

[13]  O. C. Zienkiewicz,et al.  The Finite Element Method for Fluid Dynamics , 2005 .

[14]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[15]  A. A. Amsden,et al.  The SMAC method: A numerical technique for calculating incompressible fluid flow , 1970 .

[16]  T. Tanahashi,et al.  Short Communication Finite Element Analysis of Visco-elastic Flow Under High Shear Rate Using the GSMAC-method , 2000 .

[17]  D. Drikakis,et al.  AUnified Fractional-Step , Artificial Compressibility and Pressure-Projection Formulation for Solving the Incompressible Navier-Stokes Equations , 2014 .

[18]  Takahiko Tanahashi,et al.  Numerical analysis of moving interfaces using a level set method coupled with adaptive mesh refinement , 2004 .

[19]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[20]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[21]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .