The FVM based on the artificial compressibility algorithm [27,28], which has been successfully applied to the problem of magnetohydrodynamic (MHD) convection [29,30] and nanofluid flows [31,32], is employed to solve the coupled non-linear differential equations describing the free convection of NePCMs suspension. The computational domain is discretized using a uniform staggered grid system, wherein the scalar quantities and the velocity components are stored at the center of the control volume