The characteristic-based split (CBS) scheme for viscoelastic flow past a circular cylinder

Abstract

A fully explicit, characteristic-based split (CBS) method for viscoelastic flow past a circular cylinder, placed in a rectangular channel, is presented. The pressure equation in its explicit form is employed via an artificial compressibility parameter. The constitutive equations used here are based on the Oldroyd-B model. No loss of convergence to steady state was observed in any of the results presented in this paper. Comparison of the present results with other available numerical data shows that the CBS algorithm is in excellent agreement with them at lower Deborah numbers. However, at higher Deborah numbers, the present results differ from other numerical solutions. This is due to the fact that the positive definitiveness of the conformation matrix is lost between a Deborah number of 0.6 and 0.7. However, the positive definitiveness is retained when an artificial diffusion is added to the discrete constitutive equations at higher Deborah numbers. It appears that the fractional solution stages used in the CBS scheme and the higher-order time step-based convection stabilization clearly reduce the instability at higher Deborah numbers. The Deborah number limit reached in the present work is three without artificial dissipation and two with artificial dissipation. Copyright © 2007 John Wiley & Sons, Ltd.