Home » 论文总结

2026年09月15日论文分析

共有 11 篇重点关注论文, 31 篇论文可以了解, 70 篇不相关论文

September 15, 2026 · 98 min · 20608 words · P0lyno3ial

Table of Contents

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

报告日期: 2026-09-15（补报）

检索区间（UTC）: 2026-09-11T18:00:00+00:00 至 2026-09-14T18:00:00+00:00（不含结束时刻）

实际补生成时间: 2026年09月18日 22:26:59

运行状态: 完成

计划处理论文数量: 112

已完成论文数量: 112

未完成论文数量: 0

重点关注论文数量: 11

了解领域论文数量: 31

不相关论文数量: 70

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

1. 一维可压缩 Navier–Stokes 方程在一般压强定律下强解的局部与整体存在性

Local and Global Existence of Strong Solutions to the One-Dimensional Compressible Navier–Stokes Equations with General Pressure Law

作者: Qingsong Zhao

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15317v1

摘要翻译: 我们研究一维完整可压缩 Navier–Stokes 方程的 Cauchy 问题，其中采用一般的本构关系 $$p=p(v,θ)$$ 和 $$e=e(v,θ)$$ 。假设压强和内能足够光滑，并且在热力学相容意义下满足 $$e_v=θp_θ-p$$ ，且满足 $$e_θ>0$$ 。首先，我们在不对压强施加单调性条件的情况下，建立了大初值数据下强解的局部时间存在性与唯一性。对于围绕常平衡态 $$(ar v,0,arθ)$$ 的扰动，我们进一步假设机械稳定性条件 $$p_v(ar v,arθ)<0$$ 。通过引入 Gibbs 关系中的相对热力学能量，我们推导了一个基本能量恒等式，并证明了其在平衡态附近的局部二次强制性。将此估计与高阶先验界相结合，我们获得了充分小的 $H^1(\mathbb R)$ 扰动下强解的整体存在性与唯一性。此外，比容和温度在所有时刻均保持远离零和无穷大的均匀有界性。

Comment: 无

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

1. 研究对象和背景

本文研究一维非等熵可压缩 Navier-Stokes 方程组在 Lagrangian 质量坐标下的 Cauchy 问题. 该方程组描述了粘性且导热的流体运动, 其形式为:

\begin{cases} v_t - u_x = 0, \\ u_t + p_x = \left(\mu \frac{u_x}{v}\right)_x, \\ \left(e + \frac{1}{2}u^2\right)_t + (up)_x = \left(\kappa \frac{\theta_x}{v}\right)_x + \left(\mu \frac{u u_x}{v}\right)_x. \end{cases}

其中 $$v(t,x) > 0$$ 为比容, $$u(t,x)$$ 为速度, $\theta(t,x) > 0$ 为绝对温度, $$e(t,x)$$ 为比内能, $$p(t,x)$$ 为压强. 初始条件为 $(v,u,\theta)(0,x) = (v_0(x), u_0(x), \theta_0(x))$ , 并满足远场条件 $\lim_{x \to \pm \infty}(v_0, u_0, \theta_0) = (\bar v, 0, \bar \theta)$ , 其中 $\bar v > 0, \bar \theta > 0$ 为常数. 黏度系数 $\mu > 0$ 和热导率 $\kappa > 0$ 假设为正常数.

传统的数学理论主要关注理想多温气体模型, 其状态方程为 $p(v,\theta) = R\theta/v$ 和 $e(v,\theta) = C_v \theta$ . 然而, 对于真实流体, 特别是在高压、高温或复杂分子相互作用下, 压强和内能对体积和温度的依赖关系可能更为复杂. 因此, 本文考虑更一般的本构函数 $p, e \in C^2(\mathbb{R}_+ \times \mathbb{R}_+)$ , 并假设 $e_\theta(v,\theta) > 0$ . 此外, 压强和内能需满足热力学相容性关系:

e_v(v,\theta) = \theta p_\theta(v,\theta) - p(v,\theta).

这一关系等价于存在比熵 $s(v,\theta)$ 满足 Gibbs 关系 $\theta ds = de + p dv$ . 利用上述关系, 总能量方程可以改写为温度方程:

c(v,\theta)\theta_t + \theta p_\theta(v,\theta)u_x = \left(\kappa \frac{\theta_x}{v}\right)_x + \mu \frac{u_x^2}{v},

其中 $c(v,\theta) := e_\theta(v,\theta) > 0$ 为定压比热. $e_\theta > 0$ 保证了温度方程时间导数的非退化性.

已有文献中, Kazhikhov 和 Shelukhin 以及 Matsumura 和 Nishida 等人建立了一维可压缩 Navier-Stokes 方程的全局存在性理论, 但大多依赖于理想气体状态方程或特定的非理想模型. 尽管近年来在非理想气体、真空边界条件和温度依赖性输运系数方面取得了一些进展, 但对于完全一般的热力学相容函数 $p(v,\theta)$ 和 $e(v,\theta)$ 的全局理论仍相对薄弱. 本文旨在填补这一空白, 即在不施加压强单调性条件的情况下, 建立局部存在性; 并在机械稳定性条件下, 证明小扰动下的全局存在性.

2. 主要定理或主要结果

本文的主要结果包括局部存在性定理和全局存在性定理.

定理 1.1 (局部存在性): 假设本构函数满足 $$C^2$$ 正则性和热力学相容性 (1.7), 且 $e_\theta > 0$ . 若初始数据 $(v_0 - \bar v, u_0, \theta_0 - \bar \theta) \in H^1(\mathbb{R})$ 且满足点态界 $c_0 \le v_0(x), \theta_0(x) \le C_0$ , 则存在时间 $$T > 0$$ 使得 Cauchy 问题存在唯一的强解, 满足:

(v - \bar v, u, \theta - \bar \theta) \in C(0, T; H^1(\mathbb{R})),

(u_x, \theta_x) \in L^2(0, T; H^1(\mathbb{R})),

v_t \in L^2(0, T; H^1(\mathbb{R})), \quad (u_t, \theta_t) \in L^2(0, T; L^2(\mathbb{R})).

此外, 比容和温度保持均匀有界: $c_0/2 \le v(t,x), \theta(t,x) \le 2C_0$ . 注意, 此局部结果不需要对 $$p_v$$ 或 $p_\theta$ 施加任何符号条件.

为了推导全局估计, 作者引入了相对于平衡态 $(\bar v, 0, \bar \theta)$ 的相对热力学势:

H(v,\theta | \bar v, \bar \theta) := e(v,\theta) - e(\bar v, \bar \theta) + \bar p(v - \bar v) - \bar \theta (s(v,\theta) - s(\bar v, \bar \theta)).

命题 1.3 (相对能量恒等式与局部强制性): 若机械稳定性条件 $p_v(\bar v, \bar \theta) < 0$ 成立, 则 $$H$$ 在 $(\bar v, \bar \theta)$ 附近局部非负且具有二次强制性. 具体地, 存在 $\delta > 0$ 和常数 $$0 < c_* < C_*$$ 使得当 $|v - \bar v| + |\theta - \bar \theta| \le \delta$ 时:

c_* (|v - \bar v|^2 + |\theta - \bar \theta|^2) \le H(v,\theta | \bar v, \bar \theta) \le C_* (|v - \bar v|^2 + |\theta - \bar \theta|^2).

相应的能量恒等式为:

\int_\mathbb{R} \left[ \frac{1}{2}u^2 + H(v,\theta | \bar v, \bar \theta) \right] dx + \bar \theta \int_0^t \int_\mathbb{R} \left( \mu \frac{u_x^2}{v\theta} + \kappa \frac{\theta_x^2}{v\theta^2} \right) dx d\tau = \int_\mathbb{R} \left[ \frac{1}{2}u_0^2 + H(v_0,\theta_0 | \bar v, \bar \theta) \right] dx.

定理 1.4 (小扰动下的全局存在性): 在上述假设基础上, 若进一步满足机械稳定性条件 $p_v(\bar v, \bar \theta) < 0$ , 则存在足够小的 $\varepsilon_0 > 0$ 使得当初始扰动满足 $\|(v_0 - \bar v, u_0, \theta_0 - \bar \theta)\|_{H^1(\mathbb{R})} \le \varepsilon_0$ 时, Cauchy 问题存在唯一的全局强解, 满足:

(v - \bar v, u, \theta - \bar \theta) \in C(0, \infty; H^1(\mathbb{R})),

(u_x, \theta_x) \in L^2(0, \infty; H^1(\mathbb{R})).

此外, 比容和温度在整个时间域上保持远离零和无穷大的均匀界:

0 < \underline v \le v(t,x) \le \bar v < \infty, \quad 0 < \underline \theta \le \theta(t,x) \le \bar \theta < \infty.

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

证明的核心技术障碍在于如何在没有压强单调性假设的情况下控制比容和温度, 从而防止解在有限时间内出现真空或奇异性. 作者的真正创新点在于引入基于 Gibbs 关系的相对热力学能量结构, 并利用其在平衡态附近的局部二次强制性来构建全局先验估计.

首先, 在局部存在性的证明中, 作者采用了 Picard 迭代法. 由于初始数据仅具有 $$H^1$$ 正则性, 直接取时间无关的零阶迭代会导致抛物项正则性不足. 因此, 作者使用热半群对 $$u$$ 和 $\theta$ 进行正则化, 构造了合适的近似序列. 通过归纳法, 证明了迭代序列在适当的范数下的一致有界性. 关键在于处理变系数抛物方程中的非线性项, 利用一维 Gagliardo-Nirenberg 不等式和 Young 不等式来控制交叉项, 并结合 Gronwall 不等式得到统一估计. 在收缩性证明中, 作者选择低阶度量进行收缩, 避免了在变系数差值中损失一个空间导数, 从而在 $$C^2$$ 本构假设下自然地恢复了强解的正则性.

其次, 在全局存在性的证明中, 核心工具是相对能量方法. 作者定义了相对热力学势 $H(v,\theta | \bar v, \bar \theta)$ , 并通过 Gibbs 关系将其与熵联系起来. 通过计算, 作者推导出了精确的能量耗散恒等式 (1.17). 该恒等式表明, 相对总能量的变化由耗散项决定, 而耗散项包含 $u_x^2/(v\theta)$ 和 $\theta_x^2/(v\theta^2)$ 的贡献. 为了确保该能量估计能够控制解的行为, 必须证明 $$H$$ 在平衡态附近的强制性. 作者利用机械稳定性条件 $p_v(\bar v, \bar \theta) < 0$ 和 $e_\theta > 0$ , 证明了 $$H$$ 的 Hessian 矩阵在 $(\bar v, \bar \theta)$ 处正定, 从而得到了局部二次强制性估计 (1.19). 这一步骤至关重要, 因为它将抽象的热力学势转化为具体的 $$L^2$$ 范数控制, 使得能量估计能够闭合.

最后, 结合相对能量估计和高阶先验估计, 作者使用了 Bootstrap-Continuation 论证. 通过证明解在任意有限时间区间上的范数不会爆破, 并将局部存在性延拓到全局, 最终证明了全局强解的存在性. 同时, 利用比较原理和 Sobolev 嵌入, 作者证明了比容和温度的均匀上下界, 从而排除了真空和奇异性的出现.

4. 与之前工作的比较

与既有文献相比, 本文的主要推进体现在以下几个方面:

本构函数的普遍性: 大多数现有的全局存在性结果, 如 Kawohl [12] 和 Jiang [8, 9] 的工作, 虽然放宽了理想气体假设, 但仍要求对本构函数施加特定的结构性假设, 如单调性、增长性或强制性条件, 以控制比容和温度. 相比之下, 本文在局部存在性部分完全不施加 $$p_v < 0$$ 或 $p_\theta > 0$ 等单调性条件, 仅依赖 $e_\theta > 0$ 保证温度方程的非退化性. 这代表了向“完全一般”热力学相容函数迈出的重要一步.
能量结构的创新: 传统方法往往依赖于理想气体的特殊代数结构 (如 $p = R\theta/v$ ) 来简化能量估计. 本文通过引入基于 Gibbs 关系的相对热力学能量 $$H$$ , 成功地将一般本构函数的热力学结构融入能量估计中. 这种相对能量方法不仅提供了正确的耗散结构, 而且通过局部强制性分析, 使得在小扰动下能够闭合高阶估计. 这与以往直接对原始变量进行能量估计的方法有所不同, 后者在处理非理想气体时往往面临困难.
假设的精细化: 在局部存在性中, 作者明确指出不需要压强的全局机械或热力学稳定性. 而在全球存在性中, 仅要求平衡态处的机械稳定性 $p_v(\bar v, \bar \theta) < 0$ . 这种区分使得理论适用范围更广, 允许压强在其他区域表现出不稳定行为, 只要局部平衡态是稳定的即可. 这与一些要求全局单调性的先前工作形成了对比.
与近期工作的联系: 文中提到了 Li and Xin [14], Cai et al. [1, 2], Liao et al. [16] 等近期工作, 这些工作分别处理了真空、退化热导率、电离气体等特殊情形. 本文则聚焦于本构自由度本身, 不依赖于特定的物理模型 (如 Saha 电离定律) 或特殊的边界条件. 此外, 与 Zhao [23] 研究的超双曲弛豫模型不同, 本文处理的是经典的粘性导热系统, 其奇异性机制和数学结构有本质区别.
局限性与代价: 本文的全局存在性结果仅限于小 $$H^1$$ 扰动. 对于大初值, 尽管局部存在性成立, 但无法保证全局存在性, 因为相对能量的局部强制性仅在平衡态附近有效. 此外, 结果依赖于 $$C^2$$ 正则性的本构函数, 对于更低正则性或奇异本构函数 (如相变过程中的非光滑行为), 该方法可能不再适用. 尽管如此, 本文为一般可压缩流体的数学理论提供了一个坚实的框架, 并为未来研究更大扰动或更复杂本构关系奠定了基础.

2. Navier-Stokes-Maxwell 系统整体解的最优时间衰减

Optimal Time-Decay of global solutions to the Navier-Stokes-Maxwell system

作者: Haroune Houamed, Slim Ibrahim, Belkacem Said-Houari

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15095v1

摘要翻译: 我们证明了 Navier–Stokes–Maxwell (NSM) 系统在临界 Fujita–Kato 空间 $\dot H^{\frac{d}{2}-1}(\mathbb R^d)$ 中小初值下整体时间解的存在性和唯一性，适用于任意维度 $d\geq 3$ 。在额外假设初值属于 Besov 空间 $\dot B^{-\frac{d}{2}}_{2,\infty}(\mathbb R^d)$ 的条件下，我们还建立了解在 $\dot H^s(\mathbb R^d)$ 中无穷远处的最优衰减率 $t^{-\frac{s}{2}-\frac{d}{4}}$ ，对任意 $s\in\left(-\frac{d}{2},\frac{d}{2}-1\right]$ 成立。这一结果通过构造一个与 Fourier 侧逐点能量等价的 Lyapunov 泛函，并结合 Fourier 分裂法的适应性应用来实现其最优衰减率。

Comment: 26 pages

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

1. 研究对象和背景

本文研究对象是描述粘性、导电流体与电磁场耦合演化的不可压 Navier-Stokes-Maxwell (NSM) 系统. 该系统在三维空间中的经典形式为:

\begin{cases} \partial_t u + u \cdot \nabla u = \nu \Delta u - \nabla p + j \times B, \\ \text{div } u = 0, \\ \frac{1}{c}\partial_t E - \nabla \times B = -j, \\ \frac{1}{c}\partial_t B + \nabla \times E = 0, \\ \text{div } B = 0, \end{cases}

其中 $$u$$ 为流体速度, $$E$$ 和 $$B$$ 分别为电场和磁场, $$j$$ 为电流密度. 参数 $\nu$ , $\sigma$ 和 $$c$$ 分别代表流体的粘度、电导率和光速. 电流密度 $$j$$ 通过 Ohm 定律与电磁场及流体速度耦合, 具体形式为 $j = \sigma(cE + u \times B)$ 或其无散版本 $j = \sigma(cE + P(u \times B))$ , 其中 $$P$$ 为 Leray 投影算子.

NSM 系统的分析难度源于其混合的抛物-双曲结构. 一方面, 它继承了 Navier-Stokes 方程的非线性对流项带来的正则性困难; 另一方面, Maxwell 方程的双曲性质导致电磁场缺乏直接的耗散机制, 使得能量不等式无法提供足够的先验控制. 在三维情况下, 仅凭能量不等式不足以建立适定性理论. 此外, Lorentz 力 $j \times B$ 的非线性耦合进一步增加了分析的复杂性.

为了克服这些困难, 作者将 NSM 系统推广到任意维度 $d \geq 3$ . 由于传统的向量叉积和旋度算子在高于三维时不再具有相同的代数性质, 作者利用微分形式的语言重新表述了该系统. 在此框架下, 速度和磁场被表示为 1-形式, 而电场和电流密度被表示为 2-形式. 三维的叉积 $u \times B$ 被外积 $u \wedge B$ 替代, 旋度 $\nabla \times B$ 被外微分 $$dB$$ 替代, 而 Lorentz 力 $j \times B$ 则对应于内积 $\iota_B j$ . 这种推广保持了系统在三维时的经典形式, 同时允许在更高维度中进行一致的分析.

已有文献中, Masmoudi [28] 在二维情况下建立了小初值的全局存在性和唯一性. 在三维情况下, 早期的工作通常要求速度场具有额外的 Besov 正则性, 例如 $\dot B^{d/2-1}_{2,1}$ . 随后, [18] 将这一条件放宽至临界 Sobolev 空间 $\dot H^{d/2-1}$ , 但仅限于三维. 关于长时间行为的研究更为有限, 目前仅在轴对称假设下有一些结果 [4]. 本文旨在建立一个统一的框架, 不仅在任意维度 $d \geq 3$ 下证明小初值的全局存在性和唯一性, 而且给出解的最优时间衰减率. 特别地, 所有估计均对光速 $$c$$ 保持一致性, 从而在非相对论极限 $c \to \infty$ 下可以恢复磁流体力学 (MHD) 系统的相应结果.

2. 主要定理或主要结果

本文的主要结果包括两个核心定理: 全局存在性与唯一性定理 (Theorem 3.1) 以及最优时间衰减定理 (Theorem 3.6).

Theorem 3.1 (全局存在性与唯一性): 设 $d \geq 3$ . 存在一个足够小的常数 $\delta > 0$ , 使得对于任意满足

\|(u_0, E_0, B_0)\|_{\dot H^{d/2-1}(\mathbb R^d)} < \delta

的初值 $(u_0, E_0, B_0) \in (\dot H^{d/2-1}(\mathbb R^d))^3$ , NSM 系统存在唯一的全局时间解 $$(u, E, B)$$ , 并且该解满足如下估计:

\|(u, E, B)\|_{C(\mathbb R^+; \dot H^{d/2-1})} + \|cE\|_{L^2(\mathbb R^+; \dot H^{d/2-1})} + \|u\|_{L^2(\mathbb R^+; \dot B^{d/2}_{2,1})} \leq 2 \|(u_0, E_0, B_0)\|_{\dot H^{d/2-1}}.

值得注意的是, 速度场 $$u$$ 在 $L^2(\mathbb R^+; \dot B^{d/2}_{2,1})$ 中的精细控制是本工作的关键改进之一, 这比先前工作中得到的 $L^2(\mathbb R^+; \dot H^{d/2} \cap L^\infty)$ 控制更强.

Theorem 3.6 (最优时间衰减): 在 Theorem 3.1 的基础上, 若初值额外属于 Besov 空间 $\dot B^{-d/2}_{2,\infty}(\mathbb R^d)$ , 即

(u_0, E_0, B_0) \in (\dot H^{d/2-1}(\mathbb R^d))^3 \cap (\dot B^{-d/2}_{2,\infty}(\mathbb R^d))^3,

则全局解 $$(u, E, B)$$ 在 $t \to \infty$ 时具有以下最优衰减率:

\|(u(t), E(t), B(t))\|_{\dot H^s(\mathbb R^d)} \lesssim t^{-\frac{s}{2} - \frac{d}{4}}, \quad \forall s \in \left(-\frac{d}{2}, \frac{d}{2}-1\right].

这里的衰减指数 $-\frac{s}{2} - \frac{d}{4}$ 与热方程解的衰减率一致, 表明尽管 NSM 系统包含双曲部分, 但在大时间尺度上, 其耗散行为主导了动力学. 此外, 所有结果对光速 $c \gg 1$ 保持一致, 因此在非相对论极限 $c \to \infty$ 下, 可以直接推导出对应的 MHD 系统的全局存在性和最优衰减结果.

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

本文的证明策略分为两个主要部分: 全局存在性的建立和大时间行为的分析.

1. 全局存在性与精细估计: 核心障碍在于如何在临界正则性水平下闭合非线性估计, 特别是如何处理 Lorentz 力 $j \times B$ 带来的非线性项. 作者的关键创新是利用了一个由 Arsénio 和第一作者在 [6] 中建立的尖锐抛物估计:

\|e^{t\Delta} u_0\|_{L^2(\mathbb R^+; \dot B^{d/2}_{2,1}(\mathbb R^d))} \lesssim \|u_0\|_{\dot H^{d/2-1}(\mathbb R^d)}.

这一估计允许作者将速度场 $$u$$ 控制在 $L^2(\mathbb R^+; \dot B^{d/2}_{2,1})$ 空间中, 而不是更弱的 $L^2(\mathbb R^+; \dot H^{d/2} \cap L^\infty)$ 空间. 这种更强的控制对于处理非线性项至关重要, 因为它提供了更好的可积性, 从而使得在临界 Sobolev 空间中可以闭合能量估计. 具体来说, 通过将速度场的精细估计与电磁场的能量估计相结合, 作者能够证明全局解的存在性和唯一性.

2. 最优时间衰减率的推导: 大时间行为的分析面临的主要挑战是如何从局部存在性过渡到全局衰减估计, 尤其是在缺乏直接耗散机制的情况下. 作者采用了两步法:

Lyapunov 泛函的构造: 首先, 作者构造了一个 Lyapunov 泛函, 该泛函在 Fourier 侧等价于点态能量. 这个泛函的设计巧妙地捕捉了系统中耗散和非耗散部分的相互作用, 使得可以通过单调性论证来推导衰减估计.
Fourier 分裂法的适应: 其次, 作者将经典的 Fourier 分裂方法适应于 NSM 系统. 通过将频率空间分解为低频和高频部分, 作者分别处理不同频率区间的贡献. 在低频区间, 利用 Lyapunov 泛函的单调性和抛物型半群的平滑效应; 在高频区间, 利用初始数据的 Besov 正则性 $\dot B^{-d/2}_{2,\infty}$ 来保证快速衰减. 结合这两部分, 作者得到了整体的最优衰减率 $t^{-\frac{s}{2} - \frac{d}{4}}$ .

3. 一致性与非相对论极限: 另一个重要的技术点是确保所有估计对光速 $$c$$ 保持一致. 这意味着在推导过程中, 任何依赖于 $$c$$ 的项都必须以可控的方式出现, 使得当 $c \to \infty$ 时, 估计不会退化. 通过仔细追踪每个估计中对 $$c$$ 的依赖关系, 作者证明了在非相对论极限下, NSM 系统收敛于 MHD 系统, 并且相应的全局存在性和衰减结果也成立. 这一一致性不仅增强了结果的普适性, 还为理解 NSM 系统与 MHD 系统之间的关系提供了严格的数学基础.

4. 与之前工作的比较

与现有文献相比, 本文在以下几个方面做出了显著推进:

1. 维度的推广: 先前的工作主要集中在二维或三维情况. 例如, Masmoudi [28] 的结果限于二维, 而 [18] 虽然放宽了正则性假设, 但仍限于三维. 本文通过将 NSM 系统推广到任意维度 $d \geq 3$ , 并利用微分形式的语言进行统一处理, 打破了维度的限制. 这种推广不仅丰富了理论框架, 还为研究高维等离子体模型提供了新的工具.

2. 正则性假设的优化: 在三维情况下, 早期工作如 [22] 要求速度场具有额外的 Besov 正则性 $\dot B^{d/2-1}_{2,1}$ . 虽然 [18] 将其放宽至 $\dot H^{d/2-1}$ , 但本文在此基础上进一步改进了速度场的精细控制, 从 $L^2(\mathbb R^+; \dot H^{d/2} \cap L^\infty)$ 提升到 $L^2(\mathbb R^+; \dot B^{d/2}_{2,1})$ . 这一改进不仅有助于闭合全局估计, 还为后续的大时间行为分析奠定了基础.

3. 大时间行为的完整刻画: 关于 NSM 系统大时间行为的研究此前非常有限, 仅在轴对称假设下有零星结果 [4]. 本文首次在无轴对称假设下, 在任意维度 $d \geq 3$ 中建立了全局解的最优时间衰减率. 这一结果填补了文献中的重要空白, 并为理解 NSM 系统的长期动力学提供了定量描述.

4. 一致性与非相对论极限: 许多先前工作在推导估计时可能隐含了对光速 $$c$$ 的依赖, 或者未明确讨论 $c \to \infty$ 的极限. 本文强调所有估计对 $$c$$ 的一致性, 并明确展示了如何从 NSM 系统过渡到 MHD 系统. 这种一致性不仅提高了结果的鲁棒性, 还澄清了 NSM 系统与 MHD 系统之间的数学联系, 为非相对论极限下的渐近分析提供了严格依据.

局限性与未来方向: 尽管本文取得了重要进展, 但仍存在一些局限性. 例如, 结果仅适用于小初值, 大初值情况下的全局存在性和大时间行为仍是开放问题. 此外, 虽然得到了最优衰减率, 但关于渐近剖面的精确刻画 (例如, 是否存在特定的自相似解或波包结构) 尚未深入探讨. 未来的工作可以考虑扩展到大初值情形, 或者结合谱分析方法进一步细化渐近行为.

3. 超粘性可压缩三维 Oldroyd 型模型的尖锐非唯一性：超越 Lions 指数

Sharp non-uniqueness for the hyper-viscous compressible 3D Oldroyd-type models: Beyond the Lions exponent

作者: Haobin Li, Peng Qu, Zirong Zeng, Mingxin Zhang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15073v1

摘要翻译: 我们证明了具有一般压力律和牛顿粘性的三维超粘性可压缩 Oldroyd 型模型弱解的非唯一性。利用一种间歇凸积分方案，我们在 $L_t^γW_x^{s,p}$ 中构造了非唯一的弱解，其中 $$(s,p,γ)$$ 位于相对于 Ladyzhenskaya-Prodi-Serrin 判据的超临界区域。此外，我们所考虑的牛顿粘性包含分数阶粘性 $$(-Δ)^α$$ ，其中粘性指数 $$α$$ 可以大于 Lions 指数 $$5/4$$ 。在经典粘性情形 $$α=1$$ 下，我们的构造还给出了 $$L_t^pL_x^ \infty$$ 中 $$1<p<2$$ 时的尖锐非唯一性。

Comment: 无

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

1. 研究对象和背景

本文研究三维可压缩 Oldroyd 型流体模型在弱解意义下的非唯一性问题. 该模型源自稀薄聚合物溶液的动力学理论, 是描述非牛顿弹性流体行为的经典方程组. 具体而言, 作者考虑定义在环面 $\mathbb{T}^3 := [-\pi, \pi]^3$ 上的如下耦合系统:

\begin{cases} \partial_t \rho + \text{div}(\rho u) = 0, \\ \partial_t (\rho u) + \text{div}(\rho u \otimes u) + \nabla P(\rho) = \text{div} S(u) + \text{div}(\rho F F^T), \\ \partial_t F + u \cdot \nabla F - (\nabla u) F = 0, \\ \text{div}(\rho F^T) = 0, \end{cases}

其中 $\rho(t,x)$ 为质量密度, $$u(t,x)$$ 为速度场, $$F(t,x)$$ 为变形梯度张量. $$S(u)$$ 表示牛顿粘性应力, 定义为:

S(u) = 2\nu_s \left( D(u) - \frac{1}{3}\text{div}u I \right) + \nu_b \text{div}u I,

其中 $D(u) = \frac{1}{2}(\nabla u + \nabla^T u)$ 为变形张量, $\nu_s > 0$ 和 $\nu_b + \frac{1}{3}\nu_s \geq 0$ 分别为剪切粘度和体积粘度系数. 压力函数 $P(\rho)$ 假设关于 $\rho$ 具有 $$C^2$$ 正则性, 且满足单调性条件 $P(\rho) > 0$ , $P'(\rho) > 0$ 以及 $P''(\rho) > 0$ . 约束条件 $\text{div}(\rho F^T) = 0$ 被称为 Piola 条件, 它在定义弱解及保持物理一致性中起着关键作用.

为了探究更广泛的耗散效应, 作者将 Newtonian 粘性项中的拉普拉斯算子替换为分数阶拉普拉斯算子 $(-\Delta)^\alpha$ , 其中 $\alpha \in [1, \frac{3}{2})$ . 通过引入动量 $m = \rho u$ 和加权变形梯度 $N = \rho F$ , 方程组可改写为守恒律形式:

\begin{cases} \partial_t \rho + \text{div} m = 0, \\ \partial_t m + \nu_s (-\Delta)^\alpha (\rho^{-1} m) - (\nu_b + \frac{1}{3}\nu_s) \nabla \text{div}(\rho^{-1} m) + \nabla P(\rho) + \text{div}\left( \rho^{-1} (m \otimes m - N N^T) \right) = 0, \\ \partial_t N^i + \text{div}\left( \rho^{-1} (N^i \otimes m - m \otimes N^i) \right) = 0, \quad i=1,2,3, \\ \text{div} N^i = 0, \quad i=1,2,3. \end{cases}

当粘性系数为零时, 该系统退化为可压缩基本弹性动力学系统. 此外, 在低马赫数极限下 ( $Ma \to 0$ ), 该系统收敛于不可压 Oldroyd 型模型.

尽管关于 Oldroyd 型模型强解的存在唯一性已有大量研究 (如 Guillopé-Saut, Molinet-Talhouk, Chemin-Masmoudi 等人的工作), 但关于其弱解的非唯一性 (即灵活侧) 的研究相对匮乏. 特别是对于可压缩情形, 由于压力项 $\nabla P(\rho)$ 依赖于沿输运方程传播的密度, 导致压力选择具有“相对刚性”, 这使得应用凸积分方法构造非唯一弱解变得极具挑战性. 此前, 仅在低耗散 (hypo-dissipative, $\alpha < 1$ ) 的可压缩 Navier-Stokes 或 MHD 方程中取得了部分非唯一性结果. 对于超黏性 (hyper-viscous, $\alpha \geq 1$ ) 的可压缩流体, 尤其是 Oldroyd 型模型, 尚无已知的弱解非唯一性结果. 本文旨在填补这一空白, 证明即使在粘性指数 $\alpha$ 超过 Lions 指数 $$5/4$$ 的情况下, 弱解依然不唯一.

2. 主要定理或主要结果

本文的主要结果概括为以下四个定理:

定理 1.2 (主要构造定理): 设 $\alpha \in [1, \frac{3}{2})$ . 给定任意无散度张量场 $\bar{N}$ 以及输运方程 $\partial_t \bar{\rho} + \text{div} \bar{m} = 0$ 的光滑解 $(\bar{\rho}, \bar{m})$ , 满足 $0 < c_* \leq \bar{\rho}(t,x) \leq C_*$ , 则存在 $\beta' \in (0, 1)$ , 使得对于任意 $\varepsilon_* > 0$ 以及满足以下条件的指数 $(s, p, \gamma)$ :

(s, p, \gamma) \in A := \left\{ (s, p, \gamma) \in [0, \frac{3}{2}) \times [1, +\infty) \times [1, +\infty) : 0 \leq s < \frac{2\alpha}{\gamma + 2\alpha - \frac{2}{p}} + 1 - 2\alpha \right\},

存在超黏性可压缩 Oldroyd 型模型 (1.3) 的弱解 $(\rho, m, N)$ , 满足: (i) $(\rho, m, N)$ 是分布意义下的弱解; (ii) 正则性: $\rho \in C_t C_x^1$ , $m, N \in H_{t,x}^{\beta'} \cap L_t^\gamma W_x^{s,p} \cap C_t H_x^{-3/2}$ ; (iii) 质量守恒: $\int_{\mathbb{T}^3} \rho(t,x) dx = \int_{\mathbb{T}^3} \bar{\rho}(t,x) dx$ ; (iv) 范数小偏差: $\|\rho - \bar{\rho}\|_{C_t C_x^1} \leq \varepsilon_*$ , $\|m - \bar{m}\|_{L_t^1 L_x^2} \leq \varepsilon_*$ , $\|N - \bar{N}\|_{L_t^1 L_x^2} \leq \varepsilon_*$ ; (v) 时间支撑的小偏差: 扰动主要集中在光滑背景解的时间支撑附近.

定理 1.3 (非唯一性推论): 存在光滑初始数据 $(\rho_0, m_0, N_0)$ ( $\rho_0 > 0$ ), 使得对于任意满足上述区域 $$A$$ 的指数 $(s, p, \gamma)$ , 存在无穷多个属于 $C_t C_x^1 \times (H_{t,x}^{\beta'} \cap L_t^\gamma W_x^{s,p})^2$ 的弱解具有相同的初始数据. 特别地, 在经典粘性情形 $\alpha=1$ 下, 该结果蕴含了 $L_t^p L_x^\infty$ ( $$1 < p < 2$$ ) 空间中的尖锐非唯一性, 这对应于 Ladyzhenskaya-Prodi-Serrin 判据的超临界区域.

定理 1.4 (强零粘性极限): 设 $\alpha \in [1, \frac{3}{2})$ , $(\rho, m, N) \in C_{t,x}^{\tilde{\beta}} \times (H_{t,x}^{\tilde{\beta}})^2$ ( $\tilde{\beta} > 0$ ) 是可压缩基本弹性动力学系统 (无粘) 的任意弱解, 且 $c_1 \leq \rho \leq c_2$ . 则存在 $\beta' \in (0, \tilde{\beta})$ 和一系列超黏性可压缩 Oldroyd 型模型的弱解 $(\rho^{(n)}, m^{(n)}, N^{(n)})$ , 其粘性系数为 $\kappa_n \nu_s$ 和 $\kappa_n \nu_b$ , 使得当 $\kappa_n \to 0$ 时:

\rho^{(n)} \to \rho \quad \text{strongly in } C_{t,x}, \\ (m^{(n)}, N^{(n)}) \to (m, N) \quad \text{strongly in } (H_{t,x}^{\beta'})^2.

定理 1.5 (不可压极限/低马赫数极限): 设 $\alpha \in [1, \frac{3}{2})$ , $(v, H) \in (H_{t,x}^{\tilde{\beta}})^2$ ( $\tilde{\beta} > 0$ ) 是不可压 Oldroyd 型模型的任意弱解. 则存在 $\beta' \in (0, \tilde{\beta})$ 和一系列超黏性可压缩 Oldroyd 型模型 (无量纲形式, 马赫数为 $\kappa_n Ma$ ) 的弱解 $(\rho^{(n)}, m^{(n)}, N^{(n)})$ , 使得当 $\kappa_n \to 0$ 时:

\rho^{(n)} \to 1 \quad \text{strongly in } C_{t,x}, \\ (m^{(n)}, N^{(n)}) \to (v, H) \quad \text{strongly in } (H_{t,x}^{\beta'})^2.

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

本文的核心技术挑战在于克服可压缩流体中压力项的“相对刚性”以及控制超耗散项 $(-\Delta)^\alpha$ 带来的额外导数损失. 作者采用间歇凸积分 (intermittent convex integration) 方案, 并引入了新的抵消机制 (cancellation mechanism).

新的几何引理与抵消机制: 在传统的可压缩 Navier-Stokes 凸积分中, 非线性项 $w_{q+1} \otimes w_{q+1}$ 需要抵消 Reynolds 应力 $$R_q$$ 以及由压力引起的各向同性部分 $\varrho Id$ . 然而, 在 Oldroyd 型模型中, 弹性应力项 $$N N^T$$ 提供了额外的自由度. 作者证明了新的几何引理 (Lemma 3.2), 允许构造扰动 $w_{q+1}$ (动量) 和 $f_{q+1}$ (弹性场), 使得:
$w_{q+1} \otimes w_{q+1} - \sum_{i=1}^3 f_{q+1}^i \otimes f_{q+1}^i \simeq R_q,$
而不是简单的 $w_{q+1} \otimes w_{q+1} \simeq R_q + \varrho Id$ . 这种机制利用了弹性应力的贡献来抵消动量雷诺应力, 从而在不破坏压力刚性结构的前提下实现了压力的有效控制. 这是处理可压缩 Oldroyd 模型的关键创新.
Mikado 流与双重整歇性: 为了控制 $\alpha \in [1, \frac{3}{2})$ 的超耗散项, 作者使用了具有时间和空间双重间歇性的 Mikado 流作为构建块. 频率参数 $\lambda_q$ 和振幅参数 $\delta_q$ 的选择遵循特定的尺度关系:
$\lambda_q := a(b^q), \quad \delta_q := \lambda_q^{-2\beta} \lambda_1^{3\beta},$
其中 $$a, b$$ 为足够大的常数, $\beta$ 为正则性参数. 这种尺度选择确保了高阶导数项在迭代过程中可控, 并且保证了最终解落在 $L_t^\gamma W_x^{s,p}$ 空间中.
松弛系统与模糊化 (Mollification): 为了避免凸积分过程中的导数损失, 作者对松弛系统进行了时空模糊化处理. 模糊尺度 $l := \lambda_q^{-30}$ 远小于特征尺度 $\lambda_q^{-1}$ . 模糊化后产生的交换子应力 (commutator stress) $R_{com}$ 和 $M_{com}$ 被证明是高阶小量, 可以被后续的迭代步骤吸收.
归纳估计框架: 证明基于严格的归纳假设 (2.5)-(2.9). 这些估计控制了密度、动量、弹性场、Reynolds 应力和弹性应力波动在不同范数下的行为. 通过验证从第 $$q$$ 步到第 $$q+1$$ 步的误差估计 (如 $\|\rho_{q+1} - \rho_q\|_{C_t C_x^1} \lesssim \delta_{q+2}^{1/2}$ ), 作者证明了序列的收敛性以及最终解的正则性.
奇异极限的证明策略: 对于零粘性和不可压极限, 作者利用凸积分构造的灵活性, 直接构造出一系列具有小粘性系数或小马赫数的弱解, 并证明其在适当拓扑下的强收敛性. 这避免了传统分析中对强解正则性的依赖, 从而在低正则性框架下建立了这些极限定理.

4. 与之前工作的比较

本文的工作在以下几个方面显著推进了既有文献:

首次确立超黏性可压缩流体的非唯一性: 此前的凸积分非唯一性结果主要集中在不可压流体 (如 Euler, Navier-Stokes) 或低耗散 (hypo-dissipative, $\alpha < 1$ ) 的可压缩流体 (如 Li-Qu-Zeng-Zhang [38] 针对可压缩 Navier-Stokes, Li-Qu [36] 针对可压缩 MHD). 对于超黏性 ( $\alpha \geq 1$ ) 的可压缩流体, 由于耗散更强, 通常认为更倾向于良定性. 本文首次证明了即使 $\alpha$ 超过 Lions 指数 $$5/4$$ (即 $\alpha \in [1, 1.5)$ ), 可压缩 Oldroyd 型模型的弱解依然不唯一. 这表明超黏性并不能完全消除非唯一性机制, 且 Oldroyd 模型的弹性结构削弱了粘性的刚性约束.
超越 Lions 指数的适用范围: 对于不可压超耗散 Navier-Stokes 方程, Li-Qu-Zeng-Zhang [37] 已经证明了在 $L_t^\gamma W_x^{s,p}$ 超临界空间中的非唯一性. 本文将此结果推广到了可压缩 Oldroyd 模型. 关键在于, 可压缩系统中的压力项 $\nabla P(\rho)$ 比不可压系统中的压力 $\nabla p$ 更难处理, 因为它是密度的函数而非独立的 Lagrange 乘子. 作者通过引入弹性应力抵消机制，成功克服了这一障碍，将非唯一性范围扩展到了包含经典粘性 ( $\alpha=1$ ) 甚至更高粘性的区域.
经典粘性下的尖锐非唯一性: 在 $\alpha=1$ 的经典粘性情形下，本文证明了 $L_t^p L_x^\infty$ ( $$1 < p < 2$$ ) 中的非唯一性. 这与 Ladyzhenskaya-Prodi-Serrin 判据的端点 $(s,p,\gamma)=(0,\infty,2)$ 密切相关. 相比之下，之前的可压缩 Navier-Stokes 非唯一性结果往往局限于低耗散区域或更弱的函数空间. 本文的结果表明，Oldroyd 模型在经典粘性下也表现出显著的 ill-posedness, 这为理解超弹性流体中粘性失效的物理机制提供了数学依据.
低正则性下的奇异极限: 定理 1.4 和 1.5 分别在 $H_{t,x}^{\tilde{\beta}}$ 低正则性框架下建立了零粘性极限和低马赫数极限. 以往的奇异极限研究通常依赖于强解的存在性或较高的正则性假设. 本文利用凸积分构造的弱解族，直接在弱解层面证明了这些极限，拓展了奇异极限理论适用的函数空间范围. 特别是定理 1.5 是超黏性可压缩 Oldroyd 模型的首个低马赫数极限结果.
与 Chen-Liu [12] 的比较: Chen-Liu [12] 之前仅证明了三维低耗散弹性动力学 (hypo-viscous elastodynamics) 存在有限动能的非平凡弱解. 本文不仅将耗散指数提升到了超耗散区域 ( $\alpha \geq 1$ )，而且明确证明了非唯一性 (infinitely many solutions with same initial data)，而不仅仅是存在性. 此外，本文处理的模型包含了完整的 Oldroyd 型弹性应力演化，而非简化的弹性动力学.

4. 具有流入-流出边界条件的可压缩完整Navier-Stokes方程的爆破判据

Blow-up criterion to the compressible full Navier-Stokes equations with inflow-outflow boundary conditions

作者: C. H. Arthur Cheng, Young-Sam Kwon, Cheng-Fang Su

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14727v1

摘要翻译: 我们研究了带有热传导的可压缩Navier-Stokes方程在流入/流出边界条件下的强解存在性。本文的主要贡献在于为三维可压缩流建立了一个新的正则性判据，该判据作为不可压缩流体中Beale-Kato-Majda判据的对应物。值得注意的是，我们的结果改进了Valli和Zajaczkowski所建立的强解结论，以及Basaric、Feireisl和Mizerova提出的条件正则性框架。我们的分析扩展了Basaric、Feireisl和Mizerova的框架，使我们能够成功推导出所需的均匀估计。

Comment: 47 pages

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

1. 研究对象和背景

本文研究对象为三维可压缩、粘性且导热流体在带有真实进出流 (inflow-outflow) 边界条件下的 Navier-Stokes-Fourier 方程组。该方程组描述了密度 $\varrho$ 、速度 $$u$$ 和绝对温度 $\vartheta$ 的演化，其数学形式如下：

\begin{cases} \varrho_t + \text{div}(\varrho u) = 0 & \text{in } \Omega \times (0, T), \\ (\varrho u)_t + \text{div}(\varrho u \otimes u) = -\nabla p + \Delta u + \varrho f & \text{in } \Omega \times (0, T), \\ \varrho(\vartheta_t + u \cdot \nabla \vartheta) = \kappa \Delta \vartheta + |\nabla u|^2 - p \text{div} u + \varrho g & \text{in } \Omega \times (0, T), \end{cases}

其中压力由理想气体状态方程给出 $p = \varrho \vartheta$ 。边界条件设定在光滑有界区域 $\Omega \subset \mathbb{R}^3$ 上，边界被分为进流部分 $\Gamma_{\text{in}}$ 和出流部分 $\Gamma_{\text{out}}$ 。在 $\Gamma_{\text{in}}$ 上，法向速度满足 $u_B \cdot N < 0$ ，并给定密度 $\varrho = \varrho_B$ ；在 $\Gamma_{\text{out}}$ 上， $u_B \cdot N \ge 0$ ，不指定密度。整个边界上均给定速度和温度的 Dirichlet 数据。

这一问题的研究背景在于处理开放系统 (open system) 中的强解理论。与不可穿透边界 ( $u \cdot N = 0$ ) 不同，进出流边界使得质量可以跨越边界，导致连续性方程在法向呈现双曲特征。这给高阶 Sobolev 空间中的先验估计带来了本质困难，特别是在对连续性方程进行微分后产生的边界通量项中，涉及密度的空间导数。以往关于此类系统的研究主要集中在弱解或耗散弱解的全局存在性上（如 Feireisl, Novotný 等人的工作），或者是在不可穿透边界下的条件正则性理论（如 Basarić, Feireisl, Mizerová 的工作）。对于带有真实进出流的强解，Valli 和 Zajaczkowski 曾建立过 Hilbert 空间框架下的局部存在性理论，但其要求初始数据满足由方程生成的额外高阶相容性条件。Meliani 则在 $$L^p-L^q$$ 框架下建立了包含磁流体力学耦合的强解理论，但其参数范围限制较严，且未直接给出纯 Navier-Stokes-Fourier 系统在标准 Sobolev 空间下的最优正则性。

本文试图填补的空白是：在仅假设基本初值-边界相容性的情况下，构造高阶 Sobolev 空间的强解，并推导适用于真实进出流边界的 Beale-Kato-Majda 型爆破判据。由于进出流边界破坏了传统能量估计中利用 $u \cdot N = 0$ 进行的边界消去技巧，传统的基于 $L^\infty_t L^\infty_x$ 范数的延拓论证在此失效。因此，需要开发新的分析工具来处理微分后的连续性方程中的边界项，特别是利用进流边界处法向速度的符号来控制这些项，从而确立一个新的正则性准则。

2. 主要定理或主要结果

本文的主要结果包括两个核心定理：强解的局部存在性和基于 Beale-Kato-Majda 思想的爆破判据。

定理 1.4 (强解的局部存在性): 假设 $\partial \Omega$ 为 $$C^3$$ 类，外力 $$f$$ 和热源 $$g$$ 满足适当的 $$L^2$$ 正则性，边界数据 $u_B, \vartheta_B$ 具有相应的时间正则性。设初始数据 $(\varrho_0, u_0, \vartheta_0) \in H^2(\Omega)^5$ ，且满足 $\varrho_0 > 0, \vartheta_0 > 0$ 以及基本的迹相容条件 $u_0 = u_B, \vartheta_0 = \vartheta_B$ on $\partial \Omega$ 和 $\varrho_0 = \varrho_B(0)$ on $\Gamma_{\text{in}}$ 。则存在 $$T^* > 0$$ 及强解 $(\varrho, u, \vartheta)$ 满足：

$\varrho \in C([0, T^*]; H^2(\Omega))$ , $\varrho_t \in C([0, T^*]; H^1(\Omega))$ ;
$(\partial_t^k u, \partial_t^k \vartheta) \in L^2(0, T^*; H^{3-2k}(\Omega))^4$ for $$k=0,1,2$$ and $(\partial_t^k u, \partial_t^k \vartheta) \in L^\infty(0, T^*; H^{2-2k}(\Omega))^4$ for $$k=0,1$$ ;
$\varrho(x,t) > 0$ and $\vartheta(x,t) > 0$ in $Q_{T^*}$ . 此结果放宽了 Valli-Zajaczkowski 理论中对初始时间导数的高阶相容性要求，代价是时间正则性略低。

定理 1.6 (爆破判据): 在上述假设基础上，进一步假设边界数据不依赖于时间，且 $f \in L^1(0, T; L^\infty(\Omega)^3)$ 。若解在 $$t = T^*$$ 时刻爆破，当且仅当存在 $\delta > 0$ 使得：

\lim_{T \nearrow T^*} \int_0^T \left[ \|\nabla u(t)\|_{L^\infty(\Omega)} + \|\vartheta(t)\|_{L^\infty(\Omega)}^{5+\delta} \right] dt = \infty.

这意味着，只要速度梯度在时间上是 $$L^1$$ 可积的，且最大温度在时间上具有足够的可积性（指数 $5+\delta$ ），强解就可以延拓。这是针对带有真实进出流边界的全 Navier-Stokes-Fourier 系统的首个此类判据。

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

证明的核心技术障碍在于如何处理进出流边界条件下微分后的连续性方程所产生的边界项。具体而言，当对连续性方程 $\varrho_t + \text{div}(\varrho u) = 0$ 进行高阶空间微分时，会在进流边界 $\Gamma_{\text{in}}$ 上产生形如 $\int_{\Gamma_{\text{in}}} (u \cdot N) |\nabla^k \varrho|^2 dS$ 的边界积分项。由于 $u \cdot N < 0$ ，该项符号为正，无法像不可穿透边界那样通过自然消去来处理，反而可能破坏能量估计的闭合。

作者的创新点在于发展了一套专门针对进出流结构的边界估计技术。首先，利用进流边界处法向速度的负号，结合密度 $\varrho$ 的正性，通过对连续性方程沿特征线追踪（使用反向流映射 $$X$$ ），建立了密度的上下界估计（Lemma 2.1）：

\min\{\underline{\varrho}_0, \underline{\varrho}_B\} \exp\left(-\int_0^t \|\text{div} u\|_{L^\infty} dt'\right) \le \varrho(t,x) \le \max\{\overline{\varrho}_0, \overline{\varrho}_B\} \exp\left(\int_0^t \|\text{div} u\|_{L^\infty} dt'\right).

这一估计确保了密度的正性和有界性，从而使得在处理边界项时，可以将 $\varrho$ 视为良态系数。

其次，在动量方程的处理中，作者采用了类似于 Huang-Xin 的方法，但进行了关键修改以适应非零法向速度。他们引入了修正的速度变量 $$w = u - u_B$$ ，并利用 $$u_B$$ 的调和延拓性质简化了非线性项。对于高阶估计，作者没有直接使用 $D_t u - u \cdot \nabla u_B$ 作为测试函数（这在 $u \cdot N \neq 0$ 时不再满足齐次边界条件），而是转而依赖 Beale-Kato-Majda 型条件 $\nabla u \in L^1_t L^\infty_x$ 来控制流体变形的累积效应。这种控制允许作者在存在进出流的情况下闭合高阶 Sobolev 范数的估计。

此外，温度方程的处理利用了热传导项 $\kappa \Delta \vartheta$ 的正则化作用，并结合粘性耗散项 $|\nabla u|^2$ 的控制。温度判据中出现的指数 $5+\delta$ 反映了在缺乏 $L^\infty_t L^\infty_x$ 均匀界时，所需的最弱时间可积性水平。整个证明框架结合了半群方法、不动点定理用于局部存在性，以及精细的先验估计用于延拓性。

4. 与之前工作的比较

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

边界条件的本质差异: 大多数现有的条件正则性或爆破判据研究（如 Basarić, Feireisl, Mizerová [4], Feireisl, Wen, Zhu [12]）都假设不可穿透边界 $u \cdot N = 0$ 。在这种设定下，可以利用 $u \cdot N = 0$ 在边界积分中进行消去，从而得到基于 $L^\infty_t L^\infty_x$ 范数的延拓判据。然而，本文处理的真实进出流边界 ( $u \cdot N \neq 0$ ) 彻底破坏了这一机制。正如 Remark 1.7 所指出的，传统的测试函数技巧在此失效。因此，本文不得不引入更强的 Beale-Kato-Majda 型条件 $\nabla u \in L^1_t L^\infty_x$ 来替代单纯的 $L^\infty$ 界，这是针对开放系统的一个必要且自然的推广。
与 Valli-Zajaczkowski 理论的比较: Valli 和 Zajaczkowski [27] 建立了带有进出流边界的可压缩流体强解理论，但其要求初始数据满足由方程生成的额外 $$H^1_0$$ 相容性条件（即初始时间导数必须满足特定的边界匹配条件）。本文的定理 1.4 仅需基本的迹相容条件，这使得适用范围更广，能够容纳更多物理上合理的初始数据。代价是时间正则性从 $$C^1$$ 级别降低到 $$L^2_t H^1_x$$ 级别，但这足以支持后续的爆破判据分析。
与 Meliani 的 $$L^p-L^q$$ 理论的比较: Meliani [21] 在 $$L^p-L^q$$ 框架下建立了强解的存在性，但其参数范围受到严格限制（例如 $q > 3, p > \max\{\frac{2q}{q-1}, \frac{2q}{2q-3}\}$ ），且在 $$q=6$$ 时排除了端点情形。相比之下，本文在标准的 Hilbert 空间 $$H^s$$ 框架下工作，得到了更清晰的高阶 Sobolev 正则性 $$H^2$$ 和 $$H^3$$ ，并且不需要复杂的 Besov 空间嵌入。虽然 Meliani 的理论包含了磁流体力学耦合，但就纯 Navier-Stokes-Fourier 系统而言，本文的结果提供了更直接的 Sobolev 正则性描述。
与 Abbatiello-Meliani 的延拓结果的比较: Abbatiello 和 Meliani [3] 为等温 (barotropic) 可压缩 Navier-Stokes 系统建立了带有进出流边界的延拓判据。他们的判据涉及 $\varrho^{-1}, u$ 或 $\nabla \varrho$ 的范数。然而，等温模型不包含温度方程，因此无法捕捉压力功、粘性加热和热传导之间的耦合。本文将这一思想推广到了完整的热力学系统 (Navier-Stokes-Fourier)，并首次将温度 $\vartheta$ 纳入爆破判据中，体现了热力学效应在开放系统中的重要性。
关于真空状态的假设: 与某些允许初始真空的判据（如 Huang, Li, Xin [16]）不同，本文要求初始密度严格为正 ( $\inf \varrho_0 > 0$ )。作者在 Remark 1.9 中解释，这是因为在进出流边界条件下，密度的高阶估计直接依赖于连续性方程在边界上的行为，而正性是传播密度上下界的关键。这一假设并非纯粹的技术便利，而是由进出流结构所决定的内在要求。

5. 关于洛伦兹空间 $L^{3,\infty}(\mathbb{R}^3)$ 中 Navier-Stokes 方程解的唯一性与正则性

On the uniqueness and the regularity of solutions to the Navier Stokes in the Lorentz space $L^{3,\infty}(\mathbb{R}^3)$

作者: Ramzi May

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14688v1

摘要翻译: 我们证明了在空间

Comment: 21 pages

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

1. 研究对象和背景

本文研究三维全空间 $\mathbb{R}^3$ 中的不可压缩 Navier-Stokes 方程:

\begin{cases} \partial_t u - \Delta u + P\nabla \cdot (u \otimes u) = 0 \\ \nabla \cdot u = 0 \\ u(0, x) = u_0(x) \end{cases}

其中 $$u(t,x)$$ 为速度场, $$P$$ 为 Leray 投影算子. 该方程具有标度不变性: 若 $$u$$ 是解, 则对任意 $\lambda > 0$ , $u_\lambda(t,x) = \lambda u(\lambda^2 t, \lambda x)$ 也是对应初始数据 $u_{0,\lambda}(x) = \lambda u_0(\lambda x)$ 的解. 这一性质引出了临界空间的概念: 一个函数空间 $$X$$ 被称为 Navier-Stokes 方程的临界空间, 如果其范数在标度变换下保持不变, 即 $\|u_{0,\lambda}\|_X = \|u_0\|_X$ . 重要的临界空间包括齐次 Sobolev 空间 $\dot{H}^{1/2}(\mathbb{R}^3)$ , Lebesgue 空间 $L^3(\mathbb{R}^3)$ 以及 Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ .

利用 Duhamel 原理, Navier-Stokes 方程可以改写为积分方程:

u(t) = e^{t\Delta} u_0 - B(u,u)(t),

其中 $e^{t\Delta}$ 是热半群, $$B$$ 是双线性算子, 定义为:

B(u,v)(t) = \int_0^t P\nabla e^{(t-s)\Delta} (u \otimes v)(s) ds.

通过在适当的 Banach 空间中建立双线性算子 $$B$$ 的连续性, 可以利用 Banach 不动点定理证明解的存在性.

Furioli, Lemarie-Rieusset 和 Teraneo [7] 在 1998 年证明了 $C([0,T], L^3(\mathbb{R}^3))$ 中解的唯一性. 值得注意的是, 双线性算子 $$B$$ 在 $C([0,T], L^3(\mathbb{R}^3))$ 上并不连续, 这一点后来由 Oru [20] 指出. Yves Meyer [19] 随后证明了 $$B$$ 在 $L^\infty([0,T], L^{3,\infty}(\mathbb{R}^3))$ 上是连续的, 从而提供了 Furioli-Lemarie-Rieusset-Teraneo 唯一性结果的直接证明. Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ 的一个重要特性是它包含非平凡的齐次函数, 例如 $f(x) = |x|^{-1} \in L^{3,\infty}(\mathbb{R}^3)$ , 这使得它成为研究自相似解的自然设置.

2007 年, Lemarie-Rieusset [12] 改进了 Chemin [4] 的唯一性结果, 证明了在 $L^p([0,T], L^q(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ 中当 $$p > 2$$ 且 $$q > 3$$ 时解的唯一性. 作者在前作 [17] 中将此结果推广到极限情形 $$q=3$$ . 由于 $L^3(\mathbb{R}^3)$ 是极限 Besov 空间 $B_{\infty}^{-1,\infty}(\mathbb{R}^3)$ 的子空间, 这提供了 Furioli-Lemarie-Rieusset-Teraneo 唯一性结果的另一种独立证明.

上述结果自然引出了一个问题: 当空间 $L^3(\mathbb{R}^3)$ 被替换为更大的 Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ 时, [17] 中的唯一性结果是否仍然成立? 本文旨在对此问题给出肯定的回答. 具体地, 作者在空间

L^p([0,T], \widetilde{L}^{3,\infty}(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3)), \quad p > 2,

中证明了唯一性, 其中 $\widetilde{L}^{3,\infty}(\mathbb{R}^3)$ 表示 Schwartz 空间 $\mathcal{S}(\mathbb{R}^3)$ 在 Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ 中的闭包. 此外, 作者还建立了属于这些空间的解的正则性性质.

2. 主要定理或主要结果

本文的主要结果是定理 1.1, 陈述如下:

定理 1.1. 设 $$T > 0$$ 且 $$p > 2$$ . Navier-Stokes 方程在空间

L^p([0,T], \widetilde{L}^{3,\infty}(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))

中至多有一个解. 此外, 如果 $$u$$ 是这样的一个解, 则它具有以下正则性性质: (1) $u \in C^\infty((0,T] \times \mathbb{R}^3)$ ; (2) 对任意 $$s > 0$$ , $u \in C_{0, s+1/2}([0,T], B_{\infty}^{s,\infty}(\mathbb{R}^3))$ ; (3) $u \in C_{0, 1/2}([0,T], L^\infty(\mathbb{R}^3))$ ; (4) 对任意 $\alpha > 0$ , $(-\Delta)^{\alpha/2} u \in C_{0, \alpha/2}([0,T], L^{3,\infty}(\mathbb{R}^3))$ .

这里 $C_{0,s}([0,T], X)$ 表示满足 $\lim_{t \to 0^+} t^{s/2} \|v(t)\|_X = 0$ 的函数 $v \in C((0,T], X)$ 的空间.

备注 1.1. 定理 1.1 中的唯一性断言可以从 [17, 定理 1.3] 推导出来, 因为 $L^{3,\infty}(\mathbb{R}^3) \subset M_{r,3}(\mathbb{R}^3)$ , $r \in (2,3)$ . 然而, 本文提供的证明更加透明, 避免了 [17] 论证中涉及的大量技术细节.

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

本文的证明框架主要依赖于以下几个关键技术和工具:

Littlewood-Paley 分解与 Besov 空间: 作者引入了 Littlewood-Paley 分解, 定义了算子 $$S_j$$ 和 $\Delta_j$ , 并基于此定义了非齐次 Besov 空间 $B_X^{s,q}(\mathbb{R}^n)$ . 这些空间在分析 Navier-Stokes 方程的正则性时起着核心作用. 特别是, 作者利用了 Besov 空间的插值性质 (引理 2.1), 该性质允许在不同正则性阶数的空间之间进行估计.
修改版 Bony 乘积分解: 作者引入了一种修改版的 Bony 乘积分解, 定义了算子 $\pi_a(f,g) = \sum_{j=3}^\infty S_{j-2} f \Delta_j g$ 和 $\pi_b(f,g) = \sum_{j=-1}^\infty S_{j+2} f \Delta_j g$ . 这种分解比经典版本更适合处理 Lorentz 空间中的估计. 引理 2.2 给出了这些双线性算子在 Besov 空间中的连续性, 这是证明非线性项估计的关键.
热半群的平滑效应: 作者详细讨论了热半群 $e^{t\Delta}$ 在 Besov 空间中的正则化效应 (引理 2.3). 特别是, 对于 $f \in B_X^{-s,\infty}$ , 有估计
$\sup_{0 < t \le T} t^{s/2} \|e^{t\Delta} f\|_X \lesssim \|f\|_{B_X^{-s,\infty}}.$
这一估计在处理小时间行为时至关重要.
Lorentz 空间的性质: 作者回顾了 Lorentz 空间 $L^{p,q}(\mathbb{R}^n)$ 的一系列重要性质 (命题 2.2), 包括它们的嵌入关系、卷积不等式以及点态乘积的估计. 特别是, 作者使用了改进的 Hardy-Littlewood-Sobolev 不等式 (引理 2.5), 该不等式在 Lorentz 空间中成立, 形式为:
$\|g\|_{L^{p_2,q}} \lesssim c_{n,\alpha} \|f\|_{L^{p_1,q}},$
其中 $g(x) = \int_{\mathbb{R}^n} \frac{f(y)}{|x-y|^{n-\alpha}} dy$ , $p_2 = \frac{np_1}{n-\alpha p_1}$ .
精细的 Sobolev 型不等式: 作者证明了一个精细的 Sobolev 型不等式 (引理 2.7), 该不等式可以视为 Gérard, Meyer 和 Oru [9] 改进的 Sobolev 不等式的非齐次版本. 具体地, 对于 $f \in W^{1,p}(\mathbb{R}^n) \cap B_{\infty}^{-1,\infty}(\mathbb{R}^n)$ , 有:
$\|f\|_{2p} \lesssim \|f\|_{W^{1,p}}^{1/2} \|f\|_{B_{\infty}^{-1,\infty}}^{1/2}.$
这个不等式通过将函数分解为低频和高频部分, 并利用 Hardy-Littlewood 极大函数的有界性来证明.
线性积分算子的有界性: 为了理解主定理中 $$p > 2$$ 的限制, 作者证明了线性算子
$L(v)(t) = \int_0^t \frac{v(\tau)}{\sqrt{t-\tau}\sqrt{\tau}} d\tau$
从 $L^p([0,T], \mathbb{R})$ 到自身的有界性 (引理 2.8), 其中 $2 < p < \infty$ . 这一结果通过 Minkowski 不等式得到.

通过这些工具的有机结合, 作者能够在 Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ 中建立所需的估计, 从而证明唯一性和正则性结果.

4. 与之前工作的比较

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

与 Furioli, Lemarie-Rieusset 和 Teraneo [7] 的比较: Furioli 等人证明了 $C([0,T], L^3(\mathbb{R}^3))$ 中解的唯一性. 本文的结果将此唯一性推广到了更大的空间 $L^p([0,T], \widetilde{L}^{3,\infty}(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ , $$p > 2$$ . 由于 $L^3(\mathbb{R}^3) \subset L^{3,\infty}(\mathbb{R}^3)$ , 本文的结果严格包含了 Furioli 等人的结果. 此外, 本文还建立了详细的正则性估计, 而 Furioli 等人的工作主要关注唯一性.
与 Yves Meyer [19] 的比较: Meyer 证明了双线性算子 $$B$$ 在 $L^\infty([0,T], L^{3,\infty}(\mathbb{R}^3))$ 上的连续性, 从而提供了 Furioli-Lemarie-Rieusset-Teraneo 唯一性结果的直接证明. 本文的工作可以看作是对 Meyer 结果的进一步推广和精细化. 本文不仅考虑了 $L^\infty$ 时间范数, 还考虑了更一般的 $$L^p$$ 时间范数 ( $$p > 2$$ ), 并且结合了 $C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ 的条件, 从而获得了更强的正则性结论.
与 Lemarie-Rieusset [12] 和 Chemin [4] 的比较: Lemarie-Rieusset 改进了 Chemin 的唯一性结果, 证明了在 $L^p([0,T], L^q(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ 中当 $$p > 2$$ 且 $$q > 3$$ 时解的唯一性. 作者在前作 [17] 中将此结果推广到极限情形 $$q=3$$ . 本文将空间从 $L^3(\mathbb{R}^3)$ 进一步推广到 Lorentz 空间 $L^{3,\infty}(\mathbb{R}^3)$ . 这一推广是自然的, 因为 $L^{3,\infty}(\mathbb{R}^3)$ 是 $L^3(\mathbb{R}^3)$ 的更大空间, 并且包含非平凡的齐次函数, 使得它成为研究自相似解的理想设置.
与前作 [17] 的比较: 作者在前作 [17] 中证明了 $L^p([0,T], L^3(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ , $$p > 2$$ 中解的唯一性. 本文的结果将此空间替换为 $L^p([0,T], \widetilde{L}^{3,\infty}(\mathbb{R}^3)) \cap C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ , $$p > 2$$ . 虽然唯一性断言可以从 [17, 定理 1.3] 推导出来 (因为 $L^{3,\infty}(\mathbb{R}^3) \subset M_{r,3}(\mathbb{R}^3)$ , $r \in (2,3)$ ), 但本文提供了更透明的证明, 避免了 [17] 中的大量技术细节. 此外, 本文还建立了详细的正则性估计, 这是 [17] 中没有明确给出的.
假设与适用范围: 本文的主要假设是 $$p > 2$$ 和初始数据属于 $C([0,T], B_{\infty}^{-1,\infty}(\mathbb{R}^3))$ . 这一条件比单纯要求初始数据属于 $L^{3,\infty}(\mathbb{R}^3)$ 更强, 但它允许获得更精细的正则性估计. 特别地, $B_{\infty}^{-1,\infty}(\mathbb{R}^3)$ 是一个较大的空间, 包含了 $L^3(\mathbb{R}^3)$ 和 $L^{3,\infty}(\mathbb{R}^3)$ , 因此本文的结果适用于较宽的初始数据类.
代价与局限: 本文的一个局限是它只证明了唯一性, 而没有直接证明存在性. 存在性问题需要额外的论证, 通常依赖于 Kato 的方法或其他构造性方法. 此外, 本文的正则性估计依赖于初始数据属于 $B_{\infty}^{-1,\infty}(\mathbb{R}^3)$ , 这对于某些特定的初始数据可能是一个较强的要求. 然而, 这一条件在理论分析中是常见的, 并且允许获得最优的正则性估计.

6. 三维磁流体力学边界层方程的线性不适定性

Linear ill-posedness of three-dimensional MHD boundary layer equations

作者: Mingxue Zhang, Zhonger Wu

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14532v1

摘要翻译: 我们证明了三维电阻磁流体力学（MHD）边界层方程在切向 Sobolev 空间中的线性不适定性。我们假设初始切向速度分量的适当线性组合存在一个非退化临界点，并且两个初始切向磁场分量在同一点具有二阶零点。与相应的二维 MHD 结果相比，我们的构造仅要求该二阶零点条件，而非更高阶的磁场剪切退化。证明结合了三维 Prandtl 临界层构造以及对二维磁场商的向量修正。该修正同时抵消了两个感应方程中的主导项，同时保留了磁场散度约束。由此得到的近似解的增长行为为 $\exp(\sigma t/\sqrt{\varepsilon})$ ，其残差的前因子量级为 $O(\varepsilon^{-3/8})$ ，在对数因子范围内成立。随后通过 Duhamel 论证，我们得到了对于每个切向导数损失 $\mu<\frac18$ 的线性不适定性。因此，当稳定化的切向磁场发生退化时，三维 Prandtl 不稳定性依然存在。

Comment: 无

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

1. 研究对象和背景

本文研究的是三维电阻磁流体动力学（MHD）边界层方程的线性稳定性问题。具体而言，考虑定义在区域 $\Omega = \{(t, x, y, z) : t > 0, (x, y) \in \mathbb{T}^2, z \in \mathbb{R}^+\}$ 上的如下方程组：

\begin{cases} (\partial_t + u \cdot \nabla_h + w\partial_z - \partial_z^2)u - (b \cdot \nabla_h + h\partial_z)b + \nabla_h p = 0, \\ (\partial_t + u \cdot \nabla_h + w\partial_z - \partial_z^2)b - (b \cdot \nabla_h + h\partial_z)u = 0, \\ \nabla_h \cdot u + \partial_z w = 0, \\ \nabla_h \cdot b + \partial_z h = 0, \\ (u, w)|_{z=0} = 0, \\ (\partial_z b, h)|_{z=0} = 0, \\ \lim_{z \to +\infty}(u, b) = (U_E, B_E), \end{cases}

其中 $$u = (u, v)$$ 和 $$b = (b, g)$$ 分别是切向速度和磁场分量， $$w$$ 和 $$h$$ 是法向速度及磁场分量， $\nabla_h = (\partial_x, \partial_y)$ 为切向梯度。外部流场 $$(U_E, B_E)$$ 满足伯努利定律。当磁场消失时，该方程组退化为经典的普朗特（Prandtl）方程。

普朗特方程的研究始于 Oleinik 在 1963 年的工作，她在单调性条件 $\partial_z u > 0$ 下利用 Crocco 变换建立了二维普朗特方程经典解的局部适定性。随后，Alexandre, Wang, Xu 和 Yang [1] 以及 Masmoudi 和 Wong [17] 分别在加权 Sobolev 空间中直接证明了相应结果。然而，当剪切流存在非退化临界点时，普朗特方程在 Sobolev 空间中表现出强烈的线性不适定性 [4]。对于三维普朗特方程，Liu, Wang 和 Yang [12] 指出次级流的存在导致了线性不稳定性，这种不稳定性源于切向流动方向随法向变量的变化。

与普朗特方程相比，MHD 边界层方程由于磁场中额外的导数损失而更加复杂。已有研究表明，若初始切向磁场是非退化的，则二维 MHD 边界层方程在 Sobolev 空间中的局部适定性成立 [7, 11, 15]。近期，Wu [23] 建立了三维非线性 MHD 边界层方程在 Sobolev 空间中的适定性。这些结果暗示磁场对边界层方程具有稳定作用。另一方面，文献 [14] 证明，若切向磁场在切向速度剪切的非退化临界点处充分退化，则二维 MHD 边界层方程会出现线性不稳定性。这表明当磁场退化时，其稳定效应可能会失效。然而，三维情形下是否存在类似的不稳定性此前一直是一个开放问题。

本文旨在填补这一空白，建立三维电阻 MHD 边界层系统的线性不稳定性。作者假设外层流矢量 $$U_E$$ 和 $$B_E$$ 为常数，从而压力梯度为零。围绕剪切流解 $$(U_s(t, z), 0, B_s(t, z), 0)$$ 进行线性化，得到线性化方程组 (1.5)。本文的核心目标是证明，在特定的磁场退化条件下，该线性化系统在某些 Sobolev 范数意义下是不适定的，即解算子在切向导数损失意义下是无界的。

2. 主要定理或主要结果

本文的主要结果是定理 1.1，它证明了三维电阻 MHD 边界层方程在线性化后的切向 Sobolev 空间中的不适定性。

定理 1.1. 假设剪切流解满足以下正则性条件：

U_s - U_E \in \bigcap_{i=0}^{1} C^i(\mathbb{R}^+; W^{4-2i, \infty}_\alpha(\mathbb{R}^+) \cap H^{4-2i}_\alpha(\mathbb{R}^+)),

B_s - B_E \in \bigcap_{i=0}^{2} C^i(\mathbb{R}^+; W^{6-2i, \infty}_\alpha(\mathbb{R}^+) \cap H^{6-2i}_\alpha(\mathbb{R}^+)).

假设存在 $$z_0 > 0$$ 使得以下条件成立：

V'_s(z_0) U''_s(z_0) \neq U'_s(z_0) V''_s(z_0),

B_s(z_0) = \partial_z B_s(z_0) = 0.

此外，交换 $$x$$ 和 $$y$$ 后，假设进一步有：

U'_s(z_0) \neq 0, \quad a := \frac{V'_s(z_0)}{U'_s(z_0)} \in \mathbb{Q}.

那么，存在 $\sigma > 0$ 使得对于任意 $\delta > 0$ , $m \ge 0$ 和 $\mu \in [0, \frac{1}{8})$ ，有：

\sup_{0 \le s \le t \le \delta} \| e^{-\sigma(t-s)/\sqrt{\varepsilon}} |\partial_h| T(t, s) \|_{L(H^m, H^{m-\mu})} = +\infty,

其中 $\partial_h$ 表示切向导数 $\partial_x$ 或 $\partial_y$ ， $$T(t, s)$$ 是线性化方程 (1.5) 的解算子。

技术细节解读：

增长率： 近似解的增长率为 $\exp(\sigma t / \sqrt{\varepsilon})$ 。这里的 $\sigma$ 取决于临界点处的二阶导数大小，具体地， $\sigma^*(t) = -|\frac{\partial_z^2 w^a_s(t, f(t))}{2}|^{1/2} \Im \tilde{\tau}$ 。
残差估计： 构造的近似解的残差 $R_\varepsilon$ 满足 $\|R_\varepsilon(t)\|_{L^2_\alpha} \le C_0 \varepsilon^{-3/8} L_\varepsilon^N e^{\sigma_0 t / \sqrt{\varepsilon}}$ ，其中 $L_\varepsilon = 1 + |\log \varepsilon|$ 。
导数损失： 通过 Duhamel 论证，得出的不适定性结论适用于任何切向导数损失 $\mu < \frac{1}{8}$ 。这个范围弱于二维普朗特方程的结果 ( $\mu < \frac{1}{2}$ ) 以及三维普朗特方程和二维 MHD 方程的结果 ( $\mu < \frac{1}{4}$ )。作者指出这并非最优范围，而是由残差阶 $O(\varepsilon^{-3/8})$ 和 Duhamel 公式中的 $\sqrt{\varepsilon}$ 因子共同决定的。
有理数假设： 参数 $$a = V'_s(z_0)/U'_s(z_0)$$ 必须是有理数。这是为了确保高频模态在环面 $\mathbb{T}^2$ 上的周期性。如果 $a = p/q \in \mathbb{Q}$ ，取 $\varepsilon_n = (q_n)^{-1}$ ，则模态 $e^{i\varepsilon_n^{-1}(y-ax)} = e^{i(q_n y - p_n x)}$ 在两个切向变量上都是周期的。

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

证明的核心在于构造一个具有特定增长行为的近似解，并利用该近似解的残差估计来推导解算子的无界性。主要技术障碍和创新点如下：

三维普朗特临界层构造与向量修正的结合： 在三维普朗特方程中，Liu, Wang 和 Yang [12] 利用了临界层结构来构造不稳定模态。然而，在 MHD 方程中，磁场引入了额外的耦合项。简单的标量商构造（如在二维情况 [11] 中使用的那样）无法同时消去两个感应方程中的主导项。本文的关键创新是引入了一个向量修正项。定义 $A_\varepsilon(t, z) = \lambda_\varepsilon(t) + w^a_s(t, z)$ 和 $Q_\varepsilon(t, z) = W^{reg}_\varepsilon(t, z) + W^{sl}_\varepsilon(t, z)$ ，并令 $\rho_\varepsilon(t, z) = Q_\varepsilon/A_\varepsilon$ 。尽管分子和分母在临界层内均为 $\sqrt{\varepsilon}$ 阶，但 $\rho_\varepsilon$ 是正则的，因为 Heaviside 函数的跳跃被剪切层轮廓精确抵消。切向磁场振幅定义为：
$B_\varepsilon := \rho_\varepsilon \partial_z B_s + \frac{c_s}{A_\varepsilon} (U_\varepsilon - \rho_\varepsilon \partial_z U_s),$
其中 $$c_s(t, z) = g_s(t, z) - a b_s(t, z)$$ 。这个表达式中的第二项 $\frac{c_s}{A_\varepsilon} D_\varepsilon$ （其中 $D_\varepsilon = U_\varepsilon - \rho_\varepsilon \partial_z U_s$ ）是关键的向量修正。它与 $$(-a, 1)$$ 正交，从而保证了磁场散度约束 $\nabla_h \cdot b + \partial_z h = 0$ 的自动满足。同时，直接计算表明，这个修正项精确地消去了感应方程中的频率项：
$-A_\varepsilon B_\varepsilon + Q_\varepsilon \partial_z B_s + c_s U_\varepsilon - c_s \rho_\varepsilon \partial_z U_s = 0.$
残差的精细估计： 将上述构造代入线性化方程，得到残差 $R_\varepsilon = (R^1_\varepsilon, R^2_\varepsilon, R^3_\varepsilon, R^4_\varepsilon)$ 。速度方程的残差主要由普朗特动量残差 $E_\varepsilon P_\varepsilon$ 和一个额外项组成：
$(R^1_\varepsilon, R^2_\varepsilon) = E_\varepsilon \left[ i(\partial_t - \partial_z^2) U_\varepsilon + \varepsilon^{-1} \left( -A_\varepsilon U_\varepsilon + Q_\varepsilon \partial_z U_s + \frac{c_s^2}{A_\varepsilon} D_\varepsilon \right) \right].$
由于 $$c_s$$ 在临界点 $$z_0$$ 处为零（因为 $$B_s(z_0)=0$$ 且 $$a$$ 的定义），该项的贡献受到控制。感应方程的残差简化为：
$(R^3_\varepsilon, R^4_\varepsilon) = i E_\varepsilon (\partial_t - \partial_z^2) B_\varepsilon.$
作者通过仔细分析临界层内的缩放行为（令 $Z = (z-f(t))/\varepsilon^{1/4}$ ），证明了残差的 $$L^2$$ 范数为 $O(\varepsilon^{-3/8} L_\varepsilon^N e^{\sigma_0 t / \sqrt{\varepsilon}})$ 。这里的 $\varepsilon^{-3/8}$ 来源于对 $B_\varepsilon$ 的时间导数和法向二阶导数的估计，特别是考虑到 $B_\varepsilon$ 中包含的 $\rho_\varepsilon$ 项及其导数在临界层的奇异行为。
Duhamel 论证与不适定性推导： 利用解析函数空间 $E_{\alpha, \beta}$ 中的适定性结果（命题 1.1），可以将线性化方程的解表示为齐次解加上由残差驱动的积分项。由于近似解在 $$H^m$$ 范数下保持有界（乘以指数衰减因子后），而其在 $H^{m-\mu}$ 范数下的增长由残差积分决定。具体来说，解算子 $$T(t, s)$$ 的作用可以分解为：
$T(t, s)(u_0, b_0) = \text{Homogeneous Part} + \int_s^t T(t, \tau) R_\varepsilon(\tau) d\tau.$
通过分析积分项在高频模态下的贡献，并结合残差的阶 $\varepsilon^{-3/8}$ 和时间积分产生的 $\sqrt{\varepsilon}$ 因子，最终得出解算子在 $L(H^m, H^{m-\mu})$ 中无界，只要 $\mu < 1/8$ 。

4. 与之前工作的比较

本文工作与既有文献的比较主要体现在以下几个方面：

与二维 MHD 结果的比较 ([14])： 文献 [14] 研究了二维电阻 MHD 边界层方程的线性不稳定性。在该工作中，为了克服磁场带来的稳定效应，作者要求切向磁场在速度剪切临界点处具有六阶退化（即磁场及其前五阶导数均为零）。相比之下，本文在三维情形下仅需假设两个切向磁场分量在临界点处具有二重零点（即 $B_s(z_0) = \partial_z B_s(z_0) = 0$ ）。这是一个显著的条件放宽。这是因为三维结构中，向量修正项能够更有效地处理感应方程中的耦合，从而允许更强的磁场存在而不破坏不稳定性机制。
与三维普朗特方程结果的比较 ([12])： Liu, Wang 和 Yang [12] 证明了三维普朗特方程在 Sobolev 空间中的线性不适定性，其导数损失范围为 $\mu < 1/4$ 。本文将这一结果推广到了 MHD 情形，但导数损失范围缩小至 $\mu < 1/8$ 。这一减弱是由于 MHD 方程中磁场扩散导致背景剪切流随时间变化，从而产生了额外的感应残差。此外，MHD 方程中磁场与速度的耦合使得构造更复杂的近似解成为必要，进而影响了残差的阶数。尽管如此，本文的结果确认了三维普朗特不稳定性在磁场退化时依然存在，即磁场的稳定效应在退化情况下会失效。
与非电阻 MHD 模型的比较 ([11])： 文献 [11] 研究了非电阻 MHD 边界层方程。在非电阻模型中，背景磁场剪切是定常的，这使得构造过程相对简单。而在本文考虑的电阻模型中，磁场扩散使得背景磁场剪切随时间演化，这引入了额外的复杂性。例如，背景磁场 $$B_s(t, z)$$ 不再是时间独立的，因此在构造近似解时必须考虑其时间导数对残差的影响。这也是为什么本文需要更细致的残差估计，并最终导致导数损失范围较小的原因之一。
与适定性结果的对比： 已有工作 [7, 11, 15, 23] 在磁场非退化或特定结构假设下建立了 MHD 边界层方程的适定性。本文的结果与之形成互补：它展示了当磁场在关键位置退化时，即使是在线性水平上，Sobolev 空间中的适定性也会失败。这强调了磁场结构对于边界层方程稳定性的重要性。特别地，Chen, Li 和 Yang [2] 在部分解析空间中建立了线性化 MHD 方程的适定性，但其假设排除了本文考虑的退化磁场情形。因此，本文填补了退化磁场下线性稳定性理论的空白。

综上所述，本文不仅在维度上从二维推广到三维，而且在磁场退化条件的强度上做出了实质性改进（从六阶退化降至二重零点），揭示了三维 MHD 边界层方程在磁场退化时的内在不稳定性机制。

7. 二维惯性磁流体力学中的线性与非线性不稳定性

Linear and Nonlinear Instability in Two-Dimensional Inertial Magnetohydrodynamics

作者: Hyungjun Choi

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14366v1

摘要翻译: 我们研究了环面上二维惯性磁流体力学系统某些稳态的谱不稳定性与非线性不稳定性。对于速度与磁场平行的平衡态，线性化问题可归约为傅里叶系数的矩阵递推关系，而稳定的矩阵连分式能够产生光滑的本征模。对于一个纯磁场族，线性化算子的平方在每个横向傅里叶块上具有紧自伴约化，该约化可对整个不稳定谱进行分类。在一个不变对称类中，我们确定了不稳定分支以及主导本征空间下方的一个谱间隙。该谱间隙进而允许构建近似轨迹，从而证明了在对数时间尺度上的非线性不稳定性。本研究受物理文献启发，在这些文献中，纯磁场撕裂模的非线性不稳定性被用来解释磁重联现象。

Comment: 25 pages, comments welcome

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

1. 研究对象和背景

本文研究在平坦环面 $\mathbb{T}^2 = (\mathbb{R}/2\pi\mathbb{Z})^2$ 上定义的二维惯性磁流体力学（Inertial MHD）系统的谱不稳定性与非线性不稳定性。该模型是带有电子惯性的薄层磁流体力学的二维约化，其动力学方程为：

\begin{aligned} \partial_t \omega + \nabla^\perp \psi \cdot \nabla \omega + \nabla^\perp \Phi \cdot \nabla F &= 0, \\ \partial_t F + \nabla^\perp \psi \cdot \nabla F &= 0, \\ \Delta \psi &= \omega, \\ (1 - \Delta) \Phi &= F. \end{aligned}

其中， $\nabla^\perp f = (-\partial_{x_2}f, \partial_{x_1}f)$ ， $\psi$ 为流体流函数， $\omega$ 为涡量， $\Phi$ 为磁通量函数。广义通量 $F = (1-d_e^2 \Delta)\Phi$ 引入了电子惯性效应，这里电子皮层深度 $$d_e$$ 被归一化为 1。

物理背景方面，磁重联与弛豫过程是研究惯性 MHD 的主要动机。在耗散磁流体力学中，重联通常与磁扩散相关，但在实验室和天体等离子体中，重联发生的时间尺度往往短于仅由碰撞电阻率所能解释的范围。在包含电子惯性的模型中，输运通量 $$F$$ 与磁通量 $\Phi$ 不同。虽然 $\Phi$ 的等值线描述了平面磁场线，但由于 $$F$$ 是被输运的量，这些磁场线并不一定随流体运动。Ottaviani 和 Porcelli 在该模型中研究了与磁重联相关的增长扰动（撕裂模），发现在早期非线性演化中存在加速增长和电流层变窄的现象。这些物理现象促使作者对稳态解进行严格的数学不稳定性分析。

在数学设定上，文章主要考察两类一维稳态解族。第一类是具有对齐速度和磁场的平衡态：

\psi^* = \sin(kx_1), \quad \Phi^* = a \sin(kx_1),

其中 $k \geq 2$ 为整数， $a \in \mathbb{R}$ 。第二类是纯磁正弦波族：

\psi^* = 0, \quad \Phi^* = \sin(kx_1),

其中 $k \geq 1$ 为整数。在这两种情况下， $\psi^*$ 和 $\Phi^*$ 都是拉普拉斯算子的特征函数，具有共同特征值 $$-k^2$$ 。因此， $\omega^* = \Delta \psi^*$ 和 $F^* = (1 - \Delta)\Phi^*$ 也是 $\sin(kx_1)$ 的常数倍，且所有非线性泊松括号项恒为零。

尽管存在关于理想流体稳定性的诸多文献，例如 Liao, Lin, 和 Zhu 对 Kelvin-Stuart 磁岛的研究，以及 Constantin 和 Hu 对另一种平面惯性模型的全局正则性证明，但针对本特定系统（特别是包含电子惯性的耦合结构）的严格线性及非线性不稳定性理论尚属空白。本文旨在填补这一空白，通过构造具体的不稳定特征模式和证明非线性逃逸，从数学上确立该类稳态的不稳定性，从而为理解磁重联机制提供坚实的数学基础。

2. 主要定理或主要结果

本文的主要结果分为三个部分，分别对应两个稳态族的线性不稳定性以及纯磁族的非线性不稳定性。

定理 2.1 (对齐平衡态的线性不稳定性): 设 $k \geq 2$ 为整数，且参数 $a \in \mathbb{R}$ 满足条件：

$$ a^2(1 + k^2) < 2. $$

则稳态 $\psi^* = \sin(kx_1), \Phi^* = a \sin(kx_1)$ 是线性不稳定的。具体而言，其线性化算子存在一个正实特征值 $\lambda > 0$ 和一个光滑的特征函数。该特征函数在横向傅里叶模式 $e^{ix_2}$ 下具有指数衰减的傅里叶系数。

命题 3.5, 3.7, 3.9 及定理 3.11 (纯磁族的谱分类与间隙): 对于纯磁平衡态 $\psi^* = 0, \Phi^* = \sin(kx_1)$ ，线性化算子的平方在每个横向傅里叶块上可约化为紧自伴算子。不稳定性谱完全由这些紧自伴算子的负特征值决定。

Fredholm 本质谱: 位于虚轴上的区间 $i[-k\sqrt{1+k^2}, k\sqrt{1+k^2}]$ 。
不稳定分支: 在不稳定对称类中，对于每个 $1 \leq n < k$ ，恰好存在一个简单的正特征值分支。
定量界: 给出了标量连分式公式及其定量界限。
谱间隙: 在主导特征空间下方存在正的谱间隙。特别地，当 $$k=1$$ 时，不存在增长谱值。

定理 4.2 (纯磁平衡态的非线性不稳定性): 对于任意 $k \geq 2$ ，纯磁平衡态是非线性不稳定的。更精确地，对于每个整数 $s \geq 4$ ，若初始扰动沿不变对称类中的归一化主导特征函数方向，大小为 $\epsilon$ ，则在时间 $T_\epsilon$ 内，其在 Sobolev 范数 $H^{s-1} \times H^s$ 中将达到固定的小量 $\delta$ （ $\delta$ 独立于 $\epsilon$ ）。该时间尺度为对数型：

T_\epsilon = \Lambda_k^{-1} \log(\delta/\epsilon),

其中 $\Lambda_k$ 是该对称类中最大的增长率。

需要注意的是，这里的非线性不稳定性是指 Sobolev 范数下的固定阶逃逸（fixed order-one escape），并未直接证明磁拓扑变化或电流片的形成。将这种不稳定性直接与磁重联过程联系起来仍是一个开放问题。

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

本文的证明策略结合了矩阵连分式方法、紧算子谱理论和高阶近似轨迹技术。

1. 矩阵连分式方法与稳定子空间构造 (针对对齐平衡态) 处理对齐平衡态 $\psi^* = \sin(kx_1), \Phi^* = a \sin(kx_1)$ 的核心困难在于线性化后的傅里叶系数递推关系是矩阵值的，而非标量的。当 $a \neq 0$ 时，速度场与磁场耦合导致递推矩阵 $$M_j$$ 非对角化。作者借鉴了 Meshalkin 和 Sinai 在 Kolmogorov 流中使用的连分式方法。首先，通过正常模假设 $(\omega, F)^T = e^{\lambda t} e^{ix_2} \sum v_j \cos(jkx_1)$ ，将线性化方程转化为近邻递推关系：

\mu b_j = M_{j-1}b_{j-1} - M_{j+1}b_{j+1}.

为了构造光滑特征函数，必须找到一组指数衰减的解。作者利用离散指数二重性（discrete exponential dichotomy）的粗糙性理论，证明了对于足够大的 $$j$$ ，存在唯一的二维稳定子空间。随后，通过构造收敛的矩阵连分式 $R(\mu)$ 来显式地表征这个稳定子空间。关键在于证明映射 $T_j(R) = (\mu I + M_{j+1}R)^{-1}M_{j-1}$ 在某个加权范数下是压缩映射。这依赖于矩阵 $$M_j$$ 的渐近性质 $M_j \to M_\infty$ 以及 $M_j^{-1}$ 对称部分的均匀正定性（引理 2.2）。通过比较 $\mu \to 0$ 和 $\mu \to \infty$ 时匹配行列式 $f(\mu) = \det\{\mu I + M_1 R(\mu)\}$ 的符号变化，利用介值定理证明了正特征值的存在。小 $\mu$ 估计利用了微分映射 $D\Gamma_j$ 的收缩性，而大 $\mu$ 估计则直接利用连分式的渐近展开。

2. 紧自伴约化与谱间隙 (针对纯磁平衡态) 对于纯磁平衡态 $\psi^* = 0$ ，由于背景速度为零，线性化算子 $$L_k$$ 的结构简化为：

L_k = -k \cos(kx_1) \partial_{x_2} \begin{pmatrix} 0 & A_k \\ B_k & 0 \end{pmatrix}.

固定横向频率 $n \neq 0$ 后，算子限制在 $$x_2$$ 的第 $$n$$ 个傅里叶块上。作者观察到 $L_{k,n}^2$ 的对角块分别为 $$P Q$$ 和 $$Q P$$ ，其中 $P = C_k A_{k,n}$ 和 $Q = C_k B_{k,n}$ 。利用算子乘积谱的性质 $\sigma(ST) \setminus \{0\} = \sigma(TS) \setminus \{0\}$ ，将 $L_{k,n}$ 的非零谱问题转化为紧自伴算子 $K_{k,n} = D_{k,n}^{1/2} A_{k,n} D_{k,n}^{1/2}$ 的谱问题，其中 $D_{k,n} = C_k B_{k,n} C_k$ 。由于 $B_{k,n}$ 是紧自伴正算子， $K_{k,n}$ 也是紧自伴的，其特征值序列趋于零。负特征值的个数 $N_{k,n}$ 可以通过计算 $A_{k,n}$ 的负子空间维度直接得到，即满足 $$m^2 + n^2 < k^2$$ 的整数 $$m$$ 的个数。这不仅分类了所有不稳定谱，还排除了右半平面上其他可能的谱点。通过对称性选择特定的傅里叶链，结合标量连分式估计，作者证明了在主导特征值下方存在正的谱间隙，这是后续非线性分析的关键。

3. 高阶近似轨迹与非线性不稳定性证明 从线性不稳定性过渡到非线性不稳定性面临的主要障碍是二次输运项导致的导数损失（derivative loss），这使得直接在固定正则性空间中进行扰动论证变得困难。作者采用了 Grenier 引入的高阶近似轨迹方法。具体步骤如下：

半群界: 利用线性化算子在 Sobolev 空间 $$X_s$$ 上生成的均匀连续半群，获得具有相同指数增长率的先验估计。
驯服能量估计 (Tame Energy Estimates): 控制非线性项的增长，防止解在有限时间内爆破。
递归构造近似解: 构造一个近似解，其残差（residual）是高阶小的。通过迭代修正，可以使得残差达到任意高的阶数。
Bootstrap 论证: 利用谱间隙提供的分离度，证明精确解在长时间尺度 $T_\epsilon = \Lambda_k^{-1} \log(\delta/\epsilon)$ 内保持在近似解附近，直到主导模式增长到固定大小 $\delta$ 。这种方法克服了导数损失问题，成功建立了对数时间尺度上的非线性逃逸。

4. 与之前工作的比较

本文将二维惯性 MHD 系统的不稳定性分析与现有的理想流体和磁流体力学文献进行了细致的比较，指出了其独特的贡献和局限。

1. 与 Constantin 和 Hu [11] 工作的比较 Constantin 和 Hu 研究了一种相关的平面惯性模型，但其 $$F$$ 方程左侧多了一项 $\nabla^\perp \Phi \cdot \nabla \omega$ 。在他们的模型中，两个场 $F \pm \omega$ 分别由流函数 $\psi \pm \Phi$ 输运，这种特殊的输运结构导致了全局正则性，并允许无电阻情况下的磁拓扑变化。相比之下，本文研究的系统 (1.1) 缺乏这种特定的输运结构，因此不能直接应用他们的正则性结果。相反，本文关注的是该系统固有的不稳定性，这表明不同的耦合结构会导致截然不同的动力学行为（一个是正则且拓扑变化，另一个是不稳定）。

2. 与 Liao, Lin, 和 Zhu [10] 工作的比较 Liao, Lin, 和 Zhu 研究了理想 MHD 中 Kelvin-Stuart 磁岛的稳定性。他们证明了共周期扰动下的轨道稳定性，以及双倍周期扰动下的线性不稳定性（合并不稳定性）。本文的纯磁正弦波族可以看作是某种简化的磁岛结构。然而，本文的系统包含了电子惯性项 $(1-\Delta)\Phi = F$ ，这与理想 MHD 不同。更重要的是，本文不仅证明了线性不稳定性，还进一步证明了非线性不稳定性，并且给出了明确的对数时间尺度 $T_\epsilon \sim \log(1/\epsilon)$ 。Liao 等人的工作主要集中在理想流体框架下的轨道稳定性，而本文则深入到了惯性效应引起的非线性增长机制。

3. 与物理文献 [1-6, 8, 9] 的比较 物理学家如 Ottaviani 和 Porcelli [3] 早已在数值模拟中观察到撕裂模的加速增长和电流层变窄，并将其用于解释磁重联。然而，这些通常是启发式或数值的结果。本文提供了严格的数学证明，确立了谱不稳定性（定理 2.1, 3.11）和非线性不稳定性（定理 4.2）。与物理直觉相比，数学证明揭示了不稳定性发生的精确参数范围（如 $$a^2(1+k^2) < 2$$ ）和时间尺度。但是，作者也诚实地指出，目前的非线性定理仅证明了 Sobolev 范数下的逃逸，并未直接证明物理上关注的磁拓扑变化或电流片形成。因此，虽然数学上不稳定性是确定的，但从数学不稳定性直接推导到物理重联过程之间仍存在差距，这是一个开放问题。

4. 方法论上的推进 在线性不稳定性方面，传统的 Euler 线性不稳定性理论 [13-16] 通常处理标量递推关系。本文处理的是矩阵值递推关系，这需要发展新的矩阵连分式技术来处理耦合效应。在非线性方面，Grenier [19] 的方法主要用于理想流体，本文将其推广到带有导数损失的惯性 MHD 系统，并结合了紧算子谱理论提供的谱间隙信息，这是对现有非线性不稳定性工具箱的重要补充。总体而言，本文在假设上放宽了对特殊输运结构的依赖，在结论上从线性谱分析推进到了非线性时间演化，并在适用范围上覆盖了包含电子惯性的更一般模型。

8. 临界空间中三维各向异性 Navier-Stokes 方程解的稳定性

On the stability of solutions to the 3-D anisotropic Navier-Stokes equations in the critical space

作者: Nicolas Burq, Ning Liu, Ping Zhang

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14304v1

摘要翻译: 本文研究三维不可压缩各向异性 Navier-Stokes 方程 $$(ANS)$$ 在临界空间 $$ rac12}$$ 中初值下解的连续依赖性。我们证明了在此设定下，从数据到解的映射是连续的。由于问题的临界正则性，经典的 Bona-Smith 论证 rac12} $$（另见抽象结果 rac12}$$ ）似乎不适用，因为一个关键要素（较低正则性处的 Lipschitz 连续性）似乎不成立。因此，我们在 Besov 框架内发展了一种新的直接 Fourier 加权方法。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩各向异性 Navier-Stokes 方程 (ANS) 在临界空间中的解对初值的连续依赖性。具体而言，考虑如下方程组：

\begin{cases} \partial_t u + u \cdot \nabla u - \Delta_h u = -\nabla p, & (t, x) \in \mathbb{R}^+ \times \mathbb{R}^3, \\ \text{div } u = 0, \\ u|_{t=0} = u_0, \end{cases}

其中 $\Delta_h \stackrel{\text{def}}{=} \partial_1^2 + \partial_2^2$ 是水平方向的拉普拉斯算子， $$u$$ 为流体速度场， $$p$$ 为满足无散度条件的标量压力函数。

该方程组的物理背景源于大气或海洋流体力学中垂直方向扩散远小于水平方向扩散的情形。与经典的各向同性 Navier-Stokes 方程不同，(ANS) 仅在水平方向具有耗散项 $-\Delta_h u$ ，而在垂直方向没有直接的粘性耗散。这种结构导致了正则性理论上的显著困难：由于缺乏垂直方向的耗散，建立唯一性和稳定性估计时需要处理额外的导数损失问题。

从尺度不变性的角度来看，若 $$u(t,x)$$ 是 (ANS) 在 $$[0,T]$$ 上以 $$u_0$$ 为初值的解，则对于任意 $\lambda > 0$ ，定义

u_\lambda(t,x) \stackrel{\text{def}}{=} \lambda u(\lambda^2 t, \lambda x), \quad u_{0,\lambda}(x) \stackrel{\text{def}}{=} \lambda u_0(\lambda x),

则 $u_\lambda$ 也是 (ANS) 在 $[0, T/\lambda^2]$ 上以 $u_{0,\lambda}$ 为初值的解。这一尺度变换确定了局部适定性的最小正则性要求。一个函数空间被称为“临界的”，如果其范数在该尺度变换下保持不变。对于 (ANS)，临界空间包括各向异性 Sobolev 空间 $H^{0,s}$ （其中 $$s=1/2$$ ）以及各向异性 Besov 空间 $\dot{B}^{0, 1/2}$ 。

在已有文献中，Chemin 等人 [5] 证明了当初值 $u_0 \in H^{0,s}$ 且 $s \in (1/2, 1)$ 时，(ANS) 存在局部解；若初值足够小，则解全局存在。然而，由于垂直方向缺乏耗散，他们在证明唯一性时需要比存在性多一个垂直方向的导数。Iftimie [9] 后来解决了 $$s > 1/2$$ 时的唯一性问题。Paicu [11] 进一步证明了初值在临界空间 $\dot{B}^{0, 1/2}$ 中的局部存在性，并在初值足够小时证明了全局存在性。Chemin 和 Zhang [7] 引入了负指标临界空间 $\dot{B}^{-1/2, 1/2}_4$ 并建立了相应的小初值全局存在性结果。值得注意的是，在这些工作中，存在性和唯一性通常在不同的函数空间中建立，这再次反映了垂直方向缺乏耗散带来的困难。

Hadamard 意义下的适定性要求：(1) 解的存在性与唯一性；(2) 解对初值的连续依赖性。对于经典 Navier-Stokes 方程，解映射关于初值是 Lipschitz 连续的 [8]。Zhang 和 Zhu [15] 证明了 (ANS) 的解在 $H^{0,s}$ ( $$s > 1/2$$ ) 中对初值连续依赖，但解映射在 $H^{0,s}$ 的有界子集上不是均匀连续的。在此之前，(ANS) 在临界空间 $\dot{B}^{0, 1/2}$ 中的连续性一直是一个开放问题。

本文的目标正是填补这一空白，即建立 (ANS) 解映射在临界空间 $\dot{B}^{0, 1/2}$ 中的连续性。这一结果不仅完善了 (ANS) 的适定性理论，也为理解临界尺度下非线性偏微分方程的稳定性提供了新的范例。

2. 主要定理或主要结果

本文的主要结果是定理 1.1，它确立了 (ANS) 解映射在临界空间 $\dot{B}^{0, 1/2}$ 中的连续性。

定理 1.1. 设 $u_\infty$ 是 (ANS) 在 $$[0, T]$$ 上以初值 $u_{0,\infty} \in \dot{B}^{0, 1/2}$ 的强解，且设 $\{u_{0,n}\}_{n \in \mathbb{N}}$ 是 $\dot{B}^{0, 1/2}$ 中的一列函数。若 $u_{0,n}$ 在 $\dot{B}^{0, 1/2}$ 中收敛于 $u_{0,\infty}$ （当 $n \to \infty$ ），则存在 $N \in \mathbb{N}$ ，使得对所有 $$n > N$$ ，存在 (ANS) 在 $$[0, T]$$ 上以 $u_{0,n}$ 为初值的唯一解 $$u_n$$ ，并且

\lim_{n \to \infty} \left( \|u_n - u_\infty\|_{L_T^\infty(\dot{B}^{0, 1/2})} + \|\nabla_h(u_n - u_\infty)\|_{L_T^2(\dot{B}^{0, 1/2})} \right) = 0.

为了证明该定理，作者发展了一套基于 Fourier 加权 Besov 空间的新方法。主要技术引理包括：

命题 2.1 (初值序列的改进收敛性). 若 $\{u_{0,n}\}_{n \in \mathbb{N}} \subset \dot{B}^{0, 1/2}$ 在 $\dot{B}^{0, 1/2}$ 中收敛于 $u_{0,\infty}$ ，则存在一个非递减的正算术函数 $b = \{b_j\}_{j \in \mathbb{Z}}$ 满足 $\lim_{j \to \infty} b_j = +\infty$ ，使得对所有 $n \in \mathbb{N}$ ， $u_{0,n}$ 属于加权空间 $\dot{B}^{0, 1/2}_b$ ，且 $\{u_{0,n}\}_{n \in \mathbb{N}}$ 在 $\dot{B}^{0, 1/2}_b$ 中收敛于 $u_{0,\infty}$ 。

这里，加权 Besov 空间 $\dot{B}^{0, 1/2}_b$ 定义为满足以下范数有限的 tempered 分布 $$f$$ ：

\|f\|_{\dot{B}^{0, 1/2}_b} \stackrel{\text{def}}{=} \sum_{j \in \mathbb{Z}} b_j 2^{j/2} \|\dot{\Delta}_v^j f\|_{L^2} < \infty.

命题 2.2 ( $\dot{B}^{0, 1/2}_b$ 中的 tame 估计). 设 $b = \{b_j\}_{j \in \mathbb{Z}}$ 是非递减正算术函数，满足 $b_j \leq b_{j+1} \leq 2b_j$ 。若 $$u(t)$$ 是 (ANS) 在 $$[0, T]$$ 上以 $u_0 \in \dot{B}^{0, 1/2}_b$ 为初值的解，则对所有 $t \in [0, T]$ ，有

\|u\|_{\tilde{L}_t^\infty(\dot{B}^{0, 1/2}_b)} + \|\nabla_h u\|_{\tilde{L}_t^2(\dot{B}^{0, 1/2}_b)} \leq C e^{C A(t)} \|u_0\|_{\dot{B}^{0, 1/2}_b},

其中

A(t) \stackrel{\text{def}}{=} \int_0^t \left(1 + \|u(\tau)\|_{\dot{B}^{0, 1/2}}^2 \right) \left(1 + \|\nabla_h u(\tau)\|_{\dot{B}^{0, 1/2}}^2 \right) d\tau.

这里的 $\tilde{L}_T^p$ 表示 Chemin-Lerner 型范数。

命题 2.3 ( $$X$$ 范数中的弱连续性). 设 $$u(t)$$ 和 $$v(t)$$ 是 (ANS) 在 $$[0, T]$$ 上分别以 $$u_0$$ 和 $$v_0$$ 为初值的两个解。若 $\|u - v\|_{\tilde{L}_T^\infty(X)} \leq 2^{-10}$ ，则对所有 $t \in [0, T]$ ，有

\|g(u-v)\|_{\tilde{L}_t^\infty(X)}^2 + \|g \nabla_h(u-v)\|_{\tilde{L}_t^2(X)}^2 \leq C \|u_0 - v_0\|_X^2 + C \int_0^t B'(\tau) \frac{-\ln \|(g(u-v))\|_{\tilde{L}_\tau^\infty(X)}^2 \cdot \ln\left(-\ln \|(g(u-v))\|_{\tilde{L}_\tau^\infty(X)}^2 \right)}{\|(g(u-v))\|_{\tilde{L}_\tau^\infty(X)}^2} d\tau,

其中

g(t) \stackrel{\text{def}}{=} e^{-C B(t)}, \quad B(t) \stackrel{\text{def}}{=} \int_0^t \left(1 + \|(u,v)(\tau)\|_{\dot{B}^{0, 1/2}}^2 \right) \left(1 + \|\nabla_h(u,v)(\tau)\|_{\dot{B}^{0, 1/2}}^2 \right) d\tau.

范数 $\|\cdot\|_X$ 定义为低频部分用 $\dot{B}^{0, 1/2}$ 控制，高频部分用 $H^{0, -1/2}$ 控制的混合范数：

\|f\|_X \stackrel{\text{def}}{=} \|f\|_{\dot{B}^{0, 1/2}_{lo}} + \|f\|_{\dot{H}^{0, -1/2}_{hi}},

其中

\|f\|_{\dot{B}^{0, 1/2}_{lo}} = \sum_{j \leq -1} 2^{j/2} \|\dot{\Delta}_v^j f\|_{L^2}, \quad \|f\|_{\dot{H}^{0, -1/2}_{hi}} = \left( \sum_{j \geq 0} 2^{-j} \|\dot{\Delta}_v^j f\|_{L^2}^2 \right)^{1/2}.

通过结合上述命题，作者首先将差值 $u_n - u_\infty$ 分解为高频部分 $(Id - \dot{S}_v^J)(u_n - u_\infty)$ 和低频部分 $\dot{S}_v^J(u_n - u_\infty)$ 。高频部分利用命题 2.2 和加权空间的性质进行控制，而低频部分利用命题 2.3 提供的弱连续性估计。最终，通过 Osgood 引理和标准连续性论证，得出定理 1.1 的结论。

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

本文的核心创新在于开发了一种新的直接 Fourier 加权方法，以克服在临界空间 $\dot{B}^{0, 1/2}$ 中应用经典 Bona-Smith 论证的障碍。Bona-Smith 论证（及其抽象版本 [1]）通常依赖于在较低正则性空间中的 Lipschitz 连续性以及在较高正则性空间中的 tame 估计，并通过非线性插值来推导中间空间的连续性。然而，在 (ANS) 的临界情形下，由于垂直方向缺乏耗散，在较弱范数（如 $\dot{B}^{0, -1/2}$ ）中难以获得所需的 Lipschitz 估计。已有的唯一性证明 [11] 仅给出类似 $N(t) \leq e^{\log(N(t_0))e^{-C|t-t_0|}}$ 的估计，这虽然足以证明唯一性（因为 $e^{-\infty}=0$ ），但远弱于 Lipschitz 估计 $N(t) \leq C_t N(t_0)$ ，因此无法直接应用于 Bona-Smith 框架。

作者的策略是绕过抽象插值，转而采用一种“直接初等插值”的方法。具体步骤如下：

构造改进的初值收敛性（命题 2.1）： 关键观察是，如果一个序列在 $\dot{B}^{0, 1/2}$ 中收敛，那么实际上可以在一个略平滑的空间 $\dot{B}^{0, 1/2}_b$ 中实现收敛，其中权重 $$b_j$$ 随频率 $$j$$ 缓慢增长至无穷大。作者通过 Lemma 3.1 和 Lemma 3.2 构造了这样的算术函数 $b = \{b_j\}$ 。Lemma 3.1 表明任何 $\ell^1(\mathbb{N})$ 序列都属于某个加权 $\ell^1(\mathbb{N})$ 空间，且权重可以趋于无穷大。Lemma 3.2 将此推广到一列趋于零的 $\ell^1$ 序列，确保整个序列在加权空间中一致有界且收敛。这一构造利用了 $\dot{B}^{0, 1/2}$ 范数的绝对求和性质，允许通过适当选择权重来吸收高频部分的微小误差。
传播加权正则性（命题 2.2）： 接下来，作者证明了加权正则性 $\dot{B}^{0, 1/2}_b$ 在时间演化下是稳定的。这需要建立 tame 估计，其中常数依赖于基本临界范数 $\dot{B}^{0, 1/2}$ 的时间积分。由于权重 $$b_j$$ 的增长是缓慢的（满足 $b_j \leq b_{j+1} \leq 2b_j$ ），非线性项在处理高频部分时不会导致导数损失，从而使得加权范数能够被控制。这一步骤确保了由 $u_{0,n}$ 生成的解 $$u_n$$ 在 $\dot{B}^{0, 1/2}_b$ 中一致有界。
弱连续性估计（命题 2.3）： 对于低频部分，作者引入混合范数 $$X$$ ，其中低频使用 $\dot{B}^{0, 1/2}$ ，高频使用 $H^{0, -1/2}$ 。由于 $H^{0, -1/2}$ 可以克服高频部分的垂直导数损失，而低频部分本身不需要额外正则性，这种混合范数是自然的。作者建立了带有时间权重的 Chemin-Lerner 范数下的连续性估计。该估计的形式类似于 Osgood 型不等式，包含对数修正项。通过应用 Osgood 引理，可以得到差值在 $$X$$ 范数中的衰减估计，尽管这不是 Lipschitz 型的，但对于低频部分已经足够，因为低频部分的有限维性使得弱收敛可以强化为强收敛。
高低频分解与组合： 最后，将差值 $u_n - u_\infty$ 分解为高频和低频两部分。高频部分利用加权空间的性质：由于 $b_J \to \infty$ ，可以选择足够大的 $$J$$ 使得高频尾部在 $\dot{B}^{0, 1/2}$ 中任意小，且这一估计独立于 $$n$$ 。低频部分利用命题 2.3 的弱连续性，并结合 $\dot{B}^{0, 1/2}$ 收敛性，得到低频差值在 $\dot{B}^{0, 1/2}$ 中趋于零。两部分结合，完成了定理 1.1 的证明。

这种方法的优势在于不依赖于全局的 Lipschitz 估计，而是巧妙地利用了临界空间收敛所隐含的“额外平滑性”（通过加权空间体现），并将问题分解为可分别处理的高频和低频部分。作者指出，这种方法具有普遍性，可能适用于其他 Lipschitz 收缩估计不成立的情形。

4. 与之前工作的比较

本文将 (ANS) 的适定性理论推进到了临界空间 $\dot{B}^{0, 1/2}$ 的连续性层面，填补了此前文献中的空白。以下从几个维度比较本文与既有工作的区别：

与 Bona-Smith 论证的比较： 经典的 Bona-Smith 论证 [3] 及其抽象形式 [1] 是证明非线性 PDE 解映射连续性的标准工具。它依赖于三个要素：(i) 在较光滑空间 $X_{s_2}$ 中的 tame 估计；(ii) 在较粗糙空间 $X_{s_1}$ 中的 Lipschitz 估计；(iii) 目标空间 $$X_s$$ 位于两者之间（ $$s_1 < s < s_2$$ ）。通过非线性插值，可以从这些估计推导出 $$X_s$$ 中的连续性。对于 (ANS)，Zhang 和 Zhu [15] 在 $H^{0,s}$ ( $$s > 1/2$$ ) 中成功应用了这一策略，因为他们可以在 $H^{0, 1/2}$ （对应 $$s_1$$ ）中获得足够的 Lipschitz 控制，并在更光滑的空间中获得 tame 估计。然而，在临界空间 $\dot{B}^{0, 1/2}$ 中，这一策略失效。原因是，为了避免垂直导数损失，自然选择的较弱空间是 $\dot{B}^{0, -1/2}$ 或 $H^{0, -1/2}$ 。但是，Paicu [11] 建立的唯一性估计 (2.9) 显示，差值在 $H^{0, -1/2}$ 中的演化满足一个带有对数修正的 Osgood 型不等式，而非线性的 Lipschitz 不等式。具体来说，估计形式为 $d/dt \|w\|^2 \leq C B' \|w\|^2 (1 - \ln \|w\|^2) \ln(1 - \ln \|w\|^2)$ ，这导致解的行为类似于 $e^{-e^{-t}}$ ，而不是指数衰减 $e^{-Ct}$ 。这种亚 Lipschitz 行为意味着 Bona-Smith 框架中所需的 Lipschitz 条件不成立。因此，作者必须发展一种替代方法，这正是本文的核心贡献。
与 Paicu [11] 和 Chemin-Zhang [7] 的比较： Paicu [11] 证明了 (ANS) 在 $\dot{B}^{0, 1/2}$ 中的局部存在性，并在小初值情况下证明了全局存在性。他的唯一性证明依赖于 $H^{0, -1/2}$ 估计，这表明存在性和唯一性是在不同范数下建立的。类似地，Chemin 和 Zhang [7] 在负指标空间 $\dot{B}^{-1/2, 1/2}_4$ 中建立了小初值的全局存在性，同样面临存在性与唯一性空间不一致的问题。这些工作均未涉及解映射的连续性。本文不仅确认了解在 $\dot{B}^{0, 1/2}$ 中的存在性和唯一性（作为连续性的前提），还进一步证明了连续依赖性，从而在 Hadamard 意义下完成了该临界空间中的适定性理论。
与 Zhang-Zhu [15] 的比较： Zhang 和 Zhu [15] 证明了 (ANS) 在超临界空间 $H^{0,s}$ ( $$s > 1/2$$ ) 中的连续性。他们的证明依赖于 Bona-Smith 论证，因此需要 $$s > 1/2$$ 以确保有足够的正则性余量来获得 Lipschitz 估计。此外，他们指出解映射在 $H^{0,s}$ 的有界子集上不是均匀连续的。本文的结果将连续性范围缩小到临界指数 $$s=1/2$$ ，这是理论上最小的正则性要求。虽然本文未明确讨论均匀连续性，但通过加权空间的方法，作者展示了即使在临界尺度下，连续性依然成立，只是证明技术完全不同。
方法论上的推进： 本文的方法论创新在于引入了 Fourier 加权 Besov 空间 $\dot{B}^{0, 1/2}_b$ 作为桥梁。这种加权空间允许作者在保持临界范数不变的同时，提取出初始数据收敛所隐含的额外高频信息。通过证明加权正则性的传播（命题 2.2）和混合范数下的弱连续性（命题 2.3），作者避免了全局 Lipschitz 估计的需求。这种方法不仅适用于 (ANS)，而且如作者所述，可能推广到其他具有类似亚 Lipschitz 行为的非线性 PDE。相比之下，传统方法往往依赖于全局的 Lipschitz 或 Hölder 连续性，这在临界或亚临界情形下经常失败。
局限性与代价： 本文方法的代价是需要构造特定的算术函数 $$b_j$$ ，这使得证明变得较为技术性。此外，连续性估计的最终形式 (1.4) 是逐点极限，而非显式的模连续性（modulus of continuity）。这是因为弱连续性估计 (2.8) 中的对数项使得很难提取出一个统一的、关于初值距离的连续函数。尽管如此，这仍然满足了 Hadamard 适定性的要求。未来工作或许可以探索是否能获得更精确的连续性模，或者将此方法应用于其他各向异性系统。

9. 复合状态方程的可压缩欧拉方程的波速保证界

Guaranteed wave-speed bounds for the compressible Euler equations with a composite equation of state

作者: Nicolas Favrie, Matthias Maier

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

发布日期: 2026-09-13

arXiv ID: 2609.14287v2

摘要翻译: 我们推导了具有复合状态方程的可压缩欧拉方程黎曼问题中最大波速的可计算上界，其中压力由一个凸的水动力学部分和一个正压修正项之和构成。该修正可能破坏状态方程的凸性，从而使黎曼解可包含复合激波-稀疏波。这些界通过仅针对水动力学状态方程的辅助黎曼问题进行比较而获得。无论外部波是基本波还是复合波，这些界均成立，并且无需对波的类型进行情形区分。

Comment: 无

链接: http://arxiv.org/abs/2609.14287v2

1. 研究对象和背景

本文研究一维可压缩欧拉方程在复合状态方程下的黎曼问题，旨在推导最大波速的可计算且严格的上界。研究对象为如下守恒律系统：

\partial_t U + \partial_x F(U) = 0,

其中 $U = [\rho, u, E]^T$ ，通量 $$F(U)$$ 对应于质量、动量和能量守恒。具体地，方程组写为：

\begin{cases} \partial_t \rho + \partial_x (\rho u) = 0, \\ \partial_t (\rho u) + \partial_x (\rho u^2 + p) = 0, \\ \partial_t E + \partial_x (u(E+p)) = 0. \end{cases}

初始条件为分段常数函数：

U(x,0) = \begin{cases} U_L, & x < 0, \\ U_R, & x > 0. \end{cases}

该问题的核心在于压力 $$p$$ 和总能量 $$E$$ 由一个“复合”状态方程给出：

p = p_h(\rho, e_h) + p_b(\rho),

E = \frac{1}{2}\rho u^2 + e_h + e_b(\rho),

其中 $p_h(\rho, e_h)$ 是满足经典热力学凸性条件的流体动力学部分（如理想气体或刚性气体），而 $p_b(\rho)$ 是仅依赖于密度的等熵修正项（barotropic correction）。两者通过关系 $p_b(\rho) = \rho^2 \partial_\rho e_b(\rho)$ 关联。

这类模型出现在 Mie-Grüneisen 型闭合关系中（结合冷压与热压）以及超弹性固体的分离能量模型中。关键困难在于，虽然 $$p_h$$ 保证了辅助问题的经典三波结构（激波/稀疏波-接触间断-激波/稀疏波），但 $$p_b$$ 的存在可能破坏总状态方程沿等熵线的凸性。当等熵压力关于比容的曲率改变符号时，声学特征场可能失去真正的非线性，导致黎曼解中出现复合波（composite shock-rarefaction waves），即同一特征族内包含连续的激波和稀疏波段。

在数值方法中，特别是 Godunov 型方法（如 HLL 格式、局部 Lax-Friedrichs 格式或图粘性不变域保持格式），需要估计局部最大传播速度 $\lambda_{\max} = \max(|\lambda_1^-|, |\lambda_3^+|)$ 以确定时间步长（CFL 条件）或人工粘性系数。如果低估了波速，可能导致数值耗散不足，破坏稳定性、正性或不变域性质；如果高估过多，则增加数值扩散并限制时间步长。因此，目标是在不构造完整且复杂的非线性波扇（wave fan）的情况下，获得一个既廉价又严格的波速上界。

现有文献中，对于理想气体或简单状态方程，已有基于迭代或解析公式的估计（如 Guermond-Popov, Davis）。然而，对于一般或非凸状态方程，精确求解黎曼问题需要解决非线性标量方程、进行数值积分甚至识别凸性变化和过渡态，计算成本极高且不具系统性。现有的近似求解器通常只依赖初始状态的线性化估计，无法保证对非线性复合波的严格界限。此外，基于代理模型（surrogate models）的方法虽然能保证代理模型的界限，但并不直接约束原始状态方程的真实波速。因此，本文试图填补的空白是：针对包含非凸等熵修正的复合状态方程，建立一种无需显式构造完整波扇、适用于基本波和复合波、且能严格界定极端波速的分析框架。

2. 主要定理或主要结果

本文的主要结果集中在两个定理及其推论，提供了针对不同假设条件下的波速上界。

定理 2.1 (通用波速界限): 假设左右外侧波均为压缩波或膨胀波。定义辅助流体动力学问题（仅含 $$p_h$$ ）的冲击函数 $$f_h^S$$ 。令 $$q_L, q_R$$ 和 $\underline{q}_L, \underline{q}_R$ 分别为以下方程的唯一解：

f_h^S(q_L; W_L) = \max(u_L - u_R, 0) = f_h^S(q_R; W_R),

f_h^S(\underline{q}_L; W_L) = \min(u_L - u_R, 0) = f_h^S(\underline{q}_R; W_R).

进而定义上下星压力边界：

\overline{p}_{h,L}^* := \max(q_L, p_R - \hat{p}_b(v_L)), \quad \overline{p}_{h,R}^* := \max(q_R, p_L - \hat{p}_b(v_R)),

\underline{p}_{h,L}^* := \min(q_L, p_R - \hat{p}_b(v_L)), \quad \underline{p}_{h,R}^* := \min(q_R, p_L - \hat{p}_b(v_R)).

通过 Rankine-Hugoniot 条件确定对应的比容 $\overline{v}_{L}^*, v_{L}^*$ 等。则完整问题 (2.1) 的极端波速 $\lambda_1^-$ 和 $\lambda_3^+$ 满足：

\lambda_1^- \ge u_L - \sqrt{(u_L - \overline{\lambda}_{1,h}^-)^2 + 2(v_L)^2 \max_{\overline{v}_{L}^* \le \tau \le v_{L}^*} \left( \frac{a_b^2(\tau)}{\tau^2} \right)},

\lambda_3^+ \le u_R + \sqrt{(ar{\lambda}_{3,h}^+ - u_R)^2 + 2(v_R)^2 \max_{\overline{v}_{R}^* \le \tau \le v_{R}^*} \left( \frac{a_b^2(\tau)}{\tau^2} \right)}.

此界限的关键特点是：无论外侧波是基本波还是复合波，也无论其是压缩还是膨胀，均成立，且无需区分波的类型。

定理 2.3 (正则波速界限): 在定理 2.1 的基础上，若满足特定的小性条件和凸性准则（确保外侧波为正则激波或正则稀疏波，即基本导数保持为正），则可以通过单调迭代序列 $\overline{P}_{L/R}^{(n)}$ 和 $\underline{P}_{L/R}^{(n)}$ 来锐化星压力的界限。对于任意迭代步 $n \ge 0$ ，上述波速界限依然有效，但使用了更紧的压力区间 $[\underline{P}_{L/R}^{(n)}, \overline{P}_{L/R}^{(n)}]$ 。每个迭代步骤都产生一个有效的波速上界。

推论 2.5 (足够小区间的波速界限): 如果在某次迭代后，压力区间 $[\underline{P}_{L}^{(n)}, \overline{P}_{L}^{(n)}]$ 不包含相邻初始状态的压力 $p_{h,L}$ （即区间足够小，从而确定了波的类型），则可以去掉定理 2.1 中因不确定波类型而产生的因子 2 带来的保守性，获得更精确的界限。例如，若确定为左侧激波，则：

\lambda_1^- \ge u_L - \sqrt{-1 + \frac{\hat{p}_b(\overline{v}_{L}^*) - \hat{p}_b(v_L)}{\underline{P}_{L}^{(n)} - p_{h,L}}}(u_L - \overline{\lambda}_{1,h}^-)^2 + (v_L)^2 \max_{\overline{v}_{L}^* \le \tau \le v_L} \left( \frac{a_b^2(\tau)}{\tau^2} \right).

若确定为左侧稀疏波，则界限简化为：

\lambda_1^- = u_L - \sqrt{a_h^2(\rho_L, e_h(\rho_L, p_{h,L})) + a_b^2(\rho_L)}.

这些结果表明，通过简单的迭代过程，可以在保证严格性的同时显著提高波速估计的精度。

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

本文的核心技术障碍在于如何处理非凸状态方程导致的复合波结构，同时避免昂贵的精确黎曼求解器。作者的真正创新点在于利用压力加性分解结构，将完整系统的波曲线参数化并与辅助流体动力学系统进行比较。

参数化策略与辅助问题: 作者没有直接求解完整系统的非线性匹配方程，而是引入仅由 $$p_h$$ 控制的辅助 Riemann 问题 (2.8)。由于 $$p_h$$ 满足凸性条件，辅助问题具有经典的三波结构，其星压力 $\tilde{p}_h^*$ 可通过单一标量方程 (2.10) 唯一确定。完整系统的波曲线被参数化为流体动力学星压力的函数。等熵修正项 $$p_b$$ 的影响仅通过可计算的修正项进入最终界限，从而完全避免了构建完整波扇。
Hugoniot 曲线与冲击函数的比较原理: 对于激波分支，作者建立了完整系统与辅助系统之间 Hugoniot 曲线的排序原则。通过分析 Rankine-Hugoniot 条件，证明了在特定条件下，完整系统的冲击速度可以被辅助系统的相应量加上一个由 $$a_b^2$$ 控制余项所界定。这允许使用辅助系统的冲击函数 $$f_h^S$$ 的延拓来界定完整系统的行为。
稀疏波的速度界限: 对于稀疏波，作者利用了声速的合成关系 $$a^2 = a_h^2 + a_b^2$$ 。通过积分声速沿等熵线的变化，证明了完整系统稀疏波边缘的速度可以被辅助系统的稀疏波速度加上一个涉及 $$a_b$$ 最大值的项所界定。这一技巧将非凸部分的复杂性转化为对已知函数 $$a_b$$ 的最大值计算。
复合波的处理: 这是最关键的难点。对于复合波（由交替的激波和稀疏波组成），传统的单值压力函数不再适用。作者证明，尽管波结构复杂，但其极端速度仍然受限于由辅助系统确定的某些“包络”状态。具体来说，他们证明了无论复合波内部有多少个基本段，其最外侧的速度都可以被归结为对辅助系统中某些极值压力（由 $$q_L, q_R$$ 等定义）的评估。这使得界限对所有类型的波统一成立，无需检测凸性变化或定位过渡态。
迭代锐化机制: 在正则波假设下，作者设计了一个单调收敛的迭代方案。利用 Lemma 4.6 和 4.7 提供的体积偏差界限和压力偏移界限，递归地更新星压力的上下界。每一步迭代都缩小了可能的星压力区间，从而使得最终的波速界限更加紧致。这种迭代不需要知道确切的波类型，直到区间足够小以排除另一种可能性为止。

整个分析框架依赖于严格的热力学稳定性假设（A.0-A.4）以及对 $$p_b$$ 的平滑性和有界性假设。通过这些结构性假设，作者成功地将一个全局非线性问题分解为一系列局部的、可计算的代数不等式和最大值运算。

4. 与之前工作的比较

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

与非凸状态方程精确求解器的对比: 传统方法（如 Sirianni, Guardone, Re [27] 或 Marquina et al. [20]）试图通过自动检测波曲线和解决非线性方程来构造精确或半精确的 Riemann 解。这些方法虽然能提供高精度的参考解，但计算成本高，且在多维计算中难以系统地应用于每个界面。本文提出的界限不需要构造波扇，也不需要解决任何非线性方程，仅需评估辅助系统的冲击函数和 $$p_b$$ 的导数，因此在计算效率上具有数量级的优势，适合用于实时数值模拟中的 CFL 条件计算。
与线性化或基于初始状态的估计的对比: Davis [6] 等经典公式基于初始状态的谱半径或线性化近似。Toro, Müller, Siviglia [32] 和 Toro, Tokareva [31] 指出，这些公式在某些情况下会低估最大波速，导致单调性丧失或稳定性降低。本文的界限是严格的（guaranteed upper bounds），明确考虑了非线性效应和 $$p_b$$ 的贡献，因此不会低估波速。代价是可能会比精确值稍大（overestimate），但这在数值稳定性上是可接受的，且通过迭代可以减小过估计的程度。
与代理模型方法的对比: 近期工作（如 [3, 4, 11]）使用局部热力学代理模型来保证不变域性质。这些方法保证的是代理模型的波速界限，而非原始状态方程的真实波速。虽然它们确保了物理量的合理性，但不能直接作为原始方程 Riemann 问题的信号速度界限。本文的方法直接针对原始复合状态方程，其界限严格适用于原始方程的黎曼解，因此更直接地服务于 HLL 型通量或图粘性系数的选择。
适用范围与假设差异: 本文假设 $$p_h$$ 满足标准凸性，而 $$p_b$$ 仅要求保持双曲性（即声速实数），允许其破坏总凸性。这与许多假设整体凸性的早期工作不同。本文的结果特别适用于 Mie-Grüneisen 等实际物理模型，其中冷压项 $$p_b$$ 往往是非凸的。相比之下，Guermond and Popov [10] 的工作主要针对理想气体或协体积状态方程，其凸性较强，未涵盖此类复合非凸情形。
局限性与代价: 主要局限在于，为了获得统一的界限，定理 2.1 中的公式包含一个因子 2（源于对波类型的不确定性），这可能导致一定的保守性。然而，定理 2.3 和推论 2.5 通过迭代消除了这一保守性，前提是波是正则的。如果波高度非正则或处于临界状态，迭代可能收敛较慢或需要更多步骤。此外，该方法依赖于 $$p_b$$ 的解析表达或易于求导的性质，对于纯表格化的状态方程，可能需要额外的插值处理来计算 $$a_b^2$$ 的最大值。尽管如此，相比于精确求解器的不可行性，这种计算代价是完全可接受的。

10. 三维可压缩非耗散磁流体力学的整体稳定性：隐藏耗散与有理背景场

Global Stability of 3D Compressible Non-Resistive MHD: Hidden Damping and Rational Background Fields

作者: Lin-An Li, Qian Li, Jiahong Wu, Xiaojing Xu

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14244v1

摘要翻译: 我们证明了在 $\mathbb{T}^3$ 上，对于小的 $$x_3$$ -对称扰动且对背景磁场不施加任何丢番图条件时，三维可压缩粘性、非耗散 MHD 系统经典解在平衡态 $(1,\mathbf{0},\mathbf{e}_3)$ 附近的整体适定性与非线性稳定性。核心障碍在于 $$x_3$$ -独立扇区，在该扇区中密度和磁场不具备耗散、阻尼或衰减。我们通过揭示精确总压力 $\mathcal{D}=P(1+a)+B_3+\frac12|\mathbf{B}|^2$ 的隐藏波动结构来克服这一障碍，其中平均量对 $(\overline{\mathbf{u}},\overline{\mathcal{D}})$ 满足一个封闭系统，且 $\overline{\mathcal{D}}$ 满足一个强阻尼波动方程，其主波动部分以快磁声速传播，而其非抛物分支的阻尼速率不涉及导数增益。这种隐藏耗散替代了缺失的磁耗散，并同时吸收了由可压缩性引起的主要障碍——磁压力 $\frac12\nabla|\mathbf{B}|^2$ 。结合振荡扇区的阻尼波动结构以及带有移位时间权重的时空加权能量泛函，这给出了整体存在性、一致稳定性和显式的多项式衰减率。

Comment: 无

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

1. 研究对象和背景

本文研究的对象是定义在三维环面 $\mathbb{T}^3$ 上的三维可压缩粘性、无电阻磁流体动力学 (MHD) 方程组。该方程组描述了导电流体在磁场作用下的运动，其数学形式为:

\begin{cases} \partial_t \rho + \text{div}(\rho \mathbf{u}) = 0, \\ \rho \partial_t \mathbf{u} + \rho \mathbf{u} \cdot \nabla \mathbf{u} + \nabla P - \mu \Delta \mathbf{u} - (\lambda + \mu) \nabla \text{div } \mathbf{u} = (\nabla \times \mathbf{b}) \times \mathbf{b}, \\ \partial_t \mathbf{b} + \mathbf{u} \cdot \nabla \mathbf{b} - \mathbf{b} \cdot \nabla \mathbf{u} + \mathbf{b} \text{div } \mathbf{u} = 0, \\ \nabla \cdot \mathbf{b} = 0, \end{cases}

其中 $\rho$ 为密度， $\mathbf{u}$ 为速度， $\mathbf{b}$ 为磁场， $P=P(\rho)$ 为压力函数， $\mu > 0$ 和 $\lambda$ 满足物理约束 $2\mu + 3\lambda \geq 0$ 。这里的“无电阻”意味着磁场方程中不包含磁扩散项（即没有 $-\nu_m \Delta \mathbf{b}$ ），这使得磁场方程是一个纯粹的输运-拉伸方程，缺乏自身的正则化机制。

文章关注的是平衡态 $(\rho, \mathbf{u}, \mathbf{b}) = (1, \mathbf{0}, \mathbf{e}_3)$ 附近的扰动，其中 $\mathbf{e}_3 = (0,0,1)^\top$ 是一个常向量背景磁场。引入扰动变量 $(a, \mathbf{u}, \mathbf{B}) = (\rho - 1, \mathbf{u}, \mathbf{b} - \mathbf{e}_3)$ ，线性化的耦合机制主要来自于算子 $\partial_{x_3}$ ，这是磁场稳定性的唯一来源。

在三维周期域中研究此类系统的稳定性面临巨大挑战。以往的工作（如 Chen, Zhang and Zhou [6], Wu and Zhai [34]）通常要求背景磁场方向 $\mathbf{n}$ 满足 Diophantine 条件：

|\mathbf{n} \cdot \mathbf{k}| \geq c |\mathbf{k}|^{-r}, \quad \forall \mathbf{k} \in \mathbb{Z}^3 \setminus \{0\}.

这一条件保证了高频模态能够被有效阻尼，但排除了诸如 $\mathbf{e}_3$ 这样的有理方向。对于有理背景场，存在无穷多频率使得 $\mathbf{n} \cdot \mathbf{k} = 0$ ，导致在这些频率上磁场无法通过粘性耗散获得控制。Pan, Zhou and Zhu [25] 针对不可压缩情形，利用对称性假设替代了 Diophantine 条件，证明了全球稳定性。然而，将这一策略推广到可压缩情形时，原有的机制失效，因为可压缩性引入了两个新的障碍：一是磁压项 $\frac{1}{2}\nabla |\mathbf{B}|^2$ 无法像不可压缩情形那样通过对称性抵消；二是密度依赖的非线性项 $I(a) = \frac{a}{1+a}$ 乘在高阶项上，且由于平均部分不衰减，任何将其视为微扰的估计都会失败。

因此，本文试图填补的空白是：在三维可压缩无电阻 MHD 系统中，对于有理背景场 $\mathbf{e}_3$ ，是否存在一种新的分析框架来克服 $$x_3$$ -独立扇区（即 $$k_3=0$$ 的频率分量）中缺乏耗散、阻尼和衰减的问题，从而证明经典解的全球适定性和非线性稳定性。

2. 主要定理或主要结果

本文的主要结果是定理 1.1，证明了在特定对称性假设下，三维可压缩无电阻 MHD 系统在全球范围内的适定性和稳定性。

定理 1.1: 设 $$s > 7$$ 为整数， $\sigma \in (0, 1/2)$ 。假设初始数据 $(a_0, \mathbf{u}_0, \mathbf{B}_0) \in H^s(\mathbb{T}^3)$ 满足 $\nabla \cdot \mathbf{B}_0 = 0$ 以及 $$x_3$$ -对称性条件 (1.3) 和零均值条件 (1.4)。存在 $\varepsilon > 0$ ，使得如果

\|(a_0, \mathbf{u}_0, \mathbf{B}_0)\|_{H^s} \leq \varepsilon,

则系统 (1.2) 存在唯一的全球经典解 $(a, \mathbf{u}, \mathbf{B})$ 。此外，存在常数 $$C > 0$$ ，使得对所有 $$t > 0$$ ，有如下估计：

一致有界性 (损失一阶导数): $\|(a, \mathbf{u}, \mathbf{B})(t)\|_{H^{s-1}} \leq C \varepsilon.$
振荡部分的衰减: 对于 $$i = 1, 2, 3$$ ，振荡部分 $\bar{f} = f - \int_{\mathbb{T}} f dx_3$ 满足多项式衰减率： $\|(\bar{a}, \mathbf{u}, \bar{\mathbf{B}})(t)\|_{H^{s-i}} \leq C \varepsilon (1+t)^{\frac{\sigma - i}{2}}.$ 注意：这里的速度 $\mathbf{u}$ 整体衰减，而密度 $$a$$ 和磁场 $\mathbf{B}$ 仅其振荡部分衰减，平均部分仅在时间上一致有界。

关键特征与技术细节:

有理背景场: 结果适用于背景场 $\mathbf{e}_3$ ，无需 Diophantine 条件。这是首个关于三维可压缩无电阻 MHD 系统在有理背景场下的稳定性结果。
各向异性衰减: 衰减行为反映了线性化系统的结构。 $$x_3$$ -独立扇区支持稳态解，因此平均部分 $(\bar{a}, \bar{\mathbf{B}})$ 不衰减是预期的，而非方法的缺陷。
高阶增长: 尽管 $H^{s-1}$ 范数一致有界，但最高阶范数表现出缓慢的增长： $\|(a, \mathbf{u}, \mathbf{B})(t)\|_{H^s} \leq C \varepsilon (1+t)^{\sigma/2}.$ 这种增长是由于缺乏平滑效应导致的，且增长率可以通过选择小的 $\sigma$ 任意慢。
隐藏阻尼机制: 核心发现是总压力 $\mathcal{D} = P(1+a) + B_3 + \frac{1}{2}|\mathbf{B}|^2$ 的平均值 $\overline{\mathcal{D}}$ 满足一个强阻尼波动方程： $\partial_t^2 \overline{\mathcal{D}} - 2\Delta \overline{\mathcal{D}} - \nu \Delta \partial_t \overline{\mathcal{D}} = \text{非线项},$ 其中 $\nu = \lambda + 2\mu$ 。该方程的主波部分以快磁声速 $\sqrt{2}$ 传播，其非抛物分支的阻尼率为 $-2/\nu$ ，且不涉及导数的增益。这一隐藏阻尼替代了缺失的磁耗散，并吸收了磁压项带来的障碍。

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

为了克服有理背景场下 $$x_3$$ -独立扇区的退化问题，作者发展了一套结合隐藏结构分析与加权能量方法的框架。以下是核心的技术工具和创新点：

1. 总压力 $\mathcal{D}$ 的隐藏波动结构与阻尼 这是本文最核心的创新。作者定义了精确的总压力变量：

\mathcal{D} = P(1+a) + B_3 + \frac{1}{2}|\mathbf{B}|^2.

这个变量的选取并非随意，它具有三重身份：它是热力学压力和磁压力的总和；它是动量方程中完整的压力力梯度；它是使平均扇区闭合的变量。

闭合性: 对平均变量 $\overline{\mathbf{u}}$ 和 $\overline{\mathcal{D}}$ ，线性化后的系统退化为一个双曲型可压缩 Navier-Stokes 系统： $\begin{cases} \partial_t \overline{\mathbf{u}} - \mu \Delta \overline{\mathbf{u}} - (\lambda+\mu)\nabla \text{div } \overline{\mathbf{u}} + \nabla \overline{\mathcal{D}} = \text{二次项}, \\ \partial_t \overline{\mathcal{D}} + 2 \text{div } \overline{\mathbf{u}} = \text{二次项}. \end{cases}$ 系数 $$2$$ 来源于声速平方 $$c_s^2$$ 和 Alfvén 速度平方 $$v_A^2$$ 之和 ( $$c_s^2+v_A^2 = 1+1=2$$ )。这意味着在平均扇区，磁场的影响完全消失，系统表现为具有声速 $\sqrt{2}$ 的粘性流体。
阻尼提取: 对上述系统取散度并利用 $\nabla \cdot \mathbf{B}=0$ ，得到 $\overline{\mathcal{D}}$ 满足 Kelvin-Voigt 阻尼波动方程。其色散关系 $\tau^2 + \nu |\mathbf{k}|^2 \tau + 2|\mathbf{k}|^2 = 0$ 的两个根中，一个趋于 $-\nu |\mathbf{k}|^2$ （抛物分支），另一个趋于 $-2/\nu$ （有限负极限）。后者提供了不依赖于高频的均匀阻尼，且不产生平滑效应（no gain of derivatives），这与无磁扩散的物理图像相容。在线性水平上，这种阻尼通过交叉项恒等式在能量估计中实现，具体表现为 $\| abla \overline{\mathcal{D}}\|^2_{\dot{H}^{s-1}}$ 项的正定贡献。

2. 振荡扇区的阻尼结构 对于 $$x_3$$ -振荡部分，作者沿用了类似 Diophantine 情形下的方法，但利用了 $$x_3$$ -对称性。

定义组合变量 $$d = a + B_3$$ 和 $G = Q\mathbf{u} - \frac{1}{\nu}\Delta^{-1}\nabla d$ 。
线性化显示 $$d$$ 以速率 $2/\nu$ 阻尼，且 $\partial_{x_3}\mathbf{B}$ 满足阻尼波动方程。
关键在于，由于对称性假设 (1.3)，横向磁场 $\mathbf{B}_h$ 和纵向速度 $$u_3$$ 的 $$x_3$$ -平均为零，因此它们等于其自身的振荡部分。这使得各向异性 Poincaré 不等式 $\|\bar{f}\|_{L^2} \lesssim \|\partial_{x_3}f\|_{L^2}$ 可以直接应用于完整场，而不需要额外的投影。这允许将 Lorentz 力非线性项 $\mathbf{B} \cdot \nabla \mathbf{B}$ 分解并控制在振荡扇区的阻尼范围内。

3. 移位时间权重能量泛函 由于平均部分 $(\bar{a}, \bar{\mathbf{B}})$ 不衰减，而振荡部分和速度衰减，不同扇区的自然时间尺度不同。作者构造了带有移位时间权重 $(T+t)^{i-\sigma}$ 的能量泛函 $$E_i(t)$$ 。

移位权重的作用: 使用 $$(T+t)$$ 而不是 $$(1+t)$$ 作为权重基底，其中 $$T$$ 是一个大参数。这在组合各个扇区的估计时，产生了小因子 $T^{2\sigma - 1/2}$ ，从而吸收掉由不同权重不匹配引起的线性耦合项。限制 $\sigma < 1/2$ 正是为了保证这个小因子的存在。
密度加权估计: 在处理振荡磁场 $\bar{\mathbf{B}}$ 的高阶估计时，直接能量法会损失一阶导数。作者采用密度加权范数 $\frac{1}{1+a}$ 来处理 $\bar{\mathbf{B}}$ 的估计，这导致了与速度方程中对应项的精确抵消，解决了密度依赖非线性项带来的困难。

4. 技术障碍的克服

磁压项的处理: 在不可压缩情形中，磁压项可以通过积分分部和对易关系抵消。但在可压缩情形中，这一机制失效。通过引入 $\mathcal{D}$ ，磁压项 $\frac{1}{2}\nabla |\mathbf{B}|^2$ 被吸收到总压力的结构中，不再需要单独估计，从而消除了可压缩性的主要障碍。
无 Diophantine 条件的补偿: 在没有 Diophantine 条件下，无法通过 $\|\mathbf{B}\|_{H^k} \lesssim \|\mathbf{n} \cdot \nabla \mathbf{B}\|_{H^{k+r}}$ 获得导数增益。作者通过分离平均和振荡扇区，分别利用隐藏阻尼（平均扇区）和各向异性 Poincaré 不等式（振荡扇区）来获得必要的控制，避免了全局的导数损失。

4. 与之前工作的比较

本文的结果在多个维度上推进了三维可压缩 MHD 系统的稳定性理论，并与既有文献形成了鲜明的对比：

1. 与 Diophantine 背景场工作的比较

Wu and Zhai [34] 等人在可压缩无电阻 MHD 系统中取得了全球稳定性，但其前提是背景磁场方向 $\mathbf{n}$ 满足 Diophantine 条件。这一条件虽然在数学上简化了高频分析，但在物理上排除了许多重要的情况，特别是所有有理方向（如 $\mathbf{e}_3$ ）。
本文的推进: 首次去除了 Diophantine 条件，针对有理背景场 $\mathbf{e}_3$ 证明了稳定性。代价是需要引入 $$x_3$$ -对称性假设 (1.3)，并且解的最高阶范数 $$H^s$$ 表现出缓慢的时间增长 $(1+t)^{\sigma/2}$ ，而不仅仅是 $H^{s-1}$ 的一致有界性。相比之下，Diophantine 情形下通常可以获得更高阶的一致有界性或更快的衰减。

2. 与 Pan, Zhou and Zhu [25] (不可压缩情形) 的比较

Pan et al. [25] 开创了利用对称性替代 Diophantine 条件的先例，但他们处理的是不可压缩 MHD 系统。
机制的差异: 在不可压缩情形中，对称性使得某些非线性项通过对称性精确抵消，且没有密度变化带来的额外复杂性。而在可压缩情形中，作者指出 [25] 中的机制“断裂”，原因有二：一是磁压项 $\frac{1}{2}\nabla |\mathbf{B}|^2$ 无法抵消；二是密度依赖项 $$I(a)$$ 使得微扰论证失效。
本文的创新: 作者通过发现总压力 $\mathcal{D}$ 的隐藏波动结构，成功地将不可压缩情形的对称性策略推广到了可压缩情形。 $\mathcal{D}$ 不仅吸收了磁压项，还揭示了平均扇区中存在的隐藏阻尼，这是 [25] 中不存在的新机制。

3. 与二维工作的比较

Wu and Wu [32], Wu and Zhu [36] 等在二维情形下已经证明了全球稳定性，且不需要对称性或 Diophantine 条件。
维度差异: 二维几何本质上更宽容，单一背景方向可以控制余维为一的频率集。而在三维中，除非背景场是无理的，否则存在大量未被控制的频率。本文通过引入对称性，有效地将三维问题降维处理，使其在平均扇区表现出一维或二维的特征（即 $$x_3$$ -独立），从而克服了三维的困难。

4. 与其他部分耗散 MHD 工作的比较

Chen, Zhang and Zhou [6] 等人研究了仅具有粘性或仅具有电阻率的 MHD 系统。这些工作通常也依赖于 Diophantine 条件或特定的几何结构。
本文的定位: 本文专注于无电阻但有粘性的情形，且背景场为有理方向。它展示了即使在缺乏磁扩散的情况下，只要存在适当的对称性和背景场结构，粘性耗散仍可通过耦合机制传递到磁场，从而实现全球稳定性。这丰富了我们对部分耗散流体系统稳定性的理解。

5. 局限性与开放问题

对称性假设: 本文的结果依赖于 $$x_3$$ -对称性 (1.3)。作者指出，如果没有这一对称性，对于有理背景场的三维周期 MHD 系统，全球稳定性仍然是一个开放问题。这表明对称性在克服有理方向退化中起到了至关重要的作用。
高阶增长: 虽然 $H^{s-1}$ 范数一致有界，但 $$H^s$$ 范数的缓慢增长表明解可能在极长时间尺度上失去光滑性（尽管全球存在性已证）。这与 Diophantine 情形下可能获得的更强正则性形成对比。
适用范围: 结果仅限于近平衡态的小扰动。大扰动情形下的行为尚不清楚。

11. 可压缩欧拉方程的结构稳定奇点与 Lipschitz 稳定的最优传输度量

Structurally stable singularities and Lipschitz stable optimal transport metrics for the compressible Euler equations

作者: Geng Chen, Yanbo Hu, Yannan Shen

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13943v2

摘要翻译: 众所周知，可压缩欧拉方程的解可能在有限时间内发展出奇异性。本文对一维可压缩欧拉方程在一般光滑初值下、直至首次奇性发生时刻的解的行为进行了细致分析。

Comment: 67 pages, 4 figures

链接: http://arxiv.org/abs/2609.13943v2

1. 研究对象和背景

本文研究一维等熵可压缩欧拉方程的初值问题。在拉格朗日坐标下，该方程组表示为

\begin{cases} \tau_t - u_x = 0, \\ u_t + p(\tau)_x = 0, \\ p(\tau) = \frac{1}{\gamma} \tau^{-\gamma}, \end{cases}

其中 $\tau > 0$ 是比容， $$u$$ 是速度， $p(\tau)$ 是多原子气体的状态方程，绝热指数 $\gamma > 1$ 。初始条件设为 $\tau(x, 0) = \tau_0(x)$ ， $$u(x, 0) = u_0(x)$$ 。

可压缩欧拉方程作为流体力学中的基本模型，其核心难点在于即使对于足够光滑且小的初始数据，由于系统的拟线性结构，解也可能在有限时间内发展出奇点（梯度爆破）。精确刻画首次奇点形成时刻（即“预激波”或 pre-shock）的解的行为，是研究后续激波演化的关键基础。

关于一维可压缩欧拉方程的奇点形成，已有大量文献。Riemann 早在 1850 年代就指出了这一现象。Lax 基于特征线方法研究了 $2 \times 2$ 可约齐次双曲系统的平滑解崩溃。对于 $\gamma \geq 3$ 的情况，已有完整的二分法结果：当且仅当初始数据包含压缩点时，有限时间奇点才会形成。然而，对于物理上更相关的 $1 < \gamma < 3$ 情况，由于密度消失导致的 Riccati 方程退化，需要密度的适当下界估计来解决。Chen 等人此前已针对此情形推导了最优的时间依赖密度下界，并获得了有限时间奇点形成的完整图像。

近年来，激波形成的构造性证明受到了广泛关注。Lebaud 对一维等熵欧拉方程进行了开创性的构造分析，在简单波初始数据和非退化临界点的假设下，证明了 $C^{1/3}$ 尖峰型预激波的形成。随后，Chen 和 Dong 去除了简单波假设，Kong 将该理论推广到一般 $2 \times 2$ 真非线性双曲系统。Yin 处理了非等熵情形，证明了三维球对称小扰动初始数据下的梯度崩溃和 $C^{1/3}$ 预激波结构。Christodoulou 及其合作者发展了几何框架，用于研究多维相对论无旋流和多维可压缩欧拉方程的激波形成。Buckmaster, Shkoller, Vicol 等人通过一系列工作建立了多维可压缩非等熵欧拉方程激波形成与发展的变革性构造理论，特别是利用了调制自相似分析方法。

尽管已有许多关于特定类型初始数据（如简单波、球对称、轴对称或小扰动）下奇点结构的构造性结果，但对于一般光滑初始数据下首次奇点形成时的解的通用结构和渐近行为，尚缺乏系统的描述。此外，虽然 $C^{1/3}$ 尖峰型奇点在多个特殊设定中被发现，但其是否为一类“通用”（generic）奇点，以及在何种意义下是稳定的，仍未得到充分阐明。

本文旨在填补这一空白。首先，利用 Thom 横截性定理，证明对于开稠密的 $$C^3$$ 初始数据集合，解在首次奇点形成时刻除有限个点外是二阶连续可微的。其次，对这些通用初始数据，给出每个奇点处解在半邻域内的精确渐近描述，验证其为 Hölder 指数 $$1/3$$ 的尖峰型奇点。最后，针对具有小 BV 范数的光滑初始数据，构造两种 Finsler 型最优传输度量，并证明在这些度量下，解直到首次奇点形成时刻对初始数据的依赖是 Lipschitz 连续的，且 Lipschitz 常数有界。这表明 $C^{1/3}$ 通用奇点在此意义下是稳定的，而其他指数的尖峰型奇点是不稳定的。

2. 主要定理或主要结果

本文的主要结果分为三部分，分别对应 Theorem A, B, C。

Theorem A (通用正则性): 设 $\bar{\tau}_0 > 0$ 为常数。存在一个开稠集 $D \subset (C^3(\mathbb{R}) \cap W^{1,1}(\mathbb{R}))^2$ ，满足

\frac{1}{c} \leq \tau_0(x) \leq c,

以及

\min \left\{ \inf \varpi'_{-,0}(x), \inf \varpi'_{+,0}(x) \right\} < 0,

其中 $\varpi_{\pm,0}(x) = u_0(x) \pm \frac{2}{\gamma-1}\tau_0(x)^{-(\gamma-1)/2}$ 是 Riemann 不变量。对于满足 $(\tau_0(x) - \bar{\tau}_0, u_0(x)) \in D$ 和上述条件的初始数据，欧拉方程 (1.1)-(1.2) 的解 $(\tau(x,t), u(x,t))$ 在 $\mathbb{R} \times [0, t^*]$ 上有定义，且在直线 $$t = t^*$$ 上除有限个点外是二阶连续可微的。这里的 $$t^*$$ 是首次奇点形成时间。

Theorem B (奇点的精确渐近结构): 考虑满足 Theorem A 中条件的初始数据，设 $$(x^*, t^*)$$ 是由平滑解形成的一个奇点。则 $$(x^*, t^*)$$ 是一个 Type-1 奇点（定义见 Section 5）。如果梯度 $\partial_x \varpi_-(x,t)$ 在 $$(x^*, t^*)$$ 处爆破，则另一个 Riemann 不变量 $\varpi_+(x,t)$ 在该点是平滑的。并且存在常数 $\zeta \neq 0$ ，使得对于 $t \leq t^*$ ，有如下渐近展开：

\tau(x, t) = \tau(x^*, t^*) + \zeta \cdot [c(\tau(x^*, t^*))(t - t^*) + (x - x^*)]^{1/3} + O(1) \left( |t - t^*| + |x - x^*| \right)^{4/9},

u(x, t) = u(x^*, t^*) + c(\tau(x^*, t^*))\zeta \cdot [c(\tau(x^*, t^*))(t - t^*) + (x - x^*)]^{1/3} + O(1) \left( |t - t^*| + |x - x^*| \right)^{4/9}.

这表明解在每个奇点处具有 $C^{1/3}$ Hölder 类的尖峰型预激波结构。

Theorem C (Lipschitz 稳定性与最优传输度量): 设 $\bar{\tau}_0, \bar{c}$ 和 $\eta$ 为正常数，定义集合

D_\eta = \left\{ (\tau_0(x), u_0(x)) \middle| (\tau_0 - \bar{\tau}_0, u_0) \in (C^1(\mathbb{R}) \cap W^{1,1}(\mathbb{R}))^2, \frac{1}{\bar{c}} \leq \tau_0(x) \leq \bar{c}, \int_\mathbb{R} [|\varpi'_{-,0}(x)| + |\varpi'_{+,0}(x)|] dx \leq \eta \right\}.

设 $(\tau(x,t), u(x,t))$ 和 $(\hat{\tau}(x,t), \hat{u}(x,t))$ 是两个初始数据在 $D_\eta$ 中的欧拉方程解，对应的首次奇点形成时间分别为 $$t^*$$ 和 $\hat{t}^*$ 。令 $t_d = \min\{t^*, \hat{t}^*\}$ 。若 $\eta$ 足够小，则对于任意 $t \leq t_d$ ，两个测地距离 $d_*^{(i)}(\cdot, \cdot)$ 和 $d_*^{(ii)}(\cdot, \cdot)$ （定义见 Section 7 的 Definition 5）满足：

d_*^{(\iota)}((\tau, u), (\hat{\tau}, \hat{u})) \leq {}^b C d_*^{(\iota)}((\tau_0, u_0), (\hat{\tau}_0, \hat{u}_0)), \quad (\iota = i, ii),

其中 ${}^b C$ 是某个正常数。这意味着在这些度量下，只要初始数据的 BV 范数足够小，解直到首次奇点形成时刻对初始数据的依赖是 Lipschitz 连续的。

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

本文的证明策略结合了变量变换、横截性理论和几何度量构造。

1. 半线性系统与坐标变换: 为了处理梯度的爆破，作者引入了新的因变量和自变量。首先，将原系统转化为关于 Riemann 不变量梯度 $(R, S) = (\partial_x \varpi_-, \partial_x \varpi_+)$ 的系统。然后，引入角度变量 $\alpha, \beta \in [-\pi, \pi)$ ，使得 $$R$$ 和 $$S$$ 的爆破对应于 $\alpha = -\pi$ 和 $\beta = -\pi$ 。受 Bressan 和 Zheng 以及 Chen 和 Shen 工作的启发，进一步引入新的因变量 $$(p, q)$$ 和坐标变量 $$(X, Y)$$ 。这一变换将原拟线性系统转化为在 $$(X, Y)$$ 坐标下关于未知量 $(\alpha, \beta, p, q, \tau, u)$ 的等价半线性系统 (3.14) 或 (3.15)。关键在于，这个半线性系统在边界 $$t(X,Y) = t^*$$ 处保持平滑，而奇点对应于雅可比矩阵 $\partial(x,t)/\partial(X,Y)$ 不可逆的点，即 $\alpha = -\pi$ 或 $\beta = -\pi$ 。

2. Thom 横截性定理的应用: 为了证明通用正则性（Theorem A），作者在 Section 4 中应用了 Thom 横截性定理。他们证明了如果 $$(X_0, Y_0)$$ 是边界上的一个奇点，且满足某些非退化条件（例如 $(\alpha, \alpha_X, \alpha_{XX})(X_0, Y_0) \neq (-\pi, 0, 0)$ ），则存在 3 参数族的平滑解，使得这些解在 $$(X_0, Y_0)$$ 处的相关导数张量的秩为 3。基于此，结合 Thom 横截性定理，证明了排除特定退化情况的 $$C^2$$ 解集合在 $$C^2$$ 空间中是开稠的。这直接导致了 Theorem A 的结论，即对于开稠集的初始数据，奇点只能是有限的，且没有高阶退化。

3. 奇点分类与泰勒展开: 在 Section 5 中，作者将边界上的奇点分为三类：Type-1 ( $\alpha=-\pi, \alpha_X=0, \beta \neq -\pi, \alpha_{XX} \neq 0$ )，Type-2 ( $\alpha=-\pi, \alpha_X \neq 0, \beta \neq -\pi$ )，和 Type-3 ( $\alpha=-\pi, \beta=-\pi, \alpha_X \neq 0, \beta_Y \neq 0$ )。通过分析控制方程 (3.14)，作者进行了严格的泰勒展开分析。他们证明了由于 $$t^*$$ 是首次奇点形成时间，只有 Type-1 奇点可以在边界 $$t(X,Y)=t^*$$ 上出现，而 Type-2 和 Type-3 被排除。对于 Type-1 奇点，利用 $$t(X,Y)$$ 和 $$x(X,Y)$$ 的泰勒展开将 $$X, Y$$ 表示为 $$t, x$$ 的函数，再代入 $\tau(X,Y)$ 和 $$u(X,Y)$$ 的展开式，从而得到了 Theorem B 中精确的 $$1/3$$ 次幂渐近公式。

4. Finsler 型最优传输度量的构造: 为了建立 Lipschitz 稳定性（Theorem C），作者在 Section 6 和 7 中构造了两种 Finsler 型规范。对于两个平滑解，考虑连接它们的所有可能路径，路径的长度由切向量的范数积分定义。关键在于如何定义 Finsler 范数 $\| \cdot \|_{(\tau_\vartheta, u_\vartheta)(t)}$ ，使其能够捕捉欧拉方程的动态特性，并满足不等式 $\| \Upsilon_t \| \leq {}^b C \| \Upsilon_0 \|$ 。作者分析了切向量的一阶变化，构造了依赖于解的范数，并证明了其在有界子集上的 Lipschitz 连续性。最后，将这些局部度量推广到一般情况，定义了测地距离 $d_*^{(i)}$ 和 $d_*^{(ii)}$ ，并比较了这些度量与 Sobolev 度量和 $$L^1$$ 度量的关系。

4. 与之前工作的比较

本文与既有工作的主要区别和推进体现在以下几个方面：

1. 从构造性结果到通用性结果: 此前的许多工作（如 Lebaud [49], Chen and Dong [33], Kong [46], Yin [64] 等）都是在特定的初始数据假设下（如简单波、球对称、轴对称或小扰动）构造性地证明了 $C^{1/3}$ 尖峰型奇点的形成。这些结果通常依赖于初始数据满足某些逐点条件或对称性。相比之下，本文的结果是针对一般光滑初始数据（在一个开稠集内）建立的。作者利用 Thom 横截性定理证明了这种 $C^{1/3}$ 尖峰型奇点是“通用”的，即对于绝大多数初始数据，首次奇点都会呈现这种结构。这大大扩展了结果的适用范围，不再局限于特殊的构造性场景。

2. 方法论的创新: 半线性系统与横截性 vs. 自相似分析: Buckmaster, Shkoller, Vicol 等人 [13-17, 56-58] 在处理多维可压缩欧拉方程的激波形成时，主要采用了调制自相似分析方法（modulated self-similar analysis）来建立显式爆破剖面的稳定性。这种方法非常强大，但通常适用于特定的对称设定或小扰动背景。本文采用了一种不同的方法，即通过将原系统转化为半线性系统，并结合 Thom 横截性定理来分析通用奇点的结构。这种方法类似于 Bressan 和 Chen [3, 4, 11] 在处理变分波方程时使用的框架。作者明确指出，这提供了一种不同于 Buckmaster 等人的新方法来研究预激波解，并且特别适用于一般初始数据。

3. 稳定性概念的深化: 最优传输度量下的 Lipschitz 稳定性: 虽然之前的工作（如 Neal, Shkoller, Vicol [56, 58]）也讨论了奇点形成的稳定性，但通常是在 Sobolev 空间或特定的拓扑下讨论的。本文的一个显著贡献是构造了两种 Finsler 型最优传输度量 $d_*^{(i)}$ 和 $d_*^{(ii)}$ ，并证明了在这些度量下，解对初始数据的依赖是 Lipschitz 连续的，且常数有界。这与传统的 $$L^1$$ 度量不同，后者只比较函数值，而忽略了导数的差异。由于激波形成后导数会变成 delta 函数，简单的 $$L^1$$ 度量不足以捕捉这种动态，因此需要引入“传输”或“移位”概念的最优传输度量。本文的度量包含了导数之间的差异，属于更高阶的度量，从而更精确地刻画了奇点形成过程中的稳定性。这也回应了 Remark 7 中指出的，这些结果不依赖于波速 $$c$$ 的具体形式，可以推广到一般的拟线性波动方程。

中文标题: 非散度型动力学方程的 Harnack 估计

摘要翻译: 我们建立了非散度型动力学方程的非负强上解的弱 Harnack 不等式和非负强解的 Harnack 不等式：

\partial_tu+v\cdot\nabla_xu+a(t,x,v):\nabla_v^2u +b(t,x,v)\cdot\nabla_vu+c(t,x,v)u=0,

其中系数 $$a$$ 、 $$b$$ 和 $$c$$ 仅为 Borel 可测，矩阵 $$a$$ 一致椭圆，且 $$b$$ 和 $$c$$ 有界。作为推论，我们得到了强解的内部动力学 Hölder 估计。我们的方法发展了 Krylov–Safonov 方法的概率动力学类似物。主要的新成分是对具有渐进可测系数的动力学 Itô 过程的 Krylov 估计的一个定量版本，这进而给出了与相关动力学随机微分方程（具有有界可测系数）的弱解的全局存在性。这些结果似乎是首个针对仅具可测系数的动力学方程的 Harnack 和 Hölder 正则性理论。

5. 二维可压缩微极流体方程中大体粘滞与微小角粘滞的联合极限

The combined limits of large bulk viscosity and small angular viscosity for the 2D compressible micropolar fluid equations

作者: Fangbin Lin, Shengquan Liu, Jianwen Zhang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15518v1

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

摘要翻译

中文标题: 二维可压缩微极流体方程中大体粘滞与微小角粘滞的联合极限

摘要翻译: 本文研究二维可压缩微极流体方程中大体粘滞与微小角粘滞的联合极限。我们首先建立了当初始速度在 $rac{ar u}{ u} abla u_0$ 意义下是良好准备的（即 $rac{ar u}{ u} abla u_0$ 有界）时，大体粘滞条件下柯西问题强解的全局适定性。所得到的全局先验估计关于体粘滞 $$ν$$ 和角粘滞 $$ar u$$ 是一致的。基于这些一致估计，我们随后证明了当 $$ν oar u$$ 且 $$ar u o0$$ 时的联合极限，并表明可压缩微极系统的解收敛于粘性且无扩散不可压缩微极系统的解。我们还推导了不可压缩极限和角粘滞消失极限的收敛率。具体而言，一般情形下不可压缩极限的收敛率为 $ν^{-1/4}$ 阶，而在对初始密度附加一个衰减条件时可改进为 $ν^{-1/2}$ 阶。与先前工作相比，角粘滞消失极限是在没有额外技术假设 $$ar u/ u^eta<ar u$$ 对于某个 $$1/2ar ueta<ar u$$ 的情况下建立的，这将相关结果从不可压缩情形推广到了可压缩设定。作为副产品，我们得到了二维可压缩微极方程大解的全局正则性，无论是否包含角粘滞项，前提是体粘滞足够大。

6. 半线性波动方程弱解的尖锐非唯一性

Sharp non-uniqueness of weak solutions to the semilinear wave equation

作者: Peng Qu, Mingxin Zhang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15442v1

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

摘要翻译

中文标题: 半线性波动方程弱解的尖锐非唯一性

摘要翻译: 在本文中，我们建立了半线性波动方程

中文标题: 正交各向异性、宽退化双非线性扩散方程的有界性与收缩估计

摘要翻译: 我们研究如下形式的双非线性正交各向异性演化方程弱解的正则性：

\partial_{t}(\vert u\vert^{α-1}u)-\sum_{i=1}^{N}\partial_{i}\left[a_{i}(x,t)\,(|\partial_{i}u|-δ_{i})_{+}^{p-1}\frac{\partial_{i}u}{\vert\partial_{i}u\vert}\right]=f\,\,\,\,\,\,\,\mathrm{in}\,\,Ω_{T}=Ω\times(0,T),

其中 $$Ω$$ 是 $\mathbb{R}^{N}$ （ $N\geq2$ ）中的有界开集，系数 $a_{i}$ 是可测且有界的， $$α>0$$ ， $p\in(1,\infty)$ ，且 $δ_{1},\ldots,δ_{N}$ 是非负数。我们证明，由于弱解属于一个适当的 De Giorgi 型能量类，因此它们在局部范围内是有界的。我们还获得了与上述偏微分方程相关的 Cauchy 问题解在空间上的全局有界性以及收缩估计。我们的分析扩展了文献中关于扩散方程的类似结果，这些方程要么不具备双非线性特征，要么其退化程度低于本文所考虑的方程。本文的另一主要创新在于方程右端存在源项 $$f$$ ，我们对时空变量施加了适当的可积性假设。

10. Vlasov–Yukawa 方程非线性 Landau 阻尼中的大盒子过渡

Large-Box Transition in Nonlinear Landau Damping for the Vlasov–Yukawa Equation

作者: Ling-Bing He, Yue Luo

类别: math.AP

中文标题: 半密度下的 Shubin–Vakilian–Wolff 不确定性原理

摘要翻译: Shubin、Vakilian 和 Wolff 证明了关于在倒数尺度上具有足够小局部密度的集合的傅里叶不确定性原理，并提出了是否所有低于一的都是可容许密度这一问题。我们在一维中通过构造一系列 $$1/2$$ -密度薄集对和单位向量序列来否定地回答该问题，这些序列在这些集合之外的总位置质量和傅里叶质量趋于零。这一障碍对于每个缩放倒数剖面 $ρ_κ(x)=\min\{1,κ/|x|\}$ , $$κ>0$$ 均持续存在。另一方面，对于 $0<κ\le1$ , 只要密度满足

中文标题: 粘性 Burgers 方程的交叉渐近与精确约束速率

摘要翻译: 我们研究 $$(0,L)$$ 上具有保守边界条件 $$u_x=-u^2$$ 的一维粘性 Burgers 方程。在半直线情形，解趋向一个非线性自相似剖面；而在有界区间上，解收敛到一个非平凡平衡态。我们显式地刻画了动力学在临界扩散尺度 $$t=cL^2$$ 处如何在这两种状态之间过渡。对于紧支集且质量为 $$M$$ 的初值，重标度后的解当 $L\to\infty$ 时收敛到一个显式的交叉剖面 $$Φ_c$$ 。在自相似变量下， $$Φ_c$$ 当 $c\downarrow0$ 时收敛到半直线剖面 $$f_M$$ ；经过重标度至区域变量后，它当 $c\to\infty$ 时收敛到区间平衡态。我们还确定了约束出现的精确阈值：

\lim_{c\downarrow0}-c\log|Φ_c(0)-f_M(0)|=1\qquad(M\ne0).

此外，在自相似变量的紧集上，区间解与半直线解在 $$C^k$$ 范数下的差异至多为 $C_{k,\varepsilon}\exp(-(1-\varepsilon)L^2/t)$ ，且该估计关于 $$L$$ 一致成立，指数常数是最佳的。因此，我们不仅识别出过渡尺度 $$L^2$$ ，还给出了支配交叉行为的剖面以及远端边界变得可见的精确速率。最后，我们讨论了空间依赖扩散系数下仍然成立的结论，并对三种渐近机制进行了数值说明。

不相关论文（基本信息）

1. $$C^1$$ 狭缝域中调和函数的边界正则性

Boundary regularity of harmonic functions in $$C^1$$ slit domains

作者: Joan Domingo-Pasarin, Pablo Hidalgo-Palencia, Alejandro Martínez, Clara Torres-Latorre

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15923v1

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

摘要: We establish precise upper and lower estimates for harmonic functions vanishing on the slit of a $$C^1$$ slit domain, with no assumption that the slit lies in a hyperplane. The classical $\sqrt{d}$ growth near the edge, $$d$$ being the distance to the slit, persists in this generality, up to an explicit factor

\exp\Big( \pm C \int_ρ^r ω(s)\, \frac{ds}{s} \Big)

determined by the $$C^1$$ -modulus of continuity $$ω$$ of the slit, where $$ρ< r$$ are the two scales being compared. The upper estimates allow a right-hand side and non-zero boundary data. The correction factors remain bounded above and below by positive constants as $ρ\to0$ precisely when $$ω$$ satisfies the Dini condition. Moduli of continuity beyond the Dini regime, as is the case for the logarithmic moduli arising at singular sets in relevant free boundary problems, were not covered by the previous $C^{1,α}$ theory. Previously, the $\sqrt{d}$ growth was known for slits contained in a hyperplane, which additionally have a $C^{1,α}$ edge (De Silva, Savin). For Lipschitz slits, there are boundary Harnack principles, but no growth rate is identified precisely. The main technical ingredient is a change of coordinates flattening a Lipschitz slit domain onto the model half-hyperplane slit, with quantitative estimates up to second order. The construction is geometric and does not use the equation, so we expect it to be useful for other boundary regularity problems.

2. 双幂非线性定常 Schrödinger 方程正解的唯一性：固定频率情形与固定质量情形

Uniqueness of positive solutions of double-power nonlinear stationary Schrödinger equations: the fixed frequency case and the fixed mass case

作者: Linjie Song, Wenming Zou

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15882v1

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

摘要: The existence of a positive radial solution to the problem \begin{align*} -Δu + λu = u^p + u^q, \quad u \in H^1(\mathbb{R}^N), ~~ N \ge 2, \end{align*} where $$λ> 0$$ and $$1 < q < p < 2^* - 1$$ , has been known for a long time. For the pure-power nonlinearity, it is well known that this solution is unique. In contrast, the uniqueness problem for the double-power case is substantially more delicate and depends on the dimension and the exponents. It has been known that the uniqueness result is in general not true for the equation with double-power nonlinearities in dimension three and it is a long-standing open question to prove the uniqueness of positive radial solutions throughout the full subcritical range for $N \ge 4$ since the work of H. Berestycki, P.-L. Lions [2] (Arch. Ration. Mech. Anal. 1983). In this paper we provide the uniqueness results for all $$λ> 0$$ and all $$1 < q < p < 2^* - 1$$ when $N \ge 5$ . In dimensions $$N = 2,3,4$$ , with additional conditions on $$p$$ and $$q$$ , we obtain the uniqueness results for all $$λ> 0$$ . Our results are in sharp contrast to those of J. Dávila, M. del Pino and I. Guerra [6] (Proc. London Math. Soc. 2013), where some non-uniqueness results were obtained for $$N=3$$ . At the same time, we give a positive answer and rigorous proof to the implication of their numerical simulation, clarify the open issues left by [6]. Finally, we completely resolve the uniqueness problem of the above equation with $$L^2$$ -mass constraints in the mass-subcritical and mass-critical regimes for all $N \ge 2$ .

3. k-凸超曲面的整体曲率估计

Global Curvature Estimates for k-Convex Hypersurfaces

作者: Fengrui Yang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15858v1

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

摘要: Establishing global curvature estimates for $$k$$ -convex solutions of prescribed curvature equations is a longstanding problem in fully nonlinear partial differential equations and geometric analysis. We develop a new approach that combines fractional-linear transformations with the logarithmic concavity of hyperbolic polynomials. The key step is a coercive estimate for the quadratic forms arising from the third-order terms in the maximum-principle argument. Together with a complementary concavity inequality, this estimate yields global curvature bounds for closed, strictly star-shaped $$k$$ -convex hypersurfaces in $\mathbb{R}^{n+1}$ satisfying $$σ_k(κ)=f(X,ν)>0$$ , throughout the range $3\leq k<n<2k$ .

4. Gurtin-MacCamy 种群模型中的双稳态波前

Bistable wavefronts in the Gurtin-MacCamy population model

作者: Quentin Griette, Franco Herrera, Hao Kang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15853v1

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

摘要: In this paper we prove the existence of monotone traveling wave solutions for an age-structured population model with a monotone and bistable nonlinearity. We formulate the problem in terms of a discrete-time recurrence $u_{n+1} = Q[u_n]$ , with operator $$Q$$ encoding the role of the speed into its definition. For each speed, we prove the existence of a traveling wave by the general theory of monotone evolution systems, and from appropriate estimates on the speed, we deduce the existence of a traveling wave of the original problem. Also, the sign of the speed is established under an invasion condition, and the asymptotic behavior of the profile is determined by some Tauberian results.

5. 一致椭圆偏微分方程的尖锐梯度估计、特征值渐近性质与 Landis 猜想

Sharp gradient bounds for uniformly elliptic PDE, eigenvalue asymptotics, and the Landis conjecture

作者: Boyan Sirakov, Philippe Souplet

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15821v1

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

摘要: We study three classical problems in the theory of linear second-order uniformly elliptic equations: (i) gradient estimates for the Dirichlet problem, (ii) estimates and asymptotics for the first eigenvalue of an elliptic operator, and (iii) the Landis conjecture on exponential decay. Our main results give versions of (i) in which the constant is sharply specified in terms of the norms of the coefficients of the operator and the size of the domain; and use these to strongly improve on known results for (ii) and (iii) by allowing both more general operators and weaker regularity assumptions on the coefficients, and by giving quantitative estimates. The proofs use a unified approach, relying on three ingredients: interior and boundary Harnack inequalities with sharp constants, duality arguments, and a $$C^1$$ estimate based on a rescaling procedure. The sharpness of the gradient and spectral estimates is demonstrated through various counterexamples.

6. $L^p(\mathbb{R}^N)$ 中由线性分数散度型算子生成的半群：常系数情形

The semigroup generated by linear fractional divergence form operators in $L^p(\mathbb{R}^N)$ : the constant coefficient case

作者: Alessandro Carbotti

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15809v1

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

摘要: Let $$0<s<1$$ and $$A$$ be a constant symmetric positive definite matrix. Given the quadratic form

Q_A(u)=\int_{\R^N} A\nabla^s u(x)\cdot \nabla^s u(x)\,dx,

where $\nabla^s$ denotes the Riesz fractional gradient, we establish necessary and sufficient conditions such that $$Q_A$$ defines a Dirichlet form. In strike contrast with the local case $$s=1$$ , positive definiteness of a constant matrix is not sufficient for the Markov property when $$0<s<1$$ . We also describe the associated symmetric Markov semigroup and its $$L^p$$ realizations, prove the existence of a smooth nonnegative heat kernel and establish a bounded $H^\infty$ functional calculus.

7. $\mathbb{R}^{7}$ 中的稳定 Bernstein 定理

The stable Bernstein theorem in $\mathbb{R}^{7}$

作者: Han Hong, Haizhong Li, Gaoming Wang

类别: math.DG, math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15720v1

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

摘要: We give a Green-function proof of the stable Bernstein theorem in $\mathbb R^7$ for smooth, connected, complete, two-sided minimal hypersurface, thus resolving the last case in stable Bernstein problem.

8. 抛物分数 $$p$$ -拉普拉斯方程的 Harnack 不等式

Harnack Inequality for Parabolic Fractional $$p$$ -Laplace Equations

作者: Harsh Prasad

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15694v1

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

摘要: We establish an intrinsic Harnack inequality with an optimal tail term for weak solutions of parabolic fractional $$p$$ -Laplace equations with $$p>2$$ and $s\in(0,1)$ whose kernels are symmetric, measurable, and comparable to $|x-y|^{-n-sp}$ . It constitutes the nonlocal analogue of the intrinsic Harnack inequality of DiBenedetto, Gianazza and Vespri. In contrast to the exponential change of variables and expansion-of-positivity used in the local case, our technique is nonlocal and relies on comparison arguments and barriers in a variational framework. Indeed, we obtain the stronger conclusion where the forward comparison starts immediately after the reference time, without an additional time gap, which is false in the local case. The latter leads to elliptic-type Harnack estimates for parabolic equations. Our result is new even in the linear case.

9. 具有能量界约束的Allen-Cahn方程有限Morse指数解

Finite Morse index solutions of the Allen-Cahn equation with energy bound

作者: Yong Liu, Kelei Wang, Juncheng Wei, Ke Wu

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15688v1

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

摘要: We study bounded entire solutions of the Allen-Cahn equation in dimension three with quadratic energy growth. We prove that such a solution of finite Morse index I has at most 2I ends. Furthermore, we show that Morse index one solutions have two ends and are axially symmetric, which generalizes an earlier result by Florit-Simon.

10. 用于四旋翼飞行器的体积谐波场导航

Volumetric Harmonic Field Navigation for Quadrotors

作者: Shuxiu Jia, Amartya Mukherjee, Yating Yuan

类别: cs.RO, math.AP, math.OC

发布日期: 2026-09-14

arXiv ID: 2609.15680v2

链接: http://arxiv.org/abs/2609.15680v2

摘要: Quadrotor navigation in cluttered 3-D environments requires global guidance while local motion remains subject to collision and motion limits. Harmonic potentials provide dense guidance from a global boundary value problem, but coupling a volumetric harmonic field to constrained physical quadrotor motion remains an open experimental problem. We couple a precomputed volumetric harmonic field with a constrained predictive planner that queries the field at predicted positions instead of extracting a global reference path. In Structured 3-D tests, harmonic guidance yields larger minimum clearance and lower RMS jerk than matched Dijkstra guidance, at the cost of longer paths; the same pattern remains when both methods use the same passage. Long maze tests span routes far beyond one prediction horizon, and Crazyflie trials validate physical execution. To the best of our knowledge, this is the first physical quadrotor demonstration of volumetric harmonic field navigation. The results show that globally constructed harmonic guidance can directly support local constrained motion generation on a physical quadrotor.

11. 带控制扩散的抛物型 Bellman 方程中 Howard 迭代的全局 Sobolev 收敛性

Global Sobolev convergence of Howard iteration for parabolic Bellman equations with controlled diffusion

作者: Lihua Bai, Linyu Miao

类别: math.AP, math.OC

发布日期: 2026-09-14

arXiv ID: 2609.15607v1

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

摘要: We prove global Sobolev convergence of Howard policy iteration for finite-horizon Bellman equations with control-dependent diffusion. In one dimension, the result holds for bounded measurable coefficients under uniform ellipticity, without a large discount, a short horizon, a small diffusion perturbation, or regularity assumptions on the improving policies. We also obtain multidimensional extensions under structural conditions on the diffusion matrix. The key observation is that vanishing policy improvements force the Bellman residual to converge in measure; higher integrability then yields strong convergence and constructs the unique strong solution. We extend the argument to entropy-regularized models, including zero-temperature limits, and establish quadratic convergence for special one-dimensional models.

12. 在给定面积的平面有界凸域中，三角形使 Monge-Ampère 特征值最小

Triangles minimize the Monge-Ampère eigenvalue among planar bounded convex domains of a given area

作者: Chong Gu

类别: math.AP, math.DG

发布日期: 2026-09-14

arXiv ID: 2609.15563v2

链接: http://arxiv.org/abs/2609.15563v2

摘要: We prove that triangles minimize the Monge-Ampère eigenvalue among bounded planar convex domains of a given area, thus resolving the two-dimensional case of Le’s simplex conjecture for the Monge-Ampère eigenvalue (The eigenvalue problem for the Monge-Ampère operator on general bounded convex domains, Ann. Sc. Norm. Super. Pisa Cl. Sci. (5) 18 (2018)). The proof uses a neighbor-decreasing property of the Monge-Ampère eigenvalues of convex polygons along vertex elimination parallel chord movements.

13. 归一化抛物型 $$p$$ -Laplace 方程的时间导数估计

Time derivative estimates for normalized parabolic $$p$$ -Laplace equations

作者: Se-Chan Lee, Hyungsung Yun

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15501v1

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

摘要: We establish the local boundedness of time derivatives of bounded viscosity solutions to the normalized parabolic $$p$$ -Laplace equation for every $p\in(1,\infty)$ . The proof combines a logarithmic Bernstein argument with quantitative stability to derive a recursive estimate between different regularizations. For $$p=1$$ and $n\geq2$ , we construct a bounded viscosity solution which is not Hölder continuous in time with any positive exponent.

14. 度量图上的 Lane–Emden 问题最小能量解的渐近行为

Asymptotic Behavior of Least Energy Solutions to the Lane–Emden Problem on Metric Graphs

作者: Kazuki Sato

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15486v1

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

摘要: In this paper, we study the asymptotic behavior of least energy solutions to Lane–Emden problems on compact metric graphs as $p \to \infty$ . For the problem with Dirichlet–Kirchhoff boundary conditions, we characterize the limiting variational problem in terms of the Dirichlet–Kirchhoff Green function and show that least energy positive solutions converge, up to a subsequence, to a normalized Green function centered at a maximizer of the diagonal Green function. We also show that their maximum points approach the set of such maximizers. For the problem with Kirchhoff–Neumann boundary conditions, we characterize the limiting variational problem on compact metric graphs without cycles and show that least energy solutions converge, up to a subsequence, to a limiting profile associated with a pair of points realizing the maximal distance on the graph. Moreover, the distance between their maximum and minimum points converges to the maximal distance on the graph.

15. 关于 Dirichlet 分数对数拉普拉斯算子的 Faber-Krahn 不等式

On Faber-Krahn inequality for Dirichlet Fractional-Logarithmic Laplacian

作者: Huyuan Chen, Rui Chen, Bobo Hua

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15481v1

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

摘要: We study the Faber–Krahn problem for the fractional–logarithmic Laplacian $(-Δ)^{s+\ln}$ , $$0<s<1$$ , with Fourier symbol $|ξ|^{2s}\ln|ξ|^2$ . We prove a small-volume and large-volume dichotomy. For sufficiently small volume, balls minimize the first Dirichlet eigenvalue, with rigidity up to translation.

16. 抛物型一致可整直性由有界热函数的 Carleson 测度估计刻画

Parabolic uniform rectifiability is characterized by Carleson measure estimates for bounded caloric functions

作者: Simon Bortz, Pablo Hidalgo-Palencia, Steve Hofmann, James Warta

类别: math.AP, math.CA

发布日期: 2026-09-14

arXiv ID: 2609.15374v1

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

摘要: We show that the parabolic uniform rectifiability of a set can be characterized by an interior PDE property, namely Carleson measure estimates for bounded solutions to the heat equation. In particular, under very mild background hypotheses on an open set $$Ω$$ and its boundary $\partial Ω$ , we show that the following are equivalent: (i) The quantity $|\nabla u(\cdot)|^2 \text{dist}(\cdot, \partialΩ)$ is the density of a Carleson measure on $$Ω$$ for all bounded solutions to the heat equation in $$Ω$$ . (ii) The caloric measure for $$Ω$$ admits a corona decomposition. (iii) $\partial Ω$ is parabolic uniformly rectifiable. The implication (iii) implies (i) is the main result of the cited paper [Bortz, Hoffman, Hofmann, Luna-García, Nyström, Anal. PDE 2023]. The remaining implications are new, and are direct parabolic analogues of results of Garnett, Mourgoglou and Tolsa [Duke Math. J. 2018] in the elliptic setting for the Laplacian; however, our proofs require several new ideas. In particular, we must adapt arguments to account both for time-lag, and for the fact that some characterizations of (elliptic) uniform rectifiability have either not been developed or have been shown to be untrue in the parabolic setting. Our background assumptions, as noted above, are very mild: we assume that the domain satisfies the time-symmetric capacity density condition and has interior corkscrews, and that the boundary is (non time-directed) Ahlfors-David regular. These are weaker than the background assumptions of Bortz, Hofmann, Martell and Nyström in arXiv:2510.22047, and are essentially optimal for the boundary continuity properties of caloric functions on which our arguments rely.

17. 奇异椭圆方程的径向节点 Dirichlet 解：全局分支、端点渐近与 Morse 指标

Radial Nodal Dirichlet Solutions of Singular Elliptic Equations: Global Branches, Endpoint Asymptotics, and Morse Indices

作者: Wenjing Chen, Zexi Wang

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15283v1

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

摘要: We establish sharp finite-ball shooting classifications for radial nodal Dirichlet solutions of the logarithmic equation and, within the maximal simple-zero shooting class, of the sublinear scalar-field equation. Recent whole-space uniqueness and phase-transition theorems are used as external inputs, while the finite-ball zero-curve ranges, endpoint asymptotics, and spectral consequences are proved here. The two models require different singular analyses: in the logarithmic problem the nonlinearity is not locally Lipschitz at a nodal zero and the linearized potential diverges there, whereas in the sublinear problem the limiting whole-space profile reaches a double zero at a finite support radius, beyond which continuation is nonunique. For every $$R>0$$ and $k\ge0$ , the logarithmic problem has, up to sign, a unique radial Dirichlet solution with exactly $$k$$ interior zeros. Its shooting height $β_k^{\log}(R)$ is a strictly decreasing $$C^1$$ bijection from $(0,\infty)$ onto $(α_k^{\log},\infty)$ and satisfies

β_k^{\log}(R)\longrightarrowα_k^{\log} \quad(R\to\infty),\qquad \log\bigl(β_k^{\log}(R)^2\bigr) =\frac{ρ_{k+1}^2}{R^2}+κ_{k+1,n}+o(1) \quad(R\downarrow0),

where $κ_{k+1,n}>0$ is given by an explicit Bessel integral. For the sublinear problem, the maximal simple-zero branch exists precisely for $R\in(ρ_{k+1},S_k)$ , where $ρ_{k+1}$ is the $$(k+1)$$ -st positive zero of the regular Bessel profile $$Φ_n$$ , and $$S_k$$ is the support radius of the unique compactly supported $$k$$ -node whole-space bound state.

18. Mather 测度与集合的全离散变分逼近

A Fully Discrete Variational Approximation of Mather Measures and Sets

作者: Fabio Camilli, Cristian Mendico

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

发布日期: 2026-09-14

arXiv ID: 2609.15281v1

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

摘要: We introduce a fully–discrete variational approximation of Mather measures and sets for Tonelli Lagrangians on the flat torus, together with a numerical procedure for approximating the entire Mather set. The scheme is based on a fully–discrete Lax–Oleinik operator with integer winding labels. We prove an $$O(τ+h/τ)$$ error estimate for the critical value, convergence of critical solutions, and a finite-dimensional characterization of fully–discrete Mather measures. Accumulation points of the reconstructed minimizing measures are continuous Mather measures, while the supports satisfy complementary upper and lower convergence results involving the Mañé and Mather sets. To avoid the selection of only some minimizing components by exact discrete minimizers, we introduce a mass-threshold approximation based on almost-minimizing holonomic measures. This yields a finite-dimensional constrained optimization procedure designed to recover the whole Mather set.

19. 色散谐振器阵列中的亚波长例外点

Subwavelength exceptional points in dispersive resonator arrays

作者: Konstantinos Alexopoulos

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

发布日期: 2026-09-14

arXiv ID: 2609.15200v1

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

摘要: Exceptional points provide a mechanism for producing strong spectral sensitivity in non-Hermitian wave systems. For subwavelength resonators with highly dispersive material parameters, however, the resonance problem is not a linear eigenvalue problem in the frequency. In this paper, we use the Lippmann-Schwinger formulation and the associated volume-potential reduction to study exceptional points in dispersive resonator arrays. The reduction leads to a finite-dimensional nonlinear matrix pencil whose entries encode both the local material response and the volume-mediated interaction between the resonators. We characterize exceptional points as defective characteristic values of this pencil and derive explicit conditions for dimers and higher-order systems. We show that passive damping and detuning can generate exceptional points without imposed gain, and we formulate an inverse-design condition for placing such degeneracies at prescribed frequencies. We also prove the corresponding Puiseux splitting laws and identify how material dispersion enters the leading sensitivity prefactor. Numerical experiments, including halide-perovskite-inspired dimers and trimers, illustrate the theory and show that the exceptional-point mechanism persists for explicit highly dispersive passive material laws.

20. 关于非线性离散 Calderón 问题唯一性

On the Uniqueness of a Nonlinear Discrete Calderón Problem

作者: Elena Beretta, Maolin Deng, Alberto Gandolfi, Bangti Jin

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15175v1

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

摘要: The discrete Calderón problem aims at recovering the conductivity on the edges of a graph from boundary measurements, encoded by the discrete Dirichlet-to-Neumann (DtN) map. The problem has been intensively studied in the linear case on square lattices since the seminal works of Curtis and Morrow. In this work, we investigate a nonlinear analogue of the discrete Calderón problem for a semilinear second-order elliptic equation on square lattices. We establish three uniqueness results for the conductivity recovery. First, we show that conductivity-dependent corner excitations allow a layer-by-layer reconstruction of the conductivity. Second, we study the linearization of the nonlinear DtN map at an arbitrary background boundary datum and prove that the conductivity and the background potential are uniquely determined by one pair of nonlinear Cauchy data and the linearized DtN map. Third, we show that the linearized data can be replaced by finitely many nonlinear measurements, which uniquely determine the conductivity.

21. 球域内任意给定质量下的NLS方程归一化解

Normalized solutions to the NLS equation in the ball for any prescribed mass

作者: Nicola Soave

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15141v1

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

摘要: Given $$ρ>0$$ , we consider the problem

\text{find $(λ,u) \in \mathbb{R} \times H_0^1(B)$ such that } \begin{cases} -Δu+λu = |u|^{p-1}u & \text{in } B \\ \int_B u^2\,dx = ρ, \end{cases}

where $$B$$ is a ball in $\mathbb{R}^N$ , $N \ge 1$ , and $$1<p<2^*-1$$ . Without any further restriction on $$N$$ , $$ρ$$ and $$p$$ , we prove the existence of infinitely many radial solutions, and, in dimension $N \ge 4$ , of at least one non-radial solution.

22. 具有强阻尼和延迟的非自治板方程系统的拉回吸引子

Pullback Attractors for a Non-Autonomous Plate System with Strong Damping and Delay

作者: Xuxiao Hou, Bin Yang, Yuming Qin, Alain Miranville

类别: math.AP, math.DS

发布日期: 2026-09-14

arXiv ID: 2609.15110v1

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

摘要: This paper addresses the existence of pullback attractors for a class of nonlinear non-autonomous strongly damped plate equations with delay. The model contains the biharmonic operator, the strong damping term, a nonlinear source term and a delay operator acting on the history of the solution. The interaction between the strong damping mechanism and the hereditary delay effect gives rise to several analytical difficulties in deriving uniform estimates and proving pullback asymptotic compactness. To overcome these difficulties, we construct a modified energy functional adapted to the strong damping structure and employ the contractive function method to handle the delay term in the history phase space. By establishing the existence and uniqueness of weak solutions, uniform pullback estimates and the pullback asymptotic compactness of the associated non-autonomous process, we prove the existence of pullback attractors for the strongly damped plate equation with delay.

23. 裸奇点标量场的内部引力扰动

Interior gravitational perturbations to naked singularities of a scalar field

作者: Junbin Li, Tingting Li

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

发布日期: 2026-09-14

arXiv ID: 2609.15107v1

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

摘要: For the $$k$$ -self-similar naked singularity solutions of Einstein–scalar field equations constructed by Christodoulou, we construct a family of interior gravitational perturbations leading to trapped surface formation that remain large but finite at the threshold Hölder regularity, while converging to zero in all regularities below the threshold. This extends the interior spherically symmetric result below the threshold in \cite{Li25} to non-spherically symmetric setting. We further construct a genuinely localized family of perturbations whose angular support shrinks to a single point, with smooth convergence to zero away from that point. This exploits the additional angular freedom available beyond spherical symmetry. Moreover, when measured instead in Sobolev regularity, such genuinely localized family of perturbations still converges to zero in all regularities below the same threshold.

24. $\mathbb{S}^n$ 中两相过定问题的刚性

Rigidity for two-phase overdetermined problems in $\mathbb{S}^n$

作者: Andrea Bisterzo, Shigeru Sakaguchi

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.15103v1

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

摘要: In this paper, we study rigidity phenomena associated with a class of overdetermined two-phase problems on the $$n$$ -dimensional sphere $\mathbb{S}^n$ . Specifically, we analyze the solutions of coupled elliptic equations defined on a domain $Ω\subset \mathbb{S}^n$ and its complement, subject to Dirichlet boundary conditions on the interface $\partialΩ$ and a compatibility condition on the gradient. We utilize two different methods. The method of moving planes works when $\partialΩ$ is contained inside an open hemisphere. The other method works when $$Ω$$ is simply connected in $\mathbb{S}^2$ . Both lead to the same conclusion. Namely, the existence of a solution to such overdetermined problems necessarily implies that $$Ω$$ is a geodesic ball. In the latter, we prove this rigidity property for three distinct scenarios: a baseline problem with piecewise constant source terms, which is the two-phase version of the torsion problem; a generalization including a class of positive and regular nonlinearities; and finally an eigenvalue problem involving the first Dirichlet eigenvalue of the respective phases.

25. 具有强阻尼和延迟的 Kirchhoff 波动方程拉回 $\mathcal{D}$ -吸引子的存在性

Existence of pullback $\mathcal{D}$ -attractors for Kirchhoff wave equations with strong damping and delay

作者: Bin Yang, Yuming Qin, Alain Miranville, Xuxiao Hou

类别: math.AP, math.DS

发布日期: 2026-09-14

arXiv ID: 2609.15069v1

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

摘要: This paper studies the existence of pullback $\mathcal{D}$ -attractors for non-autonomous Kirchhoff wave equations with strong damping and delay effects in the phase space $\mathcal{E} = C_{H^1_0(Ω)} \times C_{L^2(Ω)}$ . The model contains a strong damping term $-Δ\partial_t u$ , a nonlocal Kirchhoff term $Ψ(\|\nabla u\|^2)$ , a state-dependent delay term $$φ(t, u_t)$$ , and a time-dependent external force $$h(x, t)$$ . There appear to be no results on the pullback attractors of Kirchhoff wave equations involving both strong damping and delay effects in the literature. To this end, we establish delicate uniform estimates and employ the contraction function method to prove the pullback asymptotic compactness of the associated process, which yields the existence of pullback $\mathcal{D}$ -attractors.

26. Bernoulli 单相问题中 Perron 极值解的变分性质

Variational properties of Perron’s extremal solutions in the Bernoulli one-phase problem

作者: Farhan Abedin, William M. Feldman, Kerrek Stinson

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.14981v1

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

摘要: We study the Perron extremal solutions of the stationary one-phase Bernoulli problem. These solutions are important in several applications, but not much is known about the regularity of their free boundaries due to their non-variational construction. We show that, in fact, Perron extremal solutions retain some variational features; specifically, they are solutions in the sense of inner variations. We establish this property by showing that extremal solutions arise as infinite-time limits of monotone solutions of the parabolic Bernoulli problem. As a consequence, we are able to describe the structure of singularities for Perron extremal solutions in two dimensions: largest subsolutions have smooth free boundaries, while smallest supersolutions have free boundary points that are either smooth, or whose blow-up limits are two-plane wedges with equal slope.

27. 通过耦合上下解的异质 FitzHugh–Nagumo 系统中的传播失效

Propagation failure in heterogeneous FitzHugh–Nagumo systems via coupled upper and lower solutions

作者: Matías Courdurier, Esteban Paduro

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.14944v1

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

摘要: Traveling waves are fundamental objects in reaction–diffusion systems; however, spatial heterogeneities can distort or inhibit their propagation. Determining whether a given heterogeneity results in propagation failure remains challenging, particularly in systems with coupled components. In this paper, we present explicit constructions of stationary coupled upper and lower solutions tailored to specific heterogeneities for a one-dimensional FitzHugh–Nagumo system. By leveraging the mixed quasimonotone structure of the system, these barriers yield sufficient conditions for the failure of propagation for families of trapped initial conditions. The stationary constructions establish directional signal blocking in FitzHugh–Nagumo neuron models with localized depressed dynamics or abrupt geometric variations, and the resulting stationary barriers persist for sufficiently small recovery diffusion. Additionally, we introduce explicit almost stationary barriers, which provide conditional bounds on the expansion of componentwise threshold regions in homogeneous media.

28. 基于预训练非线性层的线性化 PINN

Linearized PINN with pretrained nonlinear layers

作者: Wenhao Chen, Alexandre M. Tartakovsky

类别: math.AP, cs.LG

发布日期: 2026-09-14

arXiv ID: 2609.14926v1

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

摘要: We propose a linearized Physics-Informed Neural Network (lPINN), a reduced-order neural basis method for forward and inverse differential equations. In an offline stage, lPINN learns operator-compatible continuous neural basis functions from an ensemble of numerical solutions. The basis functions are differentiable through automatic differentiation and are pretrained using solution data together with either derivative information or physics residuals. For each new problem instance, the basis functions are frozen and the solution is obtained by minimizing the governing-equation residual together with applicable initial, boundary, regularization, and observational terms. Unlike surrogate and operator-learning methods, the training data define the trial space offline, while the instance-specific solution is computed online by enforcing the governing physics. Relative to vanilla PINNs, lPINN pretrains the nonlinear hidden-layer representation offline and performs online inference only in the final linear layer. We evaluate lPINN on forward and inverse problems for the advection-diffusion equation, Burgers’ equation, and the nonlinear pendulum equation. Compared with vanilla PINNs, lPINN achieves lower solution and parameter errors while reducing online inference times by approximately one to more than three orders of magnitude, with the largest gains generally observed for limited residual or measurement data. Cross-resolution experiments show that the learned continuous representation can be evaluated on finer meshes without retraining and with nearly unchanged accuracy.

29. 五专家预测偏微分方程的显式解与COMB的精确最优集

An explicit solution of the five-expert prediction PDE and the exact optimality set of COMB

作者: Jeff Calder, Nadejda Drenska

类别: math.AP, cs.GT, cs.LG, math.OC

发布日期: 2026-09-14

arXiv ID: 2609.14892v1

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

摘要: In this paper, we derive an explicit solution of the stationary prediction with expert advice PDE for five experts. The formula is given in three regions. In the first two regions, it is the four-expert solution plus a single integral with an elementary positive density. In the third region, it is a finite sum of hyperbolic products whose coefficients are determined by one scalar quadrature. Our formula establishes that the direction $$(1,0,1,0,0)$$ is optimal throughout the ordered sector, and that the COMB strategy $$(1,0,1,0,1)$$ is optimal only on a lower dimensional subset of the sector (where $$x_1=x_2$$ and $$x_3=x_4$$ ). This disproves the COMB optimality conjecture of Gravin, Peres and Sivan (2016). The verification of the Hamiltonian inequalities is a tedious task, part of which is completed with a computer assisted proof. The verification reduces to 21 scalar inequalities, which we prove using 147 exact rational Bernstein polynomial certificates. The exact certificates and their independent arithmetic checks are included in a supplement to this paper.

30. 时间分数阶 Amari 模型中收敛于渐近行波的全局前缘解

Global front solutions converging to an asymptotic traveling wave in a time-fractional Amari model

作者: Hiroshi Ishii

类别: math.AP

发布日期: 2026-09-14

arXiv ID: 2609.14884v2

链接: http://arxiv.org/abs/2609.14884v2

摘要: We study front propagation in a time-fractional Amari neural field equation with a Caputo derivative and a Heaviside firing-rate function. For a class of nonnegative even connectivity kernels and strictly decreasing $$C^1$$ front-like initial data, we reduce the Cauchy problem to a nonlinear Volterra equation for the threshold location. We prove the existence and uniqueness of a global front solution, determine the short-time behavior of its threshold location, and identify its asymptotic speed. We also construct the associated asymptotic traveling-wave profile and establish uniqueness of the speed and profile, up to translation, in the nonzero-speed case. Using a rigidity argument, we prove the asymptotic linearity of the local increments of the threshold location and obtain uniform convergence to the profile in the threshold-centered frame. Finally, we characterize the kernel-dependent tail asymptotics of the wave profile.

31. Robin特征值不定权下小有利区域的最优位置

Optimal location of small favourable regions for Robin eigenvalues with indefinite weights

作者: Baruch Schneider, Diana Schneiderová, Yifan Zhang

类别: math.AP, math.SP

发布日期: 2026-09-13

arXiv ID: 2609.14763v1

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

摘要: We consider the positive principal eigenvalue of an elliptic problem with a Robin boundary condition and bang–bang indefinite weight $κ\mathbf 1_D-\mathbf 1_{Ω\setminus D}$ , and ask where a favourable region $$D$$ of prescribed small volume $$|D|=δ$$ should be located. Put $\varepsilon=δ^{1/N}$ and $τ_δ=α_δ\varepsilon$ . We prove that there is a finite threshold $$τ_*=τ_*(N,κ)$$ , independent of the ambient domain, such that $δ^{2/N}Λ_δ(α_δ)\toΛ_{\mathbb H}(τ)$ whenever $τ_δ\toτ<\infty$ . If $$τ<τ_*$$ , optimal small regions concentrate at the boundary. If $$τ>τ_*$$ , their concentration centres move to distances much larger than $\varepsilon$ from the boundary, and after recentring and rescaling the favourable sets converge in measure to the whole-space optimal ball. At $$τ=τ_*$$ , the half-space problem admits a compact boundary optimiser as well as minimising sequences escaping to infinity. In the finer regime $τ_δ=τ_*+σ\varepsilon+o(\varepsilon)$ , we determine the first-order competition between the boundary and interior configurations. Tangential symmetry of compact threshold optimisers reduces the geometric correction to mean curvature. In particular, every fixed finite Robin coefficient is asymptotically in the boundary regime. Numerical experiments illustrate the boundary–interior transition, the three cases in the first-order selection law, and the curvature-dependent boundary location; for $$N=2$$ and $$κ=1$$ they suggest a transition near $τ_*\approx3.2$ .

32. 迹为常数的非常数弱调和映射

Nonconstant weakly harmonic maps with constant trace

作者: David Ziener

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14760v1

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

摘要: We construct a nonconstant weakly harmonic map from a three-dimensional ball into the two-sphere with constant trace. This gives an affirmative answer to a question raised by Rivière in 1995.

33. 零子簇的有限性与 J-方程的最优失稳项

Finiteness of null subvarieties and optimal destabilizers for the J-equation

作者: Junbang Liu

类别: math.DG, math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14713v2

链接: http://arxiv.org/abs/2609.14713v2

摘要: Let $$X$$ be a compact Kähler manifold and $$(α,β)$$ be a pair of Kähler classes. We show that when $$(α,β)$$ is $$J$$ -semistable, the collection of positive-dimensional irreducible $$J$$ -null subvarieties is finite. We also obtain a uniform positive lower bound for the normalized $$J$$ -slope of every irreducible non- $$J$$ -null subvariety. This confirms the finiteness and uniform-gap expectation of Khalid-Sjöström Dyrefelt \cite[Section 1.2, Theorem 1.5; Remark 1.8]{KSD26}. As a corollary, we show that the collection of optimal destabilizers is finite when the stability threshold is nonpositive.

34. $C^{1,1}$ 正则性全局估计：带洞域上的齐次复 Hessian 方程

$C^{1,1}$ Regularity Global Estimates for Homogeneous Complex Hessian Equations on Punctured Domains

作者: Zhenghuan Gao, Xi-Nan Ma, Bao Yu, Dekai Zhang

类别: math.CV, math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14692v1

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

摘要: In this paper, we prove the $C^{1,1}$ regularity of Green functions associated with the complex $$k$$ -Hessian operator for $1\leq k<n$ , and give a new proof of the corresponding regularity of pluricomplex Green functions. We establish real Hessian estimates for approximating problems associated with homogeneous complex Hessian equations on punctured domains. The main idea is to introduce new auxiliary functions generated by complex linear vector fields to reduce the global real Hessian estimate to the boundary. We also treat the complex Monge-Ampère case. Similar real Hessian estimates also works for the exterior problem.

35. 弱竞争情形下气候变化背景下入侵物种与本土物种的自由边界模型

A free boundary model for invasive and native species under shifting climate in the weak competition case

作者: Phuoc Vinh Dinh, Phuong Le, Tien Dung Nguyen

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14679v1

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

摘要: We study a free boundary problem for a diffusive Lotka–Volterra competition system describing the invasion of a new species into the habitat of a native competitor, in a habitat that is shifted from unfavourable to favourable at a constant speed $$c>0$$ by climate change. Only the invader feels the shifting environment and only its range is governed by a Stefan-type free boundary, while the native species occupies the whole half line. We work throughout in the weak competition regime, in which the two species may coexist. We prove a spreading–vanishing dichotomy: either the invader spreads and the pair converges to the coexistence steady state $$(u^*,v^*)$$ , or the invader vanishes and the native species recovers its carrying capacity. In the vanishing case we obtain the explicit bound $\lim_{t\to\infty}h(t)\le\fracπ{2}\sqrt{d_1c_2/(a_1c_2-a_2c_1)}$ , and we give criteria guaranteeing each alternative. When spreading occurs, we determine the exact asymptotic spreading speed: $\lim_{t\to\infty}h(t)/t=\min\{c,c_0\}$ , where $$c_0$$ is the spreading speed of the corresponding homogeneous weak competition system. In particular the invasion is slowed down both by the competitor and by the climate shift, and the slower of the two mechanisms is the one that determines the speed. Numerical simulations illustrate the results.

36. 分数阶拉普拉斯算子的定量 Gidas–Ni–Nirenberg 定理

A Quantitative Gidas–Ni–Nirenberg Theorem for the Fractional Laplacian

作者: Qizhen Shen, Jiongduo Xie

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14618v1

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

摘要: We establish a quantitative Gidas–Ni–Nirenberg theorem for the pure fractional Dirichlet problem

(-Δ)^s u=k(x)g(u)\quad\text{in }B_1, \qquad u=0\quad\text{in }\mathbb R^n\setminus B_1, \qquad 0<s<1.

For positive bounded solutions under a two-sided $L^\infty$ normalization, we prove

\mathcal D(u)\le C\,\mathcal D(k)^γ,\qquad 0<γ<1,

where $\mathcal D$ measures spherical oscillation and outward radial increase. To the best of our knowledge, this provides the first quantitative Gidas–Ni–Nirenberg estimate for the pure fractional Dirichlet problem. The key new ingredient is a shift-uniform, weak Harnack inequality for shifted antisymmetric supersolutions, which overcomes the exterior-tail obstruction in the quantitative moving-plane argument. This mechanism extends to nonseparable nonlinearities.

37. 全局分析：黎曼流形上非线性分析及其变分方法导论

Global Analysis: An Introduction to Nonlinear Analysis and Its Variational Methods on Riemannian Manifolds

作者: Carlos Daniel Velázquez-Mendoza, María de los Ángeles Sandoval-Romero, Romulo Diaz Carlos

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

发布日期: 2026-09-13

arXiv ID: 2609.14580v1

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

摘要: This monograph develops an introduction to global analysis centered on the interaction between differential geometry, functional analysis, partial differential equations, and variational methods on Riemannian manifolds. Beginning with smooth and Riemannian geometry, it develops Sobolev spaces, distributions, interpolation and fractional regularity, differential and pseudodifferential operators on vector bundles, elliptic theory, heat methods, bounded geometry, and trace theorems. It then treats Fredholm and index theory, culminating in the Atiyah–Singer index theorem, followed by geometric evolution equations and Ricci flow, infinite-dimensional geometry on Banach and Hilbert manifolds, and variational methods including the direct method, Palais–Smale theory, deformation arguments, the mountain pass theorem, and the Nehari method. Particular emphasis is placed on explicit proofs, the passage from local Euclidean estimates to intrinsic global statements, and the precise geometric hypotheses required in compact, noncompact, and boundary settings. The text is intended for advanced undergraduate and graduate students, as well as readers approaching global analysis from geometry or differential equations.

38. 气候变化背景下本土竞争者领地内入侵物种扩散的自由边界问题

A free boundary problem for spreading of an invasive species in the territory of a native competitor under shifting climate

作者: Phuoc Vinh Dinh, Phuong Le, Tien Dung Nguyen

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14573v1

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

摘要: In this paper we consider a free boundary problem for the spreading of an invasive species in the habitat of a native competitor. We assume that the invasive species benefits from a climate shift that turns the environment from unfavorable to favorable at a constant speed $$c$$ . We show that if the invasive species is an inferior one, then it must vanish, while the native competitor always persists. However, if the invasive species is a superior one, then a dichotomy occurs: either the invasive species vanishes and the native one persists, or the invasive one spreads and the native one vanishes. We also show that in the latter case the asymptotic spreading speed equals $\min\{c, c_0\}$ , where $$c_0$$ is the spreading speed in the corresponding homogeneous environment. Numerical simulations are provided to illustrate our theoretical results.

39. 非线性费米子系统的基态：从幂律相互作用到对数Sobolev不等式

Ground States of Nonlinear Fermionic Systems: From Power-Law Interactions to Logarithmic Sobolev Inequality

作者: Bin Chen, Yujin Guo, Yong Luo, Zhenya Yan

类别: math.AP, math-ph

发布日期: 2026-09-13

arXiv ID: 2609.14568v1

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

摘要: We consider ground states of a two-component logarithmic fermionic system in $\mathbb{R}^d$ , where $d\ge 1$ is arbitrary. We prove that up to translations and scalings, ground states of the logarithmic system are the $L^\infty$ -limits of ground states for a two-component $$2p-1$$ power-law fermionic system as $p\searrow 1$ . As a byproduct, we also establish a sharp logarithmic Sobolev inequality for orthonormal functions in $\mathbb{R}^d$ , whose optimizers are, up to scalings, the minimizers of a constraint variational problem associated with the logarithmic system.

40. 星形域的全等参数化与半线性椭圆方程解的持久性

Conformal Parametrisation of Star-Shaped Domains and Persistence of Semilinear Elliptic Solutions

作者: Gianni Arioli

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14546v1

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

摘要: We study the conformal parametrisation of smooth star-shaped planar domains and its dependence on perturbations of the boundary. The boundary correspondence of the normalized Riemann map is described by the Theodorsen equation, which we formulate in weighted Wiener algebras. We derive an explicit representation of the inverse of the linearized Theodorsen operator through a canonical Wiener–Hopf factorisation and obtain quantitative a posteriori estimates that are stable under perturbations of the boundary function. For a collection of nonperturbative reference domains, rigorous computer-assisted estimates yield explicit neighborhoods in a weighted Wiener algebra such that every boundary function in these neighborhoods admits a uniquely determined conformal parametrisation in the validated ball, with quantitative control of the corresponding Riemann map. Thus the validation applies to open, infinite-dimensional families of star-shaped domains. As an application, we consider the semilinear Dirichlet problem

-Δv=v^3\quad\text{in }Ω, \qquad v=0\quad\text{on }\partialΩ.

Pulling the equation back by the validated conformal maps yields a problem on the unit disk. Uniform a posteriori estimates then imply persistence of the elliptic solutions throughout explicit neighborhoods of the reference domains.

41. 二维中周长与其分数阶类似量的精确比较

Sharp comparison between the perimeter and its fractional analogue in two dimensions

作者: Rupert L. Frank, Paata Ivanisvili

类别: math.AP, math.DG

发布日期: 2026-09-13

arXiv ID: 2609.14513v1

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

摘要: We show that the fractional perimeter of any order $$0<s<1$$ is maximal, among all planar sets of fixed perimeter, for the disk.

42. Kirchhoff 方程的台球、折射与爆破

Billiards, Refraction, and Blow-Up for Kirchhoff Equations

作者: Marina Ghisi, Massimo Gobbino

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14502v1

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

摘要: We start from a simple problem in geometric optics. A particle moves between two homothetic ellipses, with refraction at the inner interface and reflection at the outer one. In a suitable nearly circular regime, the small geometric anisotropy of the ellipses is amplified by the strong refraction, and the corresponding return map develops a transverse heteroclinic connection between the two axial motions. We turn this geometric mechanism into a construction for a two-mode Kirchhoff system while keeping the standard quadratic elastic variable throughout. The singular optical model is first reduced to an explicit kick–drift map, for which the heteroclinic connection is obtained by a contraction argument. The same argument also gives exponential convergence along the two tails and transversality. We then show that the connection persists for the genuine elliptic billiard and, subsequently, through a smooth regularization of the reflecting and refracting interfaces. Finally, a small positive background is added to the Kirchhoff coefficient, making it uniformly positive without destroying the hyperbolicity of the axial modes or the transverse connection. This produces a smooth and uniformly positive Kirchhoff coefficient admitting a transverse heteroclinic orbit between two simple modes. Combined with the road-map theorem developed in our previous work, the construction yields a finite-time blow-up example for a forced abstract Kirchhoff equation with a regular forcing term.

43. 任意维度下具有弱奇异敏感性和非局部项的趋化系统的整体有界性

Global boundedness of the chemotaxis system with weakly singular sensitivity and nonlocal term in any dimension

作者: Yang Cao, Qingchun Li, Jing Zhang

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14443v1

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

摘要: This paper is concerned with the parabolic-elliptic chemotaxis system involving weakly singular sensitivity and a nonlocal source: $u_t=Δu-χ\nabla\cdot\left(\frac{u}{v^k}\nabla v\right) +u^α\left(r-μ\int_Ωu^β\,dx\right)$ and $$0=Δv-v+u^γ$$ under homogeneous Neumann boundary conditions in a smooth bounded domain $Ω\subset\mathbb{R}^N$ with $N\ge1$ , where $$χ,r,μ,γ>0$$ , $k\in(0,1)$ and $α,β\ge1$ . In view of the enhanced aggregation induced by singular sensitivity and the suppression of excessive population growth by the nonlocal damping, we identify sufficient conditions under which solutions remain globally bounded. We show that classical solutions are globally bounded in the following cases: \par\smallskip \begin{center} \small \setlength{\tabcolsep}{5pt} \setlength{\arrayrulewidth}{0.3pt} \renewcommand{\arraystretch}{1.25} \begin{tabular}{c|c|c|c} \multirow{2}{}{ $$β>1$$ } & \multirow{2}{}{ $$γ=1$$ } & $$N=1$$ & $1\leqα<1+2β$ \ \cline{3-4} & & $N\geq2$ & $\begin{gathered} 1\leqα<2,\quad α+β>2+\tfrac{N}{2},\\[-0.5mm] \text{or}\quad 2\leqα<1+\tfrac{2β}{N} \end{gathered}$ \ \hline \multirow{2}{}{ $$β=1$$ } & $$γ=1$$ & $$N=1$$ & \multirow{2}{}{ $\begin{array}{l@{\quad}l@{\quad}l} \text{Case 1:} & 1\leqα<1+\tfrac{2}{N} & \text{if}\quad m_0<\tfrac{r}μ,\\[4pt] \text{Case 2:} & α\geq1 & \text{if}\quad m_0\geq\tfrac{r}μ \end{array}$ } \ \cline{2-3} & $0<γ<\tfrac{2}{N}$ & $N\geq2$ & \ \end{tabular} \end{center} \par \smallskip \noindent Here $m_0:=\int_Ωu_0$ denotes the initial total mass. It is worth noting that both the uniform-in-time $$L^p$$ -estimate and the $L^\infty$ -bound obtained by Moser iteration are derived without first establishing a uniform positive lower bound for $$v$$ .

44. 关于 $$oldsymbol{ abla}$$ 在 $oldsymbol{ abla}^2$ 中紧支光滑数据的光滑解存在性问题的一个注记

A note on existence of smooth solution to Jacobian equation for compactly supported smooth data in $\mathbb{R}^2$

作者: Animesh Jana

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14407v2

链接: http://arxiv.org/abs/2609.14407v2

摘要: In this note, we give an explicit construction of global solutions to the prescribed Jacobian equation

\det(\nabla u)=f\mbox{ in }\mathbb{R}^2,

for a class of data. For every $f\in C_c^1(\mathbb{R}^2)$ and $$p>1$$ , we construct a solution $u\in \dot W^{1,2p}(\mathbb{R}^2)\cap L^\infty(\mathbb{R}^2)$ . In particular, no sign condition or integral constraint is imposed on the compactly supported data. For $f\in C_c^\infty(\mathbb{R}^2)$ , the construction yields a smooth solution. We also consider rapidly decaying data and prove the existence of global solutions with bounded gradient for every $f\in\mathcal{S}(\mathbb{R}^2)$ . Finally, we extend the construction to compactly supported data that are measurable in one variable and $$C^1$$ in the other. Our proof relies on a similar idea as in [Moser, Trans. Amer. Math. Soc. 1965].

45. 光线变换的(不)稳定性结果

(In)stability results for the Light Ray Transform

作者: Leonard Busch, Sebastián Muñoz-Thon, Lauri Oksanen

类别: math.AP, math.DG

发布日期: 2026-09-13

arXiv ID: 2609.14387v1

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

摘要: We prove quantitative instability results for the light ray transform in the Minkowski setting, and locally in the smooth and Gevrey categories. In the smooth Sobolev setting, we show that no modulus of continuity can be better than $$t^α$$ for any $α\in(0,1)$ , while in the Gevrey setting we obtain an explicit logarithmic lower bound. On the other hand, using the relation between stationary geometries and magnetic(-potential) systems, we obtain a stability result for “time moments” of the light ray transform.

46. 初始不连续曲面上非线性扩散与平均曲率的相互作用

Interaction between nonlinear diffusion and mean curvature of a surface with initial discontinuities

作者: Shigeru Sakaguchi

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14358v1

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

摘要: We consider the Cauchy problems for a class of nonlinear diffusion equations where each initial data is given by a characteristic function of an open set $$Ω$$ in the whole Euclidean space. If the boundary $\partialΩ$ is of class of $$C^2$$ in a neighborhood of a point $$x$$ on it, then we extract the mean curvature of $\partialΩ$ at $$x$$ from the initial diffusion at $$x$$ . This corresponds to the author’s previous result for the two-phase linear diffusion equation. Applications to some overdetermined problems for nonlinear diffusion equations are given.

47. 梯度依赖型椭圆方程的Lewy型非退化性与凸环上的全局梯度微分同胚

Lewy-Type Nondegeneracy for Gradient-Dependent Elliptic Equations and Global Gradient Diffeomorphisms on Convex Rings

作者: Da Liu, Shujun Shi

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14280v1

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

摘要: We prove two complementary Lewy-type theorems for elliptic equations whose coefficients depend only on the gradient. First, let $Ω\subset\mathbb{R}^3$ and let $$u$$ solve

a^{ij}(Du)u_{ij}=0,

where $$a$$ is a smooth, symmetric, positive definite matrix field on an open set containing $$Du(Ω)$$ . We show that if the gradient map $$Du$$ is a local homeomorphism, then $\det D^2u$ never vanishes; hence $$Du$$ is a local $C^\infty$ -diffeomorphism. Second, in every dimension, we consider capacitary solutions on convex rings. If the solution has no critical points and its level hypersurfaces are strictly convex, then ellipticity alone forces the Hessian to have one positive and $$n-1$$ negative eigenvalues, that is, inertia $$(1,n-1)$$ . Consequently, the gradient is a global diffeomorphism onto a radially parametrized ring in gradient space. Both results apply to the $$p$$ -Laplace and minimal surface equations. For the global minimal-surface result, existence of a smooth solution is assumed. Classical convex-ring results supply the noncriticality and strict level-set convexity needed in the global corollaries. Since the two coefficient matrices are real analytic on the relevant gradient ranges, the corresponding local and global gradient diffeomorphisms are real analytic.

48. 用于半线性抛物方程 Dirichlet 边界控制的直接和间接物理信息神经网络：条件误差分析

Direct and Indirect Physics-Informed Neural Networks for Dirichlet Boundary Control of Semilinear Parabolic Equations: A Conditional Error Analysis

作者: Nguyen Thanh Quang, Ta Thi Thanh Mai, Bui Xuan Dieu

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14269v1

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

摘要: We study physics-informed neural networks (PINNs) for the Dirichlet boundary control of a semilinear parabolic equation with Tikhonov regularization. Two approaches are considered. A direct PINN parameterizes the state and control by separate networks and minimizes a penalized form of the tracking objective. An indirect PINN instead represents the state, adjoint, and control by unconstrained networks trained jointly to satisfy the first-order optimality system, with the state-control coupling and the homogeneous adjoint boundary and terminal conditions imposed as soft penalty terms rather than enforced architecturally. For the indirect formulation we develop an error estimation framework that decomposes the total error into approximation, optimization, quadrature, and soft boundary/terminal-constraint contributions. Under standing assumptions on optimal-solution regularity and compatibility, network approximability, uniform Hölder control of the soft-constraint residuals, and a local neighborhood of the reference optimality-system solution, we derive a quantitative linearized stability estimate and a conditional local nonlinear residual-to-error estimate, and construct a computable residual indicator with a conditional reliability bound. Numerical experiments on two manufactured test problems - a cubic reactiondiffusion equation and a linear equation with a nontrivial boundary control and adjoint - illustrate the behavior of the direct and indirect formulations.

49. 关于结构化弹性密度与初始状态的被动恢复

On passive recovery of structured elastic density and initial states

作者: Yixian Gao, Hongyu Liu, Yang Liu

类别: math.AP

发布日期: 2026-09-13

arXiv ID: 2609.14215v1

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

摘要: We study simultaneous recovery of the (variable) mass density, initial displacement, and initial velocity for the three-dimensional isotropic elastic wave equation with known constant Lamé parameters. The data are the complete displacement trace on an enclosing boundary. We first assume that the density-weighted initial displacement and velocity have fixed known profiles in one spatial direction. The $$s^0$$ and $$s^1$$ coefficients of the zero-frequency Laplace expansion identify these weighted states. The $$s^2$$ and $$s^3$$ coefficients then give static Lamé orthogonality identities for the density difference. The exact difference expansion through order $$s^3$$ has an $$O(s^4)$$ remainder, uniformly on bounded spatial sets and bounded density-contrast and weighted-state classes. Under alignment of the two initial states and a nonzero moment of the density-weighted initial velocity, the two density identities reduce to a constant-vector static transform. Two opposite elastic null phases give uniqueness for one or two fixed vertical density profiles when the profile family is two-sided Laplace nondegenerate. When $λ+μ\ne0$ , each nonzero normal root has partial multiplicities $$2$$ and $$1$$ and admits a length-two Jordan chain. The associated polynomial–exponential Lamé mode produces derivatives of the bilateral profile transforms. The resulting Hermite–Laplace system gives uniqueness for aligned classes with up to four fixed vertical profiles and independent horizontal coefficients, provided that the profile system is Hermite–Laplace nondegenerate. Distinct translations of one compactly supported profile provide an explicit four-profile class. We also prove rigidity of the alignment reduction and exhibit an infinite-dimensional kernel for the reduced transform on unrestricted densities.

50. 凸环中关于高斯 $$p$$ -容量势的拟凹性的反例

A counterexample to quasiconcavity for the Gaussian $$p$$ -Capacity potentials in convex rings

作者: Lei Qin

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.14194v1

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

摘要: We construct a smooth, strictly convex ring that admit a unique solution of Gaussian $$p$$ -Capacity problem whose super-level sets are nonconvex. This result provide a negative answer.

51. 一个与最优绝缘相关的非线性问题

A nonlinear problem related to optimal insulation

作者: Genival da Silva

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.14152v1

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

摘要: We study a generalization of the classical optimal insulation problem by replacing the standard Laplacian with the \ti{p-Laplacian}, leading to a \ti{p-Poisson equation} with Robin boundary conditions. We also investigate the corresponding eigenvalue problem. This work extends and complements some results from \cite{bu1,pietra1,pietra2} to the nonlinear setting.

52. 自由边界调和映射的抛物型缺陷测度

Parabolic defect measures for Harmonic Maps with Free Boundary

作者: Anirban Das, Ali Hyder, Yannick Sire

类别: math.AP, math.DG

发布日期: 2026-09-12

arXiv ID: 2609.14050v1

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

摘要: We introduce defect measures for the heat flow of a class of harmonic maps with free boundary and initiate their analysis. Our flow is motivated by singular perturbations on the boundary of the domain, emphasizing structural differences between harmonic maps with free boundary and their classical counterpart. The harmonic maps we consider are a reformulation of maps introduced by Da Lio and Rivière and arise in several problems in the spectral geometry of Steklov eigenvalues. Contrary to the classical case, the flow we consider exhibits features of global nature and the blow-up analysis turns out to be much more challenging.

53. 随机抛物方程中漂移与扩散点源逆恢复的稳定性估计

Stability Estimates for the Inverse Recovery of Drift and Diffusion Point Sources in Stochastic Parabolic Equations

作者: Yu Wang, Qi Lü

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.14009v1

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

摘要: This paper studies the simultaneous recovery of point-source locations and temporal strengths in both the drift and diffusion terms of a stochastic parabolic equation with Neumann boundary conditions, using only boundary observations on a nonempty portion of the boundary after the sources have become inactive. By separating the stochastic equation into its expectation and centered parts, the drift channel reduces to a deterministic source problem with scalar Laplace moments, while the diffusion channel is controlled by Itô’s isometry, which preserves the full $$L^2$$ -energy of the integrand. Under suitable nondegeneracy and separation assumptions, we establish Lipschitz stability for the location of a single source in either channel, logarithmic stability for the drift strength, and Lipschitz stability for the complete diffusion strength including its sign. For multiple sources, we obtain Lipschitz stability up to permutation for drift locations and weighted temporal moments in dimensions two and three, and for diffusion locations and strengths in dimensions one through three. The proof employs explicit adjoint probes, including complex-isotropic polynomial annihilators in dimensions two and three and spectral synthesis in one dimension, combined with boundary-to-source estimates derived from Green’s identities. These results extend and improve upon existing deterministic stability theories by revealing the fundamentally different temporal information carried by stochastic integrals, a distinction that is further illustrated by numerical experiments.

54. 旋度-Yamabe 问题

The curl-Yamabe problem

作者: Guofang Wang, Mingwei Zhang

类别: math.DG, math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13998v1

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

摘要: We introduce and study a Yamabe-type variational problem associated with the curl operator on middle-degree forms on a closed oriented Riemannian manifold $$(M^n,g)$$ with $n\equiv3\pmod4$ . The corresponding curl–Yamabe constant is defined by minimizing a conformally invariant, gauge-invariant quotient built from $\|{\rm curl}\,α\|_{\frac{2n}{n+1}}$ and $\int_M\langle{\rm curl}\,α,α\rangle$ , and its Euler–Lagrange equation is the critical curl–Yamabe equation ${\rm curl}\,α=\frac{n+1}{2}|α|^{\frac{2}{n-1}}α$ . Using sharp curl–Sobolev inequalities on the sphere and an Aubin-type cut-off construction, we prove that $0<Y_{\rm curl}(M,[g])\le Y_{\rm curl}(\mathbb{S}^n)$ for every such manifold, and that the infimum is attained whenever $Y_{\rm curl}(M,[g])<Y_{\rm curl}(\mathbb{S}^n)$ . Moreover, for $$n>3$$ we establish the strict inequality for manifolds that are not locally conformally flat, via a test-form expansion detecting the Weyl tensor, and deduce existence of a minimizer solving the curl–Yamabe equation.

55. 刺猬唯一地最小化所有 $$p$$ -能量

The hedgehog uniquely minimises every $$p$$ -energy

作者: Siran Li

类别: math.AP, math.DG

发布日期: 2026-09-12

arXiv ID: 2609.13979v1

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

摘要: We prove that for every dimension $n \geq 2$ and parameter $p \in ]1,n[$ , the hedgehog $$u_0(x) = x/|x|$$ is the unique minimiser of the $$p$$ -energy $\mathcal{E}_p [u]:=\int_{\mathbf{B}^n}|\nabla u|^p\,{\rm d}x$ in the class of $W^{1,p}$ -mappings from the $$n$$ -dimensional unit ball $\mathbf{B}^n$ to the $$(n-1)$$ -dimensional unit sphere $\mathbf{S}^{n-1}=\partial\mathbf{B}^n$ whose trace equal to the identity on $\mathbf{S}^{n-1}$ . This answers in the affirmative an open question on p.29 in [Robert M. Hardt, Singularities of harmonic maps, Bull. Amer. Math. Soc. (N.S.) 34 (1997), 15–34].

56. 噪声作用下随机热方程的不爆破

Non-blowup of stochastic heat equations by noise

作者: Yunfeng Chen, Tusheng Zhang

类别: math.PR, math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13963v1

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

摘要: Consider the heat equation: \begin{equation}\label{eq:00} \begin{cases} \displaystyle du(t,x)=\frac12Δu(t,x),dt+b(u(t,x)),dt, &t>0,\ x\in D,\[1mm] u(t,x)=0,&t\ge0,\ x\in\partial D, u(0,x)=u_0(x),&x\in D. \end{cases} \end{equation} It is well known that superlinear growth of the coefficient $$b$$ will result in a blowup of the solution $$u$$ at a finite time. In this paper, we show that a random noise will prevent the explosion of the solution. More precisely, we show that the stochastic heat equation: \begin{equation}\label{eq:01} \begin{cases} \displaystyle du(t,x)=\frac12Δu(t,x),dt+b(u(t,x)),dt+u(t,x)^n,dB_t, &t>0,\ x\in D,\[1mm] u(t,x)=0,&t\ge0,\ x\in\partial D, u(0,x)=u_0(x),&x\in D. \end{cases} \end{equation} has a global solution if \begin{equation} z,b(z)\leq C_b(1+z^2)+η|z|^{2n}, \qquad z\in\R. \end{equation} for some $η\in [0, \frac{1}{2})$ , here $$n$$ is an arbitrary, but fixed integer, $B_t, t\geq 0$ is a Brownian motion. Truncation and comparison techniques play an important role.

57. 格图 $\mathbb{Z}^d$ 上含 Hardy 算子的 Poisson 问题及其在半线性方程中的应用

Poisson problem with Hardy operator on lattice graph $\mathbb{Z}^d$ and the application to semilinear equations

作者: Huyuan Chen, Bobo Hua, Wendi Xu

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13933v1

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

摘要: In this paper, we investigate the analytic properties of the Hardy-type operator $$-Δ+ μH_0$$ on the $$d$$ -dimensional integer lattice graph $\mathbb{Z}^d$ , where $μ\geq -1$ and $$H_0$$ is a critical Hardy potential at infinity–arising naturally from discrete Hardy inequalities. We derive sharp fixed-pole Green kernel estimates on $\mathbb{Z}^d$ and obtain a necessary and sufficient weighted summability condition characterizing the solvability of the nonnegative Poisson problem associated with this operator. Building upon this classification, we provide a complete characterization of the nonexistence of positive solutions to the semilinear elliptic inequality

-Δu + μH_0 u \geq W u^p \quad \text{in } \mathbb{Z}^d,

where $$p > 0$$ , $d \geq 3$ , $μ\geq -1$ , and $H_0, W \in C(\mathbb{Z}^d)$ are strictly positive potentials satisfying the asymptotic behaviors $H_0(x) \sim |x|^{-2}$ and $W(x) \sim |x|^θ$ as $|x| \to \infty$ , with $θ\in \mathbb{R}$ .

58. 球面上凸域的 Alexandrov–Fenchel 不等式

Alexandrov–Fenchel inequalities for convex domains in the sphere

作者: Tianci Luo, Yong Wei, Rong Zhou

类别: math.DG, math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13829v1

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

摘要: We prove the Alexandrov–Fenchel inequalities between any two spherical quermassintegrals for smooth weakly convex domains contained in an open hemisphere, with equality if and only if the domain is a geodesic ball. For strictly convex hypersurfaces, we introduce a globally constrained curvature flow which preserves one quermassintegral and decreases the next one. We establish uniform curvature estimates by combining a pinching estimate, a support function argument, and spherical polarity. As a consequence, the flow exists for all time and converges smoothly and exponentially to a geodesic sphere. The monotonicity of the quermassintegrals gives the full family of Alexandrov–Fenchel inequalities. A short-time mean curvature flow approximation and a localized rigidity argument extend the result, including the equality characterization, to weakly convex domains.

59. 具有临界负指数的双调和方程正破裂解的刚性

Rigidity of positive rupture solutions to a biharmonic equation with critical negative exponent

作者: Xia Huang, Yahui Jiang, Xianmei Zhou

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13820v1

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

摘要: We establish two rigidity theorems for positive rupture solutions of the conformally invariant equation $Δ^2 u=u^{-7}$ in $\mathbb R^3\setminus\{0\}$ , which extend continuously to the origin with $$u(0)=0$$ . First, we prove that if the associated conformal metric $g=u^{-4}|dx|^2$ has nonnegative scalar curvature, then every such solution is radially symmetric and has the sharp rupture profile $u(x)\sim(\frac{4}{3})^{\frac{1}{4}}|x|^{\frac{1}{2}}$ as $x\to 0$ ; in particular, the metric is complete at the origin. The principal novelty is a global rigidity theorem requiring neither curvature nor symmetry: the single global condition $$u(x)=o(|x|)$$ at infinity forces $u(x)\equiv(\frac{4}{3})^{\frac{1}{4}}|x|^{\frac{1}{2}}$ . This result provides a sharp answer to the uniqueness question posed by McKenna and Reichel in a class with no a priori symmetry assumption. Finally, we construct complementary examples demonstrating the essential roles of the curvature and growth hypotheses.

60. Bogovskii 常数的刚性定理与稳定性

Rigidity and stability for the Bogovskii constant

作者: Glen Wheeler

类别: math.AP, math.FA

发布日期: 2026-09-12

arXiv ID: 2609.13763v1

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

摘要: We prove that balls are the unique minimisers of the Bogovskii constant among bounded connected Lipschitz domains in $\mathbb{R}^n$ , $n\geq 2$ , resolving Open Problem 2 of Gazzola, Grunau and Sperone (arXiv:2608.25861).

61. 一个非局部灌注门控血管生成系统：整体弱解、快速信号极限与入侵前沿

A Nonlocal Perfusion-Gated Angiogenesis System: Global Weak Solutions, Fast-Signal Limits, and Invasion Fronts

作者: Jiguang Yu, Louis Shuo Wang

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13753v1

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

摘要: Continuum angiogenesis systems often use local endothelial or vessel density as an oxygen-delivery proxy without distinguishing structurally formed vessels from pressure-supported vascular function. We formulate a two-dimensional model in which lumenized density determines a normalized nonlocal conductivity, a globally solved pressure field, and a flow-functional density that gates oxygen delivery and vessel regression. For fixed positive regularization parameters and signal relaxation times, we prove Lipschitz stability of the pressure–perfusion map and global bounded weak solvability; strong vessel compactness follows from a nonlocal ordinary differential equation stability estimate rather than spatial smoothing. As the oxygen and vascular endothelial growth factor timescales vanish simultaneously, weak solutions converge along a subsequence to a parabolic–elliptic–ordinary differential equation system, with both fast fields converging strongly in $$L^2(0,T;H^1)$$ . A locally frozen one-dimensional reduction admits monotone fronts precisely at or above the threshold $B+2\sqrt{DR}$ for $B\geq-\sqrt{DR}$ and $$-DR/B$$ otherwise. Finite-volume experiments verify the analytical mechanisms and show that equal-mass vessel fields can generate distinct functional masses and oxygenation levels.

62. 含非线性梯度源的半线性椭圆方程的梯度估计与 Liouville 型定理

Gradient estimates and Liouville-type theorems for the semilinear elliptic equation involving the nonlinear gradient source

作者: Wenguo Liang, Zhengce Zhang

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13738v1

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

摘要: We study local and global properties of positive solutions to the equation $-Δu=u^p+M|\nabla u|^q$ in a domain $$Ω$$ of $\mathbb R^N$ , where $$p,q$$ are parameters and $$M>0$$ . By constructing a linear operator, we establish the differential inequality containing an auxiliary function. By selecting appropriate auxiliary functions over various regions and employing the maximum principle, we derive the local gradient estimates for all $(p,q)\in \mathbb R^2$ , and further establish Liouville-type theorems. As an application, we acquire universal estimates for local solutions of elliptic equations with general nonlinearities. Our results extend partial conclusions established in Bidaut-Véron, Garcia-Huidobro and Véron [Math. Ann. 378 (1-2) (2020) 13-56].

63. 混合局部与非局部抛物方程的梯度自改进

Gradient self-improvement for mixed local and nonlocal parabolic equations

作者: Kenta Nakamura

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13722v1

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

摘要: We introduce a new method for studying gradient higher integrability for mixed local and nonlocal parabolic equations. More precisely, for $$p>2d/(d+2)$$ and $s \in (0,1)$ , we prove that if the inhomogeneity $F \in L^{p(1+σ)}_{\mathrm{loc}}$ for some $$σ>0$$ , then every weak solution satisfies $\nabla u \in L^{p(1+\eps)}_{\mathrm{loc}}$ for some $\eps>0$ , together with quantitative local estimates. The proof relies on stopping time arguments and an intrinsic Calderón-Zygmund-type covering decomposition in which the inhomogeneity and nonlocal energy contributions are treated by classical and fractional maximal estimates.

64. 一般聚焦源间接阻尼波-MGT系统的势阱动力学中的稳定与不稳定

Stable and Unstable Potential-Well Dynamics for an Indirectly Damped Wave-MGT System with General Focusing Sources

作者: Tae Gab Ha

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13708v1

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

摘要: We study a conservative semilinear wave equation coupled through a zero-order interaction to a dissipative Moore–Gibson–Thompson equation on a bounded domain. The wave component carries no direct damping and is driven by a general focusing source $$f(u)$$ . The augmented variable $$w=v+τv_t$$ reveals an exact coupled energy and a coercive potential-well geometry. The source assumptions are formulated through $H_θ(s)=\frac1θsf(s)-F(s)$ , $F(s)=\int_0^s f(r)\,d r$ , $$θ>2$$ . Under $$L^2$$ -subcritical $$C^1$$ growth, smallness at the origin, nonnegativity and radial monotonicity of $$H_θ$$ , and a nontrivial focusing condition, we establish local well-posedness, the exact energy identity, and a continuation alternative for arbitrary finite-energy data. For nonzero coupling, the linearized semigroup is strongly stable, whereas a wave-branch expansion precludes uniform exponential stability and positive-time compactness. Below the coupled well depth, the stable set is positively invariant and generates global solutions, while data with negative Nehari functional blow up in finite time without a sign condition on the initial velocities. At the critical level $$E(0)=d$$ , nonzero coupling yields a complete trichotomy into stable entry, finite-time blow-up, or a stationary Nehari ground state. Under the same coupling condition, every stable trajectory converges weakly to zero without any compactness hypothesis. A renormalized high-frequency identity shows that vanishing of the accumulated nonlinear high–low flux is equivalent to relative compactness of the orbit and to strong convergence in the natural energy space. We give several sufficient criteria, including finite total variation of the nonlinear force in $$L^2(Ω)$$ .

65. 穿孔圆盘上 Liouville 方程解的参数化

A Parametrization of Solutions to the Liouville Equation in a Punctured Disc

作者: Qing Han, Jui-Yun Hung

类别: math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13702v1

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

摘要: In this paper, we study parametrization of solutions to the Liouville equation near their isolated singular points. We demonstrate that prescribing asymptotic behaviors forms a new type of problem leading to all such solutions. The harmonic functions determining the asymptotic behaviors play the role of ``free data" as the boundary values in the boundary value problems.

66. 中间 Ricci 曲率下的 Jacobian 估计与几何不等式

Jacobian estimates and geometric inequalities under intermediate Ricci curvature

作者: Kwok-Kun Kwong, Jihye Lee, Fabio Ricci

类别: math.DG, math.AP

发布日期: 2026-09-12

arXiv ID: 2609.13690v1

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

摘要: We prove a Jacobian estimate underlying both the Heintze–Karcher comparison and the Alexandrov–Bakelman–Pucci method for submanifolds of arbitrary codimension. The estimate retains the contributions of ambient curvature along geodesics emanating from the submanifold and requires no curvature sign assumption. As applications, we obtain a quantitative Fenchel–Willmore inequality with explicit curvature remainders under nonnegative $$n$$ -intermediate Ricci curvature, Michael–Simon Sobolev and isoperimetric inequalities under nonnegative $$(n-1)$$ -intermediate Ricci curvature, and corresponding extensions under quadratic curvature decay. In each result, the required intermediate Ricci curvature condition depends only on the dimension of the submanifold and is independent of its codimension.

67. 奇异合作系统的解的最优正则性

Optimal regularity of solutions to a singular cooperative system

作者: Lili Du, Xu Tang, Cong Wang

类别: math.AP

发布日期: 2026-09-11

arXiv ID: 2609.13604v1

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

摘要: In the pioneer work of vectorial free boundary problems \cite[Adv. Math. 280 (2015)]{ASUW15}, Andersson, Shahgholian, Uraltseva, and Weiss investigated the regularity theory of solutions and the free boundary of a singular cooperative system \begin{equation*} Δ\mathbf{u} = \frac{\mathbf{u}}{|\mathbf{u}|} χ_{{|\mathbf{u}|>0}} \quad\text{in }Ω\subset\mathbb{R}^n, \qquad \mathbf{u}:Ω\longrightarrow \mathbb{R}^m, \qquad n,m\ge 2. \end{equation*} At the level of solutions, the boundedness of the right-hand side directly yields $W^{2,p}_{\mathrm{loc}}$ regularity for every $1<p<\infty$ by standard elliptic estimates, while the corresponding $W^{2,\infty}_{\mathrm{loc}}$ regularity question was explicitly left as an open problem on page~753 of that paper. In this paper, we give a positive answer by establishing optimal $W^{2,\infty}_{\mathrm{loc}}$ regularity through a scale-invariant interior Hessian estimate for every weak solution. Unlike the scalar case, the classical Alt–Caffarelli–Friedman monotonicity formula cannot be applied directly to the present system. Our proof instead combines a localized Newtonian potential decomposition with an explicit projection onto quadratic harmonic polynomials on the unit sphere, which isolates the symmetric trace-free quadratic part of the Hessian. A key new observation is that the resulting matrix coefficients satisfy an exact linear ordinary differential equations with constant coefficients on the logarithmic scale. The main novelty of this work lies in a so-called two-regime argument based on the competition between the affine and quadratic parts of the rescaled solution, together with a continuity argument that connects the two regimes and yields a uniform bound for the coefficient. The resulting a priori estimate provides an analytic basis for further study of the fine structure of the free boundary.

68. 合成洛伦兹时空中的Monge传输问题

Monge’s transport problem in synthetic Lorentzian spacetimes

作者: Afiny Akdemir

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

发布日期: 2026-09-11

arXiv ID: 2609.17592v1

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

摘要: We prove the existence of solutions to the Lorentzian Monge problem on synthetic Lorentzian spacetimes that satisfy the forward timelike measure contraction property $$TMCP^+(K,N)$$ . The cost for this problem is the time separation function $\ell(x,y)$ (the Lorentz distance), which represents the maximum amount a particle can age when traveling from $$x$$ to $$y$$ . The solution is based on the needle decomposition technique, and amounts to a reduction of the full problem to its one-dimensional counterparts. As an application of needle decomposition, we show that the Lorentzian timelike measure contraction property $$TMCP^+(K,N)$$ is equivalent to the $TMCP^+_{rLip}(K,N)$ condition, where the curvature dimension conditions are defined on gradient flow curves of reverse 1-Lipschitz functions.

69. 各向异性 Calderón 问题中欧几里得度量的局部刚性

Local rigidity of the Euclidean metric for the anisotropic Calderón problem

作者: Plamen Stefanov

类别: math.AP

发布日期: 2026-09-11

arXiv ID: 2609.13400v1

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

摘要: We prove local rigidity of the Euclidean metric for the anisotropic Calderón problem on smooth compact domains $M\subset \mathbb{R}^n$ , $n\ge3$ . A smooth Riemannian metric with the same induced boundary metric and Dirichlet-to-Neumann map as the Euclidean metric is isometric to it by a diffeomorphism fixing the boundary pointwise, provided the two metrics are sufficiently close in $$H^s(M)$$ , with $$s>n/2+1$$ an integer.

70. 史瓦西时空附近的爱因斯坦-弗拉索夫系统中的非线性相位混合

Nonlinear Phase Mixing in Einstein-Vlasov System near Schwarzschild Spacetime

作者: Saehoon Eo

类别: math.AP, gr-qc

发布日期: 2026-09-11

arXiv ID: 2609.13394v1

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

摘要: We study the Einstein-Vlasov system near the Schwarzschild spacetime in spherical symmetry. In particular, we consider the case when the matter is supported on bounded geodesics. In this regime, the decay of the matter is driven by phase mixing rather than dispersion. For the linearized problem, we obtain quantitative phase mixing; the key new ingredient proved here is a monotonicity property of the period function. For the nonlinear coupled problem, we prove phase mixing up to the timescale expected from the nonlinear echo mechanism; this is achieved by analysis on dynamical action angle variables and the integral estimates with the vector field method. This construction is made possible by our gauge choice, an isotropic coordinate system together with a foliation of prescribed mean curvature, in which the Einstein equations gain derivatives over the matter due to its elliptic nature.

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

1. 一维可压缩 Navier–Stokes 方程在一般压强定律下强解的局部与整体存在性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. Navier-Stokes-Maxwell 系统整体解的最优时间衰减#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 超粘性可压缩三维 Oldroyd 型模型的尖锐非唯一性：超越 Lions 指数#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 具有流入-流出边界条件的可压缩完整Navier-Stokes方程的爆破判据#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 关于洛伦兹空间 $L^{3,\infty}(\mathbb{R}^3)$ 中 Navier-Stokes 方程解的唯一性与正则性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 三维磁流体力学边界层方程的线性不适定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 二维惯性磁流体力学中的线性与非线性不稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. 临界空间中三维各向异性 Navier-Stokes 方程解的稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

9. 复合状态方程的可压缩欧拉方程的波速保证界#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

10. 三维可压缩非耗散磁流体力学的整体稳定性：隐藏耗散与有理背景场#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

11. 可压缩欧拉方程的结构稳定奇点与 Lipschitz 稳定的最优传输度量#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 三角翼空气动力特性的分析#

摘要翻译#

2. 加权亚强制算子的函数空间#

摘要翻译#

3. 关于正压可压缩原始方程的固有静水能量结构#

摘要翻译#

4. 非散度型动力学方程的 Harnack 估计#

摘要翻译#

5. 二维可压缩微极流体方程中大体粘滞与微小角粘滞的联合极限#

摘要翻译#

6. 半线性波动方程弱解的尖锐非唯一性#

摘要翻译#

7. 平衡态下耦合辐射-欧拉方程的希尔伯特展开#

摘要翻译#

8. 环面与波导流形上单粒子薛定谔方程及无穷多粒子薛定谔方程的逐点收敛#

摘要翻译#

9. 正交各向异性、宽退化双非线性扩散方程的有界性与收缩估计#

摘要翻译#

10. Vlasov–Yukawa 方程非线性 Landau 阻尼中的大盒子过渡#

摘要翻译#

11. 带排斥逆平方势的Schrödinger方程的双线性估计及其应用#

摘要翻译#

12. 无偶性假设的渐近双曲面的留数估计#

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

1. 一维可压缩 Navier–Stokes 方程在一般压强定律下强解的局部与整体存在性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. Navier-Stokes-Maxwell 系统整体解的最优时间衰减

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 超粘性可压缩三维 Oldroyd 型模型的尖锐非唯一性：超越 Lions 指数

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 具有流入-流出边界条件的可压缩完整Navier-Stokes方程的爆破判据

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 关于洛伦兹空间 $L^{3,\infty}(\mathbb{R}^3)$ 中 Navier-Stokes 方程解的唯一性与正则性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 三维磁流体力学边界层方程的线性不适定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 二维惯性磁流体力学中的线性与非线性不稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 临界空间中三维各向异性 Navier-Stokes 方程解的稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

9. 复合状态方程的可压缩欧拉方程的波速保证界

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

10. 三维可压缩非耗散磁流体力学的整体稳定性：隐藏耗散与有理背景场

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

11. 可压缩欧拉方程的结构稳定奇点与 Lipschitz 稳定的最优传输度量

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 三角翼空气动力特性的分析

摘要翻译

2. 加权亚强制算子的函数空间

摘要翻译

3. 关于正压可压缩原始方程的固有静水能量结构

摘要翻译

4. 非散度型动力学方程的 Harnack 估计

摘要翻译

5. 二维可压缩微极流体方程中大体粘滞与微小角粘滞的联合极限

摘要翻译

6. 半线性波动方程弱解的尖锐非唯一性

摘要翻译

7. 平衡态下耦合辐射-欧拉方程的希尔伯特展开

摘要翻译

8. 环面与波导流形上单粒子薛定谔方程及无穷多粒子薛定谔方程的逐点收敛

摘要翻译

9. 正交各向异性、宽退化双非线性扩散方程的有界性与收缩估计

摘要翻译

10. Vlasov–Yukawa 方程非线性 Landau 阻尼中的大盒子过渡

摘要翻译

11. 带排斥逆平方势的Schrödinger方程的双线性估计及其应用

摘要翻译

12. 无偶性假设的渐近双曲面的留数估计