生成时间: 2026年07月02日 06:27:34

运行状态: 完成

计划处理论文数量: 23

已完成论文数量: 23

未完成论文数量: 0

重点关注论文数量: 8

了解领域论文数量: 5

不相关论文数量: 10

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

1. n维球坐标系下梯度与散度的显式公式

Explicit formulas for gradients and the divergence in n-dimensional spherical coordinates

作者: Bern Rummler, Gudrun Thäter

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.01096v1

摘要翻译: 我们利用n维球坐标系（n>1）下的拉普拉斯算子，在向量微积分框架下写出定义在径向对称区域上的向量场的散度。我们应用向量微积分的基本公式，结合从笛卡尔坐标系到n维球极坐标系的nabla算子及坐标变换矩阵。例如，为证明向量场是n维环形区域与球体上Stokes算子的特征函数，就需要该向量场的散度。我们在n维球极坐标系下以偏导数形式给出的散度公式，是未来在这些区域上进一步验证Stokes特征函数的重要一步。

Comment: 无

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

1. 研究对象和背景

本文研究的核心对象是 $\mathbb{R}^n$ （ $n \geq 2$ ）中径向对称区域 $\Omega^{(n)}$ 上向量场的散度算子 $\operatorname{div}$ 与梯度算子 $\nabla$ 在 $$n$$ 维球坐标系下的显式表达式。具体而言，设 $\Omega^{(n)}$ 为开球 $\{x \in \mathbb{R}^n : \|x\| < R\}$ 或环形域 $\{x \in \mathbb{R}^n : R_i < \|x\| < R_o\}$ ，其边界具有旋转对称性。研究动机直接源于不可压 Navier-Stokes 方程和 Stokes 算子本征值问题的精确求解：在分析 Stokes 算子 $S = -\Upsilon \Delta_x$ 的本征函数时，需同时满足两个关键条件：(i) $\operatorname{div} v = 0$ （无源性），(ii) $-\Delta_x v + \nabla p = \lambda v$ （动量方程）。由于物理上自然出现的候选本征函数（如 toroidal 场、poloidal 场）常以球坐标分量形式给出——例如 $v = v_r e_r + v_{\vartheta_1} e_{\vartheta_1} + \cdots + v_\varphi e_\varphi$ ，其中 $v_r, v_{\vartheta_k}, v_\varphi$ 是 $r, \vartheta_1, \dots, \vartheta_{n-2}, \varphi$ 的标量函数——因此必须将散度算子本身表示为这些变量的偏微分组合，才能验证 $\operatorname{div} v = 0$ 。然而，标准数学物理教材（如 Moon & Spencer [4]）仅系统给出 $$n=2,3$$ 的散度公式；对 $n \geq 4$ ，文献中缺乏统一、可直接用于分析的显式公式。作者明确指出：在前期工作 [11] 中证明 $$n$$ 维环形域上 Poincaré 常数的精确值时，即因缺失 $$n>3$$ 的散度公式而被迫绕行，这构成一个实质性技术空白。更深层的理论需求在于：Stokes 算子 $$S$$ 定义为 Leray-Helmholtz 投影 $\Upsilon$ 与拉普拉斯算子 $-\Delta_x$ 的复合，其定义域 $D(S) = W^2_2(\Omega^{(n)}) \cap V(\Omega^{(n)})$ 要求向量场 $$v$$ 属于 Sobolev 空间且散度为零；而 $V(\Omega^{(n)})$ 的刻画本质上依赖于散度算子在所选坐标系下的精确微分结构。若使用球坐标表示 $$v$$ ，但散度公式未知，则无法严格验证 $v \in V(\Omega^{(n)})$ ，进而阻碍对 $$S$$ 本征函数基的构造与正交完备性（见 Corollary 8）的严格证明。此外，该问题亦关联到 Laplace-Beltrami 算子 $$B$$ 在单位球面 $\omega^{(n)} = \{x \in \mathbb{R}^n : \|x\| = 1\}$ 上的谱理论：由 Remark 16 与 Definition 17， $$n$$ 维拉普拉斯算子可分解为径向部分与球面部分：

\Delta_{\text{sph}} = \frac{1}{r^{n-1}} \frac{\partial}{\partial r}\left(r^{n-1} \frac{\partial}{\partial r}\right) - \frac{1}{r^2} B,

其中 $$B$$ 的本征值为 $\ell(\ell + n - 2)$ （Theorem 19）。而 $$B$$ 的显式微分表达式（Definition 17）正是从 $$n$$ 维球坐标散度与梯度的几何推导中自然导出的，故本文工作亦为高维球面调和分析提供基础工具。问题的重要性不仅在于计算便利性，更在于它是连接欧氏空间向量分析、Stokes 算子谱理论与高维球面几何分析的关键桥梁。

2. 主要定理或主要结果

本文的核心成果是给出了 $$n$$ 维球坐标系下散度算子 $\operatorname{div}$ 的统一显式公式（公式 (7)）及梯度算子 $\nabla$ 的完整表达式（Remark 10）。具体而言，在 $n \geq 3$ 维空间中，采用标准球坐标参数化（Definition 12）： $x_1 = r \sin \vartheta_1 \cdots \sin \vartheta_{n-2} \cos \varphi$ , $x_2 = r \sin \vartheta_1 \cdots \sin \vartheta_{n-2} \sin \varphi$ , …, $x_n = r \cos \vartheta_1$ ，对应的单位正交基为 $e_r, e_{\vartheta_1}, \dots, e_{\vartheta_{n-2}}, e_\varphi$ 。对任意向量场 $v = v_r e_r + \sum_{k=1}^{n-2} v_{\vartheta_k} e_{\vartheta_k} + v_\varphi e_\varphi$ ，其散度为：

\operatorname{div} v = \frac{1}{r^{n-1}} \frac{\partial}{\partial r}(r^{n-1} v_r) + \frac{1}{r \sin^{n-2} \vartheta_1} \frac{\partial}{\partial \vartheta_1}(\sin^{n-2} \vartheta_1 \, v_{\vartheta_1}) + \sum_{k=2}^{n-3} \frac{1}{r \sin \vartheta_1 \cdots \sin \vartheta_{k-1} \sin^{n-k-1} \vartheta_k} \frac{\partial}{\partial \vartheta_k}(\sin^{n-k-1} \vartheta_k \, v_{\vartheta_k}) + \frac{1}{r \sin \vartheta_1 \cdots \sin \vartheta_{n-2}} \left[ \frac{\partial}{\partial \vartheta_{n-2}}(\sin \vartheta_{n-2} \, v_{\vartheta_{n-2}}) + \frac{\partial v_\varphi}{\partial \varphi} \right].

此公式对所有整数 $n \geq 2$ 成立，且已通过 $$n=2,3,4$$ 的特例（公式 (4)-(6)）验证其一致性。梯度算子作用于标量函数 $\tilde{v}(r, \vartheta_1, \dots, \vartheta_{n-2}, \varphi)$ 的结果为：

\nabla \tilde{v} = \frac{\partial \tilde{v}}{\partial r} e_r + \sum_{k=1}^{n-2} \left( \prod_{j=1}^{k-1} \sin^{-1} \vartheta_j \right) \frac{1}{r} \frac{\partial \tilde{v}}{\partial \vartheta_k} e_{\vartheta_k} + \left( \prod_{j=1}^{n-2} \sin^{-1} \vartheta_j \right) \frac{1}{r} \frac{\partial \tilde{v}}{\partial \varphi} e_\varphi.

文中强调，该公式是“straightforward standard calculations”的结果（Remark 11），其推导基于度量张量 $g = J J^T = \operatorname{diag}\{g_{jj}\}$ 的显式计算（Remark 14），其中 $g_{11} = 1$ , $g_{22} = r^2$ , …, $g_{nn} = r^2 \sin^2 \vartheta_1 \cdots \sin^2 \vartheta_{n-2}$ ，并利用正交坐标系下散度的通用公式 $\operatorname{div} v = \frac{1}{\sqrt{\det g}} \sum_{i=1}^n \frac{\partial}{\partial u_i} \left( \sqrt{\det g} \, v^i \right)$ ，其中 $(u_1, \dots, u_n) = (r, \vartheta_1, \dots, \vartheta_{n-2}, \varphi)$ 为坐标， $$v^i$$ 为逆变分量。文中未声明任何存在性、唯一性或正则性定理，其结论纯粹是坐标变换下的恒等式，适用于所有 $$C^1$$ 类向量场与标量函数。

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

本文的方法论核心并非发展新的分析技术，而是对经典微分几何与张量分析工具进行系统性、高维推广与显式化实现。主要技术障碍在于：当维度 $$n$$ 一般化后，球坐标系的 Jacobian 矩阵 $$J$$ （Remark 14）与度量张量 $$g$$ 的结构变得高度非平凡，其行列式 $\det g$ 及其平方根 $\sqrt{\det g}$ 的表达式需被精确捕捉，否则无法代入正交坐标系下散度的通用公式。作者的真正创新点在于：(1) 对任意 $n \geq 3$ ，完整、无歧义地写出 $$n$$ 维球坐标参数化（Definition 12）及其 Jacobian 矩阵 $$J$$ 的通式，并据此严格推导出 $g = \operatorname{diag}\{g_{jj}\}$ 的对角元素序列 $g_{11}=1, g_{22}=r^2, \dots, g_{nn}=r^2 \sin^2 \vartheta_1 \cdots \sin^2 \vartheta_{n-2}$ ；(2) 基于此，显式计算 $\sqrt{\det g} = r^{n-1} \sin^{n-2} \vartheta_1 \sin^{n-3} \vartheta_2 \cdots \sin^2 \vartheta_{n-3} \sin \vartheta_{n-2}$ ，这是公式 (7) 中所有系数项的共同来源；(3) 将此 $\sqrt{\det g}$ 代入通用散度公式，并针对每个坐标方向 $$u_i$$ 进行逐项求导，最终整理出具有清晰递归结构的显式表达式。这一过程依赖于对三角函数幂次的系统性归纳：例如， $\vartheta_k$ 方向的系数总包含 $\sin^{n-k-1} \vartheta_k$ 的因子，这源于 $\sqrt{\det g}$ 中 $\sin \vartheta_k$ 的幂次为 $$n-k-1$$ 。文中还利用了梯度与散度的内在关系 $\Delta \tilde{v} = \operatorname{div}(\nabla \tilde{v})$ （公式 (3)）作为交叉验证手段，即先用 $$J$$ 和 $$g$$ 推出梯度（Remark 10），再对其应用已得的散度公式，应复现标准 $$n$$ 维拉普拉斯算子 $\Delta_{\text{sph}}$ （Remark 16）。此外，作者构建了从笛卡尔坐标到球坐标的完整线性变换矩阵 $T_{c,s}$ （Remark 15），其 $$n=2,3,4$$ 的具体形式（公式 (11)-(13)）及 $$n$$ 维通式（公式 (14)）为向量分量的坐标变换提供了严格框架，确保 $v_c = T_{c,s} v_s$ 的一致性，从而保证了 $\operatorname{div} v$ 作为标量场的坐标不变性。整个证明框架不涉及任何估计、紧性或变分论证，纯粹是微分几何层面的恒等式推演，其严谨性建立在对 $$J$$ 、 $$g$$ 及其代数性质的精确操控之上。

4. 与之前工作的比较

本文与既有工作的根本区别在于维度普适性与公式显式性的双重突破。经典文献如 Moon & Spencer [4]、Arfken & Weber 等仅详述 $$n=2,3$$ 的球坐标散度公式（即公式 (4)-(5)），对 $n \geq 4$ 仅给出模糊的‘模式类似’提示，缺乏可供严格分析使用的闭式表达。本文首次为任意有限维 $n \geq 2$ 提供了完全显式的、可直接代入计算的散度公式 (7)，其价值在于消除了高维情形下的‘黑箱’依赖。与 Rummler 等人前期工作 [11] 相比，本文是其技术前提的补全：[11] 的核心结论—— $$n$$ 维环形域上 Poincaré 常数的精确值——依赖于对特定向量场 $$v$$ 的散度为零的验证，而该验证在 [11] 中因缺乏 $$n>3$$ 的公式而未完成严格证明；本文则为此提供了不可或缺的工具。与 Laplace-Beltrami 算子 $$B$$ 的已有理论（如 Triebel [13]）相比，本文的贡献在于将 $$B$$ 的抽象定义（Definition 17）落地为具体的、可操作的微分算子，其系数完全由本文推导的 $\sqrt{\det g}$ 决定，从而使得 Theorem 19 中 $$B$$ 的谱性质（本征值 $\ell(\ell+n-2)$ ）能在具体计算中被显式验证。然而，本文亦有明确局限：其结论是纯代数恒等式，不涉及任何函数空间上的有界性、连续性或紧性估计；它未讨论奇点（如 $$r=0$$ 或 $\vartheta_k=0,\pi$ ）处的正则性问题，亦未给出在 $W^{1,p}$ 或 $$H^1$$ 空间中散度算子的延拓；其适用范围严格限定于光滑的径向对称区域 $\Omega^{(n)}$ 及其内部，对带角点或非光滑边界的区域无效。此外，文中未将此公式应用于新的存在性或唯一性定理的证明，其定位是‘prerequisite’（前提）而非‘main theorem’（主定理），正如摘要所言，它是‘a future verification of further Stokes eigenfunctions’ 的基石。因此，本文的推进是基础性、工具性的，其代价是牺牲了分析深度，换取了维度上的彻底普适与公式的绝对明晰。

2. 低正则性先验估计、爆破准则与自由边界理想磁流体动力学中方的自交奇点（含表面张力）

Low-regularity a priori estimates, blow-up criterion, and self-intersection singularities for free-boundary ideal magnetohydrodynamics with surface tension

作者: Tao Luo, Siqi Yang

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00891v1

摘要翻译: 我们研究三维不可压自由边界理想磁流体动力学（MHD）方程组，其中包含表面张力且自由边界为闭曲面。我们的第一个结果在一般有界区域中建立了 $$H^3$$ 先验估计，不要求图结构、周期性或单连通性；特别地，对于一般区域中带表面张力的理想 MHD，该结果将此前已知的正则性阈值从 $$H^6$$ 降低至 $$H^3$$ 。与不可压 Euler 方程的自由边界问题相比，其特征在于洛伦兹力参与椭圆型压力估计，且冻结磁场必须保持切向边界约束。利用这些估计，我们证明了一个精细的有限时间爆破准则，适用于 $$H^3$$ 解，该准则将拓扑自交、边界正则性丧失、法向速度爆破以及内部 MHD 爆破彼此分离。其中内部条件具有内在的磁场不对称性：除 $

Comment: 48 pages, 10 figures

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

1. 研究对象和背景

本文研究三维不可压自由边界理想磁流体动力学（MHD）方程组，其数学模型为：

\begin{cases} D_t u + \nabla p = h \cdot \nabla h, & \text{in } \Omega_t, \\ D_t h = h \cdot \nabla u, & \text{in } \Omega_t, \\ \nabla \cdot u = 0,\quad \nabla \cdot h = 0, & \text{in } \Omega_t, \\ u_n = U,\quad h_n = 0,\quad p = A, & \text{on } \partial \Omega_t, \\ u(\cdot, 0) = u_0,\quad h(\cdot, 0) = h_0, & \text{in } \Omega_0, \end{cases}

其中 $D_t := \partial_t + u \cdot \nabla$ 为物质导数； $u(x,t) \in \mathbb{R}^3$ 为速度场； $h(x,t) \in \mathbb{R}^3$ 为磁场； $p = p_{\text{fluid}} + \frac{1}{2}|h|^2$ 为总压力（含磁压项）； $\Omega_t \subset \mathbb{R}^3$ 是随时间演化的有界区域，其边界 $\partial \Omega_t$ 为闭曲面； $$ n(x,t) $$ 是单位外法向量； $A = \nabla \cdot n$ 是平均曲率（表面张力系数归一化为 1）； $u_n := u \cdot n$ ， $h_n := h \cdot n$ 。系统要求 $$ h $$ 在边界上切向（即 $$ h_n = 0 $$ ），且在内部满足 $\nabla \cdot h = 0$ ，体现“冻结磁场”（frozen-in field）结构。

该模型是 Euler 方程与 Maxwell 方程在理想（无粘、无电阻）极限下的耦合，其核心困难在于三重非线性耦合：(i) 自由边界运动由法向速度 $$ u_n $$ 驱动，而 $$ u_n $$ 又受压力 $$ p $$ 控制；(ii) 压力 $$ p $$ 满足椭圆边值问题 $-\Delta p = \nabla \cdot (h \cdot \nabla h - u \cdot \nabla u)$ ，其边界条件 $p|_{\partial \Omega_t} = A$ 将几何（曲率）与分析（压力）强耦合；(iii) 磁场演化由输运型方程 $D_t h = h \cdot \nabla u$ 描述，同时必须严格维持 $h \cdot n = 0$ 的切向约束——这一约束在自由边界发生自交时失效，因接触点处法向量不唯一，导致 $h \cdot n = 0$ 无法定义。因此，理解自交奇点是否可被表面张力或理想磁耦合抑制，本质上是检验该模型能否维持全局嵌入性（global embeddedness）的拓扑刚性问题。

已有理论中，带表面张力的自由边界 MHD 的适定性长期受限于高正则性假设：Hao–Yang [23] 在一般有界域上建立 $$ H^6 $$ 先验估计，Liu–Luo [28] 在单连通等限制下获得 $H^{9/2}$ 局部适定性，而图表示（graph）设定下 [32] 仅达 $H^{7/2}$ 。相比之下，Euler 方程已实现 $$ H^3 $$ 适定性 [37]。故 $$ H^3 $$ 水平对 MHD 是否可行，是关键开放问题。此外，爆破准则方面，[23] 给出 $$ H^6 $$ BKM 型准则，但其未分离自交与内爆，且当 $h \to 0$ 时不退化为 Euler 准则；而自交奇点构造在理想 MHD 中完全空白——此前仅见于二维粘性 MHD [22, 24] 或等离子体-真空问题 [5]，均依赖耗散或低维简化。本文旨在填补三项根本空白：(1) 在无图结构、无周期性、无单连通性假设的一般有界域上建立 $$ H^3 $$ 先验估计；(2) 构造能分离四类爆破机制（自交、边界正则性丧失、法向速度爆破、内爆）的精细准则；(3) 首次证明三维理想 MHD 带表面张力仍可产生有限时间自交奇点，从而否定表面张力或理想磁耦合具有拓扑阻断能力。

2. 主要定理或主要结果

定理 1.1（ $$ H^3 $$ 先验估计）：设 $(u,h,\Omega_t)$ 是方程组 (1.1) 在 $$ [0,T) $$ 上的解，定义能量泛函

E(t) := \|D_t^2 u\|_{L^2(\Omega_t)}^2 + \|D_t^2 h\|_{L^2(\Omega_t)}^2 + \|D_t u\|_{H^{3/2}(\Omega_t)}^2 + \|D_t h\|_{H^{3/2}(\Omega_t)}^2 + \|u\|_{H^3(\Omega_t)}^2 + \|h\|_{H^3(\Omega_t)}^2 + \|\nabla(D_t u \cdot n)\|_{L^2(\partial \Omega_t)}^2 + 1,

并设几何控制量 $P_T := R - \sup_{t\in[0,T)} \|\eta(\cdot,t)\|_{L^\infty(\Gamma)}$ （ $$ R $$ 为初始边界的内外球半径），分析控制量

Q_T := \sup_{t\in[0,T)} \left( \|\nabla u\|_{L^\infty(\Omega_t)} + \|\nabla h\|_{L^\infty(\Omega_t)} + \|\nabla^2 h\|_{L^2(\Omega_t)} + \|\eta(\cdot,t)\|_{H^{5/2}(\Gamma)} + \|u_n\|_{H^2(\partial \Omega_t)} \right).

若存在常数 $C^\dagger > 0$ 使得 $P_T^{-1} + Q_T < C^\dagger$ ，则成立

\sup_{t\in[0,T)} \left( E(t) + \|D_t p\|_{H^1(\Omega_t)} + \|B\|_{H^2(\partial \Omega_t)} \right) \leq C(C^\dagger) e^{C(C^\dagger)(1+T)} E(0),

其中 $C(C^\dagger)$ 依赖于 $C^\dagger$ ， $$ E(0) $$ 仅依赖于 $\|u_0\|_{H^3(\Omega_0)}$ 、 $\|h_0\|_{H^3(\Omega_0)}$ 和 $\|A\|_{H^2(\partial \Omega_0)}$ 。进一步，存在 $$ C_0 > 0 $$ 和 $$ T_0 > 0 $$ （仅依赖初值），使得该先验假设在 $$ [0,T_0) $$ 上成立。

定理 1.2（有限时间爆破准则）：设 $(u,h,\Omega_t)$ 是 (1.1) 的最大存在解，存在区间 $[0,T_{\max})$ 满足 $(u,h) \in C([0,T_{\max}); H^3(\Omega_t) \times H^3(\Omega_t))$ 且 $\partial \Omega_t \in C([0,T_{\max}); H^4)$ 。若 $T_{\max} < \infty$ ，则对任意充分小的 $\varepsilon > 0$ ，以下至少一者发生：(1) 自由边界 $\partial \Omega_t$ 在 $t = T_{\max}$ 自交；(2) 边界正则性丧失： $\limsup_{t \nearrow T_{\max}} \left( \|A\|_{H^{1/2}(\partial \Omega_t)} + \|\partial \Omega_t\|_{H^{2+\varepsilon}} \right) = \infty$ ；(3) 法向速度爆破： $\limsup_{t \nearrow T_{\max}} \|u_n\|_{H^2(\partial \Omega_t)} = \infty$ ；(4) 内部 MHD 爆破： $\limsup_{t \nearrow T_{\max}} \left( \|\nabla u\|_{L^\infty(\Omega_t)} + \|\nabla h\|_{L^\infty(\Omega_t)} + \|\nabla^2 h\|_{L^2(\Omega_t)} \right) = \infty$ 。此准则的关键不对称性在于：速度仅需 $L^\infty$ 控制，而磁场额外需要 $$ L^2 $$ 型二阶导数控制 $\|\nabla^2 h\|_{L^2}$ ，该量源于洛伦兹力对压力估计的贡献，且无速度对应项。

定理 1.3（自交奇点存在性）：存在初值 $(u_0,h_0,\Omega_0) \in H^3 \times H^3 \times \{\text{有界域}\}$ ，满足 $\nabla \cdot u_0 = 0$ 、 $\nabla \cdot h_0 = 0$ 、 $h_0 \cdot n_0 = 0$ ，使得对应解在 $[0,T_{\max})$ 上保持 $$ H^3 $$ 正则性，且 $\partial \Omega_t$ 在 $t = T_{\max}$ 发生有限时间拓扑自交；同时， $\|u\|_{H^3(\Omega_t)}$ 、 $\|h\|_{H^3(\Omega_t)}$ 、 $\|A\|_{H^2(\partial \Omega_t)}$ 在 $[0,T_{\max})$ 上一致有界。

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

本文的核心技术障碍在于：(i) 表面张力将压力与边界几何强耦合，导致椭圆压力估计与边界正则性估计必须同步封闭；(ii) 理想 MHD 的冻结磁场要求 $h \cdot n = 0$ 在边界上恒成立，而自交时法向量不唯一，该约束失效；(iii) $$ H^3 $$ 框架下，标准能量估计无法直接控制 $\|\nabla^2 h\|_{L^2}$ ，但该量却是压力估计所必需的——这构成 MHD 特有的结构性不对称。

作者的真正创新点体现在三层次分析框架：

第一，发展 MHD 专属的低正则压力估计。经典 Euler 方法中，压力由 $-\Delta p = \nabla \cdot (u \cdot \nabla u)$ 确定，其 $$ H^1 $$ 估计依赖 $\|u\|_{H^2}$ 。而在 MHD 中，右端变为 $\nabla \cdot (h \cdot \nabla h - u \cdot \nabla u)$ ，其 $$ L^2 $$ 范数需 $\|h \cdot \nabla h\|_{H^1}$ ，进而要求 $\|\nabla h\|_{L^\infty} \|\nabla h\|_{H^1} + \|h\|_{H^2} \|\nabla^2 h\|_{L^2}$ 。后者强制引入 $\|\nabla^2 h\|_{L^2}$ ，且无速度类比项。为此，作者设计非对称先验假设 $$ Q_T $$ ，在 $$ H^3 $$ 框架下显式分离速度与磁场的控制需求，并通过切向微分算子 $\nabla$ （定义为 $\nabla f - (\nabla f \cdot n)n$ ）精确刻画边界几何量（如第二基本形式 $B = \nabla n$ 、平均曲率 $A = \nabla \cdot n$ ）的正则性传递。

第二，构建反向时间紧致性论证以处理自交奇点。由于 $$ p = A $$ 在自交点无定义，作者采用一族光滑逼近域 $\Omega_\lambda$ （ $\lambda > 0$ 为两边界分支间距），在每个 $\Omega_\lambda$ 上施加标准表面张力条件。关键突破在于建立 $\lambda$ -一致的 Hölder 型边界速度控制：利用 $H^{3/2}$ 内部控制导出

\left\| v_n^\lambda \big(x + \eta^\lambda(x,t)n^\lambda(x,0), t\big) - v_n^\lambda(x,0) \right\|_{L^\infty_x(\partial \Xi_0^\lambda)} \leq C t^\alpha, \quad \alpha \in (0,1),

此估计避免了对 $L^\infty$ 边界加速度的依赖，从而获得与 $\lambda$ 无关的统一存在时间 $$ T^* > 0 $$ ，保障紧致性所需的共同区间。

第三，设计兼容冻结磁场的奇异态配置。在自交域 $\Omega_\ddagger$ 上，作者构造轴对称奇异速度场（具严格向内法向分量）与切向磁场，再通过局部坐标分解将其光滑延拓至 $\Omega_\lambda$ 。关键洞察在于：自交点处两分支的单侧法向量共线，故其切平面重合，磁场只需在此公共切平面上切向，无需指定相对取向（Remark 1.4）。这使 $h \cdot n = 0$ 约束在逼近过程中得以一致维持，且 $\|h\|_{H^3(\Omega_\lambda)}$ 保持 $\lambda$ -一致有界，从而保证洛伦兹力 $h \cdot \nabla h$ 在能量估计中可控。

4. 与之前工作的比较

本文在三个维度显著推进既有工作：

正则性门槛方面：此前一般有界域上的表面张力理想 MHD 仅有 Hao–Yang [23] 的 $$ H^6 $$ 先验估计，而本文降至 $$ H^3 $$ ，降幅达三阶 Sobolev 指标。该改进并非单纯技术优化，而是揭示了 MHD 结构的本质约束—— $\|\nabla^2 h\|_{L^2}$ 的引入使 $$ H^3 $$ 成为临界阈值：更低阶（如 $$ H^2 $$ ）无法闭合压力估计，更高阶（如 $$ H^6 $$ ）则过度消耗正则性冗余。对比 Euler 方程的 $$ H^3 $$ 理论 [37]，本文代价在于必须接受速度与磁场控制的不对称性（ $\|\nabla u\|_{L^\infty}$ vs $\|\nabla h\|_{L^\infty} + \|\nabla^2 h\|_{L^2}$ ），而 Euler 中二者对称。

爆破准则精细化方面：[23] 的 $$ H^6 $$ 准则未区分自交与内爆，且当 $h \to 0$ 时不退化为 Euler 准则；本文定理 1.2 首次将四类机制（自交、边界正则性丧失、法向速度爆破、内爆）纳入统一框架，并明确 $\|\nabla^2 h\|_{L^2}$ 为 MHD 特有爆破指标。与 Euler 准则 [21] 对比，本文在边界项上放宽了正则性要求（ $\|A\|_{H^{1/2}}$ vs $\|A\|_{H^{3/2}}$ ， $\|u_n\|_{H^2}$ vs $\|u_n\|_{H^{5/2}}$ ），但以强化内爆条件为代价（ $L^\infty_t L^\infty_x$ 控制 $\|\nabla u\|_{L^\infty}$ ，而非 $L^1_t L^\infty_x$ 积分型条件）。此 trade-off 表明：MHD 的磁场非线性迫使边界与内部控制更紧密耦合。

自交奇点构造方面：此前所有理想 MHD 自交结果均受限于维度或模型——[5] 限于二维等离子体-真空问题且仅得解析性瞬时丧失，[22,24] 依赖粘性正则化。本文首次在三维、无粘、无电阻、带表面张力设定下构造自交，其核心突破在于证明表面张力不具拓扑阻断能力：即使 $$ p = A $$ 在极限点失效，只要在逼近域上一致控制曲率 $\|A\|_{H^2}$ 和磁场 $\|h\|_{H^3}$ ，自交仍可发生。此结论否定了表面张力作为“几何稳定器”的直观猜想，揭示理想磁耦合与表面张力在拓扑层面的相容性。代价是构造高度依赖轴对称性与特定初值设计，尚未推广至一般情形；且未解决自交后解的延拓问题（即弱解定义），此为后续工作方向。

3. Chern-Simons 系统爆破解的局部唯一性与非退化性

Local Uniqueness and Non-degeneracy of Blow Up Solutions To A Chern-Simons System

作者: Zetao Cheng, Haoyu Li, Lei Zhang

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.01154v1

摘要翻译: 本文研究一类重要的 Chern-Simons 系统的爆破解。我们首先证明：当发生平均场型爆破（mean-field type blowup）时，在自然的几何假设下，相应的爆破解是唯一的。我们还建立了该类爆破解附近线性化系统的非退化性（non-degeneracy）。为证明这些主要结果，我们进行了精确的爆破分析，使得解的渐近描述揭示出唯一性与非退化性结论所需的曲率信息。与相关文献中类似问题的研究相比，我们的估计更为精细且技术上更为复杂。

Comment: 35 pages

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

1. 研究对象和背景

本文研究定义在单位面积平环面 $\Omega = \mathbb{R}/\mathbb{Z} \times \mathbb{R}/\mathbb{Z}$ 上的耦合Chern-Simons系统：

\begin{cases} \Delta u_1 + \frac{1}{\varepsilon^2} e^{u_2 + u_{0,2}} (1 - e^{u_1 + u_{0,1}}) = 4\pi N & \text{in } \Omega, \\ \Delta u_2 + \frac{1}{\varepsilon^2} e^{u_1 + u_{0,1}} (1 - e^{u_2 + u_{0,2}}) = 4\pi N & \text{in } \Omega, \\ u_1, u_2 \in H^1_{\mathrm{per}}(\Omega), \end{cases} \tag{1.5}

其中 $u_{0,i}(x) = -4\pi \sum_{j=1}^N G(x, p_{i,j})$ ， $$G(x,p)$$ 是环面上满足 $-\Delta_x G(x,p) = \delta_p - 1$ 且 $\int_\Omega G(x,p)\,dx = 0$ 的Green函数； $\{p_{1,j}, p_{2,j}\}_{j=1}^N \subset \Omega$ 是固定的奇异源点。该系统源于非阿贝尔规范场论、多层超导与相互作用玻色凝聚体模型，其物理意义在于描述带拓扑Chern-Simons项的双组分涡旋场——区别于经典Abel-Higgs模型中的Maxwell项，Chern-Simons项引入了非局部、高度非线性的指数型耦合机制，导致方程具有强非对称性与内在耦合结构。

数学上，当 $\varepsilon \to 0^+$ 时，系统出现“爆破”（blowup）现象：能量密度在有限个点 $\{q_1,\dots,q_k\} \subset \Omega$ 处集中，即存在序列 $\{x_{j,\varepsilon}\}_\varepsilon$ 满足 $\lim_{\varepsilon\to 0} x_{j,\varepsilon} = q_j$ 且 $u_{i,\varepsilon}(x_{j,\varepsilon}) + 2\ln \frac{1}{\varepsilon} \to +\infty$ ；而在任意紧支集 $K \Subset \Omega \setminus \{q_1,\dots,q_k\}$ 上， $u_{i,\varepsilon} + 2\ln \frac{1}{\varepsilon} \to -\infty$ 。此类解称为“平均场型爆破解”（mean-field type blowup solutions），其局部轮廓由全平面Liouville系统

\begin{cases} \Delta U_1 + e^{u_{0,2}(q_j) + U_2} = 0 & \text{in } \mathbb{R}^2, \\ \Delta U_2 + e^{u_{0,1}(q_j) + U_1} = 0 & \text{in } \mathbb{R}^2, \end{cases} \tag{2.6}

的径向对称解刻画；全局行为则通过Green函数 $$G$$ 及正则部分 $\gamma(x,y)$ （即 $G(x,y) = -\frac{1}{2\pi} \ln |x-y| + \gamma(x,y)$ ）编码几何与拓扑相互作用。爆破分析的核心目标是建立精确渐近展开：对每个 $$j$$ ，定义尺度 $\mu_{j,\varepsilon} = \varepsilon^{-1} e^{\frac{1}{2} \beta_{j,\varepsilon}}$ ，其中 $\beta_{j,\varepsilon} = \max\{u_{1,\varepsilon}(x_{1,j,\varepsilon}), u_{2,\varepsilon}(x_{2,j,\varepsilon})\}$ ，并设 $U^*_{i,j,\varepsilon}(x) = V_{i,j,\varepsilon}(\mu_{j,\varepsilon}^{-1}(x - x^*_{j,\varepsilon}))$ ，其中 $(V_{1,j,\varepsilon}, V_{2,j,\varepsilon})$ 是(2.6)的适当平移解，则误差项 $\eta_{i,j,\varepsilon}(x) = u_{i,\varepsilon}(x) - \beta_{j,\varepsilon} - U^*_{i,j,\varepsilon}(x) - \widetilde{M}_{i,j,\varepsilon} (\gamma(x,x_{j,\varepsilon}) - \gamma(x_{j,\varepsilon},x_{j,\varepsilon})) - \sum_{l \ne j} \widetilde{M}_{i,l,\varepsilon} (G(x,x_{l,\varepsilon}) - G(x_{j,\varepsilon},x_{l,\varepsilon}))$ 必须被控制至最优阶。已有工作（Huang–Zhang [16]）建立了初步估计，但其精度不足以支撑唯一性与非退化性证明：例如，局部质量 $\widetilde{M}_{i,j,\varepsilon} = \varepsilon^{-2} \int_{\Omega_j} e^{u_{i',\varepsilon} + u_{0,i'}} (1 - e^{u_{i,\varepsilon} + u_{0,i}}) dx$ 仅被控为 $|\widetilde{M}_{i,j,\varepsilon} - 8\pi| = O(\mu_{j,\varepsilon}^{-1+\theta} + e^{\beta_{j,\varepsilon}/2})$ ，而唯一性要求 $O(\mu_{j,\varepsilon}^{-2+\theta})$ ；梯度估计 $|\nabla f_{i,j,\varepsilon}(x_{j,\varepsilon})| = O(\mu_{j,\varepsilon}^{-1} + e^{\beta_{j,\varepsilon}})$ 亦需提升至 $O(\mu_{j,\varepsilon}^{-2})$ ；更关键的是，参数关系 $\varepsilon$ 与 $\beta_{j,\varepsilon}$ 的联系仅达 $O(\mu_{j,\varepsilon}^{-2+\theta})$ 阶，而Pohozaev恒等式要求 $O(\mu_{j,\varepsilon}^{-4+\theta})$ 高阶展开。本文旨在填补这一空白：在自然几何假设下，严格证明爆破解的局部唯一性与线性化算子的非退化性，从而为Lyapunov–Schmidt约化、涡旋动力学及Morse指标计算提供坚实基础。

2. 主要定理或主要结果

定理1.1（局部唯一性）：设 $u^{(1)}_\varepsilon = (u^{(1)}_{1,\varepsilon}, u^{(1)}_{2,\varepsilon})$ 与 $u^{(2)}_\varepsilon = (u^{(2)}_{1,\varepsilon}, u^{(2)}_{2,\varepsilon})$ 是问题(1.1)的两列满足假设(A)的爆破解，其公共爆破点集为 $\{q_1,\dots,q_k\}$ 。若以下条件成立：(A) $\vec{q} = (q_1,\dots,q_k)$ 是泛函 $$G^*_1 + G^*_2$$ 的非退化临界点，其中 $G^*_i(\vec{q}) = \sum_{j=1}^k u_{0,i}(q_j) + 8\pi \sum_{1\le j<l\le k} G(q_j,q_l)$ ；(B) 全局量 $D(\vec{q}) < 0$ ，其中 $D(\vec{q}) = \lim_{\delta\to 0^+} \big[ \sum_{j=1}^k \rho_{1,j} e^{u_{0,1}(q_j)} \big( \int_{\Omega_j} \frac{e^{f_{1,j}} - 1}{|x-q_j|^4} dx - \int_{\mathbb{R}^2 \setminus \Omega_j} \frac{1}{|x-q_j|^4} dx \big) + \sum_{j=1}^k \rho_{2,j} e^{u_{0,2}(q_j)} \big( \int_{\Omega_j} \frac{e^{f_{2,j}} - 1}{|x-q_j|^4} dx - \int_{\mathbb{R}^2 \setminus \Omega_j} \frac{1}{|x-q_j|^4} dx \big) \big]$ ， $f_{i,j}(x) = 8\pi \big( \gamma(x,q_j) - \gamma(q_j,q_j) + \sum_{l\ne j} (G(x,q_l) - G(q_j,q_l)) \big) + u_{0,i}(x) - u_{0,i}(q_j)$ ， $\rho_{i,j} = e^{8\pi(\gamma(q_j,q_j)+\sum_{l\ne j} G(q_l,q_j)) + u_{0,i}(q_j)}$ ；(C) 对任意 $j_1,j_2 = 1,\dots,k$ ，有 $u_{0,1}(q_{j_1}) - u_{0,2}(q_{j_1}) = u_{0,1}(q_{j_2}) - u_{0,2}(q_{j_2})$ ，则存在 $\varepsilon_0 > 0$ ，使得对所有 $\varepsilon \in (0,\varepsilon_0)$ ，有 $u^{(1)}_{1,\varepsilon} \equiv u^{(2)}_{1,\varepsilon}$ 且 $u^{(1)}_{2,\varepsilon} \equiv u^{(2)}_{2,\varepsilon}$ 。

定理1.2（非退化性）：在定理1.1的全部假设下，存在 $\varepsilon_1 > 0$ ，使得对任意 $\varepsilon \in (0,\varepsilon_1)$ ，爆破解 $(u_{1,\varepsilon}, u_{2,\varepsilon})$ 是非退化的，即线性化问题

\begin{cases} \Delta \phi_1 + \frac{1}{\varepsilon^2} \big[ e^{u_2 + u_{0,2}} \phi_2 - e^{u_1 + u_2 + u_{0,1} + u_{0,2}} (\phi_1 + \phi_2) \big] = 0 & \text{in } \Omega, \\ \Delta \phi_2 + \frac{1}{\varepsilon^2} \big[ e^{u_1 + u_{0,1}} \phi_1 - e^{u_1 + u_2 + u_{0,1} + u_{0,2}} (\phi_1 + \phi_2) \big] = 0 & \text{in } \Omega, \\ \phi_1, \phi_2 \in H^1_{\mathrm{per}}(\Omega) \end{cases} \tag{1.9}

仅有平凡解 $(\phi_1,\phi_2) = (0,0)$ 。

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

本文的核心技术障碍在于：耦合Chern-Simons系统的强非线性、双组分不对称性及环面曲率效应共同导致传统爆破分析工具失效。作者的真正创新点在于构建了一套“逐阶精炼的爆破分析框架”，其关键步骤如下：

第一，完全爆破行为的严格推导（Section 2.2）。不同于[16]将(B1)(B5)(B6)作为先验假设，本文从假设(A)出发，通过反证法结合Liouville系统分类定理（Theorem A）与Harnack/Liouville定理，严格证明：若某一分量爆破而另一分量不爆破，则极限函数 $$U_1$$ 必为常数，与积分有界性矛盾，从而确立 $\lim_{\varepsilon\to 0} x_{1,j,\varepsilon} = \lim_{\varepsilon\to 0} x_{2,j,\varepsilon} = q_j$ 及 $|u_{1,\varepsilon}(x_{1,j,\varepsilon}) - u_{2,\varepsilon}(x_{2,j,\varepsilon})| = O(1)$ ，即(B1)(B5)自动成立。此为后续所有估计的基石。

第二，局部质量的尖锐估计（Lemma 2.14）。利用改进的Pohozaev恒等式与精细的Green函数展开，作者将质量偏差从 $O(\mu_{j,\varepsilon}^{-1+\theta} + e^{\beta_{j,\varepsilon}/2})$ 提升至 $O(\mu_{j,\varepsilon}^{-2+\theta})$ 。关键在于引入“几乎等质量分割引理”（Lemma 2.10）：证明 $|\widetilde{M}_{1,j,\varepsilon} - \widetilde{M}_{2,j,\varepsilon}| = O(\mu_{j,\varepsilon}^{-2+\theta})$ ，从而消去主导误差项，使质量收敛速率翻倍。

第三，正则部分梯度的最优控制（Lemma 2.16）。作者发现 $\nabla f_{i,j,\varepsilon}(x_{j,\varepsilon})$ 的主部来自 $\sum_{l\ne j} \widetilde{M}_{i,l,\varepsilon} \nabla G(x_{j,\varepsilon}, x_{l,\varepsilon})$ ，而由前述质量估计及 $|x_{j,\varepsilon} - q_j| = O(\mu_{j,\varepsilon}^{-2})$ （由(B6)的验证），可得 $|\nabla f_{i,j,\varepsilon}(x_{j,\varepsilon})| = O(\mu_{j,\varepsilon}^{-2})$ ，达到与标量情形[20]相同的精度。

第四， $\varepsilon$ - $\beta$ 关系的四阶展开（Proposition 2.2）。这是唯一性证明的枢纽。作者将Pohozaev恒等式在爆破区域边界 $\partial B_{\theta_{j,\varepsilon}}(x_{j,\varepsilon})$ 上积分，利用前述所有尖锐估计（特别是 $O(\mu_{j,\varepsilon}^{-2+\theta})$ 质量与 $O(\mu_{j,\varepsilon}^{-2})$ 梯度），得到

D(\vec{q}) e^{u_{0,1}(q_j)} e^{-8\pi \gamma(x_{j,\varepsilon},x_{j,\varepsilon}) - 8\pi \sum_{l\ne j} G(x_{j,\varepsilon},x_{l,\varepsilon})} \mu_{j,\varepsilon}^{-2} + B e^{8\pi \gamma(x_{j,\varepsilon},x_{j,\varepsilon}) + 8\pi \sum_{l\ne j} G(x_{j,\varepsilon},x_{l,\varepsilon})} e^{\beta_{j,\varepsilon}} = O(\mu_{j,\varepsilon}^{-4+\theta}),

其中 $$B > 0$$ 为普适常数。此四阶精度足以消除Pohozaev恒等式中所有低阶项，最终导出 $u^{(1)}_\varepsilon - u^{(2)}_\varepsilon \equiv 0$ 。

第五，非退化性的谱分析（Section 4.3）。作者将线性化算子 $\mathcal{L}_\varepsilon$ 分解为内、外区域贡献：在 $B_\delta(q_j)$ 内，利用Liouville系统线性化核的完备性（Theorem B）及前述尖锐估计，证明任何非平凡零模必与核元素 $\partial_{x_1} U_j, \partial_{x_2} U_j, \frac{d}{d\delta}U_j(\delta x)|_{\delta=1} + 2$ 线性相关；在外区域，利用Green函数展开与 $O(\mu_{j,\varepsilon}^{-2})$ 误差控制，证明零模在 $\Omega \setminus \cup_j B_\delta(q_j)$ 上衰减至零。二者结合，得出 $\ker \mathcal{L}_\varepsilon = \{0\}$ 。

4. 与之前工作的比较

本文相较既有工作实现了三重实质性推进：

首先，相较于Huang–Zhang [16]，本文在假设强度上显著放宽：[16]将(B1)(B5)(B6)作为独立假设强加于解序列，而本文仅需更自然的几何假设(A)，并严格推导出(B1)(B5)（Corollary 2.2），且(B6)的验证也由 $O(\mu_{j,\varepsilon}^{-2})$ 估计自动蕴含。在技术精度上，[16]的局部质量估计为 $O(\mu_{j,\varepsilon}^{-1+\theta})$ ，本文提升至 $O(\mu_{j,\varepsilon}^{-2+\theta})$ ；其梯度估计为 $O(\mu_{j,\varepsilon}^{-1})$ ，本文达 $O(\mu_{j,\varepsilon}^{-2})$ ；其 $\varepsilon$ - $\beta$ 关系仅至 $O(\mu_{j,\varepsilon}^{-2+\theta})$ ，本文达 $O(\mu_{j,\varepsilon}^{-4+\theta})$ 。这些提升并非简单迭代，而是依赖于“几乎等质量分割”与“完全爆破推导”等原创技术。

其次，相较于标量Chern-Simons方程的唯一性工作（Lin–Yan [20]），本文首次将局部唯一性理论推广至耦合系统。[20]处理单方程 $\Delta u + \varepsilon^{-2} e^u(1-e^u) = 4\pi \sum \delta_{p_j}$ ，其线性化算子为标量，核空间维数为3；而本文系统(1.1)的线性化为2×2椭圆系统，核空间维数为6，且耦合项 $e^{u_1+u_2+u_{0,1}+u_{0,2}}(\phi_1+\phi_2)$ 引入强非对角耦合，导致标准分离变量法失效。本文通过构造耦合Liouville系统(2.6)的完整核并精细匹配内外区域，克服了这一障碍。

第三，相较于Lyapunov–Schmidt存在性工作（Huang [15]），本文在适用范围上互补：[15]在非退化临界点假设下构造解，但未讨论解的唯一性或线性稳定性；本文则在相同非退化性假设(A)下，证明所构造解的局部唯一性（即邻域内无其他解）与非退化性（即线性化无零模），为后续动力学研究（如涡旋运动方程）提供了必要前提。代价在于：本文要求 $D(\vec{q}) < 0$ 这一全局符号条件，而[15][16]未涉及；且结论仅适用于平均场型爆破（即局部质量趋近 $8\pi$ ），不覆盖其他爆破类型（如[1,2]中的多尺度爆破）。此外，本文结果局限于平环面 $\Omega$ ，未推广至一般曲面，因Green函数 $$G$$ 的曲率展开在一般情形下更为复杂。

4. Hasegawa-Mima 方程的弱解与耗散解

Weak and dissipative solutions for the Hasegawa-Mima equation

作者: Michele Gorini

类别: math.AP, math-ph

发布日期: 2026-07-01

arXiv ID: 2607.01119v1

摘要翻译: 我们考虑 Hasegawa-Mima 方程的“类 Euler”速度形式：

\partial_t(u-Δ^{-1}u)+(u\cdot\nabla)u-u^\perp\log n_0=0,

其中 $$n_0$$ 是粒子数密度 $n=n_0e^{\frac{e\varphi}{T}}$ 中出现的时间无关函数， $$u$$ 是漂移速度 $\nabla^\perp\varphi=-\nabla\varphi\times\hat z$ 。借鉴 Lions 关于 Euler 方程著作中提出的概念，我们证明：对任意 $$L^2$$ 散度自由初值 $w\in L^2(D)$ ，该方程在 $D=\mathbb T^2$ 及 $D\subset\mathbb R^2$ 为有界 $\mathcal{C}^1$ 区域的情形下均存在耗散解。

Comment: 无

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

1. 研究对象和背景

本文研究二维强磁化等离子体中描述漂移波湍流的基本模型——Hasegawa-Mima（HM）方程。该方程在无量纲形式下，势函数 $\varphi$ 满足

\partial_t(\Delta\varphi - \varphi) + (\nabla^\perp\varphi \cdot \nabla)(\Delta\varphi - \log n_0) = 0, \tag{1.1}

其中 $\nabla^\perp\varphi := (-\partial_2\varphi, \partial_1\varphi)$ 表示垂直于磁场方向的 $E \times B$ 漂移速度， $$n_0 = n_0(x,y) > 0$$ 是时间无关的平衡粒子数密度（满足 $n = n_0 e^{e\varphi/T_e}$ ），且 $\partial_t n_0 = 0$ 。物理上，(1.1) 源于电子绝热响应假设与强背景磁场 $$|B| = 1$$ 下的准中性近似，是二维等离子体物理中连接微观不稳定性与宏观输运的关键模型。

方程可等价重写为速度形式：令 $u := \nabla^\perp\varphi$ ，则 $$u$$ 是无散度向量场（ $\nabla \cdot u = 0$ ），且满足

\partial_t(u - \Delta^{-1}u) + (u \cdot \nabla)u - u^\perp \log n_0 = 0, \tag{1.2}

其中 $u^\perp = (-u_2, u_1)$ 。该形式在结构上高度类比于二维不可压 Euler 方程的涡度形式：若定义广义涡度 $q := \Delta\varphi - \varphi - \log n_0$ ，则 (1.1) 可化为活性标量输运方程

\partial_t q + (u \cdot \nabla) q = 0,

且 $u = \nabla^\perp \varphi = \nabla^\perp (-\Delta)^{-1}(q + \log n_0)$ ，即 $$u$$ 通过一个非局部算子由 $$q$$ 确定。当 $n_0 \equiv \text{const}$ 时， $q' = \Delta\varphi - \varphi$ ，此时 $u = \nabla^\perp (-\Delta)^{-1} q'$ ，故 $u \in H^s$ 当且仅当 $q' \in H^{s-1}$ ，体现出一阶正则性损失；而即使 $$n_0$$ 非常数，系统仍形式上守恒能量：

\frac{d}{dt} \int_D |\varphi|^2 + |\nabla\varphi|^2 \, dx = 0 \quad \Longleftrightarrow \quad \frac{d}{dt} \int_D \left|(-\Delta)^{-1/2} u\right|^2 + |u|^2 \, dx = 0.

此能量泛函 $\|u\|_{L^2}^2 + \|(-\Delta)^{-1/2} u\|_{L^2}^2$ 构成自然相空间 $\mathcal{H} := \{ u \in L^2(D; \mathbb{R}^2): \nabla \cdot u = 0 \}$ 上的 Hilbert 范数，其等价于 $\|u\|_{\dot{H}^{-1}}^2 + \|u\|_{L^2}^2$ 。

已有理论工作集中于弱解存在性，但存在显著局限。在周期域 $D = \mathbb{T}^2$ 上，文献 [6] 利用 Galërkin 逼近与紧性论证证明了对初始数据 $\varphi_0 \in H^2(\mathbb{T}^2)$ 的全局弱解存在性，但要求 $$n_0 = n_0(x)$$ 与 $$y$$ 无关，且未处理能量守恒或缺陷项；在有界 $$C^1$$ 域 $D \subset \mathbb{R}^2$ 上，文献 [3] 在 $\log n_0 \in L^p(D), p > 2$ 假设下建立了有限时间解，并给出先验估计 $\|\varphi - \Delta\varphi\|_{L^p} \lesssim \|\varphi - \Delta\varphi\|_{L^2} + \|\log n_0\|_{L^p}$ ，但唯一性仅在特殊情形成立；在全平面 $\mathbb{R}^2$ 上，文献 [4] 在相同 $\log n_0$ 条件下获得有限时间存在唯一性。所有这些结果均依赖于 $\varphi$ 的 $$H^2$$ 正则性，从而无法处理 $$L^2$$ 初始速度场这一更自然的低正则性初值问题。

根本障碍在于非线性项 $(u \cdot \nabla)u$ 的二次结构导致弱极限不封闭：若 $u_n \rightharpoonup u$ 在 $$L^2$$ 中弱收敛，则 $(u_n \cdot \nabla) u_n$ 不一定弱收敛到 $(u \cdot \nabla) u$ ，因缺乏紧性控制振荡。此现象与 Euler 方程完全一致，正是 Lions 在 [8] 中引入“耗散解”（dissipative solution）概念的原始动机：通过相对能量不等式定义一类对弱极限稳定的解，允许能量耗散缺陷，同时保持弱强唯一性（weak-strong uniqueness）。本文的核心空白即：尚未建立 HM 方程的耗散解理论，尤其缺乏针对 $u \in L^2$ 初值的全局存在性结果，亦未阐明耗散解与经典弱解、粘性正则化解之间的严格联系。这阻碍了对 HM 方程在湍流极限、小尺度振荡、边界层行为等物理关键场景下的数学建模能力。

2. 主要定理或主要结果

本文在周期域 $D = \mathbb{T}^2$ 和有界 $$C^1$$ 域 $D \subset \mathbb{R}^2$ 上，对 Hasegawa-Mima 方程的速度形式 (1.2) 证明了耗散解的全局存在性。主要结果如下：

定理 7.1（耗散解存在性）：设 $D = \mathbb{T}^2$ 或 $D \subset \mathbb{R}^2$ 为有界 $$C^1$$ 区域， $n_0 \in C^2(\overline{D})$ 且 $$n_0 > 0$$ ，则对任意发散自由初值 $w \in L^2(D; \mathbb{R}^2)$ ，存在一个耗散解 $u \in L^\infty_{\mathrm{loc}}([0,\infty); L^2(D))$ ，满足 $\nabla \cdot u(t) = 0$ 几乎处处，且对任意光滑测试场 $v \in C^\infty_c([0,T) \times D; \mathbb{R}^2)$ ， $\nabla \cdot v = 0$ ，定义 $g_v(t) := 2\big( \|\nabla^\perp \partial_1 v\|_{L^\infty} + \|\nabla^\perp \partial_2 v\|_{L^\infty} + \|\nabla \log n_0\|_{L^\infty} + \|\nabla v\|_{L^\infty} |\partial D| \big)$ ，并记 $K_{u,v}(t) := \int_D |u(t) - v(t)|^2 + \left|(-\Delta)^{-1/2}(u(t) - v(t))\right|^2 \, dx$ ，则 $$u$$ 满足相对能量不等式

K_{u,v}(t) \leq e^{\int_0^t g_v(s) \, ds} K_{u,v}(0) + \int_0^t e^{\int_s^t g_v(\tau) \, d\tau} \int_D (u - v) \cdot f_v(s,x) \, dx \, ds, \tag{4.15}

其中 $$f_v$$ 是 $$v$$ 满足的扰动方程 $-\partial_t(v - \Delta^{-1}v) = (v \cdot \nabla)v - v^\perp \log n_0 - f_v$ 的源项。特别地，当 $$v$$ 是 (1.2) 的经典解时， $f_v \equiv 0$ ，不等式退化为

K_{u,v}(t) \leq e^{\int_0^t g_v(s) \, ds} K_{u,v}(0),

即弱强唯一性成立。

引理 6.1（Leray 型粘性正则化解）：对粘性正则化方程 $\partial_t(u^\varepsilon - \Delta^{-1}u^\varepsilon) + (u^\varepsilon \cdot \nabla)u^\varepsilon - u^{\varepsilon,\perp} \log n_0 = \varepsilon \Delta u^\varepsilon$ ，存在全局强解 $u^\varepsilon \in L^\infty([0,T]; L^2(D)) \cap L^2([0,T]; H^1(D))$ ，且满足一致能量估计

\sup_{t \in [0,T]} \left( \|u^\varepsilon(t)\|_{L^2}^2 + \|(-\Delta)^{-1/2} u^\varepsilon(t)\|_{L^2}^2 \right) + \varepsilon \int_0^T \|\nabla u^\varepsilon(s)\|_{L^2}^2 \, ds \leq C(\|w\|_{L^2}, \|n_0\|_{C^2}),

其中常数 $$C$$ 独立于 $\varepsilon$ 。

定理 4.3（耗散解等价性）：若 $$u$$ 是 (1.2) 的耗散解，则 $\varphi := \operatorname{curl} \Delta^{-1} u$ 是 (1.1) 的耗散解；反之，若 $\varphi$ 是 (1.1) 的耗散解且 $\varphi|_{\partial D} = 0$ ，则 $u := \nabla^\perp \varphi$ 是 (1.2) 的耗散解。二者通过算子 $\nabla^\perp (-\Delta)^{-1}$ 一一对应。

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

本文的核心技术障碍在于：HM 方程 (1.2) 的非线性项 $(u \cdot \nabla)u$ 缺乏 $$L^2$$ 级别的乘子性质，导致标准弱解定义在 $$L^2$$ 初值下不具弱极限稳定性；同时，非局部算子 $\Delta^{-1}$ 与边界条件耦合，使得在有界域上构造耗散解需克服谱理论与迹算子的双重困难。作者的真正创新点在于将 Lions 的耗散解框架从 Euler 方程系统性移植至 HM 方程，并完成三项关键适配：

1. 相对能量泛函的精确构造与消去技巧：针对速度形式 (1.2)，作者定义相对能量为

K_{u,v}(t) := \int_D |u - v|^2 + \left|(-\Delta)^{-1/2}(u - v)\right|^2 \, dx,

这是 $\|u - v\|_{\mathcal{H}}^2$ 的自然选择，对应于势函数差 $\|\varphi - \psi\|_{H^1}^2$ 。关键突破在于对交叉项 $\int (u - v) \cdot [(u \cdot \nabla)u - (v \cdot \nabla)v]$ 的精细处理：利用恒等式 $\nabla^\perp f \cdot \nabla f = 0$ 和 $\nabla \cdot \nabla^\perp f = 0$ （见 (4.5)-(4.6)），通过两次分部积分彻底消去所有含三阶导数的危险项（如原文第 7 页中 III 项恒为零，II 项经两次 IBP 后边界项与体积分均被控制），最终将非线性误差归结为可控的 $L^\infty$ 范数乘积，即 $\|\nabla^\perp \partial_i v\|_{L^\infty}$ 与 $\|\nabla(u - v)\|_{L^2}^2$ 的耦合。

2. 粘性正则化与紧性框架的严格实现：证明存在性的核心是构造 $\varepsilon$ -粘性逼近解 $u^\varepsilon$ 并取 $\varepsilon \to 0$ 极限。作者并非简单引用 Navier-Stokes 的 Leray 理论，而是为 HM 方程专门建立引理 6.1：通过将方程改写为 $\partial_t u^\varepsilon + \mathcal{L} u^\varepsilon = \mathcal{N}(u^\varepsilon) + \varepsilon \Delta u^\varepsilon$ ，其中 $\mathcal{L} = -\partial_t \Delta^{-1} + \nabla \cdot$ 是线性双曲-椭圆混合算子，再结合 $$L^2$$ 能量估计与 Aubin-Lions 引理，证明 $u^\varepsilon$ 在 $L^\infty_t L^2_x$ 中弱-* 收敛、在 $$L^2_t L^2_x$$ 中强收敛。此过程的关键是验证 $\varepsilon \|\nabla u^\varepsilon\|_{L^2_{t,x}}^2$ 一致有界，从而保证耗散项在极限中消失，仅留下耗散解所需的相对能量不等式。

3. 有界域上分数阶拉普拉斯算子的谱适配：为使 $(-\Delta)^{-1/2}$ 在 $D \subset \mathbb{R}^2$ 上良定义，作者明确要求 $\Delta$ 为 Dirichlet 拉普拉斯，并指出可通过特征函数展开定义 $(-\Delta)^\alpha$ （傅里叶符号 $|\xi|^{2\alpha}$ 在 $\mathbb{T}^2$ 上自然成立，而在有界域则需利用 $\Delta$ 的离散谱 $\{\lambda_k\}$ 定义 $(-\Delta)^\alpha f := \sum_k \lambda_k^\alpha \langle f, e_k \rangle e_k$ ）。此适配确保了边界项 $\int_{\partial D} \cdots \, d\sigma$ 在 $u|_{\partial D} = v|_{\partial D} = 0$ 下严格消失（见第 10 页推导），是将周期域结果推广至有界域的基石。

4. 与之前工作的比较

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

1. 解的概念层级与初值正则性：此前所有存在性结果（[3,4,6]）均要求 $\varphi_0 \in H^2$ 或等价地 $u_0 \in H^1$ ，而本文将初值空间降至最优的 $$L^2$$ ，且无需任何额外光滑性或小初值假设。例如，[6] 的 Galërkin 解需 $\varphi_0 \in H^2$ 以保证非线性项在逼近中强收敛；[3] 的解依赖 $\log n_0 \in L^p$ 且仅在有限时间内存在。本文的耗散解对任意 $w \in L^2$ 全局存在，其代价是放弃经典弱解的逐点满足性，转而接受相对能量不等式中的耗散缺陷，这恰是处理湍流多尺度振荡的数学必然。

2. 边界处理与适用范围：[3] 在有界域上的结果受限于 $\log n_0 \in L^p$ 且未建立弱强唯一性；本文则通过明确定义 Dirichlet 分数阶拉普拉斯 $(-\Delta)^\alpha$ ，首次为 HM 方程在一般 $$C^1$$ 有界域上建立了完整的耗散解理论（见第 5 节备注），其假设 $n_0 \in C^2(\overline{D})$ 比 [3] 的 $$L^p$$ 条件更强，但换来了全局存在性与弱强唯一性的统一框架。此推广的代价是技术复杂度剧增：必须处理特征函数的边界迹、分数阶算子的谱表示，以及 $\nabla^\perp \partial_i v$ 在边界附近的 $L^\infty$ 控制（见 $$g_v(t)$$ 中的 $|\partial D|$ 项）。

3. 与 Euler 方程理论的同构性与差异性：本文直接继承 Lions [8] 的耗散解定义，但关键差异在于相对能量泛函：Euler 方程使用 $\|u - v\|_{L^2}^2$ ，而 HM 方程必须使用 $\|u - v\|_{L^2}^2 + \|(-\Delta)^{-1/2}(u - v)\|_{L^2}^2$ ，这反映了 $$u$$ 与 $\varphi$ 之间固有的非局部关系。因此，本文的 $$g_v(t)$$ 中出现了 $\|\nabla^\perp \partial_i v\|_{L^\infty}$ 项，而非 Euler 中的 $\|\nabla v\|_{L^\infty}$ ，体现了 HM 方程中 $\Delta^{-1}$ 算子对高阶导数的放大效应。从提供内容中无法确认本文是否改进了 [4] 在 $\mathbb{R}^2$ 上的有限时间结果，作者亦未显式展开与 [4] 的比较。

5. 具非线性扩散的动理学Fokker-Planck方程

Kinetic Fokker-Planck Equations with Nonlinear Diffusion

作者: Zimo Hao, Zhengyan Wu, Xicheng Zhang

类别: math.AP, math.PR

发布日期: 2026-07-01

arXiv ID: 2607.00458v1

摘要翻译: 我们在 $\mathbb R^{2d}$ 上研究非线性动理学Fokker–Planck方程

\partial_t f=\Delta_v\Psi(f)-v\cdot\nabla_x f, \qquad f|_{t=0}=f_0,

的存在性、正则性与唯一性。在模型情形 $\Psi(r)=r^s$ 下，该方程将非线性快扩散/多孔介质型扩散与动理学输运耦合起来。其显著特征在于扩散仅作用于速度变量 $$v$$ ，因此无法通过标准椭圆估计获得空间变量 $$x$$ 的紧性，而必须借助 hypoelliptic 结构予以恢复。

Comment: 无

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

1. 研究对象和背景

本文研究定义在全空间 $\mathbb{R}^{2d}$ 上的非线性动理学福克-普朗克方程

\partial_t f = \Delta_v \Psi(f) - v \cdot \nabla_x f, \quad f|_{t=0} = f_0, \tag{1.1}

其中未知函数 $f = f(t,x,v) : [0,T] \times \mathbb{R}^{2d} \to [0,\infty)$ 表示相空间密度， $x \in \mathbb{R}^d$ 为位置变量， $v \in \mathbb{R}^d$ 为速度变量；非线性扩散函数 $\Psi : [0,\infty) \to \mathbb{R}$ 满足结构条件（见下文），初始数据 $f_0 \in L^1(\mathbb{R}^{2d}) \cap L^2(\mathbb{R}^{2d})$ 非负。该方程的核心特征在于：扩散项 $\Delta_v \Psi(f)$ 仅作用于速度变量 $$v$$ ，而输运项 $-v \cdot \nabla_x f$ 耦合了 $$x$$ 与 $$v$$ ，但不提供任何 $$x$$ 方向的耗散。因此，标准椭圆或抛物型正则性理论无法直接给出 $$x$$ 变量的紧性——这是与经典多孔介质方程 $\partial_t u = \Delta (u^m)$ 或各向同性扩散-输运耦合模型的根本区别。从物理角度看，(1.1) 描述一类分布依赖型随机动力系统 $$dX = V dt$$ , $dV = \sqrt{2a(f)} dB$ 的宏观平均场极限，其中 $a(\zeta) := \Psi(\zeta)/\zeta$ ，故其非线性源于粒子间通过局部密度 $$f$$ 实现的自洽相互作用。数学上，它属于非线性动理学方程范畴，兼具三重困难：(i) 非线性且可能退化/奇异的扩散系数 $\Psi'(f)$ ；(ii) 动理学退化结构（仅 $$v$$ -方向扩散）；(iii) 缺乏全局势能或熵结构（无摩擦项 $\nabla_v \cdot (vf)$ ）。已有工作如 [BCD26] 仅处理幂律情形 $\Psi(\zeta) = \zeta^s$ 在 $s \in (1-1/d,1) \cup (1,1+1/d)$ 的基本解存在性；[BDM26] 在含摩擦项 $\nabla_v \cdot (vf)$ 的变体中建立弱解存在性，但其熵方法严重依赖局部平衡态结构，无法移植至无摩擦的 (1.1)。关键空白在于：对最自然的快扩散区间 $s \in (0,1)$ ，既缺乏适用于一般 $\Psi$ 的存在性框架，更无 $$L^1$$ -唯一性理论；尤其当 $$s < 1-1/d$$ （ $d \ge 2$ ）时，质量守恒与唯一性均失效，而该临界指数 $$1-1/d$$ 恰由相空间尺度不变性决定：若设自相似解 $f(t,x,v) = t^{-\alpha} F(x/t^{\beta}, v/t^{\gamma})$ ，由 $\partial_t f$ 与 $\Delta_v f^s$ 平衡得 $\alpha + 1 = s\alpha + 2\gamma$ ，由输运 $v \cdot \nabla_x f$ 得 $\beta = \gamma + 1$ ，由质量守恒 $\int f = \text{const}$ 得 $\alpha = d(\beta + \gamma) = d(2\gamma + 1)$ ，联立解得 $\alpha = 2d / (1 - d(1-s))$ ，故 $\alpha > 0$ 当且仅当 $$s > 1-1/d$$ ，端点 $$s = 1-1/d$$ 对应临界发散行为。本文旨在填补这一空白：构建适用于广义 $\Psi$ （含全部 $s \in (0,1)$ ）的存在性理论，并将 $$L^1$$ -唯一性推进至该质量临界阈值。

2. 主要定理或主要结果

在假设 1.1 下（即 $\Psi \in C([0,\infty))$ ， $\Psi', \Psi'' \in C(0,\infty)$ ， $\Psi(\zeta) > 0$ 、 $\Psi'(\zeta) > 0$ 对所有 $\zeta > 0$ 成立，且存在常数 $$0 < r_2 < r_1 < 2$$ 、 $\lambda \in (1,2]$ 、 $l \in (0,1/3)$ 、 $$C > 0$$ 使得 $\Psi(\zeta) \le C(\zeta^{r_2} + \zeta^{r_1})$ 且 $|\Psi''(\zeta)| / [\Psi'(\zeta)]^{1+l} \le C \zeta^{\lambda-2}$ ），本文证明：

定理 1.3（存在性）：存在非负弱解 $$f$$ （定义 2.3），且对任意 $p \in (1, 2d/(2d-l))$ 和 $\beta := 2l - 4(d - d/p)$ ，有统一先验估计

\sup_{t \in [0,T]} \|f(t)\|_{L^2}^2 + \int_0^T \int_{\mathbb{R}^{2d}} |\nabla_v H(f(t,z))|^2 dz dt + \int_0^T \|f(t)\|_{B^{\beta}_{p;a}} dt \le C,

其中 $H(\zeta) := \int_0^{\zeta} \sqrt{\Psi'(\zeta')} d\zeta'$ ， $$a = (3,1)$$ 为动理学各向异性标度， $C = C(d,p,\beta,l,T,\Psi,\|f_0\|_{L^1 \cap L^2}, c_0)$ ，且 $c_0 := \int_{\mathbb{R}^{2d}} \int_0^{f_0} |\Psi'(\zeta)|^{-l} d\zeta dz < \infty$ 。

定理 1.4（唯一性）：若额外满足 $r_2 \ge 1-1/d$ （当 $d \ge 2$ ），则定理 1.3 构造的弱解是重正规化动理学解（定义 2.2），并保持质量：

\|f(t)\|_{L^1(\mathbb{R}^{2d})} = \|f_0\|_{L^1(\mathbb{R}^{2d})}, \quad \forall t \in [0,T]. \tag{1.3}

此解在满足 (1.3) 的重正规化动理学解类中唯一。特别地，在幂律情形 $\Psi(\zeta) = \zeta^s$ 下，该条件精确等价于 $s \ge 1-1/d$ （ $d \ge 2$ ），且 $$d=1$$ 时整个区间 $s \in (0,1)$ 均被覆盖。

推论 1.5（随机表示）：在假设 1.1 且 $$r_2 > 1-1/d$$ 下，若 $$f_0$$ 为概率密度，则存在分布密度依赖型 SDE 系统 $$dX = V dt$$ , $dV = \sqrt{2\Psi(f)/f} dB$ 的非线性鞅解（定义 2.7）。

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

本文核心障碍在于：能量估计仅控制 $\nabla_v H(f)$ ，无法获得 $$x$$ 变量的任何导数信息，故标准 Aubin–Lions 紧性引理失效。作者的真正创新点在于构造并运用一个参数依赖的动理学半群光滑估计，该估计定量刻画了扩散系数 $\Psi'(\zeta)$ 及其导数 $\Psi''(\zeta)$ 对半群行为的影响。具体而言：

动理学公式化与参数化线性化：引入动理学函数 $\chi(t,x,v,\zeta) := \mathbf{1}_{\{f(t,x,v)>\zeta>0\}}$ ，形式上满足
$\partial_t \chi = \Psi'(\zeta) \Delta_v \chi - v \cdot \nabla_x \chi + \partial_\zeta q, \tag{1.4}$
其中动能测度 $$q$$ 满足 $q \ge \delta_{f=\zeta} \Psi'(\zeta) |\nabla_v f|^2$ （重正规化定义的关键）。此式将非线性问题转化为一族以 $\zeta$ 为参数的线性动理学方程，但 $\zeta$ -导数 $\partial_\zeta q$ 的出现要求估计必须对 $\zeta$ 微分稳定。
参数依赖半群估计（引理 3.1）：对每个 $\zeta > 0$ ，定义由算子 $\Psi'(\zeta) \Delta_v - v \cdot \nabla_x$ 生成的半群 $P_t(\zeta)$ ，其核为 $(\Gamma_t p_t(\zeta)) * (\Gamma_t h)(x,v)$ ，其中 $p_t(\zeta)$ 是扩散系数为 $\Psi'(\zeta)$ 的 Kolmogorov 过程高斯密度， $\Gamma_t h(x,v) := h(x-tv,v)$ 。作者证明了关键估计：
$\left\| \partial_\zeta \left( \Gamma_{t-s} p_{t-s}(\zeta) \right) * (\Gamma_{t-s} q) \right\|_{B^{\beta}_{p;a}} \lesssim (t-s)^{-\alpha} \|q\|_{L^1} \cdot \left( |\Psi''(\zeta)| [\Psi'(\zeta)]^{-1-l} + [\Psi'(\zeta)]^{-l} \right),$
其中 $\alpha, \beta$ 由各向异性标度 $$a=(3,1)$$ 决定，且常数显式依赖于 $\Psi'(\zeta), \Psi''(\zeta)$ 。该估计的 $\zeta$ -稳定性（由假设中 $$l < 1/3$$ 保证 $\partial_\zeta$ 奇异性可积）与 $\Psi'$ -依赖性是突破所在。
各向异性 Besov 紧性：将上述估计与动能测度的能量控制 $\int_0^T \int_{\mathbb{R}^{2d}} \int_0^\infty q \, d\zeta dz dt \le C$ 结合，得到 $$f$$ 的统一各向异性 Besov 估计 $\|f\|_{B^{\beta}_{p;a}}$ ，该空间嵌入紧致于 $L^1_{\text{loc}}$ ，从而恢复 $$x$$ 与 $$v$$ 的联合紧性。
重正规化唯一性论证：在 $$L^1$$ 框架下采用双变量法，其 coercivity 依赖于动能测度的支配性质 $q \ge \delta_{f=\zeta} \Psi'(\zeta) |\nabla_v f|^2$ ；而移除低动能截断的关键则来自质量守恒导出的小 $\zeta$ 衰减： $\lim_{\beta \downarrow 0} \beta^{-1} q(\mathbb{R}^{2d} \times [\beta/2, \beta] \times [0,T]) = 0$ （命题 2.5），该性质仅在 $r_2 \ge 1-1/d$ 下成立。

4. 与之前工作的比较

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

存在性范围：[BCD26] 仅对 $\Psi(\zeta) = \zeta^s$ 且 $s \in (1-1/d,1) \cup (1,1+1/d)$ 构造基本解，未涉及一般 $\Psi$ 或弱解理论；[BDM26] 处理含摩擦项 $\nabla_v \cdot (vf)$ 的方程，其熵方法要求 $$s > 1-1/d$$ 以保证熵不等式闭合，但该框架完全失效于无摩擦的 (1.1)。本文将存在性推广至所有满足假设 1.1 的 $\Psi$ ，特别覆盖整个快扩散区间 $s \in (0,1)$ ，且无需任何 $$s > 1-1/d$$ 的先验限制。

唯一性阈值与代价：[BDM26] 的唯一性依赖于熵解框架，其适用范围为 $$s > 1-1/d$$ ，但未触及端点 $$s = 1-1/d$$ ；本文首次在 $$L^1$$ 类中实现端点唯一性，即 $$s = 1-1/d$$ （ $d \ge 2$ ）。代价是必须采用重正规化动理学解概念（而非经典弱解），且唯一性类限定为质量守恒解。该限定是本质的：当 $$s < 1-1/d$$ 时，截断论证显示质量损失不可避免，故 $$L^1$$ -唯一性在此以下区间不成立（文中明确指出）。

分析工具革新：经典非线性扩散理论（如 [Váz07]）依赖于全空间抛物正则性，对 (1.1) 完全无效；线性动理学理论（如 [Bou02, HWZ20]）仅处理固定系数 $\Delta_v - v \cdot \nabla_x$ ，无法处理 $\Psi'(\zeta)$ 依赖于未知量的参数化情形；动能公式化在标量守恒律中用于处理激波，但此处其核心新作用是将非线性紧性问题转化为参数化线性半群的稳定性问题。本文的参数依赖半群估计是全新工具，其 $\zeta$ -微分稳定性设计（通过 $$l < 1/3$$ 控制奇异性）和 $\Psi'$ -显式依赖性是此前文献未见的。从提供内容中无法确认本文结果是否可推广至带外力场或边界条件的情形。

6. 具有常涡量的三维毛细液滴的刚性结果

A rigidity result for the 3D capillary liquid drop with constant vorticity

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

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00450v1

摘要翻译: 我们考虑流体动力学 Euler 方程组所描述的三维毛细液滴的自由边界问题，其中毛细系数为 $$σ_0$$ ，且液滴形状接近球形，假设其涡量恒为 $$(0, 0, α_0)$$ 。首先，我们研究常涡量条件与系统时间演化之间的相容性，证明：当 $α_0 \neq 0$ 时，任何定义在凸区域上的光滑解必须满足关于流体区域形状的强几何约束；并且，常涡量条件（不同于无旋情形 $$α_0 = 0$$ ）并不构成系统时间演化的不变集。随后，我们聚焦于该问题的稳态解，并证明一个新刚性结果：不预先假设任何对称性条件，我们证明：若比值 $$α_0^2/σ_0$$ 不过大，则任意接近球形的解必具有柱对称性，从而即为文献中已知的唯一轴对称解；此时流体区域接近但不等于球体，更精确地说，它是一个扁球面（oblate spheroid），极区压扁、赤道隆起；且每个流体质点沿水平圆周轨迹以恒定角速度运动。据我们所知，这是首个在未预设柱对称性假设下得到的关于具有常涡量的毛细液滴的结果。

Comment: 无

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

1. 研究对象和背景

本文研究三维空间中具有表面张力 $\sigma_0 > 0$ 的自由边界不可压理想流体液滴的稳态构型，其运动由欧拉方程控制，并附加常涡度约束：

\begin{cases} \langle u, \nabla \rangle u = -\nabla p & \text{in } \Omega, \\ \operatorname{div} u = 0 & \text{in } \Omega, \\ \operatorname{curl} u = \alpha_0 e_3 & \text{in } \Omega, \\ p = \sigma_0 H_{\partial\Omega} & \text{on } \partial\Omega, \\ \langle u, \nu_{\partial\Omega} \rangle = 0 & \text{on } \partial\Omega. \end{cases} \tag{1.5}

其中 $\Omega \subset \mathbb{R}^3$ 是光滑有界开集（液滴占据区域）， $u : \Omega \to \mathbb{R}^3$ 是速度场， $$p$$ 是压力， $$e_3 = (0,0,1)$$ ， $H_{\partial\Omega}$ 是边界 $\partial\Omega$ 的平均曲率（按外法向定义，凸域为正）， $\nu_{\partial\Omega}$ 是单位外法向量。未知量为 $\Omega$ 和 $$u$$ ；压力 $$p$$ 由泊松问题唯一确定： $\Delta p = -\operatorname{tr}(\nabla u \cdot \nabla u)$ 在 $\Omega$ 内， $p|_{\partial\Omega} = \sigma_0 H_{\partial\Omega}$ 。该模型描述无重力、无粘性、无表面活性剂、仅受毛细力驱动的旋转液滴平衡态——即 Plateau 实验中观察到的经典“扁球形旋转液滴”。

关键数学挑战在于：当 $\alpha_0 = 0$ （无旋情形）时，系统退化为经典的 Serrin 型自由边界问题，其唯一性由 Alexandrov 定理保证（常平均曲率 $\Rightarrow$ 球形）；但当 $\alpha_0 \neq 0$ 时，涡度非零彻底破坏了变分结构与对称性先验：(i) 常涡度条件 $\operatorname{curl} u = \alpha_0 e_3$ 不是时间演化不变量（Theorem 1.5 构造反例表明 Cauchy 问题不适定）；(ii) 即使对稳态解， $$u$$ 不再是梯度场，而必须包含刚体旋转部分 $\frac{\alpha_0}{2}(-x_2,x_1,0)$ 与调和扰动 $$v$$ 的叠加（见 (4.5)）；(iii) 边界条件 $p = \sigma_0 H_{\partial\Omega}$ 与动量方程耦合后，不再导出纯几何的常曲率方程，而是生成含速度梯度项的非线性边界关系（见 (1.9)）。因此，传统方法（如移动平面法直接作用于 $\Omega$ 、或基于能量泛函 $L(\Omega) = \sigma_0 P(\Omega) - \gamma_0 \int_\Omega |x'|^2 dx$ 的临界点分析 [24,39]）均失效：前者因 $$u$$ 的动力学耦合而无法分离几何；后者仅考虑区域 $\Omega$ 而忽略 $$u$$ 的内在约束，其 Euler-Lagrange 方程 $\sigma_0 H_{\partial\Omega} - \gamma_0 |x'|^2 = \lambda$ 缺失关键项（对比 (1.9) 中的调和函数 $$g$$ ）。本文旨在填补这一根本性空白：在不预设任何对称性（尤其是不假设柱对称）的前提下，证明只要 $\alpha_0^2 / \sigma_0$ 足够小，所有近球形稳态解必具柱对称性，且唯一对应已知的扁球面解。这要求建立全新的几何-分析桥梁——将三维流体问题降维为二维区域 $D \subset \mathbb{R}^2$ （赤道截面）上的超定椭圆问题，其核心对象是 $$D$$ 的扭转函数 $$v_D$$ （满足 $\Delta v_D = -1$ in $$D$$ , $$v_D = 0$$ on $\partial D$ ），并揭示 $$v_D$$ 的法向导数 $|\nabla v_D|$ 与边界曲率 $H_{\partial D}$ 之间由物理参数 $\alpha_0, \sigma_0$ 刻画的精确关系（见 (4.3)）。该框架首次将扭转函数的凹性理论（Makar-Limanov 定理）、梯度估计（Theorem 2.4）与流体动力学约束深度融合，从而绕过传统变分方法的局限。

2. 主要定理或主要结果

本文核心结论为 Theorem 1.7，其严格陈述如下：设 $\Omega \subset \mathbb{R}^3$ 为严格凸开集且 $\partial\Omega \in C^2$ ， $u \in C^2(\Omega; \mathbb{R}^3) \cap C(\overline{\Omega}; \mathbb{R}^3)$ ，且 $(\Omega, u)$ 满足稳态系统 (1.5)。则：

(i) 存在 $\lambda_0 \in \mathbb{R}$ 、严格凸开集 $D \subset \mathbb{R}^2$ 及函数 $f \in C^2(D) \cap C(\overline{D})$ 满足 $$f > 0$$ in $$D$$ 、 $$f = 0$$ on $\partial D$ ，使得

\Omega = \{(x', x_3) \in \mathbb{R}^3 : x' \in D,\ |x_3 - \lambda_0| < f(x')\}.

且 $$u$$ 与 $$x_3$$ 无关，第三分量 $u_3 \equiv 0$ 。

(ii) 若参数满足尺度条件

\frac{\alpha_0^2}{\sigma_0} |D|^{3/2} \leq \sqrt{2}\pi^3, \tag{1.6}

则 $$D$$ 必为圆盘， $\partial\Omega$ 为旋转曲面，且速度场必为

u = \frac{\alpha_0}{2} (-x_2, x_1, 0).

进一步，由 (i) 可推出 $\Omega$ 关于平面 $\{x_3 = \lambda_0\}$ 对称，且 $$u$$ 具柱对称性；由 (ii) 可得 $\Omega$ 为扁球面（oblate spheroid），其显式轮廓由径向函数 $$f(|x'|)$$ 通过积分确定（Remark 4.4），且每个流体质点沿水平圆周以恒定角速度 $\alpha_0/2$ 运动。该结果的临界阈值 (1.6) 本质源于扭转函数梯度估计 $\|\nabla v_D\|_{L^\infty(D)} \leq c |D|^{1/2}$ （Theorem 2.4）与曲率积分恒等式 $\int_{\partial D} H_{\partial D} \, d\sigma = 2\pi$ （Gauss-Bonnet）的耦合。值得注意的是，Theorem 1.7 的证明不依赖 $\Omega$ 的高阶正则性（ $$C^5$$ ），仅需 $$C^2$$ 凸性（Remark 1.8, 1.9），且结论中 $$D$$ 的圆盘性是刚性而非渐近性：一旦 (1.6) 成立， $$D$$ 严格为圆盘，无误差项或小参数展开。

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

证明的核心障碍在于：三维流体方程 (1.5) 与常涡度约束的强耦合，使得无法直接对 $\Omega$ 应用 Alexandrov 反射原理或 Serrin 对称性论证。作者的根本创新在于构建一个精巧的降维框架，将问题完全转化为二维赤道截面 $$D$$ 上的超定椭圆问题，其枢纽是扭转函数 $$v_D$$ 。具体技术路径分三步：

第一步：几何对称性提取与降维参数化。利用常涡度条件 $\operatorname{curl} u = \alpha_0 e_3$ 导出 $\partial_{x_3} u = 0$ （由 (3.2)），结合 $\langle u, \nu_{\partial\Omega} \rangle = 0$ 与凸性，通过 Alexandrov 反射原理（Theorem 2.1）证明 $\Omega$ 必关于某水平面 $\{x_3 = \lambda_0\}$ 对称（Theorem 1.1）。进而将 $\partial\Omega$ 参数化为 $\pm f(x')$ 在 $D = \Omega \cap \{x_3 = \lambda_0\}$ 上的图（Lemma 4.3），此参数化规避了球面坐标下梯度爆炸的困难，且天然适配柱对称性目标。

第二步：构造扭转函数超定问题。将速度场分解为刚体旋转部分与调和扰动： $u = v + \frac{\alpha_0}{2} J_3 x$ （(4.5)），其中 $$J_3 x = (-x_2, x_1, 0)$$ 。由 $\operatorname{div} v = 0$ 且 $\operatorname{curl} v = 0$ ，在单连通凸域 $$D$$ 上存在势函数 $\phi$ 使 $\tilde{v} = \nabla \phi$ （(4.7)）。关键洞察是：将动量方程 $\langle u, \nabla \rangle u = -\nabla p$ 在边界投影，并利用 $p = \sigma_0 H_{\partial\Omega}$ 及曲率公式 (2.2)，最终导出 $$D$$ 的边界 $\partial D$ 上的超定关系（Theorem 4.2）：

\frac{\alpha_0^2}{2} |\nabla v_D|^2 - \frac{\sigma_0 \beta}{|\nabla v_D|} + \sigma_0 H_{\partial D} = c \quad \text{on } \partial D, \tag{4.3}

其中 $\beta, c$ 为仅依赖 $$D$$ 的常数（由 Lemma 4.9, 4.10 显式给出）。此方程将三维流体动力学完全编码于二维扭转函数 $$v_D$$ 的法向梯度 $|\nabla v_D|$ 与边界曲率 $H_{\partial D}$ 的代数关系中，是全文最核心的技术突破。

第三步：扭转函数的尖锐分析与刚性判定。为从 (4.3) 推出 $$D$$ 为圆盘，需控制常数 $\beta$ 。作者利用凸域扭转函数的 $$1/2$$ -凹性（Theorem 2.3）及梯度估计（Theorem 2.4）： $\|\nabla v_D\|_{L^\infty} \leq c |D|^{1/2}$ ，结合 $\int_{\partial D} |\nabla v_D| \, d\sigma = |D|$ （(2.8)）与 $\int_{\partial D} H_{\partial D} \, d\sigma = 2\pi$ ，推导出 $\beta$ 的显式上界（(4.40)）。将此界代入 (4.3)，发现当 $\alpha_0^2 / \sigma_0$ 满足 (1.6) 时，方程强制 $|\nabla v_D|$ 在 $\partial D$ 上为常数。再应用 Serrin 刚性定理（Theorem 2.2）的推广形式（ $|\nabla v_D|$ 为常数 $\Rightarrow$ $$D$$ 为圆盘），完成证明。整个框架的独创性在于：未引入任何小参数摄动或隐函数定理，而是通过扭转函数的内在几何不等式（而非泛函分析工具）实现刚性判定。

4. 与之前工作的比较

本文与既有工作存在本质差异，主要体现在四个维度：

1. 对称性假设的彻底移除。此前所有关于常涡度毛细液滴的研究（如 [24, 39, 7, 32]）均预设柱对称性：即直接设 $\Omega$ 为旋转曲面、 $$u$$ 为 $\frac{\alpha_0}{2}(-x_2,x_1,0)$ 。本文首次在无需任何对称性先验下，证明该对称性是解的必然结果（Theorem 1.7 (ii)），这是概念性突破。例如，[24] 中的“稳态旋转解”仅定义为泛函 $L(\Omega)$ 的临界点，其 Euler-Lagrange 方程 (1.8) 隐含假设了轴对称，而本文从第一性原理的欧拉方程出发，导出相同结论。

2. 适定性与刚性范围的精确刻画。文献 [37, 29] 对三维常涡度流体建立了局部刚性（如局部轴对称），但依赖强正则性（ $$C^5$$ ）或小初值。本文的刚性条件 (1.6) 给出了明确的物理参数阈值： $\alpha_0^2 / \sigma_0$ 控制涡度强度与毛细力的竞争。当此比值过大时，(4.3) 不足以强制 $|\nabla v_D|$ 为常数，故结论失效——这解释了为何高速旋转液滴会失稳（如 Jacobi 椭球形），而本文给出了从球形到扁球形转变的定量判据。相比之下，[37] 仅给出存在性，无此类尺度关系。

3. 分析工具的根本革新。传统方法（[2,3] 的球面参数化、[35] 的 Serrin 方法）无法处理 $$u$$ 的动力学耦合。本文首创扭转函数超定框架：将三维问题降维为二维 $$D$$ 上的 (4.3)，其创新在于识别 $$v_D$$ 作为连接几何 ( $H_{\partial D}$ ) 与动力学 ( $\alpha_0, \sigma_0$ ) 的唯一桥梁。这超越了 [34] 中 Reilly 公式的应用，也不同于 [25,17] 对扭转函数凹性的纯几何研究，而是将其与流体方程深度耦合。

4. 局限性与代价。本文结论要求 $\Omega$ 为凸域（虽可弱化为凸性，Remark 1.8），且仅适用于稳态解；对非稳态解，Theorem 1.5 已证明常涡度条件通常不守恒，故无法推广。此外，阈值 (1.6) 中的常数 $\sqrt{2}\pi^3$ 并非最优（源于 Theorem 2.4 的非紧梯度估计），改进此常数需更精细的扭转函数分析。而 [9,8] 研究全局存在性时采用完全不同的 Lagrangian 方法，与本文的 Eulerian 刚性框架互补而非竞争。综上，本文不是对前人工作的渐进改进，而是开辟了用凸几何分析解决流体刚性问题的新范式。

7. 一个极小极大双线性传输问题与Nash-Monge-Kantorovich映射

A minimax Bilinear Transport Problem and Nash-Monge-Kantorovich Maps

作者: Rene Cabrera, Edward Huynh

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00344v1

摘要翻译: 我们研究一个源于具有二次动能与相互作用代价的二人零和博弈的极小极大双线性传输问题。从动态路径空间表述出发，我们建立了极小极大与极大极小规划的存在性，并证明了一个极小极大定理。我们证明该均衡通过传输规划上的端点代价诱导出一个有限维稳态问题：该端点代价在临界相互作用强度以下良定义，并在耦合上给出一个Nash均衡。在二次相互作用情形下，我们导出了一个显式的端点代价及一个对偶表述。所得的Nash-Monge-Kantorovich (NMK) 规划容许Monge解，从而恢复最优传输中的经典结构；当其存在时，最优映射由凸函数或凹函数的梯度给出。我们的分析凸显了非标准代价下的对偶性与循环（反）单调性，并将均衡映射与耦合的非线性PDEs相联系，从而架起了最优传输、零和博弈与Monge-Ampère型方程之间的桥梁。

Comment: 53 pages, 7 figures, comments welcome

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

1. 研究对象和背景

本文研究一个源于二人零和博弈的极小极大双线性传输问题，其核心对象是定义在路径空间上的双变量成本泛函

c(\gamma, \xi) = \int_0^1 \left[ \frac{1}{2}|\dot{\gamma}(t)|^2 - \frac{1}{2}|\dot{\xi}(t)|^2 + \alpha |\gamma(t) - \xi(t)|^2 \right] dt, \quad \gamma \in \Omega_1,\ \xi \in \Omega_2,

其中 $\Omega_i = \{ \gamma \in W^{1,2}([0,1]; \mathbb{R}^d) : \int_0^1 |\dot{\gamma}|^2 \leq K_i \}$ 是带能量约束的 Sobolev 路径空间， $\mu_1, \nu_1 \in \mathcal{P}(\mathbb{R}^d)$ 、 $\mu_2, \nu_2 \in \mathcal{P}(\mathbb{R}^d)$ 是具有紧支集的绝对连续概率测度，且 $\pi_1 \in \Pi_{\text{path}}(\mu_1, \nu_1)$ 、 $\pi_2 \in \Pi_{\text{path}}(\mu_2, \nu_2)$ 分别为满足端点约束 $(e_0, e_1)_\# \pi_i = (\mu_i, \nu_i)$ 的路径空间上概率测度（即路径型传输计划）。目标泛函为双线性作用

\sigma_0(\pi_1, \pi_2) := \iint_{\Omega_1 \times \Omega_2} c(\gamma, \xi) \, d\pi_1(\gamma) \, d\pi_2(\xi),

并考虑极小极大优化问题

\inf_{\pi_1 \in \Pi_{\text{path}}(\mu_1,\nu_1)} \sup_{\pi_2 \in \Pi_{\text{path}}(\mu_2,\nu_2)} \sigma_0(\pi_1, \pi_2). \tag{1}

此问题本质上是将经典最优传输理论推广至对抗性设置：不再是单一方最小化总传输代价，而是两方以相反符号（动能项 $\frac{1}{2}|\dot{\gamma}|^2$ 与 $-\frac{1}{2}|\dot{\xi}|^2$ ）竞争，并通过二次相互作用项 $\alpha |\gamma - \xi|^2$ 耦合其轨迹。该模型自然对应于微分博弈中的追逃问题（如图1所示），其中玩家策略为路径分布而非点态控制律。关键未知量包括：(i) 极小极大解 $(\pi_1^*, \pi_2^*)$ 的存在性与唯一性；(ii) 极小极大值与极大极小值是否相等（即是否成立 $\min\max = \max\min$ ）；(iii) 当解存在时，能否退化为确定性映射（即是否存在纯策略 Nash 均衡）；(iv) 相应的均衡映射是否具有 Monge 型结构（即是否可表为凸/凹函数梯度）。已有理论背景包含三方面：其一，经典 Monge–Kantorovich 理论（见 Ambrosio–Gigli [1], Brenier [2]）保证当成本为 $$ |x-y|^2 $$ 且源/靶测度绝对连续时，存在唯一 Monge 映射 $T = \nabla \phi$ ， $\phi$ 凸；其二，Gromov–Wasserstein 距离 [26] 引入了传输计划间的二次耦合，但其为单目标最小化而非零和博弈；其三，Cabrera [4]–[5] 等对路径依赖拥塞成本的研究提供了动态框架，但未引入对抗性优化结构。本文试图填补的关键空白在于：建立路径空间上双线性极小极大问题的完备存在性理论；揭示其与有限维端点问题的精确降维关系；证明在临界相互作用强度 $\alpha < \pi^2 / 2$ 下，Nash 均衡传输计划（NETP）诱导的映射满足 Brenier 型结构——即存在 $\phi, \psi$ 使得 $T_1 = \nabla \phi$ （当 $\alpha < 3$ ）或 $T_2 = \nabla \psi$ （当 $\alpha > 3$ ），从而将零和博弈、最优传输与 Monge–Ampère 型 PDE 系统（见后文）首次系统性地桥接起来。

2. 主要定理或主要结果

定理1.1（存在性）：设 $\mu_1, \nu_1, \mu_2, \nu_2 \in \mathcal{P}(\mathbb{R}^d)$ 具有紧支集，且成本 $$ c $$ 满足条件1（强制性与上下半连续性），则极小极大问题(1)至少存在一组解 $(\pi_1^*, \pi_2^*)$ ，即存在 $\pi_1^* \in \Pi_{\text{path}}(\mu_1,\nu_1)$ 、 $\pi_2^* \in \Pi_{\text{path}}(\mu_2,\nu_2)$ 使得

\sigma_0(\pi_1^*, \pi_2^*) = \inf_{\pi_1} \sup_{\pi_2} \sigma_0(\pi_1, \pi_2).

定理1.2（极小极大等于极大极小）：在条件2（二次相互作用）下，存在 Nash 均衡传输计划（NETP） $(\pi_1^*, \pi_2^*)$ ，满足

\sigma_0(\pi_1^*, \pi_2^*) = \inf_{\pi_1} \sup_{\pi_2} \sigma_0(\pi_1, \pi_2) = \sup_{\pi_2} \inf_{\pi_1} \sigma_0(\pi_1, \pi_2).

定理1.3（路径到端点约化）：若 $(\pi_1^*, \pi_2^*)$ 是路径空间问题(1)的解，则其支撑位于 $$ c $$ -极小极大路径上，且端点投影 $(e_0,e_1)_\# \pi_i^*$ 解决了由有效端点成本 $c_e^{\pi_1^*,\pi_2^*}(x,y,x',y')$ 定义的 Nash–Monge–Kantorovich 问题（见公式(40)）。定理1.4（Brenier 型结构）：设成本为

c(x,y,x',y') = \frac{1}{2}|x-y|^2 - \frac{1}{2}|x'-y'|^2 + \frac{\alpha}{3} \left[ (x-x')^2 + (x-x')(y-y') + (y-y')^2 \right],

其中 $\alpha < \pi^2 / 2$ ，且 $\mu_i, \nu_i \in \mathcal{P}_2(\mathbb{R}^d)$ 绝对连续，则当 $\alpha \neq 3$ 时，端点 Nash 均衡问题有唯一解，且存在传输映射 $T_1: \mathbb{R}^d \to \mathbb{R}^d$ 、 $T_2: \mathbb{R}^d \to \mathbb{R}^d$ 满足 $(T_i)_\# \mu_i = \nu_i$ ；进一步，若 $\alpha < 3$ ，则 $T_1 = \nabla \phi$ ， $\phi$ 严格凸；若 $\alpha > 3$ ，则 $T_2 = \nabla \psi$ ， $\psi$ 严格凹。此外，若 $(\pi_1^*, \pi_2^*)$ 是端点成本 $c_e^{\pi_1^*,\pi_2^*}$ 的 NETP，则原始路径问题(1)的唯一解由耦合极小极大路径映射 $\Gamma_i: \operatorname{spt}(\mu_i) \to \Omega_i$ 诱导，即 $\pi_i^* = (\Gamma_i)_\# \mu_i$ 。定理1.5（对偶性）：若 $(\pi_1^*, \pi_2^*)$ 解决(1)，则其可通过耦合对偶问题刻画：存在势函数 $\phi, \psi \in C^2$ 使得 $\pi_1^*$ 集中于集合 $\{ (\gamma,\xi) : \phi(\gamma(0)) + \psi(\xi(0)) = c(\gamma,\xi) + \phi(\gamma(1)) + \psi(\xi(1)) \}$ ，且满足循环（反）单调性条件。

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

证明的核心障碍在于：(i) 路径空间 $\Omega_i$ 非紧致（仅具弱拓扑），导致标准极值原理失效；(ii) 双线性泛函 $\sigma_0$ 关于 $\pi_1$ 下半连续但关于 $\pi_2$ 上半连续，无法直接应用 Sion 最小最大定理；(iii) 从路径解向端点解约化需克服“路径信息丢失”问题，尤其当 $\alpha$ 接近临界值时，相互作用可能破坏端点唯一性。作者的真正创新点在于构造性地建立三层分析框架：第一层为路径空间上的变分约化。利用条件1的强制性，证明 $\Omega_M$ 在乘积度量下紧致，并结合 Prokhorov 定理与弱收敛性，在 $\mathcal{P}(\Omega_1) \times \mathcal{P}(\Omega_2)$ 上建立 $\sigma_0$ 的下半连续性（对 $\pi_1$ ）与上半连续性（对 $\pi_2$ ），再通过 Ky Fan 不动点定理（而非 Sion）证得 NETP 存在（定理1.2）。第二层为端点约化的显式构造。关键引理是：对任意 $\pi_i \in \Pi_{\text{path}}(\mu_i,\nu_i)$ ，定义有效端点成本

c_e^{\pi_1,\pi_2}(x,y,x',y') := \iint c(\gamma,\xi) \, d\pi_1^{x,y}(\gamma) \, d\pi_2^{x',y'}(\xi),

其中 $\pi_i^{x,y}$ 是固定端点 $\gamma(0)=x, \gamma(1)=y$ 的条件路径测度。作者证明当 $\alpha < \pi^2/2$ 时， $c_e^{\pi_1^*,\pi_2^*}$ 是良定义的、连续的且满足拟凸性，从而将无限维路径问题精确投影至有限维端点问题（定理1.3）。第三层为 Brenier 结构的 PDE 刻画。作者推导出均衡映射 $$ (T_1,T_2) $$ 必须满足耦合非线性椭圆系统：

\det D^2 \phi = \frac{\rho_1(x)}{\rho_2(T_1(x))}, \quad -\det D^2 \psi = \frac{\rho_2(x')}{\rho_1(T_2(x'))},

其中 $\rho_i$ 为密度函数，符号差异源于成本中动能项的正负号。此系统本质为 Monge–Ampère 型方程，但因 $\alpha$ 的符号与大小调控 Hessian 的正定性（ $\alpha < 3$ 时 $D^2\phi > 0$ ， $\alpha > 3$ 时 $D^2\psi < 0$ ），故需分别应用 Caffarelli 正则性理论与 Gutiérrez–Huang 的凹 Monge–Ampère 理论。此外，作者引入“循环反单调性”概念：对 $\alpha > 3$ ，最优耦合 $\pi_2^*$ 满足 $\sum_{i=1}^k c_e(x_i,y_i,x'_i,y'_i) \geq \sum_{i=1}^k c_e(x_i,y_{\sigma(i)},x'_i,y'_{\sigma'(i)})$ 对所有置换 $\sigma,\sigma'$ 成立，这直接保证了 $$ T_2 $$ 的凹性。

4. 与之前工作的比较

与既有工作相比，本文在三个维度实现实质性推进：(1) 与经典最优传输比较：Brenier 定理 [2] 要求单目标成本 $$ |x-y|^2 $$ 且无博弈结构；本文将 Monge 映射推广至零和博弈框架，证明当 $\alpha < 3$ 时 $T_1 = \nabla \phi$ （凸），当 $\alpha > 3$ 时 $T_2 = \nabla \psi$ （凹），且临界阈值 $\alpha = 3$ 对应 Hessian 符号切换，这是经典理论完全未涉及的新现象。代价是需更强的正则性假设（ $\mu_i,\nu_i \in \mathcal{P}_2$ 绝对连续）及交互强度限制 $\alpha < \pi^2/2$ 。(2) 与 Gromov–Wasserstein (GW) 距离比较：GW [26] 使用成本 $$ |d_X(x,x') - d_Y(y,y')|^2 $$ 并最小化单个耦合 $\pi \in \Pi(\mu,\nu)$ ；本文采用双耦合 $(\pi_1,\pi_2)$ 与零和结构，且成本含时间积分项，故 GW 的三角不等式与度量性质在此不适用，但获得了博弈均衡解释。(3) 与 Cabrera [4]–[5] 的路径拥塞模型比较：[4] 研究单目标最小化 $\int L(\gamma) d\pi(\gamma) + \int V(\gamma,\xi) d\pi(\gamma)d\pi(\xi)$ ，而本文将 $$ V $$ 提升为对抗性双线性项 $\alpha|\gamma-\xi|^2$ ，并证明其诱导的端点成本 $$ c_e $$ 满足 $\partial_{xx'}c_e = -2\alpha I_d$ ，从而在 $\alpha > 0$ 时保证 $$ c_e $$ 的“扭曲单调性”，这是 [4] 中未出现的结构性特征。局限性在于：本文未处理 $\alpha \geq \pi^2/2$ 的超临界情形（此时路径振荡加剧，端点约化失效）；亦未给出 $\phi, \psi$ 的显式构造（仅证明存在性），而 Brenier 定理可显式求解 $\phi$ 为凸包络。作者明确指出，与 Jeong [17] 的单计划双线性问题不同，本文的极小极大结构天然避免了非唯一性，但代价是需同时处理两个耦合测度的正则性。

8. Onsager型能量等式与Nernst–Planck流体系统的Prodi–Serrin唯一性

Onsager-Type Energy Equality and Prodi–Serrin Uniqueness for Nernst–Planck Fluid Systems

作者: Ruimeng Hu, Quyuan Lin, Qirui Peng

类别: math.AP

发布日期: 2026-06-30

arXiv ID: 2607.00225v1

摘要翻译: 我们研究耦合Nernst–Planck方程与流体模型的电扩散系统的弱解。首先，对于三维Nernst–Planck–Euler系统，我们建立了一个Onsager型判据，以保证耦合的动能-静电势能平衡的有效性。该能量等式对满足临界Besov正则性且具有消失二进制通量条件的速度场的弱解成立。此外，在假设离子浓度亦满足相应Onsager型正则性的前提下，我们还证明了抛物正则性、浓度非负性的保持性，以及相应的电荷密度能量恒等式。其次，对于三维Nernst–Planck–Navier–Stokes系统，我们证明了Leray–Hopf解的Prodi–Serrin型唯一性判据：当速度场属于Ladyzhenskaya–Prodi–Serrin类 $$L^p_tL^q_x$$ 且满足 $$2/p+3/q=1$$ 及 $$q>3$$ 时，Leray–Hopf类中的解是唯一的。这些结果将不可压缩流体动力学中的能量等式与弱-强唯一性原理推广至包含对流、扩散及自洽静电作用力的电扩散模型。

Comment: 26 pages

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

1. 研究对象和背景

本文研究三类耦合电扩散模型的弱解理论：Nernst–Planck–Euler（NPE）系统与Nernst–Planck–Navier–Stokes（NPNS）系统。其核心方程为

abla \cdot u = 0,

\partial_t c_i + u \cdot \nabla c_i = \nabla \cdot \left( D_i (\nabla c_i + z_i c_i \nabla \Phi) \right), \quad i = 1,\dots,N,

-\epsilon \Delta \Phi = \rho := \sum_{i=1}^N z_i c_i,

\partial_t u + u \cdot \nabla u - \nu \Delta u + \nabla p = - (k_B T_k) \rho \nabla \Phi,

其中 $c_i \geq 0$ 为第 $$i$$ 种离子浓度， $z_i \in \mathbb{Z}$ 为其价态， $$D_i > 0$$ 为扩散系数， $\Phi$ 为自洽电势， $$u$$ 为不可压缩流速场， $\nu \geq 0$ 为粘性系数。当 $\nu = 0$ 时即为NPE系统；当 $\nu > 0$ 时即为NPNS系统。文中聚焦于 $$N = 2$$ 、 $$z_1 = 1$$ 、 $$z_2 = -1$$ 、 $$D_1 = D_2 = D = 1$$ 、 $\epsilon = 1$ 、 $$k_B T_k = 1$$ 的对称情形，并引入电荷密度 $\rho = c_1 - c_2$ 与总浓度 $\sigma = c_1 + c_2$ ，从而将系统化为

\partial_t \rho + u \cdot \nabla \rho = \Delta \rho + \nabla \sigma \cdot \nabla \Phi + \sigma \Delta \Phi, \tag{1.1a}

\partial_t \sigma + u \cdot \nabla \sigma = \Delta \sigma + \nabla \rho \cdot \nabla \Phi + \rho \Delta \Phi, \tag{1.1b}

-\Delta \Phi = \rho, \tag{1.1c}

\partial_t u + (u \cdot \nabla) u + \nabla p = -\rho \nabla \Phi, \tag{1.1d}

\nabla \cdot u = 0. \tag{1.1e}

该系统定义在三维环面 $\mathbb{T}^3$ 上，初值 $(\rho_0, \sigma_0, u_0) \in L^2(\mathbb{T}^3)^3$ 满足 $\int_{\mathbb{T}^3} \rho_0 \, dx = 0$ ，且 $\Phi$ 满足零均值归一化。物理上，该模型描述电解质中离子在流体输运、扩散与自洽静电场共同作用下的演化，广泛应用于离子通道、微纳流控、半导体器件及生物膜电生理建模。数学上，其关键难点在于非线性耦合结构：一方面，Euler部分继承了湍流能量级联与Onsager临界正则性的深层障碍；另一方面，Nernst–Planck部分虽具抛物正则化效应，但其漂移项 $z_i c_i \nabla \Phi$ 含双线性乘积 $\rho \nabla \Phi$ ，而由 $-\Delta \Phi = \rho$ 导出的椭圆关系又使 $\nabla \Phi$ 仅比 $\rho$ 高一阶正则性——这导致经典能量恒等式

\frac{d}{dt}\left( \|u\|_{L^2}^2 + \|\nabla \Phi\|_{L^2}^2 \right) = -2 \int_{\mathbb{T}^3} |\Delta \Phi|^2 \, dx - 2 \int_{\mathbb{T}^3} \sigma |\nabla \Phi|^2 \, dx \tag{1.2}

在弱解框架下无法直接成立：对流项 $u \cdot \nabla u$ 、静电力项 $\rho \nabla \Phi$ 及漂移项 $\sigma \nabla \Phi \cdot \nabla \psi$ 均可能产生能量缺陷（energy defect），除非施加精细的正则性控制。已有工作如[16,18,29]建立了全局弱解存在性，但未解决能量守恒与唯一性问题；[2]在解析正则性下获得强适定性；[12,13]在二维或特殊参数下证明强解唯一性；而Onsager型能量等式此前仅见于纯流体系统（如Euler方程的 $B^{1/3}_{3,\infty}$ 准则[7,11]），尚未延拓至含自洽电动力学耦合的PDE系统。本文旨在填补这一空白：一方面，建立首个适用于NPE系统的Onsager型能量等式判据，揭示流速与离子浓度在临界Besov空间中的协同正则性如何抑制多重非线性导致的能量耗散；另一方面，在NPNS框架下发展Prodi–Serrin型唯一性理论，克服静电耦合力对相对能量估计的干扰。

2. 主要定理或主要结果

本文共得三项主定理。定理3.5（NPE系统的Onsager型能量等式）：设 $(u, \rho, \sigma)$ 是系统(1.1)在 $$[0,T]$$ 上的弱解（定义3.1），若流速满足

u \in L^3([0,T]; B^{1/3}_{3,\infty}(\mathbb{T}^3))

且满足消失二进制通量条件

\lim_{q \to \infty} \int_0^T \lambda_q \|\Delta_q u(t)\|_{L^3(\mathbb{T}^3)}^3 \, dt = 0, \tag{3.12}

则能量恒等式(1.2)成立，即

\|u(t)\|_{L^2}^2 - \|u_0\|_{L^2}^2 + \|\nabla \Phi(t)\|_{L^2}^2 - \|\nabla \Phi_0\|_{L^2}^2 = -2 \int_0^t \int_{\mathbb{T}^3} |\Delta \Phi|^2 \, dx d\tau - 2 \int_0^t \int_{\mathbb{T}^3} \sigma |\nabla \Phi|^2 \, dx d\tau.

该结论不要求 $\rho, \sigma$ 具备任何额外正则性，仅依赖于Definition 3.1中的 $$L^2$$ 弱解结构。定理3.6（离子非负性与电荷-总浓度能量等式）：若在定理3.5假设基础上进一步要求

\rho, \sigma \in L^3([0,T]; B^{1/3}_{3,\infty}(\mathbb{T}^3)),

则由 $\Theta = \frac{1}{2}(\sigma + \rho), \theta = \frac{1}{2}(\sigma - \rho)$ 定义的原始离子浓度满足 $\Theta(t,x) \geq 0$ , $\theta(t,x) \geq 0$ 对所有 $(t,x) \in [0,T] \times \mathbb{T}^3$ 成立；且 $(u, \rho, \sigma)$ 为抛物型弱解（即 $\rho, \sigma \in L^2([0,T]; H^1(\mathbb{T}^3))$ ），并满足电荷-总浓度能量恒等式

\frac{1}{2} \left( \|\rho(t)\|_{L^2}^2 + \|\sigma(t)\|_{L^2}^2 \right) - \frac{1}{2} \left( \|\rho_0\|_{L^2}^2 + \|\sigma_0\|_{L^2}^2 \right) = -\int_0^t \left( \|\nabla \rho\|_{L^2}^2 + \|\nabla \sigma\|_{L^2}^2 \right) ds - \int_0^t \langle \sigma, |\rho|^2 \rangle ds. \tag{3.13}

此处 $\langle \sigma, |\rho|^2 \rangle = \int_{\mathbb{T}^3} \sigma |\rho|^2 dx$ 为非线性耗散项。定理4.3（NPNS系统的Prodi–Serrin唯一性）：设 $(u^{(1)}, \rho^{(1)}, \sigma^{(1)})$ 与 $(u^{(2)}, \rho^{(2)}, \sigma^{(2)})$ 均为NPNS系统（即 $\nu > 0$ ）的Leray–Hopf型弱解，初值相同，若 $u^{(1)} \in L^p_t L^q_x$ 满足尺度不变条件

\frac{2}{p} + \frac{3}{q} = 1, \quad q > 3,

则二者在 $$[0,T]$$ 上恒等。特别地，当 $$p = 2, q = 6$$ 时，条件退化为 $u^{(1)} \in L^2([0,T]; L^6(\mathbb{T}^3))$ ，此为经典Ladyzhenskaya–Prodi–Serrin类中最弱允许情形。

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

本文方法的核心障碍在于双重非线性耦合破坏经典能量估计的闭合性：一方面，NPE系统中Euler部分的 $u \cdot \nabla u$ 与Nernst–Planck部分的 $\rho \nabla \Phi$ 、 $\sigma \nabla \Phi$ 均含临界正则性乘积，传统Sobolev嵌入失效；另一方面，NPNS系统中静电耦合力 $\rho \nabla \Phi$ 使相对能量泛函的耗散项出现交叉项 $\int (u^{(1)} - u^{(2)}) \cdot (\rho^{(1)} - \rho^{(2)}) \nabla \Phi^{(1)}$ ，无法被粘性项 $\nu \|\nabla(u^{(1)} - u^{(2)})\|_{L^2}^2$ 直接控制。作者的真正创新点有三：第一，发展了一套针对电扩散耦合系统的Littlewood–Paley通量分解框架。不同于纯Euler方程中仅处理 $\Pi^u_q = \int \mathrm{Tr}[S_q(u \otimes u) \cdot \nabla S_q u]$ ，本文定义了三个通量项：流体通量 $\Pi^u_q$ 、电-流耦合通量 $\Pi^1_q = \int [S_q(\rho u) \cdot S_q \nabla \Phi - S_q(\rho \nabla \Phi) \cdot S_q u]$ 、离子漂移通量 $\Pi^2_q = \int [(S_q \sigma)(S_q \nabla \Phi) - S_q(\sigma \nabla \Phi)] \cdot S_q \nabla \Phi$ 。关键突破在于利用椭圆关系 $-\Delta \Phi = \rho$ 将 $\Pi^1_q$ 重写为 $\int [S_q \rho \, S_q(u \cdot \nabla \Phi) - S_q(\rho u) \cdot S_q \nabla \Phi]$ ，再通过交换子估计与Bernstein不等式（引理2.1）将其控制为 $|\Pi^1_q| \lesssim (K * d^2)^{3/2}(q)$ ，其中 $d_j(s) = \lambda_j^{1/3} \|\Delta_j u(s)\|_{L^3}$ ， $$K$$ 为衰减核（公式3.22）。第二，提出“频率卷积主导”技术（公式3.25）：将 $\int_0^t |\Pi^u_q| ds$ 的估计转化为 $\sum_{j \geq -1} K(q-j) \int_0^t \lambda_j \|\Delta_j u\|_{L^3}^3 ds$ ，再结合消失二进制通量条件（3.12）与 $\ell^1$ 核 $$K$$ 的可和性，严格证明 $\lim_{q \to \infty} \int_0^t (|\Pi^u_q| + |\Pi^1_q| + |\Pi^2_q|) ds = 0$ 。第三，对NPNS唯一性，构造修正的相对能量泛函

\mathcal{E}(t) = \frac{1}{2} \|u^{(1)} - u^{(2)}\|_{L^2}^2 + \frac{1}{2} \|\nabla(\Phi^{(1)} - \Phi^{(2)})\|_{L^2}^2 + \frac{1}{2} \|\rho^{(1)} - \rho^{(2)}\|_{L^2}^2 + \frac{1}{2} \|\sigma^{(1)} - \sigma^{(2)}\|_{L^2}^2,

并证明其时间导数满足

\frac{d}{dt} \mathcal{E}(t) \leq C \left( \|u^{(1)}\|_{L^q}^p + \|\nabla \Phi^{(1)}\|_{L^6}^2 \right) \mathcal{E}(t) + \text{低阶项}.

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

摘要翻译

中文标题: 柱状域中拟线性椭圆方程非负解的单调性

摘要翻译: 我们考虑在柱状域中满足混合齐次 Dirichlet-Neumann 边界条件的 $$p$$ -Laplace 方程的弱解。假设右端项为正且局部 Lipschitz 连续，我们证明：对任意 $$p>1$$ ，任何正解关于 $$x_N$$ 方向单调递增。作为应用，我们证明了 Allen-Cahn 型方程的解是一维的，并得到了 Lane-Emden 型方程的一个 Liouville 型结果。

不相关论文（基本信息）

1. 路径规划与速度限制对交通排放建模中最优解的影响

Influence of Routing and Speed Limits on Optimal Solutions in Traffic Emission Modeling

作者: Marc-André Bach, Simone Göttlich, Alena Ulke

类别: math.OC, math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00944v1

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

摘要: We investigate the influence of routing strategies and speed limit policies on optimal solutions in traffic emission models. Building on a first-order macroscopic traffic model coupled with an advection-diffusion model, we formulate single- and multi-objective optimization problems to simultaneously maximize traffic efficiency and minimize air pollution. We compare three control scenarios: optimizing only the routing strategy, optimizing only the speed limit policy, and optimizing both simultaneously. Numerical experiments on a small road network demonstrate that speed limit policies consistently achieve larger reductions in emissions and greater gains in traffic efficiency than routing strategies. Multi-objective optimization reveals the trade-off between the two goals and confirms that including speed limits in the control set yields Pareto-optimal solutions that are strictly superior to those obtained by routing control only. Our results provide quantitative guidance for traffic management seeking to balance mobility and environmental objectives.

2. Ashbaugh–Benguria 互隙猜想（Dirichlet 特征值）

The Ashbaugh–Benguria reciprocal-gap conjecture for Dirichlet eigenvalues

作者: Yanyang Li, Quanyu Tang, Haiqi Zhang

类别: math.AP, math.SP

发布日期: 2026-07-01

arXiv ID: 2607.01135v1

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

摘要: We prove the Ashbaugh–Benguria reciprocal-gap conjecture for the Dirichlet Laplacian in every dimension $N\ge2$ . Specifically, if $Ω\subset\mathbb R^N$ is a bounded domain and

0<λ_1(Ω)<λ_2(Ω)\leλ_3(Ω)\le\cdots

are its Dirichlet eigenvalues, then

\sum_{i=1}^{N} \frac{λ_1(Ω)} {λ_{i+1}(Ω)-λ_1(Ω)} \ge \frac{N}{j_{N/2,1}^2/j_{N/2-1,1}^2-1},

where $j_{μ,1}$ denotes the first positive zero of the Bessel function $$J_μ$$ of the first kind of order $$μ$$ . We also characterize the equality case: equality holds precisely when $$Ω$$ agrees with a Euclidean ball up to a set of Sobolev $$H^1$$ -capacity zero. In particular, among bounded Lipschitz domains, equality holds if and only if $$Ω$$ is a Euclidean ball.

3. 对余维数为二情形下一种分数质量概念的再审视

Another look at a notion of fractional mass in codimension two

作者: Michele Caselli, Mattia Freguglia, Nicola Picenni

类别: math.DG, math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00810v1

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

摘要: We study a notion of fractional $$s$$ -mass for codimension-two currents on closed Riemannian manifolds, defined via energy minimization with a prescribed Jacobian constraint. We prove equi-coercivity and $$Γ$$ -convergence, with respect to the flat topology, of the $$s$$ -mass on general codimension-two currents. We also prove several additional results for fixed $$s$$ . We establish improved regularity for $$s$$ -harmonic maps that are minimizing among competitors with vanishing Jacobian and show that their singular set has Minkowski dimension at most $$n-3$$ . Moreover, we show that the $$s$$ -mass defined via weak linking, as recently introduced by the authors, agrees with the prescribed Jacobian formulation used here, clarifying the extent to which the $$s$$ -mass depends, or ultimately does not depend, on the way singularities are prescribed.

4. 闵可夫斯基范数空间中的布朗运动

Brownian motion in Minkowski normed spaces

作者: Shin-ichi Ohta, Marco Rehmeier, Kohei Suzuki

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

发布日期: 2026-07-01

arXiv ID: 2607.00800v1

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

摘要: A Minkowski normed space is the Euclidean space equipped with a (possibly asymmetric) uniformly convex and smooth norm, forming a particular class of Finsler manifolds. We construct a stochastic process with one-dimensional time marginal densities given by the fundamental solution to the nonlinear Finsler heat equation in Minkowski normed spaces. This process is constructed as a solution to a singular McKean–Vlasov stochastic differential equation and constitutes a nonlinear Markov process in the sense of McKean. Furthermore, we show that solutions to this stochastic differential equation are pathwise unique, and thus probabilistically strong solutions, though the equation has singular coefficients beyond the subcritical regime. Since our construction is a natural extension of the construction of standard Brownian motion from the standard heat kernel, we call this process \emph{Brownian motion in Minkowski normed spaces.} To the best of our knowledge, this is the first construction of stochastic processes associated with nonlinear heat equation in Finslerian spaces.

5. 测度空间中结构化种群模型的异步指数增长

Asynchronous exponential growth for structured population models in measure space

作者: Christian Düll, József Z. Farkas, Piotr Gwiazda, Anna Marciniak-Czochra

类别: math.AP, math.FA

发布日期: 2026-07-01

arXiv ID: 2607.00658v1

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

摘要: This paper studies the asymptotic behaviour of a structured population model on the space of nonnegative Radon measures. Such formulations naturally arise when solutions develop concentration phenomena or when the population is represented by discrete cohorts. Asynchronous exponential convergence of measure solutions towards a one-dimensional global attractor is established. While such results are classical in the $$L^1$$ setting, their extension to measure spaces requires different compactness and spectral arguments. We identify conditions under which the classical asymptotic behaviour persists in the space of Radon measures endowed with the flat metric, thereby extending the theory of asynchronous exponential growth beyond the classical $$L^1$$ framework.

6. 线性弹性力学中的反密度问题：从部分可达边界上的局部测量唯一确定密度

Inverse Density Problem for Linear Elasticity: Uniqueness from Local Measurements on a Partially Accessible Boundary

作者: Yu Jia, Chengyu Wu, Jiaqing Yang

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00639v1

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

摘要: We consider the inverse boundary value problem in an elasticity system. It is proved that the density function $$ρ$$ and its derivatives at the boundary can be uniquely determined from the local Cauchy data. Furthermore, if the density function is analytic, we can uniquely determine the internal buried objects, as well as the unknown boundary and the boundary conditions imposed on it. Our methods mainly based on a precise characterization for the principal part of the difference between a special first-order singular solution and the fundamental solution in the $$H^m$$ norm, and the blow-up property for the boundary Sobolev norms of the volume potential corresponding to the fundamental solution.

7. 度量图上的薛定谔问题

The Schrödinger problem on metric graphs

作者: Juliane Krautz, Jan-F. Pietschman

类别: math.AP

发布日期: 2026-07-01

arXiv ID: 2607.00655v1

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

摘要: We study the Schrödinger problem on metric graphs and its different formulations. Starting from a static version, we introduce an equivalent reformulation as entropic optimal transport and show $$Γ$$ -convergence towards static optimal transport. We then rigorously derive a Benamou-Brenier type dynamic version of the Schrödinger problem, thereby extending known results from ${\rm RCD}^*(K,N)$ -spaces. With this equivalence at hand, we conclude that the minimum values of the dynamic Schrödinger problem converge towards the squared Wasserstein distance, and minimizers converge to Wasserstein geodesics. We also extend the dynamic formulation to a more general class of initial and final data and show existence of solutions in this setting using the direct method. Lastly, we illustrate our analytical findings by a numerical investigation.

8. 微型尺度下的多层级选择与人口衰老

Demographic senescence as multi-level selection in miniature

作者: Ananda Shikhara Bhat, Hanna Kokko

类别: q-bio.PE, cond-mat.stat-mech, math.AP, math.PR, physics.bio-ph

发布日期: 2026-06-30

arXiv ID: 2607.00262v1

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

摘要: Multi-level selection and senescence do not at first sight have much in common. Here, we demonstrate that the emergent mortality patterns generated by demographic senescence can be understood as the product of multi-level selection. We formulate a two-level Moran type process and use its scaling limits to illustrate that a simple mathematical framework that models multi-level selection in group-structured populations also models damage accumulation patterns and resultant mortality curves in ageing organisms. To verbally make the connection, observe that defectors spread within a group consisting of cooperators and defectors; when groups compete against each other, defector-rich groups suffer, and between-group selection causes such groups to be systematically under-represented. Exactly analogously, senescing individuals accumulate damage to physiological sub-systems, and `damage begets damage’; individuals who are more damaged are more likely to die, hence damage-rich individuals are systematically under-represented in later age classes. Thus, emergent senescence patterns in complex, integrated organisms are formally equivalent to the patterns generated by a within-generation multi-level selection process in which intra-organismal sub-systems play the role of particles, organisms play the role of collectives, and selective disappearance plays the role of group selection.

9. 热带几何作为物理信息神经网络的受限架构：在非线性流-固耦合例子中的应用

Tropical Geometry as a Restricted Architecture for Physics-Informed Neural Networks: Applications in Nonlinear Fluid-Structure Examples

作者: Carla Valencia-Negrete, Cristhian Garay-Lopez, Marco Favela-Rodriguez, Alonso Andapia-Viveros

类别: math.AP, math.AC

发布日期: 2026-06-30

arXiv ID: 2607.00237v1

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

摘要: Nonlinear algebraic (polynomial) differential equations that govern fluid-structure interactions, such as those modeling vortex-induced vibrations, and shock waves, often lack analytical solutions, creating significant challenges to efficient prediction and control. While Physics-Informed Neural Networks (PINNs) offer a mesh-free numerical alternative, they frequently suffer from convergence stagnation when optimizing over chaotic landscapes or stiff singularities. This paper introduces a hybrid methodology that integrates tropical differential algebraic geometry with deep learning. Using tropical algebra, we algorithmically determine a hard constraint, which we use to restrict the neural network’s hypothesis space to the exact support of the valid formal power series solution. We establish a theoretical Valuation-Support equivalence between classical Briot-Bouquet indicial analysis and the fundamental theorem of tropical differential algebraic geometry, proving that tropical methods accurately identify singularity structures. Numerical experiments on the Van der Pol and Burgers’ equations demonstrate that embedding these tropical constraints directly into the network architecture drastically reduces the search space, overcoming optimization stagnation and improving both accuracy and convergence speed in non-homogeneous physical regimes.

10. 光滑有界区域中接触角平均曲率流的极小化运动格式的收敛性

Convergence of a minimizing movement scheme for contact-angle mean curvature flow in a smooth bounded domain

作者: Tokuhiro Eto, Jiwoong Jang

类别: math.AP

发布日期: 2026-06-30

arXiv ID: 2607.00094v1

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

摘要: This paper studies a Chambolle-type minimizing movement scheme for mean curvature flow with prescribed contact angle in a smooth bounded domain. The scheme is based on the capillary functional and the geodesic signed distance relative to the container, and yields a time-discrete level-set approximation. The main result asserts that, for every $$C^1$$ boundary function prescribing a strictly nondegenerate contact angle, the approximate solutions converge locally uniformly to the unique viscosity solution of the corresponding level-set mean curvature equation with oblique derivative boundary condition. This improves a previous convergence theorem, where the container was assumed to be convex and a curvature-type condition relating the tangential derivative of the prescribed contact-angle function to the principal curvatures of the container boundary was imposed. The main new ingredient is a uniform Lipschitz estimate for the solutions of the variational problems defining the scheme. This estimate is derived by applying a Bernstein-type argument to a suitable weighted gradient, rather than to the gradient itself, which rules out boundary maxima without relying on the previous curvature-type condition.

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

1. n维球坐标系下梯度与散度的显式公式#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 低正则性先验估计、爆破准则与自由边界理想磁流体动力学中方的自交奇点（含表面张力）#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. Chern-Simons 系统爆破解的局部唯一性与非退化性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. Hasegawa-Mima 方程的弱解与耗散解#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 具非线性扩散的动理学Fokker-Planck方程#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 具有常涡量的三维毛细液滴的刚性结果#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 一个极小极大双线性传输问题与Nash-Monge-Kantorovich映射#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. Onsager型能量等式与Nernst–Planck流体系统的Prodi–Serrin唯一性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 通过除法的广义刘维尔定理#

摘要翻译#

2. 扰动环形区域中热问题能量的Schauder空间形状分析#

摘要翻译#

3. 与 Euler-Maxwell-双流体系统相关联的快极限模型#

摘要翻译#

4. 从部分边界数据同时重构 Stokes-Darcy 系统中的多个未知量#

摘要翻译#

5. 柱状域中拟线性椭圆方程非负解的单调性#

摘要翻译#

不相关论文（基本信息）#

1. 路径规划与速度限制对交通排放建模中最优解的影响#

2. Ashbaugh–Benguria 互隙猜想（Dirichlet 特征值）#

3. 对余维数为二情形下一种分数质量概念的再审视#

4. 闵可夫斯基范数空间中的布朗运动#

5. 测度空间中结构化种群模型的异步指数增长#

6. 线性弹性力学中的反密度问题：从部分可达边界上的局部测量唯一确定密度#

7. 度量图上的薛定谔问题#

8. 微型尺度下的多层级选择与人口衰老#

9. 热带几何作为物理信息神经网络的受限架构：在非线性流-固耦合例子中的应用#

10. 光滑有界区域中接触角平均曲率流的极小化运动格式的收敛性#

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

1. n维球坐标系下梯度与散度的显式公式

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 低正则性先验估计、爆破准则与自由边界理想磁流体动力学中方的自交奇点（含表面张力）

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. Chern-Simons 系统爆破解的局部唯一性与非退化性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. Hasegawa-Mima 方程的弱解与耗散解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 具非线性扩散的动理学Fokker-Planck方程

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 具有常涡量的三维毛细液滴的刚性结果

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 一个极小极大双线性传输问题与Nash-Monge-Kantorovich映射

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. Onsager型能量等式与Nernst–Planck流体系统的Prodi–Serrin唯一性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 通过除法的广义刘维尔定理

摘要翻译

2. 扰动环形区域中热问题能量的Schauder空间形状分析

摘要翻译

3. 与 Euler-Maxwell-双流体系统相关联的快极限模型

摘要翻译

4. 从部分边界数据同时重构 Stokes-Darcy 系统中的多个未知量

摘要翻译

5. 柱状域中拟线性椭圆方程非负解的单调性

摘要翻译

不相关论文（基本信息）

1. 路径规划与速度限制对交通排放建模中最优解的影响

2. Ashbaugh–Benguria 互隙猜想（Dirichlet 特征值）

3. 对余维数为二情形下一种分数质量概念的再审视

4. 闵可夫斯基范数空间中的布朗运动

5. 测度空间中结构化种群模型的异步指数增长

6. 线性弹性力学中的反密度问题：从部分可达边界上的局部测量唯一确定密度

7. 度量图上的薛定谔问题

8. 微型尺度下的多层级选择与人口衰老

9. 热带几何作为物理信息神经网络的受限架构：在非线性流-固耦合例子中的应用

10. 光滑有界区域中接触角平均曲率流的极小化运动格式的收敛性