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

Milad Tajik Jamalabad (), Cristobal Cortes, Javier Pallarés, Antonia Gil and Inmaculada Arauzo
Additional contact information
Milad Tajik Jamalabad: Institute of Energy and Resources Efficiency of Aragón (ENERGAIA), University of Zaragoza, 50018 Zaragoza, Spain
Cristobal Cortes: Institute of Energy and Resources Efficiency of Aragón (ENERGAIA), University of Zaragoza, 50018 Zaragoza, Spain
Javier Pallarés: Institute of Energy and Resources Efficiency of Aragón (ENERGAIA), University of Zaragoza, 50018 Zaragoza, Spain
Antonia Gil: Institute of Energy and Resources Efficiency of Aragón (ENERGAIA), University of Zaragoza, 50018 Zaragoza, Spain
Inmaculada Arauzo: Institute of Energy and Resources Efficiency of Aragón (ENERGAIA), University of Zaragoza, 50018 Zaragoza, Spain

Energies, 2024, vol. 17, issue 23, 1-17

Abstract: The impact of the mushy zone parameter ( A m u s h y ) and the chosen numerical model during the solidification of a commercial paraffin-type phase change material (PCM) in a vertical cylinder under T-history conditions was examined through numerical simulations. The cooling process was modeled using three methods implemented in the CFD software ANSYS Fluent 2020 R2: the enthalpy–porosity method, the apparent heat capacity (AHC) method, and a new model proposed by the authors which incorporates heat capacity directly into ANSYS Fluent. To accurately define the boundary conditions, radiative heat transfer between surfaces was taken into account. Furthermore, the influence of the mushy zone parameter on the simulation accuracy and solidification rate was investigated, with the parameter being treated as a function of the liquid fraction. The results indicate that the proposed model aligns closely with experimental data regarding cooling temperature, offering better predictions compared to the other models. It was observed that temperature varies with time but not with position, suggesting that this model more effectively satisfies the lumped system condition—an essential characteristic of the T-history experiment—compared to the other methods. Additionally, the analysis showed that a higher mushy zone parameter enhances the accuracy of simulations and predicts a shorter solidification time; approximately 11% for the E-p and 7% for the AHC model. Using a variable mushy zone parameter based on the liquid fraction also produced similar results, resulting in an increased solidification rate.

Keywords: phase change material; variable mushy zone parameter; solidification; numerical simulation; paraffin-type PCM (search for similar items in EconPapers)
JEL-codes: Q Q0 Q4 Q40 Q41 Q42 Q43 Q47 Q48 Q49 (search for similar items in EconPapers)
Date: 2024
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/1996-1073/17/23/5982/pdf (application/pdf)
https://www.mdpi.com/1996-1073/17/23/5982/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jeners:v:17:y:2024:i:23:p:5982-:d:1531606

Access Statistics for this article

Energies is currently edited by Ms. Agatha Cao

More articles in Energies from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().