生成时间: 2025年10月23日 03:39:30

重点关注论文数量: 4

了解领域论文数量: 2

不相关论文数量: 18

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

1. Classical Euler flows generate the strong Guderley imploding shock wave

经典欧拉流产生强古德莱伊内爆激波

作者: Giorgio Cialdea, Steve Shkoller, Vlad Vicol

类别: math.AP, physics.flu-dyn, 35L67, 76L05

发布日期: 2025-10-22

arXiv ID: 2510.19688v1

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

详细分析

中文标题: 经典欧拉流产生强古德莱伊内爆激波

1. 研究对象和背景

本文研究的是多维可压缩欧拉方程（compressible Euler equations）在径向对称情形下的解的演化行为，特别是从光滑、无激波的初始数据出发，如何自然地演化出著名的古德莱伊自相似内爆激波（Guderley imploding shock wave）。这一现象在物理上具有深远意义：它描述了一个强对称激波向中心点聚焦并最终坍缩的过程，在高能密度物理、惯性约束聚变以及超新星爆发等场景中扮演关键角色。

可压缩欧拉方程组

论文考虑的空间维度为 $d = 2$ 或 $3$，定义在空间区域 $\mathbb{R}^d$ 上的时间演化问题。原始形式的守恒律形式如下：

$$ \partial_t(\rho u) + \mathrm{div}(\rho u \otimes u + pI) = 0, $$
$$ \partial_t \rho + \mathrm{div}(\rho u) = 0, $$
$$ \partial_t E + \mathrm{div}((p + E)u) = 0, $$

其中：

  • $u : \mathbb{R}^d \times \mathbb{R} \to \mathbb{R}^d$ 是速度场，
  • $\rho > 0$ 是密度函数，
  • $E$ 是总能量（包含动能与内能），
  • $p$ 是压力，由理想气体状态方程闭合：
$$ p = (\gamma - 1)\left(E - \frac{1}{2}\rho |u|^2\right), \quad \gamma > 1. $$

这里的 $\gamma$ 称为绝热指数，是给定的常数，取值范围通常在 $(1,3]$ 内，对应不同的气体类型。

为了便于分析非等熵流动中的波传播结构，作者引入了伪熵平方根 $b = e^{S/2}$（$S$ 为比熵），并将声速重新标度为 $\sigma = \rho^\alpha b$，其中 $\alpha = \frac{\gamma - 1}{2}$。由此得到一组对称形式的欧拉方程：

$$ \partial_t u + u \cdot \nabla u + \alpha \sigma \nabla \sigma = \frac{1}{\gamma} \rho^{2\alpha} b \nabla b, $$
$$ \partial_t \sigma + u \cdot \nabla \sigma + \alpha \sigma \, \mathrm{div}\, u = 0, $$
$$ \partial_t b + u \cdot \nabla b = 0. $$

该系统保留了原欧拉方程的双曲结构，并更清晰地揭示了特征传播机制。

径向对称假设

由于目标是构造一个完全对称的内爆过程，作者假设所有变量仅依赖于半径 $r = |x|$ 和时间 $t$，且速度方向沿径向单位向量 $e_r = x / |x|$：

$$ u(r,t) = u(r,t)e_r, \quad \rho = \rho(r,t), \quad b = b(r,t), \quad \sigma = \sigma(r,t). $$

在此对称设定下，上述对称形式可进一步简化为一维偏微分方程组：

$$ \partial_t u + u \partial_r u + \alpha \sigma \partial_r \sigma = \frac{1}{\gamma} \rho^{2\alpha} b \partial_r b, $$
$$ \partial_t \sigma + u \partial_r \sigma + \alpha \sigma \left( \partial_r u + \frac{d-1}{r} u \right) = 0, $$
$$ \partial_t b + u \partial_r b = 0. $$

接下来，引入黎曼不变量（Riemann invariants）以分离不同特征方向的信息：

$$ z = u - \sigma, \quad w = u + \sigma, \quad b = b. $$

它们分别关联三个特征速度：

$$ \lambda_1 = u - \alpha \sigma = \frac{1 + \alpha}{2} z + \frac{1 - \alpha}{2} w, $$

$$ \lambda_2 = u = \frac{1}{2}(z + w), $$

$$ \lambda_3 = u + \alpha \sigma = \frac{1 - \alpha}{2} z + \frac{1 + \alpha}{2} w. $$

在内爆过程中 $u \leq 0$，主导的快声学波速是 $\lambda_1$，对应的变量 $z$ 成为主要动力学变量；而 $w$ 为次主导变量。

为进一步捕捉梯度奇异性的发展，作者采用微分黎曼变量（Differentiated Riemann Variables, DRV）：

$$ \mathring{z} = \partial_r z + \frac{1}{\gamma} \alpha \rho^\alpha \partial_r b, \quad \mathring{b} = \partial_r b, \quad \mathring{w} = \partial_r w - \frac{1}{\gamma} \alpha \rho^\alpha \partial_r b. $$

这些变量满足耦合的输运型方程组（见文中的 (1.9a)-(1.9c)），用于精确追踪解的正则性退化过程。

Guderley 自相似激波解

核心研究对象之一是 Guderley 解——首个被发现的“第二类自相似解”（self-similar solution of the second kind），其自相似指数 $\lambda > 1$ 不能通过量纲分析确定，而必须通过求解非线性本征值问题获得。

设自相似变量为：

$$ \xi = r (-t)^{-1/\lambda}, \quad t < 0, $$

则原始变量具有如下自相似结构：

$$ u(r,t) = (-t)^{(\lambda - 1)/\lambda} U(\xi), \quad c(r,t) = (-t)^{(\lambda - 1)/\lambda} C(\xi), \quad \rho(r,t) = R(\xi). $$

激波面固定在 $\xi = 1$，即物理空间中 $r = g(t) = (-t)^{1/\lambda}$。函数 $U(\xi), C(\xi), R(\xi)$ 满足一组常微分方程系统（ODEs），其边界条件由 Rankine-Hugoniot 跳跃条件给出：

$$ U(1) = -\frac{2}{\gamma + 1}, \quad C(1) = \sqrt{\frac{2\gamma(\gamma - 1)}{(\gamma + 1)^2}}, \quad R(1) = \frac{\gamma + 1}{\gamma - 1}. $$

此解的关键特性包括：

  • 在 $t \to 0^-$ 时发生中心坍缩（implosion singularity），速度趋于负无穷；
  • 激波强度随时间增强，局部马赫数保持 $O(1)$ 大小；
  • 数值上长期被视为强收敛激波的“普适吸引子”。

然而，长期以来，Guderley 解的存在性缺乏严格的 PDE 理论支持，直到最近 Jang-Liu-Schrecker [18] 才首次建立了其作为真实弱熵解的数学基础。

函数空间与边界条件特点

本文构造的初始数据属于 $C^{1,1/3}$ 正则性空间，即一阶导数 Hölder 连续，指数为 $1/3$。特别地，初始时刻 $t = T_{\mathrm{in}} < 0$ 的数据在核心区 $r < (-T_{\mathrm{fin}})^{1/\lambda}$ 处处于静止状态：

$$ (u_{\mathrm{in}}, \rho_{\mathrm{in}}, b_{\mathrm{in}}) \equiv (0, 1, 0), $$

表示前方为冷、静止、单位密度的理想气体。这种“静止核”（quiescent core）结构确保了激波前方介质不扰动，符合标准内爆模型的物理设定。

整个时空被划分为多个区域（如图3所示）：

  • 白色区域：光滑演化区；
  • 浅蓝色区域：静止态；
  • 黄色区域：精确匹配 Guderley 自相似解。

激波面 $\mathcal{S}$ 是分段定义的曲线：先由一段高正则曲线 $r = c(t)$ 从预激波点连接到弱激波态，再过渡到另一段 $r = s(t)$，最终与 Guderley 前端完全重合。

2. 主要定理或主要结果

定理 1.3（主定理）

对于固定的维度 $d \in \{2,3\}$ 和任意选定的匹配时间 $T_{\mathrm{fin}} < 0$，存在初始时间 $T_{\mathrm{in}} < T_*$ 和预激波时间 $T_* \in (T_{\mathrm{in}}, T_{\mathrm{fin}})$，以及相应的径向初始数据集 $(w_{\mathrm{in}}, z_{\mathrm{in}}, b_{\mathrm{in}})$ 在 $t = T_{\mathrm{in}}$ 时刻，满足以下性质：

  • 正则性：数据属于 $C^{1,1/3}(\mathbb{R}_+)$；
  • 静止核：当 $r < (-T_{\mathrm{fin}})^{1/\lambda}$ 时，$(u_{\mathrm{in}}, \rho_{\mathrm{in}}, b_{\mathrm{in}}) \equiv (0,1,0)$。

该初始数据生成唯一的径向欧拉解，其演化分为五个阶段：

  1. 光滑演化：在 $[T_{\mathrm{in}}, T_*)$ 内，解保持光滑，属于 $C^{1,1/3}$；

  2. 预激波形成：在 $t = T_*$ 时刻，于某点 $r_* > 0$ 处形成预激波（preshock），此时主导变量 $z$ 出现 $C^{1/3}$ 尖点，局部渐近行为为：

    $$ z(r, T_*) = z(r_*, T_*) + a(r - r_*)^{1/3} + b(r - r_*)^{2/3} + O(|r - r_*|), $$

    而 $w, b$ 保持更高正则性；

  3. 激波发展：在 $(T_*, T_{\mathrm{fin}}]$ 上，从预激波处诞生唯一、初始较弱的正规激波（regular shock），沿轨迹 $r = s(t)$ 传播，两侧解光滑；

  4. 与 Guderley 解匹配：在 $t = T_{\mathrm{fin}}$ 时刻，解及其激波前沿与经典的 Guderley 自相似剖面精确匹配至所有阶；此后 $[T_{\mathrm{fin}}, 0)$ 上解恒等于 Guderley 解；

  5. 最终坍缩：在 $t = 0$ 时，解在原点 $r = 0$ 发生内爆奇异性，$u \to -\infty$, $E \to +\infty$，但 $\rho$ 有界。

推论 1.4（全局存在性）

存在初始时间 $T_{\mathrm{in}} < 0$ 和属于 $C^{1,1/3}(\mathbb{R}_+)$ 的径向初始数据 $(u_{\mathrm{in}}, \rho_{\mathrm{in}}, E_{\mathrm{in}})$，使得对应的唯一欧拉解在整个时间区间 $[T_{\mathrm{in}}, \infty)$ 上存在：

  • 当 $t \in [T_{\mathrm{in}}, 0)$ 时，经历平滑→预激波→弱激波→强 Guderley 内爆；
  • 当 $t > 0$ 时，延续为反射的 Guderley 爆轰波（reflected blast wave），向外传播。

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

a) 核心分析工具

1. 反向构造法（Backward Construction / Backward-in-time Method）

这是整篇论文的核心策略：不是从前向后模拟激波演化，而是从目标状态——Guderley 解在 $t = T_{\mathrm{fin}}$ 的状态——逆时间回溯构造激波发展路径。

具体分为三步：

  • Step I：从 Guderley 强激波反推至弱激波（§4）；
  • Step II：从弱激波反推至预激波（§3）；
  • Step III：从预激波反推至光滑初值（§5）。

每一步都涉及奇异极限或非均匀尺度变换下的渐近分析。

2. 迭代映射与不动点定理

在每个反向阶段，作者构建适当的函数空间（如 $\mathcal{D}, \mathcal{B}$ 类空间，附录 A、B），并在其上定义收缩映射，利用Banach 不动点定理证明解的存在唯一性。

例如，在从弱激波反推至预激波的过程中，未知量包括激波轨迹 $c(t)$ 和内部流场，通过将控制方程重写为积分方程形式，建立映射 $\Phi: X \to X$，并验证其在加权 Hölder 空间中的压缩性。

3. 精细的能量估计与加权范数

由于在接近预激波时出现 $C^{1/3}$ 奇异性，传统 Sobolev 或 Hölder 能量不再适用。作者设计了基于分数阶导数和渐近展开系数的加权范数，例如对 $\mathring{z} = \partial_r z$ 使用 $L^\infty$ 加权：

$$ \|\mathring{z}\|_{\mathcal{X}} = \sup_{r,t} \frac{|\mathring{z}(r,t)|}{|(r - r_*)/(t - T_*)|^{-2/3}}, $$

以捕捉 $|\mathring{z}| \sim |r - r_*|^{-2/3}$ 的发散行为。

4. 特征线方法与渐近匹配

利用三条特征族 $\lambda_1, \lambda_2, \lambda_3$ 构造特征坐标系，在其上分析变量沿特征的传播规律。尤其重要的是识别出导致 $z$ 形成 $C^{1/3}$ 尖点的共振机制：$\lambda_1$ 特征线上非线性项累积效应引发梯度爆炸。

此外，使用中间渐近匹配技术（intermediate asymptotics matching）确保不同区域（弱激波区 vs 自相似区）之间的光滑衔接。

5. 正规激波理论框架

遵循 Christodoulou 等人发展的“正规激波解”框架 [6]，要求：

  • 激波前后状态满足 Rankine-Hugoniot 条件；
  • 满足 Lax 熵条件；
  • 解在激波外光滑，且在 $t = T_*$ 时连续（即从预激波“诞生”）。

这保证了所构造的激波是物理合理的、唯一的。

b) 技术难点与突破技巧

难题1：避免非物理的不对称预激波

若直接将强 Guderley 激波反向减弱至零强度，会导致激波加速度 $s''(t) \to \infty$ 当 $t \searrow T_*$，从而产生非对称的、仅一侧 Hölder 连续的预激波，无法由 $C^1$ 初值得到。

解决技巧：引入两段式激波路径：

  • 第一段 $r = c(t)$：高度正则，用于实现从弱激波到 $C^{1/3}$ 预激波的平滑过渡；
  • 第二段 $r = s(t)$：连接至 Guderley 解。

这样避免了反向演化中加速度奇异性。

难题2：处理 $C^{1/3}$ 奇异性的传播

$z$ 在 $t = T_*$ 形成 $C^{1/3}$ 尖点，意味着 $\partial_r z$ 发散。传统的经典解理论失效。

解决技巧：使用微分黎曼变量 $\mathring{z}, \mathring{w}, \mathring{b}$，它们即使在奇点附近也保持有界或可控增长。通过分析这些变量的输运方程（含非线性源项），控制其在反向演化中的行为。

难题3：实现与自相似解的“全阶匹配”

仅仅函数值匹配不足以保证后续演化一致；需要所有导数在 $t = T_{\mathrm{fin}}$ 匹配。

解决技巧：精心选择中间时间 $T_\circ$ 和连接曲线 $c(t), s(t)$，使得在 $T_{\mathrm{fin}}$ 处不仅函数值相等，而且所有时空导数通过泰勒展开逐阶逼近 Guderley 剖面。这依赖于自相似解的解析性和稳定性。

难题4：非线性项的控制

在 DRV 方程中，存在形如 $\mathring{z}^2, \mathring{z}\mathring{w}, r^{-1}(w^2 - z^2)$ 等非线性项，在 $r \to 0$ 或 $t \to T_*$ 时可能放大误差。

解决技巧：引入渐近展开 ansatz，将解分解为主项 + 余项，主项精确捕获奇异性（如 $(r - r_*)^{1/3}$），余项则在加权空间中进行迭代估计，利用小参数控制非线性反馈。

4. 与之前工作的比较

a) 突破与改进

改进1：首次实现从光滑数据到强 Guderley 激波的严格构造

此前的工作，如 Christodoulou [6]、Luo-Xin [26]、Li-Wang [9] 等，主要集中于弱激波形成（weak shock formation）的研究，即激波强度满足：

$$ [\![u]\!] \sim (t - T_*)^{1/2}, \quad [\![\rho]\!] \sim (t - T_*)^{1/2}, $$

增长缓慢。而本文成功跨越了从弱激波到完全非线性强激波的鸿沟，实现了与真实物理中观测到的 Guderley 解的精确对接。

改进2：突破线性不稳定性障碍，证实“普适吸引子”角色

尽管 Guderley 解已被广泛认为是强内爆的吸引子，但它本身是线性不稳定的（linearly unstable），任何微小扰动都会被放大。因此人们怀疑它是否真的可达。

本文通过构造一条稳定流形上的轨线，证明即使起点远离自相似轮廓，只要初始数据适当配置，仍可精确趋近并锁定 Guderley 解。这为“中间渐近性”提供了首个严格数学实现。

改进3：放宽初始数据正则性要求

以往许多全局解研究要求初始数据足够光滑（如 $C^\infty$ 或解析）。本文允许初始数据仅为 $C^{1,1/3}$，这是形成 $C^{1/3}$ 预激波所需的最低正则性，体现了对临界正则性阈值的精准把握。

改进4：推广至非等熵、多维情形

多数早期自相似解研究集中在等熵或一维情况。本文处理的是全非等熵、二维或三维径向对称的完整欧拉系统，更具物理真实性。

改进5：提供完整的物理演化链条

本文是首个完整展示如下五阶段演化的数学解：

光滑初值 → 预激波形成 → 弱激波诞生 → 向强自相似激波过渡 → 精确匹配 Guderley 解 → 坍缩 → 反射爆轰

这一完整路径填补了数学理论与物理直觉之间的空白。

总结对比

方面以往工作（如 [6,9,10,26]）本文贡献
激波强度仅限弱激波 $(\sim (t-T_*)^{1/2})$实现强激波，精确匹配 Guderley
初始数据多为抽象构造，未必物理合理明确构造具静止核的物理初值
匹配精度局部近似或渐近逼近全阶精确匹配
正则性门槛要求较高光滑性达到临界 $C^{1,1/3}$
全局演化通常止于激波形成延伸至 $t=0$ 并反射继续
数学意义局部存在性验证普适吸引子猜想

综上所述，本文不仅是技术上的重大推进，更是对非线性波动系统中普适性与稳定性辩证关系的一次深刻揭示：即便某个状态本身不稳定，它仍可作为动态演化过程中的“必经之路”和“组织中心”。

2. The local regularity theory for the Stokes and Navier–Stokes equations near the curved boundary

Stokes方程与Navier-Stokes方程在曲边界附近的局部正则性理论

作者: Hui Chen, Su Liang, Tai-Peng Tsai

类别: math.AP

发布日期: 2025-10-14

arXiv ID: 2510.13028v2

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

详细分析

中文标题: Stokes方程与Navier-Stokes方程在曲边界附近的局部正则性理论

1. 研究对象和背景

本文研究的是不可压缩粘性流体在具有光滑弯曲边界的区域中，其运动所满足的Stokes方程和Navier–Stokes方程的局部正则性理论。核心目标是：在非常弱的初始正则性假设下，证明解在靠近曲边界附近具有更高的空间正则性（即导数估计），并基于此提出一个新的“边界正则点”定义。

所研究的偏微分方程系统

设 $\Omega \subset \mathbb{R}^d$ ($d \geq 2$) 是一个有界或无界区域，其边界 $\partial\Omega \in C^{1,1}$。原点 $0 \in \partial\Omega$。记：

  • $B_R = B_R(0)$ 为以原点为中心、半径为 $R$ 的球；
  • $B_R^+ = B_R \cap \Omega$；
  • $Q_R^+ = B_R^+ \times (-R^2, 0)$ 为空间时间柱体。

考虑如下非定常Stokes方程组：

$$ (\mathrm{S}) \quad \begin{cases} \partial_t u - \Delta u + \nabla p = f + \mathrm{div}\, F, \\ \mathrm{div}\, u = 0, \end{cases} $$

其中 $u = (u_1, \dots, u_d)$ 是速度场，$p$ 是压力，$f = (f_1, \dots, f_d)$ 是体积力，$F = (F_{ij})_{1 \le i,j \le d}$ 是表面应力张量。当取 $F = -u \otimes u$ 且 $f=0$ 时，该系统退化为经典的Navier–Stokes方程：

$$ (\mathrm{NS}) \quad \begin{cases} \partial_t u + u \cdot \nabla u - \Delta u + \nabla p = 0, \\ \mathrm{div}\, u = 0. \end{cases} $$

这些方程描述了不可压缩牛顿流体的动力学行为，在数学物理中有深远意义。

边界条件

论文重点研究两种重要的边界条件，作用于边界部分 $\Sigma = \Gamma \times (-1,0)$，其中 $\Gamma = B_1 \cap \partial\Omega$。

  1. 无滑移边界条件（no-slip BC）：

    $$ u = 0 \quad \text{on } \Sigma. $$

  2. Navier滑移边界条件（Navier BC）：

    $$ [(2Du + F)n]_{\tan} + \alpha u_{\tan} = 0, \quad u \cdot n = 0 \quad \text{on } \Sigma, $$

    其中：

    • $n$ 是单位外法向量；
    • $Du = \frac{1}{2}(\nabla u + (\nabla u)^T)$ 是变形张量；
    • $v_{\tan} = v - (v \cdot n)n$ 表示切向分量；
    • $\alpha \geq 0$ 是摩擦系数。

特别地，当 $\alpha = 2\kappa_\tau$（$\kappa_\tau$ 为主曲率张量）时，Navier边界条件可简化为：

$$ \omega n = 0, \quad u \cdot n = 0, $$

其中 $\omega$ 是速度梯度的反对称部分（即涡量矩阵）。在三维情形下，这等价于 $(\mathrm{curl}\, u) \times n = 0$。

函数空间与正则性假设

文章使用混合范数空间 $L^{q,r}(\Omega \times I) = L^r(I; L^q(\Omega))$，其中 $1 < q,r < \infty$。此外引入临界指数：

$$ q_* = \begin{cases} \max\left\{1, \frac{dq}{d+q}\right\}, & q \ne d, \\ 1+, & q = d. \end{cases} $$

这是为了处理压力项对速度的非局部影响而引入的更弱的空间。

核心挑战与动机

对于抛物型方程（如热方程），弱解通常自动光滑。但Stokes方程由于压力项的存在，即使在内部区域也未必光滑；而在边界附近，情况更为复杂。

已有研究表明，在平面边界上：

  • 对于无滑移边界，若不假设压力的可积性，则无法获得高阶导数估计（反例见Kang, Seregin-Sverák等）。
  • 而对于Navier边界，作者前期工作[5]表明：仅需 $u \in L^{q,r}$ 即可推出 $\nabla u \in L^{q,r}$，无需压力假设，这是一个显著优势。

本文将这些结果从平面边界推广到曲边界，技术难度大大增加，因为必须通过坐标变换（boundary straightening）将曲边界拉直，从而引入一系列高阶扰动项，这些项破坏了原有结构。

因此，本文的关键在于：如何在极低初始正则性（例如 $u,p \in L^{q,r}$，甚至仅为非常弱解）的前提下，克服由边界曲率引起的高阶扰动，建立速度场的梯度乃至更高阶导数的估计。

2. 主要定理或主要结果

定理 1.1（无滑移边界下的导数估计）

假设：$1 < q,r < \infty$，$\Gamma \in C^{1,1}$，$(u,p)$ 是Stokes方程(S)在 $Q_1^+$ 上关于无滑移边界条件(1.1)的非常弱解对（见定义2.3），且满足：

$$ u, p, F \in L^{q,r}(Q_1^+), \quad f \in L^{q_*,r}(Q_1^+). $$

结论：则 $(u,p)$ 实际上是弱解对（定义2.4），并且有：

$$ \nabla u \in L^{q,r}(Q_{1/2}^+), $$

且满足梯度估计：

$$ \|\nabla u\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q_*,r}(Q_1^+)} + \|F\|_{L^{q,r}(Q_1^+)}. $$

进一步，若 $f + \mathrm{div}\, F \in L^{q,r}(Q_1^+)$，则：

$$ \partial_t u, \nabla^2 u, \nabla p \in L^{q,r}(Q_{1/2}^+), $$

且满足二阶导数估计：

$$ \|\partial_t u, \nabla^2 u, \nabla p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f + \mathrm{div}\, F\|_{L^{q,r}(Q_1^+)}. $$

重要性：首次在曲边界上证明非常弱解能提升为强解，且估计中不需要 $\|\nabla u\|$ 在右端，优于以往工作。

定理 1.2（Navier边界下的边界梯度估计）

假设：同上 $q,r$ 条件，$\Gamma \in C^{1,1}$，$u,F \in L^{q,r}(Q_1^+)$，$f \in L^{q_*,r}(Q_1^+)$，$\nabla u \in L^1(Q_1^+)$，且 $u$ 是Navier边界条件下Stokes方程的弱解（定义2.5）。

结论：则

$$ \nabla u \in L^{q,r}(Q_{1/2}^+), $$

且满足：

$$ \|\nabla u\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q_*,r}(Q_1^+)} + \|F\|_{L^{q,r}(Q_1^+)}. $$

改进：将[5]中平面边界的结果推广至任意维数的曲边界，且右端无需压力项。

定理 1.3（Navier边界下的高阶导数估计）

假设：$\Gamma \in C^{2,1}$，$u,p,f \in L^{q,r}(Q_1^+)$，$F = 0$，$\nabla u \in L^1(Q_1^+)$，$(u,p)$ 是Navier边界条件下的弱解对（定义2.6）。

结论：则

$$ \partial_t u, \nabla^2 u, \nabla p \in L^{q,r}(Q_{1/2}^+), $$

且满足：

$$ \|\partial_t u, \nabla^2 u, \nabla p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q,r}(Q_1^+)}. $$

进一步，若 $\Gamma \in C^{3,1}$ 且 $\nabla f \in L^{q,r}(Q_1^+)$，则：

$$ \nabla \partial_t u, \nabla^3 u, \nabla^2 p \in L^{q,r}(Q_{1/2}^+), $$

且满足三阶导数估计：

$$ \|\nabla \partial_t u, \nabla^3 u, \nabla^2 p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^r(-1,0; W^{1,q}(B_1^+))}. $$

突破：首次在曲边界上建立Navier条件下三阶导数估计，且无需先验假设高阶导数存在。

应用：新定义的边界正则点

作者提出新的边界正则点定义（Definition 1.4）：

称点 $z_0 = (x_0,t_0) \in \partial\Omega \times (0,T]$ 是Navier–Stokes方程的一个边界正则点，如果存在邻域 $Q_R^+(z_0)$ 使得：

$$ u \in L^\infty(Q_R^+(z_0)), \quad p \in L^{q_0,r_0}(Q_R^+(z_0)), \quad \text{for some } 1 < q_0,r_0 < \infty. $$

关键改进：相比经典定义（仅要求 $u$ 有界），此处显式要求压力 $p$ 属于某个 $L^{q_0,r_0}$ 空间，从而排除已知反例，确保更高正则性。

由此得到两个应用定理：

  • 定理 1.5（NS, 无滑移）：若 $z_0$ 是上述意义下的边界正则点，则存在 $\rho > 0$ 使得 $\partial_t u, \nabla^2 u, \nabla p \in L^{q_1,r_0}(Q_\rho^+(z_0))$ 对任意 $q_1 < \infty$ 成立。
  • 定理 1.6（NS, Navier）：若 $\partial\Omega \in C^{3,1}$，且 $z_0$ 是边界正则点，则更高阶导数 $\nabla \partial_t u, \nabla^3 u, \nabla^2 p \in L^{q_1,r_0}(Q_\rho^+(z_0))$。

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

a) 核心分析工具

  1. 边界拉直技术（Boundary Straightening）
    将曲边界局部映射为平面。设 $\partial\Omega$ 在原点附近由 $x_d = \gamma(x')$ 给出，$\gamma \in C^{1,1}$，$\gamma(0)=0$, $\nabla\gamma(0)=0$。定义变换：

    $$ T x = (x', x_d - \gamma(x')), \quad T^{-1}x = (x', x_d + \gamma(x')). $$

    则 $T$ 将 $\Omega$ 映射到上半空间 $\mathbb{R}_+^d$。此变换引入了高阶扰动项，如：

    $$ T^{-1} \Delta = \Delta - 2\nabla\gamma \cdot \nabla \partial_d + |\nabla\gamma|^2 \partial_d^2 - \Delta\gamma \cdot \partial_d T^{-1}, $$

    这些项在后续估计中必须精细控制。

  2. 磨光逼近（Mollification with Cut-off and Extension）
    构造速度场的逼近序列。定义截断函数 $\zeta(x,t)$ 和磨光核 $\eta^\pm_{\varepsilon_1,\varepsilon_2,\varepsilon_3}$，采用奇偶延拓结合水平/垂直/时间方向分别磨光：

    $$ E^\pm_{\varepsilon}(g)(x,t) = \int g(y,s) \eta^\pm(x,y,t-s) dy ds, $$

    特别地，$E^+$ 对应偶延拓，$E^-$ 对应奇延拓，用于保持边界条件。

  3. 正规形构造（Normal Form）
    引入中间变量 $v_k = T^{-1} u_k \cdot \zeta$，$\pi = T^{-1} p \cdot \zeta$，称为“正规形”。这些量在拉直后空间中满足修正的Stokes方程，包含由边界曲率产生的高阶扰动项。

  4. 能量型估计与吸收技巧
    通过对磨光后的方程乘以测试函数并积分，利用Young不等式、Gronwall引理等，将误差项吸收至左侧。例如，通过选择足够小的 $\rho$ 使 $\|\nabla\gamma\|_{L^\infty} \ll 1$，从而控制扰动项。

  5. 嵌入定理与迭代论证
    利用Sobolev嵌入 $W^{1,q} \hookrightarrow L^{q^*}$ 及其时间版本，结合迭代方法提升正则性。例如在Theorem 1.5中，由 $u \in L^\infty$ 推出 $\nabla u \in L^{q_0,r_0}$，再得 $p \in L^{q_1,r_0}$，循环迭代直至 $q_1 \to \infty$。

b) 关键技术难点与解决技巧

  1. 恢复垂直导数估计（Recovering Vertical Derivatives）
    在边界附近，直接估计 $\partial_d u$ 困难。作者利用方程约束（如 $\Delta u = 0$ 或 $\mathrm{div}\, u = 0$）将 $\partial_d^2 u$ 表达为 $\nabla' \nabla u$ 的组合，从而通过水平导数估计控制垂直导数。例如：

    $$ \partial_d^2 u_d = -\sum_{i=1}^{d-1} \partial_i \partial_d u_i, $$

    故 $\|\partial_d^2 u\| \sim \|\nabla' \nabla u\|$。

  2. 时间导数转为空间导数（Transferring Temporal to Spatial Derivatives）
    方程中出现 $\partial_t u$ 项难以处理。作者利用Stokes方程本身：

    $$ \partial_t u = \Delta u - \nabla p + f + \mathrm{div}\, F, $$

    将时间导数转化为空间二阶导数加压力梯度，从而统一估计。

  3. 处理高阶扰动项的“正常形”技巧
    拉直后方程右侧出现形如 $\nabla\gamma \cdot \nabla \partial_d v$ 的高阶项。作者通过设计特殊测试函数（满足边界条件），并利用交换子估计（commutator estimate）控制此类项。例如引理2.2(v)给出：

    $$ \|\partial_{x_k} E_h(ag) - a \partial_{x_k} E_h(g)\| \lesssim \|\nabla a\|_{L^\infty} \|g\|_{L^{q,r}}, $$

    用于处理系数与导数的非交换性。

  4. 吸收法控制小参数项
    所有扰动项均含 $\nabla\gamma$ 或 $\Delta\gamma$，其模长由 $C^{1,1}$ 正则性控制为 $O(\rho)$。通过选取足够小的邻域半径 $\rho$，使得这些项前的系数小于1/2，从而可被主部吸收。

  5. 避免压力项主导的策略
    在无滑移情况下，压力可能引起非局部效应。作者通过在估计中显式保留 $\|p\|$ 项，排除已知反例（如Kang构造的有界速度但无界压力情形），确保估计有效性。

4. 与之前工作的比较

a) 主要突破与改进

  1. 从平面边界到曲边界的推广
    此前关于Stokes方程边界正则性的研究（如Seregin-Shilkin [31], Chang-Kang [2], 作者前期工作[5]）大多局限于平面边界。本文首次在一般 $C^{1,1}$ 曲边界上建立完整的局部正则性理论，技术上需应对坐标变换带来的复杂扰动。

  2. 起始正则性假设大幅降低
    多数先验估计要求解已具备一定正则性（如 $\nabla u \in L^{q,r}$）。而本文从非常弱解出发（仅 $u,p \in L^{q,r}$），通过逼近法严格推导出高阶可微性，更具基础性。

  3. Navier边界下的最优高阶估计
    在Navier边界下，作者不仅推广了[5]的一阶估计，还建立了二阶乃至三阶导数估计（Theorem 1.3），且无需 $F \ne 0$（除非 $[Fn]_{\tan}|_\Sigma = 0$）。这比Dong-Kwon [9]的结果更强，因其不要求高阶导数的先验有界性。

  4. 提出新的边界正则点定义
    经典定义（如Scheffer, Seregin）仅要求速度有界，但在边界上不足以保证光滑性（反例存在）。本文提出的定义加入压力可积性条件，从而确保解在边界附近真正光滑，填补了部分正则性理论的空白。

  5. 维度推广至任意 $d \geq 2$
    许多早期工作集中于三维情形。本文结果适用于所有 $d \geq 2$，增强了普适性。

b) 结果强度对比

工作边界类型起始正则性是否需要压力最高导数阶
Seregin-Shilkin [31]平面$u,\nabla u,p \in L^{q,r}$是2
Chang-Kang [2]平面$u \in L^{q,r}$是2
本文 Theorem 1.1曲面$u,p \in L^{q,r}$（非常弱解）是2
作者 [5]平面$u \in L^{q,r}$（非常弱解）否1–2
本文 Theorem 1.2曲面$u \in L^{q,r}, \nabla u \in L^1$否1
本文 Theorem 1.3曲面$u,p \in L^{q,r}$是3

可见，本文在边界几何、起始正则性、结果阶数三个方面均实现显著推进。

c) 技术创新点总结

  • 首次系统处理边界拉直引入的高阶扰动项；
  • 发展了适用于曲边界的正规形+磨光+吸收法框架；
  • 提出“时间导数转空间导数”策略以统一估计；
  • 建立了从边界正则点定义到高阶正则性的完整链条。

综上，本文在Stokes与Navier–Stokes方程的边界正则性理论中做出了奠基性贡献，为后续研究曲边界上的流体行为提供了强有力的工具和新的视角。

3. Intermittent Solutions of the Stationary 2D Surface Quasi-Geostrophic Equation

二维定常表面准地转方程的间歇性解

作者: Nicholas Gismondi, Alexandru F. Radu

类别: math.AP

发布日期: 2025-10-18

arXiv ID: 2510.16583v2

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

详细分析

中文标题: 二维定常表面准地转方程的间歇性解

1. 研究对象和背景

本文研究的是二维定常表面准地转方程（Stationary Surface Quasi-Geostrophic Equation, 简称 SQG），其形式为：

$$ \begin{cases} \mathrm{div}(\theta u) + \Lambda^\gamma \theta = 0 \\ u = \Lambda^{-1} \nabla^\perp \theta \end{cases} \quad \text{在 } \mathbb{T}^2 \text{ 上}, $$

其中 $\theta : \mathbb{T}^2 \to \mathbb{R}$ 是一个均值为零的标量函数，代表表面势温；$u : \mathbb{T}^2 \to \mathbb{R}^2$ 是由 $\theta$ 通过非局部关系导出的速度场；$\nabla^\perp = (-\partial_2, \partial_1)$ 表示旋转梯度；$\Lambda = (-\Delta)^{1/2}$ 是分数拉普拉斯算子，其傅里叶乘子为 $|\xi|$。参数 $\gamma \in (0,2]$ 控制耗散强度。

该方程是流体力学中重要的准地转模型，用于描述大气或海洋表层在强旋转与分层条件下的动力学行为。从数学角度看，SQG 方程因其非线性结构与三维欧拉方程中涡量演化的相似性而受到广泛关注。特别地，当 $\gamma = 0$ 时称为无粘 SQG 方程，对应于保守系统；当 $\gamma > 0$ 时引入了耗散项 $\Lambda^\gamma \theta$，模拟小尺度能量耗散。

本文关注的是定常情形，即时间导数项为零的情况，这与大多数以往研究集中在非定常问题（如 $\partial_t \theta + \mathrm{div}(\theta u) + \Lambda^\gamma \theta = 0$）形成对比。定常解的存在性本身具有重要意义：它们可能作为长期动力行为的候选态，或揭示方程本身的结构性缺陷（如弱解的非唯一性）。

论文考虑的空间区域是二维环面 $\mathbb{T}^2 = [0,1]^2$，所有函数均为周期函数且均值为零。这一设定避免了边界效应，便于使用傅里叶分析工具。作者特别强调，他们构造的解属于 负索博列夫空间 $\dot{H}^{-1-\epsilon}(\mathbb{T}^2)$ 对任意 $\epsilon > 0$ 成立，但不落在 $\dot{H}^{-1/2}(\mathbb{T}^2)$ 中。

这一点极为关键。标准的弱解理论通常要求 $\theta \in \dot{H}^{-1/2}$，因为在此空间中可以良好定义非线性项 $\mathrm{div}(\theta u)$ 的对偶配对。例如，在 Cheng-Kwon-Li [11] 的工作中，正是基于此正则性假设来定义弱解：

$$ \frac{1}{2} \int_{\mathbb{T}^2} (\Lambda^{-1/2}\theta) \Lambda^{1/2}[R_\perp, \nabla \psi]\theta \, dx = - \int_{\mathbb{T}^2} (\Lambda^{-1/2}\theta) \Lambda^{\gamma+1/2}\psi \, dx, $$

其中 $[R_\perp, \nabla \psi]$ 是 Riesz 变换与梯度的交换子。然而，当 $\theta$ 的正则性进一步下降至 $\dot{H}^{-1-\epsilon}$ 时，上述积分不再有意义——乘积 $\theta u$ 失去可积性，传统的弱解框架失效。

因此，本文的核心挑战之一是如何重新定义“解”的概念，以容纳更低正则性的对象。为此，作者引入了抛物积（paraproduct）方法，借助 Littlewood-Paley 分解将非线性项分解为高频-低频相互作用，并利用双线性估计控制其分布意义下的行为。

更进一步，作者指出传统凸积分方法（convex integration）在处理此类方程时面临所谓的“奇乘子障碍”（odd multiplier obstruction）。这是因为速度 $u = \Lambda^{-1} \nabla^\perp \theta$ 与 $\theta$ 之间的关系由一个奇函数傅里叶乘子 $i\xi^\perp / |\xi|$ 给出，导致高频率振荡项在非线性相互作用中完美抵消，无法生成所需的低频误差项来进行误差消除。

为克服这一障碍，作者采用 Buckmaster-Shkoller-Vicol [8] 提出的松弛动量方程（relaxed momentum formulation）：

$$ \begin{cases} u \cdot \nabla v - (\nabla v)^T \cdot u + \Lambda^\gamma v + \nabla p = \mathrm{div}(R) \\ \mathrm{div}(v) = 0 \\ v = \Lambda^{-1} u \end{cases} $$

在这个新框架下，未知量变为 $(u,v,p,R)$，其中 $R$ 是雷诺应力张量，代表近似解与真实解之间的偏差。通过将原方程重写为此形式，可以绕过奇乘子问题，因为现在构造的关键对象 $v$ 与 $u$ 的关系为 $u = \Lambda v$，其乘子 $|\xi|$ 是偶函数，更适合进行高频振荡设计。

最终目标是构造一系列近似解 $(u_q, v_q, p_q, R_q)$，使得雷诺应力 $R_q \to 0$ 在适当拓扑下收敛，从而极限 $(u,v)$ 满足原方程。整个构造过程基于凸积分法，这是一种用于证明偏微分方程弱解非唯一性的强大技术，近年来在流体方程（如 Euler、Navier-Stokes）中取得突破性进展。

值得注意的是，由于使用了间歇性 Mikado 流作为构建块，作者能够确保所构造的解满足额外的可积性性质：$\Lambda^{-1}\theta \in L^p(\mathbb{T}^2)$ 对所有 $1 \leq p < 2$ 成立。这是非间歇性构造（如平面波或 Beltrami 流）难以实现的，后者往往只能保证 $\Lambda^{-1}\theta \in \dot{H}^{-\epsilon}$。

综上所述，本文的研究对象是一个物理上有意义、数学上极具挑战性的定常非线性非局部偏微分方程系统。其特殊之处在于：

  • 解的空间正则性极低（$\dot{H}^{-1-\epsilon}$）；
  • 非线性项需用抛物积重新解释；
  • 使用松弛动量形式规避奇乘子障碍；
  • 构造依赖高度振荡且具间歇性的 Mikado 流；
  • 目标是证明非平凡弱解的存在性，而非经典解或适定性结果。

这些特征共同构成了本文工作的创新基础。

2. 主要定理或主要结果

本文的主要定理如下：

定理 1.3（二维 SQG 方程的非平凡定常解）
对于任意给定的耗散指数 $0 < \gamma \leq 2$，存在非零函数 $u, \theta \in \dot{H}^{-1-\epsilon}(\mathbb{T}^2)$ 对所有 $\epsilon > 0$ 成立，使得：

(a) $u = \Lambda^{-1} \nabla^\perp \theta$；

(b) 抛物积 $P_{\neq 0}(\theta u)$ 在 $\dot{H}^{-5}(\mathbb{T}^2)$ 中良定义；

(c) $(u, \theta)$ 构成 SQG 方程的弱抛物积解（weak paraproduct solution），即对所有光滑测试函数 $\phi$，有

$$ \langle \theta, \Lambda^\gamma \phi \rangle_{\dot{H}^s, \dot{H}^{-s}} - \left\langle P_{\neq 0}(\theta u), \nabla \phi \right\rangle_{\dot{H}^{s'}, \dot{H}^{-s'}} = 0, $$

并且 $\Lambda^{-1}\theta \in L^p(\mathbb{T}^2)$ 对所有 $1 \leq p < 2$ 成立。

此外，作者还提出了一个归纳命题（Proposition 4.1），它是构造解序列的核心工具：

命题 4.1（归纳命题）
设 $(u_q, v_q, p_q, R_q)$ 是光滑、均值为零、满足松弛动量方程的近似解序列，且满足以下归纳假设：

  1. $u_q, v_q$ 均值为零且散度自由；
  2. $(u_q, v_q, p_q, R_q)$ 满足
    $$ u_q \cdot \nabla v_q - (\nabla v_q)^T \cdot u_q + \Lambda^\gamma v_q + \nabla p_q = \mathrm{div}(R_q); $$
  3. $\|R_q\|_{\dot{H}^{-4}} < 2^{-q}$；
  4. 对所有 $q' \leq q$，有 $\|v_{q'} - v_{q'-1}\|_{L^p} \lesssim 2^{-q'}$，$1 \leq p < 2$；
  5. 存在与 $q$ 无关的常数 $\delta > 0$，使得 $\|v_q\|_{L^1} > (1 + 2^{-q})\delta$；
  6. 存在唯一的 $j$ 使得 $P_{2j}(v_q - v_{q-1}) = v_q - v_{q-1}$；
  7. 存在常数 $C_1, C_2 > 0$ 独立于 $q$，使得
    $$ \sum_{\substack{n,m \leq q \\ n \ne m}} \|(\nabla^\perp \cdot (v_n - v_{n-1}))(u_m - u_{m-1})\|_{\dot{H}^{-5}} < C_1 - 2^{-q}, $$
    $$ \sum_{n \leq q} \|(\nabla^\perp \cdot (v_n - v_{n-1}))(u_n - u_{n-1})\|_{\dot{H}^{-5}} < C_2 - 2^{-q+100}. $$

若上述归纳假设成立，则可通过迭代构造收敛到非平凡弱解。

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

a) 核心分析学工具

本文采用了现代偏微分方程中前沿的构造性方法，主要包括以下几个关键技术：

首先是 Littlewood-Paley 理论。作者使用标准的频率投影算子 $P_{2^j}$ 将函数分解为不同频带上的成分。这种分解使得他们可以在频域中精确控制各项的正则性和相互作用。特别是，通过定义 $P_{\neq 0}(fg)$ 为抛物积：

$$ P_{\neq 0}(fg) := \sum_{j,j' \geq 0} P_{\neq 0}(P_{2^j}f \, P_{2^{j'}}g), $$

并要求其在某个 $\dot{H}^s$ 空间中绝对可和，从而赋予低正则性函数乘积以分布意义。这种方法源于调和分析中的 Bony 抛物积分解思想。

其次是 对偶配对与负索博列夫范数。由于解位于 $\dot{H}^{-1-\epsilon}$，传统的 $L^2$ 内积不再适用。作者使用对偶配对 $\langle \cdot, \cdot \rangle_{\dot{H}^s, \dot{H}^{-s}}$ 来定义弱解，并利用如下等价范数：

$$ \|f\|_{\dot{H}^{-s}} = \sup_{\|\phi\|_{\dot{H}^s} = 1} |\langle f, \phi \rangle|. $$

这允许他们在分布意义下处理方程的各项。

第三是 凸积分法（Convex Integration）。这是一种源自 Gromov 几何理论的方法，后被 De Lellis 和 Székelyhidi 引入流体力学，用于构造反例（如 Onsager 猜想的违反者）。其基本思想是构造一个快速振荡的增量 $w_{q+1}$，使其非线性部分能精确抵消前一步的雷诺应力 $R_q$。每一步都保持方程的近似满足，同时让误差趋于零。

第四是 Mikado 流（Mikado Flows）作为构建块。这类函数最早由 Daneri 和 Székelyhidi 引入，具有如下特性：

  • 支集集中在细管状区域（体现间歇性）；
  • 散度自由；
  • 其张量积的平均等于预设方向上的外积；
  • $L^1$ 范数可任意小，而 $L^2$ 范数固定。

这些性质使得 Mikado 流非常适合用于构造低正则性解，尤其是在需要精细控制非线性项时。

第五是 逆散度算子 $R$ 的使用。给定一个均值为零的向量场 $u$，定义对称矩阵 $Ru$ 满足 $\mathrm{div}(Ru) = u$。这在构造雷诺应力时至关重要，因为它允许将误差项“提升”为一个张量，从而匹配方程结构。

最后是 Kato-Ponce 型估计。文中使用的引理 2.14：

$$ \|\alpha \beta\|_{\dot{H}^{-s}} \lesssim \|\alpha\|_{C^s} \|\beta\|_{\dot{H}^{-s}}, $$

表明当一个因子足够光滑时，乘积的负索博列夫范数可由另一个因子控制。这在处理高频-低频乘积时极为有用。

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

本文最关键的突破在于成功克服了“奇乘子障碍”。如前所述，原始 SQG 中 $u = \Lambda^{-1}\nabla^\perp \theta$ 的乘子是奇函数，导致高频振荡在非线性项中完全抵消，无法生成有效误差。作者通过引入松弛动量方程彻底规避了这个问题。

另一个核心技术是间歇性的引入。传统的凸积分使用平面波或 Beltrami 波，其能量均匀分布在整个空间。而间歇性 Mikado 流的能量集中在稀疏区域内，使得 $\|W_k\|_{L^1} \ll 1$，即使 $\|W_k\|_{L^2} \sim 1$。这一机制使得作者能够控制 $\Lambda^{-1}\theta$ 的 $L^p$ 可积性，达到 $p < 2$，这是非间歇性方法无法做到的。

此外，作者巧妙利用了几何引理（Lemma 2.10）：通过对称张量空间中的单位球邻域 $B(I,\epsilon)$，找到有限个方向 $k \in \Omega \subset S^1$ 和光滑函数 $\Gamma_k(R)$，使得

$$ R = \frac{1}{2} \sum_{k \in \Omega} (\Gamma_k(R))^2 k^\perp \otimes k^\perp. $$

这使得他们可以将任意小的雷诺应力 $R_q$ 分解为若干 Mikado 流的方向叠加，从而构造出精确匹配的修正项。

最后，为了处理抛物积的收敛性，作者在归纳假设中加入了交叉项求和的控制（假设 7），确保所有历史增量之间的非线性耦合不会累积破坏整体收敛性。

4. 与之前工作的比较

a) 突破与改进

本文相较于先前工作有多项显著改进：

首先，放宽了解的正则性要求。此前 Cheng-Kwon-Li [11] 的结果要求 $\theta \in \dot{H}^{-1/2}$，这是一个临界空间。本文首次构造出仅属于 $\dot{H}^{-1-\epsilon}$ 的非平凡解，突破了这一门槛。这意味着弱解理论可以延拓到更广的函数类中。

其次，适用于任意耗散指数 $\gamma \in (0,2]$。已有工作往往对 $\gamma$ 有严格限制，如 [11] 要求 $\gamma < 3/2$。本文由于使用了间歇性构造，避免了对 $\gamma$ 的谱条件依赖，实现了完整区间覆盖。

第三，解决了奇乘子障碍的根本问题。虽然 Buckmaster 等人 [8] 已提出松弛动量方法，但他们处理的是非定常情形，拥有时间自由度。本文将其成功应用于定常情形，没有时间演化辅助，难度更高。这是该方法首次在纯空间设置中奏效。

第四，获得更强的可积性结果：$\Lambda^{-1}\theta \in L^p$ 对所有 $p < 2$ 成立。相比之下，非间歇性构造通常只能得到 $\Lambda^{-1}\theta \in \dot{H}^{-\epsilon}$，连 $L^1$ 都难保证。这一改进不仅增强了结果的实用性，也为未来研究能量守恒等问题提供了可能路径。

最后，本文展示了间歇性在凸积分中的决定性作用。过去虽有多个工作使用间歇性（如 [1], [3], [9]），但本文明确指出它是解除 $\gamma$ 限制的关键机制，并系统论证了其在控制 $L^p$ 范数方面的优势。

综上，本文不仅扩展了 SQG 方程弱解的存在范围，而且推动了凸积分技术的发展，特别是在处理非局部、奇乘子、定常系统方面提供了新范式。

4. Discretely self-similar solutions for electron MHD

电子磁流体动力学中的离散自相似解

作者: Nada Adzic Vukotic, Mimi Dai

类别: math.AP, 35Q35, 76B03, 76D09, 76E25, 76W05

发布日期: 2025-10-21

arXiv ID: 2510.19000v1

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

详细分析

中文标题: 电子磁流体动力学中的离散自相似解

研究对象和背景

本文研究的是无电阻的电子磁流体动力学（electron magnetohydrodynamics, 简称 electron MHD）方程中的一类特殊奇性解——离散自相似解（discretely self-similar solutions），并旨在证明在多种条件下，这类非平凡的奇性解不存在。

所研究的主控方程为：

$$ \partial_t B + \nabla \times ((\nabla \times B) \times B) = 0, $$

$$ \nabla \cdot B = 0. $$

其中 $ B(x,t) $ 是三维空间中的磁场向量场，满足磁高斯定律 $\nabla \cdot B = 0$。这个系统是从完整的非电阻MHD系统中引入霍尔效应（Hall effect）后，在忽略离子运动、仅考虑电子运动的理想化情形下导出的模型（参见文献[1,3]）。该模型在等离子体物理中尤为重要，因为它能够捕捉由霍尔项驱动的快速磁重联现象（rapid magnetic reconnection），这是太阳耀斑、地球磁层活动等天体物理过程中的关键机制。

与经典的欧拉方程或标准MHD不同，electron MHD中的非线性项 $\nabla \times ((\nabla \times B) \times B)$ 被称为霍尔项（Hall term），其结构高度奇异且具有更强的非局部性。这一项使得系统的数学分析极具挑战性。特别地，由于方程中没有耗散项（如拉普拉斯项 $\Delta B$），因此无法通过粘性正则化来控制解的行为，这导致人们普遍预期：若解发生爆破（blow-up），其奇性可能非常剧烈。

为了理解这种潜在的奇性行为，作者聚焦于一类具有特定对称结构的候选奇性解——自相似解与离散自相似解。这些解假设在接近奇点 $(x,t)=(0,0)$ 时，解呈现出某种尺度不变的形式。

具体来说，若 $ B(x,t) $ 是方程的一个解，则通过如下尺度变换：

$$ B_\lambda(x,t) = \lambda^\alpha B(\lambda x, \lambda^{\alpha+2} t), \quad \forall \lambda > 0, $$

得到的新函数 $ B_\lambda $ 也是原方程的解（但初始数据相应缩放）。当该关系对所有 $\lambda > 1$ 成立时，称 $ B $ 为自相似解；而若仅对某个固定的 $\lambda > 1$ 成立，则称为离散自相似解。后者意味着解在时间对数尺度下具有周期性结构。

引入自相似变量：

$$ y = \frac{x}{(-t)^{1/(\alpha+2)}}, \quad s = -\log(-t), $$

并将原解表示为：

$$ B(x,t) = \frac{1}{(-t)^{\alpha/(\alpha+2)}} H\left( \frac{x}{(-t)^{1/(\alpha+2)}}, -\log(-t) \right), $$

其中 $\alpha \ne -2$。代入原方程可得剖面函数 $ H(y,s) $ 满足如下系统：

$$ \partial_s H + \frac{\alpha}{\alpha+2} H + \frac{1}{\alpha+2}(y \cdot \nabla) H + \nabla \times ((\nabla \times H) \times H) = 0, $$

$$ \nabla \cdot H = 0. $$

此时，$ B $ 是自相似解当且仅当 $ H $ 与时间变量 $ s $ 无关；而 $ B $ 是离散自相似解当且仅当 $ H(y,s) $ 在 $ s $ 上是周期函数，周期为 $ S_0 = (\alpha+2)\log \lambda $。

因此，寻找非平凡的离散自相似爆破解，等价于寻找上述剖面方程的非零周期解 $ H(y,s) $。本文的目标正是在各种空间衰减或非衰减条件下，排除此类非平凡周期解的存在性，从而证明：在这些条件下，不存在非平凡的离散自相似爆破解。

研究区域为全空间 $ \mathbb{R}^3 \times [0,S_0] $，边界条件隐含于函数的空间行为中（如无穷远处趋于零或有界），并未施加显式的边值条件。函数空间方面，文中多次使用 $ L^q $、$ W^{1,p} $、$ C^k $ 类空间，并强调了解的光滑性假设（如 $ C^2C^1 $ 或 $ C^\infty $），以及在时空柱状域上的可积性条件（如 $ H \in L^q(\mathbb{R}^3 \times [0,S_0]) $）。这些空间设定反映了对解的整体行为与局部正则性的双重关注。

综上所述，本文的研究对象是一个描述无电阻电子MHD中磁场演化的拟线性偏微分方程组，重点在于分析其是否存在具有离散尺度不变性的爆破解，并通过严格的能量估计与渐近分析方法，在多种函数空间和衰减条件下证明了这类解的非存在性。

主要定理或主要结果

本文共陈述并证明了六个主要定理，分别针对不同的尺度参数 $\alpha$ 和剖面函数 $ H $ 与电流密度剖面 $ J = \nabla \times H $ 的空间衰减或增长条件，给出了排除非平凡离散自相似解存在的充分条件。

记 $ J = \nabla \times H $ 为电流密度剖面，$ S_0 > 0 $ 为周期长度。

定理 1.1

假设条件：

  • $\alpha \ne -2$；
  • $ H \in C^2C^1(\mathbb{R}^3 \times \mathbb{R}) $（即空间二阶连续可微，时间一阶连续可微）；
  • $ J(y) = o(|y|) $ 当 $ |y| \to \infty $；
  • 存在 $ r > 0 $，使得对所有 $ 0 < q < r $，有 $ H \in L^q(\mathbb{R}^3 \times [0,S_0]) $。

结论：
$ H \equiv 0 $，即唯一可能的解是平凡解。

此定理表明，在电流密度增长慢于线性、且剖面函数在低阶 $ L^q $ 空间可积的前提下，不存在非平凡解。

定理 1.2

假设条件：

  • $ p > \frac{6}{5} $；
  • 分两种情况：
    • (i) $\alpha > \frac{3}{2}$ 且 $ H \in L^3_s W^{1,p}_y $；
    • (ii) $-2 < \alpha < \frac{3}{2}$ 且 $ H \in L^2_s L^2_y \cap L^3_s W^{1,p}_y $。

结论：
$ H \equiv 0 $。

该定理利用了剖面函数梯度的 $ L^p $ 可积性，结合Sobolev嵌入，适用于不同尺度范围下的正则性条件。

定理 1.3

假设条件：

  • $ H \in C^2C^1(\mathbb{R}^3 \times \mathbb{R}) $，周期为 $ S_0 > 0 $；
  • $\limsup_{|y|\to\infty} (|J(y,s)| + |\nabla J(y,s)|) = 0$ 对所有 $ s \in (0,S_0) $ 成立；
  • 进一步假设：
    • (i) 若 $\alpha > 0$，存在 $ 0 < q < \frac{3}{\alpha} $；
    • (ii) 若 $-2 < \alpha \leq 0$，存在 $ q > 0 $， 使得 $ H \in L^q(\mathbb{R}^3 \times [0,S_0]) $。

结论：
$ H \equiv 0 $。

此定理要求电流及其梯度在无穷远处一致趋于零，结合适当的 $ L^q $ 可积性即可排除非平凡解。

定理 1.4

假设条件：

  • $ H \in C^2C^1(\mathbb{R}^3 \times \mathbb{R}) $，周期为 $ S_0 > 0 $，$\alpha < -2$；
  • $ H(y,s) = o(1) $，$\nabla J(y,s) = o(1)$ 当 $ |y| \to \infty $，对所有 $ s \in [0,S_0] $。

结论：
$ H \equiv 0 $。

注意此处 $\alpha < -2$ 意味着时间尺度指数为负，对应“超快”压缩情形。尽管如此，只要剖面函数和电流梯度趋于零，仍可排除非平凡解。

定理 1.5

假设条件：

  • $ H \in C^3C^1(\mathbb{R}^3 \times \mathbb{R}) $，周期为 $ S_0 > 0 $，$\alpha > 0$；
  • 原点 $ y=0 $ 是每个分量 $ J_i(\cdot, s) $ 的局部极值点；
  • $ J(0,s) = 0 $ 对所有 $ s \in (0,S_0) $。

结论：
$ H \equiv 0 $。

这是一个基于局部结构分析的结果，不依赖无穷远衰减，而是利用原点处的临界点性质推导矛盾。

定理 1.6

假设条件：

  • $\alpha > \frac{3}{2}$；
  • $ H \in C^\infty(\mathbb{R}^3 \times \mathbb{R}) $，周期为 $ S_0 > 0 $；
  • $\sup_{s \in \mathbb{R}} |H(y,s)| = o(|y|^2)$ 当 $ |y| \to \infty $。

结论：
$ H \equiv 0 $。

值得注意的是，这是一个非衰减型条件：允许 $ H $ 随 $ |y| $ 增长（最多次二次），但仍能排除非平凡解。这说明即使剖面函数不趋于零，只要增长速度低于某阈值，也无法支持离散自相似爆破。

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

本文的核心技术路线是能量估计法结合截断函数技巧（cut-off function method）与极限过程分析，通过对剖面方程进行测试函数乘积积分，建立关于 $ H $ 的积分恒等式或不等式，再通过令截断半径 $ R \to \infty $ 来提取全局信息，最终推出 $ H \equiv 0 $。

a) 核心分析工具

  1. 能量估计（Energy Estimates）
    这是全文最核心的工具。作者将剖面方程（如式 (2.1)）与测试函数（如 $ |H|^{q-2}H \phi_R $ 或 $ H \phi_R $）作内积后在时空区域上积分，利用分部积分处理对流项与非线性项，从而获得含有 $ |H|^q $ 的积分表达式。例如在定理1.1中，通过选取 $ |H|^{q-2}H \phi_R $ 作为测试函数，得到：

    $$ \frac{3b}{q} \int_0^{S_0} \int |H|^q \phi_R \,dyds + \text{nonlinear term} = \text{boundary terms on annulus}. $$

    利用周期性消去时间导数项，再估计右侧边界项随 $ R \to \infty $ 趋于零，左侧若系数发散（如 $ q \to 0 $ 时 $ |a - 3b/q| \to \infty $），即可迫使 $ \|H\|_{L^q} = 0 $。

  2. 截断函数方法（Cut-off Function Method）
    定义光滑径向截断函数 $ \phi_R(y) = \phi(|y|/R) $，满足：

    $$ \phi(r) = \begin{cases} 1, & r < 1 \\ 0, & r > 2 \end{cases}, \quad 0 \leq \phi \leq 1. $$

    其梯度 $ \nabla \phi_R $ 支集在环形区域 $ R \leq |y| \leq 2R $，且满足 $ |\nabla \phi_R| \lesssim 1/R $。这一构造使得所有边界项集中在薄环上，便于利用 $ H,J $ 的渐近行为控制其大小。

  3. Hölder 不等式与 Sobolev 嵌入
    在定理1.2的证明中，面对非线性项 $ ((\nabla \times H) \times H) \cdot (H \times \nabla \phi_R) $，作者应用 Hölder 不等式：

    $$ \int | \nabla H ||H|^2 |\nabla \phi_R| dy \lesssim \|\nabla H\|_{L^p} \|H\|_{L^{2p_1}}^2 \cdot R^{-1} \cdot R^{3(1 - 1/p_2)}, $$

    并借助 Sobolev 嵌入 $ W^{1,p} \hookrightarrow L^{p^*} $（其中 $ p^* = 3p/(3-p) $）来控制 $ H $ 的 $ L^{2p_1} $ 范数。当 $ p > 6/5 $ 时，可通过选择适当的 $ p_1, p_2 $ 使得整体指数为负，从而保证当 $ R \to \infty $ 时该项趋于零。

  4. 发散定理（Divergence Theorem）
    在定理1.3中，作者将方程改写为守恒律形式，涉及 $ \mathrm{div}(y |H|^q) $ 和 $ \mathrm{div}(J |H|^q) $，然后在环形区域 $ B_{R_2} \setminus B_{R_1} $ 上积分，转化为球面上的面积分。利用 $ H \in L^q $ 推出存在序列 $ R_j \to \infty $ 使得球面积分为零，从而让外部边界项消失。

  5. 反证法与最大值原理思想
    在定理1.4中，作者采用反证法：假设 $ H $ 不具紧支集，则 $ |H|^2 \phi_R $ 在某非紧区域取得正最大值。结合方程导出的微分不等式：

    $$ \partial_s (|H|^2 \phi_R) + \cdots \leq 0, $$

    并利用 $ a > 0 $（因 $\alpha < -2$）、$ b < 0 $，以及 $ \nabla \phi_R \geq 0 $，构造出一个矛盾，从而证明 $ H $ 必须具有紧支集，进而归约为定理1.1的情形。

b) 处理技术难题的关键技巧

  1. 处理霍尔项的代数变形
    霍尔项 $ \nabla \times ((\nabla \times H) \times H) $ 极其复杂。作者将其展开为：

    $$ (H \cdot \nabla) J - (J \cdot \nabla) H, $$

    其中 $ J = \nabla \times H $。这一分解使非线性项更易于与测试函数配对，特别是在计算 $ ((\nabla \times H) \times H) \cdot \nabla \times (H \phi_R) $ 时，能有效分离出可被控制的边界项。

  2. 利用周期性消除时间边界项
    所有能量恒等式中出现的时间导数项 $ \int \partial_s (\cdots) dy ds $，由于 $ H $ 是周期函数，其在 $ s=0 $ 与 $ s=S_0 $ 的积分相等，故总积分为零。这是整个方法得以成立的前提。

  3. 多参数协调选取技巧
    在定理1.2中，为使估计闭合，需协调选取 $ p, p_1, p_2 $ 满足 $ 1/p + 1/p_1 + 1/p_2 = 1 $，并确保 $ R $ 的幂次为负。作者指出只要 $ p > 6/5 $，总能找到这样的组合，体现了对函数空间精细结构的把握。

  4. 从局部极值出发的结构性论证（定理1.5）
    此处未用能量法，而是基于 $ J(0,s)=0 $ 且为局部极值点，结合方程在原点处的展开，推导出 $ H $ 在原点附近必须为零，再通过唯一延拓性得出全局为零。这是一种典型的局部到整体推理。

  5. 非衰减条件下的控制（定理1.6）
    即便 $ H $ 不趋于零，只要 $ |H| = o(|y|^2) $，作者仍能构造合适的估计。这表明增长率本身也蕴含信息，突破了传统“必须衰减”的思维定势。

与之前工作的比较

本文是对已有工作的显著推进，尤其相对于早期关于自相似解的研究（如[8]），实现了从“连续自相似”到“离散自相似”的跨越，并放宽了多个关键假设。

a) 突破与改进之处

  1. 研究对象的推广：从自相似到离散自相似
    以往工作（如[8]）主要排除连续自相似解（即 $ H $ 与 $ s $ 无关），这类解要求更强的对称性。而本文研究的是更一般的离散自相似解，对应 $ H(y,s) $ 为周期函数。这类解在物理上更合理（无需无限精细尺度对称），数学上也更难处理（因时间依赖性引入额外复杂性）。因此，本文的结果覆盖了更广泛的奇性解候选者。

  2. 条件更弱，适用范围更广

    • 定理1.1 中仅要求 $ J = o(|y|) $，而非 $ J \in L^p $ 或更高正则性，这是比许多现有结果更弱的渐近条件。
    • 定理1.3 和 1.4 不要求 $ H $ 本身衰减，而只限制 $ J $ 和 $ \nabla J $ 的行为，更具物理意义（电流决定电磁力）。
    • 定理1.6 完全放弃衰减假设，转而使用增长速率控制，这是全新的思路，在同类问题中罕见。

  3. 维度与尺度参数的完整覆盖
    文章系统地讨论了不同 $ \alpha $ 区间（$ \alpha > 3/2 $、$ -2 < \alpha < 3/2 $、$ \alpha < -2 $、$ \alpha > 0 $ 等），几乎涵盖了所有可能的尺度情形，展现出完整的分类思想。

  4. 方法论上的统一框架
    尽管各定理条件不同，但作者构建了一个统一的能量估计框架，灵活运用截断函数、Hölder不等式、Sobolev嵌入和极限分析，展示了强大而系统的分析能力。

b) 相较于现有结果的优势

相比经典Euler方程或Navier-Stokes中关于自相似解的研究（如[4,5,6]），本文处理的是更高阶奇异性的非线性项（霍尔项），且缺乏任何耗散机制，因此技术难度更大。作者成功将能量方法适配于这种强非线性系统，证明了即使在无粘情况下，某些类型的奇性解仍不可能存在。

此外，本文结果为electron MHD 的全局正则性猜想提供了支持：既然常见的奇性试探解均被排除，或许该系统本身是整体光滑的。这一方向未来有望结合数值模拟进一步验证。

总之，本文不仅在技术上有所创新，更在概念层面拓展了我们对奇性形成的理解，为强非线性PDE的奇性分析提供了新范例。

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

1. Joint filtration of compressible immiscible liquids

可压缩不相溶液体的联合渗流

作者: Anvarbek Meirmanov

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

发布日期: 2025-10-22

arXiv ID: 2510.19502v1

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

摘要翻译

中文标题: 可压缩不相溶液体的联合渗流

摘要翻译: 本文考虑在固体骨架孔隙空间中两种不相溶黏性液体共同运动的数学模型所导出的初边值问题。首先，我们在微观尺度下研究这一物理过程（孔隙尺度 $\varepsilon$ 约为 5–20 微米），该过程由模型 $\textbf{A}^{\varepsilon}$ 所描述，该模型包含描述固体骨架的Lamé方程和描述液体组分的Stokes方程。接着，假设对应于微观尺度下逼近模型 $\textbf{A}^{\varepsilon}$ 的模型 $\textbf{B}^{\varepsilon}$ 的初边值问题存在弱解，并通过均匀化方法，推导出在宏观尺度下描述弹性固体骨架中两种弱黏性不相溶流体共同运动的Biot模型。

2. Boundary controllability of the Korteweg-de Vries equation: The Neumann case

Korteweg-de Vries方程的边界可控制性：Neumann情形

作者: R. de A. Capistrano-Filho, J. S. da Silva

类别: math.AP, math.OC

发布日期: 2023-10-08

arXiv ID: 2310.04977v3

链接: http://arxiv.org/abs/2310.04977v3

摘要翻译

中文标题: Korteweg-de Vries方程的边界可控制性：Neumann情形

摘要翻译: 本文在区间 $[0,L]$ 上考虑具有Neumann边界条件且仅含一个控制输入的Korteweg-de Vries（KdV）方程，对其临界集现象的理解提供了必要的第一步。我们证明了在临界情形下KdV方程是可控的，即当空间区域长度 $L$ 属于集合 $\mathcal{R}_c$ 时（其中 $c \neq -1$），KdV方程在空间 $L^2(0,L)$ 中是精确可控制的，这里

$$ \mathcal{R}_c:=\left\{\frac{2\pi}{\sqrt{3(c+1)}}\sqrt{m^2+ml+l^2};\ m,l\in \mathbb{N}^*\right\}\cup\left\{\frac{m\pi}{\sqrt{c+1}};\ m\in \mathbb{N}^*\right\}. $$

该结果通过运用返回法（return method）结合不动点论证得出。

不相关论文（基本信息）

1. On the singularities of the spectral shift function for some tight-binding models

关于某些紧束缚模型谱位移函数奇点的研究

作者: Marouane Assal, Olivier Bourget, Diomba Sambou, Amal Taarabt

类别: math.SP, math.AP, math.FA, 35J10, 81Q10, 35P20, 35P25, 47A10, 47A11, 47A55, 47F05

发布日期: 2024-09-20

arXiv ID: 2409.13942v2

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

摘要: We consider perturbed discrete tight-binding models in $\ell^2(\mathbb{Z_h},\mathcal{G})$ describing union of quantum particles with localized interactions, where $\mathbb{Z_h}$ is the 1D lattice $h\mathbb{Z_h}$, $h > 0$, and $\mathcal G$ is a separable Hilbert space. The perturbations play the role of self-adjoint relatively compact (matrix-valued) electric potentials with $\mathcal B(\mathcal G)$-valued coefficients decaying polynomially at infinity. We analyze the Spectral Shift Function (SSF) associated to the pair of the perturbed and the unperturbed operators. On the one hand, we show that the SSF is bounded near the spectral thresholds of the essential spectrum if $\dim(\mathcal G) < +\infty$. On the other hand, if $\dim(\mathcal G) = +\infty$, we show that it may have singularities at some thresholds points $\mu$ of the essential spectrum. In particular, new mechanisms allowing the SSF to have singularities at the thresholds are exhibited, based on the degeneracy of the spectrum of the unperturbed operator. Moreover, we give the main terms of the asymptotic behaviors of the SSF near $\mu$ described in terms of some explicit effective Berezin-Toeplitz type operators. These results are completed by Levinson type formulas and examples of eigenvalues asymptotics for power-like and exponential decay potentials.

2. Uniqueness in the Plateau problem for calibrated currents

校准流的Plateau问题解的唯一性

作者: Bryan Dimler, Chen-Kuan Lee

类别: math.DG, math.AP, 49Q15

发布日期: 2025-10-02

arXiv ID: 2510.02299v2

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

摘要: We show that every compactly supported smoothly calibrated integral current with connected $C^{3,\alpha}$ boundary is the unique solution to the oriented Plateau problem for its boundary data. The same holds true for compactly supported $``$continuously calibrated$"$ integral flat chains. This is proved as a consequence of the boundary regularity theory for area-minimizing currents and a unique continuation argument in the spirit of Frank Morgan. In codimension one, the argument yields a sufficient condition for uniqueness in the oriented Plateau problem expressed in terms of the regularity of the calibrating form.

3. Observability of Schrödinger propagators on tori in rough settings

粗糙情形下环面上薛定谔传播子的可观测性

作者: Nicolas Burq, Hui Zhu

类别: math.AP, math-ph, math.FA, math.MP, math.OC

发布日期: 2025-09-28

arXiv ID: 2509.23965v2

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

摘要: On tori of arbitrary dimensions, Schr"odinger propagators with bounded potentials are conjectured to be observable from space-time domains of positive Lebesgue measure. We reduce this conjecture to certain integrability bounds for free Schr"odinger waves, thereby proving the conjecture on the one-dimensional torus and producing new examples of observation domains. These bounds are far weaker than Bourgain’s conjectured periodic Strichartz estimates, yet remain highly nontrivial.

4. Second domain variation for the $p$ - capacity and the $q$ - torsional rigidity

$p$-容度与$q$-扭转变形刚度的二阶区域变分

作者: Alfred Wagner

类别: math.AP, math.OC

发布日期: 2025-05-19

arXiv ID: 2505.20307v3

链接: http://arxiv.org/abs/2505.20307v3

摘要: The second domain variation of the $p$-capacity and the $q$ - torsional rigidity for compact sets in $R^d, d\geq3$ with $1<p<d$ is computed. Conditions on $p$ and $q>1$ are given such that the ball is a local minimzer or maximizer of the product.

5. Existence of multi-solitons with any parameters for the 5D energy critical wave equation

对于五维临界能量波动方程，任意参数下多孤子解的存在性

作者: Yvan Martel, Frank Merle

类别: math.AP

发布日期: 2025-10-22

arXiv ID: 2510.19609v1

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

摘要: For the focusing, energy critical wave equation in dimension 5, we construct multi-solitons with any number of solitons, any choice of signs, speeds, scaling parameters and translation parameters. This requires to revisit in depth previous constructions of multi-solitons based on a unidirectional approach, to fully take into account the dimension of the space and the possibility for solitons to move in any direction. Then, as a consequence of this more general construction and of the arguments developed in a previous article, the inelastic nature of any collision of solitons is proved under a non-cancellation assumption on the parameters.

6. Resonance-induced nonlinear bound states

由共振诱导的非线性束缚态

作者: Jackson C. Turner, Michael I. Weinstein

类别: math-ph, math.AP, math.MP, nlin.PS

发布日期: 2025-10-22

arXiv ID: 2510.19538v1

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

摘要: We study nonlinear bound states – time-harmonic and spatially decaying ($L^2$) solutions – of the nonlinear Schr"odinger / Gross–Pitaevskii equations (NLS/GP) with a compactly supported linear potential. Such solutions are known to bifurcate from the $L^2$ bound states of an underlying Schr"odinger operator $H_V=-\partial_x^2+V$. In this article we prove an extension of this result: for the 1D NLS/GP, nonlinear bound states also arise via bifurcation from the scattering resonance states and transmission resonance states of $H_V$, associated with the poles and zeros, respectively, of the reflection coefficients, $r_\pm(k)$, of $H_V$. The corresponding resonance states are non-decaying and only $L^2_{\rm loc}$. In contrast to nonlinear states arising from $L^2$ bound states of $H_V$, these resonance bifurcations initiate at a strictly positive $L^2$ threshold which is determined by the position of the complex scattering resonance pole or transmission resonance zero.

7. Analysis of quasi-periodic waves of cubic nonlinear Schr{ö}dinger equations

三次非线性薛定谔方程拟周期波的分析

作者: Perla Kfoury, Stefan Le Coz, Tai-Peng Tsai

类别: math.AP

发布日期: 2025-10-22

arXiv ID: 2510.19524v1

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

摘要: We study the quasi-periodic standing wave solutions of the focusing and defocusing cubic nonlinear Schr{"o}dinger equations in dimension one. In the defocusing case, we establish a diffeomorphic correspondence between the invariants of the ordinary differential equation of the wave profiles and the conserved quantities of the evolution equation. We introduce a numerical scheme to compute the minimizers of the energy at fixed mass and momentum for both focusing and defocusing cases. The scheme is based on a gradient flow approach with discrete renormalization at each time step. The novelty of our scheme is that the renormalization step deals at the same time with the mass and the momentum constraints. In numerical experiments, we observe that a given solution of the profile ordinary differential equation is also a minimizer of the energy at corresponding mass and momentum.

8. On the number of nodal domains of homogeneous caloric polynomials

关于齐次热多项式节点域数量的研究

作者: Matthew Badger, Cole Jeznach

类别: math.AP, Primary: 35K05. Secondary: 26C05, 26C10, 35R35

发布日期: 2024-01-14

arXiv ID: 2401.07268v2

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

摘要: We investigate the minimum and maximum number of nodal domains across all time-dependent homogeneous caloric polynomials of degree $d$ in $\mathbb{R}^{n}\times\mathbb{R}$ (space $\times$ time), i.e., polynomial solutions of the heat equation satisfying $\partial_t p\not\equiv 0$ and

$$ p(\lambda x, \lambda^2 t) = \lambda^d p(x,t)\quad\text{for all $x \in \mathbb{R}^n$, $t \in \mathbb{R}$, and $\lambda > 0$.} $$

When $n=1$, it is classically known that the number of nodal domains is precisely $2\lceil d/2\rceil$. When $n=2$, we prove that the minimum number of nodal domains is 2 if $d\not \equiv 0\pmod 4$ and is 3 if $d\equiv 0\pmod 4$. When $n\geq 3$, we prove that the minimum number of nodal domains is $2$ for all $d$. Finally, we show that the maximum number of nodal domains is $\Theta(d^n)$ as $d\rightarrow\infty$ and lies between $\lfloor \frac{d}{n}\rfloor^n$ and $\binom{n+d}{n}$ for all $n$ and $d$. As an application and motivation for counting nodal domains, we confirm existence of the singular strata in Mourgoglou and Puliatti’s two-phase free boundary regularity theorem for caloric measure.

9. On singularly perturbed $(p, N )$-Laplace Schrödinger equation with logarithmic nonlinearity

关于具对数非线性的奇异摄动$(p, N)$-Laplace Schrödinger方程

作者: Deepak Kumar Mahanta, Tuhina Mukherjee, Patrick Winkert

类别: math.AP, 35A15, 35D30, 35J10, 35J60, 35J92

发布日期: 2024-08-07

arXiv ID: 2408.03788v2

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

摘要: This article focuses on the study of the existence, multiplicity and concentration behavior of ground states as well as the qualitative aspects of positive solutions for a $(p, N)$-Laplace Schr"{o}dinger equation with logarithmic nonlinearity and critical exponential nonlinearity in the sense of Trudinger-Moser in the whole Euclidean space $\mathbb{R}^N$. Through the use of smooth variational methods, penalization techniques, and the application of the Lusternik-Schnirelmann category theory, we establish a connection between the number of positive solutions and the topological properties of the set in which the potential function achieves its minimum values.

10. $(INV)$ condition and regularity of the inverse

$(INV)$ 条件与逆的正则性

作者: Anna Doležalová, Stanislav Hencl, Jani Onninen

类别: math.FA, math.AP

发布日期: 2024-12-25

arXiv ID: 2412.18976v2

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

摘要: Let $f \colon \Omega \to \Omega' $ be a Sobolev mapping of finite distortion between planar domains $\Omega $ and $\Omega'$, satisfying the $(INV)$ condition and coinciding with a homeomorphism near $\partial\Omega $. We show that $f$ admits a generalized inverse mapping $h \colon \Omega' \to \Omega$, which is also a Sobolev mapping of finite distortion and satisfies the $(INV)$ condition. We also establish a higher-dimensional analogue of this result: if a mapping $f \colon \Omega \to \Omega' $ of finite distortion is in the Sobolev class $W^{1,p}(\Omega, \mathbb{R}^n)$ with $p > n-1$ and satisfies the $(INV)$ condition, then $f$ has an inverse in $W^{1,1}(\Omega', \mathbb{R}^n)$ that is also of finite distortion. Furthermore, we characterize Sobolev mappings satisfying $(INV)$ whose generalized inverses have finite $n$-harmonic energy.

11. Homogenisation of singular SPDEs

奇异随机偏微分方程的均匀化

作者: Martin Hairer, Harprit Singh

类别: math.AP, math.PR, 60H17, 35B27

发布日期: 2025-10-22

arXiv ID: 2510.19339v1

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

摘要: We introduce an approach to study homogenisation of a large class of singular SPDEs of the form

$$ \partial_t u_\varepsilon - \nabla\cdot {A}(x/\varepsilon,t/\varepsilon^2) \nabla u_\varepsilon = F(x/\varepsilon , t/\varepsilon^2, u_\varepsilon , \nabla u_\varepsilon , \xi ) $$

which is based on the idea of importing (classical) homogenisation results into the framework of regularity structures and the insight that one can rewrite the SPDE under consideration in terms of a model, where the correctors (from homogenisation theory) are seen as further abstract noises'. As applications, we establish periodic space-time homogenisation results for oscillatory generalisations of the 2d g-PAM and $\Phi^4_3$ equation proving that when the noise is regularised at scale $\delta\ll 1$ solutions to the equation with coefficient field ${A}(x/\varepsilon,t/\varepsilon^2)$, when appropriately renormalised, converge to solutions to the corresponding homogenised equation along any sequence $(\varepsilon,\delta)\to 0$. We make the observation that the unbounded divergences can be written as sums of two types of terms: small scale’ terms, the spatial dependence of which is an explicit local function of the coefficient field and `large scale’ terms, which for logarithmic divergences are explicit involving the homogenised matrix and correctors. Furthermore, in order to recover the same solution to the corresponding homogenised equation along any joint limit $(\varepsilon, \delta)\to 0$ one has to subtract additional bounded renormalisation constants which appear due to oscillations at mesoscopic scales, as well as due to resonances between the coefficient field and the oscillations in the nonlinearity.

12. Global Well-Posedness and Soliton Resolution for the Half-Wave Maps Equation with Rational Data

具有有理数据的半波映照方程的整体适定性与孤立子分解

作者: Patrick Gérard, Enno Lenzmann

类别: math.AP, math-ph, math.MP, nlin.SI

发布日期: 2024-12-04

arXiv ID: 2412.03351v3

链接: http://arxiv.org/abs/2412.03351v3

摘要: We study the energy-critical half-wave maps equation:

$$ \partial_t \mathbf{u} = \mathbf{u} \times |D| \mathbf{u} $$

for $\mathbf{u} : [0, T) \times \mathbb{R} \to \mathbb{S}^2$. Our main result establishes the global existence and uniqueness of solutions for all rational initial data $\mathbf{u}_0 : \mathbb{R} \to \mathbb{S}^2$. This demonstrates global well-posedness for a dense subset within the scaling-critical energy space $\dot{H}^{1/2}(\mathbb{R}; \mathbb{S}^2)$. Furthermore, we prove soliton resolution for a dense subset of initial data in the energy space, with uniform bounds for all higher Sobolev norms $\dot{H}^s$ for $s > 0$. Our analysis utilizes the Lax pair structure of the half-wave maps equation on Hardy spaces in combination with an explicit flow formula. Extending these results, we establish global well-posedness for rational initial data (along with a soliton resolution result) for a generalized class of matrix-valued half-wave maps equations with target spaces in the complex Grassmannians $\mathbf{Gr}_k(\mathbb{C}^d)$. Notably, this includes the complex projective spaces $ \mathbb{CP}^{d-1} \cong \mathbf{Gr}_1(\mathbb{C}^d)$ thereby extending the classical case of the target $\mathbb{S}^2 \cong \mathbb{CP}^1$.

13. Quantitative lower bound for solutions to the Boltzmann equation in non-convex domains

玻尔兹曼方程在非凸区域解的定量下界估计

作者: Jhe-Kuan Su

类别: math.AP, 35Q20

发布日期: 2025-05-06

arXiv ID: 2505.03396v3

链接: http://arxiv.org/abs/2505.03396v3

摘要: In this article we deal with the continuous mild solutions of the Boltzmann equation in a bounded space domain, either with angular cutoff or non cutoff assumption. Without assuming that the space domain is convex, we construct a Maxwellian lower bound in the cutoff case, and a weaker-than-Maxwellian lower bound in the non-cutoff case, thus extending results of \cite{Bri 1,Bri 2} in which the domain was assumed to be convex.

14. Stability and slow dynamics of an interior spiky pattern in a one-dimensional spatial Solow model with capital-induced labor migration

一类具资本驱动劳动力迁移的一维空间Solow模型中内部尖峰模式的稳定性与慢动力学

作者: Fanze Kong, Jiayi Sun, Shuangquan Xie

类别: math.AP, econ.TH

发布日期: 2025-10-22

arXiv ID: 2510.19204v1

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

摘要: One of the most significant findings in the study of spatial Solow-Swan models is the emergence of economic agglomeration, in which economic activities concentrate in specific regions. Such agglomeration provides a fundamental mechanism driving the spatial patterns of urbanization, labor migration, productivity growth, and resource allocation. In this paper, we consider the one-dimensional spatial Solow-Swan model with capital-induced labor migration, which captures the dynamic interaction between labor and capital through migration and accumulation. Focusing on the regime of sufficiently small capital diffusivity, we first construct an interior spike (spiky economic agglomeration) quasi-equilibrium. Next, we perform the linear stability of the corresponding spike equilibrium by using a hybrid asymptotic and numerical method. We show that a single interior spike remains stable for small reaction-time constants but undergoes a Hopf bifurcation when the constant is sufficiently large, leading to oscillations in spike height (economic fluctuation). Finally, we derive a differential-algebraic system to capture the slow drift motion of quasi-equilibrium (core-periphery shift). Numerical simulations are carried out to support our theoretical studies and reveal some intriguing yet unexplained dynamics.

15. Self-similar solutions of semilinear heat equations with positive speed

具有正速度的半线性热方程的自相似解

作者: Kyeongsu Choi, Jiuzhou Huang

类别: math.AP, math.DG

发布日期: 2024-03-25

arXiv ID: 2403.16641v2

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

摘要: We classify the smooth self-similar solutions of the semilinear heat equation $u_t=\Delta u+|u|^{p-1}u$ in $\mathbb{R}^n\times (0,T)$ satisfying an integral condition for all $p>1$ with positive speed. As a corollary, we prove that finite time blowing up solutions of this equation on a bounded convex domain with $u(\cdot,0)\geq 0$ and $u_t(\cdot,0)\geq 0$ converges to a positive constant after rescaling at the blow-up point for all $p>1$.

16. Second radial eigenfunctions to a fractional Dirichlet problem and uniqueness for a semilinear equation

分数阶狄利克雷问题的第二个径向特征函数及半线性方程的唯一性

作者: Mouhamed Moustapha Fall, Tobias Weth

类别: math.AP

发布日期: 2024-05-03

arXiv ID: 2405.02120v2

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

摘要: We analyze the shape of radial second Dirichlet eigenfunctions of fractional Schr"odinger type operators of the form $(-\Delta)^s +V$ in the unit ball $B$ in $\mathbb{R}^N$ with a nondecreasing radial potential $V$. Specifically, we show that the eigenspace corresponding to the second radial eigenvalue is simple and spanned by an eigenfunction $u$ which changes sign precisely once in the radial variable and does not have zeroes anywhere else in $B$. Moreover, by a new Hopf type lemma for supersolutions to a class of degenerate mixed boundary value problems, we show that $u$ has a nonvanishing fractional boundary derivative on $\partial B$. We apply this result to prove uniqueness and nondegeneracy of positive ground state solutions to the problem $(-\Delta)^s u+\lambda u=u^p$ on ${B}$, $\; u=0$ on $\mathbb{R}^N\setminus B$. Here $s\in (0,1)$, $\lambda\geq 0$ and $p>1$ is strictly smaller than the critical Sobolev exponent.

17. High Contrast Transmission and Fabry-Pérot-type Resonances

高对比度透射与法布里-珀罗型共振

作者: Long Li, Mourad Sini

类别: math.AP

发布日期: 2025-10-21

arXiv ID: 2510.19096v1

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

摘要: It is well known, in the acoustic model, that highly contrasting transmission leads to the so-called Minnaert subwavelength resonance. In this work, we show that such highly contrasting transmissions create not only one resonance but a family of infinite resonances located near the real axis where the first one (i.e. the smallest) is indeed the Minnaert one. This family of resonances are the shifts (in the lower complex plan) of the Neumann eigenvalues of the Laplacian. The well known Minneart resonance is nothing but the shift of the trivial (zero) Neumann eigenvalue of the bubble. These resonances, other than the Minnaert ones, are Fabry-P'erot-type resonances as the generated total fields, in the bubble, are dominated by a linear combination of the Neumann eigenfunctions which, in particular, might create interferences. In addition, we establish the following properties.

  1. We derive the asymptotic expansions, at the second order, of this family of resonances in terms of the contrasting coefficient.
  2. In the time-harmonic regime, we derive the resolvent estimates of the related Hamiltonian and the asymptotics of scattered fields that are uniform in the whole space, highlighting the contributions from this sequence of resonances.
  3. In the time domain regime, we derive the time behavior of the acoustic microresonator at large time-scales inversely proportional to powers of microresonator’s radius.
  4. The analysis shows that near Fabry-P'erot resonances, the mircoresonator exhibits pronounced anisotropy. We believe that such a feature may pave the way for designing anisotropic metamaterials from simple configurations of a single microresonator.

18. Local boundedness of weak solutions to elliptic equations under unbalanced Orlicz growth conditions

在非平衡Orlicz增长条件下的椭圆型方程弱解的局部有界性

作者: Gabriele Giannone

类别: math.AP, 35J60 (Primary) 46E30 (Secondary)

发布日期: 2025-10-21

arXiv ID: 2510.18979v1

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

摘要: We establish sufficient conditions for the local boundedness of weak solutions to a broad class of nonlinear elliptic equations in divergence form, under unbalanced growth conditions on the stress field. Our analysis is carried out in a non-variational setting, with no symmetry or structural assumptions on the operator. The ellipticity and growth are prescribed via distinct Young functions, leading to a Orlicz-type setting that captures a wide class of nonstandard behaviors. As a special case, the theory encompasses and extends the classical results on equations with $p,q$-growth.