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

Authors:
V. Prabhakar;J. N. Reddy
Affiliations:
Department of Mechanical Engineering, Texas A&M University, College Station, TX 77843-3123, USA;Department of Mechanical Engineering, Texas A&M University, College Station, TX 77843-3123, USA
Venue:
Journal of Computational Physics
Year:
2006

Citing 8
Cited 2

A discourse on the stability conditions for mixed finite element formulations

Computer Methods in Applied Mechanics and Engineering
Mixed finite element methods—reduced and selective integration techniques: a unification of concepts

Computer Methods in Applied Mechanics and Engineering - Special edition on the 20th Anniversary
A fast Schur complement method for the spectral element discretization of the incompressible Navier-Stokes equations

Journal of Computational Physics
Least-Squares Finite Element Method for the Stokes Problem with Zero Residual of Mass Conservation

SIAM Journal on Numerical Analysis
Spectral penalty methods

Proceedings of the fourth international conference on Spectral and high order methods (ICOSAHOM 1998)
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems
A spectral multidomain penalty method model for the simulation of high Reynolds number localized incompressible stratified turbulence

Journal of Computational Physics
Least-squares finite-element methods for optimization and control problems for the stokes equations

Computers & Mathematics with Applications

On the roles of minimization and linearization in least-squares finite element models of nonlinear boundary-value problems

Journal of Computational Physics
Least-squares spectral element solution of incompressible Navier-Stokes equations with adaptive refinement

Journal of Computational Physics

Quantified Score

Hi-index	31.46

Visualization

Abstract

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^-^1-@c(@?.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.