On the hydrodynamic behaviour of the immersed boundary -- lattice Boltzmann method for wetting problems

This paper evaluates the hydrodynamic behavior and validity limits of an immersed boundary–lattice Boltzmann method for wetting problems by comparing its contact-line model and thin-film formation against boundary element and volume of fluid solvers.

The Gist

Imagine you are watching a drop of water land on a table. It doesn't just sit there; it squishes, spreads out, and eventually settles into a flat puddle. This process, called wetting, is a dance between the liquid, the air, and the solid surface. It sounds simple, but for scientists trying to simulate it on a computer, it's a nightmare.

Why? Because the point where the water, air, and table meet (the "contact line") is incredibly messy. The water tries to stick to the table, but physics says it can't actually touch the table perfectly without causing mathematical chaos.

This paper is about a new way of teaching computers how to handle this messy dance. The authors are testing a specific digital tool called the Immersed Boundary - Lattice Boltzmann (IBLB) method.

Here is the breakdown of their work using simple analogies:

1. The Problem: The "Ghost" Gap

In the past, the IBLB method used a clever trick to stop the water from crashing into the table. Instead of letting the water touch the table, it created a microscopic "ghost" film of air or fluid between the drop and the surface.

Think of it like a hovercraft. The drop doesn't touch the floor; it hovers on a tiny cushion of air.

• The Good News: This prevents the computer math from breaking (crashing).
• The Bad News: Scientists were worried. Does this "hovering" drop behave like a real drop? Does it spread at the right speed? Does it bounce the right way? Or is the ghost film messing up the physics?

2. The Solution: The "Three-Way Race"

To find out if their "hovercraft" method was accurate, the authors didn't just guess. They set up a race between three different computer programs (solvers) to see who could describe the drop's behavior most accurately.

• Runner A (The IBLB Method): The new method with the "ghost film."
• Runner B (The BEM Solver): An old-school, highly precise method that works like a perfectly smooth, frictionless slide. It's great for slow, sticky movements (like honey dripping) but struggles with fast, bouncy ones.
• Runner C (Basilisk): A modern, heavy-duty engine used by engineers. It's like a high-speed camera that can handle fast crashes and splashes (inertia) very well.

3. The Experiments: Slow vs. Fast

The authors ran two main types of tests to see how the runners compared.

Test 1: The Slow Crawl (Low Inertia)
Imagine a drop of honey slowly spreading on a table.

• They compared the IBLB "hovercraft" against the BEM "smooth slide."
• Result: As the drop slowed down, the IBLB method matched the BEM perfectly. The "ghost film" didn't ruin the physics; it just acted like a very thin, invisible layer that didn't change the outcome.

Test 2: The Fast Bounce (High Inertia)
Imagine dropping a water balloon from a height. It hits the table, splats, bounces back up, and then settles. This is chaotic and fast.

• They compared the IBLB method against the heavy-duty Basilisk engine.
• Result: Even in this chaotic, bouncy scenario, the IBLB method got it right! It predicted the exact same "bouncing" shape and the speed at which the drop flattened out as the heavy-duty engine did.

4. The Big Reveal

The most important discovery is that the "ghost film" is harmless.

For a long time, scientists worried that because the IBLB method kept the drop slightly off the ground, it would produce fake results. This paper proves that it doesn't. The method is smart enough to mimic the real physics perfectly, even with that tiny gap.

Why Does This Matter?

Think of this like validating a new recipe for a cake.

• Before, people were worried that using a new type of flour (the IBLB method) might make the cake taste weird or collapse.
• This paper is the taste test. It says, "We compared our cake to the gold-standard cake (BEM) and the professional bakery cake (Basilisk). They all taste the same!"

The Takeaway:
This new computer method is a reliable, versatile tool. It can handle both slow, sticky spreads and fast, bouncy splashes. Because it works so well, scientists can now use it to study more complex things, like:

• How soft gels or pastes spread.
• How elastic materials (like rubbery droplets) interact with surfaces.
• Designing better micro-chips where tiny amounts of liquid need to move precisely.

In short, the authors have proven that their "hovercraft" drop is a real drop in disguise, and we can trust it to solve some of the trickiest fluid problems in engineering and science.

Technical Summary

1. Problem Statement

The paper addresses the hydrodynamic consistency of the Immersed Boundary - Lattice Boltzmann (IBLB) method when applied to wetting dynamics (the spreading of liquid droplets on solid surfaces).

• The Core Issue: The specific IBLB model used (previously proposed in Ref. [25]) employs a "wetting potential" (a disjoining-pressure-like term) to prevent the droplet interface from making direct contact with the solid wall. This regularization prevents abrupt curvature changes but inevitably creates a thin film of fluid between the droplet and the wall.
• The Question: Does the presence of this artificial thin film, which is resolved by only a few lattice nodes, compromise the method's ability to correctly recover the macroscopic hydrodynamic equations (Navier-Stokes) and the physical evolution of the droplet shape, particularly in the contact-line region?
• Previous Limitations: Prior validation (Ref. [25]) focused on static equilibrium shapes and scaling exponents. It lacked a rigorous comparison of the dynamic shape evolution and hydrodynamic behavior against other established solvers, especially in regimes where inertia plays a role.

2. Methodology

The authors conducted a comprehensive benchmarking study comparing their IBLB method against two independent, high-fidelity hydrodynamic solvers:

A. The IBLB Method (Subject of Study)

• Framework: Based on the Lattice Boltzmann method (D3Q19 stencil) coupled with the Immersed Boundary (IB) technique.
• Interface Handling: The droplet interface is represented by a Lagrangian triangular mesh.
• Wetting Model: A specific force density term, $\Pi(d)$, is added to the momentum equation. This term acts as a repulsive/attractive potential based on the distance $d$ from the wall, creating a thin film of thickness $\epsilon$ (where $\epsilon \to 0$ as the regularization length $\xi \to 0$).
• Viscosity: Handles variable viscosity inside and outside the droplet by assigning different relaxation times ($\tau_{LB}$) to fluid nodes.

B. Benchmark Solver 1: Boundary Element Method (BEM)

• Regime: Stokes flow (low inertia, $Re \to 0$).
• Approach: Converts volume partial differential equations into boundary integral equations. It exploits the linearity of the Stokes equations to represent velocity as a superposition of fundamental singular solutions (Green's functions).
• Adaptation: The BEM was adapted to include the same wetting interaction term $\Pi(d)$ used in IBLB to ensure a fair comparison of the contact-line physics.

C. Benchmark Solver 2: Basilisk (Volume of Fluid - VoF)

• Regime: Full Navier-Stokes equations (capable of handling both low and high inertia).
• Approach: A sharp-interface solver using a geometric VoF method on adaptive Cartesian meshes.
• Adaptation: The standard Basilisk "height function" contact angle implementation was discarded. Instead, the authors implemented the same wetting interaction term $\Pi(d)$ (Eq. 10) used in IBLB. This was crucial to replicate the thin-film formation and ensure the comparison tested the hydrodynamic recovery of the specific contact-line model, not just the contact angle imposition.

D. Simulation Scenarios

1. Spreading Dynamics (Low Inertia): Comparison of droplet height evolution ($h(t)$) during spreading from an initial contact angle to a final equilibrium angle.
2. Shear Flow Deformation: Validation of surface tension implementation by simulating a free droplet in shear flow and comparing the deformation index $D$ against analytical solutions.
3. Inertia-Dominated Spreading: Simulation of a droplet impacting/spreading on a wall where inertial forces are significant, leading to complex phenomena like "necking" and bouncing.

3. Key Contributions

• Rigorous Hydrodynamic Validation: The paper moves beyond static shape validation to provide a detailed, dynamic comparison of the IBLB method against BEM and VoF solvers.
• Fair Comparison Protocol: The authors uniquely modified the BEM and Basilisk solvers to incorporate the specific IBLB wetting potential ($\Pi$). This ensures that differences in results are due to hydrodynamic recovery capabilities, not differences in contact-angle boundary conditions.
• Analysis of the Thin Film: The study explicitly investigates whether the thin film (an artifact of the regularization) disrupts the hydrodynamic consistency.
• Inertia Effects: The paper characterizes the behavior of the IBLB method in high-inertia regimes, identifying oscillatory behaviors and "bouncing" effects that were not previously analyzed in depth.

4. Results

• Low-Inertia Regime (BEM Comparison):
  • As the Reynolds number ($Re$) decreases (approaching the Stokes limit), the IBLB simulation results for droplet height evolution ($h(t)$) converge perfectly to the BEM solution.
  • This confirms that the IBLB method correctly recovers the hydrodynamic limit when inertia is negligible, despite the presence of the thin film.
• Surface Tension Consistency (Basilisk Comparison):
  • In shear flow tests, both IBLB and Basilisk matched the analytical solution for droplet deformation ($D$) with high accuracy (relative errors < 3%).
  • This validates that the surface tension implementation in IBLB is physically consistent with established sharp-interface methods.
• Dynamic Spreading (Basilisk Comparison):
  • Shape Evolution: The IBLB and Basilisk solvers produced nearly identical droplet profiles during spreading, with a maximum relative $L_2$-norm error of 1.1%.
  • Inertia-Dominated Regime: In high-$Re$ cases, both solvers captured a "bouncing" effect where the droplet forms a neck near the wall, flattens rapidly, and then relaxes to equilibrium.
  • Height Evolution: The time evolution of the droplet height showed good agreement, with a peak relative error of 9.0% occurring during the transient inertial phase.
• Oscillations: The study noted that at higher Reynolds numbers, IBLB simulations exhibit oscillations in droplet height, which are attributed to inertial effects and are consistent with the physical "bouncing" observed in the Basilisk results.

5. Significance

• Validation of the Contact Line Model: The most critical finding is that the thin film formed by the wetting potential does not compromise the hydrodynamic consistency of the IBLB method. The method correctly recovers the Navier-Stokes dynamics and droplet shape evolution even in the presence of this regularization.
• Beyond Steady-State: Unlike previous work that only validated equilibrium shapes, this paper proves the method is robust during the transient spreading process, including complex inertial dynamics.
• Future Applications: The validated model opens the door to simulating complex wetting scenarios involving:
  • Soft particle pastes and gels.
  • Particles with complex interfacial properties (e.g., elasticity).
  • Data-driven shape prediction tools for microfluidics and materials design.
• Methodological Benchmark: The paper establishes a new standard for validating mesoscale wetting models by requiring direct, one-to-one comparisons of dynamic shapes and hydrodynamic probes against independent solvers, rather than relying solely on scaling laws.

In conclusion, the paper successfully demonstrates that the IBLB method, despite its regularization-induced thin film, is a hydrodynamically consistent and accurate tool for simulating wetting dynamics across both viscous and inertial regimes.