Let $$f(n) = 2 \sum\limits_{a = 2}^{n - 1} \sum\limits_{b = n + 1}^{n + a - 1}\frac{1}{ab} .$$
One can see that $$\lim\limits_{n \to \infty}f(n) = 2\int\limits_0^1 \frac{dx}{x} \int\limits_1^{1 + x}\frac{dy}{y} = \frac{\pi^2}{6}.$$ This suggests that the double sum defining $f(n)$ can be transformed into a single sum resembling the inverse square sum (although I didn't quite manage to do it after a few minutes with pencil and paper).
Question: Is there a simple closed form for $f(n)$ so that the limit value (or existence, for that matter) is evident?
Motivation: consider a random "wannabe-graceful-tree" graph $G = (V, E)$ with $V = [n]$, and $E$ containing a single random pair $\{x, x + d\}$ for each $d$ from $1 \to n - 1$. Then $f(n)$ is (almost) the expectation of the number of triangles in $G$.