2024/10/30 by Leonard Storcks, Tobias Buck, Storcks, Leonard +1 · 3 citations
Physics and Astronomy · #Astronomy and Astrophysical Research #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Solar and Space Plasma Dynamics #Stellar, planetary, and galactic studies
paper · pdf · doi:10.48550/arxiv.2410.23093
openalex publication_date 2024/10/30 · openalex created_date 2024/11/14 · openalex updated_date 2026/07/28
We present jf1uids, a one-dimensional fluid solver that can, by virtue of a geometric formulation of the Euler equations, model radially symmetric fluid problems in a conservative manner, i.e., without losing mass or energy. For spherical problems, such as ideal supernova explosions or stellar wind-blown bubble expansions, simulating only along a radial dimension drastically reduces compute and memory demands compared to a full three-dimensional method. This simplification also alleviates constraints on backpropagation through the solver. Written in JAX, jf1uids is a GPU-compatible and fully differentiable simulator. We demonstrate the advantages of this differentiable physics simulator by retrieving the wind's parameters for an adiabatic stellar wind expansion from the final fluid state using gradient descent. As part of a larger "stellar winds, cosmic rays and machine learning" research track, jf1uids serves as a solid foundation to be extended with additional physics modules, foremost cosmic rays and a neural-net powered gas-cooling surrogate and improved by higher order and more accurate numerical schemes. All code is available under https://github.com/leo1200/jf1uids/.