Basal melting driven by turbulent thermal convection

Babak Rabbanipour Esfahani, Silvia C. Hirata, Stefano Berti, and Enrico Calzavarini
Phys. Rev. Fluids 3, 053501 – Published 21 May 2018
PDFHTMLExport Citation

Abstract

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énard 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énard 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 of extrapolating results in numerically delicate low-Stefan setups from more convenient high-Stefan ones. Finally, we discuss the implications of our findings for the geophysically relevant problem of modeling Arctic ice melt ponds.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
3 More
  • Received 11 January 2018

DOI:https://doi.org/10.1103/PhysRevFluids.3.053501

©2018 American Physical Society

Physics Subject Headings (PhySH)

Fluid Dynamics

Authors & Affiliations

Babak Rabbanipour Esfahani*, Silvia C. Hirata, Stefano Berti, and Enrico Calzavarini

  • Université de Lille, Unité de Mécanique de Lille, UML EA 7512, F-59000 Lille, France

  • *babak.rabbanipour@polytech-lille.fr

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 3, Iss. 5 — May 2018

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review Fluids

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×