Spectral/hp penalty least-squares finite element formulation for the steady incompressible Navier-Stokes equations

In this paper, we present spectral/hp penalty least-squares finite element formulation for the numerical solution of steady incompressible Navier-Stokes equations. The continuity equation is treated as a constraint on velocity field and the constraint is enforced using the penalty method, which is equivalent to replacing the pressure variable p with p^n=p^n^-^[email protected](@?.u), where u is the velocity field, @c is the penalty parameter and superscript 'n' is non-linear iteration number. The modified Navier-Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity-dilatation or stresses as additional independent variables and the least-squares method is used to develop the finite element model. High-order basis functions are used to construct the discrete model. Unlike other penalty finite element formulations, equal order integration is used for all terms of the coefficient matrix. Schur complement method is used and all interior element degrees of freedom are condensed out and interface degrees of freedom are solved either by direct solver (banded Cholesky factorization) or iterative solver (preconditioned conjugate gradient). Present formulation produces very accurate results for even very low penalty parameters (10-40). In this way present formulation (and implementation) fixes perennial problem of bad conditioning of coefficient matrix associated with penalty formulations. Proposed formulation is verified using the exact solution of Kovasznay flow. Numerical results are also presented for a number of two-dimensional benchmark problems, e.g., 2D lid-driven cavity flow, flow over a backward-facing step, steady flow past a circular cylinder, and 2D flow past a large circular cylinder in a channel. Finally, the proposed formulation is extended to coupled velocity-temperature problems, and buoyant flow inside a square enclosure is solved to test the penalty least-squares finite element model developed herein.

[1]  John J. Nelson,et al.  Least-Squares Finite Element Method for the Stokes Problem with Zero Residual of Mass Conservation , 1997 .

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

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

[4]  D. Gartling.  A test problem for outflow boundary conditions—flow over a backward-facing step , 1990 .

[5]  B. Jiang.  The Least-Squares Finite Element Method , 1998 .

[6]  Mitutosi Kawaguti,et al.  Numerical Study of a Viscous Fluid Flow past a Circular Cylinder , 1966 .

[7]  A. S. Grove,et al.  An experimental investigation of the steady separated flow past a circular cylinder , 1964, Journal of Fluid Mechanics.

[8]  L. Kovasznay.  Laminar flow behind a two-dimensional grid , 1948 .

[9]  Jan S. Hesthaven,et al.  A spectral multidomain penalty method model for the simulation of high Reynolds number localized incompressible stratified turbulence , 2005 .

[10]  J. Reddy.  An introduction to nonlinear finite element analysis , 2004 .

[11]  F. Brezzi,et al.  A discourse on the stability conditions for mixed finite element formulations , 1990 .

[12]  Karan S. Surana,et al.  A SPACE-TIME COUPLED P-VERSION LEAST SQUARES FINITE ELEMENT FORMULATION FOR UNSTEADY TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS , 1996 .

[13]  Jan S. Hesthaven,et al.  Spectral penalty methods , 2000 .

[14]  J. P. Pontaza,et al.  Space-time coupled spectral/ hp least-squares finite element formulation for the incompressible Navier-Stokes equations , 2004 .

[15]  W. Couzy,et al.  A fast Schur complement method for the spectral element discretization of the incompressible Navier-Stokes equations , 1995 .

[16]  Gene H. Golub,et al.  Matrix computations , 1983 .

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

[18]  G. D. Davis.  Natural convection of air in a square cavity: A bench mark numerical solution , 1983 .

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

[20]  M. Gunzburger,et al.  Least-squares finite-element methods for optimization and control problems for the stokes equations , 2003 .

[21]  J. Reddy.  An introduction to the finite element method , 1989 .

[22]  D. Malkus,et al.  Mixed finite element methods—reduced and selective integration techniques: a unification of concepts , 1990 .

[23]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

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

[25]  D. Tritton.  Experiments on the flow past a circular cylinder at low Reynolds numbers , 1959, Journal of Fluid Mechanics.

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

[27]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

[28]  J. N. Reddy,et al.  Spectral/ hp least-squares finite element formulation for the Navier-Stokes equations , 2003 .