生成时间: 2026年07月21日 05:43:18

运行状态: 完成

计划处理论文数量: 51

已完成论文数量: 51

未完成论文数量: 0

重点关注论文数量: 13

了解领域论文数量: 19

不相关论文数量: 19

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

1. 具有速度依赖阻尼的 Euler-Poisson 方程

Euler-Poisson equations with velocity-dependent damping

作者: Olga S. Rozanova

类别: math.AP, math-ph

发布日期: 2026-07-20

arXiv ID: 2607.18035v1

摘要翻译: 我们考虑一维空间以及具有径向对称性的高维情形下的排斥性 Euler-Poisson 方程，并假设阻尼系数为速度的幂律函数。我们证明：在此假设下，一维情形中对应全局时间光滑解的初值集合扩大；而在其他维度（除四维外），该阻尼并不能改善解的光滑性。具体而言，任意小的稳态扰动均会发生爆破，尽管其振荡幅度衰减。

Comment: 14 pages, 4 figures

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

1. 研究对象和背景

本文研究如下具有速度依赖阻尼项的排斥型欧拉-泊松方程组（即描述冷电子等离子体的无压模型）：

其中 $N > 0$ 为电子密度，$V \in \mathbb{R}^d$ 为速度场，$\Phi$ 为排斥势（$k > 0$），$N_0 > 0$ 为背景密度（文中取 $N_0 = 1$），阻尼系数为幂律形式 $\nu(|V|) = -\nu_0 |V|^k$，$\nu_0 > 0$, $k > 0$。通过引入电场 $E = -\nabla \Phi$ 并利用约束 $N = 1 - \operatorname{div} E$（物理要求 $\operatorname{div} E < 1$），系统可等价重写为闭合一阶双曲型系统：

该模型在等离子体物理（冷电子近似）、半导体输运及天体物理中具有基本意义；其数学核心难点在于：系统是非严格双曲的，特征线存在但退化，解的奇性形成机制由梯度爆炸（即 $|\nabla V|$ 或 $|\nabla E|$ 在有限时间趋于无穷）主导，而非密度零点或真空形成。关键未知量是初始数据 $(V_0, E_0)$ 的结构如何决定全局光滑解的存在性——即刻画所谓“光滑性区域”（smoothness domain）或“次临界区域”（subcritical region）：该区域是 $C^1$ 初值空间中所有导致全局 $C^1$ 解的初值构成的集合。

已有理论表明，该区域的存在性对空间维数 $d$ 极其敏感：当 $\nu \equiv 0$ 时，仅在 $d = 1$ 和 $d = 4$ 时存在非空光滑性区域（$d=1$ 时由阈值曲线 $q_0^2 > 1 - 2 s_0$ 精确刻画，$d=4$ 时结构复杂）；其余维数下，任意小的非零扰动均导致有限时间爆破。引入常数阻尼 $\nu = \nu_0 > 0$ 后，所有维数下光滑性区域均扩张，且解呈指数衰减 $V,E \sim e^{-\nu_0 t}$。若阻尼依赖于密度 $\nu = \nu(N)$，则在一维情形下可通过积分判据保证全局光滑性。然而，速度依赖阻尼 $\nu = \nu_0 |V|^k$ 的影响此前未被系统研究：它既非常数亦非密度函数，其非线性耗散强度随流速动态变化，在湍流、非牛顿流体建模中具物理合理性，但其对奇性抑制能力的维数依赖性尚属未知。本文旨在填补这一空白：严格证明该类阻尼在 $d = 1$（及 $d = 4$）下可扩张光滑性区域，而在其他维数下完全失效——即无法阻止任意小扰动的爆破，尽管其振幅衰减。此现象挑战了“耗散必改善正则性”的直觉，根源在于幂律衰减与导数层面的 Floquet 型指数不稳定性之间的根本冲突。

2. 主要定理或主要结果

本文的核心结论按维数分述如下：

1. 一维情形（$d = 1$）：设初值 $(V_0, E_0) \in C^2(\mathbb{R})$ 满足 $s_0 = \partial_x E_0 < 1$（物理约束），且 $\nu(|V|) = -\nu_0 |V|^k$，$k = 2m$, $m \in \mathbb{N}$。则：

   • 长时间渐近行为（Theorem 1 & 2）：若解全局光滑，则沿任意特征线 $x = x(t)$，速度满足
     $$ V(x(t)) \sim t^{-1/k} \sum_{n=0}^\infty \Psi_n(t) t^{-n}, \quad t \to \infty, $$
     其中 $\Psi_n$ 为周期函数；其导数 $q = \partial_x V$, $s = \partial_x E$ 满足
     $$ q(t), s(t) \sim t^{-\nu_0 A^k a_k / 2} \left( A \sin(t + B) + O(t^{-1}) \right), \quad t \to \infty, $$
     其中 $a_k = (k+1)!!/k!!$, $|A| \le (\nu_0 k / 2)^{-1/k}$。特别地，$V$ 衰减率随 $k$ 增大而变慢，而 $q,s$ 衰减率随 $k$ 增大而加快（Remark 1）。
   • 阈值曲线精确构造（Theorem 3）：对零初速扰动 $V_0 \equiv 0$，在 $(s_0, E_0)$ 平面上，光滑性区域与爆破区域的分界曲线（threshold curve）在 $(1/2, 0)$ 附近具有展开式
     $$ \delta = s_0 - \frac{1}{2} = \nu_0 C_k \epsilon^k + o(\epsilon^k), \quad \epsilon = E_0, \quad C_k = \frac{k+1}{4} \int_0^\pi \sin^{k+2} \tau \, d\tau > 0. $$
     该曲线在 $\epsilon \to 0$ 时严格依赖于 $k$ 和 $\nu_0$，且 $C_k$ 随 $k$ 单调递增。

2. 径向对称高维情形（$d \ge 2$）：设解具径向形式 $V = F(t,r) r$, $E = G(t,r) r$，初值满足 (21)，且 $\nu(|V|) = -\nu_0 |V|^k$。则：

   • 爆破定理（Theorem 4）：对任意 $d \neq 1$, $d \neq 4$，存在 $\delta_0 > 0$，使得对所有满足 $\|V_0\|_{C^1} + \|E_0\|_{C^1} < \delta_0$ 的非平凡初值，其对应解必在有限时间 $T^* < \infty$ 内发生梯度爆破，即 $\sup_{t \in [0,T^*)} (|\nabla V| + |\nabla E|) = \infty$。爆破时间 $T^*$ 随维数 $d$ 和指数 $k$ 增大而减小。
   • 四维特例：文中明确指出 $d = 4$ 时行为类比 $d = 1$，即阻尼可扩张光滑性区域，但未给出具体定理陈述。

所有结论均建立在 $C^1$ 初值框架下，且爆破定义为一阶导数（$q,s$ 或其高维类比）在有限时间趋于无穷。

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

本文证明的核心障碍在于：速度依赖阻尼 $\nu(|V|) = -\nu_0 |V|^k$ 破坏了系统的线性耗散结构，导致标准能量估计失效；同时，高维径向情形下无法直接分离变量，必须处理耦合的非线性 ODE 系统。作者采用并发展了一套基于Radon 引理的特征线约化与线性化框架，其创新点与关键技术如下：

1. 特征线上的精确 ODE 约化：对一维系统 (5)，沿特征 $\dot{x} = V$，将 PDE 系统降为关于 $(V,E,q,s)$ 的四维非线性 ODE 系统 (6)-(7)。关键创新在于显式计算出阻尼诱导的非线性反馈项 $a(|V|) = (|V| \nu(|V|))'$，当 $\nu = -\nu_0 |V|^k$ 时，$a(|V|) = -\nu_0 (k+1) |V|^k$，此表达式直接进入 $q,s$ 的演化方程，成为后续分析的枢纽。

2. Bihari-LaSalle 型非线性 Gronwall 估计：为控制 $V,E$ 的振幅，作者对 $z(t) = V^2(x(t))$ 建立微分不等式 $\dot{z} \le -2\nu_0 z^{k/2 + 1}$，并通过变量代换 $v(t) = \alpha - 2\nu_0 \int_0^t w(z(\tau)) d\tau$（其中 $w(z) = z^{k/2 + 1}$）导出关键先验估计 (8)-(9)：

   $$ V^2(x(t)) \le \left( (V_0^2 + E_0^2)^{-k/2} + \nu_0 k t \right)^{-2/k}. $$

   此估计首次定量刻画了幂律阻尼下解振幅的代数衰减率 $t^{-2/k}$，是后续渐近分析的基础。

3. Liénard 方程的奇异摄动渐近展开：将 (6) 中 $V$ 的演化视为 Lénard 型振荡器 $\ddot{V} + \nu_0(k+1)|V|^k \dot{V} + V = 0$，通过尺度变换 $V = t^{-1/k} v(t)$ 及形式级数 $v(t) = \sum_{n=0}^\infty \Psi_n(t) t^{-n}$（$\Psi_n$ 周期），将原方程转化为递推的 ODE 序列 $\ddot{v}_n + v_n = P_2(v_j, \dot{v}_j)$（$j < n$）。关键技巧在于识别并抵消由 $\dot{v}_{n-1}$ 生成的 secular terms（如 $t \cos t$），其被 $t^{-n}$ 因子精确补偿，从而保证一致有界性，最终确立 $V$ 的 $t^{-1/k}$ 主导衰减律。

4. Radon 引理驱动的导数线性化与 Floquet 分析：对 $q,s$ 的演化，作者应用 Radon 引理将其表示为线性系统 (14) 的商：$q = p_1/Q$, $s = p_2/Q$，其中 $(Q,p_1,p_2)$ 满足三阶线性 ODE。核心创新在于对 $p_2$ 方程 $\ddot{p}_2 + a(V) \dot{p}_2 + p_2 = 0$ 进行进一步变换：令 $p_2 = u v$, $v = \exp\left(-\frac{1}{2} \int_0^t a(V) d\tau\right)$，导出 $u$ 满足的方程 $\ddot{u} + (1 - \phi(t)) \dot{u} = 0$，其中 $\phi(t) = \frac{1}{4} a^2(V) + \frac{1}{2} a'(V)$。利用 (8) 得 $a(V) \sim \nu_0 A^k (k+1) t^{-1} \sin^k(t+B)$，进而证明 $\int_\delta^\infty \phi(t) dt < \infty$，从而引用经典 ODE 理论（[2, Ch.VI, Th.1]）得出 $u$ 有界且渐近于 $A \sin(t+B)$。结合 $v \sim t^{-\nu_0 A^k a_k / 2}$，最终获得导数的精细渐近式 (12)。

5. 高维爆破的线性化判据：在径向情形下，作者将 $D = \operatorname{div} V$, $\lambda = \operatorname{div} E$ 的演化（23）通过变量替换 $u = D - dF$, $v = \lambda - dG$ 化为 (24)，再应用 Radon 引理得到线性系统 (25)。爆破判据归结为标量函数 $Q(t)$ 是否零点：$Q(t) = 0$ 当且仅当 $\nabla V$ 或 $\nabla E$ 爆破。数值实验显示，对 $d \neq 1,4$，$Q(t)$ 的振幅呈指数增长并过零，而其增长率由 Floquet 乘子主导，幂律阻尼无法压制此指数机制。

4. 与之前工作的比较

本文与既有工作的本质区别体现在三个维度：阻尼结构、维数敏感性、以及衰减-不稳定性平衡机制。以下分层比较：

1. 与零阻尼工作（[6],[9],[11],[13],[14]）比较：

   • 当 $\nu \equiv 0$ 时，$d = 1$ 存在精确阈值 $q_0^2 > 1 - 2 s_0$；$d = 4$ 存在复杂阈值；其余维数无光滑性区域。本文证明：在 $d = 1$，幂律阻尼将阈值曲线从直线 $s_0 = 1/2$ 变形为 $\delta \sim \nu_0 C_k \epsilon^k$，即扩张了光滑性区域，且扩张量由 $k$ 和 $\nu_0$ 精确量化；在 $d \neq 1,4$，则完全不扩张，延续了零阻尼的奇性必然性。这是对 [11],[14] 结论的直接推广，但揭示了阻尼类型的决定性作用。

2. 与常数阻尼工作（[6],[9],[12]）比较：

   • 常数阻尼 $\nu = \nu_0 > 0$ 下，所有维数均存在光滑性区域，且解呈指数衰减 $V,E \sim e^{-\nu_0 t}$，导数亦指数衰减。本文发现：幂律阻尼仅在 $d = 1,4$ 有效，且衰减为代数阶 $t^{-1/k}$（$V$）和 $t^{-\nu_0 A^k a_k / 2}$（导数），其代价是丧失普适性——在 $d = 2,3,5,6,...$ 完全失效。这表明：耗散的‘强度’不等于其‘有效性’；常数阻尼提供均匀压制，而幂律阻尼在振幅小时耗散弱，无法遏制导数层面的 Floquet 不稳定性。

3. 与密度依赖阻尼工作（[10]）比较：

   • 文献 [10] 证明，若 $\nu = \nu(N)$ 且 $\int_0^1 \nu(\rho) d\rho = \infty$，则一维下任意初值均得全局光滑解。本文的 $\nu(|V|)$ 不依赖密度，故不适用此判据；其结论更精细：并非所有初值都光滑，而是存在一个非平凡的、可计算的光滑性区域，其边界由 $k$ 显式参数化。这体现了速度依赖阻尼的中间性质——强于零阻尼但弱于满足发散积分条件的密度阻尼。

4. 物理机制层面的突破：

   • 本文首次明确指出并数学证实：幂律衰减与 Floquet 指数增长的尺度不匹配是高维失效的根本原因。如引言所述，“振幅按幂律衰减，而导数不稳定性按指数增长，幂因子无法压制指数增长，只能延缓”。此论断在 Theorem 4 的证明中通过 $Q(t)$ 的数值爆破行为得到验证，并在渐近分析中体现为 $\phi(t)$ 的可积性（保证 $u$ 有界）与 $v(t)$ 的代数衰减（决定 $q,s$ 衰减率）的分离。这一机制在常数阻尼下不存在，因后者使整个系统谱移至左半平面。

综上，本文的推进在于：首次建立了速度依赖阻尼对欧拉-泊松方程奇性形成的维数选择性效应，给出了 $d = 1$ 下阈值曲线的首个显式幂级数展开，并揭示了代数耗散与指数不稳定性不可调和的本质冲突。其局限在于：$d = 4$ 的严格证明未给出；多维非径向情形未涉及；且所有结果依赖于 Radon 引理框架，对更一般初值（如非径向、非小扰动）的适用性未讨论。

2. Robin 边界条件下的 Couette 流（I）：与粘性无关的摩擦

Couette Flow with Robin Boundary Condition (I): the viscosity-independent friction

作者: Siming He, Binqian Niu, Weiren Zhao

类别: math.AP, math-ph

发布日期: 2026-07-20

arXiv ID: 2607.17468v1

摘要翻译: 本文是该系列论文的第一篇。在本系列文章中，我们将考察固–液边界处摩擦因子 $α$ 对 Couette 流稳定性的影响。具体而言，我们考虑在有界周期通道 $\mathbb{T} \times [-1,1]$ 中、Robin 型边界条件下 Couette 流的稳定性（即 $u^2|_{y=\pm 1} = 0$，$[α\partial_n u^1 + u^1]|_{y=\pm1} = f $），其中 $α$ 为摩擦因子，$n$ 为单位外法向量。本文证明：对给定摩擦因子 $α$，只要流体粘性系数 $ν\ll α$ 充分小，则当初始扰动满足 $\|ω_{\rm in}\| \leq εν^{1/3}$ 时，系统是渐近稳定的。此外，无粘阻尼（inviscid damping）与增强耗散（enhanced dissipation）均成立。

Comment: 无

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

1. 研究对象和背景

本文研究二维不可压缩Navier-Stokes方程在有限周期通道 $\mathbb{T} \times [-1,1]$ 上围绕Couette基流 $(y,0)$ 的非线性稳定性问题，其核心在于刻画Robin型边界条件对临界阈值的定量影响。设扰动速度 $v = (v_1,v_2)$ 满足

其中 $\nu > 0$ 为运动粘性系数（即雷诺数倒数 $\mathrm{Re}^{-1}$），$\kappa \in (0,1)$ 为无量纲滑移参数，$\partial_n = n \cdot \nabla$ 为外法向导数，且定义摩擦因子 $\alpha := \frac{1-\kappa}{\kappa} > 0$。该Robin条件等价于

它介于经典无滑移（$\kappa=1, \alpha=0$）与Navier自由滑移（$\kappa=0, \alpha=\infty$）之间，物理上可由Boltzmann方程在壁面原子结构处的边界层推导而来。研究对象是扰动涡量 $\omega = \partial_x v_2 - \partial_y v_1$，其演化满足vorticity形式

并辅以Robin型涡量边界约束：$-\alpha \omega(\pm 1) = \mp v_1(\pm 1)$，即 $T_\pm[\omega] := \alpha \omega(\pm 1) \pm \partial_y \Delta^{-1} \omega(\pm 1) = 0$，其中 $\Delta^{-1}$ 指Dirichlet Laplacian的逆。该问题的重要性在于：一方面，Couette流是唯一已知在无限域中具有严格线性稳定性的剪切流，其非线性稳定性阈值问题自Kelvin、Orr、Rayleigh以来即为核心课题；另一方面，在有界域中，边界效应（尤其是壁面涡量堆积）严重干扰经典机制——无滑移下$\omega$在$y=\pm 1$附近产生$O(\nu^{-1/2})$尖峰，迫使文献[14]引入加权范数；而自由滑移则因$\omega|_{y=\pm 1}=0$规避此障碍，恢复$O(\nu^{1/3})$阈值。然而，真实物理系统常处于二者之间：摩擦因子$\alpha$有限且固定，此时边界既非完全粘附亦非完全滑移，传统工具失效。本文试图填补的关键空白是：在$\alpha \in (0,\infty)$固定、$\nu \ll \alpha$小参数 regime 下，建立统一的、不退化至无滑移或自由滑移极限的稳定性理论，并精确刻画临界指数$\beta$如何依赖于$\alpha$——特别地，证明$\beta = 1/3$在Robin条件下仍成立，从而揭示“粘性无关的摩擦”这一反直觉现象：当$\nu$足够小时，稳定性阈值不再受$\alpha$尺度影响，仅由$\nu^{1/3}$主导。

2. 主要定理或主要结果

在固定摩擦因子 $\alpha > 0$ 且 $0 < \nu \leq \nu_0(\alpha)$ 充分小的假设下，对初始扰动 $v_{\rm in} \in H^6(\mathbb{T} \times [-1,1])$ 满足相容性约束 $\left(\alpha \partial_y v_{1,{\rm in}} \pm v_{1,{\rm in}},\, v_{2,{\rm in}}\right)|_{y=\pm 1} \equiv (0,0)$ 及小性条件 $\|v_{\rm in}\|_{H^6} \leq \epsilon \nu^{1/3}$，本文证明了全局非线性稳定性及双重耗散机制：

（1）主稳定性估计（Theorem 1.1）：存在与 $\alpha$ 相关的常数 $\epsilon_0 \in (0,1/4)$ 及通用常数 $C \geq 1$，使得对所有 $t \geq 0$， $$ |e^{\epsilon_0 \nu^{1/3} t} (1-|y|) \operatorname{curl} v_{\neq}|{L^\infty_t L^2{x,y}} + |e^{\epsilon_0 \nu^{1/3} t} v_{\neq}|{L^2_t L^2{x,y}} \

• \nu^{1/4} |e^{\epsilon_0 \nu^{1/3} t} \operatorname{curl} v_{\neq}|{L^2_t L^2{x,y}} + |\operatorname{curl} v_0|{L^\infty_t L^2{x,y}} \leq C \epsilon \nu^{1/3}. $$ 其中 $v_{\neq} = v - v_0$ 为非零模部分，$v_0$ 为 $x$-平均分量，$(1-|y|)$ 权重精确控制壁面附近涡量衰减。

（2）双重耗散机制：对涡量扰动 $\omega = \operatorname{curl} v$，存在 $\delta > 0$ 使得

• 无粘阻尼（I.D.）：$\|\nabla \Delta^{-1} \omega_{\neq}(t)\|_{L^2} \leq C (1+t)^{-1} \|\omega_{\rm in}\|_{H^1}$，即速度扰动 $v_{\neq}$ 的 $L^2$ 范数以 $O(t^{-1})$ 代数衰减，速率与 $\nu$ 无关；
• 增强耗散（E.D.）：$\|\omega_{\neq}(t)\|_{L^2} \leq C \|\omega_{\rm in}\|_{L^2} e^{-\delta \nu^{1/3} t}$，即涡量扰动在 $O(\nu^{-1/3})$ 时间尺度上指数衰减，远快于热传导尺度 $O(\nu^{-1})$。

（3）边界涡量控制：由Proposition 2.2及Corollary 2.5，边界校正项 $\omega^{(b)}$ 满足

表明边界效应被 $\nu^{1/3}$ 抑制，且权重 $(1-|y|)$ 精确捕获壁面层衰减率。

（4）共振控制：对长时共振分量 $\omega^{(re)}_{\neq}$（满足方程(2.13)），通过共振追踪泛函框架，证明其增长被 $e^{-\epsilon_0 \nu^{1/3} t}$ 指数压制，从而避免echo chain导致的norm inflation突破 $O(\nu^{1/3})$ 邻域。

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

本文的核心技术障碍有三：一是Robin边界下涡量在壁面不消失，导致经典Fourier分析失效且需处理强边界层；二是时间依赖背景剪切流 $U(t,y)$ 使线性算子非自治，传统冻结时间法引入显著误差；三是非线性共振（echo chain）在 $O(\nu^{-1/3})$ 时间尺度上可能引发norm inflation，威胁 $O(\nu^{1/3})$ 阈值。作者的真正创新点在于构建一套层级化、物理适配的分析框架，而非简单移植既有工具。

第一，层级振幅分解（Hierarchical Amplitude-based Decomposition）：将解 $v$ 分解为不同尺度项：零模 $\tilde{U}(t,y)$（由热方程驱动）、一级近似 $\omega^{(a)}$（内部主动力学）、一级边界校正 $\omega^{(b)}$（修复Robin条件）、一级余项 $\omega^{(*)}$；再对 $\omega^{(*)}$ 进行二级分解，分离出共振分量 $\omega^{(re)}$ 与高阶误差 $\omega^{(e)}$。此分解的关键在于将 $O(\nu^{1/3})$、$O(\nu^{2/3-})$ 项严格分离，并为每类项配备专属估计工具——例如，$\omega^{(a)}$ 用显式公式估计（Proposition 2.1），$\omega^{(b)}$ 用边界校正算子 $C$ 控制（Proposition 2.3），$\omega^{(re)}$ 用共振追踪泛函，$\omega^{(e)}$ 用奇异积分算子 $J_k$。该框架成功解耦了无粘阻尼、增强耗散、Orr机制与边界效应，确保各部分相互干扰可控。

第二，鲁棒边界校正算子 $C$：针对时间依赖剪切流，作者设计了冻结时间方案：将时间轴划分为区间 $I^{[j]} = [t_j, t_{j+1})$，其中 $t_j = j(\nu^{-1/3}+1)$，并在每个区间上冻结剪切 $U^{[j]}(y) = U(t_j,y)$。边界校正算子定义为 $C[g_{\neq}] = \sum_{j'=0}^j C^{[j']}[g_{\neq}]$，其中 $C^{[j']}$ 满足

$\chi_j$ 为光滑截断函数。关键创新在于误差泛函 $E[g_{\neq}]$ 的构造（公式(2.21)），它量化了冻结剪切与真实剪切的偏差，并通过Lemma 2.4证明其 $L^2$ 范数被 $\nu^{7/3}|k|^{-1/3}$ 控制，远小于主项，从而保证整个校正过程的鲁棒性。

第三，物理侧奇异积分算子（Physical-side SIOs）：为捕捉无粘阻尼与增强耗散，作者在物理域而非频域定义SIOs：

• 无粘阻尼算子 $J_k[f_k](y) = k \, \mathrm{P.V.} \int_{-1}^1 \frac{G_k(y,y')}{2i(U(t,y)-U(t,y'))} f_k(y') \, dy'$，
• 增强耗散算子 $J^{(e)}_k[f_k](y) = k \, \mathrm{P.V.} \int_{-1}^1 \frac{\nu^{1/3}|k|^{2/3} G^{(e)}_k(y,y')}{2i(U(t,y)-U(t,y'))} f_k(y') \, dy'$， 其中 $G_k$ 和 $G^{(e)}_k$ 分别为算子 $\partial_y^2 - |k|^2$ 与 $\nu^{2/3}|k|^{4/3} \partial_y^2 - |k|^2$ 的Dirichlet Green函数。这些算子直接作用于物理空间，避免了Fourier变换在Robin边界下的困难，并通过能量泛函 $E_k(t) = \|f_k\|_{L^2}^2 + c_\alpha \langle f_k, J_k[f_k] \rangle + c_\beta \langle f_k, J^{(e)}_k[f_k] \rangle$ 实现耗散机制的定量捕获。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度：临界指数、边界处理、时间尺度与适用范围。

临界指数 $\beta$ 的推进：在无界域 $\mathbb{T} \times \mathbb{R}$，[12,33,44] 证明 $\beta = 1/3$（$H^3$ 初值）；在有界域，无滑移下[14]得 $\beta = 1/2$（加权 $H^3$），而自由滑移下[11,45]恢复 $\beta = 1/3$。本文首次证明：对任意固定 $\alpha \in (0,\infty)$，只要 $\nu \ll \alpha$，仍有 $\beta = 1/3$，即 $\|v_{\rm in}\|_{H^6} \leq \epsilon \nu^{1/3}$ 即可保证全局稳定性。这打破了“有界域必然退化至 $\beta=1/2$”的固有认知，揭示了Robin摩擦的特殊性——当粘性足够小时，摩擦效应被耗散主导，阈值回归无界域标准。

边界处理的范式革新：以往工作或回避边界（无界域）、或极端化边界（无滑移/自由滑移）。无滑移需引入 $y$-加权范数抑制壁面涡量尖峰；自由滑移则因 $\omega|_{y=\pm 1}=0$ 自然规避。本文首创“冻结时间+边界校正算子 $C$”框架，直接处理 $T_\pm[\omega]=0$ 这一非齐次约束。相比[14]的加权估计，本文的 $(1-|y|)$ 权重更精确匹配Robin层厚度 $O(\nu^{1/3})$；相比[11]的自由滑移，本文的 $C$ 算子能处理 $\alpha$ 有限情形，且其误差 $E[g]$ 的 $\nu^{7/3}$ 阶控制造成了根本性优势。

时间尺度的统一：无界域与自由滑移的增强耗散时间为 $O(\nu^{-1/3})$，无滑移下因边界层干扰，文献[14]仅得到积分意义下的耗散。本文首次在Robin条件下严格建立 $O(\nu^{-1/3})$ 指数衰减 $\|\omega_{\neq}(t)\|_{L^2} \lesssim e^{-\delta \nu^{1/3} t}$，且该尺度与 $\alpha$ 无关，证实了“粘性无关的摩擦”——摩擦因子 $\alpha$ 仅影响常数 $\nu_0(\alpha)$ 和 $\epsilon_0(\alpha)$，但不改变耗散速率的 $\nu$-幂律。

适用范围与代价：本文结论要求 $\nu \leq \nu_0(\alpha)$，即粘性必须远小于摩擦，这比自由滑移（$\alpha=\infty$）条件更苛刻，但比无滑移（$\alpha=0$）更具物理普适性。技术代价是正则性要求提升至 $H^6$（高于无界域的 $H^3$ 或自由滑移的 $H^4$），源于边界校正与共振控制的双重复杂性。此外，从提供内容中无法确认本文结果是否可推广至三维或其它剪切流，亦未显式展开与[13]中计算机辅助发现的 $O(\nu^{1/2})$ 边界不稳定性阈值的定量比较。

3. 具有常涡度的近圆形截面毛细液滴的刚性

Rigidity for capillary liquid drops of nearly circular section with constant vorticity

作者: Pietro Baldi, Domenico Angelo La Manna, Giuseppe La Scala

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17844v1

摘要翻译: 我们考虑具有毛细作用的三维液滴的自由边界 Euler 方程在时间无关情形下、具有常涡度的解。最近已有研究利用变分方法证明了该问题解的一个刚性结果：若某个涉及涡度参数、毛细系数以及液滴赤道截面面积的量低于某一阈值，则解必然具有柱对称性，液滴形状为扁椭球体（在两极处压扁、在赤道处隆起），且每个流体粒子沿水平圆周轨迹以恒定角速度运动。本文中，我们对赤道截面在 $C^2$ 范数意义下接近于圆盘的流体区域开展扰动分析，并证明：即使在变分方法所得阈值之上，刚性结果依然成立。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压理想流体在表面张力作用下的稳态自由边界问题，其物理模型为一个处于旋转平衡态的毛细液滴（capillary liquid drop）。该系统由欧拉方程、不可压条件、恒定涡度约束、动态接触角边界条件（即压力等于边界平均曲率乘以毛细系数）以及无穿透边界条件共同刻画。具体地，设 $\Omega \subset \mathbb{R}^3$ 是一个有界开集，代表液滴占据的区域；$u : \Omega \to \mathbb{R}^3$ 是速度场；$p : \Omega \to \mathbb{R}$ 是压力；$\alpha_0 > 0$ 是恒定涡度参数；$\sigma_0 > 0$ 是毛细系数；$e_3 = (0,0,1)$；$H_{\partial\Omega}$ 是 $\partial\Omega$ 的平均曲率；$\nu_{\partial\Omega}$ 是外单位法向量。所考虑的时间无关解满足如下稳态自由边界欧拉系统：

由于 $\operatorname{curl} u = \alpha_0 e_3$，速度场必为水平方向且仅依赖于水平变量（见[4]），故可将问题降维至二维截面：令 $D := \Omega \cap \{x_3 = 0\} \subset \mathbb{R}^2$ 为赤道截面（equatorial section），则 $\Omega$ 关于平面 $x_3 = 0$ 对称，且上半边界 $\partial\Omega \cap \{x_3 > 0\}$ 可表示为图 $\{(x_1,x_2,f(x_1,x_2)) : (x_1,x_2) \in D\}$，其中 $f > 0$。进一步，$D$ 上的扭转函数 $v_D$（即 $\Delta v = -1$ in $D$, $v = 0$ on $\partial D$）控制着几何结构：文献[4]证明了 $f$ 在 $v_D$ 的等值面上为常数，并导出关于 $D$ 的超定边值问题：

其中 $H_{\partial D}$ 是平面曲线 $\partial D$ 的曲率，常数 $\beta,c$ 由 $D$ 的几何量显式表达：

其中 $|D|$ 为面积，$P(D)$ 为周长，而 $Q, I_R, I_C, I_H$ 是定义在 $D$ 上的形状泛函：

将(1.7)代入(1.6)并整理，可得关键恒等式：

其中 $A(D) := (P - I_H) I_R - P^2$, $B(D) := I_R (3Q - I_C) - 2P Q$。该恒等式是整个分析的核心——它将原三维流体力学问题完全转化为一个纯二维几何变分约束。若 $D$ 是圆盘，则 $v_D(r) = (R^2 - r^2)/4$，$|\nabla v_D| = R/2$ 为常数，$H_{\partial D} = 1/R$，所有泛函均可显式计算，从而验证(1.15)成立。反之，当 $D$ 偏离圆盘时，(1.15)成为对 $D$ 形状的强限制。此前，文献[4]利用Serrin型超定问题理论证明：若参数满足

则 $D$ 必为圆盘，从而 $\Omega$ 为扁球体（oblate spheroid），$u(x) = \frac{1}{2} \alpha_0 (-x_2, x_1, 0)$，即刚性结论成立。该阈值源于要求(1.6)中关于 $|\nabla v|$ 与 $H_{\partial D}$ 的关系满足Serrin单调性条件。然而，该条件是否最优？能否在更宽泛的参数范围内建立刚性？这构成了本文的研究动机：我们不采用变分极小化或比较原理，而是发展一套精细的扰动分析框架，处理 $\partial D$ 在 $C^2$ 拓扑下接近单位圆的情形，目标是突破(1.10)的限制，揭示刚性现象的实际适用范围。

2. 主要定理或主要结果

本文的核心结果是定理1.1（近似圆形截面解的刚性）：给定 $\alpha_0, \sigma_0 > 0$，存在仅依赖于 $\alpha_0, \sigma_0$ 的 $\delta > 0$，使得对任意满足以下条件的 $C^2$ 解 $(u, \Omega)$：(i) 赤道截面 $D$ 的质心位于 $\mathbb{R}^2$ 原点；(ii) $\partial D$ 可表示为单位圆 $S^1$ 上的图 $\partial D = \{ r x (1 + h(x)) : x \in S^1 \}$，其中 $r = (|D|/\pi)^{1/2}$，且 $\|h\|_{C^2(S^1)} < \delta$；(iii) 参数满足

则必有 $h \equiv 0$，即 $D$ 是圆盘；进而 $\Omega$ 是扁球体，其轮廓函数 $f$ 由(1.11)给出：

且速度场为 $u(x) = \frac{1}{2} \alpha_0 (-x_2, x_1, 0)$。该结论的关键技术细节在于：(1.14)中的阈值 $64/3 \approx 21.33$ 远大于(1.10)中的 $\sqrt{2}\, \pi^{3/2}/3 \approx 2.61$，表明变分方法所得阈值非最优；证明中构造的泛函 $C(D') := A(D') + \lambda_0 B(D')$，其中 $\lambda_0 = \frac{\alpha_0^2}{2\sigma_0} \left( \frac{|D|}{\pi} \right)^{3/2}$，其二阶Taylor展开在 $h=0$ 处为

其中 $h_{\ell,m}$ 是 $h$ 在球谐函数基 $Y_{\ell,m}$ 下的Fourier系数。因此，当 $\lambda_0 < 32/3$（即(1.14)）时，所有二次项系数为正，且存在常数 $\mu_0 > 0$ 使得 $T(h,C) \geq \mu_0 \|a\|_{H^{3/2}(S^1)}^2$，其中 $a$ 是 $h$ 中 $\ell \geq 2$ 的模态分量。结合余项估计 $|C(D) - T(h,C)| \leq C \|h\|_{C^2}^3$ 及插值不等式 $\|a\|_{C^2} \leq C \|a\|_{H^{3/2}}^{2/3} \|h\|_{C^2}^{1/3}$ 和 $\|b\|_{C^2} \leq C' \|a\|_{L^2}^2$（$b$ 含 $\ell = 0,1$ 模态，由面积与质心归一化强制为高阶小量），最终得到二次主导下界 $C(D) \geq \frac{\mu_0}{2} \|a\|_{H^{3/2}}^2$。故 $C(D) = 0$ 推出 $a = 0$，再结合 $b=0$ 得 $h=0$。

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

本文证明的核心障碍在于：如何从一个高度非线性的超定几何约束（1.15）中精确提取出对形状扰动 $h$ 的定量控制。传统变分方法依赖于Serrin单调性，但该性质在参数超出阈值后失效；而直接对（1.6）做线性化会丢失关键的非线性耦合信息。作者的真正创新点在于构建了一套完整的、基于形状导数（shape derivative）与Dirichlet-Neumann映射的解析扰动理论框架。该框架包含四个相互支撑的技术支柱：

第一，解析性与高阶正则性提升。引理2.1通过在（1.6）基础上迭代Schauder估计，证明了只要 $\partial D$ 是 $C^2$ 的，其边界就自动是 $C^\infty$ 的，且 $\|h\|_{C^6} \leq C \|h\|_{C^2}$。这一先验估计至关重要，它允许我们在后续所有分析中自由使用高阶Sobolev范数，为解析性论证铺平道路。

第二，Dirichlet-Neumann映射的精确公式化。引理2.2给出了扭转函数梯度模长 $\phi(h) = |(\nabla v) \circ \gamma|$ 的显式表达式：$\phi(h) = (1+h)^2/(2J) - G(h)\psi$，其中 $\psi = (1+h)^2/4$，$G(h)$ 是依赖于 $h$ 的Dirichlet-Neumann算子，$J$ 是弧长元素。此公式将 $\phi(h)$ 完全转化为 $h$ 和 $G(h)$ 的代数运算，从而将几何问题转化为算子分析问题。

第三，形状导数的系统计算与Taylor展开。引理2.4–2.6通过参数化变形 $D_\varepsilon$（边界为 $\gamma_\varepsilon(x) = x(1+\varepsilon h(x))$）并计算 $g(\varepsilon) = (\nabla v_\varepsilon) \circ \gamma_\varepsilon$ 的二阶导数，得到了 $\phi^2$ 的二阶Taylor多项式（2.21）及其一般幂次 $\phi^q$ 的展开（2.22）。该计算避开了复杂的坐标变换，直接利用 $v_0(x) = (1-|x|^2)/4$ 的显式形式和调和延拓性质，确保了所有系数的精确性。

第四，泛函的谱分解与模态能量估计。引理2.7–2.8将所有形状泛函 $A, B$ 的二阶Taylor展开投影到球谐基上，关键步骤是利用 $G Y_{\ell,m} = \ell Y_{\ell,m}$ 和 $-\Delta_{S^1} Y_{\ell,m} = \ell^2 Y_{\ell,m}$ 将积分转化为模态系数的代数运算，最终得到（2.37）。此谱分解揭示了刚性机制的本质：低阶模态（$\ell=0,1$）被面积和质心条件强制为高阶小量，而高阶模态（$\ell \geq 2$）的能量由系数 $\ell(\ell-1)(\ell+2-3\lambda_0/8)$ 控制；当 $\lambda_0 < 32/3$ 时，该系数对所有 $\ell \geq 2$ 严格为正，从而构成一个正定二次型。整个证明框架的精妙之处在于，它没有引入任何新的不等式或估计，而是通过对原始方程进行极致的解析展开与精确计算，将刚性问题彻底还原为一个可计算、可验证的代数判据。

4. 与之前工作的比较

本文与既有工作的比较需从三个维度展开：假设强度、适用范围与技术本质。

假设强度方面：文献[4]的变分方法要求（1.10）成立，其推导依赖于将（1.6）重写为 $|\nabla v| = \varphi(H_{\partial D})$ 并要求 $\varphi$ 单调，这本质上是一个全局凸性条件。而本文仅要求 $\partial D$ 在 $C^2$ 拓扑下接近圆盘（即 $\|h\|_{C^2} < \delta$），且该 $\delta$ 不依赖于 $\alpha_0, \sigma_0$ 的具体值，仅依赖于其比值 $\alpha_0^2/\sigma_0$。这是一个典型的局部假设，显著弱于[4]的全局参数阈值。此外，[4]需要 $\Omega$ 为凸集以启动bootstrap，而本文的正则性引理2.1虽也引用[4]，但其扰动框架本身对凸性无额外要求——因为 $C^2$ 小扰动自动保持凸性，故该假设在本文中是‘免费’获得的。

适用范围方面：[4]的阈值（1.10）约为2.61，而本文的阈值（1.14）为64/3≈21.33，提升了近8倍。这意味着在 $\alpha_0^2/\sigma_0$ 较大（即涡度效应相对于表面张力占优）的物理情形下，刚性结论依然成立，这极大地扩展了Plateau实验中观察到的稳定旋转液滴的理论解释范围。值得注意的是，本文并未声称（1.14）是最优阈值；它只是证明了变分方法的阈值远非最优，且存在一个更大的邻域。从提供内容中无法确认是否存在绝对最优阈值，或该阈值是否可达。

技术本质方面：[4]的方法属于经典椭圆PDE的定性理论（Serrin问题、Hopf引理、最大值原理），其优势在于简洁和普适，但劣势在于难以量化。本文则属于现代几何分析与微局部分析的交叉领域：它将自由边界问题转化为一个定义在无穷维Banach空间（$C^2(S^1)$）上的非线性方程 $C(h) = 0$，并运用解析映射理论、形状导数、Dirichlet-Neumann算子的解析依赖性（引理2.3）、以及球谐分析来研究其零点集的局部结构。这是一种典型的‘隐函数定理’思路，其代价是技术复杂度剧增，且结论仅适用于小扰动；但其回报是获得了前所未有的定量精度——不仅给出了新的阈值，还揭示了刚性失效的精确机制：当 $\lambda_0 > 32/3$ 时，（2.37）中会出现负系数的模态（例如 $\ell = 2$ 时，$2(2-1)(2+2-3\lambda_0/8) < 0$），意味着存在非平凡的‘失稳模态’，这为未来研究分岔（bifurcation）和非轴对称解提供了明确的数学线索。作者未显式展开与[13,14,25]等工作的比较，但从方法论看，本文的扰动框架与Wente[25]对旋转液滴的变分构造、Lopez[14]的能量最小化有本质不同：后者关注解的存在性，而本文聚焦于解的唯一性（刚性）；前者是全局方法，后者是局部方法。

4. 不可压 Euler 方程在端点 Sobolev 空间中的整体适定性

Global well-posedness for the incompressible Euler equations in an endpoint Sobolev space

作者: Raphaël Danchin, In-Jee Jeong

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17110v1

摘要翻译: 我们考虑二维和三维情形下涡量方程在端点临界 Sobolev 空间 $W^{d,1}(\bR^{d})$ 中的初值问题，其中 $d = 2, 3$。在二维情形下，我们证明了涡量的 $W^{2,1}(\bR^{2})$ 正则性具有整体传播性。在三维情形下，对无旋轴对称流，我们证明了涡量的 $W^{3,1}(\bR^{3})$ 正则性对所有时间均成立。这与已知的在临界 Sobolev 空间 $W^{d/p,p}(\bR^{d})$（对所有 $1 < p < \ ty$）中成立的强不适定性结果形成鲜明对比；这些不适定性结果在 $d = 3$ 时正是基于无旋轴对称流构造的。

Comment: 16 pages

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

1. 研究对象和背景

本文研究 $d = 2,3$ 维欧氏空间中不可压 Euler 方程的涡量形式初值问题。在二维情形下，速度场 $u(t,x) \in \mathbb{R}^2$ 满足 $\mathrm{div}\, u = 0$，其标量涡量 $\omega(t,x) : \mathbb{R}^2 \to \mathbb{R}$ 满足纯输运方程

；在三维情形下，涡量 $\omega(t,x) : \mathbb{R}^3 \to \mathbb{R}^3$ 满足涡拉伸方程

，其中速度由 Biot–Savart 律 $u = \nabla \times \Delta^{-1} \omega$ 唯一确定。该系统等价于经典速度形式

关键未知量是涡量 $\omega$ 的正则性演化，尤其关注其在临界函数空间中的长时间行为。所谓“临界”指与 Lipschitz 类 $C^{0,1}(\mathbb{R}^d)$ 具有相同缩放不变性的空间：对 $1 \le p \le \infty$，Sobolev 空间 $W^{d/p,p}(\mathbb{R}^d)$ 是临界的，因为 $\|\omega\|_{W^{d/p,p}}$ 在缩放 $\omega_\lambda(x) = \lambda^{d/p} \omega(\lambda x)$ 下不变，且嵌入 $W^{d/p,p} \hookrightarrow C^{0,1}$ 连续但非紧。当 $1 < p < \infty$ 时，Bourgain–Li (2015) 证明了强不适定性：存在初始涡量 $\omega_0 \in W^{d/p,p}(\mathbb{R}^d)$（具足够衰减），使得任何解 $\omega \notin L^\infty(0,t; W^{d/p,p})$ 对任意 $t > 0$ 成立；这源于可构造初值使对应速度梯度无界，从而在短时间内剧烈拉伸涡量。而在端点 $p = \infty$，Yudovich 定理给出二维全局适定性（$\omega_0 \in L^\infty \cap L^1$），但三维中因拉伸项导致强不适定（$\omega_0 \in L^\infty \cap L^1$ 不足以控制 $\|\nabla u\|_{L^\infty}$）。

本文聚焦另一端点 $p = 1$，即空间 $W^{d,1}(\mathbb{R}^d)$。该空间具有根本区别：$\omega \in W^{d,1}(\mathbb{R}^d)$ 蕴含 $u \in C^{0,1}(\mathbb{R}^d)$（由 Sobolev 嵌入 $W^{1,1} \hookrightarrow L^{d/(d-1),1} \hookrightarrow L^\infty$ 及 Calderón–Zygmund 理论保证），故 Bourgain–Li 的构造失效——其核心依赖于初值使 $\|\nabla u\|_{L^\infty} = \infty$。然而，$W^{d,1}$ 的技术难度极高：需控制 $d$ 阶导数 $\nabla^d \omega \in L^1$，而 $L^1$ 缺乏弱紧性、乘子估计失效、且标准插值不适用。此前仅 Cozzi–Harrison (2024) 在二维证明了局部传播 $W^{2,1}$ 正则性，但全局性未证；三维 $W^{3,1}$ 更无结果。本文旨在填补这一空白：证明二维全局传播 $W^{2,1}$，及三维轴对称无旋情形下全局传播 $W^{3,1}$，从而确立端点 Sobolev 空间中首个非平凡的整体适定性结果。

2. 主要定理或主要结果

定理 1.1（二维全局适定性）：设初始涡量 $\omega_0 \in W^{2,1}(\mathbb{R}^2)$，则方程 (1.2) 存在唯一全局解 $\omega \in C(\mathbb{R}_+; W^{2,1}(\mathbb{R}^2))$，且满足

存在普适常数 $C > 0$，使得对所有 $t \ge 0$，

定理 1.2（三维轴对称无旋全局适定性）：设 $\omega_0 \in W^{3,1}(\mathbb{R}^3)$ 具轴对称无旋结构，即在柱坐标下 $u = u_r e_r + u_z e_z$，$\omega = \omega_\theta e_\theta$，且 $\alpha_0 := \omega_\theta / r \in L^{3,1}(\mathbb{R}^3)$。则方程 (1.1) 存在唯一全局解 $\omega$ 满足

且 $\|\alpha(t)\|_{L^{3,1}} = \|\alpha_0\|_{L^{3,1}}$ 恒定。进一步，

其中 $C_0$ 仅依赖于 $\|\omega_0\|_{W^{3,1}}$ 和 $\|\alpha_0\|_{L^{3,1}}$。

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

证明的核心障碍在于 $L^1$-型范数的脆弱性：$L^1$ 非自反、缺乏弱*紧性、标准乘子定理（如 Mikhlin）失效、且 $L^1$ 上的微分算子无界。作者突破性地构建了一个多层嵌套的 Lorentz–Besov 分析框架，而非直接处理 $W^{d,1}$。

第一层：利用 $W^{d,1}(\mathbb{R}^d) \hookrightarrow \dot{B}^0_{\infty,1}(\mathbb{R}^d)$ 的连续嵌入（Prop. A.2），将问题降维至 $\dot{B}^0_{\infty,1}$。此空间具有关键性质：$\|f\|_{\dot{B}^0_{\infty,1}} \approx \|\nabla f\|_{\dot{B}^{-1}_{\infty,1}} \lesssim \|\nabla f\|_{L^{d,1}}$（式 (2.1)），且 Vishik 的结果保证 $\dot{B}^0_{\infty,1}$ 正则性在二维全局传播。由此获得 $\|\omega(t)\|_{\dot{B}^0_{\infty,1}}$ 的双指数增长控制（式 (2.6) 及 Gronwall）。

第二层：引入 Lorentz 空间 $L^{p,q}$ 作为中间桥梁。在二维，关键估计是 $\|\nabla \omega(t)\|_{L^{2,1}} \le \|\nabla \omega_0\|_{L^{2,1}} + \int_0^t \|\nabla u\|_{L^\infty} \|\nabla \omega\|_{L^{2,1}} d\tau$（式 (2.3)），结合 $\|\nabla u\|_{L^\infty} \lesssim \|\omega\|_{\dot{B}^0_{\infty,1}}$ 导出双指数增长。在三维轴对称情形，定义 $\alpha = \omega_\theta / r$，其满足纯输运方程 $\partial_t \alpha + u \cdot \nabla \alpha = 0$（式 (1.6)），故 $\|\alpha(t)\|_{L^{3,1}}$ 守恒。此守恒律是全局控制的基石，因为 $\|\alpha_0\|_{L^{3,1}} \lesssim \|\nabla \omega_0\|_{L^{3,1}} \lesssim \|\nabla^3 \omega_0\|_{L^1}$，且 $\|r^{-1} u_r\|_{L^\infty} \lesssim \|\alpha\|_{L^{3,1}}$（式 (3.1)）。

第三层：逐阶导数的精细估计。对 $\nabla \omega$，应用 $\nabla^\perp$ 于 (1.5) 得 $\partial_t \nabla^\perp \omega_\theta + u \cdot \nabla \nabla^\perp \omega_\theta = r^{-1} u_r \nabla^\perp \omega_\theta + \nabla^\perp(r^{-1} u_r) \omega_\theta$，其中右端首项由 (3.1) 控制，次项通过分解 $\nabla^\perp(r^{-1} u_r) \omega_\theta = -r^{-1} u_r \alpha + \nabla^\perp u_r \alpha$ 并利用 $\|\nabla u\|_{L^\infty} \lesssim \|\nabla \omega_\theta\|_{L^{3,1}} + \|\alpha_0\|_{L^{3,1}}$（链式不等式）控制，最终导出单指数增长（式 (3.4)）。对 $\nabla^2 \omega$，需估计 $\nabla^\perp \alpha$ 和 $\nabla^2_{r,z} \omega_\theta$，其方程含 $\partial u \cdot \nabla \nabla^\perp \alpha$ 项，系数 $\|\nabla u\|_{L^\infty}$ 已单指数有界，故得双指数增长。对最棘手的 $\nabla^3 \omega$，作者采用分量精细化策略：对 $\partial^3_{zzz} \omega_\theta$，方程右端含 $\alpha \partial^3_{zzz} u_r$ 等七项，其中 $\|\alpha\|_{L^3} \|\partial^3_{zzz} u_r\|_{L^{3/2}}$ 等项均被前两阶估计控制；对 $\partial^3_{rzz} \omega_\theta$，关键项 $\omega_\theta \partial^3_{rzz} v$（$v = r^{-1} u_r$）通过恒等式 $r \partial^k_r (r^{-1} f) = \partial^k_r f - k \partial^{k-1}_r (r^{-1} f)$（式 (3.7)）化为 $\alpha \partial^3_{rzz} u_r - \alpha \partial^2_{zz} v$，从而落入双指数框架。整个框架依赖于附录中建立的轴对称坐标下导数等价关系（Lemma B.1）及 Lorentz 空间嵌入（Prop. A.1, A.3）。

4. 与之前工作的比较

本文与既有工作存在本质差异。与 Bourgain–Li (2015) 的强不适定性结果对比：其针对 $W^{d/p,p}$（$1<p<\infty$）构造初值使 $\|\nabla u\|_{L^\infty} = \infty$，从而触发涡量拉伸爆炸；而 $W^{d,1}$ 中 $u \in C^{0,1}$ 强制成立，故其构造完全失效。本文首次证明端点 $p=1$ 存在全局解，构成对不适定性图景的关键补全。

与 Yudovich 定理 ($L^\infty$ 涡量) 比较：Yudovich 给出二维全局解但仅 $L^\infty$ 正则性，无法控制导数；本文在 $W^{2,1}$ 中获得 $\nabla^2 \omega \in L^1$ 的显式双指数增长，提供更高阶正则性信息。与 Vishik ($B^{2/p}_{p,1}$) 结果相比：Vishik 在 Besov 空间中证明全局适定性，但 $W^{2,1}(\mathbb{R}^2) \not\subset B^{2/p}_{p,1}$ 对 $p>1$，且 $B^{0}_{\infty,1}$ 是 $W^{2,1}$ 的真子空间，故本文结果不被其覆盖。

与 Cozzi–Harrison (2024) 的局部结果对比：其证明 $W^{2,1}$ 在二维局部传播，但未解决爆破准则；本文通过建立 $\|\omega(t)\|_{W^{2,1}}$ 的先验估计（式 (2.5)）并结合 $\|\nabla u\|_{L^\infty} \lesssim \|\omega\|_{\dot{B}^0_{\infty,1}}$ 的双指数控制，首次实现全局化。代价是增长率更差（双指数 vs 单指数），但这是 $L^1$ 估计固有的最优可能。

与三维轴对称已有工作比较：Ukhovskiǐ–Yudovich 利用 $\alpha$ 输运证明 $L^\infty$ 全局存在；Cheskidov–Friedlander 等在 $L^p$ ($p>3$) 中建立适定性；Bang–Cheskidov 证明 $\alpha_0 \in L^{3,q}$ ($q>1$) 导致强不适定，本文取 $q=1$ 达到 sharpness。本文首次将正则性提升至 $W^{3,1}$，且明确给出 $\nabla^3 \omega \in L^1$ 的双指数增长，而此前工作仅处理低阶导数或 $L^p$ 空间。

局限性方面：本文结果严格依赖轴对称无旋结构，一般三维 $W^{3,1}$ 仍开放（第 4 页末段明确指出）；二维结果虽全局，但 $W^{2,1}$ 是否为最大临界空间未知；所有增长率（尤其双指数）是否可改进亦未解决。此外，证明中使用 Lorentz 空间作为中介，暗示 $W^{d,1}$ 本身可能非最优框架，未来或需发展 $L^1$-型 Littlewood–Paley 理论。

5. 基于累积分布变换的守恒型一维偏微分方程降阶建模

Reduced Order Modeling of One-Dimensional Conservative PDEs via the Cumulative Distribution Transform

作者: Harbir Antil, Rocío Díaz Martín, Ivan V. Medri, Kristofor E. Pas, Gustavo K. Rohde, Aryan Saxena, Sarswati Shah

类别: math.AP, math.NA

发布日期: 2026-07-19

arXiv ID: 2607.17066v1

摘要翻译: 我们提出一种基于累积分布变换（Cumulative Distribution Transform, CDT）的一维守恒型偏微分方程（PDE）降阶建模（Reduced Order Modeling, ROM）框架。CDT 将非负、等质量的状态映射到一个 Hilbert 空间，其中一维 Wasserstein 距离转化为加权 $L^2$ 距离，而平移则变为仿射平移。这使得该变换特别适用于以输运为主导的动力学系统，而在该类系统中，欧拉坐标下的线性子空间 ROM 常因 Kolmogorov 宽度衰减缓慢而失效。

Comment: 无

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

1. 研究对象和背景

本文研究一类一维标量守恒型偏微分方程（PDE）的降阶建模（Reduced Order Modeling, ROM）问题，核心对象为满足守恒律的非负、等质量演化密度函数 $ u = u(x,t) : (a,b) \times [0,T] \to [0,\infty) $，其动力学由如下初值问题刻画：

其中通量 $ f[u](x,t) $ 依赖于 $ u $ 及其导数（如 $ f[u] = A u - D u_x $ 对应守恒型对流-扩散方程；$ f[u] = \frac{1}{2}u^2 $ 对应无粘Burgers方程）。关键假设是：所有快照 $ u(\cdot,t) $ 均为非负概率密度（即 $ \int_a^b u(x,t)\,dx = M > 0 $ 对所有 $ t \in [0,T] $ 成立），这由零通量边界条件 $ f[u](a,t) = f[u](b,t) = 0 $ 或无穷远处充分衰减所保证。该设定使累积分布变换（Cumulative Distribution Transform, CDT）可严格定义——CDT 将每个快照 $ u(\cdot,t) $ 映射为一个单调递增的传输映射 $ \tilde{u}(\cdot,t) : K \to I $，其定义依赖于一个固定参考密度 $ r : K \to [0,\infty) $，满足

等价地，$ \tilde{u}(\xi,t) = F_u^{-1}(F_r(\xi)) $，其中 $ F_u(x) = \int_{-\infty}^x u(y,t)\,dy $ 是 $ u $ 的累积分布函数（CDF），$ F_u^{-1} $ 是其广义逆（quantile function）。此映射在最优传输理论中即为从参考测度 $ r \,dx $ 到当前测度 $ u(\cdot,t)\,dx $ 的唯一 $ W_2 $-最优传输映射。研究动机源于传统欧拉坐标系下线性子空间ROM（如POD）的根本局限：对于以平移、传播为主的守恒型动力学（如线性输运 $ u_t + A u_x = 0 $、激波解、稀疏波），其解流形 $ \mathcal{M}_{[0,T]} = \{ u(\cdot,t) : t \in [0,T] \} $ 在 $ L^2 $ 空间中 Kolmogorov $ n $-宽度衰减极慢（例如对带跳跃的平移族，$ d_n(\mathcal{M}_{[0,T]}) \sim n^{-1/2} $），导致POD需大量模态才能达到精度要求。这一“Kolmogorov障碍”源于欧拉坐标中移动结构在时间上严重失准。已有对策包括非线性流形、自适应基、Lagrangian坐标或显式传输图学习，但均面临构造复杂、依赖先验知识或难以泛化等问题。本文试图填补的空白在于：系统建立CDT作为通用、数据驱动、解析可处理的Lagrangian坐标框架，并严格证明其如何将守恒型PDE的解流形在变换空间中“拉直”为低维（甚至精确二维）结构，从而绕过Kolmogorov障碍。特别地，需回答：CDT是否能将激波后、非光滑、含扩散的保守演化也映射为具有良好逼近性质的轨迹？其理论基础是CDT的两个核心几何性质：(i) 它是 $ W_2 $ 度量到加权 $ L^2 $ 度量的等距同构，即 $ W_2(u^{(1)}, u^{(2)}) = \| \tilde{u}^{(1)} - \tilde{u}^{(2)} \|_{L^2(r)} $；(ii) 平移 $ u(x,t) \mapsto u(x - At,t) $ 在CDT空间中变为仿射平移 $ \tilde{u}(\xi,t) \mapsto \tilde{u}(\xi,t) + At $。因此，CDT空间中的解流形 $ \widetilde{\mathcal{M}}_{[0,T]} = \{ \tilde{u}(\cdot,t) : t \in [0,T] \} $ 不再是扭曲的欧拉曲线，而是具有明确正则性的Hilbert空间曲线，其逼近难度由 $ \widetilde{\mathcal{M}}_{[0,T]} $ 的Kolmogorov宽度决定，而后者可被PDE本身的守恒结构与正则性所控制。

2. 主要定理或主要结果

本文主要结果分为三类，均围绕CDT空间 $ \widetilde{\mathcal{M}}_{[0,T]} $ 的Kolmogorov $ n $-宽度 $ d_n(\widetilde{\mathcal{M}}_{[0,T]}) $ 展开，该宽度定义为 $ d_n(\widetilde{\mathcal{M}}_{[0,T]}) = \inf_{\substack{V_n \subset L^2(r) \\ \dim V_n = n}} \sup_{t \in [0,T]} \inf_{v \in V_n} \| \tilde{u}(\cdot,t) - v \|_{L^2(r)} $。

1. 线性输运方程（定理未编号，见第2.3节及引言）：对 $ u_t + A u_x = 0 $，$ u(x,0) = u_0(x) $，若 $ u_0 \geq 0 $ 且 $ \int u_0 = \int r $，则CDT轨迹满足 $ \tilde{u}(\xi,t) = \tilde{u}_0(\xi) + A t $，其中 $ \tilde{u}_0 = \text{CDT}(u_0) $。因此 $ \widetilde{\mathcal{M}}_{[0,T]} \subset \text{span}\{ \tilde{u}_0, \, 1 \} $，即完全包含于二维子空间，故 $ d_2(\widetilde{\mathcal{M}}_{[0,T]}) = 0 $。

2. 非线性双曲守恒律（定理4.2与4.5）：对熵解 $ u $ 满足 $ u_t + (f(u))_x = 0 $，$ u(x,0) = u_0(x) \geq 0 $，$ \int u_0 = \int r $。

   • 鲁棒估计（定理4.2）：仅依赖守恒结构与CDT-Wasserstein等距性，不假设光滑性或特征可逆。证明 $ \tilde{u}(\cdot,t) $ 关于 $ t $ 是全局Lipschitz连续的，即 $ \| \tilde{u}(\cdot,t_1) - \tilde{u}(\cdot,t_2) \|_{L^2(r)} \leq C |t_1 - t_2| $，其中常数 $ C $ 仅依赖于 $ f' $ 的界和 $ u_0 $ 的总变差。由此推出 $ d_n(\widetilde{\mathcal{M}}_{[0,T]}) \leq C n^{-1} $，该估计对任意 $ t \in [0,T] $ 成立，即使激波已形成。
   • 高阶估计（定理4.5）：在激波前光滑区域，若 $ u_0 \in C^2 $ 且 $ u_0 > 0 $，则 $ \tilde{u}(\xi,t) = \xi + t f'(u_0(\xi)) + O(t^2) $，且 $ \partial_t^2 \tilde{u} \in L^2(r) $。由此得到 sharper 估计 $ d_n(\widetilde{\mathcal{M}}_{[0,T]}) \leq C n^{-2} $。

3. 守恒型对流-扩散方程（定理5.1与5.6）：对 $ u_t + (A u - D u_x)_x = 0 $，$ D > 0 $，$ u(x,0) = u_0(x) \geq 0 $，$ \int u_0 = \int r $。记 $ u^{\text{adv}} $ 为对应纯对流方程 $ u_t + A u_x = 0 $ 的解，则其CDT为 $ \tilde{u}^{\text{adv}}(\xi,t) = \tilde{u}_0(\xi) + A t $。

   • 鲁棒估计（定理5.1）：证明 $ \| \tilde{u}(\cdot,t) - \tilde{u}^{\text{adv}}(\cdot,t) \|_{L^2(r)} \leq C \sqrt{D T} $，即CDT轨迹始终位于纯输运平面 $ \text{span}\{ \tilde{u}_0, \, 1 \} $ 的 $ O(\sqrt{D T}) $ 邻域内。故 $ d_2(\widetilde{\mathcal{M}}_{[0,T]}) \leq C \sqrt{D T} $，且当 $ D \to 0 $ 时恢复 $ d_2 = 0 $。
   • 高阶估计（定理5.6）：若 $ u_0 \in H^2 $ 或远离初始层（$ t \gg \sqrt{D} $），则有 $ \| \tilde{u}(\cdot,t) - \tilde{u}^{\text{adv}}(\cdot,t) \|_{L^2(r)} \leq C D^2 T^2 $，故 $ d_2(\widetilde{\mathcal{M}}_{[0,T]}) \leq C D^2 T^2 $。

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

证明的核心障碍在于：CDT空间中的解轨迹 $ \tilde{u}(\xi,t) $ 并非直接由PDE演化方程所描述，而是一个隐式定义的传输映射，其正则性（如Lipschitz性、二阶可微性）无法从原始PDE的欧拉形式直接导出。作者的真正创新点在于构建了一个双重分析框架，将CDT轨迹的逼近性质与PDE的物理结构（守恒、输运、扩散）精确耦合。

第一重工具：Wasserstein几何与守恒结构的鲁棒估计（对应定理4.2与5.1）。这是全文最根本的创新。作者并未尝试对 $ \tilde{u} $ 求导，而是利用CDT的等距性 $ W_2(u^{(1)}, u^{(2)}) = \| \tilde{u}^{(1)} - \tilde{u}^{(2)} \|_{L^2(r)} $，将问题转化为对原始解 $ u $ 的 $ W_2 $ 距离估计。对于双曲方程，利用熵解的$ L^1 $-稳定性与 $ W_2 $ 距离的收缩性（由通量函数的Lipschitz性保证），直接得到 $ W_2(u(\cdot,t_1), u(\cdot,t_2)) \leq C |t_1 - t_2| $，从而推出 $ \| \tilde{u}(\cdot,t_1) - \tilde{u}(\cdot,t_2) \|_{L^2(r)} \leq C |t_1 - t_2| $。对于对流-扩散方程，作者构造了一个精巧的比较原理：将真实解 $ u $ 与纯对流解 $ u^{\text{adv}} $ 的 $ W_2 $ 距离分解为由扩散引起的“质量位移”项，并利用扩散算子的耗散性（通过能量估计或Benamou-Brenier公式）证明该位移被 $ O(\sqrt{D T}) $ 控制。此方法完全规避了对 $ \tilde{u} $ 的微分运算，使其在激波、不连续性等奇异性下依然有效。

第二重工具：半Lagrangian坐标下的正则性提升（对应定理4.5与5.6）。当解光滑时，作者将CDT轨迹 $ \tilde{u}(\xi,t) $ 解释为广义特征曲线：它满足隐式关系 $ u(\tilde{u}(\xi,t), t) \, \partial_\xi \tilde{u}(\xi,t) = r(\xi) $（由(6)式对 $ \xi $ 求导得到）。在此基础上，作者推导出 $ \tilde{u} $ 所满足的形式PDE：对双曲方程，有 $ \partial_t \tilde{u}(\xi,t) = f'(u(\tilde{u}(\xi,t),t)) $；对对流-扩散方程，有更复杂的半Lagrangian形式（见第3.1节）。关键突破在于，作者证明了在光滑预激波区域，$ u $ 的正则性（$ C^2 $）可传递至 $ \tilde{u} $，即 $ \partial_t^2 \tilde{u} \in L^2(r) $。这一结论并非平凡，因为它要求验证 $ \partial_t u $ 和 $ \partial_t^2 u $ 在 $ \tilde{u} $ 的像集上的可积性，作者通过细致的链式法则与反函数定理完成。正是这一 $ L^2 $-二阶时间导数的存在性，结合Bochner积分理论中Hilbert空间值曲线的逼近定理（$ \| \tilde{u}(\cdot,t) - P_n \tilde{u}(\cdot,t) \| \lesssim n^{-2} \| \partial_t^2 \tilde{u} \|_{L^2} $），最终导出了 $ O(n^{-2}) $ 和 $ O(D^2 T^2) $ 的sharp估计。

第三重工具：几何视角的统一解释。作者将CDT空间中的解流形同时视为两条曲线：(i) Hilbert空间中的曲线 $ t \mapsto \tilde{u}(\cdot,t) \in L^2(r) $，用于应用泛函分析的逼近理论；(ii) 参数化的族 $ \{ t \mapsto \tilde{u}(\xi,t) \}_{\xi \in K} $，即每个质量标签 $ \xi $ 对应一条“累积质量轨迹”，用于进行逐点正则性分析。这种双重视角是连接鲁棒估计（整体 $ L^2 $ 性质）与sharp估计（逐点或高阶导数）的桥梁，也是区别于以往仅关注单一表示的创新所在。

4. 与之前工作的比较

本文与既有工作的比较需从四个维度展开：理论基础、适用范围、估计精度与计算代价。

理论基础维度：经典ROM理论（如[15,20]）将Kolmogorov宽度障碍归因于欧拉坐标中平移的非线性性，其对策多为构造特定传输图（如[25–31]）。本文则首次将CDT这一源自最优传输的数学变换系统性地嵌入ROM框架，并赋予其严格的PDE理论支撑。与[44,45]的Radon-CDT相比，本文的CDT-POD直接作用于一维密度，避免了Radon变换引入的额外维度和伪影；与[46–48]的参数化传输图方法相比，本文的CDT是数据自适应的、精确的、一维单调最优传输映射，而非需要预设参数族或学习过程的近似。

适用范围维度：现有Lagrangian方法（如[33–35]）通常要求方程为纯双曲或具有显式特征结构，且在激波后失效。本文的鲁棒估计（定理4.2, 5.1）首次证明CDT-POD对熵解（含激波）、守恒型对流-扩散（含粘性）乃至非双曲守恒律均有效，其理论根基是守恒律本身与Wasserstein几何，而非特征线的可逆性。例如，定理4.2的 $ O(n^{-1}) $ 估计在激波形成后依然成立，而[20,21]中针对欧拉坐标的 $ n^{-1/2} $ 障碍在此被彻底突破。此外，定理5.1中 $ \sqrt{D T} $ 的误差界表明，CDT天然兼容从双曲到抛物的连续过渡，这是多数Lagrangian方法无法实现的。

估计精度维度：在光滑区域，本文的 $ O(n^{-2}) $ 估计（定理4.5）优于标准POD的 $ O(n^{-1}) $（由Sobolev嵌入导出），也优于基于Fourier或小波的线性方法（其对平移的逼近仍为 $ O(n^{-1/2}) $）。然而，此优势是有代价的：它严格依赖于 $ u_0 \in C^2 $ 且 $ u_0 > 0 $，而经典POD的 $ O(n^{-1}) $ 估计对 $ L^2 $ 初始数据即成立。同样，定理5.6的 $ O(D^2 T^2) $ 界虽比 $ O(\sqrt{D T}) $ 更sharp，但要求 $ u_0 \in H^2 $ 或远离初始层，这在实际瞬态问题中可能不满足。因此，本文的贡献不是取代所有方法，而是提供了一个精度-鲁棒性权衡的新谱系：鲁棒估计保证底线性能，sharp估计在良好条件下提供超额收益。

计算代价维度：CDT-POD的预处理（CDT计算）比标准POD昂贵，因其涉及求解单调传输映射（即求CDF的逆），计算复杂度为 $ O(N \log N) $（$ N $ 为网格点数），而POD仅为 $ O(N) $。然而，作者指出，这一代价被后续POD维数的大幅降低（实验显示模态数减少5–10倍）所补偿。更重要的是，CDT-POD的重构（逆CDT）必须保证单调性以维持质量守恒，而线性截断POD在CDT空间中可能破坏单调性（见贡献(i)），这是一个明确的局限，也是未来工作（如约束优化）的方向。相比之下，[33–35]的Lagrangian ROM虽也具高精度，但其变形场的学习或求解往往比CDT更昂贵且不稳定。综上，本文推进的本质是：将ROM的理论基础从线性代数（SVD）拓展至最优传输几何，并证明该几何结构与守恒PDE的物理本质深度契合，从而在保持算法简洁性的同时，实现了理论深度与实用广度的统一。

6. 非等熵可压缩流体-粒子相互作用模型（含热效应）的增强稳定性与渐近极限

Enhanced stability and asymptotic limits to the non-isentropic compressible fluid-particle interaction model with thermal effects

作者: Fucai Li, Jinkai Ni, Zhouping Xin

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17115v1

摘要翻译: 在爱因斯坦的开创性工作 [Ann. Physik, 17 (1905), 549-560] 中，他指出流体的温度会显著影响悬浮粒子的运动。为更精确地描述该物理过程中温度的影响，Boudin 等人 [ESAIM Proc., 28 (2009), 195-210] 引入了一个新的流体-粒子相互作用模型，其中流体部分由非等熵可压缩 Euler 方程描述，粒子部分则由一类非线性 Vlasov-Fokker-Planck 型方程描述。Mu 和 Wang [Calc. Var. Partial Differential Equations, 59 (2020), Paper no. 110] 通过在该模型的流体部分添加粘性项与热传导项，建立了其在平衡态附近的经典解的整体存在性。

Comment: 76 pages

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

1. 研究对象和背景

本文研究的是一类耦合了宏观流体动力学与微观粒子统计行为的非线性偏微分方程系统，其核心对象为Boudin等人于2009年提出的非等熵可压缩Euler-Vlasov-Fokker-Planck（Euler-VFP）系统（见(1.12)），该模型旨在精确刻画温度场对悬浮粒子运动的物理影响——这一思想直接源于爱因斯坦1905年关于布朗运动的奠基性工作。系统由四组方程构成：

其中$(\varrho, u, \vartheta) : \mathbb{R}_+ \times \mathbb{R}^3 \to \mathbb{R}_+ \times \mathbb{R}^3 \times \mathbb{R}_+$为流体密度、速度与温度；$F = F(t,x,v) \ge 0$为粒子在相空间$\mathbb{R}^3_x \times \mathbb{R}^3_v$上的分布函数；压力与内能满足理想气体状态方程$P = R \varrho \vartheta$, $e = C_v \vartheta$（文中设$R = C_v = 1$）；耦合项$\int (v-u)F dv$与$\int(|v-u|^2 - 3\vartheta)F dv$分别表征动量与能量交换；而右端VFP型算子$\operatorname{div}_v((v-u)F + \vartheta \nabla_v F)$则同时包含与流体速度$u$相关的线性拖曳力及与温度$\vartheta$成正比的热扩散项——这显著区别于等熵情形中常系数扩散$\operatorname{div}_v((v-u)F + \nabla_v F)$。

该模型的数学挑战性根植于其双重退化结构：一方面，当$F \equiv 0$时，(1.12)退化为经典的三维非等熵可压缩Euler系统（见(1.13)），已知其光滑解在一般初值下必于有限时间发生激波奇点（参见Sideris, 1985; Christodoulou, 2007）；另一方面，尽管VFP方程在$v$变量上具有Fokker-Planck型耗散，但其在$x$方向完全双曲，无法单独提供对流体变量$(\varrho,u,\vartheta)$的任何空间正则化或衰减控制。因此，自2009年提出以来，(1.12)的全局适定性始终是开放问题。此前唯一可用的理论结果来自Mu-Wang（2020）与Li-Ni-Wu（2024）：他们在流体方程中引入粘性项$-\mu \Delta u - (\mu+\lambda)\nabla \operatorname{div} u$与热传导项$-\kappa \Delta \vartheta$，从而获得全局经典解，但其证明严重依赖这些人工耗散项所提供的$H^4$能量耗散。一个根本性问题是：在$\mu = \lambda = \kappa = 0$的纯双曲极限下，是否存在内在的、由流体-粒子耦合激发的新型耗散机制？本文即旨在填补这一关键空白——通过揭示并量化由宏观粒子速度$b = \int v \sqrt{M} f \, dv$与流体速度$u$之差、以及宏观粒子温度$\omega = \frac{1}{\sqrt{6}} \int (|v|^2 - 3) \sqrt{M} f \, dv$与流体温度$\vartheta$之差所诱导的“增强耗散”（enhanced dissipation），从而在不添加任何人工粘性/热导的前提下，建立(1.12)的全局经典解及其最优衰减率。

2. 主要定理或主要结果

本文在$\mathbb{R}^3$上对小初值扰动建立了如下严格结论：

定理1.1（一致全局适定性）：设$\mu, \lambda, \kappa > 0$，初始数据$(\rho_0^{\mu,\lambda,\kappa}, u_0^{\mu,\lambda,\kappa}, \theta_0^{\mu,\lambda,\kappa}, f_0^{\mu,\lambda,\kappa})$满足$\|\cdot\|_{H^4 \times H^4_{x,v}}^2 \le \varepsilon_0$，其中$\varepsilon_0 > 0$独立于$\mu,\lambda,\kappa$。则系统(1.14)（即(1.8)在平衡态附近的扰动形式）存在唯一全局经典解$(\rho^{\mu,\lambda,\kappa}, u^{\mu,\lambda,\kappa}, \theta^{\mu,\lambda,\kappa}, f^{\mu,\lambda,\kappa})$，且满足一致先验估计： $$ \sup_{t\ge0} \Big( | (\rho^{\mu,\lambda,\kappa}, u^{\mu,\lambda,\kappa}, \theta^{\mu,\lambda,\kappa})(t) |{H^4}^2 + | f^{\mu,\lambda,\kappa}(t) |{H^4_{x,v}}^2 \Big)

• \mu \int_0^t | \nabla u^{\mu,\lambda,\kappa} |_{H^4}^2 d\tau
• (\mu+\lambda) \int_0^t | \operatorname{div} u^{\mu,\lambda,\kappa} |_{H^4}^2 d\tau
• \kappa \int_0^t | \nabla \theta^{\mu,\lambda,\kappa} |_{H^4}^2 d\tau
• \int_0^t D(\rho^{\mu,\lambda,\kappa}, u^{\mu,\lambda,\kappa}, \theta^{\mu,\lambda,\kappa}, f^{\mu,\lambda,\kappa}) d\tau \le C_0 \varepsilon_0,
  $$ 其中耗散泛函$D$定义为： $$
  D(\cdot)(t) := | b^{\mu,\lambda,\kappa} - u^{\mu,\lambda,\kappa} |{H^4}^2 + | \sqrt{2} \omega^{\mu,\lambda,\kappa} - \sqrt{3} \theta^{\mu,\lambda,\kappa} |{H^4}^2 + | \nabla(a^{\mu,\lambda,\kappa}, b^{\mu,\lambda,\kappa}, \omega^{\mu,\lambda,\kappa}, \rho^{\mu,\lambda,\kappa}) |{H^3}^2 + \sum{|\alpha|\le4} | {I-P} \partial^\alpha f^{\mu,\lambda,\kappa} |\nu^2 + \sum{\substack{1\le|\beta|\le4 \ |\alpha|+|\beta|\le4}} | \partial^\alpha_\beta {I-P} f^{\mu,\lambda,\kappa} |_\nu^2. $$

定理1.2（零粘性-零热导极限）：当$\mu,\lambda,\kappa \to 0$时，上述解弱-*收敛于(1.16)（即(1.12)的扰动形式）的唯一全局解$(\rho,u,\theta,f)$，且满足：

定理1.3（收敛速率）：若初值误差满足$\| (\rho_0^{\mu,\lambda,\kappa}-\rho_0, u_0^{\mu,\lambda,\kappa}-u_0, \theta_0^{\mu,\lambda,\kappa}-\theta_0) \|_{H^1} + \| f_0^{\mu,\lambda,\kappa}-f_0 \|_{H^1_{x,v}} \le \max\{\mu,\lambda,\kappa\}$，则解误差满足： $$ \sup_{t\ge0} \Big( | (\rho^{\mu,\lambda,\kappa}-\rho, u^{\mu,\lambda,\kappa}-u, \theta^{\mu,\lambda,\kappa}-\theta)(t) |{H^1}^2 + | (f^{\mu,\lambda,\kappa}-f)(t) |{H^1_{x,v}}^2 \Big)

• \int_0^t \Big( | (b^{\mu,\lambda,\kappa}-u^{\mu,\lambda,\kappa}) - (b-u) |{H^1}^2 + | (\sqrt{2}\omega^{\mu,\lambda,\kappa}-\sqrt{3}\theta^{\mu,\lambda,\kappa}) - (\sqrt{2}\omega-\sqrt{3}\theta) |{H^1}^2 \Big) d\tau \le C_2 \max{\mu,\lambda,\kappa}^2. $$

定理1.4（最优衰减率）：对$\mu=\lambda=\kappa=0$的Euler-VFP系统(1.12)，其全局解满足$L^2$最优衰减：

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

证明的核心障碍在于：在$\mu=\lambda=\kappa=0$时，系统(1.12)完全丧失抛物正则化，标准能量方法无法控制非线性项的级联增长，尤其无法处理$u\cdot\nabla u$、$u\cdot\nabla\theta$等输运项导致的梯度爆炸风险。作者的真正创新点在于构建了一套基于“宏观-微观分解”与“有效耗散模式识别”的全新分析框架，其技术要点分述如下：

1. 宏观-微观分解与新型耗散泛函$D$的设计：不同于传统Navier-Stokes理论中依赖$\nabla u$、$\nabla\theta$的耗散，本文关键观察到VFP方程的线性化算子$L$（定义见(1.18)）在$v$变量上具有谱隙性质，其宏观投影$Pf = a\sqrt{M} + b\cdot v\sqrt{M} + \omega \frac{|v|^2-3}{\sqrt{6}}\sqrt{M}$（见(1.19)-(1.21)）与流体变量存在强耦合。作者定义两个核心“有效耗散模式”：$b-u$（宏观粒子动量与流体动量之差）与$\sqrt{2}\omega-\sqrt{3}\theta$（经尺度归一化的宏观粒子能量与流体内能之差）。通过精细计算(1.16)中$u$与$\theta$方程的演化，发现$\partial_t(b-u)$与$\partial_t(\sqrt{2}\omega-\sqrt{3}\theta)$均含负定线性项$-(b-u)$与$-\sqrt{3}(\sqrt{2}\omega-\sqrt{3}\theta)$，从而构成真正的阻尼机制。耗散泛函$D$正是围绕这两个模式构造，强制将高阶能量控制锚定于这些差量之上。

2. 低-高频分解与Lyapunov型不等式：为克服纯双曲系统的色散不足，作者在第5节引入低频主导的Lyapunov泛函$\mathcal{L}(t)$，其构造融合了流体变量的$\dot{H}^{-1}$范数（捕获质量守恒的长波行为）与粒子分布的微观部分$\{I-P\}f$的加权$L^2_\nu$范数（$\nu(v)=1+|v|^2$）。通过设计$\mathcal{L}(t)$满足$\frac{d}{dt}\mathcal{L}(t) + c \mathcal{L}(t) \lesssim \|\nabla u\|_{L^2}^2 + \|\nabla\theta\|_{L^2}^2$，再结合$D$中对$\nabla u$、$\nabla\theta$的控制，最终导出$\mathcal{L}(t) \lesssim (1+t)^{-3/2}$，进而推出$L^2$最优衰减率。

3. 一致先验估计的闭合技巧：为证明定理1.1中的$\mu,\lambda,\kappa$一致性，作者发展了一种“解耦式能量估计”。具体而言，在估计$\|\nabla^k u\|_{L^2}^2$时，并不直接对(1.14)_2求导（会引入$\mu\Delta\nabla^k u$项破坏一致性），而是利用$u$方程与$b-u$方程的差分结构，将$\mu\|\nabla^{k+1}u\|_{L^2}^2$项吸收进$\|b-u\|_{H^k}^2$的耗散中；同理，$\kappa\|\nabla^{k+1}\theta\|_{L^2}^2$被吸收进$\|\sqrt{2}\omega-\sqrt{3}\theta\|_{H^k}^2$。这种“以宏观差量代偿粘性/热导”的策略，是突破一致性的核心技术。

4. 与之前工作的比较

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

1. 与Mu-Wang (2020) 和 Li-Ni-Wu (2024) 的比较：前人工作[46,36]仅在$\mu,\lambda,\kappa>0$的粘性-热导框架下建立全局解，其结论本质是“带耗散的近似模型”的适定性，且收敛速率仅为$O(\mu^{1/2})$（见摘要提及的Wang [53]结果）。本文首次证明：即使在$\mu=\lambda=\kappa=0$的纯双曲极限下，Boudin模型(1.12)本身即具有内在全局稳定性。这不仅是假设的彻底放宽（从“必须有粘性”到“无需粘性”），更是物理机制认知的跃迁——揭示了粒子相与流体相的动量/能量交换本身即可产生等效耗散，其强度由$\|b-u\|_{H^4}^2$与$\|\sqrt{2}\omega-\sqrt{3}\theta\|_{H^4}^2$量化，这在纯流体Euler系统中完全不存在。

2. 与纯非等熵Euler系统(1.13)的比较：对于(1.13)，已知其经典解在$\mathbb{R}^3$上必然爆破（Sideris, CPAM 1985）。本文表明，一旦耦合粒子相，系统便获得全局解，其代价是引入更高维的相空间分析与更复杂的耗散结构。值得注意的是，本文的最优衰减率$(1+t)^{-3/4}$与无碰撞Boltzmann方程的流体极限一致，但其机制迥异：后者依赖于Boltzmann算子的谱隙，而本文依赖于VFP算子与流体方程的耦合共振。

3. 与等熵NS-VFP系统(1.2)的比较：等熵模型（如[44,13]）中VFP算子为$\operatorname{div}_v((v-u)F + \nabla_v F)$，其扩散系数恒为1，无法反映温度变化。本文的热扩散项$\vartheta \nabla_v F$使VFP算子成为变系数算子，导致宏观投影$Pf$的定义（见(1.21)）必须包含$\vartheta$-依赖的归一化，从而使得$\sqrt{2}\omega-\sqrt{3}\theta$这一耗散模式成为必要——这是非等熵特有现象，从提供内容中无法确认等熵模型中存在类似结构。此外，本文的收敛速率$O(\max\{\mu,\lambda,\kappa\})$优于等熵情形常见的$O(\mu^{1/2})$，表明热效应增强了极限过程的稳定性。

7. 环面上的投影型连续数据同化：共振与无核情形

Projectional continuous data assimilation on the torus: Resonant and kernel-free regimes

作者: Tsogtgerel Gantumur

类别: math.AP, math.DS

发布日期: 2026-07-18

arXiv ID: 2607.16827v1

摘要翻译: 我们在周期环面 $\mathbb{T}^2$ 上研究二维不可压 Navier-Stokes 方程的连续数据同化问题，观测方式为沿给定空间变方向的单个带符号标量速度投影。我们通过两种互补机制证明了 $L^2$ 和 $H^1$ 范数下的指数同步性。在常规共振情形下，非平凡的不可见流通过动坐标系展开、适配坐标以及 Farhat-Lunasin-Titi 对数估计加以控制，前提是由此产生的几何剪切缺陷可被粘性吸收。该机制恢复了常值有理方向下的周期单分量机制，并适用于真正非恒定的投影场。在无核情形下，定性单射性与紧性导出可观测性，且粘性泄漏可任意小，从而对每个足够大的增益（无需上界限制）均实现同步。对于足够正则的投影场，上述两种机制均可推广至满足一阶逼近性质的 $L^2$-稳定 Type-I 粗粒度标量观测，前提是满足通常的增益-分辨率条件。随后，一个共同的抛物正则化论证在不增加任何观测假设的前提下，将全部四个 $L^2$-同步结果提升为 $H^1$-同步结果。

Comment: 24 pages, 9 sections, 1 appendix

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

1. 研究对象和背景

本文研究二维不可压 Navier-Stokes 方程在周期环面 $ \mathbb{T}^2 = (\mathbb{R}/2\pi\mathbb{Z})^2 $ 上的连续数据同化（continuous data assimilation, CDA）问题，核心观测方式为单一分量的空间变向标量投影：给定一个可测单位向量场 $ m \in L^\infty(\mathbb{T}^2; \mathbb{R}^2) $，观测量为参考解 $ u $ 与待估解 $ v $ 的误差 $ e = v - u $ 在方向 $ m $ 上的标量投影 $ r = e \cdot m $。该模型统一涵盖了经典的一分量 nudging（如 $ m = e_1 $）与更一般的几何投影观测，其数学本质是通过低维（标量）观测驱动高维（二维向量）动力系统实现指数同步。

所研究的方程为标准周期 Navier-Stokes 系统：

其中 $ P_\sigma $ 是到零均值无散子空间 $ H = L^2_{\text{div},0}(\mathbb{T}^2) $ 的正交投影，$ B(u,v) = P_\sigma((u\cdot\nabla)v) $，双线性形式 $ b(e,u,e) = \int_{\mathbb{T}^2} (e \cdot \nabla u) \cdot e \, dx $。nudging 系统为

其误差 $ e = v - u $ 满足关键能量恒等式：

同步的核心障碍在于非线性项 $ -b(e,u,e) $ 的控制：粘性耗散 $ \nu \|e\|_{H^1}^2 $ 与观测耗散 $ \mu |r|_{L^2}^2 $ 必须共同压制该项。而 $ r = e \cdot m $ 仅提供 $ e $ 在 $ m $ 方向的信息，其“盲区”由无散核 $ K_m = \{ e \in H : e \cdot m = 0 \} $ 刻画。若 $ K_m \neq \{0\} $，则存在非零无散“隐形电流” $ e = \sigma n $（其中 $ n = m^\perp $），它们完全不被观测捕获，构成同步的根本障碍。例如，当 $ m $ 为常向量且 $ n $ 平行于有理格点 $ \ell \in \mathbb{Z}^2 \setminus \{0\} $ 时，$ e(x) = \psi(\ell^\perp \cdot x) n $ 属于 $ K_m $，即经典的周期横向剪切模态；此时系统处于共振情形（resonant regime）。反之，若 $ K_m = \{0\} $（如 $ n $ 无理方向），则称无核情形（kernel-free regime）。

已有工作（Farhat–Lunasin–Titi, FLT）证明：对常有理方向 $ m $，可通过 Navier-Stokes 非线性耦合将信息从 $ r $ 传递至隐形分量，从而实现同步，其技术核心是对数型相容估计（logarithmic compatibility estimate）：存在常数 $ C_{\text{FLT}} $，使得

但该估计依赖于坐标系的全局平直性，无法直接推广至非恒定 $ m $。本文旨在填补两大空白：（1）建立适用于真正非恒定投影场 $ m $ 的共振理论，超越 FLT 仅处理常向量的局限；（2）构建无核情形下无需显式小除子估计的同步机制，避免对 $ n $-流轨道闭合性的依赖，从而覆盖更广泛的 $ m $（如无理方向或一般光滑 $ m $）。这两个机制分别对应论文标题中的“resonant”与“kernel-free”两种互补范式，其统一框架揭示了投影观测下同步的几何本质——即由 $ m $ 的微分几何结构（曲率、剪切）与拓扑结构（核的存在性）共同决定。

2. 主要定理或主要结果

本文在 $ \mathbb{T}^2 $ 上建立了四类指数同步定理，分为两个几何机制（共振/无核）与两种观测精度（精确/粗粒化），并统一升级至 $ H^1 $。所有结果均要求参考解 $ u $ 一致有界：存在 $ R > 0 $ 使得 $ \sup_{t\geq0} \big( |u(t)|_{L^2} + \|u(t)\|_{H^1} \big) \leq R $。

1. 精确观测下的共振机制（Theorem 5.2）：设 $ m \in W^{2,\infty}(\mathbb{T}^2; \mathbb{R}^2) $ 为正则共振场（admits a controlled adapted covering，见 Definition 3.2），其剪切缺陷 $ A_m = D_n n - (\nabla \cdot n) n $ 满足粘性可吸收条件 $ \|A_m\|_{L^\infty} < \nu / C_\Phi $，其中 $ C_\Phi $ 为覆盖几何常数。则存在增益阈值 $ \mu_0 = C(R, \nu, C_\Phi, \|m\|_{W^{2,\infty}}) $，使得对任意 $ \mu \geq \mu_0 $，误差满足

2. 精确观测下的无核机制（Theorem 6.4）：设 $ m \in L^\infty(\mathbb{T}^2; \mathbb{R}^2) $ 满足 $ K_m = \{0\} $（即无非零无散隐形电流）。则存在 $ \mu_0 = C(R, \nu, \mathcal{C}_{\text{obs}}) $，其中 $ \mathcal{C}_{\text{obs}} $ 为紧致可观测性常数（由 $ V \hookrightarrow H $ 紧嵌入导出），使得对任意 $ \mu \geq \mu_0 $，有

且无上界增益限制（no upper bound on $ \mu $）。

3. Type-I 粗粒化观测下的共振机制（Theorem 7.3）：设 $ m \in W^{2,\infty} $ 为正则共振场，$ I_h \in \mathcal{L}(L^2(\mathbb{T}^2)) $ 为 Type-I 观测算子（满足 $ L^2 $-稳定性 $ \|I_h\| \leq C $ 及一阶逼近性 $ \|I_h \phi - \phi\|_{L^2} \leq C h \|\nabla \phi\|_{L^2} $）。若分辨率满足增益-分辨率条件 $ \mu h^2 \leq c_0 \nu $（$ c_0 $ 为绝对常数），且剪切缺陷仍满足 $ \|A_m\|_{L^\infty} < \nu / C_\Phi $，则存在 $ \mu_0 $，使得 $ \forall \mu \geq \mu_0 $，

4. Type-I 粗粒化观测下的无核机制（Theorem 7.5）：设 $ m \in W^{1,\infty} $（额外正则性要求），$ I_h $ 同上，且 $ K_m = \{0\} $。若 $ \mu h^2 \leq c_0 \nu $，则存在 $ \mu_0 $，使得 $ \forall \mu \geq \mu_0 $，

5. $ H^1 $-同步升级（Theorem 8.1）：对上述四类 $ L^2 $-同步结果，利用抛物正则化论证（parabolic smoothing），将 $ e $ 视为带 $ L^\infty $-有界外力的 Stokes 方程解，可得

其中 $ \alpha $ 同上，且无需额外观测假设。

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

证明的核心障碍在于：非线性项 $ -b(e,u,e) $ 既不能被 $ \mu|r|^2 $ 单独控制（因 $ r $ 信息不完备），也不能被 $ \nu\|e\|^2 $ 单独控制（因 $ \|e\|^2 $ 与 $ |e|^2 $ 间存在对数型间隙）。作者针对两种几何情形设计了根本不同的分析框架：

1. 共振机制中的移动坐标系与对数估计（Sections 4–5）：当 $ K_m \neq \{0\} $ 时，障碍是隐形电流 $ e = \sigma n $ 的动态。作者引入受控自适应覆盖（controlled adapted covering） $ \Phi : \widetilde{\mathbb{T}}^2 \to \mathbb{T}^2 $，使得 $ \partial_s \Phi = q n \circ \Phi $，从而将 $ n $-流直线化为 $ s $-坐标方向。在此覆盖下，误差 $ e $ 被分解为 $ e = U m + V n $，其中 $ U = r $ 为可观测分量，$ V $ 为隐形分量。关键创新在于移动坐标系下的非线性展开（Lemma 4.3）：

其中 $ b_1 $ 和 $ b_2 $ 均被拉回至 $ \widetilde{\mathbb{T}}^2 $，并利用加权对数三线性估计（Lemma 4.2）控制：

而剪切缺陷项 $ \int e \cdot A_m V $ 被直接由 $ \|A_m\|_{L^\infty} |e|_{L^2} |V|_{L^2} $ 控制。最终，通过 Young 不等式与对数函数的凸性，将 $ b_1, b_2 $ 的上界吸收进 $ \mu|r|^2 + \nu\|e\|^2 $，得到

从而导出指数衰减。此框架的创新在于：将几何缺陷 $ A_m $ 显式分离并量化，使共振理论摆脱对坐标平直性的依赖。

2. 无核机制中的紧致可观测性（Section 6）：当 $ K_m = \{0\} $ 时，障碍是 $ V \hookrightarrow H $ 紧嵌入导出的定性单射性如何转化为定量可观测性。作者的关键洞察是：无需显式构造逆算子或处理小除子，而是利用紧嵌入的抽象性质——对任意 $ \varepsilon > 0 $，存在 $ C_\varepsilon $ 使得

此即“任意小粘性泄漏”的可观测性估计。将其代入误差恒等式，得到

选择 $ \varepsilon = \nu/(2C R^2) $，$ \mu $ 足够大以吸收 $ C_\varepsilon $ 项，即得同步。此方法的创新在于：将可观测性泄漏乘以参考解界 $ R^2 $ 而非增益 $ \mu $，故无上界限制。

3. Type-I 粗粒化分析（Section 7）：对粗观测 $ I_h $，作者建立统一的标量强制性估计：存在 $ c_0 > 0 $，使得

该估计将粗观测问题还原为精确观测框架，代价仅为增益-分辨率条件 $ \mu h^2 \lesssim \nu $。

4. 与之前工作的比较

本文与既有工作的推进体现在三个维度：几何广度、技术深度与机制普适性。

1. 几何广度：从常向量到非恒定场。FLT 定理 [2] 仅处理 $ m $ 为常向量的情形，其证明依赖于 $ \mathbb{T}^2 $ 上 Fourier 模式的正交分解与 $ m $-方向的全局平移对称性。本文通过受控自适应覆盖将这一对称性几何化：对正则共振场 $ m $（如 Example 3.7 中的 $ m_\varepsilon(x) = (1,-\varepsilon \cos x_2)/\sqrt{1+\varepsilon^2 \cos^2 x_2} $），覆盖 $ \Phi_\varepsilon(y,s) = (y+\varepsilon \sin s, s) $ 将 $ n $-流直线化，使 FLT 对数估计得以移植。这显著拓展了适用范围——FLT 无法处理任何非恒定 $ m $，而本文明确给出一类 $ W^{2,\infty} $-正则、剪切缺陷 $ O(|\varepsilon|) $ 的共振场，并证明当 $ |\varepsilon| < \nu/C_\Phi $ 时同步成立。相比之下，[3,4] 中的移动观测器仍基于固定坐标系，未解决投影方向空间变化的根本问题。

2. 技术深度：无小除子的无核理论。对无理方向 $ m $，传统方法需估计 $ n $-流的轨道闭合性，导致可观测性常数 $ \mathcal{C}_{\text{obs}} $ 依赖于 Diophantine 常数 $ \inf_{k \in \mathbb{Z}^2 \setminus \{0\}} |k \cdot n| $，即小除子问题。本文 Section 6 完全规避此困难：利用 $ V \hookrightarrow H $ 的紧性，仅需 $ K_m = \{0\} $ 这一定性条件，即可导出对任意 $ \varepsilon > 0 $ 成立的 $ \|e\|_{H^1}^2 \leq \varepsilon \|e\|_{H^1}^2 + C_\varepsilon |r|^2 $。该论证不显式计算 $ C_\varepsilon $，故无小除子出现。这比 [5] 中对边界观测的谱分析更稳健，也比 [6] 中依赖显式插值常数的网格方法更本质。

3. 机制普适性：四重同步框架的统一与分工。本文首次将 CDA 分解为共振与无核两种互斥且完备的几何范式，并各自配备精确/粗粒化版本。这种分工揭示了根本规律：共振机制的成功依赖于剪切缺陷的粘性可吸收性（$ \|A_m\|_{L^\infty} \ll \nu $），而无核机制的成功依赖于核的空性（$ K_m = \{0\} $）。二者代价不同：共振机制要求 $ m \in W^{2,\infty} $ 且 $ \|A_m\|_{L^\infty} $ 小，但增益下界 $ \mu_0 $ 与 $ \nu $ 无关；无核机制仅需 $ m \in L^\infty $ 且 $ K_m = \{0\} $，但增益下界 $ \mu_0 $ 依赖于 $ \mathcal{C}_{\text{obs}} $（虽未显式计算）。此外，Type-I 粗粒化扩展要求共振情形下 $ m \in W^{2,\infty} $（保障 $ V $ 的正则性），而无核情形下仅需 $ m \in W^{1,\infty} $（保障 $ r $ 的 $ H^1 $-正则性），体现了不同机制对正则性的差异化需求。这种精细分类与代价分析，远超 AOT 算法 [1] 的粗粒化框架或后续文献对单一机制的孤立改进。

8. Giga-Kohn 型结果对全分数阶热方程

Giga-Kohn-type results for the fully fractional heat equation

作者: Yannick Sire, Juncheng Wei, Ke Wu, Zikai Ye

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16785v1

摘要翻译: 我们考虑半线性全分数阶热方程

Comment: 无

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

1. 研究对象和背景

本文研究的模型是全分数阶热方程（fully fractional heat equation）

其中 $\mathbb{R}_- := (-\infty, 0)$，非局部算子 $(\partial_t - \Delta)^\sigma$ 由奇异积分定义：

其中归一化常数 $c_{n,\sigma} = (4\pi)^{-n/2} / |\Gamma(-\sigma)|$。该算子统一了时间与空间的非局部性：当 $u$ 与 $t$ 无关时退化为分数拉普拉斯 $(-\Delta)^\sigma$；当 $u$ 与 $x$ 无关时退化为 Marchaud 分数导数 $\partial_t^\sigma$；且满足 $\lim_{\sigma \to 1^-} (\partial_t - \Delta)^\sigma = \partial_t - \Delta$（强算子收敛意义下）。

该方程具有尺度不变性：若 $u(x,t)$ 是解，则对任意 $\lambda > 0$，缩放函数 $u_\lambda(x,t) = \lambda^{2\beta} u(\lambda x, \lambda^2 t)$ 亦为解，其中 $\beta = \sigma/(p-1)$。由此自然引入自相似变量：

在该坐标下，原方程转化为关于 $v$ 的空间-时间非局部椭圆型方程（见 Lemma 2.4）：

其中 $A = \Gamma(\beta + \sigma)/\Gamma(\beta)$，核 $M_\tau(z,w)$ 显式给出且含非平凡时空耦合结构。特别地，对自相似解（即 $v = v(z)$ 独立于 $s$），方程进一步简化为非局部椭圆方程：

其中对称核 $K(z,w)$ 为 $\tau$-积分后的正则化核（见 Lemma 2.4 表达式），其渐近行为既非齐次亦非卷积型，破坏了经典 Pohozaev 恒等式所需的标度对称性。

研究动机源于经典 Fujita 型抛物方程 $u_t = \Delta u + |u|^{p-1}u$ 的爆破分类理论。Giga–Kohn [10] 在临界与次临界指数范围 $n \leq 2$ 或 $1 < p \leq \frac{n+2}{n-2}$ 下，通过构造自相似变量 $v(z,s)$ 满足 Ornstein–Uhlenbeck 型方程 $Lv + \frac{1}{p-1}v = |v|^{p-1}v$（其中 $L = -\Delta + \frac{1}{2} z \cdot \nabla$），并利用两个关键工具：（i）Pohozaev 恒等式（基于 $L$ 的微分算子结构与高斯测度 $\rho(z) = e^{-|z|^2/4}$ 下的分部积分），以及（ii）能量泛函 $E[v](s) = \frac{1}{2}\int \rho |\nabla v|^2 + \frac{1}{2(p-1)}\int \rho v^2 - \frac{1}{p+1}\int \rho |v|^{p+1}$ 的单调性，证明了所有 Type I 自相似爆破解必为常数，从而唯一性地确定为 $u(x,t) = \pm \kappa (-t)^{-1/(p-1)}$。该框架深刻依赖于局部算子的微分结构、齐次核的标度不变性及 Gauss 测度下 $L$ 的自伴性。

然而，对全分数阶情形，根本障碍在于：$(\partial_t - \Delta)^\sigma$ 是真正空间-时间耦合的非局部算子，既无标准分部积分公式，其核 $M_\tau(z,w)$ 亦丧失齐次性与平移不变性，导致经典 Pohozaev 推导中关键的抵消项消失；同时，Stinga–Torrea 扩展方法虽可将问题降维至局部高维抛物方程，但会引入额外维度并掩盖原方程的本征几何结构。因此，既有技术无法直接迁移。本文旨在填补这一空白：在保持原方程空间-时间非局部本质的前提下，建立首个不依赖扩展的、适用于 $(\partial_t - \Delta)^\sigma$ 的 Giga–Kohn 型单调性公式与 Liouville 型定理，从而系统推进非局部抛物奇性分析理论。

2. 主要定理或主要结果

本文在参数范围 $n \leq 2\sigma$ 或 $1 < p \leq \frac{n + 2\sigma}{n - 2\sigma}$（即分数阶 Sobolev 临界指数）下，建立了两个核心定理：

定理 1.1（自相似解分类）：设 $u$ 是方程 $(\partial_t - \Delta)^\sigma u = |u|^{p-1}u$ 在 $\mathbb{R}^n \times \mathbb{R}_-$ 上的弱解（定义见 Definition 2），且满足先验有界性条件

则 $u$ 必为平凡自相似解：

该结论对所有 $0 < \sigma < 1$ 成立，且 $\kappa$ 中的 Gamma 函数比精确反映了分数阶时间导数对爆破速率的修正。

定理 1.2（一般解分类）：在相同参数范围与先验有界性条件下，若进一步假设

则 $u$ 同样必为上述平凡解：$u(x,t) = \pm \kappa (-t)^{-\beta}$。

两个定理的证明均依赖于作者构建的全新分析工具：（i）基于 Hermite 展开的 Pohozaev 型恒等式族（Theorem 3.1），其核心是将自相似函数 $v(z)$ 展开为 Hermite 多项式 $h_\alpha(z)$，并利用算子 $T_{\beta,\sigma} = \Gamma(L + \beta + \sigma)/\Gamma(L + \beta)$（其中 $L$ 为 Ornstein–Uhlenbeck 算子）的谱表示，得到对每个多重指标 $\alpha$ 成立的恒等式

（ii）空间-时间非局部能量泛函的单调性公式（Proposition 4.2）：定义 $\mathcal{E}[v](s) = \frac{1}{2} \int_{\mathbb{R}^n} \rho(z) \, \mathcal{D}_{\sigma}[v](z,s) \, dz + \frac{1}{2(p-1)} \int_{\mathbb{R}^n} \rho(z) v(z,s)^2 \, dz - \frac{1}{p+1} \int_{\mathbb{R}^n} \rho(z) |v(z,s)|^{p+1} \, dz$，其中 $\mathcal{D}_{\sigma}[v]$ 为由核 $K(z,w)$ 定义的非局部 Dirichlet 型，则存在常数 $C > 0$ 使得

该单调性直接导出 $\mathcal{E}[v](s)$ 非增，结合紧性论证完成定理 1.2 的证明。

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

本文的核心创新在于突破传统非局部分析对扩展方法的依赖，构建了一套基于谱理论与特殊函数的全新技术框架，其关键点分述如下：

1. Hermite 谱表示与非局部算子的对角化：作者观察到，Ornstein–Uhlenbeck 算子 $L = -\Delta + \frac{1}{2} z \cdot \nabla$ 的 Hermite 特征函数 $h_\alpha(z)$ 不仅对 $L$ 正交对角化，而且也构成了全分数阶算子在自相似坐标下的特征基。具体而言，Lemma 2.7 严格证明了自相似方程对应的非局部算子 $T_{\beta,\sigma}$ 具有谱表示

其中 $\Gamma(\cdot)$ 为 Gamma 函数。这一公式将非局部算子完全由局部算子 $L$ 的谱函数刻画，从而将复杂积分核的分析转化为对 $L$ 的离散谱 $\{ \lambda_\alpha = |\alpha| \}_{\alpha \in \mathbb{N}^n}$ 上的函数计算。例如，$T_{\beta,\sigma} h_\alpha = \frac{\Gamma(|\alpha| + \beta + \sigma)}{\Gamma(|\alpha| + \beta)} h_\alpha$。此洞察是全文基石，它使得 Pohozaev 恒等式的构造摆脱了对核函数显式计算的依赖，转而通过 Hermite 展开系数 $\hat{v}_\alpha(s) = \langle v(\cdot, s), h_\alpha \rangle_\rho$ 的演化方程来实现。

2. Hermite 系数层面的 Pohozaev 恒等式族：利用上述谱表示，作者将自相似解 $v(z)$ 的方程 $T_{\beta,\sigma} v = |v|^{p-1} v - A v$（其中 $A = \Gamma(\beta + \sigma)/\Gamma(\beta)$）在每个 Hermite 模式 $h_\alpha$ 上投影，得到一族无穷多个恒等式。关键技巧在于选择特定的测试函数（如 $z_j h_\alpha$）并利用 Hermite 函数的升降算子关系 $\partial_j h_\alpha = \alpha_j h_{\alpha - e_j}$ 和 $z_j h_\alpha = \sum_k c_{\alpha,k} h_{\alpha + e_k}$，将非局部项 $\langle T_{\beta,\sigma} v, z_j h_\alpha \rangle_\rho$ 化为关于 $\hat{v}_\alpha$ 及其邻域系数的代数关系。通过精细调整参数（如选取 $\alpha = 0$ 对应零模），最终证明所有高阶系数 $\hat{v}_\alpha \equiv 0$（$|\alpha| \geq 1$），从而 $v$ 必为常数。该方法完全规避了经典 Pohozaev 中因核非齐次导致的不可控余项，且其灵活性允许推广至更一般的“主方程”（master equations）类算子。

3. 空间-时间非局部单调性公式的构造：为证明定理 1.2，作者定义了一个内蕴的非局部能量 $\mathcal{E}[v](s)$，其非局部 Dirichlet 型部分为

其中 $K(z,w)$ 由 Lemma 2.4 给出。证明单调性的核心在于计算 $\frac{d}{ds} \mathcal{E}[v](s)$ 时，非局部项的时间导数 $\partial_s v$ 与空间差分 $v(z) - v(w)$ 形成完美配对，产生负定耗散项 $$

• \iint \rho(z) K(z,w) , |v_s(z,s) - v_s(w,s)|^2 , dw, dz. $$ 这一结构源于 $K(z,w)$ 的对称性及 $\rho$ 测度下 $T_{\beta,\sigma}$ 的自伴性，是首个不借助 Stinga–Torrea 扩展、直接在原时空上建立的非局部单调性公式。

4. 与之前工作的比较

本文与既有工作存在三重根本性区别，分别体现在技术路径、适用范围与理论深度上：

1. 与经典 Giga–Kohn [10] 的比较：Giga–Kohn 的工作严格依赖 $\sigma = 1$ 时的局部结构——其 Pohozaev 恒等式 (1.5) 本质是 $L$ 在 Gauss 测度下的分部积分恒等式，而单调性公式则源于局部能量泛函的梯度流结构。本文首次将该框架推广至 $0 < \sigma < 1$，但并非简单替换算子，而是彻底重构证明逻辑：用 Hermite 谱理论替代微分算子的分部积分，用 Gamma 函数谱表示替代拉普拉斯算子的齐次性，用非局部 Dirichlet 型替代经典 Dirichlet 能量。值得注意的是，作者强调该新证明对 $\sigma = 1$ 依然成立，从而为经典结果提供了全新的、更具普适性的谱解释，这是对原有理论的深化而非平移。

2. 与 Stinga–Torrea 扩展方法的比较：Bustamante [2] 等人通过将 $u$ 扩展至 $\mathbb{R}^{n+1}_+ \times \mathbb{R}$ 并应用高维抛物-椭圆变换，获得了单调性公式。然而，此类方法引入额外维度，导致能量泛函依赖于扩展变量，且无法直接提取原方程的时空耦合信息。本文 Proposition 4.2 的单调性公式完全在原时空 $\mathbb{R}^n \times \mathbb{R}$ 上定义，其耗散项 $\iint \rho K |v_s(z) - v_s(w)|^2$ 清晰揭示了分数阶时间导数如何通过非局部核 $K$ 耦合空间差异，这是扩展方法无法捕捉的本征结构。文中明确指出：“据我们所知，这是首个不使用 Stinga–Torrea 扩展的空间-时间非局部单调性公式”。

3. 与 Deng–Sire–Wei–Wu [7] 的半拉普拉斯情形比较：[7] 研究方程 $\partial_t u + (-\Delta)^{1/2} u = |u|^{p-1} u$，其优势在于 $(-\Delta)^{1/2}$ 具有显式 Poisson 核，允许直接计算。但该方法受限于具体核的代数性质，难以推广至一般 $\sigma$ 或时空耦合情形。本文处理的 $(\partial_t - \Delta)^\sigma$ 无显式核，且其时空混合特性使标准分数阶技术失效。作者通过 Hermite 谱方法绕过核的显式计算，实现了对任意 $\sigma \in (0,1)$ 的统一处理，显著拓宽了适用范围。代价在于：当前结果要求 $p$ 在分数阶 Sobolev 临界值 $\frac{n+2\sigma}{n-2\sigma}$ 以下，而 [7] 在 $n=1$ 时得到稍优的 $p < \frac{n+1}{n-1}$；此外，本文未处理超临界情形或 Type II 爆破，这些仍是开放问题。从提供内容中无法确认本文是否改进了 [8] 关于 Fujita 指数 $p_{F,\sigma} = 1 + \frac{2\sigma}{n+2-2\sigma}$ 的存在性结果，作者亦未显式展开比较。

9. 非变分型自由边界系统存在性理论

Existence theory for non-variational systems with free boundaries

作者: Layan El Hajj, Seongmin Jeon, Henrik Shahgholian

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16767v1

摘要翻译: 我们研究定义在一般光滑区域上的椭圆型与抛物型偏微分方程组解的存在性，其右端项可能奇异，且具有自由边界。

Comment: 43 pages

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

1. 研究对象和背景

本文研究如下两类非变分耦合系统解的存在性：一类是椭圆型系统，另一类是抛物型系统，二者均具有潜在的奇异右端项与自由边界。核心研究对象为定义在光滑有界区域 $\Omega \subset \mathbb{R}^n$ 上的两分量弱耦合系统，其椭圆情形精确表述为：

其中 $L_1, L_2$ 是一般非变分算子（如完全非线性一致椭圆算子或 $p$-拉普拉斯算子 $\Delta_p u = \operatorname{div}(|\nabla u|^{p-2} \nabla u)$）；$f_i \in L^\infty(\Omega; \mathbb{R}_+)$，$g_i \in C(\partial\Omega; \mathbb{R}_+)$ 满足 $0 < c_1 \leq g_i \leq C_1$；非线性项 $h_i$ 属于函数类 $\mathcal{H}$，即满足：对某 $-1 < a < 1$、$C_0 > 0$ 及上界 $M_0 = \sup_{\partial\Omega} |g_i|$，有

耦合函数 $k_i$ 属于类 $\mathcal{K}$：$k_i(t) = 0$ for $t \leq 0$，且在 $\mathbb{R}$ 上非减连续。该结构天然诱导自由边界 $\partial\{u > 0\}$ 与 $\partial\{v > 0\}$，其几何性质（如是否重合或分离）亦为研究重点。

此类系统源于变分框架中能量泛函的极小化，例如标量泛函 $\int_\Omega (|\nabla u|^2 + |u|^\gamma)\,dx$（$0 < \gamma < 2$），其欧拉–拉格朗日方程为 $-\Delta u = \gamma |u|^{\gamma-2} u$，右端在 $u=0$ 处奇异（当 $\gamma < 1$ 时不可积；当 $1 \leq \gamma < 2$ 时为吸收项），导致死核 $\{u = 0\}$ 的形成，其边界即自由边界。然而，当推广至非变分耦合系统（如原型 $\Delta u = v^a$, $\Delta v = u^b$），问题本质改变：右端项依赖于另一未知量，且在自由边界附近可能同步退化或奇异性增强，致使标准迭代法、不动点定理或能量方法失效——因缺乏统一泛函结构，无法控制迭代中范数的累积爆炸；又因奇异源项恰位于待确定的自由边界邻域，经典先验估计难以建立。现有理论（如 [17]）多局限于球对称情形或特殊结构（如单调性、可积性假设），缺乏对一般光滑域、非变分算子及混合算子的统一处理框架。本文旨在填补这一空白：构建一个不依赖具体算子形式、仅需验证三项结构性质（可解性、比较原理、一致极限稳定性）的抽象存在性框架，并将其切实应用于完全非线性算子与 $p$-拉普拉斯算子两大重要族，尤其攻克 $p$-抛物情形下具奇异/退化右端项解的正则性这一文献缺失的关键环节。

2. 主要定理或主要结果

本文核心成果由两个元定理（meta-theorems）及其实例化构成。首先，椭圆元定理（Theorem 1）断言：设 $\Omega \subset \mathbb{R}^n$ 为有界光滑域，$g_i \in C(\partial\Omega; \mathbb{R}_+)$ 满足 $c_1 < g_i < C_1$，$f_i \in L^\infty(\Omega; \mathbb{R}_+)$，$h_i \in \mathcal{H}$，$k_i \in \mathcal{K}$，且算子 $L_i$ 在函数空间 $F_i(\Omega)$ 上满足三条结构性质（PE1）–（PE3）：（PE1）可解性：对任意 $\tilde{h}_i \in \widetilde{\mathcal{H}}$（正则化类）、$f \in L^\infty$、$g \in C(\partial\Omega)$，Dirichlet 问题 $L w = f \tilde{h}(w)$ 在 $\Omega$ 内存在解 $w$，且具 Hölder 连续性 $|w(x_1) - w(x_2)| \leq \rho(|x_1 - x_2|)$ 及单侧增长估计 $\sup_{B_r(x) \cap \Omega'} w \leq C(r^{1+\alpha} + w(x))$；（PE2）比较原理：若 $f_1 < f_2$，则对应解 $w_1 \geq w_2$；（PE3）收敛性：若 $f_j \to f$、$w_j \to w$ 局部一致、$\tilde{h}_j \to \hat{h}$ 在 $(0,\infty)$ 上局部一致，则 $L w = f \hat{h}(w)$ 在 $\{w > 0\}$ 内成立。在此条件下，系统 (1.1) 存在解 $(u,v) \in F_1(\Omega) \times F_2(\Omega)$。其次，抛物元定理（Theorem 7）给出平行结论，适用于时空柱 $\Omega_T = \Omega \times (-T,0]$，要求算子 $L_i$ 满足对应抛物性质（PP1）–（PP3），并引入抛物距离 $d(X,Y) = \max\{|x-y|, |t-s|^{1/2}\}$ 与抛物 Hölder 正则性。

进一步，作者将元定理实例化：对完全非线性一致椭圆/抛物算子（如 Pucci 算子 $\mathcal{M}^\pm_{\lambda,\Lambda}(D^2 u)$），验证其满足（PE1）–（PE3）及（PP1）–（PP3）；对 $p$-拉普拉斯算子 $\Delta_p u$，证明其在 $W^{1,p}(\Omega)$ 中满足（PE1）–（PE3）；对 $p$-抛物算子 $\partial_t u - \Delta_p u$，通过附录 A 严格建立其解的正则性：在奇异情形（右端 $\in L^m$，$m > 1$）下，解 $u \in L^\infty_{\rm loc}(\Omega_T)$ 且具局部 Hölder 连续性；在退化情形（右端 $\in L^1$）下，解具 $L^p$ 时间导数与梯度的可积性。此外，附加非退化性（PE2’）与 Hopf 边界引理（PE1’）后，Theorem 2 给出自由边界几何刻画：若 $f_i > 0$、$k_i > 0$ on $(0,\infty)$、$h_i(t) \geq c_0 t^b$ ($0 < b < 1$)，则零集 $\{u=0\}^\circ$ 与 $\{v=0\}^\circ$ 的任意连通分支 $C_u$、$C_v$ 必满足 $C_u = C_v$ 或 $C_u \cap C_v = \emptyset$；且存在参数 $(f_i,g_i)$ 使二者非空且互斥。

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

本文证明的核心障碍在于：如何在缺乏统一变分结构、右端项于自由边界处奇异/退化、且耦合机制导致相互依赖的三重困难下，构造一致有界的解序列并保证极限满足原始方程与自由边界条件。作者的真正创新点在于提出并实施一套“元定理驱动”的抽象框架，其技术内核并非直接估计，而是将存在性问题转化为对算子结构性质的验证问题。具体而言，证明框架包含三个关键层级：

第一层级：单调迭代方案的设计与控制。针对正则化系统 $L_1 u_\varepsilon = (f_1 k_1(v_\varepsilon) + \varepsilon) \tilde{h}_1(u_\varepsilon)$，作者构造双重单调序列 $\{u_j\}$（递增）与 $\{v_j\}$（递减），其构造依赖于（PE2）比较原理的严格不等式版本：利用 $\sigma_j \searrow 1$、$\tau_j \nearrow 1$ 的扰动序列，确保 $f_1 k_1(v_{j-1}) + \varepsilon \sigma_j > f_1 k_1(v_j) + \varepsilon \sigma_{j+1}$，从而触发 $u_j \leq u_{j+1}$；同理 $v_j \geq v_{j+1}$。此设计巧妙规避了直接处理耦合非线性的困难，将系统解的存在性归约为标量问题的可解性（PE1）与比较性（PE2）。

第二层级：一致先验估计的传递与紧性。序列 $\{u_j\}$、$\{v_j\}$ 被上下界 $u^*, v^*, u_*, v_*$ 控制（由（PE1）中的屏障函数 $w^*$ 与调和函数 $w^*$ 提供），结合（PE1）给出的统一 Hölder 模 $\rho$ 与增长估计 $\sup_{B_r(x)} w \leq C(r^{1+\alpha} + w(x))$，应用 Arzelà–Ascoli 定理获得局部一致收敛子列 $u_j \to u_\varepsilon$、$v_j \to v_\varepsilon$。此处关键在于（PE1）的增长估计（2.6）不仅保障 Hölder 连续性，更在自由边界点 $x \in \partial\{u > 0\}$ 处强制 $|\nabla u(x)| = 0$，因为若 $|\nabla u(x)| > 0$，则由隐函数定理 $u$ 在 $x$ 邻域内线性非零，与 $u(x)=0$ 矛盾；而估计本身蕴含 $u(y) \leq C |y-x|^{1+\alpha}$，故 $|\nabla u(x)| = 0$。

第三层级：极限过程与自由边界条件的恢复。利用（PE3）收敛性，由 $u_\varepsilon \to u$、$v_\varepsilon \to v$ 局部一致及 $f_i k_i(v_\varepsilon) + \varepsilon \to f_i k_i(v)$ 局部一致，推出 $L_i u = f_i k_i(v) h_i(u)$ in $\{u > 0\}$。自由边界梯度条件 $|\nabla u| = 0$ on $\partial\{u > 0\}$ 直接由（PE1）的增长估计（2.6）在极限下继承；而 Dirichlet 边界条件 $u = g_1$ on $\partial\Omega$ 则通过上下屏障 $u^* \leq u_\varepsilon \leq u^*$ 且 $u^* = u^* = g_1$ on $\partial\Omega$ 得到。对于 $p$-抛物情形，附录 A 的核心技术是建立解的局部有界性与 Hölder 连续性：在奇异情形，采用 De Giorgi 迭代结合 $L^m$ 右端的嵌入估计；在退化情形，运用 $p$-抛物 Sobolev 不等式与 intrinsic scaling 方法克服退化性，证明解属于 $L^p_{\rm loc}(\Omega_T)$ 且梯度具适当可积性，为（PP1）–（PP3）的验证奠定基础。

4. 与之前工作的比较

本文与既有工作存在多维度实质性推进。首先，相较于 Giaquinta–Giusti [20,21] 的变分框架，本文彻底摆脱能量泛函依赖，将存在性理论拓展至非变分系统（如完全非线性算子），其元定理框架可覆盖任何满足（PE1）–（PE3）的算子，远超经典变分方法的适用范围。其次，对比 [17] 的球对称设定，本文在一般光滑域 $\Omega$ 上建立理论，且通过引入 Lipschitz 边界数据与距离函数屏障，在完全非线性情形下绕过粘性解理论中构造屏障的困难；而对 $p$-拉普拉斯情形，利用变分结构仅需 $C^{0,1}$ 边界，显著弱于 [17] 所需的更高正则性。第三，在技术层面，本文首次系统处理 $p$-抛物算子与奇异/退化右端项的组合，附录 A 中建立的正则性结果（特别是退化情形下的 $L^p$ 梯度估计）是文献中缺失的关键环节，此前工作（如 DiBenedetto）多假设右端有界或 $L^q$（$q > 1$），而本文将 $q=1$ 的临界情形纳入框架。第四，关于自由边界几何，[17] 假设右端在 $v \to 0$ 时爆破（$\sim v^{-\beta}$），强制 $\{u=0\} = \{v=0\}$；本文允许 $k_i(v)$ 在 $v=0$ 处有界（仅非减连续），故 Theorem 2 揭示了更丰富的几何可能性：零集可分离，这反映了耦合机制的内在灵活性。代价在于：元定理的抽象性要求对每个具体算子严格验证（PE1）–（PE3），如 $p$-抛物情形的验证耗尽附录 A 的全部篇幅；且比较原理（PE2）需严格不等式 $f_1 < f_2$，虽在构造中可行，但弱于通常的 $\leq$ 形式，表明框架在某些极端情形下可能存在技术局限。最后，本文未涉及解的唯一性或更高阶正则性（如自由边界的 $C^{1,\alpha}$ 性），这些仍是开放问题，而 [5,6] 等文献虽讨论死核生成机制，但缺乏统一的存在性构造工具。

10. 能量超临界热方程的坍缩管型 II 类爆破

Collapsing-tube type II blow-up for the energy-supercritical heat equation

作者: Manuel del Pino, Monica Musso, Juncheng Wei, Yifu Zhou

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16733v1

摘要翻译: 我们为能量超临界热方程

Comment: 无

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

1. 研究对象和背景

本文研究能量超临界半线性热方程

在正解、单点爆破（blow-up at the origin only）、有限时间 $ t \nearrow T < \infty $ 情形下的奇性形成机制。该方程的能量泛函为

其临界性由 Sobolev 指数 $ p_S = \frac{n+2}{n-2} $ 刻画：当 $ p = 3 > p_S $ 对所有 $ n \geq 5 $ 成立，故属能量超临界（energy-supercritical）范畴。值得注意的是，$ p = 3 $ 落在 Matano–Merle 区间

此区间内，经典结果 [27,28,31] 已严格证明：在 $ \mathbb{R}^n $ 或球形有界域中，任何径向正解均不可能发生有限时间 II 类爆破——即排除了所有已知对称性下最自然的爆破路径。因此，构造一个正的、单点、有限时间 II 类爆破解，本身即构成一个根本性存在性障碍。

为突破该障碍，作者引入强各向异性对称性：将变量分组为 $ x = (x', z) \in \mathbb{R}^{n-3} \times \mathbb{R}^3 $，定义柱坐标 $ r = |x'| $，并限定解为轴对称且关于 $ z $ 偶函数，即 $ u(x,t) = \tilde{u}(r,z,t) $。在此对称性下，拉普拉斯算子化为

故原方程等价于四维空间 $ (r,z) \in \mathbb{R}^+ \times \mathbb{R}^3 $ 上带几何漂移项的方程：

关键观察在于：非线性指数 $ p = 3 $ 恰好是 $ \mathbb{R}^4 $ 的 Sobolev 临界指数 $ p_S^{(4)} = \frac{4+2}{4-2} = 3 $。因此，问题被“降维”为一个四维临界热方程，其内在奇性由 Aubin–Talenti 泡泡 $ U(y) = \alpha_0 (1 + |y|^2)^{-1} $，$ y \in \mathbb{R}^4 $，所主导，但叠加了一个由高维背景 $ n $ 引入的、奇异的几何漂移项 $ \frac{n-4}{r} \partial_r \tilde{u} $。该漂移项在 $ r \to 0 $ 处发散，其物理含义是 $ (n-4) $-维球面 $ \{ |x'| = r,\, z = 0 \} $ 的平均曲率贡献；它不改变方程的临界性，却彻底改变了爆破几何——不再是点集中，而是沿一个维度为 $ n-4 $ 的球面集中，且该球面自身以自相似速率坍缩至原点。这一设定将问题置于一个此前未被系统研究的交叉地带：一方面，它规避了 Matano–Merle 禁止定理所依赖的径向对称性假设；另一方面，它又非简单地恢复到低维临界情形，而是将临界泡泡嵌入一个动态坍缩的几何背景中，从而催生全新的调制动力学。本文试图填补的核心空白是：在能量超临界、正解、单点爆破的严格框架下，首次定量刻画一种由几何坍缩驱动的、具有明确尺度分离的 II 类爆破机制，并揭示其与四维临界现象及高维轴对称热核之间的深层耦合。

2. 主要定理或主要结果

在 $ n \geq 5 $ 下，存在初始数据 $ u_0 \in C^\infty(\mathbb{R}^n) $（或在光滑有界域 $ \Omega \subset \mathbb{R}^n $ 上满足相应对称性的边界/初值），使得方程 $ u_t = \Delta u + u^3 $ 的解 $ u(x,t) $ 在有限时间 $ T > 0 $ 发生单点爆破，且爆破仅发生在原点。更精确地，存在参数函数 $ \lambda(t) > 0 $、$ \xi_r(t) > 0 $，满足以下渐近律：

1. 主爆破轮廓：在柱坐标 $ r = |x'| $, $ z \in \mathbb{R}^3 $ 下，解渐近于一个被高度拉伸的四维 Aubin–Talenti 泡泡：

   $$ u(x,t) \sim \frac{1}{\lambda(t)} U\left( \frac{r - \xi_r(t)}{\lambda(t)}, \frac{z}{\lambda(t)} \right), \quad \text{其中 } U(y) = \alpha_0 (1 + |y|^2)^{-1},\; y \in \mathbb{R}^4. $$

2. 几何坍缩尺度：浓度球面 $ \{ |x'| = \xi_r(t),\, z = 0 \} $ 的半径以标准抛物自相似速率坍缩：

   $$ \xi_r(t) = \sqrt{2(n-4)(T-t)} + o\big(\sqrt{T-t}\big). $$

3. 横向爆破尺度：泡泡在垂直于球面的方向（即 $ (r,z) $ 平面内）的厚度遵循一个带有精确对数修正的 II 类速率：

   $$ \lambda(t) = \kappa_* \frac{T-t}{|\log(T-t)|^{\frac{n}{n-2}}} + o\left( \frac{T-t}{|\log(T-t)|^{\frac{n}{n-2}}} \right), \quad \text{其中 } \kappa_* > 0 \text{ 为一普适常数}. $$

4. 尺度分离：两个尺度满足严格的分离关系

   $$ \lambda(t) \ll \sqrt{T-t}, $$

   即横向厚度远小于球面半径，从而形成一个“薄管”（thin tube）状奇性结构。该结果被精确表述为定理 1.1，其证明要求 $ T > 0 $ 充分小（即构造适用于小时间尺度的爆破解）。此外，文中指出，当 $ n = 4 $ 时，坍缩球面退化为一点（$ n-4 = 0 $），漂移项消失，此时速率退化为 $ \lambda(t) \sim (T-t)/|\log(T-t)|^2 $，与 $ \mathbb{R}^4 $ 上临界立方热方程的稳定 $ k=1 $ 率 [16] 一致，验证了该构造的连续性。

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

证明的核心是精细的内–外匹配（inner–outer gluing）方案，其创新点远超标准能量估计或紧性方法，而在于对多重尺度、非局部耦合与几何奇性的协同处理。

首要障碍与核心创新：非局部调制方程的导出与求解。标准 II 类爆破中，尺度参数 $ \lambda(t) $ 通常满足一个局部 ODE（如 $ \dot{\lambda} \sim -C \lambda^2/(T-t) $）。本文的关键突破在于，由于四维临界泡泡 $ U $ 的长程尾部 $ U(y) \sim |y|^{-2} $ 与高维轴对称热核 $ \Gamma_n $ 的非平凡卷积，导致 $ \lambda(t) $ 的演化被一个非局部积分微分方程所支配：

该方程的来源是：将误差投影到线性化算子 $ L_0 = \Delta + 3U^2 $ 的零空间（scaling mode）上，其右端项包含一个由全局校正 $ \Psi_0 $ 引入的记忆项，该记忆项正是上述积分形式。求解此方程需进行高精度的渐近分析，最终导出对数幂律 $ \lambda(t) \sim (T-t)/|\log(T-t)|^{n/(n-2)} $。这是本文最本质的技术创新，它将爆破律的起源从局部ODE提升到了非局部动力学层面。

第二重障碍与创新：“再匹配”（re-gluing）技术。由于爆破发生在自相似区域 $ r \sim \sqrt{T-t} \to 0 $，传统的内问题（在缩放变量 $ y = (r-\xi_r, z)/\lambda $ 中）会因原始漂移项 $ \frac{n-4}{r} \partial_r $ 在 $ r \to 0 $ 处的奇性而失效。为克服此困难，作者发展了一种“再匹配”技术：在内问题中，将解进一步分解为 scaling、translation 和高阶模态，并对最棘手的 scaling 和 translation 模态，再次应用内–外匹配，引入一个“再缩放”变量和新的截断函数。这本质上是在内区内部构建了一个更精细的层次结构，以隔离并控制由 $ r \to 0 $ 奇性引发的破坏性项。

第三重障碍与创新：轴对称热核的显式构造与应用。标准泡泡 $ U_{\lambda,\xi} $ 的误差中，存在一个慢衰减项 $ E_0 \approx -\alpha_0 \dot{\lambda}/\rho^2 $，无法作为小扰动处理。为此，作者在附录 B.1 中显式构造了 $ \mathbb{R}^{n-3} \times \mathbb{R}^3 $ 上的轴对称热核 $ \Gamma_n(r,z; \tilde{r},\tilde{z}; \tau) $，并定义全局校正 $ \Psi_0 $ 为其与该慢误差的卷积：

这一构造是连接四维临界现象与 $ n $ 维背景的桥梁，其作用在第 4 节被量化：它正是非局部调制方程中关键记忆项的物理来源。

第四重障碍与创新：加权空间的精巧选取。为使不动点论证闭合，内区半径 $ R(t) $ 必须精心选择：

此选择确保了物理尺度 $ \lambda(t) R(t) \ll \sqrt{T-t} $，即内区严格位于坍缩球面的自相似区域内，同时其对数幂次 $ \theta $ 又足够大，以保证所需的时间可积性。这一细节是整个分析得以成立的基石。

4. 与之前工作的比较

本文与既有工作存在多维度的本质区别，主要体现在爆破几何、适用范围、技术深度与理论启示上。

与径向爆破工作的比较：Herrero–Velázquez [22,23] 构造了 Joseph–Lundgren 指数 $ p > p_{JL}(n) $（$ n \geq 11 $）下的径向 II 类爆破，其速率是幂律 $ \lambda(t) \sim (T-t)^{\beta} $。而本文工作处于互补的 Matano–Merle 区间 $ p_{S} < p < p_{JL} $，且完全绕开了径向对称性。正如引言 1.3 节所强调，Matano–Merle 定理 [27,28,31] 已证明在此区间内径向爆破不可能存在。本文通过引入 $ (n-3,3) $ 分组的轴对称性，将问题转化为一个带漂移的四维临界问题，从而在同一个 $ p=3 $、$ n \geq 5 $ 的参数下，实现了对禁令的突破，提供了首个该参数下的正解单点 II 类爆破实例。

与固定几何集爆破的比较：[14] 构造了在能量超临界热方程中沿一个固定圆（$ n=5 $ 时）的 II 类爆破。其几何集不随时间变化，浓度尺度为 $ \mu(t) \sim (T-t)/|\log(T-t)|^2 $，与 $ \mathbb{R}^4 $ 临界情形一致。本文则首次实现了几何集自身的坍缩，其半径 $ \xi_r(t) \sim \sqrt{T-t} $，且横向尺度 $ \lambda(t) $ 的对数幂次 $ n/(n-2) $ 显式依赖于环境维数 $ n $，这是固定几何集情形所不具备的全新自由度。这种“坍缩管”机制是几何驱动的，而非单纯由非线性强度驱动。

与各向异性爆破的比较：Collot–Merle–Raphaël [9] 发现了在孤立点处不同方向具有不同爆破速率的机制，但其各向异性源于线性化算子的各向异性谱，而非一个宏观坍缩的几何结构。本文的各向异性是几何本征的：径向 $ r $ 方向由漂移项主导，以 $ \sqrt{T-t} $ 速率坍缩；而横向 $ (r,z) $ 平面内则由四维临界泡泡主导，以更慢的 $ \lambda(t) $ 速率集中。二者通过同一个调制方程耦合，形成统一的两尺度图像。

与数值场景的比较：Hou–Huang [24] 在三维轴对称 Euler 方程的数值模拟中观察到一个“坍缩环”（collapsing ring）的潜在奇性，其环心位置 $ \xi(t) \sim \sqrt{T-t} $，厚度 $ \lambda(t) \sim T-t $。本文在标量抛物模型中，为这一数值现象提供了首个严格解析实现，不仅确认了 $ \xi_r(t) \sim \sqrt{T-t} $ 的鲁棒性（由显式调制方程导出），更揭示了其厚度尺度受对数修正的深层原因（四维临界尾部与热核的非局部相互作用）。这表明，坍缩环几何可能是轴对称流体奇性的一类普适模式。

代价与局限：该构造的代价是极高的技术复杂性——需要非局部方程、再匹配、显式热核等全套工具。其局限在于目前仅适用于 $ p=3 $ 这一特定指数（因其恰好是 $ \mathbb{R}^4 $ 的临界指数），且对称性要求严格（$ (n-3,3) $ 分组）。文中第 1.3 节的猜想探讨了推广到其他临界维度 $ k $ 的可能性，但尚未实现。从提供内容中无法确认该解是否稳定，也未显式展开与 [12,15,35] 等临界热方程工作的稳定性比较。

11. 动力学 Lamb-Oseen 涡旋的初始层不稳定性

Initial layer instability of the kinetic Lamb-Oseen Vortex

作者: Michele Dolce, Isabelle Gallagher

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16729v1

摘要翻译: 众所周知，当马赫数与克努森数趋于零时，不可压 Navier-Stokes 方程组的解在各种意义下（弱或强）是 Boltzmann 方程解的极限。特别是对于光滑解的情形，目前已有较深入的理解。近期工作致力于选取尽可能接近不可压 Navier-Stokes 方程组适定性所对应函数空间的初值。本文研究二维情形下测度值初值涡量：在初值涡量为 Dirac 质点的特殊情形下，Navier-Stokes 方程组存在唯一解，即热方程的解。我们证明，该初值涡量的动力学激发（经轻微磨光）所诱导的 Boltzmann 方程解，在一个极小的时间层内以强方式偏离预期的流体力学极限。

Comment: 无

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

1. 研究对象和背景

本文研究二维空间中 Boltzmann 方程在 Knudsen 数 $\varepsilon \to 0$ 极限下的流体动力学极限行为，聚焦于一类具有测度初值（特别是 Dirac 质点涡旋）的奇异初始数据所激发的 kinetic 解与对应 Navier-Stokes 流体解之间的偏差。核心对象是如下缩放后的 Boltzmann 方程：

其中 $L$ 是线性化玻尔兹曼算子，$\mathrm{Ker}(L) = \mathrm{span}\left\{ \sqrt{\mu},\, v_1 \sqrt{\mu},\, v_2 \sqrt{\mu},\, \left(\frac{1}{2}|v|^2 - 1\right) \sqrt{\mu} \right\}$，投影 $P$ 将 $f$ 分解为宏观部分 $P f = (\rho[f] + u[f] \cdot v + \theta[f] (\frac{1}{2}|v|^2 - 1)) \sqrt{\mu}$ 和微观部分 $P^\perp f$；宏观变量定义为

目标流体系统为二维不可压 Navier-Stokes-Fourier 系统：

其经典解由 Lamb–Oseen 涡旋给出：初始涡度 $\omega^{\mathrm{in}} = \delta_{x=0}$ 对应速度场 $u^{\mathrm{in}}(x) = \frac{1}{2\pi} \frac{x^\perp}{|x|^2}$，满足 $u^{\mathrm{LO}}(t,x) = e^{\nu t \Delta} u^{\mathrm{in}}(x)$，且 $\omega^{\mathrm{LO}}(t,x) = \frac{1}{\nu t} G\left(\frac{x}{\sqrt{\nu t}}\right)$，其中 $G(\xi) = \frac{1}{4\pi} e^{-|\xi|^2/4}$。该解是 $-1$-齐次初值驱动的自相似演化，且是二维 Navier-Stokes 全局吸引子。问题在于：能否从 Boltzmann 方程出发，以 kinetic Lamb–Oseen 初值 $f^{\varepsilon}_{\mathrm{in}}(x,v) = u^{\mathrm{in}}(x) \cdot v \sqrt{\mu}(v)$（或其光滑正则化）出发，在 $\varepsilon \to 0$ 极限下恢复该流体解？现有理论（如 [12]）要求初值属于 $X^{1+,k}$ 类（即 $H^{1+\eta}_x$ 关于 $v$ 加权 $L^\infty_k$），但 $u^{\mathrm{in}} \notin L^2_{\mathrm{loc}}$，更不满足 $H^1$ 正则性；其 $|x|^{-1}$ 奇异性使 $f^{\varepsilon}_{\mathrm{in}}$ 不属于任何已知 Boltzmann Cauchy 问题适定性框架。因此，本文旨在填补这一关键空白：当 Navier-Stokes 解允许测度初值（如 Dirac 涡旋）时，Boltzmann 动力学是否仍能实现该极限？若不能，则其失效机制、时间尺度与定量偏差为何？这触及 kinetic 理论中‘奇异性传播’与‘多尺度松弛’的根本问题——Navier-Stokes 中的瞬时压力平衡与不可压约束是否能在 kinetic 演化中被即时建立？

2. 主要定理或主要结果

在 $k > 2$ 下，存在普适常数 $C_\star > 0$、大常数 $K \gg 1$、小常数 $\delta > 0$ 及 $0 < a_0 < b_0 < \infty$、$C,c > 0$（均独立于 $\varepsilon$），使得对充分小的 $\varepsilon$，令 $T_\varepsilon := \frac{\delta}{\varepsilon^2 C_\star |\ln \varepsilon|^6}$，则对所有 $t \leq T_\varepsilon$，有： （i）$X^{1,k}$-范数偏差：$\| (f^\varepsilon - f^\varepsilon_{\mathrm{LO}})(t) \|_{X^{1,k}} \geq C \frac{t}{\varepsilon^3}$； （ii）宏观速度分量偏差：在环形区域 $A_\varepsilon := \{ x \in \mathbb{R}^2 : a_0 \varepsilon \leq |x| \leq b_0 \varepsilon \}$ 内，

（iii）不可压性破坏：$| \nabla_x \cdot u[f^\varepsilon](t,x) | \geq c \frac{t^2}{\varepsilon^6}$； （iv）密度与温度扰动：$| \rho[f^\varepsilon](t,x) | + | \theta[f^\varepsilon](t,x) | \geq c \frac{t^3}{\varepsilon^7}$。 其中 $f^\varepsilon_{\mathrm{LO}}(t) := (e^{\nu t \Delta_x} u^\varepsilon_{\mathrm{in}}) \cdot v \sqrt{\mu}$ 是期望的流体目标解，$u^\varepsilon_{\mathrm{in}}(x) = \frac{1}{2\pi} \chi(\varepsilon^\gamma |x|) \frac{x^\perp}{|x|_\varepsilon^2}$ 为带 Gevrey 截断的正则化初值，$|x|_\varepsilon := \sqrt{|x|^2 + K^2 \varepsilon^2}$，$0 < \gamma < 1$。所有下界在 $t \ll \varepsilon^2$ 的初始层内成立，且最优时间尺度 $T_\varepsilon$ 含对数修正因子 $|\ln \varepsilon|^{-6}$，源于递归展开阶数 $M \sim |\ln \varepsilon|$ 的误差控制需求。

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

证明的核心障碍在于：$-1$-齐次初值导致空间奇异性与 $\varepsilon$-缩放耦合，使标准 Chapman-Enskog 展开失效，且无法直接应用 $L^2$ 或 Sobolev 能量估计。作者的真正创新点在于构建一个精细的、逐阶追踪奇异性演化的递归近似解框架，并严格量化每阶的 $\varepsilon$-依赖性。具体分为三步： （1）奇异性分离与类 $P_\varepsilon(n)$ 函数类：定义函数类 $P_\varepsilon(n)$，其元素形如 $g(x,v) = \frac{1}{\varepsilon^n |x|^{n+1}_\varepsilon} \sum_{j=0}^n \left( \frac{|x|_\varepsilon}{\varepsilon} \right)^j \Phi_{n,j}(r(x),v,z(x))$，其中 $r(x) = x / |x|_\varepsilon$，$z(x) = \varepsilon^\gamma x$，$\Phi_{n,j}$ 具备高阶 Gevrey 正则性与 Gaussian 衰减 $\mu_n(v) = e^{-\beta_n |v|^2/4}$（$\beta_n = \frac{1}{4} - \frac{n}{32M}$）。该结构将 $|x|^{-1}$ 奇异性完全吸收到 $|x|^{-n-1}_\varepsilon$ 因子中，而 $\Phi_{n,j}$ 保持光滑有界，从而将奇异性分析转化为对 $\varepsilon$ 幂次的精确计数。 （2）递归构造与误差控制：基于 (1.3) 导出递归关系 $n f^\varepsilon_n = \varepsilon^{-2} L f^\varepsilon_{n-1} + \varepsilon^{-1} \left( -v \cdot \nabla_x f^\varepsilon_{n-1} + \sum_{k=0}^{n-1} \Gamma(f^\varepsilon_k, f^\varepsilon_{n-1-k}) \right)$，并证明若 $f^\varepsilon_{n-1} \in P_\varepsilon(n-1)$，则 $f^\varepsilon_n \in P_\varepsilon(n)$，且其 profile bound $g_n$ 满足 $g_n \lesssim n^3 g_{n-1}$。关键引理 3.4 给出 $\|g\|_{X^{1,k}} \leq C_k (n+1)^3 g_n K^{-1} \varepsilon^{-2n-1}$，揭示了 $X^{1,k}$-范数随阶数 $n$ 呈超几何增长 $\varepsilon^{-2n-1}$。 （3）初始层主导项提取：利用 $P f^\varepsilon_1 = 0$（因初值无宏观分量且碰撞守恒律保证 $P\Gamma = 0$），证明首阶修正纯微观，故宏观速度 $u[f^\varepsilon]$ 在 $O(t)$ 阶无贡献；而 $f^\varepsilon_2$ 的宏观投影 $P f^\varepsilon_2$ 生成 $O(t^2 \varepsilon^{-5})$ 径向速度，经 $\nabla_x \cdot$ 作用得 $O(t^2 \varepsilon^{-6})$ 散度，再通过质量/能量守恒方程导出 $O(t^3 \varepsilon^{-7})$ 密度/温度扰动。整个框架的成败依赖于取 $M \sim |\ln \varepsilon|$ 以平衡递归误差 $\sum_{n=0}^M \varepsilon^{-2n-1} t^n$ 在 $t \lesssim \varepsilon^2$ 下的收敛性，最终得到 $T_\varepsilon \sim \varepsilon^{-2} |\ln \varepsilon|^{-6}$。

4. 与之前工作的比较

本文与既有工作存在本质差异： （1）初值正则性维度：经典结果如 [12] 要求初值 $f^{\varepsilon}_{\mathrm{in}} \in X^{1+,k}$，隐含 $u[f^{\varepsilon}_{\mathrm{in}}] \in H^{1+\eta}_x$，而本文初值 $u^{\varepsilon}_{\mathrm{in}} \sim |x|^{-1}$ 仅属 $\mathcal{M}(\mathbb{R}^2)$（有限 Radon 测度），远低于 $L^2$ 阈值。这突破了 Boltzmann 理论中‘初值需强于流体解’的范式，首次在 kinetic 框架下处理 Navier-Stokes 允许的最奇异初值。 （2）时间尺度与失效机制：以往工作（如 [28]）通过 $O(1)$ 时间热演化预先平滑初值，避开初始层；本文则聚焦 $t \ll \varepsilon^2$ 的极端短时尺度，揭示 kinetic 系统在此尺度内无法建立流体所需的瞬时压力平衡，转而发展出可压缩的径向离心动力学。这种 $O(t^2/\varepsilon^5)$ 径向速度与 $O(t^2/\varepsilon^6)$ 散度下界，是此前所有 work 未观测到的定量失效模式。 （3）技术代价与适用范围：本文结论依赖于 $\varepsilon$-依赖正则化（$|x|_\varepsilon$ 与 $\chi(\varepsilon^\gamma |x|)$），若采用 $\varepsilon$-无关正则化则退化为 [28] 的 setting；且当前证明限定于二维硬球模型，虽 Remark 1.6 指出可推广至硬势，但软势情形因速度权重损失更小而技术上更简单，本文未覆盖。此外，$|\ln \varepsilon|^{-6}$ 时间尺度中的对数因子源于递归阶数截断，这是处理奇异性时不可避免的代价，而 [12] 的全局收敛无需此类修正。 （4）物理诠释深化：不同于 [7]（Onsager 关于无粘点涡的不可压质疑）或 [14]（Lamb-Oseen 的吸引性），本文首次从 kinetic 微观机制解释为何‘粘性正则化’在 Boltzmann 框架下并非瞬时：碰撞算子 $L$ 的弛豫时间 $O(\varepsilon^2)$ 与输运/碰撞非线性项的 $O(\varepsilon)$ 时间尺度竞争，导致在 $t \sim \varepsilon^2$ 前，宏观约束（不可压、Boussinesq）无法闭合，系统被迫经历一个 genuine compressible initial layer。

12. 投影-提升等价性与耗散紧致性：连续态FENE-Markov流体

Projection-Lift Equivalence and Dissipation-Tight Compactness for Continuum-State FENE-Markov Fluids

作者: Sai Peng

类别: math.AP, math-ph

发布日期: 2026-07-18

arXiv ID: 2607.16993v1

摘要翻译: 我们考虑不可压FENE哑铃模型，其耦合于一个定义在紧致连续内态空间上的可逆Markov算子。在质心扩散为零的情形下，我们证明态平均是一个精确的因子映射，且在自然能量解类上，Lagrangian态提升是其唯一的逆映射。因此，态分辨系统的存在性、多重性与唯一性恰好与其标量FENE投影一致；内态动力学不会引入额外的大初值非唯一性。该提升由一个迹为零的矩阵纤维演化所驱动，并允许非线性局域活性，包括对奇异Kramers应力具有线性依赖性的速率项。我们还证明了这一结构的序列形式：对于具有相对熵预置初始纤维且无粘性耗散缺陷的状态分辨正则化系统，其完整密度强收敛，无需在标量极限后重构它们。从而，完整的拖曳力、态依赖活性、奇异应力以及非原子Jeffreys产生均通过极限传递。该论证结合了正则Lagrangian流的稳定性与一个针对同时变化的矩阵漂移和跳跃率的相对熵估计。一个稳态振荡表明，该预置不能仅由自然熵界推出。当质心扩散为正时，我们独立构造了全局大初值弱解，并识别出相同的非线性项。一个显式的无限秩核证明了这些结果并非有限物种约化。

Comment: 无

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

1. 研究对象和背景

本文研究的是不可压Navier–Stokes方程与一类具有连续内部状态空间的FENE（Finite Extensible Nonlinear Elastic）哑铃模型的强耦合系统。该模型描述了聚合物溶液中溶剂流体与构型–状态联合概率密度的动力学演化，其核心难点在于：(i) FENE弹簧势在有限伸展边界∂D处产生奇异性；(ii) 内部状态空间Y为紧致度量空间（非原子、可能无限维），其上由可逆Markov核k(y, y′)驱动的跳跃动力学依赖于局部矩η_ah，而这些矩本身是弱紧致的密度h的非线性泛函；(iii) Kramers应力τ_Y[h] = ∫_{D×Y} (∇_qU ⊗ q) M h , dq dπ在h仅弱收敛时无法被识别，因∇_qU ⊗ q在∂D处发散，且M(q) ∼ dist(q, ∂D)^γ（γ > 1）呈幂律衰减。具体而言，系统由如下方程组封闭：

其中M(q) = Z^{-1} e^{-U(q)}为FENE Maxwellian，满足定义2.1的admissible finite-extension条件：存在s_0 > 0, γ > 1及正常数c_j，使得在边界邻域{s_D < s_0}内，c_1 s_D^γ ≤ M ≤ c_2 s_D^γ且c_3 ≤ s_D |∇q U| ≤ c_4；U ∈ C^2{loc}(D)，D ⊂ ℝ^d为有界C^2区域，0 ∈ D；δ ≥ 0为中心质心扩散系数；Wi > 0为Weissenberg数；非线性反应算子\mathcal{R}_Y h定义为

其中K_h(t,x,q,y,y’) = k(y,y’) a\big(q, ηh, y, y’\big)，而局部矩η_a[h]由式(6)给出：

关键未知量为速度场u(t,x) ∈ ℝ^d与联合密度h(t,x,q,y) ≥ 0，定义在时空域(0,T) × Ω × D × Y上，Ω为有界Lipschitz区域或环面，Y为紧致度量空间，π为其全支撑Borel概率测度。问题的重要性在于：传统有限种群（finite-species）FENE–Markov模型[9,10]无法刻画真实聚合物链的连续构型–状态谱；而现有理论[4–6]均建立在标量密度（即Y退化为单点）基础上，其证明严重依赖径向对称性、H^1紧致性及熵–能量估计的标量结构。当Y为连续非原子空间时，Jeffreys产额D_Y(Mh) = \frac{1}{2} ∫_{Ω×D×Y^2} M K_h (h’ - h)(\log h’ - \log h) , dπ dπ’ dq dx不再集中于离散点，其弱下半连续性需在固定参考测度π ⊗ π下验证（引理2.7）；同时，Kramers应力τ_Y[h]的识别要求对奇异张量W = ∇_qU ⊗ q进行加权Hardy估计（引理2.5）：

本文试图填补的核心空白是：在零中心质心扩散（δ = 0）这一物理上更相关但分析上更病态的极限下，能否建立连续态FENE–Markov系统的存在性、唯一性与紧致性理论？特别地，是否存在一种结构性等价——即连续态解集与标量FENE解集之间存在双射——从而将复杂连续态问题完全归约为已知的标量问题？这要求克服三大障碍：(1) δ = 0时Fokker–Planck方程退化为一阶输运–反应方程，缺乏扩散正则化；(2) 连续Y上的Markov跳跃率a依赖于弱紧致的h，导致非线性反馈循环；(3) 标量投影\bar{h}(t,x,q) = ∫_Y h , dπ与状态提升之间的可逆性需在自然能量类中严格构造，而非仅作形式推导。

2. 主要定理或主要结果

本文在三个不同扩散机制下建立了完整的理论框架，主要结果如下：

（1）零扩散（δ = 0）下的投影–提升等价性（定理4.20与推论4.21）：设Ω = 𝕋^d（d = 2,3）或有界C^2区域，Y为紧致度量空间，k满足(3)，活动率a满足(7)–(8)且为Lipschitz连续的局部有限矩函数。则状态平均映射\mathcal{P}: h ↦ \bar{h} = ∫_Y h , dπ是自然能量解类上的精确因子映射（exact factor map），其逆映射\mathcal{L}（Lagrangian状态提升）由DiPerna–Lions正则流X(t,a)与迹为零矩阵纤维演化B(t,a) = ∇_x u(t,X(t,a)) ∈ L^2(0,T)唯一确定。即：对任意标量FENE能量解\bar{h}，存在唯一连续态解h满足\mathcal{P}h = \bar{h}，且h沿几乎每条流线满足

其中B(t,a)满足tr B = 0。因此，解的存在性、多重性与唯一性完全等价于其标量投影，内部状态动力学不引入额外的大数据非唯一性。

（2）耗散紧致序列（定理4.24与推论4.25）：设{h_n}为δ = 0的正则化序列，初始纤维相对熵有界：∫_{D×Y} M h_n^0 \log(h_n^0 / \bar{h}n^0) , dq dπ ≤ C，且粘性耗散无缺陷：lim{n→∞} ∫_0^T ∫_Ω |∇_x u_n|^2 dx dt = ∫_0^T ∫_Ω |∇x u|^2 dx dt。则h_n在L^1((0,T)×Ω; L^1{M⊗π}(D×Y))中强收敛，且收敛速率为

由此，完整拖曳项τ_Y[h_n]、非线性活动率a(q,η[h_n],y,y’)、奇异Kramers应力∫ M W h_n及非原子Jeffreys产额D_Y(M h_n)均直接通过序列极限识别，无需先投影后提升。

（3）正扩散（δ > 0）下的全局弱解（定理3.4）：在d = 2,3维有界Lipschitz区域Ω上，对任意满足(7)–(8)的有界Lipschitz局部矩活动率a，及满足初值条件(49)的数据，系统存在全局大初值能量–熵弱解(u,h)。其中u ∈ L^∞(0,T; H_σ(Ω)) ∩ L^2(0,T; V_σ(Ω))，h ≥ 0满足Mh ∈ L^∞(0,T; L^1)、Mh log h ∈ L^∞(0,T; L^1)、√h ∈ L^2(0,T; H^1_x(L^2_{M⊗π}))、D_Y(Mh) ∈ L^1(0,T)。关键识别结果为：

且由加权Hardy估计(17)保证奇异应力的L^2_t,x有界性：‖∫{D×Y} M W h_n dq dπ‖{L^2_{t,x}} ≤ C_R^{1/2} \left( ∫ M(h_n + |∇_q √{h_n}|^2) \right)^{1/2}。

（4）无穷秩核的非有限种群性（命题2.3）：取Y = 𝕋^1，π为Lebesgue测度，k_∞(y,y’) = 1 + ∑{m=1}^∞ 2^{-m} cos(2πm(y-y’))，则对应生成元L∞v = ∫Y k∞(y,y’)[v(y’)-v(y)] dπ(y’)具有无穷多个互异非零特征值{−(1−2^{-|n|-1})}_{n∈ℤ{0}}，故为无穷秩算子，证伪了所有结果仅为有限种群约化的可能性。

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

本文的核心技术障碍在于δ = 0时Fokker–Planck方程丧失扩散正则化，导致标准紧致性方法失效，且连续态Y上的非线性反馈使传统逐点估计失效。作者的真正创新点在于构建了一套基于正则Lagrangian流与矩阵纤维演化的结构性分析框架，而非依赖泛泛的能量估计。

第一，迹为零矩阵纤维演化（Proposition 4.5）：这是全文最关键的分析输入。作者观察到，在δ = 0时，标量投影\bar{h}满足的输运方程中，拖曳项∇_x u q的Jacobi矩阵B = ∇_x u沿DiPerna–Lions流X(t,a)演化，且tr B = 0自动满足（由∇_x · u = 0）。作者据此构造纤维演化方程

并证明该方程在L^1_{M⊗π}中保持正性与精确L^1收缩性。其核心在于利用FENE Hardy估计（引理2.5）控制熵–拉伸平衡：

其中右端C依赖于‖B‖{L^2{t,x}}，而B的L^2有界性由标量解u的自然能量类保证。此估计结合Bielecki–Volterra不动点论证，首次在全时间区间上闭合了非线性局部活动。

第二，相对熵估计处理同步变化的漂移与跳变率（Theorem 4.24）：为证明耗散紧致序列的强收敛，作者发展了一种新型相对熵估计，用于比较两个纤维h_n与h，其漂移矩阵B_n与B及跳变核K_{h_n}与K_h均随时间同步变化。关键引理给出

其中误差项由a的Lipschitz性与η[h_n] − η[h]的L^2控制。该估计与正则Lagrangian流的稳定性（即X_n → X在C([0,T]; L^p)中）耦合，最终导出h_n → h的强L^1收敛。

第三，非原子Jeffreys产额的下半连续性（引理2.7）：针对Y上非原子测度π，作者采用Ioffe型固定系数积分下半连续性定理，而非经典凸泛函弱下半连续性。给定a_n → a in measure且(r_n,s_n) ⇀ (r,s) in L^1(μ;ℝ^2)，则

其中J(r,s) = (r−s)(log r−log s)在[0,∞)^2上联合凸且下半连续，其Hessian矩阵D^2J(r,s) = (r+s)\begin{pmatrix} r^{-2} & -(rs)^{-1} \ -(rs)^{-1} & s^{-2} \end{pmatrix}半正定。此性质是识别D_Y(Mh_n) → D_Y(Mh)的基石。

第四，无穷秩核的显式构造（命题2.3）：作者摒弃抽象存在性论证，显式构造k_∞(y,y’) = 1 + ∑{m≥1} 2^{-m} cos(2πm(y−y’))，并计算其在L^2_π(𝕋^1)上的谱：L∞ e_{2πiny} = −(1−2^{-|n|-1}) e_{2πiny}，n ≠ 0。无穷多互异非零特征值直接证伪有限种群还原，确保所有结果具有本质的连续态内涵。

4. 与之前工作的比较

本文与既有工作在多个维度上形成实质性推进，其差异与代价需分层剖析：

与标量FENE理论（[4–6]）的比较：Lions–Masmoudi [4]的corotational构造要求拖曳项为Ω(u)q（Ω为旋度），且仅适用于b > 2的Warner势；Barrett–Süli [5]的正扩散理论依赖δ > 0提供的H^1_x正则性，无法处理δ = 0的物理极限；Masmoudi [6]的缺陷传播法虽处理δ = 0，但仅针对标量密度，且其缺陷项‖∇x u_n‖^2{L^2} − ‖∇x u‖^2{L^2}的控制高度脆弱。本文彻底摆脱标量框架：在δ = 0下，通过投影–提升等价性（定理4.20）将连续态解集与标量解集一一对应，从而将[4–6]的所有标量结论（存在性、唯一性阈值、正则性）自动提升至连续态；代价是必须假设活动率a为Lipschitz连续的局部矩，而[6]允许更一般的有界 measurable a。

与有限种群FENE–Markov模型（[9,10]）的比较：文献[9,10]将Y离散化为有限点集，其反应算子\mathcal{R}Y为有限维矩阵作用，紧致性由有限维空间自动保证。本文的连续态Y导致\mathcal{R}Y作用于无限维L^1{M⊗π}空间，其非线性反馈a(q,η[h],y,y’)的弱紧致性成为核心障碍。本文通过命题2.3的显式无穷秩核k∞，严格证明结果非有限种群还原：若为有限种群，则生成元必为有限秩，而L_∞具有无穷秩，故二者线性共轭不可能。代价是必须引入相对熵初值准备（定理4.24），而[9,10]仅需L^1初值。

与一般连续态反应扩散理论（[11,12]）的比较：[11,12]建立在熵梯度流结构上，要求跳变率a仅依赖于全局平均，且Y为光滑流形以定义梯度。本文允许a依赖于局部矩η_ah，其定义涉及空间变量x，故为非局部反应；且Y仅为紧致度量空间，无需微分结构。本文的突破在于将熵梯度结构推广至状态敏感的局部反馈，并用引理2.2的连续细致平衡（continuum detailed balance）替代经典梯度流，代价是必须假设k(y,y’)对称且可逆（即π-几乎处处k > 0）。

与紧致性方法（Aubin–Lions, Kolmogorov–Riesz）的比较：标准紧致性依赖空间导数与时间导数的联合控制（引理2.4）。本文在δ = 0时，时间导数∂_t(Mh)仅在(W^{1,∞})′中有界，远弱于L^1_t(W^{-1,1}_x)，故标准Aubin–Lions失效。作者转而利用正则Lagrangian流的稳定性（[13]）与相对熵估计的双重控制，实现强收敛。代价是必须附加‘无粘性耗散缺陷’这一非平凡条件（推论4.25），而标准方法无需此假设。

与Hardy估计技术的比较：以往FENE文献（如[5]）仅使用径向Hardy估计。本文引理2.5将估计推广至非径向admissible Maxwellian（定义2.1），其核心是利用边界法坐标(p,s)与切向单位分解，将多维估计归约为一维加权Hardy型

该估计成立当且仅当γ > 1，恰好对应FENE势的物理参数范围（如Warner势b > 2 ⇒ γ = b/2 > 1）。此推广使理论适用于各向异性势U_A(q) = −b/2 log(1−q·Aq/b)，而此前工作仅限球对称情形。

13. 关于刚性球簇对粘性流体运动的影响

On the impact of clusters of rigid balls on the motion of a viscous fluid

作者: Marco Bravin, Eduard Feireisl, Arnab Roy, Arghir Zarnescu

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

发布日期: 2026-07-17

arXiv ID: 2607.16470v1

摘要翻译: 我们为大量刚性物体浸没于粘性流体中的运动问题发展了一种新方法。其核心思想是‘簇’（cluster）的概念——即一组个体刚性物体，它们可被聚类甚至连接，使得其对系统整体运动的集体影响类似于单个物体。该新方法的应用包括：

Comment: 无

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

1. 研究对象和背景

本文研究的是高密度刚性小球浸没于不可压缩牛顿粘性流体中所构成的流固耦合系统在动态均匀化（dynamic homogenization）极限下的渐近行为。核心对象是 $ d = 2 $ 或 $ 3 $ 维周期条形域

其中 $ \mathbb{T}^{d-1} $ 是 $ (d-1) $ 维平坦环面。流体占据 $ \Omega \setminus \bigcup_{n=1}^N B(r, h_n(t)) $，其中 $ B(r, h_n(t)) = \{ x \in \mathbb{R}^d \mid |x - h_n(t)| < r \} $ 表示半径为 $ r > 0 $、质心位置为 $ h_n(t) \in \Omega $ 的刚性球；$ N = N(r) \to \infty $ 随 $ r \to 0 $ 增长。流体速度 $ u_F = u_F(t,x) $ 与压力 $ p = p(t,x) $ 满足不可压缩 Navier–Stokes 方程：

其中黏性应力张量 $ S = S(D_x u_F) = \nu (\nabla_x u_F + \nabla_x^\top u_F) $，应变率张量 $ D_x u_F = \frac{1}{2}(\nabla_x u_F + \nabla_x^\top u_F) $，$ \nu > 0 $ 为动力黏度，$ \varrho_F > 0 $ 为流体密度，$ g = (0,\dots,0,-1) $ 为竖直重力场。边界条件包括：在球面 $ \partial B(r,h_n(t)) $ 上施加无滑移条件 $ u_F|_{\partial B(r,h_n)} = u_n^S $，其中刚体速度为

在上下刚性壁 $ x_d = 0,H $ 上满足不可穿透与无切向应力：$ u_F \cdot n = u_{F,d}|_{x_d = 0,H} = 0 $，且 $ [S \cdot n]_{\mathrm{tan}} = 0 $。刚性球具有恒定密度 $ 0 < \varrho_n^S \leq \varrho_\varepsilon $（允许随 $ \varepsilon \to 0 $ 缓慢增长），其运动由动量守恒控制：对 $ d = 3 $，

其中惯性张量 $ J_n(t)a \cdot b = \int_{B(r,h_n)} \varrho_n^S \big[ a \times (x - h_n) \big] \cdot \big[ b \times (x - h_n) \big] \, dx $。该系统构成强形式的流固耦合问题，其解在碰撞发生时可能失去经典正则性，故需采用 Judakov 型弱解框架（见 Definition 2.1）。

该问题的重要性在于它介于经典静态均匀化（如 Allaire [1,2] 对固定小孔的穿孔区域分析）与真实多体动力学之间：静态均匀化假设障碍物固定不动，其影响由容量（capacity）主导，临界数量满足 $ N \sim r^{-(d-2)} $（$ d \geq 3 $）或 $ N \sim |\log r|^{-1} $（$ d = 2 $）；而动态情形下，小球可自由运动、聚集甚至碰撞，其集体效应不再由几何容量决定，而是由动力学耦合强度、相对速度与局部聚集结构共同调控。已有工作 [12] 首次建立动态均匀化框架，证明当 $ N \lesssim |\log r| $ 时，小球云对极限 Navier–Stokes 系统无扰动，但该估计过于保守，且未刻画小球轨迹如何趋近流体速度场。更关键的空白在于：（i）缺乏对“小球是否跟随流体运动”的定量刻画——即 $ \| \dot{h}_n(t) - u(t, h_n(t)) \|_{L^2_t} $ 的收敛速率；（ii）无法处理高密度小球簇（cluster）内部频繁接触与粘连现象，传统单球分离假设失效；（iii）现有估计严重依赖维数特异性技术（如 $ d = 2 $ 的对数容量），缺乏统一尺度分析框架。本文旨在通过引入“簇”（cluster）这一新概念，将空间上邻近、动力学上强耦合的小球分组，将其等效为一个广义刚体，并发展适配该结构的投影算子与能量比较方法，从而系统性填补上述三项空白。

2. 主要定理或主要结果

本文在两个严格设定下建立了动态均匀化的最优渐近理论：

Theorem 2.2（弱极限）：设 $ (\varrho_\varepsilon, u_\varepsilon, (h_{n,\varepsilon})_{n=1}^{N_\varepsilon}, (O_{n,\varepsilon})_{n=1}^{N_\varepsilon}) $ 是流固耦合问题的弱解族，初始数据满足 $ |h_{i,\varepsilon}(0) - h_{j,\varepsilon}(0)| \geq 2r_\varepsilon $（无初始重叠），且 $ \sqrt{\varrho_{0,\varepsilon}} u_{0,\varepsilon} \to \sqrt{\varrho_F} u_0 $ 弱收敛于 $ L^2(\Omega; \mathbb{R}^d) $。若参数满足

则存在子列（仍记为 $ \varepsilon $）使得

且极限 $ (\varrho_F, u) $ 是不可压缩 Navier–Stokes 系统的弱解，满足

对所有 $ \varphi \in C_c^1([0,T) \times \Omega; \mathbb{R}^d) $，$ \operatorname{div}_x \varphi = 0 $，$ \varphi \cdot n|_{\partial \Omega} = 0 $ 成立。该结果将临界数量从 [12] 的 $ N \lesssim |\log r| $ 大幅提升至 $ N \ll r^{-d/5} |\log r|^{-4/5} $（$ d = 3 $），远超 Allaire 静态均匀化临界值 $ N \sim r^{-1} $（$ d = 2 $）。

Theorem 2.3（强极限与轨迹收敛）：在相同弱解族下，若初始位形满足 $ h_{n,\varepsilon}(0) \to h_{n,0} $ 一致，且初速度误差满足 $ \int_{\mathbb{T}^d} \varrho_\varepsilon |u_{0,\varepsilon} - u_0|^2 dx \leq r_\varepsilon^\alpha \beta(\varepsilon) $，其中 $ \beta(\varepsilon) \to 0 $，$ \alpha = 1 $（$ d = 3 $）或任意 $ \alpha > 0 $（$ d = 2 $）；并假设极限 Navier–Stokes 系统存在经典解 $ u \in C^1([0,T] \times \Omega) \cap C([0,T]; W^{2,\infty}(\Omega)) $，且参数满足

则有更强收敛：

且对每个 $ n $，$ h_{n,\varepsilon} \to h_n $ 在 $ W^{1,2}(0,T; \mathbb{R}^d) $ 中，且极限轨迹满足常微分方程

这是首个给出动态均匀化中 $ \| \dot{h}_n - u(\cdot, h_n) \|_{L^2_t} \to 0 $ 显式速率的结论。

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

本文的核心创新在于构建了一套基于“簇”（cluster）的动态均匀化分析框架，彻底摆脱了传统单球分离假设与静态容量判据的束缚。其技术障碍与突破可分解为以下三层次：

1. 簇的定义与投影算子（Section 4）：作者定义“簇距离” $ \mathrm{dist}_C(h_i, h_j) $ 为两球质心间沿流场特征线的加权距离，并据此将球集划分为极大连通子集 $ \mathcal{C}_k = \{ n \mid \mathrm{dist}_C(h_n, h_m) \leq \ell_\varepsilon \} $，其中 $ \ell_\varepsilon $ 是与 $ r_\varepsilon $ 和 $ N_\varepsilon $ 匹配的尺度参数。关键工具是“簇投影算子” $ \Pi_C $：对任意向量场 $ v $，定义 $ (\Pi_C v)(x) = \frac{1}{|\mathcal{C}_k|} \sum_{n \in \mathcal{C}_k} v(h_n) $ 当 $ x \in B(r_\varepsilon, h_n) $，并在流体区延拓为光滑截断。该算子满足 $ \| v - \Pi_C v \|_{L^2(\Omega)} \lesssim r_\varepsilon^{1/2} \| \nabla_x v \|_{L^2(\Omega)} + N_\varepsilon^{1/2} \ell_\varepsilon \| v \|_{W^{1,\infty}} $，精准量化了簇内速度离散对能量估计的贡献。

2. 相对能量方法的重构（Section 7）：为证明 $ \dot{h}_n \to u(t, h_n) $，作者不直接估计 $ |\dot{h}_n - u(t, h_n)| $，而是构造修正相对能量泛函

其中第二项是对每个簇 $ \mathcal{C}_k $ 内所有球的动能偏差求和。利用簇投影算子，将动量方程中的流体-球耦合项 $ \int_{\partial B(r,h_n)} (S - pI) \cdot n \cdot (\dot{h}_n - u) $ 控制为 $ \lesssim r_\varepsilon^{d-1} \| \nabla_x u \|_{L^\infty} \mathcal{E}_\varepsilon^{1/2} + \text{error terms} $，再结合 Gronwall 不等式导出 $ \mathcal{E}_\varepsilon(t) \leq C \beta(\varepsilon) e^{Ct} $，从而得到 $ \| \dot{h}_n - u(\cdot, h_n) \|_{L^2(0,T)} \to 0 $ 的显式速率。

3. 簇内碰撞的鲁棒处理（Section 4.3）：针对簇内球体可能发生的无限次碰撞，作者引入“簇惯性张量” $ J_C = \sum_{n \in \mathcal{C}_k} J_n $ 与“簇质心速度” $ V_C = \frac{1}{\sum_{n \in \mathcal{C}_k} m_n} \sum_{n \in \mathcal{C}_k} m_n \dot{h}_n $，将整个簇视为一个广义刚体。关键引理 4.5 证明：只要簇直径 $ \mathrm{diam}(\mathcal{C}_k) \leq \ell_\varepsilon $，则簇内相对速度满足 $ \max_{i,j \in \mathcal{C}_k} |\dot{h}_i - \dot{h}_j| \leq C \ell_\varepsilon^{1/2} \| \nabla_x u \|_{L^\infty} $，从而保证簇的整体运动由流场主导，内部微结构不影响宏观极限。这一框架使证明完全规避了“避免碰撞”或“瞬时分离”的强假设，实现了真正全局时间的无条件收敛。

4. 与之前工作的比较

本文与既有工作的差异体现在四个根本维度：

1. 临界数量尺度的跃升：对比 Allaire [1,2] 的静态均匀化，其临界数量由小球容量 $ \mathrm{cap}(B(r)) \sim r^{d-2} $（$ d \geq 3 $）或 $ \sim |\log r|^{-1} $（$ d = 2 $）决定，要求 $ N \cdot \mathrm{cap}(B(r)) \to 0 $，故 $ N \ll r^{-(d-2)} $（$ d \geq 3 $）或 $ N \ll |\log r| $（$ d = 2 $）。而本文 Theorem 2.2 给出 $ N \ll r^{-d/5} |\log r|^{-4/5} $（$ d = 3 $），比 Allaire 的 $ r^{-1} $ 提升整整 $ r^{-2/5} $ 阶；在 $ d = 2 $ 时，$ N \ll r^{-2/5} |\log r|^{-4/5} $ 远超 $ |\log r| $。这一进步源于放弃容量判据，转而用簇投影控制速度离散，将关键障碍从几何测度转为动力学 Lipschitz 常数。

2. 轨迹收敛的首次定量刻画：Höfer–Schubert [18] 与 Jabin–Otto [19] 研究沉降问题，但仅得 $ h_n(t) \to h_\infty $（稳态位置），未涉及 $ \dot{h}_n $ 与 $ u(t,h_n) $ 的瞬时匹配；[12] 仅证明弱极限存在，未给出任何轨迹信息。本文 Theorem 2.3 首次建立 $ \| \dot{h}_n - u(\cdot, h_n) \|_{L^2_t} \to 0 $ 的显式速率，且该速率与初值误差 $ \beta(\varepsilon) $ 和簇尺度 $ \ell_\varepsilon $ 直接相关，填补了动态均匀化理论中“运动学一致性”的根本空白。

3. 假设的实质性放宽：[12] 要求小球密度 $ \varrho_n^S \to \infty $（类固定孔洞），而本文允许 $ \varrho_n^S \leq \varrho_\varepsilon $ 且 $ \varrho_\varepsilon $ 可缓慢增长（如 $ \varrho_\varepsilon \sim |\log r_\varepsilon|^{\sigma} $）；[12] 排除碰撞，本文明确允许并处理无限次碰撞；[12] 限于 $ d = 3 $，本文统一处理 $ d = 2,3 $，且 $ d = 2 $ 结果不依赖对数容量技巧，更具普适性。

4. 代价与局限：本文的簇框架以引入额外参数 $ \ell_\varepsilon $ 为代价，其最优选择需平衡 $ r_\varepsilon $ 与 $ N_\varepsilon $（见 (2.12)(2.17)），导致临界数量表达式较复杂；目前仅适用于周期域 $ \mathbb{T}^{d-1} \times (0,H) $，尚未推广到带边界的有界域（如 $ C^2 $ 区域），因簇投影在边界附近的截断会引入额外边界层误差；此外，相对能量估计依赖极限解 $ u $ 的 $ W^{2,\infty} $ 正则性，在 $ d = 3 $ 时仅能保证局部时间存在，全局强解仍受限于 Navier–Stokes 千禧年难题。

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

1. 非局部Trudinger方程的Harnack估计

Harnack estimates for the nonlocal Trudinger equation

作者: Simone Ciani, Kenta Nakamura

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17912v1

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

摘要翻译

中文标题: 非局部Trudinger方程的Harnack估计

摘要翻译: 我们继续研究具有可测且有界奇异核的非局部Trudinger方程的逐点正则性理论。我们在最优尾部条件下建立了精细的定量上界估计、弱抛物Harnack不等式和强抛物Harnack不等式。我们的分析依赖于对De Giorgi-Moser方法论的精细化适配，该适配基于Di Benedetto式的抛物方法，专门设计用于克服非局部、双重非线性框架下的技术复杂性。

2. 若干随机色散方程的几乎必然空间衰减与几乎必然非线性光滑化

Almost sure spatial decay and almost sure nonlinear smoothing of some stochastic dispersive equations

作者: Xiangqian Yan, Yongsheng Li, Jianhua Huang, Wei Yan

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17863v1

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

摘要翻译

中文标题: 若干随机色散方程的几乎必然空间衰减与几乎必然非线性光滑化

摘要翻译: 本文研究随机mKdV方程与随机三次KdV-Benjamin-Ono方程的几乎必然非线性光滑化、几乎必然空间衰减以及几乎必然一致收敛性。首先，对初值 $g\in H^{s}(\mathbb{R})(s\geq\frac{1}{4})$ 与 $Φ_{2}\in L_{2}^{0,s}$，我们证明了随机三次KdV-Benjamin-Ono方程的局部适定性。其次，我们建立了随机mKdV方程与随机三次KdV-Benjamin-Ono方程的几乎必然非线性光滑化。最后，借助几乎必然非线性光滑化，我们得到了随机mKdV方程与随机三次KdV-Benjamin-Ono方程的路径解中积分项的几乎必然空间衰减与几乎必然一致收敛性。更精确地，我们有如下结果：对随机mKdV方程，设 $s>\frac{1}{3}$，$f\in H^{s}(\mathbb{R})$ 且 $Φ_{1}\in L_{2}^{0,s}$，则局部路径解 $u$ 满足 \begin{eqnarray*} &&\mathbb{P}\Big(\Big{ω: \lim_{t\rightarrow0}\Big|u-U(t)f-\int_{0}^{t}U(t-s)Φ_{1}dW(s)\Big|{L{x}^{\∞}}=0\Big}\Big)=1,\ &&\mathbb{P}\Big(\Big{ω: \forall t\in[0,T_ω], \lim_{|x|\rightarrow\infty}\Big(u-U(t)f-\int_{0}^{t}U(t-s)Φ_{1}dW(s)\Big)=0\Big}\Big)=1. \end{eqnarray*} 对随机三次KdV-Benjamin-Ono方程，设 $s>\frac{1}{3}$，$g\in H^{s}(\mathbb{R})$ 且 $Φ_{2}\in L_{2}^{0,s}$，则局部路径解 $v$ 满足 \begin{eqnarray*} &&\mathbb{P}\Big(\Big{ω: \lim_{t\rightarrow0}\Big|v-V(t)g-\int_{0}^{t}V(t-s)Φ_{2}dW(s)\Big|{L{x}^{\∞}}=0\Big}\Big)=1,\ &&\mathbb{P}\Big(\Big{ω: \forall t\in[0,T_ω],\lim_{|x|\rightarrow\infty}\Big(v-V(t)g-\int_{0}^{t}V(t-s)Φ_{2}dW(s)\Big)=0\Big}\Big)=1. \end{eqnarray*}

3. Sobolev 迹不等式的尖锐梯度稳定性

Sharp Gradient Stability for the Sobolev Trace Inequality

作者: Xi-Nan Ma, Yitian Zhang, Yang Zhou

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17588v1

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

摘要翻译

中文标题: Sobolev 迹不等式的尖锐梯度稳定性

摘要翻译: 设 $n\ge3$ 且 $1<p<n$。我们证明了上半空间上临界 Sobolev 迹不等式的定量稳定性估计。更精确地，Sobolev 迹亏量控制迹气泡流形的梯度距离的 $\max\{2,p\}$ 次幂。证明的核心部分是迹气泡的谱非退化性：线性化加权 Steklov 问题的前两个特征子空间恰好由振幅、缩放及切向平移模构成。

4. 分形 Turán-Nazarov 不等式与薛定谔方程的可观测性

Fractal Turán-Nazarov Inequality and Observability for Schrödinger Equations

作者: Jiaqi Yu, Shanlin Huang

类别: math.AP, math.CA

发布日期: 2026-07-20

arXiv ID: 2607.17505v1

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

摘要翻译

中文标题: 分形 Turán-Nazarov 不等式与薛定谔方程的可观测性

摘要翻译: 本文建立了薛定谔方程在分形集上可观测性不等式与唯一延拓性的限制。我们证明：与热方程不同，此类性质在分形环境下可能失效。为达成此目标，我们首先将经典的 Turán–Nazarov 不等式推广至分形情形；该经典不等式对形如 $\sum_{k=1}^nc_ke^{2πim_kt}$ 的三角多项式在正测度集合上给出下界估计。与经典情形不同，该不等式中的常数在次数 $n$ 上不再具有一致性，我们得到了依赖于 $n$ 与频率差 $m_n-m_1$ 的 sharp 界。这些精细估计进而使我们得以构造显式反例，表明当观测集为分形集时，薛定谔方程的可观测性与唯一延拓性可能失效。

5. Fefferman–Stein型估计与Fourier Sobolev空间中的分数阶非线性薛定谔方程

Fefferman–Stein-type estimates and fractional NLS in Fourier Sobolev spaces

作者: Divyang G. Bhimani, Ryosuke Hyakuna

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17546v1

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

摘要翻译

中文标题: Fefferman–Stein型估计与Fourier Sobolev空间中的分数阶非线性薛定谔方程

摘要翻译: 本文证明了分数阶薛定谔方程的一个sharp Fefferman–Stein型估计，该估计可视为定义在Fourier Lebesgue空间 $\widehat{L^p}$ 上初值数据的广义Strichartz估计。随后，作为Fefferman–Stein不等式及其非对角线推广的应用，我们在齐次与非齐次Fourier–Sobolev空间 $\widehat{\dot{H}^s_p},\widehat{H^s_p}$ 中，证明了一维分数阶非线性薛定谔方程在纯幂次非线性情形下的大初值局部适定性与小初值全局适定性结果。解被构造于 $L_x^r(\mathbb{R} ;L^q_t(I))$ 空间中，以克服标准Strichartz估计中导数损失的困难。

6. 在Lorentz流形上半线性波动方程中光滑非线性项的唯一反演

Unique inversion of smooth nonlinearity in semilinear wave equations on Lorentzian manifolds

作者: Xi Chen, Shuai Lu, Ruochong Zhang

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17501v1

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

摘要翻译

中文标题: 在Lorentz流形上半线性波动方程中光滑非线性项的唯一反演

摘要翻译: 我们在全局双曲Lorentz流形上研究半线性波动方程中光滑非线性项的反演问题。我们的方法结合高阶线性化、畸变平面波与高斯光束，利用三波、四波及 m 波相互作用来恢复非线性项的不同Taylor系数。

7. 平面Ginzburg–Landau方程整体解的有限势能

Finite Potential Energy for Entire Solutions of the Planar Ginzburg–Landau Equation

作者: Hongge Chen, Juncheng Wei, Haicheng Yan, Wen Yang

类别: math.AP, math-ph

发布日期: 2026-07-20

arXiv ID: 2607.17490v1

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

摘要翻译

中文标题: 平面Ginzburg–Landau方程整体解的有限势能

摘要翻译: 我们证明：对Ginzburg–Landau方程 $ -Δu=u(1-|u|^2) $ 的任意光滑整体解 $ u\colon\mathbb{R}^2\to\mathbb{R}^2 $，若满足 $ |u(x)|\to1 $ 当 $ |x|\to\infty $，则其势能有限，即 \begin{equation*} \int_{\mathbb{R}^2}(1-|u|^2)^2\dd x<+\infty, \end{equation*} 从而解决了 \cite{BrezisProblems} 中Brezis提出的公开问题2.5。主要困难源于可能存在一个无旋模态（curl-free mode），该模态携带非零环量且仅以 $ |x|^{-1} $ 速率衰减；此类模态不属于 $ L^2 $ 空间，亦不具有单值势函数。我们通过对 $ L^2 $ 梯度修正项取极小，构造了一个比较场，该场满足齐次方程并继承相同的环量。结合Kelvin反演与拟线性椭圆方程的De Giorgi–Nash–Moser理论，我们得到最优衰减估计 $ O(|x|^{-1}) $。对于Ginzburg–Landau解，Bernstein估计与Jacobi型的强制性（coercivity）在外部相位方程中产生一个 $ L^2 $ 型外力项。由此导出的相位场 $ L^4 $ 有界性蕴含 $1-|u|^2\in L^2(\mathbb{R}^2)$，故势能有限。

8. 超扩散型积分微分方程的临界阈值与瞬时范数膨胀

Critical thresholds and instantaneous norm inflation for super-diffusive integro-differential equations

作者: Bruno de Andrade

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17430v1

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

摘要翻译

中文标题: 超扩散型积分微分方程的临界阈值与瞬时范数膨胀

摘要翻译: 本文研究定义在 $\mathbb{R}^N$ 上的一类非线性积分微分方程的柯西问题，该方程刻画反常超扩散输运。其线性动力学由一个双尺度记忆核驱动，该核的拉普拉斯变换具有扇形性，并在高频与低频处呈现不同的幂律渐近行为。这种超扩散结构排除了经典抛物理论所具有的无穷正则化能力；因此，相应的预解算子在傅里叶空间中具有重代数尾部，作为 Hörmander 类 $S^{-2}_{1,0}$ 中的拟微分算子作用，从而限制空间光滑性。通过建立严格的 $L^q-L^p$ 乘子估计，确定了局部适定性的临界勒贝格阈值 $q_c$。为证明该阈值的尖锐性，在超临界区域 $1 < q < q_c$ 中证明了瞬时范数膨胀——并由此导出不适定性。此外，追踪向长时松弛参数的结构相变，解析了全局渐近动力学。识别出非局部 Fujita 型临界指数 $ρ_F$，并在初始数据属于交集空间且足够小时，建立了整体存在性及代数衰减；前提是该非线性项保持超临界性，并克服连接双尺度的结构性代数障碍。该一般性框架可直接应用于典型物理模型，包括 Cole-Cole 分数迟滞模型与多尺度 Prabhakar 记忆模型。

9. 具有间断通量的多维标量守恒律：无需非退化性条件的适定性

Multidimensional scalar conservation laws with discontinuous flux: well-posedness without non-degeneracy

作者: Darko Mitrovic

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17310v1

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

摘要翻译

中文标题: 具有间断通量的多维标量守恒律：无需非退化性条件的适定性

摘要翻译: 我们在不施加物理时空通量的非退化性条件的前提下，建立了具有间断异质通量的多维标量守恒律的存在性、唯一性及局部 $L^1$-稳定性。间断集可为一族局部有限的 $C^2$ 超曲面，其交集与奇点包含于一个 $(d-1)$-维 Hausdorff 测度为零的闭集之中。 主要障碍在于：当时空通量分量的一个非平凡线性组合在状态区间上为常数时，多维紧性可能失效。我们证明，对于界面问题，该障碍可归约为物理法向通量；其平坦区间通过一个保持法向通量（从而保持 Rankine–Hugoniot 关系）的法向通量商被坍缩。在非平坦区域，时间分量与切向分量被替换为法向通量的辅助多项式函数，而物理法向分量保持不变。所得辅助时空向量满足 Panov 紧性理论所要求的区间非退化性。这导出了商变量的强紧性与强单侧迹。 界面容许性条件由与两侧物理法向通量相关的 vanishing-viscosity germ 的投影给出。局域粘性比较表明，商迹属于该极大 $L^1$-耗散 germ；而 Young 测度收缩论证则恢复了物理状态的强收敛性。对弯曲界面，则通过压平、局域化与有限传播拼接进行处理。该理论首先对 $BV$ 初值数据构造，再通过连续性延拓至取值于不变区间内的 $L^1_{\mathrm{loc}}\cap L^\infty$ 初值数据。

10. 相场模型、尖锐界面极限及接触线动力学的数值格式

Phase-Field Models, Sharp Interface Limits, and Numerical Schemes for Contact Line Dynamics

作者: Guosheng Fu, Yuan Gao, Jian-Guo Liu

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17410v1

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

摘要翻译

中文标题: 相场模型、尖锐界面极限及接触线动力学的数值格式

摘要翻译: 我们在一个统一的变分框架下研究液滴在固体基底上的接触线动力学的相场模型与尖锐界面模型。接触线是液体、气体与固体三相交汇处，其运动由于无滑移流体力学的经典应力奇异性而在连续介质建模中构成根本性困难。相场模型通过引入厚度为 $δ$ 的薄过渡层来正则化该奇异性，并借助附加于基底上的壁能（wall energy）所增强的 Ginzburg-Landau 自由能来编码界面效应。从总自由能 $E = E_b + E_w$ 出发，我们分析两类相场模型：Allen-Cahn 方程与 Cahn-Hilliard 方程。利用匹配渐近展开法（matched asymptotic expansions）在 $δ\to 0$ 极限下，我们恢复了它们各自对应的尖锐界面极限。在 Allen-Cahn 情形中，极限给出平均曲率驱动的运动，并伴随一个由动态接触角偏离 Young 角所驱动的接触线律；在 Cahn-Hilliard 情形中，极限导出一个 Mullins-Sekerka 问题，其接触线动力学具有相同形式。本工作的核心结果在于识别出跨越这两类模型的一致梯度流结构：Allen-Cahn 动力学对应于 $L^2$-梯度流，而 Cahn-Hilliard 动力学对应于 $H^{-1}$-梯度流；且二者均收敛至保持相同能量耗散结构的尖锐界面演化。这为接触线运动提供了一个基于单一变分原理的统一诠释。最后，我们基于极小化运动原理（minimizing movement principle）发展了能量稳定的数值格式，并建立了全离散问题的离散能量耗散性与适定性。数值算例证实：两种格式均松弛至同一稳态尖锐界面解，而其动力学行为则分别反映了不同的耗散机制。

11. 一维双线性粗糙奇异积分的端点判据

Endpoint Criteria for One-Dimensional Bilinear Rough Singular Integrals

作者: Binwei Dan, Qingying Xue

类别: math.CA, math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17207v1

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

摘要翻译

中文标题: 一维双线性粗糙奇异积分的端点判据

摘要翻译: 我们证明了一维双线性粗糙奇异积分的端点定理。我们的出发点是对相应角乘子的尖锐结构刻画。对任意均值为零的 $Ω\in L^1(\mathbb{S}^1)$，与 $T_Ω$ 关联的有限部角乘子具有有界变差当且仅当 $Ω$ 的对径偶部属于 $H^1(\mathbb{S}^1)$。该刻画精确识别出一维双线性情形中所需的旋转正则性。它还导出一个与双线性 Hilbert 变换的一致估计相容的 Stieltjes 分解。接着，我们在临界核假设下建立了两个有界性判据。第一，若 $Ω\in L\log L(\mathbb{S}^1)$，则 $T_Ω$ 从 $ L^{p_1}(\mathbb{R})\times L^{p_2}(\mathbb{R})$ 到 $ L^p(\mathbb{R})$ 有界，只要 $1<p_1,p_2,p<\infty$ 且 $ \frac{1}{p}=\frac{1}{p_1}+\frac{1}{p_2}$。此外，对数指数 $1$ 在尺度 $L(\log L)^A$ 中是最优的。第二，在临界方向指标处，若 $Ω\in\mathcal{K}_{1/2,β}(\mathbb{S}^1)$ 且 $β>\frac{3}{2}\max\bigl\{p_1,p_1',p_2,p_2'\bigr\}-1$，则相同的有界性成立。这两个临界核类互不包含。$L\log L$ 结果是通过将乘子约化为一个有界变差的有限部角剖面而获得的；方向性结果则源于端点 Fourier 衰减、乘积小波分解及插值。

12. 二维空间中拟线性电导率的偏数据Calderón问题

Partial data Calderón problem for quasilinear conductivities in dimension 2

作者: Ruirui Wu, Ting Zhou

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17126v1

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

摘要翻译

中文标题: 二维空间中拟线性电导率的偏数据Calderón问题

摘要翻译: 本文证明了二维空间中拟线性电导率偏数据Calderón问题的唯一性结果。证明基于高阶线性化方法以及在边界一部分上消失的二维复几何光学（CGO）解。由于解的导数出现在积分恒等式中，我们需要改进余项估计；为此，我们引入相位选取方式的修正，其作用类似于极限Carleman权。此外，我们分析了相位组合如何在解的乘积中产生特定模式，从而使得我们既能应用稳相近似（stationary phase）也能应用非稳相近似（nonstationary phase）来重构电导率。

13. 双曲线一支的海森堡唯一性对：平移格点十字形的超临界非唯一性

Heisenberg Uniqueness Pairs for a Hyperbola Branch: Supercritical Nonuniqueness for Shifted Lattice Crosses

作者: Zhiqiang Wan

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

发布日期: 2026-07-19

arXiv ID: 2607.17159v1

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

摘要翻译

中文标题: 双曲线一支的海森堡唯一性对：平移格点十字形的超临界非唯一性

摘要翻译: 我们研究正双曲线一支与平移格点十字形在超临界区域 $q=\alpha\gamma>1$ 下的海森堡唯一性。我们解决了 Giri 与 Manna 提出的任意平移问题中的无限维性断言：对两支上的任意平移，归一化预零化子空间是无限维的。更精确地，每个 $v\in BV((1,q))$ 均存在一个全局 $BV$ 预零化延拓；该延拓在两个扭转相位均平凡时具有唯一性，此时其不确定性为一维。证明将零化条件约化为一个扭转 Perron–Frobenius 算子的图方程，并结合相位一致的 Lasota–Yorke 估计与谱边缘刚性。我们还给出了整个 $L^1$ 预零化子空间的精确算子论规范型，以关联 Green 级数的最大收敛域来刻画。记 $Q$ 为扭转乘积算子、$A$ 为强迫算子，我们证明：$Q$ 在 $L^1((0,1))$ 上的谱为闭单位圆盘，$\Ran(I-Q)$ 非闭，且对除可数个代数值 $q>1$ 外的所有 $q$，算子 $\sum_{j=0}^{N-1}Q^jA:L^1((1,q))\to L^1((0,1))$ 的范数恒为 $2N$（对每个 $N\ge1$）。

14. 仅在 Hessian 较大区域成立的椭圆方程估计

Estimates on elliptic equations that hold only where the Hessian is large

作者: Amit Kumar Acharya, Ram Baran Verma

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17045v1

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

摘要翻译

中文标题: 仅在 Hessian 较大区域成立的椭圆方程估计

摘要翻译: 本文中，我们建立了如下一类退化完全非线性椭圆方程的粘性解的 Hölder 正则性：

其中该算子仅在 Hessian $ D^2u $ 足够大的区域中满足椭圆性。此类方程自然出现在自由边界问题及具有部分椭圆性的模型中。证明结合了一种修正的尖点函数与点对测度论证中接触集的分解。作为推论，在自然的结构假设下可推出内 Hölder 连续性。

15. 源自负幂次 Alt-Phillips 方程的一类退化单相自由边界问题

A Degenerate One-Phase Free Boundary Problem Arising From the Alt-Phillips Equation for Negative Powers

作者: Antonio Farah

类别: math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17009v1

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

摘要翻译

中文标题: 源自负幂次 Alt-Phillips 方程的一类退化单相自由边界问题

摘要翻译: 我们研究一类形如 $wΔw = h(\nabla w)$ 的退化单相自由边界问题的粘性解。我们假设存在一个星形区域 $D$，使得 $h < 0$ 在 $D$ 内成立，$h = 0$ 在 $\partial D$ 上成立，且 $h > 0$ 在 $\bar{D}^{c}$ 中成立。此类退化单相自由边界问题出现在对半线性方程 $Δu = f(u)$ 实施标准变换时，其中 $f$ 在某 $γ\in (0,2)$ 下表现如 $-γu^{-(γ+ 1)}$。此时该方程被称为负幂次势的 Alt-Phillips 方程，且 $h(ρ) = c(|ρ|^2 - 1)$。我们证明了粘性解的存在性、Lipschitz 正则性，以及自由边界在平坦点处的正则性。此外，我们证明当 $γ$ 趋于 $2$ 时，自由边界收敛于极小曲面。

16. Dirichlet 热流在边界-储库传输度量下的尖锐连续模

Sharp Continuity Moduli for Dirichlet Heat Flow in Boundary-Reservoir Transport Metrics

作者: Maja Gwozdz

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16874v1

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

摘要翻译

中文标题: Dirichlet 热流在边界-储库传输度量下的尖锐连续模

摘要翻译: 设 $Ω\subset\mathbb R^n$ 为有界 $C^2$ 开集，$P_t$ 为被杀死的 Dirichlet 热半群。我们证明了 $P_t$ 关于 Figalli–Gigli 边界-储库传输距离 $W_{b,p}$ 的尖锐固定时间幂尺度连续模。对任意 $t>0$，$P_t$ 关于 $W_{b,1}$ 全局 Lipschitz 连续；对任意 $p>1$，在总质量子水平集 $\{μ:μ(Ω)\le m\}$ 上，$P_t$ 是 $1/p$-Hölder 连续的：

对 $p>1$，我们证明指数 $1/p$ 在幂模尺度下是最优的。在全有限测度空间上，$P_t$ 在零测度处不连续。为建立下界，我们依赖于边界层的放大效应：更精确地，初始置于距 $\partialΩ$ 距离 $\varepsilon$ 处的单位质量，其输入 $W_{b,p}$-距离为 $O(\varepsilon)$，而经过任意固定的正时间后，其 $p$ 次边界矩有下界 $c\varepsilon$。因此，在二次情形及原始有限测度 $W_{b,2}$ 度量下，不存在一个定义在 $W_{b,2}$-度量域上的标准有限-$λ$ $\mathrm{EVI}_λ$ 半群，该半群既包含仿射常边界数据类，又能限制为仿射常边界 Dirichlet 热流。最后，我们亦描述了光滑一致椭圆散度型扰动所对应的下界障碍。

17. 具有非紧支集的Cucker-Smale型模型的一致时间稳定性与平均场极限

Uniform-in-time stability and mean-field limit of the Cucker-Smale-type model with noncompact support

作者: Seung-Yeal Ha, Xinyu Wang, Wook Yoon

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16794v1

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

摘要翻译

中文标题: 具有非紧支集的Cucker-Smale型模型的一致时间稳定性与平均场极限

摘要翻译: 我们研究一个无穷维Cucker–Smale型（ICS型）模型在完全非紧设定下的一致时间稳定性与平均场极限。首先，我们阐明了当空间支集为非紧时，原始ICS模型及其对应的动理学Cucker–Smale（KCS）模型所面临的主要困难。特别地，我们证明速度直径可能在时间上保持恒定，从而使得基于位置直径与速度直径的经典方法失效。此外，我们构造了一个反例，表明原始KCS模型在非紧空间设定下不满足一致时间稳定性。为克服这些障碍，我们引入了一类新的ICS型模型，其通讯权重依赖于成对空间矩。在适当假设下，我们亦建立了该ICS型模型在非紧设定下关于初始数据的一致时间稳定性。最后，我们推导出该ICS型模型的一致时间平均场极限，并同时得到相应KCS型模型的一致时间稳定性。

18. 在洛伦兹规范下陈-西蒙斯规范化的 $O(3)$ σ 模型的几乎最优适定性

Almost optimal well-posedness for Chern–Simons gauged $O(3)$ sigma model under the Lorenz gauge

作者: Huali Zhang, Jie Zhou

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16792v1

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

摘要翻译

中文标题: 在洛伦兹规范下陈-西蒙斯规范化的 $O(3)$ σ 模型的几乎最优适定性

摘要翻译: 本文研究在洛伦兹规范下 $\mathbb{R}^{1+d}$（$d=1,2$）中陈-西蒙斯规范化的 $O(3)$ σ 模型的低正则性柯西问题。 对于 $d=1$，我们证明了初值 $(\boldsymbolφ_0,\mathbf{A}_0)\in H^{s_1}(\mathbb{R})\times H^{s_1-1}(\mathbb{R})$ 且 $s_1>\frac12$ 时的局部适定性。该结果比 Jin 与 Huh \cite{HJ} 的先前结果提升 $\frac14$ 阶导数，并在标度不变正则性 $\dot H^{1/2}(\mathbb{R})$（物质场）与 $\dot H^{-1/2}(\mathbb{R})$（规范场）意义下几乎最优。 对于 $d=2$，我们证明了初值 $(\boldsymbolφ_0,\mathbf{A}_0)\in H^{s_2}(\mathbb{R}^2)\times H^{s_2-\frac34}(\mathbb{R}^2)$ 且 $s_2>1$ 时的局部适定性。该结果比 Jin 与 Zhang \cite{JZ} 的先前结果提升 $\frac14$ 阶导数，并使正则性阈值逼近标度不变指数 $\dot H^{1}(\mathbb{R}^2)$（物质场）与 $\dot H^{0}(\mathbb{R}^2)$（规范场）。 分析依赖于两个主要要素：在二维空间中，我们识别出导数非线性项的完整零结构，使得整个系统可在统一的零形式框架内处理；在一维空间中，我们在 Keel 与 Tao 引入的函数空间中建立直接的能量估计，避免了有限传播性约化至小数据问题，并使得对一般初值的低正则性迭代成为可能。

19. 定量可 rectifiability、奇异积分与调和函数边值问题之间的相互作用

Interactions between quantitative rectifiability, singular integrals, and boundary value problems for harmonic functions

作者: Xavier Tolsa

类别: math.CA, math.AP

发布日期: 2026-07-17

arXiv ID: 2607.16457v1

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

摘要翻译

中文标题: 定量可 rectifiability、奇异积分与调和函数边值问题之间的相互作用

摘要翻译: 本文综述了定量可 rectifiability 理论起核心作用的若干不同课题。首先，它回顾了以涉及 $β$ 型系数的平方函数为工具的可 rectifiability 刻画，以及 Carleson 的 $\varepsilon^2$ 猜想。其次，本文讨论了可 rectifiability 与 Riesz 变换在 $L^2$ 上的有界性之间的深刻联系，及其在 Lipschitz 调和函数的 Painlevé 问题中的应用。最后，本文探讨了关于调和测度以及 Laplace 方程在粗糙区域上的 Dirichlet、正则性（regularity）与 Neumann 问题的 $L^p$ 可解性的若干最新重大进展，并强调定量可 rectifiability 在这些发展中所起的关键作用。

不相关论文（基本信息）

1. 自适应 Mamba 神经算子

Adaptive Mamba Neural Operators

作者: Zeyuan Song, Zheyu Jiang

类别: cs.LG, math.AP, math.NA

发布日期: 2026-07-20

arXiv ID: 2607.18043v1

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

摘要: Accurately solving partial differential equations (PDEs) on arbitrary geometries and a variety of meshes is an important task in science and engineering applications. In this paper, we propose Adaptive Mamba Neural Operators (AMO), which integrates reproducing kernels for state-space models (SSMs) rather than the kernel integral formulation of SSMs. This is achieved by constructing Takenaka-Malmquist systems for the PDEs. AMO offers new representations that align well with the adaptive Fourier decomposition (AFD) theory and can approximate the solution manifold of PDEs on a wide range of geometries and meshes. In several challenging benchmark PDE problems in the fields of fluid physics, solid physics, and finance on point clouds, structured meshes, regular grids, and irregular domains, AMO consistently outperforms state-of-the-art solvers in terms of relative $L^2$ error. Overall, this work presents a new paradigm for designing explainable neural operator frameworks.

2. 关于细长圆柱体扭转函数的梯度

On the gradient of the torsion function for elongated cylinders

作者: Michiel van den Berg

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.18055v1

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

摘要: Let $(0,L)\times E_a\subset \R^{d+1}$ denote the cylinder of length $L$, and base $E_a\subset \R^{d},$ where $E_a$ is an open ellipsoid with semi-axes $a=(a_1,...,a_d)$. Let $w_{(0,L)\times E_a}$ denote its torsion function. Heat equation tools are used to show that (i) if $E_a$ is has small eccentricity and $L$ is sufficiently large, then the maxima of $|\nabla w_{(0,L)\times E_a}|$ are located at the centres $(0,0)$ and $(L,0)$ respectively, (ii) if $E_a$ has large eccentricity and $L$ is sufficiently large, then the maxima are located on the lateral side of the cylinder.

3. 正则化最优传输的尖锐渐近性

Sharp Asymptotics for Regularized Optimal Transport

作者: Carlos Cardoso-Perelló, Alberto González-Sanz, Marcel Nutz

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

发布日期: 2026-07-20

arXiv ID: 2607.18191v1

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

摘要: We study the small-regularization limit for $L^p$-regularized optimal transport with $1<p<\infty$ and for entropically regularized optimal transport (EOT). The exact first-order (respectively, second-order) asymptotics are determined explicitly under mild assumptions on the source and target measures. Our work generalizes the existing results for quadratic and entropic regularization, and connects them by a natural interpolation via $p\in(1,2)$. We derive all these asymptotics in a unified manner by a novel approach that separates the local computation of the optimal profile from the global enforcement of the marginal constraints: convex duality leads to Gaussian profiles for entropy and Barenblatt profiles for $L^p$-regularization, while a quantization construction turns these local profiles into couplings.

4. 关于 $ abla^6$-型 $Q$-曲率在 $ abla^6$-型共形度量中的正性与非正性结果

Positivity and non-positivity results for the sixth-order $Q$-curvature of conformal metrics in $\mathbb{R}^n$

作者: Jérôme Vétois, Samuel Zeitler

类别: math.DG, math.AP

发布日期: 2026-07-20

arXiv ID: 2607.18205v1

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

摘要: Given $n,m\in\mathbb{N}$ such that $n\ge2m\ge4$, letting $g$ be a conformally Euclidean metric on $\mathbb{R}^n$, we consider the question of positivity of the lower-order $Q$-curvatures $Q_g^{(2k)}$ for $k\in\left\{1,\dotsc,m-1\right\}$ when $Q_g^{(2m)}$ is assumed to be nonnegative and not identically zero. We assume moreover that the scalar curvature of the metric $g$ is nonnegative near infinity if $n=2m$ or that $Q_g^{(2m)}$ satisfies a slow decay barrier condition near infinity if $n>2m$. Positive results for this question have been obtained by Gursky and Malchiodi for $m=2$ and $k=1$ in the context of closed manifolds with nonnegative scalar curvature and by Li and Xu and Li, Wei, and Xu for $m\ge2$ and $k\in\left\{1,\min(m-1,2)\right\}$ in the context of conformally Euclidean metrics on $\mathbb{R}^n$. These results hold for all $n\ge2m$. Considering the case where $m\ge4$ and $k=3$, we obtain a positive result for this question when $n\in\left\{2m,2m+1,\dotsc,4m-6\right\}$, namely for these dimensions, we obtain that if $Q_g^{(2m)}\ge0$ and $Q_g^{(2m)}\not\equiv0$ in $\mathbb{R}^n$, then $Q_g^{(6)}>0$. On the other hand, in surprising contrast with the results of Gursky and Malchiodi, Li and Xu, and Li, Wei, and Xu, we find that the answer to this question is negative when $k=3$ and $n\ge N_m$ for some $N_m\in\mathbb{R}$. In this case, we are able to construct examples of conformally Euclidean metrics such that $Q_g^{(2m)}$ is positive everywhere, but $Q_g^{(6)}$ is negative at some point. By stereographic projection, our examples extend to metrics conformal to the standard metric on $\mathbb{S}^n$.

5. 一类含偶非线性项的非局部椭圆方程的多重性定理

A multiplicity theorem for a class of nonlocal elliptic equations involving even nonlinearities

作者: Biagio Ricceri

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17969v1

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

摘要: In this paper, on a bounded domain $Ω\subset {\bf R}^n$, we consider a nonlocal problem of the type

where $Q:H^1_0(Ω)\times {\bf R}\to {\bf R}$, proving, under suitable assumptions, the existence of at least three weak solutions for each $λ$ running in a suitable interval.

6. 二维空间形式中某些凸区域上第二Neumann特征函数的临界点

Critical points of second Neumann eigenfunctions on some convex domains in two-dimensional space forms

作者: Haiyun Deng, Xuyong Jiang, Xiaoping Yang

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17882v1

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

摘要: In this paper, we investigate critical points of second Neumann eigenfunctions on convex domains in the two-dimensional space forms. We approach this problem from three complementary perspectives: spectral and geometric conditions; explicit quantitative location restrictions; the hot spots constant. Precisely, for the spectral and geometric conditions, we prove that if a convex domain $Ω$ is contained in the hemisphere and satisfies $μ_{2}(Ω)\leq 2$, then its second Neumann eigenfunction has no interior critical points. Beyond this, we establish a unified diameter-based criterion $μ_2(Ω)D^2\leq j_{1,1}^2$ ensuring the absence of interior critical points in $\mathbb{S}^{2}$ and $\mathbb{H}^{2}$. Moreover, when interior critical points may exist, we derive explicit quantitative location restrictions in terms of the domain’s diameter in $\mathbb{S}^{2}$ and $\mathbb{H}^{2}$. Finally, we study the hot spots constant $\mathfrak{C}(Ω)$ on convex domains using purely analytical methods. We refine the known Euclidean upper bound of $\mathfrak{C}(Ω)$ to $2.4828,$ and obtain the corresponding hot spots constants for convex domains in non-Euclidean space forms for the first time. Our proofs combine the properties of Bessel and Legendre functions, estimation of eigenvalues and Green formulas. Our results quantitatively measure ``how wrong’’ the hot spots conjecture can be.

7. 各向异性极小曲面方程与Dirichlet边界条件

Anisotropic minimal surface equation with Dirichlet boundary condition

作者: Lu Chen, Jiali Lan, Haolin Liu

类别: math.AP

发布日期: 2026-07-20

arXiv ID: 2607.17853v1

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

摘要: This paper investigates the Dirichlet problem for the anisotropic minimal surface equation in a bounded domain. Under the natural assumption of non-negative boundary anisotropic mean curvature, we establish the unique solvability of the Dirichlet problem for continuous boundary data. To achieve this, an essential a priori gradient estimate is established, which also allows us to prove a weak version of Bernstein’s theorem for entire solutions with linear growth. Moreover, using the direct method in the calculus of variations, we prove the existence and local Lipschitz regularity of generalized minimizers in $BV(Ω)$ with $L^1(\partialΩ)$ boundary data. We also find that this variational formulation naturally yields a Neumann boundary condition, geometrically explaining why the Neumann problem requires no boundary curvature constraints.

8. 图中区域的刘维尔定理

A Liouville Theorem for Domains in Graphs

作者: Philipp Hake, Matthias Keller, Felix Pogorzelski, Marcel Schmidt

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

发布日期: 2026-07-20

arXiv ID: 2607.17814v1

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

摘要: We show a Liouville theorem for domains in graphs with Dirichlet boundary conditions. More specifically, we characterize the non-existence of non-zero bounded harmonic functions. Since Dirichlet boundary conditions give rise to Laplacians with a positive killing term, we can first characterize the validity of the Liouville theorem by the fact that the Green operator applied to the killing terms is equal to 1, or in other words, that the constant function $1$ is a potential. Secondly, we derive a characterization in terms of stochastic completeness at infinity and and total loss of heat. Thirdly, we give a characterization in terms of a Green formula for superharmonic potentials. Finally, we investigate the Liouville property in terms of recurrence and transience of the graph without killing term. As an application, we consider subsets of the Euclidean space such as cones and percolation clusters, and weakly spherically symmetric graphs.

9. 奇异 Born–Infeld 型泛函的极小元与弱解

Minimizers and Weak Solutions for Singular Born–Infeld Type Functionals

作者: Tengyang Liu, Ruifeng Zhang

类别: math.AP, math-ph

发布日期: 2026-07-20

arXiv ID: 2607.17794v1

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

摘要: We investigate the relation between minimizers and weak solutions for a class of singular functionals arising from Born–Infeld type theories $\mathcal{L}(s)$. In the setting of an electrostatic field $s=\frac{1}{2}|\nablaφ|^2$, $\mathcal{L}(s)$ satisfies $\lim_{s\to(1/2)^-}\mathcal{L}(s)=+\infty$, which naturally enforces the finite gradient bound $|\nablaφ|\le 1$, also called the truncation threshold. For a prescribed extended charge density $ρ$, we consider the relation between the weak solution of the system \begin{equation} \begin{cases} -{\rm div}\left(b\left(\frac12|\nablaφ|^2\right)\nablaφ\right)=ρ,& \text{in }\mathbb{R}^N,\ b(s)=\mathcal{L}’(s),\quad\lim_{s\to\frac12^-} b(s)=+\infty,\ \lim_{|x|\to\infty}φ(x)=0 \end{cases} \end{equation} and the minimizer $φ_0$ of the singular functional. We propose a monotonic approximation method to handle the intrinsic singularities of $\mathcal{L}(s)$. We prove that the gradient of the minimizer never touches the singular boundary $|\nablaφ|^2=1$; this structural result yields a key integrability property, the existence and uniqueness of the minimizer, and the corresponding variational inequality. Under the additional assumption that $ρ$ is radially distributed, we show that the minimizer is the unique weak solution. Furthermore, we establish the $C^1$ and $C^2$ regularity of the minimizer under suitable integrability conditions on $ρ$, and provide a uniform estimate for the strict spacelikeness condition $|\nablaφ_0|\le 1-ε$, where the parameter $ε>0$ is explicitly characterized in terms of the spatial dimension, the spatial region, and $ρ$. Our results extend the classical Born–Infeld theory to a general class of singular Born–Infeld type theories, thereby providing a unified framework for the variational analysis and regularity of such singular functionals and systems.

10. Banach 和 Hilbert 空间上脉冲线性时不变系统的统一指数稳定性分析：非强制性与强制性稳定性条件

Uniform Exponential Stability Analysis of Impulsive Linear Time-Invariant Systems on Banach and Hilbert Spaces: Non-Coercive and Coercive Stability Conditions

作者: Corentin Briat, Francesco Ferrante, Christophe Prieur

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

发布日期: 2026-07-19

arXiv ID: 2607.17260v1

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

摘要: We consider the uniform exponential stability analysis of infinite-dimensional impulsive systems defined on a Banach or Hilbert space, whose flow is governed by a fixed $C_0$-semigroup generator and whose jumps occur at a prescribed time sequence. While the flow and jump maps are themselves time-invariant, the time-triggered impulses render the propagator a genuinely time-varying evolution family, which is the source of the analysis difficulty addressed here. We combine ideas from hybrid systems theory and infinite-dimensional systems to produce operator-based stability conditions, which can be analytically or numerically checked via convex programming. Necessary and sufficient conditions for the uniform exponential stability of impulsive systems on Banach spaces are obtained in the context of a fixed impulse-times sequence but also of arbitrary, constant, minimum, and range dwell-times using both non-coercive and coercive Lyapunov functionals. Some of those results are then adapted to systems on a Hilbert space and quadratic Lyapunov functionals. As an application, linear switched systems are shown to be an exact special case: reformulated as impulsive systems with unit-norm selector jumps, they inherit non-coercive and clock-dependent dwell-time stability conditions on both Banach and Hilbert spaces. Theoretical and numerical examples are given for illustration, notably on the sampled-data control of time-delay systems.

11. 平行形变下平衡包络的Hadamard公式

A Hadamard Formula for Equilibrium Envelopes under Parallel Deformation

作者: Ziyu Li

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

发布日期: 2026-07-19

arXiv ID: 2607.17187v1

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

摘要: Let $(X,ω_0)$ be a compact Kähler manifold, and let $U\Subset X$ have $C^{3,1}$ uniformly strongly pseudoconvex boundary. For the equilibrium envelope $u_0$ associated with $X\setminus U$, the normalized Monge–Ampère measure vanishes on $U$ and equals the background measure $V^{-1}ω_0^n$ on $X\setminus\overline U$; its remaining component is a singular measure supported on $\partial U$, where $V=\int_Xω_0^n$. We identify this boundary component as the negative outward Anzellotti trace of a divergence-measure flux current. We then prove a one-sided Hadamard formula for the normalized Monge–Ampère energy along the outward parallel family $U_\varepsilon=\{ρ<\varepsilon\}$. The nonlinear telescoping identity gives the mixed Bedford–Taylor boundary traces that sum up to the trace of the current

These results provide a local weak formulation of boundary flux and normal variation for regular interface problems related to Darcy/Hele–Shaw type problems and Monge–Ampère growth.

12. 一种应力-应变对偶性的切片方法

A slicing approach to stress-strain duality

作者: Virginia De Cicco, Giovanni Scilla

类别: math.FA, math.AP

发布日期: 2026-07-19

arXiv ID: 2607.17183v1

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

摘要: The classical Kohn-Temam stress-strain pairing $({\bf A}:E{\bf u})$ for symmetric tensors ${\bf A}$ and ${\bf u}\in BD$ is typically formulated under summability assumptions on the divergence of ${\bf A}$. This excludes stress fields whose divergence has singular surface contributions, as occurs at cracks and material interfaces in continuum mechanics. We define and study stress-strain pairings for bounded symmetric divergence-measure tensor fields. For general ${\bf u}\in BD$, we introduce a slicing pairing $(({\bf A}:E{\bf u}))_Ξ$ for tensor fields satisfying a directional $BV$-type condition with respect to a finite frame $Ξ$. The definition is based on a one-dimensional disintegration strategy, and despite this construction, the new pairing enjoys analogous properties of the usual pairing $({\bf A}:E{\bf u})$, such as the absolutely continuity with respect to $|E{\bf u}|$ and the Gauss-Green formulas. We also identify several situations in which the pairing is independent of the choice of frame $Ξ$, including the relevant case in which the stress field ${\bf A}$ belongs to $BV$. While a distributional stress-strain pairing can be defined naturally for bounded $BD$ functions, it cannot be extended to the unbounded setting, since the truncation techniques available in $BV$ fail in $BD$. The slicing pairing is consistent with the distributional one whenever the latter is defined, while being more general even for bounded ${\bf u}$. Indeed, its existence does not require the compatibility condition $|{\rm Div}\,{\bf A}|(S_{\bf u}\setminus J_{\bf u})=0$ which is necessary for the distributional definition. This allows the treatment of stress fields interacting with diffuse micro-cracking.

13. 几何图上的统一Sobolev不等式

Uniform Sobolev inequalities on geometric graphs

作者: Samuel Mercer, Yves van Gennip

类别: math.AP, math.CO

发布日期: 2026-07-18

arXiv ID: 2607.16948v1

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

摘要: There is significant interest in the study of calculus on graphs, especially regarding the use of gradient-based methods for applications in data driven problems such as classification, clustering and regularisation for inverse problems. Geometric graphs, whose vertices are take from from a Euclidean domain and whose edge structure is determined by the distance between the nodes in the domain, have been central in theoretical studies. Typical approaches for analysis, such as studying consistency and the existence of continuum limits, rely on $Γ$-convergence. This technique has some limitations, as it requires the typical length scale which determines the connectivity structure of the graph to be much larger than the scales frequently used for applications. Moreover, it may fail to provide quantitative results. This paper provides necessary and sufficient conditions on the asymptotic behaviour of this length scale for the existence of a uniform collection of Sobolev inequalities on a sequence of geometric graphs. Furthermore, these inequalities hold when the length scales are much smaller than what is typically assumed for $Γ$-convergence results and within the range of what is used for data-driven problems. The Sobolev inequalities provide a quantitative estimate on the $L^q$-regularisation effect of discrete gradients.

14. 来自椭圆等位面的几何梯度流：法向分解与反射动力学

Geometric Gradient Flows from Elliptic Level Sets: Normal Decomposition and Reflection Dynamics

作者: Mohammed Barkatou, Mohamed El Morsalani

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16912v1

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

摘要: We investigate the asymptotic geometry of shifting superlevel sets $Ω_t = \{x \in Ω: u(x) > t\}$ generated by solutions to the elliptic Dirichlet problem $-Δu = f$ in $Ω$, where the non-negative source $f \not\equiv 0$ is compactly supported within a strictly convex inner core $C \subset Ω$. Under a quantitative radial monotonicity condition, each boundary $\partialΩ_t$ is a smooth normal graph over $\partial C$ characterized by a thickness function $d_t \in C^{1,α}(\partial C)$ tracking $d_0$ as $t \to 0$.A central contribution is a rigorous decomposition of the inward unit normal field along the level surfaces: $\mathbf{n}_{Ω_t} = ν- \nabla_{\partial C} d_t + \mathcal{G} + \mathcal{P}$, where $ν$ is the static radial normal, $-\nabla_{\partial C} d_t$ is the kinematic driving vector, and $\mathcal{G}, \mathcal{P}$ are curvature and PDE Hessian remainder operators.In a thin-shell configuration $(\Vert{}d_0\Vert{}_{C^1} \ll 1)$, we formalize a discrete specular point-reflection mapping $F_n$ on $\partial C$. We prove the tangential displacement satisfies $F_n(p) - p = -2d_n(p)\nabla_{\partial C}d_n(p) + R_n(p)$, with quadratic control $\Vert{}R_n\Vert{}_{L^\infty} \le C\Vert{}d_n\Vert{}_{C^1}^2$. Using the energy $\mathcal{E}(t) = \int_{\partial C} d_t^2 \, d\mathcal{H}^{N-1}$, we show these orbits approximate a continuous gradient flow driven by $+\nabla_{\partial C} d_{\tilde{t}}(p)$ to first order. Finite element computations (FEniCS) validate these convergence rates.

15. 理想 Gårding 多项式

Ideal Gårding polynomials

作者: Hao Fang, Biao Ma

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

发布日期: 2026-07-18

arXiv ID: 2607.16832v1

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

摘要: We introduce ideal Gårding polynomials, a convexity-enhanced subclass of Gårding polynomials whose Gårding components are recursively convex under partial differentiation. This class strictly contains real stable polynomials and, after translation and homogenization, lies in the Lorentzian class. Our main result is that ideal Gårding polynomials still admit a robust structure theory despite this additional convexity: they are preserved under polarization, satisfy natural closure properties, and support a linear preserver theory. A key contribution of this paper is a universal model for univariate Gårding polynomials, described by monotone root sequences and equivalently by volume polynomials of Pitman–Stanley polytopes. We establish quotient concavity, and Newton–Maclaurin type inequalities, which leads to the polarization theorem, and suggests further connections with convex geometry and Lorentzian polynomials.

16. 对数型双相问题的 Hölder 正则性

Hölder regularity for logarithmic double phase problems

作者: Ky Ho, Yun-Ho Kim, Patrick Winkert

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16749v1

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

摘要: We investigate boundedness and regularity properties of weak solutions to a class of generalized logarithmic double phase equations with variable exponents. The considered operators arise from Musielak-Orlicz type energies of logarithmic double phase type and exhibit nonstandard growth features. Under general structural assumptions, we derive a priori boundedness estimates in the subcritical setting and establish boundedness of weak solutions also in the presence of critical growth terms. In addition, we prove global Hölder continuity up to the boundary by means of the De Giorgi iteration scheme, localization arguments, and the frozen functional technique. The obtained results extend several existing regularity results for double phase and related nonstandard growth problems.

17. 三维与四维带边界的Yamabe方程解的孤立奇点

Isolated Singularities of Solutions to the Yamabe Equation with boundary in Dimension 3 and 4

作者: Yuxuan Liao, Yuexiao Ma

类别: math.AP

发布日期: 2026-07-18

arXiv ID: 2607.16650v1

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

摘要: This paper studies the asymptotic behavior of positive solutions to the boundary Yamabe equation near an isolated singularity when the metric is not conformally flat. In dimensions 3 and 4, we establishes the sharp upper bound and, in the non-removable case, the matching lower bound, which also gives a necessary and sufficient condition for removability. Moreover, every solution with a non-removable singularity is shown to be asymptotically cylindrically symmetric, without relying on a global classification of Fowler-type singular solutions. These results aslo extend the flat half-space theory of Caffarelli-Jin-Sire-Xiong (2014) to non-flat boundary geometries and provide boundary analogues of the interior theories developed by Marques (2008) and Xiong-Zhang (2022).

18. 量子图上热半群的Benjamini-Schramm极限

Benjamini-Schramm limit of the heat semigroup on quantum graphs

作者: Mihály Kovács, Eszter Sikolya

类别: math.AP, math.CO

发布日期: 2026-07-17

arXiv ID: 2607.16489v1

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

摘要: We study the behaviour of heat semigroups on quantum graphs under Benjamini-Schramm convergence. For quantum graphs with uniformly bounded geometry, equipped with continuity and Kirchhoff vertex conditions, we show that the heat semigroup, transported to a common Hilbert space by a canonical breadth-first identification of the edges, depends continuously on the underlying rooted quantum graph with respect to a local Benjamini-Schramm-type metric. As a consequence, root-averaged pairings of the semigroup converge along Benjamini-Schramm convergent sequences of finite quantum graphs. Combining this with a Trotter-Kato-type approximation of the semigroup by semigroups on metric balls, we obtain a double-limit theorem interchanging the truncation radius and the graph limit.

19. 分数阶 $p$-拉普拉斯算子的伯努利问题

Bernoulli problem for the fractional $p$-Laplacian

作者: Rafael Costa, Diego Marcon

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

发布日期: 2026-07-17

arXiv ID: 2607.16444v1

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

摘要: We study regularity properties of minimizers for the one-phase Alt–Caffarelli problem associated with the fractional $p$-Laplacian in the range $p\geq2$. We consider minimizers of the fractional $p$-energy penalized by the measure of the positivity set, with prescribed nonnegative exterior datum. We prove existence of minimizers and derive their basic variational properties: minimizers are nonnegative and are weak subsolutions of the homogeneous fractional $p$-Laplace equation. The main regularity argument combines fractional $p$-harmonic replacements, energy-gap estimates, nonlocal tail bounds, and a Campanato-type iteration. This yields local Hölder continuity of minimizers and implies that they solve the homogeneous equation in their positivity set. Finally, we prove the optimal free boundary growth estimate, showing that the sharp order known in the linear fractional Bernoulli problem persists in the fractional $p$-Laplacian setting.