-
Evolution of the autoresonant plasma wave excitation in two-dimensional particle-in-cell simulations
Authors:
Mufei. Luo,
Caterina. Riconda,
Anna. Grassi,
Ning. Wang,
Jonathan S. Wurtele,
Tünde Fülöp,
István Pusztai
Abstract:
The generation of an autoresonantly phase-locked high amplitude plasma waves to the chirped beat frequency of two driving lasers is studied in two dimensions using particle-in-cell simulations. The two-dimensional plasma and laser parameters correspond to those that optimized the plasma wave amplitude in one-dimensional simulations. Near the start of autoresonant locking, the two-dimensional simul…
▽ More
The generation of an autoresonantly phase-locked high amplitude plasma waves to the chirped beat frequency of two driving lasers is studied in two dimensions using particle-in-cell simulations. The two-dimensional plasma and laser parameters correspond to those that optimized the plasma wave amplitude in one-dimensional simulations. Near the start of autoresonant locking, the two-dimensional simulations appear similar to one-dimensional particle-in-cell results [Luo et al., Phys. Rev. Res. 6, 013338 (2024)] with plasma wave amplitudes above the Rosenbluth-Liu limit. Later, just below wave-breaking, the two-dimensional simulation exhibits a Weibel-like instability and eventually laser beam filamentation. These limit the coherence of the plasma oscillation after the peak plasma wave field is obtained. In spite of the reduction of spatial coherence of the accelerating density structure, the acceleration of self-injected electrons in the case studied remains at $70\%$ to $80\%$ of that observed in one dimension. Other effects such as plasma wave bowing are discussed.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Control of autoresonant plasma beat-wave wakefield excitation
Authors:
M. Luo,
C. Riconda,
I. Pusztai,
A. Grassi,
J. S. Wurtele,
T. Fülöp
Abstract:
Autoresonant phase-locking of the plasma wakefield to the beat frequency of two driving lasers offers advantages over conventional wakefield acceleration methods, since it requires less demanding laser parameters and is robust to variations in the target plasma density. Here, we investigate the kinetic and nonlinear processes that come into play during autoresonant plasma beat-wave acceleration of…
▽ More
Autoresonant phase-locking of the plasma wakefield to the beat frequency of two driving lasers offers advantages over conventional wakefield acceleration methods, since it requires less demanding laser parameters and is robust to variations in the target plasma density. Here, we investigate the kinetic and nonlinear processes that come into play during autoresonant plasma beat-wave acceleration of electrons, their impact on the field amplitude of the accelerating structure, and on acceleration efficiency. Particle-in-Cell simulations show that the process depends on the plasma density in a non-trivial way but can be reliably modeled under specific conditions. Beside recovering previous fluid results in the deeply underdense plasma limit, we demonstrate that robust field excitation can be achieved within a fully kinetic self-consistent modeling. By adjusting the laser properties, we can amplify the electric field to the desired level, up to wave-breaking, and efficiently accelerate particles; we provide suggestions for optimized laser and plasma parameters. This versatile and efficient acceleration scheme, producing electrons from tens to hundreds of MeV energies, holds promise for a wide range of applications in research industry and medicine.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
Fluid and kinetic studies of tokamak disruptions using Bayesian optimization
Authors:
Ida Ekmark,
Mathias Hoppe,
Tünde Fülöp,
Patrik Jansson,
Liam Antonsson,
Oskar Vallhagen,
Istvan Pusztai
Abstract:
When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for…
▽ More
When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for the Compton scattering and tritium beta decay runaway electron generation mechanisms in our simulation tool DREAM [M. Hoppe et al 2021 Comp. Phys. Commun. 268, 108098]. To achieve a diverse set of disruption scenarios, Bayesian optimization is used to explore a range of massive material injection densities for deuterium and neon. The cost function is designed to distinguish between successful and unsuccessful disruption mitigation based on the runaway current, current quench time and transported fraction of the heat loss. In the non-activated scenarios, we find that fluid and kinetic disruption simulations can have significantly different runaway electron dynamics, due to an overestimation of the runaway seed by the fluid model. The primary cause of this is that the fluid hot-tail generation model neglects superthermal electron transport losses during the thermal quench. In the activated scenarios, the fluid and kinetic models give similar predictions, which can be explained by the significant influence of the activated sources on the RE dynamics and the seed.
△ Less
Submitted 5 July, 2024; v1 submitted 8 February, 2024;
originally announced February 2024.
-
Runaway Electron Dynamics in ITER Disruptions with Shattered Pellet Injections
Authors:
Oskar Vallhagen,
Lise Hanebring,
Javier Artola,
Michael Lehnen,
Eric Nardon,
Tünde Fülöp,
Mathias Hoppe,
Sarah Newton,
Istvan Pusztai
Abstract:
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well a…
▽ More
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well as including additional effects, such as ion transport and enhanced runaway electron transport during the thermal quench. The work aims to provide a fairly comprehensive coverage of experimentally feasible scenarios, considering plasmas representative of both non-activated and high-performance DT operation, different thermal quench onset criteria and transport levels, a wide range of hydrogen and neon quantities injected in one or two stages, and pellets with various characteristic shard sizes. Using a staggered injection scheme, with a pure hydrogen injection preceding a mixed hydrogen-neon injection, we find injection parameters leading to acceptable runaway electron currents in all investigated discharges without activated runaway sources. Dividing the injection into two stages is found to significantly enhance the assimilation and minimize runaway electron generation due to the hot-tail mechanism. However, while a staggered injection outperforms a single stage injection also in cases with radioactive runaway electron sources, no cases with acceptable runaway electron currents are found for a DT-plasma with a 15 MA plasma current.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Bayesian optimization of massive material injection for disruption mitigation in tokamaks
Authors:
Istvan Pusztai,
Ida Ekmark,
Hannes Bergström,
Peter Halldestam,
Patrik Jansson,
Mathias Hoppe,
Oskar Vallhagen,
Tünde Fülöp
Abstract:
A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide in…
▽ More
A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide insights into the behaviour of the disruption dynamics for representative scenarios. The simulations are conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We show that irrespective of the quantities of the material deposition, multi-megaampere runaway currents will be produced in the deuterium-tritium phase of operations, even in the optimal scenarios. However, the severity of the outcome can be influenced by tailoring the radial profile of the injected material; in particular if the injected neon is deposited at the edge region it leads to a significant reduction of both the final runaway current and the transported heat losses. The Bayesian approach allows us to map the parameter space efficiently, with more accuracy in favorable parameter regions, thereby providing us information about the robustness of the optima.
△ Less
Submitted 24 February, 2023; v1 submitted 2 February, 2023;
originally announced February 2023.
-
Drift of ablated material after pellet injection in a tokamak
Authors:
O. Vallhagen,
I. Pusztai,
P. Helander,
S. L. Newton,
T. Fülöp
Abstract:
Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are…
▽ More
Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are computationally expensive. Here, we present a one-dimensional semi-analytical model for the radial displacement of ablated material after pellet injection, taking into account both the Alfvén and ohmic currents which short-circuit the charge separation creating the drift. The model is suitable for rapid calculation of the radial drift displacement, and can be useful for e.g. modelling of disruption mitigation via pellet injection.
△ Less
Submitted 24 April, 2023; v1 submitted 30 January, 2023;
originally announced January 2023.
-
The impact of fusion-born alpha particles on runaway electron dynamics in ITER disruptions
Authors:
Andrej Lier,
Gergely Papp,
Philipp Lauber,
Istvan Pusztai,
Konsta Särkimäki,
Ola Embreus
Abstract:
In the event of a tokamak disruption in a D-T plasma, fusion-born alpha particles take several milliseconds longer to thermalise than the background. As the damping rates drop drastically following the several orders of magnitudes drop of temperature, Toroidal Alfven Eigenmodes (TAEs) can be driven by alpha particles in the collapsing plasma before the onset of the current quench. We employ kineti…
▽ More
In the event of a tokamak disruption in a D-T plasma, fusion-born alpha particles take several milliseconds longer to thermalise than the background. As the damping rates drop drastically following the several orders of magnitudes drop of temperature, Toroidal Alfven Eigenmodes (TAEs) can be driven by alpha particles in the collapsing plasma before the onset of the current quench. We employ kinetic simulations of the alpha particle distribution and show that the TAEs can reach sufficiently strong saturation amplitudes to cause significant core runaway electron transport in unmitigated ITER disruptions. As the eigenmodes do not extend to the plasma edge, this effect leads to an increase of the runaway electron plateau current. Mitigation via massive material injection however changes the Alfven frequency and can lead to mode suppression. A combination of the TAE-caused core runaway electron transport with other perturbation sources could lead to a drop of runaway current in unmitigated disruptions.
△ Less
Submitted 16 March, 2023; v1 submitted 18 January, 2023;
originally announced January 2023.
-
On the minimum transport required to passively suppress runaway electrons in SPARC disruptions
Authors:
R. A. Tinguely,
I. Pusztai,
V. A. Izzo,
K. S'"arkimäki,
T. Fülöp,
D. T. Garnier,
R. S. Granetz,
M. Hoppe,
C. Paz-Soldan,
A. Sundström,
R. Sweeney
Abstract:
In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2…
▽ More
In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1 MA compared to ~8.7 MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.3 and with kinetic energies 0.2-15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 $\mathrm{m^2/s}$, is needed in the core to fully mitigate these REs. Properly accounting for (i) the CQ electric field's effect on RE transport in islands and (ii) the contribution of significant RE currents to disruption MHD may help achieve this.
△ Less
Submitted 3 January, 2023;
originally announced January 2023.
-
Runaway dynamics in reactor-scale spherical tokamak disruptions
Authors:
E. Berger,
I. Pusztai,
S. L. Newton,
M. Hoppe,
O. Vallhagen,
A. Fil,
T. Fülöp
Abstract:
Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runawa…
▽ More
Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runaway production. The study is conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We find that, in many cases, mitigation strategies are necessary to prevent the runaway current from reaching multi-megaampere levels. Our results indicate that with a suitably chosen deuterium-neon mixture for mitigation, it is possible to achieve a tolerable runaway current and ohmic current evolution. With such parameters, however, the majority of the thermal energy loss happens through radial transport rather than radiation, which poses a risk of unacceptable localised heat loads.
△ Less
Submitted 3 August, 2022;
originally announced August 2022.
-
Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil
Authors:
V. A. Izzo,
I. Pusztai,
K. Särkimäki,
A. Sundström,
D. Garnier,
D. Weisberg,
R. A. Tinguely,
C. Paz-Soldan,
R. S. Granetz,
R. Sweeney
Abstract:
The operation of a 3D coil--passively driven by the current quench loop voltage--for the deconfinement of runaway electrons is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear MHD modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code…
▽ More
The operation of a 3D coil--passively driven by the current quench loop voltage--for the deconfinement of runaway electrons is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear MHD modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for runaway electrons based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely, et al [2021 Nucl. Fusion 61 124003], considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the $n=1$ coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the TQ MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1-2MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the current quench is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Stimulated-Raman-scattering amplification of attosecond XUV pulses with pulse-train pumps and application to local in-depth plasma-density measurement
Authors:
Andréas Sundström,
Mickael Grech,
István Pusztai,
Caterina Riconda
Abstract:
We present a scheme for amplifying an extreme-ultraviolet (XUV) seed isolated attosecond pulse via stimulated Raman scattering of a pulse-train pump. At sufficient seed and pump intensity, the amplification is nonlinear, and the amplitude of the seed pulse can reach that of the pump, one order of magnitude higher than the initial seed amplitude. In the linear amplification regime, we find that the…
▽ More
We present a scheme for amplifying an extreme-ultraviolet (XUV) seed isolated attosecond pulse via stimulated Raman scattering of a pulse-train pump. At sufficient seed and pump intensity, the amplification is nonlinear, and the amplitude of the seed pulse can reach that of the pump, one order of magnitude higher than the initial seed amplitude. In the linear amplification regime, we find that the spectral signature of the pump pulse train is imprinted on the spectrum of the amplified seed pulse. Since the spectral signature is imprinted with its frequency downshifted by the plasma frequency, it is possible to deduce the electron density in the region of interaction. This region can be of micrometer length scale longitudinally. By varying the delay between the seed and the pump, this scheme provides a local electron-density measurement inside solid-density plasmas that cannot be probed with optical frequencies, with micrometer resolution.
△ Less
Submitted 21 October, 2022; v1 submitted 14 July, 2022;
originally announced July 2022.
-
Runaway dynamics in disruptions with current relaxation
Authors:
István Pusztai,
Mathias Hoppe,
Oskar Vallhagen
Abstract:
The safe operation of tokamak reactors requires a reliable modeling capability of disruptions, and in particular the spatio-temporal dynamics of associated runaway electron currents. In a disruption, instabilities can break up magnetic surfaces into chaotic field line regions, causing current profile relaxation, as well as a rapid radial transport of heat and particles. Using a mean-field helicity…
▽ More
The safe operation of tokamak reactors requires a reliable modeling capability of disruptions, and in particular the spatio-temporal dynamics of associated runaway electron currents. In a disruption, instabilities can break up magnetic surfaces into chaotic field line regions, causing current profile relaxation, as well as a rapid radial transport of heat and particles. Using a mean-field helicity transport model implemented in the disruption runaway modeling framework DREAM, we calculate the dynamics of runaway electrons in the presence of current relaxation events. In scenarios where flux surfaces remain intact in parts of the plasma, a skin current is induced at the boundary of the intact magnetic field region. This skin current region becomes an important center concerning the subsequent dynamics: It may turn into a hot ohmic current channel, or a sizable radially localized runaway beam, depending on the heat transport. If the intact region is in the plasma edge, runaway generation in the counter-current direction can occur, which may develop into a sizable reverse runaway beam. Even when the current relaxation extends to the entire plasma, the final runaway current density profile can be significantly affected, as the induced electric field is reduced in the core and increased in the edge, thereby shifting the center of runaway generation towards the edge.
△ Less
Submitted 2 June, 2022;
originally announced June 2022.
-
Attosecond dispersion as a diagnostics tool for solid-density laser-generated plasmas
Authors:
Andréas Sundström,
István Pusztai,
Per Eng-Johnsson,
Tünde Fülöp
Abstract:
Extreme-ultraviolet pulses can propagate through ionised solid-density targets, unlike optical pulses and, thus, have the potential to probe the interior of such plasmas on sub-femtosecond timescales. We present a synthetic diagnostic method for solid-density laser-generated plasmas based on the dispersion of an extreme-ultraviolet attosecond probe pulse, in a pump--probe scheme. We demonstrate th…
▽ More
Extreme-ultraviolet pulses can propagate through ionised solid-density targets, unlike optical pulses and, thus, have the potential to probe the interior of such plasmas on sub-femtosecond timescales. We present a synthetic diagnostic method for solid-density laser-generated plasmas based on the dispersion of an extreme-ultraviolet attosecond probe pulse, in a pump--probe scheme. We demonstrate the theoretical feasibility of this approach through calculating the dispersion of an extreme-ultraviolet probe pulse propagating through a laser-generated plasma. The plasma dynamics is calculated using a particle-in-cell simulation, whereas the dispersion of the probe is calculated with an external pseudo-spectral wave solver, allowing for high accuracy when calculating the dispersion. The application of this method is illustrated on thin-film plastic and aluminium targets irradiated by a high-intensity pump pulse. By comparing the dispersion of the probe pulse at different delays relative to the pump pulse, it is possible to follow the evolution of the plasma as it disintegrates. The high-frequency end of the dispersion provides information on the line-integrated electron density, whereas lower frequencies are more affected by the highest density encountered along the path of the probe. In addition, the presence of thin-film interference could be used to study the evolution of the plasma surface.
△ Less
Submitted 26 April, 2022; v1 submitted 1 February, 2022;
originally announced February 2022.
-
Effect of Two-Stage Shattered Pellet Injection on Tokamak Disruptions
Authors:
O. Vallhagen,
I. Pusztai,
M. Hoppe,
S. L. Newton,
T. Fülöp
Abstract:
An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hopp…
▽ More
An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hoppe et al (2021) Comp. Phys. Commun. 268, 108098]. To minimize the obtained runaway currents an optimal range of injected deuterium quantities is found. This range is sensitive to the opacity of the plasma to Lyman radiation, which affects the ionization degree of deuterium, and thus avalanche runaway generation. The two-stage injection scheme, where dilution cooling is produced by deuterium before a radiative thermal quench caused by neon, reduces both the hot-tail seed and the localized transported heat load on the wall. However, during nuclear operation, additional runaway seed sources from the activated wall and tritium make it difficult to reach tolerably low runaway currents.
△ Less
Submitted 26 January, 2022; v1 submitted 25 January, 2022;
originally announced January 2022.
-
Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC
Authors:
RA Tinguely,
VA Izzo,
DT Garnier,
A Sundström,
K Särkimäki,
O Embréus,
T Fülöp,
RS Granetz,
M Hoppe,
I Pusztai,
R Sweeney
Abstract:
The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetr…
▽ More
The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium-tritium plasmas.
△ Less
Submitted 20 October, 2021;
originally announced October 2021.
-
Spatiotemporal analysis of the runaway distribution function from synchrotron images in an ASDEX Upgrade disruption
Authors:
M. Hoppe,
L. Hesslow,
O. Embreus,
L. Unnerfelt,
G. Papp,
I. Pusztai,
T. Fülöp,
O. Lexell,
T. Lunt,
E. Macusova,
P. J. McCarthy,
G. Pautasso,
G. I. Pokol,
G. Por,
P. Svensson,
the ASDEX Upgrade team,
the EUROfusion MST1 team
Abstract:
Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quenc…
▽ More
Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quench, together with an avalanche runaway tail which has an exponentially decreasing energy spectrum. We find that, although the avalanche component carries the vast majority of the current, it is the high-energy seed remnant that dominates synchrotron emission. With insights from the fluid-kinetic simulations, an analytic model for the evolution of the runaway seed component is developed and used to reconstruct the radial density profile of the RE beam. The analysis shows that the observed change of the synchrotron pattern from circular to crescent shape is caused by a rapid redistribution of the radial profile of the runaway density.
△ Less
Submitted 5 February, 2021; v1 submitted 29 May, 2020;
originally announced May 2020.
-
Runaway dynamics in the DT phase of ITER operations in the presence of massive material injection
Authors:
O. Vallhagen,
O Embreus,
I Pusztai,
L Hesslow,
T Fülöp
Abstract:
A runaway avalanche can result in a conversion of the initial plasma current into a relativistic electron beam in high current tokamak disruptions. We investigate the effect of massive material injection of deuterium-noble gas mixtures on the coupled dynamics of runaway generation, resistive diffusion of the electric field, and temperature evolution during disruptions in the DT phase of ITER opera…
▽ More
A runaway avalanche can result in a conversion of the initial plasma current into a relativistic electron beam in high current tokamak disruptions. We investigate the effect of massive material injection of deuterium-noble gas mixtures on the coupled dynamics of runaway generation, resistive diffusion of the electric field, and temperature evolution during disruptions in the DT phase of ITER operations. We explore the dynamics over a wide range of injected concentrations and find substantial runaway currents, unless the current quench time is intolerably long. The reason is that the cooling associated with the injected material leads to high induced electric fields that, in combination with a significant recombination of hydrogen isotopes, leads to a large avalanche generation. Balancing Ohmic heating and radiation losses provides qualitative insights into the dynamics, however, an accurate modeling of the temperature evolution based on energy balance appears crucial for quantitative predictions.
△ Less
Submitted 13 July, 2020; v1 submitted 27 April, 2020;
originally announced April 2020.
-
Collisional effects on the electrostatic shock dynamics in thin-foil targets driven by an ultraintense short pulse laser
Authors:
Andréas Sundström,
Laurent Gremillet,
Evangelos Siminos,
István Pusztai
Abstract:
We numerically investigate the impact of Coulomb collisions on the ion dynamics in high-$Z$, solid density caesium hydride and copper targets, irradiated by high-intensity ($I\approx2{-}5\times10^{20}{\rm\,Wcm^{-2}}$), ultrashort (${\sim}10{\rm\,fs}$), circularly polarized laser pulses, using particle-in-cell simulations. Collisions significantly enhance electron heating, thereby strongly increasi…
▽ More
We numerically investigate the impact of Coulomb collisions on the ion dynamics in high-$Z$, solid density caesium hydride and copper targets, irradiated by high-intensity ($I\approx2{-}5\times10^{20}{\rm\,Wcm^{-2}}$), ultrashort (${\sim}10{\rm\,fs}$), circularly polarized laser pulses, using particle-in-cell simulations. Collisions significantly enhance electron heating, thereby strongly increasing the speed of a shock wave launched in the laser-plasma interaction. In the caesium hydride target, collisions between the two ion species heat the protons to ${\sim}100{-}1000{\rm\,eV}$ temperatures. However, in contrast to previous work (A.E. Turrell etal., 2015 Nat. Commun. 6, 8905), this process happens in the upstream only, due to nearly total proton reflection. This difference is ascribed to distinct models used to treat collisions in dense/cold plasmas. In the case of a copper target, ion reflection can start as a self-amplifying process, bootstrapping itself. Afterwards, collisions between the reflected and upstream ions heat these two populations significantly. When increasing the pulse duration to $60{\rm\,fs}$, the shock front more clearly decouples from the laser piston, and so can be studied without direct interference from the laser. The shock wave formed at early times exhibits properties typical of both hydrodynamic and electrostatic shocks, including ion reflection. At late times, the shock is seen to evolve into a hydrodynamic blast wave.
△ Less
Submitted 30 July, 2020; v1 submitted 25 February, 2020;
originally announced February 2020.
-
Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields
Authors:
István Pusztai,
James Juno,
Axel Brandenburg,
Jason M. TenBarge,
Ammar Hakim,
Manaure Francisquez,
Andréas Sundström
Abstract:
We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on…
▽ More
We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.
△ Less
Submitted 26 June, 2020; v1 submitted 31 January, 2020;
originally announced January 2020.
-
Fast collisional electron heating and relaxation in thin foils driven by a circularly polarized ultraintense short-pulse laser
Authors:
Andréas Sundström,
Laurent Gremillet,
Evangelos Siminos,
István Pusztai
Abstract:
The creation of well-thermalized, hot and dense plasmas is attractive for warm dense matter studies. We investigate collisionally induced energy absorption of an ultraintense and ultrashort laser pulse in a solid copper target using particle-in-cell simulations. We find that, upon irradiation by a $2\times10^{20}{\rm\,W\,cm^{-2}}$ intensity, $60{\rm\,fs}$ duration, circularly polarized laser pulse…
▽ More
The creation of well-thermalized, hot and dense plasmas is attractive for warm dense matter studies. We investigate collisionally induced energy absorption of an ultraintense and ultrashort laser pulse in a solid copper target using particle-in-cell simulations. We find that, upon irradiation by a $2\times10^{20}{\rm\,W\,cm^{-2}}$ intensity, $60{\rm\,fs}$ duration, circularly polarized laser pulse, the electrons in the collisional simulation rapidly reach a well-thermalized distribution with ${\sim}3.5{\rm\,keV}$ temperature, while in the collisionless simulation the absorption is several orders of magnitude weaker. Circular polarization inhibits the generation of suprathermal electrons, while ensuring efficient bulk heating through inverse bremsstrahlung, a mechanism usually overlooked at relativistic laser intensity. An additional simulation, taking account of both collisional and field ionization, yields similar results: the bulk electrons are heated to ${\sim}2.5{\rm\,keV}$, but with a somewhat lower degree of thermalization than in the pre-set, fixed-ionization case. The collisional absorption mechanism is found to be robust against variations in the laser parameters. At fixed laser pulse energy, increasing the pulse duration rather than the intensity leads to a higher electron temperature.
△ Less
Submitted 28 April, 2020; v1 submitted 21 November, 2019;
originally announced November 2019.
-
The importance of the classical channel in the impurity transport of optimized stellarators
Authors:
S. Buller,
A. Mollén,
S. L. Newton,
H. M. Smith,
I. Pusztai
Abstract:
In toroidal magnetic confinement devices, such as tokamaks and stellarators, neoclassical transport is usually an order of magnitude larger than its classical counterpart. However, when a high-collisionality species is present in a stellarator optimized for low Pfirsch-Schlüter current, its classical transport can be comparable to the neoclassical transport. In this letter, we compare neoclassical…
▽ More
In toroidal magnetic confinement devices, such as tokamaks and stellarators, neoclassical transport is usually an order of magnitude larger than its classical counterpart. However, when a high-collisionality species is present in a stellarator optimized for low Pfirsch-Schlüter current, its classical transport can be comparable to the neoclassical transport. In this letter, we compare neoclassical and classical fluxes and transport coefficients calculated for Wendelstein 7-X (W7-X) and Large Helical Device (LHD) cases. In W7-X, we find that the classical transport of a collisional impurity is comparable to the neoclassical transport for all radii, while it is negligible in the LHD cases, except in the vicinity of radii where the neoclassical transport changes sign. In the LHD case, electrostatic potential variations on the flux-surface significantly enhance the neoclassical impurity transport, while the classical transport is largely insensitive to this effect in the cases studied.
△ Less
Submitted 9 June, 2019; v1 submitted 29 March, 2019;
originally announced March 2019.
-
Optimization of flux-surface density variation in stellarator plasmas with respect to the transport of collisional impurities
Authors:
Stefan Buller,
Håkan M. Smith,
Albert Mollén,
Sarah Louise Newton,
Istvan Pusztai
Abstract:
Avoiding impurity accumulation is a requirement for steady-state stellarator operation. The accumulation of impurities can be heavily affected by variations in their density on the flux-surface. Using recently derived semi-analytic expressions for the transport of a collisional impurity species with high-$Z$ and flux-surface density-variation in the presence of a low-collisionality bulk ion specie…
▽ More
Avoiding impurity accumulation is a requirement for steady-state stellarator operation. The accumulation of impurities can be heavily affected by variations in their density on the flux-surface. Using recently derived semi-analytic expressions for the transport of a collisional impurity species with high-$Z$ and flux-surface density-variation in the presence of a low-collisionality bulk ion species, we numerically optimize the impurity density-variation on the flux-surface to minimize the radial peaking factor of the impurities. These optimized density-variations can reduce the core impurity density by $0.75^Z$ (with $Z$ the impurity charge number) in the Large Helical Device case considered here, and by $0.89^Z$ in a Wendelstein 7-X standard configuration case. On the other hand, when the same procedure is used to find density-variations that maximize the peaking factor, it is notably increased compared to the case with no density-variation. This highlights the potential importance of measuring and controlling these variations in experiments.
△ Less
Submitted 12 April, 2019; v1 submitted 21 December, 2018;
originally announced December 2018.
-
Effect of a weak ion collisionality on the dynamics of kinetic electrostatic shocks
Authors:
Andréas Sundström,
James Juno,
Jason M. TenBarge,
István Pusztai
Abstract:
In strictly collisionless electrostatic shocks, the ion distribution function can develop discontinuities along phase-space separatrices, due to partial reflection of the ion population. In this paper, we depart from the strictly collisionless regime and present a semi-analytical model for weakly collisional kinetic shocks. The model is used to study the effect of small but finite collisionalities…
▽ More
In strictly collisionless electrostatic shocks, the ion distribution function can develop discontinuities along phase-space separatrices, due to partial reflection of the ion population. In this paper, we depart from the strictly collisionless regime and present a semi-analytical model for weakly collisional kinetic shocks. The model is used to study the effect of small but finite collisionalities on electrostatic shocks, and they are found to smooth out these discontinuities into growing boundary layers. More importantly, ions diffuse into and accumulate in the previously empty regions of phase space, and, by upsetting the charge balance, lead to growing downstream oscillations of the electrostatic potential. We find that the collisional age of the shock is the more relevant measure of the collisional effects than the collisionality, where the former can become significant during the lifetime of the shock, even for weak collisionalities.
△ Less
Submitted 6 February, 2019; v1 submitted 4 October, 2018;
originally announced October 2018.
-
Proton Acceleration in a Laser-induced Relativistic Electron Vortex
Authors:
Longqing Yi,
István Pusztai,
Alexander Pukhov,
Baifei Shen,
Tünde Fülöp
Abstract:
We show that when a solid plasma foil with a density gradient on the front surface is irradiated by an intense laser pulse at a grazing angle, around 80 degrees, a relativistic electron vortex is excited in the near-critical-density layer after the laser pulse depletion. The vortex structure and dynamics are studied using particle-in-cell simulations. Due to the asymmetry introduced by nonuniform…
▽ More
We show that when a solid plasma foil with a density gradient on the front surface is irradiated by an intense laser pulse at a grazing angle, around 80 degrees, a relativistic electron vortex is excited in the near-critical-density layer after the laser pulse depletion. The vortex structure and dynamics are studied using particle-in-cell simulations. Due to the asymmetry introduced by nonuniform background density, the vortex drifts at a constant velocity, typically 0.2 to 0.3 times the speed of light. The strong magnetic field inside the vortex leads to significant charge separation; in the corresponding electric field initially stationary protons can be captured and accelerated to twice the velocity of the vortex (100-200 MeV). A representative scenario - with laser intensity of 10^21 W/cm^2 -is discussed: two dimensional simulations suggest that a quasi-monoenergetic proton beam can be obtained with a mean energy 140 MeV and an energy spread of about 10%. We derive an analytical estimate for the vortex velocity in terms of laser and plasma parameters, demonstrating that the maximum proton energy can be controlled by the incidence angle of the laser and the plasma density gradient.
△ Less
Submitted 17 July, 2019; v1 submitted 14 August, 2018;
originally announced August 2018.
-
Collisional transport of impurities with flux-surface varying density in stellarators
Authors:
S. Buller,
H. M. Smith,
P. Helander,
A. Mollén,
S. L. Newton,
I. Pusztai
Abstract:
High-Z impurities in magnetic confinement devices are prone to develop density variations on the flux-surface, which can significantly affect their transport. In this paper, we generalize earlier analytic stellarator calculations of the neoclassical radial impurity flux in the mixed-collisionality regime (collisional impurity and low-collisionality bulk ions) to include the effect of such flux-sur…
▽ More
High-Z impurities in magnetic confinement devices are prone to develop density variations on the flux-surface, which can significantly affect their transport. In this paper, we generalize earlier analytic stellarator calculations of the neoclassical radial impurity flux in the mixed-collisionality regime (collisional impurity and low-collisionality bulk ions) to include the effect of such flux-surface variations. We find that only in the homogeneous density case is the transport of highly collisional impurities (in the Pfirsch-Schlüter regime) independent of the radial electric field. We study these effects for a Wendelstein 7-X (W7-X) vacuum field, with simple analytic models for the potential perturbation, under the assumption that the impurity density is given by a Boltzmann response to a perturbed potential. In the W7-X case studied, we find that larger amplitude potential perturbations cause the radial electric field to dominate the transport of the impurities. In addition, we find that classical impurity transport can be larger than the neoclassical transport in W7-X.
△ Less
Submitted 9 June, 2019; v1 submitted 2 May, 2018;
originally announced May 2018.
-
First principles of modelling the stabilization of microturbulence by fast ions
Authors:
George J. Wilkie,
Aylwin Iantchenko,
Ian G. Abel,
Edmund Highcock,
Istvan Pusztai
Abstract:
The observation that fast ions stabilize ion-temperature-gradient-driven microturbulence has profound implications for future fusion reactors. It is also important in optimizing the performance of present-day devices. In this work, we examine in detail the phenomenology of fast ion stabilization and present a reduced model which describes this effect. This model is derived from the high-energy lim…
▽ More
The observation that fast ions stabilize ion-temperature-gradient-driven microturbulence has profound implications for future fusion reactors. It is also important in optimizing the performance of present-day devices. In this work, we examine in detail the phenomenology of fast ion stabilization and present a reduced model which describes this effect. This model is derived from the high-energy limit of the gyrokinetic equation and extends the existing "dilution" model to account for nontrivial fast ion kinetics. Our model provides a physically-transparent explanation for the observed stabilization and makes several key qualitative predictions. Firstly, that different classes of fast ions, depending on their radial density or temperature variation, have different stabilizing properties. Secondly, that zonal flows are an important ingredient in this effect precisely because the fast ion zonal response is negligible. Finally, that in the limit of highly-energetic fast ions, their response approaches that of the "dilution" model; in particular, alpha particles are expected to have little, if any, stabilizing effect on plasma turbulence. We support these conclusions through detailed linear and nonlinear gyrokinetic simulations.
△ Less
Submitted 5 June, 2018; v1 submitted 2 January, 2018;
originally announced January 2018.
-
Low Mach-number collisionless electrostatic shocks and associated ion acceleration
Authors:
Istvan Pusztai,
Jason M. TenBarge,
Aletta N. Csapó,
James Juno,
Ammar Hakim,
Longqing Yi,
Tünde Fülöp
Abstract:
The existence and properties of low Mach-number ($M \gtrsim 1$) electrostatic collisionless shocks are investigated with a semi-analytical solution for the shock structure. We show that the properties of the shock obtained in the semi-analytical model can be well reproduced in fully kinetic Eulerian Vlasov-Poisson simulations, where the shock is generated by the decay of an initial density discont…
▽ More
The existence and properties of low Mach-number ($M \gtrsim 1$) electrostatic collisionless shocks are investigated with a semi-analytical solution for the shock structure. We show that the properties of the shock obtained in the semi-analytical model can be well reproduced in fully kinetic Eulerian Vlasov-Poisson simulations, where the shock is generated by the decay of an initial density discontinuity. Using this semi-analytical model, we study the effect of electron-to-ion temperature ratio and presence of impurities on both the maximum shock potential and Mach number. We find that even a small amount of impurities can influence the shock properties significantly, including the reflected light ion fraction, which can change several orders of magnitude. Electrostatic shocks in heavy ion plasmas reflect most of the hydrogen impurity ions.
△ Less
Submitted 19 December, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Isotope and density profile effects on pedestal neoclassical transport
Authors:
Stefan Buller,
Istvan Pusztai
Abstract:
Cross-field neoclassical transport of heat, particles and momentum is studied in sharp density pedestals, with a focus on isotope and profile effects, using a radially global approach. Global effects -- which tend to reduce the peak ion heat flux, and shift it outward -- increase with isotope mass for fixed profiles. The heat flux reduction exhibits a saturation with a favorable isotopic trend. A…
▽ More
Cross-field neoclassical transport of heat, particles and momentum is studied in sharp density pedestals, with a focus on isotope and profile effects, using a radially global approach. Global effects -- which tend to reduce the peak ion heat flux, and shift it outward -- increase with isotope mass for fixed profiles. The heat flux reduction exhibits a saturation with a favorable isotopic trend. A significant part of the heat flux can be convective even in pure plasmas, unlike in the plasma core, and it is sensitive to how momentum sources are distributed between the various species. In particular, if only ion momentum sources are allowed, in global simulations of pure plasmas the ion particle flux remains close to its local value, while this may not be the case for simulations with isotope mixtures or electron momentum sources. The radial angular momentum transport that is a finite orbit width effect, is found to be strongly correlated with heat sources.
△ Less
Submitted 9 August, 2017; v1 submitted 18 May, 2017;
originally announced May 2017.
-
Electromagnetic zonal flow residual responses
Authors:
Peter J. Catto,
Felix I. Parra,
Istvan Pusztai
Abstract:
The collisionless axisymmetric zonal flow residual calculation for a tokamak plasma is generalized to include electromagnetic perturbations. We formulate and solve the complete initial value zonal flow problem by retaining the fully self-consistent axisymmetric spatial perturbations in the electric and magnetic fields. Simple expressions for the electrostatic, shear and compressional magnetic resi…
▽ More
The collisionless axisymmetric zonal flow residual calculation for a tokamak plasma is generalized to include electromagnetic perturbations. We formulate and solve the complete initial value zonal flow problem by retaining the fully self-consistent axisymmetric spatial perturbations in the electric and magnetic fields. Simple expressions for the electrostatic, shear and compressional magnetic residual responses are derived that provide a fully electromagnetic test of the zonal flow residual in gyrokinetic codes. Unlike the electrostatic potential, the parallel vector potential and the parallel magnetic field perturbations need not relax to flux functions for all possible initial conditions.
△ Less
Submitted 26 June, 2017; v1 submitted 20 April, 2017;
originally announced April 2017.
-
Edge Momentum Transport by Neutrals: an Interpretive Numerical Framework
Authors:
J. T. Omotani,
S. L. Newton,
I. Pusztai,
T. Fülöp,
E. Viezzer,
the ASDEX Upgrade Team
Abstract:
Due to their high cross-field mobility, neutrals can contribute to momentum transport even at the low relative densities found inside the separatrix and they can generate intrinsic rotation. We use a charge-exchange dominated solution to the neutral kinetic equation, coupled to neoclassical ions, to evaluate the momentum transport due to neutrals. Numerical solutions to the drift-kinetic equation…
▽ More
Due to their high cross-field mobility, neutrals can contribute to momentum transport even at the low relative densities found inside the separatrix and they can generate intrinsic rotation. We use a charge-exchange dominated solution to the neutral kinetic equation, coupled to neoclassical ions, to evaluate the momentum transport due to neutrals. Numerical solutions to the drift-kinetic equation allow us to cover the full range of collisionality, including the intermediate levels typical of the tokamak edge. In the edge there are several processes likely to contribute to momentum transport in addition to neutrals. Therefore, we present here an interpretive framework that can evaluate the momentum transport through neutrals based on radial plasma profiles. We demonstrate its application by analysing the neutral angular momentum flux for an L-mode discharge in the ASDEX Upgrade tokamak. The magnitudes of the angular momentum fluxes we find here due to neutrals of up to $1{-}2\;\mathrm{N\, m}$ are comparable to the net torque on the plasma from neutral beam injection, indicating the importance of neutrals for rotation in the edge.
△ Less
Submitted 16 December, 2016;
originally announced December 2016.
-
Edge rotation from momentum transport by neutrals
Authors:
JT Omotani,
SL Newton,
I Pusztai,
T Fülöp
Abstract:
Due to their high cross field mobility, neutral atoms can have a strong effect on transport even at the low relative densities found inside the separatrix. We use a charge-exchange dominated model for the neutrals, coupled to neoclassical ions, to calculate momentum transport when it is dominated by the neutrals. We can then calculate self-consistently the radial electric field and predict the int…
▽ More
Due to their high cross field mobility, neutral atoms can have a strong effect on transport even at the low relative densities found inside the separatrix. We use a charge-exchange dominated model for the neutrals, coupled to neoclassical ions, to calculate momentum transport when it is dominated by the neutrals. We can then calculate self-consistently the radial electric field and predict the intrinsic rotation in an otherwise torque-free plasma. Using a numerical solver for the ion distribution to allow arbitrary collisionality, we investigate the effects of inverse aspect ratio and elongation on plasma rotation. We also calculate the rotation of a trace carbon impurity, to facilitate future comparison to experiments using charge exchange recombination spectroscopy diagnostics.
△ Less
Submitted 16 December, 2016;
originally announced December 2016.
-
Neoclassical flows in deuterium-helium plasma density pedestals
Authors:
Stefan Buller,
Istvan Pusztai,
Sarah Newton,
John Omotani
Abstract:
In tokamak transport barriers, the radial scale of profile variations can be comparable to a typical ion orbit width, which makes the coupling of the distribution function across flux surfaces important in the collisional dynamics. We use the radially global steady-state neoclassical δf code Perfect to calculate poloidal and toroidal flows, and radial fluxes, in the pedestal. In particular, we hav…
▽ More
In tokamak transport barriers, the radial scale of profile variations can be comparable to a typical ion orbit width, which makes the coupling of the distribution function across flux surfaces important in the collisional dynamics. We use the radially global steady-state neoclassical δf code Perfect to calculate poloidal and toroidal flows, and radial fluxes, in the pedestal. In particular, we have studied the changes in these quantities as the plasma composition is changed from a deuterium bulk species with a helium impurity to a helium bulk with a deuterium impurity, under specific profile similarity assumptions. The poloidally resolved radial fluxes are not divergence-free in isolation in the presence of sharp radial profile variations, which leads to the appearance of poloidal return-flows. These flows exhibit a complex radial-poloidal structure that extends several orbit widths into the core and is sensitive to abrupt radial changes in the ion temperature gradient. We find that a sizable neoclassical toroidal angular momentum transport can arise in the radially global theory, in contrast to the local.
△ Less
Submitted 12 January, 2017; v1 submitted 24 November, 2016;
originally announced November 2016.
-
Global anomalous transport of ICRH- and NBI-heated fast ions
Authors:
George J. Wilkie,
Istvan Pusztai,
Ian G. Abel,
William Dorland,
Tünde Fülöp
Abstract:
By taking advantage of the trace approximation, one can gain an enormous computational advantage when solving for the global turbulent transport of impurities. In particular, this makes feasible the study of non-Maxwellian transport coupled in radius and energy, allowing collisions and transport to be accounted for on similar time scales, as occurs for fast ions. In this work, we study the fully-n…
▽ More
By taking advantage of the trace approximation, one can gain an enormous computational advantage when solving for the global turbulent transport of impurities. In particular, this makes feasible the study of non-Maxwellian transport coupled in radius and energy, allowing collisions and transport to be accounted for on similar time scales, as occurs for fast ions. In this work, we study the fully-nonlinear ITG-driven trace turbulent transport of locally heated and injected fast ions. Previous results indicated the existence of MeV-range minorities heated by cyclotron resonance, and an associated density pinch effect. Here, we build upon this result using the t3core code to solve for the distribution of these minorities, consistently including the effects of collisions, gyrokinetic turbulence, and heating. Using the same tool to study the transport of injected fast ions, we contrast the qualitative features of their transport with that of the heated minorities. Our results indicate that heated minorities are more strongly affected by microturbulence than injected fast ions. The physical interpretation of this difference provides a possible explanation for the observed synergy when NBI heating is combined with ICRH. Furthermore, we move beyond the trace approximation to develop a model which allows one to easily account for the reduction of anomalous transport due to the presence of fast ions in electrostatic turbulence.
△ Less
Submitted 8 January, 2017; v1 submitted 25 August, 2016;
originally announced August 2016.
-
Plasma rotation from momentum transport by neutrals in tokamaks
Authors:
John Omotani,
István Pusztai,
Tünde Fülöp
Abstract:
Neutral atoms can strongly influence the intrinsic rotation and radial electric field at the tokamak edge. Here, we present a framework to investigate these effects when the neutrals dominate the momentum transport. We explore the parameter space numerically, using highly flexible model geometries and a state of the art kinetic solver. We find that the most important parameters controlling the tor…
▽ More
Neutral atoms can strongly influence the intrinsic rotation and radial electric field at the tokamak edge. Here, we present a framework to investigate these effects when the neutrals dominate the momentum transport. We explore the parameter space numerically, using highly flexible model geometries and a state of the art kinetic solver. We find that the most important parameters controlling the toroidal rotation and electric field are the major radius where the neutrals are localized and the plasma collisionality. This offers a means to influence the rotation and electric field by, for example, varying the radial position of the X-point to change the major radius of the neutral peak.
△ Less
Submitted 5 May, 2016; v1 submitted 27 April, 2016;
originally announced April 2016.
-
Global effects on neoclassical transport in the pedestal with impurities
Authors:
I. Pusztai,
S. Buller,
M. Landreman
Abstract:
We present a numerical study of collisional transport in a tokamak pedestal in the presence of non-trace impurities, using the radially global $δf$ neoclassical solver PERFECT [M. Landreman et al. 2014 Plasma Phys. Control. Fusion 56 045005]. It is known that in a tokamak core with non-trace impurities present the radial impurity flux opposes the bulk ion flux to provide an ambipolar particle tran…
▽ More
We present a numerical study of collisional transport in a tokamak pedestal in the presence of non-trace impurities, using the radially global $δf$ neoclassical solver PERFECT [M. Landreman et al. 2014 Plasma Phys. Control. Fusion 56 045005]. It is known that in a tokamak core with non-trace impurities present the radial impurity flux opposes the bulk ion flux to provide an ambipolar particle transport, with the electron transport being negligibly small. However, in a sharp density pedestal with sub-sonic ion flows the electron transport can be comparable to the ion and impurity flows. Furthermore, the neoclassical particle transport is not intrinsically ambipolar, and the non-ambipolarity of the fluxes extends outside the pedestal region by the radial coupling of the perturbations. The neoclassical momentum transport, which is finite in the presence of ion orbit-width scale profile variations, is significantly enhanced when impurities are present in non-trace quantities, even if the total parallel mass flow is dominated by the bulk ions.
△ Less
Submitted 30 June, 2016; v1 submitted 2 March, 2016;
originally announced March 2016.
-
Turbulent transport of MeV range cyclotron heated minorities as compared to alpha particles
Authors:
István Pusztai,
George J. Wilkie,
Yevgen O. Kazakov,
Tünde Fülöp
Abstract:
We study the turbulent transport of an ion cyclotron resonance heated (ICRH), MeV range minority ion species in tokamak plasmas. Such highly energetic minorities, which can be produced in the three ion minority heating scheme [Ye. O. Kazakov et al. (2015) Nucl. Fusion 55, 032001], have been proposed to be used to experimentally study the confinement properties of fast ions without the generation o…
▽ More
We study the turbulent transport of an ion cyclotron resonance heated (ICRH), MeV range minority ion species in tokamak plasmas. Such highly energetic minorities, which can be produced in the three ion minority heating scheme [Ye. O. Kazakov et al. (2015) Nucl. Fusion 55, 032001], have been proposed to be used to experimentally study the confinement properties of fast ions without the generation of fusion alphas. We compare the turbulent transport properties of ICRH ions with that of fusion born alpha particles. Our theoretical predictions indicate that care must be taken when conclusions are drawn from experimental results: While the effect of turbulence on these particles is similar in terms of transport coefficients, differences in their distribution functions - ultimately their generation processes - make the resulting turbulent fluxes different.
△ Less
Submitted 30 June, 2016; v1 submitted 22 February, 2016;
originally announced February 2016.
-
Bump formation in the runaway electron tail
Authors:
J. Decker,
E. Hirvijoki,
O. Embreus,
Y. Peysson,
A. Stahl,
I. Pusztai,
T. Fülöp
Abstract:
Runaway electrons are generated in a magnetized plasma when the parallel electric field exceeds a critical value. For such electrons with energies typically reaching tens of MeV, the Abraham-Lorentz-Dirac (ALD) radiation force, in reaction to the synchrotron emission, is significant and can be the dominant process limiting the electron acceleration. The effect of the ALD-force on runaway electron…
▽ More
Runaway electrons are generated in a magnetized plasma when the parallel electric field exceeds a critical value. For such electrons with energies typically reaching tens of MeV, the Abraham-Lorentz-Dirac (ALD) radiation force, in reaction to the synchrotron emission, is significant and can be the dominant process limiting the electron acceleration. The effect of the ALD-force on runaway electron dynamics in a homogeneous plasma is investigated using the relativistic finite-difference Fokker-Planck codes LUKE [Decker & Peysson, Report EUR-CEA-FC-1736, Euratom-CEA, (2004)] and CODE [Landreman et al, Comp. Phys. Comm. 185, 847 (2014)]. Under the action of the ALD force, we find that a bump is formed in the tail of the electron distribution function if the electric field is sufficiently large. We also observe that the energy of runaway electrons in the bump increases with the electric field amplitude, while the population increases with the bulk electron temperature. The presence of the bump divides the electron distribution into a runaway beam and a bulk population. This mechanism may give rise to beam-plasma types of instabilities that could in turn pump energy from runaway electrons and alter their confinement.
△ Less
Submitted 16 March, 2015; v1 submitted 12 March, 2015;
originally announced March 2015.
-
Radiation reaction induced non-monotonic features in runaway electron distributions
Authors:
E. Hirvijoki,
I. Pusztai,
J. Decker,
O. Embréus,
A. Stahl,
T. Fülöp
Abstract:
Runaway electrons, which are generated in a plasma where the induced electric field exceeds a certain critical value, can reach very high energies in the MeV range. For such energetic electrons, radiative losses will contribute significantly to the momentum space dynamics. Under certain conditions, due to radiative momentum losses, a non-monotonic feature - a "bump" - can form in the runaway elect…
▽ More
Runaway electrons, which are generated in a plasma where the induced electric field exceeds a certain critical value, can reach very high energies in the MeV range. For such energetic electrons, radiative losses will contribute significantly to the momentum space dynamics. Under certain conditions, due to radiative momentum losses, a non-monotonic feature - a "bump" - can form in the runaway electron tail, creating a potential for bump-on-tail-type instabilities to arise. Here we study the conditions for the existence of the bump. We derive an analytical threshold condition for bump appearance and give an approximate expression for the minimum energy at which the bump can appear. Numerical calculations are performed to support the analytical derivations.
△ Less
Submitted 23 March, 2015; v1 submitted 11 February, 2015;
originally announced February 2015.
-
Core micro-instability analysis of JET hybrid and baseline discharges with carbon wall
Authors:
S. Moradi,
I. Pusztai,
I. Voitsekhovitch,
L. Garzotti,
C. Bourdelle,
M. J. Pueschel,
I. Lupelli,
M. Romanelli,
the JET-EFDA contributors
Abstract:
The core micro-instability characteristics of hybrid and baseline plasmas in a selected set of JET plasmas with carbon wall are investigated through local linear and non-linear and global linear gyro-kinetic simulations with the GYRO code [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. In particular, we study the role of plasma pressure on the micro-instabilities, and scan the pa…
▽ More
The core micro-instability characteristics of hybrid and baseline plasmas in a selected set of JET plasmas with carbon wall are investigated through local linear and non-linear and global linear gyro-kinetic simulations with the GYRO code [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. In particular, we study the role of plasma pressure on the micro-instabilities, and scan the parameter space for the important plasma parameters responsible for the onset and stabilization of the modes under experimental conditions. We find that a good core confinement due to strong stabilization of the micro-turbulence driven transport can be expected in the hybrid plasmas due to the stabilizing effect of the fast ion pressure that is more effective at the low magnetic shear of the hybrid discharges. While parallel velocity gradient destabilization is important for the inner core, at outer radii the hybrid plasmas may benefit from a strong quench of the turbulence transport by $\mathbf{E}\times\mathbf{B}$ rotation shear.
△ Less
Submitted 18 November, 2014;
originally announced November 2014.
-
High-m Kink/Tearing Modes in Cylindrical Geometry
Authors:
J. W. Connor,
R. J. Hastie,
I. Pusztai,
P. J. Catto,
M. Barnes
Abstract:
The global ideal kink equation, for cylindrical geometry and zero beta, is simplified in the high poloidal mode number limit and used to determine the tearing stability parameter, $Δ^\prime$. In the presence of a steep monotonic current gradient, $Δ^\prime$ becomes a function of a parameter, $σ_0$, characterising the ratio of the maximum current gradient to magnetic shear, and $x_s$, characterisin…
▽ More
The global ideal kink equation, for cylindrical geometry and zero beta, is simplified in the high poloidal mode number limit and used to determine the tearing stability parameter, $Δ^\prime$. In the presence of a steep monotonic current gradient, $Δ^\prime$ becomes a function of a parameter, $σ_0$, characterising the ratio of the maximum current gradient to magnetic shear, and $x_s$, characterising the separation of the resonant surface from the maximum of the current gradient. In equilibria containing a current "spike", so that there is a non-monotonic current profile, $Δ^\prime$ also depends on two parameters: $κ$, related to the ratio of the curvature of the current density at its maximum to the magnetic shear, and $x_s$, which now represents the separation of the resonance from the point of maximum current density. The relation of our results to earlier studies of tearing modes and to recent gyro-kinetic calculations of current driven instabilities, is discussed, together with potential implications for the stability of the tokamak pedestal.
△ Less
Submitted 26 September, 2014;
originally announced September 2014.
-
Impurity transport in Alcator C-Mod in the presence of poloidal density variation induced by ion cyclotron resonance heating
Authors:
Albert Mollén,
István Pusztai,
Matthew L. Reinke,
Yevgen O. Kazakov,
Nathan T. Howard,
Emily A. Belli,
Tünde Fülöp,
the Alcator C-Mod Team
Abstract:
Impurity particle transport in an ion cyclotron resonance heated Alcator C-Mod discharge is studied with local gyrokinetic simulations and a theoretical model including the effect of poloidal asymmetries and elongation. In spite of the strong minority temperature anisotropy in the deep core region, the poloidal asymmetries are found to have a negligible effect on the turbulent impurity transport d…
▽ More
Impurity particle transport in an ion cyclotron resonance heated Alcator C-Mod discharge is studied with local gyrokinetic simulations and a theoretical model including the effect of poloidal asymmetries and elongation. In spite of the strong minority temperature anisotropy in the deep core region, the poloidal asymmetries are found to have a negligible effect on the turbulent impurity transport due to low magnetic shear in this region, in agreement with the experimental observations. According to the theoretical model, in outer core regions poloidal asymmetries may contribute to the reduction of the impurity peaking, but uncertainties in atomic physics processes prevent quantitative comparison with experiments.
△ Less
Submitted 10 November, 2014; v1 submitted 3 February, 2014;
originally announced February 2014.
-
Radially global $δf$ computation of neoclassical phenomena in a tokamak pedestal
Authors:
Matt Landreman,
Felix I Parra,
Peter J Catto,
Darin R Ernst,
Istvan Pusztai
Abstract:
Conventional radially-local neoclassical calculations become inadequate if the radial gradient scale lengths of the H-mode pedestal become as small as the poloidal ion gyroradius. Here, we describe a radially global $δf$ continuum code that generalizes neoclassical calculations to allow stronger gradients. As with conventional neoclassical calculations, the formulation is time-independent and requ…
▽ More
Conventional radially-local neoclassical calculations become inadequate if the radial gradient scale lengths of the H-mode pedestal become as small as the poloidal ion gyroradius. Here, we describe a radially global $δf$ continuum code that generalizes neoclassical calculations to allow stronger gradients. As with conventional neoclassical calculations, the formulation is time-independent and requires only the solution of a single sparse linear system. We demonstrate precise agreement with an asymptotic analytic solution of the radially global kinetic equation in the appropriate limits of aspect ratio and collisionality. This agreement depends crucially on accurate treatment of finite orbit width effects.
△ Less
Submitted 5 February, 2014; v1 submitted 7 December, 2013;
originally announced December 2013.
-
Radio frequency induced and neoclassical asymmetries and their effects on turbulent impurity transport in a tokamak
Authors:
István Pusztai,
Matt Landreman,
Albert Mollén,
Yevgen O. Kazakov,
Tünde Fülöp
Abstract:
Poloidal asymmetries in the impurity density can be generated by radio frequency heating in the core and by neoclassical effects in the edge of tokamak plasmas. In a pedestal case study, using global neoclassical simulations we find that finite orbit width effects can generate significant poloidal variation in the electrostatic potential, which varies on a small radial scale. Gyrokinetic modeling…
▽ More
Poloidal asymmetries in the impurity density can be generated by radio frequency heating in the core and by neoclassical effects in the edge of tokamak plasmas. In a pedestal case study, using global neoclassical simulations we find that finite orbit width effects can generate significant poloidal variation in the electrostatic potential, which varies on a small radial scale. Gyrokinetic modeling shows that these poloidal asymmetries can be strong enough to significantly modify turbulent impurity peaking. In the pedestal the ExB drift in the radial electric field can give a larger contribution to the poloidal motion of impurities than that of their parallel streaming. Under such circumstances we find that up-down asymmetries can also affect impurity peaking.
△ Less
Submitted 20 November, 2013; v1 submitted 16 September, 2013;
originally announced September 2013.
-
A current driven electromagnetic mode in sheared and toroidal configurations
Authors:
István Pusztai,
Peter J. Catto,
Felix I. Parra,
Michael Barnes
Abstract:
The induced electric field in a tokamak drives a parallel electron current flow. In an inhomogeneous, finite beta plasma, when this electron flow is comparable to the ion thermal speed, the Alfven mode wave solutions of the electromagnetic gyrokinetic equation can become nearly purely growing kink modes. Using the new "low-flow" version of the gyrokinetic code GS2 developed for momentum transport…
▽ More
The induced electric field in a tokamak drives a parallel electron current flow. In an inhomogeneous, finite beta plasma, when this electron flow is comparable to the ion thermal speed, the Alfven mode wave solutions of the electromagnetic gyrokinetic equation can become nearly purely growing kink modes. Using the new "low-flow" version of the gyrokinetic code GS2 developed for momentum transport studies [Barnes et al 2013 Phys. Rev. Lett. 111, 055005], we are able to model the effect of the induced parallel electric field on the electron distribution to study the destabilizing influence of current on stability. We identify high mode number kink modes in GS2 simulations and make comparisons to analytical theory in sheared magnetic geometry. We demonstrate reassuring agreement with analytical results both in terms of parametric dependences of mode frequencies and growth rates, and regarding the radial mode structure.
△ Less
Submitted 6 January, 2014; v1 submitted 4 September, 2013;
originally announced September 2013.
-
Three-dimensional modeling of beam emission spectroscopy measurements in fusion plasmas
Authors:
D. Guszejnov,
G. I. Pokol,
I. Pusztai,
D. Refy,
S. Zoletnik,
M. Lampert,
Y. U. Nam
Abstract:
One of the main diagnostic tools for measuring electron density profiles and the characteristics of long wavelength turbulent wave structures in fusion plasmas is Beam Emission Spectroscopy (BES). The increasing number of BES systems necessitated an accurate and comprehensive simulation of BES diagnostics, which in turn motivated the development of the RENATE simulation code that is the topic of t…
▽ More
One of the main diagnostic tools for measuring electron density profiles and the characteristics of long wavelength turbulent wave structures in fusion plasmas is Beam Emission Spectroscopy (BES). The increasing number of BES systems necessitated an accurate and comprehensive simulation of BES diagnostics, which in turn motivated the development of the RENATE simulation code that is the topic of this paper. RENATE is a modular, fully three-dimensional code incorporating all key features of BES systems from the atomic physics to the observation, including an advanced modeling of the optics. Thus RENATE can be used both in the interpretation of measured signals and the development of new BES systems. The most important components of the code have been successfully benchmarked against other simulation codes. The primary results have been validated against experimental data from the KSTAR tokamak.
△ Less
Submitted 10 January, 2013;
originally announced January 2013.
-
Impurity transport in trapped electron mode driven turbulence
Authors:
A. Mollén,
I. Pusztai,
T. Fülöp,
S. Moradi
Abstract:
Trapped electron mode turbulence is studied by gyrokinetic simulations with the GYRO code and an analytical model including the effect of a poloidally varying electrostatic potential. Its impact on radial transport of high-Z trace impurities close to the core is thoroughly investigated and the dependence of the zero-flux impurity density gradient (peaking factor) on local plasma parameters is pres…
▽ More
Trapped electron mode turbulence is studied by gyrokinetic simulations with the GYRO code and an analytical model including the effect of a poloidally varying electrostatic potential. Its impact on radial transport of high-Z trace impurities close to the core is thoroughly investigated and the dependence of the zero-flux impurity density gradient (peaking factor) on local plasma parameters is presented. Parameters such as ion-to-electron temperature ratio, electron temperature gradient and main species density gradient mainly affect the impurity peaking through their impact on mode characteristics. The poloidal asymmetry, the safety factor and magnetic shear have the strongest effect on impurity peaking, and it is shown that under certain scenarios where trapped electron modes are dominant, core accumulation of high-Z impurities can be avoided. We demonstrate that accounting for the momentum conservation property of the impurity-impurity collision operator can be important for an accurate evaluation of the impurity peaking factor.
△ Less
Submitted 3 April, 2013; v1 submitted 8 January, 2013;
originally announced January 2013.
-
Micro-tearing modes in spherical and conventional tokamaks
Authors:
S. Moradi,
I. Pusztai,
W. Guttenfelder,
T. Fülöp,
A. Mollén
Abstract:
The onset and characteristics of Micro-Tearing Modes (MTM) in the core of spherical (NSTX) and conventional tokamaks (ASDEX-UG and JET) are studied through local linear gyrokinetic simulations with gyro [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. For experimentally relevant core plasma parameters in the NSTX and ASDEX-UG tokamaks, in agreement with previous works, we find MTM…
▽ More
The onset and characteristics of Micro-Tearing Modes (MTM) in the core of spherical (NSTX) and conventional tokamaks (ASDEX-UG and JET) are studied through local linear gyrokinetic simulations with gyro [J. Candy and E. Belli, General Atomics Report GA-A26818 (2011)]. For experimentally relevant core plasma parameters in the NSTX and ASDEX-UG tokamaks, in agreement with previous works, we find MTMs as the dominant linear instability. Also, for JET-like core parameters considered in our study an MTM is found as the most unstable mode. In all these plasmas, finite collisionality is needed for MTMs to become unstable and the electron temperature gradient is found to be the fundamental drive. However, a significant difference is observed in the dependence of linear growth rate of MTMs on electron temperature gradient. While it varies weakly and non-monotonically in JET and ASDEX-UG plasmas, in NSTX it increases with the electron temperature gradient.
△ Less
Submitted 7 January, 2013;
originally announced January 2013.
-
Turbulent transport of impurities and their effect on energy confinement
Authors:
I. Pusztai,
A. Mollen,
T. Fulop,
J. Candy
Abstract:
By presenting linear and nonlinear gyrokinetic studies, based on a balanced neutral beam injection deuterium discharge from the DIII-D tokamak, we demonstrate that impurities alter the scaling of the transport on the charge and mass of the main species, and even more importantly, they can dramatically change the energy transport even in relatively small quantities. A poloidally varying equilibrium…
▽ More
By presenting linear and nonlinear gyrokinetic studies, based on a balanced neutral beam injection deuterium discharge from the DIII-D tokamak, we demonstrate that impurities alter the scaling of the transport on the charge and mass of the main species, and even more importantly, they can dramatically change the energy transport even in relatively small quantities. A poloidally varying equilibrium electrostatic potential can lead to a strong reduction or sign change of the impurity peaking factor due to the combined effect of the in-out impurity density asymmetry and the EXB drift of impurities. We present an approximate expression for the impurity peaking factor and demonstrate that impurity peaking is not significantly affected by impurity self-collisions.
△ Less
Submitted 29 November, 2012; v1 submitted 10 October, 2012;
originally announced October 2012.
-
Possible mechanism responsible for observed impurity outward flow under radio frequency heating
Authors:
S. Moradi,
T. Fülöp,
A. Mollén,
I. Pusztai
Abstract:
The effect of poloidal asymmetry of impurities on impurity transport driven by electrostatic turbulence in tokamak plasmas is analyzed. It is found that in the presence of in-out asymmetric impurity populations the zero-flux impurity density gradient (the so-called peaking factor) is significantly reduced. A sign change in the impurity flux may occur if the asymmetry is sufficiently large. This ma…
▽ More
The effect of poloidal asymmetry of impurities on impurity transport driven by electrostatic turbulence in tokamak plasmas is analyzed. It is found that in the presence of in-out asymmetric impurity populations the zero-flux impurity density gradient (the so-called peaking factor) is significantly reduced. A sign change in the impurity flux may occur if the asymmetry is sufficiently large. This may be a contributing reason for the observed outward convection of impurities in the presence of radio frequency heating. The effect of in-out asymmetry is most significant in regions with low temperature gradients. In the trapped electron mode dominated case also an up-down asymmetry can lead to a sign change in the peaking factor from positive to negative. The effect of ion parallel compressibility on the peaking factor is significant, and leads to positive peaking factors in regions with high temperature gradients, even in the presence of in-out asymmetry.
△ Less
Submitted 1 June, 2011;
originally announced June 2011.
-
Turbulence characteristics of electron cyclotron and ohmic heated discharges
Authors:
Istvan Pusztai,
Sara Moradi,
Tunde Fulop,
Natalia Timchenko
Abstract:
Turbulence characteristics of electron cyclotron (EC) and ohmic heated (OH) discharges has been analyzed by electrostatic gyrokinetic simulations with GYRO [J. Candy, R.E. Waltz, Journal of Computational Physics 186, 545-581 (2003)] aiming to find insights into the effect of auxiliary heating on the transport. Trapped electron modes are found to be unstable in both OH and the EC heated scenarios.…
▽ More
Turbulence characteristics of electron cyclotron (EC) and ohmic heated (OH) discharges has been analyzed by electrostatic gyrokinetic simulations with GYRO [J. Candy, R.E. Waltz, Journal of Computational Physics 186, 545-581 (2003)] aiming to find insights into the effect of auxiliary heating on the transport. Trapped electron modes are found to be unstable in both OH and the EC heated scenarios. In the OH case the main drive is from the density gradient and in the EC case from the electron temperature gradient. The growth rates and particle fluxes exhibit qualitatively different scaling with the electron-to-ion temperature ratios in the two cases. This is mainly due to the fact that the dominant drives and the collisionalities are different. The inward flow velocity of impurities and the impurity diffusion coefficient decreases when applying EC heating, which leads to lower impurity peaking, consistently with experimental observations.
△ Less
Submitted 6 April, 2011; v1 submitted 5 April, 2011;
originally announced April 2011.