Diffuse Adiabatic Flows in Thermally Coupled Grounded Shallow Ice Sheets: Modelling and Analysis

Resumen Técnico: Flujos Adiabáticos Difusos en Capas de Hielo Someras Ancladas Térmicamente Acopladas

Planteamiento del Problema
Este artículo aborda el modelado matemático de la evolución simultánea del espesor y la temperatura interna de una capa de hielo somera anclada. El desafío central radica en el acoplamiento de dos fenómenos físicos distintos: la evolución mecánica de la superficie del hielo (regida por el balance de masa y la dinámica de flujo) y la evolución termodinámica de la temperatura interna (regida por la advección-difusión).

El problema se formula bajo restricciones físicas específicas:

1. Restricción Geométrica: El espesor del hielo debe permanecer no negativo, descansando sobre un lecho plano y congelado. Esto transforma la evolución del espesor en un problema de obstáculo dependiente del tiempo.
2. Restricción Térmica: La temperatura interna está restringida para permanecer por debajo del punto de fusión por presión (régimen congelado), lo que evita la generación de calor latente y el cambio de fase. Esto introduce un segundo problema de obstáculo para el campo de temperatura.
3. Dominio Móvil: La ecuación de la temperatura se plantea sobre un dominio cuya extensión vertical está determinada por el espesor de hielo desconocido, lo que crea un problema de frontera móvil.

El autor señala que estudios rigurosos previos (p. ej., [63]) trataron la evolución del espesor por separado o asumieron un lecho plano sin acoplamiento térmico. Este trabajo extiende el análisis al sistema totalmente acoplado donde la suavidad del hielo (viscosidad) depende de su temperatura interna, la cual, a su vez, depende del espesor del hielo.

Metodología
El artículo emplea un marco analítico riguroso que combina la derivación formal, desigualdades variacionales y argumentos de compacidad.

1. Derivación Formal y Reducción:

   • El autor deriva las ecuaciones gobernantes para la elevación de la superficie del hielo y la temperatura interna basándose en la Aproximación de Hielo Somero y la ley de flujo de Glen.
   • Para manejar el dominio móvil y la degeneración del $p$-Laplaciano cerca del margen del hielo (donde el espesor se anula), se introduce un cambio de variables ($h = u^{(p-1)/2p}$), siguiendo las estrategias de Raviart.
   • Se adopta una estrategia de "interfaz difusa" para fijar el dominio computacional. Se construye un campo de temperatura auxiliar para extender la solución desde la región del hielo hacia un cilindro fijo $\Omega \times (0, L)$. Esto implica una función de corte cerca de la superficie del hielo, modelando efectivamente el sistema como si estuviera cubierto por una capa de escombros donde la temperatura se anula cerca de la superficie. Esto permite plantear el problema en un dominio fijo, aunque con una aproximación cerca de la frontera.

2. Penalización y Discretización:

   • Los problemas de obstáculo acoplados se reformulan como un sistema "penalizado" ($P_\ell$) utilizando términos de penalización para imponer las restricciones ($u \ge 0$ y $0 \le \tilde{\Theta} \le \Theta_m$).
   • La existencia de soluciones para el problema penalizado se establece mediante un esquema de paso de tiempo semidiscreto (diferencias finitas). La prueba se basa en el teorema de Minty-Browder, utilizando la monotonicidad y coercitividad de los operadores asociados.

3. Estimaciones A Priori y Compacidad:

   • Se deriva un conjunto exhaustivo de estimaciones a priori para las soluciones discretas, independientes del paso de tiempo y del parámetro de penalización. Estas estimaciones acotan las soluciones en los espacios de Sobolev y Bochner apropiados.
   • La compacidad se logra mediante el teorema de Aubin-Lions-Simon y una versión discreta del lema de Dubinskii. Se presta especial atención al manejo de las no linealidades derivadas del $p$-Laplaciano y la dependencia exponencial de la suavidad del hielo respecto a la temperatura.

4. Paso al Límite:

   • La existencia de soluciones para el problema penalizado se utiliza para pasar al límite cuando el paso de tiempo tiende a cero, recuperando una solución para el problema continuo penalizado.
   • Finalmente, se envía el parámetro de penalización $\ell$ a cero. El autor establece la existencia de soluciones "muy débiles" para el sistema acoplado original. La prueba evita supuestos de regularidad adicional sobre el límite débil que estaban presentes en trabajos anteriores, alineando el concepto de solución con el marco de las desigualdades variacionales evolucionarias gobernadas por operadores pseudo-monótonos.

Contribuciones Clave y Novedades
El artículo identifica tres novedades principales:

1. Reducción a un Dominio Fijo: La introducción de técnicas de reducción (Sección 2) transforma el modelo original, que acopla espesor y temperatura en un dominio móvil, en un modelo equivalente en un dominio fijo. Esto implica una aproximación donde la evolución de la temperatura se anula cerca de la superficie del hielo, caracterizando capas de hielo con una cubierta de escombros (flujos adiabáticos).
2. Existencia de Soluciones Muy Débiles: El artículo establece la existencia de soluciones para el problema de valor inicial-frontera donde tanto el espesor del hielo como la temperatura interna están sujetos a restricciones físicas (problemas de obstáculo). El concepto de solución es "muy débil", debido a la baja regularidad del problema que impide las definiciones de solución estándar.
3. Marco Analítico: El trabajo proporciona un nuevo marco analítico para una clase de problemas de obstáculo acoplados planteados en dominios móviles, extendiéndose más allá de aplicaciones glaciológicas específicas.

Resultados Principales
El resultado central es el Teorema 8.1, que demuestra la existencia de soluciones muy débiles para el problema de valor inicial-frontera que acopla la evolución del espesor de la capa de hielo somera con su temperatura interna, asumiendo un lecho plano y congelado.

Aspectos clave del resultado incluyen:

• La eliminación del "supuesto de disminución extra" (4.6) requerido en análisis previos (p. ej., [63]).
• La derivación de un modelo límite abstracto que no depende de supuestos de regularidad adicional para el límite débil.
• El concepto de solución se alinea con marcos no estándar para desigualdades variacionales evolucionarias que involucran operadores pseudo-monótonos.

Significancia y Reivindicaciones
El artículo afirma extender el estudio matemático riguroso de los modelos de capas de hielo iniciado en [63] al caso de acoplamiento térmico. El autor enfatiza que el acoplamiento introduce desafíos significativos debido al dominio móvil y las no linealidades asociadas a la frontera libre.

La significancia del trabajo se enmarca de la siguiente manera:

• Modelado: Ofrece un modelo termodinámico novedoso que acopla rigurosamente la elevación superficial y la temperatura interna bajo restricciones físicas.
• Analítica: Supera la dificultad de definir soluciones en dominios móviles empleando un enfoque de interfaz difusa y un cambio de variables para manejar la degeneración.
• Robustez: Al prescindir de supuestos de regularidad adicional, el resultado proporciona una teoría de existencia más robusta, consistente con la baja regularidad inherente a los modelos de hielo somero.

El artículo concluye que el modelo propuesto y la prueba de existencia representan un avance en el análisis matemático de las capas de hielo someras ancladas, particularmente en contextos donde el acoplamiento térmico y las restricciones físicas (lecho congelado) son críticos. El autor señala explícitamente que la naturaleza "adiabática" del título se refiere a la elección específica del modelo térmico justificada por la presencia de capas de escombros, más que a una afirmación general sobre todas las dinámicas de las capas de hielo.