生成时间: 2026年04月14日 05:17:59

重点关注论文数量: 7

了解领域论文数量: 28

不相关论文数量: 14

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

1. 含被动示踪剂的不可压缩流的瞬时爆破

Instantaneous blowup of incompressible flow with passive tracer

作者: Mimi Dai, Xiaotong “Dawson” Yang

类别: math.AP

发布日期: 2026-04-13

arXiv ID: 2604.11769v1

摘要翻译: 我们构造了一族含被动示踪剂的不可压缩流解 $$(u,b)$$ ，使得 $

Comment: 无

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

1. 研究对象和背景

本文研究如下不可压缩流与被动示踪剂耦合系统的爆破行为：

\begin{cases} \partial_t u - \Delta u + \mathrm{div}(u \otimes u) + \nabla p = 0, \\ \partial_t b - \Delta b + (u \cdot \nabla) b = 0, \\ \mathrm{div}\, u = 0, \end{cases} \quad \text{在 } \mathbb{T}^d \times [0,T] \text{ 或 } \mathbb{R}^d \times [0,T],\ d \geq 2.

其中 $$u(x,t)$$ 是速度场， $$p(x,t)$$ 是压力， $$b(x,t)$$ 是标量示踪剂。该系统可视为 Navier–Stokes 方程与线性对流-扩散方程的耦合：一旦 $$u$$ 给定， $$b$$ 满足一个由 $$u$$ 驱动的抛物型方程。值得注意的是，在二维情形下，若将磁流体力学（MHD）系统限制为“2.5维”结构——即速度场仅依赖于前两个空间变量且无第三分量，而磁场仅有第三分量 $$B = (0,0,b(x_1,x_2))$$ ——则 MHD 系统恰好退化为上述系统（1.1）。因此，该模型不仅具有流体力学背景，也与 MHD 的简化情形密切相关。

从分析角度看，系统（1.1）继承了 Navier–Stokes 方程的尺度不变性。具体而言，若 $$(u,p,b)(x,t)$$ 是解，则对任意 $\lambda > 0$ 和任意实数 $$a$$ ，缩放后的场

u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t),\quad p_\lambda(x,t) = \lambda^2 p(\lambda x, \lambda^2 t),\quad b_\lambda(x,t) = \lambda^a b(\lambda x, \lambda^2 t)

仍是解。这表明速度场具有标准 Navier–Stokes 尺度，而示踪剂 $$b$$ 因其线性演化性质允许额外的齐次参数 $$a$$ 。当 $$a=1$$ 时，该缩放与 MHD 系统一致。由此导出的速度场的临界空间包括 $\dot{H}^{d/2 -1}$ 、 $$L^d$$ 、 $L^{d,\infty}$ 、 $\dot{W}^{-1,\infty}$ 、 $\mathrm{BMO}^{-1}$ 和 $\dot{B}^{-1}_{\infty,\infty}$ 等。

尽管在二维情形下 Navier–Stokes 方程存在整体光滑解，但本文关注的是弱解框架下的非唯一性与爆破现象。特别地，作者旨在构造一类弱解，其在某个预定时刻 $$T_*$$ 处发生 瞬时爆破（instantaneous blowup）：即在 $$[0,T_*]$$ 上解光滑，在 $$(T_*,T]$$ 上也光滑，但在 $t \to T_*^+$ 时 $\|u(t)\|_{L^\infty}$ 与 $\|b(t)\|_{L^\infty}$ 均趋于无穷，且具有明确的爆破速率。这种现象不同于经典有限时间爆破（因解在爆破点两侧均光滑），而是通过凸积分（convex integration）方法构造的“人工”爆破解，揭示了弱解理论中的病态行为。

此前，Cheskidov–Dai–Pocovnicu（[CDP25]）已在纯 Navier–Stokes 情形中利用 逆级联机制（inverse cascade）构造了类似瞬时爆破解：高频能量被迭代注入低频，导致 $L^\infty$ 范数在极限时刻发散。然而，在耦合系统（1.1）中，新增的示踪剂 $$b$$ 引入了结构性约束：每一步迭代不仅要生成新的低频速度场，还必须同步传播示踪剂，同时保持主部速度剖面的一致性。这一兼容性要求无法通过标准几何引理满足，因为速度方程涉及对称张量（来自 $u \otimes u$ ），而示踪剂方程涉及向量场（来自 $$u b$$ ）。因此，现有凸积分框架缺乏处理此类 张量-向量耦合结构 的工具。本文试图填补这一空白，通过引入一个新的 同时分解引理（simultaneous decomposition lemma），统一处理对称张量与向量场的振荡构造，从而实现耦合系统的逆级联。

2. 主要定理或主要结果

作者证明了以下两个主要定理，分别针对二维和三维情形：

定理 1.2（二维情形）：设 $u_0 \in C^\infty(\mathbb{T}^2)$ 无散， $b_0 \in C^\infty(\mathbb{T}^2)$ 。则存在 $$T > 0$$ ，使得对任意 $T_* \in [0,T)$ ，存在系统（1.1）的弱解 $$(u(t), b(t))$$ 在 $\mathbb{T}^2 \times [0,T]$ 上满足：

在 $$[0,T_*]$$ 上， $$(u,b)$$ 是以 $$(u_0,b_0)$$ 为初值的经典解；
在 $$(T_*,T]$$ 上， $$(u,b)$$ 仍是经典解，且存在常数 $$c > 0$$ ，使得沿某序列 $t_n \to T_*^+$ ，有 $\|u(t_n)\|_{L^\infty} \geq \frac{c}{\sqrt{t_n - T_*}},\quad \|\nabla u(t_n)\|_{L^\infty} \geq \frac{c}{t_n - T_*},$ $\|b(t_n)\|_{L^\infty} \geq \frac{c}{\sqrt{t_n - T_*}},\quad \|\nabla b(t_n)\|_{L^\infty} \geq \frac{c}{t_n - T_*};$ 同时，存在 $$C > 0$$ ，使得对所有 $t \in (T_*,T]$ ， $\|u(t)\|_{L^\infty} + \|b(t)\|_{L^\infty} \leq \frac{C}{\sqrt{t - T_*}},\quad \|\nabla u(t)\|_{L^\infty} + \|\nabla b(t)\|_{L^\infty} \leq \frac{C}{t - T_*};$
轨迹 $$(u(t), b(t))$$ 在时间上关于 $\mathrm{BMO}^{-1} \times \mathrm{BMO}^{-1}$ 弱-*连续；
解满足 $u, b \in L^2([0,T]; L^p)$ 对任意 $p < \infty$ 。

定理 1.3（三维及以上）：对任意 $u_0 \in C^\infty(\mathbb{T}^3)$ （无散）和 $b_0 \in C^\infty(\mathbb{T}^3)$ ，存在 $$T > 0$$ 及弱解 $$(u,b)$$ 在 $\mathbb{T}^3 \times [0,T]$ 上，使得对任意 $T_* \in [0,T)$ ，该解满足定理 1.2 的全部性质，并额外满足 $u, b \in L^\infty_t \dot{W}^{-1,\infty}_x$ 。

关键假设包括：初值光滑、周期边界条件、以及构造中使用的频率尺度满足严格递增关系 $N_{j,k} \ll N_{j',k'}$ 当 $$k < k'$$ 或 $$k=k', j<j'$$ 。爆破速率精确为 $(t - T_*)^{-1/2}$ （函数值）和 $(t - T_*)^{-1}$ （梯度），与热核衰减尺度一致。此外，解属于 Ladyzhenskaya–Prodi–Serrin 类之外的空间，但仍在某些超临界空间中具有良好性质（如 $$L^2_t L^p_x$$ ）。这些结果表明，即使在二维（Navier–Stokes 整体正则性成立的情形），耦合系统仍可通过弱解表现出强奇异性。

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

本文的核心技术障碍在于：在凸积分迭代过程中，必须同时满足速度方程和示踪剂方程的非线性结构兼容性。具体而言，标准凸积分通过振荡剖面近似非线性项 $\mathrm{div}(u \otimes u)$ ，其主部由对称张量 $$R$$ 表示，通常借助 几何引理（Geometric Lemma）将其分解为若干方向上的秩一矩阵之和（如 $\sum \Gamma_\eta^2 \eta^\perp \otimes \eta^\perp$ ）。然而，在系统（1.1）中，示踪剂方程的非线性项为 $\mathrm{div}(u b)$ ，其主部对应一个向量场 $$g$$ ，需同时用相同的振荡方向表示为 $\sum \Gamma_\eta \eta^\perp$ 。这就要求存在一组系数 $\{\Gamma_m\}$ ，既能重构给定对称张量 $$R$$ ，又能重构给定向量 $$g$$ ，且两者共享相同的基底方向。这是标准几何引理无法处理的新约束。

作者的真正创新点是提出了 张量-向量同时分解引理（Lemma 3.2）。该引理断言：存在 $\varepsilon_* > 0$ 和有限单位向量集 $\Lambda_b \subset \mathbb{S}^1 \cap \mathbb{Q}^2$ ，以及光滑映射 $\Gamma_m: B_{\varepsilon_*}(\mathrm{Id}) \times \mathbb{R}^2 \to (0,\infty)$ 和 $p: B_{\varepsilon_*}(\mathrm{Id}) \times \mathbb{R}^2 \to \mathbb{R}$ ，使得对任意 $(R,g) \in B_{\varepsilon_*}(\mathrm{Id}) \times \mathbb{R}^2$ ，有

R - p(R,g) \mathrm{Id} = \sum_{m=1}^N \Gamma_m^2(R,g) \eta_{m,\perp} \otimes \eta_{m,\perp}, \quad g = \sum_{m=1}^N \Gamma_m(R,g) \eta_{m,\perp}.

证明的关键在于将向量 $$g$$ 的贡献“编码”进张量分解的常数项中：通过引入辅助量 $$M(g) = 1 + |g|^2$$ ，并巧妙构造 $\pm e_1, \pm e_2$ 方向上的系数 $\alpha_i(g), \beta_i(g)$ ，使得这些方向既贡献所需的向量 $$g$$ ，又产生一个各向同性的张量 $M(g) \mathrm{Id}$ 。再结合标准几何引理对 $$R$$ 的分解（使用另一组方向 $\Lambda_u$ ），最终通过设置 $$p(R,g) = -M(g)$$ ，实现张量与向量的同时表示。该引理是整个构造可行性的基石。

在分析框架上，作者沿用 [CDP25] 的 逆级联凸积分方案，但进行了适配。他们将解分解为主部 $$(v,h)$$ 与修正项 $(w,\zeta)$ ，其中主部由多尺度振荡剖面叠加而成：

v_k(x,t) \approx -\sum_{j \in \Lambda_u} N_{j,k} e^{-N_{j,k}^2 t} \mathcal{P} \Delta \psi_{j,k,u}(x) - \sum_{j \in \Lambda_b} N_{j,k} e^{-N_{j,k}^2 t} \mathcal{P} \Delta \psi_{j,k,c}(x),

h_k(x,t) \approx -\sum_{j \in \Lambda_b} N_{j,k} e^{-N_{j,k}^2 t} \Delta \psi_{j,k,b}(x),

其中 $\psi_{j,k,u}, \psi_{j,k,c}, \psi_{j,k,b}$ 是具有高频相位 $\sin(N_{j,k} \eta_j \cdot x)$ 和低频振幅的剖面函数。关键在于设计振幅 $a_{j,k+1,b}$ ，使其同时满足（3.4）中的三个近似等式——这正是同时分解引理的应用场景。通过选择频率尺度 $N_{j,k} = m_* \lceil A^{b k + (j-1)/J} \rceil$ 和时间尺度 $t_k \sim N_{1,k+1}^{-2}$ ，确保在区间 $[t_{k+1}, t_k]$ 上，第 $$k+1$$ 层的高频相互作用主导生成第 $$k$$ 层的低频主部，且交叉项可被控制。

误差估计依赖于热半群与 Leray 投影的振荡估计（如 Lemma 2.1–2.3），这些引理控制了形如 $\nabla^n e^{t\Delta} \mathcal{R} (a(x) \sin(\xi \cdot x))$ 的项，其中 $\mathcal{R}$ 是 Riesz 型算子（如 $\mathcal{R} = \Delta^{-1} \tilde{\mathcal{D}} \mathcal{P} \mathrm{div}$ ）。这些估计保证了修正项 $(w,\zeta)$ 可通过压缩映射原理在适当函数空间中求解，且其范数远小于主部，从而不破坏爆破机制。

4. 与之前工作的比较

本文与已有工作的比较可从多个维度展开：

1. 与纯 Navier–Stokes 爆破构造的比较：本文直接推广了 Cheskidov–Dai–Pocovnicu ([CDP25]) 的瞬时爆破结果。在 [CDP25] 中，作者通过逆级联凸积分在二维和三维 Navier–Stokes 方程中构造了 $L^\infty$ 爆破解，其核心是利用几何引理处理对称应力张量。本文保留了该框架的频率层次、时间尺度和误差控制策略，但面临新挑战：示踪剂方程引入了额外的向量场约束。如文中所述，“被动分量带来了新的兼容性约束：迭代必须在传播示踪剂的同时保持相同的主要速度剖面”。这一约束在纯流体情形中不存在，因此 [CDP25] 的方法无法直接套用。本文通过同时分解引理解决了这一问题，这是对凸积分技术的重要扩展。

2. 与被动示踪剂相关工作的比较：以往关于被动示踪剂的研究多集中于正则性传递（如若 $$u$$ 光滑则 $$b$$ 光滑）或耗散异常，而非爆破构造。例如，在二维情形下，由于 $$u$$ 整体光滑， $$b$$ 通常也被预期整体光滑。本文的结果表明，即使在二维，弱解框架下仍可出现 $$b$$ 与 $$u$$ 同步爆破的现象，这揭示了弱解非唯一性对耦合系统的影响。值得注意的是，作者指出该系统可视为 2.5D MHD 的简化，而 MHD 系统的爆破问题长期开放。尽管本文未直接处理完整 MHD，但其构造为理解 MHD 中磁场与速度场的耦合奇异性提供了新视角。

3. 技术假设与适用范围的推进：与 [CDP25] 相比，本文在二维情形下的证明“更为技术性”（见第4页），主要源于同时分解引理的构造复杂性以及耦合误差项的控制。然而，所得结果在爆破速率和函数空间性质上与 [CDP25] 一致：爆破速率为 $(t-T_*)^{-1/2}$ ，解属于 $$L^2_t L^p_x$$ （任意 $p<\infty$ ）和 $\mathrm{BMO}^{-1}$ 。这表明耦合并未削弱爆破强度。此外，本文明确将结果扩展至三维及以上（定理 1.3），并指出高维证明类似于 Dai 的另一工作 [Dai26]，暗示该方法具有维度普适性。

4. 代价与局限：本文构造的解是弱解，不满足能量等式，且爆破是“人为”的（依赖于凸积分的非物理振荡）。与 [CDP25] 一样，这些解不属于任何 Ladyzhenskaya–Prodi–Serrin 临界类（如 $$L^p_t L^q_x$$ 满足 $$2/p + d/q = 1, q > d$$ ），因此不违背已知的正则性准则。主要局限在于：构造依赖于周期边界条件和光滑初值，且未探讨物理初值（如有限能量）下是否存在类似现象。此外，虽然示踪剂 $$b$$ 被动（不反馈速度场），但其爆破速率与 $$u$$ 相同，这可能暗示在更一般的耦合系统中奇异性会相互增强。

综上，本文在 [CDP25] 的基础上，通过引入新的分解工具，成功将瞬时爆破构造推广至含被动示踪剂的耦合系统，揭示了弱解理论中耦合非线性带来的新结构挑战。

2. 具有不同粘性的质量守恒Navier-Stokes/Allen-Cahn系统的锐界面极限

Sharp Interface Limit for a Mass-Conserving Navier-Stokes/Allen-Cahn System with Different Viscosities

作者: Helmut Abels, Hanifah Mumtaz

类别: math.AP

发布日期: 2026-04-13

arXiv ID: 2604.11597v1

摘要翻译: 我们在二维有界光滑区域中，对耦合的Navier-Stokes方程与质量守恒Allen-Cahn系统在参数 $\varepsilon > 0$ （表示弥散界面厚度）趋于零时的锐界面极限进行了严格分析。我们证明了质量守恒Navier-Stokes/Allen-Cahn系统的解收敛于其锐界面极限问题的解。在此极限下，界面按照带有对流项的质量守恒平均曲率流演化，并与带表面张力的两相Navier-Stokes系统耦合。我们的方法通过匹配渐近展开构造极限系统的近似解，并辅以首阶项的特殊拟设。为估计该近似解与精确解之间的误差，我们利用了在线性化Allen-Cahn算子关于近似解附近的精细谱估计。

Comment: 30 pages

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

1. 研究对象和背景

本文研究二维光滑有界区域 $\Omega \subset \mathbb{R}^2$ 中耦合的、质量守恒的 Navier-Stokes/Allen-Cahn 系统在界面厚度参数 $\varepsilon \to 0$ 下的锐界面极限。该系统由以下方程组描述：

\partial_t v_\varepsilon + v_\varepsilon \cdot \nabla v_\varepsilon - \mathrm{div}(2\nu(c_\varepsilon) D v_\varepsilon) + \nabla p_\varepsilon = -\varepsilon \, \mathrm{div}(\nabla c_\varepsilon \otimes \nabla c_\varepsilon), \tag{1.1}

\mathrm{div}\, v_\varepsilon = 0, \tag{1.2}

\partial_t c_\varepsilon + v_\varepsilon \cdot \nabla c_\varepsilon - \Delta c_\varepsilon = \frac{1}{\varepsilon^2} \left( \fint_\Omega f'(c_\varepsilon)\, dx - f'(c_\varepsilon) \right), \tag{1.3}

配合齐次无滑移与 Neumann 边界条件 $(v_\varepsilon, \partial_n c_\varepsilon)|_{\partial\Omega} = (0, 0)$ ，以及适配的初值 $(v_{0,\varepsilon}, c_{0,\varepsilon})$ 。

其中 $v_\varepsilon$ 为速度场， $p_\varepsilon$ 为压力， $c_\varepsilon$ 为序参量（通常解释为两相浓度差）， $\nu(c_\varepsilon)$ 是依赖于 $c_\varepsilon$ 的粘度函数， $$f$$ 是标准双阱势（如 $f(c) = \frac{1}{8}(c^2 - 1)^2$ ），满足 $f'(\pm 1) = 0$ , $f''(\pm 1) > 0$ ，且对称。关键的是，方程 (1.3) 中的非局部项 $\fint_\Omega f'(c_\varepsilon)\, dx$ 强制系统满足总质量守恒：

\frac{d}{dt} \int_\Omega c_\varepsilon\, dx = 0.

该模型属于“扩散界面”方法，用于模拟不可混溶流体的两相流动。当 $\varepsilon \ll 1$ 时，界面层宽度为 $O(\varepsilon)$ ，而远离界面处 $c_\varepsilon \approx \pm 1$ 对应两个纯相。一个核心问题是：当 $\varepsilon \to 0$ 时，解 $(v_\varepsilon, c_\varepsilon)$ 是否收敛到某个“锐界面”模型的解？后者将界面视为零厚度曲面 $\Gamma_t$ ，并在其两侧分别求解 Navier-Stokes 方程，通过界面条件耦合。

作者证明，极限系统为如下两相 Navier-Stokes 系统耦合带对流项的质量守恒平均曲率流：

\partial_t v^{\pm}_0 + v^{\pm}_0 \cdot \nabla v^{\pm}_0 - \nu^{\pm} \Delta v^{\pm}_0 + \nabla p^{\pm}_0 = 0 \quad \text{in } \Omega^{\pm}(t), \tag{1.7}

\mathrm{div}\, v^{\pm}_0 = 0 \quad \text{in } \Omega^{\pm}(t), \tag{1.8}

界面条件包括无滑移连续性 $$[v_0] = 0$$ 、应力跳跃条件

[2\nu^{\pm} D v^{\pm}_0 - p^{\pm}_0 I] n_{\Gamma_t} = -\sigma H_{\Gamma_t} n_{\Gamma_t}, \tag{1.9}

以及界面运动学方程

V_{\Gamma_t} - n_{\Gamma_t} \cdot v^{\pm}_0 = H_{\Gamma_t} - \fint_{\Gamma_t} H_{\Gamma_t}\, d\mathcal{H}^1. \tag{1.11}

这里 $\nu^{\pm} = \nu(\pm 1)$ ， $\sigma = \int_{\mathbb{R}} (\theta_0')^2\, d\rho$ 是表面张力系数， $\theta_0$ 是 Allen-Cahn 剖面方程 $-\theta_0'' + f'(\theta_0) = 0$ 的异宿轨道。

已有工作如 Abels–Fei [1] 处理了非质量守恒情形；Chen–Hilhorst–Logak [13] 给出了质量守恒 Allen-Cahn 方程的形式渐近分析；Hensel–Liu [21] 和 Kroemer–Laux [28] 分别用相对熵方法处理了常粘度 Navier-Stokes/Allen-Cahn 和质量守恒 Allen-Cahn 的极限。但尚未有严格结果处理质量守恒、且允许不同粘度 $\nu^+ \neq \nu^-$ 的 Navier-Stokes/Allen-Cahn 系统的锐界面极限。本文旨在填补这一空白，并揭示质量守恒项如何影响界面动力学（引入体积保持项 $-\fint_{\Gamma_t} H_{\Gamma_t}$ ）及渐近展开结构。

2. 主要定理或主要结果

作者的主要成果是定理 1.1，它在强解框架下建立了从扩散模型到锐界面模型的高阶收敛性。设 $(v_0^{\pm}, \Gamma)$ 是极限系统 (1.7)–(1.13) 在时间区间 $$[0, T_0]$$ 上的光滑解。则存在一族依赖于 $\varepsilon$ 的近似解 $$(c_A, v_A)$$ ，使得：若扩散系统的初值 $(v_{0,\varepsilon}, c_{0,\varepsilon})$ 满足高阶匹配条件

\|c_{0,\varepsilon} - c_A|_{t=0}\|_{L^2(\Omega)} + \varepsilon^2 \|\nabla(c_{0,\varepsilon} - c_A|_{t=0})\|_{L^2(\Omega)} + \|v_{0,\varepsilon} - v_A|_{t=0}\|_{L^2(\Omega)} \leq C \varepsilon^{N + \frac{1}{2}}, \tag{1.16}

其中 $N \geq 3$ 为任意整数，则存在 $\varepsilon_0 \in (0,1]$ , $T_1 \in (0, T_0]$ , $$R > 0$$ ，使得对所有 $\varepsilon \in (0, \varepsilon_0]$ ，误差估计成立：

序参量 $c_\varepsilon$ 的整体与局部估计：
$\sup_{0 \leq t \leq T_1} \|c_\varepsilon(t) - c_A(t)\|_{L^2(\Omega)} + \|\nabla(c_\varepsilon - c_A)\|_{L^2((\Omega \times (0,T_1)) \setminus \Gamma(\delta))} \leq R \varepsilon^{N + \frac{1}{2}}, \tag{1.17a}$
$\|\nabla_\tau(c_\varepsilon - c_A)\|_{L^2(\Omega \times (0,T_1) \cap \Gamma(2\delta))} + \varepsilon \|\partial_n(c_\varepsilon - c_A)\|_{L^2(\Omega \times (0,T_1) \cap \Gamma(2\delta))} \leq R \varepsilon^{N + \frac{1}{2}}, \tag{1.17b}$
$\|\nabla(c_\varepsilon - c_A)\|_{L^\infty(0,T_1; L^2(\Omega))} + \|\nabla^2(c_\varepsilon - c_A)\|_{L^2(\Omega \times (0,T_1))} \leq R \varepsilon^{N - \frac{3}{2}}. \tag{1.17c}$
速度场 $v_\varepsilon$ 的收敛性：
$\|v_\varepsilon - v_A\|_{L^\infty(0,T_1; L^2(\Omega))} + \|v_\varepsilon - v_A\|_{L^2(0,T_1; H^1(\Omega))} \leq C(R) \varepsilon^{N + \frac{1}{2}}. \tag{1.18}$

此外，近似解具有明确的渐近结构： $c_A \to \pm 1$ 在 $\Omega^{\pm}$ 的紧子集上一致成立，且 $v_A = v_0^{\pm} + O(\varepsilon)$ 在 $L^\infty(\Omega \times (0, T_1))$ 中。特别地，初值近似解 $c_{A,0}$ 的构造为

c_{A,0}(x) = \zeta(d_{\Gamma_0}(x)) \theta_0\left( \frac{d_{\Gamma_0}(x)}{\varepsilon} \right) + (1 - \zeta(d_{\Gamma_0}(x))) (\chi_{\Omega^+(0)} - \chi_{\Omega^-(0)}) + O(\varepsilon),

其中 $\zeta$ 是截断函数， $d_{\Gamma_0}$ 是初始界面 $\Gamma_0$ 的符号距离函数。

这些结果表明，在适当的初值条件下，扩散模型的解在 $\varepsilon \to 0$ 时以任意高阶精度（由 $$N$$ 控制）逼近锐界面模型的解，且误差在界面内外均有精细刻画。

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

证明的核心在于构造高精度的近似解 $$(c_A, v_A)$$ 并控制真实解与近似解之间的误差。作者采用的方法融合了形式渐近展开与严格的能量估计，其创新点主要体现在以下几个方面：

1. 质量守恒约束下的匹配渐近展开：与非守恒情形 [1] 不同，质量守恒项 $\lambda_\varepsilon(t) = \varepsilon^{-1} \fint_\Omega f(c_\varepsilon)\, dx$ 引入了一个全局耦合。这导致在内层展开（使用拉伸变量 $\rho = d_\Gamma / \varepsilon - h_\varepsilon$ ）中， $O(\varepsilon)$ 阶项不再自动消失（见 [1, Remark A.6]）。作者必须引入一个新颖的 ansatz 来处理此阶项，特别是修正剖面 $\theta_1(\rho)$ ，它满足线性化方程

-\theta_1'' + f''(\theta_0) \theta_1 = 1 - \frac{2}{\sigma} \theta_0'. \tag{2.17}

这个修正项对于满足质量守恒至关重要，并保证了后续谱估计的有效性。

2. 关键的谱估计（Theorem 2.6）：为了控制误差 $u = c_\varepsilon - c_A$ ，作者需要估计线性化 Allen-Cahn 算子的能量泛函

\int_\Omega \left( |\nabla u|^2 + \frac{f''(c_A)}{\varepsilon^2} u^2 \right) dx.

由于 $$f''(c_A)$$ 在界面附近变号，该泛函并非正定。作者引用并应用了来自 [33] 的一个精细谱估计：

\int_\Omega \left( |\nabla \psi|^2 + \frac{f''(c_A)}{\varepsilon^2} \psi^2 \right) dx \geq -C_L \int_\Omega \psi^2 dx + \int_{\Omega \setminus \Gamma_t(\delta)} |\nabla \psi|^2 dx + \int_{\Gamma_t(\delta)} |\nabla_\tau \psi|^2 dx. \tag{2.20}

该估计的关键在于利用了修正剖面 $\theta_1$ 所满足的正交性条件 (2.18)，即 $\int_{\mathbb{R}} (\theta_0')^2 f'''(\theta_0) \theta_1\, d\rho = 0$ 。这一条件消除了潜在的不稳定模态，使得负特征值被有效控制。

3. 误差方程的精细分析：将真实解代入后，误差 $u = c_\varepsilon - c_A$ 和 $w = v_\varepsilon - v_A$ 满足一个复杂的耦合系统。作者通过以下步骤处理：

利用 Corollary 2.7 对 $$u$$ 进行分解，分离出沿界面法向的主导模态（与 $\theta_0'$ 相关）和剩余部分。
结合 Lemma 2.5 中关于局域积分的高阶估计。例如，该引理表明，若函数 $$a, b$$ 满足特定的正交条件 (2.16)，则 $\left| \int_{\Gamma_t(2\delta)} \zeta_\Gamma (f'(\rho_\varepsilon) a + \varepsilon f(\rho_\varepsilon) b)\, dx \right| \leq C \varepsilon^3 (\|a\|_{L^1} + \|b\|_{L^1}).$ 这类估计对于处理质量守恒项和粘度跳跃带来的低阶项至关重要，能确保残差足够小。
对速度误差 $$w$$ 使用标准的 Stokes 算子能量估计，并利用 $$u$$ 的高阶小性来闭合论证。

整个证明框架是一个迭代过程：先构造到 $O(\varepsilon^2)$ 阶的近似解（Theorem 3.1），再通过上述工具证明误差确实保持在 $O(\varepsilon^{N+1/2})$ 阶，从而完成收敛性证明。

4. 与之前工作的比较

本文的工作在多个维度上推进了已有的锐界面极限理论：

1. 与非质量守恒模型的比较： Abels 和 Fei 在 [1] 中首次严格处理了 Navier-Stokes/Allen-Cahn 系统（允许不同粘度）的锐界面极限，但其模型缺少质量守恒项。本文的核心挑战和创新正是源于此全局约束。质量守恒不仅改变了极限界面动力学（增加了体积保持项 $-\fint_{\Gamma_t} H_{\Gamma_t}$ ），更深刻地影响了渐近分析的结构。如文中所述，在 [1] 中 $O(\varepsilon)$ 阶的内层展开项会自动消失，但在本文中必须显式构造一个非平凡的修正剖面 $\theta_1$ 来处理它。这使得近似解的构造和残差估计都更为复杂。

2. 与常粘度或准静态模型的比较： Hensel 和 Liu [21] 使用相对熵方法处理了常粘度 ( $\nu^+ = \nu^-$ ) 的 Navier-Stokes/Allen-Cahn 系统。Kroemer 和 Laux [28] 将相对熵方法应用于质量守恒的 Allen-Cahn 方程（无流体耦合）。Abels 和 Liu [6] 则研究了准静态 Stokes 流（忽略惯性项 $\partial_t v + v\cdot\nabla v$ ）且假设等粘度的情形。本文同时放松了等粘度假设和准静态假设，处理了最一般的、具有不同粘度的、完全动态的 Navier-Stokes 流，并耦合了质量守恒的相场方程，显著拓宽了适用范围。

3. 与形式渐近分析的比较： Chen, Hilhorst, 和 Logak [13] 对质量守恒 Allen-Cahn 方程进行了形式上的匹配渐近展开，并推导出了正确的极限方程（即带体积保持的平均曲率流）。然而，他们的工作停留在形式层面，缺乏严格的收敛性证明。本文不仅验证了 [13] 的形式推导的正确性，更将其嵌入到一个完整的、严格的分析框架中，并进一步耦合了 Navier-Stokes 方程，处理了流体-界面的双向耦合效应。

4. 方法论上的比较：相对熵方法（如 [21], [28]）是一种强大而优雅的工具，尤其适用于处理弱解和较低正则性的情形。然而，它通常难以获得高阶收敛率，并且在处理像不同粘度这样复杂的系数跳跃时可能面临技术困难。本文采用的经典匹配渐近展开结合能量估计的方法，虽然技术细节更为繁琐，但能够提供任意高阶的收敛精度（由定理中的 $$N$$ 控制），并能对解在界面内外的行为给出非常精细的刻画（如 (1.17a)-(1.17c) 所示）。这是一种互补的、针对强解的高精度分析路径。

代价与局限：本文的主要代价是需要假设极限系统存在光滑解，并且要求初值与近似解高度匹配（条件 (1.16)）。这限制了结果在长时间或奇性发展情形下的直接应用。相比之下，相对熵方法或变分方法（如 varifold 解）可能在更弱的解框架下得到全局存在性结果，但牺牲了收敛精度和对解结构的精细描述。

3. Hall磁流体与电子磁流体湍流的决定性波数

Determining wavenumbers for Hall and electron magnetohydrodynamics turbulence

作者: Hassan Babaei, Mimi Dai

类别: math.AP

发布日期: 2026-04-12

arXiv ID: 2604.10804v1

摘要翻译: 在湍流中，Kolmogorov波数刻画了粘性效应占主导的最小尺度。Foias与Prodi [8] 首次引入了该概念的一个数学类比——决定性波数，用于量化唯一确定解长时间行为所需的最少模态集合。我们将这一框架从Navier-Stokes方程推广至磁化等离子体模型，重点研究亚离子尺度与耗散尺度范围内的Hall-MHD与电子-MHD湍流。

Comment: 无

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

1. 研究对象和背景

本文研究两类描述等离子体小尺度湍流的偏微分方程模型：霍尔磁流体动力学（Hall-MHD）系统与电子磁流体动力学（Electron-MHD, EMHD）系统。这些模型在空间尺度接近或小于离子惯性长度 $$d_i$$ 时变得重要，此时经典磁流体力学（MHD）中忽略的两流体效应（特别是霍尔电流项）不可忽略。

具体而言，三维不可压缩 Hall-MHD 系统为

\begin{aligned} \partial_t u + u \cdot \nabla u + b \cdot \nabla b + \nabla p &= \nu \Delta u, \\ \partial_t b + u \cdot \nabla b - b \cdot \nabla u + d_i \nabla \times \big( (\nabla \times b) \times b \big) &= \mu \Delta b, \\ \nabla \cdot u = 0, \quad \nabla \cdot b = 0, \end{aligned} \tag{1.1}

其中 $$u(x,t)$$ 为速度场， $$b(x,t)$$ 为磁场， $\nu>0$ 、 $\mu>0$ 分别为粘性系数与磁扩散系数， $$d_i>0$$ 为离子惯性长度。当 $$d_i=0$$ 时，系统退化为经典 MHD；当 $d_i \neq 0$ 时，霍尔项引入了色散性与更强的非线性耦合。

在更极端的小尺度极限下（离子视为静止背景），仅由电子运动主导磁场演化，得到 EMHD 方程：

\partial_t b + d_i \nabla \times \big( (\nabla \times b) \times b \big) = \mu \Delta b, \quad \nabla \cdot b = 0. \tag{1.2}

该模型刻画了亚离子尺度下的哨声波（whistler wave）相互作用主导的湍流，其能量级串机制与流体湍流显著不同。现象学理论预测 EMHD 磁能谱服从 $k^{-10/3}$ 标度律，而非 Kolmogorov 的 $k^{-5/3}$ 。

在湍流理论中，Kolmogorov 耗散波数 $\kappa_\eta = (\varepsilon / \nu^3)^{1/4}$ （三维情形）刻画了粘性开始主导的最小尺度，也隐含了系统有效自由度的有限性。Foias 与 Prodi 在二维 Navier–Stokes 方程中引入了“决定性模”（determining modes）的概念：若两个解在低频部分（低于某临界波数）一致，则其长时间行为趋于一致。这为湍流的有限维吸引子提供了数学基础。

近期，Cheskidov 与 Dai 等人在三维 Navier–Stokes 方程中建立了“决定性波数” $\Lambda_u(t)$ 的存在性，并证明其时间平均满足

\langle \Lambda_u \rangle \lesssim \lambda_0 + \kappa_d := \lambda_0 + \left( \frac{\varepsilon}{\nu^3} \right)^{\frac{1}{d+1}},

其中 $\varepsilon = \nu \langle \|\nabla u\|_{L^2}^2 \rangle$ 为平均耗散率， $d \in [0,3]$ 为间歇性维度（intermittency dimension）。当 $$d=3$$ 时即为经典 Kolmogorov 尺度。

然而，对于 Hall-MHD 与 EMHD 这类具有复杂非线性结构（如霍尔项 $\nabla \times ((\nabla \times b) \times b)$ ）和多重尺度耦合的等离子体模型，尚缺乏与现象学耗散尺度相匹配的决定性波数理论。此前工作（如文中提及的 [11]）虽尝试定义决定性波数，但未达到物理最优性，亦未建立与 $\kappa_e^{\delta_b} = (\varepsilon_b / \mu^3)^{1/(\delta_b - 1)}$ 等 EMHD 预测耗散波数的定量联系。

因此，本文旨在填补这一空白：在 Hall-MHD 与 EMHD 的弱解框架下，构造物理上最优的决定性波数 $\Lambda_{u,r}(t)$ 与 $\Lambda_{b,r}(t)$ ，并证明在明确的“尺度局域化间歇性假设”下，其时间平均被现象学耗散波数所控制。这不仅将决定性模理论从流体拓展至磁化等离子体，也为湍流有效自由度的数学刻画提供了新范式。

2. 主要定理或主要结果

作者在弱解与强解两类框架下建立了决定性波数的存在性及其与现象学耗散尺度的定量关系。主要结果如下：

定理 3.1（EMHD 决定性波数的存在性）：设 $b^{(1)}(t)$ 与 $b^{(2)}(t)$ 是 EMHD 方程 (3.1) 的两个弱解，定义其对应的磁决定性波数为

\Lambda_{b,r}(t) = \min\left\{ \lambda_q : (L \lambda_{p-q})^\delta \lambda_q^{n/r} \|b_p\|_{L^r} < c_r \mu, \, \forall p > q; \text{ 且 } \lambda_q^{-1} \|\nabla b_{\leq q}\|_{L^\infty} < c_r \mu, \, q \in \mathbb{N} \right\}, \tag{3.2}

其中 $r \in (n, 2n)$ ， $\delta \geq 0$ ， $$L>0$$ 为特征长度尺度， $$c_r>0$$ 为小常数。令 $\Lambda(t) = \max\{\Lambda_{b^{(1)}}(t), \Lambda_{b^{(2)}}(t)\}$ ，若两解在频率 $\leq \Lambda(t)$ 上完全一致，即 $(b^{(1)} - b^{(2)})|_{\leq \Lambda} = 0$ ，则对任意 $s \in (-n/r, n/r - 1)$ ，有

\lim_{t \to \infty} \|b^{(1)}(t) - b^{(2)}(t)\|_{H^s} = 0.

特别地，当 $$n=3$$ 且 $r \in (3,6)$ 时， $$s$$ 可取接近 0 的值，故收敛发生在接近 $$L^2$$ 的空间中。

定理 3.2（Hall-MHD 决定性波数的存在性）：设 $(u,b^{(1)})$ 与 $(v,b^{(2)})$ 是 Hall-MHD 系统 (1.1) 的两个弱解。定义速度与磁场的决定性波数分别为

\begin{aligned} \Lambda_{u,r}(t) &:= \min\left\{ \lambda_q : (L \lambda_{p-q})^\sigma \lambda_q^{-1} \|u_p\|_{L^\infty} < c_r \nu, \, \forall p > q; \, \lambda_q^{-2} \|\nabla u_{\leq q}\|_{L^\infty} < c_r \nu \right\}, \\ \Lambda_{b,r}(t) &:= \min\left\{ \lambda_q : (L \lambda_{p-q})^\delta \lambda_q^{n/r} \|b_p\|_{L^r} < c_r \mu, \, \forall p > q; \, \lambda_q^{-1} \|\nabla b_{\leq q}\|_{L^\infty} < c_r \mu \right\}, \end{aligned}

其中 $\sigma = (\gamma - 1)/2$ ， $0 < \gamma \leq 3$ ，故 $-1/2 \leq \sigma \leq 1$ 。令 $\Lambda_U(t) = \max\{\Lambda_u^{(1)}, \Lambda_u^{(2)}\}$ ， $\Lambda_B(t) = \max\{\Lambda_{b^{(1)}}, \Lambda_{b^{(2)}}\}$ 。若

(u - v)|_{\leq \Lambda_U} = 0, \quad (b^{(1)} - b^{(2)})|_{\leq \Lambda_B} = 0,

且满足参数条件 $\delta > 1 - n/r$ ， $\sigma > 1 - n/r$ ，则对任意

s \in \left( -\min\left\{ \frac{n}{r}, \delta, \sigma \right\}, \, \frac{n}{r} - 1 \right),

有

\lim_{t \to \infty} \left( \|u(t) - v(t)\|_{H^s} + \|b^{(1)}(t) - b^{(2)}(t)\|_{H^s} \right) = 0.

进一步，在尺度局域化间歇性假设下（见第5节），作者证明了决定性波数的时间平均被现象学耗散波数所控制：

对 EMHD，有 $\langle \Lambda_{b,r} \rangle \lesssim \lambda_0 + \kappa_e^{\delta_b}$ ，其中 $\kappa_e^{\delta_b} := \left( \frac{\varepsilon_b}{\mu^3} \right)^{\frac{1}{\delta_b - 1}}, \quad \varepsilon_b := \mu \lambda_0^{\delta_b} \left\langle \sum_{q \geq -1} \lambda_q^2 \|b_q\|_{L^2}^2 \right\rangle;$
对 Hall-MHD，有 $\langle \Lambda_{u,r} \rangle \lesssim \lambda_0 + \kappa_u^{\delta_u}, \quad \langle \Lambda_{b,r} \rangle \lesssim \lambda_0 + \kappa_e^{\delta_b},$ 其中 $\kappa_u^{\delta_u}$ 为 Navier–Stokes 型耗散波数， $\delta_u$ 为速度场间歇性维度。

此外，对强 EMHD 解（即 $b \in L^\infty(0,\infty; B_{2,\infty}^s)$ ， $$s>0$$ ），作者在第6节证明了磁决定性波数存在一致上界：

\Lambda_{b,r}(t) \leq C(\|b\|_{B_{2,\infty}^s}), \quad \forall t \geq 0,

这表明 Besov 正则性可直接控制决定性波数的大小。

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

本文的核心技术在于将 Navier–Stokes 方程中发展的决定性波数框架（基于精细的 Littlewood–Paley 分解与尺度分离）适配到 Hall-MHD 与 EMHD 的复杂非线性结构中，并克服霍尔项带来的新困难。

第一，决定性波数的物理最优定义。作者并未沿用早期文献中粗糙的定义，而是借鉴 [5,6] 中对 Navier–Stokes 的处理，对速度场采用与 NS 相同的尺度不变量（如 $\lambda_q^{-1} \|u_q\|_{L^\infty}$ ），对磁场则根据 EMHD 的标度性质引入 $\lambda_q^{n/r} \|b_q\|_{L^r}$ 。这种定义确保了决定性波数与现象学耗散尺度在量纲上一致，是后续比较的基础。

第二，霍尔项的非线性估计与对消结构利用。霍尔项 $\nabla \times ((\nabla \times b) \times b)$ 在 dyadic 分解下产生高频-低频相互作用。作者通过 Bony 仿微分分解将其拆分为三部分：低频驱动高频（transport-like）、高频驱动低频、以及可比拟频率相互作用。关键在于处理形如

\int \Delta_q \big( (\nabla \times h) \times b \big) \cdot \nabla \times h_q \, dx

的项（见 (4.4) 中的 $$J$$ 和 $$K$$ ）。为此，作者引入了霍尔交换子（Hall commutator）：

[\Delta_q, b_{\leq p-2} \times \nabla \times] h_p = \Delta_q(b_{\leq p-2} \times (\nabla \times h_p)) - b_{\leq p-2} \times (\nabla \times \Delta_q h_p), \tag{2.2}

并证明了其满足类似于输运交换子的估计（引理 2.5）：

\|[\Delta_q, b_{\leq p-2} \times \nabla \times] h_p\|_{L^{r_1}} \lesssim \|h_p\|_{L^{r_2}} \sum_{p' \leq p-2} \lambda_{p'} \|b_{p'}\|_{L^{r_3}}, \quad \frac{1}{r_1} = \frac{1}{r_2} + \frac{1}{r_3}.

这一估计依赖于 $\nabla \cdot b = 0$ 所提供的对消，是控制非线性项的关键。

第三，尺度局域化间歇性假设下的求和技巧。在估计非线性项（如 $$J_1, J_2, J_3$$ ）时，作者利用决定性波数定义 (3.2) 中的不等式（如 $\lambda_q^{n/r} \|b_q\|_{L^r} < c_r \mu \lambda_Q^\delta \lambda_q^{-\delta}$ 对 $$q>Q$$ ）将高频模的大小用 $\lambda_Q$ （即 $\Lambda(t)$ ）控制。随后通过精心选择指数（如利用 $\delta \geq 0$ 、 $$s - n/r < 0$$ ）并应用 Jensen 不等式，将双重求和

\sum_{q>Q} \sum_{Q<p'\leq q} \cdots

转化为单重求和 $\sum_q \lambda_q^{2s+2} \|h_q\|_{L^2}^2$ ，从而与耗散项 $\mu \sum_q \lambda_q^{2s+2} \|h_q\|_{L^2}^2$ 匹配。这一过程要求对频率指数进行精细平衡，是证明收敛性的核心。

第四，强解情形下 Besov 空间的应用。在第6节，作者利用强解的 Besov 正则性 $b \in B_{2,\infty}^s$ （ $$s>0$$ ）直接估计 $\|b_q\|_{L^r}$ 与 $\|\nabla b_{\leq q}\|_{L^\infty}$ ，从而给出 $\Lambda_{b,r}(t)$ 的显式上界。这展示了函数空间正则性与决定性波数之间的直接联系，为数值模拟中估计有效自由度提供了理论依据。

综上，作者通过结合精细的频率局部化分析、针对霍尔项定制的交换子估计、以及尺度不变量的巧妙运用，成功将决定性模理论拓展至磁化等离子体湍流模型。

4. 与之前工作的比较

本文相较于已有工作在多个维度上实现了实质性推进：

1. 与经典 Navier–Stokes 决定性模理论的对比：Cheskidov 与 Dai 在 [5,6] 中为三维 Navier–Stokes 方程建立了决定性波数 $\Lambda_u(t)$ ，并证明 $\langle \Lambda_u \rangle \lesssim \kappa_d$ 。本文将此框架推广至 Hall-MHD 与 EMHD，但面临两大新挑战：(i) 磁场方程包含霍尔项这一高阶非线性项，其结构不同于对流项；(ii) 系统包含两个耦合场（速度与磁场），需分别定义决定性波数。作者通过为磁场设计新的尺度不变量（基于 $$L^r$$ 范数而非 $L^\infty$ ）并发展霍尔交换子估计，成功克服了这些障碍。值得注意的是，速度场的决定性波数定义直接沿用 [5,6] 的形式，这反映了 Hall-MHD 中速度方程与 NS 方程在标度上的相似性。

2. 与先前 Hall/EMHD 决定性模工作的对比：文中明确指出，此前关于 Hall-MHD 决定性模的工作（引用 [11]）存在两大缺陷：(i) 定义的决定性波数非物理最优，即其量纲或增长速率与现象学预测不符；(ii) 未建立与耗散波数的定量联系。本文通过引入尺度局域化间歇性假设，首次严格证明了 $\langle \Lambda_{b,r} \rangle \lesssim \kappa_e^{\delta_b}$ ，填补了这一空白。这一比较凸显了本文定义的优越性：只有采用正确的尺度不变量，才能实现与物理预测的一致。

3. 与现象学湍流理论的对接：等离子体湍流的现象学理论（如 [9]）长期预测了 Kolmogorov-like 耗散波数 $\kappa_e^{\delta_b}$ 的存在，但缺乏严格的数学支撑。本文首次在弱解框架下证明，决定性波数（一个纯数学概念）的时间平均被这些现象学尺度所控制。这为“有效自由度由耗散尺度决定”这一物理直觉提供了首个针对 Hall/EMHD 的严格数学验证。特别地，当间歇性维度 $\delta_b = 3$ 时， $\kappa_e^{\delta_b} = (\varepsilon_b / \mu^3)^{1/2}$ ，与 EMHD 的经典预测一致。

4. 适用范围与假设的权衡：本文结果适用于弱解（Leray–Hopf 型），这比要求强解（如 $$H^s$$ , $$s>5/2$$ ）的工作更具普适性，因为三维 Hall-MHD 的全局强解存在性仍是公开问题。代价是需要引入尺度局域化间歇性假设来控制高频模。这一假设比全局正则性假设弱，但强于无任何先验信息的情形。作者在第6节补充了强解情形下的结果（一致有界性），展示了不同正则性假设下的结论差异。

5. 技术工具的演进：相比早期使用 Galerkin 截断或简单能量估计的方法，本文全面采用现代调和分析工具：Littlewood–Paley 分解、Bony 仿微分分解、Bernstein 不等式、以及定制的交换子估计。特别是霍尔交换子估计（引理 2.5）是本文的技术创新点，此前在类似 PDE 文献中未见报道。这使得作者能够精确追踪非线性相互作用的频率转移，从而获得最优估计。

综上，本文不仅将决定性模理论拓展至新的物理模型，更通过精细的分析框架建立了数学理论与物理现象之间的定量桥梁，代表了该方向的重要进展。

4. 一类修正NSCH-Oldroyd系统的局部适定性：基于PINN的数值验证

Local Well-Posedness of a Modified NSCH-Oldroyd System: PINN-Based Numerical Illustrations

作者: Woojeong Kim

类别: math.AP, math.NA

发布日期: 2026-04-12

arXiv ID: 2604.10616v1

摘要翻译: 受血栓建模的启发，我们研究了一类修正的Navier-Stokes-Cahn-Hilliard-Oldroyd系统，并对该修正系统进行了基于物理信息神经网络（PINN）的数值验证。为便于分析，我们在保持相应耗散能量结构的前提下，对形变变量引入了一个扩散增强系统。我们证明了该新系统的局部适定性。此外，我们还针对典型血栓情形给出了基于PINN的数值验证，并报告了基于能量衰减的Metropolis-Hastings采样所获得的残差损失与基准误差。

Comment: 58pages, 14 figures

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

1. 研究对象和背景

本文研究一个源于血栓（thrombus）建模的耦合偏微分方程系统，该系统结合了Navier–Stokes–Cahn–Hilliard（NSCH）模型与Oldroyd-B型粘弹性流体模型。具体而言，所考察的未知量包括速度场 $$u(x,t)$$ 、相场变量 $$ϵ(x,t)$$ 、形变梯度张量 $$F(x,t)$$ 和压力 $$p(x,t)$$ ，定义在有界光滑区域 $$Ω \subset ℝd$$ （ $$d=2,3$$ ）上，并满足如下控制方程组：

\begin{cases} ρ(\partial_t u + u \cdot \nabla u) + \nabla p - \nabla \cdot (\eta(ϵ)\nabla u) = -\lambda \nabla \cdot (\nabla ϵ \otimes \nabla ϵ) + \nabla \cdot \big( \nu(ϵ)(FF^\top - I) \big) - \dfrac{\eta(ϵ)(1-ϵ)}{\kappa(ϵ)} u, \\ \nabla \cdot u = 0, \\ \partial_t F + u \cdot \nabla F - k\big( \nu(ϵ)\Delta F + 2\nabla \nu(ϵ) \cdot \nabla F \big) = \nabla u\, F, \\ \partial_t ϵ + u \cdot \nabla ϵ = \tau \Delta \mu, \\ \mu = -\lambda \Delta ϵ + \lambda \gamma f'(ϵ) + \dfrac{\nu'(ϵ)}{2} \mathrm{tr}(FF^\top - I), \end{cases}

其中 $$ϵ=1$$ 表示血液区域， $$ϵ=0$$ 表示血栓区域，中间值对应混合界面；函数 $f(ϵ) = \frac{(ϵ-1)^2ϵ^2}{4h^2}$ 是标准双阱势；参数 $ρ, \eta, \kappa, \nu, \lambda, \gamma, \tau, k$ 均为正的物理常数或光滑系数函数，满足一致上下界条件 $\alpha \leq \eta(x), \nu(x), \kappa(x) \leq \beta$ 。

该系统的物理动机来自对血栓形成过程中流体-结构相互作用的连续介质描述。已有工作如 [6,25,26] 提出了类似的NSCH-Oldroyd模型，并用于数值模拟，但缺乏严格的适定性理论。特别地，文献 [6] 首次证明了不含扩散项的原始形变方程（即令 $$k=0$$ ）下强解的局部存在唯一性，并建立了能量耗散结构。然而，由于形变梯度 $$F$$ 的演化方程缺少二阶扩散项，其分析面临两大困难：一是高阶导数估计中缺乏对 $\Delta F$ 的控制，导致无法闭合先验估计；二是在数值实现中界面附近梯度剧烈变化时稳定性差。

为克服上述障碍，本文引入一个修正模型，在形变方程中添加一个小的扩散项 $-k\nu(ϵ)\Delta F$ （以及伴随的 $-2k\nabla\nu(ϵ)\cdot\nabla F$ 项以保持变分结构），同时严格保持原系统的能量耗散性质。这一修改虽看似微小，但在分析上至关重要：它使得 $$F$$ 满足一个拟线性抛物型方程，从而允许使用椭圆正则性理论控制 $$H^3$$ 范数。本文的核心目标即是针对此扩散增强型NSCH-Oldroyd系统，建立其局部强解的存在唯一性，并说明该框架为后续数据同化与PINN模拟提供更稳健的基础。

2. 主要定理或主要结果

作者在空间维数 $$d=2,3$$ 下证明了修正系统的局部强解适定性，主要结论见定理1.1。具体而言，在如下初始正则性假设下：

$u_0 \in D(A) = H^2(\Omega)^d \cap V$ （其中 $A=-P\Delta$ 为Stokes算子， $$V$$ 为无散向量场空间），
$\u03f5_0 \in H^5(\Omega)$ 且满足齐次Neumann边界条件 $\partial_n \u03f5_0 = 0$ ，
$F_0 \in H^3(\Omega)^{d\times d}$ 且满足 $\partial_n F_0^{ij} = \partial_n \Delta F_0^{ij} = 0$ （对所有 $$i,j$$ ），

则存在时间 $$T_0 > 0$$ （依赖于初值范数），使得系统 (1.2)–(1.3) 在区间 $$[0,T_0]$$ 上存在唯一解 $(u, F, \u03f5, p)$ ，且满足以下正则性：

\begin{aligned} u &\in C([0,T_0]; D(A)) \cap L^2(0,T_0; H^3(\Omega)^d), \\ \partial_t u &\in C([0,T_0]; H) \cap L^2(0,T_0; D(A)), \\ F &\in C([0,T_0]; H^2(\Omega)^{d\times d}) \cap H^2(0,T_0; H^3(\Omega)^{d\times d}) \cap H^1(0,T_0; H^2(\Omega)^{d\times d}), \\ ϵ &\in C([0,T_0]; H^3(\Omega)) \cap L^2(0,T_0; H^4(\Omega)) \cap H^1(0,T_0; H^3(\Omega)), \\ \nabla p &\in L^2(0,T_0; H^1(\Omega)). \end{aligned}

此外，作者通过能量测试推导出系统的耗散结构。定义总能量为

E(t) = \int_\Omega \left( |u|^2 + \lambda|\nabla\u03f5|^2 + 2\lambda\gamma f(\u03f5) + \nu(\u03f5)\,\mathrm{tr}(FF^\top - I) \right) dx,

则在适当正则性下，有

\frac{d}{dt} E(t) = -2\int_\Omega \left( \eta(\u03f5)|\nabla u|^2 + \tau|\nabla\mu|^2 + k\nu^2(\u03f5) \sum_{i,j} |\nabla F^{ij}|^2 + \frac{\eta(\u03f5)(1-\u03f5)}{\kappa(\u03f5)} |u|^2 \right) dx \leq 0.

这表明扩散项 $k\nu(\u03f5)\Delta F$ 的引入并未破坏原系统的能量耗散机制，反而提供了额外的耗散项 $k\nu^2(\u03f5)|\nabla F|^2$ ，这对先验估计至关重要。

作为补充，作者还展示了基于物理信息神经网络（PINN）的数值实验，针对代表性血栓情形，利用Metropolis–Hastings采样结合能量衰减构造损失函数，报告了残差损失与基准误差，验证了修正系统在界面区域的数值稳定性。

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

本文证明的核心在于构建适用于修正系统的高阶先验估计，并通过Galerkin逼近结合紧性方法完成局部适定性证明。关键创新点和工具如下：

1. 扩散增强带来的椭圆正则性提升：原始系统（ $$k=0$$ ）中 $$F$$ 满足输运型方程 $\partial_t F + u\cdot\nabla F = \nabla u\, F$ ，仅能通过沿迹估计获得有限正则性。而本文引入的扩散项 $-k\nu(\u03f5)\Delta F$ 使 $$F$$ 方程变为拟线性抛物型。作者在第4页引用并推广了文献[6]中的两个关键引理：

引理2.1（Stokes型正则性）：对含相场依赖粘性的Stokes方程 $-\nabla\cdot(\eta(\u03f5)\nabla u) + \nabla p = f$ ，若 $\u03f5 \in H^3$ ，则解 $u \in H^3$ ，且有估计 $|u|_{H^3} + \left| \frac{p}{\eta(\u03f5)} \right|_{H^1} \leq C |f|_{H^1} \left( 1 + (1 + |\u03f5|_{H^2}^2)(|\u03f5|_{H^2}^2 + |\u03f5|_{H^2}^{1/2} + |\u03f5|_{H^3}^{1/2}) \right).$
引理2.2（四阶椭圆正则性）：对双调和Neumann问题 $\Delta^2 \u03f5 = f$ ， $\partial_n \u03f5 = \partial_n \Delta \u03f5 = g$ ，若兼容性条件成立，则 $\u03f5 \in H^4$ ，且 $|\u03f5|_{H^4} \leq \tilde{C} \left( |f|_{L^2} + |g|_{H^{1/2}(\partial\Omega)} + \left| \int_\Omega \u03f5 \, dx \right| \right).$ 这两条引理使得在Galerkin逼近中可分别控制 $$u$$ 和 $\u03f5$ 的 $$H^3$$ 与 $$H^4$$ 范数。

2. 能量结构的精细保持：作者通过多重测试函数组合（ $$u$$ , $\mu$ , $F\nu(\u03f5)$ ）推导出封闭的能量恒等式 (2.37)–(2.39)。特别地，在测试形变方程时，采用内积 $(\cdot, F\nu(\u03f5))$ 并取迹，巧妙地将 $\mathrm{tr}(\nabla u\, FF^\top)\nu(\u03f5)$ 项与动量方程中的对应项抵消，从而保留了耗散结构。新增的扩散项贡献了负定项 $-k\int \nu^2(\u03f5) |\nabla F|^2 dx$ ，这是闭合高阶估计的关键。

3. 高阶先验估计的构造：在第10页，作者定义复合量

Z := |\nabla u|^2 + |\nabla \partial_t u|^2 + |\Delta \u03f5|^2 + |\nabla \partial_t \u03f5|^2 + |\Delta F|^2 + \sum_{i,j} |\nabla F^{ij}|^2 + \cdots,

并通过测试动量方程与 $$Au$$ （Stokes算子作用于 $$u$$ ）来估计 $$|Au|$$ 。利用Sobolev嵌入 $H^1 \hookrightarrow L^6$ , $H^{1/2} \hookrightarrow L^3$ 及Young不等式，将非线性项如 $(u\cdot\nabla u, Au)$ 、 $(\eta'(\u03f5)\nabla\u03f5\cdot\nabla u, Au)$ 控制为 $\alpha|Au|^2/12 + C(\text{低阶项})$ ，从而吸收主部。类似技巧用于处理 $\u03f5$ 和 $$F$$ 的时间导数方程。

4. 边界条件与平均值处理：Neumann边界条件 $\partial_n \u03f5 = \partial_n \Delta \u03f5 = 0$ 允许使用广义Poincaré不等式（如 (2.12)–(2.17)）将 $$H^1$$ 或 $$H^2$$ 范数控制为高阶梯度范数。例如， $|\nabla \u03f5|_{H^1} \leq |\Delta \u03f5|_{L^2}$ 。此外，相场平均值 $\langle \u03f5 \rangle_\Omega$ 守恒（见 (2.5)），避免了零特征值问题。

综上，作者通过引入扩散项重构系统，再结合椭圆正则性引理、精细能量估计与高阶Galerkin逼近，完成了局部强解的存在唯一性证明。

4. 与之前工作的比较

本文与既有文献的主要区别体现在模型构造、分析框架与适用范围三个层面：

1. 与原始NSCH-Oldroyd模型 [6] 的比较：

模型差异：文献 [6] 研究的系统中形变方程为纯输运型 $\partial_t F + u\cdot\nabla F = \nabla u\, F$ （即 $$k=0$$ ）。本文则引入扩散项 $-k\nu(\u03f5)\Delta F$ （ $$k>0$$ ），构成“扩散增强”系统。
分析推进：尽管 [6] 已证明局部适定性，但其证明依赖于对输运方程沿迹估计的精细处理，且无法自然获得 $F \in H^3$ 的正则性。本文通过添加扩散项，将 $$F$$ 方程转化为抛物型，从而可直接应用椭圆正则性理论（引理2.2的变体），显著简化了高阶估计的构造。更重要的是，新系统具有更强的耗散性（多出 $k\nu^2|\nabla F|^2$ 项），这在长时间行为或扰动分析中可能更具优势。
代价与局限：引入扩散项是对物理模型的微小修改。虽然作者论证其保持了能量结构，但严格来说已非原始Oldroyd-B模型。因此，本文结果不能直接推出 [6] 中系统的适定性，而是提供了一个更易处理的近似模型。不过，当 $k \to 0$ 时，修正系统应收敛于原系统，这为未来研究提供了桥梁。

2. 与一般NSCH或Oldroyd系统文献的比较：

耦合复杂性：标准NSCH系统（如 [27]）仅耦合Navier-Stokes与Cahn-Hilliard方程，而Oldroyd-B模型单独研究粘弹性应力。本文系统同时包含三者耦合，且形变梯度 $$F$$ 与相场 $\u03f5$ 通过系数 $\nu(\u03f5)$ 和化学势 $\mu$ 强耦合，远比单一模型复杂。
正则性要求：为处理耦合非线性项（如 $\nu(\u03f5)\mathrm{tr}(FF^\top)$ 出现在 $\mu$ 中），本文要求初值 $\u03f5_0 \in H^5$ 、 $F_0 \in H^3$ ，这比标准NSCH（通常 $\u03f5_0 \in H^2$ ）或Oldroyd（ $F_0 \in H^2$ ）更高。这是由耦合项的高阶导数需求决定的，属于必要代价。
边界条件处理：本文采用齐次Neumann边界条件 $\partial_n \u03f5 = \partial_n \Delta \u03f5 = 0$ 和 $\partial_n F = \partial_n \Delta F = 0$ ，这与多数文献一致，但作者通过引理2.3处理了边界迹的负范数估计（如 $|\partial_n F|_{H^{-1/2}(\partial\Omega)} \leq C(|F|_{H^1} + |\Delta F|_{H^{-1}})$ ），这是处理高阶Neumann问题的关键技术细节，未见于早期工作。

3. 数值方法层面的推进：

与依赖参考数据的模拟（如 [30,31]）不同，本文的PINN实验仅基于物理方程残差和能量衰减构造损失函数，未使用外部数据监督。这验证了修正系统本身具有更好的数值稳定性，尤其在界面梯度剧烈区域。虽然数值部分非本文重点，但它为分析结果提供了支撑：若系统不适定，PINN难以收敛。

总体而言，本文并非在经典适定性理论上取得突破（如全局存在性或弱解唯一性），而是针对一个具体生物物理模型，通过巧妙的模型修正，构建了一个兼具物理合理性与数学可处理性的新框架，填补了“严格适定性理论”与“复杂多物理场模拟”之间的空白。

5. 旋转球面上二维欧拉方程中四面体非带状流的分歧

Bifurcation of Tetrahedral Non-Zonal Flows in the 2D Euler Equations on a Rotating Sphere

作者: Yuri Cacchiò

类别: math.AP

发布日期: 2026-04-11

arXiv ID: 2604.10092v1

摘要翻译: 我们研究了球面 $\mathbb{S}^2$ 上由二维欧拉方程的定常解所引发的有限振幅非带状流的出现。通过将拉普拉斯-贝尔特拉米特征空间限制到四面体对称群 $\mathbf{T}$ 的不变子空间，我们绕过了 $$(2l+1)$$ 维的核退化，从而获得了一个标量形式的李亚普诺夫-施密特约化。我们分析了四种不同的物理非线性模型：多项式模型、正弦-戈登模型、双曲正弦-戈登模型以及指数型（刘维尔）模型。我们通过谱投影显式地导出了分歧参数，并证明了分歧拓扑（亚临界或超临界）并非几何不变量，而是由非线性的奇偶性及质量守恒所决定。

Comment: 无

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

1. 研究对象和背景

本文研究定义在单位球面 $\mathbb{S}^2$ 上的二维不可压缩欧拉方程的定常解，特别关注由对称性破缺产生的有限振幅非带状流（non-zonal flows）。设 $\psi = \psi(\varphi, \theta)$ 为流函数，其中 $\varphi \in [0, 2\pi)$ 为经度， $\theta \in (-\pi/2, \pi/2)$ 为纬度， $x_3 = \sin\theta$ 为极轴坐标。在以角速度 $\omega$ 绕 $$x_3$$ 轴旋转的参考系中，相对涡度由球面拉普拉斯–贝尔特拉米算子给出：

\Delta \psi = \frac{1}{\cos\theta} \partial_\theta (\cos\theta \, \partial_\theta \psi) + \frac{1}{\cos^2\theta} \partial_{\varphi}^2 \psi.

定常欧拉方程可写为

J(\psi, \Delta \psi + 2\omega x_3) = 0,

其中雅可比行列式 $J(f,g) = \frac{1}{\cos\theta}(\partial_\varphi f \, \partial_\theta g - \partial_\theta f \, \partial_\varphi g)$ 。在梯度非退化区域，绝对涡度是流函数的函数，从而导出半线性椭圆方程

-\Delta \psi + F(\psi) = 2\omega x_3.

当 $\omega = 0$ （无旋转极限）时，方程简化为

-\Delta \psi + F(\lambda, \psi) = 0 \quad \text{on } \mathbb{S}^2,

其中 $F(\lambda, 0) = 0$ ， $\lambda \in \mathbb{R}$ 为控制非线性的分歧参数。此模型虽忽略科里奥利力的几何不对称性，但通过变换 $\psi_\omega(\varphi, \theta, t) = \psi_0(\varphi + \omega t, \theta) + \omega x_3$ ，其定常解可映射为旋转情形下的西传罗斯贝波。

关键困难在于球面拉普拉斯算子的本征空间具有高维简并性：对球谐次数 $l \geq 1$ ，对应本征值 $$l(l+1)$$ 的核空间维数为 $$2l+1$$ 。这使得经典一维分歧理论（如 Crandall–Rabinowitz 定理）无法直接应用，因为分歧方程变为 $$2l+1$$ 个耦合非线性方程。为克服此障碍，作者采用对称性约化策略：将函数空间限制在四面体对称群 $\mathbf{T} \subset O(3)$ 的不变子空间内。该群对应于正四面体的对称操作，在 $$l=3$$ 的七维本征空间中，其表示分解后仅含一个完全不变的一维子空间，由实四面体球谐函数

Y^*(\varphi, \theta) = \sqrt{\frac{105}{16\pi}} \sin\theta \cos^2\theta \cos(2\varphi)

张成。由此，核空间降维至一维，允许进行标量 Liapunov–Schmidt 约化。

本文试图填补的关键空白是：分歧拓扑（亚临界或超临界）是否由几何对称性决定？ 已有工作（如 Constantin–Germain）证明了分歧分支的存在性，但未刻画其局部方向。作者通过分析四类物理非线性模型——广义多项式、sine-Gordon、sinh-Gordon 和 Liouville 指数模型——揭示分歧类型并非几何不变量，而是由非线性项的奇偶性与质量守恒约束共同决定。

2. 主要定理或主要结果

作者在四面体对称不变子空间 $X_{\mathbf{T}}$ 中，对四类非线性模型分别进行了三阶渐近展开，得到以下主要结果：

通用领头阶结构：对任意满足 $F(\lambda, 0)=0$ 的光滑非线性，分歧解在 $O(\varepsilon)$ 阶具有统一形式：
$\psi(\varphi, \theta) = \varepsilon Y^*(\varphi, \theta) + O(\varepsilon^2),$
其中 $$Y^*$$ 如上定义。该流在赤道和极点为零，涡旋位于纬度 $\theta = \pm \arctan(1/\sqrt{2}) \approx \pm 35.26^\circ$ ，经度 $0^\circ, 90^\circ, 180^\circ, 270^\circ$ ，形成正负交替的四面体构型，且关于赤道反对称： $Y^*(-\theta) = -Y^*(\theta)$ 。
多项式模型：取 $F(\lambda, \psi) = P(\lambda + \psi) - P(\lambda)$ ，其中 $P(\lambda) = \mu_1 \lambda^3 - (\mu + 12)\lambda$ ， $\mu, \mu_1 > 0$ 。则：
- 分歧点为 $\lambda^* = \sqrt{\mu/(3\mu_1)}$ ；
- 一阶参数修正 $\lambda_1 = 0$ （因 $$(Y^*)^3$$ 奇性导致积分消失）；
- 二阶空间修正 $\psi_2 = -3\mu_1 \lambda^* L^{-1}((Y^*)^2)$ ，其中 $L = -\Delta - 12I$ ， $$(Y^*)^2$$ 展开含 $$l=0,2,4,6$$ 偶次球谐， $$l=0$$ 分量用于满足质量约束 $\int_{\mathbb{S}^2} F \, d\sigma = 0$ ；
- 二阶参数修正为 $\lambda_2 = -\frac{1}{6\lambda^*} \left[ \int_{\mathbb{S}^2} (Y^*)^4 \, d\sigma - 18\mu_1 (\lambda^*)^2 \int_{\mathbb{S}^2} (Y^*)^2 L^{-1}((Y^*)^2) \, d\sigma \right].$ 利用 $18\mu_1 (\lambda^*)^2 = 6\mu$ ，可进一步简化。若 $\lambda_2 < 0$ ，则为超临界分歧；若 $\lambda_2 > 0$ ，则为亚临界。
sine-Gordon 模型： $F(\lambda, \psi) = \lambda (\sin \psi - \psi)$ 。由于 $$F$$ 为奇函数，质量约束自动满足（无需 $$l=0$$ 模式）。分歧点 $\lambda^* = 12$ ，且 $\lambda_1 = 0$ （因 $$(Y^*)^3$$ 积分为零）。二阶修正 $\psi_2$ 由 $L\psi_2 = -\frac{\lambda^*}{2} (Y^*)^2$ 决定，同样含偶次球谐。三阶投影给出 $\lambda_2$ 表达式，其符号决定分歧类型。
sinh-Gordon 模型： $F(\lambda, \psi) = \lambda (\sinh \psi - \psi)$ 。同为奇非线性，质量约束自动满足。分歧点 $\lambda^* = 12$ ， $\lambda_1 = 0$ ， $\psi_2$ 由 $L\psi_2 = -\frac{\lambda^*}{2} (Y^*)^2$ 给出，与 sine-Gordon 在二阶相同，但三阶非线性项不同，导致 $\lambda_2$ 表达式差异。
Liouville 指数模型： $F(\lambda, \psi) = \lambda (e^\psi - 1 - \psi - m)$ ，其中 $m = \frac{1}{4\pi} \int_{\mathbb{S}^2} (e^\psi - 1 - \psi) \, d\sigma$ 为引入的质量位移项，用于强制满足 $\int F \, d\sigma = 0$ 。此时非线性为正，必须显式修正以满足全局约束。分歧点 $\lambda^* = 12$ ，但高阶展开涉及指数函数的泰勒展开，导致 $\psi_2, \psi_3$ 的表达式更复杂，且 $\lambda_2$ 的计算需处理非局部项 $$m$$ 。

核心结论：分歧拓扑（ $\lambda_2$ 的符号）依赖于非线性项的奇偶性与质量守恒机制。奇非线性（sine/sinh-Gordon）自动满足约束，而偶或正非线性（多项式、Liouville）需激发 $$l=0$$ 模式来平衡，这显著影响高阶相互作用与 $\lambda_2$ 的符号。

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

本文的核心技术障碍在于球面拉普拉斯算子本征空间的高维简并性（ $\dim \ker = 2l+1$ ），这使得标准分歧理论失效。作者通过结合等变分歧理论与Liapunov–Schmidt 约化，构建了一个适用于高对称性问题的分析框架。具体创新点如下：

四面体对称性约化：作者没有尝试处理完整的 $$2l+1$$ 维分歧方程组，而是将 Sobolev 空间 $H^s(\mathbb{S}^2)$ 限制到四面体群 $\mathbf{T}$ 的不变子空间 $X_{\mathbf{T}}$ 。利用群表示论， $$l=3$$ 的七维表示在 $\mathbf{T}$ 作用下分解出唯一的一维平凡子表示，由显式给出的四面体球谐 $$Y^*$$ 张成。这一约化不仅将核空间降至一维，满足 Crandall–Rabinowitz 定理的前提，而且由对称临界原理保证：在 $X_{\mathbf{T}}$ 中找到的解自动是原方程在整个球面上的解。
带约束的 Liapunov–Schmidt 展开：作者系统地执行了三阶渐近展开 $\psi = \varepsilon \psi_1 + \varepsilon^2 \psi_2 + \varepsilon^3 \psi_3 + O(\varepsilon^4)$ ， $\lambda = \lambda^* + \varepsilon \lambda_1 + \varepsilon^2 \lambda_2 + O(\varepsilon^3)$ 。关键步骤是逐阶求解线性方程 $L \psi_k = \text{已知源项}$ ，其中 $L = -\Delta + \partial_\psi F(\lambda^*, 0)I$ 。由于 $$L$$ 在像空间 $\mathrm{Im}(L)$ 上是同构，只要源项与核正交，解就存在唯一。作者精确计算了各阶源项（如 $$(Y^*)^2$$ , $Y^* \psi_2$ , $$(Y^*)^3$$ 等）在核上的投影，这些投影决定了 $\lambda_k$ 。
全局质量约束的精细处理：球面上的积分约束 $\int_{\mathbb{S}^2} F(\lambda, \psi) \, d\sigma = 0$ 是一个全局条件，源于 $\int \Delta \psi \, d\sigma = 0$ 。作者敏锐地指出，该约束的满足机制因非线性奇偶性而异：
- 对奇非线性（ $F(-\psi) = -F(\psi)$ ），若解保持反对称，则约束自动成立，故 $\psi_2$ 等高阶修正不含 $$l=0$$ 模式。
- 对偶或正非线性， $$(Y^*)^2$$ 的积分为正，必须由 $\psi_2$ 中的 $$l=0$$ 分量（即常数项）来抵消，以保证总积分为零。这要求在求解 $L \psi_2 = -\text{源项}$ 时，明确分离出 $$l=0$$ 分量，并验证其与 $$L$$ 的可逆性（因 $$L$$ 在 $$l=0$$ 上的作用为 $-\partial_\psi F(\lambda^*, 0) \neq 0$ ，故可逆）。
显式谱投影计算：为确定 $\lambda_2$ 的符号，作者需要计算形如 $\int (Y^*)^2 L^{-1}((Y^*)^2) \, d\sigma$ 的非局部积分。这涉及将 $$(Y^*)^2$$ 展开为球谐级数（仅含 $$l=0,2,4,6$$ ），然后对每个分量应用 $L^{-1}$ （即除以对应的本征值偏移量 $$l(l+1) - 12$$ ），最后再与 $$(Y^*)^2$$ 做内积。虽然论文未给出最终数值，但框架清晰表明 $\lambda_2$ 的符号由这些谱系数的加权和决定，从而将分歧拓扑问题转化为具体的谱分析问题。

综上，作者的真正创新在于将抽象的对称性约化与具体的渐近分析相结合，通过精细处理全局约束，揭示了非线性结构如何通过高阶谐波相互作用影响分歧方向。

4. 与之前工作的比较

本文与已有工作的比较可从多个维度展开：

与 Constantin–Germain (2019) 的比较：Constantin 和 Germain 在其工作中首次将四面体对称性用于球面欧拉方程的分歧分析，并利用 Rabinowitz 全局分歧定理证明了非带状流分支的存在性。然而，他们的方法本质上是非构造性的，无法提供分歧点附近的局部信息，特别是无法判断分歧是亚临界还是超临界。本文则通过 Liapunov–Schmidt 约化和高阶渐近展开，首次给出了分歧拓扑的显式判据，并揭示其依赖于非线性项的具体形式。这是对 Constantin–Germain 工作的重要深化和补充。
与一般球面分歧理论的比较：传统处理球面高维核的方法包括：(a) 引入小扰动打破对称性；(b) 使用 Lyapunov–Schmidt 约化但保留多维核，导致复杂的向量值分歧方程。本文采用的对称性子空间限制法是一种更优雅且物理上合理的策略。它不破坏问题的内在对称性，而是利用对称性来简化问题。这种方法在物理学中常见（如寻找具有特定晶格对称性的解），但在严格的 PDE 分歧分析中，本文对其在四面体情形下的应用提供了完整的数学实现。
关于分歧拓扑的理解：以往许多分歧研究（尤其在流体力学中）往往默认或仅考虑单一类型的非线性（如立方非线性），并得出特定的分歧类型。本文通过系统比较四类截然不同的物理模型，颠覆了“分歧拓扑由几何或线性部分决定”的潜在假设。作者证明，即使在线性部分（即分歧点 $\lambda^*$ ）和领头阶模态（ $$Y^*$$ ）完全相同的情况下，仅改变非线性项的奇偶性或凸性，就能反转 $\lambda_2$ 的符号，从而改变分歧类型。这是一个深刻的洞察，表明在预测模式形成稳定性时，必须仔细考虑非线性相互作用的细节。
对质量守恒约束的处理：在封闭流形（如球面）上研究椭圆方程时，全局积分约束是一个普遍存在的技术难点。许多工作要么忽略此约束（在无界域或周期域中可行），要么将其作为事后修正。本文则将此约束置于分析的核心，明确展示了它如何作为不同非线性模型之间的分水岭：奇非线性天然兼容约束，而其他非线性则必须通过激发平均流（zonal mean flow, $$l=0$$ 模式）来适应。这种区分对于理解大气或海洋环流中带状流与非带状流的耦合机制具有直接的物理意义。
适用范围与局限性：本文的框架高度依赖于所选对称群能将核空间降至一维。对于更高阶的球谐（如 $$l=4$$ 对应八面体对称），虽然也可行，但不变子空间的维数可能大于一，导致仍需处理多维分歧方程。此外，分析局限于局部分歧，未涉及大振幅解或全局分支结构。然而，作为对特定对称性下局部分歧机制的精细刻画，本文为更复杂的全局分析奠定了坚实基础。

6. 通过五维提升的轴对称约化构造三维Navier–Stokes方程的有限时间稳定奇性

Stable Finite-Time Singularity Formation for 3D Navier–Stokes via 5D-Lifted Axisymmetric Reductions

作者: Rishad Shahmurov

类别: math.AP

发布日期: 2026-04-10

arXiv ID: 2604.09949v1

摘要翻译: 我们提出了一种五维提升的解析轮廓方案，用于在周期环面 $\T^3$ 上的三维不可压Navier–Stokes方程中构造有限时间奇性。该构造的核心是一个稳态重标度轮廓 $\Ombar$ ，其满足一个在解析加权Hilbert空间 $\Xspace$ 中的非线性椭圆不动点方程，并辅以基于区间算术的计算机辅助Newton–Kantorovich验证。该轮廓被重构为近似自相似的奇性演化，再通过周期延拓与精确Leray投影转移到 $\T^3$ 上。本文采用计算机辅助证明论文的组织形式，包含定理陈述、证明模块，以及关于余项、逆稳定性与Lipschitz界的显式验证常数。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩Navier–Stokes方程在周期环面 $\mathbb{T}^3$ 上是否存在有限时间奇点的问题。该问题自 Leray (1934) 提出以来，是偏微分方程与流体力学中最核心的公开问题之一。经典理论表明：若解在有限时间 $$T^*$$ 内爆破，则必须满足 Beale–Kato–Majda 准则

\int_0^{T^*} \|\nabla u(\cdot, t)\|_{L^\infty} \, dt = \infty,

等价于涡量 $\omega = \nabla \times u$ 的 $L^\infty$ 范数在 $$T^*$$ 处不可积。然而，长期以来尚无严格证明存在满足此条件的光滑初值。

作者采用一种“提升-约化”策略：首先将三维轴对称无旋（swirl-free）情形推广至包含 swirl 分量 $u_\theta$ 的情形，并引入提升变量 $\Gamma = r u_\theta$ 和 swirl 标量 $\Omega = u_\theta / r$ 。在轴对称假设下，三维Navier–Stokes方程可约化为关于 $$(r,z)$$ 的二维系统，但其自然测度并非欧氏，而是带有权重 $$r^3 dr dz$$ 的五维提升几何（5D-lifted geometry）。在此框架下，swirl 方程写作：

\partial_t \Omega + b \cdot \nabla_{r,z} \Omega - \nu B \Omega = -2 \frac{u_r}{r} \Omega,

其中 $$b = (u_r, u_z)$$ 是子午面速度，而算子

B = \partial_{rr} + \frac{3}{r}\partial_r + \partial_{zz}

是对应于五维径向拉普拉斯的 Bessel-caloric 算子。同时，子午速度满足五维不可压条件：

\partial_r u_r + \frac{3}{r} u_r + \partial_z u_z = 0.

关键观察是：若存在一个近自相似（nearly self-similar）爆破解，则其应具有如下形式：

\Omega(r,z,t) = \frac{1}{T^* - t} \bar{\Omega}\left( \rho, \zeta \right), \quad \rho = \frac{r}{\sqrt{T^* - t}}, \quad \zeta = \frac{z}{\sqrt{T^* - t}},

其中 $\bar{\Omega}(\rho, \zeta)$ 是一个稳态剖面（stationary profile），满足一个非线性椭圆型不动点方程。将上述 ansatz 代入演化方程后，得到剖面方程：

G(\bar{\Omega}, \nu) := \bar{\Omega} + \frac{1}{2}(\rho \partial_\rho + \zeta \partial_\zeta)\bar{\Omega} + \bar{b} \cdot \nabla \bar{\Omega} - \nu \left( \partial_{\rho\rho} \bar{\Omega} + \frac{3}{\rho}\partial_\rho \bar{\Omega} + \partial_{\zeta\zeta} \bar{\Omega} \right) + 2 \frac{\bar{u}_\rho}{\rho} \bar{\Omega} = 0.

此处 $\bar{b} = (\bar{u}_\rho, \bar{u}_\zeta)$ 并非自由变量，而是通过一个恢复律（recovery law）由 $\bar{\Omega}$ 隐式确定：引入五维流函数 $\bar{\psi}$ ，定义

B\left( \frac{\bar{\psi}}{\rho^4} \right) = \partial_\zeta (\bar{\Omega}^2),

并令

\bar{u}_\rho = -\frac{1}{\rho^3} \partial_\zeta \bar{\psi}, \quad \bar{u}_\zeta = \frac{1}{\rho^3} \partial_\rho \bar{\psi}.

这一构造保证了五维散度自由条件 $\partial_\rho \bar{u}_\rho + \frac{3}{\rho} \bar{u}_\rho + \partial_\zeta \bar{u}_\zeta = 0$ 成立，且在轴 $\rho=0$ 处无奇异残留（由加权格林恒等式验证）。

本文试图填补的关键空白是：在真实三维Navier–Stokes方程中，是否存在一个数学上严格验证的、稳定的、有限时间爆破的光滑解。此前工作多集中于简化模型（如 Tao 的平均化模型）、Euler 方程（如 Elgindi 的 $C^{1,\alpha}$ 奇点），或数值启发但缺乏严格验证的方案。作者提出一个完整的“解析剖面+计算机辅助证明”程序，在一个精心设计的解析加权 Hilbert 空间 $\mathcal{X}$ 中，通过 Newton–Kantorovich 定理严格验证该非线性椭圆方程存在解，并最终将其转移回 $\mathbb{T}^3$ 上。

2. 主要定理或主要结果

作者在粘性系数 $\nu \in (0, \nu_c)$ 且 $\nu_c \approx 0.00582$ 的条件下，证明了以下主要结论：

定理 2.1（ $\mathbb{T}^3$ 上有限时间奇点的存在性）：存在一个光滑、无散、平方可积的初值 $u_0 \in C^\infty(\mathbb{T}^3)$ ，使得对应的三维不可压缩 Navier–Stokes 方程

\partial_t u + (u \cdot \nabla) u + \nabla p = \nu \Delta u, \quad \nabla \cdot u = 0, \quad x \in \mathbb{T}^3,

的唯一解 $$u(x,t)$$ 在有限时间 $T^* < \infty$ 处发生爆破，即满足 Beale–Kato–Majda 准则：

\int_0^{T^*} \|\nabla u(\cdot, t)\|_{L^\infty} \, dt = \infty.

该爆破是渐近近自相似的，并且在解析加权空间 $\mathcal{X}$ 中是稳定的。

该结论的证明依赖于一系列中间结果：

定理 11.1（精确稳态剖面的存在性）：在解析剖面空间 $\mathcal{X}$ 中，存在唯一的精确解 $\bar{\Omega} \in \mathcal{X}$ 满足剖面方程 $G(\bar{\Omega}, \nu) = 0$ 。该解位于一个数值构造的近似解 $\bar{\Omega}_{\text{app}}$ 的邻域内。
定理 12.1（奇点重构）：由精确剖面 $\bar{\Omega}$ 重构的时间相关场
$\Omega(r, z, t) = \frac{1}{T^* - t} \bar{\Omega}\left( \frac{r}{\sqrt{T^* - t}}, \frac{z}{\sqrt{T^* - t}} \right)$
诱导出三维速度场和涡量场分别为
$u(x, t) = \frac{1}{\sqrt{T^* - t}} \bar{u}\left( \frac{x}{\sqrt{T^* - t}} \right), \quad \omega(x, t) = \frac{1}{T^* - t} \bar{\omega}\left( \frac{x}{\sqrt{T^* - t}} \right).$
因此，涡量的 $L^\infty$ 范数满足 $\|\omega(\cdot, t)\|_{L^\infty} \sim (T^* - t)^{-1} \to \infty$ ，其时间积分发散。
定理 13.2（到周期环面的嵌入）：通过周期延拓和精确的 Leray 投影，将局部提升解转移至 $\mathbb{T}^3$ 上，得到初值 $$u_0$$ 。转移过程的总误差被严格控制在 $\varepsilon_{\mathbb{T}^3} \leq 1.42 \times 10^{-20}$ ，远小于后续分析中的其他误差项，因此不影响爆破性质。

这些结论的成立依赖于一套严格的计算机辅助验证包（validation package），其核心常数如下：

残差界： $\|G(\bar{\Omega}_{\text{app}}, \nu)\|_{\mathcal{X}} \leq \delta = 8.421739 \times 10^{-12}$ 。
线性化算子逆的范数界： $\|DG(\bar{\Omega}_{\text{app}})^{-1}\| \leq M = 482.6$ 。
导数的 Lipschitz 常数： $\|DG(u) - DG(v)\|_{\mathcal{L}(\mathcal{Y}, \mathcal{X})} \leq K \|u - v\|_{\mathcal{Y}}$ ，其中 $K = 1.1 \times 10^4$ 。

这些常数共同满足 Newton–Kantorovich 条件： $2\delta M K \approx 8.9 \times 10^{-5} < 1$ ，从而保证了精确解的存在性和唯一性。

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

本文的核心技术障碍在于：Navier–Stokes 方程的非线性项（对流项和涡量拉伸项）与耗散项之间存在精细的平衡，任何微小的扰动都可能导致爆破场景崩溃。传统的能量方法或先验估计无法处理这种临界平衡。作者的真正创新点在于构建了一个封闭的、可计算的、解析的框架，将无限维的 PDE 问题转化为一个可在计算机上严格验证的非线性泛函方程问题。具体技术工具如下：

1. 五维提升几何与解析加权空间 $\mathcal{X}$ ：作者没有在标准 Sobolev 或 Hölder 空间中工作，而是定义了一个基于紧化提升算子（compactified lifted operator）本征函数 $\{\Phi_j\}$ 的解析 Hilbert 空间 $\mathcal{X}$ ，其范数为

\|\Omega\|_{\mathcal{X}}^2 := \sum_{j=1}^\infty (1 + j^2)^6 e^{2\tau j} |\hat{\Omega}_j|^2, \quad \tau = 0.08.

该空间的选择至关重要：指数权重 $e^{2\tau j}$ 强制剖面具有 Gevrey 正则性（强于解析性），确保了所有非线性操作（如乘积、求导）的封闭性；多项式权重 $$(1+j^2)^6$$ 则用于吸收恢复律带来的高阶导数损失。同时，作者还定义了一个更强的源空间 $\mathcal{Y}$ （ $\tau' = 0.081$ ，权重指数为7），为 Newton–Kantorovich 定理提供所需的范数对 $(\mathcal{Y}, \mathcal{X})$ 。

2. 动力学恢复律及其惩罚因子 $j^{7/2}$ ：从 $\bar{\Omega}$ 恢复速度场 $\bar{b}$ 是证明中最脆弱的环节。引理 6.1 严格量化了这一过程的“惩罚”：拉伸项 $\|u_r/r\|_{L^\infty}$ 的界依赖于谱系数的加权和 $\sum j^{7/2} E_j^{1/2}$ 。这个 $j^{7/2}$ 因子源于两部分：椭圆反演（ $B(\bar{\psi}/\rho^4) = \partial_\zeta(\bar{\Omega}^2)$ ）带来 $j^{-2}$ 的增益，而加权导数（ $\bar{u}_\rho = -\rho^{-3}\partial_\zeta \bar{\psi}$ ）带来 $j^{11/2}$ 的损失，净效应为 $j^{7/2}$ 。这个巨大的惩罚因子直接决定了 Lipschitz 常数 $$K$$ 的大小，并成为整个验证方案成败的关键。

3. 计算机辅助的 Newton–Kantorovich 验证包：这是本文方法论的核心。作者将问题分解为三个可独立验证的“包”：

残差包：通过 Arb 库的区间算术，对截断至 $$N=450$$ 模的近似解 $\bar{\Omega}_{\text{app}}$ 计算残差 $R = G(\bar{\Omega}_{\text{app}}, \nu)$ 。由于非线性项的卷积最多影响到第900模，尾部（ $$j>450$$ ）可通过解析方法严格界定，总残差 $\delta$ 极小。
Jacobian 逆包：对有限维截断 Jacobian $$J_N$$ ，通过构造近似逆 $$R$$ 并验证 $\|I - RJ_N\| < 1$ ，利用 Neumann 级数保证其可逆性，并给出逆的范数界 $$M$$ 。对于无限维尾部（ $$j>450$$ ），通过引理 9.1 的“尾部强制不等式”证明其正定性，其下界 $\gamma \geq 7182.4$ 远大于由剖面引起的扰动项 $\text{Inter}_j$ 。
Lipschitz 包：通过引理 10.1 和 10.2，分别控制恢复算子的映射常数 $C_{\text{rec}}^{\text{map}} \leq 2.5652 \times 10^7$ 和非线性相互作用的卷积常数 $C_{\text{conv}} \leq 4.2872 \times 10^{-4}$ ，从而得到整体 Lipschitz 常数 $K \approx C_{\text{rec}}^{\text{map}} \cdot C_{\text{conv}}$ 。

这三个包的数据最终汇入 Corollary 10.4，验证了 Newton–Kantorovich 条件 $2\delta M K < 1$ ，从而完成了从数值近似到严格数学证明的飞跃。

4. 到周期环面的严格转移：最后一步是将局部定义的提升解无缝嵌入到全局的 $\mathbb{T}^3$ 上。作者通过 Poincaré 和（周期延拓）加上精确的 Leray 投影来实现。引理 13.1 证明，由于剖面是高斯局域化的（尺度 $\sigma=0.05$ ），非中心周期像之间的重叠误差、投影修正误差和压力畸变误差均为指数小量，总和不超过 $10^{-20}$ ，完全可以忽略不计。

4. 与之前工作的比较

本文与已有工作的比较可从多个维度展开：

1. 与简化模型和 Euler 方程工作的比较： Elgindi (2019) 在轴对称 Euler 方程中证明了 $C^{1,\alpha}$ 初值下的有限时间奇点，但其结果依赖于 Euler 方程无耗散的特性，且初值正则性低于 $$C^2$$ 。Hou 和 Wang 等人对一维模型的分析揭示了输运与拉伸的竞争机制，但模型本身已大幅偏离真实 Navier–Stokes 方程。Tao (2016) 的平均化模型虽然保留了 Navier–Stokes 的尺度不变性，但修改了非线性项的结构。相比之下，本文处理的是原汁原味的三维不可压缩 Navier–Stokes 方程，且初值是 $C^\infty$ 光滑的，这是一个质的飞跃。

2. 与数值模拟和启发式工作的比较：过去几十年有大量高分辨率数值模拟声称观测到潜在的自相似爆破（如 Kerr, Hou 等人的工作），但这些结果始终面临网格分辨率不足、耗散效应干扰等质疑，无法构成数学证明。本文的方案虽然也始于数值构造（ $\bar{\Omega}_{\text{app}}$ ），但通过计算机辅助证明（CAP）框架，将数值结果提升为严格数学定理。CAP 方法近年来在动力系统和 PDE 领域取得成功（如文中引用的 [11,12,3,4]），但将其应用于 Navier–Stokes 奇点问题尚属首次，且本文构建的解析框架（五维提升、Gevrey 空间、恢复律）是专门为该问题定制的。

3. 与 Jia–Šverák 自相似框架的比较： Jia 和 Šverák (2013, 2015) 发展了一套关于 Navier–Stokes 方程自相似解的深刻理论，但他们主要关注全局存在性和稳定性，并未构造出爆破解。他们的框架通常假设剖面属于较弱的空间（如 $$L^p$$ 或能量空间）。本文则反其道而行之，主动寻找一个满足特定边界条件（奇性、轴对称、 $\zeta$ -奇）的剖面，并在一个非常强的解析空间 $\mathcal{X}$ 中进行验证。这种强空间的选择牺牲了普适性，但换来了非线性项的良好控制和 CAP 方法的可行性。

4. 关于粘性 $\nu$ 的范围和代价：本文的结果仅在小粘性 $\nu < \nu_c \approx 0.00582$ 下成立。这是一个明显的局限性，因为物理上更关心的是任意固定 $\nu > 0$ （特别是 $\nu=1$ ）的情形。然而，这几乎是所有试图构造爆破解的工作都面临的困境：粘性越大，耗散越强，越难维持爆破。作者通过 CAP 方法将 $\nu_c$ 显式计算出来，这本身就是一个进步。代价是，该方法目前无法处理大粘性情形，且对 $\nu$ 的依赖是通过线性化算子的强制性（ $\nu j^2$ ）体现的。

5. 时间尺度与适用范围：本文构造的解是近自相似的，而非严格的自相似解。这意味着爆破率可能有对数修正，但其主导行为仍是 $(T^*-t)^{-1}$ 。此外，该解是稳定的，即在 $\mathcal{X}$ 范数的小邻域内，所有初值都会导致爆破。这比仅仅证明存在一个孤立的爆破解要强得多。然而，这种稳定性是在一个非常特殊的函数空间中定义的，其物理意义和在更广泛空间（如 $$H^1$$ ）中的稳定性仍有待研究。

综上所述，本文并非宣称解决了 Navier–Stokes 的全局正则性猜想，而是通过一个高度创新的、结合了深刻解析洞察和前沿计算验证的方法，在一个受控的、小粘性的、强正则性的设定下，首次为三维 Navier–Stokes 方程的有限时间奇点提供了严格的数学证据。

7. 不可压缩Navier-Stokes方程含未知体力的初值反问题的时间维约化与Carleman压缩原理

Time-dimensional reduction and a Carleman contraction principle for an inverse initial data problem for the incompressible Navier-Stokes equations with unknown body force

作者: Phuong M. Nguyen, Loc H. Nguyen

类别: math.AP

发布日期: 2026-04-10

arXiv ID: 2604.09934v1

摘要翻译: 我们研究不可压缩Navier-Stokes方程组的一个初值反问题，目标是从侧边界的观测数据中重构初始速度与压力，且不假设时间无关的体力项已知。为消除该未知体力项，我们对动量方程关于时间求导，随后应用Legendre多项式-指数型时间维约化方法。该过程导出一组关于展开系数的耦合椭圆方程组。接着，我们在一个配备Carleman加权范数的适当容许集上，为此约化系统构造一个压缩映射。该映射的不动点给出时间维约化模型的一个近似解，而压缩性保证了Picard迭代的全局收敛性。最后，我们基于此框架提出一种数值算法，并通过数值实验展示了从合成边界数据中对初始速度和压力的高精度重构。

Comment: 无

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

1. 研究对象和背景

本文研究如下不可压Navier–Stokes系统的初值反问题：设 $\Omega \subset \mathbb{R}^d$ （ $d \geq 2$ ）为有界光滑区域， $$T > 0$$ 为给定时间。考虑速度场 $u: \Omega \times (0,T) \to \mathbb{R}^d$ 和压力场 $p: \Omega \times (0,T) \to \mathbb{R}$ 满足

\begin{cases} u_t + (u \cdot \nabla) u = -\nabla p + \mu \Delta u + f & \text{in } \Omega \times (0,T), \\ \operatorname{div} u = 0 & \text{in } \Omega \times (0,T), \\ u = 0 & \text{on } \partial\Omega \times (0,T), \\ u(x,0) = u_0(x) & \text{in } \Omega, \end{cases} \tag{1.1}

其中 $\mu > 0$ 为粘性系数， $f: \Omega \to \mathbb{R}^d$ 是时间无关但未知的体力项。目标是从侧边界的观测数据

g(x,t) := \partial_\nu u(x,t), \quad h(x,t) := \begin{bmatrix} p(x,t) \\ \partial_\nu p(x,t) \end{bmatrix}, \quad (x,t) \in \partial\Omega \times (0,T) \tag{1.2}

中重构初始速度 $u_0 = u(\cdot,0)$ 和初始压力 $p_0 = p(\cdot,0)$ 。

该问题在流体力学应用中具有重要意义：在航空航天、海洋环流或环境流体模拟中，初始状态常无法直接测量，但高分辨率边界传感器可提供近壁面的速度梯度和压力信息。然而，标准反问题理论通常假设模型中除目标量外其余参数（如 $$f$$ ）已知，而本文允许 $$f$$ 完全未知，这更贴近实际但显著增加数学难度。

从分析角度看，问题的核心困难在于：(i) 控制方程是非线性的；(ii) 观测仅限于边界；(iii) 未知量包括 $$u_0, p_0$$ 和 $$f$$ ，但方程仅提供 $$d+1$$ 个分量（动量方程 $$d$$ 个 + 不可压条件 1 个），自由度严重不足。作者的关键洞察是利用 $$f$$ 的时间无关性——对动量方程关于时间求导可消去 $$f$$ ，得到

u_{tt} + (u_t \cdot \nabla) u + (u \cdot \nabla) u_t = -\nabla p_t + \mu \Delta u_t \quad \text{in } \Omega \times (0,T). \tag{3.1}

再结合不可压条件 $\operatorname{div} u = 0$ ，可进一步导出压力时间导数的椭圆方程

\Delta p_t = -\nabla u_t : (\nabla u)^\top - \nabla u : (\nabla u_t)^\top \quad \text{in } \Omega \times (0,T). \tag{3.2}

至此，原抛物型系统被转化为一个不含 $$f$$ 的二阶时间导数系统，但仍是时变的偏微分方程组。

为将时间变量离散化，作者引入Legendre多项式-指数基函数：定义标准Legendre多项式 $$P_n$$ 在 $$(0,T)$$ 上的平移归一化版本

Q_n(t) = \sqrt{\frac{2n+1}{T}} P_n\left( \frac{2t}{T} - 1 \right),

并构造加权基函数 $\Psi_n(t) = e^t Q_n(t)$ 。该基在加权内积 $\langle u,v \rangle_{e^{-2t}} = \int_0^T e^{-2t} u(t)v(t) dt$ 下正交。关键性质是：时间导数算子在此基下的表示矩阵 $S_N = [s_{mn}]_{m,n=0}^N$ （其中 $s_{mn} = \int_0^T e^{-2t} \Psi_n'(t) \Psi_m(t) dt$ ）是上三角且对角元全为1（见Proposition 1），故可逆。

将 $u(x,t) = \sum_{n=0}^\infty u_n(x) \Psi_n(t)$ 、 $p(x,t) = \sum_{n=0}^\infty p_n(x) \Psi_n(t)$ 代入(3.1)-(3.2)，截断至前 $$N+1$$ 项，并投影到基函数上，得到关于系数向量 $U^N = (u_0,\dots,u_N)^\top$ 、 $P^N = (p_0,\dots,p_N)^\top$ 的耦合椭圆系统：

\begin{cases} \mu S_N (\Delta U^N) = R_N U^N + F^N(x, U^N, \nabla U^N) + S_N (\nabla P^N) & \text{in } \Omega, \\ S_N (\Delta P^N) = G^N(x, U^N, \nabla U^N) & \text{in } \Omega, \end{cases} \tag{3.13a}

其中 $R_N = [r_{mn}]$ 来自二阶时间导数项，非线性项 $$F^N, G^N$$ 由三重积分系数 $c_{mkn} = \int_0^T e^{-2t} \Psi_n'(t) \Psi_k(t) \Psi_m(t) dt$ 构成。边界条件由观测数据投影得到：

U^N|_{\partial\Omega} = 0, \quad \partial_\nu U^N|_{\partial\Omega} = \Gamma^N, \quad P^N|_{\partial\Omega} = H_1^N, \quad \partial_\nu P^N|_{\partial\Omega} = H_2^N. \tag{3.13b}

此约化模型将原时间依赖的反问题转化为有限维椭圆系统的Cauchy问题。由于Cauchy问题不适定，需借助Carleman估计构建稳定解法。本文试图填补的空白是：在未知体力项条件下，建立一种不依赖初值猜测、全局收敛的重构框架，克服传统优化方法因非凸性导致的局部极小问题。

2. 主要定理或主要结果

作者在以下假设下证明了主要结果：

正则性假设：真实解满足 $u \in H^k((0,T); H^s(\Omega))$ 、 $p \in H^k((0,T); H^s(\Omega))$ ，其中整数 $k \geq 5$ ， $$s > d/2 + 2$$ （保证嵌入 $H^s(\Omega) \hookrightarrow C^2(\bar{\Omega})$ ）。
几何假设：存在点 $x_0 \notin \bar{\Omega}$ 使得 $$r(x) = |x - x_0| > 1$$ 对所有 $x \in \Omega$ 成立（用于Carleman权重 $\mu_\beta(x) = r^{-\beta}(x)$ ）。
适定性假设：存在足够大的常数 $$M > 0$$ ，使得容许集 $\mathcal{A} = \left\{ (V,Q) \in [H^s(\Omega)]^{d(N+1)} \times [H^s(\Omega)]^{N+1} : \|V\|_{[L^\infty]^d} \leq M, \|Q\|_{L^\infty} \leq M, \text{ 并满足边界条件} \right\}$ 非空。

定理4.1（压缩映射的存在唯一性）：存在常数 $\lambda_0 > 0$ （依赖于 $x_0, \Omega, d, \beta, N, \mu, M$ ），使得对任意 Carleman 参数 $\lambda \geq \lambda_0$ ，定义在 $\mathcal{A}$ 上的算子 $\mathcal{T}_\lambda: \mathcal{A} \to \mathcal{A}$ 是压缩映射。具体而言，对任意 $(V_1, Q_1), (V_2, Q_2) \in \mathcal{A}$ ，有

\|\mathcal{T}_\lambda(V_1, Q_1) - \mathcal{T}_\lambda(V_2, Q_2)\|_{\mathcal{H}_\lambda} \leq \theta \|(V_1, Q_1) - (V_2, Q_2)\|_{\mathcal{H}_\lambda},

其中 $\theta \in (0,1)$ ， $\mathcal{H}_\lambda$ 是配备 Carleman 加权范数

\|(V,Q)\|_{\mathcal{H}_\lambda}^2 = \sum_{n=0}^N \int_\Omega e^{2\lambda \mu_\beta(x)} \left( |\Delta v_n|^2 + |\Delta q_n|^2 \right) dx

的 Hilbert 空间。

推论4.2（重构误差估计）：设 $$(U^N, P^N)$$ 为约化系统(3.13)的真实解（若存在）， $(\tilde{U}^N, \tilde{P}^N)$ 为 $\mathcal{T}_\lambda$ 的不动点，则存在常数 $$C > 0$$ 使得

\|(U^N - \tilde{U}^N, P^N - \tilde{P}^N)\|_{[H^1(\Omega)]^{d(N+1)} \times [H^1(\Omega)]^{N+1}} \leq C \lambda^{-1/2} \delta,

其中 $\delta$ 表征观测数据误差（即 $\Gamma^N, H_1^N, H_2^N$ 的扰动大小）。这表明误差随 Carleman 参数 $\lambda$ 增大而衰减。

数值算法：基于上述框架，作者提出 Picard 迭代算法：从任意初值 $(U^{(0)}, P^{(0)}) \in \mathcal{A}$ 出发，迭代求解线性椭圆系统

\begin{cases} \mu S_N \Delta U^{(k+1)} - S_N \nabla P^{(k+1)} = R_N U^{(k)} + F^N(x, U^{(k)}, \nabla U^{(k)}), \\ S_N \Delta P^{(k+1)} = G^N(x, U^{(k)}, \nabla U^{(k)}), \end{cases}

配合 Carleman 加权最小二乘正则化。数值实验显示，即使初值远离真解，迭代仍能准确重构 $$u_0$$ 和 $$p_0$$ 。

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

本文的核心技术障碍在于：约化后的椭圆系统(3.13)包含强非线性项 $$F^N, G^N$$ ，其结构复杂（涉及速度梯度的双线性耦合），无法直接套用现有Carleman压缩框架（如[24,33]中要求非线性项满足Lipschitz条件且增长阶较低）。作者的创新点在于构建了一个适配Navier–Stokes结构的Carleman压缩映射，具体分三步：

1. 力消除与时间维约化的协同设计：通过时间微分消去未知体力 $$f$$ 是关键预处理。但直接对(3.1)-(3.2)进行时间离散会引入高阶导数项。作者采用Legendre多项式-指数基 $\Psi_n(t) = e^t Q_n(t)$ ，其优势在于：(a) 时间导数矩阵 $$S_N$$ 上三角且可逆（Proposition 1），便于显式解耦；(b) 加权正交性使投影操作保持能量结构。特别地， $$S_N$$ 的对角占优性（对角元为1，上三角部分由低阶多项式正交性决定）确保了离散系统的良态性。

2. 非线性项的精细估计：非线性项 $F^N_m = \sum_{k,n} c_{mkn} [(u_n \cdot \nabla) u_k + (u_k \cdot \nabla) u_n]$ 和 $G^N_m = -\sum_{k,n} c_{mkn} [\nabla u_n : (\nabla u_k)^\top + \nabla u_k : (\nabla u_n)^\top]$ 的处理是难点。作者利用以下事实：(i) 系数 $c_{mkn}$ 由三重积分定义，但因 $\Psi_n'$ 是次数 $\leq n$ 的多项式乘以 $$e^t$$ ，结合 $$Q_n$$ 的正交性，可证 $c_{mkn}$ 在 $$m > k+n+1$$ 时为零（虽未显式写出，但隐含在截断合理性中）；(ii) 在容许集 $\mathcal{A}$ 中， $\|U^N\|_{L^\infty} \leq M$ ，故非线性项满足全局Lipschitz条件：

\|F^N(U_1) - F^N(U_2)\|_{L^2} \leq C(M,N) \|U_1 - U_2\|_{H^1},

类似对 $$G^N$$ 成立。这为后续压缩性证明奠定基础。

3. Carleman加权范数下的压缩性证明：定义算子 $\mathcal{T}_\lambda$ 为冻结非线性项后的线性系统解算子。对任意 $(V_1,Q_1), (V_2,Q_2) \in \mathcal{A}$ ，令 $(W,R) = \mathcal{T}_\lambda(V_1,Q_1) - \mathcal{T}_\lambda(V_2,Q_2)$ ，则 $$(W,R)$$ 满足齐次边界条件的线性系统：

\begin{cases} \mu S_N \Delta W - S_N \nabla R = \mathcal{F}(V_1,V_2), \\ S_N \Delta R = \mathcal{G}(V_1,V_2), \end{cases}

其中 $\mathcal{F}, \mathcal{G}$ 为非线性项差。应用Lemma 1的Carleman估计到每个分量 $$w_n, r_n$$ （注意 $$S_N$$ 可逆，可解出 $\Delta w_n, \Delta r_n$ ），得到

\int_\Omega e^{2\lambda \mu_\beta} (|\Delta W|^2 + |\Delta R|^2) dx \geq C \lambda^3 \int_\Omega e^{2\lambda \mu_\beta} (|W|^2 + |R|^2) dx + C \lambda \int_\Omega e^{2\lambda \mu_\beta} (|\nabla W|^2 + |\nabla R|^2) dx.

同时，右边非线性项差满足

\|\mathcal{F}\|_{L^2} + \|\mathcal{G}\|_{L^2} \leq C \|(V_1 - V_2, Q_1 - Q_2)\|_{H^1}.

结合Carleman权重 $e^{2\lambda \mu_\beta}$ 的指数衰减特性（在远离 $$x_0$$ 处主导），通过选择足够大的 $\lambda$ ，可使非线性项贡献被Carleman主部控制，从而导出压缩因子 $\theta < 1$ 。此步骤的关键是Carleman估计中 $\lambda^3$ 阶项对低阶项的压制能力。

4. 一致性分析：最后证明当正则化参数（隐含在 $\lambda$ 中）趋于零时，不动点收敛到约化系统的真实解。这依赖于先验正则性假设和紧性论证。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开：

1. 未知体力项的处理：绝大多数Navier–Stokes反问题文献（如[8,12,20,30,50]）均假设体力 $$f$$ 已知。例如，[48]虽研究初值反问题，但明确要求已知 $$f$$ 和 $$p$$ 。本文首次在 $$f$$ 完全未知条件下重构 $$(u_0, p_0)$$ ，通过时间微分消去 $$f$$ 是本质创新。代价是方程阶数提升（需二阶时间导数），但通过时间维约化有效化解。

2. 方法论对比：优化 vs. 压缩映射：传统方法（如[8,50]）将反问题表述为输出最小二乘优化，依赖伴随梯度下降。但Navier–Stokes的强非线性导致目标泛函高度非凸，需良好初值（实践中难获）。本文的Carleman-Picard方法不依赖初值，通过构造压缩映射保证全局收敛（Theorem 4.1）。这是对[24,33]框架的实质性拓展——此前该框架仅用于拟线性抛物方程，而本文成功处理了Navier–Stokes特有的对流非线性和不可压约束耦合。

3. 时间处理技术：时间维约化思想源于[45]，后用于反问题[28,48]。但[48]的阻尼Picard迭代未证明收敛性，且需已知 $$f$$ 。本文的Legendre-指数基选择更优： $$S_N$$ 的上三角结构（Proposition 1）比[48]中使用的其他基更利于分析。此外，[24,33]的约化针对单一方程，而本文需同时处理速度-压力耦合系统，非线性项 $$F^N, G^N$$ 的交叉结构带来新挑战。

4. Carleman估计的应用深度：Carleman估计在反问题中常见（如[3,4,16,32]），但多用于唯一性或稳定性证明。本文将其嵌入数值算法设计：Carleman加权范数直接用于定义压缩映射，使理论收敛性（Theorem 4.1）与数值实现无缝衔接。相比之下，[12]用复几何光学解证唯一性，[20,30]用边界观测证粘性唯一性，均未提供重构算法。本文的误差界（推论4.2）明确量化了Carleman参数 $\lambda$ 与重构精度的关系，这是纯理论工作不具备的。

中文标题: 一个双重非线性抛物方程的博弈论解释

摘要翻译: 我们为双重非线性抛物方程

中文标题: 具有退化势的长程相共存模型

摘要翻译: 本综述概述了非局部相变问题的近期进展，此类问题由Ginzburg–Landau型能量泛函建模，其形式为

中文标题: 流体力学中非局部拟线性演化方程的Schauder型估计与适定性

摘要翻译: 我们建立了由变系数非局部伪微分算子（阶数 $$s>0$$ ）驱动的线性抛物系统的Schauder型估计。这些估计在临界时间加权的Hölder/Besov型空间中表述，并适用于处于尺度临界正则性的拟线性方程。

中文标题: Robin 边界条件下椭圆问题解的渐近行为

摘要翻译: 本文研究半线性椭圆 Robin 问题

中文标题: 具有局部—非局部耦合算子的抛物—椭圆动力学

摘要翻译: 本文研究两类局部—非局部型抛物—椭圆演化系统。在所考虑的问题中，空间区域 $$Ω$$ 被划分为互不相交的子区域 $$A$$ 与 $$B$$ ，即 $Ω = A \cup B$ 。首先，我们在子区域 $$A$$ 上考虑一个局部抛物方程，并在补集 $$B$$ 上配以一个非局部椭圆平衡方程；随后，我们交换二者角色，在 $$A$$ 上设置局部椭圆方程，并与定义在 $$B$$ 上的非局部抛物方程耦合。在这两种模型中，两个区域之间的相互作用由一个非局部传输项驱动，该传输项通过一个核函数实现界面间的质量转移，从而产生混合的局部—非局部、椭圆—抛物动力学。两个问题均配备 Neumann 边界条件。为开展分析，我们首先利用不动点方法建立解的存在唯一性。接着，我们对解的定性行为进行详细分析。特别地，我们证明尽管系统具有椭圆—抛物结构，其耦合机制仍诱导出一个自然的能量泛函，系统的演化由该能量泛函的梯度流所支配。正如 Neumann 边界条件下所预期的那样，我们证明整个区域 $$Ω$$ 中的总质量在时间演化中保持守恒。我们还分析了长时间行为，并获得了抛物分量的衰减估计，这些估计进一步驱动椭圆部分收敛到一个常数解。最后，我们证明所研究的抛物—椭圆问题可视为一个纯抛物问题在某一控制某组分演化速度的参数趋于零时的极限情形。

不相关论文（基本信息）

1. 平均场量子电动力学中被俘获玻色子、非线性共振级联与动力学玻色-爱因斯坦凝聚的形成

Trapped bosons in mean field QED, nonlinear resonance cascades and dynamical BEC formation

作者: Thomas Chen, Ali Mezher

类别: math-ph, math.AP, quant-ph

发布日期: 2026-04-13

arXiv ID: 2604.11756v1

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

摘要: In this paper, we study a system of bosons trapped in a confining potential, interacting with a quantized field of coherent photons in the mean field description of non-relativistic Quantum Electrodynamics (QED) obtained by [N. Leopold and P. Pickl , 2017]. We derive the effective nonlinear cascade equations governing the emission and absorption of coherent photons by the boson subsystem in a combined weak-coupling and kinetic-scaling limit. We demonstrate that solutions to this nonlinear cascade describe a monotone decreasing energy flow in the boson subsystem. Thereby, we prove that a Bose-Einstein condensate (BEC) forms dynamically, under conservation of the total boson $$L^2$$ mass. We note that this process is crucially different from thermal relaxation to the ground state, and fundamentally depends on the nonlinear nature of the cascade dynamics.

2. 数学中的人工智能讲座

Lectures on AI for Mathematics

作者: Xiaoyang Chen, Xiaoyang Chen

类别: cs.AI, math.AP, math.AT, math.DG

发布日期: 2026-04-13

arXiv ID: 2604.11504v1

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

摘要: This book provides a comprehensive and accessible introduction to the emerging field of AI for mathematics. It covers the core principles and diverse applications of using artificial intelligence to advance mathematical research. Through clear explanations, the text explores how AI can discover hidden mathematical patterns, assist in proving complicated theorems, and even construct counterexamples to challenge conjectures.

3. 细胞电穿孔的随机扰动模型

Stochastically perturbed model of cell electropermeabilization

作者: Tobias Gebäck, Oleksandr Misiats, Ioanna Motschan-Armen, Irina Pettersson

类别: math.AP

发布日期: 2026-04-13

arXiv ID: 2604.11260v1

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

摘要: Reversible electropermeabilization, commonly referred to as electroporation, is a transient increase in cell membrane permeability induced by short, high-voltage electric pulses. We present a stochastically perturbed version of a phenomenological electroporation model introduced in the deterministic setting by \cite{kavian2014classical}. The deterministic model couples the electrostatic equations for the electric potential in the extra- and intracellular domains with a nonlinear evolution law for the transmembrane potential jump, itself coupled to an ordinary differential equation describing the porosity degree of the membrane. To account for various random effects, such as temperature fluctuations or uncerntainty in the applied electric field, we add noise on the cell membrane. We establish the existence and uniqueness of a variational solution to the resulting coupled SPDE-ODE system governing the membrane potential and the degree of porosity, where the stochastic perturbation is multiplicative and degenerate, acting only on the SPDE component of the coupled SPDE-ODE system. Any mixing in the ODE variables is therefore induced indirectly through the nonlinear coupling in the drift. The main technical challenge arises from the nonlinearities, which are neither Lipschitz continuous nor monotone. The result is proved by means of Galerkin method, following the methodology by Liu and Röckner \cite{liu2015stochastic} for treating equations under generalized monotonicity and coercivity conditions. Finally, we present numerical simulations of the solution and its time averages for both additive and multiplicative noise, that provide a numerical indication for existence of invariant measure.

4. 抛物方程在稀疏移动观测下的反源问题

Inverse source problem for the parabolic equation with sparse moving observations

作者: Qiling Gu, Wenlong Zhang, Zhidong Zhang

类别: math.AP, math.NA

发布日期: 2026-04-13

arXiv ID: 2604.11157v1

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

摘要: This paper considers the inverse problem of identifying the source term of parabolic equations from sparse boundary measurements. We used data from moving sensors to locate the unknown source term. This work first proves the uniqueness of the inverse problem under such measurements. Then the movement strategy of the sensor is given, from which the authors build the reconstruction algorithm. Finally, some numerical experiments are performed and the corresponding results are generated, which indicate the effectiveness of the algorithms.

5. 完备流形中余维一等距浸入的刚性

Rigidity of codimension-1 isometric immersions in complete manifolds

作者: Mert Baştuğ

类别: math.AP, math.DG

发布日期: 2026-04-13

arXiv ID: 2604.11130v1

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

摘要: We establish an asymptotic rigidity result for isometric immersions of codimension-1. Specifically, we consider a sequence of immersions from a compact $$d$$ -dimensional manifold into a complete $$(d+1)$$ -dimensional manifold whose elastic energies vanish asymptotically, where the elastic energy quantifies both stretching and bending. We show that such a sequence admits a subsequence converging to an isometric immersion. This extends a result of Alpern, Kupferman, and Maor to the case of complete target manifolds, where the lack of compactness introduces additional analytical difficulties. The proof is based on an approach using local quantitative rigidity estimates, obtained via a reduction to the Euclidean setting. This method avoids the use of Young measures and provides a flexible framework that may be of independent interest.

6. 一类Hessian商方程的 $$L_p$$ 对偶Christoffel-Minkowski型问题

The $$L_p$$ dual Christoffel-Minkowski type problem for a class of Hessian quotient equations

作者: Shasha Luo, Jiabao Gong, Qiang Tu

类别: math.AP

发布日期: 2026-04-13

arXiv ID: 2604.10924v1

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

摘要: In this paper, we investigate an $$L_p$$ dual Christoffel-Minkowski type problem for the Hessian quotient operator $\frac{σ_{k}(Λ)}{σ_{l}(Λ)}$ , where the operator $$Λ$$ has been widely studied in the literature. Exploiting the recently discovered ``inverse convexity’’ property of this class of operators, we establish a full rank theorem under suitable structural assumptions. Together with a priori estimates, this result enables us to prove the existence and uniqueness of strictly spherically convex solutions to the above $$L_p$$ dual Christoffel-Minkowski type problem.

7. 密度测度及其应用

Density measures and applications

作者: Friedemann Schuricht

类别: math.AP, math.OC

发布日期: 2026-04-12

arXiv ID: 2604.10670v1

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

摘要: The paper, that continuous some previous work of Schönherr & Schuricht, treats density measures on ${\mathbb R}^n$ that concentrate in any neighborhood of a Lebesgue null set. Such measures are typical for purely finitely additive measures. We study their basic properties and investigate related integrals. Measures taking only the values 0 and 1 are considered as special case. The results are first applied to weak convergence in $\mathcal{L}^\infty(Ω)$ . Then we derive integral representations by means of such measures for several notions of differentiability for integrable functions and we show a kind of mean value theorem for some class of Sobolev functions. Finally we provide a new approach to the generalized Jacobians in the sense of Clarke.

8. 对数加权求积域的分析

Analysis of Log-Weighted Quadrature Domains

作者: Andrew Graven

类别: math.CV, math-ph, math.AP, math.NA

发布日期: 2026-04-12

arXiv ID: 2604.10394v1

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

摘要: This paper studies plane domains satisfying a quadrature identity with respect to the singular weight $ρ_0(w)=|w|^{-2}$ . These are referred to as log-weighted quadrature domains (LQDs). The logarithmic singularity at $$w=0$$ leads to phenomena not present in the classical theory: in particular, when the domain contains the origin, the associated quadrature data are no longer unique, but are determined only up to a point charge at $$0$$ . A generalized Schwarz function characterization of LQDs is established together with a natural formulation of the inverse problem in the singular setting. In the simply connected case, it is shown that a domain is an LQD if and only if the outer factor of its Riemann map extends to the exponential of a rational function. This characterization yields explicit formulae relating the quadrature function and the Riemann map via the Faber transform, thereby extending earlier formulae from the non-singular theory. Several basic classes of LQDs are also covered, and explicit examples are computed.

9. 具有任意泊松比的非胡克弹性

Non-Hookean elasticity with arbitrary Poisson’s ratios

作者: Mikhail Itskov

类别: math.AP, cond-mat.mtrl-sci, math-ph

发布日期: 2026-04-11

arXiv ID: 2604.10153v1

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

摘要: In a previous paper \cite{Itskov-MoSM} we presented a hyperelastic isotropic material model whose stress-strain response is non-linear even at infinitesimal deformations and cannot thus be linearized. As a result values of Poisson’s ratio greater than one half were obtained. In this contribution, we further propose an isotropic strain energy function which is always positive-definite and depending on material constants delivers arbitrary values of Poisson’s ratio (except of $$-1$$ ) in agreement with the laws of thermodynamics. The model response appears plausible in various deformation states.

10. 关于一类 $$N$$ -刘维尔方程在 $\mathbb{R}^N$ 中解的分类

On the classification of solutions to a class of $$N$$ -Liouville equations in $\mathbb{R}^N$

作者: Giulio Ciraolo, Pierpaolo Esposito, Xiaoliang Li

类别: math.AP

发布日期: 2026-04-11

arXiv ID: 2604.10050v1

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

摘要: Given $N\geq 2$ and $$α>-1$$ , we consider the following weighted Liouville-type equation involving the $$N$$ -Laplacian: \begin{equation*} \left{ \begin{aligned} -& Δ_N u = |x|^{Nα} e^u \quad \text{ in } \mathbb{R}^N && , \ & \int_{\mathbb{R}^N} |x|^{Nα} e^u , dx < + \infty,. &&\end{aligned} \right. \end{equation*} Solutions have been completely classified when $$N=2$$ via complex analysis, and when $$α=0$$ using Pohozaev identities and an isoperimetric argument. In this paper, we first devise a $$P$$ -function approach to the classification result for all $$α>-1$$ when $$N=2$$ . Since it is not based on complex analysis, this alternative and more PDE-oriented approach naturally extends to $N\geq 3$ by providing the classification for any $-1<α\leq 0$ . In particular, the explicit radial solutions are the unique ones for $-1<α\leq0$ but become degenerate for special values $$α_k>0$$ , a hint that non-radial solutions might arise for $$α>0$$ as it happens when $$N=2$$ .

11. 关于（超）Ricci 流的弱形式

On weak formulations of (super) Ricci flows

作者: Sajjad Lakzian

类别: math.DG, math.AP, math.FA, math.MG

发布日期: 2026-04-11

arXiv ID: 2604.10007v1

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

摘要: We present two characterizations of smooth compact Ricci flow solutions solely in terms of metrics and measures (one of them only works under positive scalar curvature along the flow); thus, provide weak formulations that are generalized to the singular setting in a straightforward manner. These formulations are achieved by weakly formulating super Ricci flows and imposing a saturation condition (solely in terms of metric and measure) to ensure the super Ricci flow inequality is an equality.

12. 关于具有夹层型与任意增长的Kirchhoff型p(.)-Laplace问题

On Kirchhoff-type p(.)-Laplacian problems withsandwich-type and arbitrary growth

作者: Ky Ho

类别: math.AP

发布日期: 2026-04-10

arXiv ID: 2604.09929v1

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

摘要: We establish the existence of a positive bounded weak solution for a class of Kirchhoff-type $p(\cdot)$ -Laplacian problems involving an arbitrary growth and a sandwich-type growth $s(\cdot)\in (\inf p,\sup p)$ . This setting leads to substantial analytical difficulties in the variational analysis of the associated energy functional. By combining truncation arguments with a priori estimates, we prove the existence result under suitable assumptions on the data.

13. 量子环中严格关联的电子：从Kohn-Sham势到Kantorovich势

Strictly correlated electrons in a quantum ring: from Kohn-Sham to Kantorovich potentials

作者: Thiago Carvalho Corso

类别: math-ph, cond-mat.str-el, math.AP, quant-ph

发布日期: 2026-04-10

arXiv ID: 2604.09908v1

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

摘要: Our goal in this paper is twofold. First, we characterize the class of pairwise interactions for which the Seidl conjecture on the structure of optimal plans for the symmetric multimarginal optimal transport problem with one-dimensional marginal holds. This extends previous results by Colombo, De Pascale, and Di Marino [Can. Jou. Math., 67 (2015), https://doi.org/10.4153/CJM-2014-011-x], which treated the case of translation-invariant, convex and decreasing interactions. In particular, our results apply to physically relevant interactions for electrons living on a quantum ring. The second main goal of the paper is to rigorously derive the leading order asymptotics of the adiabatic connection potential for strongly interacting systems. More precisely, we show that for electrons in a quantum ring (or one-dimensional interval), not only the Lieb density functional converges to the optimal transport (or strictly correlated) functional in the semiclassical limit, but also the representing potential converges to a regular Kantorovich potential. As an intermediate step, we also extend previous results on the strongly interacting limit of the Lieb functional to periodic systems in arbitrary dimensions.

14. 一个高阶非齐次柯西-欧拉方程的特解

A particular solution of a higher-order non-homogeneous Cauchy-Euler equation

作者: Miloud assal, Skander Belhaj

类别: math.AP

发布日期: 2026-04-10

arXiv ID: 2604.09778v1

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

摘要: In this paper we introduce a new concept of atoms on discrete sets to develop an advanced method to find a particular solution for higher-order non-homogeneous Cauchy-Euler equations. The proposed method provides also an approximate solution by using approximate roots for the characteristic polynomial of the Cauchy-Euler equation.

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

1. 含被动示踪剂的不可压缩流的瞬时爆破#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 具有不同粘性的质量守恒Navier-Stokes/Allen-Cahn系统的锐界面极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. Hall磁流体与电子磁流体湍流的决定性波数#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 一类修正NSCH-Oldroyd系统的局部适定性：基于PINN的数值验证#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 旋转球面上二维欧拉方程中四面体非带状流的分歧#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 通过五维提升的轴对称约化构造三维Navier–Stokes方程的有限时间稳定奇性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 不可压缩Navier-Stokes方程含未知体力的初值反问题的时间维约化与Carleman压缩原理#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 高维随机Zakharov系统（$d \geq 4$）：局部适定性与噪声对散射的正则化效应#

摘要翻译#

2. 五阶KP方程的唯一延拓性及其在控制问题中的应用#

摘要翻译#

3. 欧氏空间中薛定谔方程的可观测性#

摘要翻译#

4. 一个双重非线性抛物方程的博弈论解释#

摘要翻译#

5. 平面非散度型椭圆方程的改进Hölder正则性#

摘要翻译#

6. 障碍条件下非局部非参数Plateau问题当 $s\searrow 0$ 时的渐近行为#

摘要翻译#

7. $p$-变分容量的内凝聚器与固定相位的几何约化#

摘要翻译#

8. 无需显式格林函数的傅里叶型位势理论#

摘要翻译#

9. 等离子体-点电荷模型：边界效应与整体适定性#

摘要翻译#

10. 欧氏凸区域上拉普拉斯算子Dirichlet特征值的一些普适不等式#

摘要翻译#

11. 具有有界核与乘性噪声的粒子系统的混沌传播的定量估计#

摘要翻译#

12. 具有退化势的长程相共存模型#

摘要翻译#

13. 关于从局部到分数阶各向异性Calderón问题的稳定性传递（含外部测量）#

摘要翻译#

14. 一类余维数为一约束下的非线性热方程的收敛速率#

摘要翻译#

15. 流体力学中非局部拟线性演化方程的Schauder型估计与适定性#

摘要翻译#

16. 具有边界奇异漂移项的凸区域上Dirichlet Green函数#

摘要翻译#

17. 由图核调制的意见驱动型流行病动力学的动理学模型#

摘要翻译#

18. 三维径向数据下Klein-Gordon-Schrödinger系统的散射理论#

摘要翻译#

19. 圆平均的加权混合范数估计与波动方程的例外集估计#

摘要翻译#

20. 对齐动力学中的拉格朗日表述与欧拉闭合#

摘要翻译#

21. p-调和函数的 Harnack 不等式：通过 tug-of-war 方法改进维数依赖性#

摘要翻译#

22. Robin 边界条件下椭圆问题解的渐近行为#

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

1. 含被动示踪剂的不可压缩流的瞬时爆破

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 具有不同粘性的质量守恒Navier-Stokes/Allen-Cahn系统的锐界面极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. Hall磁流体与电子磁流体湍流的决定性波数

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 一类修正NSCH-Oldroyd系统的局部适定性：基于PINN的数值验证

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 旋转球面上二维欧拉方程中四面体非带状流的分歧

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 通过五维提升的轴对称约化构造三维Navier–Stokes方程的有限时间稳定奇性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 不可压缩Navier-Stokes方程含未知体力的初值反问题的时间维约化与Carleman压缩原理

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 高维随机Zakharov系统（ $d \geq 4$ ）：局部适定性与噪声对散射的正则化效应

摘要翻译

2. 五阶KP方程的唯一延拓性及其在控制问题中的应用

摘要翻译

3. 欧氏空间中薛定谔方程的可观测性

摘要翻译

4. 一个双重非线性抛物方程的博弈论解释

摘要翻译

5. 平面非散度型椭圆方程的改进Hölder正则性

摘要翻译

6. 障碍条件下非局部非参数Plateau问题当 $s\searrow 0$ 时的渐近行为

摘要翻译

7. $$p$$ -变分容量的内凝聚器与固定相位的几何约化

摘要翻译

8. 无需显式格林函数的傅里叶型位势理论

摘要翻译

9. 等离子体-点电荷模型：边界效应与整体适定性

摘要翻译

10. 欧氏凸区域上拉普拉斯算子Dirichlet特征值的一些普适不等式

摘要翻译

11. 具有有界核与乘性噪声的粒子系统的混沌传播的定量估计

摘要翻译

12. 具有退化势的长程相共存模型

摘要翻译

13. 关于从局部到分数阶各向异性Calderón问题的稳定性传递（含外部测量）

摘要翻译

14. 一类余维数为一约束下的非线性热方程的收敛速率

摘要翻译

15. 流体力学中非局部拟线性演化方程的Schauder型估计与适定性

摘要翻译

16. 具有边界奇异漂移项的凸区域上Dirichlet Green函数

摘要翻译

17. 由图核调制的意见驱动型流行病动力学的动理学模型

摘要翻译

18. 三维径向数据下Klein-Gordon-Schrödinger系统的散射理论

摘要翻译

19. 圆平均的加权混合范数估计与波动方程的例外集估计

摘要翻译

20. 对齐动力学中的拉格朗日表述与欧拉闭合

摘要翻译

21. p-调和函数的 Harnack 不等式：通过 tug-of-war 方法改进维数依赖性

摘要翻译

22. Robin 边界条件下椭圆问题解的渐近行为