Showing 1–36 of 36 results for author: Garbet, X

Using Convolutional Neural Networks to detect Edge Localized Modes in DIII-D from Doppler Backscattering measurements

Authors: N. Q. X. Teo, V. H. Hall-Chen, K. Barada, R. J. H. Ng, L. Gu, A. K. Yeoh, Q. T. Pratt, X. Garbet, T. L. Rhodes

Abstract: In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods $\unicode{x2013}$ D-alpha spectr… ▽ More In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods $\unicode{x2013}$ D-alpha spectroscopy and Doppler backscattering (DBS). The former detects ELMs by measuring Balmer alpha emission while the latter uses microwave radiation to probe the plasma. DBS has the advantage of having higher temporal resolution and robustness to damage. These advantages of DBS diagnostics may be beneficial for future operational tokamaks and thus data processing techniques for DBS should be developed in preparation. In sight of this, we explore the training of neural networks to detect ELMs from DBS data, using D-alpha data as the ground truth. With shots found in the DIII-D database, the model is trained to classify each time step based on the occurrence of an ELM event. The results are promising. When tested on shots similar to those used for training, the model is capable of consistently achieving a high f1-score of 0.93. This score is a performance metric for imbalanced datasets that ranges between 0 and 1. We evaluate the performance of our neural network on a variety of ELMs $\unicode{x2013}$ grasssy, suppressed, and wide pedestal $\unicode{x2013}$ finding broad applicability. Beyond ELMs, our work demonstrates the wider feasibility of applying neural networks to data from DBS diagnostics. △ Less

Submitted 3 July, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Transition to super-diffusive transport in turbulent plasmas

Authors: Matteo Stanzani, Filippo Arlotti, Guido Ciraolo, Xavier Garbet, Cristel Chandre

Abstract: We investigate the motion of charged particles in a turbulent electrostatic potential using guiding-center theory. By increasing the Larmor radius, the dynamics exhibit close-to-ballistic transport properties. The transition from diffusive to ballistic transport is analyzed using nonlinear dynamics. It is found that twistless invariant tori in the guiding-center dynamics are responsible for this t… ▽ More We investigate the motion of charged particles in a turbulent electrostatic potential using guiding-center theory. By increasing the Larmor radius, the dynamics exhibit close-to-ballistic transport properties. The transition from diffusive to ballistic transport is analyzed using nonlinear dynamics. It is found that twistless invariant tori in the guiding-center dynamics are responsible for this transition, drastically affecting transport properties of charged particles. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Self-consistent gyrokinetic modelling of turbulent and neoclassical tungsten transport in toroidally rotating plasmas

Authors: Kyungtak Lim, Xavier Garbet, Yanick Sarazin, Etienne Gravier, Maxime Lesur, Guillaume Lo-Cascio, Timothe Rouyer

Abstract: The effect of toroidal rotation on both turbulent and neoclassical transport of tungsten (W) in tokamaks is investigated using the flux-driven, global, nonlinear 5D gyrokinetic code GYSELA. Nonlinear simulations are carried out with different levels of momentum injection that drive W to the supersonic regime, while the toroidal velocity of the main ions remains in the subsonic regime. The numerica… ▽ More The effect of toroidal rotation on both turbulent and neoclassical transport of tungsten (W) in tokamaks is investigated using the flux-driven, global, nonlinear 5D gyrokinetic code GYSELA. Nonlinear simulations are carried out with different levels of momentum injection that drive W to the supersonic regime, while the toroidal velocity of the main ions remains in the subsonic regime. The numerical simulations demonstrate that toroidal rotation induces centrifugal forces that cause W to accumulate in the outboard region, generating an in-out poloidal asymmetry. This asymmetry enhances neoclassical inward convection, which can lead to central accumulation of W in cases of strong plasma rotation. The core accumulation of W is mainly driven by inward neoclassical convection. However, as momentum injection continues, roto-diffusion, proportional to the radial gradient of the toroidal velocity, becomes significant and generate outward turbulent flux in the case of ion temperature gradient (ITG) turbulence. Overall, the numerical results from nonlinear GYSELA simulations are in qualitative agreement with the theoretical predictions for impurity transport, as well as experimental observations. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 26 pages, 10 figures, to be published

Global linear stability analysis of kinetic Trapped Ion Mode (TIM) turbulence in tokamak plasma using spectral method

Authors: D. Mandal, M. Lesur, E. Gravier, J. N. Sama, A. Guillevic, Y. Sarazin, X. Garbet

Abstract: Trapped ion modes (TIM) belong to the family of ion temperature gradient (ITG) modes, which are one of the important ingredients in heat turbulent transport at the ion scale in tokamak plasmas. It is essential to properly estimate their linear growth rate to understand their influence on ion-scale turbulent transport. A global linear analysis of a reduced gyro-bounce kinetic model for trapped part… ▽ More Trapped ion modes (TIM) belong to the family of ion temperature gradient (ITG) modes, which are one of the important ingredients in heat turbulent transport at the ion scale in tokamak plasmas. It is essential to properly estimate their linear growth rate to understand their influence on ion-scale turbulent transport. A global linear analysis of a reduced gyro-bounce kinetic model for trapped particle modes is performed, and a spectral method is proposed to solve the dispersion relation. Importantly, the radial profile of the particle drift velocity is taken into account in the linear analysis by considering the magnetic flux ψ dependency of the equilibrium Hamiltonian H_{eq}(ψ) in the quasi-neutrality equation and equilibrium gyro-bounce averaged distribution function F_{eq} . Using this spectral method, linear growth-rates of TIM instability in presence of different temperature profiles and precession frequencies of trapped ions, with an approximated constant Hamiltonian and the exact ψ dependent equilibrium Hamiltonian, are investigated. The growth-rate depends on the logarithmic gradient of temperature κ_{T} , density κ_{n} and equilibrium Hamiltonian κ_Λ . With the exact ψ dependent Hamiltonian, the growth rates and potential profiles are modified significantly, compared to the cases with approximated constant Hamiltonian. All the results from the global linear analysis agree with a semi-Lagrangian based linear Vlasov solver with a good accuracy. This spectral method is very fast and requires very less computation resources compared to a linear version of Vlasov-solver based on a semi-Lagrangian scheme. △ Less

Submitted 15 March, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Journal ref: Plasma Phys. Controlled Fusion 65 (2023), 055001

The Problem of Marginality in Model Reductions of Turbulence

Authors: C. Gillot, G. Dif-Pradalier, Y. Sarazin, C. Bourdelle, Y. Camenen, J. Citrin, X. Garbet, Ph. Ghendrih, V. Grandgirard, P. Manas, F. Widmer

Abstract: Reduced quasilinear (QL) and nonlinear (gradient-driven) models with scale separations, commonly used to interpret experiments and to forecast turbulent transport levels in magnetised plasmas are tested against nonlinear models without scale separations (flux-driven). Two distinct regimes of turbulence -- either far above threshold or near marginal stability -- are investigated with Boltzmann elec… ▽ More Reduced quasilinear (QL) and nonlinear (gradient-driven) models with scale separations, commonly used to interpret experiments and to forecast turbulent transport levels in magnetised plasmas are tested against nonlinear models without scale separations (flux-driven). Two distinct regimes of turbulence -- either far above threshold or near marginal stability -- are investigated with Boltzmann electrons. The success of reduced models especially hinges on the reproduction of nonlinear fluxes. Good agreement between models is found above threshold whilst reduced models would significantly underpredict fluxes near marginality, overlooking mesoscale flow organisation and turbulence self-advection. Constructive prescriptions whereby to improve reduced models is discussed. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Full Self-Consistent Vlasov-Maxwell Solution

Authors: Aurélien Cordonnier, Xavier Leoncini, Guilhem Dif-Pradalier, Xavier Garbet

Abstract: Full self-consistent stationary Vlasov-Maxwell solutions of magnetically confined plasmas are built for systems with cylindrical symmetries. The stationary solutions are thermodynamic equilibrium solutions. These are obtained by computing the equilibrium distribution function resulting from maximizing the entropy and closing the equations with source terms that are then computed by using the obtai… ▽ More Full self-consistent stationary Vlasov-Maxwell solutions of magnetically confined plasmas are built for systems with cylindrical symmetries. The stationary solutions are thermodynamic equilibrium solutions. These are obtained by computing the equilibrium distribution function resulting from maximizing the entropy and closing the equations with source terms that are then computed by using the obtained distribution. This leads to a self-consistent problem corresponding to solving a set of two coupled second order non-linear differential equations. Relevant plasma parameters are introduced and a bifurcation leading to an improvement of plasma confinement is shown. Conversely in the improved confinement setting, we exhibit the emergence of a separatrix in the integrable motion of a charged particles. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Transport barrier onset and edge turbulence shortfall in fusion plasmas

Authors: Guilhem Dif-Pradalier, Philippe Ghendrih, Yanick Sarazin, Yann Camenen, Peter Donnel, Xavier Garbet, Virginie Grandgirard, Yann Munschy, Laure Vermare, Fabien Widmer

Abstract: Turbulent plasmas notably self-organize to higher energy states upon application of additional free energy sources or modification of edge operating conditions. Mechanisms whereby such bifurcations occur have been actively debated for decades. Enhanced confinement occurs at the plasma edge, where a shortfall of predicted turbulence intensity has been puzzling scientists for decades. We show, from… ▽ More Turbulent plasmas notably self-organize to higher energy states upon application of additional free energy sources or modification of edge operating conditions. Mechanisms whereby such bifurcations occur have been actively debated for decades. Enhanced confinement occurs at the plasma edge, where a shortfall of predicted turbulence intensity has been puzzling scientists for decades. We show, from the primitive kinetic equations that both problems are connected and that interplay of confined plasma turbulence with its material boundaries is essential to curing the shortfall of predicted turbulence and to triggering spontaneous transport barrier onset at the plasma edge. Both problems determine access to improved confinement and are central to fusion research. A comprehensive discussion of the underlying mechanisms is proposed. These results, highly relevant to the quest for magnetic fusion may also be generic to many problems in fluids and plasmas where turbulence self-advection is active. △ Less

Submitted 19 September, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Journal ref: COMMUNICATIONS PHYSICS 5:229 (2022)

Quasilinear gyrokinetic theory: A derivation of QuaLiKiz

Authors: Cole Darin Stephens, Xavier Garbet, Jonathan Citrin, Clarisse Bourdelle, Karel Lucas van de Plassche, Frank Jenko

Abstract: In order to predict and analyze turbulent transport in tokamaks, it is important to model transport that arises from microinstabilities. For this task, quasilinear codes have been developed that seek to calculate particle, angular momentum, and heat fluxes both quickly and accurately. In this tutorial, we present a derivation of one such code known as QuaLiKiz, a quasilinear gyrokinetic transport… ▽ More In order to predict and analyze turbulent transport in tokamaks, it is important to model transport that arises from microinstabilities. For this task, quasilinear codes have been developed that seek to calculate particle, angular momentum, and heat fluxes both quickly and accurately. In this tutorial, we present a derivation of one such code known as QuaLiKiz, a quasilinear gyrokinetic transport code. The goal of this derivation is to provide a self-contained and complete description of the underlying physics and mathematics of QuaLiKiz from first principles. This work serves both as a comprehensive overview of QuaLiKiz specifically as well as an illustration for deriving quasilinear models in general. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 52 pages

Journal ref: Journal of Plasma Physics, Volume 87, Issue 4, August 2021, 905870409

Nonlinear dynamics of the fishbone-induced alpha transport on ITER

Authors: G. Brochard, R. Dumont, H. Lütjens, X. Garbet, T. Nicolas, P. Maget

Abstract: The fishbone-induced transport of alpha particles is computed for the ITER 15 MA baseline scenario, using the nonlinear hybrid Kinetic-MHD code XTOR-K. Two limit cases have been studied, in order to analyse the characteristic regimes of the fishbone instability : the weak kinetic drive limit and the strong kinetic drive limit. In both those regimes, characteristic features of the n = m = 1 fishbon… ▽ More The fishbone-induced transport of alpha particles is computed for the ITER 15 MA baseline scenario, using the nonlinear hybrid Kinetic-MHD code XTOR-K. Two limit cases have been studied, in order to analyse the characteristic regimes of the fishbone instability : the weak kinetic drive limit and the strong kinetic drive limit. In both those regimes, characteristic features of the n = m = 1 fishbone instability are recovered, such as a strong up/down-chirping of the mode frequency, associated with a resonant transport of trapped and passing alpha particles. The effects of the n = m = 0 sheared poloidal and toroidal plasma rotation are taken into account in the simulations. The shear is not negligible, which implies that the fishbone mode frequency has a radial dependency, impacting the wave-particle resonance condition. Phase space hole and clump structures are observed in both nonlinear regimes, centered around the precessional and passing resonances. These structures remains attached to the resonances as the different mode frequencies chirp up and down. In the nonlinear phase, the transport of individual resonant trapped particles is identified to be linked to mode-particle synchronization. On this basis, a partial mechanism for the nonlinear coupling between particle transport and mode dominant down-chirping is proposed. The overall transport of alpha particles inside out the q = 1 surface is of order 2-5% of the initial population between the simulations. The loss of alpha power is found to be directly equal to the loss of alpha particles. △ Less

Submitted 12 May, 2020; originally announced May 2020.

En route to realistic modeling of the kinetic-MHD interaction between macroscopic modes and fast particles induced by neutral beam injection in tokamaks

Authors: F. Orain, G. Brochard, T. Nicolas, H. Lütjens, X. Garbet, R. Dumont, P. Maget

Abstract: A new model of fast ion source induced by Neutral Beam Injection (NBI) in tokamaks has been implemented in the hybrid kinetic-magnetohydrodynamic code XTOR-K. This source, combined with the collisions also recently implemented, allows for the modeling of realistic slowing-down populations and their interplay with macroscopic modes. This paper describes the Neutral Beam injection and ionization mod… ▽ More A new model of fast ion source induced by Neutral Beam Injection (NBI) in tokamaks has been implemented in the hybrid kinetic-magnetohydrodynamic code XTOR-K. This source, combined with the collisions also recently implemented, allows for the modeling of realistic slowing-down populations and their interplay with macroscopic modes. This paper describes the Neutral Beam injection and ionization model designed to reproduce experimental configurations and its validation for a typical discharge of the ASDEX Upgrade tokamak. The application to the interaction of NBI-induced fast particles with a kink mode in ASDEX Upgrade demonstrates a resonance between passing particles and the n=1 mode. This resonance partially stabilizes the kink mode and induces a radial transport of fast particles. Preliminary results in ITER-like circular geometry show that NBI induces a toroidal torque but has little impact on the kink mode dynamics. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Linear stability of the ITER 15 MA scenario against the alpha fishbone

Authors: G. Brochard, R. Dumont, H. Lütjens, X. Garbet

Abstract: The stability of the $n=m=1$ alpha-fishbone kinetic-MHD mode on the ITER 15 MA baseline scenario \cite{ITERB} is analyzed using the nonlinear hybrid Kinetic-MHD code XTOR-K. Quantitative agreement is found between the complex frequencies $ω+ iγ$ computed with the linear model in \cite{Brochard2018} and XTOR-K's linear simulations. Identical precessional resonance positions in phase space are also… ▽ More The stability of the $n=m=1$ alpha-fishbone kinetic-MHD mode on the ITER 15 MA baseline scenario \cite{ITERB} is analyzed using the nonlinear hybrid Kinetic-MHD code XTOR-K. Quantitative agreement is found between the complex frequencies $ω+ iγ$ computed with the linear model in \cite{Brochard2018} and XTOR-K's linear simulations. Identical precessional resonance positions in phase space are also found between the linear model and XTOR-K. Linear hybrid simulations performed with XTOR-K on the ITER 15 MA scenario reveal that this configuration is likely to be unstable against the alpha fishbone mode. The fishbone thresholds for kinetic-MHD equilibria with flat q profiles with on-axis safety factor just below unity lies between $β_{α,thres}/β_{tot} = 6-10\%$, whereas the expected beta ratio on ITER is $β_α/β_{tot} = 15-20\%$ \cite{ITERB}. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 22 pages, 12 figures, submitted to Nuclear Fusion

Flux driven pedestal formation in tokamaks: turbulence simulations validated against the isotope effect

Authors: G. De Dominici, G. Fuhr, P. Beyer, C. Bourdelle, L. Chôné, F. Cianfrani, G. L. Falchetto, X. Garbet, Y. Sarazin

Abstract: The spontaneous pedestal formation above a power threshold at the edge of magnetically confined plasma is modelled for the first time in flux driven three-dimensional fluid simulations of electromagnetic turbulence with the code EMEDGE3D. The role of the collisional friction between trapped and passing particles is shown to be the key ingredient for shearing the radial electric field, hence stabil… ▽ More The spontaneous pedestal formation above a power threshold at the edge of magnetically confined plasma is modelled for the first time in flux driven three-dimensional fluid simulations of electromagnetic turbulence with the code EMEDGE3D. The role of the collisional friction between trapped and passing particles is shown to be the key ingredient for shearing the radial electric field, hence stabilizing the turbulence, rather than the Maxwell and Reynolds stresses. The isotope effect, observed in many tokamaks worldwide, is recovered in EMEDGE3D simulations: the power threshold for pedestal formation is lower for Tritium than for Deuterium. The turbulence auto-correlation time is found to increase with the ion mass easing the radial electric shear stabilization, hence the pedestal formation. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Comments: 7 pages, 7 figures

Influence of toroidal flow on stationary density of collisionless plasmas

Authors: Elias Laribi, Shun Ogawa, Guilhem Dif-Pradalier, Alexei Vasiliev, Xavier Garbet, Xavier Leoncini

Abstract: Starting from the given passive particle equilibrium particle cylindrical profiles, we built self-consistent stationary conditions of the Maxwell-Vlasov equation at thermodynamic equilibrium with non-flat density profiles. The solutions to the obtained equations are then discussed. It appears that the presence of an azimuthal (poloidal) flow in the plasma can insure radial confinement, while the p… ▽ More Starting from the given passive particle equilibrium particle cylindrical profiles, we built self-consistent stationary conditions of the Maxwell-Vlasov equation at thermodynamic equilibrium with non-flat density profiles. The solutions to the obtained equations are then discussed. It appears that the presence of an azimuthal (poloidal) flow in the plasma can insure radial confinement, while the presence of a longitudinal (toroidal) flow can enhance greatly the confinement. Moreover in the global physically reasonable situation, we find that no unstable point can emerge in the effective integrable Hamiltonian of the individual particles, hinting at some stability of the confinement when considering a toroidal geometry in the large aspect ratio limit. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Journal ref: Fluids 2019, 4, 172 (2019)

Tractable flux-driven temperature, density, and rotation profile evolution with the quasilinear gyrokinetic transport model QuaLiKiz

Authors: J Citrin, C Bourdelle, F J Casson, C Angioni, N Bonanomi, Y Camenen, X Garbet, L Garzotti, T Görler, O Gürcan, F Koechl, F Imbeaux, O Linder, K van de Plassche, P Strand, G Szepesi, JET Contributors

Abstract: Quasilinear turbulent transport models are a successful tool for prediction of core tokamak plasma profiles in many regimes. Their success hinges on the reproduction of local nonlinear gyrokinetic fluxes. We focus on significant progress in the quasilinear gyrokinetic transport model QuaLiKiz [C. Bourdelle et al. 2016 Plasma Phys. Control. Fusion 58 014036], which employs an approximated solution… ▽ More Quasilinear turbulent transport models are a successful tool for prediction of core tokamak plasma profiles in many regimes. Their success hinges on the reproduction of local nonlinear gyrokinetic fluxes. We focus on significant progress in the quasilinear gyrokinetic transport model QuaLiKiz [C. Bourdelle et al. 2016 Plasma Phys. Control. Fusion 58 014036], which employs an approximated solution of the mode structures to significantly speed up computation time compared to full linear gyrokinetic solvers. Optimization of the dispersion relation solution algorithm within integrated modelling applications leads to flux calculations $\times10^{6-7}$ faster than local nonlinear simulations. This allows tractable simulation of flux-driven dynamic profile evolution including all transport channels: ion and electron heat, main particles, impurities, and momentum. Furthermore, QuaLiKiz now includes the impact of rotation and temperature anisotropy induced poloidal asymmetry on heavy impurity transport, important for W-transport applications. Application within the JETTO integrated modelling code results in 1s of JET plasma simulation within 10 hours using 10 CPUs. Simultaneous predictions of core density, temperature, and toroidal rotation profiles for both JET hybrid and baseline experiments are presented, covering both ion and electron turbulence scales. The simulations are successfully compared to measured profiles, with agreement mostly in the 5-25% range according to standard figures of merit. △ Less

Submitted 7 August, 2017; v1 submitted 3 August, 2017; originally announced August 2017.

Comparison between measured and predicted turbulence frequency spectra in ITG and TEM regimes

Authors: J. Citrin, H. Arnichand, J. Bernardo, C. Bourdelle, X. Garbet, F. Jenko, S. Hacquin, M. J. Pueschel, R. Sabot

Abstract: The observation of distinct peaks in tokamak core reflectometry measurements - named quasi-coherent-modes (QCMs) - are identified as a signature of Trapped-Electron-Mode (TEM) turbulence [H. Arnichand et al. 2016 Plasma Phys. Control. Fusion 58 014037]. This phenomenon is investigated with detailed linear and nonlinear gyrokinetic simulations using the \gene code. A Tore-Supra density scan is stud… ▽ More The observation of distinct peaks in tokamak core reflectometry measurements - named quasi-coherent-modes (QCMs) - are identified as a signature of Trapped-Electron-Mode (TEM) turbulence [H. Arnichand et al. 2016 Plasma Phys. Control. Fusion 58 014037]. This phenomenon is investigated with detailed linear and nonlinear gyrokinetic simulations using the \gene code. A Tore-Supra density scan is studied, which traverses through a Linear (LOC) to Saturated (SOC) Ohmic Confinement transition. The LOC and SOC phases are both simulated separately. In the LOC phase, where QCMs are observed, TEMs are robustly predicted unstable in linear studies. In the later SOC phase, where QCMs are no longer observed, ITG modes are identified. In nonlinear simulations, in the ITG (SOC) phase, a broadband spectrum is seen. In the TEM (LOC) phase, a clear emergence of a peak at the TEM frequencies is seen. This is due to reduced nonlinear frequency broadening of the underlying linear modes in the TEM regime compared with the ITG regime. A synthetic diagnostic of the nonlinearly simulated frequency spectra reproduces the features observed in the reflectometry measurements. These results support the identification of core QCMs as an experimental marker for TEM turbulence △ Less

Submitted 12 July, 2017; originally announced July 2017.

Cross-code gyrokinetic verification and benchmark on the linear collisionless dynamics of the geodesic acoustic mode

Authors: A. Biancalani, A. Bottino, C. Ehrlacher, V. Grandgirard, G. Merlo, I. Novikau, Z. Qiu, E. Sonnendruecker, X. Garbet, T. Goerler, S. Leerink, F. Palermo, D. Zarzoso

Abstract: The linear properties of the geodesic acoustic modes (GAM) in tokamaks are investigated by means of the comparison of analytical theory and gyrokinetic numerical simulations. The dependence on the value of the safety factor, finite-orbit-width of the ions in relation to the radial mode width, magnetic-flux-surface shaping, and electron/ion mass ratio are considered. Nonuniformities in the plasma p… ▽ More The linear properties of the geodesic acoustic modes (GAM) in tokamaks are investigated by means of the comparison of analytical theory and gyrokinetic numerical simulations. The dependence on the value of the safety factor, finite-orbit-width of the ions in relation to the radial mode width, magnetic-flux-surface shaping, and electron/ion mass ratio are considered. Nonuniformities in the plasma profiles (such as density, temperature, or safety factor), electro-magnetic effects, collisions and presence of minority species are neglected. Also, only linear simulations are considered, focusing on the local dynamics. We use three different gyrokinetic codes: the lagrangian (particle-in-cell) code ORB5, the eulerian code GENE and semi-lagrangian code GYSELA. One of the main aims of this paper is to provide a detailed comparison of the numerical results and analytical theory, in the regimes where this is possible. This helps understanding better the behavior of the linear GAM dynamics in these different regimes, the behavior of the codes, which is crucial in the view of a future work where more physics is present, and the regimes of validity of each specific analytical dispersion relation. △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: Submitted to Phys. Plasmas

Tailoring steep density profile with unstable points

Authors: Shun Ogawa, Xavier Leoncini, Alexei Vasiliev, Xavier Garbet

Abstract: The mesoscopic properties of a plasma in a cylindrical magnetic field are investigated from the view point of test-particle dynamics. When the system has enough time and spatial symmetries, a Hamiltonian of a test particle is completely integrable and can be reduced to a single degree of freedom Hamiltonian for each initial state. The reduced Hamiltonian sometimes has unstable fixed points (saddle… ▽ More The mesoscopic properties of a plasma in a cylindrical magnetic field are investigated from the view point of test-particle dynamics. When the system has enough time and spatial symmetries, a Hamiltonian of a test particle is completely integrable and can be reduced to a single degree of freedom Hamiltonian for each initial state. The reduced Hamiltonian sometimes has unstable fixed points (saddle points) and associated separatrices. To choose among available dynamically compatible equilibrium states of the one particle density function of these systems we use a maximum entropy principle and discuss how the unstable fixed points affect the density profile or a local pressure gradient, and are able to create a steep profile that improves plasma confinement. △ Less

Submitted 5 April, 2018; v1 submitted 31 October, 2016; originally announced November 2016.

Comments: 6 pages, 5 figures, some content was added, some other removed in order to clarify the scope of the paper

Journal ref: Physics Letters A 383 (2019) 35-39

Study on creation and destruction of transport barriers via effective safety factors for energetic particles

Authors: Shun Ogawa, Xavier Leoncini, Guilhem Dif-Pradalier, Xavier Garbet

Abstract: Charged particles with low kinetic energy move along magnetic field lines, but so do not energetic particles. We investigate the topological structure changes in the phase space of energetic particles with respect to the magnetic one. For this study cylindrical magnetic fields with non-monotonic safety factors that induce the magnetic internal transport barrier are considered. We show that the top… ▽ More Charged particles with low kinetic energy move along magnetic field lines, but so do not energetic particles. We investigate the topological structure changes in the phase space of energetic particles with respect to the magnetic one. For this study cylindrical magnetic fields with non-monotonic safety factors that induce the magnetic internal transport barrier are considered. We show that the topological structure of the magnetic field line and of the particle trajectories can be quite different. We explain this difference using the concept of effective particle $q$-profile. Using this notion we can investigate the location and existence of resonances for particle orbits that are different from the magnetic ones. These are examined both numerically by integrating an equation of motion and theoretically by use of Alfvén's guiding center theory and by use of the effective reduced Hamiltonian for the integrable unperturbed system. It is clarified that, for the energetic particles, the grad $B$ drift effect shifts the resonances and the drift induced by curvature of the magnetic field line leads to the vanishing of the resonances. As a result, we give two different mechanisms that lead to the creation of transport barriers for energetic particles in the region where the magnetic field line is chaotic. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Comments: 8 pages, 10 figures

Journal ref: Physics of Plasmas 23, 122510 (2016)

Full particle orbit effects in regular and stochastic magnetic fields

Authors: Shun Ogawa, Benjamin Cambon, Xavier Leoncini, Michel Vittot, Diego del Castillo-Negrete, Guilhem Dif-Pradalier, Xavier Garbet

Abstract: We present a numerical study of charged particle motion in a time-independent magnetic field in cylindrical geometry. The magnetic field model consists of an unperturbed reversed-shear helical part and a perturbation consisting of a superposition of modes. Contrary to most of the previous studies, the particle trajectories are computed by directly solving the full Lorentz force equations of motion… ▽ More We present a numerical study of charged particle motion in a time-independent magnetic field in cylindrical geometry. The magnetic field model consists of an unperturbed reversed-shear helical part and a perturbation consisting of a superposition of modes. Contrary to most of the previous studies, the particle trajectories are computed by directly solving the full Lorentz force equations of motion in a six-dimensional phase space using a sixth-order, implicit, symplectic Gauss-Legendre method. The level of stochasticity in the particle orbits is diagnosed using averaged, effective Poincare sections. It is shown that when only one mode is present the particle orbits can be stochastic even though the magnetic field line orbits are not stochastic. The lack of integrability of the particle orbits in this case is related to separatrix crossing and the breakdown of the global conservation of the magnetic moment. Some perturbation consisting of two modes creates resonance overlapping, leading to Hamiltonian chaos in magnetic field lines. Then, the particle orbits exhibit a nontrivial dynamics depending on their energy and pitch angle. It is shown that the regions where the particle motion is stochastic decrease as the energy increases. The non-monotonicity of the $q$-profile implies the existence of magnetic ITBs which correspond to shearless flux surfaces located in the vicinity of the $q$-profile minimum. It is shown that depending on the energy, these magnetic ITBs might or might not confine particles. That is, magnetic ITBs act as an energy-dependent particle confinement filter. Magnetic field lines in reversed-shear configurations exhibit topological bifurcations due to separatrix reconnection. We show that a similar but more complex scenario appears in the case of particle orbits that depends in a non-trivial way on the energy and pitch angle of the particles. △ Less

Submitted 1 July, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 25 pages, accepted for publication in Phys. Plasmas

Journal ref: Phys. Plasmas 23, 072506 (2016)

Geodesic acoustic modes with poloidal mode couplings ad infinitum

Authors: Rameswar Singh, Ö D Gürcan, X Garbet, P Hennequin, L Vermare, P Morel, R Singh

Abstract: Geodesic acoustic modes (GAMs) are studied, for the first time, including all poloidal mode $(m)$ couplings using drift reduced fluid equations. The nearest neighbor coupling pattern, due to geodesic curvature, leads to a semi-infinite chain model of the GAM with the mode-mode coupling matrix elements proportional to the radial wave number $k_{r}$. The infinite chain can be reduced to a renormaliz… ▽ More Geodesic acoustic modes (GAMs) are studied, for the first time, including all poloidal mode $(m)$ couplings using drift reduced fluid equations. The nearest neighbor coupling pattern, due to geodesic curvature, leads to a semi-infinite chain model of the GAM with the mode-mode coupling matrix elements proportional to the radial wave number $k_{r}$. The infinite chain can be reduced to a renormalized bi-nodal chain with a matrix continued fractions. Convergence study of linear GAM dispersion with respect to $k_{r}$ and the $m$-spectra confirms that high m couplings become increasingly important with $k_{r}$. The radially sorted roots overlap with experimentally measured GAM frequency profile in low collisionality shots in Tore Supra thus explaining the reduced frequency of GAM in Tore Supra. △ Less

Submitted 18 December, 2015; originally announced December 2015.

PoPe (Projection on Proper elements) for code control: verification, numerical convergence and reduced models. Application to plasma turbulence simulations

Authors: T Cartier-Michaud, P Ghendrih, Y Sarazin, J Abiteboul, H Bufferand, G Dif-Pradalier, X Garbet, V Grandgirard, G Latu, C Norscini, C Passeron, P Tamain

Abstract: The Projection on Proper elements (PoPe) is a novel method of code control dedicated to 1) checking the correct implementation of models, 2) determining the convergence of numerical methods and 3) characterizing the residual errors of any given solution at very low cost. The basic idea is to establish a bijection between a simulation and a set of equations that generate it. Recovering equations is… ▽ More The Projection on Proper elements (PoPe) is a novel method of code control dedicated to 1) checking the correct implementation of models, 2) determining the convergence of numerical methods and 3) characterizing the residual errors of any given solution at very low cost. The basic idea is to establish a bijection between a simulation and a set of equations that generate it. Recovering equations is direct and relies on a statistical measure of the weight of the various operators. This method can be used in any dimensions and any regime, including chaotic ones. This method also provides a procedure to design reduced models and quantify the ratio costs to benefits. PoPe is applied to a kinetic and a fluid code of plasma turbulence. △ Less

Submitted 5 November, 2015; originally announced November 2015.

Transport Barrier generation at the interface of regions with different zonal flows dynamics

Authors: C. Norscini, Ph. Ghendrih, T. Cartier-Michaud, G. Dif-Pradalier, X. Garbet, N. Nace, Y. Sarazin, P. Tamain

Abstract: A novel and generic understanding of spontaneous generation of transport barriers and zonation regimes in turbulent self-organization is presented. It associates the barrier onset to the development of a spectral gap between large scale flows and turbulence modes leading to a zonation regime. A robust barrier builds-up at the interface of such a region and a neighboring one with reduced zonal flow… ▽ More A novel and generic understanding of spontaneous generation of transport barriers and zonation regimes in turbulent self-organization is presented. It associates the barrier onset to the development of a spectral gap between large scale flows and turbulence modes leading to a zonation regime. A robust barrier builds-up at the interface of such a region and a neighboring one with reduced zonal flow generation. This more complex and generic transition paradigm could fit the numerous and sometimes conflicting observations as in fusion plasma experiments. Barrier relaxation by bursts of turbulence regenerate the zonal flows that are eroded by viscous (collisional-like) damping. The duration of the quiescent phase between the quasi-periodic relaxation events is governed by this damping process, hence the barrier collision frequency for fusion plasmas. △ Less

Submitted 19 June, 2015; v1 submitted 9 June, 2015; originally announced June 2015.

Chaotic motion of charged particles in toroidal magnetic configurations

Authors: Benjamin Cambon, Xavier Leoncini, Michel Vittot, Rémi Dumont, Xavier Garbet

Abstract: We study the motion of a charged particle in a tokamak magnetic field and discuss its chaotic nature. Contrary to most of recent studies, we do not make any assumption on any constant of the motion and solve numerically the cyclotron gyration using Hamiltonian formalism. We take advantage of a symplectic integrator allowing us to make long-time simulations. First considering an idealized magnetic… ▽ More We study the motion of a charged particle in a tokamak magnetic field and discuss its chaotic nature. Contrary to most of recent studies, we do not make any assumption on any constant of the motion and solve numerically the cyclotron gyration using Hamiltonian formalism. We take advantage of a symplectic integrator allowing us to make long-time simulations. First considering an idealized magnetic configuration, we add a non generic perturbation corresponding to a magnetic ripple, breaking one of the invariant of the motion. Chaotic motion is then observed and opens questions about the link between chaos of magnetic field lines and chaos of particle trajectories. Second, we return to a axi-symmetric configuration and tune the safety factor (magnetic configuration) in order to recover chaotic motion. In this last setting with two constants of the motion, the presence of chaos implies that no third global constant exists, we highlight this fact by looking at variations of the first order of the magnetic moment in this chaotic setting. We are facing a mixed phase space with both regular and chaotic regions and point out the difficulties in performing a global reduction such as gyrokinetics. △ Less

Submitted 4 December, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Journal ref: Chaos, American Institute of Physics (AIP), 2014, 24 (3), pp.033101

Angular momentum transport modeling: achievements of a gyrokinetic quasi-linear approach

Authors: P. Cottier, C. Bourdelle, Y. Camenen, O. D. Gurcan, F. J. Casson, X. Garbet, P. Hennequin, T. Tala

Abstract: QuaLiKiz, a model based on a local gyrokinetic eigenvalue solver is expanded to include momentum flux modeling in addition to heat and particle fluxes. Essential for accurate momentum flux predictions, the parallel asymmetrization of the eigenfunctions is successfully recovered by an analytical fluid model. This is tested against self-consistent gyrokinetic calculations and allows for a correct pr… ▽ More QuaLiKiz, a model based on a local gyrokinetic eigenvalue solver is expanded to include momentum flux modeling in addition to heat and particle fluxes. Essential for accurate momentum flux predictions, the parallel asymmetrization of the eigenfunctions is successfully recovered by an analytical fluid model. This is tested against self-consistent gyrokinetic calculations and allows for a correct prediction of the ExB shear impact on the saturated potential amplitude by means of a mixing length rule. Hence, the effect of the ExB shear is recovered on all the transport channels including the induced residual stress. Including these additions, QuaLiKiz remains ~10 000 faster than non-linear gyrokinetic codes allowing for comparisons with experiments without resorting to high performance computing. The example is given of momentum pinch calculations in NBI modulation experiments. △ Less

Submitted 10 January, 2014; originally announced January 2014.

Journal ref: Plasma Physics and Controlled Fusion 56 (2013) 015011

Beam-Plasma Instability and Fast Particles: the Lynden-Bell Approach

Authors: Nakia Carlevaro, Duccio Fanelli, Xavier Garbet, Philippe Ghendrih, Giovanni Montani, Marco Pettini

Abstract: The beam-plasma instability, i.e. the response of the plasma bulk to the injection of supra thermal charged-particle beams, results to be appropriately characterized by a long-range interaction system. This physical system hosts a number of very interesting phenomena and, in particular, the emergence of long-lived quasi-stationary states. We characterize the self-consistent distribution functions… ▽ More The beam-plasma instability, i.e. the response of the plasma bulk to the injection of supra thermal charged-particle beams, results to be appropriately characterized by a long-range interaction system. This physical system hosts a number of very interesting phenomena and, in particular, the emergence of long-lived quasi-stationary states. We characterize the self-consistent distribution functions of such out-of-equilibrium states by means of the Lynden-Bell's theory. The prediction of this theory, based on the statistical mechanics of the Vlasov equation, are checked against the outcomes of numerical simulations of the discrete system. Moreover, a phenomenological study of the effective resonance band for the system response is also addressed. A threshold value has been found in the initial spread of beam-particle momenta. This threshold allows to discriminate between the resonant and non-resonant regimes. The analysis of the thermalization of a few percents of the beam population characterized by large initial momenta (with respect to the main part of the beam itself) is also performed and it confirms and deepens the understanding of the physical meaning of the mentioned threshold. △ Less

Submitted 9 February, 2014; v1 submitted 19 September, 2013; originally announced September 2013.

Comments: 21 pages, 14 figures

Journal ref: Plasma Phys. Contr. Fus. 56, 035013 (2014)

Asymptotic theory of magnetic island saturation in equilibria with a finite current gradient

Authors: Andrei Smolyakov, Alexander Poye, Olivier Agullo, Sadruddin Benkadda, Xavier Garbet

Abstract: A novel asymptotic matching procedure is developed to revisit the problem of magnetic island saturation in case of a finite current gradient at the rational surface. Nonlinear dispersion relation is derived for saturated magnetic island. It is shown that arbitrary normalization factors that were present in previous theories are fully specified with the asymptotic matching. A novel asymptotic matching procedure is developed to revisit the problem of magnetic island saturation in case of a finite current gradient at the rational surface. Nonlinear dispersion relation is derived for saturated magnetic island. It is shown that arbitrary normalization factors that were present in previous theories are fully specified with the asymptotic matching. △ Less

Submitted 28 February, 2013; v1 submitted 31 January, 2013; originally announced February 2013.

Comments: This paper has been withdrawn to correct comments and comparison with previous work

Nonlinear Dynamics of Magnetic Islands Imbedded in Small-Scale Turbulence

Authors: Magali Muraglia, Oliver Agullo, Sadruddin Benkadda, Xavier Garbet, P. Beyer, Abhijit Sen

Abstract: The nonlinear dynamics of magnetic tearing islands imbedded in a pressure gradient driven turbulence is investigated numerically in a reduced magnetohydrodynamic model. The study reveals regimes where the linear and nonlinear phases of the tearing instability are controlled by the properties of the pressure gradient. In these regimes, the interplay between the pressure and the magnetic flux determ… ▽ More The nonlinear dynamics of magnetic tearing islands imbedded in a pressure gradient driven turbulence is investigated numerically in a reduced magnetohydrodynamic model. The study reveals regimes where the linear and nonlinear phases of the tearing instability are controlled by the properties of the pressure gradient. In these regimes, the interplay between the pressure and the magnetic flux determines the dynamics of the saturated state. A secondary instability can occur and strongly modify the magnetic island dynamics by triggering a poloidal rotation. It is shown that the complex nonlinear interaction between the islands and turbulence is nonlocal and involves small scales. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Journal ref: Physical Review Letters 103 (2009) 145001

Effect of the curvature and the β parameter on the nonlinear dynamics of a drift tearing magnetic island

Authors: Magali Muraglia, Olivier Agullo, Masatoshi Yagi, Sadruddin Benkadda, Beyer Peter, Xavier Garbet, Sanae -I. Itoh, Kimitaka Itoh, Abhijit Sen

Abstract: We present numerical simulation studies of 2D reduced MHD equations investigating the impact of the electronic βparameter and of curvature effects on the nonlinear evolution of drift tearing islands. We observe a bifurcation phenomenon that leads to an amplification of the pressure energy, the generation of E \times B poloidal flow and a nonlinear diamagnetic drift that affects the rotation of the… ▽ More We present numerical simulation studies of 2D reduced MHD equations investigating the impact of the electronic βparameter and of curvature effects on the nonlinear evolution of drift tearing islands. We observe a bifurcation phenomenon that leads to an amplification of the pressure energy, the generation of E \times B poloidal flow and a nonlinear diamagnetic drift that affects the rotation of the magnetic island. These dynamical modifications arise due to quasilinear effects that generate a zonal flow at the onset point of the bifurcation. Our simulations show that the transition point is influenced by the βparameter such that the pressure gradient through a curvature effect strongly stabilizes the transition. Regarding the modified rotation of the island, a model for the frequency is derived in order to study its origin and the effect of the βparameter. It appears that after the transition, an E \times B poloidal flow as well as a nonlinear diamagnetic drift are generated due to an amplification of the stresses by pressure effects. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Journal ref: Nuclear Fusion 49 (2009) 055016

Microturbulence studies in RFX-mod

Authors: F Sattin, S C Guo, I Predebon, S F Liu, X Garbet, M Veranda

Abstract: Present-days Reversed Field Pinches (RFPs) are characterized by quasi-laminar magnetic configurations in their core, whose boundaries feature sharp internal transport barriers, in analogy with tokamaks and stellarators. The abatement of magnetic chaos leads to the reduction of associated particle and heat transport along wandering field lines. At the same time, the growth of steep temperature grad… ▽ More Present-days Reversed Field Pinches (RFPs) are characterized by quasi-laminar magnetic configurations in their core, whose boundaries feature sharp internal transport barriers, in analogy with tokamaks and stellarators. The abatement of magnetic chaos leads to the reduction of associated particle and heat transport along wandering field lines. At the same time, the growth of steep temperature gradients may trigger drift microinstabilities. In this work we summarize the work recently done in the RFP RFX-mod in order to assess the existence and the impact upon transport of such electrostatic and electromagnetic microinstabilities as Ion Temperature Gradient (ITG), Trapped Electron Modes (TEM) and microtearing modes. △ Less

Submitted 10 September, 2010; originally announced September 2010.

Comments: Work presented at the 2010 Varenna workshop "Theory of Fusion Plasmas". To appear in Journal of Physics Conference Series

Journal ref: Journal of Physics: Conference Series 260 (2010) 012018

Study of ITG modes in RFX-mod using TRB code

Authors: F. Sattin, X. Garbet, S. C. Guo

Abstract: We present here a study about the stability of Ion-Temperature-Gradient drift turbulence in the Quasi-Single-Helicity regime of RFX-mod Reversed Field Pinch (RFP) using the TRB fluid electrostatic turbulence code. Our results suggest that present-days RFP plasmas are marginally stable against this kind of turbulence. The onset of the instability may be envisaged for close future regimes, in the pr… ▽ More We present here a study about the stability of Ion-Temperature-Gradient drift turbulence in the Quasi-Single-Helicity regime of RFX-mod Reversed Field Pinch (RFP) using the TRB fluid electrostatic turbulence code. Our results suggest that present-days RFP plasmas are marginally stable against this kind of turbulence. The onset of the instability may be envisaged for close future regimes, in the presence of hotter plasmas with sharper internal transport barriers. △ Less

Submitted 29 July, 2010; originally announced July 2010.

Comments: Plasma Physics and Controlled Fusion, accepted

Journal ref: Plasma Physics and Controlled Fusion 52 (2010) 105002

Shear flow effects on double tearing mode global magnetic reconnection

Authors: Thibaut Voslion, Olivier Agullo, Peter Beyer, Masatoshi Yagi, Sadruddin Benkadda, Xavier Garbet, Kimitaka Itoh, Sanae-I. Itoh

Abstract: The dynamics of a global reconnection in the presence of a poloidal shear flow which is located in between magnetic islands is investigated. Different linear regimes are identified according to the value of the resistivity and the distance between the low-order resonant surfaces. It is found that the presence of a small shear flow affects and significantly delays the global reconnection processe… ▽ More The dynamics of a global reconnection in the presence of a poloidal shear flow which is located in between magnetic islands is investigated. Different linear regimes are identified according to the value of the resistivity and the distance between the low-order resonant surfaces. It is found that the presence of a small shear flow affects and significantly delays the global reconnection processes. It is shown that this delay is linked to a breaking of symmetry imposed by the existence of the shear flow and the generation of a mean poloidal flow in the resistive layers. △ Less

Submitted 3 December, 2009; v1 submitted 26 November, 2009; originally announced November 2009.

Comments: 6 pages, 10 figures, APFA/APPTC 2009 Conference

Universal Probability Distribution Function for Bursty Transport in Plasma Turbulence

Authors: I. Sandberg, S. Benkadda, X. Garbet, G. Ropokis, K. Hizanidis, D. del-Castillo-Negrete

Abstract: Bursty transport phenomena associated with convective motion present universal statistical characteristics among different physical systems. In this letter, a stochastic univariate model and the associated probability distribution function for the description of bursty transport in plasma turbulence is presented. The proposed stochastic process recovers the universal distribution of density fluc… ▽ More Bursty transport phenomena associated with convective motion present universal statistical characteristics among different physical systems. In this letter, a stochastic univariate model and the associated probability distribution function for the description of bursty transport in plasma turbulence is presented. The proposed stochastic process recovers the universal distribution of density fluctuations observed in plasma edge of several magnetic confinement devices and the remarkable scaling between their skewness $S$ and kurtosis $K$. Similar statistical characteristics of variabilities have been also observed in other physical systems that are characterized by convection such as the X-ray fluctuations emitted by the Cygnus X-1 accretion disc plasmas and the sea surface temperature fluctuations. △ Less

Submitted 24 September, 2009; v1 submitted 3 August, 2009; originally announced August 2009.

Comments: 10 pages, 5 figures

Control of Transport-barrier relaxations by Resonant Magnetic Perturbations

Authors: M. Leconte, P. Beyer, X. Garbet, S. Benkadda

Abstract: Transport-barrier relaxation oscillations in the presence of resonant magnetic perturbations are investigated using three-dimensional global fluid turbulence simulations from first principles at the edge of a tokamak. It is shown that resonant magnetic perturbations have a stabilizing effect on these relaxation oscillations and that this effect is due mainly to a modification of the pressure pro… ▽ More Transport-barrier relaxation oscillations in the presence of resonant magnetic perturbations are investigated using three-dimensional global fluid turbulence simulations from first principles at the edge of a tokamak. It is shown that resonant magnetic perturbations have a stabilizing effect on these relaxation oscillations and that this effect is due mainly to a modification of the pressure profile linked to the presence of both residual residual magnetic island chains and a stochastic layer. △ Less

Submitted 16 February, 2009; originally announced February 2009.

Comments: 4 pages

Journal ref: Phys. Rev. Lett. 102, 045006 (2009)

Turbulent particle transport in magnetized fusion plasma

Authors: Clarisse Bourdelle, Tuong G. Hoang, Xavier Garbet

Abstract: The understanding of the mechanisms responsible for particle transport is of the utmost importance for magnetized fusion plasmas. A peaked density profile is attractive to improve the fusion rate, which is proportional to the square of the density, and to self-generate a large fraction of non-inductive current required for continuous operation. Experiments in various tokamak devices (AUG, DIII-D… ▽ More The understanding of the mechanisms responsible for particle transport is of the utmost importance for magnetized fusion plasmas. A peaked density profile is attractive to improve the fusion rate, which is proportional to the square of the density, and to self-generate a large fraction of non-inductive current required for continuous operation. Experiments in various tokamak devices (AUG, DIII-D, JET, TCV, TEXT, TFTR) have indicated the existence of an anomalous inward particle pinch. Recently, such an anomalous pinch has been unambiguously identified in Tore Supra very long discharges, in absence of toroidal electric field and of central particle source, for more than 4 minutes [1]. This anomalous particle pinch is predicted by a quasilinear theory of particle transport [2], and confirmed by non-linear turbulence simulations [3] and general considerations based on the conservation of motion invariants [4]. Experimentally, the particle pinch is found to be sensitive to the magnetic field gradient in many cases [5, 6, 7, 8], to the temperature profile [5, 9] and also to the collisionality that changes the nature of the microturbulence [10, 11, 12]. The consistency of some of the observed dependences with the theoretical predictions gives us a clearer understanding of the particle pinch in tokamaks, allowing us to predict more accurately the density profile in ITER. △ Less

Submitted 22 October, 2004; originally announced October 2004.

Comments: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France)

Te/Ti effects on JET energy confinement properties

Authors: Elina Asp, Jan Weiland, Xavier Garbet, Paola Mantica, Vassili Parail, Wolfgang Suttrop

Abstract: Lately the question has been raised if a modification of the energy-confinement scaling law with respect to the electron to ion temperature ratio, Te/Ti, is required. Theoretically, like in e.g. the Weiland model, the confinement is thought to degrade with Te/Ti and studies of the hot-ion (Ti>/Te) mode seems to corroborate this. In this paper, it is shown that due to a number of effects that can… ▽ More Lately the question has been raised if a modification of the energy-confinement scaling law with respect to the electron to ion temperature ratio, Te/Ti, is required. Theoretically, like in e.g. the Weiland model, the confinement is thought to degrade with Te/Ti and studies of the hot-ion (Ti>/Te) mode seems to corroborate this. In this paper, it is shown that due to a number of effects that cancel each other out, the energy confinement time remains constant for Te/Ti>~1. The numerical study relies on a series of JET shots specifically designed to reveal an effect of Te/Ti in the hot-electron (Te>Ti) mode. A distinct effort was made to keep all current scaling-law parameters constant, including the total heating power. The effects that provide the constant confinement times have therefore nothing to do with the global properties of the plasma, but are rather due to variations in the temperature gradients which affects the transport locally. △ Less

Submitted 21 October, 2004; originally announced October 2004.

Comments: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France)

Non-linear temperature oscillations in the plasma centre on Tore Supra and their interplay with MHD

Authors: V. S. Udintsev, G. Giruzzi, F. Imbeaux, J. -F. Artaud, X. Garbet, G. Huysmans, P. Maget, J. -L. Segui, A. Becoulet, G. T. Hoang, E. Joffrin, X. Litaudon, B. Saoutic, The Tore Supra Team

Abstract: Regular oscillations of the central electron temperature have been observed by means of ECE and SXR diagnostics during non-inductively driven discharges on Tore Supra. These oscillations are sustained by LHCD, do not have a helical structure and, therefore, cannot be ascribed as MHD phenomena. The most probable explanation of this oscillating regime (O-regime) is the assumption that the plasma c… ▽ More Regular oscillations of the central electron temperature have been observed by means of ECE and SXR diagnostics during non-inductively driven discharges on Tore Supra. These oscillations are sustained by LHCD, do not have a helical structure and, therefore, cannot be ascribed as MHD phenomena. The most probable explanation of this oscillating regime (O-regime) is the assumption that the plasma current density (and, thus, the q-profile) and the electron temperature evolve as a non-linearly coupled predator-pray system. The integrated modelling code CRONOS has been used to demonstrate that the coupled heat transport and resistive diffusion equations admit solutions for the electron temperature and the current density which have a cyclic behaviour. Recent experimental results in which the O-regime co-exists with MHD modes will be presented. Because both phenomena are linked to details of the q-profile, some interplay between MHD and oscillations may occur. The localisation of magnetic islands allows to obtain an accurate picture of the q-profile in the plasma core. In some case, MHD-driven reconnection helps in maintaining a weakly inverted q-profile that is found to be, in the CRONOS simulations, a necessary condition to trigger the oscillations. △ Less

Submitted 21 October, 2004; originally announced October 2004.

Comments: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France)