Understanding and predicting the dynamics of complex fluid systems including liquid–liquid phase separation, relevant to both biological and engineered applications, typically uses a nonideal free energy. Introducing such a thermodynamic constraint into the Lattice-Boltzmann Method can be accomplished by altering either the equilibrium distribution function or the external force. The former requires a lengthy parameterization for a free energy of multiple independent variables which becomes cumbersome for more than three components. The latter has been done for a multicomponent compressible system, but a correction term for the force is required to recover the expected conservation equations. This work builds upon the incompressible single component forcing method from He et al. (Journal of Computational Physics, Vol. 152, No. 2, 1999) by deriving and implementing the required force needed to successfully recover the expected mass conservation from a nonideal free energy with an arbitrary number of components. This allows the simulation of more realistic phase separating fluid systems by including many interacting components, which is demonstrated here for up to five components and phases.
Read full abstract