A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology

Aquí se presenta un resumen técnico detallado del artículo "Un estudio matemático de un modelo de hielo marino elástico-viscoso-plástico con la reología de Kelvin-Voigt" de Boutros, Liu, Thomas y Titi.

1. Planteamiento del Problema y Motivación

El artículo aborda el análisis matemático del modelo de hielo marino Elástico-Viscoso-Plástico (EVP), una regularización numérica ampliamente utilizada del modelo clásico viscoso-plástico de Hibler. Si bien el modelo EVP es computacionalmente eficiente para simulaciones climáticas, su buena posición matemática rigurosa (existencia, unicidad y dependencia continua de los datos iniciales) sigue siendo un desafío, particularmente en presencia de términos de advección no lineales y sin recortes artificiales en los coeficientes de viscosidad.

Desafíos Clave:

• Reología No Newtoniana: El hielo marino se modela como un fluido no newtoniano con relaciones complejas entre tensión y deformación que involucran plasticidad.
• Término de Advección: La inclusión del término de advección no lineal ($u \cdot \nabla u$) en el balance de momento complica las estimaciones de energía y las pruebas de existencia global.
• Recorte de Viscosidad: Los modelos estándar a menudo requieren un "recorte" (parámetro de regularización $\epsilon$) en la tasa de deformación para evitar singularidades. Eliminar este recorte (tomando el límite $\epsilon \to 0$) es un problema abierto importante para modelos multidimensionales.
• Limitaciones Previas: En un trabajo anterior de los autores (arXiv:2505.03080), demostraron la buena posición para un modelo regularizado por Voigt donde la regularización se aplicaba a la relación constitutiva (evolución del tensor de tensiones). Sin embargo, ese enfoque requería alta regularidad de los datos iniciales y no podía manejar el término de advección al eliminar el recorte de viscosidad.

Objetivo:
Los autores proponen una nueva formulación donde la regularización de Kelvin-Voigt se aplica directamente al balance de momento (campo de velocidad) en lugar de a la evolución del tensor de tensiones. Su objetivo es demostrar:

1. Buena posición local para el modelo con el término de advección.
2. Buena posición global para el modelo sin el término de advección.
3. La capacidad de eliminar el recorte de viscosidad ($\epsilon \to 0$) en ambos casos.
4. Existencia de soluciones para datos iniciales menos regulares en comparación con estudios anteriores.

2. Formulación Matemática

El estudio considera el modelo EVP en un toro plano bidimensional $\mathbb{T}^2$. El sistema consiste en una ecuación de momento y una ecuación de evolución para el tensor de tensiones $\sigma$.

El Sistema Regularizado (Kelvin-Voigt en el Momento):
Los autores introducen un parámetro de regularización $\alpha > 0$ y modifican la ecuación de momento (1.6a/1.7a) para incluir un término $-\alpha^2 \partial_t \Delta u$:
$$\partial_t(u - \alpha^2 \Delta u) + (u \cdot \nabla)u = \nabla \cdot \sigma + \text{términos de forzamiento}$$
$$\frac{1}{E} \partial_t \sigma + \text{términos no lineales de tensión} = D(u)$$
donde $D(u)$ es el tensor de tasa de deformación.

Regularización de la Tasa de Deformación:
Para manejar la singularidad en la relación tensión-deformación, la magnitud de la tasa de deformación $D = |D(u)|$ se regulariza como $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$. Las pruebas implican construir soluciones para $\epsilon > 0$ y luego pasar al límite $\epsilon \to 0$.

3. Metodología

La estrategia de prueba se basa en un argumento de aproximación y compacidad de múltiples pasos:

1. Doble Regularización:

   • Tasa de Deformación: Utilizar $D_\epsilon$ para suavizar la no linealidad.
   • Momento: Introducir un término de regularización de orden superior $\beta^4 \partial_t \Delta^2 u$ (con $\beta > 0$) en el balance de momento para asegurar suficiente regularidad para el esquema de aproximación de Galerkin.
   • Datos Iniciales: Suavizar los datos iniciales $(u_0, \sigma_0)$ utilizando un mollificador estándar $\phi_\delta$.

2. Aproximación de Galerkin:

   • Construir soluciones aproximadas de dimensión finita $(u_N, \sigma_N)$ utilizando modos de Fourier.
   • Probar la existencia local para estos sistemas de EDO utilizando el teorema de Picard-Lindelöf.

3. Estimaciones A Priori:

   • Estimación $L^\infty$ para Tensión (Paso Crucial): Los autores derivan una nueva cota $L^\infty$ para el tensor de tensiones $\sigma$ (Lema 2.3). Al introducir una variable transformada $\tau = \sigma + \frac{P}{2}I_2$ y utilizar la estructura de amortiguamiento de la ecuación de evolución de tensiones, demuestran:
     $$\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$$
     Esta cota es independiente de los parámetros de regularización y es esencial para cerrar las estimaciones de orden superior.
   • Estimaciones de Energía: Derivar estimaciones $H^1$ y $H^2$ para la velocidad y estimaciones $H^1$ para la tensión. La cota $L^\infty$ sobre $\sigma$ permite controlar los términos no lineales que involucran al tensor de tensiones.

4. Paso de Límite:

   • Paso 1: Enviar $\beta, \delta \to 0$ para recuperar la solución del sistema con $D_\epsilon$ pero sin la regularización artificial $\beta$.
   • Paso 2: Enviar $\epsilon \to 0$ para eliminar el recorte de la tasa de deformación.
   • Compacidad: Utilizar el teorema de Banach-Alaoglu (para convergencia débil-*) y el lema de Aubin-Lions (para convergencia fuerte) para extraer subsucesiones convergentes.

5. Unicidad y Dependencia Continua:

   • Establecer la unicidad analizando la diferencia entre dos soluciones y aplicando la desigualdad de Grönwall, utilizando la cota $L^\infty$ sobre el tensor de tensiones.

4. Contribuciones y Resultados Clave

Teorema 1.1 (Buena Posición Local con Advección):

• Resultado: Demuestra la existencia de una solución fuerte única local en el tiempo para el modelo EVP de Kelvin-Voigt adveccionado.
• Regularidad: $u \in C([0, T]; H^2)$ y $\sigma \in C([0, T]; H^1) \cap C([0, T]; L^\infty)$.
• Significado: Este es el primer resultado riguroso para el modelo EVP multidimensional con advección donde se elimina el recorte de viscosidad ($\epsilon \to 0$). Las estimaciones anteriores no podían manejar el término de advección en el límite.

Teorema 1.2 (Buena Posición Global sin Advección):

• Resultado: Demuestra la existencia de una solución fuerte global en el tiempo única para el modelo EVP de Kelvin-Voigt cuando se omite el término de advección.
• Significado: La existencia global se mantiene para datos iniciales con menor regularidad que la requerida en el trabajo anterior de los autores (donde la regularización estaba en el tensor de tensiones). Específicamente, la velocidad se regulariza, permitiendo datos iniciales $H^2$ en lugar de requerir mayor suavidad.

Manejo del Recorte de Viscosidad:

• El artículo justifica con éxito el límite $\epsilon \to 0$ (eliminando el recorte artificial en los coeficientes de viscosidad) en presencia de advección (Teorema 1.1), un problema que permanece abierto para el modelo estándar de Hibler.

Interpretación Física:

• Se argumenta que la regularización de Kelvin-Voigt en el balance de momento es físicamente más intuitiva que regularizar la relación constitutiva, ya que modela directamente el comportamiento elástico-viscoso del hielo en el balance de fuerzas.

5. Significado e Impacto

• Rigor Matemático: El artículo proporciona la primera teoría rigurosa de buena posición para el modelo EVP con advección y sin recortes de viscosidad en múltiples dimensiones.
• Superación de Problemas Abiertos: Resuelve un problema abierto importante relacionado con el límite $\epsilon \to 0$ en presencia de advección, que previamente era intratable para los modelos de Hibler y Voigt-EVP.
• Requisitos de Regularidad Mejorados: Al desplazar la regularización a la ecuación de momento, los autores reducen la regularidad requerida de los datos iniciales, haciendo que el modelo sea aplicable a una clase más amplia de escenarios físicos.
• Nuevas Herramientas Analíticas: La derivación de la estimación $L^\infty$ para el tensor de tensiones (Lema 2.3) es una contribución técnica novedosa que permite cerrar las estimaciones de energía en presencia de fuertes no linealidades.
• Modelado Climático: Estos resultados proporcionan una base matemática sólida para los esquemas numéricos utilizados en simulaciones climáticas a gran escala, asegurando que los modelos sean matemáticamente sólidos y estables bajo condiciones realistas.

En resumen, este trabajo avanza la teoría matemática de la dinámica del hielo marino al establecer resultados robustos de buena posición para una regularización motivada físicamente del modelo EVP, superando obstáculos analíticos significativos relacionados con la no linealidad, la advección y los límites de viscosidad singular.