Search NASA⌕ Search

SEARCH · Search NASA

Results for “corona simulation heating”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

Flare Ribbon Signatures of Reconnection Plasmoids

Solar flares are explosive space weather events that rapidly convert stored magnetic energy into bulk motion, plasma heating, and particle acceleration. Understanding the structure and dynamics of the magnetic reconnection that powers flares is critical for predicting the energy release. In particular, the amount of energy transferred to energetic particles is thought to be highly dependent on whether the reconnection is primarily turbulent (e.g., plasmoid dominated) or instead laminar. We present high-resolution MHD simulations of three-dimensional reconnection in an eruptive flare and compare the results to recent data. Although flare reconnection is challenging to observe directly in the corona, highly detailed constraints on its dynamics can be obtained from observations of flare ribbons that track the chromospheric footpoints of newly reconnected field lines. The analogues of flare ribbons in our simulations are identified by tracking discontinuous changes in field-line magnetic connectivity due to the reconnection. In our highest-resolution calculations, we find that these ribbon analogues are highly structured and exhibit many ‘whorl’ patterns that are linked to turbulent plasmoids in the reconnecting current sheet. Flare ribbon fine structure therefore reveals crucial information about the fundamental turbulent vs. laminar nature of the reconnection. We discuss the implications of these results for understanding reconnection-driven energy release throughout the heliosphere.

J. T. Dahlin↗

Three-Dimensional MHD Modeling of The Solar Corona and Solar Wind: Comparison with The Wang-Sheeley Model

We present simulation results from a tilted-dipole steady-state MHD model of the solar corona and solar wind and compare the output from our model with the Wang-Sheeley model which relates the divergence rate of magnetic flux tubes near the Sun (inferred from solar magnetograms) to the solar wind speed observed near Earth and at Ulysses. The boundary conditions in our model specified at the coronal base and our simulation region extends out to 10 AU. We assumed that a flux of Alfven waves with amplitude of 35 km per second emanates from the Sun and provides additional heating and acceleration for the coronal outflow in the open field regions. The waves are treated in the WKB approximation. The incorporation of wave acceleration allows us to reproduce the fast wind measurements obtained by Ulysses, while preserving reasonable agreement with plasma densities typically found at the coronal base. We find that our simulation results agree well with Wang and Sheeley's empirical model.

Usmanov, A. V.↗

Flare Ribbon Signatures of Plasmoids

Solar flares are explosive space weather events that rapidly convert stored magnetic energy into bulk motion, plasma heating, and particle acceleration. Understanding the structure and dynamics of the magnetic reconnection that powers flares is critical for predicting the energy release. In particular, the amount of energy transferred to energetic particles is thought to be highly dependent on whether the reconnection is primarily turbulent (e.g., plasmoid dominated) or instead laminar. We present new high-resolution MHD simulations of three-dimensional reconnection in an eruptive flare and compare the results to recent data. Although flare reconnection is challenging to observe directly in the corona, highly detailed constraints on its dynamics can be obtained from observations of flare ribbons that track the chromospheric footpoints of newly reconnected field lines. The analogues of flare ribbons in our simulations are identified by tracking discontinuous changes in field-line magnetic connectivity due to the reconnection. In our highest-resolution calculations, we find that these ribbon analogues are highly structured and exhibit many ‘whorl’ patterns that are linked to turbulent plasmoids in the reconnecting current sheet. Such flare ribbon fine structure therefore reveals crucial information about the fundamental turbulent vs. laminar nature of the reconnection. Plasmoid-driven particle acceleration models would therefore predict co-location of energetic particle signatures with plasmoid-associated ribbon fine structure.

J T Dahlin↗

Coronal current-sheet formation - The effect of asymmetric and symmetric shears

A 2.5D numerical code is used to investigate the results of an asymmetric shear imposed on a potential quadrupolar magnetic field under two sets of atmospheric boundary conditions - a low-beta plasma with line tying at the base, similar to the line-tied analytic model, and a hydrostatic-equilibrium atmosphere with solar gravity, typical of the observed photosphere-chromosphere interface. The low-beta simulation confirms the crucial role of the line-tying assumption in producting current sheets. The effects of a symmetric shear on the same hydrostatic-equilibrium atmosphere is examined, using more grid points to improve the resolution of the current structures which form along the flux surfaces. It is found that true current sheets do not form in the corona when a more realistic model is considered. The amount of Ohmic dissipation in the thick currents is estimated to be two to four orders of magnitude below that required to heat the corona. It is concluded that magnetic topologies of the type examined here do not contribute significantly to coronal heating.

Karpen, Judith T.↗

Prominence Mass Supply and the Cavity

A prevalent but untested paradigm is often used to describe the prominence-cavity system; the cavity is under-dense because it it evacuated by supplying mass to the condensed prominence. The thermal non-equilibrium (TNE) model of prominence formation offers a theoretical framework to predict the thermodynamic evolutin of the prominence and the surrounding corona. We examine the evidence for a prominence-cavity connection by comparing the TNE model and diagnostics of dynamic extreme ultraviolet (EUV) emission surrounding the prominence, specifically prominence horns. Horns are correlated extensions of prminence plasma and coronal plasma which appear to connect the prominence and cavity. The TNE model predicts that large-scale brightenings will occur in the Solar Dynamics Observatory Atmospheric Imaging Assembly 171 A badpass near he prominence that are associated with the cooling phase of condensation formation. In our simulations, variations in the magnitude of footpoint heating lead to variations in the duration, spatial scale, and temporal offset between emission enhancements in the other EUV bandpasses. While these predictions match well a subset of the horn observations, the range of variations in the observed structures is not captured by the model. We discuss the implications of one-dimensional loop simulations for the three-dimensional time-averaged equilibrium in the prominence and the cavity. Evidence suggests that horns are likely caused by condensing prominence plasma, but the larger question of whether this process produces a density-depleted cavity requires a more tightly constrained model of heating and better knowledge of the associated magnetic structure.

PROMINENCES↗

Numerical simulations of a siphon mechanism for quiescent prominence formation

Quiescent prominences represent a significant challenge to our understanding of the flow of mass and energy in the outer layers of the solar atmosphere. A small number of quiescent prominences contain as much mass as the entire corona (Athay, 1976). The problem then is how to get that much material into the relatively small volume of a prominence and maintain it at a temperature of 10,000 K in close proximity to material at one million K. The thermal insulation to conduction provided by the magnetic field explains the disparate temperatures. The mass source problem is less well understood. One method for supplying mass to the prominence is to siphon it from the chromosphere. The siphon mechanism begins with a magnetic loop that evolves into a configuration with a gravitational well, such as that described by Kippenhahn and Schluter (1957). This could be formed, for example, by a twist in the magnetic field. A gravitational well could also be formed by a condensation induced sag in the field. This could further enhance the condensation process. Once this well has formed, or as it is forming, the material in the well area of the loop must cool and condense to the point where radiative losses exceed any heat input. Additional material must also flow into the well from the underlying chromosphere to supply the mass required to form the prominence. One example from a series of numerical simulations that were performed to study the formation of quiescent prominences is presented.

Poland, A. I.↗

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗

Generation of Electron Suprathermal Tails in the Upper Solar Atmosphere: Implications for Coronal Heating

We present a mechanism for the generation of non-Maxwellian electron distribution function in the upper regions of the solar atmosphere in the presence of collisional damping. It is suggested that finite amplitude, low frequency, obliquely propagating electromagnetic waves can carry a substantial electric field component parallel to the mean magnetic field that can be significantly larger than the Dreicer electric field. This long wavelength electric fluctuation is capable of generating high frequency electron plasma oscillations and low frequency ion acoustic-like waves. The analysis has been performed using 1-1/2D Vlasov and PIC numerical simulations in which both electrons and ions are treated kinetically and self consistently. The simulation results indicate that high frequency electron plasma oscillations and low frequency ion acoustic-like waves are generated. The high frequency electron plasma oscillation drives electron plasma turbulence, which subsequently is damped out by the background electrons. The turbulence damping results in electron acceleration and the generation of non-Maxwellian suprathermal tails on time scales short compared to collisional damping. Bulk heating also occurs if the fluctuating parallel electric field is strong enough. This study suggests that finite amplitude, low frequency, obliquely propagating, electromagnetic waves can play a significant role in the acceleration and heating of the solar corona electrons and in the coupling of medium and small-scale phenomena.

Vinas, Adolfo F.↗

Solar Flare X-Ray Source Motion as a Response to Electron Spectral Hardening

Context. Solar flare hard X-rays (HXRs) are thought to be produced by nonthermal coronal electrons stopping in the chromosphere, or remaining trapped in the corona. The collisional thick target model (CTTM) predicts that more energetic electrons penetrate to greater column depths along the flare loop. This requires that sources produced by harder power-law injection spectra should appear further down the legs or footpoints of a flareloop. Therefore, the frequently observed hardening of the injected power-law electron spectrum during flare onset should be concurrent with a descending hard X-ray source. Aims. To test this implication of the CTTM by comparing its predicted HXR source locations with those derived from observations of a solar flare which exhibits a nonthermally-dominated spectrum before the peak in HXRs, known as an early impulsive event. Methods. HXR images and spectra of an early impulsive C-class flare were obtained using the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). Images were reconstructed to produce HXR source height evolutions for three energy bands. Spatially-integrated spectral analysis was performed to isolate nonthermal emission, and to determine the power-law index of the electron injection spectrum. The observed height-time evolutions were then fit with CTTM-based simulated heights for each energy, using the electron spectral indices derived from the RHESSI spectra. Results. The flare emission was found to be dominantly nonthermal above approx. 7 keV, with emission of thermal and nonthermal X-rays likely to be simultaneously observable below that energy. The density structure required for a good match between model and observed source heights agreed with previous studies of flare loop densities. Conclusions. The CTTM has been used to produce a descent of model HXR source heights that compares well with observations of this event. Based on this interpretation, downward motion of nonthermal sources should indeed occur in any flare where there is spectral hardening in the electron distribution during a flare. However, this would often be masked by thermal emission associated with flare plasma pre-heating. As yet, flare models that predict transfer of energy from the corona to the chromosphere by means other than a flux of nonthermal electrons do not predict this observed source descent. Therefore, flares such as this will be key in explaining this elusive energy transfer process.

X-rays↗

Modeling solar flare conduction fronts. I - Homogeneous plasmas and ion-acoustic turbulence. II - Inhomogeneous plasmas and ambipolar electric fields

A one-dimensional, electrostatic, particle-in-cell simulation is used here to model the expansion of a heated electron population in a coronal loop during a solar flare and the characteristics of the associated X-ray emissions. The hot electrons expand outward from the localized region, creating an ambipolar electric field which accelerates a return current of cooler, ambient electrons. Ion-acoustic waves are generated by the return currents as proposed by Brown et al. (1979), but they play little or no role in containing energetic electrons and the conduction front proposed by Brown et al. does not form. The X-ray emission efficiency of the electrons is too low in the corona for them to be the source of hard X-ray bursts. The particle dynamics changes dramatically if the heated plasma is at low altitudes and expands upward into the more tenuous plasma at higher altitudes. Two important applications of this finding are the radio-frequency heating of the corona and the collisional heating of the chromosphere by precipitating energetic electrons. In both cases, the overlying plasma has a density that is too low to supply a balancing return current to the expanding hot electrons. As a result, an ambipolar electric field develops that tends to confine the energetic electrons behind a front that propagate outward at about the speed of sound.

Mckean, M. E.↗

Structure and dynamics of the coronal magnetic field

The last few years have seen a marked increase in the sophistication of models of the solar corona. This has been brought about by a confluence of three key elements. First, the collection of high-resolution observations of the Sun, both in space and time, has grown tremendously. The SOHO (Solar Heliospheric Observatory) mission is providing additional correlated high-resolution magnetic, white-light and spectroscopic observations. Second, the power and availability of supercomputers has made two- and three-dimensional modeling routine. Third, the sophistication of the models themselves, both in their geometrical realism and in the detailed physics that has been included, has improved significantly. The support from our current Space Physics Theory grant has allowed us to exploit this confluence of capabilities. We have carried out direct comparisons between observations and models of the solar corona. The agreement between simulated coronal structure and observations has verified that the models are mature enough for detailed analysis, as we will describe. The development of this capability is especially timely, since observations obtained from three space missions that are underway (Ulysses, WIND and SOHO) offer an opportunity for significant advances in our understanding of the corona and heliosphere. Through this interplay of observations and theory we can improve our understanding of the Sun. Our achievements thus far include progress modeling the large-scale structure of the solar corona, three-dimensional models of active region fields, development of emerging flux and current, formation and evolution of coronal loops, and coronal heating by current filaments.

VanHoven, Gerard↗

The Fundamental Physical Processes Producing and Controlling Stellar Coronal/Transition Region/Chromospheric Activity and Structure

Our LTSA (Long Term Space Astrophysics) research has utilized current NASA and ESA spacecraft, supporting ground-based IR, radio, and sub-mm telescopes, and the extensive archives of HST (Hubble Space Telescope), IUE (International Ultraviolet Explorer), ROSAT, EUVE (Extreme Ultraviolet Explorer), and other missions. Our research effort has included observational work (with a nonnegligible groundbased component), specialized processing techniques for imaging and spectral data, and semiempirical modelling, ranging from optically thin emission measure studies to simulations of optically thick resonance lines. In our previous LTSA efforts, we have had a number of major successes, including most recently: organizing and carrying out an extensive cool star UV survey in HST cycle eight; obtaining observing time with new instruments, such as Chandra and XMM (X-ray Multi-Mirror) in their first cycles; collaborating with the Chandra GTO program and participating with the Chandra Emission Line Project on multi-wavelength observations of HR 1099 and Capella. These are the main broad-brush themes of our previous investigation: a) Where do Coronae Occur in the Hertzsprung-Russell Diagram? b) Winds of Coronal and Noncoronal Stars; c) Activity, Age, Rotation Relations; d) Atmospheric Inhomogeneities; e) Heating Mechanisms, Subcoronal Flows, and Flares; f) Development of Analysis and Modelling Tools.

Ayres, T. R.↗

Experimental Study of Alfvén Wave Reflection from an Alfvén-speed Gradient Relevant to the Solar Coronal Holes

Abstract We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of the wavelength of the incident wave to the length scale of the gradient. Two-fluid simulations using the Gkeyll code qualitatively agree with and support the experimental findings. Our experimental results support models of wave heating that rely on wave reflection at low heights from a smooth Alfvén-speed gradient to drive turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Prominence Eruption Initiated by Helical Kink Instability of an Embedded Flux Rope

We study the triggering mechanism of a limb-prominence eruption and the associated coronal mass ejection (CME) near AR 12342 using Solar Dynamics Observatory and Large Angle and Spectrometric Coronagraph/Solar Heliospheric Observatory observations. The prominence is seen with an embedded flux thread (FT) at one end and bifurcates from the middle to a different footpoint location. The morphological evolution of the FT is similar to that of an unstable flux rope (FR), which we regard as a prominence-embedded FR. The FR twist exceeds the critical value. In addition, the morphology of the prominence plasma in 304 Å images marks the helical nature of the magnetic skeleton, with a total of 2.96 turns along arc length. The potential field extrapolation model indicates that the critical height of the background magnetic field gradient falls within the inner corona (105 Mm), which is consistent with the extent of coronal plasma loops. These results suggest that the helical kink instability in the embedded FR caused the slow rise of the prominence to the height of the torus instability domain. Moreover, the differential emission measure analysis unveils heating of the prominence plasma to coronal temperatures during an eruption, suggesting reconnection-related heating underneath the upward rising embedded FR. The prominence starts with a slow rise motion of 10 km s(exp. -1), which is followed by fast and slow acceleration phases that have an average acceleration of 28.9 m s(exp. -2) and 2.4 m s(exp. -2) in C2 and C3 field of view, respectively. As predicted by previous numerical simulations, the observed synchronous kinematic profiles of the CME leading edge and the core support the involved FR instability in the prominence initiation.

Vemareddy, P.↗

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗

Heating and acceleration of ions with Kappa distribution functions by low‐frequency Alfvén wave

Abstract Heating and acceleration of ions with Kappa distribution functions (with parameter ) in low‐beta plasmas, by a low‐frequency Alfvén wave, is investigated using test‐particle simulations, yielding interesting new results. As long as the Alfvén wave amplitude is sufficiently large, the computed net heating energy of ions becomes independent of the wave frequency and amplitude, always approaching the same value of . The eventual energy of ions is dictated only by the initial ion energy and the ratio of the magnetic field energy density to the plasma density. The heating effect of the Kappa ions increases with . During the heating process, the ions are picked up by the Alfvén wave and pitch angle scattered, forming a quasi‐isotropic spherical shell velocity distribution. The Kappa ions are accelerated in the parallel direction, reaching a bulk flow speed roughly equal to the local Alfvén speed. Higher ‐value in the initial Kappa distribution leads to faster saturation. The above results may explain certain features of the ion heating and acceleration in the solar wind and corona.

Li, Kehua↗