Source author record

T. Fülöp

T. Fülöp appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

26 published item(s)

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 generation during tokamak start-up

Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plasma. This is of particular concern in ITER where engineering limitations put restrictions on the allowable electric fields and limit the prefill densities during start-up. In this study, we present a new 0D burn-through simulation tool called STREAM (STart-up Runaway Electron Analysis Model), which self-consistently evolves the plasma density, temperature and electric field, while accounting for the generation and loss of relativistic runaway electrons. After verifying the burn-through model, we investigate conditions under which runaway electrons can form during tokamak start-up as well as their effects on the plasma initiation. We find that Dreicer generation plays a crucial role in determining whether a discharge becomes runaway-dominated or not, and that a large number of runaway electrons could limit the ohmic heating of the plasma, thus preventing successful burn-through or further ramp-up of the plasma current. The runaway generation can be suppressed by raising the density via gas fuelling, but only if done sufficiently early. Otherwise a large runaway seed may have already been built up, which can avalanche even at relatively low electric fields and high densities.

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.

preprint2020arXiv

Effect of plasma elongation on current dynamics during tokamak disruptions

Plasma terminating disruptions in tokamaks may result in relativistic runaway electron beams with potentially serious consequences for future devices with large plasma currents. In this paper we investigate the effect of plasma elongation on the coupled dynamics of runaway generation and resistive diffusion of the electric field. We find that elongated plasmas are less likely to produce large runaway currents, partly due to the lower induced electric fields associated with larger plasmas, and partly due to direct shaping effects, which mainly lead to a reduction in the runaway avalanche gain.

preprint2020arXiv

Kinetic modelling of runaway electron generation in argon-induced disruptions in ASDEX Upgrade

Massive material injection has been proposed as a way to mitigate the formation of a beam of relativistic runaway electrons that may result from a disruption in tokamak plasmas. In this paper we analyse runaway generation observed in eleven ASDEX Upgrade discharges where disruption was triggered using massive gas injection. We present numerical simulations in scenarios characteristic of on-axis plasma conditions, constrained by experimental observations, using a description of the runaway dynamics with self-consistent electric field and temperature evolution in two-dimensional momentum space and zero-dimensional real space. We describe the evolution of the electron distribution function during the disruption, and show that the runaway seed generation is dominated by hot-tail in all of the simulated discharges. We reproduce the observed dependence of the current dissipation rate on the amount of injected argon during the runaway plateau phase. Our simulations also indicate that above a threshold amount of injected argon, the current density after the current quench depends strongly on the argon densities. This trend is not observed in the experiments, which suggests that effects not captured by 0D kinetic modeling -- such as runaway seed transport -- are also important.

preprint2020arXiv

Runaway electron synchrotron radiation in a vertically translated plasma

Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron diagnostic Soft. This experimental validation lends confidence to the theory underlying the synthetic diagnostics which are used for benchmarking theoretical models of and probing runaway dynamics. We present a comparison of synchrotron measurements in TCV with predictions of kinetic theory for runaway dynamics in uniform magnetic fields. We find that to explain the detected synchrotron emission, significant non-collisional pitch angle scattering as well as radial transport of REs would be needed. Such effects could be caused by the presence of magnetic perturbations, which should be further investigated in future TCV experiments.

preprint2016arXiv

Kinetic modelling of runaway electrons in dynamic scenarios

Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concern. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.

preprint2016arXiv

NORSE: A solver for the relativistic non-linear Fokker-Planck equation for electrons in a homogeneous plasma

Energetic electrons are of interest in many types of plasmas, however previous modeling of their properties has been restricted to the use of linear Fokker-Planck collision operators or non-relativistic formulations. Here, we describe a fully non-linear kinetic-equation solver, capable of handling large electric-field strengths (compared to the Dreicer field) and relativistic temperatures. This tool allows modeling of the momentum-space dynamics of the electrons in cases where strong departures from Maxwellian distributions may arise. As an example, we consider electron runaway in magnetic-confinement fusion plasmas and describe a transition to electron slide-away at field strengths significantly lower than previously predicted.

preprint2016arXiv

Runaway-electron formation and electron slide-away in an ITER post-disruption scenario

Mitigation of runaway electrons is one of the outstanding issues for the reliable operation of ITER and other large tokamaks, and accurate estimates for the expected runaway-electron energies and current are needed. Previously, linearized tools (which assume the runaway population to be small) have been used to study the runaway dynamics, but these tools are not valid in the cases of most interest, i.e. when the runaway population becomes substantial. We study runaway-electron formation in a post-disruption ITER plasma using the newly developed non-linear code NORSE, and describe a feedback mechanism by which a transition to electron slide-away can be induced at field strengths significantly lower than previously expected. If the electric field is actively imposed using the control system, the entire electron population is quickly converted to runaways in the scenario considered. We find the time until the feedback mechanism sets in to be highly dependent on the details of the mechanisms removing heat from the thermal electron population.

preprint2016arXiv

Vlasov modelling of laser-driven collisionless shock acceleration of protons

Ion acceleration due to the interaction between a short high-intensity laser pulse and a moderately overdense plasma target is studied using Eulerian Vlasov-Maxwell simulations. The effects of variations in the plasma density profile and laser pulse parameters are investigated, and the interplay of collisionless shock and target normal sheath acceleration is analyzed. It is shown that the use of a layered-target with a combination of light and heavy ions, on the front and rear side respectively, yields a strong quasi-static sheath-field on the rear side of the heavy-ion part of the target. This sheath-field increases the energy of the shock-accelerated ions while preserving their mono-energeticity.

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

Distinguished rheological models in the framework of a thermodynamical internal variable theory

We present and analyze a thermodynamical theory of rheology with single internal variable. The universality of the model is ensured as long as the mesoscopic and/or microscopic background processes satisfy the applied thermodynamical principles, which are the second law, the basic balances and the existence of an additional-tensorial-state variable. The resulting model, which we suggest to call the Kluitenberg-Verhás body, is the Poynting-Thomson-Zener body with an additional inertial element, or, in other words, is the extension of Jeffreys model to solids. We argue that this Kluitenberg-Verhás body is the natural thermodynamical building block of rheology. An important feature of the presented methodology is that nontrivial inequality-type restrictions arise for the four parameters of the model. We compare these conditions and other aspects to those of other known thermodynamical approaches, like Extended Irreversible Thermodynamics or the original theory of Kluitenberg.

preprint2014arXiv

Quasi-linear analysis of the extraordinary electron wave destabilized by runaway electrons

Runaway electrons with strongly anisotropic distributions present in post-disruption tokamak plasmas can destabilize the extraordinary electron (EXEL) wave. The present work investigates the dynamics of the quasi-linear evolution of the EXEL instability for a range of different plasma parameters using a model runaway distribution function valid for highly relativistic runaway electron beams produced primarily by the avalanche process. Simulations show a rapid pitch-angle scattering of the runaway electrons in the high energy tail on the $100-1000\;\rm μs$ time scale. Due to the wave-particle interaction, a modification to the synchrotron radiation spectrum emitted by the runaway electron population is foreseen, exposing a possible experimental detection method for such an interaction.

preprint2014arXiv

Switching Reciprocity On and Off in a Magneto-Optical X-Ray Scattering Experiment Using Nuclear Resonance of 57Fe Foils

Reciprocity is when the scattering amplitude of wave propagation satisfies a symmetry property, connecting a scattering process with an appropriate reversed one. We report on an experiment using nuclear resonance scattering of synchrotron radiation, which demonstrates that magneto-optical materials do not necessarily violate reciprocity. The setting enables to switch easily between reciprocity and its violation. In the latter case, the exhibited reciprocity violation is orders of magnitude larger than achieved by previous wave scattering experiments.

preprint2014arXiv

Thermodynamic hierarchies of evolution equations

Non-equilibrium thermodynamics with internal variables introduces a natural hierarchical arrangement of evolution equations. Three examples are shown: a hierarchy of linear constitutive equations in thermodynamic rhelogy with a single internal variable, a hierarchy of wave equations in the theory of generalized continua with dual internal variables and a hierarchical arrangement of the Fourier equation in the theory of heat conduction with current multipliers.

preprint2013arXiv

Effect of impurities on the transition between minority ion and mode conversion ICRH heating in (3He)-H tokamak plasmas

Hydrogen majority plasmas will be used in the initial non-activated phase of ITER operation. Optimizing ion cyclotron resonance heating (ICRH) in such scenarios will help in achieving H-mode in these plasmas. Past JET experiments with the carbon wall revealed a significant impact of intrinsic impurities on the ICRH performance in (3He)-H plasmas relevant for the full-field initial ITER phase. High plasma contamination with carbon impurities resulted in the appearance of a supplementary mode conversion layer and significant reduction in the transition concentration of 3He minority ions, defined as the concentration at which the change from minority heating to mode conversion regime occurs. In view of the installation of the new ITER-like wall at JET, it is important to evaluate the effect of Be and W impurities on ICRH scenarios in (3He)-H plasmas. In this paper, an approximate analytical expression for the transition concentration of 3He minority ions is derived as a function of plasma and ICRH parameters, and accounting for typical impurity species at JET. The accompanying 1D wave modeling supports the analytical results and suggests a potential experimental method to reduce 3He level needed to achieve a specific heating regime by puffing a small amount of 4He ions additionally to (3He)-H plasma.

preprint2013arXiv

Electromagnetic waves destabilized by runaway electrons in near-critical electric fields

Runaway electron distributions are strongly anisotropic in velocity space. This anisotropy is a source of free energy that may destabilize electromagnetic waves through a resonant interaction between the waves and the energetic electrons. In this work we investigate the high-frequency electromagnetic waves that are destabilized by runaway electron beams when the electric field is close to the critical field for runaway acceleration. Using a runaway electron distribution appropriate for the near-critical case we calculate the linear instability growth rate of these waves and conclude that the obliquely propagating whistler waves are most unstable. We show that the frequencies, wave numbers and propagation angles of the most unstable waves depend strongly on the magnetic field. Taking into account collisional and convective damping of the waves, we determine the number density of runaways that is required to destabilize the waves and show its parametric dependences.

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

Synchrotron radiation from a runaway electron distribution in tokamaks

The synchrotron radiation emitted by runaway electrons in a fusion plasma provides information regarding the particle momenta and pitch-angles of the runaway electron population through the strong dependence of the synchrotron spectrum on these parameters. Information about the runaway density and its spatial distribution, as well as the time evolution of the above quantities, can also be deduced. In this paper we present the synchrotron radiation spectra for typical avalanching runaway electron distributions. Spectra obtained for a distribution of electrons are compared to the emission of mono-energetic electrons with a prescribed pitch-angle. We also examine the effects of magnetic field curvature and analyse the sensitivity of the resulting spectrum to perturbations to the runaway distribution. The implications for the deduced runaway electron parameters are discussed. We compare our calculations to experimental data from DIII-D and estimate the maximum observed runaway energy.

preprint2013arXiv

The effect of ITER-like wall on runaway electron generation in JET

This paper investigates the effect of the ITER-like wall (ILW) on runaway electron (RE) generation through a comparative study of similar slow argon injection JET disruptions, performed with different wall materials. In the carbon wall case, a runaway electron plateau is observed, while in the ITER-like wall case, the current quench is slower and the runaway current is negligibly small. The aim of the paper is to shed light on the reason for these differences by detailed numerical modelling to study which factors affected the RE formation. The post-disruption current profile is calculated by a one-dimensional model of electric field, temperature and runaway current taking into account the impurity injection. Scans of various impurity contents are performed and agreement with the experimental scenarios is obtained for reasonable argon- and wall impurity contents. Our modelling shows that the reason for the changed RE dynamics is a complex, combined effect of the differences in plasma parameter profiles, the radiation characteristics of beryllium and carbon, and the difference of the injected argon amount. These together lead to a significantly higher Dreicer generation rate in the carbon wall case, which is less prone to be suppressed by RE loss mechanisms. The results indicate that the differences are greatly reduced above ~50% argon content, suggesting that significant RE current is expected in future massive gas injection experiments on both JET and ITER.

preprint2012arXiv

Kinematic quantities of finite elastic and plastic deformation

Kinematic quantities for finite elastic and plastic deformations are defined via an approach that does not rely on auxiliary elements like reference frame and reference configuration, and that gives account of the inertial-noninertial aspects explicitly. These features are achieved by working on Galilean spacetime directly. The quantity expressing elastic deformations is introduced according to its expected role: to measure how different the current metric is from the relaxed/stressless metric. Further, the plastic kinematic quantity is the change rate of the stressless metric. The properties of both are analyzed, and their relationship to frequently used elastic and plastic kinematic quantities is discussed. One important result is that no objective elastic or plastic quantities can be defined from deformation gradient.

preprint2012arXiv

Universal heat conduction -- the thermodynamics of some weakly nonlocal theories

A linear irreversible thermodynamic framework of heat conduction in rigid conductors is introduced. The deviation from local equilibrium is characterized by a single internal variable and a current intensity factor. A general constitutive evolution equation of the current density of the internal energy is derived by introducing linear relationship between the thermodynamic forces and fluxes. The Fourier, Maxwell-Cattaneo-Vernotte, Guyer-Krumhansl, Jeffreys type and Green-Naghdi type equations of heat conduction are obtained as special cases. The universal character of the approach is demonstrated by two examples. Solutions illustrating the properties of the equation with jump boundary conditions are given.

preprint2011arXiv

Impurity flows and plateau-regime poloidal density variation in a tokamak pedestal

In the pedestal of a tokamak, the sharp radial gradients of density and temperature can give rise to poloidal variation in the density of impurities. At the same time, the flow of the impurity species is modified relative to the conventional neoclassical result. In this paper, these changes to the density and flow of a collisional impurity species are calculated for the case when the main ions are in the plateau regime. In this regime it is found that the impurity density can be higher at either the inboard or outboard side. This finding differs from earlier results for banana- or Pfirsch-Schlüter-regime main ions, in which case the impurity density is always higher at the inboard side in the absence of rotation. Finally, the modifications to the impurity flow are also given for the other regimes of main-ion collisionality.

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.