DG FEM numerical simulation of a plasma jet in a stellar atmosphere

Motivated by studying an outflow of plasma caused by magnetic reconnection in Solar and stellar flares we use a novel, Discontinuous Galerkin Finite Element Method (DG FEM) MHD code to study the behavior of the jet and of subsequent turbulence. We show conditions, under which the turbulence develops and general behavior under various initial parameters: speed of the flow, plasma beta, background density, and magnetic configuration. Furthermore, we focus on the current sheet-like configuration of the magnetic field and super alfvénic flows. Under these conditions, we present power-law fit of the energetic spectrum, energy conservation, and divergence-free nature of the simulation even without excessive artificial resistivity and viscosity.