Karsten Eppler, Helmut Harbrecht:

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 domain, especially L^2-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.

Erstveröffentlichung: Preprint WIAS-Preprint No. 1029 (WIAS Berlin)

Mathematics Subject Classification (1991): 49Q10 49M15 65N38 65K10 49K20 65T60

Bibliographical note: Numer. Math. 104, No. 1, 47-68 (2006)

Keywords: shape calculus, Newton method, boundary integral equations, finite element method, multiscale methods, sufficient second order conditions

Mail an Jens Burmeister

[Thu Feb 19 18:56:36 2009]

Impressum