A regularised one-dimensional drop formation and coalescence model using a total variation diminishing (TVD) scheme on a single Eulerian grid

The breakup of an axisymmetric viscous jet is considered in the lubrication approximation. The discretised equations are solved on a fixed equidistant one-dimensional Eulerian grid. The governing equations are implemented in a conservative second order accurate total variation diminishing (TVD) scheme, preventing the numerical diffusivity. Singularities that occur at pinchoff and coalescence are regularised by a small modification on the surface tension. The modification is of the order of the spatial step Delta x. This regularisation ensures that the solution of the presented numerical model converges to the exact solution of the breakup of a jet in the lubrication approximation. The results of the presented numerical model agree quantitatively with the analytical solution of the Rayleigh-Plateau instability, and with experimental results on the final stage of the Rayleigh-Plateau instability.