Yung Mok

Calculating the Thermal Structure of Solar Active Regions in Three Dimensions

We describe a technique to obtain the temperature and density distribution in an active region for a specified plasma heating model. The technique can be applied in general to determine the magnetic field and thermal structure self-consistently. For simplicity, we illustrate the application of this technique in the limit of small plasma β, in which the plasma dynamics decouples from that of the magnetic field, a good approximation in active regions, in which the magnetic field is strong. We select a particular active region, observed in 1996 August, to demonstrate the methodology. We apply the technique to a force-free magnetic field with a plasma heating model in which the volumetric coronal heating rate is directly proportional to the strength of the local magnetic field, and we compute the expected extreme-ultraviolet and soft X-ray emissions from the resulting thermal structure. We compare our solutions with one-dimensional loop models and analytic loop scaling laws. In the future, we plan to compare these emission images with those obtained by the SOHO EUV Imaging Telescope (EIT) and the Yohkoh Soft X-Ray Telescope (SXT) and to explore the relationship between coronal emission and various coronal heating models.

THE FORMATION OF CORONAL LOOPS BY THERMAL INSTABILITY IN THREE DIMENSIONS

Plasma loops in solar active regions have been observed in EUV and soft X-rays for decades. Their formation mechanism and properties, however, are still not fully understood. Predictions by early models, based on 1D hydrostatic equilibria with uniform plasma heating, are not consistent with high-resolution measurements. In this Letter, we demonstrate, via 3D simulations, that a class of heating models can lead to the dynamic formation of plasma loops provided the plasma is heated sufficiently to match SXT soft X-ray measurements. We show that individual flux tubes in a 3D magnetic structure tend to stand out against their neighbors. The loops have large aspect ratios and nearly uniform cross sections in the corona, similar to those observed by EIT and TRACE. The coronal EUV emission from these thermally unstable solutions is roughly consistent with EIT measurements. The solution oscillates in time through a large-amplitude, nonlinear cycle, leading to repeated brightening and fading of the loops.

The Importance of Geometric Effects in Coronal Loop Models

We systematically investigate the effects of geometrical assumptions in one-dimensional (1D) models of coronal loops. Many investigations of coronal loops have been based on restrictive assumptions, including symmetry in the loop shape and heating profile, and a uniform cross-sectional area. Starting with a solution for a symmetric uniform-area loop with uniform heating, we gradually relax these restrictive assumptions to consider the effects of nonuniform area, nonuniform heating, a nonsymmetric loop shape, and nonsymmetric heating, to show that the character of the solutions can change in important ways. We find that loops with nonuniform cross-sectional area are more likely to experience thermal nonequilibrium, and that they produce significantly enhanced coronal emission, compared with their uniform-area counterparts. We identify a process of incomplete condensation in loops experiencing thermal nonequilibrium during which the coronal parts of loops never fully cool to chromospheric temperatures. These solutions are characterized by persistent siphon flows. Their properties agree with observations (Lionello et al.) and may not suffer from the drawbacks that led Klimchuk et al. to conclude that thermal nonequilibrium is not consistent with observations. We show that our 1D results are qualitatively similar to those seen in a three-dimensional model of an active region. Our results suggest that thermal nonequilibrium may play an important role in the behavior of coronal loops, and that its dismissal by Klimchuk et al., whose model suffered from some of the restrictive assumptions we described, may have been premature.

Prominence formation in a coronal loop

A model is presented which depends on the preferential deposition of heating in the legs of a coronal loop and which produces a stable prominence-scale condensation at the loop top. Dynamic stability is attained by the subsequent adjustment of local parallel gravity by a magnetic inversion at the loop (or arcade) apex. A nonlinear numerical simulation of this process, which includes a deep chromosphere, a heating rate with a fixed dissipation length, and full solar gravity is described. 12 refs.

The thermal instability in a sheared magnetic field - Filament condensation with anisotropic heat conduction

The condensation-mode growth rate of the thermal instability in an empirically motivated sheared field is shown to depend upon the existence of perpendicular thermal conduction. This typically very small effect (perpendicular conductivity/parallel conductivity less than about 10 to the -10th for the solar corona) increases the spatial-derivative order of the compressible temperature-perturbation equation, and thereby eliminates the singularities which appear when perpendicular conductivity = 0. The resulting growth rate is less than 1.5 times the controlling constant-density radiation rate, and has a clear maximum at a cross-field length of order 100 times and a width of about 0.1 the magnetic shear scale for solar conditions. The profiles of the observable temperature and density perturbations are independent of the thermal conductivity, and thus agree with those found previously. An analytic solution to the short-wavelength incompressible case is also given.

Coronal loop formation resulting from photospheric convection

We have demonstrated the dynamic formation of coronal magnetic loops in three dimensions as a result of horizontal vortex-like convection on the photosphere. Localized plasma motions twist bipolar magnetic field lines which are tied to the dense photosphere by high electrical conductivity. The twists propagate into the corona along the field and create a narrow quasi-toroidal region where the field lines interwind. At the same time, this tubeline region rises in altitude, expands in cross section, and distorts into a slight S shape before settling into an equilibrium state. The MHD stability of such line-tied magnetic loop structures is directly exhibited by this dynamic simulation.

Thermal Non-equilibrium Revisited: A Heating Model for Coronal Loops

The location and frequency of events that heat the million-degree corona are still a matter of debate. One potential heating scenario is that the energy release is effectively steady and highly localized at the footpoints of coronal structures. Such an energy deposition drives thermal non-equilibrium solutions in the hydrodynamic equations in longer loops. This heating scenario was considered and discarded by Klimchuk et al. on the basis of their one-dimensional simulations as incapable of reproducing observational characteristics of loops. In this paper, we use three-dimensional simulations to generate synthetic emission images, from which we select and analyze six loops. The main differences between our model and that of Klimchuk et al. concern (1) dimensionality, (2) resolution, (3) geometrical properties of the loops, (4) heating function, and (5) radiative function. We find evidence, in this small set of simulated loops, that the evolution of the light curves, the variation of temperature along the loops, the density profile, and the absence of small-scale structures are compatible with the characteristics of observed loops. We conclude that quasi-steady footpoint heating that drives thermal non-equilibrium solutions cannot yet be ruled out as a viable heating scenario for EUV loops.

Interaction of Two Magnetic Loops in the Solar Corona

The solar corona is populated by a large number of semitoroidal magnetic loops, some of which are sufficiently close to each other within a neighborhood that the probability of loop-to-loop interaction is not negligible. The interaction of two coronal loops is studied using a three-dimensional numerical simulation. The first loop is an established, current-carrying magnetic loop in hydromagnetic equilibrium. The second loop dynamically emerges from the photosphere in the same neighborhood. There are a large number of possible configurations in a two-loop system regarding their relative orientation, physical size, and directions of their toroidal magnetic field and electric current. We present three representative, but characteristically different, configurations whose interactions result in releasing various amounts of energy stored in the magnetic field. Using typical coronal parameters, some of them can take place in a relatively short timescale and release sufficient energy to account for a small flare.

A THREE-DIMENSIONAL MODEL OF ACTIVE REGION 7986: COMPARISON OF SIMULATIONS WITH OBSERVATIONS

In the present study, we use a forward modeling method to construct a 3D thermal structure encompassing active region 7986 of 1996 August. The extreme ultraviolet (EUV) emissions are then computed and compared with observations. The heating mechanism is inspired by a theory on Alfven wave turbulence dissipation. The magnetic structure is built from a Solar and Heliospheric Observatory (SOHO)/MDI magnetogram and an estimated torsion parameter deduced from observations. We found that the solution to the equations in some locations is in a thermal nonequilibrium state. The time variation of the density and temperature profiles leads to time dependent emissions, which appear as thin, loop-like structures with uniform cross-section. Their timescale is consistent with the lifetime of observed coronal loops. The dynamic nature of the solution also leads to plasma flows that resemble observed coronal rain. The computed EUV emissions from the coronal part of the fan loops and the high loops compare favorably with SOHO/EIT observations in a quantitative comparison. However, the computed emission from the lower atmosphere is excessive compared to observations, a symptom common to many models. Some factors for this discrepancy are suggested, including the use of coronal abundances to compute the emissions and the neglect of atmospheric opacity effects.

VERIFICATION OF CORONAL LOOP DIAGNOSTICS USING REALISTIC THREE-DIMENSIONAL HYDRODYNAMIC MODELS

Many different techniques have been used to characterize the plasma in the solar corona: density-sensitive spectral line ratios are used to infer the density, the evolution of coronal structures in different passbands is used to infer the temperature evolution, and the simultaneous intensities measured in multiple passbands are used to determine the emission measure distributions. All these analysis techniques assume that the intensity of the structures can be isolated through background subtraction. In this paper, we use simulated observations from a three-dimensional hydrodynamic simulation of a coronal active region to verify these diagnostics. The density and temperature from the simulation are used to generate images in several passbands and spectral lines. We identify loop structures in the simulated images and calculate the background. We then determine the density, temperature, and emission measure distribution as a function of time from the observations and compare these with the true temperature and density of the loop. We find that the overall characteristics of the temperature, density, and emission measure are recovered by the analysis methods, but the details are not. For instance, the emission measure curves calculated from the simulated observations are much broader than the true emission measure distribution, though the average temperature evolution is similar. These differences are due, in part, to a limitation of the analysis methods, but also to inadequate background subtraction.