Solution behavior of the falling sphere problem in viscoelastic flows

Fan Yurun
DOI: https://doi.org/10.1007/BF02484574
IF: 3.5
2003-01-01
Acta Mechanica Sinica
Abstract:The hp -finite element method and pseudo arc-length method are combined to track solution behavior of the steady flow of a sphere falling in a tube filled with viscoelastic fluids. The computations is proved to be convergent and stable by a posteriori error analyses; the solutions smooth and with extremely steep stress gradient are obtained by using appropriate high-order interpolation distributions. The commonly used drag coefficient is proved not a reliable indicator for the approximation error. The solution curves of the upper convected Maxwell (UCM) and the Oldroyd-B fluids clearly indicate that there exist limiting points, parameterized by the Deborah number, where the magnitude of the solutions and the axial normal stress gradients behind the rear stagnation rise up rapidly, maybe go to infinity. The limiting points can be simply removed by employing the modified Chilcott and Rallison (MCR) fluid, which replaces the infinitely extendable, linear Gaussian chain with the finitely extendable nonlinear elastic spring (the FENE chain). Therefore the limiting behavior is caused by the physical model; it is not a numerical artifact.
What problem does this paper attempt to address?