Engineering topics

Dif-Pradalier, G.

Publications and source records attributed to Dif-Pradalier, G..

14 records

DOE OSTI2022

Operating a full tungsten actively cooled tokamak: overview of WEST first phase of operation

WEST is an MA class superconducting, actively cooled, full tungsten (W) tokamak, designed to operate in long pulses up to 1000 s. In support of ITER operation and DEMO conceptual activities, key missions of WEST are: (i) qualification of high heat flux plasma-facing components in integrating both technological and physics aspects in relevant heat and particle exhaust conditions, particularly for the tungsten monoblocks foreseen in ITER divertor; (ii) integrated steady-state operation at high confinement, with a focus on power exhaust issues. During the phase 1 of operation (2017–2020), a set of actively cooled ITER-grade plasma facing unit prototypes was integrated into the inertially cooled W coated startup lower divertor. Up to 8.8 MW of RF power has been coupled to the plasma and divertor heat flux of up to 6 MW m –2 were reached. Long pulse operation was started, using the upper actively cooled divertor, with a discharge of about 1 min achieved. This paper gives an overview of the results achieved in phase 1. Perspectives for phase 2, operating with the full capability of the device with the complete ITER-grade actively cooled lower divertor, are also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Towards an Emerging Understanding of Non-local Transport

In this paper, recent progress on experimental analysis and theoretical models for non-local transport (non-Fickian fluxes in real space) is reviewed. The non-locality in the heat and momentum transport observed in the plasma, the departures from linear flux-gradient proportionality, and externally triggered non-local transport phenomena are described in both L-mode and improved-mode plasmas. Ongoing evaluation of ‘fast front’ and ‘intrinsically non-local’ models, and their success in comparisons with experimental data, are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Edge localized mode rotation and the nonlinear dynamics of filaments

Edge Localized Modes (ELMs) rotating precursors were reported few milliseconds before an ELM crash in several tokamak experiments. Also, the reversal of the filaments rotation at the ELM crash is commonly observed. In this article, we present a mathematical model that reproduces the rotation of the ELM precursors as well as the reversal of the filaments rotation at the ELM crash. Linear ballooning theory is used to establish a formula estimating the rotation velocity of ELM precursors. The linear study together with nonlinear magnetohydrodynamic simulations give an explanation to the rotations observed experimentally. Unstable ballooning modes, localized at the pedestal, grow and rotate in the electron diamagnetic direction in the laboratory reference frame. Approaching the ELM crash, this rotation decreases corresponding to the moment when the magnetic reconnection occurs. During the highly nonlinear ELM crash, the ELM filaments are cut from the main plasma due to the strong sheared mean flow that is nonlinearly generated via the Maxwell stress tensor.

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

L to H mode transition: Parametric dependencies of the temperature threshold

The L to H mode transition occurs at a critical power which depends on various parameters, such as the magnetic field, the density, etc. Experimental evidence on various tokamaks (JET, ASDEX-Upgrade, DIII-D, Alcator C-Mod) points towards the existence of a critical temperature characterizing the transition. This criterion for the L-H transition is local and is therefore easier to be compared to theoretical approaches. In order to shed light on the mechanisms of the transition, simple theoretical ideas are used to derive a temperature threshold (T th ). They are based on the stabilization of the underlying turbulence by a mean radial electric field shear. The nature of the turbulence varies as the collisionality decreases, from resistive ballooning modes to ion temperature gradient and trapped electron modes. The obtained parametric dependencies of the derived T th are tested versus magnetic field, density, effective charge. Furthermore, various robust experimental observations are reproduced, in particular T th increases with magnetic field B and increases with density below the density roll-over observed on the power threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Towards an emerging understanding of non-locality phenomena and non-local transport

Here, recent progress on experimental analysis and theoretical models for non-local transport (non-Fickian fluxes in real space) is reviewed. The non-locality in the heat and momentum transport observed in the plasma, the departures from linear flux-gradient proportionality, and externally triggered non-local transport phenomena are described in both L-mode and improved-mode plasmas. Ongoing evaluation of 'fast front' and 'intrinsically non-local' models, and their success in comparisons with experimental data, are reported on.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

E × B shear pattern formation by radial propagation of heat flux waves

A novel theory to describe the formation of E×B flow patterns by radially propagating heat flux waves is presented. A model for heat avalanche dynamics is extended to include a finite delay time between the instantaneous heat flux and the mean flux, based on an analogy between heat avalanche dynamics and traffic flow dynamics. The response time introduced here is an analogue of the drivers' response time in traffic dynamics. The microscopic foundation for the time delay is the time for mixing of the phase space density. The inclusion of the finite response time changes the model equation for avalanche dynamics from Burgers equation to a nonlinear telegraph equation. Based on the telegraph equation, the formation of heat flux jams is predicted. The growth rate and typical interval of jams are calculated. The connection of the jam interval to the typical step size of the E×B staircase is discussed.

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 OSTI2013

Transport of radial heat flux and second sound in fusion plasmas

Simple flux-gradient relations that involve time delay and radial coupling are discussed. Such a formulation leads to a rather simple description of avalanches and may explain breaking of gyroBohm transport scaling. The generalization of the flux-gradient relation (i.e., constitutive relation), which involve both time delay and spatial coupling, is derived from drift-kinetic equation, leading to kinetic definitions of constitutive elements such as the flux of radial heat flux. This allows numerical simulations to compute these cubic quantities directly. The formulation introduced here can be viewed as an extension of turbulence spreading to include the effect of spreading of cross-phase as well as turbulence intensity, combined in such a way to give the flux. The link between turbulence spreading and entropy production is highlighted. An extension of this formulation to general quasi-linear theory for the distribution function in the phase space of radial position and parallel velocity is also discussed.

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↗