You do very much need "such type of [a] thing" in both cases, as you are really doing the same type of calculation, really.
For the static electron flipflop, you are normally assuming a hamiltonian with an asymmetric term like $S_z$ giving the $|\uparrow\rangle$ state a slightly higher energy eigenvalue than that of the $|\downarrow\rangle$ state. The relevant eigenvalues for the phase differences of the two states are $\epsilon_\pm$, respectively. Thus, for
$$ |\psi(0)\rangle=\cos\theta |\uparrow\rangle+\sin\theta |\downarrow\rangle ~~~~\leadsto \\ |\psi(t)\rangle=e^{-i\epsilon_+ t/\hbar } ~\cos\theta |\uparrow\rangle+e^{-i\epsilon_- t/\hbar } \sin\theta |\downarrow\rangle .
$$
Proceed to compute the oscillating transition amps for $\langle \psi(t)|\psi(0)\rangle $, etc...
For neutrino oscillations, sticking to your two-flavor model,
$$|\Psi(0)\rangle= |\nu_e \rangle, ~~~\hbox {or} ~~ |\nu_{\mu} \rangle ,$$
as you start from the production moment associated with a muon or an electron respectively. But these are not the mass, and hence hamiltonian, eigenstates.
Instead, they are linear combinations of these 1,2 eigenstates,
$$
|\nu_e \rangle = \cos\theta |\nu_1 \rangle+ \sin\theta |\nu_2 \rangle \\
|\nu_{\mu} \rangle = -\sin\theta |\nu_1 \rangle+ \cos\theta |\nu_2 \rangle .
$$
For ultra relativistic propagation of "comoving" mass eigenstates, the simple kinematics yields
$$
E_i\approx E+ {m_i^2\over 2E},
$$
so that, nondimensionalizing $\hbar$ and c as is customary in HEP,
$$
|\nu_e (t) \rangle = e^{-i (E+ m_1^2/2E)t} \cos\theta |\nu_1 \rangle+ e^{-i (E+ m_2^2/2E)t}\sin\theta |\nu_2 \rangle \\
|\nu_{\mu} (t)\rangle = -e^{-i (E+ m_1^2/2E)t}\sin\theta |\nu_1 \rangle+ e^{-i (E+ m_2^2/2E)t}\cos\theta |\nu_2 \rangle ,
$$
from each of which you compute the transition amps $\langle \nu_e (t) |\nu_e (0)\rangle$, etc, as before.
The takeaway is that states which are not eigenstates of the hamiltonian flip flop as they propagate, the common feature of oscillations, here.