The work presented here aims at better understanding the complex non linear dynamics of a large magnetization subject to efficient radiation damping. To do so, the dynamical effects of the dipolar field were investigated through numerical simulations of a simplified model, which provided results that are consistent with several observed effects, both in this and previous work. In particular, it was found that for strong dipolar fields, the predicted evolution of the magnetization (Figures 7a) was similar to the observed one and depicted in Figures , that is, a "steady state" of regularly spaced maser bursts. It is also noteworthy that in the absence of dipolar field, the model gives qualitatively the same predictions as the {3+1}-dimensional simplified BMP model introduced in Reference . This transition is illustrated in Figure . It is also interesting that we were not able to reproduce the kind of long term evolution either with the simplified BMP model, or with the present lattice model without DDF. In our view, these elements strongly suggest that the kind of induction signal observed in this work attests for the presence of an effective DDF. The noted asymmetry and the "bursts within burst" profile are clearly also due to the presence of the DDF. However, the mechanism through which the dipolar field produces these effects is somewhat unclear. The dipolar field essentially contributes an average z component that is much larger than its transverse x and y components. The average ⟨ω d,z ⟩ = -⟨B d,z (r)⟩/γ component is depicted in Figure for a trajectory of model (B). It appears that its time evolution mirrors the M z component of the magnetization, and therefore contributes a time-varying precession frequency to the magnetization. Considering its effect on an isolated burst, one would expect its effect to be that of a chirp pulse, with a time-varying offset frequency. However, the local dipolar field B d (r i ) at each site of the sample is not uniform, achieving a distribution of its z components, therefore of precession frequencies across the sample, as shown in Figure . From the dynamical viewpoint, this amounts to the existence of an inhomogeneity of B 0 , albeit a time-dependent one. This is interesting, as it points to an effect that was already invoked in the context of solution NMR, where B 0 inhomogeneity was shown to lead to such non-decaying sustained maser pulses, an effect not explained by the Bloch-Maxwell equations.