Coriolis Effect on the Stability of Centrifugally Driven Convection in a Rotating Anisotropic Porous Layer Subjected to Gravity

We investigate natural convection in a fluid saturated rotating anisotropic porous layer subjected to centrifugal gravitational and Coriolis body forces. The Darcy model (including the centrifugal, gravitational and Coriolis terms; and permeability anisotropy effects) and a modified energy equation (including the effects of thermal anisotropy) is used in the current analysis. The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection in the presence of thermal and mechanical anisotropy. It is shown that the preferred solution comprises roll cells aligned parallel to the vertical z-axis. As a result, it is found that the Coriolis acceleration (or Taylor number) and the gravitational term play no role in the stability of convection.     

Analytical Solution for Forced Convection in a Semi-Circular Channel Filled with a Porous Medium

Fully developed laminar forced convection inside a semi-circular channel filled with a Brinkman-Darcy porous medium is studied. Analytical solutions for flow and constant flux heat transfer are found using a mixture of Cartesian and cylindrical coordinates. The problem depends on a parameter s, which is proportional to the inverse square of the Darcy number. Velocity boundary layers exist when s is large. Both friction factor-Reynolds number product and Nusselt number are determined. Closed form expressions for the clear fluid () limit are found. Rare analytical solutions not only describe fundamental channel flows, but also serve as a check for more complicated numerical solutions.     

Stability of Solutal Convection in a Gravity Modulated Mushy Layer During the Solidification of Binary Alloys

The linear stability theory is used to investigate analytically the effects of gravity modulation on solutal convection in the mushy layer of solidifying binary alloys. The gravitational field consists of a constant part and a sinusoidally varying part, which is synonymous to a vertically oscillating mushy layer subjected to constant gravity. The linear stability results are presented for both the synchronous and subharmonic solutions. It is demonstrated that up to the transition point between the synchronous and subharmonic regions, increasing the frequency of vibration rapidly stabilizes the solutal convection. Beyond the transition point, further increases in the frequency tend to destabilize the solutal convection, but gradually. It is also demonstrated that the effect of increasing the ratio of the Stefan number and the solid composition (₀) is to destabilize the solutal convection.     

Destabilising A Fluid Saturated Gravity Modulated Porous Layer Heated From Above

The linear stability theory is used to investigate analytically the possibility of the motionless basic state in a porous layer heated from above and subjected to vibration. The linear stability results presented for the specific case of low amplitude vibration shows that there exists a bandwidth of frequencies for which the convection in a porous layer with a stable density gradient can be destabilized. In addition the scaled Darcy–Prandtl number is shown to influence the onset of the subharmonic and synchronous solutions.     

Nonlinear Stability of Convection in a Porous Layer with Solid Partitions

We show that for many classes of convection problem involving a porous layer, or layers, interleaved with finite but non-deformable solid layers, the global nonlinear stability threshold is exactly the same as the linear instability one. The layer(s) of porous material may be of Darcy type, Brinkman type, possess an anisotropic permeability, or even be such that they are of local thermal non-equilibrium type where the fluid and solid matrix constituting the porous material may have different temperatures. The key to the global stability result lies in proving the linear operator attached to the convection problem is a symmetric operator while the nonlinear terms must satisfy appropriate conditions.     

Finite element solution of the equations governing the flow of electrolyte in charged microporous membranes

Electrical double-layer effects are unimportant in flows through porous media except when the Debye length k^?1 is comparable in magnitude with the pore radius a. Under these conditions the equations governing the flow of electrolyte are those of Stokes, Nernst-Planck and Poisson. These equations are non-linear and require numerical solution. The finite element method provides a useful basis for solution and various algorithms are investigated. The numerical stability and errors of each scheme are analysed together with the development of an appropriate finite element mesh. The electro-osmotic flow of a typical electrolyte (barium chloride) through a uniformly charged cylindrical membrane pore is investigated and the ion fluxes are post-computed from the numerical solutions. The ion flux is shown to be strongly dependent on both zeta potential and pore radius, ka, indicating the effects of overlapping electrical double layers.     

Linear Stability and Convection in a Gravity Modulated Porous Layer Heated from Below: Transition from Synchronous to Subharmonic Solutions

The linear stability theory is used to investigate analytically the effects of gravity modulation on convection in a homogenous porous layer heated from below. The linear stability results are presented for both the synchronous and subharmonic solutions and the exact point for the transition from synchronous to subharmonic solutions is computed. It is also demonstrated that increasing the excitation frequency rapidly stabilizes the convection up to the transition point from synchronous to subharmonic convection. Beyond the transition point, the effect of increasing the frequency is to slowly destabilize the convection.     

The Onset of Darcy–Brinkman Electroconvection in a Dielectric Fluid Saturated Porous Layer

The combined effect of a vertical AC electric field and the boundaries on the onset of Darcy–Brinkman convection in a dielectric fluid saturated porous layer heated either from below or above is investigated using linear stability theory. The isothermal bounding surfaces of the porous layer are considered to be either rigid or free. It is established that the principle of exchange of stability is valid irrespective of the nature of velocity boundary conditions. The eigenvalue problem is solved exactly for free–free (F/F) boundaries and numerically using the Galerkin technique for rigid–rigid (R/R) and lower-rigid and upper-free (F/R) boundaries. It is observed that all the boundaries exhibit qualitatively similar results. The presence of electric field is emphasized on the stability of the system and it is shown that increasing the AC electric Rayleigh number R _ea is to facilitate the transfer of heat more effectively and to hasten the onset of Darcy–Brinkman convection. Whereas, increase in the ratio of viscosities Λ and the inverse Darcy number Da ⁻¹ is to delay the onset of Darcy–Brinkman electroconvection. Besides, increasing R _ea and Da ⁻¹ as well as decreasing Λ are to reduce the size of convection cells.     

Onset of convection in a horizontal porous channel with uniform heat generation using a thermal nonequilibrium model

This paper considers the onset of free convection in a horizontal fluid-saturated porous layer with uniform heat generation. Attention is focused on cases where the fluid and solid phases are not in local thermal equilibrium, and where two energy equations describe the evolution of the temperature of each phase. Standard linearized stability theory is used to determine how the criterion for the onset of convection varies with the inter-phase heat transfer coefficient, H, and the porosity-modified thermal conductivity ratio, γ. We also present asymptotic solutions for small values of H. Excellent agreement is obtained between the asymptotic and the numerical results.     

The effect of mechanical and thermal anisotropy on the stability of gravity driven convection in rotating porous media in the presence of thermal non-equilibrium

We investigate Rayleigh–Benard convection in a porous layer subjected to gravitational and Coriolis body forces, when the fluid and solid phases are not in local thermodynamic equilibrium. The Darcy model (extended to include Coriolis effects and anisotropic permeability) is used to describe the flow, whilst the two-equation model is used for the energy equation (for the solid and fluid phases separately). The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection and the effect of both thermal and mechanical anisotropy on the critical Rayleigh number is discussed.     

Fluid convection in a rotating porous layer under modulated temperature on the boundaries

The linear stability of thermal convection in a rotating horizontal layer of fluid-saturated porous medium, confined between two rigid boundaries, is studied for temperature modulation, using Brinkman’s model. In addition to a steady temperature difference between the walls of the porous layer, a time-dependent periodic perturbation is applied to the wall temperatures. Only infinitesimal disturbances are considered. The combined effect of rotation, permeability and modulation of walls’ temperature on the stability of flow through porous medium has been investigated using Galerkin method and Floquet theory. The critical Rayleigh number is calculated as function of amplitude and frequency of modulation, Taylor number, porous parameter and Prandtl number. It is found that both, rotation and permeability are having stabilizing influence on the onset of thermal instability. Further it is also found that it is possible to advance or delay the onset of convection by proper tuning of the frequency of modulation of the walls’ temperature.     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Stability of Thermosolutal Natural Convection in Superposed Fluid and Porous Layers

This article deals with the onset of thermosolutal natural convection in horizontal superposed fluid and porous layers. A linear stability analysis is performed using the one-domain approach. As in the thermal convection case, the results show a bimodal nature of the marginal stability curves where each mode corresponds to a different convective instability. At small wave numbers, the convective flow occurs in the whole cavity (“porous mode”) while perturbations of large wave numbers lead to a convective flow mainly confined in the fluid layer (“fluid mode”). Furthermore, it is shown that the onset of thermosolutal natural convection is characterized by a multi-cellular flow in the fluid region for negative thermal Rayleigh numbers. For positive thermal Rayleigh numbers, the convective flow takes place both in the fluid and porous regions. The influence of the depth ratio and thermal diffusivity ratio is also investigated for a wide range of the thermal Rayleigh numbers.     

Stability of Moderate Vadasz Number Solutal Convection in a Cylindrical Mushy Layer Subjected to Vertical Vibration

We consider vibration effects on the stability of solutal convection in a mushy layer being cast in a cylindrical geometry. The near eutectic limit is applied and moderate Vadasz numbers are considered to retain the second-order time derivative in the Darcy equation. Since small to moderate radii casting crucibles are the current area of interest, only synchronous modes are analyzed. The results indicate that the presence of vibration in solidifying mushy layers stabilizes the convection, and provides a quantification of the Rayleigh number associated with solutal convection. Of particular interest is the fact that in solidifying systems, the Rayleigh numbers are significantly smaller than that of a passive porous layer.     

Stability of a fluid-saturated porous medium contained in a vertical cylinder heated from below by forced convection

A linear stability analysis determining the critical Rayleigh number R _c for onset of convection in a bounded vertical cylinder containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top surface, and bottom surface heated by forced convection. This Newtonian heating of the bottom surface involves a Biot number Bi that allows consideration of the continuum of boundary conditions ranging from constant heat flux, with global minimum R _min=27.096 found as Bi→0, to isothermal, with global minimum R _min=4π² found as Bi→ ∞. In both cases and for most cylinder aspect ratios, incipient convection sets in as an asymmetric mode, though islands of aspect ratio exist where the onset mode is symmetric. Sample three-dimensional renderings of disturbance temperature distributions showing preferred modes at onset of convection for fixed Bi are provided and an analytical fit to R _min as a function of Bi is given.     

Stability of double-diffusive convection induced by selective absorption of radiation in a fluid layer

Linear and nonlinear stability analyses were performed on a fluid layer with a concentration-based internal heat source. Clear bimodal behaviour in the neutral curve (with stationary and oscillatory modes) is observed in the region of the onset of oscillatory convection, which is a previously unobserved phenomenon in radiation-induced convection. The numerical results for the linear instability analysis suggest a critical value γ _c of γ, a measure for the strength of the internal heat source, for which oscillatory convection is inhibited when γ > γ _c . Linear instability analyses on the effect of varying the ratio of the salt concentrations at the upper and lower boundaries conclude that the ratio has a significant effect on the stability boundary. A nonlinear analysis using an energy approach confirms that the linear theory describes the stability boundary most accurately when γ is such that the linear theory predicts the onset of mostly stationary convection. Nevertheless, the agreement between the linear and nonlinear stability thresholds deteriorates for larger values of the solute Rayleigh number for any value of γ.     

Coriolis effect on thermal convective instability of viscoelastic fluids in a rotating porous cylindrical annulus

Linear stability analysis of thermal convection is studied for a viscoelastic fluid in a rotating porous cylindrical annulus. The modified Darcy–Jeffrey model with the addition of the Coriolis term in a rotating frame of reference is applied to characterize the non-Newtonian rheology in porous media. We investigate how the interaction among the Coriolis force, viscoelasticity, and bounded sidewalls affects the preferred mode at the onset of convection. The results show that for a slowly rotating case, the oscillatory mode is always preferred for any considered cylindrical radii. However, for a moderately rotating case, the oscillatory preferred mode only arises intermittently as the outer cylindrical radius gradually increases. This result is quite different from the case for viscoelastic fluids in a rotating porous layer or in a porous cylinder without rotation. Further, we discover that for a pair of given cylindrical radii when the Taylor number exceeds a critical value depending on the viscoelastic parameters, the oscillatory convection does not occur. We also examine how the variations of the Taylor number and the viscoelastic parameters affect the patterns of temperature disturbance at the onset of convection.     

Investigation of convection in a horizontal cylindrical layer of a saturated porous medium

The structures of the convective motions and the nature of the heat transfer in a horizontal cylindrical layer are studied numerically for the Forchheimer model of a porous medium in the Boussinesq approximation. New asymmetric solutions of the equations of convection flow through a porous medium are found. Their development, domains of existence, and stability are investigated. One consists of a multivortex structure with asymmetric vortices in the near-polar region. Another asymmetric solution is realized at large Grashof numbers in the form of a convective plume deflected from the vertical. The threshold Grashof number of formation of the asymmetric motions depends on the Prandtl number and the cylindrical layer thickness.     

Variable permeability effect on convection in binary mixtures saturating a porous layer

The Darcy Model with the Boussinesq approximation is used to study natural convection in a shallow porous layer, with variable permeability, filled with a binary fluid. The permeability of the medium is assumed to vary exponentially with the depth of the layer. The two horizontal walls of the cavity are subject to constant fluxes of heat and solute while the two vertical ones are impermeable and adiabatic. The governing parameters for the problem are the thermal Rayleigh number, R _T, the Lewis number, Le, the buoyancy ratio, φ, the aspect ratio of the cavity, A, the normalized porosity, ε, the variable permeability constant, c, and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in an infinite layer, an analytical solution of the steady form of the governing equations is obtained on the basis of the parallel flow approximation. The onset of supercritical convection, or subcritical, convection are predicted by the present theory. A linear stability analysis of the parallel flow model is conducted and the critical Rayleigh number for the onset of Hopf’s bifurcation is predicted numerically. Numerical solutions of the full governing equations are found to be in excellent agreement with the analytical predictions.     

On the Effect of Anisotropy on the Stability of Convection in Rotating Porous Media

We investigate natural convection in an anisotropic porous layer subjected to centrifugal body forces. The Darcy model (including centrifugal and permeability anisotropy effects) is used to describe the flow and a modified energy equation (including the effects of thermal anisotropy) is used in the current analysis. The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection in the presence of thermal and mechanical anisotropy. It is found that the convection is stabilized when the thermal anisotropy ratio (which is a function of the thermal and mechanical anisotropy parameters) is increased in magnitude.