pith. sign in

arxiv: 2604.03089 · v1 · submitted 2026-04-03 · ⚛️ physics.plasm-ph

Finite Ion Temperature Effects on the Merging of Current-Carrying ELM Filaments in the edge region of a tokamak

Pith reviewed 2026-05-13 19:04 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords ELM filamentsfinite ion temperaturefilament mergingtokamak edge plasmafluid modelradial transportpoloidal flowsrotational motion
0
0 comments X

The pith

Increasing ion temperature delays the merging of ELM filaments by channeling energy into rotation.

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

This paper uses a three-dimensional fluid model to examine how finite ion temperature affects current-carrying ELM filaments in the tokamak edge. The key finding is that warmer ions create asymmetric electric potentials, driving poloidal flows and rotation that divert energy from radial motion. As a result, filaments merge later even though the total kinetic energy increases from stronger pressure gradients. This shift explains reduced radial transport observed in experiments where ion and electron temperatures are comparable. The work highlights why cold-ion approximations may not suffice for accurate edge plasma modeling.

Core claim

The central discovery is that finite ion temperature substantially alters filament propagation and interaction, resulting in a delay of filament merging despite an increase in total kinetic energy due to a stronger pressure-gradient drive. Examination of single-filament dynamics shows that finite ion temperature generates asymmetric potential structures, strong poloidal flows, and persistent rotational motion, which channel kinetic energy from radial propagation into vortical dynamics. A transition from radially dominated to rotation-dominated behavior occurs as the ion-to-electron temperature ratio increases.

What carries the argument

Asymmetric potential structures in the normalized three-dimensional fluid model that induce poloidal flows and rotational motion in warm-ion filaments.

If this is right

  • Filament merging is delayed with rising ion temperature.
  • Radial transport of filaments decreases as energy is redirected to rotation.
  • Total kinetic energy increases due to enhanced pressure-gradient drive.
  • Single filaments show persistent rotational motion instead of pure radial propagation.
  • A clear transition point exists from radial to rotation-dominated dynamics with increasing ion-to-electron temperature ratio.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Realistic ELM modeling may require ion temperature profiles to predict accurate transport times.
  • The rotational dynamics could influence interactions with background turbulence or magnetic fluctuations not included here.
  • Experimental measurements of filament rotation speeds versus ion temperature could test the predicted transition.
  • Extending the model to bidirectional currents or kinetic ions might reveal additional effects on merging delays.

Load-bearing premise

The chosen initial conditions and normalized three-dimensional fluid model accurately represent realistic ELM filament dynamics in the tokamak edge without needing kinetic corrections.

What would settle it

If increasing the ion temperature in the model does not produce delayed merging or if experiments show faster merging with higher ion temperatures, the claim would be falsified.

Figures

Figures reproduced from arXiv: 2604.03089 by Abhijit Sen, Indranil Bandyopadhyay, Nirmal Bisai, Souvik Mondal.

Figure 1
Figure 1. Figure 1: Time evolution of the density during the merging of two unidirectional current-carrying filaments for (a-d) the cold-ion [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the electrostatic potential during filament merging for (a-d) the cold-ion case ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the separation distance between [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of (left) peak vorticity |ω|max, (middle) total circulation R |ω|, dA, and (right) shear rate |∂vy/∂x| for two interacting current-carrying filaments in the cold-ion (τ = 0.0) and warm-ion (τ = 1.0) regimes. The warm-ion case exhibits significantly enhanced vorticity generation, leading to a strong increase in circulation and shear at later times. This indicates a transition from translation… view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of the density for an isolated filament in (a-d) the cold-ion case ( [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Time evolution of the electrostatic potential for an isolated filament in (a-d) the cold-ion case ( [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Time evolution of the density-weighted center-of-mass (a) poloidal and (b) radial velocities of an isolated filament [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Time evolution of the (a) radial electric field [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Time evolution of the radial, poloidal, and total kinetic energy of an isolated filament for (a) the cold-ion case ( [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Time evolution of the fractional contributions of radial and poloidal kinetic energy to the total kinetic energy for an [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Time evolution of the density-weighted center-of-mass velocity of an isolated filament for different [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Dependence of key driving and energetic quantities on [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Time-averaged radial and poloidal kinetic energy [PITH_FULL_IMAGE:figures/full_fig_p011_13.png] view at source ↗
read the original abstract

Edge-localized-mode (ELM) filaments are crucial for cross-field transport at the tokamak edge; yet, their dynamics are often analyzed using the cold-ion approximation, despite experimental data indicating that Ti~Te . This study employs a normalized three-dimensional fluid model to investigate the influence of finite ion temperature on the dynamics of unidirectional current-carrying ELM-like filaments. We demonstrate that increasing ion temperature substantially alters filament propagation and interaction, resulting in a delay of filament merging despite an increase in total kinetic energy due to a stronger pressure-gradient drive. The examination of single-filament dynamics indicates that finite ion temperature generates asymmetric potential structures, strong poloidal flows, and persistent rotational motion, which channel kinetic energy from radial propagation into vortical dynamics. A comprehensive examination of the ion-to-electron temperature ratio reveals a distinct transition from radially dominated to rotation-dominated behavior as ion temperature increases. These results provide a unified physical explanation for reduced radial transport and delayed merging in the warm-ion domain, emphasizing the necessity of incorporating ion temperature effects in the modeling of ELM filament dynamics and edge plasma transport.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript employs a normalized three-dimensional fluid model to study the dynamics of unidirectional current-carrying ELM-like filaments, focusing on the effects of finite ion temperature (Ti/Te ratio). It claims that increasing ion temperature delays filament merging by generating asymmetric potentials, strong poloidal flows, and persistent rotational motion that channels kinetic energy from radial propagation into vortical dynamics, despite an overall increase in total kinetic energy from a stronger pressure-gradient drive. A transition from radially dominated to rotation-dominated behavior is reported as the ion-to-electron temperature ratio rises.

Significance. If the central claims hold under broader conditions, the work would provide a mechanistic explanation for reduced radial transport in warm-ion ELM filaments and highlight limitations of the cold-ion approximation commonly used in edge plasma modeling. The direct numerical simulation approach yields clear physical insights into energy channeling and flow structures, which is a positive aspect of the study.

major comments (2)
  1. [Results section on multi-filament merging and Ti/Te scan] The reported delay in filament merging and the transition to rotation-dominated dynamics are shown only for a single choice of initial current and density profiles. No systematic scans of filament width, peak current density, or background gradients are presented, even though these directly set the pressure-gradient drive and resulting E×B flows; this leaves the headline claim conditional on the particular initialization rather than a general property of finite-Ti filaments.
  2. [Single-filament dynamics subsection] The analysis of single-filament dynamics asserts that finite ion temperature channels kinetic energy from radial to vortical motion, but provides no explicit decomposition or error bounds on the radial versus poloidal kinetic energy components, nor any grid-convergence or normalization-sensitivity tests to support the robustness of the reported transition.
minor comments (2)
  1. [Abstract] The abstract states that a 'comprehensive examination' of the ion-to-electron temperature ratio was performed; the manuscript should specify the exact range, number of values, and any convergence criteria used in that scan.
  2. [Model description] Clarify the precise normalization chosen for the fluid equations and the rationale for the unidirectional current-carrying initial conditions, including any assumptions about background profiles.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments. We address the major comments point by point below, indicating the revisions we will make to the manuscript.

read point-by-point responses
  1. Referee: The reported delay in filament merging and the transition to rotation-dominated dynamics are shown only for a single choice of initial current and density profiles. No systematic scans of filament width, peak current density, or background gradients are presented, even though these directly set the pressure-gradient drive and resulting E×B flows; this leaves the headline claim conditional on the particular initialization rather than a general property of finite-Ti filaments.

    Authors: The initial profiles were chosen to represent typical ELM filament conditions observed in tokamak experiments. The Ti/Te scan varies the effective drive strength, supporting the generality of the transition. However, to address this concern, we will revise the manuscript to include a discussion on parameter sensitivity and add results from simulations with a different filament width to demonstrate robustness. revision: partial

  2. Referee: The analysis of single-filament dynamics asserts that finite ion temperature channels kinetic energy from radial to vortical motion, but provides no explicit decomposition or error bounds on the radial versus poloidal kinetic energy components, nor any grid-convergence or normalization-sensitivity tests to support the robustness of the reported transition.

    Authors: We agree that an explicit decomposition would enhance the clarity of the results. In the revised manuscript, we will add a figure showing the time evolution of radial and poloidal kinetic energy components with error bounds derived from the simulation data. We will also include a statement on grid convergence, noting that the employed resolution has been tested for convergence in similar setups. revision: yes

Circularity Check

0 steps flagged

No circularity; results from direct numerical integration of fluid equations

full rationale

The paper reports outcomes from solving a normalized 3D fluid model initialized with chosen filament profiles. Claims of delayed merging and transition to rotation-dominated dynamics are computed results, not parameters fitted to target data or quantities defined in terms of themselves. Standard normalizations and initial conditions influence thresholds but do not create self-referential reductions. No load-bearing self-citations, ansatzes, or uniqueness theorems are invoked to force the central findings. This is the expected non-circular outcome for a simulation study.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on standard plasma fluid closure assumptions plus the specific choice of normalization and initial filament current; no new particles or forces are introduced.

free parameters (1)
  • ion-to-electron temperature ratio
    Varied parametrically to demonstrate the transition; specific values are simulation inputs rather than derived.
axioms (1)
  • domain assumption Fluid approximation remains valid for ELM filament scales and velocities
    Invoked by choice of normalized 3D fluid model without kinetic terms.

pith-pipeline@v0.9.0 · 5510 in / 1264 out tokens · 60103 ms · 2026-05-13T19:04:16.197108+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    Intermittentparticletrans- portintwo-dimensionaledgeturbulence.Physics of Plas- mas, 5(12):4214–4228, December 1998

    Y.SarazinandPh.Ghendrih. Intermittentparticletrans- portintwo-dimensionaledgeturbulence.Physics of Plas- mas, 5(12):4214–4228, December 1998

  2. [2]

    M. V. Umansky, S. I. Krasheninnikov, B. LaBombard, and J. L. Terry. Comments on particle and energy bal- ance in the edge plasma of Alcator C-Mod.Physics of Plasmas, 5(9):3373–3376, 09 1998

  3. [3]

    Krasheninnikov

    S.I. Krasheninnikov. On scrape off layer plasma trans- port.Physics Letters A, 283(5):368–370, 2001

  4. [4]

    S. J. Zweben, D. P. Stotler, J. L. Terry, B. LaBombard, M. Greenwald, M. Muterspaugh, C. S. Pitcher, K. Hal- latschek, R. J. Maqueda, B. Rogers, J. L. Lowrance, V. J. Mastrocola, and G. F. Renda. Edge turbulence imag- ing in the alcator c-mod tokamak.Physics of Plasmas, 9(5):1981–1989, 2002

  5. [5]

    J. L. Terry, S. J. Zweben, K. Hallatschek, B. LaBom- bard, R. J. Maqueda, B. Bai, C. J. Boswell, M. Green- wald, D. Kopon, W. M. Nevins, C. S. Pitcher, B. N. Rogers, D. P. Stotler, and X. Q. Xu. Observations of the turbulence in the scrape-off-layer of Alcator C-Mod and comparisons with simulation.Physics of Plasmas, 10(5):1739–1747, 05 2003

  6. [6]

    Simulation of plasma transport by coherent structures in scrape-off layer tokamak plasmas.Physics of Plasmas, 11(8):4018–4024, 2004

    Nirmal Bisai, Amita Das, Shishir Deshpande, Ratnesh- war Jha, Predhiman Kaw, Abhijit Sen, and Raghven- dra Singh. Simulation of plasma transport by coherent structures in scrape-off layer tokamak plasmas.Physics of Plasmas, 11(8):4018–4024, 2004

  7. [7]

    Simulation of plasma transport by coherent structures in scrape-off layer tokamak plasmas.Physics of Plasmas, 11(8):4018–4024, August 2004

    Nirmal Bisai, Amita Das, Shishir Deshpande, Ratnesh- war Jha, Predhiman Kaw, Abhijit Sen, and Raghven- dra Singh. Simulation of plasma transport by coherent structures in scrape-off layer tokamak plasmas.Physics of Plasmas, 11(8):4018–4024, August 2004

  8. [8]

    Bisai, A

    N. Bisai, A. Das, S. Deshpande, R. Jha, P. Kaw, A. Sen, and R. Singh. Formation of a density blob and its dy- namics in the edge and the scrape-off layer of a tokamak plasma.Physics of Plasmas, 12(10):102515, 2005

  9. [9]

    Edge localized modes (elms).Plasma Physics and Controlled Fusion, 38(2):105, feb 1996

    H Zohm. Edge localized modes (elms).Plasma Physics and Controlled Fusion, 38(2):105, feb 1996

  10. [10]

    Wang, G.S

    L. Wang, G.S. Xu, H.Y. Guo, H.Q. Wang, S.C. Liu, K.F. Gan, X.Z. Gong, Y. Liang, N. Yan, L. Chen, J.B. Liu, W. Zhang, R. Chen, L.M. Shao, H. Xiong, J.P. Qian, B. Shen, G.J. Liu, R. Ding, X.J. Zhang, C.M. Qin, S. Ding, L.Y. Xiang, G.H. Hu, Z.W. Wu, G.N. Luo, J.L. Chen, L.Q. Hu, X. Gao, B.N. Wan, J.G. Li, and the EAST Team. Characterizations of power loads o...

  11. [11]

    T. Kass, S. Gunter, M. Maraschek, W. Suttrop, H. Zohm, and ASDEX Upgrade Team. Characteristics of type i and type iii elm precursors in asdex upgrade.Nuclear Fusion, 38(1):111, jan 1998

  12. [12]

    Hironori Takahashi, E. D. Fredrickson, and M. J. Schaf- fer. Scrape-off-layer current model for filament structure observed during edge-localized modes in the diii-d toka- mak.Phys. Rev. Lett., 100:205001, May 2008

  13. [13]

    Spolaore, N

    M. Spolaore, N. Vianello, M. Agostini, R. Cavazzana, E. Martines, P. Scarin, G. Serianni, E. Spada, M. Zuin, and V. Antoni. Direct measurement of current filament structuresinamagnetic-confinementfusiondevice.Phys. Rev. Lett., 102:165001, Apr 2009

  14. [14]

    O. E. Garcia, N. H. Bian, and W. Fundamenski. Ra- dial interchange motions of plasma filaments.Physics of Plasmas, 13(8):082309, 08 2006

  15. [15]

    D. A. D’Ippolito, J. R. Myra, and S. J. Zweben. Con- vective transport by intermittent blob-filaments: Com- parison of theory and experiment.Physics of Plasmas, 18(6):060501, June 2011

  16. [16]

    S. I. KRASHENINNIKOVK, D. A. D’IPPOLITO, and J. R. MYRA. Recent theoretical progress in understand- ing coherent structures in edge and sol turbulence.Jour- nal of Plasma Physics, 74(5), 2008

  17. [17]

    J. R. Myra. Current carrying blob filaments and edge-localized-mode dynamics.Physics of Plasmas, 14(10):102314, 10 2007

  18. [18]

    Bisai, Abhijit Sen, and Indranil 13 Bandyopadhyay

    Souvik Mondal, N. Bisai, Abhijit Sen, and Indranil 13 Bandyopadhyay. Merging dynamics of unidirectional cur- rent carrying filamentary plasma blobs in the edge region of a tokamak.Physics of Plasmas, 32(3):032503, 03 2025

  19. [19]

    Ion temperature fluctuations in the asdex upgrade scrape-off layer.Plasma Physics and Controlled Fusion, 54(8):085009, jun 2012

    M Kocan, F P Gennrich, A Kendl, H W Muller, and the ASDEX Upgrade Team. Ion temperature fluctuations in the asdex upgrade scrape-off layer.Plasma Physics and Controlled Fusion, 54(8):085009, jun 2012

  20. [20]

    P. Manz, D. Carralero, G. Birkenmeier, H. W. Muller, S. H. Muller, G. Fuchert, B. D. Scott, and U. Stroth. Filament velocity scaling laws for warm ions.Physics of Plasmas, 20(10):102307, 10 2013

  21. [21]

    The influ- ence of finite ion temperature on plasma blob dynamics

    P Manz, G Birkenmeier, D Carralero, G Fuchert, H W Muller, S H Muller, B D Scott, U Stroth, T T Ribeiro, E Wolfrum, and the ASDEX Upgrade Team. The influ- ence of finite ion temperature on plasma blob dynamics. Plasma Physics and Controlled Fusion, 57(1):014012, nov 2014

  22. [22]

    Jovanovic, P

    D. Jovanovic, P. K. Shukla, and F. Pegoraro. Effects of the parallel electron dynamics and finite ion temperature on the plasma blob propagation in the scrape-off layer. Physics of Plasmas, 15(11):112305, 11 2008

  23. [23]

    Garcia, Jeppe StÊrk Larsen, Volker Naulin, Anders H

    Jens Madsen, Odd E. Garcia, Jeppe StÊrk Larsen, Volker Naulin, Anders H. Nielsen, and Jens Juul Ras- mussen. The influence of finite larmor radius effects on the radial interchange motions of plasma filaments. Physics of Plasmas, 18(11):112504, 11 2011

  24. [24]

    Wiesenberger, J

    M. Wiesenberger, J. Madsen, and A. Kendl. Radial con- vection of finite ion temperature, high amplitude plasma blobs.Physics of Plasmas, 21(9):092301, 09 2014

  25. [25]

    M. Held, M. Wiesenberger, J. Madsen, and A. Kendl. The influence of temperature dynamics and dynamic fi- nite ion larmor radius effects on seeded high amplitude plasma blobs.Nuclear Fusion, 56(12):126005, sep 2016

  26. [26]

    Halpern, Sébastien Jolliet, Joaquim Loizu, and Paolo Ricci

    Annamaria Mosetto, Federico D. Halpern, Sébastien Jolliet, Joaquim Loizu, and Paolo Ricci. Finite ion tem- perature effects on scrape-off layer turbulence.Physics of Plasmas, 22(1):012308, 01 2015

  27. [27]

    Storm modelling of scrape-off layer filament behaviour with hot ions.Plasma Physics and Controlled Fusion, 67(4):045035, apr 2025

    S Ahmed, J T Omotani, S L Newton, J Birch, and F Militello. Storm modelling of scrape-off layer filament behaviour with hot ions.Plasma Physics and Controlled Fusion, 67(4):045035, apr 2025

  28. [28]

    Transport processes in a plasma.Reviews of plasma physics, 1:205, 1965

    SI Braginskii. Transport processes in a plasma.Reviews of plasma physics, 1:205, 1965

  29. [29]

    Angus, Maxim V

    Wonjae Lee, J.R. Angus, Maxim V. Umansky, and Sergei I. Krasheninnikov. Electromagnetic effects on plasma blob-filament transport.Journal of Nuclear Ma- terials, 463:765–768, August 2015

  30. [30]

    A. A. Stepanenko. Impact of electromagnetic effects on macroscopic dynamics of blobs in hot low-edge plasma of fusion devices.Physics of Plasmas, 27(9):092301, September 2020

  31. [31]

    Umansky, J

    Wonjae Lee, Maxim V. Umansky, J. R. Angus, and Sergei I. Krasheninnikov. Electromagnetic effects on dy- namics of high-beta filamentary structures.Physics of Plasmas, 22(1):012505, 01 2015

  32. [32]

    Dudson, M.V

    B.D. Dudson, M.V. Umansky, X.Q. Xu, P.B. Snyder, and H.R. Wilson. Bout++: A framework for parallel plasma fluid simulations.Computer Physics Communications, 180(9):1467–1480, 2009