Shape optimization of harmonic helicity in toroidal domains

In this paper, we introduce a new shape functional defined for toroidal domains that we call harmonic helicity, and study its shape optimization. Given a toroidal domain, we consider its associated harmonic field. The latter is the magnetic field obtained uniquely up to normalization when imposing zero normal trace and zero electrical current inside the domain. We then study the helicity of this field, which is a quantity of interest in magneto-hydrodynamics corresponding to the L2 product of the field with its image by the Biot--Savart operator. To do so, we begin by discussing the appropriate functional framework and an equivalent PDE characterization. We then focus on shape optimization, and we identify the shape gradient of the harmonic helicity. Finally, we study and implement an efficient numerical scheme to compute harmonic helicity and its shape gradient using finite elements exterior calculus.