Abstract:The elastic viscous split stress formulation (EVSS) and the discrete EVSS formulation (DEVSS) are effective in stabilizing numerical simulations of viscoelastic flows and have been widely used. Following the concept of Galerkin least-square perturbations proposed by Hughes et al. [Comput. Meth. Appl. Mech. Eng. 73 (1989) 173–189] and Franca et al. [SIAM J. Numer. Anal. 28(6) (1991) 1680–1697; Comput . Meth. Appl. Mech. Eng. 99 (1992) 209–233; Ibid. 104 (1993) 31–48] we are able to give the DEVSS formulation a new explanation as a perturbation to the Galerkin method based on the strain-rate residual, and furthermore, introduce another stabilized formulation, here named as MIX1, based on the incompressibility residual of the finite element discretizations. The three formulations (EVSS, DEVSS, MIX1), combined with a h–p type finite element algorithm that employs the SUPG technique to solve the viscoelastic constitutive equations are then tested on three benchmark problems: the flow of the upper-convected Maxwell fluid between eccentric cylinders, the flow of the Maxwell fluid around a sphere in a tube and the flow of the Maxwell and Oldroyd-B fluids around a cylinder in a channel. The results are checked with previous published works; good agreement is observed. Our numerical experiments convincingly demonstrate that the MIX1 is an accurate algorithm and convergent in terms of the p-extension, it has the same level of stability and robustness as the DEVSS method and is superior to the EVSS method in some respects. More important is that with MIX1 method one needs not solve for the strain-rate tensor as in EVSS and DEVSS methods, therefore, the CPU time consumption in the MIX1 method especially when using a coupled iteration scheme can be radically reduced. The success of the MIX1 method presents a challenge to the widely accepted concept of making the momentum equation explicitly elliptic.

Fully Discrete Finite Element Method Based on Second-Order Crank-Nicolson/Adams-Bashforth Scheme for the Equations of Motion of Oldroyd Fluids of Order One

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

A discontinuous Galerkin finite element method for the Oldroyd model of order one

Convergence analysis and error estimates for a second order accurate finite element method for the Cahn–Hilliard–Navier–Stokes system

Second-order accurate numerical scheme with graded meshes for the nonlinear partial integrodifferential equation arising from viscoelasticity

Novel numerical analysis of multi-term time fractional viscoelastic non-Newtonian fluid models for simulating unsteady MHD Couette flow of a generalized Oldroyd-B fluid

Unconditionally stable fully‐discrete finite element numerical scheme for active fluid model

Discretisation of an Oldroyd-B viscoelastic fluid flow using a Lie derivative formulation

Positive-definiteness preserving and energy stable time-marching scheme for a diffusive Oldroyd-B electrohydrodynamic model

A Two-level Finite Element Method for Viscoelastic Fluid Flow: Non-smooth Initial Data

Numerical Simulation of Planar 4: 1 Contraction Flow of a Viscoelastic Fluid Using a Higher Upwind Finite Volume Method

Second-order energy-stable scheme and superconvergence for the finite difference method on non-uniform grids for the viscous Cahn–Hilliard equation

Linear and Energy-Stable Method with Enhanced Consistency for the Incompressible Cahn–Hilliard–Navier–Stokes Two-Phase Flow Model

Semi-Lagrangian finite element exterior calculus for incompressible flows

A two-stage fourth order time-accurate discretization for Lax--Wendroff type flow solvers II. High order numerical boundary conditions

Galerkin/least-square Finite-Element Methods for Steady Viscoelastic Flows

High Order Compact Generalized Finite Difference Methods for Solving Inviscid Compressible Flows.

A structure-preserving projection method with formal second-order accuracy for the incompressible Navier–Stokes equations

Compact High Order Finite Volume Method on Unstructured Grids II: Extension to Two-Dimensional Euler Equations.

Exponential Compact ADI Method for a Coupled System of Convection-Diffusion Equations Arising from the 2D Unsteady Magnetohydrodynamic (MHD) Flows

Numerical approximations of the Navier–Stokes equation coupled with volume-conserved multi-phase-field vesicles system: Fully-decoupled, linear, unconditionally energy stable and second-order time-accurate numerical scheme