Numerical analysis of vertical solidification process allowing solid-liquid density change is performed by a hybrid method between a winite volume method (FVM) and a finite element method (FEM). The investigation focuses on the influence of solid-liquid density change and cooling rates on the motion of solid-liquid interface, solidified mass fraction, temperatures and thermal stresses in the solid region. Due to the density change of pure aluminium, solid-liquid interface moves more slowly but the solidified mass fraction is larger. The cooling rate of the wall is shown to have a significant influence on the phase change heat transfer and thermal stresses, while the density change has a small influence on the motion of the interface, solidified mass fraction, temperature distributions and thermal stresses. As the cooling rate increases, the thermal stresses become higher at the early stage of a solidification process, but it has small influence on the final stresses as the steady state is reached.