{"id":"arxiv/2605.25282","url":"https://arxiv.org/abs/2605.25282","tags":["math.NA","cs.MS","physics.comp-ph","physics.flu-dyn"],"score":0,"title":"Computing statistical solutions of a Mach 2000 astrophysical jet","domain":"math.NA","source":"arxiv","authors":["Stephan Simonis","Gauthier Wissocq"],"reviews":0,"subject":"math.NA","summary":"The multi-dimensional compressible Euler equations admit non-unique entropy solutions in turbulent regimes, and extreme-Mach astrophysical flows are a natural setting in which this breakdown of deterministic well-posedness becomes computationally visible. We compute statistical solutions of a Mach~2000 astrophysical jet, defined as the pushforward of an initial probability measure through a vectorial lattice Boltzmann method, by Monte Carlo sampling with $M=1000$ realizations on grids of up to $12.8$ million cells. Under mesh refinement the individual realizations diverge pathwise, while the…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-05-24","reactions":{}}