An element-by-element BICGSTAB iterative method for three-dimensional steady Navier-Stokes equations

Construction of a stabilized Galerkin upwind finite element model for steady and incompressible Navier-Stokes equations in three dimensions is the main theme of this study. In the time-independent context, the weighted residuals statement is kept biased in favor of the upstream flow direction by adding an artificial damping term of physical plausibility to the Galerkin framework. This upwind approach has significant advantage of seeking solutions free from cross-stream diffusion error. Finite element solutions have been found by mixed formulation, implemented in quadratic cubic elements which are characterized as possessing the so-called LBB (Ladyzhenskaya-Babuska-Brezzi) condition. An element-by-element BICGSTAB solution solver is intended to alleviate difficulties regarding the asymmetry and indefiniteness arising from the use of a mixed formulation for incompressible fluid flows. The developed three-dimensional finite element code is first rectified by solving a problem amenable to analytic solution. A well-known lid-driven cavity flow problem in a cubical cavity is also studied.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  R. Fletcher.  Conjugate gradient methods for indefinite systems , 1976 .

[3]  A comparison study on multivariant and univariant finite elements for three‐dimensional incompressible viscous flows , 1995 .

[4]  Wing Kam Liu,et al.  Finite Element Analysis of Incompressible Viscous Flows by the Penalty Function Formulation , 1979 .

[5]  J. Douglas,et al.  Stabilized mixed methods for the Stokes problem , 1988 .

[6]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[7]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[8]  I. Babuska.  Error-bounds for finite element method , 1971 .

[9]  Max Gunzburger,et al.  Finite Element Methods for Viscous Incompressible Flows: A Guide to Theory, Practice, and Algorithms , 1989 .

[10]  R. Sani,et al.  On pressure boundary conditions for the incompressible Navier‐Stokes equations , 1987 .

[11]  Roland W. Freund,et al.  A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems , 1993, SIAM J. Sci. Comput..

[12]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[13]  J. Donea.  Generalized Galerkin Methods for Convection Dominated Transport Phenomena , 1991 .

[14]  V. Babu,et al.  Numerical solution of the incompressible three-dimensional Navier-Stokes equations , 1994 .

[15]  Lloyd N. Trefethen,et al.  How Fast are Nonsymmetric Matrix Iterations? , 1992, SIAM J. Matrix Anal. Appl..

[16]  R. A. Silverman,et al.  The Mathematical Theory of Viscous Incompressible Flow , 1972 .

[17]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[18]  P. Sonneveld.  CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear systems , 1989 .

[19]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .