Resumen de Un modelo numérico para la simulación de transporte de calor y liberación de materia en un almacenamiento profundo de residuos radiactivos
El objetivo de esta Tesis Doctoral es el desarrollo de un modelo numérico que permita simular fenómenos de liberación y transporte de isótopos radiactivos y de calor a través de la geosfera, procedentes de un almacenamiento geológico profundo de residuos radiactivos de alta actividad, Para ello, se consideran los diferentes fenómenos fisicos que intervienen en el problema: flujo de agua subterránea, desintegración radiactiva, disolución de radionúclidos en agua subterranea, generación de calor, transporte a través de la geosfera de materia y de calor. Algunos de esos fenómenos están relacionados entre si, lo que permite construir un modelo acoplado, a partir del cual crear un código escrito en el lenguaje FORTRAN.
Los modelos de flujo y transporte se desarrollan en dos dimensiones espaciales, y se resuelven mediante un método de volúmenes finitos para la integración espacial, mientras que la integración temporal se realiza mediante un theta-metodo.
Asimismo, se resuelve la ecuación de convección-difusión unidimensional mediante dos técnicas de volúmenes finitos, una con interpolación lineal y otra mediante interpolación cuadrática con descentraje en el sentido de la velocidad de convencción. También se realiza un estudio de consistencia y estabilidad de ambos métodos comparando sus zonas de estabilidad y los errores que se cometen. Para el estudio de estabilidad de emplea la técnica de von Neumann, basada en armónicos de Fourier, que, si bien es una técnica clásica en diferencias finitas, en esta Tesis se aplica a dos esquemas en volúmenes finitos.
