Buoyancy-driven flow regimes for a melting vertical ice cylinder in saline water (2410.22050v2)

Abstract: The presence of salt in seawater significantly affects the melt rate and morphological evolution of ice. This study investigates the melting process of a vertical cylinder in saline water using a combination of laboratory experiments and direct numerical simulations. The two-dimensional (2D) direct numerical simulations and 3D experiments achieve thermal Rayleigh numbers up to $\text{Ra}{T}= \mathcal{O}\left(10^{9}\right)$ and saline Rayleigh numbers up to $\text{Ra}{S}=\mathcal{O}\left(10^{12}\right)$. Some 3D simulations of the vertical ice cylinder are conducted at $\text{Ra}{T}= \mathcal{O}\left(10^{5}\right)$ to confirm that the results in 2D simulations are qualitatively similar to those in 3D simulations. The mean melt rate exhibits a non-monotonic relationship with ambient salinity, initially decreasing towards the point that thermal and saline effects compensate each other, and then increasing again. Based on the ambient salinity, the flow can be categorized into three regimes: temperature-driven flow, salinity-driven flow, and thermal-saline competing flow. In the temperature-driven and competing flow regimes, we find that the mean melt rate follows a $\text{Ra}{T_d}^{1/4}$ scaling. In contrast, in the salinity-driven flow regime, we see a transition from a $\text{Ra}{T_d}^{1/4}$ to a $\text{Ra}{T_d}^{1/3}$ scaling. Additionally, the mean melt rate follows a $\text{Ra}_{S_d}^{1/3}$ scaling in this regime. The ice cylinder develops distinct morphologies in different flow regimes. In the thermal-saline competing flow regime, distinctive scallop (dimpled) patterns emerge along the ice cylinder due to the competition between thermal buoyancy and saline buoyancy. We observe these scallop patterns to migrate downwards over time, due to local differences in the melt rate, for which we provide a qualitative explanation.