📄 chemistry

Effective diffusion in sharp-front reactive transport experiments

Este estudio presenta una solución analítica de forma cerrada para el transporte reactivo de frente agudo en medios porosos, demostrando que el coeficiente de difusión efectiva inferido experimentalmente difiere fundamentalmente del coeficiente de difusión molecular intrínseca debido a efectos acoplados de advección-difusión-termodinámica, un marco validado mediante experimentos de carbonatación de cemento Portland bajo diversas condiciones de flujo.

Autores originales: Viktoriya M. Yarushina, Reinier van Noort, Yury Y. Podladchikov, Yannik F. Schneider

Publicado 2026-09-01

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Viktoriya M. Yarushina, Reinier van Noort, Yury Y. Podladchikov, Yannik F. Schneider

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

En lo profundo del subsuelo, donde la roca y el fluido se encuentran, se libra constantemente una batalla química silenciosa. Cuando el agua o el gas ricos en dióxido de carbono fluyen a través de materiales porosos como el cemento o la roca, desencadenan una reacción que transforma la estructura sólida. Este proceso no ocurre de manera uniforme en todo el material; en su lugar, crea un límite distintivo, un frente nítido, que se desplaza hacia el interior, separando la zona alterada y reaccionada del material original intacto. Este fenómeno es crítico para la seguridad del almacenamiento geológico de carbono, donde se utilizan sellos de cemento para atrapar gases de efecto invernadero en las profundidades del subsuelo. Si estos sellos se degradan demasiado rápido, el gas almacenado podría escapar. Durante décadas, los científicos han intentado medir qué tan rápido se mueve este frente para predecir la vida útil de estos sellos, asumiendo a menudo que la velocidad es impulsada por la facilidad con la que las moléculas se desplazan a través de los diminutos poros del material, un proceso conocido como difusión.

Sin embargo, un nuevo estudio realizado por investigadores del Instituto de Tecnología de la Energía y la Universidad de Lausana desafía la forma en que interpretamos estas mediciones. Descubrieron que la velocidad a la que viaja este frente de reacción no es un simple reflejo de qué tan rápido se mueven las moléculas a través del fluido. En cambio, el movimiento está gobernado por una compleja interacción entre el flujo del fluido, la difusión de las moléculas y un efecto de "amortiguación" química donde el material sólido absorbe la sustancia reactiva. Los investigadores desarrollaron un nuevo marco matemático para describir este movimiento, demostrando que la velocidad "efectiva" medida en los experimentos es, en realidad, un resultado compuesto del movimiento del fluido y de la enorme cantidad de almacenamiento químico que el material sólido puede contener. Esta distinción es vital porque explica por qué diferentes experimentos, incluso aquellos que utilizan métodos muy distintos, suelen producir números confusamente diferentes sobre qué tan rápido se degrada el cemento.

Para probar esta teoría, el equipo recurrió al cemento Portland, el material utilizado para sellar pozos para el almacenamiento de carbono. Sometieron muestras de cemento idénticas a dos entornos muy diferentes. En un conjunto de experimentos, el cemento reposaba tranquilamente en un contenedor lleno de agua saturada de dióxido de carbono o de dióxido de carbono supercrítico húmedo. Esto se conoce como un experimento de lote (batch), donde no hay flujo forzado y el fluido simplemente permanece alrededor de la muestra. En el segundo conjunto, forzaron a estos mismos fluidos a fluir a través de las muestras de cemento bajo presión, simulando un escenario donde el fluido se mueve activamente a través de un pozo. La sabiduría convencional sugeriría que el flujo forzado empujaría el frente de reacción mucho más rápido que las condiciones estáticas y tranquilas. Sin embargo, cuando los investigadores midieron la profundidad de la carbonatación después de semanas y meses, los resultados fueron sorprendentemente similares. En ambos escenarios, tanto en el de flujo tranquilo como en el de flujo forzado, el frente de reacción se movió a un ritmo que seguía el mismo patrón: la distancia recorrida aumentaba con la raíz cuadrada del tiempo. Este patrón es la marca distintiva de un proceso controlado por la difusión, no por la velocidad del fluido circulante.

Los investigadores utilizaron esta sorprendente similitud para explicar un enigma de larga data en el campo. Cuando los científicos miden qué tan rápido se mueve un frente de reacción en un experimento, a menudo calculan un "coeficiente de difusión efectivo" y lo tratan como una propiedad fija del material, como la densidad de una roca. El nuevo estudio muestra que esto es un malentendido. La velocidad del frente no está determinada únicamente por qué tan rápido las moléculas de dióxido de carbono pueden zigzaguear a través del agua en los poros. También está determinada por cuánto de ese dióxido de carbono puede absorber y almacenar el cemento sólido al reaccionar. El cemento actúa como una esponja, absorbiendo el dióxido de carbono y convirtiéndolo en un nuevo mineral. Debido a que el sólido contiene mucho más carbono que el fluido, el frente se mueve lentamente, no porque el fluido sea lento, sino porque la reacción tiene que "llenar" la esponja sólida antes de poder avanzar. Los investigadores introdujeron una nueva forma de ver esto, utilizando números adimensionales para describir el equilibrio entre el flujo del fluido y la capacidad química del sólido. Encontraron que incluso cuando el fluido es forzado a través del cemento, la capacidad de almacenamiento químico es tan grande que el flujo apenas hace diferencia en la velocidad del frente.

Para confirmar que el frente de reacción era, de hecho, agudo y distintivo, el equipo examinó las muestras de cemento bajo potentes microscopios y utilizó luz infrarroja para mapear los cambios químicos. Observaron el límite entre el cemento reaccionado y el no reaccionado con extrema precisión. Las imágenes revelaron que la transición era increíblemente nítida, ocurriendo en una distancia de solo unos 50 a 100 micrómetros, que es aproximadamente el ancho de un cabello humano. En esta zona estrecha, la composición química cambiaba abruptamente. El silicato de calcio hidratado, el pegamento que mantiene unido al cemento, desaparecía y era reemplazado por carbonato de calcio. Los investigadores también observaron que el calcio en el cemento no se disolvía ni se lavaba; en su lugar, permanecía en su lugar, reorganizándose en la nueva estructura mineral. Esta observación respalda la idea de que el frente es una interfaz aguda donde la fase sólida actúa como un amortiguador estacionario, absorbiendo el fluido entrante y reaccionando localmente.

Las implicaciones de estos hallazgos son significativas para la forma en que diseñamos y probamos materiales para el almacenamiento de carbono. El estudio demuestra que simplemente aumentar el flujo de fluido a través de un pozo no acelerará necesariamente la degradación del sello de cemento, porque el proceso está dominado por la capacidad de almacenamiento químico del material más que por la velocidad del fluido. Esto significa que los experimentos realizados en contenedores estáticos y tranquilos pueden compararse directamente con aquellos con flujo forzado, siempre que las condiciones químicas sean similares. La velocidad "efectiva" medida en el laboratorio no es una propiedad fundamental de la capacidad del fluido para moverse, sino un resultado de la danza química específica entre el fluido y el sólido. Al comprender esto, los científicos pueden predecir mejor cuánto tiempo durará un sello de cemento en el entorno hostil de un pozo profundo, asegurando que el dióxido de carbono permanezca atrapado de forma segura bajo tierra durante miles de años. El trabajo aclara que el movimiento de estos frentes químicos es una historia de almacenamiento y equilibrio, no solo de flujo.

Resumen Técnico: Difusión Efectiva en Experimentos de Transporte Reactivo de Frente Agudo

Planteamiento del Problema
El transporte reactivo en materiales porosos, como los sistemas de cemento expuestos a CO₂, genera frecuentemente frentes de alteración agudos que separan las regiones reaccionadas de las no alteradas. Dado que estos frentes son críticos para comprender las interacciones fluido-roca en el almacenamiento geológico de CO₂ y la integridad de pozos, sus tasas de propagación se interpretan a menudo mediante leyes de difusión empíricas (p. ej., $x \propto \sqrt{t}$ ). Existe una ambigüedad significativa en la literatura con respecto al "coeficiente de difusión efectivo" ( $D_{eff}$ ) inferido de estos experimentos. Los valores reportados varían ampliamente y suelen ser órdenes de magnitud menores que los coeficientes de difusión molecular intrínsecos, lo que genera confusión sobre si estos valores representan propiedades reales del material o son artefactos de los procesos acoplados de transporte-reacción. Además, no está claro cómo comparar consistentemente experimentos realizados bajo diferentes condiciones de contorno (p. ej., lote vs. flujo forzado) o cómo la advección influye en la propagación del frente cuando el amortiguamiento químico es fuerte.

Metodología
Los autores desarrollan una solución analítica de forma cerrada para el transporte reactivo de advección-difusión que presenta un frente de reacción móvil. El modelo asume:

Conservación de la Masa: El transporte de una especie reactiva está gobernado por la advección y la difusión en la fase fluida, acoplado con el intercambio con la fase sólida mediante partición en equilibrio.
Interfaz Aguda: Una frontera móvil (tipo Stefan) separa las zonas reaccionadas de las no reaccionadas, a través de la cual la porosidad, la densidad y las concentraciones presentan saltos finitos.
Condición de Rankine-Hugoniot: La velocidad del frente se deriva de la conservación de la masa en la interfaz, considerando las discontinuidades en la concentración total.

La solución introduce dos números adimensionales clave:

Número de Péclet ( $Pe$ ): Controla la importancia relativa de la advección frente a la difusión.
Número de Stefan Químico ($Ste$): Controla la tasa de propagación del frente cuantificando la relación entre el flujo impulsado por la concentración móvil y la acumulación total de especies requerida a través del frente (incorporando el almacenamiento en la fase sólida y los saltos de concentración).

El marco es validado y aplicado a datos experimentales de muestras de cemento Portland expuestas a fluidos ricos en CO₂, tanto en condiciones de lote como de flujo forzado. Los experimentos utilizaron CO₂ supercrítico húmedo y agua saturada con CO₂. La caracterización se realizó mediante Microscopía Electrónica de Barrido con Electrones Retrodispersados (SEM-BSE), Espectroscopía de Energía Dispersiva (EDS) y microscopía de Infrarrojo por Transformada de Fourier (FTIR) para resolver la morfología del frente y los gradientes químicos.

Contribuciones Clave

Marco Analítico: El artículo proporciona una solución analítica de tipo Stefan que explica explícitamente las discontinuidades en la porosidad, la densidad y las concentraciones en una interfaz móvil. Esta solución unifica los regímenes dominados por difusión y advección a través de un parámetro compuesto $\Lambda = Pe / (Ste \cdot R)$ , donde $R$ es el factor de retardo.
Redefinición de la Difusión Efectiva: El estudio demuestra que el "coeficiente de difusión efectivo" inferido de la propagación del frente no es una propiedad intrínseca del fluido de los poros. En cambio, es un coeficiente de propagación compuesto que depende del coeficiente de difusión molecular ( $D_f$ ) escalado por el número de Stefan químico ( $D_{eff} = D_f \cdot Ste$ ).
Resolución de Discrepancias: El marco explica por qué los coeficientes de difusión inferidos experimentalmente son a menudo varios órdenes de magnitud menores que los valores moleculares: un número de Stefan pequeño implica que solo una pequeña fracción del flujo difusivo contribuye al avance del frente, ya que la mayor parte es consumida por el almacenamiento en la fase sólida.
Validación Experimental: Los autores muestran que, a pesar de las metodologías de exposición fundamentalmente diferentes (lote vs. flujo forzado), las tasas de propagación del frente pueden ser similares si los parámetros adimensionales gobernantes son comparables.

Resultados

Propagación Controlada por Difusión: Tanto en los experimentos de lote como en los de flujo forzado con cemento Portland, la propagación del frente de carbonatación siguió una escala de raíz cuadrada del tiempo ( $x \propto \sqrt{t}$ ). Los números de Péclet calculados fueron bajos ( $Pe \approx 0.2$ para CO₂ supercrítico húmedo y menores para agua saturada con CO₂), lo que indica que, incluso bajo flujo impuesto, el transporte siguió siendo dominado por la difusión.
Coeficientes Intrínsecos vs. Efectivos: Los coeficientes de difusión efectivos inferidos fueron de $1.4 \times 10^{-11}$ m²/s (agua saturada con CO₂) y $2.5 \times 10^{-11}$ m²/s (CO₂ supercrítico húmedo). Estos son órdenes de magnitud menores que la difusión molecular intrínseca del CO₂ en agua ( $\sim 10^{-8}$ m²/s). Los autores atribuyen esto a números de Stefan pequeños ( $Ste \sim 10^{-3}$ ), lo que refleja la alta capacidad de carbonatación de la fase sólida en relación con el contraste de concentración del fluido.
Morfología del Frente: Las observaciones de alta resolución mediante SEM-BSE, EDS y FTIR confirmaron que los frentes de reacción son extremadamente agudos, con zonas de transición localizadas a escalas micrométricas (aprox. 50–110 µm).
- El análisis FTIR reveló que la aparición de carbonatos y la desaparición del gel C-S-H ocurren en zonas estrechas y distintas.
- Los perfiles de EDS mostraron transiciones abruptas en las relaciones Ca/Si, lo que respalda la suposición de una interfaz aguda y un equilibrio local.
- En los experimentos de flujo forzado, se observó una estrecha zona de enriquecimiento de Ca/Si inmediatamente detrás del frente, pero el frente en general se mantuvo agudo, validando la aproximación de calcio inmóvil en el régimen de bajo- $Pe$ .

Significancia y Reivindicaciones
El artículo afirma proporcionar una base general y físicamente transparente para interpretar la propagación de frentes de alteración en medios porosos reactivos. Al distinguir entre la difusión molecular intrínseca y el coeficiente de propagación de frente efectivo, el estudio resuelve la aparente contradicción entre los diversos valores de difusión experimental y las propiedades de difusión molecular consistentes. Los autores sostienen que su marco permite la comparación consistente de experimentos realizados bajo diferentes condiciones de contorno (lote vs. paso de flujo) y composiciones de fluidos. Los resultados respaldan la aplicabilidad de las formulaciones de interfaz aguda y equilibrio local para modelar el transporte reactivo en sistemas que van desde barreras de cemento de ingeniería hasta zonas naturales de interacción fluido-roca, siempre que los números adimensionales gobernantes ( $Pe$ , $Ste $y$ \Lambda$) se caractericen adecuadamente.

Más como este