Diffusion of heat during phase transition When one substance and one phase is present, the diffusion of temperature is given by the heat equation
$$ \partial_t T = \nabla\cdot(\kappa\nabla T) $$
Can the heat equation be modified to accommodate phase transitions? If so, how?
 A: If we look at the transition from liquid to solid you can think of the problem as a Stefan Problem or moving boundary problem. Assuming a constant diffusion coefficient in the different phases and working with some dimensionless temperature we get the 2 heat equations in the bulk regions
$$ D_i \nabla^{2} T = \partial_{t} T $$ 
Where $ i = s,l $ for solid or liquid. In order to describe heat conservation due to the phase transition at the boundary you get the interface boundary condition $$  v_n = [D_i \partial_n T ]^{s}_{l} $$ 
Where $ v_n $ is the normal velocity of the interface and $ \partial_n $ is the normal derivative. The brackets means the difference of the gradients evaluated on the both sides of the interface. Finally to include surface tension and kinetic cooling the boundary condition for the temperature becomes
$$ T|_{\Gamma} = - d_0 \kappa - \beta v_n $$
Where $ d_0 $ is the capillary length, $ \kappa $ is the curvature and $ \beta $ is some kinetic coefficient. For references see for example James Langers Instabilities and pattern formation in crystal growth. For a discussion on how to numerically solve the equations using phase fields you could take a look at the relevant chapters in Phase-Field Methods in Materials Science and Engineering by Elder and Provatas.  
