We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
We numerically study the melting process of a solid layer heated from below such that a liquid melt layer develops underneath. The objective is to quantitatively describe and understand the emerging topography of the structures (characterized by the amplitude and wavelength), which evolve out of an initially smooth surface. By performing both two-dimensional (achieving Rayleigh number up to 
$Ra=10^{11}$) and three-dimensional (achieving Rayleigh number up to 
$Ra=10^9$) direct numerical simulations with an advanced finite difference solver coupled to the phase-field method, we show how the interface roughness is spontaneously generated by thermal convection. With increasing height of the melt the convective flow intensifies and eventually even becomes turbulent. As a consequence, the interface becomes rougher but is still coupled to the flow topology. The emerging structure of the interface coincides with the regions of rising hot plumes and descending cold plumes. We find that the roughness amplitude scales with the mean height of the liquid layer. We derive this scaling relation from the Stefan boundary condition and relate it to the non-uniform distribution of heat flux at the interface. For two-dimensional cases, we further quantify the horizontal length scale of the morphology, based on the theoretical upper and lower bounds given for the size of convective cells known from Wang et al. (Phys. Rev. Lett., vol. 125, 2020, 074501). These bounds agree with our simulation results. Our findings highlight the key connection between the morphology of the melting solid and the convective flow structures in the melt beneath.
We derive from first principles analytic relations for the second- and third-order moments of 
$\boldsymbol{\mathsf{m}}$, the spatial gradient of fluid velocity 
$\boldsymbol{u}$, 
$\boldsymbol{\mathsf{m}} = \nabla \boldsymbol{u}$, in compressible turbulence, which generalize known relations in incompressible flows. These relations, although derived for homogeneous flows, hold approximately for a mixing layer. We also discuss how to apply these relations to determine all the second- and third-order moments of the velocity gradient experimentally for isotropic compressible turbulence.
Enhanced fluctuations, steep gradients, and intensified heat transfer are characteristics of wall-bounded turbulence at transcritical conditions. Although such conditions are prevalent in numerous technical applications, the structure of the thermal boundary layer under realistic density gradients and heating conditions remains poorly understood. Specifically, statistical descriptions of the temperature field in such flows are provided inconsistently using existing models. To address this issue, direct numerical simulations are performed by considering fully developed transcritical turbulent channel flow at pressure and temperature conditions that cause density changes of a factor of up to 
$O(20)$ between the hot and cold walls. As a consequence of the proximity of the Widom line to the hot wall, significant asymmetries are observed when comparing regions near the cold wall and near the hot wall. Previous transformations that attempt to collapse the near-wall mean temperature profiles among different cases to a single curve are examined. By addressing model deficiencies of these transformations, a formulation for an improved mean temperature transformation is proposed, with appropriate considerations for real fluid effects that involve strong variations in thermodynamic quantities. Our proposed transformation is shown to perform well in collapsing the slope of the logarithmic region to a single universal value with reduced uncertainty. Coupled with a predictive framework to estimate the non-universal shift parameter of the logarithmic region using a priori information, our transformation provides an analytic profile to model the near-wall mean temperature. These results thus provide a framework to guide the development of models for wall-bounded transcritical turbulence.
