Lattice Boltzmann Model for Numerical Relativity

E. Ilseven,M. Mendoza
DOI: https://doi.org/10.1103/PhysRevE.93.023303
2015-07-28
Abstract:In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
Computational Physics,Cosmology and Nongalactic Astrophysics,General Relativity and Quantum Cosmology,Fluid Dynamics
What problem does this paper attempt to address?