Nonlinear alfv\'enic fast particle transport and losses

Magnetohydrodynamic instabilities like Toroidal Alfv\'en Eigenmodes or core-localized modes such as Beta Induced Alfv\'en Eigenmodes and Reversed Shear Alfv\'en 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.Comment: 7 pages, 6 figures, Reviewed Conference Proceedings (Joint Varenna - Lausanne International Workshop on the Theory of Fusion Plasmas) to be published in IOP's "Journal of Physics: Conference Series

Large-mode-number magnetohydrodynamic instability driven by sheared flows in a tokamak plasma with reversed central shear

The effect of a narrow sub-Alfvenic shear flow layer near the minimum q_min of the tokamak safety factor profile in a configuration with reversed central shear is analyzed. Sufficiently strong velocity shear gives rise to a broad spectrum of fast growing Kelvin-Helmholtz (KH)-like ideal magnetohydrodynamic (MHD) modes with dominant mode numbers m,n ~ 10. Nonlinear simulations with finite resistivity show magnetic reconnection near ripples caused by KH-like vortices, the formation of turbulent structures, and a flattening of the flow profile. The KH modes are compared to double tearing modes (DTM) which dominate at lower shearing rates. The possible application of these results in tokamaks with internal transport barrier is discussed.Comment: 4 pages, 4 figure

Multi-mode Alfv\'enic Fast Particle Transport and Losses: Numerical vs. Experimental Observation

In many discharges at ASDEX Upgrade fast particle losses can be observed due to Alfv\'enic gap modes, Reversed Shear Alfv\'en Eigenmodes or core-localized Beta Alfv\'en 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\"udgam 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\'ia-Mu\~noz'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.Comment: 16 Pages, 17 Figure

A three-dimensional reduced MHD model consistent with full MHD

Within the context of a viscoresistive magnetohydrodynamic (MHD) model with anisotropic heat transport and cross-field mass diffusion, we introduce novel three-term representations for the magnetic field (background vacuum field, field line bending and field compression) and velocity ($\vec E\times\vec B$ flow, field-aligned flow and fluid compression), which are amenable to three-dimensional treatment. Once the representations are inserted into the MHD equations, appropriate projection operators are applied to Faraday's law and the Navier-Stokes equation to obtain a system of scalar equations that is closed by the continuity and energy equations. If the background vacuum field is sufficiently strong and the $\beta$ is low, MHD waves are approximately separated by the three terms in the velocity representation, with each term containing a specific wave. Thus, by setting the appropriate term to zero, we eliminate fast magnetosonic waves, obtaining a reduced MHD model. We also show that the other two velocity terms do not compress the magnetic field, which allows us to set the field compression term to zero within the same reduced model. Dropping also the field-aligned flow, a further simplified model is obtained, leading to a fully consistent hierarchy of reduced and full MHD models for 3D plasma configurations. Finally, we discuss the conservation properties and derive the conditions under which the reduction approximation is valid. We also show that by using an ordering approach, reduced MHD equations similar to what we got from the ansatz approach can be obtained by means of a physics-based asymptotic expansion.Comment: 18 pages. This article was published in Physics of Plasma

Detection of DC currents and resistance measurements in longitudinal spin Seebeck effect experiments on Pt/YIG and Pt/NFO

In this work we investigated thin films of the ferrimagnetic insulators YIG and NFO capped with thin Pt layers in terms of the longitudinal spin Seebeck effect (LSSE). The electric response detected in the Pt layer under an out-of-plane temperature gradient can be interpreted as a pure spin current converted into a charge current via the inverse spin Hall effect. Typically, the transverse voltage is the quantity investigated in LSSE measurements (in the range of \mu V). Here, we present the directly detected DC current (in the range of nA) as an alternative quantity. Furthermore, we investigate the resistance of the Pt layer in the LSSE configuration. We found an influence of the test current on the resistance. The typical shape of the LSSE curve varies for increasing test currents.Comment: 4 pages, 2 figure

How well can VMEC predict the initial saturation of external kink modes in near circular tokamaks and l=2l=2 stellarators?

The equilibrium code, VMEC, is used to study external kinks in low $\beta$ tokamaks and $l=2$ stellarators. The applicability of the code when modelling nonlinear MHD effects is explored in an attempt to understand and predict how the initial saturation of the MHD mode depends on the external rotational transform. It is shown that helicity preserving, free boundary VMEC computations do not converge to a single perturbed solution with increasing spectral resolution. Additional constraints are therefore applied to narrow down the numerical resolution parameters appropriate for physical scans. The dependence of the modelled (4, 1) kink mode on the external rotational transform and field periodicity is then studied. While saturated states can be identified which decrease in amplitude with increasing external rotational transform, bifurcated states are found that contradict this trend. It was therefore not possible to use VMEC alone to identify the physical dependency of the nonlinear mode amplitude on the magnetic geometry. The accuracy of the VMEC solutions is nevertheless demonstrated by showing that the expected toroidal mode coupling is captured in the magnetic energy spectrum for stellarator cases. Comparing with the initial value code, JOREK, the predicted redistribution of poloidal magnetic energy from the vacuum to plasma region in VMEC is shown to be physical. This work is a first step towards using VMEC to study MHD modes in stellarator geometry.Comment: Submitted to Physics of Plasmas. The submission has been modified according to reviewer comment

MHD simulations of formation, sustainment and loss of Quiescent H-mode in the all-tungsten ASDEX Upgrade

Periodic edge localized modes (ELMs) are the non-linear consequences of pressure-gradient-driven ballooning modes and current-driven peeling modes becoming unstable in the pedestal region of high confinement fusion plasmas. In future tokamaks like ITER, large ELMs are foreseen to severely affect the lifetime of wall components as they transiently deposit large amounts of heat onto a narrow region at the divertor targets. Several strategies exist for avoidance, suppression, or mitigation of these instabilities, such as the naturally ELM-free quiescent H-mode (QH-mode). In the present article, an ASDEX Upgrade equilibrium that features a QH-mode is investigated through non-linear extended MHD simulations covering the dynamics over tens of milliseconds. The equilibrium is close to the ideal peeling limit and non-linearly develops saturated modes at the edge of the plasma. A dominant toroidal mode number of $n=1$ is found, for which the characteristic features of the edge harmonic oscillation are recovered. The saturated modes contribute to heat and particle transport preventing pedestal build-up to the ELM triggering threshold. The non-linear dynamics of the mode, in particular its interaction with the evolution of the edge safety factor is studied, which suggest a possible new saturation mechanism for the QH-mode. The simulations show good qualitative and quantitative agreement to experiments in AUG. In particular, the processes leading to the termination of QH-mode above a density threshold is studied, which results in the transition into an ELM regime. In the vicinity of this threshold, limit cycle oscillations are observed.Comment: Revised version with modifications from review process include

CAMELS-CH: hydro-meteorological time series and landscape attributes for 331 catchments in hydrologic Switzerland

We present CAMELS-CH (Catchment Attributes and MEteorology for Large-sample Studies – Switzerland), a large-sample hydro-meteorological data set for hydrologic Switzerland in central Europe. This domain covers 331 basins within Switzerland and neighboring countries. About one-third of the catchments are located in Austria, France, Germany and Italy. As an Alpine country, Switzerland covers a vast diversity of landscapes, including mountainous environments, karstic regions, and several strongly cultivated regions, along with a wide range of hydrological regimes, i.e., catchments that are glacier-, snow- or rain dominated. Similar to existing data sets, CAMELS-CH comprises dynamic hydro-meteorological variables and static catchment attributes. CAMELS-CH (Höge et al., 2023; available at https://doi.org/10.5281/zenodo.7784632) encompasses 40 years of data between 1 January 1981 and 31 December 2020, including daily time series of stream flow and water levels, and of meteorological data such as precipitation and air temperature. It also includes daily snow water equivalent data for each catchment starting from 2 September 1998. Additionally, we provide annual time series of land cover change and glacier evolution per catchment. The static catchment attributes cover location and topography, climate, hydrology, soil, hydrogeology, geology, land use, human impact and glaciers. This Swiss data set complements comparable publicly accessible data sets, providing data from the “water tower of Europe”

Nanofasern als neuartige Träger für flüchtige Signalstoffe zur biotechnischen Regulierung von Schadinsekten im integrierten und ökologischen Landbau

Using nanofibers as dispensers for pheromones and kairomones in plant protection for disrupting insect chemical communication is a novel approach aiming at popularizing this technique in organic and integrated plant production. Expected advantages of the nanofibers are highly controlled spatiotemporal release rates of pheromones / kairomones, improved climatic stability, and mechanized application. Dispenser types used so far show deficiencies in one or more of these requirements. Mechanical application of pheromones is a new approach to reduce the costs of manual labour and therefore the environmentally compatible, highly specific and efficient technique of mating disruption may become an alternative to the use of synthetic pesticides in integrated pest management. The nanofibers are highly elastic, which prevents breakoff of smaller pieces, and polymers used are biocompatible. Due to the scale of nanofibers the mass input both for pheromones and for polymeric nanofibers is very low. Major environmental benefits are high control specificity, very low concentrations of residues and reduced risk of development towards insect resistance