Stability of buoyancy-driven convection in an Oldroyd-B fluid-saturated anisotropic porous layer

The nonlinear stability of thermal convection in a layer of an Oldroyd-B fluid-saturated Darcy porous medium with anisotropic permeability and thermal diffusivity is investigated with the perturbation method. A modified Darcy-Oldroyd model is used to describe the flow in a layer of an anisotropic porous medium. The results of the linear instability theory are delineated. The thresholds for the stationary and oscillatory convection boundaries are established, and the crossover boundary between them is demarcated by identifying a codimension-two point in the viscoelastic parameter plane. The stability of the stationary and oscillatory bifurcating solutions is analyzed by deriving the cubic Landau equations. It shows that these solutions always bifurcate supercritically. The heat transfer is estimated in terms of the Nusselt number for the stationary and oscillatory modes. The result shows that, when the ratio of the thermal to mechanical anisotropy parameters increases, the heat transfer decreases.     

An extended theory to predict the onset of viscous instabilities for miscible displacements in porous media

Many enhanced oil recovery schemes involve the displacement of oil by a miscible fluid. Whether a displacement is stable or unstable has a profound effect on how efficiently a solvent displaces oil within a reservoir. That is, if viscous fingers are present, the displacement efficiency and, hence, the economic return of the recovery scheme is seriously impaired bacause of macroscopic bypassing of the oil. As a consequence, it is of interest to be able to predict the boundary which separates stable displacements from those which are unstable.This paper presents a dimensionless scaling group for predicting the onset of hydrodynamic instability of a miscible displacement in porous media. An existing linear perturbation analysis was extended in order to obtain the scaling group. The new scaling group differs from those obtained in previous studies because it takes into account a variable unperturbed concentration profile, both transverse dimensions of the porous medium, and both the longitudinal and the transverse dispersion coefficient.It has been shown that stability criteria derived in the literature are special cases of the general condition given here. Therefore, the stability criterion obtained in this study should be used for a displacement conducted under arbitrary conditions. The stability criterion is verified by comparing it with miscible displacement experiments carried out in a Hele-Shaw cell. Moreover, a comparison of the theory with some porous medium experiments from the literature also supports the validity of the theory.Nomenclature c solvent concentration - C _g fractional glycerine volume - D molecular diffusion coefficient, cm²/s - D _L longitudinal dispersion coefficient, cm²/s - D _T transverse dispersion coefficient, cm²/s - g gravitational acceleration, cm/s² - h distance between the plates, cm - I _sr dimensionless scaling group - k permeability, cm² - L _x width of the porous medium, cm - L _y height of the porous medium, cm - t time, s - u velocity in thex direction, cm/s - v velocity in they direction, cm/s - V displacement velocity, cm/s - w velocity in thez direction, cm/s - z length of the graded viscosity bank, cm - eigenvalue in thex direction - eigenvalue in they direction - wave number - viscosity, poise - density, g/cc - time constant, s^-1 - porosity     

多孔介质中二维对流传热的分叉

本文利用分叉理论研究了流体饱和的二维多孔介质从底部加热所引起的自然对流,用有限差分方法确定对流的分叉进程;揭示其模式转换机理及分叉对非正常流动图象形成的影响;同时确定了矩形截面宽高比与临界端利数的关系。还提出了一个判别分支稳定笥的简明方法。     

Vibration effect on convection onset in a system consisting of a horizontal pure liquid layer and a layer of liquid-saturated porous medium

The onset of convection in a system of two horizontal layers (a pure liquid and a porous medium saturated with the same liquid) heated from below under the action of vertical vibration is investigated. For describing the free thermal convection, in the liquid layer the Boussinesq approximation and in the porous layer the Darcy-Boussinesq approximation are used. In the limiting case of a thin liquid layer, effective boundary conditions on the upper boundary of the porous layer with account for convection in the liquid layer are obtained and it is shown that vibration has a stabilizing effect, whereas the presence of a liquid layer leads to destabilization. For an arbitrary liquid to porous layer thickness ratio the onset of convection is investigated numerically. In the case of a thin liquid layer there are two (short-and long-wave) unstable modes. In the case of thick layers the neutral curves are unimodal. Vibration has a stabilizing effect on perturbations with any wave number but affects short-wave perturbations much more strongly than long-wave ones.     

Effects of anisotropy on the free convection from a vertical plate embedded in a porous medium

The effects of anisotropy on the steady laminar boundary-layer free convection over a vertical impermeable surface are analysed by using the method of integral relations. If the permeability in the direction orthogonal to the plate is greater than the permeability along the plate, then there is an increase in the temperature field.     

The inertial effect on the natural convection flow within a fluid-saturated porous medium

The general momentum equation for fluid flow within a porous medium is supposedly valid for any fluid-porous medium configuration. One of the main concerns of using the general equations refers to the inclusion of both inertia terms, namely, the convective inertia term and the Forchheimer term. In this study, we go beyond the important discussion about the correctness of including both terms in the general momentum equations by focusing upon the effect of the convective inertia term on the heat transfer results. The fluid-porous medium system considered here is a cavity bounded by solid surfaces with vertical walls maintained at constant but different temperatures. The natural convection problem is solved numerically, and the results are compared with a general theory developed by using the method of scale analysis. It is demonstrated that the convective inertia term effect upon the heat transfer results is minor for 0.01 ≤ Pr ≤ 1, 10 ≤ Ra_D ≤ 10⁴, 10⁻⁸ ≤ Da ≤ 10⁻², and porosities 0.4 and 0.8. It is also shown that, contrary to the general belief, the convective inertial effect upon the heat transfer within the cavity is minimized when the Prandtl number is reduced.     

Penetrative convection in a horizontally isotropic porous layer

Penetrative convection in a horizontally isotropic porous layer is investigated primarily using an internal heat sink model and alternatively, a quadratic density temperature law. Employing the heat sink model, we show that the temporal growth rate for the linearised system is real, which allows us to perform a linear instability analysis. A nonlinear energy analysis is also presented, yielding a global stability threshold. We compare the heat sink model to the quadratic density model and find that the linearised systems are adjoint, implying that the instability boundaries which can be derived from the two models are the same. Received April 12, 2002 / Published online September 4, 2002 RID="a" ID="a" e-mail: Magdalen.Carr@durham.ac.uk RID="b" ID="b" e-mail: s.d.putter@tue.nl Communicated by Brian Straughan, Durham     

Nonlinear convection in a non-Darcy porous medium

In this paper, the natural convection in a non-Darcy porous medium is studied using a temperature-concentration-dependent density relation. The effect of the two parameters responsible for the nonlinear convection is analyzed for different values of the inertial parameter, dispersion parameters, Rayleigh number, Lewis number, Soret number, and Dufour number. In the aiding buoyancy, the tangential velocity increases steeply with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter when the inertial effect is zero. However, when the inertial effect is non-zero, the effect of the nonlinear temperature parameter and the nonlinear concentration parameter on the tangential velocity is marginal. The concentration distribution varies appreciably and spreads in different ranges for different values of the double dispersion parameters, the inertial effect parameter, and also for the parameters which control the nonlinear temperature and the nonlinear concentration. Heat and mass transfer varies extensively with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter depending on Dacry and non-Darcy porous media. The variation in heat and mass transfer when all the effects, i.e., the inertial effect, double dispersion ef- fects, and Soret and Dufour effects, are simultaneously zero and non-zero. The combined effects of the nonlinear temperature parameter, the nonlinear concentration parameter and buoyancy are analyzed. The effect of the nonlinear temperature parameter and the nonlinear concentration parameter and also the cross diffusion effects on heat and mass transfer are observed to be more in Darcy porous media compared with those in non- Darcy porous media. In the opposing buoyancy, the effect of the temperature parameter is to increase the heat and mass transfer rate, whereas that of the concentration parameter is to decrease.     

The onset of penetrative double-diffusive convection

The onset of double-diffusive convection in a horizontal fluid layer is studied. The density is assumed to depend quadratically on the temperature and linearly on the solute concentration. Under the Boussinesq approximation, the linear stability of the conduction state is investigated with respect to the oscillatory and steady convection modes. For steady onset, the critical thermal Rayleigh number is found to be a double-valued function of the solutal Rayleigh number as long as the relative maximum of the density profile exists within the fluid layer. Driving mechanisms of the steady convections are discussed.     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL（Karhunen-Loeve展开）-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

Radiation effect on mixed convection from a horizontal surface in a porous medium

The effect of thermal radiation on the non-Darcy mixed convection flow over a non-isothermal horizontal surface immersed in a saturated porous medium has been studied. The wall temperature is assumed to have a power-law variation with the distance measured from the leading edge of the plate. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using a finite-difference scheme. For some particular cases, the self-similar solution has also been obtained. The heat transfer is found to be strongly influenced by the radiative flux number, buoyancy parameter, variation of wall temperature, non-Darcy parameter and the nature of the free stream velocity.     

Horizontal thermal convection in a porous medium

Linearised instability and nonlinear stability bounds for thermal convection in a fluid-filled porous finite box are derived. A nonuniform temperature field in the steady state is generated by maintaining the vertical walls at different temperatures. The linearised instability threshold is shown to be well above the global stability boundary, which is strongly suggested by the lack of symmetry in the perturbed system. The numerical results are evaluated utilising a newly developed Legendre-polynomial-based spectral method.     

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.     

Linear stability of triple-diffusive convection in micropolar ferromagnetic fluid saturating porous medium

The triple-diffusive convection in a micropolar ferromagnetic fluid layer heated and soluted from below is considered in the presence of a transverse uniform magnetic field. An exact solution is obtained for a flat fluid layer contained between two free boundaries. A linear stability analysis and a normal mode analysis method are carried out to study the onset convection. For stationary convection, various parameters such as the medium permeability, the solute gradients, the non-buoyancy magnetization, and the micropolar parameters (i.e., the coupling parameter, the spin diffusion parameter, and the micropolar heat conduction parameter) are analyzed. The critical magnetic thermal Rayleigh number for the onset of instability is determined numerically for a sufficiently large value of the buoyancy magnetization parameter M ₁. The principle of exchange of stabilities is found to be true for the micropolar fluid heated from below in the absence of the micropolar viscous effect, the microinertia, and the solute gradients. The micropolar viscous effect, the microinertia, and the solute gradient introduce oscillatory modes, which are non-existent in their absence. Sufficient conditions for the non-existence of overstability are also obtained.     

Natural convection flow of a viscous incompressible fluid in a rectangular porous cavity heated from below with cold sidewalls

We consider the flow, which is induced by differential heating on the boundaries of a porous cavity heated from below. In particular we allow the sidewalls to have the same cold temperature as the upper surface, and thus the problem is a variant of the Darcy-Bénard convection problem, but one where there is flow at all non-zero Grashof numbers. Attention is focused on how the flow and heat transfer is affected by variations in the cavity aspect ratio, the Grashof number and the Darcy number. The flow becomes weaker as the Darcy number decreases from the pure fluid limit towards the Darcy-flow limit. In addition the number of cells which form in the cavity varies primarily with the aspect ratio and is always even due to the symmetry imposed by the cold sidewalls.     

TEMPERATURE PROFILES OF LOCAL THERMAL NONEQUILIBRIUM FOR THERMAL DEVELOPING FORCED CONVECTION IN POROUS MEDIUM PARALLEL PLATE CHANNEL

Based on the two-energy equation model, taking into account viscous dissipation due to the interaction between solid skeleton and pore fluid flow, temperature expressions of the solid skeleton and pore fluid flow are obtained analytically for the thermally developing forced convection in a saturated porous medium parallel plate channel, with walls being at constant temperature. It is proved that the temperatures of the two phases for the local thermal nonequilibrium approach to the temperature derived from the one-energy equation model for the local thermal equilibrium when the heat exchange coefficient goes to infinite. The temperature profiles are shown in figures for different dimensionless parameters and the effects of the parameters on the local thermal nonequilibrium are revealed by parameter study.     

Diffuse approximation method for solving natural convection in porous media

The diffuse approximation is presented and applied to natural convection problems in porous media. A comparison with the control volume-based finite-element method shows that, overall, the diffuse approximation appears to be fairly attractive.Nomenclature H height of the cavities - I functional - K permeability - p(M _i ,M) line vector of monomials - p ^T p-transpose - M current point - Nu Nusselt number - Ri inner radius - Ro outer radius - Ra Rayleigh number - x, y cartesian coordinates - u, v velocity components - T temperature - M vector of estimated derivatives - _t thermal diffusivity - coefficient of thermal expansion - practical aperture of the weighting function - scalar field - (M, M _i ) weighting function - streamfunction - kinematic viscosity     

The Prandtl number effect near the onset of Bénard convection in a porous medium

This note focuses on Kladias and Prasad's claim that the critical Rayleigh number for the onset of Bénard convection in an infinite horizontal porous layer increases as the Prandtl number decreases, and argues that the critical Rayleigh number (Ra_c) depends only on the Darcy number (Da), as linear stability analysis indicates. The two-dimensional steady-convection problem is then solved numerically to document the convection heat transfer effect of the Rayleigh number, Darcy number, Prandtl number, and porosity. The note concludes with an empirical correlation for the overall Nusselt number, which shows the effect of Prandtl number at above-critical Rayleigh numbers. The correlation is consistent with the corresponding correlation known for Bénard convection in a pure fluid.     

Measurements of thermally-induced convection in a model porous medium

The velocity field generated by thermal convection in a model porous medium is experimentally determined by means of both PIV and LDA techniques. Details of matching refraction index under non isothermal conditions are given. Fields are measured in the empty parallelepipedic cell and in a model medium made of parallel circular bundles. Results are in good agreement. Moreover, by an averaging technique, we are able to measure seeping velocity profiles.     

Convection due to the selective absorption of radiation in a porous medium

A model for convection due to the selective absorption of radiation in a fluid saturated porous medium is investigated. The model is based on a similar one introduced for a viscous fluid by Krishnamurti [x]. Employing this adapted model we show the growth rate for the linearised system is real. A linear instability analysis is performed. Global stability thresholds are also found using nonlinear energy theory. An excellent agreement is found between the linear instability and nonlinear stability Rayleigh numbers, so that the region of potential subcritical instabilities is very small, demonstrating that the linear theory accurately emulates the physics of the onset of convection. Received February 10, 2003 / Accepted February 10, 2003/ Published online May 9, 2003 / B. Straughan