Heat transfer through partitioned enclosures

A closed-form model for the computation of heat transfer rates through the rectangular-partitioned enclosures is investigated. The rectangular-partitioned enclosures may contain solid or gas with or without a constant and uniformly distributed heat generation. The conduction in the enclosures is considered as two-dimensional, whereas one-dimensional heat transfer through the fin-type partition is assumed. Dimensionless heat flux plots are parametrically studied by varying the aspect ratio (L/H) of the enclosure, the ratio of thermal conductivities of the enclosure to the fin-type partition (k _a /k _f ), the Biot number (β _a =h _a L/k _a ), and the reduced partition thickness (t ^*/L). It is demonstrated through an example problem that there is a large error in using one-dimensional analysis, particularly at lower values of k _a /k _f , and β _a .     

Transient heat transfer in a functionally graded convecting longitudinal fin

The study investigates the transient thermal performance of a constant area longitudinal fin made of a functionally graded material. Such a fin offers advantages that are not attainable with a traditional fin made of a homogeneous material. A numerical approach has been used to study the transient response of the fin with a step change in its base temperature. The fin is assumed to have an adiabatic tip. Three types of variations in the thermal conductivity with the longitudinal distance along the fin are considered: (a) linear, (b) quadratic, and (c) exponential. New analytical solutions for the steady state performance of the fin are derived in terms of the Bessel functions for cases (a) and (c) and in terms of the Legendre functions for case (b). These solutions provide a check on the accuracy of the transient numerical predictions for large times. The thermal performance of the fin is governed by the classical fin parameter, N _c, and the fin thermal conductivity grading parameter, a. Results are presented for the transient temperature distribution, base heat flow, convective heat loss, the energy stored in the fin and the fin efficiency for representative values of N _c and a. It is found that the transient, as well the steady state performance of the fin, is significantly affected by the functional grading of the fin material. The results presented are not only of fundamental interest but can also be used to design a functionally graded fin with the desirable steady and transient thermal characteristics.     

On cooling behavior of a vertical plate in a phase change material/water composite enclosure under pulsating heat load

This paper presents a numerical and experimental study concerning cooling characteristics of a pulsating heated vertical plate sandwiched between a substrate of phase change material (PCM) and an enclosure of water, forming a composite vertical rectangular enclosure. The vertical plate is assumed to have a uniform pulsating (on/off) volumetric heat source. The PCM considered in the present study is n-Octadecane. In the finite-difference simulation, the two-dimensional buoyancy-driven fluid flows developed in both the water-filled subenclosure and the molten PCM region of the PCM-filled subenclosure were modeled as laminar Newtonian fluid flow adhering to the Boussinesq approximation. Meanwhile, two-dimensional conduction is accounted for the plate heater as well as the solid PCM zone. Numerical results are presented to unveil the cooling behavior of the pulsating heat-generating plate through the PCM substrate and water-filled enclosure. Results of the parametric simulations reveal that the water layer has the better capability of heat dissipation than the PCM substrate. Heat dissipation from the plate through the PCM substrate is mainly via the latent heat absorption as associated with melting phenomenon. Moreover, numerical results obtained are compared with the corresponding experiments.     

Experimental and numerical investigation of dynamic heat transfer parameters in packed bed

Dynamic experiments in a nonadiabatic packed bed were carried out to evaluate the response to disturbances in wall temperature and inlet airflow rate and temperature. A two-dimensional, pseudo-homogeneous, axially dispersed plug-flow model was numerically solved and used to interpret the results. The model parameters were fitted in distinct stages: effective radial thermal conductivity (K _r) and wall heat transfer coefficient (h _w) were estimated from steady-state data and the characteristic packed bed time constant (τ) from transient data. A new correlation for the K _r in packed beds of cylindrical particles was proposed. It was experimentally proved that temperature measurements using radially inserted thermocouples and a ring-shaped sensor were not distorted by heat conduction across the thermocouple or by the thermal inertia effect of the temperature sensors.     

COUPLING BETWEEN LAMINAR FILM CONDENSATION AND NATURAL CONVECTION ON OPPOSITE SIDES OF A VERTICAL PLATE

Theoretical analyses which incorporate one-dimensional heat conduction along a plate and transverse heat conduction approximations are presented to predict the net heat transfer between laminar film condensation of a saturated vapour on one side of a vertical plate and boundary layer natural convection on the other side. It is assumed that countercurrent boundary layer flows are formed on the two sides. The governing boundary layer equations of this problem and their corresponding boundary conditions are all cast into dimensionless forms by using a non-similarity transformation. Thus the resulting system of equations can be solved by using the local non-similarity method for the boundary layer equations and a finite difference method for the heat conduction equation of the plate. The plate temperature and the heat flux through the plate are repetitively determined until the solutions for each side of the plate match. The predicted results show that the effect of Pr_c is not negligible for larger values of A* (thermal resistance ratio between natural convecti on side and condensing film side) and the approximation of transverse heat conduction overpredicts the plate temperature for lower values of R_t (thermal resistance ratio between plate and condensing film). However, no significant differences are observed between the two different approximations for higher values of R_t. © by 1997 John Wiley & Sons, Ltd.     

Heat transfer enhancement of PCM melting in 2D horizontal elliptical tube using metallic porous matrix

In this study, the melting process of ice as a phase-change material (PCM) saturated with a nickel–steel porous matrix inside a horizontal elliptical tube is investigated. Due to the low thermal conductivity of the PCM, it is motivated to augment the heat transfer performance of the system simultaneously by finding an optimum value of the aspect ratio and impregnating a metallic porous matrix into the base PCM. The lattice Boltzmann method with a double distribution function formulated based on the enthalpy method, is applied at the representative elementary volume scale under the local thermal equilibrium assumption between the PCM and porous matrix in the composite. While reducing or increasing the aspect ratio of the circular tubes leads to the expedited melting, the 90\(^{\circ }\) inclination of each elliptical tube in the case of the pure PCM melting does not affect the melting rate. With the reduction in the porosity, the effective thermal conductivity and melting rate in all tubes promoted. Although the natural convection is fully suppressed due to the significant flow blockage in the porous structure, the melting rates are generally increased in all cases.     

Natural convection heat transfer in enclosures with microemulsion phase change material slurry

This paper has dealt with the natural convection heat transfer characteristics of microemulsion slurry composed of water, fine particles of phase change material (PCM) in rectangular enclosures. The microemulsion slurry exhibited non-Newtonian pseudoplastic fluid behavior, and the phase changing process can show dramatically variations in both thermophysical and rheological properties with temperature. The experiments have been carried out separately in three subdivided regions in which the state of PCM in microemulsion is in only solid phase, two phases (coexistence of solid and liquid phases) or only liquid phase. The complicated heat transfer characteristics of natural convection have appeared in the phase changing region. The phase change phenomenon of the PCM enhanced the heat transfer in natural convection, and the Nusselt number was generalized by introducing a modified Stefan number. However, the Nusselt number did not show a linear output with the height of the enclosure, since a top conduction lid or stagnant layer was induced over a certain height of the enclosure. The Nusselt number increased with a decrease in aspect ratio (width/height of the rectangular enclosure) even including the side-wall effect. However, the microemulsion was more viscous while the PCM was in the solid phase, the side-wall effect on heat transfer was greater for the PCM in the solid region than that for the PCM in the liquid region. The correlation generalized for the PCM in a single phase is $ Nu = 1/3(1 - C_1 )Ra^{{1 \over {3.5n + 1}}} , $ where C ₁ = e ^–0.09AR for the PCM in solid phase and C ₁ = e ^–0.33AR for the PCM in liquid phase. For the PCM in the phase changing region, the correlation can be expressed as $ Nu = CRa^{{1 \over {7n + 2}}} Ste^{ - (1.9 - 1.65n)} , $ where C = 1.22 – 0.035AR for AR > 10 and C = 0.55 – 16.4e ^–1.1AR for AR < 10. The enclosure height used in the present experiments was varied from H = 5.5 [mm] to 30.4 [mm] at the fixed width W = 120 [mm] and depth D = 120 [mm]. The experiments were done in the range of modified Rayleigh number 7.0 × 10² ≤ Ra ≤ 3.0 × 10⁶, while the enclosure aspect ratio AR varied from 3.9 to 21.8.     

Exact solution of conductive heat transfer in cylindrical composite laminate

This paper presents an exact solution for steady-state conduction heat transfer in cylindrical composite laminates. This laminate is cylindrical shape and in each lamina, fibers have been wound around the cylinder. In this article heat transfer in composite laminates is being investigated, by using separation of variables method and an analytical relation for temperature distribution in these laminates has been obtained under specific boundary conditions. Also Fourier coefficients in each layer obtain by solving set of equations that related to thermal boundary layer conditions at inside and outside of the cylinder also thermal continuity and heat flux continuity between each layer is considered. In this research LU factorization method has been used to solve the set of equations.     

Conjugate Phase Change Heat Transfer in an Inclined Compound Cavity Partially Filled with a Porous Medium: A Deformed Mesh Approach

In this paper, the melting process of a PCM inside an inclined compound enclosure partially filled with a porous medium is theoretically addressed using a novel deformed mesh method. The sub-domain area of the compound enclosure is made of a porous layer and clear region. The right wall of the enclosure is adjacent to the clear region and is subject to a constant temperature of T_c. The left wall, which is connected to the porous layer, is thick and thermally conductive. The thick wall is partially subject to the hot temperature of T_h. The remaining borders of the enclosure are well insulated. The governing equations for flow and heat transfer, including the phase change effects and conjugate heat transfer at the thick wall, are introduced and transformed into a non-dimensional form. A deformed grid method is utilized to track the phase change front in the solid and liquid regions. The melting front movement is controlled by the Stefan condition. The finite element method, along with Arbitrary Eulerian–Lagrangian (ALE) moving grid technique, is employed to solve the non-dimensional governing equations. The modeling approach and the accuracy of the utilized numerical approach are verified by comparison of the results with several experimental and numerical studies, available in the literature. The effect of conjugate wall thickness, inclination angle, and the porous layer thickness on the phase change heat transfer of PCM is investigated. The outcomes show that the rates of melting and heat transfer are enhanced as the thickness of the porous layer increases. The melting rate is the highest when the inclination angle of the enclosure is 45°. An increase in the wall thickness improves the melting rate.

     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.     

Calculation of convective heat flux in the vicinity of the stagnation point for partially ionized gas flow past a body

From numerical solutions of the boundary layer equations for a four-component gas mixture (E, N⁺, N₂, and N) with gas injection, approximate formulas for the heat flux as a function of the variation of λρ/c_p and h^* across the boundary layer and the magnitude of the objection are obtained (λ is the thermal conductivity of the mixture,ρ is density, c_p is the specific heat, and h^* is the enthalpy of the ideal gas state of the mixture). An effective ambipolar diffusion coefficient D^(a)(i) is introduced, making possible finite formulas for the convective heat fluxes in the “frozen” boundary layer. We study the behavior of these coefficients within the boundary layer. A formula is obtained for convective heat flux to the wall from partially ionized air for a nine-component mixture (E, O⁺, N⁺, NO⁺, O, N, NO, O₂ N₂). Even for simpler four-component gas model three effective ambipolar diffusion coefficients are necessary: $$\begin{gathered} D^{(a)} (A) = D (A, M) D^{(a)} (I) = 2D (A, M), \hfill \\ D^{(a)} (M) = [ 1 + c_e (I)] D(A, M). \hfill \\ \end{gathered} $$ Here D(A, M) is the binary diffusion coefficient of the atoms into molecules, and c_e(I) is the ion concentration at the outer edge of the boundary layer. The assumption of an infinitely large charge-exchange cross section and the other simplifying assumptions used in [1] lead to overestimation of the magnitude of the dimensionless heat flux by 7–15% for the “frozen” boundary layer case.     

Effects of anisotropy in permeability on the two-phase flow and heat transfer in a porous cavity

This paper reports on the results of a numerical study of convection flow and heat transfer in a rectangular porous cavity filled with a phase change material under steady state conditions. The two vertical walls of the cavity are subject respectively to temperatures below and above the melting point of the PCM while adiabatic conditions are imposed on the horizontal walls. The porous medium is characterized by an anisotropic permeability tensor with the principal axes arbitrarily oriented with respect to the gravity vector. The problem is governed by the aspect ratioA, the Rayleigh numberRa, the anisotropy ratioR and the orientation angle θ of the permeability tensor. Attention is focused on these two latter parameters in order to investigate the effects of the anisotropic permeability on the fluid flow and heat transfer of the liquid/solid phase change process. The method of solution is based on the control volume approach in conjunction with the Landau-transformation to map the irregular flow domain into a rectangular one. The results are obtained for the flow field, temperature distribution, interface position and heat transfer rate forA=2.5,Ra=40, 0≤θ≤π, 0.25≤R≤4. It was found that the equilibrium state of the solid/liquid phase change process may be strongly influenced by the anisotropy ratioR as well as by the orientation angle θ of the permeability tensor. First, for a given set of parametersA,Ra andR, there exists an optimum orientation θ_max for which the flow strength, the liquid volume and the heat transfer rate are maximum. There also exists an orientation θ_min=θ_max+π/2 for which these quantities are minimum. Second, when an anisotropic medium is oriented along the optimum direction θ_max, an increase of the permeability component along that direction will increase the flow and heat transfer rate in a same order while an increase of the other permeability component only has a negligible effect. For the parameter ranges considered in the present study, it was found that the optimum direction is lying between the gravity vector and the dominant flow direction.     

Heat transfer to non-newtonian laminar falling liquid films with smooth wave free gas—liquid interface

The present paper investigates analytically the problem of heat transfer to a non-Newtonian laminar falling liquid film flowing along an inclined wall for the thermally developing and thermally developed regions. In the developing region of the temperature profile, the Nusselt number decreases monotonically until the thermal boundary layer touches the interface. But immediately after this point, the liquid film thickness decreases as well as the temperature difference in the film. The influence of parameters such as α (i.e. Fr/Re_mod ratio), γ (i.e. modified form of ?μ), modified Prandtl number and the flow behaviour index “n’ on heat transfer results is also presented.     

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.     

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).     

MHD flow and heat transfer of a UCM fluid over a stretching surface with variable thermophysical properties

In this paper we investigate the effects of temperature-dependent viscosity, thermal conductivity and internal heat generation/absorption on the MHD flow and heat transfer of a non-Newtonian UCM fluid over a stretching sheet. The governing partial differential equations are first transformed into coupled non-linear ordinary differential equation using a similarity transformation. The resulting intricate coupled non-linear boundary value problem is solved numerically by a second order finite difference scheme known as Keller-Box method for various values of the pertinent parameters. Numerical computations are performed for two different cases namely, zero and non-zero values of the fluid viscosity parameter. That is, 1/?? _r ??0 and 1/?? _r ??0 to get the effects of the magnetic field and the Maxwell parameter on the velocity and temperature fields, for several physical situations. Comparisons with previously published works are presented as special cases. Numerical results for the skin-friction co-efficient and the Nusselt number with changes in the Maxwell parameter and the fluid viscosity parameter are tabulated for different values of the pertinent parameters. The results obtained for the flow characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid phenomena, especially the UCM fluid phenomena. Maxwell fluid reduces the wall-shear stress.     

Transient thermal stress analysis of a functionally graded magneto-electro-thermoelastic strip due to nonuniform surface heating

This article is concerned with the theoretical analysis of the functionally graded magneto-electro-thermoelastic strip due to unsteady and nonuniform surface heating in the width direction. We analyze the transient thermal stress problem for a functionally graded strip constructed of the anisotropic and linear magneto-electro-thermoelastic materials using a laminated composite mode as one of theoretical approximation. The transient two-dimensional temperature is analyzed by the methods of Laplace and finite sine transformations. We obtain the solution for the simply supported and functionally graded magneto-electro-thermoelastic strip under a plane strain state. As an illustration, we carried out numerical calculations for a functionally graded strip composed of piezoelectric BaTiO₃ and magnetostrictive CoFe₂O₄, and examined the behaviors in the transient state for temperature change, stress, electric potential and magnetic potential distributions. Furthermore, the effects of the nonhomogeneity of material on the stresses, electric potential, and magnetic potential are investigated.     

Modeling thermal insulation of firefighting protective clothing embedded with phase change material

Experiments and research on heat transport through firefighting protective clothing when exposed to high temperature or intensive radiation are significant. Phase change material (PCM) takes energy when changes from solid to liquid thus reducing heat transmission. A numerical simulation of heat protection of the firefighting protective clothing embedded with PCM was studied. We focused on the temperature variation by comparing different thicknesses and position conditions of PCM combined in the clothing, as well as the melting state of PCM and human irreversible burns through a simplified one-dimensional model. The results showed it was superior to place PCM between water and proof layer and inner layer, in addition, greater thickness increased protection time while might adding extra burden to the firefighter.     

Constructal multi-scale structure of PCM-based heat sinks

This paper inquires the effectiveness of a PCM-based heat sink as a reliable solution to portable electronic devices. This sink is composed of a PCM with low thermal conductivity and fins to boost its conductivity. The optimization is subjected to fixed heat sink volume filled with PCM between vertical equidistant fins. New fins are installed in the unheated space existing in each enclosure which is not involved in thermal distribution from vertical fins to the PCM. Based on the same principle, new fins generations are augmented stepwise to the multi-scale structure. The steps of adding fins will continue up to the point that the objective function reaches its maximal value, i.e., maximizing the longest safe operation time without allowing the electronics to reach the critical temperature. The results indicate that in each length of the enclosure, the optimum volume fraction and the best fins distance values exist in which the heat sink performance becomes maximum, and adding more fins lowers the performance of the heat sink. Increasing the enclosure’s length by \(2^{n}\) does not change them. For an enclosure with constant length, the optimal number of steps for adding fins within the enclosure is a function of the fin thickness. The results indicate that increasing the thickness changes the optimal number of adding fins inside the enclosure (normally a decrease). As the fin thickness is lowered, there will be a higher effect by adding vertical fins in the enclosure. Numerical simulations cover the Rayleigh number range \(2\times 10^{5}\le \hbox {Ra}_{\mathrm{H}} \le 2.7\times 10^{8}\), where H is the heat sink height.     

Natural convection-radiation interaction on boundary layer flow along a thin vertical cylinder

A free convertion flow of an optically dense viscous incompressible fluid along a vertical thin circular cylinder has been studied with effect of radiation when the surface temperature is uniform. With appropriate transformations, the boundary layer equations governing the flow are reduced to local nonsimilarity equations. Solutions of the governing equations are obtained employing the implicit finite difference methods together with Keller box scheme as well the local nonsimilarity method with second order truncation for all ξ (nondimensional transverse curvature parameter) in the interval [0,10] and are expressed in terms of local Nusselt number for a range of values of the pertinent parameters. Effects of pertinent parameters, such as, the radiation parameter, R _d , the surface temperature parameter, θ _w , taking Prandtl number, Pr, equals 0.7 on the velocity and temperature field are also presented graphically. From the solution it is seen that increase of R _d , or θ _w leads to increase in the local rate of heat transfer coefficients. Results obtained by both the methods are obtained in excellent agreement between each other upto ξ = 10.