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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗
NASA NTRSDate not supplied

Nanoflare Heating of an X-Ray Bright Point

Nanoflares are thought to be one of the prime candidates that can keep the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing the nanoflare heated simulated plasma emissions with the observed emission. Here, we present a simulation of emission from an X-ray Bright Point (XBP) that was observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from SDO/AIA and Hinode/XRT. We use EBTEL hydrodynamic code to simulate the XBP loops. Length and magnetic field strength of these loops are derived from the potential field extrapolation of the observed photospheric magnetogram by HMI/SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the looplength and magnetic field strength. The simulated outputs are used to predict the intensity of spectrally pure map of Fe-18, Fe-17, Ne-9 ,O-8, O-9, Ne-9 etc, which are then compared with the derived intensity from MaGIXS observation. Further we have predicted the intensity map as observed by AIA and XRT and compared them with the observation. We also estimated the temperature distribution of the XBP from the simulation and found a good agreement with the derived distribution from MaGIXS observation.

coronal heating↗
NASA NTRS1989

Chromospheric-coronal coupling during solar flares: Current systems and particle acceleration

Two-dimensional (three velocity) electrostatic particle simulations are used to investigate the particle heating and acceleration associated with the impulsive phase of a solar flare. A crossfield current in the high corona (which is presumably driven by reconnection processes) is used to initiate the flare. Due to the differential motion of the electrons and ions, currents, and associated quasi-static electric fields are generated with the primary current and balancing return current being on adjacent field lines. These currents extend from the corona down into the chromosphere. Electrons can be accelerated to energies exceeding 100 keV on short time scales via the quasi-static fields and wave-particle interactions. The spectra of these electrons has a broken power-law distribution which hardens in time. The spatially separate primary and return currents are closed by the cross-field acceleration of the ambient ions into the primary current regions. These ions are then accelerated upwards into the corona by the same quasi-static electric field accelerating the electrons downwards. This acceleration can account for the broadened stationary and weak blue shifted component seen in soft x ray line emissions and enhancements in heavy ion abundances seen in the solar wind in associations with solar flares.

Winglee, Robert M.↗
NASA NTRS1973

Skylab and solar exploration

Review of some of the findings concerning solar structure, energy production, and heat transport obtained with the aid of the manned Skylab space station observatory launched on May 14, 1973. Among the topics discussed are the observation of thermonuclear fusion processes which cannot be simulated on earth, the observation of short-wave solar radiation not visible to observers on earth, and the investigation of energy-transport processes occurring in the photosphere, chromosphere, and corona. An apparent paradox is noted in that the cooler chromosphere is heating the hotter corona, seemingly in defiance of the second law of thermodynamics, thus suggesting that a nonthermal mechanism underlies the energy transport. Understanding of this nonthermal mechanism is regarded as an indispensable prerequisite for future development of plasma systems for terrestrial applications.

Von Puttkamer, J.↗
NASA NTRS1980

Mass ejections

Observations and model simulations of solar mass ejection phenomena are examined in an investigation of flare processes. Consideration is given to Skylab and other observations of flare-associated sprays, eruptive prominences, surges and coronal transients, and to MHD, gas dynamic and magnetic loop models developed to account for them. Magnetic forces are found to confine spray material, which originates in preexisting active-region filaments, within steadily expanding loops, while surges follow unmoving, preexisting magnetic field lines. Simulations of effects of a sudden pressure pulse at the bottom of the corona are found to exhibit many characteristics of coronal transients associated with flares, and impulsive heating low in the chromosphere is found to be able to account for surges. The importance of the magnetic field as the ultimate source of energy which drives eruptive phenomena as well as flares is pointed out.

Rust, D. M.↗
NASA NTRS2000

MACS, An Instrument, and a Methodology for Simulations and Global Measurements of the Coronal Electron Temperature and the Solar Wind Velocity on the Solar Corona

The determination of the radial and latitudinal temperature and wind profiles of the solar corona is of great importance in understanding the coronal heating mechanism and the dynamics of coronal expansion. Cram presented the theory for the formation of the K-coronal spectrum and identified two important observations. He observed the existence of temperature sensitive anti-nodes at certain wavelengths in the theoretical K-coronal spectra. The anti-nodes are separated by temperature-insensitive nodes. Remarkably, Cram showed that the wavelengths of the nodes and anti-nodes are almost independent of altitude above the solar limb. Because of these features, Cram suggested that the intensity ratios at two anti-nodes could be used as a diagnostic of the electron temperature in the K-corona. Based on this temperature diagnostic technique prescribed by Cram a slit-based spectroscopic study was performed by Ichimoto et al. on the solar corona in conjunction with the total solar eclipse of 3 Nov 1994 in Putre, Chile to determine the temperature profile of the solar corona. In this thesis Cram's theory has been extended to incorporate the role of the solar wind in the formation of the K-corona, and we have identified both temperature and wind sensitive intensity ratios. The instrument, MACS, for Multi Aperture Coronal Spectrometer, a fiber optic based spectrograph, was designed for global and simultaneous measurement of the thermal electron temperature and the solar wind velocity in the solar corona. The first ever experiment of this nature was conducted in conjunction with the total solar eclipse of 11 Aug 1999 in Elazig, Turkey. In this instrument one end of each of twenty fiber optic tips were positioned in the focal plane of the telescope in such a way that we could observe conditions simultaneously at many different latitudes and two different radial distances in the solar corona. The other ends of the fibers were vertically aligned and placed at the primary focus of the collimating lens of the spectrograph to obtain simultaneous and global spectra on the solar corona. By isolating the K-coronal spectrum from the spectrum recorded by each fiber the temperature and the wind sensitive intensity ratios were calculated to obtain simultaneous and global measurements of the thermal electron temperature and the solar wind velocity. We were successful in obtaining reliable estimates of the coronal temperature at many positions in the corona. This is the first time that simultaneous measurements of coronal temperatures have been obtained at so many points. However due to instrumental scattering encountered during observations, reliable estimates of the wind velocity turned out to be impossible to obtain. Although remedial measures were taken prior to observation, this task proved to be difficult owing to the inability to replicate the conditions expected during an eclipse in the laboratory. The full extent of the instrumental scattering was apparent only when we analyzed the observational sequence. Nevertheless the experience obtained from this very first attempt to simultaneously and globally measure both the wind velocity and the temperature on the solar corona have provided valuable information to conduct any future observations successfully.

Reginald, Nelson L.↗
NASA NTRS2014

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗
NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗
NASA NTRS2004

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗
NASA NTRS2023

Role of Small-Scale Impulsive Events in Heating the X-Ray Bright Points of the Quiet Sun

Small-scale impulsive events, known as nanoflares, are thought to be one of the prime candidates that can keep the solar corona hot at its multimillion-Kelvin temperature. Individual nanoflares are difficult to detect with the current generation of instruments; however, their presence can be inferred through indirect techniques such as Differential Emission Measure (DEM) analysis. Here, we employ this technique to investigate the possibility of nanoflare heating of the quiet corona during the minimum of solar cycle 24. We estimate the DEM of disk-integrated quiet Sun and X-ray bright points (XBP) using the observations from XSM on board the Chandrayaan-2 orbiter and AIA on board the Solar Dynamic Observatory. XBPs are found to be the dominant contributor to disk-integrated X-rays, with a radiative flux of ∼2 × 10 5 erg cm −2 s −1 . XBPs consist of small-scale loops associated with bipolar magnetic fields. We simulate such XBP loops using the EBTEL hydrodynamic code. The lengths and magnetic field strengths of these loops are obtained through a potential field extrapolation of the photospheric magnetogram. Each loop is assumed to be heated by random nanoflares having an energy that depends on the loop properties. The composite nanoflare energy distribution for all the loops has a power-law slope close to −2.5. The simulation output is then used to obtain the integrated DEM. It agrees remarkably well with the observed DEM at temperatures above 1 MK, suggesting that the nanoflare distribution, as predicted by our model, can explain the XBP heating.

Solar coronal heating↗
NASA NTRS1994

Structure and Dynamics of the Solar Corona

Advanced computational techniques were used to study solar coronal heating and coronal mass ejections. A three dimensional, time dependent resistive magnetohydrodynamic code was used to study the dynamic response of a model corona to continuous, slow, random magnetic footpoint displacements in the photosphere. Three dimensional numerical simulations of the response of the corona to simple smooth braiding flows in the photosphere were calculated to illustrate and understand the spontaneous formation of current filaments. Two dimensional steady state helmet streamer configurations were obtained by determining the time asymptotic state of the interaction of an initially one dimensinal transponic solar wind with a spherical potential dipole field. The disruption of the steady state helmet streamer configuration was studied as a response to shearing of the magnetic footpoints of the closed field lines under the helmet.

Schnack, D. D.↗
NASA NTRS2008

Explaining Warm Coronal Loops

One of the great mysteries of coronal physics that has come to light in the last few years is the discovery that warn (- 1 INK) coronal loops are much denser than expected for quasi-static equilibrium. Both the excess densities and relatively long lifetimes of the loops can be explained with bundles of unresolved strands that are heated impulsively to very high temperatures. Since neighboring strands are at different stages of cooling, the composite loop bundle is multi-thermal, with the distribution of temperatures depending on the details of the "nanoflare storm." Emission hotter than 2 MK is predicted, but it is not clear that such emission is always observed. We consider two possible explanations for the existence of over-dense warm loops without corresponding hot emission: (1) loops are bundles of nanoflare heated strands, but a significant fraction of the nanoflare energy takes the form of a nonthermal electron beam rather then direct plasma heating; (2) loops are bundles of strands that undergo thermal nonequilibrium that results when steady heating is sufficiently concentrated near the footpoints. We present numerical hydro simulations of both of these possibilities and explore the observational consequences, including the production of hard X-ray emission and absorption by cool material in the corona.

Klimchuk, James A.↗
NASA NTRS1999

What We Don't Understand About Ion Acceleration Flares

There are now strong associations between the (3)He-rich, Fe-rich ions in "impulsive" solar energetic particle (SEP) events and the similar abundances derived from gamma-ray lines from flares. Compact flares, where wave energy can predominate, are ideal sites for the study of wave-particle physics. Yet there are nagging questions about the magnetic geometry, the relation between ions that escape and those that interact, and the relative roles of cascading Alfven waves and the EMIC waves required to enhance He-3. There are also questions about the relative timing of ion and electron acceleration and of heating; these relate to the variation of ionization states before and during acceleration and during transport out of the corona. We can construct a model that addresses many of these issues, but problems do remain. Our greatest lack is realistic theoretical simulations of element abundances, spectra, and their variations. By contrast, we now have a much better idea of the acceleration at CME-driven shock waves in the rare but large "gradual" SEP events, largely because of their slow temporal evolution and great spatial extent.

Reames, Donald V.↗
NASA NTRS2011

What Dominates the Coronal Emission Spectrum During the Cycle of Impulsive Heating and Cooling?

The smoking gun of small-scale, impulsive events heating the solar corona is expected to be the presence of a hot ( > 5 MK) plasma component. Evidence for this has been scarce, but has gradually begun to accumulate due to recent studies designed to constrain the high temperature part of the emission measure distribution. However, the detected hot component is often weaker than models predict and this is due in part to the common modeling assumption that the ionization balance remains in equilibrium. The launch of the latest generation of space-based observing instrumentation aboard Hinode and the Solar Dynamics Observatory (SDO) has brought the matter of the ionization state of the plasma firmly to the forefront. It is timely to consider exactly what emission current instruments would detect when observing a corona heated impulsively on small-scales by nanoflares. Only after we understand the full effects of nonequilibrium ionization can we draw meaningful conclusions about the plasma that is (or is not) present. We have therefore performed a series of hydrodynamic simulations for a variety of different nanoflare properties and initial conditions. Our study has led to several key conclusions. 1. Deviations from equilibrium are greatest for short-duration nanoflares at low initial coronal densities. 2. Hot emission lines are the most affected and are suppressed sometimes to the point of being invisible. 3. The emission detected in all of the SDO-AIA channels is generally dominated by warm, over-dense, cooling plasma. 4. It is difficult not to create coronal loops that emit strongly at 1.5 MK and in the range 2 to 5 MK, which are the most commonly observed kind, for a broad range of nanoflare scenarios. 5. The Fe XV (284.16 ) emission in most of our models is about 10 times brighter than the Ca XVII (192.82 ) emission, consistent with observations. Our overarching conclusion is that small-scale, impulsive heating inducing a nonequilibrium ionization state leads to predictions for observable quantities that are entirely consistent with what is actually observed.

Bradshaw, Stephen J.↗
NASA NTRSDate not supplied

Constraining Nanoflare Frequency Through Observations and Models

Nanoflares are considered prime candidates for heating the non-flaring solar corona. However, direct detection of individual nanoflares with present-generation instruments remains challenging. Understanding the frequency and magnitude of nanoflares is crucial for comprehending their role in coronal heating. In this study, we employ field-aligned hydrodynamic model to simulate the emission from a non-flaring X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS). The length and magnetic field strength of the coronal loops associated with this XBP are derived from magnetic field extrapolation of observed photospheric magnetograms by SDO/HMI. Each loop is assumed to be heated by random nanoflares, whose frequency and magnitude are determined by the loop length and their field strength. The simulation results are then compared and matched against the observation to determine the nanoflare heating frequency. We will discuss the findings of this study. Furthermore, we will discuss the evolution of nanoflare heating frequency during the formation and decay of an active region (AR) on the solar disk.

coronalheating↗
NASA NTRS2020

Transition Region Contribution to AIA Observations in the Context of Coronal Heating

We investigate the ratio of coronal and transition region intensity in coronal loops observed by the AtmosphericImaging Assembly(AIA)on the Solar Dynamics Observatory(SDO). Using Enthalpy-based Thermal Evolution ofLoops(EBTEL)hydrodynamic simulations, we model loops with multiple lengths and energyfluxes heatedrandomly by events drawn from power-law distributions with different slopes and minimum delays between eventsto investigate how each of these parameters influences observable loop properties. We generate AIA intensitiesfrom the corona and transition region for each realization. The variations within and between models generatedwith these different parameters illustrate the sensitivity of narrowband imaging to the details of coronal heating.We then analyze the transition region and coronal emission from a number of observed active regions andfindbroad agreement with the trends in the models. In both models and observations, the transition region brightness issignificant, often greater than the coronal brightness in all six“coronal”AIA channels. We also identify an inverserelationship, consistent with heating theories, between the slope of the differential emission measure(DEM)coolward of the peak temperature and the observed ratio of coronal to transition region intensity. These resultshighlight the use of narrowband observations and the importance of properly considering the transition region ininvestigations of coronal heating

S J Schonfeld↗
NASA NTRS1998

The Paradox of Filamented Coronal Hole Flow but Uniform High Speed Wind

Plumes and rays in coronal holes are nearly radially aligned density striations that follow the ambient magnetic field. They have long been known, but have gained new interest with growing awareness that coronal hole flow is inherently filamentary. In retrospect, filamentary flow should have been no surprise. This is because,Beta much less than 1 in coronal holes inside approximately 10 Solar radius, allowing the flow to be filamentary down to the smallest scale of photospheric magnetic activity. While the magnetic field itself is locally smooth across any height above ca. 50,000 km, SOHO/MDI has shown that the photospheric magnetic field is a complex array of rapidly evolving small bipoles that are constantly emerging, evolving, and cancelling. The resulting activity is manifested in microflares, concentrated in the magnetic network, that produce Impulsive injections at the footpoints of coronal field lines. The uneven distribution of this activity in space and time is the source of coronal hole filamentation. What is surprising is that the radial flow speed also exhibits filamentary structure. It is not well described as smooth, spherically symmetric, diverging flow, but instead ranges from 300 to over 1000 km/s at 5.5 Solar radius among field-aligned filaments like those seen in plumes and rays [Feldman et al., JGR, Dec. 1997]. This is completely unlike the constant high speed solar wind reported beyond 0.3 AU. Consequently, plumes and filamentary structure must be strongly mixed, and the mixing must be far along by 0.3 AU to be consistent with Helios observations. The paradox is what causes the mixing? Existing models of coronal heating and solar wind acceleration hardly address this issue. One possibility we are investigating is the MHD Kelvin-Helmholtz instability, to which the shear between plumes and interplume corona is expected to become unstable at 5-10 Solar radius. This instability can be simulated and followed far into the nonlinear regime and may lead to Alfvenic fluctuations like those seen at 1 AU.

Suess, Steven T.↗
NASA NTRS2015

An MHD Code for the Study of Magnetic Structures in the Solar Wind

We have developed a 2.5D MHD code designed to study how the solar wind influences the evolution of transient events in the solar corona and inner heliosphere. The code includes thermal conduction, coronal heating and radiative cooling. Thermal conduction is assumed to be magnetic field-aligned in the inner corona and transitions to a collisionless formulation in the outer corona. We have developed a stable method to handle field-aligned conduction around magnetic null points. The inner boundary is placed in the upper transition region, and the mass flux across the boundary is determined from 1D field-aligned characteristics and a 'radiative energy balance' condition. The 2.5D nature of this code makes it ideal for parameter studies not yet possible with 3D codes. We have made this code publicly available as a tool for the community. To this end we have developed a graphical interface to aid in the selection of appropriate options and a graphical interface that can process and visualize the data produced by the simulation. As an example, we show a simulation of a dipole field stretched into a helmet streamer by the solar wind. Plasmoids periodically erupt from the streamer, and we perform a parameter study of how the frequency and location of these eruptions changed in response to different levels of coronal heating. As a further example, we show the solar wind stretching a compact multi-polar flux system. This flux system will be used to study breakout coronal mass ejections in the presence of the solar wind.

corona↗