Performance of passive scalar method in heat transfer simulation of a two-dimensional droplet focusing on the parasitic temperatures (Q6058556)

One of complicated problems in mechanics of continuous media is modelling of two-phase flows with heat transfer. Lattice-Boltzmann method (LBM) is an effective numerical approach to solve this problem. But spurious velocities arise on the phase interfaces while using LBM, and that is due to improper equations of state and also with discretization errors. To fix this, passive scalar method may be used, introducing some additional force in equation of energy distribution. In turn, some times (especially while using pseudopotential model) this force may generate parasitic temperatures, and the question arises: how to neglect such negative effect. In the paper the authors compare two approaches to solution of this problem. The first one, that of Li and Luo, is based on elimination of parasitic temperatures by adding the source term to the equation describing temperature distribution. The second one, that of Kuperstokh et al., tries to neglect temperature noise by adding two new terms to temperature distribution function. One of these summands describes temperature changes due to heat conduction and to pressure work. The second summand takes into account the temperature jump that changes internal energy during transfer from the denser phase to the phase with lower density. Several numerical simulations are made to compare two methods mentioned. The initial series of simulations is addressed to verify the entire model: Laplace test of droplet formation, one-phase Rayleigh-Benard convection and the Marangoni effect. After such validation parasitic temperatures on a droplet surface are calculated for different stationary condition. Two-phase Rayleigh-Benard convection is also simulated. Examining simulation results shows that the Kuperstokh model is more exact than that of Li and Luo.