We discuss large scale numerical simulations as a tool for obtaining qualitative understanding of the processes directly and indirectly responsible for coronal heating. The actual heating process in the low beta coronal plasma is most likely driven by transfer of magnetic energy from the subsurface high beta region, where magnetic energy is created as an energetically insignificant byproduct of solar convection and rotation. Based on the results of recent numerical experiments, we discuss some of the processes involved.