生成时间: 2026年04月30日 05:52:28

重点关注论文数量: 3

了解领域论文数量: 3

不相关论文数量: 12

重点关注论文（完整分析）

1. 具记忆效应的非局部守恒律系统及其零记忆极限

Systems of Nonlocal Conservation Laws with Memory and Their Zero Retention Limit

作者: Aekta Aggarwal, Ganesh Vaidya

类别: math.NA, math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26877v1

摘要翻译: 我们研究一类具记忆效应的非局部守恒律系统的熵解：其中对流通量关于空间与时间变量均与一个核函数做卷积。该构型建模了通量对解在其空间邻域内（空间非局部性）以及对其历史状态（时间非局部性）的依赖，从而纳入记忆效应。此外，借助一个收敛的有限体积逼近格式，我们讨论了熵解的存在性，并建立了此类熵解的唯一性。

Comment: 无

链接: http://arxiv.org/abs/2604.26877v1

1. 研究对象和背景

本文研究一类同时具有空间和时间非局部性的守恒律系统，其核心方程为：

\partial_t U^k + \partial_x \left(f^k(U^k)\nu^k((\Theta \ast U)^k)\right) = 0, \quad (t, x) \in Q_T := (0, T) \times \mathbb{R},

U^k(0, x) = U^k_0(x), \quad x \in \mathbb{R},

其中 $U = (U^k)_{k \in \mathbb{N}} : [0, \infty) \times \mathbb{R} \to \mathbb{R}^N$ 是未知函数， $\mathbb{N} = \{1, \ldots, N\}$ ， $$T$$ 为最终时间。卷积项定义为：

(\Theta^{i,k} \ast U^i)(t, x) := \int_0^t \int_{\mathbb{R}} U^i(\tau, \xi)\mu^{j,k}(x - \xi)\Gamma^{i,k}(t - \tau) d\tau d\xi, \quad i \in \mathbb{N},

且 $\Theta^{j,k}(t, x) := \mu^{j,k}(x)\Gamma^{j,k}(t)$ ，其中 $\mu$ 和 $\Gamma$ 是光滑的 $N \times N$ 矩阵。假设条件为：(H1) $f^k \in \mathrm{Lip}(\mathbb{R})$ 且 $$f^k(0) = 0 = f^k(1)$$ ；(H2) $\nu^k \in (C^2 \cap BV \cap W^{2,\infty})(\mathbb{R}^N, \mathbb{R})$ ；(H3) $\mu^{j,k} \in C^2(\mathbb{R}) \cap W^{2,\infty}(\mathbb{R})$ ， $\Gamma^{j,k} \in C^2([0, \infty); \mathbb{R}^+) \cap W^{2,\infty}([0, \infty); \mathbb{R}^+)$ 。

该模型的物理意义在于：在任意时空点 $$(t_0, x_0)$$ ，通量不仅依赖于瞬时密度 $$U^k(t_0, x_0)$$ ，还依赖于 (i) 由空间核 $\mu^{j,k}$ 确定的 $$x_0$$ 邻域内状态 $U^j(t_0, \cdot)$ 的加权平均；(ii) 由时间核 $\Gamma^{j,k}$ 调制的过去状态 $U^j(\tau, x_0)$ （ $\tau < t_0$ ）的累积影响。这使得第 $$k$$ 个方程在 $$(t_0, x_0)$$ 处的特征速度反映了系统所有分量之间的历史依赖性交互，而非仅由局部瞬时数据决定。

此类模型在诸多实际应用中出现：弹性介质中的阻尼现象 [22,23]、粘弹性 [32]、分数时间行为模型 [54]、煤层中涉及吸附和扩散的气体传输模型 [24, 55]、地下和次河带区域的传输模型 [41,43,44] 等。然而，这类时间非局部守恒律的适定性理论在很大程度上尚不完善。主要难点在于：时间卷积项通常缺乏空间光滑性，导致通量作为空间变量的函数不连续；时间卷积破坏了通常意义上的半群性质，产生新的数学对象，与局部或纯空间非局部守恒律中出现的对象不同。因此，经典的标准守恒律方法不直接适用。

相比之下，仅空间非局部的守恒律（即 $\Gamma \equiv 1$ ）：

\partial_t U^k + \partial_x \left(f^k(U^k)\nu^k((\mu \ast U)^k)\right) = 0, \quad (t, x) \in Q_T,

(\mu^{j,k} \ast U^j)(t, x) := \int_{\mathbb{R}} U^j(t, \xi)\mu^{j,k}(x - \xi) d\xi,

近十年来已得到广泛研究 [1–12,17,19–21,29,30,37–39,42,48,50]，原因有二：(1) 广泛应用于人群 [1,29]、交通 [4,5,12,17,38]、沉降模型 [11]、激光技术 [30]、颗粒材料动力学 [8] 等领域；(2) 其分析/数学特性与局部非线性守恒律接近。

然而，同时考虑空间和时间非局部性的系统研究甚少。特别地，当 $$N > 1$$ 时，熵解的存在性和唯一性问题尚未解决。本文试图填补以下空白：(1) 建立一般形式的带记忆效应非局部守恒律系统的熵解适定性理论；(2) 分析当时间卷积核支持集缩小时（记忆效应减弱）的渐近行为；(3) 研究有限体积方法在记忆到无记忆极限过渡中的渐近兼容性。值得注意的是，对于 $$N > 1$$ 的系统，到局部对应模型的过渡（即同时让空间和时间卷积核支持集趋于零）极其困难，因为对应的局部系统 $$(1.8)$$ 在 $$N > 2$$ 时仅有有限的适定性结果（如对足够小的 $$BV$$ 初值 [13, 45] 或 $$N = 2$$ 时对足够小的 $L^\infty$ 初值 [14]）。这种显著差异使得 $$N > 1$$ 时的极限问题成为悬而未决的开放问题。

2. 主要定理或主要结果

在假设 (H1)-(H3) 下，本文证明了以下主要结果：

唯一性定理（Theorem 1）：对任意 $$T > 0$$ ，设 $$U, V$$ 为系统 (1.3)-(1.4) 初值问题的熵解，对应初值分别为 $$U_0, V_0$$ ，则 $\|U(T, \cdot) - V(T, \cdot)\|_{(L^1(\mathbb{R}))^N} \leq \|U_0 - V_0\|_{(L^1(\mathbb{R}))^N} + C \sum_{k \in \mathbb{N}} \int_0^T \|U^k(\tau, \cdot) - V^k(\tau, \cdot)\|_{L^1(\mathbb{R})} d\tau,$ 其中 $C = NT|f|_{(\mathrm{Lip}(\mathbb{R}))^N} |\nu|_{(\mathrm{Lip}(\mathbb{R}^N))_N} \|\Theta\|_{(L^\infty(Q_T))^{N^2}} |U|_{(L^\infty_t BV_x)^N}$ $$

NT|f|{(\mathrm{Lip}(\mathbb{R}))^N} |U_0|{(L^1(\mathbb{R}))^N} |\nabla \nu|{(L^\infty(\mathbb{R}^N))^{N^2}} |\Gamma|{(L^\infty(\mathbb{R}^+))^{N^2}} |\mu’|_{(L^\infty(\mathbb{R}))^{N^2}} $$$ $$$
NT|f|{(\mathrm{Lip}(\mathbb{R}))^N} |U_0|{(L^1(\mathbb{R}))^N} |V|{(L^1(Q_T))^N} |\Gamma|{(L^\infty(\mathbb{R}^+))^{N^2}} |\mu’|{(L^\infty(\mathbb{R}))^{N^2}} \times |Hess \nu|{(L^\infty(\mathbb{R}^N))^{N^3}} |\Theta|_{(L^\infty(Q_T))^{N^2}}. $$

有限体积方法收敛性：在 CFL 条件
$\lambda \leq \frac{\min(1, 4 - 6\beta, 6\beta)}{1 + 6|f|_{(\mathrm{Lip}(\mathbb{R}))^N} \|\nu\|_{(L^\infty(\mathbb{R}^N))^N}}$
下，所提出的有限体积格式满足：(a) 单调性；(b) 不变区域原理；(c) 守恒性；(d) $$L^1$$ 有界性；(e) 空间有界变差估计；(f) 时间估计；(g) 离散熵条件；(h) 卷积项边界估计。特别地，空间有界变差估计为：
$\sum_{i \in \mathbb{Z}} |U^{k,n}_{i+1} - U^{k,n}_i| \leq \left(\exp(C^k_7 t) \sum_{i \in \mathbb{Z}} |U^{k,0}_{i+1} - U^{k,0}_i| + \frac{\exp(C^k_7 t) - 1}{C^k_7} C^k_8 \right),$
其中 $C^k_7 = C^k_5 |f^k|_{\mathrm{Lip}(\mathbb{R})} \|\nabla \nu^k\|_{(L^\infty(\mathbb{R}^N))^N}$ ， $C^k_8 = C^k_6 |f^k|_{\mathrm{Lip}(\mathbb{R})} \|\nabla \nu^k\|_{(L^\infty(\mathbb{R}^N))^N} \|U^k_0\|_{L^1(\mathbb{R})} + 2(C^k_5)^2|f^k|_{\mathrm{Lip}(\mathbb{R})}|\nabla \nu^k|_{\mathrm{Lip}(\mathbb{R}^N)} \|U^k_0\|_{L^1(\mathbb{R})}$ ，且 $C^k_5 = \|U^k_0\|_{L^1(\mathbb{R})} \|\dot{\mu}^{j,k}\|_{L^\infty(\mathbb{R})} \|\Gamma^{j,k}\|_{L^1(\mathbb{R}^+)}$ 。这些估计的关键在于仅依赖于 $\|\Gamma^{j,k}\|_{L^1(\mathbb{R})}$ ，而与 $\|\Gamma^{j,k}\|_{L^\infty(\mathbb{R})}$ 无关。
渐近行为分析：当时间卷积核支持集半径 $\delta \to 0$ 时，空间-时间非局部守恒律的唯一熵解以 $O(\sqrt{\delta})$ 的速率收敛到纯空间非局部守恒律的熵解。
渐近兼容性：所提出的有限体积近似在记忆到无记忆极限过渡中是渐近兼容的，误差估计为 $O(\sqrt{\delta} + \sqrt{\Delta x})$ 。具体而言，数值解 $U^{\Delta}$ 与极限解 $$U^0$$ 满足：
$\|U^{\Delta}(\cdot, T) - U^0(\cdot, T)\|_{L^1(\mathbb{R})} \leq C(\sqrt{\delta} + \sqrt{\Delta x}),$
其中 $$C$$ 为与 $\delta$ 和 $\Delta x$ 无关的常数。
熵解存在性：结合唯一性定理和有限体积方法的收敛性，建立了系统 (1.3)-(1.4) 熵解的适定性理论。

3. 研究方法、关键技术和核心工具

证明中的核心障碍在于时间卷积项缺乏空间光滑性，导致通量作为空间变量的函数不连续，且时间卷积破坏了通常的半群性质，使经典守恒律方法不直接适用。本文的创新点在于系统地将纯空间非局部守恒律的理论技术适应到同时包含空间和时间非局部效应的更广泛通量函数类，特别利用了对流通量关于空间变量 $$x$$ 的正则性。

关键估计与技术框架包括：

唯一性证明中的精细分解：将差值 $\|U^k - V^k\|_{L^1}$ 分解为两项：
$I^k_1 := \int_0^T \int_{\mathbb{R}} |U^k_x(t, x)| \left|\nu^k((\Theta \ast U)^k(t, x)) - \nu^k((\Theta \ast V)^k(t, x))\right| dx d\tau,$
$I^k_2 := \int_0^T \int_{\mathbb{R}} |U^k(t, x)| \left|\partial_x(\nu^k((\Theta \ast U)^k(t, x))) - \partial_x(\nu^k((\Theta \ast V)^k(t, x)))\right| dx d\tau.$
对 $$I^k_1$$ 的估计关键在于：
$\left|\nu^k((\Theta \ast U)^k) - \nu^k((\Theta \ast V)^k)\right| \leq |\nu^k|_{\mathrm{Lip}(\mathbb{R}^N)} \|\Theta\|_{(L^\infty(Q_T))^{N^2}} \sum_{k \in \mathbb{N}} \int_0^t \|U^k(\tau, \cdot) - V^k(\tau, \cdot)\|_{L^1(\mathbb{R})} d\tau,$
这利用了 $\nu^k$ 的 Lipschitz 连续性和卷积核的 $L^\infty$ 有界性。对 $$I^k_2$$ 的估计更为复杂，需分解为：
$\nabla \nu^k((\Theta \ast U)^k) \cdot \partial_x((\Theta \ast U)^k) - \nabla \nu^k((\Theta \ast V)^k) \cdot \partial_x((\Theta \ast V)^k) = I^k_{21} + I^k_{22},$
其中 $I^k_{21}$ 涉及 $$U - V$$ 的 $$L^1$$ 差，而 $I^k_{22}$ 涉及 $\nabla \nu^k$ 的 Lipschitz 连续性，最终得到：
$I^k_{21} \leq C_2 \sum_{k \in \mathbb{N}} \int_0^T \|U^k(\tau, \cdot) - V^k(\tau, \cdot)\|_{L^1(\mathbb{R})} d\tau,$
$I^k_{22} \leq C_3 \sum_{k \in \mathbb{N}} \int_0^T \|U^k(\tau, \cdot) - V^k(\tau, \cdot)\|_{L^1(\mathbb{R})} d\tau,$
其中 $C_2 = \|\nabla \nu\|_{(L^\infty(\mathbb{R}^N))^{N^2}} \|\Gamma\|_{(L^\infty(\mathbb{R}^+))^{N^2}} \|\mu'\|_{(L^\infty(\mathbb{R}))^{N^2}}$ ， $C_3 = \|V\|_{(L^1(Q_T))^N} \|\Gamma\|_{(L^\infty(\mathbb{R}^+))^{N^2}} \|\mu'\|_{(L^\infty(\mathbb{R}))^{N^2}} \|Hess \nu\|_{(L^\infty(\mathbb{R}^N))^{N^3}} \|\Theta\|_{(L^\infty(Q_T))^{N^2}}$ 。
有限体积方法中的关键估计：卷积项满足
$|c^{j,k,n}_{i+1/2} - c^{j,k,n}_{i-1/2}| \leq C^k_5 \Delta x, \quad |c^{j,k,n}_{i+3/2} - 2c^{j,k,n}_{i+1/2} + c^{j,k,n}_{i-1/2}| \leq C^k_6 \Delta x^2,$
其中 $C^k_5 = \|U^k_0\|_{L^1(\mathbb{R})} \|\dot{\mu}^{j,k}\|_{L^\infty(\mathbb{R})} \|\Gamma^{j,k}\|_{L^1(\mathbb{R}^+)}$ ， $C^k_6 = 2 \|U^k_0\|_{L^1(\mathbb{R})} \|\ddot{\mu}^{j,k}\|_{L^\infty(\mathbb{R})} \|\Gamma^{j,k}\|_{L^1(\mathbb{R}^+)}$ 。这些估计的关键在于仅依赖于 $\|\Gamma^{j,k}\|_{L^1(\mathbb{R})}$ 而非 $\|\Gamma^{j,k}\|_{L^\infty(\mathbb{R})}$ ，使得在渐近分析中可处理 $\delta \to 0$ 的情形。
渐近分析中的核心技巧：利用空间有界变差估计 $|U|_{(L^\infty_t BV_x)^N} \leq C$ （与 $\|\Gamma\|_{L^\infty}$ 无关），结合对时间卷积核支持集缩小的精细分析。具体而言，当 $\Gamma^{j,k}$ 的支持集为 $[0, \delta]$ 时，通过估计
$\left| \int_0^t \int_{\mathbb{R}} U(\tau, \xi) \mu(x - \xi) \Gamma(t - \tau) d\xi d\tau - \int_{\mathbb{R}} U(t, \xi) \mu(x - \xi) d\xi \right| \leq C \sqrt{\delta},$
利用 $$U$$ 的时间连续模 $\gamma(U, \tau) = \max_{k \in \mathbb{N}} \sup_{|t_1 - t_2| \leq \tau} \|U^k(t_1, \cdot) - U^k(t_2, \cdot)\|_{L^1(\mathbb{R})}$ 和空间有界变差性，得到 $\gamma(U, \tau) \leq C \tau$ ，从而导出 $O(\sqrt{\delta})$ 的收敛率。
渐近兼容性证明：将数值解与极限解的误差分解为三角不等式中的三项，其中数值方法的离散误差为 $O(\sqrt{\Delta x})$ ，而渐近过渡的误差为 $O(\sqrt{\delta})$ ，最终得到组合误差 $O(\sqrt{\delta} + \sqrt{\Delta x})$ 。特别地， $\Gamma^{j,k,s}$ 采用单元平均而非单元中心点值的离散化对证明渐近兼容性至关重要（见 Remark 1）。

4. 与之前工作的比较

与既有文献相比，本文在多个维度上推进了非局部守恒律理论：

适定性理论的扩展：

现有文献 [1–12,17,19–21,29,30,37–39,42,48,50] 主要关注仅空间非局部的守恒律（ $\Gamma \equiv 1$ ），通常假设空间核 $\mu$ 具有凸性等几何限制。
本文首次为同时具有空间和时间非局部性的系统建立完整的适定性理论，且不施加任何几何限制（如核函数的凸性），适用于更一般的 $\mu$ 和 $\Gamma$ 。
对于 $$N > 1$$ 的系统，本文解决了熵解的存在性和唯一性问题，而此前 [36] 仅针对 $$N = 1$$ 和特定 $\Theta$ 及线性 $$f$$ 使用松弛方案。

渐近行为分析的创新：

本文首次系统分析当时间卷积核支持集缩小时（ $\delta \to 0$ ）的渐近行为，证明了 $O(\sqrt{\delta})$ 的收敛率，而此前文献未涉及此类分析。
现有工作 [15,16,25–28,31,38,40,49] 仅针对 $$N = 1$$ 、特定 $\mu$ 和线性 $$f$$ 研究从纯空间非局部到局部极限的过渡，而本文关注的是从空间-时间非局部到仅空间非局部的过渡（问题 Q1）。
本文证明了有限体积方法在记忆到无记忆极限中的渐近兼容性，这是文献中首次对带记忆的非局部守恒律进行此类分析。

技术方法的突破：

本文的关键创新在于获得与 $\|\Gamma\|_{L^\infty}$ 无关的空间有界变差估计，这使渐近分析成为可能。相比之下，[36] 的方法依赖于特定 $\Theta$ 的结构，无法推广到一般情形。
在唯一性证明中，本文通过精细分解 $$I^k_1$$ 和 $$I^k_2$$ ，克服了时间卷积导致的非局部耦合困难，而 [18, 22, 23, 32–34, 51–53] 仅针对特定方程利用结构特性证明适定性。
卷积项的边界估计（Lemma 2(h)）仅依赖于 $\|\Gamma\|_{L^1}$ 而非 $\|\Gamma\|_{L^\infty}$ ，这是处理 $\delta \to 0$ 时 $\|\Gamma\|_{L^\infty} \to \infty$ 情形的关键，而现有文献未考虑此类估计。

适用范围与局限：

本文成功解决 (Q1) 和 (Q4)，即从有记忆模型到无记忆（仅空间非局部）对应模型的过渡，但未解决 (Q2)（到仅时间非局部设定的过渡），因为这会产生通量局部但在空间不连续的系统，其适定性在很大程度上开放。
对于 (Q3)（到局部对应模型的过渡），本文指出这在 $$N > 1$$ 时极为困难：局部系统 (1.8) 在 $$N > 2$$ 时仅有对小 $$BV$$ 初值 [13, 45] 或 $$N = 2$$ 时对小 $L^\infty$ 初值 [14] 的存在性结果，而本文系统 (1.3)-(1.4) 对任意 $$BV$$ 初值已证适定 [3,4]。这种差距使得 $$N > 1$$ 时的极限问题成为悬而未决的难题。
本文的收敛率 $O(\sqrt{\delta})$ 可能非最优，但现有文献未提供更优估计。数值实验表明该速率是紧的，因为当 $\delta$ 减小时，误差以 $\sqrt{\delta}$ 的速率衰减。

与数值分析文献的对比：

文献 [35,46] 研究了 $$N = 1$$ 、 $\Theta = \mu$ 且 $$f$$ 线性时的渐近兼容性，但未考虑时间非局部性。
本文首次将渐近兼容性分析扩展到带记忆效应的系统，证明了有限体积方法在 $\delta \to 0$ 时的一致收敛性，而 [36] 仅通过松弛方案间接处理特定情形。
本文的误差估计 $O(\sqrt{\delta} + \sqrt{\Delta x})$ 揭示了数值参数 $\Delta x$ 与物理参数 $\delta$ 的耦合效应，这是非局部问题数值分析中的新现象。

代价与推进的平衡：

本文放宽了对卷积核几何性质的假设，但代价是收敛率 $O(\sqrt{\delta})$ 低于某些特定情形下的 $O(\delta)$ （如当 $\Gamma$ 为 Dirac 序列且 $$U$$ 光滑时）。
通过避免对 $\Gamma$ 的特殊假设（如单调性），本文获得了更广泛适用的理论，但无法利用特定核的结构获得更优估计。
本文的分析框架可推广到更一般的非局部算子，但处理 $$N > 1$$ 时到局部极限的过渡仍需突破性进展，因为局部系统的内在困难（如共振现象）尚未解决。

2. 相对论等离子体–真空界面的稳定性与存在性

Stability and existence of relativistic plasma–vacuum interfaces

作者: Paolo Secchi, Yuri Trakhinin, Tao Wang

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26864v1

摘要翻译: 我们研究二维与三维空间中相对论等离子体–真空界面的自由边界问题。等离子体流动由理想相对论磁流体力学方程组描述，而真空中的电磁场满足麦克斯韦方程组。等离子体与真空的磁场均切于界面，且界面随等离子体流动而运动。由此得到一个具有可变重数特征边界的非线性、多维双曲型自由边界问题。我们识别出一个定量稳定性条件，并在变系数线性化问题于各向异性索博列夫空间中满足能量估计的意义下，建立了三维相对论等离子体–真空界面的线性稳定性。在估计切向导数时，我们利用一种内蕴抵消效应，将边界项转化为瞬时积分；进而将涉及空间导数的估计与涉及时间导数的估计分离，使得该瞬时积分在稳定性条件下主要被瞬时能量所吸收。此外，在初始界面上任意点处等离子体与真空磁场不同时为零的假设下，我们证明了二维空间中非线性问题的局部时间存在性与唯一性。证明结合了线性化问题的可解性与适定估计（tame estimate），以及一种适当的修正纳什–莫泽迭代法。特别地，为建立其可解性，我们将二维线性化问题约化为关于界面函数的输运方程，以及一个具有极大非负边界条件的双曲型边值问题。

Comment: 52 pages

链接: http://arxiv.org/abs/2604.26864v1

1. 研究对象和背景

本文研究空间维度 $$d = 2, 3$$ 下相对论等离子体-真空界面的自由边界问题。等离子体区域 $\Omega^{+}(t) \subset \mathbb{R}^d$ 由理想相对论磁流体动力学（RMHD）方程描述：

\partial_t(\rho\Gamma) + \nabla\cdot (\rho\Gamma v) = 0,

\partial_t\left(\rho f\Gamma^2 + \epsilon^2|H|^2 -\epsilon^2 q\right) + \nabla\cdot\left(\rho f\Gamma^2 v + \epsilon^2|H|^2 v -\epsilon^2(v \cdot H)H\right) = 0,

\partial_t\left(\rho f\Gamma^2 v + \epsilon^2|H|^2 v -\epsilon^2(v \cdot H)H\right) -\nabla\cdot\left(\Gamma^{-2}H \otimes H\right) + \nabla q + \nabla\cdot\left((\rho f\Gamma^2 + \epsilon^2|H|^2)v \otimes v -\epsilon^2(v \cdot H)(H \otimes v + v \otimes H)\right) = 0,

\partial_t H -\nabla\times (v \times H) = 0,

以及散度约束 $\nabla\cdot H = 0$ 。其中 $\rho$ 为粒子数密度， $v = (v_1, \ldots, v_d)$ 为坐标速度， $H = (H_1, \ldots, H_d)$ 为磁场， $$q$$ 为总压力， $\epsilon^{-1}$ 为光速，洛伦兹因子 $\Gamma = 1/\sqrt{1 - \epsilon^2|v|^2}$ ，流体指标 $f = 1 + \epsilon^2(e + p/\rho)$ ，总压力 $q = p + |H|^2/(2\Gamma^2) + \epsilon^2(v \cdot H)^2/2$ ， $$e$$ 为内能， $$p$$ 为压力，满足吉布斯关系 $\vartheta dS = de - p\rho^{-2} d\rho$ ， $$S$$ 为熵， $\vartheta > 0$ 为绝对温度。

真空区域 $\Omega^{-}(t)$ 由麦克斯韦方程描述：

\epsilon\partial_t h + \nabla\times e = 0, \quad \epsilon\partial_t e - \nabla\times h = 0,

以及散度约束 $\nabla\cdot h = 0$ 和 $$d = 3$$ 时 $\nabla\cdot e = 0$ 。界面 $\Sigma(t) = \Omega^{+}(t) \cap \Omega^{-}(t)$ 随等离子体流动移动，边界条件为：

H \cdot n = 0, \quad h \cdot n = 0, \quad V = v \cdot n, \quad 2q = |h|^2 - |e|^2, \quad \epsilon V h = n \times e \quad \text{on } \Sigma(t),

其中 $$n$$ 为指向 $\Omega^{+}(t)$ 的单位法向量， $$V$$ 为界面法向速度。

该问题的重要性在于：等离子体-真空界面在磁约束等离子体模型中起核心作用。非相对论MHD适用于低速等离子体，而RMHD适用于相对论流。然而，相对论情形下自由边界为特征边界且具有可变重数（边界矩阵奇异且秩非常数），导致不同界面区域需施加不同数量边界条件，取决于法向速度符号（扩张或收缩）。已有工作主要集中于非相对论情形，相对论情形仅[51]在额外小性假设下给出基本能量估计。本文旨在填补空白，建立3D线性稳定性条件并证明2D非线性解的局部存在性。

2. 主要定理或主要结果

在稳定性条件(2.53)下，3D线性化问题满足能量估计：对于基本状态 $(\hat{U}, \hat{h}, \hat{e}, \hat{\varphi})$ 满足 $\|\hat{U}\|_{W^{2,\infty}(\Omega_T)} \leq K$ ， $\|\hat{\psi}\|_{W^{2,\infty}(\Sigma_T)} \leq M$ ，以及稳定性条件

\left| \frac{1}{\epsilon}(\hat{e}_1 + \epsilon \hat{v}_2 \hat{h}_3 - \epsilon \hat{v}_3 \hat{h}_2) \right| < \sqrt{\delta^2 - \left( \frac{1}{\epsilon} \right)^2 \left| \hat{H} \times \hat{h} \right|^2},

其中 $\delta > 0$ 为常数， $m \geq 3$ ，存在常数 $C(K, M, \delta) > 0$ 使得有效线性化问题(2.48)的解满足

\sum_{|\alpha| \leq m} \|\mathcal{U}_{\alpha}\|_{L^2(\Omega_T)}^2 + \sum_{|\alpha| \leq m} \|\psi_{\alpha}\|_{L^2(\Sigma_T)}^2 \leq C(K, M, \delta) \left( \|\mathcal{U}\|_{L^2(\Omega)}^2 + \|\psi\|_{L^2(\Sigma)}^2 + \|\mathcal{F}\|_{m,*,T}^2 \right),

其中 $\mathcal{U}_{\alpha}$ 和 $\psi_{\alpha}$ 为Alinhac好未知量， $\mathcal{F}$ 为源项， $\|\cdot\|_{m,*,T}$ 为各向异性Sobolev范数。

对于2D非线性问题，定理2.3证明：若初始数据满足 $|H_0| + |h_0| \geq \delta > 0$ 在初始界面 $\Sigma(0)$ 上每一点成立，则存在 $$T > 0$$ 和唯一解 $(U, u, \varphi)$ 满足 $U \in H^m_*(\Omega_T)$ ， $u \in H^m_*(\Omega_T)$ ， $\varphi \in H^m(\Sigma_T)$ ，其中 $m \geq 15$ ， $$H^m_*$$ 为各向异性Sobolev空间（两个切向导数计为一个法向导数）。该结果通过结合线性化问题的可解性与tame估计，采用修正的Nash-Moser迭代得到。具体而言，2D线性化问题可约化为界面函数 $\psi$ 的输运方程和双曲边值问题，后者具有最大非负边界条件，其解满足tame估计

\|\mathcal{U}\|_{m,*,T} + \|\psi\|_{H^m(\Sigma_T)} \leq C \left( \|\mathcal{U}_0\|_{H^m(\Omega)} + \|\psi_0\|_{H^m(\Sigma)} + \|\mathcal{F}\|_{m,*,T} \right),

其中 $$C$$ 依赖于基本状态但不随 $$m$$ 增长而恶化。

3. 研究方法、关键技术和核心工具

主要技术障碍在于问题为特征边界且具有可变重数：边界矩阵在 $\Sigma(t)$ 上奇异且秩非常数，导致需在界面不同部分施加不同数量边界条件。作者首先通过真空散度方程(1.7)-(1.8)将原问题重新表述为具有常重数特征边界的对称化双曲问题。具体地，引入等离子体速度 $$v$$ 在真空区域的延拓 $\nu$ 满足 $\nu|_{\Sigma(t)} = v|_{\Sigma(t)}$ 且 $|\nu| < \epsilon^{-1}$ ，将麦克斯韦方程简化为

\epsilon\partial_t u + A^{-}_j(\nu)\partial_j u = 0 \quad \text{in } \Omega^{-}(t),

其中 $A^{-}_j(\nu)$ 由(2.12)-(2.13)定义，边界矩阵在 $\Sigma(t)$ 上具有恒定秩（3D情形有2个负、2个零、2个正特征值）。

对于3D线性稳定性，核心挑战是处理切向导数估计中的边界项 $Q_{\alpha}$ 。作者发现关键创新：利用边界条件(2.4)和内部方程在边界上的限制，产生内在抵消效应，将 $Q_{\alpha}$ 转换为瞬时边界积分 $R_{\alpha}$ ：

Q_{\alpha} = \int_{\Sigma_t} \left( -2\epsilon \partial_t \psi \partial_2 \psi \partial_t \psi_{\alpha} w_4 + 2\epsilon \partial_t \psi \partial_3 \psi \partial_t \psi_{\alpha} w_4 \right) d\sigma d\tau + \text{可吸收项}.

为确定定量稳定性条件，作者分离空间导数与时间导数的估计。对空间导数部分：在非共线性条件 $|\hat{H} \times \hat{h}| \geq \delta > 0$ 下，利用齐次约束(3.4)-(3.5)将 $\psi$ 的空间导数转换为法向磁场迹：

\partial_2 \psi = -\frac{W_5}{\hat{H}_1}, \quad \partial_3 \psi = -\frac{w_1}{\hat{h}_1},

再通过局部化方法将边界积分 $$T$$ 转换为涉及 $\partial_1 W_5, \partial_1 w_1, \partial_1 w_4$ 的体积积分。结合齐次约束(3.6)-(3.8)和瞬时能量显式下界，导出稳定性条件(2.53)使体积积分主要被瞬时能量吸收。

对时间导数部分：由边界条件(3.27)、恒等式(4.25)及抵消性质(4.12)， $\partial_t \psi$ 和 $\partial_t w_4$ 可表示为 $$(W_5, w_1, w_4)$$ 和 $\partial_3 \mu_2 - \partial_2 \mu_3$ 的迹，从而通过局部化论证控制 $R_{\alpha}$ 而无需额外假设。

对于2D问题，边界条件(3.27c)仅通过非特征变量 $W_{nc}$ 依赖于 $$(W, w)$$ 的迹。作者通过消除 $\psi$ 的一阶导数，将边界条件重新表述为 $\psi$ 的输运方程和两个仅含 $(W_{nc}, \psi)$ 的边界条件，使问题解耦为 $\psi$ 的输运方程和 $$(W, w)$$ 的双曲边值问题。后者在 $\psi = 0$ 时具有最大非负边界条件，通过[44]中的不动点论证建立可解性。

非线性存在性证明采用修正的Nash-Moser迭代：构造光滑中间状态（命题5.7）满足线性可解性约束，并利用2D线性化问题的tame估计（解的正则性损失固定）保证迭代收敛。关键步骤是估计中间状态对基本状态的偏离，确保每次迭代的线性化问题满足稳定性条件。

4. 与之前工作的比较

本文识别的稳定性条件(2.53)明确指定了3D线性稳定性的显式域：

\left| \frac{1}{\epsilon}(\hat{e}_1 + \epsilon \hat{v}_2 \hat{h}_3 - \epsilon \hat{v}_3 \hat{h}_2) \right| < \sqrt{\delta^2 - \left( \frac{1}{\epsilon} \right)^2 \left| \hat{H} \times \hat{h} \right|^2},

该条件不仅包含非共线性条件 $|\hat{H} \times \hat{h}| \geq \delta > 0$ ，还对电场-速度组合施加精确约束。当 $\hat{H} \times \hat{h} = 0$ 时，条件简化为 $|\hat{e}_1 + \epsilon \hat{v}_2 \hat{h}_3 - \epsilon \hat{v}_3 \hat{h}_2| < \epsilon\delta$ ，比[51]的模糊小性假设更严格且可验证。例如，若 $\hat{H} = (1,0,0)$ ， $\hat{h} = (1,0,0)$ ，则要求 $|\hat{e}_1| < \epsilon\delta$ ，而[51]仅要求 $|\hat{e}_1|$ 足够小。

本文是首个严格证明相对论等离子体-真空界面存在性的研究。2D非线性存在性定理2.3的初始条件 $|H_0| + |h_0| \geq \delta > 0$ 比[34]中非相对论2D问题的条件更弱：[34]要求 $|H| + |h| \geq \delta$ 在整个演化过程中成立，而本文仅需初始时刻满足，且允许 $$H$$ 和 $$h$$ 在部分点同时为零（只要不恒为零）。代价是仅适用于2D，3D存在性仍为开放问题，因3D边界条件涉及特征变量迹，导致2D解耦方法失效。

技术层面，本文对3D问题的处理具有原创性。先前工作如[12,14,15,26,33,34,44-46,48-50,53-56]的边界条件不涉及特征变量迹，而本文发现的内在抵消效应（将 $Q_{\alpha}$ 转换为 $R_{\alpha}$ ）和分离空间/时间导数估计的策略，为处理可变重数特征边界问题提供了新工具。对于2D问题，作者将线性化问题约化为输运方程与双曲边值问题的组合，比[34]中直接处理完整系统更精细，且通过修正的Nash-Moser迭代克服正则性损失，比Lindblad-Zhang[26]的迭代更适应各向异性空间。

适用范围上，本文结果适用于一般本构关系（满足 $0 < c_s(\rho,S) < \epsilon^{-1}$ ），但要求界面为图形式且 $|v| < \epsilon^{-1}$ 。相比非相对论情形，相对论模型中真空电场动态参与演化（而非次生变量），导致问题完全双曲化，但引入额外复杂性。本文推进了对相对论等离子体现象的数学理解，为天体物理中相对论喷流等研究提供理论基础，但3D非线性存在性仍需突破。

3. 基于Kelvin-Voigt流变律的弹性-黏性-塑性海冰模型的数学研究

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

作者: Daniel W. Boutros, Xin Liu, Marita Thomas, Edriss S. Titi

类别: math.AP, physics.ao-ph, physics.geo-ph

发布日期: 2026-04-29

arXiv ID: 2604.26295v1

摘要翻译: 受广泛应用于大尺度数值气候模拟的弹性-黏性-塑性（EVP）海冰模型［E. C. Hunke 与 J. K. Dukowicz，《物理海洋学杂志》，27 卷 9 期（1997 年），1849–1867］启发，我们在［D. W. Boutros, X. Liu, M. Thomas 与 E. S. Titi，arXiv:2505.03080（2025）］中提出了对本构关系（即应力张量关系）采用无黏性 Voigt 正则化，并证明了所得模型的整体适定性。EVP 模型将海冰视为非牛顿流体；而弹性-黏性-塑性固体通常在应变率意义上引入 Kelvin-Voigt 黏性。因此，本文构造了一个在应变率意义上带 Kelvin-Voigt 正则化的弹性-黏性-塑性海冰模型——换言之，我们将 Voigt 正则化引入动量平衡方程，而非引入本构关系（应力张量关系）。随后，我们证明了含对流项的 Kelvin-Voigt EVP 模型的局部适定性；并在不含对流项的情形下（后者沿用一种极为标准的逼近方法）证明了其整体适定性。这些结果证明中的一个关键环节，是依赖于阻尼结构所建立的一个新的应力张量 $L^\infty$ -估计。值得注意的是，无论是否含对流项，我们均能处理黏性系数无上界截断的情形——这仍是与之密切相关的 Hibler 海冰模型中一个主要未解决问题。此外，相较于我们此前关于 Voigt-EVP 模型的工作，本文还能对正则性更低的初始数据证明解的存在性。

Comment: 无

链接: http://arxiv.org/abs/2604.26295v1

1. 研究对象和背景

本研究聚焦于弹性-粘性-塑性（Elastic-Viscous-Plastic, EVP）海冰动力学模型的数学分析。该模型将海冰视为非牛顿流体，在大型数值气候模拟中具有重要应用。研究对象是定义在时间区间 $$[0, T]$$ 和二维平坦环面 $\mathbb{T}^2 = (\mathbb{R}/\mathbb{Z})^2$ 上的EVP模型：

\partial_t u = \nabla \cdot \sigma + T_a + T_w + \Omega u^\perp - g \nabla H_0, \quad (1.1a)

\frac{1}{E} \partial_t \sigma + \frac{4D}{P} \left(\sigma - \frac{1}{2} \text{Tr} \, \sigma I_2\right) + \frac{D}{2P} \text{Tr} \, \sigma I_2 + \frac{D}{2} I_2 = D(u), \quad (1.1b)

u u|_{t=0} = u_0, \quad \sigma|_{t=0} = \sigma_0, \quad (1.1c)

其中 $$P > 0$$ 为给定的内部冰强度常数， $u: \mathbb{T}^2 \times [0,T] \to \mathbb{R}^2$ 为速度场， $\sigma: \mathbb{T}^2 \times [0,T] \to \mathbb{R}^{2\times 2}_{\text{sym}}$ 为对称应力张量（关键未知量）。变形率张量定义为 $D(u) := \frac{1}{2}[\nabla u + (\nabla u)^\top]$ ，应变率 $$D = |D(u)|$$ 采用张量欧氏范数 $|A| := \sqrt{\sum_{i,j=1}^2 |A_{ij}|^2}$ 。大气和海洋拖曳力为：

T_a := c_a \rho_a |U_a| \left(U_a \cos \phi + U_a^\perp \sin \phi \right), \quad (1.4)

T_w := c_w \rho_w |U_w - u| \left[(U_w - u) \cos \theta + (U_w - u)^\perp \sin \theta \right], \quad (1.5)

科里奥利项中 $v^\perp = (-v_2, v_1)^\top$ 。EVP模型最初由Hunke和Dukowicz (1997)提出，作为粘性-塑性Hibler模型的数值正则化，后者对应形式极限 $E \to \infty$ 。先前研究[10]在应力张量本构关系中引入了Voigt正则化 $-\alpha^2 \partial_t \Delta \sigma$ ，而本文创新性地在动量平衡方程(1.1a)中引入Kelvin-Voigt正则化 $-\alpha^2 \partial_t \Delta u$ ，即研究如下修正模型：

\partial_t(u - \alpha^2 \Delta u) + u \cdot \nabla u = \nabla \cdot \sigma + T_a + T_w + \Omega u^\perp - g \nabla H_0, \quad (1.6a)

\frac{1}{E} \partial_t \sigma + \frac{4D}{P} \left(\sigma - \frac{1}{2} \text{Tr} \, \sigma I_2\right) + \frac{D}{2P} \text{Tr} \, \sigma I_2 + \frac{D}{2} I_2 = D(u), \quad (1.6b)

本文试图推进的关键空白包括：(1) 处理粘度系数无上界截断（即 $\epsilon \to 0$ 极限）的情形，这仍是Hibler模型的重大开放问题；(2) 证明更低正则性初始数据的解存在性；(3) 在含对流项的多维情况下首次严格处理 $\epsilon \to 0$ 极限，而此前[10]的估计无法处理含对流项时的此极限。

2. 主要定理或主要结果

在以下假设下证明了两个主要结果：

Theorem 1.1（含对流项的局部适定性）：设 $u_0 \in H^2(\mathbb{T}^2)$ ， $\sigma_0 \in H^1(\mathbb{T}^2) \cap L^\infty(\mathbb{T}^2)$ 且 $\sigma_0 = \sigma_0^\top$ a.e.， $U_a, U_w \in L^4((0,T); L^4(\mathbb{T}^2))$ ， $H_0 \in L^2((0,T); H^1(\mathbb{T}^2))$ 。则存在 $$T > 0$$ （依赖于初始数据及参数 $\alpha, c_a, c_w, \rho_a, \rho_w, E, g, H_0, P, \phi, \theta, U_a, U_w$ ），使得含对流项的Kelvin-Voigt EVP模型(1.6)存在唯一局部强解 $(u, \sigma)$ ，满足：

u u \in C([0, T]; H^2(\mathbb{T}^2)), \quad \sigma \in C([0, T]; H^1_w(\mathbb{T}^2)) \cap C([0, T]; L^2(\mathbb{T}^2)).

解对初始数据连续依赖（见估计式(3.25)）。

Theorem 1.2（无对流项的全局适定性）：设 $u_0 \in H^2(\mathbb{T}^2)$ ， $\sigma_0 \in H^1(\mathbb{T}^2) \cap L^\infty(\mathbb{T}^2)$ 且 $\sigma_0 = \sigma_0^\top$ a.e.， $U_a \in L^4((0,T); L^4(\mathbb{T}^2))$ ， $U_w \in L^4((0,T); L^\infty(\mathbb{T}^2))$ ， $H_0 \in L^2((0,T); H^1(\mathbb{T}^2))$ 。则对任意 $$T > 0$$ ，无对流项的Kelvin-Voigt EVP模型(1.7)存在唯一全局强解 $(u, \sigma)$ ，满足：

u u \in C([0, T]; H^2(\mathbb{T}^2)), \quad \sigma \in C([0, T]; H^1_w(\mathbb{T}^2)) \cap C([0, T]; L^2(\mathbb{T}^2)).

证明的核心是新的 $L^\infty$ 估计（Lemma 2.3）：对 $\tau = \sigma + \frac{P}{2} I_2$ ，有

\|\sigma(\cdot, t)\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P,

该估计通过在方程(2.9)上测试 $P^{-p+1}\tau|\tau|^{p-2}$ 并取 $p \to \infty$ 得到。证明中引入了应变率正则化 $D_\epsilon = \sqrt{|D(u)|^2 + \epsilon^2}$ ，先证明正则化系统(1.11)-(1.12)的解存在性，再令 $\epsilon \to 0$ 。关键技术细节包括：(1) 应用Galerkin方法时引入额外正则化 $\beta^4 \partial_t \Delta^2 u$ 以获得足够正则性；(2) 利用应力张量对称性保持性质（Lemma 2.2）；(3) 应变率差的估计 $\|\sqrt{|D(v_1)|^2 + \epsilon^2} - \sqrt{|D(v_2)|^2 + \epsilon^2}\|_{L^2} \leq \|\nabla v_1 - \nabla v_2\|_{L^2}$ （Lemma 2.1）。

3. 研究方法、关键技术和核心工具

证明中的核心障碍与创新点如下：

1. $L^\infty$ 估计的构造：最大技术突破是建立应力张量的 $L^\infty$ 估计（Lemma 2.3）。传统方法难以处理无上界截断的粘度系数，作者通过引入变换 $\tau = \sigma + \frac{P}{2} I_2$ 将方程(1.6b)转化为：

\frac{1}{E} \partial_t \tau + \frac{4D_\epsilon}{P} \left(\tau - \frac{1}{2} \text{Tr} \, \tau I_2\right) + \frac{D_\epsilon}{2P} \text{Tr} \, \tau I_2 = D(u),

然后对 $$p > 3$$ 测试函数 $P^{-p+1}\tau|\tau|^{p-2}$ 进行能量估计，得到：

\frac{P}{pE} \frac{d}{dt} \|\tau/P\|^p_{L^p} + \int_{\mathbb{T}^2} \left[ \frac{3D_\epsilon}{P^p} \left| \tau - \frac{1}{2} \text{Tr} \, \tau I_2 \right|^2 |\tau|^{p-2} + \frac{D_\epsilon}{p} \left( \frac{|\tau|}{P} \right)^p \right] dx \leq \frac{1}{p} \int_{\mathbb{T}^2} D_\epsilon dx.

关键观察是利用恒等式 $|\tau|^2 = \left| \tau - \frac{1}{2} \text{Tr} \, \tau I_2 \right|^2 + \frac{1}{2}|\text{Tr} \, \tau|^2$ 保留阻尼项，通过时间积分和 $p \to \infty$ 极限得到 $\|\tau(\cdot, t)\|_{L^\infty} \leq \|\tau_0\|_{L^\infty} + P$ ，从而导出 $\|\sigma\|_{L^\infty} \leq \|\sigma_0\|_{L^\infty} + 2P$ 。这一估计独立于 $\epsilon$ ，是处理 $\epsilon \to 0$ 极限的关键。

2. 对称性保持机制：Lemma 2.2证明应力张量的对称性在演化中保持。定义反对称部分 $A(\sigma) = \frac{1}{2}(\sigma - \sigma^\top)$ ，由方程(1.6b)导出：

\frac{1}{E} \partial_t A(\sigma) + \frac{4D}{P} A(\sigma) = 0.

在 $$L^2$$ 中测试 $A(\sigma)$ 得 $\frac{d}{dt} \|A(\sigma)\|^2_{L^2} \leq 0$ ，故若 $\sigma_0$ 对称则 $\sigma(t)$ 恒对称。这一性质对后续分析至关重要，避免了对称性破坏导致的奇异性。

3. Galerkin近似方案的改进：为克服Galerkin近似与 $L^\infty$ 估计的兼容性问题，作者在动量方程中引入双重正则化项 $\beta^4 \partial_t \Delta^2 u$ （系统(3.1)），构造解的形式为：

u u^N(x,t) = \sum_{|k|\leq N} a_k(t) e^{2\pi i k \cdot x}, \quad \sigma^N(x,t) = \sum_{|k|\leq N} b_k(t) e^{2\pi i k \cdot x}.

关键创新在于利用 $L^\infty$ 估计获得与 $\beta$ 无关的正则性界，从而在 $\beta \to 0$ 时保留解的性质。这避免了传统方法中因高阶导数估计导致的参数依赖。

4. 应变率正则化的精细处理：通过Lemma 2.1的估计 $\|\sqrt{|D(v_1)|^2 + \epsilon^2} - \sqrt{|D(v_2)|^2 + \epsilon^2}\|_{L^2} \leq \|\nabla v_1 - \nabla v_2\|_{L^2}$ ，作者控制了 $\epsilon$ 变化时解的连续性。在证明 $\epsilon \to 0$ 极限时，利用此估计和Gronwall不等式证明解序列的Cauchy性质，从而获得收敛性。这是多维情况下首次在含对流项时严格处理 $\epsilon \to 0$ 极限。

4. 与之前工作的比较

与既有文献的比较如下：

1. 与作者前期工作[10]的比较：在[10]中，作者在应力张量本构关系中引入Voigt正则化 $-\alpha^2 \partial_t \Delta \sigma$ ，证明了无对流项模型的全局适定性。本文创新点在于将正则化移至动量平衡方程 $-\alpha^2 \partial_t \Delta u$ ，这带来三重优势：(1) 物理上更合理，因Kelvin-Voigt流变学自然对应于应变率的粘性项；(2) 数学上更有效，因EVP模型的非线性涉及 $$u$$ 的导数而非 $\sigma$ ，正则化 $$u$$ 可直接控制非线性项；(3) 初始数据正则性要求显著降低——本文仅需 $u_0 \in H^2$ ，而[10]要求 $u_0 \in H^3$ 。特别地，Theorem 1.2中无需对水转向角 $\theta$ 附加限制（[10, Theorem 1.2]要求 $|\theta| < \pi/2$ ），这是因为正则化 $$u$$ 使高阶估计闭合时无需小参数假设。

2. 与Hibler模型研究的比较：Hibler模型对应 $E \to \infty$ 极限，其适定性研究存在根本困难。[48]首次证明Hibler模型局部适定性，但需修改应变率正则化（用 $D_\epsilon$ 代替最大函数）；[11]用不同方法给出证明。然而，多维Hibler模型的 $\epsilon \to 0$ 极限（即去除粘度截断）仍是重大开放问题（[48]明确指出）。本文突破在于：(1) 对含对流项的模型，首次在多维情况下严格处理 $\epsilon \to 0$ 极限（Remark 1.3）；(2) 对无对流项模型，证明全局适定性且无需 $\theta$ 限制。代价是：对含对流项模型仅得局部解（ $$T$$ 依赖初始数据），而Hibler模型即使局部解也尚未完全解决。

3. 与其他海冰模型研究的关联：[24,25]利用[10]中忽略对流项的近似，获得Hibler模型动量方程的全局存在性，但仅适用于无对流项情形。本文将此推广至Kelvin-Voigt正则化框架，且处理了更一般的初始数据。与[49]研究1D Hibler模型 $\epsilon \to 0$ 极限不同，本文解决多维情形，技术上需处理向量场耦合和更高维奇异性。值得注意的是，[52]曾对修正EVP模型进行形式 $$H^1$$ 估计，但缺乏严格证明；本文通过 $L^\infty$ 估计提供了严格基础。

4. 与非牛顿流体模型的联系：EVP模型分析与Oldroyd-B模型有相似性（均含张量演化方程），但关键区别在于EVP模型的粘性项含 $\text{Tr} \, \sigma$ 耦合。[16,18,20]等对Oldroyd-B模型的研究依赖小初值或特殊结构，而本文处理大初值且利用 $L^\infty$ 阻尼结构获得全局解（无对流项时）。Kelvin-Voigt正则化在Navier-Stokes方程中用于湍流建模（[4,19,35]），但本文首次将其应用于海冰动力学，揭示了正则化位置（动量方程vs本构关系）对解性质的决定性影响。

5. 局限性与推进方向：本文主要局限在于含对流项时仅得局部解，而物理上长期行为更重要。与Hibler模型相比，本文模型因 $E < \infty$ 保留弹性项，避免了Hibler模型的数值刚性，但牺牲了 $E \to \infty$ 极限的严格联系。未来工作需解决：(1) 含对流项模型的全局适定性；(2) $E \to \infty$ 时向Hibler模型的收敛性；(3) 三维情形的推广。值得注意的是，本文方法可能适用于其他含张量演化的地质流体模型（如冰盖动力学[55]）。

了解领域论文（摘要翻译）

1. 带源项与漂移项的Hele-Shaw流自由边界之 $$C^1$$ -正则性

$$C^1$$ -Regularity of the Free Boundary for Hele-Shaw Flow with Source and Drift

作者: Yuming Paul Zhang

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26906v1

链接: http://arxiv.org/abs/2604.26906v1

摘要翻译

中文标题: 带源项与漂移项的Hele-Shaw流自由边界之 $$C^1$$ -正则性

摘要翻译: 本文是文献 \cite{kimzhang2024} 关于带漂移项与源项的Hele-Shaw流工作的延续。我们证明：在局部邻域内，若自由边界为Lipschitz连续且其Lipschitz常数充分小，则该自由边界属于 $C^{1}$ 类。作为推论，我们还考虑二维竖直方向Hele-Shaw问题（即单相Muskat问题）带对流项的情形。我们证明：只要初始数据与对流项足够小，且传播速度足够大，则自由边界在有限时间后一致地属于 $$C^1$$ 类。

2. 逆流换热器的谱边界观测器

Spectral Boundary Observer for Counter-Flow Heat Exchangers

作者: Mohamed Camil Belhadjoudja, Mohamed Maghenem, Emmanuel Witrant

类别: eess.SY, math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26863v1

链接: http://arxiv.org/abs/2604.26863v1

摘要翻译

中文标题: 逆流换热器的谱边界观测器

摘要翻译: 我们考虑一个由两个耦合的一阶线性双曲型偏微分方程组成的系统，该系统建模逆流换热器中的热量传输：其中一个方程描述高温流体的输运，另一个方程描述低温流体在相反方向上的输运。针对该系统，我们设计了一种仅依赖于低温流体在一端边界处所测温度的边界观测器。我们的方法是谱方法：通过将支配观测误差动力学的算子的谱指定于开左半复平面内的某个预设区域，可在 $$L^2$$ 范数意义下自由调节观测误差收敛至零的速率。主要技术贡献在于证明了谱稳定性（即该算子的谱位于开左半复平面内）等价于观测误差动力学在原点处的 $$L^2$$ 指数稳定性；该等价性是通过证明该支配观测误差动力学的算子满足所谓谱映射性质（spectral mapping property）而建立的。

3. 具有不匹配密度与奇异势的拟不可压两相流扩散界面模型的整体弱解

Global weak solutions to a diffuse-interface model for quasi-incompressible two-phase flows with unmatched densities and singular potential

作者: Mingwen Fei, Xiang Fei, Yadong Liu, Hao Wu

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26660v1

链接: http://arxiv.org/abs/2604.26660v1

摘要翻译

中文标题: 具有不匹配密度与奇异势的拟不可压两相流扩散界面模型的整体弱解

摘要翻译: 我们研究一个热力学协调的扩散界面模型，该模型描述两种宏观上互不混溶、不可压且具粘性的牛顿流体在密度不匹配情形下的运动。该模型与连续介质混合物理论相容。它采用质量平均（质心）速度，从而使得两相流为拟不可压：速度场不再散度为零，且压力项出现在化学势方程中。对于定义在环面 $\mathbb{T}^3$ 上、初边值问题及一类物理上合理的奇异自由能密度，我们证明了整体时间弱解的存在性。证明依赖于将原系统适当地约化为一类 Korteweg 型流体模型，并结合双层逼近；同时借助受著名 Bresch–Desjardins 熵方法启发的关于质量密度与相场变量的精细估计。关键观察在于自由界面处的毛细效应为密度演化提供了阻尼作用。在极限过程分析中，由于压力的可积性在先验上不可得，我们推导出精细的尾部估计以排除奇异势可能出现的集中现象。本工作似乎是首例关于具有不匹配密度与质量平均速度、且未施加空间正则化项的 Navier–Stokes/Cahn–Hilliard 型系统的整体存在性结果。

不相关论文（基本信息）

1. 生长组织中黏附驱动自组织的理论

Theory of adhesion-driven self-organisation in growing tissues

作者: Carles Falcó, Samuel W. S. Johnson, Mohit P. Dalwadi, Philip K. Maini

类别: q-bio.TO, math.AP, physics.bio-ph

发布日期: 2026-04-29

arXiv ID: 2604.26928v1

链接: http://arxiv.org/abs/2604.26928v1

摘要: Cell invasion and spatial pattern formation are two distinct manifestations of cellular self-organisation in development, regeneration, and disease. Here, we develop and analyse a unified theoretical framework that links these two seemingly different behaviours within a single mechanistic model for adhesion-mediated self-organisation in growing cell populations. Using a multiscale analysis, we show that the balance between cell-cell adhesion, self-diffusion, and proliferation controls the emergence of distinct collective dynamics. We find that for weak adhesion, tissues invade through stable monotone fronts. As adhesion increases, invasion slows, fronts become unstable, leading to aggregates and spatial patterns emerging behind the advancing edge. In two spatial dimensions, these instabilities generate fingering morphologies reminiscent of dysregulated invasion in cancer. Crucially, we show that density-dependent regulation of adhesion suppresses these instabilities and restores cohesive tissue expansion. Together, our results identify adhesion strength and its regulation as key determinants of whether tissues invade cohesively or fragment into patterns, and provide a unified framework for understanding collective migration, morphogenesis, and dysregulated growth.

2. 弱奇异敏感性下斑秃趋化系统解的整体有界性与渐近行为

Global boundedness and asymptotic behavior of the chemotaxis system for Alopecia Areata with weakly singular sensitivity

作者: Pengxue Xiang, Yuebo Cao, Hongying Yang

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26905v1

链接: http://arxiv.org/abs/2604.26905v1

摘要: This paper considers the homogeneous Neumann initial-boundary value problem for Alopecia Areata chemotaxis model with weakly singular sensitivity. For any appropriately regular initial conditions,it is shown that the problem admits a global boundedness of classical solutions in two spatial dimensions. Moreover, through the explicit construction of Lyapunov functions, we establish that the globally bounded solution converges exponentially to a constant steady state. The paper concludes with numerical experiments that serve to visually illustrate and corroborate some of the theoretically derived findings.

3. 关于仿射Quermass积分的单调性

On the monotonicity of affine quermassintegrals

作者: Shibing Chen, Yuanyuan Li, Xianduo Wang

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26828v1

链接: http://arxiv.org/abs/2604.26828v1

摘要: Lutwak’s affine quermassintegral theory is a foundational component of modern affine Brunn–Minkowski theory. Developed in the 1980s, it provides affine analogues of the classical quermassintegrals and has led to a rich family of sharp affine isoperimetric inequalities. A central question in this program, going back to Lutwak’s 1988 work, is an Alexandrov–Fenchel-type monotonicity principle for the normalized $L^{-n}$ -moment quermassintegrals $I_{k,-n}$ . In one form, this principle predicts that

I_{m,-n}(K)^{1/m}\ge I_{k,-n}(K)^{1/k}, \qquad 1\le m<k\le n .

The question was recorded in Gardner’s 2006 book Geometric Tomography in its problem list, and it was later placed by Milman and Yehudayoff in their 2023 JAMS paper within a broader $$L^p$$ -moment framework. Their sharp affine quermassintegral theorem verifies the conjectured comparison in a substantial range; in particular, at the affine-invariant endpoint $$p=-n$$ , it includes the comparison with the top dimension. We show that the proposed monotonicity does not persist in the full range. More precisely, for every triple of integers $$m,k,n$$ satisfying $1\leq m<k\leq n-1$ and $$n>(m+2)(k+2)-2$$ , there exists an origin-symmetric $$C^2_+$$ convex body $K\subset\mathbb R^n$ such that

I_{m,-n}(K)^{1/m} < I_{k,-n}(K)^{1/k}.

The example is obtained from the Euclidean ball by an arbitrarily small degree-four spherical harmonic perturbation. The proof is based on a second-variation formula for negative projection-volume moments, together with an explicit Grassmannian average for the fourth zonal harmonic in arbitrary projection dimension.

4. 抛物型 $$p$$ -拉普拉斯方程解的概率表示

Probabilistic representation of solutions to the parabolic $$p$$ -Laplace equation

作者: Viorel Barbu, Michael Röckner

类别: math.AP, math.PR

发布日期: 2026-04-29

arXiv ID: 2604.26719v1

链接: http://arxiv.org/abs/2604.26719v1

摘要: This work is concerned with the probabilistic representation of solutions to the $$p$$ -Laplace evolution equation $\frac{\partial u}{\partial t}={\rm div}(|\nabla u|^{p-2}\nabla u)$ in $(0,\infty)\times\mathbb{R}^d$ , $$u(0,x)=u_0(x),$$ $x\in\mathbb{R}^d$ . One proves that, if $p\geq 4$ , and if $$u_0$$ is a probability density with compact support and $u_0\in L^2$ , $|\nabla u_0|\in L^\infty$ , then $$u$$ can be represented as $u(t,x)dx=\mathscr L_{X(t)}(dx)$ , where $\mathscr L_{X(t)}$ denotes the time marginal law of $$X$$ at time $$t$$ with $$X$$ being a probabilistically weak solution to a corresponding McKean-Vlasov stochastic differential equation. This result is based on a new second order global regularity result for the weak solutions to the parabolic $$p$$ -Laplace equation.

5. 对数障碍问题中奇性自由边界点的几何结构

Geometric structure of singular free boundary points for the logarithmic obstacle problem

作者: Lili Du, Xu Tang, Yi Zhou

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26485v1

链接: http://arxiv.org/abs/2604.26485v1

摘要: In the previous work [Interfaces Free Bound., 19, 351–369, 2017], de Queiroz and Shahgholian established the optimal $C^{1,\log}_{\mathrm{loc}}$ regularity of solutions for the obstacle problem with singular logarithmic forcing term

-Δu = \log u\,χ_{\{u>0\}} \quad \text{in } Ω,

where $Ω\subset\mathbb{R}^d$ ( $d\geq 2$ ) is a smooth bounded domain. In our earlier work [arXiv:2408.08104, 2024], we proved the $C^{1,α}$ regularity of the free boundary $Ω\cap\partial\{u>0\}$ near regular points. In this paper, we investigate the more delicate structure of the \emph{singular} free boundary. Since the nonlinearity $-\log u$ is singular near the free boundary and destroys the scaling invariance, so that neither the classical blow-up arguments nor the standard epiperimetric inequality [Weiss, Invent.\ Math., 138, 23–50, 1999] apply directly; moreover, the Weiss type monotonicity formula requires a variable-parameter correction that introduces non-integrable remainder terms into the energy estimates. Motivated by Colombo–Spolaor–Velichkov [Geom.\ Funct.\ Anal., 28, 1029–1061, 2018], we develop a new \emph{log-epiperimetric inequality} for the modified Weiss energy, also proved by the direct method. A key novelty is the introduction of an auxiliary correction term $$T$$ that absorbs the non-integrable errors. As consequences, we establish a logarithmic energy decay, uniqueness of blow-ups at singular points, and a $C^{1,\log}$ -type geometric description of the singular strata. In dimension two, the logarithmic modulus improves to a Hölder modulus.

6. 局部Calderón问题与复各向异性电导率在边界上的确定

The local Calderón problem and the determination at the boundary of a complex anisotropic admittivity

作者: Jessica Crosse, Romina Gaburro

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26458v1

链接: http://arxiv.org/abs/2604.26458v1

摘要: We address Calderón’s problem of stably determining the anisotropic complex admittivity $$σ$$ in a domain $Ω\subset\mathbb{R}^n$ , with $n\geq3$ , representing a conducting medium, in terms of a Dirichlet-to-Neumann map locally prescribed on a non-empty portion $$Σ$$ of the boundary of $$Ω$$ , $\partialΩ$ . $$σ$$ is assumed to be of type $σ(\cdot)=A(\cdot,a(\cdot))$ in $$Ω$$ , where the one-parameter family of complex-symmetric matrices $[λ^{-1},\:λ]\ni t\mapsto A(\cdot,\: t)$ is assumed to be a-priori known and the scalar function $$a$$ is unknown. We establish Lipschitz and Hölder stability estimates at the boundary for $$σ$$ and its derivatives of arbitrary order on $$Σ$$ , respectively, in terms of the local map.

7. 恰加斯病野外传播模型扩散过程生成的解析半群

Analytic semigroup generated by the dispersal process of a sylvatic transmission model of Chagas disease

作者: Narimene Benarbia, Rabah Labbas, Tewfik Mahdjoub, Alexandre Thorel

类别: math.AP, math.FA, math.SP

发布日期: 2026-04-29

arXiv ID: 2604.26436v1

链接: http://arxiv.org/abs/2604.26436v1

摘要: In this work, we develop a new biological transmission model for Chagas disease. This model, set in two juxtaposed habitats with skew Brownian motion conditions at the interface, is composed of two reaction–diffusion equations and takes into account the sylvatic transmission. We write it as an abstract perturbed Cauchy problem using operator theory. Then, we show that the main operator, which models the dispersal process, generates an analytic semigroup in an adequate Banach space.

8. 广义Hénon方程的Neumann问题：局部分析

The Neumann problem for the generalized Hénon equation. Local analysis

作者: Alexander I. Nazarov, Alexandra P. Shcheglova

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26286v1

链接: http://arxiv.org/abs/2604.26286v1

摘要: For the boundary value problem

\left\{ \begin{array}{rcll} -Δ_p u+u^{p-1}&=&|x|^αu^{q-1}&\mbox{in }Ω,\\ \frac{\displaystyle\partial u}{\displaystyle\partial{\bf n}}&=&0&\mbox{on }\partial Ω, \end{array}\right.

in the unit ball $$Ω$$ , we investigate the properties of the positive radial solution. It is known, that for $$1<p<n$$ , $\frac{(n-1)p}{n-p}<q<\frac{np}{n-p}$ and sufficiently large $$α$$ this solution does not provide global minimum to the corresponding energy functional, see [M. Gazzini, E. Serra, 2008] for $$p=2$$ and [A.P. Shcheglova, 2018] in general case. Nevertheless, it is shown in [M. Gazzini, E. Serra, 2008] that for $n\ge 4$ , $$p=2$$ , $2<q<\frac{2n}{n-2}$ and sufficiently large $$α$$ the radial solution is at least a local minimizer of the energy functional. We partially generalize this result. Namely, let $n\ge4$ and let $$p>2$$ be sufficiently close to $$2$$ . Then for all $p<q<\frac{np}{n-p}$ , for sufficiently large $$α$$ the second variation of the energy functional is positive. The same holds true for all $$2<p<n$$ if $$q>p$$ is sufficiently close to $$p$$ .

9. 二维超临界相位下拉格朗日平均曲率方程解在无穷远处的最优渐近行为

Optimal Asymptotic Behavior at Infinity of Solutions to the Lagrangian Mean Curvature Equation with Supercritical Phase in Dimension Two

作者: Jiguang Bao, Qinfeng Jiang

类别: math.AP

发布日期: 2026-04-29

arXiv ID: 2604.26246v1

链接: http://arxiv.org/abs/2604.26246v1

摘要: We employ a nonlocal method to study the asymptotic behavior at infinity ofsolutions to the two-dimensional supercritical Lagrangian mean curvature equation

\arctan λ_1(D^2u)+\arctan λ_2(D^2u) = θ+ f(x)

on exterior domains in $\mathbb{R}^2$ , where $|θ| \in (0, π)$ is a constant and $$f$$ is a Lipschitz continuous perturbation satisfying $f(x) = O(|x|^{-β})$ with decay rate $$β> 0$$ at infinity. This work generalizes the convergence results in \cite{BJ2026}, where $$f$$ is required to be at least $$C^3$$ and $$β>2$$ . Moreover, all asymptotic results established in this paper are optimal.

10. 关于四维常数量曲率收缩梯度Ricci孤立子的一个注记

A note on four dimensional Shrinking Gradient Ricci Solitons with Constant Scalar Curvature

作者: Chen Wang, Guoqiang Wu

类别: math.DG, math.AP

发布日期: 2026-04-28

arXiv ID: 2604.26163v1

链接: http://arxiv.org/abs/2604.26163v1

摘要: Let $$(M^4, g, f)$$ be a four-dimensional complete noncompact gradient shrinking Ricci soliton with the equation $Ric+\nabla^2f= \frac{1}{2}g$ . If its scalar curvature is $$1$$ , Cheng-Zhou \cite{Cheng-Zhou} proved that it is a finite quotient of $\mathbb{R}^2\times \mathbb{S}^2$ . In this note we present an alternative proof by analyzing the asymptotic geometry at infinity.

11. 分数阶亥姆霍兹方程中非均匀性支撑集的重构

Reconstruction of the Support of an Inhomogeneity for the Fractional Helmholtz Equation

作者: Dana Zilberberg

类别: math.AP

发布日期: 2026-04-28

arXiv ID: 2604.26154v1

链接: http://arxiv.org/abs/2604.26154v1

摘要: We consider the inverse scattering problem for inhomogeneous media of compact support governed by the fractional s-Helmholtz equation, with $$0<s<1$$ , in dimensions $$d=1,2,3$$ . In particular, we study the determination of the support of the inhomogeneity from the far-field pattern of the scattered field generated by plane waves for all incident directions at a fixed frequency. The far-field pattern is defined as the principal term in the asymptotic expansion of the scattered field at infinity. It is shown in \cite{zilberberg2026limiting} that, up to a multiplicative constant, this coincides with the far-field pattern corresponding to the classical Helmholtz equation with the same inhomogeneity. Our approach is based on the development of the factorization method, which not only leads to an efficient and easily implementable reconstruction algorithm, but also provides a uniqueness result for determining the support of an admissible set of inhomogeneities. A fundamental ingredient in the analysis is a new transmission eigenvalue problem, whose eigenvalues must be excluded. Therefore, we prove that they are discrete with no finite accumulation points. We present numerical examples in dimension $$d=2$$ for both the direct and inverse problems.

12. Schrödinger问题的Benamou–Brenier公式的偏微分方程方法

A PDE approach to Benamou–Brenier formula for the Schrödinger problem

作者: Mattia Garatti, Luca Nenna, Simona Rota Nodari, Luca Tamanini

类别: math.OC, math.AP

发布日期: 2026-04-28

arXiv ID: 2604.26023v1

链接: http://arxiv.org/abs/2604.26023v1

摘要: We studied the Benamou–Brenier formulation of the Schrödinger problem, focusing on a gap between theoretical results and applications, that often involve measures with unbounded support. While the existing proof in the literature relies on the compactness of the marginals’ supports to ensure the necessary regularity of the Schrödinger potentials, we extend the validity of the Benamou–Brenier formula to the larger class of sub-Gaussian probability measures. Exploiting fine estimates on the Hessian of the potentials and the entropic interpolation, we provide an almost self-contained proof that establishes the existence of a velocity field with the appropriate polynomial growth that ensures the right integrability. This result justifies the use of the dynamic formulation in more general settings, such as Gaussian and mixture-of-Gaussians models, important also for the applications.

重点关注论文（完整分析）#

1. 具记忆效应的非局部守恒律系统及其零记忆极限#

1. 研究对象和背景#

2. 主要定理或主要结果#

3. 研究方法、关键技术和核心工具#

4. 与之前工作的比较#

2. 相对论等离子体–真空界面的稳定性与存在性#

1. 研究对象和背景#

2. 主要定理或主要结果#

3. 研究方法、关键技术和核心工具#

4. 与之前工作的比较#

3. 基于Kelvin-Voigt流变律的弹性-黏性-塑性海冰模型的数学研究#

1. 研究对象和背景#

2. 主要定理或主要结果#

3. 研究方法、关键技术和核心工具#

4. 与之前工作的比较#

了解领域论文（摘要翻译）#

1. 带源项与漂移项的Hele-Shaw流自由边界之 $C^1$-正则性#

摘要翻译#

2. 逆流换热器的谱边界观测器#

摘要翻译#

3. 具有不匹配密度与奇异势的拟不可压两相流扩散界面模型的整体弱解#

摘要翻译#

不相关论文（基本信息）#

1. 生长组织中黏附驱动自组织的理论#

2. 弱奇异敏感性下斑秃趋化系统解的整体有界性与渐近行为#

3. 关于仿射Quermass积分的单调性#

4. 抛物型 $p$-拉普拉斯方程解的概率表示#

5. 对数障碍问题中奇性自由边界点的几何结构#

6. 局部Calderón问题与复各向异性电导率在边界上的确定#

7. 恰加斯病野外传播模型扩散过程生成的解析半群#

8. 广义Hénon方程的Neumann问题：局部分析#

9. 二维超临界相位下拉格朗日平均曲率方程解在无穷远处的最优渐近行为#

10. 关于四维常数量曲率收缩梯度Ricci孤立子的一个注记#

11. 分数阶亥姆霍兹方程中非均匀性支撑集的重构#

12. Schrödinger问题的Benamou–Brenier公式的偏微分方程方法#

重点关注论文（完整分析）

1. 具记忆效应的非局部守恒律系统及其零记忆极限

1. 研究对象和背景

2. 主要定理或主要结果

3. 研究方法、关键技术和核心工具

4. 与之前工作的比较

2. 相对论等离子体–真空界面的稳定性与存在性

1. 研究对象和背景

2. 主要定理或主要结果

3. 研究方法、关键技术和核心工具

4. 与之前工作的比较

3. 基于Kelvin-Voigt流变律的弹性-黏性-塑性海冰模型的数学研究

1. 研究对象和背景

2. 主要定理或主要结果

3. 研究方法、关键技术和核心工具

4. 与之前工作的比较

了解领域论文（摘要翻译）

1. 带源项与漂移项的Hele-Shaw流自由边界之 $$C^1$$ -正则性

摘要翻译

2. 逆流换热器的谱边界观测器

摘要翻译

3. 具有不匹配密度与奇异势的拟不可压两相流扩散界面模型的整体弱解

摘要翻译

不相关论文（基本信息）

1. 生长组织中黏附驱动自组织的理论

2. 弱奇异敏感性下斑秃趋化系统解的整体有界性与渐近行为

3. 关于仿射Quermass积分的单调性

4. 抛物型 $$p$$ -拉普拉斯方程解的概率表示

5. 对数障碍问题中奇性自由边界点的几何结构

6. 局部Calderón问题与复各向异性电导率在边界上的确定

7. 恰加斯病野外传播模型扩散过程生成的解析半群

8. 广义Hénon方程的Neumann问题：局部分析

9. 二维超临界相位下拉格朗日平均曲率方程解在无穷远处的最优渐近行为

10. 关于四维常数量曲率收缩梯度Ricci孤立子的一个注记

11. 分数阶亥姆霍兹方程中非均匀性支撑集的重构

12. Schrödinger问题的Benamou–Brenier公式的偏微分方程方法