Droplet impact on a superhydrophobic surface under shear airflow: Lattice Boltzmann simulations and scaling analyses

This study utilizes three-dimensional lattice Boltzmann simulations and scaling analyses to investigate droplet impact on superhydrophobic surfaces under shear airflow, revealing how aerodynamic forces enhance spreading and deflection while establishing refined scaling laws to predict the resulting contact footprint and rebound characteristics.

The Gist

Imagine a tiny water droplet falling onto a surface that is super-repellent, like a lotus leaf or a freshly waxed car. In a calm room, this droplet hits, flattens out like a pancake, and then springs back up like a rubber ball, bouncing straight back into the air.

Now, imagine that same droplet hitting that same surface, but this time, a strong wind is blowing sideways across it. What happens? The paper you shared investigates exactly this scenario using powerful computer simulations.

Here is the story of that research, broken down into simple concepts:

1. The Setup: The "Lotus Leaf" and the "Gale"

The researchers are studying superhydrophobic surfaces (think of them as "magic non-stick surfaces"). They wanted to see what happens when a raindrop hits this surface while a strong wind is blowing past it.

In the real world, this happens everywhere:

• Rain hitting an airplane wing while it's flying.
• Spray from a car washing hitting the side of a building while the wind blows.
• Cooling sprays in factories where fans are running.

Most previous studies only looked at drops hitting surfaces in still air. This study adds the "wind" variable to see how it changes the game.

2. The Method: A Digital Wind Tunnel

Instead of dropping millions of real water drops in a lab (which is messy and hard to control), the team built a virtual 3D wind tunnel inside a supercomputer.

They used a special math technique called the Lattice Boltzmann Method. You can think of this like a video game engine that doesn't just draw a picture of water; it actually simulates how every tiny molecule of water and air pushes and pulls against each other. They even added a "wind shear" effect, meaning the wind blows faster at the top of the drop and slower near the ground, just like real air does.

3. What Happens? The "Sliding Pancake"

When the wind blows, the droplet doesn't just bounce straight up. Here is the choreography:

• The Asymmetric Pancake: When the drop hits, it spreads out. Without wind, it's a perfect circle. With wind, the air pushes against the front of the spreading drop, stretching it into an oval (like a rugby ball).
• The Slide: Because the surface is so slippery, the wind doesn't just push the air; it pushes the drop. The drop starts sliding sideways while it's flattening and while it's shrinking back up.
• The "Thumb" Takeoff: When the drop finally decides to leave the surface, it doesn't pop straight up. The wind pushes the top of the drop forward, giving it a weird, thumb-like shape as it detaches. It flies off at a tilted angle, like a skier jumping off a ramp with a side wind.

4. The Big Discovery: The "Footprint"

The most surprising finding is about the contact area (the "footprint" the drop leaves behind).

• In still air: The drop hits, spreads, and retracts. The footprint is just the size of the biggest circle it made.
• In wind: Because the drop is sliding the whole time, it sweeps out a much larger area. The researchers found that the wind can make the drop's final "footprint" 80% larger than it would be in still air.

Why does this matter?
Think of a drop of water cooling a hot engine part. If the drop just sits there, it cools a small circle. If the wind blows it and makes it slide, it cools a much larger strip of the engine. This is huge for designing better cooling systems or anti-icing strategies for planes.

5. The "Magic Formulas"

The researchers didn't just watch the drops; they did the math to predict exactly what would happen.

• The "Super-Weber" Number: They created a new math formula (a modified "Weber number") that combines the speed of the falling drop and the speed of the wind. This single number tells you exactly how much the drop will spread and slide.
• The Bounce Prediction: They figured out how to predict exactly how high the drop will bounce and at what angle it will fly off. It's like having a crystal ball that tells you, "If a drop hits at this speed with this wind, it will bounce off at a 30-degree angle."

The Bottom Line

This paper is like a guidebook for understanding how water behaves when it's caught between gravity, a slippery surface, and a strong wind.

The Analogy:
Imagine a kid jumping on a trampoline.

• No wind: The kid jumps straight up and lands in the same spot.
• With wind: The kid jumps, but a strong fan blows them sideways. They stretch out, slide across the trampoline mat, and land further away, flying off at a weird angle.

The researchers figured out the exact physics of that "windy trampoline jump." This knowledge helps engineers design better surfaces for airplanes (to stop ice), cars (to stay clean), and industrial machines (to cool down faster).

Technical Summary

Here is a detailed technical summary of the paper "Droplet impact on a superhydrophobic surface under shear airflow: Lattice Boltzmann simulations and scaling analyses."

1. Problem Statement

Droplet impact in airflow environments is a critical phenomenon in natural and industrial processes, including anti-icing, spray cooling, and aircraft de-icing. While extensive research exists on droplet impact on superhydrophobic surfaces (SHPS) in quiescent (still) air, the dynamics of droplets impacting SHPS under shear airflow remain underexplored.

• The Gap: Existing studies often overlook the tangential effects of shear flow. In reality, even short-duration impacts (milliseconds) are significantly influenced by airflow due to the low contact angle hysteresis and high interfacial mobility of SHPS.
• The Challenge: Understanding the coupling between shear airflow, impinging droplets, and solid surfaces is essential for predicting droplet behavior (spreading, sliding, rebound) and optimizing surface designs for functional performance.

2. Methodology

The authors employed a three-dimensional pseudopotential multiphase Lattice Boltzmann Method (LBM) to numerically simulate the complex gas-liquid interactions.

• Numerical Model:

  • Solver: Pseudopotential multiphase LBM.
  • Collision Operator: Non-orthogonal multiple-relaxation-time (NMRT) scheme to ensure stability in high-density ratio systems (water-air).
  • Boundary Conditions:
    • Shear Flow: Implemented a piecewise parabolic velocity profile at the inlet to emulate a boundary layer shear flow.
    • Wettability: Used a geometric wetting boundary scheme with a contact angle hysteresis window (Advancing $\theta_A = 155^\circ$, Receding $\theta_R = 145^\circ$) to accurately capture dynamic wetting and contact line pinning/detachment.
  • Validation: The model was rigorously validated against experimental data for both quiescent impact and shear airflow conditions, demonstrating grid independence and high accuracy in capturing morphological evolution and key metrics (spreading diameter, contact time, restitution).

• Simulation Parameters:

  • Variables: Systematic study of the droplet impact Weber number ($We \approx 10-60$) and airflow Reynolds number ($Re \approx 0-700$).
  • Regime: Focused on the complete rebound regime (avoiding splashing/breakup).
  • Hardware: Simulations were accelerated using GPU-based parallel computing (NVIDIA 4090D), achieving speedups of 10–100x compared to CPUs.

3. Key Contributions

The paper makes three primary theoretical and practical contributions:

1. Development of a Modified Weber Number ($We^*$):

   • Introduced a composite scaling parameter, $We^*$, which incorporates the kinetic energy contribution of the horizontal airflow. This allows for the unification of vertical impact inertia and horizontal aerodynamic forcing.
   • Derived a composite scaling law for the dimensionless contact area ($A^*$) based on $We^*$, which successfully collapses simulation data across various $We$ and $Re$ conditions ($R^2 = 0.983$).

2. Scaling Laws for Sliding and Footprint:

   • Developed a theoretical relation for slide displacement ($L_x^*$) based on aerodynamic force analysis (form drag and skin friction).
   • Formulated a piecewise function to predict the final contact footprint ($A_{final}^*$), distinguishing between regimes of minor sliding (elliptical footprint) and strong sliding (elongated footprint). This model accounts for the continuous horizontal sliding unique to shear airflow.

3. Analytical Models for Rebound Dynamics:

   • Derived a refined power law for the vertical velocity restitution coefficient ($\epsilon_z$) using energy partition analysis at detachment.
   • Formulated the streamwise velocity restitution coefficient ($\epsilon_x$) using a sliding velocity approximation.
   • Combined these to accurately predict the total restitution coefficient and the take-off angle, capturing the deflection caused by aerodynamic shear.

4. Key Results

• Morphological Evolution:

  • Airflow breaks the symmetry of the impact. The contact area deforms into an oval shape during spreading and a "thumb-like" shape during detachment.
  • The droplet exhibits continuous horizontal sliding during both spreading and retraction phases, leading to a significant increase in the final contact footprint (up to 80% larger than in quiescent air at high $Re$).

• Contact-Line Characteristics:

  • Spreading: The streamwise spreading factor increases with $Re$, while the transverse factor remains largely independent of airflow.
  • Slide Displacement: Strongly correlated with $Re$ but weakly dependent on $We$. High airflow velocities drive the droplet out of the maximum spreading region, drastically altering the contact footprint.
  • Contact Time: Surprisingly, the total contact time shows only weak dependence on $Re$ and $We$, remaining relatively constant.

• Rebound Dynamics:

  • Energy Partition: Airflow replenishes kinetic energy (driving horizontal motion) but increases viscous dissipation due to irregular deformation.
  • Restitution Coefficients:
    • $\epsilon_z$ (vertical) increases with $Re$ (due to energy replenishment) but decreases with $We$.
    • $\epsilon_x$ (streamwise) increases with $Re$ but decreases with $We$ (high inertia suppresses relative tangential acceleration).
  • Take-off Angle: The rebound angle is deflected downstream. The proposed scaling laws predict the total restitution coefficient and take-off angle with relative deviations of <15% and <20%, respectively.

5. Significance

• Theoretical Insight: The study clarifies the underlying mechanisms of droplet-airflow-surface interactions, specifically how aerodynamic shear couples with impact inertia to alter wetting dynamics on low-adhesion surfaces.
• Predictive Capability: The derived scaling laws provide a robust, quantitative framework for predicting droplet behavior (spreading, sliding, rebound) under realistic aerodynamic conditions without needing full-scale simulations for every scenario.
• Engineering Application: The findings offer practical guidance for designing surfaces in applications such as anti-icing (where ice accretion depends on contact area and time), spray cooling (heat transfer efficiency), and aerospace (rain erosion and de-icing strategies). The ability to predict the enlarged contact footprint and deflected take-off angle is crucial for optimizing surface engineering in windy environments.