A new numerical method based on fictitious domain methods for shapeoptimization problems governed by the Poisson equation is proposed.The basic idea is to combine the boundary variation technique, in whichthe mesh is moving during the optimization, and efficient fictitiousdomain preconditioning in the solution of the (adjoint) state equations.Neumann boundary value problems are solved using an algebraic fictitiousdomain method. A mixed formulation based on boundary Lagrangemultipliers is used for Dirichlet boundary problems and the resultingsaddle-point problems are preconditioned with block diagonal fictitiousdomain preconditioners. Under given assumptions on the meshes, thesepreconditioners are shown to be optimal with respect to the conditionnumber. The numerical experiments demonstrate the efficiency ofthe proposed approaches.
