We investigate the dynamics of a closed corona cartesian reduced magnetohydrodynamic (MHD) model where photospheric vortices twist the coronal magnetic field lines. We consider two corotating or counter-rotating vortices localized at the center of the photospheric plate, and additionally more corotating vortices that fill the plate entirely. Our investigation is specifically devoted to study the fully nonlinear stage, after the linear stage during which the vortices create laminar and smoothly twisting flux tubes. Our main goal is to understand the dynamics of photospheric vortices twisting the field lines of a coronal magnetic field permeated by finite amplitude broadband fluctuations. We find that depending on the arrangement and handedness of the photospheric vortices an inverse cascade storing a significant amount of magnetic energy may occur or not. In the first case a reservoir of magnetic energy available to large events such as destabilization of a pre-CME configuration develops, while in the second case the outcome is a turbulent heated corona. Although our geometry is simplified our simulations are shown to have relevant implications for coronal dynamics and CME initiation.

The heating of coronal loops is investigated to understand the observational consequences in terms of the thermodynamics and radiative losses from the Sun as well as the magnetized coronae of stars with an outer convective envelope. The dynamics of the Parker coronal heating model are studied for different ratios of the photospheric forcing velocity timescale $t_p$ to the Alfv\'en crossing time along a loop $t_A$. It is shown that for $t_p/t_A \gtrsim$ 10--24 the heating rate and maximum temperature are largest and approximately independent of $t_p/t_A$, leading to a strong emission in X-rays and EUV. On the opposite decreasing $t_p/t_A$ to smaller values leads to lower heating rates and plasma temperatures, and consequently fading high-energy radiative emission once $t_p/t_A \lesssim$ 1--3. The average volumetric loop heating rate is shown to scale as $\ell_p u_p B_0^2/4\pi L^2$, where $\ell_p$ and $u_p$ are respectively the convective granule length-scale and velocity, $B_0$ is the intensity of the strong magnetic field threading the loop, and $L$ the loop length. These findings support a recent hypothesis explaining ultracool dwarf observations of stars with similar magnetic field strength but radically different topologies displaying different radiative emission.

Ponderomotive acceleration has been asserted to be a cause of the First Ionization Potential (FIP) effect, the by now well known enhancement in abundance by a factor of 3-4 over photospheric values of elements in the solar corona with FIP less than about 10 eV. It is shown here by means of numerical simulations that ponderomotive acceleration occurs in solar coronal loops, with the appropriate magnitude and direction, as a "byproduct" of coronal heating. The numerical simulations are performed with the HYPERION code, which solves the fully compressible three-dimensional magnetohydrodynamic equations including nonlinear thermal conduction and optically thin radiation. Numerical simulations of a coronal loops with an axial magnetic field from 0.005 Teslas to 0.02 Teslas and lengths from 25000 km to 75000 km are presented. In the simulations the footpoints of the axial loop magnetic field are convected by random, large-scale motions. There is a continuous formation and dissipation of field-aligned current sheets which act to heat the loop. As a consequence of coronal magnetic reconnection, small scale, high speed jets form. The familiar vortex quadrupoles form at reconnection sites. Between the magnetic footpoints and the corona the reconnection flow merges with the boundary flow. It is in this region that the ponderomotive acceleration occurs. Mirroring the character of the coronal reconnection, the ponderomotive acceleration is also found to be intermittent.

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is non-uniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales which, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multi-thermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 Tesla are presented. To connect these simulations to observations we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme ultraviolet Imaging Spectrometer (EIS) on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

Context. Photospheric motions shuffle the footpoints of the strong axial magnetic field that threads coronal loops giving rise to turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets where energy is deposited at small-scales and the heating occurs. Previous studies show that current sheets thickness is orders of magnitude smaller than current state of the art observational resolution (~700 km). Aim. In order to understand coronal heating and interpret correctly observations it is crucial to study the thermodynamics of such a system where energy is deposited at unresolved small-scales. Methods. Fully compressible three-dimensional magnetohydrodynamic simulations are carried out to understand the thermodynamics of coronal heating in the magnetically confined solar corona. Results. We show that temperature is highly structured at scales below observational resolution and nonhomogeneously distributed so that only a fraction of the coronal mass and volume gets heated at each time. Conclusions. This is a multi-thermal system where hotter and cooler plasma strands are found one next to the other also at sub-resolution scales and exhibit a temporal dynamics.

The plasma density enhancements recently observed by the Large-Angle Spectrometric Coronagraph (LASCO) instrument onboard the Solar and Heliospheric Observatory (SOHO) spacecraft have sparked considerable interest. In our previous theoretical study of the formation and initial motion of these density enhancements it is found that beyond the helmet cusp of a coronal streamer the magnetized wake configuration is resistively unstable, that a traveling magnetic island develops at the center of the streamer, and that density enhancements occur within the magnetic islands. As the massive magnetic island travels outward, both its speed and width increase. The island passively traces the acceleration of the inner part of the wake. In the present paper a few spherical geometry effects are included, taking into account either the radial divergence of the magnetic field lines and the average expansion suffered by a parcel of plasma propagating outward, using the Expanding Box Model (EBM), and the diamagnetic force due to the overall magnetic field radial gradients, the so-called melon-seed force. It is found that the values of the acceleration and density contrasts can be in good agreement with LASCO observations, provided the spherical divergence of the magnetic lines starts beyond a critical distance from the Sun and the initial stage of the formation and acceleration of the plasmoid is due to the cartesian evolution of MHD instabilities. This result provides a constraint on the topology of the magnetic field in the coronal streamer.

MHD turbulence has long been proposed as a mechanism for the heating of coronal loops in the framework of the Parker scenario for coronal heating. So far most of the studies have focused on its dynamical properties without considering its thermodynamical and radiative features, because of the very demanding computational requirements. In this paper we extend this previous research to the compressible regime, including an energy equation, by using HYPERION, a new parallelized, viscoresistive, three-dimensional compressible MHD code. HYPERION employs a Fourier collocation -- finite difference spatial discretization, and uses a third-order Runge-Kutta temporal discretization. We show that the implementation of a thermal conduction parallel to the DC magnetic field induces a radiative emission concentrated at the boundaries, with properties similar to the chromosphere--transition region--corona system.

The Parker or field line tangling model of coronal heating is studied comprehensively via long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy. In physical space this corresponds to a magnetic topology where magnetic field lines are barely entangled, nevertheless current sheets (corresponding to the original tangential discontinuities hypothesized by Parker) are continuously formed and dissipated. Current sheets are the result of the nonlinear cascade that transfers energy from the scale of convective motions ($\sim 1,000 km$) down to the dissipative scales, where it is finally converted to heat and/or particle acceleration. Current sheets constitute the dissipative structure of the system, and the associated magnetic reconnection gives rise to impulsive ``bursty'' heating events at the small scales. This picture is consistent with the slender loops observed by state-of-the-art (E)UV and X-ray imagers which, although apparently quiescent, shine bright in these wavelengths with little evidence of entangled features. The different regimes of weak and strong MHD turbulence that develop, and their influence on coronal heating scalings, are shown to depend on the loop parameters, and this dependence is quantitatively characterized: weak turbulence regimes and steeper spectra occur in {\it stronger loop fields} and lead to {\it larger heating rates} than in weak field regions.

Long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity $B_0$, which depends on the spectral index of total energy for given loop parameters, must vary from $B_0^{3/2}$ for weak fields to $B_0^{2}$ for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet Sun coronal energy losses.

We simulate the Galactic Center nonthermal filaments as magnetized wakes formed dynamically from amplification of a weak (tens of $\mu$G) global magnetic field through the interaction of molecular clouds with a Galactic Center wind. One of the key issues in this cometary model is the stability of the filament against dynamical disruption. Here we show 2-dimensional MHD simulations for interstellar conditions that are appropriate for the Galactic Center. The structures eventually disrupt through a shear driven nonlinear instability but maintain coherence for lengths up to 100 times their width as observed. The final instability, which destroys the filament through shredding and plasmoid formation, grows quickly in space (and time) and leads to an abrupt end to the structure, in accord with observations. As a by-product, the simulation shows that emission should peak well downstream from the cloud-wind interaction site.