Top: convergence study of dissipation with grid spacing and numerical resolution using 2D simulations (see Table 1 for grid parameters). The 3D run with χ = 85° is reported here for comparison with 2D runs (black solid line). Middle: dependence of dissipation with the plasma magnetization for κ_⋆ = 20. Bottom: dependence of dissipation with the plasma multiplicity for σ_⋆ = 250. The dashed lines are best-fit models to the analytical solution proposed in Sect. 4, Eq. (20). In all panels, the Poynting flux is normalized to its light-cylinder value, L/L_LC.