Abstract

The log-conformation reformulation, originally proposed by Fattal and Kupferman (2004), allows computing incompressible viscoelastic problems with high Weissenberg numbers which are impossible to solve with the typical three-field formulation. By following this approach, in this work we develop a new stabilized finite element formulation based on the logarithmic reformulation using the Variational Multiscale (VMS) method as stabilization technique, together with a modified log-conformation formulation. Our approach follows the term-by-term stabilization proposed by Castillo and Codina (2014) for the standard formulation, which is more effective when there are stress singularities. The formulation can be used when the relaxation parameter is set to zero, and permits a direct steady numerical resolution. The formulation is validated in the classical benchmark flow past a cylinder and in the well-known planar contraction 4:1, achieving very accurate, stable and mesh independent results for highly elastic fluids.