Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

En el vasto y gélido norte, el suelo no es simplemente tierra; es un cimiento estructural que sostiene las ciudades más grandes del mundo, oleoductos y carreteras. Durante siglos, este permafrost se ha mantenido sólido, actuando su estado congelado como un lecho rocoso rígido para la construcción humana. Sin embargo, a medida que el clima se calienta, este suelo congelado está comenzando a descongelarse. Cuando el hielo dentro del suelo se derrite, el terreno pierde su resistencia, lo que provoca hundimientos desiguales y la agrietación de los edificios. Predecir exactamente qué tan profundo llegará este deshielo es una cuestión de seguridad e ingeniería, pero el suelo mismo presenta una complicación masiva. El suelo real no es un bloque uniforme de hielo; es una mezcla caótica de capas, lentes y bolsas de diferentes materiales, con propiedades que cambian en distancias de apenas unos pocos centímetros. Para simular el derretimiento de un suelo tan complejo con una precisión perfecta, una computadora necesitaría rastrear cada una de las capas, una tarea que requiere más potencia de cálculo de la que está disponible actualmente para pronósticos a gran escala y largo plazo.

Para resolver esto, los investigadores recurren a una técnica llamada homogeneización numérica. Imagine intentar predecir el flujo de agua a través de una esponja. En lugar de modelar cada diminuto agujero y fibra, usted trata a la esponja como un único material con una capacidad "promedio" para dejar pasar el agua. Esta es la esencia de la homogeneización: reemplazar una realidad desordenada y detallada con un modelo simplificado y promediado que sea lo suficientemente rápido para ejecutarse en una computadora. Sin embargo, cuando el suelo se está derritiendo, este promedio se convierte en una trampa. El límite entre el suelo congelado y el descongelado —el "frente"— no se detiene ordenadamente en el borde del bloque de cálculo de la computadora. Corta directamente a través de su centro. Esto crea un dilema: si un bloque de suelo está medio congelado y medio descongelado, ¿cómo debe decidir la computadora sus propiedades? ¿Debería actuar más como el hielo, más como el agua, o algo intermedio? La elección realizada aquí determina si el modelo predice que un edificio se hundirá o se mantendrá firme.

En un estudio reciente, Timur Zakharov, de la Universidad Federal del Noreste, abordó este dilema específico probando tres formas distintas de manejar estos bloques de suelo "mixtos". La investigación se centró en un edificio climatizado situado sobre suelo congelado, un escenario donde el calor de la estructura derrite lentamente el suelo debajo de él. El equipo creó un gemelo digital de este entorno, utilizando un patrón aleatorio y realista de conductividad del suelo para imitar el caos natural del terreno. Luego, ejecutaron la simulación en una malla gruesa, donde cada bloque de suelo era mucho más grande que las capas reales de tierra y hielo, obligando al frente de cambio de fase a atravesar el medio de estos bloques. Compararon tres métodos distintos para decidir qué sucede dentro de esos bloques seccionados.

El primer método, conocido como el enfoque de tabla, precalcula una lista de propiedades para cada estado posible de un bloque de suelo y busca la respuesta basándose en la temperatura promedio del bloque. El segundo método, un esquema nodal, decide el estado del suelo en las esquinas de los bloques y promedia esas decisiones. El tercer método, llamado suavizado h (h-smoothing), calcula la fracción exacta del bloque que se ha derretido basándose en un gradiente de temperatura lineal, midiendo esencialmente la porción del bloque que se ha convertido en agua. Los investigadores realizaron estas simulaciones tanto en dos como en tres dimensiones, comparando los resultados con una solución de "referencia" altamente detallada que rastreaba cada capa de suelo.

Los hallazgos revelaron un intercambio sorprendente que desafía la forma en que los ingenieros podrían elegir sus herramientas. Cuando el objetivo es predecir la temperatura del suelo, el método de suavizado h resultó ser el más preciso. Rastreó el flujo de calor con la mayor precisión, especialmente en las mallas gruesas donde los otros métodos tuvieron dificultades. Sin embargo, cuando los investigadores observaron la posición del frente de derretimiento —el límite real entre el suelo congelado y el descongelado— la historia cambió por completo. El método de suavizado h, aunque bueno en temperatura, sobreestimó consistentemente qué tan profundo llegaría el deshielo. Predijo que el suelo se derretiría más de lo que realmente lo haría.

La razón de esta discrepancia reside en las leyes de la conservación de la energía. El método de suavizado h, en su esfuerzo por ser preciso sobre la forma de la zona de derretimiento, añadió inadvertidamente una pequeña cantidad de calor extra al sistema con cada paso de tiempo. A lo largo de un año, este excedente diminuto se acumuló, causando que el suelo simulado se derritiera más de lo que el suelo real lo haría. En contraste, los otros dos métodos, particularmente el enfoque de tabla y el esquema nodal, fueron "conservadores". No añadieron ni perdieron energía y, como resultado, predijeron la profundidad del deshielo con mucha más exactitud, incluso si sus lecturas de temperatura eran ligeramente menos precisas. El estudio demostró que si un investigador necesita saber la temperatura exacta del suelo, el suavizado h es la opción superior. Pero si el objetivo es saber qué tan profundo se hundirá el cimiento debido al deshielo, los métodos conservadores son mucho más fiables.

La investigación también exploró qué sucede cuando la temperatura a la que el suelo se derrite no es la misma en todas partes, sino que varía a través de las diferentes capas. En este escenario más complejo, el método basado en tablas, que utiliza un registro detallado de propiedades, superó a los métodos más simples de dos tensores en la predicción de la posición del frente de deshielo. Esto sugiere que, para condiciones de suelo altamente variables, es necesario un mapa de estados del suelo precalculado más detallado. En última instancia, el estudio proporciona una guía clara para ingenieros y científicos: no existe un único método "mejor". La elección depende enteramente de qué pregunta se esté haciendo. Si la prioridad es la precisión térmica, un método gana; si la prioridad es el límite físico del deshielo, otro toma la delantera. Al aclarar estas diferencias, este trabajo ayuda a asegurar que los modelos utilizados para proteger los edificios en el Ártico no sean solo rápidos, sino también confiables.