  The classical level set method, which represents the boundary of the unknown
geometry as the zero-level set of a function, has been shown to be very
effective in solving shape optimization problems. The present work addresses
the issue of using a level set representation when there are simple geometrical
and topological constraints. We propose a logarithmic barrier penalty which
acts to enforce the constraints, leading to an approximate solution to shape
design problems.
