🔬 condensed matter

Investigating frictional instability due to pressurization in granular media: insights from coupled computational fluid dynamics discrete element method

Cette étude emploie des simulations 3D couplées CFD-DEM pour démontrer que l'instabilité induite par le fluide dans les milieux granulaires n'est pas régie par la seule pression interstitielle, mais par des rétroactions hydromécaniques complexes impliquant les conditions de drainage, l'évolution de la structure granulaire et la transition d'un renforcement dominé par la dilatance vers un affaiblissement piloté par la pression interstitielle.

Auteurs originaux : Bimal Chhushyabaga (Department of Civil and Environmental Engineering, University of Houston), Behrooz Ferdowsi (Department of Civil and Environmental Engineering, University of Houston)

Publié 2026-06-05

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Bimal Chhushyabaga (Department of Civil and Environmental Engineering, University of Houston), Behrooz Ferdowsi (Department of Civil and Environmental Engineering, University of Houston)

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez un tas de sable, comme un bac à sable, mais celui-ci est imbibé d'eau et serré étroitement entre deux plaques lourdes. C'est la configuration d'une étude sur la façon dont les matériaux granulaires humides (comme le sol, le sable ou la poussière de roche) échouent soudainement ou glissent lorsque la pression de l'eau augmente à l'intérieur d'eux.

Les chercheurs ont utilisé une simulation informatique surpuissante pour observer ce phénomène particule par particule. Imaginez un film en haute définition où chaque grain de sable et chaque goutte d'eau est suivi individuellement.

Voici ce qu'ils ont découvert, expliqué par des analogies simples :

1. La configuration : Le « Sandwich » et la « Cocotte-minute »

Imaginez un sandwich où le pain est une plaque lourde et la garniture du sable mouillé.

La contrainte : La plaque supérieure appuie vers le bas (gravité) et tente de glisser latéralement (contrainte de cisaillement).
Le déclencheur : Ils ont commencé à pomper la pression de l'eau depuis le bas du sandwich, comme si l'on gonflait un ballon à l'intérieur de la garniture.
L'objectif : Ils voulaient voir à quel moment la pression de l'eau ferait perdre leur adhérence au sable, provoquant un glissement incontrôlable.

2. Les deux scénarios : La « Boîte scellée » vs La « Fenêtre ouverte »

Les chercheurs ont testé deux manières différentes dont l'eau pouvait se comporter, ce qui a radicalement changé le résultat :

La « Boîte scellée » (Non drainée) : Imaginez que le haut du sandwich est scellé hermétiquement. Lorsque l'on pompe l'eau à l'intérieur, la pression n'a nulle part où aller. Elle s'accumule rapidement, écartant les grains et les faisant flotter légèrement.
- Le résultat : Cela a provoqué un effondrement soudain et violent. Le sable a perdu sa force instantanément, et la plaque supérieure a glissé rapidement. C'était comme un tremblement de terre soudain.
La « Fenêtre ouverte » (Drainée) : Imaginez que le haut du sandwich possède un petit trou ou une éponge qui laisse l'eau s'échapper. Lorsque l'on pompe l'eau à l'intérieur, une partie s'échappe.
- Le résultat : La pression ne s'est pas accumulée aussi haut. Le sable ne s'est pas effondré d'un coup. Au lieu de cela, il s'est affaibli lentement et régulièrement, comme une descente lente sur une colline plutôt qu'une chute brutale.

3. L'ingrédient secret : Ce n'est pas seulement la pression de l'eau

Une idée reçue courante est que la pression de l'eau seule cause le glissement. L'article dit non.

Imaginez les grains de sable comme une équipe de personnes se tenant la main en cercle (une « chaîne de force »).

Le travail d'équipe : Tant qu'elles se tiennent la main fermement, l'équipe est forte.
L'eau : La pression de l'eau essaie d'écarter les personnes.
La véritable cause : Le glissement ne se produit que lorsque l'eau pousse assez fort ET que les personnes (les grains) ont déjà lâché la main les unes des autres.
- Si l'équipe est déjà faible (le sable était déjà proche de glisser), un peu de pression d'eau provoque un crash.
- Si l'équipe est forte (le sable est loin de glisser), il faut beaucoup de pression d'eau pour briser l'étreinte.

L'étude a découvert que ce « lâcher prise » se produit parce que les grains commencent à pivoter et à se réorganiser, brisant les chaînes solides qui maintiennent la structure.

4. Le « Embouteillage » de grains

À l'intérieur du sable, les choses ne bougent pas uniformément.

Dilatance et compaction : Certaines parties du sable ont tenté de s'étendre (dilatance), comme une foule de personnes essayant de s'éparpiller, tandis que d'autres parties ont été compressées (compaction).
Le changement de perméabilité : Ce réarrangement a modifié la facilité avec laquelle l'eau pouvait circuler à travers le sable. C'est comme si un embouteillage se débloquait soudainement dans une voie pendant qu'une autre voie se bouche. L'eau a trouvé de nouveaux chemins pour circuler, ce qui a modifié l'endroit où la pression était la plus élevée.

5. Le « Rampement » vs Le « Crash »

Les chercheurs ont observé la vitesse à laquelle le sable se déplaçait.

Rampement lent (Creep) : Au début, le sable bougeait très lentement, presque comme s'il hésitait. Pendant cette période, le sable est devenu temporairement plus fort car les grains se bousculaient pour trouver une position plus serrée et plus stable (comme des gens trouvant leur équilibre).
Le Crash : Finalement, la pression de l'eau a submergé cette stabilité. Les grains ont cessé de trouver un appui et ont commencé à tourner follement. Les « chaînes de force » (le fait de se tenir la main) se sont effondrées, et la friction a chuté, entraînant un glissement rapide.

6. Le « Nombre visqueux » (La limite de vitesse)

Les chercheurs ont essayé de comparer leur sable à une soupe épaisse ou une suspension dense (comme un smoothie). Ils ont utilisé un concept mathématique appelé le « nombre visqueux » pour voir si le sable se comportait comme un fluide.

La découverte : Cela ne correspondait pas tout à fait aux règles d'une « soupe » standard. Le sable était trop encombré et les grains étaient trop solides. Il se comportait davantage comme un embouteillage lent où les voitures (les grains) se cognent occasionnellement et se réorganisent, plutôt que comme le flux fluide d'un liquide. Les règles standards pour les fluides épais ne prédisaient pas parfaitement ce qui s'est passé ici.

Résumé

L'article conclut que lorsque le sable humide échoue, ce n'est pas seulement la pression de l'eau qui pousse. C'est une danse complexe :

La pression de l'eau essaie d'écarter les grains.
Les grains tentent de se réorganiser pour maintenir leur prise.
Si l'eau ne peut pas s'échapper (boîte scellée), la pression gagne, et les grains perdent leur adhérence d'un seul coup.
Si l'eau peut s'échapper (fenêtre ouverte), la pression est relâchée, et la défaillance se produit lentement.
Crucialement, la défaillance ne se produit que lorsque le réseau de contact (le fait de se tenir la main) s'est suffisamment dégradé pour que les grains ne puissent plus supporter le poids qui presse vers le bas.

Résumé technique : Étude de l'instabilité frictionnelle due à la pressurisation dans les milieux granulaires

Énoncé du problème
La pressurisation fluide est un déclencheur primaire de la déformation et de l'instabilité dans les géomatériaux saturés, notamment les gouges de failles, les pentes affectées par les précipitations et les réservoirs perturbés par l'injection. Bien que les critères de rupture macroscopiques (ex. : Mohr–Coulomb) décrivent la réduction de la contrainte normale effective due à la pression interstitielle, ils ne parviennent pas à expliquer les rétroactions à l'échelle du grain liant la diffusion de la pression, le drainage et la dégradation du réseau de contacts à la transition d'un fluage stable vers un glissement dynamique. Plus précisément, les mécanismes régissant la manière dont les couches granulaires sous contrainte subcritique se réactivent sous une contrainte de cisaillement prescrite, et comment les conditions de drainage régulent ce processus, restent difficiles à résoudre à l'échelle particulaire.

Méthodologie
L'étude utilise des simulations numériques couplées en deux étapes, CFD–DEM (méthode des éléments discrets couplée à la mécanique des fluides computationnelle), en trois dimensions, pour étudier la réactivation induite par la pression interstitielle dans une couche de cisaillement granulaire confinée et saturée de fluide.

Cadre numérique : La phase fluide est résolue sur une grille eulérienne fixe à l'aide d'équations de bilan de masse et de quantité de mouvement moyennées par volume, tandis que les grains individuels sont suivis comme des particules lagrangiennes obéissant aux lois de Newton. L'échange de quantité de mouvement entre les phases est modélisé via des forces de gradient de pression et un traînée entravée de type Syamlal–O'Brien. Les interactions entre particules utilisent une loi de contact de Hertz à sphères souples avec frottement de Coulomb.
Protocole de chargement : L'étude utilise une approche en deux étapes. Premièrement, des simulations sous contrôle de déformation sous contrainte normale asservie établissent la résistance frictionnelle de référence et l'enveloppe de rupture de Mohr–Coulomb. Deuxièmement, des simulations sous contrôle de contrainte imposent des contraintes de cisaillement subcritiques ( $\tau = 0,46\tau_{ss}$ à $0,90\tau_{ss}$ ) tout en augmentant progressivement la pression interstitielle basale.
Conditions aux limites : Les simulations comparent des conditions non drainées (limite supérieure imperméable, retenant l'excès de pression interstitielle) avec des conditions drainées (pression atmosphérique au sommet, permettant la dissipation de la pression verticale).
Diagnostics : L'analyse intègre des métriques macroscopiques (friction apparente, vitesse de glissement, pression interstitielle excédentaire) avec des diagnostics micromécaniques (rotation des particules, réorganisation des chaînes de force, nombre de coordination, redistribution de la porosité) et des métriques d'anisotropie de tissu de second ordre.

Résultats clés

Rôle du drainage : L'instabilité n'est pas gouvernée par la pression interstitielle seule, mais par son évolution couplée avec la contrainte effective et le drainage.
- Non drainé : L'excès de pression interstitielle est retenu, entraînant une réduction rapide de la contrainte effective, un affaiblissement frictionnel abrupt et une accélération du glissement dynamique.
- Drainé : Des gradients de pression verticaux sont maintenus, supprimant l'accumulation de pression excédentaire. L'affaiblissement est plus fluide et progressif, car l'échappement du fluide tamponne la conversion du forçage de pression interstitielle en perte de contrainte effective.
Proximité de la contrainte et dégradation du tissu : L'incrément de pression requis pour déclencher la rupture dépend de l'état de contrainte initial par rapport à la résistance au cisaillement.
- Les couches proches de la rupture ( $\tau \approx 0,9\tau_{ss}$ ) échouent après de faibles augmentations de pression avec une dégradation moins étendue du tissu.
- Les couches éloignées de la rupture ( $\tau \approx 0,46\tau_{ss}$ ) nécessitent des accumulations de pression plus importantes et subissent une dégradation plus significative du réseau de contacts, incluant un effondrement plus marqué des chaînes de force et une réduction du nombre de coordination.
Évolution micromécanique : La rupture est caractérisée par la transition d'une déformation distribuée vers un glissement localisé.
- Chaînes de force : L'affaiblissement macroscopique coïncide avec l'effondrement et la réorganisation de l'ossature des forces fortes. Après la rupture, le cas à faible cisaillement présente une perte substantielle de l'organisation directionnelle des chaînes de force, tandis que le cas à cisaillement élevé conserve une structure anisotrope plus cohérente.
- Porosité et perméabilité : La rupture implique des bandes alternées de dilatation et de compaction. Ces changements volumétriques réorganisent la perméabilité effective dérivée de la porosité, indiquant que les voies hydrauliques évoluent dynamiquement pendant la déformation.
Mise à l'échelle rhéologique : La réponse frictionnelle s'organise partiellement avec le nombre visqueux ( $I_v$ ) dans le régime de faible $I_v$ ( $10^{-8} \lesssim I_v \lesssim 10^{-5}$ ). Cependant, les données ne se regroupent pas sur une rhéologie locale unique, ce qui est cohérent avec un régime de fluage dominé par des réarrangements intermittents et une transmission de contrainte non locale plutôt que par une simple loi $\mu(I_v)$ locale.

Signification et affirmations
L'article affirme fournir des informations à l'échelle du grain sur la manière dont les rétroactions hydromécaniques contrôlées par le drainage et la dégradation du tissu convertissent le forçage de la pression interstitielle en instabilité frictionnelle. Les principales contributions incluent :

La démonstration que la transition vers l'instabilité est régulée par la compétition entre la génération de pression, la diffusion hydraulique et le changement de volume des pores induit par la déformation.
La démonstration que les conditions de drainage modifient fondamentalement le chemin de rupture : les systèmes non drainés favorisent un affaiblissement abrupt via la rétention de pression, tandis que les systèmes drainés échouent progressivement via des gradients de pression.
L'établissement que le degré de dégradation du réseau de contacts requis pour la rupture dépend de l'état de contrainte initial ; la proximité de la rupture dicte l'étendue de la réorganisation du tissu nécessaire avant que le squelette granulaire ne puisse plus supporter la charge de cisaillement imposée.
La clarification que l'affaiblissement macroscopique est l'expression à l'échelle du système de la dégradation du réseau de contacts, spécifiquement la perte de la structure de support de charge directionnelle, plutôt qu'une simple fonction de la pression interstitielle globale ou de la perte de contact.

Les auteurs concluent que ces résultats offrent une base mécaniste pour interpréter l'instabilité induite par la pressurisation dans les gouges de failles, les pentes déclenchées par les précipitations et les systèmes perturbés par l'injection, soulignant que les couches mal drainées sont sujettes à un affaiblissement abrupt lorsque la dégradation locale du tissu se produit dans un contexte plus large de rétention de pression.