Unsteady natural convection heat and mass transfer in a saturated porous enclosure

A detailed numerical study has been performed to investigate transient natural convection heat and mass transfer in a porous enclosure. Major dimensionless groups governing the present problem areRa,N,Le, φ andAr. Results are particular presented to illustrate the effects of the combined thermal and solutal buoyancy forces on the temporal evolution of local/average Nusselt and Sherwood numbers. The results show that with the increase in the Rayleigh number, the heat and mass transfer is enhanced as a result of greater buoyancy effect. Additionally, the increase in buoyancy ratioN results in an improvement in the heat and mass transfer rates and in the mean time causes a short time duration for the flow to approach the steady-state condition.     

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

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

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

Magnetohydrodynamic natural convection heat transfer from radiate vertical surfaces with fluid suction or injection is considered. The nonsimilarity parameter is found to be the conductive fluid injection or suction along the streamwise coordinate = V{4x/² g(T _w – T )}^1/4. Three dimensionless parameters had been found to describe the problem: the magnetic influence number N = B ² _y /V ², the radiation-conduction parameter R _d = k _R /4aT ³ , and the Gebhart number Ge _x = gx/cp to represent the effect of the viscous dissipation. It is found that increasing the magnetic field strength causes the velocity and the heat transfer rates inside the boundary layer to decrease. Its apparent that increasing the radiation-conduction parameter decreases the velocity and enhances the heat transfer rates. The Gebhart number, i.e, the viscous dissipation had no effect on the present problem.Nomenclature a Stefan-Boltzmann constant - B _y Magnetic field flux density Wb/m² - Cf _x Local skin friction factor - c _p Specific heat capacity - f Dimensionless stream function - Ge _x Gebhart number, gx/cp - g Gravitational acceleration - k Thermal Conductivity - L Length of the plate - N Magnetic influence number, B ² _y /V ² - p Pressure - Pr Prandtl number - q _r Radiative heat flux - q _w (x) Local surface heat flux - Q _w (x) Dimensionless Local surface heat flux - R _d Planck number (Radiation-Conduction parameter), k _R /4aT ³ - T Temperature - T Free stream temperature - T _w Wall temperature - u, v Velocity components in x- and y-directions - V Porous wall suction or injection velocity - V _w Porous wall suction or injection velocity - x, y Axial and normal coordinates - Thermal diffusivity Greek symbols _R Roseland mean absorption coefficient, 4/3R _d - Coefficient of thermal expansion - Nonsimilarity parameter, V{4x/² g(T _w – T )}^1/4 - Peseudo-similarity variable - Dimensionless temperature - _w Ratio of surface temperature to the ambient temperature, T _w /T - Dynamice viscosity - Kinemtic viscosity - Fluid density - Electrical conductivity - _w Local wall shear stress - Dimensional stream function     

Mixed convection heat transfer from a horizontal plate with variable surface heat flux in a porous medium

In this article nonsimilarity solution for mixed convection from a horizontal surface in a saturated porous medium was obtained for the case of variable surface heat flux. The entire mixed convection regime, ranging from pure forced convection to pure free convection, is considered by introducing a single nonsimilarity parameter. Heat transfer results are predicted by employing four different flow models, namely, Darcy's law, the Ergun model, and the Brinkman-Forchheimer-extended Darcy model with constant and variable porosity. The variable porosity effect is approximated by an exponential function. Effects of transverse thermal dispersion are taken into consideration in the energy equation, along with variable stagnant thermal conductivities. The formulation of the present problem shows that the flow and heat transfer characteristics depend on five parameters, that is, the power in the variation of surface heat flux, the nonsimilarity mixed-convection parameter, the inertia effect parameter, the boundary effect parameter, and the ratio of thermal conductivity of the fluid phase to that of the solid phase. Numerical results for the local Nusselt number variations, based on the various flow models, are presented for the entire mixed convection regime. The impacts␣of different governing parameters on the heat transfer results are thoroughly investigated. Received on 7 August 1997     

Natural convection due to heat and mass transfer in a composite system

A numerical study is performed to analyse heat and mass transfer phenomena due to natural convection in a composite cavity containing a fluid layer overlying a porous layer saturated with the same fluid. The flow in the porous region is modelled using Brinkman–Forchheimer-extended Darcy model that includes both the effect of macroscopic shear (Brinkman effect) and flow inertia (Forchheimer effect). The vertical walls of the two-dimensional enclosure are isothermal whilst the horizontal walls are adiabatic. The two regions are coupled by equating the velocity and stress components at the interface. The resulting coupled equations in non-dimensional form are solved by an alternating direction implicit method by transforming them into parabolic form by the addition of false transient terms. The numerical results show that the amount of fluid penetration into the porous layer depends strongly upon the Darcy, thermal and solutal Rayleigh numbers. Average Nusselt number decreases while average Sherwood number increases with an increase of the Lewis number. The transfer of heat and mass on the heated wall near the interface depends strongly on the Darcy number. Received on 11 May 1998     

Natural convection heat transfer in the annular region between porous confocal ellipses

The Darcy–Boussinesq equations are solved in two dimensions and in elliptical cylindrical co‐ordinates using a second‐order‐accurate finite difference code and a very fine grid. For the limiting case of a circular geometry, the results show that a hysteresis loop is possible for some values of the radius ratio, in agreement both with previous calculations using cylindrical co‐ordinates and with the available experimental data. For the general case of an annulus of elliptical cross‐section, two configurations, blunt or slender, are considered. When the major axes are horizontal (blunt case) a hysteresis loop appears for a certain range of Raleigh numbers. For the slender configuration, when the major axes are vertical, a transition from a steady to a periodic regime (Hopf bifurcation) has been evidenced. In all cases, the heat transfer rate from the slender geometry is greater than that obtained in the blunt case. Copyright © 1999 John Wiley & Sons, Ltd.     

Natural convection in partitioned enclosed spaces of solar collector

The two-dimensional, steady, laminar natural convection phenomena in partitioned enclosure of solar collector has been studied numerically. Heat conduction along the partition is considered. An iterative finite-difference scheme is employed to solve the governing equations in the flow field. The effects of Rayleigh number, thermal conductivity ratio, partition angle, tilt angle, and aspect ratio on both the local and average heat transfer coefficients of the solar collector have been discussed. The range of Rayleigh number tested was up to 5 × 10⁴, the thermal conductivity ratio of 4.5 and 30, partition angle from 10 deg to 170 deg, tilt angle from 10 deg to 90 deg, and aspect ratio varied between 0.2 and 10. The results indicate that the convective heat transfer is strongly affected with the aspect ratio of the enclosures.

---

Freie Konvektion in unterteilten Kammern von Solarkollektoren

Zusammenfassung Die zweidimensionale, stetige, laminare freie Konvektion in unterteilten Kammern von Solarkollektoren wurde numerisch untersucht. Die Wärmeübertragung entlang dieser Kammern wurde betrachtet. Ein iteratives Finite-Differenzen-Schema wurde angewandt um die Gleichungen, welche das Strömungsfeld beschreiben, zu lösen. Der Einfluß der Rayleigh-Zahl, der thermische Leitfähigkeit, des Kammerwinkels, des Neigungswinkels und der Längenverhältnisse auf die örtlichen und durchschnittlichen Wärmeübertragungskoeffizienten des Solarkollektors wurde diskutiert. Der Bereich der Rayleigh-Zahl erstreckte sich bis zu 5 × 10⁴, das Verhältnis der thermischen Leitfähigkeit betrug 4.5 und 30, der Kammerwinkel lag zwischen 10° und 170°, der Neigungswinkel zwischen 10° und 90° und das Längenverhältnis variierte zwischen 0.2 und 10. Die Ergebnisse beinhalten, daß die konvektive Wärmeübertragung sehr stark durch das Längenverhältnis der Kammern beeinflußt wird.

---

Nomenclature a slope of the partition with respect to the horizontal - A H/L=cell aspect ratio - A _w t/L=wall aspect ratio - g acceleration due to gravity - h local heat transfer coefficient - average heat transfer coefficient - H cell length - k thermal conductivity of fluid within the cell - k _w thermal conductivity of the cell wall - L plate spacing - Nu _f h L/k=local cell Nusselt number - L/k=average cell Nusselt number - overall average Nusselt number - qL/k _w t(T _h–T _c)=average wall Nusselt number - Pr /=Prandtl number - q heat transfer in the cell wall from the hot to cold plate per unit depth - Ra g L ³ T/=Rayleigh number - R _k k _w/k=ratio of wall thermal conductivity to that of the fluid - t thickness of cell wall - T _c cold plate temperature - T _f temperature in cell - T _h hot plate temperature - T _w temperature in cell wall - u, U dimension and dimensionless velocities inx-direction - v, V dimension and dimensionless velocities iny-direction - x distance measured from the bottom of the cell (Fig. 1) - X x/L=normalized distance ofx - y distance measured from hot plate (Fig. 1) - Y y/L=normalized distance ofy - x ₁ distance measured in wall (Fig. 1) - X ₁ x/L=normalized distance ofx ₁ Greek symbols thermal diffusivity of fluid - coefficient of volumetric expansion of fluid - partition angle with respect to the hot plate - _f T _f–T _c/T _h–T _c=dimensionless temperature in cell - _w T _w–T _c/T _h–T _c=dimensionless temperature in cell wall - kinematic viscosity of fluid - enclosure tilt angle from horizontal - dimensional vorticity - L ²/=dimensionless vorticity - dimensionless streamline     

Three dimensional free convection flow and heat transfer along a porous vertical plate

The free convection flow along a vertical porous plate with transverse sinusoidal suction velocity distribution is investigated. Due to this type of suction velocity at the plate the flow becomes three dimensional one. For the asymptotic flow condition, the wall shear stress in the direction of main flow for different values of buoyancy parameter G is obtained. For G=0, the skin friction in the direction of free stream and the rate of heat transfer from the plate to the fluid are given. It is found that these results differ from those obtained by Gersten and Gross.     

Unglazed selective absorber solar air collector: Heat exchange analysis

Unglazed solar air collectors show promise for applications such as ventilation air heating or crop drying. In this paper a mathematical model is developed to analyze the heat exchanges in an unglazed non-porous selective absorber air heater. It is shown that at quasi-steady state the energy balance equations of the components of the collector cascade into a single first order differential equation. The solution of this differential equation is written down as an explicit expression of the local temperature of the fluid flowing in the collector in terms of the time dependent solar intensity. The effect of various parameters such as the inlet fluid temperature, the mass flow rate, and the depth of the air channel on the thermal performances of the unglazed selective absorber collector are also studied. These performances are comparable to those of a conventional two glass covers air collector for low wind speeds. Received on 8 February 1999     

Free convection heat transfer from an isothermal wavy surface in a porous enclosure

The coupled streamfuction–temperature equations governing the Darcian flow and convection process in a fluid-saturated porous enclosure with an isothermal sinusoidal bottom sun face, has been numerically analyzed using a finite element method (FEM). No restrictions have been imposed on the geometrical non-linearity arising from the parameters like wave amplitude (a), number of waves per unit length (N), wave phase (Φ), aspect ratio (A) and also on the flow driving parameter Rayleigh number (Ra). The numerical simulations for varying values of Ra bring about interesting flow features, like the transformation of a unicellular flow to a multicellular flow. Both with increasing amplitude and increasing number of waves per unit length, owing to the shift in the separation and reattachment points, a row–column pattern of multicellular flow transforms to a simple row of multicellular flow. A cycle of n celluar and n+1 cellular flows, with the flow in adjacent cells in the opposite direction, periodically manifest with phase varying between 0 and 360°. The global heat transfer into the system has been found to decrease with increasing amplitude and increasing number of waves per unit length. Only marginal changes in the global heat flux are observed, either with increasing Ra or varying Φ. Effectively, sinusoidal bottom surface undulations of the isothermal wall of a porous enclosure reduces the heat transfer into the system. © 1998 John Wiley & Sons, Ltd.     

Forced convection heat transfer from a flat-plate model collector on roof of a model house

Presented in this abstract are the considerations and results of an investigation of forced convective heat loss from the surface of flat-plate model collectors flush mounted on the roof of a model residential house. Forced convection with air as the external fluid was treated.The study was motivated by the need to better understand heat transfer from the top surface of flat-plate solar energy collectors. The experimental work dealt with heat transfer from the top of rectangular parallelpipeds flush mounted on roof of a model residential house. The smallest collector was 7.62×10^–2 m long by 7.62×10^–2 m wide by 3.3×10^–2 m deep, and the largest collector was 15.875×10^–2 m long by 15.875×10^–2 m wide by 3.3×10^–2 m deep. These small scale units are expected to simulate heat loss pattern from a flat-plate solar collector placed over the roof of a residential dwelling. Experiments were performed to determine the average heat transfer coefficients for forced convective air flow over the surface of the model collector which was inclined and yawed relative to the incoming air flow. The wind speeds varied approximately from 2.5 m/sec to 15 m/sec, the flow was laminar for all wind speeds.

---

Konvektiver Wärmeübergang von Flachplatten-Kollektoren auf den Dächern von Wohnhäusern im Modellversuch

Zusammenfassung Es werden Vorüberlegungen und Ergebnisse bezüglich einer Untersuchung mitgeteilt, deren Ziel es war, die Wärmeverluste infolge Zwangskonvektion von Flachplatten-Kollektoren auf den Dächern von Wohnhäusern im Kleinmodell zu bestimmen. Angeregt wurde die Studie durch den Bedarf an genaueren Berechnungsunterlagen für den Wärmeübergang an Flachplatten-Solarkollektoren. Die für das Experiment verwendeten zwei Kollektoren waren quadratisch, mit den Seitenlägen 7,62 bzw. 15,9 cm, bei jeweils 3,3 cm Tiefe. Die Versuche lieferten unter Zwangskonvektion mittlere Wärmeübergangskoeffizienten an den Oberflächen der Modellkollektoren, welche zur ankommenden Luftströmung sowohl geneigt als auch gedreht orientiert waren. Bei laminarer Strömung lagen die Anströmgeschwindigkeiten zwischen 2,5 and 15 m/s.

     

Analytical and numerical study of natural convection heat transfer in an inclined composite enclosure

Laminar natural convection heat transfer in inclined fluid layers divided by a partition with finite thickness and conductivity is studied analytically and numerically. The governing equations for the fluid layers are solved analytically in the limit of a thin layered system with constant flux boundary conditions. The study covers of the range of Ra from 10³ to 10⁷, from 0° to 180° and the thermal conductivity ratio of partition to fluid ratioK from 10^–2 to 10⁶. The Prandtl number was 0.72 (for air). Results are obtained in terms of an overall Nusselt number as a function of Rayleigh number, angle of inclination of the system, mid layer thickness, and mid layer thermal conductivity. The critical Rayleigh number for the onset of convection in a bottom-heated horizontal system is predicted. The results are compared with the numerical results obtained by solving the complete system of governing equations, using SIMPLER method, as well as with the limiting cases in the literature.     

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     

Non-linear convection in a porous medium with internal heat sources

The paper studies non-linear thermal convection in a horizontal porous layer of fluid with nearly insulating boundaries and in the presence of internal heat sources. Square and hexagonal cells are found to be the only possible stable convection cells. Finite amplitude instability could exist for some particular forms of an internal heat source Q. For a uniform Q, the preferred flow pattern is that of hexagons for amplitude ε smaller than some critical value ε_c, while both squares and hexagonal cells are stable for ε ? ε_c. The convective motion is downward at the hexagonal cell's centers. For a non-uniform Q, the qualitative features of thermal convection depend on the actual form Q. In particular, a non-uniform Q can increase or decrease the cell's size and the critical Rayleigh number at the onset of convection, and stabilize or destabilize the convective motion in the form of hexagonal cells with either upward or downward motion at the cell's centers.     

Combined heat and mass transfer by mixed convection MHD flow along a porous plate with chemical reaction in presence of heat source

An exact and a numerical solutions to the problem of a steady mixed convective MHD flow of an incompressible viscous electrically conducting fluid past an infinite vertical porous plate with combined heat and mass transfer are presented.A uniform magnetic field is assumed to be applied transversely to the direction of the flow with the consideration of the induced magnetic field with viscous and magnetic dissipations of energy.The porous plate is subjected to a constant suction velocity as well as a uniform mixed stream velocity.The governing equations are solved by the perturbation technique and a numerical method.The analytical expressions for the velocity field,the temperature field,the induced magnetic field,the skin-friction,and the rate of heat transfer at the plate are obtained.The numerical results are demonstrated graphically for various values of the parameters involved in the problem.The effects of the Hartmann number,the chemical reaction parameter,the magnetic Prandtl number,and the other parameters involved in the velocity field,the temperature field,the concentration field,and the induced magnetic field from the plate to the fluid are discussed.An increase in the heat source/sink or the Eckert number is found to strongly enhance the fluid velocity values.The induced magnetic field along the x-direction increases with the increase in the Hartmann number,the magnetic Prandtl number,the heat source/sink,and the viscous dissipation.It is found that the flow velocity,the fluid temperature,and the induced magnetic field decrease with the increase in the destructive chemical reaction.Applications of the study arise in the thermal plasma reactor modelling,the electromagnetic induction,the magnetohydrodynamic transport phenomena in chromatographic systems,and the magnetic field control of materials processing.     

Natural convection momentum and heat transfer in saturated highly porous media — An asymptotic approach —

Starting from a clear flow situation with no porous matrix a regular perturbation analysis is applied to account for the influence of a highly porous matrix. The perturbation parameter is 1 –n,n being the porosity. By means of asymptotic formulae thex-dependence of the problem under consideration as well as most parameters of the problem can be separated. Thus only ordinary differential equations with the Prandtl number as a parameter have to be solved. Skin friction and heat transfer formulae are given asymptotically which compare well with literature data for highly porous media.Ausgehend von der Strömungssituation ohne poröse Matrix wird eine reguläre Störungsrechnung durchgeführt, die den Einfluß einer hoch porösen Matrix erfaßt. Der Störparameter ist 1 –n, wobein die Porosität ist. Mit Hilfe der asymptotischen Formulierung können diex-Abhängigkeit sowie eine Reihe von Parameterabbängigkeiten in dem Problem separiert werden. Auf diese Weise müssen nur gewöhnliche Differentialgleichungen gelöst werden, die als einzigen Lösungsparameter die Prandtl-Zahl besitzen. Die asymptotischen Ergebnisse für Impuls- und Wärmeübergang stimmen gut mit Literaturdaten zu hoch porösen Medien überein.     

Natural convection in composite systems with phase change materials

In this study, we carried out a numerical simulation of transient heat transfer in a composite passive system consisting of air–phase change material–air, arranged as a rectangular enclosure. The vertical boundaries of the enclosure are isothermal and the horizontal ones adiabatic. The enthalpy formulation with a fixed grid is used to study the process of phase change with liquid–solid interface zone controlled by natural convection. The flow in this zone is simulated by a model based on the Darcy porous medium. The numerical solution of the mathematical model is done using finite difference–control volume algorithm. The influence of the geometrical and thermal parameters is studied. It is found that subcooling coefficient is the most important parameter influencing heat transfer, and for a given subcooling, there is an optimum phase change partition thickness.