Source author record

I. Pusztai

I. Pusztai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph

Catalog footprint

What is connected

11works

1topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Runaway dynamics in reactor-scale spherical tokamak disruptions

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.

preprint2022arXiv

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

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.

preprint2021arXiv

Spatiotemporal analysis of the runaway distribution function from synchrotron images in an ASDEX Upgrade disruption

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.

preprint2016arXiv

Global effects on neoclassical transport in the pedestal with impurities

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.

preprint2015arXiv

Bump formation in the runaway electron tail

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.

preprint2014arXiv

Core micro-instability analysis of JET hybrid and baseline discharges with carbon wall

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.

preprint2014arXiv

High-m Kink/Tearing Modes in Cylindrical Geometry

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.

preprint2013arXiv

Micro-tearing modes in spherical and conventional tokamaks

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.

preprint2013arXiv

Three-dimensional modeling of beam emission spectroscopy measurements in fusion plasmas

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.

preprint2012arXiv

Turbulent transport of impurities and their effect on energy confinement

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.

preprint2011arXiv

Possible mechanism responsible for observed impurity outward flow under radio frequency heating

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.