3D Topology Optimization of Heat Sinks for Liquid Cooling
Total Page:16
File Type:pdf, Size:1020Kb
3D Topology Optimization of Heat Sinks for Liquid Cooling Sicheng Sun Piotr Liebersbach Xiaoping Qian Department of Mechanical Engineering University of Wisconsin-Madison TOP Webinar, 02/23/2021 S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 1/19 Introduction f c Breaks the barrier Nowadays the electronics have greater power density and smaller size Conventional heat sinks face the restriction of dissipating more heat while consuming less energy Topology optimization method combining with additive manufacturing shows the potential to solve the problem S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 2/19 Design problem Symmetric boundary on two sides Parabolic inlet condition Heat flux from bottom S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 3/19 Design problem TdΓsc min J = Γsc (1a) γ dΓ R Γsc sc s.t. ⌫ u +(u u)+ p = ↵(γ)u (1b) r·R r ·r r − u = 0 (1c) r⇢·c χ(γ) (u T )= ((γ) T ) (1d) f ·r r ·r 1 µ u 2d⌦+ ↵(γ) u 2d⌦ Φ. (1e) 2 |r | | | 1 Z⌦ Z⌦1 q(1 γ) where ↵(γ)=↵ +(↵ ↵ ) − ,with↵ =0and min max − min q + γ min 2 ↵max = ⌫/(DaL ) S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 4/19 Optimization process Initial γ PDE filter Heaviside filter Solve NS equation Solve the transport equation Solve the adjoint transport equation Solve the adjoint NS equation Compute filtered sensitivity Compute sensitivity Update γ No, update γ Converged Yes, output the result S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 5/19 Optimization process Navier-Stokes equations are discretized with equal order elements with SUPG/PSPG stabilization Continuous consistent adjoint equations are used State and adjoint Navier-Stokes equations are solved with PCD(pressure convection-di↵usion) preconditioner Solvers are implemented with open source finite element software FEniCS 96 CPUs, 6.4 million tetrahedron cells Approx. speed: 20-60 minutes for each optimization iteration, 2 weeks for the total optimization S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 6/19 TO result Inlet Symmetric wall Outlet 3D TO design and geometric model S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 7/19 TO result animation by animate[2019/08/20] S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 8/19 TO result Design viewed from the left side and cut from the middle plane S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 9/19 Cooling mechanism: flow split Disconnected piece Connected free-form fins Staggered fins Hot flow from upstream is sent to top layer Cold flow from upstream is sent to bottom layer S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 10 / 19 Cooling mechanism: flow split Disconnected piece Region A Region B Groove Hot flow goes up Hot flow goes up Middle fin Hot flow from upstream is sent to top layer Cold flow from upstream is sent to bottom layer S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 11 / 19 Cooling mechanism: re-initialization Re-initialization Re-initialization Hot flow from bottom Cold flow to bottom Cold flow 2D temperature plot shows the re-initialization e↵ect Flow split e↵ect is also shown in the 2D figure S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 12 / 19 TO design advantages f c 2.5 0.5 Compared with the TO design, the optimized PF heat sink on the Pareto front costs 40% high pressure drop when the temperature is same raises the temperature by 10% when the pressure drop is same S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 13 / 19 TO design advantages When increasing the flow rate on the inlet The temperature of the optimal PF heat sink is 10% to 40% higher than the TO design The pumping power of the optimal PF heat sink is 200% to 450% higher than the TO design S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 14 / 19 Surface area Surface area measured in Fluent Surface area: per 0.01m width Parallel plate: 1.6 10 3m2 ⇤ − TO bottom structure: 0.57 10 3m2 ⇤ − TO top structure: 0.08 10 3m2 ⇤ − TO design show less surface area than the PF heat sink TO design beat the PF heat sink by better convection This is in contrast to optimization without considering the Navier-Stokes equation, which usually gets tree-type structures by increasing the surface area S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 15 / 19 Conclusion Topology optimization method is a viable and e↵ective computational design tool for heat sink design Compare our topologically optimized design and the optimized conventional plate fin(PF) heat sink: Optimized PF heat sink has 10% higher temperature than TO design when pressure drop is same Optimized PF heat sink cost 40% higher pressure drop than TO design when temperature is same Interesting cooling mechanisms automatically emerge from TO and lead to better heat dissipation Sun, S., Liebersbach, P., & Qian, X. (2020). 3D topology optimization of heat sinks for liquid cooling. Applied Thermal Engineering, 178, 115540. S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 16 / 19 Appendix: heat sink experiments DLMS printed TO design Outside shell TO design and shell assembled S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 17 / 19 Appendix: static mixers 0.5 Γout Γin1 1.0 2.5 Γin2 1.0 x 1.5 z y 10 million tetrahedral cells 32 nodes, 2048 CPUs S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 18 / 19 Appendix: static mixers With 2.5 times pressure drop than the open channel, the mixing performance at the outlet is improved by 90% S. Sun, P. Liebersbach, X. Qian 3D TO of Heat Sinks TOP Webinar, 02/23/2021 19 / 19.