We present a series of numerical simulations of the quiet Sun plasma threaded by magnetic fields that extend from the upper convection zone into the low corona. We discuss an efficient, simplified approximation to the physics of optically thick radia
tive transport through the surface layers, and investigate the effects of convective turbulence on the magnetic structure of the Suns atmosphere in an initially unipolar (open field) region. We find that the net Poynting flux below the surface is on average directed toward the interior, while in the photosphere and chromosphere the net flow of electromagnetic energy is outward into the solar corona. Overturning convective motions between these layers driven by rapid radiative cooling appears to be the source of energy for the oppositely directed fluxes of electromagnetic energy.
The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faradays law to find an approximate solution for the electric field in the photosphere, using
a Poloidal-Toroidal Decomposition (PTD) of the magnetic field and its partial time derivative. Without additional information, however, the electric field found from this technique is under-determined -- Faradays law provides no information about the electric field that can be derived the gradient of a scalar potential. Here, we show how additional information in the form of line-of-sight Doppler flow measurements, and motions transverse to the line-of-sight determined with ad-hoc methods such as local correlation tracking, can be combined with the PTD solutions to provide much more accurate solutions for the solar electric field, and therefore the Poynting flux of electromagnetic energy in the solar photosphere. Reliable, accurate maps of the Poynting flux are essential for quantitative studies of the buildup of magnetic energy before flares and coronal mass ejections.
We compute the change in the Lorentz force integrated over the outer solar atmosphere implied by observed changes in vector magnetograms that occur during large, eruptive solar flares. This force perturbation should be balanced by an equal and opposi
te force perturbation acting on the solar photosphere and solar interior. The resulting expression for the estimated force change in the solar interior generalizes the earlier expression presented by Hudson, Fisher and Welsch (CS-383, ASP, 221, 2008), providing horizontal as well as vertical force components, and provides a more accurate result for the vertical component of the perturbed force. We show that magnetic eruptions should result in the magnetic field at the photosphere becoming more horizontal, and hence should result in a downward (towards the solar interior) force change acting on the photosphere and solar interior, as recently argued from an analysis of magnetogram data by Wang and Liu (Astrophys. J. Lett. 716, L195, 2010). We suggest the existence of an observational relationship between the force change computed from changes in the vector magnetograms, the outward momentum carried by the ejecta from the flare, and the properties of the helioseismic disturbance driven by the downward force change. We use the impulse driven by the Lorentz-force change in the outer solar atmosphere to derive an upper limit to the mass of erupting plasma that can escape from the Sun. Finally, we compare the expected Lorentz-force change at the photosphere with simple estimates from flare-driven gasdynamic disturbances and from an estimate of the perturbed pressure from radiative backwarming of the photosphere in flaring conditions.
We introduce two new methods that are designed to improve the realism and utility of large, active region-scale 3D MHD models of the solar atmosphere. We apply these methods to RADMHD, a code capable of modeling the Suns upper convection zone, photos
phere, chromosphere, transition region, and corona within a single computational volume. We first present a way to approximate the physics of optically-thick radiative transfer without having to take the computationally expensive step of solving the radiative transfer equation in detail. We then briefly describe a rudimentary assimilative technique that allows a time series of vector magnetograms to be directly incorporated into the MHD system.
Where does solar flare energy come from? More specifically, assuming that the ultimate source of flare energy is mechanical energy in the convection zone, how is this translated into energy dissipated or stored in the corona? This question appears to
have been given relatively little thought, as attention has been focussed predominantly on mechanisms for the rapid dissipation of coronal magnetic energy by way of MHD instabilities and plasma micro instabilities. We consider three types of flare theory: the steady state photospheric dynamo model in which flare power represents coronal dissipation of currents generated simultaneously by sub-photospheric flows; the magnetic energy storage model where sub-photospheric flows again induce coronal currents but which in this case are built up over a longer period before being released suddenly; and emerging flux models, in which new magnetic flux rising to the photosphere already contains free energy, and does not require subsequent stressing by photospheric motions. We conclude that photospheric dynamos can power only very minor flares; that coronal energy storage can in principle meet the requirements of a major flare, although perhaps not the very largest flares, but that difficulties in coupling efficiently to the energy source may limit this mechanism to moderate sized flares; and that emerging magnetic flux tubes, generated in the solar interior, can carry sufficient free energy to power even the largest flares ever observed.
Determining the electric field (E-field) distribution on the Suns photosphere is essential for quantitative studies of how energy flows from the Suns photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input
for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faradays Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a poloidal-toroidal decomposition (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faradays Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faradays Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcopes Minimum Energy Fit. The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulations known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data.
