Numerical studies of the optimization of the first eigenvalue of the heat diffusion in inhomogeneous media

Abstract

In this paper, we study optimization of the first eigenvalue of $$-abla \cdot (\rho (x) abla u) = \lambda u$$-∇·(ρ(x)∇u)=λu in a bounded domain $$\Omega \subset {\mathbb {R}}^n$$Ω⊂Rn under several constraints for the function $$\rho $$ρ. We consider this problem in various boundary conditions and various topologies of domains. As a result, we numerically observe several common criteria for $$\rho $$ρ for optimizing eigenvalues in terms of corresponding eigenfunctions, which are independent of topology of domains and boundary conditions. Geometric characterizations of optimizers are also numerically observed.