Consistent evaporation formulation for the phase-field lattice Boltzmann method

Thermal lattice Boltzmann model for liquid-vapor phase change

The lattice Boltzmann method is adopted to solve the liquid-vapor phase change problems in this article. By modifying the collision term for the temperature evolution equation, a thermal lattice Boltzmann model is constructed. As compared with previous studies, the most striking feature of the present approach is that it could avoid the calculations of both the Laplacian term of temperature [∇·(κ∇T)] and the gradient term of heat capacitance [∇(ρc_{v})]. In addition, since the present approach adopts a simple linear equilibrium distribution function, it is possible to use the D2Q5 lattice for the two-dimensional cases considered here. Thus, the present model is more efficient than previous models in which the lattice is usually limited to the D2Q9. The proposed model is first validated by the problems of droplet evaporation in open space and droplet evaporation on a heated surface, and the numerical results show good agreement with the analytical results and the finite difference method. Then it is used to model the nucleate boiling problem, and the relationship between detachment bubble diameter and gravitational acceleration obtained with the present approach fits well with previous works.

Droplet evaporation in finite-size systems: Theoretical analysis and mesoscopic modeling

The classical D^{2}-Law states that the square of the droplet diameter decreases linearly with time during its evaporation process, i.e., D^{2}(t)=D_{0}^{2}-Kt, where D_{0} is the droplet initial diameter and K is the evaporation constant. Though the law has been widely verified by experiments, considerable deviations are observed in many cases. In this work, a revised theoretical analysis of the single droplet evaporation in finite-size open systems is presented for both two-dimensional (2D) and 3D cases. Our analysis shows that the classical D^{2}-Law is only applicable for 3D large systems (L≫D_{0}, L is the system size), while significant deviations occur for small (L≤5D_{0}) and/or 2D systems. Theoretical solution for the temperature field is also derived. Moreover, we discuss in detail the proper numerical implementation of droplet evaporation in finite-size open systems by the mesoscopic lattice Boltzmann method (LBM). Taking into consideration shrinkage effects and an adaptive pressure boundary condition, droplet evaporation in finite-size 2D/3D systems with density ratio up to 328 within a wide parameter range (K=[0.003,0.18] in lattice units) is simulated, and remarkable agreement with the theoretical solution is achieved, in contrast to previous simulations. The present work provides insights into realistic droplet evaporation phenomena and their numerical modeling using diffuse-interface methods.