On convex domains maximizing the gradient of the torsion function
Linhang Huang
TL;DR
The paper investigates how large the gradient of the torsion function can be on convex planar domains by recasting the PDE problem via conformal maps $f:\mathbb{D}\to\Omega$ and formulating two shape-functionals $\mathbb{L}_1$ and $\mathbb{L}_2$ that are scale-invariant. Through domain-restriction and curvature-signature perturbations, the authors derive two sharp structural conclusions about extremizers: (i) if the boundary around a gradient-maximum point is $C^{1+\alpha}$, no adjacent subarcs can meet at an angle sharper than $\pi/2$, and (ii) extremizers must exhibit either non-$C^{2+\varepsilon}$ regularity or have infinitely many zero-curvature points accumulating near the maximizing boundary point. The core contribution is an Euler-Lagrange framework expressed entirely in conformal geometry, linking the PDE problem to a rich geometric-analytic structure via $v_f(z)=\mathrm{Re}(1+z f''(z)/f'(z))$ and the curvature signature $k_f$. Numerical results based on a Sweers-type family indicate near-optimal domains resemble a rounded 'D', corroborating the theoretical bounds and illustrating the practical impact for shape optimization in elasticity contexts.
Abstract
We consider the solution of $-Δu = 1$ on convex domains $Ω\subset \mathbb{R}^2$ subject to Dirichlet boundary conditions $u =0$ on $\partial Ω$. Our main concern is the behavior of $\|\nabla u\|_{L^{\infty}}$, also known as the maximum shear stress in Elasticity Theory and first investigated by Saint Venant in 1856. We consider the two shape optimization problems $\| \nabla u\|_{L^{\infty}}/ |Ω|^{1/2}$ and $\| \nabla u\|_{L^{\infty}}/ H^1( \partial Ω)$. Numerically, the extremal domain for each functional looks a bit like the rounded letter `D'. We prove that (1) either the extremal domain does not have a $C^{2 + \varepsilon}$ boundary or (2) there exists an infinite set of points on $\partial Ω$ where the curvature vanishes. Either scenario seems curious and is rarely encountered for such problems. The techniques are based on finding a representation of the functional using only conformal geometry and classic perturbation arguments.
