BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

P. Ricci

P. Ricci appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phnucl-thastro-ph.SRhep-phnucl-exastro-ph.HEnlin.CDphysics.comp-ph

Catalog footprint

What is connected

20works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

20 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
preprint2015arXiv

SMWDs as SGRs/AXPs and the lepton number violation

Possible nature of strongly magnetized white dwarfs (SMWDs) is studied. It is shown that for relatively low values of the equatorial surface magnetic field $B\,\sim\,10^9\,-\,10^{11}$ G they can be good candidates for soft gamma-ray repeaters and anomalous X-ray pulsars (SGRs/AXPs). For the case of iron SMWDs the influence of a neutrinoless electron to positron conversion on the SGRs/AXPs luminosity is estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical analysis and modeling of intermittent transport events in the tokamak SOL

The turbulence observed in the scrape-off-layer of a tokamak is often characterized by intermittent events of bursty nature, a feature which raises concerns about the prediction of heat loads on the physical boundaries of the device. It appears thus necessary to delve into the statistical properties of turbulent physical fields such as density, electrostatic potential and temperature, focusing on the mathematical expression of tails of the probability distribution functions. The method followed here is to generate statistical information from time-traces of the plasma density stemming from Braginskii-type fluid simulations, and check this against a first-principles theoretical model. The analysis of the numerical simulations indicates that the probability distribution function of the intermittent process contains strong exponential tails, as predicted by the analytical theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Calculation of Doublet Capture Rate for Muon Capture in Deuterium within Chiral Effective Field Theory

The doublet capture rate of the negative muon capture in deuterium is calculated employing the nuclear wave functions generated from accurate nucleon-nucleon potentials constructed at next-to-next-to-next-to-leading order of heavy-baryon chiral perturbation theory and the weak meson exchange current operator derived within the same formalism. All but one of the low-energy constants that enter the calculation were fixed from pion-nucleon and nucleon-nucleon scattering data. The low-energy constant d^R (c_D), which cannot be determined from the purely two-nucleon data, was extracted recently from the triton beta-decay and the binding energies of the three-nucleon systems. The calculated values of the doublet capture rates show a rather large spread for the used values of the d^R. Precise measurement of the doublet capture rate in the future will not only help to constrain the value of d^R, but also provide a highly nontrivial test of the nuclear chiral EFT framework. Besides, the precise knowledge of the constant d^R will allow for consistent calculations of other two-nucleon weak processes, such as proton-proton fusion and solar neutrino scattering on deuterons, which are important for astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nondiffusive suprathermal ion transport in simple magnetized toroidal plasmas

We investigate suprathermal ion dynamics in simple magnetized toroidal plasmas in the pres- ence of electrostatic turbulence driven by the ideal interchange instability. Turbulent fields from fluid simulations are used in the non-relativistic equation of ion motion to compute suprathermal tracer ion trajectories. Suprathermal ion dispersion starts with a brief ballistic phase, during which particles do not interact with the plasma, followed by a turbulence interaction phase. In this one simple system, we observe the entire spectrum of suprathermal ion dynamics, from subdiffusion to superdiffusion, depending on beam energy and turbulence amplitude. We estimate the duration of the ballistic phase and identify basic mechanisms during the interaction phase that determine the character of suprathermal ion dispersion upon the beam energy and turbulence fluctuation amplitude.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Muon Capture in Deuterium

Model dependence of the capture rates of the negative muon capture in deuterium is studied starting from potential models and the weak two-body meson exchange currents constructed in the tree approximation and also from an effective field theory. The tree one-boson exchange currents are derived from the hard pion chiral Lagrangians of the $N Δπρωa_1$ system. If constructed in conjunction with the one-boson exchange potentials, the capture rates can be calculated consistently. On the other hand, the effective field theory currents, constructed within the heavy baryon chiral perturbation theory, contain a low energy constant $\hat d ^R$ that cannot be extracted from data at the one-particle level nor determined from the first principles. Comparative analysis of the results for the doublet transition rate allows us to extract the constant $\hat d ^R$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Linearized model Fokker-Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests

A set of key properties for an ideal dissipation scheme in gyrokinetic simulations is proposed, and implementation of a model collision operator satisfying these properties is described. This operator is based on the exact linearized test-particle collision operator, with approximations to the field-particle terms that preserve conservation laws and an H-Theorem. It includes energy diffusion, pitch-angle scattering, and finite Larmor radius effects corresponding to classical (real-space) diffusion. The numerical implementation in the continuum gyrokinetic code GS2 is fully implicit and guarantees exact satisfaction of conservation properties. Numerical results are presented showing that the correct physics is captured over the entire range of collisionalities, from the collisionless to the strongly collisional regimes, without recourse to artificial dissipation.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Interactions of the solar neutrinos with the deuterons

Starting from chiral Lagrangians, possessing the SU(2)_L x SU(2)_R local chiral symmetry, we derive weak axial one-boson exchange currents in the leading order in the 1/M expansion (M is the nucleon mass). We apply these currents in calculations of the cross sections for the disintegration of the deuterons by the low energy neutrinos. The nuclear wave functions are derived from a variant of the OBEPQB potential and from the Nijmegen 93 and Nijmegen I nucleon-nucleon interactions. The comparison of our cross sections with those obtained within the pionless effective field theory and other potential model calculations shows that the solar neutrino-deuteron cross sections can be calculated within an accuracy of 3.3 %.

0 citations0 reviewsNo code linkedOpen access