The limit (as $s\downarrow1/2$) will not change if the double sum is replaced by the corresponding double integral, for which we find this, as desired: [![enter image description here][1]][1] --- Alternatively, the limit (as $s\downarrow1/2$) will not change if the double sum is replaced by the double integral $$I(s):=(2s-1) \iint\limits_{x,y>0,\,x+y>1}\frac{dx\,dy}{x^s y^s(x+y)} \\ =(2s-1)\int_1^\infty \frac{du}u\,\int_0^u \frac{dx}{x^s (u-x)^s} \\ =(2s-1)\int_1^\infty \frac{du}{u^{2s}}\,J(s) =J(s),$$ where $$J(s):=\int_0^1 \frac{dt}{t^s (1-t)^s}\to J(1/2)=\pi.$$ So, $I(s)\to\pi$, as desired. **Details on why the limit (as $s\downarrow1/2$) will not change if the double sum is replaced by an appropriate double integral:** The reasons for this are as follows: 1. $f(x,y):=\dfrac1{x^s y^s(x+y)}$ is decreasing in $x>0$ for each real $y>0$ and in $y>0$ for each real $x>0$. So, for any natural $i$ and $j$, $f(i,j)\ge\int_{(i,i+1]\times(j,j+1]}f\ge f(i+1,j+1)$. 2. For each natural $m$, $\sum_{n\ge1}f(m,n)<\infty$ and, similarly, for each natural $n$, $\sum_{m\ge1}f(m,n)<\infty$. Each of such ordinary sums will get multiplied by $2s-1$, and so, will not contribute to the limit as $s\downarrow1/2$. 2a. Similarly to point 2, for each real $x>0$, $\int_{\max(0,1-x)}^\infty f(x,y)\,dy\le\dfrac1{x^s}\int_0^\infty \min\Big(\dfrac1{y^s},\dfrac1{y^{s+1}}\Big)\,dy =\dfrac1{(1-s)sx^s}$, and $\int_{\max(0,1-y)}^\infty f(x,y)\,dx\le\dfrac1{(1-s)sy^s}$ for each real $y>0$. [1]: https://i.sstatic.net/l8QbX.png