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

Questo articolo presenta un nuovo modello Eulerian-VOF a griglia singola all'interno del framework open-source Basilisk che risolve completamente la dinamica accoppiata solido-gas e la deformazione anisotropa delle particelle durante la pirolisi della biomassa, dimostrando un eccellente accordo con i dati sperimentali e fornendo al contempo uno strumento robusto per lo sviluppo di processi di pirolisi sostenibili.

L'essenziale

Il quadro generale: Cucinare un tronco senza ricetta

Immaginate di cercare di cuocere perfettamente un tronco di legno in un fuoco per trasformarlo in combustibile (un processo chiamato pirolisi). Per farlo bene, è necessario sapere esattamente cosa accade all'interno del tronco mentre si scalda.

Per molto tempo, gli scienziati hanno avuto due modi separati per osservare questo problema:

1. La visione "interna": Osservavano come il legno si restringesse e cambiasse all'interno, ma ipotizzavano come l'aria calda esterna lo stesse toccando.
2. La visione "esterna": Osservavano come l'aria calda si muovesse attorno al tronco, ma trattavano il legno come una roccia statica che non cambiava mai forma.

Il problema è che il legno non è una roccia. Mentre "cuoce", si restringe, diventa spugnoso (poroso) e l'aria calda che scorre accanto ad esso cambia perché la forma del tronco sta cambiando. I vecchi metodi non coglievano la conversazione tra l' "interno" e l' "esterno".

La nuova soluzione: Una singola telecamera intelligente

Questo articolo presenta un nuovo modello informatico che agisce come una singola telecamera ad alta definizione che osserva l'intera scena contemporaneamente. Non ipotizza come l'aria e il legno interagiscano; calcola l'esatto ballo tra di loro.

Ecco come gli autori hanno costruito questa "telecamera":

1. Il trucco del "Volume di Fluido" (L'analogia del palloncino d'acqua)

Di solito, i computer faticano a tracciare un confine in movimento, come un palloncino che si restringe. Questo modello utilizza un metodo chiamato Volume di Fluido (VOF).

• L'analogia: Immaginate una griglia di minuscole scatole che copre il vostro schermo. Alcune scatole sono piene di "legno", altre di "aria" e alcune sono un mix delle due. Man mano che il legno si restringe, il modello aggiorna semplicemente la percentuale di "leganno" in ogni scatola. Traccia il bordo del legno mentre si muove, proprio come si traccia il bordo di un palloncino d'acqua che viene schiacciato.

2. L'effetto "Spugna" (Porosità e contrazione)

Il legno è come una spugna. Quando si scalda, accadono due cose simultaneamente:

• La spugna si riempie di buchi: Il materiale all'interno si decompone, creando più spazio vuoto (porosità).
• La spugna diventa più piccola: L'intero tronco si restringe nelle dimensioni.

Gli autori hanno creato una regola speciale (una funzione matematica che chiamano Z) per decidere quanta parte della reazione causi la formazione di buchi nel legno rispetto a quanta ne causi il restringimento. È come decidere se un cubetto di ghiaccio che si scioglie si sta trasformando in una pozzanghera (creando buchi) o se sta solo diventando più piccolo (restringendosi). Hanno scoperto che i risultati migliori derivano da un mix di entrambi.

3. Il "Ingorgo stradale" (Flusso di gas all'interno)

Mentre il legno cuoce, rilascia gas. Questi gas devono incastrarsi attraverso i piccoli fori all'interno del legno per uscire.

• L'analogia: Immaginate persone che cercano di uscire da uno stadio affollato. Se lo stadio è spalancato, corrono veloci. Se le uscite sono strette e affollate, si muovono lentamente. Il modello utilizza le equazioni di Darcy-Forchheimer per calcolare questo effetto di "ingorgo", assicurando che il gas non appaia magicamente all'esterno, ma che si faccia strada attraverso i pori del legno.

4. La "Venatura del legno" (Anisotropia)

Il legno non è uguale in tutte le direzioni. Il calore viaggia più velocemente lungo la venatura (come correre lungo un corridoio) rispetto a trasversalmente ad essa (come attraversare una folla).

• L'analogia: Pensate a una pila di fogli di carta. È facile far scorrere un dito lungo la pila (veloce), ma difficile spingere attraverso la pila (lento). Il modello tiene conto di questo rendendo il flusso di calore e di gas più veloce nella direzione delle fibre del legno e più lento attraverso di esse.

Cosa hanno testato?

Il team ha testato il loro modello confrontandolo con esperimenti reali con particelle di legno che vanno da piccole sfere a cilindri. Hanno controllato:

• Temperatura: Il modello prevede che il legno si scaldi alla velocità corretta? (Sì, ha corrisposto bene).
• Perdita di massa: Il modello prevede quanto legno si trasforma in gas rispetto al carbone? (Sì, entro un margine di errore molto piccolo).
• Cambio di forma: Il modello mostra correttamente il restringimento del legno? (Sì, anche se prevedere l'esatta forma finale è ancora un po' complicato, l'andamento generale era corretto).

In sintamente

Questo articolo presenta uno strumento unificato che smette di ipotizzare come il legno si restringa e come l'aria si muova intorno ad esso. Invece, simula l'intero processo in un unico passaggio.

• Perché è importante: Aiuta gli ingegneri a progettare sistemi migliori per trasformare il legno in energia rinnovabile.
• Il limite: Il modello è complesso e richiede molta potenza di calcolo, ma gli autori hanno reso il loro codice open-source (gratuito per chiunque voglia usarlo e migliorarlo).

In breve, hanno costruito un gemello digitale di un pezzo di legno che brucia, che comprende sia l'interno che l'esterno, permettendo agli scienziati di vedere i cambiamenti "invisibili" che avvengono all'interno del legno mentre si trasforma in combustibile.

Sommario tecnico

Sintesi Tecnica: Modellazione Multifase della Pirolisi di Particelle di Biomassa Anisotrope

Definizione del Problema
Gli attuali modelli numerici per la pirolisi delle particelle di biomassa soffrono tipicamente di un limite fondamentale: si concentrano o sull'evoluzione della particella solida (struttura interna, contrazione, porosità) o sulla dinamica della fase gassosa circostante, ma raramente sulle interazioni accoppiate tra le due. Gli approcci esistenti si affidano spesso a correlazioni a scala sub-griglia per stimare i coefficienti di scambio termico e di massa all'interfaccia della particella, assumendo un'interfaccia stazionaria. Questa ipotesi semplifica la complessità del modello ma non riesce a catturare l'impatto significativo della variazione del diametro della particella e del movimento dell'interfaccia sui tassi di scambio, particolarmente per la biomassa che subisce cambiamenti morfologici sostanziali (contrazione e evoluzione della porosità) durante la conversione. Inoltre, la natura anisotropa delle proprietà della biomassa, come la conducibilità termica, è spesso trascurata nei framework multidimensionali che risolvono il flusso esterno.

Metodologia
Per affrontare queste limitazioni, gli autori propongono un modello Euleriano unificato a griglia singola che risolve completamente sia la particella di biomassa solida che la fase gassosa circostante senza fare affidamento su correlazioni sub-griglia. La metodologia principale include:

• Approccio Volume-of-Fluid (VOF): Il sistema bifase è tracciato utilizzando un metodo VOF per risolvere l'interfaccia netta tra il gas esterno e la biomassa. La biomassa è trattata come una "pseudo-fase" contenente una miscela di una matrice solida e gas interno all'interno dei pori.
• Equazioni Governative Accoppiate: Il modello risolve le equazioni di conservazione per massa, quantità di moto, energia e specie chimiche sia per l'ambiente esterno che per la pseudo-fase interna.
  • Quantità di Moto: I termini di Darcy–Forchheimer sono incorporati nelle equazioni di Navier-Stokes per rendere conto degli effetti di trascinamento all'interno della matrice porosa.
  • Energia: Viene risolta un'unica equazione dell'energia per la pseudo-fase assumendo l'equilibrio termico tra il gas interno e il solido. La conducibilità termica anisotropa è modellata esplicitamente.
  • Chimica: Le reazioni di pirolisi della fase solida sono modellate utilizzando un meccanismo cinetico dettagliato (meccanismo CRECK) gestito tramite la libreria OpenSMOKE++. Le reazioni in fase gassosa sono trascurate per le specifiche condizioni di pirolisi a bassa temperatura testate, ma sono riconosciute come necessarie per futuri studi di combustione/gassificazione.
• Accoppiamento tra Contrazione e Porosità: Viene introdotto un approccio innovativo per accoppiare l'evoluzione della porosità della biomassa con la contrazione della particella. Il modello introduce una funzione, $Z(x,t)$, che ripartisce il termine sorgente della reazione chimica tra la riduzione di volume (contrazione) e i cambiamenti della porosità interna. Questa funzione viene risolta insieme a un'approssimazione di flusso potenziale per determinare la velocità di contrazione della matrice porosa.
• Implementazione Numerica: Il modello è implementato nel framework open-source Basilisk utilizzando un metodo a volumi finiti con raffinamento adattivo della mesh (griglia quad/octree). Un metodo di proiezione approssimato del secondo ordine risolve le equazioni di Navier-Stokes, mentre il trasporto scalare (temperatura, specie) è gestito tramite un approccio a due campi per minimizzare la diffusione numerica attraverso l'interfaccia.

Contributi Chiave

1. Framework Multifase Unificato: Il lavoro presenta un modello che risolve dinamicamente l'interfaccia gas-solido, eliminando la necessità di correlazioni empiriche per lo scambio termico e di massa.
2. Modellazione di Particelle Deformabili e Anisotrope: Il framework tiene conto della natura anisotropa della biomassa (conducibilità termica e permeabilità dipendenti dalla direzione) e accoppia i cambiamenti di porosità interna con la contrazione esterna della particella in una singola simulazione multidimensionale.
3. Conservazione della Massa e Convergenza: Il modello dimostra la conservazione della massa e la convergenza numerica. Sebbene l'interfaccia mobile comporti una convergenza del primo ordine per la conservazione della massa (a causa di flussi non privi di divergenza), le interfacce stazionarie mostrano una convergenza del secondo ordine.
4. Disponibilità Open-Source: Il codice e le configurazioni di simulazione sono resi pubblicamente disponibili all'interno del framework Basilisk per migliorare la riproducibilità.

Risultati e Validazione
Il modello è stato validato rispetto a estesi dati sperimentali riguardanti particelle di legno su scala centimetrica a temperature operative comprese tra 400°C e 700°C.

• Sfera Isotropa (Temperatura Imposta): Il modello ha previsto accuratamente i profili di temperatura (nucleo, raggio medio, superficie) e i tassi di rilascio delle specie volatili (CO, CO2, CH4, ecc.) per una sfera di legno di pioppo. Ha catturato con successo i picchi di decomposizione endotermica e i picchi di formazione di char esotermica.
• Sfera di Legno in Flusso: Le simulazioni di sfere di legno di Schima superba in un flusso laminare di azoto caldo hanno mostrato un eccellente accordo con i dati sperimentali di perdita di massa e profili di temperatura del nucleo. Il modello ha catturato correttamente il plateau di temperatura causato dall'evaporazione dell'umidità (circa 373 K) implementando un meccanismo di reazione reversibile per la desorbimento/riassorbimento dell'acqua.
• Profili di Contrazione: Il modello ha previsto le tendenze di contrazione con un errore medio di circa il 10% nella forma finale della particella. Sebbene il modello abbia catturato la tendenza generale dell'evoluzione della forma, ha evidenziato che è necessaria una comprensione fondamentale più profonda dell'evoluzione strutturale interna per prevedere perfettamente le forme finali.
• Cilindro Anisotropo: Per un cilindro di legno di faggio in convezione forzata, il modello ha riprodotto con successo la deformazione "simile a un osso" osservata negli esperimenti, dove gli angoli degradano più velocemente del centro a causa del trasferimento di calore anisotropo e delle condizioni al contorno. I coefficienti di scambio termico calcolati (convezione e radiazione) erano in buon accordo con le stime sperimentali.

Significato e Rivendicazioni
Gli autori affermano che questo lavoro rappresenta un passo avanti nell'aiutare lo sviluppo di processi di pirolisi sostenibili, fornendo una descrizione fisica più rigorosa del comportamento delle particelle. La capacità del modello di risolvere l'interfaccia e la dinamica interna senza correlazioni empiriche lo rende uno strumento potente per:

• Modellazione Multi-scala: Servire come generatore di modelli sub-griglia per derivare specifiche correlazioni per simulazioni di reattori su scala maggiore (ad esempio, letti fluidizzati, reattori ad elica), dove la risoluzione dell'interfaccia è computazionalmente proibitiva.
• Ottimizzazione dei Processi: Prevedere rese, tempi di degradazione e profili di temperatura sotto complesse condizioni di fluidodinamica.
• Intuizione Fondamentale: Fornire approfondimenti sull'interazione tra cambiamenti di porosità e contrazione di volume, che sono difficili da misurare sperimentalmente.

Il documento conclude con modestia, osservando che sebbene il modello catturi le principali tendenze sperimentali, la forma funzale precisa dell'accoppiamento contrazione-porosità ($Z$) rimane un'area che richiede un'ulteriore investigazione fondamentale. Si identifica come lavoro futuro la necessità di estendere il modello alle condizioni di combustione e gassificazione, incorporando reazioni in fase gassosa e gestendo meccanismi chimici più ampi.