Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient Algorithm on a Non-staggered Mesh for Simulating Rayleigh-Benard Convection in a Box
197
0
0.0
(
0
)
Authors
K.-H. Chiam
Ask ChatGPT about the research
No Arabic abstract
An efficient semi-implicit second-order-accurate finite-difference method is described for studying incompressible Rayleigh-Benard convection in a box, with sidewalls that are periodic, thermally insulated, or thermally conducting. Operator-splitting and a projection method reduce the algorithm at each time step to the solution of four Helmholtz equations and one Poisson equation, and these are are solved by fast direct methods. The method is numerically stable even though all field values are placed on a single non-staggered mesh commensurate with the boundaries. The efficiency and accuracy of the method are characterized for several representative convection problems.
rate research
Read More
We use a quantitative topological characterization of complex dynamics to measure geometric structures. This approach is used to analyze the weakly turbulent state of spiral defect chaos in experiments on Rayleigh-Benard convection. Different attractors of spiral defect chaos are distinguished by their homology. The technique reveals pattern asymmetries that are not revealed using statistical measures. In addition we observe global stochastic ergodicity for system parameter values where locally chaotic dynamics has been observed previously.
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.
Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta (non-Oberbeck-Boussinesqness) for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.
We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of non-Oberbeck-Boussinesqness or NOBness ($Delta$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
