Implementación de un Método Euleriano con Esquema TVD para la Resolución de Problemas de Transporte Reactivo en Medios Porosos con Alta Heterogeneidad

Resumen

La resolución de problema de transporte en medios porosos con alta heterogeneidad en la conductividad es comúnmente realizada a través de métodos Lagrangianos. Estos métodos son adecuados para problemas puramente advectivos, sin embargo, cuando se consideran efectos difusivos el número de partículas necesarias para describir el sistema puede ser muy alto, lo que aumenta la demanda computacional. También presentan sobrecostos de cómputo importantes a la hora de representar procesos que implican la interacción con fases no móviles. Por otro lado, los métodos Eulerianos presentan ventajas para resolver problemas que implican la interacción entre fases móviles y no móviles, como la disolución o la precipitación de minerales, y cuando se combinan con esquemas de alta resolución, como los TVD, pueden ser aplicados para resolver problemas advectivos-difusivos. En este trabajo se presenta una implementación en GPU de un método Euleriano con un esquema TVD para la resolución de problemas de transporte reactivo en medios porosos con alta heterogeneidad, considerando una reacción de disolución/precipitación en equilibro entre dos especies disueltas y una especie inmóvil la cual afecta la porosidad.