Home » 论文总结

2026年08月18日论文分析

共有 12 篇重点关注论文, 33 篇论文可以了解, 6 篇不相关论文

August 18, 2026 · 45 min · 9367 words · P0lyno3ial

Table of Contents

重点关注论文（完整分析）
了解领域论文（摘要翻译）
不相关论文（基本信息）

生成时间: 2026年08月18日 03:31:07

运行状态: 完成

计划处理论文数量: 51

已完成论文数量: 51

未完成论文数量: 0

重点关注论文数量: 12

了解领域论文数量: 33

不相关论文数量: 6

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

1. 二维定常低耗散 Navier–Stokes 系统的 Liouville 定理

A Liouville theorem for the two-dimensional stationary hypodissipative Navier–Stokes system

作者: Nicola De Nitti, Lukas Niebel, Jiaqi Yang

类别: math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16437v1

摘要翻译: 我们研究定义在 $\mathbb R^2$ 上的具有分数耗散 $$(-Δ)^s$$ 的二维定常不可压缩 Navier–Stokes 方程。对于整个范围 $s \in (0,1)$ ，我们证明了每一个满足自然能量条件 $u\in\dot{\mathrm H}^s(\mathbb R^2;\mathbb R^2)$ 的光滑解都有 $u\equiv0$ 且压力为常数。这是 Gilbarg 和 Weinberger 在 $$s=1$$ 时的平面有限 Dirichlet 定理的分式对应版本。证明在三个范围内采用了不同的论证方法。对于 $0<s<\frac13$ ，我们将从方程导出的 $\mathrm L^2$ -估计与流函数截断论证相结合。对于 $\frac13\leq s\leq\frac23$ ，我们使用局部能量估计，其边界项支撑在扩张的圆环上。对于 $\frac23<s<1$ ，我们通过 Lorentz 空间 bootstrap 建立正则性和衰减性，然后对涡度应用极大值原理。此外，我们在包含 $u\in \mathrm L^r(\mathbb R^2)$ 对所有 $1\le r\le2$ 的圆环增长条件下，结合 Bernoulli 恒等式与截断论证，处理了 $$s=0$$ 时的定常阻尼 Euler 系统。

Comment: 无

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

1. 研究对象和背景

本文研究的是定义在全平面 $\mathbb R^2$ 上的二维不可压缩定常分数阶纳维 - 斯托克斯（Navier–Stokes）方程组，其核心特征在于扩散项被分数阶拉普拉斯算子 $(-\Delta)^s$ 取代，其中参数 $s \in (0,1)$ 。该方程组的形式为：

\begin{cases} u \cdot \nabla u + (-\Delta)^s u + \nabla p = 0, & x \in \mathbb R^2, \\ \text{div } u = 0, & x \in \mathbb R^2, \end{cases}

其中 $u: \mathbb R^2 \to \mathbb R^2$ 是速度场， $p: \mathbb R^2 \to \mathbb R$ 是压力场。当 $$s=1$$ 时，该系统退化为经典的二维定常纳维 - 斯托克斯方程；当 $s \to 0$ 时，则趋近于带线性阻尼的欧拉方程（damped Euler system）。分数阶算子的定义为

(-\Delta)^s u(x) = c_{2,s} \, \text{p.v.} \int_{\mathbb R^2} \frac{u(x) - u(y)}{|x-y|^{2+2s}} \, dy,

其傅里叶表示为 $\widehat{(-\Delta)^s u}(\xi) = |\xi|^{2s} \widehat{u}(\xi)$ 。

该问题的物理与数学背景在于理解“弱耗散”流体在无限域上的刚性行为。对于经典情形 ( $$s=1$$ )，Gilbarg 和 Weinberger 证明了有限狄利克雷积分（即 $u \in \dot H^1(\mathbb R^2)$ ）的解必为零（刘维尔定理），这是二维流体力学中的经典结果。然而，当耗散减弱至 $$s < 1$$ 时，能量空间 $\dot H^s(\mathbb R^2)$ 的嵌入性质发生显著变化。尽管由索博列夫嵌入可知 $\dot H^s(\mathbb R^2) \hookrightarrow L^{q_s}(\mathbb R^2)$ ，其中 $q_s = \frac{2}{1-s}$ ，但这一嵌入并不足以直接提供 $$L^2$$ 控制或足够强的点态衰减，使得 $$u$$ 本身不能直接作为全局测试函数来验证形式上的能量恒等式 $\|(-\Delta)^{s/2}u\|_{L^2}^2 = 0$ 。

本文旨在填补的理论空白是：在仅假设自然能量条件 $u \in \dot H^s(\mathbb R^2; \mathbb R^2)$ 的前提下，证明对于全范围 $s \in (0,1)$ ，任何光滑解必须满足 $u \equiv 0$ 且压力为常数。这不仅是经典 Gilbarg-Weinberger 定理在分数阶耗散情形下的推广，也是对 $s \in (0,1)$ 区间内不同正则性层级下解的唯一性机制的系统性刻画。此前关于分数阶情形的研究多集中于三维情形，或者在二维情形中附加了额外的勒贝格可积性条件（如 $u \in L^2 \cap \dot H^s$ ），而本文试图证明仅凭 $\dot H^s$ 能量条件即可导出平凡解，从而确立了该模型在二维空间中的最优刘维尔定理。

2. 主要定理或主要结果

本文的核心成果是建立了覆盖全参数范围 $s \in (0,1)$ 的刘维尔定理，并补充处理了端点 $$s=0$$ 的情形。主要结论如下：

定理 1.1 (分数阶纳维 - 斯托克斯系统)：设 $s \in (0,1)$ ，若 $(u, p) \in C^\infty(\mathbb R^2; \mathbb R^2) \times C^\infty(\mathbb R^2)$ 是方程 (1.1) 的光滑解，且满足自然能量条件

u \in \dot H^s(\mathbb R^2; \mathbb R^2),

则必有 $u \equiv 0$ 且 $$p$$ 为常数。

定理 1.2 (零阶端点：带阻尼欧拉系统)：对于 $$s=0$$ 的情形，方程退化为定常带阻尼欧拉系统：

\begin{cases} u \cdot \nabla u + u + \nabla p = 0, \\ \text{div } u = 0. \end{cases}

若光滑解 $$(u, p)$$ 满足环形增长条件

\liminf_{R \to \infty} \frac{1}{R} \int_{B_{2R} \setminus B_R} |u(x)| \, dx = 0,

则同样有 $u \equiv 0$ 且 $$p$$ 为常数。特别地，若 $u \in L^r(\mathbb R^2; \mathbb R^2)$ 对某个 $1 \le r \le 2$ 成立，则该条件自动满足，结论亦成立。

这些定理表明，在二维空间中，无论耗散强度如何（只要 $$s>0$$ ），只要速度场具有有限的分数阶能量，就不存在非平凡的定常解。证明过程根据 $$s$$ 的取值范围分为三个互补的策略：

低耗散区 ( $$0 < s < 1/3$$ )：利用方程导出 $u \in L^2$ ，构造流函数截断，证明能量恒等式成立。
中耗散区 ( $1/3 \le s \le 2/3$ )：使用局部化能量估计，边界项支撑在扩张的圆环上，利用 $R^{3s-2}$ 的尺度关系闭合估计。
高耗散区 ( $$2/3 < s < 1$$ )：通过洛伦兹空间的自举法建立正则性与衰减，应用涡度的极大值原理导出 $\omega = 0$ 。

所有结论均为严格证明的数学定理，而非启发式推导。

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

证明的主要技术障碍在于：在仅假设 $u \in \dot H^s(\mathbb R^2)$ 的情况下，无法直接验证形式能量恒等式 $\langle (-\Delta)^s u, u \rangle = 0$ ，因为 $$u$$ 可能不具备足够的局部可积性或衰减性以支持全局积分操作。此外，非线性对流项 $u \cdot \nabla u$ 和压力项的处理在不同 $$s$$ 值下表现出截然不同的分析特性。作者的创新之处在于针对 $$s$$ 的三个不同区间开发了三种独立但互补的分析框架，并利用归一化压力和流函数截断技术克服了传统方法的局限性。

关键技术点 1：归一化压力的识别与处理 首先，作者通过 Calderón-Zygmund 理论将压力 $$p$$ 替换为归一化压力 $$p_0 = R_i R_j (u_i u_j)$$ ，其中 $$R_i$$ 为黎斯变换。引理 2.1 证明了 $$p - p_0$$ 为常数，从而可以将原方程改写为

(-\Delta)^s u + u \cdot \nabla u + \nabla p_0 = 0.

这一技巧消除了压力项在分布意义下的多项式模糊性，使得后续的能量估计和投影运算可以在 $S'(\mathbb R^2)$ 框架下严格进行。特别是，它允许作者使用 Leray 投影 $$P$$ 将方程转化为

(-\Delta)^s u + P \text{div}(u \otimes u) = 0,

从而避免了直接处理梯度项带来的困难。

关键技术点 2： $$0 < s < 1/3$$ 区域的流函数截断策略 在此区间，由于 $$s$$ 较小，直接的能量估计难以闭合。作者首先利用投影方程和傅里叶分析证明 $u \in L^2(\mathbb R^2)$ （引理 3.1）。随后，引入流函数 $\psi$ 使得 $u = \nabla^\perp \psi$ ，并证明 $\psi \in H^{1+s}(\mathbb R^2)$ 且具有连续衰减性（引理 3.2）。核心创新在于构造截断函数 $C_\varepsilon(\psi)$ （引理 2.6），生成紧支集、无散度的测试向量场 $z_\varepsilon = \nabla^\perp (\psi - C_\varepsilon(\psi))$ 。由于 $z_\varepsilon$ 的支集性质，对流项和压力项在测试过程中精确消失（或趋于零）。令 $\varepsilon \to 0$ ，利用 $$H^s$$ 空间的弱下半连续性恢复完整的能量恒等式 $\|(-\Delta)^{s/2}u\|_{L^2}^2 = 0$ 。该方法巧妙地避开了对 $$u$$ 的直接 $$L^2$$ 加权估计，转而利用流函数的几何结构。

关键技术点 3： $1/3 \le s \le 2/3$ 区域的局部化能量估计 在此中间区间，作者采用截断函数 $\chi_R$ 对能量恒等式进行局部化。测试函数取为 $u \chi_R$ ，其中 $\chi_R$ 在 $$B_R$$ 外截断。关键难点在于处理由此产生的边界项（annular boundary terms）。通过精细的尺度分析，发现边界通量项的量级约为 $R^{3s-2}$ 。当 $s \ge 1/3$ 时，非线性项的指数 $qs \ge 3$ 保证了边界项的可控性；而当 $s \le 2/3$ 时，因子 $R^{3s-2} \to 0$ 保证了随着半径扩张，边界贡献消失。这一估计依赖于 $u \in \dot H^s$ 导出的 $u \in L^{q_s}$ 以及洛伦兹空间的 Hölder 不等式，成功在无需额外 $$L^r$$ 假设的情况下闭合了能量估计。

关键技术点 4： $$2/3 < s < 1$$ 区域的洛伦兹空间自举与极大值原理 对于高耗散区域，直接的能量方法失效，因为 $R^{3s-2}$ 不再衰减。作者转而利用算子 $T = (-\Delta)^{-s} P \text{div}$ 的映射性质（引理 2.5）。通过洛伦兹空间 $L^{p,r}$ 的迭代论证（bootstrap argument），逐步提升解的正则性：从 $$L^p$$ 到 $L^\infty$ ，再到 Hölder 连续空间 $C^{0,\theta}_b$ ，最终获得 $u \in C^{1,\gamma}_b$ 及其梯度的衰减。一旦建立了 $$u$$ 和 $\nabla u$ 在无穷远处的衰减，即可定义涡度 $\omega = \partial_1 u_2 - \partial_2 u_1$ 并证明其在无穷远处为零。此时，利用分数阶微积分中的极大值原理（maximum principle for nonlocal operators），结合 $\omega$ 的连续性，直接推出 $\omega \equiv 0$ ，进而由旋度为零及无散条件导出 $u \equiv 0$ 。这一路径完全避开了能量恒等式的直接验证，转而依赖解的精细正则性和极值性质。

技术总结：整个证明框架并非单一方法的简单延伸，而是根据 $$s$$ 值导致的 Sobolev 嵌入临界指数变化（特别是 $$q_s$$ 与 3 的关系，以及 $$3s-2$$ 的符号），动态切换了流函数截断、局部化能量估计和正则性自举三种截然不同的分析工具。

4. 与之前工作的比较

本文与既有文献相比，在假设条件的宽松度和适用范围上取得了显著推进，主要体现在以下几个方面：

1. 假设条件的最小化 既往关于二维分数阶纳维 - 斯托克斯方程的研究（如 Wang 和 Xiao [29]）通常要求解同时满足 $u \in \dot H^s(\mathbb R^2)$ 和 $u \in L^2(\mathbb R^2)$ ，甚至在某些情况下需要 $u \in L^6(\mathbb R^2)$ 。例如，[29] 的结果在 $$s < 1/2$$ 时需要 $u \in \dot H^s \cap L^2 \cap L^6$ 。相比之下，本文定理 1.1 仅需假设 $u \in \dot H^s(\mathbb R^2)$ 即可证明 $u \equiv 0$ 。这是目前已知最弱的能量假设，直接对应于系统自然的能量空间，去除了人为添加的 $$L^2$$ 或高阶可积性限制。这一改进填补了从“有限能量加额外可积性”到“纯有限能量”之间的理论空白。

2. 与经典 Gilbarg-Weinberger 定理的对齐 经典结果指出，对于 $$s=1$$ 的二维定常 NS 方程，有限狄利克雷积分（ $u \in \dot H^1$ ）蕴含平凡解。本文的工作将这一结论完美推广至 $s \in (0,1)$ 的全范围，建立了分数阶版本的 Gilbarg-Weinberger 定理。此前，虽然有一些针对三维分数阶情形的结果（如 Lee and Lee [19]），但其适用范围往往受限于特定的 $$s$$ 区间（如 $s \ge n/6$ ）或需要额外的衰减假设（如 $u(x) \to 0$ 一致收敛）。本文在二维情形下不仅放宽了衰减要求（仅需能量条件隐含的衰减），而且覆盖了 $$s$$ 的整个开区间 $$(0,1)$$ ，实现了理论的完备性。

3. 对端点 $$s=0$$ 的拓展 对于 $$s=0$$ 的带阻尼欧拉系统，现有文献（如 Chepyzhov et al. [10] 和 Chae [7]）通常需要更强的正则性假设（如 $u \in H^1$ ）或更高的勒贝格可积性（如 $u \in L^q, q>6$ ）。本文定理 1.2 在 $$s=0$$ 处引入了一个基于环形平均的温和增长条件，并证明了若 $u \in L^r$ ( $1 \le r \le 2$ ) 则结论成立。这一结果将阻尼欧拉系统的唯一性理论推进到了更低的可积性层级，且无需 $$H^1$$ 正则性，展示了作者在处理零阶耗散时的独特技巧（Bernoulli 恒等式结合截断）。

4. 方法与代价的比较 本文的推进是以牺牲证明技术的统一性为代价的。不同于以往文献可能尝试寻找一种通用的能量方法来覆盖所有 $$s$$ ，本文不得不根据 $$s$$ 的取值分段处理：

在 $$0 < s < 1/3$$ 时，需引入流函数和复杂的截断构造；
在 $1/3 \le s \le 2/3$ 时，需依赖精细的圆环边界项估计；
在 $$2/3 < s < 1$$ 时，需转向洛伦兹空间自举和极大值原理。这种分段策略反映了不同 $$s$$ 值下分析性质的本质差异，但也意味着本文的方法论不如某些通用方法简洁。然而，这种代价换来了在最优假设条件下的完整证明，是理论深度的体现。

5. 与三维结果的对比 在三维情形中，即使对于 $$s=1$$ ，仅凭 $u \in \dot H^1$ 仍不足以证明刘维尔定理（这是一个著名的开放问题），通常需要额外的空间 $$X$$ （如 $L^{9/2}$ 或 BMO 类空间）。本文明确指出二维情形的特殊性：在二维中， $\dot H^s$ 的自然嵌入性质（特别是 $$q_s$$ 随 $$s$$ 变化的行为）使得纯能量假设足以排除非平凡解，而三维则不然。这一对比突显了二维问题的独特优势，也解释了为何本文能取得比三维情形更完美的结果。

综上所述，本文通过精细的分段分析和创新的工具组合，将二维弱耗散纳维 - 斯托克斯系统的刘维尔定理推进到了理论极限，即在仅保留自然能量假设的前提下，证明了全范围 $s \in (0,1]$ 及端点 $$s=0$$ 的解的唯一性。

2. 三阶对称尺度不变欧拉流的梯度增长与跳跃剖面的弛豫

Gradient growth and relaxation to jump profiles for 3-fold symmetric scale-invariant Euler flows

作者: Daomin Cao, Junhong Fan, Guolin Qin

类别: math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16755v1

摘要翻译: 我们考虑二维欧拉方程在三阶对称情形 $$m=3$$ 下的零次齐次约化，该问题在 Said、Elgindi 和 Murray extcite{EMS} 关于 $$m extgreater=4$$ 的弛豫理论中尚未解决。我们证明：对于 $$1<p extless= extinfty$$ ，每一个非常数的 $W^{1,p}$ 解都满足当 $$t o extpm extinfty$$ 时 $

Comment: 无

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

1. 研究对象和背景

二维不可压缩欧拉方程的涡度形式为 $\partial_t \omega + u \cdot \nabla \omega = 0$ ，其中速度场由流函数通过 $u = \nabla^\perp \Delta^{-1} \omega$ 恢复。尽管 Yudovich 理论保证了有界涡度的全局适定性，但长时动力学中是否产生无限时间尺度的小尺度结构（即梯度增长）以及解是否趋向于某种稳态或奇异结构，仍是分析领域的核心难题。在自由边界设置下，已知存在线性、超线性甚至双指数级的梯度增长机制，但在无边界的全平面上，对于具有特定对称性的解，其渐近行为更为微妙。

本文研究的是零齐次涡度解 $\omega(t, r, \theta) = g(t, \theta)$ ，这种解在空间无穷远处不衰减，因此需要离散旋转对称性来保证全平面公式的良定义性。考虑 $$m$$ -重对称情形，即 $g(\theta + 2\pi/m) = g(\theta)$ 。在此假设下，流函数可设为 $\psi(t, r, \theta) = r^2 G(t, \theta)$ ，代入欧拉方程得到关于角向变量 $\theta$ 的约化方程：

(\partial_{\theta\theta} + 4)G = g, \quad g_t + 2G g_\theta = 0

其中算子 $A = \partial_{\theta\theta} + 4$ 在周期区间 $$T_L$$ ( $L=2\pi/3$ ) 上定义了 $$G$$ 与 $$g$$ 的关系。该方程组描述了尺度不变涡度的角向演化。

此前 Said, Elgindi 和 Murray [33] 建立了 $m \geq 4$ 时的弛豫理论，证明了在 $m \geq 4$ 时，非平凡解的角导数范数随时间趋于无穷，且解趋向于跳跃剖面（jump profiles）。然而，当 $$m=3$$ 时，问题处于临界状态。关键在于对称化后的格林核 $K_3^{1/(\partial_{\theta\theta}+4)}$ 不再保持定号。具体而言，该核在 $[-\pi/3, \pi/3]$ 上的表达式为：

K_3^{1/(\partial_{\theta\theta}+4)}(\theta) = \frac{3\pi}{10} \left| \sin \frac{3\theta}{2} \right| - \frac{7}{20}

该核在 $$m=3$$ 时变号，导致基于核正定性的传统弛豫论证失效。此外， $$m=2$$ 是共振情形，而 $$m=3$$ 是非共振情形中最低阶的，也是第一个无需额外归一化即可唯一确定约化动力学的情况。本文旨在填补这一空白，解决 $$m=3$$ 情形下的梯度增长与渐近结构问题，回答 Drivas 和 Elgindi 综述中提出的关于小尺度产生及紧轨道生成的 Problem 5。

2. 主要定理或主要结果

本文针对 $$m=3$$ 的零齐次欧拉流，建立了完整的梯度增长理论与跳跃剖面弛豫理论。主要结果如下：

定理 1.1 (双向梯度增长)：设 $p \in (1, \infty]$ ，初始数据 $g_0 \in W^{1,p}(T_L)$ 且非常数。则对应的唯一完全解满足：

\lim_{t \to +\infty} \|g_\theta(t)\|_{L^p} = +\infty, \quad \lim_{t \to -\infty} \|g_\theta(t)\|_{L^p} = +\infty

对于端点情形 $$p=1$$ ，总变差守恒，但若初始数据满足 $L \log L$ 模有限（即 $\int_{T_L} \Phi(|D_\theta g_0|) d\theta < \infty$ ，其中 $\Phi(\xi)=\xi \log(e+\xi)$ ），则该模随时间趋于无穷：

\int_{T_L} \Phi(|g_\theta(t)|) d\theta \to \infty \quad (t \to \pm \infty)

定理 1.3 (测度分量弛豫)：若 $D_\theta g_0$ 可分解为非原子单符号分量与原子的可和级数（见定义 1.1），记正负分量的质量分别为 $$m_j^+, m_j^-$$ 。随着 $t \to +\infty$ ，正分量 $\nu_j^+(t)$ 弱*收敛至单原子测度集 $$A_1(m_j^+)$$ ，负分量 $\nu_j^-(t)$ 收敛至至多两个原子的测度集 $$A_2(m_j^-)$$ ；反之，当 $t \to -\infty$ 时角色互换。更具体地，对于任意 $1 \leq r < \infty$ 和 $0 \leq \alpha < 1/r$ ，解在 Sobolev-Slobodeckij 空间中逼近步函数集：

\text{dist}_{W^{\alpha,r}}(g(t), P_{N^+}(\|g_0\|_\infty)) \to 0 \quad (t \to +\infty)

\text{dist}_{W^{\alpha,r}}(g(t), P_{N^-}(\|g_0\|_\infty)) \to 0 \quad (t \to -\infty)

其中 $$P_N(B)$$ 表示至多有 $$N$$ 个跳跃的有界函数集合。

定理 1.4 (极限集结构)：在上述分解条件下，正向和反向的 $\omega$ -极限集 $\Omega^\pm(g_0)$ 非空、流不变且在 $W^{\alpha,r}$ 中紧致。极限集中的每个剖面 $$h$$ 的导数 $$Dh$$ 是纯原子的（purely atomic），即不含绝对连续部分和 Cantor 部分。具体形式为：

Dh = \sum_{j \in J^+} m_j^+ \delta_{x_j^+} - \sum_{j \in J^-} \rho_j^- + \sum_{k \in J_0} c_k \delta_{y_k^+}

其中 $\rho_j^- \in A_2(m_j^-)$ 表示至多两个原子的测度。

推论 1.5 (光滑数据的完备性)：对于非常数 $$C^1$$ 初值，不仅梯度范数双向发散，且任何 $$L^2$$ 弱极限生成的完整轨道在 $W^{\alpha,r}$ 中是预紧的（precompact）。

定理 1.8 (有限跳跃轨道分类)：对于具有三个非零跳跃的初值，其轨道在模去刚性旋转意义下是同宿轨道（heteroclinic），区间长度随时间严格单调变化，最终趋向于具有更少跳跃的状态（如两个跳跃或常数），不存在真正的周期轨道。

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

证明的核心障碍在于 $$m=3$$ 时对称化核 $$K_3$$ 失去正定性，使得基于能量耗散或核正定性的经典方法失效。作者引入了两种全新的分析工具来分别处理梯度增长和弛豫问题。

1. 熵恒等式与双向梯度增长：为了证明梯度增长，作者没有直接使用 $$L^p$$ 能量估计，而是构造了一个类似于熵的泛函 $H(t) = \int_{T_L} q \log |q| d\theta$ ，其中 $q = g_\theta$ 。通过对约化方程求导并利用算子 $A = \partial_{\theta\theta} + 4$ 的性质，推导出了 $$H(t)$$ 的时间导数恒等式：

H'(t) = 2 \int_{T_L} (G_\theta^2 - 4G^2) d\theta \geq 10 \|G_\theta(t)\|_{L^2}^2

这一不等式表明，只要解不是常数，熵泛函就会随时间增加。结合振荡守恒律（oscillation conservation）和固定时间区间上的紧性，排除了导数范数有界的可能性，从而确立了双向的梯度增长。对于 $$p=1$$ 的端点情形，利用 $L \log L$ 模的发散性替代了 $$L^1$$ 范数的守恒性，这是证明中的关键精细估计。

2. 射影坐标下的四点交叉比单调性：针对弛豫问题，由于核的正定性缺失，作者放弃了传统的谱分析，转而采用几何方法。将角度坐标变换为射影坐标 $z = \tan(\theta - c)$ ，并考察四个有序粒子构成的交叉比：

R = \frac{(x_2 - x_1)(x_4 - x_3)}{(x_3 - x_1)(x_4 - x_2)}

通过计算 $$R$$ 的时间导数，发现其可以表示为第三阶导数 $$D_z^3 U$$ 与一个 Peano 核 $$P_x(s)$$ 的积分：

\frac{d}{dt} \log R = \int_{[x_1, x_4]} P_x(s) d(D_z^3 U)(s)

其中 $$P_x > 0$$ 在区间内成立。这一构造巧妙地利用了 $$D_z^3 U$$ 的符号与 $D_\theta g$ 的符号一致这一事实，从而在不依赖核正定性的情况下，获得了正负分量质量的“收缩”估计。具体来说，这证明了正分量在向前时间中必须坍缩为单原子，负分量坍缩为至多两个原子。这种四点交叉比单调性是本文最核心的创新点，它替代了 $m \geq 4$ 情形下的核正定性假设。

3. 测度分解与紧性论证：作者详细分析了 $$BV$$ 空间中测度的结构，利用上述收缩估计排除了极限测度中的绝对连续部分和 Cantor 部分，证明极限必须是纯原子的。同时，利用时间平移的紧性和极限集的不变性，建立了从弱收敛到强收敛（在 $W^{\alpha,r}$ 中）的桥梁，确保了极限轨道的完整性与预紧性。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上实现了显著推进：

1. 对称性范围的突破： Said, Elgindi 和 Murray [33] 的工作主要解决了 $m \geq 4$ 的情形。他们利用了对称化核的正定性来建立弛豫理论。然而，正如文中指出的，当 $$m=3$$ 时，核 $$K_3$$ 变号，导致 [33] 中的论证框架直接失效。本文首次攻克了 $$m=3$$ 这一临界非共振情形，填补了 $m \geq 3$ 理论体系中的最后一块拼图。这使得关于尺度不变欧拉流的梯度增长与紧轨道结论能够覆盖所有非共振情形（ $m \geq 3$ ）。

2. 技术方法的革新：在 $m \geq 4$ 的研究中，核正定性是核心工具，它允许直接控制测度的分布。而在 $$m=3$$ 的情况下，作者无法沿用此路，必须发明新的几何工具。本文引入的“四点交叉比单调性”（four-point cross-ratio monotonicity）是区别于前人工作的根本创新。该方法不依赖于全局核的符号，而是利用局部粒子的相对位置关系和第三阶导数的符号性质，成功导出了原子化界限。相比之下，Alazard 和 Said [2] 在无对称限制下证明了 $$H^s$$ 范数的无穷大 limsup 增长，但并未给出像本文这样精确的双向增长结果（two-sided limit）以及具体的极限剖面结构（跳跃剖面）。

3. 结论的强度与适用范围：本文的结果比之前的稀疏增长结果更强。例如，Chen 和 Sun [5] 讨论了线性时间尺度的细丝化，Denisov [10, 11] 构造了超线性增长，但这些多为特例或仅涉及 limsup。本文证明了对于所有非常数 $W^{1,p}$ 解，梯度范数在 $t \to \pm \infty$ 时均趋于无穷，且给出了精确的 $L \log L$ 端点行为。在弛豫方面，本文不仅证明了收敛，还精确刻画了极限集的结构（纯原子测度）和收敛速率（在 $W^{\alpha,r}$ 中距离趋于零），并证明了极限轨道的预紧性。这对于理解 Euler 方程在无限时间尺度下的“粗粒化”行为至关重要。

4. 代价与局限：本文的结论依赖于 $$m=3$$ 的特殊对称性以及初值的 $$BV$$ 分解条件（尽管对 $$C^1$$ 数据自动满足）。虽然 $$m=3$$ 是最小的非共振对称，但其结构的复杂性（核变号）导致了证明技术的巨大飞跃。对于 $$m=2$$ 的共振情形，本文的方法并不适用，因为此时算子 $(\partial_{\theta\theta}+4)$ 的核非平凡，需要额外的归一化条件，这留待未来研究。此外，虽然本文处理了 $m \geq 3$ 的所有情况，但对于一般非对称或低正则性数据，其结果的普适性仍需进一步探索。

3. 具有锐角接触角的可压缩毛细 - 重力水波局部适定性

Local Well-Posedness for Compressible Capillary-Gravity Water Waves with Acute Contact Angles

作者: Jingchi Huang, Shanmu Li, Chao Wang

类别: math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16007v1

摘要翻译: 本文旨在研究二维有界角域内具有锐角接触角的可压缩欧拉方程的局部适定性。该构型描述了一个自由表面与固定底部在两点相交的情形，其中流体处于重力场中，且流体与空气之间的界面受到毛细力的影响。当接触角小于 $$π/2$$ 时，我们建立了解的局部存在理论，并在接触点处出现耗散效应。主要的分析挑战源于接触点处的奇异性，这使得处理可压缩自由边界问题的先前方法不再适用。为克服这一困难，我们首先建立可压缩欧拉方程的几何结构，该方法最初由 Shatah 和 Zeng \cite{Shatah2008} 引入用于不可压缩流体。此外，我们提供了角域内波动方程的奇异性分析，从而保证了角域附近计算的有效性。最后，基于几何结构和奇异性分析，我们获得了先验能量估计。利用这些估计，我们在几何表述下证明了该系统的局部适定性。据我们所知，这是首个针对涉及接触点的自由边界可压缩欧拉方程所得到的结果。

Comment: 无

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

1. 研究对象和背景

本文研究的是二维有界角域内可压缩欧拉方程在自由边界与固定底面相交且形成锐角接触点情形下的局部适定性问题。该物理模型描述了一个受重力场作用的流体，其自由表面与固定底面在两个接触点 $$X_l$$ 和 $$X_r$$ 处相交，且接触角 $\theta_i(t) \in (0, \pi/2)$ 。此时，流体界面不仅受重力影响，还受到毛细力（表面张力）的作用。控制方程组为可压缩欧拉方程：

\begin{cases} \rho(\partial_t v + v \cdot \nabla v) + \nabla P = -\rho g e_2 & \text{in } \Omega_t, \\ \partial_t \rho + \operatorname{div}(\rho v) = 0 & \text{in } \Omega_t, \end{cases}

其中 $$v=(v_1, v_2)^T$$ 为速度场， $$P$$ 为压强， $\rho$ 为密度， $$g>0$$ 为重力加速度常数， $$e_2=(0,1)^T$$ 。状态方程取等熵形式 $P = A\rho^\gamma - B$ （ $A,B>0, \gamma>1$ ），为简化记号设 $A=1/2, \gamma=2$ ，则 $\nabla P / \rho = \nabla \rho$ 。

边界条件包括：

自由表面运动学条件：物质导数 $D_t = \partial_t + v \cdot \nabla$ 沿自由表面 $\Gamma_t$ 切向；
动力学边界条件（杨氏方程的推广）： $P|_{\Gamma_t} = \sigma \kappa$ ，其中 $\sigma$ 为表面张力系数， $\kappa$ 为平均曲率；
固定底面无穿透条件： $v \cdot n_b = 0$ 于 $\Gamma_b$ ；
接触点处的动态接触角条件（滑移速度由未平衡的杨应力主导）： $\mp \beta_c v \cdot \tau_b = \sigma(\cos \theta_s - \cos \theta_i) \quad \text{at } X_i, i=l,r,$ 其中 $\beta_c$ 为有效摩擦系数， $\theta_s$ 为平衡接触角， $\theta_i$ 为动态接触角。

该问题的核心难点在于接触点处的奇异性。由于几何域的角点结构，传统的基于拉格朗日坐标或切向-法向分解的能量估计方法失效，因为接触点附近切向方向定义不明确且随时间变化。此外，可压缩性引入了声波与角点奇异性的耦合，使得椭圆正则性理论在角点附近退化，导致解的正则性损失。已有文献多集中于不可压缩流体或光滑边界情形，对于可压缩流体结合接触点动力学的研究尚属空白。本文旨在填补这一理论缺口，建立在该配置下的局部存在性与唯一性理论。

2. 主要定理或主要结果

本文证明了在锐角接触角条件下，可压缩欧拉方程自由边界问题的局部适定性。主要定理如下：

定理 1.1：假设初始数据 $(\rho_0, v_0, \Gamma_0) \in H^3(\Omega_0) \times H^{2+}(\Omega_0) \times H^{4.5}$ ，满足能量泛函 $E(0) < \infty$ ，初始接触角 $\theta_i(0) \in (0, \pi/2)$ ( $$i=l,r$$ )，且初始密度 $\rho_0$ 满足正性条件 $\rho_0 \ge c_0 > 0$ 及非退化条件：

-(\nabla_{\tau_b}\rho + g e_2 \cdot \tau_b)|_{(0,X_l)} \ge c_2, \quad (\nabla_{\tau_b}\rho + g e_2 \cdot \tau_b)|_{(0,X_r)} \ge c_3,

同时满足相容性条件（对 $$j=0,1,2,3$$ ）：

\begin{cases} D_t^j P(\rho)|_{t=0} = \sigma D_t^j \kappa|_{t=0} & \text{on } \Gamma_0, \\ D_t^j (v \cdot n_b)|_{t=0} = 0 & \text{on } \Gamma_b, \\ \mp \beta_c D_t^j v \cdot \tau_b|_{t=0} = \sigma D_t^j (\cos \theta_s - \cos \theta_i)|_{t=0} & \text{at } X_i(0). \end{cases}

则存在时间 $$T>0$$ ，使得系统存在唯一解 $(v, \rho, \Omega_t)$ 满足以下先验估计：

\sup_{t \in [0,T]} E(t) + \int_0^T F(t) \, dt \le P(E(0)) + \int_0^T P(E(t)) \, dt,

其中 $$P$$ 是依赖于 $\sigma, \beta_c, \|\Gamma_b\|_{H^{4.5}}$ 和 $$g$$ 的正系数多项式， $$T$$ 仅依赖于上述量及初始数据的范数。

能量泛函定义为：

E(t) = \sum_{j=0}^3 \|D_t^j \rho\|_{L^2(\Omega_t)}^2 + \sum_{j=0}^2 \|\nabla D_t^j \rho\|_{L^2(\Omega_t)}^2 + \sigma \sum_{j=0}^1 \|\nabla_{\tau_t} D_t^j \nabla_{n_t} \rho\|_{L^2(\Gamma_t)}^2 + \|\omega\|_{H^2(\Omega_t)}^2 + \|v\|_{H^{3/2}(\Omega_t)}^2 + \|\Gamma_t\|_{H^1}^2.

耗散项 $$F(t)$$ 包含接触点处的三项：

F(t) = \sum_{i=l,r} |\sin \theta_i \nabla_{\tau_t} \nabla_{n_t} \rho|_{X_i}|^2 + \sum_{i=l,r} |(\sin \theta_i)(\nabla_{\tau_t} D_t \nabla_{n_t} \rho - 2(\nabla_{\tau_t} v \cdot n_t) \Delta_{\Gamma_t} \rho)|_{X_i}|^2 + \left| \frac{1}{\sin \theta_i} (-\nabla_{\tau_b} \rho - g e_2 \cdot \tau_b) \rho |D_t^2 \rho|^2 \right|_{X_l}^{X_r}.

该结果首次确立了含接触点的可压缩欧拉方程的自由边界问题的局部适定性，并揭示了接触点处特有的耗散机制。

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

证明的核心障碍在于接触点引起的奇异性破坏了传统能量估计中的切向-法向分解的有效性，且可压缩性导致声波与角点奇异性耦合，使得椭圆正则性理论在角点附近失效。作者的创新点主要体现在以下三个方面：

可压缩欧拉方程的几何结构构建：借鉴 Shatah 和 Zeng [35] 针对不可压缩流体的几何框架，作者将其推广至可压缩情形。通过引入物质导数与曲率的演化关系，建立了密度 $\rho$ 的波动方程与曲率 $\kappa$ 的演化方程。具体地，利用运动学边界条件 $D_t \kappa = -\Delta_{\Gamma_t} v \cdot n_t - 2 \nabla_{\tau_t} n_t \cdot \nabla_{\tau_t} v$ ，结合动量方程推导出密度的高阶演化方程：
$\frac{\rho}{\sigma} D_t^2 \rho + A_{com}(\rho) = 2(\nabla_{\tau_t} v \cdot n_t) \nabla_{\tau_t}(\nabla_{\tau_t} v \cdot \tau_t) + R,$
其中 $A_{com}(f) := -\Delta_{\Gamma_t} \nabla_{n_t} f$ 为作用于内部超双曲系统的算子，区别于不可压缩情形中的椭圆算子 $A_{in}$ 。这一结构使得能够追踪自由表面的几何演化，并在能量估计中保持对称性。
角点区域的奇异性分析与正则性控制：针对角域上的椭圆边值问题，作者详细分析了解在接触点附近的渐近行为。解可分解为正则部分与奇异部分，奇异部分形如 $r^{\pi/\theta_i} \phi(\theta)$ 。通过选取 $\epsilon \in (0, \min\{1, \pi/\theta_i - 2\})$ ，证明了在锐角条件下可获得 $H^{3+\epsilon}$ 正则性。特别地，利用 Helmholtz 分解 $v = \nabla V + \nabla^\perp U$ 将旋度 $\omega$ 的控制转化为输运方程，避免了线性化过程中的导数损失。同时，通过精细的迹定理（Lemma 2.3-2.5）处理了角点处的边界积分，确保了计算在奇点附近的合法性。
迭代构造与无导数损失策略：为解决近似解构造中的正则性损失问题，作者设计了一套新的迭代方案。首先通过演化方程确定密度 $\rho_{k+1}$ 和自由面 $\Gamma_{k+1}_t$ ，再利用非退化椭圆系统求解中间变量，最后通过除-旋系统（Div-curl system）唯一确定速度场 $v_{k+1}$ 。该过程不依赖抛物正则化或 Nash-Moser 迭代，而是在低正则 Sobolev 空间下直接构造，证明了序列的收敛性且无导数损失。关键步骤在于利用几何结构将边界积分中的奇异项显式分离，并通过引入辅助变量 $\dot{v}_{k+1}$ 控制旋度演化。

这些技术共同构成了一个自洽的分析框架，成功克服了可压缩性与几何奇异性带来的双重挑战。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上实现了显著推进：

研究对象与适用范围的突破：既往关于接触点问题的工作主要集中在不可压缩流体（如 Alazard et al. [4], Ming [27], Ming and Wang [29-31]）或光滑边界的可压缩流体（如 Lindblad [21], Luo [24]）。本文首次将研究对象扩展至可压缩流体与接触点动力学的耦合系统。此前，没有任何文献处理过可压缩欧拉方程在自由边界涉及接触点的情形，这填补了理论空白。
假设条件的放宽与深化：与 Kinsey and Wu [19] 或 Wu [42] 处理的无表面张力、尖顶角情形不同，本文考虑了表面张力效应，并允许接触角在 $(0, \pi/2)$ 范围内变化（锐角）。更重要的是，本文明确刻画了接触点处的耗散机制（见公式 (1.10)），这是由可压缩性导致的独特现象（Remark 1.2），而不可压缩情形下的耗散仅源于摩擦或几何效应。相比之下，Alazard et al. [4] 仅在垂直壁面（ $\theta=\pi/2$ ）情形下讨论了 Taylor 符号条件不退化的情况，而本文处理了更一般的锐角情形，并揭示了角度越小正则性越高的规律。
方法论的革新：经典可压缩自由边界问题常采用 Lagrangian 坐标或 Nash-Moser 迭代（如 Lindblad [21], Trakhinin [38]），但这些方法在处理角点奇异性时面临巨大困难，往往导致正则性损失。本文采用的几何结构方法（Geometric Structure Approach）源自 Shatah and Zeng [35]，但进行了关键性推广：将原本用于不可压缩流体的椭圆算子 $A_{in}$ 替换为适用于可压缩流体的双曲型算子 $A_{com}$ ，从而解决了声波与角点奇异性的耦合难题。这种方法不仅避免了导数损失，还为一般自由边界问题提供了更灵活的低正则性分析框架。
代价与局限：本文的结果依赖于锐角条件 $\theta_i \in (0, \pi/2)$ 以保证 $H^{3+\epsilon}$ 正则性，当角度接近 $\pi$ 时奇异指数 $\pi/\theta_i$ 趋近于 1，可能导致正则性丧失。此外，接触点处的耗散项依赖于特定的非退化初始条件（1.7），这在物理上对应于密度梯度与重力方向的特定匹配。尽管本文建立了局部适定性，但尚未探讨大时间行为或全局稳定性，这仍是未来研究方向。总体而言，本文在保持数学严谨性的同时，极大地拓展了可压缩水波理论的适用范围。

4. 具有条件稳定性的抛物型反源问题的无伴随积分反馈方法

An adjoint-free integral feedback method for a parabolic inverse source problem with conditional stability

作者: Sedar Ngoma

类别: math.NA, math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16133v1

摘要翻译: 我们研究了一个带有加性源项 $$f(t)+ζ(t,x)$$ 的线性非自治抛物方程的反源问题，其中未知分量仅依赖于时间，并通过解的积分观测进行恢复。在将问题归约为一个等价的线性反问题后，我们建立了 Tikhonov 正则化解的存在性与唯一性，并导出了一阶最优性条件。我们表明前向算子在时间上具有 Volterra 表示，从而在 $H^{-1}(0,T)$ 中给出弱范数稳定性估计；在源项的先验 $$H^r(0,T)$$ 界下，给出 $$L^2(0,T)$$ 中的条件 Hölder 稳定性估计。受此 Volterra 结构的启发，我们引入了一种无伴随积分反馈方法，该方法仅利用前向求解即可重构源项。通过分析反馈迭代的良定性、不动点性质、精确数据的收敛性以及关于噪声数据的有限迭代稳定性，我们对反馈迭代进行了分析。进一步表明，在适当的依赖噪声的停止规则下，该方法构成了一种迭代正则化方案。针对光滑源和分段常数源的数值实验，辅以时间正则化和自动参数选择，证明了在存在噪声的情况下能够实现准确且稳定的重构。

Comment: 35 pages, 4 figures, 2 tables

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

1. 研究对象和背景

本文研究一类线性非自治抛物型方程的逆源问题。考虑定义在有界光滑区域 $\Omega \subset \mathbb{R}^d$ ( $d \geq 1$ ) 和时间区间 $$(0, T)$$ 上的偏微分方程：

\partial_t u - \nabla \cdot (a(t, x)\nabla u) + b(t, x) \cdot \nabla u + c(t, x)u = f(t) + \zeta(t, x), \quad (t, x) \in (0, T) \times \Omega,

其中 $$a(t, x)$$ 为对称一致椭圆矩阵， $$b(t, x)$$ 和 $$c(t, x)$$ 为有界系数。源项由两部分组成：未知的仅依赖于时间的分量 $$f(t)$$ 和已知的时空依赖分量 $\zeta(t, x)$ 。问题设定还包括初始条件 $$u(0, x) = u_0(x)$$ 和齐次或给定的 Neumann 边界条件。

核心目标是利用对解的全域积分观测数据恢复未知时间源函数 $$f(t)$$ 。具体而言，已知观测关系为：

\int_{\Omega} u(t, x) \, dx = \mu(t), \quad t \in [0, T].

该问题在 Hadamard 意义下是不适定的，因为从积分观测反演源函数是一个典型的病态逆问题，微小的数据扰动会导致解的巨大偏差。尽管已有文献（如 Ngoma [19]）在 Hölder 空间框架下研究了此类问题，并证明了在特定条件下（如零阶系数为零）的问题适定性，但现有理论主要局限于经典解，且往往假设算子是自洽的（即系数不显含时间）。在实际应用中，数据通常带有噪声，且物理系数可能不具备高阶正则性，这促使我们需要在弱形式（Hilbert 空间）框架下建立理论，并开发数值稳定的重构算法。

本文试图填补的空白在于：针对非自治抛物算子，系统性地构建一个基于变分法的弱解存在唯一性理论，深入分析前向算子的 Volterra 结构及其稳定性性质（包括弱范数稳定性和先验条件下的 Hölder 稳定性），并在此基础上提出一种无需计算伴随方程的“积分反馈”迭代方法，以实现对未知源的稳定、高效重构。

2. 主要定理或主要结果

本文的主要数学成果可概括为以下几个方面：

Tikhonov 正则化解的存在性与唯一性：在 $$L^2(0, T)$$ 空间中定义 Tikhonov 泛函 $J_\lambda(f) = \|F(f) - h_\delta\|^2_{L^2(0,T)} + \lambda\|f\|^2_{L^2(0,T)}$ 。证明了对于任意正则化参数 $\lambda > 0$ 和适配后的数据 $h_\delta \in L^2(0, T)$ ，该泛函存在唯一的极小元 $f_\lambda \in L^2(0, T)$ 。此结论基于变分法中的直接方法，利用了 $$L^2$$ 空间的自反性、泛函的强制性（coercivity）、弱下半连续性以及严格凸性。
一阶最优性条件：推导了极小元满足的 Euler-Lagrange 方程，即一阶最优性条件：
$\langle F^* (F(f_\lambda) - h_\delta) + \lambda f_\lambda, v \rangle_{L^2(0,T)} = 0, \quad \forall v \in L^2(0, T),$
其中 $$F^*$$ 为前向算子的伴随算子。这一条件虽然未用于后续的无伴随算法，但提供了正则化解的变分刻画，是连接伴随方法与反馈方法的理论桥梁。
前向算子的 Volterra 结构与稳定性估计：证明了前向算子 $F: L^2(0, T) \to L^2(0, T)$ （定义为 $F(f)(t) = \int_\Omega u_f(t, x) dx$ ）在时间上具有 Volterra 积分算子的结构。基于此结构，作者建立了以下稳定性结果：
- 弱范数稳定性：在 $H^{-1}(0, T)$ 范数下，存在稳定性估计，表明逆问题的不适定程度被控制在较弱拓扑中。
- 条件 Hölder 稳定性：若对未知源施加先验约束 $f \in H^r(0, T)$ 且 $\|f\|_{H^r} \leq M$ ，则在 $$L^2(0, T)$$ 范数下可获得 Hölder 型稳定性估计，即存在常数 $C, \nu \in (0, 1)$ 使得： $\|f_1 - f_2\|_{L^2(0,T)} \leq C \|F(f_1) - F(f_2)\|_{L^2(0,T)}^\nu.$ 这一结果揭示了平滑机制（smoothing mechanism）如何导致逆问题的病态性，同时也指出了通过先验正则化获得稳定解的可能性。
无伴随积分反馈方法的收敛性与正则化性质：提出了如下迭代格式：
$f^{(k+1)} = f^{(k)} - \frac{1}{\alpha} \Phi(F(f^{(k)}) - h_\delta),$
其中 $\Phi$ 为由非线性反馈函数诱导的 Nemytskii 算子， $\alpha > 0$ 为反馈参数。证明了该迭代在适当假设下是良定义的，精确数据情形下收敛于真解，且在噪声数据下具有有限步稳定性。进一步证明，配合适当的基于噪声水平的停止规则（stopping rule），该方法构成了一种迭代正则化方案。

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

本文的核心技术障碍在于处理非自治抛物算子下的逆问题不适定性，同时避免传统伴随方法高昂的计算成本（需反复求解伴随 PDE）。作者的真正创新点在于利用前向算子的内在结构（Volterra 性质）设计纯前向迭代的反馈机制。

1. 变分框架与能量估计：证明的基础建立在严格的弱解理论之上。作者首先利用标准抛物型方程理论确立了正问题在 $L^2(0, T; H^1(\Omega))$ 空间中的适定性。随后，通过对控制方程乘以解本身并积分，结合散度定理和 Young 不等式，导出了关键的能量不等式（见文中公式 4.13-4.15）：

\frac{d}{dt}\|u_f(t)\|^2_{L^2(\Omega)} \leq C_1 \|u_f(t)\|^2_{L^2(\Omega)} + C_2 |f(t)|^2.

利用 Gronwall 引理，作者得出了前向算子 $$F$$ 的有界性估计（公式 4.16），这是证明 Tikhonov 泛函 coercivity 的关键。这一系列估计不仅保证了优化问题的适定性，也为后续分析算子谱性质奠定了基础。

2. Volterra 结构的识别与稳定性分析：这是本文理论深度的体现。作者指出，由于抛物型方程的时间演化特性，解 $$u(t)$$ 仅依赖于 $s \leq t$ 时刻的源项 $$f(s)$$ ，因此前向算子 $$F$$ 本质上是一个 Volterra 算子。这种结构意味着 $$F$$ 具有紧性特征，从而导致逆问题不适定。作者利用这一结构，将逆问题的稳定性分析与 Volterra 算子的性质联系起来。通过引入先验空间 $$H^r(0, T)$$ 的约束，作者成功地将稳定性从较弱的 $H^{-1}$ 范数提升到 $$L^2$$ 范数下的 Hölder 稳定性。这一分析框架避免了直接使用谱分解等强假设，适用于更广泛的非自治系数情形。

3. 无伴随积分反馈算法的设计与分析：传统 PDE 约束优化通常依赖伴随梯度法，需要求解伴随方程 $$F^*$$ 。本文提出的核心创新是构造了一个仅依赖前向求解的迭代格式：

f^{(k+1)} = f^{(k)} - \frac{1}{\alpha} \Phi(F(f^{(k)}) - h_\delta).

这里的关键工具是 Nemytskii 算子 $\Phi$ 的性质分析。作者没有简单地将此视为数值实验技巧，而是将其作为不动点迭代进行严格的数学分析。证明过程涉及：

良定义性与 Lipschitz 连续性：验证反馈算子在 $$L^2$$ 空间中的良好行为。
不动点性质：证明精确解是该迭代的不动点。
误差传播分析：通过分析残差 $r_k = F(f^{(k)}) - h_\delta$ 的演化，建立有限步稳定性界限。
迭代正则化理论：引入基于噪声水平 $\delta$ 的停止准则（如 Morozov 不一致性原理的变体），证明当迭代次数 $k(\delta)$ 随噪声减小而趋于无穷时，近似解 $f^{k(\delta)}$ 收敛到真实解。这一部分将数值算法的提升到了理论分析的高度，填补了无伴随方法在理论保证方面的空白。

4. 数值实现的辅助技术：为了抑制数值振荡，作者在反馈迭代中引入了二阶时间正则化项，并结合曲率准则自动选择反馈参数 $\alpha$ 和正则化参数。这些技术细节虽然属于数值实现，但其设计依据正是前述的理论分析（如 Hölder 稳定性要求的先验约束）。

4. 与之前工作的比较

本文与既有 PDE/分析理论文献相比，在多个维度上实现了显著推进：

1. 问题设定的推广：从自洽到非自治：早期工作（如 Ngoma [19]）及多数经典文献主要集中在自洽（autonomous）抛物方程或系数不显含时间的情况。在这些设定下，算子谱理论更为成熟，稳定性分析相对直接。本文明确处理了非自治情形，即系数 $$a(t, x), b(t, x), c(t, x)$$ 均显含时间 $$t$$ 。这意味着不能简单利用半群理论或分离变量法，必须依赖变分法和能量估计。作者成功地将分析框架扩展到了非自治算子，这在物理建模中更具普适性（例如热传导系数随温度或时间变化的场景）。

2. 函数空间与正则性假设的放宽：从 Hölder 空间到 Hilbert 空间：既往研究多基于经典解理论，要求解和系数具有较高的正则性（如 Hölder 连续），以便应用 Schauder 估计等工具。然而，实际测量数据往往粗糙且带噪，系数也可能仅具 $L^\infty$ 正则性。本文转向弱解框架，在 Hilbert 空间 $L^2(0, T; H^1(\Omega))$ 中讨论问题。这不仅降低了对系数的光滑性要求，也使得处理带噪声数据成为可能。特别是，本文在 $$L^2$$ 空间中建立了 Tikhonov 正则化的存在唯一性，并给出了条件 Hölder 稳定性，这是对经典 Hölder 空间理论的补充而非替代，更贴合现代反问题研究的趋势。

3. 算法范式的转变：从伴随驱动到无伴随反馈：现有的数值重构方法大多基于伴随梯度下降法（Adjoint-based optimization），这需要每步迭代求解一个伴随 PDE，计算成本极高，尤其在大规模三维问题中。Ngoma [19] 使用了配置法结合 Tikhonov 正则化，但仍依赖于全局优化框架。本文提出的“无伴随积分反馈方法”彻底摒弃了伴随方程的求解，仅利用前向求解器更新源函数。这一突破不仅降低了计算复杂度（避免了双 PDE 求解循环），更重要的是提供了一种新的算法范式：将逆问题转化为一个带有特定结构的不动点迭代问题。文中明确指出，这种方法在处理非自治问题时，避免了伴随方程系数复杂变化带来的额外困难。

4. 稳定性分析的精细化：虽然许多文献讨论了逆问题的不适定性，但本文明确区分了弱范数（ $H^{-1}$ ）稳定性和条件 Hölder 稳定性。这种分层分析揭示了 Volterra 算子的平滑效应如何限制了解的不稳定性范围，并为迭代正则化提供了理论依据。相比之下，之前的工作往往笼统地讨论“稳定性”或仅在特定光滑假设下给出 Lipschitz 估计，缺乏对噪声环境下有限步稳定性的定量描述。

代价与局限：本文的方法并非完美无缺。首先，无伴随方法依赖于特定的反馈参数 $\alpha$ 和非线性函数 $\Phi$ 的选择，其收敛速度可能不如精心调优的伴随梯度法快。其次，条件 Hölder 稳定性依赖于先验信息 $f \in H^r$ ，如果先验假设过强或不准确，稳定性将无法保证。此外，虽然数值实验展示了效果，但对于极度不规则或高度振荡的源函数，该方法的表现仍需更多验证。尽管如此，本文为非自治抛物逆源问题提供了一个兼具理论严谨性和计算可行性的新路径。

5. 具有孤立边界奇点的临界半线性椭圆方程 II

The Critical Semilinear Elliptic Equation with Isolated Boundary Singularities II

作者: Hua-Yang Wang, Jingang Xiong

类别: math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16106v1

摘要翻译: 延续第二作者（2017）的工作，我们研究带有孤立边界奇点和零狄利克雷边界条件的半空间中的索伯列夫临界半线性椭圆方程。本文解决了该设定下的两个开放性问题：del~Pino–Musso–Pacard (2007) 提出的关于存在德劳内型对数周期解的问题，以及 Bidaut-Véron–Ponce–Véron (2007) 提出的关于奇异解的渐近分类问题。我们建立了一个全局正对数周期解分支，沿该分支在唯一确定的周期处发生爆破。我们还构造了相应的集中族并证明了其局部唯一性。因此，预期的平稳渐近分类失效，且在整个半空间中无法成立普适的临界尺度不变上界。这一行为与 Caffarelli–Gidas–Spruck (1989) 的经典内部奇点理论形成鲜明对比。

Comment: 52 pages

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

1. 研究对象和背景

本文研究的是在 $N \ge 3$ 维上半空间 $\mathbb{R}^N_+ = \{x = (x', x_N) \in \mathbb{R}^N : x_N > 0\}$ 中，带有原点处孤立边界奇点且满足零狄利克雷边界条件的索博列夫临界半线性椭圆方程的正解问题。具体而言，研究对象为以下边值问题：

\begin{cases} -\Delta u = u^p & \text{在 } \mathbb{R}^N_+ \text{ 中}, \\ u = 0 & \text{在 } \partial \mathbb{R}^N_+ \setminus \{0\} \text{ 上}, \end{cases}

其中指数 $$p$$ 取为索博列夫临界指数 $p = p_S(N) = \frac{N+2}{N-2}$ 。该问题在几何分析中对应于寻找具有常正标量曲率的共形平坦度量，且在原点处存在非可去奇点。

对于次临界情形 ( $$1 < q < p$$ )，Bidaut-Véron, Ponce 和 Véron [3, 4] 已经建立了完整的理论框架。他们证明了当 $$q$$ 处于特定范围时，存在唯一的正稳态奇异解 $\phi_q$ （即与径向距离无关的解），并给出了奇异解的分类：任何非可去奇点的解在趋于原点时都渐近收敛于该稳态解 $\phi_q$ ，且满足尺度不变的先验估计 $u(x) \le C |x|^{-\frac{2}{q-1}}$ 。这一结论依赖于圆柱能量（cylindrical energy）关于对数半径变量 $t = -\log|x|$ 的单调性。

然而，在索博列夫临界指数 $$q=p$$ 处，情况发生了根本性的变化。通过 Emden-Fowler 变换 $u(x) = |x|^{-a} v(t, z)$ （其中 $a = \frac{N-2}{2}$ , $t = -\log|x|$ , $$z = x/|x|$$ ），原方程转化为定义在圆柱面 $\mathbb{R} \times S^{N-1}_+$ 上的方程：

\begin{cases} -\left(\partial_{tt} + \Delta_{S^{N-1}} - a^2\right)v = v^p & \text{在 } \mathbb{R} \times S^{N-1}_+ \text{ 中}, \\ v = 0 & \text{在 } \mathbb{R} \times \partial S^{N-1}_+ \text{ 上}. \end{cases}

这里的关键障碍在于，临界情形下的一阶导数项消失，导致对应的圆柱能量不再是单调的，而是守恒的。这使得传统的单调性论证无法排除非稳态解的存在。特别是，del Pino, Musso 和 Pacard [16] 提出了两个开放性问题：是否存在类似于 Fowler 解的对数周期解（Delaunay-type log-periodic solutions）？以及临界情形下的奇异解是否仍遵循 Bidaut-Véron 等人建立的“稳态渐近分类”？

本文旨在填补这一空白。作者指出，由于能量守恒，系统可能存在非平凡的周期轨道，这将直接破坏原有的渐近分类猜想（Problem A）。此外，这种振荡行为可能导致解在原点附近的增长率超过标准的尺度不变上界，从而否定全局一致的上界估计（Problem C）。这与 Caffarelli-Gidas-Spruck [8] 在内部奇点情形下证明的刚性结果形成了鲜明对比，后者表明所有非可去奇点解均渐近于 Fowler 解且受统一常数控制。本文的研究核心在于揭示边界奇点在临界情形下的丰富动力学结构及其对经典理论的颠覆。

2. 主要定理或主要结果

本文主要建立了三个核心定理，彻底解决了上述开放问题，并揭示了临界边界奇点的新现象：

定理 1.1 (局部分岔与对数周期解的存在性)：固定 $0 < \alpha < \min\{1, p-1\}$ ，存在 $\varepsilon_0 > 0$ 和一条光滑曲线 $(k_\varepsilon, w_\varepsilon) \in (0, \infty) \times C^{2,\alpha}(T \times S^{N-1}_+)$ ，其中 $T = \mathbb{R}/(2\pi\mathbb{Z})$ 。该曲线从第一轴向共振点 $(k^*, \phi)$ 分岔而出，其中 $k^* = \sqrt{-\lambda_1}$ ， $\lambda_1$ 是线性化算子 $A_\phi = -\Delta_{S^{N-1}} + a^2 - p\phi^{p-1}$ 在轴对称类中的唯一负特征值。展开式为：

w_\varepsilon(s, z) = \phi(z) + \varepsilon \psi_1(z) \cos s + O(\varepsilon^2),

其中 $\psi_1$ 是对应的归一化主特征函数。对于 $0 < |\varepsilon| < \varepsilon_0$ ，解 $w_\varepsilon$ 是非稳态的，且具有最小周期 $2\pi$ 。这直接构造了 del Pino-Musso-Pacard 猜想中的对数周期解族。对应的原始解 $u_\varepsilon$ 满足离散自相似性： $u_\varepsilon(e^{-T}x) = e^{aT}u_\varepsilon(x)$ ，其中 $T = 2\pi/k_\varepsilon$ 。

定理 1.2 (全局分支与爆破周期)：局部分支包含在一个闭连通集 $\mathcal{C}$ 中。该集合具有以下性质：

每个非稳态解 $(k, w) \in \mathcal{C}$ 都是“鼓包”（bumpy）的，即存在平移使得 $$w$$ 关于 $$s=0$$ 偶对称且在 $$(0, T/2)$$ 上严格递减。
存在仅依赖于 $$N$$ 的正常数 $T_-, T_+, \Lambda_N$ ，使得对所有非稳态解，其周期 $T = 2\pi/k$ 满足 $T_- \le T \le T_+$ ，且能量有界： $\int_{\Sigma_T} v^{p+1} dV \le \Lambda_N$ 。
集合 $\mathcal{C}$ 是非紧致的。若序列 $(k_n, w_n) \in \mathcal{C}$ 离开任意紧集，则其周期 $T_n \to T^*$ ，且 $L^\infty$ 范数爆破 $\|w_n\|_{L^\infty} \to \infty$ 。这里 $$T^*$$ 是由级数定义的唯一点： $M_N(T) = -2^{2-N} + 2 \sum_{\ell=1}^\infty e^{-a\ell T} [(1-e^{-\ell T})^{2-N} - (1+e^{-\ell T})^{2-N}] = 0,$ 且 $$T^*$$ 是 $$M_N(T)$$ 的唯一零点。同时，重整化后的解在 $L^{p+1}$ 范数上收敛于欧氏空间中的标准气泡（bubble）积分： $\int_{\Sigma_{T_n}} v_n^{p+1} dV \to \int_{\mathbb{R}^N} U^{p+1} dy$ 。

定理 1.3 (集中族的构造与局部唯一性)：存在 $\mu_0 > 0$ 和一个参数族 $(T(\mu), v_\mu)$ ， $0 < \mu < \mu_0$ ，由中心化和缩放参数 $\mu$ 标记。当 $\mu \downarrow 0$ 时，周期 $T(\mu) \to T^*$ 且 $L^\infty$ 范数爆破。该定理进一步证明了在极限周期 $$T^*$$ 附近的局部唯一性：任何以 $$P=(0, e_N)$$ 为中心、周期趋于 $$T^*$$ 且最大值趋于无穷的正解序列，最终必然落入该构造的分支中。

推论 (渐近分类与上界的失效)：基于上述定理，作者否定了 Problem A 和 Problem C。具体而言：

渐近分类失效：存在非稳态的对数周期解，它们不收敛于稳态解 $\phi$ ，因此不存在统一的稳态渐近分类。
无全局尺度不变上界：由于解可以在有限周期 $$T^*$$ 处发生集中爆破，且爆破点处的值可以任意大，导致不存在常数 $$C=C(N)$$ 使得对所有解都有 $u(x) \le C x_N |x|^{-N/2}$ 成立。这意味着 Bidaut-Véron 等人在次临界或超临界情形下建立的锥内估计无法推广到整个上半空间。

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

证明的核心技术障碍在于处理临界情形下的能量守恒导致的非紧性以及缺乏相平面约简。作者采用了一套精密的分析框架，结合了分岔理论、全局延拓技术和精细的爆破分析。

1. 线性化分析与 Morse 指标计算：首先，作者深入研究了稳态解 $\phi$ 附近的线性化算子 $A_\phi = -\Delta_{S^{N-1}} + a^2 - p\phi^{p-1}$ 。通过变分法证明了 $\phi$ 作为约束极小值点，其 Morse 指标为 1。关键在于排除了零特征值的存在性（Lemma 2.1）。作者利用移动平面法（moving-plane method）证明了 $\phi$ 的轴对称性和单调性，进而将问题转化为一维奇异 Sturm-Liouville 问题。通过构造辅助函数 $Y = -\phi'/\phi$ 和 $$S = Y'$$ ，并利用积分因子公式分析了 Wronskian 的性质，严格证明了轴向核空间为零。这一步骤确立了第一共振频率 $k^* = \sqrt{-\lambda_1}$ 的唯一性，为 Crandall-Rabinowitz 分岔定理的应用奠定了基础。

2. Crandall-Rabinowitz 分岔定理的应用：在确认了 Fredholm 指标为零且横截性条件满足后（Lemma 2.3），作者应用 Crandall-Rabinowitz 定理构建了局部的对数周期解分支。通过引入轴向偶对称类（axial-even class）消除了时间平移的简并性，并在 Hölder 空间中进行了高阶正则性分析。这一步不仅构造了解，还确定了其最低周期为 $2\pi$ 的结构特征。

3. 全局延拓与 Leray-Schauder 度：为了将局部分支延拓至全局，作者采用了 Leray-Schauder 度理论。通过将方程重写为不动点形式 $w = (-L_{g_k})^{-1}(w^p - \phi^p)$ ，并利用了紧算子族的一致估计。关键的一步是证明了分支不会经过其他共振点 $$k_m = k^*/m$$ ( $m \ge 2$ )。作者通过能量估计和最大原理，证明了分支上的解必须保持“鼓包”（bumpy）结构，即解在半个周期内严格单调递减。这一几何约束排除了高阶共振的可能性，确保了分支的连通性和拓扑结构的稳定性。

4. 有限周期的爆破分析与精确周期确定：这是本文最具创新性的部分。当分支趋于非紧致时，解可能趋向于 $k \to 0$ （周期发散）或 $\|w\|_{L^\infty} \to \infty$ （爆破）。作者通过精细的能量估计（Proposition 4.1）排除了周期发散的情形，证明周期必须有界。随后，针对爆破情形，作者引入了“单气泡”（one-bubble）分析。通过对解进行中心化和平移，将问题转化为在圆柱面上逼近标准欧氏气泡 $$U$$ 的过程。通过计算特定的级数 $$M_N(T)$$ ，作者精确地定位了爆破发生的唯一周期 $$T^*$$ 。这个 $$T^*$$ 是一个超越数，由几何参数 $$N$$ 和临界指数决定，反映了圆柱几何与欧氏气泡之间的共振关系。最后，利用 Lyapunov-Schmidt 约简和逆映射定理，作者在 $$T^*$$ 附近构造了唯一的集中族，并证明了其局部唯一性。

5. 反对称性与唯一性论证：在整个过程中，作者充分利用了轴对称性（axial symmetry）来降维。虽然原方程在球面上是偏微分方程，但通过限制在轴对称类中，并结合 Sturm 振荡理论和 Lagrange 恒等式，成功地将复杂的谱问题简化为可处理的常微分方程性质分析。这种策略避免了直接处理高维球面上的复杂频谱结构。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上实现了显著的推进和突破：

1. 与 Caffarelli-Gidas-Spruck (1989) 的内部奇点理论对比： Caffarelli, Gidas 和 Spruck [8] 在内部奇点情形下证明了著名的刚性结果：任何非可去奇点的正解都渐近于 Fowler 解（Fowler asymptotics），且满足统一的上界估计 $u(x) \le C |x|^{-a}$ 。他们的证明依赖于径向对称性将方程约化为常微分方程，从而可以通过相平面分析完全分类解的行为。然而，本文指出，在边界奇点情形下，尽管解具有轴对称性，但由于球面变量的存在，无法进行类似的有限维相平面约简。更关键的是，临界边界情形打破了内部的刚性：本文构造的非稳态对数周期解表明，解不再收敛于单一的稳态 profile，而是呈现出周期振荡行为。这种“柔性”是内部理论所不具备的，直接导致了渐近分类的失效。

2. 与 Bidaut-Véron-Ponce-Véron (2007) 的次临界及超临界理论的对比： Bidaut-Véron 等人在次临界 ( $$q < p$$ ) 和超临界 ( $$p < q < q_S(N-1)$$ ) 情形下建立了完整的奇点分类理论。他们证明了在适当的尺度不变估计假设下，所有解都渐近于唯一的稳态解 $\phi_q$ 。在次临界情形，能量单调性保证了这种收敛；而在超临界情形，虽然需要额外假设上界，但分类结论依然成立。本文的工作揭示了临界情形 ( $$q=p$$ ) 是该理论的“临界断裂点”。由于能量守恒，单调性丧失，导致稳态解不再是吸引子，取而代之的是周期轨道。这不仅推翻了 Problem A 中的猜想，也表明 Bidaut-Véron 等人的分类定理在临界指数处无法直接推广。此外，本文还否定了 Problem C 中关于全局一致上界的猜想，指出在临界情形下，解的爆破行为可能发生在任意接近原点的区域，且没有统一的锥状限制，这与他们在超临界情形下建立的锥内估计形成了鲜明反差。

3. 与 del Pino-Musso-Pacard (2007) 的未解决问题对比： del Pino, Musso 和 Pacard 在 [16] 中明确提出了关于临界边界奇点是否存在 Delaunay 型对数周期解的问题（Problem B），并指出了其与内部问题的差异。本文通过严格的数学推导，不仅正面回答了这个问题（肯定存在），而且给出了详尽的解集结构描述。作者证明了这些解构成了一个全局连通的分支，并且精确刻画了分支端点的爆破行为（唯一周期 $$T^*$$ ）。这一成果将之前的启发式猜测和存在性猜想提升到了定量的、结构化的理论高度。

4. 代价与局限：本文的推进是以牺牲“普适性”为代价的。在内部奇点情形，所有解的行为是统一的（收敛于 Fowler 解）；而在本文研究的边界临界情形，解的行为变得极其丰富且不稳定，依赖于具体的初始数据或参数，导致无法建立统一的渐近分类或全局上界。此外，本文的结果目前仅限于轴对称类（axial class），虽然作者引用了 Xiong [42] 证明了所有正解均为轴对称，但这一依赖关系意味着如果未来发现非轴对称解的存在，本文的部分结论可能需要修正。尽管如此，本文成功地打开了理解临界边界奇点动力学的大门，指出了该领域与经典内部理论的本质区别。

6. 二维Hénon问题中的新变号解

New sign-changing solutions to the Hénon problem in dimension $$2$$

作者: Francesca Gladiali, Pingping Yang

类别: math.AP

发布日期: 2026-08-16

arXiv ID: 2608.15638v1

摘要翻译: 在本文中，我们构造了Hénon问题

Comment: 无

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

1. 研究对象和背景

本文研究的是二维有界光滑区域 $\Omega$ （包含原点）上的Hénon方程：

\begin{cases} -\Delta u = |x|^\alpha |u|^{p-1}u, & \text{in } \Omega, \\ u = 0, & \text{on } \partial\Omega, \end{cases}

其中参数 $\alpha > 0$ ，指数 $p \to +\infty$ 。该方程由M. Hénon提出，用于描述旋转恒星结构，是偏微分方程分析理论中关于加权非线性椭圆问题的经典模型。

在二维情形下，与高维（ $N \ge 3$ ）不同，Hénon问题的解在 $p \to +\infty$ 时表现出独特的性质：解是一致有界的，且不会发生边界集中现象（boundary concentration）。对于正解，已有文献（如Esposito, Pistoia, Wei [21]）建立了当 $p \to \infty$ 时，解在原点附近形成 $$k$$ 个对称尖峰的集中行为，且 $k < \alpha + 1$ 。然而，关于变号解（nodal solutions）的理论远不如正解完善。现有的主要结果集中在径向解上，例如Amadori和Gladiali [6] 证明了存在具有任意数量节点区域的径向解，并计算了其Morse指标；而针对非径向变号解的研究，主要集中在单位球域 $$B_1$$ 上，且多局限于“拟径向”（quasiradial）解或特定对称空间中的最小能量解。

尽管Lane-Emden问题（即 $\alpha=0$ 的特例）中已经构造出了具有簇状结构（cluster structures）的变号解（见[11]），其渐近形态由一个中心正气泡和围绕它的 $$k$$ 个对称负气泡组成，但将这一机制推广到具有权重项 $|x|^\alpha$ 的Hénon问题上仍面临巨大挑战。特别是，权重项改变了原点处的线性化算子谱性质，导致正解的气泡形状不再是标准的Liouville气泡，而是Hénon型气泡。此外，如何在一般满足一定对称性的区域（而非仅限于单位球）上构造此类非径向、具有复杂簇状结构的变号解，仍是未完全解决的空白。本文旨在填补这一空白，通过引入新的对称性约束和精细的渐近展开，证明在二维情形下存在一类新的非径向簇状变号解，并精确刻画其浓度速率和渐近形态。

2. 主要定理或主要结果

作者在假设区域 $\Omega \subset \mathbb{R}^2$ 为光滑有界域、包含原点、关于 $$x_1$$ 轴反射对称（(1.2)式）且关于原点在角度 $2\pi/k$ 下旋转对称（(1.3)式）的前提下，证明了如下主要结论：

定理 1.1: 设 $\alpha > 0$ ，整数 $$k$$ 满足条件

\frac{3(2+\alpha)}{2} < k \le \left\lceil \frac{2+\alpha}{2}t \right\rceil - 1,

其中常数 $t \approx 5.1869$ 是方程 $2\sqrt{e} \log t + t = 0$ 的唯一正解。则存在 $$p_k > 1$$ ，使得当 $$p > p_k$$ 时，问题 (1.1) 在 $$k$$ -对称函数空间 $H^1_{0,k}(\Omega)$ 中存在一个变号解 $u_{k,p}$ ，其满足以下渐近展开：

u_{k,p} = \tau_p (P U_1 - \eta_p P U_2) + o\left(\frac{1}{p}\right), \quad \text{在 } L^\infty_{\text{loc}}(\Omega \setminus \{0\}) \text{ 中},

其中 $$P$$ 表示向 $H^1_0(\Omega)$ 的正交投影，核心气泡项定义为：

u_1(x) := \log \frac{2(2+\alpha)^2 \delta^{2+\alpha}_p}{(|x|^{2+\alpha} + \delta^{2+\alpha}_p)^2}, \quad \text{满足 } -\Delta U_1 = |x|^\alpha e^{U_1} \text{ in } \mathbb{R}^2,

u_2(x) := \log \frac{8k^2 \beta^{2k}_p}{(|x^k - \rho^k_p|^2 + \beta^{2k}_p)^2}, \quad \text{满足 } -\Delta U_2 = |x|^{2k-2} e^{U_2} \text{ in } \mathbb{R}^2.

这里 $$x^k$$ 表示复数坐标的 $$k$$ 次幂对应的向量。随着 $p \to +\infty$ ，所有参数满足以下尺度关系：

\tau_p \sim A_k p^{-1}, \quad \eta_p \sim \eta_{\infty, k} := \frac{2(2+\alpha)}{2k - 2 - \alpha},

\delta_p \sim B_k e^{-(b_k + o(1))p}, \quad \beta_p \sim C_k e^{-(c_k + o(1))p}, \quad \rho_p \sim R_k e^{-(r_k + o(1))p},

其中 $$A_k, B_k, C_k, R_k$$ 为仅依赖于 $$k$$ 和 $\alpha$ 的正常数，具体的指数系数为：

b_k = \frac{1}{2(2+\alpha)} \left( 1 + \frac{16k(2+\alpha)}{(2k+2+\alpha)^2} \log \frac{2k-2-\alpha}{2(2+\alpha)} \right),

c_k = \frac{1}{4k} \left( 1 + \frac{2(2k-2-\alpha)^2}{(2k+2+\alpha)^2} \log \frac{2k-2-\alpha}{2(2+\alpha)} \right),

r_k = \frac{2(2k-2-\alpha)}{(2k+2+\alpha)^2} \log \frac{2k-2-\alpha}{2(2+\alpha)}.

此外，解的能量泛函极限行为为：

p \int_\Omega |x|^\alpha u_+(x)^{p+1} dx \to 4(2+\alpha)\pi e \left( \frac{2k-2-\alpha}{2(2+\alpha)} \right)^{\frac{16k(2+\alpha)}{(2k+2+\alpha)^2}},

p \int_\Omega |x|^\alpha u_-(x)^{p+1} dx \to 8k\pi e \left( \frac{2(2+\alpha)}{2k-2-\alpha} \right)^{\frac{2(2k-2-\alpha)^2}{(2k+2+\alpha)^2}}.

关键推论（Remark 1.2）指出，由于参数满足 $0 < \delta_p < \beta_p < \rho_p$ ，这意味着正气泡在原点的浓度速度最快，负气泡的浓度速度次之，而负气泡中心趋近原点的漂移速度最慢。这种层级分离的结构构成了所谓的“簇状”（clustered）解。

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

本文的核心技术障碍在于处理二维Hénon方程中权重项 $|x|^\alpha$ 导致的奇异性和非标准气泡结构，以及构造非径向变号解所需的精细参数匹配。作者的创新点主要体现在以下几个方面：

混合气泡构造与高阶修正：传统的集中现象分析通常使用单一的标准Liouville气泡。本文为了构造变号解，巧妙地组合了两种不同类型的基解：
- 中心正气泡 $$U_1$$ ：这是Hénon方程特有的径向解，形式为 $U_1(x) \sim \log \frac{\delta^{2+\alpha}}{(|x|^{2+\alpha} + \delta^{2+\alpha})^2}$ ，它解决了原点附近的加权非线性项 $|x|^\alpha e^u$ 。
- 周围负气泡 $$U_2$$ ：这些气泡的形式类似于Lane-Emden问题中的标准气泡，但在复平面坐标下表现为 $U_2(x) \sim \log \frac{\beta^{2k}}{|x^k - \rho^k|^2 + \beta^{2k}}$ ，它们满足 $-\Delta U_2 = |x|^{2k-2} e^{U_2}$ 。为了控制误差，作者并未止步于一阶近似，而是引入了二阶和三阶修正项 $$V_1, W_1$$ （对应正气泡）和 $$V_2, W_2$$ （对应负气泡）。这些修正项是全局定义的辅助函数，满足特定的线性化方程（如 $\Delta V_\alpha + |x|^\alpha e^{U_\alpha} V_\alpha = \dots$ ），并具有对数增长行为（ $V \sim C \log|x|$ ）。这种高阶展开确保了在远离峰值的区域误差足够小，从而能够应用不动点定理。
对称性空间的选取与降维：作者没有直接在全空间寻找解，而是利用区域的几何对称性，将搜索空间限制在 $$k$$ -对称子空间 $H^1_{0,k}(\Omega)$ 中。该空间要求函数在旋转 $2\pi/k$ 下不变。这一策略极大地简化了问题：首先，它保证了构造出的解具有预期的对称结构；其次，它将原本复杂的无限维优化问题转化为有限维约化问题。通过在该子空间中寻找临界点，作者成功避免了非对称扰动带来的复杂性。
隐函数定理与势函数极小化：证明过程遵循标准的Lyapunov-Schmidt约化框架，但细节处理极为精细。首先，将解写为 $u = \Upsilon + \phi$ ，其中 $\Upsilon$ 是精心构造的近似解， $\phi$ 是小扰动。线性化算子 $L = \Delta + p|x|^\alpha |\Upsilon|^{p-1}$ 在垂直于核空间的方向上是可逆的（Section 4）。通过求解 $L\phi = -E$ （ $$E$$ 为误差项），可以将原方程约化为一个有限维问题：寻找参数 $(\tau, \delta, \beta, \rho, \eta)$ 使得投影后的误差为零。作者进一步证明，这个有限维问题等价于寻找一个单变量实值函数（关于参数 $\eta$ ）的临界点。通过详细的能量估计和渐近分析，作者证明了该函数在特定区间内存在唯一的极小值点，从而确定了所有参数的具体渐近行为。
加权范数与误差控制：由于解在不同区域（原点附近、负气泡附近、远处）的行为差异巨大，作者使用了分段定义的加权范数来统一控制误差。特别是在处理 $$U_1$$ 和 $$U_2$$ 相互作用产生的交叉项时，利用了参数间的尺度分离（ $\delta_p \ll \beta_p \ll \rho_p$ ）来保证各项的可控性。这种对参数依赖关系的精确推导（Lemma 3.1）是证明成立的关键。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上取得了显著推进：

从正解到变号簇状解的跨越：之前的研究（如[21], [37]）主要关注正解的集中行为，证明了正解可以形成 $$k$$ 个对称尖峰。虽然[6]讨论了变号解的存在性，但其构造的多是基于径向对称性或简单的“塔状”结构（tower of bubbles），即所有峰值都在原点交替出现。本文首次构造了具有空间分离结构的变号解：一个位于原点的正气泡被一圈（ $$k$$ 个）位于外围的负气泡包围。这种结构更接近于Lane-Emden问题（ $\alpha=0$ ）中的簇状解（[11]），但在Hénon问题中，由于权重的影响，正气泡的形状发生了根本改变，而负气泡的分布几何主要由对称性决定。这填补了二维Hénon问题中关于非径向、非原点集中变号解的理论空白。
适用域的扩展：以往关于非径向变号解的多数结果（如[6], [27]）严格依赖于单位球域 $$B_1$$ 的旋转对称性。本文虽然也要求区域 $\Omega$ 具有 $$k$$ -旋转对称性和反射对称性，但这比单纯的球形更宽泛。只要区域满足这些对称性条件（例如某些星形域或多边形域），结论依然成立。这表明所构造的簇状结构并非球域的独有现象，而是对称性本身的产物。
参数范围的界定与深化：文献[6]证明了存在至少 $\lceil \frac{2+\alpha}{2}t \rceil - 1$ 个非径向解，但未明确其具体的簇状形态。本文通过精确的参数计算，给出了存在此类簇状解的具体 $$k$$ 值范围： $\frac{3(2+\alpha)}{2} < k \le \lceil \frac{2+\alpha}{2}t \rceil - 1$ 。特别地，文中指出了下界 $k > \frac{3(2+\alpha)}{2}$ 的来源，这与确保负气泡浓度速度慢于正气泡、从而形成稳定簇状结构的物理直觉一致。同时，文章探讨了该下界是否仅为技术限制还是本质限制，暗示可能存在更小 $$k$$ 值的解，这为未来研究指明了方向。
渐近行为的精细化描述：相比于前人工作仅给出解的存在性定性描述，本文提供了详尽的渐近公式。特别是明确了正负气泡浓度速率的层级关系（ $\delta_p \ll \beta_p \ll \rho_p$ ）以及能量泛函的精确极限值（公式1.9）。这种定量精度使得该结果不仅是一个存在性定理，更是一个强有力的工具，可用于后续研究解的稳定性、Morse指标计算以及 bifurcation 分析。相比之下，早期关于Lane-Emden问题的类似工作[11]虽然也有类似结构，但本文成功地将这些技巧移植到了具有权重奇异性的高难度Hénon方程中，展示了更强的分析鲁棒性。
代价与局限：本文的代价是必须假设区域具有特定的离散对称性（ $$k$$ -旋转对称）。目前尚不清楚在非对称的一般区域中是否存在类似的簇状解，这也是作者指出的开放问题。此外，构造过程高度依赖于 $p \to \infty$ 的渐近分析，对于固定的 $$p$$ 或较小的 $$p$$ ，解的性质仍需进一步探索。尽管如此，本文在二维Hénon问题变号解理论中的地位，正如其在Lane-Emden问题中[11]的地位一样，起到了承前启后的关键作用。

7. 可压缩欧拉方程剪切层的长波谱不稳定性

Long-Wave Spectral Instability of Shear Layers for the Compressible Euler Equations

作者: Chao Wang, Yuxi Wang, Wenzhi Wu, Zhifei Zhang

类别: math.AP

发布日期: 2026-08-16

arXiv ID: 2608.15730v1

摘要翻译: 我们研究二维可压缩欧拉方程在光滑单调剪切层周围的长波谱不稳定性。通过长波展开构造可压缩Rayleigh方程的衰减半直线解，并推导匹配Wronskian的二阶展开式。对于任意固定的马赫数 $$m>0$$ ，我们证明了对于足够小的波数存在不稳定模态。当 $$m< anh Y$$ 时，该结论对一类剖面成立，且 $$c_i$$ 随 $$α o0$$ 趋于正常数。在 $$m= anh Y$$ 处，剖面 $$U_s(Y)= anh Y$$ 具有一个不稳定的模态，其中 $$c o0$$ 且 $$c_i$$ 的量级为 $α^{1/3}$ 。当 $$m> anh Y$$ 时，同一剖面仍保持不稳定，其中 $$c o c_*(m) anh Y$$ 且 $$c_i>0$$ 的量级为 $$α$$ 。相应的时间增长率分别为 $$α$$ 、 $α^{4/3}$ 和 $$α^2$$ 量级，表明在 $$m= anh Y$$ 处长波不稳定性标度发生了变化。在零厚度极限下，超临界不稳定特征值趋近实轴，这与 extcite{CS1,CS2} 中关于超音速可压缩涡片的稳定性结果一致。

Comment: 35 pages

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

1. 研究对象和背景

本文研究二维可压缩欧拉方程在光滑单调剪切层背景下的长波谱不稳定性问题。流体运动由无量纲可压缩欧拉方程组描述：

\begin{cases} \partial_t \rho_e + \nabla \cdot (\rho_e U_e) = 0, \\ \rho_e (\partial_t U_e + U_e \cdot \nabla U_e) + \nabla P(\rho_e) = 0, \end{cases}

其中 $\rho_e$ 为密度， $$U_e = (U_e^1, U_e^2)$$ 为速度向量， $P(\rho_e)$ 为压强函数。我们考虑连接无穷远处两个均匀状态的剪切层剖面 $(\rho_s, U_s)(Y) = (1, U_s(Y), 0)$ ，满足边界条件 $\lim_{Y \to \pm \infty} U_s(Y) = \pm 1$ 。该剖面是上述方程组的稳态解。

为了分析该剪切层剖面的（不）稳定性，我们将线性化可压缩欧拉方程组在稳态解附近展开。线性化系统在整体域 $(X, Y) \in \mathbb{R} \times \mathbb{R}$ 上表示为：

\begin{cases} \partial_t \varrho + U_s \partial_X \varrho + \nabla \cdot U = 0, \\ \partial_t U + U_s \partial_X U + \partial_Y U_s(Y) (U^2, 0) + m^{-2} \nabla \varrho = 0, \end{cases}

其中 $\varrho = \rho_e - 1$ ， $$U = (U^1, U^2) = U_e - (U_s(Y), 0)$$ ， $m = \frac{1}{\sqrt{P'(1)}}$ 为马赫数。我们的目标是寻找形如 $(\varrho, U^1, U^2)(t, X, Y) = (\rho, u, v)(Y)e^{i\alpha(X-ct)} + c.c.$ 的不稳定解，其中 $\alpha > 0$ 为波数， $c = c_r + i c_i \in \mathbb{C}$ 为复相速度，且要求 $$c_i > 0$$ 以保证指数增长。

该问题的核心在于构造对应的瑞利型方程（Rayleigh-type equation）的解。通过引入广义流函数 $\phi(Y)$ ，系统可降阶为关于 $\phi$ 的二阶常微分方程：

(U_s - c) \Lambda \phi - \partial_Y (A^{-1} \partial_Y U_s) \phi = 0,

其中 $\Lambda = \partial_Y (A^{-1} \partial_Y) - \alpha^2$ ， $$A(Y) = 1 - m^2(U_s(Y) - c)^2$$ ，并需满足边界条件 $\phi(\pm \infty) = 0$ 。一旦找到非平凡解 $\phi$ ，即可还原出密度扰动 $\rho$ 和速度分量。

尽管对于不可压缩流动，Rayleigh 判据和 Fjørtoft 判据等经典理论已较为完善，且 Lin 等人建立了充分不稳定准则，但对于可压缩流动，尤其是全平面光滑剪切层在临界马赫数 $m=\sqrt{2}$ 附近的长波不稳定性机制尚不明确。Coulombel 和 Secchi 等人在涡片（vortex sheet）的稳定性方面取得了重要进展，证明了 $m>\sqrt{2}$ 时超音速可压缩涡片的线性稳定性，但光滑剪切层与间断涡片在物理机制和数学结构上存在显著差异。特别是当马赫数接近或超过临界值 $\sqrt{2}$ 时，特征方程中的系数可能退化，导致传统的摄动方法失效。本文旨在填补这一空白，证明对于任意固定的马赫数 $$m>0$$ ，只要波数 $\alpha$ 足够小，均存在不稳定的模态，并精确刻画了不同马赫数区间内增长率 $$c_i$$ 对波数 $\alpha$ 的依赖关系及标度律的变化。

2. 主要定理或主要结果

作者证明了对于任意固定的马赫数 $$m > 0$$ ，存在一个足够小的波数阈值 $\alpha_0 > 0$ ，使得当 $0 < \alpha < \alpha_0$ 时，线性化系统存在具有正虚部相速度 $$c_i > 0$$ 的不稳定本征模。具体结论根据马赫数 $$m$$ 的范围分为三种情况：

亚临界情形 ( $0 < m < \sqrt{2}$ )：对于满足结构假设（单调性、指数衰减等）的一类剖面 $$U_s(Y)$$ ，存在不稳定模态，其相速度 $c(\alpha)$ 满足：
$c(\alpha) = i\sigma(m) + O(\alpha),$
其中 $\sigma(m) = \left( \frac{\sqrt{4m^2+1} - (m^2+1)}{m^2} \right)^{1/2} > 0$ 。此时实部 $c_r(\alpha) = O(\alpha)$ ，虚部 $c_i(\alpha) = \sigma(m) + O(\alpha)$ 趋于正常数。时间增长率 $\alpha c_i \sim \alpha$ 。
临界情形 ( $m = \sqrt{2}$ )：对于特定剖面 $U_s(Y) = \tanh Y$ ，存在不稳定模态，其相速度满足：
$c(\alpha) = e^{i\pi/6} \left( \frac{\alpha}{3} \right)^{1/3} + o(\alpha^{1/3}).$
具体地， $c_r(\alpha) = \frac{\sqrt{3}}{2} (\frac{\alpha}{3})^{1/3} + o(\alpha^{1/3})$ ， $c_i(\alpha) = \frac{1}{2} (\frac{\alpha}{3})^{1/3} + o(\alpha^{1/3}) > 0$ 。此时时间增长率 $\alpha c_i \sim \alpha^{4/3}$ 。
超临界情形 ( $m > \sqrt{2}$ )：对于同一剖面 $U_s(Y) = \tanh Y$ ，仍存在不稳定模态，其相速度满足：
$c(\alpha) = c_*(m) + K_m \alpha + o(\alpha),$
其中 $c_*(m) = \left( \frac{m^2+1 - \sqrt{4m^2+1}}{m^2} \right)^{1/2} \in (0, 1)$ 为实数根， $$K_m$$ 为复常数且 $\text{Im } K_m > 0$ 。具体表达式为：
$K_m = \frac{1-c_*^2}{8c_*^2\sqrt{4m^2+1}} \left[ 2\pi c_* + i \left( 4 + 2c_* \log \frac{1+c_*}{1-c_*} \right) \right].$
此时 $c_r(\alpha) = c_*(m) + \text{Re } K_m \alpha + o(\alpha)$ ， $c_i(\alpha) = \text{Im } K_m \alpha + o(\alpha) > 0$ 。时间增长率 $\alpha c_i \sim \alpha^2$ 。

此外，在零厚度极限下（即剖面趋近于涡片），当 $m > \sqrt{2}$ 时，超临界不稳定本征值的虚部趋于零，这与 Coulombel 和 Secchi 关于超音速可压缩涡片稳定性的结果一致；而当 $m < \sqrt{2}$ 时，不稳定性依然存在。

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

证明的核心障碍在于处理可压缩瑞利方程在不同马赫数 regime 下的奇异行为，特别是当 $m \geq \sqrt{2}$ 时，特征系数 $$A(Y) = 1 - m^2(U_s - c)^2$$ 可能在某些点接近零，导致代数奇异性。作者并未直接求解全直线上的方程，而是采用了以下创新的分析框架：

半直线解的构造与长波展开：作者将问题分解为两个半直线问题，分别构造在 $Y \to +\infty$ 和 $Y \to -\infty$ 处衰减的解 $\phi_+$ 和 $\phi_-$ 。为了获得精确的渐近行为，作者引入了长波展开形式：
$\phi_\pm(Y) = e^{\mp \beta_\pm Y} \chi_\pm(Y), \quad \chi_\pm = \sum_{n=0}^\infty (\pm \beta_\pm)^n \chi_{\pm, n},$
其中 $\beta_\pm = \alpha \gamma_\pm$ ， $\gamma_\pm = A_{\pm \infty}^{1/2}$ 。通过递归求解算子 $\text{Ray}_0[\chi_{\pm, n}] = F_{\pm, n}$ ，作者成功构建了级数解。关键的技术突破在于证明了该级数在 $m \geq \sqrt{2}$ 时的收敛性，即便此时 $$A(Y)$$ 可能接近零。作者利用 $\text{Ray}_0$ 算子的发散结构（divergence structure），将明显的代数奇异性转化为对数损失（logarithmic losses），并建立了统一估计：
$\| \chi_\pm(Y) - (U_s(Y) - c) \|_{Y\theta, \pm} \lesssim \alpha^\lambda,$
其中 $\lambda = 1$ 当 $m < \sqrt{2}$ ，而 $\lambda = |\log c_i|^2$ 当 $m \geq \sqrt{2}$ 。这一估计保证了在临界马赫数附近展开的有效性。
Wronskian 匹配条件的二阶展开：全局解存在的充要条件是半直线解在 $$Y=0$$ 处的 Wronskian 为零，即色散关系 $W(\alpha, c) = 0$ 。作者推导了 $W(\alpha, c)$ 关于小参数 $\alpha$ 的二阶展开式：
$(1 - m^2 c^2)^{-1} W(\alpha, c) = \frac{4\alpha c P_m(c)}{D_m(c)} + \alpha^2 I_{m, U_s}(c) + R(\alpha, c),$
其中 $$P_m(c)$$ 仅依赖于剪切层两端的均匀状态（profile-independent），而 $I_{m, U_s}(c)$ 是包含剖面信息的二阶修正项。这一分解至关重要，因为它分离了通用部分和剖面依赖部分。
色散关系的渐近求解与收缩映射：针对不同马赫数范围，作者分析了多项式 $$P_m(c)$$ 的根结构：
- 当 $m < \sqrt{2}$ 时， $$P_m(c)$$ 在上半平面有一个简单虚根 $c=i\sigma(m)$ ，结合有界的二阶项，利用收缩映射原理得到 $c(\alpha) = i\sigma(m) + O(\alpha)$ 。
- 当 $m = \sqrt{2}$ 时，主项根退化为原点， $P_{\sqrt{2}}(c) \sim -6c^2$ ，主导平衡变为 $c^3 \sim i\alpha/3$ ，从而导出 $c \sim \alpha^{1/3}$ 的标度律。
- 当 $m > \sqrt{2}$ 时， $$P_m(c)$$ 有一个简单的实根 $$c_*(m)$$ ，此时一阶项为零，必须依靠二阶项 $\alpha^2 I_{m, U_s}(c)$ 的虚部符号来决定根是否移入上半平面。作者通过显式计算 $I_{m, \tanh}(c)$ 的极限值，证明了 $\text{Im } K_m > 0$ ，从而确认了弱不稳定性。

整个证明框架巧妙地结合了奇异摄动理论、ODE 的渐近分析和复变函数论中的零点存在性定理，克服了可压缩效应带来的系数退化难题。

4. 与之前工作的比较

本文与既有文献相比，在以下几个方面实现了显著的推进和区分：

从间断涡片到光滑剪切层的跨越： Coulombel 和 Secchi 在 [6, 7] 中主要研究了间断的可压缩涡片（vortex sheets）。他们证明了当 $m > \sqrt{2}$ 时，超音速可压缩涡片是线性稳定的（即不存在 $$c_i > 0$$ 的本征模）。然而，本文针对的是光滑剪切层（smooth shear layers）。结果显示，即使对于 $m > \sqrt{2}$ ，光滑剪切层依然表现出长波不稳定性，尽管这种不稳定性非常微弱（增长率仅为 $O(\alpha^2)$ ）。这揭示了光滑性与间断性在稳定性机制上的本质差异：光滑剪切层通过连续的速度梯度提供了额外的不稳定源，而间断涡片在超音速下则受到声波的抑制作用。此外，本文还指出，在零厚度极限下，超临界情形的不稳定本征值确实趋于实轴，与涡片稳定性结果相容，但在有限厚度下，光滑性保留了不稳定性。
临界马赫数 $m=\sqrt{2}$ 处的精细标度律：此前关于可压缩流动的稳定性研究往往关注 $m < \sqrt{2}$ 的亚临界区或 $m > \sqrt{2}$ 的超音速区，而对临界点 $m=\sqrt{2}$ 附近的过渡行为缺乏定量分析。本文填补了这一空白，精确导出了临界马赫数下相速度随波数的标度律 $c \sim \alpha^{1/3}$ ，以及相应的增长率 $\alpha c_i \sim \alpha^{4/3}$ 。这一标度律既不同于亚临界的 $O(\alpha)$ ，也不同于超临界的 $O(\alpha^2)$ ，标志着不稳定性机制在该点发生了质的改变（从代数主导变为立方根主导）。这一发现完善了可压缩剪切层稳定性理论的完整图景。
剖面依赖性与非普适性：在亚临界区 ( $m < \sqrt{2}$ )，不稳定性主要由端点状态决定，具有某种“普适性”；而在临界和超临界区，不稳定性强烈依赖于具体的剖面形状（如 $U_s(Y)=\tanh Y$ 的具体积分性质）。作者明确指出了 $I_{m, U_s}(c)$ 作为剖面依赖系数的关键作用，并给出了判断不稳定性是否存在的一般条件（ $\text{Re } I_{m, U_s}(c) < 0$ ）。这与 Lin 在不可压缩流动中的工作形成对比，后者更多关注 Rayleigh 判据的几何条件，而本文在处理可压缩问题时，深入到了具体的积分算子层面。
技术方法的改进：相比于以往处理类似奇异摄动问题的方法，本文发展了一套统一的半直线展开框架，成功处理了 $m \geq \sqrt{2}$ 时出现的对数奇异性。这种方法不仅适用于剪切层，也为研究其他涉及可压缩效应的边界层或自由界面问题提供了新的分析工具。特别是将 Wronskian 展开至二阶并分离通用项与修正项的策略，使得能够精确捕捉到超临界情形下微小的不稳定增长，这是传统能量估计或紧性方法难以达到的精度。

8. 偏微分方程中驻定周期图案尖点分岔的严格验证

Rigorous Validation of Cusp Bifurcations of Stationary Periodic Patterns in Partial Differential Equations

作者: Dominic Blanco, Jean-Philippe Lessard

类别: math.AP

发布日期: 2026-08-16

arXiv ID: 2608.15613v1

摘要翻译: 本文提出了一种计算机辅助框架，用于严格验证抛物型半线性偏微分方程中空间对称驻定周期图案的尖点分岔。我们的方法将先前在有限维情形下发展的尖点映射公式推广至无限维设定。我们将尖点条件表述为傅里叶系数希尔伯特空间上的零点查找问题，其非退化解对应于尖点分岔点。一个关键的技术要素是对所得序列空间中伴随乘法算子的细致处理，这比有限维情形更为复杂。从数值计算的近似解出发，我们发展了一个构造性的牛顿-坎托罗维奇论证，以证明存在靠近该近似值的尖点映射的零点及其局部唯一性。该解的非退化性直接导致三次正规形式系数 $$c$$ 非零。为了完成验证，我们通过适配问题对称结构的类热尔戈林估计严格界定了线性化算子的谱，确认恰好有一个特征值具有零实部。作为进一步推论，在尖点处严格正实部的特征值数量 $$k$$ 与 $$c$$ 的符号（两者均由该框架严格认证）共同决定了尖点区域内共存的三个解的稳定性结构：当 $$k=0$$ 且 $$c<0$$ 时出现双稳态，当 $$k=0$$ 且 $$c>0$$ 时出现单稳态，而当 $k

Comment: 无

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

1. 研究对象和背景

本文研究一类抛物型半线性偏微分方程（PDE）系统中，具有空间对称性的驻定周期模式的尖点（cusp）分岔现象。考虑如下形式的方程组：

\partial_t u = f(\lambda, u) \stackrel{\text{def}}{=} l_\lambda u + g(\lambda, u),

其中 $u = (u_1(x,t), \dots, u_p(x,t))$ 定义在周期性区域 $D \subset \mathbb{R}^m$ 上， $\lambda = (\lambda_1, \lambda_2) \in \mathbb{R}^2$ 为二维参数向量， $l_\lambda$ 是依赖于参数的线性微分算子， $g(\lambda, u)$ 为非线性项（可能包含空间导数）。驻定周期模式对应于时间无关解，即满足 $f(\lambda, u) = 0$ 的函数 $u \in H^{2r}_G(D)^p$ ，其中 $H^{2r}_G(D)$ 表示在有限离散对称群 $$G$$ 作用下不变的 Sobolev 空间。

尖点分岔是参数依赖动力系统全局分岔结构中的核心组织中心，属于余维二（codimension-two）分岔。在二维参数空间中，它标志着两条折叠（fold/saddle-node）分岔曲线的切向交汇与终止。其局部解结构具有普适性：在尖点邻域内，参数平面被两条折叠曲线划分为三个区域——含一个解的区域、含三个解的区域以及两条曲线本身（此处两个解碰撞）。这种结构直接导致了双稳态（bistability）或单稳态（monostability）现象，进而产生滞后回线（hysteresis loop），在生态模式形成、反应 - 扩散系统中的波钉扎（wave pinning）及化学自催化斑纹/条纹图案组织中具有广泛应用背景。

尽管尖点分岔的重要性已众所周知，但在无限维相空间中严格证明特定 PDE 系统存在尖点分岔（即数学上确证所有非简并条件）仍是一个开放难题。现有的计算机辅助证明框架多集中于有限维动力系统或一阶分岔（如鞍结分岔、Hopf 分岔等）。本文旨在填补这一空白，将有限维情形下发展的尖点映射（cusp map）形式体系推广至无限维 PDE 设定，通过构造性牛顿 - 坎托罗维奇（Newton-Kantorovich）论证，严格验证尖点分岔点的存在性、局部唯一性及其非简并性，并进一步确定分岔后共存解的稳定性结构。

2. 主要定理或主要结果

本文建立了基于傅里叶系数希尔伯特空间的零值寻找问题框架，成功证明了多个具体模型中尖点分岔的存在性与唯一性。主要结果可概括为以下定理形式：

一般性存在性定理（基于 Newton-Kantorovich 论证）：设 $\bar{x}$ 为陈氏映射 $F: X \to X$ （其中 $X = \ell^2_{G,\omega}$ 为加权对称序列空间）的一个数值近似解。若存在常数 $$Y_0, Z_1, Z_2$$ 满足特定的范数不等式约束（详见 Theorem 2.7），使得：
$Y_0 (1 - Z_1 \|A F'(\bar{x})^{-1}\|)^{-1} < \eta$
且
$Z_1 \|A F'(\bar{x})^{-1}\| < \frac{1}{2}, \quad Z_2 \|A F'(\bar{x})^{-1}\|^2 < \frac{1}{4},$
则存在唯一的真解 $\tilde{x}$ 位于以 $\bar{x}$ 为中心、半径为 $\rho = (1 - \sqrt{1 - 2Z_1\|AF'(\bar{x})^{-1}\|})/(Z_1\|AF'(\bar{x})^{-1}\|)$ 的闭球内，满足 $F(\tilde{x}) = 0$ 。该解的非简并性（即 $DF(\tilde{x})$ 可逆）直接蕴含了尖点分岔所需的三次正规形系数 $c \neq 0$ 及横截性条件 $\Delta \neq 0$ 。
稳定性结构判定（Lemma 1.2）：设 $$k$$ 为线性化算子 $D_u f(\lambda, u)$ 实部严格大于零的特征值个数（代数重数计数）， $$c$$ 为三次正规形系数。在尖点区域内，三个共存解 $\tilde{u}_1, \tilde{u}_2, \tilde{u}_3$ 的稳定性由 $(k, \text{sgn}(c))$ 完全决定：
- 若 $$k=0$$ 且 $$c<0$$ ：外两层解 $\tilde{u}_1, \tilde{u}_3$ 渐近稳定，中间解 $\tilde{u}_2$ 不稳定（双稳态）；
- 若 $$k=0$$ 且 $$c>0$$ ：仅中间解 $\tilde{u}_2$ 渐近稳定，外层解不稳定（单稳态）；
- 若 $k \geq 1$ ：三个解均不稳定。
具体应用实例（Theorems 3.6-3.8, 4.4, 4.5）：作者对以下五个独立场景给出了严格证明，所有场景中均证实 $$k=0$$ ：
- 1D Swift-Hohenberg 方程（ $$Z_2$$ 对称）： $$c < 0$$ ，呈现双稳态；
- 2D Swift-Hohenberg 方程（ $$D_4$$ 正方形对称）： $$c > 0$$ ，呈现单稳态；
- 2D Swift-Hohenberg 方程（ $$D_6$$ 六边形对称）： $$c > 0$$ ，呈现单稳态；
- 1D Gray-Scott 系统（ $$Z_2$$ 对称）： $$c > 0$$ ，呈现单稳态；
- 2D Gray-Scott 系统（ $$D_2$$ 矩形对称）：取参数 $\delta_1 \approx 0.08583846158364704, \delta_2 = 1, d = 1.5$ ，存在 $\tilde{\lambda} \approx (1, 4.183827930189122)$ 及解 $\tilde{u}$ ，半径误差界为 $3 \times 10^{-12}$ ， $$c > 0.1796$$ ，呈现单稳态。

所有结论均通过区间算术（interval arithmetic）严格封装了 $$c$$ 的符号及 $$k$$ 的计数，排除了数值舍入误差的影响。

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

本文的核心技术障碍在于如何将有限维动力系统的尖点分岔理论严谨地推广至无限维 Fourier 序列空间，同时处理对称性约束下的算子结构与谱估计问题。作者的创新点主要体现在以下三个相互关联的分析框架中：

对称性约化的傅里叶序列空间与伴随乘法算子：在无限维设定下，传统的矩阵转置不再适用，必须精确刻画序列空间中的伴随算子。由于引入了有限离散对称群 $$G$$ ，函数 $$u$$ 的傅里叶系数需满足代数约束 $u_n = \alpha_g(n) u_{\beta_g(n)}$ 。这导致 $$L^2_G(D)$$ 内积转化为带轨道权重（orbit-weighted）的序列内积 $(\cdot, \cdot)_2$ 。关键难点在于：线性化算子 $$D_u f$$ 的伴随 $D_u f^\dagger$ 在序列空间中不仅涉及共轭转置，还涉及一个复杂的“伴随乘法算子”（adjoint multiplication operator）。作者推导了该算子在 $\ell^2_{G,\omega}$ 上的显式公式（Section 2.3），指出其形式为 $(M_u w)_n = \sum_k \overline{u_{\beta_g(k)}} w_{\dots}$ 的卷积结构，并证明了其在加权空间中的有界性。这一结果是构建陈氏映射（cusp map）并验证 Fredholm 性质（如左零向量 $$w$$ 的存在性）的基础，也是有限维情形所未涉及的深层分析工具。
构造性 Newton-Kantorovich 论证与陈氏映射：作者将尖点分岔的所有代数条件（ $$f=0$$ , $$D_u f v = 0$$ , $D_u f^\dagger w = 0$ , $$b=0$$ , $c \neq 0$ 的隐式表达等）编码为一个增广系统 $$F(x)=0$$ ，其中 $x = (\lambda, u, v, w, h, z)$ 。不同于传统数值方法仅给出近似解，本文采用 Newton-Kantorovich 定理进行刚性验证。核心步骤包括：
- 构造近似逆算子 $A \approx [DF(\bar{x})]^{-1}$ ：利用线性化算子在高频模态下的对角占优性质（diagonal dominance），将 $$A$$ 分解为有限维块（数值求逆）与显式控制的无穷对角尾部；
- 计算残差界 $Y_0 = \|A F(\bar{x})\|$ 及 Lipschitz 常数 $$Z_1, Z_2$$ ：利用 $\ell^2_{G,\omega}$ 上的卷积不等式（Lemma 2.2, 2.4）控制非线性项的范数增长；
- 验证压缩条件：确保 $Z_1 \|A DF(\bar{x})\| < 1/2$ 等不等式成立，从而保证解的唯一性与存在性。该论证不仅确认了零点的存在，还通过 $DF(\tilde{x})$ 的可逆性自动保证了 $c \neq 0$ 和横截性条件，无需单独验证。
适应对称结构的 Gershgorin 谱封装：为了验证谱条件（即 $$D_u f$$ 恰有一个零特征值，其余特征值实部远离虚轴），作者发展了针对对称约束的 Gershgorin 圆盘定理。由于对称性耦合了不同傅里叶模态，线性化算子不再具有标准的卷积对角结构。作者将谱分析分解为三部分：
- 低频部分：通过显式计算有限维矩阵的包围盒（enclosure）；
- 高频部分：利用对角占优性质构造 Gershgorin 圆盘，证明其不接触虚轴；
- 中间过渡区：结合前两者的估计。此外，对于非自伴系统（如 Gray-Scott），该方法扩展至复平面上的圆盘估计，精确计数实部为正的特征值个数 $$k$$ 。这一过程完全通过区间算术实现，确保了 $$k$$ 值的数学确定性。

整个证明框架将解析分析与数值计算无缝结合，避免了传统摄动法中的小参数假设，实现了从数值近似到严格定理的逻辑闭环。

4. 与之前工作的比较

本文与既有文献相比，在理论深度、适用范围及技术细节上均有显著推进：

从有限维到无限维的维度跨越：此前关于尖点分岔的计算机辅助证明工作主要集中在有限维常微分方程（ODE）或低维流形上，特别是文献 [65] 建立了有限维情形的陈氏映射框架。然而，将此类框架直接移植到 PDE 面临本质困难：无限维算子的谱性质、伴随算子的复杂结构以及非线性项在函数空间中的正则性要求。本文首次将陈氏映射形式体系推广至无限维 Hilbert 空间，解决了序列空间中伴随乘法算子的精确刻画问题（这是有限维情形中矩阵转置所掩盖的深层结构）。这一突破使得尖点分岔的严格验证能够应用于具有空间延展性的 PDE 系统，如 Swift-Hohenberg 方程和 Gray-Scott 系统。
对称性处理的精细化：在周期性 PDE 中，平移不变性导致解通常不是孤立的，这使得零值寻找问题 ill-posed。虽然前人工作（如 [98, 99]）已利用对称性约化来消除平移简并，但本文进一步将这些技术整合进尖点分岔的验证流程中。具体而言，本文定义了带轨道权重的序列空间 $\ell^2_{G,\omega}$ ，并在此空间内重新推导了伴随算子公式。相比之下，以往工作要么忽略了对称性对伴随算子的影响，要么仅在自伴情形下讨论。本文的处理方式适用于更广泛的晶体学对称群（如 $$D_4, D_6$$ ），并能严格处理非自伴系统（如 Gray-Scott）中的谱问题。
谱估计与稳定性判据的完备性：现有文献多关注分岔点的存在性，而较少严格验证分岔后的稳定性结构。本文不仅证明了尖点分岔的存在，还通过改进的 Gershgorin 估计精确计算了正实部特征值个数 $$k$$ ，并结合 $$c$$ 的符号给出了完整的稳定性分类（双稳态 vs 单稳态）。在 [65] 等有限维工作中， $$k$$ 的计算相对简单，但在 PDE 中，由于连续谱的潜在干扰和高频模态的复杂性，这一任务极具挑战。本文通过构造对角占优的近似逆和分区间的谱包围，成功克服了这一障碍，这是以往工作未达到的精度与广度。
适用范围与代价权衡：本文的方法适用于五类不同的物理模型（1D/2D Swift-Hohenberg, 1D/2D Gray-Scott），涵盖了标量方程与系统方程，以及多种对称类型。代价在于计算复杂度随空间维数和对称群阶数增加而急剧上升，需要精心选择加权序列空间以控制卷积误差。尽管如此，作者展示了该方法在保持严格性的同时，仍能处理高维（2D）和复杂对称（六边形）情形，这在之前的计算机辅助证明中极为罕见。与纯数值模拟相比，本文提供了数学上无懈可击的保证；与纯解析摄动法相比，本文摆脱了小参数假设，适用于强非线性区域。

总体而言，本文不仅扩展了 [65] 的理论边界，还建立了一套通用的 PDE 尖点分岔验证范式，为后续研究更复杂的时空模式分岔奠定了坚实基础。

9. 具有不规则漂移的超布朗运动的紧支集性质

Compact Support Property of Super-Brownian Motion with Irregular Drift

作者: Leonid Mytnik, Johanna Weinberger

类别: math.PR, math.AP

发布日期: 2026-08-16

arXiv ID: 2608.15616v1

摘要翻译: 我们研究一维随机偏微分方程

Comment: 无

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

1. 研究对象和背景

本文研究一维随机偏微分方程（SPDE）：

d_t X_t(x) = \frac{1}{2}\Delta X_t(x) + b_1\mathbb{1}_{\{X_t(x)>0\}} + \sqrt{X_t(x)}\dot W(t,x),

其中 $$b_1 > 0$$ ， $\dot W$ 为时空白噪声，初始条件 $$X_0 = f$$ 为非负紧支连续函数。该方程描述了一类具有非零常数漂移项（当密度为正时）的超布朗运动（Super-Brownian Motion, SBM）。在经典的超布朗运动理论中，漂移项通常假设为线性或满足一定正则性条件，以便应用 Dawson-Girsanov 变换将问题转化为无漂移情形处理。然而，本文中的漂移项 $b_1\mathbb{1}_{\{x>0\}}$ 在原点处存在跳跃间断，属于高度不规则的漂移，不满足经典 Dawson-Girsanov 定理所需的绝对连续性或 Lipschitz 条件。

研究对象的核心在于证明该 SPDE 的唯一弱解 $$X_t$$ 具有紧支集性质（Compact Support Property, CSP），即若初始状态 $$f$$ 具有紧支集，则对于任意 $$t>0$$ ， $$X_t$$ 几乎必然仍具有紧支集。这一性质在物理上对应于粒子群扩散过程中不会瞬间“填充”整个空间，而是保持有限的传播速度。CSP 是分析超过程路径性质、接触概率及灭绝概率的基础。已有文献如 Iscoe [10] 通过倒向热方程的对偶性证明了经典 SBM 的 CSP；Shiga [24] 和 Englander & Pinsky [4,5] 将其推广至更广泛的系数范围；Jin & Perkowski [11] 处理了粗糙超布朗运动。然而，对于此类具有跳跃漂移的模型，传统基于测度值过程对偶性或 Girsanov 变换的方法失效。本文试图填补的空白正是针对这类非正则漂移项下超布朗运动的紧支集性质的严格证明，并进一步推导其总质量与平方 Bessel 过程的比较关系，从而得出正灭绝概率的结论。

2. 主要定理或主要结果

本文建立了以下核心结论：

解的层分解构造与收敛性（定理 2.2）：作者构造了一个序列 $(\bar{X}^n)_{n\ge 1}$ ，其中 $\bar{X}^n_t = \sum_{k=0}^n X^k_t$ 。初始项 $$X^0$$ 是无漂移的标准超布朗运动，后续项 $$X^k$$ ( $k\ge 1$ ) 是在前序层支集的正性区域上具有随机迁入率（immigration rate）的超布朗运动。具体地，迁入率由递归公式给出：
$g_k(t,x) = b_1 \mathbb{1}_{\{X^{k-1}_t(x)>0\}} \mathbb{1}_{\{\sum_{i=1}^{k-1} g_i(t,x)=0\}},$
且 $g_1(t,x) = b_1 \mathbb{1}_{\{X^0_t(x)>0\}}$ 。该序列在空间 $L^1(\Omega; C([0,T], (M^+_F, dbl)))$ 中依范数收敛到方程 (1) 的唯一弱解 $$X$$ 。
紧支集半径估计（命题 2.3）：设 $R_T = \sup\{|x|: \exists t\in[0,T], X_t(x)>0\}$ 为最大支撑半径。对于初始总质量 $m_f = \langle f, 1\rangle$ 和初始支撑半径 $$R_0$$ ，存在有限常数 $$A(T)$$ 使得期望满足：
$\mathbb{E}_f[R_T] \le A(T) \left( R_0 + \sqrt{R_0} + \beta(T, m_f) \right),$
其中 $\beta(T, m_f) \to 0$ 当 $m_f \downarrow 0$ 。这表明随着初始质量的减小，支撑半径的期望趋于初始半径。
紧支集性质成立（定理 2.4）：基于上述半径估计，证明了对于任意 $$T>0$$ 和紧支初始条件 $f \in C^+_c(\mathbb{R})$ ，解几乎必然保持紧支集：
$\mathbb{P}_f\left( X_t \in C^+_c(\mathbb{R}), \forall t \ge 0 \right) = 1.$
此外，该结果可推广至更一般的漂移项 $$h(u)$$ ，只要 $h(u) \le b_1 \mathbb{1}_{\{u>0\}}$ 且满足弱唯一性。
正灭绝概率（定理 2.6）：利用支撑半径估计将总质量 $\langle X_t, 1\rangle$ 与一个参数化的平方 Bessel 过程 $\tilde{z}_t$ 进行比较：
$d\tilde{z}_t = \delta dt + \sqrt{\tilde{z}_t} dB_t,$
其中 $\delta$ 的选择依赖于 $$b_1 R_T$$ 。分析表明，当时间 $$T$$ 、初始质量及初始支撑直径足够小时， $\tilde{z}_t$ 在时间 $$T$$ 前以正概率触及零点。由此导出原过程 $$X$$ 具有正灭绝概率：
$\mathbb{P}_f(\langle X_t, 1\rangle = 0) > 0, \quad \forall t > 0.$

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

本文的核心技术障碍在于漂移项 $b_1\mathbb{1}_{\{x>0\}}$ 的不规则性破坏了传统的 Girsanov 变换框架，使得无法直接利用已知结果。作者的创新在于构建了一种**层分解（Layer Decomposition）**策略，将复杂的非线性随机演化分解为一系列具有随机迁入率的超布朗运动的单调极限。

具体而言，证明框架分为三个关键步骤：

递归构造与随机迁入机制：作者没有尝试直接求解原方程，而是定义了一个迭代过程。第一层 $$X^0$$ 是标准超布朗运动。第二层 $$X^1$$ 的迁入率仅在 $$X^0$$ 的正支集上激活，且强度为 $$b_1$$ 。更一般地，第 $$k$$ 层 $$X^k$$ 的迁入率 $$g_k$$ 仅在前一层 $X^{k-1}$ 为正且之前所有层在该点均为零的区域生效。这种构造巧妙地将原本耦合的漂移项转化为随机的、逐层激活的源项。每一步的构造都依赖于前一层的过滤信息 $\mathcal{F}^{k-1}_\infty$ ，确保了因果性和可测性。这种分解本质上是将非线性的漂移效应“解耦”为一系列线性的（相对于当前层）随机过程叠加。
对偶方程的奇异性分析：为了控制每一层 $$X^k$$ 的支集扩张，作者深入分析了与超布朗运动对偶的 log-Laplace 方程（Log-Laplace Equation）：
$\partial_t V_t = \frac{1}{2}\Delta V_t - \frac{1}{2}V_t^2 + \psi,$
特别是边界迹问题（Boundary Trace Problems）。在处理无穷大源项（对应于支撑边界的爆炸性增长）时，作者引用并扩展了 Iscoe [10] 和 Evans [9] 的结果。通过构造特定的测试函数和势函数（如 Lemma 3.7 中的 $\psi_{m,r}$ ），导出了支撑边界附近解的精确衰减估计（如不等式 (28) 和 (29)）：
$\lim_{m\to\infty} V^{m,r}_t(x) \le c_{29}(r-|x|)t^{-3/2}\exp\left(-\frac{(r-|x|)^2}{2t}\right).$
这些估计量化了扩散前沿的传播速度，证明了即使存在不规则漂移，支集也不会无限迅速扩张。
单调收敛与误差控制：通过证明部分和序列 $\bar{X}^n$ 的单调性及在 $$L^1$$ 意义下的收敛性，作者确立了极限过程即为原方程的弱解。关键的技术细节在于证明每一层 $$X^k$$ 的支集半径增量是可求和的。利用上述对偶估计，作者控制了迁入项对支集扩张的贡献，从而导出了总支撑半径的期望界（命题 2.3）。最后，通过将总质量与平方 Bessel 过程进行鞅比较，利用 Bessel 过程在低维参数下的回归性质，严格证明了灭绝概率的正性。整个证明避免了直接处理不规则漂移带来的路径不可积性问题，转而通过级数展开和对偶分析间接控制全局行为。

4. 与之前工作的比较

本文与既有文献相比，在以下几个方面取得了显著推进：

漂移正则性的突破：以往关于超布朗运动带漂移的研究，如 Mytnik 等人在 [20] 中的工作，主要依赖 Dawson-Girsanov 变换。这要求漂移项 $$h(X_t)$$ 必须满足一定的正则性（通常是 Lipschitz 或满足线性增长条件），以保证测度变换的 Radon-Nikodym 导数存在且良好定义。本文处理的漂移项 $b_1\mathbb{1}_{\{x>0\}}$ 在原点处不连续，完全超出了 Dawson-Girsanov 定理的适用范围。这是本文最根本的理论贡献：首次在不使用 Girsanov 变换的情况下，证明了具有跳跃间断漂移的超布朗运动仍具有紧支集性质。
方法论的革新：Iscoe [10] 开创性地利用对偶热方程证明了经典 SBM 的 CSP，随后 Shiga [24] 和 Englander-Pinsky [4,5] 将其推广至更广泛的扩散系数。然而，这些方法在处理外部漂移源时，往往假设漂移是平滑的或具有特定结构。本文提出的层分解法是一种全新的构造性工具。它将复杂的随机 PDE 转化为一系列具有随机迁入率的超过程之和。这种方法不仅解决了 CSP 问题，还为研究其他具有奇异系数的 SPDE 提供了新的范式，即通过递归分解将非线性项“线性化”为可控的源项。
适用范围与代价：本文的结果适用于一类更广泛的漂移函数 $$h(u)$$ ，只要 $h(u) \le b_1\mathbb{1}_{\{u>0\}}$ 。这与 [20] 中处理的 $h(u) = \theta u - u^2$ 等光滑函数形成了鲜明对比。然而，这种推广是有代价的：证明过程变得极其复杂，需要精细地控制每一层的支集重叠和迁入率的时间依赖性，导致无法像经典情形那样得到简洁的路径模估计。此外，虽然证明了正灭绝概率，但具体的灭绝时间分布或临界参数并未像经典模型那样被完全刻画，目前仅给出了存在性结果。
与近期工作的对比：Jin & Perkowski [11] 研究了二维粗糙超布朗运动的 CSP，结合了测度值方法和内估计技巧。本文虽然也在处理不规则性，但侧重于漂移项的几何不连续性而非空间正则性不足。两者在技术路线上截然不同：前者依赖解析估计和正则化理论，后者依赖概率构造和对偶分析。本文的工作表明，即使在最原始的“开关型”漂移下，超布朗运动的物理直观（有限传播速度）依然成立，这巩固了该类模型在数学物理中的基础地位。

10. 扩散方程的随机团块方法

Random blob methods for diffusion

作者: Katy Craig, Claire Murphy

类别: math.NA, math.AP

发布日期: 2026-08-14

arXiv ID: 2608.14914v1

摘要翻译: 线性和非线性扩散方程出现在一系列具有数学意义的现象中，包括慢速和快速扩散、沙堆动力学、高度约束输运、二维Navier-Stokes方程以及概率测度采样的动力学。近年来，团块方法作为一种数值模拟这类偏微分方程（PDEs）的方法受到了广泛关注。为了解决团块方法中 $$O(N^2)$$ 的计算瓶颈，我们考虑了时空的随机离散化。我们将随机批处理方法与我们称之为随机多速率方法的新方法进行了比较。这两种方法都建立在优化文献中的经典随机方法之上，与随机梯度下降和随机坐标下降有许多相似之处。我们发现，对于线性和非线性扩散方程，在计算复杂度与精度之间的权衡中，随机多速率方法具有最佳性能。一方面，我们证明了随机多速率方法以 $$O(k Δt)$$ 的速率收敛到潜在的常微分方程（ODE）系统，这与前向欧拉法一致，并通过示例表明该速率在理论上是尖锐的。另一方面，当应用于由扩散PDEs的团块方法产生的ODE时，我们观察到随机多速率方法具有更快的收敛速率。最后，由于其能够模拟广泛的非线性扩散方程，包括高度约束输运和沙堆动力学，我们的方法成功捕捉到了那些仅有少量数值方法存在的PDE的关键特征。

Comment: 无

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

1. 研究对象和背景

本文的研究对象是广泛存在于数学物理中的线性与非线性扩散方程，其统一形式可表述为如下偏微分方程组：

\begin{cases} \partial_t \rho = \nabla \cdot (\nabla f^*(p) + \rho v), \\ p \in \partial f(\rho), \\ \rho(0) = \rho_0. \end{cases}

其中 $f: [0, \infty) \to \mathbb{R} \cup \{\infty\}$ 是一个固定的凸函数， $$f^*$$ 为其凸共轭， $\partial f(\rho)$ 表示 $$f$$ 在 $\rho$ 处的次微分， $$v$$ 为给定的速度场。该框架涵盖了多种重要的物理模型：当 $f(s) = s \ln s - s$ 时对应慢扩散（ $$m>1$$ ）与快扩散（ $$m<1$$ ）方程；当 $$f$$ 取特定分段形式时对应沙堆动力学（sandpile dynamics）；当 $$f$$ 为高度约束指示函数时对应高度受限输运问题；当 $$d=2$$ 且 $$f$$ 与 Biot-Savart 核相关时，则对应二维 Navier-Stokes 方程的粘滞形式。

这些方程的一个关键特性在于，若初始测度 $\rho_0$ 为有限支撑测度，其解 $\rho(t)$ 通常不再保持有限支撑性质。因此，直接应用粒子方法会失效。为了解决这一困难，Craig 等人 [14] 引入了正则化能量泛函 $F_\epsilon$ ，诱导出一类近似的偏微分方程：

\begin{cases} \partial_t \rho_\epsilon = \nabla \cdot (\rho_\epsilon (\nabla p_\epsilon + v)), \\ p_\epsilon = \phi_\epsilon * f'_\epsilon(\phi_\epsilon * \rho_\epsilon), \\ \rho_\epsilon(0) = \rho_{\epsilon, 0}. \end{cases}

在该正则化框架下，若初始条件为有限个粒子的加权和 $\rho_{\epsilon, 0} = \sum_{j=1}^N m_j \delta_{x_0^j}$ ，则解保持粒子结构 $\rho_\epsilon(t) = \sum_{j=1}^N m_j \delta_{x_j(t)}$ 。此时，求解 PDE 转化为求解如下的常微分方程组（ODE）：

\dot{x}_i(t) = -f''_\epsilon\left(\sum_{j=1}^N m_j \phi_\epsilon(x_i(t) - x_j(t))\right) \sum_{j=1}^N m_j \nabla \phi_\epsilon(x_i(t) - x_j(t)) - v(t, x_i(t)).

尽管“团块方法”（blob methods）在模拟此类非线性扩散方程方面具有显著优势，但其数值实现面临巨大的计算瓶颈。对于包含 $$N$$ 个粒子的系统，显式欧拉法（Forward Euler, FE）在每个时间步需要计算所有粒子对之间的相互作用，导致单步计算复杂度为 $$O(N^2)$$ 。若总时间步数为 $$M$$ ，则总计算成本高达 $$O(M N^2)$$ 。随着粒子数 $$N$$ 的增加，这成为了限制高精度模拟的主要障碍。

本文旨在填补的空白在于：如何在保持团块方法对复杂非线性扩散方程（特别是高度受限输运和沙堆动力学等缺乏现有数值方法的模型）的高保真度的同时，突破 $$O(N^2)$$ 的计算瓶颈。作者并未提出新的空间离散格式，而是专注于改进上述 ODE 系统的时间离散策略。通过引入随机化的时间步长机制，试图在计算效率与数值精度之间寻找更优的平衡点，特别是针对那些传统确定性多速率方法难以处理的动态尺度变化问题（即某些粒子运动快、某些粒子运动慢，且这种状态随时间动态转换）。

2. 主要定理或主要结果

本文的核心贡献在于提出并严格分析了两种随机离散化方法：随机批处理法（Random Batch, RB）与随机多速率法（Random Multirate, RM），并证明了后者在扩散方程背景下的优越性。

随机多速率法的收敛性定理：文章建立了随机多速率法（RM）相对于底层 ODE 系统（2.4）的强收敛性估计。设 $$G$$ 满足有界性 $||G(x, u)|| \le L_0$ 和 Lipschitz 连续性 $||G(x, u)-G(y, v)|| \le L_1 \max(||x-y||, ||u-v||)$ 。令 $$X_i(t)$$ 为精确 ODE 解， $$x_i^m$$ 为 RM 方法在第 $$m$$ 步的近似解，总时间 $T = M\Delta t$ 且 $$k$$ 为粗/细时间步长比。定理 2.1 指出：
$\max_{1 \le i \le N} ||X_i(T) - x_i^M|| \le \frac{L_0 e^{6 L_1 T}}{4} \left( L_1 (k\Delta t)^2 + k\Delta t \right).$
该结果表明，RM 方法的误差阶为 $O(k\Delta t)$ ，这与标准向前欧拉法（FE）使用步长 $k\Delta t$ 时的理论精度一致。文中通过具体算例证明，对于一般的 ODE 系统，这一收敛率是理论上最优的（sharp）。
扩散方程背景下的超收敛现象：尽管理论界限显示 RM 方法与粗步长的 FE 方法同阶，但数值实验表明，当该方法应用于由团块方法离散化得到的扩散型 ODE 系统时，其表现远优于理论预期。具体而言，在固定时间步长 $k\Delta t$ 下，RM 方法的实际误差显著小于对应的 FE 方法。特别是在慢扩散设置（slow diffusion）和高度受限输运（height-constrained transport）问题中，RM 方法展现出最佳的精度-效率权衡。
参数优化与性能对比：通过数值实验确定了最优参数配置：
- 对于 RB 方法，采用两批处理（ $$p=N/2$$ ）以平衡速度与精度。
- 对于 RM 方法，设定 $50\%$ 的粒子进行细时间步更新（ $$p=N/2$$ ），且时间步长比 $$k=2$$ 时效果最佳。在此配置下，RM 方法的计算复杂度为 $O(MN(p + \frac{N-p}{k}))$ ，即约为 $$O(MN)$$ 量级（常数因子取决于 $$k$$ 和 $$p/N$$ ），远低于 $$O(MN^2)$$ 。更重要的是，在相同的计算成本下，RM 方法的精度显著高于 RB 方法和标准 FE 方法。
适用范围验证：研究证实了该方法在处理非光滑内能密度（如高度受限输运和沙堆动力学）以及高维（ $$d=2$$ ）情形下的有效性，能够捕捉到这些复杂 PDE 的关键特征，而这些特征往往是现有数值方法难以处理的。

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

本文解决 $$O(N^2)$$ 瓶颈的核心技术障碍在于：传统的确定性多速率积分方法依赖于静态划分“快”与“慢”组件，而扩散系统中的粒子动力学行为是动态变化的——一个粒子可能在某一时刻快速移动，随后变慢，反之亦然。静态划分无法适应这种时变的尺度分离，导致效率低下或精度损失。

创新点一：随机多速率（RM）时间离散框架 作者提出的真正创新在于将随机性引入时间离散过程，而非空间交互。RM 方法在每一个时间步 $$m$$ 上：

随机选取 $$p$$ 个粒子作为“细步长”粒子（fine particles），执行一步大小为 $\Delta t$ 的演化。
剩余的 $$N-p$$ 个粒子作为“粗步长”粒子（coarse particles），执行一步大小为 $k\Delta t$ 的演化。
在细步长粒子的演化过程中，利用线性插值 $L(i, m, I_{fine}, I_{coarse})$ 来近似粗步长粒子在中间时刻的位置。这种机制本质上是一种随机时间离散，它避免了为每个粒子单独计算局部误差估计（这在大规模系统中计算不可行），而是通过随机重分配粒子角色来平均化误差。

核心工具与证明技术：

确定性误差分析框架：与文献 [24] 中基于期望（in-expectation）的随机批处理误差分析不同，本文针对 RM 方法建立了一个确定性的误差上界（Theorem 2.1）。证明过程不依赖于统计平均，而是对任意一次随机划分的轨迹进行控制。
离散 Gronwall 不等式：证明的核心步骤涉及构造误差递推关系。定义 $E_m^i = ||x_m^i - X_i(m\Delta t)||$ 为第 $$i$$ 个粒子的全局误差， $$H_m^i$$ 为由于粗步长插值引起的局部截断误差。通过 Lemma A.1（离散 Gronwall 型不等式）将局部误差累积放大效应控制住。
稳定性与一致性引理：Lemma A.2 和 A.3 分别提供了关于 ODE 向量场的 Lipschitz 连续性与数值格式的稳定性估计。特别是 Lemma 2.2 详细讨论了内能密度 $f_\epsilon$ 和速度场 $$v$$ 的局部 Lipschitz 条件，确保在紧集 $$B(0, R)$$ 内，修改后的向量场 $$G$$ 满足定理所需的有界性和 Lipschitz 条件。
线性插值误差控制：在处理粗步长粒子位置的不确定性时，利用线性插值公式 $$L$$ 将粗步长粒子的位置误差控制在 $O((k\Delta t)^2)$ 量级，从而保证了整体收敛率的稳定性。

分析框架的推进：该方法巧妙地避开了随机批处理法（RB）中因完全切断粒子间长程相互作用而导致的精度损失。RB 方法在每个时间步仅允许粒子与其所在批次内的粒子交互，这相当于人为切断了部分物理关联，导致在非均匀分布下产生较大偏差。而 RM 方法保留了全耦合的相互作用项（在计算力时仍对所有粒子求和，只是时间步长不同），仅通过调整时间步进频率来降低计算量。这种设计使得 RM 方法在保留物理机制完整性的同时，实现了类似 RB 方法的降阶计算复杂度。

4. 与之前工作的比较

本文与既有 PDE/分析理论文献的比较主要体现在以下几个方面：

与随机批处理法（RB Method）的比较：
- 假设与机制差异：Jin, Li, Liu [23, 24] 提出的 RB 方法是一种随机空间离散，即在每一步随机分组，仅计算组内相互作用。这种方法虽然将复杂度降至 $$O(MNp)$$ ，但牺牲了粒子间的长程关联信息。本文指出，在扩散方程背景下，这种空间截断往往导致显著的精度损失，即使在使用两批（ $$p=N/2$$ ）的情况下，其精度也往往不如 RM 方法。
- 收敛性理论：[24] 的结果是基于期望的强收敛性，且假设时间积分是精确的（仅考虑空间批处理的误差）。本文的 Theorem 2.1 则是确定性的误差界，同时考虑了时间离散误差和粗步长插值带来的轨迹近似误差，适用范围更广。
- 性能表现：数值实验表明，在相同的计算成本下，RM 方法在扩散型 ODE 上的收敛速度优于 RB 方法。特别是在高度受限输运等非光滑情形下，RB 方法容易丢失关键特征，而 RM 方法能有效捕捉。
与传统多速率方法（Multirate Methods）的比较：
- 静态 vs 动态划分：传统的多速率方法（如 [20, 28]）通常基于静态分区，预先识别系统中的“快”和“慢”分量。然而，在团块方法中，粒子的动力学行为随时间剧烈变化，静态划分不仅难以实施，而且会导致效率低下。Savcenco 等人 [33] 曾提出自适应方案，但需要频繁计算局部误差估计，计算代价过高。
- 创新路径：RM 方法通过随机重分配替代了复杂的误差估计过程。它不需要知道哪个粒子是“快”的，而是通过概率分布让粒子在不同时间尺度上自然流动，从而在无需额外计算开销的情况下实现了类似自适应的效果。
与团块方法（Blob Methods）文献的比较：
- 正则化策略的演进：早期工作（如 Craig et al. [14]）侧重于证明正则化方案 $f_\epsilon$ 收敛到原 PDE 的理论性质。本文在此基础上，进一步解决了数值实现的效率问题。文中提到，虽然 [14] 使用了外卷积 mollifier，但为了数值效率，本文去除了外层卷积，仅依赖 $f_\epsilon$ 自身的正则性，这在保证理论收敛的同时简化了计算。
- 适用范围扩展：现有的团块方法多集中于光滑的内能密度（如多孔介质方程）。本文明确展示了 RM 方法在处理非光滑内能（如高度受限输运 $f(s)=\infty$ for $$s>1$$ 和沙堆动力学）时的鲁棒性，填补了该类问题在高效数值算法方面的空白。
代价与局限：
- RM 方法并非在所有情况下都优于其他方法。在非扩散型的相互作用粒子系统（如公式 4.1 所示的特定吸引/排斥势）中，RM 方法的性能与标准 FE 方法相当，未显示出显著优势。这表明该方法的优势高度依赖于扩散算子带来的平滑效应和特定的相互作用结构。
- 尽管理论上是确定性的误差界，但在实际应用中，由于随机种子的选择，单次模拟的误差仍存在波动。文中建议在 $$d=2$$ 的高维情形下进行多次随机种子平均以获得稳定的误差估计，这在一定程度上增加了预处理成本，但相比 $$O(N^2)$$ 的原始成本仍可忽略不计。

总体而言，本文通过将随机性从空间域转移到时间域，成功地在保持物理机制完整性的前提下，突破了传统团块方法的计算瓶颈，为模拟广泛的非线性扩散方程提供了一套既高效又高精度的新范式。

11. 关于计算机辅助证明3D欧拉奇异性中非局部误差与轮廓残差误差区别的澄清

A clarification on the distinction between nonlocal error and profile residual error in a computer-assisted proof of 3D Euler singularity

作者: Jiajie Chen, Thomas Y. Hou

类别: math.AP

发布日期: 2026-08-15

arXiv ID: 2608.15174v1

摘要翻译: 在\cite{zhang2025dimension}中，作者讨论了自相似轮廓不存在性结果对Hou–Luo情形的适用性，并对\cite{ChenHou2023a,ChenHou2023b}中构造的近似自相似轮廓提出了质疑。在本注记中，我们澄清两个主要观点。首先，\cite{zhang2025dimension}中指出的远场差异是原始的泊松非局部误差，而非\cite{ChenHou2023a,ChenHou2023b}稳定性证明中所用的轮廓残差误差。与轮廓残差相关的非局部贡献包含乘以额外衰减因子的非局部速度误差，这些因子为远场估计提供了关键的有界小性。其次，我们考察了\cite[Remark 2.4]{zhang2025dimension}中用于支持\cite[Proposition 2.3(i)]{zhang2025dimension}假设适用性的基于图形的证据，并表明相同网格点数据并不支持该解释。因此，\cite{zhang2025dimension}中的比较和基于图形的证据既不能推翻\cite{ChenHou2023a,ChenHou2023b}中的残差估计，也不能证实所声称的角度增加假设对Hou–Luo情形的适用性。

Comment: 7 pages

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

1. 研究对象和背景

在不可压缩流体方程有限时间奇点形成的研究中，自相似方法扮演着核心角色。特别是针对三维欧拉方程（3D Euler equations）的奇点形成问题，Hou-Luo 模型及其数值构造提供了重要的理论框架。该研究旨在通过构建近似的自相似爆破轮廓（blowup profile），并证明其非线性稳定性来确立奇点的存在性。

研究对象具体为二维正半平面 $\mathbb{R}^2_+$ 上的 Boussinesq 方程的动力缩放形式，其控制方程组如下：

\begin{aligned} \omega_t + (c_l x + u) \cdot \nabla \omega &= \theta_x + c_\omega \omega, \\ \theta_t + (c_l x + u) \cdot \nabla \theta &= c_\theta \theta, \end{aligned}

其中速度场 $$u = (u, v)^T$$ 由 Biot-Savart 律确定：

-\Delta \phi = \omega, \quad u = -\phi_y, \quad v = \phi_x,

并满足边界条件 $\phi(x, 0) \equiv 0$ 。精确的自相似爆破解对应于上述系统的稳态方程：

(c_l x + u) \cdot \nabla \omega = \theta_x + c_\omega \omega, \quad (c_l x + u) \cdot \nabla \theta = c_\theta \theta.

本文试图解决的争议背景源于近期 Zhang [9] 的工作。Zhang 证明了轴对称欧拉方程在某些假设下不存在自相似剖面，并质疑这些非存在性结果是否适用于 Hou-Luo 场景。更重要的是，Zhang 指出在 Chen-Hou 构造的近似自相似剖面中存在一个“非局部误差”的潜在差异，认为该差异超出了先前文献中声称的精度范围（如 $10^{-7}$ ），从而暗示了数值构造的有效性存疑或无法支撑稳定性分析。

这一争议的实质在于对误差项定义的混淆。Chen 和 Hou 在此文中指出，Zhang 所关注的远场不匹配实际上是求解泊松方程产生的原始非局部 Poisson 误差（raw nonlocal Poisson error），而非用于非线性稳定性证明中的轮廓残差误差（profile residual error）。由于 Biot-Savart 算子的非局部性质，数值上求解流函数 $\phi$ 会引入误差 $\bar{\varepsilon} = \bar{\omega} - (-\Delta)\bar{\phi}_N$ 。然而，这个误差进入动力学方程时，并非直接以 $\bar{\varepsilon}$ 的形式出现，而是通过速度扰动 $u(\bar{\varepsilon})$ 乘以对流项的梯度因子进入残差。在远场区域，这些梯度因子具有显著的衰减特性，使得非局部误差对总残差的贡献被大幅抑制。因此，本文的研究空白在于厘清这两个不同量级的误差概念，纠正因误用原始 Poisson 误差而导致的对数值构造精度的错误解读，并验证基于图形证据的假设（如角向增长假设）是否成立。

2. 主要定理或主要结果

本文主要证明了以下两个关键结论，旨在澄清 Zhang [9] 提出的质疑：

非局部误差与轮廓残差误差的量级分离：作者严格区分了原始非局部误差 $\bar{\varepsilon}$ 与轮廓残差误差 $\bar{F}_i$ 。定义轮廓残差为稳态方程的偏差：
$\bar{F}_\omega = -(\bar{c}_l x + \bar{u}) \cdot \nabla \bar{\omega} + \bar{\theta}_x + \bar{c}_\omega \bar{\omega}, \quad \bar{F}_\theta = -(\bar{c}_l x + \bar{u}) \cdot \nabla \bar{\theta} + \bar{c}_\theta \bar{\theta}.$
将残差分解为局部部分和非局部部分：
$\bar{F}_1^{\text{loc}} = -(\bar{c}_l x + \bar{u}_N) \cdot \nabla \bar{\omega} + \bar{\theta}_x + \bar{c}_\omega \bar{\omega}, \quad \bar{F}_1^{\text{nloc}} = -u(\bar{\varepsilon}) \cdot \nabla \bar{\omega}.$
在远场区域（ $|x| \gtrsim L_2 \approx 10^{15}$ ），利用剖面的衰减估计 $|\nabla^i \bar{\omega}| \lesssim \langle x \rangle^{\bar{\alpha}_1 - i}$ 以及奇异积分算子估计 $|u(\bar{\varepsilon})| \lesssim |x| \langle x \rangle^{\bar{\alpha}_1}$ ，证明了非局部残差项的衰减速率显著快于原始误差：
$|\bar{F}_1^{\text{nloc}}| \leq C \langle x \rangle^{2\bar{\alpha}_1}, \quad |\bar{F}_{i+1}^{\text{nloc}}| \leq C \langle x \rangle^{3\bar{\alpha}_1}, \quad i=1,2.$
其中 $\bar{\alpha}_1 \in [-0.343, -0.342]$ 。这意味着即使原始非局部误差 $\bar{\varepsilon}$ 在网格点 $$(36, 720)$$ 处达到 $O(10^{-5})$ 的量级（如 Zhang 所指出的），经过衰减因子 $\langle x \rangle^{\bar{\alpha}_1} \approx 2.7 \times 10^{-5}$ 和 $\langle x \rangle^{2\bar{\alpha}_1} \approx 7.4 \times 10^{-10}$ 的加权后，其对总残差 $\bar{F}_i$ 的贡献已降至 $O(10^{-10})$ 甚至更低，完全满足稳定性分析所需的摄动小量要求。因此，Zhang 指出的 $10^{-5}$ 差异并不构成对残差估计有效性的反驳。
角向增长假设的数值证伪：针对 Zhang [9, Remark 2.4] 引用的 Figure 1 作为支持角向增长假设（Assumption (2.14)）的证据，作者利用原始数据文件 solu.w{1,1} 进行了严格的数值检验。通过将角度区间 $[0, \pi/2]$ 划分为 $$n=200$$ 个扇区，计算每个扇区内 $\bar{\omega}$ 的最大网格点值 $A_k = \max_{\beta_{ij} \in I_k} \bar{\omega}(x_i, x_j)$ 。结果显示，随着角度 $\beta$ 的增加， $$A_k$$ 呈现单调递减趋势（如图 2 所示），并未观察到任何径向最大值随角度增加的迹象。这直接推翻了 Zhang 基于可视化图形推断的“沿黄色脊线上升”的结论，表明相同的网格点数据不支持将角向增长假设应用于当前的 Hou-Luo 近似剖面。

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

本文的证明依赖于对数值构造中误差传播机制的精细分析，以及对数值数据的重新处理，其核心技术障碍和创新点如下：

1. 误差项的层级分解与衰减机制分析：核心障碍在于理解非局部算子如何影响全局残差。作者没有停留在宏观的能量估计层面，而是深入到具体的算子结构。首先，利用 Biot-Savart 律将速度场分解为 $u(\bar{\omega}) = \bar{u}_N + u(\bar{\varepsilon})$ ，其中 $\bar{u}_N$ 是基于数值流函数的近似速度， $u(\bar{\varepsilon})$ 是由泊松方程离散误差引起的非局部速度扰动。关键的创新在于识别出 $u(\bar{\varepsilon})$ 进入残差方程时的“二次衰减”效应。在远场区域，对流项 $u \cdot \nabla \omega$ 中的 $\nabla \omega$ 本身具有 $\langle x \rangle^{\bar{\alpha}_1}$ 的衰减速率。虽然奇异积分算子 $\nabla \nabla^\perp (-\Delta)^{-1}$ 是零阶算子，保持 $\bar{\varepsilon}$ 的衰减阶数不变，但速度 $u(\bar{\varepsilon})$ 在远场的表现受 $|x|^{-1}$ 因子的调节（见公式 2.10），且必须与梯度的衰减相乘。具体地，作者推导了不等式：

|\bar{F}_1^{\text{nloc}}| \leq |x|^{-1} |u(\bar{\varepsilon})| \cdot (|x| \cdot |\nabla \bar{\omega}|).

这里 $|x| \cdot |\nabla \bar{\omega}| \sim \langle x \rangle^{\bar{\alpha}_1}$ ，而 $|x|^{-1} |u(\bar{\varepsilon})| \sim \langle x \rangle^{\bar{\alpha}_1}$ 。两者相乘得到 $\langle x \rangle^{2\bar{\alpha}_1}$ 的衰减速率。这种数学结构是纯解析的，不依赖于特定的数值网格，它揭示了为什么原始误差 $\bar{\varepsilon}$ 的大小（如 $10^{-5}$ ）不能直接等同于残差的大小。作者还利用了加权范数（weighted norms） $au^2 \mu^4 |\bar{F}_1 \phi_{g1}|$ 等，确保在能量估计中这些高阶衰减项能够主导，从而保证摄动分析的封闭性。

2. 基于原始数据的扇区统计分析：针对图形证据的争议，作者摒弃了依赖低分辨率渲染图的主观视觉判断，转而采用严谨的离散数据分析方法。工具包括：

原始数据提取：直接从 MATLAB 文件 Steady_state_pertb720_Nlewcor4.mat 中提取未插值的网格点数据 solu.w{1,1}。
扇区极值统计：将角度空间离散化为 200 个扇区 $$I_k$$ ，计算每个扇区内的最大网格值 $$A_k$$ 。这种方法避免了插值带来的平滑效应，直接反映了数值解在网格上的真实分布特征。
单调性检验：通过绘制 $$A_k$$ 随 $$k$$ （即角度）变化的序列，直观且定量地展示了数据的递减趋势。这一方法不仅反驳了“角向增长”的假设，也指出了之前图形解释中可能存在的透视失真（perspective distortion）误导。

3. 数值解释的修正：作者还澄清了关于“单点消失”（one-point vanishing）的误解。利用刚性定理（rigidity property）对于精确解的结论不能直接套用到带有残差误差的数值近似解上。作者强调了计算机辅助证明（computer-assisted proof）的核心逻辑：数值解不需要在机器精度下逐点等于精确解，只需要在加权范数下的残差足够小。因此，孤立网格点的舍入误差或下溢（underflow）到零，只要被包含在残差控制的范围内，就不会导致数值解崩溃。

4. 与之前工作的比较

本文与 Zhang [9] 的工作构成了直接的对话与辨析，主要体现在以下几个维度的推进与差异：

1. 误差定义的精确性与适用范围的澄清： Zhang [9] 的主要推进在于提出了轴对称欧拉方程在非轴点附近的降维结果，并尝试将其应用于 Hou-Luo 场景。然而，其在比较数值构造时，错误地将求解泊松方程的原始误差 $\bar{\varepsilon}$ （量级约 $10^{-5}$ ）等同于用于稳定性证明的轮廓残差 $\bar{F}_i$ 。本文指出，这种混淆导致了错误的结论，即认为数值构造的精度未达到 $10^{-7}$ 的标准。本文通过解析推导证明，由于远场衰减因子的存在，实际的有效残差远小于原始误差，因此 Zhang 指出的“超出精度”并不构成对稳定性分析有效性的挑战。本文的推进在于明确了“原始误差”与“有效残差”之间的尺度关系，恢复了 Hou-Luo 数值构造在理论上的自洽性。

2. 假设适用性的实证检验： Zhang [9] 基于文献 [2] 中的可视化图形（Figure 1），主观推断 Hou-Luo 剖面满足角向增长假设（ $\sup_r W(r, \beta_2) > \sup_r W(r, \beta_1)$ ），以此作为应用其非存在性定理的前提。本文则通过重新处理原始网格数据，发现该假设在数值构造中并不成立（数据呈现角向递减）。这表明 Zhang 的论证链条在前提验证环节存在缺陷。本文的贡献在于提供了基于第一手数据的反例，否定了将该特定假设强加于当前数值剖面的合理性，同时也说明不能简单地将图形视觉特征转化为数学假设的依据。

3. 对计算机辅助证明逻辑的再阐释： Zhang [9] 似乎隐含地认为数值解必须在逐点意义上逼近精确解，否则即为无效。本文明确区分了“精确解的存在性”与“数值近似解的残差控制”。在计算机辅助证明框架下，只要加权残差范数可控，即可建立非线性稳定性，无需数值解与精确解在机器精度下逐点重合。本文澄清了这一点，指出 Zhang 关于“单点消失导致矛盾”的论点忽略了数值解仅满足弱意义方程（带残差项）的事实，从而错误地引入了不适用的刚性约束。

4. 代价与局限：本文的代价在于未能解决 Zhang 提出的其他维度缩减问题（dimension reduction results），仅专注于澄清误差定义和数值解释的具体技术细节。此外，本文虽然证伪了基于当前数值数据的角向增长假设，但并未排除在真实的物理精确解中可能存在某种形式的角向增长（尽管目前尚无证据支持）。本文的作用主要是“拨乱反正”，消除因概念混淆导致的误解，而非提出全新的理论突破。

12. 图上的莱本森方程

Leibenson’s equation on graphs

作者: Philipp Sürig, Xinrong Zhao

类别: math.AP

发布日期: 2026-08-15

arXiv ID: 2608.15168v1

摘要翻译: 本文在无限图上研究莱本森方程

\partial_t u = Δ_p u^q,

其中 $$p>1$$ , $$q>0$$ ，且 $$Δ_p$$ 表示离散 $$p$$ -拉普拉斯算子。我们证明，对于任意可积初值数据 $$u_0$$ ，存在全局解；当 $$p$$ 和 $$q$$ 处于某一特定范围时，该解是唯一的。假设费伯–克朗不等式（Faber–Krahn inequality）成立，我们获得了尖锐的 $\oldsymbol{\ell}^1$ - $\\boldsymbol{\ell}^\infty$ 正则化估计，以及关于初始支集有限解的传播量的定量界。在关于 $$p$$ 和 $$q$$ 的某些假设下，若图满足等周不等式（isoperimetric inequality），我们还证明了解在有限时间内消失的结果。特别地，在具有多项式体积增长的凯莱图上，当 $$q(p-1)>1$$ 时，我们确立了非负有限质量解的 $\\boldsymbol{\ell}^\infty$ -范数的最优长时间衰减速率，并展示了关于耗尽解有限时间消失的尖锐二分性结果。

Comment: 32 pages, no figures

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

1. 研究对象和背景

本文研究定义在无限加权图 $G=(V, E, \omega)$ 上的双非线性抛物型方程，即莱本森（Leibenson）方程：

\partial_t u = \Delta_p u^q, \quad x \in V, \ t > 0,

其中 $$p > 1$$ , $$q > 0$$ 。算子 $\Delta_p$ 为离散 $$p$$ -拉普拉斯算子，定义为：

\Delta_p f(x) = \frac{1}{\mu(x)} \sum_{y \sim x} |f(y) - f(x)|^{p-2}(f(y) - f(x)) \omega_{xy},

这里 $\mu(x) = \sum_{y \sim x} \omega_{xy}$ 是顶点 $$x$$ 的测度， $\omega_{xy}$ 为边权。该方程属于多孔介质、快扩散及双非线性扩散方程族在离散几何结构上的推广。当 $$q=1$$ 时，方程退化为抛物 $$p$$ -拉普拉斯方程；当 $$p=2$$ 且 $$q=1$$ 时，即为离散热方程。

连续情形下（如黎曼流形或欧氏空间），此类方程的正则性理论、有限传播速度、渐近行为及质量守恒性质已得到深入发展（参见 DiBenedetto, Vázquez 等经典文献）。然而，在图论背景下，由于缺乏微分结构，且图的拓扑与几何性质（如体积增长、等周不等式）对解的行为具有决定性影响，现有理论尚不完整。特别是，关于一般参数 $$p>1, q>0$$ 下的解的存在唯一性、质量守恒条件、 $\ell^1$ - $\ell^\infty$ 平滑估计以及有限时间熄灭（finite-time extinction）现象的系统性研究仍属空白。

本文旨在填补这一空白：首先建立任意可积初值下全局解的存在性，并确定解的唯一性范围；其次，引入 Faber-Krahn 不等式作为几何假设，推导 sharp 的 $\ell^1$ - $\ell^\infty$ 正则化估计及有界支撑解的传播速度界限；最后，利用等周不等式分析解的有限时间熄灭行为，并在多项式体积增长的 Cayley 图上刻画大时间衰减率及熄灭/非熄灭的二元性。这项工作将连续 PDE 中的双非线性扩散理论成功拓展至离散图结构，揭示了图几何参数与方程参数之间的精细相互作用。

2. 主要定理或主要结果

本文建立了以下核心结果：

存在性与唯一性（定理 1.3, 1.4）：对于任意 $u_0 \in \ell^r(V)$ ( $$r>1$$ )，初值问题 (1.2) 存在 exhaustion solution $u \in L^\infty([0, +\infty); \ell^r(V))$ 。若 $q(p-1) \ge 1$ 且 $u_0 \in \ell^1(V)$ ，则解满足 $u \in L^\infty([0, T); \ell^1(V))$ 。进一步，若 $q(p-1) \ge r \ge 1$ ，则比较原理成立：若 $u_0 \le v_0$ ，则 $u(t) \le v(t)$ 对所有 $$t$$ 成立，从而保证了解的唯一性。
质量守恒（定理 1.5）：设图 $$G$$ 满足体积增长条件 $\mu_V(B_R) \le C(1+R)^d$ 。若参数满足 $q(p-1) \ge \frac{d-p}{d}$ ，则对任意 $u_0 \in \ell^1(V)$ 的 exhaustion solution，总质量守恒：
$\sum_{x \in V} u(x,t)\mu(x) = \sum_{x \in V} u_0(x)\mu(x), \quad \forall t > 0.$
反之，若 $q(p-1) < \frac{d-p}{d}$ ，质量可能不守恒。
$\ell^1$ - $\ell^\infty$ 平滑估计（定理 1.7）：假设图满足 Faber-Krahn 不等式 (1.4) 且 $$q(p-1) > 1$$ 。对非负初值 $u_0 \in \ell^r(V)$ ，解满足 sharp 估计：
$\|u(t)\|_{\ell^\infty(V)} \le C \|u_0\|_{\ell^r(V)} \phi_r\left( t^{-1} \|u_0\|_{\ell^r(V)}^{-q(p-1)-1} \right)^{1/r},$
其中 $\phi_r$ 是由 Faber-Krahn 函数 $\Lambda_p$ 导出的严格递增双射函数。
传播速度界限（定理 1.8）：对有有限支撑的非负初值 $u_0 \in \ell^1(V)$ ，支撑集随时间扩张。对任意 $\varepsilon \in (0,1)$ ，当半径 $$R$$ 满足：
$R \ge c \varepsilon^{-1} t^{1/p} \|u_0\|_{\ell^1}^{\frac{q(p-1)-1}{p}} \phi_1\left( t^{-1} \|u_0\|_{\ell^1}^{-[q(p-1)-1]} \right)^{\frac{q(p-1)-1}{p}}$
时，有 $\|u(t)\|_{\ell^1(B_R)} \ge (1-\varepsilon)\|u_0\|_{\ell^1}$ 。
有限时间熄灭（定理 1.9）：若图满足 $$d$$ -等周不等式 ( $$d>p>1$$ ) 且 $q(p-1) < \frac{d-p}{d}$ ，取 $m = \frac{d(1-q(p-1))}{p} > 1$ ，则对 $u_0 \in \ell^m(V)$ ，解在有限时间 $T \le C \|u_0\|_{\ell^m}^{mp/d}$ 内熄灭（即 $u(\cdot, t) \equiv 0$ ）。特别地，在多项式体积增长为 $$N$$ 的 Cayley 图上，当 $0 < q(p-1) < \frac{N-p}{N}$ 时，所有 exhaustion solution 均有限时间熄灭；而当 $q(p-1) \ge \frac{N-p}{N}$ 时，质量守恒，永不熄灭。此即严格的二元性。
大时间衰减率（推论 1.10）：在多项式增长 $$N$$ 的 Cayley 图上，若 $$q(p-1) > 1$$ ，则对非负有限质量解，其 $\ell^\infty$ 范数满足最优衰减率：
$\|u(t)\|_{\ell^\infty} \asymp t^{-\alpha^*}, \quad \alpha^* = \frac{N}{N[q(p-1)-1] + p}.$
上界由定理 1.7 给出，下界由定理 1.8 结合质量守恒导出。

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

证明面临的核心障碍在于离散结构中缺乏光滑性，且双非线性项 $$u^q$$ 与 $$p$$ -拉普拉斯算子的耦合导致能量估计复杂化。作者通过以下关键技术推进了证明：

截断逼近与对角论证：针对无限图上的存在性问题，作者未直接构造解，而是采用 exhaustion 方法。定义在有限球 $$B_n(x_0)$$ 上的混合问题 (1.3)，利用 Peano 定理获得局部解。关键步骤是通过能量估计证明解序列 $\{u_n\}$ 在 $\ell^r$ 范数下一致有界且等度连续。具体地，利用恒等式：
$\sum_x (\Delta_p u^q) v \mu = -\frac{1}{2} \sum_{x,y} |\nabla_{xy} u^q|^{p-2} \nabla_{xy} u^q \nabla_{xy} v \omega_{xy},$
选取测试函数 $v = u^{r-1}_+$ ，导出 $\frac{d}{dt}\|u\|_{\ell^r}^r \le 0$ ，从而获得先验估计。随后通过对角线法提取收敛子列，并利用图的局部有限性验证极限函数满足原方程。此框架避免了直接处理无限维系统的紧性问题。
广义最大值原理与比较原理：为证明唯一性及非负性，作者构建了适用于双非线性情形的比较原理。难点在于 $$u^q$$ 项可能导致非线性退化或奇异性。作者引入截断函数 $\zeta_{R_1, R_2}$ ，乘以 $(u-h)_+^{r-1}$ 后求和积分。通过精细估计交叉项 $$I_2$$ （涉及 $\nabla \zeta$ 的误差项），证明当 $R \to \infty$ 时误差消失。关键在于利用条件 $q(p-1) \ge r$ 确保非线性项的可积性，使得 H"older 不等式能有效控制边界项。由此导出：若 $u_0 \le v_0$ ，则 $u(t) \le v(t)$ 。此结论不仅保证了唯一性，还隐含了非负性保持（令 $h \to 0$ ）及强极大值原理（若某点为零则全域为零）。
Faber-Krahn 不等式与变分技巧：为获得 $\ell^1$ - $\ell^\infty$ 平滑估计，作者并未使用传统的 Moser 迭代，而是利用了图上的 Faber-Krahn 不等式 (1.4)。该不等式提供了第一特征值与集合体积的关系：
$\Lambda_p(\mu_V(U)) \|f\|_{\ell^p}^p \le \mathcal{E}_p(f),$
其中 $\mathcal{E}_p$ 为 Dirichlet 能量。通过构造特定的测试函数序列并结合 $\psi_r(s)$ 函数的单调性，作者将能量衰减转化为 $L^\infty$ 范数的衰减。这一方法巧妙地将几何约束（Faber-Krahn 函数 $\Lambda_p$ ）融入解析估计，得到了依赖图几何参数的 sharp 常数 $$C$$ 。
等周不等式与熄灭机制分析：对于有限时间熄灭，作者利用 $$d$$ -等周不等式 $\mu_V(\Omega)^{\frac{d-1}{d}} \le C_d \mu_E(\partial \Omega)$ 来控制质量流失。通过选取适当的测试函数并分析能量不等式中的边界项，证明了当 $q(p-1) < \frac{d-p}{d}$ 时，能量衰减速度快于质量扩散，导致解在有限时间内归零。这一分析揭示了参数 $$p, q$$ 与图维度 $$d$$ 之间的临界阈值关系，填补了离散情形下熄灭条件的理论空白。
传播速度的定量控制：在有有限支撑初值的情况下，作者通过分析 $\ell^1$ 范数在膨胀球内的分布，结合 Faber-Krahn 估计，导出了支撑集半径 $$R(t)$$ 的增长速率公式。该公式显式依赖于 $t^{1/p}$ 及初始质量的幂次，精确刻画了信息在图上的传播速度，区别于连续情形中的波前概念，体现了离散扩散的“跳跃”特性。

4. 与之前工作的比较

本文显著推进了既有 PDE 与分析理论在图结构上的应用，具体对比如下：

与 Andreucci & Tedeev [AT20] 的比较： [AT20] 研究了 $$q=1$$ 时的抛物 $$p$$ -拉普拉斯方程在满足 Faber-Krahn 不等式图上的大时间行为，给出了 $\ell^1$ - $\ell^\infty$ 估计。本文将其推广至一般 $$q>0$$ 的双非线性情形。主要差异在于：
- 参数范围：[AT20] 仅覆盖 $$q=1$$ ，而本文处理 $$q(p-1)>1$$ 及 $$q(p-1)<1$$ 的全范围，特别是揭示了 $$q(p-1)$$ 在决定质量守恒与熄灭性质中的临界作用。
- 估计形式：本文得到的衰减率指数 $\alpha^* = \frac{N}{N[q(p-1)-1]+p}$ 在 $$q=1$$ 时退化为 [AT20] 的结果，但本文进一步证明了该指数的最优性（上下界匹配），而 [AT20] 仅给出上界。
- 技术难度：双非线性项 $$u^q$$ 使得能量估计中的插值不等式更为复杂，本文通过引入 $\psi_r$ 函数解决了这一困难，这是 [AT20] 无需处理的。
与 Hua & Mugnolo [HM15] 的比较： [HM15] 研究了无穷图上抛物 $$p$$ -拉普拉斯方程（ $$q=1$$ ）的有限时间熄灭问题，指出当 $$p$$ 较小时在等周条件下会发生熄灭。本文的贡献在于：
- 一般化熄灭条件：[HM15] 的熄灭条件依赖于 $$p$$ 的大小，而本文给出了明确的参数区域 $q(p-1) < \frac{d-p}{d}$ ，统一了 $$p$$ 与 $$q$$ 的影响。例如，即使 $$p$$ 较大，只要 $$q$$ 足够小，仍可能发生熄灭，这扩展了 [HM15] 的结论。
- 二分性证明：本文在 Cayley 图上严格证明了熄灭与非熄灭的 sharp dichotomy（当 $q(p-1) = \frac{N-p}{N}$ 时为临界状态），而 [HM15] 未涉及这种临界相变现象。
- 质量守恒的精确刻画：本文定理 1.5 给出了质量守恒的充要条件（基于体积增长 $$d$$ 和参数 $$q(p-1)$$ ），而 [HM15] 仅定性讨论了质量是否守恒。
与连续理论（DiBenedetto, Vázquez 等）的比较：连续理论中，莱本森方程的有限传播速度（finite propagation speed）是已知现象（当 $$q(p-1)>1$$ 时）。然而，本文发现一个反直觉的结论：在图结构上，即使 $$q(p-1)>1$$ ，解也不具备经典的有限传播速度（即支撑集不会以有限速度扩张，而是瞬间扩散到整个连通分支，尽管 $\ell^1$ 范数集中在某个区域）。这体现在 Remark 2.4 中：若某点为零，则全域为零。这与连续情形截然不同，揭示了离散几何对扩散性质的根本性修正。
局限性与代价：本文的进展是以更强的几何假设为代价的：必须假设图满足 Faber-Krahn 不等式或等周不等式。对于一般图（如无界曲率或异常体积增长），这些估计可能失效。此外，唯一性仅在 $q(p-1) \ge r$ 范围内得到证明，对于 $$q(p-1) < 1$$ 的情况，唯一性尚未完全解决（需进一步研究）。尽管如此，本文建立的框架为后续研究更复杂的图结构上的双非线性方程奠定了坚实基础。

摘要翻译: 我们研究奇异摄动完全非线性方程

F(D^2u_\ uarepsilon) = \frac\alpha{\varepsilon} \beta\left(\frac{u_\varepsilon}{\varepsilon}\right) \qquad\text{in }B_1\subset\mathbb R^n,

的变号粘性解，其中 $$F$$ 是一致椭圆的，且 $\beta\in C_c (-1,1)$ 是非负的。我们证明了尺度尖端的估计

中文标题: 非局部Tikhonov正则化：Hilbert尺度、显式速率与经典极限

摘要翻译: 我们研究有界Lipschitz域上线性反问题的分数阶Sobolev Tikhonov正则化。正则化惩罚项由受限Dirichlet分数阶Laplacian生成，并证明相关的变分问题存在唯一极小值点，且该极小值点对数据满足Lipschitz连续依赖关系。通过识别正自伴算子

A_s=I+(-Δ)^s, ext{ } D(A_s^{1/2})=H_0^s(Ω)，

我们将问题等距变换为具有观测算子 $B=KA_s^{-1/2}$ 的经典Hilbert空间Tikhonov问题。这导出了显式的均方误差界以及Holder型源条件下的最优先验和事后参数规则。该框架通过部分观测和反向分数阶热方程进行说明。在后一种情形下，

中文标题: 波长尺度的可选停时、临界 Feynman-Kac 规范与容量谱不等式

摘要翻译: 我们发展了两种布朗运动停时方法：一种利用击中概率和容量，另一种利用加权出口分布的矩。在闭 $$m$$ -维黎曼流形上（ $$m eq 2$$ ），第一种方法给出了支撑于零体积集合上的平衡测度的低能谱不等式。结合尺度给出了热可观测性和从具有全支撑的测度出发的零控制性，该测度由 Hausdorff 维数为 $$m-2$$ 的稠密集承载，且在受控 $$H^1$$ -迹条件下达到最小。对于每个满足 $$Φ(4E) ext{ extless} DΦ(E)$$ 的非递减无界函数 $$Φ$$ ，可以选择测度使得最优谱常数与 $$Φ(E)$$ 相当，且小时间控制代价为 $T^{-1/2}Φ(T^{-1})^{1/2}$ 。平方 Feynman-Kac 鞅使势加倍：对于径向势 $$V$$ ，若 $$r^2V(r)$$ 非递减，则在正则解 $( frac12Δ+2V)g=0$ 保持正定的区域，边界二阶矩关于 $$ r$$ 的对数 $$ r$$ 是对数凸的；当 $$V(0)>0$$ 时，因子 $$2$$ 是尖锐的。在小测地球上，加权二阶矩给出几乎单调频率和加权三半径不等式；其对数导数的界限给出倍增估计，经基态变换后，支持 Logunov-Malinnikova 的 Remez 定理。归一化一阶矩出口律给出拉普拉斯特征函数的平均边界方差恒等式。在 $r=cλ^{-1/2}$ 处，除承载 $O(c^{2/m})$ 的 $$ r^2$$ -质量集合外， $$ r(x) r>0 r rcap B(x,r)$$ 中通过 $$x$$ 的分量包含一个半径接近完整的同心球，填满几乎整个 $$B(x,r)$$ ，并几乎匹配其第一 Dirichlet 特征值。一阶矩停时也控制正超水平分量。反射布朗运动和局部时给出有界 Lipschitz 域上 Steklov 特征函数的边界到领域的 $L^ ext{ extless}∞ r$ 估计，以及 $$C^2$$ 域上的 $$L^2$$ 估计。

不相关论文（基本信息）

1. 关于梯度映射的离散性与开性的 Pogorelov 型反例

A Pogorelov-type counterexample to the discreteness and openness of gradient mappings

作者: Deguang Zhong

类别: math.AP, math.CV

发布日期: 2026-08-17

arXiv ID: 2608.16757v1

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

摘要: Let $Ω\subset\mathbb{R}^{n}$ be a domain, and suppose that $u\in W^{2,n}_{loc}(Ω)$ satisfies the following ineqiality

\det D^2u\geq δ>0\qquad\text{a.e. in }Ω.

A question of Guerra–Tione \cite[Question 5.5]{GuerraTione} asks whether the gradient mapping $$Du$$ must be open and discrete. In this paper, we give an explicit Pogorelov-type construction showing that the answer is negative in every dimension $n\geq4$ : there exists $u\in W^{2,n}_{loc}(Ω)$ satisfying the above lower bound for the Hessian determinant, with $$D^2u>0$$ a.e., such that $$Du$$ collapses an entire line segment to a single point and hence is not discrete. We also show that the same construction has a logarithmic divergence when $$n=3$$ and therefore does not directly settle the three-dimensional case.

2. 常数正 $$σ_k$$ -曲率度量的光滑奇点集：模型阈值维数界与尖锐临界端

A Model-threshold Dimension Bound and Sharp Critical Ends for Smooth Singular Sets of Constant Positive $$σ_k$$ -curvature Metrics

作者: Jiahuan Li, Yilu Liu, Xi-Nan Ma

类别: math.DG, math.AP

发布日期: 2026-08-17

arXiv ID: 2608.16571v1

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

摘要: Let $k\in\mathbb N$ satisfy $$1<k<n/2$$ , and let $Σ^p\subset\Sn^n$ be a closed smooth embedded submanifold. We prove that a complete conformal metric $g=v^{-2}g_{\Sn^n}$ on $\Sn^n\setminusΣ$ satisfying $λ(g^{-1}A_g)\in\Gk$ and $σ_k(g^{-1}A_g)=κ>0$ must obey

p\leq p_k(n),

where $$p_k$$ is the model threshold determined by $\Hh^{p+1}\times\Sn^{n-p-1}$ . When $$k=2$$ and $$n=m^2$$ , we construct a smooth complete equality example on $\Sn^n\setminus\Sn^{(m^2-m-2)/2}$ . We also prove that the strict inequality $$p<p_k(n)$$ holds under a finite positive linear-contact hypothesis.

3. 空间形式中测地三角形第一非零 Neumann 特征值的一个尖锐直径依赖下界

A Sharp Diameter-Dependent Lower Bound for the First Nonzero Neumann Eigenvalue of Geodesic Triangles in Space Forms

作者: Shoo Seto, Guofang Wei, Yusen Xia

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

发布日期: 2026-08-17

arXiv ID: 2608.16136v1

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

摘要: We prove a sharp lower bound for the first nonzero Neumann eigenvalue of geodesic triangles of given diameter in two-dimensional space forms. The bound is given by the first positive radial Neumann eigenvalue of an one dimensional model; it is approached by degenerating isosceles triangles. When $$K>0$$ and $D=π/(2\sqrt K)$ , equality is attained precisely by birectangular triangles. We also prove a hot-spots theorem for non-acute spherical triangles of diameter at most $$π/2$$ , and establish antisymmetry and eigenvalue monotonicity for isosceles spherical triangles of diameter $$π/2$$ .

4. Log-Concave Compactly Supported Targets的Brenier映射的无维Lipschitz界

Dimension-Free Lipschitz Bounds for Brenier Maps to Compactly Supported Log-Concave Targets

作者: Maja Gwozdz

类别: math.AP, math.PR

发布日期: 2026-08-16

arXiv ID: 2608.15906v1

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

摘要: We fix an integer $d\ge1$ and a symmetric positive-definite matrix $Q\in\mathbb{R}^{d\times d}$ . Let $V:\mathbb{R}^d\to\mathbb{R}$ be finite, set

Z_μ:=\int_{\mathbb{R}^d}e^{-V(x)}\,dx\in(0,\infty), \qquad dμ(x):=Z_μ^{-1}e^{-V(x)}\,dx,

and assume that $$μ$$ has finite second moment and that

x\longmapsto \frac12\langle Qx,x\rangle-V(x)

is convex. Let $$ν$$ be a compactly supported log-concave probability measure with support $$K$$ , and let $\nablaΦ$ be the Brenier map from $$μ$$ to $$ν$$ . For $v\in\mathbb{R}^d$ , define

w_K(v):= \sup_{y\in K}\langle y,v\rangle - \inf_{y\in K}\langle y,v\rangle.

We prove that

\partial_{vv}Φ\le 0.587 \sqrt{\langle Qv,v\rangle}\,w_K(v) \qquad(v\in\mathbb{R}^d)

in the sense of distributions. We show that $\nablaΦ$ has an everywhere-defined globally Lipschitz representative such that

\operatorname{Lip}(\nablaΦ) \le 0.587 \sqrt{\|Q\|_{\mathrm{op}}}\,\operatorname{diam}(K).

The directional Hessian estimate is affinely covariant, whereas the global Lipschitz estimate is dimension-free. The result also applies to singular or lower-dimensional targets. In particular, it removes the $\sqrt d$ loss in Kolesnikov’s estimate for the Brenier map from Gaussian measure to normalised Lebesgue measure on a convex body. We also prove new bounds that depend only on the support for compactly supported semi-log-concave targets, which includes targets with bounded negative curvature.

5. 各向异性 Calderón 问题中的拟解析性与几何刚性

Quasianalyticity and geometric rigidity in anisotropic Calderón’s problem

作者: Lu Chen, Yan Jiang, Hongyu Liu, Longyue Tao

类别: math.AP

发布日期: 2026-08-16

arXiv ID: 2608.15813v1

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

摘要: The anisotropic Calderón problem in dimensions $n\ge3$ remains open for general smooth metrics~\cite{Uhlmann2009}. We establish uniqueness results in two complementary regimes. In the first, the identity principle for quasianalytic functions propagates boundary information and yields uniqueness in general geometry, including a partial-boundary consequence; under a prescribed normal geometry, quasianalyticity is needed only in the distinguished direction. In the second, suitable symmetry or one-sided ordering assumptions lead to uniqueness at $C^\infty$ regularity with full or restricted boundary access. Taken together, the results exhibit a tradeoff among regularity, geometric structure, and boundary access: quasianalyticity supplies continuation in general geometry, while symmetry or one-sided order replaces that continuation at $C^\infty$ regularity.

6. $\mathbb{R}^{n}$ 上某些多值调和函数的刻画

Characterisation of some multivalued harmonic functions on $\mathbb{R}^{n}$

作者: Yang Li, Parsa Mashayekhi, Yichen Zhang

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

发布日期: 2026-08-16

arXiv ID: 2608.15722v1

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

摘要: Using different approaches, Simon Donaldson and Dashen Yan recently constructed $\mathbb Z_2$ -harmonic functions on $\mathbb R^3$ and on $\mathbb R^n$ for $n\geqslant 3$ , respectively. Their examples have quadratic growth at infinity and $\mathcal{O}(r^{3/2})$ local growth near a smooth codimension-two branching set $$Σ$$ . We show that these properties uniquely characterise such $\mathbb Z_2$ -harmonic functions up to rigid motions of $\mathbb R^n$ .

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

1. 二维定常低耗散 Navier–Stokes 系统的 Liouville 定理#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 三阶对称尺度不变欧拉流的梯度增长与跳跃剖面的弛豫#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 具有锐角接触角的可压缩毛细 - 重力水波局部适定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 具有条件稳定性的抛物型反源问题的无伴随积分反馈方法#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 具有孤立边界奇点的临界半线性椭圆方程 II#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 二维Hénon问题中的新变号解#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 可压缩欧拉方程剪切层的长波谱不稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. 偏微分方程中驻定周期图案尖点分岔的严格验证#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

9. 具有不规则漂移的超布朗运动的紧支集性质#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

10. 扩散方程的随机团块方法#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

11. 关于计算机辅助证明3D欧拉奇异性中非局部误差与轮廓残差误差区别的澄清#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

12. 图上的莱本森方程#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 接触型广义哈密顿梯度流 I：基本解的正则性#

摘要翻译#

2. 抽象希尔伯特空间中变双线性形式特征值的稳定性及应用#

摘要翻译#

3. 强竞争下对称Lotka-Volterra系统中波速符号的完全刻画#

摘要翻译#

4. 半空间和正象限上尖锐的 $L^2$-Caffarelli–Kohn–Nirenberg 与加权 Poincaré 不等式及其稳定性#

摘要翻译#

5. 几何流下热方程正解的光化估计#

摘要翻译#

6. 具有简并扩散和紧支集初始种群密度的趋化系统中避免卡尔德拉型死核#

摘要翻译#

7. 有界密度下排斥库仑流的定量平均场极限与Riesz弱强稳定性#

摘要翻译#

8. 两类相奇异摄动完全非线性椭圆方程的一致Lipschitz正则性#

摘要翻译#

9. 具有 $L^{4/3+}$ 活性标量的 SQG 方程的灵活性#

摘要翻译#

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

1. 二维定常低耗散 Navier–Stokes 系统的 Liouville 定理

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 三阶对称尺度不变欧拉流的梯度增长与跳跃剖面的弛豫

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 具有锐角接触角的可压缩毛细 - 重力水波局部适定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 具有条件稳定性的抛物型反源问题的无伴随积分反馈方法

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 具有孤立边界奇点的临界半线性椭圆方程 II

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 二维Hénon问题中的新变号解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 可压缩欧拉方程剪切层的长波谱不稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 偏微分方程中驻定周期图案尖点分岔的严格验证

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

9. 具有不规则漂移的超布朗运动的紧支集性质

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

10. 扩散方程的随机团块方法

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

11. 关于计算机辅助证明3D欧拉奇异性中非局部误差与轮廓残差误差区别的澄清

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

12. 图上的莱本森方程

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 接触型广义哈密顿梯度流 I：基本解的正则性

摘要翻译

2. 抽象希尔伯特空间中变双线性形式特征值的稳定性及应用

摘要翻译

3. 强竞争下对称Lotka-Volterra系统中波速符号的完全刻画

摘要翻译

4. 半空间和正象限上尖锐的 $$L^2$$ -Caffarelli–Kohn–Nirenberg 与加权 Poincaré 不等式及其稳定性

摘要翻译

5. 几何流下热方程正解的光化估计

摘要翻译

6. 具有简并扩散和紧支集初始种群密度的趋化系统中避免卡尔德拉型死核

摘要翻译

7. 有界密度下排斥库仑流的定量平均场极限与Riesz弱强稳定性

摘要翻译

8. 两类相奇异摄动完全非线性椭圆方程的一致Lipschitz正则性

摘要翻译

9. 具有 $L^{4/3+}$ 活性标量的 SQG 方程的灵活性

摘要翻译