  Smoothed particle hydrodynamics (SPH) is developed for modelling of melting
and solidification. Enthalpy method is used to solve heat conduction equations
which involved moving interface between phases. At first, we study the melting
of floating ice in the water for two-dimensional system. The ice objects are
assumed as solid particles floating in fluid particles. The fluid and solid
motion are governed by Navier-Stokes equation and basic rigid dynamics
equation, respectively. We also propose a strategy to separate solid particles
due to melting and solidification. Numerical results are obtained and plotted
for several initial conditions.
