Subgrid scale modeling of droplet bag breakup in VOF simulations

The mesh-dependency of the breakup of liquid films, including their breakup length scales and resulting drop size distributions, has long been an obstacle inhibiting the computational modeling of large-scale spray systems. With the aim of overcoming this barrier, this work presents a framework for the prediction and modeling of subgrid-thickness liquid film formation and breakup within two-phase simulations using the volume of fluid method. A two-plane interface reconstruction is used to capture the development of liquid films as their thickness decreases below the mesh size. The breakup of the film is predicted with a semi-analytical model that incorporates the film geometry captured through the two-plane reconstruction. The framework is validated against experiments of the bag breakup of a liquid drop at $\mathrm{We}=13.8$ through the comparison of the resulting drop size and velocity distributions. The generated distributions show good agreement with experimental results for drop resolutions as low as 25.6 cells across the initial diameter. The presented framework enables these drop breakup simulations to be performed at a computational cost three orders of magnitude lower than the cost of simulations utilizing adaptive mesh refinement.