Modelling Material Injection Into Porous Structures Under Non-isothermal Conditions

Dit artikel presenteert een thermodynamisch consistent model op basis van de theorie van poreuze media voor de injectie van acrylcement in wervelbot onder niet-isotherme omstandigheden, waarbij lokale thermische niet-evenwichtstoestanden worden beschouwd om de temperatuurverschillen tussen het cement en het menselijk lichaam nauwkeuriger te simuleren.

De kern

Het Verwarmen van een Spons: Een Simpel Verhaal over Ruggegraatsinjecties

Stel je voor dat je een oude, beschadigde rugwervel moet repareren. De wervel is als een brokachtige spons die niet meer stevig genoeg is om het lichaam te dragen. Artsen lossen dit op met een procedure genaamd vertebroplastie: ze spuiten een vloeibare, acryl-cement in de wervel, die daar uithardt en de wervel weer verstevigt.

Maar hier zit een probleem: het cement dat de arts inspuit is vaak kouder dan het menselijk lichaam (ongeveer 37 graden). Als je koude vloeistof in een warme spons spuit, ontstaat er een klein temperatuurverschil. Bovendien reageert het cement chemisch en wordt het iets warmer tijdens het harden.

De auteurs van dit paper (Jan-Sören, Zainab en hun collega's) hebben een wiskundig model gemaakt om precies te begrijpen wat er gebeurt in die wervel tijdens zo'n injectie. Ze willen weten: Hoe stroomt het cement? Hoe verspreidt de warmte zich? En is het veilig?

Hier is hoe ze dat hebben aangepakt, vertaald in alledaagse taal:

1. De "Geest van de Spons" (De Theorie)

Stel je de wervel niet voor als een steen, maar als een drijvende stad in een dichte mist.

• De gebouwen: Dat is het bot (de vaste stof).
• De lucht: Dat is het beenmerg (een vloeistof die al in de wervel zit).
• De nieuwe bewoners: Dat is het cement dat erin wordt gepompt.

Vroeger maakten wetenschappers modellen die aannamen dat alles even warm was (als een kamer met een perfecte verwarming). Maar in werkelijkheid is het cement kouder dan het bot, en het bot is warmer dan het merel. Het model van deze auteurs is speciaal gemaakt om deze temperatuurverschillen mee te nemen. Ze kijken naar hoe warmte "huilt" van het ene materiaal naar het andere, net zoals warmte van een hete kop koffie naar je koude hand gaat.

2. De Drie Spelregels (De Wiskunde)

Om dit te simuleren, gebruiken ze drie hoofdregels, alsof ze een spelletje spelen met drie spelers:

1. Massa: Hoeveel cement en merel is er? (Niets verdwijnt zomaar).
2. Kracht: Hoe duwt het cement het merel opzij? (Het cement moet de merel verdringen om ruimte te maken).
3. Energie (Warmte): Hoe warmt het cement op en koelt het bot af?

Ze hebben een slimme truc bedacht: ze nemen aan dat de temperatuur van het cement, het merel en het bot op elk punt evenwichtig kan zijn, maar dat ze ook even anders kunnen zijn. Dit noemen ze "lokale thermische niet-evenwicht". Klinkt ingewikkeld, maar het is alsof je in een drukke trein zit: iedereen heeft een eigen temperatuur, maar ze wisselen warmte uit met elkaar.

3. De Twee Experimenten (De Simulaties)

Om hun model te testen, hebben ze twee scenario's doorgespeeld in een virtuele buis (die de wervel voorstelt):

• Scenario 1: De "Lekker Warme" Injectie.
  Het cement is precies even warm als het bot (35°C). Hierdoor is er geen temperatuurverschil. Het model laat zien hoe het cement de merel verdringt. Het resultaat? Het cement stroomt als een golf vooruit. Als ze een bepaalde wiskundige formule gebruiken (de "Brooks-Corey" formule), ontstaat er een scherpe grens (een "schokgolf") tussen het cement en het merel.

• Scenario 2: De "Koude" Injectie.
  Het cement is kouder (25°C) dan het bot. Nu gebeurt er iets interessants: de koude golf van het cement beweegt door de warme wervel.

  • Het cement koelt het bot direct om zich heen iets af.
  • Het merel koelt ook iets af, maar iets later (want het bot fungeert als een buffer).
  • Uiteindelijk warmt alles weer op door de warmte van het lichaam.

4. Wat Leerden Ze?

De resultaten waren geruststellend:

• Geen ramp: De temperatuurverschillen die ze zagen waren heel klein (slechts een paar tienden van een graad). In de echte wereld is dit waarschijnlijk verwaarloosbaar, maar het is goed om te weten dat hun model dit wel kan zien.
• De druk: Als het cement kouder is, moet de arts iets meer druk uitoefenen om het erin te krijgen, maar het verschil is klein.
• De veiligheid: Het model bevestigt dat hun wiskundige aanpak "logisch" is. Het gedraagt zich zoals de natuur het zou doen (thermodynamisch consistent).

Conclusie in Eén Zin

Deze wetenschappers hebben een digitale proefomgeving gebouwd die niet alleen kijkt naar hoe cement door bot stroomt, maar ook hoe de temperatuur daarbij een rol speelt. Het is alsof ze een simulator hebben gemaakt voor een chirurg, zodat die in de toekomst nog beter kan voorspellen wat er in een patiënt gebeurt, zonder dat ze het echt hoeven te proberen.

Hoewel ze in dit specifieke artikel zeggen dat de temperatuurverschillen klein zijn, is het belangrijk dat ze hun model zo hebben gebouwd dat het wel die verschillen kan zien als het nodig is. Dat maakt hun werk een sterke basis voor toekomstige, nog veiligere operaties.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het paper "Modelling Material Injection Into Porous Structures Under Non-isothermal Conditions" in het Nederlands.

Titel: Modellering van Materiaalinjectie in Poreuze Structuren onder Niet-isotherme Omstandigheden

Auteurs: J.-S. L. Völter et al. (Universiteit van Stuttgart & Ansbach University of Applied Sciences)

---

1. Probleemstelling

Het paper richt zich op de percutane vertebroplastie (PV), een medische procedure waarbij acrylcement wordt geïnjecteerd in beschadigde wervels om deze te stabiliseren. Hoewel er reeds macroscopische modellen bestaan die gebaseerd zijn op de Theorie van Poreuze Media (TPM), zijn deze tot nu toe beperkt tot isotherme aannames (constante temperatuur).

Dit is een beperking omdat:

1. De temperatuur van het geïnjecteerde cement doorgaans lager is dan de lichaamstemperatuur.
2. Het uithardingsproces van het cement exotherm is (warmte genereert).
3. Het negeren van temperatuurverschillen kan leiden tot onnauwkeurige voorspellingen van de injectie-dynamiek en mechanische interacties.

Het doel van dit werk is het uitbreiden van bestaande TPM-modellen naar een niet-isotherme case, waarbij rekening wordt gehouden met lokale thermische niet-evenwichtstoestanden (LTNE).

2. Methodologie

De auteurs hanteren de principes van de Theorie van Poreuze Media (TPM) voor een meerfasig continuüm.

• Fasen: Het model beschouwt drie fasen:
  • Vaste fase ($\phi_S$): Trabeculair bot.
  • Vloeibare fase 1 ($\phi_M$): Botmerg (natmakend fluïdum).
  • Vloeibare fase 2 ($\phi_C$): Botcement (niet-natmakend fluïdum).
• Thermodynamische Raamwerk:
  • Er worden drie energiebalansen afgeleid (één per fase) in plaats van één gezamenlijke balans.
  • De afleiding volgt de Coleman en Noll-procedure om thermodynamische consistentie te garanderen.
  • Er wordt uitgegaan van lokale thermische niet-evenwicht (LTNE), wat betekent dat de temperaturen van de vaste stof, het merg en het cement op elk punt kunnen verschillen ($T_S \neq T_M \neq T_C$).
• Constitutieve Relaties:
  • Mechanica: Het bot wordt gemodelleerd als een hyperelastisch materiaal (neo-Hookean model).
  • Stroming: De stroming wordt beschreven door filterwetten (Darcy-wet) met inachtneming van capillaire druk en relatieve permeabiliteit.
  • Warmtetransport: Er wordt gebruikgemaakt van de wet van Fourier voor warmtegeleiding en een term voor warmteoverdracht tussen de fasen (interface-uitwisseling).
  • Capillaire Druk: Hoewel capillaire krachten in de afleiding worden meegenomen, worden ze in de simulaties verwaarloosbaar klein gehouden om de focus op thermische effecten te houden.
• Numerieke Implementatie:
  • De vergelijkingen worden opgelost met een monolithische Finite Element Methode (FEM) (solver: PANDAS).
  • Ruimtelijke discretisatie: Petrov-Galerkin methode met Box-discretisatie (upwinding).
  • Tijdsdiscretisatie: Crank-Nicolson schema.

3. Belangrijkste Bijdragen

1. Thermodynamisch Consistent Model: Het paper presenteert een nieuw, thermodynamisch consistent model voor materiaalinjectie in poreuze media onder LTNE-omstandigheden, specifiek toegepast op vertebroplastie.
2. Afleiding van Energiebalansen: De auteurs leiden drie gekoppelde energiebalansen af die de warmtegeleiding binnen elke fase en de warmteoverdracht tussen de fasen beschrijven.
3. Validatie van Aannames: Het model toont aan dat, indien beperkt tot lokale thermische evenwicht (LTE), de vergelijkingen overeenkomen met bestaande TPM-modellen.
4. Vergelijking van Permeabiliteitsmodellen: Er wordt een vergelijking gemaakt tussen het gebruik van Brooks-Corey relatieve permeabiliteitsfactoren en lineaire permeabiliteitsfactoren, en hoe deze de injectiedruk en temperatuurverdeling beïnvloeden.

4. Resultaten

Twee simulatiescenario's werden uitgevoerd in een tubulair domein:

• Scenario 1: Injectie van cement op lichaamstemperatuur (308,15 K) in een omgeving met dezelfde temperatuur (alleen dissipatie-effecten).
• Scenario 2: Injectie van kouder cement (298,15 K) in een warmer bot (308,15 K).

Kernbevindingen:

• Temperatuurverdeling: In Scenario 2 ontstaat een niet-uniforme temperatuurverdeling. De temperatuur van het cement verandert het snelst, gevolgd door het bot en vervolgens het merg. De temperatuurverschillen zijn echter klein (maximaal ~0,93 K na 120 seconden), wat suggereert dat voor deze specifieke injectiesnelheid de LTE-aanname (lokale thermische evenwicht) in de praktijk mogelijk voldoende is, hoewel het LTNE-model fysisch correcter is.
• Drukgedrag:
  • Bij gebruik van Brooks-Corey permeabiliteit: De injectiedruk neemt significant toe over de tijd, zelfs bij constante viscositeit.
  • Bij gebruik van lineaire permeabiliteit: De injectiedruk blijft constant over de tijd.
  • Opmerking: Experimentele data suggereert geen sterke drukstijging, maar het model negeert niet-Newtonse viscositeitseigenschappen van cement, wat de vergelijking bemoeilijkt.
• Stromingsdynamiek: Met Brooks-Corey factoren ontstaat een schokfront gevolgd door een verdunningswaaier in de verzadiging. Met lineaire factoren ontstaat een meer sigmoïde stap.

5. Betekenis en Conclusie

Dit werk is significant omdat het een fundamentele stap zet in het realistisch modelleren van medische ingrepen waarbij thermische effecten een rol spelen.

• Thermodynamische consistentie: Het model is strikt afgeleid volgens thermodynamische principes, wat het betrouwbaar maakt voor toekomstige uitbreidingen.
• Toekomstperspectief: De auteurs benadrukken dat de huidige resultaten gebaseerd zijn op verwaarloosbare capillaire krachten en constante viscositeit. In toekomstig werk zullen ze het uithardingsproces van het cement (als warmtebron) integreren en de noodzaak van de LTNE-aanname verder onderzoeken, aangezien snellere injecties of exotherme reacties grotere temperatuurverschillen kunnen veroorzaken.
• Toepassing: Het model biedt een basis voor pre-operatieve planning en het voorspellen van complicaties (zoals cementlekkage) door rekening te houden met thermische en mechanische koppelingen.

Samenvattend biedt dit paper een robuust theoretisch raamwerk voor niet-isotherme stroming in poreuze media, specifiek gericht op het verbeteren van de veiligheid en effectiviteit van vertebroplastie-procedures.