Natural convection in inclined rectangular enclosures with perfectly conducting fins attached on the heated wall

The natural convection heat transfer in inclined rectangular enclosures with perfectly conducting fins attached to the heated wall is numerically studied. The parameters governing this problem are the Rayleigh number (10²≤Ra≤2×10⁵), the aspect ratio of the enclosures (2.5≤A=H′/L′≤∞), the dimensionless lengths of the partitions (0≤B=?′/L′≤1), the aspect ratio of micro-cavities (A≤C=h′/L′≤0.33), the inclination angle (0≤φ≤60^°) and the Prandtl number (Pr=0.72). The results indicate that the heat transfer through the cover is considerably affected by the presence of the fins. At low Rayleigh numbers, the heat transfer regime is dominated by conduction. When B≈0.75 and C≈0.33, the heat transfer through the cold wall decreases considerably. This trend is enhanced when the enclosure is inclined. Useful engineering correlations are derived for practical applications.     

Mixed convection heat transfer from a horizontal channel with protruding heat sources

A numerical investigation is carried out to study fluid flow and heat transfer characteristics of conjugate mixed convection from a two dimensional horizontal channel with four protruding heat sources mounted on one of the finite thick channel walls. The flow is assumed as laminar, hydrodynamically and thermally developing. Water and FC70 are the fluids under consideration. The geometric parameters such as spacing between the channel walls (S), size of protruding heat sources (L_h×t_h), thickness of substrate (t) and spacing between heat sources (b) are fixed. Results are presented to show the effect of parameters such as Re_S, Gr_S^*, Pr, k_p/k_f and k_s/k_f on fluid flow and heat transfer characteristics. Using the method of asymptotic expansions, correlations are also presented for the maximum temperature of heat source.     

Mixed convection and conduction heat transfer in open cavities

A numerical study has been carried out on mixed convection and conduction in open cavities. The study covers the Rayleigh number from 10 to 10⁶, the Reynolds number from 1 to 10³ andk _r =k _w /k _a from 1 to 100 forPr=0.72 (air) and cavity aspect ratioB from 0.8 to 1.3. The flow is assumed to be laminar and two-dimensional. The density variation is taken into account by the Boussinesq approximation. The controlvolume approach is used for solving the governing equations of conjugate heat transfer involving conduction in the walls. Streamlines and isotherms in the system are produced and the results are represented in terms of the Nusselt number as function of other parameters. Correlations are derived to calculate heat transfer through the cavity openings.     

Unsteady entrance heat transfer in laminar pipe flows with convection from the ambient

The characteristics of unsteady entrance heat transfer in the combined entrance heat transfer region of laminar pipe flows resulting from time-varying inlet temperature are numerically investigated. Three non-dimensional parameters,Nu ₀, a^*, andf are identified in the study. Also, their effects on the non-dimensional duct wall temperature, fluid bulk temperature, and duct wall heat flux are discussed in great detail. Comparisons are made with the zero thermal capacity wall solution.     

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     

Soret and Dufour effects on Double-Diffusive Free Convection From a Corrugated Vertical Surface in a Non-Darcy Porous Medium

Combined heat and mass transfer process by natural convection from a wavy vertical surface immersed in a fluid-saturated semi-infinite porous medium due to Soret and Dufour effects for Forchheimer extended non-Darcy model has been analyzed. A similarity transformation followed by a wavy to flat surface transformation is applied to the governing coupled non-linear partial differential equations, and they are reduced to boundary layer equations. The obtained boundary layer equations are solved by finite difference scheme based on the Keller-Box approach in conjunction with block-tridiagonal solver. Detailed simulations are carried out for a wide range of parameters like Groshof number (Gr*), Lewis number (Le), Buoyancy ratio (B), Wavy wall amplitude (a), Soret number (S _r ), and Dufour number (D _f ). Comparison tables local and average Nusselt (Nu) number, local and average Sherwood (Sh) number plots are presented.     

Thermal protection characteristics of a vertical rectangular cell filled with PCM/air layer

The present paper presents a numerical analysis concerning thermal protection characteristics of a vertical rectangular composite cell filled with a solid-liquid phase change material (PCM) and air layer. Inside the composite cell the PCM layer is separated from air layer by a solid partition of negligible thickness. The buoyancy-induced flows developed in both the air-filled layer and the molten PCM zone inside the PCM layer were modeled as two-dimensional laminar Newtonian fluid flow adhering to the Boussinesq approximation. Meanwhile, two-dimensional conduction heat transfer was accounted for the unmelted solid PCM region. Delineation is made via a parametric simulation of the effects of the pertinent parameters:Ste (Stefan number),Sc (subcooling factor),Ra (Rayleigh number), aspect ratio of composite cell,A, and relative thickness ratioA _p /A _a , on the transient thermal protection performance of the composite cell. Results demonstrate that by means of the latent-heat absorption inside the PCM layer, heat penetration across the composite cell can be greatly retarded over an effective duration until a critical instant until the melting front of PCM reaches the partition wall. Such an effective thermal protection duration is found to be a strong function ofRa, Ste, A _p /A _a , andA. In addition, the results of the transient heat transfer rate penetrating through the composite cell are examined as a function of the pertinent parameters of the problem.     

Transient behavior of heat removal from a cylindrical heat storage vessel packed with spherical porous particles

The transient heat transfer behavior in the case of heat removal from a cylindrical heat storage vessel packed with spherical particles was investigated experimentally for various factors (flow rate, diameter of spherical particles packed, temperature difference between flowing cold air and spherical particles accumulating heat, and physical properties of spherical particles). The experiments were covered in ranges of Reynolds number based on the mean diameter of spherical particles packed Re_d = 10.3–2200, porosity?=0.310 to 0.475, ratio of spherical particle diameter to cylinder diameterd/D = 0.0075–0.177 and ratio of length of the cylinder to cylinder diameterL/D=2.5–10. It was found that especially the flow rate and the dimension of spherical particles played an important role in estimating the transient local heat transfer characteristics near the wall of the cylindrical vessel in the present heat storage system. As flow rate and diameter of spherical particles were increased under a given diameter of the cylinder heat storage vessel, the mean heat transfer coefficient between the flow cold air and the hot spherical particles increased and the time period to finish removing heat from the vessel reduced. In addition, the useful experimental correlation equations of mean heat transfer coefficient between both phases and the time period to finish removing heat from the vessel were derived with the functional relationship of Nusselt numberNu _d=f [modified Prandtl numberPr ^* (d/D), Re_d) and Fourier numberFo = f(d/D, L/D, Pr^*, Re_d).     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.     

Heat and mass transfer characteristics of the self-similar boundary-layer flows induced by continuous surfaces stretched with rapidly decreasing velocities

The mechanical and thermal characteristics of the self-similar boundary-layer flows induced by continuous surfaces stretched with rapidly decreasing power-law velocities U _w∝x ^m , m<?1 are considered. Comparing to the well studied cases of the increasing stretching velocities (m>0) several new features of basic significance have been found. Thus: (i) for m<?1 the boundary layer equations admit self-similar solutions only if a lateral suction is applied; (ii) the dimensionless suction velocity f _w<0 must be strong enough, i.e. f _w<f _w,max(m) where f _w,max(m) depends on m so that its absolute maximum max (f _w,max(m))=?2.279 is reached for m→?∞, while for m→?1, f _w,max(m)→?∞; (iii) the case {m→?∞, f _w,max(m)=?2.279} of the flow boundary value problem is isomorphic to the stretching problems with exponentially decreasing velocities U _w∝e ^ax with arbitrary a<0; (iv) for any fixed m<?1 and f _w<f _w,max(m) the flow problem admits a non-denumerable infinity of multiple solutions corresponding to the values of the dimensionless skin friction f ^″(0)≡s belonging to a finite interval s∈ [s _min(f _w,m), s _max(f _w,m)]; (v) the solution is only unique for f _w=f _w,max(m) where s=s _min(f _w,m)= s _max(f _w,m) holds; (vi) to every one of the multiple solutions of the flow problem there corresponds a unique solution of the heat transfer problem with a wall temperature distribution T _w?T _∞∝x ⁿ and a well defined and distinct value of the dimensionless wall temperature gradient ?^′(0), except for the cases n=(|m|?1)/2 where ?^′(0) has the same value ?^′(0)=Pr·f _w for the whole class of flow solutions with s∈[s _min(f _w,m), s _max(f _w,m)]; (vii) for f _w→?∞ one obtains the `asymptotic suction profiles' corresponding to s=s _min(f _w,m)?f _w and ?^′(0)?Pr·f _w in an explicit analytic form. The paper includes several examples which illustrate the dependence of the heat and fluid flows induced by surfaces stretching with rapidly decreasing velocities on the physical parameters f _w, m, n and Pr.     

Natural convection heat transfer in enclosures with conducting multiple partitions and side walls

The flow and heat transfer in enclosures with conducting multiple partitions and side walls were numerically analyzed. Side walls were kept at isothermal conditions, while top and bottom walls were insulated. Employing control volume approach, a computer program based on SIMPLE algorithm was developed. Computations were carried out to investigate the effects of Rayleigh number, number of partitions and cavity aspect ratios on the heat transfer rate. The mean Nusselt numbers were calculated from computed temperature fields. It was observed that, the mean Nusselt number decreases with increasing partition number. It is inversely proportional to (1+N) for N≤4. For all partition numbers, the mean Nusselt number increases with increasing Rayleigh number. On the other hand, the cavity aspect ratio does not affect the mean Nusselt number to a considerable extent for considered aspect ratios in this study.     

Heat transfer and skin-friction in a turbulent boundary layer under a non-equilibrium longitudinal adverse pressure gradient

The experimental data on the effect of weak and moderate non-equilibrium adverse pressure gradients (APG) on the parameters of dynamic and thermal boundary layers are presented. The Reynolds number based on the momentum thickness at the beginning of the APG region was Re^** = 5500. The APG region was a slot channel with upper wall expansion angles from 0 to 14°. The profiles of the mean and fluctuation velocity components were measured using a single-component hot-wire anemometer. The friction coefficients were determined using two methods, namely, the indirect Clauser method and the direct method of weighting the lower wall region on a single-component strain-gage balance. The heat transfer coefficients were determined by a transient method using an IR camera. It is noticed that in the pressure gradient range realized the universal logarithmic region in the boundary layer profile is conserved. The values of the relative (divided by the parameters in zero gradient flow at the same value of Re^**) friction and heat transfer coefficients, together with the Reynolds analogy factor, are determined as functions of the longitudinal pressure gradient. The values of the relative friction coefficient reduced to c_f/c_f0 = 0.7 and those of the heat transfer to St/St₀ = 0.9. A maximum value of the Reynolds analogy factor (St/St₀)/(c_f/c_f0) = 1.16 was reached for the pressure gradient parameter β = 2.9.     

Natural convection in shallow enclosures with multiple conducting partitions

Laminar natural convection and conduction in shallow enclosures having multiple partitions with finite thickness and conductivity have been studied. An approximate analytical solution is obtained by using the parallel flow approximation in horizontal shallow enclosures heated isothermally at two vertical ends while adiabatic on horizontal end walls. The same problem is solved also using a finite difference formulation and the control volume method. The study covers the range ofRa from 105 to 107,A=H/L0.2, C=1/L from 0 to 0.15, and the thermal conductivity ratio of partition to fluidk _r from 10^–4 to 10¹¹. The partition numberN was varied from 0 to 5. The Prandtl number was 0.72 (for air). The results are reduced in terms ofNu as a function ofRa, k, and various geometrical parameters (A, C). The streamlines and isotherms are produced to visualize the flow and temperature fields.Es wird der kombinierte Einfluß von laminarer Naturkonvektion und Leitung in flachen Behältern mit mehreren Trennwänden endlicher Dicke und Leitfähigkeit untersucht. Eine analytische Näherungslösung läßt sich über die Parallelstromapproximation bezüglich horizontaler flacher Behälter finden, deren zwei vertikale Begrenzungswände isotherm beheizt sind, während die Horizontalflächen adiabat sein sollen. Das selbe Problem wird unter Verwendung eines Differenzverfahrens und der Kontrollvolumen-Methode gelöst und zwar für die Parameterbereiche 105 Ra 10⁷;A=H/L<0.2;>C=1/L 0.15; 10^–4k_r 1011, wobei der letzte Parameter das Verhältnis der Leitfähigkeit von Trennwand und Fluid bezeichnet. Die Zahl der TrennwändeN variierte Zwischen 0 und 5, die Prandtl-Zahl betrug 0.72 (Luft). Die Ergebnisse werden in dimensionsloser Form gemäß der BeziehungNu =f (Ra, k _r ,A, C) mitgeteilt bzw. durch Diagrammdarstellungen der Stromlinien- und Isothermenfelder veranschaulicht.Financial support from the Natural Sciences and Engineering Council Canada is acknowledged. Financial support to A. Kangni from Canadian Fellowship Program For French Speaking Countries is also acknowledged.     

The comparison model method: A new arithmetic approach to the discrete inverse problem of groundwater hydrology

Let the steady-state pressure z(·) of a fluid in a one-dimensional domain be governed by the equation d _x (a d _x z) = f subject to Dirichlet boundary conditions. We consider the identification of the transmissivity a (·), given f(·), and measured pressure z(·) by the comparison model method, a direct method which has been known and applied for some time but lacked theoretical background. With reference to a domain in one spatial dimension, we examine both the infinite-(‘continuous’) and finite-(discrete) dimensional cases. In the former, the method is based on the solution p(·) of an auxiliary flow equation, where f(·) and the two-point boundary conditions are unchanged with respect to the original or z(·) equation, whereas a tentative constant value b is assigned to the auxiliary transmissivity. The ratio of the first derivatives of p(·) and z(·) multiplied by b yields a solution ã(·) to the inverse problem. We examine in detail the nonuniqueness of ã(·) as a function of b, some cases where existence implies uniqueness, the role of positivity constraints, and a special feature: self-identifiability. We then translate all available results into the discrete case, where the good unknowns for the inverse problem are the internode coefficients. Several algebraic and numerical examples are presented.     

Homogenization of a degenerate parabolic problem in a highly heteregeneous mediumHomogénéisation d'un problème parabolique dégénéré dans un milieu fortement hétérogène

We consider the homogenization of a time-dependent heat transfer problem in a highly heteregeneous periodic medium made of two connected components having finite heat capacities c_α(x) and heat conductivities a_α(x), α=1,2, of order one, separated by a third material with thickness of order ε the size of the basic periodicity cell, but with conductivity λa₃(x) where a₃=O(1) and λ tends to zero with ε. Assuming only that c_i(x)?0 a.e., such that the problem can degenerate (parabolic-elliptic), we identify the homogenized problem following the values of δ=lim_ε→0ε²/λ. To cite this article: M. Mabrouk, A. Boughammoura, C. R. Mecanique 331 (2003).     

Theoretical prediction of the mass transfer coefficients in bubble columns operating in churn-turbulent flow regime. Study in Newtonian and non-Newtonian fluids under different operation conditions

 A comprehensive experimental study of the volumetric transfer coefficient k _L a with Newtonian and non-Newtonian fluids in bubble columns using CO₂ as gas phase is the objective of this work. The evaluation of the hydrodynamic characteristics of the bubble columns and delineated the different hydrodynamic regimes considering column geometry, gas flow, liquid height and type of fluid (Newtonian and non-Newtonian) suggest a general applicability of the proposed model. An explanation about of the k _L a values in non-Newtonian fluid is offered take into account shear rate, column geometry, viscosity and results reported in the literature previously. Received on 31 July 1999     

Interaction of heat transfer by conduction and radiation in a nongray planar medium

This paper considers the problem of the steady energy transfer due to combined effects of conduction and radiation in a medium with frequency dependent properties. The particular problem studied is the one-dimensional energy transfer in an absorbing, emitting, and conducting medium capable of generating heat which is bounded by two black plates. The analysis is restricted to an absorption coefficient x_v of theMilne-Eddington type, i.e., x_v(T)=αv β(T), and the frequency dependence of α(v) is approximated by a rectangular model. Temperature distribution and heat fluxes are reported in the paper for two models of spectral absorption coefficient and are compared with those for the gray case.     

The effect of suction or injection on the boundary layer flows induced by continuous surfaces stretched with prescribed skin friction

The boundary layer flow and heat transfer on a stretched surface moving with prescribed skin friction is studied for permeable surface. Three major cases are studied for isothermal surface (n=0) stretched corresponding to different dimensional skin friction boundary conditions namely; skin friction at the surface scales as (x ^?1/2) at m=0, constant skin friction at m=1/3 and skin friction scales as (x) at m=1. The constants m and n are the indices of the power law velocity and temperature exponent respectively. Similarity solutions are obtained for the boundary layer equations subject to power law temperature and velocity variation. The effect of various governing parameters, such as Prandtl number Pr, suction/injection parameter f _w , m and n are studied. The results show that for isothermal surface increasing m enhances the dimensionless heat transfer coefficient for fixed f _w at the suction case and the reverse is true at the injection case. Furthermore, for fixed m, as f _w increases the dimensionless heat transfer coefficient increases. Large enhancements are observed in the heat transfer coefficient as the temperature boundary condition along the surface changes from uniform to linear where the dimensional skin friction is of order (x) at m=1. This enhancement decreases as the suction increases.     

Thermofluiddynamic experiments with a heated and rotating circular cylinder in crossflow

Heat transfer experiments with a heated and rotating circular cylinder in an air crossflow were performed for subcritical freestream-Reynolds-number 8.3 · 10³ < Re _f < 7.1 · 10⁴ and for Nusselt-numbers up to 300. The experiments took into consideration important practical parameters such as freestream-turbulence intensity, aspect ratio, blockage ratio and temperature loading. The measured values were standardized to the ‘idealized case’ to allow adequate comparison with pertinent literature. The results encompass the range of velocity ratios α (circular/freestream-velocity) between zero and 2.5. For values α > 0.5, the typical Nusselt-number behaviour noted by some authors has been confirmed. For small velocity ratios α < 0.5, however, there are significant deviations from the heat transfer coefficients previously believed to be constant.     

A modified Lin equation for the energy balance in isotropic turbulence

At sufficiently large Reynolds numbers, turbulence is expected to exhibit scale-invariance in an intermediate (“inertial”) range of wavenumbers, as shown by power law behavior of the energy spectrum and also by a constant rate of energy transfer through wavenumber. However, there is an apparent contradiction between the definition of the energy flux (i.e., the integral of the transfer spectrum) and the observed behavior of the transfer spectrum itself. This is because the transfer spectrum T(k) is invariably found to have a zero-crossing at a single point (at k = k_*), implying that the corresponding energy flux cannot have an extended plateau but must instead have a maximum value at k = k_*. This behavior was formulated as a paradox and resolved by the introduction of filtered/partitioned transfer spectra, which exploited the symmetries of the triadic interactions (J. Phys. A: Math. Theor., 2008). In this paper we consider the more general implications of that procedure for the spectral energy balance equation, also known as the Lin equation. It is argued that the resulting modified Lin equations (and their corresponding Navier–Stokes equations) offer a new starting point for both numerical and theoretical methods, which may lead to a better understanding of the underlying energy transfer processes in turbulence. In particular the filtered partitioned transfer spectra could provide a basis for a hybrid approach to the statistical closure problem, with the different spectra being tackled using different methods.