Melting driven by rotating Rayleigh–Bénard convection

Ravichandran, S.; Wettlaufer, J. S.

doi:10.1017/jfm.2021.223

Cited by 12 publications

(35 citation statements)

References 53 publications

(139 reference statements)

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…Furthermore, the results presented here were ascertained to be grid-independent by doubling the number of grid points in the vertical direction for the smallest Ri b simulations. The penalization parameter is fixed at η ≈ 1.4 × 10 −3 , which is sufficiently small for results to be independent of this parameter [41], and the time-step used is ∆t ≤ 7 × 10 −4 .…”

Section: Resultsmentioning

confidence: 99%

“…Equations 11 -13, along with the boundary conditions, are solved using the finite-volume solver Megha-5, with a volume penalization method for the melting solid. The solver uses second-order central differences in space and a second-order Adams-Bashforth time-stepping scheme, and has been extensively validated for free-shear flows [44,45] and Rayleigh-Bénard convection [41,46,47]. Details of the volume-penalization algorithm, including validation against the analytical solution for the purely conductive Stefan problem, tests of sensitivity to the penalization parameter, and the convergence of the solution under grid-refinement may be found in Ravichandran and Wettlaufer [41].…”

Section: Methodsmentioning

confidence: 99%

“…The solver uses second-order central differences in space and a second-order Adams-Bashforth time-stepping scheme, and has been extensively validated for free-shear flows [44,45] and Rayleigh-Bénard convection [41,46,47]. Details of the volume-penalization algorithm, including validation against the analytical solution for the purely conductive Stefan problem, tests of sensitivity to the penalization parameter, and the convergence of the solution under grid-refinement may be found in Ravichandran and Wettlaufer [41]. We have further verified that the body force prescribed in equation 1 produces the known analytical plane-Poiseuille velocity profile with an exponential relaxation timescale of (π 2 ν) −1 for the peak velocity.…”

Section: Methodsmentioning

confidence: 99%

“…An important means of quantifying the response of the system is in terms of the Nusselt number, which is defined as the ratio of total heat flux to the heat flux due only to conduction. Following [41], we define the Nusselt number of melting as…”

Section: F Non-dimensional Equations Of Motionmentioning

confidence: 99%

“…Ravichandran and Wettlaufer [41] explored the combined effects of rotation and thermal convection on the evolution of a phase boundary. They found that the rotation-dominated convection of a pure liquid, which takes the form of columnar vortices [42], melts voids into its overlying solid phase.…”

Section: Introductionmentioning

confidence: 99%

See 4 more Smart Citations

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

Ravichandran,

Toppaladoddi,

Wettlaufer

2021

Preprint

Self Cite

View full text Add to dashboard Cite

We use well-resolved numerical simulations to study the combined effects of buoyancy, pressuredriven shear and rotation on the melt rate and morphology of a layer of pure solid overlying its liquid phase in three dimensions at a Rayleigh number Ra = 1.25×10 5 . During thermal convection we find that the rate of melting of the solid phase varies non-monotonically with the strength of the imposed shear flow. In the absence of rotation, depending on whether buoyancy or shear is dominant, we observe either domes or ridges aligned in the direction of the shear flow respectively. Furthermore, we show that the geometry of the phase boundary has important effects on the magnitude and evolution of the heat flux in the liquid layer. In the presence of rotation, the strength of which is characterized by the Rossby number (Ro), we observe that for Ro = O(1) the mean flow in the interior is perpendicular to the direction of the constant applied pressure gradient. As the magnitude of the applied horizontal pressure gradient increases, the geometry of solid-liquid interface evolves from the voids characteristic of melting by rotating convection, to grooves oriented perpendicular to or aligned at an oblique angle to the applied pressure gradient.

show abstract

Section: Resultsmentioning

confidence: 99%

Section: Methodsmentioning

confidence: 99%

Section: Methodsmentioning

confidence: 99%

Section: F Non-dimensional Equations Of Motionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

Ravichandran,

Toppaladoddi,

Wettlaufer

2021

Preprint

Self Cite

View full text Add to dashboard Cite

show abstract

Experimental study on heat transfer characteristics of ice melting processes under point-source bubble flows

Liu,

Zhang,

Song

et al. 2024

International Communications in Heat and Mass Transfer

View full text Add to dashboard Cite

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

2022

Self Cite

View full text Add to dashboard Cite

We use well-resolved numerical simulations to study the combined effects of buoyancy, pressure-driven shear and rotation on the melt rate and morphology of a layer of pure solid overlying its liquid phase in three dimensions at a Rayleigh number $Ra=1.25\times 10^5$ . During thermal convection, we find that the rate of melting of the solid phase varies non-monotonically with the strength of the imposed shear flow. In the absence of rotation, depending on whether buoyancy or shear dominates the flow, we observe either domes or ridges aligned in the direction of the shear flow, respectively. Furthermore, we show that the geometry of the phase boundary has important effects on the magnitude and evolution of the heat flux in the liquid layer. In the presence of rotation, the strength of which is characterized by the Rossby number, $Ro$ , we observe that for $Ro={O}(1)$ , the mean flow in the interior is perpendicular to the direction of the constant horizontal applied pressure gradient. As the magnitude of this pressure gradient increases, the geometry of solid–liquid interface evolves from the voids characteristic of melting by rotating convection, to grooves oriented perpendicular or obliquely to the direction of the pressure gradient.

show abstract

Melting driven by rotating Rayleigh–Bénard convection

Abstract: Abstract

Cited by 12 publications

References 53 publications

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

Experimental study on heat transfer characteristics of ice melting processes under point-source bubble flows

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

Contact Info

Product

Resources

About