An open-source computational framework for immersed fluid-structure interaction modeling using FEBio and MFEM

Deze paper introduceert een nieuw open-source computerramwerk voor immersieve fluid-structure interaction-simulaties dat de schaalbaarheid van MFEM en de geavanceerde biomechanische modellering van FEBio combineert om complexe biologische systemen, zoals hartkleppen, nauwkeurig te modelleren.

De kern

Stel je voor dat je een heel complexe dans wilt filmen tussen twee partners: een stromende rivier (het bloed) en een flexibele, dansende jurk (een hartklep). De uitdaging is dat de jurk enorm beweegt, soms zelfs in contact komt met zichzelf, en de rivier moet meebewegen met elke beweging.

In de wereld van computersimulaties is dit een enorme hoofdpijn. Traditionele methoden proberen de rivier en de jurk in één groot, flexibel net te vangen. Maar als de jurk te veel beweegt of in contact komt met zichzelf, wordt dat net zo verdraaid dat de computer het niet meer kan volgen. Het is alsof je probeert een elastiekje te rekken tot het kapot gaat; je moet dan het hele net opnieuw maken, wat veel tijd kost en fouten introduceert.

De Oplossing: Een Nieuwe Manier van Kijken

De auteurs van dit paper hebben een slimme, open-source oplossing bedacht. Ze hebben twee krachtige gereedschappen samengevoegd die eerder apart werkten:

1. MFEM: Dit is als een superkrachtige, razendsnelle vrachtwagen (de "computermotor") die gebouwd is voor enorme taken. Hij kan rekenen op duizenden processoren tegelijk en is klaar voor de toekomstige kracht van grafische kaarten (GPUs). Hij is geweldig in het simuleren van stromend water.
2. FEBio: Dit is een meester-architect die alles weet over zachte materialen. Hij begrijpt precies hoe menselijk weefsel (zoals hartkleppen) reageert, rekken, en vervormen. Hij is de expert in de "dansende jurk".

De Magische Koppeling: Het "Geestelijke" Net

In plaats van de rivier en de jurk in één verwarrend net te stoppen, gebruiken ze een techniek die "geïmmersion" (onderdompeling) heet.

• De Analogie: Stel je voor dat je een zwembad hebt met een vast rooster op de bodem (het water-net). Je gooit nu een zachte, beweeglijke bal (de hartklep) in het zwembad. De bal beweegt door het rooster heen, maar het rooster zelf verandert niet.
• Het Probleem: Hoe zorg je dat het water weet dat de bal er is en dat het water niet door de bal heen stroomt?
• De Oplossing: De auteurs gebruiken een slimme wiskundige truc (een "distributie Lagrange-multiplicator"). Het is alsof er onzichtbare, elastische bandjes zijn die de bal en het water aan elkaar koppelen. Als het water tegen de bal duwt, voelt de bal het. Als de bal beweegt, duwt hij het water opzij.

Waarom is dit zo speciaal?

1. Geen Remeshing: Omdat het water-net vaststaat, hoeft de computer nooit te stoppen om het net opnieuw te bouwen, zelfs niet als de hartklep helemaal dichtklapt en weer opent. Het is alsof je een film maakt zonder dat de camera ooit hoeft te stoppen om de belichting aan te passen.
2. Krachtige Samenwerking: De "vrachtwagen" (MFEM) doet het zware rekenwerk voor het water, terwijl de "architect" (FEBio) precies weet hoe het weefsel zich gedraagt. Samen kunnen ze een hartklep simuleren die niet alleen beweegt, maar ook spanningen en rekken in het weefsel berekent. Dit is cruciaal om te begrijpen waarom een kunstklep na een paar jaar kan breken of waarom een aangeboren hartklep misvormd is.
3. Open Source: Dit is geen geheim recept voor een rijke farmaceutische firma. Het is gratis code die iedereen kan gebruiken, aanpassen en verbeteren. Dit helpt artsen en ingenieurs om betere operaties te plannen en betere medische hulpmiddelen te ontwerpen.

Wat hebben ze getest?

Ze hebben hun nieuwe systeem laten zien aan de hand van verschillende proeven:

• Een ring die in water drijft (om te zien of de wiskunde klopt).
• Een simpele klep die opent en sluit.
• Een vlag die wappert in de wind (om te zien of het systeem grote vervormingen aankan).
• Een bal die in water valt (om de zwaartekracht en weerstand te testen).
• De grote test: Een volledige, driedimensionale simulatie van een menselijk hart met een aortaklep, inclusief de bloedstroom en de spanningen in het weefsel.

Conclusie

Kortom, deze paper introduceert een nieuw, gratis en krachtig computerprogramma. Het combineert de snelheid van moderne supercomputers met de diepgaande kennis van weefselmechanica. Hierdoor kunnen wetenschappers voor het eerst heel gedetailleerd kijken naar hoe hartkleppen werken, hoe ze falen, en hoe we ze het beste kunnen repareren of vervangen, zonder dat de computer vastloopt door de complexe bewegingen. Het is een belangrijke stap naar betere zorg voor mensen met hartproblemen.

Technische samenvatting

Probleemstelling

De simulatie van vloeistof-structuurinteractie (FSI) in biologische systemen, met name voor hartkleppen, vormt een aanzienlijke computationele uitdaging. Traditionele methoden, zoals de Arbitrary Lagrangian-Eulerian (ALE) benadering, vereisen dat het vloeistofnetwerk zich aanpast aan de beweging van de vaste structuur. Bij grote vervormingen, contactmechanica (zoals het sluiten van hartklepbladen) en topologische veranderingen in het vloeistofgebied, leidt dit tot ernstige vervorming van het mesh. Dit vereist frequente herschikking van het mesh (remeshing), wat computatief duur is en fouten kan introduceren.

Bestaande open-source FSI-oplossers missen vaak geavanceerde solid-mechanica-modellering (zoals hyperelastische en visco-elastische constitutieve modellen voor zacht weefsel) of vereisen aanzienlijke gebruikersinspanning om deze te implementeren. Daarnaast zijn er beperkingen in de schaalbaarheid voor high-performance computing (HPC), met name met betrekking tot GPU-versnelling, die noodzakelijk is voor realistische, patiëntspecifieke simulaties.

Methodologie

Het artikel presenteert een nieuw, open-source framework dat twee gevestigde finite element-bibliotheken strategisch koppelt:

1. MFEM: Een schaalbare, GPU-bereide library ontwikkeld bij Lawrence Livermore National Laboratory, die zorgt voor de vloeistofoplosser en parallelle verwerking.
2. FEBio: Een niet-lineaire finite element-oplosser ontwikkeld bij de Universiteit van Utah en Columbia University, gespecialiseerd in biomechanica en geavanceerde solid-mechanica.

De Immersed Boundary Methode:
In plaats van een bewegend mesh te gebruiken, wordt de vaste structuur (bijv. de hartklep) "geïmmisceerd" in een vast achtergrondvloeistofdomein. De koppeling wordt bereikt via een fictitious domain/distributed Lagrange multiplier methode:

• De vloeistofvergelijkingen worden uitgebreid over het gehele domein, inclusief het gebied waar de vaste stof zich bevindt.
• Een Lagrange-multiplicator wordt gebruikt om de snelheid van de vloeistof en de vaste stof in het overlappende gebied gelijk te houden.
• Variational Multiscale (VMS) Stabilisatie: Om de nauwkeurigheid te verbeteren op onderopgeloste meshes (wat typisch is voor niet-aangepaste meshes), wordt VMS gebruikt. Een specifieke aanpassing van de stabilisatiecoëfficiënten (via een schalingsfactor $s$) wordt ingevoerd om grote drukstappen over de vloeistof-structuur-interface (zoals bij het sluiten van een klep) nauwkeurig te modelleren.
• Monolithische Koppeling: Het systeem wordt volledig impliciet en monolithisch opgelost (alle onbekenden tegelijk), wat robuustheid biedt voor sterk gekoppelde interacties. Door een speciale keuze voor de Lagrange-multiplicator basisfuncties, kan de multiplicator worden geëlimineerd, waardoor het lineaire systeem dezelfde grootte heeft als een puur vloeistofprobleem, wat de rekentijd verlaagt.

Implementatie:
De koppeling is geïmplementeerd als een plugin voor FEBio (de MFEMiFSI plugin). Deze plugin fungeert als een brug die datastructuren converteert en de monolithische matrix opbouwt met bijdragen van zowel de MFEM vloeistofoplosser als de FEBio solid-oplosser.

Belangrijkste Bijdragen

1. Unieke Synergie: Voor het eerst worden de schaalbare parallelle capaciteiten van MFEM (inclusief GPU-potentieel) gekoppeld aan de uitgebreide biomechanische modellering van FEBio. Dit stelt onderzoekers in staat om complexe weefselmechanica te combineren met hoge prestaties.
2. Open Source Framework: Het biedt een volledig open-source platform dat specifiek is ontworpen voor grote vervormingen en contactproblemen in de cardiovasculaire biomechanica, een gat in de huidige markt van open-source software.
3. Modulaire Architectuur: De plugin-architectuur maakt het eenvoudig om aanvullende fysica (zoals transport of chemische reacties) en geavanceerde elementtechnieken (zoals de "mean dilatation" methode voor bijna-incompressible materialen) toe te voegen zonder de kern van het framework aan te passen.
4. Geavanceerde Stabilisatie: De implementatie van aangepaste VMS-stabilisatiecoëfficiënten verbetert de nauwkeurigheid bij het modelleren van drukstappen over dunne structuren, wat cruciaal is voor hartklepsimulaties.

Resultaten

De auteurs hebben het framework gevalideerd en getest via een reeks numerieke benchmarks:

• Statisch evenwicht van een ring: Toonde optimale convergentiepercentages voor snelheid en druk, vergelijkbaar met analytische oplossingen.
• Idealiserde hartklep (Open en Gesloten):
  • Open klep: Toonde convergentie van de bladverplaatsing naar een referentie-ALE-oplossing.
  • Gesloten klep: Demonstreerde het vermogen om drukdiscontinuïteiten correct te vangen, wat essentieel is voor het sluiten van kleppen.
• Oscillerend flexibel blad: Simuleerde grote vervormingen gedreven door vloeistofkrachten, een scenario dat met ALE-methoden vaak onmogelijk is zonder herhaaldelijk remeshing.
• Vrij vallende bol: Bevestigde de nauwkeurigheid van de 3D-implementatie door de terminale snelheid te vergelijken met een analytische oplossing (met correctie voor wandeffecten).
• 3D Semilunaire Hartklep: De meest complexe test simuleerde een volledige hartcyclus van een tri-leaflet klep in een geïdealiseerde aorta. De simulatie toonde complexe vortex-dynamica, jet-vorming bij opening en het sluitmechanisme, inclusief de spanningen in het weefsel.

Betekenis en Toekomstperspectief

Dit werk vult een kritieke leemte in de computationele biomechanica. Het stelt onderzoekers in staat om patiëntspecifieke simulaties uit te voeren die zowel de hemodynamiek als de mechanische spanningen in het weefsel tegelijkertijd oplossen. Dit is essentieel voor:

• Het optimaliseren van reparatietechnieken voor aangeboren hartafwijkingen bij kinderen.
• Het voorspellen van de duurzaamheid en vermoeiing van bioprosthetische kleppen.
• Het evalueren van transcatheter apparaten.

De modulariteit en de integratie met MFEM leggen de basis voor toekomstige ontwikkelingen, waaronder volledige GPU-versnelling, adaptief mesh-verfijning (AMR) rond de interface, en de integratie van contactmodellen en schaalmodellen (shell formulations) voor nog efficiëntere simulaties. De code is beschikbaar via de FEBio Plugin Repository, wat de adoptie en uitbreiding door de gemeenschap faciliteert.