生成时间: 2026年10月01日 09:40:29

运行状态: 完成

计划处理论文数量: 41

已完成论文数量: 41

未完成论文数量: 0

重点关注论文数量: 7

了解领域论文数量: 11

不相关论文数量: 23

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

1. 具有自由表面的可压缩液体的哈达玛局部适定性

Hadamard Local Well-Posedness for Compressible Liquids with a Free Surface

作者: Tao Luo, Ruixi Zhang

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.40163v1

摘要翻译: 我们证明了描述具有自由表面且无表面张力液体的三维可压缩 Euler 方程在哈达玛意义下的局部适定性。在泰勒符号条件和边界相容条件下，对于每个 $$s>3$$ ，我们在一个混合 Sobolev 类中建立了存在性、唯一性和强连续依赖性。速度和界面具有 $$H^s$$ 正则性，而焓和速度散度分别属于 $H^{s+1/2}$ 和 $H^{s-1/2}$ 。这种正则性划分反映了非退化声学方程与自由表面和涡量动力学之间的耦合。证明使用正则化算子对移动区域进行平滑，并扩展非线性相容性层级；同时，对定义在不同区域上的解建立 $$L^2$$ 和部分 $$H^2$$ 距离估计。通过正则化前向 Euler 格式构造正则解。随后，高阶能量估计和频率包络给出光滑近似解的强收敛，并证明数据到解映射的连续性。整数 Sobolev 指数则通过 Lions–Magenes 端点相容条件处理。

Comment: 71 pages

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

1. 研究对象和背景

本文研究三维可压缩 Euler 液体自由边界问题, 无表面张力. 流体占据随时间变化的有界区域 $\Omega_t\subset\mathbb{R}^3$ , 边界为 $\Gamma_t=\partial\Omega_t$ , 速度为 $$u$$ , 密度为 $\rho$ . 原始模型是

\begin{cases} (\partial_t+u\cdot\nabla)\rho+\rho\,\operatorname{div}u=0, & x\in\Omega_t,\\ \rho(\partial_t+u\cdot\nabla)u+\nabla P(\rho)=0, & x\in\Omega_t. \end{cases}

自由面随流体运动, 即物质导数 $D_t=\partial_t+u\cdot\nabla$ 在时空自由边界上切向. 边界压强取常数, 并将边界密度归一化为 $$1$$ . 对文中采用的线性压强律 $P(\rho)=\rho$ , 令 $h=\log\rho$ , 方程化为

\begin{cases} D_t h+\operatorname{div}u=0, & x\in\Omega_t,\\ D_t u+\nabla h=0, & x\in\Omega_t,\\ h=0, & x\in\Gamma_t. \end{cases}

因此边界密度保持正值, 声速在边界不退化. 这与 physical vacuum 情形不同: 后者密度和声速在边界消失, 需要适配退化算子的加权正则性; 本文处理的是非退化液体, 但 Dirichlet 边界条件 $h|_{\Gamma_t}=0$ 及其物质时间导数会产生一串非线性相容条件.

核心未知量不仅包括体内的 $$h,u$$ , 还包括自由面几何 $\Gamma_t$ . 边界由速度的法向分量推进, 而边界稳定性受 Taylor 系数控制:

a=-\nabla h\cdot n_{\Gamma_t},\qquad a|_{t=0}\ge c>0.

该符号条件使线性化能量中的边界项具有正性; 若其失效, 会出现 Rayleigh–Taylor 型不稳定. 因为 $D_t h=-\operatorname{div}u$ 且 $$h$$ 在边界为零, 正则解还必须满足 $\operatorname{div}u|_{\Gamma_t}=0$ , 继续作物质微分则得到更高阶条件. 记

z_j(h_0,u_0)=(-D_t)^j h\big|_{t=0},

其中各时间导数借助方程改写成初始数据的空间微分与非线性组合, 相容性要求相应的 $$z_j$$ 在边界消失.

该问题的重要分析结构是可压缩声学演化与自由面、涡量动力学同时存在. 对第一条方程再施加 $$D_t$$ , 可得到

D_t^2h-\Delta h=\operatorname{tr}\big((\nabla u)^2\big),

所以 $$h$$ 与 $\operatorname{div}u=-D_th$ 构成带移动 Dirichlet 边界的非退化声学部分; 自由面法向、切向速度与涡量则具有类似不可压缩自由边界问题的几何演化. 两部分通过 Euler 方程相互耦合, 不能像不可压缩情形那样仅从瞬时速度解一个椭圆压力方程来消去声学变量. 尤其 Taylor 系数的演化含有法向导数 $\nabla_n\operatorname{div}u$ , 这使得控制自由面能量需要额外的声学正则性.

论文采用的混合 Sobolev 标度反映了这种不同的正则性需求: 自由面和速度处于 $$H^s$$ , 而 $$h$$ 比 $$u$$ 多半阶, $\operatorname{div}u$ 比 $$u$$ 少半阶, 即 $h\in H^{s+1/2}$ 、 $\operatorname{div}u\in H^{s-1/2}$ . 该结构也对应于声学部分与自由面部分各自的尺度, 而非整个耦合系统的单一共同尺度. 问题的目标不止是构造局部解, 而是在自然的移动区域拓扑中证明存在性、唯一性以及数据到解映射的强连续性, 即 Hadamard 意义下的局部适定性. 既有压缩液体自由面理论已发展出局部存在与先验估计, 但如何在低正则混合空间中同时处理相容性保持的逼近、不同移动区域上的解差以及强连续依赖, 是本文试图弥补的缺口。对一般光滑压强律, 作者说明可用归一化焓变量及变系数能量作相应推广; 论文详细证明则限于线性压强律。

2. 主要定理或主要结果

主要结论是定理 1.4. 令参考边界 $\Gamma_*$ 为光滑有界曲面, 初始自由面是其附近的法向图, 且初值 $(h_0,u_0,\Omega_0)$ 属于文中定义的混合空间 $H^s_{1/2}$ , $$s>3$$ . 假设初始 Taylor 条件成立, 并满足有意义的边界相容性条件

z_j(h_0,u_0)\big|_{\Gamma_0}=0,\qquad 0\le j<s.

若 $s\in\mathbb{N}$ , 还须满足端点条件 $z_s\in H^{1/2}_{00}(\Omega_0)$ ; 此时 $H^{1/2}_{00}$ 表示零延拓属于 $H^{1/2}(\mathbb{R}^3)$ 的 Lions–Magenes 空间, 不要求对只有 $H^{1/2}$ 正则性的 $$z_s$$ 取经典边界迹. 则存在 $$T>0$$ , 其取值仅依赖于初始混合范数及 $\inf_{\Gamma_0}a_0$ , 使得存在唯一解

(h,u,\Omega_t)\in C([0,T],H^s_{1/2}).

这里自由面 $\Gamma_t$ 属于 $$H^s$$ , $u\in H^s(\Omega_t)$ , 而声学变量满足

h\in H^{s+1/2}(\Omega_t),\qquad \operatorname{div}u\in H^{s-1/2}(\Omega_t).

论文不要求无旋性, 涡量由速度的 $$H^s$$ 控制, 因而处于 $H^{s-1}$ 正则性. 结论针对任意光滑参考边界, 并非仅限于平底或图形自由面.

定理还给出强连续依赖. 若一列初始状态在 $H^s_{1/2}$ 中收敛到给定初值, 则在共同存在时间区间上, 对应解在 $C([0,T],H^s_{1/2})$ 中收敛. 因区域随数据变化, 此处收敛按文中定义理解: 将自由面写为 $\Gamma_*$ 上的法向图, 要求图函数在 $H^s(\Gamma_*)$ 收敛, 并选取到 $\mathbb{R}^3$ 的延拓, 使 $$u$$ 在 $$H^s$$ 中以及 $(h,\operatorname{div}u)$ 在 $H^{s+1/2}\times H^{s-1/2}$ 中收敛. 当 $$s$$ 为整数时, 还要求 $z_s\mathbf{1}_{\Omega}$ 的零延拓在 $H^{1/2}(\mathbb{R}^3)$ 中收敛. 这比仅在较弱拓扑下得到紧性极限更强.

证明中使用的另一组结论服务于上述定理: 作者构造相容性保持的正则化算子 $\Psi_m$ , 将粗糙初值平滑为更高正则状态, 同时满足预先指定的有限阶相容条件. 命题 2.2 给出相应正则性增长与逼近误差界: 平滑后 $H^{s+r}$ 型范数至多按 $2^{mr}$ 增长乘以原有控制量, 而低阶状态误差为 $O(2^{-ms})$ ; 对 $(\operatorname{div}u,\nabla h)$ 的误差则为 $O(2^{-m(s'-1)})$ , 其中 $$s'=s$$ 或 $$s+1/2$$ . 该构造还使 $z_{j,m}|_{\Gamma_m}=0$ 对预设的所有有限阶 $j\le K$ 成立.

作者进一步建立不同自由面解之间的 $$L^2$$ 距离估计及部分 $$H^2$$ 距离估计, 并在高阶能量 $$H^k$$ 、 $$H^k_1$$ 与频率包络控制的基础上, 通过正则化前向 Euler 离散方案构造正则解. 这些是证明存在、唯一以及混合拓扑强收敛的关键中间结论, 不是对定理适用范围的额外放宽。作者明确指出 $$s>3$$ 的阈值未证明为最优; 从给出的结果不能推出该阈值可进一步降低。

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

证明的首要障碍是声学变量与自由面几何的闭合耦合. 对线性化扰动 $$(z,w)$$ 和边界法向位移 $\vartheta$ , 线性化边界条件为 $z=a\vartheta$ , 相应能量包含体积分和带 Taylor 权的边界项:

E_{\mathrm{lin}}=\int_{\Omega_t}(z^2+|w|^2)\,dx+\int_{\Gamma_t}a\vartheta^2\,dS.

命题 2.1 对带源项的线性化系统导出能量恒等式. 能量导数中出现 $$D_ta$$ , 而其主项含 $\nabla_n\operatorname{div}u$ ; 仅有 $$u$$ 的 Lipschitz 控制不足以控制该项. 因此作者不是给所有未知量强加同一 Sobolev 阶数, 而是利用 $\operatorname{div}u\in H^{s-1/2}$ (在三维且 $$s>3$$ 时可提供所需点态控制) 及 $h\in H^{s+1/2}$ , 配合速度、自由面的 $$H^s$$ 估计. 高阶估计把物质导数与边界切向导数结合: 物质导数捕捉声学演化, 切向能量与曲率演化相配合以消去边界项, 随后借助 div–curl 和椭圆估计恢复完整空间正则性. 作者还传播一个额外声学导数层级 $$H^k_1$$ , 而不要求速度和自由面整体多一阶正则性; 这两层能量共同恢复最终的半阶正则性分离。

第二个核心困难是逼近不能破坏非线性边界相容性. 低阶条件包括 $h|_\Gamma=0$ , $\operatorname{div}u|_\Gamma=0$ 以及

\left(\Delta h+\operatorname{tr}((\nabla u)^2)\right)\big|_\Gamma=0,

更高阶条件涉及更高导数与非线性项. 普通卷积既改变边界, 也不保留这些迹恒等式; 同时粗糙初值只满足其正则性所允许的有限条件, 但构造正则演化又需要更高阶条件. 作者以分解 $$z_j=Z_j+N_j(h,u)$$ 将主线性边界算子与非线性余项分开, 再通过椭圆边界修正逐阶消除误差, 构造 $\Psi_m$ . 该算子同时平滑区域和流体变量、保持已有相容性并补足指定的高阶相容性, 还给出依赖频率尺度 $$2^m$$ 的正则性增长和逼近误差. 整数 $$s$$ 的端点上, $$z_s$$ 只有 $H^{1/2}$ 正则性, 未必有经典迹; 采用 $H^{1/2}_{00}$ 条件使零延拓具备所需边界控制, 从而纳入同一正则化与收敛框架。

第三个困难在于离散时间推进也必须定量控制相容性. 直接前向 Euler 步会损失空间导数并破坏高阶边界条件; 每一步重新平滑虽可恢复正则性, 却会扰动边界迹和能量. 固定时间内步数约为 $\varepsilon^{-1}$ , 因而单步能量误差必须为 $O(\varepsilon)$ , 不能只要求每步能量有界. 作者将区域、声学部分及切向速度分阶段用椭圆 resolvent 正则化, 再用具有指定法向导数的边界提升算子削弱相容性误差, 并利用 resolvent 的耗散吸收最高阶误差. 定理 5.1 传播的是量化的近似相容性, 而非要求每个中间离散状态严格满足全部相容条件. 文中给出 $O(\varepsilon)$ 单步能量增量与 $O(\varepsilon^2)$ 一致性误差, 由此获得对步长一致的控制并通过极限构造正则解。

第四个障碍是比较不同移动区域上的解. 作者沿用 Eulerian 距离思路, 在解区域交集上建立 $$L^2$$ 距离, 并依据线性化边界能量给声学差配置 Taylor 系数倒数权重. Taylor 符号使焓差在自由面附近与两边界的法向间距联系起来, 因而该距离同时控制体变量与界面位移; 其传播用于唯一性及稳定性. 但 $$L^2$$ 距离不足以得到混合空间要求的声学逼近速率, 所以另建基于二阶物质导数的部分 $$H^2$$ 距离, 控制 $\Delta h$ 与 $\nabla\operatorname{div}u$ 的差及低阶余项。两解各自的焓在各自边界为零, 其差在交集边界上却一般不满足齐次 Dirichlet 条件, 因而不能直接把 Laplacian 估计转为完整 Hessian 估计. 作者先用微分同胚对齐区域, 估计椭圆算子因几何变换产生的误差, 再将声学估计转回公共区域。

最后, 紧性或弱收敛本身不能给出 $H^s_{1/2}$ 中的强收敛. 作者将上述距离估计与 $$H^k,H^k_1$$ 高阶能量结合, 对自由面、速度、声学变量分别使用频率包络控制高频尾部. 对收敛初值, 高频尾部的一致控制与固定正则化尺度下的连续性相结合, 得到光滑近似在完整混合拓扑中的强收敛, 并进一步证明数据到解映射连续; 整数端点的零延拓相容性量也包含在这一收敛论证中。

4. 与之前工作的比较

压缩液体自由面局部存在理论方面, 作者将 Lindblad [20,21] 的线性化估计及 Nash–Moser 非线性构造、Trakhinin [33] 的 Euler 坐标 Nash–Moser 理论, 以及 Wang、Zhang、Zhao [38] 关于 Taylor 条件与自由面演化方程双曲性的工作作为背景. 这些工作奠定了正边界密度液体模型的局部存在理论; 本文着重推进的是 Hadamard 意义下的完整低正则适定性论证, 特别是混合 Sobolev 空间中的强连续依赖、相容性保持的正则化和跨移动区域距离估计. 从摘录内容无法确认上述每篇先前工作的精确正则指数与连续依赖拓扑, 因而不宜将本文的结论表述为对所有既有存在结果的统一低正则改进。

最直接的对照是不可压缩自由边界 Euler 理论. Ifrim、Pineau、Tataru、Taylor [14] 在三维得到速度和自由面 $$H^s$$ , $$s>5/2$$ 的 Hadamard 适定性, 其框架使用几何能量、不同区域上的解距离及正则解逼近. Wang、Zhang、Zhao、Zheng [35] 在有限深度、平底及图形自由面等设定下也达到 $$r>5/2$$ , 并以 paradifferential 方法、边界速度 good unknown、Dirichlet–Neumann 算子估计和修正能量处理涡量项表面上的半阶损失, 同时得到 Beale–Kato–Majda 型判据. 相比之下, 本文不要求无旋性, 允许一般光滑参考边界, 但需要同时控制非退化声学演化及其高阶非线性边界相容条件, 得到的阈值为 $$s>3$$ , 且声学量采用 $h\in H^{s+1/2}$ 、 $\operatorname{div}u\in H^{s-1/2}$ 的混合标度. 作者明确表示 $$s>3$$ 未知是否最优, 并指出若要达到不可压缩理论的 $H^{5/2+}$ 阈值, 还需要更精细的声学距离估计以及能在更低正则性下保持相容性的正则化程序。因此本文的代价是更高的正则性门槛和额外的声学相容性结构, 而非把不可压缩阈值直接移植到可压缩液体。

在压缩液体模型自身的工作中, Lindblad 和 Luo [22] 建立允许非零涡量的先验能量估计, 并利用对声速一致的估计研究不可压缩极限; Luo [24] 处理无界区域中的可压缩重力水波先验估计和不可压缩极限, Luo 和 Zhang [25] 则对该水波问题用切向平滑证明局部适定性. 本文的问题是有界液体区域、无表面张力的自由边界问题, 主要目标不是不可压缩极限或重力水波长时间行为, 而是建立一般初值混合拓扑下的 Hadamard 局部理论. 从摘录无法确认 [25] 与本文在全部函数空间及连续依赖结论上的逐项差别, 因此只能确定模型与证明侧重点不同, 不应推断其结论被本文全面包含。

物理真空理论构成另一类重要对照. Coutand–Shkoller [4]、Jang–Masmoudi [17] 等研究密度和声速在边界退化的三维物理真空, Ifrim–Tataru [13] 及 Liu–Luo [23] 发展适配该退化结构的 Hadamard 理论. 本文的边界密度为正, 声学算子到边界仍非退化, 因而使用普通 Sobolev 空间, 而非物理真空所需的加权空间. 但非退化并不意味着边界条件更简单: 移动 Dirichlet 条件强迫 $\operatorname{div}u$ 及其进一步物质导数满足边界约束, 使相容性保持成为这里特有的技术负担. 论文关于一般压力律的讨论也明确限于密度处于正值紧区间且平方声速有正上下界的液体情形, 不适用于 physical vacuum。

其他相邻模型显示本文方法的范围边界. Coutand、Hole、Shkoller [3] 研究有表面张力情形并论证零表面张力极限; Huang、Li、Wang [12] 处理二维有界角域中的可压缩毛细–重力水波, 还需分析角点奇性、接触角与耗散. Ginsberg、Lindblad、Luo [7] 的自引力液体, Luo、Trivisa、Zeng [26] 的非等熵 Euler 和 Euler–Poisson, 以及磁流体自由边界工作 [34] 均含有与本文不同的额外场或稳定性机制. 本文不含表面张力、重力、角点或额外热力学变量, 也不声称建立长时间理论; 其结果是 Taylor 符号假设下的局部适定性。

方法论上, 本文从 [14,13,23] 的 Eulerian 低正则框架吸收了跨区域距离、正则解逼近和频率包络等思路, 但必须为压缩液体增加声学部分的部分 $$H^2$$ 距离、相容性修正及混合正则性的高频尾部控制. 作者声称对足够光滑且 $P'(\rho)>0$ 的一般压力律, 可改用归一化焓 $$h_P$$ 与平方声速 $$q(h_P)$$ , 在密度范围紧含于 $(0,\infty)$ 且 $$q$$ 有正上下界时沿用同一 Sobolev 指数; 不过详细证明仅针对线性压强律, 变系数项和相容性表达式的实际推导并未在本文摘录中展开. 因此这一推广是作者说明可按相同估计处理的结论, 不应误作本文已逐项呈现的完整变系数证明。

2. $\mathbb{R}^3$ 中两相流模型的低马赫数极限

Low Mach number limit of a two-phase flow model in $\mathbb{R}^3$

作者: Fucai Li, Jinkai Ni, Zhipeng Zhang, Zhu Zhang

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39900v1

摘要翻译: 我们研究 $\mathbb R^3$ 中可压缩 Navier–Stokes–Euler 两相系统的同时低马赫数极限，其中两相的压力项均按 $\varepsilon^{-2}$ 缩放。对于常值平衡态附近充分小的 $$H^3$$ 扰动，我们建立缩放可压缩系统的整体适定性，并推导出关于马赫数 $\varepsilon$ 一致的整体能量耗散估计。分析的一个关键特征是耗散结构的退化性：黏性仅作用于 Navier–Stokes 相，而阻力耦合通过弛豫模态将耗散传递至无黏 Euler 相。随后，我们证明极限不可压缩两相系统的整体适定性和大时间衰减，并表明相对速度 $$u-ω$$ 比完整速度对衰减更快。最后，对于良备初值，我们引入压力修正的声学变量，以消除奇异压力失配，并建立一个整体时间范围内的 $$H^2$$ 误差估计。因此，缩放可压缩解关于时间一致地收敛到相应的极限不可压缩系统解。

Comment: 32 pages. All comments are welcome

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

1. 研究对象和背景

本文研究的对象是三维全空间 $\mathbb{R}^3$ 上的可压缩 Navier–Stokes–Euler 两相流模型. 一相为具有黏性的可压缩流体, 密度和速度分别记为 $\rho$ 、 $$u$$ ; 另一相为无黏的可压缩 Euler 流体, 密度和速度分别记为 $$n$$ 、 $\omega$ . 两相通过与相对速度 $u-\omega$ 成正比的阻力耦合. 原模型为

\begin{cases} \partial_t\rho+\operatorname{div}(\rho u)=0,\\ \partial_t(\rho u)+\operatorname{div}(\rho u\otimes u)+\nabla P(\rho) =\mu\Delta u+(\mu+\lambda)\nabla\operatorname{div}u-\kappa n(u-\omega),\\ \partial_t n+\operatorname{div}(n\omega)=0,\\ \partial_t(n\omega)+\operatorname{div}(n\omega\otimes\omega)+\nabla n =\kappa n(u-\omega), \end{cases} \qquad P(\rho)=A\rho^\gamma,

其中 $$A>0$$ , $\gamma\geq 1$ , $\mu>0$ , $2\mu+\lambda>0$ , $\kappa>0$ . 这一模型可作为颗粒与流体相互作用的宏观描述, 与喷雾、沉降等问题相关; 文中提及其动力学背景可来自 Vlasov–Fokker–Planck/Navier–Stokes 系统在强噪声及局部对齐条件下的流体极限.

研究的奇异极限是同时令两相的压力都处于低马赫数尺度. 经过时间、速度及系数缩放后, 压力项具有 $\varepsilon^{-2}$ 系数, 其中 $\varepsilon\in(0,1)$ 为马赫数. 以常密度平衡态 $$(1,0,1,0)$$ 为中心, 写成 $\rho^\varepsilon=1+\varepsilon q^\varepsilon$ , $n^\varepsilon=1+\varepsilon r^\varepsilon$ , 两个速度记作 $u^\varepsilon,\omega^\varepsilon$ , 则扰动方程包含 $\varepsilon^{-1}$ 阶声学耦合项. 例如连续性方程的主部为

\partial_tq^\varepsilon+u^\varepsilon\cdot\nabla q^\varepsilon +\frac{1+\varepsilon q^\varepsilon}{\varepsilon}\operatorname{div}u^\varepsilon=0, \qquad \partial_tr^\varepsilon+\omega^\varepsilon\cdot\nabla r^\varepsilon +\frac{1+\varepsilon r^\varepsilon}{\varepsilon}\operatorname{div}\omega^\varepsilon=0.

形式上, 若密度扰动及其导数有界, $\varepsilon\to0$ 时两相速度都趋于无散场. 对动量方程施加 Leray 投影, 压力梯度消失, 得到两相均不可压缩、但各自具有压力的极限系统:

\begin{cases} \partial_tu+u\cdot\nabla u+\nabla\pi_1-\mu\Delta u=\omega-u,\\ \operatorname{div}u=0,\\ \partial_t\omega+\omega\cdot\nabla\omega+\nabla\pi_2=u-\omega,\\ \operatorname{div}\omega=0. \end{cases}

这里 $\pi_1$ 、 $\pi_2$ 分别是维持两相无散约束的压力. 问题的重要分析困难在于耗散是退化的: 黏性仅作用于 $$u$$ , Euler 相的 $\omega$ 没有自身的黏性项; 而低马赫数缩放又引入奇异声学项, 不能简单视为高阶方程中的小扰动. 因此需要同时处理声学振荡、非线性高阶估计以及由阻力耦合传递的间接耗散.

已有理论分别涉及可压缩 Navier–Stokes–Euler 系统的小扰动整体适定性与大时间衰减, 以及单相可压缩流体的低马赫数极限. 但本文所考虑的双重压力缩放使两相都进入不可压缩极限, 与只对 Navier–Stokes 相作低马赫数缩放、而让 Euler 相保持可压缩的设定不同. 本文试图补足的空白是: 在 Euler 相无内禀黏性的退化结构下, 对该双相同时低马赫数问题建立对 $\varepsilon$ 一致的整体解估计, 并在良备初值下给出全时间的定量误差估计, 而不只依赖紧性论证得到定性收敛.

2. 主要定理或主要结果

以下结论均在 $\mathbb{R}^3$ 上成立, 常数可取为与 $\varepsilon$ 无关. 记 $P(\rho)=A\rho^\gamma$ , 并令平衡态密度为 $$1$$ .

定理 1.1 给出缩放可压缩系统的整体适定性及一致能量耗散. 若初值属于 $H^3(\mathbb{R}^3)$ , 满足 $1+\varepsilon q_0^\varepsilon>0$ , $1+\varepsilon r_0^\varepsilon>0$ , 且

X_0^\varepsilon=\|(q_0^\varepsilon,u_0^\varepsilon,r_0^\varepsilon,\omega_0^\varepsilon)\|_{H^3}^2\leq\delta_0^2,

其中 $\delta_0$ 与 $\varepsilon$ 无关且充分小, 则存在唯一整体强解 $(q^\varepsilon,u^\varepsilon,r^\varepsilon,\omega^\varepsilon)\in C([0,\infty);H^3)$ , 且两密度保持正. 并有

\sup_{t\geq0}\|(q^\varepsilon,u^\varepsilon,r^\varepsilon,\omega^\varepsilon)(t)\|_{H^3}^2 +\int_0^\infty\!\left(\|u^\varepsilon-\omega^\varepsilon\|_{H^3}^2 +\|\nabla(q^\varepsilon,r^\varepsilon,\omega^\varepsilon)\|_{H^2}^2 +\|\nabla u^\varepsilon\|_{H^3}^2\right)dt \leq C_0X_0^\varepsilon.

该估计明确体现耗散的层级差异: $u^\varepsilon$ 有完整的 $$H^3$$ 梯度耗散, $\omega^\varepsilon$ 仅得到 $$H^2$$ 阶梯度耗散, 但相对速度在 $$H^3$$ 中可积平方.

定理 1.2 对极限不可压缩系统给出整体适定性. 若 $U_0=(u_0,\omega_0)\in H^3$ , $\operatorname{div}u_0=\operatorname{div}\omega_0=0$ , 且 $X_0=\|U_0\|_{H^3}^2\leq\delta_1^2$ , 则存在唯一整体强解 $U=(u,\omega)\in C([0,\infty);H^3)$ , 满足

\sup_{t\geq0}\|U(t)\|_{H^3}^2 +\int_0^\infty\left(\|u-\omega\|_{H^3}^2+\|\nabla u\|_{H^3}^2+\|\nabla\omega\|_{H^2}^2\right)dt \leq C_1X_0.

若进一步 $U_0\in L^1$ 且 $\|U_0\|_{H^3\cap L^1}\leq\delta_1$ , 则有大时间衰减

\|U(t)\|_{H^3}\leq C(1+t)^{-3/4}\|U_0\|_{H^3\cap L^1},

\|\nabla U(t)\|_{H^2}+\|u(t)-\omega(t)\|_{H^3} \leq C(1+t)^{-5/4}\|U_0\|_{H^3\cap L^1}.

特别地, 相对速度比完整速度对 $$U$$ 多衰减一个 $(1+t)^{-1/2}$ 因子.

定理 1.3 给出良备初值下的整体定量低马赫数极限. 极限压力由 Riesz 变换表示为

\pi_1=\sum_{i,j=1}^3R_iR_j(u_i u_j),\qquad \pi_2=\sum_{i,j=1}^3R_iR_j(\omega_i\omega_j).

要求初值满足

\|u_0^\varepsilon-u_0\|_{H^2} +\|q_0^\varepsilon-\varepsilon[P'(1)]^{-1}\pi_1(0)\|_{H^2} +\|\omega_0^\varepsilon-\omega_0\|_{H^2} +\|r_0^\varepsilon-\varepsilon\pi_2(0)\|_{H^2}\leq\varepsilon.

则压力修正后的误差满足

\begin{aligned} &\sup_{t\geq0}\left(\|u^\varepsilon-u\|_{H^2}^2 +\|q^\varepsilon-\varepsilon[P'(1)]^{-1}\pi_1\|_{H^2}^2 +\|\omega^\varepsilon-\omega\|_{H^2}^2 +\|r^\varepsilon-\varepsilon\pi_2\|_{H^2}^2\right)\\ &\quad+\int_0^\infty\left(\|\nabla(u^\varepsilon-u)\|_{H^2}^2 +\|\nabla(\omega^\varepsilon-\omega)\|_{H^1}^2 +\|\nabla(q^\varepsilon-\varepsilon[P'(1)]^{-1}\pi_1)\|_{H^1}^2\right.\\ &\qquad\qquad\left.+\|\nabla(r^\varepsilon-\varepsilon\pi_2)\|_{H^1}^2 +\|(u^\varepsilon-\omega^\varepsilon)-(u-\omega)\|_{H^2}^2\right)dt \leq C_2\varepsilon^2. \end{aligned}

因此速度误差在 $L^\infty_tH^2_x$ 中为 $O(\varepsilon)$ ; 同时 $q^\varepsilon,r^\varepsilon$ 的主阶不是零, 而分别由 $\varepsilon[P'(1)]^{-1}\pi_1$ 、 $\varepsilon\pi_2$ 给出. 对原密度而言, $\rho^\varepsilon-1=\varepsilon q^\varepsilon$ 、 $n^\varepsilon-1=\varepsilon r^\varepsilon$ , 所以密度的首个非零压力修正是 $O(\varepsilon^2)$ . 文中进一步给出 $(q^\varepsilon,u^\varepsilon,r^\varepsilon,\omega^\varepsilon)\to(0,u,0,\omega)$ 于 $C([0,\infty);H^2)$ , 以及 $\varepsilon^{-2}(\rho^\varepsilon-1)\to[P'(1)]^{-1}\pi_1$ 、 $\varepsilon^{-2}(n^\varepsilon-1)\to\pi_2$ ; 两个速度的非无散部分也在同一时空连续范数中趋于零.

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

证明的第一重障碍是扰动系统中 $\varepsilon^{-1}$ 阶声学项. 直接把这些项当作高阶估计中的源项会产生随 $\varepsilon\to0$ 发散的常数. 作者改用对称化能量: 在密度方程和动量方程配对时, 选择依赖密度的权重, 使声学主部在积分分部后相消. 高阶能量中相应出现类似

\left\|\sqrt{\frac{P'(1+\varepsilon q^\varepsilon)}{1+\varepsilon q^\varepsilon}}\,\partial^\alpha q^\varepsilon\right\|_2^2, \qquad \left\|\frac{\partial^\alpha r^\varepsilon}{\sqrt{1+\varepsilon r^\varepsilon}}\right\|_2^2,

的加权密度项, 与速度的 $$L^2$$ 能量组成对称双曲型主能量. 关键不是单纯使用能量法, 而是让权重与非线性压力系数相匹配, 从而避免显式控制大小为 $\varepsilon^{-1}$ 的项.

第二个关键是高阶交换子中奇异系数的抵消. 对 $1\leq|\alpha|\leq3$ 微分后, 连续性方程中的系数可利用

\frac{1}{\varepsilon}[\partial^\alpha,(1+\varepsilon q^\varepsilon)\operatorname{div}]u^\varepsilon =[\partial^\alpha,q^\varepsilon\operatorname{div}]u^\varepsilon

以及对应的 $r^\varepsilon,\omega^\varepsilon$ 恒等式改写. 右侧不再带 $\varepsilon^{-1}$ . 对压力系数的交换子, 密度系数的每次微分会带出 $\varepsilon$ , 这抵消了表面上的奇异因子. 再以 Sobolev 乘积和交换子估计控制非线性项, 得到与马赫数无关的 $$H^3$$ 能量界. 小扰动假设同时保证 $1+\varepsilon q^\varepsilon$ 、 $1+\varepsilon r^\varepsilon$ 一致远离零, 使这些加权能量等价于标准 Sobolev 范数.

第三个障碍是耗散退化. 黏性只控制 $u^\varepsilon$ 的梯度, 不能直接吸收 Euler 相的高阶非线性. 阻力项给出相对速度耗散 $\|u^\varepsilon-\omega^\varepsilon\|^2$ , 并与 $u^\varepsilon$ 的黏性耗散结合, 通过耦合或交互型估计把耗散传递给 $\omega^\varepsilon$ 及密度梯度. 所得耗散并非两相完全对称: 对 $\omega^\varepsilon$ 只恢复较低阶的空间导数, 但作者证明这一较弱耗散足以控制非线性项并闭合整体 $$H^3$$ 估计. 这一结构也解释了定理 1.1 中 $\nabla\omega^\varepsilon$ 只在 $$H^2$$ 阶可积平方.

对不可压缩极限系统, 作者结合线性化系统的频谱分析和带时间权的非线性能量方法. 低频部分决定三维扩散型的代数衰减, 而阻力耦合的松弛模态为相对速度 $u-\omega$ 提供额外衰减因子, 从而得到 $$H^3$$ 速度对的 $(1+t)^{-3/4}$ 衰减与相对速度的 $(1+t)^{-5/4}$ 衰减. 同时, 对两条动量方程取散度得到压力的椭圆方程; 由此建立的压力梯度及压力物质导数估计, 是处理后续压力修正误差时不可缺少的输入. 摘录未给出完整频谱分解公式及相关命题的具体证明细节, 因而不能进一步确认各频段符号的精确表达.

最后, 低马赫数误差不能直接用 $q^\varepsilon,r^\varepsilon$ 与零比较: 它们包含由不可压缩压力决定的 $O(\varepsilon)$ 平衡剖面, 若不先扣除该项, 误差方程仍会保留奇异压力失配. 作者定义

\widetilde q=q^\varepsilon-\varepsilon[P'(1)]^{-1}\pi_1,\quad \widetilde r=r^\varepsilon-\varepsilon\pi_2,\quad \widetilde u=u^\varepsilon-u,\quad \widetilde\omega=\omega^\varepsilon-\omega.

压力修正把可压缩密度扰动中的准平衡声学部分提取出来, 使误差系统的主压力项能够与极限方程的压力相消. 随后沿用对称化能量、阻力松弛耗散及交互估计, 并将极限压力的梯度和物质导数界用于控制修正项. 文中得到的误差泛函满足示意不等式

\widetilde X(t)\lesssim\widetilde X(0)+\varepsilon^2+(\delta_0+\delta_1)\widetilde X(t).

结合良备初值给出的初始误差以及足够小的扰动, 可吸收最后一项并推出全时间 $$H^2$$ 误差平方为 $O(\varepsilon^2)$ . 因而该论证的关键组织方式是先构造正确的压力修正变量, 再在退化耗散框架中对误差闭合, 而非仅由一致有界性抽取子列.

4. 与之前工作的比较

本文在既有两相模型整体适定性理论与单相低马赫数理论的交界处推进. 对模型 (1.1) 的 Cauchy 问题, 文中回顾 Choi [7] 在全空间及周期区域建立小光滑解整体存在, 并在周期情形证明速度对齐; 该工作揭示黏性相与 Euler 相的耦合可抑制小扰动 Euler 动力学可能出现的奇性. Wu、Zhang、Zou [34] 以及 Tang、Zhang [30] 研究大时间行为并在 $$L^1$$ 型假设下得到最优代数衰减. Zhang 等 [37] 在 $\mathbb{R}^3$ 中以 $$H^3$$ 小性证明整体存在, 不要求更高阶 Sobolev 范数小. Li、Shou [19] 则在多维临界 Besov 框架中研究该模型. 本文与这些工作相衔接, 但核心目标不是单独改善未缩放方程的适定性, 而是获得对低马赫数参数一致的估计并定量识别不可压缩极限.

与 Hong、Jong [14] 的低马赫数研究相比, 本文最明确的模型差别是压力缩放范围. [14] 只在 Navier–Stokes 相引入低马赫数尺度, Euler 相仍保持可压缩; 本文则同时缩放两相压力, 因而极限中同时出现 $\operatorname{div}u=0$ 和 $\operatorname{div}\omega=0$ , 并需要两个压力 $\pi_1,\pi_2$ . 这不是把单相结论直接并列应用: Euler 相没有黏性, 新增的奇异声学结构必须依靠与黏性相的阻力耦合恢复耗散. 摘录没有给出 [14] 的完整定理假设和收敛范数, 因此不能据此作更细的定量比较.

与经典可压缩 Euler 低马赫数理论 (Ebin [11]、Klainerman–Majda [17,18]、Ukai [31]、Schochet [27,28]、Métivier–Schochet [24]) 以及可压缩 Navier–Stokes 低马赫数理论 (Lions–Masmoudi [20]、Danchin [9]、Alazard [1]) 相比, 本文的结构重点是两相耦合与部分黏性并存. 经典单相结果不能直接处理“一个速度有黏性, 另一个速度无黏性, 两者只通过松弛项联系”的耗散退化. 本文在小 $$H^3$$ 扰动范围内先建立统一于 $\varepsilon$ 的整体强解估计, 再将该机制用于两相不可压缩极限, 从而给出特定框架下的全局定量收敛.

结论形式上也超出仅由紧性得到的定性极限. 定理 1.3 在良备数据条件下控制整个时间半轴上的 $$H^2$$ 误差平方为 $C\varepsilon^2$ , 即速度误差为 $O(\varepsilon)$ ; 并且密度扰动经过压力校正后同样有 $$H^2$$ 控制. 这同时给出密度首阶剖面 $\varepsilon^{-2}(\rho^\varepsilon-1)\to[P'(1)]^{-1}\pi_1$ 、 $\varepsilon^{-2}(n^\varepsilon-1)\to\pi_2$ , 以及速度非无散部分趋零. 这里的全时间估计依赖小数据整体耗散, 与只在有限时间区间内进行的误差控制有所不同; 但摘录未逐一列举或比较先前两相极限工作的具体时间区间和收敛率, 故不宜断言它在所有相关文献中首次达到某个普遍最优尺度.

本文的代价和适用边界同样清楚. 结论针对 $\mathbb{R}^3$ 、常密度平衡态附近的小 $$H^3$$ 扰动; 衰减结论另需 $$L^1$$ 初值, 定量极限则需密度初值匹配极限压力, 且速度初值与极限速度相差不超过 $O(\varepsilon)$ 的 $$H^2$$ 量级. 因此定理 1.3 属于良备初值结果, 不能据此推出对一般未准备声学初值也有同样的全时间强收敛或相同误差率. 文中没有在所给摘录中证明大数据适定性、边界区域结果、临界空间中的低马赫数极限, 或允许两相压力采用不同马赫数尺度的统一定理. 这些均不应从现有结论外推. 此外, 文中把新结果与既有研究的主要差异概括为双相同时不可压缩化、退化耗散处理及定量全局误差; 对部分相关文献的详细假设和最佳性比较, 从提供内容中无法确认.

3. Boussinesq 系统存在时间的改进

Improved existence time for a Boussinesq system

作者: Achenef Tesfahun, Sigmund Selberg

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39891v1

摘要翻译: 我们考虑一类强色散 Boussinesq 型系统，该系统是弱非线性水波传播的一个模型。该系统由非线性参数 $0<ε\le 1$ 和浅水参数 $0<μ\le 1$ 表征。在 Boussinesq 区间 $μ\sim ε$ 内，Saut、Wang 和 Xu \cite{SWX2017} 在非空化条件下采用双曲型方法，证明了该系统在一维和二维中的适定性，其时间尺度为 $$1/ε$$ 。本文假设带有参数 $$h_0>0$$ 的非空化条件，建立了一维情形下时间尺度为 $( μ^{ 1/4} ε^{-1} h_0 )^{4/3}$ 、二维情形下时间尺度为 $(μ^{ 1/4}ε^{-1} h_0)^{2-}$ 的适定性结果。我们的结果定量地揭示了在 $ε\lesssim μ$ 和 $μ\llε$ 区间内，存在时间如何取决于非线性与色散之间的平衡。结果还明确说明了当 $h_0\to0$ 时，即波趋近空化时，存在时间如何退化。我们的证明结合了色散估计、Strichartz 估计与能量估计。

Comment: 24 pages

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

1. 研究对象和背景

论文研究一类具有强色散效应的 Boussinesq 型水波方程. 未知量 η(x,t) 表示自由液面的高程, v(x,t) 是与流体速度相关的速度场; 空间维数 d 为 1 或 2. 方程为

\begin{cases} \eta_t+\nabla\cdot v-\mu\nabla\cdot\Delta v+\varepsilon\nabla\cdot(\eta v)=0,\\ v_t+\nabla\eta+\dfrac{\varepsilon}{2}\nabla |v|^2=0. \end{cases}

这里 $0<\varepsilon,\mu\leq 1$ 分别衡量非线性和浅水色散效应. 对二维问题还要求初始速度场无旋, 即 $\operatorname{curl}v_0=0$ ; 这一条件由第二个方程沿演化保持. 该系统是经典 $$(abcd)$$ Boussinesq 系统的一种特例, 对应 $$a=-1,b=c=d=0$$ , 并可作为具有强表面张力的毛细-重力波模型. 其线性部分具有高频 Schrödinger 型色散, 但方程仍是拟线性系统, 非线性项会使能量增长并影响解的存在时间.

适定性问题不仅关乎解的存在、唯一性和对初值的连续依赖, 也关系到该类近似模型能否在水波渐近分析所需的长时间区间内有效. 对 Boussinesq 模型进行 Euler 方程的渐近近似时, 通常希望在 $O(1/\varepsilon)$ 时间内获得适当 Sobolev 范数的统一控制; 与模型一致性结合后, 可导出相对于 Euler 系统的 $O(\varepsilon^2t)$ 误差估计. 然而, Boussinesq 系统一般缺少足以直接给出所需统一控制的守恒律, 长时间估计因而需要利用具体方程结构.

本文特别关注非空化条件. 它要求水深因子满足

\inf_{x\in\mathbb R^d}\bigl(1+\varepsilon\eta_0(x)\bigr)\geq h_0>0.

$$h_0$$ 越小, 初始波形越接近水深退化的空化状态, 能量的强制性越弱. 因而, 除了证明解存在, 还需定量追踪存在时间对 $\varepsilon,\mu,h_0$ 以及初值大小的依赖. 这也使得 $\varepsilon$ 与 $\mu$ 不预先绑定的分析具有意义: 浅水色散与非线性并非总处于同一量级, 需分别考察 $\varepsilon\lesssim\mu$ 和 $\mu\ll\varepsilon$ 等平衡情形.

已有结果中, Saut、Wang 和 Xu 对该强色散系统先在 $H^s\times H^{s+1/2}$ 、 $$s>d/2$$ 的框架下建立了 $O(1/\sqrt{\varepsilon})$ 时间的适定性; 其能量中的速度高阶导数带有 $\sqrt{\varepsilon}$ 权重, 而拟线性项对 $\nabla v$ 的控制造成了寿命损失. 他们随后利用依赖解的对称化能量, 在 Boussinesq 区域 $\mu\sim\varepsilon$ 、较高正则性及非空化假设下得到 $O(1/\varepsilon)$ 的寿命. 本文试图弥补的空白是: 在将 $\mu$ 和 $\varepsilon$ 作为独立参数的同时, 把线性色散纳入拟线性能量增长控制, 获得更长且对参数及空化余量显式的存在时间.

2. 主要定理或主要结果

主要结论是定理 1.1. 设 $d\in\{1,2\}$ , $0<\varepsilon,\mu\leq1$ , $$s>d/2+2$$ , 且初值 $u_0=(\eta_0,v_0)^\top$ 属于

X_\mu^s(\mathbb R^d),\qquad \|u\|_{X_\mu^s}^2=\|\eta\|_{H^s}^2+\|v\|_{H^s}^2+\mu\|\nabla v\|_{H^s}^2.

假设 $\inf_x(1+\varepsilon\eta_0(x))\geq h_0>0$ , 且当 $$d=2$$ 时 $\operatorname{curl}v_0=0$ . 定义

p_0=1+\varepsilon\|u_0\|_{L^\infty},\qquad \delta_1=0,\quad \delta_2=\delta>0,

其中二维的 $\delta$ 可取充分小. 则初值问题存在唯一解

u=(\eta,v)\in C([0,T_d];X_\mu^s)\cap C^1([0,T_d];X_\mu^{s-1}),

且流映射对初值连续. 论文给出的定量时间及范数控制为

T_d=\left(c\,\mu^{-(1+\delta_d)/4}\varepsilon\,h_0^{-1}p_0\|u_0\|_{X_\mu^s}\right)^{-\frac{4}{4-d+2\delta_d}}, \qquad \sup_{0\leq t\leq T_d}\|u(t)\|_{X_\mu^s} \leq 4\sqrt{\frac{p_0}{h_0}}\,\|u_0\|_{X_\mu^s},

其中 $$c$$ 仅依赖于维数、正则性指数及相关小参数选择. 因此, 在省略由 $$p_0$$ 和初值范数产生的因子后, 一维寿命尺度为

T_1\sim(\mu^{1/4}\varepsilon^{-1}h_0)^{4/3},

二维寿命尺度为

T_2\sim(\mu^{1/4}\varepsilon^{-1}h_0)^{2-},

这里 $$2-$$ 表示通过充分小的 $\delta>0$ 得到略低于二次幂的指数, 并非严格的指数 $$2$$ .

作为证明中的先验估计, 引入修正能量

E_s(u)=\frac12\int_{\mathbb R^d}\left((J^s\eta)^2+|J^sv|^2+\mu|J^s\nabla v|^2+\varepsilon\eta|J^sv|^2\right)\,dx, \qquad J=(1-\Delta)^{1/2}.

若 $1+\varepsilon\eta\geq h_0$ , 则该能量满足与 $X_\mu^s$ 范数相联系的强制性估计, 且其增长受

q_d(t)=\|\nabla\eta(t)\|_{L^\infty}+\|\nabla v(t)\|_{L^\infty} +\sqrt{\mu}\|\Delta v(t)\|_{L^\infty}

控制. 对 $$s>2$$ , 色散估计给出

\int_0^Tq_d(t)\,dt \lesssim \sigma_d(\mu,T)\|u_0\|_{X_\mu^s} +\varepsilon\sigma_d(\mu,T)^2\|u\|_{L_T^\infty X_\mu^s}^2, \quad \sigma_d(\mu,T)=\mu^{-(1+\delta_d)/4}T^{(4-d+2\delta_d)/4}.

与此相配合, $$s>3$$ 时还可对包含 $\mu\|\nabla\cdot\Delta v\|_{L^\infty}$ 的量得到同型积分控制, 用以保持非空化条件. 这些估计经自举论证导出定理中的寿命与解范数界. 在 Boussinesq 区域 $\mu\sim\varepsilon$ 中, 一维尺度化为 $\varepsilon^{-1}h_0^{4/3}$ , 与既有 $O(1/\varepsilon)$ 结果同阶且显示了 $$h_0$$ 依赖; 二维则给出 $\varepsilon^{-3/2+}h_0^{2-}$ 型改进尺度. 这些尺度描述需结合定理中显式的初值范数和 $$p_0$$ 因子理解.

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

证明的关键困难是同时处理两个彼此不同的控制需求. 一方面, 能量必须控制 $\eta,v$ 的 $$H^s$$ 范数以及带 $\sqrt{\mu}$ 权重的速度导数 $\sqrt{\mu}\nabla v$ ; 另一方面, 拟线性能量不等式需要时间积分控制 $\nabla\eta,\nabla v,\sqrt{\mu}\Delta v$ 的 $L^\infty$ 范数. 单靠能量估计无法直接得到这些 Lipschitz 型量的时间可积性, 而若忽略 $\mu$ 权重, 也无法在 $\mu$ 与 $\varepsilon$ 独立变化时保持估计尺度正确.

作者为独立参数设计了修正能量 $$E_s$$ . 其中 $\mu|J^s\nabla v|^2$ 与 $X_\mu^s$ 的定义匹配, 非线性权重项 $\varepsilon\eta|J^sv|^2$ 则将水深因子 $1+\varepsilon\eta$ 纳入能量主部. 在非空化条件下, 该能量满足形如

h_0\|u(t)\|_{X_\mu^s}^2\lesssim E_s(u(t)) \lesssim p_d(t)\|u(t)\|_{X_\mu^s}^2, \qquad \frac{d}{dt}E_s(u)\lesssim\varepsilon h_0^{-1}p_d(t)q_d(t)E_s(u),

的强制性及增长估计, 其中 $p_d(t)=1+\varepsilon(\|\eta(t)\|_{L^\infty}+\|v(t)\|_{L^\infty})$ . 通过分部积分、乘积估计及 Kato–Ponce 型交换子控制高阶微分算子与非线性乘积的交换误差, 得以将能量增长归结为 $$q_d$$ 的时间积分. 同一结构也用于解差的方程, 构造 $$E_0(z)$$ 并建立唯一性所需的低阶稳定性估计.

补足 Lipschitz 积分控制的核心是线性化系统的频率局部化色散估计. 相关传播群可用符号

m_\mu(D)=|D|J_\mu(D),\qquad J_\mu(D)=(1-\mu\Delta)^{1/2}

表示. 高频时 $m_\mu(\xi)\sim\sqrt{\mu}|\xi|^2$ , 因而传播具有 Schrödinger 型振荡, 但色散常数显式依赖 $\mu$ 和频率. 作者先证明频率局部化的 $L^1\to L^\infty$ 衰减, 再以 $$TT^*$$ 推导 Strichartz 估计. 一维使用

\|S_\mu(t)P_\lambda f\|_{L_t^4L_x^\infty} \lesssim \mu^{-1/4}\lambda^{-1/4}\langle\sqrt{\mu}\lambda\rangle^{1/4}\|f\|_{L^2},

二维使用 $0<\delta<1$ 时的

\|S_\mu(t)P_\lambda f\|_{L_t^{2/(1-\delta)}L_x^{2/\delta}} \lesssim\mu^{-(1+\delta)/4}\|f\|_{L^2}.

随后将线性传播作用于 Duhamel 公式的各频率块, 并结合 Littlewood–Paley 分解、Bernstein 不等式和双线性 Besov 乘积估计, 控制非线性强迫项. 由此得到 $$q_d$$ 的积分界, 其时间增长指数正是 $\sigma_d(\mu,T)$ 中的指数; 二维额外引入小的 $\delta$ 以适配相应 Strichartz 指数与求和, 因而最终只能给出 $$2-$$ 型寿命幂次.

最后的自举不仅要控制 Sobolev 范数, 还须防止水深因子在演化中降至零. 作者同时控制 $$p_d(t)$$ 和 $1+\varepsilon\eta(x,t)$ : 由原方程估计 $\eta_t,v_t$ , 并使用对 $Q_d=q_d+\mu\|\nabla\cdot\Delta v\|_{L^\infty}$ 的积分界, 使 $$p_d$$ 保持在初始值常数倍以内, 且在所选区间上有 $1+\varepsilon\eta\geq h_0/2$ . 时间 $$T_d$$ 的选取使 $\varepsilon\ell_0^2\sigma_d(\mu,T_d)\|u_0\|_{X_\mu^s}$ 足够小, 其中 $\ell_0=(p_0/h_0)^{1/2}$ . 因此能量指数增长可以被吸收到自举假设中. 对适定性其余部分, 文中说明通过正则化与紧性获得解, 以 Bona–Smith 论证连续依赖, 并用解差能量估计证明唯一性; 细节在摘录中被略去.

4. 与之前工作的比较

与直接能量法相比, Saut、Wang 和 Xu 在文中讨论的早期结果适用于 $u_0\in H^s\times H^{s+1/2}$ 、 $$s>d/2$$ , 得到 $O(1/\sqrt{\varepsilon})$ 的存在时间. 其使用常系数权重 $\operatorname{diag}(1,1-\varepsilon\Delta)$ 使线性部分呈反自伴结构, 再配合 Kato–Ponce 交换子估计及 Bona–Smith 方法. 该方法的寿命损失来自拟线性项 $\varepsilon\eta\nabla\cdot v$ : 估计中需要控制 $\nabla v$ 于 $$H^s$$ , 但能量只控制 $\sqrt{\varepsilon}\nabla v$ . 本文不是简单加强该能量估计, 而是额外利用色散产生的时空可积性来控制相关 $L^\infty$ 导数, 同时将高阶速度导数的权重按独立的色散参数改为 $\sqrt{\mu}$ .

更直接的比较对象是 Saut、Wang 和 Xu 对同一强色散系统的长时间定理. 在 $\mu\sim\varepsilon$ 的 Boussinesq 区域, 他们利用依赖解的对称化能量, 在 $$s>2+d/2$$ 、初值具有额外速度正则性条件

\|\eta_0\|_{H^s}+\|v_0\|_{H^s}+\sqrt{\varepsilon}\|v_0\|_{H^{s+1}}<\infty

并满足非空化条件时, 得到 $[0,T/\varepsilon]$ 上的适定性, 其中 $$T$$ 与 $\varepsilon$ 无关. 本文将 $\mu,\varepsilon$ 作为独立参数, 初值空间写为 $X_\mu^s$ , 即以 $\sqrt{\mu}\nabla v$ 加权控制速度的额外导数, 并由色散估计给出显式参数寿命. 它因此不仅覆盖 $\mu\sim\varepsilon$ 的量级比较, 也量化 $\varepsilon\lesssim\mu$ 和 $\mu\ll\varepsilon$ 等非平衡情形下非线性与色散的竞争. 对这些区域具体取值后的比较, 应依据定理中的完整公式而非仅依据省略初值因子的摘要尺度.

在 $\mu\sim\varepsilon$ 时, 一维结果的参数尺度为 $\varepsilon^{-1}h_0^{4/3}$ , 所以对固定 $$h_0$$ 与有界初值, 与先前 $O(1/\varepsilon)$ 结果同阶, 并未宣称在该维数、该参数区间获得更高的 $\varepsilon$ 幂次. 新增信息是对空化余量的显式依赖, 以及对独立参数的量化. 二维结果则为 $\varepsilon^{-3/2+}h_0^{2-}$ 型尺度, 摘录明确称其改进了 [28] 的相应寿命. 此处的 $$-$$ 表示存在一个小的指数损失, 因而不能将结果表述为无损的 $\varepsilon^{-3/2}$ 或精确的二次 $$h_0$$ 幂次.

对空化条件的处理也更定量. 既有长时间结论要求 $1+\varepsilon\eta_0\geq h_0$ , 但摘要所述先前结果没有在寿命尺度中展示其如何随 $$h_0$$ 退化. 本文的能量强制性含 $$h_0$$ , 自举时间中相应出现 $$h_0$$ 及 $$p_0$$ 、初值范数因子; 解在短时间内保有 $1+\varepsilon\eta\geq h_0/2$ . 这既明确显示接近空化时估计变差, 也表明结论并未覆盖 $$h_0=0$$ 或空化发生的情形.

改进的代价和适用边界同样重要. 定理要求 $$s>d/2+2$$ , 初值属于参数依赖空间 $X_\mu^s$ , 二维需无旋条件, 且寿命常数依赖初值大小、 $$p_0$$ 及非空化下界. 二维还保留任意小的 $\delta>0$ 所造成的指数损失. 论文以色散估计为关键, 其结论针对所给强色散系统及 $$d=1,2$$ , 不能仅凭摘录推广到整个 $$(abcd)$$ 家族或更一般的底部几何. 文中列举了其他 Boussinesq 与 Whitham–Boussinesq 长时间研究, 以及 SQG 和 Whitham–Boussinesq 中类似色散-能量思路, 但提供内容未展开这些工作的具体定理假设和定量比较; 因此更细致的跨模型优劣判断从提供内容中无法确认.

4. 广义 SQG 方程的适定性与极限行为

On wellposedness and limiting behavior of generalized SQG equations

作者: Anuj Kumar

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39580v1

摘要翻译: 我们研究定义在 $\mathbb{R}^2$ 上的二维广义表面准地转方程：

Comment: 20 pages

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

1. 研究对象和背景

本文研究二维广义表面准地转方程 (generalized surface quasi-geostrophic equation, gSQG) 在全空间 $\mathbb{R}^2$ 上的光滑解，以及解对 constitutive law 参数 $\beta$ 的稳定性。方程可写成输运形式

\begin{cases} \partial_t\theta+u\cdot\nabla\theta=0,\\ u=-\nabla^\perp\Lambda^{\beta-2}\theta, \end{cases} \qquad x\in\mathbb{R}^2,\quad t>0,\quad \beta\in[1,2),

其中 $\theta(x,t)$ 是被输运的标量， $$u$$ 是由 $\theta$ 非局部决定的速度， $\Lambda=(-\Delta)^{1/2}$ ，且 $\nabla^\perp=(-\partial_{x_2},\partial_{x_1})$ 。等价地，引入流函数 $\psi$ ，有 $u=\nabla^\perp\psi$ 、 $\Delta\psi=\Lambda^\beta\theta$ 。由于 $\nabla\cdot u=0$ ，标量沿速度场输运；但速度与标量之间的非局部关系具有奇异性，且奇异程度随 $\beta$ 改变。 $\beta=1$ 时即为 SQG 方程， $\beta>1$ 时 constitutive law 比 SQG 更奇异。本文考察给定初值 $\theta(x,0)=\theta_0(x)$ 时，方程的光滑解在固定时间区间内是否能在参数微扰下继续存在，以及对应解是否定量接近。

这类方程一方面来自地球物理中的旋转、稳定分层流体模型。SQG 标量可解释为表面温度或浮力，其速度反演包含非局部算子；另一方面，SQG 的涡旋拉伸特征与三维 Euler 涡量方程有相似之处，因此也常被用作研究输运方程正则性问题的模型。需要区分的是，本文讨论的是无耗散 gSQG 的光滑解及参数依赖，并非证明耗散 SQG 或无耗散 gSQG 的全局正则性。

在 $\beta$ 的其他区间，gSQG 与熟知模型相连： $\beta=0$ 对应二维 Euler 涡量形式， $\beta=1$ 对应 SQG。本文只处理 $\beta\in[1,2)$ ，其中特别关注从 SQG 参数向更奇异的 $\beta>1$ 一侧扰动。对 $\beta>1$ 的 Cauchy 问题，既有工作已给出局部适定性结果，例如在 $$H^4$$ 中建立局部适定性，后来将初值正则性改进到 $H^{\beta+1+\varepsilon}$ 一类条件；文中还回顾了耗散版本及临界空间中的相关研究。与此同时，临界 Sobolev 正则性附近的适定性问题已有若干病态性结果，说明控制非线性并不能简单依赖常规乘积估计，而需要利用方程结构。

本文的核心问题可以精确表述为：固定同一个初值 $\theta_0$ ，假设参数为 $\beta_0\in[1,2)$ 时存在光滑解 $\theta_{\beta_0}$ 于 $$[0,T]$$ ，是否存在 $\delta>0$ ，使所有足够接近 $\beta_0$ 的允许参数 $\beta$ 对应的解也存在于整个 $$[0,T]$$ ？并且，能否给出 $\theta_\beta-\theta_{\beta_0}$ 的定量 Sobolev 估计？这里的“存在区间不缩短”是相对于已知的参考解时间区间而言，并不表示存在统一的全局时间，也不推出任意初值的全局正则性。

参数比较的困难在于，两解之差的速度并不能在 $\beta>1$ 时直接通过简单的分数积分估计从标量差控制；特别是常用的对易子估计常数在 $\beta\to1^+$ 时发散。因而，单纯把方程视为关于参数连续变化的常规拟线性输运方程，无法覆盖 SQG 端点附近的统一控制。本文试图弥补的空白，是在更奇异的 $\beta>1$ 区域建立“已知光滑解可延续到邻近参数”的结论，并以对 $\beta\to1^+$ 一致的估计处理 SQG 端点。作者还将此前关于 $\beta\to1^-$ 的结果与本工作相结合，指出 SQG 两侧的小参数扰动稳定性由此得到补全；关于先前结果的具体定理条件仍应以相应文献为准。

2. 主要定理或主要结果

主要结论是定理 1.1 和定理 1.2，分别处理参考参数 $\beta_0\in(1,2)$ 与 SQG 端点 $\beta_0=1$ 。两者讨论相同初值下参数变化导致的解的延续与差值估计，结论均为局部时间区间 $$[0,T]$$ 上的光滑解存在，而非全局适定性。

定理 1.1：设 $\beta_0\in(1,2)$ ， $$s>3$$ ，且 $\theta_{\beta_0}$ 是初值 $\theta_0\in H^{s+1}(\mathbb{R}^2)$ 所对应的 gSQG 光滑解，并在 $$[0,T]$$ 上存在。则存在 $\delta>0$ ，其依赖于 $$T$$ 和参考解范数积分

\int_0^T\|\theta_{\beta_0}(t)\|_{H^{s+1}}\,dt,

使得对满足 $\beta\in(1,2)$ 且 $|\beta-\beta_0|<\delta$ 的参数，使用同一初值 $\theta_0$ 的解 $\theta_\beta$ 也在 $$[0,T]$$ 上光滑存在。并且，对应解满足

\|\theta_\beta(t)-\theta_{\beta_0}(t)\|_{H^s} \le C\left( |\beta-\beta_0|^{2-\beta} +|\beta-\beta_0|^{2-\beta_0} +|\beta-\beta_0|\,|\ln|\beta-\beta_0|| \right).

这里 $$C>0$$ 依赖于 $$T$$ 及 $\int_0^T\|\theta_{\beta_0}\|_{H^{s+1}}\,dt$ 。该估计给出固定区间上的 $$H^s$$ 连续依赖定量模；其中幂次含有变化参数 $\beta$ ，不能在未指定参数邻域及常数控制的情况下将其改写为单一、与参数无关的幂次阶。

定理 1.2：设 $\beta_0=1$ ， $1<\beta<2$ ，且参考 SQG 解以初值 $\theta_0\in H^{s+1}(\mathbb{R}^2)\cap L^1(\mathbb{R}^2)$ 、 $$s>3$$ 在 $$[0,T]$$ 上存在。则存在 $\delta>0$ ，依赖于 $$T$$ 和

\int_0^T\|\theta_{1}(t)\|_{H^{s+1}\cap L^1}\,dt,

使得当 $0<\beta-1\le\delta$ 时，使用同一初值的 $\theta_\beta$ 在 $$[0,T]$$ 上也光滑存在，并有

\|\theta_\beta(t)-\theta_1(t)\|_{H^s} \le C\left((\beta-1)^{2-\beta} +(\beta-1)|\ln|\beta-1||\right).

此处 $$C$$ 依赖于 $$T$$ 和 $\int_0^T\|\theta_1\|_{H^{s+1}\cap L^1}\,dt$ 。与定理 1.1 相比，端点结论需要额外的 $$L^1$$ 控制；参数只从 $$1$$ 的右侧趋近，因为本文处理的是 $\beta>1$ 。当 $\beta\downarrow1$ 时，右侧两项均趋于零，因而给出该方向的定量收敛。

证明中还建立了对易子引理 3.2。对 $\beta\in(1,2)$ 、 $p\in(1,\infty)$ ，若 $$1/q_1+1/r_1=1/p$$ ，则

\|[\nabla\Lambda^{\beta-2}\cdot,g]f\|_{L^p} \le C_2(\beta)\|\nabla g\|_{L^{q_1}}\|f\|_{L^{r_1}} +C_3(\beta,p)\|g\|_{L^\infty}\|f\|_{L^1}.

在另一个指数条件 $$1/q_2+1/r_2>1/p$$ 下，还证明

\begin{aligned} \|[\nabla\Lambda^{\beta-2}\cdot,g]f\|_{L^p} \le {}& C_2(\beta)\|\nabla g\|_{L^{q_1}}\|f\|_{L^{r_1}} +C_4(\beta)\|g\|_{L^\infty}\|f\|_{\dot W^{-1,p}}\\ &+C_5(\beta,p,q_2,r_2)\|\nabla g\|_{L^{q_2}}\|f\|_{\dot W^{-1,r_2}}. \end{aligned}

上述常数在 $\beta\to1^+$ 时一致有界，这是端点证明的重要技术结论。另一个关键辅助结果是引理 4.1：对 $\sigma>3$ ，修改通量方程的非线性满足

\left|\left\langle\operatorname{div}F_q^\beta(\theta),(-\Delta)^\sigma\theta\right\rangle\right| \le C_6(\beta,\sigma)\|q\|_{H^\sigma}\|\theta\|_{H^\sigma}^2,

且 $$C_6$$ 在 $\beta\to1^+$ 时一致有界。该估计用于高阶能量控制；引理 4.2 则提供相应双线性项的连续性估计。提供的 PDF 摘录未展示定理证明的全部后续细节，故此处不补写未出现在摘录中的中间估计或额外定理。

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

证明的主要障碍不是单独获得每个参数方程的局部解，而是比较两个不同参数的解时，速度反演算子阶数也发生变化。记 $\theta=\theta_\beta-\theta_{\beta_0}$ 、 $u=u_\beta-u_{\beta_0}$ ，差值方程中既有由 $\theta$ 引起的非线性项，也有由算子 $\nabla^\perp\Lambda^{\beta-2}$ 与 $\nabla^\perp\Lambda^{\beta_0-2}$ 不同产生的项。对于 $\beta>1$ ，速度反演具有更强奇异性，直接以 $\|\theta\|_{H^s}$ 控制速度差并不足够；而文中回顾的已有对易子估计，其常数在 $\beta\to1^+$ 发散，不能直接用于端点参数扰动。

作者的核心结构性处理，是将 gSQG 改写为带修改通量的守恒律。对足够光滑的 $$q$$ ，定义

F_q^\beta(\theta) =(\nabla^\perp\Lambda^{\beta-2}q)\theta +\Lambda^{\beta-2}((\nabla^\perp\theta)q), \qquad \partial_t\theta+\operatorname{div}F_q^\beta(\theta)=0.

取 $q=-\theta$ 时，形式上有

\operatorname{div}F_{-\theta}^\beta(\theta) =-(\nabla^\perp\Lambda^{\beta-2}\theta)\cdot\nabla\theta =u\cdot\nabla\theta,

从而恢复原方程。此重写并非仅是记号变化：它把原来难以直接估计的输运非线性组织成具有抵消结构的通量，使能量配对时可以把关键部分改写为对易子。作者借鉴并发展了先前用于局部适定性证明的修改通量思路，将其用于不同参数解之间的比较。

具体而言，高阶能量估计考察 $\operatorname{div}F_q^\beta(\theta)$ 与 $(-\Delta)^\sigma\theta$ 的配对。若直接对乘积逐项估计，导数会落在非局部速度和标量上，可能产生无法由给定 Sobolev 正则性吸收的项。作者将配对拆成若干部分：利用散度自由结构，使某些主项积分为零；将剩余项写成形如 $[\Lambda^\alpha,a]b$ 的对易子；再通过 Kato–Ponce 型对易子估计以及 Sobolev 嵌入分配导数。引理 4.1 总结了这一处理的能量结果，给出 $\|q\|_{H^\sigma}\|\theta\|_{H^\sigma}^2$ 型控制，并确保相应常数在 $\beta\to1^+$ 时不退化。摘录中展示的 $$I_2$$ 、 $$I_4$$ 项正是借助积分分部和散度结构消失的例子；其他项则通过逐次减去主导乘积项形成对易子。

为解决端点附近传统估计的常数退化，作者另证引理 3.2。一般已有估计将对易子控制为 $\|f\|_{L^p}\|\Lambda^{\beta-1}g\|_{L^\infty}$ ，但相关常数随 $\beta\to1^+$ 失去一致性。新证明在空间变量 $$|x-y|$$ 上作近场与远场分解：近场使用 $$g(y)-g(x)$$ 的一阶差分及中值定理，将核的奇异性转化为 $\nabla g$ ；二维局部核的积分给出包含 $(2-\beta)^{-1}$ 的量，在 $\beta\to1^+$ 时保持有界。远场则不再沿用同一估计，而分别采用 $f\in L^1$ 的卷积控制，或将导数转移到 $\Lambda^{-2}f$ 上，以 $\dot W^{-1,p}$ 和 $\dot W^{-1,r_2}$ 范数控制。后一处理要求相应指数满足 $$1/q_2+1/r_2>1/p$$ 。这两种远场界对应引理 3.2 的两条不等式，便于在不同能量项中选用适配的可积性或负阶正则性。

从证明架构看，作者先建立对参数端点一致的对易子控制与修改通量非线性能量界，再将其用于比较方程的能量估计。参数变化产生的误差项需要结合算子阶数变化进行估计；获得差值控制后，再借助关于非负函数的微分不等式引理 2.3，将形如“参数小量源项加二次能量项”的估计闭合。该引理说明，当小参数 $\nu$ 足够小时，若

y'(t)\le \nu F(t)+G y(t)^2,\qquad y(0)=0,

则在规定的小性条件下有 $y(t)\lesssim \nu\int_0^T F(t)\,dt$ 。这种闭合机制把参数扰动误差与参考解的时间积分范数联系起来，并用于保证差值在整个 $$[0,T]$$ 上受控，从而延续邻近参数的光滑解。就提供的摘录而言，差值方程所有项如何逐项对应至定理中的三种误差率并未完整出现，因此不进一步推断其每一项的来源或隐藏常数的精细参数依赖。

端点 $\beta_0=1$ 还需要额外的 $$L^1$$ 信息。引理 3.2 的第一种远场估计直接含有 $\|f\|_{L^1}$ ，而定理 1.2 假设参考初值属于 $H^{s+1}\cap L^1$ ，并在时间积分控制中保留相应范数。这表明端点所需的额外可积性与克服对易子估计退化的远场控制相配合；但摘录不足以断言 $$L^1$$ 假设在所有可能的证明方案中都是必要条件。方法上的实质推进是将“参数扰动导致的奇异算子差”纳入带抵消的通量框架，并用端点一致对易子估计处理常规估计无法覆盖的 $\beta\to1^+$ 极限。

4. 与之前工作的比较

与既有工作相比，本文的主要推进是从“固定参数下的适定性”转向“光滑解存在时间对参数扰动的稳定性”。文中回顾， $\beta\in(1,2)$ 的 Cauchy 问题先在 $$H^4$$ 中建立局部适定性，之后相关结果将初值要求改进至 $H^{\beta+1+\varepsilon}$ 。这些结果回答给定 $\beta$ 时是否局部存在解；本文则假设某一参考参数 $\beta_0$ 的解已在给定 $$[0,T]$$ 上存在，证明相邻参数的同初值解也覆盖该区间，并给出 $$H^s$$ 差值的显式上界。它不取代各固定参数下的适定性理论，也不把有限时间结论升级为全局正则性。

与文中明确提到的 [37] 相比，本文处理了此前结果没有覆盖的奇异方向。作者称 [37] 建立了 $\beta_0\in[0,1]$ 、 $\beta\in(0,1)$ 的相应参数扰动结论，即在其所研究范围内处理向较低参数一侧的扰动；本文给出 $\beta_0\in(1,2)$ 邻域内的结果，并特别证明从 SQG 的 $\beta_0=1$ 向 $\beta>1$ 的右侧稳定性。因此，按作者的表述，本文与 [37] 合起来补全了 SQG 在参数 $\beta=1$ 两侧的小扰动分析。两侧结果的具体正则性假设、差值范数与误差率是否完全对称，不能仅由本文摘录确认；尤其本文的端点定理明确要求 $H^{s+1}\cap L^1$ ，不能据此把右侧结论说成与左侧条件完全一致。

证明技术也与 [37] 所采用的策略不同。本文指出，先前工作在某些较弱奇异性区间可将速度差分解后直接使用 Hardy–Littlewood–Sobolev 不等式，而在 $\beta_0\in(1,2)$ 或 $\beta_0=1$ 且 $\beta\to1^+$ 时，constitutive law 更奇异，这种直接控制不再可用。本文改为把非线性重写成修改通量并使用对易子抵消；同时针对传统对易子界在 $\beta\to1^+$ 时常数发散的问题，建立了常数一致有界的新估计。这一改变使参数稳定性论证能覆盖更奇异的一侧，但端点结果付出的条件代价是加入 $$L^1$$ 控制，并且初值仍需有较高 Sobolev 正则性 $H^{s+1}$ 、 $$s>3$$ 。

本文还与 [10] 的思路有概念上的类比：[10] 研究若三维不可压 Euler 方程有光滑解，则同初值的三维 Navier–Stokes 方程在黏性足够小时也在相同时间区间存在。本文把“参考模型有光滑解，则相邻模型保持该存在区间”的稳定性问题用于 gSQG 参数 $\beta$ 。类比只在问题结构上成立：本文扰动的是速度反演算子的阶数，而不是加入黏性耗散；不能由此将两种结果的估计机制或假设视为相同。

此外，文中回顾了修改通量方法在作者与合作者先前工作 [25] 中的应用：该方法用于带分数耗散的 gSQG，在临界 Sobolev 空间中建立大数据局部适定性、小数据全局适定性及 Gevrey 正则性；[24] 则用于对数正则化版本的局部适定性。本文的推进不在于首次提出修改通量，而在于将此结构用于参数连续性，并与新的端点一致对易子界结合。作者也提及 $\beta>1$ 的无耗散方程已有局部适定性研究，以及临界正则性附近的病态性结果；这些工作与本文问题不同：本文在较高正则性和已存在参考光滑解的前提下证明有限时间稳定，不解决临界空间的适定性或病态性。

适用范围与限制应据定理准确理解。对内部参数 $\beta_0\in(1,2)$ ，可在其邻域内对两侧参数作扰动，但始终要求新参数 $\beta\in(1,2)$ ；对端点 $\beta_0=1$ ，本文只证明 $\beta>1$ 的右侧邻域。参数上界 $$2$$ 不包含在结论内。时间区间由参考解给定，阈值 $\delta$ 依赖于 $$T$$ 及参考解的时间积分范数；因此这不是对所有初值统一的稳定性结论。定量结论是在 $$H^s$$ 中控制解差，且端点需额外 $$L^1$$ 假设。文中摘要和定理未声称全局存在、低正则性延拓、解对参数的可微性，或在 $\beta=2$ 极限处的结论，这些都不应从本文结果外推。

5. 密度-涡量反馈下稳定二维 Boussinesq 方程组在 $$H^2$$ 中的非存在性

Non Existence of the Stable 2D Boussinesq Equations in $$H^2$$ via Density-Vorticity Feedback

作者: Roberta Bianchini, Luis Martínez-Zoroa

类别: math.AP

发布日期: 2026-09-29

arXiv ID: 2609.38584v1

摘要翻译: 我们证明了二维 Boussinesq 方程在关于线性稳定背景剖面 $$-x_2$$ 的情形下发生瞬时 $$H^2$$ 范数爆炸。以往关于该系统的不适定性结果依赖于 $L^\infty$ 型空间，而本文在能量临界空间中给出了首个强不适定性结果；在该空间中，稳定背景分层会抑制形变增长。

Comment: 无

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

1. 研究对象和背景

本文研究二维无粘 Boussinesq 方程在稳定背景密度剖面 $\bar\rho(x)=-x_2$ 附近的强适定性问题。取重力加速度为 $$1$$ ，原始变量形式为

\begin{cases} \partial_t\rho+u\cdot\nabla\rho=0,\\ \partial_tu+u\cdot\nabla u+\nabla P=(0,-1)\rho,\\ \nabla\cdot u=0. \end{cases}

令 $\omega=\partial_1u_2-\partial_2u_1$ ，并用流函数表示速度： $u=\nabla^\perp\psi=(-\partial_2\psi,\partial_1\psi)$ 、 $\Delta\psi=\omega$ 。相应的涡量形式为

\begin{cases} \partial_t\rho+u[\omega]\cdot\nabla\rho=0,\\ \partial_t\omega+u[\omega]\cdot\nabla\omega=-\partial_1\rho,\\ u[\omega]=K*\omega, \end{cases} \qquad K(x)=\dfrac{1}{2\pi}\left(-\dfrac{x_2}{|x|^2},\dfrac{x_1}{|x|^2}\right).

背景态为 $(\bar\rho,\bar u)=(-x_2,0)$ 。作者把总密度写为 $\rho=-x_2+\rho_{\mathrm{pert}}$ ，研究的是扰动 $\rho_{\mathrm{pert}}$ 和涡量，而不是总密度本身的有限能量性质。稳定背景的线性作用通常抑制速度形变增长，因此，能否在这一稳定化效应下产生临界正则性失控，是问题的关键。本文所说的主要失控表现为：初值扰动可以任意小地属于 $$H^2$$ ，但任意正时间的密度扰动不再属于 $$H^2$$ ，且涡量不再属于 $$H^1$$ 。

这一问题与 Boussinesq 方程的整体正则性及小尺度形成研究相关，但逻辑上不同于证明光滑解有限时间爆破。作者构造的是从小扰动出发的强解，并证明其在每个正时间都发生临界 Sobolev 范数爆炸；定理没有声称初始数据本身不光滑，也没有证明解在某个有限正时间之前不存在。解仍有较低阶正则性，并在作者给出的能量解类中唯一。因而“ $$H^2$$ 中的非存在性”指无法在所要求的临界正则性类中持续存在，而不是完全没有解。

临界正则性失适定在流体方程中已有重要先例。Bourgain 与 Li 的 Euler 构造利用速度形变产生 Sobolev 范数失控；二维 Boussinesq 方程在 $L^\infty$ 型临界空间中的失适定性也已有研究，Elgindi 与 Masmoudi 的结果与 Riesz 变换在 $L^\infty$ 上无界有关，后续工作进一步得到范数膨胀。本文强调的空白是：这些 Boussinesq 结果依赖 $L^\infty$ 型框架，而本文试图在能量临界的 $$H^2$$ 扰动尺度上，在稳定背景的抑制作用下证明更强的瞬时失控，并且不把不稳定性归结为 Euler 子系统。

稳定背景带来一个直接障碍：若密度扰动在形变发生的中心附近保留背景的稳定梯度，背景效应会削弱所需的形变。作者据此安排初值，使 $\rho_{\mathrm{pert}}=x_2$ 在原点邻域成立，从而总密度 $-x_2+\rho_{\mathrm{pert}}$ 在该邻域为零，局部中和稳定分层；与此同时，扰动仍可具有任意小的 $$H^2$$ 范数。由于密度方程是输运方程，初始密度扰动的支集与原点保持正距离一段时间，这使局部中和结构可以被利用。

本构造的第二个核心特征是取初始涡量恒为零。这样一来，初始速度为零，形变并非来自预置的 Euler 型速度场，而是先由密度梯度通过涡量方程 $\partial_t\omega=-\partial_1\rho$ 生成涡量，再由 Biot-Savart 算子生成速度形变；该形变又影响密度输运并改变后续涡量生成。要证明失控，必须控制这一密度—涡量反馈在时间演化中的符号与积累，而不能仅将密度视作被动标量。

技术上，论文先通过对称性确保 $$u(0,t)=0$$ ，把原点的速度形变 $\partial_1u_1(0,t)$ 作为可追踪量。随后利用多尺度、局部化的密度扰动制造负的形变，并逐层叠加以累积任意强的原点形变。由于各尺度会通过非局部 Biot-Savart 速度相互作用，简单叠加并不是精确解；作者因此引入内外近似及误差控制，再以迭代构造建立极限强解。论文试图补足的正是稳定背景、临界 $$H^2$$ 正则性、零初始涡量和非线性反馈这几项同时存在时的强失适定性结果。所给摘录未包含论文后半部分的完整估计和迭代参数条件，故无法据此进一步确认所有中间误差界的参数依赖。

2. 主要定理或主要结果

主要结论是定理 1.1。对任意 $\epsilon>0$ ，存在时间 $T_\epsilon>0$ 和 $\rho_{\mathrm{in}}\in H^2(\mathbb R^2)$ ，满足

\|\rho_{\mathrm{in}}\|_{H^2}\leq\epsilon,

使得以

\rho(x,0)=-x_2+\rho_{\mathrm{in}}(x),\qquad \omega(x,0)=0

为初值的稳定二维 Boussinesq 方程在 $[0,T_\epsilon]$ 上存在强解。对每个 $t\in(0,T_\epsilon]$ ，该解满足

\|\rho(t)+x_2\|_{H^2(\mathbb R^2)}=\infty, \qquad \|\omega(t)\|_{H^1(\mathbb R^2)}=\infty.

因此，尽管初始密度扰动可任意小地属于 $$H^2$$ 且初始涡量为零，任意正时间的解都不属于密度扰动 $$H^2$$ 、涡量 $$H^1$$ 的临界正则性类。这里结论是对所有 $t\in(0,T_\epsilon]$ 成立，而不是仅在某个趋于零的时间序列上成立；不过摘录没有给出 $T_\epsilon$ 关于 $\epsilon$ 的明确下界或渐近关系，不能据此声称统一的爆炸时间尺度。

解仍保有较低阶强正则性。定理具体给出

(\rho+x_2,\omega) \in C\!\left([0,T_\epsilon];(H^1\cap C^1)\times(L^2\cap C^0)\right) \cap C^1\!\left([0,T_\epsilon];L^2\times\dot H^{-1}\right).

因此， $H^2\times H^1$ 的瞬时范数爆炸并非解在所有意义下失去正则性。作者还指出解在原点之外具有更高的经典性：密度方程处处按经典意义成立，涡量方程在 $\mathbb R^2\setminus\{0\}$ 上经典成立；障碍集中在原点附近涡量的空间可微性。由对称性，原点是静止点，且文中指出 $\partial_t\omega(0,t)=0$ 、 $u[\omega](0,t)=0$ ，所以沿拉格朗日轨道理解时仍有相应方程关系。这里应区分强解结论与处处经典解结论，作者并未声称构造出全空间经典解。

定理还给出唯一性：解在

C([0,T_\epsilon];H^1\times L^2)\cap C^1([0,T_\epsilon];L^2\times\dot H^{-1})

这一能量解类中唯一。因而范数爆炸不是通过选择非唯一分支获得的，而发生在从给定初值出发的唯一能量类解上。

构造稳定背景初值时，作者使用一个局部抵消性质：对任意足够小的 $\epsilon_0>0$ ，可取对称光滑函数 $\rho_{\mathrm{in}}$ ，使其支集包含于 $$B_1(0)$$ ，在某个 $B_{\delta_0}(0)$ 中满足 $\rho_{\mathrm{in}}(x)=x_2$ ，且 $\|\rho_{\mathrm{in}}\|_{H^2}\leq\epsilon_0$ 。于是总密度在该邻域满足 $-x_2+\rho_{\mathrm{in}}=0$ 。后续构造进一步安排初始扰动支集避开原点，使这一空穴结构在短时间内保持。

引理 2.3 的关键初始形变结论是：可在保持 $\|\rho_{\mathrm{pert}}(0)\|_{H^2}\leq\epsilon_0$ 、 $\omega(0)=0$ 的条件下，使某个短时间区间上

\partial_t\partial_1u_1[\omega](0,t)<0.

其机制可以从原点形变的 Biot-Savart 公式看出：

\partial_1u_1[\omega](0,t) =\frac{1}{\pi}\int_{\mathbb R^2}\frac{y_1y_2}{|y|^4}\,\omega(t,y)\,dy.

初始涡量为零时，形变的时间导数由涡量源项 $-\partial_1\rho$ 决定。通过选取不同尺度、符号相容的密度层，可使该时间导数为负且绝对值任意大；摘录中的具体多尺度形式为

\rho_{K,\epsilon}(x)=\epsilon\sum_{i=1}^K\frac{f(2^ix)}{i2^i},

其中 $$f$$ 光滑、具指定对称性且支集位于环带 $B_2(0)\setminus B_1(0)$ 。各层支集互不相交，其 $$H^2$$ 范数平方受 $\epsilon^2\sum_{i=1}^K i^{-2}$ 控制，而对原点形变导数的贡献包含调和和 $\sum_{i=1}^K i^{-1}$ ，故可在小 $$H^2$$ 代价下产生大形变。

引理 2.4 针对线性形变场 $$k(t)(x_1,-x_2)$$ 给出密度—涡量反馈的定量估计。若 $$k(t)<0$$ 、 $-\int_0^T k(s)\,ds\leq M$ ，且初始密度支集满足 $|x_2|\geq2|x_1|$ 以及 $x_1x_2\partial_1\widetilde\rho_{\mathrm{in}}\geq0$ ，则生成的线性模型解满足

-\partial_1u_1[\widetilde\omega](0,t) \asymp_M \partial_1u_1[\partial_1\widetilde\rho_{\mathrm{in}}](0,0)\,t, \qquad 0\leq t\leq T,

其中隐含常数依赖于 $$M$$ 。更明确地，证明中的流映射为

\Phi(x,t)=\left(x_1e^{\int_0^t k(\tau)\,d\tau},\ x_2e^{-\int_0^t k(\tau)\,d\tau}\right),

并且涡量沿流由输运密度的 $$x_1$$ 导数产生。由扇形支集条件，Biot-Savart 核在形变下的分母可作统一比较，从而得到形变的符号保持与时间线性增长估计。引理 2.5 又把早期区间的形变积分与全区间积分联系起来，避免每次迭代新增的形变只集中在极短时间内而对总形变贡献不足。

论文使用内外近似的粘合引理 3.1 控制把局部高频层加入已有解时产生的误差。摘录所给结论表明，在适当有界性假设及大尺度参数 $\lambda$ 下，密度的 $$H^3$$ 误差和涡量的 $H^2\cap\dot H^{-1}$ 误差为 $O(t\lambda^{-1})$ ，相应时间导数误差为 $O(\lambda^{-1})$ 。这些是中间近似的误差阶，不应误读为最终范数爆炸的速率。完整定理的迭代参数、最终增长率以及各层误差如何统一求和，无法从所提供的第 1 至 10 页摘录中确认。

作者在备注 1.4 中说明，预期构造的解满足对任意 $\delta\in(0,1)$ ， $\rho+x_2\in H^{2-\delta}$ 、 $\omega\in H^{1-\delta}$ ；但论文为简化估计没有追踪这些次临界范数。因此这属于作者说明可依证明步骤补出的预期性质，而不是定理 1.1 中严格陈述并证明的结论。另一个备注指出，全密度取任意小 $$H^2$$ 范数的真空态附近也可作类似构造，初值会更简单；该扩展在摘录中只作说明，未作为定理 1.1 的正式结论展开。

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

证明的实质障碍不是单纯找一个令 $\partial_1u_1(0,t)$ 变大的初值，而是在稳定背景中同时做到：抵消背景稳定化、让初始涡量严格为零、保持密度—涡量反馈所需的形变符号，并将局部尺度上的构造嵌入真实非线性解。速度由涡量经 Biot-Savart 算子非局部决定，新增局部层会影响其他位置和后续输运；因此，线性化模型中的形变增长不能未经误差分析就直接视为非线性方程的结论。

首先，作者通过对称性和局部空穴设计处理稳定背景。密度取偶于 $$x_1$$ 、奇于 $$x_2$$ ，涡量取关于两个坐标分别为奇函数。对称性在解的演化中保持，并推出 $$u(0,t)=0$$ ，所以原点沿流不动，且可将原点速度梯度作为稳定的追踪对象。再选 $\rho_{\mathrm{in}}=x_2$ 于原点附近，使总密度 $-x_2+\rho_{\mathrm{in}}$ 在那里消失。与此同时，参与反馈的密度扰动支集安排在原点之外。由于 $$u(0,t)=0$$ 且速度在短时间内有 $$C^1$$ 控制，流对原点附近距离至多产生指数型变化，故支集与原点的正距离在短时间内保留。局部适定性引理 2.2 提供构造初始阶段所需的高阶正则性控制。

其次，作者把多尺度密度层转换为形变增长。原点处速度梯度可写为

\partial_1u_1[\omega](0,t) =\frac1\pi\int_{\mathbb R^2}\frac{y_1y_2}{|y|^4}\omega(t,y)\,dy,

而涡量方程给出

\partial_t\partial_1u_1[\omega](0,t) =-\frac1\pi\int_{\mathbb R^2}\frac{y_1y_2}{|y|^4} \big(u\cdot\nabla\omega+\partial_1\rho\big)(t,y)\,dy.

在初始时刻 $\omega=0$ ，非线性输运项消失，形变的初始变化由 $\partial_1\rho$ 决定。这使作者能先精确安排密度源项的符号。对满足 $x_1x_2\partial_1f\geq0$ 的环带函数 $$f$$ ，缩放 $$f(2^ix)/(i2^i)$$ 的 $$H^2$$ 代价平方可求和，而其对形变导数的贡献按 $$1/i$$ 累积。利用调和级数发散，作者能在固定的小 $$H^2$$ 预算下把初始形变导数推得足够负。该处的关键不是泛泛的“多尺度方法”，而是 $$H^2$$ 缩放下二阶导数代价与 Biot-Savart 形变泛函对不同尺度的贡献具有不同求和行为。

第三，作者证明负形变在演化中可以持续，并控制由该形变诱发的新涡量。其模型把局部速度近似为 $$k(t)(x_1,-x_2)$$ ，此时流映射有显式表达

\Phi(x,t)=\left(x_1e^{\int_0^t k},x_2e^{-\int_0^t k}\right).

沿流密度保持常数，而涡量满足源项 $-\partial_1\rho$ ，故可显式积分得到涡量。将其代入原点形变公式，并用初始支集的扇形限制 $|x_2|\geq2|x_1|$ 控制形变分母，可避免分母估计因各向异性拉伸而退化。引理 2.4 因而给出负号保持以及形变与初始源强、时间长度之间的定量比较。引理 2.5 进一步保证较早时间内已经积累的形变积分，不会相对于全时间区间的形变积分小得不可用；这对逐层迭代至关重要。

第四，作者需要把显式线性模型嵌入非线性方程。引理 3.1 的做法是并行构造内外解：内层在局部尺度上看到已有流的线性形变 $$k(t)(x_1,-x_2)$$ ，外层则由已有背景及内层产生的速度驱动；此外还构造相应的真实耦合系统解。精确解是内外分量的耦合演化，而“旧解加上线性近似内层”仅是形式近似。引理通过高阶 Sobolev 控制以及两尺度分离，证明两者之间在 $$H^3$$ （密度）和 $H^2\cap\dot H^{-1}$ （涡量）中有 $O(t\lambda^{-1})$ 误差，时间导数误差为 $O(\lambda^{-1})$ 。这些误差估计的作用是确保真实非线性解的原点形变仍继承近似模型中的符号和增长，而不是把近似解误当作精确解。

第五，整体构造采取迭代累积形变。若第 $$j$$ 步已有

\int_0^T\partial_1u_1[\omega_j](0,t)\,dt=-M_j,

则在更小尺度加入新密度扰动，目标是使下一步的积分进一步变负，同时保存此前已累积的贡献。每一步都需要：选择新层的支集和符号；利用线性模型估算其自身反馈；用内外粘合误差保证新层与旧解相互作用可控；选择尺度参数使误差小于新增形变。迭代极限产生速度梯度在原点不受界的解，由此导致高阶 Sobolev 正则性在每个正时间失效。摘录只展示了迭代目标、若干稳定性引理及粘合引理的开头和部分结论，并未展示第 4 节的完整收敛论证，故不能据此补写具体尺度递推式或最终误差求和条件。

最后，作者需要区分“高阶范数爆炸”与“低阶强解存在及唯一”。局部适定性估计以扰动的 $H^{s+1}\times(H^s\cap\dot H^{-1})$ 范数控制密度与涡量，并通过速度梯度的 Sobolev 控制建立短时解；最终则在 $H^1\times L^2$ 的能量解类中证明唯一性。换言之，构造通过把失控放在高阶导数上，避免了把解的存在性或能量类唯一性一并破坏。根据摘录，作者还指出解除原点外近于经典解，但原点速度梯度的极限失控妨碍全局经典正则性。整体创新是将可精确计算的密度驱动涡量反馈模型、多尺度 $$H^2$$ 低代价形变、以及内外近似的定量粘合结合起来，而不是仅依赖一般性的能量估计或紧性论证。

4. 与之前工作的比较

与 Boussinesq 方程既有失适定性结果相比，本文的主要推进在正则性框架和失控机制两方面。作者指出，Elgindi 与 Masmoudi 的相关结果处于 $L^\infty$ 型临界空间，并与 Riesz 变换在 $L^\infty$ 上无界有关；后续工作 [4] 将这一方向推进到强失适定性或范数膨胀，并通过受非线性效应启发的主导阶模型分析。本文则在扰动 $$H^2$$ 范数这一能量临界尺度上证明任意正时间的 $H^2\times H^1$ 爆炸，同时保留能量类唯一性。就所给摘录而言，作者将其称为该系统在这一能量临界空间中的首个强失适定结果。摘录没有给出上述文献的完整定理假设和范数定义，因而无法逐项比较它们的精确参数范围或定量范数膨胀速度。

与 Bourgain-Li 的 Euler 型构造相比，本文不是把二维 Euler 的速度形变机制直接嵌入 Boussinesq 系统。作者明确指出，若只需展示一般 Boussinesq 系统附近的 $$H^2$$ 失适定性，可以尝试令初始密度为零并利用 Euler 构造；但这不能直接解决稳定背景问题，因为稳定分层会抑制速度形变。即使能够改造 Euler 构造，所得不稳定性也会继承 Euler 动力学，而不是揭示稳定 Boussinesq 系统自身的密度—涡量作用。本文通过令 $\omega(0)=0$ ，排除了初始速度涡量所触发的 Euler 型形变；先由密度梯度生成涡量，再由涡量反馈生成形变，故不稳定性的因果链条不同。摘要将相关 Euler 失适定性表述为 $$H^1$$ 机制，而正文提到 Bourgain-Li 的 $$H^2$$ Euler 构造；两者对应的具体比较口径在摘录中并未进一步澄清，因此宜只作上述机制层面的区分。

与作者及 Córdoba 合作者此前的稳定 IPM 方程工作 [3] 相比，本文继承了“在原点附近设计局部空穴以中和稳定化效应”的初值思想，并沿用由局部形变导致 Sobolev 失控的总体策略。但 IPM 的演化结构不同，本文额外面对 Boussinesq 的密度—涡量耦合：密度生成涡量，涡量通过 Biot-Savart 作用于速度，再反过来改变密度输运。因此速度形变的时间演化不再具有 IPM 情形中可直接利用的单调性，作者需要另行证明负形变持续性、早期形变积分的有效性，并控制内外尺度间的非线性误差。作者还指出一个解正则性的差别：其构造序列在原点的 $\partial_1u_1$ 发散，所以极限解速度梯度不有界，最终得到的是强解而非 IPM 工作中的经典解。这里的比较反映了机制和正则性代价的变化，而不是声称本文的所有证明步骤都无法从 IPM 工作借鉴。

与一般稳定背景下的正则性研究相比，本文不解决光滑有限能量初值是否全局正则，也不构造总密度本身有限能量的解。定理的初始总密度为 $-x_2+\rho_{\mathrm{in}}$ ，由于背景 $$-x_2$$ 在整个平面上不属于有限能量空间，作者明确说明总密度具有无限能量；所控制并任意缩小的是扰动 $\rho_{\mathrm{in}}$ 的 $$H^2$$ 范数。相应地，爆炸陈述写作 $\|\rho(t)+x_2\|_{H^2}=\infty$ ，而非总密度的 $$H^2$$ 范数从有限值变为无限值。作者备注称，若改在真空态附近并要求全密度初值具有任意小的 $$H^2$$ 范数，类似机制仍可构造，且无需空穴引理 2.1；但这一版本在摘录中只作说明，未以主定理形式完整陈述。

本文的时间结论也应准确理解。与只证明某个序列时间上的范数变大或短时间范数膨胀不同，定理对固定构造所得的解断言所有 $t\in(0,T_\epsilon]$ 均有 $H^2\times H^1$ 范数无穷大；然而它只在一个依赖于初值的正时间区间内给出构造。摘录没有提供 $T_\epsilon$ 关于 $\epsilon$ 的显式尺度，也没有给出范数从有限值增长到某个大但有限阈值的速率，因为定理结论直接是每个正时间属于性失效。不能据此推断任意给定正时间都存在相应构造，或推断 $T_\epsilon$ 可统一取正下界。

代价和局限主要体现在解的正则性、初值背景以及证明适用范围。初值密度扰动是光滑且可小于任意给定 $$H^2$$ 阈值的，但极限解在每个正时间不属于临界高阶空间；它仍属于 $H^1\times L^2$ 能量类并具有定理所列的时间连续性和唯一性。构造需要对称性、局部空穴以及特殊符号的多尺度数据，不是对一般小扰动的定量增长定理。作者备注预期解在 $H^{2-\delta}\times H^{1-\delta}$ 中保持次临界正则性，但没有在主结果中追踪证明，故这不能作为已正式陈述的定理结论。论文还提出该密度—涡量反馈或可用于三维轴对称带旋 Euler 方程，但这属于未来方向；所给内容没有相应的三维定理或证明。总体而言，本文把稳定 Boussinesq 的临界失适定推进到强烈的瞬时 Sobolev 正则性丧失，同时以受限的特殊初值构造、非经典极限正则性及有限时间区间为代价；更广泛初值稳定性、爆炸时间的定量下界以及三维推广，均无法从所提供内容中确认。

6. 偶-奇扰动下广义 Boussinesq 方程孤立波附近的全局动力学

Global dynamics near solitary waves for the generalized Boussinesq equation under even-odd perturbations

作者: Xiaoguang Li, Jun Wang

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39356v1

摘要翻译: 我们研究广义 Boussinesq 方程在驻波孤立波附近的全局动力学。对于偶-奇扰动，Maulen [J. Math. Pures Appl. 177 (2023)] 在孤立波 Q 附近构造了一个渐近稳定的中心-稳定流形 M。我们通过证明：任何位于 M 之外、且充分小的 Q 的偶-奇扰动，要么在 t -> +infinity 时于能量空间中散射，要么在有限时间内爆破，从而补全这一局部描述。证明依赖于两个主要要素。首先，我们通过识别两个不变区域，对能量低于基态能量的偶-奇解进行分类，这两个区域分别对应散射和爆破。其次，我们在一维情形下建立单次通过定理：非散射解一旦离开适当选取的孤立波邻域，就不能再次进入该邻域。因此，解始终局限于两个不变区域之一。结合上述分类，我们得到 M 之外解的散射与爆破二分结论。

Comment: 无

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

1. 研究对象和背景

本文研究一维聚焦广义 Boussinesq 方程在静止孤立波附近的正向时间动力学. 标量方程及其一阶 Hamilton 系统形式为

\partial_t^2u-\partial_x^2u+\partial_x^4u+\partial_x^2(|u|^pu)=0,

\partial_tu=\partial_xv,\qquad \partial_tv=\partial_x(-\partial_x^2u+u-|u|^pu),\qquad (u,v)|_{t=0}=(u_0,v_0).

状态变量为 $\vec u=(u,v)^T$ , 自然能量空间是 $X_1=H^1(\mathbb R)\times L^2(\mathbb R)$ . 系统守恒能量和另一守恒量分别为

E(u,v)=\frac12\int_{\mathbb R}\big(v^2+u^2+|\partial_xu|^2\big)\,dx-\frac{1}{p+2}\int_{\mathbb R}|u|^{p+2}\,dx, \qquad I(u,v)=\int_{\mathbb R}uv\,dx.

方程存在行波族 $\vec Q_c(x-ct)=(Q_c(x-ct),-cQ_c(x-ct))^T$ , 其中 $$c^2<1$$ 且

-\partial_x^2Q_c+(1-c^2)Q_c-Q_c^{p+1}=0.

本文聚焦 $$c=0$$ 的静止基态 $\vec Q=(Q,0)^T$ , 即 $-Q''+Q-Q^{p+1}=0$ , 并限制在偶奇对称子空间 $X_1^{eo}=H^1_e(\mathbb R)\times L^2_o(\mathbb R)$ . 这是与反射对称相容的演化子空间, 但该对称性并不自动保证紧性: 一对互为镜像、彼此向空间无穷远移动的局域波包仍可能保持整体偶奇对称.

这一问题处于孤立波稳定性与聚焦色散方程阈值动力学的交界. 已有局部理论给出基态附近的中心稳定流形, 其上解可长期保持靠近 $\vec Q$ ; 但这并未说明离开流形的扰动最终会发生什么. 另一方面, 低于基态能量时常见的势阱分解以变分泛函的符号区分散射与爆破, 但将这种次临界能量结论用于基态附近的超临界能量轨道, 还需要控制轨道离开孤立波邻域后的行为. 本文的目标是补足这一环节: 对充分小的偶奇扰动, 证明中心稳定流形之外不存在其他长期滞留或反复返回的动力学, 并将离开后的演化归入散射或有限时爆破两类.

作者采用的关键势阱泛函为

K_1(\vec u)=\int_{\mathbb R}\big(|\partial_xu|^2+|u|^2+|v|^2-|u|^{p+2}\big)\,dx.

它与能量约束下的振幅变分相关, 其正负号在相应能量区域内区分不同动力学通道. 本文的主要定理要求 $$p>6$$ . 这一限制不只是孤立波局部谱分析的要求; 它还用于能量空间上的色散与时空估计、非线性插值以及严格散射阈值的建立. 因而结果不是对所有幂次的完整分类, 而是在该参数范围和偶奇对称设定下完成基态附近的正向整体分类.

2. 主要定理或主要结果

主要结论首先是定理 1.5 的基态能量以下分类. 设 $$p>6$$ , 并令

\mathcal E_{<Q}=\{\vec u\in X_1^{eo}:E(\vec u)<E(\vec Q)\},\qquad K_+=\{\vec u\in\mathcal E_{<Q}:K_1(\vec u)>0\},\quad K_- =\{\vec u\in\mathcal E_{<Q}:K_1(\vec u)<0\}.

则 $$K_+$$ 与 $$K_-$$ 在流下不变. 若初值属于 $$K_+$$ , 解对所有正负时间全局存在, 并在 $t\to\pm\infty$ 时于能量空间散射, 即分别趋近线性流 $e^{tJA}\vec u_\pm$ ; 若初值属于 $$K_-$$ , 解在正、负两个时间方向均有限时爆破. 这里能量空间为 $H^1\times L^2$ , 线性群由 $$J A$$ 生成.

定理 1.2 给出孤立波附近的正向整体分类. 对每个 $$p>6$$ , 存在充分小的 $\varepsilon>0$ , 使得任意 $\vec u_0\in U_\varepsilon(\vec Q)\cap X_1^{eo}$ 对应的解满足下列二分. 若 $\vec u_0\in\mathcal M$ , 其中 $\mathcal M$ 是 Maulén 构造的中心稳定流形, 则解正向全局存在且一致保持近邻：

\|\vec u(t)-\vec Q\|_{X_1}\le C\|\vec u_0-\vec Q\|_{X_1},\qquad t\ge0.

若 $\vec u_0\notin\mathcal M$ , 则解或者正向散射, 即存在 $\vec u_+\in X_1$ 使

\lim_{t\to+\infty}\|\vec u(t)-e^{tJA}\vec u_+\|_{X_1}=0,

或者在某个有限正时间 $$T$$ 爆破, 且 $\lim_{t\to T^-}\|\vec u(t)\|_{X_1}=+\infty$ . 这是正向结论; 定理 1.2 并未声称中心稳定流形外的近孤立波数据也具有相同的负向二分. 流形上局部趋近孤立波的结论来自 Maulén 的先前结果, 本文引用并结合该局部描述.

为处理不能仅凭势阱分类的基态附近超临界能量轨道, 论文还建立了单次穿越定理. 在合适选取的尺度 $0<\varepsilon\ll R\ll1$ 下, 非散射解离开 $\pm\vec Q$ 的指定邻域后不能再次进入该邻域; 离开后其轨道保持在由 $$K_1$$ 符号确定的正、负通道之一. 因而后续分别可调用正通道的散射机制和负通道的爆破机制. 这一结论是由定理 1.5 推到定理 1.2 的关键桥梁.

另一个定量结果是命题 3.4 的一般数据严格散射阈值. 若 $E_c^{gen}$ 表示一般能量空间数据在 $E<E_c^{gen}$ 、 $$K_1>0$$ 条件下具有一致散射范数界的阈值, 则存在与 $$p>6$$ 无关的 $\varepsilon_*>0$ , 使

E_c^{gen}\ge\left(\frac12+\varepsilon_*\right)E(\vec Q)>\frac12E(\vec Q).

该结论不是一般数据低于完整基态能量的散射定理, 而是严格超过半基态能量的统一阈值改进. 它用于排除偶奇对称剖分中可能逃逸到无穷远的一对反射分量, 也用于单次穿越论证中可能恰有半基态能量的逃逸分量. 作者未在所给内容中给出 $\varepsilon_*$ 的显式数值.

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

证明的第一项困难是次临界散射中的紧性失败. 对偶函数作线性剖分时, 除了有固定中心的剖面, 还可能出现空间中心趋于无穷的镜像剖面对. 一维偶对称并不像高维径向情形那样由 Strauss 型尾部估计直接压制空间逃逸; 本方程又没有 Galilean 或 Lorentz 变换, 无法通过改变参考系消除移动中心. 作者因此没有把偶奇对称误当作紧性, 而是先建立一个适用于一般数据的严格阈值, 再利用该阈值处理逃逸对.

这个阈值改进的核心是构造对线性流具有单调性结构的有界色散泛函. 令 $D=|\partial_x|$ , $\Lambda=(1+D^2)^{1/2}$ , $$H$$ 为 Hilbert 变换, 并作实线性变换

z=u-i\Lambda^{-1}Hv,\qquad \|z\|_{H^1}=\|\vec u\|_{X_1}.

在该变量下系统对角化为

i\partial_tz=D\Lambda z-D\Lambda^{-1}(|u|^pu).

记自由群 $S_0(t)=e^{-itD\Lambda}$ , 并设 $\mathcal D(z)=\|\partial_x|z|^2\|_{L^2}^2$ . 作者通过 Fourier 相位变量变换得到

\int_{\mathbb R}\mathcal D(S_0(s)z)\,ds\le\frac23\|\Lambda^{1/2}z\|_{L^2}^4.

于是可定义有界泛函 $\mathcal M(z)=\frac12\int\operatorname{sgn}(s)\mathcal D(S_0(s)z)\,ds$ 及其归一化 $\mathcal V(z)=\mathcal M(z)/\|\Lambda^{1/2}z\|_{L^2}^4$ . 对非线性项 $D\Lambda^{-1/2}(|u|^pu)$ , 作者分离其沿 $\Lambda^{1/2}z$ 的复数内积平行分量与正交分量. 平行分量恰与归一化分母求导产生的项抵消, 从而导数化为

\partial_t\mathcal V(z(t))=\frac{-\mathcal D(z(t))+\mathcal R_\perp(t)}{\|\Lambda^{1/2}z(t)\|_{L^2}^4}.

复正交性消去平方展开中的交叉项; 再结合密度估计与分数阶链式估计, 得到 $|\mathcal R_\perp|\le a_p\|z\|_{H^1}^p\mathcal D(z)$ . 能量陷获在半基态能量处提供严格吸收余量, 并通过关于 $$p$$ 的一致连续性将余量延伸到 $E<(1/2+\varepsilon_*)E(\vec Q)$ , 即 $|\mathcal R_\perp|\le(1-\delta)\mathcal D(z)$ , 其中 $\delta>0$ 可对 $$p>6$$ 一致选取. 积分单调性控制 $L_t^6L_x^\infty$ 范数, 再结合能量界获得非线性项 $JN(\vec u)\in L_t^1X_1$ , 最终推出双向散射.

偶奇次临界散射部分采用剖分、临界元与刚性论证, 但关键在于如何使用上述一般数据阈值. 若偶奇散射阈值低于 $E(\vec Q)$ , 逃逸镜像对中的每个组成部分能量至多为偶奇阈值的一半, 因而严格落在一般数据的改进阈值之下, 各自散射. 这排除了对称逃逸对作为临界障碍. 剩余临界剖面可固定在原点, 其轨道在 $$X_1$$ 中预紧; 正势阱及紧性给出 $$K_2$$ 的一致正下界, 局部 virial 恒等式遂与预紧性矛盾. 负势阱的有限时爆破则由正则化凹性论证得到.

局部单次穿越分析利用线性化算子 $$JL$$ 在偶奇子空间中的一对实特征值 $\pm\mu$ 及正交余项的能量相干性. 写 $\vec u=\vec Q+\vec\eta$ , 分解为稳定、不稳定模态系数和辐射部分 $\vec\gamma$ , 其中 $\langle L\vec\gamma,\vec\gamma\rangle\gtrsim\|\vec\gamma\|_{X_1}^2$ . 能量展开控制一个与到 $\pm\vec Q$ 距离等价的量 $$d_Q$$ , 并推出解离开小邻域后以指数速度到达较大邻域. 振幅和空间缩放泛函 $$K_1,K_2$$ 在满足 $|E(\vec u)-E(\vec Q)|\ll d_Q(\vec u)^2$ 时具有相同符号, 因而可确定离开后的通道.

单次穿越证明按泛函符号分别构造 virial 控制. 当 $$K_1<0$$ 时, 使用带低频截断及时间端点修正的双线性泛函 $$T_R(u,v)$$ , 建立 $\partial_tT_R\ge c[-K_1]-C\sqrt R$ ; 该增长与假设返回时由指数离开带来的大小矛盾. 当 $$K_1>0$$ 时, 使用 $I_\phi=\int\phi uv\,dx$ 的局部 virial 恒等式, 其主项为 $$-K_2$$ . 正通道缺少直接的径向尾部控制, 作者在假设首次返回的有限时间区间上使用紧性论证控制截断误差; 末段靠近孤立波的时间长度趋于无穷, 与非散射解的时空范数行为冲突. 两种符号分别得到矛盾, 从而排除返回. 最后, 对离开后的正通道通过次临界散射与扰动论证完成散射分类, 对负通道通过正则化凹性推出有限时爆破.

4. 与之前工作的比较

本文的直接局部参照是 Maulén [29]. 该工作在 $$p>2$$ 时构造了静止孤立波附近的偶奇中心稳定流形 $\mathcal M$ , 并说明流形上的解正向全局、保持邻近且局部趋于孤立波; 若一个全局解始终留在足够小邻域, 则它必须属于 $\mathcal M$ . 这些结论识别了不离开邻域的解, 但并未分类流形外解离开之后的命运. 本文在更强的幂次假设 $$p>6$$ 下补上这一点: 小邻域内但不在 $\mathcal M$ 上的偶奇扰动正向散射或有限时爆破. 因而推进的是局部流形描述到正向全局二分, 不是对 Maulén 的局部流形构造本身作替代.

在次临界能量动力学方面, 本文的思路承接 Kenig–Merle [15,16] 的浓缩紧性与刚性框架, 以及其后在 NLS 等模型中的发展. 定理 1.5 给出与聚焦问题势阱结构相应的低于基态能量分类: $E<E(\vec Q)$ 时, $$K_1$$ 的正号导向双向散射, 负号导向双向有限时爆破. 文中将其与 Killip–Visan [17]、Holmer–Roudenko [14] 等阈值结果联系起来, 但本方程是一维广义 Boussinesq 系统, 不能直接照搬具有 Galilean 或 Lorentz 不变性的模型中的移动剖面处理.

相对于高维径向框架, 本文面临的具体差异是: 一维偶奇对称仍允许两份反射剖面向两端逃逸, 而 Nakanishi–Schlag [32,33] 的一遍通过方法所用的三维径向衰减估计在这里不可用. 作者并非声称偶奇对称本身恢复了紧性, 而是先证明一般数据散射阈值满足 $E_c^{gen}\ge(1/2+\varepsilon_*)E(\vec Q)$ , 再借此排除逃逸对. 该严格超过半基态能量的间隙也是处理单次穿越误差时可能出现的半能量逃逸分量所必需的. 因此这里的贡献是针对一维对称逃逸构造补偿机制, 而不是将高维径向论证原样移植.

作者还专门指出 Chen–Shao [6] 的相关线性剖分没有纳入对称逃逸对, 并认为其将所有空间平移参数置零的断言未充分处理一维偶奇情形. 本文据此区分一般数据阈值与偶奇阈值, 以严格阈值改进支撑相应的偶奇散射证明. 这一讨论针对的是所提供文本明确指出的技术缺口; 对 Chen–Shao 工作其他结论的有效性, 从提供内容中无法确认, 不宜作更广泛判断.

与已有结果相比, 本文也清楚保留了适用范围和代价. 第一, 基态以下的分类及基态附近的整体二分均要求 $$p>6$$ , 而中心稳定流形先前可在 $$p>2$$ 构造; 因而全局分类的幂次范围更窄. 第二, 散射/爆破二分只针对偶奇扰动; 作者明确说明一般非对称数据低于完整基态能量的正势阱散射尚未解决, 因为仍需处理带移动空间中心的临界元. 第三, 阈值改进仅给出一个与 $$p>6$$ 无关的正间隙存在性, 并不等同于证明所有一般数据在 $E<E(\vec Q)$ 、 $$K_1>0$$ 时散射. 第四, 定理 1.2 是正向时间分类, 流形外的负向动力学并未由该定理分类; 双向结论只在定理 1.5 的严格次基态能量区域成立. 这些限制与本文证明中对对称性、能量阈值和单次穿越机制的依赖相对应.

7. 玻尔兹曼方程黏性 BGK 模型的流体力学极限

Hydrodynamic Limit of a Viscous BGK Model of the Boltzmann Equation

作者: Gabriella Puppo, Thomas Rey, Tommaso Tenna

类别: math.AP, math-ph

发布日期: 2026-09-30

arXiv ID: 2609.39283v1

摘要翻译: 本文研究玻尔兹曼方程 BGK 模型的一种黏性正则化及其流体力学极限。首先，我们利用由空间 Fisher 信息导出的界，严格证明了向黏性 Euler 方程的收敛。随后，我们分析了单组分和多组分情形下带熵不等式的 Euler 方程的形式推导，并将我们的结果与对称双曲热力学相容理论进行比较。

Comment: 无

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

1. 研究对象和背景

本文研究带空间扩散的 BGK 动力学方程在碰撞频率趋于无穷时的流体力学极限。未知量是非负粒子分布函数 $f^\varepsilon(t,x,v)$ ，其中 $t\geq 0$ ， $x,v\in\mathbb{R}^d$ ， $\varepsilon>0$ 为 Knudsen 数。方程为

\partial_t f^\varepsilon+v\cdot\nabla_x f^\varepsilon =\frac{1}{\varepsilon}\bigl(M[f^\varepsilon]-f^\varepsilon\bigr)+\tau\Delta_x f^\varepsilon, \qquad f^\varepsilon(0,x,v)=f_0^\varepsilon(x,v),

其中 $\tau>0$ 固定， $M[f^\varepsilon]$ 是与 $f^\varepsilon$ 具有相同质量、动量和能量矩的局部 Maxwell 分布：

M[f](t,x,v)=\frac{\rho}{(2\pi T)^{d/2}} \exp\!\left(-\frac{|v-u|^2}{2T}\right),

\rho=\int_{\mathbb{R}^d}f\,dv,\qquad \rho u=\int_{\mathbb{R}^d}vf\,dv,\qquad d\rho T=\int_{\mathbb{R}^d}|v-u|^2f\,dv.

BGK 松弛项 $\varepsilon^{-1}(M[f]-f)$ 保持质量、动量和动能；当 $\varepsilon\ll1$ 时，碰撞使分布趋向局部 Maxwell 分布。本文额外加入的 $\tau\Delta_x f$ 是作用于位置变量而非速度变量的人工空间黏性。它对宏观守恒量产生同样的扩散项，并通过熵耗散提供空间正则性。作者将这种空间扩散解释为一种模型正则化；文中也指出，它与速度扩散所描述的热浴机制不同。

由碰撞不变量得到的宏观变量为

U^\varepsilon=(\rho^\varepsilon,\rho^\varepsilon u^\varepsilon, \rho^\varepsilon E^\varepsilon)^T, \qquad E^\varepsilon=\frac12|u^\varepsilon|^2+\frac d2T^\varepsilon,

其守恒律具有形式 $\partial_tU^\varepsilon+\nabla_x\cdot F^\varepsilon=\tau\Delta_xU^\varepsilon$ 。预期的极限系统是可压缩 Euler 方程的黏性正则化：

\begin{aligned} \partial_t\rho+\nabla_x\cdot(\rho u)&=\tau\Delta_x\rho,\\ \partial_t(\rho u)+\nabla_x\cdot(\rho u\otimes u+pI)&=\tau\Delta_x(\rho u),\\ \partial_t(\rho E)+\nabla_x\cdot\bigl(u(\rho E+p)\bigr)&=\tau\Delta_x(\rho E), \end{aligned} \qquad p=\rho T.

这里的目标不是直接证明 $\tau\to0$ 的无黏极限，而是在固定 $\tau>0$ 时令碰撞尺度 $\varepsilon\to0$ ，并严格识别上述宏观黏性系统。问题的重要性在于：经典 BGK 方程的分布函数缺乏空间紧性，宏观矩通常只能借助速度平均引理获得分数阶正则性；而人工空间扩散可能提供更直接的矩紧性，从而在不依赖重整化解框架的情形下识别流体极限。

论文还提出第二条线索：利用 Chapman-Enskog 展开形式讨论带熵不等式的 Euler 极限，并将其与对称双曲、热力学相容理论（SHTC）比较。单组分和多组分情形均在摘要及第三节开头被提及。其动机是说明宏观熵产生可以由动力学模型的结构导出，而不必像 SHTC 构造中那样额外指定熵产生以保证总能量守恒。不过，所提供的 PDF 内容只到第三节开头，未包含该形式推导的具体方程、假设和结论，因而无法进一步确认其精确模型与定量结果。

2. 主要定理或主要结果

严格证明的核心结果是固定 $\tau>0$ 时， $\varepsilon\to0$ 的黏性流体力学极限。定理 2.6 假设非负初值满足

\int_{\mathbb{R}^d\times\mathbb{R}^d} f_0^\varepsilon(x,v)\bigl(1+|x|^2+|v|^q+|\log f_0^\varepsilon(x,v)|\bigr)\,dx\,dv<\infty, \qquad q\leq4,

并取满足文中相应一致有界条件的初值族。则存在子列使得

(\rho^\varepsilon,\rho^\varepsilon u^\varepsilon, \rho^\varepsilon E^\varepsilon) \longrightarrow (\rho,\rho u,\rho E) \quad\text{强收敛于 }L^1_{\mathrm{loc}}((0,T)\times\mathbb{R}^d).

在 $\{\rho>0\}$ 上， $u^\varepsilon\to u$ 、 $T^\varepsilon\to T$ 几乎处处；极限 $(\rho,u,T)$ 是上述黏性 Euler 系统的弱解。该定理给出的是子列紧性及弱解识别，并未在所提供内容中声称极限解唯一，故不能据此将子列收敛加强为整个族收敛。定理针对固定正黏性系数 $\tau$ ，不涵盖 $\tau\to0$ 的进一步极限。

证明所依赖的熵估计是命题 2.1。对足够正则的非负解，Boltzmann 熵

H(f)=\int_{\mathbb{R}^d\times\mathbb{R}^d}f\log f\,dx\,dv

满足

H(f^\varepsilon(t))+ \frac1\varepsilon\int_0^t\!\int_{\mathbb{R}^d} D_{\mathrm{BGK}}(s,x)\,dx\,ds+ 4\tau\int_0^t\!\int_{\mathbb{R}^d\times\mathbb{R}^d} |\nabla_x\sqrt{f^\varepsilon}|^2\,dx\,dv\,ds \leq H(f_0^\varepsilon),

其中

D_{\mathrm{BGK}}(t,x)=\int_{\mathbb{R}^d} (f^\varepsilon-M[f^\varepsilon]) \log\!\left(\frac{f^\varepsilon}{M[f^\varepsilon]}\right)\,dv\geq0.

因此，松弛耗散具有 $\varepsilon^{-1}$ 权重，空间 Fisher 信息具有 $4\tau$ 权重。结合熵下界、质量和总动能守恒，作者得到 Fisher 信息在有限时间区间上的一致积分控制，以及 $f^\varepsilon$ 的局部弱紧性。

命题 2.2 还证明了速度高阶矩的传播：若初值的 $$q$$ 阶速度矩一致有界（ $q\geq2$ ），则对有限 $$T$$ ，有

\sup_{\varepsilon>0}\sup_{0\leq t\leq T} \int f^\varepsilon(t,x,v)|v|^q\,dx\,dv\leq C_T.

其常数依赖初始矩界和 $$T$$ ，但不依赖 $\varepsilon$ 。在极限通量识别中，本文具体使用到四阶速度矩控制。

命题 2.3 将微观 Fisher 信息转移为宏观矩的空间梯度控制。若 $\psi(v)$ 至多按 $$1+|v|^q$$ 增长，且相应 $$2q$$ 阶矩有界，则对 $\mu_\psi^\varepsilon=\int\psi(v)f^\varepsilon\,dv$ 有

|\nabla_x\mu_\psi^\varepsilon(t,x)| \leq2\left(\int |\psi(v)|^2f^\varepsilon\,dv\right)^{1/2} \left(\int|\nabla_x\sqrt{f^\varepsilon}|^2\,dv\right)^{1/2}.

特别地，质量、动量和总能量密度的空间梯度一致有界于 $L^2(0,T;L^1(\mathbb{R}^d))$ 。命题 2.4 随后利用 Aubin-Lions-Simon 紧性框架，推出宏观守恒变量在 $L^1_{\mathrm{loc}}$ 中强相对紧。

命题 2.5 对每个固定 $\varepsilon>0$ 、 $\tau>0$ ，在

f_0^\varepsilon\geq0,\qquad \int f_0^\varepsilon(1+|v|^2+|\log f_0^\varepsilon|)\,dx\,dv<\infty

下构造全局非负弱解，并保留质量、总能量守恒及熵耗散不等式。该存在性结论与定理 2.6 的流体极限结论是不同层次的结果，后者还要求空间二阶矩和最高至四阶的速度矩等假设。

所提供内容还声称第三节将对单组分和多组分模型作带熵产生的形式流体极限分析，并讨论 SHTC 理论；但具体展开、定理或公式未出现在摘录中。因此，这部分只能视为作者在摘要和章节开头说明的形式分析目标，不能据现有文本列出其严格结论。

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

主要技术障碍是 $\varepsilon\to0$ 时既要获得宏观变量的强紧性，又要处理守恒律通量中的非 Maxwell 偏差和高阶速度矩。仅有熵界通常只能给出分布函数的弱紧性；若无空间正则性，不能直接对相空间变量应用标准 Rellich 紧性。本文的关键做法是利用 $\tau\Delta_x f^\varepsilon$ 在熵估计中产生的 Fisher 信息，将微观空间正则性传递到宏观矩，再以此建立宏观强紧性。

第一步是严谨利用熵结构。将方程乘以 $1+\log f^\varepsilon$ 并对相空间积分，输运项由空间衰减消失。扩散项经分部积分成为

\tau\int\frac{|\nabla_x f^\varepsilon|^2}{f^\varepsilon}\,dx\,dv =4\tau\int|\nabla_x\sqrt{f^\varepsilon}|^2\,dx\,dv.

松弛项则利用 $f^\varepsilon$ 与 $M[f^\varepsilon]$ 具有相同碰撞不变量矩，使 $\log M[f^\varepsilon]$ 对碰撞算子的贡献抵消，留下非负耗散 $D_{\mathrm{BGK}}$ 。因此两种耗散同时可用：空间 Fisher 信息负责紧性，BGK 熵耗散负责极限分布接近 Maxwellian。熵本身在无穷测度相空间上并非自动有统一下界；作者以质量、动能和二阶空间矩，通过广义 Young 不等式控制 $\int f|\log f|$ ，从而获得熵下界及一致可积性。

第二步是把 Fisher 信息有效地传递到宏观量。对任意多项式增长的 $\psi(v)$ ，写成

\nabla_x\int\psi(v)f^\varepsilon\,dv =2\int\psi(v)\sqrt{f^\varepsilon}\, \nabla_x\sqrt{f^\varepsilon}\,dv,

再在速度变量中使用 Cauchy-Schwarz。将结果对 $$x$$ 积分时，第一因子由 $$2q$$ 阶速度矩控制，第二因子由 Fisher 信息控制，最终得到宏观矩梯度的 $$L^2_tL^1_x$$ 界。这一步是本文紧性论证的实质性桥梁：扩散只直接正则化 $$f$$ 的 $$x$$ 变量，而速度方向没有相应导数控制；作者不追求相空间的全正则性，而是针对宏观矩使用矩条件与平方根变换。

第三步是时间紧性与宏观闭合。将动力学方程对 $$1,v,|v|^2$$ 积分，得到

\partial_tU^\varepsilon+\nabla_x\cdot F^\varepsilon =\tau\Delta_xU^\varepsilon.

由第三阶速度矩控制，宏观通量在局部可积；由上一节的空间梯度估计，扩散项可写成散度并置于相应的负 Sobolev 空间。于是 $U^\varepsilon$ 在 $L^2_tW^{1,1}_{\mathrm{loc}}$ 有界，时间导数在 $L^1_tW^{-1,1}_{\mathrm{loc}}$ 有界。结合

W^{1,1}_{\mathrm{loc}}\Subset L^1_{\mathrm{loc}} \hookrightarrow W^{-1,1}_{\mathrm{loc}},

应用 Aubin-Lions-Simon 定理，获得 $U^\varepsilon$ 在时空局部 $$L^1$$ 的强相对紧性。空间扩散的作用不只是形式上给宏观方程加黏性，而是提供了获得该强紧性的微观来源。

第四步是识别局部 Maxwell 平衡及极限通量。由熵耗散有

\frac1\varepsilon\int_0^T\!\int D_{\mathrm{BGK}}\,dx\,dt\leq C_0.

使用 Csiszár-Kullback-Pinsker 不等式及质量守恒，得到

\int_0^T\|f^\varepsilon-M[f^\varepsilon]\|_{L^1_{x,v}}^2\,dt \leq 2M_0C_0\varepsilon\longrightarrow0.

因此极限中分布趋于由宏观矩确定的 Maxwellian。对动量通量，将其拆成 Maxwell 部分与偏差部分：

\int v\otimes v f^\varepsilon\,dv =\int v\otimes v M[f^\varepsilon] \,dv +\int v\otimes v(f^\varepsilon-M[f^\varepsilon])\,dv.

偏差项通过 $$L^1$$ 小量与四阶速度矩界的插值估计消失；能量通量中的三阶速度权重也以相同思路处理。Maxwell 部分则由显式矩公式给出

\int v\otimes vM[f^\varepsilon] \,dv =\rho^\varepsilon u^\varepsilon\otimes u^\varepsilon +\rho^\varepsilon T^\varepsilon I,

\int v\frac{|v|^2}{2}M[f^\varepsilon] \,dv =(\rho^\varepsilon E^\varepsilon+\rho^\varepsilon T^\varepsilon)u^\varepsilon.

再结合宏观变量的强收敛，便可在弱形式中识别 Euler 型通量。由此，证明不是单靠“弱收敛到平衡态”，而是依次建立宏观强紧性、松弛偏差消失，并用高阶矩控制完成非线性通量的闭合。

全局存在性部分采用 Perthame 的 BGK 构造思路。作者先将温度下截断为 $T_\delta=\max(T,\delta)$ ，并截断宏观速度以避免 Maxwellian 退化；再利用输运-扩散-衰减算子的正基本解

\Gamma_\tau(t,x,v)=\frac1{(4\pi\tau t)^{d/2}} \exp\!\left(-\frac{|x-vt|^2}{4\tau t}\right)

写出 Duhamel 映射，在 $C([0,T];L^1_2\cap L^2)$ 中作局部压缩映射，并迭代延拓。最后令 $\delta\to0$ ，借助宏观矩紧性识别正则化 Maxwellian 的极限。此处所提供文本给出了框架与步骤，但关于固定点映射的完整定量常数及各函数空间细节未进一步展开，不能补充摘录以外的技术结论。

4. 与之前工作的比较

本文与经典 BGK 流体极限工作的主要差异，是在动力学方程中加入正的空间扩散 $\tau\Delta_x f$ ，并以此取得宏观矩的空间紧性。引言指出，经典 BGK 缺少分布函数的空间正则性，宏观矩的紧性通常依靠速度平均引理，所得正则性一般只有分数阶 Sobolev 控制，相关背景引文为 [24,25]。本文的路线不同：熵耗散给出 $\nabla_x\sqrt f$ 的 Fisher 信息，再借助速度矩界直接控制质量、动量和能量密度的 $$L^2_tL^1_x$$ 空间梯度。作者据此强调该模型无需依靠重整化解来获得所需的紧性。需要注意，这不是对无扩散经典 BGK 方程的同条件推广，而是对加入空间正则化后的方程建立结果；所需的 $\tau>0$ 是关键结构假设。

就极限范围而言，定理 2.6 处理的是 $\varepsilon\to0$ 、 $\tau>0$ 固定的情形，极限是仍带 $\tau\Delta_x$ 的黏性 Euler 系统。它没有证明 $\tau\to0$ ，也没有据此得到无黏 Euler 的熵解或激波选择结果。引言将 DiPerna [17] 关于一维等熵 Euler 黏性近似的工作，以及 Bianchini-Bressan [5] 关于一维严格双曲小总变差系统的消失黏性收敛结果，作为宏观 vanishing-viscosity 理论背景；这些工作的目标涉及黏性趋零及物理解的选择，而本文严格证明的核心是微观 BGK 碰撞极限与固定黏性宏观系统之间的对应。两类极限所取参数不同，不能将本文定理解释成对上述消失黏性结论的替代或直接延伸。

本文结论的适用代价体现在初值矩与熵条件上。定理 2.6 要求有限二阶空间矩、有限熵以及至多四阶的速度矩；高阶矩传播和四阶矩插值用于控制动量、能量通量的非 Maxwell 偏差。存在性命题 2.5 的假设较侧重质量、动能与熵，但极限定理还需要额外的空间衰减及高阶速度矩。结果是局部时空 $$L^1$$ 强收敛的子列结论，加上正密度区域上速度和温度的几乎处处收敛；摘录没有给出极限唯一性，也未证明全序列收敛。文章还指出初始二阶空间矩用于保证空间衰减，使分部积分和守恒量计算成立。

与 SHTC 理论的比较是本文另一项动机。按引言及第三节开头的说明，SHTC 框架把熵产生方程作为变量演化结构的一部分，并通过对熵产生项施加条件来确保总能量守恒；本文则从保持质量、动量和能量的 BGK 碰撞结构出发，希望在流体极限中自动得到与能量守恒相容的熵产生关系。摘要称作者讨论单组分和多组分情形，并比较其形式极限与 SHTC 理论。然而，提供的 PDF 仅包含第三节开头，未展示具体 Chapman-Enskog 展开、熵产生公式、单组分或多组分定理，也未说明与 SHTC 各项结构条件的逐式对应。因此，关于这部分究竟放宽了哪些假设、导出了何种具体熵不等式或适用于哪些多组分参数范围，从提供内容中无法确认；作者的比较细节在所给摘录中未显式展开。

中文标题: 标度临界空间中的 Calogero–Moser 导数非线性薛定谔方程的孤子分解猜想

中文标题: 具有时空依赖材料性质的石墨烯等离激元建模近似方法

摘要翻译: 本文提出了用于模拟具有时空依赖材料性质的石墨烯片上表面等离激元传播的近似计算方法。从 Maxwell 方程出发，将空间和时间调制的 Drude 权重下电流密度的演化表述为偏积分微分方程（PIDE）。在假设 Drude 权重受到小振幅扰动的条件下，本文构造了正则摄动展开，并在适当的 Sobolev 空间中证明了所得逐阶积分方程的可解性与傅里叶可逆性。本文还介绍了一种数值求解器，该求解器结合了 Volterra 积分方程表述与快速傅里叶变换（FFT），并给出了逐点误差估计和计算复杂度分析。此外，本文针对行波型 Drude 权重调制，利用围道积分和 Laplace 反演推导出一种基于变换的解析解。数值模拟展示了多种波动现象，包括行波、驻波以及指数增长的等离激元电流密度。

不相关论文（基本信息）

1. 由一类非线性随机偏微分方程得到的半空间 KPZ 方程

Half-space KPZ equation from a class of nonlinear SPDEs

作者: Kevin Yang

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

发布日期: 2026-09-30

arXiv ID: 2609.40352v1

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

摘要: We derive the half-space KPZ equation as a continuum limit from a general class of Ginzburg-Landau SPDEs under weak nonlinearity scaling. To our knowledge, this is the first half-space universality result for the KPZ equation beyond integrable models (and their perturbations).

2. 有限测量下反演测地线射线变换：稳定性与重建

Inverting the geodesic ray transform with finite measurements: stability and reconstruction

作者: Yernat M. Assylbekov

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

发布日期: 2026-09-30

arXiv ID: 2609.40015v1

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

摘要: We study the inversion of the geodesic ray transform from finitely many local averages of its data. Under suitable geometric assumptions, we prove Lipschitz stability on any fixed finite-dimensional reconstruction space when the measurement partition is sufficiently fine. We develop convergent reconstruction algorithms based on Steepest Gradient Descent and Conjugate Gradients and implement them using piecewise constant finite element spaces. Numerical experiments in 2D and 3D illustrate the influence of geometry and measurement discretization on reconstruction quality.

3. 关于真空初始数据集空间的拓扑

On the topology of the space of vacuum initial data sets

作者: Romain Gicquaud, Jonathan Glöckle

类别: gr-qc, math.AP, math.AT, math.DG

发布日期: 2026-09-30

arXiv ID: 2609.39987v1

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

摘要: We show that the space of vacuum initial data sets on a closed manifold often has many non-trivial homotopy groups. The starting point is a result of the second named author, which constructs non-trivial elements in the homotopy groups of the space of initial data sets satisfying the strict dominant energy condition, together with a later result in joint work with Bernd Ammann showing that these elements often persist when the strictness assumption is dropped. In this work, we use a parametrized version of the conformal method to show that these elements may also be represented by maps into the space of vacuum initial data sets. This requires two results that may be of independent interest: metrics admitting conformal Killing vectors can be removed from the space of metrics without changing its weak homotopy type, and over the remaining metrics the York decomposition can be carried out in families, the TT-tensors forming a trivial Hilbert bundle. To our knowledge, this is the first result on the global topology of this space.

4. 异质介质中多个小陷阱下的捕获时间最小化

Minimizing capture time with many small traps in heterogeneous media

作者: Denis S. Grebenkov, Theodore Kolokolnikov

类别: cond-mat.stat-mech, math.AP, math.PR

发布日期: 2026-09-30

arXiv ID: 2609.39990v1

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

摘要: We study the problem of optimally placing a large number of small absorbing traps to minimize the mean first-passage time (MFPT) of particles diffusing in a heterogeneous medium with space-dependent diffusivity $$D(x)$$ and prescribed initial particle distribution $$ω(x)$$ . In two and three dimensions we identify two distinct regimes, depending on how strongly the trap ensemble depletes the particles. In the weak trapping regime (fewer traps) the optimal trap density $$ρ$$ is the arithmetic average $ρ= \frac12(μ+ω)$ of the normalized inverse diffusivity $μ\propto 1/D$ and the initial particle distribution $$ω$$ . In the strong trapping regime (more traps), it is proportional to their geometric average, $ρ\propto \sqrt{μω}$ . We revisit the classical Green’s function approach and show that it applies to weak trapping only; in 2D this requires the trap size to be exponentially small in the number of traps, so that most applications of interest fall in the strong trapping regime instead. We therefore develop a homogenization approach that captures both regimes. A special choice of initial particle distribution is $ω\propto 1/D$ , which corresponds to the equilibrium particle distribution in the absence of traps according to the Itô interpretation. In this case we show that the optimal trap distribution is the stationary particle distribution itself, regardless of the trapping strength, and this makes the MFPT constant throughout the domain; the same trap distribution also minimizes the worst-case capture time. Finally, in one dimension we find instead $ρ\propto\left(μω\right) ^{1/3}$ , and a similar analysis covers space-dependent drift and thin domains of variable cross-section. Direct numerical optimization confirms the analytical results.

5. 关于边界 Yamabe 问题中紧性的临界维数，I

On critical dimensions for compactness in the boundary Yamabe problem, I

作者: Liuwei Gong, Seunghyeok Kim, Monica Musso, Juncheng Wei

类别: math.DG, math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39945v1

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

摘要: We construct smooth, non-locally-conformally-flat metrics on the closed ball for which the boundary Yamabe equation admits $L^\infty$ -unbounded sequences of positive solutions. The background metric and the prescribed scalar and boundary mean curvatures are fixed along each sequence. For zero scalar curvature and positive boundary mean curvature, such examples exist with umbilic boundary in every dimension $N\ge22$ and with nonumbilic boundary in every $N\ge15$ . For positive scalar curvature and minimal boundary, the corresponding ranges are $N\ge21$ and $N\ge15$ . For every $N\ge9$ , examples with either umbilic or nonumbilic boundary also exist when the scalar curvature is $$N(N-1)$$ and the prescribed boundary mean curvature is below a negative threshold depending on $$N$$ . The construction uses polynomial metric perturbations and corrected bubbles on the half-space, followed by conformal compactification to the ball. We evaluate or estimate the correction solving the linearized Neumann or Robin boundary problem and its contribution to the quadratic term of the reduced energy. We prove that this quadratic term has negative, nondegenerate local extrema with respect to tangential translations and scale in every stated dimension.

6. Hermitian流形上全非线性椭圆方程的稳定性与梯度估计

Stability and gradient estimates for fully nonlinear elliptic equations on hermitian manifolds

作者: Bin Guo, Jian Song

类别: math.AP, math.DG

发布日期: 2026-09-30

arXiv ID: 2609.39726v1

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

摘要: The known proofs of the gradient estimate for complex Hessian equations rely on the second order estimate of Hou-Ma-Wu and on a Liouville theorem of Dinew-Kolodziej. In this paper, we give a direct proof of the gradient estimate for concave fully nonlinear elliptic equations on compact Hermitian manifolds, which uses neither the Hou-Ma-Wu estimate nor the Liouville theorem. Instead, the proof combines a stability estimate with a comparison argument against nearby smooth admissible functions. We also give an independent proof of the complex Hessian estimate on compact Kähler manifolds, based on Dirichlet Green’s functions, under additional assumptions on the operator.

7. 平面各向同性 $$L_p$$ 对偶闵可夫斯基问题非常值正解的确切个数

The exact number of nonconstant positive solutions to the planar isotropic $$L_p$$ dual Minkowski problem

作者: Guanyu Tao

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

发布日期: 2026-09-30

arXiv ID: 2609.39636v1

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

摘要: We determine the exact number, up to rotations, of nonconstant positive $$C^2$$ solutions to the planar isotropic $$L_p$$ dual Minkowski problem for every $(p, q) \in \mathbb{R}^2$ . We obtain a new parametrization of the associated period integral, characterize the regions where the period is strictly increasing or strictly decreasing, and prove that in the remaining nonmonotone region it has a unique nondegenerate maximum. As a consequence, we have a complete classification.

8. 凯勒-里奇流的I型估计 I

Type I estimates for the Kahler-Ricci flow I

作者: Wangjian Jian, Jian Song

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

发布日期: 2026-09-30

arXiv ID: 2609.39620v1

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

摘要: We establish the Type I estimate for the scalar curvature of the Kahler-Ricci flow on compact Kahler manifolds, for arbitrary smooth initial metrics, whenever the flow develops finite time singularities. This estimate is a building block for the analytic minimal model program with Ricci flow.

9. $�

Regularity of Critical Points of Scale-Invariant Geometric Energies for Even-Dimensional Submanifolds of $\mathbb R^m$

作者: Tian Lan

类别: math.AP, math.DG

发布日期: 2026-09-30

arXiv ID: 2609.39606v1

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

摘要: We consider scale-invariant curvature energies for immersions of closed manifolds of even dimension $$n=2h$$ into $\mathbb R^m$ , with principal term $\int_Σ \big|\nabla^{(h-1)} \vec{\mathrm{I\!I}}\big|_g^2\,d\text{vol}_g$ and arbitrary lower-order polynomial extrinsic invariants of the same scaling. Following the four-dimensional approach developed in joint work with Bernard, Martino, and Rivière, we prove that every weak critical immersion in the natural Sobolev class $W^{h+1,2}$ , whose induced metric and its inverse have $L^\infty$ coefficients, is real-analytic in harmonic coordinates. The proof combines geometric conservation laws, additional structural identities, and elliptic estimates with critical Sobolev coefficients to obtain Morrey decay and bootstrap to full regularity.

10. 双稳态脉动波与周期平流：均匀化及波速的尖锐渐近

Bistable pulsating waves with periodic advection: homogenization and sharp speed asymptotics

作者: Weiwei Ding, Linfeng Xu

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39444v1

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

摘要: We study bistable pulsating waves for reaction-diffusion equations with general periodic advection in arbitrary space dimension, allowing the diffusion matrix to be nonsymmetric. Assuming that the homogenized equation admits a traveling wave with nonzero speed in a given direction, we construct moving pulsating waves for all sufficiently small spatial periods $$L$$ and prove their convergence to the homogenized wave as $L\to0^+$ . The existence range and convergence are uniform in the propagation direction when the homogenized speeds never vanish. We also prove uniqueness of the wave speed for arbitrary periods, profile uniqueness for moving waves, and stationary-wave uniqueness in the standing case under continuity of the competing profile. For spatially homogeneous reactions, we derive the expansion $$c_L=c_0+Lc_1+O(L^2)$$ and an explicit formula for $$c_1$$ . Examples with constant diffusion and zero-mean periodic advection show that $$c_1$$ can have either sign, demonstrating that the heterogeneity in advection may accelerate or decelerate propagation relative to the homogenized limit and that the general $$O(L)$$ speed estimate is sharp.

11. 各向异性拟线性椭圆系统的拓扑不变量：一个庞加莱-霍普夫公式

Topological invariants for anisotropic quasilinear elliptic systems: A Poincaré-Hopf formula

作者: Natalino Borgia

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.40145v1

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

摘要: We consider the functional $I_{δ,Ψ_1,Ψ_2}:X \to \mathbb{R}$ defined for any $z=(u,v) \in X$ as \begin{align*} I_{δ,Ψ_1,Ψ_2}(z) & = \int_Ω Ψ_1(\nabla u) , dx + \int_Ω Ψ_2(\nabla v ) , dx & - \int_Ω H(δ,x,u(x),v(x)) ,dx, \end{align*} where $$Ω$$ is a smooth bounded domain of $\mathbb{R}^N$ , $Ψ_1, Ψ_2: \mathbb{R}^N \to \mathbb{R}$ are convex functions satisfying suitable conditions for $$1 < p, q < N$$ , and the nonlinearity $$H$$ is allowed to exhibit both subcritical and critical growth and may depend on a parameter $δ\in I \subseteq \mathbb{R}$ . Here, the functional space setting is given by the product space $X :=W_0^{1,p}(Ω)\times W_0^{1,q}(Ω)$ equipped with the norm $\|z\|= \|u\|_{1,p} + \|v\|_{1,q}$ for any $z=(u,v)\in X$ , where $\| \cdot \|_{1,s}$ is the usual norm in $W^{1,s}_0(Ω)$ . In this paper we prove that $I_{δ,Ψ_1,Ψ_2}'$ is of class $$(S)_+$$ and we infer that each isolated critical point of $I_{δ,Ψ_1,Ψ_2}$ has critical groups of finite type and a Poincaré-Hopf formula holds.

12. 多项式解映射的自然重参数化与连续性

Natural reparameterization and continuity of the solution map for polynomials

作者: Adam Parusiński, Armin Rainer

类别: math.CA, math.AP, math.FA, math.SP

发布日期: 2026-09-30

arXiv ID: 2609.39389v1

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

摘要: The optimal Sobolev regularity of the roots of a smooth curve of monic complex polynomials is $W^{1,q}$ , for $q \in [1,\frac{d}{d-1})$ , where $$d$$ is the degree. This result is stable in the sense that the solution map from $$C^d$$ coefficients to $W^{1,q}$ roots is continuous. We prove that, after a natural Lipschitz reparameterization, the roots are Lipschitz and the map from $$C^d$$ coefficients to reparameterization and reparameterized roots is continuous with respect to the $W^{1,q}$ topology on the target spaces, for all $q \in [1,\infty)$ . The result is based on a convergence and reparameterization theorem for Almgren’s $$Q$$ -valued Sobolev functions.

13. 五维和六维中的平面性估计

The planarity estimate in dimensions five and six

作者: Stephen Lynch

类别: math.DG, math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39330v1

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

摘要: Naff’s planarity estimate is a crucial tool for analysing singularities of quadratically pinched mean curvature flows in high codimension, showing that high-curvature regions become asymptotically codimension one. Naff established the estimate throughout the full Andrews–Baker pinching range in dimensions at least $$7$$ , while in dimensions $$5$$ and $$6$$ a stronger pinching assumption was required. We close this remaining gap, proving the planarity estimate in dimensions $$5$$ and $$6$$ throughout the full Andrews–Baker range. The key new ingredient is a sharp Kato-type inequality relating the twisting of the principal normal to the gradient of the second fundamental form.

14. 图形曲率商方程的内部曲率估计

Interior curvature estimates for graphical curvature quotient equations

作者: Fei Han, Yan Ma

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39276v1

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

摘要: We prove interior curvature estimates for admissible graphical solutions of curvature quotient equations in the range $3\leq k<n$ on the full Gårding cone $$Γ_k$$ . No convexity, semiconvexity, or $Γ_{k+1}$ -admissibility assumption is required. The proof combines a new quantitative compression inequality, a two-surface comparison and doubling argument, and a Pogorelov estimate.

15. 抛物型均匀化中的奇异集与倍增性质

Singular Sets and Doubling Properties in Parabolic Homogenization

作者: Xiujin Chen

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.38952v1

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

摘要: We study the singular sets of real-valued solutions to periodic parabolic equations with rapidly oscillating coefficients. Under a nondegeneracy condition on the corrector matrix, we obtain a local parabolic $$n$$ -dimensional Hausdorff measure bound and a codimension-two Minkowski estimate. The Minkowski estimate holds at every scale, including scales below the oscillation scale $\varepsilon$ , with a constant independent of $\varepsilon$ . The bounds depend on a coarse-scale ratio of space-time $$L^2$$ mass to terminal-slice $$L^2$$ mass. Independently of corrector nondegeneracy, we prove that this ratio controls doubling at all smaller scales uniformly in $\varepsilon$ . The singular-set argument combines summable errors in truncated Gaussian approximations, comparison of caloric polynomials at different time centers, and a stopping cover with microscopic singular-set estimates.

16. 通过热核正则化的康托罗维奇势的稳定性

Stability of Kantorovich potentials via heat kernel regularization

作者: Bang-Xian Han, Zhuo-Nan Zhu

类别: math.AP, math.MG

发布日期: 2026-09-30

arXiv ID: 2609.39213v1

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

摘要: This survey presents heat kernel regularization as a method for quantitative stability of Kantorovich potentials for the quadratic transport cost. We give a complete new heat kernel proof on Heisenberg groups equipped with the Carnot–Carathéodory distance. For source densities bounded above and away from zero on bounded John domains, we obtain $$L^2$$ stability with rate $W_1^{1/2}$ for arbitrary targets in a fixed compact set, including atomic measures. We also establish a new $$L^2$$ stability estimate on finite-dimensional RCD spaces, with the optimal dimension-independent exponent $$1/2$$ .

17. 随机 Landau–Lifshitz–Baryakhtar 方程在 $\mathbb{R}^d$ 中的适定性

Well-posedness of the Stochastic Landau–Lifshitz–Baryakhtar Equation in $\mathbb{R}^d$

作者: Bin Liu, Fan Xu

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.38888v1

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

摘要: This paper studies the Cauchy problem for the stochastic Landau–Lifshitz–Baryakhtar equation on $\mathbb{R}^d$ , $$d=1,2,3$$ , subject to Stratonovich Gaussian perturbations. We establish a low-regularity well-posedness theory for initial data in $\mathbb{L}^2$ , $\mathbb{H}^1$ , and $\mathbb{H}^2$ . For $\mathbb{L}^2$ initial data, we prove the existence and pathwise uniqueness of global very weak solutions in dimensions one and two, while in dimension three we construct martingale very weak solutions on every finite time interval. For arbitrary $\mathbb{H}^1$ initial data, we prove the existence and uniqueness of global pathwise weak solutions in all dimensions $d\leq3$ . Furthermore, for $\mathbb{H}^2$ initial data, we establish the existence and uniqueness of global pathwise strong solutions in the same range of dimensions. The analysis is based on a frequency-truncation approximation scheme combined with stochastic compactness arguments in the whole-space setting. The key step in global continuation is to justify the stochastic effective-field balance at $\mathbb{H}^1$ regularity, where the drift is available only in a negative Sobolev space. Spatial mollification and passage to the limit yield the energy identity for local weak solutions. A coercive modification then gives global energy moments and excludes finite-time blow-up.

18. 临界椭圆性矩条件下非散度型方程的定量均匀化与大尺度正则性

Quantitative homogenization and large-scale regularity for nondivergence-form equations under a critical ellipticity moment

作者: Jizu Huang, Yong Ma

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.38849v1

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

摘要: We prove quantitative homogenization estimates for linear elliptic equations in nondivergence form with stationary, symmetric coefficients, finite range of dependence, and a deterministic upper ellipticity bound. No deterministic positive lower bound is imposed. We assume that the reciprocal of the infimum of the smallest eigenvalue on a unit ball has a finite moment of order $$d$$ , together with the common continuity condition of Armstrong and Smart. Under these assumptions, we obtain algebraic probability bounds for finite-cell errors and for Dirichlet homogenization errors with nonzero sources. We construct quadratic correctors on the whole space, modulo affine functions, and prove first- and second-order large-scale regularity and the corresponding Liouville theorems. We also identify the effective matrix through the invariant density and quantify smooth spatial averages of the density and the weighted coefficients. The proof separates a unit-trace diffusion from its physical clock. A reverse Hölder estimate for the Green function of a stopped coarse process yields the integrability gain needed to control the clock at the critical moment. Applications include weighted gradient convergence and a finite-domain approximation of the effective matrix.

19. 最优传输的超定问题

Overdetermined problem for optimal transportation

作者: Qing Zhao, Feida Jiang

类别: math.AP

发布日期: 2026-09-30

arXiv ID: 2609.39782v1

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

摘要: In this paper, we establish symmetry results for solutions of overdetermined problems arising in optimal transportation. These problems involve the Monge-Ampère equation with a Dirichlet boundary condition $$u=0$$ on $\partial Ω$ and the natural boundary condition $Du(Ω)=Ω^{*}$ . We show that symmetry holds when the target is the unit ball $$B$$ . For more general target domains, including the cases $$Ω^*=Ω$$ and arbitrary $Ω^{*}$ , symmetry is retained under an additional volume constraint on $$Ω$$ and a boundary condition on $$|Du|$$ . Finally, we introduce a curvature-type overdetermined problem for the Monge-Ampère equation and obtain ellipsoidal or spherical symmetry under an additional integral normalization condition. Our proofs rely on distinct techniques in different contexts: optimal transport and convex analysis for some cases, integral identities and isoperimetric inequalities for others, and P-function methods for the remaining cases. As a byproduct, in the $$τ=2$$ case, a new maximum principle for the $$P$$ -function $φ(x)=\sum_{k,l=1}^{n}{\frac{\partial{S_τ(D^2{u})}}{\partial{u_{kl}}}u_{k}u_{l}}-2\binom{n-1}{τ-1}\int_{0}^{u}{f^{\fracτ{n}}(t)\,dt}$ is established for arbitrary positive and nondecreasing $$f$$ , which is of independent interest.

20. 细长杆中格林希尔不稳定性的磁稳定化与滞回切换

Magnetic Stabilization of the Greenhill Instability and Hysteretic Switching in Slender Rods

作者: Vivekanand Dabade, Likhit Ganedi, Mohd Tahseen, Raghavendra Venkatraman

类别: math.AP, cond-mat.mtrl-sci, math.OC, math.SP, physics.app-ph

发布日期: 2026-09-30

arXiv ID: 2609.38676v1

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

摘要: In a bulk ferromagnet, magnetic energy is six orders of magnitude weaker than elastic energy. Slenderness removes this disparity. In a thin ferromagnetic rod, the two energies scale differently with aspect ratio, and an ordinary nickel wire deforms nonlinearly even under modest remote fields. We develop a three-dimensional variational theory of a ferromagnetic Kirchhoff rod under gravity and magnetic loading and study it through bifurcation analysis, computation, and experiment. A uniform field can stabilize a vertical rod against buckling under its own weight, the classical Greenhill instability. We give sharp criteria for clamped and pinned supports in terms of an Airy-type principal eigenvalue, and give a fine description of the character of the bifurcation landscape, which has qualitatively different features than its non-magnetic counterpart. When subject to the field generated by a permanent magnet, we identify a rich bifurcation landscape as the magnet is rotated quasistatically along a circle of radius slightly larger than the rest length of the rod. We prove that instability occurs only through a planar mode and that each snap is a nondegenerate saddle-node. We also locate the two folds that produce hysteresis. Experiments on $50\,μ$ m nickel wires confirm the Greenhill threshold, restabilization at ~ $$30$$ mT for both clamped and pinned supports, as well as the salient features of the hysteresis loop such as switching angles, width, and geometry of folds accurately. Our analysis and experiments offer principled design rules for field-stabilized filaments and remotely switched bistable actuators. Our analysis also offers a well-conditioned numerical method is valid near the switching events, where the direct numerical minimization is strongly ill-conditioned.

21. 非线性扩散问题中对称性诱导分岔点的验证

Validation of Symmetry Induced Bifurcation Points in Nonlinear Diffusion Problems

作者: Maxime Breden, Evelyn Sander, Thomas Wanner

类别: math.DS, math.AP

发布日期: 2026-09-29

arXiv ID: 2609.38629v1

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

摘要: Nonlinear parabolic partial differential equations frequently occur in the modeling of applied phenomena. They often exhibit complicated dynamics which is organized by their associated equilibrium solutions. As a function of the involved model parameters, the equilibria can be detected using a bifurcation analysis, which in turn is centered around a deeper understanding of bifurcation points. In many applied systems establishing these organizing centers rigorously is difficult to achieve using traditional mathematical techniques. In this paper, we propose a theoretical framework for computer-assisted proofs of symmetry-induced transcritical and pitchfork bifurcations in nonlinear problems. This s accomplished by reformulating the existence question via a zero finding problem for an associated extended nonlinear system. The approach is general and allows for a wide variety of underlying symmetries. We demonstrate how these results can be applied in the setting of parabolic systems and higher-order equations, on both one- and two-dimensional domains.

22. 受迫非线性弹性单指向矢曲面的全局解分支

Global Solution Branches of Forced Nonlinearly Elastic Single-Director Surfaces

作者: Timothy J. Healey

类别: math.AP

发布日期: 2026-09-29

arXiv ID: 2609.38498v1

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

摘要: We prove the existence of classical solutions for nonlinearly elastic, single-director surfaces or shells characterized by finite mid-plane strains and curvature. We employ topological methods, yielding solution continua far (in norm) from an unloaded, stress-free reference configuration, i.e., a global implicit function theorem of sorts. We first establish the existence of a local solution path, which is subsequently shown to be part of a global solution branch, viz., a maximal connected, locally compact set, emanating from the unloaded state.

23. 一般二分量反应扩散系统中的转变动力学

Transition Dynamics in General Two-Component Reaction-Diffusion Systems

作者: Taylan Şengül, Burhan Tiryakioglu, Esmanur Yıldız Akıl

类别: math.DS, math.AP

发布日期: 2026-09-29

arXiv ID: 2609.38450v1

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

摘要: This work provides a comprehensive characterization of the first dynamic transition of a general two-component reaction-diffusion system on a bounded interval in one spatial dimension. By rigorously addressing the difficulties arising from the non-linear interactions between the two components, this study significantly extends the scope of previous studies on dynamic transitions of reaction-diffusion systems, which focused on single-equation cases. The analysis is conducted by employing center manifold reduction within the standard classification of dynamic transition theory. A key contribution of this study is that the nonlinear terms are treated to involve arbitrary powers of the unknowns $$u$$ and $$v$$ . Furthermore, the theoretical framework is complemented by the analysis of a representative and simplified system. As a result of this generalization, the dynamics of a broad class of fundamental reaction-diffusion systems can be systematically classified.

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

1. 具有自由表面的可压缩液体的哈达玛局部适定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. $\mathbb{R}^3$ 中两相流模型的低马赫数极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. Boussinesq 系统存在时间的改进#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 广义 SQG 方程的适定性与极限行为#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 密度-涡量反馈下稳定二维 Boussinesq 方程组在 $H^2$ 中的非存在性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 偶-奇扰动下广义 Boussinesq 方程孤立波附近的全局动力学#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 玻尔兹曼方程黏性 BGK 模型的流体力学极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 一般严格凸区域内 Neumann 边界条件下波动方程的色散#

摘要翻译#

2. Kerr扰动中的Teukolsky方程#

摘要翻译#

3. 整个次极端范围内的克尔稳定性#

摘要翻译#

4. 图像分割的基于共识的动力学模型的整体存在性#

摘要翻译#

5. 基于产物通量的反应 Boltzmann 核层析成像#

摘要翻译#

6. 环域与球壳上的显式罗宾格林函数、共振谱与阻抗反演#

摘要翻译#

7. 不满足等波速型假设的带间接一般记忆阻尼的耦合二阶演化方程的稳定性#

摘要翻译#

8. 标度临界空间中的 Calogero–Moser 导数非线性薛定谔方程的孤子分解猜想#

摘要翻译#

9. 归一化线性随机薛定谔方程的全局能量解#

摘要翻译#

10. 从粒子系统推导 Enskog–Enskog 系统#

摘要翻译#

11. 具有时空依赖材料性质的石墨烯等离激元建模近似方法#

摘要翻译#

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

1. 由一类非线性随机偏微分方程得到的半空间 KPZ 方程#

2. 有限测量下反演测地线射线变换：稳定性与重建#

3. 关于真空初始数据集空间的拓扑#

4. 异质介质中多个小陷阱下的捕获时间最小化#

5. 关于边界 Yamabe 问题中紧性的临界维数，I#

6. Hermitian流形上全非线性椭圆方程的稳定性与梯度估计#

7. 平面各向同性 $L_p$ 对偶闵可夫斯基问题非常值正解的确切个数#

8. 凯勒-里奇流的I型估计 I#

9. $�#

10. 双稳态脉动波与周期平流：均匀化及波速的尖锐渐近#

11. 各向异性拟线性椭圆系统的拓扑不变量：一个庞加莱-霍普夫公式#

12. 多项式解映射的自然重参数化与连续性#

13. 五维和六维中的平面性估计#

14. 图形曲率商方程的内部曲率估计#

15. 抛物型均匀化中的奇异集与倍增性质#

16. 通过热核正则化的康托罗维奇势的稳定性#

17. 随机 Landau–Lifshitz–Baryakhtar 方程在 $\mathbb{R}^d$ 中的适定性#

18. 临界椭圆性矩条件下非散度型方程的定量均匀化与大尺度正则性#

19. 最优传输的超定问题#

20. 细长杆中格林希尔不稳定性的磁稳定化与滞回切换#

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

1. 具有自由表面的可压缩液体的哈达玛局部适定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. $\mathbb{R}^3$ 中两相流模型的低马赫数极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. Boussinesq 系统存在时间的改进

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 广义 SQG 方程的适定性与极限行为

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 密度-涡量反馈下稳定二维 Boussinesq 方程组在 $$H^2$$ 中的非存在性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 偶-奇扰动下广义 Boussinesq 方程孤立波附近的全局动力学

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 玻尔兹曼方程黏性 BGK 模型的流体力学极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 一般严格凸区域内 Neumann 边界条件下波动方程的色散

摘要翻译

2. Kerr扰动中的Teukolsky方程

摘要翻译

3. 整个次极端范围内的克尔稳定性

摘要翻译

4. 图像分割的基于共识的动力学模型的整体存在性

摘要翻译

5. 基于产物通量的反应 Boltzmann 核层析成像

摘要翻译

6. 环域与球壳上的显式罗宾格林函数、共振谱与阻抗反演

摘要翻译

7. 不满足等波速型假设的带间接一般记忆阻尼的耦合二阶演化方程的稳定性

摘要翻译

8. 标度临界空间中的 Calogero–Moser 导数非线性薛定谔方程的孤子分解猜想

摘要翻译

9. 归一化线性随机薛定谔方程的全局能量解

摘要翻译

10. 从粒子系统推导 Enskog–Enskog 系统

摘要翻译

11. 具有时空依赖材料性质的石墨烯等离激元建模近似方法

摘要翻译

不相关论文（基本信息）

1. 由一类非线性随机偏微分方程得到的半空间 KPZ 方程

2. 有限测量下反演测地线射线变换：稳定性与重建

3. 关于真空初始数据集空间的拓扑

4. 异质介质中多个小陷阱下的捕获时间最小化

5. 关于边界 Yamabe 问题中紧性的临界维数，I

6. Hermitian流形上全非线性椭圆方程的稳定性与梯度估计

7. 平面各向同性 $$L_p$$ 对偶闵可夫斯基问题非常值正解的确切个数

8. 凯勒-里奇流的I型估计 I

9. $�

10. 双稳态脉动波与周期平流：均匀化及波速的尖锐渐近

11. 各向异性拟线性椭圆系统的拓扑不变量：一个庞加莱-霍普夫公式

12. 多项式解映射的自然重参数化与连续性

13. 五维和六维中的平面性估计

14. 图形曲率商方程的内部曲率估计

15. 抛物型均匀化中的奇异集与倍增性质

16. 通过热核正则化的康托罗维奇势的稳定性

17. 随机 Landau–Lifshitz–Baryakhtar 方程在 $\mathbb{R}^d$ 中的适定性

18. 临界椭圆性矩条件下非散度型方程的定量均匀化与大尺度正则性

19. 最优传输的超定问题

20. 细长杆中格林希尔不稳定性的磁稳定化与滞回切换