Modeling and simulation of viscoelastic film retraction

In this paper, we investigate the retraction of a circular viscoelastic liquid film with a hole initially present in its center by means of finite element numerical simulations. We study the whole retraction process, aiming at understanding the hole opening dynamics both when the hole does not feel any confinement and when it interacts with the solid wall bounding the film. The retraction behavior is also interpreted through a simple toy model, that highlights the physical mechanism underlying the process.We consider three different viscoelastic constitutive equations, namely, Oldroyd-B, Giesekus (Gsk), and Phan Thien-Tanner (PTT) models, and several system geometries, in terms of the film initial radius and thickness. For each given geometry, we investigate the effects of liquid inertia, elasticity, and flow-dependent viscosity on the dynamics of the hole opening. Depending on the relative strength of such parameters, qualitatively different features can appear in the retracting film shape and dynamics.When inertia is relevant, as far as the opening hole does not interact withthe wall bounding the film, the influence of liquid elasticity is very moderate,and the retraction dynamics tends to the one of Newtonian sheets; whenthe hole starts to interact with the solid wall, hole radius/opening velocityoscillations are detected. Such oscillations enhance at increasing elasticity.From the morphological point of view, the formation of a rim at the edge ofthe retracting film is observed. If inertial forces become less relevant withrespect to viscous forces, R-oscillations disappear, the hole opening velocitygoes through a maximum and then monotonically decays to zero, and norim forms during the film retraction. Geometrical changes have the effect ofenlarging or reducing the portion of the retraction dynamics not influencedby the presence of the solid wall with respect to the one governed by thehole-wall interactions.