Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.