Local Thermal Non-Equilibrium Models in Porous Media: A Comparative Study of Conduction Effects

Deze studie vergelijkt lokale thermische niet-evenwichtsmodellen (LTNE) met een porie-opgeloste referentie voor puur geleidende poreuze media en toont aan dat REV-schaalmodellen die gebruikmaken van op homogenisatie gebaseerde effectieve parameters interfaciële weerstandseffecten nauwkeurig vastleggen, terwijl dual-netwerkmodellen met vaste ruimtelijke resolutie grotere afwijkingen vertonen.

De kern

Stel je een poreus materiaal voor, zoals een spons of een rots, als een drukke stad. Deze stad heeft twee soorten bewoners: de vaste korrels (de gebouwen) en het vloeistof (het water of de lucht die door de straten stroomt).

Wanneer je één kant van deze stad opwarmt, wil je weten hoe de temperatuur zich verspreidt. De grote vraag die dit artikel stelt is: Verwarmen de gebouwen en het water erin zich met exact dezelfde snelheid, of lopen ze achter op elkaar?

Hier volgt een uiteenzetting van wat de onderzoekers deden, met behulp van eenvoudige analogieën.

1. De Twee Manieren om over Warmte na te Denken

In het verkeer gingen wetenschappers meestal uit van Lokale Thermische Evenwicht (LTE).

• De Analogie: Stel je een kamer vol mensen voor die hand in hand houden. Als één persoon heet wordt, voelt iedereen anders dit direct. In dit model zijn de "gebouwen" en het "water" zo perfect verbonden dat ze op elk gegeven moment exact dezelfde temperatuur hebben op dezelfde plek. Het is alsof ze één enkel brein delen.

Echter, de onderzoekers wisten dat dit niet altijd waar is. Soms is de verbinding tussen het gebouw en het water "kleverig" of traag. Dit is Lokale Thermische Niet-Evenwicht (LTNE).

• De Analogie: Stel je voor dat de mensen in aparte kamers zitten met dikke, geïsoleerde deuren ertussen. Als je het water in de hal verwarmt, kunnen de gebouwen een tijdje koel blijven omdat de warmte moeite moet doen om door de deur te komen. Het water wordt heet, maar het gebouw blijft even koud. Ze hebben verschillende temperaturen op dezelfde locatie.

2. De Drie "Kaarten" die Gebruikt Worden om Warmte te Voorspellen

Om uit te zoeken wanneer deze "vertraging" optreedt en hoe je deze kunt voorspellen, vergeleek het team drie verschillende manieren om een kaart van deze stad te tekenen:

• Kaart A: Het "Straatniveau"-beeld (Pore-Resolved Model)

  • Wat het is: Dit is de meest gedetailleerde kaart. Het tekent elk enkel gebouw en elke enkele straat. Het ziet de exacte vorm van de rots en het water.
  • Het nadeel: Het is ongelooflijk traag en computergewijs duur, alsof je probeert elk zandkorreltje op een strand te simuleren. De onderzoekers gebruikten dit als hun "Gouden Standaard" of referentie om te zien of de andere kaarten juist waren.

• Kaart B: Het "Buurt"-beeld (Dual-Network Model)

  • Wat het is: In plaats van elke straat te tekenen, vereenvoudigt deze kaart de stad tot een netwerk van punten (die de gebouwen en waterzakken vertegenwoordigen) verbonden door lijnen (die de verbindingen ertussen vertegenwoordigen).
  • Het nadeel: Het is sneller, maar het heeft een vaste resolutie. Het is alsof je naar een stad kijkt door een raamrooster; je kunt niet dichter inzoomen dan de grootte van het raam. Het artikel vond dat omdat dit rooster vaststaat, het soms de scherpe temperatuurveranderingen die direct aan de randen plaatsvinden, mist.

• Kaart C: Het "Luchtfoto"-beeld (REV-Schaal Model)

  • Wat het is: Dit is een hoog niveau, een gemiddelde kaart. Het ziet geen individuele gebouwen; het ziet "blokken" van de stad. Het gebruikt wiskunde om het gemiddelde gedrag van het hele blok te raden.
  • Het nadeel: Om dit werk te laten doen, moet je de "gemiddelde eigenschappen" van het blok raden. Als je verkeerd raadt, is de hele kaart verkeerd.

3. Het Grote Experiment

De onderzoekers voerden simulaties op een computer uit om te zien hoe warmte zich door deze "stad" bewoog onder twee verschillende omstandigheden:

Scenario 1: De Open Deur (Lage Weerstand)

• De Opzet: De verbinding tussen het water en de rots was perfect (zoals een wijd open deur). Warmte stroomde vrij.
• Het Resultaat: De "Open Deur" betekende dat het water en de rots direct samen opwarmden. De LTE-aanname (het enkele brein) werkte perfect. Alle drie de kaarten gaven bijna hetzelfde antwoord. De "vertraging" bestond niet.

Scenario 2: De Geïsoleerde Deur (Hoge Weerstand)

• De Opzet: De verbinding was geblokkeerd of "kleverig" (zoals een dikke, geïsoleerde deur). Warmte had het moeilijk om van het water naar de rots te springen.
• Het Resultaat: Nu werd het water heet, maar bleef de rots een tijdje koel. De LTE-aanname faalde volledig.
  • De Straatniveau-kaart toonde de exacte vertraging.
  • De Luchtfoto-kaart (als deze correct werd berekend met een specifieke wiskundige methode genaamd homogenisatie) kwam zeer goed overeen met de Straatniveau-kaart.
  • De Buurt-kaart was oké, maar omdat zijn "ramen" een vaste grootte hadden, gladde het de scherpe verschillen iets te veel af.

4. De Belangrijkste Conclusie

De belangrijkste bevinding gaat over hoe je de "Luchtfoto"-kaart berekent.

• Sommige oude manieren om de gemiddelde eigenschappen voor de Luchtfoto-kaart te berekenen, negeerden de "kleverige deur". Ze gingen ervan uit dat warmteoverdracht altijd perfect was. Toen de onderzoekers deze oude formules gebruikten, faalde de Luchtfoto-kaart om de vertraging tussen het water en de rots te tonen.
• Echter, toen ze een specifieke, geavanceerdere wiskundige methode (homogenisatie) gebruikten die wel rekening hield met de "kleverige deur" (de grensvlakweerstand), werd de Luchtfoto-kaart ongelooflijk nauwkeurig. Het kwam bijna perfect overeen met het gedetailleerde Straatniveau-beeld, hoewel het veel eenvoudiger was.

Samenvatting

• Als de verbinding perfect is: Je kunt eenvoudige modellen gebruiken; alles warmt samen op.
• Als de verbinding traag/kleverig is: Je moet modellen gebruiken die toestaan dat water en rots verschillende temperaturen hebben.
• De Beste Afkorting: Als je een enorm systeem moet modelleren (zoals een hele aquifer of een brandstofcel) en niet elke korrel kunt simuleren, gebruik dan het "Luchtfoto"-model, maar zorg ervoor dat je de specifieke wiskunde gebruikt die rekening houdt met de weerstand tussen de materialen. Als je dat doet, zal je eenvoudige model net zo nauwkeurig zijn als het supergedetailleerde model.

Opmerking: Het artikel stelt expliciet dat deze studie alleen keek naar warmte die door stilstaande materialen beweegt (geleiding). Ze keken niet naar warmte die met stromend water beweegt (convectie), wat ze zeggen in een toekomstig artikel te zullen bestuderen.

Technische samenvatting

1. Probleemstelling

Het modelleren van warmteoverdracht in poreuze media vertrouwt vaak op de aanname van Lokale Thermische Evenwicht (LTE), die stelt dat de vloeistof- en vaste fasen binnen een controlevolume onmiddellijk dezelfde temperatuur delen. Hoewel dit voor veel scenario's geldig is, faalt deze aanname in systemen die gekenmerkt worden door:

• Grote temperatuurgradiënten.
• Beduidende verschillen in thermische eigenschappen tussen de fasen.
• Hoge thermische weerstanden aan het interface (lage warmteoverdrachtscoëfficiënten).
• Snelle dynamiek.

Wanneer LTE faalt, zijn Lokale Thermische Niet-Evenwicht (LTNE)-modellen vereist om onderscheiden temperaturen voor de vloeistof ($T_f$) en vaste ($T_s$) fasen en de warmte-uitwisseling tussen hen op te lossen. Er is echter een gebrek aan uitgebreide vergelijkende studies die continuum-schaal LTNE-modellen (Dual-Network en REV-schaal) valideren tegen referenties met hoge precisie die de poriën oplossen, met name wat betreft de nauwkeurigheid van formuleringen voor effectieve parameters (bijv. homogenisatie versus empirische correlaties).

2. Methodologie

De studie vergelijkt drie benaderingen op continuum-schaal voor een volledig verzadigd systeem (één vloeistof, één vast stof) waarbij pure warmtegeleiding wordt beschouwd (convectie en massatransport worden verwaarloosd).

A. De Drie Modelklassen

1. Poriën-opgelost Model (Referentie):

   • Benadering: Lost energiebalansvergelijkingen afzonderlijk op voor de vloeistof ($\Omega_f$) en het vaste ($\Omega_s$) domein op een rooster dat de exacte geometrie van de poriën-schaal oplost.
   • Interface: Lost expliciet de scherpe vloeistof-vaste interface ($\Gamma_{fs}$) op.
   • LTNE-implementatie: Gebruikt een eindige warmteoverdrachtscoëfficiënt ($h$) om een temperatuursprong over de interface toe te staan ($\lambda_f \nabla T_f \cdot n = h(T_s - T_f)$).
   • Oplosser: Finite Element Method (FEM) met behulp van Netgen/NGSolve.

2. Dual-Network Model:

   • Benadering: Stelt het poreuze medium voor als twee onderling verbonden netwerken (vloeistof en vast) bestaande uit geïdealiseerde lichamen (bollen) en keelstukken (cilinders).
   • Interface: Lost de scherpe interface niet op; in plaats daarvan gebruikt het interfaciale keelstukken met thermische transmissibiliteit gebaseerd op vormfactoren en effectieve oppervlakken.
   • Oplosser: Control Volume Finite Element (CVFE) methode met behulp van DuMux.

3. REV-Schaal Model (Representative Elementary Volume):

   • Benadering: Gebruikt gemiddelde balansvergelijkingen over een controlevolume waar vloeistof- en vastestoftemperaturen overlappen.
   • LTNE-implementatie: Lost twee gekoppelde energie-vergelijkingen op met een term voor interfaciale warmte-uitwisseling ($q_{cond} = \lambda_{eff,I} a_{fs} (T_s - T_f)$).
   • Effectieve Parameters: De studie evalueert drie verschillende formuleringen voor effectieve thermische geleidbaarheden en interfaciale termen:
     • Nuske et al.: Eenvoudige schaling op basis van volumefractie.
     • Nakayama et al.: Omvat tortuositeit en verhoudingen van geleidbaarheid.
     • Homogenisatietheorie: Leidt effectieve parameters wiskundig af uit vergelijkingen op poriën-schaal, waarbij expliciet de interfaciale warmteoverdrachtscoëfficiënt ($h$) wordt opgenomen.
   • Oplosser: Finite Volume Method (FV) met behulp van PorePy.

B. Simulatieopstelling

• Geometrie: Een periodiek 3D-domein bestaande uit kubische referentiecellen die een bol bevatten die verbonden is met zes buizen.
• Materialen: Graniet (vast) en Water (vloeistof).
• Randvoorwaarden: Kwaasi-1D opstelling met een temperatuurstap aan de linkerrand (vloeistof verwarmd met 10K, vast stof geïsoleerd).
• Testgevallen:
  • Geval 1 (Lage Weerstand): Hoge warmteoverdrachtscoëfficiënt ($h \approx 8.3 \times 10^7$ W/m²K), die standaard water-graniet interfaces nabootst.
  • Geval 2 (Hoge Weerstand): Lage warmteoverdrachtscoëfficiënt ($h = 100$ W/m²K), die interfaces met hoge weerstand nabootst (bijv. koper-helium) om sterke LTNE-effecten te induceren.

3. Belangrijkste Bijdragen

• Systematische Vergelijking: Biedt een rigoureuze benchmark van Dual-Network- en REV-schaalmodellen tegen een poriën-opgeloste referentie voor LTNE-omstandigheden.
• Analyse van Parameterformuleringen: Evalueert de doeltreffendheid van verschillende formuleringen voor effectieve parameters. Het toont aan dat homogenisatie-gebaseerde benaderingen superieur zijn voor het vastleggen van LTNE-gedrag omdat ze de interfaciale warmteoverdrachtscoëfficiënt wiskundief incorporeren, terwijl empirische formuleringen (Nuske, Nakayama) vaak falen in het voorspellen van LTNE in puur geleidende systemen zonder expliciete $h$-afhankelijkheid.
• Invloed van Resolutie: Benadrukt de kritieke rol van ruimtelijke resolutie. Het Dual-Network-model vertoonde afwijkingen van de referentie vanwege zijn vaste ruimtelijke resolutie, die niet kan worden verfijnd om steile temperatuurgradiënten nabij grenzen vast te leggen, in tegenstelling tot het REV-schaalmodel waar roosterverfijning mogelijk is.

4. Belangrijkste Resultaten

• Lage Interfaciale Weerstand (Hoge $h$):

  • Het systeem gedraagt zich als LTE; vloeistof- en vastestoftemperaturen zijn nagenoeg identiek in alle modellen.
  • Het REV-schaalmodel (met gebruik van gehomogeniseerde parameters) komt nauw overeen met de poriën-opgeloste referentie.
  • Het Dual-Network-model vertoonde iets lagere temperaturen en minder steile gradiënten in vergelijking met de referentie, toegeschreven aan zijn vaste ruimtelijke resolutie en numerieke diffusie.

• Hoge Interfaciale Weerstand (Lage $h$):

  • LTNE-effecten: Beduidende temperatuurverschillen tussen vloeistof- en vastestoffasen werden waargenomen, met name nabij de verwarmde grens.
  • Modelprestaties:
    • Het REV-schaalmodel met gehomogeniseerde parameters legde het LTNE-gedrag en de temperatuursprong nauwkeurig vast en kwam dicht in de buurt van de poriën-opgeloste referentie (hoewel iets sneller vanwege numerieke diffusie).
    • Het Dual-Network-model vertoonde opnieuw een minder steile gradiënt en lagere temperaturen vanwege het onvermogen om het rooster bij de interface te verfijnen.
    • LTE vs. LTNE: Modellen die interfaciale weerstand negeren (LTE) faalden in het vastleggen van de temperatuursprong, wat de noodzaak van LTNE-modellen voor scenario's met hoge weerstand bevestigt.

• Formuleringen voor Effectieve Parameters:

  • In het geval van hoge weerstand faalden formuleringen van Nuske en Nakayama (die $h$ inherent niet bevatten in de geleidende limiet) om LTNE-effecten te produceren en voorspelden ze in plaats daarvan LTE.
  • Alleen de homogenisatie-gebaseerde benadering slaagde erin het LTNE-gedrag vast te leggen omdat de opgeschaalde interfaciale term ($H = h \cdot a_{fs}$) expliciet afhankelijk is van de warmteoverdrachtscoëfficiënt.

5. Betekenis en Conclusie

• Validatie van Homogenisatie: De studie bevestigt dat homogenisatietheorie de meest betrouwbare effectieve parameters biedt voor REV-schaal LTNE-modellen, omdat het de fysica op poriën-schaal (inclusief interfaciale weerstand) rigoureus overbrugt naar de continuum-schaal.
• Beperkingen van Dual-Network-modellen: Hoewel nuttig voor inzichten op poriën-schaal, worden Dual-Network-modellen beperkt door hun vaste ruimtelijke resolutie. Ze worstelen om scherpe thermische gradiënten nabij grenzen op te lossen in vergelijking met continuummodellen die roosterverfijning toestaan.
• Praktische Implicaties: Voor toepassingen met hoge interfaciale weerstand (bijv. specifieke geothermische of cryogene systemen) kan het vertrouwen op empirische effectieve parameters zonder rekening te houden met de warmteoverdrachtscoëfficiënt leiden tot aanzienlijke fouten (LTE voorspellen terwijl LTNE bestaat).
• Toekomstig Werk: De auteurs merken op dat convectie een dominante factor is in de meeste praktische toepassingen en plannen deze vergelijkende studie uit te breiden om massatransport en convectieve warmteoverdracht op te nemen in een vervolgartikel.

Kortom, het artikel stelt vast dat REV-schaalmodellen die gebruikmaken van gehomogeniseerde effectieve parameters de meest robuuste continuumbenadering zijn voor het simuleren van LTNE in poreuze media, mits de interfaciale warmteoverdrachtscoëfficiënt expliciet wordt meegenomen in de modelformulering.