Influence Coefficient Method for Calculating Discrete Heat Source Temperature on Finite Convectively Cooled Substrates. Y.S.
|
|
- Emory Booth
- 6 years ago
- Views:
Transcription
1 Influence Coefficient Method for Calculating Discrete Heat Source Temperature on Finite Convectively Cooled Substrates Y.S. Muzychka Faculty of Engineering and Applied Science Memorial University of Newfoundland St. John s, NF, Canada, A1B 3X5 Phone: (709) Fax: (709) ABSTRACT A simple method is developed for predicting discrete heat source temperatures on a finite convectively cooled substrate. The method is based on the principle of superposition using a single source solution for the mean or maximum contact temperature of an eccentric uniform heat source on a rectangular substrate. By means of influence coefficients, the effect of neighboring source strength and location may be assessed. It is shown that the influence coefficients represent localized thermal resistances, which must be weighted according to source strength. For a system of N heat sources, there exists N effects of source strength and position on any one heat source. This includes a self effect (source of interest) and N-1 influence effects (neighboring sources). The method is developed for isotropic, orthotropic, and compound systems. Convection in the source plane is addressed for isotropic and orthotropic systems. Expressions are developed for both mean and centroidal temperature. Keywords: Conduction, Electronics Cooling, Heat Spreaders, Compound Systems, Heat Sinks, Spreading Resistance, Orthotropic Properties NOMENCLATURE a, b, c, d = linear dimensions, m A 0,A m,a n,a mn = Fourier coefficients B 0,B m,b n,b mn = Modified Fourier coefficients ˆf,f = influence coefficient, K/W h,h 1,h = contact conductance or film coefficient, W/m K k,k 1,k = thermal conductivities, W/m K m,n = indices for summations N = number of heat sources Q = heat flow rate, W R = thermal resistance, K/W t,t 1,t = total and layer thicknesses, m Assistant Professor T s = surface temperature, K T f = sink temperature, K X c,y c = heat source centroid, m Greek Symbols β mn = eigenvalues, λ m + δn δ n = eigenvalues, nπ/b γ = transform variable k z /k xy θ = temperature excess, T T f, K θ = mean temperature excess, T T f, K ˆθ = centroidal temperature excess, ˆT T f, K θ s = surface temperature excess, T s T f, K κ = relative conductivity, k /k 1 λ m = eigenvalues, mπ/a φ = spreading function ζ = dummy variable, m 1 Subscripts eff = effective value i,j = denotes the i th and j th sources xy = xy-plane z = z-plane Superscripts ( ) = mean value ( ) ˆ = centroid value INTRODUCTION Calculation of the mean or centroidal value of discrete heat sources on a rectangular subststrate is of interest in electronic packages, circuit boards, and heat sink systems. In the simplest level of analysis, the total heat dissipated from all sources may be lumped together and evenly distributed over the substrate giving rise to the lowest system temperature. Further refinements in analysis may be made by treating this lumped heat source as a single discrete heat source with an area equivalent to the total area of all heat sources. In this case, the
2 discrete lumped formulation will usually give rise to the highest system temperature, depending upon the distribution and size of individual heat sources. In most cases, these two approaches will yield useful information for preliminary sizing of cooling systems. A more refined discrete heat source analysis is often desired to minimize hot spots and evenly distribute heat flows. This paper presents a simplified method of analysis for systems with multiple discrete heat sources, which enables the effects of neighboring source location and strength to be determined. Presently, a number of methods are widely used for examining systems with multiple packages or heat sources, see Fig. 1. These include the TAMS method of Ellison [1], the GENPaK model of Culham et al. [], and full numerical solutions using Finite Element Methods (FEM), among others. The analytic methods are based on a Fourier series solutions to Laplace s equation in isotropic or multilayered systems. They differ primarily in how the local heat source is introduced into the analysis. In the TAMS method, the source is specified through the governing partial differential equation, whereas in the GENPaK method, the heat source is specified through the die plane boundary condition. The present method is based on the latter approach with a significant simplification of the remaining boundary conditions and a limitation on the number of layers. This simplified system was recently addressed by Muzychka et al. [3], who obtained a managable Fourier series based solution for a single eccentrically located heat source on an isotropic or compound substrate, which is convectively cooled with a uniform film coefficient or contact conductance. Using the principles of superposition, this solution may be used for multiple discrete sources. The present analysis considers both heat loss through the sink plane and die plane. The method of Ellison [1] and Culham et al. [] also allow for convective cooling in the die plane in addition to the sink plane. Heat Sink Heat Sources k n, t n k n-1, t n-1. k, t k 1, t 1 Fig. 1 - Multi-Source System []. In this work, the solution of Muzychka et al. [3] is re-cast in terms of influence coefficients, which allow the effects of neighbouring heat sources to be easily assessed. A temperature for each heat source may be computed in terms of these influence coefficients which shows that the total temperature excess of any given heat source is comprised of a self effect in addition to the sum of all induced effects due to neighboring sources. Using the influence coefficients it is shown that a unique thermal resistance for each heat source cannot be defined in the presence of other heat sources. The present approach does allow for more efficient computation of heat source temperatures. LITERATURE REVIEW A review of the literature reveals that several approaches for computing the thermal spreading resistance and/or heat source temperature have been developed for a rectangular substrate with single or multiple discrete sources. In the case of multiple heat sources several approaches are found in the open literature. Hein and Lenzi [4] obtained a solution for an IC package using Fourier transforms. In their development, the heat source is specified by means of a Poisson equation using a piecewise function to model discrete heat sources. Both the die plane and sink plane are convectively cooled using uniform heat transfer coefficients. Later, Kokkas [5] obtained a Fourier/Laplace transform solution for a multi-layer substrate containing discrete heat sources. The substrate base was assumed to be attached to a heat sink of fixed temperture. Discrete heat sources were dealt with using the die plane boundary condition. Ellison [1] developed a method refered to as TAMS. This method is similar to that of Hein and Lenzi [4], but considers multiple layers. More recently, Culham et al. [] developed a three dimensional Fourier series model for an electronic packaging system. There model is very general and allows for the specification of a mixed boundary condition in the die plane. Heat sources are specified through the boundary condition in the die plane. Due to the complex nature of the die plane boundary condition, numerical analysis is required to complete the solution. In all of the above methods significant effort is required to code the analysis. In the case of a single discrete heat source, several approaches are readily found in the open literature. Kadambi and Abuaf [6] obtained steady and transient solutions for a central heat source on an isotropic rectangular substrate which was convectively cooled in the sink plane. Later, Krane [7] obtained a steady solution for a similar system in which the sink plane is at a constant temperature. More recently, Yovanovich et al. [8] and Muzychka et al. [3] obtained solutions for a compound convectively cooled rectangular substrate, containing a central and eccentric heat source, respectively.
3 Finally, Muzychka et al. [9] extended these solutions to orthotropic systems, while Muzychka et al. [10] obtained a solution for a central heat source on an isotropic convectively cooled rectangular substrate with edge cooling. MATHEMATICAL MODELLING The system of interest in the present work is idealized as a rectangular substrate which may be either isotropic, orthotropic, or compound in nature. For the time being we will only consider an isotropic system, see Fig.. Later, the effects of adding a conductive layer to promote the spreading of heat, and system orthotropy will be examined. In the present system all of the edges are assumed be adiabatic, a reasonable assumption in many electronics applications where edge area is significantly less than the area of the source and sink planes. Finally, there is no heat loss through the source plane, such that all heat is dissipated through the sink plane by means of a uniform film coefficient, i.e., thermal wake effects are neglected. The addition of convection in the source plane is dealt with in a separate section. Single Source Solution The single source solution of Muzychka et al. [3] for a single eccentric uniform heat source on an isotropic substrate has the following form: θ(x,y,z) = A 0 + C 0 z+ A m cos(λ m x)[cosh(λ m z) φ m (λ m )sinh(λ m z)]+ A n cos(δ n y)[cosh(δ n z) φ n (δ n )sinh(δ n z)] + A mn cos(λ m x)cos(δ n y) [cosh(β mn z) φ mn (β mn )sinh(β mn z)] (1) where λ m = mπ/a, δ n = nπ/b, and β mn = λ m + δ n are the eigenvalues. The origin of the coordinate system is taken to be the lower left corner of the substrate. The general solution contains four components, a uniform flow solution and three spreading (or constriction) solutions which vanish when the heat flux is uniformily distributed over the entire source plane, z = 0. The genral solution is a linear superposition of each component. Application of the boundary conditions in the through plane direction yields solutions for one half of the unknown constants and gives rise to the following expression for the spreading parameter φ: φ(ζ) = ζ sinh(ζt 1) + h/k 1 cosh(ζt 1 ) ζ cosh(ζt 1 ) + h/k 1 sinh(ζt 1 ) () where ζ is replaced by λ m, δ n, or β mn, accordingly. The spreading parameter accounts for the effects of conductivity, thickness, and convection cooling. Fig. - Single Eccentric Heat Source [3]. The final Fourier coefficients A m, A n, and A mn were obtained by taking Fourier series expansions of the boundary condition in the source plane, z = 0. This yielded the following expressions for the Fourier coefficients: A m = A n = [ ( ) ( )] Q sin (Xc+c) λ m sin (Xc c) λ m abck 1 λ mφ(λ m ) [ ( ) ( )] Q sin (Yc+d) δ n sin (Yc d) δ n abdk 1 δ nφ(δ n ) A mn = 16Qcos(λ mx c )sin( 1 λ mc)cos(δ n Y c )sin( 1 δ nd) abcdk 1 β mn λ m δ n φ(β mn ) (3) where X c and Y c are the coordinates of the centroid of an arbitrarily placed heat source with respect to the lower left corner of the substrate as shown in Fig.. Finally, values for the coefficients in the uniform flow solution are given by and A 0 = Q ab ( t1 + 1 ) k 1 h C 0 = Q k 1 ab (4) (5) Centroidal Source Temperature The maximum or centroidal heat source temperature may be determined from Eq. (1) when x = X c, y = Y c,
4 and z = 0. This gives: ˆθ = A 0 + A m cos(λ m X c ) + A n cos(δ n Y c )+ A mn cos(λ m X c )cos(δ n Y c ) (6) Mean Source Temperature The mean heat source temperature is obtained by integrating the local source temperature over the source area, i.e., θ = 1 A θ(x, y, 0)dA (7) This leads to the following result for the mean temperature excess of a single eccentric heat source: may be computed using Eq. (1) evaluated at the surface θ i (x,y,0) = A i 0 + A i m cos(λ m x) + A i n cos(δ n y)+ A i mn cos(λ m x)cos(δ n y) (10) The Fourier coefficients are now evaluated at each of the i th heat source characteristics, i.e. c i,d i,x c,i,y c,i and Q i. Centroidal Source Temperature The maximum or centroidal temperature is now the sum of all heat source contributions at the point of interest. Thus, using Eq. (10) evaluated at the centroid of the j th heat source, we may write θ = A o + cos(λ m X c )sin( 1 A λ mc) m + λ m c cos(δ n Y c )sin( 1 A δ nd) n + δ n d where ˆθ j = ˆθ i (X c,j,y c,j,0) = ˆθ ij (11) A mn 4cos(δ n Y c )sin( 1 δ nd)cos(λ m X c )sin( 1 λ mc) λ m cδ n d (8) The results given by Eqs. (6) and (8) may now be used to analyze systems containing multiple heat sources. These expressions may also be used as a fundamental surface element for analyzing irregularly shaped heat sources, by discretizing the region into several rectangular strip sources. MULTIPLE DISCRETE HEAT SOURCES If more than one heat source is present (see Fig. 3), the solution for the temperature distribution on the surface of the circuit board, heat sink, or chip substrate may be obtained using superposition. Both the centroidal and mean heat source temperatures will be obtained for each heat source. Surface Temperature Distribution For N discrete heat sources the maximum temperatures occur in the source plane. The surface temperature distribution is obtained from T s (x,y,0) T f = θ s = θ i (x,y,0) (9) where θ i is the temperature excess for each heat source by itself. The temperature excess of each heat source ˆθ ij = A i 0 + A i m cos(λ m X c,j ) + A i n cos(δ n Y c,j )+ A i mn cos(λ m X c,j )cos(δ n Y c,j ) (1) The present notation θ ij, denotes the effect of the i th heat source in the region of the j th heat source. Mean Source Temperature The mean heat source temperature of an arbitrary rectangular patch of dimensions c j and d j, i.e. the j th heat source, located at X c,j and Y c,j, may be computed by integrating Eq. (7) over the region A j = c j d j, i.e., θ j = 1 A j A j θda j = 1 A j which may be written as θ j = 1 A j A j A j θ i (x,y,0)da j = θ i (x,y,0)da j (13) θ ij (14) Using Eqs. (10) results in the following expression for the mean temperature excess contribution of the i th heat
5 source in the region of the j th heat source: θ ij = A i o + A i cos(λ m X c,j )sin( 1 λ mc j ) m + λ m c j A i cos(δ n Y c,j )sin( 1 δ nd j ) n + 4 δ n d j A i mn cos(δ n Y c,j )sin( 1 δ nd j )cos(λ m X c,j )sin( 1 λ mc j ) or which may be written as: ˆθ j = Q 1 ˆf1j + Q ˆfj + + Q N ˆfNj (17) ˆθ j = Q i ˆfij (18) λ m c j δ n d j (15) Equation (14) represents the sum of the effects of all sources over an arbitrary region c j d j. Equation (15) is evaluated over the region of interest c j d j located at X c,j,y c,j. The coefficients A i 0,A i m,a i n and A i mn are then evaluated at each of the i th source parameters. Fig. 3 - Multiple Heat Sources [3]. INFLUENCE COEFFICIENT METHOD The present results may now be used to define an influence coefficient. Influence coefficients were first proposed by Negus and Yovanovich [11] for semi-infinite domains and later applied by Negus et al. [1] and Negus and Yovanovich [13,14] for multiple sources on a half space. The concept of an influence coefficient for a finite substrate was partially addressed by Hein and Lenzi [4]. Influence coefficients offer an insightful assessment of the effect of neighbouring heat sources on thermal resistance, and hence the mean or centroidal temperature excess of each discrete heat source. We begin by examining Eq. (1) and Eq. (15), for the centroidal and mean temperature excess θ. Beginning first with Eq. (11) we may write ˆθ j = ˆθ 1j + ˆθ j ˆθ Nj (16) where ˆf ij = B 0 + where Bm i cos(λ m X c,j ) + Bn i cos(δ n Y c,j )+ Bmn i cos(λ m X c,j )cos(δ n Y c,j ) B 0 = 1 ab B i m = B i n = ( t1 + 1 ) k 1 h [ ( ) ( )] sin (Xc+c) λ m sin (Xc c) λ m abck 1 λ mφ(λ m ) [ ( ) ( )] sin (Yc+d) δ n sin (Yc d) δ n abdk 1 δ nφ(δ n ) (19) B i mn = 16cos(λ mx c )sin( 1 λ mc)cos(δ n Y c )sin( 1 δ nd) abcdk 1 β mn λ m δ n φ(β mn ) (0) are modified Fourier coefficients, since Q i has now been factored out. Once again it is noted that the coefficients are evaluated at each of the i th heat source characteristics, i.e. c i,d i,x c,i, an Y c,i. Thus, the influence coefficients are only functions of the substrate properties and dimensions and of heat source geometry and location. Similarly, we may obtain an expression for the mean temperature excess of the j th heat source using Eq. (14): or which may be written as: θ j = θ 1j + θ j θ Nj (1) θ j = Q 1 f 1j + Q f j + + Q N f Nj () θ j = Q i f ij (3)
6 where f ij = B o + Bm i cos(λ m X c,j )sin( 1 λ mc j ) + λ m c j Bn i cos(δ n Y c,j )sin( 1 δ nd j ) + δ n d j B i mn 4cos(δ n Y c,j )sin( 1 δ nd j )cos(λ m X c,j )sin( 1 λ mc j ) λ m c j δ n d j (4) When i = j, the contribution is a self effect, i.e., the effect of the source acting alone. When i j the contribution to the temperature excess is an influence effect. The self effect f ii, is merely the single source thermal resistance. The influence effects f ij, are affected by two factors: source strength and the location and size of neighboring sources, i.e. a geometry effect. The influence coefficients f ij are clearly functions only of the location of the neighboring sources. Finally, we may write the temperature excess in the following matrix form: or θ 1 θ θ 3.. θ N f 11 f 1 f 1N f 1 f f N = f 31 f 3 f 3N.... f N1 f N f NN Q 1 Q Q 3.. Q N (5) {θ} = [F ij ][Q] (6) where F ij is the matrix of influence coefficients. The influence coefficient method offers a number of advantages. First, it becomes obvious what the effect a neighboring heat source has on the thermal resistance of a particular heat source. Examination of Eqs. (18) or (3) reveals that an influence effect arises by virtue of proximity and strength. In otherwords, a remote and/or weak heat source has little influence on another heat source. Second, it can be shown that the influence coefficients also possess reciprocity for the case when i j, f ji = f ij (7) This property significantly reduces computation for systems where more than five sources are present. In general, for a system of N sources, a symmetric N N matrix results for the influence coefficients. As a result of this symmetry only (N +N)/ coefficients need be computed. An upper triangular matrix is all that is needed to compute the temperature excesses. Thus the influence method offers a substantial savings in computation over the use of Eq. (9). The reciprocity is a result of the property of Greens functions [15], i.e. the potential at X c,i,y c,i due to a unit heat input at X c,j,y c,j is the same as the potential at X c,j,y c,j due to a unit heat input at X c,i,y c,i. This property also holds upon integration over a finite region. The reciprocity of the influence coefficients was also observed by Negus and Yovanovich [11] for semi-infinite regions. Thermal Resistance Finally, if we consider defining a thermal resistance R j = θ j /Q j, for each heat source, it can be shown that or ˆR j = R j = Q i Q j ˆfij (8) Q i Q j f ij (9) The above equations clearly demonstrate that the concept of thermal resistance is not strictly applicable in multiple source systems, since the total resistance of any given source depends on both proximity of the neighboring heat source, i.e. f ij, and the relative strength ratio, i.e. Q i /Q j. Changing location or strength of any source leads to a new value of thermal resistance. CONVECTION IN THE SOURCE PLANE Convection in the source plane may now be dealt with using results of Hein and Lenzi [4]. Comparison of the solution of Muzychka et al. [3] with that of Hein and Lenzi [4] shows that coefficient B o becomes: ( 1 t1 + 1 ) ab k B 0 = 1 h 1 + h 1 + h (30) 1t 1 h k 1 where h 1 denotes the film coefficient in the source plane and h denotes the film coefficient in the sink plane. Further, the spreading function φ becomes: φ(ζ) = ( h1 k 1 ζ + k 1ζ h ) ( sinh(ζt 1 ) h 1 h k 1 ζ h cosh(ζt 1 ) + sinh(ζt 1 ) ) cosh(ζt 1 ) (31) Both Eqs. (30) and (31) reduce to Eqs. () and (0), when h 1 = 0, i.e. adiabatic source plane. COMPOUND AND ORTHOTROPIC SYSTEMS The results developed earlier may be easily adapted to compound and orthotropic systems with little effort. In a recent paper, Muzychka et al. [9], applied the necessary transformations to show the relationship between isotropic and orthotropic systems. Further, using the results of Yovanovich et al. [8], one may modify the
7 isotropic model to effectively model a resistive or conductive layer placed on a rectangular substrate. Each modification is discussed below. Orthotropic Systems If the rectangular flux channel is orthotropic such that the in plane and through plane conductivities are different, i.e. k xy k z, then the following transformations may be made to apply the present method to such systems (Muzychka et al. [9]): k k eff = k xy k z (3) where, k xy and k z represent the in-plane and throughplane thermal conductivity, and t t eff = t γ (33) where γ = k z /k xy is the conductivity ratio of the orthotropic system. The orthotropic transformation is also valid for a substrate which is convectively cooled in the source plane. Compound Systems The effect of an additional layer was also examined by Muzychka et al. [3]. It was shown that the effect of an additional layer (see Fig. 4) may be handled by means of the modified spreading parameter given by: ( ) ( αe 4ζt 1 e ζt1 + ϱ e ζ(t 1+t ) αe ζ(t1+t)) φ(ζ) = (αe 4ζt1 + e ζt1 ) + ϱ ( e ζ(t1+t) + αe ζ(t1+t)) where ϱ = ζ + h/k ζ h/k and α = 1 κ 1 + κ (34) with κ = k /k 1, and ζ is replaced by λ m, δ n, or β mn, accordingly. Further, the coefficient B 0 is now given by: B 0 = 1 ( t1 + t + 1 ) (35) ab k 1 k h This modification can only be applied to the case when there is no convection in the source plane. APPLICATION OF RESULTS The results may now be applied to a simple system. Three cases will be examined: isotropic, compound, and orthotropic. The thermal property and component thicknesses are given in Table 1. In all three cases, the heat source layout summarized in Table is used along with the following substrate properties: a = 00 [mm], b = 100 [mm], and h = 100 [W/m K]. In the first case, an isotropic substrate which is cooled in teh sink plane is considered. Next, the effect of a heat spreader is examined through the addition of a conductive layer. Finally, the effect of orthotropic properties is examined. This gives rise to t = t eff = [mm] and k = k eff = [W/mK] using the properties in Table 1. Maple V Release 8 [16] was used to perform the necessary calculations. The simple code is given in the Appendix for the isotropic case. To ensure convergence, 100 terms were used in each of the single summations and 50 terms in the double summation. The results of each of the six runs are summarized in Tables 3-5, which report the centroidal and mean temperature excess for each case. In general, convergence is much slower for the centroidal temperature excess due to an alternating series. In the case of the mean temperature excess, convergence is much more rapid since all terms are positive. The results illustrate the effect that a conductive layer or layers have on the temperature. The addition of a thin conductive layer in Case B, reduces the overall temperature level in addition to flattening the temperature distribution. Similar results are also obtained for the orthotropic case where the in-plane conductivity is higher than the through plane. In both cases thermal spreading is promoted due to the presence of a higher conductivity material. The results illustrate the ease with which discrete heat source temperatures may be determined. Typical computation times ranged between 100 and 400 seconds depending on whether an isotropic or compound system was considered and whether the centroid or mean value of temperature was computed. Table 1 - Case Studies t 1 t k 1 k k xy k z [mm] [W/mK] Case A Case B Case C Fig. 4 - Compound System [3].
8 Table - Source Layout Q c d X c Y c [W] [mm] [mm] [mm] [mm] Source Source Source Source Table 3 - Isotropic Substrate Results [ C] ˆθ θ Source Source Source Source Table 4 - Compound Substrate Results [ C] ˆθ θ Source Source Source Source Table 5 - Orthotropic Substrate Results [ C] ˆθ θ Source Source Source Source SUMMARY AND CONCLUSIONS A simple method for predicting mean and centroidal heat source temperature was developed by means of an influence coefficient. It was shown that this coefficient is only a function of source location and size. It was also shown that discrete heat source thermal resistance is weighted according to the relative source strength ratios. Further, it was also shown that the influence coefficients lend themselves to more efficient computation due to the reciprocity property. Several examples were computed to demonstrate the ease of application. Finally, the method was developed for isotropic, compound, and orthotropic systems. Modification of the basic equations for the case where heat is dissipated in the source plane was also discussed. ACKNOWLEDGMENTS The author acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC). The author also thanks Prof. M.M. Yovanovich for comments given during manuscript preparation. REFERENCES [1] Ellison, G., Thermal Computations for Electronic Equipment, Krieger Publishing, Malabar, FL, [] Culham, J.R., Yovanovich, M.M., and Lemczyk, T.F., Thermal Characterization of Electronic Packages Using a Three-Dimensional Fourier Series Solution, Journal of Electronic Packaging, 000, Vol. 1, pp [3] Muzychka, Y.S., Yovanovich, M.M., and Culham, J.R., Thermal Spreading Resistance of Eccentric Heat Sources on Rectangular Flux Channels, Journal of Electronic Packaging, Vol. 15, pp , 003. [4] Hein, V.L. and Lenzi, V.D., Thermal Analysis of Substrates and Integrated Circuits, pp , Electronics Components Conference, [5] Kokkas, A., Thermal Analysis of Multiple-Layer Structures, IEEE Transactions on Electron Devices, Vol. Ed-1, No. 14, pp , [6] Kadambi, V. and Abuaf, N., Analysis of Thermal Response for Power Chip Packages, IEEE Trans. Elec. Dev., Vol. ED-3, No. 6, [7] Krane, M.J.H., Constriction Resistance in Rectangular Bodies, Journal of Electronic Packaging, Vol. 113, 1991, pp [8] Yovanovich, M.M., Muzychka, Y.S., and Culham, J.R., Spreading Resistance of Isoflux Rectangles and Strips on Compound Flux Channels, Journal of Thermophysics and Heat Transfer, Vol. 13, 1999, pp [9] Muzychka, Y.S., Yovanovich, M.M., and Culham, J.R., Thermal Spreading Resistances in Compound and Orthotropic Systems, Journal of Thermophysics and Heat Transfer, In Press, 003. [10] Muzychka, Y.S., Culham, J.R., and Yovanovich, M.M., Thermal Spreading Resistances of Rectangular Flux Channels: Part II Edge Cooling, 36th AIAA Thermophysics Conference, Orlando, FL, 003. [11] Negus, K.J. and Yovanovich, M.M., Thermal Resistance of Arbitrarily Shaped Contacts, Numerical Methods in Thermal Problems, Proceedings of the 3rd International Conference, Seattle, WA, 1983, pp [1] Negus, K.J., Yovanovich, M.M., and DeVaal, J.W., Development of Thermal Constriction Resistance for
9 Anisotropic Rough Surfaces by the Method of Infinite Images, National Heat Trnasfer Conference, Denver, CO, [13] Negus, K.J. and Yovanovich, M.M., Transient Temperature Rise at Surface Due to Arbitrary Contacts on Half Spaces, Transactions of the CSME, 1987, Vol. 13, pp [14] Negus, K.J. and Yovanovich, M.M., Thermal Computations in a Semiconductor Die Using Surface Elements and Infinite Images, International Symposium on Cooling Technology in Electronic Equipment, Honolulu, HI, 1987, pp [15] Morse, P.M., and Feshbach, H., Methods of Theoretical Physics, Part I, McGraw-Hill, New York, [16] Maple TM Release 8, Waterloo Maple Inc., Waterloo, ON, 00. APPENDIX Simple Maple Release 8 code for Case A results. Define Influence Coefficient > restart; > lambda:=m*pi/a; > delta:=n*pi/b; > beta:=sqrt(lambda^+delta^); > phi:=zeta->(zeta*sinh(zeta*t)+h/k* cosh(zeta*t))/(zeta*cosh(zeta*t)+ h/k*sinh(zeta*t)); > B[0]:=1/(a*b)*(t/k+1/h); > Bm[i]:=*(sin((*X[i]+c[i])*lambda/)- sin((*x[i]-c[i])*lambda/)) /(a*b*c[i]*k*lambda^*phi(lambda)); > Bn[i]:=*(sin((*Y[i]+d[i])*delta/)- sin((*y[i]-d[i])*delta/)) /(a*b*d[i]*k*delta^*phi(delta)); > Bmn[i]:=16*(cos(lambda*X[i])*sin(1/* lambda*c[i])*cos(delta*y[i])*sin(1/*delta*d[i]) )/(a*b*c[i]*d[i]*k*lambda*delta*beta*phi(beta)); > f[i]:=value(b[0]+*add(bm[i]*cos(lambda*x[j]) *sin(1/*lambda*c[j])/(lambda*c[j]),..100) +*add(bn[i]*cos(delta*y[j])*sin(1/*delta*d[j]) /(delta*d[j]),..100)+4*add(add(bmn[i]*cos( lambda*x[j])*sin(1/*lambda*c[j])*cos(delta*y[j] )*sin(1/*delta*d[j])/(lambda*c[j]*delta*d[j]),..100),..100)): Input System Parameters > baseparameters:={ a=0.,b=0.1,k=10,h=100,t=0.01}; > sourceparameters:={ c[1]=0.0,d[1]=0.0,x[1]=0.04,y[1]=0.03,q[1]=10, c[]=0.03,d[]=0.04,x[]=0.095,y[]=0.03,q[]=15, c[3]=0.03,d[3]=0.07,x[3]=0.155,y[3]=0.045,q[3]=5, c[4]=0.05,d[4]=0.01,x[4]=0.055,y[4]=0.075,q[4]=10}; Calculate Influence Coefficients > f1s:=[seq(evalf(subs(j=1,i=n,baseparameters, sourceparameters,f[i])),..4)]; > fs:=[seq(evalf(subs(j=,i=n,baseparameters, sourceparameters,f[i])),..4)]; > f3s:=[seq(evalf(subs(j=3,i=n,baseparameters, sourceparameters,f[i])),..4)]; > f4s:=[seq(evalf(subs(j=4,i=n,baseparameters, sourceparameters,f[i])),..4)]; Calculate Source Temperature Excesses > Source1Theta:=subs(j=1,sourceparameters, add(q[i]*f1s[i],..4)); > SourceTheta:=subs(j=,sourceparameters, add(q[i]*fs[i],..4)); > Source3Theta:=subs(j=3,sourceparameters, add(q[i]*f3s[i],..4)); > Source4Theta:=subs(j=4,sourceparameters, add(q[i]*f4s[i],..4));
CALCULATION of the mean or centroidal value of discrete
636 IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL 29, NO 3, SEPTEMBER 2006 Influence Coefficient Method for Calculating Discrete Heat Source Temperature on Finite Convectively Cooled
More informationInfluence of Geometry and Edge Cooling on Thermal Spreading Resistance
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 20, No. 2, April June 2006 Influence of Geometry Edge Cooling on Thermal Spreading Resistance Y. S. Muzychka Memorial University of Newfoundl, St. John s,
More informationThermal Spreading Resistance in Compound and Orthotropic Systems
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 8, No., January March 004 Thermal Spreading Resistance in Compound and Orthotropic Systems Y. S. Muzychka Memorial University of Newfoundland, St. John s,
More informationOptimization of Plate Fin Heat Sinks Using Entropy Generation Minimization
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL 24, NO 2, JUNE 2001 159 Optimization of Plate Fin Heat Sinks Using Entropy Generation Minimization J Richard Culham, Member, IEEE, and Yuri
More informationTHE generalsolutionfor the spreadingresistanceof a ux speci-
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 5, No. 3, July September 00 Thermal Spreading Resistances in Compound Annular Sectors Y. S. Muzychka Memorial University of Newfoundland, St. John s, Newfoundland
More informationTHE recently published paper of Hui and Tan [1] was
168 IEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY PART A, VOL. 21, NO. 1, MARCH 1998 Analytical Modeling of Spreading Resistance in Flux Tubes, Half Spaces, Compound Disks M.
More informationExperimental Investigation of Heat Transfer in Impingement Air Cooled Plate Fin Heat Sinks
Zhipeng Duan Graduate Research Assistant e-mail: zpduan@engr.mun.ca Y. S. Muzychka Associate Professor Member ASME e-mail: yuri@engr.mun.ca Faculty of Engineering and Applied Science, Memorial University
More informationMicroelectronics Heat Transfer Laboratory
Microelectronics Heat Transfer Laboratory Department of Mechanical Engineering University of Waterloo Waterloo, Ontario, Canada http://www.mhtl.uwaterloo.ca Outline Personnel Capabilities Facilities Research
More informationModeling Transient Conduction in Enclosed Regions Between Isothermal Boundaries of Arbitrary Shape
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 19, No., July September 2005 Modeling Transient Conduction in Enclosed Regions Between Isothermal Boundaries of Arbitrary Shape Peter Teertstra, M. Michael
More informationModeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL. 24, NO. 2, JUNE 2001 207 Modeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate Mirko Stevanović, M. Michael
More informationA Simple Closed Form Solution to Single Layer Heat Spreading Angle Appropriate for Microwave Hybrid Modules
Journal of Electronics Cooling and Thermal Control, 2016, 6, 52-61 Published Online June 2016 in SciRes. http://www.scirp.org/journal/jectc http://dx.doi.org/10.4236/jectc.2016.62005 A Simple Closed Form
More informationTwo-Dimensional Numerical Investigation on Applicability of 45 Heat Spreading Angle
Journal of Electronics Cooling and Thermal Control, 24, 4, - Published Online March 24 in SciRes. http://www.scirp.org/journal/jectc http://dx.doi.org/.4236/jectc.24.4 Two-Dimensional Numerical Investigation
More informationTHERMAL spreading resistance occurs whenever heat leaves a
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 3, No. 4, October December 999 Spreading Resistance of Iso ux Rectangles and Strips on Compound Flux Channels M. M. Yovanovich, Y. S. Muzychka, and J. R.
More informationUsing FLOTHERM and the Command Center to Exploit the Principle of Superposition
Using FLOTHERM and the Command Center to Exploit the Principle of Superposition Paul Gauché Flomerics Inc. 257 Turnpike Road, Suite 100 Southborough, MA 01772 Phone: (508) 357-2012 Fax: (508) 357-2013
More informationModeling of Natural Convection in Electronic Enclosures
Peter M. Teertstra e-mail: pmt@mhtlab.uwaterloo.ca M. Michael Yovanovich J. Richard Culham Microelectronics Heat Transfer Laboratory, Department of Mechanical Engineering, University of Waterloo, Waterloo,
More informationJournal of Solid and Fluid Mechanics. An approximate model for slug flow heat transfer in channels of arbitrary cross section
Vol. 2, No. 3, 2012, 1 7 Journal of Solid and Fluid Mechanics Shahrood University of Technology An approximate model for slug flow heat transfer in channels of arbitrary cross section M. Kalteh 1,*, A.
More informationReview of Thermal Joint Resistance Models for Non-Conforming Rough Surfaces in a Vacuum
Review of Thermal Joint Resistance Models for Non-Conforming Rough Surfaces in a Vacuum M. Bahrami J. R. Culham M. M. Yovanovich G. E. Schneider Department of Mechanical Engineering Microelectronics Heat
More informationConstructal multi-scale design of compact micro-tube heat sinks and heat exchangers
JID:THESCI AID:2493 /FLA [m5+; v 1.60; Prn:29/06/2006; 9:31] P.1 (1-8) International Journal of Thermal Sciences ( ) www.elsevier.com/locate/ijts Constructal multi-scale design of compact micro-tube heat
More informationTHERMAL DESIGN OF POWER SEMICONDUCTOR MODULES FOR MOBILE COMMNICATION SYSYTEMS. Yasuo Osone*
Nice, Côte d Azur, France, 27-29 September 26 THERMAL DESIGN OF POWER SEMICONDUCTOR MODULES FOR MOBILE COMMNICATION SYSYTEMS Yasuo Osone* *Mechanical Engineering Research Laboratory, Hitachi, Ltd., 832-2
More informationThermal Characterization of Electronic Packages Using a Three-Dimensional Fourier Series Solution
J. R. Culham Associate Professor and Director, Mem. ASME, M. M. Yovanovich Professor Emeritus and Principal Scientific Advisor, Fellow ASME, Microelectronics Heat Transfer Laboratory, Department of Mechanical
More informationConjugate heat transfer from an electronic module package cooled by air in a rectangular duct
Conjugate heat transfer from an electronic module package cooled by air in a rectangular duct Hideo Yoshino a, Motoo Fujii b, Xing Zhang b, Takuji Takeuchi a, and Souichi Toyomasu a a) Fujitsu Kyushu System
More informationTHERMAL PERFORMANCE EVALUATION AND METHODOLOGY FOR PYRAMID STACK DIE PACKAGES
THERMAL PERFORMANCE EVALUATION AND METHODOLOGY FOR PYRAMID STACK DIE PACKAGES Krishnamoorthi.S, *W.H. Zhu, C.K.Wang, Siew Hoon Ore, H.B. Tan and Anthony Y.S. Sun. Package Analysis and Design Center United
More informationconstriction/spreading RESISTANCE MODEL FOR ELECTRONICS PACKAGING
ASME/JSME Thermal Engineering Conference: Volume 4 ASME 1995 constrcton/spreadng RESSTANCE MODEL FOR ELECTRONCS PACKAGNG Seri Lee Aavid Engineering, nc. Laconia, New Hampshire Seaho Song, Van Au Bell-Northern
More informationMOVING HEAT SOURCES IN A HALF SPACE: EFFECT OF SOURCE GEOMETRY
Proceedings of the ASME 2009 Heat Transfer Summer Conference HT2009 July 19-23, 2009, San Francisco, California, USA HT2009-88562 MOVING HEAT SOURCES IN A HALF SPACE: EFFECT OF SOURCE GEOMETRY Mohsen Akbari*,
More informationStacked Chip Thermal Model Validation using Thermal Test Chips
Stacked Chip Thermal Model Validation using Thermal Test Chips Thomas Tarter Package Science Services ttarter@pkgscience.com Bernie Siegal Thermal Engineering Associates, Inc. bsiegal@thermengr.net INTRODUCTION
More informationELEC9712 High Voltage Systems. 1.2 Heat transfer from electrical equipment
ELEC9712 High Voltage Systems 1.2 Heat transfer from electrical equipment The basic equation governing heat transfer in an item of electrical equipment is the following incremental balance equation, with
More informationDeveloping an Empirical Correlation for the Thermal Spreading Resistance of a Heat Sink
San Jose State University SJSU ScholarWorks Master's Theses Master's Theses and Graduate Research Spring 2016 Developing an Empirical Correlation for the Thermal Spreading Resistance of a Heat Sink Andrew
More informationExperimental Analysis of Wire Sandwiched Micro Heat Pipes
Experimental Analysis of Wire Sandwiched Micro Heat Pipes Rag, R. L. Department of Mechanical Engineering, John Cox Memorial CSI Institute of Technology, Thiruvananthapuram 695 011, India Abstract Micro
More informationOptimization of Heat Spreader. A thesis presented to. the faculty of. In partial fulfillment. of the requirements for the degree.
Optimization of Heat Spreader A thesis presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillment of the requirements for the degree Master of
More informationTools for Thermal Analysis: Thermal Test Chips Thomas Tarter Package Science Services LLC
Tools for Thermal Analysis: Thermal Test Chips Thomas Tarter Package Science Services LLC ttarter@pkgscience.com INTRODUCTION Irrespective of if a device gets smaller, larger, hotter or cooler, some method
More informationPressure Drop of Impingement Air Cooled Plate Fin Heat Sinks
Zhipeng Duan Graduate Research Assistant e-mail: zpduan@engr.mun.ca Y. S. Muzychka Associate Professor Mem. ASME e-mail: yuri@engr.mun.ca Faculty of Engineering and Applied Science, Memorial University
More informationIEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY PART A, VOL. 20, NO. 4, DECEMBER
IEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY PART A, VOL. 20, NO. 4, DECEMBER 1997 463 Pressure Loss Modeling for Surface Mounted Cuboid-Shaped Packages in Channel Flow Pete
More informationThermal Resistance (measurements & simulations) In Electronic Devices
Thermal Resistance (measurements & simulations) In Electronic Devices A short online course PART 3 Eric Pop Electrical Engineering, Stanford University 1 Topics 1) Basics of Joule Heating 2) Heating in
More informationBoundary Condition Dependency
Boundary Condition Dependency of Junction to Case Thermal Resistance Introduction The junction to case ( ) thermal resistance of a semiconductor package is a useful and frequently utilized metric in thermal
More informationThermal Resistance Measurement
Optotherm, Inc. 2591 Wexford-Bayne Rd Suite 304 Sewickley, PA 15143 USA phone +1 (724) 940-7600 fax +1 (724) 940-7611 www.optotherm.com Optotherm Sentris/Micro Application Note Thermal Resistance Measurement
More informationChapter 3: Transient Heat Conduction
3-1 Lumped System Analysis 3- Nondimensional Heat Conduction Equation 3-3 Transient Heat Conduction in Semi-Infinite Solids 3-4 Periodic Heating Y.C. Shih Spring 009 3-1 Lumped System Analysis (1) In heat
More informationSPREADING resistance, also sometimes referred to as constriction
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL. 33, NO., JUNE 00 67 Thermal Spreading Resistance of Arbitrary-Shape Heat Sources on a Half-Space: A Unified Approach Ehsan Sadeghi, Majid
More information7-9 October 2009, Leuven, Belgium Electro-Thermal Simulation of Multi-channel Power Devices on PCB with SPICE
Electro-Thermal Simulation of Multi-channel Power Devices on PCB with SPICE Torsten Hauck*, Wim Teulings*, Evgenii Rudnyi ** * Freescale Semiconductor Inc. ** CADFEM GmbH Abstract In this paper we will
More informationResistance Post-Trim Drift Index for Film Resistors to be Trimmed Abstract Introduction
Resistance Post-Trim Drift Index for Film Resistors to be Trimmed (K. Schimmanz, Numerical and Applied Mathematics, Technical University Cottbus, P.O. Box 10 13 44, Cottbus 03013 Germany schimm@math.tu-cottbus.de
More informationThe Increasing Importance of the Thermal Management for Modern Electronic Packages B. Psota 1, I. Szendiuch 1
Ročník 2012 Číslo VI The Increasing Importance of the Thermal Management for Modern Electronic Packages B. Psota 1, I. Szendiuch 1 1 Department of Microelectronics, Faculty of Electrical Engineering and
More informationNUMERICAL SIMULATION OF CONJUGATE HEAT TRANSFER FROM MULTIPLE ELECTRONIC MODULE PACKAGES COOLED BY AIR
Proceedings of IPACK03 International Electronic Packaging Technical Conference and Exhibition July 6 11 2003 Maui Hawaii USA InterPack2003-35144 NUMERICAL SIMULATION OF CONJUGATE HEAT TRANSFER FROM MULTIPLE
More informationSEMICONDUCTOR THERMAL MEASUREMENT PROCEDURE
SEMICONDUCTOR TERMAL MEASUREMENT PROCEDURE The following general procedure is equally applicable to either JEDEC or SEMI thermal measurement standards for integrated circuits and thermal test die. 1. Determine
More informationAnalytical solutions of heat transfer for laminar flow in rectangular channels
archives of thermodynamics Vol. 35(2014), No. 4, 29 42 DOI: 10.2478/aoter-2014-0031 Analytical solutions of heat transfer for laminar flow in rectangular channels WITOLD RYBIŃSKI 1 JAROSŁAW MIKIELEWICZ
More informationLocalized TIM Characterization Using Deconstructive Analysis
Localized TIM Characterization Using Deconstructive Analysis By Phillip Fosnot and Jesse Galloway Amkor Technology, Inc. 1900 South Price Road Chandler, AZ 85286 Phillip.Fosnot@amkor.com Abstract Characterizing
More informationAN ANALYTICAL THERMAL MODEL FOR THREE-DIMENSIONAL INTEGRATED CIRCUITS WITH INTEGRATED MICRO-CHANNEL COOLING
THERMAL SCIENCE, Year 2017, Vol. 21, No. 4, pp. 1601-1606 1601 AN ANALYTICAL THERMAL MODEL FOR THREE-DIMENSIONAL INTEGRATED CIRCUITS WITH INTEGRATED MICRO-CHANNEL COOLING by Kang-Jia WANG a,b, Hong-Chang
More information3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationLecture 5 - Electromagnetic Waves IV 19
Lecture 5 - Electromagnetic Waves IV 9 5. Electromagnetic Waves IV 5.. Symmetry in EM In applications, we often have symmetry in the structures we are interested in. For example, the slab waveguide we
More informationUsing Excel to Implement the Finite Difference Method for 2-D Heat Transfer in a Mechanical Engineering Technology Course
Paper ID #9196 Using Excel to Implement the Finite Difference Method for -D Heat ransfer in a Mechanical Engineering echnology Course Mr. Robert Edwards, Pennsylvania State University, Erie Bob Edwards
More informationHEAT TRANSFER THERMAL MANAGEMENT OF ELECTRONICS YOUNES SHABANY. C\ CRC Press W / Taylor Si Francis Group Boca Raton London New York
HEAT TRANSFER THERMAL MANAGEMENT OF ELECTRONICS YOUNES SHABANY C\ CRC Press W / Taylor Si Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group, an informa business
More informationSimulation of the Temperature Profile During Welding with COMSOL Multiphysics Software Using Rosenthal s Approach
Simulation of the Temperature Profile During Welding with COMSOL Multiphysics Software Using Rosenthal s Approach A. Lecoanet*, D. G. Ivey, H. Henein Department of Chemical & Materials Engineering, University
More informationFast field solver programs for thermal and electrostatic analysis of microsystem elements
Fast field solver programs for thermal and electrostatic analysis of microsystem elements Vladimir Székely Márta Rencz szekely@eet.bme.hu rencz@eet.bme.hu Technical University of Budapest, H-1521 Budapest,
More informationConduction Heat Transfer. Fourier Law of Heat Conduction. x=l Q x+ Dx. insulated x+ Dx. x x. x=0 Q x A
Conduction Heat Transfer Reading Problems 10-1 10-6 10-20, 10-48, 10-59, 10-70, 10-75, 10-92 10-117, 10-123, 10-151, 10-156, 10-162 11-1 11-2 11-14, 11-20, 11-36, 11-41, 11-46, 11-53, 11-104 Fourier Law
More informationChapter 3 Three Dimensional Finite Difference Modeling
Chapter 3 Three Dimensional Finite Difference Modeling As has been shown in previous chapters, the thermal impedance of microbolometers is an important property affecting device performance. In chapter
More informationAustralian Journal of Basic and Applied Sciences. Numerical Investigation of Flow Boiling in Double-Layer Microchannel Heat Sink
AENSI Journals Australian Journal of Basic and Applied Sciences ISSN:1991-8178 Journal home page: www.ajbasweb.com Numerical Investigation of Flow Boiling in Double-Layer Microchannel Heat Sink Shugata
More informationTRANSIENT HEAT CONDUCTION
TRANSIENT HEAT CONDUCTION Many heat conduction problems encountered in engineering applications involve time as in independent variable. This is transient or Unsteady State Heat Conduction. The goal of
More informationABSTRACT. 1. Introduction
Journal of Electronics Cooling and Thermal Control, 03, 3, 35-4 http://dx.doi.org/0.436/jectc.03.3005 Published Online March 03 (http://www.scirp.org/journal/jectc) Evaluation of Inherent Uncertainties
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationPower Stage Thermal Design for DDX Amplifiers
Power Stage Thermal Design for DDX Amplifiers For Applications Assistance Contact: Apogee Technical Support e-mail: support@apogeeddx.com CONTROLLED DOCUMENT: P_901-000002_Rev04 Power Stage Thermal Design
More informationEffects of Chrome Pattern Characteristics on Image Placement due to the Thermomechanical Distortion of Optical Reticles During Exposure
Effects of Chrome Pattern Characteristics on Image Placement due to the Thermomechanical Distortion of Optical Reticles During Exposure A. Abdo, ab L. Capodieci, a I. Lalovic, a and R. Engelstad b a Advanced
More informationAnalysisofElectroThermalCharacteristicsofaConductiveLayerwithCracksandHoles
Global Journal of Researches in Engineering Mechanical and Mechanics Engineering Volume 14 Issue 1 Version 1.0 Type: Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
More informationTable of Contents. Foreword... Introduction...
Table of Contents Foreword.... Introduction.... xi xiii Chapter 1. Fundamentals of Heat Transfer... 1 1.1. Introduction... 1 1.2. A review of the principal modes of heat transfer... 1 1.2.1. Diffusion...
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationPROBLEM Node 5: ( ) ( ) ( ) ( )
PROBLEM 4.78 KNOWN: Nodal network and boundary conditions for a water-cooled cold plate. FIND: (a) Steady-state temperature distribution for prescribed conditions, (b) Means by which operation may be extended
More informationIntroduction: Plate Fin Heat Sinks Heat transfer enhancement for air cooled applications: { increase eective surface area { decrease thermal resistanc
Heat Transfer Laoratory Microelectronics ofwaterloo University Analytical Forced Convection Modeling of Plate Fin Heat Sinks P. Teertstra, M.M. Yovanovich and J.R. Culham Department of Mechanical Engineering
More informationC ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127
C ONTENTS Preface xviii Nomenclature xxvi CHAPTER ONE BASICS OF HEAT TRANSFER 1 1-1 Thermodynamics and Heat Transfer 2 Application Areas of Heat Transfer 3 Historical Background 3 1-2 Engineering Heat
More informationDesign and Optimization of Horizontally-located Plate Fin Heat Sink for High Power LED Street Lamps
Design and Optimization of Horizontally-located Plate Fin Heat Sink for High Power LED Street Lamps Xiaobing Luo 1,2*, Wei Xiong 1, Ting Cheng 1 and Sheng Liu 2, 3 1 School of Energy and Power Engineering,
More information6.2 Modeling of Systems and Components
Chapter 6 Modelling of Equipment, Processes, and Systems 61 Introduction Modeling is one of the most important elements of thermal system design Most systems are analyzed by considering equations which
More informationOptimal Design Methodology of Plate-Fin Heat Sinks for Electronic Cooling Using Entropy Generation Strategy
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL. 27, NO. 3, SEPTEMBER 2004 551 Optimal Design Methodology of Plate-Fin Heat Sinks for Electronic Cooling Using Entropy Generation Strategy
More informationHEAT TRANSFER FROM FINNED SURFACES
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 HEAT TRANSFER FROM FINNED SURFACES Mehmet Kanoglu Copyright The McGraw-Hill Companies,
More informationReview: Conduction. Breaking News
CH EN 3453 Heat Transfer Review: Conduction Breaking News No more homework (yay!) Final project reports due today by 8:00 PM Email PDF version to report@chen3453.com Review grading rubric on Project page
More informationSIMULATION AND ASSESSMENT OF AIR IMPINGEMENT COOLING ON SQUARED PIN-FIN HEAT SINKS APPLIED IN PERSONAL COMPUTERS
20 Journal of Marine Science and Technology, Vol. 13, No. 1, pp. 20-27 (2005) SIMULATION AND ASSESSMENT OF AIR IMPINGEMENT COOLING ON SQUARED PIN-FIN HEAT SINKS APPLIED IN PERSONAL COMPUTERS Hwa-Chong
More informationEXPERIMENTAL INVESTIGATION OF HIGH TEMPERATURE THERMAL CONTACT RESISTANCE WITH INTERFACE MATERIAL
HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 26 July 214 Orlando, Florida EXPERIMENTAL INVESTIGATION OF HIGH TEMPERATURE THERMAL CONTACT RESISTANCE WITH
More informationThe Teaching of Unsteady Heat Conduction Using the Thermo-electric Analogy and the Code pspice. Nonlineal Models
International Conference on Engineering Education and Research "Progress Through Partnership" VSB-TUO, Ostrava, ISSN 56-58 The Teaching of Unsteady Heat Conduction Using the Thermo-electric Analogy and
More informationChapter 9 NATURAL CONVECTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition in SI Units Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 9 NATURAL CONVECTION PM Dr Mazlan Abdul Wahid Universiti Teknologi
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationOptimizing Diamond Heat Spreaders for Thermal Management of Hotspots for GaN Devices
Optimizing Diamond Heat Spreaders for Thermal Management of Hotspots for GaN Devices Thomas Obeloer*, Bruce Bolliger Element Six Technologies 3901 Burton Drive Santa Clara, CA 95054 *thomas.obeloer@e6.com
More informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
More informationCOMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS
HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 214 Orlando, Florida COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES
More informationOptimization of Microchannel
Optimization of Microchannel Heat Exchangers It s no secret that power densities in electronics have continued to rise, and researchers have been forced to explore new thermal management technologies to
More informationChapter 4: Transient Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University
Chapter 4: Transient Heat Conduction Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Assess when the spatial
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationSupplementary Information for On-chip cooling by superlattice based thin-film thermoelectrics
Supplementary Information for On-chip cooling by superlattice based thin-film thermoelectrics Table S1 Comparison of cooling performance of various thermoelectric (TE) materials and device architectures
More informationan alternative approach to junction-to-case thermal resistance measurements
an alternative approach to junction-to-case thermal resistance measurements Bernie Siegal Thermal Engineering Associates, Inc. Introduction As more and more integrated circuits dissipate power at levels
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationUsing LBM to Investigate the Effects of Solid-Porous Block in Channel
International Journal of Modern Physics and Applications Vol. 1, No. 3, 015, pp. 45-51 http://www.aiscience.org/journal/ijmpa Using LBM to Investigate the Effects of Solid-Porous Bloc in Channel Neda Janzadeh,
More informationQUESTION ANSWER. . e. Fourier number:
QUESTION 1. (0 pts) The Lumped Capacitance Method (a) List and describe the implications of the two major assumptions of the lumped capacitance method. (6 pts) (b) Define the Biot number by equations and
More informationChapter 10: Steady Heat Conduction
Chapter 0: Steady Heat Conduction In thermodynamics, we considered the amount of heat transfer as a system undergoes a process from one equilibrium state to another hermodynamics gives no indication of
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
More informationLab 5: Post Processing and Solving Conduction Problems. Objective:
Lab 5: Post Processing and Solving Conduction Problems Objective: The objective of this lab is to use the tools we have developed in MATLAB and SolidWorks to solve conduction heat transfer problems that
More informationGeneral Properties for Determining Power Loss and Efficiency of Passive Multi-Port Microwave Networks
University of Massachusetts Amherst From the SelectedWorks of Ramakrishna Janaswamy 015 General Properties for Determining Power Loss and Efficiency of Passive Multi-Port Microwave Networks Ramakrishna
More informationANALYSIS OF A ONE-DIMENSIONAL FIN USING THE ANALYTIC METHOD AND THE FINITE DIFFERENCE METHOD
J. KSIAM Vol.9, No.1, 91-98, 2005 ANALYSIS OF A ONE-DIMENSIONAL FIN USING THE ANALYTIC METHOD AND THE FINITE DIFFERENCE METHOD Young Min Han* Joo Suk Cho* Hyung Suk Kang** ABSTRACT The straight rectangular
More informationThermal Systems. What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance
Introduction to Heat Transfer What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance Thermal Resistance Thermal Capacitance Thermal
More informationAnalysis of Thermal Behavior of High Frequency Transformers Using Finite Element Method
J. Electromagnetic Analysis & Applications, 010,, 67-63 doi:10.436/emaa.010.1108 ublished Online November 010 (http://www.scirp.org/ournal/emaa) Analysis of Thermal Behavior of High Frequency Transformers
More informationCircle one: School of Mechanical Engineering Purdue University ME315 Heat and Mass Transfer. Exam #1. February 20, 2014
Circle one: Div. 1 (Prof. Choi) Div. 2 (Prof. Xu) School of Mechanical Engineering Purdue University ME315 Heat and Mass Transfer Exam #1 February 20, 2014 Instructions: Write your name on each page Write
More informationUsing Computational Fluid Dynamics And Analysis Of Microchannel Heat Sink
International Journal of Engineering Inventions e-issn: 2278-7461, p-issn: 2319-6491 Volume 4, Issue 12 [Aug. 2015] PP: 67-74 Using Computational Fluid Dynamics And Analysis Of Microchannel Heat Sink M.
More informationThermal aspects of 3D and 2.5D integration
Thermal aspects of 3D and 2.5D integration Herman Oprins Sr. Researcher Thermal Management - imec Co-authors: Vladimir Cherman, Geert Van der Plas, Eric Beyne European 3D Summit 23-25 January 2017 Grenoble,
More informationLinear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More informationMECH 375, Heat Transfer Handout #5: Unsteady Conduction
1 MECH 375, Heat Transfer Handout #5: Unsteady Conduction Amir Maleki, Fall 2018 2 T H I S PA P E R P R O P O S E D A C A N C E R T R E AT M E N T T H AT U S E S N A N O PA R T I - C L E S W I T H T U
More informationThermal Characterization of Packaged RFIC, Modeled vs. Measured Junction to Ambient Thermal Resistance
Thermal Characterization of Packaged RFIC, Modeled vs. Measured Junction to Ambient Thermal Resistance Steven Brinser IBM Microelectronics Abstract Thermal characterization of a semiconductor device is
More informationFinite Element Stress Evaluation Of a Composite Board using 2-D and 3-D Convection Models
Finite Element Stress Evaluation Of a Composite Board using 2-D and 3-D Convection Models Amir Khalilollahii Pennsylvania State University, The Behrend College Abstract Thermal fatigue and high levels
More information