PDE constrained shape optimisation with first-order and Newton-type methods in the $W^{1,\infty}$ topology