Following Gerald Edgar, $$W_3(1)\equiv\mathbb{E}[|Y|]=\int dx\int dy \,\sqrt{x^2+y^2}\,p(x,y)=2\pi\int_0^\infty r^2p(r)dr$$ with $p(x,y)dxdy=p(r)\,rdrd\phi$ the rotationally invariant distribution of the complex variable $Y=x+iy=r\cos\phi+ir\sin\phi$. From here $$\mathbb{E}[|X|]=\int dx\int dy \,|x|\,p(x,y)=\int_0^\infty \int_0^{2\pi} |\cos\phi|r^2 p(r)\,drd\phi=\frac{2}{\pi}W_3(1).$$ We thus arrive at $$\int\limits_0^{2\pi}\int\limits_0^{2\pi}\int\limits_0^{2\pi}|\cos x+\cos y+\cos z|\ dx\ dy\ dz=(2\pi)^3\mathbb{E}[|X|]=16\pi^2 W_3(1)$$ $$\qquad=3\;\frac{2^{1/3}}{\pi^2}\Gamma({\textstyle \frac{1}{3}})^6 + 108\;\frac{2^{2/3}}{\pi^2}\Gamma({\textstyle \frac{2}{3}})^6 \approx 248.65$$ --- curiosities: - since $(2\pi)^3\approx 248.05$, the average of the absolute value of the sum of three cosines equals unity within 0.2%. - the integral $\int_0^{2\pi}\cdots\int_0^{2\pi}|\cos x_1+\cdots + \cos x_n|\,dx_1\cdots dx_n$ equals a power of 2 for $n=1$ and $n=2$, but not for $n=3$.