The Influence of Real Gas Radiation On the Stability and Development of Benard Convection in a Two-Dimensional Layer

2021 ◽  
Author(s):  
Brent W. Webb ◽  
Vladimir Solovjov
Keyword(s):  
Rayleigh Number ◽  
Critical Rayleigh Number ◽  
Radiative Heating ◽  
Two Dimensional ◽  
Gas Temperature ◽  
Layer Depth ◽  
Real Gas ◽  
Gas Radiation ◽  
Induced Flow ◽  
Buoyancy Induced Flow

Abstract The influence of real gas radiation on the thermal and hydrodynamic stability is investigated in a two-dimensional layer of radiatively participating H2O and/or CO2 heated from below. The non-gray radiation effects of the two species are treated rigorously using a global spectral approach, the Spectral Line Weighted-sum-of-gray-gases model. The phenomena are explored by solving the full coupled laminar equations of motion, energy, and radiative transfer from the low-Rayleigh number, pure conduction-radiation regime through the onset of buoyancy-induced flow to the developed Bénard convection regime. The evolution of the thermal, velocity, and radiative heating fields is studied, and the critical Rayleigh number is characterized as a function of species mole fraction, average layer gas temperature, layer depth, wall emissivity, and the total gas pressure. It is found that participating radiation in the medium has the effect of stabilizing the layer, delaying transition to buoyancy-induced flow. The development of buoyancy-induced flow and temperature, along with the radiative heating are presented. It is found that the critical Rayleigh number in the radiatively participating gas layer can be more than an order of magnitude higher than the classical convection-only scenario. The onset of instability is found to depend on the species mole fractions, average gas temperature in the layer, wall emissivity, layer depth, and total pressure. Generally, all other variables being the same, H2O has a greater stabilizing influence on the layer than CO2.

scholarly journals Linear biglobal analysis of Rayleigh–Bénard instabilities in binary fluids with and without throughflow

2012 ◽  
Vol 713 ◽  
pp. 216-242 ◽  
Author(s):  
Jun Hu ◽  
Daniel Henry ◽  
Xie-Yuan Yin ◽  
Hamda BenHadid
Keyword(s):  
Reynolds Number ◽  
Rayleigh Number ◽  
Aspect Ratio ◽  
Critical Rayleigh Number ◽  
Two Dimensional ◽  
Binary Fluids ◽  
Transverse Modes ◽  
Aspect Ratios ◽  
Rayleigh Benard ◽  
Mode A

AbstractThree-dimensional Rayleigh–Bénard instabilities in binary fluids with Soret effect are studied by linear biglobal stability analysis. The fluid is confined transversally in a duct and a longitudinal throughflow may exist or not. A negative separation factor $\psi = \ensuremath{-} 0. 01$, giving rise to oscillatory transitions, has been considered. The numerical dispersion relation associated with this stability problem is obtained with a two-dimensional Chebyshev collocation method. Symmetry considerations are used in the analysis of the results, which allow the classification of the perturbation modes as ${S}_{l} $ modes (those which keep the left–right symmetry) or ${R}_{x} $ modes (those which keep the symmetry of rotation of $\lrm{\pi} $ about the longitudinal mid-axis). Without throughflow, four dominant pairs of travelling transverse modes with finite wavenumbers $k$ have been found. Each pair corresponds to two symmetry degenerate left and right travelling modes which have the same critical Rayleigh number ${\mathit{Ra}}_{c} $. With the increase of the duct aspect ratio $A$, the critical Rayleigh numbers for these four pairs of modes decrease and closely approach the critical value ${\mathit{Ra}}_{c} = 1743. 894$ obtained in a two-dimensional situation, one of the mode (a ${S}_{l} $ mode called mode A) always remaining the dominant mode. Oscillatory longitudinal instabilities ($k\approx 0$) corresponding to either ${S}_{l} $ or ${R}_{x} $ modes have also been found. Their critical curves, globally decreasing, present oscillatory variations when the duct aspect ratio $A$ is increased, associated with an increasing number of longitudinal rolls. When a throughflow is applied, the symmetry degeneracy of the pairs of travelling transverse modes is broken, giving distinct upstream and downstream modes. For small and moderate aspect ratios $A$, the overall critical Rayleigh number in the small Reynolds number range studied is only determined by the upstream transverse mode A. In contrast, for larger aspect ratios as $A= 7$, different modes are successively dominant as the Reynolds number is increased, involving both upstream and downstream transverse modes A and even the longitudinal mode.

Transition and turbulence in horizontal convection: linear stability analysis

2017 ◽  
Vol 821 ◽  
pp. 31-58 ◽  
Author(s):  
Pierre-Yves Passaggia ◽  
Alberto Scotti ◽  
Brian White
Keyword(s):  
Stability Analysis ◽  
Rayleigh Number ◽  
Boundary Conditions ◽  
Stokes Equations ◽  
Critical Rayleigh Number ◽  
Base Flow ◽  
Two Dimensional ◽  
Slip Boundary Conditions ◽  
Slip Boundary ◽  
Horizontal Convection

The linear instability mechanisms of horizontal convection in a rectangular cavity forced by a horizontal buoyancy gradient along its surface are investigated using local and global stability analyses for a Prandtl number equal to unity. The results show that the stability of the base flow, a steady circulation characterized by a narrow descending plume and a broad upwelling region, depends on the Rayleigh number, $Ra$. For free-slip boundary conditions at a critical value of $Ra\approx 2\times 10^{7}$, the steady base flow becomes unstable to three-dimensional perturbations, characterized by counter-rotating vortices originating within the plume region. A Wentzel–Kramers–Brillouin (WKB) method applied along closed streamlines demonstrates that this instability is of a Rayleigh–Taylor type and can be used to accurately reconstruct the global instability mode. In the case of no-slip boundary conditions, the base flow also becomes unstable to a self-sustained two-dimensional instability whose critical Rayleigh number is $Ra\approx 1.7\times 10^{8}$. Beyond this critical $Ra$, two-dimensional equilibrium stationary states of the Navier–Stokes equations are computed using the selective frequency damping method. The two-dimensional onset of instability is shown to be characterized by a family of modes also originating within the plume. A local spatio-temporal stability analysis shows that the flow becomes absolutely unstable at the origin of the plume. Taken together, these results illustrate the mechanisms behind the onset of turbulence that has been observed in horizontal convection.

Numerical Simulations of an Unsteady Confined Natural Convection Flow

2009 ◽  
Author(s):  
Gillian Leplat ◽  
Emmanuel Laroche ◽  
Philippe Reulet ◽  
Pierre Millan
Keyword(s):  
Natural Convection ◽  
Rayleigh Number ◽  
Large Scale ◽  
Critical Rayleigh Number ◽  
Three Dimensional ◽  
High Amplitude ◽  
Joint Fluid ◽  
Convection Flow ◽  
Two Dimensional ◽  
Natural Convection Flow

A two-dimensional numerical analysis of a laminar natural convection flow within an air-filled enclosure is proposed in this paper from an unstable configuration previously studied experimentally. The flow is driven by a heated square-section cylinder located at the center of a square-section enclosure. Instabilities are observed for an aspect ratio (height of the cylinder over the height of the cavity) of 0.4 and cause the flow to turn into a three-dimensional and unsteady regime characterized by a symmetry breaking and large scale high amplitude flappings around the cylinder. The multi-physic computational software CEDRE, developed at the ONERA, is used to study this unstable behavior and a time-dependent compressible flow solver is used to perform the two-dimensional simulations under the low Mach number approximation, corresponding to the mid-depth cross-section of the enclosure from the experimental configuration. The first results on the investigation of the first unstable modes confirm the onset of the instabilities at the Rayleigh number of the experiment with asymmetrical motions of the fluid around the cylinder. Further analyses highlight the critical Rayleigh number that defines the instability threshold of the first bifurcation which origin and nature could have been identified. Finally, joint fluid-solid simulations are performed to determine more precisely the role of boundary conditions in the onset of instabilities.

scholarly journals LARGE EDDY SIMULATIONS OF HIGH ROSSBY NUMBER FLOW IN THE HIGH PRESSURE COMPRESSOR INTER-DISK CAVITY

2021 ◽  
pp. 1-22
Author(s):  
Deepak Saini ◽  
Richard Sandberg
Keyword(s):  
Rayleigh Number ◽  
Outer Region ◽  
Large Eddy Simulations ◽  
Rossby Number ◽  
Air Cooling ◽  
Large Eddy ◽  
Number Flow ◽  
Induced Flow ◽  
Inner Region ◽  
Buoyancy Induced Flow

Abstract The focus of the present study is to understand the effect of Rayleigh number on a high Rossby number flow in a high pres- sure compressor (HPC) inter-disk cavity. These cavities form between the compressor disks of a gas turbine engine, and they are an integral part of the internal air cooling system. We perform highly resolved large eddy simulations for two Rayleigh numbers of 0.76 × 108 and 1.54 × 108 at a fixed Rossby number of 4.5 by solving the compressible Navier–Stokes equations. The results show a flow structure dominated by a toroidal vortex in the inner region of the cavity. In the outer region, the flow is observed to move radially outwards by Ekman layers formed on the side disks and to move radially inwards through the central core region of the cavity. An enhancement in the intensity of the radial flares is observed in the outer region of the cavity for the high Rayleigh number case with no perceivable effect in the inner region. The near shroud region is mostly dominated by the centrifugal buoyancy-induced flow and the wall Nusselt number calculated at the shroud is in close agreement with centrifugal buoyancy-induced flow without an axial bore flow.

Two Dimensional-Steady State-Natural Convection During the Buoyancy-Induced Flow of CuO-Water Nanofluid Along a Vertical Channel

2014 ◽  
Author(s):  
I. P. Koronaki ◽  
M. T. Nitsas ◽  
Ch. Vallianos
Keyword(s):  
Heat Transfer ◽  
Steady State ◽  
Natural Convection ◽  
Base Fluid ◽  
Dynamic Models ◽  
Effective Conductivity ◽  
Vertical Channel ◽  
Two Dimensional ◽  
Induced Flow ◽  
Buoyancy Induced Flow

In many engineering applications, heat and mass transfer is of vital importance. Therefore a lot of research has been done trying to maximize the heat transfer rate. It is proved, mostly through experimental processes that nanofluids i.e., liquid suspensions of nanometer size particles, have the required capability to augment heat transfer since their efficacy is based on their improved properties compared to those of the base fluid. The present paper examines the two dimensional-steady state-natural convection during the buoyancy-induced flow of the incompressible CuO-water nanofluid along a vertical channel whose walls are uniformly heated. The available literature suggests static and dynamic models for calculating the effective conductivity and viscosity of nanofluids. In this work, both models are assumed so as the Brownian motion of nanoparticles to be considered. The governing equations of continuity, momentum and energy have been solved numerically with the Finite Difference Method (FDM) by using suitable dimensionless variables. The results of the aforementioned analysis prove that the convection coefficient is enhanced due to the presence of nanofluids and it increases further by changing the volume concentration of the nanoparticles. Finally, the effect of the nanoparticles size on heat transfer and the type of the base fluid is investigated.

Interferometric study of two-dimensional Bénard convection cells

1974 ◽  
Vol 66 (4) ◽  
pp. 739-752 ◽  
Author(s):  
R. Farhadieh ◽  
R. S. Tankin
Keyword(s):  
Rayleigh Number ◽  
Sea Water ◽  
Critical Rayleigh Number ◽  
Width Ratio ◽  
Two Dimensional ◽  
Free Boundaries ◽  
Cell Height ◽  
Bénard Convection ◽  
Benard Convection ◽  
Convection Rolls

A Mach-Zehnder interferometer was used to stud two-dimensional Bénard convection cells. The experiments were performed with distilled water and sea water in the region where density is a linear function of temperature. Two-dimensional convection rolls were formed with Rayleigh numbers as great as 23400. Reversal in the temperature profile was obtained for R/Rc ≥ 3·8, and an overshoot of about 6% was observed at R/Rc = 9·2 and 13·8. This agrees with the values predicted theoretically by Veronis (1966) for stress-free boundaries and Royal (1969) for rigid boundaries. This disagrees with the experimental results of Gille (1967), who reports an overshoot of only 1 ½% at R/Rc = 16. Many of the other results agree with those of other experimenters, such as the relation between the cell height-to-width ratio and Rayleigh number, the relation between the Nusselt number and Rayleigh number, and the value of the critical Rayleigh number.

The smooth transition to a convective régime in a two-dimensional box

1978 ◽  
Vol 358 (1693) ◽  
pp. 199-221 ◽  
Keyword(s):  
Rayleigh Number ◽  
Critical Rayleigh Number ◽  
Fluid Motion ◽  
Finite Amplitude ◽  
Linear Stability Theory ◽  
Smooth Transition ◽  
Two Dimensional ◽  
Onset Of Convection ◽  
Free Mode ◽  
Limiting Case

The fluid motion in a two-dimensional box heated from below is considered. The horizontal surfaces are taken to be free and isothermal while the sidewalls are first taken to be rigid and perfect insulators. Linear stability theory shows that the critical Rayleigh number for the onset of convection is higher than that when no side walls are present and the eigenvalue spectrum is discrete. Finite amplitude theory shows that the onset of convection is sudden, that is, bifurcation occurs. The effect of allowing the sidewalls to be slightly imperfect insulators is also investigated. It is found that if the boundary conditions of the sidewalls depart only slightly from those given above, there is a significant change in the response of the fluid. In the most general circumstances a resonance of the free mode is excited as the Rayleigh number approaches its critical value and finite amplitude effects become important. Then it is shown that the onset of convection is quite smooth and the concept of a sharp bifurcation at a critical Rayleigh number is no longer tenable. For a particular class of imperfections it is shown that a ‘transcritical’ bifurcation as described by Benjamin (1976) is possible. The limiting case of a very long box is given special consideration.

Thermoconvective instabilities in a porous medium bounded by two concentric horizontal cylinders

1976 ◽  
Vol 76 (2) ◽  
pp. 337-362 ◽  
Author(s):  
Jean-Paul Caltagirone
Keyword(s):  
Porous Medium ◽  
Finite Elements ◽  
Rayleigh Number ◽  
Critical Rayleigh Number ◽  
Three Dimensional ◽  
Two Dimensional ◽  
Dimensional Model ◽  
Different Types ◽  
Rayleigh Numbers ◽  
Steady Convection

The study of natural convection in a saturated porous medium bounded by two concentric, horizontal, isothermal cylinders reveals different types of evolution according to the experimental conditions and the geometrical configuration of the model. At small Rayleigh numbers the state of the system corresponds to a regime of pseudo-conduction. The isotherms are coaxial with the cylinders. At larger Rayleigh numbers a regime of steady two-dimensional convection sets in between the two cylinders. Finally, for Rayleigh numbers above the critical Rayleigh number Ra*c the phenomena become three-dimensional and fluctuating. The appearance of these different regimes depends, moreover, on the geometry considered and, in particular, on two numbers: R, the ratio of the radii of the cylinders, and A, the ratio of the length of the cylinders to the radius of the inner one. In order to approach these experimental observations and to obtain realistic theoretical models, several methods of solving the equations have been used.The perturbation method yields information about the thermal field and the heat transfer between the cylinders under conditions close to the equilibrium state.A numerical two-dimensional model enables us to extend the range of investigation and to represent properly the phenomena when steady convection appreciably modifies the temperature distribution and the velocities within the porous layer.Neither of these models allows account to be taken of the instabilities observed experimentally above a critical Rayleigh number Ra*c. For this reason, a study of stability has been carried out using a Galerkin method based on equations corresponding to an initial state of steady convection. The results obtained show the importance of three-dimensional effects for the onset of fluctuating convection. The critical transition Rayleigh number Ra*c is thus determined in terms of the ratio of the radii R by solving an eigenvalue problem.A numerical three-dimensional model based on the method of finite elements has thus been developed in order to point out the different types of evolution with time. Steady two-dimensional convection and fluctuating three-dimensional convection have been actually found by calculation. The solution of the system of equations by the method of finite elements is briefly described.The experimental and theoretical results are then compared and a general physical interpretation is given.

Roll-pattern evolution in finite-amplitude Rayleigh–Beńard convection in a two-dimensional fluid layer bounded by distant sidewalls

1984 ◽  
Vol 143 ◽  
pp. 125-152 ◽  
Author(s):  
P. G. Daniels
Keyword(s):  
Rayleigh Number ◽  
Temporal Evolution ◽  
Fluid Layer ◽  
Critical Rayleigh Number ◽  
Finite Amplitude ◽  
Two Dimensional ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

This paper considers the temporal evolution of two-dimensional Rayleigh–Bénard convection in a shallow fluid layer of aspect ratio 2L ([Gt ] 1) confined laterally by rigid sidewalls. Recent studies by Cross et al. (1980, 1983) have shown that for Rayleigh numbers in the range R = R0 + O(L−1) (where R0 is the critical Rayleigh number for the corresponding infinite layer) there exists a class of finite-amplitude steady-state ‘phase-winding’ solutions which correspond physically to the possibility of an adjustment in the number of rolls in the container as the local value of the Rayleigh number is varied. It has been shown (Daniels 1981) that in the temporal evolution of the system the final lateral positioning of the rolls occurs on the long timescale t = O(L2) when the phase function which determines the number of rolls in the system satisfies a one-dimensional diffusion equation but with novel boundary conditions that represent the effect of the sidewalls. In the present paper this system is solved numerically in order to determine the precise way in which the roll pattern adjusts after a change in the Rayleigh number of the system. There is an interesting balance between, on the one hand, a tendency for the number of rolls to change by the least number possible and, on the other, a tendency for the even or odd nature of the initial configuration to be preserved during the transition. In some cases this second property renders the natural evolution susceptible to arbitrarily small external disturbances, which then dictate the form of the final roll pattern.The complete transition involves an analysis of the motion on three timescales, a conductive scale t = O(1), a convective growth scale t = O(L) and a convective diffusion scale t = O(L2).

Sign in / Sign up

Export Citation Format

Share Document