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

Este artículo presenta un novedoso modelo de una sola malla, de tipo Euleriano-VOF, dentro del marco de código abierto Basilisk que resuelve plenamente la dinámica acoplada sólido-gas y la deformación anisotrópica de partículas durante la pirólisis de biomasa, demostrando una excelente concordancia con los datos experimentales al tiempo que proporciona una herramienta robusta para el desarrollo de procesos de pirólisis sostenibles.

Lo esencial

La visión general: Cocinar un tronco sin una receta

Imagina que estás intentando cocinar perfectamente un tronco de madera en un fuego para convertirlo en combustible (un proceso llamado pirólisis). Para hacer esto bien, necesitas saber exactamente qué está pasando dentro del tronco a medida que se calienta.

Durante mucho tiempo, los científicos tuvieron dos formas separadas de ver este problema:

1. La visión del "interior": Observaban cómo la madera se encoge y cambia por dentro, pero suponían cómo el aire caliente exterior la tocaba.
2. La visión del "exterior": Observaban cómo el aire caliente se movía alrededor del tronco, pero trataban al tronco como una roca estática que nunca cambiaba de forma.

El problema es que la madera no es una roca. Mientras se cocina, se encoge, se vuelve esponjosa (porosa) y el aire caliente que pasa a su alrededor cambia porque la forma del tronco está cambiando. Los métodos antiguos no captaban la conversación entre el "interior" y el "exterior".

La nueva solución: Una única cámara inteligente

Este artículo presenta un nuevo modelo computacional que actúa como una única cámara de alta definición que observa toda la escena a la vez. No supone cómo interactúan el aire y la madera; calcula la danza exacta entre ambos.

Así es como los autores construyeron esta "cámara":

1. El truco del "Volumen de Fluido" (La analogía del globo de agua)

Normalmente, las computadoras tienen dificultades para rastrear un límite en movimiento, como un globo que se encoge. Este modelo utiliza un método llamado Volumen de Fluido (VOF).

• La analogía: Imagina una cuadrícula de pequeñas cajas cubriendo tu pantalla. Algunas cajas están llenas de "madera", otras de "aire" y algunas son una mezcla. A medida que la madera se encoge, el modelo simplemente actualiza el porcentaje de "madera" en cada caja. Rastrea el borde de la madera a medida que se mueve, tal como se rastrea el borde de un globo de agua que es apretado.

2. El efecto "Esponja" (Porosidad y contracción)

La madera es como una esponja. Cuando se calienta, ocurren dos cosas simultáneamente:

• La esponja desarrolla agujeros: El material en el interior se descompone, creando más espacio vacío (porosidad).
• La esponja se hace más pequeña: Todo el tronco se encoge de tamaño.

Los autores crearon una regla especial (una función matemática que llaman Z) para decidir cuánto de la reacción causa que la madera desarrolle agujeros frente a cuánto causa que se encoja. Es como decidir si un cubo de hielo derritiéndose se convierte en un charco (desarrollando agujeros) o simplemente se hace más pequeño (encogiéndose). Descubrieron que los mejores resultados provienen de una mezcla de ambos.

3. El "Atasco de tráfico" (Flujo de gas en el interior)

A medida que la madera se cocina, libera gases. Estos gases tienen que escurrirse a través de los diminutos agujeros dentro de la madera para salir.

• La analogía: Imagina personas intentando salir de un estadio abarrotado. Si las salidas del estadio están muy abiertas, corren rápido. Si las salidas son estrechas y están congestionadas, se mueven lentamente. El modelo utiliza las ecuaciones de Darcy-Forchheimer para calcular este efecto de "atasco de tráfico", asegurando que el gas no aparezca mágicamente afuera, sino que realmente se abra paso a través de los poros de la madera.

4. La "Veta de la madera" (Anisotropía)

La madera no es igual en todas las direcciones. El calor viaja más rápido a lo largo de la veta (como correr por un pasillo) que a través de ella (como atravesar una multitud).

• La analogía: Piensa en una pila de papeles. Es fácil deslizar un dedo a lo largo de la pila (rápido), pero difícil empujar a través de la pila (lento). El modelo tiene en cuenta esto haciendo que el calor y el flujo de gas sean más rápidos en la dirección de las fibras de la madera y más lentos a través de ellas.

¿Qué probaron?

El equipo probó su modelo con experimentos del mundo real utilizando partículas de madera que iban desde pequeñas esferas hasta cilindros. Comprobaron:

• Temperatura: ¿Predice el modelo que la madera se calienta a la velocidad correcta? (Sí, coincidió bien).
• Pérdida de masa: ¿Predice el modelo cuánta madera se convierte en gas frente a carbón? (Sí, dentro de un margen de error muy pequeño).
• Cambio de forma: ¿Muestra el modelo que la madera se encoge correctamente? (Sí, aunque predecir la forma final exacta sigue siendo un poco complicado, la tendencia general fue correcta).

La conclusión fundamental

Este artículo presenta una nueva herramienta unificada que deja de suponer cómo la madera se encoge y cómo el aire se mueve a su alrededor. En su lugar, simula todo el proceso de una sola vez.

• Por qué es importante: Ayuda a los ingenieros a diseñar mejores sistemas para convertir la madera en energía renovable.
• El inconveniente: El modelo es complejo y requiere mucha potencia de cómputo, pero los autores han puesto su código como código abierto (gratis para que cualquiera lo use y lo mejore).

En resumen, construyeron un gemelo digital de un trozo de madera ardiendo que entiende tanto el interior como el exterior, permitiendo a los científicos ver los cambios "invisibles" que ocurren dentro de la madera mientras se convierte en combustible.

Resumen técnico

Resumen Técnico: Modelado Multifásico de la Pirólisis de Partículas de Biomasa Anisotrópicas

Planteamiento del Problema
Los modelos numéricos actuales para la pirólisis de partículas de biomasa sufren habitualmente de una limitación fundamental: se centran ya sea en la evolución de la partícula sólida (estructura interna, contracción, porosidad) o en la dinámica de la fase gaseosa circundante, pero rara vez en las interacciones acopladas entre ambas. Los enfoques existentes suelen depender de correlaciones de sub-escala para estimar los coeficientes de transferencia de calor y masa en el límite de la partícula, asumiendo una interfaz estacionaria. Esta suposición simplifica la complejidad del modelo, pero no logra capturar el impacto significativo de la variación del diámetro de la partícula y el movimiento de la interfaz en las tasas de intercambio, particularmente para la biomasa, que experimenta cambios morfológicos sustanciales (contracción y evolución de la porosidad) durante la conversión. Además, la naturaleza anisotrópica de las propiedades de la biomasa, como la conductividad térmica, suele ser omitida en los marcos multidimensionales que resuelven el flujo externo.

Metodología
Para abordar estas limitaciones, los autores proponen un modelo euleriano unificado de rejilla única que resuelve plenamente tanto la partícula de biomasa sólida como la fase gaseosa circundante sin depender de correlaciones de sub-escala. La metodología central incluye:

• Enfoque de Volumen de Fluido (VOF): El sistema bifásico se rastrea utilizando un método VOF para resolver la interfaz nítida entre el gas externo y la biomasa. La biomasa es tratada como una "pseudo-fase" que contiene una mezcla de una matriz sólida y gas interno dentro de los poros.
• Ecuaciones de Gobierno Acopladas: El modelo resuelve las ecuaciones de conservación de masa, cantidad de movimiento, energía y especies químicas tanto para el entorno externo como para la pseudo-fase interna.
  • Cantidad de Movimiento: Se incorporan términos de Darcy–Forchheimer en las ecuaciones de Navier-Stokes para dar cuenta de los efectos de arrastre dentro de la matriz porosa.
  • Energía: Se resuelve una única ecuación de energía para la pseudo-fase asumiendo un equilibrio térmico entre el gas interno y el sólido. La conductividad térmica anisotrópica se modela explícitamente.
  • Química: Las reacciones de pirólisis de la fase sólida se modelan mediante un mecanismo cinético detallado (mecanismo CRECK) gestionado a través de la biblioteca OpenSMOKE++. Las reacciones de fase gaseosa se omiten para las condiciones específicas de pirólisis de baja temperatura probadas, pero se reconocen como necesarias para futuros estudios de combustión/gasificación.
• Acoplamiento de Contracción y Porosidad: Se introduce un enfoque novedoso para acoplar la evolución de la porosidad de la biomasa con la contracción de la partícula. El modelo introduce una función, $Z(x,t)$, que particiona el término de fuente de la reacción química entre la reducción de volumen (contracción) y los cambios de porosidad interna. Esta función se resuelve junto con una aproximación de flujo potencial para determinar la velocidad de contracción de la matriz porosa.
• Implementación Numérica: El modelo se implementa en el marco de código abierto Basilisk mediante un método de volúmenes finitos con refinamiento de malla adaptativo (rejilla quad/octree). Un método de proyección aproximado de segundo orden resuelve las ecuaciones de Navier-Stokes, mientras que el transporte escalar (temperatura, especies) se gestiona mediante un enfoque de dos campos para minimizar la difusión numérica a través de la interfaz.

Contribuciones Clave

1. Marco Multifásico Unificado: El trabajo presenta un modelo que resuelve la interfaz gas-sólido de forma dinámica, eliminando la necesidad de correlaciones empíricas de transferencia de calor y masa.
2. Modelado de Partículas Deformables y Anisotrópicas: El marco tiene en cuenta la naturaleza anisotrópica de la biomasa (conductividad térmica y permeabilidad dependientes de la dirección) y acopla los cambios de porosidad interna con la contracción externa de la partícula en una única simulación multidimensional.
3. Conservación de Masa y Convergencia: El modelo demuestra la conservación de la masa y la convergencia numérica. Aunque la interfaz móvil resulta en una convergencia de primer orden para la conservación de la masa (debido a flujos no divergentes libres), las interfaces estacionarias muestran una convergencia de segundo orden.
4. Disponibilidad de Código Abierto: El código y las configuraciones de simulación se ponen a disposición pública dentro del marco Basilisk para mejorar la reproducibilidad.

Resultados y Validación
El modelo fue validado mediante extensos datos experimentales que involucran partículas de madera a escala de centímetros, operando a temperaturas entre 400 °C y 700 °C.

• Esfera Isotrópica (Temperatura Impuesta): El modelo predijo con precisión los perfiles de temperatura (núcleo, radio medio, superficie) y las tasas de liberación de especies volátiles (CO, CO2, CH4, etc.) para una esfera de madera de álamo. Capturó con éxito los picos de descomposición endotérmica y los picos de formación de carbón exotérmico.
• Esfera de Madera en Flujo: Las simulaciones de esferas de madera de Schima superba en un flujo laminar de nitrógeno caliente mostraron un excelente acuerdo con los datos experimentales de pérdida de masa y perfiles de temperatura del núcleo. El modelo capturó correctamente la meseta de temperatura causada por la evaporación de la humedad (aprox. 373 K) mediante la implementación de un mecanismo de reacción reversible para la desorción/resorción de agua.
• Perfiles de Contracción: El modelo predijo las tendencias de contracción con un error promedio de aproximadamente el 10% en la forma final de la partícula. Aunque el modelo capturó la tendencia general de la evolución de la forma, destacó que aún se requiere una comprensión fundamental más profunda de la evolución estructural interna para predecir perfectamente las formas finales.
• Cilindro Anisotrópico: Para un cilindro de madera de haya en convección forzada, el modelo reprodujo con éxito la deformación "tipo hueso" observada en los experimentos, donde las esquinas se degradan más rápido que el centro debido a la transferencia de calor anisotrópica y las condiciones de contorno. Los coeficientes de intercambio de calor calculados (convección y radiación) se alinearon bien con las estimaciones experimentales.

Significancia y Reivindicaciones
Los autores afirman que este trabajo representa un paso adelante en el apoyo al desarrollo de procesos de pirólisis sostenibles al proporcionar una descripción física más riguro de la conducta de las partículas. La capacidad del modelo para resolver la interfaz y la dinámica interna sin correlaciones empíricas lo convierte en una herramienta poderosa para:

• Modelado Multiescala: Servir como generador de modelos de sub-escala para derivar correlaciones específicas para simulaciones de reactores a mayor escala (por ejemplo, lechos burbujeantes, reactores de tornillo sin fin) donde la resolución de la interfaz es computacionalmente prohibitiva.
• Optimización de Procesos: Predecir rendimientos, tiempos de degradación y perfiles de temperatura bajo condiciones de dinámica de fluidos complejas.
• Perspectiva Fundamental: Proporcionar información sobre la interacción entre los cambios de porosidad y la contracción de volumen, que son difíciles de medir experimentalmente.

El artículo concluye con modestia, señalando que aunque el modelo captura las principales tendencias experimentales, la forma funcional precisa del acoplamiento entre la contracción y la porosidad ($Z$) sigue siendo un área que requiere mayor investigación fundamental. Se identifica como trabajo futuro la necesidad de extender el modelo a condiciones de combustión y gasificación mediante la incorporación de reacciones de fase gaseosa y el manejo de mecanismos químicos más complejos.