Source author record

Philipp Lauber

Philipp Lauber 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

7works
2topics
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

7 published item(s)

preprint2022arXiv

Gyrokinetic modelling of anisotropic energetic particle driven instabilities in tokamak plasmas

Energetic particles produced by neutral beams are observed to excite energetic-particle-driven geodesic acoustic modes (EGAMs) in tokamaks. We study the effects of anisotropy of distribution function of the energetic particles on the excitation of such instabilities with ORB5, a gyrokinetic particle-in-cell code. Numerical results are shown for linear electrostatic simulations with ORB5. The growth rate is found to be sensitively dependent on the phase-space shape of the distribution function. The behavior of the instability is qualitatively compared to the theoretical analysis of dispersion relations. Realistic neutral beam energetic particle anisotropic distributions are obtained from the heating solver RABBIT and are introduced into ORB5 as input distribution function. Results show a dependence of the growth rate on the injection angle. A qualitative comparison to experimental measurements is presented and few disagreements between them are found, being the growth rate in the simulations much lower than that from experiments. An explanation for the difference is advanced.

preprint2022arXiv

One dimensional reduced model for ITER relevant energetic particle transport

We set up a mapping procedure able to translate the evolution of the radial profile of fast ions, interacting with Toroidal Alfvén Eigenmodes, into the dynamics of an equivalent one dimensional bump-on-tail system. We apply this mapping technique to reproduce ITER relevant simulations, which clearly outlined deviations from the diffusive quasi-linear (QL) model. Our analysis demonstrates the capability of the one-dimensional beam-plasma dynamics to predict the relevant features of the non-linear hybrid LIGKA/HAGIS simulations. In particular, we clearly identify how the deviation from the QL evolutive profiles is due to the presence of avalanche processes. A detailed analysis regarding the reduced dimensionality is also addressed, by means of phase-space slicing based on constants of motion. In the conclusions, we outline the main criticalities and outcomes of the procedure, which must be satisfactorily addressed to make quantitative prediction on the observed outgoing fluxes in a Tokamak device.

preprint2022arXiv

Representation and modeling of charged particle distributions in tokamaks

Experimental diagnostics, analysis tools and simulations represent particle distributions in various forms and coordinates. Algorithms to manage these data are needed on platforms like the ITER Integrated Modelling & Analysis Suite (IMAS), performing tasks such as archiving, modeling, conversion and visualization. A method that accomplishes some of the required tasks for distributions of charged particles with arbitrarily large magnetic drifts in axisymmetric tokamak geometry is described here. Given a magnetic configuration, we first construct a database of guiding center orbits, which serves as a basis for representing particle distributions. The orbit database contains the geometric information needed to perform conversions between arbitrary coordinates, modeling tasks, and resonance analyses. Using that database, an imported or newly modeled distribution is mapped to an exact equilibrium, where the dimensionality is reduced to three constants of motion (CoM). The orbit weight is uniquely given when the input is a true distribution: one that measures the number of physical particles per unit of phase space volume. Less ideal inputs, such as distributions estimated without drifts, or models of particle sources, can also be processed. As an application example, we reconstruct the drift-induced features of a distribution of fusion-born alpha particles in a large tokamak, given only a birth profile, which is not a function of the alpha's CoM. Repeated back-and-forth transformations between CoM space and energy-pitch-cylinder coordinates are performed for verification and as a proof of principle for IMAS.

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

Numerical validation of the electromagnetic gyrokinetic code NEMORB on global axisymmetric modes

This is a report about a comparison of collisionless simulations on global modes (i.e. low poloidal mode number) with gyrokinetic code NEMORB against analytical theory and other codes. Only axisymmetric modes, i.e. with toroidal mode number n=0, are considered, and flat equilibrium profiles. Benchmarks are performed for GAMs against local analytical theory. In the presence of energetic ions, local benchmarks of NEMORB are performed against semilagrangian gyrokinetic code GYSELA. The models of adiabatic vs trapped-kinetic- vs fully-kinetic-electrons and of electrostatic vs electromagnetic at very low beta are compared. Scalings of Alfven modes are also presented.

preprint2012arXiv

Double-resonant fast particle-wave interaction

In future fusion devices fast particles must be well confined in order to transfer their energy to the background plasma. Magnetohydrodynamic instabilities like Toroidal Alfvén Eigenmodes or core-localized modes such as Beta Induced Alfvén Eigenmodes and Reversed Shear Alfvén Eigenmodes, both driven by fast particles, can lead to significant losses. This is observed in many ASDEX Upgrade discharges. The present study applies the drift-kinetic HAGIS code with the aim of understanding the underlying resonance mechanisms, especially in the presence of multiple modes with different frequencies. Of particular interest is the resonant interaction of particles simultaneously with two different modes, referred to as 'double-resonance'. Various mode overlapping scenarios with different q profiles are considered. It is found that, depending on the radial mode distance, double-resonance is able to enhance growth rates as well as mode amplitudes significantly. Surprisingly, no radial mode overlap is necessary for this effect. Quite the contrary is found: small radial mode distances can lead to strong nonlinear mode stabilization of a linearly dominant mode.

preprint2012arXiv

Nonlinear alfvénic fast particle transport and losses

Magnetohydrodynamic instabilities like Toroidal Alfvén Eigenmodes or core-localized modes such as Beta Induced Alfvén Eigenmodes and Reversed Shear Alfvén Eigenmodes driven by fast particles can lead to significant redistribution and losses in fusion devices. This is observed in many ASDEX Upgrade discharges. The present work aims to understand the underlying resonance mechanisms, especially in the presence of multiple modes with different frequencies. Resonant mode coupling mechanisms are investigated using the drift kinetic HAGIS code [Pinches 1998]. Simulations were performed for different plasma equilibria, in particular for different q profiles, employing the availability of improved experimental data. A study was carried out, investigating double-resonant mode coupling with respect to various overlapping scenarios. It was found that, depending on the radial mode distance, double-resonance is able to enhance growth rates as well as mode amplitudes significantly. Small radial mode distances, however can also lead to strong nonlinear mode stabilization of a linear dominant mode. With the extended version of HAGIS, losses were simulated and directly compared with experimental loss measurements. The losses' phase space distribution as well as their ejection signal is consistent with experimental data. Furthermore, it allowed to characterize them as prompt, resonant or stochastic. It was found that especially in multiple mode scenarios (with different mode frequencies), abundant incoherent losses occur in the lower energy range, due to a broad phase-space stochastization. The incoherent higher energetic losses are "prompt", i.e. their initial energy is too large for confined orbits.