Conjugate Natural Convection in a Porous Enclosure with Non-Uniform Heat Generation

Effects of a conductive wall on natural convection in a square porous enclosure having internal heating at a rate proportional to a power of temperature difference is studied numerically in this article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (0 ???Ra ???1000), the internal heating and the local exponent parameters (0 ????? ???5), (1 ????? ???3), the wall to porous thermal conductivity ratio (0.44 ???Kr ???9.9) and the ratio of wall thickness to its width (0.02 ???D ???0.5). The results are presented to show the effect of these parameters on the fluid flow and heat transfer characteristics. It is found a strong internal heating can generate significant maximum fluid temperature more than the conductive solid wall. Increasing value thermal conductivity ratio and/or decreasing the thickness of solid wall can increase the maximum fluid temperature. It is also found that at very low Rayleigh number, the heat transfer across the porous enclosure remain stable for any values of the thermal conductivity ratio.     

Effect of Conduction in Bottom Wall on Darcy–Bénard Convection in a Porous Enclosure

Conjugate natural convection-conduction heat transfer in a square porous enclosure with a finite-wall thickness is studied numerically in this article. The bottom wall is heated and the upper wall is cooled while the verticals walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and the COMSOL Multiphysics software is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (100 ≤ Ra ≤ 1000), the wall to porous thermal conductivity ratio (0.44 ≤ K _r ≤ 9.90) and the ratio of wall thickness to its height (0.02 ≤ D ≤ 0.4). The results are presented to show the effect of these parameters on the heat transfer and fluid flow characteristics. It is found that the number of contrarotative cells and the strength circulation of each cell can be controlled by the thickness of the bottom wall, the thermal conductivity ratio and the Rayleigh number. It is also observed that increasing either the Rayleigh number or the thermal conductivity ratio or both, and decreasing the thickness of the bounded wall can increase the average Nusselt number for the porous enclosure.     

Effect of an Inserted Porous Layer Located at a Wall of a Parallel Plate Channel on Forced Convection Heat Transfer

A theoretical study is performed on heat and fluid flow in partially porous medium filled parallel plate channel. A uniform symmetrical heat flux is imposed onto the boundaries of the channel partially filled with porous medium. The dimensional forms of the governing equations are solved numerically for different permeability and effective thermal conductivity ratios. Then, the governing equations are made dimensionless and solved analytically. The results of two approaches are compared and an excellent agreement is observed, indicating correctness of the both solutions. An overall Nusselt number is defined based on overall thermal conductivity and difference between the average temperature of walls and mean temperature to compare heat transfer in different channels with different porous layer thickness, Darcy number, and thermal conductivity ratio. Moreover, individual Nusselt numbers for upper and lower walls are also defined and obtained. The obtained results show that the maximum overall Nusselt number is achieved for thermal conductivity ratio of 1. At specific values of Darcy number and thermal conductivity ratio, individual Nusselt numbers approach to infinity since the value of wall temperatures approaches to mean temperature.     

Stationary Fingering Instability in a Non-equilibrium Porous Medium with Coupled Molecular Diffusion

Finger type double diffusive convective instability in a fluid-saturated porous medium is studied in the presence of coupled heat-solute diffusion. A local thermal non-equilibrium (LTNE) condition is invoked to model the Darcian porous medium which takes into account the energy transfer between the fluid and solid phases. Linear stability theory is implemented to compute the critical thermal Rayleigh number and the corresponding wavenumber exactly for the onset of stationary convection. The effects of Soret and Dufour cross-diffusion parameters, inter-phase heat transfer coefficient and porosity modified conductivity ratio on the instability of the system are investigated. The analysis shows that positive Soret mass flux triggers instability and positive Dufour energy flux enhances stability whereas their combined influence depends on the product of solutal Rayleigh number and Lewis number. It also reveals that cell width at the convection threshold gets affected only in the presence of both the cross-diffusion fluxes. Besides, asymptotic solutions for both small and large values of the inter-phase heat transfer coefficient and porosity modified conductivity ratio are found. An excellent agreement is found between the exact and asymptotic solutions.     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

Laminar natural convection in inclined enclosures bounded by a solid wall

Laminar natural convection has been studied in enclosures bounded by a solid wall with its outer boundary at constant temperature while the opposing side has a constant heat flux. Two-dimensional equations of conservation of mass, momentum and energy, with the Boussinesq approximation are solved using a finite difference method. The numerical procedure adopted is based on the SIMPLER algorithm. Various parameters were: Rayleigh number (from 10³ to 10⁶), dimensionless conductivity of bounding wall (from 1 to 10) and dimensionless wall width (from 0.15 to 0.5), aspect ratio (from 0.5 to 1) and the inclination angle (from 30^° to 180^°). The results are reduced in terms of the normalized Nusselt number as a function of the Rayleigh number, and other dimensionless parameters. The isotherms and streamlines are produced for various Rayleigh numbers and geometrical conditions. It is found that the heat transfer is an increasing function of the Rayleigh number, wall to fluid conductivity ratio, enclosure aspect ratio and a decreasing function of the wall thickness. It passes from a maximum for the inclination angle of about 80^°.     

Natural convection in an inclined rectangular porous cavity with uniform heat flux from the side

The problem of natural convection in an inclined rectangular porous layer enclosure is studied numerically. The enclosure is heated from one side and cooled from the other by a constant heat flux while the two other walls are insulated. The effect of aspect ratio, inclination angle and Rayleigh number on heat transfer is studied. It is found that the enclosure orientation has a considerable effect on the heat transfer. The negative orientation sharply inhibits the convection and consequently the heat transfer and a positive orientation maximizes the energy transfer. The maximum temperature within the porous medium can be considerably higher than that induced by pure conduction when the cavity is negatively oriented. The peak of the average Nusselt number depends on the Rayleigh number and the aspect ratio. The heat transfer between the two thermally active boundaries is sensitive to the effect of aspect ratio. For an enclosure at high or low aspect ratio, the convection is considerably decreased and the heat transfer depends mainly on conduction.     

Lattice Boltzmann Pore Scale Simulation of Natural Convection in a Differentially Heated Enclosure Filled with a Detached or Attached Bidisperse Porous Medium

In this research, pore scale simulation of natural convection in a differentially heated enclosure filled with a conducting bidisperse porous medium is investigated using the thermal lattice Boltzmann method. For the first time, the effect of connection of the bidisperse porous medium to the enclosure walls is studied by considering the attached geometry in addition to the detached one. Effect of most relevant parameters on the streamlines and isotherms as well as hot wall average Nusselt number is studied for two of the bidisperse porous medium configurations. It is observed that effect of geometrical and thermo-physical parameters of the bidisperse porous medium on the heat transfer characteristics is more complicated for the attached configuration. To assess the validity of the local thermal equilibrium condition in the micro-porous media, the pore scale results are used to compute the percentage of the local thermal non-equilibrium for two of the bidisperse porous medium configurations. It is concluded that for the detached configuration, the local thermal equilibrium condition is confirmed in the entire micro-porous media for the ranges of the parameters studied here. However, for the attached geometry, it is shown that departure from the local thermal equilibrium condition is observed for the higher values of the Rayleigh number, micro-porous porosity, solid–fluid thermal conductivity ratio, and the smaller values of the macro-pores volume fraction.     

Free Convection Heat Transfer From A Sphere In A Porous Medium Using A Thermal Non-equilibrium Model

This work studies the free convection heat transfer from a sphere with constant wall temperature embedded in a fluid-saturated porous medium using a thermal non-equilibrium model. The governing equations are transformed into boundary-layer partial differential equations by the coordinate transform, and the obtained governing equations are then solved by the cubic spline collocation method. The temperature distributions for fluid and solid phases are shown for different values of the porosity scaled thermal conductivity ratio, the interphase heat transfer parameter, and the streamwise coordinate. The effects of the porosity scaled thermal conductivity ratio and the interphase heat transfer parameter between solid and fluid phases on the local Nusselt numbers for fluid and solid phases are examined. Results show the local Nusset number for the porous medium can be increased by increasing the porosity scaled thermal conductivity ratio. Moreover, the thermal non-equilibrium effect is more significant for low values of the porosity scaled thermal conductivity ratio or the interphase heat transfer parameter.     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

The Use of Porous Fins for Heat Transfer Augmentation in Parallel-Plate Channels

In this study, a steady, fully developed laminar forced convection heat augmentation via porous fins in isothermal parallel-plate duct is numerically investigated. High-thermal conductivity porous fins are attached to the inner walls of two parallel-plate channels to enhance the heat transfer characteristics of the flow under consideration. The Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous fins. This study reports the effect of several operating parameters on the flow hydrodynamics and thermal characteristics. This study demonstrates, mainly, the effects of porous fin thickness, Darcy number, thermal conductivity ratio, Reynolds number, and microscopic inertial coefficient on the thermal performance of the present flow. It is found that the highest Nusselt number is achieved at fully filled porous duct which requires the highest pumping pressure. The results show that using porous fins requires less pumping pressure with comparable high heat augmentation weight against fully filled porous duct. It is found that higher Nusselt numbers are achieved by increasing the microscopic inertial coefficient (A), the Reynolds number (Re), and the thermal conductivity of the porous substrate k ₂. The results show that heat transfer can be enhanced (1) with the use of high thermal conductivity fins, (2) by decreasing the Darcy number, and (3) by increasing microscopic inertial coefficient.     

多层平行裂隙型多孔介质通道内流体流动传热特性

基于Brinkman-extended Darcy模型和局部热平衡模型,对多层平行裂隙型多孔介质通道内的流动传热特性进行研究.获得了多层平行裂隙型多孔介质通道内各区域的速度场、温度场、摩擦系数及努塞尔数解析解,并分析了裂隙层数、达西数、空心率、有效热导率之比等对通道内流动传热特性的影响.结果表明:达西数较小时,通道多孔介质层内会出现不随高度变化的达西速度,此达西速度会随裂隙层数的增加而增大,但却不受各裂隙层下多孔介质层位置变化的影响.增加裂隙层数会减弱空心率对压降的影响,会使通道内流体压降升高,但升高程度会逐渐降低.增大热导率之比或减小空心率会使多裂隙通道内出现阶梯式温度分布,而在较小热导率之比或较大空心率时多裂隙情况下的温度分布曲线会趋于一致.此外,当热导率之比较小时,多层裂隙通道内的传热效果在任何空心率下都要优于单裂隙情况,当热导率之比较大时,存在临界空心率使各裂隙层数通道内的传热效果相同,且多裂隙通道内继续增加裂隙层数对传热强度影响不大.     

Non-Darcian Effects on Buoyancy-Induced Heat Transfer in a Partially Divided Square Enclosure with Internal Heat Generation

The present investigation addresses non-Darcian effects on the buoyancy-induced heat transfer in a partially divided square enclosure with internal heat generation. The generalized model of the momentum equation, which is also known as the Forchheimer–Brinkman extended Darcy model, which takes into account boundary and inertia effects, was used in representing the fluid motion inside the porous layer. The local thermal equilibrium condition was assumed to be valid for the range of the thermophysical parameters considered in the present investigation. The transport equations were solved using the finite element formulation based on the Galerkin method of weighted residuals. The validity of the numerical code used was ascertained by comparing our results with previously published results. Results were obtained in terms of streamlines, isotherms, and Nusselt number for various geometrical parameters specifying the height and width of the partition. In addition, the effects of external and internal Rayleigh numbers and Darcy number were highlighted in the proposed study.     

Linear and nonlinear double diffusive convection in a rotating sparsely packed porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and nonlinear stability analyses. The Brinkman model that includes the Coriolis term is employed as the momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for the energy equation. The onset criterion for stationary, oscillatory, and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal diffusion, solute diffusion, and rotation that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number, and Taylor number on the stability of the system is investigated. The nonlinear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

Conjugate natural convection of non-Newtonian hybrid nanofluid in wavy-shaped enclosure

The present study concerns the modelization and numerical simulation for the heat and flow exchange characteristics in a novel configuration saturated with a nonNewtonian Ag-MgO hybrid nanofluid. The wavy shaped enclosure is equipped with onequarter of a conducting solid cylinder. The system of equations resulting from the mathematical modeling of the physical problem in its dimensionless form is discretized via the higher-order Galerkin-based finite element method(GFEM). The dependency of vario...     

MHD natural convection with Joule and viscous heating effects in iso-flux porous medium-filled enclosures

In this study, the magnetohydrodynamics （MHD） natural convection heat transfer with Joule and viscous heating effects inside an iso-flux porous medium-filled inclined rectangular enclosure is studied numerically. An iso-heat flux is applied for heating and cooling the two opposing walls of the enclosure while the other walls are adiabatic. The Forchheimer extension of Darcy-Oberbeck-Boussinesq and energy equations is transformed into a dimensionless form using a set of suitable variables instead of a finite difference scheme. The governing parameters are the magnetic influence number, the modified Rayleigh number, the inclination angle, and the aspect ratio of the enclosure. The results show that viscous and Joule heating effects decrease heat transfer rates.     

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.     

Natural convection and conduction heat transfer in open shallow cavities with bounding walls

A numerical study has been carried out on inclined open shallow cavities, which are formed by a wall and horizontal fins. Constant heat flux is applied on the surface of the wall inside the cavity while its other surface was kept isothermal. The wall and the fins are conductive. Conjugate heat transfer by natural convection and conduction is studied by numerically solving equations of mass, momentum and energy. Streamlines and isotherms are produced, heat and mass transfer is calculated. A parametric study is carried out using following parameters: Rayleigh number from 10⁶ to 10¹², conductivity ratio from 1 to 60, open cavity aspect ratio from 1 to 0.125, dimensionless end wall thickness from 0.05 to 0.20, horizontal walls from 0.01 to 0.15 and inclination of the end wall from 90° to 45°. It is found that the volume flow rate and Nusselt number are a decreasing function of the cavity aspect ratio, horizontal fin thickness and conductivity ratio. They are an increasing function of end wall thickness and inclination angle, except in the latter case optima exist at high Rayleigh numbers.     

Numerical Simulation of Natural Convection in Wavy Porous Cavities Under the Influence of Thermal Radiation Using a Thermal Non-equilibrium Model

Continuum equations governing thermal non-equilibrium modeling of steady natural convection inside wavy enclosures with the effect of thermal radiation are developed. These equations account for such effects as the inter-phase heat transfer coefficient effect, the thermal radiation effect, the modified conductivity ratio effect and the Rayleigh number effect. Finite difference method is employed to solve these equations and comparisons between previous published works are presented. Numerical results for the flow and heat transfer for the fluid and solid phases are obtained for various combinations of the physical parameters. Graphical and tabular results illustrating interesting features of the physics of the problem are presented and discussed.