High-order finite volume methods for viscoelastic flow problems

Abstract

In this paper accurate and stable finite volume schemes for solving viscoelastic flow problems are presented. Two contrasting finite volume schemes are described: a hybrid cell-vertex scheme and a pure cell-centred counterpart. Both schemes employ a time-splitting algorithm to evolve the solution through time towards steady state. In the case of the hybrid scheme, a semi-implicit formulation is employed in the momentum equation, based on the Taylor–Galerkin approach with a pressure-correction step to enforce incompressibility. The basis of the pure finite volume approach is a backward Euler scheme with a semi-Lagrangian step to treat the convection terms in the momentum and constitutive equations. Two distinct finite volume schemes are presented for solving the systems of partial differential equations describing the flow of viscoelastic fluids. The schemes are constructed to be second-order accurate in space. The issue of stability is also addressed with respect to the treatment of convection. Numerical examples are presented illustrating the performance of these schemes on some steady and transient problems that possess analytical solutions.