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

Dieser Beitrag stellt eine patientenspezifische Gitter-Boltzmann-Methode vor, die nicht-newtonsche Blutströmung in gekapselten intrakraniellen Aneurysmen modelliert, indem die Spulmasse als inhomogenes poröses Medium behandelt wird, wodurch ein Arbeitsablauf zur Bewertung der posttherapeutischen Hämodynamik validiert wird.

Das Wesentliche

Das große Ganze: Ein undichter Luftballon reparieren

Stellen Sie sich einen winzigen, schwachen Punkt auf einem Wasserballon in Ihrem Kopf vor. Medizinisch ausgedrückt handelt es sich dabei um ein intrakranielles Aneurysma. Wenn es platzt, verursacht es eine gefährliche Art von Schlaganfall.

Um dies zu beheben, wenden Ärzte häufig eine Technik namens Coiling an. Sie führen einen winzigen Draht in den Ballon ein und füllen ihn mit einem Geflecht aus Metalldrähten auf. Denken Sie dabei daran, wie man einen Schwamm in einen undichten Eimer stopft. Das Ziel ist es, das rasende Wasser (Blut) im Inneren des Ballons zu verlangsamen, damit es nicht mehr gegen die schwachen Wände drückt, und dem Ballon so eine Chance zu geben, zu heilen oder sich zu vernarben.

Das Problem: Das Ergebnis erraten

Der knifflige Teil besteht darin, dass der „undichte Ballon" jedes Patienten eine andere Form hat und jeder Arzt den „Schwamm" (die Coils) etwas anders einfüllt. Vor der Operation ist es sehr schwierig, genau vorherzusagen, wie sich der Wasserfluss ändern wird, sobald die Coils im Inneren sind. Ärzte verlassen sich derzeit stark auf ihre früheren Erfahrungen, um zu erraten, ob die Behandlung funktionieren wird.

Die Lösung: Ein digitaler „Flugsimulator"

Die Autoren dieses Papiers haben ein Computerprogramm entwickelt, das wie ein Flugsimulator für den Blutfluss funktioniert. Anstatt zu raten, können sie einen virtuellen Test durchführen, um genau zu sehen, wie sich das Blut verhalten wird, nachdem die Coils eingesetzt wurden.

So haben sie diesen Simulator aufgebaut:

1. Der Trick mit dem „porösen Schwamm"
Normalerweise muss ein Computer, um Wasser simulieren zu wollen, das durch ein chaotisches Haufen von Metalldrähten (den Coils) fließt, jeden einzelnen Draht zeichnen. Das ist wie der Versuch, jedes einzelne Sandkorn in einem Eimer zu zählen; es dauert ewig und erfordert einen Supercomputer.

Die Autoren fanden einen schlaueren Weg. Anstatt jeden Draht zu zeichnen, behandelten sie den Haufen Coils als einen Schwamm.

• Die Analogie: Stellen Sie sich vor, Sie gießen Wasser durch einen dichten Wald. Sie könnten versuchen zu berechnen, wie das Wasser um jeden einzelnen Baumstamm fließt (sehr schwierig). Oder Sie könnten einfach sagen: „Dieses ganze Gebiet ist ein dicker, langsam fließender Sumpf" (viel einfacher).
• Die Wissenschaft: Sie verwendeten ein mathematisches Modell namens Volume-Averaged Navier-Stokes (VANSE). Dies behandelt den mit Coils gefüllten Bereich als ein „poröses Medium" (einen Schwamm), bei dem das Wasser sich verlangsamt, je dichter der Schwamm ist.

2. Der Faktor „intelligentes Blut"
Blut ist nicht wie Wasser; es ist „zäh" und ändert sein Fließverhalten je nach Geschwindigkeit (wie Ketchup oder Honig). Das Computermodell berücksichtigt dieses „nicht-newtonsche" Verhalten, sodass die Simulation sich wie echtes Blut anfühlt und nicht nur wie einfaches Wasser.

3. Der „schnelle Motor" (Lattice-Boltzmann-Methode)
Um diese Berechnungen schnell genug laufen zu lassen, damit sie nützlich sind, verwendeten sie einen speziellen mathematischen Motor namens Lattice-Boltzmann-Methode (LBM).

• Die Analogie: Denken Sie daran wie an eine Hochgeschwindigkeits-Game-Engine. Während andere Methoden versuchen könnten, die Physik des gesamten Ozeans auf einmal zu lösen, zerlegt LBM den Ozean in winzige, handhabbare Kacheln und simuliert, wie Partikel zwischen ihnen hin und her springen. Dies ermöglicht es, die Simulation auf modernen Grafikkarten (GPUs) unglaublich schnell laufen zu lassen.

Was sie in der Studie taten

Das Team nahm einen echten CT-Scan des Gehirns eines Patienten mit einem Aneurysma und erstellte ein 3D-Modell davon. Anschließend führten sie zwei Arten von Simulationen durch:

1. Die „super detaillierte" Version: Sie modellierten jeden einzelnen Draht der Coils (wie das Zählen jedes Baumes im Wald).
2. Die „Schwamm"-Version: Sie verwendeten ihr neues „poröses Medium"-Modell (wie die Behandlung des Waldes als Sumpf).

Sie testeten dies mit drei verschiedenen Stufen der „Verstopfung" (Packungsdichten von 15 %, 20 % und 25 %), was den Mengen entspricht, die Ärzte tatsächlich in der Operation verwenden.

Die Ergebnisse: Der Schwamm funktioniert!

Die Ergebnisse waren ermutigend:

• Genauigkeit: Das „Schwamm"-Modell lieferte fast die gleichen Ergebnisse wie das „super detaillierte" Modell. Die „Schwamm"-Methode war viel schneller, verlor aber nicht den Überblick über das große Ganze.
• Flussreduktion: Als sie mehr Coils hinzufügten (den Schwamm dichter machten), verlangsamte sich der Blutfluss im Aneurysma erheblich.
• Sicherheit: Die „Wandschubspannung" (die Kraft, mit der das Blut gegen die schwache Wand reibt) sank drastisch. Beim unbehandelten Aneurysma wurde die Wand stark bombardiert. Mit den Coils sank die Kraft um etwa 40 %, was auf ein deutlich geringeres Risiko für das Platzen des Ballons hindeutet.

Das Fazit

Dieses Papier stellt eine neue, schnellere Methode zur Simulation von Behandlungen bei Gehirnaneurysmen vor. Indem sie die Metalldrähte als einen „Schwamm" behandeln, anstatt jeden Draht zu zählen, könnten Ärzte patientenspezifische Simulationen schnell durchführen. Dieser Arbeitsablauf ermöglicht eine bessere Einschätzung, ob ein bestimmter Behandlungsplan den Blutfluss erfolgreich verlangsamen und den Patienten schützen wird, und bewegt sich weg von reinem Raten hin zu datengestützten Entscheidungen.

Technische Zusammenfassung

Technischer Überblick: Eine Gitter-Boltzmann-Methode für nicht-newtonsche Blutströmung in gewickelten intrakraniellen Aneurysmen

Problemstellung
Intrakranielle Aneurysmen sind eine der Hauptursachen für hämorrhagische Schlaganfälle und werden häufig durch eine Coil-Embolisation behandelt, um eine Strömungsreduktion und eine Verringerung der Wandschubspannung (WSS) herbeizuführen. Die Vorhersage der hämodynamischen Ergebnisse solcher Behandlungen bleibt jedoch eine offene Herausforderung. Eine genaue Vorhersage erfordert patientenspezifische Geometrien und Parameter, insbesondere die Darstellung der inhomogenen Coil-Verteilung innerhalb des Aneurysmas. Traditionelle Ansätze sehen sich einem Zielkonflikt gegenüber: Vollständig aufgelöste Simulationen diskreter Coil-Geometrien sind rechenintensiv, während vereinfachte Modelle oft die Heterogenität der Coil-Packung nicht erfassen, was für korrekte Strömungsvorhersagen entscheidend ist. Darüber hinaus ist Blut eine nicht-newtonsche Flüssigkeit, was die Simulation rheologisch komplexer macht.

Methodik
Die Autoren schlagen einen Workflow vor, der patientenspezifische Geometrie mit einem volumenmittelnden Ansatz kombiniert, um das gewickelte Aneurysma als inhomogenes poröses Medium zu modellieren. Die Kernmethodik umfasst:

1. Governing Equations (Grundgleichungen): Das Modell verwendet die volumenmittelten Navier-Stokes-Gleichungen (VANSE) für nicht-newtonsche Strömung. Das Fluidgebiet wird durch ein Porositätsfeld ($\phi$) charakterisiert, das den Volumenanteil der Flüssigkeit darstellt. Die Impulsgleichung enthält einen Volumenkraftterm ($f(u)$), der den Darcy- und Forchheimer-Widerstand modelliert, um das poröse Medium (die Coils) zu berücksichtigen.
2. Rheologie: Die Blutviskosität wird mittels des Carreau-Yasuda-Modells modelliert, das das scherverdünnende Verhalten berücksichtigt. Die dynamische Viskosität hängt von der Scherrate ab.
3. Numerisches Schema: Die Gleichungen werden mit einer problemangepassten Gitter-Boltzmann-Methode (LBM) auf einem D3Q27-Gitter gelöst.
   • Handhabung poröser Medien: Die LBM integriert räumlich variierende Porosität und eine variable Relaxationsrate, um die nicht-newtonsche Viskosität und die porösen Widerstandskräfte zu behandeln.
   • Kraftimplementierung: Das Guo-Kraftschema wird verwendet, um die Widerstandskräfte einzubinden, die quadratisch von der Geschwindigkeit abhängen.
   • Geschwindigkeitslösung: Für isotrope Permeabilität wird die Geschwindigkeit explizit gelöst. Für anisotrope Permeabilität wird eine Fixpunktiteration vorgeschlagen, wobei die Autoren anmerken, dass dies die Rechenzeit erheblich erhöht.
4. Geometrie und Simulationseinrichtung:
   • Geometrie: Eine patientenspezifische Aneurysmageometrie, abgeleitet aus einem anonymisierten CT-Scan, wurde verwendet.
   • Coil-Darstellung: Coils wurden mit einem diskreten elastischen Stabmodell mit Packungsdichten von 15 %, 20 % und 25 % simuliert. Die resultierenden diskreten Coil-Strukturen wurden mittels Fensterfaltung in kontinuierliche Porositätsfelder umgewandelt.
   • Randbedingungen: Haftbedingungen wurden an den Wänden angewendet (Half-Way-Bounce-Back). Der Zufluss wurde durch ein zeit aufgelöstes Geschwindigkeitsprofil angetrieben, das aus einem 1D-0D-arteriellen Netzwerkmodell abgeleitet wurde, während der Abfluss eine Extrapolationsmethode verwendete.
   • Implementierung: Die Simulationen wurden mit dem XLB-Framework implementiert und nutzten GPU-Beschleunigung (Nvidia RTX A6000) für die parallele Berechnung.

Hauptbeiträge

• Hybrider Modellierungsansatz: Die Arbeit integriert LBM für VANSE mit nicht-newtonscher Rheologie, um die Blutströmung in gewickelten Aneurysmen zu simulieren, wobei die Coil-Masse als poröses Medium behandelt wird, anstatt einzelne Drähte aufzulösen.
• Validierung der Mittelung: Die Studie bietet einen direkten Vergleich zwischen vollständig aufgelösten Coil-Simulationen und dem volumenmittelnden Ansatz für poröse Medien.
• Patientenspezifischer Workflow: Die Autoren demonstrieren eine vollständige Pipeline von der medizinischen Bildgebung (CT-Scan) bis zur hämodynamischen Bewertung, einschließlich der Generierung von Porositätsfeldern aus simulierten Coil-Implantationen.

Ergebnisse

• Strömungsreduktion: Die Simulationen bestätigen, dass die Coil-Einführung die Strömungsgeschwindigkeit innerhalb des Aneurysmasacks signifikant reduziert, während die Strömung im Muttergefäß erhalten bleibt. Die Reduktion ist bei höheren Packungsdichten (bis zu 25 %) ausgeprägter.
• Modellgültigkeit: Die volumenmittelten (porösen Medium) Simulationen zeigen eine gute Übereinstimmung mit vollständig aufgelösten Simulationen. Während der Mittelungsansatz lokale Feinskaleneffekte verliert, sind das gesamte Strömungsverhalten und makroskopische Größen konsistent.
• Wandschubspannung (WSS): Die Behandlung reduziert die WSS-Stärke drastisch. Im unbehandelten Aneurysma erreichte die WSS etwa 2 Pa, während gewickelte Aneurysma Werte unter 1,5 Pa aufwiesen, was auf ein verringertes Rupturrisiko hindeutet. Der Unterschied in der WSS zwischen den drei Packungsdichten (15 %, 20 %, 25 %) erwies sich als gering, was darauf hindeutet, dass die Verteilung der Coil-Platzierung möglicherweise einflussreicher ist als die Packungsdichte allein.
• Quantitative Metriken: Über den Aneurysmabereich hinweg sagte das volumenmittelnde Modell eine Reduktion der Strömungsstärke von etwa 50 % und eine WSS-Reduktion von etwa 40 % im Vergleich zum unbehandelten Fall voraus.

Bedeutung und Behauptungen
Die Autoren behaupten, dass dieser Workflow eine patientenspezifische Bewertung von Blutströmungsänderungen aufgrund potenzieller Coil-Behandlungen ermöglicht. Die primäre Bedeutung liegt in der Gültigkeit des volumenmittelnden Ansatzes für poröse Medien, der eine rechen-effiziente Alternative zu vollständig aufgelösten Simulationen bietet, ohne die Genauigkeit globaler Strömungsmetriken zu beeinträchtigen.

Die Arbeit schließt bescheiden, dass die Berücksichtigung anisotroper Permeabilität theoretisch die Genauigkeit erhöhen könnte, die aktuelle isotrope Annahme jedoch für die klinische Praxis wahrscheinlich ausreicht, angesichts der Unsicherheiten bei der genauen Coil-Platzierung und der Notwendigkeit schneller Entscheidungsfindung. Die Autoren identifizieren zukünftige Arbeiten als notwendig für eine detaillierte Parameterkalibrierung, Unsicherheitsquantifizierung und die Analyse von Erweiterungen für anisotrope Permeabilität.