Numerical approximation of phase field based shape and topology optimization for fluids

We consider the problem of finding optimal shapes of fluid domains. The fluid obeys the Navier--Stokes equations. Inside a holdall container we use a phase field approach using diffuse interfaces to describe the domain of free flow. We formulate a corresponding optimization problem where flow outside the fluid domain is penalized. The resulting formulation of the shape optimization problem is shown to be well-posed, hence there exists a minimizer, and first order optimality conditions are derived. For the numerical realization we introduce a mass conserving gradient flow and obtain a Cahn--Hilliard type system, which is integrated numerically using the finite element method. An adaptive concept using reliable, residual based error estimation is exploited for the resolution of the spatial mesh. The overall concept is numerically investigated and comparison values are provided

Topology optimization of multiple anisotropic materials, with application to self-assembling diblock copolymers

We propose a solution strategy for a multimaterial minimum compliance topology optimization problem, which consists in finding the optimal allocation of a finite number of candidate (possibly anisotropic) materials inside a reference domain, with the aim of maximizing the stiffness of the body. As a relevant and novel application we consider the optimization of self-assembled structures obtained by means of diblock copolymers. Such polymers are a class of self-assembling materials which spontaneously synthesize periodic microstructures at the nanoscale, whose anisotropic features can be exploited to build structures with optimal elastic response, resembling biological tissues exhibiting microstructures, such as bones and wood. For this purpose we present a new generalization of the classical Optimality Criteria algorithm to encompass a wider class of problems, where multiple candidate materials are considered, the orientation of the anisotropic materials is optimized, and the elastic properties of the materials are assumed to depend on a scalar parameter, which is optimized simultaneously to the material allocation and orientation. Well-posedness of the optimization problem and well-definition of the presented algorithm are narrowly treated and proved. The capabilities of the proposed method are assessed through several numerical tests

Non-local control in the conduction coefficients: well posedness and convergence to the local limit

We consider a problem of optimal distribution of conductivities in a system governed by a non-local diffusion law. The problem stems from applications in optimal design and more specifically topology optimization. We propose a novel parametrization of non-local material properties. With this parametrization the non-local diffusion law in the limit of vanishing non-local interaction horizons converges to the famous and ubiquitously used generalized Laplacian with SIMP (Solid Isotropic Material with Penalization) material model. The optimal control problem for the limiting local model is typically ill-posed and does not attain its infimum without additional regularization. Surprisingly, its non-local counterpart attains its global minima in many practical situations, as we demonstrate in this work. In spite of this qualitatively different behaviour, we are able to partially characterize the relationship between the non-local and the local optimal control problems. We also complement our theoretical findings with numerical examples, which illustrate the viability of our approach to optimal design practitioners

Revisiting topology optimization with buckling constraints

We review some features of topology optimization with a lower bound on the critical load factor, as computed by linearized buckling analysis. The change of the optimized design, the competition between stiffness and stability requirements and the activation of several buckling modes, depending on the value of such lower bound, are studied. We also discuss some specific issues which are of particular interest for this problem, as the use of non-conforming finite elements for the analysis, the use of inconsistent sensitivities in the optimization and the replacement of the single eigenvalue constraints with an aggregated measure. We discuss the influence of these practices on the optimization result, giving some recommendations.Comment: 15 pages, 12 figures, 2 table