  New method to simulate heat transport in multiphase lattice Boltzmann (LB)
method is proposed. The energy transport equation needs to be solved when phase
boundaries are present. Internal energy is represented by an additional set of
distribution functions, which evolve according to a LB-like equation simulating
the transport of a passive scalar. Parasitic heat diffusion near boundaries
with large density gradient is supressed by using the interparticle
"pseudoforces" which prevent the spreading of energy. The compression work and
heat diffusion are calculated by finite differences. The latent heat of a phase
transition is released or absorbed in the inner side of a thin transition layer
between liquid and vapor. This allows one to avoide the interface tracking.
Several tests were carried out concerning all aspects of the processes. It was
shown that the Galilean invariance and the scaling of thermal conduction
process hold as well as the correct dependence of sound speed on the heat
capacity ratio. The method proposed has low scheme diffusion of the internal
energy, and it can be applied to modeling a wide range of multiphase flows with
heat and mass transfer.
