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
All sourcesNASA NTRSDOE OSTI
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)↗
NASA NTRS1999

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

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

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

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

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

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

Cluster Analysis of IRIS Spectroscopic Line Profiles and SDO/AIA EUV Emission in Observations and RMHD Simulations of the Solar Atmosphere

Spatially-resolved observations from the IRIS and SDO/AIA satellites, especially when coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar chromosphere, transition region, and corona. However, the complexity of the data makes understanding the observations and modeling results difficult. In this work, we apply unsupervised clustering algorithms for analysis of observational and synthetic chromospheric Mg II h&k 2796Å&2803Å and transition region C II 1334Å&1335Å line profiles observed by IRIS, and extreme ultraviolet (EUV) emission observed by SDO/AIA, for various types of problems. The synthetic line profiles are computed for simulations of the quiescent solar atmosphere (using the StellarBox and RH1.5 codes). The K-Means clustering algorithm is applied, and the selection of an optimal number of clusters is supported by the average silhouette width technique. We discuss applications of the line profile clustering method to 1) visualization of computational and observational spectroscopic imaging data; 2) understanding of evolutionary trends and behavior patterns of quiet Sun emission and during solar flares; and 3) recognition of heating events and shock waves.

Sadykov, Viacheslav↗
NASA NTRSDate not supplied

Nanoflare heating in coronal X-ray Bright Points

Understanding the heating of the non-flaring solar corona is an active topic of research in heliophysics. It is well accepted that magnetic fields are mainly responsible for coronal heating. The photospheric driver randomly moves the foot-points of the magnetic field lines, and either generates waves or the quasi-static buildup of magnetic energy, depending on the timescale of motion. Heating by the dissipation of the magnetic energy is termed as DC heating while the dissipation of wave is known as AC heating mechanism. Both the AC and DC heating mechanisms can lead to impulsive heating events, termed nanoflares. The magnitude and frequency of these nanoflares determine whether they can adequately satisfy the coronal heating budget. Thus it is of great importance to study the nanoflares properties to validate their role in coronal heating. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred through indirect observational techniques combined with simulated nanoflare-heated plasma emissions. We combined the observed emission properties of coronal X-ray bright points (XBP) in EUV and X-ray wavelengths with the simulated nanoflare heated plasma to understand nanoflare properties and their contribution in coronal heating. Our results indicate that the heating of coronal XBPs could be explained by frequent nanoflares. Further we demonstrate that the sensitive spectroscopic observations in X-ray wavelength are crucial to diagnose the nanoflare properties.

Biswajit Mondal↗
NASA NTRS2016

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗
NASA NTRS2024

Direct Imaging of Magnetohydrodynamic Wave Mode Conversion Near A 3D Null Point on the Sun

Mutual conversion of various kinds of magnetohydrodynamic (MHD) waves can have profound impacts on wave propagation, energy transfer, and heating of the solar chromosphere and corona. Mode conversion occurs when an MHD wave travels through a region where the Alfvén and sound speeds are equal (e.g., a 3D magnetic null point). Here we report the direct extreme ultraviolet (EUV) imaging of mode conversion from a fast-mode to a slow-mode MHD wave near a 3D null point using Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) observations. An incident fast EUV wavefront associated with an adjacent eruptive flare propagates laterally through a neighboring pseudostreamer. Shortly after the passage of the fast EUV wave through the null point, a slow-mode wave appears near the null that propagates upward along the open structures and simultaneously downward along the separatrix encompassing the fan loops of the pseudostreamer base. These observations suggest the existence of mode conversion near 3D nulls in the solar corona, as predicted by theory and MHD simulations. Moreover, we observe decaying transverse oscillations in both the open and closed structures of the pseudostreamer, along with quasiperiodic type III radio bursts indicative of repetitive episodes of electron acceleration.

Pankaj Kumar↗
NASA NTRS2011

Widespread Nanoflare Variability Detected with Hinode/XRT in a Solar Active Region

It is generally agreed that small impulsive energy bursts called nanoflares are responsible for at least some of the Sun s hot corona, but whether they are the explanation for most of the multi-million degree plasma has been a matter of ongoing debate. We here present evidence that nanoflares are widespread in an active region observed by the X-Ray Telescope on-board the Hinode mission. The distributions of intensity fluctuations have small but important asymmetries, whether taken from individual pixels, multi-pixel subregions, or the entire active region. Negative fluctuations (corresponding to reduced intensity) are greater in number but weaker in amplitude, so that the median fluctuation is negative compared to a mean of zero. Using Monte Carlo simulations, we show that only part of this asymmetry can be explained by Poisson photon statistics. The remainder is explainable with a tendency for exponentially decreasing intensity, such as would be expected from a cooling plasma produced, e.g., from a nanoflare. We suggest that nanoflares are a universal heating process within active regions.

Reale, Fabio↗
NASA NTRS1999

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar Corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and Momentum sources are required to produce high speed wind from Corona] holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature above the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UltraViolet Coronagraph Spectrometer instrument (UVCS), and with the Ulysses/Solar Wind Observations Over the Poles of the Sun instrument (SWOOPS) proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 and 5 solar radii (2 and 5 R(sub S)) is similar to the density reported from SPARTAN 201.-01 measurements by Fisher and Guhathakurta [19941. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer the temperature and density are similar to those reported empirically by Li et al. [1998], and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub S), as it is in all other MHD coronal streamer models [e.g., Steinolfson et al., 1982; also G. A. Gary and D. Alexander, Constructing the coronal magnetic field, submitted to Solar Physics, 1998].

Suess, S. T.↗
NASA NTRS2018

Investigating the Response of Loop Plasma to Nanoflare Heating Using RADYN Simulations

We present the results of 1D hydrodynamic simulations of coronal loops that are subject to nanoflares, caused by either in situ thermal heating or nonthermal electron (NTE) beams. The synthesized intensity and Doppler shifts can be directly compared with Interface Region Imaging Spectrograph (IRIS) and Atmospheric Imaging Assembly (AIA) observations of rapid variability in the transition region (TR) of coronal loops, associated with transient coronal heating. We find that NTEs with high enough low-energy cutoff (EC) deposit energy in the lower TR and chromosphere, causing blueshifts (up to approximately 20 kilometers per second) in the IRIS Si IV lines, which thermal conduction cannot reproduce. The EC threshold value for the blueshifts depends on the total energy of the events (approximately 5 kiloelectronvolts for 1024 ergs, up to 15 kiloelectronvolts for 1025 ergs). The observed footpoint emission intensity and flows, combined with the simulations, can provide constraints on both the energy of the heating event and EC. The response of the loop plasma to nanoflares depends crucially on the electron density: significant Si IV intensity enhancements and flows are observed only for initially low-density loops (less than 10 (sup 9) per cubic centimeter). This provides a possible explanation of the relative scarcity of observations of significant moss variability. While the TR response to single heating episodes can be clearly observed, the predicted coronal emission (AIA 94 Angstroms) for single strands is below current detectability and can only be observed when several strands are heated closely in time. Finally, we show that the analysis of the IRIS Mg II chromospheric lines can help further constrain the properties of the heating mechanisms.

Sun: Corona↗
NASA NTRS2022

Solar Flare Heating with Turbulent Suppression of Thermal Conduction

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function

Solar flares↗
NASA NTRS1998

Beta in Streamers

Streamers are generally described as regions of the corona in which the density is higher than in coronal holes because the plasma is trapped by closed loops of magnetic flux. In contrast, MHD models of the global corona show the plasma beta in streamers above siml.2R_S (heliocentric). There are three recent contributions to this topic. The first is that heating near the cusp further drives beta up and results in release of new slow solar wind plasma from the top of the streamer. The second is SOHO/UVCS observations, in combination with a potential field/source surface model of the magnetic field, that show beta above 1.2R_S in a streamer observed near solar sunspot minimum. The third is a magnetic field reconstruction technique that uses field deforming algorithms and is more versatile for local fields than potential field models . The field reconstruction algorithm was applied to an isolated active region (AR 7999) and to the Pneuman Kopp global MHD model (beta has never been published for their model). In the active region, beta becomes larger than unity at siml.2 R_S. In the Pneuman & Kopp model, beta at the base of the streamer and rises with increasing height, becoming 15-20 at 1.6R_S and 35- 50 at 1.7R_S. Global simulations go on to show that the reason streamers do not simply explode under such high beta conditions is that they are held down by pressure from the sides due to the magnetic fields (and low beta) in adjacent coronal holes. The main role of the closed magnetic loop near the cusp is to keep the steamer from continuously leaking plasma, as otherwise happens in a magnetic pinch which is similar but has no closed loops. Awareness of MHD physical conditions in streamers is causing us to focus more attention on the details of the heating. On obvious suggestion is that heating is at much lower heights in streamers than in coronal holes. Also, energy which is directly delivered as momentum in coronal holes might all be deposited as heat in streamers.

Suess, Steven T.↗
NASA NTRS1998

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and momentum sources are required to produce high speed wind from coronal holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature in the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UVCS, and with the Ulysses/SWOOPS proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 solar radii and 5 solar radii (2RS and 5RS) is similar to the density reported from SPARTAN 201-01 measurements by Fisher and Guhathakurta. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer, the temperature and density are similar to those reported empirically by Li et al and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub s), as it is in all other MHD coronal streamer models.

Suess, Steven T.↗
NASA NTRS1992

Diagnostics of electron-heated solar flare models. III - Effects of tapered loop geometry and preheating

A series of hydrodynamic numerical simulations of nonthermal electron-heated solar flare atmospheres and their corresponding soft X-ray Ca XIX emission-line profiles, under the conditions of tapered flare loop geometry and/or a preheated atmosphere, is presented. The degree of tapering is parameterized by the magnetic mirror ratio, while the preheated atmosphere is parameterized by the initial upper chromospheric pressure. In a tapered flare loop, it is found that the upward motion of evaporated material is faster compared with the case where the flare loop is uniform. This is due to the diverging nozzle seen by the upflowing material. In the case where the flare atmosphere is preheated and the flare geometry is uniform, the response of the atmosphere to the electron collisional heating is slow. The upward velocity of the hydrodynamic gas is reduced due not only to the large coronal column depth, but also to the increased inertia of the overlying material. It is concluded that the only possible electron-heated scenario in which the predicted Ca XIX line profiles agree with the BCS observations is when the impulsive flare starts in a preheated dense corona.

Emslie, A. G.↗