Abstract:Several finite element methods for simulating incompressible flows rely on the streamline upwind Petrov-Galerkin stabilization (SUPG) term, which is weighted by tau_SUPG. The conventional formulation of tau_SUPG includes a constant that depends on the time step size, producing an overall method that becomes exceedingly less accurate as the time step size approaches zero. In practice, such method inconsistency introduces significant error in the solution, especially in cardiovascular simulations, where small time step sizes may be required. To overcome this issue, we propose a consistent method that is based on a new definition of tau_SUPG. This method, which can be easily implemented on top of an existing streamline upwind Petrov-Galerkin and pressure stabilizing Petrov-Galerkin method, involves the replacement of the time step size in tau_SUPG with a physical time scale. This time scale is calculated in a simple operation once every time step for the entire computational domain from the ratio of the L2-norm of the acceleration and the velocity. The proposed method is compared against the conventional method using four cases: a steady pipe flow, a blood flow through vascular anatomy, an external flow over a square obstacle, and a fluid-structure interaction case. These numerical experiments show that the proposed formulation eliminates the inconsistency issue associated with the conventional formulation in all cases. While the proposed method is slightly more costly than the conventional method, it significantly reduces the error, particularly at small time step sizes. For the pipe flow where an exact solution is available, we show the conventional method can over-predict the pressure drop by a factor of three. This large error is almost completely eliminated by the proposed formulation, dropping to approximately 1% for all time step sizes and Reynolds numbers considered.

Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations

Development of a Finite Element/discontinuous Galerkin/level Set Approach for the Simulation of Incompressible Two Phase Flow.

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier-Stokes: Applications to fluid-structure interaction

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation

A Less Time-Consuming Upwind Compact Difference Method with Adjusted Dissipation Property for Solving the Unsteady Incompressible Navier-Stokes Equations

An Improved Finite Element Technique for the Simulation of Unsteady Incompressible Flows

A New Combined Finite Element-Upwind Finite Volume Method for Convection-Dominated Diffusion Problems

Finite Difference Streamline Diffusion Method Using Nonconforming Space for Incompressible Time-Dependent Navier-Stokes Equations

Finite element methods for multicomponent convection-diffusion

A Higher Order Compact Finite Difference Algorithm for Solving the Incompressible Navier–Stokes Equations

A time-consistent stabilized finite element method for fluids with applications to hemodynamics

A Parallel High-Order Compact Scheme for the Pure Streamfunction Formulation of the 3D Unsteady Incompressible Navier–Stokes Equation

Least-squares Finite Element Method for Unsteady Stress Formulation of Navier-Stokes Equations

High-order finite element methods for three-dimensional multicomponent convection-diffusion

A fully conservative and shift-invariant formulation for Galerkin discretizations of incompressible variable density flow

A high‐order compact scheme for solving the 2D steady incompressible Navier‐Stokes equations in general curvilinear coordinates

An advanced algorithm for solving incompressible fluid dynamics: from Navier–Stokes to Poisson equations

A super-convergent quasi-discontinuous Galerkin method for approximating inertia-gravity and Rossby waves in geophysical flows

A Meshless Local Petrov–Galerkin Approach for Solving the Convection-Dominated Problems Based on the Streamline Upwind Idea and the Variational Multiscale Concept

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

A gas-kinetic theory based multidimensional high-order method for the compressible Navier–Stokes solutions