The onset of thermo-bioconvection in a shallow fluid saturated porous layer heated from below in a suspension of oxytactic microorganisms

The aim of this paper is to present a continuum model for thermo-bioconvection of oxytactic bacteria in a porous medium and investigate the combined effects of microorganisms' upswimming and heating from below on the stability of bioconvection in a horizontal layer filled with a fluid saturated porous medium. Different from traditional bioconvection, thermo-bioconvection has two destabilizing mechanisms that contribute to creating the unstable density stratification. This problem may be relevant to a number of geophysical applications, such as the investigation of the dynamics of oxytactic species of thermophiles (heat loving microorganisms) living in hot springs, microbial-enhanced oil recovery, and modeling oil- and gas-bearing sedimentary basins. The utilization of the Galerkin method to solve a linear stability problem leads to a correlation between the critical value of the bioconvection Rayleigh number and the traditional “thermal” Rayleigh number.     

Overstability Analysis of Thermo-bioconvection Saturating a Porous Medium in a Suspension of Gyrotactic Microorganisms

A linear stability analysis is performed to analyze bioconvection in a dilute suspension of gyrotactic microorganisms in horizontal shallow fluid layer cooling from below and saturated by a porous medium, in the rigid boundary case. It is established that due to cooling from below thermally stratified layer is stabilized, which opposes the development of bioconvection and the situations for oscillatory convection may take place. The stability criterion is obtained in terms of thermal Rayleigh number, bioconvection Rayleigh number, gyrotactic number, bioconvection Peclet number, measure of cell eccentricity, Prandtl number, and Lewis number. It is observed that the presence of porous medium results in decrease of the magnitude of critical bioconvection Rayleigh number in comparison with its non-existence; hence due to porous effect, the system becomes less stable.     

Linear Instability Analysis of a Suspension of Oxytactic Bacteria in Superimposed Fluid and Porous Layers

The purpose of this paper is to perform a pioneering investigation of the stability of bioconvection of oxytactic bacteria in superimposed fluid and porous layers. A dilute suspension of oxytactic bacteria in a shallow system that consists of superimposed fluid and porous layers is considered. A linear instability analysis of this problem is performed and the Galerkin method is utilized to solve the eigenvalue problem. The analysis leads to an equation for the critical Rayleigh number.*Author for correspondence: Tel.: +1-919-515-5292; Fax: +1-919-515-7968; e-mail: avkuznet@eos.ncsu.edu     

Thermo-Bioconvection in a Square Porous Cavity Filled by Oxytactic Microorganisms

This paper studies the thermo-bioconvection in a square porous cavity filled by oxytactic microorganisms. The Darcy model with Boussinesq approximation has been used to solve the flow and heat and mass transfer in the porous region. The governing equations formulated in terms of the dimensionless stream function, temperature and concentration have been solved using the finite difference method. Comparison with results from the open literature of the mean Nusselt number for a square cavity filled with a regular porous medium is made. It is shown that the results are in very good agreement. The main objective was to investigate the influence of the traditional Rayleigh number Ra = 10, 100, bioconvection Rayleigh number Rb = 10, 100, Lewis number Le = 1, 10, and Péclet number Pe = 0.1, 1 on the fluid flow and heat and mass transfer. Comprehensive analysis of an effect of these key parameters on the Nusselt and Sherwood numbers at the vertical walls has been conducted.     

Instability Analysis of Gyrotactic Microorganisms: A Combined Effect of High-Frequency Vertical Vibration and Porous Media

A linear stability analysis is carried out to predict the instability analysis in a dilute suspension of gyrotactic microorganisms in horizontal fluid-saturated porous layer influenced by high-frequency vertical vibration. The governing equations, describing the mean flow, are the time-averaged Boussinesq equations and the analytical solution of the problem has been obtained using Galerkin method. A secular relation involving bioconvection Rayleigh number and its vibrational analogs and other parameters have been established. The graphical interpretations for dependence of bioconvection Rayleigh–Darcy number and corresponding wave number, on gyrotactic number and bioconvection Péclet number in the presence of vibration are utilized to understand the problem.     

The Onset of Bioconvection in a Horizontal Porous-Medium Layer

In this note the problem of the onset of bioconvection in a horizontal layer occupied by a saturated porous medium is analyzed. Gyrotactic effects are incorporated in the analysis. The Darcy flow model is employed, and it is assumed that the bioconvection Péclet number is not greater than unity. Critical values of the bioconvection Rayleigh number and the corresponding critical Rayleigh number are obtained for various values of the bioconvection Péclet number, the gyrotaxis number and the cell eccentricity.     

Stability Analysis of Bioconvection of Gyrotactic Motile Microorganisms in a Fluid Saturated Porous Medium

Despite a large number of publications on bioconvection in suspensions of motile microorganisms, bioconvection in a fluid saturated porous medium is a relatively new area of research. This paper is motivated by experimental research by Kessler (1986) who established that a porous medium prevents the development of convection instability in algal suspensions. This suggests that there may exist a critical value of the permeability of a porous medium. If the permeability is smaller than critical, the system is stable and bioconvection does not develop. If the permeability is larger than critical, bioconvection may develop. This paper presents a model of bioconvection of gyrotactic motile microorganisms in a fluid saturated porous medium. The focus of this research is the determination of the critical value of permeability of a porous medium by a linear stability analysis. A simple but elegant analytical solution for the critical Darcy number is obtained.     

Thermal Instability in a Horizontal Porous Channel with Horizontal Through Flow and Symmetric Wall Heat Fluxes

The linear thermoconvective instability of the basic parallel flow in a plane and horizontal porous channel is investigated. The boundary walls are assumed to be impermeable and subject to symmetric and uniform heat fluxes. The wall heat fluxes produce either a net heating or a net cooling of the fluid saturated porous medium. A horizontal mass flow rate is externally impressed leading to a stationary basic state with a temperature gradient inclined to the vertical. A region of possibly unstable thermal stratification exists either in the lower half-channel (boundary heating), or in the upper half-channel (boundary cooling). The convective instability of the basic flow is governed by the Rayleigh number and by the Péclet number. In the case of boundary heating, the thermal instability arises when the Rayleigh number exceeds its critical value, that depends on the Péclet number. The change of the critical Rayleigh number as a function of the Péclet number is determined numerically for arbitrary normal modes oblique to the basic flow direction. The most dangerous modes are the longitudinal rolls, with a wave vector perpendicular to the basic velocity. There exists a minimum value of the Péclet number, 19.1971, below which no linear instability is detected.     

Investigation of the onset of bioconvection in a suspension of oxytactic microorganisms subjected to high-frequency vertical vibration

This paper investigates the effect of vertical vibration on the stability of a dilute suspension of oxytactic microorganisms in a shallow horizontal fluid layer. For the case of high-frequency vibration, an averaging method is utilized to derive the equations describing the mean flow by decomposing the solutions of governing equations into two components: one that varies slowly with time, and a second that varies rapidly with time. Linear stability analysis is used to investigate the stability of the obtained averaged equations. It is predicted that high-frequency, low-amplitude vertical vibration has a stabilizing effect on a suspension of oxytactic microorganisms confined in a shallow horizontal layer. PACS 47.27 Te; 68.60 Dv     

The Onset of Convection in a Suspension of Gyrotactic Microorganisms in Superimposed Fluid and Porous Layers: Effect of Vertical Throughflow

The purpose of this paper is to investigate the effect of vertical throughflow on the onset of bioconvection in a suspension of gyrotactic microorganisms. A dilute suspension of gyrotactic microorganisms in a shallow system that consists of superimposed fluid and porous layers is considered. A linear instability analysis of this problem is performed and the Galerkin method is utilized to solve the eigenvalue problem. The analysis leads to an equation for the critical Rayleigh number. It is shown that the vertical throughflow stabilizes the system.     

Anomalous reactivity of thermo-bioconvective nanofluid towards oxytactic microorganisms

The peristaltic flow of a non-Newtonian nanofluid with swimming oxytactic microorganisms through a space between two infinite coaxial conduits is investigated. A variable magnetic field is applied on the flow. The bioconvection flow and heat transfer in the porous annulus are formulated, and appropriate transformations are used, leading to the non-dimensionalized ruling partial differential equation model. The model is then solved by using the homotopy perturbation scheme. The effects of the germane parameters on the velocity profile, temperature distribution, concentration distribution, motile microorganism profile, oxytactic profile, pressure rise, and outer and inner tube friction forces for the blood clot and endoscopic effects are analyzed and presented graphically.It is noticed that the pressure rise and friction forces attain smaller values for the endoscopic model than for the blood clot model. The present analysis is believed to aid applications constituting hemodynamic structures playing indispensable roles inside the human body since some blood clotting disorders, e.g., haemophilia, occur when some blood constituents on the artery wall get confined away from the wall joining the circulation system.     

Penetrative convection in sublayer of water including density inversion

Numerical solutions of stability and convective flow in an infinite horizontal water layer, including density inversion, have been obtained using a finite element code. The evolution of the temperature field and flow pattern near the onset of convection are studied in detail. It is known that natural convection develops primarily in the lower unstably stratified layer. Of interest is the penetration of the convection rolls into the upper stably stratified layer and concurrent liquid entrainment as a function of the increasing Rayleigh number at different aspect ratios. Individual convection rolls may grow and expand before splitting up into two roll cells. It is shown that changing the aspect ratio influences critical Rayleigh number, flow symmetry, flow pattern, and transitions between flow patterns. Numerical results on heating from above or from below, agree well with available results in the literature. A correlation to predict critical Rayleigh numbers is given for the case of heating from above.     

Development and instability of steady convective flows in a square cavity heated from below and a field of vertically directed vibrational forces

The present paper is devoted to numerical investigation of the spatial structure and stability of secondary vibrational convective flows resulting from instability of the equilibrium of a fluid heated from below. Vibrations parallel to the vector of the gravitational force (vertical vibrations) are considered. As in earlier work [7–9], a region of finite size is used — a square cavity heated from below. It is shown that enhancement of the vibrational disturbance of the natural convective flow may either stabilize or destabilize flows with different spatial structures; it may also stabilize certain solutions of the system of convection equations that are unstable in the absence of vibrational forces. In addition, increase of the vibrational Rayleigh number can lead to a change of the mechanisms responsible for equilibrium instability and oscillatory instability of the secondary steady flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–18, March–April, 1991.I thank G. Z. Gershuni for assistance and extremely fruitful discussions of the results of the paper.     

Convective stability of a weakly conducting liquid in an electric field

In an inhomogeneously heated weakly conductive liquid (electrical conductivity 10^–12–1 cm^–1) located in a constant electric field a volume charge is induced because of thermal inhomogeneity of electrical conductivity and dielectric permittivity. The ponderomotive forces which develop set the liquid into intense motion [1–6]. However, under certain conditions equilibrium proves possible, and in that case the question of its stability may be considered. A theoretical analysis of liquid equilibrium stability in a planar horizontal condenser was performed in [2, 4]. Critical problem parameters were found for the case where Archimedean forces are absent [2]. Charge perturbation relaxation was considered instantaneous. It was shown that instability is of an oscillatory character. In [4] only heating from above was considered. Basic results were obtained in the limiting case of disappearingly small thermal diffusivity in the liquid (infinitely high Prandtl numbers). In the present study a more general formulation will be used to examine convective stability of equilibrium of a vertical liquid layer heated from above or below and located in an electric field. For the case of a layer with free thermally insulated boundaries, an exact solution is obtained. Values of critical Rayleigh number and neutral oscillation frequency for heating from above and below are found Neutral curves are constructed. It is demonstrated that with heating from below instability of both the oscillatory and monotonic types is possible, while with heating from above the instability has an oscillatory character. Values are found for the dimensionless field parameter at which the form of instability changes for heating from below and at which instability becomes possible for heating from above.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 16–23, September–October, 1976.In conclusion, the author thanks E. M. Zhukhovitskii for this interest in the study and valuable advice.     

Effect of a permeable partition on convective instability in a rectangular cavity

The equilibrium stability of a fluid, heated from below, in a rectangular cavity with a vertical permeable partition is investigated. The small perturbation problem is solved by the Galerkin-Kantorovich method. The relations obtained for the dependence of the critical Rayleigh numbers on the partition parameters and the cavity dimensions make it possible to identify regions in which either even or odd perturbations, sensitive to only the normal or only the tangential resistance of the partition, respectively, are responsible for equilibrium crisis. The effect of a permeable partition on the convective instability of a horizontal layer of fluid under various heating conditions was considered in [1–3], where a significant dependence of the critical Rayleigh numbers on the properties of the partition was established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 6–10, May–June, 1989.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Stability of a plane-parallel convective flow of a liquid in a horizontal layer

We consider the stationary plane-parallel convective flow, studied in [1], which appears in a two-dimensional horizontal layer of a liquid in the presence of a longitudinal temperature gradient. In the present paper we examine the stability of this flow relative to small perturbations. To solve the spectral amplitude problem and to determine the stability boundaries we apply a version of the Galerkin method, which was used earlier for studying the stability of convective flows in vertical and inclined layers in the presence of a transverse temperature difference or of internal heat sources (see [2]). A horizontal plane-parallel flow is found to be unstable relative to two critical modes of perturbations. For small Prandtl numbers the instability has a hydrodynamic character and is associated with the development of vortices on the boundary of counterflows. For moderate and for large Prandtl numbers the instability has a Rayleigh character and is due to a thermal stratification arising in the stationary flow.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 95–100, January–February, 1974.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

Experimental determination of the boundaries of the region of existence of anomalous convection flow in an tilted cube

This paper presents an experimentally study of the bifurcation of steady-state air convection in a cubic cavity heated from below under controlled deviations from equilibrium heating conditions due to a slow quasisteady-state tilt of the cavity at a predetermined angle α. It is found that in the supercritical range of Rayleigh numbers Ra at a tilt of the cavity not exceeding 7°, the existence of two stable steady-state convection regimes (normal and anomalous) with circulation in opposite directions is possible. A study is made of the transformations of the temperature distribution in the middle (with respect to the planes in which heat exchangers are located) plane during transition from the anomalous flow regime to the normal regime by instantaneous rotation of the entire mass of air in the cavity around the vertical axis by an angle of 90 to 135°. It is shown that this rotation occurs when the tilt of the cavity exceeds a critical value α _cr(Ra), which was determined experimentally for Rayleigh numbers 0 < Ra < 25Racr, where Racr is the critical Rayleigh number for stability of mechanical equilibrium for heating from below.     

Convective Instability in a Horizontal Porous Channel with Permeable and Conducting Side Boundaries

The stability analysis of the motionless state of a horizontal porous channel with rectangular cross-section and saturated by a fluid is developed. The heating from below is modelled by a uniform flux, while the top wall is assumed to be isothermal. The side boundaries are considered as permeable and perfectly conducting. The linear stability of the basic state is studied for the normal mode perturbations. The principle of exchange of stabilities is proved, so that only stationary normal modes need to be considered in the stability analysis. The eigenvalue problem for the neutral stability condition is solved analytically, and a closed-form dispersion relation is obtained for the neutral stability. The Darcy–Rayleigh number is expressed as an implicit function of the longitudinal wave number and of the aspect ratio. The critical wave number and the critical Darcy–Rayleigh number are evaluated for different aspect ratios. The preferred modes under critical conditions are detected. It is found that the selected patterns of instability at the critical Rayleigh number are two-dimensional, for slender or square cross-sections of the channel. On the other hand, instability is three dimensional when the critical width-to-height ratio, 1.350517, is exceeded. Eventually, the effects of a finite longitudinal length of the channel are discussed.