生成时间: 2026年05月26日 06:40:38

运行状态: 完成

计划处理论文数量: 49

已完成论文数量: 49

未完成论文数量: 0

重点关注论文数量: 11

了解领域论文数量: 20

不相关论文数量: 18

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

1. 拟周期行波电子层

Quasi-periodic traveling electron layers

作者: Luca Franzoi, Emeric Roulley

类别: math.AP, physics.plasm-ph

发布日期: 2026-05-25

arXiv ID: 2605.25885v1

摘要翻译: 我们考虑一维空间周期性 Vlasov-Poisson 方程，并在对称平坦速度条带附近，构造出小振幅的拟周期行波电子层解，即相空间中条带状的电子斑块。这些解对条带面积的绝大多数取值存在。证明采用 Nash-Moser 构造方法，并结合基于拟微分齐次展开与 KAM 约化的可约性论证。借助于未知函数的适当线性变换，我们揭示了电子斑块问题与经典的、具有物理意义的电子 Euler-Poisson 系统（压强律为 $p(\rho) = \rho^3$）之间的联系。作为该联系的一个直接而非平凡的推论，我们亦得到该模型的小振幅拟周期行波解。据我们所知，本工作提供了动力学模型中拟周期解构造的第一个严格例子。

Comment: 81 pages

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

1. 研究对象和背景

本文研究一维空间周期性 Vlasov-Poisson 方程

$$ \begin{cases} \partial_t f(t,x,v) + v \partial_x f(t,x,v) - E(t,x) \partial_v f(t,x,v) = 0, \\ E(t,x) = \partial_x \varphi(t,x), \\ \partial_x^2 \varphi(t,x) = 1 - \int_{\mathbb{R}} f(t,x,v)\,dv, \\ (t,x,v) \in \mathbb{R}_+ \times \mathbb{T} \times \mathbb{R}, \end{cases} \tag{1.1} $$

其中 $\mathbb{T} = \mathbb{R}/2\pi\mathbb{Z}$。该方程描述无碰撞等离子体中电子在均匀离子背景下的动力学演化；$f(t,x,v) \ge 0$ 表示相空间 $(x,v)$ 上电子的概率密度，满足中性条件 $\frac{1}{2\pi}\int_0^{2\pi} \int_{\mathbb{R}} f(t,x,v)\,dv\,dx = 1$。核心研究对象是‘电子层’（electron layers）——即形如

$$ f(t,x,v) = \frac{2\pi}{|S_t|} \mathbf{1}_{S_t}(x,v), \quad S_t = \{ (x,v) \in \mathbb{T} \times \mathbb{R} : v_-(t,x) < v < v_+(t,x) \} $$

的弱解，其中 $v_\pm(t,\cdot) \in C^0(\mathbb{T})$ 是 $2\pi$-周期剖面函数，且面积守恒 $|S_t| = \int_0^{2\pi} [v_+(t,x)-v_-(t,x)]\,dx$ 恒成立。特别地，对称平坦条带 $S_{\text{flat}}(a) = \mathbb{T} \times [-a,a]$ 对应平衡态

$$ f_a(x,v) = \frac{1}{2a} \mathbf{1}_{S_{\text{flat}}(a)}(x,v), \quad a > 0, $$

其扰动形式为 $v_\pm(t,x) = \pm a + \varepsilon r_\pm(t,x)$，其中 $\varepsilon \in (0,1)$ 为小参数，$\int_0^{2\pi} r_\pm(t,x)\,dx = 0$。代入 (1.1) 并利用 Poisson 方程消去电势后，可导出变形函数 $(r_+,r_-)$ 满足的封闭系统：

$$ \begin{cases} \partial_t r_+(t,x) = -\partial_x \left[ \frac{1}{2\varepsilon} (a + \varepsilon r_+)^2 - \frac{1}{2a} \partial_x^{-1} r_+ + \frac{1}{2a} \partial_x^{-1} r_- \right], \\ \partial_t r_-(t,x) = -\partial_x \left[ \frac{1}{2\varepsilon} (-a + \varepsilon r_-)^2 - \frac{1}{2a} \partial_x^{-1} r_+ + \frac{1}{2a} \partial_x^{-1} r_- \right]. \end{cases} \tag{1.7} $$

此系统具有哈密顿结构：$\partial_t r = \frac{1}{\varepsilon^2} J \nabla E(r)$，其中 $r = (r_+,r_-)$，$J = \operatorname{diag}(-\partial_x, \partial_x)$，能量泛函 $E$ 包含动能与势能之和。问题的重要性在于：一方面，经典 Landau 阻尼理论表明解析/Gevery 类初值会衰减至齐次平衡态，而电子层属于低正则性解，其长期行为是否仍服从阻尼尚不明确；另一方面，BGK 波与周期性 E-态（如 [87] 中构造的行波）已被证实存在，但它们仅填充环面 $\mathbb{T}^1$，动力学维度受限。本文试图填补的关键空白是：首次严格构造高维不变环面 $\mathbb{T}^d$（$d \ge 2$）上拟周期行波解，即寻找形如 $r(t,x) = \check{r}(\omega t - \vec{\jmath} x)$ 的解，其中 $\omega \in \mathbb{R}^d$ 为频率向量，$\vec{\jmath} \in \mathbb{Z}^d$ 为速度向量，且满足非共振条件 $\omega \cdot \ell \neq 0$（$\forall \ell \in \mathbb{Z}^d \setminus \{0\}$）。这要求突破传统 bifurcation 方法的局限（仅适用于单频），直面由线性算子谱 $\{ \Omega_j(a) = j/|j| \sqrt{a^2 j^2 + 1} \}_{j \in \mathbb{Z}^*}$ 导致的严重小分母问题：当 $\ell \in \mathbb{Z}^d$ 满足 $\sum_{k=1}^d \ell_k \Omega_{j_k}(a) \approx 0$ 时，线性化算子 $L_\omega = \omega \cdot \partial_\varphi - \mathcal{L}$ 在 Fourier 模式 $\ell$ 上的本征值趋于零，导致标准 Newton 迭代失效。已有工作（如 [87]）仅处理了 $d=1$ 的周期情形，而拟周期情形需同时控制无限多个共振集，并保证解保持可逆性（reversibility）与动量守恒（momentum-preserving），这是此前动能模型中从未实现的。

2. 主要定理或主要结果

定理 1.1（拟周期行波电子层）：设 $0 < a_0 < a_1$，$(d_+, d_-) \in \mathbb{N}^2 \setminus \{(0,0)\}$，并固定有限 Fourier 模式集 $S_+ = \{j_1^+, \dots, j_{d_+}^+\} \subset \mathbb{N}^*$、$S_- = \{j_1^-, \dots, j_{d_-}^-\} \subset \mathbb{N}^*$，满足 $S_+ \cap S_- = \emptyset$。给定振幅向量 $(a_j^+)_{j \in S_+} \in \mathbb{R}^{d_+}$、$(a_j^-)_{j \in S_-} \in \mathbb{R}^{d_-}$，则存在 $s_0 = s_0(d, q_0) > 0$ 与 $\varepsilon_0 > 0$，使得对任意 $\varepsilon \le \varepsilon_0$，存在一个 Cantor 型参数集 $C(\varepsilon) \subset (a_0, a_1)$ 满足 $\lim_{\varepsilon \to 0} |C(\varepsilon)| = a_1 - a_0$；对每个 $a \in C(\varepsilon)$，Vlasov-Poisson 方程 (1.1) 存在可逆拟周期行波电子层解

$$ f(t,x,v) = \frac{1}{2a} \mathbf{1}_{S_t}(x,v), \quad S_t = \{ (x,v) \in \mathbb{T} \times \mathbb{R} : -a + \varepsilon r_-(t,x) < v < a + \varepsilon r_+(t,x) \}, $$

其中变形函数 $r = (r_+, r_-)$ 具有如下结构：

$$ r(t,x) = \sum_{j \in S_+} \frac{a_j^+ \cos\big(\Omega_j^+(a,\varepsilon)t - j x\big)}{\sqrt{1 - b_j^2(a)}} \begin{pmatrix} 1 \\ b_j(a) \end{pmatrix} + \sum_{j \in S_-} \frac{a_j^- \cos\big(\Omega_j^-(a,\varepsilon)t - j x\big)}{\sqrt{1 - b_j^2(a)}} \begin{pmatrix} b_j(a) \\ 1 \end{pmatrix} + p(t,x), $$

其中 $b_j(a) = \frac{1}{2a|j|} \left( \sqrt{a^2 j^2 + 1} + a|j| + \frac{1}{2a|j|} \right) \in (0,1)$，扰动项 $p(t,x) = \check{p}(\omega(a,\varepsilon)t - \vec{\jmath} x)$ 属于 Sobolev 空间 $H^{s_0}(\mathbb{T}^d; \mathbb{R}^2)$，满足偶性 $\check{p}(-\varphi) = \check{p}(\varphi)$ 及小性估计 $\|\check{p}\|_{H^{s_0}} = O(\varepsilon)$。频率向量为 $\omega(a,\varepsilon) = \big( \Omega_{j_1^+}^+(a,\varepsilon), \dots, \Omega_{j_{d_+}^+}^+(a,\varepsilon), \Omega_{j_1^-}^-(a,\varepsilon), \dots, \Omega_{j_{d_-}^-}^-(a,\varepsilon) \big)$，速度向量为 $\vec{\jmath} = (j_1^+, \dots, j_{d_+}^+, j_1^-, \dots, j_{d_-}^-)$，且扰动频率满足极限 $\lim_{\varepsilon \to 0} \Omega_j^\kappa(a,\varepsilon) = \kappa \sqrt{a^2 j^2 + 1}$（$\kappa = \pm$）。此外，该解是线性稳定的：线性化算子在拟周期行波解上的 Floquet 指数纯虚。

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

证明的核心障碍是线性化算子 $L_\omega$ 的严重小分母问题与拟周期结构的耦合：在 Nash-Moser 迭代第 $n$ 步，需构造近似右逆 $T_n$ 满足 tame 估计 $\|T_n u\|_{s+\mu} \lesssim \|u\|_{s} + \|i_n\|_{s+\mu'} \|u\|_{s_0}$（$\mu, \mu'$ 为固定损失阶），而直接求逆 $L_\omega^{-1}$ 会导致无穷阶导数损失。作者的真正创新点在于将 KAM 约化与 pseudo-differential 算子理论深度耦合，构建了一个多尺度约化框架：首先，通过 symplectic 变换 $Q$（由 (1.13) 定义的 Fourier 乘子矩阵族 $Q_j(a)$ 组成）将原系统对角化为 $\partial_t \tilde{r} = -i \Omega(a,D) \tilde{r} + \varepsilon J \nabla P H(\tilde{r})$，其中 $\Omega(a,D)$ 是主部 Fourier 乘子；其次，在 action-angle-normal 坐标系下（见 (1.22)-(1.23)），将线性化算子 $L_\omega$ 分解为切向（tangential）与法向（normal）两部分，关键约化集中于法向算子 $L_\omega = \Pi_\perp^S (L^{(0)} - \varepsilon \partial_x R) \Pi_\perp^S$（见 (1.27)）。此处，$L^{(0)}$ 具有展开式 $L^{(0)} = \omega \cdot \partial_\varphi I_2 + i \Omega(a,D) + \varepsilon \partial_x F + R^{(d)}_0 + R^{(o)}_0 + S_{0,M}$（见 (1.28)），其中 $F = \operatorname{diag}(f_+, f_-)$ 是对角乘法算子，$R^{(d)}_0$ 是 diagonal smoothing 余项，$R^{(o)}_0$ 是 off-diagonal 余项，$S_{0,M}$ 是 $M$-阶 smoothing 算子。作者采用三阶段约化：(i) Block-decoupling（6.2 节）：利用 $Q$ 变换的结构，将 $L^{(0)}$ 分块对角化至 $O(\varepsilon^2)$ 精度，消除非对角耦合；(ii) Transport reduction & Egorov argument（6.3 节）：引入 symplectic 微分同胚 $\Phi$ 将 transport 项 $\omega \cdot \partial_\varphi + i \Omega(a,D)$ 共轭为常系数算子 $\omega \cdot \partial_\varphi + i \widetilde{\Omega}(a)$，其中 $\widetilde{\Omega}(a)$ 是 $\Omega(a,D)$ 的平均值，此步依赖 Egorov 定理处理 pseudo-differential 算子的共轭变换；(iii) Diagonalization of $L_\omega$（6.5 节）：对剩余的 $\varepsilon$-阶扰动项 $\varepsilon \partial_x F + R^{(d)}_0$，应用 KAM 迭代型约化，构造一系列近恒同 symplectic 变换 $\Psi^{(k)}$，使得 $\big(\Psi^{(k)}\big)^{-1} L_\omega \Psi^{(k)} = \omega \cdot \partial_\varphi I_2 + i D^{(k)} + R^{(k)}$，其中 $D^{(k)}$ 为对角矩阵，$R^{(k)}$ 的范数以几何速率 $\rho^k$ 衰减。最终，$L_\omega$ 被约化为一个可显式求逆的 diagonal operator 加上一个足够光滑的余项，从而获得满足 tame 估计的右逆 $T_n$。整个框架的枢纽是 pseudo-differential homogeneous expansion：所有算子（如 $F, R^{(d)}_0$）均被表示为齐次符号 $a(x,\xi;\varepsilon)$ 的 Weyl 量化，其符号满足 $a(x,\xi;\varepsilon) \sim \sum_{k \ge 0} \varepsilon^k a_k(x,\xi)$，且各阶 $a_k$ 具备明确的奇偶性与对称性，以兼容可逆性与动量守恒约束。

4. 与之前工作的比较

与既有工作相比，本文实现了三重实质性推进：(1) 解的类型与维度：此前关于 Vlasov-Poisson 方程的周期解构造（如 [87]）仅限于 $d=1$ 的行波，其频率集合为单元素 $\{\Omega_j(a)\}$；而本文首次构造了 $d \ge 2$ 的拟周期解，频率向量 $\omega \in \mathbb{R}^d$ 由 $d$ 个独立模态 $\{\Omega_{j_k}^\kappa(a,\varepsilon)\}$ 组成，解填充 $\mathbb{T}^d$ 不变环面，动力学复杂度呈指数增长。代价是必须将参数 $a$ 限制在满测度 Cantor 集 $C(\varepsilon)$ 上，以满足 Diophantine 条件 $|\omega \cdot \ell| \gtrsim \gamma |\ell|^{-\tau}$（$\forall \ell \in \mathbb{Z}^d \setminus \{0\}$），而 [87] 中的周期情形无需此类限制。(2) 技术工具与假设：[87] 使用局部/全局分歧理论，依赖于系统在平衡态处的 Fredholm 性质及 Crandall-Rabinowitz 定理，要求非退化条件 $\partial_c (c - \Omega_j(a)) \neq 0$；本文则采用 Nash-Moser-KAM 框架，完全绕过 Fredholm 假设，转而要求初始数据处于高正则 Sobolev 空间 $H^s$（$s > s_0$），并利用面积参数 $a$ 作为外部调谐参数来打破共振，而非依赖波速 $c$（如 [87] 中的 $c_m^\pm(a)$）。(3) 物理意义与适用范围：本文揭示了 Vlasov-Poisson 系统与电子 Euler-Poisson 系统（具三次压强律 $p(\rho) = \rho^3$）的精确等价性（见第 2 节），因此定理 1.1 直接蕴含后者的小振幅拟周期行波解（定理 2.1），这是流体力学中 Euler-Poisson 模型的首个拟周期解结果。相比之下，水波或二维 Euler 方程的拟周期解（如 [17, 13, 58]）均基于几何参数（如涡旋半径、水深）调谐，而本文首次将相空间几何参数 $a$（条带半宽） 用作 KAM 调谐参数，建立了动能论与连续介质力学间的深层桥梁。局限性在于：目前结果仅适用于一维空间周期情形；Dirichlet 或反射边界下的 patch 解尚未定义，且三维 Vlasov-Poisson 的拟周期构造仍是开放问题。从提供内容中无法确认本文方法能否推广至非对称条带（即 $v_+$ 与 $v_-$ 的均值不为 $\pm a$），因作者明确声明‘为简化研究，仅考虑对称情形’。

2. 一维标量守恒律色散边界层的稳定性

Stability of dispersive boundary layers for scalar conservation laws in one space dimension

作者: Paolo Antonelli, Pierangelo Marcati, Laura V. Spinolo

类别: math.AP, math-ph

发布日期: 2026-05-25

arXiv ID: 2605.25853v1

摘要翻译: 我们研究一类定义在半直线上的 Korteweg–de Vries (KdV) 型初边值问题的零色散极限，其在 $x=0$ 处施加 Dirichlet 边界条件。我们聚焦于流出（outflow）情形，即极限标量守恒律的解不满足色散问题所施加的边界条件。我们构造了一个依赖于快变量的边界层剖面，该剖面由相应的定常三阶边界层方程唯一确定，其确定依据为边界条件之间的失配，以及该剖面在快变量趋于无穷时的指数衰减性。我们的主要结果表明：在数据满足适当的正则性与相容性假设下，色散解可被一个 WKB 展开式良好逼近，该展开式由守恒律的光滑解与该边界层剖面之和构成。特别地，我们通过在加权能量范数下对余项建立定量估计，证明了该边界层剖面的稳定性，并进一步证明该余项在 $H^1$ 中收敛于 $0$，且该收敛在时间上是一致的，持续至守恒律光滑解的生命史（lifespan）结束为止。该证明基于线性化能量泛函的分析，不依赖于完全可积性或逆散射方法。该结果适用于一般通量函数，且对边界层振幅无小性假设。据我们所知，这是首个关于半直线上 KdV 型方程边界层的稳定性结果。

Comment: 27 pages

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

1. 研究对象和背景

本文研究半直线 $\mathbb{R}_+ = (0, +\infty)$ 上一类 Korteweg–de Vries (KdV) 型色散方程在零色散极限 $\varepsilon \to 0^+$ 下产生的边界层的稳定性问题。核心对象是如下初边值问题：

$$ \begin{cases} \partial_t u^\varepsilon + \partial_x [f(u^\varepsilon)] + \varepsilon^2 \partial_{xxx} u^\varepsilon = 0, & (t,x) \in \mathbb{R}_+ \times \mathbb{R}_+, \\ u^\varepsilon|_{t=0} = u^{\mathrm{in}} \in H^4(\mathbb{R}_+), \\ u^\varepsilon|_{x=0} = u^b \in H^4([0,T_0]). \end{cases} \tag{1.1} $$

其中非线性通量函数 $f \in C^4(\mathbb{R})$ 满足严格递减性假设

$$ f'(u) \le -c < 0, \quad \forall u \in J \subseteq \mathbb{R}, \tag{1.2} $$

即特征速度为负，所有特征线向左传播（斜率为负），从而边界 $x=0$ 是出流边界（outflow boundary）。该假设排除了入流或混合边界情形，是本文分析成立的关键几何前提。

当 $\varepsilon \to 0^+$ 时，(1.1) 形式上趋近于标量守恒律初值问题

$$ \begin{cases} \partial_t u + \partial_x [f(u)] = 0, & (t,x) \in \mathbb{R}_+ \times \mathbb{R}_+, \\ u|_{t=0} = u^{\mathrm{in}}, \end{cases} \tag{1.3} $$

但注意：(1.3) 在半直线上的标准初值问题不配备边界条件；由于 $f'<0$，其解 $u(t,x)$ 在 $x=0$ 处的迹 $u_0(t) := u(t,0)$ 完全由初始数据 $u^{\mathrm{in}}$ 通过特征线决定（参见 Dafermos [11]），一般与给定的色散边界条件 $u^b(t)$ 不一致，即 $u_0(t) \neq u^b(t)$。这种不匹配（mismatch）导致在 $x=0$ 附近出现尺度为 $\mathcal{O}(\varepsilon)$ 的快速过渡层——即色散边界层（dispersive boundary layer）。

为刻画该层，作者引入快变量 $y = x/\varepsilon$，并构造边界层剖面 $V(t,y)$，它满足静止的三阶常微分方程（即边界层方程）：

$$ \begin{cases} \partial_y \big[ f(u_0(t) + V) - f(u_0(t)) \big] + \partial_{yyy} V = 0, & y > 0, \\ V(t,0) = u^b(t) - u_0(t), \\ \lim_{y \to +\infty} V(t,y) = 0. \end{cases} \tag{1.5} $$

该方程是将原PDE在边界附近作WKB展开 $u^\varepsilon \approx u(t,x) + V(t,x/\varepsilon)$ 后，提取 $\mathcal{O}(\varepsilon^0)$ 项所得的主导平衡。关键在于：(1.5) 是一个非自治、非线性、三阶ODE边值问题，其适定性本身即构成一个非平凡分析问题。已有文献（如 Xin [24] 对粘性边界层的研究）表明，此类问题的解存在唯一且指数衰减，但其证明严重依赖于二阶扩散算子 $-\partial_{xx}$ 的正定性与最大值原理；而此处的三阶导数 $\partial_{yyy}$ 既非椭圆亦非抛物，不具备类似结构，传统方法失效。因此，本文首要空白在于：在缺乏完全可积性（如KdV逆散射）和强凸性（如 $f''>0$）假设下，建立 (1.5) 解的存在性、唯一性及定量指数衰减估计（见引理2.1–2.3），并进而证明该剖面在色散扰动下的非线性稳定性——即验证WKB展开的余项 $w^\varepsilon := u^\varepsilon - u - V(\cdot, x/\varepsilon)$ 在 $\varepsilon \to 0^+$ 时以明确定量阶趋于零。这一问题在色散PDE边界层理论中长期悬而未决：此前工作 [9,18] 针对非线性Schrödinger方程，而KdV型方程的半直线边界层稳定性尚属空白（作者明确声明“This is the first stability result for boundary layers of KdV-type equation on the half line”）。其困难本质在于：三阶色散项破坏能量守恒结构，标准能量法直接应用会导致符号不定项；同时，出流边界下无法施加传统Neumann或Robin型控制，必须仅凭Dirichlet条件与指数衰减渐近条件实现全局控制。

2. 主要定理或主要结果

定理1.1是全文核心结果，其结论分为两部分：

  1. 适定性：在假设(1.2)、(1.4)、(1.8)、(1.9)下，存在 $\varepsilon_* > 0$，使得对任意 $\varepsilon \in (0, \varepsilon_*]$，问题(1.1)在时间区间 $[0,T_0]$ 上存在唯一解 $u^\varepsilon \in C^0([0,T_0]; H^2(\mathbb{R}_+)) \cap L^\infty([0,T_0]; H^3(\mathbb{R}_+))$。

  2. 边界层稳定性（主结论）：定义WKB展开余项 $w^\varepsilon(t,x) := u^\varepsilon(t,x) - u(t,x) - V(t,x/\varepsilon)$，其中 $u$ 满足(1.3)，$V$ 满足(1.5)。则存在常数 $K_* > 0$（仅依赖于 $T_0$, $\|u^b - u_0\|_{C^3}$, $\|f\|_{C^3}$, $c$, $\|u\|_{C^0(H^4)}$），使得对任意 $t \in [0,T_0]$，有加权能量估计

    $$ \int_{\mathbb{R}_+} \left[ (w^\varepsilon)^2 + \frac{\varepsilon^2}{2} (\partial_x w^\varepsilon)^2 \right](t,x) \, dx \le K_* \varepsilon^3. \tag{1.10} $$

    由此立即推出 $w^\varepsilon \to 0$ 在 $C^0([0,T_0]; H^1(\mathbb{R}_+))$ 中一致成立。特别地，该估计的误差阶为 $\mathcal{O}(\varepsilon^{3/2})$ 在 $H^1$ 范数下（因 $\|w^\varepsilon\|_{H^1}^2 \lesssim \varepsilon^3$），且时间范围达光滑解 $u$ 的整个生命期 $T_0$，无需小初值假设。此外，引理2.1–2.3给出了边界层剖面 $V$ 的精细定量性质：对任意 $t \in [0,T_0]$，有

    $$ |V(t,y)| \le \|u^b - u_0\|_{C^0} e^{-\sqrt{c}\, y}, \quad |\partial_t V(t,y)| \le \hat{C} e^{-(\sqrt{c}/2) y}, \quad |\partial_{tt} V(t,y)| \le \hat{C} e^{-(\sqrt{c}/4) y}, \quad |\partial_{ttt} V(t,y)| \le \hat{C} e^{-(\sqrt{c}/8) y}, $$

    其中 $\hat{C} = C(T_0, \|u^b - u_0\|_{C^3}, \|f\|_{C^3}, c, \|u\|_{C^0(H^4)})$，这些指数衰减率随时间导数阶数加倍而减半，是后续能量估计中控制源项 $E_b$ 和 $E_{\mathrm{inn}}$ 的关键。

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

证明的核心障碍在于：标准能量法应用于余项 $w^\varepsilon$ 的方程(3.3)时，三阶色散项 $\varepsilon^2 \partial_{xxx} w^\varepsilon$ 与一阶对流项 $\partial_x[f'(u^\varepsilon_a) w^\varepsilon]$ 耦合，导致直接乘以 $w^\varepsilon$ 或 $\partial_x w^\varepsilon$ 均无法获得符号确定的能量泛函。作者的根本创新在于构造并分析一个线性化能量泛函

$$ E(t) := \int_{\mathbb{R}_+} \left[ \frac{\varepsilon^2}{2} (\partial_x w^\varepsilon)^2 - f'(u^\varepsilon_a) (w^\varepsilon)^2 \right] dx, \tag{3.6} $$

该泛函并非原始系统的哈密顿量，而是围绕近似解 $u^\varepsilon_a = u + V(\cdot,x/\varepsilon)$ 的线性化能量。其关键性质是：由(1.2)知 $-f'(u^\varepsilon_a) \ge c > 0$，故第二项正定，整体 $E(t)$ 等价于 $\|w^\varepsilon\|_{H^1}^2$ 的加权范数（即(1.10)左端），从而控制 $E(t)$ 即得所需估计。

为控制 $E(t)$，作者采用反导数技巧（anti-derivative method）：定义 $W^\varepsilon(t,x) := \int_0^x w^\varepsilon(t,x') dx'$，则 $\partial_x W^\varepsilon = w^\varepsilon$，且由(3.3)的齐次初边值条件得 $W^\varepsilon|_{t=0}=0$、$W^\varepsilon|_{x=0}=0$。将余项方程(3.3)两边乘以 $\partial_t W^\varepsilon$ 并在 $[0,T] \times \mathbb{R}_+$ 上积分，得到

$$ 0 = \int_0^T \int_{\mathbb{R}_+} \partial_t W^\varepsilon \left[ \partial_t w^\varepsilon + \partial_x(f'(u^\varepsilon_a) w^\varepsilon + g(w^\varepsilon)^2) + \varepsilon^2 \partial_{xxx} w^\varepsilon + E_{\mathrm{inn}} + E_b \right] dx dt. $$

其中，$I_2 = \int_0^T \int_{\mathbb{R}_+} \partial_t W^\varepsilon \partial_x(f'(u^\varepsilon_a) w^\varepsilon) dx dt$ 和 $I_4 = \int_0^T \int_{\mathbb{R}_+} \partial_t W^\varepsilon \varepsilon^2 \partial_{xxx} w^\varepsilon dx dt$ 是主导项，经分部积分后精确对应 $\frac{d}{dt} E(t)$ 的主要部分。其余项（如 $I_1, I_3, I_5, I_6$）需被吸收：$I_1$ 给出 $\frac{1}{2} \|w^\varepsilon(T)\|_{L^2}^2$；$I_5, I_6$ 涉及源项 $E_{\mathrm{inn}} = \varepsilon^2 \partial_{xxx} u$ 和 $E_b$，其大小由 $V$ 的指数衰减（引理2.2–2.3）与 $u$ 的 $H^4$ 正则性控制，最终贡献 $\mathcal{O}(\varepsilon^3)$；$I_3$ 含非线性项 $g(w^\varepsilon)^2$，利用连续性归纳法（continuity argument）定义 $T^* := \sup\{ t \in [0,T_0] : \|w^\varepsilon(\tau)\|_{H^1}^2 \le \varepsilon^2, \, \forall \tau \in [0,t] \}$，结合Sobolev嵌入 $\|w^\varepsilon\|_{L^\infty} \le \sqrt{2} \varepsilon^{1/2}$（由(3.11)），可证 $g(w^\varepsilon)^2$ 项为高阶小量 $\mathcal{O}(\varepsilon^{5/2})$，可被主导项吸收。最终，通过Grönwall型不等式与归纳法闭合，得出 $T^* = T_0$，从而完成(1.10)的证明。整个框架不依赖逆散射或完全可积性，适用于一般通量 $f$，且未对边界层振幅作小性假设，这是区别于以往工作的本质技术突破。

4. 与之前工作的比较

本文与既有工作存在多维度根本性差异：

  1. 与Xin [24]（粘性边界层）的比较：Xin针对方程 $\partial_t u^\varepsilon + \partial_x f(u^\varepsilon) = \varepsilon \partial_{xx} u^\varepsilon$ 建立了粘性边界层稳定性，其证明依赖二阶算子 $-\partial_{xx}$ 的正定性及最大值原理，可自然导出 $L^\infty$ 衰减和 $H^1$ 能量耗散。而本文处理三阶色散项 $\varepsilon^2 \partial_{xxx}$，其符号不定性迫使作者放弃经典能量法，转而构造线性化泛函(3.6)并采用反导数技巧。Xin的结果分两种情形：i) 小振幅层在 $[0,\hat{T}]$（$\hat{T} < T_0$）上稳定；ii) 一般振幅层在 $[0,T_0]$ 上稳定但需 $f$ 凸（$f''>0$）。本文则在无小性假设、无凸性要求下，直接在完整生命期 $[0,T_0]$ 上建立稳定性，代价是仅适用于出流边界（$f'<0$），且指数衰减率随导数阶数降低（引理2.3）。

  2. 与KdV方程整体线研究（如 [20]）的比较：在全直线 $\mathbb{R}$ 上，KdV的零色散极限研究高度依赖完全可积性（逆散射、Riemann-Hilbert问题），能精确刻画解的渐近行为（如破裂前的振荡区）。本文则刻意避开可积结构（强调“does not rely on complete integrability”），在半直线非可积框架下建立鲁棒稳定性，适用性更广但精度较低（仅给出 $\mathcal{O}(\varepsilon^{3/2})$ 余项估计，而非点态渐近）。

  3. 与非线性Schrödinger边界层工作（如 [9,18]）的比较：NLS方程含二阶色散 $i\varepsilon \partial_t \psi + \varepsilon^2 \partial_{xx} \psi + |\psi|^2 \psi = 0$，其色散项是二阶且纯虚，能量结构（$L^2$ 范数守恒）与相位特性不同。本文的三阶实色散项导致能量不守恒且无自然守恒律，技术路径完全不同。

  4. 与KdV半直线适定性工作（如 [2,3]）的比较：文献[3]利用调和分析高级工具（Bourgain空间）建立了KdV在 $H^s$（$s \ge 0$）上的全局适定性，但未涉及零色散极限。本文的适定性证明（§4）采用更初等的抛物逼近与Aubin-Lions引理，虽仅得局部解（至 $T_0$），却为边界层分析提供了恰好的正则性框架（$H^3$ 解保证 $E_{\mathrm{inn}}$ 有界），体现了目标驱动的精巧设计。综上，本文填补了色散PDE边界层理论的关键空白，其方法论对其他高阶色散方程（如BBM、Rosenau）的边界层研究具有范式意义。

3. 非局部 Whitham 型方程中波破碎的尖锐阈值

A sharp threshold for wave breaking in nonlocal Whitham-type equations

作者: Yongki Lee

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.24839v1

摘要翻译: 本文建立了一类非线性非局部 Whitham 型方程的尖锐且扩展的波破碎判据，显著推广了 Seliger 所引入的经典阈值。尽管控制流体偏差空间极值的不等式组已被广泛研究，但建立一个尖锐条件仍是一项挑战。这一困难主要源于该系统的非合作性（non-cooperative nature），从而排除了标准比较原理（comparison principles）的应用。我们通过严格识别出精确的非线性阈值——该阈值将不等式组的解限制在其中——并分析其时间演化，从而严谨地克服了这一分析障碍。

Comment: 无

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

1. 研究对象和背景

本文研究一类非局部Whitham型浅水波方程的波破碎（wave breaking）现象，即解 $u(t,x)$ 本身在有限时间内保持有界，但其空间导数 $u_x(t,x)$ 在某点趋于负无穷：

$$ \lim_{t \to T^-} \inf_{x \in \mathbb{R}} u_x(t,x) = -\infty, \quad \text{而} \quad \sup_{(t,x) \in [0,T) \times \mathbb{R}} |u(t,x)| < \infty. $$

该现象是水波理论中核心的奇性形成机制，物理上对应激波前缘的陡峭化与翻卷，无法被KdV等弱非线性色散模型捕获。所研究方程为Whitham于1974年提出的非局部拟线性演化方程：

$$ \partial_t u + u \partial_x u + \int_{\mathbb{R}} K(x-\xi) u_\xi(t,\xi)\, d\xi = 0, \quad t>0,\ x\in\mathbb{R}, $$

其中核 $K$ 是相速度 $c(\kappa)$ 的傅里叶逆变换：$K(x) = \frac{1}{2\pi} \int_{\mathbb{R}} c(\kappa) e^{i\kappa x}\, d\kappa$；典型物理核如Whitham原始建议的奇异核

$$ K_0(x) = \int_{\mathbb{R}} \left(\frac{\tanh \xi}{\xi}\right)^{1/2} e^{i\xi x}\, d\xi, $$

其傅里叶乘子 $\widehat{K_0}(\kappa) = \left(\frac{\tanh \kappa}{\kappa}\right)^{1/2}$ 精确刻画了完全线性色散关系。函数 $u(t,x)$ 表示流体自由表面相对于静止位置的偏差。问题的关键数学对象是极值函数

$$ m_1(t) := \inf_{x\in\mathbb{R}} u_x(t,x), \quad m_2(t) := \sup_{x\in\mathbb{R}} u_x(t,x), $$

它们在经典分析中被证明（在适当正则性假设下，如 $u_0 \in H^2(\mathbb{R})$）是局部Lipschitz连续的，且几乎处处满足耦合微分不等式系统：

$$ \begin{cases} m_1'(t) \leq -m_1^2(t) + K(0)(m_2(t) - m_1(t)), \\ m_2'(t) \leq -m_2^2(t) + K(0)(m_2(t) - m_1(t)), \end{cases} \quad \text{a.e. } t>0, $$

其中 $K(0)>0$ 为核在原点的值（归一化后取 $K(0)=1$）。该系统源于对原方程关于 $x$ 求导并在极值点 $x=\xi_i(t)$ 处求值，其结构本质是非合作的（non-cooperative）：右端向量场 $F(m_1,m_2) = (-m_1^2 + m_2 - m_1,\ -m_2^2 + m_2 - m_1)$ 的雅可比矩阵的非对角元 $\partial F_2/\partial m_1 = -1 < 0$，违反Kamke比较定理所需的准单调性（quasi-monotonicity）条件，故标准比较原理完全失效。Seliger（1968）最早提出形式化判据 $m_1(0)+m_2(0) \leq -2K(0)$，其几何意义是初始点位于相平面中连接鞍点 $(-2K(0),2K(0))$ 的稳定流形下方的半平面内；Constantin–Escher（1998）首次将其严格化。然而，该判据远非最优：它仅利用了线性组合 $m_1+m_2$ 的单调性，而忽略了系统内在的非线性几何结构——特别是与鞍点相连的精确稳定流形（separatrix）。长期以来，如何从非合作不等式系统中提取出精确的、不可改进的初值域 $\Omega \subset \mathbb{R}_- \times \mathbb{R}_+$，使得 $ (m_1(0),m_2(0)) \in \Omega $ 当且仅当波破碎发生，构成该领域的根本性空白。本文旨在严格构造并证明这一精确阈值，其核心挑战在于：（1）不等式系统的轨迹无法用等式系统的解直接控制；（2）必须在缺乏全局单调性的前提下，建立一个“单向闸门”（one-way trapdoor），确保一旦初值落入 $\Omega$，轨迹便永不能逃逸；（3）需精确量化 $G(m_1,m_2)$ 沿轨迹的衰减速率以获得爆破时间上界。

2. 主要定理或主要结果

在核 $K$ 光滑、偶、在 $[0,\infty)$ 上单调递减且可积，并满足 $K(0)=1$ 的假设下，对初值 $u_0 \in H^2(\mathbb{R})$，定义

$$ m_1(0) = \inf_{x\in\mathbb{R}} u_0'(x), \quad m_2(0) = \sup_{x\in\mathbb{R}} u_0'(x). $$

令

$$ G(x,y) := x + y + \operatorname{sgn}(x-y+4) \sqrt{(x-y)\big((x-y) + 4 \ln(y-x) + 4 - 4 \ln 4\big)}. $$

定义区域

$$ \Omega := \{(x,y) \in \mathbb{R}_- \times \mathbb{R}_+ \mid G(x,y) < 0\}. $$

主要结果为：定理1.1——若 $(m_1(0), m_2(0)) \in \Omega$，则解 $u$ 在有限时间 $T$ 内发生波破碎，且

$$ \lim_{t \to T^-} \inf_{x\in\mathbb{R}} u_x(t,x) = -\infty, $$

其中爆破时间满足显式上界

$$ T \leq -\frac{2}{G(m_1(0), m_2(0))}. $$

此外，该阈值是sharp的：若 $(m_1(0), m_2(0))$ 位于曲线 $G(x,y)=0$ 上，则存在全局有界解（由辅助等式系统 $(2.1)$ 的相轨分析保证）；若初值在 $\Omega$ 外，则波破碎未必发生。作为推论，Seliger判据 $m_1(0)+m_2(0) \leq -2$ 被严格包含于 $\Omega$ 内，即 $\{x+y \leq -2\} \subsetneq \Omega$，且边界 $G(x,y)=0$ 恰好是辅助系统 $(2.1)$ 在鞍点 $(-2,2)$ 处的稳定流形。

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

本文的核心障碍是系统 $(1.6)$ 的严格非合作性导致所有标准比较原理（如Kamke定理、Müller比较定理）完全失效，无法将不等式系统的轨迹与等式系统的轨迹进行逐点控制。作者的真正创新点在于摒弃传统比较思路，转而构建一个基于相平面几何的“单向闸门”（one-way trapdoor）框架，其技术支柱由三个精密引理构成：（1）辅助等式系统的首次积分与稳定流形精确构造（Proposition 2.1）：通过变量代换 $u=x-y$, $v=x+y$ 将系统 $(2.1)$ 化为可积形式，导出首次积分

$$ \frac{(x+y)^2}{x-y} = x-y + 4 \ln|y-x| + C, $$

并代入鞍点 $(-2,2)$ 确定常数 $C=4-4\ln 4$，从而显式解出稳定流形 $G(x,y)=0$。该曲线是唯一能将相平面分割为爆破区 $\Omega$ 与全局存在区的精确分界。（2）不等式系统的单向不变性证明（Lemma 2.2）：定义垂直距离函数 $V(t) = y(t) - g(x(t))$，其中 $y=g(x)$ 是稳定流形。关键观察是，在接触点 $V(t)=0$ 处，有 $g'(x) = H/F$，代入 $V'$ 表达式后得到 $V'(t) = -b(t) + a(t)g'(x(t))$。由于 $g'<0$ 且 $a,b\geq0$，立即得 $V'(t) \leq 0$，从而证明轨迹一旦进入 $\Omega$，便永远无法穿越 $G=0$ 边界——这是克服非合作性的几何钥匙。（3）阈值函数的时间演化控制（Lemma 2.3）：对 $G(x(t),y(t))$ 求导，利用链式法则及 $G$ 的显式结构，严格证明

$$ \frac{d}{dt}G(x(t),y(t)) \leq -\frac{1}{2} G^2(x(t),y(t)). $$

此估计的取得依赖于两个精细计算：首先验证 $G_x = 1-S'(z), G_y = 1+S'(z)$ 在 $\Omega$ 内非负（即 $|S'(z)| \leq 1$），这等价于不等式 $z - (\ln(z/4)+2)^2 \geq 0$，而该式在 $z>4/e$ 时恒成立（最小值在 $z=4$ 处为0）；其次，通过代入 $v=G-S(z)$ 并利用恒等式

$$ -\frac{1}{2}S^2 - \frac{1}{2}z^2 + 2z + z S S' \equiv 0, $$

将导数上界简化为 $-\frac{1}{2}G^2 + G(S - zS')$，最终证明修正项 $G(S - zS') \leq 0$，从而获得决定性微分不等式。此不等式经Grönwall积分即得爆破时间上界。整个框架不依赖任何紧性或能量泛函，纯粹基于ODE相平面分析与不等式系统的几何控制，是处理非合作动力学的范式性突破。

4. 与之前工作的比较

本文与既有工作的推进体现在四个根本维度：（1）判据精度与sharpness：Seliger（1968）和Constantin–Escher（1998）的判据 $m_1(0)+m_2(0) \leq -2$ 是线性阈值，仅覆盖相平面中一个半平面；本文的 $G(m_1,m_2)<0$ 是非线性精确阈值，其边界 $G=0$ 是辅助系统稳定流形的解析表达，因此是不可改进的（sharp）。图1明确显示 $\Omega$ 严格大于蓝色半平面，且包含大量满足 $m_1+m_2 > -2$ 但仍会破碎的初值，例如 $(-3.5, 0.1) \in \Omega$ 但 $-3.5+0.1 = -3.4 > -2$。（2）方法论范式转移：此前所有工作（包括[1,2,3,6,7,11]）均依赖对 $m_1,m_2$ 的线性组合（如 $m_1+m_2$, $m_2-m_1$）构造Lyapunov型函数，其本质是降维投影，必然损失信息；本文首次将不等式系统提升至二维相平面，通过构造几何不变集 $\Omega$ 实现全状态控制，是方法论上的质变。（3）适用范围与代价：本文要求核 $K$ 可积且 $K(0)=1$，而Hur（2017）[4]等处理了非可积核（如 $K(x)\sim |x|^{-1/2}$），但其判据更粗糙；本文的sharpness是以牺牲核的奇异性为代价的。此外，本文结论针对 $H^2$ 初值，而[1]可推广至 $H^{3/2}$，但本文的阈值结构在更高正则性下依然sharp。（4）对后续工作的奠基性：文中明确指出，该几何框架可系统性推广至其他非合作非局部方程，如Fornberg-Whitham方程（见[2,3,11]）或更一般的非局部守恒律（见作者前作[6]）。与[7]（Ma–Liu–Qu）相比，后者虽也考虑Whitham型方程，但其判据仍为线性组合形式 $\alpha m_1 + \beta m_2 \leq C$，未触及精确流形；而本文证明了任何线性判据必为 $\Omega$ 的真子集。从提供内容中无法确认本文结果是否可推广至周期域或高维情形，作者亦未显式展开与[4]（Hur）在非可积核下的比较，仅说明‘for non-integrable kernels, we refer to [4] and references therein’。

4. 二维磁流体方程弱解的Sharp非唯一性

Sharp non-uniqueness of weak solutions to 2D magnetohydrodynamic equations

作者: Changxing Miao, Yao Nie, Weikui Ye

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.25097v1

摘要翻译: 本文证明：对任意给定的 $1\le p<\infty$，二维粘性、电阻性磁流体动力学（MHD）方程的弱解在函数空间 $L^2_t L^p(\mathbb{R}^2) \cap L^1_t W^{1,p}(\mathbb{R}^2)$ 中是非唯一的，从而揭示了Ladyzhenskaya–Prodi–Serrin条件在端点 $(2,\infty)$ 处的sharp性，且所构造的解位于Beale–Kato–Majda准则的临界边界上。据我们所知，这是首个关于二维粘性、电阻性MHD方程组的非唯一性结果。作为副产品，我们还得到了Navier–Stokes方程在 $L^2_t L^p$（其中 $1\le p<\infty$）中的非唯一性，以及具有大 $\mathrm{BMO}^{-1}$ 初值的MHD方程组的非唯一性。

Comment: 无

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

1. 研究对象和背景

本文研究二维粘性与电阻性不可压缩磁流体动力学（MHD）方程的弱解非唯一性问题，其控制方程为

$$ \begin{cases} \partial_t u - \Delta u + (u \cdot \nabla) u + \nabla p = B \cdot \nabla B, \\ \partial_t B - \Delta B + u \cdot \nabla B = B \cdot \nabla u, \\ \nabla \cdot u = \nabla \cdot B = 0, \end{cases} \quad (x,t) \in \mathbb{R}^2 \times [0,T], $$

其中 $ u : \mathbb{R}^2 \times [0,T] \to \mathbb{R}^2 $ 为速度场，$ B : \mathbb{R}^2 \times [0,T] \to \mathbb{R}^2 $ 为磁场，$ p $ 为标量压力。该系统是 Navier–Stokes 方程与麦克斯韦–欧姆定律在低马赫、高雷诺与高磁雷诺极限下的耦合模型，广泛应用于天体物理（如太阳耀斑、星际介质湍流）、聚变等离子体约束及液态金属流动建模。关键数学挑战在于：尽管二维情形通常具有更强的先验控制（如 $ L^2 $ 能量守恒、涡度 $ \omega = \nabla^\perp \cdot u $ 满足 $ L^1 \cap L^\infty $ 控制），但粘性与电阻项 $ -\Delta u $、$ -\Delta B $ 的存在使得标准能量法无法排除多个 Leray–Hopf 型弱解共存的可能性；而经典唯一性理论——Ladyzhenskaya–Prodi–Serrin（LPS）准则——指出：若弱解 $ (u,B) \in L^q_t L^p_x $ 满足

$$ \frac{2}{q} + \frac{2}{p} \le 1, \quad 2 \le p \le \infty, $$

则解唯一且正则（对 $ p > 2 $ 光滑，对 $ p = 2 $ 至少属于 $ C_t L^2 $）。该准则在端点 $ (q,p) = (2,\infty) $ 处恰好处于临界尺度：其齐次范数 $ \|u\|_{L^2_t L^\infty_x} $ 在自然缩放 $ u(t,x) \mapsto \lambda u(\lambda^2 t, \lambda x) $ 下不变，且与 Beale–Kato–Majda 准则 $ \int_0^T \|\nabla u(s)\|_{L^\infty} ds < \infty $ 的边界重合。然而，此前所有关于 MHD 方程非唯一性的构造均局限于三维（如 [34,44] 在 $ L^q_t L^\infty_x $、$ 1 \le q < 2 $ 中构造），其技术基础是凸积分（convex integration）框架，依赖于三维空间中丰富的振荡方向（即 $ \Lambda \subset \mathbb{S}^2 $ 具有非空内部的 $ \Lambda $-凸包），而在二维中，由于 $ \Lambda \subset \mathbb{S}^1 $ 导致 $ \Lambda $-凸包无内点，标准凸积分失效。此外，针对二维理想 MHD 的构造（如 [38]）利用有限传播速性质，但无法处理粘性与电阻项。因此，一个根本性空白长期存在：是否存在二维粘性–电阻 MHD 的两个不同 Leray–Hopf 弱解，共享同一初值，且同时落入 LPS 准则的临界补集 $ \frac{2}{q} + \frac{2}{p} > 1 $？特别是，能否实现 $ q = 2 $、任意 $ 1 \le p < \infty $（即 $ L^2_t L^p_x \cap L^1_t W^{1,p}_x $）这一此前完全未被触及的函数类？该问题不仅关乎 LPS 准则在二维的尖锐性验证，更触及二维 MHD 是否真比三维“更良态”的深层结构认知——若连 $ L^2_t L^p_x $（$ p < \infty $）都不足以保证唯一性，则二维系统的病态性远超传统预期。

2. 主要定理或主要结果

定理 1.2（尖锐非唯一性）：对任意给定 $ 1 \le p < \infty $，存在两个不同的 Leray–Hopf 弱解 $ (u,B) $ 与 $ (\tilde{u},\tilde{B}) $，满足相同的散度自由初值 $ (u^{\text{in}}, B^{\text{in}}) \in \mathcal{D}'(\mathbb{R}^2) $，且

$$ (u,B),\, (\tilde{u},\tilde{B}) \in L^2([0,T]; L^p(\mathbb{R}^2)) \cap L^1([0,T]; W^{1,p}(\mathbb{R}^2)). $$

更精确地，对 $ 2 < p < \infty $，存在参数 $ r = r(p) \in \left(0, \frac{1}{10} \min\left\{1 - \frac{2}{p}, \frac{2}{p}\right\}\right) $，整数列 $ \lambda_q = N_\Lambda^{a b^q} $（其中 $ N_\Lambda \in \mathbb{N} $ 由几何引理 2.3 确定，$ a,b \gg 1 $ 依赖于 $ T,p,\|f\|_{\dot{B}^0_{2,1} \cap \dot{B}^1_{2,1} \cap B^{-1}_{1,1}} $），以及截断尺度 $ \ell_q = \min\left\{ \lambda_q^{-50/(1-2/p)}, \lambda_q^{-25p} \right\} $，使得迭代构造的解序列 $ \{(u_n,B_n)\}_{n\ge1} $ 在 $ L^2_T L^p \cap L^1_T W^{1,p} $ 中收敛到极限解 $ (u,B) $，并满足误差估计

$$ \|(u_n,B_n) - (u_{n-1},B_{n-1})\|_{L^2_T L^p \cap L^1_T W^{1,p}} \lesssim C_0^{3/4} 2^{-n/p}, $$

其中常数 $ C_0 $ 为绝对常数。推论 1.3 进一步指出：（i）该结果亦适用于周期域 $ \mathbb{T}^2 $，且方法可直接推广至任意 $ d \ge 2 $ 维，从而确立 LPS 准则在端点 $ (q,p) = (2,\infty) $ 对所有 $ d \ge 2 $ 的尖锐性；（ii）对任意 $ d \ge 2 $ 及 $ 1 \le q < 2 $，存在两个不同弱解属于 $ L^q_T L^\infty $，统一覆盖了 [34,44] 的三维结果并扩展至新区域；（iii）所构造解还满足 $ (u,B) \in L^\infty(\mathbb{R}_+; \mathrm{BMO}^{-1}(\mathbb{R}^2)) $，从而证明大初值 $ \mathrm{BMO}^{-1} $ 数据下解的非唯一性，揭示了 [42] 中小初值 $ \mathrm{BMO}^{-1} $ 全局适定性结果的尖锐性；（iv）作为副产品，Navier–Stokes 方程在 $ L^2_T L^p(\mathbb{R}^d) $（$ 1 \le p < \infty $）中的非唯一性亦被证实，推广了 Cheskidov–Luo [9] 在 $ L^q_T L^\infty $（$ 1 \le q < 2 $）的结果。

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

本文的核心障碍在于：二维空间中标准凸积分因几何限制（$ \Lambda \subset \mathbb{S}^1 $ 导致 $ \Lambda $-凸包无内点）而彻底失效；同时，粘性项 $ -\Delta u $ 和电阻项 $ -\Delta B $ 不能像在欧拉方程中那样被忽略，必须被精确吸收而非视为低阶误差。作者提出一种全新的迭代机制，其创新点在于构造精确解序列而非近似解，且每一迭代步均严格满足原 MHD 方程（1.1），这与 [33,34,36,44] 中求解松弛系统的策略有本质区别。该机制由三部分增量组成：热主导 Mikado 流 $ w^{(h)}_n $、反级联流 $ w^{(i)}_n $ 与隐式修正流 $ (w^{(m)}_n, d^{(m)}_n) $。首先，$ w^{(h)}_n $ 定义为

$$ w^{(h)}_n \sim \lambda_n e^{-\lambda_n^2 t} \sum_{k \in \Lambda} a_k \left( \mathrm{Id} + \frac{\mathcal{R} w_{n-1}}{1000 \|\mathcal{R} w_{n-1}\|_{L^\infty}} \right) \eta_{n-1} \chi_{n-1} \phi(\lambda_n^{r} k \cdot x) e^{i \lambda_n k \cdot x} k^\perp + \text{l.o.t.}, $$

其中 $ \phi(\lambda_n^{r} k \cdot x) k^\perp $ 是二维 Mikado 流（由引理 2.2 提供），指数 $ r = r(p) $ 精确调控振荡频率以匹配 $ L^p $ 范数衰减；光滑截断 $ \chi_{n-1} $ 保证支集 $ \operatorname{spt} w^{(h)}_n \subset O_n $ 单调收缩，这是实现 $ L^2_t L^p $ 收敛的关键。其次，为处理 Mikado 流二次项 $ \operatorname{div}(w^{(h)}_n \otimes w^{(h)}_n) $ 产生的几何交叠误差，作者引入高度振荡流 $ e^{-\lambda_n^2 t} \sum_{k \in \Lambda} e^{i \lambda_n k \cdot x} k^\perp $，它同时是二维稳态 Euler 方程与自由热方程的精确解：前者允许将振荡部分吸收到压力项 $ \nabla p $ 中，后者则直接抵消粘性项 $ -\Delta w^{(h)}_n $ 的主部，避免将其降为误差。第三，反级联流 $ w^{(i)}_n = w^{(h)}_{n-1} e^{-2 \lambda_{n-1}^2 t} $ 的设计旨在精确抵消 $ \partial_t w^{(i)}_n $ 与 $ \operatorname{div}(w^{(h)}_n \otimes w^{(h)}_n) $ 中的 $ \lambda_n $-级振荡误差，仅留下低阶余项。最后，隐式部分 $ (w^{(m)}_n, d^{(m)}_n) $ 通过求解强制 MHD 系统获得：利用 $ w^{(h)}_n $、$ w^{(i)}_n $ 的指数衰减 $ e^{-\lambda_n^2 t} $，可证残差在 $ L^1_t \dot{B}^{-1+2/p}_{p,1} $ 中以双指数速率 $ \lambda_{n-1}^{-30} $ 衰减，而 $ (u_{n-2},B_{n-2}) $ 的 $ L^1_t C^1 $ 范数仅线性增长 $ O(n) $，故残差可视为小强迫项；再结合引理 2.6（带 $ v $-依赖指数的输运–扩散估计）与命题 2.7（振荡因子的 $ \dot{B}^{-1}_{\infty,\infty} $ 衰减），可证强制系统存在唯一小解 $ (w^{(m)}_n, d^{(m)}_n) $ 满足

$$ \|(w^{(m)}_n, d^{(m)}_n)\|_{\widetilde{L}^\infty_t \dot{B}^{-1+2/p}_{p,1} \cap L^1_t \dot{B}^{1+2/p}_{p,1}} \le \lambda_{n-1}^{-20}. $$

此框架彻底规避了二维凸积分的几何障碍，并将粘性/电阻项从‘误差源’转化为‘构造工具’，构成方法论的根本突破。

4. 与之前工作的比较

与既有工作相比，本文在多个维度实现实质性推进：（1）维度与模型适用性：此前所有 MHD 非唯一性结果均限于三维（[33,34,36,44]），依赖三维凸积分所需的丰富振荡方向；而本文首次攻克二维粘性–电阻 MHD，且方法可无缝推广至任意 $ d \ge 2 $，填补了该领域的根本性空白。相比之下，[38] 虽在二维理想 MHD 中构造了 Hölder 解，但其依赖有限传播速，不适用于含耗散项的系统。（2）函数空间尺度：本文目标空间 $ L^2_t L^p_x $（$ 1 \le p < \infty $）位于 LPS 准则 $ \frac{2}{q} + \frac{2}{p} > 1 $ 的临界线上，而此前结果均落在更‘远离’临界的区域：[34,44] 在 $ L^q_t L^\infty_x $（$ 1 \le q < 2 $）中满足 $ \frac{2}{q} + \frac{2}{\infty} = \frac{2}{q} > 1 $，但 $ p = \infty $ 无法控制 $ L^p $ 范数；[33] 的 $ \frac{1}{q} + \frac{2}{p} > 3 $、[36] 的 $ q = \infty, p = 2 $ 均显著偏离临界线。本文将临界尺度精确推进至 $ (q,p) = (2,\infty) $ 的‘对偶’区域 $ q = 2, p < \infty $，实现真正尖锐性验证。（3）技术路径：本文摒弃凸积分，开创性采用‘精确解迭代’框架，其核心是将粘性/电阻项主动融入构造（如用热方程解抵消 $ -\Delta w^{(h)}_n $），而非视作扰动。这与 [34,44] 将 $ -\Delta u $ 归入误差项、再用 Mikado 流平衡的策略截然不同，也避免了 [12] 在 $ \mathrm{BMO}^{-1} $ 构造中对周期域的依赖。（4）代价与局限：本文构造的解虽属于 $ L^2_t L^p_x $，但未达到 $ L^\infty_t L^2_x $（即标准 Leray–Hopf 类），其能量可能随时间增长；且解的正则性止步于 $ L^1_t W^{1,p}_x $，未涉及更高阶 Sobolev 或 Besov 空间。此外，构造依赖于特定的二维 Mikado 流与指数衰减设计，目前尚未推广至三维无粘性情形或更一般的各向异性耗散模型。从提供内容中无法确认该方法是否能处理 $ p = \infty $ 的端点情形，作者亦未显式展开与 [25] 计算机辅助证明的比较。

5. 关于满足全滑移边界条件并流过圆锥的轴对称Navier-Stokes流

On the axisymmetric Navier-Stokes flow passing a cone with the total-slip boundary condition

作者: Zijin Li, Xin Yang, Qi S. Zhang

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.25137v1

摘要翻译: 对于满足Navier全滑移边界条件、流过某类圆锥区域 $D$ 的不可压Navier-Stokes流，我们证明存在一个绝对常数 $C_* > 0$，使得若

Comment: 77 pages. All comments are welcome

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

1. 研究对象和背景

本文研究三维轴对称不可压Navier-Stokes方程（ASNS）在锥形外区域 $ D \subset \mathbb{R}^3 $ 上初边值问题的全局适定性，其中 $ D $ 是以 $ x_3 $-轴为对称轴、顶点在原点、半顶角为 $ \alpha \in (0, \pi/6] $ 的开圆锥体：

$$ D = \left\{ (r, \theta, x_3) \in \mathbb{R}^3 :\, 0 < r^2 + x_3^2 < 1,\; -r \tan\alpha < x_3 < r \tan\alpha,\; \theta \in [0,2\pi) \right\}. $$

在柱坐标系 $ (r,\theta,x_3) $ 下，速度场 $ \boldsymbol{v} = v_r e_r + v_\theta e_\theta + v_3 e_3 $ 满足轴对称性，即所有分量 $ v_r, v_\theta, v_3 $ 均与 $ \theta $ 无关；其控制方程为

$$ \begin{cases} \left( \Delta - \frac{1}{r^2} \right) v_r - (v_r \partial_r + v_3 \partial_{x_3}) v_r + \frac{v_\theta^2}{r} - \partial_r P - \partial_t v_r = 0, \\ \left( \Delta - \frac{1}{r^2} \right) v_\theta - (v_r \partial_r + v_3 \partial_{x_3}) v_\theta - \frac{v_\theta v_r}{r} - \partial_t v_\theta = 0, \\ \Delta v_3 - (v_r \partial_r + v_3 \partial_{x_3}) v_3 - \partial_{x_3} P - \partial_t v_3 = 0, \\ \frac{1}{r} \partial_r (r v_r) + \partial_{x_3} v_3 = 0, \end{cases} \tag{1.1} $$

其中 $ \Delta = \partial_r^2 + \frac{1}{r} \partial_r + \partial_{x_3}^2 $，$ P $ 为压力。该系统等价于经典三维不可压Navier-Stokes方程 $ \partial_t \boldsymbol{v} + (\boldsymbol{v} \cdot \nabla) \boldsymbol{v} - \Delta \boldsymbol{v} + \nabla P = 0 $，$ \nabla \cdot \boldsymbol{v} = 0 $ 在轴对称约束下的约化形式。

关键未知量是方位角分量 $ v_\theta $（即“swirl”），因其直接耦合涡量拉伸项并主导奇异性形成机制。定义 $ \Gamma := r v_\theta $，则 $ \Gamma $ 满足输运-扩散型方程

$$ \Delta \Gamma - \boldsymbol{b} \cdot \nabla \Gamma - \frac{2}{r} \partial_r \Gamma - \partial_t \Gamma = 0, \quad \boldsymbol{b} = v_r e_r + v_3 e_3, $$

该方程无源项且满足最大值原理，但仅在特定边界条件下可导出 $ L^\infty $ 控制。本文的核心困难在于：所施加的Navier全滑移（Navier total-slip, NTS）边界条件

$$ \boldsymbol{v} \cdot \boldsymbol{n} = 0, \quad (S \boldsymbol{v} \cdot \boldsymbol{n})_{\text{tan}} = 0 \quad \text{on } \partial D, $$

其中 $ S \boldsymbol{v} = \frac{1}{2}(\nabla \boldsymbol{v} + \nabla \boldsymbol{v}^\top) $ 为应变率张量，$ \boldsymbol{n} $ 为单位外法向，$ (\cdot)_{\text{tan}} $ 表示切向分量——该条件在锥顶 $ \rho = 0 $ 处退化，且不提供 $ v_\theta $ 的Robin型边界控制（如 $ \partial_n v_\theta + \beta v_\theta = 0 $），导致标准能量估计中出现无法吸收的坏符号边界项：例如在射线边界 $ R_1 \cup R_2 $ 上，由分部积分产生的 $ \int_{R_1 \cup R_2} F \partial_\rho F \, dS $ 项因缺乏正定性而无法被内域耗散 $ \|\nabla F\|_{L^2}^2 $ 控制。相比之下，此前工作 \cite{LPYZZZ24} 所采用的Navier-Hodge-Lions（NHL）边界条件 $ \omega \times \boldsymbol{n} = 0 $ 隐含 $ \partial_n v_\theta = 0 $，从而允许直接应用最大值原理证明 $ \|\Gamma(\cdot,t)\|_{L^\infty} \leq \|\Gamma_0\|_{L^\infty} $；而NTS条件下此性质失效，迫使作者发展全新技术框架。

此外，锥域几何本身引入双重奇异性：一是 $ r=0 $ 处的柱坐标奇异性（$ \partial_r $、$ 1/r $ 系数发散），二是锥面 $ \partial D $ 在原点处的非光滑性（曲率无界）。这使得经典Sobolev嵌入、迹定理及椭圆正则性均失效，传统方法（如基于 $ \Gamma $ 的De Giorgi迭代或Biot-Savart核的 $ L^p $ 估计）无法直接移植。已有理论中，\cite{Zhang2023} 在 $ \mathbb{R}^3 $ 上证明了 $ \|r v_\theta\|_{L^\infty} \leq C |\ln r|^{-2} $ 可排除奇点；\cite{LPYZZZ24} 在对称锥域上结合NHL条件与初值奇偶性假设 $ v_\theta(x_1,x_2,-x_3) = -v_\theta(x_1,x_2,x_3) $，得到小初值全局解；但二者均未处理NTS条件与非对称初值。本文旨在填补这一空白：在不依赖任何奇偶性假设、仅要求 $ \int_D r v_{0,\theta} \, dx = 0 $ 及 $ \sup_D r|v_{0,\theta}| \leq C_* $ 的前提下，建立NTS条件下ASNS的全局有界强解存在唯一性。

2. 主要定理或主要结果

在锥域 $ D $（定义见(1.4)，$ \alpha \in (0, \pi/6] $）上，对初值 $ \boldsymbol{v}_0 \in H^2(D) $ 满足 $ \nabla \cdot \boldsymbol{v}_0 = 0 $ 及NTS边界条件，若 (i) $ \int_D r v_{0,\theta}(x) \, dx = 0 $， (ii) 存在仅依赖于 $ D $ 上Sobolev嵌入常数的绝对常数 $ C_* > 0 $，使得 $ \sup_{x \in D} r |v_{0,\theta}(x)| \leq C_* $， 则对任意 $ T > 0 $，问题(1.11)存在唯一强解 $ (\boldsymbol{v}, P) $ 满足

$$ \|\boldsymbol{v}\|_{L^\infty_{t,x}(D \times [0,T])} + \|\boldsymbol{v}\|_{H^1_t L^2_x(D \times [0,T])} + \|\boldsymbol{v}\|_{L^2_t H^2_x(D \times [0,T])} + \|P\|_{L^2_t H^1_x(D \times [0,T])} \leq C, $$

其中 $ C $ 仅依赖于 $ \|\boldsymbol{v}_0\|_{C^2(\overline{D})} $；且解保持零平均性质：

$$ \int_D r v_\theta(x,t) \, dx = 0, \quad \forall t \in [0,T]. $$

能量不等式成立：

$$ \int_D |\boldsymbol{v}(x,T)|^2 \, dx + \frac{3}{4} \int_0^T \int_D |\nabla \boldsymbol{v}(x,t)|^2 \, dx dt \leq \int_D |\boldsymbol{v}_0(x)|^2 \, dx. $$

进一步，该解属于能量空间 $ \mathcal{E} = L^2_t H^1_x \cap L^\infty_t L^2_x $，且 $ \boldsymbol{v} $ 全局有界：$ \|\boldsymbol{v}\|_{L^\infty_{t,x}} \leq C_*' $，其中 $ C_*' $ 由 $ C_* $ 和 $ \|\boldsymbol{v}_0\|_{C^2} $ 控制。定理1.3的关键突破在于：完全取消了前作 \cite{LPYZZZ24} 中对初值速度场的奇偶性限制，并将适用域推广至非对称锥域 $ \widetilde{D} = \{ (r,\theta,x_3) : -r \tan \alpha_1 < x_3 < r \tan \alpha_2 \} $（$ 0 \leq \alpha_1, \alpha_2 \leq \pi/6 $，不同时为零），包括珊瑚状外域 $ \widetilde{D}_0 = \{ 0 < x_3 < r \tan \alpha_2 \} $。

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

证明的核心障碍是：NTS边界条件在锥顶处退化，导致标准能量估计中出现无法控制的坏符号边界项；同时，$ \Gamma = r v_\theta $ 不再满足最大值原理，需另寻途径控制其 $ L^\infty $ 范数。作者构建了一个全新的三元“好未知量”系统 $ (\mathcal{K}, \mathcal{F}, \mathcal{O}) $，其定义为

$$ \mathcal{K} \coloneqq \frac{\sin \phi}{\rho^2} \partial_\phi \left( \frac{v_\theta}{\sin \phi} \right), \quad \mathcal{F} \coloneqq -\partial_\rho \left( \frac{v_\theta}{\rho} \right), \quad \mathcal{O} \coloneqq \frac{1}{\rho \sin \phi} \left( \omega_\theta - \frac{2 v_\phi \eta(\rho)}{\rho} \right), $$

其中 $ \eta $ 是截断函数（$ \eta(\rho) = 0 $ 当 $ \rho \leq 1/3 $，$ \eta(\rho) = 1 $ 当 $ \rho \geq 2/3 $）。$ \mathcal{K} $ 和 $ \mathcal{F} $ 的创新在于其边界行为：在射线边界 $ R_1 \cup R_2 $ 上，$ \mathcal{K} $ 满足齐次Dirichlet条件（$ \mathcal{K}|_{R_1 \cup R_2} = 0 $），而 $ \mathcal{F} $ 满足Robin型条件（$ \partial_\rho \mathcal{F} \sim \text{lower order} $），从而彻底消除了分部积分中坏符号边界积分；在弧形边界 $ A $ 上，二者均自然满足Robin条件。$ \mathcal{O} $ 则是连接经典未知量 $ \Omega = \omega_\theta / (\rho \sin \phi) $ 与NTS适配量 $ \tilde{\Omega} = \partial_\rho(v_\phi / (\rho \sin \phi)) - \partial_\phi(v_\rho / \rho^2) $ 的桥梁，通过截断函数 $ \eta $ 实现局部切换，确保其在全边界上具备可控的边界通量。

为闭合 $ (\mathcal{K}, \mathcal{F}, \mathcal{O}) $ 的能量估计，必须先控制压力 $ P $。作者提出一个非标准的边界值问题：在近似域 $ D_m $ 上，$ P $ 满足 $ \Delta P = \nabla \cdot (\boldsymbol{v} \cdot \nabla \boldsymbol{v}) $，且边界条件由NTS推导出的动量方程切向投影给出（引理3.1）。由此建立关键椭圆估计：$ \|\nabla P\|_{L^2(D_m)} \lesssim \|\boldsymbol{v}\|_{H^1}^2 + \|\mathcal{O}\|_{L^2}^2 $，该估计依赖于对非线性项 $ \|v_\phi \partial_\rho v_\rho\|_{L^2} $ 的精细处理（引理9.4），其证明使用了各向异性插值不等式而非标准Gagliardo-Nirenberg。

对 $ \Gamma $ 的 $ L^\infty $ 控制，作者摒弃了在NTS下失效的Moser迭代，转而采用De Giorgi迭代：定义 $ \Gamma_k = (\Gamma - k)_+ $，通过构造层级能量不等式 $ \|\Gamma_{k_{j+1}}\|_{L^2}^2 \leq C 2^{-j} \|\Gamma_{k_j}\|_{L^2}^2 $ 并利用 $ \int_D r v_\theta \, dx = 0 $ 导出的零平均性质，最终获得 $ \|\Gamma\|_{L^\infty} \leq C_* $。该迭代成功的关键在于：NTS下坏边界项的增长率被精确计算为 $ \lesssim k \|\nabla \Gamma_k\|_{L^2} $，可被耗散项吸收。

最根本的工具是新证明的各向异性Hardy不等式（推论5.1）：对任意满足 $ \int_D r v_\theta \, dx = 0 $ 的函数，

$$ \int_{D_m} \left| \frac{1}{\rho} \frac{v_\theta}{\rho} \right|^2 dx \leq \frac{8}{\sqrt{3}} \int_{D_m} \left| \partial_\rho \left( \frac{v_\theta}{\rho} \right) \right|^2 dx + \frac{1}{5\sqrt{3}} \int_{D_m} \left| \frac{1}{\rho} \partial_\phi \left( \frac{v_\theta}{\rho} \right) \right|^2 dx. $$

该不等式右侧系数 $ 8/\sqrt{3} \approx 4.618 $ 远小于标准三维Hardy常数4？此处需修正：实际是远小于导致估计破裂的临界值（原文指出若为4则方法失效），其微小性源于零平均假设，使 $ v_\theta/\rho $ 在 $ \phi $-方向具有Poincaré型控制，从而将 $ L^2 $ 范数控制为一阶导数的加权组合，为闭合 $ (\mathcal{K}, \mathcal{F}, \mathcal{O}) $ 的能量估计提供了不可或缺的尺度关系。

4. 与之前工作的比较

与前人工作相比，本文在四个维度实现实质性推进：

  1. 边界条件维度：此前唯一处理锥域ASNS全局解的工作是 \cite{LPYZZZ24}，其采用NHL边界条件 $ \omega \times \boldsymbol{n} = 0 $，该条件隐含 $ \partial_n v_\theta = 0 $，从而保证 $ \Gamma $ 满足最大值原理，可直接导出 $ L^\infty $ 控制。而本文首次攻克更物理自然但数学上更棘手的NTS条件 $ (S \boldsymbol{v} \cdot \boldsymbol{n})_{\text{tan}} = 0 $，该条件不控制 $ v_\theta $ 的法向导数，导致经典方法全面失效。代价是必须引入截断函数 $ \eta $ 构造混合未知量 $ \mathcal{O} $，并额外建立压力的椭圆估计以处理新增的非线性耦合项。

  2. 初值假设维度：\cite{LPYZZZ24} 要求 $ v_\theta $ 关于 $ x_3 $ 为奇函数（即 $ v_\theta(r,-x_3) = -v_\theta(r,x_3) $），该对称性用于在 $ \phi $-方向应用Poincaré不等式，导出关键的Hardy型估计 $ \|v_\theta/\rho^2\|_{L^2} \lesssim \|\nabla(v_\theta/\rho)\|_{L^2} $。本文以更弱的零平均条件 $ \int_D r v_{0,\theta} \, dx = 0 $ 替代，该条件不依赖任何空间对称性，从而允许初值在锥域两侧不对称分布，极大拓展了适用场景（如珊瑚状域 $ \widetilde{D}_0 $）。从技术上看，零平均条件导出的各向异性Hardy不等式系数 $ 8/\sqrt{3} $ 虽略大于奇函数假设下的 $ 2/19 $，但作者通过精巧设计 $ \mathcal{K}, \mathcal{F} $ 的边界行为，成功补偿了这一代价。

  3. 几何适用性维度：\cite{LPYZZZ24} 的结果严格限定于对称锥域 $ D $（$ \alpha_1 = \alpha_2 $），其证明中多次使用 $ \phi $-方向的对称性。本文证明全程未使用对称性，故结论自动推广至非对称锥域 $ \widetilde{D} = \{ -r \tan \alpha_1 < x_3 < r \tan \alpha_2 \} $，且 $ \alpha_1, \alpha_2 $ 可独立选取（只要 $ \leq \pi/6 $）。这是首个在非对称几何上建立ASNS全局有界解的工作。

  4. 技术框架维度：\cite{LPYZZZ24} 基于经典未知量 $ (K,F,\Omega) = (\omega_\rho/\rho, \omega_\phi/\rho, \omega_\theta/(\rho \sin \phi)) $，其边界行为在NHL下良好；本文则完全重构未知量体系，$ \mathcal{K}, \mathcal{F} $ 直接以 $ v_\theta $ 为基本变量构造，确保在NTS下具备理想边界性质；$ \mathcal{O} $ 则是首个为NTS定制的涡量相关未知量。此外，De Giorgi迭代替代Moser迭代，避免了NTS下迭代幂次增长导致的系数爆炸问题。从提供内容中无法确认本文结果是否改进了时间尺度（如从有限时间到全局）或能量类（如从Leray-Hopf到强解）的严格性，因 \cite{LPYZZZ24} 同样构造了强解；但本文在更弱假设下达到同等正则性，实为实质强化。

6. 共振微结构作为狄拉克型驱动器用于声波控制

Resonant Microstructures as Dirac-type Actuators for Acoustic Wave Control

作者: Arpan Mukherjee, Mourad Sini

类别: math.AP, math.OC

发布日期: 2026-05-25

arXiv ID: 2605.25978v1

摘要翻译: 我们研究通过有限个共振扰动（建模为声学亚波长气泡）产生的有效点源，对声波方程进行内域控制。在抽象层面，先将全空间动力学局部化至一个大的辅助观测区域，再考虑带有限个位于预设内点的点驱动器的狄利克雷谱形式波方程。限制在狄利克雷本征频率的一个有限谱带内，我们在关联耦合矩阵满足自然满秩条件下，证明了该谱子空间上的任意轨迹均可被实现，并给出了控制代价的定量估计，其依赖于谱带几何结构与驱动器位置。随后我们证明：这些理想点驱动器可由小尺寸、高对比度气泡簇实现。利用时域渐近展开，散射场被表示为一系列延迟单极子的叠加，其幅度满足一个有限维延迟双曲系统。在拉普拉斯域中，该系统诱导出一个传输算子，其极点结构编码了具有集体衰减效应的Minnaert共振。我们证明：当远离共振时，相应的驱动映射是病态的；而在一种簇层级换能器可达性条件（将入射场与主导簇通道相联系）下，该驱动映射在合适的Minnaert谱带上存在有界右逆。因此，当气泡尺寸 $\varepsilon \to 0$ 时，可实现波场的谱跟踪，误差为 $\mathcal{O}(\varepsilon^\gamma)$。

Comment: 无

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

1. 研究对象和背景

本文研究三维全空间中声波方程的内点控制问题，其物理背景是利用亚波长尺度、高对比度气泡簇（resonant subwavelength bubbles）作为可编程的声学执行器。核心控制对象是齐次背景介质 $\mathbb{R}^3$ 中满足声波方程的总压力场 $u(x,t)$，其动力学由如下传输问题刻画：

$$ \begin{cases} \kappa^{-1}(x) \partial_{tt} u - \operatorname{div}\big( \rho^{-1}(x) \nabla u \big) = \sum_{m=1}^{M_{\mathrm{tr}}} \lambda_m(t) \delta_{x^{\mathrm{tr}}_m}(x), & (x,t) \in (\mathbb{R}^3 \setminus \partial D) \times (0,T), \\ u|_+ = u|_-, & \text{on } \partial D \times (0,T), \\ \rho_c^{-1} \partial_\nu u|_+ = \rho_b^{-1} \partial_\nu u|_-, & \text{on } \partial D \times (0,T), \\ u(\cdot,0) = 0,\quad \partial_t u(\cdot,0) = 0, & \text{in } \mathbb{R}^3, \end{cases} $$

其中 $D = \bigcup_{i=1}^M D_i$ 是 $M$ 个半径为 $\varepsilon \ll 1$ 的球形气泡的并集，其材料参数按 Assumption 1 缩放：$[\kappa(x), \rho(x)] = [\kappa_c, \rho_c]$（背景），$[\kappa_b^i \varepsilon^2, \rho_b^i \varepsilon^2]$（气泡内部），从而保证局部波速 $c = \sqrt{\kappa/\rho} = O(1)$ 一致有界。未知量包括外加激励 $\lambda = (\lambda_1,\dots,\lambda_{M_{\mathrm{tr}}}) \in U_{\mathrm{tr}} := \{ \lambda \in H^2(0,T;\mathbb{R}^{M_{\mathrm{tr}}}) : \lambda(0)=\dot{\lambda}(0)=0 \}$，以及由多散射生成的微观气泡振幅向量 $q^\varepsilon_{\mathrm{mic}}(t) = (q^\varepsilon_1(t),\dots,q^\varepsilon_M(t))$，它满足一个具有时滞项的有限维代数-微分系统（见(1.1)与(4.4)）。关键建模步骤是将微观响应 $q^\varepsilon_{\mathrm{mic}}$ 聚合为 $N$ 个宏观簇级源强度 $Q^\varepsilon_\alpha(t) = \sum_{i\in I_\alpha} q^\varepsilon_i(t)$，进而定义有效散射场

$$ p^{\mathrm{eff},\varepsilon}(x,t) = \sum_{\alpha=1}^N \frac{1}{4\pi |x - y_\alpha|} Q^\varepsilon_\alpha\big(t - c_0^{-1}|x - y_\alpha|\big), $$

其中 $y_\alpha$ 是第 $\alpha$ 个气泡簇的宏观中心。该场在大而光滑的观测域 $\Omega \subset \mathbb{R}^3$（满足 $\operatorname{dist}(D,\partial\Omega) > c_0 T$）内被截断，并通过齐次 Dirichlet 边界条件引入谱框架。由此导出的理想化控制模型为

$$ \left( \frac{1}{c_0^2} \partial_{tt} - \Delta \right) p(t,x) = \sum_{\alpha=1}^N q_\alpha(t) \delta_{y_\alpha}(x), \quad (t,x) \in (0,T) \times \Omega, $$

其中 $p|_{\partial\Omega} = 0$，初始条件为零。由于 $\delta_{y_\alpha} \notin H^{-1}(\Omega)$，该方程仅在投影到有限维 Dirichlet 谱子空间 $H_J = \operatorname{span}\{\phi_k : \omega_k = c_0 \sqrt{\lambda_k} \in J\}$ 后良定；投影后得到耦合振子系统：

$$ \ddot{p}_k(t) + \omega_k^2 p_k(t) = c_0^2 \sum_{\alpha=1}^N q_\alpha(t) \phi_k(y_\alpha), \quad \forall k \in K_J. $$

本文试图填补的核心空白在于：如何在严格数学框架下，将“理想狄拉克点源”的可控性理论（即对任意给定轨迹 $p_r(t) \in H_J$ 实现精确跟踪）与真实物理实现（即通过亚波长气泡簇的共振响应近似生成该点源）统一起来。已有工作如 [6,14,17,22,24] 主要处理边界或体分布控制，依赖几何控制条件（GCC）；而点源控制 [22,23] 多停留在抽象分布意义，缺乏物理可实现性分析；气泡渐近理论 [2,3,4,5,7,15,18,19] 则聚焦于单频散射或稳态响应，未建立与时间域轨迹跟踪、谱带可控性及控制代价的定量联系。本文首次构建了从微观多散射动力学到宏观谱控目标的双重尺度分析链：一方面需刻画气泡簇在 Laplace 域的传递矩阵 $H_b(s)$ 的极点结构（特别是 Minnaert 共振极点 $s^*_{\alpha,1} = -\eta_{\alpha,1} + i \omega_{\alpha,1}(\varepsilon)$ 的分离性与残差非退化性），另一方面需证明该物理系统诱导的“执行器映射” $T^\varepsilon: U_{\mathrm{tr}} \to L^2(0,T;\mathbb{R}^N)$ 在适当 Minnaert 频带 $I_M = \bigcup_{\alpha=1}^N I_\alpha(\varepsilon)$ 上具有有界右逆，从而将抽象可控性转化为可工程实现的声学控制协议。

2. 主要定理或主要结果

本文主要结果分为两个层次：理想点源可控性与气泡物理实现性。首先，在理想层面，设 $\Omega \subset \mathbb{R}^3$ 为 $C^{1,1}$ 光滑有界区域，$\{\phi_k, \lambda_k\}_{k\ge1}$ 为其 Dirichlet 特征对，$\omega_k = c_0 \sqrt{\lambda_k}$。对任意有限频带 $J \subset (0,\infty)$，记 $K_J = \{k : \omega_k \in J\}$，$H_J = \operatorname{span}\{\phi_k : k \in K_J\}$，$P_J$ 为其正交投影。令 $y = (y_1,\dots,y_N) \in \Omega^N$ 为 $N$ 个宏观簇中心，定义空间耦合矩阵 $C_M(y) = \big( \phi_k(y_\alpha) \big)_{k \in K_J,\, \alpha=1}^N \in \mathbb{R}^{|K_J| \times N}$。定理 3.2（隐含于 Section 3.2）断言：若 $\operatorname{rank} C_M(y) = |K_J|$（即满秩条件），则对任意目标轨迹 $p_r \in C^2([0,T]; H_J)$ 满足 $p_r(0) = \dot{p}_r(0) = 0$，存在控制 $q \in L^2(0,T;\mathbb{R}^N)$ 使得对应解 $p$ 满足 $P_J p(t) = p_r(t)$ 对所有 $t \in [0,T]$ 成立；且控制代价满足定量估计 $\|q\|_{L^2(0,T;\mathbb{R}^N)} \le C_J \|\ddot{p}_r + \Lambda_J p_r\|_{L^2(0,T;H_J)}$，其中 $\Lambda_J = \operatorname{diag}(\omega_k^2)_{k \in K_J}$，常数 $C_J$ 依赖于 $C_M(y)$ 的最小奇异值 $\sigma_{\min}(C_M(y))$ 及频带几何（如 $\max J - \min J$）。其次，在物理实现层面，设气泡簇满足 Assumption 1–3：材料缩放、双尺度几何（簇内间距 $\sim \varepsilon^p$，簇间间距 $O(1)$）、宏观失谐（$\omega_{M,\alpha} \neq \omega_{M,\beta}$）及换能器可及性（$\inf_{\omega \in K_{\mathrm{tr}}} \sigma_{\min}(\widetilde{G}_{\mathrm{tr}}(i\omega)) \ge c_{\mathrm{tr}} > 0$）。令 $I_\alpha(\varepsilon) = [\hat{\omega}_{\tau(\alpha)} - \delta(\varepsilon), \hat{\omega}_{\tau(\alpha)} + \delta(\varepsilon)]$，其中 $\delta(\varepsilon) \in (0, \frac{g_0}{2} \varepsilon^{1-p})$，$\hat{\omega}_\ell$ 为选定的 Dirichlet 频率，$\tau: \{1,\dots,N\} \twoheadrightarrow \{1,\dots,L\}$ 为满射。定义复合 Minnaert 带 $I_M = \bigcup_{\alpha=1}^N I_\alpha(\varepsilon)$。主定理（Theorem 6.3 与 Corollary 7.1 综合）断言：存在 $\varepsilon_0 > 0$ 及 $\gamma > 0$（显式依赖于 $p$ 与阻尼率 $\eta_{\alpha,1} \sim \varepsilon$），使得对所有 $\varepsilon \in (0,\varepsilon_0]$，存在有界线性算子 $\mathcal{R}^\varepsilon: L^2(0,T; H_J) \to U_{\mathrm{tr}}$ 满足

$$ \left\| P_J u^\varepsilon(\cdot, t) - p_r(t) \right\|_{L^2(\Omega)} \le C \varepsilon^\gamma \|p_r\|_{C^2([0,T]; H_J)}, \quad \forall t \in [0,T], $$

其中 $u^\varepsilon$ 是由 $\lambda = \mathcal{R}^\varepsilon(p_r)$ 激发的真实气泡系统解，且 $\|\mathcal{R}^\varepsilon\|_{\mathcal{L}(L^2(0,T;H_J), U_{\mathrm{tr}})} \le C' \varepsilon^{-1}$（控制增益发散由共振带宽 $\sim \varepsilon^{1-p}$ 所致）。误差阶 $\gamma$ 源于两层渐近：一是簇级展开余项 $R^\varepsilon_{\mathrm{cl}} = O(\varepsilon^{1-p})$（来自(1.2)），二是 Laplace 域极点逼近误差（Proposition 2.1(iii) 的残差展开截断）。此外，Proposition 2.1(ii) 给出极点隔离尺度：$\operatorname{dist}(s^*_{\alpha,1}, P \setminus \{s^*_{\alpha,1}\}) \ge g_0 \varepsilon^{1-p}$，且跨簇极点间距 $|s^*_{M,\alpha}(\varepsilon) - s^*_{M,\beta}(\varepsilon)| \ge \Delta_\omega > 0$，其中 $\Delta_\omega = \frac{1}{2} \min_{\alpha \neq \beta} |\omega_{M,\alpha} - \omega_{M,\beta}|$。

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

证明的核心障碍在于跨越三个尺度的耦合：微观气泡动力学（$\varepsilon$-尺度延迟系统）、宏观簇级输出（$O(1)$ 空间尺度但含 $\varepsilon$-依赖时滞）、以及目标谱子空间（Dirichlet 模态的离散频谱）。首要技术难点是建立 Laplace 域传递矩阵 $H_b(s)$ 的精细极点分析，因为其谱性质直接决定物理执行器的可控性。作者并未采用标准的 Krein–Rutman 或抽象算子扰动理论，而是发展了一套针对高对比度多体共振系统的“双尺度极点分离技术”。具体而言，在 Proposition 2.1 的证明（Section 5）中，关键工具是将动态矩阵 $K_{\mathrm{dyn}}(s) = I + D(s)^{-1} s^2 Q(s)$ 的行列式零点问题分解为块对角主导项与交互项：$D(s) = \operatorname{diag}(\omega_{M,i}^{-2} s^2 + 1)$ 是解耦的 Minnaert 近似，而 $Q(s)$ 的非对角元 $Q_{ij}(s) = \varepsilon \frac{C_j e^{-s \tau_{ij}}}{4\pi d_{ij}}$（$i \neq j$）携带簇内时滞 $\tau_{ij} = c_0^{-1} |z_i - z_j|$ 与几何衰减。利用 Assumption 2(ii) 的双尺度间距，作者证明 $Q(s)$ 的范数在 $s$ 接近 $i\omega_{M,\alpha}$ 时为 $O(\varepsilon^{1-p})$，从而将 $\det(K_{\mathrm{dyn}}(s)) = 0$ 的求解转化为对每个 $\alpha$ 的标量方程 $\omega_{M,\alpha}^{-2} s^2 + 1 + O(\varepsilon^{1-p}) = 0$，由此导出极点 $s^*_{\alpha,1} = -\eta_{\alpha,1} + i\omega_{\alpha,1}(\varepsilon)$ 及其 $\sim \varepsilon^{1-p}$ 的隔离间隙。第二个创新点是“簇级投影与换能器可及性”的协同分析。传统方法难以处理外加激励 $\lambda$ 到簇输出 $Q^\varepsilon_\alpha$ 的映射，因为 $G_{\mathrm{tr}}(s) = (G_s(z_i, x^{\mathrm{tr}}_m))$ 是 $M \times M_{\mathrm{tr}}$ 矩阵，而 $B_{\mathrm{out}}$ 是 $N \times M$ 投影。作者引入宏观迹矩阵 $\widetilde{G}_{\mathrm{tr}}(s) = (G_s(y_\alpha, x^{\mathrm{tr}}_m))$ 并利用其远场渐近 $G_{i\omega}(y_\alpha, x^{\mathrm{tr}}_m) = c_m(\omega) e^{i(\omega/c_0) \theta_m \cdot y_\alpha} + O(R_m^{-1})$，结合 Rellich 引理与唯一延拓性，证明 $\operatorname{rank}(\widetilde{G}_{\mathrm{tr}}) = N$ 且 $\sigma_{\min}(\widetilde{G}_{\mathrm{tr}}(i\omega)) \ge c_{\mathrm{tr}} > 0$ 在紧集 $K_{\mathrm{tr}}$ 上一致成立。再通过 $G_{\mathrm{tr}}(s) = \widetilde{G}_{\mathrm{tr}}(s) + O(r_\varepsilon)$（$r_\varepsilon = \max_{\alpha,i\in I_\alpha} |z_i - y_\alpha|$）及极点处的残差展开（(2.19)），最终获得 $\sigma_{\min}(H^\varepsilon_{\mathrm{ext}}(i\omega)) \ge c_{\mathrm{ext}} \varepsilon / \eta_{\max} \ge c > 0$（(2.20)），这是构造有界右逆 $\mathcal{R}^\varepsilon$ 的基石。第三个关键技术是“谱带滤波与误差传播控制”。为避免非共振频段的病态放大，作者在 Section 6.1 明确定义带通滤波算子 $\Pi_{I_M}: L^2(0,T;\mathbb{R}^N) \to L^2(0,T;\mathbb{R}^N)$，其 Fourier 乘子为 $\mathbf{1}_{I_M}(\omega)$。结合能量估计（Section 7），将时间域误差 $\|P_J u^\varepsilon - p_r\|_{L^2}$ 控制为 $\|\Pi_{I_M} Q^\varepsilon_{\mathrm{cl}} - q_{\mathrm{ideal}}\|_{L^2} + \|R^\varepsilon_{\mathrm{cl}}\|_{L^2} + \|\text{non-resonant modes}\|$，其中前两项由极点分析与滤波保证为 $O(\varepsilon^\gamma)$，最后一项由非共振带的指数衰减 $e^{-\sigma_0 t}$ 控制。整个框架回避了对无限维系统做全局可控性分析，转而聚焦于可工程实现的有限谱带，体现了控制论与多尺度 PDE 分析的深度融合。

4. 与之前工作的比较

本文与既有工作的推进体现在四个维度。第一，与经典边界/体分布控制理论（如 Bardos–Lebeau–Rauch GCC [6] 或 Lions’ Hilbert uniqueness method [14]）相比，本文将控制自由度从无限维函数空间压缩至有限维 $N$ 个点源，并证明其在物理上可由亚波长结构实现，从而绕开了 GCC 对观测区域几何的苛刻要求（如非俘获性、测地线覆盖等），代价是控制目标局限于有限谱带 $H_J$ 且代价常数 $C_J$ 随 $\sigma_{\min}(C_M(y))^{-1}$ 发散。第二，相较于纯抽象点源控制研究（如 [22,23]），本文首次建立了点源“理想可控性”与“物理可实现性”的定量桥梁。以往工作仅证明 $\delta_{y_\alpha}$ 在 $H^{-1}(\Omega)$ 中的稠密性或弱可控性，未涉及实现机制；而本文通过 Minnaert 共振极点的精细渐近（$\operatorname{dist}(s^*_{\alpha,1}, \text{其他极点}) \sim \varepsilon^{1-p}$）与换能器可及性条件（(2.17)），给出了物理执行器映射 $T^\varepsilon$ 具有有界右逆的充分必要条件，且误差阶 $O(\varepsilon^\gamma)$ 显式依赖于气泡尺寸与几何参数。第三，相对于气泡散射的渐近理论（如 [2,3,4,5,7,15,18,19]），本文将静态/单频散射结果升级为时间域轨迹跟踪控制。例如，[18,19] 证明了散射场的单极子展开 (1.1)，但未分析其作为执行器的可控性；本文则进一步将 $q^\varepsilon_i$ 视为受控输入，通过延迟系统 (4.4) 与 Laplace 域极点结构（Proposition 2.1），揭示了共振带宽 $\sim \varepsilon^{1-p}$ 与控制带宽 $\delta(\varepsilon)$ 的匹配关系，这是前人工作中完全缺失的控制论视角。第四，与近期基于超材料的波前整形工作（如 [20] 中的等离子体热生成）相比，本文的“狄拉克型执行器”具有明确的物理起源（Minnaert 共振）与严格的多尺度分析框架，而非经验性等效参数拟合。其局限在于：目前仅适用于三维全空间中的稀疏气泡簇（Assumption 2 要求簇间距离 $O(1)$），尚未推广至高密度随机介质或二维情形；且控制代价 $\|\mathcal{R}^\varepsilon\| \sim \varepsilon^{-1}$ 表明小尺寸虽提升空间分辨率，却以更高驱动功率为代价，这与实际换能器的功率饱和限制构成潜在冲突。此外，文中未显式展开与 Kato 解析扰动理论 [12] 的比较，仅引用其用于证明 Proposition 2.3 中特征函数行列式的零测集性质；亦未将 Perron–Frobenius 谱（关键词提及）在证明中实质性使用，从提供内容中无法确认其具体作用。

7. Hardy-Sobolev 抛物方程在 $\mathbb{R}^N$ 中的整体存在性与爆破

Global existence and blow-up for the Hardy-Sobolev parabolic equation in RN

作者: Fei Fang, Zhong Tan

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.24802v1

摘要翻译: 本文应用自相似变换，将含 Hardy 项的抛物方程 \begin{equation*} \begin{cases}u_t-\Delta u-\mu\frac{u}{|x|^2}=|u|^{2^-2} u & \text { in } \mathbb{R}^N \times(0, T), \ u(x, 0)=u_0(x) & \text { in } \mathbb{R}^N , \end{cases} \end{equation} 转化为如下抛物方程 \begin{equation*} \begin{cases} v_s-\Delta v-\frac{1}{2} y \cdot \nabla v=\beta v+\frac{\mu v}{|y|^2}+|v|^{2^-2} v &\text { in } \mathbb{R}^N \times(0, S), \ \left.v\right|_{s=0}=v_0 & \text { in } \mathbb{R}^N, \end{cases} \end{equation} 其中 $N \geqslant 3$, $\mu\in [0,(N-2)^2 /8]$ 且 $2^{\ast}=2N /(N-2)$。对此方程，我们建立了一个加权 Hardy 不等式。此外，借助修正的势井方法与 Palais-Smale 序列分析，我们研究了该抛物方程解的长时间行为与有限时间爆破性质。

Comment: 无

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

1. 研究对象和背景

本文研究如下带Hardy位势和临界Sobolev非线性的半线性抛物方程初值问题：

$$ \begin{cases} u_t - \Delta u - \mu \dfrac{u}{|x|^2} = |u|^{2^*-2} u & \text{in } \mathbb{R}^N \times (0, T), \\ u(x, 0) = u_0(x) & \text{in } \mathbb{R}^N, \end{cases} \tag{1.1} $$

其中 $N \geqslant 3$，$2^* = 2N/(N-2)$ 为Sobolev临界指数，参数 $\mu \in [0, (N-2)^2/8]$。该方程是经典热方程 $u_t - \Delta u = |u|^{2^*-2}u$ 的加权推广，其Hardy项 $-\mu u/|x|^2$ 在原点处产生强奇异性，显著改变解的局部行为与全局适定性结构。从物理角度看，该模型源于量子力学中电子在分子偶极场下的运动（此时 $1/|x|^2$ 势对应长程偶极-电荷相互作用），亦通过球面共形几何（经球极投影 $\mathbb{R}^N \cong \mathbb{S}^N \setminus \{\text{N-pole}\}$）与Yamabe问题紧密关联：寻找标量曲率为1的共形度量等价于求解(1.1)。数学上，该方程处于两个核心理论交界处：一是Hardy不等式与加权Sobolev空间理论，二是临界非线性抛物方程的全局动力学分类（整体存在 vs 有限时间爆破）。关键难点在于：当 $\mu > 0$ 时，标准Sobolev嵌入 $H^1(\mathbb{R}^N) \hookrightarrow L^{2^*}(\mathbb{R}^N)$ 失效，且 $L^2$ 能量泛函 $E(u) = \frac{1}{2}\int_{\mathbb{R}^N} \left(|\nabla u|^2 - \mu \frac{|u|^2}{|x|^2}\right) dx - \frac{1}{2^*} \int_{\mathbb{R}^N} |u|^{2^*} dx$ 不再正定——其二次型部分 $\int |\nabla u|^2 - \mu |u|^2/|x|^2$ 在 $\mu$ 接近临界值 $(N-2)^2/4$ 时可负，导致经典Nehari流形方法失效；同时，临界非线性 $|u|^{2^*-2}u$ 破坏尺度不变性与紧性，使能量控制不足以排除振荡或集中现象。已有工作主要集中于三类情形：(i) $\mu = 0$ 时的Fujita型爆破理论（临界指数 $q^* = 1 + 2/N$）；(ii) 次临界非线性 $|u|^{p-1}u$ ($p < 2^*$) 下带位势的抛物方程（如[9,10,24,25,33]），其Fujita指数依赖于位势衰减率；(iii) 非临界Hardy位势（$\mu < (N-2)^2/4$）下椭圆或抛物问题的正则性与存在性（如[3,12]）。然而，对 $\mu \in (0, (N-2)^2/8]$ 且非线性取Sobolev临界指数 $2^*$ 的抛物初值问题，既无适用于加权空间的精确Hardy-Sobolev不等式，也缺乏能处理临界增长与原点奇异性耦合的统一动力学框架。本文旨在填补这一空白：建立适配方程(1.1)的加权Hardy不等式，构造修正的势井方法，并在加权Gaussian空间中刻画解的三重动力学分岔（低、临界、高初始能量），从而首次给出Hardy-Sobolev抛物方程在临界非线性下的完整存在性-爆破判据。

2. 主要定理或主要结果

本文在加权Gaussian空间 $X = \overline{C_c^\infty(\mathbb{R}^N)}^{\|\nabla v\|_{K,2}}$ 中（其中权重 $K(y) = e^{|y|^2/4}$），通过自相似变换 $v(y,s) = (1+t)^{1/(p-1)} u(x,t)$, $t = e^s -1$, $x = (1+t)^{1/2} y$ 将(1.1)化为等价问题(2.1)，并证明以下核心结果：

定理3.1（低能量整体存在性）：设 $\mu \in [0, (N-2)^2/8]$, $v_0 \in X$ 满足 $E_K(v_0) < d$ 且 $v_0 \in W$（即 $D_K(v_0) > 0$），其中 $d = \inf_{v \in \mathcal{N}} E_K(v) > 0$ 为势井深度，$E_K(v) = \frac{1}{2} A(v) - \frac{1}{2^*} B(v)$, $A(v) = \int_{\mathbb{R}^N} \left(|\nabla v|^2 - \mu \frac{|v|^2}{|y|^2} - \beta |v|^2\right) K(y) dy$, $B(v) = \int_{\mathbb{R}^N} |v|^{2^*} K(y) dy$, $\beta = (N-2)/4$。则问题(2.1)存在唯一弱解 $v \in L^\infty(0,\infty; X) \cap L^2_{\text{loc}}(0,\infty; L^2_K(\mathbb{R}^N))$，且满足 $\|v(s)\|_{K,2} \leq C e^{-\lambda s}$，其中 $\lambda = \frac{1}{4} \min\{1, N/4\}$。

定理4.1（临界能量整体存在性）：若 $E_K(v_0) = d$ 且 $v_0 \in \overline{W}$，则问题(2.1)存在整体弱解 $v$，且 $\lim_{s \to \infty} \|v(s)\|_{K,2} = 0$；进一步，若 $v_0 \in \mathcal{N}$，则 $v(s) = e^{-\beta s} v_0$ 是稳态解。

定理5.1（高能量爆破）：若 $E_K(v_0) > d$ 且 $v_0 \in V$（即 $D_K(v_0) < 0$），则问题(2.1)的弱解在有限时间 $S^* < \infty$ 内爆破，即 $\limsup_{s \to S^*} \|v(s)\|_{K,2} = \infty$；更精确地，存在常数 $C > 0$ 使得 $\|v(s)\|_{K,2} \geq C (S^* - s)^{-1/(2^*-2)}$。

引理2.2（加权Hardy不等式）：对任意 $u \in X$，有

$$ \int_{\mathbb{R}^N} \frac{|u|^2}{|x|^2} K(x) dx \leq C_N \int_{\mathbb{R}^N} |\nabla u|^2 K(x) dx, $$

其中最优常数 $C_N = \left(\frac{2}{N-2}\right)^2$，且当 $\mu \leq (N-2)^2/8 < (N-2)^2/4$ 时，二次型 $A(v)$ 正定（引理2.5：$A(v) \geq \frac{1}{4} \|v\|_{K,2}^2$）。

引理2.7（势井深度性质）：定义修正Nehari流形 $\mathcal{N}_\delta = \{v \in X \setminus \{0\} : D_{K,\delta}(v) = 0\}$ 及 $d_\delta = \inf_{v \in \mathcal{N}_\delta} E_K(v)$，其中 $D_{K,\delta}(v) = \delta A(v) - B(v)$。则 $d_\delta$ 在 $\delta = 1$ 处取最大值 $d = d_1 > 0$，且 $d_\delta$ 关于 $\delta$ 单调：在 $(0,1]$ 上递增，在 $[1, (2^*)/2]$ 上递减。

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

本文的核心技术障碍在于：(1) 标准Sobolev嵌入在Hardy位势下失效，需构造适配方程(1.1)的加权嵌入与不等式；(2) 临界非线性破坏尺度不变性，使能量泛函缺乏紧性控制；(3) 原点奇异性导致经典势井方法中Nehari流形 $\mathcal{N} = \{v: A(v) = B(v)\}$ 不再是自然不变流形，且 $A(v)$ 可能非正。作者的真正创新点在于构建了一个三层分析框架：

第一层：加权Gaussian空间与最优Hardy不等式。区别于常规加权空间（如 $L^2(|x|^{-2}dx)$），作者引入Gaussian权重 $K(y) = e^{|y|^2/4}$，将自相似变换后的算子 $L v = -\Delta v - \frac{1}{2} y \cdot \nabla v$ 表达为散度形式 $L v = -K^{-1} \nabla \cdot (K \nabla v)$，从而赋予 $X$ 以Hilbert空间结构 $\|\nabla v\|_{K,2} = \left(\int |\nabla v|^2 K dy\right)^{1/2}$。在此空间中，作者给出了引理2.2的加权Hardy不等式，其最优常数 $C_N = (2/(N-2))^2$ 严格小于经典Hardy常数 $(N-2)^2/4$，这得益于权重 $K$ 对原点奇异性与无穷远衰减的双重平衡。该不等式是后续所有估计的基础，确保了 $A(v)$ 的强制性（引理2.5）。

第二层：修正势井方法与参数化Nehari流形。为克服 $A(v)$ 符号不定性，作者未采用传统 $\mathcal{N}$，而是定义一族参数化Nehari流形 $\mathcal{N}_\delta = \{v: \delta A(v) = B(v)\}$（见(2.3)），并系统研究 $d_\delta$ 的单调性（引理2.7）。关键洞察是：当 $\delta = 1$ 时，$\mathcal{N}_1 = \mathcal{N}$，而 $d_\delta$ 在 $\delta = 1$ 处取最大值 $d$，这表明 $\mathcal{N}$ 是“最深”的势井边界。作者进而定义修正势井 $W_\delta = \{v: D_{K,\delta}(v) > 0, E_K(v) < d_\delta\} \cup \{0\}$，并利用引理2.6建立 $\sqrt{A(v)}$ 与 $r(\delta) = \delta^{(N-2)/2} S_K^{N/2}$ 的阈值关系：$\sqrt{A(v)} < r(\delta)$ 当且仅当 $D_{K,\delta}(v) > 0$。此关系将能量条件 $E_K(v) < d$ 与梯度范数 $\sqrt{A(v)}$ 直接关联，为不变集论证提供定量依据。

第三层：Palais-Smale序列分析与渐近行为刻画。针对整体解的长时间行为（第6节），作者将解 $v(s)$ 视为参数 $s$ 的曲线，并证明其满足 $\lim_{s \to \infty} \|v_s(s)\|_{L^2_K} = 0$，即 $v(s)$ 是相应椭圆问题 $Lv = \beta v + \mu v/|y|^2 + |v|^{2^*-2}v$ 的Palais-Smale序列。结合引理2.1的紧嵌入 $X \hookrightarrow L^q_K$ ($q < 2^*$)，作者提取收敛子列并证明极限必为稳态解，从而确立 $v(s) \to 0$ 的强收敛性。此方法绕开了对解的显式先验估计，直接利用变分结构刻画渐近态。

4. 与之前工作的比较

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

1. 空间框架与不等式精度：此前Hardy抛物方程研究（如[9,10,24]）多基于标准 $L^p$ 或加权 $L^p(|x|^\alpha dx)$ 空间，其Hardy不等式常采用非最优常数（如 $C < (N-2)^2/4$），导致能量泛函 $A(v)$ 的下界估计粗糙。本文首创Gaussian加权空间 $X$，并证明引理2.2的最优加权Hardy不等式，其常数 $C_N = (2/(N-2))^2$ 与 $\mu$ 的允许范围 $[0, (N-2)^2/8]$ 精确匹配（因 $(2/(N-2))^2 \cdot (N-2)^2/8 = 1/2 < 1$），从而保证 $A(v)$ 的强制性（引理2.5）。相较之下，文献[3]在 $\mu < (N-2)^2/4$ 下使用标准Hardy不等式，但未处理临界非线性；文献[25]假设位势 $V(x)=1$，回避了原点奇异性。

2. 动力学分类的完整性：Fujita型工作（如[4,30]）仅针对次临界非线性 $|u|^{q-1}u$ ($q < 2^*$)，其临界指数 $q^*$ 与空间维数相关，但无法覆盖 $q = 2^*$。本文首次在临界情形下实现三重能量分类（低/临界/高），且对每类均给出定量结论：低能量解的指数衰减率 $e^{-\lambda s}$（$\lambda > 0$）比文献[7]中仅有的代数衰减更精确；临界能量解的收敛性证明（定理4.1）补充了[15,16]中未涉及的稳态解刻画；高能量爆破的速率下界 $(S^*-s)^{-1/(2^*-2)}$ 细化了[18,19]中仅有的存在性断言。

3. 方法论的结构性改进：经典势井方法（[23,27]）依赖固定Nehari流形 $\mathcal{N}$，而本文的修正方法通过参数化 $\mathcal{N}_\delta$ 揭示了势井深度 $d_\delta$ 的内在结构（引理2.7），证明 $\delta=1$ 是临界参数。这一框架比[15,16]中仅引入单一家族 $\mathcal{N}_\delta$ 更深入，因其定量刻画了 $d_\delta$ 的单调性与极值点，为能量阈值 $d$ 的唯一性提供了变分依据。此外，将自相似变换与Palais-Smale分析结合（第6节），是首次将椭圆变分技巧系统应用于临界抛物方程的长时间行为研究。

4. 适用范围与代价：本文优势在于处理了最难的临界非线性与中等强度Hardy位势（$\mu \leq (N-2)^2/8$），但代价是结果局限于自相似坐标系下的加权空间 $X$，未给出原变量 $u(x,t)$ 的显式爆破模式（如爆破点位置、轮廓函数）。相较文献[18,19]在 $\mathbb{R}^1$ 中构造自相似爆破解，本文未显式构造爆破解，而是通过不变集论证确立爆破必然性。从提供内容中无法确认本文是否推广至变系数位势或非径向初值，亦未讨论 $\mu > (N-2)^2/8$ 时的情形（此时 $A(v)$ 可能非正，需全新框架）。

8. 与RN中Hardy-Sobolev热方程自相似解相关的变分问题

Variational problems related to self-similar solutions of Hardy-Sobolev heat equation in RN

作者: Fei Fang, Zhong Tan

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.24804v1

摘要翻译: 本文应用自相似变换，将带Sobolev-Hardy项的抛物型方程\begin{align*} u_t-\Delta u= \frac{|u|^{q-2}u}{\left|x\right|^s} & \text { in } \mathbb{R}^N \times(0, \infty), \end{align*}转化为如下椭圆型方程\begin{equation*} -\Delta v-\frac{1}{2} y \cdot \nabla v=\alpha v+ \frac{|v|^{q-2} v}{|y|^s}, \end{equation*}其中$2 < q \leq 2^*(s)=\frac{2 N-2 s}{N-2}, 0 \leq s < 2, \alpha=\frac{2-s}{2q-4}$。对此方程，我们建立了加权Hardy不等式与Sobolev不等式。进一步，借助变分方法，我们在次临界情形下得到了无穷多个解，并在临界情形下证明了解的存在性。此外，我们还利用Pohozaev恒等式，在某些条件下建立了解的不存在性。

Comment: 无

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

1. 研究对象和背景

本文研究的原始模型是定义在 $\mathbb{R}^N \times (0, \infty)$ 上的带Hardy-Sobolev型奇异势的半线性抛物方程：

$$ u_t - \Delta u = \frac{|u|^{q-2}u}{|x|^s}, \quad (x,t) \in \mathbb{R}^N \times (0,\infty), \tag{1.1} $$

其中参数满足 $2 < q \leq 2^*(s) := \frac{2N-2s}{N-2}$，$0 \leq s < 2$，$N \geq 3$。该方程是经典热方程的奇异扰动版本，其右端项 $|x|^{-s}$ 在原点处具有非可积奇性（当 $s \geq 2$ 时 $|x|^{-s} \notin L^1_{\mathrm{loc}}(\mathbb{R}^N)$），故不能通过标准的 $L^p$-理论或光滑势方法处理；而 $q$ 的上界 $2^*(s)$ 恰为加权Sobolev嵌入 $D^{1,2}(\mathbb{R}^N; |x|^{-s}dx) \hookrightarrow L^{2^*(s)}(\mathbb{R}^N; |x|^{-s}dx)$ 的临界指数，标志着变分泛函从强制性到缺乏紧性的相变阈值。该模型物理上刻画了在存在中心对称奇点（如引力源、点缺陷、几何奇点）约束下的异常热传导与浓度演化过程，例如多孔介质中带中心吸附势的扩散、量子点附近的载流子局域化、以及具奇点的宇宙学标量场动力学。

为提取长时间渐近结构，作者采用自相似变换：

$$ u(x,t) = t^{-\alpha} v(y), \quad y = \frac{x}{\sqrt{t}}, \quad \alpha = \frac{2-s}{2q-4} > 0. \tag{1.2'} $$

代入 (1.1) 后，$v = v(y)$ 满足如下椭圆型Fokker–Planck型方程：

$$ -\Delta v - \frac{1}{2} y \cdot \nabla v = \alpha v + \frac{|v|^{q-2}v}{|y|^s}, \quad y \in \mathbb{R}^N. \tag{1.2} $$

该方程等价于加权算子 $L v := -K^{-1} \nabla \cdot (K \nabla v)$ 作用下的特征值问题，其中权重函数 $K(y) = e^{|y|^2/4}$ 是 Ornstein–Uhlenbeck 算子的不变测度密度。因此，(1.2) 并非标准的欧氏拉普拉斯方程，而是定义在加权空间 $L^2_K(\mathbb{R}^N) := L^2(\mathbb{R}^N; K(y)dy)$ 上的非均匀椭圆问题，其主部具有强退化性（当 $|y| \to \infty$ 时 $K(y) \to \infty$ 指数增长），且非线性项受双重奇性耦合：一是 $|y|^{-s}$ 在原点的Hardy奇性，二是 $|v|^{q-2}v$ 在 $q = 2^*(s)$ 处的Sobolev临界增长。研究的核心困难在于：(i) 加权空间 $X := \overline{C_c^\infty(\mathbb{R}^N)}^{\|\nabla \cdot\|_{K,2}}$ 中缺乏标准 Sobolev 嵌入的紧性（因 $K(y)$ 非常数）；(ii) $|y|^{-s}$ 导致 $0$ 点附近能量集中，破坏 Poincaré 不等式与 Rellich 型紧嵌入；(iii) 当 $q = 2^*(s)$ 时，泛函 $E_K(v)$ 在 $X$ 上既不满足 Palais–Smale 条件，也不具备山路引理所需的几何结构。已有工作如 Escobedo–Kavian [5]（$s=0$）、Tan [20]（$s=0$，低能初值）及 [12]（$s=0$，衰减率）均未处理 $s>0$ 的Hardy–Sobolev耦合情形。本文试图填补的关键空白是：建立适用于 $s \in [0,2)$ 的加权Hardy–Sobolev不等式体系，并在此基础上，首次系统解决加权Fokker–Planck–Hardy–Sobolev方程 (1.2) 在亚临界 ($q < 2^*(s)$) 和临界 ($q = 2^*(s)$) 两种情形下的变分适定性——即无穷多解的存在性、临界能量水平下的解存在性与Pohozaev型非存在性。

2. 主要定理或主要结果

本文在加权Hilbert空间 $X = \{ v \in L^2_K(\mathbb{R}^N) : \|\nabla v\|_{K,2} < \infty \}$ 上建立了以下核心结论：

  1. 【加权不等式】证明了最优加权Hardy不等式：对任意 $u \in X$，

    $$ \int_{\mathbb{R}^N} \frac{|u|^2}{|x|^2} K(x)\,dx \leq C_N \int_{\mathbb{R}^N} |\nabla u|^2 K(x)\,dx, \quad C_N = \left( \frac{2}{N-2} \right)^2, \tag{2.1} $$

    且 $C_N$ 为最佳常数；并建立了加权Sobolev–Hardy不等式：对任意 $2 \leq q \leq 2^*(s)$，存在 $C > 0$ 使得

    $$ \left( \int_{\mathbb{R}^N} \frac{|u|^q}{|x|^s} K(x)\,dx \right)^{2/q} \leq C \int_{\mathbb{R}^N} |\nabla u|^2 K(x)\,dx. \tag{2.4} $$

  2. 【亚临界情形】（Theorem 3.1）：若 $2 < q < 2^*(s)$，则问题 (1.2) 存在无穷多个变分解 $\{\pm v_k\}_{k=1}^\infty \subset X$，满足 $E_K(\pm v_k) \to +\infty$，其中能量泛函为

    $$ E_K(v) = \frac{1}{2} \int_{\mathbb{R}^N} |\nabla v|^2 K(y)\,dy - \frac{\alpha}{2} \int_{\mathbb{R}^N} |v|^2 K(y)\,dy - \frac{1}{q} \int_{\mathbb{R}^N} \frac{|v|^q}{|y|^s} K(y)\,dy. $$

    该结论由Fountain定理保证，其关键尺度关系为：令 $Z_k = \bigoplus_{j=k}^\infty \mathbb{R} e_j$（$\{e_j\}$ 为 $X$ 的标准正交基），定义

    $$ \beta_k = \sup_{\substack{v \in Z_k \\ \|\nabla v\|_{K,2} = 1}} \left( \int_{\mathbb{R}^N} \frac{|v|^q}{|y|^s} K(y)\,dy \right)^{1/q}, \quad \gamma_k = \sup_{\substack{v \in Z_k \\ \|\nabla v\|_{K,2} = 1}} \left( \int_{\mathbb{R}^N} |v|^2 K(y)\,dy \right)^{1/2}, $$

    则 $\beta_k \to 0$，$\gamma_k \to 0$，且取 $r_k = (2c_2 q \beta_k^q)^{1/(2-q)}$，有 $\inf_{\substack{v \in Z_k \\ \|\nabla v\|_{K,2} = r_k}} E_K(v) \geq \frac{q-2}{4q} \left( \frac{1}{2q c_2 \beta_k^q} \right)^{2/(q-2)} \to +\infty$。

  3. 【临界情形】（Theorems 4.1 & 4.2）：定义加权最佳常数

    $$ S_{K,\alpha} = \inf_{v \in X \setminus \{0\}} \frac{\int_{\mathbb{R}^N} K|\nabla v|^2 - \alpha K v^2}{\left( \int_{\mathbb{R}^N} K |v|^{2^*(s)} |y|^{-s} \right)^{2/2^*(s)}}, $$

    则：(i) 若 $\alpha \leq N/4$（$N \geq 3$），则 (1.2) 仅有零解（Theorem 4.1）；(ii) 若 $N \geq 4$ 且 $\alpha \in (N/4, N/2)$，则 (1.2) 至少存在一个非平凡解；若 $\alpha \leq N/4$ 或 $\alpha \geq N/2$，则无正解；(iii) 若 $N = 3$，则当 $\alpha \in (1, 3/2)$ 时存在非平凡解，而当 $\alpha \leq 3/4$ 或 $\alpha \geq 3/2$ 时无正解。

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

本文的核心技术障碍在于：加权空间 $X$ 中，由于权重 $K(y) = e^{|y|^2/4}$ 的指数增长性，标准的 Rellich–Kondrachov 紧嵌入完全失效；同时，Hardy项 $|y|^{-s}$ 在原点引入第二重奇性，导致 $0$ 点邻域内能量无法控制，使得传统的截断–估计策略（如 Brezis–Lieb 引理直接应用）不可行。作者的真正创新点在于构建了一套适配双重奇性的加权不等式链，并设计了与之匹配的变分框架。

第一，加权Hardy–Sobolev不等式的构造是全文基石。Lemma 2.3 并非简单引用文献[5]，而是给出了 $C_N = (2/(N-2))^2$ 的严格最优性证明：通过构造径向试探函数 $V(r) = r^{-\gamma} e^{-r^2/4}$，其中 $\gamma = (N-2+2\varepsilon)/2$，精确计算其梯度与 $|x|^{-2}$ 加权范数的比值，证明当 $\varepsilon \to 0^+$ 时，比值趋近于 $C_N$，从而确立最优性。这一构造克服了经典Hardy不等式在加权测度下常数漂移的难题。进而，Lemma 2.4 将此与加权Sobolev嵌入耦合：利用Hölder不等式将 $\int |u|^{2^*(s)} |x|^{-s} K$ 拆分为 $\int |u|^2 |x|^{-2} K$ 与 $\int |u|^{2^*} K$ 的乘积形式，再分别用 Lemma 2.3 和 Lemma 2.1 控制，最终得到 $\left( \int |u|^{2^*(s)} |x|^{-s} K \right)^{2/2^*(s)} \leq C \int |\nabla u|^2 K$。该不等式是后续所有紧性论证的前提。

第二，在亚临界情形（Section 3），作者未使用常规的山路引理，而是采用 Fountain 定理，其关键在于证明 $\beta_k \to 0$ 和 $\gamma_k \to 0$。此处的精妙之处在于：利用 Lemma 2.4 的紧嵌入 $X \hookrightarrow L^q_K(\mathbb{R}^N; |y|^{-s})$（$q < 2^*(s)$），对任意 $v_k \in Z_k$ 满足 $\|\nabla v_k\|_{K,2} = 1$，由 $v_k \rightharpoonup 0$ 在 $X$ 中弱收敛，可推出 $v_k \to 0$ 在 $L^q_K(\mathbb{R}^N; |y|^{-s})$ 中强收敛，从而 $\beta_k \to 0$。这绕开了对 $L^\infty$ 或 Hölder 连续性的需求，是处理加权空间中‘高频振荡’解的典型技巧。

第三，在临界情形（Section 4），作者引入加权最佳常数 $S_{K,\alpha}$ 并建立其与欧氏最佳常数 $S_\alpha$ 的比较。非存在性（Theorem 4.1）的证明依赖于Pohozaev恒等式在加权空间中的修正形式：对任意解 $v \in X$，有

$$ \frac{N-2}{2} \int_{\mathbb{R}^N} |\nabla v|^2 K - \frac{N}{2} \alpha \int_{\mathbb{R}^N} v^2 K = \frac{N-s}{2^*(s)} \int_{\mathbb{R}^N} \frac{|v|^{2^*(s)}}{|y|^s} K. $$

结合 $S_{K,\alpha}$ 的定义与Poincaré不等式 $\lambda_1 = N/2$，导出矛盾的条件 $\alpha \leq N/4$。存在性则通过极小化序列的集中紧性分析完成：利用 $S_{K,\alpha} < S_{K,0}$（当 $\alpha > 0$）及 $S_{K,\alpha}$ 的可达性，证明极小化序列在平移后不逃逸至无穷远，从而强收敛到非零解。

4. 与之前工作的比较

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

  1. 【假设放宽】Escobedo–Kavian [5] 及后续工作 [2,6,8] 均仅处理 $s = 0$ 的经典自相似热方程，其对应的椭圆方程为 $-\Delta v - \frac{1}{2} y \cdot \nabla v = \alpha v + |v|^{q-2}v$，此时无需Hardy不等式，且 $L^2_K$ 嵌入 $X \hookrightarrow L^q_K$ 对所有 $q \in [2,2^*)$ 紧致（$2^* = 2N/(N-2)$）。本文将 $s$ 放宽至 $[0,2)$，首次建立了 $s > 0$ 下的加权Hardy–Sobolev不等式 (2.1) 与 (2.4)，这是处理 $|x|^{-s}$ 奇性的唯一工具。特别地，当 $s \to 2^-$ 时，$2^*(s) \to 2$，临界情形退化为 $L^2$ 边界，此时本文的 $C_N$ 最优性证明（Lemma 2.3）比文献[15]中仅针对 $s=0$ 的论证更具普适性。

  2. 【时间尺度与适用范围】Tan [20] 与 [12] 关注原始抛物方程 (1.1) 的初值问题，其结论依赖于初始数据的能量大小（低能全局存在，高能可能爆破），且仅给出 $s=0$ 下的衰减率 $\|u(t)\|_{L^2} \sim t^{-\sigma}$。本文则通过自相似变换，将问题完全转化为稳态椭圆问题 (1.2)，其解 $v(y)$ 直接对应原始方程的自相似解 $u(x,t) = t^{-\alpha} v(x/\sqrt{t})$，因此结论刻画的是 $t \to \infty$ 时所有自相似轨道的结构，而非单个初值的长时间行为。Theorem 3.1 给出的无穷多解意味着原始热方程存在无穷多个不同渐近剖面的自相似演化模式，这是初值问题框架无法揭示的全局动力学图景。

  3. 【代价与局限】本文的代价是：所有结果均建立在加权空间 $X$ 上，其范数 $\|\nabla v\|_{K,2}$ 包含指数权重，导致解 $v(y)$ 在 $|y| \to \infty$ 时必须以 $e^{-|y|^2/8}$ 速率衰减，这比欧氏空间中 $|y|^{-(N-2)/2}$ 的衰减更强。因此，本文解对应的原始解 $u(x,t) = t^{-\alpha} v(x/\sqrt{t})$ 在固定 $t$ 下具有高斯型衰减，而无法覆盖具有代数衰减的更广义解。此外，Theorem 4.2 中对 $N=3$ 与 $N\geq4$ 的区分，以及 $\alpha$ 范围的精细划分（如 $N=3$ 时 $\alpha \in (1,3/2)$），表明临界情形的解存在性高度依赖于维数与参数的共振效应，这比 $s=0$ 时的单一阈值 $\alpha < N/2$ 更为复杂。作者未显式展开与 [3,7,15] 中动力学分类的对接，亦未讨论解的稳定性或奇点类型（如 type I/II blow-up），这些是未来工作的重要方向。

9. 含真空的可压 Navier-Stokes/Allen-Cahn 系统强解的适定性与爆破准则

Well-posedness and blow-up criterion for strong solutions of the compressible Navier-Stokes/Allen-Cahn system with vacuum

作者: Yinghua Li, Haoran Zheng, Junquan Zhou

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24568v1

摘要翻译: 本文研究定义在有界区域 $Ω\subset\mathbb R^3$ 上的可压 Navier-Stokes/Allen-Cahn 系统的强解，允许初值出现真空。该系统的一个特征是密度与 Allen-Cahn 方程之间的强耦合，这导致在真空区域出现强退化性。在初值相场变量满足相容性条件的前提下，我们建立了局部强解的存在唯一性：要求 $0\leρ_0\in W^{1,q}$（其中 $q\in(3,6)$）、$u_0\in H_0^1$ 且 $χ_0\in H^2$。借助时间加权估计，速度初值无需相容性条件，但此类估计在证明唯一性时引入了奇异性。随后，我们建立了一个刻画该局部强解可能在有限时刻发生破裂的准则，即当量 $\|\nabla u\|_{L_t^{1} L_x^{\infty}}$、$\|u\|_{L_t^{2} L_x^{\infty}}$ 与 $\|\nabla χ\|_{L_t^{2} L_x^{\infty}}$ 发生爆破时，解将不再延拓。

Comment: 无

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

1. 研究对象和背景

本文研究三维有界区域 $\Omega \subset \mathbb{R}^3$ 上带初始真空的可压Navier-Stokes/Allen-Cahn (NSAC) 耦合系统的强解理论。该系统由Blesgen [3] 提出，其数学形式为：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0, \\ \rho \partial_t u + \rho (u \cdot \nabla) u + \nabla P(\rho) = \mathcal{L} u - \mathrm{div}\left( \nabla \chi \otimes \nabla \chi - \frac{|\nabla \chi|^2}{2} I \right), \\ \rho \partial_t \chi + \rho (u \cdot \nabla) \chi = -\mu, \\ \rho \mu = -\Delta \chi + \rho f(\chi), \end{cases} \tag{1.1} $$

其中 $\rho = \rho(x,t) \geq 0$ 为总密度，$u = u(x,t)$ 为平均速度，$\chi = \chi(x,t) \in \mathbb{R}$ 为相场变量（表征两相浓度差），$\mu = \mu(x,t)$ 为化学势；压力 $P(\rho) = A \rho^\gamma$，$A > 0$, $\gamma > 1$；Lamé算子 $\mathcal{L} u = \nu \Delta u + (\nu + \lambda) \nabla (\mathrm{div}\, u)$，粘性系数满足 $\nu > 0$, $\lambda = \eta - \frac{2}{3}\nu$, $2\nu + 3\lambda \geq 0$；自由能密度 $f(s) = F'(s)$，Landau型势能 $F(s) = \frac{1}{4}(s^2 - 1)^2$，故 $f(s) = s(s^2 - 1)$。边界条件为 $u|_{\partial\Omega} = 0$, $\partial_n \chi|_{\partial\Omega} = 0$；初始条件为 $(\rho, \rho u, \chi)|_{t=0} = (\rho_0, \rho_0 u_0, \chi_0)$，且允许 $\rho_0$ 在 $\Omega$ 的非空开子集上恒为零（即初始真空）。

该系统的核心困难在于强耦合导致的退化性：Allen-Cahn方程中出现 $\rho^2 \chi_t$ 项，当 $\rho \to 0$ 时，时间导数项完全消失，使方程在真空区退化为椭圆型 $-\Delta \chi + \rho f(\chi) = 0$，而右端 $\rho f(\chi)$ 又随 $\rho$ 消失，物理上意味着真空区无相分离驱动力，故要求 $\Delta \chi_0 = 0$ 当 $\rho_0 = 0$。这一结构远超经典可压Navier-Stokes方程（仅含 $\rho u_t$）的退化程度，亦不同于Cahn-Hilliard型（含高阶 $\Delta \mu$）的正则化机制。已有工作多限于一维情形：Ding-Li-Luo [13]、Chen-Guo [10] 等在无真空或强兼容性条件下建立强解；Chen-Zhu [12] 引入密度-相场依赖粘性并施加兼容性条件；Su [41] 利用时间加权估计去除速度兼容性但未处理三维真空。三维情形下，Feireisl等 [21]、Chen-Wen-Zhu [11] 仅得弱解；Kotschote [27] 仅对非等熵情形给出局部适定性，且未允许真空。因此，三维带真空NSAC强解的存在性是长期悬而未决的公开问题。本文旨在填补这一空白：在仅对相场变量 $\chi_0$ 施加轻度兼容性条件（见(1.4)）、不对速度 $u_0$ 作任何兼容性假设的前提下，建立局部强解的存在唯一性，并给出首个三维爆破准则。

2. 主要定理或主要结果

定理1.1（局部存在唯一性）：设 $\Omega \subset \mathbb{R}^3$ 为光滑有界区域，初始数据满足 $0 \leq \rho_0 \in W^{1,q}(\Omega)$，$q \in (3,6)$，$u_0 \in H_0^1(\Omega)$，$\chi_0 \in H^2(\Omega)$，且满足相场兼容性条件

$$ \begin{cases} \Delta \chi_0 = \rho_0 h, & \text{a.e. in } \Omega, \\ \partial_n \chi_0 = 0, & \text{on } \partial\Omega, \end{cases} \quad \text{其中 } h \in L^2(\Omega). \tag{1.4} $$

则存在依赖于 $A, \gamma, \nu, \lambda, q, \Omega$ 和 $\Phi_0 = \|\rho_0\|_{W^{1,q}} + \|\nabla u_0\|_{L^2} + \|\chi_0\|_{H^2} + \|h\|_{L^2}$ 的 $T_0 > 0$，使得系统(1.1)–(1.3)存在唯一强解 $(\rho, u, \chi, \mu)$，满足定义1.1中的正则性：

  • $0 \leq \rho \in L^\infty(0,T_0; W^{1,q}) \cap C([0,T_0]; C(\Omega))$，$\rho_t \in L^\infty(0,T_0; L^2)$；
  • $u \in L^\infty(0,T_0; H_0^1) \cap L^2(0,T_0; H^2) \cap L^1(0,T_0; W^{2,q})$，$\sqrt{\rho} u \in C([0,T_0]; L^2)$，$\sqrt{t} \nabla^2 u \in L^\infty(0,T_0; L^2) \cap L^2(0,T_0; L^q)$，$\sqrt{\rho} u_t \in L^2(0,T_0; L^2)$，$\sqrt{t} u_t \in L^2(0,T_0; H_0^1)$；
  • $\chi \in C([0,T_0]; H^1) \cap L^\infty(0,T_0; H^2) \cap L^2(0,T_0; H^3)$，$\sqrt{t} \nabla^3 \chi \in L^\infty(0,T_0; L^2)$，$\chi_t \in L^2(0,T_0; H^1)$，$\rho \chi_t \in L^\infty(0,T_0; L^2)$，$\sqrt{t} \nabla \chi_t \in L^\infty(0,T_0; L^2)$；
  • $\mu \in L^\infty(0,T_0; L^2)$。

定理1.2（爆破准则）：设 $(\rho, u, \chi, \mu)$ 为上述强解，$T^* \in (0,\infty)$ 为其最大存在时间，则必有

$$ \lim_{T \to T^*} \left( \|\nabla u\|_{L^1(0,T; L^\infty)} + \|u\|_{L^2(0,T; L^\infty)} + \|\nabla \chi\|_{L^2(0,T; L^\infty)} \right) = \infty. \tag{1.5} $$

该准则表明，若解在有限时间 $T^*$ 处爆破，则上述三项范数至少有一项趋向无穷；其形式与Li-Ding-Huang [33] 对三维NSAC的准则一致，但本文首次在真空情形下证明。

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

证明的核心障碍源于 $\rho^2 \chi_t$ 导致的强退化：在真空区，Allen-Cahn方程失去抛物性，标准能量估计无法控制 $\chi_t$ 的初始奇异性。作者的真正创新点在于构建了一套分层的时间加权估计框架，并严格区分速度与相场的处理策略。

首先，针对速度场 $u$，作者采用 $\sqrt{t}$ 加权的 $H^2$-估计以规避兼容性条件。关键引理2.3通过Lamé方程的 $W^{2,q}$-正则性给出：

$$ \int_0^T \|\nabla^2 u\|_{L^q} \, dt \leq C T^{\frac{6-q}{4q}} \Phi^{\gamma+1}(T), $$

其中 $\Phi(t)$ 是包含 $\|\rho\|_{W^{1,q}}$、$\|\nabla u\|_{L^2}^2$、$\|\sqrt{t} \sqrt{\rho} u_t\|_{L^2}^2$、$\|\chi\|_{H^2}^2$ 等项的复合能量泛函。此估计结合Gronwall不等式（引理2.4）导出密度的一致有界性 $\sup_{t\in[0,T]} \|\rho\|_{W^{1,q}} \leq C$，前提是小时间尺度 $T^{\frac{6-q}{4q}} \Phi^{\gamma+1}(T) \leq \varepsilon_0$。该步骤成功绕开了对 $u_0$ 的兼容性要求，但代价是引入了 $\sqrt{t}$ 奇异性，使唯一性证明需额外处理。

其次，针对相场 $\chi$，作者必须直面退化性。关键突破在于对(1.1)₃–(1.1)₄进行重写与高阶微分：将原方程组合并为 $\rho^2 \chi_t + \rho^2 u \cdot \nabla \chi = \Delta \chi - \rho f(\chi)$（即(2.6)），再对时间求导得(2.7)，乘以 $\chi_t$ 并积分。此时，$\rho^2 \chi_{tt}$ 项产生 $\frac{1}{2} \frac{d}{dt} \|\rho \chi_t\|_{L^2}^2$，而退化项 $2\rho \rho_t \chi_t$ 等通过 $\|\rho\|_{L^\infty}$ 有界性及 $\|\nabla u\|_{L^2} \|\nabla^2 u\|_{L^2}$ 控制（引理2.5）。特别地，初始值 $\lim_{\tau \to 0} \|\rho \chi_t\|_{L^2}^2(\tau)$ 的估计直接依赖兼容性条件(1.4)：

$$ \lim_{\tau \to 0} \|\rho \chi_t\|_{L^2}^2(\tau) \leq C \left( \|\rho_0 (u_0 \cdot \nabla \chi_0)\|_{L^2}^2 + \left\| \frac{\Delta \chi_0}{\rho_0} \right\|_{L^2}^2 + \|f(\chi_0)\|_{L^2}^2 \right) \leq C \left( \|\rho_0\|_{L^\infty}^2 \|u_0\|_{L^6}^2 \|\nabla \chi_0\|_{L^3}^2 + \|h\|_{L^2}^2 + 1 \right), $$

其中 $\| \Delta \chi_0 / \rho_0 \|_{L^2}$ 的有界性正是由 $\Delta \chi_0 = \rho_0 h$ 保证。这解释了为何兼容性条件对 $\chi_0$ 不可省略。

最后，为获得 $\chi$ 的高阶正则性，作者对(2.6)作用 $\nabla$ 得 $\|\nabla \Delta \chi\|_{L^2}$ 的估计（引理2.6），关键项 $\|\rho^2 \nabla \chi_t\|_{L^2} \leq C \|\rho\|_{L^\infty}^2 \|\nabla \chi_t\|_{L^2}$ 与 $\|\nabla \Delta \chi\|_{L^2}^{1/2}$ 构成吸收不等式，最终导出 $\int_0^T \|\nabla^3 \chi\|_{L^2}^2 \, dt \leq C$。整个框架的枢纽是时间加权范数与退化结构的精确匹配：$\sqrt{t}$ 加权用于 $u$ 以消兼容性，而 $\rho$-加权（如 $\rho \chi_t$）用于 $\chi$ 以驯服退化。

4. 与之前工作的比较

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

1. 维度一：真空处理与兼容性条件。此前所有三维NSAC结果（如Feireisl [21], Chen-Wen-Zhu [11]）均要求 $\rho_0 \geq \delta > 0$，而一维工作（Chen-Guo [10], Chen-Zhu [12]）虽允许真空，但强制施加 $u_0$ 的兼容性条件（如 $[\nu(\rho_0,\chi_0) u_{0x}]_x - (A \rho_0^\gamma)_x - \frac{1}{2}(\chi_{0x}^2)_x = \rho_0 f(x)$）。本文首次在三维实现“半兼容性”：仅对 $\chi_0$ 要求(1.4)，对 $u_0$ 完全免去兼容性。这一推进的代价是引入 $\sqrt{t}$ 时间加权，导致 $\sqrt{t} \nabla^2 u$ 等范数出现，而经典框架（如Huang-Li-Xin [23]）中对应范数为无权的 $\nabla^2 u$。

2. 维度二：相场边界条件。所有已知一维真空结果均采用Dirichlet条件 $\chi|_{x=0,1} = 0$，而本文采用Neumann条件 $\partial_n \chi|_{\partial\Omega} = 0$，更符合物理实际（无通量边界）。该选择使 $\chi$ 的能量估计（引理2.2）中 $\|\nabla \chi\|_{L^2}$ 项自然出现，但增加了 $\nabla^3 \chi$ 估计的难度（引理2.6需处理边界项）。

3. 维度三：爆破准则的适用性。Li-Ding-Huang [33] 的准则针对无真空三维NSAC，其证明依赖 $\rho \geq \delta > 0$ 以保证 $\mu$ 的正则性；Zhang [43] 将 $\|\nabla \chi\|_{L^2(0,T;L^\infty)}$ 放宽至 $\|\nabla \chi\|_{L^2(0,T;\mathrm{BMO})}$，但仍未处理真空。本文准则(1.5)是首个在真空情形下成立的三维结果，其证明关键在于利用 $\rho \chi_t \in L^\infty(0,T;L^2)$ 和 $\nabla \chi_t \in L^2(0,T;L^2)$ 推出 $\nabla \chi \in L^2(0,T;L^\infty)$ 的Sobolev嵌入（因 $\Omega$ 有界，$H^2 \hookrightarrow L^\infty$），而该嵌入的常数不依赖 $\rho$ 的下界。

4. 维度四：方法论差异。与热传导可压NS方程（需Lagrangian坐标处理 $\rho_0 \theta_0$ 初始数据 [31]）不同，本文因直接给定 $\chi_0 \in H^2$ 且无热效应，得以全程在Eulerian坐标下完成存在性与唯一性证明（Remark 1.3）。这是结构优势的体现，但亦意味着无法借用Lagrangian方法中成熟的粒子轨迹技术来处理真空奇点。综上，本文在真空NSAC领域实现了存在性、兼容性、边界条件、爆破准则四重突破，但尚未解决全局存在性（仍为开放问题），且时间加权范数的存在限制了对长时间行为的直接分析。

10. 具有非局部动力学的标量守恒律的障碍问题

The obstacle problem for scalar conservation laws with nonlocal dynamics

作者: Paulo Amorim, Alexander Keimer, Lukas Pflug, Jakob Rodestock

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24677v1

摘要翻译: 本文中，我们提出一种方法来求解一维非局部守恒律方程，使其满足一个依赖于空间的映射——称为障碍（obstacle）。该方法通过推广已有的局部守恒律结果实现：我们考虑速度的一个松弛化（relaxation），该松弛化显式地依赖于障碍。我们证明了该松弛化问题解的存在性，并证明：当松弛映射收敛于 Heaviside 型函数时，相应的解收敛于一个间断非局部守恒律方程的弱解。此外，在若干情形下，我们可对极限通量（limiting flux）加以刻画。\ 本文最后给出一项数值研究，用以说明上述收敛性。

Comment: 17 pages, 5 figures

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

1. 研究对象和背景

本文研究一维非局部标量守恒律在空间依赖障碍约束下的解的存在性与极限行为。核心模型是如下非局部守恒律：

$$ \partial_t q(t, x) + \partial_x \big( q(t, x)\, U(W[q](t, x)) \big) = 0, \quad (t,x) \in (0,T) \times \mathbb{R}, $$

其中非局部加权平均算子定义为

$$ W[q](t, x) := \int_{x}^{\infty} \gamma(y - x)\, q(t, y)\, dy, $$

核函数 $\gamma \in L^1(\mathbb{R}_{\ge 0}; \mathbb{R}_{\ge 0})$ 满足 $\int_0^\infty \gamma(s)\, ds = 1$，且通常取为单调递减（如指数衰减），从而体现下游交通密度对当前速度的抑制效应。速度函数 $U \in C^2(\mathbb{R})$ 被假定为单调减且正（即 $U' \le 0$，且在障碍值域内 $U \ge U_{\min} > 0$），这是交通流建模中典型的“密度越高、速度越低”的Folklore本构关系。

障碍问题要求解 $q$ 在任意时刻满足点态约束 $q(t,x) \le o(x)$，其中 $o \in W^{2,\infty}(\mathbb{R}; \mathbb{R}_{>0})$ 是给定的空间依赖障碍函数（例如路障、限速区或瓶颈位置）。该约束不能通过标准初边值条件实现，而需嵌入方程结构本身——这正是本文的核心挑战。与椭圆/抛物型障碍问题（如Signorini问题、膜问题）不同，双曲型障碍问题缺乏自然的能量泛函和最大值原理；其解通常不连续，且障碍激活区域（即 $q = o$ 的集合，称为 coincidence set）随时间演化，导致通量 $q U(W[q])$ 出现不连续依赖性，破坏经典解的存在性。

已有工作主要分三类：(i) 局部守恒律障碍问题：Lévi [47] 首次建立一维框架，后续工作采用源项法 [5,10,11]、投影法 [58] 或逼近法 [25]；Amorim 等 [4] 提出关键的光滑罚函数方法：引入 $V_\varepsilon(o - q)$ 作为速度乘子，其中 $V_\varepsilon(x) = 1 - e^{-x/\varepsilon}$ 是 $\chi_{\mathbb{R}_{>0}}$ 的单调光滑逼近，使 $q$ 接近 $o$ 时速度衰减，从而“自动减速”以避免超障。(ii) 非局部守恒律理论：已知方程 $\partial_t q + \partial_x(q U(W[q])) = 0$ 具有优异的正则保持性——若 $q_0 \in BV(\mathbb{R})$，则解 $q \in C([0,T]; BV(\mathbb{R}))$ 且唯一，无需熵条件 [40,44,3]；这与局部方程 $\partial_t q + \partial_x(q U(q)) = 0$ 形成鲜明对比（后者即使光滑初值也产生激波，需Kruzhkov熵解）[46,31]。(iii) 非局部-局部奇异性极限：当核 $\gamma$ 逼近 Dirac 测度时，非局部解收敛于局部熵解 [19,3,26]；但本文明确声明不研究此极限下的障碍问题耦合。

本文试图填补的关键空白在于：将 [4] 的罚函数策略系统性地推广至非局部-非线性耦合情形。这引发本质新困难：(1) 通量 $V_\varepsilon(o - q)\, q\, U(W[q])$ 同时含 $q$ 的点态非线性（罚项）、乘性非线性（$q U$）及非局部依赖（$W[q]$），导致标准能量估计失效；(2) $V_\varepsilon$ 的导数 $V_\varepsilon'$ 在 $\varepsilon \to 0$ 时无界（$\|V_\varepsilon'\|_{L^\infty} = \varepsilon^{-1} e^{-1} \to \infty$），使经典紧性论证（如Arzelà–Ascoli）无法直接应用；(3) 非局部算子 $W$ 与罚函数 $V_\varepsilon(o - q)$ 的复合破坏了比较原理的直接可得性，需精细控制 $\partial_x W[q]$ 等高阶项。因此，本文并非简单复刻局部结果，而是构建了一套适配非局部结构的、带参数一致性的分析框架，旨在严格证明：当 $\varepsilon \to 0$ 时，罚函数解 $q_\varepsilon$ 收敛到一个满足障碍约束、质量守恒、有限传播速度的弱解，并精确刻画其极限通量结构。

2. 主要定理或主要结果

本文在以下假设下建立了完整存在性与收敛性理论：

(Assumption 1) 设 $U \in C^2(\mathbb{R})$, $U' \le 0$, 且 $\exists U_{\min} > 0$ 使得 $U \ge U_{\min}$ 在 $[-\|o\|_{L^\infty} - 1, \|o\|_{L^\infty} + 1]$ 上成立；核 $\gamma \in W^{2,1}(\mathbb{R}_{\ge 0})$ 单调递减且 $\int_0^\infty \gamma = 1$；障碍 $o \in W^{2,\infty}(\mathbb{R}; \mathbb{R}_{>0})$, $o' \in W^{1,1}(\mathbb{R})$, $\inf o > 0$；初值 $q_0 \in BV(\mathbb{R}; \mathbb{R}_{\ge 0})$；且 $d_o := \operatorname{ess\,inf}_{x\in\mathbb{R}} (o(x) - q_0(x)) > 0$。

Theorem 8 (a)：对每个固定 $\varepsilon > 0$，松弛问题

$$ \partial_t q_\varepsilon + \partial_x \big( V_\varepsilon(o(x) - q_\varepsilon)\, q_\varepsilon\, U(W[q_\varepsilon](t,x)) \big) = 0, \quad q_\varepsilon(0,\cdot) = q_0 $$

存在唯一熵解 $q_\varepsilon \in C([0,T]; L^1_{\mathrm{loc}}(\mathbb{R}))$，且满足 $0 \le q_\varepsilon(t,x) < o(x)$ 几乎处处成立。

Theorem 8 (b)：若 $q_0$ 还满足单侧Lipschitz下界 $\frac{q_0(x)-q_0(y)}{x-y} \ge c_{q_0}$，则存在子列 $\{q_{\varepsilon_k}\}$ 使得 $q_{\varepsilon_k} \to q$ 在 $C([0,T]; L^1_{\mathrm{loc}}(\mathbb{R}))$ 中强收敛，且 $q$ 满足障碍约束 $q(t,x) \le o(x)$，质量守恒 $\|q(t,\cdot)\|_{L^1(\mathbb{R})} = \|q_0\|_{L^1(\mathbb{R})}$，以及有限传播速度（由Theorem 6的BV估计保证）。

Theorem 9：设 $q_{\varepsilon_k} \to q$ 如上，且 $V_{\varepsilon_k}(o - q_{\varepsilon_k}) \overset{*}{\rightharpoonup} \tilde{V}$ 在 $L^\infty((0,T) \times \mathbb{R})$ 中弱*收敛，则极限 $q$ 是如下间断非局部守恒律的解：

$$ \partial_t q + \partial_x \big( \tilde{V}(t,x)\, q(t,x)\, U(W[q](t,x)) \big) = 0, \quad q(0,\cdot) = q_0, $$

其中 $\tilde{V} \in L^\infty((0,T) \times \mathbb{R}; [0,1])$，且进一步有 $\tilde{V}(t,x) = 1$ 几乎处处在 $\{(t,x) : q(t,x) \ne o(x)\}$ 上成立。这意味着在非接触区域，极限通量恢复为原始非局部通量 $q U(W[q])$；而在接触区域 $q = o$，$\tilde{V}$ 可能取 $[0,1]$ 内任意值，反映障碍激活时速度的退化不确定性。

Theorem 10：在相同假设下，罚函数复合项 $V_\varepsilon(o - q_{\nu,\varepsilon})$ 满足单侧Lipschitz上界：存在常数 $\hat{C} = \hat{C}(U_{\min}, \|U\|_{W^{2,\infty}(J)}, \|\gamma\|_{W^{2,1}}, c_{q_0}, \|o\|_{W^{2,1}})$ 使得

$$ \frac{V_\varepsilon(o(x)-q_{\nu,\varepsilon}(t,x)) - V_\varepsilon(o(y)-q_{\nu,\varepsilon}(t,y))}{x-y} \le \hat{C}, \quad \forall t \in (0,T),\, x \ne y. $$

该估计是获得 $V_\varepsilon(o - q_\varepsilon)$ 的 $L^\infty$ 弱*紧性的关键，其证明依赖于对 $V_\varepsilon'$ 的精细控制（见Methods部分）。

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

本文证明的核心障碍在于：非局部算子 $W$ 与奇异性罚函数 $V_\varepsilon$ 的耦合彻底破坏了标准双曲方程的分析工具链。作者的真正创新点在于构建了一个三层嵌套的正则化-紧性框架，每一层均针对特定障碍设计，且各层估计必须关于参数 $\nu$（粘性）和 $\varepsilon$（罚尺度）保持一致。

第一层：粘性近似与收缩映射（Theorems 2 & 3） 障碍：非局部-非线性通量 $F(q) = V_\varepsilon(o - q) q U(W[q])$ 的弗雷歇导数难以显式计算，标准不动点论证（如Banach压缩映射）失效。作者的创新在于引入加权时间变换 $r_{\nu,\varepsilon}[\lambda](t,x) = e^{-\lambda t} q_{\nu,\varepsilon}(t,x)$（Remark 2），将原粘性方程 $\partial_t q_{\nu,\varepsilon} + \partial_x F(q_{\nu,\varepsilon}) = \nu \partial_{xx} q_{\nu,\varepsilon}$ 转化为关于 $r$ 的新方程，其中非线性项 $F$ 的依赖被解耦为 $r$ 和 $\varrho$（用于构造不动点算子 $F$）。关键技巧是：在估计 $\|F[\varrho] - F[\hat{\varrho}]\|_{L^2(L^2)}$ 时，利用Young卷积不等式控制 $\|W[\varrho] - W[\hat{\varrho}]\|_{L^\infty}$，并选择足够大的 $\lambda$（依赖于 $\nu$）使 $F$ 成为压缩映射（Theorem 2）。这避免了直接处理 $W[q]$ 对 $q$ 的Gâteaux导数，是处理非局部性的核心技巧。

第二层：比较原理与障碍保持（Theorem 4） 障碍：标准最大值原理不适用于含 $W[q]$ 的方程，因为 $W[q]$ 依赖全局信息。作者的创新在于构造障碍差函数 $m_{\nu,\varepsilon}(t) = \sup_{x} (q_{\nu,\varepsilon}(t,x) - o(x))$，并利用粘性解的最优性条件（在极大值点 $x^*(t)$ 处 $\partial_x q = o'$, $\partial_{xx} q \le o''$）对 $m'_{\nu,\varepsilon}$ 进行精细估计。关键步骤是将 $\partial_x W[q]$ 控制为 $\|q\|_{L^\infty} |\gamma|_{TV}$（总变差），从而得到微分不等式

$$ m'_{\nu,\varepsilon}(t) \le V_\varepsilon(-m_{\nu,\varepsilon}(t)) \left( \|o'\|_{L^\infty} \|U_\nu\|_{L^\infty} - \|q_{\nu,\varepsilon}\|_{L^\infty}^2 \|U'_\nu\|_{L^\infty} |\gamma|_{TV} \right). $$

由于 $V_\varepsilon(-s) > 0$ 当 $s < 0$，且初始值 $m_{\nu,\varepsilon}(0) < 0$（由Assumption (D)保证），ODE比较原理直接推出 $m_{\nu,\varepsilon}(t) < 0$，即 $q_{\nu,\varepsilon} < o$ 严格成立。此论证完全规避了非局部算子的复杂性，仅依赖其 $L^\infty$-有界性。

第三层：单侧Lipschitz紧性与极限通量刻画（Theorems 5 & 10） 障碍：$V_\varepsilon'$ 的爆炸性（$\|V_\varepsilon'\|_{L^\infty} \sim \varepsilon^{-1}$）导致 $\partial_x q_{\nu,\varepsilon}$ 的经典估计失效。作者的创新在于对 $m_{\nu,\varepsilon}(t) = \inf_x \partial_x q_{\nu,\varepsilon}(t,x)$ 建立二次多项式型微分不等式（Theorem 5, Eq. (10)）：

$$ m'_{\nu,\varepsilon}(t) \ge a_1 m_{\nu,\varepsilon}(t)^2 + a_2 m_{\nu,\varepsilon}(t) + a_3, $$

其中系数 $a_i$ 由 $\|o\|_{W^{2,\infty}}$, $\|U\|_{W^{2,\infty}}$, $|\gamma|_{TV}$ 等确定，且 $a_1 = U_{\min} > 0$。此结构源于对 $\partial_{xx} q$ 和 $\partial_{xxx} q$ 的最优性条件的充分利用，并巧妙利用 $V_\varepsilon'' \le 0$ 和 $q \le o$ 的符号。由此，$m_{\nu,\varepsilon}$ 被迫大于某个负常数 $a_-$，从而获得一致的OSL下界。类似地，Theorem 10 对 $V_\varepsilon(o - q_{\nu,\varepsilon})$ 的导数进行估计，关键观察是：当 $q_{\nu,\varepsilon} \le o_\nu / 2$ 时，$V_\varepsilon'(o_\nu - q_{\nu,\varepsilon})$ 被指数衰减项 $e^{-o_{\min}/(2\varepsilon)}$ 控制，使其贡献可忽略；而在 $q_{\nu,\varepsilon} \approx o_\nu$ 区域，$o' - \partial_x q$ 很小，故 $V_\varepsilon'(o_\nu - q_{\nu,\varepsilon})(o' - \partial_x q)$ 仍可控。这体现了作者对奇异性与正则性之间精细平衡的深刻把握。

4. 与之前工作的比较

本文与既有工作的推进主要体现在三个维度：问题设定、技术假设与极限解释，且均与文中明确引用的工作（[4,38,40,47]等）形成清晰对比。

1. 与局部障碍问题 [4] 的比较：从局部到非局部的本质跃迁 Amorim et al. [4] 研究的是局部方程 $\partial_t q + \partial_x(q U(q)) = 0$ 的罚函数方法，其核心估计（如OSL）仅依赖 $U$ 和 $o$ 的局部性质。本文将此框架推广至非局部方程，代价是引入了更严格的核正则性假设：[4] 仅需 $\gamma \in L^1$，而本文要求 $\gamma \in W^{2,1}(\mathbb{R}_{\ge 0})$ 且单调递减（Assumption (K)）。这是因为 $\partial_x W[q] = -\gamma(0) q + \int_x^\infty \gamma'(y-x) q(t,y) dy$ 的估计需要 $\gamma'$ 的可积性（见Theorem 4的推导），而 $\partial_{xx} W[q]$ 的出现则要求 $\gamma'' \in L^1$（Theorem 5）。此外，[4] 的收敛极限是标准局部熵解，而本文极限是间断非局部解，其通量包含未知因子 $\tilde{V}$，反映了非局部性与奇异性耦合产生的新物理现象（如下游密度突变对上游速度的瞬时反馈）。

2. 与一般非局部方程理论 [38] 的比较：从存在性到参数一致性的深化 Very recently, existence and uniqueness of entropy solutions to the penalized equation have been shown in [38] using fixed point methods. However, [38] works in a very general setting without tracking parameter dependence. 本文的 Theorem 8 (a) 虽然也得到存在性，但其证明（基于粘性近似和Theorem 2的收缩映射）刻意构造了关于 $\varepsilon$ 一致的先验估计（如Theorem 6的BV界中 $\|V'_\varepsilon\|_{L^\infty}$ 项被显式吸收），这是 [38] 未提供的。这种一致性是进行 $\varepsilon \to 0$ 极限的先决条件，也是本文真正的技术贡献所在。换言之，[38] 解决了“是否存在解”，而本文解决了“解族是否紧致并收敛”。

3. 与非局部方程基础理论 [40,44,3] 的比较：障碍引入的正则性损失 经典非局部理论 [40,44,3] 证明：若 $q_0 \in BV$，则 $\partial_t q + \partial_x(q U(W[q])) = 0$ 的解 $q \in C([0,T]; BV)$，即BV正则性被完美保持。本文的Theorem 8 (b) 表明，在加入障碍约束后，极限解 $q$ 仍属于 $C([0,T]; L^1_{\mathrm{loc}})$，但无法保证 $q \in BV$。这是因为罚函数 $V_\varepsilon$ 的奇异性在极限中可能诱导出新的不连续性（如接触前沿的尖峰）。Theorem 9 中 $\tilde{V}$ 的弱*极限而非强极限，正是这种正则性损失的数学体现。相比之下，[47] 的局部障碍问题虽也产生不连续性，但其极限通量结构（如 $\tilde{V} = \chi_{\{q < o\}}$）更为明确；而本文的 $\tilde{V}$ 在接触集上的值未被唯一确定，这是非局部依赖带来的根本不确定性，也是未来工作（如添加微观模型或更高阶修正）的起点。

11. 围绕小 Landau 解的定常 Navier-Stokes 方程的精细渐近行为

Refined asymptotics of the steady Navier Stokes equation around small Landau solutions

作者: Hao Jia, Vladimír Šverák

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.24200v1

摘要翻译: 本文研究三维 Navier-Stokes 方程在外部区域上小定常解的大距离渐近行为。Korolev 与第二作者 \cite{SverakKorolev} 已证明其主导项由 Landau 解给出，并猜想下一项应为 $O(1/|x|^2)$（当 $x\to\infty$）。我们证实该猜想成立，并将下一项渐近式用某个在单位球面上围绕 Landau 解构造的线性化算子的特征值显式表达出来。其中部分项的衰减阶恰好为 $O(1/|x|^2)$，而其余项的衰减则略有加速。

Comment: 50 pages, comments welcome

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

1. 研究对象和背景

本文研究三维外区域 $\mathbb{R}^3 \setminus B_1$ 上小定常 Navier-Stokes 方程解在无穷远处的精细渐近结构。所研究方程为

$$ -\Delta u + u \cdot \nabla u + \nabla p = 0, \quad \mathrm{div}\, u = 0, \qquad x \in \mathbb{R}^3 \setminus B_1, \tag{1.1} $$

其中 $B_1 = \{x \in \mathbb{R}^3 : |x| < 1\}$，速度场 $u : \mathbb{R}^3 \setminus B_1 \to \mathbb{R}^3$ 与压力 $p : \mathbb{R}^3 \setminus B_1 \to \mathbb{R}$ 满足衰减条件 $|u(x)| = O(|x|^{-1})$、$|p(x)| = O(|x|^{-2})$，并附加自然无通量条件

$$ \int_{\partial B_R} u \cdot n \, d\sigma = 0, \quad \forall R > 1. \tag{1.2} $$

该问题的核心对象是 Landau 解族 $U_b$, $P_b$（$b \in \mathbb{R}^3$），其定义于全空间 $\mathbb{R}^3$ 上，满足分布意义下方程

$$ -\Delta U_b + U_b \cdot \nabla U_b + \nabla P_b = b\, \delta(x), \quad \mathrm{div}\, U_b = 0, \tag{1.4} $$

并具有 $-1$-齐次性：$U_b(\lambda x) = \lambda^{-1} U_b(x)$。在球坐标 $(r,\theta,\phi)$ 下，其显式形式为

$$ U_b(x) = \frac{1}{|x|} V^b\left(\frac{x}{|x|}\right) e_\theta + \frac{1}{|x|} F^b\left(\frac{x}{|x|}\right) e_r, \tag{1.7} $$

其中 $V^b(\omega) = -2\epsilon \sin\theta / (1 - \epsilon \cos\theta)$，$F^b(\omega) = 2[(1-\epsilon^2)/(1-\epsilon \cos\theta)^2 - 1]$，且 $|b| = f(\epsilon) = 16\pi[\epsilon^{-1} + \frac{1}{2}\epsilon^{-2} \log\frac{1-\epsilon}{1+\epsilon} + \frac{4\epsilon}{3(1-\epsilon^2)}]$，$\epsilon \in (0,1)$。Landau 解是唯一光滑、尺度不变的 $\mathbb{R}^3 \setminus \{0\}$ 上定常 Navier-Stokes 解（见 [21]）。Korolev–Šverák [14] 已证明：对满足 $|u(x)| \le \delta_0 / |x|$ 的小解，存在 $b \in \mathbb{R}^3$（由通量张量 $T_{ij}(x) = -(\partial_{x_j} u_i + \partial_{x_i} u_j) + u_i u_j + p \delta_{ij}$ 在球面上积分定义）使得

$$ |u(x) - U_b(x)| = O(|x|^{-\gamma}), \quad \gamma \in (1,2). \tag{1.13} $$

然而，该估计中 $\gamma < 2$ 是严格成立的；端点 $\gamma = 2$ 被排除，且当 $\gamma \to 2^-$ 时，允许的初值大小 $\delta_0$ 收缩至零。文献 [14] 明确提出猜想：$\gamma = 2$ 应可实现，即次主导项应为 $O(|x|^{-2})$。但标准扰动论失效——因非线性耦合 $U_b$ 与 $O(|x|^{-2})$ 项生成 $O(|x|^{-4})$ 项，其在 Stokes 算子逆作用下产生对数发散积分

$$ D(x) = \int_{|y|>1} \frac{1}{|x-y|^2} \frac{1}{|y|^3} \, dy \sim |x|^{-2} \log|x|, $$

这正是 $\gamma = 2$ 障碍的根源。因此，本文试图填补的关键空白在于：严格证明 $\gamma = 2$ 可达，并精确刻画 $O(|x|^{-2})$ 项的结构——它并非任意 $C^\infty(S^2)$ 函数，而是受 Landau 解线性化谱约束的特定子空间上的球面函数，且其衰减速率可能被进一步加速（如 $|x|^{-(\alpha_j+1)}$ 形式，$\alpha_j > 1$）。该问题的重要性在于：它揭示了非线性对远场渐近的深层影响——即使在小数据情形，非线性亦不可视为微扰，而必须通过谱分析在临界尺度上进行精细分解。

2. 主要定理或主要结果

定理 1.1：存在 $\delta_0 \in (0,1)$，使得对任意满足 $|u(x)| \le \delta_0 / |x|$ 及 (1.2) 的光滑解 $(u,p)$，令 $b \in \mathbb{R}^3$ 由 (1.10)–(1.11) 定义，$\epsilon$ 由 (1.8) 确定，则存在一个八维空间 $E_b \subset C^\infty(S^2; \mathbb{R}^3)$ 及其正交直和分解 $E_b = E^1_b \oplus E^2_b \oplus E^3_b$，其中 $\dim E^1_b = 4$，$\dim E^2_b = \dim E^3_b = 2$。进而有如下渐近展开：

$$ u(x) = U_b(x) + \frac{\Phi_1(x/|x|)}{|x|^2} + \frac{1}{|x|^{\alpha_2+1}} \left[ \cos(\beta_2 \log |x|) \Phi_{2,1}(x/|x|) + \sin(\beta_2 \log |x|) \Phi_{2,2}(x/|x|) \right] + \frac{1}{|x|^{\alpha_3+1}} \left[ \cos(\beta_3 \log |x|) \Phi_{3,1}(x/|x|) + \sin(\beta_3 \log |x|) \Phi_{3,2}(x/|x|) \right] + O(|x|^{-5/2}), \tag{1.16} $$

其中 $\Phi_1 \in E^1_b$，$\{\Phi_{2,1}, \Phi_{2,2}\} \subset E^2_b$，$\{\Phi_{3,1}, \Phi_{3,2}\} \subset E^3_b$，且参数满足

$$ \alpha_2 = 1 + \frac{1}{15\epsilon^2} + O(\epsilon^3), \quad \beta_2 = O(\epsilon^3), \quad \alpha_3 = 1 + \frac{4}{15\epsilon^2} + O(\epsilon^3), \quad \beta_3 = O(\epsilon^3). \tag{1.17} $$

特别地，$\alpha_2, \alpha_3 > 1$，故后两项衰减快于 $|x|^{-2}$（即 $|x|^{-(\alpha_j+1)} = o(|x|^{-2})$）。注 1.2 进一步指出：(i) 展开式中所有项均以 $|x|^{-2}$ 或更快速率衰减，这与线性 Stokes 方程的 $|x|^{-2}$ 项形成对比，表明 Landau 解的非线性效应加速了次主导项衰减；(ii) 数值计算强烈暗示 $\beta_2 = \beta_3 = 0$，即无真实振荡项，但严格证明尚未完成；(iii) 余项 $O(|x|^{-5/2})$ 可改进为 $O(|x|^{-\mu})$，其中 $\mu < 3$ 且 $\mu \to 3^-$ 当 $|b| \to 0^+$，但 $\mu = 3$ 可能因 $b=0$ 时潜在的对数项而不可达。

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

核心障碍在于：标准扰动框架无法处理 $U_b$ 与 $O(|x|^{-2})$ 项的非线性共振，其导致的 $|x|^{-2} \log|x|$ 型误差破坏 $\gamma = 2$ 的可能性。作者的根本创新是将远场渐近问题完全重铸为球面上线性化算子的谱问题。具体而言：(1) 引入球坐标 $(r,\omega) \in (1,\infty) \times S^2$，并设 $s = -\log r$，将原方程转化为关于 $s$ 的演化系统（见 (2.6)–(2.9)）。定义六分量态向量 $\mathbf{X} = (\xi, \xi', \vartheta, \vartheta^*)^\top$，其中 $\xi$ 为切向速度分量，$\vartheta$ 为径向分量，$\xi', \vartheta^*$ 为其对数导数。则线性化系统写为 $\partial_s \mathbf{X} = L \mathbf{X} + \mathcal{N}(\mathbf{X})$，其中非线性项 $\mathcal{N}$ 在 $\mathbf{X} = O(e^{-s})$ 时为高阶小量。(2) 关键算子 $L$（定义于 (2.10)）及其未扰动对应物 $L_0$（定义于 (2.11)）作用于球面函数空间 $D$（见 (2.12)），其谱决定渐近行为。$L_0$ 对应线性 Stokes 系统，其特征值为 $\lambda = 0$（四重）和 $\lambda = 1$（八重），后者源于 $1/|x|^2$ 项的球谐展开自由度。(3) 精细谱分析聚焦于 $\lambda = 1$ 的分裂：利用 $U_b$ 的轴对称性，对傅里叶模 $m = 0,1,2$ 分别处理（见第 6–7 节）。证明 $\lambda = 1$ 在扰动下保持为四重特征值（对应平移对称性及一个非平凡‘纯涡旋’模），而剩余四维空间分裂为两个二重特征值簇，其本征值 $\lambda_{2,3} = -\alpha_{2,3} \pm i \beta_{2,3}$ 直接给出 (1.16) 中的指数与振荡频率（因 $e^{\lambda s} = e^{-\alpha s} e^{\pm i \beta s} = |x|^{\alpha} \cos(\beta \log|x|)$）。(4) 计算扰动特征值的核心工具是‘不变子空间图表示法’：将扰动后的特征子空间表示为未扰动子空间 $\mathrm{ker}(L_0 - \lambda I)$ 上的图 $\Gamma : \mathrm{ker}(L_0 - \lambda I) \to \mathrm{ran}(L_0 - \lambda I)$，再将特征值问题约化为有限维矩阵摄动计算（类似量子力学中的非简并/简并微扰论），从而获得 (1.17) 的渐近展开。此方法绕开了对数发散的直接处理，将问题彻底代数化。

4. 与之前工作的比较

与 Korolev–Šverák [14] 相比，本文实现了根本性推进：[14] 仅得到 $|u - U_b| = O(|x|^{-\gamma})$（$\gamma < 2$），且其证明依赖于对 $U_b$ 的粗略控制和能量估计，无法触及 $\gamma = 2$ 的临界情形；本文则通过谱重铸，严格证明了 $\gamma = 2$ 可达，并给出了精确到 $O(|x|^{-5/2})$ 的完整渐近级数。与线性 Stokes 理论（即 $U_b \equiv 0$）相比，本文揭示了非线性导致的本质差异：线性理论中 $1/|x|^2$ 项由八维空间 $E_0$ 参数化，而本文中该空间分裂为三个子空间 $E^1_b, E^2_b, E^3_b$，且其中两部分对应加速衰减项（$|x|^{-(\alpha_j+1)}$，$\alpha_j > 1$），这是纯线性框架完全无法预测的。与以往基于势论或多重极展开的工作（如 [1, 2, 4, 6, 7, 8, 9, 10, 12, 15, 17, 18, 19]）相比，本文不依赖于特定坐标系下的分离变量假设，而是建立了一个普适的、基于球面谱的分析框架，适用于任意小 $b$。与 [21]（刻画 Landau 解的唯一性）相比，本文不关注解的存在唯一性，而专注于解的精细结构，其工具（$s = -\log r$ 演化、算子 $L$ 的谱分析）是全新的。本文的代价在于：结果仅适用于小解（$\delta_0$ 依赖于 $b$），且对 $b=0$ 的极限情形，$\mu = 3$ 的余项尚不能确认（可能存在 $|x|^{-3} \log|x|$ 项）；此外，振荡频率 $\beta_2, \beta_3$ 的精确零值尚未证明，数值证据虽强，但理论缺口仍存。最后，本文明确指出其方法可推广至其他具临界齐次性的非线性椭圆/稳态问题，只要能构造出相应的球面线性化算子并分析其谱分裂。

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

1. Schrödinger 算子在 Bessel 位势空间中的逐点收敛性

Pointwise Convergence of Schrödinger Operators in Bessel Potential Spaces

作者: Yucheng Pan, Wenchang Sun, Jiheng Tan

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25833v1

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

摘要翻译

中文标题: Schrödinger 算子在 Bessel 位势空间中的逐点收敛性

摘要翻译: 我们研究初值属于 Bessel 位势空间 $L_s^p(\mathbb{R}^n)$ 的自由 Schrödinger 方程解的逐点收敛性。我们建立了覆盖全范围 $1 \leq p < \infty$ 的逐点收敛性的新充分正则性指标，并通过反例证明：在一维情形下，这些指标对所有 $1 \leq p \leq 2$ 均是 sharp 的；在高维情形下，这些指标在端点 $p=1$ 及 $p>2+p_n$ 处亦是 sharp 的。证明依赖于高维平稳相位法。

2. Orlicz 空间中高阶椭圆系统解的内部先验估计

Interior a priori estimate for higher order elliptic systems in Orlicz spaces

作者: Amiran Gogatishvili, Pia Salerno, Lubomira Softova

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25780v1

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

摘要翻译

中文标题: Orlicz 空间中高阶椭圆系统解的内部先验估计

摘要翻译: 我们研究具有可变 Calderón–Zygmund 核的奇异积分算子及其与 $VMO$ 函数的交换子在 Orlicz 空间框架下的性质。在回顾经典 $L^p$ 理论之后，我们在 Young 函数满足标准 $Δ_2$ 与 $ abla_2$ 条件的前提下，建立了其在 $L^Φ$ 中的有界性结果。证明依赖于分解技巧与弱型估计。作为应用，这些结果为具不连续系数的高阶椭圆算子解的先验估计与内部正则性提供了泛函分析基础。

3. Bochner空间中Laplace方程的非局部问题

Nonlocal problem for Laplace equation in Bochner spaces

作者: Bilal Bilalov, Sabina Sadigova, Pia Salerno, Lubomira Softova

类别: math.AP, math.FA

发布日期: 2026-05-25

arXiv ID: 2605.25761v1

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

摘要翻译

中文标题: Bochner空间中Laplace方程的非局部问题

摘要翻译: 我们在无界矩形区域 $Π= I \times (0,\infty)$ 上研究Laplace方程，其中 $I= (0,2π)$，并施加在边界 $\partial Π$ 上以迹意义理解的非局部边界条件。分析在Bochner-Sobolev空间 $W^2_{p,1}(Π;X)$ 中进行，该空间与Bochner空间 $L^{p,1}(Π;X)$ 相关联，其中 $ p \in (1,\infty)$ 且 $X$ 是一个适当的Banach空间。为求解该问题，我们采用广义谱方法。特别地，我们引入由张量积生成的 $\otimes$-基的概念，并将经典标量情形中已知的谱方法框架推广至当前向量值设定。 此外，我们证明：相应非局部谱问题的根函数系在 $L^p(I;X)$ 中构成一个 $\otimes$-基。

4. 关于一类带权的尖锐Sobolev型估计

On a class of sharp Sobolev type estimates with weights

作者: Raul Hindov, Evgeniy Lokharu

类别: math.AP, math-ph

发布日期: 2026-05-25

arXiv ID: 2605.25637v1

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

摘要翻译

中文标题: 关于一类带权的尖锐Sobolev型估计

摘要翻译: 我们研究如下形式的尖锐带权Sobolev型不等式：

$$ \int_{0}^{1}|u(x)|ρ(x) \diff x \leqslant Λ \Bigl(\int_{0}^{1}|u^{(k)}(x)|^2 \diff x\Bigr)^{1/2}, \qquad u\in H_0^k(0,1), $$

其中 $ρ$ 是一个非负权函数。我们刻画了极小化子，并通过证明每个极小化子具有恒定符号、从而满足一类多调和型非线性特征值问题，确定了最优常数 $Λ(k,ρ)$。这给出了大类权函数下极值函数的显式刻画。此外，我们甚至给出了最优常数关于权函数 $ρ$ 的显式计算公式。这些新的带权估计非常有用；除其他应用外，它们使我们得以复现若干此前已知的尖锐估计以及有限区间上的Hardy型不等式。

5. 加权Weyl-Sonine算子的半群方法：绕过Bernstein函数实现适定性

A Semigroup Approach to Weighted Weyl-Sonine Operators: Bypassing Bernstein Functions for Well-Posedness

作者: Gustavo Dorrego

类别: math.FA, math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25630v1

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

摘要翻译

中文标题: 加权Weyl-Sonine算子的半群方法：绕过Bernstein函数实现适定性

摘要翻译: 本文为半无限区间上的加权Weyl-Sonine分数阶算子建立了严格的泛函分析框架。尽管经典的Phillips函数演算严格依赖于完全单调的Bernstein函数来对线性半群进行从属化（subordination），我们证明这一结构性限制可被完全绕过。通过利用容许的失忆权（$ω$）所诱导的内蕴拓扑衰减，我们证明抽象积分表示对一般类别的缓增Sonine核绝对收敛，从而在一个统一框架下同时涵盖扩散型松弛与振荡型非局部动力学。借助强连续半群理论，我们将这些算子刻画为一类形变平移半群的无穷小生成元，该半群同时纳入主观老化尺度（$ψ$）与历史衰减型失忆效应。这一基于相似性的方法不仅为这些算子的超奇异Marchaud表示提供了坚实的抽象基础，而且保证了其闭性以及相应抽象分数阶Cauchy问题在一般Banach空间中的解析适定性（analytical well-posedness）。此外，我们引入相应的广义分数阶积分，并证明基本反演定理，揭示该理论的完整解析对称性。最后，我们通过加权Sobolev嵌入建立解的逐点衰减估计。

6. Grushin 算子的统一热核渐近性

Uniform heat kernel asymptotics for the Grushin operator

作者: Yimeng Chen, Hong-Quan Li, Jun-Cheng Tang, Jia-Yu Yang

类别: math.AP, math.CA

发布日期: 2026-05-25

arXiv ID: 2605.25622v1

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

摘要翻译

中文标题: Grushin 算子的统一热核渐近性

摘要翻译: 采用文献 \cite{li2021carnotcaratheodory, LZ2025} 中引入的强有力方法，我们研究 Grushin 算子 $Δ_G = Δ_x + |x|^2 Δ_u$ 在 $ \mathbb{R}_x ^n \times \mathbb{R}_u ^{n'} $ 上对应的热核在无穷远处的渐近行为，其中 $ n , \, n' \ge 1 $。我们进一步建立了其空间导数的尖锐估计。作为副产品，我们将给出 Grushin 热核的精确估计及其小时间渐近行为，其显式程度达到所能期望的最高水平。

7. 非交换环面上的高阶热估计与半线性热方程

Higher-Order Heat Estimates and Semilinear Heat Equations on the Noncommutative Torus

作者: Fulin Yang, Zhipeng Yang

类别: math.AP, math.FA

发布日期: 2026-05-25

arXiv ID: 2605.25513v1

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

摘要翻译

中文标题: 非交换环面上的高阶热估计与半线性热方程

摘要翻译: 我们在非交换环面 $\mathbb{T}_θ^{n}$ 上建立了具有完全界控的尖锐高阶热估计，并证明了其在小时间尺度下的阶最优性。作为其在 $\mathbb{T}_θ^{n}$ 上多项式半线性热方程中的应用，我们给出了局部适定性、爆破替代准则、高阶正则性的持续性，以及 Sobolev 代数尺度 $H^k_θ$（其中 $k>n/2$）上的瞬时光滑性。

8. 关于能量临界非均匀Hartree方程适定性的若干注记

Remarks on the well-posedness of the energy-critical inhomogeneous Hartree equation

作者: Seongyeon Kim, Tarek Saanouni, Ihyeok Seo

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25341v1

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

摘要翻译

中文标题: 关于能量临界非均匀Hartree方程适定性的若干注记

摘要翻译: 我们研究三维及更高维空间中的能量临界非均匀Hartree方程。此前的局部适定性结果未能覆盖非均匀性指数较小、且Riesz位势指数既可较小又可较大的参数区域。我们建立了新参数范围下的局部适定性，从而实质性地填补了剩余的未解决问题的参数区域。特别地，我们的结果完全解决了五维和六维情形中所遗留的间隙。

9. Benjamin-Ono 方程的孤子分解猜想：显式的 $L^\infty$ 渐近误差公式

Soliton resolution conjecture for the Benjamin-Ono equation: Explicit $L^\infty$ asymptotic error formula

作者: Hong-Yu Pan, Shou-Fu Tian

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25327v1

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

摘要翻译

中文标题: Benjamin-Ono 方程的孤子分解猜想：显式的 $L^\infty$ 渐近误差公式

摘要翻译: 我们证明了 Benjamin-Ono (BO) 方程的孤子分解猜想，并在 $L^\infty$-范数下给出显式误差界。对于有限阶多孤子情形，当初始数据 $u_0 \in H^{s,\alpha}(\mathbb{R})$ 且任意 $s>1/2$、$\alpha\geqslant 1$ 时，显式的 $L^\infty$-范数误差被控制在 $\mathcal{O}(|t|^{-\frac{1}{4}(1-\frac{1}{2s})})$。对于无限阶多孤子情形，当 $u_0$ 可表为孤子剖面的无穷级数时，显式的 $L^\infty$-范数误差被控制在 $\mathcal{O}(|t|^{-1/3})$。最近，Gassot、Gérard 与 Miller（arXiv:2601.10488, 2026）在 $u_0 \in H^{1,1}\left( \mathbb{R} \right)$ 的有限阶多孤子情形下，证明了孤子分解在 $H^1$-范数下的隐式误差界，但需额外条件 $x^2u_0(x) = c_0 + v_0(x)$，其中 $c_0\in \mathbb{R}$、$v_0(x) \in L^2(\mathbb{R})$。在无限阶多孤子情形下，Gassot 与 Gérard（arXiv:2603.15419, 2026）证明了当 $u_0$ 可表为孤子剖面的无穷级数时，孤子分解在 $L^\infty$-范数下的隐式误差界。值得注意的是，他们指出 Lax 算子的反谱问题是一个有趣的公开问题。为了建立有限阶/无限阶多孤子情形下具有显式 $L^\infty$-范数误差的孤子分解，目前仍存在诸多关于初始条件、误差精度及其他相关问题的公开问题。解决这些公开问题是本文工作的核心目标。为扩大初始数据空间并去除额外条件，我们运用 Kato-Rellich 定理将孤子分解猜想转化为序列与解之间的误差估计问题。值得指出的是，我们通过基于离散谱构造一个迹类算子，解决了 Lax 算子的公开反谱问题。

10. 任意维空间中具有非捕获特征的首阶双曲方程的反推控制

Backstepping Control of First-Order Hyperbolic Equations in Arbitrary Dimensions with Non-Trapping Characteristics

作者: Mohamed Camil Belhadjoudja

类别: eess.SY, math.AP

发布日期: 2026-05-24

arXiv ID: 2605.25217v1

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

摘要翻译

中文标题: 任意维空间中具有非捕获特征的首阶双曲方程的反推控制

摘要翻译: 本文提出了一种针对定义在任意维数区域上的、具空间变系数的首阶双曲方程的边界控制反推方法。该方法基于由时不变输运算子 $\partial_t + \mathbf{a}(x) \cdot \nabla$ 的特征流所诱导的变量变换，将原多维系统转化为沿各条特征曲线演化的连续统形式的一维双曲方程组，各成员方程彼此解耦。随后，对分解所得的每个一维方程分别设计反推控制器，并将所得控制律在原始坐标系下重新组合，从而实现整个系统的有限时间稳定化。该框架依赖于空间区域上存在一族叶状分布的特征曲线，且其穿越时间一致有界（即非捕获条件）

11. 带超线性扰动的混合局部-非局部Brezis-Nirenberg问题：紧性与应用

Brezis-Nirenberg problems for mixed local-nonlocal operators with superlinear perturbations: compactness and applications

作者: Mousomi Bhakta, Nirjan Biswas, Paramananda Das

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.25157v1

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

摘要翻译

中文标题: 带超线性扰动的混合局部-非局部Brezis-Nirenberg问题：紧性与应用

摘要翻译: 本文考虑如下混合局部-非局部Brezis-Nirenberg问题\begin{equation}\label{crit_pro_abstract}\tag{$\mathcal{P}_{2^*}$} -Δu+(-Δ)^s u=λ|u|^{p-2}u+|u|^{2^*-2}u\text{ in }Ω,\quad u=0\text{ in }\mathbb{R}^N \setminus Ω, \end{equation}其中$Ω\subset\mathbb{R}^N$为有界区域，$N\geq3$，$s\in(0,1)$，$λ>0$，且$2\leq p<2^*=\frac{2N}{N-2}$。我们为如下子临界/临界问题族建立紧性结果\begin{equation}\label{sub_pro_abstract}\tag{$\mathcal{P}_{p_n}$} -Δu+(-Δ)^s u=λ|u|^{p-2}u+|u|^{p_n-2}u\text{ in }Ω,\quad u=0\text{ in }\mathbb{R}^N \setminus Ω, \end{equation}其中$p_n \in (p,2^* ]$且$p_n\to 2^*$。具体而言，当$N>6-4s$时对$p \in (2+\frac{4s}{N-2},2^*)$，以及当$N\leq6-4s$时对$p \in (2^*-1,2^*)$，我们证明：问题\eqref{sub_pro_abstract}的任意有界解序列$\{u_n\}$在能量空间中相对紧，且强收敛于问题\eqref{crit_pro_abstract}的一个非平凡解。这是首篇针对非齐次算子建立此类紧性结果的工作。由于非齐次算子的存在，我们的证明需对Devillanova与Solimini（Adv. Differential Equations, 2002）以及Yan、Yang与Yu（J. Funct. Anal., 2015）所发展的方法进行非平凡的适应性改进。作为该紧性结果的应用，在相同的$N$与$p$取值范围内，我们证明问题\eqref{crit_pro_abstract}存在无穷多个变号解。我们预期该方法论可推广至更广泛的一类相关问题。

12. SU(1,1) 非线性傅里叶变换的逐点行为

Pointwise behavior of SU(1,1) nonlinear Fourier transform

作者: Sergey A. Denisov

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

发布日期: 2026-05-24

arXiv ID: 2605.25108v1

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

摘要翻译

中文标题: SU(1,1) 非线性傅里叶变换的逐点行为

摘要翻译: 我们证明：SU(1,1) 非线性傅里叶变换（NLFT）对平方可和系数可能逐点发散。作为推论，我们证明：单位圆上正交多项式的经典逐点渐近行为在 Szegö 类测度上可能失效。我们还讨论了逐点收敛成立的一些特殊情况。

13. 具有变号系数的多层球中 Helmholtz 方程的极限吸收原理

Limiting Absorption Principle for the Helmholtz Equation with Sign-Changing Coefficients in Multilayer Spheres

作者: Wenjing Zhang, Yixian Gao

类别: math.AP

发布日期: 2026-05-24

arXiv ID: 2605.24826v1

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

摘要翻译

中文标题: 具有变号系数的多层球中 Helmholtz 方程的极限吸收原理

摘要翻译: 本文研究 $\mathbb{R}^d$ ($d \ge 2$) 中一类多层 Helmholtz 模型，其特征为具有交替正负折射率的同心层状材料。为克服变号材料参数所导致的强制性（coercivity）丧失，我们构造了一类特制的 $\mathbb T$-强制性算子以恢复问题的强制结构。此外，为应对无界区域上固有的紧性缺失，我们将一个复波数 Dirichlet-to-Neumann (DtN) 算子嵌入该框架。通过将这一变分合成与尖锐的先验估计 $\textit{a priori}$ estimates 相结合，我们严格建立了极限吸收原理，并在适当函数空间中证明了相应传输问题的适定性（well-posedness）。关键地，我们通过显式分析最优迹常数（optimal trace constant），量化了唯一性对区域几何的依赖关系，从而为多层超材料的设计提供了严格的数学判据。

14. 线性动力学方程的自相似性与扩散极限：Wild 级数方法

Self-similarity and diffusive limits for linear kinetic equations: a Wild sum approach

作者: José A. Cañizo, Stéphane Mischler, Niccolò Tassi

类别: math.AP, math.PR

发布日期: 2026-05-23

arXiv ID: 2605.24685v1

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

摘要翻译

中文标题: 线性动力学方程的自相似性与扩散极限：Wild 级数方法

摘要翻译: 我们证明，在全空间中且无约束机制的线性碰撞型动力学方程，其解在长时间尺度下呈现自相似行为。 这一结果极大地改进了该背景下近期已知的衰减估计结果，并首次确立了此类自相似行为的存在性。作为推论，我们还建立了适当缩放后的解关于时间的一致收敛性——即一致地收敛至其扩散极限，该结果亦属首次。 所考虑的方程类包括若干 BGK 型方程、若干动力学非局部 Fokker–Planck 型方程，以及若干动力学（可能为分数阶）Fokker–Planck 方程；对这些方程，我们或能通过 Wild 级数（亦称 Dyson 级数）显式写出解，或能在 Fourier 侧进行精确计算。

15. 散度型椭圆方程带低阶项的Schauder–Orlicz型估计

Schauder–Orlicz-Type Estimates for Divergence-Form Elliptic Equations with Lower-Order Terms

作者: Jaouad Bourabiaa, Youssef Elmadani, Abdelouahab Hanine

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24596v1

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

摘要翻译

中文标题: 散度型椭圆方程带低阶项的Schauder–Orlicz型估计

摘要翻译: 针对具低阶项的二阶线性散度型椭圆方程的弱解，推导了Schauder Orlicz型估计。首先处理Orlicz空间 $X=L^ψ$ 情形：在Young函数 $ψ$ 及方程系数的适当假设下，确定了低阶项数据的最优相伴空间；进而得到 $W^{1,ψ}$ 中的先验估计。随后将讨论推广至重排不变Banach函数空间；为此引入一类空间 $(\mathcal C)$，用以刻画所有满足存在相应相伴空间 $Y$ 从而可获得Schauder型估计的空间 $X$；最后以Lorentz空间作为具体实例加以考察。

16. 多重极Aharonov–Bohm哈密顿量下薛定谔方程的Strichartz估计

Strichartz estimates for Schrödinger equations with the multipole Aharonov–Bohm Hamiltonian

作者: Mengxuan Yang, Junyong Zhang

类别: math.AP, math-ph

发布日期: 2026-05-23

arXiv ID: 2605.24361v1

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

摘要翻译

中文标题: 多重极Aharonov–Bohm哈密顿量下薛定谔方程的Strichartz估计

摘要翻译: 我们在 $\mathbb{R}^2$ 上证明了具有多重极Aharonov–Bohm哈密顿量的薛定谔方程的整体时间Strichartz估计。作为中间步骤，我们证明了多重极Aharonov–Bohm哈密顿量的整体时间局部光滑性估计。

17. Boltzmann方程朝向含激波的通用Riemann解的流体力学极限

Hydrodynamic Limit of the Boltzmann Equation toward Generic Riemann Solutions with Shocks

作者: Mingi Choe, Moon-jin Kang, Chanwoo Kim

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24392v1

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

摘要翻译

中文标题: Boltzmann方程朝向含激波的通用Riemann解的流体力学极限

摘要翻译: 我们建立了具硬球碰撞的一维Boltzmann方程朝向可压Euler系统Riemann解的流体力学极限。本结果所覆盖的Riemann解包括通用的基本波叠加：或为两道激波与一道接触间断，或为一道稀疏波、一道接触间断与一道激波。对适当预处理的初始数据及足够小的波强度，我们证明相应的Boltzmann解全局存在，并且当Knudsen数趋于零时，在 $L^2([0,T]\times\mathbb R_x\times\mathbb R^3_ξ)$ 中收敛于与该Riemann解相关联的局域Maxwellian分布，其中 $T>0$ 任意。激波位置由动力学未知量——即所谓“位移”（Shifts）——进行调制，这些位移被构造为 $[0,T]$ 上的有界变差（BV）函数。本结果的一个显著特点是：该收敛性证明未去除激波层亦未去除初始层。在单激波的特殊情形下，我们的分析给出了直至动力学选取之位移为止的动能激波层的sharp定量刻画。证明结合了宏观-微观分解、激波的 $a$-contraction 方法的动能适配、层分析以及关于位移的紧性论证。

18. 四维空间中具有大指数的双调和问题解的精确渐近估计与非退化性

Precise asymptotic estimates and non-degeneracy of solutions to a biharmonic problem with large exponents in dimension four

作者: Xiuda Liang, Wenjie Wang

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24338v1

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

摘要翻译

中文标题: 四维空间中具有大指数的双调和问题解的精确渐近估计与非退化性

摘要翻译: 我们研究定义在光滑有界区域 $Ω\subset \mathbb{R}^4$ 上、满足 Dirichlet 边界条件的半线性双调和问题：\begin{equation*} \begin{cases} Δ^2 u=(u^+)^{p} &{\text{in}~Ω},\[0.5mm] u \not\equiv 0 &{\text{in}Ω},\[0.5mm] u=\partial u / \partial ν= 0 &{\text{on}\partial Ω}, \end{cases} \end{equation*} 其中 $p>1$ 充分大。 该问题解的基本渐近行为与集中现象已在文献中建立。本文旨在改进若干已知的解的渐近估计，使其更为显式，从而证明一般区域上多峰解的非退化性。主要方法包括常微分方程理论、爆破分析、局部 Pohozaev 恒等式，以及 Green 函数与 Green 表示公式的应用。

19. 利用 Hirota 方法研究可积离散聚焦 Manakov 系统中的孤子与呼吸子相互作用

Soliton and breather interactions in the integrable discrete focusing Manakov system via Hirota’s method

作者: Uyen Le, Alexander Chernyavsky, Barbara Prinari

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

发布日期: 2026-05-23

arXiv ID: 2605.24333v1

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

摘要翻译

中文标题: 利用 Hirota 方法研究可积离散聚焦 Manakov 系统中的孤子与呼吸子相互作用

摘要翻译: 本文旨在将 Hirota 双线性方法应用于聚焦色散情形下的可积离散 Manakov 系统，以构造并分析孤子与呼吸子解。在导出该系统的通用双线性形式后，我们展示了如何获得基本孤子解，以及基本呼吸子与复合呼吸子解。随后，我们得到包含 2 个孤子与 2 个呼吸子的解，以及孤子与呼吸子组合的解，并系统讨论所有“两体”相互作用性质，特别强调显式公式、可视化图像与长时间渐近行为，从而严格证实了这些相干结构所具有的高度非平凡相互作用性质。

20. 双曲守恒律的变分公式

Variational formulation of hyperbolic conservation laws

作者: Eitan Tadmor

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.24214v1

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

摘要翻译

中文标题: 双曲守恒律的变分公式

摘要翻译: 熵函数在双曲守恒律数学理论的发展中起着关键作用。熵的概念与对称性密切相关，是对非线性守恒律系统所施加的一种扩展。在此背景下，Friedrichs 提出了一个问题：所假设的对称性是否也可被导出？我们引入一种变分公式以回应 Friedrichs 的问题：熵函数被导出为一个极值对象，由此我们推导出（而非施加）Dafermos 的最大熵产生原理。

不相关论文（基本信息）

1. 随机环境中可逆扩散过程的定量爱因斯坦关系

Quantitative Einstein relation for reversible diffusions in a random environment

作者: Ahmed Bou-Rabee, Ruizhe Xu

类别: math.PR, math.AP

发布日期: 2026-05-25

arXiv ID: 2605.26082v1

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

摘要: The Einstein relation describes the response of a diffusing particle to a small constant external force. It states that, as the force tends to zero, the ratio of the limiting velocity to the force magnitude converges to the diffusivity matrix of the unforced particle, evaluated in the force direction. Gantert, Mathieu, and Piatnitski (2012) proved this identity for reversible diffusions in random environments. We prove a quantitative version, with an explicit quenched algebraic rate.

2. 从Bogoliubov-de Gennes到Ginzburg-Landau：非磁情形下临界温度 $T_{ m c}$ 附近的临界点

From Bogoliubov-de Gennes to Ginzburg-Landau: Critical Points Near $T_{\rm c}$ in the Non-Magnetic Case

作者: Rupert L. Frank, Christian Hainzl, Dong Hao Ou Yang

类别: math.AP, math-ph

发布日期: 2026-05-25

arXiv ID: 2605.26008v1

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

摘要: We study the relation between the Bogoliubov-de Gennes equation and the Ginzburg-Landau equation for a BCS model without external fields. While previous rigorous derivations of Ginzburg-Landau theory from BCS theory have focused on energies and minimizers, here we consider arbitrary critical points in the relevant energy regime. For temperatures close to the critical temperature, we prove that every sufficiently small solution of the BdG equation admits an asymptotic factorization into a microscopic Cooper-pair profile and a macroscopic order parameter. The latter satisfies the Ginzburg-Landau equation up to an error that vanishes in the scaling limit. Our analysis relies on a Birman-Schwinger reformulation of the BdG equation, a Lyapunov-Schmidt type reduction, and semiclassical estimates at low regularity.

3. 具有给定质量的由局部-非局部算子驱动的拟线性临界Choquard方程解的多重性

Multiplicity of solutions with prescribed mass for a quasilinear critical Choquard equation driven by a local-nonlocal operator

作者: J. Giacomoni, Nidhi Nidhi, K. Sreenadh

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25787v1

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

摘要: In this paper we study the normalized solutions of the following critical growth Choquard equation with mixed local and non-local operators: \begin{equation*} \begin{array}{rcl} -Δ_p u +(-Δ_p)^s u & = & λ|u|^{p-2}u +μ|u|^{q-2}u +(I_α*|u|^{p^_α})|u|^{p^_α-2}u \text{ in } \mathbb{R}^N; \left| u \right|_p & = & τ. \end{array} \end{equation*} Here, $N\geq 3$, $2 \le p<N$, $τ>0$, $I_α$ is the Riesz potential of order $α\in (\max\{0,N-2p\}, N)$, $p^*_α=\frac{p}{2}\left(\frac{N+α}{N-p}\right)$ is the critical exponent corresponding to the Hardy Littlewood Sobolev inequality, $(-Δ_p)^s$ is the non-local fractional p-Laplacian operator with $s\in (0,1)$, $μ>0$ is a parameter and $λ$ appears as a Lagrange multiplier. We show the existence of at least two distinct solutions in the presence of a mass subcritical perturbation, $μ|u|^{q-2}u$ with $p<q<p+\frac{sp^2}{N}$ under some conditions on $p,N$ and $s$.

4. 从多体量子Gibbs态导出最优质量区域内聚焦的 $Φ^6_1$ 测度

Derivation of the focusing $Φ^6_1$ measure in the optimal mass regime from many-body quantum Gibbs states

作者: Lin Lü, Phan Thành Nam, Rongchan Zhu

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

发布日期: 2026-05-25

arXiv ID: 2605.25755v1

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

摘要: We derive the focusing $Φ^6_1$ measure on the torus $\mathbb{T}$ as the high-temperature/mean-field limit of many-body quantum Gibbs states with an attractive three-body interaction. The main difficulty in the focusing setting is to relate the classical mass cutoff to the quantum particle-number cutoff. Our result reaches the optimal mass threshold for the classical field identified by Oh, Sosoe, and Tolomeo (2022), and thereby extends the earlier work of Rout and Sohinger (2025). At the critical threshold, the short-range interaction is allowed to shrink to a Dirac delta function on a logarithmic scale in the temperature parameter. Strictly below the threshold, the same convergence holds with a polynomial dependence on the temperature. Moreover, we establish a quantum-level phase transition at the same mass threshold. The proof develops the variational framework of Lewin, Nam, and Rougerie (2015) in the focusing setting and relies on two new ingredients: a non-factorized trial state construction and a delicate tail estimate for the interacting lower symbol. These allow us to control the localization and relative entropy errors caused by the particle-number cutoff, as well as the contribution of the focusing exponential weight.

5. 加权 1-拉普拉斯算子的首个狄利克雷特征值

First Dirichlet eigenvalue of the weighted 1-Laplacian operator

作者: R. Barbato, J. C. Sabina de Lis, S. Segura de León

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25642v1

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

摘要: In this paper, we study the eigenvalue problem

$$ \left\{\begin{array}{cl}-\hbox{div}\left(a(x)\frac{Du}{|Du|}\right)=Λ\, b(x)\frac{u}{|u|} & \text{in }Ω\\u=0 & \text{on }\partialΩ,\end{array}\right. $$

where $a(x)$ and $b(x)$ are suitable nonnegative functions. We prove that the first eigenvalue coincides with the weighted Cheeger constant. To see this identity, we analyze the behavior of the first Dirichlet eigenvalue of the weighted $p$-Laplacian operator as $p$ goes to $1$. In the case that the weight $a(x)$ is Lipschitz-continuous, we show that the limit of eigenvalues of the weighted $p$-Laplacian exists, and it is the weighted Cheeger constant. In addition, we check that the sequence of normalized $p$-eigenfunctions converges to the normalized eigenfunction of our limiting problem, which turns out to be bounded. For more general weights, we identify the first eigenvalue of the weighted 1-Laplacian operator with the weighted Cheeger constant and prove that the associated eigenfunction is bounded.

6. 具有非线性边界条件的麦克斯韦方程组的形状导数

Shape Derivatives for Maxwell’s Equations with Nonlinear Boundary Conditions

作者: Chao Deng, Yixian Gao

类别: math.AP

发布日期: 2026-05-25

arXiv ID: 2605.25579v1

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

摘要: This paper develops a trace-regular variational framework for time-harmonic Maxwell scattering problems involving pointwise nonlinear boundary and interface responses. We investigate three canonical classes of models: nonlinear impedance, nonlinear perfect electric conductor, and nonlinear transmission conditions. Since the standard Maxwell tangential trace belongs to a space of negative order, the nonlinearities are formulated in refined functional spaces where the tangential electric field admits an $L^2(Γ)$-trace. Under the assumption of a sufficiently small Lipschitz constant for the nonlinear response, we establish the well-posedness of the direct problems via fixed-point arguments leveraging the mapping properties of the associated linear Maxwell operators. Within this framework, we perform a rigorous sensitivity analysis of the electromagnetic fields with respect to perturbations of the scattering interface. By employing the covariant Piola transform, we prove the continuity and Fréchet differentiability of the pulled-back solutions with respect to domain variations. The material derivative is characterized as the unique solution to a corresponding $\mathbb{R}$-linearized Maxwell system, and the shape derivative is shown to satisfy explicit boundary or interface conditions for each of the three nonlinear models. We further demonstrate that the resulting sensitivity expressions possess the Hadamard structure, depending exclusively on the normal component of the boundary deformation. The resulting derivative characterizations provide a mathematical basis for subsequent adjoint-based sensitivity analysis, shape optimization, and gradient-driven inverse reconstruction.

7. 具尖点的Poincaré-Einstein 4维流形

Poincaré-Einstein 4-manifolds with cusps

作者: Mingyang Li, Hongyi Liu

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

发布日期: 2026-05-25

arXiv ID: 2605.25462v1

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

摘要: In this paper, we construct Poincaré-Einstein 4-manifolds with various kinds of cusps. In particular, we construct: (1) Infinite families of Einstein metrics on $(0,\infty)\times \mathscr{N}$, where $\mathscr{N}\to T^2$ is a principal $\mathbb{S}^1$-bundle over $T^2$, with one Poincaré-Einstein end and one end asymptotic to a real or complex hyperbolic cusp. (2) Infinite families of Einstein metrics on $(0,\infty)\times P$, where $P\to Σ_{\mathtt{g}}$ is a principal $\mathbb{S}^1$-bundle over a closed Riemann surface $Σ_{\mathtt{g}}$ of genus $\mathtt{g}\geq 2$, with one Poincaré-Einstein end and one end asymptotic to a bundle of two-dimensional hyperbolic cusps over hyperbolic $Σ_{\mathtt{g}}$. Universal covers of (1) and (2) provide new complete negative Einstein metrics on $\mathbb{R}^4$. These Einstein metrics also exhibit interesting degeneration phenomena. With this construction, we give a negative answer to a question of Anderson concerning cusp formation for Poincaré-Einstein 4-manifolds.

8. 一种具有宿主内病毒载量与抗体动力学的媒介传播疾病结构化模型

A structured model of vector-borne disease with within-host viral load and antibody dynamics

作者: Paulo Amorim, Maria Soledad Aronna, Debora de Oliveira Medeiros

类别: math.AP, math.NA

发布日期: 2026-05-25

arXiv ID: 2605.25370v1

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

摘要: We present an epidemiological model for vector-borne diseases that includes within-host viral load and antibody dynamics using structured transport equations. By incorporating the internal dynamics into the infected and recovered host compartments, the formulation introduces nonlinearities and nonlocalities. We establish analytical properties, including well-posedness and mass conservation, and characterize its characteristic curves. Furthermore, we derive a simplified Uniform Host Response (UHR) model featuring delay-type terms. For both the full and UHR frameworks, the basic reproduction number is determined and shown to serve as a threshold for the existence of an endemic equilibrium, and is related to the linear stability of the disease-free state. Finally, numerical experiments, parameterized specifically for Dengue fever, demonstrate how within-host mechanisms influence population-level epidemiological outcomes.

9. 关于高维Hersch-Weinberger不等式

On the Hersch-Weinberger inequality in higher dimensions

作者: T. V. Anoop, Vladimir Bobkov, Mrityunjoy Ghosh, Olga Pochinka

类别: math.AP, math.SP

发布日期: 2026-05-24

arXiv ID: 2605.25182v1

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

摘要: We investigate a reverse Faber-Krahn type inequality for the Robin Laplacian in a bounded smooth domain $Ω\subset \mathbb{R}^N$ whose boundary has two connected components. We prove that a concentric spherical shell maximizes the first eigenvalue over a class of such domains under perimeter and volume constraints, and under an additional convexity assumption when $N \geq 3$. This result generalizes to a wider class, and extends to higher dimensions, the inequality of Hersch [20], whose approach was substantially based on a construction of the so-called effectless cut by Weinberger [35], so that we call it the Hersch-Weinberger inequality. Our method is based on the analysis of the gradient flow of the first eigenfunction and several approximation procedures, without relying on the effectless cut itself. The effectless cut being a complicated object related to the attractor of the gradient flow, we describe its most fundamental topological properties. In particular, we show that it does not necessarily have to be a hypersurface.

10. 时空对数Sobolev不等式与欠阻尼Langevin动力学的弱强制超收缩性

Space-Time Log-Sobolev Inequality and Hypocoercive Hypercontractivity for Underdamped Langevin Dynamics

作者: Bowen Li, Jianfeng Lu

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

发布日期: 2026-05-24

arXiv ID: 2605.25083v1

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

摘要: We study hypercontractivity for the underdamped Langevin dynamics with a convex confining potential. Unlike in the overdamped case, the noise acts only on the velocity variable, so the usual argument based on the logarithmic Sobolev inequality (LSI) does not apply. Nevertheless, we prove that, when the spatial marginal satisfies an LSI with constant $ρ$ and the friction parameter is of order $\sqrtρ$, the underdamped Langevin semigroup improves integrability over the kinetic time scale $t \sim ρ^{-1/2}$. The proof rests on a novel space-time logarithmic Sobolev inequality, derived from a controlled version of the entropy decay estimate, which captures how dissipation in the velocity variable is transferred to the position variable. Combining this space-time LSI with a duality argument based on a forward/backward interpolation of the underdamped Langevin semigroup yields the desired hypocoercive hypercontractivity estimate. As a corollary, we obtain decay of the Rényi divergence at the sharp hypocoercive rate $\mathcal{O}(\sqrtρ)$.

11. 非凸优化的最优漂移优化器

Optimal drift optimizer for non-convex optimization

作者: Qin Li, Sixu Li, Eitan Tadmor, Emmanuel Trélat

类别: math.OC, math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24485v1

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

摘要: We study a finite-horizon stochastic control criterion for non-convex optimization in which Brownian exploration is balanced against a quadratic control cost. Rather than emphasizing the classical Hopf–Cole representation, we isolate the exact drift selected by the criterion and reorganize it in a form adapted to optimization. The key object is the conditional terminal law of the optimal process. We show that this law is a Gibbs measure for a proximally penalized energy, yielding three exact representations of the drift: potential, averaged-gradient, and barycentric. We then analyze two asymptotic regimes relevant for optimization. As terminal time is approached, the drift recovers a scaled gradient-descent field. In the low-temperature regime, assuming a unique global minimizer, the conditional terminal law concentrates on it even in the presence of nonglobal local minima, and the drift converges to an affine attraction field toward it. In the nondegenerate case we also derive Laplace asymptotics for the drift, the value function, and the covariance of the conditional terminal law. Finally, we record a simple gradient-free discretization suggested by the barycentric formula.

12. 关于边值问题 $-\varepsilon \Delta_p u = a|u|^{q-2}u - b|u|^{\gamma-2}u$ 的基态的渐近性

On the asymptotics of ground states for a boundary value problem for the equation $-\varepsilon Δ_p u = a|u|^{q-2}u - b|u|^{γ-2}u$

作者: Yavdat Sh. Il’yasov, Elvira I. Turianova

类别: math.AP, math-ph

发布日期: 2026-05-23

arXiv ID: 2605.24482v1

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

摘要: We study a singularly perturbed Dirichlet problem for the $p$-Laplacian with competing superlinear terms,

$$ -\varepsilon Δ_p u = a(x)|u|^{q-2}u - b(x)|u|^{γ-2}u, \qquad u|_{\partialΩ}=0, $$

where $1<p<q<γ<p^*$, $a\geq 0$, $b\geqσ_b>0$, and $\varepsilon>0$ is small. By means of the nonlinear Rayleigh quotient method, we introduce two critical parameter values, $\varepsilon^*$ and $\varepsilon_e^*$, related respectively to the Nehari manifold and to the zero energy level. We prove the nonexistence of nontrivial weak solutions for $\varepsilon>\varepsilon^*$, and the existence of at least two positive weak solutions for $0<\varepsilon<\varepsilon_e^*$; one of them is a ground state. The main result describes the asymptotic behaviour of ground states as $\varepsilon\to 0^+$. If, in addition, $a\geqσ_a>0$, then every family of positive ground states $u_\varepsilon$ converges in measure in $Ω$ to the explicit profile

$$ \bar u_0(x) = \left(\frac{a(x)}{b(x)}\right)^{1/(γ-q)}. $$

Moreover,

$$ u_\varepsilon\to\bar u_0 \quad\text{strongly in }L^r(Ω), \qquad 1\le r<γ, $$

and

$$ u_\varepsilon\rightharpoonup\bar u_0 \quad\text{weakly in }L^r(Ω), \qquad 1<r\leγ. $$

13. Caffarelli 收缩定理的无维数插值

A dimension-free interpolation of Caffarelli’s contraction theorem

作者: Bader Ammari, Alessio Figalli

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24443v1

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

摘要: We prove global Lipschitz estimates for Brenier maps between probability measures on $\mathbb{R}^n$ whose densities belong to the family

$$ ρ_{U,\,p}=Z_{U,\, p}^{-1}\exp(-Θ_p(U)), \qquad Θ_p(t)=p\log\Bigl(1+\frac{t}{p}\Bigr), \qquad p\in[n,+\infty], $$

with finite normalization constant $Z_{U,\, p}$, and with the convention $Θ_{\infty}(t)=t$. We allow different parameters for source and target, $d,D\in[n,+\infty]$, with $d\le D$. Our global estimate is uniform in $n,d,D$, and in the case $d=D<+\infty$, it improves the bounds of arXiv:2404.05456 by removing their exponential dependence on the dimension. We also prove localized estimates inside fixed balls $B_R$ whose constants are stable under the limits $d,D\to+\infty$ and they allow us to recover Caffarelli’s celebrated contraction theorem with sharp constants.

14. 广义 $p$-面积极小曲面的定量稳定性

Quantitative Stability of Generalized $p$-Area Minimizing Surfaces

作者: Amir Moradifam, Gerardo Orozco-Fernandez

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24317v1

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

摘要: We study the stability of $p$-area minimizing surfaces in the Heisenberg group under perturbations of the weight function and the drift vector field in generalized least gradient problems of the form

$$ \inf_{w\in BV_0(Ω)} \int_Ω\left(a(x)|Dw+F(x)|+H(x)w\right)\,dx. $$

Owing to the lack of strict convexity, establishing stability of minimizers is challenging. We derive quantitative stability estimates for minimizers. In particular, under suitable nondegeneracy and geometric assumptions, we obtain $L^1$ and $W^{1,1}$ stability estimates with respect to perturbations of the weight function $a$ and the drift vector field $F$. We further establish unified quantitative stability estimates under simultaneous perturbations of all principal parameters, namely $a$, $F$, and $H$. Numerical simulations illustrating the stability theory are also presented.

15. 最优传输映射的稳定性与2-Monge-Kantorovich距离的二阶变分

Stability of optimal transport maps and second variation of the 2-Monge-Kantorovich distance

作者: F. -U. Caja-Lopez, Matias G. Delgadino, Jun Kitagawa

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.24232v1

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

摘要: We establish several quantitative stability estimates for optimal transport maps between non-degenerate densities on uniformly convex domains for the quadratic cost. Under Hölder regularity assumptions, we prove Lipschitz $L^2$ (respectively $C^{1,α}$) stability estimates for optimal transport maps in terms of the 2-Monge-Kantorovich distance (respectively $L^{p}$ distances) between pairs of source and target densities. When the continuity assumption is removed, we obtain a Lipschitz $L^2$ stability estimate for the Brenier potentials in terms of the $L^2$ distance between the source and target densities. The proofs rely on a precise characterization of the linear response of the Brenier potential along smooth interpolations of the data, obtained by linearizing the Monge-Ampère equation in divergence form. As a further application of this approach, we derive an explicit formula for the second variation of the quadratic Monge-Kantorovich distance.

16. 具有轴向力和内部时滞的 Euler-Bernoulli 梁的边界稳定化

Boundary stabilization of an Euler-Bernoulli beam with axial force and internal delay

作者: Ben Bakary Junior Siriki, Adama Coulibaly

类别: math.AP

发布日期: 2026-05-23

arXiv ID: 2605.24315v1

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

摘要: This paper investigates the boundary stabilization of an Euler-Bernoulli beam under constant axial tension and subject to an internal time-delay. First, the well-posedness of the system is established using semigroup of linear operators theory. Then, by constructing an adequate Lyapunov functional, we establish the exponential stability of the system within a region of the $(α, ξ)$-plane, explicitly characterized by a system of inequalities involving the physical parameters. This work extends the results of Han et al. (Journal of Difference Equations, 2016) by incorporating the effect of axial tension and an internal delay term without sign constraints.

17. 一类时间切换的波动方程与扩散方程系统的新型反问题

New inverse problems for a time-switched system of wave and diffusion equations

作者: E. T. Karimov, N. A. Murolimova

类别: math.AP, math-ph

发布日期: 2026-05-22

arXiv ID: 2605.24189v1

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

摘要: We study two new classes of inverse problems for a time-switched system in which a fractional wave equation (with Caputo derivative of order $α\in (1,2)$) governs the dynamics on the interval $[0,a)$, and a fractional diffusion equation (with Caputo derivative of order $β\in (0,1)$ taken with respect to the switching point $t=a$) governs the dynamics on $(a,b]$. The two problems differ in which part of the transmitting condition at the interface $t=a$ is regarded as unknown. In both cases the overdetermination data consist of a single spatial measurement of the solution at a fixed time $ξ\in (a,b)$. Using the spectral expansion method with respect to the classical Sturm-Liouville eigensystem on $[0,1]$, we reduce each problem to a sequence of coupled scalar Cauchy problems involving the two-parameter Mittag-Leffler function. Explicit series representations for the solution $u(t,x)$ and the unknown interface functions $h(x)$ and $\bar{h}(x)$ are derived. Uniform convergence of the resulting infinite series and their relevant derivatives is established through four auxiliary lemmas, using the decay estimates for the Mittag-Leffler function, integration-by-parts arguments, the Cauchy–Schwarz inequality, and the Weierstrass $M$-test. A uniqueness and existence theorem is stated for Problem1 under explicit Sobolev-type regularity conditions on the data, with an analogous result for Problem2.

18. 极化诱导的光束弯曲：数学模型、离散化与算法

Polarization-Induced Beam Bending: Mathematical Model, Discretization, and Algorithm

作者: Harbir Antil, Rainald Löhner, Sarswati Shah

类别: physics.optics, math.AP, math.NA

发布日期: 2026-05-22

arXiv ID: 2605.24151v1

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

摘要: We study a reduced hydrodynamic formulation of paraxial vector beam propagation in which the beam intensity, optical phase, and spatially-dependent polarization are coupled through a nonlinear dispersive system. While prior analytical work derived a solution for the beam path valid for short propagation distances, a fully resolved numerical treatment of the model over long ranges has not previously been available. Here we present a conservative numerical scheme for the coupled system, combining a finite-volume discretization of the intensity equation with monotone Hamilton–Jacobi (H-J) solvers for the phase dynamics and upwind transport of polarization. The method preserves the nonnegativity of the intensity and remains stable under long-distance propagation. We perform large-scale simulations over propagation distances of tens of meters, while resolving millimeter-scale transverse structure. The numerical results reproduce the analytically predicted and experimentally observed quadratic beam bending at short distances and reveal systematic deviations beyond the asymptotic regime. These deviations arise from nonlinear phase accumulation and dispersive effects captured by the full model but are neglected in the short-distance approximation.