September 15, 2026 | New peer-reviewed paper published
In a new article published (open access) in the Journal of Non-Newtonian Fluid Mechanics, we present a fully coupled implicit finite-volume method for the simulation of viscoelastic interfacial flows. Such flows are particularly challenging to compute at high elasticity, where numerical instabilities often limit conventional solution approaches. The new method solves the continuity, momentum, and polymer extra-stress equations together in a single coupled system and combines this formulation with a front-tracking method for the fluid interface. The approach retains the second-order accuracy of the underlying finite-volume discretisation and enables robust simulations over a wide range of elastic conditions. For a droplet in a shear-thinning viscoelastic liquid, simulations could be performed with a Weissenberg number of 10,000 without any numerical tricks, while simulations of a rising bubble reproduce the experimentally observed jump in terminal rise velocity and the associated negative wake. The results demonstrate the potential of fully coupled solution strategies for robust numerical simulations of strongly elastic interfacial flows.