The pseudopotential lattice Boltzmann model (PP-LBM) is a very popular model for simulating multiphase systems. In this model, phase separation occurs via a short-range attraction between different phases when the interaction potential term is properly chosen. Therefore, the potential term is expected to play a significant role in the model and to affect the accuracy and the stability of the computations. The original PP-LBM suffers from some drawbacks such as being capable of dealing with low density ratios only, thermodynamic inconsistency, and spurious velocities. In this paper, we aim to analyze the PP-LBM with the view to simulate single-component (non-)isothermal multiphase systems at large density ratios and in spite of the presence of spurious velocities. For this purpose, the performance of two popular potential terms and of various implementation schemes for these potential terms is examined. Furthermore, the effects of different parameters (i.e., equation of state, viscosity, etc.) on the simulations are evaluated, and, finally, recommendations for a proper simulation of (non-)isothermal multiphase systems are presented.

To simulate the hydrodynamics and mixing characteristics of chemical reactors by means of a lattice Boltzmann method (LBM), it is essential to consider components with varying molecular weights (and therefore speeds of sound). This option requires modication of the standard equilibrium distribution function and the use of an extended velocity set. In this paper, we show that, for isothermal incompressible single-component non-reactive ows, tuning the speed of sound with a modied equilibrium distribution and an extended velocity set allows for reproducing the proper ow characteristics with strongly reduced errors compared to LBM simulations on standard lattices. This is done for two isothermal benchmarks, viz. a damped standing pressure wave and a decaying viscous Taylor-Green vortex. The convergence as a function of the number of lattice nodes used improves substantially for varying values of the speed of sound.

A surface reaction boundary condition in multicomponent lattice Boltzmann simulations is developed. The method is applied to a test case with nonlinear reaction rates and nonlinear density profiles. The results are compared to the corresponding analytical solution, which shows that the error of the method scales with the square of the lattice spacing.