The origin of cosmic magnetic fields is an unsolved problem and magnetogenesis could have occurred in the early Universe. We study the evolution of such primordial magnetic fields across the cosmological recombination epoch via 3D magnetohydrodynamic numerical simulations. We compute the effective or net heating rate of baryons due to decaying magnetic fields and its dependence on the magnetic field strength and spectral index. In the drag-dominated regime ($z \gtrsim 1500$), prior to recombination, we find no real heating is produced. Our simulations allow us to smoothly trace a new transition regime ($600 \lesssim z \lesssim 1500$), where magnetic energy decays, at first, into the kinetic energy of baryons. A turbulent velocity field is built up until it saturates, as the net heating rate rises from a low value at recombination to its peak towards the end of the transition regime. This is followed by a turbulent decay regime ($z \lesssim 600$) where magnetic energy dissipates via turbulent decay of both magnetic and velocity fields while net heating remains appreciable and declines slowly. Both the peak of the net heating rate and the onset of turbulent decay are delayed significantly beyond recombination, by up to 0.5 Myr (until $z\simeq 600-700$), for scale-invariant magnetic fields. We provide analytic approximations and present numerical results for a range of field strengths and spectral indices, illustrating the redshift-dependence of dissipation and net heating rates. These can be used to study cosmic microwave background constraints on primordial magnetic fields.

We explore the phenomenon of the recently discovered inverse transfer of energy from small to large scales in decaying magnetohydrodynamical turbulence by Brandenburg et al. (2015) even for nonhelical magnetic fields. For this investigation we mainly employ the Pencil-Code performing a parameter study, where we vary the Prandtl number, the kinematic viscosity and the initial spectrum. We find that in order to get a decay which exhibits this inverse transfer, large Reynolds numbers ($\mathcal{O}\sim 10^{3}$) are needed and low Prandtl numbers of the order unity $Pr = 1$ are preferred. Compared to helical MHD turbulence, though, the inverse transfer is much less efficient in transferring magnetic energy to larger scales than the well-known effect of the inverse cascade. Hence, applying the inverse transfer to the magnetic field evolution in the Early Universe, we question whether the nonhelical inverse transfer is effective enough to explain the observed void magnetic fields if a magneto- genesis scenario during the electroweak phase transition is assumed.