🔬 física

Dimensional effects on breakthrough behavior and pore-scale residence time distributions

Este estudio utiliza simulaciones de Lattice Boltzmann para demostrar que, si bien los modelos a escala de poros en 2D subestiman sistemáticamente la permeabilidad en 3D debido a la reducción de la conectividad fuera del plano —una discrepancia mejor predicha por la característica de Euler—, el impacto de la dimensionalidad en las métricas de transporte, como las curvas de ruptura y los tiempos de residencia, varía significativamente dependiendo de la geometría de poros específica y la porosidad del medio.

Autores originales: Shahram Asgari, Christof Meile

Publicado 2026-10-16

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Shahram Asgari, Christof Meile

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

El agua y las sustancias disueltas se mueven a través del mundo oculto e intrincado del suelo, la roca y los materiales de ingeniería de formas que a menudo son imposibles de ver. Este movimiento, conocido como transporte de solutos, es crítico para todo, desde la limpieza de aguas subterráneas contaminadas hasta el almacenamiento de carbono en las profundidades del subsuelo. Para comprender cómo estos fluidos navegan por los túneles diminutos y serpenteantes entre los granos de arena o las fibras de un filtro, los científicos construyen modelos computacionales que imitan el flujo a la escala de los poros individuales. Durante décadas, un dilema práctico ha enfrentado a los investigadores: si bien los modelos tridimensionales ofrecen la imagen más realista de estas redes complejas, exigen una potencia de cálculo inmensa. Los modelos bidimensionales, que son esencialmente secciones transversales planas del material, son mucho más económicos y rápidos de ejecutar, lo que ha llevado a muchos a utilizarlos como sustitutos de la realidad tridimensional completa. La pregunta ha sido durante mucho tiempo si este atajo sacrifica demasiada precisión para ser útil.

Un equipo de investigadores de la Universidad de Georgia se propuso probar exactamente cuánto pierden estos modelos planos. Crearon versiones digitales de tres tipos distintos de materiales porosos: uno hecho de fibras enredadas, otro que recordaba a esferas compactadas como una pila de canicas, y un tercero que imitaba la estructura aleatoria y grumosa de las formaciones geológicas naturales. Para cada material, construyeron tanto un bloque tridimensional completo como una serie de secciones bidimensionales cortadas de este. Utilizando un método de simulación potente que rastrea el movimiento de las partículas de fluido y las sustancias disueltas, compararon cómo el agua y los trazadores se movían a través de las secciones planas frente a los bloques completos. Observaron específicamente con qué facilidad podía pasar el fluido a través del material, cuánto tiempo tardaba una sustancia en viajar de un lado al otro y cuánto se dispersaba la sustancia en el camino.

Los resultados confirmaron que los modelos planos subestiman consistentemente la facilidad con la que el fluido puede fluir a través del material. En los bloques tridimensionales, el fluido tiene la libertad de moverse alrededor de los obstáculos en todas las direcciones, encontrando atajos que simplemente no existen en una sección plana. Cuando los investigadores observaron las secciones bidimensionales, descubrieron que el fluido se veía forzado a seguir caminos más largos y sinuosos porque no podía sortear los obstáculos moviéndose hacia "arriba" o hacia "abajo", fuera del plano. Este efecto fue más dramático en materiales con baja porosidad, donde el espacio disponible para el flujo ya es reducido. En estos casos, los modelos planos predijeron valores de permeabilidad casi diez veces menores que la verdadera realidad tridimensional para los materiales fibrosos, y aproximadamente cuatro veces menores para las estructuras de tipo geológico. Los investigadores descubrieron que esta discrepancia no es solo una cuestión de cuánto espacio vacío existe, sino más bien de una pérdida fundamental de conexión. Al aplanar el mundo, los modelos cortan los vitales enlaces fuera del plano que permiten al fluido tejerse a través de la red de manera eficiente.

Para predecir cuándo un modelo plano podría ser un sustituto seguro, el equipo recurrió a una forma matemática de contar las conexiones dentro de la red de poros. Utilizaron una medida que esencialmente cuenta el número de bolsas separadas, bucles y espacios cerrados en la estructura. Descubrieron que cuando la sección plana y el bloque completo tenían conteos de estas conexiones muy diferentes, las predicciones de flujo divergían dramente. Esta medida topológica actúa como una herramienta de diagnóstico rápido: si una sección plana captura el mismo patrón de bucles y conexiones que el bloque completo, es probable que proporcione una imagen precisa del flujo. Si los conteos difieren significamente, el modelo plano está perdiendo vías críticas, y es necesaria una simulación tridimensional completa.

Más allá de la velocidad a la que fluye el agua, el estudio examinó cuánto tiempo tarda una sustancia en atravesar el material y cuánto se dispersa. El comportamiento de estas métricas de transporte dependía fuertemente del tipo de material. En el medio fibroso, la diferencia entre los modelos plano y completo se mantuvo constante independientemente de cuán permeable fuera el material, probablemente debido a que el flujo está dominado por unos pocos canales principales que se ven similares en ambas dimensiones. En las estructuras de tipo geológico, la brecha entre los dos modelos se reducía a medida que el material se volvía más permeable, sugiriendo que en materiales más sueltos y abiertos, las secciones planas capturan bien las rutas de flujo. Sin embargo, en los modelos de esferas compactadas, las secciones planas a veces predecían tiempos de viaje más rápidos y otras veces más lentos, sin un patrón claro, lo que indica que la disposición específica de las esferas en una sola sección puede ser engañosa.

Los investigadores también probaron si la naturaleza química de la sustancia rastreada cambiaba estos resultados. Compararon un trazador conservativo, que simplemente se mueve con el agua sin reaccionar, contra un trazador reactivo que interactúa con el material a medida que se desplaza. Sorprendentemente, el tipo de trazador apenas influyó en la comparación general. Ya fuera que la sustancia fuera inerte o reactiva, los errores introducidos por el uso de un modelo plano estaban gobernados principalmente por la forma y la conectividad de los poros, no por la química del fluido. Esto sugiere que las limitaciones estructurales de los modelos bidimensionales son un problema geomético fundamental, no uno químico.

En última instancia, el estudio proporciona una guía clara para científicos e ingenieros que deben elegir entre velocidad y precisión. Si bien los modelos bidimensionales son útiles para estimaciones rápidas, omiten sistemáticamente la facilidad de flujo que se encuentra en la realidad tridimensional porque no pueden contabilizar las rutas adicionales disponibles en la tercera dimensión. La investigación muestra que, al verificar una medida simple de cómo están conectados los poros, se puede determinar antes de ejecutar una simulación si una sección plana es un sustituto fiel o una simplificación engañosa. Este conocimiento permite un uso más eficiente de los recursos computacionales, asegurando que la naturaleza compleja y tridimensional del mundo natural sea respetada cuando más importa.

Resumen Técnico: Efectos Dimensionales en el Comportamiento de Ruptura y la Distribución del Tiempo de Residencia a Escala de Poros

Planteamiento del Problema
Los modelos de poros a escala de dos dimensiones (2D) se emplean frecuentemente como sustitutos computacionalmente eficientes de las simulaciones tridimensionales (3D) en el estudio del transporte de solutos dentro de medios porosos. Sin embargo, las estructuras porosas 2D y 3D exhiben comportamientos de flujo y transporte fundamentalmente diferentes debido a las diferencias en conectividad, tortuosidad y regímenes de mezcla. Mientras que las simulaciones 3D ofrecen la representación más realista de las redes de poros, son computacionalmente intensivas. En consecuencia, los investigadores suelen recurrir a cortes 2D derivados de imágenes de micro-CT o secciones delgadas. El problema central abordado en este estudio es la medida en que estos análogos 2D pueden predecir con precisión la permeabilidad 3D, las curvas de ruptura y las distribuciones de la edad del fluido, y si indicadores topológicos o geométricos específicos pueden determinar cuándo un modelo 2D es un sustituto adecuado para una simulación 3D completa.

Metodología
Los autores utilizaron simulaciones de Lattice Boltzmann (LB) para resolver las ecuaciones de Navier–Stokes y de advección–difusión–reacción (ADR) en rejillas de vóxeles discretas. El estudio se centró en tres tipos distintos de medios porosos sintéticos:

Fibroso: Representa tejidos biológicos y materiales de ingeniería, generado mediante la colocación de fibras cilíndricas.
Gaussiano: Se asemeja a materiales geológicos como agregados de suelo y arenisca, generado mediante un campo aleatorio filtrado.
Esferas Superpuestas: Representa materiales granulares y lechos empacados, construido mediante la colocación aleatoria de esferas.

Para cada tipo de medio, se generaron dominios computacionales 3D (100 × 100 × 100 vóxeles) con porosidades baja, media y alta (0.4, 0.6 y 0.8). Los dominios 2D se extrajeron como cortes en el plano xz de estas estructuras 3D, asegurando la percolación del flujo.

Las simulaciones compararon:

Permeabilidad: Calculada mediante la ley de Darcy tanto para los cortes 2D como para sus dominios parentales 3D.
Métricas de Transporte: Se analizaron curvas de ruptura ponderadas por flujo para derivar el primer momento temporal ( $M_1$ , tiempo de llegada promedio) y el coeficiente de variación ( $CV$ , dispersión de la ruptura).
Tipos de Trazadores: Se simularon trazadores tanto conservativos (no reactivos) como reactivos (que representan la edad del fluido con una tasa de reacción de primer orden) bajo condiciones de advección dominante (números de Péclet de 4.1, 12.5 y 37.5).
Análisis Topológico: Se computó la característica de Euler normalizada por poros ( $\tilde{\chi}$ ), derivada de la homología persistente (números de Betti $\beta_0, \beta_1, \beta_2$ ), para cuantificar la conectividad de la red de poros como una métrica puramente geométrica.

Resultos Clave

Subestimación Sistemática de la Permeabilidad:
Las simulaciones 2D subestimaron sistemáticamente la verdadera permeabilidad 3D en todos los tipos de medios y porosidades. La magnitud de esta discrepancia fue más pronunciada en las porosidades bajas. Por ejemplo, en el medio Fibroso con $\phi=0.4$ , la permeabilidad 2D fue casi 10 veces menor que los valores 3D. Este sesgo se atribuyó a la pérdida de conectividad fuera del plano en los cortes 2D, lo que aumenta la tortuosidad efectiva y reduce la fracción de espacio poroso que contribuye al flujo activo. La relación de Kozeny–Carman confirmó que las diferencias de porosidad por sí solas no podían explicar esta discrepancia; el desconecte topológico fue el principal motor.
Papel de la Característica de Euler ( $\tilde{\chi}$ ):
La característica de Euler normalizada por poros ( $\tilde{\chi}$ ), una medida topológica de conectividad, rastreó en gran medida la diferencia entre la permeabilidad 2D y 3D. Los valores negativos de $\tilde{\chi}$ (que indican redes bien conectadas dominadas por bucles) se correlacionaron con permeabilidades más altas y menores discrepancias 2D–3D, mientras que los valores positivos de $\tilde{\chi}$ (redes fragmentadas) se correlacionaron con permeabilidades más bajas y mayores discrepancias. Los autores proponen $\tilde{\chi}$ como un proxy de bajo costo, basado únicamente en la geometría, para filtrar los cortes 2D antes de realizar simulaciones de flujo.
Tendencias de Transporte Específicas del Medio:
Los efectos dimensionales en las métricas de transporte ( $M_1$ y $CV$ ) variaron significamente según el tipo de medio:

Medio Fibroso: La diferencia relativa entre el tiempo de viaje promedio y la dispersión de la ruptura 2D y 3D se mantuvo casi constante a través del rango de permeabilidad. Esto se atribuye a que el flujo se organiza en canales dominantes que se preservan en 2D.
Medio Gaussiano: Las diferencias 2D–3D disminuyeron a medida que la permeabilidad aumentó. A baja permeabilidad, los cuellos estrechos restringen el flujo 2D, pero estos pueden ser evitados mediante rutas fuera del plano en 3D; a medida que la permeabilidad aumenta, estas restricciones se relajan, reduciendo la brecha dimensional.
Medio de Esferas Superpuestas: Las diferencias en el tiempo de llegada promedio fueron grandes y bidireccionales (algunos cortes 2D sobreestimaron, otros subestimaron los resultados 3D) sin una tendencia clara respecto a la permeabilidad. Esto refleja la compleja interacción entre los poros grandes e interconectados y la porosidad específica de los cortes 2D individuales.

Impacto de la Reactividad del Trazador:
Los trazadores conservativos y reactivos exhibieron sesgos dimensionales similares. Aunque los trazadores reactivos mostraron una diferencia relativa en el tiempo de viaje promedio ( $M_1$ ) entre un 20 y 40% mayor en comparación con los trazadores conservativos, el patrón general estuvo gobernado por la geometría de los poros más que por la reactividad del trazador. La dispersión de la ruptura ( $CV$ ) mostró más variabilidad, donde el medio Gaussiano presentó un sesgo distintivo en el que los trazadores conservativos exhibieron mayores diferencias 2D–3D que los reactivos.

Significancia y Reivindicaciones
El artículo concluye que la geometría y la conectividad de los poros son los principales reguladores del sesgo dimensional en el modelado a escala de poros, mientras que la reactividad del trazador desempeña un papel secundario. El estudio demuestra que, si bien los modelos 2D son computacionalmente eficientes, introducen errores sistemáticos en las predicciones de permeabilidad y transporte que no pueden corregirse únicamente igualando la porosidad.

Los autores afirman que la característica de Euler normalizada por poros ( $\tilde{\chi}$ ) ofrece una herramienta de cribado práctica previa a la simulación. Al comparar el $\tilde{\chi}$ de los candidatos a cortes 2D frente al sistema 3D parental, los investigadores pueden identificar qué cortes preservan la conectividad topológica del sistema 3D. Si los valores de $\tilde{\chi}$ coinciden en signo y magnitud, un corte 2D puede servir como un análogo válido; de lo contrario, se requiere una simulación 3D completa. Este enfoque sugiere que, al construir micromodelos 2D (por ejemplo, dispositivos microfluídicos) para representar sistemas 3D, igualar la conectividad a escala de poros mediante descriptores topológicos es tan crítico como igualar la porosidad o el área superficial.

Más como este