Search

The dynamics of heat transfer in a model system of Rayleigh-B\'enard (RB) convection reduced to its essential, here dubbed Burgers-Rayleigh-B\'enard (BRB), is studied. The system is spatially one-dimensional, the flow field is compressible and its evolution is described by the Burgers equation forced by an active temperature field. The BRB dynamics shares some remarkable similarities with realistic RB thermal convection in higher spatial dimensions: i) it has a supercritical pitchfork instability for the onset of convection which solely depends on the Rayleigh number $(Ra)$ and not on Prandlt $(Pr)$, occurring at the critical value $Ra_c = (2\pi)^4$ ii) the convective regime is spatially organized in distinct boundary-layers and bulk regions, iii) the asymptotic high $Ra$ limit displays the Nusselt and Reynolds numbers scaling regime $Nu = \sqrt{RaPr}/4$ for $Pr\ll 1$, $Nu=\sqrt{Ra}/(4\sqrt{\pi})$ for $Pr\gg1$ and $Re = \sqrt{Ra/Pr}/\sqrt{12}$, thus making BRB the simplest wall-bounded convective system exhibiting the so called ultimate regime of convection. These scaling laws, derived analytically through a matched asymptotic analysis are fully supported by the results of the accompanying numerical simulations. A major difference with realistic natural convection is the absence of turbulence. The BRB dynamics is stationary at any $Ra$ number above the onset of convection. This feature results from a nonlinear saturation mechanism whose existence is grasped by means of a two-mode truncated equation system and via a stability analysis of the convective regime., Comment: 17 pages, 9 figures

Melting and, conversely, solidification processes in the presence of convection are key to many geophysical problems. An essential question related to these phenomena concerns the estimation of the (time-evolving) melting rate, which is tightly connected to the turbulent convective dynamics in the bulk of the melt fluid and the heat transfer at the liquid-solid interface. In this work, we consider a convective-melting model, constructed as a generalization of the Rayleigh-B\'enard system, accounting for the basal melting of a solid. As the change of phase proceeds, a fluid layer grows at the heated bottom of the system and eventually reaches a turbulent convection state. By means of extensive Lattice-Boltzmann numerical simulations employing an enthalpy formulation of the governing equations, we explore the model dynamics in two and three-dimensional configurations. The focus of the analysis is on the scaling of global quantities like the heat flux and the kinetic energy with the Rayleigh number, as well as on the interface morphology and the effects of space dimensionality. Independently of dimensionality, we find that the convective-melting system behavior shares strong resemblances with that of the Rayleigh-B\'enard one, and that the heat flux is only weakly enhanced with respect to that case. Such similarities are understood, at least to some extent, considering the resulting slow motion of the melting front (with respect to the turbulent fluid velocity fluctuations) and its generally little roughness (compared to the height of the fluid layer). Varying the Stefan number, accounting for the thermodynamical properties of the material, also seems to have only a mild effect, which implies the possibility to extrapolate results in numerically delicate low-Stefan setups from more convenient high-Stefan ones., Comment: 22 pages, 10 figures (flow visualizations in low resolution), videos are not included