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

重点关注论文数量: 5

了解领域论文数量: 10

不相关论文数量: 18

重点关注论文(完整分析)

1. 建模血管网络形成的具密度依赖黏性的趋化-可压Navier-Stokes方程

Chemotaxis compressible Navier-Stokes equations with density-dependent viscosity modeling vascular network formation

作者: Ansgar Jüngel, Flora Philipp

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24725v1

摘要翻译: 本文在三维环面 $\mathbb{T}^3$ 上,对具有能量有限初值的内皮细胞密度及其速度所满足的可压Navier-Stokes方程,耦合一个趋化因子浓度所满足的反应-扩散方程,证明了全局弱解的存在性。该耦合通过趋化力(出现在动量平衡方程中)以及细胞对趋化信号的产生(出现在趋化方程中)实现。该模型刻画了血管生成早期阶段内皮细胞的自组装过程。存在性结果适用于绝热压力指数 $γ>4/3$,该指数与退化Keller-Segel方程存在性分析中所得指数一致。证明借助于含Korteweg项与拖曳项的逼近、BD熵不等式,以及在速度变量上被重正规化的弱解的构造。

Comment: 无

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

1. 研究对象和背景

本文研究内皮细胞自组装过程中血管网络形成的数学模型,该模型由描述内皮细胞密度$\rho$和速度场$v$的可压缩Navier-Stokes方程,与描述趋化剂浓度$c$的反应-扩散方程耦合而成,研究区域为三维环面$\Omega$。方程系统为:

$$ \partial_t \rho + \text{div}(\rho v) = 0, $$

$$ \partial_t(\rho v) + \text{div}(\rho v \otimes v) + \nabla p(\rho) = \nu \text{div}(\rho D(v)) + \rho \nabla c - \frac{\rho v}{\zeta}, $$

$$ \partial_t c = \Delta c - c + \rho, $$

其中$p(\rho) = \rho^\gamma$为绝热压力,$D(v) = \frac{1}{2}(\nabla v + \nabla v^T)$,$\nu > 0$为粘性系数,$\zeta > 0$为松弛时间,初始条件为$\sqrt{\rho}(0, \cdot) = \sqrt{\rho_0}$,$(\sqrt{\rho}v)(0, \cdot) = \sqrt{\rho_0}v_0$,$c(0, \cdot) = c_0$。该模型源于血管发生过程,当忽略惯性和粘性时退化为退化Keller-Segel方程$\partial_t \rho + \zeta \text{div}(\nabla p(\rho) - \rho \nabla c) = 0$,其全局经典解存在当且仅当$\gamma > 4/3$[21],且当$\gamma \geq 4/3$时存在爆破[24]。在常数粘性假设下,文献[1]证明了$\gamma \geq 4$时全局弱解存在,[12]扩展至$\gamma > 8/5$,其关键依赖自由能$E_0(\rho, v, c) = \int_\Omega \left(\frac{1}{2}\rho|v|^2 + h(\rho) + \frac{1}{2}(|\nabla c|^2 + c^2) - \rho c\right) dx$与Sugiyama不等式。当粘性依赖于密度时,Bresch-Desjardins的BD熵$E_{BD}(\rho, v) = \frac{1}{2} \int_\Omega \rho|v + \nu \nabla \log \rho|^2 dx$提供$\nabla \sqrt{\rho}$的$L^2$一致界,但趋化力项$\rho \nabla c$引入交叉项$K_1 = \nu \int_0^t \int_\Omega \nabla \rho \cdot \nabla c \, dx d\tau$,其估计需控制$\int_0^t \int_\Omega \rho^2 \, dx d\tau$。本文旨在填补的空白是:在密度依赖粘性框架下,将$\gamma$的可接受范围推进至$\gamma > 4/3$(与退化Keller-Segel方程临界指数一致),并探究引入Korteweg项是否可覆盖$\gamma > 1$的完整范围。

2. 主要定理或主要结果

在$\gamma > 4/3$的假设下,本文证明了趋化性可压缩Navier-Stokes系统(1)-(4)全局弱解的存在性,且该结果与退化Keller-Segel方程的临界指数相匹配。具体而言:

定理1(趋化性Navier-Stokes-Korteweg系统的存在性):设$\kappa, r_0, r_1 > 0$,$T > 0$,$\gamma > 1$,且初始条件满足

$$ \rho_0 \in L^\gamma(\Omega), \quad \sqrt{\rho_0}, c_0 \in H^1(\Omega), \quad \sqrt{\rho_0}v_0 \in L^2(\Omega; \mathbb{R}^3), \quad r_0 \log \rho_0 \in L^1(\Omega). $$

则存在系统(4),(12)-(13)在$[0,T]$上的弱解。若$\gamma > 4/3$,则正则性结果(16)与参数$(\kappa, r_0, r_1)$无关,即:

$$ \rho \in L^\infty(0, T; L^\gamma(\Omega)), \quad \rho^{\gamma/2} \in L^2(0, T; H^1(\Omega)), $$

$$ \nabla\sqrt{\rho}, \sqrt{\rho}v \in L^\infty(0, T; L^2(\Omega; \mathbb{R}^3)), \quad T_\nu \in L^2(\Omega_T; \mathbb{R}^{3\times 3}), $$

$$ c \in L^\infty(0, T; H^1(\Omega)) \cap H^1(0, T; L^2(\Omega)), \quad \sqrt{r_0}v \in L^2(\Omega_T), $$

$$ \sqrt{\kappa}\Delta\sqrt{\rho} \in L^2(\Omega_T), \quad \nabla^4\sqrt{\kappa\rho}, \quad ^4\sqrt{r_1\rho}v \in L^4(\Omega_T; \mathbb{R}^3). $$

定理2(趋化性Navier-Stokes系统的存在性):设$T > 0$,$\gamma > 4/3$,且初始条件满足(21),则存在系统(1)-(4)在$[0,T]$上的弱解。

关键估计来自引理3:对任意$\sigma_1, \sigma_2 > 0$,存在$C > 0$使得:

$$ \int_\Omega \rho c \, dx \leq \begin{cases} C + \sigma_1 \int_\Omega \rho^\gamma \, dx + \sigma_2 \int_\Omega |\nabla c|^2 \, dx & \text{if } \gamma > 8/5, \\ C + \sigma_1 \int_\Omega |\nabla\sqrt{\rho}|^2 \, dx + \sigma_2 \int_\Omega |\nabla c|^2 \, dx. \end{cases} $$

对任意$\sigma > 0$,存在$C > 0$使得:

$$ \int_0^t \int_\Omega \rho^2 \, dx d\tau \leq \begin{cases} Ct + \int_0^t \int_\Omega \rho^\gamma \, dx d\tau & \text{if } \gamma \geq 2, \\ Ct + \sigma \int_0^t \int_\Omega |\nabla \rho^{\gamma/2}|^2 \, dx d\tau & \text{if } 4/3 < \gamma < 2, \\ Ct + \sigma \int_0^t \int_\Omega |\nabla^4 \sqrt{\rho}|^4 \, dx d\tau. \end{cases} $$

当$4/3 < \gamma < 2$时,关键不等式(9)为:

$$ \int_0^t \int_\Omega \rho^2 \, dx d\tau \leq C(\delta) + \delta \int_0^t \int_\Omega |\nabla \rho^{\gamma/2}|^2 \, dx d\tau, $$

其中$\delta > 0$任意小,右边积分可由耗散项$D_1 = \frac{4}{\gamma} \nu \int_0^t \int_\Omega |\nabla \rho^{\gamma/2}|^2 \, dx d\tau$吸收。当$\gamma > 1$时,需更强不等式(10):

$$ \int_0^t \int_\Omega \rho^2 \, dx d\tau \leq C(\delta) + \delta \int_0^t \int_\Omega |\nabla^4 \sqrt{\rho}|^4 \, dx d\tau, $$

该估计通过Korteweg项提供的额外耗散$D_2 = \kappa \int_0^t \int_\Omega \rho |\nabla^2 \log \rho|^2 \, dx d\tau$实现,且满足下界估计(11):

$$ D_2 \geq 2C_K \kappa \int_\Omega \left(|\Delta \sqrt{\rho}|^2 + |\nabla^4 \sqrt{\rho}|^4 \right) dx, \quad C_K = \frac{1}{16}. $$

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

本文核心障碍在于趋化力项$\rho \nabla c$与密度依赖粘性项$\nu \text{div}(\rho D(v))$的耦合,尤其在$\gamma \to 1^+$时,标准能量估计无法控制交叉项$K_1 = \nu \int_0^t \int_\Omega \nabla \rho \cdot \nabla c \, dx d\tau$及其导致的$\int_0^t \int_\Omega \rho^2 \, dx d\tau$。作者创新性结合BD熵不等式与自由能估计,并引入Korteweg项作为近似方案。具体技术包括:

BD熵不等式的应用:利用$E_{BD}(\rho, v) = \frac{1}{2} \int_\Omega \rho |v + \nu \nabla \log \rho|^2 dx$,结合Gagliardo-Nirenberg不等式,得到关键估计(8):

$$ \int_\Omega \rho c \, dx \leq C(\delta) + \delta \int_\Omega (|\nabla \sqrt{\rho}|^2 + |\nabla c|^2) dx, $$

其中$\delta > 0$任意小,该式由$E_0$和$E_{BD}$控制。但$K_1$的估计需控制$\rho^2$,故需对$\int_0^t \int_\Omega \rho^2 \, dx d\tau$进行尺度分析。

关键不等式与尺度分析:对$4/3 < \gamma < 2$,利用Gagliardo-Nirenberg不等式:

$$ \int_\Omega \rho^2 dx = \|\rho^{\gamma/2}\|^{4/\gamma}_{L^{4/\gamma}(\Omega)} \leq C + C\|\nabla \rho^{\gamma/2}\|^{4\theta/\gamma}_{L^2(\Omega)}, \quad \theta = \frac{3\gamma}{6\gamma-2} < 1, $$

且$4\theta/\gamma < 2$当且仅当$\gamma > 4/3$。通过Young不等式,得到:

$$ \int_\Omega \rho^2 dx \leq C(\sigma) + \sigma \|\nabla \rho^{\gamma/2}\|^2_{L^2(\Omega)}, $$

使该估计可被耗散项$D_1 = \frac{4}{\gamma} \nu \int_0^t \int_\Omega |\nabla \rho^{\gamma/2}|^2 \, dx d\tau$吸收。

Korteweg项的引入:为覆盖$\gamma > 1$,引入Korteweg项$\text{div} K = \kappa \rho \nabla (\Delta \sqrt{\rho}/\sqrt{\rho})$,其耗散为:

$$ D_2 = \kappa \int_0^t \int_\Omega \rho |\nabla^2 \log \rho|^2 \, dx d\tau. $$

利用引理[14, Appendix]与[22, Lemma 2.1],得到下界估计(11):

$$ D_2 \geq 2C_K \kappa \int_\Omega \left(|\Delta \sqrt{\rho}|^2 + |\nabla^4 \sqrt{\rho}|^4 \right) dx, \quad C_K = \frac{1}{16}. $$

结合Gagliardo-Nirenberg不等式($\theta = 3/8$),得到:

$$ \int_\Omega \rho^2 dx = \|^4\sqrt{\rho}\|^8_{L^8(\Omega)} \leq C(\sigma) + \sigma \|\nabla^4 \sqrt{\rho}\|^4_{L^4(\Omega)}, $$

其中$8\theta = 3 < 4$确保估计成立。

两层正则化方案:第一层添加Korteweg项与拖曳项$(r_0, r_1)$,构造系统(12)-(13);第二层添加高阶正则化项$(\delta, \varepsilon, \eta, \mu)$如$\varepsilon \Delta \rho$与$\mu \Delta^2 v$。证明框架为:

  1. 通过Faedo-Galerkin近似证明正则化系统(24)-(25)解的存在性,利用修改能量
    $$ E_K(\rho, v, c) = \int_\Omega \left(\frac{1}{2}\rho|v|^2 + h(\rho) + \kappa|\nabla \sqrt{\rho}|^2 + \frac{1}{2}(|\nabla c|^2 + c^2) \right) dx $$
    获得与Galerkin维数$n$无关的一致估计;
  2. 证明BD熵不等式,获得额外一致估计,允许$(\delta, \varepsilon, \eta, \mu) \to 0$,得到带Korteweg项系统的弱解;
  3. 采用Lacroix-Violet和Vasseur[15]的方法,证明该解为速度变量中重归一化的弱解,从而允许$(\kappa, r_0, r_1) \to 0$,得到原始系统(1)-(4)的弱解。

重归一化解的构造:定义张量$\sqrt{\kappa \rho} S_\kappa = \kappa \sqrt{\rho} (\nabla^2 \sqrt{\rho} - 4 \nabla^4 \sqrt{\rho} \otimes \nabla^4 \sqrt{\rho})$与$\sqrt{\nu \rho} T_\nu = \nu \nabla(\rho v) - 2 \nu \sqrt{\rho} v \otimes \nabla \sqrt{\rho}$,其在密度为零处仍有定义。通过定义3中的重归一化解框架,建立弱解与重归一化解的等价性,这是处理$\kappa \to 0$极限的关键。

4. 与之前工作的比较

本文与既有PDE/分析理论文献的比较如下:

与标准可压缩Navier-Stokes方程的比较:Lions[17]首次证明常数粘性下$\gamma > 9/5$时全局弱解存在,Feireisl等[11]扩展至$\gamma > 3/2$,依赖粘性系数有正下界。当粘性依赖密度时,Bresch等[7]通过Korteweg项克服密度为零处的困难,Bresch和Desjardins[5]移除Korteweg项,依赖BD熵提供$\nabla \sqrt{\rho}$控制。Mellet和Vasseur[18]进一步移除拖曳项。本文首次将趋化力项纳入密度依赖粘性框架,证明$\gamma > 4/3$时解存在。

与趋化性模型的比较:文献[1]在常数粘性下证得$\gamma \geq 4$时解存在,[12]利用Sugiyama不等式扩展至$\gamma > 8/5$,其关键为能量$E_0$与不等式:

$$ \int_\Omega \rho c \, dx \leq \frac{1}{2} \|h(\rho)\|_{L^1(\Omega)} + \frac{1}{4} \|\nabla c\|^2_{L^2(\Omega)} + C_1(\gamma) \|c\|^{C_2(\gamma)}_{L^1(\Omega)}, $$

要求$\gamma > 8/5$。本文通过BD熵不等式将范围推进至$\gamma > 4/3$,与退化Keller-Segel方程临界指数一致。特别地,当$\gamma > 4/3$时,引理3中不等式(23)第二式成立,而$\gamma = 4/3$时退化Keller-Segel方程存在爆破解[24],表明本文结果在某种意义上最优。

与Korteweg模型的比较:Lacroix-Violet和Vasseur[15]证明带Korteweg项的可压缩Navier-Stokes方程对$\gamma > 1$存在弱解,并通过$\kappa \to 0$极限得到$\kappa = 0$时的解,但未考虑趋化力项。本文将此框架扩展至趋化性系统,证明定理1:当$\gamma > 1$时,系统(12)-(13)存在弱解,且当$\gamma > 4/3$时正则性与$(\kappa, r_0, r_1)$无关,允许$\kappa \to 0$。关键区别在于趋化力项引入$K_1$,需更精细估计:当$4/3 < \gamma < 2$时,$\rho^2$由$|\nabla \rho^{\gamma/2}|^2$控制;当$\gamma > 1$时,需Korteweg项提供的$|\nabla^4 \sqrt{\rho}|^4$控制。

时间尺度与适用范围:本文结果适用于任意有限时间$T > 0$,相比[1,12]将$\gamma$范围从$\gamma > 8/5$推进至$\gamma > 4/3$,引入Korteweg项后覆盖$\gamma > 1$。特别地,$\gamma > 4/3$与退化Keller-Segel方程的全局存在性条件完全匹配,表明趋化力项是限制$\gamma$范围的主要因素,而非粘性项的密度依赖性。

代价与局限:1. 为处理趋化力项,必须结合BD熵与自由能估计,分析复杂度显著增加,需对$\gamma$分段讨论($\gamma \geq 2$,$4/3 < \gamma < 2$,$\gamma > 1$);2. 当$\gamma \leq 4/3$时,无法证明原始系统弱解存在,与退化Keller-Segel方程爆破结果一致,表明$\gamma > 4/3$可能为最优条件;3. 引入Korteweg项虽允许覆盖$\gamma > 1$,但最终需取$\kappa \to 0$极限,限制实际应用中$\gamma$范围为$\gamma > 4/3$;且Korteweg项在生物模型中物理意义不明,可能仅为数学工具。

与退化Keller-Segel方程的联系:本文明确指出,当惯性和粘性为零时,系统退化为退化Keller-Segel方程$\partial_t \rho + \zeta \text{div}(\nabla p(\rho) - \rho \nabla c) = 0$。文献[21]证明$\gamma > 4/3$时全局经典解存在,[24]证明$\gamma = 4/3$时存在爆破解。本文结果表明,在完整流体模型中,$\gamma > 4/3$同样是全局弱解存在的临界条件,揭示趋化力项主导解行为,而非粘性项的密度依赖性,为理解血管形成等生物过程的数学模型提供新视角。


2. 接触型Hamilton-Jacobi方程Lyapunov稳定性的偏微分方程刻画

A PDE formulation of Lyapunov stability for contact-type Hamilton-Jacobi equations

作者: Panrui Ni, Jun Yan

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24329v1

摘要翻译: 我们研究紧致流形上接触型Hamilton-Jacobi方程稳态解的Lyapunov稳定性。以往工作通常假设Hamilton量为$C^3$ Tonelli型,并借助Mather测度来刻画稳定性。本文考虑连续、凸且强制的Hamilton量,建立了可验证的偏微分方程型判据,分别判定稳定性与不稳定性。特别地,原先涉及Mather测度的动力学条件被替换为仅依赖于Hamilton量临界值及粘性次解的条件。该结果为稳定性分析提供了基于偏微分方程的框架,并揭示了其与粘性解各类渐近行为之间的联系。

Comment: 33 pages

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

1. 研究对象和背景

本文研究紧致流形上接触型Hamilton-Jacobi方程稳态解的李亚普诺夫稳定性问题。研究对象是以下初值问题:

$$ \begin{cases} u_t(x, t) + H(x, Du(x, t), u(x, t)) = 0 & \text{for } (x, t) \in M \times (0, +\infty), \\ u(x, 0) = \varphi(x) & \text{for } x \in M, \end{cases} \tag{eHJ} $$

其中 $M$ 是紧致、连通、无边界的光滑流形,$H: T^*M \times \mathbb{R} \to \mathbb{R}$ 是Hamiltonian,$\varphi \in C(M)$。所有解均在粘性解意义下理解。稳态解满足:

$$ H(x, Du(x), u(x)) = 0 \quad \text{for } x \in M. \tag{sHJ} $$

定义向后半群 $T^-_t: C(M) \to C(M)$ 使得 $(x, t) \mapsto T^-_t \varphi(x)$ 给出(eHJ)的解,则稳态解对应于 $T^-_t$ 的不动点。李亚普诺夫稳定性的定义如下:若对任意 $\varepsilon > 0$,存在 $\delta > 0$ 使得对所有满足 $\|\varphi - u\|_\infty \leq \delta$ 的 $\varphi \in C(M)$,有 $\|T^-_t \varphi - u\|_\infty < \varepsilon$ 对所有 $t > 0$ 成立,则称解 $u$ 是李亚普诺夫稳定的;若此外还有 $\lim_{t \to +\infty} \|T^-_t \varphi - u\|_\infty = 0$,则称 $u$ 是渐近稳定的。

已有工作(如文献[27])通常假设 $C^3$ Tonelli Hamiltonian(即关于 $p$ 严格凸且超线性的Hamiltonian),并用Mather测度来刻画稳定性。当 $H$ 关于 $u$ 非单调时,解的大时间行为变得复杂,例如在位错动力学模型中会出现这种情况。本文考虑更一般的连续、凸且强制的Hamiltonian,具体形式为:

$$ H(x, p, u) = G(x, p) + W(x, u), $$

满足条件:(A1) $G(x, p)$ 连续、关于 $p$ 凸,且 $\lim_{|p| \to +\infty} \inf_{x \in M} G(x, p) = +\infty$;(A2) $W(x, u)$ 及其偏导数 $\partial_u W(x, u)$ 连续,且存在常数 $\Lambda > 0$ 使得 $|\partial_u W(x, u)| \leq \Lambda$ 对所有 $(x, u) \in M \times \mathbb{R}$ 成立。这种假设自然地出现在常微分方程视角中,其中 $\partial_u W$ 的连续性和有界性类似于标准Lipschitz条件,确保解的唯一性。

本文试图填补的空白是:建立基于PDE的稳定性判据,避免使用接触Hamiltonian流和Mather测度等动力系统工具。具体而言,用Hamiltonian的临界值 $c(G)$ 和粘性次解替代涉及Mather测度的动力学条件,从而提供纯PDE视角的稳定性分析框架。这不仅推广了[27]的结果,还揭示了Lyapunov稳定性与粘性解大时间行为、vanishing discount问题以及均匀化问题之间的深刻联系。

2. 主要定理或主要结果

在假设(A1)-(A2)下,定义 $c(\varepsilon) := c(H(x, p, u_-(x) + \varepsilon))$,其中 $u_- \in C(M)$ 是(sHJ)的解,$c(\cdot)$ 表示相应Hamiltonian的临界值。主要结果如下:

Proposition 1.2:$c(\varepsilon)$ 是Lipschitz连续的,且左导数 $D^-c(\varepsilon)|_{\varepsilon=0}$ 存在。考虑以下条件:

  • (A3)₁:存在 $\zeta > 0$ 使得 $c(H(x, p, u_-(x)) - \zeta \cdot \partial_u W(x, u_-(x))) < 0$;
  • (A3)₂:存在 $\zeta > 0$ 和Lipschitz连续函数 $w(x)$ 使得 $H(x, Dw, u_-(x)) - \zeta \cdot \partial_u W(x, u_-(x)) < 0$ a.e. $x \in M$;
  • (A3)₃:$D^-c(\varepsilon)|_{\varepsilon=0} > 0$。

则(A3)₁与(A3)₂等价,且均蕴含(A3)₃。

Theorem 1.3:若 $u_-$ 满足(A1)-(A2)和(A3)₁-(A3)₃之一,则 $u_-$ 是局部渐近稳定的。具体而言,存在 $\Delta > 0$,使得对任意 $\varphi \in C(M)$ 满足 $\|\varphi - u_-\\_\infty \leq \Delta$,有

$$ \limsup_{t \to \infty} \frac{\ln \|u(\cdot, t) - u_-(\cdot)\|_\infty}{t} \leq -A $$

对某常数 $A > 0$(依赖于 $H$)成立,其中 $u(x, t)$ 是(eHJ)的解。特别地,$\|u(\cdot, t) - u_-(\cdot)\|_\infty = O(e^{-At})$。

Example 1.4:对 $\zeta, \theta > 0$ 和 $\phi \in C^1(M)$,考虑方程

$$ -(|D\phi|^2 - \theta)u + \zeta|Du|^2 + (|D\phi|^2 - \theta)\phi - \zeta|D\phi|^2 = 0, $$

则 $u_- = \phi$ 是解且局部渐近稳定。

Proposition 1.5:考虑条件:

  • (A4)₁:存在 $\zeta > 0$ 使得 $c(H(x, p, u_-(x)) + \zeta \cdot \partial_u W(x, u_-(x))) < 0$;
  • (A4)₂:存在 $\zeta > 0$ 和Lipschitz连续函数 $w(x)$ 使得 $H(x, Dw, u_-(x)) + \zeta \cdot \partial_u W(x, u_-(x)) < 0$ a.e. $x \in M$;
  • (A4)₃:$D^+c(\varepsilon)|_{\varepsilon=0} < 0$。

则(A4)₁与(A4)₂等价,且均蕴含(A4)₃。

Theorem 1.6:若 $u_-$ 满足(A1)-(A2)和(A4)₁-(A4)₃之一,则 $u_-$ 是李亚普诺夫不稳定的。具体而言,存在 $\Delta > 0$,使得对任意 $\varepsilon \in (0, 1)$,存在 $\varphi_\varepsilon \in C(M)$ 满足 $\|\varphi_\varepsilon - u_-\\_\infty \leq \varepsilon$ 且

$$ \limsup_{t \to +\infty} \|T^-_t \varphi_\varepsilon - u_-\\_\infty \geq \Delta. $$

Proposition 1.7:在(A1)-(A2)和(A1)’($p \mapsto G(x, p)$ 严格凸且 $G(x, p) \geq \Theta(|p|)$,$\Theta$ 超线性)下,考虑条件:

  • (A5)₁:对(sHJ)的任意次解 $w(x)$,存在 $\zeta > 0$ 使得 $c(H(x, p, w(x)) - \zeta \cdot \partial_u W(x, w(x))) < 0$;
  • (A5)₂:对(sHJ)的任意次解 $w(x)$,存在 $\zeta > 0$ 和Lipschitz连续函数 $\tilde{w}(x)$ 使得 $H(x, D\tilde{w}, w(x)) - \zeta \cdot \partial_u W(x, w(x)) < 0$ a.e. $x \in M$。

则(A5)₁与(A5)₂等价。

Theorem 1.8:若满足(A1)-(A2)、(A1)‘和(A5)₁-(A5)₂之一,则(sHJ)至多有一个解。若解 $u_-$ 存在,则它是全局渐近稳定的,即对所有 $\varphi \in C(M)$,有 $\lim_{t \to +\infty} u(x, t) = u_-(x)$ 一致成立。

Corollary 1.9:设 $a(x) \in C(M)$,$a(x) \geq 0$ 且 $a(x) > 0$ 在投影Aubry集 $A$ 上成立,则方程

$$ a(x)u + G(x, Du) = c(G) $$

有唯一全局渐近稳定的解。

Proposition 1.10:对满足条件(1)-(4)的 $H \in BUC(\mathbb{R}^n \times \mathbb{T}^n \times B(0, R) \times \mathbb{R})$,设 $u_\varepsilon$ 是

$$ H\left(x, \frac{x}{\varepsilon}, Du_\varepsilon, u_\varepsilon\right) = 0 $$

在 $BUC(\mathbb{R}^n)$ 中的唯一解。定义有效Hamiltonian $\overline{H}(x, p, c)$ 为方程 $H(x, y, p + D_y w, c) = \overline{H}(x, p, c)$ 的解,则存在唯一解 $u \in BUC(\mathbb{R}^n)$ 满足 $\overline{H}(x, Du, u) = 0$,且

$$ \|u_\varepsilon - u\|_\infty \leq C\sqrt{\varepsilon} $$

对某常数 $C > 0$(仅依赖于 $H$)成立。

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

主要技术障碍源于Hamiltonian的低正则性(仅假设连续而非 $C^3$),这阻碍了接触Hamiltonian流的构造,特别是无法得到最小化曲线速度的估计。当 $H$ 关于 $u$ 非单调时,解的大时间行为更加复杂,无法像单调情况那样简单收敛到稳态解。

作者的核心创新在于将涉及Mather测度的动力学条件替换为基于Hamiltonian临界值和粘性次解的PDE类型条件。具体技术包括:

  1. 临界值函数 $c(\varepsilon)$ 的微分分析: 引入 $c(\varepsilon) := c(H(x, p, u_-(x) + \varepsilon))$,证明其Lipschitz连续性及单侧导数的存在性(Proposition 1.2和1.5)。关键引理Lemma 3.1表明,在(A3)₁下,对与 $H_-(x, p) = H(x, p, u_-(x))$ 相关的任意Mather测度 $\tilde{\mu}_-$,有

    $$ \int_{TM} \partial_u W(x, u_-(x)) d\tilde{\mu}_- > A > 0. $$

    证明利用了Mather测度的最小化性质(2.2):

    $$ \min_{\tilde{\mu}} \int_{TM} L(x, v) d\tilde{\mu} = -c(G), $$

    以及 $u_-$ 作为(sHJ)解的性质 $\int_{TM} (L(x, v) - W(x, u_-(x))) d\tilde{\mu}_- = 0$。这建立了Mather测度与PDE条件之间的联系,避免了显式构造Hamiltonian流。

  2. 解半群的精细估计:

  • 向后半群的表示:
    $$ T^-_t \varphi(x) = \inf \left\{ \varphi(\gamma(0)) + \int_0^t L_H(\gamma(s), \dot{\gamma}(s), T^-_s \varphi(\gamma(s))) ds \right\}, $$
    其中下确界取遍满足 $\gamma(t) = x$ 的绝对连续曲线(Proposition 2.6)。
  • 半群的比较原理:若 $\varphi_1 < \varphi_2$,则 $T^-_t \varphi_1 < T^-_t \varphi_2$(Proposition 2.7(1))。
  • 指数型Lipschitz估计:$\|T^-_t \varphi_1 - T^-_t \varphi_2\|_\infty \leq e^{\Lambda t} \|\varphi_1 - \varphi_2\|_\infty$(Proposition 2.7(2)),这是处理非单调项 $W(x, u)$ 的关键。
  • 等Lipschitz连续性:在(A1)‘下,若 $T^-_t \varphi$ 一致有界,则族 $\{T^-_t \varphi\}_{t \geq \delta}$ 等Lipschitz连续(Proposition 2.13),这对全局稳定性证明至关重要。
  1. 稳定性证明框架:
  • 对Theorem 1.3,利用(A3)₃和半群的指数估计,证明存在 $\Delta > 0$ 和 $A > 0$ 使得对 $\|\varphi - u_-\\_\infty \leq \Delta$,有
    $$ \|T^-_t \varphi - u_-\\_\infty \leq \|\varphi - u_-\\_\infty e^{-At}. $$
    这里关键步骤是将 $c(\varepsilon)$ 的导数与解的衰减速率联系起来,通过构造适当的次解和超解来控制解的演化。
  • 对Theorem 1.6,利用(A4)₁-(A4)₃构造初始扰动 $\varphi_\varepsilon$,使得 $T^-_t \varphi_\varepsilon$ 在长时间后偏离 $u_-$ 至少 $\Delta$。证明中利用了半群的单调性和临界值条件的符号性质。
  1. 全局稳定性与唯一性:
  • Theorem 1.8的证明结合了等Lipschitz连续性、半群的单调性以及条件(A5)₁-(A5)₂。关键观察是:对任意两个次解 $w_1, w_2$,条件(A5)保证了 $T^-_t w_1$ 和 $T^-_t w_2$ 会收敛到同一解。这里利用了向前半群 $T^+_t$ 与向后半群的关系 $T^+_t \circ T^-_t \varphi \leq \varphi \leq T^-_t \circ T^+_t \varphi$(Proposition 2.12)。
  • Corollary 1.9通过将条件转化为投影Aubry集上的正性条件,简化了稳定性判据。
  1. 均匀化问题的收敛率:
  • Proposition 1.10将大时间行为结果应用于多尺度问题。证明中结合了曲线切割方法(curve-cutting approach)和大时间行为的指数衰减估计。关键创新是利用 $\|u(\cdot, t) - u_-\\_\infty = O(e^{-At})$ 来控制多尺度解 $u_\varepsilon$ 与有效解 $u$ 之间的差,最终得到 $\|u_\varepsilon - u\|_\infty \leq C\sqrt{\varepsilon}$。与线性情况不同,这里无法使用重缩放参数,因此 $\sqrt{\varepsilon}$ 的收敛率是最佳可能。

4. 与之前工作的比较

与既有工作的比较主要体现在以下几个方面:

  1. 与[27]的对比:
  • [27]假设 $C^3$ Tonelli Hamiltonian(即 $H$ 关于 $p$ 严格凸且超线性),并用Mather测度刻画稳定性。具体而言,[27]证明当 $\int_{TM} \partial_u H d\tilde{\mu} > 0$ 对所有Mather测度 $\tilde{\mu}$ 成立时,稳态解是李亚普诺夫稳定的。
  • 本文将假设放宽至连续、凸且强制的Hamiltonian(仅需(A1)-(A2)),并用PDE条件替代Mather测度条件。例如,Theorem 1.3中的(A3)₁等价于 $c(H(x, p, u_-(x)) - \zeta \cdot \partial_u W(x, u_-(x))) < 0$,这可以通过计算临界值直接验证,无需构造接触Hamiltonian流。
  • 适用范围:本文结果适用于更广泛的Hamiltonian类(如 $G(x, p)$ 可能不是 $C^2$ 的情况),但代价是仅处理了可分离形式 $H(x, p, u) = G(x, p) + W(x, u)$。一般Hamiltonian的处理留待未来工作。
  1. 与[7, 9]的联系:
  • [7, 9]研究vanishing discount问题,使用条件 $\int_{TM} \partial_u H d\tilde{\mu} > 0$(对所有Mather测度 $\tilde{\mu}$)来保证解的收敛性和唯一性([7, Theorems 1.7, 1.8]和[9, Theorem 3.1])。
  • 本文揭示了这些条件与Lyapunov稳定性的等价性:Lemma 3.1表明(A3)₁蕴含 $\int_{TM} \partial_u W d\tilde{\mu}_- > A > 0$,这正是[7, 9]中的关键假设。但本文通过临界值 $c(\varepsilon)$ 提供了PDE视角的解释,避免了对Mather测度的显式依赖。
  • 本文进一步将条件推广至局部稳定性(Theorem 1.3)和全局稳定性(Theorem 1.8),而[7, 9]主要关注vanishing discount问题。
  1. 与单调情形的比较:
  • 当 $H$ 关于 $u$ 非递减时(即 $\partial_u H \geq 0$),[28]证明解 $u(x, t)$ 一致收敛到稳态解。
  • 本文处理了更一般的非单调情形($\partial_u H$ 可正可负),并给出了稳定性/不稳定性判据。例如,Corollary 1.9中 $a(x) > 0$ 在Aubry集上成立的条件,比[28]中 $a(x) \geq 0$ 更弱(仅需在Aubry集上严格正)。
  • 代价是:当 $\partial_u H$ 变号时,解可能不具有等连续性(文中指出这是自然问题,但未解决)。
  1. 与均匀化理论的联系:
  • [10]对非线性均匀化问题仅获得收敛性(无速率);[14]对线性问题获得 $O(\varepsilon)$ 收敛率;[19]处理了单尺度非线性问题。
  • 本文Proposition 1.10对多尺度非线性问题获得 $O(\sqrt{\varepsilon})$ 收敛率,这是现有技术下的最佳结果。关键区别在于:多尺度性质阻止了使用重缩放参数和比较原理,因此无法达到 $O(\varepsilon)$ 速率。
  • 创新点在于结合大时间行为的指数衰减估计(Theorem 1.3)与曲线切割方法,这比[21, Theorem 1.1]的纯曲线切割方法更精细。
  1. 与弱KAM理论的融合:
  • 传统弱KAM理论(如[11, 12])主要处理 $H = H(x, p)$ 的情况,通过Aubry集和Mather测度研究解的大时间行为。
  • 本文将弱KAM框架扩展至接触型Hamilton-Jacobi方程($H = H(x, p, u)$),定义了投影Aubry集 $A$(公式(2.1))和Mather测度(Proposition 2.4),并建立了它们与Lyapunov稳定性的联系。
  • 特别地,Theorem 1.8表明当 $a(x) > 0$ 在Aubry集上成立时,方程 $a(x)u + G(x, Du) = c(G)$ 有唯一解,这推广了弱KAM理论中关于解唯一性的经典结果。

局限性与推进:

  • 本文仅处理了可分离形式的Hamiltonian,而[27]处理更一般的接触Hamiltonian系统。作者明确指出一般Hamiltonian的处理是未来工作方向。
  • 当 $H$ 关于 $u$ 非单调时,解的正则性问题(如等连续性)仍未解决,文中指出这是自然问题但未提供答案。
  • 二阶方程和弱耦合系统的情形未被涵盖,作者在引言中提到这些扩展将是未来工作。
  • 收敛率 $O(\sqrt{\varepsilon})$ 可能不是最优的,但文中论证了在多尺度非线性情形下难以改进。

总体而言,本文通过PDE方法重构了Lyapunov稳定性理论,避免了动力系统工具的复杂性,同时揭示了与大时间行为、vanishing discount问题和均匀化问题的深刻联系。这种PDE视角不仅放宽了正则性假设,还提供了更直接的稳定性判据,对应用数学和动力系统理论均有重要价值。


3. 关于小阿尔芬数下不可压磁流体动力学方程的整体解

On Global-in-time Solutions of Incompressible MHD Equations with Small Alfvén Numbers

作者: Fei Jiang, Xiao Ren, Yi Zhou

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23735v1

摘要翻译: 1965年,Kraichnan指出:足够强的背景磁场(即小阿尔芬数情形)将削弱非线性相互作用,并抑制磁流体动力学(简称MHD)系统中强梯度的形成;该系统满足 $\mu = \nu \geqslant 0$,其中 $\mu$ 与 $\nu$ 分别为运动粘性系数与电阻率系数。这意味着:当 $\mu = \nu \geqslant 0$ 时,MHD系统存在具有小阿尔芬数的整体大扰动解。此类大扰动解的存在性最早于1988年由Bardos–Sulem–Sulem在Hölder空间中对 $\mu = \nu = 0$ 情形给出了严格数学验证;而对 $\mu = \nu > 0$ 情形,Cai–Cui–Jiang–Liu最近在Sobolev空间中完成了相应验证。本文进一步将该结论推广至更一般的“$\mu > 0$ 且 $\nu > 0$”情形,并通过发展一种新方法给出严格证明;该方法的核心是一个用于处理非线性相互作用项的关键双线性估计。此外,本文还同时建立了两个附加结果:非线性相互作用项的消失行为,以及解在小阿尔芬数极限下的收敛性。

Comment: 无

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

1. 研究对象和背景

本文研究不可压缩磁流体动力学(MHD)方程在小Alfvén数条件下的全局解问题。原始MHD系统为:

$$ \begin{cases} \rho v_t + \rho v \cdot \nabla v - \tilde{\mu}\Delta v + \nabla P = \lambda M \cdot \nabla M/(4\pi), \\ M_t + v \cdot \nabla M - \nu \Delta M = M \cdot \nabla v, \\ \text{div}\, v = \text{div}\, M = 0. \end{cases} $$

其中$v$、$M$和$P$分别表示速度、磁场和压力,$\rho > 0$和$\lambda$为常数,$\tilde{\mu} \geqslant 0$和$\nu \geqslant 0$分别表示运动粘性和电阻系数。考虑静止状态$(0, \bar{M})$,其中$\bar{M} = \varpi e_1$为背景磁场,$e_1 \in \mathbb{R}^n$为单位向量,通过变换:

$$ H = \sqrt{\frac{\lambda}{4\pi\rho}}(M - \bar{M}), \quad p = \frac{P}{\rho}, \quad \mu = \frac{\tilde{\mu}}{\rho}, \quad \varepsilon = \sqrt{\frac{4\pi\rho}{\lambda\varpi^2}}, $$

得到扰动系统:

$$ \begin{cases} v_t + v \cdot \nabla v - \mu\Delta v + \nabla p = (\varepsilon^{-1}e_1 + H) \cdot \nabla H, \\ H_t + v \cdot \nabla H - \nu\Delta H = (\varepsilon^{-1}e_1 + H) \cdot \nabla v, \\ \text{div}\, v = \text{div}\, H = 0. \end{cases} $$

Kraichnan于1965年指出,当$\mu = \nu \geqslant 0$时,足够强的背景磁场(即小$\varepsilon$)会减弱非线性相互作用并抑制强梯度形成。Bardos-Sulem-Sulem(1988)在Hölder空间中首次数学验证了$\mu = \nu = 0$的情况,Cai-Cui-Jiang-Liu近期在Sobolev空间中验证了$\mu = \nu > 0$的情况。本文试图填补$\mu > 0$且$\nu > 0$这一一般情况下的理论空白,解决如下断言:若初始扰动满足$\|(v_0, H_0) \|_X \ll \varepsilon^{-s}$($s > 0$),则系统存在全局强解。该问题的重要性在于:一方面,它验证了Kraichnan关于强磁场抑制湍流的物理直觉;另一方面,它揭示了磁场对MHD系统非线性结构的深层影响,特别是当粘性与电阻系数不同时的数学机制。已有工作主要集中在$\mu = \nu$的特殊情况,而一般情形$\mu \neq \nu$的全局解问题仍是开放的,这正是本文试图推进的理论空白。

2. 主要定理或主要结果

在$n = 2, 3$、$\mu > 0$、$\nu > 0$、$(v_0, H_0) \in H^m_\sigma$($m \geqslant 3$)的假设下,若小Alfvén数条件$\varepsilon \max\{\|(v_0, H_0) \|_m^2, \|(v_0, H_0) \|_m^3\} \leqslant \delta$成立($\delta > 0$为仅依赖$\mu, \nu, n, m$的小常数),则:(1)定理2.1证明Cauchy问题(1.6)存在唯一全局解$(u, h) \in C^0(\mathbb{R}^+_0, H^m \times H^m)$,且满足$\|(u, h)(t) \|_m + \sqrt{\varepsilon}\|\nabla(u, h) \|_{L^2_t H^m} \lesssim \varepsilon\|(v_0, H_0) \|_m$;(2)定理2.2给出非线性解与线性解$(u^L, h^L)$的误差估计:$\|\nabla(u^d, h^d)(t) \|_{m-2} + \sqrt{\varepsilon}\|\nabla^2(u^d, h^d) \|_{L^2_t H^{m-2}} \lesssim \varepsilon^{5/4}(\|(v_0, H_0) \|_m^{3/2} + \|(v_0, H_0) \|_m^{7/4})$,其中$(u^d, h^d) = (u, h) - (u^L, h^L)$;(3)定理2.3表明原始系统(1.3)的全局解满足稳定性估计$\|(v, H)(t) \|_m + \|\nabla(v, H) \|_{L^2_t H^m} \lesssim \|(v_0, H_0) \|_m$,且非线性相互作用消失现象:$\|\nabla(v^d, H^d)(t) \|_{m-2} + \|\nabla^2(v^d, H^d) \|_{L^2_t H^{m-2}} \lesssim \varepsilon^{1/4}(\|(v_0, H_0) \|_m^{3/2} + \|(v_0, H_0) \|_m^{7/4})$;(4)定理2.5证明小Alfvén数极限:当$\varepsilon \to 0$时,解$(v^\varepsilon, H^\varepsilon) \to 0$,且在平板$S_R$上满足$\|(v^\varepsilon, H^\varepsilon)(t) \|_{L^\infty(S^l_R)} \lesssim_{\theta,R,t_0,l} \varepsilon^{n/16}(\|(v_0, H_0) \|_m^{1+n/8} + \|(v_0, H_0) \|_m^{1+3n/16} + \|(v_0, H_0) \|_{L^1})$。特别地,当$\|(v_0, H_0) \|_m \geqslant 1$时,小Alfvén数条件简化为$\varepsilon\|(v_0, H_0) \|_m^3 \leqslant \delta$,即断言1.1中$s = 1/3$。

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

核心障碍在于标准能量方法无法处理非线性积分。例如,对$|\alpha| = m$,需估计$I := \int_0^T \int \partial^\alpha h \cdot \nabla h \cdot \partial^\alpha u \, dx dt$。标准估计给出$|I| \lesssim \varepsilon K \|\nabla(u, h) \|_{L^2_T H^{m-1}}^2$,但建立稳定性估计(2.4)需要将$\varepsilon$替换为$\varepsilon^\vartheta$($\vartheta > 1$)。作者的创新点是引入$\Lambda^\pm := u \pm h$重写非线性项,并利用div-curl引理(Lemma 3.4)推导关键双线性估计:

$$ \left\|\|\nabla\Lambda^+\|_{H^{m-2}_{x_v}} \|\nabla\Lambda^-\|_{H^{m-1}_{x_v}} + \|\nabla\Lambda^-\|_{H^{m-2}_{x_v}} \|\nabla\Lambda^+\|_{H^{m-1}_{x_v}}\right\|_{L^2_{T,x_1}} \lesssim \sup_{0 \leqslant t \leqslant T} \{\|(u, h)(t) \|_m^2\} + \sup_{0 \leqslant t \leqslant T} \{\|(u, h)(t) \|_m\} (\sup_{0 \leqslant t \leqslant T} \{\|(u, h)(t) \|_m\} + \varepsilon)^{1/2} \|\nabla(u, h) \|_{L^2_T H^m} $$

该估计将非线性积分控制为$\sup_{0 \leqslant t \leqslant T} \{\|(u, h)(t) \|_m\} \left( \sup_{0 \leqslant t \leqslant T} \{\|(u, h)(t) \|_m\} + \varepsilon^{1/2}(K^{1/2} + 1)\|\nabla(u, h) \|_{L^2_T H^m} \right) \|\nabla(u, h) \|_{L^2_T H^m}$。在小Alfvén数条件$\varepsilon^{1/2} \max\{K, K^{3/2}\} \leqslant \delta \ll 1$下,非线性项可被能量项$E(T)$和耗散项$D(T)$吸收,从而得到稳定性估计$E(T) + D(T) \leqslant (c\varepsilon\|(v_0, H_0) \|_m)^2$。此外,作者通过分析线性基本解的振荡结构(定理2.4),结合误差估计(2.16)和Nirenberg插值不等式$\|f\|_{L^\infty} \lesssim \|f\|_0^{1-n/4} \|\nabla^2 f\|_0^{n/4}$,严格证明了小Alfvén数极限行为。证明框架的关键在于:(1) 通过$\Lambda^\pm$分解揭示非线性项的特殊结构;(2) 利用横截方向$H^{m-k}_{x_v}$范数的乘积估计克服标准方法中$\varepsilon$阶数不足的问题;(3) 将全局解存在性问题转化为对非线性相互作用的精细估计。

4. 与之前工作的比较

与Bardos-Sulem-Sulem(1988)相比:本文处理的是$\mu > 0$且$\nu > 0$的情形,而Bardos等仅研究$\mu = \nu = 0$的理想MHD系统;Bardos等利用双曲结构在Hölder空间中证明全局适定性,而本文在Sobolev空间中建立解的稳定性估计$\|(v, H)(t) \|_m + \|\nabla(v, H) \|_{L^2_t H^m} \lesssim \|(v_0, H_0) \|_m$,且无需额外正则性条件(如$\||\nabla|^{-1}(v_0, H_0) \|_0 < \infty$)。与Cai-Cui-Jiang-Liu相比:Cai等仅适用于$\mu = \nu > 0$的对称情形,其方法依赖Alinhac ghost weight技术,而本文通过$\Lambda^\pm$分解和双线性估计处理$\mu \neq \nu$的非对称情形;Cai等在证明中需要$\||\nabla|^{-1}(v_0, H_0) \|_0 < \infty$的附加条件,而本文定理2.3完全移除了这一限制;在小Alfvén数极限方面,Cai等仅证明$\mu = \nu$时解收敛到零,但无法推广到$\mu \neq \nu$,而本文通过误差估计(2.16)和线性解的振荡分析,严格证明了一般情形的极限行为。与Zhang、Jiang-Jiang等研究$\mu > 0$且$\nu = 0$的工作相比:这些工作利用波动方程或Lagrange坐标下的磁抑制理论,而本文方法同时依赖粘性和电阻效应,故无法处理$\mu = 0$且$\nu > 0$的情形(需进一步研究);本文的误差阶$\varepsilon^{1/4}$在$\mu, \nu > 0$时优于$\nu = 0$情形的估计。适用范围方面:本文结果适用于$\mathbb{R}^n$($n = 2, 3$)上的大扰动解,而Cai等的结果仅限于周期域;代价是本文方法要求$\mu, \nu > 0$,无法覆盖理想MHD($\mu = \nu = 0$)或仅电阻耗散的情形($\mu = 0, \nu > 0$)。时间尺度上,本文证明了任意时间$T > 0$的全局解存在性,且解在$\varepsilon \to 0$时以$\varepsilon^{n/16}$速率收敛,这比Cai等的收敛速率更精确(后者仅给出逐点收敛)。


4. 关于多维可压Navier-Stokes-Korteweg方程组的Cauchy问题:具任意大初值的全局强解

On the Cauchy problem for the multi-dimensional compressible Navier-Stokes-Korteweg system: Global strong solutions with arbitrarily large initial data

作者: Xiangdi Huang, Muxi Lei, Huitao Zhou

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23652v1

摘要翻译: 自Korteweg(1901)的开创性工作及Dunn与Serrin(1985)对毛细流体模型的后续完善以来,多维可压Navier-Stokes-Korteweg(NSK)方程组在初值任意大的条件下全局强解的存在性,一直是流体力学中一个极具挑战性的公开问题。该问题最近由[Gu-Huang-Meng-Zhou, arXiv:2603.11762]部分解决:他们在周期域$\mathbb{T}^N$($N=2,3$)上,在粘性系数满足BD型代数关系($μ(ρ) = νρ^α$, $λ(ρ) = 2ν(α-1)ρ^α$)且Korteweg应力张量满足广义Bohm恒等式($κ(ρ) = \varepsilon^2 α^2 ρ^{2α-3}$)的条件下,建立了具任意大初值的全局强解的存在性。然而,在相同条件下,Cauchy问题(即全空间$\mathbb{R}^N$上的初值问题)的全局强解存在性仍为未解决问题。本文解决了这一问题:我们证明了在$\mathbb{R}^N$($N=2,3$)上、初值任意大且远场密度非真空的条件下,Cauchy问题全局强解的存在性。通过结合精细的截断分析与原创的修正型Nash-Moser迭代格式,我们克服了全空间中密度缺乏可积性所引发的困难。该结果将可压Navier-Stokes-Korteweg方程的大初值理论从有界环面$\mathbb{T}^N$推广至无界全空间$\mathbb{R}^N$,因而适用于更一般的物理背景。

Comment: 61 pages

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

1. 研究对象和背景

本文研究$\mathbb{R}^N$($N=2,3$)上可压缩Navier-Stokes-Korteweg(NSK)系统的Cauchy问题:

$$ \begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla \rho^\gamma = \text{div}\left(2\mu(\rho)Du\right) + \nabla\left(\lambda(\rho) \text{div} u\right) + \text{div} K, \end{cases} $$

其中$\gamma \geq 1$,$\rho = \rho(x,t)$、$u = (u_1,\ldots,u_N)^\top = u(x,t)$和$P = \rho^\gamma$分别表示流体密度、速度和压力。变形张量$Du = \frac{\nabla u + (\nabla u)^\top}{2}$,Korteweg应力张量$K$定义为:

$$ K = \left(\rho \text{div}(\kappa(\rho)\nabla\rho) - \frac{\rho\kappa'(\rho) - \kappa(\rho)}{2}|\nabla\rho|^2\right)I - \kappa(\rho)\nabla\rho \otimes \nabla\rho, $$

$I$表示$N \times N$单位矩阵。研究对象的特殊性在于其系数结构:

$$ \mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha-1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}, \quad \nu \geq \varepsilon > 0. $$

特别地,$\lambda(\rho) = 2\rho\mu'(\rho) - 2\mu(\rho)$,满足BD代数关系,且物理约束$\mu(\rho) \geq 0$和$2\mu(\rho) + N\lambda(\rho) \geq 0$要求$\alpha \geq \frac{N-1}{N}$。远场平衡态设定为$\rho(x,t) \to 1$,$u(x,t) \to 0$当$|x| \to \infty$,初值满足$0 < \underline{\rho}_0 \leq \rho_0 \leq \overline{\rho}_0$,$\rho_0 - 1 \in H^3(\mathbb{R}^N)$,$u_0 \in H^2(\mathbb{R}^N)$。

该问题的重要性在于,自Korteweg(1901)开创性工作和Dunn-Serrin(1985)对毛细管流体模型的完善以来,具有任意大初值的多维可压缩NSK系统强解的全局存在性一直是流体力学中的重大开放问题。Gu-Huang-Meng-Zhou(arXiv:2603.11762)近期在周期域$\mathbb{T}^N$($N=2,3$)上解决了此问题,但Cauchy问题在这些条件下的全局强解存在性仍属开放。本文填补的关键空白是:在$\mathbb{R}^N$($N=2,3$)上,对任意大初值和非真空远场密度的Cauchy问题建立全局强解理论,将大初值理论从有界环面$\mathbb{T}^N$扩展到无界全空间$\mathbb{R}^N$,适用于更广泛的物理场景。

2. 主要定理或主要结果

在假设1.1条件下,Theorem 1.1证明:对$N \in \{2,3\}$,Cauchy问题(1),(5)-(8)存在唯一全局强解$(\rho,u)$。具体而言,对任意$0 < T < \infty$,存在常数$C(T) > 0$使得:

$$ (C(T))^{-1} \leq \rho(x,t) \leq C(T) \quad \forall (x,t) \in \mathbb{R}^N \times [0,T], $$

且解满足:

$$ \begin{cases} \rho - 1 \in C([0,T];H^3(\mathbb{R}^N)) \cap L^2(0,T;H^4(\mathbb{R}^N)), \\ \rho_t \in C([0,T];H^1(\mathbb{R}^N)) \cap L^2(0,T;H^2(\mathbb{R}^N)), \\ u \in C([0,T];H^2(\mathbb{R}^N)) \cap L^2(0,T;H^3(\mathbb{R}^N)), \\ u_t \in L^\infty(0,T;L^2(\mathbb{R}^N)) \cap L^2(0,T;H^1(\mathbb{R}^N)). \end{cases} $$

其中$C(T)$仅依赖于$T$、系统参数及初值范数。假设1.1要求$(\alpha,\gamma,\beta)$满足:

  • $N=2$时,$\alpha \in \left(\frac{\sqrt{5}-1}{2},1\right)$,$\gamma \in [1,\infty)$,$\beta \in [0,\beta_2^+(\alpha))$;
  • $N=3$时,$\alpha \in \left(\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}},1\right)$,$\gamma \in \left[1,\frac{15\alpha-7}{3}\right)$,$\beta \in [0,\beta_3^+(\alpha))$。

这里$\beta = \sqrt{1 - \frac{\varepsilon^2}{\nu^2}} \in [0,1)$,$\beta_N^+(\alpha)$是二次方程:

$$ \left(\beta_N^+(\alpha)\right)^2 + \frac{2(1-\alpha)(N\alpha^2 + \sqrt{N}\alpha + 1 - N)}{\alpha(1 - \sqrt{N}(1-\alpha))^2}\beta_N^+(\alpha) + \frac{(1-\alpha)(-N\alpha^2 - N\alpha + 2N - 2)}{\alpha(1 - \sqrt{N}(1-\alpha))^2} = 0 $$

的唯一正解。值得注意的是,当$\alpha \to 1$时,$\frac{15\alpha-7}{3} \to \frac{8}{3}$,与Huang-Gu-Lei[13]中量子NSK系统($\alpha=1$)的$\gamma$范围一致。

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

证明中的核心障碍源于全空间$\mathbb{R}^N$上密度缺乏自然可积性,导致周期域方法失效。作者的创新点主要体现在四个关键技术:

§2.3.1 $\gamma$的分段分析($N=3$):为获得有效速度的加权可积性,在$\gamma \in [\frac{3}{2}\alpha - \frac{1}{6}, \frac{15\alpha-7}{3})$范围内分两段处理:(1)当$\gamma \in [\frac{3}{2}\alpha - \frac{1}{6}, \frac{15\alpha-7}{3})$时,取$1 < q < 2$满足$\frac{3}{2}\alpha - \frac{1}{6} \leq \gamma < 3\alpha - 1 + \frac{6\alpha-4}{q+2}$,利用$\|\rho^{\frac{3}{2}\alpha-1 -1}\|_{L^2(0,T;H^2(\mathbb{R}^3))} \leq C$,将压力项按区域$\Omega_1(t) = \{x: \rho(x,t) \leq 4\}$和$\Omega_2(t) = \{x: \rho(x,t) > 4\}$分解;(2)当$\gamma \in [1, \frac{9}{2}\alpha - 2)$时,通过定义$\zeta' = \frac{(5\alpha-3)q+2}{2\alpha-12}$,利用附加可积性条件$\int_0^T \left(\int_{\Omega_2(t)} \rho^{\zeta'}dx\right)^{\frac{2\lambda}{(q+2)\zeta'}} dt \leq C$处理高密度区域。

§2.3.2 全空间截断过程:针对$\rho \to 1$在无穷远的特性,引入截断函数$\eta \in C^\infty([0,\infty))$满足$0 \leq \eta \leq 1$,$\eta(s)=0$当$s \leq 2$,$\eta(s)=1$当$s \geq 4$,定义$f := \eta(\rho)\rho$。相对熵估计保证其支撑集具有有限测度,使迭代可行。例如在建立密度上界时,通过测试适当选择的截断密度函数,仅在有限测度的高密度区域改进$L^p$可积性。

§2.3.3 密度下界的精细Nash-Moser迭代:这是全空间分析最精细部分。令$\tau := \rho^{-1}$,对逆密度截断$f := \eta(\tau)\tau$,将分析定位到真空区域。关键创新在于结合相对熵界(支撑集一致有限测度)与密度上界后的加权估计,设计新的Nash-Moser方案:在$\tau$的奇异区域($\rho \to 0$)中,即使扩散与传输强耦合,迭代仍能闭合。具体地,通过构造反向Hölder不等式,确保迭代收敛,最终得到$\rho(x,t) \geq (C(T))^{-1} > 0$。

§2.3.4 一阶高阶估计中强耦合的解耦:构造大规模Gronwall不等式:

$$ \frac{d}{dt}\left(\|\nabla^2\rho(t)\|_{L^2}^2 + \|\nabla v(t)\|_{L^2}^2 + \|\nabla\rho(t)\|_{L^4}^4\right) + \left(\|\nabla^2 v\|_{L^2}^2 + \|\nabla\Delta\rho\|_{L^2}^2 + \||\nabla\rho||\nabla^2\rho|\|_{L^2}^2\right) \leq R(t). $$

其中$\||\nabla\rho||\nabla^2\rho|\|_{L^2}^2$项提供的耗散抵消右侧非线性项。假设1.1中$\alpha$下界$\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}$恰是确保该耗散为正的阈值。若缺少此项,$\|\nabla^2\rho\|_{L^2}$与$\|\nabla v\|_{L^2}$的系统将强耦合无法解耦,这是三维情况的核心技术难点。

4. 与之前工作的比较

本文与既有工作的比较主要体现在以下维度:

• 与Gu-Huang-Meng-Zhou(arXiv:2603.11762)的周期域结果比较:该工作首次在$\mathbb{T}^N$上对一般系数结构NSK系统证明了任意大初值的全局强解存在性,但Cauchy问题仍是开放问题。本文将相同参数范围($N=2$时$\alpha \in (\frac{\sqrt{5}-1}{2},1)$,$\gamma \in [1,\infty)$;$N=3$时$\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}},1)$,$\gamma \in [1,\frac{15\alpha-7}{3})$)推广到$\mathbb{R}^N$,推进的关键在于克服全空间特有的困难:(1)无穷远点远场行为导致密度缺乏自然可积性;(2)域的非紧致性使周期域的迭代方法失效。代价是技术复杂度显著提升,需开发截断分析与改进的Nash-Moser方案。

• 与Huang-Gu-Lei[13]的量子NSK系统比较:该工作研究了$\alpha=1$的特殊情况(量子NSK系统),在$\mathbb{R}^N$上得到$\gamma \in [1,\frac{8}{3})$。本文中当$\alpha \to 1$时,$\frac{15\alpha-7}{3} \to \frac{8}{3}$,验证了参数连续性。但一般$\alpha$情况更复杂,需处理$\alpha$依赖的耗散结构,特别是$\||\nabla\rho||\nabla^2\rho|\|_{L^2}^2$项的阈值$\alpha > \frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}} \approx 0.865$($N=3$),而量子系统($\alpha=1$)无此限制。

• 与周期域方法的比较:周期域中密度上下界估计依赖域的有限测度和密度可积性,全空间中这些条件不成立。本文通过四个技术突破实现推广:(1)$\gamma$分段分析实现有效速度加权估计;(2)全空间截断过程处理无穷远行为;(3)密度下界的反向Hölder不等式确保迭代收敛;(4)利用$\|\nabla\rho\|_{L^4}^4$生成的耗散解耦高阶估计。特别是密度下界证明,周期域中可直接迭代$\rho^{-1}$,而全空间需先通过截断定位真空区域,再结合支撑集控制与加权估计。

• 与弱解结果的比较:Jüngel[16]、Antonelli-Spirito[3]、Lacroix-Violet-Vasseur[18]等建立了允许真空的弱解,但要求$\nu_0 < \kappa_0$或特定系数结构。本文关注强解(更高正则性),且通过证明密度正下界排除真空,代价是要求更严格的系数关系(BD型代数关系与广义Bohm恒等式)和非真空初值。

• 与一般系数结构结果的比较:Danchin-Desjardins[8]、Bresch-Desjardins-Lin[4]、Haspot[12]等对一般系数结构的研究多依赖初值小性假设。本文在无小性假设下证明全局强解,推进的关键在于利用BD熵结构和有效速度变量,但代价是要求特定系数关系($\mu(\rho)=\nu\rho^\alpha$,$\lambda(\rho)=2\nu(\alpha-1)\rho^\alpha$,$\kappa(\rho)=\varepsilon^2\alpha^2\rho^{2\alpha-3}$)。

• 与一维结果的比较:Jüngel[15]、Chen-Zhao[7]、Tsyganov[19]等在一维情况下取得进展,但高维问题更复杂。特别是当$\alpha \neq 1$时,系统失去有效速度等分析结构,使得全局强解存在性成为长期开放问题。本文首次在$\mathbb{R}^N$($N=2,3$)上解决此问题,适用范围扩展至任意大初值,但参数范围受限(如$N=3$时$\gamma < \frac{15\alpha-7}{3}$)。


5. 陆缘冰模型的分析与数值模拟

Analysis and numerical simulations of a landfast ice model

作者: Felix Brandt, Carolin Mehlmann

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23596v1

摘要翻译: 本文考虑了基于Lemieux等人[27]工作的北极陆缘冰(landfast ice)通用建模框架,该框架专为大规模气候模型设计。该方法扩展了Hibler[18]提出的经典黏塑性海冰模型,后者至今仍是气候科学中模拟大尺度海冰动力学最广泛使用的模型。特别地,陆缘冰指附着于海岸线或搁浅的海冰,因而其运动近乎消失。我们对该陆缘冰模型开展了严格的解析与数值研究。主要解析贡献包括:局部强适定性;在无外力作用且初始数据接近常值平衡解情形下的全局强适定性;以及时间周期解的存在性。作为解析研究的补充,我们进行了数值模拟,揭示了陆缘冰模型与经典黏塑性海冰模型之间若干关键的定性差异。特别是,模拟结果展现出以冰速恒为零为特征的稳态平衡态的形成。这些现象与定理4.1所建立的、在平衡态邻域内的整体时间存在性结果,以及定理5.2所给出的时间周期性结果相一致。解析与数值结果的结合,为陆缘冰动力学的结构、稳定性及长期行为提供了新的理论认识。

Comment: 无

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

1. 研究对象和背景

本文研究北极固定冰(landfast ice)的数学模型,这是附着在海岸线或搁浅的海冰,其运动几乎可以忽略。研究对象是基于Lemieux等人工作的固定冰模型,该模型是对Hibler粘塑性海冰模型的扩展,通过调整海冰拉伸强度和引入基底应力参数化来考虑浅水区域的锚定效应。考虑有界区域 $G \ abla v^ op) = \dot{\varepsilon}' + \frac{1}{2} \text{tr}(\dot{\varepsilon})I$ 是应变率张量,$\zeta$ 是粘度,$P'$ 是压力项,具体表达式为:

$$ \zeta := \frac{P + T}{2 \max(\Delta_P(v), \Delta_{\text{min}})}, \quad \Delta_P(v) := \frac{1}{2} \dot{\varepsilon}' : \dot{\varepsilon}', \quad \Delta_{\text{min}} = 2 \times 10^{-9}, $$

$$ P' := P - T. $$

冰强度 $P$ 和拉伸强度 $T$ 建模为:

$$ P = h P^*\exp(-c^*(1-A)) - T, \quad T = k_t h P^*\exp(-c^*(1-A)) = k_t P, $$

其中 $P^* > 0$ 是冰强度参数,$c^* > 0$ 是冰浓度参数,$k_t \in [0,1]$ 是各向同性拉伸强度参数。为确保粘性闭合 $\Delta_{\text{min}}$ 和塑性区域 $\Delta_P(v)$ 之间的平滑过渡,定义正则化应变率不变量:

$$ \Delta(\dot{\varepsilon}) := \sqrt{\Delta_P(v) + \Delta_{\text{min}}^2}. $$

基底应力 $f_b$ 定义为:

$$ f_b = \\ \begin{cases} 0, & \text{if } h \leq h_{\text{crit}}, \\ k_2\left(-\frac{v}{|v|} + v_0\right)(h - h_{\text{crit}}) \exp(-\alpha_b(1-A)), & \text{if } h > h_{\text{crit}}, \end{cases} $$

其中 $h_{\text{crit}}$ 是临界平均厚度,$k_2$ 是确定最大基底应力的自由参数,$\alpha_b$ 是基底应力冰浓度参数,$v_0$ 是小速度参数。初始条件和边界条件为:

$$ v(0) = v_0, \quad h(0) = h_0, \quad A(0) = A_0, \quad \text{on } G, $$

$$ v = 0, \quad \partial_\nu h = \partial_\nu A = 0, \quad \text{on } \partial G \times (0, T). $$

固定冰在北极海岸线附近扮演重要角色,影响海洋和大气之间的热量和动量交换,对气候系统有重要影响。已有理论主要集中在Hibler的粘塑性海冰模型,这是气候科学中模拟大尺度海冰动力学最常用的模型。Lemieux等人的工作引入了基底应力参数化和海冰拉伸强度的修改,成功再现了北极海岸区域的固定冰空间分布。然而,对固定冰模型的严格数学分析几乎空白,本文首次填补了这一空白,特别关注模型的适定性、稳定性及长期行为。

2. 主要定理或主要结果

在 $p, q \in (1, \infty)$ 和 $\mu \in (\frac{1}{p}, 1]$ 满足 $\frac{1}{2} + \frac{1}{p} + \frac{1}{q} < \mu$ 的条件下,对初始数据 $u_0 = (v_0, h_0, A_0) \in V_\mu$(其中 $V_\mu := \{u = (v, h, A) \in X_{\gamma,\mu} : h > \kappa \text{ 且 } A \in (0, 1+\delta)\}$,$X_{\gamma,\mu}$ 是时间迹空间),证明了以下主要结果:

  1. 局部强适定性(Theorem 3.1):存在 $T' = T'(u_0) \in (0, T]$ 和 $r = r(u_0) > 0$,使得系统 (3.4) 在 $[0, T']$ 上对任意初始值 $u_1 \in B_{X_{\gamma,\mu}}(u_0, r)$ 有唯一解 $u = (v, h, A) \in W^{1,p}_\mu(0, T'; X_0) \cap L^p_\mu(0, T'; X_1) \cap C([0, T']; V_\mu)$,其中 $X_0 := L^q(G)^2 \times L^q(G) \times L^q(G)$,$X_1 := W^{2,q}(G)^2 \cap W^{1,q}_0(G)^2 \times W^{2,q}_N(G) \times W^{2,q}_N(G)$。解具有以下性质:(a) 依赖于初始数据的连续性;(b) 瞬时正则化;(c) 最大存在区间由全局存在、趋近 $V_\mu$ 边界或极限不存在三种情况之一确定。

  2. 全局强适定性(Theorem 4.1):在无外部力(即 $v_a = v_o = 0$)且初始数据接近常数平衡解的条件下,系统存在全局强解。具体而言,对于足够小的 $\|u_0 - \bar{u}\|_{X_{\gamma,\mu}}$(其中 $\bar{u} = (\bar{v}, \bar{h}, \bar{A})$ 是常数平衡解),解 $u(t)$ 在 $t \to \infty$ 时收敛到平衡解,且满足 $\|u(t) - \bar{u}\|_{X_\gamma} \leq C e^{-\omega t} \|u_0 - \bar{u}\|_{X_{\gamma,\mu}}$,其中 $\omega > 0$ 是衰减率。

  3. 时间周期解的存在性(Theorem 5.2):对于足够小的时间周期外力(即 $v_a, v_o \in L^{2p}(\mathbb{R}; L^{2q}(G)^2)$ 且周期为 $T$),系统存在时间周期强解 $u(t+T) = u(t)$。特别地,当 $k_t \equiv 1$ 时,需要正则化参数满足小条件 $\Delta_{\text{min}} < \delta$($\delta > 0$ 足够小),解满足 $\|u(t) - \bar{u}\|_{X_\gamma} \leq \varepsilon$($\varepsilon > 0$ 足够小)。

  4. 数值验证:数值模拟显示固定冰模型会形成以消失冰速为特征的静止平衡状态,这与Theorem 4.1的全局存在结果一致。对于时间周期风力,模拟结果验证了Theorem 5.2中时间周期解的存在性。

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

本文的主要技术障碍和创新点如下:

  1. 快冰算子(Fastice operator)的椭圆性分析:与Hibler模型不同,固定冰模型的流变学通过拉伸强度参数 $k_t$ 进行了调整,导致应力张量形式变化。作者定义了快冰算子 $AFI$:

    $$ [AFI v]_i := -\sum_{j,k,l=1}^2 \frac{\tilde{P}}{2\rho h} \frac{1}{\Delta(\dot{\varepsilon})} \left( S^{kl}_{ij} - \frac{1}{\Delta^2(\dot{\varepsilon})}(S\dot{\varepsilon})_{ik}(S\dot{\varepsilon})_{jl} \right) D_k D_l v_j + \frac{1}{2\rho h \Delta(\dot{\varepsilon})} \sum_{j=1}^2 (\partial_j \tilde{P})(S\dot{\varepsilon})_{ij}, $$

    其中 $\tilde{P} = P + T = (1 + k_t)h P^*\exp(-c^*(1-A))$。通过分析其主部系数 $a^{kl}_{ij}(\dot{\varepsilon}, \tilde{P})$,证明了负线性化快冰算子 $-AFI(u_0)$ 的强椭圆性和参数椭圆性(Proposition 3.2)。关键估计是利用 $\triangle^2(d) = d_I^2 + \frac{1}{e^2}(d_{II}^2 + 4d_{III}^2)$(其中 $d_I = d_{11} + d_{22}$,$d_{II} = d_{11} - d_{22}$,$d_{III} = \frac{d_{12} + d_{21}}{2}$)证明:

    $$ \sum_{i,j,k,l=1}^2 -a^{kl}_{ij} d_{ik} d_{jl} = \frac{\tilde{P}(h_0,A_0)}{2\rho h_0 \Delta^3(\dot{\varepsilon})} \left( \Delta^2(\dot{\varepsilon}) d^\top S d - (d^\top S \dot{\varepsilon})^2 \right) \geq c \Delta^2_{\text{min}} \Delta^2(d), $$

    从而得到强椭圆性条件。

  2. 强法向椭圆性的证明:通过分析边界上的行为,证明了负快冰算子的强法向椭圆性(Lemma 3.3)。关键步骤是验证对 $x \in \partial G$,$|\xi| = |\nu| = 1$ 且 $(\xi|\nu) = 0$,以及 $u, v \in \mathbb{C}^2$,有:

    $$ \text{Re} \left( \sum_{i,j,k,l=1}^2 -a^{kl}_{ij}(u_0)(\xi_l u_j - \nu_l v_j)(\xi_k u_i - \nu_k v_i) \right) \geq 0, $$

    且当 $\text{Im}(u|v) \neq 0$ 时严格大于零。这保证了Lopatinskii-Shapiro条件的满足。

  3. 最大 $L^p$-正则性理论的应用:基于快冰算子的椭圆性性质,证明了其 $L^q$ 实现 $-AFI_D(u_0) + \omega$ 具有有界的 $H^\infty$-演算(Proposition 3.4),从而具有最大 $L^p$-正则性。进一步,考虑完整线性化算子矩阵:

    $$ ALFI(u_0) = \\ \begin{pmatrix} -AFI_D(u_0) & \frac{\partial_h P'(h_0,A_0)}{2\rho h_0} \nabla & \frac{\partial_A P'(h_0,A_0)}{2\rho h_0} \nabla \\ 0 & -d_h \Delta_N & 0 \\ 0 & 0 & -d_A \Delta_N \end{pmatrix}, $$

    证明了 $ALFI(u_0) + \omega$ 具有最大 $L^p$-正则性(Proposition 3.6)。关键技巧是利用对角算子 $A_{\text{diag}} = \text{diag}(-AFI_D(u_0), -d_h \Delta_N, -d_A \Delta_N) + \omega$ 的 $H^\infty$-演算性质,并通过分数幂域的等价性 $D(A^{1/2}_{\text{diag}}) \simeq [X_0, X_1]_{1/2}$ 证明非对角项的相对有界性。

  4. 基底应力 $f_b$ 的非线性估计:对 $f_b$ 的关键估计是:

    $$ \left\| \frac{1}{\rho h_1} f_b(v_1, h_1, A_1) - \frac{1}{\rho h_2} f_b(v_2, h_2, A_2) \right\|_{L^q(G)} \leq C \left( \|u_1\|_{X_{\gamma,\mu}} + \|u_2\|_{X_{\gamma,\mu}} + 1 \right) \|u_1 - u_2\|_{X_{\gamma,\mu}}, $$

    其中利用了 $h_i > \kappa$ 和 $B^{2(\mu-1/p)}_{qp}(G) \hookrightarrow L^\infty(G)$ 的嵌入(由 $\frac{1}{2} + \frac{1}{p} + \frac{1}{q} < \mu$ 保证)。这证明了非线性项的Lipschitz连续性(Lemma 3.7)。

  5. 广义线性化稳定性原理:用于处理全局适定性,特别解决了由于Neumann拉普拉斯算子导致算子矩阵谱包含零点的问题。通过引入新基空间(对 $h$ 和 $A$ 分量考虑空间平均为零的函数空间),处理了 $k_t \equiv 1$ 的特殊情况,此时需要正则化参数满足小条件 $\Delta_{\text{min}} < \delta$。

  6. Arendt-Bu定理的应用:用于证明时间周期解的存在性,该定理将最大周期 $L^p$-正则性刻画为最大 $L^p$-正则性和半群的谱条件。通过处理算子不可逆性,证明了对于足够小的时间周期外力,系统存在时间周期强解。

4. 与之前工作的比较

本文与既有PDE/分析理论文献的比较如下:

  1. 与Hibler模型分析的比较:
  • Hibler的粘塑性海冰模型是气候科学中模拟大尺度海冰动力学的标准模型,但不考虑固定冰的特殊性质。Brandt和Mehlmann [9] 首次对Hibler模型的完全抛物正则化版本进行了局部强适定性和全局强适定性(接近常数平衡解)的严格分析。本文将这一分析框架扩展到固定冰模型,关键区别在于流变学的修改:固定冰模型通过拉伸强度参数 $k_t$ 调整了应力张量,特别是当 $k_t \equiv 1$ 时,分析更为复杂,需要对正则化参数施加小条件 $\Delta_{\text{min}} < \delta$。
  • 与 [9] 相比,本文的线性化算子矩阵 $ALFI(u_0)$ 包含额外的基底应力项 $f_b$,这导致非线性估计更为复杂。特别地,$f_b$ 在 $h = h_{\text{crit}}$ 处不光滑,作者通过仔细分析证明了其Lipschitz连续性(Lemma 3.7),而 [9] 中的非线性项相对更规则。
  • 在全局适定性方面,本文处理了由于Neumann拉普拉斯算子导致算子矩阵谱包含零点的情况,而 [9] 中的Dirichlet边界条件避免了这一问题。本文采用广义线性化稳定性原理 [40] 来处理这一困难,扩展了分析框架的适用范围。
  1. 与其他海冰模型分析的比较:
  • [29] 建立了Hibler模型不同正则化的局部强适定性,但使用了不同的应力张量正则化。本文的快冰算子分析基于 [9] 的方法,但需要处理 $k_t$ 参数带来的额外复杂性。
  • [8] 使用拉格朗日方法研究了Hibler的抛物-双曲模型,强调了平衡定律的双曲性。本文主要关注抛物正则化模型,但也讨论了抛物-双曲设置(Section 3.3),指出在无扩散项情况下,局部适定性证明需要额外的技术处理。
  • [6] 研究了带有Voigt正则化的EVP海冰模型的全局适定性。与本文相比,EVP模型通过引入弹性项降低了计算成本,但数学结构不同:EVP模型是二阶时间导数系统,而本文是准线性抛物系统。
  • [16, 14] 研究了Hibler模型动量方程的弱解存在性,特别是处理了应力张量的奇异极限。本文关注强解,要求更高的正则性,但提供了更详细的解结构和稳定性信息。
  1. 时间尺度和长期行为的推进:
  • 本文证明了在无外部力且初始数据接近常数平衡解的情况下,固定冰模型存在全局解,并且解指数收敛到平衡解(Theorem 4.1),衰减速率 $\omega > 0$ 由线性化算子的谱间隙决定。相比之下,[9] 仅证明了Hibler模型的局部适定性和接近平衡解的全局适定性,但未详细分析长期行为。
  • 本文首次证明了固定冰模型时间周期解的存在性(Theorem 5.2),适用于足够小的时间周期外力(如风力)。这扩展了 [10] 中对Hibler模型的周期解分析,特别处理了 $k_t \equiv 1$ 的情况,此时需要正则化参数满足小条件。
  • 数值模拟显示固定冰模型会形成以消失冰速为特征的静止平衡状态,这与Theorem 4.1的全局存在结果一致,而经典粘塑性海冰模型通常不会形成这种状态。这揭示了固定冰与普通海冰在长期行为上的本质区别。
  1. 假设差异与适用范围:
  • 本文假设初始数据在 $C^1$ 空间中(通过 $X_{\gamma,\mu} \hookrightarrow C^{1,\alpha}(G)^4$ 保证),而 [16] 中的弱解分析允许更低的正则性。但本文的强解提供了更详细的解结构信息,适用于需要高精度模拟的气候模型。
  • 与 [9] 相比,本文需要额外假设 $h > \kappa > 0$(避免冰厚为零的奇点),这是固定冰模型的物理要求,因为基底应力 $f_b$ 在 $h = 0$ 附近行为特殊。
  • 对于 $k_t \equiv 1$ 的情况(即拉伸强度等于冰强度),本文需要正则化参数 $\Delta_{\text{min}}$ 足够小,这限制了结果的适用范围。相比之下,当 $k_t < 1$ 时,这一条件可以放宽。
  1. 代价与局限:
  • 为了进行严格的数学分析,本文引入了人工扩散项 $d_h\Delta h$ 和 $d_A\Delta A$,而实际气候模型中可能不需要这些项。虽然作者在Section 3.3讨论了无扩散项的情况,但局部适定性证明需要更复杂的分析。
  • 与数值模拟相比,理论分析需要更严格的假设,如初始数据的正则性要求和外部力的小性条件。特别是时间周期解的存在性仅对足够小的周期外力成立,而实际气候系统中的外力可能不满足这一条件。
  • 本文主要关注固定冰模型的确定性分析,未考虑随机扰动,而实际气候系统中存在显著的随机性。未来工作可能需要扩展到随机PDE框架。

总体而言,本文首次填补了固定冰模型严格数学分析的空白,将海冰模型的PDE理论从Hibler模型扩展到更复杂的固定冰模型,揭示了固定冰在长期行为上的独特性质,并为气候模型的数学基础提供了新的见解。


了解领域论文(摘要翻译)

1. Gagliardo–Nirenberg不等式的一个推广及其在质量临界与质量次临界椭圆方程中的应用

A generalisation of the Gagliardo–Nirenberg Inequality with applications to mass-critical and mass-subcritical elliptic equations

作者: Bartosz Bieganowski, Jacopo Schino

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24381v1

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

摘要翻译

中文标题: Gagliardo–Nirenberg不等式的一个推广及其在质量临界与质量次临界椭圆方程中的应用

摘要翻译: 通过一个新型的类Gagliardo–Nirenberg不等式,我们在质量临界与质量次临界情形下证明了如下问题解的存在性与不存在性:\begin{equation*} \begin{cases} (-Δ)^s u + \frac{μ}{|y|^{2s}} u + λu = f(u), \quad \mathbb{R}^N \ni x = (y,z) \in \mathbb{R}^K \times \mathbb{R}^{N-K}, \ \int_{\mathbb{R}^N} u^2 , \mathrm{d}x = ρ\end{cases} \end{equation*}其中 $s>0$,$N \ge K \ge 2$,$μ\in \mathbb{R}$ 属于某一特定范围,$ρ>0$ 是预先给定的,而 $λ\in \mathbb{R}$ 是未知拉格朗日乘子。此外,我们还得到了当 $μ=0$ 且 $N \ge 1$ 时上述问题的类似结果,以及一个相关的curl-curl型方程的结果。最后,我们对划分负最小能量与零最小能量两种情形的质量阈值 $ρ$ 进行了系统深入的分析。


2. 含拉普拉斯算子与双拉普拉斯算子的某些积分-微分方程组的定态解存在性

Existence of stationary solutions for some systems of integro-differential equations with Laplace and bi-Laplace operators

作者: Vitali Vougalter, Vitaly Volpert

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24327v1

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

摘要翻译

中文标题: 含拉普拉斯算子与双拉普拉斯算子的某些积分-微分方程组的定态解存在性

摘要翻译: 本文研究一类积分-微分方程组的可解性,其中扩散项由标准拉普拉斯算子与双拉普拉斯算子之差构成。解的存在性证明基于不动点方法;我们运用了无Fredholm性质的椭圆型算子在无界区域上的可解性条件。


3. 非对称快速反应极限中的界面消失现象

Vanishing interfaces in an asymmetric fast reaction limit

作者: Yuki Tsukamoto

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24318v1

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

摘要翻译

中文标题: 非对称快速反应极限中的界面消失现象

摘要翻译: 我们研究一类具有非对称反应项的双组分反应-扩散系统的快速反应极限,其中仅有一个组分具有扩散性。对于非负且互不重叠的初始数据,我们证明初始界面在时刻 t = 0+ 处瞬时消失。更精确地,扩散组分一致收敛于热方程的解,而非扩散组分在初始时刻以外的所有时刻均趋于零。该证明基于显式上下界函数(barriers)及比较原理,并同时适用于 Dirichlet 边界条件与 Neumann 边界条件。


4. 带二维控制输入的Benjamin-Ono方程:近似可控性及其应用

The Benjamin-Ono equation with 2D control input: approximate controllability and its application

作者: Jia-Cheng Zhao

类别: math.OC, math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24130v1

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

摘要翻译

中文标题: 带二维控制输入的Benjamin-Ono方程:近似可控性及其应用

摘要翻译: 我们证明了在环面上,通过二维控制输入实现非线性Benjamin-Ono方程在 $L^2$ 空间中的近似可控性。我们的证明基于Agrachev与Sarychev所提出的几何控制方法,并对其进行了适配。作为该可控性结果的应用,我们研究了一类随机强迫方程的长时间动力学行为。当随机外力非退化且在时间上统计周期时,我们证明其轨道在Sobolev范数下几乎必然无界。


5. 具有各向异性奇异时变扩散系数的热方程的极弱解

Very weak solutions of the heat equation with anisotropically singular time-dependent diffusivity

作者: Zhirayr Avetisyan, Zahra Keyshams, Monire Mikaeili Nia, Michael Ruzhansky

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23856v1

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

摘要翻译

中文标题: 具有各向异性奇异时变扩散系数的热方程的极弱解

摘要翻译: 我们研究具时变、各向异性且可能奇异的扩散系数张量的热方程。由于在此设定下,Sobolev意义下的弱解或分布意义下的解可能并不存在,我们采用极弱解框架,证明了具奇异、各向异性、时变扩散系数的热方程的解的存在性与唯一性。


6. 具Hartree型非线性的半线性阻尼波动方程在齐次Besov空间初值下的一个新临界指数

A new critical exponent for the semilinear damped wave equation with Hartree-type nonlinearity and initial data from homogeneous Besov spaces

作者: Duc An Phan

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23764v1

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

摘要翻译

中文标题: 具Hartree型非线性的半线性阻尼波动方程在齐次Besov空间初值下的一个新临界指数

摘要翻译: 本文研究具Hartree型非线性 $\mathcal{I}_γ\left(|u|^{p_1}\right)|u|^{p_2}$(其中 $p_1, p_2>0$,$γ\in[0, n)$)的半线性阻尼波动方程的临界指数问题,初值取于齐次Besov空间 $\dot{B}_{2, \infty}^{-β}$,其中 $β\in\left[0, \frac{n}{2}\right)$。我们的方法基于对初值属于 $\dot{B}_{2, \infty}^{-β}$ 的相应线性阻尼波动方程解导出衰减估计,并结合调和分析中的精细工具。由此,我们确定了一个新的临界指数

$$ p_1+p_2:=p_{\mathrm{Fuji}}\left(\tfrac{n+2β}{2+γ}\right):=1+\tfrac{4+2γ}{n+2β} \quad \text{ 对 } β\in\left[0, \tfrac{n}{2}\right) \text{ 且 } γ\in [0, n). $$

更精确地,在超临界与临界区域 $p_1+p_2 \geq p_{\mathrm{Fuji}}\left(\frac{n+2 β}{2+γ}\right)$ 中,我们建立了小初值解的整体(时间)存在性;而在次临界区域 $2<p_1+p_2<p_{\mathrm{Fuji}}\left(\frac{n+2 β}{2+γ}\right)$ 中,我们证明了弱解的有限时间爆破现象,即使初值任意小亦然。


7. 黎曼流形上Leibenson方程解的长时间上界估计

Long time upper bounds for solutions of Leibenson’s equation on Riemannian manifolds

作者: Philipp Sürig

类别: math.AP, math.DG

发布日期: 2026-04-25

arXiv ID: 2604.23227v1

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

摘要翻译

中文标题: 黎曼流形上Leibenson方程解的长时间上界估计

摘要翻译: 我们在黎曼流形上考虑Leibenson方程

$$ \partial _{t}u=Δ_{p}u^{q}. $$

我们证明:该方程弱解的某一类长时间上界估计等价于欧氏型Sobolev不等式。


8. 三维Zakharov-Kuznetsov方程的散射性

Scattering of the 3D Zakharov-Kuznetsov equation

作者: Philippe Anjolras

类别: math.AP

发布日期: 2026-04-24

arXiv ID: 2604.22957v1

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

摘要翻译

中文标题: 三维Zakharov-Kuznetsov方程的散射性

摘要翻译: 我们考虑三维空间中的Zakharov-Kuznetsov方程:

$$ \left\{ \begin{array}{l} \partial_t u + \partial_x \Delta u + \partial_x \frac{u^2}{2} = 0 \\ u(t = 0) = u_0 \end{array} \right. $$

其中 $u : (t, x, y) \in \mathbb{R} \times \mathbb{R} \times \mathbb{R}^2 \mapsto u(t, x, y) \in \mathbb{R}$,且 $\Delta = \partial_x^2 + \Delta_y$ 为全拉普拉斯算子。我们证明:对任意初值 $u_0$ 满足

$$ \Vert (1 + x^2 + |y|^2) u_0 \Vert_{H^1} \ll 1, $$

其整体解在 $H^1$ 空间中具有散射性。该结论通过时空共振方法——更精确地说,采用为处理各向异性而提出的偏对称性方法 [GPW23]——得以建立。我们引入适配的各向异性加权范数,证明了适配于这些范数的色散衰减估计,并建立了一致先验估计,从而可通过Bootstrap论证完成闭合。


9. 周期性动力学 Fokker-Planck 方程的大尺度正则性

Large-Scale Regularity for the Periodic Kinetic Fokker-Planck equation

作者: Philip Gaddy

类别: math.AP

发布日期: 2026-04-24

arXiv ID: 2604.22956v1

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

摘要翻译

中文标题: 周期性动力学 Fokker-Planck 方程的大尺度正则性

摘要翻译: 我们首先证明线性动力学 Fokker-Planck 方程基本解的均匀化结果:该基本解在平均 $L^2$ 意义下收敛于一个有效热方程的基本解,其有效扩散系数由定义在环面 $\mathbb{T}^d$ 上的相伴胞腔问题的校正函数所确定。证明的一个关键特征是必须引入二阶校正函数以控制速度变量的平均过程,并需处理由解的空间正则性受限所导致的非散度型误差项。

此外,在上述均匀化结果基础上,我们建立了该 Fokker-Planck 方程解的大尺度正则性结果:更具体地,我们证明解可在大尺度区域上被一类异质多项式(类比于 Taylor 多项式)一致逼近,且具有显式的误差估计;进一步,我们在更大的尺度范围内证明该逼近多项式本身精确满足该 Fokker-Planck 方程。


10. 线性动理学Fokker-Planck方程的Nash-Aronson型估计

Nash-Aronson Estimate for the Linear Kinetic Fokker-Planck equation

作者: Philip Gaddy

类别: math.AP

发布日期: 2026-04-24

arXiv ID: 2604.22944v1

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

摘要翻译

中文标题: 线性动理学Fokker-Planck方程的Nash-Aronson型估计

摘要翻译: 我们证明了带摩擦项的线性动理学Fokker-Planck方程基本解的Nash-Aronson型上界估计,并区分两种时间尺度情形:在长时间尺度下,我们导出一个与经典抛物型估计相匹配的高斯型上界,该估计反映了此情形下速度变量所发生的平均化效应;在短时间尺度下,基本解由常系数Kolmogorov方程的基本解主导,此时摩擦项与势能项可忽略不计。


不相关论文(基本信息)

1. 任意维的双曲面型与时空正质量定理

The Hyperboloidal and Spacetime Positive Mass Theorem in All Dimensions

作者: Sven Hirsch, Marcus Khuri, Martin Lesourd, Yiyue Zhang

类别: math.DG, gr-qc, math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24746v1

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

摘要: Using the recent work of Brendle–Wang on the Riemannian positive mass theorem, we prove the spacetime positive mass theorem for asymptotically flat and asymptotically hyperboloidal initial data sets in arbitrary dimension $n$.


2. 具有给定质量的耦合薛定谔系统多重正解:质量混合与索博列夫临界耦合情形

Multiple positive solutions with prescribed masses for a coupled Schrödinger system: mass mixed and Sobolev critical coupled case

作者: Qing Guo, Qihan He, Wei Shuai, Xuexiu Zhong

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24438v1

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

摘要: The aim of this paper is to establish multiple positive normalized solutions $(u,v,λ_1,λ_2)\in H^1(\mathbb{R}^N,\mathbb{R}^2)\times \mathbb{R}^2$ to the following coupled Schrödinger system involving Sobolev critical exponent:

$$ \begin{cases} -Δu+λ_1 u=μ_1|u|^{p-2}u+να|u|^{α-2}u|v|^β, x\in \mathbb{R}^N,\\ -Δv+λ_2 v=μ_2|v|^{q-2}v+νβ|v|^{β-2}v|u|^α, x\in \mathbb{R}^N,\\ \int_{\mathbb{R}^N}|u|^2\mathrm{d}x=a, \int_{\mathbb{R}^N}|v|^2\mathrm{d}x=b, \end{cases} N\geq 3, $$

where $μ_1,μ_2, ν, a, b>0$. We are particularly interested in the mass mixed case that $2<p, q<2+\frac{4}{N}, α>1, β>1$, and $α+β=2^*:=\frac{2N}{N-2}$. For sufficiently small $ν>0$, we demonstrate that the above system admits two positive solutions, one of which serves as a local minimizer, and the other as a mountain pass solution. By developing some new technical lemmas on the interaction estimates, we are managed to resolves Soave’s open problem [{\it J. Funct. Anal.}, 2020, Remark 1.1] within the context of the system case. Notably, our existence result holds true for all dimensions $N\geq 3$. Our results also significantly extend the result of Gou and Jeanjean [{\it Nonlinearity}, 2018, Theorem 1.1] to the Sobolev critical coupled case and removing the hypothesis ``either $p,q\leq α+β-\frac{2}{N}$ or $|p-q|\leq \frac{2}{N}$" for $N\geq 5$. Additionally, we also establish a sequence of properties for the local minimizer, including local uniqueness, continuity with respect to the small parameter $ν$, and the limiting profiles for $ν\rightarrow 0^+$.


3. 变截面几何下非线性广义热方程的李对称性分析

Lie symmetry analysis of the nonlinear generalized heat equation for varying cross-section geometry

作者: Targyn A. Nauryz

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24418v1

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

摘要: We study the nonlinear generalized heat equation $C(u)u_t=\frac{1}{z^ν}\left(K(u)z^νu_z\right)_z$, where $C(u)$ and $K(u)$ are temperature-dependent thermal coefficients and $ν>0$ is a geometric parameter describing the varying cross-section geometry. By applying the classical Lie symmetry method, we derive the determining equations and perform a complete classification of the admitted Lie point symmetries according to the functional dependence between $C(u)$ and $K(u)$. The analysis shows that the symmetry structure splits naturally into two principal cases: $C(u)/K(u)$ non-constant and $C(u)/K(u)=β$ constant. In the first case, only the basic symmetries are admitted for arbitrary coefficients, whereas additional generators appear under special compatibility relations. In the second case, the equation can be transformed to a linear radial heat equation by the substitution $v=\int K(u)du$, yielding an extended symmetry algebra. For each case, we construct the infinitesimal generators, commutator tables, one-parameter transformation groups, and corresponding invariant reductions. Invariant and similarity solutions are obtained and then specialized to several physically relevant subclasses, including power-law, exponential-type, and linear constitutive coefficients. The results provide a unified symmetry-based model for the analysis of generalized nonlinear heat equations in non-Cartesian geometries.


4. 从标准区域上的偏微分方程到分形上的自相似粒子系统

From PDEs on standard domains to self-similar particle systems on fractals

作者: Georgi Medvedev, Emmanuel Trélat

类别: math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24388v1

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

摘要: We construct transported PDEs on self-similar fractal domains from reference equations posed on the unit interval, and derive explicit self-similar interacting particle systems that approximate the resulting dynamics. The construction combines a measure-preserving isometry between $L^2$-spaces on $[0,1]$ and on the fractal \cite{Med2026}, a nonlocal-to-local approximation of differential operators \cite{PauTre2025}, and a Galerkin discretization on the canonical self-similar partitions. This yields a two-parameter approximation scheme whose error separates a nonlocal consistency term from a Galerkin network term. We work out the transport, Burgers, and heat equations, discuss the relation with intrinsic operators on fractals, and outline extensions to local charts and to pullbacks of nonlocal equations on fractal domains. Moreover, the reverse mapping transforms a nonlocal evolution equation on a fractal domains into the evolution equation on the unit interval, where the methods of classical numerical analysis can be applied. This suggests a promising direction for the development of numerical methods for nonlocal models on fractals, including fractional heat equation, fractal scattering, and related models.


5. 具有一般余维奇点的各向异性磁算子的局部正则性

Local regularity for anisotropic magnetic operators with general codimension singularities

作者: Giovanni Siclari, Stefano Vita

类别: math.AP, math-ph

发布日期: 2026-04-27

arXiv ID: 2604.24314v1

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

摘要: We study local regularity properties of solutions to stationary anisotropic magnetic Schrödinger equations in $\mathbb{R}^d$, $d \ge 2$, arising from singular magnetic potentials concentrated along manifolds of general codimension $2 \le n \le d$. The magnetic interaction is modeled through a covariant gradient of the form

$$ \nabla_m u = (iM\nabla + A)u, $$

where $M^T M$ is a uniformly elliptic matrix encoding anisotropy and $A$ is a magnetic potential with critical Hardy-type scaling along the $n$-codimensional singular set $Σ_0$; that is, $A\sim \mathrm{dist}(\cdot,Σ_0)^{-1}$. We establish local Hölder $C^{0,α}$ and Schauder $C^{1,α}$ estimates for weak solutions via a blow-up analysis adapted to the magnetic structure. The regularity is deeply influenced by the combined effect of anisotropy and the singular magnetic potential, which determines the spectrum of the limiting spherical Laplace-Beltrami operator arising in the blow-up at the singular set. Our model is motivated by the study of magnetic potentials generated by shrinking solenoids onto an axis $Σ_0$, in the three-dimensional setting $d=3$, $n=2$, leading to Aharonov-Bohm-type (AB) models. In this framework, we show that the geometry of the solenoidal loops plays a crucial role: in particular, any deviation from planar cross-sections orthogonal to $Σ_0$ induces a twofold effect. On the one hand, it breaks the ideal AB configuration, in the sense that the magnetic field outside the solenoid is no longer vanishing. On the other hand, it yields an unexpected regularizing mechanism on the wave functions, through a positive shift in the eigenvalues of the asymptotic spectral problem. This purely three-dimensional effect is consistent with our $C^{1,α}$ regularity.


6. Aubry集与Mather集的半离散逼近

Semi-Discrete Approximation of Aubry and Mather sets

作者: Fabio Camilli, Cristian Mendico

类别: math.DS, math.AP, math.OC

发布日期: 2026-04-27

arXiv ID: 2604.24148v1

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

摘要: We study the semi-discrete approximation of Aubry and Mather sets for Tonelli Lagrangians on the flat torus. Starting from the discrete Lax–Oleinik equation, we introduce natural discrete analogues of these sets and analyze their convergence, as the time step tends to zero, in the sense of Kuratowski. Our results show that the semi-discrete variational framework captures not only the ergodic constant, but also the minimizing invariant geometry of the continuous dynamics. In full generality, we prove upper Kuratowski limit inclusions for both the Aubry and Mather sets. For the Aubry set, we establish full convergence under a hyperbolicity assumption on the continuous Aubry set. For the Mather set, we prove full convergence under a genericity assumption ensuring that the Lagrangian admits finitely many ergodic Mather measures. This provides a first rigorous step toward a structure-preserving approximation theory for Aubry and Mather sets in the Tonelli setting, and clarifies how discrete variational models recover the central geometric objects of weak KAM and Aubry–Mather theory.


7. Caginalp型四阶非等温肿瘤生长模型的最优控制

Optimal control for a fourth-order nonisothermal tumor growth model of Caginalp type

作者: Cavalleri Giulia, Pierluigi Colli, Elisabetta Rocca

类别: math.OC, math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24108v1

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

摘要: We study a distributed optimal control problem for a nonisothermal Caginalp-type phase-field model that describes tumour growth under thermal therapy. The PDE system couples a possibly viscous Cahn-Hilliard equation, governing the evolution of the healthy and tumor phases, with an equation for the heat balance, and a reaction-diffusion equation for the nutrient concentration. Chemotaxis and active transport effects are taken into account, and hyperthermia appears as a control variable. We introduce a suitable tracking-type cost functional and show the existence of optimal controls. Then, we analyse the differentiability of the control-to-state operator and establish necessary first-order conditions expressed through a variational inequality involving the adjoint state variables.


8. QED:一个面向开放数学问题证明生成的开源多智能体系统

QED: An Open-Source Multi-Agent System for Generating Mathematical Proofs on Open Problems

作者: Chenyang An, Qihao Ye, Minghao Pan, Jiayaun Zhang

类别: cs.AI, math.AP

发布日期: 2026-04-27

arXiv ID: 2604.24021v1

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

摘要: We explore a central question in AI for mathematics: can AI systems produce original, nontrivial proofs for open research problems? Despite strong benchmark performance, producing genuinely novel proofs remains an outstanding challenge for LLMs. Through systematic experiments with frontier LLMs on research-level proof tasks, we identify seven failure modes that prevent reliable proof generation, including context contamination, citation hallucination, hand-waving on key steps and misallocation of proof effort, unstable proof plans, unfocused verification, problem modification and single-model bottleneck. We argue that the gap between benchmark success and research-level proving is primarily one of system design, due to those failure modes. We present QED, an open-source multi-agent proof system in which each architectural decision directly addresses a specific failure mode. Evaluated on five open problems in applied analysis and PDEs contributed by domain experts, QED produces correct proofs for three problems, each verified by the contributing experts as original and nontrivial. QED is released as open-source software at https://github.com/proofQED/QED.


9. 关于具有常数量曲率的收缩梯度Ricci孤立子刚性的注记

A note on Rigidity of Shrinking Gradient Ricci Solitons with Constant Scalar Curvature

作者: Chen Wang, Guoqiang Wu

类别: math.DG, math.AP

发布日期: 2026-04-27

arXiv ID: 2604.23939v1

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

摘要: Let $(M^n, g, f)$ be an $n$-dimensional complete noncompact gradient shrinking Ricci soliton with the equation $Ric+\nabla^2f= \frac{1}{2}g$.

  1. If its scalar curvature is $\frac{k}{2}$, Ricci curvature is nonnegative and sectional curvature has upper bound $\frac{1}{2(k-1)}$, we prove that the Ricci shrinker is isometric to a finite quotient of $\mathbb{R}^{n-k}\times \mathbb{S}^k$.
  2. If $M$ has constant scalar curvature $R=\frac{n-2}{2}$, and each level set of $f$ has vanishing Weyl curvature, we prove that it is a finite quotient of $\mathbb{R}^2\times \mathbb{S}^{n-2}$. This can be seen a generalization of Cheng-Zhou’s four dimensional result \cite{Cheng-Zhou} to high dimension, since the level set of the potential function $f$ has vanishing Weyl curvature automatically when $n=4$.

10. Allen–Cahn方程及系统中拓扑方法综述

Survey on topological methods for Allen–Cahn equations and systems

作者: João Henrique Andrade, Stefano Nardulli, Raoní Ponciano

类别: math.AP, math.DG

发布日期: 2026-04-27

arXiv ID: 2604.23920v1

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

摘要: We present a survey on multiplicity results for the Allen–Cahn equation and systems in the singular perturbation regime, emphasizing their geometric interpretation through $Γ$-convergence and isoperimetric theory. In the scalar case, the Allen–Cahn functional converges to perimeter, giving rise to minimal and constant-mean-curvature hypersurfaces, while vectorial Allen–Cahn systems lead to multi-phase isoperimetric clusters. The main methodological tool discussed is the photography method, a variational-topological approach based on localized approximate solutions and barycenter maps, which enables one to encode the topology of the ambient manifold into multiplicity results. We compare problems posed on closed manifolds with those on manifolds with boundary, describing the distinct geometric effects induced by Neumann and Dirichlet boundary conditions. The survey highlights both the effectiveness and the limitations of this framework, particularly in the vectorial case, where the lack of a full classification of isoperimetric clusters creates fundamental analytical challenges.


11. 二维两相Bernoulli自由边界的精细结构

Fine structure of the two-phase Bernoulli free boundaries in 2D

作者: Lorenzo Ferreri, Luca Spolaor, Bozhidar Velichkov

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23843v1

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

摘要: We prove that the branching set of a solution to a two-dimensional two-phase Bernoulli problem with constant coefficients is locally finite. We do this via a Weierstrass representation formula, which allows to transform the problem into a new geometric two-phase problem for capillary minimal surfaces. We also apply this method to the obstacle problem establishing a connection between the directional derivatives of solutions to the obstacle problem and the linear thin two-membrane problem.


12. 超导自由边界问题爆破解的唯一性

Uniqueness of Blow-ups for the Superconductivity Free Boundary Problem

作者: Shibing Chen, Yuanyuan Li, Xianduo Wang

类别: math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23682v1

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

摘要: We study the free-boundary equation

$$ Δu=χ_{\{|\nabla u|>0\}} $$

near the origin. We prove that, at a singular point of $\partial\{|\nabla u|>0\}$, the quadratic blow-up is unique. As noted in \cite[Notes to Chapter 7]{PSU2012}, little is known about the singular set for this problem. The usual Weiss–Monneau monotonicity argument does not seem to apply directly, because the inactive set is determined by the vanishing of $\nabla u$, rather than by a sign condition on $u$. The proof follows the quadratic part of the rescalings. Projecting onto the trace-free quadratic harmonics yields a finite-dimensional differential equation for the quadratic coefficient. Together with a Lyapunov identity and estimates on dyadic annuli, this implies convergence of the quadratic coefficient, and hence uniqueness of the blow-up.


13. 关于黎曼流形上微分不等式 $Δ_pu+ Δ_q u+V(x)u^s\leq 0$ 的非负解

On nonnegative solutions of the differential inequality $Δ_pu+ Δ_q u+V(x)u^s\leq 0$ on Riemannian manifolds

作者: Biqiang Zhao

类别: math.AP, math.DG

发布日期: 2026-04-26

arXiv ID: 2604.23624v1

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

摘要: In this paper, we are concerned with differential inequalities with $(p,q)$-Laplacian operator on Riemannian manifolds. Using a test function argument, we establish Liouville-type theorems under the manifold’s geometry and the potential’s behavior at infinity.


14. 单位分解高斯Kolmogorov-Arnold网络

Partition-of-Unity Gaussian Kolmogorov-Arnold Networks

作者: Amir Nooeizadegan

类别: cs.CE, cs.AI, math.AP

发布日期: 2026-04-26

arXiv ID: 2604.23599v1

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

摘要: Gaussian basis functions provide an efficient and flexible alternative to spline activations in KANs. In this work, we introduce the partition-of-unity Gaussian KAN (PU-GKAN), a Shepard-type normalized Gaussian KAN in which the Gaussian basis values on each edge are divided by their local sum over fixed centers. This produces a partition-of-unity feature map with trainable coefficients, while preserving the standard edge-based KAN structure. The normalized construction gives exact constant reproduction at the edge level and admits an explicit finite-feature kernel interpretation. We formulate both the standard Gaussian KAN (GKAN) and PU-GKAN from a finite-feature and additive-kernel viewpoint, making the induced layer kernels and empirical feature matrices explicit. Using the first-layer feature matrix as the reference object, we adopt a practical scale-selection interval for $ε$, with the lower endpoint determined by adjacent-center overlap and the upper endpoint determined by a conservative conditioning threshold. Numerical experiments show that PU-GKAN reduces sensitivity to $ε$, improves validation accuracy for most smooth and moderately non-smooth targets, and gives more stable training behavior. The benefit persists across sample-size and center-number sweeps, higher-dimensional architectures, Matérn RBF bases, and physics-informed examples involving Helmholtz and wave equations. These results indicate that Shepard-type partition-of-unity normalization is a simple and effective stabilization mechanism for RBF-based KANs.


15. 具混合奇异-退化结构的拟线性抛物方程粘性解的梯度正则性

Gradient regularity for viscosity solutions to quasilinear parabolic equations with mixed singular-degenerate structure

作者: Junior da Silva Bessa, João Vitor da Silva, Ginaldo de Santana Sá

类别: math.AP

发布日期: 2026-04-25

arXiv ID: 2604.23421v1

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

摘要: We establish regularity results for viscosity solutions to a class of quasilinear parabolic equations exhibiting nonhomogeneous degeneracy or singularity (a double phase regime) of the form

$$ u_t - \big(|Du|^{\mathfrak{p}} + \mathfrak{a}(x,t)|Du|^{\mathfrak{q}}\big)Δ_p^{\mathrm N} u = f(x,t) \quad \text{in } Q_1, $$

where $-1 < \mathfrak{p} < 0$, $\mathfrak{p} \leq \mathfrak{q}$, and $\mathfrak{a}, f : Q_1 \to \mathbb{R}$ are prescribed functions. Using the Jensen–Ishii method, we prove Lipschitz regularity for appropriately translated solutions. Moreover, combining this approach with intrinsic scaling techniques, we establish interior Hölder continuity estimates for the gradient. Our results extend recent work of Fang and Zhang on the homogeneous case via a different approach.


16. 任意维数下一类Hessian商方程的半凸解的内 $C^{2}$ 估计

Interior $C^{2}$ estimate for semi-convex solutions to a class of Hessian quotient equations in arbitrary dimensions

作者: Xinqun Mei, Jin Yan

类别: math.AP

发布日期: 2026-04-25

arXiv ID: 2604.23349v1

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

摘要: In this paper, we study the interior $C^{2}$ estimates for Hessian quotient equations $\frac{σ_{3}(D^{2}u)}{σ_{l}(D^{2}u)}=1$ for $l=1, 2$, in arbitrary dimensions, under the natural ellipticity and semi-convexity conditions. We further derive analogous results for the corresponding sum Hessian equations. In addition, we establish several rigidity results.


17. 略低于临界情形的狄利克雷问题中具有剩余质量的爆破解

Blowing-up Solutions with Residual Mass in a Slightly Subcritical Dirichlet Problem

作者: Rufaidah Alharbi, Mohamed Ben Ayed, Khalil El Mehdi

类别: math.AP

发布日期: 2026-04-25

arXiv ID: 2604.23339v1

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

摘要: In this paper, we study the Dirichlet elliptic problem $(\mathcal{P}_\varepsilon)$: $-Δu +V\,u = u^{p-\varepsilon}$, $u>0$ in $Ω$, $u=0$ on $\partialΩ$, where $Ω\subset \R^n$ ( $n\geq 3$) is a bounded domain, $V$ is a smooth positive function on $\overlineΩ$, $p+1= 2n/(n-2)$ is the critical Sobolev exponent, and $\varepsilon >0$ is a small parameter. First, we show that, unlike the case of weak convergence to zero, interior bubbling solutions with a nonzero weak limit cannot occur in low dimensions. We then treat the general setting by removing the restriction that blow-up points are confined to the interior. Using delicate asymptotic expansions of the gradient of the associated functional, we prove that in dimensions $n=4$ and $n=5$, single blow-up point cannot coexist with residual mass.\ We further elucidate the role of the sign of the normal derivative of the potential $V$ on the boundary: if it is positive, any single blow-up solution with residual mass must occur in the interior; if it is negative at some boundary point, boundary blow-up solutions with residual mass can be constructed. Finally, we construct both simple and non-simple interior blow-up solutions exhibiting residual mass, without any assumption on the sign of the normal derivative of $V$. These results provide new insights into the interaction between the potential, the geometry of the domain, and the critical nonlinearity.


18. 一种拓扑固定型形状梯度框架:用于带排斥能的CIE Lab彩色图像的非简单边界提取

A Topology fixated Shape Gradient Framework for Non Simple Boundary Extraction for CIE Lab color images with Repulsive Energy

作者: Shafeequdheen Palengara, Jyotiranjan Nayak, Vijayakrishna Rowthu

类别: cs.CV, math.AP

发布日期: 2026-04-25

arXiv ID: 2604.23167v1

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

摘要: A levelset free but a hybrid image segmentation approach based on a modified version of the piece wise constant shape gradient of an Mumford Shah shape functional and a repulsive function is considered. The segmentation is performed a non-local shape based through an evolution of discrete curves driven by a non local shape based energy to segment images containing disjoint regions and multiple boundaries. This formulation has a novel additional component as a multivariable function dependent on a few sampled points of the curves that handles the occurrence of self intersection during boundary curves evolution. The method is applied to a few gray scale and color images, including images with nested structures and astronomical objects. The results indicate effective segmentation in complex scenarios with absolute control on the topology of the segments and self-intersections of the boundaries