Lattice Boltzmann Model for Numerical Relativity

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.

9 pages, 11 figures

References in corpus (6)

• The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics
• Conformal and covariant formulation of the Z4 system with constraint-violation damping
• Solving Einstein's Equations With Dual Coordinate Frames
• Implementation of standard testbeds for numerical relativity
• Kinetic Formulation of the Kohn-Sham Equations for ab initio Electronic Structure Calculations
• Modelling the influence of photospheric turbulence on solar flare statistics

Cited by in corpus (2)

• Lattice Boltzmann Methods and Active Fluids
• Lattice Wigner equation