I came across this integral when learning the Ising model. Without external field Onsager's solution of a 2D square lattice with $J_2=0$ should equal the solution of a 1D Ising model, which leads to this
$$ \int_0^{2\pi}\ln \left(\cosh 2K-\sinh 2K\cos\theta\right)\,\mathrm{d}\theta=4\pi\ln\cosh K $$ where $K\in\mathbb{R}$. I tested this integral numerically using Mathematica and it holds. But I want to prove it.
I tried splitting the LHS to
$$ 2\pi\ln\cosh2K+\int_0^{2\pi}\ln(1-\tanh 2K\cos\theta)\,\mathrm{d}\theta $$ and expanding the log function using Taylor series, but then the integration result became a complicated series
$$ -\sum_{n=1}^\infty \frac{1}{2n}\frac{2\pi}{4^n}\binom{2n}{n}\tanh^{2n}2K $$ which does not easily relate to the RHS. Anyone has better idea? Ising model is famous so this may have been done long ago.