Rate of convergence of the $p$-torsion function to the distance function

We prove that the Dirichlet $p$-torsion function $u_p$ converges to the distance to the boundary $d$ at the rate $O(1/p)$ in the uniform norm, on every bounded domain and with a constant depending only on the dimension and the diameter. The upper bound comes from a radial barrier with its pole at a