A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

En el mundo oculto bajo nuestros pies, el agua y los productos químicos disueltos se mueven a través de un complejo laberinto de roca. Esta roca rara vez es uniforme; a menudo es un bloque sólido, llamado matriz, atravesado por finas grietas conocidas como fracturas. Estas fracturas son las autopistas del subsuelo. Pueden actuar como canales superrápidos que hacen pasar fluidos a través de la roca, o como barreras que bloquean el movimiento por completo, dependiendo de su naturaleza. Debido a que estas grietas suelen ser de un ancho microscópico pero pueden extenderse por kilómetros de longitud, modelar cómo viajan las sustancias a través de ellas es un desafío masivo para los científicos. Si intentan construir un modelo computacional que capture cada pequeño detalle de la apertura de la grieta, el cálculo se vuelve tan pesado que incluso las supercomputadoras más potentes no pueden resolverlo. Para evitar esto, los investigadores utilizan un atajo ingenioso llamado modelo de dimensión híbrida. En lugar de tratar la grieta como una fina rebanada tridimensional, la aplanan en una línea bidimensional o una línea unidimensional que se asienta dentro de la roca. Esto permite simular el flujo sin necesidad de resolver el diminuto grosor de la grieta en sí.

La dificultad surge cuando los científicos intentan conectar las matemáticas de la roca sólida con las matemáticas de estas grietas aplanadas. En el mundo real, la roca y la grieta están en constante conversación, intercambiando fluidos y productos químicos a través de su frontera compartida. Para simular esto, las computadoras suelen dividir la roca y la grieta en piezas pequeñas, como un mosaico, para calcular el flujo. Durante años, el enfoque estándar requería que las piezas de la roca y las piezas de la grieta encajaran perfectamente a lo largo de su borde compartido, como dos piezas de un rompecabezas que encajan. Si la roca se dividía en cuadrados grandes y la grieta en segmentos diminutos, la computadora tenía que forzarlos a coincidir, a menudo creando una capa compleja y artificial de datos para cerrar la brecha. Este requisito hacía que los modelos fueran rígidos y difíciles de usar cuando diferentes partes de un sistema necesitaban ser refinadas a diferentes escalas.

Un nuevo estudio de Alvaro H. Salas y Carlos D. Acosta, de la Universidad Nacional de Colombia, aborda esta rigidez planteando una pregunta simple pero profunda: ¿qué pasaría si la roca y la grieta no tuvieran que alinearse en absoluto? Desarrollaron un método donde la roca y la fractura pueden dividirse en piezas de tamaños y formas completamente diferentes, generadas de manera independiente entre sí, y aun así intercambiar información perfectamente. Imagine intentar verter agua de un cubo grande y tosco en un tamiz fino y detallado. Normalmente, necesitaría un embudo especial para hacer que el agua se transfiera suavemente sin derramarse. Los investigadores encontraron una forma de verter el agua directamente, calculando la cantidad exacta que pasa por cada punto donde el borde del cubo toca el tamiz, independientemente de cuán desajustados estén los agujeros.

El equipo creó un marco matemático que trata el intercambio de productos químicos entre la roca y la fractura como una interacción física directa. En lugar de forzar a las dos mallas a coincidir o usar un intermediario para traducir los datos, construyeron un sistema que calcula la transferencia en las secciones diminutas y superpuestas donde el borde de la roca y la superficie de la fractura realmente se encuentran. Este enfoque se llama un método de malla no coincidente (nonmatching-grid method). En sus simulaciones, probaron este sistema con reacciones químicas complejas, donde las sustancias cambian a medida que se mueven, un escenario que es notoriamente difícil de modelar. Demostraron matemáticamente que esta conexión directa es estable y precisa, incluso cuando la roca se divide en bloques muy grandes y la fractura en segmentos muy finos, o viceversa. La estabilidad de su método no depende de la relación entre los tamaños de estos bloques, lo que significa que la computadora no se bloquea ni produce errores erráticos solo porque una malla sea mucho más fina que la otra.

Uno de los hallazgos más significativos de su trabajo es que este método directo es superior a los enfoques indirectos más antiguos que intentan proyectar los datos de una malla sobre la otra. Cuando el intercambio de productos químicos es fuerte, los métodos de proyección antiguos pierden información. Esencialmente, suavizan los detalles de la fractura, descartando las variaciones de escala fina que son cruciales para una predicción precisa. Los investigadores demostraron que su método directo conserva toda la información necesaria. En pruebas donde la transferencia de productos químicos era intensa, el método antiguo produjo errores un veintiuno por ciento mayores que su nuevo método. Además, cuando refinaron la malla de la fractura para que fuera aún más detallada manteniendo la malla de la roca fija, el método antiguo dejó de mejorar, estancándose en un estado de parálisis. Su nuevo método, sin embargo, continuó volviéndose más preciso, demostrando que captura la verdadera física de la interacción sin limitaciones artificiales.

Los investigadores también mostraron que este método funciona eficientemente en las computadoras. Escribieron un programa de software para probar su teoría, ejecutando simulaciones en una computadora de escritorio estándar. Compararon el tiempo que tomó configurar los cálculos para mallas coincidentes frente a mallas no coincidentes. Los resultados mostraron que el trabajo adicional requerido para manejar las piezas desajustadas fue insignificante. El tiempo para ensamblar las ecuaciones fue casi idéntico, lo que significa que los científicos ahora pueden usar mallas mucho más flexibles y realistas sin pagar un penalizador en tiempo de cómputo. Esta flexibilidad es vital para aplicaciones del mundo real, donde diferentes partes de una formación geológica podrían ser modeladas por diferentes equipos usando diferentes software, o donde la red de fracturas es tan compleja que una malla única y perfectamente alineada es imposible de crear.

El estudio también exploró cómo se comporta el método bajo condiciones extremas. Probaron escenarios donde las reacciones químicas eran muy fuertes y la transferencia de fluidos era rápida. En estos casos difíciles, el método se mantuvo robusto, convergiendo a una solución sin fallar. Incluso probaron la geometría de las fracturas, simulando grietas que se doblan o se intersectan entre sí. Aunque las pruebas principales utilizaron una grieta recta, la lógica subyacente de su método se mantuvo al aplicarse a estas formas más complejas, lo que sugiere que el enfoque podría extenderse a las redes enredadas que se encuentran en la naturaleza. Los investigadores fueron cuidadosos en señalar que su trabajo actual se centra en tipos específicos de rocas y geometrías de fractura, y que extender el método a redes tridimensionales completas, curvas e intersectadas es una tarea para el futuro.

Lo que hace que este trabajo sea particularmente importante es que elimina una barrera de larga data en la simulación de medios fracturados. Durante décadas, la necesidad de mallas coincidentes ha limitado cómo los científicos podían modelar estos sistemas, obligándolos a menudo a comprometer el detalle de la fractura o de la roca para que las matemáticas funcionaran. Al demostrar que las dos pueden acoplarse directamente y de forma conservativa sin una coincidencia perfecta, los investigadores han abierto la puerta a simulaciones más realistas y flexibles. Esto significa que, en el futuro, los modelos de contaminación de aguas subterráneas, almacenamiento de carbono o extracción de energía geotérmica pueden construirse con mayor libertad, permitiendo que la roca y las grietas se representen exactamente como son, en lugar de como deben ser para encajar en una restricción matemática rígida. El método proporciona una nueva y estable base para comprender cómo los fluidos y los productos químicos se mueven a través de la tierra fracturada, asegurando que la conexión entre la roca sólida y las grietas sea tan fluida en la computadora como lo es en el suelo.