Numerical modelling of phase change in freezing and thawing unsaturated soil