Transient solutions for three‐dimensional lid‐driven cavity flows by a least‐squares finite element method

SUMMARY A time-accurate least-squares finite element method is used to simulate three-dimensional flows in a cubic cavity with a uniform moving top. The time- accurate solutions are obtained by the Crank-Nicolson method for time integration and Newton linearization for the convective terms with extensive linearization steps. A matrix-free algorithm of the Jacobi conjugate gradient method is used to solve the symmetric, positive definite linear system of equations. To show that the least-squares finite element method with the Jacobi conjugate gradient technique has promising potential to provide implicit, fully coupled and time-accurate solutions to large-scale three-dimensional fluid flows, we present results for three-dimensional lid-driven flows in a cubic cavity for Reynolds numbers up to 3200.

[1]  Giuseppe Gambolati,et al.  Is a simple diagonal scaling the best preconditioner for conjugate gradients on supercomputers , 1990 .

[2]  Robert L. Lee,et al.  A MODIFIED FINITE ELEMENT METHOD FOR SOLVING THE TIME-DEPENDENT, INCOMPRESSIBLE NAVIER-STOKES EQUATIONS. PART 1: THEORY* , 1984 .

[3]  T. Taylor,et al.  A Pseudospectral method for solution of the three-dimensional incompressible Navier-Stokes equations , 1987 .

[4]  Jae Min Hyun,et al.  Numerical simulation of three-dimensional flow structure in a driven cavity , 1989 .

[5]  Wagdi G. Habashi,et al.  Finite Element Solution of the 3D Compressible Navier-Stokes Equations by a Velocity-Vorticity Method , 1993 .

[6]  T. H. Tsang,et al.  AN EFFICIENT LEAST-SQUARES FINITE ELEMENT METHOD FOR INCOMPRESSIBLE FLOWS AND TRANSPORT PROCESSES , 1995 .

[7]  A pseudospectral matrix element method for solution of three-dimensional incompressible flows and its implementation on a parallel computer , 1989 .

[8]  H. Kawai,et al.  Numerical Flow Analysis in a Cubic Cavity by the GSMAC Finite-Element Method : In the Case Where the Reynolds Numbers are 1000 and 3200 , 1990 .

[9]  O. C. Zienkiewicz,et al.  Least square-finite element for elasto-static problems. Use of `reduced' integration , 1974 .

[10]  Fujima Shoichi,et al.  Extension to three-dimensional problems of the upwind finite element scheme based on the choice of up- and downwind points , 1994 .

[11]  Robert L. Street,et al.  Numerical simulation of three‐dimensional flow in a cavity , 1985 .

[12]  Li Q. Tang,et al.  A least‐squares finite element method for time‐dependent incompressible flows with thermal convection , 1993 .

[13]  Louis A. Povinelli,et al.  Optimal least-squares finite element method for elliptic problems , 1993 .

[14]  F. Stella,et al.  A vorticity-velocity method for the numerical solution of 3D incompressible flows , 1993 .

[15]  Robert L. Street,et al.  Three‐dimensional unsteady flow simulations: Alternative strategies for a volume‐averaged calculation , 1989 .

[16]  R. B. Kellogg,et al.  Least Squares Methods for Elliptic Systems , 1985 .

[17]  G. de Vahl Davis,et al.  An evaluation of upwind and central difference approximations by a study of recirculating flow , 1976 .

[18]  Louis A. Povinelli,et al.  Large-scale computation of incompressible viscous flow by least-squares finite element method , 1994 .

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

[20]  Urmila Ghia,et al.  A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations , 1987 .

[21]  Dennis C. Jespersen,et al.  A least squares decomposition method for solving elliptic equations , 1977 .

[22]  P. P. Lynn,et al.  Use of the least squares criterion in the finite element formulation , 1973 .

[23]  B. Jiang.  A least‐squares finite element method for incompressible Navier‐Stokes problems , 1992 .