Landau was the first to advance hydrodynamic concepts such as density and velocity to describe the superfluidity of liquid He
4. Due to the recent spectacular success of experiments demonstrating Bose Einstein condensation in dilute Bose atomic gases, interest has been revitalized in the theoretical description of Bose Einstein condensates. Many of the properties of these gases were obtained by using the Gross-Pitaevskii equation (GPE) to derive the hydrodynamic equations for the gases. However, it is interesting to apply the hydrodynamic equations obtained by Yee for bosons. Many of the properties obtained for the dilute Bose gases are also consequences of Yee's hydrodynamic equations, which derive from a formalism distinct from that of the GPE.
We introduce a technique to simulate the dynamics of a thermal dilute Bose gas. With the method of the nonequilibrium statistical operator, a correction to the Gross-Pitaevskii equation (GPE) is derived. This modification of the GPE describes spontaneous scattering processes. It regularizes the ultraviolet catastrophe and captures the dynamics of the particles at all energies, including the tail of the momentum distribution. We apply our model to the simulation of an evaporative cooling quench of a two-dimensional Bose gas.