生成时间: 2026年01月07日 03:53:58

重点关注论文数量: 5

了解领域论文数量: 4

不相关论文数量: 12

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

1. 二维 Navier-Stokes 方程的向前自相似解

Forward self-similar solutions to the 2D Navier–Stokes equations

作者: Dallas Albritton, Julien Guillod, Mikhail Korobkov, Xiao Ren

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

发布日期: 2026-01-06

arXiv ID: 2601.03161v1

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

详细分析

中文标题: 二维 Navier-Stokes 方程的向前自相似解

研究对象和背景

本文主要研究二维不可压缩 Navier-Stokes 方程（NS）在全空间 $\mathbb{R}^2$ 上的向前自相似解 (forward self-similar solutions)。Navier-Stokes 方程描述了粘性流体的运动，其标准形式如下：

其中 $u(x,t)$ 是速度场，$p(x,t)$ 是压力，$\Delta$ 是拉普拉斯算子。该方程组在全空间 $\mathbb{R}^2$ 上定义。

背景与动机： Navier-Stokes 方程的适定性理论中，解的唯一性是一个核心问题。在三维情形下，Leray-Hopf 弱解已被证明是不唯一的（例如 Jia 和 Sverák 的工作）。然而，在二维情形下，经典的 Leray-Hopf 解在有限能量空间中是唯一的。本文关注的是一个特殊的解类——自相似解，它们在 Navier-Stokes 方程的缩放对称性下保持不变。这种缩放变换定义为：

如果解 $(u,p)$ 满足上述变换下的不变性，则它具有如下形式：

其中 $(U, P)$ 被称为自相似剖面（profile），仅依赖于变量 $y = x/\sqrt{t}$。将此形式代入原方程组，得到关于剖面的椭圆型方程组，称为 Leray 方程：

这里 $y = x/\sqrt{t}$。

研究对象的特殊性：

1. 初始数据的性质：自相似解对应于 $-1$-齐次（1-homogeneous）的初始速度场 $u_0(x) \sim 1/|x|$。这类数据在 $d=2$ 时具有局部无限的动能（Infinite energy locally），甚至在 $d=2$ 时，其涡度 $\omega_0 = \nabla^\perp \cdot u_0$ 可能是狄拉克测度或柯西主值分布，不具有局部可积性。这与 $d=3$ 的情形不同（3D 中 $1/|x|$ 数据局部能量有限），使得标准的能量估计方法失效。
2. 临界空间：标准的局部适定性理论通常在临界空间（如 $L^d$）中成立，但 $-1$-齐次数据位于临界空间的边缘，属于弱 Lebesgue 空间 $L^{d,\infty}$。在 2D 中，这类数据甚至比 $L^{2,\infty}$ 更差，因为它导致动能无限。
3. 非唯一性现象：本文旨在解决两个问题（Q1）和（Q2）：
   • 对于任意大的 $-1$-齐次初始数据，是否存在自相似解？
   • 这些解是否唯一？

文章通过构造性和数值模拟指出，在 2D 中也存在类似 3D 的解的分岔（Bifurcation）和非唯一性现象。

主要定理或主要结果

论文的主要结果包括存在性定理和数值观察到的非唯一性现象。

1. 存在性定理 (Theorem 1.1): 假设 $\alpha \in (0,1)$，且 $u_0 \in C^\alpha(\mathbb{R}^2 \setminus \{0\})$ 是一个 $-1$-齐次、散度为零的向量场。那么，存在一个自相似解 $u$ 满足：

• $u$ 是 2D Navier-Stokes 方程的解，且初始数据为 $u_0$。
• 其剖面 $U$ 是光滑的。
• 剖面 $U$ 与热半群作用下的初始数据之间的差满足如下衰减估计：

其中 $M = \|u_0\|_{C^\alpha(S^1)}$，$\langle y \rangle = (1+|y|^2)^{1/2}$。这意味着解在无穷远处很好地逼近其初始数据的齐次形式。

2. 非唯一性的数值观察 (Observations 1.2): 考虑特定的 $-1$-齐次初始速度场：

其中 $y^\perp$ 表示将向量旋转 $90^\circ$（即 $(y_1, y_2)^\perp = (-y_2, y_1)$），$\sigma$ 是振幅参数。数值模拟表明，当 $\sigma$ 在区间 $39.2 \lesssim \sigma \le 80$ 内时，对于给定的初始数据，存在三个不同的自相似解。这揭示了 2D Navier-Stokes 方程在无限能量类中解的非唯一性。

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

本文的证明策略主要依赖于** Leray-Schauder 不动点定理**，其核心在于建立一致的先验估计 (A priori estimates) 和证明非线性算子的紧性 (Compactness)。由于 2D 自相似问题中初始数据局部能量无限，标准的流体力学估计不适用，作者发展了一套针对 Leray 方程结构的特定分析工具。

a) 核心分析学工具

1. Leray-Schauder 不动点定理： 作者将求解 Leray 方程 (1.4) 的问题转化为 Banach 空间 $X_\alpha$ 上的不动点问题。

   • 空间设置：定义空间 $X_\alpha = \{ V \in L^\infty(\mathbb{R}^2) : \langle y \rangle^{1+\alpha} V \in L^\infty, \text{div } V = 0 \}$，范数为 $\|V\|_{X_\alpha} = \| \langle y \rangle^{1+\alpha} V \|_{L^\infty}$。该空间专门用于捕捉剖面 $U$ 在无穷远处衰减至 $O(|y|^{-1})$ 的行为，同时允许 $V = U - a$（$a$ 是初始数据的热半群近似）有更快的衰减。
   • 算子设定：将方程改写为 $V - \lambda \mathcal{K}(V) = 0$，其中 $\mathcal{K}$ 是非线性算子，$\lambda \in [0,1]$ 是参数。通过证明对于任意 $\lambda \in [0,1]$，解 $V$ 在 $X_\alpha$ 中一致有界（先验估计），并结合算子 $\mathcal{K}$ 的紧性和连续性，即可由不动点定理保证解的存在性。

2. 加权 $L^2$ 能量估计与 Hardy 不等式： 由于 $U$ 本身在 $L^2$ 中可能发散（衰减仅为 $|y|^{-1}$），作者使用加权积分估计来控制解的行为。这涉及到处理如下形式的积分：

   $$ \int_{\mathbb{R}^2} \langle y \rangle^{-2\gamma} |U(y)|^2 dy $$

   这种加权估计有助于在保持椭圆结构的同时，获得比纯 $L^\infty$ 估计更精细的控制。

3. 椭圆正则性理论： Leray 方程 (1.4) 是一个椭圆系统。作者利用了标准的椭圆理论（如 Calderón-Zygmund 估计）来提升正则性。例如，通过估计涡度 $\Omega = \nabla^\perp \cdot U$ 的梯度（即 $\Delta U$）的 $L^2$ 范数。

4. 抛物型局部平滑效应 (Local-in-space smoothing)： 在证明的最后阶段，为了将静态剖面 $U$ 的性质反馈回时间相关的解 $u(x,t)$，作者引用了 [41] 中发展的“局部在空间上的平滑效应”估计。这一工具用于控制解在环形区域（远离原点，但时间靠近 $t=0$）的行为，证明 $U$ 与 $u_0$ 之间的差值衰减更快（对应于 Hlder 连续性）。

b) 解决技术性难题的主要技巧

难题 1：2D 无穷能量与低频控制 在 3D 中，$H^1$ 嵌入到 $L^6$，利用能量估计可以控制解的大小。但在 2D 中，$\dot{H}^1$ 函数允许对数增长甚至 $1/|y|$ 的衰减，这使得仅靠 $\nabla U$ 的 $L^2$ 范数无法控制 $U$ 的点态行为（即 $|U|$）。此外，初始数据的涡度可能包含狄拉克质量，导致标准涡度估计失效。

应对技巧：

• 自相似 Bernoulli 压力 (Self-similar Bernoulli Pressure) $\Phi$： 这是本文最关键的新工具。作者定义了量：
  $$ \Phi := \frac{1}{2}|U|^2 - \frac{1}{2}y \cdot U + P \quad (1.15) $$
  它满足方程：
  $$ -(\Delta + \frac{1}{2} U \cdot \nabla) \Phi = |\Omega|^2 \quad (1.16) $$
  利用最大值原理，作者证明了 $\Phi(y) \le C(\|u_0\|_{C^\alpha(S^1)})$。这个结果非常强，因为它表明 $|U|^2$ 被常数控制（减去有界的 $y \cdot U$ 和 $P$）。由此可以推断出 $|U|$ 在某个半径的圆上是有界的，从而克服了 $H^1$ 嵌入失效的问题。
• 利用大附加项 $y \cdot \nabla_y U$： 在 3D 的分析中，$y \cdot \nabla_y U$ 是一个大项，但在 2D 中，作者巧妙地利用这个项来获得额外的衰减。作者观察到 $\partial_r(rU)$ 几乎属于 $L^2$，这使得可以通过加权积分 (1.18) 获得 $U$ 的强衰减。本质上，算子 $1 + y \cdot \nabla_y$ 的主导作用迫使 $U$ 在无穷远处表现出良好的衰减特性，即使初始数据仅属于临界空间。

难题 2：压力项的点态控制 如何在没有 $L^2$ 能量控制的情况下估计压力 $P$？

应对技巧：

• Poisson 方程与 Hardy 空间： 利用压力满足的 Poisson 方程 $\Delta P = - \text{div}(U \cdot \nabla U)$ 或 $\Delta P = U : \nabla \nabla U$（来自 (1.14)）。作者使用 Div-Curl 引理 和 补偿紧性 (Compensated Compactness) 理论，指出右端项属于 Hardy 空间 $\mathcal{H}^1$。由此得出结论，$\nabla P$ 在 $L^2$ 上有界，$P$ 在 $L^\infty$ 上有界。这是处理低正则性流体方程中压力项的经典且强力的工具。

难题 3：从静态剖面到时间演化的过渡 椭圆方程 (1.4) 的解 $U$ 是否真的对应一个时间演化 $u(x,t)$？即需要验证当 $t \to 0$ 时，$u(x,t) \to u_0(x)$ 的收敛率。

应对技巧：

• 自相似分解与平滑估计： 文章将解分解为 $U = a + V$，其中 $a = e^\Delta u_0$。对于 $V$，通过前面的椭圆估计得到其在 $X_\alpha$ 中的有界性（即 $|V(y)| \le C \langle y \rangle^{-1-\alpha}$）。这正是 (1.6) 的来源。最后，利用 [41] 中的局部平滑估计（Parabolic regularity）来处理 $a$ 在 $t \to 0$ 时的行为，将 $U$ 的性质提升为时间演化 $u$ 的性质，完成证明。

与之前工作的比较

本文的工作在 2D Navier-Stokes 理论中填补了重要空白，并与 3D 及凸积分方法的结果进行了对比。

1. 与 3D 自相似解构造 ([41]) 的对比：

• 维度的挑战：之前 Jia 和 Sverák 在 3D ([41]) 构造大自相似解时，关键步骤是利用 $u_0$ 的局部有限能量（因为 $1/|x|$ 在 $d=3$ 时局部可积）来获得先验估计。在 2D 中，这一条件完全失效（$1/|x|$ 局部不可积，动能无穷）。因此，本文不能照搬 3D 的证明。
• 方法的革新：本文引入的 自相似 Bernoulli 压力 $\Phi$ 是 3D 文献中未使用过的技巧。3D 的证明更依赖于 $H^1$ 嵌入和局部能量估计，而 2D 的证明必须依赖于这种特定的结构恒等式来“绕过”能量无穷的障碍。

2. 与临界空间适定性理论的对比：

• 通常，NS 方程在 $L^d$ 或 $L^{d,\infty}$ 空间有小数据全局存在性或大数据局部存在性。但 $-1$-齐次数据甚至不属于这些标准的临界空间（在 2D 中，它比 $L^{2,\infty}$ 更奇异，因为 $\frac{1}{|x|} \notin L^{2,\infty}_{loc}$）。
• 本文的结果表明，虽然在标准的流体框架下这些数据是“病态”的，但在自相似结构下，解仍然存在。这揭示了 Navier-Stokes 方程在临界边缘行为的复杂性。

3. 与非唯一性及凸积分 (Convex Integration) 结果的对比：

• 凸积分法：Buckmaster, Vicol, Cheskidov 等人利用凸积分证明了 Navier-Stokes 方程在各种空间中的非唯一性（包括 2D）。这些方法通常产生具有非常低正则性（如 $L^\infty_t L^p_x, p<2$）的弱解。
• 本文结果的区别：本文的非唯一性是通过分岔 (Bifurcation) 现象观察到的，发生在特定的自相似解族中。数值模拟显示，对于一个给定的“大”初始数据，随着参数 $\sigma$ 的变化，解的数量发生变化（从 1 个变为 3 个）。这与凸积分构造的大量互不相关的解不同，分岔现象揭示了解空间结构的内在不稳定性。
• 与 3D 数值证据的呼应：本文的数值策略复现了 3D [36] 中的发现（即 Pitchfork 分岔），首次在 2D 中提供了类似的证据。这表明非唯一性在 2D 无限能量类中也是存在的，打破了“2D Navier-Stokes 解总是唯一”的简单认知（注：经典 2D 有限能量解确实是唯一的，但本文处理的是更宽泛的无限能量类）。

总结： 本文的主要贡献在于克服了 2D 几何与无限能量初始数据带来的本质困难，通过引入创新的椭圆估计工具（Bernoulli 压力），证明了大自相似解的存在性，并通过数值模拟揭示了 2D 情形下的非唯一性分岔现象，将非唯一性理论从 3D 推广到了 2D 的临界情形。

2. 轴对称无旋情形下由Boussinesq方程控制的三维流体上的自由漂浮圆柱体

Freely floating cylinder on a 3D fluid governed by the Boussinesq equations in the axisymmetric without swirl case

作者: Geoffrey Beck, Ewan Contentin, Ludovic Martaud

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03133v1

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

详细分析

中文标题: 轴对称无旋情形下由Boussinesq方程控制的三维流体上的自由漂浮圆柱体

研究对象和背景

本文主要研究浅水波与漂浮结构物相互作用的一类数学模型。具体而言，研究对象是一个部分浸没在三维流体中的刚性圆柱体，该圆柱体被约束仅在垂直方向上运动（heave motion）。流体的运动由带有粘性项的Boussinesq型方程组控制，该方程组是对自由表面欧拉方程在非线性参数 $\varepsilon$ 和浅水参数（即色散参数）$\kappa$ 均为小量时的近似。

文章关注的物理场景是波浪能量转换器（Wave Energy Converters）附近的流体动力学，特别是“衰减测试”（decay test）情形，即结构物从初始位移释放后，在波浪的作用下回归平衡位置的过程。

1. 数学模型与方程组

在忽略粘性（$\nu=0$）和色散（$\kappa=0$）的情况下，浅水方程退化为双曲型的圣维南方程组。然而，为了捕捉沿海区域中波浪的色散效应，必须引入Boussinesq近似。文章考虑的是经过粘性修正的Boussinesq-Abbott方程组。在轴对称无旋（axisymmetric without swirl）的假设下，流体变量仅依赖于径向坐标 $r$ 和时间 $t$。设 $\zeta(t,r)$ 为水面升高，$Q(t,r)=q(t,r)e_r$ 为截面流量（discharge），$h(t,r)=1+\varepsilon\zeta(t,r)$ 为总水深，$P(t,r)$ 为表面压力。方程组写为：

在极坐标下，对于 $r>0$，上述方程化为径向Boussinesq-Abbott系统（BAr）：

其中 $d_r q = \partial_r q + \frac{q}{r}$ 是径向散度算子。

2. 边界条件与流固耦合

系统分为两个区域：

• 内部区域 $I$ ($0 \le r \le R$)：圆柱体下方。此处水深受圆柱体底部约束，满足 $h(t,r) = \delta(t) + h_{i,eq}$。这导致该区域的水面升高 $\zeta_i$ 与空间无关，且流量的散度满足 $\text{div} Q_i = -\dot{\delta}(t)$。此外，该区域存在非零压力 $P_i$，作为施加在圆柱体上的水动力升力。
• 外部区域 $E$ ($r > R$)：自由水面。此处压力 $P_e = 0$，水面升高 $\zeta_e$ 自由演化。

两区域在接触线 $\Gamma$（半径为 $R$ 的圆）处通过传输条件耦合：

1. 体积守恒：$[Q \cdot e_r]_\Gamma = 0$。
2. 能量/压力条件：$[\frac{P}{\varepsilon} + \frac{|Q|^2}{2h^2} + \zeta + \nu \dot{\zeta} + \kappa^2 \ddot{\zeta}]_\Gamma = 0$。
3. 无旋假设：$[\frac{Q}{h} \times n]_\Gamma = 0$。

圆柱体的垂直运动由牛顿第二定律描述（加去量纲化）：

这里的 $\tau_{\text{buoy}}^2$ 代表浮力因子，$\frac{1}{\varepsilon |I|} \int_I P_i \, dx$ 是水动力升力项，它依赖于流体的压力分布，从而构成了复杂的流固耦合系统。

3. 研究区域与函数空间

由于问题具有径向对称性，空间区域被简化为 $[0, \infty)$ 上的径向坐标。文章的主要分析是在无界域 $r \in [R, \infty)$ 上进行的。为了处理色散算子 $(1-\kappa^2 \partial_r d_r)$ 引入的非局部效应和边界层，文章引入了特定的函数空间，如 $W^{1,\infty}_{r,\kappa}$（包含 $\kappa \partial_r$ 项）和加权 Sobolev 空间 $H^1_{\kappa}$（包含 $\kappa^2 \|\partial_r d_r f\|_{L^2}$）。这些空间对于捕捉色散波在无穷远处的衰减行为至关重要。

主要定理或主要结果

文章提出了三个核心定理，分别对应于模型的重构、适定性分析和渐近行为研究。

定理 1.1 (增广公式, Augmented formulation) 在假设 1 成立的前提下，原流固耦合问题等价于以下的增广系统。设 $u = (\zeta_e, q_e)$ 为外部波场变量，$Z = (\delta, \dot{\delta}, \zeta_e, \dot{\zeta}_e)$ 为内部变量（注意 $\zeta_e$ 在 $r=R$ 处的迹）。系统如下：

• PDE部分：
  $$ \partial_t u + \partial_r \widetilde{F}_\kappa[u, Z] = S_r(u) + \nu \varepsilon S_\nu(u), \quad r > R $$
• ODE部分：
  $$ \partial_t Z + D_{\widetilde{}}(Z, I_{\text{hyd}}, F_{\text{ext}}) = 0, \quad t > 0 $$
• 耦合边界条件（在 $r=R$）：
  $$ u_1|_{r=R} = Z_3, \quad u_2|_{r=R} = -\frac{1}{2}Z_2 $$
  其中 $\widetilde{F}_\kappa$ 包含了非局部色散项，$S_r$ 是辐射阻尼项（2D效应），$S_\nu$ 是粘性耗散项。该定理的意义在于将原本依赖于全流场解的边界条件，转化为仅依赖于边界迹 $Z$ 和历史变量的局部耦合，从而使得问题适定性分析变为可能。

定理 1.2 (增广公式的适定性)

• (i) 一般存在性 ($\nu \ge 0$)：在 Banach 空间 $L^\infty_r \times W^{1,\infty}_{r,\kappa} \times \mathbb{R}^4$ 中，系统存在局部唯一解，存在时间为 $T_{\text{ODE}} \sim \frac{\kappa^2}{\varepsilon (1+\nu\kappa^{-1})^2}$。这里 $L^\infty_r$ 指 $L^\infty([R, \infty), rdr)$。
• (ii) 正则性更高的存在性 ($\nu=0$)：在 Sobolev 空间 $H^1_r \times H^2_{\kappa} \times \mathbb{R}^4$ 中，系统存在局部唯一解，存在时间为 $T_{\varepsilon, \kappa, R} \sim \frac{1}{2} \left( -(\alpha + \dots) + \sqrt{(\dots)^2 + 4\varepsilon^{-1}} \right)$，其中 $\alpha$ 依赖于 $R$。
• 爆破准则：如果最大存在时间 $T^* < \infty$，则必须满足 $\limsup_{t \to T^*} \left( \|\zeta_e(t)\|_{L^\infty} + \frac{|\dot{\delta}(t)|}{h_i(t)} + |F_{\text{ext}}(t)| \right) = \infty$。这表明只要水深 $h_i$ 不趋于零且外部力有界，解就不会爆破。

定理 1.3 (回归平衡与衰减估计) 考虑线性情形（$\varepsilon=0$），初始状态为 $(\zeta, q, \delta, \dot{\delta}) = (0, 0, \delta_0 \neq 0, 0)$。

• 无粘色散情形 ($\nu=0, \kappa>0$)：圆柱体中心质量的回归速度不能快于 $O(t^{-1/2+})$。
• 有粘或无色散情形 ($\nu>0$ 或 $\kappa=0$)：圆柱体中心质量的回归速度不能快于 $O(t^{-3/2+})$（在2D情形下）或 $O(t^{-1+})$（在1D情形下）。 该定理揭示了色散效应（$\kappa>0$）显著减缓了系统回归平衡的速度，阻止了指数衰减的发生。

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

该论文在处理非局部色散算子与自由边界问题的耦合时，展示了一系列高超的分析技巧。证明过程主要分为三个阶段：模型重构、局部适定性证明、线性衰减分析。

a) 核心分析学工具

1. 傅里叶变换与象征分析 (Fourier Transform & Symbol Analysis)： 在证明增广公式（Theorem 1.1）时，作者利用了Boussinesq算子的傅里叶变换特性。算子 $(1-\kappa^2 \partial_r d_r) \partial_t$ 在频域对应于乘子 $(1+\kappa^2 k^2) i \omega$（在全空间情形）。这种符号分析帮助作者从原始的三阶方程（关于时间）中解出水面升高 $\zeta$ 的二阶隐藏ODE，从而重构变量。这是处理高阶色散项的标准且有效的手段。

2. 加权 Sobolev 空间与能量估计 (Weighted Sobolev Spaces & Energy Estimates)： 为了证明Theorem 1.2，文章构造了特定的能量空间 $H^1_r$ 和 $H^2_{\kappa}$。$H^2_{\kappa}$ 包含了 $\kappa \|\partial_r d_r u\|_{L^2}$ 项，这与色散算子的能量守恒密切相关。在无粘 ($\nu=0$) 情况下的证明 (ii) 依赖于精细的能量估计。由于存在对流项（$\varepsilon$ 项）和色散项（$\kappa$ 项），能量估计需要平衡非线性增长和色散正则化。作者通过估计 $L^\infty_t H^1_x$ 范数和 $L^2_t H^2_{\kappa, x}$ 范数，结合Gronwall不等式来确定解的存在时间 $T$。

3. 不动点定理 (Fixed Point Theorems)： 局部适定性的证明通常遵循标准的Picard迭代或Banach不动点定理框架。作者定义了一个映射 $\Phi$，将一个候选解映射为由增广方程确定的新解。为了证明该映射是收缩的，必须利用上述的能量估计来控制解的增量。定理1.2中两种不同的存在时间结果（$T_{\text{ODE}}$ 和 $T_{\varepsilon, \kappa, R}$）对应了两种不同的范数控制策略。

4. 拉普拉斯变换与复分析 (Laplace Transform & Complex Analysis)： 在研究定理1.3（线性衰减）时，作者引入了关于时间变量的拉普拉斯变换 $\hat{\delta}(s)$。耦合的ODE/积分方程被转化为复平面右半平面内的全纯函数方程。通过分析该方程在 $s \to 0$（长时行为）和 $s \to \infty$ 处的渐近展开，可以确定衰减率。特别是，利用Hardy空间理论或复围道积分可以得到 $t \to \infty$ 时的 $O(t^{-\alpha})$ 行为。

b) 解决技术性难题的关键技巧

1. 增广公式 (The Augmented Formulation)： 这是本文最核心的技术创新。原始的流固耦合问题是“非局部”的，即边界上的压力条件依赖于整个流体域的历史（通过 $\zeta$ 和 $\ddot{\zeta}$ 的迹）。

   • 难点：标准的双曲方程初边值问题理论（如Kreiss对称化方法）依赖于边界条件的局部性和兼容性条件，但对于色散方程，边界层的存在使得兼容性条件变得复杂，且Kreiss方法通常要求色散项为小扰动或具有特殊结构，此处 $\kappa^2 \nabla \text{div}$ 算子破坏了标准对称化所需的对易性质。
   • 技巧：作者利用色散算子的性质，推导出一个关于 $\zeta$ 在边界 $r=R$ 处迹的二阶常微分方程（ODE）。这使得原本由流体方程决定的边界项，变成了由局部变量 $Z$ 及其导数控制的代数关系。这实际上是将无限维的边界控制问题约化为有限维的ODE-ODE/PDE耦合系统。这种“增广”处理使得系统在没有粘性时也能利用色散算子的正则化效应来获得局部适定性。

2. 色散边界层与非局部流 (Dispersive Boundary Layers & Non-local Flux)： Boussinesq方程可以重写为带有非局部通量和源项的守恒律形式。在证明Theorem 1.2的(i)时，作者利用了这种结构。虽然没有粘性，但色散项 $\kappa^2 \partial_r d_r \partial_t q$ 产生了一个随距离 $r$ 指数衰减的源项（边界层效应）。这种非局部效应使得通量 $\widetilde{F}_\kappa$ 即使在无粘情形下也表现出阻尼特性，从而允许在 $L^\infty$ 类型的空间中获得解，尽管存在时间 $T_{\text{ODE}}$ 随着 $\kappa \to 0$ 而趋于0，这与纯双曲情形（$\kappa=0$）的适定性理论不兼容，暗示了极限的奇异行为。

3. 结合复分析的渐近展开： 在证明衰减定理1.3时，关键技术在于处理积分方程中的记忆核（memory kernel）。作者通过拉普拉斯变换将积分核转化为解析函数，并通过分析该函数在原点附近的奇异性（极点和分支点）来提取时间渐近行为。文章指出，无粘色散情形下，$\delta$ 的衰减率受限于 $t^{-1/2}$，这与1D情形下的 $t^{-3/2}$ 不同，是因为2D几何导致的波浪扩散特性更强，且色散算子在2D径向坐标下的格林函数衰减特性决定了这一下界。

与之前工作的比较

该论文在多个方面扩展并改进了现有的关于流体-结构相互作用（FSI）及Boussinesq方程的数学理论。

a) 维度的推广与几何结构的引入

• 从1D到2D的跨越：此前的许多工作（如作者团队在 [5, 32] 中的研究）主要集中在一维（1D）设置，即沿水平轴的一维波导。本文将问题推广到了三维流体中的轴对称情形（2D空间）。这一推广并非简单的形式上的变化，而是带来了本质的数学困难。
• 几何效应：在2D径向几何中，波浪方程的色散项和辐射项表现出与1D截然不同的性质。例如，Theorem 1.3指出，在2D无粘色散情形下，衰减率为 $O(t^{-1/2})$，而在1D中通常认为是 $O(t^{-3/2})$。这是因为2D波浪具有圆柱形扩散（Cylindrical spreading），导致波幅随距离自然衰减（几何扩散），这种扩散效应与色散效应耦合，改变了回归动力学的长时间行为。论文证明了即使在有粘性的情况下，2D情形下也不能获得指数衰减，这与某些1D结果不同（虽然1D也有慢衰减，但几何源项的结构不同）。

b) 模型理论与近似精度的提升

• 超越圣维南方程：许多关于漂浮体的经典工作（如 [11, 23]）基于非线性浅水方程（Saint-Venant equations），即忽略了色散项 ($\kappa=0$)。本文明确引入了色散参数 $\kappa$，并刻画了其对适定性和衰减的影响。特别是，Theorem 1.2 详细分析了存在时间 $T$ 对 $\kappa$ 的依赖关系，揭示了当 $\kappa \to 0$ 时，基于增广公式的适定性结果会失效（时间趋于0），从而指出了双曲极限的奇异性。
• 包含粘性项：在 [5] 的基础上，本文引入了粘性项 $\nu$。虽然物理上这代表了一种简化的湍流或边界层模型，但数学上它为系统提供了额外的耗散，使得在 $\nu > 0$ 时能获得不同的存在时间估计（Theorem 1.2-(i) 中的 $T_{\text{ODE}}$ 依赖于 $\nu$）和更快的衰减率（Theorem 1.3）。

c) 适定性理论的创新：增广公式 vs. 对称化方法

• 避免Kreiss对称化器：对于双曲型方程的边界值问题，Kreiss对称化器理论是证明适定性的标准工具。然而，对于高阶的Boussinesq方程，特别是当色散算子不可对称化时，该方法难以直接应用。
• 增广公式的普适性：与之前的工作（如 [11] 对浅水方程使用对称化）相比，本文提出的“增广公式”提供了一种替代路径。它通过引入额外的边界变量 $Z$（包含 $\zeta_e$ 和 $\dot{\delta}$ 等），将问题转化为一个定义在外部区域的PDE与一个ODE在边界上的耦合。这种结构不仅绕过了对称化矩阵构造的困难，还自然地捕捉了由浮体运动引起的“记忆效应”（Added Mass Effect）。相比 John 模型（通常忽略接触面变化），本文的模型严格处理了接触面运动（通过 $\dot{\delta}$ 项），且不需要像 John 模型那样假设小振幅运动。

d) 衰减估计的精确化与改进

• 打破指数衰减的幻想：Theorem 1.3 的一个重要结论是，在无粘或有粘但有色散的情况下，漂浮体的回归运动不可能是指数快的。这在2D情形下尤为显著。之前的文献可能侧重于存在性，或者在1D情形下给出衰减估计。本文通过严格的复分析论证，给出了2D情形下的多项式衰减下界 $O(t^{-1/2})$。
• 对比无色散情形：文章指出，如果忽略色散（$\kappa=0$），衰减率会是 $O(t^{-3/2})$（2D有粘/无粘，或1D有粘）。这意味着色散效应实际上阻碍了能量的耗散，使得结构物回归平衡的过程更加漫长。这是物理直觉上的反直觉之处（通常认为色散会分散能量），但在数学上，这是由于色散波的守恒性质造成的，色散波包不像耗散波那样迅速消散，而是维持着与结构物的长期相互作用。

总结来说，这篇论文通过引入增广公式这一核心技巧，成功地将高阶色散项和复杂的流固耦合边界纳入了一个适定的数学框架，并在2D轴对称情形下得到了比前人更精细的存在性结果和衰减估计，是该领域理论分析的一个重要进展。

3. 三维无D条件等熵Navier-Stokes方程组的Liouville型定理

On Liouiville Type Theorem for the 3D Isentropic Navier-Stokes System without D-condition

作者: Quansen Jiu, Jie Tan, Zhihong Yan

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.02791v1

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

详细分析

中文标题: 三维无D条件等熵Navier-Stokes方程组的Liouville型定理

1. 研究对象和背景

本论文研究的对象是三维稳态等熵可压缩Navier-Stokes方程组（Compressible Navier-Stokes Equations）。该方程组描述了粘性流体（如液体或气体）在空间中随时间演化的运动状态，而“稳态”意味着流体的状态不随时间变化，即速度场 $u$ 和密度 $\rho$ 仅依赖于空间变量 $x$。

在流体力学中，Navier-Stokes方程组是描述粘性流体运动的核心方程。对于可压缩流体，方程组不仅包含动量守恒方程，还包含质量守恒方程（连续性方程）以及状态方程（联系压力与密度）。在 $\mathbb{R}^3$ 上，该系统的数学形式如下：

其中：

• $\rho = \rho(x) \ge 0$ 表示流体的密度。
• $u = (u_1, u_2, u_3)$ 表示流体的速度场。
• $P(\rho)$ 表示压力函数，满足状态方程 $P(\rho) = A \rho^\gamma$，其中 $A > 0$ 且 $\gamma > 1$（通常对于物理上的绝热气体，$\gamma \ge 1$，这里 $\gamma > 1$ 符合等熵假设）。
• $\mu > 0$ 和 $\mu + \lambda > 0$ 是粘性系数，分别对应第一和第二粘性系数。

Liouville型定理在偏微分方程理论中通常指：如果一个定义在整个无界空间 $\mathbb{R}^n$ 上的解具有某种衰减性或正则性（例如属于某个特定的 $L^p$ 空间），那么这个解必须是平凡的（即零解或常数解）。对于流体方程，这意味着流体处于静止状态（$u \equiv 0$）或者密度为常数。

背景与D条件： 该论文的一个核心背景在于对“D-condition”（Dirichlet-type regularity condition）的处理。在之前的许多Liouville型定理研究中（特别是对于可压缩情况），为了保证解的某些积分估计成立，研究者通常假设解满足 $\int_{\mathbb{R}^3} |\nabla u|^2 dx < \infty$，这被称为 D-condition 或 D-solution。这个条件本质上要求解的梯度属于 $L^2$ 空间，是一个较强的正则性衰减条件。

然而，对于无界区域上的稳态问题，直接验证 $u \in \dot{H}^1$ 或满足D条件往往很困难。本文的目标正是移除这一先验假设，即研究“无D条件”下的Liouville型定理。这意味着研究者必须在更弱的函数空间框架下，利用解的结构和不等式技巧来推导出解的平凡性。

函数空间的特殊性： 论文涉及的函数空间不仅包括经典的勒贝格空间 $L^p(\mathbb{R}^3)$，还深入探讨了Besov空间 $\dot{B}^{s}_{\infty, \infty}$ 和 Morrey型空间 $\dot{M}^{s,p}$。

• Besov空间 $\dot{B}^{s}_{\infty, \infty}$：用于描述解的“振荡”性质。如果一个函数不属于 $L^p$ 但具有某种精细的正则性，Besov范数可以捕捉到这种特征。
• Morrey空间 $\dot{M}^{s,p}$：这是一种局部化程度介于 $L^p$ 和 $L^\infty$ 之间的空间，用于控制函数在不同尺度下的积分行为。

本文研究的区域是全空间 $\mathbb{R}^3$，没有边界条件（边值条件为零），主要关注的是解在无穷远处的衰减行为如何强迫解在整个空间上消失。

2. 主要定理或主要结果

论文建立了稳态可压缩Navier-Stokes方程组的一系列Liouville型定理，主要结果概括如下：

定理 1 (基本衰减结果): 假设 $(u, \rho)$ 是方程组的光滑解，满足 $\rho \in L^\infty(\mathbb{R}^3)$（密度有界），且速度场满足积分衰减条件 $u \in L^p(\mathbb{R}^3)$，其中 $3 \le p \le \frac{9}{2}$。在此条件下，可以推出 $u \equiv 0$ 且 $\rho$ 为常数。

定理 3 (衰减-振荡条件): 假设 $(u, \rho)$ 是方程组的光滑解。

1. (端点情形) 如果 $u \in L^6(\mathbb{R}^3)$ 且动量 $\rho u \in BMO^{-1}(\mathbb{R}^3)$，则 $u \equiv 0$ 且 $\rho$ 为常数。
2. (中间范围) 如果 $u \in L^p(\mathbb{R}^3)$ 且动量 $\rho u \in \dot{B}^{\frac{3}{p} - \frac{3}{2}}_{\infty, \infty}(\mathbb{R}^3)$，其中 $\frac{9}{2} < p < 6$，则 $u \equiv 0$ 且 $\rho$ 为常数。

定理 5 (Morrey空间结果): 假设 $(u, \rho)$ 是方程组的光滑解，且满足以下Morrey型空间条件：

• 速度 $u \in \dot{M}^{s,6}(\mathbb{R}^3)$，参数 $2 \le s \le 6$。
• 动量 $\rho u \in \dot{M}^{q,3}_{w}(\mathbb{R}^3)$，参数 $\frac{3}{2} < q \le 3$。 在此条件下，推出 $u \equiv 0$ 且 $\rho$ 为常数。

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

本文的证明方法主要基于能量估计（Energy Estimates）和Caccioppoli型不等式的迭代技术，辅以针对不同函数空间特性的精细插值和交换子估计。核心难点在于处理动量方程中的非线性项（$\operatorname{div}(\rho u \otimes u)$ 和 $\nabla P(\rho)$），并在缺乏 $L^2$ 梯度控制的情况下，利用解的衰减性和振荡性来补偿。

a) 核心分析学工具

1. 截断函数与辅助函数构造 (Cut-off Functions and Auxiliary Functions): 所有的证明都始于选取一个径向截断函数 $\phi$。在全空间问题中，直接对整个空间积分会导致积分发散。因此，必须引入截断函数 $\phi$，其支撑集限制在大球 $B_r$ 上，且在小球 $B_\rho$ 上恒为1。

   • 关键技巧：为了处理动量方程中的压力项 $\nabla(\rho^\gamma)$ 和散度项 $\operatorname{div}u$，论文引入了一个辅助向量场 $w$。这个 $w$ 不是随意选取的，而是满足方程 $\operatorname{div}w = \gamma \nabla \phi \cdot u$。通过选择这样的 $w$，可以使得压力项的积分贡献精确抵消或转化，从而将能量估计转化为仅涉及 $u$ 的梯度和 $w$ 的积分不等式。

2. Caccioppoli型不等式 (Caccioppoli Inequality): 这是椭圆型和抛物型PDE正则性理论中的核心工具。它通常表示为：解在小球 $B_{\rho}$ 上的能量（如梯度的 $L^2$ 范数）可以被控制在稍大一点的球 $B_r$ 上的解本身（$L^p$ 范数）所控制。 论文的证明结构通常是推导如下形式的不等式：

   $$ \Phi_u(\rho) \le \frac{1}{4} \Phi_u(r) + C \cdot \mathcal{N}(u) $$

   其中 $\Phi_u(\rho) = \int_{B_\rho} |\nabla u|^2 dx$，$\mathcal{N}(u)$ 是包含 $u$ 的 $L^p$ 范数（或Morrey范数）的项。通过迭代这一不等式（令 $\rho = R/2, r=R$ 等），可以消去 $\Phi_u(r)$ 项，得到 $\Phi_u(R/2) \le \epsilon(R) \cdot \|u\|^k_{L^p}$。当 $R \to \infty$ 时，若 $\epsilon(R) \to 0$，则梯度的积分必须为零。

3. 函数空间理论与插值 (Function Space Theory & Interpolation):

   • Morrey空间与Lebesgue空间的联系：在定理5中，利用了Morrey范数和Lebesgue范数之间的关系。例如，$\|u\|_{L^6(B_R)} \le C R^{1/3 - 1/s} \|u\|_{\dot{M}^{s,6}}$。这使得可以通过 $R$ 的幂次增长来控制局部的 $L^6$ 范数。
   • 弱型估计与对偶性：处理 $BMO^{-1}$ 和 Besov 空间时，经常利用这些空间的对偶性质或其与弱 $L^p$ 空间 (Lorentz spaces) 的关系。例如，$BMO^{-1}$ 的对偶是 Hardy 空间 $\mathcal{H}^1$，但在实际估计中，往往将其转化为对分布的估计。

b) 解决技术性难题的技巧

1. 处理非线性项 $\operatorname{div}(\rho u \otimes u)$ 与 D-Condition 的移除: 这是本文最大的技术挑战。在经典的 $L^2$ 能量估计中，这一项通过分部积分会产生 $\int \rho |u|^2 u \cdot \nabla \phi$ 和 $\int \rho u \cdot (\nabla u \cdot u)$ 等项。这些项通常需要 $u$ 在 $L^3$ 或 $L^6$ 中有界，且最好有 $u \in \dot{H}^1$ 来控制余项。

   • 技巧：吸收与迭代。作者没有假设 $\int |\nabla u|^2 < \infty$，而是利用迭代法。通过估计 $|I_1|, |I_2|$ 等非线性项，将它们表示为 $C(R) \cdot \|\nabla u\|_{L^2(B_\rho)}^\alpha$ 的形式。随后利用能量不等式的结构，将 $\nabla u$ 项吸收回左边的主项中。
   • 特别是在定理1中，当 $p=3$ 时，$I_1 \sim \|u\|_{L^3}^3 / R$。如果 $u \in L^3$，随着 $R \to \infty$，这一项趋于0，从而迫使梯度项为0。对于 $3 < p \le 4.5$，通过 $L^p$ 到 $L^3$ 的插值或直接利用 $L^p$ 在大球上的衰减性质，同样使得余项随 $R \to \infty$ 消失。

2. 处理振荡条件：BMO 和 Besov 空间 (Theorem 3): 当 $u$ 的积分衰减较弱（例如 $p > 4.5$ 或 $p=6$）时，仅靠 $L^p$ 范数不足以抵消非线性项。此时必须引入对动量 $\rho u$ 的振荡控制。

   • BMO^{-1} 的表示：文章利用了 $BMO^{-1}$ 空间的结构特征。若 $\rho u \in BMO^{-1}$，则存在一个反称矩阵 $d \in BMO$ 使得 $\rho u = \operatorname{div} d$。利用 $d$ 的 BMO 性质，可以对含有 $\rho u$ 的项进行精细积分估计（利用高斯平均值性质），从而控制非线性项 $I_6$。
   • Lorentz 空间与弱范数：在 Morrey 空间结果（定理5）中，作者大量使用了 Lorentz 空间（弱 $L^p$）的插值不等式。例如，利用 $\|f\|_{L^{q,\infty}}$ 来控制由 Morrey 范数转化而来的局部弱范数。这种技术允许处理比经典 $L^p$ 更粗糙的函数。

3. 处理 $L^6$ 边界情况 (Endpoint Case): 当 $p=6$ 时，能量估计中的某些系数不再随 $R \to \infty$ 衰减（出现 $R^0$ 项）。此时，必须利用更强的振荡条件 $BMO^{-1}$。论文的证明展示了如何利用 $BMO$ 范数的对数性质来弥补 $L^6$ 估计中的亏缺，或者通过更精细的变量替换（如利用 $\rho u = \operatorname{div} d$ 并结合 Stokes 算子的估计）来获得所需的衰减。

4. 与之前工作的比较

该论文在可压缩Navier-Stokes方程的Liouville型定理领域做出了显著的推进，主要体现在对先验假设的放宽和函数空间的拓展上。

a) 对假设条件的改进：D-Condition 的移除

• 先前的工作 (Li et al., Chae):

  • D. Chae (2012) 在 $L^{3/2}$ 衰减条件下证明了Liouville定理，但要求解是 D-solution，即满足 $\int |\nabla u|^2 dx < \infty$。
  • D. Li 等人 (后续工作) 将衰减条件推广到 $L^{9/2}$，但这仍然基于 D-solution 框架。
  • 这意味着之前的结论依赖于解不仅在无穷远处衰减，而且其导数也有足够的能量衰减。这在物理上是一个较强的约束。

• 本文的突破:

  • 本文的 Theorem 1 明确指出，只需 $u \in L^p$ ($3 \le p \le 9/2$) 和 $\rho$ 有界，无需任何关于 $\nabla u$ 的积分假设，即可推出 $u \equiv 0$。
  • 技术原因分析：之前的 D-condition 往往是为了保证 $L^2$ 能量估计的收敛性。本文通过巧妙的截断函数构造和迭代法，将对 $\nabla u$ 的依赖转化为对 $u$ 本身的依赖。当 $R \to \infty$ 时，这种依赖项衰减得足够快，从而无需先验的 $L^2$ 梯度控制。

b) 结果的精细度：衰减-振荡条件的对应

• Incompressible vs. Compressible:

  • 对于不可压缩流体，Seregin (2012) 等人建立了著名的 $L^6 \cap BMO^{-1}$ 判据。
  • 对于可压缩流体，物理量是动量 $\rho u$ 而非速度 $u$。本文 Theorem 3 将不可压缩的 $BMO^{-1}$ 条件完美平移到了动量 $\rho u$ 上。
  • 这种迁移不仅是形式上的类比，而且在数学上是深刻的：因为连续性方程 $\operatorname{div}(\rho u) = 0$，动量天然具有散度为零的性质，这使得 $\rho u \in BMO^{-1}$ 成为一个自然的假设（因为 $BMO^{-1}$ 可以表示为散度的像）。

• Morrey 空间的推广:

  • 之前的定理主要集中在 $L^p$ 或 Besov 空间。
  • 本文引入 Morrey 空间 $\dot{M}^{s,6}$ 和 $\dot{M}^{q,3}_w$。这是一个重要的泛化。
  • 比较：$L^6$ 空间控制的是整体积分，而 Morrey 空间控制的是局部积分的加权平均。如果 $u \in \dot{M}^{s,6}$，意味着在无穷远处，$u$ 可以比 $L^6$ 允许的衰减更慢，但在局部尺度上仍受到控制。
  • 意义：定理5表明，即使解不具备全局的 $L^6$ 可积性，只要其局部振荡和衰减由 Morrey 范数控制，且动量具有适当的 Morrey 振荡，解依然是平凡的。这比单纯推广到 Lorentz 空间（如 $L^{6, \infty}$）更具一般性，因为 Morrey 空间包含了更多关于解在不同尺度下行为的信息。

c) 边界情况 ($p=6$) 和 $p \in (4.5, 6)$ 的处理

• 先前的局限：在 Li 等人的工作中，主要集中在 $p \le 9/2$ 或特定的 $L^3$ 情况。对于 $p$ 接近 6 的情况，由于临界指数的变化，标准的能量方法往往失效。
• 本文的提升：Theorem 3 (2) 填补了 $\frac{9}{2} < p < 6$ 的空白，并给出了临界指数 $\frac{3}{p} - \frac{3}{2}$ 下的确切 Besov 空间要求。这显示了作者对临界空间（Critical Spaces）理论的熟练运用。在流体力学中，这类结果通常被称为“临界正则性结果”，即在不损失正则性的前提下所能得到的最弱衰减条件。

综上所述，该论文通过引入更精细的函数空间工具（Morrey, BMO, Besov）和改进的能量估计技术，成功移除了可压缩Navier-Stokes方程Liouville定理中长期存在的D-condition限制，显著推广了前人的结果。

4. 时离散广义欧拉方程中的质量分裂与多边际最优传输中的非蒙日解

Mass splitting in the time-discrete generalized Euler equations and non-Monge solutions in multi-marginal optimal transport

作者: Gero Friesecke

类别: math.AP, math.OC

发布日期: 2026-01-06

arXiv ID: 2601.02616v1

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

详细分析

中文标题: 时离散广义欧拉方程中的质量分裂与多边际最优传输中的非蒙日解

1. 研究对象和背景

这篇论文的核心研究对象是时离散、空间连续的广义欧拉方程（Generalized Euler Equations）及其与多边际最优传输（Multi-Marginal Optimal Transport, MMOT）之间的深刻联系。作者Gero Friesecke试图解决一个长期悬而未决的问题：在时离散化的框架下，广义欧拉方程是否存在质量分裂（Mass Splitting）现象，或者等价地，其对应的最优传输问题是否存在非蒙日（Monge）形式的解。

背景始于欧拉流体动力学方程。经典的不可压欧拉方程由以下方程组描述：

这里 $\Omega \subset \mathbb{R}^d$ 是流体占据的空间区域，$u$ 是速度场，$p$ 是压力。Arnold (1966) 提出了一个著名的几何观点，认为欧拉方程的解对应于**保体积微分同胚群（group of volume-preserving diffeomorphisms）**中的测地线。这可以表述为一个变分原理：

然而，Arnold 的变分原理在一般情况下是不适定的（ill-posed），即可能不存在最小值。为了解决这个问题，Brenier (1989) 引入了一个松弛形式，即广义欧拉方程。在这个松弛框架下，流体粒子不再必须遵循单一的轨迹，而是可以遵循一组轨迹构成的概率分布（即质量分裂）。物理上，这意味着允许流体质量在移动过程中“分裂”并在终点重新组合。

为了分析的方便，论文将连续时间问题离散化。设时间点为 $t_0, \dots, t_N$，将路径空间 $C([0, T]; \Omega)$ 替换为 $\Omega^{N+1}$。这导出了一个多边际最优传输问题。具体来说，论文关注的是如下离散变分问题（经过端点条件简化）：

约束条件为：

其中 $g_*$ 是给定的保体积映射（例如 $g_*(x) = -x$）。

论文的研究背景在于理解多边际最优传输中解的结构。如果最优测度 $\gamma$ 可以写成 $\gamma = (T_1, \dots, T_N)_{\sharp} \mu$，则称为蒙日解（Monge solution），对应于无质量分裂。如果 $\gamma$ 具有更复杂的结构（例如是多个蒙日映射的凸组合），则称为质量分裂。此前已知在连续时间极限下（即步长趋于0）存在质量分裂，但在固定的时离散设置下，由于边际固定为均匀测度 $\mathbf{1}_\Omega$，问题变得非常困难，长期未解。

2. 主要定理或主要结果

论文的核心贡献是解决了上述开放问题，证明了即使在时离散情形下，广义欧拉方程依然存在质量分裂。

主要结果（Proposition 4.1 及其推论）：

论文通过构造一个具体的反例（在空间维度 $d=1$，空间网格 $\Omega = \{-1, 0, 1\}$），证明了以下结论：

1. 存在性：对于广义欧拉方程的时离散多边际最优传输问题，存在**质量分裂（Mass Splitting）**的最优解。这意味着最优测度 $\gamma$ 不是蒙日形式的，即流体粒子在离散时间步之间分叉。
2. 机制揭示：在离散空间和时间的简单例子中，质量分裂的发生是因为存在一种“交换路径”的策略，该策略在能量上严格优于任何保持粒子一一对应的蒙日策略。

具体来说，对于端点映射 $g_*(x) = -x$（将流体倒置），在离散空间 $\{-1, 0, 1\}$ 上，论文证明了任何试图避免质量分裂的策略（蒙日策略）都必须携带“动能惩罚”，而允许质量分裂的策略可以通过让 $0$ 处的粒子在两步内完成跳跃并交换位置，从而避免这种惩罚。

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

这篇论文的研究方法主要依赖于变分法（Variational Methods）和构造性证明（Constructive Proof），而不是高深的泛函分析技巧。其核心在于通过极度简化的离散模型（Minimalistic Model）来揭示复杂的物理现象背后的机制。

a) 核心分析学工具

论文中最核心的工具是有限维变分原理和**组合优化（Combinatorial Optimization）**的直觉。由于作者将空间离散化为有限集 $\Omega = \{-1, 0, 1\}$，无限维测度空间的问题转化为了有限维概率单纯形上的最优化问题。研究者不需要使用复杂的紧性定理（如弱收敛下的Prokhorov定理）或不动点定理，而是直接计算目标函数值并比较不同策略（Plan）的优劣。

使用的具体工具包括：

1. 测度分解与延拓：利用 Lemma 2.1，将涉及 $N+1$ 个边际的路径空间测度 $\gamma$ 简化为涉及 $N$ 个边际的中间状态测度 $\bar{\gamma}$。这使得分析集中在粒子在每一步的分布上。
2. 边际约束（Marginal Constraints）：约束条件 $(\pi_i)_{\sharp} \bar{\gamma} = \mathbf{1}_\Omega$ 意味着在每个时间步，粒子必须均匀分布。在有限离散空间中，这意味着概率质量必须均匀分配。例如，在 3 个点的空间中，每个点的总概率为 $1/3$。这使得可行解空间结构非常清晰。
3. 能量泛函的直接比较：通过计算具体的动能 $\sum \frac{|x_{i+1}-x_i|^2}{\Delta t}$ 来比较不同策略。这种方法在处理非凸、非光滑问题时非常直观有效。

b) 关键技巧：离散机制的透明化

论文的关键技巧在于构造了一个极其简单的例子，使得质量分裂的原因一目了然。

技巧细节（基于 Proposition 4.1）：

考虑空间 $\Omega = \{-1, 0, 1\}$，时间步长 $\Delta t = 1$，目标是将初始分布（均匀分布）在两步内变为 $g_*(x) = -x$（即 $1 \to -1, -1 \to 1, 0 \to 0$）。

• 策略 A：无质量分裂（蒙日策略） 我们要求每个粒子在每一步只能走一条确定的路径。

  • 初始状态：粒子位于 $-1, 0, 1$（概率各 $1/3$）。
  • 我们需要 $1 \to -1$ 和 $-1 \to 1$，这通常需要两步。
  • 但是，$0$ 必须保持在 $0$。
  • 如果我们要在两步内完成 $1 \to -1$ 和 $-1 \to 1$，路径可能是 $1 \to 0 \to -1$ 和 $-1 \to 0 \to 1$。
  • 问题来了：在 $t=1$ 时刻，三个粒子 $1, 0, -1$ 都在 $0$ 点汇合？这违反了边际约束（$0$ 点必须有总质量 $1/3$，而这里汇集了 3 份 $1/3$ 的质量，即总质量为 1）。所以我们不能让所有粒子在中间时刻都聚集在 0。
  • 另一种尝试：$1 \to 0 \to -1$ (粒子A), $-1 \to 0 \to 1$ (粒子B), $0 \to 0 \to 0$ (粒子C)。这在 $t=1$ 时，粒子 A, B, C 都在 0。总质量为 1。不行。
  • 必须存在一种“交换”或者“分裂”。如果我们不允许分裂，粒子 $1$ 和 $-1$ 无法在保持 $0$ 粒子不动的情况下完成位置交换。如果 $0$ 粒子移动（例如 $0 \to 1$），虽然可能满足边际，但 $0$ 粒子的目标是 $0$，移动它会产生额外的动能，且终点不对。

  论文中的论证指出：任何非分裂的策略都面临着无法同时满足端点映射 $g_*$ 和中间边际约束的困境，或者被迫引入巨大的动能（让粒子跑得更远）。

• 策略 B：质量分裂（最优策略） 论文展示了如下机制是严格优越的：

  • 让 $0$ 点的粒子产生“分裂”：在 $t=1$ 时刻，一半质量去 $1$，一半去 $-1$。
  • 让 $1$ 点的粒子去 $0$，让 $-1$ 点的粒子去 $0$。
  • 这样在 $t=1$：
    • 来自 $1$ 的质量在 $0$。
    • 来自 $-1$ 的质量在 $0$。
    • 来自 $0$ 的质量分成了两半，分别在 $1$ 和 $-1$。
  • 结果是：$0$ 点聚集了 $2/3$ 的质量（来自 $1$ 和 $-1$），而 $1$ 和 $-1$ 各聚集了 $1/6$ 的质量（来自 $0$ 的分裂）。这满足了均匀分布的边际约束（$1/3$ 在每一点）！
  • 第二步：来自 $1$ 的粒子（在 $0$）走到 $-1$，来自 $-1$ 的粒子（在 $0$）走到 $1$。来自 $0$ 的两部分粒子分别停留在 $1$ 和 $-1$。
  • 最终，起点 $1$ 的粒子到了 $-1$，起点 $-1$ 的粒子到了 $1$，起点 $0$ 的粒子分列 $1$ 和 $-1$。咦？这终点不对，起点 $0$ 的粒子应该在 $0$。
  • 修正后的分裂机制（论文实际含义）： 质量分裂允许路径的交换。例如，让起点 $1$ 的粒子走路径 $1 \to 0 \to -1$，起点 $-1$ 的粒子走路径 $-1 \to 0 \to 1$。 但是中间时刻 $t=1$ 都在 $0$，这就要求起点 $0$ 的粒子必须离开 $0$ 去填补 $1$ 和 $-1$ 的空缺，以维持 $t=1$ 时的均匀分布。 起点 $0$ 的粒子可以走 $0 \to 1 \to 0$ 和 $0 \to -1 \to 0$。 这样，$t=1$ 时，$0$ 处有来自 $1$ 和 $-1$ 的粒子，$1$ 处有来自 $0$ 的粒子，$-1$ 处有来自 $0$ 的粒子。这满足了边际约束！ 这个策略产生了动能。论文证明存在某种分裂策略（实际上是允许质量重排的策略，或者真正的多对多映射）比任何一对一的策略动能更低。

核心洞察：该机制揭示了在离散设置下，质量分裂允许粒子在中间时刻“交换身份”，从而利用空间上的对称性降低总动能。这种机制在连续极限下通常被平滑化了，但在离散极限下是显而易见且占主导地位的。

4. 与之前工作的比较

这篇论文在多边际最优传输和广义欧拉方程的研究中占据了一个独特的位置，主要体现在它连接了连续极限与离散计算，并解决了特定的开放问题。

a) 突破与改进

1. 从连续极限到固定离散：

   • 先前工作：Brenier (1989) 和 BFS (2009) 已经证明了在连续时间模型（$N \to \infty$ 或 $\Delta t \to 0$）中存在质量分裂。
   • 本论文：指出对于时离散问题，情况并非显然。因为边际被严格固定为均匀测度 $\mathbf{1}_\Omega$，这在多边际传输中是一个非常强的约束（通常导致非紧性或非存在性）。先前的文献虽然研究了MMOT，但没有专门解决“固定均匀边际的时离散广义欧拉方程是否存在非蒙日解”这一具体问题。本论文证明了即使在这些刚性约束下，质量分裂依然发生。

2. 透明度的提升：

   • 先前工作：Brenier 和 BFS 的证明依赖于对连续路径空间的深刻分析和对偶性理论（Kantorovich Duality）。这种证明是存在性的，但难以直观理解为什么粒子“想要”分裂。
   • 本论文：提供了一个完全离散、有限维的例子。这使得证明从复杂的测度论分析转变为简单的组合优化和变分计算。这种“透明机制”是之前文献所没有的。

b) 结果的强度与适用范围

1. 存在性 vs. 计算：

   • 之前的结果主要保证了连续模型中解的存在性。
   • 本论文的结果更贴近数值计算。因为现代最优传输算法（如作者在参考文献 [FP23] 中提到的算法）通常是在离散网格上运行的。本论文的结果表明，即使在粗糙的网格上，数值模拟捕捉到的非蒙日（质量分裂）解并非数值误差，而是方程本身的本质特征。

2. 对蒙日猜想的补充：

   • 在最优传输领域，一个中心问题是“在什么条件下最优传输是蒙日型的”（即没有质量分裂）。对于二次成本，当边际绝对连续时，$N=2$ 的情况是 Monge 的（Brenier 定理）。
   • 对于 $N \ge 3$，情况复杂。Agueh 和 Carlier 证明了对于 Wasserstein 巴黎中心成本（$L^2$ 巴黎中心），在一般条件下是 Monge 的。
   • 本论文展示了广义欧拉方程的动能成本（也是一种类似于 $L^2$ 的成本，但具有时间依赖性和端点约束）却会导致非 Monge 解。这丰富了我们对多边际传输解结构的理解，说明即使在简单的二次成本下，边际的特殊几何结构（这里是恒等式保持）和端点约束的耦合会导致复杂的解结构。

总结来说，这篇论文并非在寻找更强的正则性估计或解的存在性（这方面已有Brenier的奠基性工作），而是通过构造精巧的反例，澄清了离散化后模型的行为，证明了质量分裂在计算相关的时间离散层次上依然成立，并揭示了其背后的变分机制。这对于理解流体方程的数值模拟和多边际最优传输的非凸性具有重要意义。

5. 基于薛定谔方程的一维浅水方程数值模拟色散正则化方法

A Schrödinger-Based Dispersive Regularization Approach for Numerical Simulation of One-Dimensional Shallow Water Equations

作者: Guosheng Fu, Chun Liu

类别: math.NA, math.AP

发布日期: 2026-01-05

arXiv ID: 2601.02561v1

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

详细分析

中文标题: 基于薛定谔方程的一维浅水方程数值模拟色散正则化方法

研究对象和背景

本文主要研究了一维浅水方程（Shallow Water Equations, SWE）的数值模拟问题，特别是在处理包含干湿界面（wetting-drying interfaces）和真空状态（vacuum states）等复杂流体动力学场景时的鲁棒性问题。浅水方程是描述自由表面流动（如河流、海岸工程和地球物理流体动力学）的基础模型。尽管其形式相对简单，但在数值模拟中面临巨大挑战，特别是在水深 $h$ 趋于零时，方程会发生退化，导致标准的求解器失效。

论文首先介绍了标准的一维带底地形 $b(x)$ 的浅水方程组：

其中 $g > 0$ 为重力常数，$h(x,t) \ge 0$ 为水深，$u(x,t)$ 为垂向平均流速。

为了克服标准数值方法在处理干床或真空状态时的困难，本文引入了一种色散正则化（dispersive regularization）框架。该研究的核心在于通过在动量方程中引入三阶色散项，将原双曲型系统转化为一个可以通过马德隆变换（Madelung transform）与带有漂移项的非线性薛定谔方程（Nonlinear Schrödinger Equation, NLS）建立等价关系的系统。

文中从能量变分方法（Energetic Variational Approach, EnVarA）的角度推导了这一正则化过程。通过在总能量泛函中添加费雪信息（Fisher-information）能量项：

得到修正的量子势（量子波姆势）。结合质量守恒和动量平衡，作者导出了正则化后的浅水方程组：

其中 $\mu = -\frac{\varepsilon^2}{2} \frac{(\sqrt{h})_{xx}}{\sqrt{h}}$ 为化学势（或量子势）。

该系统的一个关键性质是其与非线性薛定谔方程的等价性。引入复值波函数 $\psi(x,t) = \sqrt{h} e^{i\phi/\varepsilon}$，正则化系统等价于：

该方程是一个带有底地形诱导漂移项的**散焦型（defocusing）**立方非线性薛定谔方程。

研究区域通常设定在有界域 $\Omega = [-L, L]$，并处理边界条件问题。由于研究关注小参数 $\varepsilon$ 下的半经典极限（semiclassical limit），解的振荡特性非常显著。为了处理开放边界，文中引入了海绵层（sponge layer）和复吸收势（Complex Absorbing Potential, CAP）技术。

主要定理或主要结果

这篇文章主要是一篇数值分析与建模的论文，其核心成果不是传统的解析存在性定理，而是关于数值逼近的准确性和鲁棒性的论证以及方程等价性的推导。文中没有以 “Theorem X.X” 命名的严格数学定理陈述，而是通过数值实验和理论推导展示以下核心结果：

1. 等价性结果：通过马德隆变换，证明了所提出的色散正则化浅水方程系统与带有底地形漂移项的散焦立方非线性薛定谔方程在数学上是形式等价的。这建立了流体力学模型与量子力学模型之间的桥梁。

2. 逼近误差估计：在亚临界（subcritical）且无激波形成的区域，当 $\varepsilon$ 较小时，基于薛定谔方程的正则化解与经典的浅水方程解之间的逼近误差为 $O(\varepsilon)$。这意味着随着正则化参数 $\varepsilon \to 0$，该方法能一致收敛于原问题的解。

3. 干湿界面的鲁棒性：一个非常重要的发现是，该方法在处理移动的湿干界面时表现出极强的鲁棒性。由于波函数的模平方 $h = |\psi|^2$ 自动保证非负性，该方法不需要像传统方法那样引入特殊的通量限制器（limiters）或正性保持技术，即可自然处理真空状态（$h=0$）。

4. 激波形成的局限性：文中指出（通过数值实验 3.1.2），当浅水方程产生激波时，薛定谔正则化方法不再收敛于经典的熵解，而是产生色散激波波列（dispersive shock waves）。因此，该方法适用于无激波（如稀疏波）或涉及干床的复杂界面问题，但不适用于强激波占主导的流场。

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

本文融合了偏微分方程理论、变分法、数值线性代数和谱方法。以下是其主要的技术细节：

a) 核心分析学工具

1. 能量变分方法 (EnVarA)： 这是本文推导正则化项的核心物理和数学工具。作者利用 Lagrangian 坐标系下的动能作用量和势能泛函，通过变分原理导出了力平衡方程。

   • 力的推导：通过引入流映射 $x(X,t)$ 和质量守恒关系 $\delta h = -\partial_x(h\xi)$，分别推导了惯性力 $f_{inertial} = -h(u_t + u u_x)$ 和保守力 $f_{conservative} = h \partial_x g(h+b)$。
   • 色散项的引入：通过泛函增广，引入费雪信息能量 $F_{disp}$，并计算其变分导数得到量子势 $\mu$，从而在力平衡中引入色散项 $h \mu_x$。这一过程保证了正则化后的系统仍然具有 Hamilton 结构和能量守恒律（$d/dt(E + F_{disp}) = 0$）。

2. 马德隆变换 (Madelung Transform)： 这是连接浅水动力学与薛定谔方程的数学桥梁。通过变换 $\psi = \sqrt{h} e^{i\phi/\varepsilon}$，将流体变量 $(h, u)$ 转化为复波函数 $(\text{Re} \psi, \text{Im} \psi)$。这一工具使得原本的非线性双曲型系统（可能面临奇异性）转化为半线性抛物型（色散型）方程，便于使用谱方法求解。

3. 半经典极限理论 (Semiclassical Limit Theory)： 文中提及了关于非线性薛定谔方程在 $\varepsilon \to 0$ 时的收敛性理论。这是理解该方法有效性的理论基础：在解保持光滑且远离真空的区域，NLS 的解会收敛于经典的双曲方程解（即欧拉-马德隆方程）。

b) 解决技术性难题的技巧

1. 数值振荡的处理 (Stiffness & Oscillations)： 在小 $\varepsilon$ 下，波函数 $\psi$ 具有 $O(1/\varepsilon)$ 频率的快速振荡。直接使用显式时间步进会导致严重的 CFL 限制。

   • Strang 分裂法 (Strang Splitting)：将时间积分分为两步：
     1. 非线性势项演化：$i\varepsilon \psi_t = (g|\psi|^2 + gb)\psi$。由于这是一个对角化的 ODE 系统（在 Gauss-Lobatto 节点上），可以精确积分（Exact integration），得到解析更新公式：
        $$ \psi(t + \Delta t/2) = e^{-i \frac{\Delta t}{2\varepsilon} (g|\psi|^2 + gb)} \psi(t). $$
        这一步不仅高效，而且在离散层面保持了 $|\psi|^2$ 的守恒性（即水深的非负性）。
     2. 线性色散项演化：$i\varepsilon \psi_t = -\frac{\varepsilon^2}{2} \psi_{xx}$。这一步使用隐式格式（如 Crank-Nicolson）求解，保证无条件稳定，并通过频域或预条件迭代解法高效处理。

2. 空间离散：谱元法 (Spectral Element Method, SEM)：

   • 为了在小 $\varepsilon$ 下准确捕捉波函数的高频振荡，必须使用高阶精度的空间离散。本文采用基于 Gauss-Lobatto 节点的谱元法（$k=1$ 即二次插值）。
   • 这种方法结合了有限元的几何灵活性和谱方法的高精度收敛性。质量矩阵是对角化的（ lumped mass matrix），进一步提高了计算效率。

3. 开放边界条件：复吸收势 (CAP) 与海绵层：

   • 为了在有限计算域 $[-L-\ell, L+\ell]$ 上模拟无限域问题，防止边界反射干扰内部解，作者引入了复吸收势 $-i\sigma(x)\psi$。这相当于在边界层添加人工阻尼。
   • 光滑截断函数：使用五次平滑阶梯（quintic smoothstep ramp）函数定义 $\sigma(x)$，使其在边界处平滑增长，避免引入非物理的反射或数值不稳定性。
   • 边界处理：在海绵层内，非线性步的更新变为 $e^{-i...}e^{-\sigma(x)\Delta t/2}$，引入了指数衰减。

4. 干湿界面的自然处理：

   • 传统浅水求解器在 $h \to 0$ 时面临 $u = q/h$ 的奇异性。
   • 本文方法通过 $h = |\psi|^2$ 恢复水深。即使 $h$ 在离散点上由于数值误差接近零，波函数 $\psi$ 只要模很小即可，不会出现除以零的错误。
   • 初始条件的构造：为了启动薛定谔求解器，作者将不连续的 Riemann 初始数据正则化为光滑的 tanh 轮廓，并据此构造相函数 $\phi_0$，使得 $u_0 = \phi_{0,x}$ 平滑逼近目标流速。这解决了从非光滑初始场启动高精度谱方法的问题。

与之前工作的比较

本文提出的框架相对于传统的浅水方程数值方法（如有限体积法、间断 Galerkin 法）以及标准的正则化方法有显著的改进和不同视角：

1. 克服干床处理的复杂性：

   • 先前工作：传统方法（如 [3, 7, 18, 26, 27]）在处理 $h=0$ 时，通常需要复杂的人工干预，例如：
     • 通量限制器 (Flux Limiters)：防止负水深出现。
     • 底地形重构 (Hydrostatic Reconstruction)：调整底地形以保持干湿界面处的稳定性。
     • 正性保持限制器 (Positivity-preserving Limiters)：强制离散解的非负性。
     • 这些技术通常具有“问题依赖性”（problem-dependent tuning），算法复杂且难以保证一般情况下的稳定性。
   • 本文工作：利用薛定谔框架，水深 $h$ 作为波函数模的平方自然非负 ($h = |\psi|^2 \ge 0$)。真空状态对应于 $\psi \approx 0$。因此，该方法完全消除了对特殊干湿处理算法的需求。这是方法学上的一个重大简化，特别是在处理移动干湿界面时。

2. 从“人工粘性”到“物理色散”的转变：

   • 先前工作：许多正则化方法通过添加人工粘性项（Artificial Viscosity）来平滑激波或防止数值爆炸，但这往往会抹平物理细节，且在干床处难以平衡。
   • 本文工作：引入的 $\varepsilon^2$ 色散项源于物理上的费雪信息能量，具有 Hamilton 结构和能量守恒性质。这种色散机制在光滑区域提供高精度逼近，同时在宏观尺度上保持了物理守恒律。

3. 数值算法的效率与稳定性：

   • 先前工作：对于高频振荡问题，显式方法通常受限于严格的 CFL 条件。
   • 本文工作：采用 Strang 分裂法，将非线性项精确求解，线性色散项隐式求解。这种分离使得时间步长 $\Delta t$ 不受 $\varepsilon$ 的严格限制（仅受空间分辨率限制），且利用谱元法实现了高空间精度。

4. 应用场景的特殊性：

   • 与传统的激波捕捉方法不同，本文的方法在激波存在时会产生**色散激波（振荡）**而非熵激波（单调跳跃）。这看似是缺点，但在某些应用中（如量子流体动力学或某些无激波的自然流动模拟），这可能提供更丰富的内部结构信息。
   • 文章明确界定了该方法的优势领域：亚临界、无激波、含干湿界面的流动。这为数值方法的选择提供了新的指导原则。

综上所述，本文提供了一种基于量子力学视角的流体力学求解策略，其核心优势在于利用复波函数的结构特性，优雅地解决了传统流体求解器中最棘手的非负性和奇异性问题，同时提供了一套高效的数值实现方案（谱元+Strang分裂）。

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

1. 平面上退化椭圆型方程的平均梯度局域化

Average gradient localisation for degenerate elliptic equations in the plane

作者: Thibault Lacombe

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03078v1

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

摘要翻译

中文标题: 平面上退化椭圆型方程的平均梯度局域化

摘要翻译: 我们考虑 $B_1\subset\R^2$ 上可能高度退化的椭圆型方程 $ \dv G(\nabla u)=0 $ 的 Lipschitz 解，其中 $ G\colon\R^2\to\R^2 $ 是任意连续严格单调向量场。我们证明了 $u$ 要么在 $0$ 点处为 $C^1$ 类，要么沿序列 $δ\to 0$ 的任意爆破极限 $v(x)=\lim \frac{u(δx)-u(0)}δ $ 满足 $ \nabla v\in \mathcal{D}\cap \mathcal{S} \text{ a.e} $。这里，$\mathcal{D}$ 和 $\mathcal{S}$ 可以粗略地解释为椭圆性从下方和上方退化的集合，即 $ \nabla G$ 和 $(\nabla G)^{-1}$ 的对称部分具有零特征值。这为 $H(\nabla u)$ 的连续性（对于任意在 $\mathcal{D}\cap \mathcal{S}$ 上消失的连续函数 $H$）提供了强有力的佐证。与前人同类研究相比，除了连续性和严格单调性外，我们没有对 $G$ 的结构做任何假设。

2. 具有几乎扇形算子的有理核分数阶演化方程：统一ABC与W动力学的预解框架

Rational-Kernel Fractional Evolution Equations with Almost Sectorial Operators: A Resolvent Framework Unifying ABC and W Dynamics

作者: Mohamed Wakrim

类别: math.AP, math.FA

发布日期: 2026-01-06

arXiv ID: 2601.02894v1

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

摘要翻译

中文标题: 具有几乎扇形算子的有理核分数阶演化方程：统一ABC与W动力学的预解框架

摘要翻译: 我们研究了由具有非奇异记忆的有理核时间算子驱动的分数阶演化方程，其中包括 Atangana-Baleanu-Caputo (ABC) 算子和一种广义 W-算子。这些算子的拉普拉斯符号特征在于并不一定属于经典的 Bernstein 类。分析是在几乎扇形算子的框架下进行的，这使得能够获得超出标准解析半群理论的预解估计。通过拉普拉斯变换技术和围道积分，确立了温和解的存在性、唯一性和时间正则性，从而导出了相关预解族的构造。我们开发了一个统一的预解框架，使得能够在 ABC 和 W 动力学之间进行精确比较，并阐明有理记忆核对衰减和光滑化性质的影响。包括分数阶扩散型方程在内的几个例子阐释了该抽象理论，并突出了非奇异记忆对长时间行为的影响。

3. 黎曼流形背景下的拟凸性

Quasiconvexity in the Riemannian setting

作者: Aurora Corbisiero, Chiara Leone, Carlo Mantegazza

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.02642v1

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

摘要翻译

中文标题: 黎曼流形背景下的拟凸性

摘要翻译: 我们针对光滑黎曼流形 $(M,g)$ 的切空间之间的线性映射构成的向量丛上的连续函数 $f$，引入了一种拟凸性概念。该概念自然地推广了经典的欧几里得定义。我们证明了，对于 $M$ 的任意有界开子集 $Ω$，该条件刻画了如下关联积分泛函在空间 $W^{1,\infty}(Ω, \mathbb{R}^m)$ 的弱$*$ 拓扑下的序列下半连续性：

4. 再论输运-交换子估计中的正则性

Another look at regularity in transport-commutator estimates

作者: Elias Hess-Childs, Matthew Rosenzweig, Sylvia Serfaty

类别: math.AP, math-ph, math.CA

发布日期: 2026-01-05

arXiv ID: 2601.02326v1

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

摘要翻译

好的，这是您要求的摘要翻译，严格遵循了学术性和准确性的要求：

中文标题: 再论输运-交换子估计中的正则性

摘要翻译: 我们关注输运速度场需要何种正则性才能控制Riesz型交换子。这类交换子的估计在分析具有成对Riesz相互作用的粒子系统的平均场极限和涨落时起着核心作用，我们首先对此进行回顾。我们的第一个新结果表明，通常对速度场梯度所做的 $L^\infty$ 假设一般不能放宽为BMO假设。我们对所有维数和所有满足 $-2< s<d$ 的Riesz奇异性（除了一维对数端点 $s=0$ 外）都构造了反例。在这一例外的端点处，这样的放宽是可能的，这一事实与Hilbert变换的经典Coifman-Rochberg-Weiss交换子界有关。我们的第二个结果揭示了相互作用势的奇异性与速度场所需正则性之间的权衡关系。粗略地说，如果希望得到由势能决定的自然能量半范数下的交换子估计，那么更平滑（奇异性更弱）的相互作用需要更强的速度场控制。我们为一类广泛的势能阐述了这一原理，并表明在亚库仑Riesz机制下，已知交换子不等式中出现的速度场正则性是最佳的。尽管有这些否定性的发现，我们的第三个结果表明，对于近李普希茨输运场，一个“有缺陷的”交换子估计仍然成立。这个有缺陷的估计是著名的Brezis-Wainger-Hansson不等式的一个推论，它使我们能够在平均场密度属于临界缩放Sobolev空间时证明收敛速率。

不相关论文（基本信息）

1. 次主控拟模与谱不稳定性

Subprincipal Controlled Quasimodes and Spectral Instability

作者: Pelle Brooke Borgeke

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03188v1

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

摘要: Here we explore, in a series of articles, semiclassical quasimodes u(h,b), approximative solutions P(h)u(h,b)\sim 0, depending on $0<h<1$, and on b, the subprincipal symbol. We study a pseudodifferential operator with transversal intersections of bicharacteristics, where the principal symbol has double multiplicity, $p=dp=0$, in a small neigborhood $Ω$. Because of this fact, we instead study the subprincipal symbol b, and we can conclude that we get transport equations depending on b where sign changes for the imaginary part of b give approximative solutions with small support. These modes are used to estimate spectral instability, or the pseudospectrum. We also investigate the possibility that we can factorize the model operator as $P(h)=h^2P_1P_2,$ in this way actually annihilating the subprincipal symbol, thus there is no condition for the imaginary part of b. In a follow-up article, we examine different cases for more complex operators with tangential intersections of bicharacteristics, thereby generalizing the findings here.

2. 论朗道-麦克斯韦方程熵产生的单调性

On the monotonicity of the entropy production in the Landau-Maxwell equation

作者: Côme Tabary

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03107v1

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

摘要: We study the homogeneous Landau equation with Maxwell molecules and prove that the entropy production is non-increasing provided the directional temperatures are well-distributed and the solution admits a moment of order $\ell$, for some $\ell$ arbitrarily close to $2$. It implies that for an initial condition with finite moment of order $\ell$, the entropy production is guaranteed to be non-increasing after a certain time, that we explicitly compute. This is the first partial answer to a conjecture made by Henry P. McKean in 1966 on the sign of the time-derivatives of the entropy. We also obtain algebraic decay estimates for the entropy production for large time; as well as a short-time estimate without moment assumptions.

3. 具有无穷守恒律的基尔霍夫方程

A Kirchhoff equation with infinite conservation laws

作者: Chiara Boiti, Renato Manfrin

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03095v1

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

摘要: We show here that the quasilinear Kirchhoff-Pokhozhaev equation

with $a\neq0$, admits conservation laws of all orders.

4. 涉及混合局部-非局部算子的指数临界 Choquard 方程基态解的存在性与集中性

Existence and concentration of ground state solutions for an exponentially critical Choquard equation involving mixed local-nonlocal operators

作者: Shaoxiong Chen, Min Yang, Zhipeng Yang

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.03060v1

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

摘要: We study the Choquard equation involving mixed local and nonlocal operators

where $\varepsilon>0$, $s\in(0,1)$, $0<μ<2$, $f$ has Trudinger–Moser critical exponential growth, and $F(t)=\int_{0}^{t}f(τ)\,dτ$. By variational methods, combined with the Trudinger–Moser inequality and compactness arguments adapted to the critical growth and the nonlocal interaction term, we prove the existence of ground state solutions and describe their concentration behavior as $\varepsilon\to0^{+}$.

5. 具有竞争性奇异非线性的非齐次非线性薛定谔方程

Inhomogeneous nonlinear Schrödinger equations with competing singular nonlinearities

作者: Elisandra Gloss, Kanishka Perera, Bruno Ribeiro

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.02909v1

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

摘要: We study nonlinear elliptic equations arising as stationary states of inhomogeneous nonlinear Schrödinger equations with competing singular nonlinearities. Working in a weighted Sobolev space that combines the homogeneous Sobolev space with a weighted Lebesgue term, we establish continuous and compact embeddings of Caffarelli–Kohn–Nirenberg type. These embeddings, together with a model that displays a natural scaling, allow us to apply the abstract critical point framework of Mercuri and Perera (2025), yielding a sequence of nonlinear eigenvalues for the associated problem. This scaling property leads to a classification of weighted power-type nonlinearities into subscaled, scaled, and superscaled regimes. Within this variational setting, we obtain broad existence and multiplicity results for equations driven by sums of weighted power nonlinearities, covering superscaled, scaled, and subscaled interactions, both in the subcritical and critical cases. We also provide a nonexistence result as a consequence of a Pohozaev-type identity. Finally, in the radial setting we employ improved radial CKN inequalities to enlarge the admissible embedding ranges. This yields strengthened radial versions of all our main results, including two-dimensional configurations with more singular weights, where no compact embeddings are available in the nonradial case.

6. W-算子：一个具有非伯恩斯坦符号的沃尔特拉分数时间算子

The W-Operator: A Volterra Fractional Time Operator with Non-Bernstein Symbol

作者: Mohamed Wakrim

类别: math.AP, math.FA

发布日期: 2026-01-06

arXiv ID: 2601.02876v1

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

摘要: We introduce a new two-parameter fractional time operator with Volterra structure, denoted by the W-operator, defined through a generalized Laplace symbol. The operator preserves the Caputo-type high-frequency behavior while allowing a controlled modification of the low-frequency regime through an additional parameter, leading to regularized memory effects. We develop a complete symbolic and Volterra theory, including explicit Prabhakar-type kernels, a left-inverse Volterra integral, and a fractional fundamental theorem of calculus. We show that the natural factorization of the Laplace symbol does not fit the classical Bernstein product mechanism and that the symbol is not a Bernstein function in general. Despite this non-Bernstein character, we establish well-posedness of abstract fractional Cauchy problems with sectorial generators by resolvent estimates and Laplace inversion, yielding a W-resolvent family with temporal regularity and smoothing properties. As an illustration, we apply the theory to a W-fractional diffusion model and discuss the influence of the modulation parameter on the relaxation of spectral modes.

7. 具有容量密度条件的区域上抛物型方程的整体 Hölder 可解性

Global Hölder Solvability of parabolic equations on domains with capacity density conditions

作者: Takanobu Hara

类别: math.AP

发布日期: 2026-01-06

arXiv ID: 2601.02863v1

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

摘要: We investigate the Cauchy-Dirichlet problem for linear parabolic equations in divergence form. Under mild assumptions on the source term and the domain, we prove the existence of globally Hölder continuous solutions. Notably, our results accommodate data exhibiting singularities nearly as critical as the inverse square of the distance from the boundary.

8. 极小极大平稳积分变分线的 Morse 指数

Morse index of min-max stationary integral varifolds

作者: Mitchell Gaudet, Talant Talipov

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

发布日期: 2026-01-06

arXiv ID: 2601.02860v1

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

摘要: We prove an upper bound for the Morse index of min-max stationary integral varifolds realizing the $d$-dimensional $p$-width of a closed Riemannian manifold.

9. 二阶椭圆方程解的高阶Hölder逼近

Higher order H{ö}lder approximation by solutions of second order elliptic equations

作者: Grigori Rozenblum, Nikolay Shirokov

类别: math.AP, math.CA

发布日期: 2026-01-06

arXiv ID: 2601.02859v1

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

摘要: For a given second order elliptic operation $\mathcal{L}$ in a domain $Ω\subset{\mathbb{R}}^\mathbf{N}$, $\mathbf{N}\ $, and a compact set $\mathbf{K}\subsetΩ$, order $\mathbf{N}$-$2$-Ahlfors-David regular, we define the space $\mathcal{H}^{\mathbf{r}+ω}_{\mathcal{L}}(\mathbf{K})$ of continuous functions $f(x),\, x\in\mathbf{K}$, admitting, for any $δ>0$, a local approximation in the $δ$-neighborhood of any point $x\in\mathbf{K}$, with $δ^{\mathbf{r}}ω(δ)$-error estimate, by solutions of the equation $\mathcal{L} u=0$. For such functions, we prove the existence of a global approximation $v_δ$ on $\mathbf{K}$ with the same order of error estimate, by a solution of the same equation in a $δ$-neighborhood of $\mathbf{K}$. A number of properties of these functions $v_δ$ and their derivatives are established.

10. 具有奇异应力函数的弹性碰撞新方法

New approach for elastic collisions with singular stress functions

作者: Toyohiko Aiki, Chiharu Kosugi

类别: math.AP, math-ph

发布日期: 2026-01-06

arXiv ID: 2601.02639v1

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

摘要: A collision of a rubber rod to a hard floor is regarded as a simple example of obstacle problems for elastic material. In this article we have proposed a new mathematical model for the collision phenomenon by applying beam equations with singular stress functions, which is investigated in our recent works. As in the works we have established a mathematical method to deal with the singular stress function. Here, we demonstrate the validity of our modeling through observation to the numerical results. Also, we present existence and uniqueness results of the model given as initial boundary value problems.

11. 溶瘤病毒疗法数学模型中的行波

Travelling Waves in a Mathematical Model for Oncolytic Virotherapy

作者: Negar Mohammadnejad, Thomas Hillen

类别: math.AP

发布日期: 2026-01-05

arXiv ID: 2601.02568v1

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

摘要: Oncolytic virotherapy (OVT) is a promising cancer treatment strategy in which engineered viruses selectively infect and destroy tumor cells. Motivated by the biological mechanisms underlying viral spread and tumor invasion into the tissue, we analyze a non-cooperative reaction-diffusion model capturing the invasion of tumor tissue by oncolytic viruses. Using carefully constructed upper and lower solutions together with Schauder’s fixed point theorem, we establish the existence of positive travelling-wave solutions. In particular, we identify a minimal wave speed value $\bar c$ such that positive travelling waves exist for all $c \geq \bar c$ . Our analysis also highlights parameter regions where the existence of travelling waves remains ambiguous, suggesting new mathematical questions about the propagation of viral treatments through tumor environments.

12. 快速扩散型趋化系统弱解的Hölder估计

Hölder estimates of weak solutions to chemotaxis systems of fast diffusion type

作者: M. Marras, F. Ragnedda, S. Vernier-Piro, V. Vespri

类别: math.AP

发布日期: 2026-01-05

arXiv ID: 2601.02528v1

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

摘要: We study a quasilinear chemotaxis system of singular type, where the diffusion operator is given by $Δu^m$ with $0<m<1$, corresponding to the fast diffusion regime, and where the chemotactic drift is nonlinear. Since Hölder continuity constitutes the optimal regularity class for weak solutions to the porous medium equation, we establish analogous regularity results for bounded solutions of parabolic–parabolic chemotaxis systems in this setting. The proof is based on a refined De Giorgi–Di Benedetto iteration scheme adapted to the coupled structure of the system. These results advance the understanding of the fine regularity properties of chemotaxis models with nonlinear diffusion, and demonstrate that the interplay between singular diffusion and aggregation exhibits a regularizing mechanism consistent with the porous medium paradigm.