Saved in:
| Main Authors: | , , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2511.04489 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866913028287496192 |
|---|---|
| author | Lappi, Oskar Leggate, Huw Marandet, Yannick Åström, Jan Heljanko, Keijo Borodin, Dmitriy V. |
| author_facet | Lappi, Oskar Leggate, Huw Marandet, Yannick Åström, Jan Heljanko, Keijo Borodin, Dmitriy V. |
| contents | EIRENE [1] is a Monte Carlo neutral transport solver heavily used in the fusion community. EIRENE does not implement domain decomposition, making it impossible to use for simulations where the grid data does not fit on one compute node (see e.g. [2]). This paper presents a domain-decomposed Monte Carlo (DDMC) algorithm implemented in a new open source Monte Carlo code, Eiron. Two parallel algorithms currently used in EIRENE are also implemented in Eiron, and the three algorithms are compared by running strong scaling tests, with DDMC performing better than the other two algorithms in nearly all cases. On the supercomputer Mahti [3], DDMC strong scaling is superlinear for grids that do not fit into an L3 cache slice (4 MiB). The DDMC algorithm is also scaled up to 16384 cores in weak scaling tests, with a weak scaling efficiency of 45% in a high-collisional (heavier compute load) case, and 26% in a low-collisional (lighter compute load) case. We conclude that implementing this domain decomposition algorithm in EIRENE would improve performance and enable simulations that are currently impossible due to memory constraints. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2511_04489 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Scalable Domain-decomposed Monte Carlo Neutral Transport for Nuclear Fusion Lappi, Oskar Leggate, Huw Marandet, Yannick Åström, Jan Heljanko, Keijo Borodin, Dmitriy V. Computational Physics Distributed, Parallel, and Cluster Computing Performance 68W10 (Primary), 68W15, 65C05 (Secondary) D.1.3; G.3; I.6.8; J.2 EIRENE [1] is a Monte Carlo neutral transport solver heavily used in the fusion community. EIRENE does not implement domain decomposition, making it impossible to use for simulations where the grid data does not fit on one compute node (see e.g. [2]). This paper presents a domain-decomposed Monte Carlo (DDMC) algorithm implemented in a new open source Monte Carlo code, Eiron. Two parallel algorithms currently used in EIRENE are also implemented in Eiron, and the three algorithms are compared by running strong scaling tests, with DDMC performing better than the other two algorithms in nearly all cases. On the supercomputer Mahti [3], DDMC strong scaling is superlinear for grids that do not fit into an L3 cache slice (4 MiB). The DDMC algorithm is also scaled up to 16384 cores in weak scaling tests, with a weak scaling efficiency of 45% in a high-collisional (heavier compute load) case, and 26% in a low-collisional (lighter compute load) case. We conclude that implementing this domain decomposition algorithm in EIRENE would improve performance and enable simulations that are currently impossible due to memory constraints. |
| title | Scalable Domain-decomposed Monte Carlo Neutral Transport for Nuclear Fusion |
| topic | Computational Physics Distributed, Parallel, and Cluster Computing Performance 68W10 (Primary), 68W15, 65C05 (Secondary) D.1.3; G.3; I.6.8; J.2 |
| url | https://arxiv.org/abs/2511.04489 |