I'm going to answer how the Bohr calculation works. The potential is $V(r)=-Z/r$. The force is $F=Z/r^2$. This has to match the centripetal force $mv^2/r$. So $v=\sqrt\frac{Z}{mr}$. The momentum is $p=mv=\sqrt\frac{mZ}{r}$. This has, so far, been classical.
Quantum mechanical duality tells us that a particle of momentum $p$ is a wave of wavelength $\lambda=2\pi\hbar/p$ (the De-Broglie relation). But the wave has to "fit""fit nicely" (i.e. no discontinuities in the wave or its derivatives) into the circular orbit which has length $2\pi r$. So the wavelength has to be quantized: $\lambda_n=2\pi r/n$ (with $n=1,...,\infty$). So $p$ has to be quantized: $p_n=n\hbar/r$.
Equating this with the classical expression for $p$ we obtained earlier, we then have $\sqrt\frac{mZ}{r}=n\hbar/r$ which gives the quantized radii: $r_n=\frac{\hbar^2n^2}{mZ}$. Substituting this radius into the potential gives the quantized energy levels: $V(r_n)=E_n=-\frac{mZ^2}{\hbar^2n^2}$. The actual result is $E_n=-\frac{mZ^2}{2\hbar^2n^2}$ so we're off by a factor of 2.