It seems that the following functions $$G(R,s)=(1-R^2)\int_0^1\int_0^{2\pi} Adr da,$$ where $$A=\frac{ \sqrt{\left(\left(1+r^2\right) \cos(a+s)-2 r R \cos s\right)^2+\left(1-r^2\right)^2 \sin^2(a+s)}}{\left(1+r^2 R^2-2 r R \cos a\right)^2}$$ attains its maximum for $R=0$, if $R\in[0,1]$, however I dont have the proof.
Geometricaly, it is equivalent to find $$\sup_{z\in U}\int_{U}\left|\frac{1}{z-w}+\frac{w}{1-\bar z w}\right|du dv, \ \ w=u+iv,$$ where $U$ is the unit disk.