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

以下は、Boutros、Liu、Thomas、Titi による論文「A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology（ケルビン・フォイトレオロジーを伴う弾性・粘性・塑性海氷モデルの数学的研究）」の詳細な技術的要約である。

1. 問題提起と動機

本論文は、古典的な粘性・塑性 Hibler モデルの広範に用いられる数値的正則化である弾性・粘性・塑性（EVP）海氷モデルの数学的解析を取り扱っている。EVP モデルは気候シミュレーションにおいて計算効率が優れているが、非線形移流項の存在下および粘性係数に対する人工的なカットオフなしでは、その厳密な数学的適切性（存在、一意性、初期データに対する連続依存性）は依然として困難な課題である。

主要な課題：

• 非ニュートンレオロジー： 海氷は、塑性を含む複雑な応力 - ひずみ関係を持つ非ニュートン流体としてモデル化される。
• 移流項： 運動量バランスに非線形移流項（$u \cdot \nabla u$）を含めることは、エネルギー評価および大域的存在証明を複雑にする。
• 粘性カットオフ： 標準的なモデルでは、特異性を回避するためにひずみ速度に「カットオフ」（正則化パラメータ $\epsilon$）を課すことがよく行われる。このカットオフを除去すること（極限 $\epsilon \to 0$ を取ることは）は、多次元モデルにおける主要な未解決問題である。
• 以前の限界： 著者らの以前の研究（arXiv:2505.03080）では、正則化が構成則（応力テンソルの進化）に適用される Voigt 正則化モデルの適切性を証明した。しかし、そのアプローチは初期データの高度な正則性を要求し、粘性カットオフを除去する際に移流項を扱うことができなかった。

目的：
著者らは、ケルビン・フォイト正則化を応力テンソルの進化ではなく、直接運動量バランス（速度場）に適用する新しい定式化を提案する。彼らは以下のことを証明することを目的としている：

1. 移流項を含むモデルの局所適切性。
2. 移流項を含まないモデルの大域適切性。
3. 両ケースにおいて粘性カットオフ（$\epsilon \to 0$）を除去する能力。
4. 以前の研究と比較して、より低い正則性を持つ初期データに対する解の存在。

2. 数学的定式化

本研究は、2 次元平坦トーラス $\mathbb{T}^2$ 上の EVP モデルを検討する。この系は、運動量方程式と応力テンソル $\sigma$ の進化方程式から構成される。

正則化された系（運動量におけるケルビン・フォイト）：
著者らは正則化パラメータ $\alpha > 0$ を導入し、運動量方程式（1.6a/1.7a）に項 $-\alpha^2 \partial_t \Delta u$ を含むように修正する：
$$\partial_t(u - \alpha^2 \Delta u) + (u \cdot \nabla)u = \nabla \cdot \sigma + \text{強制項}$$
$$\frac{1}{E} \partial_t \sigma + \text{非線形応力項} = D(u)$$
ここで、$D(u)$ はひずみ速度テンソルである。

ひずみ速度の正則化：
応力 - ひずみ関係における特異性を処理するために、ひずみ速度の大きさ $D = |D(u)|$ は $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$ として正則化される。証明には、$\epsilon > 0$ に対して解を構成し、その後極限 $\epsilon \to 0$ を通じる過程が含まれる。

3. 手法

証明戦略は、多段階の近似とコンパクト性の議論に依存している：

1. 二重正則化：

   • ひずみ速度： 非線形性を滑らかにするために $D_\epsilon$ を使用する。
   • 運動量： ガレルキン近似スキームに十分な正則性を確保するために、運動量バランスに高次正則化項 $\beta^4 \partial_t \Delta^2 u$（$\beta > 0$）を導入する。
   • 初期データ： 標準的なモルフィファイア $\phi_\delta$ を用いて初期データ $(u_0, \sigma_0)$ を滑らかにする。

2. ガレルキン近似：

   • フーリエモードを用いて有限次元近似解 $(u_N, \sigma_N)$ を構成する。
   • ピカール・リンデレフの定理を用いて、これらの常微分方程式系の局所存在を証明する。

3. 事前評価：

   • 応力に対する $L^\infty$ 評価（決定的なステップ）： 著者らは応力テンソル $\sigma$ に対する新しい $L^\infty$ 評価（補題 2.3）を導出する。変換された変数 $\tau = \sigma + \frac{P}{2}I_2$ を導入し、応力進化方程式の減衰構造を利用することで、以下を証明する：
     $$\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$$
     この評価は正則化パラメータに依存せず、高次評価を閉じるために不可欠である。
   • エネルギー評価： 速度に対する $H^1$ および $H^2$ 評価、および応力に対する $H^1$ 評価を導出する。$\sigma$ に対する $L^\infty$ 評価により、応力テンソルを含む非線形項の制御が可能となる。

4. 極限通過：

   • ステップ 1： $\beta, \delta \to 0$ を行って、人工的な $\beta$ 正則化なしで $D_\epsilon$ を含む系の解を取り戻す。
   • ステップ 2： $\epsilon \to 0$ を行って、ひずみ速度のカットオフを除去する。
   • コンパクト性： 弱*収束のための Banach-Alaoglu の定理と、強収束のための Aubin-Lions の補題を用いて、収束する部分列を抽出する。

5. 一意性と連続依存性：

   • 2 つの解の差を解析し、応力テンソルに対する $L^\infty$ 評価を利用してグロンワールの不等式を適用することで、一意性を確立する。

4. 主要な貢献と結果

定理 1.1（移流を含む局所適切性）：

• 結果： 移流を含むケルビン・フォイト EVP モデルに対する、時間局所的な強解の一意存在を証明する。
• 正則性： $u \in C([0, T]; H^2)$ および $\sigma \in C([0, T]; H^1) \cap C([0, T]; L^\infty)$。
• 意義： これは、粘性カットオフが除去された（$\epsilon \to 0$）移流を伴う多次元 EVP モデルに対する最初の厳密な結果である。以前の評価は極限において移流項を扱うことができなかった。

定理 1.2（移流なしの大域適切性）：

• 結果： 移流項を省略した場合のケルビン・フォイト EVP モデルに対する、一意な時間大域的強解の存在を証明する。
• 意義： 大域的存在は、著者らの以前の研究（正則化が応力テンソルに行われていた場合）で要求されたものよりも低い正則性を持つ初期データに対して成り立つ。具体的には、速度が正則化されるため、より高い滑らかさを要求するのではなく、$H^2$ 初期データで十分である。

粘性カットオフの処理：

• 本論文は、移流の存在下（定理 1.1）において極限 $\epsilon \to 0$（粘性係数に対する人工的カットオフの除去）を正当化することに成功している。これは、標準的な Hibler モデルにおいては依然として未解決の問題である。

物理的解釈：

• ケルビン・フォイト正則化を運動量バランスに適用することは、構成則を正則化することよりも物理的に直感的であると論じられている。なぜなら、それは力バランスにおける氷の弾性・粘性挙動を直接モデル化するためである。

5. 意義と影響

• 数学的厳密性： 本論文は、多次元において移流を伴い粘性カットオフなしの EVP モデルに対する、最初の厳密な適切性理論を提供する。
• 未解決問題の克服： 以前は Hibler モデルおよび Voigt-EVP モデルに対して扱いが困難であった、移流の存在下における極限 $\epsilon \to 0$ に関する主要な未解決問題を解決する。
• 改善された正則性要件： 正則化を運動量方程式に移すことで、著者らは初期データに要求される正則性を低下させ、モデルをより広範な物理的シナリオに適用可能にした。
• 新しい解析ツール： 応力テンソルに対する $L^\infty$ 評価（補題 2.3）の導出は、強い非線形性の存在下でエネルギー評価を閉じることを可能にする新規の技術的貢献である。
• 気候モデリング： これらの結果は、大規模気候シミュレーションで使用される数値スキームに堅固な数学的基盤を提供し、現実的な条件下でもモデルが数学的に健全で安定であることを保証する。

要約すると、この研究は、非線形性、移流、および特異な粘性極限に関連する重要な解析的障壁を克服し、EVP モデルの物理的に動機付けられた正則化に対する堅牢な適切性結果を確立することにより、海氷力学の数学的理論を前進させた。