Evolution of solid–liquid interface in bottom heated cavity for low Prandtl number using lattice Boltzmann method

Abstract

The influence of natural convection cells on heat transfer and the evolution of melt interface is studied for low Prandtl number fluid (Pr = 0.025) in phase-change Rayleigh–Benard convection using the lattice Boltzmann method. The thermal lattice Boltzmann model is used to evaluate the effect of Rayleigh number (Ra = 6708, 11 708, and 21 708) and cavity aspect ratio (γ = 0.062 5, 0.125, 0.25, 0.5, and 1) on the onset of convection, number of convection cells, and Nusselt number in the classical Rayleigh–Benard convection. The modified equilibrium distribution function-based thermal lattice Boltzmann model is applied to evaluate the effect of Stefan number (Ste = 0.025, 0.05, and 0.1) in the phase change Rayleigh–Benard convection. Distinct flow configurations depend on the Rayleigh number, aspect ratio, and Stefan number. The number of convection cells follows an inverse relation with the aspect ratio. Nusselt number increases with decreasing cavity aspect ratio and increasing Rayleigh number in the classical Rayleigh–Benard convection. With the variation in the aspect ratio based on the melt layer height during melting of phase change material, the number of convection cells changes resulting in the change in the evolution of the melt interface and convective heat transfer. Melting in a cavity of aspect ratio less than 0.5, the evolution of melt interface remains symmetrical. For an aspect ratio greater than 0.5, the interface evolution becomes unsymmetrical depending on the transition to single convection cell-dominated heat transfer.