skip to main content

Title: Level-Set-Based Shape & Topology Optimization of Thermal Cloaks
Thermal metamaterials are gaining increasing popularity, especially for heat flux manipulation purposes. However, due to the high anisotropy of the structures resulting from the transformation thermotics or scattering cancellation methods, researchers are resorting to topology optimization as an alternative to find the optimal distribution of constituent bulk materials to realize a specific thermal function. This paper proposes to design a thermal cloak using the level-set-based shape and topology optimization. The thermal cloak design is considered in the context of pure heat conduction. The cloaking effect is achieved by reproducing the reference temperature field through the optimal distribution of two thermally conductive materials. The structural boundary is evolved by solving the Hamilton-Jacobi equation. The feasibility and validity of the proposed method to design thermal meta-devices with cloaking functionality are demonstrated through two numerical examples. The optimized structures have clear boundaries between constituent materials and do not exhibit thermal anisotropy, making it easier for physical realization. The first example deals with a circular cloaking region as a benchmark design. The robustness of the proposed method against various cloaking regions is illustrated by the second example concerning a human-shaped cloaking area. This work can inspire a broader exploration of the thermal meta-device in more » the heat flux manipulation regime. « less
Award ID(s):
Publication Date:
Journal Name:
ASME Design Engineering Technical Conferences
Sponsoring Org:
National Science Foundation
More Like this
  1. Mechanical cloaks are materials engineered to manipulate the elastic response around objects to make them indistinguishable from their homogeneous surroundings. Typically, methods based on material-parameter transformations are used to design optical, thermal, and electric cloaks. However, they are not applicable in designing mechanical cloaks, since continuum-mechanics equations are not form invariant under general coordinate transformations. As a result, existing design methods for mechanical cloaks have so far been limited to a narrow selection of voids with simple shapes. To address this challenge, we present a systematic, data-driven design approach to create mechanical cloaks composed of aperiodic metamaterials using a large precomputed unit cell database. Our method is flexible to allow the design of cloaks with various boundary conditions, multiple loadings, different shapes and numbers of voids, and different homogeneous surroundings. It enables a concurrent optimization of both topology and properties distribution of the cloak. Compared to conventional fixed-shape solutions, this results in an overall better cloaking performance and offers unparalleled versatility. Experimental measurements on additively manufactured structures further confirm the validity of the proposed approach. Our research illustrates the benefits of data-driven approaches in quickly responding to new design scenarios and resolving the computational challenge associated with multiscale designs ofmore »functional structures. It could be generalized to accommodate other applications that require heterogeneous property distribution, such as soft robots and implants design.« less
  2. In this paper, the authors propose a new dimension reduction method for level-set-based topology optimization of conforming thermal structures on free-form surfaces. Both the Hamilton-Jacobi equation and the Laplace equation, which are the two governing PDEs for boundary evolution and thermal conduction, are transformed from the 3D manifold to the 2D rectangular domain using conformal parameterization. The new method can significantly simplify the computation of topology optimization on a manifold without loss of accuracy. This is achieved due to the fact that the covariant derivatives on the manifold can be represented by the Euclidean gradient operators multiplied by a scalar with the conformal mapping. The original governing equations defined on the 3D manifold can now be properly modified and solved on a 2D domain. The objective function, constraint, and velocity field are also equivalently computed with the FEA on the 2D parameter domain with the properly modified form. In this sense, we are solving a 3D topology optimization problem equivalently on the 2D parameter domain. This reduction in dimension can greatly reduce the computing cost and complexity of the algorithm. The proposed concept is proved through two examples of heat conduction on manifolds.
  3. Generators are considered as the core application of electromagnetic machines, which require high-cost rare-earth-based permanent magnets. The development of generators is moving toward high efficiency and increased environmental friendliness. Minimizing the use of rare earth materials such as magnetic materials under the premise of machine performance emerges as a challenging task. Topology optimization has been promisingly applied to many application areas as a powerful generative design tool. It can identify the optimal distribution of magnetic material in the defined design space. This paper employs the level-set-based topology optimization method to design the permanent magnet for generators. The machine under study is a simplified 2D outer rotor direct-drive wind power generator. The dynamic and static models of this generator are studied, and the magnetostatic system is adopted to conduct the topology optimization. The optimization goals in this study mainly focused on two aspects, namely the maximization of the system magnetic energy and the generation of a target magnetic field in the region of the air gap. The continuum shape sensitivity analysis is derived by using the material time derivative, the Lagrange multiplier method, and the adjoint variable method. Two numerical examples are investigated, and the effectiveness of the proposed design frameworkmore »is validated by comparing the performance of the original design against the optimized design.« less
  4. Abstract

    The goal of this work is to predict the effect of part geometry and process parameters on the instantaneous spatial distribution of heat, called the heat flux or thermal history, in metal parts as they are being built layer-by-layer using additive manufacturing (AM) processes. In pursuit of this goal, the objective of this work is to develop and verify a graph theory-based approach for predicting the heat flux in metal AM parts. This objective is consequential to overcome the current poor process consistency and part quality in AM. One of the main reasons for poor part quality in metal AM processes is ascribed to the heat flux in the part. For instance, constrained heat flux because of ill-considered part design leads to defects, such as warping and thermal stress-induced cracking. Existing non-proprietary approaches to predict the heat flux in AM at the part-level predominantly use mesh-based finite element analyses that are computationally tortuous — the simulation of a few layers typically requires several hours, if not days. Hence, to alleviate these challenges in metal AM processes, there is a need for efficient computational thermal models to predict the heat flux, and thereby guide part design and selection of processmore »parameters instead of expensive empirical testing. Compared to finite element analysis techniques, the proposed mesh-free graph theory-based approach facilitates layer-by-layer simulation of the heat flux within a few minutes on a desktop computer. To explore these assertions we conducted the following two studies: (1) comparing the heat diffusion trends predicted using the graph theory approach, with finite element analysis and analytical heat transfer calculations based on Green’s functions for an elementary cuboid geometry which is subjected to an impulse heat input in a certain part of its volume, and (2) simulating the layer-by-layer deposition of three part geometries in a laser powder bed fusion metal AM process with: (a) Goldak’s moving heat source finite element method, (b) the proposed graph theory approach, and (c) further comparing the heat flux predictions from the last two approaches with a commercial solution. From the first study we report that the heat flux trend approximated by the graph theory approach is found to be accurate within 5% of the Green’s functions-based analytical solution (in terms of the symmetric mean absolute percentage error). Results from the second study show that the heat flux trends predicted for the AM parts using graph theory approach agrees with finite element analysis with error less than 15%. More pertinently, the computational time for predicting the heat flux was significantly reduced with graph theory, for instance, in one of the AM case studies the time taken to predict the heat flux in a part was less than 3 minutes using the graph theory approach compared to over 3 hours with finite element analysis. While this paper is restricted to theoretical development and verification of the graph theory approach for heat flux prediction, our forthcoming research will focus on experimental validation through in-process sensor-based heat flux measurements.

    « less
  5. Abstract Heterogeneous and complex electronic packages may require unique thermomechanical structures to provide optimal heat guiding. In particular, when a heat source and a heat sink are not aligned and do not allow a direct path, conventional thermal management methods providing uniform heat dissipation may not be appropriate. Here we present a topology optimization method to find thermally conductive and mechanically stable structures for optimal heat guiding under various heat source-sink arrangements. To exploit the capabilities, we consider complex heat guiding scenarios and three-dimensional (3D) serpentine structures to carry the heat with corner angles ranging from 30 deg to 90 deg. While the thermal objective function is defined to minimize the temperature gradient, the mechanical objective function is defined to maximize the stiffness with a volume constraint. Our simulations show that the optimized structures can have a thermal resistance of less than 32% and stiffness greater than 43% compared to reference structures with no topology optimization at an identical volume fraction. The significant difference in thermal resistance is attributed to a thermally dead volume near the sharp corners. As a proof-of-concept experiment, we have created 3D heat guiding structures using a selective laser melting technique and characterized their thermal propertiesmore »using an infrared thermography technique. The experiment shows the thermal resistance of the thermally optimized structure is 29% less than that of the reference structure. These results present the unique capabilities of topology optimization and 3D manufacturing in enabling optimal heat guiding for heterogeneous systems and advancing the state-of-the-art in electronics packaging.« less