  Starting from the Maxwell-Juettner equilibrium distribution, we develop a
relativistic lattice Boltzmann (LB) algorithm capable of handling
ultrarelativistic systems with flat, but expanding, spacetimes. The algorithm
is validated through simulations of quark-gluon plasma, yielding excellent
agreement with hydrodynamic simulations. The present scheme opens the
possibility of transferring the recognized computational advantages of lattice
kinetic theory to the context of both weakly and ultra-relativistic systems.
