A Lattice Boltzmann Method for Non-Newtonian Blood Flow in Coiled Intracranial Aneurysms

Cet article présente une méthode de Boltzmann sur réseau spécifique au patient qui modélise l'écoulement sanguin non newtonien dans des anévrismes intracrâniens embolisés en traitant la masse de coils comme un milieu poreux hétérogène, validant ainsi une démarche d'évaluation de l'hémodynamique post-traitement.

L'essentiel

La vue d'ensemble : Réparer un ballon percé

Imaginez un tout petit point faible sur un ballon d'eau à l'intérieur de votre tête. En termes médicaux, il s'agit d'un anévrisme intracrânien. S'il éclate, cela provoque un type d'AVC dangereux.

Pour réparer cela, les médecins utilisent souvent une technique appelée coiffage (ou embolisation par coils). Ils enfilent un fil fin dans le ballon et le remplissent d'une maille de coils métalliques. Pensez-y comme à l'action de bourrer une éponge dans un seau qui fuit. L'objectif est de ralentir l'eau qui s'écoule rapidement (le sang) à l'intérieur du ballon afin qu'elle cesse de frapper les parois fragiles, donnant ainsi au ballon une chance de cicatriser ou de se recouvrir de tissu cicatriciel.

Le problème : Deviner le résultat

La partie délicate est que le « ballon percé » de chaque patient a une forme différente, et chaque médecin bourre l'« éponge » (les coils) d'une manière légèrement différente. Avant la chirurgie, il est très difficile de prédire exactement comment l'écoulement de l'eau changera une fois les coils à l'intérieur. Les médecins s'appuient actuellement fortement sur leur expérience passée pour deviner si le traitement fonctionnera.

La solution : Un « simulateur de vol » numérique

Les auteurs de cet article ont créé un programme informatique qui agit comme un simulateur de vol pour la circulation sanguine. Au lieu de deviner, ils peuvent effectuer un test virtuel pour voir exactement comment le sang se comportera après l'insertion des coils.

Voici comment ils ont construit ce simulateur :

1. L'astuce de l'« éponge poreuse »
Habituellement, pour simuler l'eau s'écoulant à travers un amas désordonné de fils métalliques (les coils), un ordinateur doit dessiner chaque fil individuellement. C'est comme essayer de compter chaque grain de sable dans un seau ; cela prend une éternité et nécessite un super-ordinateur.

Les auteurs ont trouvé une méthode plus intelligente. Au lieu de dessiner chaque fil, ils ont traité l'amoncellement de coils comme une éponge.

• L'analogie : Imaginez que vous versez de l'eau à travers une forêt dense. Vous pourriez essayer de calculer comment l'eau s'écoule autour de chaque tronc d'arbre (très difficile). Ou, vous pourriez simplement dire : « Toute cette zone est un marais épais et lent » (beaucoup plus facile).
• La science : Ils ont utilisé un modèle mathématique appelé Navier-Stokes moyennés en volume (VANSE). Cela traite la zone remplie de coils comme un « milieu poreux » (une éponge) où l'eau ralentit en fonction de la densité de l'éponge.

2. Le facteur « sang intelligent »
Le sang n'est pas comme l'eau ; il est « épais » et change de comportement d'écoulement selon sa vitesse (comme le ketchup ou le miel). Le modèle informatique prend en compte ce comportement « non newtonien », garantissant que la simulation ressemble à du sang réel et non à de l'eau ordinaire.

3. Le « moteur rapide » (Méthode Lattice Boltzmann)
Pour que ces calculs s'exécutent assez rapidement pour être utiles, ils ont utilisé un moteur mathématique spécifique appelé la Méthode Lattice Boltzmann (LBM).

• L'analogie : Pensez-y comme à un moteur de jeu vidéo haute vitesse. Tandis que d'autres méthodes pourraient essayer de résoudre la physique de tout l'océan en une seule fois, la LBM décompose l'océan en de minuscules tuiles gérables et simule comment les particules rebondissent entre elles. Cela permet à la simulation de s'exécuter incroyablement vite sur les cartes graphiques modernes (GPU).

Ce qu'ils ont fait dans l'étude

L'équipe a pris une véritable IRM d'un anévrisme cérébral de patient et a construit un modèle 3D de celui-ci. Ils ont ensuite exécuté deux types de simulations :

1. La version « Super détaillée » : Ils ont modélisé chaque fil individuel du coil (comme compter chaque arbre dans la forêt).
2. La version « Éponge » : Ils ont utilisé leur nouveau modèle de « milieu poreux » (comme traiter la forêt comme un marais).

Ils ont testé cela avec trois niveaux différents de « densité de bourrage » (densités de remplissage de 15 %, 20 % et 25 %), qui sont les quantités réellement utilisées par les médecins en chirurgie.

Les résultats : L'éponge fonctionne !

Les résultats étaient encourageants :

• Précision : Le modèle « Éponge » a donné des résultats presque identiques au modèle « Super détaillé ». La méthode « Éponge » était beaucoup plus rapide sans perdre la vue d'ensemble.
• Réduction du flux : À mesure qu'ils ajoutaient plus de coils (rendant l'éponge plus dense), l'écoulement sanguin à l'intérieur de l'anévrisme ralentissait considérablement.
• Sécurité : La « contrainte de cisaillement pariétal » (la force du sang frottant contre la paroi fragile) chutait drastiquement. Dans l'anévrisme non traité, la paroi était violemment battue. Avec les coils, la force chutait d'environ 40 %, suggérant un risque beaucoup plus faible d'éclatement du ballon.

La conclusion

Cet article présente une nouvelle méthode plus rapide pour simuler les traitements d'anévrismes cérébraux. En traitant les coils métalliques comme une « éponge » plutôt que de compter chaque fil, les médecins pourraient potentiellement exécuter rapidement des simulations spécifiques au patient. Ce flux de travail permet une meilleure évaluation de savoir si un plan de traitement spécifique ralentira avec succès l'écoulement sanguin et protégera le patient, passant d'une pure devinette à des décisions fondées sur les données.

Résumé technique

Résumé Technique : Une Méthode de Boltzmann sur Réseau pour l'Écoulement Sanguin Non Newtonien dans des Aneurysmes Intracrâniens Enroulés

Énoncé du Problème
Les anévrismes intracrâniens sont une cause majeure d'accident vasculaire cérébral hémorragique, souvent traités par embolisation par coils pour induire une réduction de l'écoulement et atténuer la contrainte de cisaillement pariétal (WSS). Cependant, la prédiction des résultats hémodynamiques de tels traitements reste un défi ouvert. Une prédiction précise nécessite une géométrie et des paramètres spécifiques au patient, en particulier la représentation de la distribution hétérogène des coils à l'intérieur de l'anévrisme. Les approches traditionnelles font face à un compromis : les simulations entièrement résolues de géométries de coils discrètes sont coûteuses en calcul, tandis que les modèles simplifiés échouent souvent à capturer l'hétérogénéité de l'empilement des coils, ce qui est crucial pour des prédictions d'écoulement correctes. De plus, le sang est un fluide non newtonien, ajoutant une complexité rhéologique à la simulation.

Méthodologie
Les auteurs proposent un flux de travail combinant une géométrie spécifique au patient avec une approche moyennée en volume pour modéliser l'anévrisme enroulé comme un milieu poreux hétérogène. La méthodologie centrale comprend :

1. Équations Gouvernantes : Le modèle utilise les Équations de Navier-Stokes Moyennées en Volume (VANSE) pour un écoulement non newtonien. Le domaine fluide est caractérisé par un champ de porosité ($\phi$) représentant la fraction volumique de fluide. L'équation de quantité de mouvement inclut un terme de force volumique ($f(u)$) modélisant la traînée de Darcy et de Forchheimer pour tenir compte du milieu poreux (les coils).
2. Rhéologie : La viscosité du sang est modélisée à l'aide du modèle de Carreau–Yasuda, qui rend compte du comportement de fluidité par cisaillement. La viscosité dynamique dépend du taux de déformation de cisaillement.
3. Schéma Numérique : Les équations sont résolues à l'aide d'une Méthode de Boltzmann sur Réseau (LBM) adaptée au problème sur un réseau D3Q27.
   • Gestion des Milieux Poreux : La LBM intègre une porosité spatialement variable et un taux de relaxation variable pour gérer la viscosité non newtonienne et les forces de traînée poreuse.
   • Implémentation de la Force : Le schéma de forçage de Guo est utilisé pour intégrer les forces de traînée, qui dépendent quadratiquement de la vitesse.
   • Solution de Vitesse : Pour une perméabilité isotrope, la vitesse est résolue explicitement. Pour une perméabilité anisotrope, une itération de point fixe est proposée, bien que les auteurs notent que cela augmente considérablement le temps de calcul.
4. Géométrie et Configuration de Simulation :
   • Géométrie : Une géométrie d'anévrisme spécifique au patient, dérivée d'un scanner CT anonymisé, a été utilisée.
   • Représentation des Coils : Les coils ont été simulés à l'aide d'un modèle de tige élastique discrète avec des densités de remplissage de 15 %, 20 % et 25 %. Les structures de coils discrètes résultantes ont été converties en champs de porosité continus par convolution fenêtrée.
   • Conditions aux Limites : Des conditions d'adhérence (no-slip) ont été appliquées aux parois (rebond à mi-chemin). L'entrée était pilotée par un profil de vitesse résolu dans le temps dérivé d'un modèle de réseau artériel 1D-0D, tandis que la sortie utilisait une méthode d'extrapolation.
   • Implémentation : Les simulations ont été implémentées en utilisant le cadre XLB, exploitant l'accélération GPU (Nvidia RTX A6000) pour le calcul parallèle.

Contributions Clés

• Approche de Modélisation Hybride : L'article intègre la LBM pour les VANSE avec une rhéologie non newtonienne pour simuler l'écoulement sanguin dans des anévrismes enroulés, traitant la masse de coils comme un milieu poreux plutôt que de résoudre les fils individuels.
• Validation du Moyennage : L'étude fournit une comparaison directe entre les simulations de coils entièrement résolues et l'approche de milieu poreux moyennée en volume.
• Flux de Travail Spécifique au Patient : Les auteurs démontrent un pipeline complet, allant de l'imagerie médicale (scanner CT) à l'évaluation hémodynamique, incluant la génération de champs de porosité à partir de déploiements de coils simulés.

Résultats

• Réduction de l'Écoulement : Les simulations confirment que l'insertion de coils réduit considérablement la vitesse d'écoulement à l'intérieur du sac anévrismal tout en maintenant l'écoulement dans le vaisseau parent. La réduction est plus prononcée avec des densités de remplissage plus élevées (jusqu'à 25 %).
• Validité du Modèle : Les simulations de milieu poreux (moyennées en volume) montrent un bon accord avec les simulations entièrement résolues. Bien que l'approche de moyennage perde les effets locaux à petite échelle, le comportement global de l'écoulement et les grandeurs macroscopiques sont cohérents.
• Contrainte de Cisaillement Pariétal (WSS) : Le traitement réduit drastiquement l'amplitude du WSS. Dans l'anévrisme non traité, le WSS atteignait environ 2 Pa, tandis que les anévrismes enroulés présentaient des amplitudes inférieures à 1,5 Pa, indiquant un risque réduit de rupture. La différence de WSS entre les trois densités de remplissage (15 %, 20 %, 25 %) s'est avérée légère, suggérant que la distribution du placement des coils pourrait être plus influente que la densité de remplissage seule.
• Métriques Quantitatives : Sur toute la région de l'anévrisme, le modèle moyenné en volume a prédit une réduction de l'amplitude de l'écoulement d'environ 50 % et une réduction du WSS d'environ 40 % par rapport au cas non traité.

Signification et Revendications
Les auteurs affirment que ce flux de travail permet une évaluation spécifique au patient des changements d'écoulement sanguin dus à des traitements potentiels par coils. La signification principale réside dans la validité de l'approche de milieu poreux moyennée en volume, qui offre une alternative efficace en calcul aux simulations entièrement résolues sans sacrifier la précision des métriques d'écoulement globales.

L'article conclut modestement que, bien que la prise en compte de la perméabilité anisotrope pourrait théoriquement augmenter la précision, l'hypothèse isotrope actuelle est probablement suffisante pour la pratique clinique étant donné les incertitudes sur le placement exact des coils et la nécessité d'une prise de décision rapide. Les auteurs identifient un travail futur nécessaire pour l'étalonnage détaillé des paramètres, la quantification des incertitudes et l'analyse des extensions de perméabilité anisotrope.