An efficient and accurate fully discrete finite element method for unsteady incompressible Oldroyd fluids with large time step

This paper proposes a second‐order accuracy in time fully discrete finite element method for the Oldroyd fluids of order one. This new approach is based on a finite element approximation for the space discretization, the Crank–Nicolson/Adams–Bashforth scheme for the time discretization and the trapezoid rule for the integral term discretization. It reduces the nonlinear equations to almost unconditionally stable and convergent systems of linear equations that can be solved efficiently and accurately. Here, the numerical simulations for L², H¹ error estimates of the velocity and L² error estimates of the pressure at different values of viscoelastic viscosities α, different values of relaxation time λ₁, different values of null viscosity coefficient μ₀ are shown. In addition, two benchmark problems of Oldroyd fluids with different solvent viscosity μ and different relaxation time λ₁ are simulated. All numerical results perfectly match with the theoretical analysis and show that the developed approach gives a high accuracy to simulate the Oldroyd fluids under a large time step. Furthermore, the difference and the connection between the Newton fluids and the viscoelastic Oldroyd fluids are displayed. Copyright © 2015 John Wiley & Sons, Ltd.