Journal of Non-Newtonian Fluid Mechanics, Vol.126, No.2-3, 123-143, 2005
Computation of incompressible and weakly-compressible viscoelastic liquids flow: finite element/volume schemes
In this study, we analyse viscoelastic numerical solution for an Oldroyd-B model under incompressible and weakly-compressible liquid flow conditions. We consider flow through a planar four-to-one contraction, as a standard benchmark, throughout a range of Weissenberg numbers up to critical levels. At the same time, inertial and creeping flow settings are also addressed. Within our scheme, we compare and contrast, two forms of stress discretisation, both embedded within a high-order pressure-correction time-marching formulation based on triangles. This encompasses a parent-cell finite element/SUPG scheme, with quadratic stress interpolation and recovery of velocity gradients. The second scheme involves a sub-cell finite volume implementation, a hybrid fe/fv-scheme for the full system. A new feature of this study is that both numerical configurations are able to accommodate incompressible, and low to vanishing Mach number compressible liquid flows. This is of some interest within industrial application areas. We are able to provide parity between the numerical solutions across schemes for any given flow setting. Close examination of flow patterns and vortex trends indicates the broad differences anticipated between incompressible and weakly-compressible solutions. Vortex reduction with increasing Weissenberg number is a common feature throughout. Compressible solutions provide larger vortices (salient and lip) than their incompressible counterparts, and larger stress patterns in the re-entrant corner neighbourhood. Inertia tends to reduce such phenomena in all instances. The hybrid fe/fv-scheme proves more robust, in that it captures the stress singularity more tightly than the fe-form at comparable Weissenberg numbers, reaching higher critical levels. The sub-cell structure, the handling of cross-stream numerical diffusion, and corner discontinuity capturing features of the hybrid fe/fv-scheme, are all perceived as attractive additional benefits that give preference to this choice of scheme. (c) 2004 Elsevier B.V. All rights reserved.
Keywords:compressible liquid flow;high-order pressure-correction;low Mach number;Oldroyd-B;planar contraction;critical Weissenberg number;vortex activity