Initially assuming that the input current is sinusoidal.
and
Why does the presence of the third harmonic in the flux, alter the
current (by adding a third harmonic out of phase by a 180∘)?
That's the bottom line - no matter what the flux does (in terms of non-linearity), the current is sinusoidal by definition.
However, if the winding were driven by a sine wave voltage then that's a different story and, current can appear very oddly affected by the diminishing flux at higher peak levels: -
GIF from here.
Reason
In a normal inductor this formula applies \$V = L\dfrac{di}{dt}\$
This gives us the sine-cosine relationship between voltage and current but, only if inductance is unchanging. As the core saturates, inductance drops dramatically and this, in turn, leads to a much greater \$\frac{di}{dt}\$.
Then, knowing\$^1\$ that the definition of inductance per turn is \$\dfrac{\Phi}{I}\$, if \$L\$ decreases and \$I\$ increases, then flux remains unaffected.
\$^1\$ We are taught that \$V = L\dfrac{di}{dt}\$ and this is taken from Faraday's law: \$V = N\dfrac{d\Phi}{dt}\$.
So, if we equate them: \$\dfrac{L}{N} = \dfrac{d\Phi}{di}\$ or, inductance per turn equals flux per amp.