Source author record

the ASDEX Upgrade team

the ASDEX Upgrade team 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

22works
6topics
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

22 published item(s)

preprint2022arXiv

Implementation of synthetic fast-ion loss detector and imaging heavy ion beam probe diagnostics in the 3D hybrid kinetic-MHD code MEGA

A synthetic Fast-Ion Loss Detector (FILD) and an imaging Heavy Ion Beam Probe (i-HIBP) have been implemented in the 3D hybrid kinetic-magnetohydrodynamic code MEGA. First synthetic measurements from these two diagnostics have been obtained for neutral beam injection (NBI) driven Alfvén Eigenmode (AE) simulated with MEGA. The synthetic fast-ion losses show a strong correlation with the AE amplitude. This correlation is observed in the phase-space, represented in coordinates toroidal canonical momentum and energy. Fast-ion losses and the energy exchange diagrams of the confined population are connected with lines of constant E' , a linear combination of E and Pϕ . First i-HIBP synthetic signals also have been computed for the simulated AE, showing displacements in the strikeline of the order of around 1 mm, above the expected resolution in the i-HIBP scintillator of approximately 100 μm.

preprint2021arXiv

I-mode pedestal relaxation events in the Alcator C-Mod and ASDEX Upgrade tokamaks

In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light into their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks $-$ Alcator C-Mod and ASDEX Upgrade $-$ is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, the nature of the triggering instability is discussed by comparing measurements close to the separatrix in both devices. The PRE relative energy loss from the confined region increases with decreasing pedestal top collisionality $ν_{\mathrm{ped}}^*$. In addition, the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing $ν_{\mathrm{ped}}^*$. Finally, the peak parallel energy fluence due to the PRE measured on the divertor in both devices is compared to the model introduced in [1] for type-I ELMs. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections to future devices such as DEMO and ARC.

preprint2021arXiv

MHD simulations of small ELMs at low triangularity in ASDEX Upgrade

The development of small- and no-ELM regimes for ITER is a high priority topic due to the risks associated to type-I ELMs. By considering non-linear extended MHD simulations of the ASDEX Upgrade tokamak with the JOREK code, we probe a regime that avoids type-I ELMs completely provided that the separatrix density is high enough. The dynamics of the pedestal in this regime are observed to be qualitatively similar to the so-called quasi-continuous exhaust (QCE) regime in several ways. Repetitive type-I ELMs are substituted by roughly constant levels of outwards transport caused by peeling-ballooning modes (with dominant ballooning characteristics) which are localised in the last 5\% of the confined region (in normalised poloidal flux). The simulated low triangularity plasma transitions to a type-I ELMy H-mode if the separatrix density is sufficiently reduced or if the input heating power is sufficiently increased. The stabilising factors that play a role in the suppression of the small ELMs are also investigated by analysing the simulations, and the importance of including diamagnetic effects in the simulations is highlighted. By considering a scan in the pedestal resistivity and by measuring the poloidal velocity of the modes (and comparing to theoretical estimates for ideal and resistive modes), we identify the underlying instabilities as resistive peeling-ballooning modes. Decreasing the resistivity below experimentally-relevant conditions (i.e., going towards ideal MHD), the peeling-ballooning modes that constrain the pedestal below the type-I ELM stability boundary display sharply decreasing growth rates.

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

Flux-driven integrated modelling of main ion pressure and trace tungsten transport in ASDEX Upgrade

Neoclassical and turbulent heavy impurity transport in tokamak core plasmas are determined by main ion temperature, density and toroidal rotation profiles. Thus, in order to understand and prevent experimental behaviour of W accumulation, flux-driven integrated modelling of main ion heat and particle transport over multiple confinement times is a vital prerequisite. For the first time, the quasilinear gyrokinetic code QuaLiKiz has been applied for successful predictions of core kinetic profiles in an ASDEX Upgrade H-mode discharge in the turbulence dominated region within the integrated modelling suite JETTO. Neoclassical contributions are calculated by NCLASS; auxiliary heat and particle deposition profiles due to NBI and ECRH prescribed from previous analysis with TRANSP. Turbulent and neoclassical contributions are insufficient in explaining main ion heat and particle transport inside the $q=1$ surface, necessitating the prescription of further transport coefficients to mimic the impact of MHD activity on central transport. The ion to electron temperature ratio at the simulation boundary at $ρ_\mathrm{tor} = 0.85$ stabilizes ion scale modes while destabilizing ETG modes when significantly exceeding unity. Careful analysis of experimental measurements using Gaussian process regression techniques is carried out to explore reasonable uncertainties. In following trace W impurity transport simulations performed with additionally NEO, neoclassical transport under consideration of poloidal asymmetries alone is found to be insufficient to establish hollow central W density profiles. Reproduction of these conditions measured experimentally is found possible only when assuming the direct impact of a saturated $(m,n)=(1,1)$ MHD mode on heavy impurity transport.

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

On the role of filaments in perpendicular heat transport at the Scrape-off Layer

In this work we carry out quantitative measurements of particle and heat transport associated to SOL filaments in a tokamak, and relate density shoulder formation to the advection of energy in the far SOL. For the first time, this attempt includes direct measurements of ion and electron temperatures for background and filaments. With this aim, we combine data from a number of equivalent L-mode discharges from the ASDEX Upgrade tokamak in which different probe heads were installed on the midplane manipulator. This approach is validated by a comparison with independent diagnostics. Results indicate an increase of heat transport associated to filaments after the shoulder formation. Several centimeters into the SOL, filaments are still found to carry a substantial fraction (up to one fifth) of the power ejected at the separatrix.

preprint2020arXiv

Self-consistent modeling of runaway electron generation in massive gas injection scenarios in ASDEX Upgrade

We present the first successful simulation of a induced disruption in ASDEX Upgrade from massive material injection (MMI) up to established runaway electron (RE) beam, thus covering pre-thermal quench, thermal quench and current quench (CQ) of the discharge. For future high-current fusion devices such as ITER, the successful suppression of REs through MMI is of critical importance to ensure the structural integrity of the vessel. To computationally study the interplay between MMI, background plasma response, and RE generation, a toolkit based on the 1.5D transport code coupling ASTRA-STRAHL is developed. Electron runaway is described by state-of-the-art reduced kinetic models in the presence of partially ionized impurities. Applied to argon MMI in ASDEX Upgrade discharge #33108, key plasma parameters measured experimentally, such as temporal evolution of the line averaged electron density, plasma current decay rate and post-CQ RE current, are well reproduced by the simulation presented. Impurity ions are transported into the central plasma by the combined effect of neoclassical processes and additional effects prescribed inside the $q = 2$ rational surface to explain experimental time scales. Thus, a thermal collapse is induced through strong impurity radiation, giving rise to a substantial RE population as observed experimentally.

preprint2016arXiv

GPEC, a real-time capable Tokamak equilibrium code

A new parallel equilibrium reconstruction code for tokamak plasmas is presented. GPEC allows to compute equilibrium flux distributions sufficiently accurate to derive parameters for plasma control within 1 ms of runtime which enables real-time applications at the ASDEX Upgrade experiment (AUG) and other machines with a control cycle of at least this size. The underlying algorithms are based on the well-established offline-analysis code CLISTE, following the classical concept of iteratively solving the Grad-Shafranov equation and feeding in diagnostic signals from the experiment. The new code adopts a hybrid parallelization scheme for computing the equilibrium flux distribution and extends the fast, shared-memory-parallel Poisson solver which we have described previously by a distributed computation of the individual Poisson problems corresponding to different basis functions. The code is based entirely on open-source software components and runs on standard server hardware and software environments. The real-time capability of GPEC is demonstrated by performing an offline-computation of a sequence of 1000 flux distributions which are taken from one second of operation of a typical AUG discharge and deriving the relevant control parameters with a time resolution of a millisecond. On current server hardware the new code allows employing a grid size of 32x64 zones for the spatial discretization and up to 15 basis functions. It takes into account about 90 diagnostic signals while using up to 4 equilibrium iterations and computing more than 20 plasma-control parameters, including the computationally expensive safety-factor q on at least 4 different levels of the normalized flux.

preprint2016arXiv

Non-linear modeling of the plasma response to RMPs in ASDEX Upgrade

The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces $q = m/n$, the coupling between the m + 2 kink component and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.

preprint2015arXiv

Doppler tomography in fusion plasmas and astrophysics

Doppler tomography is a well-known method in astrophysics to image the accretion flow, often in the shape of thin discs, in compact binary stars. As accretion discs rotate, all emitted line radiation is Doppler-shifted. In fast-ion D-alpha (FIDA) spectroscopy measurements in magnetically confined plasma, the D-alpha-photons are likewise Doppler-shifted ultimately due to gyration of the fast ions. In either case, spectra of Doppler-shifted line emission are sensitive to the velocity distribution of the emitters. Astrophysical Doppler tomography has lead to images of accretion discs of binaries revealing bright spots, spiral structures, and flow patterns. Fusion plasma Doppler tomography has lead to an image of the fast-ion velocity distribution function in the tokamak ASDEX Upgrade. This image matched numerical simulations very well. Here we discuss achievements of the Doppler tomography approach, its promise and limits, analogies and differences in astrophysical and fusion plasma Doppler tomography, and what can be learned by comparison of these applications.

preprint2015arXiv

Effect of sonic poloidal flows in determining flow and density asymmetries for trace impurities in the tokamak edge pedestal

The structure of poloidal and toroidal flows of trace impurities in the edge pedestal of tokamak plasmas is studied analytically and numerically. Parallel momentum balance is analysed upon retaining the following terms: poloidal and toroidal centrifugal forces (inertia), pressure force, electric force, and the friction force. It is shown that, when the poloidal flow is such to produce a properly defined Mach number of order unity somewhere on the flux surface, shock fronts can form. The shock fronts can modify the predicted asymmetry structures in both the flow and the density profile along the poloidal arc. Predictions of the theory are shown against experimental observations in the ASDEX Upgrade tokamak, showing good qualitative and quantitative agreement if the inertia term associated with the poloidal flow is retained.

preprint2015arXiv

Gyrokinetic studies of core turbulence features in ASDEX Upgrade H-mode plasmas

Gyrokinetic validation studies are crucial in developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

preprint2015arXiv

Measurement of a 2D fast-ion velocity distribution function by tomographic inversion of fast-ion D-alpha spectra

We present the first measurement of a local fast-ion 2D velocity distribution function $f(v_\parallel, v_\perp)$. To this end, we heated a plasma in ASDEX Upgrade by neutral beam injection and measured spectra of fast-ion D-alpha (FIDA) light from the plasma center in three views simultaneously. The measured spectra agree very well with synthetic spectra calculated from a TRANSP/NUBEAM simulation. Based on the measured FIDA spectra alone, we infer $f(v_\parallel, v_\perp)$ by tomographic inversion. Salient features of our measurement of $f(v_\parallel, v_\perp)$ agree reasonably well with the simulation: the measured as well as the simulated $f(v_\parallel, v_\perp)$ are lopsided towards negative velocities parallel to the magnetic field, and they have similar shapes. Further, the peaks in the simulation of $f(v_\parallel, v_\perp)$ at full and half injection energies of the neutral beam also appear in the measurement at similar velocity-space locations. We expect that we can measure spectra in up to seven views simultaneously in the next ASDEX Upgrade campaign which would further improve measurements of $f(v_\parallel, v_\perp)$ by tomographic inversion.

preprint2015arXiv

Pellet refuelling of particle loss due to ELM mitigation with RMPs in the ASDEX Upgrade tokamak at low collisionality

The complete refuelling of the plasma density loss (pump-out) caused by mitigation of Edge Localised Modes (ELMs) is demonstrated on the ASDEX Upgrade tokamak. The plasma is refuelled by injection of frozen deuterium pellets and ELMs are mitigated by external resonant magnetic perturbations (RMPs). In this experiment relevant dimensionless parameters, such as relative pellet size, relative RMP amplitude and pedestal collisionality are kept at the ITER like values. Refuelling of density pump out requires a factor of two increase of nominal fuelling rate. Energy confinement and pedestal temperatures are not restored to pre-RMP values by pellet refuelling.

preprint2014arXiv

An experimental investigation on the high density transition of the Scrape-off Layer transport in ASDEX Upgrade

A multidiagnostic approach, utilizing Langmuir probes in the midplane, X-point and divertor walls, along with Lithium beam and infrared measurements is employed to evaluate the evolution of the Scrape-off Layer (SOL) of ASDEX Upgrade across the L-mode density transition leading to the formation of a density shoulder. The flattening of the SOL density profiles is linked to a regime change of filaments, which become faster and larger, and to a similar flattening of the $q_{\parallel}$ profile. This transition is related to the beginning of outer divertor detachment and leads to the onset of a velocity shear layer in the SOL. Experimental measurements are in good agreement with several filament models which describe the process as a transition from conduction to convection-dominated SOL perpendicular transport caused by an increase of parallel collisionality. These results could be of great relevance since both ITER and DEMO will feature detached divertors and densities largely over the transition values, and might therefore exhibit convective transport levels different to those observed typically in present-day devices.

preprint2013arXiv

Applications of Large Eddy Simulation methods to gyrokinetic turbulence

The Large Eddy Simulation (LES) approach - solving numerically the large scales of a turbulent system and accounting for the small-scale influence through a model - is applied to nonlinear gyrokinetic systems that are driven by a number of different microinstabilities. Comparisons between modeled, lower resolution, and higher resolution simulations are performed for an experimental measurable quantity, the electron density fluctuation spectrum. Moreover, the validation and applicability of LES is demonstrated through a series of diagnostics based on the free energetics of the system.

preprint2013arXiv

Improved Collective Thomson Scattering measurements of fast ions at ASDEX Upgrade

Understanding the behaviour of the confined fast ions is important in both current and future fusion experiments. These ions play a key role in heating the plasma and will be crucial for achieving conditions for burning plasma in next-step fusion devices. Microwave-based Collective Thomson Scattering (CTS) is well suited for reactor conditions and offers such an opportunity by providing measurements of the confined fast-ion distribution function resolved in space, time and 1D velocity space. We currently operate a CTS system at ASDEX Upgrade using a gyrotron which generates probing radiation at 105 GHz. A new setup using two independent receiver systems has enabled improved subtraction of the background signal, and hence the first accurate characterization of fast-ion properties. Here we review this new dual-receiver CTS setup and present results on fast-ion measurements based on the improved background characterization. These results have been obtained both with and without NBI heating, and with the measurement volume located close to the centre of the plasma. The measurements agree quantitatively with predictions of numerical simulations. Hence, CTS studies of fast-ion dynamics at ASDEX Upgrade are now feasible. The new background subtraction technique could be important for the design of CTS systems in other fusion experiments.

preprint2013arXiv

Multi-mode Alfvénic Fast Particle Transport and Losses: Numerical vs. Experimental Observation

In many discharges at ASDEX Upgrade fast particle losses can be observed due to Alfvénic gap modes, Reversed Shear Alfvén Eigenmodes or core-localized Beta Alfvén Eigenmodes. For the first time, simulations of experimental conditions in the ASDEX Upgrade fusion device are performed for different plasma equilibria (particularly for different, also non-monotonic q profiles). The numerical tool is the extended version of the HAGIS code [Pinches'98, Brüdgam PhD Thesis, 2010], which also computes the particle motion in the vacuum region between vessel wall in addition to the internal plasma volume. For this work, a consistent fast particle distribution function was implemented to represent the strongly anisotropic fast particle population as generated by ICRH minority heating. Furthermore, HAGIS was extended to use more realistic eigenfunctions, calculated by the gyrokinetic eigenvalue solver LIGKA [Lauber'07]. The main aim of these simulations is to allow fast ion loss measurements to be interpreted with a theoretical basis. Fast particle losses are modeled and directly compared with experimental measurements [García-Muñoz'10]. The phase space distribution and the mode-correlation signature of the fast particle losses allows them to be characterized as prompt, resonant or diffusive (non-resonant). The experimental findings are reproduced numerically. It is found that a large number of diffuse losses occur in the lower energy range (at around 1/3 of the birth energy) particularly in multiple mode scenarios (with different mode frequencies), due to a phase space overlap of resonances leading to a so-called domino [Berk'95] transport process. In inverted q profile equilibria, the combination of radially extended global modes and large particle orbits leads to losses with energies down to 1/10th of the birth energy.

preprint2013arXiv

Sawtooth control using electron cyclotron current drive in the presence of energetic particles in high performance ASDEX Upgrade plasmas

Sawtooth control using steerable electron cyclotron current drive (ECCD) has been demonstrated in ASDEX Upgrade plasmas with a significant population of energetic ions in the plasma core and long uncontrolled sawtooth periods. The sawtooth period is found to be minimised when the ECCD resonance is swept to just inside the q = 1 surface. By utilising ECCD inside q = 1 for sawtooth control, it is possible to avoid the triggering of neoclassical tearing modes, even at significnatly higher pressure than anticipated in the ITER baseline scenario. Operation at 25% higher normalised pressure has been achieved when only modest ECCD power is used for sawtooth control compared to identical discharges without sawtooth control when neo-classical tearing modes are triggered by the sawteeth. Modelling suggests that the destabilisation arising from the change in the local magnetic shear caused by the ECCD is able to compete with the stabilising influence of the energetic particles inside the q = 1 surface.

preprint2012arXiv

Reduced-MHD Simulations of Toroidally and Poloidally Localized ELMs

We use the non-linear reduced-MHD code JOREK to study ELMs in the geometry of the ASDEX Upgrade tokamak. Toroidal mode numbers, poloidal filament sizes, and radial propagation speeds of filaments into the scrape-off layer are in good agreement with observations for type-I ELMs in ASDEX Upgrade. The observed instabilities exhibit a localization of perturbations which is compatible with the "solitary magnetic perturbations" recently discovered in ASDEX Upgrade [R.Wenninger et.al., Solitary Magnetic Perturbations at the ELM Onset, Nucl.Fusion, submitted]. This localization can only be described in numerical simulations with high toroidal resolution.

preprint2012arXiv

Solitary magnetic perturbations at the ELM onset

Edge localised modes (ELMs) allow maintaining sufficient purity of tokamak H-mode plasmas and thus enable stationary H-mode. On the other hand in a future device ELMs may cause divertor power flux densities far in excess of tolerable material limits. The size of the energy loss per ELM is determined by saturation effects in the non-linear phase of the ELM, which at present is hardly understood. Solitary magnetic perturbations (SMPs) are identified as dominant features in the radial magnetic fluctuations below 100kHz. They are typically observed close (+-0.1ms) to the onset of pedestal erosion. SMPs are field aligned structures rotating in the electron diamagnetic drift direction with perpendicular velocities of about 10km/s. A comparison of perpendicular velocities suggests that the perturbation evoking SMPs is located at or inside the separatrix. Analysis of very pronounced examples showed that the number of peaks per toroidal turn is 1 or 2, which is clearly lower than corresponding numbers in linear stability calculations. In combination with strong peaking of the magnetic signals this results in a solitary appearance resembling modes like palm tree modes, edge snakes or outer modes. This behavior has been quantified as solitariness and correlated to main plasma parameters. SMPs may be considered as a signature of the non-linear ELM-phase originating at the separatrix or further inside. Thus they provide a handle to investigate the transition from linear to non-linear ELM phase. By comparison with data from gas puff imaging processes in the non-linear phase at or inside the separatrix and in the scrape-off-layer (SOL) can be correlated. A connection between the passing of an SMP and the onset of radial filament propagation has been found. Eventually the findings related to SMPs may contribute to a future quantitative understanding of the non-linear ELM evolution.