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

Cet article présente un nouveau modèle VOF eulérien à grille unique au sein du cadre open-source Basilisk qui résout pleinement la dynamique couplée solide-gaz et la déformation anisotrope des particules lors de la pyrolyse de la biomasse, démontrant un excellent accord avec les données expérimentales tout en fournissant un outil robuste pour le développement de processus de pyrolyse durables.

L'essentiel

La vue d'ensemble : Cuire une bûche sans recette

Imaginez que vous essayiez de cuire parfaitement une bûche de bois dans un feu pour la transformer en combustible (un processus appelé pyrolyse). Pour y parvenir, vous devez savoir exactement ce qui se passe à l'intérieur de la bûche pendant qu'elle chauffe.

Pendant longtemps, les scientifiques avaient deux manières distinctes d'aborder ce problème :

1. La vue de « l'intérieur » : Ils observaient comment le bois rétrécissait et changeait à l'intérieur, mais ils devinaient comment l'air chaud extérieur le touchait.
2. La vue de « l'extérieur » : Ils observaient comment l'air chaud se déplaçait autour de la bûche, mais ils traitaient la bûche comme un rocher statique qui ne change jamais de forme.

Le problème est que le bois n'est pas un rocher. À mesure qu'il cuit, il rétrécit, devient spongieux (poreux) et l'air chaud qui passe à côté de lui change car la forme de la bûche évolue. Les anciennes méthodes manquaient la conversation entre « l'intérieur » et « l'extérieur ».

La nouvelle solution : Une caméra intelligente unique

Ce document présente un nouveau modèle informatique qui agit comme une caméra haute définition unique observant toute la scène à la fois. Il ne devine pas comment l'air et le bois interagissent ; il calcule la danse exacte entre eux.

Voici comment les auteurs ont construit cette « caméra » :

1. L'astuce du « Volume de Fluide » (L'analogie du ballon d'eau)

D'ordinaire, les ordinateurs ont du mal à suivre une limite mobile, comme un ballon qui rétrécit. Ce modèle utilise une méthode appelée Volume de Fluide (VOF).

• L'analogie : Imaginez une grille de petites boîtes couvrant votre écran. Certaines boîtes sont remplies de « bois », d'autres d'« air », et certaines d'un mélange des deux. À mesure que le bois rétrécit, le modèle met simplement à jour le pourcentage de « bois » dans chaque boîte. Il suit le bord du bois au fur et à mesure qu'il se déplace, tout comme on suivrait le bord d'un ballon d'eau qui est pressé.

2. L'effet « Éponge » (Porosité et rétrécissement)

Le bois est comme une éponge. Lorsqu'il chauffe, deux choses se produisent simultanément :

• L'éponge se creuse : Le matériau à l'intérieur se décompose, créant plus d'espaces vides (porosité).
• L'éponge rétrécit : La taille globale de la bûche diminue.

Les auteurs ont créé une règle spéciale (une fonction mathématique qu'ils appellent Z) pour décider quelle part de la réaction provoque la création de trous dans le bois par rapport à la part qui le fait rétrécir. C'est comme décider si un glaçon qui fond se transforme en une flaque (création de trous) ou s'il devient simplement plus petit (rétrécissement). Ils ont découvert que les meilleurs résultats proviennent d'un mélange des deux.

3. Le « Embouteillage » (Flux de gaz à l'intérieur)

À mesure que le bois cuit, il libère des gaz. Ces gaz doivent se frayer un chemin à travers les minuscules trous à l'intérieur du bois pour sortir.

• L'analogie : Imaginez des gens essayant de sortir d'un stade bondé. Si le stade est grand ouvert, ils courent vite. Si les sorties sont étroites et encombrées, ils avancent lentement. Le modèle utilise les équations de Darcy-Forchheimer pour calculer cet effet d'« embouteillage », garantissant que le gaz n'apparaît pas par magie à l'extérieur, mais qu'il pousse réellement à travers les pores du bois.

4. Le « Grain du bois » (Anisotropie)

Le bois n'est pas le même dans toutes les directions. La chaleur se propage plus rapidement le long du grain (comme courir dans un couloir) que perpendiculairement à celui-ci (comme traverser une foule).

• L'analogie : Pensez à une pile de feuilles de papier. Il est facile de faire glisser un doigt le long de la pile (rapide), mais difficile de pousser à travers la pile (lent). Le modèle tient compte de cela en faisant circuler la chaleur et les gaz plus rapidement dans la direction des fibres du bois et plus lentement à travers elles.

Qu'ont-ils testé ?

L'équipe a testé son modèle par rapport à des expériences réelles avec des particules de bois allant de petites sphères à des cylindres. Ils ont vérifié :

• La température : Le modèle prédit-il que le bois chauffe à la bonne vitesse ? (Oui, la correspondance était bonne).
• La perte de masse : Le modèle prédit-il quelle quantité de bois se transforme en gaz par rapport au charbon ? (Oui, avec une marge d'erreur très faible).
• Le changement de forme : Le modèle montre-t-il que le bois rétrécit correctement ? (Oui, bien que prédire la forme finale exacte soit encore un peu complexe, la tendance générale était correcte).

L'essentiel

Ce document présente un nouvel outil unifié qui cesse de deviner comment le bois rétrécit et comment l'air circule autour de lui. Au lieu de cela, il simule l'ensemble du processus en une seule étape.

• Pourquoi c'est important : Cela aide les ingénieurs à concevoir de meilleurs systèmes pour transformer le bois en énergie renouvelable.
• Le bémol : Le modèle est complexe et nécessite une grande puissance de calcul, mais les auteurs ont rendu leur code open-source (libre d'accès pour que quiconque puisse l'utiliser et l'améliorer).

En résumé, ils ont construit un jumeau numérique d'un morceau de bois en train de brûler qui comprend à la fois l'intérieur et l'extérieur, permettant aux scientifiques de voir les changements « invisibles » qui se produisent à l'intérieur du bois lorsqu'il se transforme en combustible.

Résumé technique

Résumé Technique : Modélisation multiphasique de la pyrolyse de particules de biomasse anisotropes

Énoncé du problème
Les modèles numériques actuels de la pyrolyse des particules de biomasse souffrent d'une limitation fondamentale : ils se concentrent soit sur l'évolution de la particule solide (structure interne, rétrécissement, porosité), soit sur la dynamique de la phase gazeuse environnante, mais rarement sur les interactions couplées entre les deux. Les approches existantes reposent souvent sur des corrélations à l'échelle sous-maille pour estimer les coefficients de transfert de chaleur et de masse à la limite de la particule, en supposant une interface stationnaire. Cette hypothèse simplifie la complexité du modèle mais ne parvient pas à capturer l'impact significatif de la variation du diamètre de la particule et du mouvement de l'interface sur les taux d'échange, particulièrement pour la biomasse qui subit des changements morphologiques substantiels (rétrécissement et évolution de la porosité) lors de la conversion. De plus, la nature anisotrope des propriétés de la biomasse, telle que la conductivité thermique, est souvent négligée dans les cadres multidimensionnels qui résolvent l'écoulement externe.

Méthodologie
Pour remédier à ces limitations, les auteurs proposent un modèle eulérien unifié à grille unique qui résout pleinement la particule de biomasse solide et la phase gazeuse environnante sans recourir à des corrélations sous-maille. La méthodologie centrale comprend :

• Approche Volume de Fluide (VOF) : Le système biphasique est suivi à l'aide d'une méthode VOF pour résoudre l'interface nette entre le gaz externe et la biomasse. La biomasse est traitée comme une « pseudo-phase » contenant un mélange de matrice solide et de gaz interne au sein des pores.
• Équations de Conservation Couplées : Le modèle résout les équations de conservation de la masse, de la quantité de mouvement, de l'énergie et des espèces chimiques pour l'environnement externe et pour la pseudo-phase interne.
  • Quantité de mouvement : Des termes de Darcy–Forchheimer sont incorporés dans les équations de Navier-Stokes pour rendre compte des effets de traînée au sein de la matrice poreuse.
  • Énergie : Une seule équation d'énergie est résolue pour la pseudo-phase en supposant un équilibre thermique entre le gaz interne et le solide. La conductivité thermique anisotrope est explicitement modélisée.
  • Chimie : Les réactions de pyrolyse de la phase solide sont modélisées à l'aide d'un mécanisme cinétique détaillé (mécanisme CRECK) géré via la bibliothèque OpenSMOKE++. Les réactions en phase gazeuse sont négligées pour les conditions spécifiques de pyrolyse à basse température testées, mais sont reconnues comme nécessaires pour de futures études de combustion/gazéification.
• Couplage Rétrécissement et Porosité : Une approche novatrice est introduite pour coupler l'évolution de la porosité de la biomasse avec le rétrécissement de la particule. Le modèle introduit une fonction, $Z(x,t)$, qui répartit le terme source de la réaction chimique entre la réduction de volume (rétrécissement) et les changements de porosité interne. Cette fonction est résolue parallèlement à une approximation d'écoulement potentiel pour déterminer la vitesse de rétrécissement de la matrice poreuse.
• Implémentation Numérique : Le modèle est implémenté dans le cadre open-source Basilisk en utilisant une méthode de volumes finis avec un raffinement de maillage adaptatif (grille quad/octree). Une méthode de projection approchée de second ordre résout les équations de Navier-Stokes, tandis que le transport scalaire (température, espèces) est géré via une approche à deux champs pour minimiser la diffusion numérique à travers l'interface.

Principales Contributions

1. Cadre Multiphasique Unifié : Ce travail présente un modèle qui résout dynamiquement l'interface gaz-solide, éliminant ainsi le besoin de corrélations empiriques de transfert de chaleur et de masse.
2. Modélisation de Particules Anisotropes et Déformables : Le cadre rend compte de la nature anisotrope de la biomasse (conductivité thermique et perméabilité dépendantes de la direction) et couple les changements de porosité interne avec le rétrécissement externe de la particule dans une seule simulation multidimensionnelle.
3. Conservation de la Masse et Convergence : Le modèle démontre la conservation de la masse et la convergence numérique. Bien que l'interface mobile entraîne une convergence de premier ordre pour la conservation de la masse (en raison de flux non divergence-free), les interfaces stationnaires présentent une convergence de second ordre.
4. Disponibilité en Open-Source : Le code et les configurations de simulation sont rendus publics au sein du cadre Basilisk afin de favoriser la reproductibilité.

Résultats et Validation
Le modèle a été validé par rapport à des données expérimentales étendues impliquant des particules de bois à l'échelle du centimètre, à des températures de fonctionnement comprises entre 400 °C et 700 °C.

• Sphère Isotrope (Température imposée) : Le modèle a prédit avec précision les profils de température (cœur, demi-rayon, surface) et les taux de libération des espèces volatiles (CO, CO2, CH4, etc.) pour une sphère de peuplier. Il a capturé avec succès les pics de décomposition endothermique et les pics de formation de charbon exothermique.
• Sphère de Bois en Écoulement : Les simulations de sphères de bois de Schima superba dans un écoulement laminaire d'azote chaud ont montré un excellent accord avec les données expérimentales de perte de masse et les profils de température du cœur. Le modèle a correctement capturé le plateau de température causé par l'évaporation de l'humidité (environ 373 K) en implémentant un mécanisme de réaction réversible pour la désorption/resorption de l'eau.
• Profils de Rétrécissement : Le modèle a prédit les tendances de rétrécissement avec une erreur moyenne d'environ 10 % sur la forme finale de la particule. Bien que le modèle ait capturé la tendance générale de l'évolution de la forme, il a souligné qu'une compréhension fondamentale plus approfondie de l'évolution structurelle interne est encore nécessaire pour prédire parfaitement les formes finales.
• Cylindre Anisotrope : Pour un cylindre de hêtre en convection forcée, le modèle a reproduit avec succès la déformation « en forme d'os » observée dans les expériences, où les coins se dégradent plus rapidement que le centre en raison de l'anisotropie du transfert de chaleur et des conditions limites. Les coefficients d'échange de chaleur calculés (convection et rayonnement) étaient en adéquation avec les estimations expérimentales.

Signification et Revendications
Les auteurs affirment que ce travail constitue une avancée pour aider au développement de processus de pyrolyse durables en fournissant une description physique plus rigoureuse du comportement des particules. La capacité du modèle à résoudre l'interface et la dynamique interne sans corrélations empiriques en fait un outil puissant pour :

• La Modélisation Multi-échelle : Servir de générateur de modèles sous-maille pour dériver des corrélations spécifiques pour des simulations de réacteurs à plus grande échelle (ex: lits bouillonnants, réacteurs à vis sans fin) où la résolution de l'interface est trop coûteuse en calcul.
• L'Optimisation des Procédés : Prédire les rendements, les temps de dégradation et les profils de température sous des conditions de dynamique des fluides complexes.
• L'Apport de Compréhension Fondamentale : Fournir des informations sur l'interaction entre les changements de porosité et le rétrécissement de volume, qui sont difficiles à mesurer expérimentalement.

L'article conclut avec modestie, notant que bien que le modèle capture les principales tendances expérimentales, la forme fonctionnelle précise du couplage rétrécissement-porosité ($Z$) reste un domaine nécessitant des investigations fondamentales supplémentaires. Des travaux futurs sont identifiés comme nécessaires pour étendre le modèle aux conditions de combustion et de gazéification en incorporant des réactions en phase gazeuse et en gérant des mécanismes chimiques plus vastes.