BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

P. Ricci

P. Ricci contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.plasm-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

Development of Advanced Linearized Gyrokinetic Collision Operators Using a Moment Approach

The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge R. et al., J. Plasma Phys. 85, 905850604 (2019)) and of the widely-used linearized GK Sugama collision operator (Sugama H. et al., Phys. Plasmas 16, 112503 (2009)) is reported. An approach based on a Hermite-Laguerre moment expansion of the perturbed gyrocenter distribution function is used, referred to as gyro-moment expansion. This approach allows considering arbitrary perpendicular wavenumber and expressing the two linearized GK operators as a linear combination of gyro-moments where the expansion coefficients are given by closed analytical expressions that depend on the perpendicular wavenumber and on the temperature and mass ratios of the colliding species. The drift-kinetic (DK) limits of the GK linearized Coulomb and Sugama operators are also obtained. Comparisons between the gyro-moment approach with the GK continuum code GENE are reported focusing on the ion-temperature-gradient (ITG) instability and zonal flow (ZF) damping, finding an excellent agreement. In particular, we demonstrate that the GK linearized Sugama yields a stronger collisional damping of the ZF residual compared to the GK linearized Coulomb. Finally, we show that the numerical efficiency of the gyro-moment approach increases with collisionality, a desired property for boundary plasma applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First-principles density limit scaling in tokamaks based on edge turbulent transport and implications for ITER

A first-principles scaling law, based on turbulent transport considerations, and a multi-machine database of density limit discharges from the ASDEX Upgrade, JET and TCV tokamaks, show that the increase of the boundary turbulent transport with the plasma collisionality sets the maximum density achievable in tokamaks. This scaling law shows a strong dependence on the heating power, therefore predicting for ITER a significantly larger safety margin than the Greenwald empirical scaling (Greenwald et al, Nucl. Fusion, 28(12), 1988) in case of unintentional H-L transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Global fluid simulation of plasma turbulence in a stellarator with an island divertor

Results of a three-dimensional, flux-driven, electrostatic, global, two-fluid turbulence simulation for a 5-field period stellarator with an island divertor are presented. The numerical simulation is carried out with the GBS code, recently extended to simulate plasma turbulence in non-axisymmetric magnetic equilibria. The vacuum magnetic field used in the simulation is generated with the theory of Dommaschk potentials, and describes a configuration with a central region of nested flux surfaces, surrounded by a chain of magnetic islands, similarly to the diverted configurations of W7-X. The heat outflowing from the core reaches the island region and is transported along the magnetic islands, striking the vessel walls, which correspond to the boundary of the simulation domain. The radial transport of particles and heat is found to be mainly driven by a field-aligned coherent mode with poloidal number $m=4$. The analysis of this mode, based on non-local linear theory considerations, shows its ballooning nature. In contrast to tokamak simulations and experiments, where blobs often contribute to transport, we do not observe the presence of intermittent transport events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local Gyrokinetic Collisional Theory of the Ion-Temperature Gradient Mode

We present a study of the linear properties of ion temperature gradient (ITG) modes with collisions modelled by the linearized gyrokinetic (GK) Coulomb collision operator (Frei et al. 2021) in the local limit. The study is based on a Hermite-Laguerre polynomial expansion of the perturbed ion distribution function applied to the linearized GK Boltzmann equation, yielding a hierarchy of coupled equations for the expansion coefficients, referred to as gyro-moments. We explore analytically the collisionless and high-collisional limits of the gyro-moment hierarchy. Parameter scans revealing the dependence of the ITG growth rate on the collisionality are reported, showing strong damping at small scales as the collisionality increases. These properties are compared with the predictions based on the Sugama, the momentum-conserving pitch-angle scattering, the Hirshman- Sigmar-Clarke, and the Daugherty collision operators. The importance of finite Larmor radius (FLR) terms in the collision operators is pointed out by the appearance of a short wavelength (SW) ITG branch when collisional FLR terms are neglected, this branch being completely suppressed by collisional FLR effects. We demonstrate that energy diffusion is important at high collisionality and small scale lengths and that, among the collision operators considered in this work, the GK Sugama collision operator yields, in general, the smallest deviation on the ITG growth rate compared to the GK Coulomb collision operator. Convergence studies of the gyro-moment method are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent transport regimes in the tokamak boundary and operational limits

Two-fluid, three-dimensional, flux-driven, global, electromagnetic turbulence simulations carried out by using the GBS code are used to identify the main parameters controlling turbulent transport in the tokamak boundary and to delineate an electromagnetic phase space of edge turbulence. Four turbulent transport regimes are identified: (i) a regime of fully developed turbulence appearing at intermediate values of collisionality and $β$, with turbulence driven by resistive ballooning modes, related to the L-mode operation of tokamaks, (ii) a regime of reduced turbulent transport at low collisionality and large heat source, with turbulence driven by drift-waves, related to a high-density H-mode regime, (iii) a regime of extremely large turbulent transport at high collisionality, which is associated with the crossing of the density limit, and (iv) a regime above the ideal ballooning limit at high $β$, with global modes affecting the dynamics of the entire confined region, which can be associated with the crossing of the $β$ limit. The transition from the reduced to the developed turbulent transport regime is associated here with the H-mode density limit and an analytical scaling law for maximum edge density achievable in H-mode is obtained. Analogously, analytical scaling laws for the crossing of the L-mode density and $β$ limits are provided and compared to the results of GBS simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A self-consistent multi-component model of plasma turbulence and kinetic neutral dynamics for the simulation of the tokamak boundary

A self-consistent model is presented for the simulation of a multi-component plasma in the tokamak boundary. A deuterium plasma is considered, with the plasma species that include electrons, deuterium atomic ions and deuterium molecular ions, while the deuterium atoms and molecules constitute the neutral species. The plasma and neutral models are coupled via a number of collisional interactions, which include dissociation, ionization, charge-exchange and recombination processes. The derivation of the three-fluid drift-reduced Braginskii equations used to describe the turbulent plasma dynamics is presented, including its boundary conditions. The kinetic advection equations for the neutral species are also derived, and their numerical implementation discussed. The first results of multi-component plasma simulations carried out by using the GBS code are then presented and analyzed, being compared with results obtained with the single-component plasma model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The GBS code for the self-consistent simulation of plasma turbulence and kinetic neutral dynamics in the tokamak boundary

A new version of GBS (Ricci et al. Plasma Phys. Control. Fusion 54, 124047, 2012; Halpern et al. J. Comput. Phys. 315, 388-408, 2016; Paruta et al. Phys. Plasmas 25, 112301, 2018) is described. GBS is a three-dimensional, flux-driven, global, two-fluid turbulence code developed for the self-consistent simulation of plasma turbulence and kinetic neutral dynamics in the tokamak boundary. In the new version presented here, the simulation domain is extended to encompass the whole plasma volume, avoiding an artificial boundary with the core, hence retaining the core-edge-SOL interplay. A new toroidal coordinate system is introduced to increase the code flexibility, allowing for the simulation of arbitrary magnetic configurations (e.g. single-null, double-null and snowflake configurations), which can also be the result of the equilibrium reconstruction of an experimental discharge. The implementation of a new iterative solver for the Poisson and Ampère equations is presented, leading to a remarkable speed-up of the code with respect to the use of direct solvers, therefore allowing for efficient electromagnetic simulations that avoid the use of the Boussinesq approximation. The self-consistent kinetic neutral model, initially developed for limited configurations, is ported to the magnetic configurations considered by the present version of GBS and carefully optimized. A new MPI parallelisation is implemented to evolve the plasma and neutral models in parallel, thus improving the code scalability. The numerical implementation of the plasma and neutral models is verified by means of the method of manufactured solutions. As an example of the simulation capabilities of the new version of GBS, a simulation of a TCV tokamak discharge is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and SOL of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted Ohmic L-mode discharges on the TCV tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from 5 diagnostic systems -- giving a total of 45 1- and 2-D comparison observables in two toroidal magnetic field directions. The dataset is used to validate three flux-driven 3D fluid-turbulence models: GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP -- both in terms of profile shape and absolute magnitude -- while a poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A gyrokinetic model for the plasma periphery of tokamak devices

A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts, assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behavior of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function and includes a non-linear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a velocity-space Hermite-Laguerre polynomial basis and the gyrokinetic equation is projected onto the same basis, obtaining the spatial and temporal evolution of the Hermite-Laguerre expansion coefficients. The Hermite-Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwell's equations derived from a variational principle, with the velocity integrals of the gyroaveraged distribution function explicitly evaluated.

0 citations0 reviewsNo code linkedOpen access