Abstract
A numerical method for computing the motion of bubbles undergoing liquid?vapor phase change is presented. The method is based on a level set technique for capturing the phase interface, which is modified to include the effect of phase change at the interface as well as to achieve mass conservation during the whole calculation procedure. The modified level set method is applied for numerical simulation of bubble rise and growth in a stationary liquid. The numerical results are found to compare well with the data reported in the literature and the analytical solutions.