  A novel simulation framework has been developed in this study for the direct
numerical simulation of the aerodynamic breakup of a vaporizing drop. The
interfacial multiphase flow with phase change is resolved using a consistent
geometric volume-of-fluid method. The bulk fluids are viscous and
incompressible with surface tension at the interface. The newly-developed
numerical methods have been implemented in the Basilisk solver, in which the
adaptive octree/quadtree mesh is used for spatial discretization, allowing
flexibility in dynamically refining the mesh in a user-defined region. The
simulation framework is extensively validated by a series of benchmark cases,
including the 1D Stefan and sucking problems, the growth of a 3D spherical
bubble in a superheated liquid, and a 2D film boiling problem. The simulation
results agree very well with the exact solution and previous numerical studies.
2D axisymmetric simulations were performed to resolve the vaporization of a
moving drop with a low Weber number in a high-temperature free stream. The
computed rate of volume loss agrees well with the empirical model of drop
evaporation. Finally, the validated solver is used to simulate the aerodynamic
breakup of an acetone drop at a high Weber number. A fully 3D simulation is
performed and the morphological evolution of the drop is accurately resolved.
The rate of vaporization is found to be significantly enhanced due to the drop
deformation and breakup. The drop volume decreases nonlinearly in time and at a
much higher rate than the empirical correlation for a spherical drop.
