High-order quadrature-based lattice Boltzmann models for the flow of ultrarelativistic rarefied gases

TL;DR

This paper introduces relativistic lattice Boltzmann models using quadrature rules for simulating ultrarelativistic gas flows, demonstrating high accuracy in shock and expansion problems across different regimes.

Contribution

It develops a systematic quadrature-based approach for relativistic lattice Boltzmann models tailored for massless particle flows, optimized for one-dimensional applications.

Findings

Models accurately reproduce exact solutions in inviscid and ballistic limits.

Convergence tests determine quadrature order needed for desired accuracy.

Velocity set size increases in ballistic regime for better profile reproduction.

Abstract

We present a systematic procedure for the construction of relativistic lattice Boltzmann models (R-SLB) appropriate for the simulation of flows of massless particles. Quadrature rules are used for the discretization of the momentum space in spherical coordinates. The models are optimized for one-dimensional flows. The applications considered in this paper are the Sod shock tube and the one-dimensional boost invariant expansion (Bjorken flow). Our models are tested against exact solutions in the inviscid and ballistic limits. At intermediate relaxation times (finite viscosity), we compare with the results obtained using the Boltzmann approach of multiparton scattering model for the Sod shock tube problem, as well as with a semi-analytic solution for the non-ideal Bjorken flow. In all cases our models give remarkably good results. We define a convergence test in order to find the quadrature order necessary to obtain convergence at a predefined accuracy. We show that, while in the hydrodynamic regime the number of velocities is comparable to that required for the more popular collision-streaming type models, as we go towards the ballistic regime, the size of the velocity set must be substantially increased in order to accurately reproduce the analytic profiles.