Temporal evolution of the variations in total kinetic, magnetic, and internal energies (ΔE_cin, ΔE_mag, ΔE_int), average current density squared ⟨ j² ⟩, and average linear momentum ⟨ ρv ⟩ in the standard case (in dimensionless units). The grey lines in the top panel represent the same quantities for a case without thermal conduction for comparison. The total kinetic energy is always lower than the magnetic and internal energies. The thermal conductive flux then begins to grow fast as the plasma quickly heats up locally (as a result of the ohmic dissipation) and is responsible for the decay in internal energy during the relaxation phase (note that in our setup the conductive flux can transport heat outwards through the loop footpoints).