versão impressa ISSN 1678-5878
J. Braz. Soc. Mech. Sci. & Eng. vol.33 no.3 Rio de Janeiro jul./set. 2011
Conjugate Cooling of a Discrete Heater in Laminar Channel Flow
Thiago Antonini AlvesI; Carlos A.C. AltemaniII
IUniversidade Tecnológica Federal do Paraná, Engenharia Mecânica, 84016-210 Ponta Grossa, PR, Brasil, firstname.lastname@example.org
IIUniversidade Estadual de Campinas, Faculdade de Engenharia Mecânica, Departamento de Energia, 13083-970 Campinas, SP, Brazil, email@example.com
Electronic components are usually assembled on printed circuit boards cooled by forced airflow. When the spacing between the boards is small, there is no room to employ a heat sink on critical components. Under these conditions, the components thermal control may depend on the conductive path from the heater to the board in addition to the direct convective heat transfer to the airflow.The conjugate forced convection-conduction heat transfer from a two-dimensional strip heater flush mounted to a finite thickness wall of a parallel plates channel cooled by a laminar airflow was investigated numerically. A uniform heat flux was generated along the strip heater surface. Under steady state conditions, a fraction of the heat generation was transferred by direct convection to the airflow in the channel and the remaining fraction was transferred by conduction to the channel wall. The lower surface of the channel wall was adiabatic, so that the heat conducted from the heater to the plate eventually returned to the airflow. A portion of it returned upstream of the heater, preheating the airflow before it reached the heater surface. Due to this, it was convenient to treat the direct convection from the heater surface to the airflow by the adiabatic heat transfer coefficient. The flow was developed from the channel entrance, with constant properties.The conjugate problem was solved numerically within a single solution domain comprising both the airflow region and the solid wall of the channel. The results were obtained for the channel flow Reynolds number ranging from about 600 to 1900, corresponding to average airflow velocities from 0.5 m/s to 1.5 m/s. The effects of the solid wall to air thermal conductivities ratio were investigated in the range from 10 to 80, typical of circuit board materials. The wall thickness influence was verified from 1 mm to 5 mm. The results indicated that within these ranges, the conductive substrate wall provided a substantial enhancement of the heat transfer from the heater, accomplished by an increase of its average adiabatic surface temperature.
Keywords: conjugate heat transfer, adiabatic heat transfer coefficient, laminar channel flow, numerical analysis
The purpose of the present work was to perform an analysis of the conjugate forced convection-conduction heat transfer from a small 2D foil heater flush mounted to the lower plate (substrate) of a horizontal channel, as indicated in Fig. 1. The substrate thickness (t) and its thermal conductivity (ks) were known and a laminar developed airflow was forced into the channel. The lower surface and both ends of the substrate plate were adiabatic, so that only the upper face, in contact with the airflow, could exchange heat by forced convection. There were two thermal paths available for heat transfer from the heater to the airflow. One was by convection, directly from the heater upper surface to the airflow. The other was from the heater lower surface by conduction and spreading in the substrate wall. The heat transfer through this path was transferred back to the airflow by convection at the upper substrate surface, both upstream and downstream of the heater.This substrate conduction presents two opposing effects to the heater cooling by the airflow. First, a thermal boundary layer development upstream the heater reduces the direct convective heat transfer from the heater to the airflow. On the other hand, the effective heat transfer area to the airflow increases due to the conductive spreading upstream and downstream the heater. If this effect prevails over the first one, there will be a heat transfer enhancement, as compared to the case of an adiabatic substrate plate.
The effects of the laminar airflow rate, the fluid to substrate thermal conductivities ratio and the substrate thickness were considered in the present analysis. The total heat dissipation rate in the heater was assumed known, but its distribution into the direct convection to the airflow and the conduction to the substrate wall was obtained from analysis.The convective heat transfer was characterized by the adiabatic heat transfer coefficient had, due to its independence on the thermal boundary conditions (Moffat, 1998; Alves and Altemani, 2008). The heat transfer coefficient based on the inlet flow temperature (Tin) was also evaluated for comparison with the adiabatic coefficient.This problem depends also on the heater length and position in the channel. In the present work, the heater length was equal to the channel height and its upstream edge was centered on the substrate.
A comprehensive review of the work on conjugate forced convection-conduction heat transfer from electronic components was presented by Nakayama (1997), considering distinct heat source geometries and their arrangement on the wiring board. Here, the emphasis will be directed to works associated to flush mounted heat sources.Ramadhyaniet al.(1985) presented results of a 2D numerical analysis of the conjugate heat transfer from two discrete heat sources flush mounted on one wall of a channel. The heaters were isothermal and the wall on which they were flush mounted was a thick conductive substrate. The coolant flow was laminar and fully developed, and the analysis encompassed a range of the wall to fluid thermal conductivities ratio.The reference temperature for the Nusselt number was the inlet flow temperature.Their results indicated that the fraction of heat transfer from the heat source to the substrate can be a major contribution to the total heat transfer it increased with the substrate to fluid thermal conductivity ratio and for low fluid Peclet numbers. Incroperaet al. (1986) performed experiments and numerical analysis of the problem of flush mounted isothermal heat sources on a thick insulated wall of a horizontal channel. The flow was fully developed either in the laminar or the turbulent regime. The reference temperature for the Nusselt number was the fluid inlet temperature in the channel.The calculation domain for the numerical analysis of the conjugate problem involved the solid and fluid regions.Their numerical model predictions were in good agreement with measurements for turbulent flow, but for laminar flow the predictions were lower than the measurements.Anderson (1994) presented a technique to decouple the conjugate conductive and convective heat transfer from electronic modules on a circuit board. The adiabatic heat transfer coefficient (an invariant descriptor of the local convective heat transfer) was used with a superposition method to interface between a convection solver and a conductor solver.The method was recommended for the cooling of a module on a circuit board whenever a Biot number based on the ratio of the modules resistance to conduction over that to convection was greater than 1. Sugavanamet al. (1995) presented a numerical analysis of the conjugate heat transfer from 2D uniformly powered strip heaters flush mounted on the lower wall of a horizontal parallel plates channel. The substrate wall upper surface was cooled by laminar forced convection and the lower surface was either adiabatic or subject to laminar forced convection. The dependence of the heat transfer from the heat source was investigated for a wide range of parameters, including the substrate to coolant fluid thermal conductivities ratio (ks/k), the heat source position and the channel Reynolds number. The results indicated that preheating of the flow by board conduction increased with the ratio (ks/k), decreasing the Nusselt number on the heat source. The reference temperature for the Nusselt number was the fluid inlet temperature. It was verified that the effects of substrate thickness were important only for the conductivity ratio (ks/k) > 10. The average Nusselt on the heater surface decreased with an increase in (ks/k) for fully developed flow, due solely to substrate conduction effects. This decrease was more pronounced for smaller values of the Reynolds number, due to the relative increase of the substrate conduction effect. Considering a uniform flow at the channel entrance, developing along the channel length, the average Nusselt on the heater surface decreased as the source moved to downstream positions in the channel. Cole (1997) presented an analysis for the conjugate convective-conductive heat transfer from a 2D strip heater flush mounted on the surface of a conductive plate of finite thickness. The plate was cooled on the heater side by a fluid flow with a linear velocity profile, while the opposite surface was adiabatic. The reference temperature for the Nusselt number was that of the coolant fluid far from the plate surface, typical for external flows. The heat transfer numerical results were presented in terms of a conjugate Peclet number.It was found that for large values of this parameter, the heat transfer is dominated by the fluid flow and the strip heater size is the appropriate length scale. For small values, the appropriate scale was the solid thickness. Wang and Jaluria (2004) considered the effects of the three-dimensional conjugate mixed convection and conduction heat transfer from two heaters flush mounted on the lower wall of a horizontal rectangular duct. The two heaters were deployed with both streamwise and spanwise separation on the wall. The effects of the wall to fluid thermal conductivities ratio and the Reynolds number in the laminar regime, considering a single value of the Grashoff number, were presented in the numerical results.Among their conclusions, the heaters spanwise distribution may result in lower average temperature for the two sources. Considering a duct with an aspect ratio of 10, they reported that the numerical results of temperature distribution on the solid-fluid interface were slightly lower than those of Sugavanamet al. (1995). This was attributed to conduction in the spanwise direction, associated to the three-dimensional simulations.
In the present investigation, a numerical solution of the conjugate problem was performed for a single domain comprising both the solid and fluid regions. The reference temperature selected to describe convective heat transfer from the heater was its average adiabatic surface temperature. Thus, the corresponding average adiabatic heat transfer coefficient is independent of the amount of airflow preheating upstream of the heater. This preheating is due to the thermal wake generated by conduction upstream of the heater through the substrate plate.
Natural convection effects were considered negligible in the present work, a procedure adopted in similar investigations (e.g., Alves and Altemani, 2010;Zeng and Vafai, 2009;Davalath and Bayazitoglu, 1987 and Ramadhyani et al., 1985). In order to present a perspective for this consideration, reference will be made to the work of Kang et al. (1990).They performed experiments to investigate the mixed convection from a single heat source module mounted on a thin horizontal plate well insulated at the back surface.Their results indicated that the heat transfer from the module is dominated by forced convection when (Gr/Re5/2) < 0.9. The Reynolds number Re was based on the heater length and the forced flow average velocity. The Grashoff number wasGr = gβq"(L+2H)L3/(k υ2), where L and H indicate respectively the heater length and height, and q" is the convective heat flux from the module to the fluid. When the heater is flush mounted on the substrate plate, the height is H = 0. For any specified value of Re, their correlation gives the upper limit of Gr for which natural convection effects may be neglected.
The results of the present investigation are important for the thermal design of electronic components assembled on circuit boards cooled by forced convection.The conjugate forced convection and conduction heat transfer must be accounted for, and this analysis is most important under conditions of limited available spacing between the boards. In this case, heat sinks may not be allowed for cooling purposes and the heat transfer enhancement due to heat spreading along the substrate wall may be an important contribution to the total heat transfer (Alves, 2010).
|=||specific heat, J/(kg.K)|
|=||convective heat transfer coefficient, W/(m2.K)|
|=||channel height, m|
|=||air thermal conductivity, W/(m.K)|
|=||substrate thermal conductivity, W/(m.K)|
|=||total channel length, m|
|=||downstream length, m|
|=||heater length, m|
|=||upstream length, m|
|=||mass flow rate, kg/s|
|=||convective heat transfer rate per unit heater depth,W/m|
|=||heat flux, W/m2|
|=||Reynolds number, Eq. (2)|
|=||substrate thickness, m|
|=||velocity component along the plates, m/s|
|=||dimensionless u velocity component, Eq. (6)|
|=||Cartesian coordinates, m|
|X, Y||=||dimensionless Cartesian coordinates, Eq. (6)|
|=||dynamic viscosity, Pa.s|
|=||dimensionless temperature, Eq. (6)|
Problem formulation and heat transfer parameters
A small foil heater (length Lh) is flush mounted on a conductive substrate with thickness t and thermal conductivity ks in a horizontal channel of length L and height H, as indicated in Fig. 1. The substrate lower and side surfaces, as well as the channel upper surface, are all adiabatic.
A laminar airflow, with constant properties evaluated at 300 K, was developed from the channel entrance, with uniform velocity and the analytical parabolic velocity profile given by
The average airflow velocity and the channel hydraulic diameter defined the Reynolds number
The heater upstream edge was positioned at Lu = L/2 on the upper substrate surface. A heat dissipation rate equal to qh per unit depth normal to Fig. 1 was uniformly distributed along the heater length, with a heat flux = qh/Lh. This heat flux was distributed either directly to the airflow above the heater, with local flux (x), or to the substrate under the heater, with a local flux (x). Under steady state conditions, a local energy balance along the heater requires that
The local heat fluxes (x) and (x) were obtained by an iterative numerical solution of the energy equation in the conjugate domain indicated in Fig.1, comprising the solid substrate, the heater and the flow regions. In dimensionless form, the energy conservation equations respectively for the solid and fluid regions were expressed as follows:
The dimensionless variables were
Equations (4) and (5) were solved simultaneously, subject to the following boundary conditions. The upper and lower surfaces of the domain, as well as the upstream and downstream ends of the substrate plate were adiabatic. At the channel inlet, the fluid temperature was uniform at Tin, which corresponds to θin= 0. The outflow boundary was treated with negligible diffusion. At the substrate-fluid interface, the temperature and heat flux continuity were imposed as follows.
Equation (9) is the dimensionless form of the local energy balance expressed by Eq. (3). The heater surface temperature Th(x) and the heat fluxes (x) and (x) were not uniform along the heater length Lh. They were integrated along the heater length to obtain, respectively, the average heater temperature and the heat transfer rates qf from the heater directly to the airflow and qs to the solid substrate.
Obviously, the local energy balance expressed by Eq. (3), after integration over the heater length, results in the overall energy balance qh = qf + qs.
Two average convective heat transfer coefficients were defined on the heater surface. One was , based on the heater average adiabatic surface temperature, . This coefficient is important because it is independent of the thermal conditions upstream the heater (Alves and Altemani, 2008). The heater adiabatic surface temperature may be obtained mostly, simply considering an adiabatic substrate wall, because in this case it is equal to the airflow inlet temperature. In this case, all dissipated heat is transferred directly to the airflow (qf = qh), so that
The average adiabatic Nusselt number over the heater was expressed by
The other average heat transfer coefficient was , based on the uniform inlet flow temperature Tin. It was obtained considering a conductive substrate (ks ≠ 0), for which a fraction (qf/qh) of the total dissipation rate qh was transferred directly to the airflow, defined by:
The corresponding average Nusselt number in this case was
For a conductive substrate, the heater average temperature rise above the inlet flow temperature is due to two effects. The direct convective heat transfer rate qf from the heater to the airflow is responsible for the heater temperature rise above its adiabatic temperature, . The heat transfer rate qu conducted through the substrate wall upstream the heater returns by convection to the airflow. This effect gives rise to a thermal wake responsible for the heater adiabatic temperature rise above the inlet flow temperature, . These two effects may be added as follows:
Due to imperfect fluid mixing in the channel, the adiabatic temperature rise is always greater than the corresponding fluid mixed mean temperature rise (Tm,u Tin). The ratio of these two temperature differences defined a coefficient. An energy balance in the flow region upstream the heater related (Tm,u Tin) to qu and the airflow rate .
Similarly, the heater average temperature rise was related, by an influence coefficient , to the corresponding fluid mixed mean temperature rise (ΔTm)h, which was expressed in terms of qf and the flow rate .
Substituting Eqs. (16) and (17) into Eq. (15),
Defining the Peclet number Pe = RePr, Eq. (18) was obtained in dimensionless form.
The fractions (qf/qh) and (qu /qh) and the upstream coefficient depend on the conduction through the substrate and were obtained for each test condition. The convective coefficient , on the other hand, depends solely on the flow conditions and it was obtained from numerical tests considering an adiabatic substrate wall.
Problem formulation and heat transfer parameters
The energy equations (4) and (5) were solved simultaneously within the conjugate domain indicated in Fig. 1, comprising the solid substrate and the flow regions, using the control volumes method (Patankar, 1980). The convection and diffusion in the flow region, described by Eq. (5), with fully developed flow in the channel, were treated by the power-law scheme. The linear system of algebraic equations obtained from the discretization process was solved iteratively using the line-by-line TDMA method. The numerical results were obtained imposing the stated boundary conditions and the temperature and heat flux continuity expressed by Eqs. (7) to (9) along the substrate-fluid interface. Equations (7) and (8) were automatically satisfied by the use of the harmonic mean for evaluation of the diffusion coefficients at the control volumes interfaces. Equation (9) indicates the distribution of the uniform heat flux into the local heat fluxes (x) and (x) and required an iterative solution procedure, described as follows. Along the heater length, the link between the adjacent solid and fluid control volumes, respectively below and above the heater, was removed. An initial guess of the heat fluxes (x) and (x) was assumed and used to apply source terms to the referred solid and fluid control volumes. Then, the energy equations (4) and (5) were solved and the numerically obtained temperature distributions, together with Eq.(9), were used to evaluate a new heater temperature distribution Th(x) on the solid substrate-fluid interface. In the present work, a perfect thermal contact was assumed between the heater and the substrate. From the evaluated Th(x), new distributions for (x) and (x) were obtained and the process was repeated until convergence. Then, the upstream heating rate qu was obtained from the substrate temperature distribution and Eq. (19) was employed to determine the coefficient .
The results were obtained with a non-uniform two-dimensional numerical grid deployed on the solution domain, comprising 300 control volumes in the flow direction and 24 to 40 control volumes in the transversal direction. In the flow direction, the grid was most refined over the heater, which contained 80 uniformly distributed control volumes. Along the upstream length Lu, 170 control volumes were uniformly deployed, and 50 others were distributed along the downstream length Ld. Along the y-direction, the grid consisted of 20 non-uniform control volumes along the height H and 4 control volumes uniformly distributed along each mm of the substrate thicknesst. In the flow region, the grid size was finer near the upper and lower solid surfaces and increased with a geometric progression toward the center of the channel.
Several grids were tested before the final distribution was selected to obtain the numerical results. The initial numerical grid tests were performed under the conditions of Re = 1260 and an adiabatic substrate. The number of grid points uniformly distributed along the upstream and downstream lengths was increased until 170 points along Lu and 50 points along Ld indicated that further grid refinement would not change the results. Uniform grids were tested along the heater length Lh employing 10 to 100 grid points in the x-direction, while in the fluid region uniform grids were tested from 10 to 80 grid points in the y-direction. The Richardson extrapolation technique (De Vahl Davis, 1983) was employed and indicated that the increase of the numerical error with the grid spacing was quadratic. The exact extrapolated values of the average Nusselt number over the heater were within 0.10% of the results obtained with a uniform grid of 80 x 80 grid points over the heater. Additional tests were performed employing non-uniform grids along the y-direction in the fluid region, with a geometric progression increase of the grid spacing from the top and bottom boundarysurfaces of the domain. In this case, for the same uniform grid with 80 points along the heater and a non-uniform grid with only 20 points in the y-direction, the average Nusselt number was within 0.05% of the value previously obtained by extrapolation. When a conductive substrate was considered, additional tests verified that 4 grid points per mm along the y-direction in the solid region were enough to obtain results independent of further grid refinement. The grid used to obtain the results for a substrate thickness t/H = 0.5is presented in Fig. 2. It is non-uniform and contains 300 control volumes along the flow and 40 control volumes in the transversal direction.
The iterative solution process was interrupted when the absolute changes of the heater average temperature , the average Nusselt number , and the heat flow ratio(qs/qh) between two consecutive iterations were smaller than 10-5. The numerical results were obtained in a microcomputer (Intel® Pentium® D processor 2.8 GHz and 1GB RAM), in about 10 minutes for a typical solution considering a conductive substrate.
The numerical results were obtained considering a channel length L = 0.2 m and the plates spacing H = 0.01 m. The heater length was Lh = 0.01 m, with its upstream edge at Lu = 0.1m from the channel entrance. Five values were considered for the substrate thickness t, from 1 to 5mm.The effect of the substrate conductivity was investigated for five values of (ks/k) in the range from 10 to 80. The flow was always in the laminar regime and five values of the Reynolds number were considered, in the range from about 600 to 1900 (average air velocities from 0.5 m/s to 1.5 m/s).The air properties were obtained from tabulated values at 300 K (Incropera et al., 2006).
Shah and London (1978) presented a laminar flow correlation for the local Nusselt number along a parallel plates channel with uniform heat flux on the walls, defined as follows.
It is based on a dimensionless coordinate ξ with origin at the heater upstream length, ξ = (xLu)/(2H RePr), on the local mixed mean airflow temperatureTm(ξ) and on the channel hydraulic diameter. Their correlation is given by
This correlation can be associated with that along the heater on the adiabatic substrate in the present work. A local Nusselt number Nuad (ξ) based on the adiabatic heater temperature (which is equal to Tin when ks = 0), and on the heater length, can be related to that defined in Eq. (20) by
After substituting the correlation Eq. (21) or (22) on the left side of this equation, the resulting Nuad(ξ) from the right side was compared with the numerical distribution obtained in the present work for a particular value of Re. The results obtained for Re = 1890 are shown in Fig.3.They indicated that the predictions fromEq.(21) were 0.5% above the present numerical results,whileEq. (22) predicted values 1.4% below the numerical values.
The influence coefficient , defined by Eq. (17), was obtained from simulations for an adiabatic substrate wall, because this is the simplest procedure. For the adiabatic substrate, = Tinand all the heat is transferred directly to the airflow by convection, qf = qh. The results presented in Fig. 4 are, however, valid for either an adiabatic or a conductive substrate. They indicate that the coefficient increases with the Reynolds number, due to larger mass flow rates. This coefficient was correlated to the Reynolds number by
The heater average adiabatic Nusseltnumber , as defined by Eq. (12), was obtained considering an adiabatic substrate (qf = qh and qu = 0). The heater average adiabatic temperature was obtained from Eqs. (19) and (24) and the result was replaced on the right side of Eq. (12), giving rise to the correlation
This Nusselt number does not depend on the thermal conditions it changes only with the flow conditions and the results for air (Pr = 0.707) are presented in Fig. 5 as a function of the Reynolds number.
For a conductive substrate, an upstream airflow heating qu gives rise to a temperature increase . This temperature rise was expressed with the help of the upstream coefficient defined by Eq. (16), obtained from the simulations of the conjugate problem. Considering for example t/H= 0.5, the dependence of the coefficient on the Reynolds number is presented in Fig. 6, for distinct values of (ks/k). It increases with the Reynolds number, mainly due to larger mass flow rates. For a given Re, this coefficient decreases as (ks/k) increases, due to larger upstream heating (qu) through the substrate wall. The dependence on Re was similar for thinner substrates, but the values of increased due to larger conductive spreading resistances associated with thinner walls. Table 1 presents the numerical results of for the ratio(ks/k)=80, encompassing the five substrate thicknesses and the five Reynolds numbers considered in the present work.
The heat transfer fractions (qf/qh), (qs/qh) and(qu/qh) were also obtained from the simulations of the conjugate problem. The first two fractions obviously add to unity, as can be checked by integration of Eq. (3) along the heater length, with the definitions of Eq. (10). The fraction (qs/qh) is presented in Fig. 7 for the substrate thickness t/H= 0.5, indicating that most of the heat transfer occurs through the substrate plate, except for the lowest (ks/k). As expected, the ratio (qs/qh) increased with the ratio(ks/k), due to smaller conductive spreading resistance, and it decreased as Re increased, due to larger direct convective heat transfer. The results of (qs/qh) for (ks/k) = 80 are presented in Table 2, showing that for any Re it increases with the substrate thickness, due also to smaller conductive spreading resistance.
The fraction (qu/qh) conducted upstream of the heater for t/H = 0.5 is presented in Fig. 8, considering the effects of Re and the ratio (ks/k). The effect of the substrate thickness t on (qu/qh) is presented in Table 3 for (ks/k) = 80. They follow a trend similar to (qs/qh) with respect to the effects of t, (ks/k) and Re. Comparing the data in Table 2 and Table 3, it is seen that about 60 per cent of the conduction heat transfer from the heater to the substrate wall is released to the airflow upstream of the heater.
The heater average dimensionless temperature was obtained from the previous results and the Peclet number (Pe), using Eq.(19). The distribution of for t/H=0.5 is presented in Fig. 9, showing the expected temperature decrease as the mass flow rates increases with the Reynolds number. It also shows the heater temperature decrease as the substrate thermal conductivity increases. The results for the thermal conductivities ratio (ks/k)=80 are presented in Table 4, showing that the heater average temperature decreases with the Reynolds number and the substrate thickness. The temperatures presented in Table 4 were obtained from Eq.(19) and they matched those obtained directly from the numerical simulations within 10-4.
The heater average Nusselt number based on Tin depends on the ratio (ks/k) and the Reynolds number. It may be obtained either directly from the numerical simulations, or from the previous results for (qf/qh) and , as indicated by Eq. (14). The distribution obtained for t/H = 0.5is presented in Fig. 10. It increases with Re due to larger airflow rates, but decreases with (ks/k) due to a larger conductance of the substrate wall. As the substrate thickness t decreases, larger conductive wall resistances increase , as indicated by the numerical results presented in Table 5 for the case with (ks/k) = 80.
Compared to the adiabatic substrate, the conductive substrate provides additionally a conductive path for heat transfer from the heater. Considering the same inlet flow and heater average temperatures and flow rate in the channel, the conductive substrate causes an enhancement of heat transfer from the heater when compared to the adiabatic substrate. It was evaluated as follows. An adiabatic substrate transfers heat to the airflow only by convection, at a rate qad. For a conductive substrate, the direct convective heat transfer from the heater to the airflow is equal to qf. From the definitions of and , Eq. (11) and (13), the ratio of these two heat transfer rates was expressed by
When the heat transfer rate qf obtained from this equation is substituted into the heater overall energy balance, qh = (qf +qs), an expression for the heat transfer enhancement due to the conductive substrate is obtained in the form
This result is presented in Fig. 11 for t/H =0.5 and also in Table 6, for (ks/k) = 80. The enhancement increases with both (ks/k) and t, due to a greater conductance of the substrate wall. The effect of the Reynolds number is not obvious from Eq. (27) because, as can be seen from the previous results, the ratio increases with the Reynolds number, while the ratio (qs/qh) decreases, causing opposing trends to (qh/qad). In the present investigation, the effect of (qs/qh) was slightly dominant the heat transfer enhancement decreased slightly with the Reynolds number. The presented results indicate a significant heat transfer enhancement due to substrate conduction the heat transfer ratio (qh/qad) ranged from about 150% to 280%.
Due to the substrate adiabatic lower surface, the heat conducted upstream of the heater through the substrate wall eventually returns to the airflow and causes an increase of the heater average adiabatic temperature above the inlet flow temperature. This undesirable temperature rise was compared to the total heater average temperature rise above the inlet flow temperature, . From Eqs. (16) and (18):
The results obtained for the ratio of these temperature differences are presented in Fig. 12 for t/H = 0.5 and in Table 7 for (ks/k) = 80, within the investigated range of the Reynolds number. They show that from 25% to 45% of the total heater average temperature rise is due to the thermal wake originating from the upstream heated floor of the conductive substrate.This ratio increases with the substrate conductivity and thickness, while it decreases slightly with the Reynolds number, as expected. It should be kept in mind, however, that the average adiabatic temperature rise is the effect of the enhanced total heat transfer from the heater. Under the same inlet flow conditions and heater average temperature, the direct convective heat transfer from the heater to the airflow is smaller for the conductive substrate than for the adiabatic substrate.This decrease is, however, more than compensated by the total heat transfer rate from the heater in the case of the conductive substrate, as it has been shown by the results presented in Fig. 11 and Table 6.
The conjugate forced convection and conduction heat transfer from a strip heater flush mounted to a finite thickness wall (substrate) of a parallel plates channel were investigated numerically, using the control volumes method. The investigation was performed considering laminar airflow fully developed from the channel entrance. A uniform heat flux was released along the heater length and it was transferred to the airflow and to the conductive substrate wall. The local distribution into the direct convective heat flux (x) from the heater to the airflow and the conductive heat flux (x) from the heater to the substrate wall were not known a priori and they were obtained by an iterative procedure. These two heat fluxes were integrated along the heater length to evaluate its convective and conductive heat losses. The convective heat loss qf was expressed by means of the adiabatic heat transfer coefficient, because it is independent of the thermal conditions. The heater average temperature was expressed by means of two influence coefficients: the upstream influence coefficient , and the self-heating influence coefficient . The first coefficient evaluated the heater average adiabatic temperature rise above the inlet flow temperature in the channel, and the second coefficient was used to evaluate the average heater temperature rise above its average adiabatic temperature. Considering a dissipation rate qh in the heater, the fraction (qs/qh) conducted through the substrate and the portion (qu/qh) conducted upstream of the heater were obtained numerically as functions of the channel Reynolds number, the thermal conductivities ratio (ks/k) and the substrate thickness t.The results indicated that a substantial fraction of heat transfer occurs by conduction through the substrate. The average Nusselt numbers and were related to evaluate the heat transfer enhancement due to a conductive substrate in comparison to an adiabatic substrate, indicating values from 150% to 280%. The heater average adiabatic temperature rise due to preheating of the airflow by substrate conduction was related to its total temperature rise above Tin in dimensionless form, indicating that it represents a substantial fraction of the total, ranging from 25%to 45% in the present investigation.
The support of CNPq (Brazilian Research Council) to the first author in the form of a Doctorate Program Scholarship is gratefully acknowledged.
Alves, T.A., 2010, "Conjugate Cooling of Discrete Heaters in Channels" (In Portuguese), Ph.D. Thesis, State University of Campinas, Campinas, SP, Brazil, 129 p. [ Links ]
Alves, T.A. and Altemani, C.A.C., 2008, "Convective Cooling of Three Discrete Heat Sources in Channel Flow", J. Braz. Soc. Mech. Sci. & Eng., Vol.XXX, pp. 245-252. [ Links ]
Alves, T.A. and Altemani, C.A.C., 2010, "Thermal Design of a Protruding Heater in Laminar Channel Flow", Proceedings of the 14th International Heat Transfer Conference, Washington, USA, pp.691-700. [ Links ]
Anderson, A.M., 1994, "Decoupling Convective and Conductive Heat Transfer Using the Adiabatic Heat Transfer Coefficient", ASME J. Electronic Packaging, Vol.116, pp. 310-316. [ Links ]
Cole, K.D., 1997, "Conjugate Heat Transfer from a Small Heated Strip", I. J. Heat Mass Transfer, Vol.40, pp. 2709-2719. [ Links ]
Davalath, J. and Bayazitoglu, Y., 1987, "Forced Convection Cooling Across Rectangular Blocks", ASME J. Heat Transfer, Vol.109, pp. 321-328. [ Links ]
De Vahl Davis, G., 1983, "Natural Convection of Air in a Square Cavity: A Benchmark Numerical Solution", I. J. Numerical Methods Fluids, Vol.3, pp. 249-264. [ Links ]
Incropera, F.P., DeWitt, D.P., Bergman, T.L. and Lavine, A.S., 2006, "Fundamentals of Heat and Mass Transfer", John Wiley & Sons, Hoboken, USA, 1024 p. [ Links ]
Incropera, F.P., Kerby, J.S., Moffatt, D.F. and Ramadhyani, S., 1986, "Convection Heat Transfer from Discrete Heat Sources in a Rectangular Channel", I. J. Heat Mass Transfer, Vol.29, pp. 1051-1057. [ Links ]
Kang, B.H., Jaluria, Y. and Tewari, S.S., 1990, "Mixed Convection Transport from an Isolated Heat Source Module on a Horizontal Plate", ASME J. Heat Transfer, Vol.112, pp. 653-661. [ Links ]
Moffat, R.J., 1998, "Whats New in Convective Heat Transfer?", I. J. Heat Fluid Flow, Vol.19, pp. 90-101. [ Links ]
Nakayama, W., 1997, "Forced Convective/Conductive Conjugate Heat Transfer in Microelectronic Equipment", Annual Review Heat Transfer, Vol.8, pp. 1-45. [ Links ]
Patankar, S.V., 1980, "Numerical Heat Transfer and Fluid Flow", Hemisphere Publishing Corporation, New York, USA, 197 p. [ Links ]
Ramadhyani, S., Moffat, D.F. and Incropera, F.P., 1985, "Conjugate Heat Transfer from Small Isothermal Heat Sources Embedded in a Large Substrate", I. J. Heat Mass Transfer, Vol.28, pp. 1945-1952. [ Links ]
Shah, R.K. and London, A.L., 1978, "Laminar Flow Forced Convection in Ducts", Advances in Heat Transfer, Supplement No. 1, Academic Press, New York. [ Links ]
Sugavanam, R., Ortega, A. and Choi, C.Y., 1995, "A Numerical Investigation of Conjugate Heat Transfer from a Flush Heat Source on a Conductive Board in Laminar Channel Flow", I. J. Heat Mass Transfer, Vol.38, pp. 2969-2984. [ Links ]
Wang, Q. and Jaluria, Y., 2004, "Three-Dimensional Conjugate Heat Transfer in a Horizontal Channel with Discrete Heating", ASME J. Heat Transfer, Vol.126, pp. 642-647. [ Links ]
Zeng, Y. and Vafai, K., 2009, "An Investigation of Convective Cooling of an Array of Channel-Mounted Obstacles", Numerical Heat Transfer, Part A, Vol.55, pp. 967-982. [ Links ]
Paper received 2 February 2011.
Paper accepted 23 May 2011
Technical Editor: Horácio Vielmo