Melting and vaporization of solids occur very often in nature and in man-made
processes. Many analytical and numerical solutions exist for solving the temperature
field in the liquid and solid regions, but inaccuracies persist in tracking the
phase change interfaces and the numerical solution of the temperature field is
usually cumbersome. The Boundary Element Method is proposed as an accurate,
efficient way to solve for the temperature field and the interface positions in a
phase change problem involving combined melting and vaporization. When
comparing to specific one-dimensional test cases, accurate results arc obtained
when using a sufficiently small time step. A comparison is made to existing data
from a laser drilling experiment. The anticipated physical effects which occur on
semi-infinite and finite domains arc confirmed. Consequently, this method can
be used to model natural and industrial phenomena involving phase change.