Numerical simulation of viscous filament stretching flows

A numerical study on the stretching of a Newtonian fluid filament is analysed. Stretching is performed between two retracting plates, moving under constant extension rate. A semi‐implicit Taylor‐Galerkin/pressure‐correction finite element formulation is employed on variable‐structure triangular meshes. Stability and accuracy of the scheme is maintained up to large Hencky‐strain levels. A non‐uniform radius profile, minimum at the filament mid‐plane, is observed along the filament‐length at all times. We have found maintenance of a suitable mesh aspect‐ratio around the mid‐plane region (maximum stretch zone) to restrict early filament break‐up and consequently solution divergence. As such, true transient flow evolution is traced and the numerical results bear close agreement with the literature.