Color-gradient lattice Boltzmann modeling of immiscible two-phase flows on partially wetting surface

A zero-interfacial-force condition is derived and implemented to improve the wetting boundary scheme for a lattice Boltzmann color-gradient model. This new wetting boundary scheme is validated by two static problems, i.e. a droplet resting on a flat surface and a cylindrical surface, and one dynamic problem, i.e. the capillary filling in a 2 dimensional (2D) channel. In these simulations, we observe that non-physical mass transfer is suppressed and spurious velocities become smaller. Meanwhile, accurate results including dynamic contact line movement are achieved on a broad range of contact angles. The model is then applied to study displacement of immiscible fluids in a 2D channel. Both the displacement velocity and the change rate of finger length are found to exhibit a linear dependence on the contact angle at the viscosity ratio of unity. The displacement velocity decreases but the change rate of finger length increases with increasing capillary number, while the displacement velocity tends to be constant, i.e. two-third of the maximum inlet velocity, at high viscosity ratios or low capillary numbers. In contrast to the displacement velocity, the change rate of finger length is negligible at high viscosity ratios or low capillary numbers, where the finger length is in an equilibrium state, while the equilibrium finger length itself is smaller at a higher viscosity ratio or a lower capillary number

Lattice Boltzmann methods for multiphase flow and phase-change heat transfer

Over the past few decades, tremendous progress has been made in the development of particle-based discrete simulation methods versus the conventional continuum-based methods. In particular, the lattice Boltzmann (LB) method has evolved from a theoretical novelty to a ubiquitous, versatile and powerful computational methodology for both fundamental research and engineering applications. It is a kinetic-based mesoscopic approach that bridges the microscales and macroscales, which offers distinctive advantages in simulation fidelity and computational efficiency. Applications of the LB method are now found in a wide range of disciplines including physics, chemistry, materials, biomedicine and various branches of engineering. The present work provides a comprehensive review of the LB method for thermofluids and energy applications, focusing on multiphase flows, thermal flows and thermal multiphase flows with phase change. The review first covers the theoretical framework of the LB method, revealing certain inconsistencies and defects as well as common features of multiphase and thermal LB models. Recent developments in improving the thermodynamic and hydrodynamic consistency, reducing spurious currents, enhancing the numerical stability, etc., are highlighted. These efforts have put the LB method on a firmer theoretical foundation with enhanced LB models that can achieve larger liquid-gas density ratio, higher Reynolds number and flexible surface tension. Examples of applications are provided in fuel cells and batteries, droplet collision, boiling heat transfer and evaporation, and energy storage. Finally, further developments and future prospect of the LB method are outlined for thermofluids and energy applications