Adaptive thermodynamic consistency control via interface thickness in pseudopotential lattice Boltzmann method across wide temperature ranges

Multiphase flow phenomena are ubiquitous in real-life applications, and the pseudopotential-based lattice Boltzmann multiphase simulation method has gained extensive usage in fields such as fuel cells, energy storage materials, and boiling heat transfer. Over the past two decades, significant improvements have been made to the pseudopotential model. These advancements have greatly enhanced its accuracy in simulating various processes. In this paper, we employ a numerical stability-based unit conversion scheme to ensure an accurate representation of real-world material properties. Additionally, we introduce a more precise non-ideal two-parameter gas state equation Modified Peng-Robinson (MPR) that closely aligns with experimental data, surpassing the commonly used single-parameter Peng-Robinson state equations. Furthermore, we compare the accuracy of the two-state equations for polar and non-polar substances, finding improved accuracy for non-polar substances and a higher degree of fidelity for polar substances when using the MPR equation. We analyze the constant correction coefficients chosen for thermodynamic consistency regulation from the perspectives of vapor-liquid interface thickness and numerical stability, both when held constant and when employing a dynamic correction approach that balances vapor-liquid interface width, numerical stability, and thermodynamic consistency. Finally, we validate to ensure its practical applicability.