A Modified Suspension-Balance Model for Deformable Particle Suspensions: Application to Blood Flows with Cell-Free Layer

Dit artikel stelt een aangepast suspensie-balansmodel voor dat hydrodynamische liftkrachten incorporeert om bloedstromen in microvasculaire kanalen efficiënt te simuleren, waarbij het succesvol kernfenomenen zoals de vorming van een celvrije laag, hematocriet- en snelheidsprofielen, en de Fahraeus- en Fahraeus-Lindqvist-effecten vastlegt.

De kern

Stel je een drukke snelweg voor waar de auto's eigenlijk piepkleine, zachte rode bloedcellen zijn, en de weg een microscopisch bloedvat is. In een normale verkeerssituatie zou je verwachten dat het verkeer gelijkmatig verdeeld is. Maar in ons lichaam hebben deze "auto's" een vreemde gewoonte: ze haten het om dicht bij de wanden te zijn. Ze geven er de voorkeur aan om in het midden van de weg bij elkaar te kruipen, waardoor er een duidelijke, lege rijstrook direct naast het wegdek overblijft.

Deze lege rijstrook wordt de Cell-Free Layer (CFL) genoemd. Dit is een cruciaal kenmerk van onze bloedstroom dat helpt om ons bloed sneller te laten stromen en met minder wrijving.

Het Probleem: De Oude Kaart Mist een Afslag

Wetenschappers proberen al jaren computermodellen te bouken om te simuleren hoe bloed stroomt. Ze gebruiken hiervoor iets dat een "Suspension Balance Model" (SBM) wordt genoemd. Denk aan dit model als een verkeerssimulatiesoftware.

De oude versie van deze software was goed in het voorspellen dat auto's naar het midden van de weg zouden bewegen vanwege de manier waarop ze tegen elkaar aan botsen. Echter slaagde het er niet in uit te leggen waarom de auto's zo graag de wanden willen verlaten. Het kon die lege "cell-free lane" naast de rand niet creëren. Het was als een GPS die wel wist dat de auto's bewogen, maar niet wist dat ze actief de stoeprand probeerden te vermijden.

De Oplossing: Een Nieuwe "Push"-knop

De auteurs van dit artikel, onder leiding van Hugo Castillo-Sánchez en Leonardo Liu, besloten de software te repareren. Ze realiseerden zich dat omdat rode bloedcellen zacht zijn (vervormbaar), ze een speciale soort onzichtbare kracht genereren wanneer ze te dicht bij een wand komen.

Ze noemen dit de Lift Force (opliftkracht).

• De Analogie: Stel je voor dat je zwemt vlak bij de zijkant van een zwembad. Terwijl je beweegt, duwt het water je een beetje weg van de wand. Voor rode bloedcellen is deze "duw" veel sterker omdat ze zacht zijn en van vorm veranderen terwijl ze zich langs de wand perst.
• De Fix: Het team heeft deze "Lift Force" toegevoegd aan hun computermodel. Ze hebben een Modified Suspension Balance Model (MSBM) gecreëerd. Nu kijkt de software niet alleen naar de auto's; de software duwt ze actief weg van de wand, net zoals het water een zwemmer wegduwt.

Wat gebeurde er toen ze de simulatie draaiden?

Toen ze deze nieuwe "Lift Force" inschakelden in hun computer, veranderden de resultaten drastisch:

1. De lege rijstrook verscheen: De simulatie creëerde succesvol die duidelijke zone nabij de wand (de CFL) die we in het echte leven zien.
2. De verkeersopstopping in het midden: De rode bloedcellen stapelden zich op in het centrum, wat een dichte kern creëerde.
3. De vorm van de stroming: Omdat de cellen in het midden gebundeld waren en de randen vrij waren, stroomde het bloed niet in een vloeiende, gebogen boog (zoals een normale rivier). In plaats daarvan stroomde het als een solide plug of een zuiger, met een platte bovenkant. Dit is precies wat er gebeurt in echte microvaten.

Het testen van het nieuwe model

Het team heeft niet alleen gegokt; ze hebben hun nieuwe model getest tegen echte wereldgegevens en andere complexe simulaties:

• Tijdreizen: Ze keken hoe de "lege rijstrook" in de loop van de tijd ontstond. Het begon met cellen overal, en langzaam duwde de "Lift Force" ze weg van de wanden totdat de rijstrook leeg was. Dit kwam overeen met de snelheid en het gedrag dat wordt gezien in experimenten met hogesnelheidscamera's.
• Het "Fåhræus-effect": Dit is een chique term voor een eenvoudige observatie: bloed stroomt sneller in kleine buisjes dan je zou verwachten, en de concentratie van cellen in het midden is anders dan de concentratie bij de uitgang. Hun nieuwe model voorspelde dit perfect.
• Het "Fåhræus-Lindqvist-effect": Dit is de observatie dat bloed "dunner" (minder plakkerig) wordt wanneer het door zeer kleine buisjes stroomt. Hun model legde dit ook vast, waarbij het liet zien dat de lege laag nabij de wand de wrijving vermindert, waardoor het bloed makkelijker stroomt.

De Kern van het Verhaal

Het artikel beweert dat door een eenvoudige "duw" (de liftkracht) toe te voegen aan hun computermodel, ze nu nauwkeurig kunnen simuleren hoe bloed zich in minuscule vaten gedraagt.

• Wat het doet: Het legt de vorming van de cell-free layer vast, de plug-achtige stroming en de beroemde "Fåhræus"-effecten die ervoor zorgen dat de bloedstroom efficiënt is in ons lichaam.
• Wat het (nog) niet doet: De auteurs geven toe dat voor zeer grote buisjes (groter dan 40 micrometer) het model de cellen een beetje te ver weg duwt. Ze vermoeden dat dit komt omdat hun model nog niet rekening houdt met de manier waarop cellen elkaar "afschermen" van de wand wanneer ze dicht op elkaar gepakt zitten. Ze zijn van plan dit in toekomstig werk te herstellen.

Kortom, ze hebben een betere digitale tweeling gebouwd voor bloedstroom die begrijpt dat rode bloedcellen niet alleen passieve passagiers zijn, maar actieve zwemmers die zichzelf weg van de wanden duwen om de snelweg vrij te houden.

Technische samenvatting

Technische Samenvatting: Een Gemodificeerd Suspensiebalansmodel voor Deformeerbare Deeltjessuspensies

Probleemstelling
Het modelleren van de stroming van rode bloedcellen (RBC's) en andere deformeerbare deeltjessuspensies in beperkte geometrieën vormt aanzienlijke uitdagingen vanwege inherente niet-lineariteit, multi-body interacties en niet-Newtonia true rheologische gedragingen. Een cruciaal fenomeen in de microcirculatie is de vorming van een celvrije laag (CFL), een nabij de wand gelegen regio die uitgearmd is van RBC's. Deze laag is fysiologisch essentieel voor het handhaven van een lage schijnbare viscositeit en het bevorderen van plaatjesmargenatie. Hoewel computationele modellen op celniveau de CFL-vorming kunnen vastleggen, zijn deze vaak computationeel te kostbaar voor simulaties van vaatnetwerken of orgaan-schaal. Daartegenover staan bestaande continuümmodellen, zoals het conventionele Suspensiebalansmodel (SBM), die vaak niet in staat zijn de CFL te vangen omdat ze de door deformabiliteit geïnduceerde hydrodynamische liftkrachten, die deeltjes weg van de wanden drijven, niet adequaat rekening houden.

Methodologie
De auteurs stellen een Gemodificeerd Suspensiebalansmodel (MSBM) voor dat het conventionele SBM-raamwerk aanvult met een specifieke liftkrachtterm om rekening te houden met hydrodynamische interacties tussen wand en RBC.

• Theoretisch Raamwerk: Het model behoudt de standaard massa- en momentumconservatievergelijkingen voor de suspensie-bulk en de deeltelfases. De kernmodificatie ligt in de momentumbalans van de deeltelfase. Het conventionele SBM houdt rekening met deeltjesmigratie gedreven door de divergentie van de deeltelfase-spanning ($\nabla \cdot \langle \Sigma \rangle_p$) en gehinderde bezinkingsweerstand. De MSBM introduceert een extra volumegemiddelde liftkrachtterm, $\phi L_\perp$, afgeleid van experimentele observaties van deformeerbare vesikels en RBC's.
• Formulering van de Liftkracht: De liftkracht wordt gemodelleerd als een functie van de lokale afschuifsnelheid ($\dot{\gamma}$), de afstand tot de wand ($h$), en een gereduceerde volumeparameter ($\nu$). De term heeft de vorm $L_\perp \propto \frac{\dot{\gamma}}{(h+h_0)^\beta}$, waarbij $\beta$ een machtswet-coëfficiënt is die wordt aangepast om de deeltjesafhankelijkheid te reflecteren (specifiek voor RBC's), en $h_0$ een geometrische factor is om singulariteit bij de wand te vermijden.
• Numerieke Implementatie: De MSBM werd geïmplementeerd in OpenFOAM, waarbij de bestaande SbmFoam solver is uitgebreid tot SbmBloodFoam. De solver maakt gebruik van een eindige-volumemethode (FVM). Om het continuümmodel te valideren, hebben de auteurs een experimenteel gevalideerde multiscale bloedstroomsolver gebruikt, gebaseerd op de Lattice Boltzmann-methode (LB) en grofkorrelige deformeerbare capsulemodellen voor RBC's.

Belangrijkste Resultaten
De studie demonstreert het vermogen van de MSBM om bloedstromingen in microvasculaire kanalen en buizen te simuleren, waarbij zowel transiënte als stationaire gedragingen worden gevangen:

1. CFL-vorming en Snelheidsprofielen: In tegen tegenstelling tot het conventionele SBM, dat een parabolisch snelheidsprofiel en een niet-nul concentratie bij de wand voorspelt, legt de MSBM succesvol de vorming van een CFL en het resulterende "plug-flow"-achtige snelheidsprofiel dat kenmerkend is voor bloed vast. Het model laat zien dat het verhogen van de liftkrachtparameter ($\beta$) de dikte van de CFL vergroot.
2. Parametrische Gevoeligheid:
   • Bulk Hematocriet ($\phi_b$): Naarmate $\phi_b$ toeneemt, neemt de dikte van de CFL af, wat consistent is met de literatuur.
   • Maximale Pakkingsconcentratie ($\phi_m$): Hogere $\phi_m$-waarden (die een grotere RBC-aggregatie of pakkingsneiging vertegenwoordigen) leiden tot dikkere CFL's en hogere centrale hematocrietpieken.
3. Temporele Evolutie: Het model legt de transiënte ontwikkeling van de CFL vast. Aanvankelijk is de hematocriet uniform; na verloop van tijd migreren RBC's naar het centrum, waardoor een cel-gearmeerde laag nabij de wand en een geconcentreerde kern ontstaat, om uiteindelijk een stationaire toestand te bereiken.
4. Validatie tegenover Theorie en Experimenten:
   • Eenvoudige Afschuivingsstroom: Het model reproduceert de lineaire schaling van de kubische CFL-dikte ($\delta_{CFL}^3$) met dimensieloze tijd ($t\dot{\gamma}$) op korte tijdschalen, wat overeenkomt met de schaleringswetten waargenomen in experimenten met enkelvoudige deeltjesmigratie.
   • Kanaalstromingen: De MSBM-voorspellingen voor snelheid en CFL-dikte komen goed overeen met experimentele gegevens van Losserand et al. en Salame over diverse kanaalhoogtes en bulkhematocrieten.
   • Buisstromingen: Het model legt succesvol het Fåhræus-effect (waarbij de discharge hematocriet de buishematocriet overtreft) en het Fåhræus-Lindqvist-effect (afname van de schijnbare viscositeit bij afnemende vaatdiameter) vast. De voorspelde schijnbare viscositeit en hematocrietratio's vertonen een goede overeenstemming met empirische correlaties (bijv. Pries et al.) en experimentele gegevens, met name voor vaatdiameters onder de 40 µm.

Betekenis en Claims
Het artikel vestigt een nieuw continuüm computationeel raamwerk dat efficiënt de microstructurele heterogeniteit en het niet-Newtonia true stromingsgedrag van geconcentreerde deformeerbare deeltjessuspensies onder opsluiting kan vastleggen. De auteurs beweren dat door een deformabiliteit-geïnduceerde liftkracht op te nemen, de MSBM de kloof overbrugt tussen computationeel dure discrete deeltjessimulaties en minder nauwkeurige continuümmodellen.

Het werk wordt gepresenteerd als een stap naar efficiënte, grootschalige biomedische en industriële toepassingen. De auteurs merken bescheiden op dat hoewel het model goed presteert voor kleinere vaten, het de CFL en het Fåhræus-Lindqvist-effect voor grotere buizen (>40 µm) licht overschat. Zij schrijven dit toe aan het gebrek aan een "deeltjes-afschermingseffect" in de huidige liftkrachtformulering, wat zou kunnen verklaren waarom liftkrachten afnemen in deeltjesrijke regio's door lokale deeltje-wand smeringsinteracties. Deze beperking wordt geïdentificeerd als een onderwerp voor toekomstig onderzoek in plaats van een huidig falen van het model.