Large-eddy simulation of the FDA benchmark blood pump: validation against experiments and implications for turbulent flow mechanisms

Deze studie valideert dat Large Eddy Simulation (LES) met tijdsafhankelijke rotor-statorkoppeling de experimentele stromingsmetingen in de FDA-benchmarkbloedpomp nauwkeuriger voorspelt dan RANS-methoden, en bevestigt dat deze schaaloplossende benadering essentieel is voor het correct modelleren van de complexe turbulente stromingsmechanismen in hartventrikelhulpmiddelen.

De kern

Hartpomp-simulaties: Waarom we de "microscopische chaos" moeten zien om bloed veilig te houden

Stel je voor dat je een hartpomp ontwerpt die een menselijk hart kan helpen. Dit is geen gewone waterpomp; het is een ultra-geavanceerde machine die bloed door het lichaam moet pompen zonder het bloed te beschadigen. Als de stroming binnenin te chaotisch is, kan het bloedcellen kapotmaken, wat dodelijk kan zijn.

De onderzoekers in dit paper hebben gekeken naar hoe computersimulaties (CFD) deze stroming voorspellen. Ze hebben een beroemde "testcase" gebruikt, ontworpen door de Amerikaanse gezondheidsautoriteit (FDA), om te zien welke rekenmethode het beste werkt.

Hier is de uitleg in gewone taal, met een paar creatieve vergelijkingen:

1. Het probleem: De computer ziet niet alles

Om te weten of een pomp veilig is, moeten we weten hoe het bloed stroomt. Computers proberen dit na te bootsen.

• De oude methode (RANS): Stel je voor dat je probeert het weer te voorspellen door alleen naar het gemiddelde weer van de afgelopen maand te kijken. Je ziet dat het vaak regent, maar je mist de plotselinge onweersbuien. Zo werkt de oude methode: hij berekent het gemiddelde, maar mist de snelle, chaotische wervelingen (turbulentie) die echt belangrijk zijn voor het beschadigen van bloedcellen.
• De nieuwe methode (LES): Dit is alsof je een super-snelheidscamera gebruikt die elke seconde van de onweersbui vastlegt. Je ziet precies hoe de wind draait, waar de turbulentie zit en hoe de krachten werken. Dit is wat de onderzoekers hebben gebruikt: Large Eddy Simulation (LES).

2. De test: De "FDA-bloeddrukmeter"

De onderzoekers hebben hun computermodel getest tegen echte experimenten. Ze gebruikten een doorzichtig model van een hartpomp en lieten er vloeistof doorheen stromen terwijl ze met lasers (een techniek genaamd PIV) de snelheid van de deeltjes maten.

• Het resultaat: De oude methode (RANS) gaf in de "uitlaat" van de pomp (de diffuser) een heel verkeerd beeld. Het was alsof ze dachten dat de vloeistof rustig stroomde, terwijl de camera liet zien dat het daar een wild, draaiend kolkend water was.
• De nieuwe methode (LES) met een dynamische koppeling (waarbij de draaiende delen echt bewegen in de simulatie) gaf een beeld dat bijna perfect overeenkwam met de echte metingen. Het zag precies waar de vloeistof versnelde en waar de wervelingen ontstonden.

3. De puzzel: Hoe groot moet het raster zijn?

Om zo'n gedetailleerde simulatie te doen, moet je het binnenste van de pomp opdelven in miljoenen kleine blokjes (een "mesh").

• Te grof (10 miljoen blokjes): Dit is alsof je een foto maakt met een oude, korrelige camera. Je ziet de grote lijnen, maar de fijne details van de turbulentie gaan verloren. De simulatie is dan net op de rand van bruikbaar.
• Net goed (51 miljoen blokjes): Hier begin je de chaos goed te zien. Het is een goede balans tussen kosten en kwaliteit.
• Zeer fijn (80 miljoen blokjes): Dit is een 8K-foto. Je ziet elke kleine wervel. De onderzoekers ontdekten dat dit nodig is om de kleinste details te zien die bloedcellen kunnen beschadigen, maar voor het gemiddelde stromingsbeeld is dit misschien al iets te veel werk (te duur in rekentijd).

4. Wat gebeurt er eigenlijk binnenin?

Met de goede simulatie konden ze zien wat er echt gebeurt:

• De wervel-dans: Het bloed wordt niet rustig rondgeduwd. De lepelvormige bladen van de pomp creëren enorme, draaiende wervels (zoals kleine tornado's).
• De gevaarlijke plekken: Waar deze wervels botsen of waar het bloed heel snel versnelt, ontstaat er veel turbulentie. Dit is waar het bloed het meeste risico loopt om beschadigd te worden.
• De verrassing: Zelfs bij snelheden waar men dacht dat het nog "rustig" was, bleek er al veel turbulentie te zijn. De oude methodes zagen dit niet, maar de nieuwe methode wel.

5. Waarom is dit belangrijk?

Als artsen en ingenieurs hartpompen ontwerpen, willen ze voorkomen dat het bloed kapotgaat (hemolyse).

• Als je de oude methode gebruikt, denk je misschien: "Ah, de stroming is veilig," terwijl de computer eigenlijk de gevaarlijke wervels heeft gemist.
• Met deze nieuwe, gedetailleerde methode (LES) kunnen ze zien precies waar de gevaarlijke plekken zitten. Hierdoor kunnen ze de pomp zo ontwerpen dat het bloed veilig blijft, wat levens kan redden.

Kortom:
Dit onderzoek zegt: "Stop met het kijken naar het gemiddelde. Om hartpompen veilig te maken, moeten we de computer laten kijken naar de echte, chaotische dans van het bloed, en daarvoor hebben we heel fijne rekenmodellen nodig." Het is een stap van "ruwe schatting" naar "precieze voorspelling" voor de geneeskunde van de toekomst.

Technische samenvatting

Probleemstelling

Ventriculaire hulpmiddelen (VAD's), zoals centrifugale bloedpompen, zijn cruciaal voor de behandeling van hartfalen. Het ontwerp en de optimalisatie van deze pompen vereisen een nauwkeurige voorspelling van de hemodynamiek, aangezien turbulentie en schuifspanningen direct gerelateerd zijn aan bloedschade (hemolyse en trombose).

• Uitdaging: Centrifugale pompen werken vaak in een Reynolds-getalbereik ($10^5$–$10^6$) dat dicht bij de turbulentiedrempel ligt, wat leidt tot complexe, overgangs- en turbulentestromingen.
• Huidige beperkingen: De meeste bestaande CFD-studies (Computational Fluid Dynamics) voor bloedpompen vertrouwen op Reynolds-averaged Navier-Stokes (RANS) benaderingen. Deze methoden modelleren alle turbulentie en kunnen de tijdsafhankelijke interacties tussen rotor en stator (bijv. bladpassage-effecten) niet adequaat vastleggen.
• FDA Benchmark: Ondanks uitgebreide experimentele data van de FDA (Particle Image Velocimetry of PIV) voor een gestandaardiseerd centrifugaal pompmodel, hebben eerdere CFD-studies moeite gehad om de snelheidsvelden consistent en nauwkeurig te voorspellen, vooral in de diffuser-regio waar sterke turbulentie en intermittentie optreden.

Methodologie

De auteurs voeren een systematische validatie en vergelijking uit van verschillende CFD-strategieën voor het FDA-benchmarkmodel:

1. Simulatiebenaderingen:
   • LES (Large Eddy Simulation): Een schaal-resolverende methode die de energiebevattende turbulentestructuren direct oplost en alleen de subgrid-schaal modelleert. Er wordt gebruikgemaakt van het WALE (Wall-Adapting Local Eddy-viscosity) subgrid-model.
   • RANS & URANS: Vergelijking met traditionele RANS (steady) en Unsteady RANS (URANS) methoden, gebruikmakend van het k-ω SST turbulentiemodel.
2. Rotor-Stator Koppeling:
   • SI (Sliding Interface): Een tijdsafhankelijke methode waarbij het roterende en stationaire domein fysiek bewegen ten opzichte van elkaar, waardoor tijdsafhankelijke interacties expliciet worden opgelost.
   • MRF (Multiple Reference Frame): Een steady-state benadering waarbij het roterende domein in een roterend referentiekader wordt opgelost zonder expliciete beweging van het rooster.
3. Netwerk (Mesh) en Validatie:
   • Drie mesh-resoluties werden gebruikt: ~10 miljoen, ~51 miljoen en ~80 miljoen cellen.
   • De simulaties werden gevalideerd tegen PIV-metingen onder twee operationele condities (Conditie 2 en 5), met name in het bladvlak en de diffuser.
   • Kwaliteitsmetriek: Drie complementaire metrieken werden toegepast om de kwaliteit van de LES te beoordelen:
     1. Opgelost TKE-quotiënt ($I_{Q_k}$): De fractie van de totale turbulente kinetische energie die direct wordt opgelost.
     2. Subgrid-activiteitsparameter ($s$): Vergelijking van subgrid-dissipatie met viskeuze dissipatie.
     3. $I_{Q_\nu}$ index: Een alternatieve index gebaseerd op de verhouding van turbulentieviscositeit tot viskeuze viscositeit.

Belangrijkste Resultaten

1. Validatie tegen Experimenten:

   • LES met SI toont een aanzienlijk betere overeenkomst met experimentele PIV-data dan RANS-methoden, vooral in de diffuser-regio waar sterke intermittentie en wandgebonden turbulentie optreden.
   • RANS-MRF vertoont aanzienlijke afwijkingen, waaronder niet-fysische oscillaties in het bladvlak en een onderschatting van de snelheid in de tip-clearance.
   • URANS-SI presteert beter dan RANS-MRF, maar faalt nog steeds om de complexe stromingsprofielen in de diffuser volledig vast te leggen vergeleken met LES.

2. Mesh-gevoeligheid en Kwaliteit:

   • Een mesh van ~80 miljoen cellen (Mesh 3) bereikt een $I_{Q_k}$ van >95%, wat aangeeft dat het een goed opgeloste LES-regime is.
   • De fijnere mesh (Mesh 2, 51 miljoen) levert ook voldoende kwaliteit (>91%), terwijl de grove mesh (Mesh 1, ~10 miljoen) slechts marginaal opgelost is (80%) en tekortschiet in kritieke gebieden met hoge turbulentie.
   • De resultaten suggereren dat voor betrouwbare voorspellingen van gemiddelde stromingsgrootheden een mesh van ~50 miljoen cellen voldoende kan zijn, maar voor het vastleggen van fijne schaalstructuren (relevant voor bloedschade) de 80 miljoen cellen nodig zijn.

3. Fysische Inzichten in Turbulentie:

   • De analyse van wervelstructuren (via de Q-criterium) toont aan dat turbulentie in centrifugale pompen wordt gedreven door een complex samenspel van blad-induceren wervels, wake-structuren, wandwervels en jet-stromingen.
   • Er is een sterke correlatie tussen gebieden met hoge turbulente kinetische energie (TKE) en intense werveldynamica.
   • Spectrale analyse: De snelheidsspectra tonen een duidelijk $-5/3$ bereik (inertie-subbereik), wat bevestigt dat de LES de grootschalige turbulentie correct oplost. De kleinste dissipatieschalen blijven onopgelost, wat consistent is met de LES-aanpak.

Bijdragen en Betekenis

• Validatie van LES: Het artikel biedt de eerste uitgebreide validatie van tijdsafhankelijke LES met Sliding Interface voor het FDA-benchmarkpompmodel, bewijzend dat deze methode superieur is aan RANS voor het vastleggen van de complexe, tijdsafhankelijke stromingsfysica in bloedpompen.
• Richtlijnen voor Mesh-resolutie: De studie biedt kwantitatieve richtlijnen voor toekomstige studies: een mesh van ~10 miljoen cellen is onvoldoende voor betrouwbare turbulentie-analyse, terwijl ~50-80 miljoen cellen nodig zijn voor hoogwaardige hemodynamische studies.
• Implicaties voor Bloedschade: Omdat bloedschade sterk afhankelijk is van lokale schuifspanningen en tijdsafhankelijke turbulentie, suggereert de studie dat RANS-benaderingen onvoldoende zijn voor nauwkeurige hemolyse-voorspellingen. Transiënte, schaal-resolverende simulaties zijn essentieel voor het ontwikkelen van veiligere VAD's.
• Fysisch Mechanisme: De studie onthult dat werveldynamica de centrale rol speelt in de turbulentiegeneratie binnen centrifugale pompen, wat een fundamenteel inzicht biedt voor het optimaliseren van pompopstellingen om bloedschade te minimaliseren.

Kortom, dit onderzoek onderstreept de noodzaak van hoogwaardige, tijdsafhankelijke CFD-simulaties (LES) in plaats van traditionele RANS-methoden voor het nauwkeurig modelleren van bloedpompen, en biedt een blauwdruk voor het balanceren van rekenkosten en voorspellingsnauwkeurigheid.