Shape Optimization and Fictitious Domain Approach for Solving Free Boundary Problems of Bernoulli Type

This contribution deals with an efficient method for the numerical realization of the exterior and interior Bernoulli free boundary problems. It is based on a shape optimization approach. The state problems are solved by a fictitious domain solver using boundary Lagrange multipliers.     

Analytical and numerical methods in shape optimization

This paper is intended to overview on analytical and numerical methods in shape optimization. We compute and analyse the shape Hessian in order to distinguish well‐posed and ill‐posed shape optimization problems. We introduce different discretization techniques of the shape and present existence and convergence results of approximate solutions in case of well posedness. Finally, we survey on the efficient numerical solution of the state equation, including finite and boundary element methods as well as fictitious domain methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Shape and topology optimization of contact problems by level set methods

This paper deals with the numerical solution of a topology and shape optimization problems of an elastic body in unilateral contact with a rigid foundation. The contact problem with the prescribed friction is considered. The structural optimization problem consists in finding such shape of the boundary of the domain occupied by the body that the normal contact stress along the contact boundary of the body is minimized. In the paper shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization. Level set based numerical algorithm for the solution of the shape optimization problem is proposed. Level set method is used to describe the position of the boundary of the body and its evolution on a fixed mesh. This evolution is governed by Hamilton – Jacobi equation. The speed vector field driving the propagation of the boundary of the body is given by the shape derivative of a cost functional with respect to the free boundary. Numerical examples are provided. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Efficient treatment of stationary free boundary problems

In the present paper we consider the efficient treatment of free boundary problems by shape optimization. We reformulate the free boundary problem as shape optimization problem. A second order shape calculus enables us to analyze the shape problem under consideration and to prove convergence of a Ritz–Galerkin approximation of the shape. We show that Newton's method requires only access to the underlying state function on the boundary of the domain. We compute these data by boundary integral equations which are numerically solved by a fast wavelet Galerkin scheme. Numerical results prove that we succeeded in finding a fast and robust algorithm for solving the considered class of problems.     

Shape optimization for Stokes problem with threshold slip

We study the Stokes problems in a bounded planar domain Ω with a friction type boundary condition that switches between a slip and no-slip stage. Our main goal is to determine under which conditions concerning the smoothness of Ω solutions to the Stokes system with the slip boundary conditions depend continuously on variations of Ω. Having this result at our disposal, we easily prove the existence of a solution to optimal shape design problems for a large class of cost functionals. In order to release the impermeability condition, whose numerical treatment could be troublesome, we use a penalty approach. We introduce a family of shape optimization problems with the penalized state relations. Finally we establish convergence properties between solutions to the original and modified shape optimization problems when the penalty parameter tends to zero.     

A new fictitious domain method in shape optimization

The present paper is concerned with investigating the capability of the smoothness preserving fictitious domain method from Mommer (IMA J. Numer. Anal. 26:503–524, 2006) to shape optimization problems. We consider the problem of maximizing the Dirichlet energy functional in the class of all simply connected domains with fixed volume, where the state equation involves an elliptic second order differential operator with non-constant coefficients. Numerical experiments in two dimensions validate that we arrive at a fast and robust algorithm for the solution of the considered class of problems. The proposed method can be applied to three dimensional shape optimization problems.     

Hybrid optimization procedure applied to optimal location finding for piezoelectric actuators and sensors for active vibration control

This paper presents an efficient hybrid optimization approach using a new coupling technique for solving the constrained optimization problems. This methodology is based on genetic algorithm, sequential quadratic programming and particle swarm optimization combined with a projected gradient techniques in order to correct the solutions out of domain and send them to the domain’s border. The established procedures have been successfully tested with some well known mathematical and engineering optimization problems, also the obtained results are compared with the existing approaches. It is clearly demonstrated that the solutions obtained by the proposed approach are superior to those of existing best solutions reported in the literature. The main application of this procedure is the location optimization of piezoelectric sensors and actuators for active control, the vibration of plates with some piezoelectric patches is considered. Optimization criteria ensuring good observability and controllability based on some main eigenmodes and residual ones are considered. Various rectangular piezoelectric actuators and sensors are used and two optimization variables are considered for each piezoelectric device: the location of its center and shape orientation. The applicability and effectiveness of the present methodological approach are demonstrated and the location optimization of multiple sensors and actuators are successfully obtained with some main modes and residual ones. The shape orientation optimization of sensors observing various modes as well as the local optimization of multiple sensors and actuators are numerically investigated. The effect of residual modes and the spillover reduction can be easily analyzed for a large number of modes and multiple actuators and sensors.     

Shape Optimization for Semi-Linear Elliptic Equations Based on an Embedding Domain Method

We study a class of shape optimization problems for semi-linear elliptic equations with Dirichlet boundary conditions in smooth domains in ℝ². A part of the boundary of the domain is variable as the graph of a smooth function. The problem is equivalently reformulated on a fixed domain. Continuity of the solution to the state equation with respect to domain variations is shown. This is used to obtain differentiability in the general case, and moreover a useful formula for the gradient of the cost functional in the case where the principal part of the differential operator is the Laplacian. Online publication 23 January 2004.     

An Approach to Shape Optimization with State Constraints

We introduce a new approach to solve shape optimization problems with state constraints. The problem is reformulated on a fixed reference domain using conformal pull-back. The shape dependence is then hidden in the conformal parameter. The problem is discretized using FEM and solved as an NLP. Finally the optimal shape can be reconstructed from the optimal conformal parameter. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Gradients for Shape Optimization Based on Embedding Domain Techniques

An embedding domain technique is proposed to characterize the gradient of shape optimization problems. A discussion of the numerical realization of the arising saddle point problems is given and numerical feasibility of the gradient information is discussed.     

Coupling of FEM and BEM in Shape Optimization

In the present paper we consider the numerical solution of shape optimization problems which arise from shape functionals of integral type over a compact region of the unknown shape, especially L ²-tracking type functionals. The underlying state equation is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the optimization problem under consideration. Since the adjoint state depends directly on the state, we propose a coupling of finite element methods (FEM) and boundary element methods (BEM) to realize an efficient first order shape optimization algorithm. FEM is applied in the compact region while the rest is treated by BEM. The coupling of FEM and BEM essentially retains all the structural and computational advantages of treating the free boundary by boundary integral equations.This research has been carried out when the second author stayed at the Department of Mathematics, Utrecht University, The Netherlands, supported by the EU-IHP project Nonlinear Approximation and Adaptivity: Breaking Complexity in Numerical Modelling and Data Representation     

Non-linear eigenvalue problems for <Emphasis Type="Italic">p</Emphasis>-Laplacian with variable domain

In this work we consider some eigenvalue problems for p-Laplacian with variable domain. Eigenvalues of this operator are taken as a functional of the domain. We calculate the first variation of this functional, using the obtained formula investigate behavior of the eigenvalues when the domain varies. Then we consider one shape optimization problem for the first eigenvalue, prove the necessary condition of optimality relatively domain, offer an algorithm for the numerical solution of this problem.     

基于形状优化框架下的Steklov特征值问题研究

<正>1引言特征值问题在应用数学分支和工程中,尤其是在最优设计问题中,有很多的应用,所以特征值问题的最优化已经有了较为深入的研究,见在我们的研究当中,最优设计问题常常以一种指定载荷的设计下、能量的极小化问题的形式出现.在大多数关于最优设计的文章里面,我们更重视在一个固定载荷下条件下结构的最     

Vortex control in channel flows using translational invariant cost functionals

The use of translation invariant cost functionals for the reduction of vortices in the context of shape optimization of fluid flow domain is investigated. Analytical expressions for the shape design sensitivity involving different cost functionals are derived. Channel flow problems with a bump and an obstacle as possible control boundaries are taken as test examples. Numerical results are provided in various graphical forms for relatively low Reynolds numbers. Striking differences are found for the optimal shapes corresponding to the different cost functionals, which constitute different quantification of a vortex.     

On optimization of the RBF shape parameter in a grid-free local scheme for convection dominated problems over non-uniform centers

Global optimization techniques exist in the literature for finding the optimal shape parameter of the infinitely smooth radial basis functions (RBF) if they are used to solve partial differential equations. However these global collocation methods, applied directly, suffer from severe ill-conditioning when the number of centers is large. To circumvent this, we have used a local optimization algorithm, in the optimization of the RBF shape parameter which is then used to develop a grid-free local (LRBF) scheme for solving convection–diffusion equations. The developed algorithm is based on the re-construction of the forcing term of the governing partial differential equation over the centers in a local support domain. The variable (optimal) shape parameter in this process is obtained by minimizing the local Cost function at each center (node) of the computational domain. It has been observed that for convection dominated problems, the local optimization scheme over uniform centers has produced oscillatory solutions, therefore, in this work the local optimization algorithm has been experimented over Chebyshev and non-uniform distribution of the centers. The numerical experiments presented in this work have shown that the LRBF scheme with the local optimization produced accurate and stable solutions over the non-uniform points even for convection dominant convection–diffusion equations.     

Modelling of topological derivatives for contact problems

The problem of topology optimization is considered for free boundary problems of thin obstacle types. The formulae for the first term of asymptotics for energy functionals are derived. The precision of obtained terms is verified numerically. The topological differentiability of solutions to variational inequalities is established. In particular, the so-called outer asymptotic expansion for solutions of contact problems in elasticity with respect to singular perturbation of geometrical domain depending on small parameter are derived by an application of nonsmooth analysis. Such results lead to the topological derivatives of shape functionals for contact problems. The topological derivatives are used in numerical methods of simultaneous shape and topology optimization. Partially supported by the grant 4 T11A 01524 of the State Committee for the Scientific Research of the Republic of Poland     

Necessary conditions of optimality in constrained planar domain optimization

The article considers optimization of the shape of a planar domain with an integral functional subject to integral inequality constraints and constraints on domain location. The boundary-value problem is a Dirichlet problem for a quasilinear elliptical equation. A method is proposed for deriving necessary conditions of optimality for problems of this kind.Translated from Nelineinye Dinamicheskie Sistemy: Kachestvennyi Analiz i Upravlenie — Sbornik Trudov, No. 3, pp. 22–29, 1993.     

Solving a free boundary problem with nonconstant coefficients

The present article is concerned with the numerical solution of a free boundary problem for an elliptic state equation with nonconstant coefficients. We maximize the Dirichlet energy functional over all domains of fixed volume. The domain under consideration is represented by a level set function, which is driven by the objective's shape gradient. The state is computed by the finite element method where the underlying triangulation is constructed by means of a marching cubes algorithm. We show that the combination of these tools lead to an efficient solver for general shape optimization problems.     

Finding the optimum domain of a nonlinear wave optimal control system by measures

We will explain a new method for obtaining the nearly optimal domain for optimal shape design problems associated with the solution of a nonlinear wave equation. Taking into account the boundary and terminal conditions of the system, a new approach is applied to determine the optimal domain and its related optimal control function with respect to the integral performance criteria, by use of positive Radon measures. The approach, say shape-measure, consists of two steps; first for a fixed domain, the optimal control will be identified by the use of measures. This function and the optimal value of the objective function depend on the geometrical variables of the domain. In the second step, based on the results of the previous one and by applying some convenient optimization techniques, the optimal domain and its related optimal control function will be identified at the same time. The existence of the optimal solution is considered and a numerical example is also given.     

Topology optimization methods for guided flow – comparison of optimality criteria vs. phase-field method

Optimization of guided flow problems is an important task for industrial applications especially those with high Reynolds numbers. There exist several optimization methods to increase the energy efficiency of these problems. Different optimization methods are shown bei Klimetzek [1], Hinterberger [2] and Pingen [3]. In recent years the phase-field method has been shown to be an applicable method for different kinds of topology optimization [4, 5]. We present results of topology optimization methods with optimality criterion and by using a phase-field model in the area of guided fluid flow problems. The two methods aim on the same main target reducing the pressure drop between the inlet and outlet of the flow domain. The first method is based on local optimality criterion, preventing the backflow in the flow domain [1, 6, 7]. The second method is based on a phase field model, which describes a minimization problem and uses a specially constructed driving force to minimize the total energy of the system [4, 5]. We investigate the capabilities and limits of both methods and present examples of different resulting geometries. The initial configurations are prepared in a way that the same optimization problem is solved with both methods. We discuss these results regarding the shape of the improved flow geometry. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)