Text this: Comparison of Different Numerical Models for Solidification of Paraffin-Based PCMs and Evaluation of the Effect of Variable Mushy Zone Constant