Studies of spin dynamics of highly polarized dilute mixtures of 3He in superfluid 4He have been performed by various researchers over the past three decades. One series of experiments performed at Cornell University in the early 1990’s revealed a novel long timescale excitation. We present the numerical solution of the non-linear Leggett spin dynamics equation in one spatial dimension subject to boundary conditions consistent with the Cornell experiments. Experimentally observed phenomena are composed of trains of bursts in the transverse magnetization lasting several seconds. The simulations capture the time evolution of the individual bursts localized in time. Preliminary results of two dimensional simulations are also presented.