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

1. Formulação do Problema e Motivação

O artigo aborda a análise matemática do modelo de gelo marinho Elástico-Viscoso-Plástico (EVP), uma regularização numérica amplamente utilizada do clássico modelo viscoso-plástico de Hibler. Embora o modelo EVP seja computacionalmente eficiente para simulações climáticas, sua boa colocação matemática rigorosa (existência, unicidade e dependência contínua em relação aos dados iniciais) permanece desafiadora, particularmente na presença de termos de advecção não lineares e sem cortes artificiais nos coeficientes de viscosidade.

Desafios Principais:

• Reologia Não Newtoniana: O gelo marinho é modelado como um fluido não newtoniano com relações complexas entre tensão e deformação envolvendo plasticidade.
• Termo de Advecção: A inclusão do termo de advecção não linear ($u \cdot \nabla u$) no balanço de momento complica as estimativas de energia e as provas de existência global.
• Corte de Viscosidade: Modelos padrão frequentemente exigem um "corte" (parâmetro de regularização $\epsilon$) na taxa de deformação para evitar singularidades. Remover esse corte (tomando o limite $\epsilon \to 0$) é um problema em aberto importante para modelos multidimensionais.
• Limitações Anteriores: Em um trabalho anterior dos autores (arXiv:2505.03080), eles provaram a boa colocação para um modelo regularizado por Voigt, onde a regularização era aplicada à relação constitutiva (evolução do tensor de tensões). No entanto, essa abordagem exigia alta regularidade dos dados iniciais e não conseguia lidar com o termo de advecção ao remover o corte de viscosidade.

Objetivo:
Os autores propõem uma nova formulação onde a regularização de Kelvin-Voigt é aplicada diretamente ao balanço de momento (campo de velocidade) em vez da evolução do tensor de tensões. Eles visam provar:

1. Boa colocação local para o modelo com o termo de advecção.
2. Boa colocação global para o modelo sem o termo de advecção.
3. A capacidade de remover o corte de viscosidade ($\epsilon \to 0$) em ambos os casos.
4. Existência de soluções para dados iniciais menos regulares em comparação com estudos anteriores.

2. Formulação Matemática

O estudo considera o modelo EVP em um toro plano 2D $\mathbb{T}^2$. O sistema consiste em uma equação de momento e uma equação de evolução para o tensor de tensões $\sigma$.

O Sistema Regularizado (Kelvin-Voigt no Momento):
Os autores introduzem um parâmetro de regularização $\alpha > 0$ e modificam a equação de momento (1.6a/1.7a) para incluir um termo $-\alpha^2 \partial_t \Delta u$:
$$\partial_t(u - \alpha^2 \Delta u) + (u \cdot \nabla)u = \nabla \cdot \sigma + \text{termos de forçagem}$$
$$\frac{1}{E} \partial_t \sigma + \text{termos não lineares de tensão} = D(u)$$
onde $D(u)$ é o tensor de taxa de deformação.

Regularização da Taxa de Deformação:
Para lidar com a singularidade na relação tensão-deformação, a magnitude da taxa de deformação $D = |D(u)|$ é regularizada como $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$. As provas envolvem a construção de soluções para $\epsilon > 0$ e, em seguida, a passagem ao limite $\epsilon \to 0$.

3. Metodologia

A estratégia de prova baseia-se em uma aproximação multi-etapa e um argumento de compacidade:

1. Dupla Regularização:

   • Taxa de Deformação: Uso de $D_\epsilon$ para suavizar a não linearidade.
   • Momento: Introdução de um termo de regularização de ordem superior $\beta^4 \partial_t \Delta^2 u$ (com $\beta > 0$) no balanço de momento para garantir regularidade suficiente para o esquema de aproximação de Galerkin.
   • Dados Iniciais: Suavização dos dados iniciais $(u_0, \sigma_0)$ usando um mollificador padrão $\phi_\delta$.

2. Aproximação de Galerkin:

   • Construção de soluções aproximadas de dimensão finita $(u_N, \sigma_N)$ usando modos de Fourier.
   • Prova de existência local para esses sistemas de EDOs usando o teorema de Picard-Lindelöf.

3. Estimativas A Priori:

   • Estimativa $L^\infty$ para Tensão (Etapa Crucial): Os autores derivam uma nova limitação $L^\infty$ para o tensor de tensões $\sigma$ (Lema 2.3). Ao introduzir uma variável transformada $\tau = \sigma + \frac{P}{2}I_2$ e utilizar a estrutura de amortecimento da equação de evolução da tensão, eles provam:
     $$\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$$
     Essa limitação é independente dos parâmetros de regularização e é essencial para fechar as estimativas de ordem superior.
   • Estimativas de Energia: Derivação de estimativas $H^1$ e $H^2$ para a velocidade e estimativas $H^1$ para a tensão. A limitação $L^\infty$ sobre $\sigma$ permite o controle de termos não lineares envolvendo o tensor de tensões.

4. Passagem ao Limite:

   • Etapa 1: Enviar $\beta, \delta \to 0$ para recuperar a solução para o sistema com $D_\epsilon$, mas sem a regularização artificial $\beta$.
   • Etapa 2: Enviar $\epsilon \to 0$ para remover o corte da taxa de deformação.
   • Compacidade: Uso do teorema de Banach-Alaoglu (para convergência fraca-*) e do lema de Aubin-Lions (para convergência forte) para extrair subsequências convergentes.

5. Unicidade e Dependência Contínua:

   • Estabelecimento da unicidade analisando a diferença entre duas soluções e aplicando a desigualdade de Grönwall, utilizando a limitação $L^\infty$ sobre o tensor de tensões.

4. Contribuições e Resultados Principais

Teorema 1.1 (Boa Colocação Local com Advecção):

• Resultado: Prova a existência de uma solução forte única local no tempo para o modelo EVP de Kelvin-Voigt advectivo.
• Regularidade: $u \in C([0, T]; H^2)$ e $\sigma \in C([0, T]; H^1) \cap C([0, T]; L^\infty)$.
• Significado: Este é o primeiro resultado rigoroso para o modelo EVP multidimensional com advecção onde o corte de viscosidade é removido ($\epsilon \to 0$). Estimativas anteriores falharam em lidar com o termo de advecção no limite.

Teorema 1.2 (Boa Colocação Global sem Advecção):

• Resultado: Prova a existência de uma solução forte única global no tempo para o modelo EVP de Kelvin-Voigt quando o termo de advecção é omitido.
• Significado: A existência global vale para dados iniciais com menor regularidade do que a exigida no trabalho anterior dos autores (onde a regularização estava no tensor de tensões). Especificamente, a velocidade é regularizada, permitindo dados iniciais $H^2$ em vez de exigir maior suavidade.

Tratamento do Corte de Viscosidade:

• O artigo justifica com sucesso o limite $\epsilon \to 0$ (removendo o corte artificial nos coeficientes de viscosidade) na presença de advecção (Teorema 1.1), um problema que permanece em aberto para o modelo Hibler padrão.

Interpretação Física:

• Argumenta-se que a regularização de Kelvin-Voigt no balanço de momento é fisicamente mais intuitiva do que regularizar a relação constitutiva, pois modela diretamente o comportamento elástico-viscoso do gelo no balanço de forças.

5. Significado e Impacto

• Rigor Matemático: O artigo fornece a primeira teoria rigorosa de boa colocação para o modelo EVP com advecção e sem cortes de viscosidade em múltiplas dimensões.
• Superação de Problemas Abertos: Resolve um grande problema em aberto relativo ao limite $\epsilon \to 0$ na presença de advecção, que anteriormente era intratável para os modelos Hibler e Voigt-EVP.
• Requisitos de Regularidade Melhorados: Ao deslocar a regularização para a equação de momento, os autores reduzem a regularidade exigida dos dados iniciais, tornando o modelo aplicável a uma classe mais ampla de cenários físicos.
• Novas Ferramentas Analíticas: A derivação da estimativa $L^\infty$ para o tensor de tensões (Lema 2.3) é uma contribuição técnica novel que permite o fechamento das estimativas de energia na presença de fortes não linearidades.
• Modelagem Climática: Estes resultados fornecem uma base matemática sólida para os esquemas numéricos usados em simulações climáticas em grande escala, garantindo que os modelos sejam matematicamente sólidos e estáveis sob condições realistas.

Em resumo, este trabalho avança a teoria matemática da dinâmica do gelo marinho ao estabelecer resultados robustos de boa colocação para uma regularização motivada fisicamente do modelo EVP, superando obstáculos analíticos significativos relacionados à não linearidade, advecção e limites de viscosidade singular.