Modelling and understanding fast particle transport in non-axisymmetric tokamak plasmas

Thumbnail Image
Files
isbn9789526088211.pdf (5 MB)
Journal Title
Journal ISSN
Volume Title
School of Science | Doctoral thesis (article-based) | Defence date: 2019-11-27
Unless otherwise stated, all rights belong to the author. You may download, display and print this publication for Your own personal use. Commercial use is prohibited.
Author
Date
2019
Department
Teknillisen fysiikan laitos
Department of Applied Physics
Major/Subject
Mcode
Degree programme
Language
en
Pages
98 + app. 67
Series
Aalto University publication series DOCTORAL DISSERTATIONS, 211/2019
Abstract
Understanding fast particle confinement is essential for successful operation of ITER and other advanced fusion experiments. Designing and interpreting experiments require that fast particle transport is studied computationally. The test-particle orbit-following Monte Carlo method is well suited for this purpose as long as the simulations are accompanied with rigorous analysis to ensure the results are valid. This thesis presents a high-performance orbit-following code ASCOT5 for the studies of fast particle transport. The code is verified with respect to known analytical results. Furthermore, this thesis introduces tools to supplement the orbit-following simulations that aid in interpreting the results of the fast particle studies. These tools give credibility to the results and they can also be used to decrease the time required to execute the simulations. ASCOT5 was designed to take full advantage of the modern CPU hardware. The code supports MPI, multithreading, and SIMD instructions resulting in a parallelisation on multiple levels. ASCOT5 uses computational resources more efficiently than its predecessor, but physics-wise they are similar with the exception of the adaptive collision operator. The tools to supplement the orbit-following simulations are demonstrated for ITER. This thesis demonstrates how to construct loss maps which are used to connect fast particle losses to the collisionless transport mechanisms. Constructing loss maps directly from the magnetic field structure provides an alternative way to estimate the fast particle losses without orbit-following simulations. ASCOT5 and the loss-map analysis are used to study fast ion confinement in the presence of various magnetic field perturbations. It was found that the plasma response to the ELM control coils introduces a new loss channel which explains the previously observed shift in fast ion divertor loads. Furthermore, it is shown that the radial transport of runaway electrons and fast ions in an externally perturbed field can be modelled as an advection-diffusion process. This result can be utilized to provide the transport coefficients to the orbit-averaged codes and to perform fast estimates on fast ion losses.

ITERin ja muiden edistyneiden fuusiokoelaitteiden kannalta on oleellista, että nopeiden hiukkasten koossapitoon vaikuttavat seikat ymmärretään. Nopeiden hiukkasten kulkeutumista tulee mallintaa laskennallisin keinoin, jotta tulevia fuusioreaktoreita voidaan suunnitella ja kehittää. Monte Carlo testihiukkas- ja radanseurantamenetelmä sopii hyvin tähän tarkoitukseen. Toisaalta ei ole suoraviivaista varmentaa mallin antamien tulosten oikeellisuutta, jotta voitaisiin poissulkea mahdollisuus koodin tai sen syötteen sisältämille virheille. Työssä esitellään suurteholaskentaan soveltuva radanseurantakoodi ASCOT5, jolla voidaan tutkia nopeiden hiukkasten kulkeutumista. Koodin toimivuus varmennetaan osoittamalla sen toistavan tunnetut analyyttiset tulokset. Tämän lisäksi työssä esitellään menetelmiä, joita voidaan käyttää radanseurantasimulaatioiden tueksi ja tulosten tulkitsemista varten. Näillä menetelmillä lisätään tulosten uskottavuutta sekä lyhennetään simulaatioaikaa. ASCOT5 on suunniteltu erityisesti nykyaikaista CPU-teknologiaa varten. Koodi on rinnakkaistettu usealla tasolla käyttäen MPI-rajapintaa, säikeitä, sekä SIMD-käskykantaa. ASCOT5 on edeltäjäänsä suorituskykyisempi, mutta mallinnetulta fysiikaltaan koodit ovat samankaltaiset. Poikkeuksena on mukautuvan törmäysoperaattorin toteutus. Työssä esitellään radanseurantalaskujen tueksi kehitettyjä menetelmiä ITERille tehdyissä mallinnuksissa. Esittelyssä näytetään miten häviökarttoja rakentamalla voidaan kytkeä nopeiden hiukkasten häviöt törmäyksettömiin kulkeutumisprosesseihin. Häviökarttoja voidaan myös rakentaa suoraan magneettikentän perusteella, jolloin niiden avulla voidaan arvioida häviöitä ilman radanseurantalaskuja. ASCOT5:sta sekä häviökarttoja käytetään työssä nopeiden ionien koossapidon selvittämiseksi, jolloin havaitaan plasman vasteen ELM-kelojen synnyttämään häiriöön aiheuttavan uudenlaisen häviökanavan. Tämä selittää aiemmin havaitun siirtymän diverttorikuormissa. Lisäksi työssä näytetään, että nopeiden ionien sekä karkurielektronien kulkeutumista häirityssä kentässä on mahdollista mallintaa advektio-diffuusio prosessina. Tämä tulos mahdollistaa kulkeutumiskertoimien käytön ratakeskiarvoistetuissa koodeissa sekä nopean tavan arvioida nopeiden ionien häviöitä.
Description
Supervising professor
Groth, Mathias, Prof., Aalto University, Department of Applied Physics, Finland
Thesis advisor
Kurki-Suonio, Taina, Dr., Aalto University, Department of Applied Physics, Finland
Keywords
tokamak, plasma, orbit-following, fast ions, runaway electrons, radanseuranta, nopeat ionit, karkurielektronit
Other note
Parts
  • [Publication 1]: K. Särkimäki, J. Varje, M. Bécoulet, Y. Liu, T. Kurki-Suonio. Mechanics of ELM control coil induced fast particle transport in ITER. Nuclear Fusion, 2018, 58, 7.
    DOI: 10.1088/1741-4326/aac393 View at publisher
  • [Publication 2]: K. Särkimäki, E. Hirvijoki, J. Decker, J. Varje, T. Kurki-Suonio. An advection-diffusion model for cross-field runaway electron transport in perturbed magnetic fields. Plasma Physics and Controlled Fusion, 2016, 58, 12, 125017.
    DOI: 10.1088/0741-3335/58/12/125017 View at publisher
  • [Publication 3]: K. Särkimäki, E. Hirvijoki, J. Terävä. Adaptive time-stepping Monte Carlo integration of Coulomb collisions. Computer Physics Communications, 2018, 222, pp.374-383.
    DOI: 10.1016/j.cpc.2017.09.025 View at publisher
  • [Publication 4]: J. Varje, K. Särkimäki, J. Kontula, P. Ollus, T. Kurki-Suonio, A. Snicker, E. Hirvijoki, S. Äkäslompolo. High-performance orbit-following code ASCOT5 for Monte Carlo simulations in fusion plasmas. Submitted toComputer Physics Communications, Aug. 2019
  • [Publication 5]: K. Särkimäki. Efficient and rigorous evaluation of fast particle losses in non-axisymmetric tokamak plasmas. Submitted to Nuclear Fusion, Aug. 2019
Citation
Permanent link to this item
https://urn.fi/URN:ISBN:978-952-60-8821-1
Collections
[diss] Perustieteiden korkeakoulu / SCI