Thermal radiation effect on mixed convection from horizontal surfaces in saturated porous media

An analysis is presented to investigate the effect of radiation on mixed convection from a horizontal flat plate in a saturated porous medium. Both a hot surface facing upward and a cold surface facing downward are considered in the analysis. The conservation equations that govern the problem are reduced to a system of nonlinear ordinary different equations. The important parameters of this problem are the radiation parameter R, the buoyancy parameter B, and the freestream to wall temperature ratio T _∞/T _w for the case of a hot surface or the wall to freestream to wall temperature T _w /T _∞ for the case of a cold surface.     

Inertia effects on mixed convection for vertical plates with variable wall temperature in saturated porous media

Mixed convection in a porous medium from a vertical plate with variable wall temperature is investigated. The entire mixed convection regime is divided into two regions and two sets of transformations are used. The first region is for the forced convection dominated regime (FCDR) and the other one is for the natural convection dominated regime (NCDR). The dimensionless parameterK′ U _∞/ν is found to characterize the effect of inertia resistance in the first region and the dimensionless parameter (K′ U _∞/ν (Ra _x /Pe _x ) is found to characterize the effect of inertia resistance in the second region. The solution of the first region is carried out from ξ _f =0.0 to ξ _f =1.0, and the solution of the second region is carried out from ξ _n =0.0 to ξ _n =1.0. Velocity and temperature profiles are calculated at different values of the governing parameters. Local Nusslet number variation for the entire mixed convection regime is also calculated and presented.     

Non-darcy natural convection on a vertical cylinder in a saturated porous medium

The steady free convection boundary layer flow of non-Darcy fluid along an isothermal vertical cylinder embedded in a saturated porous medium using the Ergun model has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme developed by Keller. It is found that the heat transfer is strongly affected by the modified Grashof number which characterizes the non-Darcy fluid, and the curvature parameter. Also the heat transfer is found to be more than that of the flat plate.     

Non-Darcy double-diffusive natural convection in axisymmetric fluid saturated porous cavities

Double-diffusive natural convection in a fluid saturated porous medium has been investigated using the finite element method. A generalised porous medium model is used to study both Darcy and non-Darcy flow regimes in an axisymmetric cavity. Results indicate that the Darcy number should be a separate parameter to understand flow characteristics in non-Darcy regime. The influence of porosity on heat and mass transfer is significant and the transport rates may differ by 25% or more, at higher Darcy and Rayleigh numbers. When compared with the Darcy and other specialised models of Brinkman and Forchheimer, the present generalised model predicts the least heat and mass transfer rates. It is also observed that an increase in radius ratio leads to higher Nusselt and Sherwood numbers along the inner wall.     

No slip boundary effects in non-Darcian mixed convection from a vertical wall in saturated porous media. Analytical solution

An analytical study is made for wall effects in non-Darcy mixed convection from vertical impermeable surfaces embedded in a saturated porous medium. The governing equations are transformed into a dimensionless form by non-similar transformation to cover both forced and natural convection dominated regimes. Two different dimensionless parameters that measure the strength of mixed convection were found in both regimes. The parameters of forced convection dominated regime can be related to those of natural convection dominated regime. An approximate analytical solution for the governing equations was obtained. Temperature and velocity profiles for both regimes are presented. Received on 9 September 1997     

Non-Darcy mixed convection from a vertical plate in saturated porous media-variable surface heat flux

Nonsimilarity solutions for non-Darcy mixed convection from a vertical impermeable surface embedded in a saturated porous medium are presented for variable surface heat flux (VHF) of the power-law form. The entire mixed convection region is divided into two regimes. One region covers the forced convection dominated regime and the other one covers the natural convection dominated regime. The governing equations are first transformed into a dimensionless form by the nonsimilar transformation and then solved by a finite-difference scheme. Computations are based on Keller Box method and a tolerance of iteration of 10⁻⁵ as a criterion for convergence. Three physical aspects are introduced. One measures the strength of mixed convection where the dimensionless parameter Ra^* _x /Pe^3/2 _x characterizes the effect of buoyancy forces on the forced convection; while the parameter Pe _x /Ra^*2/3 _x characterizes the effect of forced flow on the natural convection. The second aspect represents the effect of the inertial resistance where the parameter K′U _∞/ν is found to characterize the effect of inertial force in the forced convection dominated regime, while the parameter (K′U _∞/ν)(Ra^*2/3 _x /Pe _x ) characterizes the effect of inertial force in the natural convection dominated regime. The third aspect is the effect of the heating condition at the wall on the mixed convection, which is presented by m, the power index of the power-law form heating condition. Numerical results for both heating conditions are carried out. Distributions of dimensionless temperature and velocity profiles for both Darcy and non-Darcy models are presented. Received on 26 May 1997     

Comprehensive correlations for laminar mixed convection on vertical and horizontal flat plates

This paper presents comprehensive correlation equations of the local Nusselt numbers and surface shear stresses for laminar forced convection, natural convection, and mixed convection on vertical and horizontal flat plates which are maintained with uniform wall temperature or uniform surface heat flux. The correlation for pure forced convection and pure natural convection are very accurate for any Prandtl number between 0.001 and infinity. The correlations for mixed convection coincide very well with the numerical results over the entire regimes of mixed convection intensity and Prandtl number for the eight cases of the two plates with distinct thermal boundary conditions and buoyancy-assisting and-opposing flow conditions.     

Methodical investigation of free convection from vertical and horizontal plates

This paper attempts of review the literature on free convection from flat surfaces transferring heat in an unlimited space. Data from the references have been presented as final criterial relations both in tabular form and in graphic system, separately for the vertical plates and the horizontal ones. About 45 results of investigations on free convection have been examined and the following mean criterial relations have been calculated:Nu=0.550·(Ra)^0.250 for vertical plates,Nu=0.310·(Ra)^0.290 for horizontal plates. It has been observed for the vertical plates that divergence of these results relatively to the mean values is constant within the whole variability range of the Rayleigh number and amounts to about ± 15%. In the case of horizontal plates the divergencies are more considerable and vary within ± 50% for the laminar, and within ± 25% for the turbulent range. In the next step the attempt was made to determine the causes of these divergencies and errors. The work was performed because of two main reasons: imperfection of measuring methods and a lack of uniformity in the choice of the characteristic linear dimension.     

Darcian mixed convection plumes along vertical adiabatic surfaces in a saturated porous medium

The mixed convection flow due to a line thermal source embedded at the leading edge of an adiabatic vertical plane surface immersed in a saturated porous medium has been studied. Both weakly and strongly buoyant plume regimes have been considered. The cases of buoyancy assisting and buoyancy opposing flow conditions have been incorporated in the analysis. The results are presented for the entire range of buoyancy parameter from the pure forced convection (ξ=0) to the pure free convection (ξ → ∞@#@) regimes. For buoyancy-assisting flow, the wall temperature and the velocity at the wall increase as the plume strength increases. However, they all decrease as the free-stream velocity increases. For buoyancyopposing flow, the temperature at the wall increases as the strength of the plume increases but velocity at the wall decreases.     

Transient non-Darcy free convection between parallel vertical plates in a fluid saturated porous medium

Transient non-Darcy free convection between two parallel vertical plates in a fluid saturated porous medium is investigated using the generalized momentum equation proposed by Vafai and Tien. The effects of porous inertia and solid boundary are considered in addition to the Darcy flow resistance. Exact solutions are found for the asymptotic states at small and large times. The large time solutions reveal that the velocity profiles are rather sensitive to the Darcy number Da when Da<1. It has also been found that boundary friction alters the velocity distribution near the wall, considerably. Finite difference calculations have also been carried out to investigate the transient behaviour at the intermediate times in which no similarity solutions are possible. This analytical and numerical study reveals that the transient free convection between the parallel plates may well be described by matching the two distinct asymptotic solutions obtained at small and large times.Nomenclature C empirical constant for the Forchheimer term - f velocity function for the small time solution - F velocity function for the large time solution - g acceleration due to gravity - Gr* micro-scale Grashof number - H a half distance between two infinite plates - K permeability - Nu Nusselt number - Pr Prandtl number - t time - T temperature - u, v Darcian velocity components - x, y Cartesian coordinates - effective thermal diffusivity - coefficient of thermal expansion - porosity - dimensionless time - similarity variable - dimensionless temperature - viscosity - kinematic viscosity - density - the ratio of heat capacities     

Heat dispersion in stationary mixed convection flow about horizontal surfaces in porous media

An analysis has been performed to study the influence of velocity dependent dispersion on transverse heat transfer in mixed convection flow above a horizontal wall of prescribed temperature in a saturated porous medium. The Boussinesq approximation and boundary layer analysis were used to numerically obtain gravity affected temperature and velocity distributions within the frames of Darcy's law and a total thermal diffusivity tensor comprising both of constant coefficient heat conduction and velocity proportional mechanical heat dispersion. Dependending on Pe_∞, the molecular Peclét number basing on the effective thermal diffusivity and the velocity of the oncoming flow, density coupling has distinct influences on heat transfer rates between the wall surface and the porous medium flow region. For small Peclét numbers, when heat conduction is the prevailing mechanism, wall heat fluxes are the higher the larger the density difference between the oncoming and the near wall fluid is. The opposite is true for larger Peclét numbers, when mechanical heat dispersion is the main cause of heat spreading. For Pe_∞ tending to infinity these wall heat fluxes approach finite maximum values in the total heat diffusivity model, they grow beyond any limit if only constant coefficient heat conduction is considered. Thus, the inclusion of mechanical heat dispersion effects yields physically more realistic predictions. Received on 18 September 1996     

Onset of double-diffusive convection in a saturated porous layer with time-periodic surface heating

The stability of a fluid saturated, horizontal porous layer in the presence of a solute concentration gradient and time-periodic thermal gradient is examined. The modulated gradient is the result of a sinusoidal upper surface temperature which models the effect of variable solar radiation heating of the layer. Darcy's law and the Boussinesq approximation are employed, and we assume an equation of state linear in temperature and concentration. A linear stability analysis is carried out to obtain predictions for the onset of convection and critical wavenumbers for the system. The critical conditions are obtained via the Galerkin method and Floquet theory. The effects of variable concentration gradient, temperature modulation amplitude and frequency are examined, and compared with the results obtained analytically from the corresponding unmodulated problem. It is shown that instabilities can occur as convective motions which are synchronous or subharmonic with the surface heating, or can be identified via complex conjugate Floquet exponents. The neutral stability curves at the transitions between instabilities are found to be bimodal when the temperature is time-periodic, and are characterized by jumps in the critical wavenumbers. Received February 5, 1998     

Laminar natural convection between partially heated vertical parallel plates

A numerical analysis has been performed to examine the characteristics of laminar natural convection in vertical channel with unheated entry and unheated exit. The heated section is subjected to uniform wall temperature (UWT) or uniform heat flux (UHF). Theoretical results for average Nusselt number and induced volume flow rateQ were derived under fully developed condition. Particular attention is paid to investigating the effects of the partially heated section on the induced volume flow rate and Nusselt number for various conditions. Results show that for UWT the induced volume flow rateQ increases with decreasing unheated entry length or increasing total length of heated section and unheated exit. For a fixed unheated entry length, the channel with a longer heated section length causes a greaterQ. Additionally, for both UWT and UHF, the average Nusselt number under fully-developed condition increases with increasing value ofE ₁/E ₂.

---

Laminare Naturkonvektion zwischen teilbeheizten senkrechten Parallelplatten

Zusammenfassung Die numerische Untersuchung soll das Verhalten einer laminaren Strömung bei natürlicher Konvektion in einem senkrechten Kanal mit unbeheiztem Ein- und Austritt klären. Der beheizte Abschnitt wird entweder mit gleichförmiger Wandtemperatur (UWT) oder gleichförmigem Wärmefluß (UHF) beaufschlagt. Bezüglich voll ausgebildeter Strömung ließen sich theoretische Ergebnisse für die mittlere Nußelt-Zahl und den induzierten VolumenstromQ gewinnen. Besonderes Interesse galt der Untersuchung des Einflusses des teilbeaufschlagten Abschnittes auf Volumenstrom und Nußelt-Zahl unter verschiedenen Nebenbedingungen. Die Ergebnisse zeigen, daß im UWT-Fall der induzierte VolumenstromQ mit abnehmender unbeheizter Einlauflänge oder zunehmender Gesamtlänge des Heizabschnittes anwächst. Bei fester unbeheizter Einlauflänge erzeugt der Kanal mit längerem Heizabschnitt einen höheren StromQ.

---

Nomenclature b half channel width - E ₁ ratio of unheated exit length to channel length,l ₁/l - E ₂ ratio of heated section length to channel length,l ₂/l - E ₃ ratio of unheated entry length to channel length,l ₃/l - Gr _L Grashof number, Eq. (4) - g gravitational acceleration - k thermal conductivity - l channel length - l ₁ unheated exit length - l ₂ heated section length - l ₃ unheated entry length - L dimensionless channel length, Eq. (4) - L ₁ dimensionless unheated exit length - L ₂ dimensionless heated section length - L ₃ dimensionless unheated entry length - average Nusselt number - Nu _x local Nusselt number - p pressure defect - P dimensionless pressure defect, Eq. (4) - Pr Prandtl number,/ - q _w wall heat flux of heated section - Q dimensionless induced volume flow rate, Eq. (8) - Ra _E Rayleigh number based on the heated section length,Ra _L/E₂ - Ra _L Rayleigh number based on the full channel length,Gr _LPr - T temperature - T ₀ inlet temperature - T _w wall temperature - u, v velocity components in thex andy directions, respectively - U, V dimensionless velocity components in thex andy directions, respectively, Eq. (4) - u ₀,U ₀ dimensional and dimensionless inlet velocity, respectively - x, y coordinates in thex andy directions, respectively - X, Y dimensionless coordinates in thex andy directions, respectively, Eq. (4) - X ratio of longitudinal distance from the entrance of heated section to the heated section length,X=[X–(L–L ₁–L ₂)]/L ₂ Greek symbols thermal diffusivity - thermal expansion coefficient - kinematic viscosity - dimensionless temperature, Eq. (4) - ₀ fluid density at ambient temperature     

Free-forced convection from a heated longitudinal horizontal cylinder embedded in a porous medium

The flow and heat transfer from a heated semi-infinite horizontal circular cylinder which is moving with a constant speed into a porous medium is considered. It is assumed that the Grashof and Reynolds numbers are large so that the governing equations are the three dimensional boundary-layer equations. A numerical procedure for solving these equations is described and the asymptotic solutions which are valid both near and distant from the leading edge of the cylinder are presented. The range of validity of these asymptotic solutions is discussed and the results are compared in detail with the full numerical solution. The problem is of practical importance, for example in the drilling of pipes into a geothermal reservoir.

---

Freie erzwungene Konvektion von einem beheizten schlanken horizontalen Zylinder, eingebettet in ein poröses Medium

Zusammenfassung Es wird die Strömung und der Wärmeübergang an einem beheizten, halbunendlichen horizontalen Kreiszylinder betrachtet, der mit konstanter Geschwindigkeit sich in ein poröses Medium bewegt. Dabei wird angenommen, daß die Grashof- und Reynolds-Zahlen groß sind, so daß die Bestimmungsgleichungen von den dreidimensionalen Grenzschichtgleichungen gebildet werden. Es wird ein numerisches Verfahren zur Lösung dieser Gleichungen beschrieben und eine asymptotische Lösung präsentiert, die sowohl in der Nähe als auch in großem Abstand von dem vorderen Ende des Zylinders gültig ist. Der Gültigkeitsbereich dieser asymptotischen Lösungen wird diskutiert und die Ergebnisse werden im Detail mit vollständigen numerischen Lösungen verglichen. Das Problem ist z.B. beim Eindringen von Rohrleitungen in geothermische Reservoire von praktischer Wichtigkeit.

---

Nomenclature a radius of cylinder - Gr Grashof number (=g(T_w-Ta/²) - g acceleration due to gravity - permeability in the porous medium - Nu local Nusselt number - total heat flux from cylinder - q _w heat flux from cylinder - r radial co ordinate - Ra Rayleigh number (=g (T_w - T_t8) a/ ) - Re Reynolds number (=U _t8 a/) - T temperature - u, v, w speeds inx, , r directions - x axial co ordinate - equivalent thermal diffusivity - thermal expansion coefficient - ratioGr/Re - similarity variable - dimensionless temperature (=(T- T)/(T _w- T) - kinematic viscosity - azimuthal co ordinate - w cylinder surface - free stream     

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

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

Free convection from horizontal screened plates

Experimental investigations of laminar free convection from horizontal isothermal surfaces screened by cylindrical vertical walls are presented. Screen diameters (D) were equal to the diameter of a heating plate while their heights (H) were varied. Results obtained for several heating fluxes and two different liquids indicate that depending on the (H/D) ratio of cylindrical screens, the heat flux transferred from the heating plate to the liquid varies. Three subranges of screening effect have been distinguished: primary inhibition, intensification and secondary inhibition. For primary inhibition, a theoretical model has been proposed and critical ratio of (H/D) has been predicted. The results of flow visualization are also discussed.

---

Freie Konvektion von horizontalen, seitlich abgeschirmten Platten

Zusammenfassung Es werden experimentelle Untersuchungen der laminaren freien Konvektion an horizontalen, isothermen Flächen mit vertikalen zylindrischen Abschirmflächen beschrieben. Die Schirmdurchmesser (D) waren gleich dem einer Heizplatte, während die Schirmhöhen (H) variiert wurden. Die Resultate, welche für mehrere Wärmeflüsse und mit zwei verschiedenen Flüssigkeiten gewonnen wurden, zeigen, daß — abhängig von VerhältnisH/D der zylindrischen Schirme — der von der Heizplatte an das Fluid übergehende Wärmestrom stark variiert. Drei Bereiche konnten bezüglich des Abschirmeffektes unterschieden werden: erste Hemmung, Verstärkung und zweite Hemmung. Für den Fall der ersten Hemmung wird ein theoretisches Modell vorgeschlagen und das kritischeH/D-Verhältnis vorausberechnet. Abschließend erfolgt die Diskussion von Ergebnissen, welche aus der sichtbaren Darstellung des Strömungsfeldes gewonnen wurde.

---

Nomenclature a thermal diffusivity [m²/s] - D external diameter of the plate, internal diameter of screens [m] - g gravitational acceleration [m/s²] - H height of screen [m] - Nu= D/ Nusselt number - Nu _(H/D)=0 Nusselt number for plate without screens - Nu _H/D Nusselt number for plate with screen of heightH - Ra=g (T _w1–T _F)D ³/( a) Rayleigh number - Ra _H=g (T _w–T _w1)H ³/( a) Rayleigh number inside screen - q heat flux density [W/m²] - T _w surface temperature of the heating plate [°C] - T _w1 boundary temperature of the motionless liquid [°C] - T _F bulk temperature of the liquid [°C] - heat transfer coefficient [W/(m²·K)] - coefficient of volumetric expansion [1/K] - kinematic viscosity - thermal conductivity [W/(m·K)]     

Nonsimilar solutions for free convection from a vertical plate in porous media

A nonsimilar boundary layer analysis has been presented for the free convection along a vertical plate embedded in a fluid-saturated porous medium in the presence of surface mass transfer and internal heat generation. The transformed conservation laws are solved numerically for the cases of variable wall temperature and variable wall heat flux boundary conditions. Results are presented for the details of the velocity and temperature fields as well as Nusselt number. Received on 13 December 1996     

A correlation for laminar mixed convection from vertical plates using transient experiments

The applicability of the transient cooling experimental technique is studied for developing a correlation for Nusselt number with other pertinent parameters for laminar aiding mixed convection flows over an isothermal vertical plate. A heated aluminum plate, modeled as a lumped system, is allowed to cool in the mixed convection environment during which the transient response of the plate is recorded using a data acquisition system. The experimentally known transient response of the plate is then compared with the numerically computed transient response to estimate the square of the residual. The minimization of this residual using Levenberg–Marquardt iterative procedure retrieves the values of relevant parameters in the correlation, whose form is assumed a priori. The experiments were conducted for a temperature range of 410–305 K and a corresponding Richardson number range 20–0.04. The retrieved values of the parameters compare well with those available in literature.     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

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.