Lattice Boltzmann method for relativistic hydrodynamics: Issues on conservation law of particle number and discontinuities

In this paper, we aim to address several important issues about the recently developed lattice Boltzmann (LB) model for relativistic hydrodynamics [M. Mendoza et al., Phys. Rev. Lett. 105, 014502 (2010); Phys. Rev. D 82, 105008 (2010)]. First, we study the conservation law of particle number in the relativistic LB model. Through the Chapman-Enskog analysis, it is shown that in the relativistic LB model the conservation equation of particle number is a convection-diffusion equation rather than a continuity equation, which makes the evolution of particle number dependent on the relaxation time. Furthermore, we investigate the origin of the discontinuities appeared in the relativistic problems with high viscosities, which were reported in a recent study [D. Hupp et al., Phys. Rev. D 84, 125015 (2011)]. A multiple-relaxation-time (MRT) relativistic LB model is presented to examine the influences of different relaxation times on the discontinuities. Numerical experiments show the discontinuities can be eliminated by setting the relaxation time $\tau_e$ (related to the bulk viscosity) to be sufficiently smaller than the relaxation time $\tau_v$ (related to the shear viscosity). Meanwhile, it is found that the relaxation time $\tau_\varepsilon$, which has no effect on the conservation equations at the Navier-Stokes level, will affect the numerical accuracy of the relativistic LB model. Moreover, the accuracy of the relativistic LB model for simulating moderately relativistic problems is also investigated.Comment: 7 figure

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Li et al., Phys. Rev. E. 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/sqrt(a). Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.Comment: 9 figures, 2 tables, accepted by Physical Review E (in press