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

1. Probleemstelling en Motivatie

Het artikel behandelt de wiskundige analyse van het Elastic-Viscous-Plastic (EVP) zee-ijsmodel, een veelgebruikte numerieke regularisatie van het klassieke visco-plastische Hibler-model. Hoewel het EVP-model rekenkundig efficiënt is voor klimaatsimulaties, blijft de strikte wiskundige goedgesteldheid (existentie, uniciteit en continue afhankelijkheid van beginvoorwaarden) een uitdaging, met name in aanwezigheid van niet-lineaire advectietermen en zonder kunstmatige afsnijdingen van viscositeitscoëfficiënten.

Belangrijkste Uitdagingen:

• Niet-Newtoniaanse Reologie: Zee-ijs wordt gemodelleerd als een niet-Newtoniaanse vloeistof met complexe spanning-rekrelaties die plasticiteit omvatten.
• Advectieterm: De opname van de niet-lineaire advectieterm ($u \cdot \nabla u$) in de impulsbalans bemoeilijkt energie-estimaten en bewijzen van globale existentie.
• Viscositeitsafsnijding: Standaardmodellen vereisen vaak een "afsnijding" (regularisatieparameter $\epsilon$) op de rekrate om singulariteiten te vermijden. Het verwijderen van deze afsnijding (door de limiet $\epsilon \to 0$ te nemen) is een groot open probleem voor meerdimensionale modellen.
• Eerdere Beperkingen: In een eerder werk van de auteurs (arXiv:2505.03080) bewezen zij goedgesteldheid voor een door Voigt-geregulariseerd model waarbij de regularisatie werd toegepast op de constitutieve relatie (evolutie van de spannings tensor). Deze aanpak vereiste echter hoge regulariteit van de beginvoorwaarden en kon de advectieterm niet hanteren bij het verwijderen van de viscositeitsafsnijding.

Doel:
De auteurs stellen een nieuwe formulering voor waarbij de Kelvin-Voigt-regularisatie direct wordt toegepast op de impulsbalans (snelheidsveld) in plaats van op de evolutie van de spannings tensor. Zij beogen het volgende te bewijzen:

1. Lokale goedgesteldheid voor het model met de advectieterm.
2. Globale goedgesteldheid voor het model zonder de advectieterm.
3. De mogelijkheid om de viscositeitsafsnijding te verwijderen ($\epsilon \to 0$) in beide gevallen.
4. Existentie van oplossingen voor minder reguliere beginvoorwaarden in vergelijking met eerdere studies.

2. Wiskundige Formulering

De studie beschouwt het EVP-model op een 2D platte torus $\mathbb{T}^2$. Het systeem bestaat uit een impulsvergelijking en een evolutievergelijking voor de spannings tensor $\sigma$.

Het Geregulariseerde Systeem (Kelvin-Voigt in Impuls):
De auteurs introduceren een regularisatieparameter $\alpha > 0$ en modificeren de impulsvergelijking (1.6a/1.7a) om een term $-\alpha^2 \partial_t \Delta u$ op te nemen:
$$\partial_t(u - \alpha^2 \Delta u) + (u \cdot \nabla)u = \nabla \cdot \sigma + \text{forcerende termen}$$
$$\frac{1}{E} \partial_t \sigma + \text{niet-lineaire spannings termen} = D(u)$$
waarbij $D(u)$ de rekrate tensor is.

Regularisatie van Rekrate:
Om de singulariteit in de spanning-rekrelatie te hanteren, wordt de grootte van de rekrate $D = |D(u)|$ geregulariseerd als $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$. De bewijzen omvatten het construeren van oplossingen voor $\epsilon > 0$ en vervolgens het nemen van de limiet $\epsilon \to 0$.

3. Methodologie

De bewijsstrategie rust op een meerstapsbenadering en compactheidsargument:

1. Dubbele Regularisatie:

   • Rekrate: Gebruik $D_\epsilon$ om de niet-lineariteit te gladstrijken.
   • Impuls: Introduceer een regularisateterm van hogere orde $\beta^4 \partial_t \Delta^2 u$ (met $\beta > 0$) in de impulsbalans om voldoende regulariteit te garanderen voor het Galerkin-benaderingsschema.
   • Beginvoorwaarden: Mollificeer beginvoorwaarden $(u_0, \sigma_0)$ met een standaard mollifier $\phi_\delta$.

2. Galerkin-benadering:

   • Construeer eindig-dimensionale benaderende oplossingen $(u_N, \sigma_N)$ met behulp van Fourier-modes.
   • Bewijs lokale existentie voor deze ODE-systemen met behulp van de stelling van Picard-Lindelöf.

3. A priori Estimaten:

   • $L^\infty$-estimatie voor Spanning (Cruciaal Stap): De auteurs leiden een nieuwe $L^\infty$-grens af voor de spannings tensor $\sigma$ (Lemma 2.3). Door een getransformeerde variabele $\tau = \sigma + \frac{P}{2}I_2$ in te voeren en de dempingsstructuur van de spannings-evolutievergelijking te benutten, bewijzen zij:
     $$\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$$
     Deze grens is onafhankelijk van de regularisatieparameters en is essentieel om hogere-orde estimaten te sluiten.
   • Energie-estimaten: Leid $H^1$- en $H^2$-estimaten af voor snelheid en $H^1$-estimaten voor spanning. De $L^\infty$-grens op $\sigma$ maakt het mogelijk niet-lineaire termen die de spannings tensor betreffen, te beheersen.

4. Limietovergang:

   • Stap 1: Stuur $\beta, \delta \to 0$ om de oplossing te herstellen voor het systeem met $D_\epsilon$ maar zonder de kunstmatige $\beta$-regularisatie.
   • Stap 2: Stuur $\epsilon \to 0$ om de rekrate-afsnijding te verwijderen.
   • Compactheid: Gebruik de stelling van Banach-Alaoglu (voor zwak-* convergentie) en het lemma van Aubin-Lions (voor sterke convergentie) om convergente deelrijen te extraheren.

5. Uniciteit en Continue Afhankelijkheid:

   • Stel uniciteit vast door de verschillen tussen twee oplossingen te analyseren en de ongelijkheid van Grönwall toe te passen, gebruikmakend van de $L^\infty$-grens op de spannings tensor.

4. Belangrijkste Bijdragen en Resultaten

Stelling 1.1 (Lokale Goedgesteldheid met Advectie):

• Resultaat: Bewijst de existentie van een unieke lokaal-in-de-tijd sterke oplossing voor het advectieve Kelvin-Voigt EVP-model.
• Regulariteit: $u \in C([0, T]; H^2)$ en $\sigma \in C([0, T]; H^1) \cap C([0, T]; L^\infty)$.
• Betekenis: Dit is het eerste strikte resultaat voor het meerdimensionale EVP-model met advectie waarbij de viscositeitsafsnijding is verwijderd ($\epsilon \to 0$). Eerdere estimaten faalden om de advectieterm in de limiet te hanteren.

Stelling 1.2 (Globale Goedgesteldheid zonder Advectie):

• Resultaat: Bewijst de existentie van een unieke globaal-in-de-tijd sterke oplossing voor het Kelvin-Voigt EVP-model wanneer de advectieterm wordt weggelaten.
• Betekenis: De globale existentie geldt voor beginvoorwaarden met lagere regulariteit dan vereist in het eerdere werk van de auteurs (waar de regularisatie op de spannings tensor lag). Specifiek wordt de snelheid geregulariseerd, wat $H^2$ beginvoorwaarden toestaat in plaats van hogere gladheid te vereisen.

Hanteer van Viscositeitsafsnijding:

• Het artikel rechtvaardigt succesvol de limiet $\epsilon \to 0$ (verwijdering van de kunstmatige afsnijding van viscositeitscoëfficiënten) in aanwezigheid van advectie (Stelling 1.1), een probleem dat voor het standaard Hibler-model open blijft.

Fysische Interpretatie:

• De Kelvin-Voigt-regularisatie in de impulsbalans wordt fysisch intuïtiever geacht dan het regulariseren van de constitutieve relatie, omdat deze direct het elastisch-viskeuze gedrag van het ijs in de krachtenbalans modelleert.

5. Betekenis en Impact

• Wiskundige Strenge: Het artikel biedt de eerste strikte theorie van goedgesteldheid voor het EVP-model met advectie en zonder viscositeitsafsnijdingen in meerdere dimensies.
• Oplossen van Open Problemen: Het lost een groot open probleem op betreffende de limiet $\epsilon \to 0$ in aanwezigheid van advectie, wat voorheen onberekenbaar was voor het Hibler- en Voigt-EVP-model.
• Verbeterde Regulariteitseisen: Door de regularisatie naar de impulsvergelijking te verschuiven, verlagen de auteurs de vereiste regulariteit van beginvoorwaarden, waardoor het model toepasbaar wordt op een bredere klasse van fysische scenario's.
• Nieuwe Analytische Hulpmiddelen: De afleiding van de $L^\infty$-estimatie voor de spannings tensor (Lemma 2.3) is een nieuwe technische bijdrage die het sluiten van energie-estimaten in aanwezigheid van sterke niet-lineariteiten mogelijk maakt.
• Klimaatmodellering: Deze resultaten bieden een solide wiskundige basis voor de numerieke schema's die worden gebruikt in grootschalige klimaatsimulaties, waarbij wordt gegarandeerd dat de modellen wiskundig sound en stabiel zijn onder realistische omstandigheden.

Kortom, dit werk vordert de wiskundige theorie van zee-ijsdynamica door robuuste resultaten van goedgesteldheid vast te stellen voor een fysisch gemotiveerde regularisatie van het EVP-model, waarbij aanzienlijke analytische hindernissen met betrekking tot niet-lineariteit, advectie en singulariteitsviscositeitslimieten worden overwonnen.