Second-order finite difference approximations of the upper-convected time derivative

TL;DR

This paper introduces second-order finite difference schemes for the upper-convected time derivative, improving accuracy in viscoelastic models and demonstrating their effectiveness through numerical experiments in one and two dimensions.

Contribution

The work develops new second-order finite difference schemes for the upper-convected derivative, with theoretical error analysis and practical application to viscoelastic constitutive equations.

Findings

Second-order accuracy achieved in one and two dimensions.

Numerical results confirm theoretical error estimates.

Effective application to Oldroyd-B model at various Weissenberg numbers.

Abstract

In this work, new finite difference schemes are presented for dealing with the upper-convected time derivative in the context of the generalized Lie derivative. The upper-convected time derivative, which is usually encountered in the constitutive equation of the popular viscoelastic models, is reformulated in order to obtain approximations of second-order in time for solving a simplified constitutive equation in one and two dimensions. The theoretical analysis of the truncation errors of the methods takes into account the linear and quadratic interpolation operators based on a Lagrangian framework. Numerical experiments illustrating the theoretical results for the model equation defined in one and two dimensions are included. Finally, the finite difference approximations of second-order in time are also applied for solving a two-dimensional Oldroyd-B constitutive equation subjected to a prescribed velocity field at different Weissenberg numbers.