A consistent numerical scheme for self‐gravitating fluid dynamics

In this article, we present a numerical scheme for the 3‐D system of self‐gravitating fluid dynamics in the collisional case as well as in the non‐collisional case. Consistency in the sense of distributions is proved in 1‐D and in absence of pressure. In the other cases consistency is proved under the numerical assumptions of boundedness of the velocity field in the CFL condition and of boundedness of the gradient of the gravitation potential. In 2‐D and 3‐D, concentrations of matter in strings and points can cause a theoretical difficulty in the pressureless case although one observes that the scheme still works. The initial data are L ∞ functions in velocity and L 1 functions in density. Applications are given to numerical simulations of the role of dark matter and gravitational collapse in cosmology as well as Jeans theory. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013