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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Two dimensional PIC simulations of plasma heating by the dissipation of Alfven waves

Two dimensional plasma particle simulations of the evolution of large amplitude circularly polarized Alfven waves propagating parallel to the magnetic field show that the waves decay via both one- and two- dimensional parametric decay instabilities. For parameters studied, one-dimensional processes dominate the simulations, but two-dimensional decay processes, including the recently predicted filamentation instability are also observed. The daughter waves generated by the parametric decay are primarily damped by the ions, leading to ion heating. The parametric decay processes efficiently convert the ordered fluid ion motion in the Alfven wave into ion thermal energy. These processes may be important for the dissipation of Alfven waves in the solar wind, the corona and other space plasma environments. The computations were performed on the Intel Touchstone parallel supercomputer.

Liewer, P. C.↗
NASA NTRS1993

A model of an X-ray-illuminated accretion disk and corona

The X-ray-illuminated surface of the accretion disk in a low-mass X-ray Binary (LMXRB) and the X-ray-heated corona above the disk produce optical, UV, and soft X-ray emission lines. This paper presents 1D models of the emission line spectra and the vertical temperature and density structures at different radii. The models include a detailed treatment of the important atomic processes and an escape probability treatment of radiative transfer. Soker and Raymond (1993) use the density structure predicted by these models for a 2D Monte Carlo simulation of the photon scattering in the accretion disk corona (ADC) to examine the effects of the ADC on the angular distribution of X-rays and the flux of X-rays incident on the outer disk. This paper concentrates on the emission line fluxes for various elemental abundances and disk parameters. The UV lines of the classic LMXRBs are consistent with the model predictions. Some CNO processing is necessary to account for the nitrogen and helium abundances in Sco X-1 and other LMXRBs. Comparison of the models with observed spectra also points to a soft X-ray component with luminosity comparable to the hard X-rays. The models predict a substantial luminosity in the group of highly ionized iron lines near 100 A.

Raymond, John C.↗
NASA NTRS2008

The Role of Magnetic Reconnection in Solar Activity

The central challenge in solar/heliospheric physics is to understand how the emergence and transport of magnetic flux at the photosphere drives the structure and dynamics that we observe in the corona and heliosphere. This presentation focuses on the role of magnetic reconnection in determining solar/heliospheric activity. We demonstrate that two generic properties of the photospheric magnetic and velocity fields are responsible for the ubiquitous reconnection in the corona. First, the photospheric velocities are complex, which leads to the injection of energy and helicity into the coronal magnetic fields and to the efficient, formation of small-scale structure. Second, the flux distribution at the photosphere is multi-polar, which implies that topological discontinuities and, consequently, current sheets, must be present in the coronal magnetic field. We: present numerical simulations showing that photospherically-driven reconnection is responsible for the heating and dynamics of coronal plasma, and for the topology of the coronal/heliospheric magnetic field.

Antiochos, Spiro↗
NASA NTRS2022

Flows in Enthalpy-based Thermal Evolution of Loops

Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space-integrated 0D simulations. The enthalpy-based thermal evolution of loops (EBTEL) framework is a commonly used method to study the exchange of mass and energy between the corona and transition region. EBTEL solves for density, temperature, and pressure, averaged over the coronal part of the loop, velocity at coronal base, and the instantaneous differential emission measure distribution in the transition region. The current single-fluid version of the code, EBTEL2, assumes that at all stages the flows are subsonic. However, sometimes the solutions show the presence of supersonic flows during the impulsive phase of heat input. It is thus necessary to account for this effect. Here, we upgrade EBTEL2 to EBTEL3 by including the kinetic energy term in the Navier–Stokes equation. We compare the solutions from EBTEL3 with those obtained using EBTEL2, as well as the state-of-the-art field-aligned hydrodynamics code HYDRAD. We find that the match in pressure between EBTEL3 and HYDRAD is better than that between EBTEL2 and HYDRAD. Additionally, the velocities predicted by EBTEL3 are in close agreement with those obtained with HYDRAD when the flows are subsonic. However, EBTEL3 solutions deviate substantially from HYDRAD’s when the latter predicts supersonic flows. Using the mismatches in the solution, we propose a criterion to determine the conditions under which EBTEL can be used to study flows in the system.

Solar coronal heating↗
NASA NTRS2022

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗
NASA NTRS1992

Volcanism by melt-driven Rayleigh-Taylor instabilities and possible consequences of melting for admittance ratios on Venus

A large number of volcanic features exist on Venus, ranging from tens of thousands of small domes to large shields and coronae. It is difficult to reconcile all these with an explanation involving deep mantle plumes, since a number of separate arguments lead to the conclusion that deep mantle plumes reaching the base of the lithosphere must exceed a certain size. In addition, the fraction of basal heating in Venus' mantle may be significantly lower than in Earth's mantle reducing the number of strong plumes from the core-mantle boundary. In three-dimensional convection simulations with mainly internal heating, weak, distributed upwellings are usually observed. We present an alternative mechanism for such volcanism, originally proposed for the Earth and for Venus, involving Rayleigh-Taylor instabilities driven by melt buoyancy, occurring spontaneously in partially or incipiently molten regions.

Tackley, P. J.↗
NASA NTRS2003

Microphysics of Waves and Instabilities in the Solar Wind and their Macro Manifestations in the Corona and Interplanetary Space

Investigations of the physical processes responsible for the acceleration of the solar wind were pursued with the development of two new solar wind codes: a hybrid code and a 2-D MHD code. Hybrid simulations were performed to investigate the interaction between ions and parallel propagating low frequency ion cyclotron waves in a homogeneous plasma. In a low-beta plasma such as the solar wind plasma in the inner corona, the proton thermal speed is much smaller than the Alfven speed. Vlasov linear theory predicts that protons are not in resonance with low frequency ion cyclotron waves. However, non-linear effect makes it possible that these waves can strongly heat and accelerate protons. This study has important implications for study of the corona and the solar wind. Low frequency ion cyclotron waves or Alfven waves are commonly observed in the solar wind. Until now, it is believed that these waves are not able to heat the solar wind plasma unless some cascading processes transfer the energy of these waves to high frequency part. However, this study shows that these waves may directly heat and accelerate protons non-linearly. This process may play an important role in the coronal heating and the solar wind acceleration, at least in some parameter space.

Habbal, Shadia R.↗
NASA NTRS1994

Accretion disk coronae in high-luminosity systems

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer n two dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for L/L(sub Edd) approx. greater than 0.1, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.

Murray, Stephen D.↗
NASA NTRS2010

Can Thermal Nonequilibrium Explain Coronal Loops?

Any successful model of coronal loops must explain a number of observed properties. For warm (approx. 1 MK) loops, these include: 1. excess density, 2. flat temperature profile, 3. super-hydrostatic scale height, 4. unstructured intensity profile, and 5. 1000-5000 s lifetime. We examine whether thermal nonequilibrium can reproduce the observations by performing hydrodynamic simulations based on steady coronal heating that decreases exponentially with height. We consider both monolithic and multi-stranded loops. The simulations successfully reproduce certain aspects of the observations, including the excess density, but each of them fails in at least one critical way. -Xonolithic models have far too much intensity structure, while multi-strand models are either too structured or too long-lived. Storms of nanoflares remain the only viable explanation for warm loops that has been proposed so far. Our results appear to rule out the widespread existence of heating that is both highly concentrated low in the corona and steady or quasi-steady (slowly varying or impulsive with a rapid cadence). Active regions would have a very different appearance if the dominant heating mechanism had these properties. Thermal nonequilibrium may nonetheless play an important role in prominences and catastrophic cooling e(veen.gts..,coronal rain) that occupy a small fraction of the coronal volume. However, apparent inconsistencies between the models and observations of cooling events have yet to be understood.

Klimchuk, James A.↗
NASA NTRS2017

Kinetic Scale Structure of Low-Frequency Waves and Fluctuations

The dissipation of solar wind turbulence at kinetic scales is believed to be important for the heating of the corona and for accelerating the wind. The linear Vlasov kinetic theory is a useful tool for identifying various wave modes, including kinetic Alfven, fast magnetosonic/whistler, and ion-acoustic (or kinetic slow), and their possible roles in the dissipation. However, the kinetic mode structure in the vicinity of ion-cyclotron modes is not clearly understood. The present paper aims to further elucidate the structure of these low-frequency waves by introducing discrete particle effects through hybrid simulations and Klimontovich formalism of spontaneous emission theory. The theory and simulation of spontaneously emitted low-frequency fluctuations are employed to identify and distinguish the detailed mode structures associated with ion-Bernstein modes versus quasi-modes. The spontaneous emission theory and simulation also confirm the findings of the Vlasov theory in that the kinetic Alfven waves can be defined over a wide range of frequencies, including the proton cyclotron frequency and its harmonics, especially for high-beta plasmas. This implies that these low-frequency modes may play predominant roles even in the fully kinetic description of kinetic scale turbulence and dissipation despite the fact that cyclotron harmonic and Bernstein modes may also play important roles in wave-particle interactions.

waves↗
NASA NTRS1998

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

Our LTSA grant supports a long-term collaborative investigation of stellar activity. The project involves current NASA spacecraft and supporting ground-based telescopes, will make use of future missions, and utilizes the extensive archives of IUE, ROSAT, HST, and EUVE. Our interests include observational work (with a nonnegligible groundbased component); specialized processing techniques for imaging and spectral data; and semiempirical modeling, ranging from optically-thin emission measure studies to simulations of optically-thick resonance lines. Collaborations with our cool-star colleagues here in Boulder (at JILA and the High Altitude Observatory) provide access to even broader expertise, particularly on the solar corona, convection, and magnetohydrodynamic phenomena (including "dynamo" theories). The broad-brush of our investigation include the following: (1) where do coronae occur in the Hertzsprung-Russell diagram? (2) the winds of coronal stars: hot, cool, or both? (3) age, activity, rotation relations; (4) atmospheric inhomogeneities; and (5) heating mechanisms, subcoronal flows and flares. Our observation task has been to map the global properties of chromospheres and coronae in the H-R diagram and conduct detailed studies of key objects.

Ayres, Thomas R.↗
NASA NTRS2009

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2003

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRS1991

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.↗
NASA NTRS2013

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↗
NASA NTRS2014

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)↗