生成时间: 2026年05月14日 10:21:50

运行状态: 完成

计划处理论文数量: 49

已完成论文数量: 49

未完成论文数量: 0

重点关注论文数量: 16

了解领域论文数量: 17

不相关论文数量: 16

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

1. 三维Navier-Stokes方程一个初等模型中的有限时间爆破

Finite-time blow-up in an elementary model of the 3D Navier-Stokes equations

作者: Stan Palasek

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13827v1

摘要翻译: 我们证明了一个简单但符合物理实际的三维Navier-Stokes方程壳模型中存在有限时间爆破，该模型配备“光滑”（即在频率空间中快速衰减）的初始数据与外力项。此前研究的模型要么展现出一种能正则化三维粘性动力学的湍流级串，要么依赖于在真实Euler非线性中并不透明实现的高度人工相互作用。我们还处理了无粘、无外力情形，并在能量水平略上方得到了奇点形成结果。最后，我们讨论了将该二进制模型（dyadic model）的行为嵌入完整Euler方程与Navier-Stokes方程之中的前景。

Comment: 20 pages, 2 figures

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

1. 研究对象和背景

本文研究的对象是三维不可压Navier-Stokes方程的一个壳模型（shell model），即Obukhov型 dyadic 模型，其核心方程为

其中 $ X_k(t) \ge 0 $ 表示速度场在第 $ k $ 个频段（中心频率尺度为 $ N_k $）上的 $ L^2 $-模长近似，$ \nu \ge 0 $ 为运动粘性系数，$ \alpha \in [1,5/2] $ 是表征湍流间歇性的物理参数，$ f_k(t) $ 是外力在该频段的投影。频尺度序列 $ (N_k)_{k\ge0} $ 不再取经典几何序列 $ N_k = \lambda^k $，而是采用超指数增长形式

这是本文建模的关键设定。该模型源于对全PDE解的频域分解：设 $ u = \sum_{k\ge0} u_k $，其中 $ \operatorname{supp} \widehat{u_k} \subset \{ \xi \in \mathbb{Z}^3 : |\xi| \sim N_k \} $，并仅保留两类关键非线性相互作用——低频驱动项 $ u_k \cdot \nabla u_{k-1} $（对应 $ N_{k-1}^\alpha X_{k-1} X_k $）与高频自相互作用 $ u_{k+1} \cdot \nabla u_{k+1} $ 投影至第 $ k $ 频段的贡献（对应 $ -N_k^\alpha X_{k+1}^2 $）。后者经Littlewood-Paley投影及支持集体积比 $ \sim N_k^{-2(\alpha-1)} / N_{k+1}^{-2(\alpha-1)} $ 修正后，结合散度结构导出精确的幂律衰减因子。模型中函数空间采用Besov型加权 $ \ell^\infty $ 空间：

此空间对应于PDE中 $ B^{s+(2/p-1)(\alpha-1)}_{p,\infty}(\mathbb{T}^3) $，尤其 $ C^\infty $ 等价于光滑函数类。问题的重要性在于：三维Navier-Stokes全局正则性这一千禧年难题，其核心障碍是无法排除有限时间奇点；而直接处理PDE过于困难，故学界长期依赖简化模型探路。已有模型存在根本局限：Katz–Pavlović（KP）模型虽具物理直观性，但其能量级联（Kolmogorov-5/3律）具有内在正则化效应，当 $ \alpha < 3 $ 时导致全局适定性；标准Obukhov模型在 $ N_k = \lambda^k $ 下因高-高-低项主导而全局正则；Tao模型虽实现粘性爆破，但依赖四模/转子门等高度人工结构，与真实Euler非线性缺乏透明对应。因此，本文试图填补的关键空白是：构造一个既保持Obukhov型非线性本质、又具备有限时间爆破能力的“现实”壳模型，且其爆破机制可望通过嵌入技术（embedding）推广至全PDE系统。该空白的核心在于突破传统频尺度限制——必须让 $ N_k $ 增长快于指数，以压制高-高-低项的累积破坏效应，并允许 $ C^\infty $ 初值在爆破时刻退出所有 $ C^s $（$ s>0 $）空间，从而模拟PDE中 $ B^0_{2+\epsilon,\infty} $ 或 $ B^\epsilon_{2,\infty} $ 类型的Type II爆破。

2. 主要定理或主要结果

本文在超指数频尺度 $ N_k = N_0^{b^k} $（$ N_0 > 1, b > 1 $）下，严格证明了Obukhov模型的两类有限时间爆破：（1）粘性情形（Theorem 1.3）：对任意 $ \nu > 0 $、$ \alpha > 2 $，存在正的初值 $ X_0 \in C^\infty $ 及光滑外力 $ f \in C^\infty_t([0,\infty); C^\infty) $，使得解 $ X(t) $ 的最大存在时间 $ T^* < \infty $，且对任意 $ s > 0 $，有 $ \|X(t)\|_s \to \infty $ 当 $ t \to T^* $；进一步，该外力满足 $ \|f(t)\|_s \to 0 $（$ \forall s > 0 $）当 $ t \to T^* $，即其在爆破时刻完全退耦。（2）无粘情形（Theorem 1.8）：对 $ \nu = 0 $、任意 $ \alpha \ge 1 $，存在正的初值 $ X_0 \in C^\infty $（$ f \equiv 0 $），使得解发生有限时间爆破，且同样对任意 $ s > 0 $，$ \|X(t)\|_s \to \infty $ 当 $ t \to T^* $。两个定理的关键技术条件为：爆破剖面渐近于 $ X_k(T^*) \sim N_k^{\beta - \alpha} $，其中 $ \beta > 2 $，故 $ C^s $-范数爆破当且仅当 $ s > 0 $；该剖面对应PDE中略高于 $ L^2 $ 能量空间的临界性（Remark 1.7）。此外，定理1.3中 $ \alpha > 2 $ 是尖锐的：当 $ \alpha \le 2 $ 时，系统能量临界或超临界，已知全局适定性成立；而 $ \alpha \in (2, 5/2] $ 正是三维Navier-Stokes相关间歇性参数的物理窗口。证明中隐含的尺度关系为：第 $ k $ 模的激活时间窗 $ I_k $ 满足 $ |I_k| \sim (N_{k-1}^\alpha A_{k-1})^{-1} $，其中 $ A_{k-1} \sim N_{k-1}^{\beta - \alpha} $ 为前一模的目标幅值；为保证级联可控，必须满足 $ |I_k| \ll (N_{k-2}^\alpha A_{k-2})^{-1} $，这要求 $ N_k $ 的双指数增长以产生足够大的分离比 $ N_{k-1}^\alpha A_{k-1} / N_{k-2}^\alpha A_{k-2} \sim N_k^c $（$ c > 0 $），从而允许 $ C^\infty $ 初值构造。

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

证明的核心障碍在于：如何在保留Obukhov模型物理结构的前提下，同时压制两类破坏性效应——无粘情形下的高-高-低项 $ -N_k^\alpha X_{k+1}^2 $ 的累积负反馈，以及粘性情形下耗散项 $ -\nu N_k^2 X_k $ 对级联增长的阻滞。作者的真正创新点在于将频尺度 $ N_k $ 的增长速率本身转化为控制参数，并据此重构整个分析框架。首先，引入关键的重标度变量 $ x_k := N_k^\alpha X_k $，将原系统化为

其中 $ \delta_k $ 因 $ N_{k+1} \gg N_k $ 而成为极小无量纲参数（例如双指数下 $ \delta_k \lesssim e^{-c N_k} $），这使高-高-低项在动力学中可被一致控制。其次，采用反向构造法（backward construction）：不从初值出发正向积分，而是设定爆破时刻 $ t=0 $ 的奇异剖面 $ x_k(0) = A_k $（满足 $ A_k \sim N_k^\beta $），并考虑截断系统 $ \{x_k\}_{k=0}^K $ 在 $ [-T,0] $ 上的反向演化。利用大参数 $ N_0 \gg 1 $ 和 $ b > 1 $ 构造一个不变的“捕获区域”（trapping region）：对每个 $ k $，定义 $ x_k $ 的上下界为 $ \frac{1}{2} A_k \le x_k(t) \le 2 A_k $（$ t \in [-T_k, 0] $），其中 $ T_k $ 由尺度关系 $ T_k \sim (A_{k-1})^{-1} $ 确定。通过精细的Gronwall型估计，证明反向解在此区域内满足一致先验界，且 $ \delta_k $ 的指数小性确保高阶项扰动可忽略。随后令 $ K \to \infty $ 取极限，得到全局反向解，并验证其在 $ t = -T $ 处光滑（$ C^\infty $）。最后，将此反向解在 $ t = -T $ 处的值作为正向初值，由唯一性即得正向爆破解。对于粘性情形，额外引入外力 $ f_k $ 的设计是另一关键技术：其作用是在每个模 $ k $ 的激活时间窗 $ I_k $ 之前，精确抵消耗散项 $ -\nu N_k^2 x_k $，即令 $ g_k(t) \approx \nu N_k^2 x_k(t) $ 当 $ t \in [-T, -T_k] $，而 $ g_k(t) \equiv 0 $ 当 $ t \in [-T_k/2, 0] $；由于初值 $ x_k(-T) $ 本身随 $ k $ 指数衰减（因 $ x \in C^\infty $），该外力可被构造为 $ C^\infty_t C^\infty $ 类且支集远离 $ t=0 $。整个框架摒弃了传统能量方法，转而依赖尺度分离引理（scale-separation lemma）、反向捕获区域论证（backward trapping argument）以及超指数频尺度诱导的几何级数控制（geometric control from double-exponential separation）。

4. 与之前工作的比较

本文工作与既有文献的差异体现在四个维度：（1）模型真实性：相较于Tao [20] 模型中人为引入的“转子门”（rotor gates）和四模耦合，本文模型严格继承Obukhov原始物理动机——仅保留 $ u_k \cdot \nabla u_{k-1} $ 与 $ u_{k+1} \cdot \nabla u_{k+1} $ 两类真实存在于Euler非线性中的相互作用，且其系数 $ N_{k-1}^\alpha $、$ N_k^\alpha $ 直接由频域投影和间歇性体积比导出，无任何自由参数。因此，其嵌入至全PDE的可行性远高于Tao模型。（2）假设放宽性：KP模型与标准Obukhov模型均要求 $ N_k = \lambda^k $，而本文首次系统研究并利用 $ N_k $ 超指数增长（如 $ N_k = N_0^{b^k} $）作为正则性破缺的主动工具。此前Looi [16] 已指出，在几何尺度下粘性Obukhov模型全局正则，而本文证明只要 $ N_k $ 增长足够快，同一模型即可爆破，从而将“全局适定性”的适用范围从 $ N_k = \lambda^k $ 精确限定为“非超指数增长”。（3）时间尺度与适用范围：KP模型的爆破需初始数据非光滑（如 $ X_0 \notin \ell^2 $），而本文爆破发生在 $ C^\infty $ 初值下，对应PDE中光滑解的爆破，更具物理意义；其爆破类型为Type II（退出所有 $ C^s $，$ s>0 $），位于能量空间 $ C^0 $（等价于 $ \ell^2 $）正上方，比KP模型的 $ \ell^1 $-型爆破更接近Navier-Stokes临界空间 $ B^{1/3}_{3,\infty} $ 的尺度。（4）代价与局限：本文的代价是牺牲了模型的普适湍流描述能力——超指数尺度 $ N_k $ 不再对应各向同性湍流，而专为构造奇点设计；且目前仅证明存在性，未给出爆破时间的显式估计或稳定性分析（Remark 1.10指出爆破在 $ C^s $ 中不稳定）。与Jeong-Li [11] 的混合KP-Obukhov模型相比，本文纯Obukhov结构避免了KP主导的级联机制，揭示了Obukhov型相互作用自身蕴含的爆破潜力。从提供内容中无法确认本文爆破解是否具有自相似结构，亦未显式展开与Onsager猜想或Buckmaster-Vicol非唯一性工作的定量比较。

2. Kuramoto 均场博弈中同步稳态平衡的唯一性

Uniqueness of synchronized stationary equilibria in the Kuramoto mean field game

作者: Sebastian Munoz

类别: math.AP, math.OC

发布日期: 2026-05-13

arXiv ID: 2605.13783v1

摘要翻译: 稳态 Kuramoto 均场博弈刻画了一类相位振子群体，当相互作用强度超过某一临界值时，该群体形成同步化的 Nash 平衡。我们证明：从均匀态在分岔点处产生的同步分支是唯一的光滑 Nash 平衡族；即对每个超临界相互作用强度，同步平衡在环面旋转意义下唯一，并且当相互作用参数递减趋近于临界阈值时，该平衡光滑地收敛于均匀分布。这两个结论均源于我们的主要技术结果：标量自洽映射 $\Phi(\rho)$ 是严格凹的，从而证实了 Carmona、Cormier 与 Soner 提出的一个猜想。证明将自洽映射的二阶导数分解为平衡态的两个符号不定矩——一个三次矩与一个梯度矩——并通过对值函数的尖锐形状估计、一种逐点几何平均单调性（该单调性通过余弦偏度不等式确定三次矩的符号），以及结合相关不等式的反射论证（用于控制梯度矩的符号）来分别控制这两项的符号。

Comment: 16 pages

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

1. 研究对象和背景

本文研究的是Kuramoto均场博弈（Kuramoto mean field game）的稳态Nash均衡结构，其核心模型为一个受控扩散过程在环面 $\mathbb{T} = \mathbb{R}/2\pi\mathbb{Z}$ 上的动力学系统。代表性智能体的状态 $X_t$ 满足

其中控制 $\alpha_t$ 是智能体为最小化折扣成本而自主选择的漂移项，$\sigma > 0$ 为噪声强度，$B_t$ 为标准布朗运动。目标泛函为

其中 $\beta > 0$ 为贴现率，$\kappa > 0$ 为耦合强度，$\mu_t \in \mathcal{P}(\mathbb{T})$ 是群体在时刻 $t$ 的相位分布，而Kuramoto耦合函数定义为

其中 $\widehat{\mu}(k)$ 为傅里叶系数；特别地，当 $\mu$ 关于原点对称时（如后文旋转归一化后），$c(\theta, \mu) = 1 - \mu(\cos)\cos\theta$，即仅依赖于一阶矩 $\gamma := \mu(\cos)$。该耦合机制体现相位对齐倾向：$c$ 在 $\theta$ 与群体平均相位一致时最小。

在稳态设定下，$\mu_t \equiv \mu$，价值函数 $\Phi(t,x) = \phi(x)$ 满足定常Hamilton–Jacobi–Bellman（HJB）方程

最优反馈律为 $\alpha^*_t = -\phi_x(X_t)$，对应扩散过程 $dX_t = -\phi_x(X_t)\,dt + \sigma\,dB_t$。稳态Nash均衡要求 $\mu$ 是该扩散过程的不变测度，其密度满足Gibbs形式

由于 $c$ 具有旋转不变性，所有旋转等价类构成对称轨道；通过旋转可设 $\mu(\sin) = 0$ 且 $\mu(\cos) \geq 0$，此时 $\mu$ 完全由标量序参量 $\gamma = \mu(\cos)$ 决定。将HJB方程无量纲化：令 $\lambda = 2\beta/\sigma^2$，$\zeta = 4\gamma/\sigma^4$，$v_\zeta = 2\phi/\sigma^2$，则方程化为

且 $\mu_\zeta(dx) = f_\zeta(x)\,dx$，其中 $f_\zeta(x) = e^{-v_\zeta(x)} / Z_\zeta$，$Z_\zeta = \int_{\mathbb{T}} e^{-v_\zeta}$. 此时自洽映射定义为 $A(\zeta) := \mu_\zeta(\cos) = \int_{\mathbb{T}} \cos x\, f_\zeta(x)\,dx$，原始自洽映射为 $F_\kappa(\gamma) = \kappa A(4\gamma/\sigma^4)$。稳态均衡即解方程 $\gamma = F_\kappa(\gamma)$，等价于寻找 $A(\zeta)$ 的不动点。已知 $\gamma = 0$（即均匀分布 $\mu \equiv 1/(2\pi)$）恒为解，对应非同步相；而当 $\kappa > \kappa_c := (\beta\sigma^2 + \sigma^4)/2$ 时，Carmona–Cormier–Soner [CCS23] 已证存在非平凡同步解。但该分支是否唯一、是否光滑连续地从临界点 $\kappa = \kappa_c$ 分岔而出，仍是未解问题。尤其关键的是，[CCS23] 基于数值实验提出核心猜想：$F_\kappa$ 在 $[0,\infty)$ 上严格凹，此性质若成立，则直接推出同步解在旋转意义下的唯一性及临界连续性。本文旨在严格证明该凹性猜想，并由此确立同步稳态均衡的整体结构稳定性。

2. 主要定理或主要结果

定理1.1给出全文核心结论：对任意 $\beta, \sigma, \kappa > 0$，(i) 标量自洽映射 $F_\kappa$ 在 $[0,\infty)$ 上严格凹；(ii) 当 $\kappa \leq \kappa_c$ 时，不存在非平凡同步稳态均衡；当 $\kappa > \kappa_c$ 时，存在唯一（模环面旋转）的同步稳态均衡 $\mu_\kappa$，满足归一化条件 $\mu_\kappa(\sin) = 0$ 且 $\mu_\kappa(\cos) > 0$；进一步，$\mu_\kappa$ 的密度作为 $\kappa$ 的函数，在区间 $(\kappa_c, \infty)$ 上 $C^\infty$-光滑地取值于 $C^\infty(\mathbb{T})$，且当 $\kappa \downarrow \kappa_c$ 时，该密度在 $C^\infty(\mathbb{T})$ 中收敛于均匀密度 $1/(2\pi)$。具体而言，无量纲化后，$F_\kappa(\gamma) = \kappa A(4\gamma/\sigma^4)$，故 $F_\kappa$ 的严格凹性等价于 $A(\zeta)$ 在 $\zeta > 0$ 上严格凹，即 $A''(\zeta) < 0$ 对所有 $\lambda > 0$ 和 $\zeta > 0$ 成立。引理4.1给出精确二阶导数分解公式：令 $w := A(\zeta)/\lambda + v_\zeta$，$z := \partial_\zeta v_{\zeta,x}$，则 $\int_{\mathbb{T}} w f_\zeta = 0$，且

该式将凹性问题转化为两个积分项的符号判定：立方矩 $\int w^3 f$ 和梯度矩 $\int w z^2 f$。命题4.4证明前者 $\geq 0$，命题4.6证明后者 $> 0$，从而 $A''(\zeta) < 0$。此外，临界阈值 $\kappa_c$ 由线性化分析确定：$F_\kappa'(0) = \kappa/\kappa_c$，故分岔发生在 $\kappa = \kappa_c$；且 $|F_\kappa(\gamma)| \leq \kappa$，保证 $g_\kappa(\gamma) := F_\kappa(\gamma) - \gamma$ 在 $\gamma \to \infty$ 时趋向 $-\infty$，结合 $g_\kappa(0) = 0$ 及 $g_\kappa'(0) > 0$（当 $\kappa > \kappa_c$），严格凹性确保 $g_\kappa$ 在 $(0,\infty)$ 上恰有一个正零点 $\gamma^*(\kappa)$，且 $\gamma^*(\kappa) \to 0$ 当 $\kappa \downarrow \kappa_c$。

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

证明的核心障碍在于：引理4.1给出的 $A''(\zeta)$ 表达式中两项——立方矩 $\int w^3 f$ 和梯度矩 $\int w z^2 f$——其被积函数均无显见符号，标准能量估计或分部积分无法直接判定其正负。作者的真正创新在于发展了一套精细的形状分析框架，将高维积分不等式降维为一维单调性与凸性问题，并引入若干原创性不等式工具。

第一，针对 $u = v_{\zeta,x}$ 在 $(0,\pi)$ 上的几何结构，引理3.1建立四重形状估计：$u > 0$、$u_{xx} \leq 0$（凹性）、$(u/\sin x)_x \geq 0$（正弦缩放单调性）、$|u_x| \leq u/\sin x$（弦界）。其证明采用“Wronskian–曲率”耦合系统方法：定义 $W = \sin x\, u_x - \cos x\, u$（编码 $u/\sin x$ 单调性）与 $K = \zeta \sin x / u - \lambda - u_x$（编码 $u_{xx} = -uK$ 的凹性），二者满足封闭的一阶ODE系统

通过最大值原理与渐近展开（在 $x=\pi$ 处用 $\xi = \pi - x$ 展开 $u = a\xi + b\xi^3 + o(\xi^3)$ 并利用 $u>0$ 推出 $a>0$），联合证明 $W \geq 0$ 与 $K \geq 0$，从而同时确立单调性与凹性。

第二，针对 $z = \partial_\zeta v_{\zeta,x}$，引理3.2引入正弦缩放商 $Z = z/\sin x$，并证明 $0 \leq Z_x \leq u Z$。该上界等价于通量 $-(f Z)_x \geq 0$，其证明依赖于引理3.1中辅助函数 $B = 1 + \lambda + u_x + u \cot x$ 的单调递减性 $B_x \leq 0$，这是处理 $Z$ 方程 $-Z_{xx} + (u - 2\cot x) Z_x + B Z = 1$ 的关键。

第三，立方矩的非负性（命题4.4）依赖于原创的“几何平均单调性”（引理4.2）：定义 $\Delta_0(x) = v_\zeta(x) - v_\zeta(0)$，$\Delta_\pi(x) = v_\zeta(\pi) - v_\zeta(x)$，则函数

在 $(0,\pi)$ 上非增。其证明将 $H_x \leq 0$ 转化为关于 $\Delta_0, \Delta_\pi$ 的积分下界/上界（式(29)），再利用 $e^{-v} Z$ 的单调性（由 $Z_x \leq u Z$ 导出）将其与初值问题 $I_0' + u I_0 = \sin x$、$I_\pi' + u I_\pi = -\sin x$ 的解联系起来，最终通过比较 $I_0, I_\pi$ 与三角函数组合 $\sin x/(u+r), \sin x/(t-u)$ 完成。随后，通过变量替换 $C(x) = \cos \theta(x)$ 将立方矩化为关于 $\cos \theta$ 的中心化三阶矩，并应用引理4.3（任意 $[0,\pi]$ 上非增密度 $\eta$ 满足 $\int_0^\pi (\cos \theta - \bar{c})^3 \eta(\theta)\,d\theta \geq 0$，其中 $\bar{c}$ 为均值），该不等式由中位数对称化技巧证明。

第四，梯度矩的正性（命题4.6）基于反射论证与相关不等式：首先，引理4.5证明 $w$ 关于 $\pi/2$ 具有“中心正性”，即对任意在 $[0,\pi/2]$ 上非减、关于 $\pi/2$ 对称的测试函数 $\psi$，有 $\int_0^\pi w \psi f \geq 0$；其次，由 $w$ 与 $Z^2$ 在 $(0,\pi)$ 上同为非减函数，结合Chebyshev型相关不等式，得出 $\int w z^2 f = \int w Z^2 \sin^2 x f > 0$。

4. 与之前工作的比较

本文与既有工作存在本质差异，主要体现在问题设定、技术路径与适用范围三个维度。

首先，与Carmona–Cormier–Soner [CCS23] 相比：[CCS23] 首次建立了Kuramoto均场博弈中同步相变的存在性，但其分析止步于存在性与大耦合极限 $\kappa \to \infty$ 下的完全同步 $\mu_\kappa \to \delta_0$，对同步分支的唯一性、局部结构及临界行为未作刻画。本文直接解决其明确提出的凹性猜想（Remark 7.4），并首次证明同步均衡在旋转意义下的全局唯一性及临界连续性 $\mu_\kappa \to 1/(2\pi)$，填补了从存在性到结构性的关键空白。数值实验虽支持凹性，但缺乏严格证明；本文的证明完全脱离数值，基于解析形状估计与不等式，具有理论完备性。

其次，与Cesaroni–Cirant [CC24] 的工作相比：[CC24] 研究的是遍历型（ergodic）Kuramoto均场博弈（即以遍历常数替代贴现率 $\beta$），其方法为变分能量估计，在强耦合 $\kappa \gg 1$ 区间证明了同步均衡的局部唯一性与Turnpike性质。然而，该方法天然受限于大参数渐近，无法触及分岔点附近的弱超临界区域 $\kappa \in (\kappa_c, \kappa_c + \varepsilon)$，而这正是相变物理发生的核心区域。本文的凹性证明对所有 $\kappa > \kappa_c$ 一致有效，特别是覆盖了 $\kappa$ 紧邻 $\kappa_c$ 的整个区间，从而完整刻画了相变本身的数学结构。文中亦指出，本文方法可经微调推广至遍历情形，且因缺少立方矩项而更为简洁，印证了当前技术对临界区域的普适优势。

第三，在方法论层面，本文与经典均场博弈文献（如Lasry–Lions [LL07]、Cardaliaguet [CD18]）形成鲜明对比：后者多依赖单调性条件（monotonicity condition）保证解的唯一性，而Kuramoto模型天然属于非单调（non-monotone）情形，传统框架失效。本文绕过单调性假设，转而挖掘HJB方程解 $v_\zeta$ 及其导数 $v_{\zeta,x}, \partial_\zeta v_{\zeta,x}$ 的深层几何结构（凹性、正弦缩放单调性、几何平均单调性），并将这些点态性质升华为全局积分不等式。这种“从点态形状到全局符号”的范式，为处理其他非单调均场博弈的唯一性问题提供了新模板。代价在于技术高度特化：所有估计均紧密耦合于Kuramoto的特定非线性（$-\zeta \cos x$ 强迫项）与环面对称性，尚未推广至一般耦合函数或高维流形。文中未显式展开与其他非单调MFG（如[YMMS12], [Cir19]）的比较，但从方法复杂度看，本文的证明远超这些工作的形式分析或数值探索，代表了该方向严格的解析突破。

3. 不可压Oldroyd-B模型在临界 $L^p$ 空间中的尖锐衰减刻画

Sharp decay characterization for the incompressible Oldroyd-B model in critical $L^p$ spaces

作者: Zhi Chen, Mingwen Fei, Lvqiao Liu, Jiahong Wu

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13598v1

摘要翻译: 本文在临界 $L^p$ 框架下建立了不可压Oldroyd-B模型时间衰减速率的尖锐刻画，覆盖了流体粘性与应力张量阻尼均缺失这一物理上相关且数学上微妙的情形。我们证明：初始数据 $(u_0, τ_0)$ 的低频部分满足一个 $L^2$-型条件，该条件对于在临界Besov空间中获得解的时间衰减的最优上下界而言，几乎是必要且充分的。一项关键贡献在于将应力张量新分解为其不可压部分与可压部分，并引入一个有效张量以处理高频速度场中正则性的损失。这是首个揭示无粘性、无阻尼不可压Oldroyd-B模型此类精确双边渐近行为的结果。

Comment: All comments are welcome

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

1. 研究对象和背景

本文研究的是无粘性、无应力阻尼的不可压Oldroyd-B模型，其控制方程为

其中 $u: \mathbb{R}^+ \times \mathbb{R}^d \to \mathbb{R}^d$ 是速度场，$\tau: \mathbb{R}^+ \times \mathbb{R}^d \to \mathbb{R}^{d\times d}_{\mathrm{sym}}$ 是非牛顿应力张量（对称矩阵），$P$ 是压力，$D(u) = \frac{1}{2}(\nabla u + (\nabla u)^\top)$ 是应变率张量，$Q(\tau, \nabla u) = \tau \Omega(u) - \Omega(u) \tau + b\big(D(u)\tau + \tau D(u)\big)$ 是双线性耦合项，$\Omega(u) = \frac{1}{2}(\nabla u - (\nabla u)^\top)$ 是旋度部分，$b \in [-1,1]$，$\nu_1, \nu_2 > 0$ 为耦合常数，$\eta > 0$ 为应力扩散系数。关键特征在于：流体粘性系数 $\nu = 0$（即速度方程无拉普拉斯耗散），且应力方程中阻尼项系数 $a = 0$（即无 $a\tau$ 项）。这使得系统既非完全抛物型，亦非纯双曲型，而是呈现典型的“部分耗散”（partially dissipative）结构：仅通过 $\nu_1 \mathrm{div}\, \tau$ 和 $-\eta \Delta \tau$ 实现间接耦合耗散。物理上，该情形对应于极高弹性、极低粘性极限下的聚合物流动；数学上，它回避了经典能量方法所依赖的显式耗散，导致全局适定性与大时间行为分析异常困难。已有工作如 [12, 27] 要求 $\nu>0, a>0$ 以保证全局解存在；[15, 16, 23] 在 $\nu=0, a>0$ 下借助阻尼项构造Lyapunov泛函；而 $a=0$ 时，线性化系统（1.3）的傅里叶符号显示其特征根为

当 $|\xi| \ll 1$ 时，$\mathrm{Re}\, \lambda_\pm \sim -\frac{\eta}{2}|\xi|^2$，表现为热核型耗散；但当 $|\xi| \gg 1$ 时，$\lambda_+ \to 0$，$\lambda_- \sim -\eta |\xi|^2$，意味着高频速度 $u$ 几乎无耗散，退化为输运主导行为，从而导致 $u$ 相对于 $\tau$ 在高频处损失整整一阶正则性。这一根本不对称性使标准临界Besov框架 $\dot{B}^{s}_{p,r}$ 难以闭合估计。现有衰减研究（如 [32, 10]）均依赖 $L^1$ 初值假设或额外小性条件，无法揭示初值低频正则性与解衰减速率之间的本质对应关系。本文旨在填补这一空白：在临界 $L^p$ 框架下，建立解的上下界衰减率与初值低频部分所属函数空间之间的几乎充要（almost necessary and sufficient）对应关系，即精确刻画何种低频正则性是获得最优衰减的最小要求。

2. 主要定理或主要结果

在 $d \geq 2$ 且 $p$ 满足 $2 \leq p \leq \min\{4, 2d/(d-2)\}$（当 $d=2$ 时 $p \neq 4$）的条件下，定义 $\sigma_0 = d/2 - 2d/p$。设初值 $(u_0, \tau_0)$ 满足小性条件

其中 $c_0 > 0$ 充分小，$\ell/h$ 分别表示低/高频截断（阈值 $k_0$ 待定）。则对Cauchy问题(1.1)的解 $(u, \tau)$，有如下精确衰减刻画：

(i) 充分性与必要性（Theorem 1.1第一部分）：对任意给定 $t_0 > 0$，解满足

当且仅当

这里 $\langle t \rangle = (1+t)$。

(ii) 双向衰减率（Theorem 1.1第二部分）：若存在 $t_0 > 0$ 使得

则必有

反之，若

则上述双向衰减成立。其中 $P = I + \nabla(-\Delta)^{-1}\mathrm{div}$ 为Leray投影，$\dot{B}^{\sigma_1}_{2,\infty}$ 是一个带有特定频率集中性质的子空间（见原文Definition后描述）。

(iii) 特例 $p=2$（Corollary 1.1）：当 $p=2$ 时，衰减率简化为

其成立当且仅当 $(u_0, \tau_0)^\ell \in \dot{B}^{\sigma_1}_{2,\infty}$，且反向蕴含亦成立。

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

证明的核心障碍源于系统(1.1)在高频处的严重正则性损失：由于 $\nu=0, a=0$，线性化系统的特征根 $\lambda_+ \to 0$，导致高频速度 $u^h$ 几乎不耗散，其行为类似于输运方程，从而 $\|u^h\|_{\dot{B}^{s}_{p,1}}$ 的估计无法从 $\|\tau^h\|_{\dot{B}^{s+1}_{p,1}}$ 直接导出，标准能量法失效。作者的真正创新点在于构建了一套全新的、针对部分耗散结构的精细分析框架，包含四个相互支撑的关键技术：

1. 应力张量的不可压/可压分解与潜在函数引入：为克服 $\tau$ 缺乏散度自由性的困难，在低频分析中，作者将 $\tau$ 分解为 $v = \Lambda^{-1}\mathrm{div}\, \tau$（不可压部分，对应速度的涡量类变量）和 $w = \Lambda^{-1}\mathrm{curl}\, \tau$（可压部分）。进而考虑三元组 $(u, v, w)$，将原系统(1.3)重写为

   $$ \begin{cases} \partial_t u - \nu_1 \Lambda P v = 0, \\ \partial_t P v - \eta \Delta P v + \frac{\nu_2}{2} \Lambda u = 0, \end{cases} $$

   此形式揭示了隐藏的耗散结构（Lemma 2.3），导出关键微分不等式

   $$ \frac{d}{dt} \left( \|u^\ell\|_{\dot{B}^{s}_{2,1}}^2 + \|v^\ell\|_{\dot{B}^{s}_{2,1}}^2 \right) + c \eta \left( \|\nabla u^\ell\|_{\dot{B}^{s}_{2,1}}^2 + \|\nabla v^\ell\|_{\dot{B}^{s}_{2,1}}^2 \right) \leq 0, $$

   从而在低频区域建立了与热方程同阶的耗散控制。

2. 有效张量（Effective Tensor）的构造：为处理高频速度的正则性损失，作者受[21]启发，引入有效张量

   $$ W = \nabla(-\Delta)^{-1} \left( \frac{1}{2} u + \mathrm{div}\, \tau \right). $$

   该构造的精妙之处在于，$W$ 的傅里叶变换 $\widehat{W}(\xi)$ 与 $\widehat{u}(\xi)$ 和 $\widehat{\tau}(\xi)$ 的低阶组合相关，且其高频部分 $\widehat{W}^h$ 满足 $|\widehat{W}^h(\xi)| \lesssim |\xi|^{-1}(|\widehat{u}^h(\xi)| + |\widehat{\tau}^h(\xi)|)$，从而将 $u^h$ 的估计转化为对 $W^h$ 和 $\tau^h$ 的联合估计（Lemma 2.4），成功绕开了直接估计 $u^h$ 所需的额外导数。

3. 点态下界估计与取消性质的挖掘：为证明衰减率的必要性（即下界），作者并未直接对 $(u, \tau)$ 估计，而是转而分析 $(u, \Lambda^{-1}\mathrm{div}\, \tau)$。在此新变量下，线性系统(2.1)展现出关键的取消性质（Lemma 2.1）：其解的傅里叶模满足

   $$ |\widehat{u}(t, \xi)| + |\widehat{\tau}(t, \xi)| \geq C e^{-\eta |\xi|^2 t} \left( |\widehat{u}_0(\xi)| + |\widehat{P\Lambda^{-1}\mathrm{div}\, \tau_0}(\xi)| \right), \quad |\xi| \leq c. $$

   此下界是整个必要性论证的基石，它严格排除了衰减快于 $e^{-\eta |\xi|^2 t}$ 的可能性，从而将衰减速率的下界精确锚定在初值低频的 $\dot{B}^{\sigma_1}_{2,\infty}$ 正则性上。

4. 非线性项 $Q(\tau, \nabla u)$ 的精细化处理：在误差系统(1.4)的估计中，对 $Q(\tau, \nabla u)$ 的低频贡献，作者采用分裂技巧

   $$ \tau^\ell_{ik} \partial_k \tilde{u}^h_j = \partial_k(\tau^\ell_{ik} \tilde{u}^h_j) - (\partial_k \tau^\ell_{ik}) \tilde{u}^h_j. $$

   第一项利用 $\dot{B}^{1}_{p,1}$ 嵌入获得导数增益；第二项则利用 $\|\nabla \tilde{\tau}^\ell\|_{\dot{B}^{d/2-1}_{2,1}}$ 的快速衰减（因 $\sigma'$ 可取极小），从而分别控制了两项的时间衰减阶，避免了标准乘积估计的失效。

4. 与之前工作的比较

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

1. 与带阻尼/粘性工作的比较：经典结果如 [12, 27]（$\nu>0, a>0$）、[15, 16, 23]（$\nu=0, a>0$）均依赖显式阻尼项 $a\tau$ 或 $\nu \Delta u$ 构造Lyapunov泛函，其衰减率由这些项主导（如 $e^{-at}$ 或 $e^{-\nu |\xi|^2 t}$）。本文首次在 $a=\nu=0$ 的极端情形下，仅依靠耦合耗散 $\nu_1 \mathrm{div}\, \tau$ 和 $-\eta \Delta \tau$，就获得了与热方程同阶的 $e^{-c|\xi|^2 t}$ 衰减，且该速率被证明是最优且不可改进的。代价是必须放弃全局 $L^2$ 能量守恒，转而采用高度异质的高低频分离框架（$\ell \dot{B}^{d/2-1}_{2,1} + h \dot{B}^{d/p+1}_{p,1}$），这是对传统临界框架的重大突破。

2. 与先前衰减工作的比较：[32] 和 [10] 均假设初值 $(u_0, \tau_0) \in L^1$，并由此推出 $\|u(t)\|_{L^2} \lesssim (1+t)^{-d/4}$。本文彻底摆脱了 $L^1$ 假设，代之以更本质、更弱的低频Besov空间条件 $(u_0, \tau_0)^\ell \in \dot{B}^{\sigma_1}_{2,\infty}$。由于 $L^1 \hookrightarrow \dot{B}^0_{1,\infty} \hookrightarrow \dot{B}^{-d/2}_{2,\infty}$，本文的初值类比 [32] 更广（例如允许初值在低频有缓慢振荡），且结论更强（给出了双向衰减，而非单边）。

3. 与Navier-Stokes衰减理论的比较：本文借鉴了 [3, 4] 中用于NS方程的 $\dot{B}^{\sigma_1}_{2,\infty}$ 空间，但其应用背景和证明技术截然不同。NS方程的 $\dot{B}^{\sigma_1}_{2,\infty}$ 条件源于其线性化算子（Stokes算子）的纯抛物性；而本文的 $\dot{B}^{\sigma_1}_{2,\infty}$ 条件则源于线性化Oldroyd-B算子的双曲-抛物混合性（damped wave structure），其必要性证明（Lemma 2.1）依赖于对特征方程 $\partial_{tt} u + \eta \Lambda^2 \partial_t u - \frac{\nu_1 \nu_2}{2} \Delta u = 0$ 的精细分析，这是NS理论中不存在的。

4. 局限性与未解决问题：本文的结论局限于 $p$ 的特定范围（$2 \leq p \leq \min\{4, 2d/(d-2)\}$），其上限由临界Besov乘积估计的适用性决定，而 $p=4$ 在 $d=2$ 时的排除则是由于嵌入 $\dot{B}^{0}_{4,1} \hookrightarrow L^4$ 在端点退化。此外，本文未处理 $d=1$ 的情形（此时模型退化），也未讨论解的渐近剖面（asymptotic profile）或高阶展开。从提供内容中无法确认本文结果是否可推广至更一般的 $Q$-tensor 模型或具有分数阶扩散的变体。

4. 含 Biot-Savart 与 Keller-Segel 相互作用的 Dean-Kawasaki 方程：存在性与大偏差

Dean-Kawasaki Equation with Biot-Savart and Keller-Segel Interactions: Existence and Large Deviations

作者: Xiaohao Ji, Yue Sun, Zhengyan Wu

类别: math.PR, math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13479v1

摘要翻译: 我们建立了含奇异相互作用核（包括 Biot–Savart 型与 Keller–Segel 型）的 Dean–Kawasaki 方程的概率弱解与重正化动力学解的存在性。在对平方根噪声系数施加适当正则化后，我们进一步证明了正则化 Dean–Kawasaki 方程的概率弱解满足一个受限的大偏差原理。Biot–Savart 型与 Keller–Segel 型相互作用在 Dean–Kawasaki 方程及其关联骨架方程的 $L^1$ 框架内引入了尺度临界性，从而带来了显著的新挑战。与 [Fehrman, Gess; Invent. Math., 2023] 不同，我们的大偏差分析依赖于一种新颖的指数紧性论证，该论证专门适配于 Dean–Kawasaki 噪声。这一方法结合关联骨架方程的弱-强唯一性结果，使我们得以部分克服奇异相互作用核所引发的临界性。

Comment: 35 pages

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

1. 研究对象和背景

本文研究二维环面 $\mathbb{T}^2$ 上带奇异非局部漂移项的Dean-Kawasaki型保守随机偏微分方程（SPDE）：

其中噪声 $\xi = \sum_{j \in \mathbb{Z}^2} \dot{B}_j e_j$ 是向量值时空白噪声，$\xi_K = \sum_{|j| \leq K} \dot{B}_j e_j$ 是紫外截断，$\circ$ 表示Stratonovich积分；相互作用核 $V[\rho]$ 具有Biot-Savart与Keller-Segel叠加结构：

且 $G$ 是 $\mathbb{T}^2$ 上满足 $-\Delta G * \rho = \rho - \int_{\mathbb{T}^2} \rho$ 的Green函数。该核满足 $V \in L^p(\mathbb{T}^2; \mathbb{R}^2)$ 对所有 $1 \leq p < 2$ 成立，但不满足Ladyzhenskaya–Prodi–Serrin条件 $\|V\|_{L^r_t L^p_x} < \infty$ 且 $\frac{2}{r} + \frac{d}{p} \leq 1$（此处 $d=2$，故要求 $p > 2$），因此完全脱离了[WWZ24, WZ24, FG23]等工作的适用框架。

该方程源自涨落流体力学对布朗粒子系统经验测度的宏观极限建模，其确定性部分 $\partial_t \rho = \Delta \rho - \nabla \cdot (\rho V[\rho])$ 是经典Patlak–Keller–Segel（PKS）与二维涡量输运（Biot-Savart律）的耦合模型，而噪声项 $-\sqrt{\varepsilon}\, \nabla \cdot (\sqrt{\rho} \circ \xi_K)$ 则精确编码了粒子数守恒下的微观涨落——这是Dean-Kawasaki方程区别于加性噪声模型（如[MM25]）的本质特征。从分析角度看，该问题的核心挑战在于三重临界性：（i）空间临界性：$V \in L^{2-}$ 导致 $\nabla \cdot (\rho V[\rho])$ 在 $L^1$ 框架下无法通过标准Sobolev嵌入控制，例如由Gagliardo-Nirenberg不等式 $\|f\|_{L^4} \leq C_{\mathrm{GN}} \|\nabla f\|_{L^2}^{1/2} \|f\|_{L^2}^{1/2} + C \|f\|_{L^2}$ 可得 $\|\rho V[\rho]\|_{L^1} \lesssim \|\nabla \sqrt{\rho}\|_{L^2}^2 \|\rho_0\|_{L^1} + \|\rho_0\|_{L^1}^2$，但系数 $C_{\mathrm{GN}}$ 在临界尺度下无法吸收吸引项 $-\kappa_1 > 0$ 的发散贡献；（ii）时间-空间紧性临界性：为在 $L^1_{[0,T]} L^1_{\mathbb{T}^2}$ 中应用Aubin-Lions引理，需同时控制 $\|\nabla \sqrt{\rho}\|_{L^2_{[0,T]} L^2_{\mathbb{T}^2}}$（空间正则性）与 $\|\rho\|_{W^{\alpha,1}_{[0,T]} H^{-\gamma}_{\mathbb{T}^2}}$（时间正则性），但乘性噪声破坏了高斯型Fernique估计，使得[GHW24]中基于加性噪声的 $W^{\alpha,1} H^{-\gamma}$ 指数矩方法失效；（iii）骨架方程（skeleton equation）的解唯一性临界性：对应零噪声极限的确定性PDE为

其尺度分析显示：对伸缩 $\tilde{\rho}(x,t) = \lambda \rho(\iota x, \tau t)$，要求质量守恒 $\lambda = \iota^2$ 与控制 $g$ 的Cameron-Martin范数不变，导出 $V$ 的临界可积指数为 $p = d = 2$；而实际 $V \in L^{2-}$，故骨架方程处于超临界状态，导致在自然熵空间 $E_{\mathrm{fin}} = \{ \rho \in L^\infty_{[0,T]} L^1_{\mathbb{T}^2} : \|\nabla \sqrt{\rho}\|_{L^2_{[0,T]} L^2_{\mathbb{T}^2}} < \infty, \rho \geq 0 \}$ 中无法建立全局解唯一性。本文旨在填补这一空白：在保留Dean-Kawasaki方程原始保守乘性噪声结构的前提下，首次为Biot-Savart/Keller-Segel型奇异核建立（弱）解存在性，并克服超临界骨架方程带来的大偏差原理障碍，给出首个受限但严格的零噪声大偏差原理。

2. 主要定理或主要结果

本文核心结果分为存在性与大偏差两部分，均建立在初始数据 $\rho_0 \in \mathrm{Ent}(\mathbb{T}^2) := \left\{ \rho \in L^1(\mathbb{T}^2) : \rho \geq 0, \int_{\mathbb{T}^2} \Psi(\rho) < \infty \right\}$ 及吸引强度小性条件

之上，其中 $\Psi(\zeta) = \int_0^\zeta \log(\zeta') \, d\zeta'$，$C_{\mathrm{GN}}$ 为Gagliardo-Nirenberg不等式最优常数。

存在性定理（Theorem 1.2）：对任意固定 $\varepsilon > 0$，设 $\rho_0 \in \mathrm{Ent}(\mathbb{T}^2)$，则（1）对正则化方程（1.4）（即 $\eta > 0$），存在概率弱的重正规化动能解 $\rho^\eta$，且该解亦为（1.4）的经典弱解；（2）对原始方程（1.3）（即 $\eta = 0$），存在概率弱的重正规化动能解 $\rho^\varepsilon$。该解满足质量守恒 $\|\rho^\varepsilon(\cdot, t)\|_{L^1_{\mathbb{T}^2}} = \|\rho_0\|_{L^1_{\mathbb{T}^2}}$ 几乎必然成立，并满足路径熵耗散不等式（2.2）及动能测度 $q$ 的正则性条件（2.3）–（2.4）。

受限大偏差原理（Theorem 1.3）：考虑正则化族 $\{\rho^{\eta(\varepsilon)}\}$，其参数满足渐近标度

其中 $N_K \sim K^{4}$（因 $d=2$）。令 $\mu^\varepsilon = \mathbb{P} \circ (\rho^{\eta(\varepsilon)})^{-1}$ 为其在 $L^1_{[0,T]} L^1_{\mathbb{T}^2}$ 上的分布，则 $\{\mu^\varepsilon\}$ 满足受限大偏差原理，速率泛函为

定义域为 $E_{\mathrm{fin}} = \{ \rho \in L^\infty_{[0,T]} L^1_{\mathbb{T}^2} : \|\nabla \sqrt{\rho}\|_{L^2_{[0,T]} L^2_{\mathbb{T}^2}} < \infty, \rho \geq 0 \text{ a.e.} \}$，其余处 $I(\rho) = +\infty$。具体而言，对任意开集 $G \subset L^1_{[0,T]} L^1_{\mathbb{T}^2}$，有下界

其中限制类 $C_0 = \bigcup_{p' > 2} \{ \rho \in L^\infty_{[0,T]} L^1_{\mathbb{T}^2} : \rho \in L^\infty_{[0,T]} W^{1,p'}_{\mathbb{T}^2} \}$；对任意闭集 $F \subset L^1_{[0,T]} L^1_{\mathbb{T}^2}$，有上界

该受限性源于骨架方程弱强唯一性仅在 $C_0$ 上成立（Proposition 5.7），且 $C_0$ 非最优，将下界延拓至全空间 $L^1_{[0,T]} L^1_{\mathbb{T}^2}$ 仍是开放问题。

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

证明的核心障碍在于：（1）乘性平方根噪声导致标准能量估计失效，尤其在处理 $\nabla \cdot (\rho V[\rho])$ 时，$V \in L^{2-}$ 与 $\rho \in L^2$ 的乘积缺乏足够可积性；（2）超临界骨架方程使弱收敛法（weak convergence method）不可行，因其依赖于骨架方程解的全局稳定性；（3）时间正则性估计受乘性噪声破坏，无法直接应用加性噪声下的Fernique型指数矩技术。

创新点一：路径熵耗散的指数紧性（Exponential Tightness）。作者未采用[FG23]的弱收敛框架，而是发展了一套专为Dean-Kawasaki噪声定制的指数紧性论证。关键在于利用路径熵耗散不等式（2.2）：

通过将右侧鞅项 $\sqrt{\varepsilon} \int_0^t \int_{\mathbb{T}^2} 2 \nabla \sqrt{\rho} \cdot \xi_K$ 视为指数超鞅，并利用其二次变差 $\varepsilon \int_0^t \int_{\mathbb{T}^2} |\nabla \sqrt{\rho}|^2$ 可被左侧拉普拉斯耗散项吸收（源于涨落-耗散结构），导出

此估计是[FG23, FG24]中未出现的新结果，它将解约束在熵耗散有限子集 $E_{\mathrm{fin}}$ 内，并结合质量守恒，得到 $L^2_{[0,T]} W^{1,1}_{\mathbb{T}^2}$ 中的指数矩。

创新点二：时间正则性的Taylor展开技巧。为应用Aubin-Lions引理，需控制 $\|\rho^\varepsilon\|_{W^{\alpha,1}_{[0,T]} H^{-\gamma}_{\mathbb{T}^2}}$ 的指数矩。由于乘性噪声破坏了鞅结构，作者放弃[GHW24]的Gaussian方法，转而对 $\exp\left( \frac{1}{\varepsilon} \|\rho^\varepsilon\|_{W^{\alpha,1}_{[0,T]} H^{-\gamma}_{\mathbb{T}^2}} \right)$ 进行Taylor展开，并对每一阶 $p$-矩 $\mathbb{E}\left[ \|\rho^\varepsilon\|_{W^{\alpha,1}_{[0,T]} H^{-\gamma}_{\mathbb{T}^2}}^p \right]$ 单独估计（$p \geq 1$）。该技术绕开了对噪声高斯性的依赖，是处理保守乘性SPDE时间正则性的新工具。

创新点三：熵方法（Entropy Method）与弱强唯一性的耦合。大偏差下界不采用标准弱收敛法，而是采用熵方法（Lemma 5.1）：对任意 $\rho \in C_0$ 满足 $I(\rho) < \infty$，需构造控制 $g \in L^2_{[0,T]} L^2_{\mathbb{T}^2}$ 使得 $\rho$ 是骨架方程 $\partial_t \rho = \Delta \rho - \nabla \cdot (\rho V[\rho]) - \nabla \cdot (\sqrt{\rho}\, g)$ 的解。由于超临界性导致骨架方程在 $E_{\mathrm{fin}}$ 中无全局唯一性，作者证明：若 $\rho \in C_0$，则所有重正规化动能解必与之重合（Proposition 5.7）。这使得熵方法可在 $C_0$ 上验证——通过选取 $\rho \in C_0$，利用其正则性保证骨架方程解的紧性与极限唯一性，从而规避了无唯一性时的测度收敛困境。此策略将大偏差分析的成败系于弱强唯一性类 $C_0$ 的刻画，是本文最本质的技术突破。

4. 与之前工作的比较

本文与既有工作存在三重根本性差异：

与[FG23, WZ24]的弱收敛法对比：[FG23]为满足Ladyzhenskaya–Prodi–Serrin条件（即 $V \in L^r_t L^p_x$ 且 $\frac{2}{r} + \frac{d}{p} \leq 1$, $p > d$）的Dean-Kawasaki方程建立了完整大偏差原理，其核心是证明骨架方程在熵空间 $E_{\mathrm{fin}}$ 中的稳定性。而本文的Biot-Savart/Keller-Segel核 $V \in L^{2-}$ 严格违反该条件（因 $p \leq 2 = d$），导致[FG23]的弱收敛框架完全失效。本文放弃该路径，转而采用熵方法并引入受限类 $C_0$，代价是下界仅在 $C_0$ 上成立，而非全空间。

与[WWZ24]的存在性工作对比：[WWZ24]在相同Assumption (A1)–(A2)下证明了重正规化动能解存在性，但其假设要求 $V \in L^r_t L^p_x$（$p > d$），仍不覆盖本文核。本文通过精细控制吸引项（利用Gagliardo-Nirenberg不等式与小性条件（1.6））和Biot-Savart项的散度自由性（在（3.7）中 vanish），在超临界框架下重建了熵估计（3.5），这是存在性证明的关键推进。

与[MM25, MM24]的加性噪声模型对比：Martini-Mayorcas研究的是加性噪声Keller-Segel-Dean-Kawasaki模型 $\partial_t \rho = \Delta \rho - \nabla \cdot (\rho V[\rho]) + \nabla \cdot (\sigma \xi)$，其噪声结构与本文保守乘性噪声 $\nabla \cdot (\sqrt{\rho} \circ \xi)$ 本质不同。加性噪声允许使用paracontrolled分析与Fernique估计，而本文必须处理 $\sqrt{\rho}$ 与噪声的耦合，导致时间正则性估计需全新Taylor展开技巧（第6页），且骨架方程形式含 $\nabla \cdot (\sqrt{\rho}\, g)$ 而非 $\nabla \cdot (\sigma g)$，这使得弱强唯一性证明必须针对该非线性控制项设计。

与[GHW24, WW25]的受限LDP对比：[GHW24]在Landau-Lifshitz-Navier-Stokes方程中建立受限LDP，其限制类源于能量等式，而本文的 $C_0$ 源于弱强唯一性，且 $C_0$ 的Sobolev正则性 $W^{1,p'}$（$p'>2$）比[GHW24]的负Sobolev空间更‘强’，反映了Biot-Savart/Keller-Segel核对解正则性的更高要求。此外，[WW25]在准地转方程中使用负Sobolev空间，而本文 $C_0$ 是正则Sobolev空间的并集，体现了保守SPDE特有的正则性提升机制。

综上，本文的推进在于：首次将Dean-Kawasaki理论拓展至物理上重要的Biot-Savart/Keller-Segel奇异核；提出专用于保守乘性噪声的指数紧性与Taylor展开时间正则性技术；并将大偏差下界与弱强唯一性类 $C_0$ 精确绑定。代价是受限性，且 $C_0$ 的最优性（即能否取 $p'=2$）及全空间LDP的可行性仍是开放问题。

5. 零粘性微极流体简化模型的全局 Yudovich 型解

Global Yudovich-type solutions to a reduced model for micropolar fluids with zero viscosity

作者: Francesco Fanelli, Pedro Gabriel Fernández Dalgo

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13478v1

摘要翻译: 本文研究一个二维系统的低正则性适定性问题，该系统作为微极流体动力学的简化模型而获得。在数学层面，该系统呈现了流体二维速度场所满足的类 Euler 方程与标量微旋转场所满足的对流-扩散方程之间的耦合。

Comment: Submitted

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

1. 研究对象和背景

本文研究如下二维简化微极流体模型的低正则性适定性问题：

其中 $u = (u_1, u_2) \in \mathbb{R}^2$ 为速度场，$\Pi \in \mathbb{R}$ 为压力，$m \in \mathbb{R}$ 为标量微旋度（microrotation）场，$\omega := \operatorname{curl} u = \partial_1 u_2 - \partial_2 u_1$ 为其涡量，$\nabla^\perp m := (-\partial_2 m, \partial_1 m)$，常数 $\alpha > 0$, $\kappa > 0$。系统定义在 $(t,x) \in \mathbb{R}_+ \times \Omega$ 上，$\Omega = \mathbb{R}^2$ 或 $\mathbb{T}^2$，初始条件为 $(u,m)|_{t=0} = (u_0,m_0)$，满足 $\operatorname{div} u_0 = 0$。

该模型源自Eringen三维微极流体方程（2）：

通过四步约化获得：(i) 限制为二维水平速度 $U = (u_h, 0)$ 与垂直微旋 $M = (0,0,m)$，且 $\operatorname{div}_h u_h = 0$；(ii) 设 $\kappa = c_a + c_d$, $\alpha = 2\mu_r$；(iii) 人为置零粘性：令 $\mu + \mu_r = 0$；(iv) 人为抹除阻尼项：删去 $+4\mu_r M$。这两步物理上不成立（因 $\mu, \mu_r \geq 0$），但动机在于考察经典微极系统在粘性消失极限下的稳定性——而该极限破坏能量平衡：对（2）作 $L^2$ 能量估计时，右侧 $2\mu_r \operatorname{curl} U \cdot M$ 与 $2\mu_r \operatorname{curl} M \cdot U$ 同号，需粘性与阻尼共同抵消其能量产生效应。约化后系统（1）继承此病态结构：虽保留 $-\kappa \Delta m$ 提供对 $m$ 的耗散，但缺失对 $u$ 的耗散及对 $m$ 的阻尼，导致总能量呈指数增长趋势，构成长期适定性的根本障碍。

问题的重要性在于：它代表了Yudovich理论向“非齐次扰动Euler方程”的首次实质性拓展。经典Yudovich解要求 $\omega_0 \in L^1 \cap L^\infty$，保证 $u$ 为log-Lipschitz连续，从而存在唯一正则流，支撑全局解的存在唯一性。然而，现有工作均未覆盖（1）类耦合结构：(a) 理想MHD、无扩散Boussinesq、非齐次Euler等仅获局部Yudovich型适定性，且依赖于额外的共形正则性假设（如vortex patch正则性）；(b) 含温度扩散的无粘Boussinesq系统虽得全局解，但其正则性要求高于纯Yudovich类（如要求 $\theta_0 \in W^{1,p}$）；(c) Hmidi–Keraani–Rousset对临界扩散Boussinesq的处理虽结构相近，但因涉及非局部算子（Riesz变换），其解空间为比 $L^\infty$ 更精细的 $B^0_{\infty,1}$ 类，无法退化至纯Yudovich框架。因此，本文试图填补的关键空白是：在完全缺失速度耗散与微旋阻尼的极端设定下，能否仅凭 $-\kappa \Delta m$ 的局部耗散机制，结合Yudovich型初值（即 $\omega_0, m_0 \in L^{p_0} \cap L^\infty$, $p_0 \in (1,2)$），建立全局存在唯一性？这本质上是检验Yudovich范式的鲁棒性边界：当Euler方程被一个具有自身耗散但无阻尼的标量场以一阶微分耦合方式扰动时，其低正则性理论是否仍能维持。

2. 主要定理或主要结果

定理1.2（全局Yudovich型适定性）：设 $p_0 \in (1,2)$ 固定，取初值 $(u_0, m_0)$ 满足 $\operatorname{div} u_0 = 0$，且其涡量 $\omega_0 := \operatorname{curl} u_0$ 与微旋 $m_0$ 满足

则系统（1）存在唯一全局弱解 $(u,m)$，满足

其中 $\omega := \operatorname{curl} u$。进一步，令 $q_0 \in (2, +\infty)$ 满足 $1/q_0 = 1/p_0 - 1/2$，则有

此处 $LL(\mathbb{R}^2)$ 为log-Lipschitz空间，其半范数定义为

该解在分布意义下满足（1），且压力梯度 $\nabla \Pi$ 满足 $\nabla \Pi \in L^\infty_{\mathrm{loc}}\big(\mathbb{R}_+;\, L^{p_0}(\mathbb{R}^2) \cap L^\infty(\mathbb{R}^2)\big)$（引理2.3）。值得注意的是，解的正则性完全由初值 $L^{p_0} \cap L^\infty$ 控制，不依赖更高阶范数；时间局部性体现在 $L^\infty_{\mathrm{loc}}$ 和 $L^2_{\mathrm{loc}}$ 中，但全局性由先验估计（第2节）严格保障，无爆破可能。

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

证明的核心障碍在于：系统（1）缺乏对速度场 $u$ 的耗散，且微旋方程中缺失阻尼项，导致标准能量方法失效；同时，$u$ 仅属 $L^{q_0} \cap LL$（$q_0 > 2$），不满足 $L^2$ 可积性，使得对速度差作能量估计成为技术难点。作者的真正创新点在于构建了一个双变量耦合先验控制框架，并发展出适配Yudovich类的精细估计技术。

第一，涡量-微旋耦合控制（Section 2.3）：引入辅助变量 $\Gamma := \omega + \alpha m$，将原系统（7）重写为

关键观察是：尽管右侧含 $\Delta m$ 和 $\omega$，但利用 $\Gamma$ 的 $L^{p_0}$ 范数的极大值原理（因 $u$ 无散，对流项不增范数）及 $-\kappa \Delta$ 的耗散性，可导出全局 $L^{p_0}$ 界：

再结合 $m$ 方程的 $L^{p_0}$ 估计 $\|m(t)\|_{L^{p_0}} \leq e^{C \alpha t} \|m_0\|_{L^{p_0}}$，最终推出 $\|\omega(t)\|_{L^{p_0}}$ 和 $\|m(t)\|_{L^{p_0}}$ 的全局控制。此步骤绕开了对 $\omega$ 单独做 $L^\infty$ 估计的困难，是处理无粘耦合系统的关键突破。

第二，Yudovich流的构造与紧性（Section 3）：采用磨光初值 $u_0^\varepsilon, m_0^\varepsilon$ 求解光滑解 $(u^\varepsilon, m^\varepsilon)$。核心工具是Proposition 2.1给出的Biot-Savart律的精细估计：若 $\omega^\varepsilon \in L^{p_0} \cap L^\infty$，则对应速度 $u^\varepsilon$ 满足

且该估计在 $\varepsilon \to 0$ 下一致成立。由此，$\{u^\varepsilon\}$ 在 $L^\infty_{\mathrm{loc}}(L^{q_0})$ 中有界，$\{m^\varepsilon\}$ 在 $L^\infty_{\mathrm{loc}}(L^{p_0}) \cap L^2_{\mathrm{loc}}(H^1)$ 中有界，应用Aubin-Lions引理可提取收敛子列，其极限 $(u,m)$ 满足（1）的弱形式。

第三，Yudovich类唯一性（Section 4）：借鉴Yudovich原始论证，考虑两解 $(u^{(1)}, m^{(1)})$, $(u^{(2)}, m^{(2)})$ 的差 $\delta u = u^{(1)} - u^{(2)}$, $\delta m = m^{(1)} - m^{(2)}$。标准能量法需 $\delta u \in L^2$，但此处不成立。作者的创新在于利用Biot-Savart律将 $\delta u$ 表示为 $\delta \omega$ 的卷积，并借助log-Lipschitz性控制流映射差。关键引理是：若 $u^{(i)} \in LL$，则其诱导的流 $\Phi^{(i)}$ 满足

而 $\|\nabla u^{(i)}(s)\|_{L^\infty}$ 由 $\|\omega^{(i)}(s)\|_{L^{p_0} \cap L^\infty}$ 控制（通过Calderón-Zygmund）。最终，通过估计 $\|\delta \omega(t)\|_{L^{p_0}}$ 和 $\|\delta m(t)\|_{L^{p_0}}$ 的Grönwall型不等式完成唯一性证明。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度：

（1）与经典微极流体理论比较：Galdi–Rionero [15] 及Chen [6] 均要求 $\mu + \mu_r > 0$ 和 $c_0, c_d, c_a > 0$，从而获得Navier-Stokes型能量恒等式 $\frac{d}{dt}(\|U\|_{L^2}^2 + \|M\|_{L^2}^2) + \text{dissipation} = 0$，支撑有限能量弱解。本文则彻底放弃所有粘性与阻尼，使能量不再守恒甚至不耗散，迫使作者发展出基于 $L^{p_0} \cap L^\infty$ 的新控制机制，这是对微极理论适用边界的实质性推进。

（2）与非齐次Euler类系统比较：文献[11,16,17]处理理想MHD或Boussinesq时，仅获局部适定性，且依赖J.-Y. Chemin的共形正则性（如 $\omega_0$ 在曲线附近有界变差）。本文则在纯Yudovich初值（仅 $L^{p_0} \cap L^\infty$）下实现全局性，代价是必须引入微旋的扩散项 $-\kappa \Delta m$ 作为唯一的正则化源。这表明：对Euler方程的扰动，若扰动项本身具备耗散性，则可补偿缺失的粘性，换取全局性。

（3）与扩散Boussinesq比较：文献[18,20,24]中，温度方程 $\partial_t \theta + u \cdot \nabla \theta - \kappa \Delta \theta = 0$ 无源项，动量方程右侧为 $-\nabla \theta$（无导数）。而本文系统（1）的耦合更奇异：$m$ 方程含源项 $\alpha \omega$（涡量直接驱动微旋），动量方程右侧为 $-\alpha \nabla^\perp m$（微旋梯度直接驱动速度），形成双向一阶微分耦合。这导致标准对流-扩散估计失效，必须引入 $\Gamma = \omega + \alpha m$ 进行解耦，这是本文区别于所有前述工作的核心技术特征。

（4）与Hmidi–Keraani–Rousset [19] 比较：其临界扩散Boussinesq系统结构最接近（1），但因使用Riesz变换表示速度，导致 $\nabla u$ 为 $\theta$ 的非局部函数，迫使解空间取 $B^0_{\infty,1}$（比 $L^\infty$ 更强）。本文因 $\nabla^\perp m$ 是局部算子，且 $u$ 由Biot-Savart律局部表示，故能退化至纯 $L^\infty$ 涡量框架。但代价是：其唯一性证明无法直接套用[19]的非局部估计，而必须重构基于log-Lipschitz流的Yudovich型论证。从提供内容中无法确认本文结果是否可推广至 $\omega_0 \in L^\infty$（即 $p_0 \to 1^+$ 极限），作者亦未显式展开此比较。

6. 非凸有界区域中加权Wasserstein度量的演化变分不等式

The evolution variational inequality for weighted Wasserstein metrics in non-convex bounded domains

作者: Kyogo Murai

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13420v1

摘要翻译: 本文在不假设区域凸性的前提下，建立了加权Wasserstein距离的演化变分不等式。得益于该演化变分不等式，我们得以在不仅凸而且非凸的区域中运用加权Wasserstein度量进行若干论证。最终，我们将该演化变分不等式应用于加权Wasserstein度量下的极小化运动，从而在非凸区域中获得Keller–Segel系统与Cahn–Hilliard型方程的弱解。去除凸性假设的关键在于对边界积分的控制。为处理边界积分，我们使用了边界上函数的估计、Sobolev迹嵌入以及Kato不等式的变体。由此，边界积分可被已知的良好项所吸收。

Comment: 23 pages

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

1. 研究对象和背景

本文研究的核心对象是加权Wasserstein度量 $ W_m $ 下的概率测度演化问题，其定义基于连续性方程（Continuity Equation, CE）的极小作用原理：

其中 $ \Omega \subset \mathbb{R}^d $ 是具有光滑边界的有界区域（不假设凸性），$ m : [0,M] \to [0,\infty) $ 是满足特定正则性与凹性的迁移率函数（mobility），$ \mathrm{CE}(0,1; \mu_0 \to \mu_1) $ 表示所有满足弱意义下 $ \partial_t \mu + \nabla \cdot \nu = 0 $、且初终值为 $ \mu(\cdot,0)=\mu_0 $、$ \mu(\cdot,1)=\mu_1 $ 的非负Radon测度对 $ (\mu,\nu) $。当 $ m(r)=r $ 时，该度量退化为经典 $ L^2 $-Wasserstein距离 $ W_2 $，由Benamou–Brenier定理保证其与Monge–Kantorovich公式等价；但一般情形下，$ W_m $ 无矩表示（no moment formulation），因而无法直接构造平移扰动 $ \mu_a = T_{a\xi}\#\mu $ 来推导Euler–Lagrange方程——这是JKO方案在加权Wasserstein空间中实施的根本障碍。为克服此障碍，现有理论（如[12,17,19,20]）依赖于流交换引理（flow interchange lemma），其关键前提是演化变分不等式（Evolution Variational Inequality, EVI）成立：存在 $ \lambda \in \mathbb{R} $，使得对所有 $ u,v \in \mathrm{Dom}(F) $ 满足 $ W_m(u,v) < \infty $，有

其中 $ S_t $ 是由泛函 $ F $ 生成的梯度流半群。EVI的证明归结为对辅助能量 $ A_h(z) := \int_\Omega m_\varepsilon(\rho_h(z)) |\nabla \phi_h(z)|^2 dx $ 的微分估计，其中 $ \rho_h(z) = S_{hz} \rho(z) $ 是由热-迁移率耦合方程 $ \partial_t \rho = \delta \Delta \rho + \nabla \cdot (m_\varepsilon(\rho) \nabla \phi) $ 驱动的流，$ \phi_h(z) $ 是对应连续性方程 $ \partial_z \rho_h = -\nabla \cdot (m_\varepsilon(\rho_h) \nabla \phi_h) $ 的势函数，满足Neumann边界条件 $ \nabla \phi_h \cdot n = 0 $。在[12]中，EVI的证明依赖于对边界项 $ J_5 := \int_{\partial \Omega} m_\varepsilon(\rho_h(z)) \nabla(|\nabla \phi_h(z)|^2) \cdot n \, dS $ 的控制；当 $ \Omega $ 凸时，由凸性导致的二阶导数符号性质可得 $ J_5 \leq 0 $，从而被忽略。然而，非凸区域中该边界项可正可负，且无法被几何条件支配，成为阻碍EVI推广至一般有界区域的根本技术瓶颈。本文旨在填补这一空白：在完全放弃区域凸性假设的前提下，建立适用于非凸有界区域的EVI，从而为Keller–Segel系统与Cahn–Hilliard型方程在非凸域中的弱解存在性提供统一的变分框架。问题的重要性在于：大量实际模型（如细胞趋化、相分离）天然定义于非凸几何（如生物组织褶皱、多孔介质通道），而既有加权Wasserstein梯度流理论因凸性限制无法覆盖此类情形，导致变分结构与解的存在性之间存在断裂。

2. 主要定理或主要结果

在区域 $ \Omega \subset \mathbb{R}^d $ 为具有光滑边界的有界区域（不假设凸性）、迁移率 $ m $ 满足条件 (M-LSC) 与 (M-A)（见原文第7页）的前提下，本文证明了核心微分不等式（Theorem 3.3）：对任意 $ \phi \in C^\infty(\bar{\Omega}) $ 满足 $ \nabla \phi \cdot n = 0 $ 在 $ \partial \Omega $ 上，以及任意 $ \delta > 0 $、$ \varepsilon > 0 $，存在仅依赖于 $ \Omega $ 和 $ \varepsilon $ 的充分小常数 $ \sigma > 0 $，使得对所有 $ z \in [0,1] $、$ h \in (0,1) $，有

其中 $ A_h(z) = \int_\Omega m_\varepsilon(\rho_h(z)) |\nabla \phi_h(z)|^2 dx $，$ F(u) = \int_\Omega u \phi \, dx + \delta \int_\Omega U_\varepsilon(u) \, dx $，且 $ \lambda_\delta = -\frac{1}{\delta} \|\nabla \phi\|_{L^\infty(\Omega)}^2 \sup_{0<r<M} (m_\varepsilon(r)|m_\varepsilon''(r)|) - \|\nabla^2 \phi\|_{L^\infty(\Omega)} \sup_{0<r<M} |m_\varepsilon'(r)| - \delta C / \sigma \leq 0 $，$ C > 0 $ 为仅依赖于 $ \Omega $ 的常数。由此导出的EVI（1.3）中参数 $ \lambda $ 由 $ \lambda_\delta $ 给出，且其负性由 $ \delta $ 的选择与 $ \sigma $ 的小性共同保障。作为应用，作者在Section 4中证明：对Keller–Segel系统 $ \partial_t u = \Delta u - \nabla \cdot (u \nabla v), \, \tau \partial_t v = \Delta v + u - v $（其中迁移率 $ m(u)=u $）与Cahn–Hilliard型方程 $ \partial_t u = \nabla \cdot (m(u) \nabla (-\Delta u + f'(u))) $（其中 $ m(u)=u $ 或 $ m(u)=u(1-u) $），在非凸有界区域 $ \Omega $ 上，通过加权Wasserstein空间中的JKO方案可构造出满足质量守恒、能量耗散与弱连续性方程的弱解 $ u \in L^\infty(0,T; L^1(\Omega)) \cap L^2(0,T; H^1(\Omega)) $，且解序列在 $ C([0,T]; \mathcal{P}_2(\Omega)) $ 中收敛。

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

本文的核心技术障碍是：在非凸区域中控制边界积分 $ J_5 = \int_{\partial \Omega} m_\varepsilon(\rho_h) \nabla(|\nabla \phi_h|^2) \cdot n \, dS $，该积分在凸情形下自动非正，但在非凸情形下可能主导整个估计并破坏EVI。作者的真正创新点在于构建了一套三重嵌套的边界-体积分析框架，将 $ J_5 $ 分解、估计并最终吸收进已有的负项中。第一层是边界法向导数的几何控制（Lemma 2.1）：利用 $ \nabla \phi_h \cdot n = 0 $ 在 $ \partial \Omega $ 上，结合文献[16]的估计，得到 $ \nabla(|\nabla \phi_h|^2) \cdot n \leq C_\Omega |\nabla \phi_h|^2 $，将 $ J_5 $ 控制为 $ \int_{\partial \Omega} m_\varepsilon(\rho_h) |\nabla \phi_h|^2 dS $。第二层是Sobolev迹嵌入与分数阶Gagliardo–Nirenberg不等式（Lemmas 2.2 & 2.3）：将边界积分上提至体积分，即

其中 $ s \in (1/2,1) $。第三层是Kato型不等式与Hessian负项的协同吸收（Lemma 2.4）：展开梯度项后，出现需控制的新增体积分 $ \int_\Omega m_\varepsilon(\rho_h) |\nabla |\nabla \phi_h||^2 dx $，而经典证明中被忽略的负项 $ -\int_\Omega m_\varepsilon(\rho_h) |\nabla^2 \phi_h|^2 dx $ 此刻成为关键资源；由Kato不等式 $ |\nabla |\nabla \phi_h|| \leq |\nabla^2 \phi_h| $ a.e.，该新增项被完全吸收进负Hessian项中，即

这一吸收机制是全文最精妙之处：它将原本因放弃凸性而“释放”出的边界正贡献，通过分数阶嵌入转化为体积分，并借由Kato不等式将其“锁回”到已知的负Hessian结构内，从而在不损失能量耗散率的前提下，彻底摆脱了对区域凸性的依赖。整个框架的代价是引入了更严格的迁移率条件 (M-A)，以确保 $ m_\varepsilon $ 的二阶导数与一阶导数的比值有界，从而在吸收过程中系数可控。

4. 与之前工作的比较

本文与先前工作（特别是[12,17,19,20]）的本质区别在于凸性假设的完全移除，这带来四个维度的实质性推进与相应代价。第一，适用范围维度：[12]等文要求 $ \Omega $ 为凸区域（甚至需 $ C^2 $ 边界曲率非负），而本文仅需 $ \Omega $ 为具有光滑边界的有界区域，涵盖所有星形域、环形域及任意复杂拓扑的非凸区域，极大扩展了加权Wasserstein梯度流的应用场景。第二，技术假设维度：[12]中EVI的证明仅需迁移率 $ m $ 满足 (M-LSC)（Lipschitz且 $ m^2 $ 半凸），而本文因需控制额外边界衍生项，必须追加条件 (M-A)：$ \sup_{r\geq 0} |m_\varepsilon'(r)|^2 / |m_\varepsilon''(r) m_\varepsilon(r)| < \infty $。该条件排除了过于平坦的迁移率（如 $ m(r) = r $ 在 $ r=0 $ 附近），但保留了所有物理相关模型（如 $ m(r)=r^\alpha $, $ m(r)=\log(1+r) $, $ m(r)=[r(M-r)]^\alpha $）。第三，时间尺度与正则性维度：在[12]中，EVI的 $ \lambda $ 仅依赖于 $ \phi $ 和 $ m $，而本文的 $ \lambda_\delta $ 显含参数 $ \delta $ 与 $ \sigma $，且需 $ \sigma $ 充分小以压倒边界项；这意味着在非凸情形下，EVI的耗散率 $ \lambda $ 对扩散参数 $ \delta $ 更为敏感，且其最优值需通过精细平衡获得，而非像凸情形下可直接取几何下界。第四，解的构造维度：在[14,15,11]中，经典Wasserstein空间的JKO方案已推广至非凸区域，因其依赖矩表示与平移扰动，边界处理相对直接；而本文首次将同等广义性赋予加权Wasserstein空间，使得Keller–Segel与Cahn–Hilliard型方程的弱解可在非凸域中通过纯变分方法（而非PDE先验估计）获得，其解满足与凸情形完全一致的能量耗散律 $ \frac{d}{dt} E(u_t) \leq -\lambda W_m(u_t, u_s)^2 $。从提供内容中无法确认本文结果是否可进一步推广至仅具Lipschitz边界的区域，亦未显式展开与[9]（其处理的是不同方程的边界控制）在函数空间和尺度上的定量比较。

7. Keller-Segel 系统的谱不稳定性与温和解的非唯一性

Spectral instability and non-uniqueness of mild solutions for the Keller-Segel system

作者: Eliseo Luongo, Umberto Pappalettera

类别: math.AP, math.SP

发布日期: 2026-05-13

arXiv ID: 2605.13592v1

摘要翻译: 我们证明：对于维数 $n \in \{3,\dots,9\}$ 及超临界范围 $q\in [1,\frac{n}{2})$，抛物-椭圆型 Keller-Segel 模型所对应的柯西问题在 $L^q(\mathbb{R}^n)$ 中局部不适定。该非唯一性由自相似变量下的不稳定性机制所驱动，其精神与 Jia 和 Šverák 为三维 Navier-Stokes 方程所提出的纲领一致。

Comment: 48 pages, comments welcome!

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

1. 研究对象和背景

本文研究经典抛物-椭圆型Keller-Segel系统在超临界Lebesgue空间中的适定性问题。该系统描述细菌群体在自分泌化学信号驱动下的趋化运动，其数学模型为：

其中未知量 $c : [0,T] \times \mathbb{R}^n \to \mathbb{R}$ 表示细胞密度，$v$ 是由浓度梯度诱导的漂移速度场，通过Riesz变换 $\nabla(-\Delta)^{-1}$ 非局部耦合于 $c$。该系统具有缩放不变性：若 $c(t,x)$ 是解，则对任意 $\lambda > 0$，函数 $c_\lambda(t,x) := \lambda^2 c(\lambda t, \lambda^2 x)$ 也是解。由此导出临界空间为 $C([0,T]; L^{n/2}(\mathbb{R}^n))$：当 $q > n/2$（亚临界）时，扩散项主导，已知存在局部唯一温和解；当 $q = n/2$（临界）时，存在有限时间爆破准则（如 $\|c_0\|_{L^{n/2}}$ 过大则爆破）；而当 $q < n/2$（超临界）时，非线性项主导，理论上预期不适定，但此前尚无严格证明。

本文聚焦于超临界范围 $q \in [1, n/2)$、维数 $n \in \{3,\dots,9\}$ 下的局部不适定性，核心目标是填补以下关键空白：(1) 证明存在非零初值 $c_0 \in L^q(\mathbb{R}^n)$，使得Cauchy问题存在至少两个不同的温和 $L^q$-解；(2) 揭示该非唯一性的深层机制——并非源于弱解框架下的能量缺陷或缺乏紧性，而是源于自相似变量下线性化算子的谱不稳定性；(3) 在无外部强迫、无计算机辅助、纯解析的条件下，首次为Keller-Segel系统实现Jia–Šverák纲领（原为三维Navier-Stokes方程设计），从而将谱不稳定性这一分析工具从流体力学成功移植至趋化模型。该问题的重要性在于：它挑战了“超临界不适定即意味着解不存在”的朴素直觉，揭示了一类新型病态性——解存在但不唯一，且其根源可精确追溯至自相似结构中某个正实部特征值所激发的指数增长模态；同时，它为理解生物趋化聚集、恒星动力学中引力坍缩等物理过程的内在不稳定性提供了严格的PDE理论基础。特别地，由于Keller-Segel系统与Navier-Stokes vorticity方程共享相同的缩放律、非局部结构及梯度非线性项 $\nabla c \cdot \nabla(-\Delta)^{-1}c$，其谱不稳定性分析面临与后者同等严峻的技术障碍，尤其是线性化算子 $L_{\bar{\Xi}}[\Phi] = \Delta \Phi + \frac{\theta}{2} \cdot \nabla \Phi + \Phi - \nabla \bar{\Xi} \cdot \nabla(-\Delta)^{-1}\Phi - \nabla \Phi \cdot \nabla(-\Delta)^{-1}\bar{\Xi} + 2\bar{\Xi} \Phi$ 的非自伴性、强奇异性与高维退化性，使得标准半群理论与谱扰动方法完全失效。

2. 主要定理或主要结果

本文在 $n \in \{3,\dots,9\}$ 维欧氏空间中，对超临界指数 $q \in [1, n/2)$，严格证明了Keller-Segel系统温和解的非唯一性。主要结果如下：

定理1.2（非唯一性）：对任意 $n \in \{3,\dots,9\}$ 和任意 $q \in [1, n/2)$，存在非零初值 $c_0 \in L^q(\mathbb{R}^n)$、时间 $T > 0$，以及两个互异的温和 $L^q$-解 $c_1, c_2 \in C([0,T]; L^q(\mathbb{R}^n))$，均满足定义1.1中的积分方程

其中 $v_i = \nabla(-\Delta)^{-1} c_i$。

定理1.3（谱不稳定性）：对任意 $n \in \{3,\dots,9\}$，固定 $1 < \eta < n/2 < \gamma$ 及 $\bar{\lambda} > 0$，存在一个与 $\tau$ 无关的径向对称自相似剖面 $\bar{\Xi} \in C^2_{\mathrm{rad}}(\mathbb{R}^n)$，使得其对应的线性化算子 $L_{\bar{\Xi}}$ 在加权空间 $X_{\eta,\gamma} := L^\eta_{\mathrm{rad}}(\mathbb{R}^n) \cap L^\gamma_{\mathrm{rad}}(\mathbb{R}^n)$ 上存在一个正实部特征值 $\lambda \in (0, \bar{\lambda}]$，且谱 $\sigma(L_{\bar{\Xi}}) \subset \{ \mu \in \mathbb{C} : \operatorname{Re}(\mu) \leq \lambda \}$。相应特征函数 $\phi_\lambda$ 满足指数衰减估计 $|\phi_\lambda(\theta)| \lesssim e^{-c|\theta|}$，$c>0$。

定理3.13（最小化谱假设）：存在仅依赖于 $n$ 的隐式阈值 $\alpha_0 = \alpha_0(n) > 0$，使得对所有 $\alpha > \alpha_0$，对应剖面 $\bar{\Omega}(\xi) = u_\alpha(|\xi|)$ 的线性化算子 $L_\alpha$ 在空间 $Y^{\eta,\gamma}_\nabla$ 上具有正特征值；而当 $\alpha \in [0, 1/(16n))$ 时，$L_\alpha$ 稳定（即谱位于左半平面）。该定理允许覆盖端点 $q=1$，因映射 $A: X_{1,\gamma} \to Y^{1,\gamma}_\nabla$ 非同构，但其逆 $A^{-1}$ 连续有界。

非唯一性构造细节：所构造的两个解均为径向对称，其初值 $c_0$ 在原点附近具有 $|x|^{-2}$ 奇异性：$c_0(x) \sim \frac{2\ell_\alpha (n-2)}{|x|^2}$，$\ell_\alpha \in (0,2]$；在空间无穷远处被截断以保证 $L^q$ 可积性。非唯一性程度可推广：对任意 $k \in \mathbb{N}$ 及互异正参数 $\{a_j\}_{j=1}^k$，可构造 $k$ 个互异的温和 $L^q$-解。

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

本文的核心创新在于将Jia–Šverák纲领完整、严格地应用于Keller-Segel系统，并克服了三大根本性技术障碍：

障碍一：非自伴线性化算子 $L_{\bar{\Xi}}$ 的谱分析。直接分析 $L_{\bar{\Xi}}$ 极其困难，因其既非自伴也非正常算子，且系数含强增长项 $\frac{\theta}{2} \cdot \nabla$。作者的关键突破是引入径向约化与加权变换：利用线性同构 $A[c](t,y) = \frac{1}{2\omega_{n-1}|y|^n} \int_{B(0,|y|)} c(t,x) \, dx$ 将原系统（1.1）映射为 $\mathbb{R}^{n+2}$ 上的局部PDE（1.5）：$\partial_t w = \Delta w + y \cdot \nabla(w^2) + 2n w^2$。此变换将 $L_{\bar{\Xi}}$ 的谱问题等价转化为 $\mathbb{R}^{n+2}$ 上算子 $L_{\bar{\Omega}}$ 的谱问题，并进一步通过加权 $L^2_\pi$ 空间 $\pi(\xi) = |\xi|^2$ 将 $L_{\bar{\Omega}}$ 重写为自伴Sturm-Liouville形式 $L^\pi_\alpha$。由此，谱不稳定性问题被归约为研究一个二阶ODE的零点个数：对剖面 $u_\alpha$ 对应的线性化方程

其中 $f(0)=1$, $f'(0)=0$，证明存在 $\rho_0 > 0$ 使得 $f(\rho_0) = 0$。这正是Sturm振荡理论的应用前提。

障碍二：非显式剖面 $u_\alpha$ 的渐近行为精确定量。由于 $u_\alpha$ 无闭式解，无法直接求解上述ODE。作者发展了一套融合动力系统与比较原理的新方法：首先对 $\rho$ 进行尺度变换 $\tilde{u}_\alpha(\rho) = u_\alpha(\rho / \sqrt{\alpha}) / \alpha$，使方程在 $\alpha \to \infty$ 极限下退化为自治系统（1.9）：$\ddot{z} + (n-4 + 2z) \dot{z} - 2(n-2)(z - z^2) = 0$，其中 $z(t) = e^{2t} \tilde{u}(e^t)$。利用Poincaré-Bendixson定理与Bendixson-Dulac判据，严格证明了当 $n \in \{3,\dots,9\}$ 时，$z(t)$ 收敛到平衡点 $(z, \dot{z}) = (1,0)$，从而得到关键极限 $\lim_{\rho \to \infty} \rho^2 \tilde{u}(\rho) = 1$ 和 $\lim_{\rho \to \infty} \rho^3 \tilde{u}'(\rho) = -2$。这些极限是应用Picone比较原理 [Pic10] 的基石。

障碍三：Picone比较原理的精准应用。作者将变换后的 $f$ 方程与标准振荡方程 $y'' + \frac{A}{\rho} y' + \frac{B}{\rho^2} y = 0$ 比较，其中 $A = n+1$, $B = -2(n-2)$。Picone原理指出，若 $(A-1)^2 - 4B < 0$，则解振荡。代入得判据 $n^2 - 10n + 17 < 0$，其正根为 $5 \pm 2\sqrt{2} \approx 2.17, 7.83$，故 $n \in \{3,\dots,7\}$ 满足；结合 $n=8,9$ 的单独验证（见原文Subsubsection 2.3.1），最终确定维数范围 $n \in \{3,\dots,9\}$。此判据是全文维数限制的终极来源，且 $n=2$ 和 $n \geq 10$ 时判据失效，作者推测此时系统稳定。

障碍四：非唯一解的定位与正则性恢复。谱不稳定性产生的是定义在 $(-\infty, \log T]$ 上的“古代解”，其空间衰减不足（如 $\sim |\xi|^{-2}$），不属 $Y^{q,\gamma}_\nabla$。作者设计了一个精巧的两步定位：(i) 将初值分解为 $\tilde{w}_0 = w_0 + h_0$，其中 $w_0$ 在原点附近保留 $|y|^{-2}$ 奇异性，在无穷远快速衰减以保证 $w_0 \in Y^{q,\gamma}_\nabla$；(ii) 求解一个带奇异自相似势的线性化方程（1.13）：$\partial_t h = \Delta h + 4n \bar{w} h + 2y \cdot \nabla(\bar{w} h) + f(h)$，其适定性依赖于谱不稳定强度 $\lambda_\alpha < 1 - n/(2r)$。该条件可通过选择足够大的 $\alpha$ 满足（由Theorem 3.13保证）。整个框架的关键在于，通过控制 $|y|h$ 而非仅 $h$，绕开了 $L^1$ 中Hardy不等式失效的困境，从而将结果延拓至 $q=1$。

4. 与之前工作的比较

本文与既有工作存在本质性推进，主要体现在五个维度：

1. 与Jia–Šverák纲领的原始工作 [JŠ15] 比较：[JŠ15] 为三维Navier-Stokes方程提出该纲领，但未给出任何谱不稳定性证明。本文是首个在无外部强迫（unforced）、无计算机辅助（computer-assisted）、纯解析（analytic）条件下，对该纲领在另一重要PDE（Keller-Segel）上实现的范例。与[JŠ15]不同，本文的线性化算子 $L_{\bar{\Xi}}$ 含有梯度非线性项 $\nabla \Phi \cdot \nabla(-\Delta)^{-1}\bar{\Xi}$，导致其比Navier-Stokes的vorticity线性化更复杂，且无法直接套用Vishik [Vis18a,b] 或Albritton–Bruè–Colombo [ABC22] 的方法（后者需人为添加强迫项）。本文的纯解析证明为该纲领的普适性提供了强有力的新证据。

2. 与Navier-Stokes相关工作比较：[GŠ23] 仅提供数值证据；[ABC22] 的证明依赖于精心设计的外部强迫；[HWY26] 宣称的无强迫证明仍为计算机辅助。本文则完全规避了这些手段，代价是维数受限于 $n \leq 9$，而Navier-Stokes的谱不稳定性在三维已被证实。这表明Keller-Segel的谱结构对维数更敏感，其不稳定性是一种“高维脆弱现象”。

3. 与Keller-Segel自身适定性文献比较：经典结果如[Kar99, BHN94] 仅处理亚临界 $q > n/2$；[BM14] 处理二维临界情形；[CPZ04] 给出临界爆破准则。本文首次在超临界 $q < n/2$ 范围内建立不适定性，且明确指出其机制是“存在性+非唯一性”而非“不存在性”。这修正了以往认为超临界必然导致解不存在的误解。

4. 与非线性热方程工作 [GHLL25] 比较：[GHLL25] 是此前唯一无强迫、无计算机的Jia–Šverák实现，但其方程 $\partial_t u = \Delta u + u^p$ 缺少梯度非线性项，其ODE分析可直接引用[HW82, Nai06]。本文的梯度项 $2\rho u_\alpha f' + 2\rho u'_\alpha f$ 彻底破坏了这种可引用性，迫使作者发展全新的动力系统-Picone联合分析法，这是本文最核心的方法论贡献。

5. 与定位技术文献比较：[JŠ15] 指出定位的难点在于“余项方程含临界奇异的低阶项”，本文对此给出了首个完整的解析解法：通过引入加权空间 $Y^p_1$ 和控制 $|y|h$，建立了带奇异自相似势的线性化方程的适定性理论（引理5.2），其关键条件 $\lambda_\alpha < 1 - n/(2r)$ 明确量化了谱不稳定强度与定位可行性的尺度关系。这一技术为未来处理其他具奇异势的PDE提供了新范式。从提供内容中无法确认本文结果是否可推广至 $n \geq 10$ 或 $q \geq n/2$；作者亦未显式展开与[CC57, Wol92] 等天体物理文献的比较。

8. 精确守恒与Onsager阈值：不可压缩Navier–Stokes方程的离散外微积分理论

Exact conservation and the Onsager threshold: a discrete exterior calculus theory for incompressible Navier–Stokes

作者: Peter Korn

类别: math.AP, math.NA

发布日期: 2026-05-13

arXiv ID: 2605.13048v1

摘要翻译: 我们基于闭Riemann流形上棱柱型Delaunay–Voronoi网格的离散外微积分，发展了一套针对不可压缩Euler方程与Navier–Stokes方程的保结构离散化严格理论。核心结果是一个选择原理：离散层面的精确代数守恒不仅是一种保真性性质，而且无条件地排除了其他离散方法所能达到的整类弱解。我们在四个正则性框架下确立该原理。\emph{光滑解}：收敛速率为 $\mathcal{O}(h^{\min(r_{\rm rec},\,r_\star)}\,|\log h|^{β_d})$，且在粘性系数 $ν\ge 0$ 上一致成立，其中 $β_3 = 0$、$β_2 = 1$；在一般网格上为一阶收敛，在具有形心邻近性与重构对称性的网格上为二阶收敛。\emph{Leray–Hopf弱解框架}：子列的 $L^2$ 极限是粘性系统的弱解。\emph{无粘测度值框架}：极限是守恒型测度值Euler解；其浓度缺陷在Onsager阈值 $α> 1/3$ 之上消失，\emph{前提是离散解在该区域具有一致的 $C^{0,α}$ 模界}。\emph{耗散框架}：不存在任何子列收敛至任一正则性下的能量耗散型Euler解，这一结构性排除源于离散能量的精确守恒，并使该格式区别于所有Galerkin方法与有限体积法。区间 $1/3 < α< 1$ 构成一个间隙——在此区间内能量守恒与无缺陷收敛均成立，但唯一性问题仍未解决，从而凸显了无粘流体动力学中最核心的公开问题。

Comment: 无

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

1. 研究对象和背景

本文研究对象是定义在闭定向黎曼流形 $\Omega$（维数 $d = 2$ 或 $3$）上的不可压缩Euler与Navier–Stokes方程的向量不变形式：

其中 $v : [0,T] \times \Omega \to T\Omega$ 是速度场，$p$ 是压力，$\nu \ge 0$ 为运动粘性系数。该系统具有深层几何结构：其动力学由Lie导数 $\mathcal{L}_v \omega = d(i_v \omega) + i_v (d\omega)$ 控制，能量守恒源于Lamb矢量 $Q(v,v) = i_v \omega$ 的反对称性 $\langle v, Q(v,v) \rangle = 0$，而Kelvin环量守恒则源于 $d$ 与 $i_v$ 的交换关系。在连续层面，这些结构对应于de Rham复形 $\cdots \xrightarrow{d} \Lambda^k(\Omega) \xrightarrow{d} \Lambda^{k+1}(\Omega) \xrightarrow{d} \cdots$ 上的微分算子、Hodge星算子 $\star$ 及内积 $i_v$ 的协调作用。

然而，标准数值方法（如Galerkin型有限元、有限体积法）在离散化过程中不可避免地破坏这些结构：其离散能量仅满足近似恒等式 $\frac{d}{dt} E_{\rm kin}^h \approx -\nu \|\nabla v_h\|_{L^2}^2 + \mathcal{O}(h^r)$，导致高波数模态的能量耗散或反向级联失真；更严重的是，它们无法区分物理上允许与禁止的弱解——例如，经典格式可收敛至任意Onsager正则性 $C^{0,\alpha}$ 下的能量耗散Euler解（$\alpha \le 1/3$），而此类解虽数学上满足分布意义下的方程，却违反热力学第二定律所要求的能量守恒。这一缺陷在长时间积分、湍流模拟及气候模型中引发不可控的偏差。

现有结构保持格式亦存在根本局限：FEEC框架[1,2]严格构造了离散de Rham复形并保证拓扑精确性，但未处理非线性对流项 $v \cdot \nabla v$ 的几何实现；Elcott等[22]与Mullen–Pavlov[43,44]在单纯形网格上实现了离散能量与环量的精确守恒，但缺乏收敛性分析；EMAC格式[12]通过混合对流与旋转形式实现能量守恒，却牺牲了Kelvin环量；Zhang等[55]以双场变量代价获得螺度守恒；mimetic有限差分[9]满足离散Green恒等式，但非线性问题的收敛理论缺失。所有前述工作均未建立覆盖全正则性谱（从光滑解到测度值解）的统一收敛理论，亦未揭示离散守恒律作为连续极限选择原理的本质作用。

本文试图填补的核心空白在于：构建一个兼具代数精确守恒性与全正则性谱收敛性的离散框架，并严格证明：精确的离散能量恒等式 $\frac{d}{dt} E_{\rm kin}^h = -\nu \|\delta_h \tilde{d}_1 v_h\|_{L^2_h}^2$ 不仅是精度指标，更是对连续极限的强约束——它系统性地排除整类非物理弱解，从而在数学上实现类似熵条件对激波的筛选功能。具体而言，需回答：（1）在何种网格与算子设计下，能将所有近似误差严格集中于Hodge星算子 $M_k$，使外微分 $d$ 与内积 $i_v$ 的离散实现保持组合精确性？（2）如何基于Delaunay–Voronoi正交网格的几何特性，构造满足 $\langle v_h, I_v(\tilde{D}_1 v_h) \rangle_1 = 0$ 的离散Lamb矢量，从而继承反对称性？（3）在此框架下，能否建立统一证明：光滑解的最优收敛率、Leray–Hopf弱解的存在性、Onsager阈值 $\alpha > 1/3$ 处测度值解的浓度缺陷消失、以及能量耗散Euler解的结构性排除？这些问题直指计算流体力学基础——如何使离散模型不仅‘看起来像’连续方程，而且‘逻辑上强制’其极限满足物理第一性原理。

2. 主要定理或主要结果

本文在满足Assumption 2.1（拟均匀性、形状正则性、Delaunay性质、有界价）的prismatic Delaunay–Voronoi网格族 $\{K_h\}_{h>0}$ 上，对不可压缩Navier–Stokes/Euler系统建立如下四重收敛模式，统称为Main Theorem 1.1：

(I) 光滑解正则性：设 $u \in C^1([0,T]; W^{r_{\rm rec},\infty}(\Omega))$ 为参考光滑解，$v_h$ 为其对应的半离散DEC解，则

其中 $\beta_3 = 0$, $\beta_2 = 1$，且 $r_\star = 1$（一般网格）或 $r_\star = 2$（满足Property 2.2的质心邻近网格）。该速率关于 $\nu \ge 0$ 一致成立（Theorems 5.7 & 5.13）。

(II) Leray–Hopf弱解正则性：对任意 $L^2$ 初值，Whitney重构序列 $W_h v_h$ 在 $L^2(0,T; L^2(\Omega))$ 中存在子列收敛至Navier–Stokes方程的Leray–Hopf弱解 $u$（Theorem 5.22）。若存在强解，则弱–强唯一性定理（Theorem 5.24）恢复(I)中的收敛速率。

(III) 无粘测度值解正则性：当 $\nu = 0$ 时，$W_h v_h$ 存在子列在测度值意义下收敛至保守测度值Euler解 $(w, \sigma)$（Theorem 5.28），即对任意测试函数 $\phi \in C_c^\infty((0,T) \times \Omega; \mathbb{R}^d)$，

若进一步假设 $W_h v_h$ 在 $C^{0,\alpha}$ 中具有一致界且 $\alpha > 1/3$，则浓度缺陷 $\sigma$ 消失（Theorem 5.35）；若 $u \in L^\infty(0,T; W^{1,\infty}(\Omega))$ 是Lipschitz弱Euler解，则全序列收敛且 $\sigma = 0$（Theorem 5.32）。

(IV) 耗散解排除正则性：对任意 $\alpha \in (0,1]$，不存在任何子列 $W_h v_h$ 收敛至能量耗散的Euler解（即满足 $\frac{d}{dt} \int_\Omega \frac{|w|^2}{2} dx < 0$ 的解），无论网格尺寸 $h$、时间步长或初值如何（Proposition 5.33）。此排除由离散能量恒等式 $\frac{d}{dt} E_{\rm kin}^h = 0$（当 $\nu=0$）严格保证，且对Galerkin或有限体积法不成立。

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

本文核心障碍在于：如何在非线性、非均匀、非结构化网格上，同时实现（a）组合精确的拓扑算子（外微分 $d$）、（b）几何精确的度量算子（Hodge星 $\star$）、（c）代数精确的物理算子（内积 $i_v$ 与Lie导数 $\mathcal{L}_v$），三者缺一则无法导出全局守恒律。作者的真正创新点并非引入新算子，而是通过几何驱动的算子重构与误差源的严格隔离，使三者协同生效。

首先，利用Delaunay–Voronoi网格的正交性（Primal edge $e_k \perp$ Dual face $f_k^*$），将Hodge星 $M_k$ 严格对角化为 $ (M_k)_{ii} = |\sigma_i^*| / |\sigma_i| $，从而将全部近似误差集中于单一算子 $M_k$，而外微分 $D_k, \tilde{D}_k$ 作为纯组合coboundary算子完全精确（Section 2.1）。此设计使de Rham插值 $R_h$ 与离散外微分严格交换：$\tilde{D}_1 R_h = R_h d$（Lemma 2.8），消除了FEM中因质量矩阵稠密性导致的误差耦合。

其次，为实现Lamb矢量的精确反对称性 $\langle v_h, I_v(\tilde{D}_1 v_h) \rangle_1 = 0$，作者构造了基于extrusion的离散内积 $I_v$（Definition 2.13），其矩阵元 $ (I_v)_{jk} $ 显式依赖于对偶边 $e_j^*$ 与对偶面 $f_k^*$ 的几何权重，并被设计为满足 $I_v^T M_1 = -M_2 \tilde{D}_1$（Proposition 2.19）。此恒等式直接推出 $\langle v_h, I_v(\tilde{D}_1 v_h) \rangle_1 = -\langle \tilde{D}_1 v_h, \tilde{D}_1 v_h \rangle_2$，将非线性对流项的离散能量输运完全转化为耗散项（在粘性情形）或零（在无粘情形），这是所有守恒律（Theorems 3.4, 3.6, 3.8）的基石。

第三，在收敛性证明中，作者发展了一套针对DEC框架的稳定性分析工具：（i）Lamb反对称性消除立方自相互作用项 $\langle v_h, I_v(\tilde{D}_1 v_h) \rangle_1$ 的增长；（ii）利用DV正交性，将离散Hodge-Laplacian $\Delta_h = \delta_h \tilde{d}_1$ 精确识别为cell-centred有限体积Laplacian，从而获得尖锐的逐点重构界 $\|v_h - \bar{u}_j(v_h)\|_{L^\infty} \le C h^{r_{\rm rec}}$（Lemma 5.10(ii)）；（iii）极化反对称恒等式 $\langle v_h, I_u(w_h) \rangle_1 + \langle w_h, I_u(v_h) \rangle_1 = 0$ 与复形关系 $\tilde{D}_1 \tilde{D}_0 = 0$ 联用，彻底消除Leray投影余项中可能降低收敛阶的双线性项（Proof of Theorem 5.7）。这些工具共同构成‘一致性–稳定性’范式，使光滑收敛率不依赖于inf-sup条件，且关于 $\nu$ 一致。

4. 与之前工作的比较

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

1. 守恒律的严格性与选择能力：此前所有结构保持格式（包括FEEC[1,2]、EMAC[12]、Mullen–Pavlov[43,44]）均未证明其离散守恒律能作为连续极限的选择原理。FEEC仅处理线性Hodge-Laplace，未涉及非线性对流；EMAC的守恒依赖于特定形式的非线性混合，其离散能量恒等式含 $\mathcal{O}(h)$ 误差项，无法排除耗散解；Mullen–Pavlov虽在单纯形网格上实现精确守恒，但未建立任何收敛理论。本文首次证明：精确的离散能量恒等式 $\frac{d}{dt} E_{\rm kin}^h = 0$（$\nu=0$）在取极限时强制 $\sigma \equiv 0$，从而无条件排除所有能量耗散Euler解（Proposition 5.33），此为Galerkin/有限体积法绝对不具备的结构性优势。

2. 收敛范围的完整性：Guermond[25,26]对广义格式证明了收敛至CKN意义下的适定弱解，但其框架未利用守恒结构，且未覆盖Onsager阈值。本文在相同弱假设下，不仅得到Leray–Hopf解（II），更前延至测度值解（III）并后延至光滑解（I），形成首套覆盖 $C^{0,\alpha}$ 全谱（$0<\alpha\le1$）的统一理论。特别地，对 $\alpha > 1/3$ 的一致 $C^{0,\alpha}$ 界假设下缺陷消失（Theorem 5.35），直接对接Onsager正向证明[14]，而此前无任何离散格式能刻画此阈值。

3. 网格适用性与实用性：FEEC与多数DEC工作限于单纯形网格，而本文采用prismatic Delaunay–Voronoi网格（如ICON[34]、MPAS[47]所用），其正交性使Hodge星对角化，误差集中且可控（Lemma 2.10）。Property 2.2（质心邻近）将 $r_\star$ 从1提升至2，使收敛率从一阶升至二阶，此改进在球面大气模型中具直接工程价值，而FEM需复杂inf-sup稳定对才能达二阶。

4. 代价与局限：本文的严格守恒以放弃局部守恒通量为代价——离散动量方程不显式含通量 $F_{ij}$，而是通过extrusion算子隐式定义，故无法直接用于激波捕捉。此外，压力未作为演化变量求解，而是后处理自Bernoulli函数（Section 5），其收敛性弱于速度（仅 $L^{q_d}_t L^2_x$，$q_2=4/3$, $q_3=8/7$），而FEM可自然获得 $L^2_t L^2_x$ 压力收敛。最后，全文未处理边界条件下的Onsager阈值（仅在Section 6提及Dirichlet情形的扩展），而无界域或周期边界下的完整Onsager负向证明（凸积分构造）尚未离散化。

9. 超临界广义SQG方程在临界Sobolev空间中的渐近研究

Asymptotic study of supercritical generalized SQG equations in critical Sobolev spaces

作者: Anuj Kumar

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13176v1

摘要翻译: 我们研究超临界广义SQG（gSQG）方程在完全非线性区域中正则解的长时间行为。更准确地说，在初始数据关于临界Sobolev范数 $||{θ_0}||_{H^{1+β-2α}}$ 充分小的假设下，我们证明了唯一全局解的存在性；该解满足能量不等式 (1.3)，且其临界范数 $||{θ(t)}||_{H^{1+β-2α}}$ 随时间趋于无穷而衰减至 0。

Comment: 18 pages

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

1. 研究对象和背景

本文研究二维欧几里得空间 $\mathbb{R}^2$ 上超临界广义表面准地转（generalized surface quasi-geostrophic, gSQG）方程的长期行为，其具体形式为

其中 $\theta = \theta(x,t): \mathbb{R}^2 \times [0,\infty) \to \mathbb{R}$ 是演化标量场（如位势温度或浮力），$\Lambda = (-\Delta)^{1/2}$ 为Riesz位势算子，$\alpha \in (0,1)$ 为耗散阶数，$\beta \in (0,2)$ 为构型律参数。速度场 $u_\theta$ 由 $\theta$ 通过奇异积分算子非局部地确定，其傅里叶符号为 $\widehat{u_\theta}(\xi) = i \, \mathrm{sgn}(\xi_1,\xi_2)^\perp |\xi|^{\beta-2} \widehat{\theta}(\xi)$，其中 $\mathrm{sgn}(\xi_1,\xi_2)^\perp = (-\xi_2, \xi_1)/|\xi|$。

该方程族在 $\beta=0$ 时退化为二维欧拉方程的涡度形式，$\beta=1$ 时即经典SQG方程，而 $\beta \in (1,2)$ 则对应更强奇异性（比SQG更奇异）的构型律。物理上，SQG模型描述强分层、快速旋转流体中温度/浮力的演化，其数学重要性在于：它既是二维方程，又展现出类似三维欧拉方程的涡拉伸机制与湍流特征；同时，它作为“玩具模型”深刻揭示了非线性输运与分数阶耗散之间的微妙平衡——这使得其全局正则性问题成为偏微分方程分析领域的核心开放问题之一。

按耗散与非线性导数损失的尺度补偿关系，(1.1)被划分为三类：亚临界（$\beta < 2\alpha$）、临界（$\beta = 2\alpha$）与超临界（$\beta > 2\alpha$）。在亚临界与临界情形下，全局正则性已由多项工作确立（如[12,13,27,10,16,24]）；而在超临界情形，即使对小初值，全局解的存在性也极为困难。尤其当 $\beta \in [2\alpha+1, 2)$ 时，方程呈现“完全非线性”（fully nonlinear）结构：非线性项 $u_\theta \cdot \nabla \theta$ 中的速度 $u_\theta$ 本身含 $\Lambda^{\beta-2}$，与梯度 $\nabla \theta$ 结合后整体导数阶达 $\beta-1$，远超耗散项 $(-\Delta)^\alpha$ 提供的 $2\alpha$ 阶正则增益。此时，标准能量方法失效，因非线性项无法被控制在关键空间 $\dot{H}^{1+\beta-2\alpha}$ 中——该空间正是方程 (1.1) 的尺度不变临界空间：若 $\theta(t,x)$ 是解，则尺度变换 $\theta_\lambda(t,x) = \lambda^{2\alpha-\beta} \theta(\lambda^{2\alpha} t, \lambda x)$ 仍为解，且 $\|\theta_\lambda(t)\|_{\dot{H}^{1+\beta-2\alpha}} = \|\theta(\lambda^{2\alpha}t)\|_{\dot{H}^{1+\beta-2\alpha}}$，故 $\dot{H}^{1+\beta-2\alpha}$ 在尺度变换下不变。

已有工作表明，在 $\beta \in (1,2\alpha+1)$（拟线性区）和 $\beta = 1$（SQG）情形下，小初值可保证全局解存在并满足临界范数衰减（如[26,6,5,4,3,1]）；在 $\beta \in (2\alpha, \alpha+1)$ 且初值属于 Besov 空间 $\dot{B}^{\beta-2\alpha}_{\infty,\infty}$ 小范数时，亦有 $\|\theta(t)\|_{H^s} \to 0$（[31]）。然而，对最奇异的完全非线性区 $\beta \in [2\alpha+1, 2)$，此前尚无关于临界 Sobolev 范数 $\|\theta(t)\|_{H^{1+\beta-2\alpha}}$ 全局存在性与渐近衰减的严格结果。本文旨在填补这一空白：在 $\alpha \in (0,1/2)$、$\beta \in [2\alpha+1,2)$ 的完全非线性超临界设定下，建立小初值在临界空间 $H^{1+\beta-2\alpha}$ 中的全局适定性，并证明该临界范数随时间趋于零。

2. 主要定理或主要结果

定理 1.2 给出三项严格结论，全部在假设 $\alpha \in (0,1/2)$、$\beta \in [2\alpha+1,2)$、$\theta_0 \in H^{1+\beta-2\alpha}(\mathbb{R}^2)$ 下成立。设 $\theta \in C([0,T^*); H^{1+\beta-2\alpha})$ 为定理 1.1 所给局部解的最大存在解，则：

i) 爆破准则：若最大存在时间 $T^* < \infty$，则必有

该准则表明，$H^{1+\beta-2\alpha}$ 范数本身的有界性不足以防止爆破，而耗散项 $\Lambda^\alpha \theta$ 在该临界空间中的 $L^2_t$ 可积性是全局存在的充要条件。

ii) 全局存在性与能量不等式：存在仅依赖于 $\alpha, \beta$ 的常数 $C(\alpha,\beta) > 0$，使得当 $\|\theta_0\|_{H^{1+\beta-2\alpha}} < C(\alpha,\beta)$ 时，存在唯一全局解 $\theta$ 满足

且对任意 $t \geq 0$，能量不等式成立：

注意此处耗散空间为 $\dot{H}^{1+\beta-\alpha}$，其导数阶 $1+\beta-\alpha$ 严格大于临界空间导数阶 $1+\beta-2\alpha$（因 $\alpha > 0$），体现了耗散的正则化效应。

iii) 临界范数衰减：设 $\theta$ 为 ii) 中所得全局解，则

该衰减是强收敛意义下的，即在 $H^{1+\beta-2\alpha}$ 范数下收敛至零，而非弱收敛。证明中未给出衰减速率，但由能量不等式 (1.3) 可知 $\|\Lambda^\alpha \theta\|_{H^{1+\beta-2\alpha}} \in L^2_t$，结合插值与Gronwall型论证，隐含了某种（可能慢速的）代数衰减结构。

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

本文的核心技术障碍在于：在完全非线性区 $\beta \in [2\alpha+1,2)$，非线性项 $u_\theta \cdot \nabla \theta = \nabla^\perp \Lambda^{\beta-2} \theta \cdot \nabla \theta$ 的导数阶高达 $\beta-1 \geq 2\alpha$，而临界空间 $H^{1+\beta-2\alpha}$ 的嵌入性质不足以直接控制该高阶非线性。标准能量估计会遭遇 $\|\nabla^\perp \Lambda^{\beta-2} \theta\|_{L^\infty}$ 或 $\|\nabla \theta\|_{L^\infty}$ 类的失控项，而这些范数在临界空间中无定义。因此，必须绕过点态估计，发展一种能精细分离高低频相互作用、并精确量化奇异算子间交换误差的分析框架。

作者的真正创新点体现在引理 2.4 的构造与证明中。该引理断言：对 $\sigma \in (2-\alpha, 3-\alpha)$，有

此不等式的关键在于其右侧仅含临界范数 $\|\theta\|_{\dot{H}^{1+\beta-2\alpha}}$ 与一个“稍高阶”范数 $\|\theta\|_{\dot{H}^{\sigma+\alpha}}$ 的乘积，且系数 $C$ 不依赖于 $j$ 或 $t$。为证此式，作者将非线性项分解为三部分 $J_1 + J_2 + J_3$，每部分均通过引入合适的交换子（commutator）来剥离奇异算子作用顺序。例如，$J_1$ 被写为 $-\langle [\partial^\perp_\ell \Lambda^{\beta-2}, \partial_\ell \theta] \Lambda^\sigma_j \theta, \Lambda^\sigma_j \theta \rangle$，从而可直接应用引理 2.2（来自[23]）的尖锐交换子估计：对 $s=2-\beta, \rho_1=2\alpha, \rho_2=\alpha$，有

而 $J_2, J_3$ 则需更精细的处理，作者将其进一步拆解为 $J_{21}, J_{22}, J_{31}, J_{32}$，每一项均形如 $\langle [\Lambda^{s+\rho} \partial_\ell \Delta_i, g] f, h \rangle$，并应用引理 2.3（推广自[22]）：对 $|i-j| \leq k$，有

其中 $c_j \in \ell^2$ 且 $\|\{c_j\}\|_{\ell^2} \leq 1$。此估计的精妙之处在于，它允许将 $g$（此处为 $\theta$ 或其导数）的正则性要求降至 $\dot{H}^{\beta+1-2\alpha}$（即临界空间导数阶加1），而 $f,h$ 的正则性则通过 Littlewood-Paley 块 $\Delta_j$ 的频率局部化得以优化。最终，通过对所有 $j$ 求和并利用 Plancherel 和 $\ell^2$-正交性 (2.6)，将高频振荡抵消，得到所需的统一控制。

整个证明框架依赖于三个支柱：1) Littlewood-Paley 分解：将方程投影到 dyadic 频率块 $\Delta_j$，使非线性相互作用的频率几何（如低×高→高）得以显式追踪；2) 尖锐交换子估计（引理 2.2 与 2.3）：提供奇异算子与函数乘积间误差的最优 $\dot{H}^s$ 范数控制；3) 迭代能量法：在证明 i) 时，先假设 (1.2) 不成立，利用引理 2.5（引理 2.4 的特例，取 $\sigma = 1+\beta-2\alpha$）导出 $\|\theta(t)\|_{H^{1+\beta-2\alpha}}$ 在 $[0,T^*)$ 上一致有界；再结合 $L^2$ 强连续性（引理 2.6）与插值，推出 $\theta(t)$ 在 $t \uparrow T^*$ 时强收敛于某 $\tilde{\theta} \in H^{1+\beta-2\alpha}$，从而可将解延拓过 $T^*$，矛盾。ii) 与 iii) 的证明则基于对 $\gamma_n \uparrow 1+\beta-2\alpha$ 的一族逼近能量不等式取极限，并利用单调收敛与支配收敛定理。

4. 与之前工作的比较

本文与既有工作的核心推进体现在问题设定、技术工具与结论强度三个维度。

问题设定维度：此前关于超临界 gSQG 渐近性的成果主要集中于两个区域：(a) 拟线性区 $\beta \in (1,2\alpha+1)$，如 [26] 所证；(b) SQG 本体 $\beta = 1$，如 [6,5,4,3,1] 所证；(c) Besov 小初值情形 $\theta_0 \in \dot{B}^{\beta-2\alpha}_{\infty,\infty}$，如 [31] 所证。本文首次将结论拓展至最困难的完全非线性区 $\beta \in [2\alpha+1,2)$。该区域的标志性特征是 $u_\theta$ 的奇异性 $\Lambda^{\beta-2}$ 与 $\nabla \theta$ 的结合导致非线性项的整体导数阶 $\beta-1 \geq 2\alpha$，使得 [26] 中用于拟线性情形的“导数分裂”技巧（如将 $u_\theta$ 视为低阶系数）彻底失效。本文必须直面完全非线性结构，其技术难度远超前人。

技术工具维度：[26] 的核心工具是针对拟线性情形设计的交换子估计，其假设 $\beta < 2\alpha+1$ 以保证 $\Lambda^{\beta-2}$ 的奇异性可控；而本文引理 2.4 是专为 $\beta \geq 2\alpha+1$ 设计的全新估计，其证明中 $J_2,J_3$ 的分解与引理 2.3 的多参数应用（特别是对 $\nu$ 的灵活选取，如 (2.18) 中取 $\nu = 2-\beta+2\alpha$）是此前文献未见的。此外，[31] 依赖于 Besov 空间的极大函数刻画与 $L^p$ 插值，其结论 $\|\theta(t)\|_{H^s} \to 0$ 仅对 $s > 0$ 成立，且初值需在 $\dot{B}^{\beta-2\alpha}_{\infty,\infty}$ 中小，该空间与 $H^{1+\beta-2\alpha}$ 互不嵌入（前者侧重高频振荡，后者侧重整体平方可积性），故 [31] 的结果与本文不可比较。本文坚持在纯 Sobolev 框架内工作，结论更贴近经典正则性理论的语言。

结论强度维度：本文定理 1.2 ii) 给出了严格的小初值全局存在性，其小性阈值 $C(\alpha,\beta)$ 是明确的（尽管未给出显式公式），且解自动满足能量不等式 (1.3)，这是强解（strong solution）的标志。而 [26] 与 [6] 的结论虽也包含小初值全局性，但其证明中隐含的常数依赖于更高阶范数，本文则通过引理 2.4 实现了仅依赖临界范数的闭环估计。更重要的是，iii) 中的强范数衰减 (1.4) 是全新结果：此前工作或仅给出弱收敛（如 [26] 的注释），或仅在更低阶范数中衰减（如 [31] 的 $H^s$）。本文证明了临界空间自身范数的强衰减，这蕴含了系统在临界尺度上的彻底弛豫。代价在于：本文要求 $\alpha < 1/2$（而 [26] 允许 $\alpha \in (0,1)$），且 $\beta$ 的下界被锁定在 $2\alpha+1$，这反映了完全非线性区固有的技术限制；此外，衰减速率未被量化，而 [6] 对 SQG 给出了具体的 $t^{-1/2}$ 速率。从提供内容中无法确认本文是否可推广至周期域或有界区域，亦无法确认其是否蕴含解的解析性或 Gevrey 光滑性，这些是未来可拓展的方向。

10. 退化粘性Hamilton-Jacobi方程的一个新选择问题

A new selection problem for degenerate viscous Hamilton-Jacobi equations

作者: Qinbo Chen, Zhi-Xiang Zhu

类别: math.AP, math.DS

发布日期: 2026-05-13

arXiv ID: 2605.12996v1

摘要翻译: 我们研究具有凸Hamiltonian的退化粘性Hamilton–Jacobi方程的选择问题，其中逼近过程结合了一种非线性折扣逼近与一个小势能扰动。一个关键问题是：二者的联合作用如何影响相应遍历问题的粘性解的渐近选择。基于非线性伴随方法，我们证明了逼近解一致收敛于遍历问题的一个特定解，并推导出所选极限的表达式，该表达式以广义Mather测度和势能为变量。作为应用，我们证明该选择原理具有足够灵活性，可实现遍历问题的任意给定解，且具显式的收敛速率。

Comment: 无

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

1. 研究对象和背景

本文研究的是环面 $\mathbb{T}^n = \mathbb{R}^n / \mathbb{Z}^n$ 上退化粘性Hamilton-Jacobi方程的遍历问题及其逼近解的渐近选择机制。核心模型为如下退化二阶偏微分方程：

其中 $H : \mathbb{T}^n \times \mathbb{R}^n \to \mathbb{R}$ 是 $C^2$ 光滑、关于 $p$ 严格凸且超线性的Hamiltonian；扩散矩阵 $A(x) = (a_{ij}(x))_{n\times n} \in C^2(\mathbb{T}^n)$ 是半正定的（即 $A(x) \succeq 0$ 对所有 $x$ 成立），但允许在某些点或子集上秩亏缺（degenerate），例如 $A(x) \equiv 0$（一阶情形）或 $A(x) = a(x) I_n$ 且 $a(x) \ge 0$ 在某处为零；$c_H \in \mathbb{R}$ 是唯一确定的遍历常数（ergodic constant），其存在性与唯一性由 $H$ 和 $A$ 决定，参见 [42, 5]。

由于退化性，方程 (E) 的粘性解在加法常数意义下不唯一：即使固定 $c_H$，解集 $\mathcal{S}_H := \{ u \in C(\mathbb{T}^n) : u \text{ 是 } (E) \text{ 的粘性解} \}$ 通常是一个无限维仿射空间（模常数）。这一非唯一性引出一个根本性问题：给定某种自然逼近程序 $\{u_\lambda\}_{\lambda > 0}$，当 $\lambda \to 0^+$ 时，$u_\lambda$ 是否一致收敛？若收敛，极限 $u_0$ 是 $\mathcal{S}_H$ 中哪一个特定解？该问题称为“选择问题”（selection problem），其物理含义在于：不同逼近路径对应不同物理正则化（如不同时间尺度的折扣、不同强度的噪声、不同外力扰动），它们在零极限下“挑选”出遍历系统中具有特定动力学/变分意义的临界解。

本文考虑的逼近方程为：

该方程包含两个关键扰动项：(i) 非线性折扣项 $f(x, \lambda u_\lambda)$，满足 $\partial_r f > 0$ 且 $f(x,0)\equiv 0$，涵盖标准折扣 $f(x,r)=r$ 和空间变折扣 $f(x,r)=\sigma(x)r$（$\sigma>0$）；(ii) 小势能扰动 $\lambda V(x)$，其中 $V \in \mathrm{Lip}(\mathbb{T}^n)$。二者同时作用，且均以 $O(\lambda)$ 尺度 vanish，这构成一个全新的耦合选择机制。当 $V \equiv 0$ 时，(E$_\lambda$) 退化为经典 vanishing discount 问题；当 $A \equiv 0$ 时，则为一阶动力系统中的弱KAM选择问题。然而，已有工作表明：仅靠改变折扣结构（如 [12, Thm 1.5] 中的非线性折扣）无法实现对任意 $\bar{u}_0 \in \mathcal{S}_H$ 的精确选择。例如，对机械Hamiltonian $H(x,p) = \frac{1}{2} p^2 + \cos 4\pi x$ 在 $\mathbb{T} = \mathbb{R}/\mathbb{Z}$ 上，其Mather测度支撑于双曲平衡点 $(0,0)$ 和 $(\frac{1}{2}, 0)$，vanishing discount 极限恒选同一解，与折扣形式无关。因此，引入 $\lambda V(x)$ 并非技术冗余，而是突破选择刚性的本质手段：它将选择机制从单参数（折扣率）推广至函数空间参数（势能 $V$），从而有望覆盖整个解集 $\mathcal{S}_H$。本文旨在严格建立这一新机制的数学基础，核心挑战在于：(1) 退化扩散 $A(x)$ 破坏了标准椭圆估计的统一性；(2) 非线性折扣 $f(x,\lambda u)$ 与势能扰动 $\lambda V$ 的耦合导致传统线性对偶方法失效；(3) 粘性解缺乏光滑性，阻碍对广义Mather测度的直接积分操作；(4) 必须精确刻画 $V$ 的尺度效应——作者强调 $\lambda V$ 是临界尺度（Remark 1.3），若换为 $\lambda^\alpha V$，则 $\alpha > 1$ 时扰动无效，$\alpha < 1$ 时破坏收敛性。

2. 主要定理或主要结果

在假设 (H1)–(H3) 下（即 $H$ 严格凸超线性、$A(x) \succeq 0$ 且 $C^2$、$V \in \mathrm{Lip}$、$f \in C^1$ 且 $\partial_r f > 0$，并归一化 $f(x,0)\equiv 0$），本文证明了以下三个主定理：

定理 A（选择极限的刻画）：存在 $\lambda_0 > 0$，使得对每个 $\lambda \in (0, \lambda_0]$，方程 (E$_\lambda$) 存在唯一 Lipschitz 连续粘性解 $u_\lambda : \mathbb{T}^n \to \mathbb{R}$；且当 $\lambda \to 0^+$ 时，$u_\lambda$ 一致收敛于遍历问题 (E) 的某个特定解 $u_0 \in C(\mathbb{T}^n)$，其表达式为

其中 $\mathcal{E}$ 是所有满足如下积分约束的连续粘性解 $w$ 的集合：

此处 $\widetilde{\mathcal{M}}_g$ 是由非线性伴随方法构造的广义Mather测度集（定义见原文 (39) 式），其元素是 $\mathbb{T}^n \times \mathbb{R}^n$ 上的概率测度，满足动力学相容性与变分最小性。

定理 B（基于广义Mather测度的比较原理）：在 (H1)–(H2) 下，对任意正函数 $\sigma \in C(\mathbb{T}^n)$ 及 (E) 的任意两个连续粘性解 $u_1, u_2$，若

则必有 $u_1 \le u_2$ 在 $\mathbb{T}^n$ 上成立。特别地，取 $\sigma \equiv 1$，即得：$u_1, u_2$ 由其在全体广义Mather测度上的积分唯一确定，此为 Mitake–Tran [44] 结果的推广（假设更弱）。

定理 C（全解可实现性与显式收敛率）：在 (H1)–(H2) 及 $f$ 满足 (H3) 下，对遍历问题 (E) 的任一粘性解 $\bar{u}_0 \in C(\mathbb{T}^n)$，存在势能 $\bar{V}_0 \in \mathrm{Lip}(\mathbb{T}^n)$，使得对充分小的 $\lambda > 0$，扰动方程

有唯一 Lipschitz 解 $u_\lambda$，且满足

其中常数 $C > 0$ 与 $\lambda$ 无关。该 $O(\lambda)$ 收敛率是尖锐的（sharp），由 $\bar{V}_0$ 的特定构造保证。

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

本文的核心技术障碍在于：退化扩散 $A(x)$ 导致标准椭圆正则性理论（如 Schauder 或 $W^{2,p}$ 估计）完全失效；粘性解 $u_\lambda$ 仅为 Lipschitz 连续，无法直接对其应用分部积分或构造测度；而广义Mather测度 $\widetilde{\mathcal{M}}_g$ 的定义依赖于解的二阶导数信息，需在缺乏光滑性的前提下建立其存在性与变分性质。作者发展了一套精密的分析框架，其创新点分述如下：

1. 退化方程的统一Lipschitz估计（Theorem 2.1）：针对形如 $-\operatorname{tr}(A(x)D^2 u) + f(x, \delta u) + H(x, Du) + \delta V(x) = 0$ 的方程，作者采用 doubling variables 方法，但关键改进在于：利用 $A(x) = \Sigma(x)^T \Sigma(x)$ 的Lipschitz平方根分解（[47, Theorem 5.2.3]），并构造辅助矩阵 $G_t = \begin{pmatrix} t^2 \Sigma(x_0)^T \Sigma(x_0) & t \Sigma(x_0)^T \Sigma(y_0) \\ t \Sigma(y_0)^T \Sigma(x_0) & \Sigma(y_0)^T \Sigma(y_0) \end{pmatrix}$。通过将Crandall-Ishii引理所得矩阵不等式与 $G_t$ 做迹运算，成功将退化项 $\operatorname{tr}(A(x_0) X_\varepsilon) - \operatorname{tr}(A(y_0) Y_\varepsilon)$ 控制在 $O(|x_0-y_0|)$ 尺度内，最终导出显式Lipschitz常数 $K = e^{2\gamma_1}(\gamma_2 + 2\|f\|_{C^1(U)}) + \delta(e^{2\gamma_1} \|V\|_{L^\infty} + \mathrm{Lip}(V)) + C$，其中 $C$ 仅依赖于 $\gamma_1$ 和 $\Lambda_0 = (\mathrm{Lip}(\Sigma))^2 + \sup_x \|\Sigma(x)\|_F^2$。此估计是后续所有分析（存在性、紧性、渐近分析）的基石。

2. 广义Mather测度的非线性伴随构造（Section 3）：为克服 $u_\lambda$ 不光滑的困难，作者引入椭圆正则化 $u_\lambda^\eta$（添加 $\eta^2 \Delta u$ 项），使其成为 $C^\infty$ 解。对正则化方程 (E$_\lambda^\eta$)，定义其伴随方程（adjoint equation）：

其中 $\mathbf{b}_\lambda^\eta$ 是由 $A(x)$ 和 $f$ 诱导的漂移场。解 $m_\lambda^\eta$ 是概率密度，其弱极限（当 $\eta \to 0$）定义广义Mather测度 $\tilde{\mu}_\lambda$。关键创新在于：作者证明了 $\{\tilde{\mu}_\lambda\}_{\lambda > 0}$ 在 $\mathcal{P}(\mathbb{T}^n \times \mathbb{R}^n)$ 中紧，并提取极限 $\tilde{\mu} \in \widetilde{\mathcal{M}}_g$，该测度自动满足 $\int \phi(x) \, d\tilde{\mu} = \lim_{\lambda \to 0} \int \phi(x) \, m_\lambda^\eta(x) \, dx$ 对所有 $\phi \in C(\mathbb{T}^n)$ 成立，从而将解的宏观行为编码于测度中。

3. 针对粘性解的弱光滑逼近（Theorem 4.1）：为处理 $u_\lambda$ 的不可微性，作者未采用文献 [43] 中依赖于标量扩散 $A(x)=a(x)I_n$ 的交换引理（commutation lemma），而是设计了一种新的逼近方案：先对 $u_\lambda$ 应用 Lasry–Lions 正则化（得到 $u_\lambda^{\varepsilon}$），再进行额外的 mollification（磨光）。该组合确保 $u_\lambda^{\varepsilon, \delta} \in C^\infty$，且其梯度和Hessian在 $L^1$ 意义下收敛于 $Du_\lambda$ 和 $D^2 u_\lambda$ 的Radon–Nikodym导数，从而允许在极限过程中合法地执行 $\int \psi(x, Du_\lambda^{\varepsilon, \delta}) \, d\tilde{\mu}_\lambda \to \int \psi(x, Du_\lambda) \, d\tilde{\mu}_\lambda$。

4. 渐近约束的双重估计（Lemmas 5.1 & 5.2）：这是选择机制的核心。作者证明：对任意 $\tilde{\mu} \in \widetilde{\mathcal{M}}_g$，极限解 $u_0$ 必须满足

且 $u_0$ 是满足此不等式的最大解。该估计源于对正则化解 $u_\lambda^\eta$ 与伴随解 $m_\lambda^\eta$ 的乘积做分部积分，并利用 $H$ 的凸性和 $f$ 的单调性控制误差项，最终令 $\eta, \lambda \to 0$。

4. 与之前工作的比较

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

1. 选择机制的灵活性维度：此前所有工作（[17, 43, 31, 32, 53, 12]）均局限于单一扰动源（discount 或 viscosity），其选择结果被证明是固定的：vanishing discount 总收敛到以经典Mather测度为权重的“最平坦”解；vanishing viscosity 则倾向于最小化扩散耗散的解。本文首次引入双扰动耦合 $f(x, \lambda u) + \lambda V(x)$，将选择参数从标量 $\lambda$ 扩展至函数空间 $V \in \mathrm{Lip}(\mathbb{T}^n)$。定理 C 表明，对任意预设解 $\bar{u}_0$，总存在 $V$ 实现 $u_\lambda \to \bar{u}_0$，这在以往任何单一扰动框架下均不可能（如 [12] 明确指出其选择与折扣结构无关）。代价是：必须精确调控 $V$ 的大小（$O(1)$ 量级），且收敛率 $O(\lambda)$ 依赖于 $V$ 的特定构造（Remark 1.4）。

2. 退化性处理的技术维度：Mitake–Tran [43] 仅处理标量扩散 $A(x) = a(x) I_n$，其证明依赖于 $a(x)$ 的标量性以应用交换引理；Ishii–Mitake–Tran [31, 32] 的duality方法要求 $A(x)$ 为正定（uniformly elliptic）。本文通过 $A(x) = \Sigma(x)^T \Sigma(x)$ 的Lipschitz平方根分解与矩阵迹估计（Theorem 2.1），首次将选择理论推广至一般半正定矩阵 $A(x)$，覆盖了从完全退化 ($A\equiv 0$) 到部分退化 ($\mathrm{rank}(A(x)) < n$) 的全谱情形。此推广的代价是：Lipschitz常数 $K$ 显式依赖于 $\mathrm{Lip}(\Sigma)$ 和 $\|\Sigma\|_F$，而非仅 $\|A\|_{C^2}$。

3. 广义测度的构造维度：经典Mather测度 [23] 定义于一阶系统，要求 $A\equiv 0$；随机Mather测度 [26, 27] 仅适用于 $A(x) = a I_n$（常数扩散）。本文定义的 $\widetilde{\mathcal{M}}_g$ 是首个适用于一般 $A(x)$ 的广义测度，其构造完全基于非线性伴随方法（而非变分公式），且其元素 $\tilde{\mu}$ 满足 $\int \phi(x) \, d\tilde{\mu} = \lim \int \phi \, m_\lambda^\eta \, dx$，天然兼容退化性。与 [44] 相比，本文的测度集 $\widetilde{\mathcal{M}}_g$ 更大（包含更多约束），故定理 B 的比较原理在更弱的 $H$ 假设下成立。

4. 收敛率的显式性维度：以往工作（如 [41, 43]）仅证明收敛性，未给出速率。本文定理 C 首次为退化粘性HJ方程建立了尖锐的 $O(\lambda)$ 一致收敛率，且明确指出该速率依赖于 $V$ 的精心选取。从提供内容中无法确认是否存在对任意 $V$ 的通用收敛率估计。

11. 从多项式稳定性到部分耗散系统中的周期适定性

From Polynomial Stability to Periodic Well-posedness in Partially Dissipative Systems

作者: Giovanni P. Galdi, Boris Muha, Justin T. Webster

类别: math.AP, math.DS

发布日期: 2026-05-13

arXiv ID: 2605.12892v1

摘要翻译: 在应用中，研究复杂系统在时间周期载荷下的共振（及适定性）问题具有广泛意义。Galdi 等人（2014）将无外力 Cauchy 问题解的渐近稳定性与时间周期强迫问题的可解性联系起来。解半群的一致稳定性保证了对自然 mild 强迫类中所有外力均成立的周期适定性，而强稳定性仅能推出存在一个稠密的外力集合，使得其上可排除共振。我们处理一种介于二者之间的中间情形：多项式（亦称：有理型或半一致）稳定性。借助 Hilbert 空间中的 Fourier 分解，我们证明：半群的多项式稳定性导出周期适定性成立的稠密强迫集的一个显式刻画。更精确地，预解估计直接转化为为保证适定性而在外力上所要求的时间导数损失；即若 $\| (i\tau - A)^{-1} \|_{\mathcal{L}(H)} \lesssim |\tau|^\alpha$ 对 $|\tau| \gg 1$ 成立，则需外力 $f \in L^2_{\mathrm{per}}(0,2\pi; D(A^{-\beta}))$，其中 $\beta > \alpha/2$。本结果受部分耗散模型驱动——包括著名的理想化流体-结构相互作用的热-波耦合问题，以及某些热弹性、粘弹性及弱阻尼双曲系统——在这些模型中，多项式衰减是自然的衰减机制。

Comment: 无

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

1. 研究对象和背景

本文研究的核心对象是一类由希尔伯特空间 $ H $ 上的线性无界算子 $ A : D(A) \subset H \to H $ 生成的 $ C_0 $-半群 $ S(t) = e^{At} $ 所刻画的部分耗散系统（Partially Dissipative Systems, PDS），其演化由抽象发展方程

描述。该方程建模了大量物理耦合系统：典型例子包括热-波耦合系统（如流体-结构相互作用的理想化模型）、热弹性系统、黏弹性系统、弱阻尼双曲系统等。这些系统的共同特征是耗散仅作用于系统的一部分变量或子区域，例如在热-波系统中，热方程部分提供强耗散（$ L^2 $-能量指数衰减），而波方程部分本身无耗散（仅具守恒能量），二者通过界面条件强耦合。因此，整体系统的长期行为既非均匀稳定（uniformly stable），亦非仅强稳定（strongly stable），而是处于中间状态——即存在多项式衰减率：

这一衰减率由Borichev–Tomilov定理（Proposition 3）等价于虚轴上的多项式型预解估计：

此类估计在PDE中自然出现：例如，在界面耦合热-波系统中，若几何满足经典几何控制条件（GCC），则可证得 $ \alpha = 2 $（即 $ \|S(t)A^{-1}\| \lesssim t^{-1/2} $）；在更优分析下（如[14]），甚至可达 $ \alpha = 1 $。问题的关键在于：当对系统施加时间周期外力 $ F(t+T) = F(t) $ 时，系统是否具有唯一的 $ T $-周期解？此即周期适定性（periodic well-posedness）问题，其本质是排除共振（resonance）——即避免解随时间无界增长。共振的发生直接关联于预解集 $ \rho(A) $ 在虚轴 $ i\mathbb{R} $ 上的行为：若 $ in\omega \in \sigma(A) $ 对某个 $ n \in \mathbb{Z} $ 成立（$ \omega = 2\pi/T $），则第 $ n $ 阶傅里叶模态无法求逆，导致解不存在或不唯一。已有理论框架（Galdi et al. 2014 [1]）建立了稳定性与周期适定性的宏观联系：若 $ S(t) $ 一致稳定，则对任意 $ F \in L^1(0,T;H) $，方程 (1) 存在唯一 $ T $-周期解（无条件适定）；若仅强稳定，则仅能保证存在一个抽象稠密集 $ Q \subset L^1(0,T;H) $，使得对所有 $ F \in Q $，解存在唯一。然而，强稳定情形下 $ Q $ 的结构完全未知，无法用于实际建模或控制设计。本文旨在填补这一关键空白：在多项式稳定这一PDS中最自然、最典型的中间稳定性范畴内，给出 $ Q $ 的显式、定量刻画，并揭示其与预解增长指数 $ \alpha $ 的精确对应关系。这不仅是理论上的补全，更是为工程中预测和规避共振提供可计算的正则性门槛。

2. 主要定理或主要结果

本文核心结果为Theorem 4：设 $ A : D(A) \subset H \to H $ 生成有界 $ C_0 $-半群 $ S(t) $，且 $ S(t) $ 是阶为 $ 1/\alpha $ 的多项式稳定半群（即满足定义1(iii)及Proposition 3等价条件）。则对任意 $ T > 0 $ 及任意整数 $ m \geq 1 $，方程 (1) 在强迫项 $ F \in H^{m+\alpha}_\#(0,T; H) $ 下是周期适定的，即存在唯一 $ T $-周期温和解 $ U \in C_\#([0,T]; H) $，且满足数据到解的连续映射估计

其中 $ H^{s}_\#(0,T; H) $ 表示 $ T $-周期分布意义下时间导数阶数为 $ s $ 的索博列夫空间，其范数由傅里叶系数定义：

定理要求 $ m \geq 1 $ 是为了确保解 $ U $ 连续（$ H^1_\# \hookrightarrow C_\# $），从而明确定义初始/终值条件 $ U(0) = U(T) $。作为推论，在热-波耦合模型中，若已知 $ \alpha = 2 $（标准几何控制下），则取 $ m = 1 $ 得：

即有限能量周期解的存在性要求强迫项至少属于三阶时间索博列夫空间；若采用分数阶嵌入，亦可取 $ F \in H^{5/2+}_\#(0,T; H) \implies U \in H^{1/2+}_\#(0,T; H) $。若模型达到最优 $ \alpha = 1 $（如[14]），则结论强化为 $ F \in H^{m+1}_\#(0,T; H) \implies U \in H^{m}_\#(0,T; H) $。

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

证明的核心障碍在于：多项式稳定性（即预解在 $ i\mathbb{R} $ 上的多项式增长）如何转化为周期解存在性所需的定量正则性提升？作者摒弃了传统能量方法或紧性论证，转而采用希尔伯特空间中的傅里叶模态分解这一纯粹谱分析框架。其真正创新点在于将Borichev–Tomilov定理（Proposition 3）的等价性从时间域直接“采样”到周期问题的离散频谱 $ \{in\omega : n \in \mathbb{Z}\} $ 上：由 (3) 立即导出关键的离散预解估计

其中 $ M_T $ 依赖于 $ T $ 但与 $ n $ 无关。此估计是全文的基石，它将抽象的谱性质 $ \alpha $ 直接翻译为傅里叶系数层面的乘子增长阶。随后，作者构造模态解 $ U_n = (in\omega I - A)^{-1} F_n $，并利用 (5) 对 $ U_n $ 的索博列夫范数进行逐项控制：

对 $ n $ 求和即得解 $ U(t) = \sum_n U_n e^{in\omega t} $ 属于 $ H^m_\# $，且满足估计 (4)。为验证 $ U $ 确为温和解，作者引入傅里叶部分和 $ F^N $ 和 $ U^N $，并利用 $ m \geq 1 $ 保证 $ U^N \to U $ 在 $ C([0,T]; H) $ 中一致收敛，从而可在变分参数公式

中取极限，严格确立 $ U $ 的温和解身份及周期性。整个框架的精妙之处在于：它完全绕开了PDE的具体结构（如椭圆性、特征线），仅依赖半群生成元 $ A $ 的谱性质，因而具有高度的抽象普适性；同时，它将预解增长这一“谱障碍”直接量化为解空间所需的“正则性代价”，揭示了共振规避的本质是对强迫项施加可计算的时间光滑性要求。

4. 与之前工作的比较

本文与既有工作存在三个维度的明确推进。首先，相较于Galdi et al. (2014) [1] 的Proposition 2：(i) 其一致稳定性结论给出无条件适定性（$ F \in L^1 $），而本文在更弱的多项式稳定性下，虽不能覆盖全部 $ L^1 $，却给出了显式、稠密、可计算的子集 $ H^{m+\alpha}_\# $，其结构由 $ \alpha $ 完全决定；(ii) 其强稳定性结论仅断言存在某个抽象稠密集 $ Q $，而本文将 $ Q $ 显式构造为 $ H^{m+\alpha}_\# $，并提供了定量的导数损失律（$ \alpha $ 阶时间导数损失），这是根本性突破。其次，相较于Borichev–Tomilov (2010) [4]：该文建立了多项式衰减与预解增长的等价性，但未涉及周期问题；本文则将此等价性直接应用于周期框架，实现了从渐近行为到稳态响应的桥梁构建。第三，相较于Mosný–Muha–Schwarzacher–Webster (2024) [22] 的直接时域构造：(i) [22] 针对热-波系统，通过观测性估计得到 $ F \in H^4_\# $ 的适定性，而本文由 $ \alpha = 2 $ 推出 $ F \in H^3_\# $，表明本文结果在正则性门槛上更优，且不依赖于[22]中可能存在的非最优常数（如迹嵌入损失）；(ii) 几何假设上，[22] 允许更广的几何（因其不依赖GCC），而本文需GCC以保证 $ \alpha $ 存在，故适用范围更窄但正则性要求更低；(iii) [22] 的估计是分量式的（不同方程分量有不同导数损失），而本文是全局式的（统一施加 $ H^{m+\alpha}_\# $），二者视角互补。本文的代价在于：它要求 $ A $ 生成有界半群且 $ 0 \in \rho(A) $（即 $ A^{-1} $ 有界），这在某些无界耗散模型中可能不成立；此外，其结论限于有限能量解（$ H^m_\# $），对更弱解（如[22]中的超弱解）未作讨论。从提供内容中无法确认本文结果是否可推广至非自伴或非耗散主导的算子 $ A $。

12. 无散度被动向量场中的混合与小尺度生成

Mixing and Small-Scale Formation in a Passive Divergence-Free Vector Field

作者: Anuj Kumar, Franziska Weber

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.12883v1

摘要翻译: 我们研究一个由另一无散度向量场 $U$ 输运的无散度被动向量场 $u$ 的混合问题，其中 $u$ 满足方程 $ \partial_t u + (U \cdot \nabla) u + \nabla p = 0.$ 近年来，关于标量情形下的最优混合问题受到广泛关注，此时对输运速度场施加了 Sobolev 范数约束。然而，在本文所考虑的向量情形中，压力项引入了实质性困难，因为标量情形下可用的简单拉格朗日视角已不再适用。

Comment: 无

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

1. 研究对象和背景

本文研究定义在 $ d $-维环面 $ \mathbb{T}^d $ 上的无散度向量场 $ u : \mathbb{R}_+ \times \mathbb{T}^d \to \mathbb{R}^d $ 在另一给定无散度背景速度场 $ U : \mathbb{R}_+ \times \mathbb{T}^d \to \mathbb{R}^d $ 作用下的线性输运演化，其动力学由如下系统刻画：

其中压力项 $ \nabla p $ 由 $ \mathrm{div}\, u = 0 $ 和 $ \mathrm{div}\, U = 0 $ 确定，即 $ p $ 是满足 $ -\Delta p = \mathrm{div}\, ((U \cdot \nabla) u) $ 的唯一零均值解（在 $ \mathbb{T}^d $ 上）。该模型本质上是线性化欧拉方程中将流体速度 $ u $ 与驱动场 $ U $ 解耦后的最简结构：当 $ u = U $ 时，(1.1) 恢复为不可压欧拉方程；而当 $ u $ 为标量 $ \theta $ 时，则退化为经典被动标量输运方程 $ \partial_t \theta + (U \cdot \nabla) \theta = 0 $。研究对象 $ u $ 是一个“被动”向量场——它不反作用于 $ U $，但其几何结构（如旋度、剪切、拉伸）受 $ U $ 的输运和压力投影共同调制。

该问题的核心动机源于湍流理论中两个根本性开放问题：(i) 小尺度（即高频 Fourier 模态）如何在无粘/弱粘流体中自发生成？(ii) 压力项在尺度级联与能量转移中扮演何种角色？在标量情形下，混合速率可通过 $ H^{-\alpha} $ 范数衰减量化，且已有深刻结论：若 $ U \in W^{1,q}(\mathbb{T}^d) $，$ q > 1 $，则对任意初始 $ \theta_0 \in L^2 $，存在指数下界 $ \|\theta(t)\|_{H^{-1}} \geq C e^{-Ct} $（Crippa–De Lellis），且该下界是紧的（Alberti–Crippa–Mazzucato 构造达优场）。然而，向量情形因压力项破坏了 Lagrangian 不变量（如标量沿迹线守恒），导致经典方法失效：$ u $ 不再沿 $ U $ 的轨迹平移，而是经历非局部投影修正 $ u \mapsto \mathbb{P}((U \cdot \nabla) u) $，其中 $ \mathbb{P} = I - \nabla \Delta^{-1} \mathrm{div} $ 为 Leray 投影。因此，$ u $ 的 Fourier 系数演化为耦合非局部 ODE 系统：

其中 $ \mathbb{P}_k v = v - (v \cdot k) k / |k|^2 $。该系统中，压力通过正交投影强制 $ \hat{u}_k \perp k $，从而将 $ u $ 的能量在不同波数间重分配，而非简单拉伸。现有理论空白在于：缺乏对 $ H^{-\alpha} $ 范数衰减率的定量下界估计；未厘清 $ W^{1,q} $ 类 $ U $ 对 $ u $ 小尺度生成能力的临界依赖；更关键的是，尚无框架能分离并量化压力项对混合的独立贡献——这正是本文试图填补的空白：将 (1.1) 视为一个“压力隔离平台”，在剥离非线性惯性项（如 $ (u \cdot \nabla) u $）的前提下，系统考察 $ \nabla p $ 如何与输运项协同或竞争以驱动 $ u $ 的频谱展宽。

2. 主要定理或主要结果

本文在环面 $ \mathbb{T}^d $ 上建立以下严格结果：

1. 解的存在性与正则性（引理 2.1–2.4）：对任意 $ u_0 \in L^2_{\mathrm{div}}(\mathbb{T}^d) $，若 $ U \in L^\infty(\mathbb{R}_+; W^{1,q}(\mathbb{T}^d)) $ 且 $ \mathrm{div}\, U = 0 $，则 (1.1) 存在弱解 $ u \in L^\infty(\mathbb{R}_+; L^2_{\mathrm{div}}(\mathbb{T}^d)) $，且满足能量不等式 $ \|u(t)\|_{L^2} \leq \|u_0\|_{L^2} $。进一步，时间导数正则性为：

   $$ \partial_t u \in \begin{cases} L^\infty(\mathbb{R}_+; H^{-1}(\mathbb{T}^d)), & q > d, \\ L^\infty(\mathbb{R}_+; W^{-1, \frac{2dq}{dq + 2d - 2q}}(\mathbb{T}^d)), & q < d. \end{cases} $$

   若额外假设 $ U \in L^1(\mathbb{R}_+; W^{1,\infty}(\mathbb{T}^d)) $，则解唯一（引理 2.5），且对初值稳定：$ \|u(t) - v(t)\|_{L^2}^2 = \|u_0 - v_0\|_{L^2}^2 $。更高正则性要求更强条件：若 $ u_0 \in H^1_{\mathrm{div}} $，则 $ \|\nabla u(t)\|_{L^2} \leq \|\nabla u_0\|_{L^2} \exp\left( \int_0^t \|\nabla U(s)\|_{L^\infty} \, ds \right) $（引理 2.3）；若 $ u_0 \in H^2 \cap H^1_{\mathrm{div}} $ 且 $ U \in L^1(\mathbb{R}_+; W^{2,3}(\mathbb{T}^d)) $，则 $ \|\nabla^2 u(t)\|_{L^2} \leq \|\nabla^2 u_0\|_{L^2} \exp\left( \int_0^t (\|\Delta U(s)\|_{L^3} + \|\nabla U(s)\|_{L^\infty}) \, ds \right) $（引理 2.4）。

2. 混合速率下界（第 3 节核心定理，文中虽未编号但明确推导）：混合程度由齐次 Sobolev 范数 $ \|u(t)\|_{H^{-\alpha}} $ 的衰减刻画，其中 $ \alpha \in [1/2, 1] $。对 $ U $ 施加 $ \|U(t)\|_{W^{1,q}} \leq C $ 的约束，作者证明了以下最优下界：

   • 当 $ q > d $ 时，存在常数 $ C_{q,\alpha} > 0 $，使得
     $$ \inf_{\substack{\|U\|_{W^{1,q}} \leq 1 \\ \mathrm{div}\, U = 0}} \|u(t)\|_{H^{-\alpha}} \geq C_{q,\alpha} \, t^{-\frac{\alpha q}{q - d}} \quad \text{（代数衰减）}. $$
   • 当 $ q = d $ 时，出现临界对数修正：
     $$ \inf_{\substack{\|U\|_{W^{1,d}} \leq 1 \\ \mathrm{div}\, U = 0}} \|u(t)\|_{H^{-\alpha}} \geq C_{d,\alpha} \, (\log t)^{-\alpha} \quad \text{（对数衰减）}. $$
   • 当 $ 1 < q < d $ 时，衰减率由 Sobolev 共轭指数 $ q^* = dq/(d-q) $ 决定：
     $$ \inf_{\substack{\|U\|_{W^{1,q}} \leq 1 \\ \mathrm{div}\, U = 0}} \|u(t)\|_{H^{-\alpha}} \geq C_{q,\alpha} \, t^{-\frac{\alpha q^*}{2}} \quad \text{（更慢代数衰减）}. $$
     这些下界表明：$ W^{1,q} $ 正则性越弱（即 $ q $ 越小），$ U $ 驱动 $ u $ 生成小尺度的能力越受限；而 $ \alpha $ 越大，对高频模态的探测越敏感，衰减下界越强。

3. 数值验证的指数速率（第 4 节）：在 $ W^{1,2} $ 约束下，作者设计瞬时最优控制策略：在每个时刻 $ t $，选取 $ U(t) $ 以最大化 $ \frac{d}{dt} \|u(t)\|_{H^{-\alpha}}^2 $ 的负值。数值模拟显示，该策略产生的 $ \|u(t)\|_{H^{-\alpha}} $ 至少以 $ e^{-\gamma t} $ 指数衰减（$ \gamma > 0 $ 依赖于 $ \alpha $ 和初始数据），为理论下界提供了强有力证据，暗示在非临界 $ q $ 下，真实最优速率可能远超代数阶。

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

本文证明的核心障碍在于：压力项 $ \nabla p $ 破坏了 Lagrangian 结构，使标准特征线方法、流映射分析及标量混合中成熟的双曲估计（如 $ \|\theta(t)\|_{H^{-1}} \gtrsim e^{-Ct} \|\theta_0\|_{H^{-1}} $）完全失效。作者的真正创新点在于构建了一套基于 Fourier 截断、频谱能量重分配与 Sobolev 嵌入的混合分析框架，具体包含以下关键技术：

1. Leray 投影的精细频谱分解：作者没有回避压力，而是将其显式纳入 Fourier 空间分析。利用 $ \mathbb{P}_k v = v - (v \cdot k) k / |k|^2 $，将 $ u $ 的演化分解为平行与垂直于波矢 $ k $ 的分量。关键观察是：压力仅影响 $ u $ 的纵向分量（$ \propto k $），而横向分量（$ \perp k $）直接受输运项驱动。这一分解使作者能分离出压力对 $ \|u\|_{H^{-\alpha}} = \left( \sum_{k \neq 0} |k|^{-2\alpha} |\hat{u}_k|^2 \right)^{1/2} $ 的贡献，并证明其本质是将能量从低波数 $ |k| \sim 1 $ 向高波数 $ |k| \gg 1 $ 转移，从而驱动范数衰减。

2. 加权能量估计与 Sobolev 嵌入耦合：为获得 $ H^{-\alpha} $ 下界，作者构造加权能量泛函 $ \mathcal{E}_\alpha(t) = \sum_k |k|^{-2\alpha} |\hat{u}_k(t)|^2 $。对其求导并利用 $ \mathbb{P}_k $ 的正交性，得到

   $$ \frac{d}{dt} \mathcal{E}_\alpha(t) = -2 \Re \sum_{\ell,m} |k|^{-2\alpha} \, \hat{u}_k \cdot \mathbb{P}_k \big( \widehat{U}_\ell \cdot i m \, \hat{u}_m \big), \quad k = \ell + m. $$

   关键一步是应用 Hausdorff–Young 不等式与 Sobolev 嵌入 $ W^{1,q} \hookrightarrow L^{q^*} $（$ q^* = dq/(d-q) $），将三重和控制为

   $$ \left| \frac{d}{dt} \mathcal{E}_\alpha(t) \right| \leq C \|U(t)\|_{W^{1,q}} \, \mathcal{E}_\alpha(t)^{1/2} \, \left( \sum_k |k|^{-2\alpha + 2s} |\hat{u}_k|^2 \right)^{1/2}, $$

   其中 $ s $ 由嵌入指数决定。通过选择 $ s $ 使 $ -2\alpha + 2s = 0 $（即 $ s = \alpha $），并利用 $ \|U\|_{W^{1,q}} \leq C $，最终导出微分不等式 $ \frac{d}{dt} \mathcal{E}_\alpha(t) \geq -C \mathcal{E}_\alpha(t)^{\beta} $，其解给出前述代数/对数下界。

3. 谱截断法的收敛性与稳定性分析：为处理弱解，作者采用谱方法（2.3）构造近似解 $ u_N $，并证明其在 $ L^\infty_t L^2_x $ 中弱*收敛到真解 $ u $。关键工具是引理 2.6 给出的收敛率估计：

   $$ \|u - u_N\|_{L^2} \leq C \|U - P_N U\|_{L^1_t L^\infty_x} + C N^{-1} \|u_0\|_{H^1} \exp\left( \|U\|_{L^1_t W^{1,\infty}_x} \right), $$

   其中 $ P_N $ 为 Fourier 截断算子。该估计确保了数值模拟的数学可靠性，并为后续混合速率分析提供了坚实基础。

4. 与之前工作的比较

本文与既有工作存在本质差异，主要体现在问题设定、技术路径与结论深度三个维度：

1. 与标量混合工作的根本区别：Crippa–De Lellis [7] 及其后继者（Seis [27], Iyer–Kiselev–Xu [17]）针对标量方程 $ \partial_t \theta + (U \cdot \nabla) \theta = 0 $，利用 Lagrangian 观点证明了 $ \|\theta(t)\|_{H^{-1}} \geq C e^{-Ct} $。但该方法完全依赖于 $ \theta $ 沿 $ U $ 的迹线守恒（即 $ \theta(t, X(t;x)) = \theta_0(x) $），而本文的向量系统 (1.1) 因压力项彻底丧失此性质。作者明确指出：“the simple Lagrangian perspective available in the scalar case is no longer applicable”，并首次系统论证了压力如何迫使分析转向 Fourier 频谱重分配框架。因此，本文不是对标量结果的简单推广，而是开辟了向量混合的新范式。

2. 与欧拉/NS 方程研究的定位差异：Bardos–Titi [3]、Buckmaster–Vicol [4] 等关于欧拉方程弱解非唯一性的工作聚焦于非线性相互作用； Constantin–Fefferman–Majda [6] 关于涡量拉伸的研究强调 $ (u \cdot \nabla) u $ 项。而本文刻意剥离非线性，将 (1.1) 定位为“压力隔离测试平台”。如引言所述：“we view model (1.1) as a more flexible version… one may first analyze the more flexible model (1.1) and then, in future work, try to close the gap between $ U $ and $ u $.” 这种策略使本文能单独量化压力对小尺度形成的贡献，这是欧拉方程研究中无法做到的。

3. 与向量场输运的早期工作相比的推进：Lin–Thiffeault–Doering [22] 的数值工作虽涉及向量场，但未建立理论下界；Elgindi–Zlatos [14] 的构造性工作针对标量。本文首次为向量混合建立了严格的、依赖于 $ q $ 和 $ \alpha $ 的代数/对数下界，其形式 $ t^{-\alpha q/(q-d)} $ 揭示了 Sobolev 指数 $ q $ 的临界作用——这在标量工作中仅表现为 $ q > 1 $ 的门槛，而本文则刻画了 $ q $ 连续变化时的精细衰减速率谱。代价在于：本文下界是“inf over $ U $”意义上的最坏情形保证，而非构造性达优；且未解决 $ q = 1 $ 的端点情形（此时 $ W^{1,1} $ 不嵌入连续函数空间，技术上更棘手）。作者坦承：“We do not know whether these weak solutions are unique or satisfy an energy equality instead of the inequality above”，这表明在最低正则性下，解的适定性仍是开放问题。

13. Viterbo 平面中 ABP 不等式辛几何证明的一个经典分析对应版本

A Classical Analysis Counterpart of Viterbo’s Symplectic Geometry Proof of ABP in the Plane

作者: Daniel Maienshein, Juan J. Manfredi

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.12712v1

摘要翻译: 我们首先基于 Viterbo 的辛几何证明方法，给出二维紧支集 $C^2$ 函数上 Alexandroff-Bakelman-Pucci 不等式（ABP）一个版本的经典分析证明；该方法避开了凸性或接触集的使用。随后，我们说明如何修改该证明以去除紧支集假设，并恢复 ABP 的通常陈述形式——其中包含边界项。我们还讨论了将纯经典分析证明推广至三维及更高维的可能性（及其困难）。

Comment: 无

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

1. 研究对象和背景

本文研究二维空间中Alexandroff–Bakelman–Pucci（ABP）不等式的经典分析重构，核心对象是紧支集 $ C^2 $ 函数 $ f \in C^2(\mathbb{R}^2) $ 在紧集 $ K \subset \mathbb{R}^2 $ 上的振幅 $ \operatorname{osc}_K f := \sup_K f - \inf_K f $，其控制目标为Hessian行列式绝对值的 $ L^1 $ 范数：

这一不等式是椭圆偏微分方程理论中弱极值原理的基石，尤其对非散度型、不可逆系数或退化椭圆方程的解的 $ L^\infty $ 估计至关重要。标准ABP估计（Theorem 1.1）针对有界区域 $ \Omega \subset \mathbb{R}^n $ 中满足 $ Lu \geq f $ 的强上解 $ u $，给出

其证明依赖于上接触集 $ \Gamma_+ := \{ x \in \Omega : u(x) = \Gamma_u(x) \} $（其中 $ \Gamma_u $ 是 $ u $ 的凹包络），并利用面积公式建立 $ \mathcal{L}^n(Du(\Gamma_+)) \leq \int_{\Gamma_+} |\det D^2 u| $，再通过几何比较 $ (\sup_\Omega u - \sup_{\partial \Omega} u)^n \lesssim \mathcal{L}^n(Du(\Gamma_+)) $ 完成。该框架天然要求凸性结构（接触集定义依赖于凹包络）与区域几何（如直径控制），且难以脱离测度论工具。2000年，Viterbo [Vit00] 利用辛几何方法给出了一个全新的ABP版本（Theorem 1.3）：对紧支集 $ f \in C^2(\mathbb{R}^n) $，存在仅依赖维数 $ n $ 的常数 $ \rho_n $，使得

其证明基于余切丛 $ T^*K $ 中拉格朗日子流形的哈密顿同痕不变性，核心是构造生成函数 $ S_z(x) = f(x) - z x_2 $，并将 $ \det D^2 f $ 解释为拉格朗日子流形的雅可比行列式。该方法完全规避了凸包络、接触集与区域边界等经典分析中的几何障碍，但代价是引入高阶辛拓扑语言，使其难以被PDE分析社群直接吸收与推广。本文旨在填补这一空白：在 $ n = 2 $ 情形下，将Viterbo的辛几何思想纯化为经典分析语言，即不使用任何辛流形、拉格朗日子流形、哈密顿同痕或生成函数的抽象定义，而仅依赖微积分基本定理、共面积公式（co-area formula）、一维总变差（total variation）与平面拓扑（如环绕数、二色染色）。关键挑战在于：如何在不预设凸性、不引入接触集、不依赖区域边界形状（如不假设 $ K $ 为凸集或有光滑边界）的前提下，将二维振幅控制还原为一系列可处理的一维估计？本文将此问题转化为对水平集 $ \Sigma_z = \{ x \in \mathbb{R}^2 : \partial_{x_2} f(x) = z \} $ 的精细分析，并通过构造特定的“$ (\Sigma, c) $-容许路径”来连接极值点与边界，从而绕过传统框架中对全局凸性的隐含依赖。这一重构不仅为ABP提供了一个新的、更初等的证明视角，更揭示了振幅控制本质上是一种拓扑-分析耦合现象：水平集的连通分支结构（由环绕数刻画）与函数沿垂直方向的变差共同决定了整体振幅。

2. 主要定理或主要结果

本文在 $ n = 2 $ 维建立了两个严格证明的ABP型不等式，其结论均具有明确的常数依赖、尺度关系与函数空间假设。

Theorem 2.1（紧支集情形）: 设 $ K \subset \mathbb{R}^2 $ 为具有非空内部的紧集，则对任意 $ f \in C^2(\mathbb{R}^2) $ 满足 $ \operatorname{supp} f \subset K $，有

常数 $ C_K = 16 \, \mathrm{diam}(K)^2 $ 显式给出，仅依赖于 $ K $ 的直径平方，而非体积；该结果适用于所有紧支集 $ C^2 $ 函数，无需对 $ K $ 的正则性（如凸性、光滑性）作任何假设。

Theorem 2.6（一般情形）: 设 $ K \subset \mathbb{R}^2 $ 为具有非空内部的紧集，则对任意 $ f \in C^2(\mathbb{R}^2) $，有

其中常数 $ C_K $ 同样仅依赖于 $ \mathrm{diam}(K)^2 $。此处 $ \operatorname{osc}_{\partial K} f := \sup_{\partial K} f - \inf_{\partial K} f $ 是边界上的振幅，该形式恢复了经典ABP中不可或缺的边界项，且证明过程通过修改紧支集情形中的染色方案与容许路径构造（引理2.4、2.7、2.8）实现，未引入额外的正则性假设（如 $ \partial K $ 光滑）。

此外，论文在引理2.3中确立了一维垂直ABP的核心引理：对光滑紧连通曲线 $ M \subset \mathbb{R}^2 $ 及 $ u \in C^2(\mathbb{R}^2) $ 满足 $ M \subset u_{x_2}^{-1}(0) $，其垂直振幅 $ \operatorname{vosc}_M(u) := \sup\{ u(x) - u(y) : x, y \in M, \, p_1(x) = p_1(y) \} $ 满足

其中 $ \gamma $ 是 $ M $ 的单位速参数化，$ \mathcal{H}^1 $ 为一维Hausdorff测度。该引理的常数仅为投影长度 $ \mathcal{H}^1(p_1(M)) $，是二维证明中所有振幅控制的原子单元。

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

本文证明的核心障碍在于：如何在完全放弃凸包络与接触集的前提下，将二维函数的全局振幅 $ \operatorname{osc}_K f $ 控制为Hessian行列式的积分？作者的真正创新点并非发明新工具，而是对已有经典工具进行结构性重组与语义重赋，形成一个自洽的、纯分析的证明框架。该框架包含三个相互支撑的关键技术层。

第一层：共面积公式驱动的维度约化。证明始于将 $ \int_K |\det D^2 f| $ 通过共面积公式（Fed69, Theorem 3.2.22）分解为对 $ x_2 $-导数水平集 $ \Sigma_z = f_{x_2}^{-1}(z) $ 的积分：

其中 $ \mathrm{TV}_z(f_{x_1}) $ 是 $ f_{x_1} $ 在 $ \Sigma_z $ 上的总变差。由于 $ n = 2 $，有 $ |\det D^2 f| = |\nabla f_{x_1} \cdot X_{f_{x_2}}| $，且 $ X_{f_{x_2}} = (-f_{x_2 x_2}, f_{x_1 x_2}) $ 是 $ \Sigma_z $ 的单位切向量场，故 $ \mathrm{TV}_z(f_{x_1}) = \int_{\Sigma_z} |\nabla_{z} f_{x_1}| \, d\mathcal{H}^1 $。这一步将二维问题彻底降维为一族一维问题，但代价是需精细控制水平集 $ \Sigma_z $ 的拓扑结构。

第二层：平面拓扑引导的振幅合成。面对紧支集假设，$ \Sigma_z $ 对几乎每个 $ z $ 是 $ K $ 内部的有限个互不相交的光滑简单闭曲线 $ \Omega_{z,i} $。作者引入环绕数 $ X_\Sigma(x,y) $ 和二元染色 $ c_{0,x_0}, c_{1,x_0} $（定义在 $ K \setminus \Sigma_z $ 上），将平面划分为由 $ \Sigma_z $ 分隔的、具有确定符号的区域。在此基础上，定义**$ (\Sigma_z, c) $-容许路径**：包括在 $ c = +1 $ 区域内向下、在 $ c = -1 $ 区域内向上的垂直线段，以及在闭曲线 $ \Omega_{z,i} $ 上从某点 $ x $ 到其“最高/最低同 $ x_1 $-坐标点” $ G(x) $ 的弧段。引理2.4保证：对任意内点 $ x^* \in K \setminus \Sigma_z $，存在一条从 $ \partial K $ 到 $ x^* $ 的分段光滑容许路径 $ \Gamma $，且该路径至多访问每个 $ \Omega_{z,i} $ 一次。这一拓扑引理是整个框架的枢纽，它替代了传统证明中“接触集必与梯度像相交”的测度论论证，将振幅控制转化为沿特定路径的增量累加。

第三层：生成函数与垂直ABP的耦合估计。为在每条闭曲线 $ \Omega_{z,i} $ 上应用一维估计，作者构造辅助函数 $ S_z(x) = f(x) - z x_2 $，其关键性质是 $ (S_z)_{x_2} = 0 $ 在 $ \Sigma_z $ 上恒成立，从而满足引理2.3的条件。由此得到 $ \operatorname{vosc}_{\Omega_{z,i}} S_z \leq \mathcal{H}^1(p_1(\Omega_{z,i})) \int_{\Omega_{z,i}} |\frac{d}{ds} f_{x_1}| $。沿容许路径 $ \Gamma $ 积分时，垂直线段上的增量由 $ z $ 控制（引理2.5：$ f(y) - f(x) \leq z \, p_2(y-x) $），而闭曲线弧段上的增量则由 $ \operatorname{vosc}_{\Omega_{z,i}} S_z $ 控制。最终，通过Markov不等式选择最优 $ z_0 \notin E_\lambda $ 并优化参数 $ \lambda $，完成 $ \operatorname{osc}_K f $ 的平方控制。整个框架的创新本质在于：用环绕数染色代替凸包络几何，用容许路径积分代替接触集覆盖，用生成函数 $ S_z $ 作为纯粹的分析工具（而非辛几何对象）来桥接二维与一维估计。

4. 与之前工作的比较

本文与既有文献的比较需从证明范式、假设强度、适用范围及技术代价四个维度展开。

与经典ABP证明（Cab95, GT98）的比较：经典证明（Theorem 1.1, Lemma 1.2）严重依赖区域 $ \Omega $ 的有界性与函数 $ u $ 在 $ \partial \Omega $ 上的取值，其核心引理1.2的右边为 $ \int_{\Gamma_+} |\det D^2 u| $，积分域 $ \Gamma_+ $ 是由凹包络定义的、依赖于 $ u $ 全局行为的集合。本文Theorem 2.1则完全摆脱了边界与凸性的束缚——它仅需 $ f $ 紧支于 $ K $，对 $ K $ 的形状（凸、星形、甚至非连通）无任何要求，且积分域是整个 $ K $。代价是：经典证明可直接推广至任意维数 $ n $ 及一般椭圆算子 $ Lu $，而本文方法目前仅严格成立於 $ n = 2 $；且经典证明的常数 $ C(n, \lambda, B \, \mathrm{diam}(\Omega)) $ 隐含对系数的控制，本文常数仅显式依赖直径。

与Viterbo辛几何证明（Vit00）的比较：本文是Viterbo工作的直接“翻译”，但绝非字面转译。Viterbo的Theorem 1.3给出常数 $ \rho_n \mathcal{L}^n(K) $，其证明依赖于拉格朗日子流形的哈密顿同痕不变性与生成函数的辛容量估计。本文Theorem 2.1的常数为 $ 16 \, \mathrm{diam}(K)^2 $，明确显示其与经典ABP的直径依赖一致，证明中所有概念（环绕数、染色、容许路径）均为初等平面拓扑，不涉及余切丛、辛形式或哈密顿流。作者明确指出：“our proof is not a one-to-one translation”，因为Viterbo的体积依赖在本文中被直径依赖取代，这反映了两种框架对“尺度”的不同敏感度。此外，Viterbo未探讨非紧支集情形，而本文Theorem 2.6首次在纯分析框架下实现了边界项的自然引入，这是对Viterbo原始构想的实质性推进。

与后续分析工作（如Wang09, Imb16）的比较：近年来有工作尝试用调和分析或概率方法重证ABP，但多聚焦于特殊算子（如$ p $-拉普拉斯）或加权空间。本文的独特性在于其问题导向的纯粹性：它不追求最一般的算子或最优常数，而是执着于复现Viterbo思想的“最小分析模型”。这种克制反而凸显了二维特有结构（如环绕数、简单闭曲线的定向）的不可替代性。论文末尾坦承：“the possibility (and difficulties) of extending to dimension 3 and above” 仍开放，因为三维中水平集 $ \Sigma_z $ 是曲面而非曲线，环绕数概念失效，且无法构造类似引理2.4的拓扑路径。因此，本文的贡献不是普适性推广，而是精确地刻画了ABP不等式在二维所蕴含的、可被经典分析完全捕获的深层拓扑-分析结构。

14. Hilbert 空间中具退化扩散的均值场博弈：粘性解方法

Mean Field Games in Hilbert Spaces with Degenerate Diffusion: A Viscosity Solution Approach

作者: Andrzej Święch, Lukas Wessels

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

发布日期: 2026-05-12

arXiv ID: 2605.12690v1

摘要翻译: 我们在 Hilbert 空间 $H$ 中研究一类退化的二阶均值场博弈（MFG）系统，该系统耦合了一个描述 $H$ 上概率测度演化的 Fokker–Planck 方程与一个刻画值函数的 Hamilton–Jacobi–Bellman（HJB）方程。我们的主要结果建立了该耦合系统的解的存在性与唯一性。HJB 方程的解以粘性解意义进行解释。为证明存在性，我们将基于 Tikhonov 定理的经典不动点方法推广至本设定。该方法中的一个核心难点在于证明相应线性退化 Fokker–Planck 方程的唯一性。为应对这一问题，我们引入一类适当的伴随方程，并运用粘性解技术构造出足够正则的解。随后，通过适配 Lasry–Lions 单调性方法，得到完整 MFG 系统的唯一性。

Comment: 无

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

1. 研究对象和背景

本文研究定义在实可分希尔伯特空间 $H$ 上的一类退化二阶均值场博弈（Mean Field Game, MFG）系统，其核心对象是如下耦合的偏微分方程组：

其中 $v : [0,T] \times H \to \mathbb{R}$ 是代表性智能体的价值函数，$m : [0,T] \to \mathcal{P}_1(H)$ 是群体概率测度流，$\mathcal{P}_1(H)$ 表示关于 $H$-范数具有一阶矩的Wasserstein空间。算子 $L$ 形式定义为

其中 $A : D(A) \subset H \to H$ 是一个闭的、稠定的、极大耗散线性算子（如椭圆型或双曲型无穷维微分算子），$\sigma : H \to \mathcal{L}_2(\Xi, H)$ 是取值于Hilbert–Schmidt算子空间的退化扩散系数（即 $\sigma(x)$ 不可逆，甚至可能秩有限），且满足 $\|\sigma(x) - \sigma(y)\|_{\mathcal{L}_2} \le C \|x-y\|_{-1}$；$L^*$ 是其形式伴随，$\operatorname{div}$ 定义为负梯度的伴随。关键未知量包括：(i) 在退化、非一致椭圆、无穷维框架下，HJB方程的粘性解 $v$ 的存在性与正则性；(ii) 对应的线性Fokker–Planck方程（FP）在给定漂移场 $w(t,x) = H_p(x,Dv,m(t))$ 下弱解 $m$ 的唯一性；(iii) 整个耦合系统的全局适定性。

该问题的重要性源于其作为大规模随机微分博弈 $N \to \infty$ 极限的严格数学基础：当每个智能体的状态由含无界算子 $A$ 和退化噪声 $\sigma$ 的无穷维SDE驱动时（例如随机延迟微分方程、随机偏微分方程），经典有限维MFG理论失效。已有理论在此方向存在显著空白：一方面，[26] 在mild解框架下处理了类似系统，但要求 $\sigma \equiv I$ 且 $A$ 满足强正则性条件（如 $(-A)^{-1+\delta}$ 迹类），以保障Ornstein–Uhlenbeck半群的光滑性与紧性，这排除了大量物理模型（如退化扩散、非自伴传输算子）；另一方面，无穷维粘性解理论虽在[24]中建立，但尚未用于MFG耦合系统，尤其缺乏对退化 $\sigma$ 与无界 $A$ 同时存在时线性FP唯一性的处理工具。本文试图填补的核心空白是：在不假设 $\sigma$ 非退化、不依赖Ornstein–Uhlenbeck半群光滑性、不预设 $A$ 自伴或谱条件的前提下，构建一套适用于退化无穷维MFG的粘性解分析框架，并确立其解的存在唯一性。

2. 主要定理或主要结果

本文在Assumptions 2.1(A1)–(A6)与3.1下，证明了以下主要结果：

定理 4.2（存在性）：设 $H$ 满足(A1)–(A6)，$G$ 满足(A5)，$m_0 \in \mathcal{P}_1(H)$。则MFG系统(HJB)–(FP)存在一对解 $(v,m)$，其中 $v$ 是(HJB)的$B$-连续粘性解，且 $v(t,\cdot) \in C^{1,1}_b(H_{-1})$ 对所有 $t \in [0,T]$ 成立；$m \in \mathcal{C}([0,T]; \mathcal{P}_1(H_{-1}))$ 是(FP)的弱解（Definition 3.2），满足 $\sup_{t \in [0,T]} \int_H \|x\|\, m(t,dx) < \infty$。证明基于Tikhonov不动点定理：定义映射 $\Phi : S \to S$，$\Phi(m) = m'$，其中 $m'$ 是线性FP方程 $\partial_t m' - L^* m' - \operatorname{div}(H_p(x,Dv,m)m') = 0$ 的弱解，$v$ 是对应HJB方程的粘性解。紧性由 $S \subset \mathcal{C}([0,T]; \mathcal{P}_1(H_{-1}))$ 中Wasserstein度量 $d_{1,-1}$ 的完备性及 $H_{-1}$ 的紧嵌入保证。

定理 3.8（线性FP唯一性）：设 $w$ 满足Assumption 3.1及额外正则性(3.7)–(3.8)，则对任意 $f \in C_b([0,T]\times H)$ 满足 $|f(t,x)-f(t,y)| \le C_2 \|x-y\|_{-1}$，线性FP方程 $\partial_t m - L^* m - \operatorname{div}(w m) = 0$，$m(0)=m_0$ 在 $\mathcal{SF}_P$ 中至多有一个弱解。该结果独立成立，是本文关键新贡献。

定理 5.1（唯一性）：在相同假设下，MFG系统(HJB)–(FP)的解 $(v,m)$ 在如下意义下唯一：若 $(v_1,m_1)$ 与 $(v_2,m_2)$ 均为解，则 $v_1 = v_2$ 在 $[0,T]\times H$ 上一致成立，且 $m_1(t) = m_2(t)$ 在 $\mathcal{P}_1(H_{-1})$ 中对所有 $t$ 成立。证明采用Lasry–Lions单调性方法的适应版本：构造正则化序列 $v_{1,n,N}, v_{2,n,N}$（通过Yosida逼近 $A_n$、投影 $P_N$ 与磨光 $\eta$），验证 $\int_0^T \int_H (v_{1,n,N} - v_{2,n,N})\, d(m_{1,n,N} - m_{2,n,N})$ 趋于零，最终利用 $v_{i,n,N} \to v_i$ 及 $m_{i,n,N} \to m_i$ 的收敛性得出结论。误差控制依赖于估计 $\|Dv_i(t,\cdot) - Dv_{i,n,N}(t,\cdot)\|_{\infty} \le C/n + C/N + C\eta$ 与 $\|v_i - v_{i,n,N}\|_{\infty} \le C T \|f_i - f_{i,n,N}\|_{\infty}$。

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

本文的核心技术障碍有三重：(i) 无穷维退化FP的唯一性缺失：经典方法依赖伴随方程的光滑解，但在 $\sigma$ 退化且 $A$ 无界时，$\partial_t u + \langle Ax, Du\rangle + \frac{1}{2} \operatorname{Tr}[\sigma\sigma^* D^2 u] - \langle w, Du\rangle = f$ 通常无经典解；(ii) 粘性解正则性不足：HJB粘性解仅保证 $v \in C^{1,1}_b(H_{-1})$，而标准单调性论证需 $Dv$ 的Lipschitz连续性及 $v$ 的二阶导数存在；(iii) 无穷维空间紧性丧失：$\mathcal{P}_1(H)$ 在Wasserstein拓扑下非紧，无法直接应用Tikhonov定理。

作者的真正创新点在于构建了一套三层正则化-紧化-极限框架：

第一层：伴随方程的粘性解正则化。为克服(i)，作者不寻求经典解，而是系统性地使用粘性解理论处理伴随方程(3.2)。关键工具是三个嵌套逼近：(a) Yosida逼近 $A_n = nA(nI-A)^{-1}$ 将无界 $A$ 替换为有界算子，确保 $e^{tA_n}$ 一致有界；(b) 有限维投影 $P_N$ 将 $H$ 截断至 $H_N = \operatorname{span}\{e_1,\dots,e_N\}$，使问题降维；(c) 标准磨光 $\rho_\eta$ 在 $[0,T]\times \mathbb{R}^N$ 上作用于 $\tilde{u}_{n,N}(t,P_N x)$，生成 $C^\infty$ 函数 $u^\eta_{n,N}$。引理3.7证明 $u^\eta_{n,N} \in D_T$（测试函数类），且满足 $\partial_t u^\eta_{n,N} + L_0 u^\eta_{n,N} - \langle w, D u^\eta_{n,N}\rangle = f^\eta_{n,N}$，其中 $f^\eta_{n,N} \to f$ 点态收敛。此构造将密度论证从经典函数空间转移到粘性解生成的正则化序列上。

第二层：$H_{-1}$ 空间的紧化枢纽。为克服(iii)，作者引入辅助空间 $H_{-1}$，其范数定义为 $\|x\|_{-1}^2 = \langle B x, x\rangle$，其中 $B$ 是紧的、正定的自伴算子，满足弱$B$-条件 $-A^*B + c_0 B \ge 0$（Assumption 2.1(A1)）。由于 $B$ 紧，$H \hookrightarrow H_{-1}$ 紧嵌入，从而 $\mathcal{P}_1(H) \subset \mathcal{P}_1(H_{-1})$ 在 $d_{1,-1}$ 度量下具备相对紧性（Arzelà–Ascoli）。所有解被强制取值于 $\mathcal{C}([0,T]; \mathcal{P}_1(H_{-1}))$，Tikhonov定理得以应用。

第三层：单调性论证的粘性解适配。为克服(ii)，作者放弃对原始 $v$ 的直接计算，而是在正则化层级进行：对两组解 $(v_1,m_1)$、$(v_2,m_2)$，分别构造 $v_{1,n,N}, v_{2,n,N}$，并定义差值 $\delta v_{n,N} = v_{1,n,N} - v_{2,n,N}$。利用 $v_{i,n,N}$ 满足的近似HJB方程及 $m_i$ 的弱解定义，推导出

再结合 $H$ 的单调性假设（如 $\langle H_p(x,p,\mu) - H_p(x,p',\mu'), p-p' \rangle \ge \lambda \|p-p'\|^2$）与 $v_{i,n,N}$ 的 $C^{1,1}_b(H_{-1})$ 估计（引理3.5给出 $\|v_{i,n,N}(t,\cdot)\|_{C^{1,1}_b(H_{-1})} \le L$），最终令 $n,N,\eta \to \infty$ 取极限完成唯一性证明。

4. 与之前工作的比较

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

1. 与[26]的对比（最直接参照）：[26]在mild解框架下处理相同形式的MFG系统，但其假设与本文本质互斥。[26]要求 $\sigma \equiv I$ 且 $A$ 是负自伴算子，满足 $(-A)^{-1+\delta}$ 迹类（$\delta>0$），以调用[6]中FP唯一性结果并保障Ornstein–Uhlenbeck半群 $e^{tL}$ 的强连续性与紧性。本文则允许 $\sigma$ 退化（如 $\operatorname{rank}(\sigma(x)) < \infty$）且 $A$ 仅为极大耗散（如运输算子 $A = -\partial_\theta$），适用范围更广。代价是：本文需 $H$ 关于 $H_{-1}$-范数连续（Assumption 2.1(A2)），且依赖 $B$-紧性构造 $H_{-1}$；而[26]对 $H$ 的正则性要求极低（仅Lipschitz）。二者结果互补：[26]覆盖强抛物情形，本文覆盖退化/双曲主导情形。

2. 与有限维退化MFG工作的对比（如[10,28,29]）：[10]利用MFG系统作为凸优化一阶条件的结构，要求 $H$ 关于 $p$ 严格凸；[28,29]采用高阶椭圆正则化 $\varepsilon \Delta$，需 $\varepsilon > 0$ 并令其趋于零，本质上仍依赖非退化扩散。本文完全摒弃椭圆正则化，直接在退化极限下工作，技术路径根本不同。代价是：本文对 $H$ 的结构要求更强（需满足Assumption 2.1(A2)–(A3)的Lipschitz与可微性），而[10]等可处理更一般的Hamiltonian。

3. 与无穷维线性FP文献的对比（如[4,6,57]）：这些工作证明退化FP唯一性，但均要求 $\sigma(x) \equiv \sigma$ 常数且满足 $e^{tA}(H) \subset \left(\int_0^t e^{sA} \sigma \sigma^* e^{sA^*} ds\right)^{1/2}(H)$，此条件在 $A$ 无界且 $\sigma$ 退化时几乎无法满足（如 $A = \Delta$，$\sigma$ 为乘子算子）。本文通过引入 $H_{-1}$ 框架与粘性解正则化，完全绕过该几何条件，为无穷维退化FP唯一性提供了新范式。

4. 与粘性解理论文献的对比（如[24,23]）：[24]建立了无穷维粘性解基本理论，[23]证明了HJB粘性解的 $C^{1,1}_b(H_{-1})$ 正则性，但二者均未涉及MFG耦合。本文首次将 $B$-连续粘性解、Yosida逼近、$H_{-1}$ 紧化三者整合为MFG分析工具链。局限在于：当前结果依赖Assumption 2.1(A6)（即预设HJB解存在且具 $C^{1,1}_b(H_{-1})$ 正则性），尚未将此正则性作为定理结论推导出来；且对Hamiltonian的局部Lipschitz性（如 $H$ 在 $p$ 上仅局部Lipschitz）尚未处理，而[26]对此更具鲁棒性。

15. 关于半线性热方程整体解的存在性与非存在性

On the existence and nonexistence of global solutions of the semilinear heat equation

作者: Kaiqiang Zhang, Zhiyu Li

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11933v1

摘要翻译: 我们考虑半线性热方程

该问题的著名难点在于：由于尺度不变性导致势阱深度为零，势阱方法无法直接应用。我们采用前向相似变换将原方程转化为一个新的抛物型方程，从而可在加权Sobolev空间中应用势阱方法。由此得到一个新判据，用以判定热方程的解是否在有限时间内爆破或整体存在。本工作将Ikehata等人（\textit{Ann. Inst. H. Poincaré Anal. Non Linéaire}, \textbf{27} (2010) 877–900）在临界Sobolev指数情形下获得的部分结果推广至 $p_F<p<p_S$ 的情形，其中 $p_F=1+2/n$ 为Fujita指数，$p_S=(n+2)/(n-2)$（当 $n\ge3$ 时）为临界Sobolev指数。

Comment: 无

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

1. 研究对象和背景

本文研究全空间 $\mathbb{R}^n$（$n \ge 3$）上初值问题

的全局适定性与有限时间爆破行为。该方程是反应扩散型抛物方程的典范模型，其非线性项 $|u|^{p-1}u$ 具有临界尺度不变性：若 $u(t,x)$ 是解，则对任意 $\lambda > 0$，函数

也是解。这一对称性诱导出能量泛函

的尺度不变性 $E(u_\lambda) = E(u)$，并导致经典势井方法在标准 Sobolev 空间 $H^1(\mathbb{R}^n)$ 中失效——因可构造 Nehari 流形 $\mathcal{N} = \{ u \in H^1 \setminus \{0\} : I(u) = 0 \}$ 上一族函数满足 $E(u_j) \to 0$，故势井深度 $d := \inf\{ E(u) : u \in \mathcal{N} \}$ 为零。此障碍使得传统势井理论无法区分低能量非负初值（即 $0 < E(u_0) < d = 0$ 不可能）下的解行为，仅能复现 Levine 型结果：当 $E(u_0) < 0$ 时必爆破；而 $E(u_0) \ge 0$ 时的判据长期缺失。已知两个关键临界指数：Fujita 指数 $p_F = 1 + \frac{2}{n}$，它刻画了小初值全局解存在的阈值——当 $p \le p_F$ 时，任意非平凡非负初值均导致爆破；Sobolev 临界指数 $p_S = \frac{n+2}{n-2}$（$n \ge 3$），对应能量临界情形，此时方程具有共形不变性，且存在非平凡稳态解。Ikehata–Ishiwata–Suzuki（2010）在 $p = p_S$ 下通过修正势井法建立了低能量阈值理论：存在正深度势井 $W = \{ w \in H^1(K) : E(w) < d,\, I(w) > 0 \} \cup \{0\}$，其中 $K(y) = e^{|y|^2/4}$ 为高斯权函数，使得初值落入 $W$ 则全局存在，落入 $Z = \{ w \in H^1(K) : E(w) < d,\, I(w) < 0 \}$ 则有限时间爆破。但该方法依赖于 $p = p_S$ 时加权空间 $H^1(K)$ 中嵌入 $H^1(K) \hookrightarrow L^{p+1}(K)$ 的紧性与临界性，无法直接推广至亚临界区域 $p_F < p < p_S$。本文旨在填补这一空白：在能量亚临界但超 Fujita 的参数范围内，建立首个基于势井框架的、具有正深度的初值分类准则，从而将势井方法的有效性从临界情形系统性地延拓至更广阔的物理相关参数域。

2. 主要定理或主要结果

设 $n \ge 3$，$p_F = 1 + \frac{2}{n} < p < \frac{n+2}{n-2} = p_S$，初值 $u_0 \in H^1(K)$，其中加权空间 $H^1(K)$ 定义为

定义前向相似变换：$s = \log(1+t)$，$y = x / \sqrt{t+1}$，$w(y,s) = e^{s/(p-1)} u(e^{s/2} y, e^s - 1)$，则原方程化为加权抛物方程

相应加权能量泛函与 Nehari 泛函为

势井深度 $d := \inf\{ E(w) : w \in H^1(K) \setminus \{0\},\, I(w) = 0 \}$ 严格为正（见 Proposition 3.1）。定义势井 $W = \{ w \in H^1(K) : E(w) < d,\, I(w) > 0 \} \cup \{0\}$ 与外势井 $Z = \{ w \in H^1(K) : E(w) < d,\, I(w) < 0 \}$。则主定理断言：\textbf{Theorem 2.2.} （1）若 $u_0 \in W$，则原方程 (1.1) 存在唯一全局解 $u \in C([0,\infty); H^1(K))$，且满足 $\|u(t)\|_{H^1(K)} \le C$ 一致有界；（2）若 $u_0 \in Z$，则解在有限时间 $T < \infty$ 内爆破，即 $\lim_{t \to T^-} \|u(t)\|_{L^\infty(\mathbb{R}^n)} = +\infty$。进一步，由变换关系 $T = e^{T_s} - 1$，$T_s$ 为 $w$-方程的最大存在时间，故 $T < \infty$ 当且仅当 $T_s < \infty$。

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

核心障碍在于：原始方程 (1.1) 的尺度不变性导致标准 $H^1$ 势井深度为零，使经典势井方法完全失效。作者的根本创新在于将前向相似变换从单纯分析工具升格为构造正深度势井的结构性手段，并在加权空间 $H^1(K)$ 中重建完整的势井理论框架。具体而言：（1）加权空间的选择与关键不等式：采用 Kavian 引入的高斯权 $K(y) = e^{|y|^2/4}$，使得算子 $L = -\Delta - y\cdot\nabla/2$ 成为 $L^2(K)$ 上的自伴正定算子，且引理 2.1(iii) 给出加权 Poincaré 型不等式 $\frac{n}{2} \int |w|^2 K \, dy \le \int |\nabla w|^2 K \, dy$，这是保证势井深度 $d > 0$ 的基石。该不等式在 $p > p_F$ 时提供必要的控制，避免能量泛函退化。（2）势井深度的正性证明（Proposition 3.1）：对任意 $w \in H^1(K) \setminus \{0\}$ 满足 $I(w) = 0$，利用 $I(w)=0$ 将 $\int |\nabla w|^2 K$ 表示为 $\int |w|^2 K$ 与 $\int |w|^{p+1} K$ 的组合，再结合嵌入 $H^1(K) \hookrightarrow L^{p+1}(K)$（引理 2.1(ii)）及加权 Poincaré 不等式，导出 $\int |w|^{p+1} K \, dy \ge C_0 > 0$，进而 $E(w) \ge (p-1)C_0 / (2(p+1)) > 0$，故 $d > 0$。此论证精确揭示了 $p > p_F$ 与 $p < p_S$ 的双重必要性：前者保障加权 Poincaré 控制，后者保障嵌入紧性以闭合估计链。（3）不变集论证与全局存在性（Theorem 2.2(i)）：利用 $E(w(s))$ 单调递减性及 $I(w(s)) > 0$ 的保持性，由 $I(w) > 0$ 推出 $\int |\nabla w|^2 K \, dy$ 主导其他项，再结合引理 2.1(iii) 得到 $\|w(s)\|_{H^1(K)}$ 一致有界，从而排除有限时间爆破。（4）凹性论证的加权适配（Theorem 2.2(ii)）：定义辅助函数 $M(s) = \frac{1}{2} \int_0^s \int_{\mathbb{R}^n} |w|^2 K \, dy d\tau$，则 $M''(s) = -I(w(s))$。当 $I(w(s)) < -\varepsilon_1 < 0$ 时，$M''(s) > \varepsilon_1$，结合能量恒等式 $E(u_0) - E(w(s)) = \int_0^s \int |w_s|^2 K \, dy d\tau$，导出 $M(s) M''(s) \ge (1+\alpha) (M'(s))^2$ 对充分大的 $s$ 成立，从而 $Y(s) = M^{-\alpha}(s)$ 凹且正，矛盾。此步骤的关键在于将 Levine 的凹性法完全移植到加权框架，并利用 $p < p_S$ 保证所需嵌入成立。

4. 与之前工作的比较

本文与既有工作的比较需从四个维度展开：（1）适用参数范围：Ikehata–Ishiwata–Suzuki（2010）[8] 的结果严格限定于 $p = p_S$（能量临界），其证明严重依赖 $H^1(K) \hookrightarrow L^{p_S+1}(K)$ 的临界 Sobolev 嵌入及相应极小化问题的紧性；本文则首次将势井方法成功拓展至整个亚临界区间 $p_F < p < p_S$，覆盖了燃烧理论中更典型的超线性但非临界反应速率。（2）初值能量门槛：Levine（1973）[12] 仅处理 $E(u_0) < 0$ 的爆破，Lee–Ni（1992）[11] 给出基于初值衰减率的爆破判据，但均未构建统一的势井结构；本文通过加权势井 $W/Z$ 实现了对 $0 \le E(u_0) < d$ 这一长期悬而未决的‘低正能’区域的精细划分，$d > 0$ 的存在本身即是突破。（3）技术代价与局限：为获得正深度势井，本文必须工作在加权空间 $H^1(K)$ 而非标准 $H^1$，这导致初值需满足 $u_0 \in H^1(K)$，即要求 $\int_{\mathbb{R}^n} |u_0|^2 e^{|x|^2/4} \, dx < \infty$，比通常的 $L^2$ 或 $H^1$ 初值条件更强（例如 Gaussian 衰减初值满足，但紧支初值不满足）。此为获得正深度所付出的必要代价。此外，本文未处理 $p = p_F$ 或 $p \le p_F$ 的情形，因此时加权 Poincaré 不等式失效，$d$ 可能退化。（4）与 Kavian 方法的差异：Kavian（1987）[10] 首次使用前向相似变换研究大时间行为，但聚焦于加权能量单调性及渐近形态；本文则将该变换作为构造新动力系统的工具，彻底重构了势井的几何结构，使其深度从零变为正，并赋予其明确的动力学含义（不变集）。综上，本文并非对 [8] 的简单推广，而是通过加权空间与变换的协同设计，在参数、初值、方法三重维度上实现了对势井理论边界的实质性拓展。

16. 大气动力学中相变模型的整体存在性与强解唯一性

Global Existence and Uniqueness of Strong Solutions for a Phase Transition Model in Atmospheric Dynamics

作者: Giada Cianfarani Carnevale, Donatella Donatelli, Stefano Spirito

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11754v1

摘要翻译: 本文研究大气动力学中的一类相变模型，该模型受文献[6,14,15]的启发，后者分析了刻画速度、温度与比湿度演化的原始方程（primitive equations）。主要困难源于温度方程与比湿度方程中出现的一个多值不连续非线性项，该项描述了在过饱和条件下降水的形成过程。为克服这一困难，我们引入一种正则化形式，从而保证近似解的存在性与唯一性。随后借助经典的紧性论证，我们建立了原模型强解的存在性。此外，在一个条件性且具有物理意义的假设下，我们还建立了强解的唯一性。该方法使得我们能在全空间 $\mathbb{R}^2$ 上对热带气候模型给出严格的数学论证，同时避免在比湿度方程中引入粘性项。

Comment: 无

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

1. 研究对象和背景

本文研究一个描述热带大气中水汽相变过程的偏微分方程系统，其物理基础是原始方程（primitive equations）在小纵横比极限下的垂直模态截断——即热带气候模型（Tropical Climate Model, TCM）。该模型耦合了水平二维空间 $\mathbb{R}^2$ 上的 barotropic（正压）速度场 $u = u(x,y,t) \in \mathbb{R}^2$、baroclinic（斜压）速度场 $v = v(x,y,t) \in \mathbb{R}^2$、位温（或简记为温度）$T = T(x,y,t) \in \mathbb{R}$ 以及比湿（specific humidity）$q = q(x,y,t) \in [0,1]$。系统由如下五维非线性抛物-双曲混合型方程组给出（见原文(1.2)及(3.16)）：

其中 $p = p(x,y,t)$ 为等效压力；$H > 0$ 为对流层厚度；$g > 0$ 为重力加速度；$\theta_0 > 0$ 为参考位温；$N > 0$ 为 Brunt–Väisälä 频率（表征层结稳定性）；$R > 0$, $R_\nu > 0$, $c_p > 0$, $L > 0$ 分别为空气气体常数、水汽气体常数、定压比热容与蒸发潜热；$q_s \in (0,1)$ 为常数饱和比湿；$\overline{Q} > 0$ 为粗略湿度层结参数（gross moisture stratification），物理上反映湿度随高度指数衰减的特征尺度。

关键物理机制在于相变项：当空气上升冷却至露点以下时，水汽过饱和（$q > q_s$），触发凝结并释放潜热，同时产生降水 $P$。该过程在数学上体现为一个多值不连续非线性源项，其原型为 $-H_g/\pi \, \mathrm{div}\, v \, G^+(T) H(q - q_s)$，其中 $H(\cdot)$ 是 Heaviside 多值函数（即 $H(s) = 0$ 若 $s < 0$, $H(0) = [0,1]$, $H(s) = 1$ 若 $s > 0$），而 $G^+(T)$ 是依赖于温度的正则化因子（见原文(3.17)–(3.18)）。该非线性项仅在 $q > q_s$ 且 $\mathrm{div}\, v < 0$（对应上升运动）时激活，导致系统在 $q = q_s$ 处发生本质奇异性——既非 Lipschitz 连续，亦非单调（因符号可变），更非标准的 $L^2$ 可控源项。此即全文核心困难：它破坏了经典能量估计的闭合性，阻碍了标准 Galerkin 或紧性方法对强解的存在性与唯一性证明。

已有理论框架对此类问题处理有限：[6] 和 [14] 分别采用压力坐标系与微分包含（differential inclusions）证明弱解存在，但未获得强解唯一性；[15] 在无相变项时建立了热带气候模型强解的整体适定性；而 [10] 和 [12] 的模型虽含降水项，但或假设其为光滑函数，或仅处理浅水近似，无法刻画真实大气中水汽饱和阈值引发的突变行为。因此，本文旨在填补这一空白：在全空间 $\mathbb{R}^2$ 上，对包含物理真实相变机制（即 $H(q-q_s)$ 型多值项）、同时保留温度扩散（$-\Delta T$）但刻意避免在湿度方程中引入人工粘性（即无 $\Delta q$ 项） 的热带气候模型，建立强解的整体存在性与条件唯一性。该设定更具物理真实性（因水汽输运以平流主导，扩散效应远弱于温度），但数学难度陡增——缺失 $\Delta q$ 意味着无法通过抛物正则性控制 $q$ 的振荡，使 $H(q-q_s)$ 的奇异性无法被常规 Sobolev 嵌入所压制。

2. 主要定理或主要结果

本文在初始数据满足 $u_0 \in H^2(\mathbb{R}^2)^2$, $v_0 \in H^2(\mathbb{R}^2)^2$, $T_0 \in H^2(\mathbb{R}^2)$, $q_0 \in H^1(\mathbb{R}^2) \cap L^\infty(\mathbb{R}^2)$ 且 $0 \leq q_0 \leq 1$ a.e. 的前提下，证明了以下主要结果：

定理 4.1（整体强解存在性）：存在一个全局强解 $(u,v,T,q,p)$，满足

• $u, v \in L^\infty_{\mathrm{loc}}([0,\infty); H^2(\mathbb{R}^2)^2) \cap L^2_{\mathrm{loc}}([0,\infty); H^3(\mathbb{R}^2)^2)$,
• $T \in L^\infty_{\mathrm{loc}}([0,\infty); H^2(\mathbb{R}^2)) \cap L^2_{\mathrm{loc}}([0,\infty); H^3(\mathbb{R}^2))$,
• $q \in L^\infty_{\mathrm{loc}}([0,\infty); H^1(\mathbb{R}^2)) \cap L^2_{\mathrm{loc}}([0,\infty); H^2(\mathbb{R}^2))$,
• $p \in L^2_{\mathrm{loc}}([0,\infty); \dot{H}^1(\mathbb{R}^2))$, 且该解在分布意义下满足系统 (1.2)，其中多值项 $H(q - q_s)$ 被理解为一个具有取值于 $[0,1]$ 的选择函数 $\chi = \chi(x,y,t) \in L^\infty(\mathbb{R}^2 \times (0,\infty))$，满足 $\chi(x,y,t) \in H(q(x,y,t) - q_s)$ a.e.

定理 4.2（强解唯一性）：若两个强解 $(u^{(1)}, v^{(1)}, T^{(1)}, q^{(1)})$ 与 $(u^{(2)}, v^{(2)}, T^{(2)}, q^{(2)})$ 满足相同的初值，且存在常数 $C_T > 0$ 使得 $\|T^{(i)}(t)\|_{L^\infty(\mathbb{R}^2)} \leq C_T$ 对所有 $t \in [0,T]$ 成立（即温度有界），则二者在 $[0,T]$ 上恒等。该条件物理上自然成立，因热带大气中位温被限制在有限范围内（如 $200\,\mathrm{K} \leq T \leq 320\,\mathrm{K}$）。

推论 4.3（温度有界性）：对上述整体强解，存在仅依赖于初值 $\|u_0\|_{H^2}, \|v_0\|_{H^2}, \|T_0\|_{H^2}, \|q_0\|_{H^1}$ 的常数 $C > 0$，使得

该估计是证明唯一性的关键支撑，且本身具有独立物理意义——排除了数学解中出现非物理的无限升温现象。

所有结论均建立在全空间 $\mathbb{R}^2$ 上，无需周期性或有界区域假设；且湿度方程中严格不引入 $\Delta q$ 项，即 $\eta = 0$，这是与既有工作（如 [6] 中需 $\eta > 0$）的本质区别。

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

证明的核心障碍在于：多值不连续项 $H(q - q_s)$ 破坏了系统的单调性与 Lipschitz 连续性，而缺失 $\Delta q$ 项又剥夺了用抛物正则性压制其奇性的标准手段。作者为此构建了一个三层正则化-紧性框架，并发展出三项关键技术：

1. 双参数正则化策略与一致先验估计：首先引入两个正则化参数 $\epsilon > 0$（Heaviside 函数的 mollifier 宽度）和 $\eta > 0$（人工湿度粘性系数），构造正则化系统（见原文 Section 5）：

关键创新在于：对每个固定 $\epsilon, \eta > 0$，利用标准能量法结合 $\eta \Delta q$ 项提供的 $H^2$ 控制，可证得全局强解 $(u^{\epsilon,\eta}, v^{\epsilon,\eta}, T^{\epsilon,\eta}, q^{\epsilon,\eta})$ 存在且满足 $\|q^{\epsilon,\eta}(t)\|_{L^\infty} \leq 1$（由最大值原理保证 $q$ 保持在 $[0,1]$ 内）。更重要的是，作者证明了一组**$\eta$-一致先验估计**：$\|u^{\epsilon,\eta}\|_{L^\infty_t H^2_x}$, $\|v^{\epsilon,\eta}\|_{L^\infty_t H^2_x}$, $\|T^{\epsilon,\eta}\|_{L^\infty_t H^2_x}$, $\|q^{\epsilon,\eta}\|_{L^\infty_t H^1_x}$ 的上界均与 $\eta$ 无关（见原文 Section 6）。这使得后续 $\eta \to 0$ 极限可行，且极限函数 $q^\epsilon$ 自动继承 $H^1$ 正则性，为处理 $H(q-q_s)$ 提供了关键的 Sobolev 可微性基础。

2. 温度 $L^\infty$ 估计与 De Giorgi 迭代技术：为建立唯一性所需的温度有界性，作者并未依赖经典的 Moser 迭代（因其要求右端源项在 $L^p$ 中具有高阶可积性，而此处 $\mathrm{div}\, v$ 仅属 $L^2_t L^2_x$），而是将温度方程重写为

并观察到右侧第二项在 $T$ 较大时为负（因 $G^+(T) > 0$ 当 $T < \xi_0 = LR/(c_p R_\nu)$，而 $\xi_0$ 是典型热带温度上限），从而构成一个内在的“阻尼”效应。作者据此设计了一个精巧的 De Giorgi 迭代方案：定义截断函数 $T_k = (T - k)_+$，乘以方程并积分，利用 $G^+(T)$ 的符号性质导出形如

的微分不等式。再结合 Gronwall 不等式与迭代技巧，最终得到 $\|T(t)\|_{L^\infty} \leq C$，该常数仅依赖于初值范数与 $\|\mathrm{div}\, v\|_{L^2_t L^2_x}$ 的积分，后者由 $v$ 的 $H^2$ 范数控制。

3. 基于 Lemma 2.1 的唯一性分析框架：为证明两个强解的差 $(\delta u, \delta v, \delta T, \delta q)$ 恒为零，作者构造能量泛函 $f(t) = \|\delta u\|_{L^2}^2 + \|\delta v\|_{L^2}^2 + \|\delta T\|_{L^2}^2 + \|\delta q\|_{L^2}^2$。难点在于处理 $\delta q$ 方程中由 $H(q^{(1)} - q_s) - H(q^{(2)} - q_s)$ 引起的奇异性。作者利用 $H$ 的多值性及 $q^{(i)} \in [0,1]$ 的事实，将该项控制为 $|\delta q| \cdot \mathbf{1}_{\{|q^{(1)} - q_s| < \epsilon\} \cup \{|q^{(2)} - q_s| < \epsilon\}}$，再结合温度有界性 $\|T^{(i)}\|_{L^\infty} \leq C_T$ 所隐含的 $G^+(T^{(i)})$ 的 Lipschitz 性质，最终导出符合 Lemma 2.1 前提的微分不等式：

其中 $G(t) = \|\nabla \delta T\|_{L^2}^2 + \|\nabla \delta q\|_{L^2}^2$，$m_1 \in L^1_t$, $m_2 \in L^2_t$, $S \in L^1_t$。应用引理即得 $f \equiv 0$，从而唯一性得证。

4. 与之前工作的比较

本文与既有文献的推进主要体现在四个维度：

1. 相变建模的物理保真度与数学严格性：相较于 [10, 12] 中将降水项 $P$ 视为光滑函数或浅水近似，本文严格采用 $H(q - q_s)$ 多值结构，忠实反映了大气中水汽饱和阈值引发的突变物理。与 [6] 相比，[6] 在压力坐标系下使用微分包含证明弱解存在，但未获得强解；而本文在 $\mathbb{R}^2$ 上证明了强解的整体存在性，并给出了明确的正则性指标（$H^2$ for $u,v,T$; $H^1$ for $q$）。与 [14] 相比，[14] 的模型虽含 $H(q-q_s)$，但完全忽略了温度扩散（即 $-\Delta T = 0$），其解的正则性严重不足；本文则恢复了物理必需的 $-\Delta T$ 项，使温度方程成为标准抛物型，为 $L^\infty$ 估计提供了可能。

2. 湿度方程的无粘性设定与技术突破：[6] 和 [2] 均需在湿度方程中添加 $\eta \Delta q$ 项以保证正则性，其解的存在性依赖于 $\eta > 0$，且 $\eta \to 0$ 极限未被考察。本文则首次在不引入任何 $\Delta q$ 项的前提下，通过双参数正则化与 $\eta$-一致估计，成功实现了 $\eta \to 0$ 极限，从而获得了真正物理意义上的无粘性湿度输运模型的强解。这一代价是必须依赖温度有界性来启动唯一性论证，而 [6] 的粘性模型则可直接利用 $q$ 的 $H^2$ 正则性进行标准能量估计。

3. 解的适定性时空尺度：[15] 和 [20] 建立了无相变热带气候模型在 $\mathbb{R}^2$ 或有界区域上的整体强解适定性，但其时间尺度受限于初值正则性（如 $H^1$ 初值仅得局部解）。本文将此尺度推进至任意 $H^2 \times H^2 \times H^2 \times H^1$ 初值下的整体解，且空间域为全平面 $\mathbb{R}^2$，无需边界条件或周期性假设。这得益于对 $H(q-q_s)$ 奇异性的精细控制，而非简单回避。

4. 唯一性条件的物理可验证性：[16] 中关于原始方程弱解的唯一性通常依赖于 $\nabla u \in L^1_t L^\infty_x$ 等难以验证的条件；[14] 的唯一性要求 $\epsilon > 0$（convective relaxation parameter）且收敛速率为 $O(\sqrt{\epsilon})$。本文的唯一性条件 $\|T(t)\|_{L^\infty} \leq C_T$ 是物理上天然满足且可由观测数据直接验证的（热带大气位温范围已知），且不依赖于任何人为小参数。从提供内容中无法确认本文是否与 [26] 或 [27] 中的更精细热力学模型进行定量误差比较；作者亦未显式展开与 [1, 17] 关于小纵横比极限严格性的对比，故该维度暂不作断言。

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

1. 非零无穷远条件下方程 $u_t + i \Delta u + i |u|^2 u = 0$ 的黑孤子的轨道稳定性

Orbital stability of black solitons for quasilinear Schrödinger equations with nonzero conditions at infinity

作者: Erwan Le Quiniou

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13629v1

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

摘要翻译

中文标题: 非零无穷远条件下方程 $u_t + i \Delta u + i |u|^2 u = 0$ 的黑孤子的轨道稳定性

摘要翻译: 我们研究一类一维拟线性薛定谔方程在无穷远处具有非零边界条件情形下黑孤子的轨道稳定性。具体而言，我们的框架涵盖一般的散焦半线性非线性项，以及聚焦或散焦的拟线性非线性项。首先，我们建立关于拟线性非线性项的充分条件，以保证存在一条由其速度参数化的有限能量孤子局部分支；在此分支中，黑孤子（亦称扭结）对应于稳态解。我们的主要结果是：若满足 Vakhitov-Kolokolov (VK) 斜率条件（即动量关于速度的导数在零处为负），则黑孤子在能量空间中是轨道稳定的。此外，我们推导出一个用于验证该 VK 条件的显式公式。证明依赖于对一个精心设计的变分问题的分析，该分析使得我们能够用能量与动量（二者均由流保持守恒）来控制扭结扰动演化的上确界范数。论证中一个精细的部分在于对该变分问题极小化序列的分析，因为该变分问题的下确界并未达到。

2. 动力学方程正则性传递的一个新证明

A new proof of the transfer of regularity for kinetic equations

作者: Lukas Niebel

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13582v1

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

摘要翻译

中文标题: 动力学方程正则性传递的一个新证明

摘要翻译: 我们提出了一种基于轨迹的新方法，用于在局部扩散的弱尺度上建立动力学方程的正则性传递估计（à la Bouchut-Hörmander）。该方法避免了在傅里叶变量中的显式计算，也不依赖于基本解，但仍能导出尖锐的、尺度不变的齐次估计。

3. 不定权函数下对数拉普拉斯算子的谱性质

Spectral Properties of the Logarithmic Laplacian with Indefinite Weights

作者: Rakesh Arora, Tuhina Mukherjee, Arshi Vaishnavi

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13513v1

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

摘要翻译

中文标题: 不定权函数下对数拉普拉斯算子的谱性质

摘要翻译: 本文研究由带不定权函数的对数拉普拉斯算子驱动的加权特征值问题。我们证明了 Lusternik-Schnirelman 特征值存在一个无界序列，并且表明第一特征值是单重的，其对应的特征函数在定义域内具有恒定符号。相比之下，对应于更高阶特征值的特征函数必然变号。我们进一步建立了一类节点域型不等式，将高阶特征值与相应特征函数的正部和负部的测度联系起来，该不等式本身具有独立兴趣。作为应用，我们证明了第一特征值是孤立的。此外，我们获得了第一与第二特征值的替代性变分刻画，并建立了特征值关于权函数和定义域的单调性性质。

4. 与 $N$ Finsler 拉普拉斯算子相关的奇异 Liouville 方程解的分类

Classification of solutions to the singular Liouville’s equation associated with the $N$ Finsler Laplacian

作者: Jianwei Xue, Maochun Zhu

类别: math.AP, math.FA

发布日期: 2026-05-13

arXiv ID: 2605.13447v1

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

摘要翻译

中文标题: 与 $N$ Finsler 拉普拉斯算子相关的奇异 Liouville 方程解的分类

摘要翻译: 本文对任意 $β\in (0,N)$，对如下与 Finsler-$N$-拉普拉斯算子相关的奇异 Liouville 方程的一类解进行分类： \begin{align*} -\mathrm{div}\left(F^{N-1}(\nabla u)DF(\nabla u)\right)=\hat{F}^{o}(x)^{-β}e^u\ \ \text{in } \mathbb{R}^{N}\backslash {0}, \end{align*} 在有限质量条件 $\int_{\mathbb{R}^{N}}\hat{F}^{o}(x)^{-β}e^u dx<+\infty$ 下成立。此处 $F$ 是一个凸函数，且是 1 次正齐次函数；其极函数 $F^{o}$ 在 $\mathbb{R}^{N}$ 上定义了一个 Finsler 度量，而 $\hat{F}^{o}(x)=F^{o}(-x)$。我们的结果弱化了文献 [39] 中分类定理所要求的质量条件。

5. 利用Nash不等式证明带非局部Robin边界条件的热半群在 $\mathrm{L}^{2}\left( Ω\right)$ 中的超收缩性

Ultracontractivity of Heat semigroups in $\mathrm{L}^{2}\left( Ω\right)$ with non-local Robin boundary conditions using Nash’s inequality

作者: Christoph Schwerdt

类别: math.AP, math.FA

发布日期: 2026-05-13

arXiv ID: 2605.13413v1

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

摘要翻译

中文标题: 利用Nash不等式证明带非局部Robin边界条件的热半群在 $\mathrm{L}^{2}\left( Ω\right)$ 中的超收缩性

摘要翻译: 我们在 $\mathbb{R}^{d}$ 中的有界Lipschitz区域 $Ω$（其中 $d>2$）上研究热方程 $\frac{\partial u}{\partial t} - \operatorname{div} \left( A \nabla u \right) = 0$，其中 $-\operatorname{div} \left( A \nabla \cdot \right)$ 是一个具有广义Robin边界条件的二阶一致椭圆算子。这些边界条件形式上由 $ν\cdot A \nabla u + Bu = 0$ 给出，其中 $ν$ 是 $\partialΩ$ 上的外单位法向量，且 $B \in \mathcal{L} \left( \mathrm{L}^{2}\left( \partial Ω\right) \right)$ 是一个一般算子，允许破坏解半群的保正性。在关于 $B$ 的相当温和的假设下，我们利用Sobolev空间 $H^{1}( Ω)$ 上的Nash不等式证明了解半群的超收缩性。

6. p-莱维算子补值问题的最优稳定性

Optimal stability of complement value problems for p-Lévy operators

作者: Guy Foghem

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13389v1

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

摘要翻译

中文标题: p-莱维算子补值问题的最优稳定性

摘要翻译: 我们建立了由对称积分微分 p-莱维算子 $1 < p < \infty$ 所支配的积分微分方程（IDEs）在非局部 Dirichlet 或 Neumann 边界条件下解的最优收敛性。为说明起见，考虑（分数阶）p-拉普拉斯算子 $(-\Delta)^s_p$（其中 $0 < s \le 1$）的特例：若 $(-\Delta)^s_p u_s = f_s $ 在 $\Omega\subset \mathbb{R}^d$ 中成立，并辅以 Dirichlet 或 Neumann 数据 $g_s$，则在关于 $\Omega$、$f_s$ 和 $g_s$ 的适当假设下，我们证明 $(u_s)_s$ 在 $s \to 1^-$ 时强收敛于最优意义下，即 $\|u_s - u_1\|_{W^{s,p}(\Omega)} \to 0$。\smallskip 我们方法所支撑的另一后续目标是非局部迹空间的鲁棒性；具体而言，我们还证明了非局部迹空间以适当意义收敛于局部迹空间。

7. 燃烧型多孔介质方程中传播速度的精细估计

Refined estimates of the propagation speed in porous medium equation of combustion type

作者: Suying Liu, Fan Wu

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13361v1

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

摘要翻译

中文标题: 燃烧型多孔介质方程中传播速度的精细估计

摘要翻译: 我们研究初值问题 $u_{t}=(u^{m})_{xx}+f(u)$，其中非线性项 $f(u)$ 为燃烧型，且初始数据具有紧支集。在 \cite{lou2024convergence} 中，作者证明：对给定初始数据作缩放（即考虑其某个倍数），存在一个临界值，使得对应的过渡解以渐近速度 $2y_{0}\sqrt{t}[1+o(1)]\ \text{as} \ t\rightarrow\infty$ 传播，但低阶项 $o(1)$ 的具体形式仍未知。本文中，针对一族燃烧型函数，我们改进了过渡解渐近传播速度的估计，并对低阶项 $o(1)$ 给出了精确刻画。我们的结果还表明：对一般的燃烧型函数 $f$，并不存在统一的低阶项 $o(1)$ 刻画。

8. 分数阶Sobolev常数的Galerkin逼近

Galerkin Approximation of the Fractional Sobolev Constant

作者: Andreea Dima, Liviu I. Ignat

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

发布日期: 2026-05-13

arXiv ID: 2605.13347v1

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

摘要翻译

中文标题: 分数阶Sobolev常数的Galerkin逼近

摘要翻译: 我们在维数 $N\geq 1$、分数阶指数 $s\in (0,\min\{1,N/2\})$ 的情形下，建立了分数阶Sobolev不等式离散最优常数的尖锐估计。我们所建立的收敛速率适用于单位球上的分片线性元Galerkin逼近；在该区域上，我们采用拟一致且正则的网格进行计算。

9. 变指数 Muckenhoupt $A_\infty$ 权

Variable Muckenhoupt $A_\infty$ Weights

作者: Dachun Yang, Wen Yuan, Zongze Zeng

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

发布日期: 2026-05-13

arXiv ID: 2605.12941v1

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

摘要翻译

中文标题: 变指数 Muckenhoupt $A_\infty$ 权

摘要翻译: 本文通过引入变标量权 $\mathcal{A}_{p(\cdot),\infty}$ 与变矩阵权 $\mathscr{A}_{p(\cdot),\infty}$ 的概念，旨在建立变指数空间框架下 $A_\infty$ 权的完整理论。我们首先证明：一个权 $w$ 属于 $\mathcal{A}_{p(\cdot),\infty}$ 当且仅当其 $p(\cdot)$ 次幂 $w^{p(\cdot)}$ 是一个经典的 $A_\infty$ 权。利用该等价性，我们以极小算子刻画了 $\mathcal{A}_{p(\cdot),\infty}$ 条件。接着，我们在变指数 Lebesgue 空间中为 $\mathcal{A}_{p(\cdot),\infty}$ 权建立了带显式常数的反向 Hölder 不等式；并结合前述 $\mathcal{A}_{p(\cdot),\infty}$ 权与经典 $A_\infty$ 权之间的关系，证明：对任意权 $w$，其在变指数 Lebesgue 空间中满足反向 Hölder 不等式当且仅当 $w \in \mathcal{A}_{p(\cdot),\infty}$。对于矩阵权 $\mathscr{A}_{p(\cdot),\infty}$，我们首先证明了其约化算子的存在性；然后将矩阵权 $\mathscr{A}_{p(\cdot),\infty}$ 与标量权 $\mathcal{A}_{p(\cdot),\infty}$ 相结合，在变指数 Lebesgue 空间中建立了 $\mathscr{A}_{p(\cdot),\infty}$ 权的反向 Hölder 不等式。最后，为进一步应用于变矩阵加权函数空间，我们为 $\mathscr{A}_{p(\cdot),\infty}$ 权引入了上维数与下维数的概念，并利用这些概念建立了涉及约化算子的尖锐估计。

10. Boltzmann 方程在无穷远处 Mach 数为零的半空间问题

Half-space problem on the Boltzmann equation with zero Mach number at infinity

作者: Hongxu Chen, Jun-ling Chen, Renjun Duan

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.12914v1

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

摘要翻译

中文标题: Boltzmann 方程在无穷远处 Mach 数为零的半空间问题

摘要翻译: 我们在三维半空间 $ \mathbb{R}^2 \times \mathbb{R}^+ $ 中研究具有硬球碰撞的时变 Boltzmann 方程的长时间动力学行为，该方程满足漫反射边界条件，并围绕全局 Maxwellian 平衡态作小扰动。假设远场速度静止；即取无穷远处 Mach 数为零。第一个目标是构造关于 Maxwellian 的整体时间低正则性解。我们利用沿二维切向方向的时间衰减性质，建立与二维热方程匹配的解的多项式衰减速率。第二个目标是进一步证明 Gevrey 正则性的传播：在切向空间变量 $x_\parallel$ 上具有解析性（Gevrey 指标为 1），在切向速度变量 $v_\parallel$ 上属于 Gevrey 类（指标为 2），前提是初始数据具有适当的正则性。证明结合了基于傅里叶空间的 $L^1_k \cap L^p_k$ 方法以获得衰减估计、适用于无界区域的宏观-微观分解与 $L^2 - L^\infty$ 框架，以及加权 Gevrey 范数以控制正则性传播，从而克服了边界效应和非线性相互作用带来的困难。

11. 结构化连续最大正则性的推论

Implications of structured continuous maximal regularity

作者: Philip Preußler, Felix L. Schwenninger

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

发布日期: 2026-05-12

arXiv ID: 2605.12121v1

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

摘要翻译

中文标题: 结构化连续最大正则性的推论

摘要翻译: 我们研究了在空间范数本质上不同于上确界范数的情形下，关于连续函数的最大正则性估计如何自动改进。更准确地说，我们利用与底层线性发展方程的温和解相关的卷积型算子的弱紧性等性质，来强化先验估计。这些结果具有若干应用：例如，为Guerre-Delabriere关于 $\mathrm{L}^1$-最大正则性的结果提供了一个新证明，并推广了Baillon定理；简化了关于 $\mathrm{C}_0$-半群生成的著名扰动定理；并且我们解决了一个控制理论中关于输入到状态稳定性（input-to-state stability）的公开问题，该问题针对一大类抽象系统。

12. 气泡流约化模型推导的一种变分方法

A variational approach to the derivation of reduced models for bubbly flows

作者: Cosmin Burtea, David Gérard-Varet

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11938v1

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

摘要翻译

中文标题: 气泡流约化模型推导的一种变分方法

摘要翻译: 本文利用哈密顿最小作用原理，推导环境无粘液体中气泡运动的约化模型。我们首先说明如何从该原理恢复经典的尖锐界面模型，其中压力在气泡表面处连续。随后，我们通过简化这些表面的描述来降低模型复杂度：即限定界面演化于由有限个参数刻画的一类超曲面子类之中（最简单的例子是球面，亦即忽略气泡偏离球形的形变）。从数学与建模角度看，难点在于确定可替代压力连续性条件的界面条件。最后，在液体流动为无旋且气泡内压力均匀的情形下，我们通过适定性分析完成约化模型的推导。

13. 基于H-收敛的隐函数定理：非线性非光滑椭圆系统的均匀化

An H-convergence-based implicit function theorem for homogenization of nonlinear non-smooth elliptic systems

作者: Lutz Recke

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11976v1

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

摘要翻译

中文标题: 基于H-收敛的隐函数定理：非线性非光滑椭圆系统的均匀化

摘要翻译: 我们考虑具有非光滑数据的半线性椭圆系统Dirichlet问题的均匀化。我们假设当均匀化参数趋于零时，扩散张量满足H-收敛。我们的结果属于隐函数定理类型：对充分小的均匀化参数，存在唯一一个弱解，该弱解位于给定的非退化均匀化问题弱解附近。证明中，我们使用Meyers型梯度估计（当空间维数为二时）或Morrey型梯度估计（当扩散张量为三角形时）来处理线性椭圆系统的解。

14. Burgers方程的统一变换法：在有限区间非饱和流中的应用

The unified transform for Burgers’ equation: Application to unsaturated flow in finite interval

作者: Konstantinos Kalimeris, Leonidas Mindrinos, Athanasios Paraskevopoulos

类别: math.AP, math-ph

发布日期: 2026-05-12

arXiv ID: 2605.11788v1

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

摘要翻译

中文标题: Burgers方程的统一变换法：在有限区间非饱和流中的应用

摘要翻译: 本文聚焦于一维竖向入渗问题，假设扩散系数为常数，且水力传导度与含水量之间满足二次关系。在此假设下，Richards方程约化为Burgers方程，我们再通过Hopf-Cole变换将其线性化。由此，初边值问题转化为一个具有混合边界条件的有限区间上的扩散方程。为求解该问题，我们采用统一变换法（亦称Fokas方法）。该方法给出了解的显式积分表示；当进行数值计算时，所得结果与经典的Fourier级数解完全一致，但具有更优的收敛性与稳定性。文中考察了两个水文应用实例。

15. 非线性薛定谔流对非线性项的依赖性

Dependence of the nonlinear Schr{ö}dinger flow upon the nonlinearity

作者: Rémi Carles, Quentin Chauleur, Guillaume Ferriere

类别: math.AP, math-ph

发布日期: 2026-05-12

arXiv ID: 2605.11643v1

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

摘要翻译

中文标题: 非线性薛定谔流对非线性项的依赖性

摘要翻译: 我们考虑能量次临界情形下的散焦非线性薛定谔方程，并研究解对非线性项幂次的依赖性。特别关注其整体时间描述。本文所处理的三个主要方面按难度递减排序如下：(1) 当总幂次趋于 1 时的极限，以及该极限与对数型薛定谔方程 $i\partial_t u + \Delta u = u \log|u|^2$ 的联系；(2) 可能存在长程效应的情形下的描述；(3) 短程情形下散射算子的连续性。本文整理自第一作者于 2026 年 5 月在巴黎综合理工学院为 Laurent Schwartz 讨论班所作的报告。

16. 移动的局部观测与保守型偏微分方程的Cesàro渐近可观测性

Moving localized observations and Ces{à}ro asymptotic observability for conservative PDEs

作者: Maarten V. de Hoop, Antti Kykkänen, Emmanuel Trélat

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11624v1

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

摘要翻译

中文标题: 移动的局部观测与保守型偏微分方程的Cesàro渐近可观测性

摘要翻译: 我们构建了一种确定性的大时间机制，该机制可从保守演化系统的移动局部观测中导出Cesàro渐近可观测性不等式。在每个观测时间区间上，紧致的、具有测度的齐性空间上的精确凸化操作，将整个观测流形上的完全观测替换为一个原型子集的有限个平移的凸组合。随后，一个切换实现定理（switching realization theorem）将该静态设计转化为真正运动的观测器；而一个Hilbert空间框架下的尾部约化命题（Hilbertian tail-reduction proposition）则表明：即使区间估计仅在增长的谱窗（growing spectral windows）上被证明，经Cesàro平均后仍可恢复全部守恒能量。该设计到可观测性的完整链条适用于如下情形：(i) 紧致的、具有测度的齐性流形上波动方程、Klein-Gordon方程与Schrödinger方程的内部观测；(ii) 欧氏球体上的移动边界帽（moving boundary caps）；(iii) 一类奇异的、近乎分离的气态巨行星边界模型（singular almost-separated gas-giant boundary model）。当每个瞬时观测集过小、以致无法期望成立有限时间GCC（Geometric Control Condition）或依赖于时间的GCC陈述时，该框架尤为相关。

17. 慢-快随机非线性分数阶薛定谔方程的平均化原理

Averaging principle for a slow-fast stochastic nonlinear fractional Schrödinger equation

作者: Manil T. Mohan, Debopriya Mukherjee, Sandip Roy

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11561v1

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

摘要翻译

中文标题: 慢-快随机非线性分数阶薛定谔方程的平均化原理

摘要翻译: 我们在一维环面上建立了一个由乘性Wiener噪声驱动的结构多尺度随机非线性分数阶薛定谔系统的平均化原理。慢变量分量由带有一般多项式非线性的分数阶薛定谔算子所支配，而快变量分量在更短的时间尺度上演化，并表现出耗散扩散、非线性相互作用及随机强迫。 在适当的耗散性假设下，我们证明：当尺度分离参数趋于零时，慢变量分量强收敛于一个有效的随机分数阶薛定谔方程 $i\partial_t u + (-\Delta)^\alpha u = \overline{F}(u) + \sigma(u)\dot{W}$，其中有效漂移项 $\overline{F}(u)$ 通过对冻结快动力学的唯一不变测度取耦合项的平均而得到。 证明依赖于一致先验估计、快方程的遍历性、通过粘性消失法（vanishing viscosity method）获得的慢变量分量的Hölder时间正则性，以及一种适配于分数阶色散算子的Khasminskii型时间离散化论证。由于分数阶薛定谔半群的光滑性有限，且存在一般多项式非线性，该分析在技术上具有挑战性；我们通过精细估计与粘性逼近予以处理。

不相关论文（基本信息）

1. 网络上非自治 Hamilton-Jacobi 方程的 Lax-Oleinik 公式

Lax-Oleinik formula for nonautonomous Hamilton-Jacobi equations on networks

作者: Marco Pozza

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13704v1

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

摘要: We provide a Lax-Oleinik-type representation formula for solutions to nonautonomous Hamilton-Jacobi equations posed on networks with a rather general geometry. The networks may possess countably many arcs and allow for the presence of loops. We consider Hamiltonians that are convex and superlinear in the momentum variable, and satisfy a Lipschitz-type condition in the time variable. The representation formula is constructed via an overall Lagrangian that accounts for both the arc-specific dynamics and vertex-based constraints, called flux limiters, which ensure the well-posedness of the problem. We prove that the corresponding action functional admits Lipschitz continuous minimizers without needing to rule out the Zeno phenomenon. Furthermore, we demonstrate that the formula yields the unique solution to the problem even when the flux limiters exceed standard upper bounds.

2. 临界双曲悬链面在 $\mathbb{H}^3$ 中的 Medvedev Morse 指标猜想的解析局部解

Analytic local resolution of Medvedev’s Morse index conjecture for the critical hyperbolic catenoid in $\mathbb{H}^3$

作者: Alexander Pigazzini

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

发布日期: 2026-05-13

arXiv ID: 2605.13562v1

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

摘要: Let $Σ_a\subset B^3(r(a))\subset\mathbb{H}^3$ ($a>1/2$) be the critical hyperbolic catenoid of the Mori family, a free boundary minimal surface in the geodesic ball. The Medvedev conjecture [8] states ind$(Σ_a)=4$ for all $a>1/2$. We study its strong form: ind$(Σ_a)=4$ and nul$(Σ_a)=2$. The nullity condition nul$(Σ_a)=2$ combines the mode-$|k|=1$ result $\text{nul}_R(Σ_a)|_{|k|=1}=2$ of [10, Cor. 4.4] with vanishing kernel in modes $|k|=0,|k|\ge2$; the latter, not in [10], is established here for $a\in(1/2,1/2+δ_0)$. The main result is the analytic local resolution of the strong Medvedev conjecture: $\existsδ_0>0$ s.t. ind$(Σ_a)=4$, nul$(Σ_a)=2$ for all $a\in(1/2,1/2+δ_0)$. This follows from the expansion $H(a):=\sinh r(a)/K(a)=σ_*\coshσ_*+C_0(a-\frac12)+O((a-\frac12)^2)$ as $a\to(1/2)^+$, with $C_0=\frac{σ_*\coshσ_*(\sinh^2σ_*-1)(3\sinh^2σ_*-2)}{12\sinh^2σ_*}$, where $σ_*>0$ the unique positive root of $σ=\cothσ$, and $C_0>0$ by $σ_*>\log(1+\sqrt2)$. The proof proceeds via three reductions: $(i)$ the Medvedev conjecture is equivalent to $μ_0^{\mathrm{even}}(2)>0$ $(E)$ and $μ_2(0)>0$ with non-degeneracy in mode $0$ $(F)$; $(ii)$ $μ_2(0)>0$ reduces, via a Sturm shooting-count argument, to $φ_a>0$ of the parametric Jacobi field on the principal branch; $(iii)$ $φ_a>0$ reduces, under $\sinh r(a)>2K(a)$ $(G)$, to $H'(a)>0$ via a constant Wronskian and Sturm separation. Auxiliary results: a Picone identity (base $f_*$) closing unconditionally the odd radial sector for $|k|\ge2$; a second Picone identity (base $B$) proving $(E)$ unconditionally on $(1/2,1]$ and, via Hardy estimates, on $(1/2,A_*]$ ($A_*>1$); analytic closure of $(G)$ on $(1/2,1]$ via strict concavity of a transcendental function; an alternative proof of ind$(Σ_a)\ge4$ via Lorentz ambient coordinates.

3. 非线性Hartree系统存在无穷多个多峰解

Infinitely many multi-peaks solutions for a nonlinear Hartree system

作者: Qihan He, Qingfang Wang

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13531v1

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

摘要: In this paper, we study the following nonlinear Hartree system: $-Δu_i + V_i(x)u_i = μ_i φ_{u_i}u_i + \sum_{j\neq i}β_{ij}φ_{u_j}u_i$ for $x\in\mathbb{R}^3$, with $u_i\in H^1(\mathbb{R}^3)$ ($i=1,2,3$), where $φ_u(x):=\int_{\mathbb{R}^3}\frac{u^2(y)}{|x-y|}dy$ for any $u\in H^1(\mathbb{R}^3)$, $V_i(x)$ ($i=1,2,3$) are continuous bounded radial functions, and $β_{ij}$ are coupling constants. We mainly investigate the effects of the potentials and the nonlinear coupling terms on the structure of solutions. Applying the Lyapunov-Schmidt reduction method, we prove the existence of infinitely many solutions to the system. Specifically, the solutions we obtain satisfy that some components are synchronized with each other but segregated from the others, and that some components are positive while others are sign-changing. To the best of our knowledge, it is the first time that solutions possessing some positive components and some sign-changing ones have been constructed using Lyapunov-Schmidt reduction methods. Moreover, it is also the first attempt to investigate systems consisting of three Hartree equations with mixed couplings.

4. 部分观测随机过程的结构可识别性：从单粒子轨迹到总粒子密度数据

Structural identifiability of partially-observed stochastic processes: from single-particle trajectories to total particle density data

作者: Arianna Ceccarelli, Alexander P. Browning, Ruth E. Baker

类别: stat.ME, math.AP, math.DS, math.PR, q-bio.QM

发布日期: 2026-05-13

arXiv ID: 2605.13504v1

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

摘要: The increasing availability of experimental data has intensified interest in calibrating stochastic models, raising fundamental questions about parameter identifiability. Structural identifiability determines whether parameters can be uniquely recovered from idealised, noise-free data, a prerequisite to allow for parameter estimation. However, existing methods to assess structural identifiability are not generally applicable to stochastic processes. We develop a methodology to analyse structural identifiability for a class of spatio-temporal stochastic processes. We investigate how identifiability depends on the type of available data, distinguishing between single-particle trajectories and total particle density measurements. For trajectory data, we use the individual-based model description that explicitly represents single-particle dynamics. For population-level data, we derive a partial differential equation model representation, that describes the evolution of total particle density, and apply a differential algebra approach, common to ordinary differential equations analysis. We further introduce a novel method to study the initial condition, based on characteristic equations to construct a Taylor expansion of the density evolution, enabling identification of additional identifiable parameter combinations. We apply our methodology to a model, and show it is identifiable with trajectory data but only locally identifiable with density data, and demonstrate the critical role of initial conditions in the identifiability analysis.

5. 具变系数的时间分数阶拟抛物方程的正向与反向问题

Forward and inverse problems for a time-fractional pseudo-parabolic equation with variable coefficients

作者: Ravshan Ashurov, Elbek Husanov

类别: math.AP

发布日期: 2026-05-13

arXiv ID: 2605.13285v1

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

摘要: In this work, forward and inverse problems for a time-fractional pseudo-parabolic equation $D_t^ρ [u(t) + μAu(t)] + σ(t) Au(t) = r(t)g$ are investigated in a Hilbert space, where $A$ is an unbounded, positive, self-adjoint operator. According to the known papers, the forward problem has been studied only in the case $σ(t) = const$. The main novelty of the forward problem in this work is that the model is further generalized and investigated for a time-dependent coefficient $σ(t)$. To determine the solution of the forward problem, the Fourier method is employed, and the global existence and uniqueness of the solution are proved. Moreover, when the operator $A$ is a second-order differential operator, a numerical scheme and an efficient computational algorithm are developed. The inverse problem of determining a time-dependent source function is considered under the overdetermination condition of the form $F[u(t)] = Φ(t)$. The functional $F$ is taken in a general form, and such an inverse problem has not been considered before. The global existence of the solution to the inverse problem is proved by applying Schauder’s fixed point theorem, and its uniqueness is established. Furthermore, several examples related to the operator $A$ and the functional $F$ are provided.

6. 物理信息神经网络的统一广义性分析

Unified generalization analysis for physics informed neural networks

作者: Yuka Hashimoto, Tomoharu Iwata

类别: cs.LG, math.AP, math.FA, stat.ML

发布日期: 2026-05-13

arXiv ID: 2605.13260v1

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

摘要: Physics-Informed Neural Networks (PINNs) and their variational counterparts (VPINNs) are neural networks that incorporate physical laws, making them useful for scientific problems. Existing generalization analyses for PINNs and VPINNs remain limited, often requiring restrictive assumptions such as stability conditions or linear ellipticity. In this paper, we derive generalization bounds for neural networks that involve differentiation with respect to input variables, covering PINNs and VPINNs under a unified framework. We apply Taylor expansion to represent nonlinear differential operators as linear operators on a high-dimensional space, enabling the use of Koopman-based analysis and showing that high-rank networks can generalize well even in settings involving differential operators. We also show that the nonlinearity of the differential operator exponentially enlarges the bound, highlighting its significant impact on generalization.

7. Nesterov 加速法用于位移凸自由能的 Wasserstein 最小化

Nesterov acceleration for the Wasserstein minimization of displacement-convex free energies

作者: Pierre Monmarché

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

发布日期: 2026-05-13

arXiv ID: 2605.13186v1

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

摘要: We show that the mean-field underdamped Langevin process (associated to the non-linear Vlasov-Fokker-Planck equation) achieves a Nesterov acceleration with respect to the Wasserstein gradient flow of a displacement-convex free energy, in the sense that it converges at a rate of order given by the square-root of the Polyak-Łojasiewicz constant of the free energy (which is the optimal convergence rate for the corresponding gradient flow). This result has been made possible by the recent breakthrough [42] by Jianfeng Lu, which establishes such a \emph{diffusive-to-ballistic} improvement in term of entropy in the linear case.

8. 环面上极小叶状结构的极小-极大间隙刻画

A min-max gap characterization of minimal foliations on the torus

作者: Hoan Nguyen

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

发布日期: 2026-05-12

arXiv ID: 2605.12428v1

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

摘要: We extend an energy introduced by Mather to the setting of Almgren-Pitts min-max theory and obtain a parametric, higher-dimensional analogue of Mather’s variational barrier theory for twist maps and geodesics on tori. We use this energy to establish several criteria for the existence of foliations of the $n$-torus by minimal hypersurfaces. We show that for a generic metric, whenever a lamination by area-minimizing hypersurfaces of the $n$-torus contains a gap, there exists a minimal hypersurface inside the gap that is not area-minimizing. As an application, we derive a recurrence property for totally irrational minimal foliations.

9. 逆双调和散射问题中扩展采样方法的新实现

Novel implementation of the extended sampling method for inverse biharmonic scattering

作者: Isaac Harris, General Ozochiawaeze

类别: math.AP, math.NA

发布日期: 2026-05-12

arXiv ID: 2605.12367v1

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

摘要: This paper considers an inverse shape problem for recovering an unknown clamped obstacle in two dimensions from far–field measurements generated by a single incident wave or just a few incident waves for the biharmonic (flexural) wave equation. Here we will develop a new extended sampling method (ESM) that is derived using the analysis of the well–known factorization method. We will also consider an ESM using both sound–soft and sound–hard sampling disks to identify sampling points where the reference disk intersects the unknown cavity. The use of a sound–hard sampling disk has not been studied in the literature whereas the sound–soft sampling disk has been used in most recent works. Traditionally the ESM seeks to find the location of the scatterer from limited incident directional data. Here, our method acts more like the factorization method to obtain the location as well as the size (and possibly the shape) of the obstacle. We present numerical experiments with synthetic data that demonstrate how effective this new implementation is with respect to noisy data and illustrate the influence of the reference disk radius on the reconstruction.

10. 关于利用可导型拉普拉斯变换求解扩散方程 $\partial_t u = \Delta u$ 的解

On solution of Diffusion Equation using Conformable Laplace Transform

作者: Somnath Sarate, Anil Khairnar, Krishnat Masalkar

类别: math.DS, math.AP, math.FA, math.OA

发布日期: 2026-05-12

arXiv ID: 2605.12157v1

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

摘要: The inversion theorem and convolution theorem of the conformable fractional Laplace transforms are developed. All the elementary properties of the classical Laplace transform are extended to the conformable fractional transform, and using these properties, we found analytical solutions to the initial-boundary value problems of the diffusion equation.

11. 基于测度源对流-扩散模型的泄漏定位

Leak localisation with a measure source convection-diffusion model

作者: Thi Tam Dang, Tuomo Valkonen

类别: math.AP, math.OC

发布日期: 2026-05-12

arXiv ID: 2605.12095v1

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

摘要: We study the inverse problem of locating gas leaks from line-of-sight concentration measurements using a convection-diffusion model with the source term a Radon measure. By imposing sparsity-promoting regularisation on this measure, we recover point sources - identifying both their locations and intensities - rather than diffuse approximations. We jointly estimate the underlying physical convection (wind) and diffusion parameters. Our main theoretical contribution is the stability analysis of the convection-diffusion equation with respect to its parameters: the measure, and the convection and diffusion fields. Numerically, we employ a semi-grid-free optimisation approach for reconstructing the source measure. Our experiments demonstrate accurate localisation, highlighting the potential of the method for practical gas emission detection.

12. 关于分数阶对数 $p$-拉普拉斯算子

On the fractional logarithmic $p$-Laplacian

作者: Anouar Bahrouni, Abdelhamid Gouasmia, Hichem Hajaiej, Anass Ouannasser

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.12050v1

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

摘要: In this paper, we introduce and investigate the fractional logarithmic $p$-Laplacian $(-Δ)_{p}^{s+\log}$, defined as the first-order derivative with respect to the parameter $t$ of the fractional $p$-Laplacian $(-Δ)_{p}^{t}$ evaluated at $t=s$. We establish that this operator admits the following integral representation

where $C(N,s,p)$ denotes the standard normalization constant associated with the fractional $p$-Laplacian, and $B(N,s,p)=\frac{d}{ds}\left(\ln C(N,s,p)\right)$. As a consequence of this representation, it follows that the operator is nonlocal and of logarithmic type, and may be viewed as a nonlinear analogue of the fractional logarithmic Laplace operator recently introduced by Chen et al. \cite{Chen-Chen-Hauer}. We further develop the associated functional framework in both $\mathbb{R}^{N}$ and bounded Lipschitz domains by introducing the natural energy spaces adapted to problems driven by $(-Δ)_{p}^{s+\log}$. Within this framework, fundamental functional inequalities are established, in particular Pohozaev-type identities and Díaz-Saa inequalities, which are of independent interest and applicable to a broader class of problems. Moreover, we derive results concerning density, continuity, and compact embedding properties. We emphasize that the compactness of the embedding is proved at the critical exponent $p^{*}_{s}=\frac{Np}{N-sp}$, which distinguishes the present setting from the classical Sobolev and fractional Sobolev frameworks. Finally, as an application, we investigate the associated Dirichlet eigenvalue problem and derive existence, uniqueness, and boundedness results for the corresponding solutions.

13. 来自 Ginzburg-Landau 模型的一类分支输运问题的尖锐上界

Sharp upper bound for a branched transport problem coming from Ginzburg-Landau models

作者: Alessandro Cosenza, Michael Goldman, Felix Otto

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11834v1

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

摘要: We consider a branched transport type problem with weakly imposed boundary conditions, which can be seen as a blown-up version of a reduced model for type-I superconductors in the regime of vanishing external magnetic field. We prove that if the irrigated measure is (locally) Ahlfors regular then it is of dimension at most $8/5$ in agreement with the conjecture by Conti, the third author and Serfaty.

14. 互联区域中的退化反应-扩散 SIR 模型

A degenerate reaction-diffusion SIR model in interconnected regions

作者: Omar Elamraoui, Jawad Salhi, Abderrahim Zafrar

类别: math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11925v1

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

摘要: This paper presents a novel time-space SIR (Susceptible-Infected-Recovered) model for simulating infectious disease dynamics in two interconnected regions. The model is formulated as a coupled reaction-diffusion system with boundary conditions that dynamically switch from Robin to Neumann types, effectively modelling policy-driven interventions such as lockdowns. A key innovation lies in the incorporation of degenerate diffusion, arising from vanishing population density, which significantly influences transmission behaviour near regional borders. The wellposedness of the model is rigorously established using the Faedo-Galerkin method, ensuring the existence, uniqueness, and positivity of weak solutions. Numerical simulations, performed using the Finite Volume Method, validate the theoretical findings and demonstrate the impact of migration and mobility restrictions on epidemic progression. This framework offers valuable insights for understanding and controlling disease spread in spatially heterogeneous and interconnected settings.

15. 假设流体密度非常数的牛顿问题

The Newton’s problem assuming non-constant density of the fluid

作者: Rafael López

类别: math.AP, math.DG

发布日期: 2026-05-12

arXiv ID: 2605.11641v1

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

摘要: This paper investigates the Newton’s problem of minimal resistance for a body moving through a fluid whose density decreases exponentially with altitude. We prove the local existence and regularity of radial solutions $u(r)$ satisfying the initial conditions $u(0)=u'(0)=0$ using a fixed-point theorem. We show that the maximal domain of the solution is finite, $[0, r_M)$, terminating at a critical slope $u'(r_M) = \frac{1}{\sqrt{3}}$.

16. 在谱Ricci下界与平均凸边界的约束下的刚性与柔性

Rigidity and flexibility under spectral Ricci lower bounds and mean-convex boundary

作者: Gioacchino Antonelli, Yangyang Li, Paul Sweeney

类别: math.DG, math.AP

发布日期: 2026-05-12

arXiv ID: 2605.11384v1

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

摘要: We study Riemannian manifolds $(M^n,g)$ with mean-convex boundary whose Ricci curvature is nonnegative in a spectral sense. Our first main result is a sharp spectral extension of a rigidity theorem by Kasue: we prove that under the conditions

and in the sharp range $0\leq γ<4$ if $n=2$, and $0\leqγ<\frac{n-1}{n-2}$ if $n\geq3$, a (possibly noncompact) complete manifold with disconnected boundary, with at least one compact boundary component, must split isometrically as a product $[0,L]\times Σ$. Our second main contribution is a topological rigidity result for the relative fundamental group $π_1(M,\partial M)$, combined with a deep theorem of Lawson–Michelsohn. We prove that, in dimensions $n\neq4$, any compact manifold with boundary satisfying the two inequalities above, with at least one of them strict, admits a metric with positive sectional curvature and strictly mean-convex boundary, provided $γ\geq0$ if $n=2$, and $0\leqγ\leq\frac{n-1}{n-2}$ if $n\geq3$. This range of $γ$ is sharp for the latter result to hold.