Sum of $(1/n^2)\sum_{k=2}^{n/2}\frac{1}{1-\cos\left(\frac{2\pi(k-1)}{n}\right)}$ in the large $n$ limit. This summation appears in studying circulant matrices. Consider an $n-$dimensional circulant matrix with with entries $c_1=c_{n-1}=1$ and $c_0=-2$ (all other coefficients beign $0$). This represents, for example, a random walk on a closed ring of $n$ states. The following calculation can be viewed as the average decay time for all modes (notice that the denominator of the series' terms correspond to the negative eigenvalues of the circulant matrix, $-\lambda_{k}$. Prefactor $1/n^2$ is a normalization constant). From the symmetry of the eigenvalues set one can distinguish two cases: $n$ odd or even.
If $n$ is odd,
$$\frac{1}{n^2}\sum_{k=2}^{(n+1)/2}\frac{1}{1-\cos\left(\frac{2\pi(k-1)}{n}\right)}\,.$$
If $n$ is even,
$$\frac{1}{n^2}\sum_{k=2}^{n/2}\frac{1}{1-\cos\left(\frac{2\pi(k-1)}{n}\right)}\,.$$
I've attempted to approximate this summation by using the continuous limit but I get the wrong values when comparing to numerical evaluation.
The goal is to find an expression for the large $n$ behavior.

Let's denote
$$S(n)=\frac{1}{n^2}\sum_{k=2}^{(n+1)/2}\frac{1}{1-\cos\left(\frac{2\pi(k-1)}{n}\right)}=\frac{1}{2n^2}\sum_{k=1}^{(n-1)/2}\frac{1}{\sin^2\left(\frac{\pi k}{n}\right)}$$
It is not difficult to check that $\displaystyle \phi\geqslant\sin\phi\geqslant\phi-\frac{\phi^3}{6}$ for any $\phi\in \big[0;\frac{\pi}{2}\big]$. Therefore, we have
$$\frac{1}{2\pi^2}\sum_{k=1}^{(n-1)/2}\frac{1}{k^2}\frac{1}{\Big(1-\frac{1}{6}\big(\frac{\pi k}{n}\big)^2\Big)^2}\geqslant S(n)\geqslant\frac{1}{2\pi^2}\sum_{k=1}^{(n-1)/2}\frac{1}{k^2}$$
and, opening the brackets in the first term,
$$\frac{1}{2\pi^2}\sum_{k=1}^{(n-1)/2}\frac{1}{k^2}\frac{1}{1-\frac{1}{3}\big(\frac{\pi k}{n}\big)^2}=\frac{1}{2\pi^2}\sum_{k=1}^{(n-1)/2}\frac{1}{k^2}\bigg(1+\frac{\frac{1}{3}\big(\frac{\pi k}{n}\big)^2}{1-\frac{1}{3}\big(\frac{\pi k}{n}\big)^2}\bigg)\geqslant S(n)\geqslant\frac{1}{2\pi^2}\sum_{k=1}^{(n-1)/2}\frac{1}{k^2}$$
$$\frac{1}{2\pi^2}\bigg(\sum_{k=1}^\infty\frac{1}{k^2}-\sum_{k=(n+1)/2}^\infty\frac{1}{k^2}+\sum_{k=1}^{(n-1)/2}\frac{\frac{1}{3n^2}\big(\frac{\pi }{2}\big)^2}{1-\frac{1}{3}\big(\frac{\pi }{2}\big)^2}\bigg)\geqslant S(n)\geqslant\frac{1}{2\pi^2}\bigg(\sum_{k=1}^\infty\frac{1}{k^2}-\sum_{k=(n+1)/2}^\infty\frac{1}{k^2}\bigg)$$
$$\frac{1}{2\pi^2}\zeta(2)+O\Big(\frac{1}{n}\Big)\geqslant S(n)\geqslant\frac{1}{2\pi^2}\zeta(2)+O\Big(\frac{1}{n}\Big)$$
$$S(n)=\frac{1}{12}+O\Big(\frac{1}{n}\Big)$$
