Is there a possibility to analytically determine how often a Lissajous curve on the surface of a sphere crosses itself?
The Lissajous figure (on the unit sphere) is given by:
x = $ \sin(n_{\theta} t) \cos(n_{\phi} t) $
y = $ \sin(n_{\theta} t) \sin(n_{\phi} t) $
z = $ \cos(n_{\theta} t) $
where $n_{\theta}$ and $n_{\phi}$ are arbitrary integers and $t \in [0, 2\pi)$.
An example for such a trajectory on a sphere with $n_{\phi}=21$ and $n_{\theta}=20$ is given in this image: 
The poles are obviously two crossing points. If you want to consider each pair of arcs separately, there are obviously $n_\theta(n_\theta+1)$ polar crossings.
Having made some sketches, I think there are simply $n_\theta n_\phi$ non-polar crossings, for a total of $n_\theta (n_\phi + n_\theta + 1)$.