ترغب بنشر مسار تعليمي؟ اضغط هنا

Melting driven by rotating Rayleigh-Benard convection

118   0   0.0 ( 0 )
 نشر من قبل Sivaramakrishnan Ravichandran
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We study numerically the melting of a horizontal layer of a pure solid above a convecting layer of its fluid rotating about the vertical axis. In the rotating regime studied here, with Rayleigh numbers of order $10^7$, convection takes the form of columnar vortices, the number and size of which depend upon the Ekman and Prandtl numbers, as well as the geometry -- periodic or confined. As the Ekman and Rayleigh numbers vary, the number and average area of vortices vary in inverse proportion, becoming thinner and more numerous with decreasing Ekman number. The vortices transport heat to the phase boundary thereby controlling its morphology, characterized by the number and size of the voids formed in the solid, and the overall melt rate, which increases when the lower boundary is governed by a no-slip rather than a stress-free velocity boundary condition. Moreover, the number and size of voids formed are relatively insensitive to the Stefan number, here inversely proportional to the latent heat of fusion. For small values of the Stefan number, the convection in the fluid reaches a slowly evolving geostrophic state wherein columnar vortices transport nearly all the heat from the lower boundary to melt the solid at an approximately constant rate. In this quasi-steady state, we find that the Nusselt number, characterizing the heat flux, co-varies with the interfacial roughness, for all the flow parameters and Stefan numbers considered here. This confluence of processes should influence the treatment of moving boundary problems, particularly those in astrophysical and geophysical problems where rotational effects are important.



قيم البحث

اقرأ أيضاً

Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh-Benard convection. In the low-latitude region $mathrm{I}$ vertical (parallel to the axis of rotation) convective columns are ge nerated between the hot inner and the cold outer sphere. The mid-latitude region $mathrm{II}$ is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion-free scaling, which indicates bulk-dominated convection, originates from this mid-latitude region. In the equator region $mathrm{III}$ the vortices are affected by the outer spherical boundary and are much shorter than in region $mathrm{II}$. Thermally driven turbulence with background rotation in spherical Rayleigh-Benard convection is found to be characterized by three distinctly different flow regions. The diffusion-free scaling, which indicates the heat transfer is bulk-dominated, originates from the mid-latitude region in which vertically aligned vortices are stretched between the Northern and Southern hemispheres of the outer sphere. These results show that the flow physics in rotating convection are qualitatively different in planar and spherical geometries. This finding underlines that it is crucial to study convection in spherical geometries to better understand geophysical and astrophysical flow phenomena.
The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a relatively low Rayleigh number, i.e. $Ra = 4times10^7$, is studied for different Pr by direct numerical simulations and the results are compared with laminar BL theory. In this regime we find a smooth onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer enhancement as function of the rotation rate for this relatively low Ra number regime and weak background rotation $Rogtrsim 1$. The smooth onset of heat transfer enhancement observed here is in contrast to the sharp onset observed at larger $Ra gtrsim 10^8$ by Stevens {it{et al.}} [Phys. Rev. Lett. {bf{103}}, 024503, 2009], although only a small shift in the Ra-Ro-Pr phase space is involved.
For rapidly rotating turbulent Rayleigh--Benard convection in a slender cylindrical cell, experiments and direct numerical simulations reveal a boundary zonal flow (BZF) that replaces the classical large-scale circulation. The BZF is located near the vertical side wall and enables enhanced heat transport there. Although the azimuthal velocity of the BZF is cyclonic (in the rotating frame), the temperature is an anticyclonic traveling wave of mode one whose signature is a bimodal temperature distribution near the radial boundary. The BZF width is found to scale like $Ra^{1/4}Ek^{2/3}$ where the Ekman number $Ek$ decreases with increasing rotation rate.
Using direct numerical simulations, we study rotating Rayleigh-Benard convection in a cylindrical cell for a broad range of Rayleigh, Ekman, and Prandtl numbers from the onset of wall modes to the geostrophic regime, an extremely important one in geo physical and astrophysical contexts. We connect linear wall-mode states that occur prior to the onset of bulk convection with the boundary zonal flow that coexists with turbulent bulk convection in the geostrophic regime through the continuity of length and time scales and of convective heat transport. We quantitatively collapse drift frequency, boundary length, and heat transport data from numerous sources over many orders of magnitude in Rayleigh and Ekman numbers. Elucidating the heat transport contributions of wall modes and of the boundary zonal flow are critical for characterizing the properties of the geostrophic regime of rotating convection in finite, physical containers and is crucial for connecting the geostrophic regime of laboratory convection with geophysical and astrophysical systems.
Steady flows that optimize heat transport are obtained for two-dimensional Rayleigh-Benard convection with no-slip horizontal walls for a variety of Prandtl numbers $Pr$ and Rayleigh number up to $Rasim 10^9$. Power law scalings of $Nusim Ra^{gamma}$ are observed with $gammaapprox 0.31$, where the Nusselt number $Nu$ is a non-dimensional measure of the vertical heat transport. Any dependence of the scaling exponent on $Pr$ is found to be extremely weak. On the other hand, the presence of two local maxima of $Nu$ with different horizontal wavenumbers at the same $Ra$ leads to the emergence of two different flow structures as candidates for optimizing the heat transport. For $Pr lesssim 7$, optimal transport is achieved at the smaller maximal wavenumber. In these fluids, the optimal structure is a plume of warm rising fluid which spawns left/right horizontal arms near the top of the channel, leading to downdrafts adjacent to the central updraft. For $Pr > 7$ at high-enough Ra, the optimal structure is a single updraft absent significant horizontal structure, and characterized by the larger maximal wavenumber.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا