Recrystallization and its related phenomena play an important role in the evolution of the microstructure-property relationship of ferritic steels during annealing after cold-rolling. The kinetics of recrystallization phenomena is heavily influenced by the mobility of high-angle grain boundaries (HAGBs). Here, alloying elements, such as Nb, V, and Ti show strong segregation tendencies to HAGBs and impede their movement, which is known as solute drag. A better understanding of this effect is required for innovative steel production processes for $\mathrm{CO_2}$ reduction like the electric arc furnace (EAF) route, where impurity elements enter from steel scrap. \newline

In this work, we employ a state-parameter based mean-field model, which describes static recrystallization in terms of nucleation and growth of recrystallized grains, and incorporates the solute drag effect as a function of alloying elements in work-hardened ferritic steels. The solute drag effect from different elements is usually used as a fitting parameter in recrystallization models. We go beyond this approach by computing solute binding energies of characteristic $\mathrm{\Sigma5}$ HAGB for a wide range of elements by \emph{ab-initio} simulations and consider their temperature and composition dependent influence on HAGB mobility via the Cahn approach. The prevalent dislocation density is obtained from the ABC-dislocation density model, which is parameterised on flow curves of relevant steels. With this model framework at hand, we can show good agreement to experimental data.