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

Bifurcation analysis of two-dimensional Rayleigh--Benard convection using deflation

127   0   0.0 ( 0 )
 نشر من قبل Nicolas Boull\\'e
 تاريخ النشر 2021
والبحث باللغة English




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

We perform a bifurcation analysis of the steady state solutions of Rayleigh--Benard convection with no-slip boundary conditions in two dimensions using a numerical method called deflated continuation. By combining this method with an initialisation strategy based on the eigenmodes of the conducting state, we are able to discover multiple solutions to this non-linear problem, including disconnected branches of the bifurcation diagram, without the need of any prior knowledge of the dynamics. One of the disconnected branches we find contains a s-shape bifurcation with hysteresis, which is the origin of the flow pattern that may be related to the dynamics of flow reversals in the turbulent regime. Linear stability analysis is also performed to analyse the steady and unsteady regimes of the solutions in the parameter space and to characterise the type of instabilities.



قيم البحث

اقرأ أيضاً

We analyse the nonlinear dynamics of the large scale flow in Rayleigh-Benard convection in a two-dimensional, rectangular geometry of aspect ratio $Gamma$. We impose periodic and free-slip boundary conditions in the streamwise and spanwise directions , respectively. As Rayleigh number Ra increases, a large scale zonal flow dominates the dynamics of a moderate Prandtl number fluid. At high Ra, in the turbulent regime, transitions are seen in the probability density function (PDF) of the largest scale mode. For $Gamma = 2$, the PDF first transitions from a Gaussian to a trimodal behaviour, signifying the emergence of reversals of the zonal flow where the flow fluctuates between three distinct turbulent states: two states in which the zonal flow travels in opposite directions and one state with no zonal mean flow. Further increase in Ra leads to a transition from a trimodal to a unimodal PDF which demonstrates the disappearance of the zonal flow reversals. On the other hand, for $Gamma = 1$ the zonal flow reversals are characterised by a bimodal PDF of the largest scale mode, where the flow fluctuates only between two distinct turbulent states with zonal flow travelling in opposite directions.
We study the stability of steady convection rolls in 2D Rayleigh--Benard convection with free-slip boundaries and horizontal periodicity over twelve orders of magnitude in the Prandtl number $(10^{-6} leq Pr leq 10^6)$ and five orders of magnitude in the Rayleigh number $(8pi^4 < Ra leq 3 times 10^7)$. The analysis is facilitated by partitioning our modal expansion into so-called even and odd modes. With aspect ratio $Gamma = 2$, we observe that zonal modes (with horizontal wavenumber equal to zero) can emerge only once the steady convection roll state consisting of even modes only becomes unstable to odd perturbations. We determine the stability boundary in the $(Pr,Ra)$-plane and observe remarkably intricate features corresponding to qualitative changes in the solution, as well as three regions where the steady convection rolls lose and subsequently regain stability as the Rayleigh number is increased. We study the asymptotic limit $Pr to 0$ and find that the steady convection rolls become unstable almost instantaneously, eventually leading to non-linear relaxation osculations and bursts, which we can explain with a weakly non-linear analysis. In the complementary large-$Pr$ limit, we observe that the stability boundary reaches an asymptotic value $Ra = 2.54 times 10^7$ and that the zonal modes at the instability switch off abruptly at a large, but finite, Prandtl number.
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$.
In this numerical study on two-dimensional Rayleigh-Benard convection we consider $10^7 leq Ra leq 10^{12}$ in aspect ratio $0.23 leq Gamma leq 13$ samples. We focus on several cases. First we consider small aspect ratio cells, where at high Ra numbe r we find a sharp transition from a low Ra number branch towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show that the influence of the aspect ratio on the heat transport decreases with increasing aspect ratio, although even at very large aspect ratio of $Gammaapprox10$ variations up to 2.5% in the heat transport as a function of Gamma are observed. Finally, we observe long-lived transients up to at least $Ra=10^9$, as in certain aspect ratio cells we observe different flow states that are stable for thousands of turnover times.
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

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