Graphite material in the form of expanded has proven to be an effective material to remedy the low thermal conductivity of paraffin storage medium. The present work focuses on the solidification performance of carbon-based approach via graphite matrix composite with phase change to eliminate the bottleneck of low thermal conductivity issue of phase change materials for thermal energy storage applications beyond the current literature. Paraffin/graphite matrix composites with various bulk density values of 0 kg/m3 (Case 1), 23 kg/m3 (Case 2), 50 kg/m3 (Case 3), 100 kg/m3 (Case 4), and 143 kg/m3 (Case 5) under various operating conditions of wall temperature (Twall: 15 °C, 25 °C, 35 °C, and 45 °C), which mimic real applications' fluid temperature, are performed numerically in detail. The numerical model was validated with experimental outputs performed by the authors. Numerical results show that paraffin/graphite matrix composites have a great potential for the solidification/discharging process to recover the heat stored in the charging/melting period. Uniform and rapid solidification behaviour is achieved with the paraffin/graphite matrix composite in terms of conduction-dominated heat transfer via abundant thermal bridges compared to the pure paraffin case. Higher bulk density values and lower cold surface temperatures present superior performance in solidification rate and thermal energy recovery rate. However, the effect of the bulk density is slowed down with the beyond of certain value of bulk density (100 kg/m3). For 100 kg/m3 (Case 4) at Twall = 15 °C, the heat recovery rate is enhanced by 31.17 times and 3.84 times compared to the 0 kg/m3 and 23 kg/m3. The heat recovery rate of 100 kg/m3 is improved by 350.8 % for Twall = 15 °C in comparison with Twall = 45 °C. The total solidification rate (t = 17.5 min) increased by 90.94 % and 62.94 % for 100 kg/m3 at Twall = 35 °C compared to 0 kg/m3 and 23 kg/m3, respectively.