Salvato in:
Dettagli Bibliografici
Autori principali: García-Morillo, José María, Menchiari, Stefano, López-Coto, Rubén
Natura: Preprint
Pubblicazione: 2026
Soggetti:
Accesso online:https://arxiv.org/abs/2604.18703
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866910149980979200
author García-Morillo, José María
Menchiari, Stefano
López-Coto, Rubén
author_facet García-Morillo, José María
Menchiari, Stefano
López-Coto, Rubén
contents In order to model astrophysical environments characterized by radial stratification, such as supernova remnants or expanding superbubbles; correctly understanding the transport of non-thermal particles in astrophysical plasmas is essential. While large-scale Galactic propagation codes exist, they are often optimized for Cartesian or cylindrical geometries and lack the efficiency of one-dimensional spherically symmetric problems. In this work, we present SAETASS (Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry), a novel, open-source numerical tool designed to solve the time-dependent transport equation for astroparticles. The solver is built upon a conservative finite-volume framework that ensures exact particle conservation and numerical stability. To manage the interplay between diverse physical processes, SAETASS employs a modular operator-splitting architecture. Radial advection and continuous momentum losses are treated using a second-order, shock-capturing MUSCL-Hancock scheme, while the diffusive operator is integrated via an implicit, batched Crank-Nicolson algorithm. This approach allows for the robust handling of steep gradients, spatial discontinuities and regularity conditions at the origin. We rigorously validate the code through a suite of tests for pure advection, diffusion and losses. Finally, we demonstrate the solver's capabilities by modelling cosmic-ray proton transport in a real astrophysical scenario. Our results successfully recover established steady-state limits while revealing relevant pre-equilibrium temporal dynamics across Kolmogorov, Kraichnan and Bohm diffusion regimes. SAETASS provides the community with a lightweight, flexible tool for investigating particle acceleration and propagation in complex, radially dependent astrophysical environments.
format Preprint
id arxiv_https___arxiv_org_abs_2604_18703
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle SAETASS: Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry
García-Morillo, José María
Menchiari, Stefano
López-Coto, Rubén
High Energy Astrophysical Phenomena
Instrumentation and Methods for Astrophysics
In order to model astrophysical environments characterized by radial stratification, such as supernova remnants or expanding superbubbles; correctly understanding the transport of non-thermal particles in astrophysical plasmas is essential. While large-scale Galactic propagation codes exist, they are often optimized for Cartesian or cylindrical geometries and lack the efficiency of one-dimensional spherically symmetric problems. In this work, we present SAETASS (Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry), a novel, open-source numerical tool designed to solve the time-dependent transport equation for astroparticles. The solver is built upon a conservative finite-volume framework that ensures exact particle conservation and numerical stability. To manage the interplay between diverse physical processes, SAETASS employs a modular operator-splitting architecture. Radial advection and continuous momentum losses are treated using a second-order, shock-capturing MUSCL-Hancock scheme, while the diffusive operator is integrated via an implicit, batched Crank-Nicolson algorithm. This approach allows for the robust handling of steep gradients, spatial discontinuities and regularity conditions at the origin. We rigorously validate the code through a suite of tests for pure advection, diffusion and losses. Finally, we demonstrate the solver's capabilities by modelling cosmic-ray proton transport in a real astrophysical scenario. Our results successfully recover established steady-state limits while revealing relevant pre-equilibrium temporal dynamics across Kolmogorov, Kraichnan and Bohm diffusion regimes. SAETASS provides the community with a lightweight, flexible tool for investigating particle acceleration and propagation in complex, radially dependent astrophysical environments.
title SAETASS: Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry
topic High Energy Astrophysical Phenomena
Instrumentation and Methods for Astrophysics
url https://arxiv.org/abs/2604.18703