Asymmetric equilibrium states for melting and freezing in thermal convection

A block of ice in a box heated from below and cooled from above can (partially) melt. Vice versa, a box of water with less heating from below or more cooling from above can (partially) re-solidify. This study investigates the asymmetric behaviours between such melting and freezing processes in this Rayleigh–Bénard geometry, focusing on differences in equilibrium flow structures, solid–liquid interface morphology, and equilibrium mean interface height. Our findings reveal a robust asymmetry across a range of Rayleigh numbers and top cooling temperature (i.e. hysteretic behaviour), where the evolution of freezing shows a unique ‘splitting event’ of convection cells that leads to a non-monotonic height evolution trend. To characterise the differences between melting and freezing, we introduce an effective Rayleigh number and the aspect ratio for the cellular structures, and apply the heat flux balance and the Grossmann–Lohse theory. Based on this, we develop a unifying model for the melting and freezing behaviour across various conditions, accurately predicting equilibrium states for both phase-change processes. This work provides insights into the role of convective dynamics in phase-change symmetry-breaking, offering a framework applicable to diverse systems involving melting and freezing.