Here's a way to think about quantum numbers that halfway explains them, but leaves a different mystery about them.
For a typical one-dimensional problem (and if there is spherical symmetry, we can factor the wavefunction into a radial wave function and an angular wave function--the radial wave function is essentially a one-dimensional problem), we can write Schrodinger's equation this way:
\frac{\partial^2 \psi}{\partial r^2} = \frac{-2m}{\hbar^2} (E - V) \psi
In a typical problem we have two regions:
- Oscillatory region, where E - V \gt 0
- Non-oscillatory region, where E - V \lt 0
When E - V \gt 0, we can define k(r) to be \sqrt{\frac{2m}{\hbar^2} (E-V)}. In that region, the equation looks like:
\frac{\partial^2 \psi}{\partial r^2} = -k^2(r) \psi
That equation has oscillatory solutions, like sines and cosines (it would be exactly sines and cosines if k were constant). An oscillatory function goes up and down, and can be characterized by a nonnegative integer, namely, the number of zeros--the number of points where \psi(r) = 0. That integer is the radial quantum number, n_r.
That's an explanation for why the radial quantum number must be an integer. It doesn't really explain why energy levels are quantized, though. If you tried to numerically solve Schrodinger's equation for an arbitrary energy E, you would definitely get an integer number of zeros, no matter what the energy was. But what's special about energy eigenvalues can only be understood by looking at the other region--the non-oscillatory region.
In this region, we have E - V \lt 0. We can define \kappa(r) to be \sqrt{\frac{2m}{\hbar^2} (V-E)} In that region, the equation looks like:
\frac{\partial^2 \psi}{\partial r^2} = +\kappa^2(r) \psi
Rather than oscillatory solutions, this equation has solutions that are like exponentials, they either blow up as r \rightarrow \infty, or they go to zero. If \kappa were constant, then these two solutions would be e^{+\kappa r} and e^{-\kappa r}. This is where the quantization of energy levels comes from. If you try to solve the Schrodinger equation for an arbitrary energy level, what you'll find is that if you make sure that it is finite at r=0, then it will blow up as r \rightarrow \infty. If you try to make it go to zero as r \rightarrow \infty, then the wave function will blow up at r=0. Only for very special energies will it be possible to have the wave function be well behaved both at r=0 and at r \rightarrow \infty.
What's not obvious from my explanation, and I think that the mathematics of proving this is beyond me, is the connection between those two different integral numbers. Why is there exactly one energy E_0 that has one zero and has nice behavior as r \rightarrow \infty and r \rightarrow 0?