The secondary flow and its stability for natural convection in a tall vertical enclosure

1989 ◽  
Vol 200 ◽  
pp. 189-216 ◽  
Author(s):  
Arnon Chait ◽  
Seppo A. Korpela

The multicellular flow between two vertical parallel plates is numerically simulated using a time-splitting pseudospectral method. The steady flow of air, and the time-periodic flow of oil (Prandtl numbers of 0.71 and 1000, respectively) are investigated and descriptions of these flows using both physical and spectral approaches are presented. The details of the time dependency of the flow and temperature fields of oil are shown, and the dynamics of the process is discussed. The spectral transfer of energy among the axial modes comprising the flow is explored. The spectra of kinetic energy and thermal variance for air are found to be smooth and viscously dominated. Similar spectra for oil are bumpier, and the dynamics of the time-dependent flow are determined to be confined to the lower end of the spectrum alone.The three-dimensional linear stability of the multicellular flow of air is parametrically studied. The domain of stable two-dimensional cellular motion was found to be constrained by the Eckhaus instability and by two types of monotone instability. The two-dimensional multicellular flow is unstable above a Grashof number of about 8550 (with the critical Grashof number for the base flow being 8037). Therefore the flow of air in a sufficiently tall vertical enclosure should be considered to be three-dimensional for most practical applications.

2018 ◽  
Vol 851 ◽  
Author(s):  
Shingo Motoki ◽  
Genta Kawahara ◽  
Masaki Shimizu

The divergence-free time-independent velocity field has been determined so as to maximise heat transfer between two parallel plates with a constant temperature difference under the constraint of fixed total enstrophy. The present variational problem is the same as that first formulated by Hassanzadeh et al. (J. Fluid Mech., vol. 751, 2014, pp. 627–662); however, the search range for optimal states has been extended to a three-dimensional velocity field. A scaling of the Nusselt number $Nu$ with the Péclet number $Pe$ (i.e., the square root of the non-dimensionalised enstrophy with thermal diffusion time scale), $Nu\sim Pe^{2/3}$, has been found in the three-dimensional optimal states, corresponding to the asymptotic scaling with the Rayleigh number $Ra$, $Nu\sim Ra^{1/2}$, expected to appear in an ultimate state, and thus to the Taylor energy dissipation law in high-Reynolds-number turbulence. At $Pe\sim 10^{0}$, a two-dimensional array of large-scale convection rolls provides maximal heat transfer. A three-dimensional optimal solution emerges from bifurcation on the two-dimensional solution branch at $Pe\sim 10^{1}$, and the three-dimensional solution branch has been tracked up to $Pe\sim 10^{4}$ (corresponding to $Ra\approx 2.7\times 10^{6}$). At $Pe\gtrsim 10^{3}$, the optimised velocity fields consist of convection cells with hierarchical self-similar vortical structures, and the temperature fields exhibit a logarithmic-like mean profile near the walls.


2019 ◽  
Vol 2019 ◽  
pp. 1-10 ◽  
Author(s):  
Xuan Huang ◽  
Lingfeng Liu ◽  
Xiangjun Li ◽  
Minrong Yu ◽  
Zijie Wu

Given that the sequences generated by logistic map are unsecure with a number of weaknesses, including its relatively small key space, uneven distribution, and vulnerability to attack by phase space reconstruction, this paper proposes a new two-dimensional mutual coupled logistic map, which can overcome these weaknesses. Our two-dimensional chaotic map model is simpler than the recently proposed three-dimensional coupled logistic map, whereas the sequence generated by our system is more complex. Furthermore, a new kind of pseudorandom number generator (PRNG) based on the mutual coupled logistic maps is proposed for application. Both statistical tests and security analysis show that our proposed PRNG has good randomness and that it can resist all kinds of attacks. The algorithm speed analysis indicates that PRNG is valuable to practical applications.


2008 ◽  
Vol 15 (4) ◽  
pp. 695-699 ◽  
Author(s):  
F. Maggi

Abstract. The need to assess the three-dimensional fractal dimension of fractal aggregates from the fractal dimension of two-dimensional projections is very frequent in geophysics, soil, and atmospheric sciences. However, a generally valid approach to relate the two- and three-dimensional fractal dimensions is missing, thus questioning the accuracy of the method used until now in practical applications. A mathematical approach developed for application to suspended aggregates made of cohesive sediment is investigated and applied here more generally to Diffusion-Limited Aggregates (DLA) and Cluster-Cluster Aggregates (CCA), showing higher accuracy in determining the three-dimensional fractal dimension compared to the method currently used.


1998 ◽  
Vol 374 ◽  
pp. 145-171 ◽  
Author(s):  
DANIEL HENRY ◽  
MARC BUFFAT

The convective flows which arise in shallow cavities filled with low-Prandtl-number fluids when subjected to a horizontal temperature gradient are studied numerically with a finite element method. Attention is focused on a rigid cavity with dimensions 4×2×1, for which experimental data are available. The three-dimensional results indicate that, after a relative concentration of the initial Hadley circulation, a transition to time-dependent flows occurs in the form of a roll oscillation with a purely dynamical origin. This transition corresponds to a Hopf bifurcation with a breaking of symmetry that gives some specific properties to the time evolution of the flow: these properties are shown to be the result of the general behaviour of the dynamical systems. Calculations performed in the case of mercury compare well with the experiments with similar power spectra of the temperature, and this validates the analysis of the nature of the global flow performed in the limiting case Pr=0. All these results are discussed with respect to the linear and nonlinear analyses and to other computational experiments. Numerical results obtained in the corresponding two-dimensional situation give a different transition to the time-dependent flow: it is shown that in the three-dimensional cavity this type of two-dimensional transition is less probable than the observed transition with breaking of symmetry.


1991 ◽  
Vol 113 (3) ◽  
pp. 627-634 ◽  
Author(s):  
K. Vafai ◽  
J. Ettefagh

The present work centers around a numerical three-dimensional transient investigation of the effects of axial convection on flow and temperature fields inside an open-ended annulus. The transient behavior of the flow field through the formation of a three-dimensional flow field and its subsequent effect on the temperature distribution at different axial locations within the annulus were analyzed by both finite difference and finite element methods. The results show that the axial convection has a distinctly different influence on the temperature and velocity fields. It is found that in the midportion of the annulus a two-dimensional assumption with respect to the temperature distribution can lead to satisfactory results for Ra<10,000. However, such an assumption is improper with respect to the flow field. Furthermore, it is shown that generally the errors for a two-dimensional assumption in the midportion of the annulus are less at earlier times (t<50Δt) during the transient development of the flow and temperature fields.


2018 ◽  
Vol 841 ◽  
pp. 636-653
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang

A three-dimensional linear stability analysis is carried out for a convecting layer in which both the temperature and solute distributions are linear in the horizontal direction. The three-dimensional results show that, for $Le=3$ and 100, the most unstable mode occurs invariably as the longitudinal mode, a vortex roll with its axis perpendicular to the longitudinal plane, suggesting that the two-dimensional results are sufficient to illustrate the stability characteristics of the convecting layer. Two-dimensional results show that the stability boundaries of the transverse mode (a vortex roll with its axis perpendicular to the transverse plane) and the longitudinal modes are virtually overlapped in the regime dominated by thermal diffusion and the regime dominated by solute diffusion, while these two modes hold a significant difference in the regime the salt-finger instability prevails. More precisely, the instability area in terms of thermal Grashof number $Gr$ and solute Grashof number $Gs$ is larger for the longitudinal mode than the transverse mode, implying that, under any circumstance, the longitudinal mode is always more unstable than the transverse mode.


1996 ◽  
Vol 326 ◽  
pp. 399-415 ◽  
Author(s):  
M. Wanschura ◽  
H. C. Kuhlmann ◽  
H. J. Rath

The stability of steady axisymmetric convection in cylinders heated from below and insulated laterally is investigated numerically using a mixed finite-difference/Chebyshev collocation method to solve the base flow and the linear stability equations. Linear stability boundaries are given for radius to height ratios γ from 0.9 to 1.56 and for Prandtl numbers Pr = 0.02 and Pr = 1. Depending on γ and Pr, the azimuthal wavenumber of the critical mode may be m = 1, 2, 3, or 4. The dependence of the critical Rayleigh number on the aspect ratio and the instability mechanisms are explained by analysing the energy transfer to the critical modes for selected cases. In addition to these results the onset of buoyant convection in liquid bridges with stress-free conditions on the cylindrical surface is considered. For insulating thermal boundary conditions, the onset of convection is never axisymmetric and the critical azimuthal wavenumber increases monotonically with γ. The critical Rayleigh number is less then 1708 for most aspect ratios.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


2021 ◽  
Author(s):  
David Naylor

An introduction is given to the optical setup and principle of operation of classical and holographic interferometers that are used for convective he at transfer measurements. The equations for the evaluation of the temperature field are derived and methods of analysis are discussed for both two-dimensional and three-dimensional temperature fields. Emphasis is given to techniques for measuring local heat transfer rates. For two-dimensional fields, a method is presented for measuring the surface temperature gradient directly from a finite (wedge) fringe interferogram. This “direct gradient method” is shown to be most useful for the measurement of low convective heat transfer rates. For three-dimensional fields, the equations for calculating the beam-averaged local heat flux are presented. The measurement of the fluid temperature averaged along the light beam is shown to be approximate. However, an analysis is presented showing that for most cases the error associated with temperature variations in the light beam direction is small. Digital image analysis of interferograms to obtain fringe spacings is also discussed briefly.


Sign in / Sign up

Export Citation Format

Share Document