Shape optimization with inradius constraint and emergence of honeycomb structures

We maximize the average torsional rigidity among open subsets of $\mathbb{R}^2$ with fixed inradius whose complements consist of pairwise disjoint balls of fixed radius $\varepsilon$ with mutual distances bounded below. We prove that, provided $\varepsilon$ is sufficiently small, the optimal value i