Approximate energy conservation of the beam processes. Results are shown for the B₀ = 35 G and B₀ = 50 G simulations assuming an arcade depth (third-dimension depth) of 100 Mm. The solid line shows the Ohmic heating rate Q_e (integrated over the whole domain) that is subtracted from the MHD model as per Eq. 8 in units of erg s⁻¹. The dashed line shows the beam heating rate (re-interpolated onto the adaptive mesh grid and integrated over the whole domain) that is the source term H_e in the MHD energy equation in units of erg s⁻¹. The dotted line shows, at each time, the total energy stored in the field lines that is carried over from the previous time step, in units of erg.