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

以下是 Boutros、Liu、Thomas 和 Titi 所著论文《具有开尔文 - 沃伊特流变学的弹性 - 粘性 - 塑性海冰模型的数学研究》的详细技术总结。

1. 问题陈述与动机

本文针对**弹性 - 粘性 - 塑性（EVP）**海冰模型的数学分析展开研究，该模型是经典粘性 - 塑性 Hibber 模型的一种广泛使用的数值正则化形式。尽管 EVP 模型在气候模拟中具有计算效率，但其严格的数学适定性（存在性、唯一性以及对初始数据的连续依赖性）仍然具有挑战性，特别是在存在非线性平流项且未对粘性系数进行人为截断的情况下。

关键挑战：

• 非牛顿流变学： 海冰被建模为具有复杂应力 - 应变关系（涉及塑性）的非牛顿流体。
• 平流项： 动量平衡方程中包含非线性平流项（$u \cdot \nabla u$），这增加了能量估计和全局存在性证明的复杂性。
• 粘性截断： 标准模型通常需要在应变率上设置“截断”（正则化参数 $\epsilon$）以避免奇点。在多维模型中移除该截断（取极限 $\epsilon \to 0$）是一个主要的未解决问题。
• 先前局限性： 作者之前的工作（arXiv:2505.03080）证明了针对本构关系（应力张量演化）施加正则化的沃伊特正则化模型的适定性。然而，该方法要求初始数据具有高正则性，并且在移除粘性截断时无法处理平流项。

目标：
作者提出了一种新公式，将开尔文 - 沃伊特正则化直接应用于动量平衡（速度场），而非应力张量演化。他们旨在证明：

1. 包含平流项的模型的局部适定性。
2. 不包含平流项的模型的全局适定性。
3. 在两种情况下均能移除粘性截断（$\epsilon \to 0$）。
4. 相较于以往研究，能够处理正则性较低的初始数据的解的存在性。

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$，恢复具有 $D_\epsilon$ 但无人为 $\beta$-正则化的系统解。
   • 步骤 2： 令 $\epsilon \to 0$，移除应变率截断。
   • 紧性： 利用巴拿赫 - 阿劳格卢定理（用于弱*收敛）和奥宾 - 利翁斯引理（用于强收敛）提取收敛子序列。

5. 唯一性与连续依赖性：

   • 通过分析两个解之间的差异并应用格朗沃尔不等式，利用应力张量的 $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)$。
• 意义： 这是针对多维 EVP 模型在移除粘性截断（$\epsilon \to 0$）且包含平流情况下的首个严格结果。先前的估计无法在极限情况下处理平流项。

定理 1.2（不含平流的全局适定性）：

• 结果： 证明了当省略平流项时，开尔文 - 沃伊特 EVP 模型存在唯一的全局时间强解。
• 意义： 全局存在性对初始数据的要求正则性更低，低于作者先前工作（当时正则化施加于应力张量）的要求。具体而言，速度被正则化，允许使用 $H^2$ 初始数据，而无需更高的光滑性。

处理粘性截断：

• 本文成功论证了在存在平流的情况下取极限 $\epsilon \to 0$（移除粘性系数的人为截断）的合理性（定理 1.1），而这一问题对于标准 Hibber 模型而言仍然是未解决的。

物理解释：

• 论证表明，在动量平衡中应用开尔文 - 沃伊特正则化比在本构关系中应用正则化更符合物理直觉，因为它直接在力平衡中模拟了冰的弹性 - 粘性行为。

5. 意义与影响

• 数学严谨性： 本文提供了首个针对多维 EVP 模型（含平流且无粘性截断）的严格适定性理论。
• 克服未解决问题： 它解决了在存在平流情况下关于极限 $\epsilon \to 0$ 的一个主要未解决问题，该问题此前对于 Hibber 和沃伊特 -EVP 模型而言是不可处理的。
• 改进的正则性要求： 通过将正则化移至动量方程，作者降低了对初始数据正则性的要求，使该模型适用于更广泛的物理场景。
• 新的分析工具： 推导应力张量的 $L^\infty$ 估计（引理 2.3）是一项新颖的技术贡献，使得在强非线性存在的情况下能够闭合能量估计。
• 气候建模： 这些结果为大规模气候模拟中使用的数值方案提供了坚实的数学基础，确保模型在现实条件下具有数学合理性和稳定性。

总之，这项工作通过为 EVP 模型的物理动机正则化建立稳健的适定性结果，推进了海冰动力学的数学理论，克服了与非线性、平流和奇异粘性极限相关的重大分析障碍。