A numerical study on the coefficient of restitution of wet collisions

En utilisant des simulations SPH validées expérimentalement, cette étude numérique établit une loi d'échelle reliant le coefficient de restitution aux nombres de Stokes et d'épaisseur de film dans le régime de collisions humides à nombre de Weber modéré-élevé, où les effets de tension de surface sont négligeables.

L'essentiel

🌧️ Le Grand Jeu de la Balle Humide : Quand l'eau freine le rebond

Imaginez que vous jouez avec une balle de tennis. Si vous la lancez sur un sol sec, elle rebondit haut et joyeusement. C'est ce qu'on appelle un choc "sec". Mais que se passe-t-il si vous lancez cette même balle sur un sol recouvert d'une fine couche d'eau ? Soudain, la balle semble "coller", elle perd de son énergie et ne remonte pas aussi haut.

C'est exactement ce que les chercheurs de cet article ont voulu comprendre : comment l'eau change le rebond d'un objet solide.

1. La Méthode : Une Simulation Virtuelle de Précision

Au lieu de lancer des milliers de billes en verre dans des bains d'eau (ce qui serait long et coûteux), les scientifiques ont utilisé un super-ordinateur. Ils ont créé un monde virtuel où ils peuvent voir chaque goutte d'eau bouger.

• L'analogie du "Jeu de Lego liquide" : Imaginez que l'eau n'est pas un fluide continu, mais une foule de petites billes invisibles qui se poussent les unes les autres. Les chercheurs ont simulé des milliards de ces interactions pour voir comment l'eau réagit quand une bille de verre tombe dessus. C'est comme si on filmait l'eau à l'intérieur d'une caméra ultra-rapide, mais en 3D et en détail mathématique.

2. Le Problème : Ce n'est pas aussi simple qu'il y paraît

Avant cette étude, les scientifiques pensaient qu'il suffisait de connaître la vitesse de la bille et la viscosité (l'épaisseur collante) de l'eau pour prédire le rebond. Ils utilisaient une formule magique appelée "Nombre de Stokes".

• L'analogie du "Manteau trop grand" : C'était comme si on pensait que la taille d'un manteau dépendait uniquement de la température. Or, cette étude a découvert que la taille de la bille par rapport à l'épaisseur de la couche d'eau est tout aussi importante !
  • Si la bille est énorme et l'eau très fine, l'eau s'écrase comme une fine pellicule de beurre.
  • Si la bille est petite et l'eau relativement épaisse, l'eau se comporte comme une piscine profonde où des tourbillons peuvent se former.

3. La Découverte : Deux Règles pour Deux Mondes

En analysant leurs simulations, les chercheurs ont découvert qu'il n'y a pas une seule règle, mais deux règlements différents selon la situation :

• Règle n°1 (La Balle Géante) : Quand la bille est grosse par rapport à la fine couche d'eau, le rebond dépend fortement de la vitesse et de la viscosité. C'est comme essayer de traverser une flaque d'eau très fine : l'eau résiste, mais elle ne fait pas beaucoup de vagues complexes.

• Règle n°2 (La Balle Petite) : Quand la bille est petite par rapport à l'épaisseur de l'eau, la vitesse compte beaucoup moins. Pourquoi ? Parce que l'eau a assez d'espace pour créer des tourbillons (des petits tourbillons d'eau qui tournent sur eux-mêmes). Ces tourbillons volent l'énergie de la bille, comme un frein à main invisible.

• L'analogie du "Trafic Routier" :

  • Dans le premier cas, c'est comme une voiture sur une route très étroite : elle doit avancer lentement à cause de la friction (l'eau).
  • Dans le deuxième cas, c'est comme une voiture dans un grand parking : elle peut faire des dérapages et des virages serrés (les tourbillons) qui la ralentissent différemment.

4. Pourquoi est-ce utile ?

Comprendre ces règles, c'est comme avoir un manuel d'instructions pour prédire comment les choses bougent dans la nature et l'industrie.

• Dans la vraie vie : Cela aide à concevoir de meilleurs procédés pour l'industrie pharmaceutique (mélanger des poudres humides), pour le transport de sédiments dans les rivières, ou même pour comprendre comment les glissements de terrain (boue) se déplacent.
• Le résultat final : Les chercheurs ont créé une nouvelle "formule mathématique" qui combine la vitesse, la taille de la bille et l'épaisseur de l'eau. Cette formule permet de prédire exactement à quelle hauteur une bille rebondira, qu'elle soit grosse, petite, rapide ou lente.

En résumé

Cette étude nous apprend que l'eau n'est pas un simple frein uniforme. Selon la taille de l'objet qui tombe et l'épaisseur de la couche d'eau, l'eau peut soit résister par friction (comme du miel), soit créer des tourbillons qui volent l'énergie (comme un courant marin). Grâce à des simulations informatiques très avancées, les scientifiques ont enfin pu écrire les deux règles du jeu pour prédire ces rebonds humides.

Résumé technique

1. Problématique et Contexte

L'estimation de la dissipation d'énergie lors de collisions de particules humides est cruciale pour modéliser divers phénomènes d'écoulement granulaire (glissements de terrain, mélange acoustique, granulation humide, transport sédimentaire).

• Limites des modèles actuels : Pour les collisions sèches, le coefficient de restitution (COR) est souvent traité comme une constante matérielle. Cependant, la présence d'un film liquide interstitiel introduit des mécanismes de dissipation supplémentaires (écoulement visqueux, accumulation de pression inertielle, déformation de la surface libre).
• Hypothèse défaillante : La plupart des modèles macroscopiques actuels supposent que le COR humide est gouverné uniquement par le nombre de Stokes (St), qui compare l'inertie de la particule à la résistance visqueuse.
• Constat expérimental précédent : Des études antérieures (notamment Gollwitzer et al.) ont montré que le COR augmente avec la vitesse d'impact avant de saturer, mais elles n'ont pas exploré systématiquement l'effet indépendant de l'épaisseur du film liquide par rapport au diamètre de la particule.
• Objectif de l'étude : Combler ce manque en déterminant si le nombre de Stokes seul suffit à caractériser la dynamique de rebond, ou si l'épaisseur relative du film joue un rôle déterminant, en particulier dans le régime de nombre de Weber modéré à élevé (où les effets inertiels dominent la tension de surface).

2. Méthodologie

Les auteurs ont utilisé une approche numérique basée sur l'Hydrodynamique des Partules Lissées Incompressible (ISPH) couplée à la dynamique des corps rigides.

• Modélisation ISPH :

  • Utilisation d'une formulation basée sur la projection pour garantir des champs de pression stables et physiquement cohérents, essentiels pour capturer les pertes d'énergie.
  • Intégration temporelle via le schéma de Cummins et modèle de viscosité de Cleary–Monaghan.
  • Résolution de l'équation de Poisson pour la pression via la méthode BiCGSTAB avec préconditionnement de Jacobi.
  • Négligence de la tension de surface (nombre de Weber $O(10) - O(10^2)$), car les conditions d'impact correspondent à un régime où les effets inertiels prédominent.

• Modélisation du Corps Rigide :

  • La bille est modélisée comme un corps rigide libre composé d'un nuage de particules lagrangiennes.
  • Optimisation : Pour réduire les coûts de calcul, seule une coque sphérique tronquée (la partie inférieure interagissant avec le fluide) est discrétisée, tandis que les propriétés dynamiques (masse, centre de gravité, moment d'inertie) sont calculées pour la sphère complète.
  • Couplage complet entre les forces hydrodynamiques (ISPH) et la dynamique du corps rigide (équations d'Euler pour la rotation, bilan de quantité de mouvement pour la translation).

• Modélisation de la Collision :

  • Détection de collision précise par interpolation linéaire du temps d'impact.
  • Le rebond sur le substrat solide est traité avec un coefficient de restitution sec ($e_r \approx 0.98$), supposé constant et élevé pour isoler les pertes d'énergie dues au liquide.
  • Adaptation dynamique du pas de temps ($\Delta t$) lors de l'approche du mur pour assurer la stabilité numérique.

• Validation :

  • Les simulations ont été validées contre les données expérimentales de Gollwitzer et al. (impact de billes de verre sur un film de silicone M5).
  • Une étude de convergence spatiale a été menée, montrant une erreur relative $L_2$ inférieure à 1% pour une résolution de 178 220 particules liquides.

3. Contributions Clés

1. Identification d'une loi d'échelle à deux paramètres : L'étude démontre que le COR ne peut pas être décrit uniquement par le nombre de Stokes ($St$). Une nouvelle variable adimensionnelle, l'épaisseur de film adimensionnelle ($\gamma = \delta/D$, où $\delta$ est l'épaisseur du film et $D$ le diamètre de la bille), est introduite comme paramètre de contrôle essentiel.
2. Découverte de deux régimes distincts : L'analyse révèle l'existence de deux régimes de collision avec des exposants de loi de puissance différents, séparés par une frontière dans l'espace des paramètres $(St, \gamma)$.
3. Algorithme de réduction de coût : Développement d'une méthode efficace pour simuler l'interaction fluide-solide en ne discrétisant qu'une fraction de la particule rigide, réduisant considérablement le coût computationnel sans perte de précision sur la cinématique de rebond.

4. Résultats Principaux

Les simulations couvrent une large gamme de vitesses d'impact, de diamètres de billes (2,8 à 8 mm) et d'épaisseurs de film (eau et huile silicone).

• Échec du modèle à un seul paramètre : À nombre de Stokes constant, le COR varie systématiquement avec le diamètre de la particule (et donc avec $\gamma$). Les courbes ne s'effondrent pas sur une seule tendance, invalidant l'hypothèse d'une dépendance exclusive à $St$.
• Lois d'échelle obtenues : Une régression linéaire dans l'espace logarithmique ($\ln(COR)$ en fonction de $\ln(St)$ et $\ln(\gamma)$) a permis d'identifier deux régimes :
  • Régime R1 (Faible $\gamma$, grosses particules) : Le COR dépend fortement à la fois de l'inertie et de la géométrie.
    $$COR \approx 0.17 \cdot St^{0.17} \cdot \gamma^{-0.16}$$
    Ici, les effets inertiels et visqueux sont tous deux significatifs.
  • Régime R2 (Fort $\gamma$, petites particules) : Le COR dépend principalement de la géométrie du film, avec une dépendance très faible au nombre de Stokes.
    $$COR \approx 0.4 \cdot St^{0.01} \cdot \gamma^{-0.22}$$
    La dépendance à $St$ est négligeable ($p \approx 0.053$).
• Analyse énergétique :
  • Dans le Régime R1, l'énergie cinétique du fluide diminue de manière monotone après le contact avec le substrat, indiquant une dissipation dominée par la viscosité et l'écoulement primaire.
  • Dans le Régime R2, une seconde augmentation de l'énergie cinétique (un pic secondaire) est observée, suggérant la présence de phénomènes d'écoulement secondaires (tourbillons, ondes) qui jouent un rôle majeur dans la dissipation d'énergie, expliquant la faible sensibilité au nombre de Stokes.

5. Signification et Conclusion

Cette étude apporte une avancée significative dans la compréhension fondamentale des collisions humides :

• Correction des modèles existants : Elle réfute l'idée que le nombre de Stokes suffit à prédire le comportement de rebond, en montrant que le rapport épaisseur/diamètre ($\gamma$) est un paramètre critique, surtout pour les petites particules ou les films épais.
• Application pratique : La loi d'échelle proposée (avec ses deux régimes) fournit des modèles robustes pour les simulations numériques de dynamique des particules discrètes (DEM) et d'autres approches statistiques, permettant de mieux prédire les pertes d'énergie dans les écoulements granulaires humides industriels.
• Mécanisme physique : L'identification de la transition entre un régime de dissipation visqueuse pure et un régime dominé par les effets inertiels secondaires (tourbillons) ouvre la voie à une modélisation plus fine des interactions fluide-structure à l'échelle microscopique.

En résumé, les auteurs établissent que la dynamique de rebond d'une particule sur un film liquide est gouvernée par une interaction complexe entre l'inertie de la particule et la géométrie du confinement liquide, nécessitant une description à deux paramètres pour une prédiction précise.