Premium
On the theory of semi‐implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 2: Implementation
Author(s) -
Gresho Philip M.,
Chan Stevens T.
Publication year - 1990
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1650110510
Subject(s) - mass matrix , finite element method , compressibility , flow (mathematics) , projection (relational algebra) , projection method , conservation of mass , matrix (chemical analysis) , incompressible flow , mathematics , pressure correction method , coupling (piping) , mechanics , mathematical optimization , computer science , geometry , algorithm , physics , mechanical engineering , engineering , structural engineering , materials science , dykstra's projection algorithm , composite material , nuclear physics , neutrino
Ever since the expansion of the finite element method (FEM) into unsteady fluid mechanics, the ‘consistent mass matrix’ has been a relevant issue. Applied to the time‐dependent incompressible Navier–Stokes equations, it virtually demands the use of implicit time integration methods in which full ‘velocity–pressure coupling’ is also inherent. The high cost of such (high‐quality) FEM calculations led to the development of simpler but ad hoc methods in which the ‘lumped’ mass matrix is employed and the velocity and pressure are uncoupled to the maximum extent possible. Resulting computer codes were less expensive to use but suffered a significant loss of accuracy, caused by lumping the mass when the flow was advection‐dominated and accurate transport of ‘information’ was important. In the second part of this paper we re‐introduce the consistent mass matrix into some semi‐implicit projection methods in such a way that the cost advantage of lumped mass and the accuracy advantage of consistent mass are simultaneously realized.