Gyrokinetic simulations on zonal flow-turbulence spreading coupling
Pith reviewed 2026-05-09 23:38 UTC · model grok-4.3
The pith
Global nonlinear gyrokinetic simulations demonstrate that turbulence spreading transports zonal flows radially into linearly stable regions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using global nonlinear gyrokinetic simulations, turbulence spreading is shown to transport zonal flow radially, extending into the linearly stable regions after local nonlinear saturation of turbulence. This is understood through analysis with a momentum theorem in toroidal plasmas, indicating a direct relation between turbulence-driven enstrophy transport and perpendicular momentum generation.
What carries the argument
The momentum theorem that relates turbulence-driven enstrophy transport to perpendicular momentum generation, applied to interpret the radial transport by turbulence spreading.
If this is right
- Zonal flows can extend their influence into linearly stable plasma regions via turbulence spreading.
- Radial transport of perpendicular momentum is tied to enstrophy carried by spreading turbulence.
- Local nonlinear saturation of turbulence does not confine its effects to unstable regions.
- The coupling between zonal flows and turbulence spreading affects how momentum is generated and transported in toroidal plasmas.
Where Pith is reading between the lines
- This radial transport could lead to broader regulation of turbulence across the plasma profile than local analysis predicts.
- Experimental measurements of fluctuations in stable regions might detect signatures of this spreading-induced zonal flow presence.
- Similar mechanisms may apply to other types of plasma instabilities or different magnetic configurations.
Load-bearing premise
The momentum theorem applies directly and without extra assumptions to the global simulation results to establish the radial transport mechanism.
What would settle it
Observation of turbulence spreading into stable regions without corresponding radial transport of zonal flow or without the expected correlation between enstrophy transport and momentum generation would contradict the claim.
Figures
read the original abstract
Zonal flows and turbulence spreading play important roles in magnetic fusion plasma confinement, yet their coupling mechanisms remain elusive. Using global nonlinear gyrokinetic simulations, we show that turbulence spreading transports zonal flow radially, extending into the linearly stable regions after local nonlinear saturation of turbulence. Theoretical understanding has been gained by analyzing the simulation results in the context of a momentum theorem in toroidal plasmas [T.S. Hahm \textit{et al.}, Phys. Plasmas \textbf{31}, 032310 (2024)] which is an extension of the Charney-Drazin non-acceleration theorem [J.G. Charney and P.G. Drazin, J. Geophys. Res. \textbf{66}, 83 (1961)]. It indicates a direct relation between turbulence-driven enstrophy transport and perpendicular momentum generation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses global nonlinear gyrokinetic simulations to demonstrate that turbulence spreading transports zonal flows radially into linearly stable regions after local nonlinear saturation. Theoretical interpretation is provided by analyzing the results in the context of the momentum theorem from Hahm et al. (Phys. Plasmas 31, 032310, 2024), an extension of the Charney-Drazin non-acceleration theorem, which is claimed to establish a direct relation between turbulence-driven enstrophy transport and perpendicular momentum generation.
Significance. If the central claim holds after validation, the work would provide a concrete mechanistic link between zonal-flow spreading and enstrophy transport in toroidal plasmas, extending non-acceleration theorems to the nonlinear gyrokinetic regime with potential implications for understanding confinement and transport barriers in fusion devices. The global simulation approach to probe spreading beyond linearly unstable regions is a strength, though the result's robustness depends on confirming the theorem's applicability.
major comments (2)
- [Interpretation / Discussion] The interpretation section (likely §4 or equivalent): the claimed direct relation between observed zonal-flow spreading and the momentum theorem rests on the unverified assumption that the Hahm et al. (2024) theorem applies without modification to the global nonlinear simulation regime. Global simulations include radial profile variations, finite domain boundaries, and possible residual sources that can introduce additional terms or violate the theorem's underlying assumptions (e.g., specific averaging, absence of external torques, or quasilinear ordering); no dedicated check confirming these conditions hold unmodified is described.
- [Results] Results section (likely §3): the central finding is stated without accompanying quantitative profiles of zonal-flow transport, error estimates, resolution convergence checks, or direct numerical comparison to the theorem equations (e.g., enstrophy flux terms). This weakens the data-to-claim linkage, as the abstract and summary provide no simulation parameters (e.g., grid resolution, time steps, or plasma parameters) to allow independent assessment of the spreading mechanism.
minor comments (2)
- [Abstract] The abstract could be expanded to include at least one key quantitative indicator (e.g., radial extent of spreading or a normalized transport coefficient) to better convey the result's magnitude.
- [Theoretical background] Notation for the momentum theorem terms should be explicitly cross-referenced to the equations in Hahm et al. (2024) for clarity when applying them to the simulation diagnostics.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments help clarify how to strengthen the linkage between the simulation results and the momentum theorem. We address each major comment below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Interpretation / Discussion] The interpretation section (likely §4 or equivalent): the claimed direct relation between observed zonal-flow spreading and the momentum theorem rests on the unverified assumption that the Hahm et al. (2024) theorem applies without modification to the global nonlinear simulation regime. Global simulations include radial profile variations, finite domain boundaries, and possible residual sources that can introduce additional terms or violate the theorem's underlying assumptions (e.g., specific averaging, absence of external torques, or quasilinear ordering); no dedicated check confirming these conditions hold unmodified is described.
Authors: We agree that explicit verification of the theorem's assumptions is necessary for a robust claim. Our simulations are performed without external torques and employ flux-surface averaging consistent with the theorem's derivation. However, the original manuscript does not contain a dedicated subsection confirming that radial profile variations and boundary effects introduce no significant extra terms. In the revision we will add a paragraph in the interpretation section that (i) restates the key assumptions of Hahm et al. (2024), (ii) shows that the residual sources remain negligible in our runs, and (iii) demonstrates that the observed enstrophy flux and perpendicular momentum generation remain directly correlated as predicted by the theorem even when modest profile variations are present. This will make the applicability explicit rather than implicit. revision: yes
-
Referee: [Results] Results section (likely §3): the central finding is stated without accompanying quantitative profiles of zonal-flow transport, error estimates, resolution convergence checks, or direct numerical comparison to the theorem equations (e.g., enstrophy flux terms). This weakens the data-to-claim linkage, as the abstract and summary provide no simulation parameters (e.g., grid resolution, time steps, or plasma parameters) to allow independent assessment of the spreading mechanism.
Authors: We accept that the results section would be strengthened by additional quantitative material. The full manuscript already reports the essential simulation parameters (grid resolution, time step, and plasma parameters) in the methods section, but they are not summarized in a single table. In the revision we will (i) add a table of key numerical parameters, (ii) include radial profiles of the zonal-flow transport and the relevant enstrophy flux terms with direct overlay of the theorem predictions, (iii) report statistical error estimates obtained from multiple independent realizations, and (iv) present resolution-convergence tests for the spreading distance and the enstrophy-momentum correlation. These additions will make the data-to-claim connection quantitative and reproducible. revision: yes
Circularity Check
Central mechanistic interpretation of simulation results as radial zonal-flow transport rests on self-cited momentum theorem without re-derivation or validation
specific steps
-
self citation load bearing
[Abstract]
"Theoretical understanding has been gained by analyzing the simulation results in the context of a momentum theorem in toroidal plasmas [T.S. Hahm et al., Phys. Plasmas 31, 032310 (2024)] which is an extension of the Charney-Drazin non-acceleration theorem [J.G. Charney and P.G. Drazin, J. Geophys. Res. 66, 83 (1961)]. It indicates a direct relation between turbulence-driven enstrophy transport and perpendicular momentum generation."
The paper claims the simulations reveal a direct relation between enstrophy transport and perpendicular momentum (i.e., radial zonal-flow transport via spreading). This relation is not extracted from the simulation data alone but is supplied by interpreting the data through the momentum theorem of Hahm et al. (2024). Because the cited paper shares an author with the present work and the current manuscript does not re-derive or validate the theorem's assumptions in the global regime, the central interpretive step is justified only by self-citation.
full rationale
The paper's simulations demonstrate turbulence spreading into linearly stable regions post-saturation, but the abstract explicitly states that theoretical understanding and the 'direct relation between turbulence-driven enstrophy transport and perpendicular momentum generation' are gained solely by analyzing results 'in the context of' the momentum theorem from Hahm et al. (2024). With T.S. Hahm as co-author on both works, this self-citation supplies the load-bearing link between observed spreading and the claimed radial transport mechanism. The provided text gives no indication that the theorem's assumptions (e.g., specific averaging, absence of external torques, or quasilinear ordering) are re-derived or explicitly checked against the global nonlinear gyrokinetic setup, so the unification claim reduces to application of prior self-work rather than independent derivation from the current simulations.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Momentum theorem in toroidal plasmas (Hahm et al. 2024) that extends the Charney-Drazin non-acceleration theorem and relates turbulence-driven enstrophy transport to perpendicular momentum generation
Reference graph
Works this paper leans on
-
[1]
M. Cao and P. H. Diamond. Physics of edge-core coupling by inward turbulence propagation. Physical Review Letters. 2025. doi:10.1103/PhysRevLett.134.235101
-
[2]
J. G. Charney and P. G. Drazin. Propagation of planetary-scale disturbances from the lower into the upper atmosphere. Journal of Geophysical Research. 1961. doi:10.1029/JZ066i001p00083
-
[3]
C.-C. Chen and P. H. Diamond. Potential vorticity mixing in a tangled magnetic field. The Astrophysical Journal. 2020. doi:10.3847/1538-4357/ab774f
-
[4]
L. Chen and Z. Lin and R. White. Excitation of zonal flow by drift waves in toroidal plasmas. Physics of Plasmas. 2000. doi:10.1063/1.874222
-
[5]
L. Chen and R. B. White and F. Zonca. Zonal-flow dynamics and size scaling of anomalous transport. Physical Review Letters. 2004. doi:10.1103/PhysRevLett.92.075004
-
[6]
M. J. Choi and L. Bard \=o czi and J.-M. Kwon and T. S. Hahm and H. K. Park and J. Kim and M. Woo and B.-H. Park and G. S. Yun and E. Yoon and G. McKee. Effects of plasma turbulence on the nonlinear evolution of magnetic island in tokamak. Nature Communications. 2021. doi:10.1038/s41467-020-20652-9
-
[7]
W. Deng and Z. Lin. Properties of microturbulence in toroidal plasmas with reversed magnetic shear. Physics of Plasmas. 2009. doi:10.1063/1.3243918
-
[8]
P. H. Diamond and M. N. Rosenbluth and F. L. Hinton and M. Malkov and J. Fleischer and A. Smolyakov. Dynamics of zonal flows and self-regulating drift-wave turbulence. 17th IAEA Fusion Energy Conference, Yokohama, Japan. 1998
work page 1998
-
[9]
P. H. Diamond and S.-I. Itoh and K. Itoh and T. S. Hahm. Zonal flows in plamsa - a review. Plasma Physics and Controlled Fusion. 2005. doi:10.1088/0741-3335/47/5/R01
-
[10]
P. H. Diamond and \"O . D. G \"u rcan and T. S. Hahm and K. Miki and Y. Kosuga and X. Garbet. Momentum theorems and the structure of atmospheric jets and zonal flows in plasmas. Plasma Physics and Controlled Fusion. 2008. doi:10.1088/0741-3335/50/12/124018
-
[11]
A. M. Dimits and G. Bateman and M. A. Beer and B. I. Cohen and W. Dorland and G. W. Hammett and C. Kim and J. E. Kinsey and M. Kotschenreuther and A. H. Kritz and L. L. Lao and J. Mandrekas and W. M. Nevins and S. E. Parker and A. J. Redd and D. E. Shumaker and R. Sydora and J. Weiland. Comparisons and physics basis of tokamak transport models and turbule...
-
[12]
D. G. Dritschel and M. E. McIntyre. Multiple jets as PV staircase: The Phillips effect and the resilience of eddy-transport barriers. Journal of the Atmospheric Sciences. 2008. doi:10.1175/2007JAS2227.1
-
[13]
T. Estrada and C. Hidalgo and T. Happel. Signatures of turbulence spreading during the H–L back-transition in TJ-II plasmas. Nuclear Fusion. 2011. doi:10.1088/0029-5515/51/3/032001
-
[14]
T. Estrada and E. Ascas \'i bar and E. Blanco and A. Cappa and C. Hidalgo and K. Ida and A. L \'o pez-Fraguas and B. Ph van Milligen. Plasma flow, turbulence and magnetic islands in TJ-II. Nuclear Fusion. 2016. doi:10.1088/0029-5515/56/2/026011
-
[15]
R. A. Fisher. The wave of advance of advantageous genes. Annals of Eugenics. 1937. doi:10.1111/j.1469-1809.1937.tb02153.x
-
[16]
X. Garbet and L. Laurent and A. Samain and J. Chinardet. Radial propagation of turbulence in tokamaks. Nuclear Fusion. 1994. doi:10.1088/0029-5515/34/7/I04
-
[17]
X. Garbet and Y. Idomura and L. Villard and T. H. Watanabe. Gyrokinetic simulations of turbulent transport. Nuclear Fusion. 2010. doi:10.1088/0029-5515/50/4/043002
-
[18]
G. Grenfell and M. Spolaore and D. Abate and L. Carraro and L. Marrelli and I. Predebon and S. Spagnolo and M. Veranda and M. Agostini and B. Ph. van Milligen. Turbulent filament properties in L and H-mode regime in the RFX-mod operating as a tokamak. Nuclear Fusion. 2020. doi:10.1088/1741-4326/abaf32
-
[19]
\"O . D. G \"u rcan and P. H. Diamond and T. S. Hahm and Z. Lin. Dynamics of turbulence spreading in magnetically confined plasmas. Physics of Plasmas. 2005. doi:10.1063/1.1853385
-
[20]
\"O . D. G \"u rcan and P. H. Diamond and T. S. Hahm. Radial transport of fluctuation energy in a two-field model of drift-wave turbulence. Physics of Plasmas. 2006. doi:10.1063/1.2180668
-
[21]
\"O . D. G \"u rcan and P. H. Diamond and T. S. Hahm. Nonlinear triad interactions and the mechanism of spreading in drift-wave turbulence. Physical Review Letters. 2006. doi:10.1103/PhysRevLett.97.024502
-
[22]
\"O . D. G \"u rcan and P. H. Diamond. Zonal flows and pattern formation. Journal of Physics A: Mathematical and Theoretical. 2015. doi:10.1088/1751-8113/48/29/293001
-
[23]
P. N. Guzdar and R. G. Kleva and L. Chen. Shear flow generation by drift waves revisited. Physics of Plasmas. 2001. doi:10.1063/1.1340618
-
[24]
T. S. Hahm. Nonlinear gyrokinetic equations for tokamak microturbulence. Physics of Fluids. 1988. doi:10.1063/1.866544
-
[25]
T. S. Hahm and K. H. Burrell. Flow shear induced fluctuation suppression in finite aspect ratio shaped tokamak plasma. Physics of Plasmas. 1995. doi:10.1063/1.871313
-
[26]
T. S. Hahm and P. H. Diamond and Z. Lin and K. Itoh and S.-I. Itoh. Turbulence spreading into the linearly stable zone and transport scaling. Plasma Physics and Controlled Fusion. 2004. doi:10.1088/0741-3335/46/5A/036
-
[27]
T. S. Hahm and P. H. Diamond and Z. Lin and G. Rewoldt and \"O . G \"u rcan and S. Ethier. On the dynamics of edge-core coupling. Physics of Plasmas. 2005. doi:10.1063/1.2034307
-
[28]
T. S. Hahm and L. Wang and J. Madsen. Fully electromagnetic nonlinear gyrokinetic equations for tokamak edge turbulence. Physics of Plasmas. 2009. doi:10.1063/1.3073671
-
[29]
T. S. Hahm and P. H. Diamond. Mesoscopic transport events and the breakdown of Fick's law for turbulent fluxes. Journal of the Korean Physical Society. 2018. doi:10.3938/jkps.73.747
-
[30]
T. S. Hahm and G. J. Choi and S. J. Park and Y.-S. Na. Fast ion effects on zonal flow generation: A simple model. Physics of Plasmas. 2023. doi:10.1063/5.0151466
-
[31]
T. S. Hahm and P. H. Diamond and S. J. Park and Y.-S. Na. Potential vorticity conservation for plasma turbulence in an inhomogeneous magnetic field: Theory and implications. Physics of Plasmas. 2024. doi:10.1063/5.0189855
-
[32]
A. Hasegawa and M. Wakatani. Plasma edge turbulence. Physical Review Letters. 1983. doi:10.1103/PhysRevLett.50.682
-
[33]
P.-C. Hsu and P. H. Diamond. On calculating the potential vorticity flux. Physics of Plasmas. 2015. doi:10.1063/1.4916401
-
[34]
K. Ida and T. Kobayashi and M. Ono and T. E. Evans and G. R. McKee and M. E. Austin. Hysteresis relation between turbulence and temperature modulation during the heat pulse propagation into a magnetic island in DIII-D. Physical Review Letters. 2018. doi:10.1103/PhysRevLett.120.245001
-
[35]
F.O. Khabanov and R. Hong and P. H. Diamond and G. R. Tynan and Z. Yan and G. R. McKee and C. Chrystal and F. Scotti and G. Yu and S. A. Zamperini and Y. Zhu. Density fluctuation statistics and turbulence spreading at the edge of L–mode plasmas. Nuclear Fusion. 2024. doi:10.1088/1741-4326/ad820d
-
[36]
M. Kobayashi and K. Tanaka and K. Ida and Y. Hayashi and Y. Takemura and T. Kinoshita. Turbulence spreading into an edge stochastic magnetic layer induced by magnetic fluctuation and its impact on divertor heat load. Physical Review Letters. 2022. doi:10.1103/PhysRevLett.128.125001
-
[37]
A. N. Kolmogoroff and I. G. Petrovsky and N. S. Piscounoff. Clinical Cancer Research. 1937
work page 1937
-
[38]
Y. Kosuga and P. H. Diamond and L. Wang and \"O . D. G \"u rcan and T. S. Hahm. Progress on theoretical issues in modelling turbulent transport. Nuclear Fusion. 2013. doi:10.1088/0029-5515/53/4/043008
-
[39]
J. M. Kwon and S. Yi and T. Rhee and P. H. Diamond and K. Miki and T. S. Hahm and J. Y. Kim and \"O . D. G \"u rcan and C. McDevitt. Analysis of symmetry breaking mechanisms and the role of turbulence self-regulation in intrinsic rotation. Nuclear Fusion. 2012. doi:10.1088/0029-5515/52/1/013004
-
[40]
J.-M. Kwon and Lei Qi and S. Yi and T. S. Hahm. ITG-TEM turbulence simulation with bounce-averaged kinetic electrons in tokamak geometry. Computer Physics Communications. 2017. doi:10.1016/j.cpc.2017.02.009
-
[41]
W. W. Lee. Gyrokinetic approach in particle simulation. Physics of Fluids. 1983. doi:10.1063/1.864140
-
[42]
N. Li and X. Q. Xu and P. H. Diamond and T. Zhang and X. Liu and Y. F. Wang and N. Yan and G. S. Xu. How fluctuation intensity flux drives SOL expansion. Nuclear Fusion. 2023. doi:10.1088/1741-4326/ad0599
-
[43]
Z. Lin and T. S. Hahm and W. W. Lee and W. M. Tang and R. B. White. Turbulent transport reduction by zonal flows: Massively parallel simulations. Science. 1998. doi:10.1126/science.281.5384.1835
-
[44]
Z. Lin and S. Ethier and T. S. Hahm and W. M. Tang. Size scaling of turbulent transport in magnetically confined plasmas. Physical Review Letters. 2002. doi:10.1103/PhysRevLett.88.195004
-
[45]
J. Madsen. Full-F gyrofluid mode. Physics of Plasmas. 2013. doi:10.1063/1.4813241
-
[46]
H. Masui and A. Ishizawa and K. Imadera and Y. Kishimoto and Y. Nakamura. Global saturation physics of ion temperature gradient turbulence in finite normalized pressure tokamaks. Nuclear Fusion. 2022. doi:10.1088/1741-4326/ac667f
-
[47]
C. J. McDevitt and P. H. Diamond and \"O . D. G \"u rcan and T. S. Hahm. Poloidal rotation and its relation to the potential vorticity flux. Physics of Plasmas. 2010. doi:10.1063/1.3490253
-
[48]
M. E. McIntyre. How Well do we Understand the Dynamics of Stratospheric Warmings?. Journal of the Meteorological Society of Japan. Ser. II. 1982. doi:10.2151/jmsj1965.60.1_37
- [49]
-
[50]
P. K. Shukla and L. Stenflo. Nonlinear interactions between drift waves and zonal flows. The European Physical Journal D - Atomic, Molecular, Optical and Plasma Physics. 2002. doi:10.1140/epjd/e2002-00119-x
-
[51]
J. H. Seo and J.-M. Kwon and S. Yi and L. Qi. A new hybrid simulation model for tokamak plasma turbulence. Computer Physics Communications. 2021. doi:10.1016/j.cpc.2020.107626
-
[52]
D. Strintzi and B. Scott. Nonlocal nonlinear electrostatic gyrofluid equations. Physics of Plasmas. 2004. doi:10.1063/1.1807850
-
[53]
G. I. Taylor. I. Eddy motion in the atmosphere. Philosophical Transactions of the Royal Society of London. Series A. 1915. doi:10.1098/rsta.1915.0001
-
[54]
G. K. Vallis. Atmospheric and Oceanic Fluid Dynamics, 2nd ed. 2017
work page 2017
-
[55]
L. Wang and T. S. Hahm. Generalized expression for polarization density. Physics of Plasmas. 2009. doi:10.1063/1.3152601
-
[56]
L. Wang and P. H. Diamond and T. S. Hahm. How does drift wave turbulence convert parallel compression into perpendicular flows?. Plasma Physics and Controlled Fusion. 2012. doi:10.1088/0741-3335/54/9/095015
-
[57]
W. X. Wang and T. S. Hahm and W. W. Lee and G. Rewoldt and J. Manickam and W. M. Tang. Nonlocal properties of gyrokinetic turbulence and the role of E B flow shear. Physics of Plasmas. 2007. doi:10.1063/1.2750647
-
[58]
M. Yagi and T. Ueda and S.-I. Itoh and M. Azumi and K. Itoh and P. H. Diamond and T. S. Hahm. Turbulence spreading in reversed shear plasmas. Plasma Physics and Controlled Fusion. 2006. doi:10.1088/0741-3335/48/5A/S42
-
[59]
S. Yi and H. Jhang and J.-M. Kwon and S. S. Kim. Role of parallel compression in potential vorticity mixing and zonal flow generation: a gyrokinetic simulation study. Nuclear Fusion. 2019. doi:10.1088/1741-4326/aaf9d3
-
[60]
S. Yi and C. Sung and E. S. Yoon and J.-M. Kwon and T. S. Hahm and D. Kim and J. Kang and J. Seo and Y. W. Cho and L. Qi. A validation study of a bounce-averaged kinetic electron model in a KSTAR L-mode plasma. Physics of Plasmas. 2024. doi:10.1063/5.0178350
-
[61]
E. S. Yoon and T. S. Hahm and G. J. Choi and Y. W. Cho and A. Ishizawa and M. J. Choi and J. M. Kwon. Turbulence spreading induced E B vortex flow generation in a magnetic island. Nuclear Fusion. 2024. doi:10.1088/1741-4326/ad82fa
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.