FEM Techniques for the LCR Reformulation of Viscoelastic Flow Problems

We present special numerical techniques for viscoelastic fluid flow utilizing a fully coupled monolithic multigrid finite element approach with consistent edge-oriented stabilization technique. The governing equations arise from the Navier–Stokes for the Oldroyd-B type of fluid with the help of the log-conformation reformulation to allow a wide range of Weissenberg numbers. The resulting nonlinear system consists of 6 variables for velocity, pressure and the logarithm of the conformation stress tensor in 2D. The system is discretized in time by using a fully implicit second order accurate time integrator. In each time step, we have to solve a discretized system in space employing the high order finite element triple Q 2∕P 1 disc ∕Q 2. We utilize the discrete damped Newton method with divided differences for handling the Jacobian, and apply a geometrical multigrid solver with a special Vanka smoother to handle the linear subproblems. Local refinement can be assigned at regions of interest to reduce the computational cost. The presented methodology is implemented on the open source software package FeatFlow (www.featflow.de) and validated for several well-known benchmark problems.