---

In this paper, simultaneous impact of two parallel drops on a thin liquid film is investigated using the lattice Boltzmann method. The purpose of this study is to investigate the effects of surface tension (characterized by Weber number), distance between two drops, and gas kinematic viscosity on the impact. The developed numerical model in this paper which is based on the Shan and Chen pseudo-potential two-phase model makes it possible to access large density ratios, low viscosities, and tunable values of surface tension independent of the density ratio. The model is validated by comparing the coexistence densities with those of Maxwell analytical solution, evaluating the Laplace law for a droplet, and simulating single droplet impact on a thin liquid film. Simulation results of two drops simultaneous impact show that after impact, two jets raised between the drops join each other and form a central jet. Height of this jet increases with time leading to separation of secondary droplets from its tip. When the surface tension value is decreased, the central jet height is increased, but size of the separated droplets is reduced. The crown shape observed in single drop impact is also seen in simultaneous impact of two drops. Increasing distance between two drops leads to a smaller central jet height and an increase in the crown radius. The crown height, however, was found to be independent of the distance. Finally, increasing gas kinematic viscosity reduces the central jet rising speed and delays separation of secondary droplets from the jet.

---

In this paper, a new method is proposed to reach high density ratios and low viscosities based on the Shan-Chen multiphase model in the lattice Boltzmann Method. In this new method the interaction force and as a result the pressure tensor is modified purposefully so that the density of the phases can be adjusted to coincide the corresponding values from the Maxwell equal area rule in thermodynamics. This leads to higher stability and therefore the mentioned purposes are achieved. This new method takes advantage of simplicity and the same implementing procedure in 2D and 3D problems with single or multi relaxation time collision operators. In order to validate the new method, first the coexistence densities of the phases at different subcritical temperatures are compared with those of the Maxwell rule, then the validity of the Laplace law for a droplet is evaluated, after that the spurious velocities around the droplet are evaluated, and finally the broken dam problem is simulated and its results are compared with an experimental data. Results show that the developed model is properly stable and is capable of simulating different multiphase flows at a wide range of density ratios and viscosities.