An analytical-numerical coupled model of liquid droplet impact on solid material surfaces

This study presents an analytical-numerical coupled model that combines a closed-form analytical solution for droplet impact pressure with finite-element simulations of solid response, achieving over 97% computational cost reduction compared to traditional SPH methods while accurately predicting erosion-relevant quantities like peak pressure and impact force.

The Gist

Imagine a raindrop hitting a wind turbine blade at high speed. To the naked eye, it looks like a tiny splash. But to the material of the blade, it's like a microscopic hammer strike. Over thousands of these hits, the blade gets chipped, eroded, and eventually fails. This is a massive problem for wind energy, especially for those giant blades spinning offshore where the tips move incredibly fast.

For decades, engineers have tried to simulate these crashes on computers to design better, tougher blades. But there's a catch: simulating the liquid part (the water drop) is computationally expensive and messy. It's like trying to film a splash in slow motion with a camera that requires a supercomputer just to render a single drop of water.

This paper introduces a clever new way to solve this problem. Here is the breakdown in simple terms:

1. The Old Way: The "Pixelated Splash" (SPH)

Traditionally, to see how a drop hits a solid, scientists use a method called SPH (Smoothed Particle Hydrodynamics).

• The Analogy: Imagine trying to model a water balloon hitting a wall by breaking the water into millions of tiny, invisible marbles. The computer has to track every single marble's position, speed, and collision with its neighbors.
• The Problem: It's incredibly heavy on the computer. It's like trying to count every grain of sand on a beach to predict how a wave will hit the shore. It takes days of computing time, and the results often look "noisy" or jittery, like a bad video game with low graphics settings.

2. The New Way: The "Magic Formula" (ANCM)

The authors of this paper said, "Wait a minute. Do we really need to simulate every single water marble to know how hard the wall gets hit?"

They developed a new method called ANCM (Analytical-Numerical Coupled Model).

• The Analogy: Instead of simulating the water marbles, they used a mathematical shortcut. Think of it like a weather forecaster. You don't need to track every single air molecule to know if it's going to rain; you use a formula based on pressure and temperature to predict the storm.
• How it works: They derived a precise mathematical formula (an "analytical solution") that predicts exactly how the water pressure spreads across the surface of the blade over time. They then fed this "magic formula" directly into the computer model of the solid blade.
• The Result: The computer doesn't waste time simulating the water. It only focuses on how the solid reacts to the pressure.

3. The "Ring of Fire" Discovery

One of the coolest findings in the paper is about where the damage happens.

• The Analogy: When a drop hits a surface, the pressure doesn't just hit the center like a bullseye. It creates a ring-shaped shockwave, like the ripples you see when you drop a stone in a pond, but happening in a split second.
• The Insight: This ring of high pressure is actually what causes the most damage (erosion) on the material. The old, messy computer simulations often missed the clear shape of this ring because of the "noise." The new "Magic Formula" method shows this ring perfectly, like a high-definition photo compared to a blurry sketch.

4. Why This Matters: The 97% Savings

The biggest win here is speed and efficiency.

• The Analogy: If the old method (SPH) was like driving a heavy, gas-guzzling truck to deliver a letter, the new method (ANCM) is like sending a high-speed drone.
• The Stats: The new method is 97% faster and requires 97% less computing power to get the same (or better) answer.
  • Old Way: Takes 100 hours of computer time.
  • New Way: Takes about 3 hours.

Summary

This paper is a game-changer for engineers designing wind turbines, airplane wings, and steam turbines. They realized that to understand how a solid breaks, they don't need to simulate the entire fluid dance in high definition. They just need the "score" (the pressure formula) to tell the solid how to react.

By swapping a heavy, messy simulation for a clean, mathematical prediction, they can now design stronger, longer-lasting materials much faster and cheaper. It's a perfect example of using a little bit of smart math to save a whole lot of computer power.

Technical Summary

Here is a detailed technical summary of the paper "An Analytical-Numerical Coupled Model of Liquid Droplet Impact on Solid Material Surfaces."

1. Problem Statement

Liquid droplet impacts (e.g., rain, sea spray) on solid surfaces, particularly wind turbine blades, cause significant material damage through erosion. This leads to leading-edge erosion (LEE), aerofoil geometry degradation, and reduced turbine lifespan.

• Current Limitations: Traditional numerical approaches for Fluid-Structure Interaction (FSI), such as Smoothed Particle Hydrodynamics (SPH) or Coupled Eulerian-Lagrangian (CEL) methods, are computationally expensive. They require extremely fine mesh resolutions to capture the complex, small-scale interface evolution and large deformations of the liquid phase.
• The Gap: There is a lack of explicit, closed-form analytical solutions that cover the entire duration of droplet impact (from initial contact to spreading) to bridge the gap between early-time self-similar flow and late-time inertia-driven spreading. Furthermore, existing studies rarely couple analytical liquid models directly with numerical solid response analyses to reduce computational costs while maintaining accuracy.

2. Methodology

The authors propose a novel Analytical-Numerical Coupled Method (ANCM) that replaces the explicit simulation of fluid dynamics with an analytical approximation, coupled with Finite Element (FE) analysis for the solid response.

A. Theoretical Framework (Liquid Phase)

The study extends the classical Wagner (1932) theory of inviscid potential flow for a rising expanding disk.

• Baseline: The authors build upon a previous framework (Hao et al., 2026) which provided solutions for early-to-intermediate times but broke down at $t \approx 0.27$ (near peak impact force) due to overestimation of the wet radius and Bernoulli constant.
• Key Modifications for Full-Duration Closure:
  1. Truncation at $r_{max}$: Instead of using the separation point ($r_{sep}$) or the full wet radius ($a(t)$) as the integration boundary, the authors truncate the analytical solution at the pressure-maximum radius ($r_{max}$). This radius corresponds to the location of the ring-shaped pressure distribution observed in experiments. Beyond $r_{max}$, the contribution to the impact force is negligible or overestimated in the potential flow model.
  2. Bernoulli Constant Adjustment: The time-dependent Bernoulli constant $b(t)$, which was previously assumed constant ($1/2$), is set to zero for the entire duration. This corrects the overprediction of pressure at late times when the flow field is no longer infinite.
• Result: These modifications yield an explicit, closed-form analytical solution for pressure $p(r, t)$ and impact force $f(t)$ that remains valid from initial impact through the spreading phase, capturing the transition from self-similar impact to inertia-driven spreading.

B. Numerical Implementation (Solid Phase)

• Coupling Strategy: The derived analytical pressure profile is applied as a time- and spatial-dependent boundary load on the solid surface using a user-defined subroutine (VDLOAD) in the commercial FE software ABAQUS.
• Comparison: The ANCM results are compared against:
  1. Experimental data (Zhang et al., 2019).
  2. Conventional SPH simulations (using ABAQUS) where the liquid is explicitly modeled as particles.

3. Key Contributions

1. Extended Analytical Solution: Derivation of a full-duration, closed-form analytical approximation for droplet impact that bridges the gap between early-time self-similar regimes and late-time spreading, specifically addressing the limitations of previous Wagner-based models.
2. ANCM Framework: Development of a coupled method that eliminates the need for explicit fluid dynamics simulation (SPH/CEL) by using the analytical pressure field as a direct load on the solid.
3. Physical Insight: Demonstration that truncating the flow solution at the pressure-maximum radius ($r_{max}$) effectively suppresses the overestimation of impact loads during the spreading phase, providing a physically meaningful "effective impact area."

4. Results

The study validates the ANCM against SPH simulations and experimental data using a water droplet impacting an aluminum alloy plate.

• Accuracy:
  • Impact Force: The ANCM-predicted force history matches the analytical solution exactly and shows good agreement with experimental data, particularly in capturing the peak force and the decay tail.
  • Pressure & Stress Fields: ANCM successfully reproduces the experimentally observed ring-shaped pressure distribution on the surface. Unlike SPH, which exhibits significant numerical noise (spurious oscillations) due to particle discretization, ANCM provides smooth, high-resolution spatial and temporal distributions of pressure, stress, and displacement.
  • Subsurface Behavior: The model reveals that while surface pressure is ring-shaped, the subsurface stress field transitions to a center-peaked profile, a detail obscured by noise in SPH simulations.
• Grid Independence:
  • SPH simulations required very fine meshes (up to $1.25 \times 10^6$ elements) to converge and reduce oscillations.
  • ANCM achieved grid independence at much coarser resolutions (e.g., $2.56 \times 10^5$ elements) because the loading is analytical and noise-free.
• Computational Efficiency:
  • For the same mesh resolution, ANCM is ~6x faster (16% of the time) than SPH.
  • When accounting for the fact that ANCM requires fewer elements to achieve convergence, the total computational cost is reduced by >97% (ANCM takes only ~2.8% of the time required by SPH).
  • Memory usage is reduced by more than 40%.

5. Significance and Conclusion

This paper presents a paradigm shift in simulating droplet impact erosion:

• Engineering Utility: The ANCM provides a highly efficient tool for analyzing material response (erosion, fatigue, failure) without the prohibitive computational cost of fully coupled fluid-structure simulations.
• Applicability: The method is specifically tailored for inertia-dominated regimes (high Reynolds and Weber numbers), which are typical for wind turbine leading-edge erosion (50–100 m/s impact speeds). While it neglects viscous and capillary effects (valid for the target regime), it captures the critical loading phases responsible for material damage.
• Future Impact: The approach enables reliable numerical predictions for erosion studies on wind turbines, aircraft, and steam turbines using standard computational resources, making high-fidelity FSI analysis accessible for design and maintenance optimization.

In summary, the authors successfully replaced the computationally intensive fluid simulation with a robust analytical approximation, achieving a 97% reduction in computational cost while maintaining or improving the accuracy of solid material response predictions.