24-26 September 2008, Rome, ItalyInitial attempts to achieve the goals set by MACE werewith constant diameter staggered pin fins. This was notsuccessful. Another parameter was introduced by allowing thepins to be tapered. The results are shown in Figs. 5 and 6,Figure 6. Pin Fin Geometry with pin tip diameter of0.958mm and pitch to pin base diameter of 2.0A Fin base ThicknessB Fin tip thicknessC Pitch/Fin base thicknessE Reynolds numberFigure 7. Model Term Ranking Pareto Chart For Plane Finsoutput from Fusion Pro. Figure 5 is a Pareto Diagram showingthe relative importance of the different parameters. Theparameter that is conspicuously absent is the base thickness. Itis only hwne it becomes very thin that it becomes important.Figure 6 is the trajectory of the temperature at the midplane ofthe heatsink (TMID), the temperature at the heatsink exit(TMAX) and the effectiveness as a function of pin basediameter and Reynolds number. These results were achievedwith a staggered array of pins having a pin base diameter of2.2mm, a pin tip diameter of 1.0mm and a pitch to pin basediameter of 2.0 at a Reynolds number of 43,000. This result isnot fully optimized because the author is not yet skilledenough at the use of DoE. It is, nevertheless, remarkable.A similar effort was made to design an optimum heatsinkusing plane fins. This was not as successful as for the pin fins.The Pareto Diagram. Fig. 7, shows that the first threeparameters are almost equal in importance making the searchfor an optimum even more difficult. An effectiveness of 15was achieved, however, with tip thickness of 0.1mm, pitch tobase diameter of 4.0 and base thickness of 1mm at a Reynoldsnumber ranging from 1150 to 2000, see Fig. 8. Note the strongdependence on the base diameter. The reason is the need formore metal to conduct heat to the fluid-solid interface as theheat transfer coefficient increases with increasing velocity.CONCLUSIONSA preliminary study of micro heat sink optimization wascarried out. The problem was kept manageable by assumingthe morphologies evaluated were relatively simple. The heatflux to the heat sink was assumed to be uniform.The design of a micro channel heat sink for 4 inch by 4inch base was optimized for a uniform 1000W heat load withtwo constraints; first, the midpoint of the base temperaturemust be below 65 0 C and the overall heat sink height is limitedto 1 inch. Statistical design of experiment was used to find theoptimum performance variables. The VAT governingequations were solved with the proper closure to describe fluidflow and heat transfer in micro channels in lieu of experiment.DOE is shown to be a useful tool in both achieving an optimaldesign and in determining the importance of the geometricparameters describing the micro channel.Table 1. Optimal DesignsStaggered Pin Fins Plane finsBase Dimension 2.2 mm 0.2 mmTip Dimension 1.0 mm 0.1 mmPitch/BaseDimension2.0 4.0Reynolds Number 41,000 1000EffectivenessQ/PP50 20Mass Flow 0.027CFM 57.40Pressure Drop -13300The major conclusion that can be drawn from this study isthat optimization is difficult to accomplish. With each value ofthe Reynolds number having a different "best geometry",©<strong>EDA</strong> <strong>Publishing</strong>/THERMINIC 2008 181ISBN: 978-2-35500-008-9
24-26 September 2008, Rome, ItalyFigure 8. Plane Fin Geometry with a fin base of 0.2mm, atip of 0.1 mm and a pitch to base dimension of 4.0another approach should be taken. An optimum configurationshould have been sought for a given mass flow rate. Theconfigurations that yielded a temperature that was withinallowable limits could then be evaluated to see which wasmost suitable based on the maximum effectiveness.Nevertheless, the results demonstrate that significantimprovements can be made by solving the conjugate problem.The axially uniform morphology and heat fluxassumption made in this work is easily relaxed so that themorphology variation can become part of the optimizationprocess. This would enable one to try and make the heattransfer coefficient high where temperature differences aresmall leading to an even more optimal configuration. It shouldbe noted that the staggered pin fin may not be the best whenthe heat flux to the heat sink is non-uniform. Both morphologyand its variation would be different.As far as we know this is the first complete optimizationof a micro heat sink.ACKNOWLEDGEMENTSThe support of a Department of Energy NERI grant, Award NumberDE-FC07-07ID14827,is gratefully acknowledged. The CFDcomputations were performed using a commercial version ofSC/Tetra.Hwang, G., J., and Cao, C., H., 1994``Heat Transfer Measurementand Analysis for sintered porous channels.``, Journal of HeatTransfer, Vol.116, pp.456-464.Jiang, P., Si, G., Li, M. Ren, Z., 2004, ``Experimental andNumerical Investigation of Forced Convection Heat Transfer of Airin Non-sintered Porous Media``, Experimental Thermal and FluidScience, Vol.28, pp.545-555.Kawano, K., Minakami, K., Iwasaki, H. and Ishizuka, M., 1998,``Micro Channel Heat Exchanger for Cooling Electrical Equipment``,in Proceedings of ASME Heat Transfer Division – volume 3, HTDvol.361-3/PID-vol.3, pp.173-180.Koşar, A., Peles, Y ,2006, ``Thermal-Hydraulic Performance ofMEMS-based Pin Fin Heat Sink`` Journal of Heat Transfer, Vol.128, No.2, pp. 121-131.Lin, Y., Y., Semenic, T., Catton, I., 2005,``ThermophysicalProperties of Biporous Sintered Copper``, ASME InternationalMechanical Engineering Congress and Exposition, Orlando, Florida,November 5-11.Peles, Y., Kou, C., Koşar, A., Mishra, C. and Schneider, B.,``Forced Convective Heat Transfer Across a Pin Fin Micro HeatEchanger``, Int. J. Heat Mass Transfer, Vol.48, No.17, pp.3615-3627.Rizzi, M., and Catton, Ivan, 2002, ``Experimental Results forEndwall and Pin Fin Heat Transfer Coefficients``, Proceedings of the12th INHTC, 2002Travkin, V.S. and I. Catton, 1999a, ''Compact Heat ExchangerOptimization Tools Based on Volume Averaging Theory,'' in Proc.33rd ASME NHTC, NHTC99-246. ASME, New Mexico.Travkin, V.S. and Catton, I., (1998), "Porous Media TransportDescriptions - Non-Local, Linear and Nonlinear Against EffectiveThermal/Fluid Properties", Advances in Colloid and InterfaceScience, Vol. 76-77, pp. 389-443.Travkin, V.S. and I. Catton, 1999b, “Turbulent Flow and HeatTransfer Modeling in a Flat Channel with Regular Highly RoughWalls,” International Journal of Fluid Mechanics Research, Vol 26,No. 2, pp 110-135.Travkin, V.S. and Catton, I. (2001), "Transport Phenomena inHeterogeneous Media Based on Volume Averaging Theory",Advances in Heat Transfer, Vol. 34, pp.1-144.Wakao, N., Kaguei, S, 1982, Heat and Mass Transfer in PackedBeds, Taylor & FrancisWatanabe H, 1989, `` Drag Coefficient And Voidage Function OnFluid-Flow Through Granular Packed-Beds ``, International JournalOf Engineering Fluid Mechanics, Vol.2, No.1:, pp.93-108.REFERENCESBejan, A. and Morega, A.M. (1993), "Optimal Arrays of Pin Finsand Plate Fins in Laminar Forced Convection," Journal of HeatTransfer, Vol. 115, pp. 75-81.Bird, R., B., Stewart, W., E., Lightfoot E., N., 2001 TransportPhenomena, 2nd ed., Wiley, pp.178-179.Catton, I. And Hu, K., 2003, ``VAT Based Optimization of HeatTransfer in a Flat Channel Filled with a Porous Media``, inProceedings of Summer Heat Transfer Conference Summer HeatTransfer conference Las Vegas, Nevada, USAGratton, L., Travkin, V.S., and Catton, I. (1996), "The Influenceof Morphology upon Two- Temperature Statements for ConvectiveTransport in Porous Media," Journal of Enhanced Heat Transfer,Vol. 3, No. 2, pp.129-145.©<strong>EDA</strong> <strong>Publishing</strong>/THERMINIC 2008 182ISBN: 978-2-35500-008-9
- Page 1 and 2:
http://cmp.imag.fr/conferences/ther
- Page 3 and 4:
24-26 September 2008, Rome, Italy©
- Page 5 and 6:
24-26 September 2008, Rome, Italy©
- Page 7 and 8:
24-26 September 2008, Rome, Italy©
- Page 9 and 10:
24-26 September 2008, Rome, ItalyBL
- Page 11 and 12:
24-26 September 2008, Rome, ItalySe
- Page 13 and 14:
24-26 September 2008, Rome, ItalyED
- Page 15 and 16:
Staggered Cartesian meshes tolerate
- Page 17 and 18:
VIII. THE CASE FOR MCAD-EMBEDDED EC
- Page 19 and 20:
24-26 September 2008, Rome, ItalyTr
- Page 21 and 22:
solution to the eikonal equation(
- Page 23 and 24:
24-26 September 2008, Rome, ItalyS
- Page 25 and 26:
24-26 September 2008, Rome, ItalyTr
- Page 27 and 28:
24-26 September 2008, Rome, Italyda
- Page 29 and 30:
q (RthΣ, t)q (R th , t )10.80.60.4
- Page 31 and 32:
24-26 September 2008, Rome, ItalyEv
- Page 33 and 34:
24-26 September 2008, Rome, Italy*R
- Page 35 and 36:
24-26 September 2008, Rome, Italy[7
- Page 37 and 38:
24-26 September 2008, Rome, ItalyNe
- Page 39 and 40:
24-26 September 2008, Rome, ItalyS=
- Page 41 and 42:
24-26 September 2008, Rome, ItalyVI
- Page 43 and 44:
24-26 September 2008, Rome, Italyth
- Page 45 and 46:
24-26 September 2008, Rome, Italypr
- Page 47 and 48:
24-26 September 2008, Rome, ItalyFi
- Page 49 and 50:
24-26 September 2008, Rome, ItalyPo
- Page 51 and 52:
24-26 September 2008, Rome, ItalyWL
- Page 53 and 54:
24-26 September 2008, Rome, ItalyTa
- Page 55 and 56:
24-26 September 2008, Rome, ItalyII
- Page 57 and 58:
Thermal resistance (K/W)0.250.200.1
- Page 59 and 60:
24-26 September 2008, Rome, ItalyB.
- Page 61 and 62:
temperatrue rise [K]76543210Sim0,00
- Page 63 and 64:
24-26 September 2008, Rome, ItalyCo
- Page 65 and 66:
24-26 September 2008, Rome, ItalyPr
- Page 67 and 68:
24-26 September 2008, Rome, ItalyPr
- Page 69 and 70:
24-26 September 2008, Rome, ItalyBl
- Page 71 and 72:
conclusion, every tile is contacted
- Page 73 and 74:
and B10), three space blocks (W4, W
- Page 75 and 76:
24-26 September 2008, Rome, ItalyMu
- Page 77 and 78:
After having determined the k eff v
- Page 79 and 80:
24-26 September 2008, Rome, ItalyVI
- Page 81 and 82:
24-26 September 2008, Rome, ItalyTh
- Page 83 and 84:
PTjunc24-26 September 2008, Rome, I
- Page 85 and 86:
tures - has been shown previously [
- Page 87 and 88:
24-26 September 2008, Rome, ItalyAu
- Page 89 and 90:
IV.CONCLUSIONSMentor Graphics Exped
- Page 91 and 92:
24-26 September 2008, Rome, ItalyIn
- Page 93 and 94:
III. HEDORIS SIMULATIONSThermal sim
- Page 95 and 96:
24-26 September 2008, Rome, Italytr
- Page 97 and 98:
24-26 September 2008, Rome, ItalyDe
- Page 99 and 100:
24-26 September 2008, Rome, ItalyII
- Page 101 and 102:
datasheets, these values are usuall
- Page 103 and 104:
24-26 September 2008, Rome, Italyte
- Page 105 and 106:
24-26 September 2008, Rome, ItalyCo
- Page 107 and 108:
24-26 September 2008, Rome, Italy20
- Page 109 and 110:
material to handle difficulties suc
- Page 111 and 112:
24-26 September 2008, Rome, ItalyAt
- Page 113 and 114:
24-26 September 2008, Rome, Italymo
- Page 115 and 116:
heatspreaders deposited during the
- Page 117 and 118:
24-26 September 2008, Rome, ItalyTh
- Page 119 and 120:
24-26 September 2008, Rome, Italyla
- Page 121 and 122:
Fig. 13 Transient of the output dio
- Page 123 and 124:
24-26 September 2008, Rome, ItalyIn
- Page 125 and 126:
chamber to calibrate the thermal te
- Page 127 and 128:
V.Table 2: FE-Simulation test-matri
- Page 129 and 130:
24-26 September 2008, Rome, ItalyCo
- Page 131 and 132:
24-26 September 2008, Rome, Italy(c
- Page 133 and 134:
Case coverSupercapacitor n°2Electr
- Page 135 and 136:
24-26 September 2008, Rome, Italyfr
- Page 137 and 138:
24-26 September 2008, Rome, ItalyTA
- Page 139 and 140:
24-26 September 2008, Rome, ItalyTh
- Page 141 and 142: 24-26 September 2008, Rome, Italych
- Page 143 and 144: 24-26 September 2008, Rome, ItalyDe
- Page 145 and 146: Peltier control units (PCUs). The t
- Page 147 and 148: 24-26 September 2008, Rome, ItalyTh
- Page 149 and 150: every step: the fewer cells the mor
- Page 151 and 152: TABLE IIEffect of multiplication-bl
- Page 153 and 154: 24-26 September 2008, Rome, ItalyTh
- Page 155 and 156: 24-26 September 2008, Rome, ItalyRe
- Page 157 and 158: Drain source leakage current (A)Rev
- Page 159 and 160: 24-26 September 2008, Rome, ItalyFP
- Page 161 and 162: V IN =V DD /2Fig. 5. Short circuit
- Page 163 and 164: 24-26 September 2008, Rome, ItalyTh
- Page 165 and 166: The overall objective of the NANOPA
- Page 167 and 168: 24-26 September 2008, Rome, ItalyRe
- Page 169 and 170: predictions show that carbon nanotu
- Page 171 and 172: 24-26 September 2008, Rome, ItalyFl
- Page 173 and 174: 55. Campbell, R.C., S.E. Smith, and
- Page 175 and 176: 24-26 September 2008, Rome, ItalyCo
- Page 177 and 178: 24-26 September 2008, Rome, ItalyTh
- Page 179 and 180: 24-26 September 2008, Rome, ItalyEX
- Page 181 and 182: 24-26 September 2008, Rome, Italyen
- Page 183 and 184: 24-26 September 2008, Rome, ItalyKE
- Page 185 and 186: An original three-step etch process
- Page 187 and 188: with a network of built-in sensors
- Page 189 and 190: 24-26 September 2008, Rome, Italyan
- Page 191: 24-26 September 2008, Rome, ItalyFu
- Page 195 and 196: 24-26 September 2008, Rome, ItalyFi
- Page 197 and 198: 24-26 September 2008, Rome, ItalyWe
- Page 199 and 200: semiconductors is of the order of 1
- Page 201 and 202: 24-26 September 2008, Rome, ItalyDu
- Page 203 and 204: transistor, through the bulk and be
- Page 205 and 206: impedance behaviour and on the othe
- Page 207 and 208: Tungstenmicro-heaterGas SensingMate
- Page 209 and 210: Tungsten micro-heater resistance (
- Page 211 and 212: 24-26 September 2008, Rome, ItalyPo
- Page 213 and 214: 24-26 September 2008, Rome, ItalyT3
- Page 215 and 216: 24-26 September 2008, Rome, ItalyEv
- Page 217 and 218: 24-26 September 2008, Rome, ItalyD.
- Page 219 and 220: 24-26 September 2008, Rome, ItalyOn
- Page 221 and 222: 24-26 September 2008, Rome, ItalyAc
- Page 223 and 224: 24-26 September 2008, Rome, Italy
- Page 225 and 226: 24-26 September 2008, Rome, ItalyLO
- Page 227 and 228: 24-26 September 2008, Rome, Italymu
- Page 229 and 230: 24-26 September 2008, Rome, Italy(R
- Page 231 and 232: 24-26 September 2008, Rome, ItalyPr
- Page 233 and 234: 24-26 September 2008, Rome, ItalyPa
- Page 235 and 236: 24-26 September 2008, Rome, ItalyEl
- Page 237 and 238: Temperature rise (°C)1,81,61,41,21
- Page 239 and 240: i résistances carbone (A)i résist
- Page 241 and 242: calculates the dissipation distribu
- Page 243 and 244:
in the next. It can also be used to
- Page 245 and 246:
allows creating special methods of
- Page 247 and 248:
OLED device (see Fig. 2.) provided
- Page 249 and 250:
24-26 September 2008, Rome, ItalyFi
- Page 251 and 252:
24-26 September 2008, Rome, ItalyIn
- Page 253:
24-26 September 2008, Rome, ItalyWa