Accelerated Solvers for Neutral Particle Dynamics in Plasma Simulation
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , , |
|---|---|
| Format: | Preprint |
| Publié: |
2025
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866908491113824256 |
|---|---|
| author | Guido, Margherita Kressner, Daniel Mancini, Davide Ricci, Paolo |
| author_facet | Guido, Margherita Kressner, Daniel Mancini, Davide Ricci, Paolo |
| contents | The simulation of turbulence in the boundary region of a tokamak is crucial for understanding and optimizing the performance of fusion reactors. In this work, the use of low-rank linear algebra techniques is shown to enhance the efficiency of boundary simulations, specifically by accelerating the solution of a kinetic model for the neutral particles. Solving the kinetic model deterministically using the method of characteristics requires the solution of integral equations, which typically result in dense linear systems upon discretization. We employ hierarchical matrix approximations to significantly reduce the computational cost of assembling and solving the linear systems, leading to substantial savings in both time and memory. The hierarchical matrix method is implemented and tested within the GBS simulation code for boundary simulations, achieving over 90\% reduction in computation time and memory, and enabling simulations with unprecedented spatial resolution for neutral particles. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_11564 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Accelerated Solvers for Neutral Particle Dynamics in Plasma Simulation Guido, Margherita Kressner, Daniel Mancini, Davide Ricci, Paolo Plasma Physics The simulation of turbulence in the boundary region of a tokamak is crucial for understanding and optimizing the performance of fusion reactors. In this work, the use of low-rank linear algebra techniques is shown to enhance the efficiency of boundary simulations, specifically by accelerating the solution of a kinetic model for the neutral particles. Solving the kinetic model deterministically using the method of characteristics requires the solution of integral equations, which typically result in dense linear systems upon discretization. We employ hierarchical matrix approximations to significantly reduce the computational cost of assembling and solving the linear systems, leading to substantial savings in both time and memory. The hierarchical matrix method is implemented and tested within the GBS simulation code for boundary simulations, achieving over 90\% reduction in computation time and memory, and enabling simulations with unprecedented spatial resolution for neutral particles. |
| title | Accelerated Solvers for Neutral Particle Dynamics in Plasma Simulation |
| topic | Plasma Physics |
| url | https://arxiv.org/abs/2508.11564 |