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

In this article, we present a numerical method for the Stokes flow of an Oldroyd-B fluid. The viscoelastic stress evolves according to a constitutive law formulated in terms of the upper convected time derivative. A finite difference method is used to discretise along fluid trajectories to approximate the advection and deformation terms of the upper convected derivative in a simple, cheap and cohesive manner, as well as ensuring that the discrete conformation tensor is positive definite. A full implementation with coupling to the fluid flow is presented, along with a detailed discussion of the issues that arise with such schemes. We demonstrate the performance of this method with detailed numerical experiments in a lid-driven cavity setup. Numerical results are benchmarked against published data, and the method is shown to perform well in this challenging case.