Boundary element based multiresolution shape optimisation in electrostatics

We consider the shape optimisation of high-voltage devices subject to electrostatic field equations by combining fast boundary elements with multiresolution subdivision surfaces. The geometry of the domain is described with subdivision surfaces and different resolutions of the same geometry are used for optimisation and analysis. The primal and adjoint problems are discretised with the boundary element method using a sufficiently fine control mesh. For shape optimisation the geometry is updated starting from the coarsest control mesh with increasingly finer control meshes. The multiresolution approach effectively prevents the appearance of non-physical geometry oscillations in the optimised shapes. Moreover, there is no need for mesh regeneration or smoothing during the optimisation due to the absence of a volume mesh. We present several numerical experiments and one industrial application to demonstrate the robustness and versatility of the developed approach.We gratefully acknowledge the support provided by the EU commission through the FP7 Marie Curie IAPP project CASOPT (PIAP-GA-2008-230224). K.B. and F.C. thank for the additional support provided by EPSRC through #EP/G008531/1. J.Z. thanks for the support provided by the European Regional Development Fund in the IT4Innovations Centre of Excellence project (CZ.1.05/1.1.00/02.0070) and by the project SPOMECH – Creating a Multidisciplinary R&D Team for Reliable Solution of Mechanical Problems, reg. no. CZ.1.07/2.3.00/20.0070 within the Operational Programme ‘Education for Competitiveness’ funded by the Structural Funds of the European Union and the state budget of the Czech Republic. Special thanks to Andreas Blaszczyk from the ABB Corporate Research Center Switzerland for fruitful discussions and for providing the industrial applications.This is the final version of the article. It first appeared from Elsevier via http://dx.doi.org/10.1016/j.jcp.2015.05.01