Pith. sign in

REVIEW 5 major objections 6 minor 1 cited by

Spectrally Accelerated Edge and Scrape-Off Layer Gyrokinetic Turbulence Simulations

T0 review · 5 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read A velocity-space trick speeds gyrokinetic edge turbulence runs 50x

desk verdict A solid numerical-methods paper: the Hermite-Laguerre velocity-space implementation in GENE-X gives a credible ~50x speed-up for TCV-X21 and reproduces grid profiles, but the single-scale-temperature basis limits how far that speed-up transfers to steep H-mode or ITER cases. read the letter →

arxiv 2411.09232 v1 pith:P2ZQGZBH submitted 2024-11-14 physics.plasm-ph

B. J. Frei , P. Ulbl , J. Trilaksono , F. Jenko This is my paper · ORCID
classification physics.plasm-ph PACS 52.65.-y52.35.Ra52.25.Fi
keywords gyrokineticsspectralvelocity-spaceexpansionHermite-Laguerrebasisedgeandscrape-offlayerturbulencefull-fgyrokineticcodetrappedelectronmodescomputationalspeed-upmethodofmanufacturedsolutions
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

The paper presents the first full-f gyrokinetic simulations of edge and scrape-off layer turbulence that accelerate the velocity-space discretization with a spectral (Hermite-Laguerre) expansion, implemented in the GENE-X code. The central claim is that this spectral formulation reproduces the outboard midplane profiles—density, electron and ion temperatures, and radial electric field—of the TCV-X21 reference case with excellent agreement against previously validated grid simulations, while requiring vastly fewer velocity-space degrees of freedom. The quantitative payoff claimed is a speed-up of approximately 50 for TCV-X21, reducing the cost from millions of CPU-hours and weeks of wall time to about 0.05–0.1 MCPUh and one to two days on CPU-based supercomputers. If correct, this would make routine high-fidelity gyrokinetic edge/SOL studies feasible for medium-sized devices and move reactor-relevant devices such as ITER within reach of this modeling approach.

What carries the argument

The central object is the spectral expansion of the full-f distribution function onto scaled Hermite-Laguerre polynomials, $F_\alpha = \sum_{p,j} N^{pj}_\alpha \hat H_p(\hat v_{\|\alpha}) L_j(\hat\mu_\alpha) F_{M\alpha}$, with a constant per-species scaled temperature $\tau_\alpha$ chosen a priori ($\tau_e = 114$ eV, $\tau_i = 102.5$ eV for TCV-X21). The Hermite-Laguerre basis carries the argument: it converts the 2D velocity-space grid into a small set of coupled spectral coefficients, with the coupling expressed through sparse recurrence identities (Eq. 15). The spectral formulation yields closed expressions for the GK Vlasov equation (Eq. 12), the quasineutrality, Ampere, and Ohm equations (Eq. 16), and the LBD collision operator (Eq. 17), with the fluid moments of interest expressed directly in terms of low-order coefficients (Eq. 21). The spectral approach also exactly conserves the total energy of the GK system, and a Landau-damping-like diagonal dissipation term is added at the highest orders to stabilize the truncation.

What would settle it

A gyrokinetic simulation of a different edge/SOL scenario with a wider temperature range or a pronounced H-mode pedestal, run with the spectral approach at the claimed low resolution (e.g., $N_{v\|} \approx 6$, $N_\mu \approx 4$) and with the single-constant $\tau_\alpha$ recipe, would falsify the generality claim if its OMP profiles deviated systematically from grid or experimental results and if increasing the resolution into the range $N_{v\|} \approx 16$ or beyond failed to recover the grid result. More narrowly, a scan over the trial-and-error parameter $\tau_\alpha$ in TCV-X21 would falsify the robustness claim if small variations of $\tau_\alpha$ around the chosen values (114 eV, 102.5 eV) changed the computed OMP density or temperature profiles by more than the quoted agreement with the grid data.

Watch

Extended reading notes

Core claim

The paper claims that a velocity-space spectral formulation based on scaled Hermite-Laguerre polynomials can replace the grid discretization in the full-f gyrokinetic code GENE-X without loss of fidelity for edge and SOL turbulence, at a fraction of the computational cost. The authors state that the spectral approach reproduces the OMP profiles of density, temperature, and radial electric field of the TCV-X21 case—turbulence dominated by trapped electron modes—with excellent agreement and significantly lower velocity-space resolution, and that a speed-up of approximately 50 is achieved for that case. They further argue that the spectral discretization is particularly advantageous at high collisionality, because the LBD collision operator has a sparse spectral representation and its CFL constraint scales linearly rather than quadratically with resolution.

Load-bearing premise

The spectral convergence and the claimed low-resolution fidelity rely on the existence of a single constant scaled temperature per species, chosen by trial and error, such that the Hermite-Laguerre basis is wide enough to represent the whole temperature range across the edge and scrape-off layer; if no such constant exists for a different scenario, the spectral resolution and speed-up will not transfer.

Editorial extensions

If this is right

  • A spectral resolution of about $(N_{v\|}, N_\mu) \approx (6,4)$ suffices to reproduce the grid simulation results for TCV-X21, compared to (80, 24) for the optimized grid case, reducing the velocity-space degrees of freedom by roughly two orders of magnitude.
  • TCV-X21-class L-mode edge/SOL gyrokinetic simulations can be completed within a few days (about 0.05–0.1 MCPUh) on current CPU-based supercomputers, versus several million CPU-hours and several weeks for the grid approach.
  • The spectral approach captures the key kinetic mechanism of the TCV-X21 validation—the collisional cooling of trapped electrons—and reproduces the electron temperature OMP profile without resolving the fine structure of the trapped-passing boundary.
  • The approach is claimed to extend to larger devices: the authors hypothesize that a similar spectral resolution is adequate for L-mode scenarios in ASDEX Upgrade-class machines, though they note the required resolution likely depends on temperature gradients and instabilities.
  • The spectral formulation exactly conserves the total energy of the GK Vlasov-Maxwell system, providing a global energy-consistency guarantee that is demonstrated analytically in the paper.

Reading between the lines

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

  • The approximately 50x speed-up and the low-resolution fidelity are demonstrated for a single scenario (TCV-X21) with a particular temperature profile; the extrapolation to other devices relies on the assumption that a single constant $\tau_\alpha$ can be found that resolves the temperature range, which may be more difficult for H-mode pedestals or detached/divertor scenarios with wider temperature
  • The spectral formulation's advantage is expected to be larger at high collisionality, because collisions damp the higher-order spectral coefficients and accelerate convergence; the collisionless limit requires higher spectral resolution, as the paper itself shows for TCV-X21.
  • The success of the approach for TEM-dominated turbulence suggests that the Hermite-Laguerre basis (which is well-suited to drift-kinetic and collisional dynamics) may also perform well for other drift-wave instabilities, but its performance for strongly electromagnetic or kinetic-ballooning regimes, where fine velocity-space structure matters, remains an open question the paper does not address.
  • The Dirichlet boundary conditions that set flux variables to zero and pin the distribution at the boundary to a local Maxwellian are acknowledged to be non-physical; improving them (e.g., with sheath boundary conditions) could change the far-SOL predictions, where the paper attributes the agreement to the boundary condition rather than to the physics.
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

5 major / 6 minor

Summary. The manuscript reports the first implementation of a velocity-space spectral (Hermite-Laguerre) discretization in the full-f gyrokinetic edge/SOL code GENE-X, replacing the grid discretization in v_parallel and mu. The authors derive the spectral form of the Vlasov-Maxwell system and the LBD collision operator, verify the implementation with the method of manufactured solutions in three geometries, and then apply the method to the TCV-X21 reference case. They compare outboard-midplane profiles of density, electron/ion temperature, and radial electric field against existing GENE-X grid results, and they report a speed-up of approximately 50 for a spectral resolution of (Nvpar,Nmu) = (6,4) relative to an optimized grid simulation with (80,24). The paper concludes that spectral acceleration enables high-fidelity edge/SOL gyrokinetic simulations on CPU-based machines within a few days and argues that this moves the method toward reactor-relevant devices.

Significance. If the central claim is correct, the paper is a substantial methodological contribution: it demonstrates that a full-f gyrokinetic edge/SOL code with X-point geometry can use a spectral velocity-space basis in place of a high-resolution grid, while preserving the physics of TEM-dominated turbulence and the kinetic collisional cooling of trapped electrons. The MMS verification in slab, circular, and toroidal geometries is a concrete strength, as is the use of openly archived grid data for comparison. The claimed order-of-magnitude reduction in velocity-space resolution and the factor-of-50 speed-up are significant for the feasibility of edge/SOL gyrokinetics on current CPU clusters. The paper is also honest about several limitations: the constant scaled temperature tau_alpha is chosen by trial and error, the far-SOL agreement is partly a boundary-condition artifact, and the required spectral resolution will vary with temperature gradients and instability type. Those caveats, however, mean that the headline transferability to ITER-like conditions is not established by the presented evidence.

major comments (5)
  1. [Sec. 6.3, Figs. 9-13] The central claim of 'excellent agreement' with the grid simulations rests exclusively on visual profile comparisons. No quantitative error metric is reported for the OMP density, Te, Ti, or Er profiles, nor are the spectral and grid profiles compared against their respective standard deviations in a quantitative way. Please compute relative L2 or pointwise errors between the spectral and grid profiles over well-defined radial ranges (e.g., edge region, near-separatrix region, and far SOL), and report how these errors vary with spectral resolution. Without such a metric, the central claim of the abstract cannot be assessed independently of visual inspection.
  2. [Sec. 7.3, Table 1] The speed-up factor of approximately 50 is based on a spectral simulation with (Nvpar,Nmu) = (6,4), but the physics validation in Section 6 uses (4,2), (6,2), (8,4), and (16,8); no OMP profile or other physical observable is shown for (6,4). Since the speed-up is the paper's headline quantitative result, the (6,4) resolution must be demonstrated to reproduce the grid profiles with the same fidelity as the resolutions studied in Section 6. Please add the (6,4) case to the profile comparison at the stated Delta_RZ = 3.7 rho_ref, or otherwise justify that interpolating between (6,2) and (8,4) is sufficient.
  3. [Sec. 4.3, Eq. (23)] The artificial Landau-damping-like sink, with coefficient K_parallel = 0.1, is stated to 'not affect the saturated turbulent state', but no sensitivity scan over K_parallel is presented. The damping term is applied precisely to the flux variables used at the low resolutions for which the speed-up is claimed, so a dependence of the saturated profiles or fluctuation level on K_parallel would directly affect the validity of the low-resolution results. Please include a K_parallel scan (at least at the resolutions used in the speed-up comparison) showing that time-averaged OMP profiles and ideally fluctuation amplitudes are unchanged.
  4. [Appendix C, Sec. 6.1, Sec. 8] The transferability of the method to H-mode or ITER-like conditions is not supported by the present evidence. Equation (C.9) gives only a lower bound on tau_alpha, and the convergence of the expansion for a Maxwellian at temperature T_alpha depends on |1 - T_alpha/tau_alpha|. Figure 8 shows that the spectral coefficients decay slowly where T_alpha deviates from tau_alpha, and Section 8 concedes that the required resolution depends on temperature gradients. For a case with T_edge/T_SOL > 10, no single constant tau_alpha will simultaneously provide fast convergence at the hot edge and at the cold SOL with only 6-8 parallel and 4 perpendicular modes. The authors should either demonstrate the approach on a synthetic or real case with a larger temperature ratio, or explicitly restrict the claim to L-mode-like cases with moderate temperature variation and remove or soften the ITER statement from the abstract.
  5. [Sec. 4.5, Sec. 6.3] The configuration-space boundary conditions used for the spectral coefficients are described as 'not based on physical principles' (projected Maxwellian with initial profiles), and the text states that the far-SOL agreement in Ti is due to the Dirichlet boundary condition rather than physical dynamics. This means that part of the claimed agreement in the SOL is a boundary artifact, not a test of the spectral method. The paper should state explicitly over which radial range the spectral/grid comparison is considered physical, and should flag the same caveat for the density and Te profiles, not only for Ti.
minor comments (6)
  1. [Table 1] The CPU/node counts appear inconsistent for the (6,2) Marconi row: 32 nodes with 48 cores per node (A3 partition) gives 1536 CPUs, not 3072. Please verify all node/core counts and the resulting CPU-hour estimates.
  2. [Sec. 6.3] The time-averaging procedure is described as a 'toroidal and time average' over 0.1 ms, but the details are incomplete: specify how the toroidal average is performed over the Nphi = 32 poloidal planes, whether the interpolation to the OMP line precedes or follows the average, and how many instantaneous samples enter the time average.
  3. [Appendix C] The sentence 'values of tau_alpha close to (but larger than) tau_alpha,c ensure a faster convergence in the SOL region but a slower convergence in the edge region' appears difficult to reconcile with Eq. (C.7): when T_alpha << tau_alpha, the factor (1 - T_alpha/tau_alpha) is near unity, which gives slow decay, not fast convergence. Please check the wording and clarify the qualitative behavior of the convergence as a function of T_alpha/tau_alpha.
  4. [Sec. 4.4, Eq. (24)] The CFL estimate in Eq. (24) uses min(R) without defining the normalization or whether R is the major radius in the normalized units of Appendix A; please define all symbols in the equation.
  5. [Sec. 2] The collision operator is called the 'Lernard-Bernstein/Dougherty' operator; the correct spelling is Lenard-Bernstein. Please correct this typo.
  6. [Sec. 4.2, Fig. 1] The description of the spectral stencil, including two-sided ghosts in p, one-sided ghosts in j, and edge ghosts for the magnetic-pumping terms, is terse. A few sentences clarifying how the truncation is applied and how ghost coefficients are set (zero? extrapolated?) would help the reader reproduce the implementation.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: the spectral formulation is derived from the GK model, verified by MMS, and benchmarked against the same group's grid GENE-X only as a consistency check, not as a fitted-input prediction.

full rationale

The central claim—that a Hermite-Laguerre spectral velocity-space discretization reproduces grid GENE-X OMP profiles at much lower resolution—does not reduce to its inputs by construction. The spectral coefficients are defined by projection (Eq. 10), the spectral Vlasov equation is derived from the conservative GK Vlasov equation (Eq. 12), and the spectral Maxwell equations are exact spectral moments (Eq. 16). The implementation is verified against analytic manufactured solutions in slab, circular, and toroidal geometry (Sec. 5, Fig. 2), providing an external functional check that does not depend on the later physics benchmarking. The scaled temperatures tau_e=114 eV and tau_i=102.5 eV are selected by trial and error to satisfy the convergence criterion Eq. (C.9) and to keep the initial simulations stable; they are not regressed against the grid OMP profiles, so the reported agreement is not a statistically forced fit. Similarly, the damping coefficient K_parallel is a numerical stabilization parameter, acknowledged to be negligible at high Nvpar (Sec. 4.3). The benchmark is the authors' own grid GENE-X simulation [25], making the comparison a same-code consistency check rather than an independent-code validation; however, the grid simulation was previously validated against TCV-X21 experimental data and its data are publicly available (Refs. 46, 47), and the spectral results are also compared with experimental TS/FHRP data. The paper explicitly flags that the far-SOL agreement is due to the Dirichlet boundary condition and should not be considered physical (Sec. 6.3), so no over-claim is hidden. The limited transferability of the ~50x speed-up to H-mode or ITER-like scenarios is a scope limitation, not circularity. No derivation step in the paper is equivalent by definition to its own inputs.

Assumptions & free parameters 3 free parameters · 7 assumptions · 0 invented entities

The central claim rests on the gyrokinetic model used by GENE-X, on the LBD collision operator, and on a set of numerical choices specific to the spectral formulation: the constant scaled temperatures tau_alpha (trial and error, tau_e=114 eV, tau_i=102.5 eV), the artificial Landau-damping coefficient K_parallel=0.1, the approximation 1/B*_parallel approximately 1/B in Eq. (11), and the unphysical Dirichlet boundary conditions. No new physical entities are introduced. The tau and K_parallel parameters are the main free parameters; neither is fitted to the grid profiles, but their adequacy is not demonstrated by sensitivity scans.

free parameters (3)
  • tau_e (electron scaled temperature) = 114 eV
    Chosen by trial and error to satisfy Eq. C.9 and stabilize the simulation; sets the width of the Hermite-Laguerre basis. Convergence degrades where local Te deviates from this value (Fig. 8).
  • tau_i (ion scaled temperature) = 102.5 eV
    Same as tau_e for ions; chosen to satisfy Eq. C.9 based on the initial Ti profile; affects spectral convergence and stability.
  • K_parallel (parallel Landau damping coefficient) = 0.1
    Artificial damping term Eq. (23) applied to flux variables with p+2j >= Nv_parallel-1. The paper states the value does not affect the saturated turbulent state, but no sensitivity scan is shown.
assumptions (7)
  • domain assumption The full-f electromagnetic and collisional gyrokinetic Vlasov-Maxwell system (Eqs. 1-5) is the correct model for edge and SOL turbulence.
    Basis of GENE-X; inherited from Ref. [18].
  • domain assumption The Lenard-Bernstein/Dougherty collision operator sufficiently models collisions for the TCV-X21 case.
    Used in both grid and spectral simulations; Ref. [31].
  • ad hoc to paper 1/B*_parallel is approximated by 1/B, neglecting the guiding-center correction in the denominator (Eq. 11).
    Needed to avoid velocity-space divisions that couple all spectral coefficients. The paper argues energy conservation is preserved, but this changes the model relative to the grid solver.
  • domain assumption The spectral expansion is truncated and boundary coefficients set to zero (Section 4.5), equivalent to truncating Eq. (8).
    Standard numerical truncation; convergence shown with increasing resolution.
  • ad hoc to paper Configuration-space Dirichlet boundary conditions use a projected Maxwellian with initial profiles (Eq. C.7); the paper states these are not based on physical principles.
    Used to ensure numerical robustness; the paper notes far-SOL agreement is due to these boundary conditions and is not physical.
  • ad hoc to paper The artificial energy sink Eq. (23) mimics Landau damping and does not affect the saturated state.
    Introduced to suppress recurrence and energy accumulation; no sensitivity study is provided.
  • ad hoc to paper The scaled temperature tau_alpha is constant in space and time.
    Chosen to avoid time-derivative and gradient terms of the basis in Eq. (12); differs from local-temperature gyro-moment models.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spectrally Accelerated Edge and Scrape-Off Layer Gyrokinetic Turbulence Simulations." pith.science (2026). https://pith.science/paper/P2ZQGZBH

@misc{pith2026241109232,
  author       = {Pith},
  title        = {Pith review of: Spectrally Accelerated Edge and Scrape-Off Layer Gyrokinetic Turbulence Simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/P2ZQGZBH}},
  note         = {Machine review of arXiv:2411.09232}
}
read the original abstract

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case [Oliveira D. S. et al., Nucl. Fusion 62, 096001 (2022)] and compare these results with previously validated grid simulations [Ulbl P. et al., Phys. Plasmas 30, 107986 (2023)]. The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. Thus, the spectral approach reduces the computational cost by at least an order of magnitude, achieving a speed-up of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

Figures

Figures reproduced from arXiv: 2411.09232 by the authors.

Figure 1
Figure 1. Velocity-space stencils associated with the grid (a) and spectral (b) implemen [PITH_FULL_IMAGE:figures/full_fig_p018_1.png] view at source ↗
Figure 2
Figure 2. MMS Lp-errors of the ion density ˆni = Nˆ 00 i for p = 2 (star markers) and p = ∞ (circular markers) as a function of the number of poloidal planes Nφ (as a proxy for the 4D increase of resolution) in the slab (blue), circular (orange), and toroidal (green) geometry. The second-order of reference is shown by the gray lines. The results of the MMS verification are shown in [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗
Figure 3
Figure 3. Normalized ion density ˆni (top) and MMS absolute error ˆni − nˆ (M) i (bottom) plotted in the φ = π poloidal plane and at tˆ = 0.625 in the toroidal geometry with Nφ = 192. the MMS errors decrease with second-order accuracy in all geometries at sufficiently high resolution (for Nφ ≳ 48), which is consistent with the dom￾inant error arising from the second-order accuracy of the Arakawa scheme used in GENE-X [18]. Th… view at source ↗
Figures from the paper (13 more)
Figure 4
Figure 4. Figure 4: Initial density n (red) and temperature Te/i (blue) profiles for the electrons (dashed) and ions (solid) as a function of the poloidal flux surface label ρpol. The density profiles are the same for both species. The horizontal blue lines represent the values of the sca…
Figure 5
Figure 5. Figure 5: Time traces of the toroidally averaged density close to the separatrix at the [PITH_FULL_IMAGE:figures/full_fig_p028_5.png]
Figure 6
Figure 6. Figure 6: Snapshots of the ion density ni in the poloidal plane φ = π at t = 0.5 ms, obtained with increasing spectral resolutions (from left to right). The green and red solid lines indicate the position of the separatrix and the OMP measurement line, respectively [PITH_FULL_I…
Figure 7
Figure 7. Figure 7: Poloidal variation of the electrostatic potential, [PITH_FULL_IMAGE:figures/full_fig_p031_7.png]
Figure 8
Figure 8. Figure 8: Normalized amplitude of the electron spectral coefficients without (left column) [PITH_FULL_IMAGE:figures/full_fig_p032_8.png]
Figure 9
Figure 9. Figure 9: The time average is calculated over a 0.1 ms period in all simulations at quasi-steady state. It is worth noticing that the OMP profiles from the grid simulations are obtained similarly. In [PITH_FULL_IMAGE:figures/full_fig_p032_9.png]
Figure 9
Figure 9. Figure 9: OMP ion density ni profiles obtained using the spectral simulations (colored lines) with different spectral resolution (Nv∥ , Nµ). The OMP profiles are computed by performing a toroidal and time average on the data interpolated on the OMP line of measurements (see red …
Figure 10
Figure 10. Figure 10: Same as Fig. 9 for the OMP electron temperature [PITH_FULL_IMAGE:figures/full_fig_p034_10.png]
Figure 11
Figure 11. Figure 11: Same as Fig. 10, but in the collisionless limit. In this case, only the spectral [PITH_FULL_IMAGE:figures/full_fig_p035_11.png]
Figure 12
Figure 12. Figure 12: Same as Fig. 10 for the OMP ion temperature [PITH_FULL_IMAGE:figures/full_fig_p036_12.png]
Figure 13
Figure 13. Figure 13: Same as Fig. 9 for the OMP radial electric field, [PITH_FULL_IMAGE:figures/full_fig_p037_13.png]
Figure 14
Figure 14. Figure 14: Roofline analysis, showing the performance (GFlops/s) as a function of the [PITH_FULL_IMAGE:figures/full_fig_p039_14.png]
Figure 15
Figure 15. Figure 15: Comparison of computational costs (MCPUh) of the spectral simulations per [PITH_FULL_IMAGE:figures/full_fig_p042_15.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Turbulence and Transport in Spectrally Accelerated full-f Gyrokinetic Simulations

    physics.plasm-ph 2025-05 conditional novelty 5.0 of 10

    Spectrally accelerated gyrokinetic simulations reproduce TEM-driven edge turbulence, particle and heat fluxes, and separatrix power in the TCV-X21 scenario, matching grid-based GENE-X and the experimental power within...

Reference graph

Works this paper leans on

52 extracted references · 51 canonical work pages · cited by 1 Pith paper

  1. [1]

    Ikeda, Progress in the iter physics basis, Nuclear Fusion 47 (6) (2007) E01

    K. Ikeda, Progress in the iter physics basis, Nuclear Fusion 47 (6) (2007) E01

    work page 2007

  2. [2]

    H. Zohm, C. Angioni, E. Fable, G. Federici, G. Gantenbein, T. Hart- mann, K. Lackner, E. Poli, L. Porte, O. Sauter, et al., On the physics guidelines for a tokamak demo, Nuclear Fusion 53 (7) (2013) 073019

    work page 2013

  3. [3]

    Litaudon, F

    X. Litaudon, F. Jenko, D. Borba, D. V. Borodin, B. J. Braams, S. Brezinsek, I. Calvo, R. Coelho, A. J. H. Donn´ e, O. Embr´ eus, et al., Eurofusion-theory and advanced simulation coordination (e-tasc): pro- gramme and the role of high performance computing, Plasma Physics and Controlled Fusion 64 (3) (2022) 034005

    work page 2022

  4. [4]

    Goerler, X

    T. Goerler, X. Lapillonne, S. Brunner, T. Dannert, F. Jenko, F. Merz, D. Told, The global version of the gyrokinetic turbulence code GENE, Journal of Computational Physics 230 (18) (2011) 7053

    work page 2011

  5. [5]

    Grandgirard, J

    V. Grandgirard, J. Abiteboul, J. Bigot, T. Cartier-Michaud, N. Crou- seilles, G. Dif-Pradalier, C. Ehrlacher, D. Esteve, X. Garbet, P. Ghen- drih, G. Latu, M. Mehrenberger, C. Norscini, C. Passeron, F. Rozar, Y. Sarazin, E. Sonnendr¨ ucker, A. Strugarek, D. Zarzoso, A 5D gy- rokinetic full-f global semi-lagrangian code for flux-driven ion turbu- lence sim...

    work page doi:10.1016/j.cpc.2016.05.007 2016

  6. [6]

    Lanti, N

    E. Lanti, N. Ohana, N. Tronko, T. Hayward-Schneider, A. Bottino, B. F. McMillan, A. Mishchenko, A. Scheinberg, A. Biancalani, P. Angelino, et al., Orb5: a global electromagnetic gyrokinetic code using the pic approach in toroidal geometry, Computer Physics Communications 251 (2020) 107072

    work page 2020

  7. [7]

    Jenko, W

    F. Jenko, W. Dorland, M. Kotschenreuther, B. N. Rogers, Electron tem- perature gradient driven turbulence, Physics of plasmas 7 (5) (2000) 1904

    work page 2000

  8. [8]

    Candy, E

    J. Candy, E. A. Belli, R. V. Bravenec, A high-accuracy eulerian gyroki- netic solver for collisional plasmas, Journal of Computational Physics 324 (2016) 73

    work page 2016

Show all 52 references
  1. [9]

    A. G. Peeters, Y. Camenen, F. J. Casson, W. A. Hornsby, A. P. Snodin, D. Strintzi, G. Szepesi, The nonlinear gyro-kinetic flux tube code gkw, Computer physics communications 180 (12) (2009) 2650

    2009

  2. [10]

    Giacomin, P

    M. Giacomin, P. Ricci, A. Coroado, G. Fourestey, D. Galassi, E. Lanti, D. Mancini, N. Richart, L. Stenger, N. Varini, The GBS code for the self- consistent simulation of plasma turbulence and kinetic neutral dynamics in the tokamak boundary, J. Comput. Phys. 463 (2021) 111294

    2021

  3. [11]

    Tamain, H

    P. Tamain, H. Bufferand, G. Ciraolo, C. Colin, D. Galassi, P. Ghendrih, F. Schwander, E. Serre, The tokam3x code for edge turbulence fluid sim- ulations of tokamak plasmas in versatile magnetic geometries, Journal of Computational Physics 321 (2016) 606

    2016

  4. [12]

    Stegmeir, A

    A. Stegmeir, A. Ross, T. Body, M. Francisquez, W. Zholobenko, D. Coster, O. Maj, P. Manz, F. Jenko, B. N. Rogers, et al., Global tur- bulence simulations of the tokamak edge region with GRILLIX, Physics of Plasmas 26 (5) (2019) 052517

    2019

  5. [13]

    Hariri, M

    F. Hariri, M. Ottaviani, A flux-coordinate independent field-aligned ap- proach to plasma turbulence simulations, Computer Physics Communi- cations 184 (11) (2013) 2419

    2013

  6. [14]

    Stegmeir, D

    A. Stegmeir, D. Coster, O. Maj, K. Hallatschek, K. Lackner, The field line map approach for simulations of magnetically confined plasmas, Computer Physics Communications 198 (2016) 139. 49

    2016

  7. [15]

    Zholobenko, T

    W. Zholobenko, T. Body, P. Manz, A. Stegmeir, B. Zhu, M. Griener, G. D. Conway, D. Coster, F. Jenko, A. U. Team, et al., Electric field and turbulence in global braginskii simulations across the asdex upgrade edge and scrape-off layer, Plasma Physics and Controlled Fusion 63 (...

    2021

  8. [16]

    Stegmeir, T

    A. Stegmeir, T. Body, W. Zholobenko, Analysis of locally-aligned and non-aligned discretisation schemes for reactor-scale tokamak edge tur- bulence simulations, Computer Physics Communications 290 (2023) 108801

    2023

  9. [17]

    T. Body, A. Stegmeir, W. Zholobenko, D. Coster, F. Jenko, Treatment of advanced divertor configurations in the flux-coordinate independent tur- bulence code GRILLIX, Contributions to Plasma Physics 60 (5) (2020) e201900139

    2020

  10. [18]

    Michels, A

    D. Michels, A. Stegmeir, P. Ulbl, D. Jarema, F. Jenko, GENE-X: A full- f gyrokinetic turbulence code based on the flux-coordinate independent approach, Comput. Phys. Commun. 264 (2021) 107986

    2021

  11. [19]

    E. L. Shi, G. W. Hammett, T. Stoltzfus-Dueck, A. Hakim, Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma, Journal of Plasma Physics 83 (3) (2017) 905830304

    2017

  12. [20]

    Boesl, A

    M. Boesl, A. Bergmann, A. Bottino, D. Coster, E. Lanti, N. Ohana, F. Jenko, Gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS, Physics of Plasmas 26 (12) (2019) 122302

    2019

  13. [21]

    Ku, C.-S

    S. Ku, C.-S. Chang, P. H. Diamond, Full-f gyrokinetic particle simula- tion of centrally heated global itg turbulence from magnetic axis to edge pedestal top in a realistic tokamak geometry, Nuclear Fusion 49 (11) (2009) 115021

    2009

  14. [22]

    Hager, S

    R. Hager, S. Ku, A. Y. Sharma, C. S. Chang, R. M. Churchill, A. Schein- berg, Electromagnetic total-f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC, Physics of Plasmas 29 (11) (2022) 112308

    2022

  15. [23]

    M. Dorf, M. Dorr, Continuum gyrokinetic simulations of edge plasmas in single-null geometries, Physics of Plasmas 28 (3) (2021) 032508. 50

    2021

  16. [24]

    Michels, P

    D. Michels, P. Ulbl, W. Zholobenko, T. Body, A. Stegmeir, T. Eich, M. Griener, G. Conway, F. Jenko, Full-f electromagnetic gyrokinetic turbulence simulations of the edge and scrape-off layer of ASDEX up- grade with GENE-X, Physics of Plasmas 29 (3) (2022) 032307

    2022

  17. [25]

    P. Ulbl, T. Body, W. Zholobenko, A. Stegmeir, J. Pfennig, F. Jenko, Influence of collisions on the validation of global gyrokinetic simulations in the edge and scrape-off layer of TCV, Physics of Plasmas 30 (5) (2023) 107986

    2023

  18. [26]

    D. S. Oliveira, T. Body, D. Galassi, C. Theiler, E. Laribi, P. Tamain, A. Stegmeir, M. Giacomin, W. Zholobenko, P. Ricci, et al., Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case, Nuclear Fusion 62 (9) (2022) 096001

    2022

  19. [27]

    P. J. Roache, Code verification by the method of manufactured solutions, J. Fluids Eng. 124 (1) (2002) 4

    2002

  20. [28]

    B. J. Frei, A. Hoffmann, P. Ricci, S. Brunner, Z. Tecchioll, Moment- based approach to the flux-tube linear gyrokinetic model, Journal of Plasma Physics 89 (4) (2023) 905890414

    2023

  21. [29]

    Lenard, I

    A. Lenard, I. B. Bernstein, Plasma oscillations with diffusion in velocity space, Physical Review 112 (5) (1958) 1456

    1958

  22. [30]

    J. P. Dougherty, Model Fokker-Planck equation for a plasma and its solution, The Physics of Fluids 7 (11) (1964) 1788

    1964

  23. [31]

    P. Ulbl, D. Michels, F. Jenko, Implementation and verifica- tion of a conservative, multi-species, gyro-averaged, full-f, Lenard- Bernstein/Dougherty collision operator in the gyrokinetic code GENE- X, Contributions to Plasma Physics 62 (5) (2022) e202100180

    2022

  24. [32]

    N. R. Mandell, A. Hakim, G. W. Hammett, M. Francisquez, Electromag- netic full-gyrokinetics in the tokamak edge with discontinuous galerkin methods, Journal of Plasma Physics 86 (1) (2020) 905860109

    2020

  25. [33]

    Scott, Turbulence and Instabilities in Magnetised Plasmas, Volume 2: Gyrokinetic theory and gyrofluid turbulence, IOP Publishing, 2021

    B. Scott, Turbulence and Instabilities in Magnetised Plasmas, Volume 2: Gyrokinetic theory and gyrofluid turbulence, IOP Publishing, 2021. 51

    2021

  26. [34]

    B. J. Frei, J. Mencke, P. Ricci, Full-f turbulent simulation in a linear plasma device using a gyro-moment approach, Physics of Plasmas 31 (2023) 012301

    2023

  27. [35]

    B. J. Frei, R. Jorge, P. Ricci, A gyrokinetic model for the plasma pe- riphery of tokamak devices, Journal of Plasma Physics 86 (2) (2020) 905860205

    2020

  28. [36]

    I. S. Gradshteyn, I. M. Ryzhik, Table of integrals, series, and products, Academic press, 2014

    2014

  29. [37]

    L. D. Landau, Kinetic equation for the coulomb effect, Phys. Z. Sowje- tunion 10 (1936) 154

    1936

  30. [38]

    Jorge, B

    R. Jorge, B. J. Frei, P. Ricci, Nonlinear gyrokinetic coulomb collision operator, Journal of Plasma Physics 85 (6) (2019) 905850604

    2019

  31. [39]

    B. J. Frei, J. Ball, A. C. D. Hoffmann, R. Jorge, P. Ricci, L. Stenger, Development of advanced linearized gyrokinetic collision operators using a moment approach, Journal of Plasma Physics 87 (5) (2021) 905870501

    2021

  32. [40]

    Arakawa, Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow

    A. Arakawa, Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow. Part I, Journal of computational physics 135 (2) (1997) 103

    1997

  33. [41]

    Michels, Development of a high-performance gyrokinetic turbulence code for the edge and scrape-off layer of magnetic confinement fusion devices, Ph.D

    D. Michels, Development of a high-performance gyrokinetic turbulence code for the edge and scrape-off layer of magnetic confinement fusion devices, Ph.D. thesis, Technische Universit¨ at M¨ unchen (2021)

    2021

  34. [42]

    Courant, K

    R. Courant, K. Friedrichs, H. Lewy, ¨Uber die partiellen differenzengle- ichungen der mathematischen physik, Mathematische annalen 100 (1) (1928) 32

    1928

  35. [43]

    Loizu, P

    J. Loizu, P. Ricci, F. D. Halpern, S. Jolliet, Boundary conditions for plasma fluid models at the magnetic presheath entrance, Physics of plas- mas 19 (12) (2012)

    2012

  36. [44]

    F. Riva, P. Ricci, F. D. Halpern, S. Jolliet, J. Loizu, A. Mosetto, Verifica- tion methodology for plasma simulations and application to a scrape-off layer turbulence code, Physics of Plasmas 21 (6) (2014). 52

    2014

  37. [45]

    S. Coda, M. Agostini, R. Albanese, S. Alberti, E. Alessi, S. Allan, J. All- cock, R. Ambrosino, H. Anand, Y. Andrebe, et al., Physics research on the TCV tokamak facility: from conventional to alternative scenarios and beyond, Nuclear Fusion 59 (11) (2019) 112023

    2019

  38. [46]

    P. Ulbl, T. Body, W. Zholobenko, A. Stegmeir, J. Pfennig, F. Jenko, TCV-X21-GENEX: influence of collisions on the validation of global gyrokinetic simulations (2023). doi:10.5281/ZENODO.7894731. URL https://zenodo.org/record/7894731

    2023

  39. [47]

    Ulbl, Collision models for gyrokinetic simulations of edge turbulence in fusion plasmas, Ph.D

    P. Ulbl, Collision models for gyrokinetic simulations of edge turbulence in fusion plasmas, Ph.D. thesis, Technische Universit¨ at M¨ unchen (2023)

    2023

  40. [48]

    A. Ross, A. Stegmeir, D. Coster, Effect of the boussinesq approximation: Turbulence studies with GRILLIX in slab geometry, Contributions to Plasma Physics 58 (6) (2018) 478

    2018

  41. [49]

    D. S. D. Oliveira, T. Body, D. Galassi, C. Theiler, TCV-X21: an open dataset for the validation of edge turbulence models (Mar. 2024). doi: 10.5281/zenodo.10841179. URL https://doi.org/10.5281/zenodo.10841179

    2024 doi

  42. [50]

    Zholobenko, A

    W. Zholobenko, A. Stegmeir, M. Griener, G. D. Conway, T. Body, D. Coster, F. Jenko, A. U. Team, et al., The role of neutral gas in validated global edge turbulence simulations, Nuclear Fusion 61 (11) (2021) 116015

    2021

  43. [51]

    Williams, A

    S. Williams, A. Waterman, D. Patterson, Roofline: an insightful visual performance model for multicore architectures, Communications of the ACM 52 (4) (2009) 65

    2009

  44. [52]

    Madsen, Full-F gyrofluid model, Physics of Plasmas 20 (7) (2013) 072301

    J. Madsen, Full-F gyrofluid model, Physics of Plasmas 20 (7) (2013) 072301. 53

    2013

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.