Double-diffusive natural convection in inclined porous cavities with various aspect ratios and temperature-dependent heat source or sink

Laminar double-diffusive natural convective flow of a binary fluid mixture in inclined square and rectangular cavities filled with a uniform porous medium in the presence of temperature-difference dependent heat generation (source) or absorption (sink) is considered. Transverse gradients of heat and mass are applied on two opposing walls of the cavity while the other two walls are kept adiabatic and impermeable to mass transfer. The problem is put in terms of the stream function-vorticity formulation. A numerical solution based on the finite-difference methodology is obtained for relatively high Lewis numbers. Representative results illustrating the effects of the inclination angle of the cavity, buoyancy ratio, Darcy number, heat generation or absorption coefficient and the cavity aspect ratio on the contour maps of the streamline, temperature, and concentration as well as the profiles of velocity, temperature and concentration at mid-section of the cavity are reported. In addition, numerical results for the average Nusselt and Sherwood numbers are presented for various parametric conditions and discussed.     

Numerical investigation of thermosolutal natural convection in a rectangular enclosure of an aspect ratio four with heat and solute sources

A numerical study of double-diffusive natural convection in an enclosure with a partial vertical heat and mass sources for an aspect ratio Ar = 4 has been carried out. The influence of various dimensionless parameters (Rayleigh number, buoyancy ratio, source location, Lewis number, and source length) on the flow behavior are investigated. Correlations of average Nusselt and Sherwood numbers are obtained as function of two parameters (Ra, d) and (Le, d), respectively.     

A Modified Nested-Gridding for Upscaling–Downscaling in Reservoir Simulation

This article reports a numerical study of double-diffusive convection in a fluid-saturated vertical porous annulus subjected to discrete heat and mass fluxes from a portion of the inner wall. The outer wall is maintained at uniform temperature and concentration, while the top and bottom walls are adiabatic and impermeable to mass transfer. The physical model for the momentum equation is formulated using the Darcy law, and the resulting governing equations are solved using an implicit finite difference technique. The influence of physical and geometrical parameters on the streamlines, isotherms, isoconcentrations, average Nusselt and Sherwood numbers has been numerically investigated in detail. The location of heat and solute source has a profound influence on the flow pattern, heat and mass transfer rates in the porous annulus. For the segment located at the bottom portion of inner wall, the flow rate is found to be higher, whereas the heat and mass transfer rates are higher when the source is placed near the middle of the inner wall. Further, the average Sherwood number increases with Lewis number, while for the average Nusselt number the effect is opposite. The average Nusselt number increases with radius ratio (λ); however, the average Sherwood number increases with radius ratio only up to λ = 5, and for λ > 5 , the average Sherwood number does not increase significantly.     

A numerical study based on a weakly compressible formulation for thermosolutal convection in vertical cavities

Laminar thermosolutal convection in cavities with uniform, constant temperature and mass fraction profiles at the vertical side is studied numerically. The study is conducted in the case where an inert carrier gas (species “1”) present in the cavity is not soluble in species “2”, and do not diffuse into the walls. A mass flux of species “2” into the cavity occurs at the hot vertical wall and a mass flux out of the cavity occurs at the opposite cold wall. The weakly compressible model proposed in this work was used to investigate the flow fields, and heat and mass transfer in cavities filled with binary mixtures of ideal gases. The dimensionless form of the seven governing equations for constant thermophysical properties, except density, show that the problem formulation involves ten dimensionless parameters. The results were validated against numerical results published in the literature for purely thermal convection, and thermodynamic predictions for transient thermosolutal flows. A parametric study has been performed to investigate the effects of the initial conditions, molecular weight ratio, Lewis number, and aspect ratio of the cavity for aiding or opposing buoyancy forces. For the range of parameters considered, the results show that variations in the density field have larger effects on mass transfer than on heat transfer. For opposing buoyancy forces, the numerical simulations predict complex flow structures and possible chaotic behavior for rectangular vertical cavities according to the value of the molecular weight ratio.     

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.     

Simultaneous heat and mass transfer in unsteady free convection from two-dimensional and axisymmetric bodies

The unsteady laminar free convection boundary layer flows around two-dimensional and axisymmetric bodies placed in an ambient fluid of infinite extent have been studied when the flow is driven by thermal buoyancy forces and buoyancy forces from species diffusion. The unsteadiness in the flow field is caused by both temperature and concentration at the wall which vary arbitrarily with time. The coupled nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. Computations have been performed for a circular cylinder and a sphere. The skin friction, heat transfer and mass transfer are strongly dependent on the variation of the wall temperature and concentration with time. Also the skin friction and heat transfer increase or decrease as the buoyancy forces from species diffusion assist and oppose, respectively, the thermal buoyancy force, whereas the mass transfer rate is higher for small values of the ratio of the buoyancy parameters than for large values. The local heat and mass transfer rates are maximum at the stagnation point and they decrease progressively with increase of the angular position from the stagnation point.     

Numerical investigation of flow and combined natural-forced convection from an isothermal square cylinder in cross flow

Unsteady flow and heat transfer from a horizontal isothermal square cylinder is studied numerically using a three-dimensional computational model to investigate the influence of buoyancy on the forced flow and heat transfer characteristics. The numerical model is based on a horizontal square cylinder subjected to laminar fluid flow in an unconfined channel. The governing equations in 3D form are solved using a fractional step method based on the finite difference discretization in addition to a Crank–Nicholson scheme employed to the convective and the viscous terms. Two working fluids–air (Pr = 0.7) and water (Pr = 7)–are considered, and the flow and heat transfer simulations were carried out for the Reynolds and Richardson numbers in the intervals 55 ≤ Re ≤ 250 and 0 ≤ Ri ≤ 2, respectively. The flow characteristics such as time-averaged drag/lift, rms drag/rms lift coefficients as well as Strouhal number were computed. The heat transfer from the cylinder is assessed by mean Nusselt number (and rms Nusselt number) over the total heated cylinder walls. As the buoyancy increases, the mass and the velocity of the fluid flowing underneath the cylinder increases. The fluid is injected into the near wake region with an upward motion which significantly alters the flow field in the downstream as well as upstream regions. The effects of Reynolds, Richardson and Prandtl numbers on the flow field and temperature distributions are discussed in detail. It is shown that the flow and heat transfer characteristics are influenced more for air than water. To fill the void in the literature, useful empirical correlations of practical importance are derived for pure forced and pure natural as well as mixed convection. The mixed convection correlations, in terms of the ratio of pure forced convection, are also developed, and their implications are discussed.     

Transient natural convection along vertical cylinder with Heat and Mass transfer

A numerical solution for the transient natural convection flow over a vertical cylinder under the combined buoyancy effect of heat and mass transfer is presented. The velocity, temperature and concentration profiles, local and average skin-friction, Nusselt number and Sherwood number are shown graphically. It is observed that time taken to reach steady state increases with Schmidt number and decreases as combined buoyancy ratio parameter N increases. Stability and convergence of the finite difference scheme are established. Received on 8 July 1997     

Effect of thermo‐solutal stratification on recirculation flow patterns in a backward‐facing step channel flow

This paper presents results on the combined effect of thermo‐solutal buoyancy forces on the recirculatory flow behavior in a horizontal channel with backward‐facing step and the ensuing impact on heat and mass transfer phenomena. The governing equations for double diffusive mixed convection are represented in velocity–vorticity form of momentum equations, velocity Poisson equations, energy and concentration equations. Galerkin's finite‐element method has been employed to solve the governing equations. Recirculatory flow fields with heat and mass transfer are simulated for opposing and aiding thermo‐solutal buoyancy forces by assuming suitable boundary conditions for energy and concentration equations. The effect of Richardson number (0.1?Ri?10) and buoyancy ratio (?10?N?10) on the recirculation bubble and Nusselt and Sherwood numbers are studied in detail. For Richardson number greater than unity, distinct variations in the gradients of Nusselt number and Sherwood number with buoyancy ratio are observed for flow regimes with opposing and aiding buoyancy forces. Copyright © 2009 John Wiley & Sons, Ltd.     

Transient free convection flow between long vertical parallel plates with ramped wall temperature at one boundary in the presence of thermal radiation and constant mass diffusion

The unsteady laminar free convection flow between two long vertical parallel plates with ramped wall temperature at one boundary has been investigated in the presence of thermal radiation and chemical species concentration. The exact solutions of the momentum, energy and concentration equations have been obtained using the Laplace transform technique. The velocity and temperature profiles, skin-friction and Nusselt number variations are shown graphically and the numerical values of the volume flow rate, the total heat rate and species rate added to the fluid are presented in a table. The influence of different system parameters such as the radiation parameter (R), buoyancy ratio parameter (N), Schmidt number (Sc) and time (t) has been analyzed carefully. A critical analysis of the coupled heat and mass transfer phenomena is provided. The free convective flow due to ramped wall temperature has also been compared with the baseline case of flow due to constant wall temperature.     

Coupled heat and mass transfer in mixed convection over a vertical flat plate embedded in saturated porous media: PST/PSC or PHF/PMF

A boundary layer analysis is used to investigate the heat and mass transfer characteristics of mixed convection about a vertical flat plate embedded in a saturated porous medium under the coupled effects of thermal and mass diffusion. The plate is maintained at prescribed surface temperature/concentration (PST/PSC) or prescribed heat/mass flux (PHF/PMF). The nonsimilar governing equations are obtained by using a suitable transformation and solved by Keller box method. Numerical results for the local heat transfer rate and the local mass transfer rate are presented for various parameters. The local heat and mass transfer rates increase with increasing n and m and buoyancy parameter ξ. When buoyancy parameter ξ is very small (large) the value of local Nusselt and the local Sherwood number correspond with the pure forced (free) convection, respectively. Increasing buoyancy ratio N (or N ^*) increases the local heat and mass transfer rates. It is apparent that Lewis number has a pronounced effect on the local mass transfer rate than it does on the local heat transfer rate. Furthermore, increasing Lewis number decreases (increases) the local heat (mass) transfer rate. Received on 8 December 1997     

Double-diffusive natural convection in inclined finned triangular porous enclosures in the presence of heat generation/absorption effects

The problem of double-diffusive convection in inclined finned triangular porous enclosures for various thermal and concentration boundary conditions and in the presence of heat source or sink was studied. The finite difference method was employed to solve the dimensionless governing equations of the problem. The effects of the governing parameters, namely the dimensionless time parameter, the inclination angle, Darcy number, heat generation/absorption parameter, the buoyancy parameter and the Rayleigh number on the streamlines, temperature and concentration contours as well as selected velocity component in the x-direction, local and average Nusselt numbers and local and average Sherwood number at the heated and concentrated wall for various values of the aspect ratio and the position of the fin were considered. The present results are validated by favorable comparisons with previously published results. All the results of the problem were presented in graphical and tabular forms and discussed.     

State space approach to unsteady magnetohydrodynamics natural convection heat and mass transfer through a porous medium saturated with a viscoelastic fluid

A technique of the state space approach and the inversion of the Laplace transform method are applied to dimensionless equations of an unsteady one-dimensional boundary-layer flow due to heat and mass transfer through a porous medium saturated with a viscoelastic fluid bounded by an infinite vertical plate in the presence of a uniform magnetic field is described. Complete analytical solutions for the temperature, concentration, velocity, and induced magnetic and electric fields are presented. The inversion of the Laplace transforms is carried out by using a numerical approach. The proposed method is used to solve two problems: boundary-layer flow in a viscoelastic fluid near a vertical wall subjected to the initial conditions of a stepwise temperature and concentration and viscoelastic fluid flow between two vertical walls. The solutions are found to be dependent on the governing parameters including the Prandtl number, the Schmidt number, the Grashof number, reaction rate coefficient, viscoelastic parameter, and permeability of the porous medium. Effects of these major parameters on the transport behavior are investigated methodically, and typical results are illustrated to reveal the tendency of the solutions. Representative results are presented for the velocity, temperature, concentration, and induced magnetic and electric field distributions, as well as the local skin-friction coefficient and the local Nusselt and Sherwood numbers.     

Thermo-Bioconvection in a Square Porous Cavity Filled by Oxytactic Microorganisms

This paper studies the thermo-bioconvection in a square porous cavity filled by oxytactic microorganisms. The Darcy model with Boussinesq approximation has been used to solve the flow and heat and mass transfer in the porous region. The governing equations formulated in terms of the dimensionless stream function, temperature and concentration have been solved using the finite difference method. Comparison with results from the open literature of the mean Nusselt number for a square cavity filled with a regular porous medium is made. It is shown that the results are in very good agreement. The main objective was to investigate the influence of the traditional Rayleigh number Ra = 10, 100, bioconvection Rayleigh number Rb = 10, 100, Lewis number Le = 1, 10, and Péclet number Pe = 0.1, 1 on the fluid flow and heat and mass transfer. Comprehensive analysis of an effect of these key parameters on the Nusselt and Sherwood numbers at the vertical walls has been conducted.     

Combined heat and mass transfer by laminar natural convection from a vertical plate

Simultaneous heat and mass transfer in buoyancy-induced laminar boundary-layer flow along a vertical plate is studied for any ratio of the solutal buoyancy force to the thermal buoyancy force by using a new similarity transformation. The effects of the buoyancy ratio and Lewis number on the rates of heat and mass transfer are presented explicitly for most practical gaseous solutions (Pr=0.7, 0.21≤Sc≤2.1) and aqueous solutions (Pr=7, 140≤Sc≤1400). Very accurate correlations of the mass transfer and heat transfer rates are developed for the cases of single and combined buoyancy forces.     

Fully Developed Forced Convection Heat Transfer in a Porous Channel with Asymmetric Heat Flux Boundary Conditions

An analytical study is performed on steady, laminar, and fully developed forced convection heat transfer in a parallel plate channel with asymmetric uniform heat flux boundary conditions. The channel is filled with a saturated porous medium, and the lower and upper walls are subjected to different uniform heat fluxes. The dimensionless form of the Darcy–Brinkman momentum equation is solved to determine the dimensionless velocity profile, while the dimensionless energy equation is solved to obtain temperature profile for a hydrodynamically and thermally fully developed flow in the channel. Nusselt numbers for the lower and upper walls and an overall Nusselt number are defined. Analytical expressions for determination of the Nusselt numbers and critical heat flux ratio, at which singularities are observed for individual Nusselt numbers, are obtained. Based on the values of critical heat flux ratio and Darcy number, a diagram is provided to determine the direction of heat transfer between the lower or upper walls while the fluid is flowing in the channel.     

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.     

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.     

Local nonsimilarity solution for the impact of the buoyancy force on heat and mass transfer in a flow over a porous wedge with a heat source in the presence of suction/injection

Combined heat and mass transfer in free, forced and mixed convection flows along a porous wedge with internal heat generation in the presence of uniform suction or injection is investigated. The boundary-layer analysis is formulated in terms of the combined thermal and solute buoyancy effect. The flow field characteristics are analyzed using the Runge-Kutta-Gill method, the shooting method, and the local nonsimilarity method. Due to the effect of the buoyancy force, power law of temperature and concentration, and suction/injection on the wall of the wedge, the flow field is locally nonsimilar. Numerical calculations up to third-order level of truncation are carried out for different values of dimensionless parameters as a special case. The effects of the buoyancy force, suction, heat generation, and variable wall temperature and concentration on the dimensionless velocity, temperature, and concentration profiles are studied. The results obtained are found to be in good agreement with previously published works.     

Unsteady Natural Convection,Heat and Mass Transfer in Inclined Triangular Porous Enclosures in the Presence of Heat Source or Sink: Effect of Sinusoidal Variation of Boundary Conditions

The problem of double-diffusive convection in inclined triangular porous enclosures with sinusoidal variation of boundary conditions in the presence of heat source or sink was discussed numerically. The dimensionless governing equations of the problem were solved numerically by using finite difference method. The effects of governing parameters, namely, the dimensionless time parameter, various values of the inclination angle, Darcy number, the heat generation/absorption parameter, the buoyancy parameter and the amplitude wave length ratio on the streamlines, temperature and concentration contours as well as the velocity component in the x-direction at the triangle mid-section, the average Nusselt and Sherwood numbers at the bottom wall of the triangle for various values of aspect ratio were considered. The present results are validated by favorable comparisons with previously published results. All the results of the problem were presented in graphical and tabular forms and discussed.