Simulation of high-viscosity-ratio multicomponent fluid flow using a pseudopotential model based on the nonorthogonal central-moments lattice Boltzmann method

In this research, the development of a pseudopotential multicomponent model with the capability of simulating high-viscosity-ratio flows is discussed and examined. The proposed method is developed based on the non-orthogonal central moments model in the lattice Boltzmann method, and the exact difference model (EDM) is used to apply the intercomponent interaction force. In contrast to the original Shan-Chen model, in which the applying force has the viscosity-dependent error term, the error term of this model does not depend on the viscosity. A GPU parallel cuda code has been developed and is employed to study the proposed method. Different cases are considered to evaluate the ability of the model, including the Laplace test, a static droplet, and a two-component concurrent channel flow. Also, wetting and nonwetting relative permeabilities for flows with dynamic viscosity ratios between 0.0002 and 5000 are predicted. Numerical results are compared with those of available analytical solutions. Very good agreement between these results are observed. The model has the capability of simulating multicomponent flows with very low kinematic viscosities of the order of ${10}^{\ensuremath{-}5}$ and dynamic viscosity ratios of up to an order of ${10}^{4}$, which is a much wider range compared with that of existing pseudopotential models. Furthermore, the results showed that the parallel processing on GPU significantly accelerated computations. The present parallel performance evaluation shows that the cuda parallel can achieve about 41 times improvement than the CPU serial implementation. The aforementioned enhancement increases the flexibility of the multicomponent lattice Boltzmann method and its applicability to a broader spectrum of engineering applications.