Engineering topics

Grandgirard, V.

Publications and source records attributed to Grandgirard, V..

10 records

DOE OSTI2018

Particle transport due to energetic-particle-driven geodesic acoustic modes

The transport of particles in the presence of energetic geodesic acoustic modes (EGAMs) is analysed by means of full-f global gyro-kinetic simulations, using the multi-species version of the code and a test particle tracking post-treatment, which solves the equations of motion of passive gyro-centres embedded in the self-consistent EGAM potential obtained from the simulations. It is found that EGAMs induce the transport of particles, which eventually results in counter-passing particle losses modulated at the EGAM frequency. A detailed analysis of the trajectories of the test gyro-centres is performed and evidence of the interaction between the EGAM island and the region of magnetically trapped particles (trapping cone) is provided. In particular, we report for the first time on the complex interaction between the stochastic separatrix of the EGAM island and the X-point of the trapping cone, creating a channel for the transport of particles from the $v_{||} < 0$ to the $v_{||} > 0$ regions. Therefore, for the cases analysed in this work, the co-injection of energetic particles might lead to a significant reduction of EGAM-induced losses. This result opens up new perspectives for further studies of the selection of energetic particle injection in order to minimize the losses due to EGAMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Intercode comparison of gyrokinetic global electromagnetic modes

Aiming to fill a corresponding lack of sophisticated test cases for global electromagnetic gyrokinetic codes, a new hierarchical benchmark is proposed. Starting from established test sets with adiabatic electrons, fully gyrokinetic electrons, and electrostatic fluctuations are taken into account before finally studying the global electromagnetic micro-instabilities. Results from up to five codes involving representatives from different numerical approaches as particle-in-cell methods, Eulerian and Semi-Lagrangian are shown. By means of spectrally resolved growth rates and frequencies and mode structure comparisons, agreement can be confirmed on ion-gyro-radius scales, thus providing confidence in the correct implementation of the underlying equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

A multi-species collisional operator for full-F gyrokinetics

A linearized multi-species collision operator has been developed for an efficient implementation in gyrokinetic codes. This operator satisfies the main expected properties: particle, momentum, and energy conservation, and existence of an H-theorem. A gyrokinetic version is then calculated, which involves derivatives with respect to the gyrocenter position, parallel velocity, and magnetic momentum. An isotropic version in the velocity space can be constructed for the specific problem of trace impurities colliding with a main species. A simpler version that involves derivatives with parallel velocity only has been developed. This reduced version has been implemented in the GYSELA gyrokinetic code, and is shown to comply with particle, momentum, and energy conservation laws. Moreover, the interspecies relaxation rates for momentum and energy agree very well with the theoretical values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Comprehensive comparisons of geodesic acoustic mode characteristics and dynamics between Tore Supra experiments and gyrokinetic simulations

In a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Shear flow instabilities induced by trapped ion modes in collisionless temperature gradient turbulence

One important issue in turbulence self-organization is the interplay between the Kelvin–Helmholtz (KH) instability and streamers and/or zonal flows. This question has been debated for a long time. The effects of the KH instability and its position in the sequence of events between streamers, turbulence, and zonal flows have been investigated with a reduced gyro-bounce averaged kinetic code devoted to study the primary ion temperature gradient (ITG) instability linked to trapped ion modes (TIM). In toroidal geometry, the specific dynamics of TIM linked to trapped particles becomes important when the frequency of ITG modes falls below the ion bounce frequency, allowing one to average on both the cyclotron and bounce motion fast time scales. This reduction of the number of degrees of freedom leads to a strong reduction of computer resources (memory and computation time). Bounce-averaged gyrokinetic code can be considered as a toy model able to describe basic structures of turbulent transport in tokamak devices. In particular, by means of this code, we have observed that the energy injected in the system by the TIM instability is exchanged between streamers and zonal flows by means of KH vortices that grow along these structures in the nonlinear phase. The energy transfer occurs throughout the relaxation phase of the streamer growth leading to a modification of the KH modes and to the generation of the zonal flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Non-linear magnetohydrodynamic modeling of plasma response to resonant magnetic perturbations

The interaction of static Resonant Magnetic Perturbations (RMPs) with the plasma flows is modeled in toroidal geometry, using the non-linear resistive MHD code JOREK, which includes the X-point and the scrape-off-layer. Two-fluid diamagnetic effects, the neoclassical poloidal friction and a source of toroidal rotation are introduced in the model to describe realistic plasma flows. RMP penetration is studied taking self-consistently into account the effects of these flows and the radial electric field evolution. JET-like, MAST, and ITER parameters are used in modeling. For JET-like parameters, three regimes of plasma response are found depending on the plasma resistivity and the diamagnetic rotation: at high resistivity and slow rotation, the islands generated by the RMPs at the edge resonant surfaces rotate in the ion diamagnetic direction and their size oscillates. At faster rotation, the generated islands are static and are more screened by the plasma. An intermediate regime with static islands which slightly oscillate is found at lower resistivity. In ITER simulations, the RMPs generate static islands, which forms an ergodic layer at the very edge (ψ≥0.96) characterized by lobe structures near the X-point and results in a small strike point splitting on the divertor targets. In MAST Double Null Divertor geometry, lobes are also found near the X-point and the 3D-deformation of the density and temperature profiles is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Turbulent acceleration and heating in toroidal magnetized plasmas

It is shown that turbulence is responsible for a source of momentum, which cannot be recast as a divergence of a momentum flux. This process is similar to turbulent heating, with similar properties. The sum over all species vanishes up to polarization contributions. Hence, toroidal momentum is transferred from species to species, mediated by turbulence. As for momentum flux, symmetry breaking is needed. Flow shear is investigated as a source of symmetry breaking, leading to a source of momentum proportional to the shear rate. Turbulent acceleration is significant for ion species. It is found that it is proportional to the charge number Z, while turbulent heating scales as Z{sup 2}/A, where A is the mass number. It is maximum in the edge, where the E × B flow shear rate and turbulence intensity are maximum. When both are large enough, the turbulent torque may overcome the collisional friction between impurities and main ions, thus leading to different toroidal velocities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

Physics of intrinsic rotation in flux-driven ITG turbulence

Global, heat flux-driven ITG gyrokinetic simulations which manifest the formation of macroscopic, mean toroidal flow profiles with peak thermal Mach number 0.05, are reported. Both a particle-in-cell (XGC1p) and a semi-Lagrangian (GYSELA) approach are utilized without a priori assumptions of scale separation between turbulence and mean fields. Flux-driven ITG simulations with different edge flow boundary conditions show in both approaches the development of net unidirectional intrinsic rotation in the co-current direction. Intrinsic torque is shown to scale approximately linearly with the inverse scale length of the ion temperature gradient. External momentum input is shown to effectively cancel the intrinsic rotation profile, thus confirming the existence of a local residual stress and intrinsic torque. Fluctuation intensity, intrinsic torque and mean flow are demonstrated to develop inwards from the boundary. The measured correlations between residual stress and two fluctuation spectrum symmetry breakers, namely E × B shear and intensity gradient, are similar. Avalanches of (positive) heat flux, which propagate either outwards or inwards, are correlated with avalanches of (negative) parallel momentum flux, so that outward transport of heat and inward transport of parallel momentum are correlated and mediated by avalanches. The probability distribution functions of the outward heat flux and the inward momentum flux show strong structural similarity.

Physics↗

DOE OSTI2012

Physics of Intrinsic Rotation in Flux-Driven ITG Turbulence

Global, heat flux-driven ITG gyrokinetic simulations which manifest the formation of macroscopic, mean toroidal flow profiles with peak thermal Mach number 0.05, are reported. Both a particle-in-cell (XGC1p) and a semi-Lagrangian (GYSELA) approach are utilized without a priori assumptions of scale-separation between turbulence and mean fields. Flux-driven ITG simulations with different edge flow boundary conditions show in both approaches the development of net unidirectional intrinsic rotation in the co-current direction. Intrinsic torque is shown to scale approximately linearly with the inverse scale length of the ion temperature gradient. External momentum input is shown to effectively cancel the intrinsic rotation profile, thus confirming the existence of a local residual stress and intrinsic torque. Fluctuation intensity, intrinsic torque and mean flow are demonstrated to develop inwards from the boundary. The measured correlations between residual stress and two fluctuation spectrum symmetry breakers, namely E x B shear and intensity gradient, are similar. Avalanches of (positive) heat flux, which propagate either outwards or inwards, are correlated with avalanches of (negative) parallel momentum flux, so that outward transport of heat and inward transport of parallel momentum are correlated and mediated by avalanches. The probability distribution functions of the outward heat flux and the inward momentum flux show strong structural similarity

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

Neoclassical physics in full distribution function gyrokinetics

Treatment of binary Coulomb collisions when the full gyrokinetic distribution function is evolved is reported here. A spectrum of different collision operators is presented, differing through both the physics that can be addressed and the numerics they are based on. Eulerian-like (semi-Lagrangian) and particle in cell (PIC) (Monte-Carlo) schemes are successfully cross-compared, and a detailed confrontation to neoclassical theory is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗