The simulations presented in this work qualitatively explain a number of features of the magnetization dynamics under DNP and radiation feedback conditions. The models introduced in this work involve a large number of parameters and time scales (γ n,ee , γ rd , γ st , γ 2 , shape of the sample and spin density) that control the kind of dynamics that the magnetization undergoes, and exploring the complete parameter space is a formidable task. We have restricted the study to ranges of parameters that were able to reproduce experimental observations and, as it turns out, some of them may seem remote from the experimental ones. Although these values were arbitrary to some extent, the important point made in our study are the roles played by the relative values of the ratios γ rd /γ 2 and ω dip /γ rd : the former determines the maser threshold, whereas the latter controls the manifestation of DDF effects. Moreover, we have assumed a time scale separation between γ rd , γ 2 and ω dip , on the one hand, and on the much slower rates γ n,ee and γ st , on the other hand. This corresponds to the much longer polarization buildup times.