Level set topology optimization of cooling and heating devices using a simplified convection model

This paper studies topology optimization of convective heat transfer problems in two and three dimensions. The convective fluxes are approximated by Newton's Law of Cooling (NLC). The geometry is described by a Level Set Method (LSM) and the temperature field is predicted by the eXtended Finite Element Method (XFEM). A constraint on the spatial gradient of the level set field is introduced to penalize small, sub-element-size geometric features. Numerical studies show that the LSM-XFEM provides improved accuracy over previously studied density methods and LSMs using Ersatz material models. It is shown that the NLC model with an iso-thermal fluid phase may over predict the convective heat flux and thus promote the formation of very thin fluid channels, depending on the Biot number characterizing the heat transfer problem. Approximating the temperature field in the fluid phase by a diffusive model mitigates this issue but an explicit feature size control is still necessary to prevent the formation of small solid members, in particular at low Biot numbers. The proposed constraint on the gradient of the level set field is shown to suppress sub-element-size features but necessitates a continuation strategy to prevent the optimization process from stagnating as geometric features merge.