Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConnecting Lyapunov Vectors with the Pattern Dynamics of Chaotic Rayleigh-Benard Convection
158
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We explore the chaotic dynamics of Rayleigh-Benard convection using large-scale, parallel numerical simulations for experimentally accessible conditions. We quantify the connections between the spatiotemporal dynamics of the leading-order Lyapunov vector and different measures of the flow field patterns topology and dynamics. We use a range of pattern diagnostics to describe the spatiotemporal features of the flow field structures which includes many of the traditional diagnostics used to describe convection as well as some diagnostics tailored to capture the dynamics of the patterns. Using precision-recall curves, we quantify the complex relationship between the pattern diagnostics and the regions where the magnitude of the leading-order Lyapunov vector is significant.
rate research
Read More
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 geophysical 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.
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.
Using direct numerical simulations, we study the statistical properties of reversals in two-dimensional Rayleigh-Benard convection for infinite Prandtl number. We find that the large-scale circulation reverses irregularly, with the waiting time between two consecutive genuine reversals exhibiting a Poisson distribution on long time scales, while the interval between successive crossings on short time scales shows a power law distribution. We observe that the vertical velocities near the sidewall and at the center show different statistical properties. The velocity near the sidewall shows a longer autocorrelation and $1/f^2$ power spectrum for a wide range of frequencies, compared to shorter autocorrelation and a narrower scaling range for the velocity at the center. The probability distribution of the velocity near the sidewall is bimodal, indicating a reversing velocity field. We also find that the dominant Fourier modes capture the dynamics at the sidewall and at the center very well. Moreover, we show a signature of weak intermittency in the fluctuations of velocity near the sidewall by computing temporal structure functions.
Rayleigh-Benard (RB) convection with free-slip plates and horizontally periodic boundary conditions is investigated using direct numerical simulations. Two configurations are considered, one is two-dimension (2D) RB convection and the other one three-dimension (3D) RB convection with a rotating axis parallel to the plate. We explore the parameter range of Rayleigh numbers Ra from $10^7 to $10^9$ and Prandtl numbers $Pr$ from $1$ to $100$. We show that zonal flow, which was observed, for example, by Goluskin emph{et al}. emph{J. Fluid. Mech.} 759, 360-385 (2014) for $Gamma=2$, is only stable when $Gamma$ is smaller than a critical value, which depends on $Ra$ and $Pr$. With increasing $Gamma$, we find a second regime in which both zonal flow and different convection roll states can be statistically stable. For even larger $Gamma$, in a third regime, only convection roll states are statistically stable and zonal flow is not sustained. For the 3D simulations, we fix $Ra=10^7$ and $Pr=0.71$, and compare the flow for $Gamma=8$ and $Gamma = 16$. We demonstrate that with increasing aspect ratio $Gamma$, zonal flow, which was observed for small $Gamma=2pi$ by von Hardenberg emph{et al}. emph{Phys. Rev. Lett.} 15, 134501 (2015), completely disappears for $Gamma=16$. For such large $Gamma$ only convection roll states are statistically stable. In between, here for medium aspect ratio $Gamma = 8$, the convection roll state and the zonal flow state are both statistically stable. What state is taken depends on the initial conditions, similarly as we found for the 2D case.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
