scholarly journals Cessation and reversals of large-scale structures in square Rayleigh–Bénard cells

2019 ◽  
Vol 877 ◽  
pp. 922-954 ◽  
Author(s):  
Andrés Castillo-Castellanos ◽  
Anne Sergent ◽  
Bérengère Podvin ◽  
Maurice Rossi
Keyword(s):  
Phase Space ◽  
Large Scale ◽  
Large Scale Structures ◽  
Multiple Paths ◽  
Data Driven Approach ◽  
Scale Structures ◽  
Prandtl Numbers ◽  
Square Cells ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

We consider direct numerical simulations of turbulent Rayleigh–Bénard convection inside two-dimensional square cells. For Rayleigh numbers $Ra=10^{6}$ to $Ra=5\times 10^{8}$ and Prandtl numbers $Pr=3$ and $Pr=4.3$, two types of flow regimes are observed intermittently: consecutive flow reversals (CR), and extended cessations (EC). For each regime, we combine proper orthogonal decomposition (POD) and statistical tools on long-term data to characterise the dynamics of large-scale structures. For the CR regime, centrosymmetric modes are dominant and display a coherent dynamics, while non-centrosymmetric modes fluctuate randomly. For the EC regime, all POD modes follow Poissonian statistics and a non-centrosymmetric mode is dominant. To explore further the differences between the CR and EC regimes, an analysis based on a cluster partition of the POD phase space is proposed. This data-driven approach confirms the successive mechanisms of the generic reversal cycle in CR as proposed in Castillo-Castellanos et al. (J. Fluid Mech., vol. 808, 2016, pp. 614–640). However, these mechanisms may take one of multiple paths in the POD phase space. Inside the EC regime, this approach reveals the presence of two types of coherent time sequences (weak reversals and actual cessations) and more rarely intense plume crossings. Finally, we analyse within a range of Rayleigh numbers up to turbulent flow, the relation between dynamical regimes and the POD energetic contents as well as the residence time in each cluster.

Numerical simulations of Rayleigh–Bénard convection for Prandtl numbers between 10−1 and 104 and Rayleigh numbers between 105 and 109

2010 ◽  
Vol 662 ◽  
pp. 409-446 ◽  
Author(s):  
G. SILANO ◽  
K. R. SREENIVASAN ◽  
R. VERZICCO
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Nusselt Number ◽  
Numerical Simulations ◽  
Multiple Solutions ◽  
Large Scale ◽  
Prandtl Numbers ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

We summarize the results of an extensive campaign of direct numerical simulations of Rayleigh–Bénard convection at moderate and high Prandtl numbers (10−1 ≤ Pr ≤ 104) and moderate Rayleigh numbers (105 ≤ Ra ≤ 109). The computational domain is a cylindrical cell of aspect ratio Γ = 1/2, with the no-slip condition imposed on all boundaries. By scaling the numerical results, we find that the free-fall velocity should be multiplied by $1/\sqrt{{\it Pr}}$ in order to obtain a more appropriate representation of the large-scale velocity at high Pr. We investigate the Nusselt and the Reynolds number dependences on Ra and Pr, comparing the outcome with previous numerical and experimental results. Depending on Pr, we obtain different power laws of the Nusselt number with respect to Ra, ranging from Ra2/7 for Pr = 1 up to Ra0.31 for Pr = 103. The Nusselt number is independent of Pr. The Reynolds number scales as ${\it Re}\,{\sim}\,\sqrt{{\it Ra}}/{\it Pr}$, neglecting logarithmic corrections. We analyse the global and local features of viscous and thermal boundary layers and their scaling behaviours with respect to Ra and Pr, and with respect to the Reynolds and Péclet numbers. We find that the flow approaches a saturation state when Reynolds number decreases below the critical value, Res ≃ 40. The thermal-boundary-layer thickness increases slightly (instead of decreasing) when the Péclet number increases, because of the moderating influence of the viscous boundary layer. The simulated ranges of Ra and Pr contain steady, periodic and turbulent solutions. A rough estimate of the transition from the steady to the unsteady state is obtained by monitoring the time evolution of the system until it reaches stationary solutions. We find multiple solutions as long-term phenomena at Ra = 108 and Pr = 103, which, however, do not result in significantly different Nusselt numbers. One of these multiple solutions, even if stable over a long time interval, shows a break in the mid-plane symmetry of the temperature profile. We analyse the flow structures through the transitional phases by direct visualizations of the temperature and velocity fields. A wide variety of large-scale circulation and plume structures has been found. The single-roll circulation is characteristic only of the steady and periodic solutions. For other regimes at lower Pr, the mean flow generally consists of two opposite toroidal structures; at higher Pr, the flow is organized in the form of multi-jet structures, extending mostly in the vertical direction. At high Pr, plumes mainly detach from sheet-like structures. The signatures of different large-scale structures are generally well reflected in the data trends with respect to Ra, less in those with respect to Pr.

Thermal convection in inclined cylindrical containers

2016 ◽  
Vol 790 ◽  
Author(s):  
Olga Shishkina ◽  
Susanne Horn
Keyword(s):  
Reynolds Number ◽  
Nusselt Number ◽  
Large Scale ◽  
Direct Numerical Simulations ◽  
Large Scale Circulation ◽  
Prandtl Numbers ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard ◽  
Small Inclination

By means of direct numerical simulations (DNS) we investigate the effect of a tilt angle ${\it\beta}$, $0\leqslant {\it\beta}\leqslant {\rm\pi}/2$, of a Rayleigh–Bénard convection (RBC) cell of aspect ratio 1, on the Nusselt number $\mathit{Nu}$ and Reynolds number $\mathit{Re}$. The considered Rayleigh numbers $\mathit{Ra}$ range from $10^{6}$ to $10^{8}$, the Prandtl numbers range from 0.1 to 100 and the total number of the studied cases is 108. We show that the $\mathit{Nu}\,({\it\beta})/\mathit{Nu}(0)$ dependence is not universal and is strongly influenced by a combination of $\mathit{Ra}$ and $\mathit{Pr}$. Thus, with a small inclination ${\it\beta}$ of the RBC cell, the Nusselt number can decrease or increase, compared to that in the RBC case, for large and small $\mathit{Pr}$, respectively. A slight cell tilt may not only stabilize the plane of the large-scale circulation (LSC) but can also enforce an LSC for cases when the preferred state in the perfect RBC case is not an LSC but a more complicated multiple-roll state. Close to ${\it\beta}={\rm\pi}/2$, $\mathit{Nu}$ and $\mathit{Re}$ decrease with increasing ${\it\beta}$ in all considered cases. Generally, the $\mathit{Nu}({\it\beta})/\mathit{Nu}(0)$ dependence is a complicated, non-monotonic function of ${\it\beta}$.

Turbulent Rotating Rayleigh–Benard Convection: Spatiotemporal and Statistical Study

2008 ◽  
Vol 131 (2) ◽  
Author(s):  
A. Husain ◽  
M. F. Baig ◽  
H. Varshney
Keyword(s):  
Rayleigh Number ◽  
Kinetic Energy ◽  
Rotation Rate ◽  
Large Scale ◽  
Coriolis Forces ◽  
Large Scale Structures ◽  
Vertical Heat ◽  
Scale Structures ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The present study involves a 3D numerical investigation of rotating Rayleigh–Benard convection in a large aspect-ratio (8:8:1) rectangular enclosure. The rectangular cavity is rotated about a vertical axis passing through the center of the cavity. The governing equations of mass, momentum, and energy for a frame rotating with the enclosure, subject to generalized Boussinesq approximation applied to the body and centrifugal force terms, have been solved on a collocated grid using a semi-implicit finite difference technique. The simulations have been carried out for liquid metal flows having a fixed Prandtl number Pr=0.01 and fixed Rayleigh number Ra=107 while rotational Rayleigh number Raw and Taylor number Ta are varied through nondimensional rotation rate (Ω) ranging from 0 to 104. Generation of large-scale structures is observed at low-rotation (Ω=10) rates though at higher-rotation rates (Ω=104) the increase in magnitude of Coriolis forces leads to redistribution of buoyancy-induced vertical kinetic energy to horizontal kinetic energy. This brings about inhibition of vertical fluid transport, thereby leading to reduced vertical heat transfer. The magnitude of rms velocities remains unaffected with an increase in Coriolis forces from Ω=0 to 104. An increase in rotational buoyancy (Raw), at constant rotation rate (Ω=104), on variation in Raw/Ta from 10−3 to 10−2 results in enhanced breakup of large-scale structures with a consequent decrease in rms velocities but with negligible reduction in vertical heat transport.

scholarly journals Large-scale mean patterns in turbulent convection

2015 ◽  
Vol 776 ◽  
pp. 96-108 ◽  
Author(s):  
Mohammad S. Emran ◽  
Jörg Schumacher
Keyword(s):  
Large Scale ◽  
Three Dimensional ◽  
Finite Amplitude ◽  
Mean Flow ◽  
Turbulent Fluctuations ◽  
The Mean ◽  
Prandtl Numbers ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

Large-scale patterns, which are well-known from the spiral defect chaos (SDC) regime of thermal convection at Rayleigh numbers $\mathit{Ra}<10^{4}$, continue to exist in three-dimensional numerical simulations of turbulent Rayleigh–Bénard convection in extended cylindrical cells with an aspect ratio ${\it\Gamma}=50$ and $\mathit{Ra}>10^{5}$. They are revealed when the turbulent fields are averaged in time and turbulent fluctuations are thus removed. We apply the Boussinesq closure to estimate turbulent viscosities and diffusivities, respectively. The resulting turbulent Rayleigh number $\mathit{Ra}_{\ast }$, that describes the convection of the mean patterns, is indeed in the SDC range. The turbulent Prandtl numbers are smaller than one, with $0.2\leqslant \mathit{Pr}_{\ast }\leqslant 0.4$ for Prandtl numbers $0.7\leqslant \mathit{Pr}\leqslant 10$. Finally, we demonstrate that these mean flow patterns are robust to an additional finite-amplitude sidewall forcing when the level of turbulent fluctuations in the flow is sufficiently high.

Some comments about observations and image processing of comet 29P/Schwassmann-Wachmann 1

1999 ◽  
Vol 173 ◽  
pp. 243-248
Author(s):  
D. Kubáček ◽  
A. Galád ◽  
A. Pravda
Keyword(s):  
Image Processing ◽  
Large Scale ◽  
Harvard College ◽  
Oak Ridge ◽  
Large Scale Structures ◽  
Short Period Comet ◽  
Short Period ◽  
Scale Structures ◽  
Harvard College Observatory ◽  
Period Comet

AbstractUnusual short-period comet 29P/Schwassmann-Wachmann 1 inspired many observers to explain its unpredictable outbursts. In this paper large scale structures and features from the inner part of the coma in time periods around outbursts are studied. CCD images were taken at Whipple Observatory, Mt. Hopkins, in 1989 and at Astronomical Observatory, Modra, from 1995 to 1998. Photographic plates of the comet were taken at Harvard College Observatory, Oak Ridge, from 1974 to 1982. The latter were digitized at first to apply the same techniques of image processing for optimizing the visibility of features in the coma during outbursts. Outbursts and coma structures show various shapes.

Sign in / Sign up

Export Citation Format

Share Document