Home » 论文总结

2026年10月07日论文分析

共有 9 篇重点关注论文, 14 篇论文可以了解, 34 篇不相关论文

October 7, 2026 · 59 min · 12452 words · P0lyno3ial

Table of Contents

重点关注论文（完整分析）
了解领域论文（摘要翻译）
不相关论文（基本信息）

生成时间: 2026年10月07日 09:39:15

运行状态: 完成

计划处理论文数量: 57

已完成论文数量: 57

未完成论文数量: 0

重点关注论文数量: 9

了解领域论文数量: 14

不相关论文数量: 34

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

1. 单连通有界区域中二维欧拉方程涡片结构持续性的几何效应

Geometric effects on the persistence of patch structures for 2D Euler equations in simply-connected bounded domains

作者: Emeric Roulley

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

发布日期: 2026-10-06

arXiv ID: 2610.08679v1

摘要翻译: 我们首先证明有界单连通平面区域中欧拉方程均匀旋转涡片解的一个刚性结果。具体而言，若一个边界为实解析的有界单连通区域内含有一个非平凡涡片，且该涡片在作匀速刚性旋转时始终严格位于区域内部，则外围区域必为以旋转中心为圆心的圆盘。证明结合了旋转参考系下的表述、调和函数唯一延拓，以及跨越解析边界弧的广义施瓦茨反射原理。这些工具表明，除非外部区域具有旋转不变性，否则速度场必在涡片周围的某个开区域内消失，这与涡片产生的非零环量相矛盾。沿着同样的思路，我们证明只有圆形区域才存在圆形点涡轨道。

Comment: 50 pages

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

1. 研究对象和背景

本文研究二维不可压缩欧拉方程在有界单连通平面区域中的涡斑块动力学, 特别关注环境区域的几何形状如何限制斑块的刚性旋转, 以及在圆盘的小非径向扰动中是否仍可构造准周期斑块。速度-涡量形式为

\partial_t\omega+v\cdot\nabla\omega=0,\qquad v=\nabla^\perp\Psi,\qquad \Delta\Psi=\omega\quad\text{于 }D,

并在边界施加不可穿透条件 $v\cdot n_{\partial D}=0$ , 等价地取流函数边界条件 $\Psi|_{\partial D}=0$ 。对涡斑块, 涡量为特征函数 $\omega(t)=\mathbf 1_{D_t}$ , 因而其演化可由边界 $\partial D_t$ 的轮廓动力学描述。有限能量涡量的全局弱解理论由 Yudovich 理论提供; 但对特定时间结构, 例如刚性旋转或准周期运动, 还需处理边界演化方程中的非线性与频率共振问题。

在有界区域中, 流函数由 Dirichlet Green 函数决定, 可写成

\Psi(t,x)=\frac{1}{2\pi}\int_D G_D(x,y)\omega(t,y)\,dy, \qquad G_D(x,y)=\log|x-y|+K_D(x,y),

其中 $$K_D$$ 是由边界诱导的调和修正项。点涡量的极限也受边界几何支配: 若 $\omega=\Gamma\delta_{\xi(t)}$ , 则点涡位置满足由 Robin 函数 $$R_D$$ 生成的哈密顿系统

\dot\xi(t)=\frac{\Gamma}{4\pi}\nabla^\perp R_D(\xi(t)), \qquad \Delta R_D=4e^{2R_D}.

因此, 即使局部涡斑块方程与平面情形相似, 全局 Green 函数和边界几何仍会改变运动结构。

第一类问题是反问题: 若一个非平凡涡斑块始终严格位于环境区域内部, 并以非零角速度绕固定中心作刚性旋转, 环境区域必须具有什么形状? 在平面或圆盘中, 已有 Kirchhoff 椭圆、Burbea 分支以及圆盘内从 Rankine 涡斑块分岔的多重对称 V 态等刚性旋转解; 圆盘中的环形斑块与点涡圆轨道也有相关构造。相反, 对一般有界区域, 已知可存在稳态斑块, 但非零角速度的刚性旋转是否与非径向边界相容, 是本文要回答的刚性问题。作者指出, 结论依赖实解析边界, 并不覆盖一般光滑边界。

第二类问题是准周期斑块的存在性。此前准周期涡斑块构造包括浅水准地转、广义准地转与 Euler- $\alpha$ 模型中的工作, 常借助方程中的外部参数避开共振; 对二维 Euler 方程本身, 没有此类独立参数, 研究者转而利用平衡态族的几何参数, 例如 Kirchhoff 椭圆长短轴比或圆盘内 Rankine 斑块半径。此前这些 Euler 准周期构造的环境区域为平面或圆盘, 均具有径向对称性。本文试图填补的空白, 是在轻微非径向的有界单连通区域中建立准周期斑块, 并厘清该运动能否仍只是整体刚性旋转。

本文的第二部分以单位圆盘内半径为 $$b$$ 的 Rankine 斑块为基态, 将环境区域的形变幅度设为 $\varepsilon^2$ , 同时将斑块边界扰动设为 $\varepsilon$ 量级。这个匹配使未扰动圆盘问题的线性频谱成为主导阶, 而区域形变作为较小扰动进入 Hamiltonian 轮廓动力学。该构造还要求环境区域关于一个轴对称, 以保留可逆性并实施 KAM 可约化。由此, 本文同时研究了几何对刚性运动的排斥作用, 以及几何微扰下非刚性准周期结构的持续性。

2. 主要定理或主要结果

定理 1.1 给出涡斑块的几何刚性。设 $D\subset\mathbb R^2$ 是有界单连通区域, $\partial D$ 为实解析曲线。若 Euler 方程存在非平凡刚性旋转斑块

\omega(t,x)=\mathbf 1_{R_{\Omega t}D_0}(x),\qquad \Omega\ne0,

其中 $D_0\Subset D$ 是非空单连通开集, $R_\theta$ 表示绕原点旋转角度 $\theta$ , 则 $$D$$ 必须是以原点为中心的圆盘。因而在非圆形实解析环境区域中, 不存在满足这些假设的非平凡非定常 V 态。定理不排除 $\Omega=0$ 的稳态斑块; 文中明确指出非零角速度是证明中利用连续旋转族的关键假设。

定理 1.2 是点涡版本。若 $D\subset\mathbb R^2$ 是有界单连通区域且边界为 $$C^1$$ , 并且区域内存在一个圆形点涡轨道, 则 $$D$$ 必为圆盘。该结论对边界正则性的要求低于斑块刚性定理, 但证明所用的 Robin 函数边界行为需要 $$C^1$$ 条件。摘录没有给出该圆轨道的额外定量参数范围。
定理 1.3 构造非径向区域内的准周期涡斑块。固定 $$0<b_0<b_1<1$$ , 有限模态集 $S\subset\mathbb N^*$ 及振幅 $(a_j)_{j\in S}$ 。存在充分小的 $\varepsilon_0>0$ , 使得当 $0<\varepsilon<\varepsilon_0$ 且环境区域由共形映射

\Phi(x)=x+\varepsilon^2f_\varepsilon(x),\qquad f_\varepsilon\text{ 在 }\mathbb D\text{ 上解析},\qquad f_\varepsilon(0)=0,

给出, 并满足文中要求的实对称条件 $f_\varepsilon(x)=\overline{f_\varepsilon(\bar x)}$ 及 $\varepsilon\mapsto f_\varepsilon$ 连续性时, 存在 Cantor 型集合 $C_\infty^\varepsilon\subset(b_0,b_1)$ , 且

\lim_{\varepsilon\to0}|C_\infty^\varepsilon|=b_1-b_0,

使每个 $b\in C_\infty^\varepsilon$ 对应一个时间准周期斑块解。其区域写为

D_t^\varepsilon=\left\{\ell e^{i\theta}:0\leq\ell<R^\varepsilon(b,t,\theta)\right\},\qquad R^\varepsilon(b,t,\theta)=\sqrt{b^2+2\varepsilon r(t,\theta)},

其中

r(t,\theta)=\sum_{j\in S}a_j\cos\big(j\theta+\omega_j(b,\varepsilon)t\big) +p\big(\omega_p(b,\varepsilon)t,\theta\big).

频率向量 $\omega_p=(\omega_j)_{j\in S}$ 满足文中所需的 Diophantine 非共振条件, 且逐分量有

\omega_j(b,\varepsilon)\longrightarrow-\Omega_j(b),\qquad \Omega_j(b)=\frac12\big(j-1+b^{2j}\big),\qquad \varepsilon\to0.

修正项 $$p$$ 为偶函数, 对某个由 $$S$$ 决定的充分大正则性指标 $$s=s(S)$$ , 满足

\|p\|_{H^s(\mathbb T^{|S|+1})}=o(\varepsilon).

所得解线性稳定; 文中对“线性稳定”的含义是, 该准周期解线性化系统的 Floquet 指数全为纯虚数。该结论由构造中的可约化论证得到, 并非另行建立的独立稳定性估计。

当 $$|S|=1$$ 且环境边界实解析时, 定理 1.1 与定理 1.3 合用可知, 该周期解不能来自纯粹的刚性旋转, 因而斑块边界必有随时间变化的振荡。定理 1.3 的环境边界本身不必实解析; 但作者要求的对称性是构造假设的一部分。摘要与定理摘录没有提供 Diophantine 常数、 $\varepsilon_0$ 对参数的显式依赖, 也没有给出可直接计算的误差常数。

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

对定理 1.1, 核心不是直接求解 V 态轮廓方程, 而是把刚性旋转转化为不同旋转时刻的流函数之间的比较。令 $\varphi_\theta(x)=\Psi(\theta/\Omega,R_\theta x)$ , 并在 $D\cap D_\theta$ 中考察 $h_\theta=\varphi_0-\varphi_\theta$ , 其中 $D_\theta=R_{-\theta}D$ 。旋转协变性和斑块随 $R_{\Omega t}$ 移动给出

\Delta\varphi_\theta=\mathbf 1_{D_0},\qquad \Delta h_\theta=0.

另一方面, 轮廓动力学积分后, 流函数在随斑块旋转的边界上满足常数相对旋转势的关系, 从而 $h_\theta$ 在 $\partial D_0$ 上为常数。最大值原理先将该常值关系传入 $$D_0$$ ; 调和函数的解析性与唯一延拓再将其传到 $D\cap D_\theta$ 的相应连通分支。这一步把边界运动信息转化为一族调和函数恒等式。

反证的几何环节利用非圆形边界必存在横截方向: 可找到边界点, 使旋转方向不与边界相切。对该点附近的边界弧作小角度旋转, 可得到位于 $$D$$ 内的开弧 $I_\theta\subset\partial D_\theta\cap D$ 。作者进一步证明初始流函数可延拓到外边界邻域, 并利用边界实解析性及跨解析弧的广义 Schwarz 反射原理, 将 $I_\theta$ 上的常值延拓到整条 $\partial D_\theta$ 。原外边界上 $\Psi=0$ , 而两条围住相同面积且相交的 Jordan 曲线有公共点, 因而上述常值只能为零。随 $\theta$ 变化的边界弧扫出斑块外部的一个开集, 流函数在该开集为零; 再由它在 $D\setminus\overline{D_0}$ 中调和及唯一延拓, 得到斑块外部流函数恒为零。于是任取围住斑块的 Jordan 曲线 $\Gamma$ 都有 $\oint_\Gamma v\cdot dp=0$ , 但 Stokes 定理给出

\oint_\Gamma v\cdot dp=\int_{\operatorname{int}(\Gamma)}\omega\,dx=|D_0|>0,

矛盾。真正的障碍是如何从局部边界弧上的信息扩展到整个旋转边界; 实解析反射与调和唯一延拓正是完成这一扩展的关键, 也解释了为何当前论证要求解析边界。

对定理 1.2, 圆形点涡轨道是 Robin 函数的一个等值线。作者使用 $\Delta R_D=4e^{2R_D}$ 和最大值原理, 先推出 Robin 函数在轨道所围区域内具有径向性, 再利用 $$R_D$$ 的实解析性将径向性质延拓到整个 $$D$$ 。最后结合 Robin 函数靠近边界时趋于 $+\infty$ 的行为, 推出环境区域必须为圆盘。摘录未展示该证明的全部边界估计细节, 因此不能据此补充更精细的正则性论证。

定理 1.3 的困难在于, 对固定环境区域的轮廓动力学做准周期解构造时, 线性化算子同时含有切向输运、无限维法向频谱和小除数; 还必须控制非径向几何扰动对共振条件及线性可约化的影响。作者从单位圆盘构造出发, 利用

G_D(x,y)=G_{\mathbb D}(x,y)+\varepsilon^2G_p^\varepsilon(x,y),

其中远离边界时的扰动核 $G_p^\varepsilon$ 光滑, 将几何变化纳入小扰动。采用 $R^\varepsilon=\sqrt{b^2+2\varepsilon r}$ 的缩放, 使斑块形变为 $\varepsilon$ 阶而域形变为 $\varepsilon^2$ 阶; 因而主阶线性化保留圆盘 Rankine 斑块的频率

\Omega_j(b)=\frac12(j-1+b^{2j}).

固定有限切向模态集 $$S$$ , 将相空间拆成切向与法向部分, 对切向 Fourier 系数作辛极坐标变换, 引入作用量 $$I$$ , 角变量 $\vartheta$ 和法向变量 $$z$$ 。寻找承载准周期轨道的不变环面, 等价于求嵌入 $i=(\vartheta,I,z)$ 的非线性方程

F(i,\omega)=\omega\cdot\partial_\phi i-X_{H\circ A}(i)=0.

以圆盘问题的准周期解为零阶近似, 再实施 Nash-Moser 迭代。每步需要构造线性化算子的近似右逆: 切向作用角变量部分按 Berti-Bolle 框架处理, 法向部分通过 KAM 可约化对角化并求逆。域形变产生的项被组织为小范数的光滑积分算子或输运项中小的变系数局部扰动, 使未扰动问题的可逆估计能够在扰动下延续。这里的技术要点不是单独使用 Nash-Moser 或 KAM, 而是将作用-角-法向分解、法向可约化和可逆性估计拼接成适用于几何微扰轮廓方程的近似右逆。

环境映射的实对称条件使 Hamiltonian 轮廓动力学具有可逆性, 即 $r(t,\theta)$ 为解时 $r(-t,-\theta)$ 也为解。作者将这种对称结构用于 KAM 可约化及 Nash-Moser 方案; 文中说明它对当前证明是重要条件。在线性化系统经可约化后, 得到常系数系统且谱为纯虚数, 因而线性稳定性随构造一并得到。文中还指出可在多重对称函数空间内实施方案以保留 $$m$$ 重对称性, 但为避免增加复杂性并未在本文完成该版本。

4. 与之前工作的比较

本文的刚性结论与圆盘、全平面中的 V 态构造形成互补。全平面中的 Rankine 圆盘是平凡旋转态, Kirchhoff 椭圆给出经典二重对称旋转族; Burbea 通过局部分岔构造了由圆盘产生的多重对称 V 态, 后续工作研究其全局延拓。圆盘中的对应理论包括 de la Hoz、Hassainia、Hmidi 与 Mateu 关于内外半径参数化 Rankine 斑块分岔的结果, 以及点涡轨道的斑块正则化工作。本文并不改进这些圆对称环境中的构造, 而是反过来证明: 在有界单连通实解析环境区域中, 若一个严格内含的非平凡斑块以非零角速度刚性旋转, 则环境只能是圆盘。因此, 圆盘中已有 V 态不是可直接移植到一般边界的结构; 一般非圆域的周期运动若存在, 必须具有更复杂的时间演化。此结论不覆盖静态斑块, 也未证明一般非圆域中不存在任何周期解。

与非圆域中的既有稳态斑块研究相比, 文中明确提到 [14, 82, 94] 给出了平面区域中的非平凡稳态斑块。因此定理 1.1 的限制针对 $\Omega\ne0$ 的刚性旋转, 并非对所有时间独立或时间依赖解的普遍刚性结论。作者也强调环境边界实解析性在当前证明中不可缺少, 尽管相信较低正则性下结论可能仍成立; 这仍是未解决方向。对斑块本身单连通的假设, 作者认为可尝试推广到有孔斑块, 但需要处理不同边界上的相对流函数常数; 对多连通环境区域, 作者则指出同一结论未必成立。相比之下, 点涡定理只要求 $$C^1$$ 边界, 但它处理的是 Robin 函数产生的单点轨道, 并非有限面积斑块的轮廓方程。

准周期构造方面, 早期工作 [75]、[59] 及 Euler- $\alpha$ 相关研究 [86] 借助模型中的外部参数避开共振; 后续 SQG 工作 [49] 处理了相应参数范围中的剩余情形。二维 Euler 方程没有可供调节的外部参数, 因而 [9] 使用 Kirchhoff 椭圆的形状参数, Hassainia 与 Roulley 使用圆盘内 Rankine 斑块半径, [61] 则使用环形斑块内半径。本文延续以平衡态几何参数代替外部参数的思路, 但改变了环境区域: 从全平面或圆盘推进到圆盘的小非径向扰动, 构成作者所称的首个非径向环境中的准周期涡斑块构造。其基态频率仍来自圆盘 Rankine 斑块的线性化, 半径 $$b$$ 在满测度极限的 Cantor 参数集中承担非共振参数的作用。

该推进有明确代价。结论只在环境形变为 $\varepsilon^2$ 阶、斑块振幅为 $\varepsilon$ 阶的相容缩放下建立, 并限于固定有限模态集 $$S$$ 、半径区间 $[b_0,b_1]\Subset(0,1)$ 及充分小的扰动。作者称尚不清楚若不采用这一相容尺度, 能否运行证明。构造还要求共形映射的实对称条件, 以保证可逆性; 因此它不是对任意小的非径向形变都适用的无条件结果。虽然参数集合在 $\varepsilon\to0$ 时渐近占满 $$(b_0,b_1)$$ , 对每个正 $\varepsilon$ 仍是 Cantor 型集合, 而不是所有半径均可用。文中未给出可供比较的显式共振排除常数或统一误差常数。

本文的准周期论证以 [64] 的圆盘构造为直接基础, 采取扰动延续而非从头建立一套不同的准周期理论。所增加的工作是把环境几何引起的 Green 函数修正纳入作用-角-法向 Nash-Moser/KAM 框架, 并利用其小而光滑的结构维持近似右逆与可约化。相较于 [64], 主要推进是环境区域从圆盘扩展到特定的小非径向扰动; 摘录并未表明它扩大了模态数、半径区间或振幅范围。关于周期情形, 定理 1.1 使本文构造进一步具有非刚性结论: 当 $$|S|=1$$ 且外边界实解析时, 所得周期斑块不能只是整体旋转。这与 [62] 等关于点涡周期轨道正则化后出现丰富周期动力学的研究相呼应, 但二者证明对象和构造机制不同。除上述明确内容外, 文中没有在提供的摘录中对所有相关先前工作的定量常数、最优正则性或最优尺度作系统比较, 因此这些方面无法确认。

2. 旋转对称环面区域中的 Seifert 型磁流体静力平衡态

Seifert type magnetohydrostatic equilibria in rotationally symmetric toroidal domains

作者: Daniel Peralta-Salas, Radu Slobodeanu

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

发布日期: 2026-10-06

arXiv ID: 2610.08519v1

摘要翻译: 我们构造了旋转对称环面区域中的磁流体静力平衡态，其磁力线定义了任意互素整数 $$m,n$$ 对应的 $$(m,n)$$ Seifert 纤维化。

Comment: 8 pages. Comments are welcome!

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

1. 研究对象和背景

本文研究三维欧氏空间中轴对称环形区域内的磁流体静力平衡, 重点不是单纯构造满足平衡方程的磁场, 而是要求磁力线具有指定的闭轨道拓扑: 环形区域中的磁力线构成具有给定绕行数的塞弗特纤维化, 且磁轴是唯一可能的特殊纤维。

磁流体静力平衡的未知量是磁场 $$B$$ 和等离子体压力 $$P$$ . 在光滑有界区域 $\Omega\subset\mathbb{R}^3$ 中, 它们满足

B\times\operatorname{curl}B+\nabla P=0,\qquad \operatorname{div}B=0\quad\text{in }\Omega,

以及边界条件

B\cdot N=0,\qquad P|_{\partial\Omega}=\text{常数}.

其中 $$N$$ 是外法向量. 第一式表示磁力与压力梯度平衡, 第二式表示磁场无散, 边界条件则要求磁场沿区域边界切向并使边界成为等压面。

在轴对称情形, 取柱坐标 $(r,\varphi,z)$ , 磁场可由流函数 $\psi(r,z)$ 和函数 $F(\psi)$ 表示为

B=-\frac{\partial_z\psi}{r}e_r+\frac{F(\psi)}{r}e_\varphi+\frac{\partial_r\psi}{r}e_z.

相应的平衡条件归结为 Grad-Shafranov 方程

\partial_{rr}\psi+\partial_{zz}\psi-\frac1r\partial_r\psi =r^2H'(\psi)-\frac12(F^2)'(\psi),

其中 $$H$$ 是流函数的函数, 按本文符号约定 $$H=-P$$ . 这一降维使平衡场构造转化为二维半线性椭圆方程问题; 同时, $\psi$ 的等值线绕旋转轴生成磁通面, 磁场轨道在这些面上运动。因此, 方程解的几何结构直接控制磁力线的闭合性、绕行类别和相互链接方式。

作者考虑以圆 $\{r=1,z=0\}$ 为中心线的细环形区域。对给定的互素非零整数 $$m,n$$ , 目标是使每条非磁轴磁力线都闭合, 并在经向和环向分别绕行 $$m$$ 次和 $$n$$ 次, 从而具有同调类 $m\mu+n\lambda$ . 这里 $\mu$ 是经线, $\lambda$ 是沿旋转方向的纬线。相应的轨道叶状结构是塞弗特纤维化; 中心磁轴为唯一的核心纤维, 其链接数满足 $\operatorname{lk}(\Gamma,C)=-m$ 。

既有轴对称构造已经覆盖若干特殊轨道类型, 但任意指定经向、环向绕行数并在欧氏环形区域中同时满足平衡方程, 并非只靠给磁场指定一个旋转变换数即可完成: 还需要先找到具有合适等值线几何的 Grad-Shafranov 解, 并严格控制其光滑性和边界。文章要补足的正是这一构造性空白。其策略是先在单位圆盘上构造一个近似径向、且可由保面积微分同胚与标准径向函数共轭的椭圆方程解, 再将其嵌入细环形区域, 并以常数环向分量调节轨道斜率。该方法给出的不仅是某个孤立绕行类型的例子, 而是任意互素非零 $$(m,n)$$ 的闭轨道构造; 同一解析构造还允许非有理斜率, 产生每个磁通面上的稠密轨道情形。

问题的重要性在于, 理想磁流体中的磁力线组织方式是磁约束与平衡几何的核心信息, 而闭合磁力线的拓扑类型会影响对磁轴、磁通面及链接结构的理解。本文的结论是数学上的存在性结果, 并不提供稳定性、动力学吸引性或实验可实现性的判据。文中也没有证明一般环形区域、任意给定压力剖面或任意给定磁通面几何下都存在此类平衡; 它建立的是一族细环形区域中的明确构造。

2. 主要定理或主要结果

主要结论是定理 1.1. 对任意整数 $k\ge2$ 和 $0<\alpha<1$ , 存在 $\varepsilon_0=\varepsilon_0(k,\alpha)>0$ , 使得对所有 $0<\varepsilon<\varepsilon_0$ 及所有互素非零整数 $$m,n$$ , 存在光滑旋转对称环形区域 $\Omega_\varepsilon$ 和轴对称平衡场 $B\in C^{k,\alpha}(\Omega_\varepsilon)$ , 满足磁流体静力平衡方程及边界条件。区域取为

\Omega_\varepsilon=\{(r\cos\varphi,r\sin\varphi,z):(\log r)^2+z^2<\varepsilon^2,\ \varphi\in\mathbb{R}/(2\pi\mathbb{Z})\},

它是圆 $\{r=1,z=0\}$ 的厚度为 $O(\varepsilon)$ 的管状邻域。定理所用的 $\varepsilon_0$ 、区域和压力构造不依赖于 $$m,n$$ 。

所得磁场处处非零, 压力 $P\in C^{k,\alpha}(\Omega_\varepsilon)$ 的等值面构成包围一条临界核心曲线的嵌套环面叶状结构。核心曲线 $$C$$ 本身是磁力线, 且可写为

C=\{(r_\varepsilon\cos\varphi,r_\varepsilon\sin\varphi,z_\varepsilon):\varphi\in\mathbb{R}/(2\pi\mathbb{Z})\},

其中

r_\varepsilon=1+\varepsilon^2+O(\varepsilon^3),\qquad z_\varepsilon=O(\varepsilon^3).

每条非核心磁力线 $\Gamma$ 都闭合, 位于某个压力等值环面上, 并具有绕行类

[\Gamma]=m\mu+n\lambda.

按文中经线、纬线的取向约定, 其与核心的链接数为

\operatorname{lk}(\Gamma,C)=-m.

因此轨道叶状结构是塞弗特纤维化, 核心 $$C$$ 是唯一特殊纤维; 正则纤维沿经向和环向分别绕行 $$m$$ 次和 $$n$$ 次。由于 $$m,n$$ 可为负数, 绕行方向由其符号决定; 核心纤维的轨道周期在文中归一化的圆作用下为 $2\pi/|n|$ , 其特殊纤维重数对应 $$|n|$$ 。特别地, $$|n|=1$$ 时核心不具有非平凡的特殊重数。

构造中令 $\kappa=n/m$ , 并取常数 $F(\psi)=\kappa$ . 设 $\Psi_\varepsilon$ 是单位圆盘上的函数, 令

r=e^{\varepsilon x},\qquad z=\varepsilon y,\qquad \psi(r,z)=\varepsilon^2\Psi_\varepsilon(x,y).

作者构造 $\Psi_\varepsilon$ 使其满足

L_\varepsilon\Psi_\varepsilon=f_\varepsilon(\Psi_\varepsilon),\qquad L_\varepsilon=e^{-4\varepsilon x}(\partial_{xx}-2\varepsilon\partial_x)+e^{-2\varepsilon x}\partial_{yy}.

定义

H(s)=\varepsilon^2\int_0^{s/\varepsilon^2}f_\varepsilon(c)\,dc,

则 $H'(\psi)=f_\varepsilon(\Psi_\varepsilon)$ , 从而得到 $\Delta^*\psi=r^2H'(\psi)$ , 即 $$F$$ 为常数时的 Grad-Shafranov 方程。相应磁场为

B=-\frac{\partial_z\psi}{r}e_r+\frac{\kappa}{r}e_\varphi+\frac{\partial_r\psi}{r}e_z.

由于 $\kappa\ne0$ , 环向分量处处非零, 因而 $$B$$ 不消失。边界上 $\psi=\varepsilon^2/2$ , 故磁场与边界相切; 压力按文中约定为 $P=-H(\psi)$ 。

更具体地, 命题 2.2 给出小 $\varepsilon$ 下保面积微分同胚 $\gamma_\varepsilon:D\to D$ 和径向函数 $f_\varepsilon$ , 使得

\det D\gamma_\varepsilon=1,\qquad \Psi_\varepsilon=\Psi_0\circ\gamma_\varepsilon^{-1},\qquad \Psi_0(x,y)=\frac{x^2+y^2}{2},\qquad L_\varepsilon\Psi_\varepsilon=f_\varepsilon(\Psi_\varepsilon).

$\Psi_\varepsilon$ 有唯一非退化极小点, 且 Hamilton 向量场 $X_{\Psi_\varepsilon}=J\nabla\Psi_\varepsilon$ 的非定常轨道最小周期均为 $2\pi$ . 这使轨道闭合性不只是近似结论, 而是由精确的保面积共轭得到。

在坐标 $(x,y,\varphi)$ 中, 磁场轨道满足

r^2B=X_{\Psi_\varepsilon}+\kappa\partial_\varphi.

再经 $\gamma_\varepsilon^{-1}\times\operatorname{id}_{S^1}$ 共轭, 归一化向量场变成

m(-Y\partial_X+X\partial_Y)+n\partial_\varphi.

其流是圆作用

(X+iY,e^{i\varphi})\longmapsto(e^{imt}(X+iY),e^{int}e^{i\varphi}).

互素条件保证非核心轨道在归一化作用下的最小周期为 $2\pi$ , 并给出所需的两种绕行数; 圆盘中心对应核心轨道。定理中的核心位置估计来自隐函数定理解对 $\varepsilon$ 的一阶展开, 文中得到 $\gamma_\varepsilon(0)=(\varepsilon,0)+O(\varepsilon^2)$ 。

文章另指出同一构造可取任意实数 $\kappa$ . 当 $\kappa\notin\mathbb{Q}$ 时, 每个非核心不变环面上的轨道为稠密轨道; 当 $\kappa=0$ 时, 可得所有轨道闭合且为经向的情形。这两项是定理构造的直接延伸, 但不属于定理 1.1 中“ $$m,n$$ 均非零”的表述。文中还提醒, 归一化向量场的周期不等于物理磁场 $$B$$ 的周期, 后者可以随不变环面变化。

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

证明的主要困难在于, 需要同时满足两个性质不同的要求: 一方面, 流函数必须解一个系数依赖空间变量的半线性椭圆方程; 另一方面, 为了精确控制磁力线, 它的等值线 Hamilton 流必须保留标准圆周运动的轨道结构。单独对椭圆方程作小扰动估计, 通常只能控制解的接近程度, 并不能保证所有轨道闭合或绕行数精确。本文的核心做法是把“构造椭圆方程解”和“保持轨道周期结构”合并为一个保面积坐标变换问题, 再以 Banach 空间隐函数定理解决。

首先将环形区域经 $r=e^{\varepsilon x}, z=\varepsilon y$ 拉回单位圆盘. 在这些变量下, 归一化的 Grad-Shafranov 算子成为

L_\varepsilon=e^{-4\varepsilon x}(\partial_{xx}-2\varepsilon\partial_x)+e^{-2\varepsilon x}\partial_{yy},

它是拉普拉斯算子的系数扰动。直接寻找 $\Psi_\varepsilon$ 可能无法同时控制其等值线的周期。于是作者从标准函数 $\Psi_0(q)=|q|^2/2$ 出发, 寻找面积保持映射 $\gamma(q)=q+v(q)$ , 并令 $\Psi=\Psi_0\circ\gamma^{-1}$ . 面积保持条件 $\det D\gamma=1$ 使 $\gamma$ 成为标准辛形式下的辛映射, 因而 Hamilton 向量场满足

\gamma_*X_{\Psi_0}=X_\Psi.

所以, 一旦找到这样的 $\gamma$ , $X_\Psi$ 的非定常轨道就与标准旋转场共轭, 其最小周期精确为 $2\pi$ . 这一步是得到全部磁力线闭合而非仅近似闭合的关键。

具体地, 定义 $A=D\gamma$ , $W=-JR\gamma$ , $M=(DW)A^{-1}$ , 并构造误差函数

E_\varepsilon(v)=e^{-4\varepsilon\gamma_1}(M_{11}-2\varepsilon W_1)+e^{-2\varepsilon\gamma_1}M_{22}.

在 $\det D\gamma=1$ 时, 有 $W=(\nabla\Psi)\circ\gamma$ 、 $M=(D^2\Psi)\circ\gamma$ , 因而

E_\varepsilon(v)=(L_\varepsilon\Psi)\circ\gamma.

这将非线性椭圆方程转化为对 $$v$$ 的条件: $E_\varepsilon(v)$ 必须在圆盘上为径向函数。作者把面积保持、径向性及边界落在单位圆上编码为映射

F(v,\varepsilon)=\left(\det(I+Dv)-1,\ QE_\varepsilon(v),\ \Pi_{\partial}|q+v|^2-1\right),

其中 $$Q=I-P$$ , $$P$$ 是对角向旋转取平均的投影, $\Pi_{\partial}$ 去掉边界平均值。三个分量分别施加面积保持约束、消除误差函数的非径向部分、以及固定像域的边界形状。边界条件只去掉均值而非强行逐点固定, 是因为映射的边界像先由该条件成为以原点为中心的圆, 再由面积保持条件推出圆半径为一。

技术上的关键不是笼统应用隐函数定理, 而是证明该非线性映射在基点的线性化是 Banach 空间同构。作者使用带图范数的函数空间 $$X_s,A_s,G_s,H_s$$ , 其中控制旋转导数的正则性, 并对旋转平均为零的函数定义角向积分算子 $$S$$ , 满足 $$RSg=g$$ . 在 $(v,\varepsilon)=(0,0)$ 处, 线性化为

T v=(\operatorname{div}v,\ R\operatorname{curl}_2v,\ \Pi_{\partial}(v\cdot\nu)).

给定目标 $$(a,g,h)$$ , 先令 $$w=Sg$$ , $\lambda=(2\pi)^{-1}\int_Da$ , 再分别解 Dirichlet 问题 $\Delta b=w$ , $b|_{\partial D}=0$ 和具有相容 Neumann 数据的 Neumann 问题 $\Delta e=a$ , $\partial_\nu e=h+\lambda$ , $\int_D e=0$ . 取

v=\nabla e+J\nabla b,

即可构造 $$Tv=(a,g,h)$$ . Schauder 估计以及对旋转导数所满足的相应椭圆问题给出 $v\in X_s$ 和有界逆估计。单射性则由 $$Rv$$ 的约束、散度为零、旋度为零以及边界法向分量的均值条件推出 $$v$$ 是具有零 Neumann 数据的调和势梯度, 因而只能为零。得到 $$T$$ 可逆后, 隐函数定理产生小 $\varepsilon$ 的解 $v_\varepsilon$ 。

该论证还处理了正则性与几何闭合之间的衔接。取 $$s=2l$$ 后, $E_\varepsilon(v_\varepsilon)$ 是径向函数; 将其限制到直径得到偶函数 $h_\varepsilon(t)$ , 再用 Whitney 关于偶函数的结论写成 $h_\varepsilon(t)=f_\varepsilon(t^2/2)$ , 从而得到所需的单变量非线性项 $f_\varepsilon(\Psi_\varepsilon)$ . 这确保椭圆方程确实具有半线性形式, 而不是仅有某个随空间变化的右端项。

最后, 作者把圆盘解缩放回欧氏环形区域, 将 $$F$$ 取为常数 $\kappa=n/m$ . 轨道分析利用 $r^2B=X_{\Psi_\varepsilon}+\kappa\partial_\varphi$ ; 保面积共轭将其化成圆盘旋转与环向平移的直积流。整数权重 $$m,n$$ 的互素性恰好确定非核心轨道的共同最小周期及其绕行数, 而圆盘中心的稳定子给出核心特殊纤维。由此, PDE 解的构造、环面纤维化和链接数计算形成一条明确的证明链条。核心曲线的 $O(\varepsilon^2)$ 径向偏移则另由隐函数解的一阶展开计算, 而非仅凭对称性断言。

文中并未建立该平衡的稳定性、唯一性或对一般几何扰动的鲁棒性。其分析框架依赖细管缩放后算子接近圆盘上的拉普拉斯算子, 因而核心创新是适用于小 $\varepsilon$ 的构造与保面积共轭, 不是对任意环形区域上 Grad-Shafranov 方程的一般求解理论。

4. 与之前工作的比较

本文在轨道拓扑范围上的主要推进, 是把轴对称欧氏环形区域中已知的若干特殊闭轨道类型扩展到任意互素非零整数 $$(m,n)$$ . 引言明确提到, Lortz [4] 在非轴对称环形区域构造了所有轨道闭合、类型为 $$(0,1)$$ 的磁流体静力平衡; 轴对称情形下, 已有许多所有轨道均为经向闭合、类型为 $$(1,0)$$ 的例子, 相关工作见 [7]. 本文的构造同时允许非零经向和环向绕行数, 并规定正则轨道的同调类为 $m\mu+n\lambda$ , 因而不是仅对已有纯经向或单一特殊类型的重述。

对于具体例子, 引言指出 Solov’ev 解 [10] 可显式给出 $$(m,n)=(2,1)$$ 的情形, 近期 Landreman [3] 又利用该例构造了 Grad 猜想的反例。相较之下, 本文的结论不是将一个显式解推广为另一个孤立的显式参数值, 而是对全部互素非零 $$m,n$$ 给出统一存在性构造。代价是结果属于小尺度、隐函数定理型构造: 文中没有给出任意 $$m,n$$ 下的闭式流函数, 也没有提供可直接计算的完整压力剖面公式; $f_\varepsilon$ 和相应解是通过 Banach 空间隐函数定理得到的。

已有工作 [5,9] 在赋予标准度量的 $$S^3$$ 上构造了具有任意塞弗特纤维化的解析磁流体静力平衡。作者强调, 欧氏空间中的证明要复杂得多。本文针对的是旋转对称的欧氏细环形区域, 通过 Grad-Shafranov 方程和保面积微分同胚完成构造, 因此把拓扑类型的广泛可实现性从 $$S^3$$ 上的设置带到特定的欧氏环形几何中。两种背景的度量、方程表达及区域几何不同, 本文并未声称将 $$S^3$$ 上的构造逐项搬运, 也未证明任意三维流形或任意欧氏环域中的相应结论。

适用范围方面, 定理要求 $\varepsilon$ 足够小, 区域是 $\{(\log r)^2+z^2<\varepsilon^2\}$ 这种围绕单位圆的薄管。小性阈值依赖正则性指标 $$k$$ 和 Hölder 指数 $\alpha$ , 但作者说明它不依赖 $$m,n$$ . 这种与绕行整数无关的阈值来自解析构造先生成一个与 $$m,n$$ 无关的流函数, 再独立选择常数斜率 $\kappa=n/m$ . 但这并不意味着所有场量的范数对 $$m,n$$ 一致有界: $\kappa$ 可随整数比值变大, 且环向分量为 $\kappa/r$ . 文中也没有给出对 $$m,n$$ 的统一场强估计。

结论还扩展到有理性以外的轨道动力学。作者指出同一构造可令 $\kappa\in\mathbb{R}$ ; 取无理 $\kappa$ 时, 各非核心不变环面上的轨道稠密, 取 $\kappa=0$ 时得到纯经向闭轨道。这说明所需椭圆解并非只适配某一组整数绕行数, 而可与连续变化的环向漂移组合。不过, 定理 1.1 本身只陈述 $$m,n$$ 为互素非零整数的闭轨道情形; 对无理斜率的表述见正文说明, 不应误写成定理 1.1 的整数型结论。

与 Constantin、Drivas 和 Ginsberg [1] 的关系主要是技术启发。作者明确说明, 保面积变形径向函数以适配半线性椭圆方程的思路受其定理 3.1 启发; 本文进一步将该策略用于具体的 $L_\varepsilon$ 算子, 并通过线性算子 $$T$$ 的同构性质实现隐函数构造。相关闭轨道判据方面, 文中引用 Perrella [6], 但并未在所提供内容中逐条比较其假设或给出本文相对于该判据的定量改进, 因而无法确认更细的优劣差别。

需要注意文中一个关于量级的表述与显式公式之间存在表面不一致。第 2 页 Remark 1.1 称平衡场及其导数为 $O(\varepsilon^2)$ , 但第 6 页的构造公式含有环向分量 $B_\varphi=\kappa/r$ , 对固定非零 $\kappa=n/m$ 通常是 $$O(1)$$ , 而经向分量按缩放为 $O(\varepsilon)$ . 因此, 从提供的公式无法确认“整个磁场为 $O(\varepsilon^2)$ ”如何成立; 严谨比较时不宜把该量级断言当作与场公式一致的统一估计。可以确定的是区域厚度为 $O(\varepsilon)$ , 核心位置满足文中给出的展开, 以及构造阈值不依赖 $$m,n$$ 。

最后, 本文提供的是存在性和精确轨道拓扑结论, 其代价与局限也应明确: 只处理轴对称的特定薄环域; 压力剖面由构造中的 $f_\varepsilon$ 产生, 并非任意预先指定; 没有证明平衡的稳定性、唯一性、对区域形状扰动的持续性, 也没有讨论非轴对称扰动下磁力线是否仍闭合。文中关于 Grad 猜想反例的联系是 Solov’ev 的 $$(2,1)$$ 例子已被 [3] 使用, 并非本文声称自行构造新的反例。

3. 临界横向低频控制下 $$$$

Nonlinear stability of 3D Couette flow on $\mathbb T\times\mathbb R^2$ under critical transverse low-frequency control

作者: Yubo Chen, Guoxu Yang, Jianbo Yu

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08288v1

摘要翻译: 我们证明了 $$$$ 上不可压 Navier–Stokes 方程三维库埃特流的非线性稳定性。满足 $$$$ 的无散扰动生成唯一的全局 $$H^2$$ 解，且这些解始终保持在库埃特流附近；其中， $$$$ 表示横向速度的流向平均。尽管不存在跨向谱隙，该结果仍保留了 $$O(ν)$$ 的稳定性阈值。非零流向模态表现出增强耗散，其壁法向分量满足一个与黏度无关的时间积分型无粘阻尼估计，而横向零模态在 $$L^2$$ 中以 $(1+νt)^{-1}$ 的速率衰减。端点空间 $$$$ 控制横向低频，容纳散度形式的非线性强迫，并且对于壁法向零模态所产生的线性 lift-up 效应是临界的。证明结合了该端点范数的传播、Fourier 分裂，以及一个双层 roll-to-streak 表示，并对 Wei–Zhang (CPAM, 2021) 的非零模态框架作了调整。

Comment: 无

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

1. 研究对象和背景

论文研究无界横向区域 $\Omega=\mathbb T_x\times\mathbb R_y\times\mathbb R_z$ 上三维不可压 Navier–Stokes 方程在平面 Couette 流 $$U=(y,0,0)$$ 附近的长期稳定性。令 $$v=U+u$$ ，其中 $$u=(u^1,u^2,u^3)$$ 是扰动， $$p$$ 是相应压力，则扰动满足

\partial_t u-\nu\Delta u+y\partial_xu+u^2e_1+u\cdot\nabla u+\nabla p=0,\qquad \nabla\cdot u=0,

其中 $0<\nu\leq1$ 是黏性系数， $$e_1=(1,0,0)$$ 。问题是：当 $\nu\to0$ 时，初始扰动需要多小，才能使解对所有时间保持在 Couette 流附近？这里的空间结构尤其重要：流向 $$x$$ 周期，而壁法向 $$y$$ 和展向 $$z$$ 均无界。对任意场 $$f$$ ，记流向平均和非零流向模态为

f_0(y,z)=\int_{\mathbb T}f(x,y,z)\,dx,\qquad f_{\ne}=f-f_0,

并写 $u^\perp=(u^2,u^3)$ 、 $\nabla_\perp=(\partial_y,\partial_z)$ 。流向平均的横向速度 $u_0^\perp$ 称为 roll，流向速度 $$u_0^1$$ 称为 streak。本文的关键并非只研究线性谱稳定性，而是要同时处理非正规线性增长、剪切混合和黏性耗散对不同模态产生的不同效应。

对 $k\ne0$ 的流向模态，剪切输运 $y\partial_x$ 使频率随时间倾斜，黏性因而在比普通扩散更短的时间尺度上起效。被动标量的线性耗散乘子满足

\exp\!\left[-\nu\int_0^t\big(k^2+(\eta-ks)^2+\ell^2\big)\,ds\right] \leq \exp\!\left(-\frac{\nu k^2t^3}{12}\right),

对应增强耗散尺度 $t\sim\nu^{-1/3}$ 。同时，Euler 线性化中壁法向非零模态满足 $(\partial_t+y\partial_x)\Delta u^2_{\ne}=0$ ；从涡量恢复速度时出现随时间变化的椭圆乘子，产生 inviscid damping。非零模态的这些机制已有相当成熟的分析框架，但流向零模态不受剪切混合，仍受 lift-up 作用支配。

具体地，零模态线性系统为

\partial_tu_0^1-\nu\Delta_\perp u_0^1=-u_0^2,\qquad \partial_tu_0^\perp-\nu\Delta_\perp u_0^\perp=0,

从而

$$ u u $$

更正地，streak 的线性公式为

$$ u $$

$$ u $$

即

$$ u $$

并可直接写为

$$ u $$

其中实际公式是

u_0^1(t)=e^{\nu t\Delta_\perp}u^1_{0,\mathrm{in}}-t e^{\nu t\Delta_\perp}u^2_{0,\mathrm{in}}.

对横向频率 $\xi_\perp=(\eta,\ell)\ne0$ ，有 $\sup_{t>0}t e^{-\nu t|\xi_\perp|^2}=(e\nu|\xi_\perp|^2)^{-1}$ 。在展向周期的几何中，非零展向频率有谱隙；而在 $\mathbb R_z$ 上，横向频率可任意接近零，低频扩散任意缓慢，普通正阶 Sobolev 控制无法统一限制线性 lift-up。因此，已有周期展向问题中依靠谱隙得到的零模态估计不能直接移植到本区域。

本文针对的空白是：在没有展向谱隙、也没有额外色散机制的普通 Navier–Stokes 方程中，能否仍以自然的 $O(\nu)$ 初值尺度获得全局非线性稳定性。作者以临界齐次 Besov 条件控制初始 roll 的横向低频，并证明该条件可在非线性演化中传播。该范数为

\|f\|_{\dot B^{-2}_{2,\infty}(\mathbb R^2)}= \sup_{j\in\mathbb Z}2^{-2j}\|\dot\Delta_j f\|_{L^2(\mathbb R^2)},

且热半群刻画给出 $\sup_{t>0}\|t e^{\nu t\Delta_\perp}f\|_2\simeq\nu^{-1}\|f\|_{\dot B^{-2}_{2,\infty}}$ 。它既是线性 lift-up 的端点低频控制，也能容纳非线性方程中的散度型应力项，这两方面共同构成本文建立稳定性理论的出发点。

2. 主要定理或主要结果

定理 1.1 给出全局稳定性。若 $u_{\mathrm{in}}\in H^2(\Omega)$ 无散，且

\varepsilon:=\|u_{\mathrm{in}}\|_{H^2(\Omega)}+ \|u^\perp_{0,\mathrm{in}}\|_{\dot B^{-2}_{2,\infty}(\mathbb R^2)} \leq c_*\nu,

则扰动方程有唯一全局解

$$ u $$

也就是

$$ u $$

其准确正则性结论为

$$ u $$

$$ u $$

即

$$ u $$

具体地，

$$ u $$

可简洁表述为

u\in C([0,\infty);H^2(\Omega))\cap L^2_{\mathrm{loc}}([0,\infty);H^3(\Omega)),\qquad \sup_{t\geq0}\|u(t)\|_{H^2(\Omega)}\leq Cc_*.

并且对任意 $\tau_0>0$ ，有 $u\in C([\tau_0,\infty);H^4(\Omega))\cap L^2_{\mathrm{loc}}([\tau_0,\infty);H^5(\Omega))$ 。这里定理原文将一致 $$H^2$$ 界写为 $$Cc_*$$ ；其他定量模态估计以下述定理 1.5 为准。

定理 1.5 在相同假设下分别描述不同模态。存在与 $\nu$ 无关的 $\lambda_*>0,C>0$ ，使非零流向模态满足增强耗散

\|u_{\ne}(t)\|_{L^2(\Omega)}\leq C\varepsilon e^{-\lambda_*\nu^{1/3}t},

以及

\nu^{1/6}\|e^{\lambda_*\nu^{1/3}t}u_{\ne}\|_{L^2_tL^2} +\nu^{1/2}\|e^{\lambda_*\nu^{1/3}t}\nabla u_{\ne}\|_{L^2_tL^2} \leq C\varepsilon.

壁法向非零模态还满足黏性一致的时间积分 inviscid-damping 估计

\|e^{\lambda_*\nu^{1/3}t}u^2_{\ne}\|_{L^2_tH^1} +\|e^{\lambda_*\nu^{1/3}t}\partial_x\nabla u^2_{\ne}\|_{L^2_tL^2} \leq C\varepsilon.

其中没有在这些范数前损失 $\nu$ 的负幂；作者明确指出，这比单凭增强耗散积分得到的带 $\nu^{-1/6}$ 损失的 $$L^2_t$$ 界更强，但并未声称非线性问题有逐点 $t^{-2}$ 衰减率。

roll 的结论是

\|u^\perp_0(t)\|_{L^2(\mathbb R^2)}\leq\frac{C\varepsilon}{1+\nu t},\qquad \sup_{t\geq0}\|u^\perp_0(t)\|_{\dot B^{-2}_{2,\infty}} +\nu^{1/2}\|u^\perp_0\|_{L^2_tL^2}\leq C\varepsilon.

streak 在正时间平滑后满足

\sup_{t\geq1}\|u^1_0(t)\|_{H^4} +\nu^{1/2}\|\nabla_\perp\Delta_\perp^2u^1_0\|_{L^2(1,\infty;L^2)} +\nu^{-1}\sup_{t\geq1}\|\partial_tu^1_0(t)\|_{H^2} \leq C\frac{\varepsilon}{\nu}.

这是 uniform boundedness 控制，不是 streak 衰减结论；当 $\varepsilon\sim\nu$ 时，streak 可以达到常数量级，而 roll 和低阶非零模态仍为 $O(\varepsilon)$ 。

推论 1.7、1.9 给出便于验证的充分条件。可以将初始 roll 的 $\dot B^{-2}_{2,\infty}$ 控制替换为 $\dot H^{-2}(\mathbb R^2)$ 控制；也可以要求 $\||x_\perp|u^\perp_{0,\mathrm{in}}\|_{L^1}\lesssim\nu$ 。特别地，对任意固定 $\sigma>2$ ，条件 $\|\langle x_\perp\rangle^\sigma u^\perp_{0,\mathrm{in}}\|_2+\|u_{\mathrm{in}}\|_{H^2}\lesssim_\sigma\nu$ 也足够。加权条件只施加于初始 roll，不要求传播空间权重。

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

证明的首要障碍是零模态的低频累积，而不是单纯缺少耗散估计。roll 方程在 $\mathbb R^2$ 上没有 Poincaré 谱隙，故不能像横向周期情形那样由耗散直接得到统一的时间积分控制。另一方面，线性 lift-up 把 $$u^2_0$$ 经过时间因子 $$t$$ 积累到 streak；而由定理得到的 roll 衰减只有 $(1+\nu t)^{-1}$ ，时间上并不可积。因此，直接估计 $\int_0^t\|u^2_0(s)\|_2\,ds$ 会产生不允许的增长。作者的核心处理是将临界低频传播、roll 的衰减、非线性 lift-up 的表示和非零模态能量估计组成闭合框架，而不是把它们分开作粗略估计。

第一项关键设计是采用端点空间 $\dot B^{-2}_{2,\infty}$ 。线性热半群公式表明它恰好控制 $t e^{\nu t\Delta_\perp}u^2_{0,\mathrm{in}}$ 的 $$L^2$$ 范数，因而与 lift-up 的临界尺度匹配。更重要的是，二维 Leray 投影后的散度型应力满足

\|P_\perp\nabla_\perp\cdot F\|_{\dot B^{-2}_{2,\infty}} \lesssim \|F\|_{L^1(\mathbb R^2)}.

对 $F=u^\perp_0\otimes u^\perp_0+(u^\perp_{\ne}\otimes u^\perp_{\ne})_0$ ，右侧可由速度的 $$L^2$$ 控制。这使端点范数能够容纳 roll 自作用及非零模态平均应力。相比之下，较强的 $\dot H^{-2}=\dot B^{-2}_{2,2}$ 虽然控制线性 lift-up，却一般不包含 $$L^1$$ 张量的投影散度；因此作者不是在演化中维持 $\dot H^{-2}$ ，而是在初值层面允许该条件，并在证明中传播更适配非线性的 $\dot B^{-2}_{2,\infty}$ 。

第二项关键步骤是用 Fourier splitting 同时闭合端点范数传播与 roll 的时间积分。低频部分由 Besov 范数控制，高频部分由黏性耗散控制；roll 自相互作用在其 $$L^2$$ 能量恒等式中抵消，剩余外力主要是非零模态应力。文中得到的结构性估计为

\sup_{0\leq t\leq T}\|u^\perp_0(t)\|_{\dot B^{-2}_{2,\infty}} \lesssim \varepsilon+\int_0^T\|u^\perp_0(t)\|_2^2\,dt+\nu^{-1/3}\varepsilon^2,

以及

\int_0^T\|u^\perp_0(t)\|_2^2\,dt \lesssim \nu^{-1}\left(\varepsilon^2+ \sup_{0\leq t\leq T}\|u^\perp_0(t)\|_{\dot B^{-2}_{2,\infty}}^2\right) +\nu^{-3}\varepsilon^4.

在 $\varepsilon\leq c\nu$ 下联立吸收，得到端点范数一致有界及 $\nu^{1/2}L^2_tL^2$ roll 控制。再结合对波应力的加权估计和 Fourier splitting，推出 $(1+\nu t)^{-1}$ 的 roll 衰减。证明中还量化了 roll 自作用造成的端点范数变化至多为 $O(\varepsilon^2/\nu)$ ，在所取阈值下不会破坏初始控制。

第三项创新是先代入 roll 方程，再估计 streak 方程，避免直接积累不可积的 roll 衰减。由两个零模态方程组合得到 streak 的 Duhamel 表示，其中 roll 应力项带有时间权重 $$t-s$$ ：

\int_0^t(t-s)e^{\nu(t-s)\Delta_\perp} \left[P_\perp\nabla_\perp\cdot \left(u^\perp_0\otimes u^\perp_0+(u^\perp_{\ne}\otimes u^\perp_{\ne})_0\right)\right]_y(s)\,ds.

这个权重并非额外损失，因为它与热半群对散度源项的衰减相配合：

\sup_{\tau>0}\|\tau e^{\nu\tau\Delta_\perp}P_\perp\nabla_\perp\cdot F\|_2 \lesssim \nu^{-1}\|F\|_1.

于是非线性 lift-up 可由二次应力的时间积分控制；roll 自作用对 streak 的贡献为 $O(\varepsilon^2/\nu^2)$ 。作者还构造相应的高阶版本，以控制正时间平滑后所需的导数，而不破坏这套时间权重结构。

第四项是将上述零模态估计接入 Wei–Zhang 非零模态框架。作者把 streak 吸收到有效剪切速度 $$V=y+u^1_0$$ 中，在 $\|\nabla_\perp u^1_0\|_\infty$ 足够小时保证 $V_y\geq1/2$ ，再使用随 streak 调整的坐标、切向导数和 good unknowns，以保留剪切混合并处理速度与涡量间的非局部关系。相应能量同时记录 $$L^2$$ 控制、类似 $\nabla\Delta^{-1}\partial_x$ 的 inviscid-damping 控制、增强耗散和黏性耗散。线性频率估计对展向频率趋近零保持一致，这是适配无界展向方向的关键。非线性混合项通过高阶因子与低阶因子的配对估计控制，避免所有项都按同一个最坏尺度估算。作者最终保持两组不同大小的能量层级：一组为 $O(\varepsilon)$ ，另一组为 $O(\varepsilon/\nu)$ ，与 roll/nonzero 模态和 lift-up 放大的 streak 尺度相匹配。初值只有 $$H^2$$ 正则性时，先用正时间平滑在 $$t=1$$ 初始化较高阶能量，再按低阶、高阶、streak 的顺序完成 bootstrap 闭合。

4. 与之前工作的比较

最直接的比较对象是 Wei–Zhang 在 $\mathbb T\times\mathbb R\times\mathbb T$ 上的三维 Couette 稳定性结果，以及 Chen–Wei–Zhang 在有限通道 $\mathbb T\times[-1,1]\times\mathbb T$ 、无滑移边界条件下的结果。两者均达到 $O(\nu)$ 的扰动阈值；本文保持相同的黏性指数，但把展向从周期方向改为整条实线。这个几何变化移除了横向频谱间隙，使零模态低频扩散可能任意慢，因而不能直接沿用周期展向下的零模态耗散论证。本文的推进在于明确添加初始 roll 的临界 $\dot B^{-2}_{2,\infty}$ 控制，并证明该低频条件在非线性演化中可传播，从而在没有谱隙时仍达到 $O(\nu)$ 阈值。代价是稳定性假设不再只有全三维 $$H^2$$ 小性，而需额外控制流向平均横向速度的低频；该条件是充分的非线性判据，作者没有声称它对所有非线性稳定数据都是必要条件。

阈值尺度与线性 lift-up 相符。初始大小为 $\varepsilon$ 的适当 roll，在频率远离零的固定频带上经过 $t\sim\nu^{-1}$ 的演化，可产生大小约 $\varepsilon/\nu$ 的 streak，因此要在固定小邻域内稳定， $\varepsilon\lesssim\nu$ 是自然尺度。本文的贡献不是改善这个线性限制，而是在无界展向低频存在时通过端点条件达到该自然尺度。需要区分的是，临界性直接针对 $$u^2_0$$ 的线性 lift-up；作者将 $\dot B^{-2}_{2,\infty}$ 条件施加于整个 $u^\perp_0$ ，是为了控制完整 roll 方程及其二次应力，而不只是因为 $$u^3_0$$ 在线性系统中直接驱动 streak。

与更早的三维周期展向结果相比，Bedrossian–Germain–Masmoudi 的 Sobolev 稳定性阈值为 $O(\nu^{3/2})$ （在相应高阶 Sobolev 正则性下），其 Gevrey 结果达到 $O(\nu)$ ；Wei–Zhang 随后在 $$H^2$$ 框架下达到 $O(\nu)$ 。本文同样以 $$H^2$$ 控制三维初始扰动，保持 $O(\nu)$ 阈值，但额外处理了整条展向轴导致的零模态低频问题。Chen–Wei–Zhang 的有限通道结果也有边界和谱结构上的差别，不能简单视为本文低频条件的替代。本文没有改变这些周期或有限通道结果中的阈值指数，而是扩展了适用几何，并以额外的端点低频假设弥补谱隙缺失。

论文还提及整平面或其他无界区域上的相关工作。例如摘要与引言指出， $\mathbb R^2$ 相关方向上的无界域结果需要附加低频控制；但从提供的文本无法确认所引工作的精确假设、阈值及其与本文 Besov 条件的严格包含关系，因此不宜据此声称本文在所有无界域结果中更强。本文的明确特点是区域为 $\mathbb T\times\mathbb R^2$ ，流向仍周期，因此非零流向模态保有离散频率间隙；困难集中在流向零模态的横向低频，而不是三个方向都无界时的全面低频退化。作者还特别指出，旋转、稳定分层等能改变零模态动力学的色散或耦合机制不在本文模型中；本文在没有这些额外物理机制时直接控制低频 lift-up。

结论的范围也应准确理解。非零流向模态具有 $e^{-c\nu^{1/3}t}$ 的增强耗散，而 roll 只有 $\|u^\perp_0(t)\|_2\lesssim\varepsilon(1+\nu t)^{-1}$ 的代数衰减；streak 给出正时间高阶有界性而不保证趋零。尤其作者明确说明，临界 Besov 条件不能推出线性 lift-up 项必然衰减，因而本文证明的是全局稳定和分模态衰减/可积估计，不是整个扰动最终衰减到零的渐近稳定结论。对壁法向非零模态，结果是黏性一致的时间积分 inviscid-damping 估计，而非所未证明的逐点非线性幂率。文中进一步给出 $\dot H^{-2}$ 和空间加权初值条件作为 Besov 假设的充分替代，但一般 $\dot H^{-2}$ 并不能承载演化中出现的任意 $$L^1$$ 张量散度；所以这些替代条件主要用于验证初值，演化估计仍依赖端点 Besov 空间。

4. 三维环面 $\mathbb T^3$ 上 Navier–Stokes 方程的 Hartman–Grobman 定理

A Hartman-Grobman Theorem for the Navier-Stokes Equation on $\mathbb T^3$

作者: Rongchang Liu, Kening Lu, Lin Shi

类别: math.AP, math.DS

发布日期: 2026-10-06

arXiv ID: 2610.08058v1

摘要翻译: 我们证明了三维环面上无外力 Navier–Stokes 方程在零解处的局部 Hartman–Grobman 定理：对每个 $$s>7/2$$ ，该定理在 $$H^s$$ 的强拓扑下成立。更准确地说，我们构造了零点的两个正向不变开邻域之间的一个同胚，使 Navier–Stokes 半流与热半群共轭。

Comment: 无

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

1. 研究对象和背景

本文研究三维环面上无外力、不可压缩 Navier–Stokes 方程在零平衡点附近的非线性动力学, 问题是能否通过一个局部拓扑坐标变换, 将完整的非线性半流化为其线性化热半群, 且保留原有时间参数。方程为

\partial_t u+(u\cdot\nabla)u+\nabla p=\nu\Delta u,\qquad \nabla\cdot u=0,\qquad \nu>0,

其中 $$u(t,x)$$ 是速度场, $$p(t,x)$$ 是压力。空间为 $\mathbb T^3=(\mathbb R/2\pi\mathbb Z)^3$ , 初值取自零均值、散度为零的 Sobolev 空间

X=H^s_\sigma(\mathbb T^3;\mathbb R^3),\qquad s>\frac72.

零均值条件排除了 Stokes 算子的常数模, 因而线性化在该空间上严格耗散。令 $A=-\Delta$ , $$P$$ 为 Leray 投影, 则方程可写为

\dot u=-\nu Au-B(u,u),\qquad B(f,g)=P((f\cdot\nabla)g),

而在 $$u=0$$ 处的线性化是 $\dot u=-\nu Au$ , 对应热半群 $e^{-\nu At}$ 。对于充分小的初值, 文中使用标准小数据理论得到全局强解及连续半流 $$S_t$$ 。这里的目标并非只证明解趋于零或估计其衰减, 而是寻找零点邻域间的同胚 $$H$$ , 使每条局部非线性轨道在一个统一坐标变换下与热方程轨道逐时对应。

这一问题与有限维经典 Hartman–Grobman 定理有本质差别。有限维双曲平衡点附近的流可逆, 但抛物型方程的正时间演化具有平滑效应, 通常不能在原相空间局部可逆。因此不能直接把非线性轨道向过去延拓, 也不能直接照搬可逆动力系统中的线性化证明。另一方面, Navier–Stokes 非线性含有一阶导数损失: 在本文记号下 $B:X\times X\to H^{s-1}_\sigma$ , 而非 $$X$$ 。简单的 Fourier 截断也不能给出精确的有限维因子, 因为被截去的高频模会通过非线性影响未来低频模。三维 Stokes 谱由 $$|k|^2$$ 给出, 相邻谱层并无本文所需的统一大间隙, 因而以谱间隙吸收导数损失的既有策略不能直接适用。

作者转而利用 Foias–Saut 渐近正规化所编码的长时间渐近信息。不同初值的解可以在某一衰减尺度以下具有相同渐近行为, 由此形成嵌套的不变等价类和有限维渐近因子。对衰减级别 $$L$$ , 文中将关系刻画为

u\sim_L v\quad\Longleftrightarrow\quad e^{\nu Lt}\|S_tu-S_tv\|_{H^s}\longrightarrow0\quad(t\to\infty),

其含义是两条轨道在衰减速率不慢于 $e^{-\nu Lt}$ 的渐近层面一致, 差异只留在更快衰减的部分。该关系在半流下保持, 因此相应商空间具有精确的因子动力学。它与普通低频投影不同: 后者一般不与非线性演化相容。

然而, 渐近系数本身即使能区分初值, 也不自动给出强 $$H^s$$ 拓扑下的坐标图, 更不控制逆映射和无限多个高阶因子的尾部。本文试图填补的空白正是: 在有平滑、不可逆性和一阶导数损失的三维 Navier–Stokes 半流中, 如何从精确渐近因子恢复 Sobolev 几何, 在有限因子上构造共轭, 并将其组装成完整相空间的局部同胚。摘要和引言称最终组装依靠散度自由输运结构的谱估计, 且不要求额外谱间隙条件。

2. 主要定理或主要结果

主要结论是定理 2.1。固定 $\nu>0$ 及 $$s>7/2$$ , 在零均值散度自由空间 $X=H^s_\sigma(\mathbb T^3;\mathbb R^3)$ 中, 存在零点的开邻域 $U,V\subset X$ , 使得 $$U$$ 对 Navier–Stokes 半流 $$S_t$$ 正向不变, $$V$$ 对热半群 $e^{-\nu At}$ 正向不变, 并存在强 $$H^s$$ 拓扑下的同胚 $H:U\to V$ , 满足

H(0)=0,\qquad H(S_tu)=e^{-\nu At}H(u),\qquad u\in U,\ t\ge0.

这里 $A=-\Delta$ 限制在零均值散度自由场上。结论是完整半流的局部拓扑共轭, 不是仅有轨道衰减率相同、渐近等价或某个截断系统的共轭; 等式对所有 $t\ge0$ 成立, 因而共轭不改变时间参数。邻域只要求是零点附近的开集且正向不变, 文中摘要与定理陈述没有给出可显式计算的半径或共轭的定量模连续性界。

定理前后的分析还给出构造所需的小数据估计。对文中选取的充分小半径 $r_s\le \rho_s\nu$ , 实解满足指数衰减及耗散积分控制, 典型形式为

\|S_tu_0\|_s\le C_s e^{-\nu t}\|u_0\|_s,

以及

\int_t^\infty\bigl(\|S_qu_0\|_s^2+\|S_qu_0\|_{s+1}^2\bigr)\,dq \le C_s\nu^{-1}\|S_tu_0\|_s^2.

还证明了对应的一阶时间积分估计, 并在复化的小球上建立解映射的全纯性及正时间平滑性。这些是渐近正规化、纤维几何和变分方程估计的分析基础, 不是定理 2.1 的替代结论。文中另用 Dirichlet 商方法处理线性化变分方程: 在扰动系数满足相应 $$L^2$$ 可积条件时, 得到有限时间反向唯一性; 对非零解则有形如 $\|v(t)\|_s\ge c e^{-C(t-t_0)}$ 的指数下界。这些结论用于证明变分演化的单射和稠密值域性质。

构造层面, Foias–Saut 渐近系数 $$W_n(u)$$ 定义嵌套的有限维精确因子。作者证明这些因子在共同小邻域上具有所需的子浸没和分离切向量性质, 再以各纤维上的最小 $$H^s$$ 范数代表元恢复 Hilbert 空间几何。由此形成有限维 Sobolev 块坐标 $$z_j$$ , 使每个有限前缀具有自治动力学, 形式上为

\dot z_j=V_j(z_1,\ldots,z_j).

在每个有限维因子上构造边界—时间共轭, 再以谱估计控制有限层共轭及其逆的块尾, 最终组装为强 $$H^s$$ 拓扑下的全空间同胚。所给摘录没有提供后续定理的完整陈述, 因此无法从中核实各层谱估计的精确常数、误差阶或显式尾部衰减率。

定理 2.1 还蕴含一个关于长期动力学的结论: 对任意具有零空间均值的全局强解, 若已知其经典地满足 $u(t)\to0$ 于 $$H^s$$ , 则存在依赖于该解的 $T\ge0$ , 使 $u(T)\in U$ 。此后

u(t)=H^{-1}\!\left(e^{-\nu A(t-T)}H(u(T))\right),\qquad t\ge T.

这表示同一个局部同胚识别所有进入 $$U$$ 的全局强解的最终演化; 进入时间 $$T$$ 依解而变, 共轭本身则是统一的。

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

证明的主要困难是三个要求必须同时满足: 正时间抛物半流不可逆, 非线性输运项造成一阶导数损失, 而有限维共轭在维数增加时可能严重扭曲高频坐标, 使映射或逆映射不再保持 $$H^s$$ 的平方可和性。作者的策略不是对 Fourier 模直接截断, 而是先把长时间渐近信息变成精确因子, 再逐层恢复相空间几何并控制无穷维组装。

第一步是建立小解、变分演化及渐近因子的共同分析框架。投影后的方程为 $\dot u=-\nu Au-B(u,u)$ 。输运项的主部具有散度自由结构, 在 Sobolev 能量估计中可利用输运抵消; 通过交换子估计控制 $J_s=(I+A)^{s/2}$ 与输运算子的对易子, 得到 $$H^s$$ 能量及线性化算子的估计。小数据条件按黏性尺度取为 $r_s\lesssim\nu$ , 从而线性耗散吸收非线性项, 并得到解的指数衰减与时间积分控制。后者保证沿轨道的系数可积, 是处理渐近展开及伴随方程的关键, 而不仅是常规的局部适定性结论。

对线性化变分方程, 作者采用 Dirichlet 商

q(t)=\frac{\|A^{1/2}v(t)\|_s^2}{\|v(t)\|_s^2}

控制解的频率尺度。若扰动满足 $\|D(t)h\|_s\le\alpha(t)\|A^{1/2}h\|_s+\beta(t)\|h\|_s$ , 则导出 $$q$$ 的微分不等式以及 $\log\|v\|_s$ 的下界。结合 $\alpha,\beta\in L^2$ , 可排除非零变分解在有限时间内消失, 并得到其非零解的指数下界。该论证进一步给出变分演化的单射性; 对伴随演化使用同类反向唯一性, 则推出演化算子的值域稠密。这些性质支撑渐近因子的正则几何与纤维切空间分析。

第二步是以 Foias–Saut 展开构造精确而非截断的渐近因子。渐近系数 $$W_n$$ 记录解在对应衰减层次上的信息, 其前缀 $W_{\le L}$ 决定商掉更快衰减差异后的状态。因子之间有自然投影, 等价类构成嵌套的不变纤维。普通投影 $P_{\le N}u$ 不具备这种自治性, 因为高频状态可以影响未来低频演化; 渐近因子的定义则直接使用轨道未来, 因而保留精确的半流因子关系。文中以 Foias–Saut 的递归多项式展开表示这些系数, 并为余项及其对初值的导数建立共同复球上的控制, 以支持共同邻域上的纤维构造。

第三步是解决渐近标签与 Sobolev 拓扑不等价的问题。对每个渐近纤维 $$[u]_L$$ , 作者选取最小范数代表

\Pi_{\le L}u=\operatorname*{argmin}_{v\in[u]_L}\|v\|_{H^s}.

纤维是弯曲的, 所以最小点的存在还不足够, 需要唯一性和不同层之间的兼容性。核心几何估计是控制纤维的第二基本形式: 将差分方程中的二次余项与相应伴随法向量配对时, 线性演化项相消, 剩余项呈现两个向前衰减因子与一个向后增长的伴随因子的结构。该抵消使纤维曲率可控, 在充分小邻域内保证 $$H^s$$ 范数平方沿纤维严格凸, 从而得到唯一最小代表。作者随后使用依赖参数的 Morse 引理, 把逐层增加的最小范数代价写成块范数:

\|z_j\|^2=\|\Pi_{\le L_j}u\|_{H^s}^2-\|\Pi_{\le L_{j-1}}u\|_{H^s}^2, \qquad L_j\uparrow\infty.

因此有限块范数望远镜相加趋于 $\|u\|_{H^s}^2$ , 而投影误差由后续块的尾和控制。这个步骤的实质是证明渐近因子不仅区分状态, 还以可逆方式编码强 Sobolev 几何。

第四步是在这些坐标中利用因子结构实施边界—时间方案。每个有限前缀 $z_{\le j}$ 对应一个自治有限维因子系统。小球边界上的径向估计保证向量场严格向内, 因而非零点的有限维反向轨道可追溯到边界交点。记录边界交点及其到当前点的经过时间, 再从相同边界数据按参考流演化相同时间, 可定义该层的拓扑共轭。对第 $$j$$ 块使用前 $$j$$ 个因子的边界交点和时间, 使增加更高块时不改变已经定义的前缀坐标。Navier–Stokes 因子与热流因子都分别如此共轭到共同的对角参考流, 最后将两套有限层坐标复合。

最后也是无穷维中真正困难的一步, 是证明有限层同胚及其逆不会随层数积累过大的畸变。作者利用散度自由输运结构, 对归一化变分方程中的斜对称输运部分作非共振与近共振分解: 非共振部分通过共轭或消去处理, 对易子估计控制其代价; 剩余近共振部分对径向运动的影响较小。进一步的谱定位估计控制因子纤维法向方向以及快速衰减切向方向在谱块中的分布, 从而估计各块共轭和逆共轭的尾部贡献, 保持平方可和性。摘要明确强调该组装不依赖额外谱间隙假设; 但所给摘录尚未展示第 5 至第 7 节的具体谱不等式, 因而无法进一步可靠列出其精确指数、常数及误差阶。

4. 与之前工作的比较

与经典有限维 Hartman–Grobman 理论相比, 本文处理的是抛物型半流而非局部可逆流。在有限维双曲平衡点理论以及 Banach 空间微分同胚理论中, 反向轨道和双曲分解是基本工具; 但 Navier–Stokes 的正时间映射因平滑效应不局部可逆, 因而本文不能直接应用这些结果。本文的结论是局部拓扑共轭, 作者也明确没有证明共轭为 $$C^1$$ ; 引言将共轭能否具有 $$C^1$$ 正则性列为开放问题。

就抛物方程而言, 文中回顾 Aulbach–Garay [3] 对非可逆映射的部分线性化, 指出其方法需要额外谱条件; Lu [27] 与 Bates–Lu [5] 对若干一维抛物方程建立完整相空间共轭, 其方法利用不断增大的谱间隙控制非线性耦合。本文的推进在于针对三维 Navier–Stokes 这类有导数损失的问题构造完整局部共轭, 而非仅得到部分线性化, 并避免把相邻谱间隙作为额外假设。相应代价是证明需引入嵌套渐近因子、非线性纤维几何、谱块估计和边界—时间组装, 技术复杂度显著提高。

作者还将本文与其关于二维半线性抛物方程的工作 [25]、以及无导数损失反应扩散方程的工作 [26] 对照。[25] 利用平面谱渐近选择谱带, 由带边界间隙补偿导数损失; 本文指出 Stokes 谱缺少该工作所依赖的分离结构, 因而原有选带策略不适用于三维 Navier–Stokes。[26] 针对无导数损失非线性, 使用固定 dyadic 谱带, 不要求相邻谱间隙下界或谱重数一致有界, 并可处理任意维环面上的反应扩散方程。但 Navier–Stokes 的对流项损失一个导数, 不能直接沿用该构造。本文的区别不只是更换谱带, 而是依靠散度自由输运的结构性抵消与纤维法向、快速衰减切向的谱定位来控制尾部畸变。

相较于直接 Fourier 截断, 本文的有限维对象也有不同性质。低频投影通常不是半流因子, 因为未来低频分量会受到省略高频分量的影响; Foias–Saut 渐近因子则按轨道的衰减层次取商, 是与半流严格相容的有限维动力系统。其优势是获得精确自治的有限维因子, 代价是必须证明这些长时间渐近数据的纤维在共同邻域中具有良好几何, 并恢复其对强 $$H^s$$ 拓扑的控制, 否则因子坐标的单射性本身并不能保证连续逆映射。

结果范围方面, 定理对每个 $$s>7/2$$ 及 $\nu>0$ 成立, 相空间是三维环面上的零均值散度自由速度场, 结论限于零解附近的局部拓扑共轭。作者指出所有零空间均值的全局强解经典地趋于零, 所以每条此类轨道最终进入共轭邻域, 但这并不把定理提升为全相空间全局共轭, 也不意味着初始时刻就处于该邻域。文中未声称对有外力、非零均值、边界域或其他流体模型成立。也未证明共轭的可微性、定量 Hölder/Lipschitz 界或显式邻域大小。

关于前人工作的具体定理假设及其最优性, 所给 PDF 摘录主要提供引言中的定性比较, 未展示这些工作的完整陈述; 对它们之间更细的参数范围和技术常数, 从提供内容中无法确认。类似地, 本文谱估计的精确形式在所给内容中尚未出现, 因此可以确认作者声称不需要额外谱间隙条件, 但不能据此补写其具体定量代价或最优性结论。

5. 纳维–斯托克斯方程的翻译失真：为何对 AI 自动形式化结果进行 Lean 验证并不能保证自然语言证明正确

Navier-Stokes lost in translation: Why Lean verification of AI autoformalisation does not guarantee correct natural language proofs

作者: Alexander Bastounis, Fabian Circelli, Anders C. Hansen

类别: math.AP, cs.AI, math.LO

发布日期: 2026-10-06

arXiv ID: 2610.08144v1

摘要翻译: 自动形式化日益被用于验证数学文本，包括 AI 生成的文本，例如 OpenAI 宣布的纳维–斯托克斯方程解爆破解。在这一过程中，AI 系统将文本从自然语言（NL）翻译成 Lean 等形式语言。完成翻译后，形式语言所表达的论证便很容易通过机器进行验证。本文旨在说明，为何这一过程并不能让人对原始 NL 论证有任何信心，原因在于，要实现语义忠实的翻译，存在种种困难。特别是，我们强调，为了提供语义忠实的翻译，必须解决数学 NL 文本中的歧义，而这一问题在可解性复杂度指标（SCI）层级/算术层级中可任意高（SCI $= \infty$ ）。因此，非正式地说，提供语义忠实的 AI 自动形式化，比包括停机问题在内的任何计算问题都更难（停机问题的 SCI $$= 1$$ ）。为展示这一结果的影响，我们给出了若干实际案例，说明 AI 将 NL 陈述和证明错误地翻译成 Lean，导致 NL 证明与其 Lean“验证”之间不匹配。其中包括 OpenAI 宣布的纳维–斯托克斯证明。特别地，我们表明，形式化后的 Lean 证明并不对应于关于纳维–斯托克斯方程解爆破的 NL 证明。

Comment: 25 pages, 4 Figures

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

1. 研究对象和背景

本文讨论的核心对象不是 Navier-Stokes 方程的解本身, 而是数学自然语言文本向 Lean 等证明助理语言的自动形式化, 以及形式化结果能否忠实表达原文语义。相关 PDE 背景是三维不可压缩 Navier-Stokes 方程, 可写为

\begin{cases} \partial_t u+(u\cdot\nabla)u-\nu\Delta u+\nabla p=0,\\ \nabla\cdot u=0,\\ u(0,x)=u_0(x), \end{cases}

其中 $$u(t,x)$$ 是速度场, $$p(t,x)$$ 是压力, $\nu>0$ 是黏性系数。该方程的三维正则性问题历史悠久, 与 Leray 的弱解理论及 Clay 千禧年问题相关。本文所针对的具体背景是 OpenAI 宣布证明三维 Navier-Stokes 解有限时间爆破, 并在 GitHub 仓库中提供相应 Lean 形式化。这里的爆破问题属于 PDE 中极具挑战性的全局正则性问题; 但本文作者明确限定讨论范围: 他们不判断 OpenAI 的自然语言证明是否正确, 而是检查所提供的 Lean 形式化是否是该自然语言文本的语义忠实翻译。

这一限定很重要。形式证明助理验证的是某个精确定义的命题及其形式推导。若自然语言证明被转换为 Lean 代码, Lean 编译成功只能推出 Lean 中实际陈述的命题可由该代码证明; 它本身不能推出原文的命题、引理、推理步骤均被保留。作者区分两类自动形式化问题。第一类是在目标定理陈述已被正确翻译的前提下, 把一份给定的自然语言证明转为该定理的 Lean 证明, 并以代码编译通过作为成功标准。第二类则要求对文本中的定义、命题、证明及其依赖关系作语义保真的逐项翻译; 若要借此验证原文证明, 还必须确认形式证明与书面论证相对应。前一种任务允许系统在不改变最终定理的条件下改用另一套证明, 甚至修补原证明中的错误; 后一种任务则必须保存源文本的数学内容和论证关系。

因此, 论文研究的问题可概括为: 给定数学自然语言文本 $$S$$ 和形式化结果 $$F$$ , 何种条件足以保证 $$F$$ 忠实表达 $$S$$ , 而不是仅仅证明一个相同、较弱或相关的命题? 数学语言中省略的量词范围、术语定义、隐含假设、证明中“因此”的逻辑连接, 以及不同符号和范数的约定, 都可能使翻译出现语义偏移。要检验忠实性, 不能只检查 $$F$$ 是否通过 Lean, 还要逐一核对自然语言陈述、形式化陈述、假设、结论与证明步骤之间的对应。本文指出, 消解数学自然语言歧义这一忠实翻译所必需的工作, 在 Solvability Complexity Index (SCI) / 算术层级意义下可处于任意高层, 即摘要所称 $\mathrm{SCI}=\infty$ ; 作者据此强调, 该问题在一般情形下比停机问题更难以算法化解决。提供的节选没有给出这一复杂性结论的完整定义、定理陈述和证明, 因而其精确技术假设与层级构造无法从当前材料中进一步确认。

论文试图填补的空白不是再提供一个 PDE 正则性定理, 而是澄清“形式化代码可机械验证”与“自然语言证明已获验证”之间的逻辑鸿沟。作者用初等例子说明这种鸿沟如何产生, 再逐项比较 OpenAI Navier-Stokes 论文中的两个估计及 Lean 对应声明, 展示结论强度、正则性损失、证明路线和时间量词均可能发生变化。就 PDE 阅读者而言, 关键警示是: Lean 机器内核的可靠性针对形式系统内部的证明对象, 而数学文本到该对象的语义对应本身是另一个独立且可能极困难的数学问题。

2. 主要定理或主要结果

文章的主要结论是关于自动形式化可靠性的分析结论和具体反例, 而非 Navier-Stokes 方程的存在、正则性或爆破定理。

作者区分两种任务: (i) 在目标定理已经忠实形式化的前提下, 让 AI 将一份自然语言证明转成 Lean 代码并使之通过编译; (ii) 将原始数学文本中的定义、定理、证明及其依赖整体作语义保真翻译。本文指出, (i) 的成功不蕴含 (ii), 也不蕴含原自然语言证明正确。Lean 代码即使不使用 sorry 或额外公理而成功编译, 仍可能证明了另一条论证所得的结论, 或只证明了比原文更弱的结论。
两个初等例子展示上述逻辑差别。对 $$p(x)=x^3-x^2-x+1$$ , 原自然语言证明错误地声称根及重数给出 $$p(x)=(x-1)(x+1)^2$$ ; 聊天机器人却返回了一个利用正确因式分解的可编译 Lean 证明。因此, 形式证明正确不能反推所提交的自然语言证明正确。对实对称矩阵 $$A$$ , 原文以正交对角化证明 $\operatorname{tr}(A^2)\geq0$ , 而 Lean 代码采用 $a_{ij}=a_{ji}$ , 得到 $\sum_{i,j}a_{ij}a_{ji}=\sum_{i,j}a_{ij}^2\geq0$ 。两个证明都正确, 但形式证明并未翻译原文的证明路线。
对 OpenAI Navier-Stokes 论文的 Lemma 8.6, 原文声称对零均值光滑函数 $$F$$ 有反演估计

\|N^{-1}F\|_{C_y^m}\leq C_m\|F\|_{C_y^{m+4}},

其中 $y\in\mathbb T^2$ , $N=v_t\cdot\nabla_y$ , $v_t=(\sqrt2-1,1)$ , 且 $N^{-1}F$ 指满足 $N\phi=F$ 的唯一零均值光滑解。对应 Lean 结果所需输入导数为 $$m+5$$ , 其可比较形式为

\|N^{-1}F\|_{C_y^m}\leq \widehat K_m\|F\|_{C_y^{m+5}}.

Lean 代码中的导数假设以 $$w.length+5$$ 控制输出中长度为 $$w.length$$ 的导数; 对所有 $|w|\leq m$ 取上确界后, 仍有额外一个输入导数。作者指出, 原文使用除数估计 $|v_t\cdot k|^{-1}\lesssim 1+|k|$ 并结合二维频率级数中 $(1+|k|)^{-3}$ 的可和性, 因而只需 $$m+4$$ 阶导数; Lean 证明通过另一套系数半范数估计得到 $\sum_{k\in\mathbb Z^2}(1+|k_1|+|k_2|)^{-4}$ , 从而要求 $$m+5$$ 阶导数。该 Lean 结论较弱; 文中还提醒, 部分相关 Lean 声明采用混合导数, 与原文只涉及环面变量导数的范数并非完全同型。

对原文式 (10.19) 的压力通量估计, 作者指出 Lean 对应结果在结论和论证上均有差异。令 $R\geq1$ , 原文的估计为

\left|\int \pi w\cdot\nabla\chi_R\right| \leq C_T\left[\frac{(B_R+1)B_R^{1/2}}{R}+R^{-3/4}B_R^{3/4}\right], \quad\text{对几乎处处 }t\in(0,T),

而被比较的 Lean 声明给出

\left|\int \pi w\cdot\nabla\chi_R\right| \leq C_T\left[(B_R^{1/2}+1)\left(\frac{A_R}{R}+\frac{1}{R^2}\right)+R^{-7/4}B_R^{3/4}\right], \quad\forall t\in(0,T).

这里 $A_R=(\int\chi_R|\nabla w|^2)^{1/2}$ , $B_R=\|\phi_R^4w\|_{L^6}$ 。Lean 估计额外依赖 $$A_R$$ , 并且 $$R$$ 的幂次及 $$B_R$$ 项的组合也不同; 原文中另一条关于 $$B_R$$ 的估计不能据此消去 $$A_R$$ 并恢复原式, 所以两条界不能简单视为等价。Lean 证明还把原文的“几乎处处时间”结论升级为每个 $t\in(0,T)$ , 这是其显式建立时间连续性后的结果, 但不消除空间尺度表达式上的差异。

摘要给出的复杂性结论是: 为实现语义保真的自动形式化, 必须处理自然语言数学文本中的歧义; 这一消歧问题在 SCI / 算术层级中可任意高, 摘要以 $\mathrm{SCI}=\infty$ 表述, 并将其与停机问题的 $\mathrm{SCI}=1$ 作比较。由于当前提供的 PDF 节选未包含相应章节和定理证明, 本文这里只能报告作者摘要中的结论, 无法补充其精确定义、归约构造或适用模型。作者明确不据此断言 OpenAI 的自然语言 Navier-Stokes 证明错误; 他们证明和展示的是 Lean 形式化与该自然语言证明之间存在未被编译检查排除的语义不一致。

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

本文的技术重点在于把“证明对象正确”与“翻译关系忠实”分开处理, 再以可逐行核验的数学差异说明后一问题不能由 Lean 编译自动解决。

首先, 作者建立了两层验证框架。形式系统内的检查是: Lean 是否接受给定声明及其证明项; 语义层面的检查则是: 该声明是否对应自然语言文本的原命题, 证明项是否保留原文假设、结论和论证步骤。只做前一层检查时, 系统只需找到任何能证明形式化目标的项, 不必忠实复现用户提交的论证。若原文推理有错, 系统可以悄悄替换成正确证明; 若原文论证正确, 系统也可以省略关键中间步骤并另走一条证明路线。这个区分解释了为什么“编译通过”不是自然语言证明的正确性证书。

其次, 作者使用初等反例把抽象的语义偏移具体化。多项式例中, 自然语言证明使用错误根重数和错误因式分解, 但生成的 Lean 代码实际依据正确的因式分解证明目标; 这不是对错误论证的形式化, 而是对目标命题另行证明。矩阵例中, 自然语言论证先对称矩阵作正交对角化, 再将迹写成特征值平方和; Lean 代码则留在原矩阵坐标, 利用对称性把 $a_{ij}a_{ji}$ 化为 $a_{ij}^2$ 。逐项比较证明依赖关系, 能区分“同一结论的不同证明”和“源证明已被验证”这两个不同主张。

第三, 对环面方向算子反演估计, 作者不是只比较定理名称, 而是追踪输入正则性如何进入输出范数。对 $N=v_t\cdot\nabla_y$ 作 Fourier 展开时, 零均值条件排除 $$k=0$$ 模式并固定反演解的常数自由度; 除数估计控制非零频率上的乘子 $|v_t\cdot k|^{-1}$ 。原文依靠二维格点级数的可和性, 将输入导数损失控制在四阶。作者进一步检查 Lean 定理 norm_derivativeWord_inverse_le 的假设, 发现其对长度为 $$m$$ 的导数词要求 $$m+5$$ 阶的坐标导数界; 追踪到 coefficient_seminorm_bound 后, 识别出相应估计使用的权重级数是 $\sum_{k\in\mathbb Z^2}(1+|k_1|+|k_2|)^{-4}$ 。因此差别不只是常数记号, 而是证明所需的正则性损失由四阶变为五阶。作者也谨慎指出, 某些 Lean 声明使用混合导数, 不能未经检查就与原文 $$C_y^m$$ 范数直接等同。

第四, 对压力通量估计, 作者沿着形式化代码追溯其调用链, 将顶层 Lean 定理与 NL 原文式 (10.19) 并列, 再辨析估计依赖的分析工具。原文通过 Riesz 变换在 $L^{3/2}$ 上的有界性处理压力项; Lean 路线则先借 Sobolev 嵌入移至 $$L^2$$ , 再使用 Riesz 变换在 $$L^2$$ 上的有界性。该替换引入额外导数, 反映为 Lean 估计中出现局部梯度量 $$A_R$$ 。作者还指出两者使用不同的 Hölder 指数, 因而最终的 $$B_R$$ 幂次及 $$R$$ 的衰减率不同。除此之外, Lean 证明显式证明通量积分关于时间连续, 从而把几乎处处成立的结论升级为所有内部时间成立。这里作者没有把后一项视为单纯错误, 而是说明它属于证明策略造成的真实差异; 与此同时, 更强的时间量词并不补偿右端表达式和依赖量与原文不一致。

第五, 论文摘要提出以 SCI / 算术层级分析语义消歧的计算复杂性, 并称其层级可任意高。其论证应当用于说明一般语义忠实翻译并非一个可由简单算法保证解决的任务, 而不是声称某个具体 AI 系统在某一次翻译中必然失败。当前节选止于第 10 页, 未包含复杂性章节的正式技术发展; 因此无法可靠复述具体编码、归约或定理条件, 也不应仅凭摘要扩展出更强的复杂性命题。

总体而言, 作者的实际分析框架是“声明对齐—假设核查—范数与量词比较—证明依赖追踪”。它使差异落在可检验的数学对象上: 导数损失阶数、所用频率级数、压力估计中的函数空间、局部梯度依赖、空间尺度幂次及时间量词。论文的创新点不是提出新的 Navier-Stokes 估计, 而是把自动形式化可能发生的语义替换落实为经过 Lean 声明和源文公式对照的案例, 并将其与一般消歧问题的复杂性联系起来。

4. 与之前工作的比较

就 PDE 文献而言, 文章把 OpenAI 所声称的三维 Navier-Stokes 有限时间爆破工作放在 Leray 1934 年以来的研究背景及 Clay 千禧年问题背景中, 但提供的节选没有展开这一 PDE 问题的技术文献综述, 也没有比较各类弱解、局部正则性或爆破判据。因此, 不能据此判断 OpenAI 的 PDE 论证相对于既有 Navier-Stokes 理论究竟推进了什么, 更不能把本文误读为对爆破定理的证明或反证。作者在第 2 页明确声明, 不对自然语言证明的正确性作结论, 只陈述 Lean 翻译存在不匹配。

本文真正的比较对象是两种自动形式化目标。相对于仅要求 Lean 代码通过编译的做法, 作者要求语义忠实翻译还必须维持定义、定理陈述、证明和依赖关系的对应。前者的验收标准是形式系统内可验证; 后者还要求审查形式化对象与自然语言源文之间的语义映射。因而本文并非否认 Lean 内核验证的价值, 而是限定其保证范围: 它能确保形式化声明在所用公理和定义下有形式证明, 不能独立确保声明就是原文所声称的命题, 也不能确保代码中的推理就是原文的推理。

与常见的“形式化一个正确定理”案例相比, 初等例子揭示两类不同风险。多项式例中, 原证明错误而 Lean 证明正确, 形式化过程事实上修正了论证; 矩阵例中, 两套证明均正确, 但 Lean 证明没有忠实复现给定论证。OpenAI 案例又增加第三类风险: 自然语言文本声称的结论比对应 Lean 声明强, 或两者的估计依赖和证明方法不同。故单纯报告形式化项目规模、编译状态或最终定理名称, 都不能替代逐项语义核对。

在具体分析推进上, 论文指出环面反演估计从原文的 $$m+4$$ 阶输入导数变为 Lean 的 $$m+5$$ 阶要求。这不是形式化后的常数变化, 而是正则性损失的实质变化; Lean 版本在该点更弱, 也不能据其成立反推原文较强估计已被验证。压力通量估计则不宜概括成简单的强弱关系: 原文界依赖 $$B_R$$ , Lean 界另含 $$A_R$$ , 且 $$R$$ 的幂和 $$B_R$$ 指数结构有别; 原文中关于 $$B_R$$ 的辅助控制并不能直接把 Lean 界化为原文界。另一方面, Lean 路线在时间变量上给出对所有 $t\in(0,T)$ 的陈述, 原文是几乎处处成立。作者把这种差异归因于时间连续性论证, 同时指出它与空间估计的差异是不同维度, 不应混为一谈。

这些例子也显示形式化证明的潜在代价。为了使证明助理中的估计易于实现, 证明可以改用不同的函数空间或不等式路线; 例如以 Sobolev 嵌入转到 $$L^2$$ 并使用该空间中的 Riesz 变换有界性, 可能引入额外导数及 $$A_R$$ 依赖。此类替代证明在 Lean 内部可以完全有效, 但若项目宣称是在形式化原文的特定估计与证明, 就必须明确标示这些变化, 并证明替代结果足以支撑原论证所需结论。本文没有主张所有更弱估计都无用, 而是强调它们不能未经额外推导就冒充原文的强估计。

在计算复杂性层面, 摘要将语义歧义消解与 SCI / 算术层级联系, 并称忠实自动形式化问题可达任意高层, 与 SCI 为 $$1$$ 的停机问题作比较。若该结论按作者所述成立, 它提示不存在对一般数学文本都能有效保证语义忠实的简单计算程序; 这与“对某个给定形式命题进行机械检查”并不矛盾, 因为后者的输入语义已经由人或其他过程固定。由于所给 PDF 内容没有包括该结果的正式陈述和证明, 具体与既有 SCI 文献、Hilbert 第十问题或停机问题的技术比较, 从提供内容中无法确认。

作者提出的实践立场是: 对 AI 生成并自动形式化的证明, 仍需像其他数学证明一样接受同行审查; 审查不仅要跑 Lean, 还要核对原文命题、形式化声明、假设、常数依赖、范数、时间量词和每个关键推理步骤。本文明确举出 OpenAI 仓库在特定提交 f9e8bc5 中的对应声明作为比较对象, 所以其代码观察依赖该版本; 后续代码若有修改, 需要重新核验。文章未证明 AI 自动形式化在所有实例中都不可信, 也未证明 OpenAI 的自然语言 Navier-Stokes 证明错误; 它证明并展示的是, 编译成功本身不足以排除翻译失真。

6. 二维 Navier–Stokes-Boussinesq 方程中 Poiseuille 流的稳定性阈值

Stability threshold of the Poiseuille flow in 2D Navier-Stokes Boussinesq equations

作者: Nangao Zhang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07598v1

摘要翻译: 我们研究二维 Navier–Stokes Boussinesq 方程中 Poiseuille 流在三个区域 $$$$ 上的转变阈值问题： $\mathbb{T} \times \mathbb{R}$ 、 $\mathbb{T} \times [0, \infty)$ 和 $\mathbb{T} \times [-1, 1]$ 。更准确地说，我们证明，若涡量和温度的初始扰动在各向异性 Sobolev 空间中的大小分别为 $O(\min\{ν^{\fr34},μ^{\fr34},ν^{-\fr14}μ\})$ 和 $O(\min\{ν^{\fr94},μ^{\fr94}\})$ ，则该系统的解对所有时间的大小均保持如此，并且非线性增强耗散以与 $O(\min\{ν^{\fr12},μ^{\fr12}\})$ 成正比的速率持续存在。

Comment: 无

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

1. 研究对象和背景

本文研究二维 Navier-Stokes-Boussinesq (NSB) 方程中 Poiseuille 剪切流的扰动稳定性及转捩阈值. 速度场 $$v=(v^1,v^2)$$ 满足不可压缩条件, 温度为 $\vartheta$ , 浮力沿竖直方向作用; 黏性系数为 $\nu>0$ , 热扩散系数为 $\mu>0$ . 研究区域为周期水平变量与竖直变量组成的三个条带型区域 $\mathbb T\times D$ , 其中 $D=\mathbb R$ , $[0,\infty)$ 或 $$[-1,1]$$ . 将重力常数归一化为 $$1$$ , 方程为

\begin{cases} \partial_t v+v\cdot\nabla v+\nabla\Pi-\nu\Delta v=\vartheta e_2,\\ \partial_t\vartheta+v\cdot\nabla\vartheta-\mu\Delta\vartheta=0,\\ \nabla\cdot v=0. \end{cases}

边界条件取为 $\partial_y v^1=2y$ , $$v^2=0$$ , $\vartheta=a$ 于 $y\in\partial D$ . 该系统描述温度差异造成的浮力与流体运动的耦合, 是对流、热输运等问题的基本模型; 从分析角度看, 它同时包含剪切混合、黏性耗散、热扩散以及速度与温度之间的耦合, 因而不能简单地将两个标量扩散方程分别处理.

相应的稳态为 Poiseuille 流 $$v_s=(y^2,0)$$ 、常温 $\vartheta_s=a$ , 以及相应压力. 令速度、温度和涡量扰动分别为 $$u=v-v_s$$ , $\theta=\vartheta-a$ 和 $\omega=\nabla^\perp\cdot u$ , 其中 $u=\nabla^\perp\phi$ , $\Delta\phi=\omega$ . 扰动方程为

\begin{cases} \partial_t\omega+y^2\partial_x\omega-2\partial_x\phi-\nu\Delta\omega+u\cdot\nabla\omega=\partial_x\theta,\\ \partial_t\theta+y^2\partial_x\theta+u\cdot\nabla\theta-\mu\Delta\theta=0,\\ u=\nabla^\perp\phi,\quad \Delta\phi=\omega. \end{cases}

这里第一行中的 $-2\partial_x\phi$ 来自基流涡量的梯度与扰动速度相互作用, $\partial_x\theta$ 则是浮力在涡量方程中的贡献. 扰动满足 $\omega=0$ , $\theta=0$ 于竖直边界; 速度边界条件对应于基流的边界条件. 初值记为 $(\omega_{\rm in},\theta_{\rm in})$ . 核心问题是在小黏性、小热扩散极限下, 初始扰动需要多小, 才能保证解全局保持受控, 并且剪切流产生的增强耗散不被非线性作用破坏.

研究此问题的重要性在于, 普通扩散给出的时间尺度通常分别为 $\nu^{-1}$ 与 $\mu^{-1}$ , 而剪切使不同水平频率持续产生竖直方向的细尺度, 可以带来快于纯扩散的增强耗散. Poiseuille 剪切率随 $$y$$ 变化, 在 $$y=0$$ 处退化, 其线性机制和 Couette 流的常剪切情形不同. 同时, 涡量与温度通过浮力耦合, 非线性速度又会输运二者, 因而既要捕捉剪切混合产生的耗散, 又要控制零水平频率与非零频率间的能量交换.

既有研究已分析了 NSB 方程在常态或 Couette 流附近的适定性、线性稳定性和非线性稳定性, 但 Poiseuille 流附近的结果较少, 尤其是整个区域 $\mathbb T\times\mathbb R$ 和半区域 $\mathbb T\times[0,\infty)$ . 对不含温度耦合的二维 Navier-Stokes 方程, 已有工作在全直条带及有限通道中建立了 Poiseuille 流的稳定性和增强耗散结果; 而对本文的 NSB 模型, Wang 在有限通道中证明过非线性稳定性, 但据本文介绍未得到增强耗散估计. 本文试图补足的正是耦合 NSB 模型在三个区域上的定量稳定性与增强耗散理论, 并为此给出基于能量估计而非直接套用半群估计的分析框架. 需要区分的是, 论文称其为转捩阈值问题, 但所提供的主要定理是初值充分小条件下的稳定性结论; 未证明大于该尺度必然不稳定, 因而不能将所得充分稳定尺度解释为已确定的尖锐阈值.

2. 主要定理或主要结果

主要结果是定理 1.1 和定理 1.2. 设 $$m>1/2$$ , $0<\nu,\mu\leq 1$ , 并记 $\kappa=\min\{\nu,\mu\}$ . 在 $D=\mathbb R$ 或 $[0,\infty)$ 时, 若初值满足下述各向异性能量小性条件, 则相应扰动解对所有 $t\geq0$ 保持控制, 且所有非零水平频率具有增强耗散衰减.

具体地, 涡量初值的条件为

\|u^1_{0,\rm in}\|+\|y\partial_y\omega_{0,\rm in}\|+\sum_{j=0}^1\|(\nu^{1/4}\partial_y)^j\langle\partial_x\rangle^{m-j/4}\omega_{\rm in}\|+\nu^{-1/4}\||\partial_x|^{m+1/4}(y\omega_{\rm in},\nabla\phi_{\rm in})\| \leq \varepsilon_0\min\{\nu^{3/4},\mu^{3/4},\nu^{-1/4}\mu\},

温度初值条件为

\|y\partial_y\theta_{0,\rm in}\|+\sum_{j=0}^1\|(\mu^{1/4}\partial_y)^j\langle\partial_x\rangle^{m-j/4}\theta_{\rm in}\|+\mu^{-1/4}\||\partial_x|^{m+1/4}y\theta_{\rm in}\| \leq \varepsilon_0\min\{\nu^{9/4},\mu^{9/4}\}.

其中 $\varepsilon_0>0$ 与 $\nu,\mu$ 无关, $\langle\partial_x\rangle$ 表示水平频率权重. 这些条件并非单一 Sobolev 范数的小性: 它们还要求零模的加权竖直正则性, 并对非零水平频率施加分数阶水平频率权重及含 $$y$$ 的矩条件. 由此, 对某个仅依赖普适常数的 $\delta_0>0$ , 定理给出

\sum_{j=0}^1\|(\nu^{1/4}\partial_y)^j\langle\partial_x\rangle^{m-j/4}\omega_{\ne0}(t)\|+\nu^{-1/4}\||\partial_x|^{m+1/4}(y\omega_{\ne0}(t),\nabla\phi_{\ne0}(t))\| \lesssim \varepsilon_0 e^{-\delta_0\kappa^{1/2}t},

以及

\sum_{j=0}^1\|(\mu^{1/4}\partial_y)^j\langle\partial_x\rangle^{m-j/4}\theta_{\ne0}(t)\|+\mu^{-1/4}\||\partial_x|^{m+1/4}y\theta_{\ne0}(t)\| \lesssim \varepsilon_0 e^{-\delta_0\kappa^{1/2}t}.

这里下标 $\ne0$ 表示去掉水平零频率. 因而增强耗散时间尺度为 $O(\kappa^{-1/2})$ , 衰减率为 $e^{-c\kappa^{1/2}t}$ ; 论文明确称该速率与相关 Poiseuille 流 Navier-Stokes 结果一致. 在无界区域和半无界区域, 定理 1.1 没有声称零模指数衰减, 这与竖直方向无 Poincare 不等式相符. 零模通过能量估计保持受控, 但不能据此改述为全解所有频率均按 $e^{-c\kappa^{1/2}t}$ 衰减.

在有限通道 $$D=[-1,1]$$ 上, 定理 1.2 的初值小性条件与定理 1.1 相同. 除上述非零模估计外, 由于可使用 Poincare 不等式, 零模也满足

\sum_{j=0}^1\|(\nu^{1/4}\partial_y)^j\omega_0(t)\|+\|y\partial_y\omega_0(t)\| \lesssim \varepsilon_0 e^{-\delta_0\kappa t},

\sum_{j=0}^1\|(\mu^{1/4}\partial_y)^j\theta_0(t)\|+\|y\partial_y\theta_0(t)\| \lesssim \varepsilon_0 e^{-\delta_0\kappa t}.

因此有限通道中零模的耗散率为 $O(\kappa)$ , 非零模仍具有 $O(\kappa^{1/2})$ 的增强耗散率. 论文的能量闭合还给出相应时间加权耗散积分有限, 但摘录中作为主定理陈述的结论是上述点态指数估计.

这些定理证明的是一组对 $\nu,\mu$ 一致的充分小性条件及其稳定性后果. 涡量阈值尺度为 $\min\{\nu^{3/4},\mu^{3/4},\nu^{-1/4}\mu\}$ , 温度阈值尺度为 $\min\{\nu^{9/4},\mu^{9/4}\}$ . 由于尚无相应不稳定性结论, 也未证明这些指数最优, 作者在备注 1.4 中明确指出转捩阈值可能并非尖锐, 最优阈值仍是开放问题.

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

证明的首要障碍是剪切输运、退化剪切点和涡量-温度耦合同时存在. 对单个线性剪切方程可以考虑半群或频率乘子, 但在这里涡量方程含有 $\partial_x\theta$ , 温度又通过非线性速度输运, 两个方程耗散系数分别为 $\nu$ 和 $\mu$ . 作者指出, 这种耦合使直接仿照纯 Navier-Stokes 方程的半群估计变得困难. 证明的主线因此不是把线性解算子分开估计, 而是构造同时反映剪切混合、扩散和空间矩的模态能量, 再将线性耗散余量用于吸收非线性项.

首先对水平变量 $$x$$ 作 Fourier 展开, 分别处理零模 $$k=0$$ 与非零模 $k\neq0$ . 这种拆分具有实质作用: 非零模感受到 $y^2\partial_x$ 的剪切输运并产生增强耗散; 零模不受该输运项影响, 在无界或半无界竖直区域也没有 Poincare 衰减, 其控制只能依靠扩散及非零模对零模的非线性强迫. 若将所有模态放在一个统一估计中, 便难以同时识别这两种不同的耗散机制.

对非零模, 作者为每个 $k\neq0$ 构造涡量能量 $E_\omega[\omega_k]$ 和温度能量 $E_\theta[\theta_k]$ . 其结构包括通常的 $$L^2$$ 项、带有 $\nu^{1/2}|k|^{-1/2}$ 或 $\mu^{1/2}|k|^{-1/2}$ 系数的竖直导数项、交叉项 $|k|^{-1}\operatorname{Re}\langle iky f_k,\partial_y f_k\rangle$ , 以及分别控制 $y\omega_k$ , $y\theta_k$ 和流函数梯度的矩项. 例如, 涡量能量中包含

\|\omega_k\|^2+\nu^{1/2}|k|^{-1/2}\|\partial_y\omega_k\|^2+5|k|^{-1}\operatorname{Re}\langle iky\omega_k,\partial_y\omega_k\rangle+100\nu^{-1/2}|k|^{1/2}(\|y\omega_k\|^2+2\|\nabla_k\phi_k\|^2).

温度能量使用对应的 $\mu$ 权重, 交叉项系数也相应调整. 交叉项不是形式上的修饰: 它把乘以 $$y$$ 的矩与竖直导数联系起来, 从而在对能量求导时提取剪切混合所产生的耗散. 对应耗散泛函含有 $\nu^{1/2}|k|^{1/2}\|\omega_k\|^2$ , $|k|\|y\omega_k\|^2$ 等项, 温度方程则有热扩散的对应项. 这类能量安排体现了 hypocoercivity 思路: 方程本身的扩散并不直接控制所有所需范数, 但剪切输运与导数/矩之间的耦合使能量获得额外衰减.

随后按 $|k|^{2m}$ 对模态能量求和, 并引入时间权重 $e^{2\delta_0\kappa^{1/2}t}$ . 该权重将目标增强耗散率直接编码进能量, 通过选择充分小的 $\delta_0$ 保证权重求导产生的正项可由模态耗散吸收. 证明中得到非零模能量不等式 (2.16)-(2.17); 其右端逐项记录涡量、温度、零模速度与各耗散量之间的耦合代价. 这些系数也解释了初值阈值为何同时依赖 $\nu$ 和 $\mu$ : 闭合非线性估计时会出现如 $\nu^{-1/2}\mu^{-5/6}$ 、 $\nu^{-1}\mu^{-1/2}$ 等损失, 必须由初值小性补偿, 不能只凭单个方程的线性耗散率推断阈值.

零模估计则直接针对水平平均方程. 作者为 $\omega_0$ 和 $\theta_0$ 的能量加入 $\|y\partial_y f_0\|^2$ 及带扩散参数的竖直导数项, 同时控制零模速度分量 $$u^1_0$$ . 关键处理是利用 Fourier 模态相互作用的结构, 对零模输运项做分部积分, 将其改写成含 $\partial_y u^1_0$ 、 $\partial_y\theta_0$ 或 $\partial_y\omega_0$ 的表达式, 再以非零模的流函数和耗散控制. 例如零模速度满足形如

\frac{d}{dt}\|u^1_0\|^2+\nu\|\partial_yu^1_0\|^2\lesssim E_{\omega,0}^{1/2}D_{\omega,\ne0},

的估计. 温度与涡量零模方程则以 $y\partial_y f_0$ 的加权能量控制空间矩项, 避免无界区域上缺少 Poincare 不等式所导致的直接衰减困难. 零模的这些加权控制并非附带正则性, 而是吸收非零模方程中带 $$y$$ 的非线性项所必需.

非线性闭合依靠将零模与非零模的能量不等式合并, 并用初值小性使所有高阶能量项可吸收. 作者定义总能量为零模能量、带指数权重的非零模能量以及时间积分耗散项之和. 非零模的耗散控制零模的非线性强迫, 而零模速度的 $L^\infty$ 控制又用于估计非零模输运项, 因此形成一个相互依赖的闭合系统, 而不是单向的线性估计. 关键不等式 (2.21)-(2.22) 将其约化为关于总涡量能量与总温度能量的代数控制, 最终由初值阈值中的不同参数幂次保证统一有界.

有限通道的差别在于可使用 Poincare 不等式控制零模的 $$L^2$$ 能量, 因而可在零模能量上加入 $e^{2\delta_0\kappa t}$ , 并获得 $e^{-c\kappa t}$ 的零模衰减. 相反, 对 $\mathbb R$ 与 $[0,\infty)$ , 证明只对非零模加增强耗散时间权重, 不声称零模指数衰减. 从提供的摘录看, 作者将上述能量构造作为主要新方法, 但未给出对某个单一乘子估计或抽象 hypocoercivity 定理的依赖; 因而更准确地说, 创新在于针对耦合 Poiseuille 系统设计含交叉项、空间矩和频率权重的能量, 并据此处理零模-非零模间的非线性反馈.

4. 与之前工作的比较

本文与既有工作的比较主要体现在模型、区域、衰减结论和分析方法四个方面. 首先, 在 Couette 流附近的 NSB 方程已有多种非线性稳定性结果. 文中回顾了 Deng 等人关于 $\nu=\mu$ 全耗散情形的结果, 其条件涉及 $\|\omega_{\rm in}\|_{H^b}\ll\nu^{1/2}$ 、 $\|\theta_{\rm in}\|_{H^b}\ll\nu$ 以及额外的温度水平分数阶正则性; 对仅竖直耗散也有相应结果. Zhang 和 Zi、Niu 和 Zhao 后续改进了 Couette 情形的初值条件. 这些工作说明剪切增强耗散与非线性稳定性已在常剪切背景下得到深入研究, 但本文的 Poiseuille 基流为 $$y^2$$ 型, 剪切率依赖 $$y$$ 且在中心位置退化, 并且论文讨论的焦点是同时存在的涡量-温度耦合. 因此其频率能量与空间矩控制不能直接由 Couette 结果替代.

其次, 与有限通道中的 Poiseuille NSB 结果相比, Wang 的工作在文中被引用为通过 resolvent estimate 证明非线性稳定性, 初始速度要求约为 $\min\{\mu,\nu\}^{2/3}$ 阶的 $H^{7/2+}$ 小量, 温度条件还包含 $$H^1$$ 与水平分数阶导数控制, 指数为 $$31/24$$ . 本文在有限通道中给出涡量和温度的各向异性初值条件, 并进一步证明非零模以 $e^{-c\min\{\nu,\mu\}^{1/2}t}$ 衰减, 零模则以 $e^{-c\min\{\nu,\mu\}t}$ 衰减. 因而相对于 Wang 文中报告的结论, 明确推进之处是建立增强耗散估计, 并将讨论扩展至整条带 $\mathbb T\times\mathbb R$ 和半条带 $\mathbb T\times[0,\infty)$ . 由于两篇工作的未知量范数和初值假设表达不同, 仅凭摘录不能断言本文的初值阈值在所有参数区间都严格优于 Wang 的阈值.

第三, 本文的增强耗散速率与不含温度耦合的 Poiseuille 二维 Navier-Stokes 结果相比较. 文中提到 Coti Zelati 等人在 $\mathbb T\times\mathbb R$ 上得到稳定性、增强耗散及转捩指数 $\alpha\leq3/4+$ , Del Zotto 将该指数推进到 $\alpha\leq2/3+$ ; Ding 等人在有限通道 Navier-slip 边界下得到阈值 $\alpha\leq2/3$ . 本文对 NSB 系统得到的非零模衰减率 $e^{-c\kappa^{1/2}t}$ 与这些 Navier-Stokes 工作的增强耗散率一致, 但由于耦合存在, 其初值条件同时约束涡量和温度, 并含有额外的各向异性正则性与 $$y$$ 加权矩. 论文备注 1.4 明确承认转捩阈值可能不是尖锐的; 它没有证明类似上述 $\alpha$ 指数意义下的最优性, 也没有给出超过所列小性尺度后的失稳构造. 因此本文的贡献应表述为耦合模型下的充分稳定尺度和增强耗散, 而不是最优阈值的确定.

第四, 方法上本文选择直接构造能量, 与文中提到的纯 Navier-Stokes Poiseuille 研究使用半群估计或 resolvent estimate 的路线有所区别. 作者认为涡量和温度的耦合使得直接使用相同的半群估计较困难, 因而以 hypocoercivity 型能量处理线性剪切耗散, 再由能量法处理非线性项. 这种做法的代价反映在较强的初值条件及对 $y\partial_y\omega_0$ , $y\partial_y\theta_0$ 和零模速度的控制上; 论文备注 2.5 特别指出, Poiseuille 流较 Couette 流更复杂, 非线性项更多, 零模速度估计不可缺少. 这也意味着方法虽直接, 但尚不能据此声称阈值最优或正则性要求最弱.

区域与边界条件方面, 本文覆盖无界竖直方向、半无界竖直方向和有限通道, 但三者结论并不完全相同. 无界及半无界情形没有零模指数衰减, 只有零模能量控制与非零模增强耗散; 有限通道依靠 Poincare 不等式获得零模 $O(\kappa)$ 衰减. 此外, 本文所处理的边界条件是使扰动涡量和温度满足齐次边界条件的设置. 作者备注 1.5 指出, 若采用无滑移速度边界条件, 边界层会要求比纯能量估计更复杂的处理, 并将其留待未来研究. 因此结果不能不加修改地推广到无滑移边界情形.

总的来说, 本文填补的是 Poiseuille 背景下耦合 NSB 系统在多种条带区域中的定量非线性稳定性与增强耗散分析空白, 尤其补上了全空间条带和半空间条带的研究, 并在有限通道中补充了增强耗散速率. 但从提供内容中无法确认其初值阈值是否与其他结果在统一范数下可直接比较, 也无法确认最优阈值、失稳侧结论或更弱正则性假设; 作者亦明确指出阈值尖锐性尚未解决.

7. 具有流出边界条件的正压流的弱 BV 稳定性与无黏极限

Weak-BV Stability and Inviscid Limits for Barotropic Flows with Outflow Boundary Conditions

作者: Jiayun Meng, Alexis F. Vasseur

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07574v1

摘要翻译: 守恒律系统（包括无黏可压缩流）的边值问题自然出现在物理模型中，例如运河中的浅水流；在这类模型中，流体可能经由边界流出区域。与柯西问题不同，无黏极限对其所依赖的黏性模型高度敏感：不同的黏性正则化可能为极限系统选取不同的有效边界条件。

Comment: 116 pages

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

1. 研究对象和背景

论文研究半直线 $\mathbb{R}_{+}=(0,\infty)$ 上的一维压强型可压缩流体。未知量为密度 $\rho$ 和速度 $$u$$ , 压强满足 $p(\rho)=\rho^{\gamma}$ , $\gamma>1$ . 无黏模型是 Euler 方程

\begin{cases} \rho_t+(\rho u)_x=0,\\ (\rho u)_t+(\rho u^2)_x+p(\rho)_x=0, \end{cases} \qquad x>0,\ t>0.

其物理能量熵及熵通量分别为

\eta(\rho,\rho u)=\frac{\rho u^2}{2}+\frac{\rho^\gamma}{\gamma-1}, \qquad q(\rho,\rho u)=(\eta(\rho,\rho u)+\rho^\gamma)u,

并要求满足熵不等式 $\partial_t\eta+\partial_xq\leq 0$ . 作者考察的黏性正则化是具有密度依赖黏性的 Navier-Stokes 方程

\begin{cases} \rho_t^\nu+(\rho^\nu u^\nu)_x=0,\\ (\rho^\nu u^\nu)_t+(\rho^\nu (u^\nu)^2)_x+p(\rho^\nu)_x =\nu\bigl(\mu(\rho^\nu)u_x^\nu\bigr)_x, \end{cases} \qquad x>0,

其中 $\nu>0$ 是黏性参数, $\mu$ 可退化或随密度增长。边界是流出边界, 黏性系统规定

$$ u u(0,t)=u_B<0. $$

此处原文边界条件为 $u^\nu(0,t)=u_B<0$ ; 即速度在边界取给定的负值。远场状态为 $U_*=(\rho_*,\rho_*u_*)$ , 且位于负速度的亚音速区域 $\Omega^-_{\mathrm{sub}}$ : 两个特征速度 $u-c(\rho)$ 和 $u+c(\rho)$ 异号, 其中 $c(\rho)=\sqrt{\gamma\rho^{\gamma-1}}$ . 因而 Euler 边界问题恰需一个边界条件。本文采用由 Navier-Stokes 边界层选出的有效条件 $$u(0,t)=u_B$$ , 而不是任意指定一个与黏性极限无关的边界关系。

问题的重要性在于, 全空间 Cauchy 问题中不同黏性近似通常选择同一熵解, 但存在边界时黏性边界层会影响极限边界条件; 换用人工黏性与物理 Navier-Stokes 黏性, 可能得到不同的有效边界条件。作者希望确定物理黏性模型是否确实收敛到具有相应边界条件的 Euler 小 $$BV$$ 解, 并证明该识别对边界诱发的复杂波相互作用仍稳定。初始 Euler 数据取于 $$E_0$$ , 要求 $$U^0$$ 局部有界、 $\inf\rho^0>0$ , 且相对于远场状态有可积扰动, 即

\int_{\mathbb R_+}\bigl(\rho^0|u^0-u_*|+|\rho^0-\rho_*|\bigr)\,dx<\infty.

既有理论留下的关键空白是: 物理黏性具有退化性, 通常不能得到对 $\nu$ 一致的总变差估计; 因而经典小 $$BV$$ 黏性极限论证不能直接使用。此前针对物理 Navier-Stokes 系统的有限能量极限结果也未识别小 $$BV$$ 极限必为唯一的小 $$BV$$ Euler 解。本文进一步把相关弱- $$BV$$ 稳定性框架从全空间推进到流出边界问题, 同时处理边界反射激波及一般波形。

2. 主要定理或主要结果

论文的主要结论是定理 2.3 与定理 2.5, 依赖黏性系数假设 2.1. 对大密度, 存在常数 $$C_1,C_2>0$$ 使 $C_1\rho^\alpha\leq\mu(\rho)\leq C_2\rho^\alpha$ ; 指数满足: 若 $1<\gamma\leq5/3$ , 则 $\alpha=\gamma-1$ ; 若 $\gamma>5/3$ , 则 $\gamma-1\leq\alpha\leq\gamma$ . 对 $0<\rho<1$ , 要求 $C_1\rho^{\widetilde\alpha}\leq\mu(\rho)\leq C_2\rho^{\widetilde\alpha}$ , 其中 $\widetilde\alpha=\alpha$ 若 $\alpha<1/2$ , 否则 $0<\widetilde\alpha<1/2$ . 文中指出低密度条件中较强的限制用于现有边界初边值 Navier-Stokes 解的存在理论。

定理 2.3. 给定负速度亚音速远场状态 $U_*=(\rho_*,\rho_*u_*)\in\Omega^-_{\mathrm{sub}}$ , 存在 $\varepsilon>0$ . 若 $U^0=(\rho^0,\rho^0u^0)\in E_0$ , $$u_B<0$$ , 且

\|U^0\|_{BV(\mathbb R_+)}+|u_*-u_B|\leq\varepsilon,

则对任一与 $$U^0$$ 相适配的 Navier-Stokes 初值族 $U_0^\nu$ , 相应 Navier-Stokes 解 $U^\nu$ 当 $\nu\to0$ 时在 $L^1_{\mathrm{loc}}(\mathbb R_+\times\mathbb R_+)$ 收敛到以 $$U^0$$ 为初值、满足流出边界条件 $$u(0,t)=u_B$$ 的 Euler 小 $$BV$$ 解。这里相适配初值是光滑的, 满足 $u_0^\nu(0)=u_B$ , 并且在 $\nu\to0$ 时满足文中给出的加权收敛条件

\int_{\mathbb R_+}|U_0^\nu-U^0|\,dx +\int_{\mathbb R_+}\rho_0^\nu|u_0^\nu-u^0|^2\,dx +\int_{\mathbb R_+}|\rho_0^\nu-\rho^0|^\gamma\,dx +\nu^2\int_{\mathbb R_+}|\partial_x\phi(\rho_0^\nu)|^2\,dx\longrightarrow0,

其中 $\phi'(\rho)=\mu(\rho)/\rho^{3/2}$ . 引理 2.2 另证明每个 $$E_0$$ 数据都能构造这样的相适配初值族。

定理 2.5. 在相同小 $$BV$$ 与边界小失配假设下, 设 $U_n^0\in E_0$ 收敛到 $$U^0$$ 的条件为

\int_{\mathbb R_+}|U_n^0-U^0|\,dx +\int_{\mathbb R_+}\rho_n^0|u_n^0-u^0|^2\,dx +\int_{\mathbb R_+}|\rho_n^0-\rho^0|^\gamma\,dx\longrightarrow0.

对每个 $$U_n^0$$ 取其 Navier-Stokes 无黏极限 $$U_n$$ , 则 $$U_n$$ 在时空局部 $$L^1$$ 中收敛到以 $$U^0$$ 为初值的 Euler 小 $$BV$$ 解。该结论给出无黏极限类中的定量初值稳定性框架, 但本摘录未显示定理所用稳定估计的显式模函数或收敛速率, 因而不能据此宣称有特定幂次误差阶。引理 2.4 还说明任意 $$E_0$$ 初值至少存在一个相应的无黏极限子列。定理适用于一般小 $$BV$$ 波形及边界相互作用, 不局限于光滑解或预先指定的一类 Riemann 解。

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

证明的核心障碍不是仅仅令黏性项形式上消失, 而是在缺少对 $\nu$ 一致的 $$BV$$ 控制、也没有已知的 Navier-Stokes 解一致 $L^\infty$ 界时, 仍须将黏性解与正确的 Euler 小 $$BV$$ 解稳定地比较; 此外, 近边界波与边界本身的相互作用不能通过全空间延拓消除。

首先, 作者在质量 Lagrange 坐标中组织相对熵论证。坐标变换为

y^\nu(x,t)=\int_0^x\rho^\nu(z,t)\,dz-\int_0^t\rho^\nu(0,s)u^\nu(0,s)\,ds,

并以比容 $v^\nu=1/\rho^\nu$ 和速度为未知量, 得到

\begin{cases} v_t^\nu-u_y^\nu=0,\\ u u u_t^\nu+(v^\nu)^{-\gamma}_y =\nu\left(\dfrac{\gamma\mu_L(v^\nu)}{v^\nu}u_y^\nu\right)_y, \end{cases}

原文第二式左端为 $u_t^\nu+p_L(v^\nu)_y$ , 其中 $p_L(v)=v^{-\gamma}$ ; 上式的黏性项即 $\nu(\gamma\mu_L(v^\nu)u_y^\nu/v^\nu)_y$ , $\mu_L(v)=\mu(1/v)/\gamma$ . 采用这组变量的目的, 是使黏性作用于具有非线性通量的变量, 并利用 Bresch-Désjardins 熵与熵耗散得到一致控制。文中进一步在 Eulerian 与 Lagrangian 表述中建立熵和耗散的统一估计。由坐标变换产生的 Lagrange 边界 $y_b^\nu(t)$ 满足

(y_b^\nu)'(t)=-\rho^\nu(0,t)u_B=-\frac{u_B}{v^\nu(y_b^\nu(t),t)}.

这使边界成为随黏性解运动的曲线, 而不是固定位置。

其次, 作者把 $$a$$ -收缩与移位的相对熵方法用于该边界问题。普通相对熵不一定对激波收缩: 微小初值扰动可改变激波速度并累积成明显的激波位置偏移。带权相对熵与依赖扰动的激波移位用于补偿这种非收缩性。本文目录显示证明中具体构造权函数、激波移位, 并分别估计稀疏波导数、非物理波、权函数导数与激波附近的主项、抛物项及边界项。提供的正文摘录没有这些命题的完整表达式和常数依赖, 因而无法可靠复述逐项不等式或精确参数选择。

再次, 为构造可与黏性解比较的参照解, 作者在 Euler 层面采用经修改的 front tracking 算法: 以有限个基本波近似小 $$BV$$ 解, 对基本波作平滑处理, 并将波及相互作用组织为可追踪的 elements 与 events。与仅证明 Euler 解存在的通常 front tracking 不同, 此处近似结构要同时兼容黏性波形与相对熵比较。小总变差使近似波保持亚音速传播, 但黏性解对应的 Lagrange 边界可能移动很快, 甚至可能以超音速方式越过近邻波, 导致近似边界条件反复失效。作者据引言所述设计特殊 resetting procedure, 在重置时恢复边界条件; 每次重置造成的增长由两次重置之间累积的边界贡献补偿。这样避免了依赖重置次数的统一上界, 这是把全空间稳定框架移植到边界情形的一项实质性设计。

最后, 证明把上述一致熵耗散、带权相对熵稳定性及近似解的全局控制结合起来, 先控制黏性解到适配近似解的差异, 再识别极限为小 $$BV$$ Euler 解; 对初值序列则利用相同稳定性传递收敛。作者还修改此前方案以去除低密度处黏性系数非退化的要求。具体 stopping times、events 以及估计的技术细节在提供的 PDF 摘录中没有完整出现, 因此关于这些环节的精细证明机制只能确认其章节安排与引言概述, 不能补写未展示的引理内容。

4. 与之前工作的比较

与全空间小 $$BV$$ 无黏极限理论相比, Bianchini 与 Bressan [2] 证明严格双曲系统的小 $$BV$$ 解可由非退化人工黏性近似得到, 但其方法依托小 $$BV$$ 结构, 不能直接覆盖物理 Navier-Stokes 系统的密度退化黏性, 因为后者通常没有关于 $\nu$ 的一致 $$BV$$ 界。本文处理的是半直线边界问题和物理黏性, 并强调黏性模型会选择有效边界条件; 因此其结论不是将 [2] 的全空间结果简单加上边界, 而是需要与 Navier-Stokes 边界层所选条件相容。

Chen 与 Perepelitsa [17] 对 $\mu=1$ 的一维 Euler 系统证明了有限能量数据下无黏极限的存在, 采用 compensated compactness. 该类结果能在较弱数据条件下获得某种极限, 但按本文引言的说明, 对小 $$BV$$ 初值并未判定极限是否仍为小 $$BV$$ 解, 也未解决其唯一性或与指定小 $$BV$$ 解的稳定识别。本文走的是不同方向: 假设初始扰动小 $$BV$$ , 换取极限的明确识别、边界条件确认及初值稳定性。

最直接的前身是 Chen、Kang 与 Vasseur [15] 的物理黏性弱- $$BV$$ 稳定性框架。本文沿用 $$a$$ -收缩与移位及适配黏性结构的近似思想, 但从全空间推进到流出边界。边界问题不能像不透流壁问题那样通过特殊对称延拓还原为全空间; 亚音速流入与流出也没有这种简化。已有边界 Euler 小 $$BV$$ 稳定性工作 [28] 通过控制 front tracking 近似的时间总变差及跟踪波元素来确保极限满足边界条件。本文面对的要求更强: 近似解还必须与 Navier-Stokes 黏性解可比较, 并处理移动 Lagrange 边界、波被边界追越以及重置导致的误差累积。

假设方面, 作者明确指出对 [15] 的算法作了实质修改, 以去掉其低密度处对黏性系数非退化的要求; 本文允许 $\mu(\rho)$ 在低密度区与 $\rho^{\widetilde\alpha}$ 同阶退化, 并称改进也适用于 [15] 的全空间情形。因此推进不仅是加入边界, 也扩大了可处理的黏性系数范围。不过, 定理仍要求正密度初值、规定的多项式黏性增长、远场负速度亚音速态以及足够小的总变差和 $$|u_*-u_B|$$ . 低密度指数限制在本文的边界初值存在理论中仍被需要; 文中指出这项更强限制只用于 Navier-Stokes 解的存在性, 在边界情形下若仅假设 $\mu$ 对所有 $\rho>0$ 与 $\rho^\alpha$ 可比, 相应存在唯一性理论尚未建立。

适用范围上, 论文强调并非只对平滑流或固定 Riemann 波形成立, 而覆盖小 $$BV$$ 解的一般波模式、波与边界的复杂相互作用以及反射激波。结果给出 $L^1_{\mathrm{loc}}$ 时空收敛及无黏极限类内的初值稳定性, 但摘录未提供收敛速率, 也不能据此推断对任意大总变差、真空初态、超音速边界态或一般非小扰动成立。作者未在所给内容中显式量化相对既有结果的最优性或具体常数代价; 这些比较无法进一步确认。总体而言, 本文的推进以小 $$BV$$ 条件和黏性系数结构为代价, 将物理 Navier-Stokes 到 Euler 的稳定极限识别扩展到具有有效流出边界条件的半直线问题, 并处理全空间理论无法覆盖的边界波相互作用。

8. 小外力作用下三维矩形喷管中常值状态稳态跨声速激波的稳定性

Stability of steady transonic shocks with constant states in 3D rectangular nozzles under small external forces

作者: Shangkun Weng, Wengang Yang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07586v1

摘要翻译: 本文研究矩形喷管中常值状态稳态跨声速激波在外力和出口压力受到小的三维扰动时的稳定性。从数学上看，该问题可表述为一个双曲—椭圆耦合偏微分方程组的非线性自由边界问题。一个关键难点在于背景激波的位置任意；在扰动作用下，这可能导致相同外力和出口压力对应多个激波解，从而揭示背景激波的不稳定性。我们用变形量和涡量重新表述稳态 Euler 方程，以有效分解双曲模态和椭圆模态。对重新表述后的方程组以及背景跨声速激波处的 Rankine–Hugoniot 条件进行线性化后，我们发现，非零水平外力这一类外力有助于唯一确定激波位置，其中质量守恒起着关键作用。我们发展了一种精细的迭代方法，以建立跨声速激波解的存在性和结构稳定性。我们的分析表明，平直喷管与扩张喷管中稳态跨声速激波的稳定机制在许多方面存在显著差异。

Comment: 无

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

1. 研究对象和背景

本文研究三维矩形喷管内完整可压缩 Euler 方程的定常跨声速激波，在小的保守外力和出口压强扰动下是否存在、是否稳定，以及扰动问题中的激波位置能否由给定数据确定。喷管为

N=[0,L]\times C,\qquad C=[-1,1]\times[-1,1],\qquad x'=(x_2,x_3),

入口为 $$x_1=0$$ ，出口为 $$x_1=L$$ ，四个侧壁固定并满足滑移边界条件。气体满足多方状态方程

p=S\rho^\gamma,\qquad e=\frac{p}{(\gamma-1)\rho},\qquad \gamma>1,

其中 $\rho$ 、 $$u=(u_1,u_2,u_3)$$ 、 $$p$$ 、 $$S$$ 分别是密度、速度、压强和熵变量。控制方程为带体力的定常质量、动量和总能量守恒律

\begin{aligned} &\nabla\cdot(\rho u)=0,\\ &\nabla\cdot(\rho u\otimes u)+\nabla p=\rho f,\\ &\nabla\cdot((E+p)u)=\rho f\cdot u, \end{aligned} \qquad E=\rho\left(\frac12|u|^2+e\right),

外力取保守形式 $f=-\sigma\nabla\Phi$ ，其中 $\sigma>0$ 小，且 $\Phi\in C^{3,\alpha}(N)$ 、入口上 $\Phi(0,x')=0$ 。势函数在侧壁满足相应的 Neumann 型兼容条件及文中给出的三阶兼容条件。入口流态固定为上游常状态，出口压强规定为 $p_+(L,x')=\bar p_++\sigma p_{\rm ex}(x')$ ，其中 $p_{\rm ex}\in C^{2,\alpha}(C)$ 并满足侧边兼容条件。

待求解的自由边界是激波面

S=\{x_1=\xi(x'),\ x'\in C\},

激波上两侧流场须满足质量、动量及广义 Bernoulli 量的 Rankine-Hugoniot 条件。背景解是无外力时的一维法向激波，它连接两常状态 $\bar U_-$ 与 $\bar U_+$ ，满足上游超声速、下游亚声速，且压强与熵均增加：

\bar M_->1>\bar M_+>0,\qquad [\bar p]>0,\qquad [\bar S]>0.

由于背景流在激波两侧均为常数，对任意 $L_s\in(0,L)$ ，把激波放在 $$x_1=L_s$$ 都得到同一类背景解。因此背景状态本身不包含激波位置的信息；这是本问题与许多已有稳定性问题的关键差异，也使仅以小扰动估计流场并不能自动确定自由边界。

这一位置退化具有实际的分析后果。对于给定入口状态和出口压强，若没有额外机制选定位置，扰动问题可能有多个不同激波解，因而不能直接把某个任意背景激波视作唯一稳定的参照。本文考察的选位机制来自小外力在横截面上的积分效应：外力改变流动，并通过质量守恒与下游边值问题的可解性条件，将激波位置与出口压强联系起来。作者要求每个内部截面上外力的轴向分量积分非零，即

\int_C\partial_{x_1}\Phi(t,x')\,dx'\ne0,\qquad 0<t<L.

同时，出口压强扰动的截面平均值须落在由外力势确定的位置函数的取值范围内。由此可先选出一个候选背景位置，再证明其附近存在唯一的小扰动解。

分析上，定常 Euler 系统在超声速区具有双曲型特征，在亚声速区则包含椭圆型约束；激波又使两区域通过未知边界上的跳跃条件耦合。作者使用形变-旋度分解，将速度方程与熵、广义 Bernoulli 量及第一旋度分量的输运方程结合，从而把不同类型的变量拆分处理。广义 Bernoulli 量为

B=\frac12|u|^2+\frac{\gamma}{\gamma-1}S\rho^{\gamma-1}+\sigma\Phi.

问题的核心不仅是获得常规的局部正则性，而是处理背景激波位置任意所导致的线性可解性退化，并在非线性迭代中持续选定位置。论文试图填补的是：在平直三维矩形喷管、完整三维扰动和常状态背景激波的框架下，利用小保守外力及出口压强数据确定激波位置，并建立相应自由边界解的结构稳定性。摘要和正文还强调，平直喷管中的位置选择机制不同于发散喷管中由几何膨胀直接进入线性化跳跃条件的机制。

2. 主要定理或主要结果

主要结论是定理 1.3。设背景常状态 $\bar U_\pm$ 满足法向激波关系及

\bar M_->1>\bar M_+>0,\qquad [\bar p]>0,\qquad [\bar S]>0.

设 $\Phi\in C^{3,\alpha}(N)$ 满足入口归一化、文中关于侧壁的一阶及三阶兼容条件，且 $p_{\rm ex}\in C^{2,\alpha}(C)$ 满足出口侧边兼容条件。定义

\bar K=\bar\rho_++\frac{\gamma-1}{\gamma+1}\bar\rho_->0,

以及

F(t)=\bar K\int_C\Phi(t,x')\,dx' -\frac{\bar\rho_+}{1-\bar M_+^2}\int_C\Phi(L,x')\,dx',\qquad 0\le t\le L.

假设出口扰动的截面平均满足严格范围条件

\inf_{0<t<L}F(t)<\int_Cp_{\rm ex}(x')\,dx'<\sup_{0<t<L}F(t),

并且轴向外力具有非退化的截面积分

\int_C\partial_{x_1}\Phi(t,x')\,dx'\ne0,\qquad 0<t<L.

则存在唯一 $\bar\xi_*\in(0,L)$ 满足位置选择方程

F(\bar\xi_*)=\int_Cp_{\rm ex}(x')\,dx'.

这里 $F'(t)=\bar K\int_C\partial_{x_1}\Phi(t,x')\,dx'$ ，第二个假设使 $$F'$$ 在连通区间 $$(0,L)$$ 上不为零，因而 $$F$$ 严格单调；严格范围条件保证所选位置在喷管内部。

进一步，存在 $\sigma_*>0$ ，使每个 $0<\sigma\le\sigma_*$ 对应一个唯一的、处于定理所述小扰动类中的跨声速激波解，且

U_\pm\in C^{2,\alpha}(N_\pm),\qquad \xi\in C^{3,\alpha}(C),

并有定量估计

\sum_{\pm}\|U_\pm-\bar U_\pm\|_{2,\alpha;N_\pm} +\|\xi-\bar\xi_*\|_{3,\alpha;C}\le C\sigma.

其中 $$C$$ 与 $\sigma_*$ 依赖于背景状态、喷管、 $\alpha$ 及给定数据。这里的 $N_\pm$ 是由实际激波面 $x_1=\xi(x')$ 划分出的上下游区域。所得流态保持上游超声速、下游亚声速：

u_{1-}^2-c_-^2>0\quad\text{于 }N_-,\qquad c_+^2-|u_+|^2>0\quad\text{于 }N_+,

并且激波熵增严格成立：

S_+(\xi(x'),x')-S_-(\xi(x'),x')>0,\qquad x'\in C.

位置选择不是预先固定激波面上一点，而由积分可解性条件确定。非线性迭代中，作者以 $\xi(1,1)=\bar\xi_*+\ddot\xi_*$ 参数化位置位移，并用下游形变-旋度问题的积分条件构造标量方程 $$I=0$$ 。在 $|\ddot\xi_*|\le C\sigma$ 时，其对位置位移的导数满足

\partial_{\ddot\xi_*}I =\frac{1-\bar M_+^2}{\bar\rho\bar u}\, \frac{1}{\sigma\bar K}\int_C\partial_{x_1}\Phi(\bar\xi_*,y')\,dy' +o(\sigma),

按文中公式 (1.44) 的记法，非零积分使主导系数非退化；同时 $I(0)=O(\sigma^2)$ 。对 $\sigma^{-1}I=0$ 使用隐函数定理，得到唯一 $\ddot\xi_*=O(\sigma)$ 。随后由切向 Rankine-Hugoniot 条件更新激波斜率，并恢复激波面。需要注意，以上导数公式是文中用于位置隐函数论证的展开；其精确系数按原文 (1.44) 书写时为

\partial_{\ddot\xi_*}I =\frac{1-\bar M_+^2}{\bar\rho\bar u}\,\frac{1}{\sigma\bar K} \int_C\partial_{x_1}\Phi(\bar\xi_*,y')\,dy'+o(\sigma).

作者还给出不满足非退化条件时可能出现多位置的明确例子。取仅依赖 $$x_1$$ 的势函数

\Phi(x_1)=\frac{L}{2\pi}\left(\cos\frac{2\pi x_1}{L}-1\right),

则 $f=\sigma(\sin(2\pi x_1/L),0,0)$ ，且 $F(t)=\frac{4\bar K L}{2\pi}(\cos(2\pi t/L)-1)$ 。对于文中选取的常出口扰动 $p_{\rm ex}=-\bar K L/(2\pi)$ ，位置方程有两个内部根 $$t_1=L/4$$ 和 $$t_2=3L/4$$ ；两根处的 $\int_C\partial_{x_1}\Phi$ 均非零，因而可分别在两处进行非线性构造。这说明定理的唯一性依赖其外力非退化假设，且外力方向改变时，同一入口状态、外力和出口压强可对应不同激波解。

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

证明的主要困难是背景激波位置自由导致线性化问题存在位置退化，而不是单纯地在已固定边界上建立椭圆估计。对常状态背景而言，线性化 Rankine-Hugoniot 条件没有发散喷管情形中那种直接含有激波位移的主导代数项。因此，不能仅从局部跳跃关系解出激波位置。作者的关键做法是把位置选择转化为下游速度边值问题的积分可解性条件，并追踪该条件在非线性迭代中的非退化性。

首先，使用形变-旋度分解改写定常 Euler 方程。以 $$U=(u,S,B)$$ 为变量， $$S$$ 与 $$B$$ 沿流线输运；第一旋度分量 $\omega_1$ 满足由 $\operatorname{div}\omega=0$ 导出的输运方程；横向旋度分量 $\omega_2,\omega_3$ 可由横向动量关系表达为 $u,S,B,\omega_1$ 的函数。速度则满足形变方程及旋度方程的耦合系统，其中形变方程为

\sum_{i,j=1}^3(c^2\delta_{ij}-u_i u_j)\partial_{x_i}u_j =\sigma u\cdot\nabla\Phi,

旋度部分为 $\operatorname{curl}u=\omega$ 。这种改写把熵、Bernoulli 量及部分涡量的信息放入输运方程，把速度的椭圆型部分留在形变-旋度系统中。下游亚声速区可据此建立速度边值问题；上游超声速区则利用入口数据沿流向构造解。三维情形中，这一分解比二维压力与流角的椭圆系统更适合处理三个速度分量及完整横向扰动。

其次，作者在线性层面先确定候选激波位置。将形变-旋度系统和激波条件在背景常状态激波处线性化，求解输运方程后形成上下游速度问题。通过对线性化修正散度方程在上下游区域积分，并结合质量 Rankine-Hugoniot 条件，得到位置可解性条件 (1.41)，即 $F(\bar\xi_*)=\int_Cp_{\rm ex}$ 。其本质是总质量通量平衡，而非人为固定激波面某点。外力通过修正散度方程中的源项进入该平衡，出口压力则通过下游边界数据进入。因此，截面平均的轴向外力不为零保证位置方程有非零导数，出口压力的严格范围条件则保证根位于喷管内部。

再次，非线性构造把位置确定步骤嵌入迭代，而非在初始线性化后就假定位置不变。先构造上游超声速流，再把下游亚声速区域通过坐标变换化到固定区域，以便在同一空间中建立迭代映射。给定迭代流场及激波切向斜率后，下游问题仍须满足一个整体积分相容条件；作者将其写成标量方程 $I(\cdots;\ddot\xi_*)=0$ 。利用导数展开 (1.44) 和 $I(0)=O(\sigma^2)$ ，通过隐函数定理逐步更新唯一的 $O(\sigma)$ 激波位移。激波的两个切向导数再由切向 Rankine-Hugoniot 条件确定；其兼容关系提供第一旋度分量的边界数据。这个安排同时处理位置、激波斜率和涡量边界数据之间的耦合。

下游形变-旋度线性问题还有一个细致的相容性障碍：迭代产生的旋度源项一般并非无散场，不能直接作为旋度方程的右端数据。作者依照 Weng-Xin [34] 的思路引入辅助标量函数，并为它设置齐次混合边界条件；用其梯度修正旋度源项，使修正后的数据满足无散条件。证明该辅助函数在迭代不动点处恒为零后，扩大的系统才还原为原有形变-旋度系统。这一步不是形式上的变量增加，而是用于补足旋度方程的相容性，使固定域线性问题可解。

最后，估计须覆盖激波与平壁相交的棱边以及侧壁彼此相交的角点。势函数和出口数据的兼容条件允许关于侧壁作偶、奇延拓；输运方程和激波关系须保持这种延拓结构。作者据此获得流场与激波斜率的一致 $C^{2,\alpha}$ 控制，并在 $C^{1,\alpha}$ 范数下建立压缩映射，最终得到 $\xi\in C^{3,\alpha}$ 。因此，证明中的链条是：分解方程以隔离双曲输运与椭圆速度部分，线性积分条件选位，辅助势函数恢复旋度数据相容性，固定域迭代控制非线性自由边界，边角延拓完成全域正则性。文中给出的材料未展开各线性边值问题的全部具体估计常数及逐项证明细节，不能据此补写未展示的估计公式。

4. 与之前工作的比较

与平直喷管中常状态背景激波的既有研究相比，本文的主要推进是直接处理激波位置不预先指定的问题，并使用外力与出口压强共同选位。文中回顾 Chen [8]、Chen-Yuan [9] 对三维管道问题的结构稳定性结果时指出，其设置要求激波前沿通过一个预先指定点，并将出口压强规定到一个加法常数的意义下。本文则不给激波面指定通过点，而由积分可解性条件确定候选位置；但相应地，作者对外力积分及出口压强平均值施加了明确的非退化与范围假设。这是边界设定和位置机制上的差别，不宜表述为无条件放宽所有数据限制。

最接近的位置选择先例是 Fang-Xin [16] 的二维完整 Euler 近乎平直有限喷管问题。该工作通过线性化自由边界问题的可解性条件筛选初始激波位置，再在可接受位置附近构造非线性解，并指出同一接收压强可能对应多个位置，其多重性依赖喷管几何。本文沿用“先解线性位置条件，再作非线性构造”的思路，但对象从二维推广到三维矩形喷管，并将位置选择因素从喷管几何转为保守外力。Fang-Xin 使用压力和流角的一阶椭圆系统与熵、Bernoulli 量输运方程相耦合；作者指出该椭圆分解依赖二维结构，不能直接扩展到一般三维流动。本文改用形变-旋度分解，建立适用于三个速度分量的三维框架。三维推广并非只增加一个横向变量，而是需要处理旋度数据相容性、固定域下游问题和边角正则性。

与 Fang-Gao [14] 的三维轴对称圆柱喷管位置研究、Fang-Gao-Xiang-Zhao [15] 的大旋流位置选择相比，本文考察非轴对称三维扰动下的矩形喷管及小保守外力机制。文中提到 [15] 中旋流对选定激波位置起关键作用；本文则由截面轴向外力积分通过质量守恒进入位置方程。与 Weng-Yang [37] 的二维等熵外力问题、Weng-Zhang-Zhou [39] 的三维矩形圆柱中的相应结果相比，本文关注完整 Euler 系统中的三维矩形喷管和常状态背景激波。作者特别区分：前两项相关工作在背景问题中已有非零轴向外力，且背景流态非恒定；本文从无外力的常状态法向激波出发，再施加小外力和出口压力扰动。两类背景的选位结构因而不同。

与发散喷管中的结果相比，本文的稳定机制也不相同。相关工作包括 Li-Xin-Yin [20,19,22] 对二维发散喷管、Xin-Yin [41] 对完整 Euler 方程相关问题、Li-Xin-Yin [21] 对圆锥喷管、Weng-Xin [34] 对三维圆柱喷管，以及 Liu-Xu-Yuan [23] 对球形激波的研究。作者的 Remark 1.6 指出，在 [34] 的发散几何中，膨胀项进入线性化 Rankine-Hugoniot 条件并产生激波位移的主导代数项；熵输运把位置信息传递到下游形变-旋度系统，后续归结为带非局部项的椭圆问题，几何符号还提供斜导数边界条件所需的可解性符号。本文的平直喷管没有这一几何主项：位置必须从修正散度方程的积分相容条件中恢复，外力通过质量守恒提供非退化性。因此，不能把发散喷管中由几何提供的稳定机制直接套用到常状态平直背景。

本文也明确说明假设失效时的代价。若去掉 $\int_C\partial_{x_1}\Phi\ne0$ ，位置函数可能不单调；文中的振荡外力例子给出两个可行初始位置，并可分别构造两个非线性解。因此定理的唯一性是带条件的局部结构稳定性，而非对任意小外力、任意出口扰动的全局唯一性。定理还要求出口平均压强扰动严格落在 $$F$$ 的内部值域中，保证位置在 $$(0,L)$$ ；端点情形及退化根没有由所给定理覆盖。证明所得 $O(\sigma)$ 估计是在背景状态与数据固定、扰动幅度足够小时的 Hölder 范数控制，并未声称大扰动稳定、时间演化稳定或激波位置对一般数据的全局连续性质。

从所提供的 PDF 摘录可以确认上述比较以及文中点名的前人工作；摘录没有给出这些先前论文的完整定理条件、常数和证明细节，因此不能据此逐项比较其精确函数空间、最优正则性或所有假设是否严格强弱。作者明确强调的结论限于：三维矩形喷管中利用小外力选定常状态背景激波位置，并在所列非退化条件下建立小扰动解的存在、局部唯一性及结构稳定性。

9. 非线性反问题中依赖误差下的后验一致性及其在时间演化方程中的应用

Posterior Consistency under dependent errors in non-linear inverse problems with applications to time evolution equations

作者: Aditya Jambhale

类别: math.ST, math.AP

发布日期: 2026-10-05

arXiv ID: 2610.07240v1

摘要翻译: 设 $θ\to \mathcal G(θ)$ 是一个映射，它将参数映射到某个刻画动力系统的非线性偏微分方程的解。我们给出了一些条件，在观测 $$4b$$ 中含有相关噪声时，这些条件能够保证关于 $$θ$$ 的后验一致性；这里的相关噪声用于刻画测量所用传感器可能存在的时间相关性。我们给出两个例子：具有未知黏性系数和初始条件的 Navier-Stokes 方程，以及扩散系数在空间中变化（非均匀）的反应-扩散型方程。后一个例子需要针对具有时间累积效应的 Darcy 输运算子建立新的稳定性估计。

Comment: 无

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

1. 研究对象和背景

本文研究由非线性偏微分方程解所定义的贝叶斯逆问题. 设参数为 $\theta\in\Theta$ , 前向映射 $\mathcal G$ 将参数送到时间区间 $$[0,1]$$ 与空间区域 $\Omega$ 上的解, 即

\mathcal G: \Theta\longrightarrow L^2(Q,V),\qquad Q=[0,1]\times\Omega.

目标是在观测解而非直接观测参数的情形下, 由数据恢复真实参数 $\theta_0$ . 论文采用空间随机、时间确定的混合设计: 在 $$T$$ 个时间点 $$t_i=(2i-1)/(2T)$$ 上, 每个时间点有 $$n$$ 个空间观测, 总样本量 $$N=nT$$ , 且

Y_{ij}=\mathcal G(\theta_0)(t_i,X_{ij})+\epsilon_{ij},\qquad X_{ij}\stackrel{\mathrm{i.i.d.}}{\sim}\operatorname{Unif}(\Omega).

同一传感器在不同时间的误差可以相关, 不同传感器的误差相互独立; 将误差整体记作 $\epsilon\sim N(0,\Gamma_N)$ , 并假设协方差矩阵可逆. 这种设定意在描述重复使用同一传感器所产生的时间相关误差, 而非通常的独立同分布噪声. 论文关注后验是否随 $$N$$ 增大而集中到真值附近, 以及在何种协方差、先验和前向映射条件下可以得到定量收缩结论.

具体应用有两个时间演化问题. 第一个是环面 $\Omega=\mathbb T^d$ 上的异质反应扩散方程. 已知反应项 $$f$$ 和初值 $\varphi$ , 未知空间扩散率 $$D(x)$$ , 方程为

\partial_t u-\nabla\cdot(D\nabla u)=f(u)\quad\text{于 }(0,1]\times\Omega,\qquad u(0,\cdot)=\varphi\quad\text{于 }\Omega,

其中 $D(x)\ge D_{\min}>0$ , $$f$$ 为光滑且紧支撑的时间无关函数. 解记为 $$u_D$$ , 前向映射为 $\mathcal G(D)=u_D$ . 与常系数反应扩散模型相比, 这里的 $$D$$ 随空间变化, 反演难点在于从一个时空解中识别空间变系数. 论文强调, 时间演化可使解的梯度信息在时间上累积, 因而在适当非退化条件下, 可以建立比单一时刻观测更有力的反向稳定性. 其关键条件之一是存在 $\delta>0$ 使

\inf_{x\in\Omega}\int_0^1|\nabla u_{D_1}(t,x)|^2\,dt\ge\delta.

直观上, 该条件排除某些空间位置在整个观测时间内都缺乏梯度信息的情形; 文中还声称这类初值在适当拓扑中是泛性的, 并在实验者可选择初值时给出依赖已知量的选择方法. 但所给摘录没有包含相应命题的证明及完整函数空间假设, 因而这些细节不能在此进一步核实.

第二个应用是二维环面上的周期 Navier-Stokes 方程, 速度场 $$u$$ 无散、均值为零, 外力 $$f$$ 已知, 黏度 $\nu$ 与初值 $\varphi$ 未知:

\partial_tu+\nu Au+B(u,u)=f,\qquad \nabla\cdot u=0,\qquad u(0,\cdot)=\varphi,

其中 $A=-P\Delta$ , $B(u,u)=P[(u\cdot\nabla)u]$ , $$P$$ 是投影到无散、周期、零均值向量场空间的算子, 且 $\nu\ge\nu_{\min}>0$ . 此处前向映射为 $\mathcal G(\nu,\varphi)=u_{\nu,\varphi}$ . 反演同时涉及有限维黏度与函数值初值, 所需稳定性须将两者差异由时空解差异控制.

上述问题的重要性在于, 实际测量往往按时间重复采样, 相关噪声会改变似然几何与可用检验方法, 不能不加论证地套用独立误差下的 Hellinger 检验理论. 本文试图补足两方面空白: 一是为具有时间依赖高斯误差的非线性回归逆问题建立一般后验收缩框架; 二是针对两个抛物型方程验证前向映射的正则性及反向稳定性, 尤其为异质反应扩散问题发展带时间累积机制的 Darcy 型输运算子估计. 摘录提及相关贝叶斯 PDE 研究及 i.i.d. 误差文献, 但没有提供所有应用定理的完整陈述, 因此各 PDE 例子的精确参数范围与收缩率需以全文后续部分为准.

2. 主要定理或主要结果

一般统计结果以定理 1.12 和定理 1.14 为核心. 观测误差为协方差 $\Gamma_N$ 的联合高斯向量, 首先要求最大方差方向有界:

\|\Gamma_N\|_{2\to2}\le C_1.

此外, 还需满足以下两种条件之一: 精度矩阵最大绝对行和一致有界,

\max_i\sum_j| (\Gamma_N^{-1})_{ij}|\le C_2,

或前向解具有参数无关的空间平均 $\int_\Omega\mathcal G(\theta)(t,x)\,dx=A(t)$ , 且 $\operatorname{tr}(\Gamma_N^{-2})\le C_2N$ . 命题 1.6 说明, 若每个传感器的时间误差是平稳高斯 $\operatorname{VAR}(p_j)$ 过程, 阶数及自回归矩阵逐项一致有界, 创新协方差的最小特征值一致远离零, 则精度矩阵行和条件成立. 命题 1.7 在这些协方差条件下, 若两个候选前向解一致有界, 则其数据分布满足

\operatorname{KL}(p_\theta^N,p_\vartheta^N)+V^{2,0}(p_\theta^N,p_\vartheta^N) \le C N\|\mathcal G(\theta)-\mathcal G(\vartheta)\|_{L^\infty_tL^2_x}^2.

先验取中心高斯先验 $\Pi'$ 的缩放版本. 若参数定义域维数为 $\ell$ , 高斯先验再生核希尔伯特空间嵌入适当的 $H^\alpha$ 空间, 则采用

\theta=N^{-\ell/(4(\alpha+\kappa)+2\ell)}\theta',\qquad \theta'\sim\Pi',

其中 $\kappa$ 控制前向映射对参数弱范数的稳定性. 定理 1.12 假设前向映射满足条件 1.9 的正则性以及条件 1.10 的抛物型 Sobolev 正则性, 其中 $a\ge(d+1)\vee6$ , 并令

\delta_N=N^{-a/(2a+d+2)},\qquad \varepsilon_N=N^{-(\alpha+\kappa)/(2\alpha+2\kappa+\ell)},\qquad \varepsilon_N\lesssim\delta_N.

若设计分配满足

N^{(2d+3)/(2a+d+2)}\lesssim n\lesssim N\delta_N, \qquad \delta_N^{-1}\lesssim T\lesssim N^{(2a-d-1)/(2a+d+2)},

且真参数属于相应正则化集合 $\Theta_N$ , 则可构造检验 $\phi_N$ , 使第一类错误趋于零, 且对任意固定 $K\in\mathbb N$ 和足够大的常数 $$C$$ ,

\sup_{\substack{\theta\in\Theta_N:\,\|\mathcal G(\theta)-\mathcal G(\theta_0)\|_{L^2(Q)}>C\delta_N}} \mathbb E_\theta^N[1-\phi_N]\le e^{-KN\varepsilon_N^2}.

定理 1.14 在上述条件上再假设 $\theta_0\in R\cap H$ , 条件 1.3 及条件 1.4 或 1.5 成立, 得到后验以 $\delta_N$ 的尺度在前向解度量中收缩. 按定理原文的概率表述, 对每个 $$b>0$$ , 取充分大的 $$m$$ 后有

P_{\theta_0}^N\!\left(\Pi_N\!\left(\theta\in\Theta,\ \|\mathcal G(\theta)-\mathcal G(\theta_0)\|_{L^2(Q)}\le m\delta_N\mid D_N\right)\le1-e^{-bN\varepsilon_N^2}\right)\longrightarrow0.

若进一步满足条件 1.13 的反向稳定性, 即在有界正则化集合上

\|\theta-\theta_0\|_{L^2}\lesssim\|\mathcal G(\theta)-\mathcal G(\theta_0)\|_{L^2(Q)}^\eta,

则参数后验在 $$L^2$$ 距离下以 $\delta_N^\eta$ 的尺度收缩, 概率界同样由 $e^{-bN\varepsilon_N^2}$ 控制. 这里 $\varepsilon_N$ 是先验小球与熵论证中的尺度, $\delta_N$ 是解空间检验的分辨尺度; 不应将两者混同. 摘录明确指出依赖误差下该直接论证所得速率劣于 i.i.d. 情形, 并说明若能另行构造足够快的估计量, 可通过线性化前向映射及惩罚最小二乘构造 plug-in 检验, 将后验推进到 i.i.d. 速率; 该更优速率不是定理 1.14 直接证明的结论.

PDE 应用方面, 对异质反应扩散方程, 摘录给出的引理 3.8 结论形式为: 在扩散率及解满足适当光滑性、椭圆性和上述时间累积非退化条件时, 存在由平滑度决定的指数 $\eta$ , 使

\|D_1-D_2\|_{L^2(\Omega)}\lesssim \|\nabla\cdot((D_1-D_2)\nabla u_{D_1})\|_{L^2(Q)}^\eta.

这为系数反演提供反向稳定性基础. 对二维 Navier-Stokes 方程, 命题 3.13 建立的目标是

|\nu_1-\nu_2|+\|\varphi_1-\varphi_2\|_{L^2(\Omega)} \lesssim\|u_{\nu_1,\varphi_1}-u_{\nu_2,\varphi_2}\|_{L^2(Q)}^\eta,

其中 $\eta$ 依赖初值光滑度. 但所提供 PDF 摘录止于定理 1.14, 未包含引理 3.8、命题 3.13 的完整假设与指数公式, 因而无法确认应用定理中的精确参数范围和收缩率. 另有一项计算结果: 在二维 Navier-Stokes 初值估计问题中, 若 $\varphi_0\in H^\alpha(\Omega)$ , $\alpha>d/2+1$ , 且黏度和初值 $H^\alpha$ 范数有界, 则存在计算量为 $O(N^\kappa)$ 、 $\kappa<2$ 的估计量 $\hat\varphi_N$ , 满足

P_{\nu_0,\varphi_0}^N\!\left(\|\hat\varphi_N-\varphi_0\|_\infty\le C_1N^{-\beta/2}\right) \gtrsim1-\exp\!\left(-cN^{(d+2)/(2(\alpha+1)+d+2)}\right), \quad \beta=\frac{2(\alpha+1)-d-2}{2(\alpha+1)+d+2}.

作者说明该初始化结果在 i.i.d. 误差情形也成立, 并将其用于讨论多项式时间贝叶斯计算的条件.

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

证明的统计核心困难是相关噪声下观测似然不再按独立样本分解, 因而经典独立误差下基于 Hellinger 距离的检验论证不能直接照搬. 作者没有仅把相关性隐藏在一个方差常数中, 而是明确控制协方差和精度矩阵, 再将相关高斯似然的复杂度转化为前向均值差的函数范数控制. 条件 1.3 控制最大方差方向, 用于高斯集中; 条件 1.4 对精度矩阵施加最大绝对行和界, 或条件 1.5 利用解的已知空间均值结构并控制 $\operatorname{tr}(\Gamma_N^{-2})$ . 由此得到命题 1.7 的 KL 散度和对数似然方差界, 使小球概率及后验分母估计仍可在依赖数据下实施. 命题 1.6 将抽象矩阵条件落实到有限阶平稳 VAR 模型, 说明时间相关并不要求误差独立, 但仍需相关结构具有一致可控性.

第二个关键环节是构造检验. 作者采用基于小波的函数估计器, 利用前向解的光滑度和集中不等式, 区分前向解在 $$L^2(Q)$$ 中相距至少 $C\delta_N$ 的参数. 这一做法沿用 Giné 和 Nickl 的小波检验思路, 但需要处理两项结构问题: 其一, 数据在时间上是确定网格而空间位置随机, 检验直接看到的是离散时间经验范数, 而目标度量是连续时空 $$L^2(Q)$$ 范数; 其二, 同一传感器的噪声在不同时间相关, 需要将高斯向量协方差带入误差概率控制. 作者通过命题 2.14 控制时间设计离散化误差, 并在上述协方差条件下得到指数型第二类错误界. 定理 1.12 所要求的 $$n,T$$ 区间正是保证空间随机采样误差、时间网格近似误差与小波分辨尺度相容的条件; 不能任意选择时间点数和每时刻传感器数.

先验部分采取缩放高斯先验与正则化集合相结合的框架. RKHS 嵌入 $H^\alpha$ 空间用于量化先验集中和逼近能力, 而正则化空间 $$R$$ 同时约束真值与候选参数的有界性. 条件 1.9 要求在 $$R$$ 有界集上解一致有界, 且前向映射对参数的 $(H^\kappa)^*$ 距离局部 Lipschitz; 条件 1.10 规定先验方向及正则化集合上的解具有抛物型 Sobolev 正则性. 这里使用各向异性空间 $H^{a,a/2}(Q)$ , 是为了反映抛物方程空间与时间导数的不同尺度, 而不是将时空正则性错误地视为各向同性. 这些条件共同决定检验尺度 $\delta_N$ 、先验尺度 $\varepsilon_N$ 以及样本分配约束. 后验收缩再通过后验收缩的一般定理完成, 并借助条件 1.13 将解空间收缩传递为参数空间收缩. 该最后一步本质上依赖逆映射的局部 Hölder 稳定性, 没有该条件时只能直接得到前向解度量下的结论.

PDE 证明的主要分析障碍是反演稳定性, 而不是仅证明正问题解存在. 对异质反应扩散方程, 将两组扩散率之差记为 $$D_1-D_2$$ , 两个方程相减后会出现由 $\nabla\cdot((D_1-D_2)\nabla u_{D_1})$ 生成的 Darcy 型输运项. 若固定时刻的梯度存在临界点, 直接除以梯度或点态反演便会退化. 作者利用整个时间区间上的梯度能量累积条件

\inf_x\int_0^1|\nabla u_{D_1}(t,x)|^2\,dt\ge\delta

避免要求每个时刻处处非退化, 并据此建立扩散率差的 Hölder 型控制. 论文摘要称此处需要关于带时间累积的 Darcy 输运算子的新稳定性估计. 此外, 附录 B 目录显示作者研究了散度形式算子的抛物方程稳定性, 附录 E 汇总 PDE 正则性. 但摘录未给出这些估计的具体插值步骤、指数推导或引理细节, 因而不能进一步断言其精确技术形式.

对 Navier-Stokes 方程, 证明目标是从全时空速度差中同时恢复 $\nu$ 和 $\varphi$ . 非线性项 $$B(u,u)$$ 、黏性系数与初值的耦合意味着不能简单把初值估计与黏度估计独立处理. 摘录表明作者通过命题 3.13 建立联合反向稳定性, 并通过抛物正则性验证一般收缩定理的前向条件. 另外, 对计算初始化器, 作者把椭圆问题中的 Wang (2026) 方法改写到抛物问题, 采用更适合抛物方程的惩罚, 得到初值估计的 $L^\infty$ 误差概率界和低于二次的浮点运算复杂度. 摘录没有给出离散化矩阵、算法迭代细节或非线性求解步骤, 因此这些不宜进一步臆测.

4. 与之前工作的比较

在统计模型层面, 作者将贝叶斯非线性逆问题从常见的独立同分布观测误差推进到按传感器分组、允许跨时间相关的高斯误差. 引言把 Stuart (2010) 列为贝叶斯 PDE 反问题的奠基工作, 并提及 Nickl (2020, 2023)、Monard 等 (2021a, 2021b) 等关于理论基础的进展; 时间演化问题相关文献包括 Nickl and Titi (2024)、Konen and Nickl (2025, 2026)、Nickl 等 (2025)、Castre and Nickl (2026)、Nickl (2024)、Kekkonen (2022) 及 Magra 等 (2026). 文中明确指出, 依赖观测误差的既有研究相对有限, 并点名 Nickl and Seizilles (2026) 与 Seizilles and Siebel (2026). 但摘录没有逐篇说明这些工作的模型、误差结构与定理差别, 因而不能据此作更细的文献优先权判断.

与 i.i.d. 情形相比, 本文的主要代价是额外的协方差结构假设和更复杂的样本设计条件. 所需时间相关性不能任意增长: 一方面要求 $\|\Gamma_N\|$ 一致有界, 另一方面还需精度矩阵行和受控, 或利用解的已知空间平均并满足迹条件. 这些条件覆盖有限阶、参数一致受控的平稳 VAR 误差, 文中还称可覆盖某些长程依赖时间序列, 但摘录没有展示后一类模型的完整条件. 当解的空间均值已知时, 条件 1.5 可替代较强的精度矩阵行和条件; Navier-Stokes 的零均值结构是作者给出的例子. 这表明本文结果利用了模型结构缓解相关噪声, 并非对所有可能协方差矩阵都成立.

在后验速率上, 作者明确承认其直接相关误差证明得到的速率劣于 i.i.d. 情形. 定理 1.14 的速率由 $\delta_N=N^{-a/(2a+d+2)}$ 控制, 并受到设计约束、前向正则性 $$a$$ 、参数先验平滑度 $\alpha$ 与弱稳定指数 $\kappa$ 的共同限制. 作者进一步指出, 若先有一个足够快的估计量, 可在线性化前向映射后构造凸的惩罚最小二乘问题, 再以 plug-in 检验恢复 i.i.d. 情形下的速率; 这属于额外策略的说明, 并非当前直接检验定理已实现的结果. 因此本文贡献更集中于证明依赖误差下后验仍可一致收缩, 而非在所有情形下达到最优速率.

在设计结构上, 时间点固定而空间位置随机, 是为了让误差相关矩阵非随机, 同时避免完全确定空间网格带来的额外近似误差. 作者认为其他设计也可考虑, 前提是相应经验范数能适当集中到 $$L^2(Q)$$ 范数. 该选择便于理论处理, 但实际采样若时间点随机或传感器之间也相关, 需重新验证经验范数近似与协方差界, 不能直接援引定理 1.14.

在 PDE 方面, 对异质反应扩散方程的推进在于将时间积累纳入 Darcy 型输运稳定性: 其非退化假设作用于 $\int_0^1|\nabla u(t,x)|^2dt$ , 而非要求单个时刻的梯度处处不消失. 这利用演化过程的信息积累, 是相对于只看静态或单时刻反演的重要机制. 作者还报告初值可选时可构造满足条件的初值, 且适当拓扑中这类初值泛性成立. 不过, 与此前异质扩散系数反演结果相比具体放宽了哪些系数正则性或可辨识性假设, 从提供内容中无法确认; 相关文献 Bonito 等 (2017)、Alberti 等 (2025)、Nickl (2024)、Nickl 等 (2020) 和 Nickl and Seizilles (2026) 虽被提及, 摘录没有展开逐项比较.

对二维 Navier-Stokes, 本文把未知黏度与未知初值作为联合参数, 目标是以全时空速度观测控制二者的差异. 这与仅恢复初值或仅估计方程系数的问题不同, 但命题 3.13 的精确假设、指数及相对于 Konen and Nickl (2025) 等前作的具体增量均未出现在所给页面, 作者也未在摘录中显式展开这些比较. 因而只能确认其在本文框架中提供了所需的 Hölder 型反向稳定性, 不能据此声称其对已知文献的某项假设作了确定放宽.

最后, 论文还将理论与计算联系起来: 在 Navier-Stokes 设定下构造多项式浮点复杂度的初值初始化器, 指数 $\kappa<2$ , 并说明该结果在 i.i.d. 误差情形也成立. 作者将其与 Wang (2026) 的椭圆情形方法联系, 称抛物惩罚选择可带来更合适的误差率, 并将初始化器视为 Nickl and Wang (2024) 多项式时间贝叶斯计算条件的验证工具; 同时提及缺少此类初始化器时 MCMC 可能需要指数运行时间, 参照 Bandeira 等 (2023). 该计算结论使用的成本定义只计矩阵-向量乘法和线性系统求解涉及的浮点操作, 与 MCMC 文献中的成本概念不同. 摘录没有提供完整算法实现及其数值条件数分析, 因此不应把该复杂度陈述泛化为实际墙钟时间保证.

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

1. 具有物理碰撞频率的 BGK 方程的大数据全局解

Global large-data solutions for the BGK equation with physical collision frequency in a slab

作者: Young-Pil Choi, Sihyun Song

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08658v1

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

摘要翻译

中文标题: 具有物理碰撞频率的 BGK 方程的大数据全局解

摘要翻译: 我们为薄层区域中的 Bhatnagar–Gross–Krook 方程建立了大数据全局存在理论，其中物理碰撞频率由局部密度给出。我们考虑入流边界条件和 Maxwell 边界条件，并且对初始数据仅假设具有有限质量、动能和熵，以及相应的入流通量界。数据无需小，也无需点态有界；我们既不要求密度或温度具有正下界，也不要求给定二阶以上的速度矩。据我们所知，在这些假设下，这是针对该初边值问题的首个非微扰全局存在性结果。核心障碍在于能量紧性与未经重整化的碰撞算子的可积性同时丧失。有限能量本身不能保证二阶速度矩的一致可积性，而依赖于密度的碰撞频率又阻碍了对增益项和损失项分别获得局部 $$L^1$$ 控制。我们通过关于碰撞频率截断一致的补偿估计、由数据确定的 de la Vallée Poussin 权函数，以及利用熵产生进行的两次连续重整化，克服了这些困难。由此得到宏观能量、局部 Maxwellian 以及非线性碰撞项的紧性，进而得到全局重整化解。

2. 关于带扰动与不带扰动的 Hardy 算子的波算子

On wave operators for Hardy operators with and without perturbations

作者: Atsuhide Ishida, Konstantin Merz

类别: math.AP, math-ph

发布日期: 2026-10-06

arXiv ID: 2610.08594v1

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

摘要翻译

中文标题: 关于带扰动与不带扰动的 Hardy 算子的波算子

摘要翻译: 我们证明了时间依赖散射理论中波算子的渐近完备性；这些波算子交织通常拉普拉斯算子或分数阶拉普拉斯算子与相应的 Hardy 算子，即加上次临界或临界 Hardy 势的分数阶拉普拉斯算子，条件是 Hardy 势为短程势。此外，我们证明了交织带有与不带外部短程扰动的 Hardy 算子的波算子的渐近完备性。我们还证明，当 Hardy 势或外部势为长程势时，所考虑的波算子不存在。我们通过三种不同的方法得到这些结果。第一种方法采用抽象半群框架，研究带短程势的 Fourier 乘子所对应波算子的渐近完备性。第二种方法采用 Kato 光滑性理论。第三种方法中，我们为带扰动与不带扰动的 Hardy 算子建立并应用新的传播估计。

3. 利用 Harnack 型方法获得退化非局部方程的瞬时正则化

Instantaneous regularisation via Harnack-type methods for a degenerate non-local equation

作者: Simon M. Schulz

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08434v1

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

摘要翻译

中文标题: 利用 Harnack 型方法获得退化非局部方程的瞬时正则化

摘要翻译: 我们证明了一种退化非局部演化方程的瞬时平滑效应（该方程最初由活性布朗粒子系统的多粒子极限导出）。其内在机制是：通过为某个平均量导出 Harnack 型不等式，我们能够证明该方程在正时间时变为强抛物型，因而具有自正则化性质。本文显著改进了此前的一项工作；在该工作中，对初始数据施加了额外的（弱）非退化性假设。本文完全去除了这些假设，并且除正则性迭代提升外，还通过一种 Bombieri-Giusti 式的方法，对上述平均量证明了定量估计。

4. $$L^p$$ 估计：Dunkl–薛定谔算子

$$L^p$$ estimates for Dunkl–Schrödinger operators

作者: Agnieszka Hejna-Łyżwa, Suman Mukherjee

类别: math.AP, math.FA

发布日期: 2026-10-06

arXiv ID: 2610.08371v1

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

摘要翻译

中文标题: $$L^p$$ 估计：Dunkl–薛定谔算子

摘要翻译: 我们研究 $$L=-Δ_k+V$$ 在 $\mathbb R^N$ 上的 Dunkl–薛定谔算子，其中 $$N>2$$ ， $$V$$ 是反向 Hölder 类 ${\rm RH}^q(dw)$ 中的非负势， $q>\mathbf N/2$ ，而 $\mathbf N$ 是齐次维数。我们建立了关于临界半径和轨道距离的预解核衰减估计。对于反射不变势，我们得到了 $$L^p(dw)$$ 上 $VL^{-1}$ 、 $L^{-1}$ 的二阶 Dunkl 导数以及 Riesz 变换 $\nabla_kL^{-1/2}$ 的有界性，并获得了分数阶势估计和混合估计。证明结合了弱解的局部估计、Dunkl Fefferman–Phong 不等式，以及针对欧氏距离和轨道距离的核估计。

5. 涨落动力学理论的一项研究纲领

A Program for Fluctuating Kinetic Theory

作者: Zhengyan Wu

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

发布日期: 2026-10-06

arXiv ID: 2610.08247v1

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

摘要翻译

中文标题: 涨落动力学理论的一项研究纲领

摘要翻译: 本文提出一项涨落动力学理论研究纲领，以统一研究动力学偏微分方程的非线性 Lévy 型涨落修正。Vlasov–Fokker–Planck (VFP)、Landau 和 Boltzmann 方程均可分解为一个动力学输运项，以及由一对非线性对偶耗散势生成的广义梯度结构。我们超越线性 Onsager 结构和经典涨落–耗散关系，提出一种新的非线性 Lévy 型涨落–耗散关系，用以识别与该广义梯度结构相关的热涨落。在一个正则倾斜类上，广义梯度结构的典范作用量与相应 SPDE 所预言的 Freidlin–Wentzell 型作用量一致。这种一致性为新的非线性涨落–耗散关系提供热力学支持，并导出一个主涨落动力学方程。我们进一步为非线性对偶耗散势建立一个充分条件，在该条件下，其中心化累积量具有 Lévy–Khintchine 表示；并证明所得的 Lévy 流可由高斯白噪声和补偿 Poisson 随机测度实现。作为这一框架的应用，我们导出一个 α-stable Lévy 型 Dean-Kawasaki 方程。此外，通过在噪声强度与流体力学缩放参数满足特定关系时对涨落 Boltzmann 方程进行重标度，我们证明其 Hamiltonian 收敛到 Landau–Lifshitz–Navier–Stokes 方程的 Hamiltonian。这揭示了动力学涨落与涨落流体力学之间的一种缩放关系。

6. 变指数空间中的分数阶梯度与非局部偏微分方程

Fractional gradients and nonlocal PDEs in variable exponent spaces

作者: Susana Moura, Cornelia Schneider

类别: math.AP, math.FA

发布日期: 2026-10-06

arXiv ID: 2610.08253v1

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

摘要翻译

中文标题: 变指数空间中的分数阶梯度与非局部偏微分方程

摘要翻译: 我们建立了变指数 Lebesgue 空间中分布意义下的 Riesz 分数阶梯度理论，并引入相应的非局部 Sobolev 空间。此外，我们建立了 Sobolev、Hardy、Morrey 和 Trudinger 型不等式，并推导出分数阶微积分基本定理。作为应用，我们研究具有非标准增长的非线性变分问题，并证明了极小化子的存在性及其 Euler–Lagrange 方程。

7. 界面阻尼下线性麦克斯韦方程组的适定性、正则性与稳定性

Well-posedness, regularity and stability for linear Maxwell’s equations with interface damping

作者: Julian Dörner, Serge Nicaise, Roland Schnaubelt

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08148v1

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

摘要翻译

中文标题: 界面阻尼下线性麦克斯韦方程组的适定性、正则性与稳定性

摘要翻译: 我们研究一个长方体区域上的线性麦克斯韦方程组，该区域被一个平面界面分为两部分，切向磁场跨越该界面发生跃变，且跃变量与切向电场成正比。这类条件出现在原子级薄超材料的建模中，其中薄片被替换为等效面电流。在允许界面上存在表面电荷的状态空间中，我们证明该问题由一个压缩半群生成。为此，我们给出底层函数空间的稠密性结果、韦伯型不等式以及迹估计。我们还证明了一个隐藏正则性结果，即麦克斯韦算子定义域中的场具有分片一阶 Sobolev 正则性。由于阻尼仅作用于界面，并非所有初始场的能量都会衰减。我们将动力学分解为两部分：一部分由对阻尼不可见的特征模态承载，且具有酉演化；另一部分是强稳定的。在无折射且界面位置非共振的情形下，我们通过两个矩条件刻画能衰减的初始场。

8. 欧几里得空间 R^n 中有界区域的庞加莱常数估计

Estimates of the Poincare constants for bounded domains in Euclidian space R^n

作者: Bernd Rummler, Gudrun Thäter

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08028v1

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

摘要翻译

中文标题: 欧几里得空间 R^n 中有界区域的庞加莱常数估计

摘要翻译: 我们研究有限维数 n 大于 1 的空间 R^n 中有界、单连通、凸且 C2-光滑的区域。对于这类区域，我们利用直径寻找包含该区域的球，并使球的直径取最小值。通过一个简单的变换，我们能够将所研究的区域映射到开单位球内的一个区域。我们推导了这些区域上标量函数和向量函数的庞加莱常数估计。主要工具分别是具有齐次狄利克雷迹的拉普拉斯算子或斯托克斯算子的第一特征值（我们采用 [21] 中的结果）。此外，我们还应用了标量拉普拉斯算子与斯托克斯算子的第一特征值及一个特征函数之间的关系，参见 [21]。我们得到的庞加莱常数估计是新的。空间维数 n 对估计的精度有显著影响。我们通过一个例子对此作了说明。

9. $\mathbb{T}^6$ 上的二次薛定谔系统

The quadratic Schrödinger system on $\mathbb{T}^6$

作者: Yilin Song, Yan Xue, Ruixiao Zhang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07821v1

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

摘要翻译

中文标题: $\mathbb{T}^6$ 上的二次薛定谔系统

摘要翻译: 本文研究 $\mathbb{T}^6$ 上二次薛定谔系统的全局适定性。为此，我们采用 Ionescu 和 Pausader \cite{IPT3} 开创的集中紧致性/刚性论证。我们的结果也可视为 Gao-Meng-Xu-Zheng \cite{GMXZ} 结果的周期性对应。

10. 五阶KdV方程中快速振荡的出现

The emergence of rapid oscillations in a fifth-order KdV equation

作者: Christopher J. Lustri, S. Jonathan Chapman, Richard C. P. Nicotra

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07797v1

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

摘要翻译

中文标题: 五阶KdV方程中快速振荡的出现

摘要翻译: 我们研究瞬时启动的五阶Korteweg–de Vries方程（KdV5）中快速传播、振幅呈指数小量级的振荡的出现，并证明这些波源于小时间边界层中的Stokes现象。我们对问题进行时间重标度，并构造一个辅助线性偏微分方程（lKdV5），以再现相关的指数渐近结构；随后采用两种主要方法进行分析。首先，通过最速下降法分析该辅助偏微分方程的傅里叶变换解，以提取相关的渐近贡献并确定Stokes曲线的位置。其次，我们通过直接的阶乘-幂次分析再现这些结果。在这一分析中，我们利用初始条件和远场条件，避免了指数渐近方法通常所需、而对于非线性偏微分方程又往往无法进行的局部渐近匹配。随后，我们将类似的阶乘-幂次分析应用于完整的非线性偏微分方程；该方程无法通过最速下降法研究。所得结果与数值模拟吻合良好，既展示了不衰减振荡的快速传播，也展示了振幅在远场衰减的瞬态涟漪，并揭示了支配其出现的复杂Stokes结构（包括常规Stokes曲线和高阶Stokes曲线）。这种发生在时间边界层中的复杂Stokes结构是一种新现象。

11. 相对论 Vlasov–Poisson 系统动量矩的局部传播

Local Propagation of Momentum Moments for the Relativistic Vlasov–Poisson System

作者: Zhaopeng Wang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07738v1

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

摘要翻译

中文标题: 相对论 Vlasov–Poisson 系统动量矩的局部传播

摘要翻译: 本文证明了三维相对论 Vlasov–Poisson 系统中阶数为 $k\ge5/2$ 的动量矩的局部传播。在库仑情形下，这改进了此前的结果；此前的结果要求阶数为 $k\ge3$ 的整数阶矩。证明采用了我们近期关于经典 Vlasov–Poisson 方程的工作中发展的临界密度方法。我们引入了一个与相对论速度相适配的加权动能张量，并得到临界密度估计。随后，该估计使我们能够通过短时间 bootstrap 闭合矩的论证。

12. 局部吸引作用下聚集-扩散的定量稳定性

Quantitative Stability Of Aggregation-Diffusion Under Local Attraction

作者: Yunhong Deng

类别: math.AP, math-ph

发布日期: 2026-10-05

arXiv ID: 2610.07264v1

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

摘要翻译

中文标题: 局部吸引作用下聚集-扩散的定量稳定性

摘要翻译: 本文证明了紧支撑初值下聚集-扩散能量的定量稳定性。我们重点考察具有局部吸引作用、且在无穷远处渐近平坦的退化扩散情形，例如 Riesz 核。我们刻画了稳态处 Wasserstein Hessian 谱下界。所需假设仅限于保证唯一性所必需的条件。

13. 对流–扩散方程的大自相似解

Large self-similar solutions of a convection–diffusion equation

作者: Xiaoxin Zheng, Dominika Pilarczyk, Grzegorz Karch, Marco Cannone

类别: math.AP

发布日期: 2026-10-05

arXiv ID: 2610.07259v1

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

摘要翻译

中文标题: 对流–扩散方程的大自相似解

摘要翻译: 我们为定义在整个欧几里得空间上的一类非线性对流–扩散方程构造自相似解，其初始数据属于一类齐次函数。相应的自相似剖面满足 $$$$ 上的非线性椭圆方程。我们通过带有截断非线性项的有界区域上的 Dirichlet 问题的解的极限，得到这些剖面；所用先验估计关于区域大小和截断水平二者均一致。我们的方法受到近期为 Navier–Stokes 系统构造大自相似解工作的启发。这里考虑的对流–扩散方程与 Navier–Stokes 方程具有某些结构特征，但其标量抛物型性质使我们能够利用最大值原理等工具，从而更直接地构造自相似解。

14. Ahlfors 正则边界区域中的分数阶光滑边界数据 Dirichlet 问题

The Dirichlet problem with boundary data of fractional smoothness in domains with Ahlfors regular boundary

作者: Ariel Barton, Svitlana Mayboroda, Alberto Pacati

类别: math.AP

发布日期: 2026-10-05

arXiv ID: 2610.07225v1

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

摘要翻译

中文标题: Ahlfors 正则边界区域中的分数阶光滑边界数据 Dirichlet 问题

摘要翻译: 我们证明，在边界具有 Ahlfors 正则性且维数可能非整数、满足 $0<\bdmn<\dmn$ 的弱局部 John 域 $Ω\subset\R^\dmn$ 中，对于系数为实数的线性二阶椭圆方程，当数据属于 Besov 空间且参数取值范围极广时，Dirichlet 问题是适定的。

不相关论文（基本信息）

1. 具有极小边界的正 Yamabe–Escobar 问题的紧性与非紧性

Compactness and Noncompactness for the Positive Yamabe–Escobar Problem with Minimal Boundary

作者: Santiago Cordero-Misteli

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08759v1

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

摘要: We study the Yamabe–Escobar equation prescribing positive constant scalar curvature and minimal boundary. The compactness conjecture asks whether being conformally equivalent to the round hemisphere is the only obstruction for its set of positive solutions to be compact. In this paper we show that this is true for dimensions $3\leq n\leq 14$ . We show that this dimensional threshold $$n_*=14$$ is sharp by constructing in $n\geq 15$ nonumbilic and non-locally-conformally-flat examples of metrics in the hemisphere for which a noncompact sequence of solutions exists. We do the same for the special classes of locally conformally flat metrics, $n_*^{LCF}=18$ , and the metrics with umbilic boundary, $n_*^{umb} =20$ . In all three settings we also prove that the Weyl–umbilicity vanishing conjecture holds for $4\leq n\leq n_*$ and construct counterexamples in $$n>n_*$$ by perturbing the noncompactness examples. To the best of the author’s knowledge, these are the first locally conformally flat noncompactness constructions for the positive-scalar-curvature, minimal-boundary case of the Yamabe–Escobar equation, outside the conformal class of the round hemisphere.

2. 含有不受限制的空间变反扩散项的 Kuramoto–Sivashinsky 方程的快速 Fredholm 稳定化

Rapid Fredholm stabilization of the Kuramoto–Sivashinsky equation with unrestricted, spatially-varying anti-diffusion

作者: Luke Bhan, Miroslav Krstic, Yuanyuan Shi

类别: eess.SY, cs.LG, math.AP, math.OC

发布日期: 2026-10-06

arXiv ID: 2610.08764v1

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

摘要: We develop the first feedback design for rapid stabilization of the Kuramoto–Sivashinsky equation with a spatially varying anti-diffusion coefficient. For constant coefficients, the single-input Fredholm design of Coron and Lü (2015) excludes a discrete set of values at which repeated unstable eigenvalues cause a loss of controllability. We overcome this obstruction by introducing a second boundary input and assigning the two inputs distinct roles. The key idea, inspired by Heymann’s Lemma, is to use the boundary value $$u(0,t)$$ entirely for a pre-feedback that renders the modified plant controllable through the curvature input $u_{xx}(0,t)$ . The latter input then stabilizes the plant through a Fredholm backstepping transformation. We show that two inputs suffice for controllability and are necessary when the plant has an unstable double eigenvalue. However, the Fredholm kernel still must be approximated for implementation. Hence, to enable kernel and gain approximation, we prove continuity of the coefficient-to-gain design map on compact admissible design classes. Unlike Volterra-based continuity proofs using successive approximations, our proof uses the modal representation to control the spectral data, the inverse coefficient system, and the tails of the kernel and gain series. This yields a single neural operator approximation of the gain to any prescribed $$L^2$$ accuracy across the class. Finally, we establish rapid local stabilization of the nonlinear closed-loop system under both the exact gains and sufficiently accurate approximations. We conclude with numerical results that illustrate prescribed decay rates and the computational cost of the approximations. In particular, we train a Fourier neural operator that achieves typical relative gain errors of approximately $0.1\%$ and stabilizes all held-out cases tested, including a plant with an unstable double eigenvalue.

3. 关于狄拉克算子的法贝尔–克拉恩猜想

On the Dirac Faber-Krahn Conjecture

作者: Habib Ammari, Jiayu Qiu

类别: math.SP, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08665v1

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

摘要: In this paper, we prove the strict Faber-Krahn inequality for non-negative-mass Dirac operators on simply connected domains with infinite-mass boundary condition. Denoting the first positive eigenvalue of the corresponding Dirac operator with mass $m\geq 0$ as $λ_1^{+}(Ω;m)$ , our main result shows that \begin{equation*} λ_1^{+}(Ω;m)\geq λ_1^{+}(B_{|Ω|};m), \end{equation*} where $B_{|Ω|}$ is the disk with the same area as the domain $$Ω$$ , and the equality is strictly valid if and only if $Ω=B_{|Ω|}$ . The key observation leading to the proof is that the upper component of the Dirac spinor associated with $λ_1^{+}(Ω;m)$ is non-vanishing over $$Ω$$ , which is proved by exploiting its Dirac current and the associated stream function. Based on this global non-vanishing property, we construct a sphere-valued spinor map of degree one, which allows a sharp comparison with the radial disk profile using the planar isoperimetry.

4. 非线性杆件连杆模型：零能量链与连续体一致性

A nonlinear rod linkage model: zero-energy chains and continuum consistency

作者: Dennis Kriventsov, Philippe G. LeFloch

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08705v1

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

摘要: We study a discrete Lagrangian model for a one-dimensional chain of rigid rods coupled by elastic springs. We characterize the geometry of its zero-energy constraint and distinguish the spinner and flipper regimes. Under graph-type hypotheses, selected zero-energy chains converge, after rescaling, to trajectories of an autonomous ordinary differential equation. For the physical linkage, we identify a simple alternating selection criterion under which chains form a unique monotone heteroclinic family in the spinner regime. In the flipper regime, we determine the topology of the constraint set and the number of successors, but leave global branch selection open. We also derive one- and two-component continuum residual expansions, distinguishing consistency from convergence of time-dependent solutions. Reproducible numerical experiments test the selected spinner recurrence and its asymptotic multiplier; for a simplified one-component model with smooth periodic data, they recover first- and second-order discrepancies for the leading and corrected continuum equations, respectively.

5. 三维立方体中的相对等周问题：分类与定量稳定性

Relative isoperimetry in the three-dimensional cube: classification and quantitative stability

作者: Bin Deng, Jiahuan Li, Yilu Liu

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08644v1

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

摘要: We determine the complete relative isoperimetric profile of the three-dimensional unit cube and classify all minimizing sets, including the equality cases at both transition volumes. Up to cube isometries and complementation, the minimizers are corner eighth-balls, edge quarter-cylinders, and coordinate slabs, with transitions at $$4π/81$$ and $$1/π$$ . Building on this classification, we prove that the perimeter deficit at each fixed volume controls the squared $$L^1$$ distance to the full minimizing family, with optimal exponent two. Near either transition, a uniform estimate combines the distance to each adjacent phase with its excess perimeter and identifies the linear degeneration, at neighboring volumes, of the best constant for distance to the actual minimizers. As consequences, we resolve the isoperimetric conjecture in Milman’s formulation (Conjecture~~1.1) for all cuboids $(0,β)\times(0,1)^2$ , $0<β\leq1$ , and the three-dimensional case of Brezis’s Open Problem~~10.1.

6. 关于巴赫张量与二次曲率泛函

On the Bach tensor and quadratic curvature functionals

作者: Letizia Branca, Davide Dameno

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08628v1

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

摘要: We study the critical points of a quadratic functional depending on the gradient of the Bach tensor on Riemannian four-manifolds, which generalize the Bach-flat condition. We show that, on every closed four-manifold, there exists a weak Bach-parallel metric, i.e. a critical point for this functional with respect to conformal variations: in particular, we prove that there exist infinitely many conformal classes which contain a unique minimizer for the functional, up to constant positive rescaling. Next, we analyze the global minima of the functional, i.e. metrics with parallel Bach tensor, relating these metrics to well-known variational problems. Using a version of de Rham’s splitting theorem on complete four-manifolds, we provide a classification result for products of surfaces, exploiting the theory of conformal gradient solitons; we also construct a new explicit example of a Bach-flat metric which is neither locally conformally flat nor conformally Einstein and we characterize HCMU metrics on complete surfaces. Finally, we prove an equivalence between the Bach-parallel condition on 4D cylinders and the existence of critical metrics for a well-known quadratic curvature functional in dimension three: in this direction, we also prove a characterization of flat three-manifolds, under some curvature and finite energy assumptions.

7. 闭流形与欧氏空间乘积上的大体积分数阶等周不等式

Large-volume fractional isoperimetry on products of closed manifolds and Euclidean spaces

作者: Juneyoung Seo

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08556v1

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

摘要: Let $$M$$ be a closed connected Riemannian manifold and let $\mathcal P_s$ be the fractional perimeter on $M\times\mathbb R^k$ defined through its heat semigroup, where $k\ge1$ and $$0<s<1$$ . We prove that the straight cylinder $M\times B_R$ , with $B_R\subset\mathbb R^k$ a ball of radius $$R$$ , uniquely minimizes $\mathcal P_s$ among all measurable sets of the same volume, up to translations in $\mathbb R^k$ and null sets, whenever $R\sqrt{λ_1(M)}$ is at least an explicit threshold depending only on $$k$$ and $$s$$ . Here $$λ_1(M)$$ is the first positive eigenvalue of the nonnegative Laplace–Beltrami operator $$-Δ_M$$ on $$M$$ . Under a larger threshold, the perimeter deficit controls both the deviation of the slice volumes from their mean and the square of the asymmetry relative to straight cylinders. For flat tori this gives large-volume uniqueness of straight cylinders for the periodic fractional perimeter, including singly periodic sets in every ambient dimension. The proof combines spectral reduction on $$M$$ with rearrangement of the Euclidean slices and a one-variable comparison between an interaction gain and a concavity loss.

8. 归一化哈密顿椭圆系统：关于双线性质量的 Gagliardo-Nirenberg 不等式

Normalised Hamiltonian Elliptic Systems: a Gagliardo-Nirenberg inequality for bilinear mass

作者: Daniele Cassani, Giulio Romani

类别: math.AP, math.FA

发布日期: 2026-10-06

arXiv ID: 2610.08774v1

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

摘要: We study Hamiltonian elliptic systems with prescribed bilinear mass $\int_{\mathbb{R}^N}uv=a>0$ . We establish a Gagliardo–Nirenberg inequality in the crossed gradient pairing and bilinear mass, on scaling-invariant component-bounded classes. It identifies the critical curve $$1/p+1/q=N/(N+2)$$ and the restricted energy trichotomy. For power nonlinearities, we establish uniqueness up to common translations and nondegeneracy of positive profiles throughout the Sobolev-subcritical hyperbola. Exact scaling then classifies positive normalised solutions and identifies the unique critical mass. In dimension two, we develop a bilinear exponential Gagliardo–Nirenberg estimate with a sharp gradient threshold and a truncated scalar refinement with the optimal quartic leading coefficient. Adapting variational ideas from [Cassani-Tarsi, Calc.Var.PDE (2015)] we combine a mass-normalising quotient, reduction over the full negative fibres and compactness below the concentration threshold to obtain least-energy positive radial solutions for pure exponential nonlinearities at every prescribed mass below the cubic limiting mass. Local perturbative branches and a small-frequency exponential branch yield explicit mass-response formulas.

9. Allen–Cahn方程稳定解的层片性质与刚性

Sheeting property and rigidity for stable solutions to the Allen–Cahn equation

作者: Yong Liu, Tianci Luo, Kelei Wang, Juncheng Wei, Yong Wei, Ke Wu

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08499v1

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

摘要: We prove a sheeting theorem for stable solutions to the Allen–Cahn equation in ambient dimensions $$2$$ through $$10$$ . We show that weak convergence of the measures $\varepsilon|Du_\varepsilon|^2\,dX$ to a hyperplane with integer multiplicity implies that the zero sets are disjoint smooth graphs. As a consequence, stable entire solutions in $\mathbb R^N$ , $N\le7$ , with $O(R^{N-1})$ energy growth are one-dimensional. When $$W''(-1)=W''(1)$$ , every blowdown of an entire solution with finite Morse index and this energy growth has multiplicity one for $4\le N\le10$ . For even potentials, we also obtain multiplicity one for limits with bounded energy and Morse index on closed manifolds of dimensions $$3$$ through $$7$$ , under bumpiness or positive Ricci curvature.

10. 具有 $H_\infty$ 误差界的偏微分方程传递函数稳定有理逼近

Stable Rational Approximation of PDE Transfer Functions with $H_\infty$ Error Bounds

作者: Aleksandr Talitckii, Matthew M. Peet

类别: math.AP, eess.SY

发布日期: 2026-10-06

arXiv ID: 2610.08497v1

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

摘要: Transfer functions of Partial Differential Equations (PDEs) are irrational and difficult to obtain. Unlike for rational transfer functions, only a few methods are available for the control and analysis of irrational transfer functions. Thus, for robust analysis and control of PDEs, we need to use a rational approximation of the transfer function, preferably with provable $H_\infty$ error bounds. In this paper, we extend Krylov subspace model reduction to infinite-dimensional systems modeled by Partial Integral Equations (PIEs) – an equivalent representation of PDEs. First, we propose a state transformation of the PIE system that allows us to compute a controllability matrix. Then, to obtain a rational approximation, we project the PIE states onto the vector space spanned by the truncated controllability matrix. Next, we show that sufficient conditions for exponential stability and well-posedness of the PIE system guarantee the exponential stability of the reduced-order system and provide constructive finite-frequency $H_\infty$ error bounds. Finally, we provide four illustrative numerical examples to demonstrate the performance of the proposed method.

11. 凸集上扭转函数梯度最大值的形状优化

Shape optimisation of the maximum of the gradient of the torsion function on convex sets

作者: Michiel van den Berg, Dorin Bucur, Mickael Nahon

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08408v1

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

摘要: We analyse the boundary behaviour of the gradient of the torsion function $w_\Om$ in convex domains $\Om\subset\mathbb{R}^d$ . We prove that its $L^\infty$ -norm is attained at a boundary point and establish the stability of this norm under general (convex) perturbations. For problems of the form

\sup\{\|\nabla w_{\Om}\|_\infty: \Om \subset \R^d,\textup{open, bounded and convex}, G(\Om)=m\},

where $$G$$ is either volume or perimeter, we show the existence of maximisers which are $$C^1$$ -smooth, rotationally symmetric about an axis and contain a $$(d-1)$$ - dimensional ball (of estimated size) in their boundary. Other constraints, whether geometric (such as diameter or circumradius) or variational (such as torsion or first Dirichlet eigenvalue), are also briefly discussed, and some of their characteristic behaviours are highlighted.

12. 一个奇异椭圆问题的极值参数、多重性与非存在性

Extremal parameters, multiplicity and nonexistence for a singular elliptic problem

作者: J. C. de Albuquerque, M. L. M. Carvalho, Edcarlos D. Silva, M. L. Silva

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08389v1

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

摘要: We consider the following class of singular elliptic problems

\begin{cases} -Δu = -\dfrac{u}{|u|^{β+1}}χ_{\{|u|>0\}} +λ|u|^{p-1}u, & \text{in }Ω,\\[1mm] u=0, & \text{on }\partialΩ, \end{cases}

where $Ω\subset\mathbb R^N$ is a bounded smooth domain, $N\geq3$ , $$0<β,p<1$$ , and $$λ>0$$ . The energy functional is not of class $$C^1$$ in $$H_0^1(Ω)$$ , and the standard critical point theory does not apply. We prove that a nontrivial nonnegative weak solution exists if and only if $$λ$$ exceeds a critical value $\barλ$ , which we identify and locate strictly between two Rayleigh-type parameters $0<λ_*<\barλ<λ^*<+\infty$ . This gives a sharp existence threshold. We also prove that for every $$λ$$ in a neighborhood on the left of $$λ^*$$ , there exist at least two distinct nontrivial nonnegative weak solutions. The first is a minimizer on the Nehari manifold, with an energy changing sign at $$λ^*$$ . The second has a strictly positive energy. For $$λ<λ^*$$ sufficiently close to $$λ^*$$ , both solutions have positive energy. This contrasts with previous multiplicity results for the same problem, in which the two solutions obtained have energies of opposite signs.

13. 连续中值滤波格式的强收敛

Strong Convergence of Continuous Median Filter Scheme

作者: Fabius Krämer, Tim Laux

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

发布日期: 2026-10-06

arXiv ID: 2610.08484v1

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

摘要: We introduce a continuous-time counterpart of the median filter that gradually denoises a given image. In the limit of vanishing stencil size, our results show that the evolution converges to level-set mean curvature flow. Surprisingly, and for the first time for any median-type filter scheme, we can prove the convergence of energies in this limit. Such strong convergence results are of major interest as they often have to be assumed in the literature to prove convergence of related schemes. Moreover, this result carries crucial geometric information, namely the unit multiplicity of interfaces. The proof relies on a compensated compactness argument inspired by the work of Evans and Spruck.

14. 各向同性 Ohta–Kawasaki 能量与 Yukawa 能量中层状结构的全局极小性

Global minimality of lamellae for isotropic Ohta–Kawasaki and Yukawa energies

作者: Dario Reggiani

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08485v1

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

摘要: We establish a sharp comparison principle for the sharp-interface Ohta–Kawasaki energy and its Yukawa-screened counterpart in dimensions $2\le d\le5$ , for $$m=0$$ , and under periodic and zero-flux Neumann boundary conditions. Whenever the finite optimal stripe width is compatible with the domain, we prove global minimality of periodic lamellae and characterize all minimizers. The proof combines a positive radial comparison measure, a sharp Fourier bound, and rigidity for a directional translation inequality in $$BV$$ . We also identify an obstruction to this measure construction in dimensions $d\ge6$ .

15. Monge-Ampère障碍问题自由边界直至边界处的正则性

Regularity of Free Boundary for Monge-Ampère Obstacle Problem up to the Boundary

作者: Hanfei Han, Genggeng Huang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08351v1

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

摘要: In this paper, we study the regularity of free boundaries for Monge-Ampère obstacle problem up to the boundary. We mainly focus on how the growth order of boundary value near the contact point impact the regularity. We also find an interesting phenomena that under quadratic growth boundary value, the free boundary enjoys better regularity in high dimensions than in two dimension.

16. 平均曲率流过程中拓扑复杂性的增加

Increase in topological complexity along the mean curvature flow

作者: Natasha Diederen

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08357v1

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

摘要: We construct an example of a smooth embedding of $\mathbb{S}^{p+q-1}$ in $\mathbb{R}^{p+q}$ whose evolution under the mean curvature flow forms an isolated singularity at one point, after which it is a smoothly embedded copy of $\mathbb{S}^{p-1}\times \mathbb{S}^q$ . In doing so, we prove some general results about $SO(p)\times SO(q)$ -invariant mean curvature flow whose initial data is generated by the rotation of a graphical profile curve, including an analysis of bubblesheet singularities. The most involved part of our analysis is an instant smoothness result, which requires a delicate approximation and pseudolocality.

17. 薄膜-体域中的热传输：有限厚度薄膜与薄膜极限

Thermal Transmission in a Film-Bulk Domain: Finite-Thickness Film and Thin-Film Limit

作者: Ferdinando Auricchio, Mihaela Chiappetta, Andrea Poiatti, Alessandro Reali

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08281v1

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

摘要: We investigate a thermal transmission problem in a domain composed of a film region and a bulk region coupled through an interface where temperature continuity and heat-flux balance are imposed. The analysis addresses both the finite-thickness film configuration and the corresponding thin-film limit, with anisotropic, non-homogeneous thermal conductivity tensors. Letting the film thickness vanish directly leads to a degenerate problem due to the vanishing measure of the film region. To preserve the contribution of the film in the limit, we introduce suitable scaling assumptions for the normal heat fluxes, namely for the flux balance across the film-bulk interface and for the flux prescribed on the external boundary of the film region, together with the requirement that the volumetric source in the film be uniformly bounded with respect to the film thickness. Under these scaling assumptions, we derive a non-trivial thin-film limit formulation that takes the form of a coupled bulk-interface problem, where the effect of the collapsing film is represented through an additional balance relation posed on the interface. Well-posedness results are established for both the finite-thickness film and thin-film limit configurations, and a rigorous convergence analysis is carried out to justify the asymptotic limit. The main mathematical achievement consists in finding explicit convergence rates in suitable functional norms, relying only on the weak regularity of the finite-thickness problem due to the transmission conditions at the film interface. To complement the study, analytical solutions are derived for a simplified film-bulk geometry and are used to compare the thermal responses of the two configurations. These solutions also provide a quantitative assessment of the convergence rates established in the theoretical analysis.

18. 利用紧性方法获得 Hessian 商方程凸解的 Hessian 估计

Hessian Estimates for Convex Solutions of Hessian Quotient Equations via Compactness

作者: Jianxiang Liu

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08442v1

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

摘要: We give a compactness proof of interior Hessian estimates for convex solutions of $σ_n(D^2u)/σ_{n-2}(D^2u)=1$ and $σ_n(D^2u)/σ_{n-1}(D^2u)=f(x,u,Du)>0$ , with $$1/f$$ concave in the gradient variable in the latter case, revisiting results of Lu and Fung.

19. Robin边界条件下平面条带中受迫曲率流的渐近动力学

Asymptotic Dynamics of Forced Curvature Flow in a Planar Strip with Robin Boundary Conditions

作者: Tianhong Pu, Xiaoqian Xin, Lixia Yuan, Xiaoyu Zhang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08270v1

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

摘要: We study the graphical forced mean curvature flow in a planar strip subject to Robin boundary conditions. The normal velocity is prescribed by $$V=H+b$$ , where $$b>0$$ is constant and $$H$$ denotes the curvature of the graph. We consider positive even solutions and reduce the problem to the half interval. The main results describe the degeneration produced by the Robin condition: the height, the instantaneous velocity, and the spatial gradient diverge, while the normalized slope $$u_x/u$$ converges to one away from the symmetry axis. The normalized velocity $$u_t/u$$ converges to the forcing parameter $$b$$ . The curvature converges to zero in the interior in an integral sense, whereas its horizontal curvature measure concentrates at the symmetry axis.

20. 带有逐点状态约束的一类反应-扩散系统的最优控制

Optimal Control of a Class of Reaction-Diffusion Systems with Pointwise State Constraints

作者: Eduardo Casas, Konstantinos Chrysafinos, Mariano Mateos

类别: math.OC, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08114v1

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

摘要: We study optimal control problems governed by a system of semilinear reaction-diffusion equations. The assumptions on the governing system are formulated so that our results apply to several models, including the Brusselator, Schnakenberg, and Gray–Scott systems. We prove the existence and uniqueness of a solution to the governing system without imposing any non-negativity assumption on the data. This allows us to study unconstrained, control-constrained, and state-constrained optimal control problems. For all these problems, we derive first-order necessary and second-order sufficient optimality conditions. The control acts on the activator equation through the boundary. Numerical examples demonstrate that both the inhibitor and the activator can be efficiently controlled.

21. Robin边界条件下的径向幂次平均曲率流

Radial Powered Mean Curvature Flow with Robin Boundary Conditions

作者: Tianhong Pu, Xiaoqian Xin, Lixia Yuan

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08283v1

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

摘要: We study a positive radial graph evolving by $$V=H^α+b$$ , with $$α>0$$ and $$b<0$$ , in the unit cylinder and with $$u_r(1,t)=u(1,t)$$ . We separate finite-time continuation estimates from the time-uniform interior estimates needed for convergence of normalized translates. On existing classical intervals we prove height and gradient bounds and an explicit conversion of velocity bounds into estimates for $$H$$ , $$u_r/r$$ , and $u_{rr}$ . We derive the correct boundary defect for comparisons with finite-slope translators. Under independent interior velocity bounds and two-sided slope trapping, we prove compactness and identify every translated limit with the cup translator. We identify the radial zero-number comparison still required to establish that trapping. Stationary caps and power-law Robin laws are included.

22. 任意快增长非线性项下半线性分数阶 Hermite 热方程的禁闭效应

Confinement effects for semilinear fractional Hermite heat equations with nonlinearities of arbitrarily rapid growth

作者: Divyang G. Bhimani, Mohamed Majdoub, Ramesh Manna

类别: math.AP, math.FA

发布日期: 2026-10-06

arXiv ID: 2610.08074v1

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

摘要: We study the Cauchy problem for a semilinear heat equation driven by the fractional harmonic oscillator. We establish local and global well-posedness in Orlicz and Morse–Transue spaces, based on sharp semigroup estimates and an integral compatibility condition between the Young function and the nonlinearity. We identify a critical Young function separating local well-posedness from instantaneous nonexistence for suitable data. Finally, we prove a confinement effect: under a Dini condition at the origin, small data yield global solutions decaying at the sharp exponential rate, with no Fujita-type restriction, in contrast with the corresponding free equation.

23. 无穷多个四维流形上的正 Einstein 度量

Positive Einstein metrics on infinitely many four-manifolds

作者: Tristan Ozuch

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08052v1

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

摘要: We construct smooth Einstein metrics with positive scalar curvature on $\#_k({\mathbb S^2}\times {\mathbb S^2})$ and on $k\mathbb{CP}^2\#(2k+1)\overline{\mathbb{CP}}{}^{2}$ , for $k\geqslant2$ . In particular, we show that positive Einstein metrics exist on infinitely many homeomorphism types of closed $$4$$ -manifolds. The first family shows that any simply connected spin $$4$$ -manifold carrying metrics with positive scalar curvature is homeomorphic to an Einstein $$4$$ -manifold. Wick rotation gives associated Lorentzian static axisymmetric Einstein metrics with positive cosmological constant and arbitrarily many horizons with the same surface gravity. The metrics are obtained by varying the cone angle $$2πβ$$ of a prescribed singularity along a torus, starting from the product metric on $\mathbb S^2\times\mathbb S^2$ and from the symmetric Kähler-Einstein metric on $\mathbb{CP}^2\#3\overline{\mathbb{CP}}{}^2$ . At $β=\frac{2}{k+1}$ , changing the period of the angle around the torus removes the cone singularity and changes the topology. For both families, as $k\to+\infty$ , the metrics collapse along the torus orbits to an explicit $$2$$ -dimensional disk, whose boundary has Hausdorff dimension $\frac{4}{3}$ . At the scale of curvature, they converge to complete periodic Ricci-flat metrics with Kasner asymptotics: the Euclidean Myers-Korotkin-Nicolai metric for the first family, and a new one for the second.

24. 一种用于具有淬灭现象的混合局部-非局部抛物型系统的移动平面法

A method of moving planes for a mixed local-nonlocal parabolic system with quenching

作者: Sergio Junquera, Leandro M. Del Pezzo

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08014v1

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

摘要: In this paper we develop a parabolic version of the method of moving planes on systems with a mixed local and nonlocal diffusion operator, that being $\mathcal{L} = -(-Δ)^s + Δ$ , where $$(-Δ)^s$$ is the integral fractional laplacian. We then use this result to study the quenching phenomena occurring in said systems with a singular absorption term of the type $(1-u)^{-p}$ . We prove that there exists a range of multiplicative parameters for which every solution presents quenching, and that for any initial datum $$u_0$$ , there exists a small enough multiplicative parameter such as the solution is globally defined. We also prove that the quenching phenomena cannot happen arbitrarily close to the boundary, but they are contained in a compact subset of the domain.

25. 紧致 Kähler 流形上复 Hessian 方程的插值方法

The interpolation method for the complex Hessian equation on compact Kähler manifolds

作者: Jiaogen Zhang, Xi Zhang

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07852v1

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

摘要: In current work, we delve into the study of complex Hessian equations on compact Kähler manifolds. By imposing a constraint on both the bisectional curvature and the oscillation of the solution, we demonstrate that the Laplacian estimate can be linearly controlled by its gradient estimate. Consequently, we are able to derive gradient estimate for the solution using a standard interpolation argument, without the need to employ the blow-up technique. As an application, we demonstrate that when the right-hand side of the complex Hessian equations approximate 1 in the $$L^1$$ norm, the desired gradient estimate is obtained.

26. Lévy符号收敛下半线性外域狄利克雷问题解的稳定性

Stability of Solutions to Semilinear Exterior Dirichlet Problems under Convergence of Lévy Symbols

作者: Tomasz Klimsiak, Andrzej Rozkosz, Leszek Słomiński

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07963v1

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

摘要: We study stability of semilinear exterior Dirichlet problems driven by possibly nonsymmetric Levy operators. Let $$L^n$$ and $$L$$ have Lévy-Khintchine symbols $$ψ_n$$ and $$ψ$$ , respectively. Our basic assumption on the operators is that $ψ_n(ξ)\toψ(ξ)$ for every $ξ\in\mathbb R^d$ . Under mild regularity assumptions on the domain and an absolute-continuity condition for the resolvent of the limiting process, we prove pointwise and $$L^p$$ -stability of solutions. No symmetry, uniform ellipticity, comparability of Lévy kernels, or common energy space is required; nonzero exterior data are allowed; and the monotone nonlinearity is subject to no growth restriction in the solution variable. The framework therefore covers local, nonlocal and mixed operators and, in particular, nonlocal-to-local limits, as well as fractional, relativistic, anisotropic and nonsymmetric stable-type operators. The proofs combine convergence of Lévy processes with a detailed analysis of exit times and exit positions and a stability theorem for monotone backward stochastic differential equations with varying terminal times. In the symmetric case, we also obtain Mosco convergence of the corresponding part forms on $$C^0$$ domains and give a counterexample showing that whole-space Mosco convergence need not pass to part forms on arbitrary open sets.

27. 基于后向散射数据的内部逆散射问题：谱比较与唯一性结果

The Interior Inverse Scattering Problem from Backscattering Data: Spectral Comparison and Uniqueness Results

作者: Xiaodong Liu

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.08011v1

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

摘要: We study an interior inverse scattering problem for the two-dimensional Helmholtz equation with backscattering data, where the source and receiver coincide. Such data are common in applications, yet existing uniqueness results are almost exclusively established for full multistatic data with separated sources and receivers. We consider two main measurement regimes. In the first regime, data are collected on a closed curve inside the unknown domain. Exploiting the meromorphic dependence of the backscattering field on the wavenumber, we show that the residues at the poles determine the diagonal of the spectral projection kernel on the curve, and hence the \emph{visible} part of the Dirichlet spectrum. This yields an unconditional spectral comparison: the symmetric difference of the two Dirichlet spectra is contained in the spectrum of the interior domain enclosed by the curve. If, in addition, the common visible eigenvalues have either matching multiplicities or matching restriction ranks, we obtain an area estimate bounding the discrepancy of the two domain areas by the area of this interior domain. In the second regime, data are given on a nonempty open subset of the domain, and a single non-resonant frequency already uniquely determines a bounded connected Lipschitz domain. The proof combines interior real-analyticity of the regular part of the Helmholtz Green function with its boundary blow-up.

28. 平坦环面上多个格林函数的临界点结构及其应用

Structure of critical points of multiple Green functions on flat torus and applications

作者: Zhijie Chen, Chang-Shou Lin

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07830v1

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

摘要: Let $E_τ:=\mathbb{C}/(\mathbb{Z}+\mathbb{Z}τ)$ be a flat torus and $$G(z)=G(z;τ)$$ be the Green function on $$E_τ$$ with the singularity at $$0$$ . Consider the multiple Green function $G_{n}$ on $(E_τ)^{n}$ :

G_{n}(z_{1},\cdots,z_{n}):=\sum_{i<j}G(z_{i}-z_{j})-n\sum_{i=1}^{n}G(z_{i}).

A critical point is called trivial if $\{z_{1},\cdots,z_{n}\}=\{-z_{1},\cdots,-z_{n}\}$ in $$E_τ$$ . Lin and Wang (Ann. Math. 2010) proved that $$G_1=-G$$ has exactly three trivial critical points and at most one pair of nontrivial critical points (depends on the choice of $$τ$$ ). For general $n\geq 2$ , Chai, Lin ang Wang (Camb. J. Math. 2015) proved that $$G_n$$ has at most $$2n+1$$ trivial critical points, but how many nontrivial critical points might exist remained completely open there. In this paper, we prove that $$G_n$$ has at most $$n$$ pairs of nontrivial critical points, and this upper bound is optimal at least for $$n=2,3$$ . The key idea is to show that nontrivial critical points (if exist) are always non-degenerate and contribute the same degree $(-1)^{n+1}$ . Applications to the curvature equation $Δu+e^{u}=8πn δ_{0}$ and the integral Lamé equation $y''=[n(n+1)\wp(z)+B]y$ will be given. In particular, we show that the kernel of the linearized operator is one-dimensional for any solution of the curvature equation.

29. 非线性分数阶演化方程解的渐近伪周期性

Asymptotic pseudo-periodicity of solutions for nonlinear fractional evolution equations

作者: Jin Liang, Yunyi Mu, Gaston Mandata N’Guérékata, Ti-Jun Xiao

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07914v1

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

摘要: This work focuses on the dynamic behaviors, particularly the asymptotic periodicity, of solutions to nonlinear fractional-order dynamical systems modeled by Caputo fractional evolution equations in Banach spaces. We propose a new class of weighted pseudo $$S$$ -asymptotically $$(ξ,τ,a)$$ -periodic functions, which unifies well-established function types such as standard periodic functions and $$S$$ -asymptotically periodic functions as its special instances. We then systematically derive the fundamental properties of this new function family, and rigorously prove its completeness, corresponding convolution theorems, and nonlinear composition rules. Building on these newly established results, we demonstrate that under reasonable hypotheses, the mild solutions to a broad class of nonlinear fractional-order dynamical systems — both delay-free and with finite delay — in Banach spaces admit the (weighted pseudo) $$S$$ -asymptotically $$(ξ,τ,a)$$ -periodic property, which further implies that no bifurcation or chaotic behavior emerges in these fractional-order systems. The theoretical conclusions are finally validated by illustrative examples describing memory-dependent complex dynamic processes, supplemented with supporting numerical simulations. Additionally, the presented weighted pseudo $$S$$ -asymptotically $$(ξ,τ,a)$$ -periodic framework can be extended to functions with first-kind discrete discontinuities, so that systems with discontinuous nonlinear terms can be handled via the same analytical paradigm.

30. 基于 Monge-增长对的 Hellinger-Kantorovich 梯度流神经网络 JKO 格式

A Neural JKO Scheme for Hellinger-Kantorovich Gradient Flows via Monge-Growth Pairs

作者: Geuntaek Seo, Cheolhyeong Kim, Hwijae Son, Hyung Ju Hwang

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

发布日期: 2026-10-06

arXiv ID: 2610.07602v1

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

摘要: We develop a mesh-free neural JKO scheme for advection-reaction-diffusion equations with a gradient-flow structure in the Hellinger-Kantorovich (HK) geometry of unbalanced optimal transport. Each update is parametrized by a spatial map and a mass-changing factor, allowing spatial redistribution and local mass creation or loss to be treated jointly within a single variational step. Their cone action bounds the squared HK distance from above, yielding a sufficient condition for discrete energy dissipation through comparison with the identity pair. Minimizing the pair objective over all admissible pairs recovers the exact JKO minimum when the source and a minimizer have positive densities. We establish existence and mass bounds for JKO minimizers and, under additional assumptions, obtain positivity and regularity together with a discrete Euler-Lagrange equation and a metric-dissipation identity. The self-consistent chemical potential is then nonincreasing along an optimal map. There exist parametric pairs whose endpoint densities and objective values converge to those of an exact JKO minimizer, provided a regular-pair approximation hypothesis holds. Finally, we show that a primal-dual gap controls objective suboptimality and, for Boltzmann entropy, the $$L^1$$ density error, assuming exact-step regularity, positive-semidefinite interactions, and global dual feasibility. Numerical experiments examine pointwise agreement with the PDE, energy dissipation, and the roles of transport, reaction, and fully implicit interactions.

31. 常 $Q_{2N}$ -曲率问题的非紧性

Noncompactness for the constant $Q_{2N}$ -curvature problem

作者: Liuwei Gong, Seunghyeok Kim, Juncheng Wei

类别: math.DG, math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07992v1

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

摘要: For every integer $N\ge4$ , we construct a fixed smooth, non-locally-conformally-flat metric on the $$n$$ -dimensional unit sphere $\mathbb{S}^n$ for which the constant $Q_{2N}$ -curvature equation admits an $L^\infty$ -unbounded sequence of positive solutions. The construction applies for $n\ge2N+20$ when $$N=4,5$$ , $n\ge2N+19$ when $6\le N\le8$ , $n\ge2N+18$ when $9\le N\le17$ , and $n\ge2N+17$ when $N\ge18$ . For each $L \in \mathbb{N} \cup \{0\}$ , the metric can be chosen arbitrarily close to the round metric in the $$C^L$$ norm. Together with the known cases $$N=1,2,3$$ , this establishes noncompactness at every order, with dimension bounds that we expect to be optimal. We derive an explicit formula for the fixed-volume Hessian of total $Q_{2N}$ -curvature in transverse-traceless (TT) directions at closed Einstein metrics satisfying $\operatorname{Ric}=(n-1)g$ . Building on Juhl’s formulas, we identify this Hessian, for every $N\in\mathbb{N}$ and $$n>2N$$ , as a degree- $$N$$ polynomial in the Lichnerowicz Laplacian. For the metric perturbations generated by algebraic Weyl tensors in the construction, the quadratic reduced energy with fixed bubble center is a positive multiple of this Hessian restricted to a four-dimensional TT space. Differentiation with respect to the bubble scale then yields a finite matrix whose exact sign analysis establishes the required indefiniteness in the stated dimension ranges. After continuation to real $$n$$ , this matrix changes from negative definite to indefinite at $n_N=2N+a_*+c_*N^{-1}+O(N^{-2})$ as $N\to\infty$ , with $a_*\approx16.201871$ and $c_*\approx13.512880$ , explaining the eventual bound $n\ge2N+17$ .

32. 关于非紧流形上的调和形式空间与薛定谔热方程

On the space of harmonic forms and Schrödinger heat equations on noncompact manifolds

作者: Jiuyuan Cui

类别: math.AP, math.DG

发布日期: 2026-10-06

arXiv ID: 2610.07617v1

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

摘要: It is known that on complete Riemannian manifolds with nonnegative Ricci curvature, the space of harmonic 1-forms of polynomial growth is finite-dimensional. In this paper, we extend this finiteness phenomenon to some manifolds whose Ricci curvature is allowed to be negative. More precisely, under a scale-invariant $$L^2$$ -Sobolev inequality, a weak Hardy inequality, and quadratic Ricci curvature decay with a sufficiently small coefficient, we prove that the space of polynomial-growth harmonic 1-forms is finite-dimensional. We also construct a counterexample showing the necessity of the smallness condition. Our result is closely related to an earlier theorem of L.F.Tam \cite{tam1998note}, where the finite-dimensionality of polynomial-growth harmonic 1-forms was established under a faster curvature decay assumption, without any smallness condition on the coefficient. In comparison, our result treats the critical quadratic decay regime, at the expense of requiring the corresponding coefficient to be sufficiently small. As a byproduct, we obtain an analogous finite-dimensionality result for polynomial-growth solutions to the Schrödinger heat equation, which generalizes the finite-dimensionality results in \cite{lin2019ancient},\cite{colding2021optimal} and \cite{lin2026space} to some extent.

33. 给定运动界面下两相调和映射热流的短时存在唯一性

Short-time existence and uniqueness for two-phase harmonic map heat flow with a prescribed moving interface

作者: Xingyu Wang

类别: math.AP

发布日期: 2026-10-06

arXiv ID: 2610.07736v1

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

摘要: Let $$N^+$$ and $$N^-$$ be disjoint compact smooth submanifolds of a Euclidean space, and let $Λ\subset N^+\times N^-$ be a compact smooth embedded submanifold prescribing admissible pairs of one-sided traces. Given a smooth family of separating hypersurfaces $Γ_t\subset\mathbb T^d$ , we study harmonic map heat flows into $N^\pm$ in the two moving phases, subject to

(u^+,u^-)\inΛ, \qquad \bigl(-\partial_{ν_t}u^+,\partial_{ν_t}u^-\bigr) \perp T_{(u^+,u^-)}Λ \quad\text{on }Γ_t.

For $$q>d+2$$ , initial data in $W_q^{2-2/q}$ satisfying these conditions at $$t=0$$ generate a unique short-time strong solution in $W_q^{2,1}$ , which is smooth up to the moving interface for every positive time; smooth initial data satisfying the compatibility conditions of every order yield solutions that are smooth up to the initial corner. We double the paired map across $$Λ$$ by normal reflection, which turns the interface conditions into a Dirichlet problem near the interface, and match this problem with the two bulk phases by an overlapping space-time Schwarz map inside a Schauder fixed-point argument. When the interface evolves independently by mean curvature, we also obtain a blow-up alternative in terms of its curvature and the two phase gradients. As an application, we prove local smooth solvability for the periodic matrix-valued sharp-interface system appearing in the convergence theory of Fei, Lin, Wang, and Zhang (\emph{Invent. Math.} 233 (2023), 1–80).

34. 具有临界指数增长的 FitzHugh-Nagumo 系统的极小极大解

Minimax solutions for FitzHugh-Nagumo systems with critical exponential growth

作者: João Marcos do Ó, José Carlos de Albuquerque, Hugo Henryque Silva

类别: math.AP

发布日期: 2026-10-05

arXiv ID: 2610.07411v1

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

摘要: We study a class of FitzHugh-Nagumo type systems with critical exponential growth in a smooth bounded domain of $\mathbb{R}^{2}$ . The classical reduction method turns the strongly indefinite system into a scalar nonlocal problem whose quadratic part induces a norm equivalent to, but different from the standard Sobolev norm, so that the Trudinger-Moser threshold must be handled in this framework. We prove the existence of a ground state solution via a mountain pass argument combined with a Nehari level comparison. We also prove the existence, without imposing the monotonicity condition required in the ground state argument, of a nontrivial solution obtained through a saddle-point linking argument based on the spectral decomposition of $\operatorname{T}=-Δ+B$ .

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

1. 单连通有界区域中二维欧拉方程涡片结构持续性的几何效应#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 旋转对称环面区域中的 Seifert 型磁流体静力平衡态#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 临界横向低频控制下 $$#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 三维环面 $\mathbb T^3$ 上 Navier–Stokes 方程的 Hartman–Grobman 定理#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 纳维–斯托克斯方程的翻译失真：为何对 AI 自动形式化结果进行 Lean 验证并不能保证自然语言证明正确#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 二维 Navier–Stokes-Boussinesq 方程中 Poiseuille 流的稳定性阈值#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 具有流出边界条件的正压流的弱 BV 稳定性与无黏极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. 小外力作用下三维矩形喷管中常值状态稳态跨声速激波的稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

9. 非线性反问题中依赖误差下的后验一致性及其在时间演化方程中的应用#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 具有物理碰撞频率的 BGK 方程的大数据全局解#

摘要翻译#

2. 关于带扰动与不带扰动的 Hardy 算子的波算子#

摘要翻译#

3. 利用 Harnack 型方法获得退化非局部方程的瞬时正则化#

摘要翻译#

4. $L^p$ 估计：Dunkl–薛定谔算子#

摘要翻译#

5. 涨落动力学理论的一项研究纲领#

摘要翻译#

6. 变指数空间中的分数阶梯度与非局部偏微分方程#

摘要翻译#

7. 界面阻尼下线性麦克斯韦方程组的适定性、正则性与稳定性#

摘要翻译#

8. 欧几里得空间 R^n 中有界区域的庞加莱常数估计#

摘要翻译#

9. $\mathbb{T}^6$ 上的二次薛定谔系统#

摘要翻译#

10. 五阶KdV方程中快速振荡的出现#

摘要翻译#

11. 相对论 Vlasov–Poisson 系统动量矩的局部传播#

摘要翻译#

12. 局部吸引作用下聚集-扩散的定量稳定性#

摘要翻译#

13. 对流–扩散方程的大自相似解#

摘要翻译#

14. Ahlfors 正则边界区域中的分数阶光滑边界数据 Dirichlet 问题#

摘要翻译#

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

1. 具有极小边界的正 Yamabe–Escobar 问题的紧性与非紧性#

2. 含有不受限制的空间变反扩散项的 Kuramoto–Sivashinsky 方程的快速 Fredholm 稳定化#

3. 关于狄拉克算子的法贝尔–克拉恩猜想#

4. 非线性杆件连杆模型：零能量链与连续体一致性#

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

1. 单连通有界区域中二维欧拉方程涡片结构持续性的几何效应

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 旋转对称环面区域中的 Seifert 型磁流体静力平衡态

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 临界横向低频控制下 $$$$

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 三维环面 $\mathbb T^3$ 上 Navier–Stokes 方程的 Hartman–Grobman 定理

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 纳维–斯托克斯方程的翻译失真：为何对 AI 自动形式化结果进行 Lean 验证并不能保证自然语言证明正确

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 二维 Navier–Stokes-Boussinesq 方程中 Poiseuille 流的稳定性阈值

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 具有流出边界条件的正压流的弱 BV 稳定性与无黏极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 小外力作用下三维矩形喷管中常值状态稳态跨声速激波的稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

9. 非线性反问题中依赖误差下的后验一致性及其在时间演化方程中的应用

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 具有物理碰撞频率的 BGK 方程的大数据全局解

摘要翻译

2. 关于带扰动与不带扰动的 Hardy 算子的波算子

摘要翻译

3. 利用 Harnack 型方法获得退化非局部方程的瞬时正则化

摘要翻译

4. $$L^p$$ 估计：Dunkl–薛定谔算子

摘要翻译

5. 涨落动力学理论的一项研究纲领

摘要翻译

6. 变指数空间中的分数阶梯度与非局部偏微分方程

摘要翻译

7. 界面阻尼下线性麦克斯韦方程组的适定性、正则性与稳定性

摘要翻译

8. 欧几里得空间 R^n 中有界区域的庞加莱常数估计

摘要翻译

9. $\mathbb{T}^6$ 上的二次薛定谔系统

摘要翻译

10. 五阶KdV方程中快速振荡的出现

摘要翻译

11. 相对论 Vlasov–Poisson 系统动量矩的局部传播

摘要翻译

12. 局部吸引作用下聚集-扩散的定量稳定性

摘要翻译

13. 对流–扩散方程的大自相似解

摘要翻译

14. Ahlfors 正则边界区域中的分数阶光滑边界数据 Dirichlet 问题

摘要翻译

不相关论文（基本信息）

1. 具有极小边界的正 Yamabe–Escobar 问题的紧性与非紧性

2. 含有不受限制的空间变反扩散项的 Kuramoto–Sivashinsky 方程的快速 Fredholm 稳定化

3. 关于狄拉克算子的法贝尔–克拉恩猜想

4. 非线性杆件连杆模型：零能量链与连续体一致性