Heat conduction is an unsteady process during which heater pads are turned on and off sequentially and heat flows through different solid materials. Heat conduction can be expressed as a partial differential equation with an added term to model phase changes:

where is the density of the solid, its specific heat, its conductivity and H is the enthalpy. The unknown is the variation of the nodal temperature between iterations (equals 0 at convergence).

At each time step, starting from an initial temperature and enthalpy , the heat conduction equation is solved iteratively until the new temperature

matches the temperature corresponding to the new enthalpy

A table that correlates the enthalpy with temperature T for each material must be provided by you.