Numerical solution of an extended White–Metzner model for eccentric Taylor–Couette flow

In this study, we have developed a new numerical approach to solve differential-type viscoelastic fluid models for a commonly used benchmark problem, namely, the steady Taylor—Couette flow between eccentric cylinders. The proposed numerical approach is special in that the nonlinear system of discretized algebraic flow equations is solved iteratively using a Newton–Krylov method along with an inverse-based incomplete lower-upper preconditioner. The numerical approach has been validated by solving the benchmark problem for the upper-convected Maxwell model at a large Deborah number. Excellent agreement with the numerical data reported in the literature has been found. In addition, a parameter study was performed for an extended White–Metzner model. A large eccentricity ratio was chosen for the cylinder system in order to allow flow recirculation to occur. We detected several interesting phenomena caused by the large eccentricity ratio of the cylinder system and by the viscoelastic nature of the fluid. Encouraged by the results of this study, we intend to investigate other polymeric fluids having a more complex microstructure in an eccentric annular flow field.

► New preconditioned Newton–Krylov method for solving highly nonlinear systems of discretized algebraic flow equations.
► Application of inverse-based ILU preconditioner makes iterative solution approach robust.
► Results provide new insight into rheology-structure relationships of polymeric liquids.