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.     

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.     

The horizontal intrusion layer of melt in a saturated porous medium

This article presents a theory of how the melt region advances as an intrusion layer along the top boundary of a solid phase-change material that is heated from the side. The phase-change material fills the pores of a solid matrix. We show that the thickness of the horizontal melt layer increases as x^3/5, where x is the horizontal distance measured by from the leading edge of the layer. The total length of the intrusion layer increases as t^3/4, and as T_max^5/4. Finite-difference simulations of convection melting in the Darcy-Rayleigh number range of 200–800 agree with the theoretical results. We also show that in a rectangular porous medium heated from the side, the size of the entire melt region is dominated by the melting contributed by the horizontal intrusion layer, if the time is great enough so that the group (Ste Fo)^3/4 is greater than 1.     

MHD natural convection in porous media-filled enclosures

In this work, the magnetohydrodynamics （MHD） natural convection heat transfer problem inside a porous medium filled with inclined rectangular enclosures is investigated numerically. The boundary conditions selected on the enclosure are two adiabatic and two isothermal walls. The governing equations, continuity, and Forchheimer extension of the Darcy law and energy are transformed into dimensionless forms by using a set of suitable variables, and then solved by using a finite difference scheme. The governing parameters are the magnetic influence number, the Darcy Rayleigh number, the inclination angle, and the aspect ratio of the enclosure. It is found that the magnetic influence number and the inclination angle have pronounced effects on the fluid flow and heat transfer in porous media-filled enclosures.     

Magnetohydrodynamics convective flow in a vertical slot through a porous medium

An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman model. The fully-developed two-dimensional flow from capped to open ends is considered for which a continuum of solutions is obtained. The influence of pertinent parameters on the flow is delineated and appropriate conclusions are drawn. The asymptotic behaviour and the volume flux are analyzed and incorporated graphically for the three-parameter family of solution.     

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.     

Numerical investigation of Dufour and Soret effects on unsteady MHD natural convection flow past vertical plate embedded in non-Darcy porous medium

The Dufour and Soret effects on the unsteady two-dimensional magnetohydro-dynamics(MHD) double-diffusive free convective flow of an electrically conducting fluidpast a vertical plate embedded in a non-Darcy porous medium are investigated numeri-cally.The governing non-linear dimensionless equations are solved by an implicit finitedifference scheme of the Crank-Nicolson type with a tridiagonal matrix manipulation.The effects of various parameters entering into the problem on the unsteady dimension-less velocity,temperature,and concentration profiles are studied in detail.Furthermore,the time variation of the skin friction coefficient,the Nusselt number,and the Sherwoodnumber is presented and analyzed.The results show that the unsteady velocity,tem-perature,and concentration profiles are substantially influenced by the Dufour and Soreteffects.When the Dufour number increases or the Soret number decreases,both the skinfriction and the Sherwood number decrease,while the Nusselt number increases.It isfound that,when the magnetic parameter increases,the velocity and the temperaturedecrease in the boundary layer.     

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.     

Numerical convergence studies of the mixed finite element method for natural convection flow in a fluid-saturated porous medium

This paper describes a numerical approximation scheme for the natural convection (NC) flow in a fluid-saturated porous medium. Our formulation of the problem is based on the mixed finite element method (FEM). Using the so-called consistent splitting scheme, a second-order accuracy in time and in space is verified by the numerical calculation. The resulting flow patterns and heat transfer for different Rayleigh numbers, Darcy numbers and porosities are numerically studied by the proposed scheme.     

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.     

Conduction-radiation effect on natural convection flow in fluid-saturated non-Darcy porous medium enclosed by non-isothermal walls

The combined effect of conduction-convection-radiation on natural convection flow of an optically thick Newtonian fluid with gray radiant properties, confined in a porous media square cavity with Darcy-Brinkman-Forchheimer drag is studied numerically. For a gray fluid, Rosseland diffusion approximation is considered. It is assumed that (i) the temperature of the left vertical wall varies linearly with height, (ii) the right vertical and top walls are at a lower temperature, and (iii) the bottom wall is uniformly-heated. The governing equations are solved using the alternate direct implicit method together with the successive over relaxation technique. The investigation of the effect of governing parameters, namely, the Forschheimer resistance (Γ), the temperature difference (Δ), and the Plank number (Rd), on the flow pattern and heat transfer characteristics is carried out. It can be seen that the reduction of flow and heat transfer occur as the Forschheimer resistance is increased. On the other hand, both the flow strength and heat transfer increase as the temperature ratio Δ is increased.     

Darcy-Forchheimer natural,forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media

The governing equation for Darcy-Forchheimer flow of non-Newtonian inelastic power-law fluid through porous media has been derived from first principles. Using this equation, the problem of Darcy-Forchheimer natural, forced, and mixed convection within the porous media saturated with a power-law fluid has been solved using the approximate integral method. It is observed that a similarity solution exists specifically for only the case of an isothermal vertical flat plate embedded in the porous media. The results based on the approximate method, when compared with existing exact solutions show an agreement of within a maximum error bound of 2.5%.Nomenclature A cross-sectional area - b _i coefficient in the chosen temperature profile - B ₁ coefficient in the profile for the dimensionless boundary layer thickness - C coefficient in the modified Forchheimer term for power-law fluids - C ₁ coefficient in the Oseen approximation which depends essentially on pore geometry - C _i coefficient depending essentially on pore geometry - C _D drag coefficient - C _t coefficient in the expression forK ^* - d particle diameter (for irregular shaped particles, it is characteristic length for average-size particle) - f _p resistance or drag on a single particle - F _R total resistance to flow offered byN particles in the porous media - g acceleration due to gravity - g _x component of the acceleration due to gravity in thex-direction - Grashof number based on permeability for power-law fluids - K intrinsic permeability of the porous media - K ^* modified permeability of the porous media for flow of power-law fluids - l _c characteristic length - m exponent in the gravity field - n power-law index of the inelastic non-Newtonian fluid - N total number of particles - Nux,_D,F local Nusselt number for Darcy forced convection flow - Nux,_D-F,F local Nusselt number for Darcy-Forchheimer forced convection flow - Nux,_D,M local Nusselt number for Darcy mixed convection flow - Nux,_D-F,M local Nusselt number for Darcy-Forchheimer mixed convection flow - Nux,_D,N local Nusselt number for Darcy natural convection flow - Nux,_D-F,N local Nusselt number for Darcy-Forchheimer natural convection flow - pressure - p exponent in the wall temperature variation - _Pe c characteristic Péclet number - _Pe x local Péclet number for forced convection flow - _Pe x modified local Péclet number for mixed convection flow - _Ra c characteristic Rayleigh number - _Ra x local Rayleigh number for Darcy natural convection flow - _Ra x local Rayleigh number for Darcy-Forchheimer natural convection flow - Re convectional Reynolds number for power-law fluids - Reynolds number based on permeability for power-law fluids - T temperature - T _e ambient constant temperature - T w,_ref constant reference wall surface temperature - T _w(X) variable wall surface temperature - T _w temperature difference equal toT w,_ref–T _e - T ₁ term in the Darcy-Forchheimer natural convection regime for Newtonian fluids - T ₂ term in the Darcy-Forchheimer natural convection regime for non-Newtonian fluids (first approximation) - T _N term in the Darcy/Forchheimer natural convection regime for non-Newtonian fluids (second approximation) - u Darcian or superficial velocity - u ₁ dimensionless velocity profile - u _e external forced convection flow velocity - u _s seepage velocity (local average velocity of flow around the particle) - u _w wall slip velocity - U c _M characteristic velocity for mixed convection - U c _N characteristic velocity for natural convection - x, y boundary-layer coordinates - x ₁,y ₁ dimensionless boundary layer coordinates - X coefficient which is a function of flow behaviour indexn for power-law fluids - effective thermal diffusivity of the porous medium - shape factor which takes a value of/4 for spheres - shape factor which takes a value of/6 for spheres - ₀ expansion coefficient of the fluid - _T boundary-layer thickness - T ₁ dimensionless boundary layer thickness - porosity of the medium - similarity variable - dimensionless temperature difference - coefficient which is a function of the geometry of the porous media (it takes a value of 3 for a single sphere in an infinite fluid) - ₀ viscosity of Newtonian fluid - ^* fluid consistency of the inelastic non-Newtonian power-law fluid - constant equal toX(2 ^2–n )/ - density of the fluid - dimensionless wall temperature difference     

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.     

Non-Darcy natural convection from a vertical wavy surface in a porous medium

We examine the combined effect of spatially stationary surface waves and the presence of fluid inertia on the free convection induced by a vertical heated surface embedded in a fluid-saturated porous medium. We consider the boundary-layer regime where the Darcy-Rayleigh number, Ra, is very large, and assume that the surface waves have O(1) amplitude and wavelength. The resulting boundary-layer equations are found to be nonsimilar only when the surface is nonuniform and inertia effects are present; self-similarity results when either or both effects are absent. Detailed results for the local and global rates of heat transfer are presented for a range of values of the inertia parameter and the surface wave amplitude.     

Natural convection of nanofluid over vertical plate embedded in porous medium: prescribed surface heat flux

The aim of the present paper is to analyze the natural convection heat and mass transfer of nanofluids over a vertical plate embedded in a saturated Darcy porous medium subjected to surface heat and nanoparticle fluxes. To carry out the numerical solution, two steps are performed. The governing partial differential equations are firstly simplified into a set of highly coupled nonlinear ordinary differential equations by appropriate similarity variables, and then numerically solved by the finite difference method. The obtained similarity solution depends on four non-dimensional parameters, i.e., the Brownian motion parameter (N _b), the Buoyancy ratio (N _r), the thermophoresis parameter (N _t), and the Lewis number (Le). The variations of the reduced Nusselt number and the reduced Sherwood number with N _b and N _t for various values of Le and N _r are discussed in detail. Simulation results depict that the increase in N _b, N _t, or N _r decreases the reduced Nusselt number. An increase in the Lewis number increases both of the reduced Nusselt number and the Sherwood number. The results also reveal that the nanoparticle concentration boundary layer thickness is much thinner than those of the thermal and hydrodynamic boundary layers.     

A boundary layer analysis of combined heat and mass transfer by natural convection from a concentrated source in a saturated porous medium

A boundary layer analysis was carried out to investigate the coupled phenomena of heat and mass transfer by natural convection from concentrated heat and mass sources embedded in saturated porous media. Both line and point source problems were treated. The boundary layer equations based on Darcy's law and Boussinesq approximation were solved by means of similarity transformation to obtain the details of velocity, temperature and concentration distributions above a concentrated heat source. Two important parameters, namely the Lewis number Le and the buoyancy ratioN were identified to conduct a series of numerical integrations. For the case of small Le, a substance diffuses further away from the plume centerline, such that the mass transfer influences both velocity and temperature profiles over a wide range. For large Le, on the other hand, the substance diffuses within a narrow range along the centerline. Naturally, the influence of mass transfer is limited to the level of the centerline velocity, so that a peaky velocity profile appears for positiveN whereas a velocity defect emerges along the centerline for negativeN. For such cases of large Le, the temperature profiles are found to be fairly insensitive to Le.     

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

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

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.     

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.     

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.