We develop new numerical methods for solving shape optimization problems with several geometrical constraints. Some of these problems are classical problems in spectral theory and have several physical applications, including musical acoustics, civil and mechanical engineering. Previous studies in the literature were restricted to 2D domains and we have been working in the development of robust numerical methods for solving such kind of problems among 3D and 4D domains.