生成时间: 2026年05月07日 06:24:54

重点关注论文数量: 8

了解领域论文数量: 9

不相关论文数量: 12

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

1. 二维 Kaup-Broer-Kuperschmidt Boussinesq 系统的数值研究

Numerical study of the 2D Kaup-Broer-Kuperschmidt Boussinesq system

作者: Théo Gaudry, Christian Klein, Jean-Claude Saut, Nikola Stoilov

类别: math.AP, math.NA

发布日期: 2026-05-06

arXiv ID: 2605.05183v1

摘要翻译: 本文考虑二维情形下适定化的 Kaup-Broer-Kuperschmidt 系统。我们数值构造了孤子型解，并证明其既在色散作用下不稳定，也在奇点形成意义下不稳定。进一步，我们研究了线孤子及其稳定性，以及一般具有局部支集的初值数据。在上述任一情形中，我们都未能发现稳定的结构。

Comment: 无

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

1. 研究对象和背景

研究对象为二维Kaup-Broer-Kuperschmidt (KBK) Boussinesq系统，这是Boussinesq ‘abcd’系统在二维的特例。系统分为线性不适定的’bad’版本和适定的’good’版本，本文聚焦后者：

其中$\eta$为自由表面高度，$\mathbf{v}$为速度场。该系统具有双哈密顿结构，其哈密顿量为：

尽管能量守恒，但不足以保证Cauchy问题的全局解存在。目前已知最好结果是长时间存在性（见[33]），即对小参数$\epsilon$系统：

存在时间尺度为$O(1/\epsilon)$的解。一维KBK系统完全可积，但其解的定性性质远未完善，如紧致初始数据分解为孤子和辐射的长时间行为（孤子分解猜想）尚未解决。二维推广后系统不可积，主要开放问题包括全局解的长时间行为和有限时间爆破可能性。作为[27]一般研究的结果，一维和二维中存在线性色散估计（如Strichartz和局部Kato平滑估计）。本文通过数值方法研究二维KBK系统的解性质，特别是寻找稳定结构，包括：1. 二维孤子及其稳定性：通过引入速度势$V$（$\mathbf{v} = \nabla V$），系统重写为：

2. 线孤子稳定性：一维孤波解（$C \in \mathbb{R}$，$|C| < 1$）：

在二维中简单横向扩展后的稳定性；3. 静态解研究：时间无关解满足：

径向对称情形下简化为：

4. 爆破机制分析：通过动态重标度$X = (x - x_s)/L$，$Y = (y - y_s)/L$，$\eta = U/L^2$，$d\tau/dt = 1/L^2$，将系统转化为：

其中$a = L_\tau/L$。$L^2$超临界特性暗示可能的自相似爆破$L \propto (t^* - t)^{1/2}$，$\eta$的$L^\infty$范数增长率为$(t^* - t)^{-1}$。本文填补了二维KBK系统数值研究的空白，此前仅有少量一维理论结果（如Angulo[2]利用高阶守恒量$I_3$证明$H^s(\mathbb{R}) \times H^{s-1}(\mathbb{R})$，$s \geq 2$上的全局适定性和轨道稳定性），而二维理论几乎空白。

2. 主要定理或主要结果

本文主要通过数值实验得出以下结论，未提供严格证明但具有高精度验证（相对能量误差$\Delta E < 10^{-3}$）：1. 静态解构造与扰动稳定性：- 成功构造径向对称静态解，满足方程$(4sV_{ss} + 4V_s - V)_s + 2sV_s^3 = 0$（$s = r^2$），使用$N = 150$阶Chebyshev离散化和松弛参数$\mu = 0.1$，残差低于$10^{-9}$。$\eta$的$L^2$范数$\|\eta\|_2^2 \sim 96.5967$。- 扰动$\eta(x,y,0) = \lambda \eta_{\text{static}}$，$V(x,y,0) = V_{\text{static}}$：* $\lambda = 0.99$（低于临界能量）：使用$L_x = L_y = 5$，$N_x = N_y = 2^9$，$N_t = 10^4$步模拟$t \leq 10$，解完全色散，$\|\eta\|_\infty$和$\|\eta\|_2^2$单调递减至零；* $\lambda = 1.01$（高于临界能量）：使用$L_x = L_y = 5$，$N_x = N_y = 2^{11}$，$N_t = 2 \times 10^4$步模拟$t \leq 5.2$，在$t \approx 5.1740$发生爆破，$\|\eta\|_\infty \propto (t^* - t)^{-1}$，$\|\eta\|_2^2$同步增长，符合$L^2$超临界系统自相似爆破特征$L \propto (t^* - t)^{1/2}$。2. 线孤子横向稳定性：- 以一维孤波（$C = 0.8$）为初始数据，在$L_x = 20$，$L_y = 5$，$N_x = 2^{11}$，$N_y = 2^7$网格上模拟$t \leq 1$（$N_t = 2000$步）。- 数值解与精确解差异$\sim 10^{-11}$，能量守恒精度$\sim 10^{-12}$，但线孤子呈现强不稳定性，横向扰动导致快速失稳。3. 高斯初始数据演化：- 对一般紧致初始数据（如高斯型）的数值实验均未发现稳定结构，解要么完全色散，要么在有限时间内爆破。4. 爆破机制验证：- 对$\lambda = 1.01$情形的动态重标度分析确认自相似爆破，$L(\tau) \propto \exp(-a^\infty \tau)$对应$L \propto (t^* - t)^{1/2}$，$a^\infty \approx 0.5$（通过拟合$\log L$与$\tau$线性关系确定）。- $\eta$的$L^\infty$范数增长符合$(t^* - t)^{-1}$标度，与理论预测一致，数值爆破时间$t^* = 5.1740 \pm 10^{-5}$通过外推法精确确定。

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

本文核心创新在于开发了高精度数值框架处理二维KBK系统的三大挑战：静态解构造、刚性时间演化和爆破检测。关键技术如下：1. 静态解的数值构造：- 极坐标奇异性处理：将径向方程$V_r - (V_{rr} + r^{-1}V_r)_r - \frac{1}{2}V_r^3 = 0$通过变量替换$s = r^2$转化为$(4sV_{ss} + 4V_s - V)_s + 2sV_s^3 = 0$，显著降低$r=0$处的奇异性。在$s \in [0, s_0]$（$s_0 = 100 \gg 1$）上应用$s = s_0(1+l)/2$变换，$l \in [-1,1]$采用$N=150$阶Chebyshev离散化$l_j = \cos(j\pi/N)$。导数通过Chebyshev微分矩阵近似，边界条件用$\tau$-方法实现。- 非平凡解收敛技术：针对牛顿迭代易收敛到平凡解$V=0$的问题，引入松弛参数$\mu=0.1$，迭代步更新为$V^{(n+1)} \leftarrow \mu V^{(n+1)} + (1-\mu)V^{(n)}$，初始猜测$V=4\exp(-s)$，使残差降至$10^{-9}$以下。- 二维傅里叶网格插值：将极坐标解插值到$\mathbb{T}^{2\pi L_x} \times \mathbb{T}^{2\pi L_y}$环面，通过牛顿-Krylov方法求解非线性系统$(-\frac{1}{2}|\nabla V|^2 \nabla V) + \nabla (V - \Delta V) = 0$，避免了对称性限制。2. 刚性时间演化的高效求解：- 分数阶NLS形式转化：将系统写为$(\hat{u}_\pm)_t = L_\pm \hat{u}_\pm + N_\pm$，其中$L_\pm = \mp i|k|\sqrt{1 + |k|^2}$，$N_\pm$为非线性项。此形式色散关系$|k|^2$（大$|k|$）使系统具有强刚性，需特殊处理。- 指数时间差分(ETD)方案：采用Cox-Matthews 4阶ETD格式，将线性项精确积分，非线性项通过多项式插值处理。对于$U_{n+1} = e^{L\Delta t}U_n + \int_0^{\Delta t} e^{L(\Delta t-\tau)}N(U(\tau))d\tau$，使用4阶精度近似积分项。该方法避免了显式格式的CFL限制和隐式格式的高计算成本。- 误差控制策略：通过精确守恒能量$E(t)$监控数值误差，定义相对能量误差$\Delta E = |E(t)/E(0)-1|$。当$\Delta E > 10^{-3}$时自动减小时间步长，确保数值解的可靠性。对于爆破情形，采用自适应时间步长$\Delta t \propto L^2$保持动态重标度下的精度。3. 爆破机制的定量分析：- 动态重标度实现：在检测到$\|\eta\|_\infty$快速增长后，自动切换到重标度变量$(X,Y,U,\tau)$，通过监测$L(\tau) = \|\eta\|_\infty^{-1/2}$验证自相似性。使用最小二乘法拟合$\log L$与$\tau$的线性关系，确定$a^\infty = -d(\log L)/d\tau$。- 爆破时间外推：基于$\|\eta\|_\infty(t) \sim (t^*-t)^{-1}$标度，对$1/\|\eta\|_\infty$进行多项式外推，结合Richardson外推法将爆破时间精度提升至$\pm 10^{-5}$。- 谱分析验证：通过计算$\hat{\eta}(k,t)$在爆破前的谱分布，验证$|\hat{\eta}| \sim k^{-\alpha}$的幂律行为（$\alpha \approx 3$），符合$L^2$超临界系统的理论预测。

4. 与之前工作的比较

本文与既有工作的比较揭示了二维KBK系统的独特性质和数值研究的推进：1. 与一维KBK理论的对比：- Angulo[2]利用高阶守恒量$I_3$证明了1D KBK在$H^s(\mathbb{R})\times H^{s-1}(\mathbb{R})$（$s\geq 2$）上的全局适定性和孤波轨道稳定性，但仅限于初始数据靠近孤波平移的情形。本文发现2D中即使存在静态解，其稳定性窗口极窄（$\lambda=0.99$色散，$\lambda=1.01$爆破），且无轨道稳定性证据。- [38]证明了1D孤波在$C\in(-1,-1/2)$时的渐近稳定性，而本文显示2D线孤子（$C=0.8$）横向扰动下迅速失稳，与KP-I方程行为类似，但不同于KP-II方程的横向稳定性。- 1D中孤子分解猜想（紧致数据分解为孤子+辐射）尚未解决，而2D数值实验表明解要么完全色散（低能量），要么爆破（高能量），未观察到孤子-辐射分解现象。2. 与Boussinesq ‘abcd’系统的比较：- [33]对小参数$\epsilon$系统证明了$O(1/\epsilon)$时间尺度的存在性，但未涉及2D全局行为。本文揭示即使’good’版本在2D也可能有限时间爆破，突破了$O(1/\epsilon)$框架的局限性。- 其他Boussinesq系统（如’bad’版本）需非空化条件$\eta+1>0$保证正则性，而KBK系统不受此限（见PDF第2页）。但本文发现KBK在2D的稳定性更差，表明适定性不等价于结构稳定性。- 与[27]的线性色散估计对比：Strichartz估计保证线性解的衰减，但本文显示非线性效应（$\lambda=1.01$）可克服色散导致爆破，揭示了非线性对线性理论的突破。3. 与NLS方程的类比与差异：- KBK的分数阶NLS形式（$L_\pm = \mp i|k|\sqrt{1+|k|^2}$）色散关系$\sim |k|^2$（大$|k|$），与$L^2$超临界NLS（如3D cubic NLS）类似。本文爆破标度$L\propto (t^*-t)^{1/2}$、$\|\eta\|_\infty\propto (t^*-t)^{-1}$与超临界NLS一致，但机制不同：NLS爆破源于非线性自聚焦，而KBK爆破涉及速度场$\mathbf{v}$与高度$\eta$的耦合。- 不同于NLS的$H^1$临界性，KBK的$L^2$超临界性由$\eta$的$L^2$范数非不变性决定（见PDF第5页），这是Boussinesq系统特有的几何结构。4. 与数值方法的演进对比：- [20]对1D KBK使用谱方法，但未处理2D问题。本文将ETD方案扩展到2D，通过分数阶NLS形式避免了传统Boussinesq系统数值方法的刚性问题，时间步长比显式格式提高2-3个数量级。- 与[12]的静态解求解相比，本文引入$s=r^2$变换和松弛技术，使非平凡解的构造残差从$10^{-6}$提升至$10^{-9}$，并首次实现2D环面上的非对称静态解计算。5. 局限性与代价：- 计算资源限制：高精度爆破模拟需$N_x=N_y=2^{11}$网格，单次模拟耗时>24小时，限制了参数扫描范围。可能遗漏狭窄稳定区域（如$|\lambda-1|<10^{-3}$）。- 数值方法固有局限：周期边界条件可能影响色散行为，但通过$L_x=L_y=5$（远大于解支撑）减轻了边界效应。自适应时间步长在爆破点附近可能引入人为平滑。- 与理论的差距：数值结果强烈暗示2D KBK无稳定孤子结构，但缺乏严格证明。特别地，无法确定是否存在指数脆弱的稳定解（如[38]中1D的渐近稳定情形）。6. 对领域的影响：- 本文否定性结果（无稳定结构）挑战了’可积系统推广应保留孤子稳定性’的直觉，表明二维化可能彻底改变系统动力学。- 爆破标度的精确验证为$L^2$超临界Boussinesq系统的分类提供了新基准，后续理论工作需解释为何KBK比类似系统（如KP-I）更不稳定。- 开发的数值框架可推广至其他二维非线性波系统，特别是具有分数阶色散关系的模型。

2. 关于三维磁流体动力学方程在速度-涡量形式下的部分Voigt正则化

On a Partial Voigt Regularization of the 3D Magnetohydrodynamic Equations in Velocity-Vorticity Form

作者: Adam Larios, Yuan Pei

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.05139v1

摘要翻译: 不可压Navier-Stokes方程的速度-涡量（VV）形式近年来广受关注，尤其在数值研究中，因其具有结构性优势。近期，我们与L. Rebholz合作，在该形式的动量方程中引入了Voigt正则化，建立了该正则化系统在三维情形下的整体适定性，并给出了收敛性结果及一个爆破准则。本文中，我们将这些思想推广至三维磁流体动力学（MHD）方程。尽管表面上看磁方程似乎需要一种“VV型”分解，但我们证明此类修改实非必要：仅对动量方程施加Voigt正则化即可保证整体适定性，同时保持涡量方程与磁方程原有的结构。我们还证明了该正则化系统在原系统可能爆破的时间之前收敛于原系统，并为原三维MHD系统的解建立了一个爆破准则。

Comment: 无

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

1. 研究对象和背景

本研究聚焦于三维周期域 $\mathbb{T}^3 = [0,1]^3$ 上磁流体动力学（MHD）方程在速度-涡度（Velocity-Vorticity, VV）形式下的部分Voigt正则化问题。原始3D MHD方程系统（1.1）为：

其中 $U$ 是速度场，$B$ 是磁场，$P$ 是压力，$\nu > 0$ 是流体粘度，$\eta > 0$ 是磁阻率。通过向量恒等式 $v \cdot \nabla v = (\nabla \times v) \times v + \frac{1}{2}\nabla|v|^2$ 和记号 $J := \nabla \times B$、$\Omega := \nabla \times U$，可将动量方程改写为：

其中 $\Pi := p + \frac{1}{2}|U|^2 - \frac{1}{2}|B|^2$。取动量方程的旋度得到涡度方程：

本文研究的正则化系统（1.5）为：

其中 $\alpha > 0$ 是正则化参数，$I$ 是恒等算子，$j = \nabla \times b$ 表示电流密度。与标准Voigt正则化不同，本系统仅在速度方程中引入Voigt项 $-\alpha^2\Delta \frac{\partial u}{\partial t}$，而保持涡度方程和磁场方程的原始结构。

研究该系统的动机源于MHD方程在数值模拟中的困难：非线性相互作用在大小尺度间形成反馈，需要极高分辨率的直接数值模拟，计算成本高昂。更重要的是，目前尚不清楚MHD方程的解是否会产生任意小的长度尺度导致激波形成。Voigt正则化因其保留原始系统稳态解、在物理边界条件下全局适定等优势而受到关注。特别地，本工作解决了关键问题：在MHD系统中是否需要对磁场方程也进行VV-type分裂？作者证明仅需在动量方程中添加Voigt正则化即可保证全局适定性，无需修改磁场方程结构，这与直观预期不同。本文填补了VV形式下MHD方程正则化理论的空白，扩展了作者先前在Navier-Stokes方程上的工作（Larios与Rebholz的VVV系统），并建立了原始MHD系统的爆破准则。

2. 主要定理或主要结果

定理2.4（弱解）：对任意初始数据 $u_0 \in V$ 和 $w_0, b_0 \in H$，系统（1.5）在定义2.1意义下存在唯一的全局弱解，满足 $\nabla \cdot w = 0$。且满足能量恒等式：

定理2.5（强解）：对初始数据 $u_0, w_0, b_0 \in V$，系统（1.5）在定义2.2意义下存在唯一的全局强解。若进一步假设 $u_0, w_0, b_0 \in H^s \cap V$（$s \geq 2$ 为整数），则解满足：

定理2.6（$w$ 到 $\omega = \nabla \times u$ 的收敛性）：对初始数据 $u_0, b_0 \in H^3 \cap V$ 且 $w_0 = \nabla \times u_0$，存在常数 $K_0, C > 0$ 使得：

其中时间区间 $[0, T]$ 满足 $u(t) \in L^\infty(0, T; H^4) \cap L^2(0, T; H^4)$ 和 $w(t) \in L^\infty(0, T; H^2) \cap L^2(0, T; H^3)$。当 $\alpha \to 0^+$ 时，$\|\omega(t) - w(t)\|_{L^2} \to 0^+$ 对几乎所有 $0 < t < T$ 成立。

定理2.7（正则化系统到原始系统的收敛性）：对初始条件 $U_0 = u_0 \in H^3 \cap V$，$B_0 = b_0 \in H^3 \cap V$ 且 $w_0 = \nabla \times u_0$，记 $\Omega = \nabla \times U$，$J = \nabla \times B$，$j = \nabla \times b$，则在时间区间 $[0, T]$ 上（其中 $u, U, b, B \in L^\infty(0, T; H^3) \cap L^2(0, T; H^4)$ 且 $w \in L^\infty(0, T; H^2) \cap L^2(0, T; H^3)$），存在仅依赖于 $U, B$ 的 $H^3$ 范数而与 $\alpha$ 无关的常数 $C$，使得：

当 $\alpha \to 0^+$ 时，$\|u(t) - U(t)\|_{L^2} \to 0^+$，$\|w(t) - \Omega(t)\|_{L^2} \to 0^+$，$\|b(t) - B(t)\|_{L^2} \to 0^+$，且 $\|j(t) - J(t)\|_{L^2} \to 0^+$ 对几乎所有 $0 < t < T$ 成立。

定理2.8（爆破准则I）：若存在 $0 < T < \infty$ 和 $\epsilon > 0$ 使得

则原始MHD系统（1.1）的解在 $[0,T]$ 上发展出奇点。

推论2.9（爆破准则II）：若存在 $0 < T < \infty$ 和 $\epsilon > 0$ 使得

则原始MHD系统（1.1）的解在 $[0,T]$ 上发展出奇点。

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

证明的核心技术障碍在于：如何在仅对速度方程添加Voigt正则化的情况下，建立全局适定性并控制涡度与磁场方程的非线性项。作者的创新点在于发现无需对磁场方程进行额外正则化，即可保持系统结构并实现全局适定性。关键估计如下：

1. 能量估计框架：对系统（1.5a）和（1.5c）分别与 $u$ 和 $b$ 做内积，利用散度自由条件得到能量恒等式（定理2.4）：

   $$ \frac{d}{dt} \left( \|u\|_{L^2}^2 + \alpha^2 \|\nabla u\|_{L^2}^2 + \|b\|_{L^2}^2 \right) + 2\nu \|\nabla u\|_{L^2}^2 + 2\eta \|\nabla b\|_{L^2}^2 = 0, $$

   此恒等式直接导出 $u \in L^\infty(0,T;H^1) \cap L^2(0,T;H^2)$ 和 $b \in L^\infty(0,T;L^2) \cap L^2(0,T;H^1)$ 的先验估计。

2. 涡度与电流密度的高阶估计（定理2.5）：对磁场方程（1.5c）取旋度得 $j = \nabla \times b$ 的演化方程：

   $$ \frac{\partial j}{\partial t} - \eta \Delta j = \nabla \times \left( (b \cdot \nabla)u - (u \cdot \nabla)b \right), $$

   与 $j$ 做内积后，利用Agmon不等式 $\|u\|_{L^\infty} \leq C \|u\|_{H^1}^{1/2} \|u\|_{H^2}^{1/2}$ 和Young不等式估计非线性项：

   $$ \left| \int_{\mathbb{T}^3} (b \cdot \nabla)u \cdot (\nabla \times j) \right| \leq C \|j\|_{L^2}^{1/2} \|\nabla j\|_{L^2}^{3/2} \|\nabla u\|_{L^2} \leq C_{\alpha,\eta} \|j\|_{L^2}^2 + \frac{3\eta}{4} \|\nabla j\|_{L^2}^2, $$
   $$ \left| \int_{\mathbb{T}^3} (u \cdot \nabla)b \cdot (\nabla \times j) \right| \leq C \|\nabla u\|_{L^2} \|j\|_{L^2}^{1/2} \|\nabla j\|_{L^2}^{3/2} \leq C_{\alpha,\eta} \|j\|_{L^2}^2 + \frac{3\eta}{4} \|\nabla j\|_{L^2}^2, $$

   从而得到：

   $$ \frac{1}{2} \frac{d}{dt} \|j\|_{L^2}^2 + \frac{\eta}{4} \|\nabla j\|_{L^2}^2 \leq C_{\alpha,\eta} \|j\|_{L^2}^2, $$

   通过Gronwall不等式导出 $b \in L^\infty(0,T;H^1) \cap L^2(0,T;H^2)$。

3. 收敛性证明的关键技术（定理2.6-2.7）：引入误差变量 $\delta u = u - U$，$\delta w = w - \Omega$，$\delta b = b - B$，建立误差方程并估计。对定理2.6，考虑 $\omega = \nabla \times u$ 与 $w$ 的差，利用以下恒等式：

   $$ \frac{d}{dt} \left( \|\omega - w\|_{L^2}^2 + \alpha^2 \|\nabla(\omega - w)\|_{L^2}^2 \right) + \|\nabla(\omega - w)\|_{L^2}^2 \leq C \left( \|\omega - w\|_{L^2}^2 + \alpha^2 \|\nabla(\omega - w)\|_{L^2}^2 \right), $$

   结合Gronwall不等式得到指数收敛率。对定理2.7，需同时处理速度、涡度、磁场和电流密度的误差，关键步骤是利用原始系统解的正则性（$U,B \in H^3$）控制交叉项，最终导出 $O(\alpha^2)$ 的收敛速率。

4. 爆破准则的构造：基于 $\alpha$-能量平衡

   $$ \alpha^2 \|\nabla u(t)\|_{L^2}^2 + \|u(t)\|_{L^2}^2 \leq \alpha^2 \|\nabla u_0\|_{L^2}^2 + \|u_0\|_{L^2}^2, $$

   若 $\limsup_{\alpha \to 0^+} \alpha \sup_{t \in [0,T]} \|\nabla u(t)\|_{L^2} \geq \epsilon$，则原始解 $U$ 的 $H^1$ 范数必在 $[0,T]$ 内爆破。这推广了[35]中SQG方程的爆破准则思想到MHD系统。

4. 与之前工作的比较

1. 与VVV系统（Navier-Stokes）的比较： 作者先前与Rebholz在[41]中研究了Navier-Stokes方程的速度-涡度-Voigt（VVV）系统，其中仅在速度方程添加Voigt正则化。本文将该框架扩展到MHD系统，但发现关键区别：对于MHD，无需对磁场方程进行类似涡度方程的VV-type分裂（即无需为电流密度 $j$ 添加额外演化方程）。在VVV-MHD系统（1.5）中，仅需三个演化方程（速度、涡度、磁场），而直观预期可能需要四个（速度、涡度、磁场、电流密度）。这表明MHD系统的结构对正则化的需求弱于Navier-Stokes系统，因为磁场方程的非线性项 $(u \cdot \nabla)b - (b \cdot \nabla)u$ 本身已提供足够正则性。

2. 与Voigt-MHD（V-MHD）系统的比较： [44]提出的V-MHD系统仅在速度方程添加Voigt项，但将速度方程取旋度后，涡度方程会继承Voigt项 $-\alpha^2 \Delta \omega_t$。而本文的VVV-MHD系统（1.5b）中涡度方程无Voigt项，这意味着VVV-MHD的正则化程度低于V-MHD。具体而言，V-MHD在 $\nu = 0, \eta > 0$ 时仍全局适定（[44]），而VVV-MHD要求 $\nu > 0, \eta > 0$（见Remark 1.1）。但VVV-MHD的优势在于保持了原始涡度方程的结构，更适合数值模拟中涡度动力学的精确捕捉。

3. 与经典Voigt正则化的比较： 传统Voigt正则化（如[43,58]）通常在所有演化方程中添加正则化项，而本文证明对MHD系统仅需在速度方程正则化即可。代价是：收敛速率从通常的 $O(\alpha)$ 提升为 $O(\alpha^2)$（定理2.7），但需假设原始解在 $[0,T]$ 上属于 $H^3$（而非一般Voigt文献中的 $H^2$）。此外，爆破准则中的条件更严格（要求 $\alpha \|\nabla u\|_{L^2}$ 有下界而非仅 $\|\nabla u\|_{L^2}$ 爆破）。

4. 与无粘情形的对比： [44]研究了V-MHD在 $\nu = 0, \eta > 0$ 的情况，但本文指出若在VVV框架下设 $\nu = 0$，则涡度方程将缺乏物理粘性（Remark 1.1），这与MHD系统的物理一致性矛盾。因此本文严格要求 $\nu > 0, \eta > 0$，确保速度-涡度方程的粘性匹配。这反映了VVV形式对物理模型一致性的敏感性：若涡度方程有粘性，速度方程必须有相应粘性。

5. 收敛性与爆破准则的推进： 本文的收敛结果（定理2.6-2.7）将[41]中VVV系统的收敛理论推广到MHD，但技术难点在于处理磁场与速度的双向耦合。特别是定理2.7同时控制了 $u-U$、$w-\Omega$、$b-B$、$j-J$ 四重误差，而Navier-Stokes情形仅需处理两重误差。爆破准则（定理2.8-2.9）则首次将[35]中SQG方程的 $\alpha$-能量方法应用于MHD系统，比经典Beale-Kato-Majda准则更易数值验证，因为只需监控正则化解的 $\alpha \|\nabla u\|_{L^2}$ 而非原始解的 $\|\nabla U\|_{L^\infty}$。

6. 局限与未来方向： 本文未解决 $\eta = 0$ 的情形（完全导电极限），而[43]证明在标准V-MHD中 $\eta = 0$ 时仍全局适定。作者在1.1节指出，若在磁场方程添加Voigt项 $-\alpha_M^2 \Delta b_t$，可能实现 $\eta = 0$ 时的全局适定性，这将是重要拓展。此外，与[61]中Boussinesq系统的对比表明，热耦合的VVV-MHD模型值得研究，但本文为简化未涉及。

3. 三维不可压活性液晶的全局适定性与衰减率

Global well-posedness and decay rates for the three dimensional incompressible active liquid crystals

作者: Fan Yang, Xiongfeng Yang

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04625v1

摘要翻译: 本文研究在恒定活性条件下三维不可压活性液晶的全局适定性与大时间行为，其模型由一个受迫不可压Navier-Stokes方程组（描述速度场）与一个关于 $Q$-张量序参量的抛物型方程组耦合而成。通过运用精细的交换子估计，我们证明了当初始数据 $(Q_0,u_0)\in H^{s+1}\times H^s$ ($s\geq 2$) 足够小、且活性 $c>c_\star$ 时，全局强解的存在性与唯一性；该结果改进了文献 \cite{active-limit} 中的先前结论。此外，若初始数据进一步属于 $L^1$ 且 $s\geq 4$，则对 $k\leq s-1$，我们得到关于 $\|\partial^kQ(t)\|_{L^2}$ 的混合衰减估计：该估计同时包含一个额外的指数衰减因子，其衰减速率为与 $(c-c_\star)\Gamma$ 成正比；以及一个最优代数衰减率，该速率与热核衰减速率一致。该结果表明，在高活性区域，活性向列相将趋于各向同性，且收敛具有依赖于活性的指数速率；并且该估计在无穷旋转粘度极限 $\Gamma\rightarrow 0$ 下保持稳定。与此同时，在附加初始假设下，我们亦导出了对 $\|\partial^ku(t)\|_{L^2}$ 的尖锐衰减估计，其中 $k\leq s-2$。证明结合了Green函数方法与时间加权能量方法。据我们所知，这些结果是首个在Beris-Edwards框架下针对活性/被动向列相液晶所建立的结果；而取向场增强的衰减效应本质上源于自由能。此外，在被动情形下，我们的结果蕴含热致向列相液晶在高温下的相变。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩活性液晶的Beris-Edwards模型，该模型由耦合的Navier-Stokes方程与Q-张量动力学方程构成：

其中$Q$为无迹对称张量序参量，描述向列相取向；$u$为速度场；$c > 0$为恒定活性浓度，$c_\star$为临界活性阈值；$a = \frac{1}{2}(c - c_\star)$，$\mathrm{tr}(Q^2)$项体现非线性自由能结构。当$c = c_\star$时，系统退化为被动液晶模型。已有工作：Paicu与Zarnescu[53]证明了二维被动模型的全局弱解存在性；Xiao[62]和Hieber等[32]在三维被动情形下建立了小初值全局强解；Chen等[15]处理了简化活性系统，Jiang等[37]首次研究全活性系统在环面上的强适定性，但要求高正则性且未分析大时间行为。本文旨在填补以下空白：(1) 在$\mathbb{R}^3$中将初始正则性降至$s \geq 2$；(2) 在$\mathbb{R}^3$中，揭示高活性下活性向列相向各向同性相的指数收敛机制；(3) 严格论证自由能项$aQ - cQ \mathrm{tr}(Q^2)$对衰减行为的增强作用，这是被动模型所不具备的特征。

2. 主要定理或主要结果

在$s \geq 2$且初始数据$(Q_0, u_0) \in H^{s+1} \times H^s$充分小时，系统在$\mathbb{R}^3$中存在唯一全局强解：

若进一步假设初值$(Q_0, u_0) \in L^1$且$s \geq 4$，则对$0 \leq k \leq s-1$，有混合衰减估计：

其中$a = \frac{1}{2}(c - c_\star)$，且该估计在$\Gamma \to 0$极限下稳定。对速度场，若初值满足额外条件，则对$0 \leq k \leq s-2$，有最优代数衰减：

$

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

证明的核心障碍在于三维非线性耦合系统中非线性项的高阶导数损失与活性项引发的指数衰减机制的数学刻画。创新点包括：(1) 低正则性全局适定性：通过精炼交换子估计，证明了

其中$\Lambda^s = (-\Delta)^{s/2}$，从而在$s \geq 2$下绕过[37]对$\|Q_0\|_{L^\infty}$小性的依赖；利用Sobolev范数与Besov范数等价性（引理A.6），将环面估计推广至$\mathbb{R}^3$，克服了$|Q|$项的非光滑性障碍。(2) 线性化系统的格林函数分析：引入辅助系统$\sigma := \mathrm{div} Q$，其线性化形式为

通过傅里叶变换显式求解格林函数$\hat{G}_{\sigma,u}(t, \xi)$，并利用傅里叶分裂法证明

从而导出线性解的混合衰减$\|\partial^k_x Q^L(t)\|_{L^2} \lesssim e^{-a\Gamma t}(1 + t)^{-3/4 - k/2}$，首次揭示自由能项$a\Gamma\sigma$产生的指数衰减因子。(3) 时间加权能量方法：为处理非线性项$u \cdot \nabla Q$与$Q\Delta Q$的导数损失，将积分分解为$[0, t/2]$与$[t/2, t]$，前者利用线性衰减，后者结合时间权重$(1+t)^\gamma$进行能量估计；对$\|B^{s-1}Q\|_{L^2}$采用卷积估计

以获得非线性解的混合衰减；尽管最高阶$\|B^s Q(t)\|_{L^2}$仅获$(1+t)^{-s/2}$衰减，但足以支撑低阶导数的最优估计，最终实现自由能驱动的指数增强衰减。

4. 与之前工作的比较

本文在全局适定性与衰减行为两方面显著推进既有工作：(1) 全局适定性：与Jiang等[37]相比，将三维初始正则性从隐含$s \geq 3$降至$s \geq 2$，并去除对$\|Q_0\|_{L^\infty}$小性的要求；与Xiao[62]的被动模型相比，本文处理活性项$\kappa\nabla\cdot Q$，且使用能量方法而非最大正则性，适用性更广；与Paicu-Zarnescu[53]的二维结果相比，本文克服了三维Navier-Stokes的困难并推广至全空间$\mathbb{R}^3$。(2) 衰减估计：与[19]、[47]的被动模型结果相比，本文首次发现活性模型中$\|\partial^k Q(t)\|_{L^2}$包含指数因子$e^{-(c-c_\star)\Gamma t/4}$，实现比纯代数衰减更优的混合衰减；与Schonbek-Shibata[47]的$L^p$衰减相比，本文定量揭示了衰减速率与$c, \Gamma$的依赖关系，且在$\Gamma \to 0$极限下稳定保留代数部分；定理1.4首次证明速度场衰减上界为最优，填补[19,47]中下界缺失的空白。(3) 物理意义：在被动情形下（$c = c_\star$），$a > 0$对应高温相变，本文数学结果解释了热致液晶从向列相向各向同性相的指数衰减过程，与实验[9,31,39,61]一致；同时，从数学上阐明了活性应力$\kappa\nabla\cdot Q$对缺陷动力学的类电场驱动机制，支持活性材料中信息传输的理论构想[69]。(4) 局限性：结果仅适用于$a > 0$（即$c > c_\star$），对$a \leq 0$无法保证全局存在；衰减估计要求初值在$L^1$中，且对$Q$的最高阶导数（$k = s$）衰减速率降为$(1+t)^{-s/2}$，非最优；无法处理$H^2 \times H^1$初值（二维可实现）；仅适用于常数活性$c = \mathrm{const}$，非常数活性系统未被涵盖。

4. Euler 奇异性 II：全空间 $\mathbb R^3$ 中带旋涡的光滑轴对称 Euler 方程的内部四极爆破

Euler Singularities II: Interior Quadrupole Blow-Up for Smooth Axisymmetric Euler with Swirl in \texorpdfstring{$\mathbb R^3$}

作者: Rishad Shahmurov

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04526v1

摘要翻译: 我们为全空间中带旋涡的三维轴对称不可压 Euler 方程构造了一种自洽的、位于内部区域的四极型有限时间奇点形成机制。该构造远离对称轴。在以被追踪的径向点为中心的局部变量

Comment: 无

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

1. 研究对象和背景

研究对象为三维不可压缩Euler方程的轴对称带旋流解：

在柱坐标系$(r, \theta, z)$下，轴对称速度场表示为：

引入关键变量$\Gamma = r u^\theta$和$G = \omega^\theta/r$，其中$\omega = \nabla \times u$为涡度场。对于光滑轴对称Euler流，满足：

这里$D_t = \partial_t + u^r \partial_r + u^z \partial_z$为物质导数，$\omega^\theta$为方位角涡度分量。轴对称散度自由条件为：

方位角涡度满足$\omega^\theta = \partial_z u^r - \partial_r u^z$。通过五维提升技术，引入标量势$\phi$满足$-\Delta_5 \phi = G$，其中$\Delta_5 = \partial_r^2 + 3r^{-1}\partial_r + \partial_z^2$为五维拉普拉斯算子，速度场可表示为：

这种提升源于将三维轴对称问题视为$\mathbb{R}^4 \times \mathbb{R}$中五维拉普拉斯方程的轴对称商空间。

本文考虑局部变量$x = r - r_*(t)$, $y = z$，其中$r_*(t) > 0$为追踪的内部径向中心，满足$r'_*(t) = u^r(r_*(t), 0, t)$。在该坐标系下，主动涡度和旋流剖面为：

第一个剖面产生正的内部Biot-Savart双曲应变，第二个剖面使$G$方程中的源项再生相同四极形状。全空间问题与边界问题有本质区别：缺乏壁面反射和单侧Dirichlet核符号，使得构造爆破机制更加困难。本文试图填补的理论空白是：在远离轴的全空间区域构造有限时间奇点，证明光滑解在轴对称带旋流类中也可能发生爆破，从而挑战光滑解全局存在的猜想。该研究对理解Beale-Kato-Majda准则和Constantin-Fefferman-Majda几何观点（涡旋拉伸与对齐机制的竞争）具有重要意义。

2. 主要定理或主要结果

在主引导假设$E(t) = \delta_{\text{jet}}(t) + \mu(t) + R_{\text{prof}}(t) + \rho(t) + \varepsilon_{\text{strain}}(t) \ll 1$下，其中$\mu(t) = b(t)\lambda(t)^3/\Gamma_*(t)$, $\rho(t) = \lambda(t)/r_*(t)$，证明了以下主要结果：

1. 应变下界（引理6.1）：当$\lambda(t) \ll r_*(t)$时，

   $$ \sigma(t) \geq cQ(t) - C\left(R_{\text{prof}}(t) + \frac{\lambda(t)}{r_*(t)}\right)Q(t), \quad \sigma(t) = -\partial_z u^z(r_*(t), 0, t) $$

   特别地，当主缺陷足够小时，$\sigma(t) \geq cQ(t)$。

2. 比较系统（引理7.1-7.2）：在主引导过程中，

   $$ c_0 C(t) \leq Q'(t) \leq C_0 C(t), \quad C'(t) \geq c Q(t)C(t), \quad C(t) \geq \kappa Q(t)^2 $$

   其中$Q(t) = H^{\text{full}}_\lambda[G](t)$为完整四象限四极分数，$C(t) = \lambda(t)^2 b(t)$，$b(t) = \partial_x \partial_y^2 \Gamma(0,0,t)$。若初始满足$C(0) \geq \kappa Q(0)^2$且$0 < \kappa < c/(4C_0)$，则源优势$C(t) \geq \kappa Q(t)^2$在整个过程中保持。

3. 显式初始数据（定理9.1）：当$b_0 = A_b a_0^2 \lambda_0^2$（$A_b \gg 1$）且$\lambda_0$足够小满足$A_b a_0^2 \lambda_0^5 \leq \varepsilon_0 \Gamma_{*,0}$, $\lambda_0/r_0 \leq \varepsilon_0$时，显式构造的光滑初始数据

   $$ G_0(r,z) = a_0(r-r_0)z \chi_{\lambda_0}(r-r_0,z), \quad \Gamma_0(r,z) = \chi_{\lambda_0}(r-r_0,z)\left(\Gamma_{*,0} + \frac{1}{2}b_0(r-r_0)z^2\right) $$

   满足$E(0) \leq C\varepsilon_0$且$C(0) \geq \kappa Q(0)^2$，从而进入四极引导区域。

4. 有限时间爆破：由比较系统可得$Q'(t) \geq c\kappa Q(t)^2$，导致$Q(t)$在有限比较时间内爆破。结合应变下界$\sigma(t) \geq cQ(t)$，得到$\|\nabla u(t)\|_{L^\infty}$的爆破，即存在有限时间$T_* < \infty$使得$\lim_{t \to T_*^-} \|\nabla u(t)\|_{L^\infty} = \infty$。

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

主要技术障碍包括：全空间问题缺乏边界反射机制和单侧Dirichlet核符号；涡度与速度场的强耦合导致$G$方程中源项$r^{-4}\partial_z(\Gamma^2)$的非线性再生；需保持四极剖面形状的稳定性以防止扰动破坏反馈循环；移动坐标系中的中心漂移、尺度变化和截断误差的精确控制。

作者的真正创新点：

1. 五维提升技术：将三维轴对称问题嵌入五维空间，通过$-\Delta_5 \phi = G$（$\Delta_5 = \partial_r^2 + 3r^{-1}\partial_r + \partial_z^2$）建立椭圆恢复公式。在$(Y,z) \in \mathbb{R}^4 \times \mathbb{R}$中，$r = |Y|$，五维牛顿势$|(Y,z)|^{-3}$的混合导数给出内部四极核

   $$ K_{\text{int}}(x,y) = C_0 \frac{xy}{(x^2 + y^2)^2} + K_{\text{rem}}(x,y), \quad C_0 > 0 $$

   其中余项$|K_{\text{rem}}| \leq C(\lambda/r_*) |xy|/(x^2 + y^2)^2 + K_{\text{smooth}}$。该符号结构是证明应变下界$\sigma(t) \geq cQ(t)$的关键。

2. 完整四象限四极分数与强制子分数：定义完整分数$Q(t) = \int_{P_\lambda} K_Q G \, dxdy$作为主动量，$P_\lambda = \{|x|,|y| < \lambda\}$，同时引入窄对角扇区$D_{\lambda,\delta_c} = \{|x|,|y| < \lambda, \left||y|/|x| - 1\right| < \delta_c\}$的子分数作为强制估计。通过引理4.1建立比较：

   $$ H^{\text{diag}}_{\lambda,\delta_c}[G] \geq c_{\delta_c} H^{\text{full}}_\lambda[G] - C D_{\text{ang},\lambda}[G] $$

   其中$D_{\text{ang},\lambda}[G] = \inf_a \int_{P_\lambda} W_Q |G - axy| \, dxdy$为角剖面缺陷，$W_Q = |xy|/(x^2 + y^2)^2$。

3. 移动包加权投影：为避免点态系数的局限性，定义加权内积$\langle F, H \rangle_\lambda = \int F H w_\lambda \, dxdy$，$w_\lambda = \chi_\lambda^2$，$\chi_\lambda(x,y) = \psi(x/\lambda)\psi(y/\lambda)$。投影振幅为：

   $$ a_\lambda(t) = \frac{\langle G, xy \rangle_\lambda}{\langle xy, xy \rangle_\lambda}, \quad b_\lambda(t) = \frac{2\langle \Gamma - \Gamma_*, xy^2 \rangle_\lambda}{\langle xy^2, xy^2 \rangle_\lambda} $$

   对应归一化量$Q_\lambda = c_Q a_\lambda \lambda^2$, $C_\lambda = \lambda^2 b_\lambda$。移动投影交换子（引理5.1）精确量化了移动中心、收缩尺度和截断的误差：

   $$ \left|\text{Comm}_{\lambda,r_*}[F, \Phi]\right| \leq C\sigma(t)(\varepsilon_{\text{strain}} + \rho) \int_{P_{2\lambda}} |F| |\Phi| \, dxdy $$

4. 关键引理的精细估计：

• 引理6.2（源兼容性）：对$\Gamma = \Gamma_* + \frac{1}{2}bxy^2 + R_\Gamma$，证明
  $$ r^{-4}\partial_y(\Gamma^2) = 2r^{-4}_* \Gamma_* b xy + \text{Err}_{\text{src}}, \quad |\text{Err}_{\text{src}}| \leq C\left(\delta_{\text{jet}} + \frac{\lambda}{r_*} + \frac{b\lambda^3}{\Gamma_*}\right) r^{-4}_* \Gamma_* b |xy| $$
• 引理6.3（旋流喷流放大）：证明$b'(t) \geq (1 - C[\delta_{\text{jet}} + \varepsilon_{\text{strain}} + \lambda/r_*]) \sigma(t)b(t)$，表明$b(t)$随应变率指数增长。
• 引理7.3（角剖面缺陷控制）：导出微分不等式
  $$ D^+ R_{\text{prof}} \leq -c\frac{C}{Q} R_{\text{prof}} + C\varepsilon_{\text{adm}}\frac{C}{Q} + C\sigma R_{\text{prof}}^2, \quad \varepsilon_{\text{adm}} = \delta_{\text{jet}} + \mu + \lambda/r_* + \varepsilon_{\text{strain}} $$
  确保剖面缺陷在引导过程中受控。

5. 显式数据构造与自进入验证：构造远离轴的光滑初始数据$G_0, \Gamma_0$，通过直接计算验证当$b_0 = A_b a_0^2 \lambda_0^2$（$A_b \gg 1$）时，初始满足$C(0) \geq \kappa Q(0)^2$且$E(0) \ll 1$，从而进入四极引导区域。证明框架为：构造初始数据→验证自进入→证明主引导传播→推导比较系统→得出有限时间爆破。

4. 与之前工作的比较

与前人工作的比较：

1. 与边界爆破工作的对比：Luo-Hou[6,7]和Chen-Hou[9]在周期圆柱体中证明了边界奇点形成，依赖壁面反射和单侧Dirichlet核符号。本文处理全空间问题，缺乏这些机制，开发了全新的’内部四极机制’。边界爆破中，奇点发生在轴上或近轴区域；而本文构造的奇点完全在远离轴的内部区域（$\lambda_0 \ll r_0$）。技术上，边界问题利用一维符号结构，而全空间需要处理二维四极符号结构，通过五维提升揭示$K_{\text{int}}(x,y) = C_0 xy/(x^2 + y^2)^2$的符号性质。适用范围上，本文机制不依赖边界，但要求初始数据远离轴；边界爆破机制适用于有界区域，但无法推广到全空间。

2. 与Elgindi工作的对比：Elgindi[10]证明了$C^{1,\alpha}$解的有限时间奇点，但解在轴上不光滑（$\omega^\theta \sim r^{-\beta}$）。本文处理全空间光滑解（$C^\infty$且远离轴），证明了在真正光滑类中也可能形成奇点。Elgindi的工作依赖轴上奇异性，而本文构造完全避免轴上区域，技术路线完全不同：Elgindi使用自相似解和谱分析，本文采用移动包引导论证。代价是：Elgindi的结果适用于更广的非轴对称扰动，而本文仅限于轴对称带旋流类；但本文结果更符合物理直觉，因为真实流体通常具有光滑初始数据。

3. 与作者第一部分论文的对比：第一部分处理周期圆柱体中的边界爆破，使用’边界双曲包’机制；本文是第二部分，开发’内部四极包’机制。两种机制的反馈循环不同：边界情况下，$G \sim a y \Rightarrow \sigma \gtrsim Q$, $\Gamma \sim \Gamma_* + b x y \Rightarrow Q' \gtrsim C$, $C' \gtrsim Q C$；而内部情况下，$G \sim a xy \Rightarrow \sigma \gtrsim Q$, $\Gamma \sim \Gamma_* + b x y^2 \Rightarrow Q' \gtrsim C$, $C' \gtrsim Q C$。内部机制需要额外控制径向反射误差$\lambda/r_*$，且源项再生要求更高阶的$xy^2$剖面。时间尺度上，内部爆破可能比边界爆破更慢，因为$C' \gtrsim Q C$与$C \geq \kappa Q^2$导致$Q' \gtrsim \kappa Q^2$，而边界情况下$C' \gtrsim Q^2$导致$Q' \gtrsim Q^{3/2}$。

4. 与Navier-Stokes工作的联系：作者提到两篇Navier-Stokes手稿[1,2]，使用相同五维提升技术但发展’第一阈值连续性机制’来证明全局存在性。这两类工作互补：Navier-Stokes中粘性抑制包放大，而Euler中无粘几何允许包放大驱动奇点形成。关键差异在于：Navier-Stokes利用$\nu \Delta u$项的强制性，而本文利用无粘反馈循环。适用范围上，Navier-Stokes结果要求初始能量低于阈值，而本文构造要求初始满足$b_0 = A_b a_0^2 \lambda_0^2$；代价是Navier-Stokes结果适用于一般光滑数据，而本文仅适用于特定构造的初始数据。

5. 技术推进与局限：本文引入’命名估计’系统（内部四极核符号展开、源兼容性等），将复杂引导论证模块化。相比前人，本文：

• 首次在全空间光滑类中构造有限时间奇点
• 开发五维提升技术处理内部符号结构
• 使用完整四象限分数而非窄扇区，增强鲁棒性
• 通过移动包投影精确控制交换子误差

局限包括：仅适用于轴对称带旋流类，尚未推广到一般三维Euler；构造的奇点是’内部’但仍在有限距离处；证明依赖精细引导论证，技术复杂度高。与边界爆破相比，本文机制的物理意义可能较不直观，但理论价值在于证明全空间光滑解也可能爆破，挑战了’光滑解全局存在’的猜想。

6. 时间尺度与爆破率：由比较系统$Q' \geq c\kappa Q^2$，得$Q(t) \geq (T_* - t)^{-1}$，$T_* < \infty$。应变率$\sigma(t) \geq cQ(t)$导致$\|\nabla u\|_{L^\infty} \gtrsim (T_* - t)^{-1}$，符合Beale-Kato-Majda准则。与Elgindi的$C^{1,\alpha}$解（爆破率$\sim (T_* - t)^{-\alpha}$）相比，本文光滑解的爆破率更接近临界情况。与边界爆破相比，内部机制的时间尺度更长，因为$C \geq \kappa Q^2$要求更强的初始条件（$b_0 \sim a_0^2 \lambda_0^2$ vs 边界情况$b_0 \sim a_0 \lambda_0$）。

7. 理论意义：本文为三维Euler方程的奇点形成问题提供了新视角，表明：

• 轴对称带旋流类中，全空间光滑解也可能有限时间爆破
• 无粘几何中的反馈循环可独立于边界机制驱动奇点形成
• 五维提升技术为分析轴对称问题提供了新工具

这些结果支持了涡旋拉伸在缺乏对齐机制时可能导致奇点的观点，对理解流体奇点形成的普适机制具有重要意义。

5. 背景磁场的稳定化作用：带速度阻尼的可压等熵理想磁流体动力学方程组的整体适定性

Stabilization by a background magnetic field: global well-posedness of the compressible isentropic ideal MHD equations with velocity damping

作者: Liening Qiao, Jiahong Wu, Fuyi Xu, Xiaoping Zhai

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04462v1

摘要翻译: 我们研究定义在周期环面 $\mathbb{T}^3$ 上的三维等熵可压理想（无粘、无电阻）磁流体动力学（MHD）方程组在速度阻尼项下的柯西问题。该系统存在一个稳态平衡态，由常数密度 $\barρ$ 和均匀背景磁场 $ω\in\mathbb{R}^3$ 构成。我们证明该平衡态是非线性稳定的。更精确地，若初始数据在 Sobolev 空间 $H^N(\mathbb{T}^3)$ 中为 $(\barρ,\mathbf{0},ω)$ 的充分小扰动，其中 $N\geq 6r+4$，且 $ω$ 满足 Diophantine 条件，则该系统存在唯一的整体光滑解。此外，扰动关于时间代数衰减。据我们所知，这是首个关于多维等熵可压理想 MHD 系统的整体适定性结果。证明揭示了一种隐含的耗散机制：尽管密度方程与磁场方程均不显含扩散或阻尼项，但通过背景磁场 $ω$ 耦合的速度与磁场之间的相互作用，结合 Diophantine–Poincaré 不等式，为密度扰动与磁场扰动共同生成了有效耗散；该耗散与速度阻尼协同作用，保证了解的整体正则性与时间衰减性。

Comment: 无

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

1. 研究对象和背景

本研究聚焦于三维周期环面 $\mathbb{T}^3$ 上的可压缩等熵理想磁流体力学（MHD）方程组，其速度场带有阻尼项。系统方程为

其中 $\rho$ 为密度，$u$ 为速度场，$H$ 为磁场，$\kappa > 0$ 为阻尼系数，等熵压力 $P = P(\rho)$ 是光滑严格递增函数。该系统描述了无粘（$\mu = \lambda = 0$）且无电阻（$\nu = 0$）的理想导电流体动力学，物理背景涵盖等离子体物理、天体物理吸积盘、液态金属流动等领域。

系统存在稳态解 $(\rho, u, H) \equiv (\bar{\rho}, 0, \omega)$，其中 $\bar{\rho} > 0$ 为常数密度，$\omega \in \mathbb{R}^3 \setminus \{0\}$ 为均匀背景磁场。研究核心问题是：当初始数据为 $(\bar{\rho}, 0, \omega)$ 的小扰动时，系统是否存在唯一的全局光滑解，且解是否保持接近稳态？

引入扰动变量 $a := \rho - \bar{\rho}$，$B := H - \omega$（并设 $\bar{\rho} = 1$，$\kappa = 1$），系统可重写为

其中 $\beta := P'(1) > 0$，非线性项为

线性化系统揭示了关键物理机制：密度扰动满足阻尼声波方程

而磁场扰动满足含阻尼的波动方程

此处 $(\omega \cdot \nabla)^2$ 项表征沿背景场 $\omega$ 以阿尔文速度 $|\omega|$ 传播的阿尔文波，而 $\partial_t B$ 项则体现了速度耦合产生的间接阻尼，这是物理稳定机制的数学体现。

在周期环面 $\mathbb{T}^3$ 上，波算子 $(\omega \cdot \nabla)$ 在傅里叶模 $e^{i k \cdot x}$ 上的作用为乘以 $i(\omega \cdot k)$。为使相应的庞加莱型不等式有效，需 $\omega \cdot k$ 对所有非零格点向量 $k$ 远离零，即要求 $\omega$ 满足 Diophantine 条件：存在 $c > 0$ 和 $r > 2$ 使得

由此导出的关键不等式为

其中 $f$ 为零均值函数。该不等式将波算子的控制转化为 Sobolev 范数估计，是揭示隐藏耗散的核心工具。

研究此问题的动机源于：一方面，无粘无电阻理想 MHD 是数学上最具挑战性的模型；另一方面，物理实验和数值模拟均表明，足够强的均匀背景磁场可抑制湍流并稳定导电流体。然而，在无粘无电阻设置下，该系统的全局适定性在多维情形长期悬而未决。本文填补了这一空白，首次证明了多维可压缩等熵理想 MHD 系统的全局适定性，揭示了背景磁场与速度阻尼协同作用产生的隐藏耗散机制。

2. 主要定理或主要结果

在以下假设下证明了主要结论：设 $N \geq 6r + 4$ 且 $r > 2$，背景磁场 $\omega \in \mathbb{R}^3$ 满足 Diophantine 条件 $|\omega \cdot k| \geq c / |k|^r$（$\forall k \in \mathbb{Z}^3 \setminus \{0\}$，$c > 0$ 为常数）。初始扰动 $(\rho_0 - \bar{\rho}, u_0, B_0)$ 满足均值条件 $\int_{\mathbb{T}^3} \rho_0 \, dx = \bar{\rho}$，$\int_{\mathbb{T}^3} \rho_0 u_0 \, dx = \int_{\mathbb{T}^3} B_0 \, dx = 0$，且 $c_0 \leq \rho_0 \leq c_0^{-1}$（$c_0 > 0$）。则存在 $\varepsilon > 0$，使得当 $\|\rho_0 - \bar{\rho}\|_{H^N} + \|u_0\|_{H^N} + \|B_0\|_{H^N} \leq \varepsilon$ 时，系统 (1.5) 存在唯一全局解 $(\rho - \bar{\rho}, u, B) \in C([0, \infty); H^N)$。此外，对任意 $t \geq 0$ 和 $r + 4 \leq \gamma \leq N$，扰动满足代数衰减估计

关键参数关系：正则性要求 $N \geq 6r + 4$ 源于 Diophantine 条件导致的 $r$ 阶导数损失，衰减速率 $-3(N-\gamma)/(2(2r+2))$ 通过低阶能量衰减与高阶能量不等式的插值得到。证明中构造的 Lyapunov 泛函为

其中 $A > 0$ 为充分大的常数，满足 $\mathcal{E}(t) \sim \|(a, u, B)\|_{H^{4r+2}}^2$。该泛函满足耗散不等式 $d\mathcal{E}(t)/dt + c \|(a, u, B)\|_{H^{2r}}^2 \leq 0$，结合插值不等式导出 $\mathcal{E}(t) \lesssim (1 + t)^{-(N-4r-2)/(2r+2)}$，最终通过 Gronwall 不等式控制最高阶范数。

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

证明面临的核心障碍在于：(1) 连续性方程 $\partial_t a + \text{div}(\rho u) = 0$ 无直接作用于 $a$ 的扩散或阻尼项；(2) 感应方程为纯传输-拉伸方程，对 $B$ 无平滑性；(3) 阻尼仅作用于 $u$，导致 $\frac{d}{dt} \mathcal{E}(t) + \|u\|^2 \leq 0$ 中 $a$ 与 $B$ 的能量无法被控制；(4) 非线性项 $B \cdot \nabla B$、$B \cdot \nabla u$ 和 $u \cdot \nabla B$ 的耦合使标准能量方法失效。

作者的真正创新点是揭示了隐藏耗散机制：通过背景场 $\omega$ 与速度-磁场耦合，结合 Diophantine-Poincaré 不等式，为 $a$ 和 $B$ 生成有效耗散。具体技术路径如下：

首先，利用 Diophantine-Poincaré 不等式（引理 2.1）：对零均值函数 $f$，有 $\|f\|_{H^s} \leq C \|\omega \cdot \nabla f\|_{H^{s+r}}$。该不等式将波算子控制转化为 Sobolev 范数估计，但损失 $r$ 阶导数，导致高正则性要求 $N \geq 6r + 4$。

其次，通过投影技术提取隐藏耗散（引理 3.2-3.3）：

• 对密度 $a$，将动量方程投影到 $\omega$ 上分离 $\omega \cdot \nabla a$，应用 $\Lambda^s$ 并与 $\omega \cdot \nabla \Lambda^s a$ 做内积，得到
  $$ \|( \omega \cdot \nabla) \Lambda^s a\|_0^2 + \frac{d}{dt} \int \omega \cdot \nabla \Lambda^s a \cdot \Lambda^s (u \cdot \omega) \, dx \lesssim \|u\|_{s+1}^2 + \|B\|_{s+1}^2. $$
  结合 Diophantine-Poincaré 不等式，$\|( \omega \cdot \nabla) \Lambda^s a\|_0$ 转化为 $\|a\|_{s-r}$ 的控制，损失 $r$ 阶导数。
• 对磁场 $B$，将 $\Lambda^s$ 作用于动量方程并与 $(\omega \cdot \nabla) \Lambda^s B$ 做内积，得到
  $$ \|( \omega \cdot \nabla) \Lambda^s B\|_0^2 - \frac{d}{dt} \int \Lambda^s u \cdot (\omega \cdot \nabla) \Lambda^s B \, dx \lesssim \|u\|_{s+1}^2. $$
  同样应用 Diophantine-Poincaré 不等式，获得 $B$ 的 $H^{s-r}$ 控制。

第三，处理非线性项需精细估计：

• 使用 Kato-Ponce 交换子估计（引理 2.4）处理 $[\Lambda^s, f \cdot \nabla]g$ 型项，关键估计为
  $$ \|[ \Lambda^s, f \cdot \nabla]g\|_{L^2} \leq C \left( \|\nabla f\|_{L^\infty} \|\Lambda^s g\|_{L^2} + \|\Lambda^s f\|_{L^2} \|\nabla g\|_{L^\infty} \right). $$
• 利用散度自由条件 $\text{div} B = 0$ 发现结构抵消：对 $(I_5, I_8)$ 和 $(I_6, I_9)$ 对进行分部积分，消除高阶导数项。
• 通过加权 Poincaré 不等式（引理 2.6）处理 $\|u\|_{L^2}$ 的控制，利用 $\int_{\mathbb{T}^3} \rho u \, dx = 0$ 得到
  $$ \|u\|_{L^2} \leq K (1 + \|a\|_{L^2}) \|\nabla u\|_{L^2}, $$
  在小扰动假设下简化为标准 Poincaré 不等式。

最后，构建 Lyapunov 泛函 $\mathcal{E}(t)$ 并闭合能量层次：结合基本能量估计 $\frac{d}{dt} \int (\rho |u|^2 + |a|^2 + |B|^2) \, dx + \int \rho |u|^2 \, dx \leq 0$、隐藏耗散估计（$4r+2$ 阶）和高阶能量不等式（$N$ 阶），通过交叉项 $\int \omega \cdot \nabla \Lambda^s a \cdot \Lambda^s (u \cdot \omega) \, dx$ 和 $\int \Lambda^s u \cdot (\omega \cdot \nabla) \Lambda^s B \, dx$ 耦合各阶导数，最终导出 $\frac{d}{dt} \mathcal{E}(t) + c \|(a, u, B)\|_{H^{2r}}^2 \leq 0$。此框架将双波结构（声波与阿尔文波）的耦合效应转化为定量的 Sobolev 控制，突破了无粘无电阻设置下的技术瓶颈。

4. 与之前工作的比较

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

与粘性/电阻 MHD 系统的比较：对于完整粘性（$\mu, \lambda > 0$）和电阻（$\nu > 0$）的可压缩 MHD 系统 (1.13)，全局适定性已被充分研究 [10, 12, 17, 18, 25, 26, 27, 35, 29, 36, 37, 46, 47, 50]。然而，本文处理的是无粘无电阻的理想情形（$\mu = \lambda = \nu = 0$），其数学难度显著更高。在非电阻情形（$\nu = 0, \mu, \lambda > 0$），环面上已有工作 [9, 54, 56, 57, 31, 23] 利用背景磁场和 Diophantine 条件生成磁场有效耗散，但这些结果依赖于粘性提供的密度耗散机制（通过有效粘性通量 $F = \text{div} u - P$）。而本文完全摒弃粘性，揭示了仅靠背景磁场与速度阻尼即可为密度和磁场同时生成隐藏耗散的新机制。

与不可压缩 MHD 的比较：在不可压缩框架下，Lin–Zhang [32] 和 Lin–Xu–Zhang [33] 证明了粘性非电阻 MHD（$\mu > 0, \nu = 0, \text{div} u = 0$）在背景磁场附近的全局适定性；Bardos–Sulem–Sulem [2] 对无粘无电阻理想 MHD 使用 Elsässer 变量 $z^\pm = u \pm B$ 获得对称双曲系统，从而得到一致能量估计。然而，可压缩情形中 $\text{div} u \neq 0$ 破坏了对称性，Elsässer 方法失效。本文通过线性化系统揭示的双波结构（声波与阿尔文波）耦合机制，为可压缩情形提供了新视角。值得注意的是，Wu–Wu–Xu [55] 和 Jiang–Jiang–Zhao [22] 研究了带速度阻尼的不可压缩 MHD，但可压缩性引入的密度扰动使问题复杂化，本文首次将此类结果推广至可压缩情形。

与可压缩无粘 MHD 的比较：Wu–Xu–Zhai [53] 研究了非等熵方程的全局光滑解，Li–Qiao [30] 获得了无粘电阻等熵情形的结果，但二者均未覆盖完全理想（$\mu = \lambda = \nu = 0$）的等熵设置。本文首次填补了这一空白，证明了 $\mathbb{T}^3$ 上小初值解的全局存在性。特别地，当 $B_0 \equiv 0$ 时，系统退化为带速度阻尼的可压缩欧拉方程，Sideris–Thomases–Wang [45] 已证明其全局存在性和指数衰减。本文表明，添加非零背景磁场 $\omega$ 不仅与阻尼相容，还增强了稳定效果，使衰减率从指数型变为代数型（$\sim (1+t)^{-\alpha}$），这是磁场耦合导致的物理效应。

推进与代价：本文的主要推进在于首次建立多维可压缩理想 MHD 的全局适定性理论，并揭示了背景磁场诱导的隐藏耗散机制。关键假设差异在于：(1) 要求 $\omega$ 满足 Diophantine 条件（尽管对 Lebesgue 几乎所有 $\omega \in \mathbb{R}^3$ 成立）；(2) 需高正则性 $N \geq 6r + 4$（源于 Diophantine-Poincaré 不等式的 $r$ 阶导数损失）。相比粘性情形，本文的代价是衰减率从指数型降为代数型，且仅适用于小初值扰动。然而，这一代价换取了对无粘无电阻物理模型的严格数学理解，为等离子体物理中的磁约束稳定性问题提供了理论基础。

适用范围与局限：本文结果适用于周期域 $\mathbb{T}^3$ 上的小扰动解，但无法处理 $\mathbb{R}^3$ 上的情形（因缺乏紧性）。Diophantine 条件要求 $r > 2$，这与 $\mathbb{T}^3$ 上的维数相关；在更高维或非周期域上可能需要调整参数。此外，结果仅针对等熵情形，非等熵情形需额外处理温度方程。尽管如此，本文建立的隐藏耗散框架为后续研究提供了新工具，例如可推广至带边界问题或更复杂的 MHD 模型。

6. 具有大初值的多维等熵可压Navier–Stokes–Korteweg系统的全局球对称解

Global spherically symmetric solutions to the multidimensional isentropic compressible Navier–Stokes–Korteweg system with large initial data

作者: Zhengzheng Chen, Fanfan Jiang

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04386v1

摘要翻译: 本文研究无界外区域上多维等熵可压Navier-Stokes-Korteweg系统初边值问题的大初值球对称强解的全局存在性。我们考虑压力 $p(ρ)=ρ^γ$、黏性系数 $μ(ρ)$ 与 $ λ(ρ)$ 满足如下两种情形之一：$μ(ρ)=\tildeμ, λ(ρ)=\tildeλρ^α$ 或 $μ(ρ)=\tildeμρ^α, λ(ρ)=\tildeλρ^α$，以及毛细系数 $κ(ρ)=\tildeκρ^β$，其中 $α,β,γ\in \mathbb{R}$ 为参数，$\tildeμ,\tildeλ,\tildeκ$ 为给定实常数。在对参数 $α,β$ 与 $γ$ 施加适当限制的条件下，我们建立了球对称强解的全局存在性与唯一性。证明依赖于径向加权能量方法，并结合Y.~Kanel’\cite{28}发展的技术。

Comment: 无

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

1. 研究对象和背景

本文研究多维等熵可压缩Navier-Stokes-Korteweg系统在外域中的球对称解全局存在性问题，针对大初值情形。考虑无界外域 $\Omega := \{\xi \in \mathbb{R}^d : |\xi| \geq a\}$（$d \geq 2$ 为空间维数，$a > 0$ 为给定常数）中的初边值问题，系统方程为：

其中应力张量 $T$ 和毛细管张量 $K$ 分别为：

主要未知量为密度 $\rho(t, \xi) > 0$ 和速度场 $u(t, \xi) \in \mathbb{R}^d$。初边值条件为：

其中 $n$ 为 $\partial \Omega$ 的单位外法向量。考虑压力函数 $p(\rho) = \rho^\gamma$，粘性系数 $\mu(\rho)$、$\lambda(\rho)$ 和毛细管系数 $\kappa(\rho)$ 满足以下两种情形之一： (1) $\mu(\rho) = \tilde{\mu}$, $\lambda(\rho) = \tilde{\lambda} \rho^\alpha$, $\kappa(\rho) = \tilde{\kappa} \rho^\beta$； (2) $\mu(\rho) = \tilde{\mu} \rho^\alpha$, $\lambda(\rho) = \tilde{\lambda} \rho^\alpha$, $\kappa(\rho) = \tilde{\kappa} \rho^\beta$， 其中 $\alpha, \beta, \gamma \in \mathbb{R}$ 为参数，$\tilde{\mu}, \tilde{\lambda}, \tilde{\kappa}$ 为常数。

在球对称假设下，解可表示为 $(\rho(t, \xi), u(t, \xi)) = \left(\hat{\rho}(t, r), \frac{\xi}{r} \hat{u}(t, r)\right)$（$r := |\xi| \geq a$），相应系统转化为：

其中 $m = d - 1$。通过拉格朗日坐标变换，系统进一步简化为：

初边值条件为 $(v, u)|_{t=0} = (v_0, u_0)(x)$，$(v_x, u)|_{x=0} = 0$，且 $\lim_{x \to +\infty} (v_0, u_0) = (1, 0)$。

Navier-Stokes-Korteweg系统描述了具有内部毛细现象的粘性可压缩流体运动，源于Van der Waals和Korteweg的工作，并由Dunn和Serrin在二阶梯度理论框架下严格推导。与经典Navier-Stokes方程相比，该系统因Korteweg应力张量引入了高度非线性项，动量方程包含密度的三阶导数，在非等熵情况下还包含温度的二阶导数，分析难度显著增加。该系统在相变、接触角问题和量子流体模型中有重要应用。

已有研究主要集中在小初值情形：Danchin和Desjardins在临界空间中建立全局解；Haspot扩展到非等熵情况；Hattori和Li证明Sobolev空间中光滑解的全局存在性；Bresch等获得三维周期情形下的全局弱解；Tsyganov等研究一维模型。然而，球对称解的研究几乎空白，特别是大初值、多维外域情形。本文填补了这一空白，首次研究多维外域中大初值的球对称解全局存在性问题。

2. 主要定理或主要结果

论文证明了两个主要定理：

定理1.1（Kazhikhov模型）：设粘性系数 $\mu(v), \lambda(v)$、热传导系数 $\kappa(v)$ 和压力函数 $p(v)$ 满足(1.18)，参数 $\alpha, \beta, \gamma$ 满足以下条件之一： (i) $\alpha \in \mathbb{R}$, $-3 \leq \beta \leq -2$, $\gamma \geq 1$； (ii) $\frac{\beta + 2 - \sqrt{-2(\beta + 2)(\beta + 5)}}{2} < \alpha \leq \beta + 3$, $-\frac{7 + \sqrt{3}}{2} \leq \beta < -3$, $\gamma > -\beta - 2$； (iii) $-\frac{3 + \sqrt{-2\beta^2 - 22\beta - 59}}{2} < \alpha \leq \beta + 3$, $-5 \leq \beta \leq -\frac{14}{3}$, $\gamma > -\beta - 2$。

若初值 $(v_0, u_0)(x)$ 满足：

其中 $V, \bar{V} > 0$ 为任意常数，则初边值问题(1.15)-(1.17)存在唯一全局强解 $(v, u)(t, x)$，满足对任意 $T > 0$：

其中 $C_1$ 仅依赖于 $\alpha, \beta, \gamma, V, \bar{V}, \|v_0 - 1\|, \|u_0\|, \|v_{0x}\|_{0,r}$；$C_2(T)$ 仅依赖于 $T, \alpha, \beta, \gamma, V, \bar{V}, \|v_0 - 1\|_1, \|u_0\|_1, \|v_{0x}\|_{1,r}$。

定理1.2（密度依赖粘性模型）：设参数满足(1.19)及以下条件之一： (i) $\alpha \in \mathbb{R}$, $-3 \leq \beta \leq -2$, $\gamma \geq 1$； (ii) $\alpha = \frac{\beta + 3}{2}$, $\beta < -3$, $\gamma > -\beta - 2$； (iii) $\frac{\beta + 2 - \sqrt{-2(\beta + 2)(\beta + 5)}}{2} < \alpha \leq \frac{\beta + 4}{3}$, $-4 \leq \beta < -3$, $\gamma > -\beta - 2$； (iv) $\frac{\beta + 2 - \sqrt{-2(\beta + 2)(\beta + 5)}}{2} < \alpha < \frac{\beta + 2 + \sqrt{-2(\beta + 2)(\beta + 5)}}{2}$, $-5 \leq \beta < -4$, $\gamma > -\beta - 2$； (v) $-\frac{3 + \sqrt{-2\beta^2 - 22\beta - 59}}{2} < \alpha < -\frac{3 - \sqrt{-2\beta^2 - 22\beta - 59}}{2}$, $-\frac{11 + \sqrt{3}}{2} \leq \beta \leq -\frac{11 - \sqrt{3}}{2}$, $\gamma > -\beta - 2$。

则初边值问题(1.15)-(1.17)存在唯一全局强解 $(v, u)(t, x)$，满足：

其中 $C_3, C_4(T)$ 的依赖关系类似定理1.1。值得注意的是，定理1.2包含Bresch-Desjardins (BD)熵结构情形（$\tilde{\mu} = 1$, $\tilde{\lambda} = 2(\alpha - 1)$, $\alpha > 0$）。

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

主要技术难点在于推导密度的均匀上下界及处理非线性高阶项。作者的核心创新点如下：

1. 基本能量估计：通过乘以 $p(1) - p(v)$ 和 $u$ 并相加，得到关键能量等式：

   $$ \left(\Phi(v) + \frac{u^2}{2} + \frac{r^{2m}}{2v^{\beta+5}} v_x^2 \right)_t + \left(2\tilde{\mu} + \frac{\tilde{\lambda}}{v^\alpha}\right) \frac{[(r^m u)_x]^2}{v} = R_{1x}, $$

   其中 $\Phi(v) = \int_1^v (p(1) - p(s)) ds$，$R_1$ 为边界项。积分后得：

   $$ \int_{\mathbb{R}^+} \left(\Phi(v) + \frac{u^2}{2} + \frac{r^{2m}}{2v^{\beta+5}} v_x^2 \right) dx + \int_0^t \int_{\mathbb{R}^+} \left(2\tilde{\mu} + \frac{\tilde{\lambda}}{v^\alpha}\right) \frac{[(r^m u)_x]^2}{v} dx d\tau = \epsilon_0. $$

   此估计中 $\int_{\mathbb{R}^+} \frac{r^{2m}}{v^{\beta+5}} v_x^2 dx$ 项对后续分析至关重要。

2. Kanel’方法的应用与改进：定义 $\Upsilon(v) = \int_1^v \sqrt{\Psi(s)} s^{\frac{\beta+5}{2}} ds$（$\Psi(v) = \frac{(1-v)^2}{1+v}$），利用基本能量估计和Cauchy-Schwarz不等式：

   $$ |\Upsilon(v)(t, x)| \leq \left\|\sqrt{\Psi(v)(t)}\right\| \cdot \left\|\frac{v_x}{v^{\frac{\beta+5}{2}}}(t)\right\| \leq C \left\|\sqrt{\Phi(v)(t)}\right\| \cdot \left\|\frac{v_x}{v^{\frac{\beta+5}{2}}}(t)\right\| \leq C. $$

   结合 $\Upsilon(v)$ 在 $v \to 0$ 和 $v \to \infty$ 时的渐近行为：

   $$ |\Upsilon(v)| \geq \begin{cases} A_1 \left(v^{-\frac{\beta+3}{2}} + \ln \frac{1}{v}\right) - A_2, & v \to 0, \\ A_1 \left(v^{-\frac{\beta+2}{2}} + \ln v\right) - A_2, & v \to \infty, \end{cases} $$

   证明当 $-3 \leq \beta \leq -2$ 时，$C_5^{-1} \leq v(t, x) \leq C_5$。该方法将经典Kanel’技术扩展到含毛细项的高维系统。

3. 关键高阶估计：通过乘以 $r^{-m} \frac{v_x}{v^{\alpha+1}}$ 得到：

   $$ \left(\frac{\tilde{\lambda} v_x^2}{2v^{2\alpha+2}} + \frac{2\tilde{\mu} v_x^2}{v^{\alpha+2}} - \frac{u v_x}{r^m v^{\alpha+1}}\right)_t + \frac{\gamma v_x^2}{v^{\alpha+\gamma+2}} + \frac{r^{2m}}{v^{\alpha+\beta+6}} v_{xx}^2 = R_{2x} + H_1, $$

   其中 $H_1$ 包含关键非线性项：

   $$ H_1 = m(m-1)\frac{v_x^2}{v^{\alpha+\beta+4} r^2} + \frac{[(r^m u)_x]^2}{r^{2m} v^{\alpha+1}} - \frac{2m u (r^m u)_x}{r^{2m+1} v^\alpha} + \frac{m u^2 v_x}{r^{m+1} v^{\alpha+1}} + \cdots - f_1(\alpha, \beta) \frac{r^{2m}}{v^{\alpha+\beta+8}} v_x^4, $$

   $f_1(\alpha, \beta) = \frac{(\beta + 5)(\alpha + \beta + 7)}{6} + \frac{(\alpha + 1)(\beta + 5)}{2} - \frac{(\alpha + \beta + 6)(\alpha + \beta + 7)}{3}$。对于不同参数范围，作者通过Sobolev不等式和Gronwall不等式控制 $\int_0^t \int_{\mathbb{R}^+} \frac{r^{2m}}{v^{\alpha+\beta+8}} v_x^4 dxd\tau$ 项，从而获得关键估计：

   $$ \int_{\mathbb{R}^+} \left(\frac{v_x^2}{v^{2\alpha+2}} + \frac{v_x^2}{v^{\alpha+2}}\right) dx + \int_0^t \int_{\mathbb{R}^+} \left(\frac{v_x^2}{v^{\alpha+\gamma+2}} + \frac{r^{2m}}{v^{\alpha+\beta+6}} v_{xx}^2 \right) dxd\tau \leq C_6(T), $$

   或

   $$ \int_{\mathbb{R}^+} \left(\frac{v_x^2}{v^{2\alpha+2}} + \frac{v_x^2}{v^{\alpha+2}}\right) dx + \int_0^t \int_{\mathbb{R}^+} \left(\frac{v_x^2}{v^{\alpha+\gamma+2}} + \left[\left(\frac{r^m}{v^{\frac{\alpha+\beta+6}{2}}} v_x \right)_x \right]^2 \right) dxd\tau \leq C_7(T). $$

4. 参数范围的精细分析：作者发现不同模型中 $\frac{v_x}{v^{\alpha+1}}$ 估计的差异导致参数范围不同。在Kazhikhov模型中，需处理 $\int_0^t \int_{\mathbb{R}^+} \left(\tilde{\mu} \frac{(r^m u)_x v_x^2}{v^{\alpha+3}} - \tilde{\mu} \frac{(r^m u)_x v_{xx}}{v^{\alpha+2}}\right) dxd\tau$，要求 $\alpha \leq \beta + 3$；而密度依赖模型中相应项 $\int_0^t \int_{\mathbb{R}^+} \left(m\tilde{\mu} \frac{u}{v^\alpha r} \left(\frac{v_x}{v^{\alpha+1}}\right)_x - m\tilde{\mu} \frac{(r^m u)_x v_x}{r^{m+1} v^{2\alpha+1}}\right) dxd\tau$ 导致 $\alpha \leq \frac{\beta+4}{3}$，从而使得 $\alpha = \frac{\beta+3}{2}$ 在密度依赖模型中可行。

5. 径向加权能量方法：作者发展了径向加权能量方法，结合Kanel’技术处理外域中球对称解的全局存在性。通过引入加权空间 $L^p_r(\mathbb{R}^+) = \left\{f \in L^p(\mathbb{R}^+) \mid \|f\|_{L^p_r} < \infty \right\}$，其范数为 $\|f\|_{L^p_r} = \left(\int_{\mathbb{R}^+} r^{2m}|f(x)|^p dx \right)^{1/p}$，以及加权Sobolev空间 $H^k_r(\mathbb{R}^+)$，有效控制了高维几何效应。

4. 与之前工作的比较

本文与已有工作的比较主要体现在以下几个方面：

1. 研究对象的创新性：如Remark 1.2指出，这是首次研究Navier-Stokes-Korteweg系统的球对称解。已有研究主要集中在小初值情形（Danchin-Desjardins [6]、Haspot [7,8]、Hattori-Li [9,10]）、一维模型（Tsyganov [20]、Charve-Haspot [21]、Germain-LeFloch [22]）或周期域（Bresch-Desjardins-Lin [5]）。本文首次处理大初值、多维外域中的球对称解，填补了研究空白。

2. 与小初值结果的对比：已有工作多要求初值足够小，而本文允许任意大的初值。例如，Danchin-Desjardins在临界空间中建立小初值的全局解，Haspot扩展到非等熵情况，但都限于小初值。本文通过Kanel’方法和径向加权能量估计，克服了大初值带来的困难。

3. 与经典Navier-Stokes方程结果的对比：对于经典可压缩Navier-Stokes方程，球对称解已有较多研究（Nikolaev [29]、Chen-Kratka [30]、Frid-Shelukhin [31]、Hoff-Jenssen [32]），但这些结果都依赖于域的有界性。本文考虑的是无界外域，且包含Korteweg应力张量，分析难度更大。特别是，本文需要处理动量方程中的三阶导数项，这是经典Navier-Stokes方程所没有的。

4. 与Chen等人[23]的结果比较：在[23]中，Chen等人建立了大初值扰动下的全局经典解，但仅获得时间依赖的上下界 $C(T)^{-1} \leq v(t, x) \leq C(T)$。相比之下，本文不仅建立了外域中球对称解的全局存在性，还推导了均匀（时间）上下界 $C_0^{-1} \leq v(t, x) \leq C_0$，这是显著改进。尽管本文对参数范围有轻微限制，但获得了更强的估计结果。

5. 参数范围的差异：定理1.1和1.2的参数范围差异源于模型结构的不同。在Kazhikhov模型中，估计 $\left\|\frac{v_x}{v^{\alpha+1}}\right\|$ 时需处理 $\int_0^t \int_{\mathbb{R}^+} \left(\tilde{\mu} \frac{(r^m u)_x v_x^2}{v^{\alpha+3}} - \tilde{\mu} \frac{(r^m u)_x v_{xx}}{v^{\alpha+2}}\right) dxd\tau$，要求 $\alpha \leq \beta + 3$，导致当 $\alpha = \frac{\beta+3}{2}$ 时无可行参数范围；而在密度依赖模型中，相应项导致条件 $\alpha \leq \frac{\beta+4}{3}$，使 $\alpha = \frac{\beta+3}{2}$ 成为可能。这表明密度依赖粘性模型在参数适应性上更具优势。

6. 与BD熵结构的联系：定理1.2包含Bresch-Desjardins (BD)熵结构情形（Remark 1.3）。取 $\tilde{\mu} = 1$ 和 $\tilde{\lambda} = 2(\alpha - 1)$ 时，满足BD约束 $2\mu(\rho) + d\lambda(\rho) > 0$（需 $\alpha > 0$）。此时参数范围为：(i) $\alpha > 0$, $-3 \leq \beta \leq -2$, $\gamma \geq 1$；(ii) $0 < \alpha < \frac{\beta+4}{3}$, $-4 < \beta < -3$, $\gamma > -\beta - 2$。BD熵结构提供了额外的先验估计，有助于处理真空状态，但本文通过Kanel’方法在更广泛的情形下获得了结果。

7. 大时间行为的局限性：如Remark 1.7所述，本文无法建立多维等熵情况下的大时间行为，主要受限于高维项 $\int_0^t \int_{\mathbb{R}^+} \frac{m(m-1) v_x^2}{v^{\alpha+\beta+4} r^2} dxd\tau$、$\int_0^t \int_{\mathbb{R}^+} \frac{2m u (r^m u)_x}{r^{2m+1} v^\alpha} dxd\tau$ 和 $\int_0^t \int_{\mathbb{R}^+} \frac{m u^2 v_x}{r^{m+1} v^{\alpha+1}} dxd\tau$。相比之下，一维情况（$d = 1$，即 $m = 0$）可获得大时间渐近行为。这表明高维几何效应是大时间行为分析的主要障碍。

8. 域的限制：本文专注于外域（$r \geq a > 0$）中的解，而包含原点的域（$r \geq 0$）中的问题仍开放（Remark 1.9）。这与经典Navier-Stokes方程的研究形成对比，后者在外域和包含原点的域中都有较多结果。Korteweg系统在原点附近的奇性更强，需要新的技术处理。

9. 技术方法的推进：本文将Kanel’方法[28]与径向加权能量估计相结合，这是对已有技术的重要拓展。Kanel’方法最初用于一维可压缩Navier-Stokes方程，本文将其推广到含毛细项的高维球对称系统，处理了额外的三阶导数项和几何效应。同时，通过精心设计的加权函数和参数范围分析，克服了外域中衰减估计的困难。

10. 与Cao-Li-Zhu工作的对比：Cao、Li和Zhu[44]研究了球外域中三维球对称正则解的全局存在性，但仅限于经典Navier-Stokes方程（无Korteweg项）和真空初值。本文将结果推广到含Korteweg项的系统，并处理了更一般的粘性系数和毛细管系数，同时允许大初值但不假设真空状态。

7. 可压缩流体–黏弹性壳耦合系统的爆破准则与弱–强唯一性

Blow-Up Criteria and Weak–Strong Uniqueness for Compressible Fluid–Viscoelastic Shell Interactions

作者: Pierre Marie Ngougoue Ngougoue, Prince Romeo Mensah

类别: math.AP

发布日期: 2026-05-05

arXiv ID: 2605.04314v1

摘要翻译: 近期已在有限时间区间上建立了等熵可压缩流体–黏弹性壳耦合系统强解的存在性与唯一性。一个自然的问题是：此类解是否可全局延拓？本文中，我们为该耦合系统导出了一条延拓准则。我们的分析基于在物质加速度层次上的能量估计，该估计在 Serrin 型与 Beale–Kato–Majda 型控制假设下得以建立。尽管在不可压情形下，此类控制已足以防止有限时间爆破，但在可压情形下，它本身并不能保证强解所需的全部正则性的传播。为获得真正的延拓准则，我们对密度与速度梯度施加了 Beale–Kato–Majda 型 Lipschitz 控制，并要求其具有更强的时间可积性。结合支撑加速度估计的控制框架，我们闭合了一个高阶能量估计，从而防止强解正则性的丧失。因此，只要相应的控制范数保持有限，解便可延拓至潜在爆破时刻之后。此外，我们在上述条件正则性准则下进一步建立了该系统的弱–强唯一性原理。

Comment: 无

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

1. 研究对象和背景

本文研究可压缩流体与粘弹性壳体相互作用系统的爆破准则与弱-强唯一性。研究对象为耦合系统：流体部分由可压缩Navier-Stokes方程描述，壳体部分由粘弹性壳方程控制，两者通过运动学边界条件耦合。具体方程如下：

流体方程（在移动域 $\Omega_\eta(t) \subset \mathbb{R}^3$ 上）：

其中 $v: (t, x) \in I \times \Omega_\eta \mapsto v(t, x) \in \mathbb{R}^3$ 和 $\rho: (t, x) \in I \times \Omega_\eta \mapsto \rho(t, x) \in \mathbb{R}$ 分别表示流体速度和密度，压力 $p(\rho) = a\rho^\gamma$（$a > 0, \gamma > 1$）为等熵状态方程。粘性系数满足物理约束 $\mu > 0$, $\lambda + \frac{2}{3}\mu \geq 0$。

壳体方程（在参考曲面 $\omega \subset \mathbb{R}^2$ 上）：

其中 $\phi_\eta: \omega \to \partial \Omega_\eta$ 是变形边界的参数化，$\det(\nabla_y \phi_\eta) = |\partial_1 \phi_\eta \times \partial_2 \phi_\eta|$ 为曲面Jacobi行列式，$n$ 和 $n_\eta$ 分别为参考边界和变形边界的外法向量，$\tau = S(\nabla v) - p(\rho)I_{3\times 3}$ 为流体Cauchy应力张量，$S(\nabla v)$ 为粘性应力张量。

耦合条件（运动学边界条件）：

该问题的重要性在于：流体占据的域 $\Omega_\eta(t)$ 依赖于壳体位移 $\eta(t, \cdot)$，导致方程在移动域上求解，增加了分析难度。[21]最近建立了该系统在有限时间区间上的局部强解存在唯一性，但自然的问题是：这些解能否全局延拓？在固定域的可压缩Navier-Stokes方程中，已知爆破通常与密度或速度梯度的奇性相关（如[27, 11]），但在移动边界配置中，界面几何的退化也可能导致解的爆破（如[3]在不可压缩情形所示）。本文试图填补的空白是：在可压缩FSI系统中，确定哪些分析或几何量控制可能的正则性损失，建立严格的延拓准则，并证明弱-强唯一性原理。特别地，与不可压缩情形不同，可压缩流体的连续性方程和非线性压力需要重新规范化处理，且Serrin型控制本身不足以保证强解正则性的传播，这构成了本文的主要挑战。

2. 主要定理或主要结果

在以下假设下证明了主要结论：

(A1) (Serrin型动量控制) 速度满足可积条件 $\rho^{1/2} v \in L^s(I^*; L^r(\Omega_\eta))$，其中 $\frac{2}{s} + \frac{3}{r} \leq 1$，$r \in (3, \infty]$，$s \in [2, \infty)$。

(A2) (压缩性效应控制) 满足以下至少一个条件：(a) $\|\text{div} v\|_{L^1(I^*; L^\infty(\Omega_\eta))} < \infty$；(b) $\|\rho\|_{L^\infty(I^*; L^\infty(\Omega_\eta))} < \infty$。

(A3) (结构几何正则性) 壳体位移满足 $\eta \in L^\infty(I^*; C^1(\omega))$，且流体-结构界面保持非退化。

(B) (Beale-Kato-Majda型Lipschitz控制) 满足：$\int_{I^*} \|\nabla v\|_{L^\infty(\Omega_\eta)}^2 dt < \infty$ 和 $\int_{I^*} \|\nabla \rho\|_{L^\infty(\Omega_\eta)}^2 dt < \infty$。

(CC) (兼容性条件) $[\partial_t v + v \cdot \nabla v] \circ \phi_\eta|_{t=0} = (\partial_t^2 \eta) n|_{t=0}$ on $I \times \omega$。

主要结果如下：

定理1.2（加速度估计）：设 $(\rho, v, \eta)$ 为定义1.1意义下的强解，若假设(A1)-(A3)成立，则有加速度能量估计：

其中隐含常数仅依赖于 $\|\rho^{1/2} v\|_{L^s(I; L^r(\Omega_\eta))}$, $\|\text{div} v\|_{L^1(I^*; L^\infty(\Omega_\eta))}$, $\|\eta\|_{L^\infty(I^*; C^1(\omega))}$ 和弱解先验能量估计常数 $C_0 > 0$（见(2.1)式）。

定理1.3（高阶能量估计）：在假设(A1)-(A3)、(B)和(CC)下，有高阶能量估计：

其中隐含常数仅依赖于(B)、(A1)-(A3)、$T^*$和加速度能量 $E_{\text{acc}}$。

定理1.4（延拓准则）：在假设(A1)-(A3)、(B)和(CC)下，有：

其中常数 $C > 0$ 仅依赖于(A1)-(A3)、(B)、$T^*$和 $C_0$。

推论1.1（延拓准则）：设 $(\rho, v, \eta)$ 为定义1.1意义下的最大强解，若兼容性条件(CC)成立，且假设(A1)-(A3)和(B)在 $I^*$ 上保持有效，则解可延拓到 $T^*$ 之后。

定理1.5（弱-强唯一性）：设初始数据 $(\rho_0, v_0, \eta_0, \eta_*)$ 满足(1.4)-(1.6)，且 $(\rho, v, \eta)$ 为[15, Section 2, Definition 2.1]意义下的有限能量弱解（绝热指数 $\gamma > \frac{12}{7}$）。若在每个有限区间 $I = (0, T)$ 上假设(A1)-(A3)、(B)和(CC)成立，则 $(\rho, v, \eta)$ 在 $I$ 上是定义1.1意义下的强解，且在具有 $\eta \in L^\infty(I; C^1(\omega))$ 的弱解类中唯一。

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

证明的核心障碍源于可压缩流体的连续性方程和非线性压力与移动边界的耦合效应。在不可压缩情形[3]中，Serrin型控制足以防止有限时间爆破，但在可压缩情形下，仅靠(A1)-(A3)无法保证强解正则性的传播，因为密度变化引入了额外的非线性项。具体而言，加速度估计(1.8)虽然提供了系统加速度的自然能量控制，但不足以恢复定义1.1所需的全部正则性（如 $\rho \in W^{1,\infty}(I^*; W^{2,2}(\Omega_\eta))$ 和 $v \in L^2(I^*; W^{4,2}(\Omega_\eta))$）。

作者的真正创新点在于：识别出在可压缩FSI系统中需要额外的Beale-Kato-Majda型Lipschitz控制（假设(B)），并证明该控制与加速度估计框架结合可闭合高阶能量估计。具体技术包括：

1. 材料加速度水平的能量估计：测试流体动量方程(1.1)₂与边界相容的材料导数型测试函数 $\psi := \partial_t v + E_\eta(\partial_t \eta n) \cdot \nabla v$，其中 $E_\eta$ 为通过Hanzawa变换构造的边界延拓算子。这一选择确保 $\psi \circ \phi_\eta = (\partial_t^2 \eta) n$，与壳体方程的测试函数匹配。关键步骤是将 $\psi$ 重写为 $\dot{v} + (E_\eta(\partial_t \eta n) - v) \cdot \nabla v$（其中 $\dot{v} = \partial_t v + v \cdot \nabla v$ 为材料导数），从而分离出 $\|\rho^{1/2} \dot{v}\|_{L^2(\Omega_\eta)}$，该量在后续分析中用于控制 $\nabla^2 v$。

2. 非线性项的精细估计：对(3.4)式中的15项 $R_k$ 进行逐项分析。以 $R_1 = \int_{I^*} \int_{\Omega_\eta} \rho |v \cdot \nabla v|^2 dx dt$ 为例，通过Hölder不等式与Gagliardo-Nirenberg插值 $\|\nabla v\|_{L^q(\Omega_\eta)} \lesssim \|v\|_{W^{2,2}(\Omega_\eta)}^{3/r} \|v\|_{W^{1,2}(\Omega_\eta)}^{(2r-6)/r}$（其中 $q = 2r/(r-2)$），得到：

   $$ R_1 \lesssim \int_{I^*} \|\rho^{1/2} v\|_{L^r(\Omega_\eta)}^2 \|v\|_{W^{2,2}(\Omega_\eta)}^{6/r} \|v\|_{W^{1,2}(\Omega_\eta)}^{(2r-6)/r} dt. $$

   应用Young不等式（指数 $s = 2r/(r-3) \in [2, \infty)$）后：

   $$ R_1 \lesssim \varepsilon \int_{I^*} \|v\|_{W^{2,2}(\Omega_\eta)}^2 dt + c(\varepsilon) \int_{I^*} \|\rho^{1/2} v\|_{L^r(\Omega_\eta)}^s \|v\|_{W^{1,2}(\Omega_\eta)}^2 dt, $$

   其中第一项可被左边吸收，第二项通过Grönwall不等式由Serrin型控制控制。类似地，边界项 $R_2$ 和 $R_7$ 利用迹嵌入 $W^{1/2,2}(\partial \Omega_\eta) \hookrightarrow L^4(\partial \Omega_\eta)$ 和分数阶Sobolev空间的插值性质 $W^{3/4,2}(\Omega_\eta) = [L^2(\Omega_\eta), W^{1,2}(\Omega_\eta)]_{3/4}$，结合 $\|\partial_t \eta \circ \phi_\eta^{-1}\|_{W^{1/4,2}(\partial \Omega_\eta)} \lesssim \|\partial_t \eta\|_{W^{1/4,2}(\omega)}$，得到：

   $$ R_2 \lesssim \varepsilon \int_{I^*} \|\nabla v\|_{W^{1,2}(\Omega_\eta)}^2 dt + c(\varepsilon) \int_{I^*} \|\partial_t \eta\|_{W^{1,2}(\omega)}^2 \|\nabla v\|_{L^2(\Omega_\eta)}^2 dt. $$

3. 高阶正则性的传播：在假设(B)下，通过控制 $\|\nabla v\|_{L^\infty(\Omega_\eta)}$ 和 $\|\nabla \rho\|_{L^\infty(\Omega_\eta)}$，作者闭合了高阶能量估计。关键步骤是利用 $\int_{I^*} \|\nabla v\|_{L^\infty(\Omega_\eta)}^2 dt < \infty$ 保证特征线稳定性（见[1, 17]），从而传播密度的高阶正则性。例如，在估计压力项 $R_9$ 时，应用Sobolev嵌入 $W^{3/2,2}(\Omega_\eta) \hookrightarrow L^4(\Omega_\eta)$ 得到：

   $$ R_{9,1} \lesssim \int_{I^*} \|\nabla v\|_{W^{1,2}(\Omega_\eta)} \|p(\rho)\|_{L^4(\Omega_\eta)} \|E_\eta(\partial_t \eta n)\|_{W^{3/2,2}(\Omega_\eta)} dt, $$

   其中 $\|p(\rho)\|_{L^4(\Omega_\eta)} \lesssim \|\rho\|_{L^{4\gamma}(\Omega_\eta)}^\gamma$，而 $\|\rho\|_{L^{4\gamma}(\Omega_\eta)}$ 可通过控制 $\|\nabla \rho\|_{L^\infty(\Omega_\eta)}$ 与初始数据估计。

4. 弱-强唯一性的证明：通过构造相对熵泛函，利用延拓准则中的控制条件，证明任何有限能量弱解在强解存在区间上必须与强解重合。核心是处理 $\int_{I^*} \int_{\Omega_\eta} \rho |v - \tilde{v}|^2 dx dt$ 项，其中 $\tilde{v}$ 为弱解速度，利用Serrin型控制和Grönwall不等式导出唯一性。

4. 与之前工作的比较

与既有PDE/分析理论文献的比较如下：

1. 与不可压缩FSI的比较：

• [3]在不可压缩流体耦合粘弹性壳的情形下，证明了Ladyzhenskaya-Prodi-Serrin型条件正则性结果，仅需Serrin型控制 $v \in L^s(I; L^r(\Omega))$（$\frac{2}{s} + \frac{3}{r} \leq 1$）和几何控制 $\eta \in L^\infty(I; C^1(\omega))$ 即可延拓强解并建立弱-强唯一性。
• 本文将结果推广到可压缩情形，但发现仅靠(A1)-(A3)不足以保证强解正则性的传播。关键差异在于：在可压缩情形下，加速度估计(1.8)无法直接控制密度的高阶导数，需要额外的Beale-Kato-Majda型Lipschitz控制（假设(B)）。具体而言，在不可压缩情形中 $\text{div} v = 0$ 简化了压力项的处理，而可压缩情形中 $\nabla p(\rho)$ 依赖于 $\nabla \rho$，必须控制 $\int_{I^*} \|\nabla \rho\|_{L^\infty(\Omega_\eta)}^2 dt < \infty$ 才能闭合高阶估计。
• 适用范围：本文结果要求更强的正则性条件（如 $\rho \in W^{3,2}(\Omega_{\eta_0})$），但允许 $\gamma > 1$ 的更广绝热指数范围（对比[4]的 $\gamma > \frac{12}{7}$）。

2. 与固定域可压缩Navier-Stokes的比较：

• [27]对三维Navier-Stokes系统证明了Beale-Kato-Majda爆破准则：在 $7\mu > \lambda$ 下，$\sup_{0 \leq t < T} \|\rho\|_{L^\infty(\Omega)} < \infty$ 足以延拓强解。其论证结合了有效粘性通量估计与Lamé算子的对数界。
• [11]将Serrin条件扩展到粘性可压缩情形：若 $v$ 满足(1.7)且 $\rho$ 有界或 $\text{div} v \in L^1(I; L^\infty(\Omega))$，则无爆破发生。在 $7\mu > \lambda$ 或无真空（$\min_{\mathbb{R}^3} \rho_0 > 0$）下，Serrin条件可省略。
• 本文推进：在移动边界配置中，不仅需要控制 $\rho$ 和 $v$，还需确保界面几何正则性（$\eta \in L^\infty(I; C^1(\omega))$）。特别地，假设(B)中的 $\int_{I^*} \|\nabla v\|_{L^\infty(\Omega_\eta)}^2 dt < \infty$ 比[11]的 $\text{div} v \in L^1(I; L^\infty(\Omega))$ 更强，因为移动域中 $\text{div} v$ 无法直接控制 $\nabla v$ 的全部分量。
• 代价：本文需要同时控制 $\nabla v$ 和 $\nabla \rho$ 的 $L^\infty$ 范数（假设(B)），而固定域情形中 $\rho$ 有界（假设(A2)(b)）可能已足够（见[11]）。

3. 与弱解理论的比较：

• [4]建立了可压缩FSI系统的弱解理论（$\gamma > \frac{12}{7}$），使用四层逼近方案（人工压力、人工粘性、边界正则化、Galerkin离散化），但仅得到有限能量解，未解决高阶正则性和唯一性问题。
• [21]证明了局部强解存在唯一性，但未提供延拓准则。
• 本文推进：在[21]的强解框架下，建立了严格的延拓准则（推论1.1）和弱-强唯一性原理（定理1.5）。关键突破是识别出可压缩FSI系统中爆破的充分条件：当控制范数 $\|\rho^{1/2} v\|_{L^s L^r}$、$\|\text{div} v\|_{L^1 L^\infty}$、$\|\eta\|_{L^\infty C^1}$、$\int \|\nabla v\|_{L^\infty}^2 dt$、$\int \|\nabla \rho\|_{L^\infty}^2 dt$ 保持有限时，解可全局延拓。
• 局限：本文结果要求初始数据具有高Sobolev正则性（$\rho_0 \in W^{3,2}$, $v_0 \in W^{3,2}$, $\eta_0 \in W^{5,2}$, $\eta_* \in W^{3,2}$），而弱解理论[4]仅需有限能量（$\rho_0 \in L^\gamma$, $\rho_0^{1/2} v_0 \in L^2$, $\eta_0 \in W^{2,2}$, $\eta_* \in L^2$）。

4. 时间尺度与适用范围：

• 在不可压缩FSI[3]中，Serrin型控制的时间尺度为 $\frac{2}{s} + \frac{3}{r} \leq 1$，与经典Navier-Stokes一致。
• 本文中，假设(A1)保持相同时间尺度，但假设(B)要求 $\int_{I^*} \|\nabla v\|_{L^\infty}^2 dt < \infty$（即 $L^2$-in-time控制），比固定域可压缩Navier-Stokes中常见的 $\int_{I^*} \|\nabla v\|_{L^\infty} dt < \infty$（$L^1$-in-time）更强。这是因为移动域中界面变形放大了速度梯度的效应，需要更严格的控制来防止几何退化。
• 适用范围：本文结果适用于 $\gamma > 1$ 的等熵流体，而弱解理论[4]要求 $\gamma > \frac{12}{7}$。然而，本文的延拓准则仅对强解有效，无法直接应用于弱解的全局存在性问题。

5. 与相关模型的比较：

• 对于结构由板或波方程描述的FSI模型（如[16]），正则性要求类似，但壳体方程(1.2)中的 $-\partial_t \Delta_y \eta$ 项引入了额外的耗散，简化了某些估计。
• 在[5]的FSI问题中（文中提及但未详细说明），可能需要不同的分析框架，因为其结构模型可能不包含粘性项 $-\partial_t \Delta_y \eta$。
• 与[9]的近期结果相比，本文处理了更复杂的移动边界效应，而[9]可能仅考虑固定域或简化几何。

8. Euler 奇异性 I：具涡旋的光滑精确奇对称 Euler 方程在边界的爆破

Euler Singularities I: Boundary Blow-Up for Smooth Exact-Odd Axisymmetric Euler with Swirl

作者: Rishad Shahmurov

类别: math.AP

发布日期: 2026-05-05

arXiv ID: 2605.04181v1

摘要翻译: 我们在周期圆柱区域中构造了光滑的轴对称带涡旋（axisymmetric-with-swirl）初始数据，使得三维不可压 Euler 方程的演化在有限时间内于边界处产生奇异性。该构造在动力学不变的精确奇类

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩Euler方程在周期圆柱域中的边界奇异性问题。考虑柱坐标系$(r, \theta, z)$下的轴对称带旋流速度场$u(r, z, t) = u^r(r, z, t)e_r + u^\theta(r, z, t)e_\theta + u^z(r, z, t)e_z$，其满足不可压缩条件$\nabla \cdot u = 0$和动量方程$\partial_t u + (u \cdot \nabla)u + \nabla p = 0$。定义关键变量：角动量$\Gamma = r u^\theta$和涡度比$G = \omega^\theta/r$，其中$\omega = \nabla \times u$为涡度场。这些变量满足动力学方程：

其中$D_t = \partial_t + u^r \partial_r + u^z \partial_z$为子午面物质导数。本文聚焦于精确奇对称类（exact-odd class）：

该类在Euler演化下保持不变。问题设定在周期圆柱$\Omega = \{(r, \theta, z) : 0 < r < 1, z \in \mathbb{T}\}$中，侧壁满足不可渗透边界条件$u^r(1, z, t) = 0$。在侧壁点$(r, z) = (1, 0)$处，精确奇对称性给出逐点恒等式：

其中$\sigma(t) = -\partial_z u^z(1, 0, t)$。研究该问题的重要性在于：有限时间奇点的存在性是数学流体动力学的核心问题，Beale-Kato-Majda准则指出最大涡度的时间可积性是决定性判据。Luo和Hou的数值工作[12,13]提供了边界奇点的数值证据，Chen和Hou最近给出了计算机辅助的边界奇点结果[4]，但缺乏对原始轴对称变量中直接核-包机制的结构化理解。本文试图填补的空白是：构建一个自包含的簇不变证明框架，展示超压缩二进展开角簇如何导致光滑初始数据在有限时间内形成边界奇点，而不依赖数值模拟或简化模型假设。

2. 主要定理或主要结果

在精确奇对称光滑初始数据假设下，证明了存在有限时间边界奇点。定理1.1（光滑精确奇对称侧壁Euler爆破）：存在满足以下条件的光滑轴对称带旋流初始数据：在周期圆柱$\Omega = \{0 < r < 1, z \in \mathbb{T}\}$中，满足不可渗透边界条件$u^r(1, z, 0) = 0$，且满足精确奇对称性$\Gamma(r, -z, 0) = -\Gamma(r, z, 0)$, $G(r, -z, 0) = -G(r, z, 0)$，使得相应的经典Euler解在有限时间$T^* < \infty$内满足

关键证明步骤包括：(i) 侧壁Dirichlet参数化给出主导压缩核$K_0(x, y) = C_0 \frac{xy}{(x^2 + y^2)^2}$ ($C_0 > 0$)及余项控制；(ii) 精确奇对称性导出核心ODE系统：

其中$a(t) = \partial_z G(1, 0, t)$, $b(t) = \partial_z \Gamma(1, 0, t)$；(iii) 定义尺度$\lambda(t)$满足$\lambda'(t) = -\sigma(t)\lambda(t)$，并令$A(t) = \lambda(t)a(t)$, $B(t) = \lambda(t)^{1/2}b(t)$，则有

该比较系统在有限时间内爆破，且$B(t) \to \infty$蕴含$\|\nabla u(t)\|_{L^\infty} \to \infty$。更精细地，通过窄对角极大分值论证得到

其中$M_{\delta c}$为锥分值。最终规范化引擎为

其中$A_*$和$B_*$是选定的均匀相干分量的振幅，该系统在有限时间内爆破。证明中关键尺度关系为：在锥$K_{\kappa\lambda} = \{0 < x < \kappa\lambda, mx < y < Mx\}$上，若满足锥非退化条件$c_0 a(t)y \leq G(1-x, y, t) \leq C_0 a(t)y$，则局部双曲质量$H_\lambda[G]$满足$c\lambda(t)a(t) \leq H_\lambda(t)[G(t)] \leq C\lambda(t)a(t)$，即$H_\lambda(t)[G(t)] \sim A(t)$。

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

核心障碍包括：(1) 非局部Biot-Savart律的复杂性，需要精确估计侧壁附近的Green核；(2) 碎片化问题，局部包可能分裂成同尺度的角碎片；(3) 远场尾部控制，远离活跃包的贡献可能破坏局部机制；(4) 剪切抑制，径向-旋流缺陷$R = \partial_r \Gamma / \partial_z \Gamma$可能旋转旋流梯度；(5) 曲率控制，$a'(t) = \sigma(t)a(t) + 2b(t)^2 + 2\Gamma_* c(t)$中可能的负项$2\Gamma_* c(t)$。作者的创新点与关键技术：(i) 五维提升与侧壁Dirichlet参数化：将轴对称恢复写成五维提升形式$-\Delta_5 \phi = G$，其中$\Delta_5 = \partial_r^2 + 3r^{-1}\partial_r + \partial_z^2$，侧壁边界条件$\phi|_{r=1} = 0$。在展平坐标$x = 1-r$, $y = z$下，Green核展开为$K_\Omega(1-x, y) = C_0 \frac{xy}{(x^2 + y^2)^2} + K_{\text{rem}}^\lambda(x, y)$，其中余项在固定锥$mx \leq |y| \leq Mx$上满足$|K_{\text{rem}}^\lambda(x, y)| \leq C\lambda \frac{|xy|}{(x^2 + y^2)^2} + K_{\text{smooth}}^\lambda(x, y)$。关键观察是：对奇符号一致包$G(1-x, -y) = -G(1-x, y)$, $G(1-x, y) \geq 0$ ($y > 0$)，局部包贡献产生正双曲压缩。(ii) 超压缩二进展开角簇功能：定义二进展开尾质量$A_j = P_{2^{j+1}\lambda} \setminus P_{2^j \lambda}$和$H_j[G] = \left| \int_{A_j} K_\Omega G \, d\mu_5 \right|$，其中$d\mu_5 = r^3 dr dz$为五维提升测度。证明超双曲尾支配或重新选择定理：若活跃边界包$P_\lambda$通过首次阈值最大性选择，则$\text{Tail}_\lambda[G] \leq \eta H_\lambda[G]$，否则存在尺度可比的竞争包，这与首次阈值最大性矛盾。(iii) 包-核耦合：在锥$K_{\kappa\lambda}$上，若$G$满足$c_0 a(t)y \leq G(1-x, y, t) \leq C_0 a(t)y$，则$H_\lambda(t)[G(t)] \sim A(t)$，其中$A(t) = \lambda(t)a(t)$。证明利用了在锥$y \sim x$上$K_0(x, y) \sim x^{-2}$, $G \sim a(t)x$，故$K_0 G \sim a(t)/x$，积分得$\lambda a(t)$。(iv) 离对角剪切抑制：证明在内核$x + |y| \leq \kappa\lambda$中，$|\partial_r u^z(1-x, y)| \leq C\kappa\sigma + \text{Tail}_{\text{shear}}[G] + \text{Smooth}_{\text{shear}}[G] + O(\lambda)\sigma$。利用Green核表示和奇对称性导致的奇偶性抵消，得到$|\partial_r u^z| \leq C \frac{x + |y|}{\lambda} H_\lambda[G] + \cdots$，结合$\sigma \gtrsim H_\lambda[G]$即得结果。进而证明若$|R(0)| \leq \varepsilon$且$\delta \ll \varepsilon \ll 1$，则$|R(t)| \leq 2\varepsilon$持久成立。(v) 分化源放大与曲率控制：利用边界速度$u^z(1, y, t)$的奇对称性（$v(y, t) = u^z(1, y, t)$为奇函数），证明曲率比$Q(t) = \Gamma_* \Gamma_{yy}(1, 0, t) / \Gamma_y(1, 0, t)^2$为常数。若$Q(0) > -1 + \varepsilon$，则$a'(t) \geq \sigma(t)a(t) + 2\varepsilon b(t)^2$，消除了可能的负项影响。(vi) Dini函数与仿射Campanato控制：超压缩二进展开角簇功能吸收同尺度角碎片、增长的二进展开窗口、动态分离的远尾和固定距离外部场，整合为可积小的仿射Campanato缺陷。结果不变簇包含振幅为$A_*(t)$和$B_*(t)$的均匀相干分量，满足Dini比较系统$D^+ A_*(t) \gtrsim B_*(t)^2$, $D^+ B_*(t) \gtrsim A_*(t) B_*(t)$，该系统在有限时间内爆破。

4. 与之前工作的比较

与Luo和Hou的数值工作[12,13]比较：Luo-Hou提供了边界奇点的数值证据，但缺乏严格数学证明；本文构建了光滑初始数据导致有限时间边界奇点的严格证明，机制理解更深入——通过精确核-包机制解释了物理循环：奇G包→边界压缩→$\partial_z \Gamma$变陡→$\partial_z(\Gamma^2)$源→更大的G包。优势在于提供自包含的结构性证明，不依赖数值模拟；局限在于需精确奇对称假设，而Luo-Hou研究更一般的轴对称情形。与Chen和Hou的边界奇点结果[4]比较：Chen-Hou使用计算机辅助的近自相似框架，本文采用直接核-包机制在原始轴对称变量中构建证明；本文更注重结构性，开发了自包含的簇不变证明方法；Chen-Hou结果可能更具普适性，但本文提供更清晰的机制理解，且证明不依赖计算机验证。与Hou-Luo型模型[7,6]比较：Hou-Luo型模型处理一维或简化模型，本文处理完整三维Euler方程；本文在原始Euler方程中构建奇点，无需简化假设；优势是直接适用于完整Euler方程，代价是需要精确奇对称假设。与Elgindi的有限时间奇点形成[10]比较：Elgindi处理$C^{1,\alpha}$解而非光滑解，本文处理光滑初始数据；Elgindi机制是内部奇点，本文是边界奇点；适用范围不同：Elgindi适用于无界域，本文限于周期圆柱边界；本文奇点机制更接近物理边界层现象。与Navier-Stokes相关工作[1,2]比较：作者提及两篇Navier-Stokes配套论文，其中五维提升、全Dirichlet可见性、粘性强制性和首次阈值延续机制用于抑制包配置集中；本文中因缺乏耗散，侧壁Dirichlet核创建正双曲反馈环；两系列互补：Navier-Stokes描述粘性如何防止包放大，Euler隔离无粘几何中包放大导致奇点；机制对比凸显粘性与无粘动力学的根本差异。与边界爆破构造的其他工作[5]比较：本文贡献在于超压缩二进展开角簇方法，提供新结构化工具；特别关注精确奇对称类机制；优势是方法可能更易推广。与小尺度创建机制[11]比较：[11]研究有界域中Euler的小尺度创建，本文专注边界奇点；机制差异在于本文依赖侧壁双曲压缩；本文提供更详细的Dini比较系统分析。本文推进：(i) 提供光滑精确奇对称轴对称Euler数据中边界奇点的直接构造；(ii) 开发超压缩二进展开角簇方法吸收各种误差；(iii) 证明Dini比较系统有限时间爆破；(iv) 明确连接Beale-Kato-Majda准则与局部机制。放宽假设：无需自相似假设（对比Chen-Hou）；使用精确奇对称类简化分析；侧壁条件比无滑移条件更弱（Euler不可渗透vs. Navier-Stokes无滑移）。适用范围：限于周期圆柱侧壁点$(r, z) = (1, 0)$；需精确奇对称初始数据；适用于光滑初始数据。代价与局限：需精确奇对称假设（物理上可能不自然）；奇点仅在特定点形成；方法依赖侧壁特殊几何，难推广到其他边界类型；作为系列第一部分仅处理边界奇点（第二部分处理全空间内部四极机制）。

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

1. b-诺维科夫方程在临界Sobolev空间中的不适定性

Ill-posedness in the critical Sobolev space for the $b$-Novikov equation

作者: Dan-Andrei Geba, A. Alexandrou Himonas, Curtis Holliman

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.05033v1

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

摘要翻译

中文标题: b-诺维科夫方程在临界Sobolev空间中的不适定性

摘要翻译: 本文证明了$b$-诺维科夫方程（一类具三次非线性的含单参数Camassa-Holm型方程）在临界Sobolev空间$H^{3/2}(\mathbb{R})$中存在范数膨胀现象。该结果完善了该方程的适定性理论：此前已知该方程在$H^{s}(\mathbb{R})$中对$s>3/2$局部适定，而在$H^{s}(\mathbb{R})$中对$s<3/2$不适定。

2. 一维拟线性波动方程经典解的生命跨度

Lifespan of Classical Solutions to One-Dimensional Quasilinear Wave Equations

作者: Yuusuke Sugiyama, Taro Yamanoi

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04976v1

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

摘要翻译

中文标题: 一维拟线性波动方程经典解的生命跨度

摘要翻译: 本文研究了一维拟线性波动方程 $u_{tt}-c(u_x)^2u_{xx}=0$ 的柯西问题的经典解生命跨度的上界与下界，其中 $c(θ)$ 的导数在原点附近趋于 $0$。特别地，我们的结果表明：解的生命跨度随初始数据的小性代数增长。此外，我们还证明：当 $c(θ)$ 在原点处平坦（flat）时，解的生命跨度随初始数据的小性指数增长。我们的证明基于 Lax 特征线方法与黎曼不变量。

3. Strichartz 极值函数的神经发现

Neural Discovery of Strichartz Extremizers

作者: Nicolás Valenzuela, Ricardo Freire, Claudio Muñoz

类别: math.AP, cs.LG, math.NA

发布日期: 2026-05-06

arXiv ID: 2605.04918v1

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

摘要翻译

中文标题: Strichartz 极值函数的神经发现

摘要翻译: Strichartz 不等式是现代色散型偏微分方程理论的基石，但其极值函数仅在少数几个 sharp 情形下被显式知晓。底层泛函的非凸性使该问题极为困难，且据我们所知，尚无系统性的数值尝试。我们提出一种基于神经网络的简单流程，将极值函数搜索转化为 Strichartz 比值 $\frac{\|e^{it\Delta}f\|_{L^q_t L^r_x}}{\|f\|_{L^2_x}}$ 的临界点搜寻，并将其应用于三种情形。第一，在 Schrödinger 群上，我们在 $d=1,2$ 维中以 $10^{-3}$ 量级的相对误差恢复了 Foschi 以及 Hundertmark–Zharnitsky 所得的 Gaussian 极值函数，且未引入任何解析先验信息。第二，在 $d=1$ 中其余 $59$ 个 admissible 对上（其极值函数形式仅为 conjectural），该方法始终收敛至 Gaussian 函数，从而支持如下 conjecture：Gaussian 函数是整个 admissible 范围内的 universal 极值函数。第三，在临界 Airy–Strichartz 不等式 $\gamma = 1/q$ 处（其极值函数的存在性仍 open），优化过程不收敛于任何 $L^2$ 剖面；相反，迭代序列自组织为 mKdV breathers $B(0,\cdot;\alpha,1,0,0)$，其内部频率 $\alpha$ 持续增长，且所发现的比值从下方趋近 Frank–Sabin universal 下界 $\widetilde A_{q,r}$，并呈现幂律间隙 $\sim \alpha^{-0.9}$。我们通过一个独立的 Hermite 基底 ansatz 验证了相同图像。我们提出一个精确 conjecture：该上确界等于 $\widetilde A_{q,r}$，且沿 breather 族被趋近但不可达。因此，该流程既可在已知情形下作为验证器（validator），亦可在极值函数不存在时作为发现工具（discovery tool）。

4. Zeldovich-von Neumann-Döring 燃烧模型活塞问题中爆轰波解的存在性

Existence of detonation wave solutions to the piston problem for the Zeldovich-von Neumann-D{ö}ring combustion model

作者: Xiaomin Zhang, Huimin Yu

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04836v1

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

摘要翻译

中文标题: Zeldovich-von Neumann-Döring 燃烧模型活塞问题中爆轰波解的存在性

摘要翻译: 本文研究具单步放热化学反应的一维 Zeldovich-von Neumann-Döring (ZND) 燃烧模型的活塞问题中的爆轰波解。作为一类特殊的激波，爆轰波的位置未知，因而该模型构成一个自由边界问题。我们证明了该自由边界问题的爆轰波解的整体存在性。与先前结果~\cite{Lai} 相比，我们未对控制方程及边界施加任何耗散性条件。

5. 具有变号非线性的双曲空间上的标量场方程

A scalar field equation on hyperbolic space with indefinite sign nonlinearity

作者: Debabrata Karmakar, Atanu Manna, Bhakti Bhusan Manna

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04687v1

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

摘要翻译

中文标题: 具有变号非线性的双曲空间上的标量场方程

摘要翻译: 本文研究双曲空间 $\mathbb{B}^N$（其中 $N \ge 3$）上半线性双幂椭圆方程

的阈值现象。对参数 $1 < p \le 2^*-1$（尽管我们偶尔允许超临界指数）与 $q > 0$，我们旨在刻画实参数 $λ\in \mathbb{R}$ 的最优谱区域，以区分正能量解的存在性与不存在性。我们在所有指数构型下——即 $p < q$、$0 < q < 1 < p$ 以及 $1 < q < p$——完全确定了这些阈值。我们的结果表明：分隔这些区域的边界由一个显式临界谱参数所控制；在 $p < q$ 区域中，该参数依赖于 $p$、$q$ 与 $N$；而在其余情形中，它仅依赖于 $N$。

6. 双曲空间中 Fisher-KPP 方程的行波行为

Traveling-wave behavior for Fisher-KPP equations in the hyperbolic space

作者: María del Mar González, Irene Gonzálvez, Fernando Quirós

类别: math.AP, math.DG

发布日期: 2026-05-06

arXiv ID: 2605.04661v1

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

摘要翻译

中文标题: 双曲空间中 Fisher-KPP 方程的行波行为

摘要翻译: 我们在双曲空间中研究带有 Fisher-KPP 型强迫项的热方程的柯西问题。解的长时间行为（传播或消失）取决于扩散与强迫项所诱导增长之间的相对强度，其中扩散强度由拉普拉斯–贝尔特拉米算子谱的下确界刻画。我们证明了关于这一二分性（dichotomy）的新结果，特别涵盖了扩散与反应项处于同阶的临界情形。若初始数据具有某种对称性（即在双曲空间等距群的一个余齐性一子群下不变），则该问题可约化为一维问题。在传播情形下，该一维问题的解在适当选取的移动坐标系中，其形状收敛于一个欧氏最小速度行波 $u(t,x) = U(x - ct)$；该坐标系的选择依赖于初始数据所保持不变的等距子群类型（椭圆型、双曲型或抛物型）。与欧氏情形不同，在适当坐标下的渐近传播速度依赖于空间维数，而波前位置中的对数修正项系数则与底层等距变换无关。

7. 多重特征值情形下的Riesz性质

Riesz property in the case of multiple eigenvalues

作者: Boris Mityagin, Petr Siegl

类别: math.SP, math-ph, math.AP, math.FA

发布日期: 2026-05-06

arXiv ID: 2605.04659v1

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

摘要翻译

中文标题: 多重特征值情形下的Riesz性质

摘要翻译: 我们分析了具有任意重数（包括无限重）特征值的自伴算子之非对称扰动的谱及其谱投影的Riesz性质。特别地，当 $d/2 < r < \infty$ 时，我们证明了复值 $L^r$-位势对多维谐振子、朗道哈密顿量以及球面上的拉普拉斯–贝尔特拉米算子所作扰动具有Riesz性质。

8. 高维 Ginzburg–Landau 型泛函极小元的渐近性

Asymptotics of Minimizers for Ginzburg–Landau-type Functionals in High Dimensions

作者: Giacomo Canevari, Haotong Fu, Wei Wang

类别: math.AP, math.DG

发布日期: 2026-05-06

arXiv ID: 2605.04442v1

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

摘要翻译

中文标题: 高维 Ginzburg–Landau 型泛函极小元的渐近性

摘要翻译: 我们研究高维 $n\geq 3$ 中满足对数型能量估计的 Ginzburg–Landau 型泛函的局部极小元，假设势函数的真空流形具有有限的基本群。我们证明：归一化能量测度收敛于一个与平稳 varifold 相关联的 $(n-2)$-可 rectifiable 测度，其密度是量子化的，由真空流形的同伦类决定。在该 $(n-2)$-可 rectifiable 测度的支集之外，极小元在 $H^1_{\text{loc}}$ 中强收敛于一个极小调和映射；该调和映射在某个 $(n-3)$-可 rectifiable 奇异集之外是光滑的。

9. 分数阶向量微积分与分数阶麦克斯韦方程组

Fractional Vector Calculus and the Fractional Maxwell’s Equations

作者: Giovanni Covi, Ruirui Wu

类别: math.AP

发布日期: 2026-05-05

arXiv ID: 2605.04241v1

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

摘要翻译

中文标题: 分数阶向量微积分与分数阶麦克斯韦方程组

摘要翻译: 我们考虑麦克斯韦方程组的一个分数阶变体，其中电场与磁场被建模为两点场（two-point fields）。为建立该系统，我们引入一个与分数阶散度算子 $\mathrm{div}^\alpha$ 相容的分数阶旋度算子 $\mathrm{curl}^\alpha$，从而保证无散条件。一个关键工具是一个投影映射 $Π$，它将两点场约化为单点场（one-point fields）。我们还定义了一个新的分数阶索博列夫空间，其元素具有分数阶亥姆霍兹分解；并观察到投影 $Π$ 在该空间中是双射，这使得我们能将整个问题完全重述为仅关于单点场的形式。随后，我们在加权分数阶索博列夫空间中证明了单点场形式方程组的适定性，并由此导出两点分数阶麦克斯韦系统的相应适定性结果。这构成了求解分数阶麦克斯韦方程组散射反问题的第一步必要工作，该反问题将是后续研究的主题。

不相关论文（基本信息）

1. 五个不等式中的可理解性

Grokability in five inequalities

作者: Paata Ivanisvili, Xinyuan Xie

类别: math.PR, cs.AI, math.AP, math.CA, math.FA

发布日期: 2026-05-06

arXiv ID: 2605.05193v1

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

摘要: In this note, we report five mathematical discoveries made in collaboration with Grok, all of which have been subsequently verified by the authors. These include an improved lower bound on the maximal Gaussian perimeter of convex sets in $\mathbb{R}^n$, sharper $L_2$-$L_1$ moment comparison inequalities on the Hamming cube $\{-1,1\}^n$, a strengthened autoconvolution inequality, improved asymptotic bounds on the size of the largest $g$-Sidon sets in $\{1,\dots,n\}$, and an optimal balanced Szarek’s inequality.

2. 部分观测下 Vlasov–Poisson 方程不稳定性控制的可证明模仿学习

Provable imitation learning for control of instability in partially-observed Vlasov–Poisson equations

作者: Xiaofan Xia, Qin Li, Wenlong Mou

类别: cs.LG, math.AP, math.OC, physics.plasm-ph

发布日期: 2026-05-06

arXiv ID: 2605.05081v1

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

摘要: We consider the stabilization of Vlasov–Poisson plasma dynamics, a central control problem in nuclear fusion. Our focus is the gap between what an ideal controller would use and what experiments can actually observe: while optimal policy may rely on the full phase-space state, practical feedback is typically limited to sparse macroscopic diagnostics. We therefore study imitation learning methods that distill a fully observed expert policy into controllers operating only on macroscopic measurements. We show the stability guarantees of the learned policy, where the error floor depends on the minimal behavior cloning loss achievable under the observation constraints. We further characterize this minimal loss in terms of a notion of entropy that quantifies the complexity of the initial distribution. Our results demonstrates the theoretical feasibility of learning stabilizing feedback policies for kinetic plasma dynamics from macroscopic observations, and exhibits the adaptivity of the learning approach to low-complexity structures. Through extensive numerical experiments, we validate our theory and show that the learned policies can stabilize the system using only macroscopic observations, within a significantly longer time horizon than non-adaptive baseline controllers.

3. 具有小奇异集的浸入稳定极小超曲面的最优正则性理论

An Optimal Regularity Theory for Immersed Stable Minimal Hypersurfaces with Small Singular Set

作者: Paul Minter, Zhengyi Xiao

类别: math.DG, math.AP

发布日期: 2026-05-06

arXiv ID: 2605.05041v1

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

摘要: We show that if $M^n$ is a properly immersed, two-sided, stable minimal hypersurface in $B^{n+1}_1(0)\setminus S$, where $S$ is closed with $\mathcal{H}^{n-2}(S)=0$, then $\text{dim}_{\mathcal{H}}\text{sing}(M)\leq n-7$, namely $\overline{M}\cap B^{n+1}_1(0)$ is represented by a smooth minimal immersion outside a closed set of generally unavoidable singularities which has Hausdorff dimension at most $n-7$. This provides the optimal a priori size assumption on the non-immersed singular set in order to guarantee optimal regularity. Consequently, such objects form a compact class under mass upper bounds.

4. 投影演化提升与无限维定常Hamilton-Jacobi-Bellman方程的适定性

Projected Evolutionary Lifting and Well-Posedness of Stationary Hamilton-Jacobi-Bellman Equations in Infinite Dimensions

作者: Gabriele Bolli, Fabian Fuchs

类别: math.OC, math.AP

发布日期: 2026-05-06

arXiv ID: 2605.05028v1

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

摘要: This paper establishes the existence and uniqueness of mild solutions to stationary Hamilton-Jacobi-Bellman (HJB) equations associated with infinite-horizon stochastic optimal control problems in separable Hilbert spaces. Our framework includes settings with a lack of global smoothing properties of the transition semigroup, singular dynamics involving unbounded control operators, and state-dependent running costs. We overcome these challenges by lifting the state space using the Projected Evolutionary Lifting technique. This work is an extension of G. Bolli and F. Gozzi, Lifting and partial smoothing for stationary HJB equations and related control problems in infinite dimensions, 2025, in which existence and uniqueness is proved via a contraction mapping argument and is consequently restricted to sufficiently large discount factors. We remove this restriction, proving existence and uniqueness for any discount rate $λ> 0$ using tools from the theory of maximally monotone operators.

5. Calderón 问题中非自伴各向异性及复值包含体的重构

Reconstruction of non-self-adjoint anisotropic and complex inclusions in the Calderón problem

作者: Henrik Garde, David Johansson, Thanasis Zacharopoulos

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.05021v1

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

摘要: We generalize recent results on the monotonicity method, for inclusion detection in the partial data anisotropic Calderón problem, to very general non-self-adjoint perturbations. This involves a forward model that accounts for both the anisotropic real conductivity and the anisotropic permittivity, and the results hold in any spatial dimension $d \geq 2$. We assume that the inclusion boundaries can be reached from the domain boundary via a set on which the background conductivity is self-adjoint, and that a definiteness condition holds near the inclusion boundaries. Away from the inclusion boundaries we allow general $L^\infty$ non-self-adjoint perturbations. We only require unique continuation based on the self-adjoint part of the background conductivity, thus making the methods compatible with generic unique continuation results.

6. 含低阶项的非线性抛物问题的重正规化解

Renormalized Solution for the Nonlinear Parabolic Problem with Lower Order Terms

作者: Shijun Li, Boai Huang, Shaopeng Xu

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04820v1

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

摘要: In this paper, we consider the following problem:

in a bounded open set $ Ω\subset \mathbb{R}^N $. We have established certain gradient estimates and proved the existence of a renormalized solution for the equation.

7. 一个微分方程的非唯一性及ChatGPT给出的证明

Non-uniqueness for a differential equation and a proof by ChatGPT

作者: Brian Street

类别: math.AP, math.CA

发布日期: 2026-05-06

arXiv ID: 2605.04810v1

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

摘要: Let $f(t,x),M(t,x)\in C([0,1]^2)$ with $M(t,x)>0$. We consider differential equations of the form

For a fixed positive weight $M$, we ask whether the condition $f(0,x)=0$ forces $f\equiv 0$. We show the answer is negative for smooth functions: there exist $f(t,x),M(t,x)\in C^{\infty}([0,1]^2)$ with $f(0,x)=0$, $f(t,0)\not\equiv 0$, and $M(t,x)>0$ satisfying the above equation. However, we show that for a large class of $M(t,x)$, the equation does have uniqueness. We relate this to uniqueness/non-uniqueness theorems for weighted Laplace transforms. A key example originated in an output by ChatGPT-5.5-Pro, and we include a discussion of its output as well as a complete proof.

8. 一类带对数扰动的Kirchhoff型方程的归一化解

Normalized solutions to a class of Kirchhoff type equations with a logarithmic perturbation

作者: Qi Li, Wenshu Zhou, Yuzhu Han

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04766v1

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

摘要: This paper is devoted to the study of normalized solutions to the Kirchhoff type equation with a logarithmic perturbation

under the normalized constraint $\int_{\mathbb{R}^3} u^2 \,\mathrm{d}x = c^2$, where $a,b>0$, $2<p\leq 6$, $c>0$ is a constant, and $λ\in\mathbb{R}$ emerges as a Lagrange multiplier which is not a priori known. A unified variational framework is developed based on Orlicz spaces together with the Pohozaev constraint method and refined fiber map analysis. For $2<p<\frac{14}{3}$ or $p=\frac{14}{3}$ with small mass, the energy functional is bounded from below and admits a positive radial ground state minimizer. For $\frac{14}{3}<p<6$, where the energy functional is unbounded from below, we establish the existence of two normalized solutions for small mass: a ground state $u_c^+$ obtained via local minimization, and a second solution $u_c^-$ obtained via minimization on the negative component of the Pohozaev manifold. For the Sobolev critical case $p=6$, we construct a ground state solution and, under a technical condition on the parameters, a second solution by introducing a proper auxiliary functional and precise energy estimates with Aubin-Talenti bubbles. Asymptotically as $c\to0^+$, the $L^{2}$ norm of the gradient of ground state solution vanishes for $2<p\le6$. Surprisingly, for $\frac{14}{3}<p<6$, the $L^{2}$ norm of the gradient of the second solution diverges to infinity as $c\to 0^+$, while for $p=6$ it concentrates around the Aubin-Talenti bubble with energy converging to the energy level of the corresponding critical Kirchhoff equation.

9. 广义哈密顿梯度流的近期进展：奇点

Recent progress in generalized Hamiltonian gradient flow: Singularities

作者: Wei Cheng, Jiahui Hong

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04658v1

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

摘要: This paper is a survey of the generalized Hamiltonian gradient flow (GHGF) framework for Hamilton-Jacobi equations, with an emphasis on the propagation of singularities and its connections to weak KAM theory, optimal transport and mean field control. In addition to reviewing the main ideas and known results, we present two new contributions. First, we provide a variational construction of generalized characteristics via a minimizing movement scheme; by taking the weak limit of approximate solutions and using Young measure compactness, we show that the limiting curve satisfies the generalized characteristic differential inclusion. Second, we lift the forward dynamics of strict singular characteristics to a semi-flow on the space of trajectories and study its invariant probability measures. We prove that the only invariant measures of the GHGF semi-flow that attain the critical value $c[H]$ are precisely the projected Mather measures, thereby giving a new dynamical characterization of Mather’s minimal measures as well as Mañé’s critical value. Finally, we discuss a number of open problems that arise from the GHGF perspective, including questions on uniqueness of strict singular characteristics, rectifiability of cut loci, stability under perturbations, contact Hamiltonian systems, vanishing noise limits, and extensions to non-convex or low-regularity Hamiltonians. These problems highlight the deeper connections between singular dynamics, ergodic theory, optimal transport, and geometric analysis, and indicate directions for future research.

10. 通过加权 $L^2$ 上同调研究梯度 Ricci 收缩子的拓扑

Topology of gradient Ricci shrinkers via weighted $L^2$ cohomology

作者: Fei He

类别: math.DG, math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04476v1

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

摘要: This paper proves several topological results for smooth gradient Ricci shrinkers. We establish upper bounds for the Betti numbers, a vanishing theorem for cohomology, and a dichotomy for the number of ends. We also prove a full Hodge theorem for a large class of shrinkers. The methods are based on weighted $L^2$ cohomology and extend to self-shrinkers of the mean curvature flow.

11. 排斥/吸引趋化模型带Logistic型源项的行波解的存在性、唯一性、稳定性与单调性

Existence, uniqueness, stability, and monotonicity of traveling waves for repulsion/attraction chemotaxis models with logistic type source

作者: Wenxian Shen

类别: math.AP

发布日期: 2026-05-06

arXiv ID: 2605.04401v1

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

摘要: This paper is devoted to the study of existence, uniqueness, stability, and monotonicity of traveling wave solutions to the following parabolic-elliptic chemotaxis system with logistic type source \begin{equation}\label{E:main-abstract-eq}\tag{CM} \begin{cases} u_t=u_{xx}-χ(u^m v_x)x +u(1-u^α),\quad &x\in\mathbb{R}\cr 0=v{xx}-v+u^γ,\quad&x\in\mathbb{R}, \end{cases} \end{equation} where $m,α,γ\ge 1$ and $χ\in\mathbb{R}$. System (CM) can be used to describe the evolution of a biological species influenced by a chemical substance produced by the species itself. In this context, the function $u$ denotes the population density of the biological species, and $v$ denotes the concentration of the chemical agent. Traveling wave solutions of (CM) connecting the two constant solutions $(1,1)$ and $(0,0)$ are among important types of solutions, which characterize the front propagation phenomena in (CM). The existence of such traveling wave solutions to (CM) with $m=α=γ=1$ has been studied in several papers. However, there is little study on the uniqueness, stability, and monotonicity of traveling wave solutions of (CM) in literature and there is also no study on the existence of traveling wave solutions of (CM) connecting $(1,1)$ and $(0,0)$ for general $m,α,γ\ge 1$. In the current paper, we prove the existence of traveling wave solutions of (CM) connecting $(1,1)$ and $(0,0)$ for any $χ\le 0$ with speed $c$ large than some number $c^*_{χ,m,γ}$, or for $0<χ<1/2$ with any speed $c>2$. We prove that the traveling wave solutions are monotone when $χ\le 0$. We also prove the uniqueness and stability of traveling wave solutions of (CM) connecting $(1,1)$ and $(0,0)$ when the speed $c$ is larger than some number $c^{**}_{χ,m,α, γ}(\ge c^*_{χ,m,γ})$.

12. 三角形上的Sharp Dirichlet特征值不等式

Sharp Dirichlet eigenvalue inequalities on triangles

作者: Ryoki Endo, Xuefeng Liu, Phanuel Mariano

类别: math.SP, math.AP

发布日期: 2026-05-05

arXiv ID: 2605.04331v1

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

摘要: We prove sharp Dirichlet eigenvalue inequalities for planar triangles. We settle a conjecture of Laugesen and Siudeja by showing that the equilateral triangle uniquely minimizes a scale-invariant functional of the first Dirichlet eigenvalue, area, and perimeter. Consequences include an optimal two-term lower bound for the first Dirichlet eigenvalue in terms of area and perimeter. We also prove a Cheeger-type inequality with an explicit best constant considered by Parini. To prove these conjectures we propose a new method for proving Dirichlet eigenvalue inequalities on triangles. Our method is based on a new computable lower bound for second-order directional shape derivatives under vertex perturbations. It also uses validated finite-element error estimates and recently developed analytic estimates for eigenvalues of nearly degenerate triangles. The method is not specific to the functionals considered in this paper and it can be used to prove various other eigenvalue inequalities on triangles.