•A conventional level set method is applied to design thermal cloaks in the Euclidean spaces (2D planar surfaces and 3D solid).•A conformal manifold thermal cloak is devised for the first time using ...the extended level set method (X-LSM) incorporating conformal geometry theory.•The topologically optimized thermal cloaks do not exhibit material anisotropy and non-homogeneity.•The robustness and validity of the proposed method are demonstrated for a variety of cloaking regions.
Thermal cloaks are devices designed to shield an object against thermal detection, which have attracted growing interest in research. This paper proposes to design thermal cloaks using the level-set-based shape and topology optimization in the context of pure heat conduction. The cloaking effect is achieved by optimizing the distribution of two bulk heat conductive materials to eliminate the temperature disturbance induced by the introduction of the insulator (cloaking region) into a homogeneous thermal conduction medium. The optimized thermal cloaks are free of high anisotropy and nonhomogeneity commonly seen in the popular transformation thermotics or scattering cancellation methods. Due to the clear boundary characteristic of the level set representation, no sophisticated filtering techniques are required to suppress the appearance of ”gray regions” as opposed to the density-based topology optimization methods. Considering the fact that the device components that need to be thermally cloaked, e.g., sensors, can take an arbitrary free-form shape, a conformal thermal cloak on the manifold is also topologically optimized using the extended level set method (X-LSM), which has not been reported in the literature. The structural boundary is evolved by solving the (modified) Hamilton-Jacobi equation. The feasibility and robustness of the proposed method to design thermal meta-devices with cloaking functionality are demonstrated through a number of 2D and 3D (solid and shell) numerical examples with different cloaking regions (circular, human-shaped, spherical, and curved circular). This work may shed light on further exploration of the thermal meta-devices in the heat flux manipulation regime.
Passive heat sinks cooled by natural convection are reliable, compact, and low‐noise. They are widely used in telecommunication devices, LEDs, and so forth. This work builds upon the recent ...advancements in fluid topology optimization (TO) to present a case study of two‐ and three‐dimensional optimum design and thermal modeling for the natural convection problems using a reaction–diffusion equation (RDE)‐based level‐set method. To this end, first, a high‐fidelity thermal‐fluid model is constructed where the full Navier–Stokes equations are strongly coupled with the energy equation through the Boussinesq approximation. We benchmark our simulation solver against the experimental analysis and other numerical analysis methods. Next, we carefully investigate the flow behavior under different Grashof numbers using a fully transient simulation solver. Then, we revisit the RDE‐based level‐set TO methodology and construct an open‐source parallel TO framework. The main findings reveal that using the body‐fitted mesh adaptation, the proposed methodology can capture the explicit fluid–solid boundary and we are free of the continuation approach to penalize the design variable to the binary structure. A moderately large‐scale TO problem with 3.56×106 DOFs can be solved in parallel on a standard multi‐process platform. Our numerical implementation uses FreeFEM for finite element analysis, PETSc for distributed linear algebra, and Mmg for mesh adaptation. For comparison and for accessing our various techniques, a variety of 2D and 3D benchmarks are presented to support these remarkable features.
Intensity inhomogeneity often appears in medical images and causes great difficulties in image segmentation. Most active contour models perform poorly when applied to intensity inhomogeneous images ...because their energy functions use local intensity information in a fixed-size domain, causing the contour to evolve in the wrong direction. To overcome the difficulties caused by intensity inhomogeneity, we propose an adaptive multi-scale Gaussian kernel function based on image local entropy, which can determine the appropriate scale for each local region. Choosing the small scale and large scale for inhomogeneous and homogeneous areas respectively make the contour move toward the target boundary accurately. We also propose three adaptive multi-scale (AMS) models, AMS-region scalable fitting (AMS-RSF) model, AMS-local image fitting (AMS-LIF) model, AMS-local and global intensity fitting (AMS-LGIF) model, to segment medical images with intensity inhomogeneity and noise, including left atrial MR images and breast ultrasound images. The experimental results show that the adaptive multi-scale Gaussian kernel function enables the active contour model to effectively segment intensity inhomogeneous images and has a certain robustness to the initial contour and noise, which achieves good performance on MR left atrial images and ultrasound images of breast cancer. The AMS-LGIF model obtained the highest DICE coefficient of 0.9532, which was better than the 0.9429 obtained by the second-ranked LGIF model to segment left atrial MR images. For segmenting breast ultrasound images, the DICE coefficient is increased by 16% than that of the U-Net++ model.
•An adaptive multi-scale Gaussian kernel function based on image local entropy is proposed.•Adaptive multi-scale RSF, LIF, and LGIF models are proposed.•The proposed adaptive multi-scale model is robust to initial contours and noise.•Numerical results such as DICE, JS, Precision, and Recall values are given.•The experimental results show that our adaptive multi-scale model is better than the other methods.
Additive manufacturing (AM) enables fabrication of multiscale cellular structures as a whole part, whose features can span several dimensional scales. Both the configurations and layout pattern of ...the cellular lattices have great impact on the overall performance of the lattice structure. In this paper, we propose a novel design method to optimize cellular lattice structures to be fabricated by AM. The method enables an optimized load-bearing solution through optimization of geometries of global structures and downscale mesostructures, as well as global distributions of spatially-varying graded mesostructures. A shape metamorphosis technology is incorporated to construct the graded mesostructures with essential interconnections. Experimental testing is undertaken to verify the superior stiffness properties of the optimized graded lattice structure compared to the baseline design with uniform mesostructures.
Display omitted
•A design method is developed for graded lattice structures.•Geometries of global structure, lattice meso-structures and their spatial layout are optimized.•Experimental testing shows 24% improved stiffness of optimized graded lattice structures than uniform ones.
Summary
In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented ...algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid.
We develop an efficient solver for the direct numerical simulations of viscoelastic and elastoviscoplastic multiphase flows. The solver is validated for various benchmark single‐phase and two‐phase elastoviscoplastic flow problems. We study a Newtonian droplet rising in an elastoviscoplastic fluid and show the appearance of a negative wake.