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.     

Nonlinear cellular convection and heat transport in a porous medium

Nonlinear steady cellular convection in a fluid-saturated porous medium is investigated using the technique of spectral analysis. The effect of permeability is shown to contract the cell and to damp the convection process. The influence of Prandtl number, though small, is seen only in the fourth order term. The cross-interactions of the higher modes caused by nonlinear effects are considered through the modal Rayleigh number R . The possibility of the existence of a steady solution with two self-excited modes in certain regions is predicted. A detailed discussion of the heat transport is made. The theoretical value of the Nusselt number is found to be in good agreement with the experimental results. The similarities and qualitative differences between the present analysis and that of the power integral technique are brought out.     

A numerical method for forced convection in porous and homogeneous fluid domains coupled at interface by stress jump

A numerical method was developed for flows involving an interface between a homogeneous fluid and a porous medium. It is based on the finite volume method with body‐fitted and multi‐block grids. The Brinkman–Forcheimmer extended model was used to govern the flow in the porous medium region. At its interface, the flow boundary condition imposed is a shear stress jump, which includes the inertial effect, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The forced convection through a porous insert over a backward‐facing step is investigated. The results are presented with flow configurations for different Darcy numbers, 10^?2 to 10^?5, porosity from 0.2 to 0.8, Reynolds number from 10 to 800, and the ratio of insert length to channel height from 0.1 to 0.3. The heat transfer is improved by using porous insert. To enhance the heat transfer with minimal frictional losses, it is preferable to have a medium length of insert with medium Darcy number, and larger Reynolds number. The interfacial stress jump coefficients β and β₁ were varied from ?1 to 1, and within this range the average and local lower‐wall Nusselt numbers are not sensitive to the parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

Onset of Convection in a Porous Layer with Continuous Periodic Horizontal Stratification. Part I. Two-Dimensional Convection

The onset of convection in a horizontal porous layer is investigated theoretically. The permeability of the porous medium is a continuous periodic function of the horizontal x coordinate. Floquet theory has been employed to determine the favoured two-dimensional mode of convection. For a wide range of periods of the permeability variation, a matrix eigenvalue technique with eighth order accuracy has been employed to find the critical Darcy– Rayleigh number. This is supplemented by a multiple-scales analysis of the large-period limit, and a brief consideration of the anisotropic limit for very short periods.     

Numerical simulation of the flow around a porous covering square cylinder in a channel via lattice Boltzmann method

In this paper, a detailed investigation on the flow past a porous covering cylinder is presented through the lattice Boltzmann method. The Brinkman‐Forchheimer‐extended Darcy model is adopted for the entire flow field with the solid, fluid, and porous medium. The effects of several parameters, such as porous layer thickness, Darcy number, porosity, and Reynolds number on flow field are discussed. Compared with the case of a solid cylinder, the present work shows that the porous layer may play an important role on the flow, the lift and drag force exerted on the cylinder. The numerical results indicate that the maximal drag coefficient C_d and maximal amplitude of lift coefficient C_l exist at certain Darcy number which is in the range of 10^?6–10^?2. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite Element Modelling of Flow Through a Porous Medium Between Two Parallel Plates Using The Brinkman Equation

Despite the widespread use of the Darcy equation to model porous flow, it is well known that this equation is inconsistent with commonly prescribed no slip conditions at flow domain walls or interfaces between different sections. Therefore, in cases where the wall effects on the flow regime are expected to be significant, the Darcy equation which is only consistent with perfect slip at solid boundaries, cannot predict velocity and pressure profiles properly and alternative models such as the Brinkman equation need to be considered. This paper is devoted to the study of the flow of a Newtonian fluid in a porous medium between two impermeable parallel walls at different Darcy parameters (Da). The flow regime is considered to be isothermal and steady. Three different flow regimes can be considered using the Brinkman equation: free flow (Da > 1), porous flow (high permeability, 1 > Da > 10⁻⁶) and porous flow (low permeability Da < 10⁻⁶). In the present work the described bench mark problem is used to study the effects of solid walls for a range of low to high Darcy parameters. Both no-slip and slip conditions are considered and the results of these two cases are compared. The range of the applicability of the Brinkman equation and simulated results for different cases are shown.     

Unsteady natural convection in a liquid-saturated porous enclosure with local thermal non-equilibrium effect

Stability analysis of free convection in a liquid-saturated sparsely-packed porous medium with local-thermal-non-equilibrium (LTNE) effect is presented. For the vertical boundaries free–free, adiabatic and rigid–rigid, adiabatic are considered while for horizontal boundaries it is the stress-free, isothermal and rigid–rigid, isothermal boundary combinations we consider. From the linear theory, it is apparent that there is advanced onset of convection in a shallow enclosure followed by that in square and tall enclosures. Asymptotic analysis of the thermal Rayleigh number for small and large values of the inter-phase heat transfer coefficient is reported. Results of Darcy–Bénard convection (DBC) and Rayleigh–Bénard convection can be obtained as limiting cases of the study. LTNE effect is prominent in the case of Brinkman–Bénard convection compared to that in DBC. Using a multi-scale method and by performing a non-linear stability analysis the Ginzburg–Landau equation is derived from the five-mode Lorenz modal. Heat transport is estimated at the lower plate of the channel. The effect of the Brinkman number, the porous parameter and the inter-phase heat transfer coefficient is to favour delayed onset of convection and thereby enhanced heat transport while the porosity-modified ratio of thermal conductivities shows the opposite effect.

     

Instability of the Benard–Marangoni Convection in a Porous Layer Affected by a Non-Vertical Magnetic Field

The onset of the Benard–Marangoni convection in a horizontal porous layer permeated by a magnetohydrodynamic fluid with a nonlinear magnetic permeability is examined. The porous layer is assumed to be governed by the Brinkman model; it is bounded by a rigid surface from below and by a non-deformable free surface from above and subjected to a non-vertical magnetic field. The critical effective Marangoni number and the critical Rayleigh number are obtained for different values of the effective Darcy number, Biot number, Chandrasekhar number, nonlinear magnetic parameter, and angle from the vertical axis for the cases of stationary convection and overstability. The related eigenvalue problem is solved by using the first-order Chebyshev polynomial method.     

Effects of Brownian Motion and Thermophoresis on the Mixed Convection of Nanofluid in a Porous Channel Including Flow Reversal

This article is devoted to combined convection heat transfer of nanofluids through a vertical channel filled with a homogeneous and isotropic porous medium. The flow is assumed to be fully developed and the “Brinkman extended Darcy” model is used for the flow in the porous media and “clear compatible” viscous dissipation model is considered. Also the model utilized for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The governing momentum, energy, and nanopartices volume fraction equations are solved both analytically and numerically. The effects of the influential dimensionless parameters such as Brownian and thermophoresis parameters, mixed convection parameter (Gr/Re), Brinkman, Darcy and Lewis numbers on dimensionless velocity and temperature distributions and pressure drop are studied. Also, the results of the Nusselt number for the both left and right walls are presented and discussed.     

On the Chebyshev collocation spectral approach to stability of fluid flow in a porous medium

In this paper, the temporal development of small disturbances in a pressure‐driven fluid flow through a channel filled with a saturated porous medium is investigated. The Brinkman flow model is employed in order to obtain the basic flow velocity distribution. Under normal mode assumption, the linearized governing equations for disturbances yield a fourth‐order eigenvalue problem, which reduces to the well‐known Orr–Sommerfeld equation in some limiting cases solved numerically by a spectral collocation technique with expansions in Chebyshev polynomials. The critical Reynolds number Re_c, the critical wave number α_c, and the critical wave speed c_c are obtained for a wide range of the porous medium shape factor parameter S. It is found that a decrease in porous medium permeability has a stabilizing effect on the fluid flow. Copyright © 2008 John Wiley & Sons, Ltd.     

The Effect of Magnetic Field Dependent Viscosity on Thermosolutal Convection in a Ferromagnetic Fluid Saturating a Porous Medium

The effect of magnetic field dependent viscosity on thermosolutal convection in a ferromagnetic fluid saturating a porous medium is considered for a fluid layer heated and soluted from below in the presence of uniform magnetic field. Using linearized stability theory and normal mode analysis, an exact solution is obtained for the case of two free boundaries. For case of stationary convection, medium permeability has a destabilizing effect, whereas a stable solute gradient and magnetic field dependent viscosity have a stabilizing effect on the system. In the absence of magnetic field dependent viscosity, the destabilizing effect of non-buoyancy magnetization is depicted but in the presence of magnetic field dependent viscosity non-buoyancy magnetization may have a destabilizing or stabilizing effect on the onset of instability. The critical wave number and the critical magnetic thermal Rayleigh number for the onset of instability are also determined numerically for sufficiently large values of buoyancy magnetization parameter M₁ and the results are depicted graphically. The principle of exchange of stabilities is found to hold true for the ferromagnetic fluid saturating a porous medium heated from below in the absence of stable solute gradient. The oscillatory modes are introduced due to the presence of the stable solute gradient, which were non-existent in its absence. A sufficient condition for the non-existence of overstability is also obtained. The paper also reaffirms the qualitative findings of earlier investigations which are, in fact, limiting cases of the present study.     

Forced Convection with Laminar Pulsating Counterflow in a Saturated Porous Circular Tube

An analytical solution is obtained for forced convection in a circular tube occupied by a core–sheath-layered saturated porous medium with counterflow produced by pulsating pressure gradients. The case of the constant heat-flux boundary conditions is considered, and the Brinkman model is employed for the porous medium. A perturbation approach is used to obtain analytical expressions for the velocity, temperature distribution, and transient Nusselt number for convection produced by an applied pressure gradient that fluctuates with small amplitude harmonically in time about a non-zero mean. It is shown that the fluctuating part of the Nusselt number alters in magnitude and phase as the dimensionless frequency increases. The magnitude increases from zero, goes through a peak, and then decreases to zero. The height of the peak depends on the values of various parameters. The phase (relative to that of the steady component) decreases from π/2 to − π/2 as the frequency increases.     

A Note on a Brinkman–Brinkman Forced Convection Problem

The problem of fully developed forced convection in a parallel-plates channel filled with a saturated porous medium (involving a Brinkman model for the momentum equation), with the effect of viscous dissipation (involving a Brinkman number), is discussed. Some general matters relating to the possibility of fully developed convection are also discussed.     

Double diffusive convection in a porous medium with modulated temperature on the boundaries

The effect of temperature modulation on the onset of double diffusive convection in a sparsely packed porous medium is studied by making linear stability analysis, and using Brinkman-Forchheimer extended Darcy model. The temperature field between the walls of the porous layer consists of a steady part and a time dependent periodic part that oscillates with time. Only infinitesimal disturbances are considered. The effect of permeability and thermal modulation on the onset of double diffusive convection has been studied using Galerkin method and Floquet theory. The critical Rayleigh number is calculated as a function of frequency and amplitude of modulation, Vadasz number, Darcy number, diffusivity ratio, and solute Rayleigh number. Stabilizing and destabilizing effects of modulation on the onset of double diffusive convection have been obtained. The effects of other parameters are also discussed on the stability of the system. Some results as the particular cases of the present study have also been obtained. Also the results corresponding to the Brinkman model and Darcy model have been compared.     

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.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

Mixed convection in the presence of horizontal impermeable surfaces in saturated porous media with variable permeability

Boundary layer approximation is applied for mixed convection about a horizontal flat plate in a saturated porous medium with aiding external flows. Similarity solutions are obtained, incorporating the variation of permeabilty, for 1) horizontal flat plate at zero angle of attack with constant heat flux; 2) stagnation point flows about horizontal flat plates with wall temperature varying asT _w ∝x ². The temperature and velocity profiles for different values of Ra/(RePr)^3/2 and the parameters governing the flow are obtained. The heat transfer rate is calculated and its implications in a geothermal application is discussed. Further, the criteria for pure mixed convection about horizontal flat plates in a porous media are established.     

Effect of Vertical Heterogeneity on Long-Term Migration of CO<Subscript>2</Subscript> in Saline Formations

For deep injection of CO₂ in thick saline formations, the movements of both the free gas phase and dissolved CO₂ are sensitive to variations in vertical permeability. A simple model for vertical heterogeneity was studied, consisting of a random distribution of horizontal impermeable barriers with a given overall volume fraction and distribution of lengths. Analytical results were obtained for the distribution of values for the permeability, and compared to numerical simulations of deep CO₂ injection and convection in heterogeneous formations, using multiple realizations for the permeability distribution. It is shown that for a formation of thickness H, the breakthrough times in two dimensions for deep injection scale as H ² for moderate injection rates. In comparison to heterogeneous shale distributions, a homogeneous medium with equivalent effective vertical permeability has a longer breakthrough time for deep injection, and a longer onset time for convection.     

On the flow of fluids through inhomogeneous porous media due to high pressure gradients

Most porous solids are inhomogeneous and anisotropic, and the flows of fluids taking place through such porous solids may show features very different from that of flow through a porous medium with constant porosity and permeability. In this short paper we allow for the possibility that the medium is inhomogeneous and that the viscosity and drag are dependent on the pressure (there is considerable experimental evidence to support the fact that the viscosity of a fluid depends on the pressure). We then investigate the flow through a rectangular slab for two different permeability distributions, considering both the generalized Darcy and Brinkman models. We observe that the solutions using the Darcy and Brinkman models could be drastically different or practically identical, depending on the inhomogeneity, that is, the permeability and hence the Darcy number.     

Convection in a porous medium with velocity slip and temperature jump boundary conditions

The onset of convection in a rarefield gas saturating a horizontal layer of a porous medium has been investigated using both Darcy and Brinkman models. It is assumed that due to rarefaction both velocity slip and temperature jump exist at the boundaries. The results show that (i) when the degree of rarefaction increases the critical Rayleigh number as well as the critical wave number for the onset of convection increases, (ii) stabilizing effect of temperature jump is more than that of velocity slip, (iii) Darcy model is seen to be the most stable one when compared to Brinkman model or the pure gaseous layer (i.e. in the absence of porous medium).