Force approach for the pseudopotential lattice Boltzmann method

Author(s)

Czelusniak, L. E. ; Mapelli, V. P. ; Guzella, M. S. ; Cabezas-Gomez, L. ; Wagner, Alexander J.

• Document type: Journal article
• Source: PHYSICAL REVIEW E; v. 102, n. 3, p. 19-pg., 2020-09-10.

Abstract

One attractive feature of the original pseudopotential method consists on its simplicity of adding a force dependent on a nearest-neighbor potential function. In order to improve the method, regarding thermodynamic consistency and control of surface tension, different approaches were developed in the literature, such as multirange interactions potential and modified forcing schemes. In this work, a strategy to combine these enhancements with an appropriate interaction force field using only nearest-neighbor interactions is devised, starting from the desired pressure tensor. The final step of our procedure is implementing this external force by using the classical Guo forcing scheme. Numerical tests regarding static and dynamic flow conditions were performed. Static tests showed that current procedure is suitable to control the surface tension and phase densities. Based on thermodynamic principles, it is devised a solution for phase densities in a droplet, which states explicitly dependence on the surface tension and interface curvature. A comparison with numerical results suggest a physical inconsistency in the pseudopotential method. This fact is not commonly discussed in the literature, since most of studies are limited to the Maxwell equal area rule. However, this inconsistency is shown to be dependent on the equation of state (EOS), and its effects can be mitigated by an appropriate choice of Carnahan-Starling EOS parameters. Also, a droplet oscillation test was performed, and the most divergent solution under certain flow conditions deviated 7.5% from the expected analytical result. At the end, a droplet impact test against a solid wall was performed to verify the method stability, and it was possible to reach stable simulation results with density ratio of almost 2400 and Reynolds number of Re = 373. The observed results corroborate that the proposed method is able to replicate the desired macroscopic multiphase behavior. (AU)

• FAPESP's process: 16/09509-1 - Phase change heat transfer processes of high performance applied to solar energy recovery
• Grantee: Gherhardt Ribatski
• Support Opportunities: Research Projects - Thematic Grants
• FAPESP's process: 18/09041-5 - Study of Lattice Boltzmann Method for modeling flows in microchannels of solar absorbers
• Grantee: Luiz Eduardo Czelusniak
• Support Opportunities: Scholarships in Brazil - Doctorate (Direct)