Numerical study on melting of phase change material in an enclosure subject to Neumann boundary condition in the presence of Rayleigh-Bénard convection

• Published in: International journal of heat and mass transfer Vol. 171; p. 121103
• Main Authors: Parsazadeh, Mohammad, Malik, Mehtab, Duan, Xili, McDonald, André
• Format: Journal Article
• Language: English
• Published: Oxford Elsevier Ltd 01.06.2021; Elsevier BV
• Subjects: Boundary conditions; Domains; Electronic devices; Enclosures; Enthalpy; Fluid flow; Heat exchangers; Heat flux; Heat transfer; Liquid-solid interfaces; Melting; Neumann boundary condition; Nusselt number; Phase change heat transfer; Phase change material (PCM); Phase change materials; Rayleigh number; Rayleigh-Benard convection; Rayleigh-Bénard convection; Phase change material (PCM); Neumann boundary condition; Rayleigh-Bénard convection; Melting; Phase change heat transfer
• ISSN: 0017-9310; 1879-2189
• DOI: 10.1016/j.ijheatmasstransfer.2021.121103

•The effect of domain size on the phase change heat transfer process was identified.•The formation of different heat transfer subregimes was evaluated during melting.•A new correlation was developed for Nu as a function of Ste, Fo and Ras.•A new correlation was developed for Save as a function of Ste, Fo, Ras, and H. This research study investigated the melting process of a phase change material (PCM), heated from the bottom, under Neumann boundary conditions, in the presence of Rayleigh-Bénard convection. The problem was numerically simulated for a range of domain sizes (1×1−8×8cm2) and heat fluxes (0.5−2W/cm2) using the enthalpy porosity technique, which is mostly applicable in cooling heat exchangers of electronic devices. Scaling analysis and the numerical results were employed to find the relationship between the Nusselt number and the solid-liquid interface location with other dimensionless parameters and develop correlations to predict the Nusselt number and the solid-liquid interface location for this type of melting problem. The results of this research could be used to better understand the heat transfer regimes formed during melting in an enclosure heated from the bottom and predict the Nusselt number and the solid-liquid interface location. It was found that the temperature and Nusselt number fluctuations are initiated in the coarsening subregime when the Rayleigh number was not larger than 106 and thus may not occur for small-sized domains. It can be concluded the coarsening and turbulent subregimes may not occur in small enclosures, while the melting process mostly occurs during turbulent subregiem in large enclosures. The numerical results revealed the fastest advance of heat transfer rate and consequently, solid-liquid interface occur during the formation of Bénard cells when the Rayleigh number was on the order of 104. The Nusselt number and solid-liquid interface location correlations developed in this study were later validated using two new cases of numerical data. The results revealed that these correlations could accurately predict the Nusselt number and solid-liquid interface location.