I sometimes had similar problems, which I solved as follows:
\documentclass{article}
\usepackage{youngtab,young}
\usepackage{amsmath,cancel}
\newcommand{\CenterObject}[1]{\ensuremath{\vcenter{\hbox{#1}}}}
\begin{document}
\section*{Multiplying Young tableaux}
\begin{enumerate}\renewcommand{\labelenumi}{step \arabic{enumi}.}
\item In the first tableau, label all boxes of the first row with an $a$, the
boxes of the second row with a $b$ etc.\label{EnumYoungStep1}
\item
\begin{enumerate}\renewcommand{\labelenumii}{(\alph{enumii})}
\item Sum all schemes with decreasing skyline which may be obtained by
combining the second tableau with boxes of type $a$. Make sure that no
column contains more than $N$ boxes and no two $a$s appear in the same
column.\label{EnumYoungStep2}
\item Continue in the same way with boxes of type $b$.
\label{EnumYoungStep3}
\item Etc.
\end{enumerate}
\item Drop all columns with $N$ boxes (as long as the scheme is not just such
a column).\label{EnumYoungStep4}
\item For each of the resulting schemes, build a string of characters by
reading the first row from the right to the left, then the second row from
the right to the left, and so on. If a given string contains left of an
arbitrary character more $b$s than $a$s or more $c$s than $b$s etc., drop
this string.\label{EnumYoungStep5}
\end{enumerate}\renewcommand{\labelenumi}{\arabic{enumi}.}
\paragraph{Example.}
Consider $\text{SU}(3)$. The gauge bosons transform in the adjoint representation. We
reduce the tensor product of the adjoint representation with itself:
\begin{eqnarray*}
\lefteqn{
\CenterObject{\yng(2,1)}\otimes \CenterObject{\yng(2,1)}
~ \xrightarrow{\mathrm{step}\:\ref{EnumYoungStep1}} ~
\CenterObject{\young(aa,b)} \otimes \CenterObject{\yng(2,1)}} \\
& \xrightarrow{\mathrm{step}\:\mathrm{\ref{EnumYoungStep2}}} &
\CenterObject{\young(\hfil\hfil aa,\hfil)} \oplus \CenterObject{\young(\hfil\hfil a,\hfil a)}
\oplus \CenterObject{\young(\hfil\hfil a,\hfil,a)} \oplus
\CenterObject{\young(\hfil\hfil,\hfil a,a)}\\
& \xrightarrow{\mathrm{step}\:\mathrm{\ref{EnumYoungStep3}}} &
\CenterObject{
\young(\hfil\hfil aab,\hfil)}
\oplus
\CenterObject{\young(\hfil\hfil aa,\hfil b)}
\oplus
\CenterObject{\young(\hfil\hfil aa,\hfil,b)}
\oplus
\CenterObject{\young(\hfil\hfil ab,\hfil a)}
\oplus
\CenterObject{\young(\hfil\hfil a,\hfil ab)}\\
&& {} \oplus
\CenterObject{\young(\hfil\hfil a,\hfil a,b)}
\oplus
\CenterObject{\young(\hfil\hfil ab,\hfil,a)}
\oplus
\CenterObject{\young(\hfil\hfil a,\hfil b,a)}
\oplus
\CenterObject{\young(\hfil\hfil b,\hfil a,a)}
\oplus
\CenterObject{\young(\hfil\hfil,\hfil a,ab)}
\\
& \xrightarrow[\mathrm{step}\:\ref{EnumYoungStep5}]{\mathrm{step}\:\ref{EnumYoungStep4}} &
\cancel{\CenterObject{
\young(\hfil\hfil aab,\hfil)}}
\oplus
\CenterObject{\young(\hfil\hfil aa,\hfil b)}
\oplus \dots
\\
& = &
\CenterObject{
\young(\hfil\hfil\hfil\hfil,\hfil\hfil)
}
\oplus
\CenterObject{
\young(\hfil\hfil\hfil)
}
\oplus
\CenterObject{
\young(\hfil\hfil\hfil,\hfil\hfil\hfil)}
\oplus 2\cdot
\CenterObject{
\young(\hfil\hfil,\hfil)}
\oplus
\CenterObject{
\young(\hfil,\hfil,\hfil)}
\\
& = & \boldsymbol{27} \oplus \boldsymbol{10} \oplus \overline{\boldsymbol{10}}
\oplus 2\cdot \boldsymbol{8} \oplus \boldsymbol{1}\;.
\end{eqnarray*}
\end{document}

UPDATE: Here is an application to your code:
\documentclass[11pt]{scrartcl}
\usepackage{amsmath}
\usepackage{IEEEtrantools}
\usepackage{commath}
\usepackage{lipsum}
\begin{document}
\section{Docendo discimus}
\label{sec:docendo-discimus}
\lipsum[2]
\begin{IEEEeqnarray*}{rCl}
F_{u}(u) &=& \int_{0}^{u}\int_{0}^{y_{1}}3y_{1} \dif y_{2} \dif y_{1} + \int_{u}^{1}\int_{y_{1}-u}^{y_{1}}3y_{1}\dif y_{2} \dif y_{1} \\[0.5em]
&\stackrel{(\ref{step1})}{=}& \int_{0}^{u}\left[\eval{3y_{1}y_{2}}_{0}^{y_{1}}\right] \dif y_{1} + \int_{u}^{1}\left[\eval{3y_{1}y_{2}}_{y_{1}-u}^{y_{1}}\right] \dif y_{1} \\[0.5em]
&\stackrel{(\ref{step2})}{=}& \int_{0}^{u}3y_{1}^{2}\dif y_{1} + \int_{u}^{1}3y_{1}u \dif y_{1}u \\[0.5em]
&\stackrel{(\ref{step3})}{=}& \left[\eval{3 \frac{1}{3}y^{3}}_{0}^{u}\right] + \left[\eval{3 \frac{1}{2}y_{1}^{2}u}_{u}^{1}\right] \\[0.5em]
&\stackrel{(\ref{step4})}{=}& u^{3} + \frac{3}{2}u - \frac{3}{2}u^{3} \\[0.5em]
\IEEEyesnumber
&\stackrel{(\ref{step5})}{=}& \frac{1}{2}(3u - u^{3})
\end{IEEEeqnarray*}
\begin{enumerate}\renewcommand{\labelenumi}{(\arabic{enumi})}
\item\label{step1} In the first step, we perform the $y_2$ integrals.
\item\label{step2} In the second step, we evaluate the inner integrals.
\item\label{step3} In the first step, we perform the $y_1$ integrals.
\item\label{step4} \dots
\item\label{step5} \dots
\end{enumerate}
\lipsum[4]
\end{document}

I'm assuming that eventually your equation numbers will become (section.number)
, otherwise I recommend to label the steps differently.
\marginpar
or\marginnote
. Or use subequations like (1.a) for the lines you want to reference. – Skillmon Feb 12 '18 at 21:28\stackrel{(1)}{=}
,\stackrel{(2)}{=}
and so on, and then refer to these steps. You can even use labels instead of hard-coded numbers. – user121799 Feb 12 '18 at 21:46