Rayleigh–Bénard instability in the presence of phase boundary and shear

• Published in: Journal of fluid mechanics Vol. 948
• Main Authors: Lu, Cailei, Zhang, Mengqi, Luo, Kang, Wu, Jian, Yi, Hongliang
• Format: Journal Article
• Language: English
• Published: Cambridge, UK Cambridge University Press 10.10.2022
• Subjects: Couette flow; Energy; Energy absorption; Flow stability; Fluid flow; Heat; Heat flux; Heat transfer; Ice shelves; Instability; Investigations; JFM Papers; Laminar flow; Liquid phases; Melt temperature; Melting; Phase boundaries; Rayleigh number; Reynolds number; Shear flow; Stability analysis; Temperature; Topography; absolute/convective instability; Bénard convection; solidification/melting
• ISSN: 0022-1120; 1469-7645
• DOI: 10.1017/jfm.2022.723

We study the two-dimensional Rayleigh–Bénard instability subject to the combined effects of a solid–liquid phase boundary and shear using linear stability theory and energy analysis. We consider two thermal states of the solid (isothermal and conducting), and two types of shear that can arise in different contexts. When the melting temperature is equal to that held at the top boundary, three instability modes can arise with the increase of Reynolds number Re that characterizes the shear intensity: a boundary mode, a mixed boundary–bulk mode and a bulk flow mode. When the melting temperature lies between the top and the bottom boundaries, the introduction of Couette flow, independent of its intensity, always leads to the mixed mode, whereas the instability with Poiseuille flow is dominated by the bulk flow mode once Re exceeds a critical value, below which the mixed mode dominates. The energy analysis suggests that there exist two mechanisms by which the shear flow affects the system: one is by inhibiting the upward heat flux and another is by absorbing energy from the perturbed hydrodynamic field. These two mechanisms can play totally different roles in different cases. Results in the high-Re regime indicate that, when Re exceeds its classical threshold, i.e. Re = 5772.2 for Poiseuille flow, the Tollmien–Schlichting instability will be dominant in the present system.