Pith. sign in

REVIEW 3 major objections 6 minor 24 references

Efficient and robust evaluation of fast particle losses in non-axisymmetric tokamak plasmas

T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read Fast-particle losses in non-axisymmetric tokamak plasmas can be reproduced by a one-dimensional advection-diffusion model whose coefficients come from short orbit-following runs, cutting scan cost a hundredfold.

desk verdict Loss maps are genuinely useful and the advection-diffusion idea is promising, but the coefficients are validated with the same orbit-following code used to make them, and the missing sensitivity analysis leaves the 100x speedup claim undermoored. read the letter →

arxiv 1908.02494 v1 pith:EBEHQIAW submitted 2019-08-07 physics.plasm-ph

classification physics.plasm-ph
keywords fastionsparticlelosseslossmapsadvection-diffusionmodelorbit-followingsimulationsITERmagneticfieldrippleELMcontrolcoils
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

Fast particle losses from a magnetically perturbed tokamak are usually estimated with expensive orbit-following simulations, but this paper argues that the losses can be captured with a much cheaper description. It constructs loss maps in a two-dimensional space of orbit constants — the radial coordinate $\rho'$ and pitch $\xi'$ at the outer mid-plane — and shows that distinct loss channels (first-orbit, ripple-induced, stochastic field-line, perturbed banana) appear as separate regions. It then treats the collisionless radial motion as a one-dimensional advection-diffusion process, evaluates the two transport coefficients with a short orbit-following run, and solves the resulting Fokker-Planck equation to predict losses on longer time scales. In ITER benchmarks the model reproduces the overall loss landscape, and a scan over ELM control coil phases runs one hundred times faster than full slowing-down simulations while preserving the locations of high-loss bands. If the approach holds, fast particle parameter scans and wall-load estimates become much cheaper and tie directly to underlying transport mechanisms.

What carries the argument

The load-bearing object is the loss map: a histogram of lost-particle fraction in $(\rho',\xi')$ space, where $\rho'$ is the normalized poloidal flux at the outer mid-plane crossing and $\xi'$ is the pitch there. Because every orbit's topology is fixed by these two coordinates at fixed energy, loss channels appear as separated regions identifiable with known mechanisms. The second mechanism is a coefficient-estimation recipe for the advection-diffusion equation $$\frac{\partial f}{\partial t}=-\frac{\partial}{\partial\rho'}(Kf)+\frac{\$partial^{2}$}{\partial{\rho'}^2}(Df),$$ using short orbit-following data: Gaussian displacement statistics for confined markers and inverse-Gaussian first-passage times for lost markers. The loss map supplies both the diagnostic that validates the model and the weighting that lets transport coefficients be evaluated only where transport actually occurs.

What would settle it

Follow a collisionless orbit ensemble in a field with strong stochastic-ripple transport and test whether the radial displacement has Gaussian statistics with variance growing linearly in time and whether loss times follow the inverse-Gaussian distribution; if the inferred coefficients depend on the averaging window $N$ or on the toroidal launch phase of otherwise identical markers, the Markovian one-dimensional description is wrong. A concrete target is the paper's own 16-37\% disagreement: a case where that gap grows with simulation time rather than staying bounded would show the model is not a reliable predictor for that regime.

Watch

Extended reading notes

Core claim

The central claim is that collisionless fast ion transport in a three-dimensional perturbed magnetic field can be modelled to good accuracy as a one-dimensional Fokker-Planck process in the radial coordinate $\rho'$ (flux-surface coordinate at the outer mid-plane), with magnetic moment $\mu$ and energy $E$ treated as parameters. The advection coefficient $K(\rho';\mu,E)$ and diffusion coefficient $D(\rho';\mu,E)$ are not derived from first principles but measured from a very short orbit-following simulation, about a millisecond in these test cases, in which markers are followed for only tens of poloidal orbits. Confined markers yield $K$ and $D$ from the mean and variance of the radial displacement after averaging over ten outer mid-plane crossings; lost markers contribute through the inverse-Gaussian distribution of first-passage times. In the benchmark with the most complete magnetic field, the model predicted 2.46 MW of total lost $\alpha$ power against 1.79 MW from the collisionless full-orbit simulation (1.42 MW versus 1.22 MW when only collisionless-time-scale losses are counted), and the $8\times8$ scan over ELM control coil phases reproduced the two high-loss bands and their crossing point at one hundredth of the computational cost.

Load-bearing premise

The load-bearing premise is that collisionless radial transport is Markovian and local, so that a particle's random walk in $\rho'$ is fully described by one advection and one diffusion coefficient at each point; if the motion has memory, depends on toroidal or poloidal phase, or needs more than two coefficients, the model fails.

Editorial extensions

If this is right

  • ELM control coil phase scans that would normally require millions of markers per configuration can be reduced to short coefficient runs plus a cheap one-dimensional solve, making it practical to map the full parameter space of coil phases and currents.
  • Because loss maps can be projected from magnetic field structure alone, one can cross-check an orbit-following result or estimate losses without running a dedicated simulation.
  • Marker initialization can be concentrated on the loss channels identified in the map, so peak power load estimates converge with one or two orders of magnitude fewer markers.
  • Orbit-averaged transport codes can incorporate three-dimensional-field fast-ion transport by adding the measured advection and diffusion coefficients rather than following orbits in full geometry.
  • The collisionless approximation captures most alpha losses at birth energy, so collisions matter mainly near thermal energies where neoclassical transport takes over.

Reading between the lines

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

  • Beyond the paper, the same coefficient-measurement recipe could be applied to runaway electrons or other fast species; the one-dimensional Markov assumption is most plausible for strongly passing populations and would need re-testing there.
  • Beyond the paper, if transport coefficients could be estimated from magnetic-field-based projections rather than short orbit runs, the approach might eventually skip orbit-following entirely, although the field-based loss-map projections shown here are close but not yet accurate enough for that.
  • Beyond the paper, the systematic overestimate of lost power in the benchmark suggests the model may be collapsing several mechanisms into a single effective diffusion coefficient; checking whether the inferred coefficients are independent of the averaging window and simulation time would test whether the Markov description is genuinely valid.
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

3 major / 6 minor

Summary. The paper introduces and demonstrates loss-map techniques for fast-particle loss studies in non-axisymmetric tokamaks, applied mainly to ITER. The loss map is a representation of particle birth location in (rho-prime, xi-prime) space that allows loss channels (first-orbit, ripple-induced, stochastic field-line, perturbed banana) to be identified and connected to magnetic field structure. The paper presents optimized marker initialization for improved convergence of total and peak power loads, magnetic-field-based projections for estimating losses without orbit-following, and an analysis of collisional effects. The central new claim is that collisionless fast-ion transport can be described as a one-dimensional advection-diffusion process with coefficients K(rho',mu,E) and D(rho',mu,E) evaluated from a short (1 ms) orbit-following simulation, and that the resulting model can replace full slowing-down simulations for parameter scans, specifically a 64-case ELM control coil phase scan that is claimed to be 100 times faster than full simulations while reproducing the overall loss landscape.

Significance. If the central claim holds, the advection-diffusion model would be a practically useful tool for fast design scans of fast-ion losses, a problem of direct relevance to ITER and future devices. The loss-map framework itself is a useful interpretive and diagnostic tool: the magnetic-field-based loss projection in Section 5.1 provides a partially independent cross-check (the +PR estimate of 1.24 MW versus 1.18 MW from orbit-following is encouraging), and the optimized marker initialization demonstrably improves peak-load convergence in Figure 7. The paper is also honest in reporting the disagreements (2.46 MW versus 1.79 MW in Section 6.1 and underestimation of high-loss scan points in Section 6.2). However, the validation of the advection-diffusion model is currently self-referential, and the extraction of the transport coefficients uses manually tuned parameters without a sensitivity study.

major comments (3)
  1. [Section 6, Eqs. (13)-(15)] The transport coefficients K and D are evaluated as finite-time estimates from a single 1 ms collisionless simulation with N=10 OMP-crossing averaging, and the manuscript states that choosing N and t is critical and that they were deduced with trial and error. No sensitivity analysis with respect to t, N, or the reflecting boundary location is reported. This is load-bearing because Eqs. (11)-(14) yield genuine Fokker-Planck coefficients only if K and D are independent of the estimation window; if they are not, the advection-diffusion model is a tuned fitting procedure rather than a physics-based transport model. The point is reinforced by Figure 9, which shows that after 1 ms only the first-orbit loss channel is fully developed while the stochastic-ripple, stochastic-field-line, and perturbed-banana channels contribute only through small sub-threshold displacements. The authors should show that K and D converge as t and N are varied, and ideally that the model reproduces the time history of losses, not just the final loss map.
  2. [Section 6.1 and Fig. 10; Section 6.2 and Fig. 11] The benchmark accuracy is quantified only through total lost power: the advection-diffusion model gives 2.46 MW versus 1.79 MW from the collisionless orbit-following simulation (a 37% overestimate), and in the short-time comparison 1.42 MW versus 1.22 MW (a 16% overestimate). In the ECC phase scan the model systematically underestimates the losses in exactly the high-loss cases that a design scan is intended to identify. Since the stated purpose of the model is to find interesting regions in parameter space and provide rough estimates, the paper should report a quantitative accuracy metric for the scan, such as rank correlation or per-case relative error, and demonstrate that the high-loss cases are not systematically missed or suppressed. As it stands, the claim that the model is suitable for fast parameter scans is supported mainly by the qualitative visual agreement of the loss contours.
  3. [Section 6.1 and Section 6.2] The transport coefficients are evaluated with the same orbit-following code (ASCOT5) that provides the reference losses, so the advection-diffusion benchmarks test only whether the reduced 1D model can reproduce the output of the same code; they do not test the physical fidelity of the coefficients. The genuine 1 ms-to-100 ms time extrapolation is a meaningful test of time-independence, but it remains a self-consistency check. To support the physical claim that collisionless fast-ion transport is advection-diffusive, the authors should compare the inferred coefficients with the analytic expectations for at least one channel, e.g., the stochastic-ripple diffusion coefficient of Eq. (7) or the stochastic-field-line estimate of Eq. (9), or benchmark against an independent orbit-following implementation.
minor comments (6)
  1. [Eq. (2)] The definition of the magnetic moment appears to be missing the factor m/2 and should be clarified; as written, mu has units of m^2/s^2 rather than the usual J/T.
  2. [Section 5.1] The text says the stochastic-ripple diffusion coefficient is projected using Eq. (9), but Eq. (9) is the stochastic-field-line diffusion coefficient; the ripple diffusion coefficient is Eq. (7).
  3. [Section 6, first paragraph] The sentence 'comes from Ref. where this was shown' has an empty citation; the reference number is missing.
  4. [Eq. (16)] The inverse-Gaussian first-passage-time formula is rendered ambiguously: the denominator in the exponential should be 2 c_1^2 t, and c_1, c_2 should be defined more explicitly. The current notation is hard to parse.
  5. [Section 6, coefficient evaluation paragraph] The averaging in Eq. (15) is not fully specified: it should be stated explicitly that the rho'_j are OMP-crossing averages for a single marker and that the final K and D are weighted averages over markers in each (rho',mu) bin, including how the weights are defined.
  6. [Throughout] There are numerous typos and misspellings, including 'leves' for 'levels', 'quaranteed' for 'guaranteed', 'ploidal' for 'poloidal', 'extent' for 'extend', and 'absent' for 'absence'. A careful proofreading pass is needed.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the advection-diffusion prediction is a long-time extrapolation from short-trajectory statistics, and the ECC phase scan is out-of-sample.

full rationale

The central claim in Section 6 is that short 1 ms collisionless orbit-following statistics yield Fokker-Planck coefficients K and D (Eqs. 13-14) which, when inserted into Eq. (11), reproduce the 100 ms loss map and enable fast ECC-phase scans. This is not circular by construction: the final 100 ms loss fractions are not used to evaluate K and D; the coefficients come from displacement moments and first-passage times of the short simulation. The benchmark in Section 6.1 reports a genuine mismatch (2.46 MW model vs 1.79 MW collisionless orbit-following), and the 64-case scan in Section 6.2 systematically underpredicts high-loss cases while reproducing only the overall shape. Those discrepancies are exactly what an out-of-sample test should show and are not consistent with a forced fit. The sentence that choosing N and t 'can be deduced with some trial and error by evaluating coefficients with different values and carrying out comparisons' discloses hyperparameter tuning of the smoothing/averaging window; this is a robustness concern, not a circular reduction, because the predicted scan is not statistically forced by those hyperparameters. The magnetic-field-based loss projection in Section 5.1 provides an independent cross-check for one case (1.24 MW vs 1.18 MW) without using the advection-diffusion coefficients. Self-citations, mainly to Ref. [7] for loss-map concepts and perturbed-banana transport, are not load-bearing for the 100x speedup or scan suitability claims. The dangling citation 'Ref. where this was shown to be true for runaway electrons' is an editorial omission, not a circular dependency. Overall, the derivation is self-contained as a reduced-model validation: a short-time surrogate is built and then tested against a full simulation at long times and across new configurations.

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

The central claim rests on several modeling choices: the use of P_phi for orbit categorization in 3D fields, the 1D Fokker-Planck representation of transport, and the neglect of collisions. The averaging parameters N and t are explicitly tuned by trial and error. No new physical entities are introduced; the 'perturbed banana transport' is a newly identified mechanism, not a new particle or force.

free parameters (3)
  • Averaging window N and evaluation time t for transport coefficients = N = 10 subsequent OMP crossings, t = 1 ms
    Section 6: 'Choosing N and t is critical for success, and they can be deduced with some trial and error by evaluating coefficients with different values and carrying out comparisons between the advection-diffusion model and orbit-following simulation.' These are tuned to the benchmark.
  • Reflecting boundary location in rho-prime = Not specified; 'or at rho-prime value below which there is believed to be no significant transport'
    Section 6: the outer boundary is absorbing at rho'=1 and the inner boundary is reflecting at rho'=0 or at a selected rho' value where transport is negligible. The choice is not quantified.
  • Reduced-field temperature and density assumptions = T = 10 keV, n = 0.5 n_G (Greenwald density)
    Section 4.1: these values are assumed for the ICRH reduced-field scenarios. They affect the collision operator and the loss-map results, and the conclusion is explicitly tentative pending a known birth distribution.
assumptions (4)
  • domain assumption P_phi can be used to categorize particle orbit topology even when it is not conserved in non-axisymmetric fields.
    Section 2: 'Even in non-axisymmetric plasma, P_phi can be used to categorize particles with respect to their orbit topology.' This is load-bearing for the loss-map coordinates (rho', xi').
  • ad hoc to paper Collisionless fast ion transport can be represented as a 1D Fokker-Planck process in rho-prime with mu and E as parameters.
    Section 6: 'We begin by assuming that the transport can be modelled as one-dimensional process in (rho', mu, E)-space.' This is the central modeling assumption of the advection-diffusion model.
  • domain assumption Collisions can be neglected for alpha particle losses on the relevant 100 ms timescale.
    Section 5.3 shows collisionless loss maps are close to slowing-down maps for 3.5 MeV alphas, but collisions matter at 500 keV and dominate at 10 keV. The advection-diffusion model uses collisionless coefficients.
  • standard math Standard Fokker-Planck and inverse Gaussian first-passage time statistics.
    Equations (11)-(18) rely on standard stochastic process results for Gaussian spreading and first-passage time distributions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Efficient and robust evaluation of fast particle losses in non-axisymmetric tokamak plasmas." pith.science (2026). https://pith.science/paper/EBEHQIAW

@misc{pith2026190802494,
  author       = {Pith},
  title        = {Pith review of: Efficient and robust evaluation of fast particle losses in non-axisymmetric tokamak plasmas},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EBEHQIAW}},
  note         = {Machine review of arXiv:1908.02494}
}
read the original abstract

We present various techniques that make orbit-following Monte Carlo simulations faster and more reliable when assessing collisionless fast particle losses due to magnetic field perturbations. These techniques are based on identifying various loss channels in constants of motion space using the so-called loss maps. We demonstrate that this allows one to attribute losses quantitatively to different transport mechanisms, increase signal-to-noise ratio when estimating FILD signal and peak power loads, and connect magnetic field structure directly to fast particle losses. Furthermore, we show that collisionless fast particle transport can be treated as an advection-diffusion process where the transport coefficients can be evaluated with the orbit-following method. Applying these techniques has the potential to make orbit-following simulations faster to perform, or to avoid them completely, while making the results more reliable as they become more clearly connected to underlying physics. We demonstrate these techniques for ITER by showing how alpha particle losses are affected by various magnetic field perturbations, estimating ICRH losses in reduced field scenarios, and performing a scan on alpha particle losses as a function of ELM control coil current phases.

Figures

Figures reproduced from arXiv: 1908.02494 by the authors.

Figure 1
Figure 1. Loss channels sketched on (ρ, ξ)-plane. Each colored region corresponds to a different loss or transport mechanism: (orange) open-field-line losses, (yellow) gradient￾drift losses, (blue) stochastic ripple transport, (light blue) ripple-trapping, (green) stochastic field-line transport, (purple) perturbed banana transport. Passing-trapped boundaries are shown with dashed black lines. loss maps. Assuming low-collisio… view at source ↗
Figure 2
Figure 2. Distributions related to marker initialization, optimized for loss calculations, for fusion alphas. (a) Physical particle distribution which, in this case, is the alpha particle birth rate. (b) Marker distribution chosen so that it is uniform in (ρ 0 , ξ0 )-space and only represents the edge population. (c) Weighted marker distribution. All these distributions are given in (R, z, ξ) but here only the Rz-profile is s… view at source ↗
Figure 3
Figure 3. Loss maps for alpha particles during slowing-down in different magnetic configurations. The plot shows particle birth position in (ρ 0 , ξ0 ) space, and the fraction of markers lost from that region and their mean loss time. The color shows the mean loss time and the color lightness is varied according to what is the fraction of particles lost locally; no losses occur on regions that appear white. The black contours… view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Loss maps for 1 MeV hydrogen ions, representing ICRH ions, in (a) full-field, (b) half-field, and (c) third-field scenarios. Orange lines marks the location of particles whose banana tip is at R = R0 where the resonance surface is assumed to be. The meaning of other cu…
Figure 3
Figure 3. Figure 3: These cases cover all transport mechanisms [PITH_FULL_IMAGE:figures/full_fig_p008_3.png]
Figure 5
Figure 5. Figure 5: Comparison between losses evaluated with orbit-following method and loss regimes deduced from magnetic field structure for cases (a) TF (b) +ECC, and (c) +PR. Thick and thin black lines are the 10% and 90% loss contours from orbit￾following calculations. The ripple-tra…
Figure 6
Figure 6. Figure 6: Connection between particles’ birth locations in (ρ 0 , ξ0 )-space and their final location on the wall or the divertor. The wall projection is shown in (i) and the close-up of the divertor in (ii). The four plots on the upper-left corner shows the loss maps for partic…
Figure 7
Figure 7. Figure 7: Convergence of (a) total power lost and (b) peak power load as a function of orbit-following calculation simulation time with different marker initialization procedures. The solid vertical line indicates simulation has used 105 markers while the dashed line is for 104 …
Figure 8
Figure 8. Figure 8: Loss maps for alpha particles of different initial energies at +PR case. Simulations were done twice: once without any collisions and once with just pitch collisions, i.e., with energy collisions disabled. ion collisions begin to dominate. At thermal energies, scatteri…
Figure 9
Figure 9. Figure 9: Loss map for the 1 ms simulation that was used to evaluate the transport coefficients. The meaning of the different curves is the same as in [PITH_FULL_IMAGE:figures/full_fig_p012_9.png]
Figure 10
Figure 10. Figure 10: Transport coefficients and model benchmark using +PR case as a testbed. (a) Advection and (b) diffusion coefficients from an orbit-following simulation. (c) P´eclet number P = KL/D which indicates whether transport is dominated by advection, P > 1 (red regions), or di…
Figure 11
Figure 11. Figure 11: Scan of alpha particle losses in ITER baseline scenario with ECC coil configuration n = 3, I = 45 kAt. (a) Results orbit-following simulations where markers were simulated for the full slowing-down process. (b) Results for advection-diffusion model where the transport…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 18 canonical work pages

  1. [1]

    Varje J, Asunta O, Cavinato M, Gagliardi M, Hirvijoki E, Koskela T, Kurki-Suonio T, Liu Y, Parail V, Saibene G, Sipil¨ a S, Snicker A, S¨ arkim¨ aki K and¨Ak¨ aslompolo S 2016 Nuclear Fusion 56 046014 doi: 10.1088/0029- 5515/56/4/046014

  2. [2]

    Varje J, S¨ arkim¨ aki K, Kontula J, Ollus P, Kurki-Suonio T, Snicker A, Hirvijoki E and ¨Ak¨ aslompolo S In press Journal

  3. [3]

    39th EPS Conference on Plasma Physics p P5 URL http://ocs.ciemat.es/ epsicpp2012pap/pdf/P5.088.pdf

    Akers R, Verwichte E, Martin T, Pinches S and Lake R 2012 GPGPU Monte Carlo calculation of gyro- phase resolved fast ion and n-state resolved neutral deuterium distributions Proc. 39th EPS Conference on Plasma Physics p P5 URL http://ocs.ciemat.es/ epsicpp2012pap/pdf/P5.088.pdf

  4. [4]

    26th IAEA Fusion Energy Conference pp TH/4–1 URL https://nucleus

    Akers R, Colling B, Hess J, Liu Y, Turner A, ¨Ak¨ aslompolo S, Varje J, S¨ arkim¨ aki K, Pinches S and Singh M 2018 High fidelity simulations of fast ion power flux driven by 3d field perturbations on iter Proc. 26th IAEA Fusion Energy Conference pp TH/4–1 URL https://nucleus. iaea.org/sites/fusionportal/Shared%20Documents/ FEC%202016/fec2016-preprints/prepr...

  5. [5]

    Kurki-Suonio T, Asunta O, Hellsten T, Hyn¨ onen V, Johnson T, Koskela T, L¨ onnroth J, Parail V, Roccella M, Saibene G, Salmi A and Sipil¨ a S 2009Nuclear Fusion 49 095001 doi: 10.1088/0029-5515/49/9/095001

  6. [6]

    Tani K, Shinohara K, Oikawa T, Tsutsui H, Miyamoto S, Kusama Y and Sugie T 2011 Nuclear Fusion 52 013012 doi: 10.1088/0029-5515/52/1/013012

  7. [7]

    S¨ arkim¨ aki K, Varje J, B´ ecoulet M, Liu Y and Kurki-Suonio T 2018 Nuclear Fusion 58 076021 doi: 10.1088/1741- 4326/aac393

  8. [8]

    White R B, Goldston R J, Redi M H and Budny R V 1996 Physics of Plasmas 3 3043–3054 doi: 10.1063/1.871641

Show all 24 references
  1. [9]

    Hsu C T and Sigmar D J 1992 Physics of Fluids B: Plasma Physics 4 1492–1505 doi: 10.1063/1.860060

  2. [10]

    Goldston R J, White R B and Boozer A H 1981 Physical Review Letters 47 647–649 doi: 10.1103/phys- revlett.47.647

  3. [11]

    Tobita K, Tani K, Neyatani Y, van Blokland A A E, Miura S, Fujita T, Takeuchi H, Nishitani T, Matsuoka M and Takechi S 1992 Physical Review Letters 69 3060–3063 doi: 10.1103/physrevlett.69.3060

  4. [12]

    Rechester A B and Rosenbluth M N 1978 Physical Review Letters 40 38–41 doi: 10.1103/physrevlett.40.38

  5. [13]

    ¨Ak¨ aslompolo S, Kurki-Suonio T, Asunta O, Cavinato M, Gagliardi M, Hirvijoki E, Saibene G, Sipil¨ a S, Snicker A, S¨ arkim¨ aki K and Varje J 2015Nuclear Fusion 55 093010 doi: 10.1088/0029-5515/55/9/093010

  6. [14]

    Kurki-Suonio T, S¨ arkim¨ aki K,¨Ak¨ aslompolo S, Varje J, Liu Y, Sipil¨ a S, Asunta O, Hirvijoki E, Snicker A, Ter¨ av¨ a J, Cavinato M, Gagliardi M, Parail V and Saibene G 2016 Plasma Physics and Controlled Fusion 59 014013 doi: 10.1088/0741-3335/59/1/014013

  7. [15]

    Kurki-Suonio T, ¨Ak¨ aslompolo S, S¨ arkim¨ aki K, Varje J, Asunta O, Cavinato M, Gagliardi M, Hirvijoki E, Parail V, Saibene G, Sipil¨ a S and Snicker A 2016Nuclear Fusion 56 112024 doi: 10.1088/0029-5515/56/11/112024

  8. [16]

    Parail V, Albanese R, Ambrosino R, Artaud J F, Besseghir K, Cavinato M, Corrigan G, Garcia J, Garzotti L, Gribov Y, Imbeaux F, Koechl F, Labate C, Lister J, Litaudon X, Loarte A, Maget P, Mattei M, McDonald D, Nardon E, Saibene G, Sartori R and Urban J 2013 Nuclear Fusion 53 1...

  9. [17]

    ¨Ak¨ aslompolo S, Asunta O, Bergmans T, Gagliardi M, Galabert J, Hirvijoki E, Kurki-Suonio T, Sipil¨ a S, Snicker A and S¨ arkim¨ aki K 2015Fusion Engineering and Design 98-99 1039–1043 doi: 10.1016/j.fusengdes.2015.05.038

  10. [18]

    Liu Y, ¨Ak¨ aslompolo S, Cavinato M, Koechl F, Kurki- Suonio T, Li L, Parail V, Saibene G, S¨ arkim¨ aki K, Sipil¨ a S and Varje J 2016Nuclear Fusion 56 066001 doi: 10.1088/0029-5515/56/6/066001

  11. [19]

    Illustration of the toroidal magnetic field ripple in different magnetic configurations

    Evans T, Orlov D, Wingen A, Wu W, Loarte A, Casper Efficient and robust evaluation of fast particle losses 16 a) Ripple at OMP separatrix b) Ripple well c) Ripple magnitude TF +FI +TBM Figure A1. Illustration of the toroidal magnetic field ripple in different magnetic configuration...

  12. [20]

    2018 ITER Research Plan within the Staged Approach (Level III - Provisional Version) Tech. rep. ITER Organi- zation URL https://www.iter.org/technical-reports

  13. [21]

    Shinohara K, Oikawa T, Urano H, Oyama N, Lon- nroth J, Saibene G, Parail V and Kamada Y 2009 Fusion Engineering and Design 84 24–32 doi: 10.1016/j.fusengdes.2008.08.040

  14. [22]

    Kurki-Suonio T, S¨ arkim¨ aki K, Snicker S, Schneider M, Polevoi A and Mitteau R 2018 Beam Ion Performance and Power Loads in the ITER Pre-Fusion Power Operating Scenarios (PFPO) with Reduced Field and Current 27th IAEA Fusion Energy Conference URL https://conferences.iaea.org...

  15. [23]

    Schneider M, POLEVOI A, Kim S, Loarte A, ARTAUD J, Beaumont B, Bilato R, Boilson D, Campbell D, Dumortier P et al. 2018

  16. [24]

    Illustration of the toroidal magnetic field ripple in reduced field scenarios

    Boozer A H 1981 Physics of Fluids 24 851 doi: 10.1063/1.863445 Efficient and robust evaluation of fast particle losses 18 a) Ripple at OMP separatrix b) Ripple well c) Ripple magnitude Full Half Third Figure A3. Illustration of the toroidal magnetic field ripple in reduced field s...

Pith tools

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