Pith. sign in

REVIEW 1 major objections 6 minor 29 references

Development and validation of a local neoclassical transport module in NLT with applications to EAST-relevant impurity transport and trapped-electron-mode stability

T0 review · 1 major / 6 minor · reviewed 2026-07-14 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A local neoclassical module in NLT matches NEO multi-species fluxes and flows, and under EAST-like conditions higher Zeff can cut TEM growth while tungsten transport flips with local gradients.

desk verdict Solid NLT methods paper: dual local neo solvers with multi-species Sugama match NEO well; EAST TEM/W bits are caveated illustrations, not the load-bearing claim. read the letter →

arxiv 2607.11103 v1 pith:WCZQUYPH submitted 2026-07-13 physics.plasm-ph

Dandan Liao , Lei Ye , Guoxu Wang , Xiaotao Xiao , Chengkang Pan , Qilong Ren , Wei Zhang , Nong Xiang This is my paper · ORCID
classification physics.plasm-ph
keywords neoclassicaltransportgyrokineticsimulationmulti-speciescollisionSugamaoperatorimpuritytrapped-electronmode
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper develops and validates a local neoclassical transport module inside the semi-Lagrangian gyrokinetic code NLT for multi-species collisional plasmas. It embeds a linearized multi-species Sugama collision operator and offers two complementary solvers: an initial-value path that uses composite orbit substeps to integrate the neoclassical drive accurately along unperturbed trajectories while keeping large macroscopic time steps, and a direct steady-state linear solve that yields the stationary response without long relaxation. Against the Eulerian code NEO, the module recovers particle and heat fluxes, B-weighted parallel flows, and bootstrap current from banana through Pfirsch–Schlüter regimes for electron–ion plasmas and three-species plasmas with carbon impurities. Applied to EAST-relevant conditions, it shows that neoclassical tungsten particle flux is sensitive to local profile gradients—including a temperature-screening sign change—and that the higher effective collisionality tied to larger Zeff can reduce linear trapped-electron-mode growth rates. The result extends NLT toward realistic multi-species collisional transport while treating neoclassical and TEM calculations independently and leaving nonlinear saturation and self-consistent impurity evolution for future work.

What carries the argument

Composite substep source integration of the neoclassical drive along unperturbed characteristics (Strang-split with semi-implicit multi-species Sugama collisions), plus a complementary direct steady-state discretization of the local drift-kinetic balance solved as a linear system; together they deliver accurate orbit-integrated sources and stationary multi-species responses.

What would settle it

At the paper’s EAST r/a≈0.6 parameters, a linear TEM scan that does not reduce growth rate when Zeff rises from ~2.4 to ~3.8 at fixed profiles, or an NEO comparison where NLT particle fluxes break ambipolarity or deviate strongly across the banana-to-Pfirsch–Schlüter collisionality scan.

Watch

Extended reading notes

Core claim

With a linearized multi-species Sugama operator, NLT’s local neoclassical module—via composite-substep initial-value integration or a direct steady-state solve—reproduces NEO particle and heat fluxes, parallel flows, and bootstrap current over a broad collisionality range for electron–ion and carbon-impurity plasmas; under the considered EAST-relevant conditions, higher Zeff-linked collisionality reduces linear TEM growth rates and tungsten neoclassical transport is sensitive to local gradients.

Load-bearing premise

That independent linear TEM growth-rate drops with higher Zeff, plus local neoclassical tungsten fluxes on fixed profiles, are enough to explain the measured reduction of TEM-like density fluctuations after impurity increase—if nonlinear saturation or profile evolution dominate, that link fails.

Editorial extensions

If this is right

  • NLT can compute local multi-species neoclassical fluxes, parallel flows, and bootstrap current at NEO-level accuracy for collisionality scans and code verification.
  • Composite orbit substeps let large macroscopic time steps remain usable for fast electrons without spoiling source accuracy or the semi-implicit collision step.
  • Under EAST-like profiles, increasing tungsten-related Zeff can reduce linear TEM growth through higher effective collisionality.
  • Neoclassical tungsten particle flux can reverse from inward to outward near ρ_pol≈0.7 when main-ion temperature-gradient screening overcomes the inward pinch.
  • The same collisional framework is positioned for future coupled neoclassical–turbulent multi-species simulations inside NLT.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the linear Zeff–collisionality reduction tracks the observed fluctuation drop, ICRF-driven impurity increases may partly stabilize core density fluctuations without requiring temperature degradation.
  • Temperature-screening sign flips of tungsten flux suggest ion-temperature-gradient control could mitigate or reverse accumulation even before turbulent transport is fully coupled.
  • Intrinsic ambipolarity recovery under refined orbit substeps is a practical numerical diagnostic other semi-Lagrangian codes could adopt for electron source accuracy.
  • The dual initial-value and direct-solver path makes systematic multi-impurity and collisionality scans cheap enough to map where neoclassical and turbulent impurity channels compete.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 6 minor

Summary. The paper develops and validates a local neoclassical transport module in the semi-Lagrangian gyrokinetic code NLT for multi-species collisional plasmas. It implements the linearized multi-species Sugama collision operator and two complementary solvers: an initial-value method with a composite substep source-integration scheme for accurate neoclassical drive along unperturbed orbits at large macroscopic time steps, and a direct steady-state solver of the stationary local drift-kinetic equation. Both are benchmarked against the Eulerian code NEO for electron–ion plasmas and three-species plasmas with carbon impurities, reproducing particle and heat fluxes, B-weighted parallel flows, and bootstrap current over a broad collisionality range. The validated framework is then applied independently to EAST-relevant tungsten neoclassical transport and linear TEM stability, showing sensitivity of W fluxes to local profile gradients (including temperature-screening sign reversal) and a reduction of linear TEM growth rates with increased Zeff-associated collisionality.

Significance. If the validation holds, this is a useful and timely extension of NLT toward realistic multi-species collisional physics. Dual independent solution routes (initial-value with composite orbit substeps and direct steady-state PETSc/MUMPS solve), recovery of intrinsic ambipolarity with Norb, and quantitative agreement with NEO full Fokker–Planck and Hirshman–Sigmar operators on fluxes, flows, and bootstrap current from banana to Pfirsch–Schlüter regimes constitute a solid, falsifiable verification package. The composite substep source-integration scheme is a concrete algorithmic contribution that preserves the large-Δt advantage of semi-implicit collisions. The EAST applications are appropriately framed as mechanism-oriented illustrations rather than self-consistent predictions, and the work positions NLT for future coupled neoclassical–turbulent multi-species studies.

major comments (1)
  1. The central validation claim (Secs. 3.2–3.3, Figs. 3–6) is well supported by dual solvers, collisionality scans, ambipolarity recovery (Fig. 2), and NEO comparisons including full FP vs HS. No load-bearing technical inconsistency is evident in the local drift-kinetic ordering, source integration, or multi-species conservation. The EAST TEM/W sections (Sec. 4) are already caveated as independent applications on prescribed profiles; they do not undermine the validation core. I therefore have no major technical objections that would require reworking the central claim.
minor comments (6)
  1. Sec. 2.2 / Fig. 1: State more explicitly the practical upper bound on macroscopic Δt once Norb is scaled (global Strang-splitting and time-discretization errors), so users know when the composite scheme ceases to help.
  2. Sec. 3.2: Briefly justify the different Nμ choices (128 for initial-value vs 32 for direct solver) in the main text rather than only in the discussion of error mechanisms, for reproducibility.
  3. Sec. 4.1 / Table 2: Clarify that neff_W is inferred from Zeff under quasineutrality (already stated) and note any sensitivity of the TEM/W results to that inference or to multi-impurity mix not retained in the model.
  4. Fig. 9: The non-monotonic growth-rate upturn at high collisionality is interesting; a short remark on whether this is a dissipative TEM-like branch or a numerical/mode-identification issue would help readers interpret the experimental markers.
  5. Presentation: Fix minor typesetting issues (e.g., “NL T” spacing in the title/header, “Pfirsch-Schluter” spelling, and ensure figure axis labels and line styles are fully legible in grayscale).
  6. References: Ref. [27] is listed as “Under review”; if possible, update status or provide a stable preprint identifier for reproducibility of the EAST context.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: validation rests on external NEO benchmarks and independent dual solvers, not on self-definitional or fitted constructions.

full rationale

The paper's central claim is a local multi-species neoclassical module (linearized Sugama operator + composite-substep initial-value scheme + direct steady-state solver) that reproduces NEO particle/heat fluxes, B-weighted parallel flows, and bootstrap current for electron-ion and three-species carbon plasmas over banana-to-Pfirsch-Schlüter collisionality (Secs. 3.2–3.3, Figs. 3–6). NEO is an independent Eulerian code (Belli & Candy); the dual NLT routes (time-relaxation vs. stationary linear solve) and ambipolarity recovery with Norb provide internal cross-checks without redefining the target moments. The Sugama operator and semi-implicit collision scheme cite prior work (including author-overlapping Ref. [21]), but that citation supplies an already-implemented operator whose multi-species conservation properties are then externally verified against NEO full-FP and Hirshman-Sigmar results; it is not a uniqueness theorem or load-bearing premise that forces the fluxes. EAST applications use prescribed experimental profiles/Zeff as inputs and compute linear TEM growth rates plus local neoclassical W fluxes independently (explicitly caveated as non-self-consistent); no parameter is fitted to the observed fluctuation drop and then re-presented as a prediction. No self-definitional identities, fitted-input-as-prediction, uniqueness import, or ansatz smuggling appear in the derivation chain. The work is therefore self-contained against external benchmarks.

Assumptions & free parameters 4 free parameters · 6 assumptions · 1 invented entities

The validation claim rests on standard local neoclassical ordering and the established linearized multi-species Sugama model, discretized and solved with author-chosen grids and orbit substeps, then checked against NEO. Application claims add experimental EAST profile/Zeff inputs and the modeling choice that linear, decoupled TEM and local neo fluxes are informative. No new physical particles or forces are invented; free parameters are numerical resolutions and experimental effective impurity content, not fitted transport coefficients.

free parameters (4)
  • Norb (orbit substeps per macroscopic Δt; Norb,e ≈ 5 Norb,i)
    Chosen for source-integration accuracy vs electron transit/bounce time; convergence shown but value is numerical free choice affecting ambipolarity recovery at low collisionality.
  • Phase-space resolution (e.g. Nz×Nv∥×Nμ = 32×32×128 IV vs 32×32×32 direct)
    Hand-chosen grids and |v|≤4.2 vT domain; different Nμ between solvers justified by error mechanisms but not uniquely determined.
  • EAST effective tungsten density neff_W and Zeff (≈2.4 / 3.8)
    Inferred from prescribed Zeff under quasineutrality rather than direct W density measurement (Table 2, Sec. 4.1); sets impurity content for TEM and neo scans.
  • Carbon impurity charge fraction fI = 0.1 in three-species benchmark
    Chosen benchmark impurity level for multi-species validation, not fitted to data but free scenario parameter.
assumptions (6)
  • domain assumption Local neoclassical drift-kinetic ordering on a fixed flux surface with prescribed n,T and radial gradients; response independent of field-line label y.
    Sec. 2.1; standard local neo framework used to drop global profile evolution and restrict characteristics to the surface.
  • domain assumption Linearized multi-species Sugama collision operator conserves particles, pair momentum/energy, and self-adjointness sufficiently for neo fluxes and bootstrap current.
    Sec. 2.1 and Refs. [22,21]; central model choice vs full Fokker–Planck (NEO FP used as external check).
  • domain assumption δf decomposition about a fixed local Maxwellian with thermodynamic drive only through radial magnetic drift across ∇n,∇T.
    Eqs. (5)–(14); excludes self-consistent Er evolution and profile relaxation in the local module.
  • ad hoc to paper Strang splitting of source advection and semi-implicit collisions with composite trapezoidal orbit substeps yields accurate large-Δt neo drive.
    Sec. 2.2; numerical axiom of the initial-value scheme, supported by convergence tests but specific to this implementation.
  • domain assumption Linear TEM growth rates at r/a=0.6 with scaled collision frequency capture the collisional stabilization trend relevant to EAST fluctuation reduction.
    Sec. 4.2; assumes linear drive reduction is the right observable for the experimental TEM-like suppression narrative.
  • standard math Axisymmetric tokamak equilibrium in field-aligned coordinates with standard GA and EAST geometric parameters.
    Sec. 2.1 and Tables 1–2; conventional geometry setup.
invented entities (1)
  • Composite substep source-integration scheme (precomputed orbit-substep trapezoidal neo source increments)
    purpose: Evaluate neoclassical drive along unperturbed characteristics accurately while retaining large macroscopic time steps under semi-implicit collisions.
    Numerical construct introduced in Sec. 2.2; not a physical entity. Independent evidence is internal convergence and NEO agreement rather than external measurement.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Development and validation of a local neoclassical transport module in NLT with applications to EAST-relevant impurity transport and trapped-electron-mode stability." pith.science (2026). https://pith.science/paper/WCZQUYPH

@misc{pith2026260711103,
  author       = {Pith},
  title        = {Pith review of: Development and validation of a local neoclassical transport module in NLT with applications to EAST-relevant impurity transport and trapped-electron-mode stability},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WCZQUYPH}},
  note         = {Machine review of arXiv:2607.11103}
}
abstract

A local neoclassical transport module has been developed and validated in the semi-Lagrangian gyrokinetic code NLT for multi-species collisional plasmas. The module incorporates a linearized multi-species Sugama collision operator and provides two complementary solution strategies. In the initial-value formulation, a composite substep source-integration scheme is introduced to accurately evaluate the neoclassical drive along unperturbed particle trajectories while retaining large macroscopic time steps. A direct steady-state solver is also implemented to obtain the stationary neoclassical response without long-time relaxation. The two approaches are benchmarked against the Eulerian neoclassical code NEO for electron-ion plasmas and three-species plasmas with carbon impurities. The NLT results reproduce the NEO particle and heat fluxes, parallel flows, and bootstrap current over a broad collisionality range. As representative applications, the validated framework is applied to EAST-relevant tungsten impurity transport and core trapped-electron-mode stability. The results show that tungsten neoclassical transport is sensitive to local profile gradients, while the increased effective collisionality associated with larger \(Z_{\rm eff}\) can reduce the linear TEM growth rate under the considered EAST-relevant conditions. These developments extend NLT toward realistic multi-species collisional transport simulations.

Figures

Figures reproduced from arXiv: 2607.11103 by the authors.

Figure 1
Figure 1. Numerical convergence of the composite substep source-integration scheme at (a/vti)τ −1 ii = 3 × 10−2 . The left panel shows the convergence with Norb,e at fixed ∆t, while the right panel shows the dependence on ∆t with nearly fixed electron orbit substep size. 102 103 -0.04 -0.03 -0.02 -0.01 0 0.01 [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Dependence of the electron and ion particle fluxes on the substep number Norb,e in the composite substep source-integration scheme at (a/vti)τ −1 ii = 3×10−4 and ∆tmax = 250Ω−1 ci . The dashed horizontal line denotes the NEO reference value. Increasing Norb,e reduces the mismatch between Γi and Γe, indicating improved recovery of intrinsic ambipolarity. A more stringent test is performed in the low-collisionality re… view at source ↗
Figure 3
Figure 3. Benchmark of the initial-value relaxation solver against NEO. The left panel shows the electron and ion particle and heat fluxes, and the right panel shows the corresponding flux￾surface-averaged parallel flows. Solid lines denote NEO results, and dashed lines denote NLT results. The same benchmark is then repeated using the direct steady-state solver. Unlike the initial-value method, the direct solver obtains the s… view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: Benchmark of the direct steady-state solver against NEO. The left panel shows the electron and ion particle and heat fluxes, and the right panel shows the corresponding flux￾surface-averaged parallel flows. Solid lines denote NEO results, and dashed lines denote NLT re…
Figure 5
Figure 5. Figure 5: Collisionality dependence of the particle and heat fluxes for electrons and carbon impurities. Blue solid lines denote NEO with the full linearized Fokker-Planck operator, green dashed lines denote NEO with the Hirshman-Sigmar operator, and red dashed lines with symbol…
Figure 6
Figure 6. Figure 6: Collisionality dependence of the flux-surface-averaged parallel flows for electrons, deuterium ions, and carbon impurities, together with the total bootstrap current. The line styles are the same as in [PITH_FULL_IMAGE:figures/full_fig_p013_6.png]
Figure 7
Figure 7. Figure 7: Linear TEM spectrum at r/a = 0.6 for impurity-free and W-containing plasmas. The left panel shows the linear growth rate γ, and the right panel shows the real frequency ω. The higher-Zeff tungsten-containing case gives the strongest reduction of the TEM growth rate. Th…
Figure 8
Figure 8. Figure 8: Collisionality scaling of the TEM growth rate for the W, Zeff = 3.8 case at r/a = 0.6. The collision frequency is scaled from the collisionless limit to 0.1νc, 0.5νc, and the experimen￾tal value νc. The decreasing growth rate with increasing collision frequency demonst…
Figure 9
Figure 9. Figure 9: Dependence of the linear growth rate and real frequency on the effective collision frequency for tungsten-containing plasmas with Zeff = 2.4 and Zeff = 3.8. The dashed horizontal lines denote the collisionless reference values, and the green markers indicate the experi…
Figure 10
Figure 10. Figure 10: Radial profile of the neoclassical tungsten particle flux normalized by the gyro￾Bohm particle flux for W, Zeff = 2.4 and W, Zeff = 3.8. The dashed horizontal line denotes the zero-flux level. Negative values correspond to inward tungsten transport, while positive val…
Figure 11
Figure 11. Figure 11: Dependence of the neoclassical tungsten particle flux on the main-ion temperature gradient a/LTi at Zeff = 3.8 and ρpol = 0.7. The green marker denotes the experimental value. Negative values correspond to inward tungsten transport, while positive values indicate outw…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references

  1. [1]

    F. L. Hinton and S. K. Wong. Neoclassical ion transport in rotating axisymmetric plasmas. The Physics of Fluids , 28(10):3082–3098, 1985

    1985

  2. [2]

    C. N. Lashmore-Davies. Collisional transport in magnet ized plasmas, by per helander and dieter j. sigmar, cambridge monograph on plasma physics , cambridge university press (2002). Journal of Plasma Physics , 70(4):501–502, 2004

    2002

  3. [3]

    S. P. Hirshman and D. J. Sigmar. Neoclassical transport o f impurities in tokamak plasmas. Nuclear Fusion, 21(9):1079, 1981

    1981

  4. [4]

    The transport and radiation loss of tungsten impurities in c ore plasma for the CFEDR conventional H-mode scenario

    Chengxi Zhou, Ling Zhang, Jiale Chen, Xiang Jian, Rui Din g, and Vincent Chan. The transport and radiation loss of tungsten impurities in c ore plasma for the CFEDR conventional H-mode scenario. Plasma Science and Technology , 27(10):104008, 2025

    2025

  5. [5]

    F. L. Hinton and M. N. Rosenbluth. Dynamics of axisymmetr ic and poloidal flows in tokamaks. Plasma Physics and Controlled Fusion , 41(3A):A653, 1999

    1999

  6. [6]

    M. N. Rosenbluth, R. D. Hazeltine, and F. L. Hinton. Plasm a Transport in Toroidal Confinement Systems. The Physics of Fluids , 15(1):116–140, 1972

    1972

  7. [7]

    E. A. Belli and J. Candy. Kinetic calculation of neoclass ical transport including self- consistent electron and impurity dynamics. Plasma Physics and Controlled Fusion , 50(9):095010, 2008

    2008

  8. [8]

    E. A. Belli and J. Candy. Full linearized Fokker–Planck c ollisions in neoclassical transport simulations. Plasma Physics and Controlled Fusion , 54(1):015015, 2012

    2012

Show all 29 references
  1. [9]

    E. A. Belli and J. Candy. Neoclassical transport in toroi dal plasmas with nonaxisym- metric flux surfaces. Plasma Physics and Controlled Fusion , 57(5):054012, 2015

    2015

  2. [10]

    Candy, E

    J. Candy, E. A. Belli, and R. V. Bravenec. A high-accurac y Eulerian gyrokinetic solver for collisional plasmas. Journal of Computational Physics , 324:73–93, 2016

    2016

  3. [11]

    Robert Hager, E. S. Yoon, S. Ku, E. F. D’Azevedo, P. H. Wor ley, and C. S. Chang. A fully non-linear multi-species Fokker–Planck–Landau col lision operator for simulation of fusion plasma. Journal of Computational Physics , 315:644–660, 2016

    2016

  4. [12]

    Robert Hager and C. S. Chang. Gyrokinetic neoclassical study of the bootstrap cur- rent in the tokamak edge pedestal with fully non-linear Coul omb collisions. Physics of Plasmas , 23(4):042503, 2016

    2016

  5. [13]

    Dominski, C

    J. Dominski, C. S. Chang, R. Hager, S. Ku, E. S. Yoon, and V . Parail. Neoclassical transport of tungsten ion bundles in total-f neoclassical g yrokinetic simulations of a whole-volume JET-like plasma. Physics of Plasmas , 31(3):032303, 2024

    2024

  6. [14]

    Vernay, S

    T. Vernay, S. Brunner, L. Villard, B. F. McMillan, S. Jol liet, T. M. Tran, A. Bottino, and J. P. Graves. Neoclassical equilibria as starting point for global gyrokinetic microturbulence simulations. Physics of Plasmas , 17(12):122301, 2010

    2010

  7. [15]

    Vernay, S

    T. Vernay, S. Brunner, L. Villard, B. F. McMillan, S. Jol liet, T. M. Tran, and A. Bot- tino. Synergy between ion temperature gradient turbulence and neoclassical processes in global gyrokinetic particle-in-cell simulations. Physics of Plasmas , 19(4):042301, 2012. 19

    2012

  8. [16]

    Lanti, N

    E. Lanti, N. Ohana, N. Tronko, T. Hayward-Schneider, A. Bottino, B. F. McMillan, A. Mishchenko, A. Scheinberg, A. Biancalani, P. Angelino, S . Brunner, J. Dominski, P. Donnel, C. Gheller, R. Hatzky, A. Jocksch, S. Jolliet, Z. X . Lu, J. P. Martin Collar, I. Novikau, E. Sonnend...

    2020

  9. [17]

    Idomura, K

    Y. Idomura, K. Obrejan, Y. Asahi, and M. Honda. Dynamics of enhanced neoclas- sical particle transport of tracer impurity ions in ion temp erature gradient driven turbulence. Physics of Plasmas , 28(1):012501, 2021

    2021

  10. [18]

    Donnel, X

    P. Donnel, X. Garbet, Y. Sarazin, V. Grandgirard, Y. Asa hi, N. Bouzat, E. Caschera, G. Dif-Pradalier, C. Ehrlacher, P. Ghendrih, C. Gillot, G. L atu, and C. Passeron. A multi-species collisional operator for full-f global gyr okinetics codes: Numerical aspects and verification...

    2019

  11. [19]

    K. Lim, X. Garbet, Y. Sarazin, E. Gravier, M. Lesur, G. Lo -Cascio, and T. Rouyer. Self-consistent gyrokinetic modeling of turbulent and neo classical tungsten transport in toroidally rotating plasmas. Physics of Plasmas , 30(8):082501, 2023

    2023

  12. [20]

    Im- plementation of field-aligned coordinates in a semi-Lagran gian gyrokinetic code for tokamak turbulence simulation

    Lei YE, Xiaotao XIAO, Yingfeng XU, Zongliang DAI, and Sh aojie W ANG. Im- plementation of field-aligned coordinates in a semi-Lagran gian gyrokinetic code for tokamak turbulence simulation. Plasma Science and Technology , 20(7):074008, jun 2018

    2018

  13. [21]

    Semi-implicit scheme for multi-species collision operators in tokamak pl asma simulations

    Dandan Liao, Lei Ye, Pengfei Zhao, Qilong Ren, and Nong X iang. Semi-implicit scheme for multi-species collision operators in tokamak pl asma simulations. Com- puter Physics Communications , 320:110012, 2026

    2026

  14. [22]

    Sugama, T

    H. Sugama, T. H. Watanabe, and M. Nunami. Linearized mod el collision operators for multiple ion species plasmas and gyrokinetic entropy ba lance equations. Physics of Plasmas , 16(11):112503, 2009

    2009

  15. [23]

    On the construction and comparison of d ifference schemes

    Gilbert Strang. On the construction and comparison of d ifference schemes. SIAM Journal on Numerical Analysis , 5(3):506–517, 1968

    1968

  16. [24]

    R. E. Waltz, G. M. Staebler, W. Dorland, G. W. Hammett, M. Kotschenreuther, and J. A. Konings. A gyro-Landau-fluid transport model. Physics of Plasmas , 4(7):2482– 2496, 1997

    1997

  17. [25]

    Gyrokinetic Simulation of Microtearing Turbulence

    Hauke Doerk. Gyrokinetic Simulation of Microtearing Turbulence . PhD thesis, Uni- versitaet Ulm, 2012

    2012

  18. [26]

    Crandall, D

    P. Crandall, D. Jarema, H. Doerk, Q. Pan, G. Merlo, T. Gör ler, A. Bañón Navarro, D. Told, M. Maurer, and F. Jenko. Multi-species collisions f or delta-f gyrokinetic simulations: Implementation and verification with gene. Computer Physics Commu- nications, 255:107360, 2020

    2020

  19. [27]

    X Wang, L Ye, P Li, F Chen, H

    W Zhang, G. X Wang, L Ye, P Li, F Chen, H. W Yang, X. M Zhai, L. N Liu, H Yang, T Jin, L. Q Xu, D. D Liao, N Gorelenkov, S. Mazzi, J. P Qia n, Y. T Song, and the EAST team. High-Z impurity polarization shielding o f short-wavelength trapped-electron-mode turbulence in EAST. Un...

    2026

  20. [28]

    Manas, Y

    P. Manas, Y. Camenen, S. Benkadda, W. A. Hornsby, and A. G . Peeters. Enhanced stabilisation of trapped electron modes by collisional ene rgy scattering in tokamaks. Physics of Plasmas , 22(6):062302, 2015

    2015

  21. [29]

    Application of Lorentz collision operator in the gyrokinetic code NLT

    Shiqi Xiao, Haochuan Wang, Zihao Wang, and Shaojie Wang . Application of Lorentz collision operator in the gyrokinetic code NLT. Plasma Science and Technology , 27(11):115102, 2025. 21

    2025

Pith tools

Reviewed July 14, 2026 · model on record in the stance chip above.