On convex domains maximizing the gradient of the torsion function

We consider the solution of $-\Delta u = 1$ on convex domains $\Omega \subset \mathbb{R}^2$ subject to Dirichlet boundary conditions $u =0$ on $\partial \Omega$. 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}}/ |\Omega|^{1/2}$ and $\| \nabla u\|_{L^{\infty}}/ H^1( \partial \Omega)$. 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 \Omega$ 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.