生成时间: 2026年07月29日 05:41:02

运行状态: 完成

计划处理论文数量: 27

已完成论文数量: 27

未完成论文数量: 0

重点关注论文数量: 8

了解领域论文数量: 7

不相关论文数量: 12

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

1. 非负Ricci曲率流形上拟线性Liouville方程的Sharp刚性

Sharp rigidity for quasilinear Liouville equation on manifolds with nonnegative Ricci curvature

作者: Xiaohan Cai

类别: math.AP, math.DG

发布日期: 2026-07-28

arXiv ID: 2607.25981v1

摘要翻译: 我们在具有非负Ricci曲率的完备非紧Riemann流形上研究拟线性Liouville方程

Comment: 18 pages. Comments are welcome

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

1. 研究对象和背景

本文研究完备非紧黎曼流形 $(M^n, g)$ 上的拟线性 Liouville 方程

该方程是经典二维 Liouville 方程 $-\Delta u = e^u$ 的 conformally invariant 高维推广，几何上对应于寻找具有常数量曲率 $R_g = n(n-1)$ 的共形度量 $\tilde{g} = e^{\frac{2}{n-2}u} g$（当 $n>2$）或 Gaussian 曲率 $K_{e^u g_0} = \frac{1}{2}$（当 $n=2$）。在欧氏空间 $\mathbb{R}^n$ 中，Esposito [23] 在有限体积假设 $\int_{\mathbb{R}^n} e^u < +\infty$ 下完全分类了 $C^3$ 解：所有解必为标准气泡解

此结果本质依赖于 Pohozaev 恒等式、等周不等式与精细渐近分析。然而，在一般非负 Ricci 曲率流形上，该分类问题呈现深刻张力：一方面，存在大量非平坦流形（如圆锥、抛物型 warped product） admitting满足 $-\Delta_n u = e^u$ 的解，其渐近行为可显著偏离欧氏气泡；另一方面，Bishop–Gromov 比较定理、Sharp Isoperimetric Inequality（Brendle [3]）等几何不等式提供了强有力的 rigidity 工具。关键科学空白在于：是否存在仅依赖于解 $u$ 的渐近下界（而非全局可积性）的尖锐判据，足以强制流形平坦化并唯一确定解？ 特别地，已知在 $n=2$ 时，Catino–Monticelli [8] 要求 $u(x) \ge -4 \log r(x) - 2\gamma \log \log r(x), \gamma < 1$；Ciraolo–Farina–Polvara [20] 利用 $P$-函数法将余项放宽至 $u(x) \ge -4 \log r(x) - \log F(r)$，其中 $\int_2^{+\infty} (t F(t))^{-1} dt = +\infty$；Ou [34] 与 Ciraolo–Farina–Gatti [19] 进一步允许 $F^\beta$ 型余项。对 $n \ge 3$，Sun–Wang [35] 和 Catino–Monticelli–Roncoroni [9] 分别要求 $u(x) \ge -\frac{n^2}{n-1} \log r(x) - \frac{n}{2(n-1)} \log F(r)$ 或 $u(x) \ge -\frac{n^2}{n-1} \log r(x) + o(\log r(x))$ 且附加 $\int_M e^u < +\infty$。但所有前述工作均未回答：若完全放弃 $\int_M e^u < +\infty$ 假设，仅凭最优对数下界能否实现刚性？该下界系数 $-\frac{n^2}{n-1}$ 是否不可改进？余项 $o(\log r(x))$ 是否为自然边界？ 本文旨在填补这一核心空白，建立首个不依赖有限体积假设的尖锐刚性定理，并系统刻画解的亚对数衰减与流形渐近几何（通过渐近体积比 $AVR(M,g)$）之间的双向定量联系。

2. 主要定理或主要结果

定理 1.2（尖锐刚性）：设 $(M^n, g)$ 为 $n \ge 2$ 维完备非紧流形，满足 $\operatorname{Ric}_g \ge 0$；$u \in C^3(M^n)$ 满足 $-\Delta_n u = e^u$。若存在某点 $p \in M$，使得距离函数 $r(x) = d_g(x,p)$ 满足

则 $(M^n, g)$ 必等距同构于欧氏空间 $\mathbb{R}^n$，且 $u$ 必为标准气泡解（即形如 (1.3)）。此外，系数 $-\frac{n^2}{n-1}$ 是最优的：对任意 $\gamma > \frac{n^2}{n-1}$，存在非平坦完备流形 $(M^n_\gamma, g_\gamma)$ 满足 $\operatorname{Ric}_{g_\gamma} \ge 0$ 及解 $u_\gamma$，使得 $\lim_{r(x)\to+\infty} u_\gamma(x)/\log r(x) = -\gamma$，且 $AVR(M^n_\gamma, g_\gamma) \in (0,1)$，$\int_{M_\gamma} e^{u_\gamma} = |S^{n-1}| AVR(M^n_\gamma, g_\gamma) \gamma^{n-1}$。命题 1.3（尖锐体积估计）：在相同几何假设下，若 $u$ 满足 $u(x) \ge -\gamma \log r(x) + o(\log r(x))$，则

引理 2.2（推广形式）：若 $-\Delta_n u = f$（$f \ge 0$），且 $\liminf_{r(x)\to+\infty} u(x)/\log r(x) > -\infty$，则

定理 3.3 与 3.4（逆命题部分成立）：定义 $L := \liminf_{r(x)\to+\infty} u(x)/\log r(x)$。当 $n=2$ 时，若 $L = -\infty$，则 $\int_M e^u / AVR(M^2,g) = +\infty$（Conjecture 3.1 成立）。当 $n \ge 3$ 时，若加强假设为 $\lim_{r(x)\to+\infty} u(x)/\log r(x) = -\infty$（即极限存在且为 $-\infty$），则同样有 $\int_M e^u / AVR(M^n,g) = +\infty$。

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

证明的核心障碍在于：在无 $\int_M e^u < +\infty$ 前提下，如何从纯渐近下界导出解的可积性及流形刚性。作者的真正创新点在于构建了一个全新的、基于 $n$-容量（$n$-capacity）与 Bishop–Gromov 定理耦合的尖锐体积估计框架。具体而言：(1) 关键引理 2.1（源自 [9]）提供了一个适用于 $n$-拉普拉斯算子的广义容量控制：对光滑区域 $K \Subset \Omega$，有

其中 $\operatorname{Cap}_n(K, \Omega) = \inf\left\{ \int_\Omega |\nabla \phi|^n : \phi \in W^{1,n}_0(\Omega), \phi \ge 1 \text{ on } K \right\}$。此不等式将源项积分与解的振幅及几何容量关联。(2) 作者设计了一个精巧的径向测试函数 $\phi(x) = \frac{\int_{r(x)}^{R_2} A^{-1/(n-1)}(t) dt}{\int_{R_1}^{R_2} A^{-1/(n-1)}(t) dt}$（其中 $A(t) = \operatorname{Area}(\partial B(p,t))$），并利用 Bishop–Gromov 定理 $A(t) \le |S^{n-1}| t^{n-1}$ 导出容量上界

(3) 将此容量上界与渐近下界 $u(x) \ge -\gamma \log r(x) + o(\log r(x))$ 结合，代入引理 2.1，取极限 $R_2 \to +\infty$ 得到命题 1.3 的尖锐上界。此步骤是技术核心：它首次将对数衰减率 $\gamma$ 直接、定量地嵌入体积上界，且系数 $\gamma^{n-1}$ 与 $AVR$ 的乘积不可改进。(4) 刚性证明则依赖于 Brendle 的尖锐等周不等式：对水平集 $\Omega_t = \{u > t\}$，有

作者将此不等式与由方程导出的恒等式 $\int_{\Omega_t} e^u = \int_{\partial \Omega_t} |\nabla u|^{n-1}$ 及 Coarea 公式结合，推导出

对 $t$ 积分后，结合命题 1.3 的上界与该下界，迫使所有不等式取等号，从而触发 Brendle 等周不等式的刚性情形——即 $(M^n,g)$ 必为 $\mathbb{R}^n$。整个框架的创新性在于：它绕开了传统 Pohozaev 或 $P$-函数方法对有限体积或额外正则性的依赖，转而以 $n$-容量为桥梁，将解的渐近行为、几何比较定理与尖锐等周不等式三者精密耦合，实现了从纯渐近信息到全局刚性的跃迁。

4. 与之前工作的比较

与既有工作相比，本文的推进体现在三个根本维度：(1) 假设的彻底弱化：此前所有 $n \ge 3$ 的刚性结果（[35], [9]）均强制要求 $\int_M e^u < +\infty$，而本文定理 1.2 完全取消此条件，仅需 $u$ 的对数下界。即使在 $n=2$，本文亦改进了 [5] 的结果（其仍需有限体积），成为首个无 $L^1$ 假设的二维刚性定理。(2) 尖锐性的实质性突破：系数 $-\frac{n^2}{n-1}$ 的最优性由 [9, Example 3.4] 构造的 warped product 流形 $(M^n_\gamma, g_\gamma)$ 严格验证，该构造显示：当 $\gamma > \frac{n^2}{n-1}$ 时，存在非平坦流形承载满足 $u(x) \sim -\gamma \log r(x)$ 的解，且 $\int_{M_\gamma} e^{u_\gamma} = |S^{n-1}| AVR \gamma^{n-1}$，这表明本文系数不可改进。余项 $o(\log r(x))$ 同样是自然边界：若替换为 $O(1)$，则 [9, Example 4.7] 给出反例（$\liminf u/\log r = -\infty$ 但 $\int_M e^u < +\infty$ 且 $AVR=0$），故 $o(\log r)$ 是保证刚性成立的最小余项。(3) 方法论的范式转换：以往工作（[20], [19], [35]）主要依赖 $P$-函数法或 Pohozaev 恒等式，其本质是微分恒等式技巧，对解的全局可积性或特定对称性敏感；而本文首创 $n$-容量与 Bishop–Gromov–Brendle 三位一体的分析框架，将问题转化为几何测度论与非线性位势理论的交叉。代价在于：该方法目前尚不能处理 $L = -\infty$ 但极限不存在的情形（即 Conjecture 3.1 的完整证明仍开放），且对 Ricci 曲率负下界的流形尚未推广。相较 [23] 的欧氏分类，本文将尖锐性从 $\mathbb{R}^n$ 推广至非负 Ricci 流形，揭示了 $AVR$ 作为几何‘尺度因子’的核心作用；相较 [3] 的等周刚性，本文首次将其与非线性椭圆方程的解的渐近行为直接绑定，建立了 $AVR$ 与 $\int e^u$ 的定量桥梁。

2. 用微分方程模拟扩散过程

Mimicking diffusion processes with differential equations

作者: Rama Cont

类别: math.PR, math.AP, math.CA

发布日期: 2026-07-28

arXiv ID: 2607.25685v1

摘要翻译: 与扩散过程相关的概率流常微分方程（PF-ODE）在基于分数的生成建模中被广泛用作确定性采样器，以复现该扩散过程的边缘分布。尽管这一边缘匹配性质有时被视为理所当然，但其有效性依赖于一个常微分方程的适定性，该方程的速度场由扩散过程的分数函数（score function）构造而成。我们在漂移项与分数函数的弱正则性假设下，考察扩散过程、Fokker–Planck 方程及相应 PF-ODE 之间的精确数学关系。

Comment: 41 pages including appendix, 2 figures

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

1. 研究对象和背景

本文研究的核心对象是如下 $d$ 维前向扩散过程 $X = (X_t)_{t \in [0,T]}$ 满足的随机微分方程（SDE）：

其中 $W_t$ 是标准布朗运动，扩散系数 $\sigma : [0,T] \to (0,\infty)$ 为标量且满足 $0 < \epsilon \le \sigma(t) \le \epsilon^{-1}$（假设 (E1)），漂移场 $f : \mathbb{R}^d \times [0,T] \to \mathbb{R}^d$ 允许高度不规则——仅要求其散度 $\nabla \cdot f$ 属于 $L^1([0,T]; L^1_{\mathrm{loc}}(\mathbb{R}^d))$ 且负部 $[\nabla \cdot f]_- \in L^1([0,T]; L^\infty(\mathbb{R}^d))$（即假设 (D2)），不要求 $f$ 本身可弱微分或满足 Lipschitz 条件。初始分布 $p_0 \in L^1 \cap L^\infty(\mathbb{R}^d)$，无对数凹性、正则性或紧支集等先验限制（假设 (I1)）。该 SDE 的边际密度 $p_t(x) := \mathrm{Law}(X_t)(dx)/dx$（若存在）满足 Fokker–Planck 方程（Kolmogorov 前向方程）：

在 $p_t > 0$ 的区域上，利用恒等式 $\Delta p_t = \nabla \cdot (p_t \nabla \log p_t)$，可将 (5) 重写为连续性方程（transport equation）：

其中速度场 $v$ 定义为概率流速度（probability-flow velocity）：

对应地，确定性常微分方程（PF-ODE）为：

其核心动机源于 score-based 生成建模：若 PF-ODE (8) 良适（well-posed），则由初始随机变量 $Z_0 \sim p_0$ 出发生成的轨迹 $Z_t = Z(t,Z_0)$ 满足 $\mathrm{Law}(Z_t) = \mathrm{Law}(X_t) = p_t(x)\, dx$，从而以确定性方式“模拟”（mimic）原始扩散的边际流。这一性质被广泛用于高效采样与精确似然计算，但其数学基础却常被默认成立。本文旨在严格厘清三者之间的逻辑依赖：SDE 的解 $X_t$、Fokker–Planck PDE (5) 的解 $p_t$、以及 PF-ODE (8) 的流 $Z(t,\cdot)$。关键科学空白在于：在生成模型实际使用的弱正则性场景下（如非对数凹数据、奇异初值、神经网络学习的粗糙 score），经典 Cauchy–Lipschitz 理论失效，而现有工作（如 Mooney et al. [26], Stéphanovitch [33]）均建立在 $\nabla \log p_t$ 一致 Lipschitz 或 Hessian 有界的强假设之上，无法覆盖真实训练中 score 的低正则性。本文填补的空白正是：在仅假设 $\nabla \log p_t$ 属于 Sobolev 空间 $W^{1,1}_{\mathrm{loc}}$ 或有界变差空间 $BV_{\mathrm{loc}}$ 的极弱条件下，系统建立 PF-ODE 的 Lagrangian 描述（即正则 Lagrangian 流）的存在性、唯一性与稳定性，并揭示其与 Eulerian 描述（Fokker–Planck 解）的本质非等价性——前者要求额外的散度控制与 score 增长性条件，而后者仅依赖漂移散度的单侧界。

2. 主要定理或主要结果

本文主要结果分为四个层次，均在明确假设下给出严格证明：

1. Eulerian 层面（定理 3.1）：在假设 (E1)、(DE)（即 (D0)+(D2)+(D3)）和 (I1) 下，Fokker–Planck 方程 (5) 在能量类 $X := \{ p \in L^\infty([0,T]; L^1 \cap L^\infty(\mathbb{R}^d)) : \nabla p \in L^2([0,T]; L^2(\mathbb{R}^d)) \}$ 中存在唯一弱解 $p_t$，且该解恰为 SDE (4) 的边际密度流。关键估计为：

   $$ \sup_{s \le T} \|p_s\|_{L^2}^2 \le \|p_0\|_{L^2}^2 e^{\theta}, \quad \epsilon^2 \int_0^T \|\nabla p_s\|_{L^2}^2 \, ds \le \|p_0\|_{L^2}^2 (1 + \theta e^{\theta}), $$

   其中 $\theta := \int_0^T \|[\nabla \cdot f(\cdot,s)]_-\|_{L^\infty} \, ds$。特别地，该唯一性证明仅依赖 $\nabla \cdot f$ 的负部，无需 $f$ 的任何梯度信息或 score 正则性。

2. Lagrangian 层面（定理 4.5）：在 Eulerian 结果基础上，附加 score 正则性假设 (S)（即 (S1)-(S4)），其中最关键的是 (S4)：$|\nabla \log p_s(x)| \le A(s)(1+|x|)$，$A \in L^1(0,T)$，以及散度控制 $\nabla \cdot v \in L^1([0,T]; L^\infty(\mathbb{R}^d))$，则 PF-ODE (8) 存在唯一的正则 Lagrangian 流 $Z(t,\cdot) : \mathbb{R}^d \to \mathbb{R}^d$，满足：对任意可测集 $E \subset \mathbb{R}^d$，有 $(Z(t,\cdot))_\# (p_0 dx) = p_t dx$，即 $\mathrm{Law}(Z(t,Z_0)) = p_t dx$。推论 4.12 进一步指出，若加强为双向散度界 $|\nabla \cdot v| \in L^1([0,T]; L^\infty(\mathbb{R}^d))$，则流 $Z(t,\cdot)$ 可逆，且其逆亦为正则 Lagrangian 流。

3. 反例构造（第 6 节）：存在具体反例（构造于一维），使得定理 3.1 的条件全部满足，Fokker–Planck 方程 (5) 在 $X$ 中有唯一解 $p_t$，但 PF-ODE (8) 在 $t=0$ 处不存在正则 Lagrangian 流。其机制是 score 在初始支撑集外出现 $1/t$ 奇异性，导致 (S4) 失败而 Fisher 信息积分 $\int_0^T I(p_s)\, ds < \infty$（即 (S4′)）仍成立，从而严格区分了两种正则性条件。

4. 学习 score 的稳定性（定理 7.1）：设学习得到的 score $s_\theta$ 诱导速度场 $v_\theta(x,t) := f(x,t) - \frac{1}{2}\sigma(t)^2 s_\theta(x,t)$，若 $v_\theta$ 满足 (S) 类似条件，则对对应流 $Z_\theta$ 与真流 $Z$，有定量对数 Lipschitz 稳定性估计：

   $$ \left| \log \frac{|Z_\theta(t,x) - Z_\theta(t,y)|}{|Z(t,x) - Z(t,y)|} \right| \le C \int_0^t \|v_\theta(\cdot,s) - v(\cdot,s)\|_{L^1(p_s)} \, ds, $$

   其中 $C$ 为仅依赖 $\sigma$ 和 $A$ 的常数。命题 7.4 将此误差与 score-matching 损失 $\mathbb{E}_{p_t}[|s_\theta - \nabla \log p_t|^2]$ 关联，引入逆密度因子 $1/p_t$；命题 7.8 在 $s_\theta$ Lipschitz 时给出线性速率估计 $\|Z_\theta(t,\cdot) - Z(t,\cdot)\|_{L^\infty} \le C' \int_0^t \|s_\theta - \nabla \log p_s\|_{L^\infty} \, ds$。

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

本文的核心技术障碍在于：PF-ODE (8) 的速度场 $v$ 显式依赖于未知密度 $p_t$ 的对数梯度 $\nabla \log p_t$，而后者在生成模型实践中常由神经网络粗糙估计，导致 $v$ 缺乏经典 ODE 理论所需的 Lipschitz 连续性。作者的真正创新点是系统性地将 DiPerna–Lions–Ambrosio 正则 Lagrangian 流理论嫁接到扩散模型分析中，并据此构建一个精细的假设层级，精确刻画每个条件的作用。

第一，针对 Eulerian 问题（定理 3.1），作者设计了一个精巧的能量估计框架，其关键在于：在测试函数选取 $\phi_R w^2$（$w = p_1 - p_2$）后，通过分布散度恒等式 $\int f \cdot \nabla \psi = -\int (\nabla \cdot f) \psi$ 将漂移项 $\int f \cdot w \nabla w$ 完全转化为 $\int (\nabla \cdot f) w^2$。这一步骤仅需 $\nabla \cdot f$ 的绝对连续性（假设 (9)）及其负部的 $L^\infty$ 可积性（(D2)），彻底规避了对 $\nabla f$ 或 $\nabla \log p_t$ 的任何要求。该估计的鲁棒性体现在：即使 $f$ 仅为 $L^2_{\mathrm{loc}}$ 且无弱导数，只要其散度作为 Radon 测度的奇异部分为零，就能保证唯一性。

第二，针对 Lagrangian 问题（定理 4.5），作者严格应用 Ambrosio 的 BV 理论框架。其核心工具是正则 Lagrangian 流的存在唯一性定理，该定理要求速度场 $v$ 满足：(i) $v \in L^1([0,T]; W^{1,1}_{\mathrm{loc}}(\mathbb{R}^d; \mathbb{R}^d))$，(ii) $\nabla \cdot v \in L^1([0,T]; L^\infty(\mathbb{R}^d))$。作者将这些条件翻译为对 score $\nabla \log p_t$ 的可验证假设：(S2) 保证 $v$ 的 Sobolev 正则性，(S4) 与 (D4) 共同确保 $v$ 的线性增长性以满足流的全局存在，而 (S3) 与 (D2) 的组合则导出 $\nabla \cdot v$ 的 $L^\infty$ 控制。特别是，作者证明 $\nabla \cdot v = \nabla \cdot f - \frac{1}{2} \sigma(t)^2 \Delta \log p_t$，其中 $\Delta \log p_t$ 的正部由 (S3) 控制，负部则由 $\nabla \cdot f$ 的负部吸收，这是连接 Eulerian 与 Lagrangian 的关键桥梁。

第三，反例构造（第 6 节）并非简单否定，而是精密的病理学分析。作者在一维构造初值 $p_0 = \mathbf{1}_{[0,1]}$ 与漂移 $f(x,t) = -x/(2t)$，导出显式解 $p_t(x) = \frac{1}{\sqrt{2\pi t}} e^{-(x-1)^2/(2t)}$，其 score $\nabla \log p_t(x) = -(x-1)/t$ 在 $t \to 0^+$ 时于 $x \neq 1$ 处爆炸为 $1/t$，违反 (S4)，但 Fisher 信息 $I(p_t) = 1/t$ 可积，故 (S4′) 成立。此反例严格证明：Fokker–Planck 解的唯一性绝不蕴含 PF-ODE 流的存在性，从而确立了 Eulerian 与 Lagrangian 描述的根本分野。

第四，稳定性分析（第 7 节）采用双尺度论证。定理 7.1 的证明基于正则 Lagrangian 流的比较原理，其核心引理（引理 7.5）建立了流差的对数 Lipschitz 控制：$\frac{d}{dt} \log |Z_\theta - Z| \le C |v_\theta - v| \circ Z$，再通过对 $p_t$-加权的 $L^1$ 误差进行时间积分。命题 7.4 的创新在于将 $L^1(p_t)$ 误差与 score 的 $L^2(p_t)$ 误差通过 Hölder 不等式与逆密度 $1/p_t$ 关联，这直接反映了学习 score 在低密度区域的高敏感性，为实际训练提供了理论警示。

4. 与之前工作的比较

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

1. 正则性假设的维度：Mooney et al. [26] 和 Stéphanovitch [33] 的工作均建立在 $\nabla \log p_t$ 一致 Lipschitz 的强假设上，即 $|\nabla \log p_t(x) - \nabla \log p_t(y)| \le L |x-y|$，这隐含了 $\nabla \log p_t$ 的有界性及 $\Delta \log p_t$ 的 $L^\infty$ 有界性，本质上要求数据分布高度光滑且对数凹。本文则将正则性门槛大幅降低至 Sobolev 类 $W^{1,1}_{\mathrm{loc}}$ 或 BV 类，仅需 $\nabla \log p_t$ 的一阶弱导数局部可积，且允许其在空间上无界（仅受线性增长约束 (S4)）。例如，对于紧支集初值 $p_0$，Proposition 5.3 证明早期停止（early stopping）后 $p_t$ 自动获得所需正则性，这在图像生成中对应噪声注入初期，而经典 Lipschitz 分析在此阶段失效。因此，本文适用范围显著拓宽，覆盖了更真实的生成场景。

2. 数学框架的维度：Han et al. [17] 和 Huang et al. [19] 的工作聚焦于“假设 PF-ODE 已良适”的前提下，分析 score 估计误差与离散化误差的传播，其结论是误差的 $L^2$ 或总变差范数界。本文则直击根基，回答“PF-ODE 是否能定义”的本体论问题。通过引入 DiPerna–Lions–Ambrosio 理论，本文首次将生成模型的确定性采样器置于现代偏微分方程正则性理论的严格框架下，其贡献不仅是技术推广，更是范式转换：将生成建模的数学基础从经典 ODE 理论提升至非光滑动力系统理论。

3. 结论强度与代价的维度：本文的推进以明确的代价为前提。定理 4.5 要求 (S4) 即 $|\nabla \log p_t(x)| \le A(t)(1+|x|)$，这是一个比 Fisher 信息有限（(S4′)）更强的条件，反例 6.1 证实了此差距的实质性。相比之下，[26, 33] 的 Lipschitz 假设自动蕴含 (S4)，但本文证明 (S4) 已是最小充分条件。此外，本文的稳定性估计（定理 7.1）依赖于 $L^1(p_t)$ 速度误差，而 [19] 的 TV 收敛率依赖于 $L^\infty$ score 误差，二者不可直接比较：前者对低密度区域的误差更敏感（因 $p_t$ 小），后者则对高密度区域更敏感。这种差异揭示了不同误差度量背后的物理含义——本文的估计更贴近实际采样中“流线分离”的几何直观，而 [19] 的估计更侧重于分布整体的接近程度。综上，本文并非取代先前工作，而是为其提供底层支撑：当 score 学习足够好以满足 (S) 时，[17, 19] 的误差传播分析才真正有效；而当 (S) 失败时，本文的反例警告我们，任何基于 PF-ODE 的确定性采样都可能从根源上崩溃。

3. 湍流火焰速度可在曲率平滑作用下增大

Turbulent Flame Speed Can Increase under Curvature Smoothing

作者: Hung V. Tran, Jack Xin, Yifeng Yu

类别: math.AP, physics.flu-dyn

发布日期: 2026-07-28

arXiv ID: 2607.25185v1

摘要翻译: 曲率效应预期会平滑火焰锋面褶皱，从而降低湍流火焰速度。我们构造了一种光滑的三维周期剪切流，使得在该流中引入 Markstein 曲率扩散项反而增大了由水平集 G-方程 $G_t + \mathbf{V} \cdot \nabla G = s_L |\nabla G|$ 预测的有效火焰速度。这是该模型中首个关于曲率平滑导致单调减速的反例，并与二维剪切流中已建立的严格单调性结果形成对比。该例子揭示了一种真正的多维机制：局部曲率平滑可增强而非抑制大尺度锋面传播。

Comment: 14 pages, 1 figure

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

1. 研究对象和背景

本文研究三维周期剪切流中带 Markstein 曲率修正的 G-方程所描述的火焰前缘传播问题。G-方程是湍流燃烧中刻画火焰面运动的经典水平集模型，其无曲率版本为一阶哈密顿–雅可比型偏微分方程：

其中 $V(x) = (0,0,f(x_1,x_2))$ 是光滑 $\mathbb{T}^2$-周期剪切流，$s_L > 0$ 为层流火焰速度，$G(x,t)$ 的零水平集 $\{G=0\}$ 表示瞬时火焰面。当计入曲率效应时，局部燃烧速度被修正为 $s_L(1 - d \delta_L \kappa)_+$，其中 $d > 0$ 为 Markstein 数，$\delta_L$ 为火焰厚度，$\kappa$ 为火焰面前缘的平均曲率（即 $\kappa = \operatorname{div}\left( \frac{DG}{|DG|} \right)$），$(\cdot)_+$ 表示正部截断以避免负速度。在均匀化尺度下，令 $V(x) \mapsto V(x/\varepsilon)$，并设 $\delta_L \sim \varepsilon$，则曲率项贡献为 $O(d\varepsilon)$，其宏观效应由对应胞问题的遍历常数 $c(d)$ 刻画。

具体地，对传播方向 $P = (p,1) \in \mathbb{R}^3$（其中 $p = (2,2) \in \mathbb{R}^2$），经变量变换后，有效燃烧速度 $c(d)$ 由如下二维周期胞问题的遍历常数给出（非截断情形）：

该方程等价于一个退化二阶拟线性椭圆型方程，其主部算子在 $d=0$ 时完全退化为一阶输运算子，而在 $d>0$ 时引入了沿法向的各向异性扩散；但该扩散并非各向同性拉普拉斯型，而是仅作用于曲率方向，且其系数依赖于梯度 $Dv$，导致其为强非线性、非一致椭圆的退化结构。特别地，当 $d=0$ 时，方程退化为 $\frac{p}{\sqrt{1 + |p + Dv|^2}} + f(x) = c(0)$，其解可能不唯一且不可微（见 [Tra21, Section 4.8]），而 $c(0)$ 对应于无曲率时的湍流火焰速度。

该问题的核心物理与数学动机在于：燃烧文献中长期认为正 Markstein 数 $d>0$ 总会通过平滑火焰褶皱而抑制湍流增强，即 $c(d)$ 应关于 $d$ 单调递减。这一猜想已在二维剪切流中被严格证明 [LXY18]：此时胞问题可约化为一维 ODE，曲率项表现为真实粘性，强制解更光滑，从而降低有效速度。然而，在三维设置下，胞问题本质是二维的，无法降维，曲率扩散不再具有全局正则化能力，其退化性导致其对宏观传播的影响可能反转。本文旨在填补这一关键空白：是否存在一个光滑三维剪切流 $V$，使得 $c(d)$ 在 $d=0$ 附近非单调——即存在 $d_0 > 0$ 使得 $c(d_0) > c(0)$？这直接挑战了曲率总是抑制湍流传播的传统范式，并触及哈密顿–雅可比方程受小退化二阶扰动时的宏观选择原理这一深层分析问题。

2. 主要定理或主要结果

定理 1（主结果）：取 $p = (2,2) \in \mathbb{R}^2$，存在 $f \in C^\infty(\mathbb{T}^2)$ 满足： (i) $\max_{\mathbb{T}^2} f < 0$； (ii) $c(0) = 0$； (iii) 右导数存在且为正：

其中 $\lambda$ 由变分公式（引理 1）给出：

且 $u(x) = p \cdot x + v(x)$，$v(x) = a(x_1) + a(x_2)$，$a \in C^\infty(\mathbb{R})$ 为 1-周期函数，满足 $\|a'\|_{L^\infty} \leq 1$。进一步，$f$ 被显式构造为 $f(x) = -\sqrt{1 + |Du(x)|^2}$，从而保证 $c(0) = 0$。推论：存在 $d_0 > 0$，使得对所有 $d \in (0,d_0)$，有 $c(d) > 0 = c(0)$，且因 $c(d) - \max f > 0$，截断无效，故 $c(d) = \tilde{c}(d)$ 为真实湍流火焰速度。因此，$c(d)$ 在 $d=0$ 邻域内非单调递减，否定 [XYR24, Question 8]。此外，附录 A 构造了另一例子，表明非单调性不仅限于 $d \to 0^+$，在 $d \approx 1$ 处亦可发生。注：该结果对 $f$ 在 $C^0$-范数下是鲁棒的——若 $f$ 被足够小的 $C^0$ 扰动，则 $c(d^*) > c(0)$ 仍成立。

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

证明的核心障碍在于：$c(d)$ 关于 $d$ 的单调性无法通过对胞问题（2）形式求导获得，因为 $d=0$ 处的退化导致线性化方程（3）不具确定符号的结构；且二维情形下缺乏 ODE 约化，无法像 [LXY18] 那样获得 $c'(d)$ 的显式积分表达式。作者的真正创新点在于建立了一个基于轨道平均的变分表征与渐近展开框架，将 $c(d)$ 在 $d=0^+$ 的行为精确归结为曲率项 $g(x)$ 沿特征轨道的极值平均。

第一，构造特殊解 ansatz：选取 $v(x) = a(x_1) + a(x_2)$，使 $Du(x) = (2 + a'(x_1), 2 + a'(x_2))$，并令 $f(x) = -\sqrt{1 + |Du(x)|^2}$，从而强制 $c(0) = 0$。此构造确保 $|Du| \geq \sqrt{2}$，避免分母退化，并使特征场 $W(x) = Du / \sqrt{1 + |Du|^2}$ 具有明确的几何结构——其积分曲线 $\xi(t)$ 满足 $\dot{\xi} = W(\xi)$，且投影到环面 $\mathbb{T}^2$ 上具有旋转向量 $(1,1)$。

第二，建立变分公式（引理 1）：定义泛函 $\lambda = \inf_\psi \sup_x \{ W \cdot D\psi - g \}$，并证明其等于 $-\min_{\xi \in \Gamma} \frac{1}{T_\xi} \int_0^{T_\xi} g(\xi(t)) \, dt$，其中 $\Gamma$ 为所有 $W$-轨道构成的集合。该公式的关键在于引入坐标变换 $(r,\theta)$，其中 $r = F(x_1) - F(x_2)$，$F(s) = \int_0^s d\tau / P(\tau)$，$P(t) = 2 + a'(t)$。在此坐标下，$W \cdot D = A(r,\theta) \partial_\theta$，从而将泛函极值问题转化为沿 $r$-参数族轨道的平均值优化问题。

第三，证明 $c(d)/d \to \lambda$（引理 2）：采用精巧的扰动检验函数法。上界部分取 $\phi = v + d\psi$，利用一致展开得到 $H_d[\phi] \leq d \sup(W \cdot D\psi - g) + O(d^2)$，从而 $\limsup c(d)/d \leq \lambda$。下界部分需克服 $W \cdot D\psi - g = \lambda$ 可能无全局解的困难，作者构造 $\phi_d = v + \sqrt{d} \chi(r) + d\psi$，其中 $\chi(r)$ 为 $L$-周期函数，用以“抬升” $m(r)$ 在其最大值点附近的值，确保 $H_d[\phi_d] \geq d(\lambda - \varepsilon) + O(d^{3/2})$，最终证得 $\liminf c(d)/d \geq \lambda$。此两步共同确立 $c(d)$ 的一阶渐近：$c(d) = \lambda d + o(d)$。

第四，显式验证 $\lambda > 0$：选取双模小扰动 $a_\delta(t) = -\frac{\delta}{2\pi} \cos(2\pi t) - \frac{\delta}{4\pi} \cos(4\pi t)$，计算得 $g(x)$ 的显式表达式（第 7 页公式 (9)），并证明存在轨道 $\xi$ 使得 $\frac{1}{T_\xi} \int_0^{T_\xi} g(\xi(t)) \, dt < 0$，故 $\lambda > 0$。整个框架揭示了退化曲率扩散的宏观效应由特征轨道上的极值平均主导，而非点态平滑，这是对 vanishing-viscosity 选择原理在退化二阶扰动下的重要推广。

4. 与之前工作的比较

本文与既有工作存在三重根本性对比：

1. 维度与模型适用性：[LXY18] 在二维剪切流中证明了 $c'(d) < 0$，其核心在于胞问题可约化为一维 ODE：$\frac{p}{\sqrt{1 + (p + v')^2}} + f(x) = c(d)$，曲率项 $-d \frac{d}{dx}\left( \frac{v'}{\sqrt{1 + (v')^2}} \right)$ 是标准粘性项，强制 $v'$ 更光滑，从而抑制由剪切 $f$ 引起的速度增强。本文的三维剪切流对应二维胞问题，无此类 ODE 约化，曲率算子 $-d \operatorname{div}\left( \frac{p + Dv}{\sqrt{1 + |p + Dv|^2}} \right)$ 是退化、各向异性的，其扩散仅沿法向，且系数依赖于梯度，导致其无法全局正则化解，反而在特定几何配置下产生加速效应。这表明二维单调性不能外推至三维，凸显了多维性的本质作用。

2. 渐近极限的非交换性：[LXY18, Section 2.2] 的弱流摄动理论假设流动为 $V_\delta = V_0 + \delta V_1$，并在 Diophantine 方向 $P$ 下证明：对固定 $d > 0$，当 $\delta \to 0^+$ 时，$c_\delta(d,P) < c_\delta(0,P)$。本文构造的流虽也是弱扰动，但其方向 $P = (2,2,1)$ 是有理的，且关键在于揭示了两个小参数极限的非交换性：对固定小 $\delta$，当 $d \to 0^+$ 时，$c_\delta(d,P) > c_\delta(0,P)$；而对固定小 $d$，当 $\delta \to 0^+$ 时，$c_\delta(d,P) < c_\delta(0,P)$。因此，$\lim_{d \to 0^+} \lim_{\delta \to 0^+} c_\delta(d,P) \neq \lim_{\delta \to 0^+} \lim_{d \to 0^+} c_\delta(d,P)$，说明弱流与 vanishing-curvature 极限不可交换，这是一个 KAM 型小除数问题，此前未被系统探讨。

3. 宏观选择机制的深化：vanishing-viscosity Hamilton–Jacobi 方程的选择原理（如 [AIPSM05]）指出，当 Aubry 集由周期轨道构成时，小粘性主导的宏观效应由曲率项在这些轨道上的极值平均决定。本文首次将此原理应用于退化曲率扩散，并证明 $c'(0^+)$ 由 $-g(x)$ 的极小轨道平均给出，而非点态值。这推广了选择原理的适用范围：它不仅适用于非退化粘性，也适用于退化二阶扰动，且其选择准则从 Aubry 集的几何性质转为特征流的轨道动力学。代价是：该机制高度依赖于流的精细构造（如双模 $a_\delta$），目前尚无一般性判据预测何时出现加速；且结果局限于 $d \to 0^+$，对大 $d$ 的全局单调性仍未解决。从提供内容中无法确认该反例是否适用于更一般的三维流（如各向同性湍流），作者亦未显式展开与 [MMT+25] 中三维有效速度不存在性结果的比较。

4. 二维电阻性磁流体方程中具常磁场的周期剪切流的长波不稳定性

Long-wave instability of periodic shear flows with constant magnetic field for the 2D resistive MHD equations

作者: Roberto Feola, Luca Franzoi, Riccardo Montalto, Claudia Peña

类别: math.AP, math-ph

发布日期: 2026-07-28

arXiv ID: 2607.25836v1

摘要翻译: 我们在周期区域 ${\mathbb T}_α\times {\mathbb T} = \Big( {\mathbb R}/(\frac{2 π}α {\mathbb Z}) \times {\mathbb R}/(2 π{\mathbb Z}) \Big)$ 上，研究二维粘性、电阻性磁流体动力学（MHD）方程在涡量-电流形式下的长波线性稳定性与不稳定性，扰动围绕周期剪切流 $(U(y),0)$ 及其耦合的常背景磁场 ${\bf b}=({\rm b}_1,{\rm b}_2)$ 展开。这是 Colombo, Dolce, Montalto & Ventura 关于 Navier-Stokes 方程的近期工作在 MHD 框架下的非平凡推广，其精神承袭 Kolmogorov、Meshalkin、Sinai 与 Yudovich 的经典工作。我们建立了关于剪切流剖面 $U(y)$ 的显式条件，该条件涉及粘性系数 $ν$、电阻率 $η$ 以及背景磁场 ${\bf b}$ 的分量，从而在 $α\ll 1$ 区间内获得长波线性稳定性与不稳定性。证明结合了非微扰的正规形变换（将零 Fourier 模式与非零模式解耦）以及关于参数 $α$ 的、从零未扰动本征值分岔出的本征值的尖锐渐近展开。作为动力学后果，我们得到相空间的分裂：分解为不稳定子空间与稳定子空间，在其上解分别以 Sobolev 范数指数增长或指数衰减。

Comment: 无

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

1. 研究对象和背景

本文研究二维电阻性、粘性不可压缩磁流体（MHD）方程在涡度-电流（vorticity-current）形式下的线性长波稳定性问题。所考虑的方程组为

其中 $ (x,y) \in \mathbb{T}_\alpha \times \mathbb{T} = \left( \mathbb{R} / \left( \frac{2\pi}{\alpha} \mathbb{Z} \right) \times \mathbb{R} / (2\pi \mathbb{Z}) \right) $，$ \alpha \ll 1 $ 为长波参数（即逆纵横比），$ \nabla = (\partial_x, \partial_y) $，$ \nabla^\perp = (-\partial_y, \partial_x) $，$ \Delta = \partial_x^2 + \partial_y^2 $；$ u $ 为速度场，$ b $ 为磁场，$ w = \nabla^\perp \cdot u $ 为标量涡度，$ j = \nabla^\perp \cdot b $ 为标量电流；$ \nu > 0 $ 为运动粘性系数，$ \eta > 0 $ 为电阻率；非线性项 $ H(u,b) = \partial_x b \cdot \nabla u_2 - \partial_y b \cdot \nabla u_1 $。该模型描述强碰撞等离子体（如天体物理中日冕、吸积盘）的动力学行为，其数学难点在于双扩散耦合结构：涡度演化受粘性与磁场剪切共同驱动，而电流演化则受电阻性耗散与洛伦兹力耦合支配。

研究对象是围绕周期剪切流 $ u(y) = (U(y), 0) $ 与恒定背景磁场 $ \mathbf{b} = (b_1, b_2) \in \mathbb{R}^2 $ 的稳态解 $ (w,j) = (\omega(y), 0) $，其中 $ \omega(y) = -U'(y) $，且 $ U \in C^{S+2}(\mathbb{T}, \mathbb{R}) $ 满足零均值条件 $ \int_{\mathbb{T}} U(y)\,dy = 0 $。需强调的是，该稳态对仅在 $ b_2 = 0 $ 时严格满足无外力诱导方程（见 Remark 1.5）；当 $ b_2 \neq 0 $ 时，感应方程中产生误差项 $ b_2 U''(y) $，无法通过物理允许的电磁力源完全抵消——这使问题兼具数学挑战性与物理真实性。因此，本文同时处理 $ b_2 = 0 $（物理相关情形）与 $ b_2 \neq 0 $（数学拓展情形），以揭示磁场分量对不稳定性机制的本质影响。

核心问题是：在长波极限 $ \alpha \to 0^+ $ 下，扰动 $ (\Omega(t,x,y), J(t,x,y)) $ 满足缩放后的线性化系统

其中线性算子 $ L_{\nu,\eta,\varepsilon} : H^{s+2}(\mathbb{T}) \times H^{s+2}(\mathbb{T}) \to H^s(\mathbb{T}) \times H^s(\mathbb{T}) $ 显式给出为

该算子具有非平凡的非局部结构（含 $ (\varepsilon^2 - \partial_y^2)^{-1} $），且其频谱在 $ \varepsilon \to 0 $ 时发生退化：零阶近似 $ D_{KK} = \begin{pmatrix} \nu \partial_y^2 & b_2 \partial_y \\ b_2 \partial_y & \eta \partial_y^2 \end{pmatrix} $ 的谱为 $ \lambda_\pm(j) = -\frac{(\nu+\eta)j^2}{2} \pm \frac{1}{2} \sqrt{(\nu-\eta)^2 j^4 - 4b_2^2 j^2} $，在 $ j \neq 0 $ 处具有负实部，但 $ j = 0 $ 模式（即 $ x $-平均扰动）的谱完全坍缩至零，成为主导不稳定性的关键。已有工作如 Kolmogorov [1]、Meshalkin–Sinai [41]、Yudovich [48] 奠定了剪切流长波不稳定的物理图景，而 Colombo–Dolce–Montalto–Ventura [15] 首次对 Navier–Stokes 方程实现了严格的谱分析。本文的空白在于：如何将该框架推广至 MHD 系统，精确刻画粘性 $ \nu $、电阻率 $ \eta $、磁场分量 $ b_1, b_2 $ 三者之间的竞争关系，并建立可验证的显式判据，以区分稳定/不稳定区域。

2. 主要定理或主要结果

本文在 Sobolev 空间 $ H^s(\mathbb{T}) \times H^s(\mathbb{T}) $（$ 0 \le s \le S $，$ S \ge 1 $）上建立了如下主要结果：

定理 1.1（谱长波稳定性/不稳定性）：设 $ \nu, \eta > 0 $，$ U \in C^{S+2}(\mathbb{T}, \mathbb{R}) $，$ \varepsilon = \alpha k \in (0, \varepsilon_0) $，其中 $ \varepsilon_0 = \varepsilon_0(S, \nu, \eta, \mathbf{b}) > 0 $ 充分小。则 $ L_{\nu,\eta,\varepsilon} $ 的谱 $ \sigma(L_{\nu,\eta,\varepsilon}) $ 是离散的，且满足：

(i) 当 $ b_1 \neq 0 $ 且 $ \eta > \nu $ 时，若

则存在两个单重特征值 $ \mu_-, \mu_+ \in \mathbb{C} \setminus \{0\} $ 满足 $ \Re(\mu_\pm) > 0 $，其余谱点实部均为负；

(ii) 当 $ b_1 = 0 $ 时，若

则存在唯一单重特征值 $ \mu_+ \in \mathbb{C} \setminus \{0\} $ 满足 $ \Re(\mu_+) > 0 $，其余谱点实部均为负；

(iii) 当 $ b_1 \neq 0 $ 且 $ \nu \ge \eta $ 时，$ \sigma(L_{\nu,\eta,\varepsilon}) $ 中所有特征值实部均为负。

进一步，由引理 4.9 给出的渐近展开表明：不稳定特征值具有尺度 $ \mu_\pm = \varepsilon^2 \lambda_0 + O(\varepsilon^3) $，其中 $ \lambda_0 $ 是 $ 2 \times 2 $ 矩阵 $ M^{(0)}_{\nu,\eta} $ 的特征值，其显式依赖于 $ \nu, \eta, b_1, b_2 $ 和 $ U $ 的 Fourier 系数。

定理 1.3（动力学行为）：在定理 1.1 条件下，Cauchy 问题 $ \partial_t u(t) = L_{\nu,\eta,\varepsilon} u(t),\; u(0) = \varphi \in H^s \times H^s $ 的解满足相空间分裂：

(i) 在 (i) 条件下，$ H^s \times H^s = U_\varepsilon \oplus S_\varepsilon $，其中 $ \dim U_\varepsilon = 2 $，$ \operatorname{codim} S_\varepsilon = 2 $，且对任意 $ \varphi \in U_\varepsilon $，有 $ \|u(t)\|_s \gtrsim e^{\mu \varepsilon^2 t} \|\varphi\|_s $；对任意 $ \varphi \in S_\varepsilon $，有 $ \|u(t)\|_s \lesssim e^{-\mu_0 t} \|\varphi\|_s $；

(ii) 在 (ii) 条件下，$ \dim U_\varepsilon = 1 $，且对任意 $ \varphi \in U_\varepsilon $，有 $ \|u(t)\|_s \gtrsim e^{\mu_1 \varepsilon^2 t} \|\varphi\|_s $；对任意 $ \varphi \in S_\varepsilon $，有 $ \|u(t)\|_s \lesssim \varepsilon^{-1} e^{-\mu_2 \varepsilon^2 t} \|\varphi\|_s $；

(iii) 在 (iii) 条件下，对任意 $ \varphi \in H^s \times H^s $，有 $ \|u(t)\|_s \lesssim \varepsilon^{-1} e^{-\mu_3 \varepsilon^2 t} \|\varphi\|_s $。

所有常数 $ \mu, \mu_0, \mu_1, \mu_2, \mu_3 > 0 $ 仅依赖于 $ \nu, \eta, b_1, b_2, S $，且与 $ \varepsilon $ 无关。

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

证明的核心障碍源于 $ L_{\nu,\eta,\varepsilon} $ 在长波极限 $ \varepsilon \to 0 $ 下的非微扰性（non-perturbative nature）：其零模（$ j = 0 $）谱完全坍缩至零，且 $ \varepsilon $-展开中出现 $ \varepsilon^{-1} $ 和 $ \varepsilon $ 的混合项（见 (1.19) 式），标准摄动理论失效。作者的真正创新点在于构建了一个非微扰正规形变换（non-perturbative normal form transformation），其作用不是将算子化为对角形式，而是实现零模与非零模的精确解耦，从而分离出决定稳定性的主导部分。

具体而言，方法分为两步：

第一步：零模解耦（Decoupling of the zero mode）。利用投影 $ \Pi_0 $（取 $ x $-平均）和 $ \Pi_K^0 $（零均值部分），将 $ L_{\nu,\eta,\varepsilon} $ 写为块矩阵 $ \begin{bmatrix} \Pi_0 L \Pi_0 & \Pi_0 L \Pi_K^0 \\ \Pi_K^0 L \Pi_0 & \Pi_K^0 L \Pi_K^0 \end{bmatrix} $。目标是构造一个可逆变换 $ T = \begin{bmatrix} \mathrm{Id}_0 & T_2 \\ T_3 & \mathrm{Id}_K \end{bmatrix} $，使得共轭后 $ T L_{\nu,\eta,\varepsilon} T^{-1} = \begin{bmatrix} M^{(0)}_{\nu,\eta} & 0 \\ 0 & L^{(0)}_{\nu,\eta,\varepsilon} \end{bmatrix} $。关键在于求解两个算子方程：$ T_2 $ 满足 $ \Pi_0 L \Pi_K^0 + \Pi_0 L \Pi_0 T_2 - T_2 \Pi_K^0 L \Pi_K^0 = 0 $，$ T_3 $ 满足 $ \Pi_K^0 L \Pi_0 + \Pi_K^0 L \Pi_K^0 T_3 - T_3 \Pi_0 L \Pi_0 = 0 $。作者证明了 $ T_2 = O(\varepsilon^3) $，$ T_3 = O(\varepsilon^{-1}) $，故 $ T = \mathrm{Id} + O(\varepsilon^{-1}) $ “远离恒等映射”，这正是非微扰性的体现。该变换成功地将 $ j = 0 $ 模式编码进 $ 2 \times 2 $ 矩阵 $ M^{(0)}_{\nu,\eta} $，其特征值 $ \mu_\pm $ 的 $ \varepsilon^2 $-阶渐近展开（Proposition 4.9）直接导出了不稳定性判据 (1.9)–(1.10)。

第二步：非零模稳定性（Stability of non-zero modes）。解耦后的算子 $ L^{(0)}_{\nu,\eta,\varepsilon} $ 作用于 $ H^s_0 \times H^s_0 $，被分解为 $ L^{(0)}_{\nu,\eta,\varepsilon} = D_{KK} + Q^{(0)}_\varepsilon $，其中 $ D_{KK} = \begin{bmatrix} \nu \partial_y^2 & b_2 \partial_y \\ b_2 \partial_y & \eta \partial_y^2 \end{bmatrix} $ 是傅里叶乘子，其谱 $ \lambda_\pm(j) = -\frac{(\nu+\eta)j^2}{2} \pm \frac{1}{2} \sqrt{(\nu-\eta)^2 j^4 - 4b_2^2 j^2} $ 对所有 $ j \neq 0 $ 满足 $ \Re(\lambda_\pm(j)) \le -\sigma |j|^2 $（$ \sigma > 0 $）。关键引理 5.2 证明了扰动 $ Q^{(0)}_\varepsilon $ 是 $ O(\varepsilon) $-阶有界算子，且其范数足够小，从而通过 Neumann 级数保证 $ L^{(0)}_{\nu,\eta,\varepsilon} $ 的谱仍位于左半平面。此步骤的精妙之处在于，它绕开了对 $ L_{\nu,\eta,\varepsilon} $ 整体进行复杂估计的困难，转而利用解耦后各部分的结构性质分别处理。

4. 与之前工作的比较

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

1. 与 Navier-Stokes 框架 [15] 的比较：当 $ b_1 = b_2 = 0 $，系统退化为二维 Navier-Stokes 方程，此时 $ L_{\nu,\eta,\varepsilon} $ 解耦为两个独立的 Orr-Sommerfeld 型算子。定理 1.1(ii) 的条件 (1.10) 退化为 $ \|\partial_y^{-1} U\|_{L^2} > \nu $，与 [15] 中的临界条件完全一致，验证了本框架的兼容性。然而，当 $ b_1 \neq 0 $，MHD 系统展现出根本不同的物理机制：不稳定性仅在电阻率主导区域 $ \eta > \nu $ 出现（定理 1.1(i)），而在粘性主导区域 $ \nu \ge \eta $（定理 1.1(iii)）则全局稳定。这与经典水动力学中“粘性总是起稳定作用”的直觉相反，凸显了磁场 $ b_1 $ 作为耦合项引入的新自由度——它在 $ \eta > \nu $ 时能放大扰动，而在 $ \nu \ge \eta $ 时被强粘性抑制。

2. 与理想 MHD 及其他 MHD 稳定性工作的比较：Tatsuno & Dorland [47] 的数值模拟指出，理想 MHD 中无需拐点即可发生不稳定性，但缺乏严格的谱分析。本文在非理想（$ \nu, \eta > 0 $）设定下，首次给出了长波不稳定的显式、可计算的谱判据（(1.9)–(1.10)），并精确量化了 $ \nu, \eta, b_1, b_2 $ 的竞争关系。相比 Dolce [21]、Knobel & Zillinger [31] 等关于 Couette 流的 Sobolev 稳定性结果，本文聚焦于更一般的周期剪切流 $ U(y) $，且不依赖于小初值假设，而是直接分析线性算子的谱结构，适用范围更广但代价是仅得到线性结论。

3. 与长波不稳定性经典理论的比较：Kolmogorov [1] 和 Meshalkin–Sinai [41] 的工作基于特定剪切流 $ U(y) = \sin y $ 进行形式分析，而 Yudovich [48] 的启发式推导未给出严格证明。本文将 [15] 的严格方法成功移植到 MHD，其推进在于：(a) 将单一方程推广至耦合方程组；(b) 引入并精确处理了磁场分量 $ b_2 $ 对 $ D_{KK} $ 谱的影响（见 (1.22) 式中的 $ b_2 \partial_y $ 交叉项）；(c) 揭示了 $ b_1 $ 分量通过 $ \varepsilon $-尺度耦合项 $ i\varepsilon b_1 $ 控制零模不稳定性阈值的全新机制。代价是：分析高度依赖于周期边界和长波假设 $ \alpha \ll 1 $，目前尚无法推广至短波或无界域情形。

5. 时间离散梯度流的均匀化

Homogenization of Time-Discrete Gradient Flows

作者: Yunhong Deng, Yuan Gao, Li Wang

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

发布日期: 2026-07-28

arXiv ID: 2607.25206v1

摘要翻译: 在本工作中，我们在时间离散层面上研究平坦环面 $\mathbb{T}^d$ 上具有快速振荡扩散系数的 Fokker–Planck 方程的均匀化极限。我们采用极小化运动格式（minimizing movement scheme）对时间进行离散，并比较两种度量选择：一种是传输型度量（transport-type metric），另一种是加权 $L^2$ 度量（weighted $L^2$ metric）。我们证明该度量选择至关重要：当尺度参数趋于零时，这两种格式导致不同的均匀化极限；进而，当时间步长趋于零时，它们导出不同的演化方程。特别地，仅基于加权 $L^2$ 度量的格式能够恢复正确的有效动力学。

Comment: 无

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

1. 研究对象和背景

本文研究定义在 d 维平环面 $\mathbb{T}^d := \mathbb{R}^d / \mathbb{Z}^d$ 上、具有快速振荡扩散系数的 Fokker–Planck 方程在时间离散层面的均匀化极限。核心对象是如下含小参数 $\varepsilon > 0$ 的非自治概率密度演化方程：

其中 $A : \mathbb{R}^d \to \mathbb{R}^{d \times d}$ 是光滑、1-周期、一致椭圆的对称正定矩阵场，满足 $\Lambda_0 I_{d\times d} \preceq A(x) \preceq \Lambda_1 I_{d\times d}$；$V^\varepsilon \in C^1(\mathbb{T}^d)$ 是有界、一致 Lipschitz 的外势，且 $\sup_{x\in\mathbb{T}^d} |V^\varepsilon(x) - V(x)| \to 0$（$\varepsilon \to 0$），$V \in C^1(\mathbb{T}^d)$ 为极限势。该方程描述一个在周期微结构中受外势驱动的扩散-漂移过程，其稳态解为 $\pi^\varepsilon(x) = e^{-V^\varepsilon(x)}$ 归一化的 Gibbs 测度。经典均匀化理论（见 [1], [20, Appendix A]）指出：当 $\varepsilon \to 0$ 时，强解 $\mu^\varepsilon$ 在 $L^2$ 拓扑下收敛于有效方程

的解，其中有效扩散张量 $A^* \in \mathbb{R}^{d\times d}$ 由周期校正器 $\chi = (\chi_1,\dots,\chi_d) \in C^2(\mathbb{T}^d;\mathbb{R}^d)$ 给出：

且 $\chi_k$ 满足 $-\nabla \cdot (A(x) \nabla \chi_k) = \nabla \cdot A_{\cdot k}(x)$（模常数）。该极限的物理意义在于：宏观尺度上，粒子运动等效于在均匀介质 $A^*$ 中遵循相同势 $V$ 的扩散-漂移。然而，这一结论建立在连续时间框架下，其本质依赖于空间均匀化与时间演化算子的可交换性——即先对每个固定时刻求解椭圆型均匀化问题，再将结果代入时间导数。本文关注的根本问题是：当采用时间离散化（特别是基于变分原理的最小移动格式）逼近原方程时，该可交换性是否仍成立？换言之，若先取 $\varepsilon \to 0$（均匀化）再取 $\tau \to 0$（时间连续极限），与先取 $\tau \to 0$ 再取 $\varepsilon \to 0$，所得极限是否一致？现有文献（如 Gao–Yip [20]）已证实：在连续时间 Wasserstein 梯度流框架下，$\varepsilon$-Wasserstein 距离 $W^\varepsilon$ 的均匀化极限并非 $W^*$，而是 Gromov–Hausdorff 极限 $W^{\mathrm{GH}}$，其诱导的成本函数 $c^{\mathrm{hom}}(x,y)$ 来源于 Riemann 度量 $\langle A(x/\varepsilon)^{-1} v, v \rangle$ 的均匀化，即

其中 $A^{\mathrm{hom}}(\xi) = \liminf_{T\to\infty} \inf_{\phi \in H^1_0([0,T];\mathbb{R}^d)} \frac{1}{T} \int_0^T \left\langle A(t\xi + \phi(t))^{-1} (\xi + \dot{\phi}(t)), \xi + \dot{\phi}(t) \right\rangle \, dt$。这揭示了均匀化操作与最优传输几何之间的深刻张力：$W^\varepsilon$ 的极限不反映有效扩散张量 $A^*$，而反映 $A^{\mathrm{hom}}$。本文旨在系统检验两种主流时间离散化方案对此张力的响应：一种是将方程重写为加权 $L^2$ 梯度流（对应热方程型结构），另一种是标准 JKO 格式（对应 Wasserstein 梯度流）。关键科学空白在于：缺乏对“时间离散化是否破坏均匀化一致性”这一根本问题的严格判别准则；尤其缺乏对不同度量选择如何导致不同有效动力学的机制性解释。本文填补此空白，证明度量选择绝非技术细节，而是决定宏观有效模型本质的结构性因素。

2. 主要定理或主要结果

本文建立了两个时间离散梯度流格式在 $\varepsilon \to 0$ 和 $\tau \to 0$ 双重极限下的精确渐近行为。主要结果分为两部分：

定理 1.4（加权 $L^2$ 格式）：设 $\mu^{\varepsilon,\tau}$ 为定义 1.2 所述加权 $L^2$ 格式的解，初始值 $\mu^0 \in H^1(\mathbb{T}^d)$。则对任意固定 $t > 0$，有 $\mu^{\varepsilon,\tau}(t) \to \mu^\tau(t)$ 在 $L^2(\mathbb{T}^d)$ 中强收敛，其中 $\mu^\tau$ 是同构于 $\mu^{\varepsilon,\tau}$ 的分段常数函数 $\mu^\tau(t) = \mu^n$（$(n-1)\tau < t \leq n\tau$），且序列 $\{\mu^n\}_{n\geq 0}$ 由以下迭代极小化问题定义：

其中 $Q(\mu) = \frac{1}{2} \int_{\mathbb{T}^d} \left\langle \pi^*(x) A^* \nabla \varrho, \nabla \varrho \right\rangle dx$，$K(\mu, \nu) = \left( \int_{\mathbb{T}^d} \frac{(\mu - \nu)^2}{\pi^*} dx \right)^{1/2}$，$\varrho = \mu / \pi^*$，$\pi^* = e^{-V}$，$A^*$ 为公式 (1.3) 定义的有效张量。

定理 1.5（$\varepsilon$-JKO 格式）：设 $\mu^{\varepsilon,\tau}$ 为定义 1.3 所述 $\varepsilon$-JKO 格式的解，初始值 $\mu^0 \in \mathcal{P}(\mathbb{T}^d)$ 且熵有限。则对任意固定 $t > 0$，有 $\mu^{\varepsilon,\tau}(t) \to \mu^\tau(t)$ 在 2-Wasserstein 距离 $W_2$ 下收敛，其中 $\mu^\tau$ 是分段常数函数 $\mu^\tau(t) = \mu^n$，且 $\{\mu^n\}_{n\geq 0}$ 由以下迭代极小化问题定义：

其中 $H(\mu) = S(\mu) + \int_{\mathbb{T}^d} V(x) d\mu$，$S(\mu)$ 为对数熵，$W^{\mathrm{GH}}$ 为由成本函数 $c^{\mathrm{hom}}$ 诱导的 Wasserstein 距离（公式 (1.20)），其对应的均匀化作用 $A^{\mathrm{hom}}$ 满足 $\Lambda_1^{-1} |\xi|^2 \leq A^{\mathrm{hom}}(\xi) \leq \Lambda_0^{-1} |\xi|^2$（公式 (1.22)）。

定理 1.7（加权 $L^2$ 格式的连续极限）：设 $\mu^\tau$ 如定理 1.4 所定义，则序列 $\{\mu^n\}$ 在 $H^1$ 弱意义下满足隐式欧拉格式：

且当 $\tau \to 0$ 时，$\mu^\tau(t) \to \mu(t)$ 在 $L^2(\mathbb{T}^d)$ 中强收敛，其中 $\mu(t)$ 是有效 Fokker–Planck 方程 (1.2) 的解，初值为 $\mu^0$。

定理 1.8（$\varepsilon$-JKO 格式的连续极限）：若 $A^{\mathrm{hom}}$ 是一致强凸的，且 $\mu^\tau$ 如定理 1.5 所定义，则存在 $\mu(t)$ 使得 $\mu^\tau(t) \to \mu(t)$ 在 $W_2$ 中收敛。进一步，若 $\mu^\tau(t) \to \mu(t)$ 在 $W^{1,1}(\mathbb{T}^d)$ 中强收敛（此为额外紧性假设），则 $\mu(t)$ 满足如下拟线性抛物方程：

特别地，在 $d=1$ 情形下，文中明确计算出：$\left( \frac{1}{2} \nabla A^{\mathrm{hom}} \right)^{-1}(\xi) = \left( \int_0^1 A(x)^{-1/2} dx \right)^{-2} \xi$，而 $A^* = \left( \int_0^1 A(x)^{-1} dx \right)^{-1}$，二者一般不等。

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

本文的核心障碍在于：如何严格分离并量化两种不同度量结构对均匀化极限的独立贡献。作者的真正创新点并非引入新工具，而是对经典工具进行精微重构与组合，构建了针对“双参数极限非交换性”的全新分析框架。

第一，加权 $L^2$ 格式的均匀化：基于 $\Gamma$-收敛的逐点能量控制。证明的关键在于将极小化问题 (1.8) 的能量泛函 $\mathcal{E}_\varepsilon(\mu) = Q_\varepsilon(\mu) + \frac{1}{2\tau} K_\varepsilon(\mu, \nu)^2$ 分解为“主导项”与“扰动项”。作者首先利用 Assumption 1.1 中 $A$ 的周期性与椭圆性，结合标准周期均匀化理论，证明加权 Dirichlet 能量 $Q_\varepsilon(\mu)$ 的 $\Gamma$-极限为 $Q(\mu)$（定理 3.1）。难点在于处理距离项 $K_\varepsilon(\mu, \nu)^2 = \int_{\mathbb{T}^d} \frac{(\mu - \nu)^2}{\pi^\varepsilon} dx$：由于 $\pi^\varepsilon = e^{-V^\varepsilon}$ 本身依赖于 $\varepsilon$，直接应用 $\Gamma$-收敛需克服 $L^\infty$ 收敛性（$V^\varepsilon \to V$）与 $L^2$ 拓扑的耦合。作者通过构造显式恢复序列 $\mu_\varepsilon = \mu$ 并利用 $\|\pi^\varepsilon - \pi^*\|_{L^\infty} \to 0$，证明 $K_\varepsilon(\cdot, \nu)^2$ 连续收敛于 $K(\cdot, \nu)^2$（引理 3.3）。最终，依据引理 2.13（$\Gamma$-收敛在连续扰动下的稳定性），得到 $\mathcal{E}_\varepsilon \xrightarrow{\Gamma} \mathcal{E} = Q + \frac{1}{2\tau} K(\cdot, \nu)^2$。此框架的威力在于：它完全规避了对 Wasserstein 几何的依赖，仅需 $H^1$ 空间上的经典椭圆估计，从而自然继承了 $A^*$ 的结构。

第二，$\varepsilon$-JKO 格式的均匀化：基于最优传输的路径空间均匀化。此处障碍截然不同：$W^\varepsilon$ 的定义涉及所有概率测度间的耦合，其均匀化本质上是对无限维路径空间 $H^1([0,1];\mathbb{T}^d)$ 上作用泛函 $\int_0^1 \langle A(x/\varepsilon)^{-1} \dot{z}, \dot{z} \rangle dt$ 的均匀化。作者没有直接处理 Monge-Kantorovich 对偶，而是回归到 Benamou-Brenier 公式，将 $W^\varepsilon$ 表示为速度场 $v_t$ 的极小化：$W^\varepsilon(\mu, \nu)^2 = \min_{\partial_t \mu_t + \nabla \cdot (\mu_t v_t) = 0} \int_0^1 \int_{\mathbb{T}^d} \langle A(x/\varepsilon)^{-1} v_t, v_t \rangle d\mu_t dt$。证明的核心是引理 4.4：对于任意固定曲线 $(\mu_t, v_t)$，其能量 $\int_0^1 \int_{\mathbb{T}^d} \langle A(x/\varepsilon)^{-1} v_t, v_t \rangle d\mu_t dt$ 的 $\varepsilon \to 0$ 极限等于 $\int_0^1 \int_{\mathbb{T}^d} A^{\mathrm{hom}}(v_t) d\mu_t dt$。这依赖于 Riemann 度量均匀化的经典结果（[9, Theorem 2.6]），但作者将其嵌入到测度流框架中，并利用 Lemma 2.6（Wasserstein 距离的等价性）和 Prokhorov 紧性（Lemma 2.10）来保证极小化序列的收敛性。最终，通过 $\Gamma$-收敛的 Urysohn 性质（Lemma 2.14），证明自由能 $H_\varepsilon$ 连续收敛于 $H$，而 $W^\varepsilon$ 的平方 $\Gamma$-收敛于 $W^{\mathrm{GH}}$ 的平方，从而得到定理 1.5。

第三，连续极限的识别：隐式欧拉格式的弱形式解析。对于加权 $L^2$ 格式，作者将极小化条件 (1.8) 转化为变分不等式：对任意测试函数 $\phi \in H^1$，有 $\int_{\mathbb{T}^d} \left\langle \pi^\varepsilon A(x/\varepsilon) \nabla \varrho^\varepsilon_n, \nabla \phi \right\rangle dx + \frac{1}{\tau} \int_{\mathbb{T}^d} \frac{(\mu^\varepsilon_n - \mu^\varepsilon_{n-1})}{\pi^\varepsilon} \phi dx = 0$。再利用 $\Gamma$-极限的性质，将此不等式在 $\varepsilon \to 0$ 后取极限，直接得到定理 1.7 中的隐式欧拉格式。而对于 $\varepsilon$-JKO 格式，作者指出其连续极限的识别高度依赖于 $\mu^\tau$ 在 $W^{1,1}$ 中的紧性（Remark 1.9），这反映了 Wasserstein 梯度流在高正则性要求上的内在困难，也解释了为何其极限方程 (1.23) 的系数 $\left( \frac{1}{2} \nabla A^{\mathrm{hom}} \right)^{-1}$ 无法简化为 $A^*$。

4. 与之前工作的比较

本文与既有工作的根本性推进体现在三个维度：

1. 与 Gao–Yip [20] 的比较：从连续到离散的范式跃迁。[20] 首次将 Fokker–Planck 方程的均匀化置于 Wasserstein 梯度流框架下，证明了连续时间极限的一致性：$\varepsilon$-Wasserstein 梯度流收敛于以 $W^*$ 为度量的有效梯度流。本文则揭示了该结论在离散层面的脆弱性。关键差异在于：[20] 的证明依赖于能量耗散不等式 (EDI) 和梯度流的 $\Gamma$-收敛理论，其核心是“先均匀化度量，再构造梯度流”。而本文证明 $\varepsilon$-JKO 格式却走向反面：它强制执行“先构造离散梯度流，再均匀化”，导致度量 $W^\varepsilon$ 的均匀化被降维至路径空间，产生 $W^{\mathrm{GH}}$ 而非 $W^*$。因此，本文不是 [20] 的简单推广，而是对其结论适用边界的严格刻画——它证明了连续时间框架的优雅一致性在时间离散化下会崩塌，除非主动选择与均匀化兼容的度量。

2. 与标准 JKO 文献（如 [24, 31, 12]）的比较：度量选择的结构性后果。经典 JKO 理论（如 Jordan–Kinderlehrer–Otto [24]）确立了 Fokker–Planck 方程作为 $W_2$ 梯度流的普适性，其核心是度量与扩散算子的内在匹配。本文颠覆了这一认知：当扩散系数快速振荡时，“匹配”不再是唯一选项。本文明确指出，$\varepsilon$-JKO 格式所依赖的 $W^\varepsilon$ 度量，其几何本质是粒子轨迹的最短路径长度，而均匀化后的 $W^{\mathrm{GH}}$ 反映的是宏观平均路径成本 $A^{\mathrm{hom}}$；而加权 $L^2$ 格式所依赖的 $K_\varepsilon$ 度量，其几何本质是加权 $L^2$ 距离，其均匀化自然继承 $A^*$ 的张量结构。因此，本文首次为“何种度量能保持均匀化一致性”提供了可判定的数学准则：当且仅当度量的定义不显式依赖于 $A(x/\varepsilon)$ 的逆（即不涉及最优传输的几何），其均匀化极限才与 $A^*$ 一致。

3. 与数值分析文献（如 [18, 21]）的比较：非交换性的统一解释。[18, 21] 在有限体积法中观察到空间与时间参数的非交换性，但未给出其与变分结构的关联。本文则从第一性原理出发，将非交换性归因于两种不同的均匀化机制：加权 $L^2$ 格式对应“空间均匀化”，即对每个固定时间点的椭圆问题进行均匀化；而 $\varepsilon$-JKO 格式对应“路径均匀化”，即对连接两点的所有可能路径进行均匀化。这种机制差异导致了根本性的结果差异：前者产生 $A^*$（源自椭圆校正器），后者产生 $A^{\mathrm{hom}}$（源自 Riemann 度量均匀化）。文中 Remark 4.12 明确指出，这正是 $\lim_{\varepsilon \to 0} \lim_{\tau \to 0}$ 与 $\lim_{\tau \to 0} \lim_{\varepsilon \to 0}$ 不可交换的深层原因。代价是：加权 $L^2$ 格式虽具有一致性，但其物理诠释较弱（不直接对应最优传输）；而 $\varepsilon$-JKO 格式虽具物理直观性，却牺牲了与宏观有效模型的一致性。从提供内容中无法确认本文是否讨论了更高阶精度或自适应时间步长的推广。

6. 退化可压缩Navier-Stokes方程在扩张区域上的无粘不可压极限

Inviscid Incompressible Limit for Degenerate Compressible Navier-Stokes Equations on Expanding Domains

作者: Tong Tang, Emil Wiedemann, Lu Zhu

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.25028v1

摘要翻译: 我们同时研究三维退化可压缩Navier-Stokes方程在扩张区域上的无粘极限与低马赫数极限，其中粘性系数依赖于密度。从不适定初值出发，我们证明极限系统为不可压Euler系统。

Comment: 18 pages, 1 figure

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

1. 研究对象和背景

本文研究三维退化可压缩Navier-Stokes方程在随参数ε→0而不断扩张的有界区域Ω_M ⊂ ℝ³上的双重奇异极限：粘性系数ε^α ρ（即剪切粘性为ρ，体积粘性为0）趋于零，同时Mach数ε→0，且空间域Ω_M满足|∂Ω_M|₂ ≤ β M²与{x ∈ ℝ³ : |x| < ω M} ⊂ Ω_M，即Ω_M以速率M线性扩张并最终覆盖全空间ℝ³。所研究的方程组为：

其中密度ρ ≥ 0，速度u ∈ ℝ³，压力状态方程取p(ρ) = ρ^γ（γ > 1），粘性张量D_x(u) = \frac{1}{2}(∇x u + ∇x u^\top)，r₁(ε) > 0且r₁(ε) → 0，α > 0为固定参数。边界条件为ρu|{∂Ω_M} = 0与∇x ρ × n|{∂Ω_M} = 0，对应无滑移-零通量混合型物理边界。初始数据为非适定（ill-prepared）：ρ|{t=0} = ρ_{0,ε} = 1 + ε ρ^{(1)}{0,ε}，u|{t=0} = u_{0,ε}，满足ρ^{(1)}_{0,ε} → ρ^{(1)}0 ∈ L²(ℝ³)，u{0,ε} → u_0 ∈ L²(ℝ³)，但u₀未必满足div u₀ = 0——这导致初值含强声波振荡分量，是低Mach极限的核心障碍。

该模型属于“退化”情形：因粘性系数ρ在真空（ρ = 0）处消失，动量方程在ρ ≈ 0区域失去椭圆正则性，∇_x u的信息严重缺失；经典Lions–Feireisl理论（要求常数粘性、γ > 3/2）失效，需依赖Bresch–Desjardins熵结构（如[3,4,24]）构造弱解。本文中，作者采用Lacroix-Violet–Vasseur型重写：div_x(ρ D_x(u)) = div_x(√ρ S_μ)，其中S_μ定义为√ρ T_μ = ∇x(ρu) − √ρ u ⊗ √ρ ∇x ln ρ，从而将粘性项转化为L²可控形式（见(2.5)中∥S_μ∥{L²{t,x}} ≤ C ε^{-α/2}）。此重写对后续相对能量估计至关重要。

问题的物理与数学重要性在于：它统一刻画了三重尺度分离——小粘性（ε^α）、弱可压缩性（ε²作为压力缩放因子）、大尺度几何（Ω_M ∼ B_{ωM}）——如何协同驱动流体从可压缩粘性行为向理想不可压无粘动力学过渡。若仅考虑常数粘性，Feireisl等[12]已建立ε→0下向Euler方程的弱强收敛；但密度依赖粘性引入根本新困难：标准相对熵方法中，粘性耗散项ε^α ∫ |∇_x u|² dx在ρ→0处无法控制，传统能量估计失效。此外，扩张区域设置虽规避了边界层（因声波传播速度有限且Ω_M扩张快于声速ε⁻¹），却使目标Euler解v(t,x)无法直接作为测试函数（因v定义于ℝ³，而积分域为Ω_M），必须构造精巧校正项w_M = −η_M v − η_M ∇_x Φ以弥合域不匹配。因此，本文试图填补的关键空白是：在退化粘性框架下，首次实现无粘极限与低Mach极限的同步严格耦合，并克服扩张域带来的非紧致性与粘性退化导致的相对能量估计失稳。

2. 主要定理或主要结果

定理2.1（主收敛结果）：设α > 0，初始数据满足ρ^{(1)}_0, v_0 = H[u_0] ∈ C^m(ℝ³)（m > 4）且具有紧支集，T_max > 0为Euler系统∂_t v + v·∇_x v + ∇_x Π = 0, div_x v = 0在ℝ³上光滑解的最大存在时间。对任意0 < T < T_max及满足(H1)–(H4)的扩张域族{Ω_M}，若(ρ_ε,u_ε)是(1.1)的弱解，则当ε→0, M→∞且εM(ε)→∞时，有：

1. 速度强收敛：√ρ_ε u_ε → v 强收敛于L²(0,T; L²(K))，对任意紧集K ⊂ ℝ³；
2. 密度强收敛：ρ_ε → 1 强收敛于L^∞(0,T; L^{min{2,γ}}(K))；
3. 收敛速率由相对能量控制：定义相对能量泛函
   $$ E(ρ,u|r,U) = \int_{Ω_M} \left[ \frac{1}{2} ρ |u−U|^2 + \frac{1}{ε^2} \big( H(ρ) − H'(r)(ρ−r) − H(r) \big) \right] dx, $$
   其中H(ρ) = ∫₁^ρ p(z)/z² dz，取参考态r = 1 + ε s，U = v + ∇_x Φ + w_M（Φ为声势，w_M为截断校正项），则存在常数C > 0使得
   $$ \sup_{t∈[0,T]} E(ρ_ε,u_ε|1+ε s,v+∇_xΦ+w_M)(t) + ε^α \int_0^T \int_{Ω_M} |S_μ|^2 dx dt ≤ C \left( ε^2 + \frac{1}{M^2} \right), $$
   其中误差项ε²来自声波Strichartz衰减（见(2.12)中τ/ε因子），1/M²来自w_M的W^{2,p}估计（见(2.13)）。特别地，当α ≥ 1时，粘性耗散项主导，收敛无需额外小ness假设；当0 < α < 1时，需ε^α ≫ ε²即ε^{2−α} → 0以保证主导误差项为ε²。

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

核心障碍有三：(i) 退化粘性导致标准∇x u估计失效，无法直接控制ε^α ∫ ρ |D_x(u)|² dx；(ii) 扩张域使Euler解v不能作为测试函数，需处理∫{Ω_M} […] dx − ∫_{ℝ³} […] dx的几何误差；(iii) ill-prepared初值激发O(1)声波振荡s, ∇_x Φ，其能量不衰减但高频振荡需被相对能量吸收。

创新点一：采用纯相对能量框架，完全规避Bresch–Desjardins结构依赖。作者将粘性项重写为ε^α ∫ √ρ S_μ : ∇_x U dx，并利用S_μ ∈ L²（由(2.5)保证）及U = v + ∇_x Φ + w_M的正则性，通过Cauchy–Schwarz得

此估计仅依赖S_μ的L²有界性，不涉及∇_x u本身，是处理退化粘性的关键突破。

创新点二：构造双尺度校正项w_M = −η_M v − η_M ∇x Φ。其中η_M满足η_M|{∂Ω_M} = 1且dist(x,∂Ω_M) > 1时η_M = 0，故w_M在∂Ω_M上精确补偿v与∇x Φ的边界贡献。由(2.13)，∥w_M∥{W^{2,p}(Ω_M)} ≤ C M^{2(1/p−1)}，特别取p = 6/5得∥∇x² w_M∥{L^{6/5}} ≤ C M^{−1/3}，此衰减率与Ω_M体积增长M³平衡后，使w_M相关项R₄,R₅的贡献为O(M⁻¹)（见(3.3)中1/ε² ∫ (r−ρ)H″(r) div_x w_M dx）。

创新点三：声波-涡量解耦与Strichartz估计嵌入。将U分解为v（涡量主导）+ ∇_x Φ（声波主导）+ w_M（几何校正），并利用声波系统(2.10)的守恒律(2.11)与Strichartz估计(2.12)：对k=0, p=6, q=6/5, d>3/2，有

结合假设εM→∞，得∫₀^T ∥∇x Φ∥{L^6}² dt ≤ C ε^{2/3}，此衰减被相对能量中的1/ε²因子放大后仍为O(ε⁻⁴⁄³)，但通过选取合适s（如s = Δ⁻¹ div(v·∇_x v)）可将其吸收进主导误差ε²中。

最后，Gronwall引理应用需验证所有余项R_i满足|R_i| ≤ C E + o(1)，其中o(1)为ε² + M⁻²。例如R₂项

利用H′(r)−H″(1)(r−1)−H′(1) = O((r−1)²) = O(ε² s²)，故被1/ε²抵消后为O(s²|U−u|)，再由s∈L^∞(L⁶)与|U−u|∈L²(L³)插值得到L¹控制。

4. 与之前工作的比较

与Feireisl et al. [12]比较：后者处理常数粘性μ = const的Navier-Stokes方程，在固定有界域上证明ε→0时向Euler方程的弱强收敛，其相对熵中粘性项为ε ∫ |∇_x u|² dx，直接由能量不等式控制；而本文粘性项为ε^α ∫ ρ |D_x(u)|² dx，退化性迫使作者放弃∇_x u估计，转而发展S_μ-L²框架，这是本质方法论跃迁。代价是：[12]允许γ > 3/2，本文因需S_μ ∈ L²及H″(ρ) = (γ−1)ρ^{γ−2}有界性，隐含要求γ ≥ 2（否则H″(ρ)在ρ→0时爆炸）。

与Bisconti & Caggio [1]比较：后者在固定Mach数下证退化NS的无粘极限，未涉及Mach数趋零，故无需处理声波振荡；本文必须同时压制O(1)声波能量与O(ε⁻²)压力梯度，为此引入声势Φ并利用其Strichartz衰减，这是[1]未涉及的新技术维度。

与Chaudhuri et al. [7]比较：后者研究退化NS的anelastic极限（保留浮力效应），使用修正相对熵但依赖Bresch–Desjardins结构来控制粘性；本文明确声明“不依赖Bresch–Desjardins结构”，证明更简洁且适用范围更广（如可推广至γ > 1的更广范围），但代价是丧失对真空附近精细结构的描述能力。

与Lions [22]和Feireisl [10]的存在性理论比较：本文不重建解的存在性，而是条件性地假设对每个ε,M存在弱解（由[4]保证），聚焦于极限过程本身；因此结论是‘weak-to-strong’而非‘weak-to-weak’，要求目标Euler解v光滑，这是此类相对熵方法的固有局限，无法处理Euler解奇点形成后的激波极限。

与扩张域先驱工作[15]比较：[15]处理常数粘性热传导方程，其校正项w_M仅需补偿v，而本文需同时补偿v与∇_x Φ，导致w_M构造更复杂且误差分析中出现额外的1/ε ∫ (r−ρ)H″(r) div_x(sU) dx项，此为退化粘性与声波耦合独有现象。

7. 扩散稳定守恒律与松弛振荡器耦合系统中的时间周期解

Diffusion stabilises time-periodic solutions in conservation laws coupled to a relaxation oscillator

作者: Julien Barré, Nils Berglund, Hiroshi Horii

类别: math.AP, math.DS

发布日期: 2026-07-27

arXiv ID: 2607.24994v1

摘要翻译: 我们研究一个粘性一维守恒律，其与一个快速常微分方程耦合。当粘性系数趋于零时，该系统具有无限维的时间周期解族，对应于松弛振荡。我们证明：对称初值条件下，任意小的正粘性选择出唯一的周期解，并证明该解是线性稳定的。证明利用系统的慢-快结构，通过平均化策略以及谱理论方法完成。结果由数值模拟予以验证。

Comment: 54 pages, 13 figures

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

1. 研究对象和背景

本文研究如下一维粘性守恒律与快速常微分方程耦合的慢-快系统：

其中 $n(t,x) \in [0,1]$ 表示粒子密度，$E(t) \in \mathbb{R}$ 是全局电场类变量，$M(t)$ 是左半轴总质量（即负方向累积量），$\varepsilon > 0$ 是扩散系数，$\eta > 0$ 是快慢时间尺度分离参数，且满足 $0 < \eta \ll \varepsilon \ll 1$。通量函数取 $f(n) = n(1-n)$，保证双稳结构与激波候选性；函数 $a$ 和 $g$ 满足 Assumption 2.4 所列几何条件：$g(E,M)=0$ 定义一条 S 形临界流形 $M = h(E)$，其具有两个非退化折点 $(E_\pm, h(E_\pm))$，且 $h'(E) < 0$ 在 $(E_-,E_+)$ 内成立；$a(E,M)=0$ 定义慢零线 $E = k(M)$，二者在 $(E_-,E_+)$ 内唯一横截相交。该结构确保当 $\varepsilon = 0$ 时，$(E,M)$ 子系统呈现典型的范德波尔型松弛振荡——轨迹沿临界流形稳定支快速滑移，在折点处跳跃，再沿另一稳定支返回，形成大振幅周期轨道 $\Gamma$，其周期 $T$ 由快子系统主导，且对初值敏感度极低。

关键物理动机源于尘埃等离子体建模：$n$ 描述微米级尘埃粒子密度，$E$ 近似代表由尘埃、离子与电子共同产生的平均电场，长程耦合通过 $M$ 实现（理想化 Poisson 方程效应）。边界条件为

并假设初始数据 $n_0(x)$ 满足中心对称性 $n_0(-x) = 1 - n_0(x)$、严格单调递增性 $\partial_x n_0 > 0$ 及正则性 $n_0 - n^* \in H^2(\mathbb{R}) \cap L^1(\mathbb{R})$，其中 $n^*(x) = \frac{1}{2}(1 + \tanh x)$ 是参考渐近剖面。此时 $n(t,0) \equiv \frac{1}{2}$ 恒成立，故 $M(t) = -\int_{-\infty}^0 n(t,x)\,dx$ 直接刻画剖面陡峭程度：$M(t) \to 0$ 对应 $x=0$ 处激波；$|M(t)|$ 增大则剖面展宽。

当 $\varepsilon = 0$ 时，系统退化为纯输运方程 $\partial_t n = -a(E,M) \partial_x [n(1-n)]$，可由特征线法显式求解。定义反函数 $X(t,n)$ 满足 $X(t,n(t,x)) = x$，则其演化满足 $\partial_t X(t,n) = a(E(t),M(t))(1-2n)$，积分得

由于 $(E,M)$ 指数趋近于极限环 $\Gamma$，且 $\int_0^T a(E_p(s),M_p(s))\,ds = 0$（因 $\dot{M}_p = \frac{1}{4} a(E_p,M_p)$ 且 $M_p$ 周期），故 $\int_0^t a(\cdot)\,ds$ 有界震荡，导致 $X(t,\cdot)$ 收敛于 $T$-周期函数，从而 $n(t,x)$ 亦收敛于某周期函数。更精确地，$n(t,x) = n_0(Y(t,x))$，其中 $Y$ 满足 $y + [1-2n_0(y)] \int_0^t a(\cdot)\,ds = x$，表明每个初值剖面 $n_0$ 均生成一个独立的极限环——构成无限维周期解族，其共同周期为 $T$。此退化性源于无扩散时系统缺乏耗散机制，无法区分不同剖面形状，仅依赖全局质量 $M$ 的动力学。

本文旨在填补这一理论空白：在 $\varepsilon > 0$ 小但非零时，扩散项 $\varepsilon \partial_{xx} n$ 如何打破无限维简并，选择并稳定唯一周期解？该问题本质是研究奇异摄动下慢-快偏微分系统的周期轨道分支与稳定性，涉及输运-扩散竞争、全局耦合非局部性、以及对称性约束下的谱分析，此前未见严格处理此类守恒律-ODE 耦合周期解选择机制的工作。

2. 主要定理或主要结果

定理 2.7 给出核心结论：在 Assumptions 2.3（初值对称、单调、正则）与 2.4（S 形临界流形、横截慢零线）下，存在常数 $\varepsilon_0, c, C, M_0, M, \Delta_0, \Delta, \delta > 0$，使得当 $0 < \varepsilon < \varepsilon_0$ 且 $0 \leq \eta < c / \log(\varepsilon^{-1})$ 时，系统 (1.1)–(1.3) 存在唯一的孤立周期轨道 $(n^*(t,x), E^*(t))$，且该轨道线性渐近稳定。具体而言，对任意满足附加斜率控制 $\partial_x n_0(x) \leq C$ 及初值接近条件

的初始数据，解 $(n(t,x), E(t))$ 保持严格单调递增（即 $\partial_x n(t,x) > 0$），且指数收敛于该周期轨道：

此外，$0 < n(t,x) < 1$ 对所有 $t \geq 0, x \in \mathbb{R}$ 成立。周期 $T_\varepsilon$ 满足 $\lim_{\varepsilon \to 0} T_\varepsilon = T$，其中 $T$ 是 $\varepsilon = 0$ 时无限维周期族的公共周期。该结果表明扩散以 $O(\varepsilon)$ 的衰减率主导稳定性，且稳定性在 $H^2$ 范数下成立，而非仅 $L^2$ 或逐点意义。定理 2.1 保证全局适定性：对 Lipschitz 连续的 $f,a,g$ 及 $a$ 有界，若 $n_0 - n^* \in H^2 \cap L^1$，则解全局存在唯一，且保持 $n(t,\cdot) - n^* \in H^2 \cap L^1$。Proposition 2.2 证明对称性守恒：若 $f(n) = f(1-n)$ 且 $n_0$ 中心对称，则 $n(t,-x) = 1 - n(t,x)$ 恒成立，从而 $M(t)$ 与剖面陡峭度直接关联。

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

证明的核心障碍在于：(i) 系统兼具慢-快结构（$\eta \ll 1$）、非局部全局耦合（$M$ 依赖 $n$ 全空间积分）、及强非线性（双稳通量 $f(n)=n(1-n)$）；(ii) $\varepsilon$-小扩散项与输运项量级相当，传统边界层或匹配渐近展开失效；(iii) 周期解的谱稳定性需处理非自伴算子 $\mathcal{L}_\varepsilon$，其在 $H^2$ 上的本征值分布受慢变量 $E(t)$ 时变调制影响。作者的创新点在于构建一个严格的 averaging strategy，并将其与 spectral-theoretic 方法深度耦合。

首先，Section 4 发展了精细的 averaging procedure：将原系统在快速时间尺度 $\tau = t/\eta$ 上平均，但关键在于不直接平均 PDE，而是利用反函数 $X(t,n)$ 的演化方程 (2.7)，将 $n$ 的动力学映射到 $X$ 空间，其中 $\partial_t X = a(E,M)(1-2n)$。由于 $n$ 由 $X$ 隐式定义，且 $M$ 通过 $n$ 依赖于 $X$，作者引入平均通量 $\bar{a}(E) = \frac{1}{T} \int_0^T a(E_p(s), M_p(s))\,ds$（注意 $\bar{a}(E)$ 实际为常数，因 $\int_0^T a = 0$），并构造 averaged equation for $X$：$\partial_t X_{\text{av}} = \bar{a}(1-2n_{\text{av}}) + \mathcal{R}_\varepsilon$，其中余项 $\mathcal{R}_\varepsilon$ 含 $\varepsilon$-阶扩散修正与 $\eta$-阶快慢脱耦误差。该 averaging 不同于标准 Kolmogorov-Arnold-Moser 类方法，而是针对反函数空间设计，天然保持单调性与边界条件。

其次，Section 5 将 averaged system 的 stationary solution（即 $\partial_t X_{\text{av}} = 0$）对应于原系统的周期解，并通过线性化分析证明其稳定性。关键工具是构造一个 time-periodic linearized operator $\mathcal{L}_\varepsilon(t)$，作用于扰动 $v = n - n^*$：

其中 $v_E, v_M$ 由耦合关系导出。作者证明 $\mathcal{L}_\varepsilon(t)$ 的 Floquet 谱中，除单重零本征值（对应时间平移对称性）外，其余本征值实部均 $\leq -\delta \varepsilon$。此结论依赖两个核心估计：(1) 对称性诱导的正交分解：将 $v$ 分解为偶函数与奇函数分量，利用 $n^*(t,-x) = 1-n^*(t,x)$ 使奇部满足齐次 Dirichlet 条件，从而应用加权 Sobolev 不等式 $\|v_{\text{odd}}\|_{H^2} \leq C \|\partial_{xx} v_{\text{odd}}\|_{L^2}$；(2) 针对偶部的 energy estimate：定义能量泛函 $\mathcal{E}(t) = \frac{1}{2} \int (\partial_x v)^2 + \alpha v^2 \,dx$，选取 $\alpha > 0$ 使耗散项 $\varepsilon \int (\partial_{xx} v)^2$ 主导输运项的不稳定贡献，得到 $\frac{d}{dt} \mathcal{E} \leq -\delta \varepsilon \mathcal{E}$。最后，Section 6 通过 contraction mapping argument 证明余项 $\mathcal{R}_\varepsilon$ 不破坏存在性与稳定性：在 Banach 空间 $C^0([0,T]; H^2) \times C^0([0,T])$ 上定义 Poincaré 映射 $\mathcal{P}_\varepsilon$，证明其为压缩映射，且不动点唯一，从而确立周期解的孤立性。

4. 与之前工作的比较

本文与既有工作存在本质差异。经典输运-扩散行波研究（如 [26,30,1,20]）聚焦于 traveling front 的存在性与 $\varepsilon$-渐近稳定性，其解形式为 $n(t,x) = q(x-ct)$，稳定性由 Sturm-Liouville 算子谱决定，且不涉及 ODE 耦合。而本文解非行波、非空间周期，其周期完全由快 ODE 子系统驱动，PDE 仅提供形状调制。与 [36]（Hopf 分岔产生周期解）相比，本文不依赖分岔参数穿越临界值，而是通过 $\varepsilon > 0$ 的奇异摄动实现解的选择，且稳定性速率 $e^{-\delta \varepsilon t}$ 显式依赖扩散强度，而非固定指数。与 FitzHugh-Nagumo 类慢-快系统（[34,11,24]）对比，后者周期解源于内部反馈，而本文周期性完全由外部 ODE 强制，且耦合为非局部（$M$ 为积分），导致标准几何奇异摄动理论（如 Fenichel 定理）无法直接应用，因中心流形维度无限。

在技术层面，[29,38,19] 对周期行波的稳定性分析依赖 Bloch 波分解与 Evans 函数，适用于空间平移不变系统；本文系统因非局部耦合失去空间平移不变性，迫使作者发展基于反函数 $X(t,n)$ 的 averaging 框架，这是对标准 slow-fast PDE 方法的重要拓展。与 [13]（非粘性双守恒律的空间-时间周期解）不同，本文解的周期性纯粹由快变量 $E$ 的振荡引发，且 $\varepsilon$-扩散是稳定性的必要条件，而非仅正则化工具。数值上（Section 3），当 $\eta$ 增大至 $O(1)$ 时，系统仍保持周期性，表明 $\eta \ll \varepsilon$ 并非绝对必要，但理论证明中该假设用于控制 averaging 余项 $\mathcal{R}_\varepsilon = O(\eta^{1/2})$。局限在于：结果仅对称初值成立（非对称情形可能有多解，见 Figure 8）；未处理 $\varepsilon$ 与 $\eta$ 同阶情形（如 $\varepsilon \sim \eta$）；且未证明全局吸引域，仅给出 $H^2$-邻域内的指数吸引。从提供内容中无法确认是否推广至高维或多组分情形，亦未显式展开与等离子体文献 [31,42,43] 的定量模型比对。

8. 具有双尺度遗传粘性的Navier-Stokes方程：超临界范数膨胀与临界空间中的整体适定性

The Navier-Stokes equations with dual-scale hereditary viscosity: supercritical norm inflation and global well-posedness in critical spaces

作者: Bruno de Andrade

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.25009v1

摘要翻译: 本文研究由双尺度遗传记忆所支配的不可压流体流动的Cauchy问题，该模型为经典Navier-Stokes方程的一种粘弹性变体，可刻画反常动量输运。其非局部耗散破坏了精确的整体尺度不变性，从而要求在Hörmander类 $S^{-2}_{1,0}$ 中进行拟微分分析，其中空间梯度诱导出分数阶时间惩罚项。我们严格建立了 $L^q-L^p$ 衰减估计，并确定了临界Lebesgue阈值 $p_c = N (\frac{1+α_\infty}{1-α_\infty})$。在超临界区域 $1 < p < p_c$ 内，我们通过将频率调制的双线性流直接映射至Fourier空间，绕过了空间局域化损失；借助Bernstein不等式，我们证明了原点处的瞬时范数膨胀，并确认其内在不适定性。反之，在拓扑极限 $p \to \infty$ 下，我们证明双尺度记忆结构上阻止了经典流体在极大临界Besov空间 $\dot{B}^{-1}_{\infty, \infty}$ 中通常观测到的坍缩现象。通过在Bony的拟微分算子演算中采用非对称插值，我们证明时间光滑性压倒了高-高对流共振级联。这一精细的分析平衡将不适定性限制于Besov拓扑中不可分的高频尾部，从而在 $\dot{B}^{-κ}_{\infty, \infty}(\mathbb{R}^N)$ 中对具有高频依附性的初值建立了小初值的整体Hadamard适定性，其中 $κ= \frac{1-α_\infty}{1+α_\infty}$；该适定性区域严格宽于可分的小Besov闭包 $\dot{b}^{-κ}_{\infty, \infty}(\mathbb{R}^N)$。

Comment: 28 pages

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

1. 研究对象和背景

本文研究如下不可压流体的初值问题，其动量方程被替换为一个具有双尺度遗传记忆的非局部耗散模型，即纳维-斯托克斯方程的粘弹性变体：

其中 $g \in L^1_{\mathrm{loc}}(\mathbb{R}_+)$ 是一个双尺度容许核，其拉普拉斯变换 $\hat{g}(\lambda) = \mathcal{L}\{g\}(\lambda)$ 满足三类结构性假设（H1）–（H3）：(H1) 高频扇形正则性，即 $\hat{g}$ 在 $\mathbb{C}\setminus(-\infty,0]$ 上解析非零，且特征方程 $\lambda + |\xi|^2 \hat{g}(\lambda) = 0$ 的根满足 $|\arg \lambda| \le \pi - \theta(|\lambda|)$，其中 $\theta(|\lambda|) \to \theta_\infty > 0$ 且 $\theta_\infty < \frac{\pi \alpha_\infty}{1+\alpha_\infty}$；(H2) 高频渐近行为：存在常数 $C_1,C_2>0$，使得当 $|\lambda|$ 充分大时，$C_1 |\lambda|^{-\alpha_\infty} \le |\hat{g}(\lambda)| \le C_2 |\lambda|^{-\alpha_\infty}$，其中 $\alpha_\infty \in (0,1)$；(H3) 低频渐近行为：存在 $\alpha_0 \in [0,1)$ 及 $c_1,c_2>0$，使得当 $|\lambda|$ 充分小时，$c_1 |\lambda|^{-\alpha_0} \le |\hat{g}(\lambda)| \le c_2 |\lambda|^{-\alpha_0}$。该模型物理上刻画了聚合物熔体、胶束网络等复杂介质中异常动量输运的双阶段松弛机制——短时（高频）呈现强弹性响应（指数 $\alpha_\infty$），长时（低频）趋于牛顿型或弱粘弹性终端松弛（指数 $\alpha_0$）。其线性部分的傅里叶-拉普拉斯联合符号为 $\widehat{S}(\lambda,\xi) = (\lambda + |\xi|^2 \hat{g}(\lambda))^{-1}$，对应的线性算子 $S(t)$ 是一个伪微分算子，其符号在 Hörmander 类 $S^{-2}_{1,0}$ 中全局一致有界。关键数学难点在于：经典 Navier-Stokes 方程的精确尺度不变性 $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$ 被彻底破坏，因 $\hat{g}$ 的双尺度结构导致系统不再具有单一自相似标度；由此，线性解算子 $S(t)$ 的衰减行为在短时与长时呈现不同代数速率：短时由 $\alpha_\infty$ 主导，长时由 $\alpha_0$ 主导。具体而言，对 $\nabla S(t)P$（含 Leray-Hopf 投影的梯度算子），短时估计为 $\|\nabla S(t)P f\|_{L^p} \lesssim t^{-\tilde{\gamma}_{q,p}} \|f\|_{L^q}$，其中 $\tilde{\gamma}_{q,p} = \frac{1+\alpha_\infty}{2} N\left(\frac{1}{q}-\frac{1}{p}\right) + \frac{1+\alpha_\infty}{2}$；长时估计为 $\|\nabla S(t)P f\|_{L^p} \lesssim t^{-\tilde{\beta}_{q,p}} \|f\|_{L^q}$，其中 $\tilde{\beta}_{q,p} = \frac{1+\alpha_0}{2} N\left(\frac{1}{q}-\frac{1}{p}\right) + \frac{1+\alpha_0}{2}$。这一双重尺度依赖性直接导致传统基于单一齐次尺度的调和分析工具（如标准子序理论、Wright 函数表示）失效，迫使作者发展一套精细的拟微分-小波复合框架。研究的核心空白在于：既有工作（如 Barbu-Sritharan [3]、de Andrade et al. [10,11,13]）均局限于单尺度分数核或次临界空间；而双尺度情形下，临界指标、范数爆炸机制及整体适定性边界均未被刻画。尤其，Bourgain-Pavlović 在 $\dot{B}^{-1}_{\infty,\infty}$ 中发现的经典范数爆炸是否延续至本模型？若否，其拓扑障碍被何种新机制所克服？这构成了本文的根本动机。

2. 主要定理或主要结果

本文建立了三类严格结论：

（1）线性估计与临界阈值：在假设（H1）–（H3）下，证明了 $L^q-L^p$ 衰减估计。特别地，对梯度-投影组合算子 $\nabla S(t)P$，当 $1<q\le p\le \infty$ 且 $\frac{1}{q}-\frac{1}{p}<\frac{1}{N}$ 时，存在常数 $C_T>0$ 使得对任意 $0<t\le T$，

由此导出 Lebesgue 临界指数 $p_c = N \frac{1+\alpha_\infty}{1-\alpha_\infty}$，即当 $1<p<p_c$ 时，空间 $L^p$ 不足以吸收对流项引起的正则性损失。

（2）超临界范数爆炸：对任意 $1<p<p_c$，存在初始数据序列 $\{u_0^{(n)}\} \subset \mathcal{S}(\mathbb{R}^N)$ 满足 $\|u_0^{(n)}\|_{L^p} \to 0$，但对应 mild 解 $u^{(n)}$ 满足 $\liminf_{t\to 0^+} t^{\varepsilon} \|u^{(n)}(t)\|_{L^p} = \infty$ 对任意 $\varepsilon>0$ 成立；更精确地，利用 Christ-Colliander-Tao 频率调制技巧与 Bernstein 不等式，证明了在原点处的瞬时范数爆炸：存在 $\delta>0$ 及 $t_n \to 0^+$ 使得 $\|u^{(n)}(t_n)\|_{L^p} \ge \delta t_n^{-\sigma}$，其中 $\sigma = \frac{1}{2}\left(\frac{N}{p} - \frac{1+\alpha_\infty}{1-\alpha_\infty}\right)>0$。

（3）临界 Besov 空间中的整体适定性：令 $\kappa = \frac{1-\alpha_\infty}{1+\alpha_\infty}$，则对任意小初值 $u_0 \in \dot{B}^{-\kappa}_{\infty,\infty}(\mathbb{R}^N)$ 满足 $\|u_0\|_{\dot{B}^{-\kappa}_{\infty,\infty}} < \varepsilon_0$（$\varepsilon_0$ 依赖于 $\alpha_\infty,N$），存在唯一全局 mild 解 $u \in C([0,\infty); \dot{B}^{-\kappa}_{\infty,\infty})$，且满足 Hadamard 连续性：若 $u_0^{(n)} \to u_0$ 在 $\dot{B}^{-\kappa}_{\infty,\infty}$ 中，则对应解 $u^{(n)} \to u$ 在 $C_{\mathrm{loc}}([0,\infty); \dot{B}^{-\kappa}_{\infty,\infty})$ 中。该结果严格强于 little Besov 空间 $\dot{b}^{-\kappa}_{\infty,\infty}$ 上的适定性，因 $\dot{b}^{-\kappa}_{\infty,\infty}$ 要求 $\lim_{|j|\to\infty} 2^{-j\kappa} \|\Delta_j u_0\|_{L^\infty} = 0$，而 $\dot{B}^{-\kappa}_{\infty,\infty}$ 允许高频分量 $\Delta_j u_0$ 仅满足 $\sup_j 2^{-j\kappa} \|\Delta_j u_0\|_{L^\infty} < \infty$，即允许非可分离的高频尾部。

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

本文的核心技术障碍是双重的：其一，双尺度记忆核 $g$ 彻底破坏全局尺度不变性，致使经典自相似方法与子序理论失效；其二，Leray-Hopf 投影 $P$ 作为 Calderón-Zygmund 奇异积分，在 $L^\infty$ 上无界，而临界 Besov 空间 $\dot{B}^{-\kappa}_{\infty,\infty}$ 的范数本质上是 $L^\infty$ 型的，故必须绕过 $P$ 的 $L^\infty$ 无界性。作者的真正创新点体现在三个相互支撑的分析框架中：

第一，拟微分-尺度分解耦合框架。作者将线性算子 $S(t)$ 的符号 $\widehat{S}(t,\xi)$ 在高频区 $|\xi| \gtrsim t^{-\frac{1+\alpha_\infty}{2}}$ 和低频区 $|\xi| \lesssim t^{-\frac{1+\alpha_\infty}{2}}$ 分别处理。对高频区，引入时空缩放 $\eta = \xi t^{\frac{1+\alpha_\infty}{2}}$ 和 $\mu = t^{\frac{1+\alpha_\infty}{2}}$，使缩放后符号 $\widetilde{\Psi}(t,\eta) = \widehat{S}(t,\mu^{-1}\eta)$ 一致属于 $S^{-2}_{1,0}$ 类；梯度作用 $\nabla$ 对应乘子 $i\xi = \mu^{-1} i\eta$，从而 $\nabla S(t)$ 的符号落入 $S^{-1}_{1,0}$，其卷积核在原点附近具有 $O(|w|^{-(N-1)})$ 奇异性，由此导出几何约束 $\frac{1}{q}-\frac{1}{p} < \frac{1}{N}$。此步骤将物理参数 $\alpha_\infty$ 精确嵌入到衰减指数 $\tilde{\gamma}_{q,p}$ 中，是临界阈值 $p_c$ 的根源。

第二，频率调制与 Bernstein 不等式的协同运用。为证明超临界范数爆炸，作者摒弃了依赖空间局域化的传统方法（因 $P$ 会破坏局域性），转而构造紧支集谱初值 $u_0^{(n)}$，其傅里叶变换 $\widehat{u_0^{(n)}}$ 支集于 $\{ |\xi| \sim 2^{j_n} \}$，并利用 Bernstein 不等式 $\|\Delta_j f\|_{L^p} \lesssim 2^{jN(\frac{1}{r}-\frac{1}{p})} \|\Delta_j f\|_{L^r}$（$r<p$）控制各频段。关键一步是将双线性项 $P\nabla\cdot(u\otimes u)$ 的频率调制形式 $\sum_{k,m} \Delta_k u \Delta_m u$ 直接映射至傅里叶空间，通过相位共振分析揭示 $k\sim m\sim j_n$ 的高-高交互主导能量向低频级联，从而触发 $\|u(t_n)\|_{L^p} \gtrsim 2^{j_n(\frac{N}{p} - \frac{1+\alpha_\infty}{1-\alpha_\infty})} t_n^{-1/2}$ 的爆炸。

第三，Bony 分解中的非对称插值。在临界 Besov 情形，作者对非线性项 $P\nabla\cdot(u\otimes u)$ 应用 Bony 分解：$uv = T_u v + T_v u + R(u,v)$。核心突破在于对高-高残差项 $R(u,u)$ 施加非对称插值：将空间基拓扑 $\dot{B}^{-\kappa}_{\infty,\infty}$ 与分数时间权 $t^\gamma$（$\gamma>0$）耦合，即估计 $\|\Delta_j R(u,u)\|_{L^\infty} \lesssim \sum_{k\ge j-3} 2^{-j\kappa} 2^{k\kappa} \|\Delta_k u\|_{L^\infty} \cdot 2^{-k\kappa} \|\Delta_k u\|_{L^\infty} \cdot t^{-\gamma}$。此插值避免了 $t=0$ 处的时间奇异性，并证明遗传记忆的平滑效应 $t^{-\gamma}$ 严格压倒高-高共振级联 $2^{k\kappa}$ 的增长，从而将不适定性严格限制在 Besov 拓扑的非可分高频尾部 $\{j: 2^{-j\kappa} \|\Delta_j u\|_{L^\infty} \not\to 0\}$，而非整个空间。

4. 与之前工作的比较

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

（1）与单尺度分数模型的比较：此前所有关于 NSHV 的工作（如 [3,10,11,13]）均假设 $\hat{g}(\lambda) = \lambda^{-\alpha}$ 或其有限和，此时 $\alpha_\infty = \alpha_0 = \alpha$，系统保持全局尺度不变性 $u_\lambda(x,t) = \lambda^{\alpha} u(\lambda x, \lambda^{1+\alpha} t)$，线性算子可用 Wright 函数显式表示，临界空间为 $\dot{B}^{-\frac{2}{1+\alpha}}_{\infty,\infty}$。本文则彻底放弃该假设，引入 $\alpha_\infty \neq \alpha_0$ 的双尺度结构，导致临界指数 $\kappa = \frac{1-\alpha_\infty}{1+\alpha_\infty}$ 仅由高频参数决定，而长时行为由 $\alpha_0$ 控制；这使模型能刻画真实复杂流体的终端弛豫，但代价是丧失显式解表达式与经典子序理论。

（2）与经典 Navier-Stokes 的比较：Koch-Tataru 在 $\mathrm{BMO}^{-1}$ 中建立的适定性要求 $\alpha_\infty = 0$（即 $p_c = N$），而 Bourgain-Pavlović 证明 $\dot{B}^{-1}_{\infty,\infty}$ 中的范数爆炸源于低频共振。本文表明，当 $\alpha_\infty > 0$ 时，双尺度记忆在 $\dot{B}^{-\kappa}_{\infty,\infty}$（$\kappa < 1$）中结构性地抑制了该爆炸，其机制并非增强耗散，而是通过高频弹性响应 $\alpha_\infty$ 改变线性算子的符号阶数，使 $\nabla S(t)$ 的时间衰减 $t^{-\frac{1+\alpha_\infty}{2}}$ 足够快以压制高-高级联。

（3）与次临界适定性工作的比较：文献 [11] 在 Besov-Morrey 空间 $\dot{\mathcal{N}}^{s}_{p,q,\lambda}$ 中建立次临界适定性，其条件 $s > -1 + \frac{N}{p}$ 本质是 $L^p$ 型；本文首次将适定性推进至真正的临界空间 $\dot{B}^{-\kappa}_{\infty,\infty}$，且证明其范围严格大于 little Besov 空间 $\dot{b}^{-\kappa}_{\infty,\infty}$，这是首个在非局部流体模型中实现的非可分离临界适定性。

（4）与范数爆炸研究的比较：Christ-Colliander-Tao [8] 在 $L^p$（$p<N$）中对经典 NS 构造范数爆炸，但依赖空间局域化；本文的爆炸构造完全在傅里叶空间进行，不依赖 $L^p$ 局部性质，因而适用于任何破坏空间局域性的算子（如 $P$）。然而，从提供内容中无法确认本文是否改进了爆炸速率的最优性，亦未显式展开与 Cheskidov-Shvydkoy [7] 几何化爆炸机制的比较。

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

1. 与平面上双抛物型Keller-Segel方程相关的两物种粒子系统的渐近行为

Asymptotics of a two-species particle system associated to the doubly parabolic Keller-Segel equation in the plane

作者: Nicolas Fournier, Milica Tomašević

类别: math.PR, math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25986v1

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

摘要翻译

中文标题: 与平面上双抛物型Keller-Segel方程相关的两物种粒子系统的渐近行为

摘要翻译: 我们考虑Stevens (2000)提出的与双抛物型Keller-Segel方程相关的两物种粒子系统。该系统由 $N$ 个细胞及数量可变的趋化因子粒子组成。细胞在平面上扩散，并沿趋化因子浓度的经验梯度（经磨光后）运动；趋化因子粒子由细胞以某一常速率产生，同时以另一常速率扩散并衰减。我们证明：当细胞对趋化因子的敏感性足够小时，在磨光子族满足某较弱条件下，该系统在 $N\to \infty$ 时逼近抛物-抛物型Keller-Segel方程。此外，我们还证明：当 $N$ 固定时，若趋化因子粒子的产生速率趋于无穷，则该系统逼近Talay-Tomašević (2020)提出的（非马尔可夫型）单物种系统，该系统随后被作者们进一步研究 (2023)。

2. 无限柱形区域上分数阶拉普拉斯算子的大解

Large Solutions for Fractional Laplacian on Infinite Cylindrical Domains

作者: Indranil Chowdhury, N N Dattatreya

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25906v1

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

摘要翻译

中文标题: 无限柱形区域上分数阶拉普拉斯算子的大解

摘要翻译: 我们研究定义在某些方向（但非全部方向）趋于无界的区域上的线性与半线性方程的大解，其中涉及分数阶拉普拉斯算子 $(-\Delta)^s$。在区域边界上爆破的解通常被称为大解。对于局部算子，此类行为仅在存在低阶非线性项时出现；然而，对于非局部算子，此类解的存在性更为微妙。与局部情形不同，算子的非局部性导致即使在线性方程情形下也产生不同的边界爆破现象。本文研究有限柱形区域上分数阶线性与半线性方程的边界爆破解的存在性与定性行为，并利用这些性质构造无限柱体上的大解。

3. 乘积紧李群上的双奇异 $(ρ,δ)$ 伪微分算子

Bisingular $(ρ,δ)$ pseudodifferential operators on products of compact Lie groups

作者: G. Drei, S. Federico, A. Parmeggiani

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25821v1

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

摘要翻译

中文标题: 乘积紧李群上的双奇异 $(ρ,δ)$ 伪微分算子

摘要翻译: 本文发展了乘积紧李群上的 $(ρ,δ)$ 型双奇异伪微分算子演算，该演算推广了第二与第三作者早先在 $(1,0)$ 情形下引入的演算。$(ρ,δ)$ 情形是适用于在 $G_1\times G_2$ 上构造 hypoelliptic 算子的基本解（parametrices）的演算。

4. 与Minkowski泛函相关的尖锐各向异性 $L^2$-Caffarelli-Kohn-Nirenberg 不等式

Sharp anisotropic $L^2$-Caffarelli-Kohn-Nirenberg inequalities associated with the Minkowski functional

作者: Zhenzhen Wei

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25876v1

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

摘要翻译

中文标题: 与Minkowski泛函相关的尖锐各向异性 $L^2$-Caffarelli-Kohn-Nirenberg 不等式

摘要翻译: 设 $K\subset \mathbb{R}^N$ 为一个包含原点在其内部的凸体，记其 Minkowski 泛函为 $\hK{\cdot}$。本文建立了一个基于恒等式的框架，用于证明与各向异性径向导数

相关的尖锐各向异性 $L^2$-Caffarelli-Kohn-Nirenberg 不等式。本工作的关键在于：未假设 $K$ 关于原点对称。因此，Minkowski 泛函 $\hK{\cdot}$ 不必为偶函数，标准的基于范数的各向异性论证无法直接应用。主要工具是带有显式非负余项的各向异性 $L^2$-Hardy 恒等式与各向异性 $L^2$-Caffarelli-Kohn-Nirenberg 恒等式。这些恒等式导出尖锐的各向异性 $L^2$-Caffarelli-Kohn-Nirenberg 不等式，其最优常数依赖于参数区域 $(a,b)\in\mathbb{R}^2$。我们还研究了该最优常数在自然完备空间中的可达性，并得到了相应的极值函数。作为进一步推论，我们导出了尖锐的各向异性 Heisenberg 型不确定性原理与极大值型各向异性梯度不等式。当 $K$ 为欧氏单位球时，我们的结果恢复经典欧氏 $L^2$ 理论；当 $K$ 关于原点对称时，结果与标准的基于范数的各向异性框架一致。特别地，本文结果将尖锐的 $L^2$-Caffarelli-Kohn-Nirenberg 理论推广至任意包含原点在其内部的凸体——此时 Minkowski 泛函可能非偶。

5. 有限质量孤子型刚性与三维非径向聚焦能量临界非线性薛定谔方程的四通道约化

Finite-mass soliton-type rigidity and four-channel reduction for the three-dimensional nonradial focusing energy-critical nonlinear Schrödinger equation

作者: Pang-Hung Chung, Dan Han

类别: math.AP, math.DS

发布日期: 2026-07-28

arXiv ID: 2607.25735v1

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

摘要翻译

中文标题: 有限质量孤子型刚性与三维非径向聚焦能量临界非线性薛定谔方程的四通道约化

摘要翻译: 本文研究三维非径向聚焦能量临界非线性薛定谔方程在基态阈值以下的集中-紧致性通道。在一侧归一化后，极小临界元落入四类通道：有限时间通道、快速级联通道、有界尺度有限质量通道以及剩余拟孤子通道。前三个通道被严格排除。主要结果表明：在无需径向对称性、零动量或缓慢空间中心的假设下，每一个有限质量、有界尺度的几乎周期解恒等于零。因此，任何低于阈值散射的极小反例必位于剩余拟孤子通道；若其尺度有界，则其在每一时刻具有无限质量。

6. 带非自治强迫项的紧李群上的半线性阻尼波动方程

Semilinear damped wave equation on a compact Lie group with a non-autonomous forcing term

作者: Wenhui Chen, Sandra Lucente, Alessandro Palmieri

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25571v1

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

摘要翻译

中文标题: 带非自治强迫项的紧李群上的半线性阻尼波动方程

摘要翻译: 在本注记中，我们考虑带非自治非线性项 $\varphi(t)|u|^p$ 的紧李群上的半线性阻尼波动方程。我们关注非负时变因子 $\varphi$ 如何影响局部解的整体延拓性（global in time prolongability）。特别地，函数 $\varphi$ 的可和性（summability）提供了一个判据，用以区分小初值解的有限时间爆破（blow-up in finite time）与整体存在性（global existence）。最后，当 $\varphi\notin L^1([0,+\infty))$ 且满足某缩放条件（我们称之为一致上缩放条件，uniform upper scaling condition）时，我们推导出局部解的精确生命跨度估计（sharp lifespan estimates）。

7. 无凝胶化与具有正则变化质量交换速率的Boltzmann方程的整体存在性

No Gelation and Global Existence for a Boltzmann Equation with Regularly Varying Mass-Exchange Rates

作者: Siwei Luo, Jian-Guo Liu

类别: math.AP, math-ph

发布日期: 2026-07-27

arXiv ID: 2607.25112v1

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

摘要翻译

中文标题: 无凝胶化与具有正则变化质量交换速率的Boltzmann方程的整体存在性

摘要翻译: 我们发展了一种无凝胶化机制，由此得到具有质量交换的均匀空间Boltzmann方程对Grad截断硬势 $B=E^γb(ξ), 0<γ<1$ 及正则变化质量交换速率的整体解。在不假设任何高阶质量矩的前提下，我们构造了一个由二进制铰链组装而成的凸超线性权函数。大质量且质量相近粒子间碰撞所产生的加权矩，被与一个质量有界粒子均匀填充的储库之间的碰撞所耗散吸收。正则变化性使得这两类碰撞构型中的带符号增量可比，而质量守恒及 $γ<1$ 则给出消失因子 $L^{γ-1}$。由此，我们在任意有限时间区间上获得了质量截断逼近解的一致矩估计，并排除了有限时间内质量逃逸至无穷大的可能性。因此，对每个具有有限物理矩的非负初始密度，我们得到了一个整体积分弱解。

不相关论文（基本信息）

1. Pólya 猜想在欧几里得球上的 Neumann Laplacian 上的成立性

Pólya’s Conjecture for the Neumann Laplacian on Euclidean Balls

作者: Yutian Li

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

发布日期: 2026-07-28

arXiv ID: 2607.25958v1

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

摘要: We prove Pólya’s conjectured lower bound for the Neumann counting function of every Euclidean ball. If $d\ge2$, $R>0$, and $E\ge0$, then \begin{equation*} N_{B_R^d}^{<}(E) \ge \frac{ω_d}{(2π)^d}|B_R^d|E^{d/2} = \frac{(R\sqrt E)^d}{2^dΓ(\frac d2+1)^2}. \end{equation*} The radial boundary condition in dimensions $d\ge3$ is a Dini condition, not a derivative-zero condition. A strict Robin comparison first reduces it to a Bessel phase estimate. Variational bounds handle low frequencies, while estimates based on finitely many radial levels and on beta moments cover the intermediate range uniformly in the dimension. The remaining high-frequency estimate is explicit. The finite rational calculations form part of the proof appendix.

2. 关于具有次线性竞争相互作用的椭圆系统中死核现象的出现

On the emergence of dead cores in elliptic systems with sublinear competitive interactions

作者: Hugo Tavares, Gianmaria Verzini, Junwei Yu

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25943v1

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

摘要: In this paper, we study semilinear elliptic systems with sublinear coupling terms, showing that, under competition-type interactions, solutions typically have dead cores, i.e., they vanish on certain open subsets of the domain. We apply our results to a large class of solutions treated in the literature, for instance to ground states and least energy sign-changing solutions, under Dirichlet boundary conditions or in the whole space. The proofs are based on a general result stating that subsolutions to a certain sublinear equation have dead cores, and on uniform Hölder bounds.

3. Barron-Lipschitz 能量间隙与科学机器学习中的深度分离现象

The Barron-Lipschitz Energy Gap and Depth Separation Phenomena in Scientific Machine Learning

作者: Nima Rezaei, Stephan Wojtowytsch

类别: math.AP, stat.ML

发布日期: 2026-07-28

arXiv ID: 2607.25905v1

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

摘要: We illustrate in several examples that even neural networks of infinite width (specifically, Barron functions) may encounter substantial obstacles when used as a model class for problems in the calculus of variations. An instance of practical relevance concerns the bending, stretching and folding of a thin elastic shell with anchored or clamped boundary conditions where elastic energy could be reduced by folding along a circular line, but the neural networks can only describe straight folds along entire lines. Conversely, we show that there is no gap between the energy that Barron functions and Lipschitz functions can achieve for a large class of integral first-order functionals.

4. 多点热带松弛与Monge–Ampère方程

Many-point tropical relaxation and the Monge–Ampère equation

作者: Nikita Kalinin, Ernesto Lupercio, Higinio Serrano, Mikhail Shkolnikov

类别: math.AP, math.AG, math.MG

发布日期: 2026-07-28

arXiv ID: 2607.25878v1

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

摘要: We construct an incidence-driven tropical approximation of the planar Aleksandrov Monge–Ampère equation. Let $Ω\subset\mathbb R^2$ be a bounded open convex domain, let $K\SubsetΩ$, and let $F_N=G_{P_N}0_Ω$ be the minimal nonnegative concave tropical series with integral slopes, zero boundary values, and corner locus containing an $N$-point set $P_N\subset K$. For universally generic configurations whose empirical measures converge to $μ$, $N^{-1/2}F_N\longrightarrow F_{μ,Ω}$ uniformly on $\overlineΩ$, where $F_{μ,Ω}$ is the unique continuous concave Aleksandrov solution of $\mathrm{MA}(F)=μ$ with zero boundary values. On every compact $L\SubsetΩ$ we prove an $O(N^{-1/2})$ bounded-Lipschitz-type estimate for the curvature discrepancy. No regularity or strict-convexity assumption is imposed on $\partialΩ$. For rational polygons, strong genericity holds on an open dense full-measure locus. The tropical curve has exactly $N$ bounded cells, the marked dual edges form a spanning tree, and every compact internal edge has weight one. These finite statements yield global weak curvature convergence and the exact identity $\mathrm{MA}(F_N)(Ω^\circ)=N-1+\frac{1}{2}D_{\mathrm{term}}(F_N)$, with $D_{\mathrm{term}}(F_N)=O(\sqrt N)$. The proof combines tropical interpolation, semilinear marked topology, an Euler–Pick curvature formula, minimal coefficient deformations, tangential coarea, and a weighted Crofton estimate uniform over rational polygonal exhaustions. We also obtain almost-sure limits for random point clouds, full affine covariance of the continuum solution, and a configuration-dependent Abelian-sandpile diagonal.

5. 三维临界Hartree型Brézis–Nirenberg问题的爆破渐近性

Blow-up asymptotics for a critical Hartree-type Brézis–Nirenberg problem in dimension three

作者: Yue Wu, Minbo Yang

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25759v1

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

摘要: In this paper, we study the following critical Hartree problem \begin{equation}\label{equationabstract} \begin{cases} \displaystyle-Δu+(Q+\varepsilon V)u= A_{3,μ} \left(\int_Ω\frac{u^{6-μ}(y)}{|x-y|^μ}dy\right) u^{5-μ} &\mathrm{in}Ω,\ \displaystyle u>0&\mathrm{in}Ω,\u=0&\mathrm{on}\partialΩ,\end{cases} \end{equation} where $Ω\subset\mathbb{R}^3$ is a bounded open set, $0<μ<2$, the exponent $6-μ$ is the upper critical exponent in the sense of the Hardy–Littlewood–Sobolev inequality and $A_{3,μ}>0$ is a normalization constant. The function $Q$ is assumed to be critical in the sense of Hebey and Vaugon, and the solutions $u_{\varepsilon}$ of \eqref{equationabstract} are assumed to be an optimizing sequence for the Hardy–Littlewood–Sobolev inequality. Under a natural nondegeneracy assumption, we derive a precise asymptotic expansion of $u_{\varepsilon}$, determine the exact blow-up rate, and identify the concentration point. We also obtain the pointwise blow-up behavior both near and away from the concentration point.

6. 任意周期介质中线性弹动力学与弹静力学的良好适定均匀化应变梯度模型

Well-posed homogenized strain-gradient models for linear elastodynamics and elastostatics in arbitrary periodic media

作者: Rémi Cornaggia, Marc Bonnet, Giuseppe Rosi, Saad El Ouafa, Nicolas Auffray

类别: math.AP, math-ph, physics.class-ph

发布日期: 2026-07-28

arXiv ID: 2607.25602v1

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

摘要: This work develops well-posed homogenized strain-gradient models for linear elastostatics and elastodynamics in periodic media, with a primary focus on elastic wave propagation. Using the classical two-scale asymptotic expansion method, we carry out second-order periodic homogenization for media in $\mathbb{R}^d$ ($d = 2, 3$), with no restriction on the periodicity cell geometry or material distribution. Reciprocity identities applied to suitably chosen pairs of cell solutions provide alternative expressions for the effective stiffness and inertia tensors arising at the leading, first and second orders, substantially reducing the number of cell problems that must actually be solved. Since direct two-scale homogenization beyond leading order generically yields ill-posed effective operators, a Boussinesq-trick procedure is introduced, involving a tunable scalar weight, to recast the resulting fourth-order partial differential equation as a valid strain-gradient elasticity (SGE) model possessing the requisite symmetry, sign-definiteness and coercivity properties. These properties are then used, via the Hille-Yosida theorem, to establish the well-posedness of the corresponding transient initial-value and forced-response problems. Several practically relevant special cases are examined, including centrosymmetric cells, homogeneous mass density and homogeneous elasticity, each yielding simplified model structures. Numerical illustrations on three two-dimensional periodicity cells (square, hexagonal and a non-centrosymmetric chiral lattice) compare the resulting dispersion relations against reference Floquet-Bloch computations and assess transient wave propagation, demonstrating the model’s capacity to capture anisotropic and dispersive effects beyond classical elasticity while preserving mathematical well-posedness.

7. 具有性状无关相互作用的性状结构化广义Lotka-Volterra系统的全局动力学

Global Dynamics of Trait-Structured Generalised Lotka-Volterra Systems with Trait-Independent Interactions

作者: Manh Hong Duong, Nataliya Balabanova, Blaine van Rensburg

类别: math.AP, q-bio.PE

发布日期: 2026-07-28

arXiv ID: 2607.25573v1

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

摘要: We study the long time dynamics of a selection-mutation integro-differential Lotka-Volterra system of $N$ populations. In our model, fitness depends on a continuous phenotypic trait, but the effect of one population on another is independent of this trait. We establish that, under some usual assumptions on the interactions between populations, the long time behaviour of solutions is exactly determined by the well studied Generalised Lotka-Volterra (ordinary differential) Equation. We also show that, all else being equal, the minimal mutation rate is selected for in a static environment, whereas in an a changing environment, an intermediate mutation rate could be selected instead. The key step in our proofs is establishing that the total population sizes are asymptotically governed by a Generalised Lotka-Volterra Equation, which then allows the use of a general result on asymptotically autonomous dynamical systems.

8. Lavrentiev 现象下的双相流

Double phase flow under Lavrentiev phenomenon

作者: Yoshiki Kaiho

类别: math.AP

发布日期: 2026-07-28

arXiv ID: 2607.25296v1

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

摘要: This paper deals with parabolic equations associated with double phase functionals. It is known that double phase functionals may exhibit the Lavrentiev phenomenon, which indicates an existence of a singular minimizer. Our aim is to investigate the process of the associated double phase flow evolving toward the singular minimizer. For this purpose, we study whether solutions can be approximated by smooth functions. We first observe a phenomenon that we call finite-time loss of smooth approximability; more precisely, we prove that the flow eventually ceases to be smoothly approximable. We also establish quantitative estimates on the time of loss. On the other hand, we investigate a phenomenon that we call short-time persistence of smooth approximability; we prove that the smooth approximability persists for a short time provided that the initial datum is regular and that the functional is nondegenerate with respect to the gradient variable. The results concerning finite-time loss of smooth approximability are derived from the evolution variational inequality, whereas the short-time persistence result is obtained by applying analytic semigroup theory. The novelty of this paper lies in studying the dynamical aspect of the Lavrentiev phenomenon, which is usually regarded as a stationary phenomenon.

9. 来自穿孔流形 Steklov 谱的不定 Coulomb 相互作用

An indefinite Coulomb interaction from the Steklov spectrum of perforated manifolds

作者: Zhonggan Huang, Raghavendra Venkatraman

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

发布日期: 2026-07-28

arXiv ID: 2607.25211v1

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

摘要: We establish optimal convergence rates for Steklov eigenvalues and harmonically extended eigenfunctions toward their weighted Laplace–Beltrami counterparts on a closed manifold perforated by many small geodesic balls. The holes have radii that scale critically with respect to their spacing in the sense that the boundary area of each hole balances with the weighted volume of its Voronoi cell. We then derive a higher-order expansion of the Steklov eigenvalues; in dimensions two and three, we identify two correction scales. The expansion is governed by an indefinite Coulomb-type energy of the discrepancy between the boundary and bulk measures, mediated by the reduced Green function of the limiting operator. The proof relies on sharp estimates for certain auxiliary functions that we introduce in order to quantify the discrepancy measure between the surface measure on the holes and their background density. Our paper serves to bridge an emerging literature in spectral geometry with one on systems of points interacting via Coulomb-type energies.

10. 退化植物-消费者反应-扩散系统中的传播窗口估计

Propagation-Window Bounds in Degenerate Plant-Consumer Reaction-Diffusion Systems

作者: Chiun-Chuan Chen, Ting-Yang Hsiao, Li-Chang Hung, Haoyuan Li, Shun-Chieh Wang

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.25149v1

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

摘要: We study traveling waves in a class of high-dimensional plant-consumer reaction-diffusion systems with $N$ competing plant species and $N$ associated consumer populations. The consumer equations are linearly degenerate at the extinction state. Under the weak-interaction condition $\|K\|_{\infty}<1$, the system admits a unique positive coexistence equilibrium. Define $ s*{\mathrm P}:=\max_{1\le i\le N}2\sqrt{d_i r_i},\, s_{\mathrm C}^{\mathrm{ex}}:=\min_{1\le i\le N}2\sqrt{D_i e_iσ_i},\, s_{\mathrm C}^{\mathrm{nec}}:=\min_{1\le i\le N}2\sqrt{D_i e_i}$, where $σ_i:=1-\sum_{j=1}^Nκ_{ij}$, and let $\mathcal S$ be the set of speeds admitting a strictly positive traveling wave connecting extinction to coexistence. If $s_{\mathrm P}<s_{\mathrm C}^{\mathrm{ex}}$, we prove $ (s_{\mathrm P},s_{\mathrm C}^{\mathrm{ex}})\subseteq \mathcal S \subseteq [s_{\mathrm P},s_{\mathrm C}^{\mathrm{nec}})$. The existence result is based on a new hybrid lower-solution construction that removes the diffusion-separation and auxiliary compatibility assumptions from the earlier four-component theory. Plant lower profiles combine a polynomial-exponential leading-edge bump with a logistic bridge, while consumer lower profiles are Gaussian near the leading edge and constant behind the front. Nonexistence below $s_{\mathrm P}$ follows from a Sturm-type oscillation argument, whereas leading-edge asymptotics and a Riccati crossing argument yield the universal upper obstruction $s<s_{\mathrm C}^{\mathrm{nec}}$. Finally, for a fixed two-species system, we construct an explicit wave with $ s_{\mathrm P}<s_{\mathrm C}^{\mathrm{ex}}<s_{\mathrm{exact}}<s_{\mathrm C}^{\mathrm{nec}}$, showing that $s_{\mathrm C}^{\mathrm{ex}}$ is only a constructive threshold, not an intrinsic maximal wave speed.

11. Barron 空间中的椭圆正则性理论及其在深度 Ritz 方法中的应用

Elliptic Regularity Theory in Barron Spaces and Applications to the Deep Ritz Method

作者: Stephan Wojtowytsch

类别: math.AP, math.NA, stat.ML

发布日期: 2026-07-27

arXiv ID: 2607.25100v1

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

摘要: We prove that harmonic functions with Dirichlet boundary data in Barron space, a function class tailored to wide ReLU networks with a single hidden layer and suitably bounded weights, are generally neither Lipschitz continuous nor in the Sobolev class $H^2$. A fortiori, they are not in any function class in which the norm controls the Lipschitz constant, which rules out not only Barron space regularity, but also regularity in function classes for deeper ReLU networks with bounded coefficients. They can, however, be approximated to accuracy $\sim \varepsilon$ by Barron functions of low norm $\sim |\log\varepsilon|$ in various Lebesgue and Sobolev norms (with at most two derivatives). The positive result holds on very simple domains: Half-spaces in arbitrary dimension and rectangular domains in two dimensions. As an application of this regularity theory, we obtain a priori error estimates for Deep Ritz neural PDE solvers.

12. 一种面向多物种趋化性中亚临界与超临界动力学的混合物理信息神经网络框架

A Hybrid Physics-Informed Neural Network Framework for Subcritical and Supercritical Dynamics in Multi-Species Chemotaxis

作者: Hamid El Bahja, Jan C. Riedel, Peter Jung

类别: math.NA, math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24949v1

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

摘要: We study a two-species chemotaxis system that exhibits two qualitatively different regimes: subcritical dynamics, in which solutions remain smooth, and supercritical dynamics, in which strong aggregation may lead to finite-time blow-up. In the subcritical regime, we use a standard continuous-time physics-informed neural network (PINN) with alternating training for the coupled fields and show that it provides accurate and efficient approximations. In the supercritical regime, however, this formulation is not sufficiently robust to resolve the highly localized structures and steep gradients associated with blow-up. To address this limitation, we introduce a discrete-time PINN based on backward Euler time stepping, combined with a logarithmic transformation to stabilize large solution values and a residual-adaptive collocation strategy that concentrates training points near regions of high activity. This hybrid framework allows us to use a simple and computationally inexpensive PINN when the solution is smooth, while switching to a more robust discrete formulation when singular dynamics emerge. Numerical experiments confirm that the continuous-time formulation accurately resolves subcritical dynamics, while the discrete-time formulation reliably captures singular aggregation and finite-time blow-up in the supercritical regime.