scholarly journals A proof that convection in a porous vertical slab may be unstable

2015 ◽  
Vol 770 ◽  
pp. 273-288 ◽  
Author(s):  
A. Barletta

The stability of the natural convection parallel flow in a vertical porous slab is reconsidered, by reformulating the problem originally solved in Gill’s classical paper of 1969 (J. Fluid Mech., vol. 35, pp. 545–547). A comparison is made between the set of boundary conditions where the slab boundaries are assumed to be isothermal and impermeable (Model A), and the set of boundary conditions where the boundaries are modelled as isothermal and permeable (Model B). It is shown that Gill’s proof of linear stability for Model A cannot be extended to Model B. The question about the stability/instability of the basic flow is examined by carrying out a numerical solution of the stability eigenvalue problem. It is shown that, with Model B, the natural convection parallel flow in the basic state becomes unstable when the Darcy–Rayleigh number is larger than 197.081. The normal modes selected at onset of instability are transverse rolls. Direct numerical simulations of the nonlinear regime of instability are carried out.

2014 ◽  
Vol 742 ◽  
pp. 636-663 ◽  
Author(s):  
P. Ripesi ◽  
L. Biferale ◽  
M. Sbragaglia ◽  
A. Wirth

AbstractWe investigate the stability and dynamics of natural convection in two dimensions, subject to inhomogeneous boundary conditions. In particular, we consider a Rayleigh–Bénard (RB) cell, where the horizontal top boundary contains a periodic sequence of alternating thermal insulating and conducting patches, and we study the effects of the heterogeneous pattern on the global heat exchange, at both low and high Rayleigh numbers. At low Rayleigh numbers, we determine numerically the transition from a regime characterized by the presence of small convective cells localized at the inhomogeneous boundary to the onset of ‘bulk’ convective rolls spanning the entire domain. Such a transition is also controlled analytically in the limit when the boundary pattern length is small compared with the cell vertical size. At higher Rayleigh number, we use numerical simulations based on a lattice Boltzmann method to assess the impact of boundary inhomogeneities on the fully turbulent regime up to $\mathit{Ra} \sim 10^{10}$.


2016 ◽  
Vol 64 (1) ◽  
pp. 31-37 ◽  
Author(s):  
Roushanara Begum ◽  
MZI Bangalee

Effects of different boundary conditions at the surfaces of the extended computational domain on buoyancy driven natural convection flow in a three dimensional open cavity are studied numerically. This study is carried out for turbulent flow where Rayleigh number is greater than 108. Air is used as working fluid having properties at 25°C temperature and 1atm pressure. To capture the turbulent nature of the flow k - ? model is used. ANSYS CFX software is used to solve the governing equations subject to the corresponding boundary conditions. The methodology is verified through a satisfactory comparison with some published results. Average mass flow, temperature, stream line, contour velocity and velocity profile are studied at different height. An extended computational domain around the physical domain of the cavity at different surrounding conditions is considered to investigate the effect of its existence on the computation. Effects of different surrounding boundary conditions on the physical domain of the cavity are studied and reported.A relation among non-dimensional parameters such as Nusselt number, Rayleigh number, Prandlt number and Reynolds number is also reported.Dhaka Univ. J. Sci. 64(1): 31-37, 2016 (January)


1998 ◽  
Vol 120 (3) ◽  
pp. 776-783 ◽  
Author(s):  
J. Melanson ◽  
J. W. Zu

Vibration analysis of an internally damped rotating shaft, modeled using Timoshenko beam theory, with general boundary conditions is performed analytically. The equations of motion including the effects of internal viscous and hysteretic damping are derived. Exact solutions for the complex natural frequencies and complex normal modes are provided for each of the six classical boundary conditions. Numerical simulations show the effect of the internal damping on the stability of the rotor system.


Author(s):  
Weilong Wang ◽  
Jilian Wu ◽  
Xinlong Feng

Purpose The purpose of this paper is to propose a new method to solve the incompressible natural convection problem with variable density. The main novel ideas of this work are to overcome the stability issue due to the nonlinear inertial term and the hyperbolic term for conventional finite element methods and to deal with high Rayleigh number for the natural convection problem. Design/methodology/approach The paper introduces a novel characteristic variational multiscale (C-VMS) finite element method which combines advantages of both the characteristic and variational multiscale methods within a variational framework for solving the incompressible natural convection problem with variable density. The authors chose the conforming finite element pair (P2, P2, P1, P2) to approximate the density, velocity, pressure and temperature field. Findings The paper gives the stability analysis of the C-VMS method. Extensive two-dimensional/three-dimensional numerical tests demonstrated that the C-VMS method not only can deal with the incompressible natural convection problem with variable density but also with high Rayleigh number very well. Originality/value Extensive 2D/3D numerical tests demonstrated that the C-VMS method not only can deal with the incompressible natural convection problem with variable density but also with high Rayleigh number very well.


2007 ◽  
Vol 586 ◽  
pp. 259-293 ◽  
Author(s):  
F. X. TRIAS ◽  
M. SORIA ◽  
A. OLIVA ◽  
C. D. PÉREZ-SEGARRA

A set of complete two- and three-dimensional direct numerical simulations (DNS) in a differentially heated air-filled cavity of aspect ratio 4 with adiabatic horizontal walls is presented in this paper. Although the physical phenomenon is three-dimensional, owing to its prohibitive computational costs the majority of the previous DNS of turbulent and transition natural convection flows in enclosed cavities assumed a two-dimensional behaviour. The configurations selected here (Rayleigh number based on the cavity height 6.4 × 108, 2 × 109 and 1010, Pr = 0.71) are an extension to three dimensions of previous two-dimensional problems.An overview of the numerical algorithm and the methodology used to verify the code and the simulations is presented. The main features of the flow, including the time-averaged flow structure, the power spectra and probability density distributions of a set of selected monitoring points, the turbulent statistics, the global kinetic energy balances and the internal waves motion phenomenon are described and discussed.As expected, significant differences are observed between two- and three-dimensional results. For two-dimensional simulations the oscillations at the downstream part of the vertical boundary layer are clearly stronger, ejecting large eddies to the cavity core. In the three-dimensional simulations these large eddies do not persist and their energy is rapidly passed down to smaller scales of motion. It yields on a reduction of the large-scale mixing effect at the hot upper and cold lower regions and consequently the cavity core still remains almost motionless even for the highest Rayleigh number. The boundary layers remain laminar in their upstream parts up to the point where these eddies are ejected. The point where this phenomenon occurs clearly moves upstream for the three-dimensional simulations. It is also shown that, even for the three-dimensional simulations, these eddies are large enough to permanently excite an internal wave motion in the stratified core region. All these differences become more marked for the highest Rayleigh number.


2011 ◽  
Vol 133 (4) ◽  
Author(s):  
Suvash C. Saha ◽  
John C. Patterson ◽  
Chengwang Lei

The natural convection boundary layer adjacent to an inclined plate subject to sudden cooling boundary condition has been studied. It is found that the cold boundary layer adjacent to the plate is potentially unstable to Rayleigh–Bénard instability if the Rayleigh number exceeds a certain critical value. A scaling relation for the onset of instability of the boundary layer is achieved. The scaling relations have been developed by equating important terms of the governing equations based on the development of the boundary layer with time. The flow adjacent to the plate can be classified broadly into a conductive, a stable convective, or an unstable convective regime determined by the Rayleigh number. Proper scales have been established to quantify the flow properties in each of these flow regimes. An appropriate identification of the time when the instability may set in is discussed. A numerical verification of the time for the onset of instability is also presented in this study. Different flow regimes based on the stability of the boundary layer have been discussed with numerical results.


2006 ◽  
Vol 63 (11) ◽  
pp. 2795-2812 ◽  
Author(s):  
Christopher L. Wolfe ◽  
Roger M. Samelson

Abstract The stability of a time-periodic baroclinic wave-mean oscillation in a high-dimensional two-layer quasigeostrophic spectral model is examined by computing a full set of time-dependent normal modes (Floquet vectors) for the oscillation. The model has 72 × 62 horizontal resolution and there are 8928 Floquet vectors in the complete set. The Floquet vectors fall into two classes that have direct physical interpretations: wave-dynamical (WD) modes and damped-advective (DA) modes. The WD modes (which include two neutral modes related to continuous symmetries of the underlying system) have large scales and can efficiently exchange energy and vorticity with the basic flow; thus, the dynamics of the WD modes reflects the dynamics of the wave-mean oscillation. These modes are analogous to the normal modes of steady parallel flow. On the other hand, the DA modes have fine scales and dynamics that reduce, to first order, to damped advection of the potential vorticity by the basic flow. While individual WD modes have immediate physical interpretations as discrete normal modes, the DA modes are best viewed, in sum, as a generalized solution to the damped advection problem. The asymptotic stability of the time-periodic basic flow is determined by a small number of discrete WD modes and, thus, the number of independent initial disturbances, which may destabilize the basic flow, is likewise small. Comparison of the Floquet exponent spectrum of the wave-mean oscillation to the Lyapunov exponent spectrum of a nearby aperiodic trajectory suggests that this result will still be obtained when the restriction to time periodicity is relaxed.


2011 ◽  
Vol 21 (06) ◽  
pp. 1539-1582 ◽  
Author(s):  
TASSOS BOUNTIS ◽  
GEORGE CHECHIN ◽  
VLADIMIR SAKHNENKO

In the present tutorial we address a problem with a long history, which remains of great interest to date due to its many important applications: It concerns the existence and stability of periodic and quasiperiodic orbits in N-degree of freedom Hamiltonian systems and their connection with discrete symmetries. Of primary importance in our study is what we call nonlinear normal modes (NNMs), i.e. periodic solutions which represent continuations of the system's linear normal modes in the nonlinear regime. We examine questions concerning the existence of such solutions and discuss different methods for constructing them and studying their stability under fixed and periodic boundary conditions. In the periodic case, we find it particularly useful to approach the problem through the discrete symmetries of many models, employing group theoretical concepts to identify a special type of NNMs which we call one-dimensional "bushes". We then describe how to use linear combinations of s ≥ 2 such NNMs to construct s-dimensional bushes of quasiperiodic orbits, for a wide variety of Hamiltonian systems including particle chains, a square molecule and octahedral crystals in 1, 2 and 3 dimensions. Next, we exploit the symmetries of the linearized equations of motion about these bushes to demonstrate how they may be simplified to study the destabilization of these orbits, as a result of their interaction with NNMs not belonging to the same bush. Applying this theory to the famous Fermi Pasta Ulam (FPU) chain, we review a number of interesting results concerning the stability of NNMs and higher-dimensional bushes, which have appeared in the recent literature. We then turn to a newly developed approach to the analytical and numerical construction of quasiperiodic orbits, which does not depend on the symmetries or boundary conditions of our system. Using this approach, we demonstrate that the well-known "paradox" of FPU recurrences may in fact be explained in terms of the exponential localization of the energies Eq of NNM's being excited at the low part of the frequency spectrum, i.e. q = 1, 2, 3, …. These results indicate that it is the stability of these low-dimensional compact manifolds called q-tori, that is related to the persistence or FPU recurrences at low energies. Finally, we discuss a novel approach to the stability of orbits of conservative systems, expressed by a spectrum of indices called GALI k, k = 2, …, 2N, by means of which one can determine accurately and efficiently the destabilization of q-tori, leading, after very long times, to the breakdown of recurrences and, ultimately, to the equipartition of energy, at high enough values of the total energy E.


1970 ◽  
Vol 42 (1) ◽  
pp. 125-127 ◽  
Author(s):  
A. E. Gill ◽  
C. C. Kirkham

The study of natural convection in a rectangular cavity whose side walls are maintained at different fixed temperatures can be regarded as one of the classical problems of thermal convection (see Batchelor 1954; Elder 1965 and Gill 1966). One aspect of this problem is the question of stability of the laminar flow solution for given values of the three governing parameters, the Rayleigh numberA= γgΔTL3/κν, the aspect ratioh=H/Land the Prandtl number σ = ν/κ. HereHis the height andLthe width of the cavity,Tis the temperature difference between the two walls andgthe acceleration due to gravity. The fluid filling the container has coefficient of expansion γ, thermal diffusivity κ and kinematic viscosity ν. Instead of studying the stability of the exact laminar flow solution, which is not a parallel flow, the stability of the following solution of the governing eauations is examined instead:


Sign in / Sign up

Export Citation Format

Share Document