A Corrected Open Boundary Framework for Lattice Boltzmann Immiscible Pseudopotential Models

This paper proposes a corrected open boundary framework based on the Multiple-relaxation-time (MRT) lattice Boltzmann method for immiscible pseudopotential models, which significantly reduces spurious currents and ensures global mass conservation through improved inlet reconstruction, dynamic outlet velocity adjustment, and viscosity-based stability tuning, as validated by four benchmark simulations showing high accuracy in droplet morphology and flow behavior.

The Gist

Imagine you are trying to simulate a busy highway where two different types of cars (let's say, red sedans and blue trucks) are driving together but never mix. They stay in their own lanes, forming distinct groups. In the world of computer science, this is called simulating multiphase flow—like oil and water, or bubbles in soda.

The paper you shared introduces a new, smarter way to run these simulations using a method called the Lattice Boltzmann Method (LBM). Think of LBM as a giant grid of tiny tiles where every tile holds a little bit of "traffic data" (density and speed). The computer updates these tiles step-by-step to see how the traffic moves.

However, the authors found that the old way of handling the entrances and exits of this digital highway was broken. Here is the problem and their solution, explained with simple analogies:

The Three Big Problems (The "Traffic Jams")

1. The "Ghost" Entrance:

   • The Issue: When cars enter the simulation, the old method was like a bouncer who guessed the speed of the incoming cars based on the people standing behind them. Sometimes, the guess was wrong, and the "real" number of cars entering didn't match the "planned" number. This caused the simulation to drift off course.
   • The Fix: The authors added a "Correction Coefficient." Imagine giving the bouncer a clipboard and a calculator. Instead of guessing, he now double-checks the math to ensure the exact number of cars enters at the exact speed required. This makes the entrance perfectly accurate.

2. The "Leaky" Exit:

   • The Issue: When cars leave the simulation, the old method was like a revolving door that sometimes spun too fast or too slow. Over time, this meant cars were either piling up inside the simulation or disappearing entirely. The total number of cars wasn't staying the same (mass wasn't conserved).
   • The Fix: They introduced a "Velocity Correction." Now, the exit door is smart. It constantly counts how many cars entered and adjusts the speed of the exit door to match. If 100 cars entered, it ensures exactly 100 leave, keeping the total count perfectly stable.

3. The "Shaky" Interface (Spurious Currents):

   • The Issue: Where the red sedans and blue trucks meet, there is a boundary. In the old simulations, this boundary would jitter and shake violently, like a shaky camera. These "spurious currents" are fake movements that don't exist in real life but ruin the computer model, making the cars look like they are vibrating apart.
   • The Fix: The authors tweaked a "Relaxation Coefficient" (think of this as a "stability dial"). They realized that if the fluid is thick (viscous), the dial needs to be set differently. By tuning this dial based on how "thick" the fluids are, they smoothed out the shaking. The boundary between the fluids became calm and stable, like a smooth glass surface.

The Results: A Smoother Ride

The authors tested their new "Smart Highway" system with four different scenarios:

• The Soap Bubble Test: Checking if a bubble stays round and stable.
• The Stretch Test: Seeing how a droplet stretches when pulled.
• The Pipe Flow: Watching two fluids flow side-by-side in a pipe.
• The Droplet Factory: Simulating how droplets are made in tiny channels (like in medical devices).

The Outcome:

• The fake shaking (spurious currents) was reduced by 65%.
• The simulation kept the total amount of fluid almost perfectly constant (only a tiny 3.5% error).
• The shapes of the droplets matched real-life experiments and math formulas with less than 5% error.

Why Does This Matter?

Imagine you are designing a new drug delivery system that uses tiny droplets, or a fuel cell that needs to manage water and gas perfectly. If your computer simulation is "leaky" or "shaky," your design might fail in the real world.

This paper provides a more reliable, stable, and accurate toolkit for scientists and engineers. It's like upgrading from a sketchy, hand-drawn map to a high-definition GPS that knows exactly where every car is, ensuring that when you build a real-world device, it works exactly as the computer predicted.

In short: They fixed the front door, the back door, and the shaky fence in between, making the simulation of mixing fluids much more trustworthy.

Technical Summary

1. Problem Statement

The pseudopotential Lattice Boltzmann Method (LBM) is a powerful tool for simulating multiphase flows due to its physical intuitiveness and ability to handle complex interfaces without explicit tracking. However, simulating dynamic multi-component immiscible fluid systems with open boundaries (inlets and outlets) presents three persistent challenges:

• Boundary Accuracy: Existing non-equilibrium (NEQ) extrapolation schemes (e.g., Zou-He or Zhao-Li) often fail to accurately recover macroscopic quantities (density and velocity) at the inlet boundary in 3D or complex scenarios, leading to interpolation errors and reduced simulation accuracy.
• Mass Conservation: Traditional outflow schemes often suffer from global mass non-conservation, particularly during prolonged interfacial disturbances (e.g., droplet generation). This occurs because the mass flow rates at the inlet and outlet become imbalanced, causing the total mass in the domain to drift over time, eventually leading to simulation instability or crashes.
• Spurious Currents: Non-physical parasitic currents (spurious currents) arise at phase interfaces with curvature. In immiscible pseudopotential models, these currents can be excessive, causing unrealistic interface deformation, breakup, or instability, especially when viscosity ratios are high.

2. Methodology

The authors propose a corrected open boundary framework based on the Multi-Relaxation-Time (MRT) collision operator and an 8th-order interaction force format. The methodology introduces three specific improvements:

A. Inlet Boundary Correction (Distribution Function Reconstruction)

• Approach: Building on the standard NEQ extrapolation scheme, the authors introduce a correction coefficient ($\alpha$) to reconstruct the distribution function ($f_i$) at the inlet boundary.
• Mechanism: The unknown distribution functions pointing into the domain are adjusted using $\alpha$ to ensure that the reconstructed macroscopic velocity and density exactly match the prescribed boundary conditions. This eliminates interpolation errors inherent in standard extrapolation methods.
• Ghost Nodes: Two layers of ghost nodes are added at the boundary to calculate interaction forces accurately, reducing anisotropy effects near the wall.

B. Outlet Boundary Correction (Velocity Adjustment for Mass Conservation)

• Approach: A velocity correction coefficient ($\nu$) is introduced to adjust the outlet velocity in real-time.
• Mechanism: The method calculates the mass flow rate at the inlet and the outlet at every time step. The outlet velocity is scaled by $\nu$ such that:
  $$\nu = \frac{\text{Mass Flow}_{inlet}}{\text{Mass Flow}_{outlet}}$$
  This ensures the net mass flow rate in the system is zero, strictly enforcing global mass conservation even during dynamic events like droplet detachment.

C. Spurious Current Suppression (Viscosity-Dependent Relaxation)

• Approach: The relaxation parameter $\tau_e$ (associated with the bulk viscosity/stability) is dynamically adjusted based on the kinematic viscosity of the two-phase fluids.
• Mechanism: By tuning $\tau_e$ according to the viscosity of the more viscous phase, the method improves the isotropy at the phase interface. This significantly attenuates spurious currents without compromising the physical evolution of the flow.

3. Key Contributions

1. Accurate Macroscopic Recovery: The introduction of correction coefficients at the inlet ensures that boundary conditions are met with high precision, overcoming the limitations of standard NEQ schemes in immiscible systems.
2. Robust Mass Conservation: The real-time velocity correction at the outlet solves the critical issue of mass drift in open-boundary multiphase simulations, enabling long-duration simulations of droplet generation and migration.
3. Enhanced Stability: The viscosity-based adjustment of the relaxation parameter suppresses spurious currents by two orders of magnitude, allowing for stable simulations across a wide range of viscosity ratios (up to 150:1).
4. 3D Applicability: The framework is validated in both 2D and 3D, including complex coaxial flow structures, demonstrating scalability for high-fidelity simulations.

4. Results and Validation

The proposed method was validated against four benchmark cases:

• Laplace Test & Taylor Deformation:
  • The model accurately reproduced Laplace's law ($\Delta P = \gamma/R$).
  • Spurious currents were reduced to the order of $10^{-4}$ (a significant improvement over the $10^{-2}$ to $10^{-3}$ seen in uncorrected models).
  • The method successfully simulated Taylor deformation across various viscosity ratios ($M=1$ to $150$) and capillary numbers, confirming that the parameter adjustments do not alter the physical deformation behavior.
• Two-Phase Poiseuille Flow:
  • Numerical velocity profiles matched analytical solutions with relative errors below 1.52% for viscosity ratios of 1/20, 1, and 20.
  • Mass Conservation: The corrected method maintained the total system mass with an average deviation of ~3.5%, whereas the uncorrected method showed significant mass drift and eventual simulation failure.
• Droplet Migration in Microchannels:
  • Simulations of droplets moving through channels (with and without obstacles) showed that the corrected method preserves the original scheme's advantages while ensuring mass conservation.
  • Deviations in droplet position and interface deformation compared to the original scheme were less than 5%.
• Droplet Generation (T-shaped & Co-flow):
  • The method successfully simulated droplet formation in T-shaped and co-flow devices, matching experimental data and previous numerical studies (Shi et al., Chen et al.).
  • In 3D focused coaxial flow simulations, the droplet diameter-to-tube width ratio matched experimental results with an average absolute deviation < 5%.

5. Significance

This work addresses critical bottlenecks in applying pseudopotential LBM to practical engineering problems involving open boundaries.

• Scientific Impact: It provides a theoretically sound and numerically stable framework for simulating complex multiphase dynamics where mass conservation and interface integrity are paramount.
• Practical Application: The method enables high-fidelity simulations of microscale droplet dynamics, which are crucial for applications in:
  • Biomedical Engineering: Lab-on-a-chip devices, drug delivery, and PCR microfluidics.
  • Energy Systems: Fuel cells and boiling heat transfer.
  • Material Science: Synthesis of monodisperse microspheres.
• Efficiency: By reducing spurious currents and preventing mass drift, the method reduces the need for excessive numerical damping or extremely small time steps, making long-term 3D simulations computationally feasible.

In conclusion, the proposed corrected open boundary framework significantly enhances the reliability, accuracy, and stability of immiscible pseudopotential LBM, making it a robust tool for next-generation multiphase flow simulations.