We study the evolution of magnetic fields in freely decaying magnetohydrodynamic turbulence. By quasi-linearizing the Navier-Stokes equation, we solve analytically the induction equation in quasinormal approximation. We find that, if the magnetic field is not helical, the magnetic energy and correlation length evolve in time respectively as EB ∝ t −2(1+p)/(3+p) and ξB ∝ t 2/(3+p) , where p is the index of initial power-law spectrum. In the helical case, the magnetic helicity is an almost conserved quantity and forces the magnetic energy and correlation length to scale as EB ∝ (log t) 1/3 t −2/3 and ξB ∝ (log t) −1/3 t 2/3 .PACS numbers: 52.30.Cv, 98.62.EnThe origin of presently-observed large scale magnetic fields throughout the universe is still unclear [1]. Essentially, there are two possible classes of mechanisms to produce cosmic fields depending on when they are generated: Astrophysical mechanisms acting during or after large-scale structure formation, and mechanisms acting in the primordial universe. Magnetic fields created in the early universe (except those generated during inflation), usually suffer from a "small-scale problem", that is their comoving correlation length is much smaller then the characteristic scale of the observed cosmic fields. However, if magnetohydrodynamic (MHD) turbulence operates during their evolution, an enhancement of correlation length can occur, especially if the magnetic field is helical. As pointed out by Banerjee and Jedamzik [2], the evolution of a magnetic field in the early universe goes through different phases depending on the particular conditions of the primordial plasma. In this paper, we are interested in the case of magnetic fields evolving in the turbulent primordial universe well before recombination epoch and when kinematic dissipative effects are due to diffusing particles. Therefore, we are concerned with the so-called phase of "turbulent MHD". In other phases, such as "viscous MHD" and "MHD with ambipolar diffusion" described in Ref.[2], the dynamics of the magnetic field is very different from that studied here. The problem of determining the evolution properties of magnetic fields in MHD turbulence has been deeply and widely discussed in the literature using different methods and approximations. A direct integration of the full set of MHD equations would allow us to deeply understand the dynamics of freely decaying MHD turbulence. However, MHD equations are quite difficult to handle due to their high non-linearity and it has not been yet brought in a definitive verdict for the evolution laws of magnetic energy and correlation length (for recent numerical studies of freely decaying magnetohydrodynamic turbulence see, e.g., Ref. [2,3,4]).The turbulent MHD equations for incompressible fluids, in the case of non-expanding universe, are [5]:and ∇·v = ∇·B = 0. Here, v is the velocity of bulk fluid motion, B the magnetic field, J = ∇ × B the magnetic current, ν the kinematic viscosity, η the resistivity. The thermal pressure of the fluid, p, is not an independen...