Iterative Solvers and Preconditioners for Fully-Coupled Finite Element Formulations of Incompressible Fluid Mechanics and Related Transport Problems

Finite element discretization of fully-coupled, incompressible flow problems with the classic mixed velocity-pressure interpolation produces matrix systems that render the best and most robust iterative solvers and preconditioners ineffective. The indefinite nature of the discretized continuity equation is the root cause and is one reason for the advancement of pressure penalty formulations, least-squares pressure stabilization techniques, and pressure projection methods. These alternatives have served as admirable expedients and have enabled routine use of iterative matrix solution techniques; but all remain plagued by exceedingly slow convergence in the corresponding nonlinear problem, lack of robustness, or limited range of accuracy. The purpose of this paper is to revisit matrix systems produced by this old mixed velocity-pressure formulation with two approaches: (1) deploying well-established tools consisting of matrix system reordering, GMRES, and ILU preconditioning on modern architectures with substantial distributed or shared memory, and (2) tuning the preconditioner by managing the condition number using knowledge of the physical causes leading to the large condition number. Results obtained thus far using these simple techniques are very encouraging when measured against the reliability (not efficiency) of a direct matrix solver. Here we demonstrate routine solution for an incompressible flow problem using the Galerkin finite element method, Newton-Raphson iteration, and the robust and accurate LBB element. We also critique via an historical survey the limitations of pressure-stabilization strategies and all other commonly used alternatives to the mixed formulation for acceleration of iterative solver convergence. The performance of the new iterative solver approaches on other classes of problems, including fluid-structural interaction, multi-mode viscoelasticity, and free surface flow is also demonstrated.

[1]  R. Codina,et al.  A finite element formulation for the Stokes problem allowing equal velocity-pressure interpolation , 1997 .

[2]  O. C. Zienkiewicz,et al.  Constrained variational principles and penalty function methods in finite element analysis , 1974 .

[3]  Yitzhak Hasbani,et al.  Out-of-core solution of linear equations with non-symmetric coefficient matrix , 1979 .

[4]  Michele Benzi,et al.  A Sparse Approximate Inverse Preconditioner for Nonsymmetric Linear Systems , 1998, SIAM J. Sci. Comput..

[5]  D. Knoll,et al.  A NEW NONLINEAR SOLUTION METHOD FOR PHASE-CHANGE PROBLEMS , 1999 .

[6]  H. Elman,et al.  Efficient preconditioning of the linearized Navier-Stokes , 1999 .

[7]  Graham F. Carey,et al.  Performance of iterative methods for Newtonian and generalized Newtonian flows , 1989 .

[8]  Edmond Chow,et al.  An object-oriented framework for block preconditioning , 1998, TOMS.

[9]  J. Reddy,et al.  The Finite Element Method in Heat Transfer and Fluid Dynamics , 1994 .

[10]  Jun Zhang,et al.  Preconditioned Krylov subspace methods for solving nonsymmetric matrices from CFD applications , 2000 .

[11]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[12]  H. V. D. Vorst,et al.  Iterative solution methods for certain sparse linear systems with a non-symmetric matrix arising from PDE-problems☆ , 1981 .

[13]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[14]  Thomas J. R. Hughes,et al.  The Stokes problem with various well-posed boundary conditions - Symmetric formulations that converge for all velocity/pressure spaces , 1987 .

[15]  Jeffrey J. Derby,et al.  Massively parallel finite element computations of three-dimensional, time-dependent, incompressible flows in materials processing systems , 1994 .

[16]  T. Manteuffel,et al.  Necessary and Sufficient Conditions for the Existence of a Conjugate Gradient Method , 1984 .

[17]  P. A. Sackinger,et al.  A finite element method for free surface flows of incompressible fluids in three dimensions. Part I. Boundary fitted mesh motion , 2000 .

[18]  F. Baaijens,et al.  An iterative solver for the DEVSS/DG method with application to smooth and non-smooth flows of the upper convected Maxwell fluid , 1998 .

[19]  K. Jensen,et al.  A finite element solutionof three‐dimensional mixed convection gas flows in horizontal chnnels using preconditioned iterative metrix methods , 1992 .

[20]  YereminA. Yu.,et al.  Factorized sparse approximate inverse preconditionings I , 1993 .

[21]  E. Cuthill,et al.  Reducing the bandwidth of sparse symmetric matrices , 1969, ACM '69.

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

[23]  T. Manteuffel.  An incomplete factorization technique for positive definite linear systems , 1980 .

[24]  P. Hood,et al.  Frontal solution program for unsymmetric matrices , 1976 .

[25]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[26]  Yousef Saad,et al.  Preconditioned Krylov Subspace Methods For CFD Applications , 1995 .

[27]  Ray S. Tuminaro,et al.  Parallel Smoothed Aggregation Multigrid : Aggregation Strategies on Massively Parallel Machines , 2000, ACM/IEEE SC 2000 Conference (SC'00).

[28]  Michel Fortin,et al.  A preconditioned generalized minimal residual algorithm for the numerical solution of viscoelastic fluid flows , 1990 .

[29]  Michael S. Engelman,et al.  Segregated finite element algorithms for the numerical solution of large‐scale incompressible flow problems , 1993 .

[30]  C. F. Curtiss,et al.  Dynamics of Polymeric Liquids, Volume 1: Fluid Mechanics , 1987 .

[31]  F. Baaijens,et al.  Application of low-order Discontinuous Galerkin methods to the analysis of viscoelastic flows , 1994 .

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

[33]  Comparison of three solvers for viscoelastic fluid flow problems , 1995 .

[34]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[35]  Unsymmetric conjugate gradient methods and sparse direct methods in finite element flow simulation , 1990 .

[36]  J. A. George.  Computer implementation of the finite element method , 1971 .

[37]  Homer F. Walker,et al.  An Inexact Newton Method for Fully Coupled Solution of the Navier-Stokes Equations with Heat and Mass Transport , 1997 .

[38]  Thomas J. R. Hughes,et al.  A boundary integral modification of the Galerkin least squares formulation for the Stokes problem , 1994 .

[39]  P. Gresho.  Some current CFD issues relevant to the incompressible Navier-Stokes equations , 1991 .

[40]  M. Fortin,et al.  Some special purpose preconditioners for conjugate gradient‐like methods applied to CFD , 1993 .

[41]  Yousef Saad,et al.  High-order ILU preconditioners for CFD problems , 2000 .

[42]  Y. Saad.  Preconditioning techniques for nonsymmetric and indefinite linear systems , 1988 .

[43]  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..

[44]  Yousef Saad,et al.  Robust preconditioning for sparse linear systems , 1997 .

[45]  Rekha Ranjana Rao,et al.  A Newton-Raphson Pseudo-Solid Domain Mapping Technique for Free and Moving Boundary Problems , 1996 .

[46]  M. Fortin,et al.  A new mixed finite element method for computing viscoelastic flows , 1995 .

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

[48]  Abdelmalek Zine,et al.  An improved GMRES method for solving viscoelastic fluid flow problems , 1992 .

[49]  Antonio Huerta,et al.  Viscous flow with large free surface motion , 1988 .