Multiphase modeling of anisotropic biomass particle pyrolysis accounting for particle deformation and coupled gas-phase dynamics

Dit artikel presenteert een nieuw single-grid, Euler-VOF-model binnen het open-source Basilisk-framework dat de gekoppelde vaste stof-gasdynamica en anisotrope deformatie van deeltjes tijdens biomassa-pyrolyse volledig oplost, waarbij een uitstekende overeenstemming met experimentele gegevens wordt aangetoond terwijl een robuuste tool wordt geboden voor de ontwikkeling van duurzame pyrolyseprocessen.

De kern

Het Grote Plaatje: Een Boomstam Koken Zonder Recept

Stel je voor dat je een boomstam perfect probeert te koken in een vuur om er brandstof van te maken (een proces dat pyrolyse wordt genoemd). Om dit goed te doen, moet je precies weten wat er binnenin de boomstam gebeurt terwijl deze opwarmt.

Lange tijd hadden wetenschappers twee aparte manieren om naar dit probleem te kijken:

1. Het "Binnenste" Perspectief: Ze keken naar hoe het hout krimpt en van binnenuit verandert, maar gokten hoe de hete lucht aan de buitenkant het aanraakte.
2. Het "Buitenste" Perspectief: Ze keken naar hoe de hete lucht rond de boomstam bewoog, maar behandelden de boomstam als een statische rots die nooit van vorm verandert.

Het probleem is dat hout geen rots is. Terwijl het "kookt", krimpt het, wordt het sponsachtig (poreus), en de hete lucht die erlangs stroomt verandert omdat de vorm van de stam verandert. De oude methoden misten het gesprek tussen het "binnenste" en het "buitenste".

De Nieuwe Oplossing: Eén Slimme Camera

Dit artikel introduceert een nieuw computermodel dat werkt als een enkele, high-definition camera die de hele scène tegelijkertijd bekijkt. Het model raadt niet hoe de lucht en het hout met elkaar interageren; het berekent de exacte dans tussen hen.

Zo hebben de auteurs deze "camera" gebouwd:

1. De "Volume-of-Fluid"-truc (De Waterballon-analogie)

Computers hebben meestal moeite met het volgen van een bewegende grens, zoals een krimpendende ballon. Dit model gebruikt een methode genaamd Volume-of-Fluid (VOF).

• De Analogie: Stel je een raster van kleine vakjes voor die je scherm bedekken. Sommige vakjes zijn gevuld met "hout", sommige met "lucht", en sommige met een mengsel. Naarmate het hout krimpt, werkt het model simpelweg het percentage "hout" in elk vakje bij. Het volgt de rand van het hout terwijl deze beweegt, net zoals het volgen van de rand van een waterballon die wordt ingedrukt.

2. Het "Spons"-effect (Porositeit en Krimping)

Hout is als een spons. Wanneer het opwarmt, gebeuren er twee dingen tegelijkertijd:

• De spons krijgt gaatjes: Het materiaal binnenin breekt af, waardoor er meer lege ruimte ontstaat (porositeit).
• De spons wordt kleiner: De hele stam krimpt in omvang.

De auteurs creëerden een speciale regel (een wiskundige functie die ze Z noemen) om te bepalen hoeveel van de reactie ervoor zorgt dat het hout gaatjes krijgt versus hoeveel het ervoor zorgt dat het krimpt. Het is alsoك beslissen of een smeltend ijsblokje verandert in een plas water (gaatjes krijgen) of gewoon kleiner wordt (krimpen). Ze ontdekten dat de beste resultaten voortkomen uit een mix van beide.

3. De "Verkeersopstopping" (Gasstroom binnenin)

Terwijl het hout kookt, komen er gassen vrij. Deze gassen moeten door de kleine gaatjes in het hout naar buiten sijpelen.

• De Analogie: Stel je mensen voor die een overvolle stadion uit willen rennen. Als het stadion wijd open is, rennen ze snel. Als de uitgangen smal en druk zijn, bewegen ze langzaam. Het model gebruikt Darcy-Forchheimer-vergelijkingen om dit "verkeersopstopping"-effect te berekenen, zodat het gas niet zomaar maglijk buiten verschijnen, maar zich echt een weg door de poriën van het hout moet banen.

4. De "Houtnerf" (Anisotropie)

Hout is niet in elke richting hetzelfde. Warmte reist sneller langs de nerf (zoals hardlopen door een gang) dan dwars door de nerf (zoals door een menigte lopen).

• De Analogie: Denk aan een stapel papier. Het is makkelijk om met je vinger langs de stapel te glijden (snel), maar moeilijk om door de stapel heen te duwen (langzaam). Het model houdt hier rekening mee door warmte en gasstroom sneller te laten verlopen in de richting van de houtvezels en langzamer er dwars op.

Wat hebben ze getest?

Het team heeft hun model getest tegen echte experimenten met houtdeeltjes variërend van kleine bollen tot cilinders. Ze controleerden:

• Temperatuur: Voorspelt het model dat het hout met de juiste snelheid warm wordt? (Ja, het kwam goed overeen).
• Massaverlies: Voorspelt het model hoeveel hout verandert in gas versus houtskool? (Ja, binnen een zeer kleine foutmarge).
• Vormverandering: Laat het model zien dat het hout correct krimpt? (Ja, hoewel het voorspellen van de exacte uiteindelijke vorm nog een beetje lastig is, was de algemene trend correct).

De Kern van het Verhaal

Dit artikel presenteert een nieuw, verenigd hulpmiddel dat stopt met gokken hoe hout krimpt en hoe lucht eromheen beweegt. In plaats daarvan simuleert het het hele proces in één keer.

• Waarom het belangrijk is: Het helpt ingenieurs bij het ontwerpen van betere systemen om hout om te zetten in hernieuwbare energie.
• De Kanttekening: Het model is complex en vereist veel rekenkracht, maar de auteurs hebben hun code open-source gemaakt (gratis voor iedereen om te gebruiken en te verbeteren).

Kortom, ze hebben een digitale tweeling gebouwd van een brandend stuk hout die zowel het binnenste als het buitenste begrijpt, waardoor wetenschappers de "onzichtbare" veranderingen kunnen zien die binnenin het hout plaatsvinden terwijl het verandert in brandstof.

Technische samenvatting

Technische Samenvatting: Multifasische Modellering van Anisotrope Biomassa-partikelpyrolyse

Probleemstelling
Huidige numerieke modellen voor de pyrolyse van biomassa-deeltjes lijden doorgaans aan een fundamentele beperking: ze richten zich ofwel op de evolutie van het vaste deeltje (interne structuur, krimp, porositeit), ofwel op de omringende gasfase-dynamica, maar zelden op de gekoppelde interacties tussen beide. Bestaande benaderingen vertrouwen vaak op sub-grid-schaal correlaties om warmte- en massaoverdrachtscoëfficiënten aan de deeltelgrens te schatten, uitgaande van een stationair grensvlak. Deze aanname vereenvoudigt de modelcomplexiteit, maar faalt in het vastleggen van de significante impact van variaties in de deeltjesdiameter en de beweging van het grensvlak op de uitwisselingssnelheden, met name voor biomassa die aanzienlijke morfologische veranderingen ondergaat (krimp en porositeïtevolutie) tijdens de conversie. Bovendien wordt het anisotrope karakter van biomassa-eigenschappen, zoals thermische geleidbaarheid, vaak verwaarloosd in multidimensionale kaders die de externe stroming oplossen.

Methodologie
Om deze beperkingen aan te pakken, stellen de auteurs een verenigde, single-grid Euleriaanse model voor dat zowel het vaste biomassa-deeltje als de omringende gasfase volledig oplost zonder afhankelijk te zijn van sub-grid correlaties. De kernmethodologie omvat:

• Volume-of-Fluid (VOF) Benadering: Het twee-fasen systeem wordt gevolgd met behulp van een VOF-methode om het scherpe grensvlak tussen het externe gas en de biomassa op te lossen. De biomassa wordt behandeld als een "pseudo-fase" die een mengsel bevat van een vaste matrix en intern gas binnen de poriën.
• Gekoppelde Besturende Vergelijkingen: Het model lost behoudvergelijkingen op voor massa, impuls, energie en chemische soorten voor zowel de externe omgeving als de interne pseudo-fase.
  • Impuls: Darcy–Forchheimer-termen worden geïntegreerd in de Navier-Stokes-vergelijkingen om effecten van weerstand binnen de poreuze matrix te verklaren.
  • Energie: Een enkele energievergelijking wordt opgelost voor de pseudo-fase, waarbij wordt uitgegaan van thermisch evenwicht tussen het interne gas en de vaste stof. Anisotrope thermische geleidbaarheid wordt expliciet gemodelleerd.
  • Chemie: Reacties van de vaste fase worden gemodelleerd met behulp van een gedetailleerd kinetisch mechanisme (CRECK-mechanisme) dat wordt beheerd via de OpenSMOKE++ bibliotheek. Gasfase-reacties worden genegeerd voor de specifieke laagtemperatuurpyrolysecondities die zijn getest, maar worden erkend als noodzakelijk voor toekomstige verbrandings-/vergassingsstudies.
• Krimp- en Porositeitskoppeling: Een nieuwe benadering wordt geïntroduceerd om de koppeling tussen de evolutie van de biomassa-porositeit en de deeltelskrimp te realiseren. Het model introduceert een functie, $Z(x,t)$, die de bronterm van de chemische reactie verdeelt tussen volumevermindering (krimp) en interne porositeitsveranderingen. Deze functie wordt opgelost naast een potentiaalstroombenadering om de krimpende snelheid van de poreuze matrix te bepalen.
• Numerieke Implementatie: Het model is geïmplementeerd in het open-source Basilisk-framework met behulp van een eindige volume-methode met adaptive mesh refinement (quad/octree grid). Een tweede-orde benaderde projectiemethode lost de Navier-Stokes-vergelijkingen op, terwijl scalaire transport (temperatuur, soorten) wordt afgehandeld via een two-field benadering om numerieke diffusie over het grensvlak te minimaliseren.

Belangrijkste Bijdragen

1. Verenigd Multifasisch Kader: Het werk presenteert een model dat het gas-vast grensvlak dynamisch oplost, waardoor de noodzaak voor empirische warmte- en massaoverdrachtcorrelaties vervalt.
2. Anisotrope en Deformeerbare Deeltjesmodellering: Het kader houdt rekening met het anisotrope karakter van biomassa (richtingsafhankelijke thermische geleidbaarheid en permeabiliteit) en koppelt interne porositeitsveranderingen aan externe deeltelskrimp in één enkele multidimensionale simulatie.
3. Massabehoud en Convergentie: Het model demonstreert massabehoud en numerieke convergentie. Hoewel het bewegende grensvlak resulteert in een eerste-orde convergentie voor massabehoud (vanwege niet-divergentievrije stromingen), vertonen stationaire grensvlakken een tweede-orde convergentie.
4. Open-Source Beschikbaarheid: De code en simulatie-opstellingen zijn publiekelijk beschikbaar gemaakt binnen het Basilisk-framework om reproduceerbaarheid te bevorderen.

Resultaten en Validatie
Het model werd gevalideerd tegen uitgebreide experimentele gegevens betreffende houtdeeltjes op centimeterschaal bij bedrijfstemperaturen tussen 400°C en 700°C.

• Isotroop Bolvormig Deeltje (Opgelegde Temperatuur): Het model voorspelde nauwkeurig de temperatuurprofielen (kern, halve straal, oppervlak) en de vrijgavecijfers van vluchtige soorten (CO, CO2, CH4, etc.) voor een populierenhoutbol. Het legde succesvol de endotherme ontledingspieken en exotherme koolstofvormingspieken vast.
• Houten Bol in Stroming: Simulaties van Schima superba houtbollen in een laminaire hete stikstofstroom vertoonden een uitstekende overeenkomst met experimentele massaverlies- en kerntemperatuurprofielen. Het model legde het temperatuurplateau correct vast dat wordt veroorzaakt door vochtverdamping (ca. 373 K) door een reversibel reactiemechanisme voor waterdesorptie/resorptie te implementeren.
• Krimpprofielen: Het model voorspelde krimptrends met een gemiddelde fout van ongeveer 10% in de uiteindelijke deeltelvorm. Hoewel het model de algemene trend van vormevolutie vastlegde, benadrukte het dat een dieper fundamenteel begrip van de interne structurele evolutie nog vereist is om de uiteindelijke vormen perfect te voorspellen.
• Anisotrope Cilinder: Voor een beuken houtcilinder onder gedwongen convectie reproduceerde het model succesvol de "botachtige" deformatie waargenomen in experimenten, waarbij hoeken sneller degraderen dan het centrum door anisotrope warmteoverdracht en randvoorwaarden. De berekende warmteoverdrachtscoëfficiënten (convectie en straling) kwamen goed overeen met experimentele schattingen.

Betekenis en Claims
De auteurs stellen dat dit werk een stap voorwaarts is in het ondersteunen van de ontwikkeling van duurzame pyrolyseprocessen door een rigoureuzere fysieke beschrijving van deeltjesgedrag te bieden. Het vermogen van het model om het grensvlak en de interne dynamiek op te lossen zonder empirische correlaties maakt het een krachtig instrument voor:

• Multi-schaal Modellering: Het dienen als een sub-grid modelgenerator om specifieke correlaties af te leiden voor grootschalige reactor-simulaties (bijv. bubbling beds, auger reactoren) waar resolutie van het grensvlak computationeel onhaalbaar is.
• Procesoptimalisatie: Het voorspellen van opbrengsten, degradatietijden en temperatuurprofielen onder complexe fluïdumdynamische condities.
• Fundamenteel Inzicht: Het bieden van inzichten in de wisselwerking tussen porositeitsveranderingen en volumekrimp, die moeilijk experimenteel te meten zijn.

Het artikel sluit bescheiden af met de opmerking dat hoewel het model de belangrijkste experimentele trends vastlegt, de exacte functionele vorm van de krimp-porositeitskoppeling ($Z$) nog steeds een gebied is dat verdere fundamentele onderzoek vereist. Toekomstig werk wordt geïdentificeerd als noodzakelijk om het model uit te breiden naar verbranding en vergassing door gasfase-reacties op te nemen en grotere chemische mechanismen te hanteren.