Stochastic relativistic viscous hydrodynamics from the Metropolis algorithm

by Rajeev Singh (West University of Timișoara)

Monday 16 Jun 2025, 16:00 → 18:00 Europe/Zurich

4/2-037 - TH meeting room (CERN)

We describe a unique formulation of relativistic viscous fluid dynamics based on the density frame. The equations of motion are strictly first order in time and have no non-hydrodynamic modes. The only fields are the energy and momentum densities $T^{0\mu}$ and the only parameters are the equilibrium equation of state and the shear and bulk viscosities. We compare numerical results to QCD kinetic theory simulations for a set of one dimensional tests -- we find excellent agreement within the regime of applicability of hydrodynamics, and the simulation is well behaved even outside of this regime. The results are also compared to a hydrodynamic formulation developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK) as well as to second-order hydrodynamics as implemented in MUSIC. In the density frame we propose a strikingly simple algorithm to simulate stochastic relativistic fluid dynamics based on the Metropolis updates. Each step of the algorithm begins with an update based on ideal hydrodynamics. Then random (three) momentum transfers between fluid cells are proposed. These proposals are then accepted or rejected using the change in entropy as the Boltzmann weight. On average the algorithm reproduces viscous hydrodynamics in the density-frame. The algorithm can be generalized to general coordinates and space-time foliations by parallel transporting the proposed momentum transfers.

 

Based on

https://arxiv.org/abs/2403.04185 

https://arxiv.org/abs/2412.10303

 https://arxiv.org/abs/2412.10306