A numerical method for computing the motion of bubbles undergoing liquid-va
por 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 d
uring the whole calculation procedure. The modified level set method is app
lied for numerical simulation of bubble rise and growth in a stationary liq
uid. The numerical results are found to compare well with the data reported
in the literature and the analytical solutions.