生成时间: 2026年09月28日 12:38:50

运行状态: 完成

计划处理论文数量: 46

已完成论文数量: 46

未完成论文数量: 0

重点关注论文数量: 8

了解领域论文数量: 21

不相关论文数量: 17

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

1. 具有 Kraichnan 输运噪声的随机二维 Euler 方程全离散近似的收敛性

Convergence of a fully discrete approximation for the stochastic 2D Euler equations with Kraichnan transport noise

作者: Abhishek Chaudhary, Ujjwal Koley, Andreas Prohl

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

发布日期: 2026-09-25

arXiv ID: 2609.31369v1

摘要翻译: 我们研究二维不可压 Euler 方程涡量形式的全离散近似理论，该方程定义在 $\mathbb{T}^{2}$ 上，受 Kraichnan 型输运噪声驱动，且初始涡量取值于 $\mathbb{H}^{-1}(\mathbb{T}^2)$。本文的主要新意在于针对环面 Kraichnan 模型建立了一个有限维强制性估计，该估计在离散层面刻画了粗糙输运噪声所诱导的正则化机制。结合一种新的紧性策略——该策略基于连续的方程型插值以及局部化的离散估计——这一强制性估计给出了近似序列的紧性，并证明其某个子列极限是 Euler 方程的弱鞅解。

Comment: 58 pages

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

1. 研究对象和背景

论文研究二维不可压 Euler 方程的涡量形式在粗糙输运噪声驱动下的全离散逼近。空间为二维环面 $\mathbb{T}^2=(\mathbb{R}/2\pi\mathbb{Z})^2$，涡量 $\omega$ 与速度 $u$ 满足 Biot-Savart 关系 $u=\nabla^\perp(-\Delta)^{-1}\omega$。确定性涡量方程为

$$ \partial_t\omega+(u\cdot\nabla)\omega=0, \qquad u=\nabla^\perp(-\Delta)^{-1}\omega. $$

本文考虑的随机模型是在输运项中加入 Kraichnan 型、无散度的空间随机场：

$$ d\omega+(u\cdot\nabla)\omega\,dt+\sum_{k\in I}(\sigma_k\cdot\nabla)\omega\circ dW_k(t)=0, \qquad u=\nabla^\perp(-\Delta)^{-1}\omega. $$

这里 $W_k$ 是相互独立的标准 Brownian motions，噪声向量场由环面上的 Fourier 模态显式给出。具体地，对 $m\in\mathbb{Z}^2\setminus\{0\}$，令 $e_m=m^\perp/|m|$、$q(m)=(1+|m|^2)^{-1-\alpha}$，其中 $0<\alpha<1$，并取

$$ \sigma_{m,c}(x)=\sqrt{2q(m)}\,e_m\cos(m\cdot x),\qquad \sigma_{m,s}(x)=\sqrt{2q(m)}\,e_m\sin(m\cdot x). $$

这些噪声不是作用于涡量的加性外力，而是改变涡量的输运；Stratonovich 形式保留了随机输运的几何结构。对本文采用的余弦、正弦模态，噪声的 Itô 修正产生与 $\Delta\omega$ 成比例的项。因而噪声一方面代表未解析小尺度对大尺度流动的影响，另一方面在负阶 Sobolev 能量估计中提供额外的耗散或正则化机制。论文的初值允许仅属于 $H^{-1}(\mathbb{T}^2)$，并要求零平均；这明显低于经典 Yudovich 理论中的有界涡量假设，也覆盖相当粗糙的涡量数据。

低正则性下，二维 Euler 方程的困难不只是构造形式上的解，还包括控制涡量的集中、获得紧性并识别非线性极限。既有理论分别处理了不同类别：Yudovich 在有界涡量类中建立适定性；DiPerna-Majda、Delort 与 Schochet 等工作发展了粗糙涡量弱解及极限识别中的紧性和集中抵消理论；确定性数值分析中，谱方法的收敛通常要求较高正则性，而谱黏性方法可处理包括 Delort 类在内的较粗初值。随机方面，粗糙 Kraichnan 输运噪声已被证明可以带来负阶 Sobolev 正则化，并已有随机弱解存在及特定数据类别下的唯一性结果。不过，连续方程的存在理论并不自动给出可实施的空间、噪声和时间同时离散的方案，也不自动保证离散解极限满足相应的随机弱解定义。

本文试图补足的正是这一数值分析与随机弱解构造之间的缺口。作者在环面上选取 Fourier 截断，使噪声模态能显式计算；同时截断空间频率、噪声频率并离散时间，构造全离散谱格式。该问题的关键在于：离散噪声必须在有限维空间中仍能提供足够的负阶 Sobolev 强制性，而格式还必须具有可用于紧性证明的稳定性。为此，论文证明环面 Kraichnan 模型的有限维一致强制性估计，并据此建立离散能量控制；随后用适应该格式的紧性策略，证明离散近似的子列极限是随机 Euler 方程的弱鞅解。论文明确强调，工作的结论是子列收敛，而非全序列收敛或弱解唯一性。

2. 主要定理或主要结果

  1. 有限维噪声强制性估计。令 $P_N$ 为去除零模并截断到 $|n|_\infty\leq N$ 的 Fourier 投影，令 $I_K$ 表示保留 $|m|_\infty\leq K$ 的余弦、正弦噪声模态。论文定义离散二次型
$$ D_{N,K}(\omega)=\sum_{k\in I_K}\left\|(-\Delta)^{-1/2}P_N[(\sigma_k\cdot\nabla)\omega]\right\|_{L^2(\mathbb{T}^2)}^2. $$

在 $0<\alpha<1$ 且尺度分离条件

$$ K\geq \sqrt{2}\,\gamma N,\qquad \gamma>1+\sqrt{2}, $$

下，定理式结论 (4.58) 给出该有限维二次型对相应 $H^{-\alpha}$ 量的统一强制性控制。也就是说，保留的噪声频率必须显著高于已解析的空间频率，才能在离散层面捕捉 Kraichnan 噪声诱导的负阶正则化。所提供的 PDF 摘录没有呈现 (4.58) 的完整不等式及其常数，因此不能可靠补写精确系数、范数权重的具体等价形式或常数对 $\alpha,\gamma$ 的依赖。

  1. 全离散格式与稳定性。格式 (3.2) 在 $V_N$ 中演化涡量，使用时间步长 $\Delta t$、空间截断 $N$ 和噪声截断 $K$：
$$ \omega_N^{n+1}-\omega_N^n+\Delta t\,P_N[(u_N^{n+1}\cdot\nabla)\omega_N^n] =\frac{\Delta t\,c_K}{2}\Delta\omega_N^{n+1} -\sum_{k\in I_K}P_N[(\sigma_k\cdot\nabla)\omega_N^n]\Delta^{n+1}W_k, $$
$$ u_N^{n+1}=\nabla^\perp(-\Delta)^{-1}\omega_N^{n+1}, \qquad c_K=\sum_{m\in M_K}q(m). $$

这里 $M_K=\{m\ne0:|m|_\infty\leq K\}$。混合时间层 $(u_N^{n+1},\omega_N^n)$ 使以 $(-\Delta)^{-1}\omega_N^{n+1}$ 测试时输运项恰好抵消；隐式 Laplacian 项则提供耗散控制。作者在 (5.2)、(5.3) 建立离散 $H^{-1}$ 能量估计及稳定性控制。由于摘录未给出这两个估计的全文，具体期望值上界、时间积分项和常数依赖无法从现有内容确认。

  1. 子列收敛及极限识别。定理 6.2 的结论是：对论文所构造的全离散近似，在相应截断及尺度条件下，可抽取子列，使近似收敛到 Itô 形式随机 Euler 方程的弱鞅解。该 Itô 形式包含噪声项及其由截断极限产生的 Laplacian 修正；极限弱解的定义见 Definition 6.1，弱形式由 (6.29) 给出。摘录未列出定理 6.2 的完整参数趋近假设，特别是 $\Delta t,N,K$ 是否需要满足除强制性条件以外的定量耦合关系，因此不宜声称某个未见于摘录的收敛速率或网格关系。

紧性证明使用连续方程型插值

$$ \omega_j^\sharp(t)=\omega_j^0+\int_0^tF_j(s)\,ds+M_j(t),\qquad t\in[0,T], $$

它在时间网格点与离散解一致。由局部化的 Hölder 控制，作者得到 $\omega_j^\sharp$ 在 $C([0,T];H^{-5}(\mathbb{T}^2))$ 中的紧性，并在包含连续插值、分段常值插值和 Brownian motions 的路径空间上建立联合紧性。文中进一步证明

$$ \omega_j^\sharp-\omega_j\longrightarrow0, \qquad \omega_j-\overline{\omega}_j\longrightarrow0 \quad\text{in probability in }L^2(0,T;H^{-4}(\mathbb{T}^2)), $$

其中不同插值符号对应论文中的连续和离散时间插值。经 Skorokhod-Jakubowski 表示后，相关插值收敛到同一极限过程。离散速度还沿同一子列强收敛到极限 Biot-Savart 速度，见 (6.24)；这一点用于识别非线性输运项。最后将极限传入离散弱形式 (6.21)，得到极限弱形式 (6.29)。所给摘录未提供速度强收敛的具体函数空间，故这里只能保留其强收敛事实，不能补充未显示的拓扑。

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

证明的首要障碍是如何把连续 Kraichnan 模型的正则化机制转化为有限维、可计算的离散估计。离散二次型 $D_{N,K}$ 包含 Fourier 投影、有限噪声模态以及 $(-\Delta)^{-1/2}$ 权重；其强制性不能仅由连续协方差结构直接推出。尤其是空间截断会改变频率相互作用，若噪声频率没有高于解析频率足够多，离散和式未必能一致控制所需的负阶 Sobolev 量。论文的核心新技术是对环面 Fourier 模态直接分析这一有限维二次型，而不是将连续结论形式上移植到离散格式。

该强制性证明由三部分构成。第一，Lemma 4.1 给出一族随方向变化的区域上的一致 Riemann 和逼近。其边界层测度满足统一估计，且被积函数具有统一 $C^1$ 控制，因此格点求和与连续积分之间的误差随尺度趋于零；证明中通过区域内部方格的求积误差和边界方格数量估计分离处理。第二，Lemma 4.2 证明连续缺陷积分为正。对方向 $e\in S^1$，论文构造 $\tau_e$ 和 $F_e(z)=|z|^{-2-2\alpha}(1-\tau_e(z))$，并证明主值积分与 $e$ 无关，且

$$ \operatorname{p.v.}\int_{|z|\leq\gamma}F_e(z)\,dz =j_{\alpha,\gamma} =2\pi\int_1^\gamma \rho^{-1-2\alpha}(1-\rho^{-2})\,d\rho>0. $$

这一步提供连续极限中的正量，而不是仅凭形式上的各向同性推断正性。第三，Lemma 4.5 的离散格点分析把该连续正性转回有限格点和式，并与一致 Riemann 和估计结合，得到 (4.23)，进而推出关键强制性估计 (4.58)。因此，真正起作用的是连续缺陷的严格正性、随方向一致的求积控制及格点频率几何的共同配合。环面几何与有限噪声截断在此处不可忽略；作者也明确将这一论证视为环面模型的特定分析。

第二个障碍是离散能量结构与随机紧性之间的衔接。格式把速度取在新时间层、涡量取在旧时间层，即采用 $(u_N^{n+1},\omega_N^n)$ 的混合输运项。这个选择并非一般性的时间离散便利：以 $\psi_N^{n+1}=(-\Delta)^{-1}\omega_N^{n+1}$ 测试时，作者证明

$$ \left\langle P_N[(u_N^{n+1}\cdot\nabla)\omega_N^n],\psi_N^{n+1}\right\rangle=0. $$

因此非线性项不会破坏 $H^{-1}$ 能量估计；隐式项 $\frac{c_K}{2}\Delta\omega_N^{n+1}$ 则以正确符号提供稳定性贡献。再结合噪声强制性，得到离散估计 (5.2)、(5.3)，为后续时间正则性、局部化和紧性提供基础。这里的设计同时服务于格式可实施性与证明中的结构抵消。

紧性策略是另一项实质创新。作者指出，已有 Ondreját、Prohl、Walkington 的抽象框架针对不连续 Skorokhod 型路径空间的分段常值插值，并要求类似全局一致界；本问题可获得的控制却是由离散停时产生的局部化控制，不能直接满足该框架的假设。因此作者不直接把分段常值近似作为紧性对象，而构造满足积分方程的连续插值 $\omega_j^\sharp$。停时局部化估计导出其在 $C([0,T];H^{-5})$ 中的局部 Hölder 控制，继而获得紧性。分段常值插值仍用于写离散弱形式，但作者通过在概率意义下的 $L^2(0,T;H^{-4})$ 渐近等价关系，保证三种插值最终具有同一极限。

具体极限论证还需要在非可分或弱拓扑相关的路径空间上处理联合紧性，因此使用 Skorokhod-Jakubowski 表示定理，而不是未经说明地套用仅适用于标准 Polish 空间的表示结论。表示后，极限过程对由极限涡量及 Brownian motions 生成的通常增广滤过适应。通过离散弱形式逐项识别漂移、噪声和修正项；其中涡量只在较弱的负阶空间中紧，而 Biot-Savart 算子额外改善速度的紧性，速度的强收敛使 $(u\cdot\nabla)\omega$ 的非线性项能够通过极限。摘录给出了这一论证路线及结论，但没有提供各项一致可积性和乘积极限的全部技术引理，故不宜据此虚构更精细的逐项估计。

4. 与之前工作的比较

与连续随机 Euler 理论相比，本文不是首先证明 Kraichnan 输运噪声能够正则化粗糙涡量，而是把这一机制落实到一个同时离散空间、时间和噪声的可实施方案中。Coghi 与 Maurelli 在 $\mathbb{R}^2$ 上对粗糙 Kraichnan 噪声驱动的随机二维 Euler 方程证明了 $\dot H^{-1}$ 初值下的弱鞅解存在，并得到额外的 $L^2([0,T];H^{-\alpha})$ 正则性；在 $L^1\cap L^p$ 初值下还建立了受参数条件限制的路径唯一性。Jiao 与 Luo 后来把该唯一性范围由 $p>3/2$ 推进到 $p>1$。Brzeźniak、Flandoli 与 Maurelli 处理有界涡量情形的存在唯一性，Brzeźniak 与 Maurelli 则证明了正 $H^{-1}$ 涡量的随机 Delort 型存在结果。上述工作主要是连续分析理论；本文新增的是可由数值离散逼近产生弱鞅解的收敛理论，而非宣称改进这些连续问题中的唯一性范围。

空间选择也体现了与既有工作不同的取舍。作者解释，在 $\mathbb{R}^2$ 上，噪声向量场缺少便于直接实现的显式表达，因而构造全离散方案面临困难；环面上的 Fourier 模态则允许显式定义、有限截断并直接计算。代价是结论针对 $\mathbb{T}^2$ 的具体周期模型，不能据此认为已经建立了 $\mathbb{R}^2$ 上同类可实施算法的收敛性。环面上显式噪声结构还使作者能够用格点和及连续缺陷积分证明新强制性估计；这既是技术优势，也意味着证明依赖特定几何与频率结构，推广到一般区域或非 Fourier 噪声并非自动成立。

与确定性粗糙 Euler 数值分析相比，经典 Fourier 谱方法在解足够光滑时可获得谱收敛，但这类正则性前提不能覆盖本文仅取 $H^{-1}$ 的初始涡量。Lopes Filho、Nussenzveig Lopes 与 Tadmor 研究了无集中效应的近似解；Lanthaler 与 Mishra 证明谱黏性方法可收敛至包括 Delort 类初值在内的粗糙 Euler 解。本文的目标并非给出光滑解情形的高阶误差界，也没有在摘录中给出 $N,K,\Delta t$ 的定量收敛速率，而是证明离散随机解存在收敛子列并识别其极限为弱鞅解。其贡献更接近构造性存在理论，而非经典意义下的误差分析定理。

相对于 Ondreját、Prohl 与 Walkington 的抽象离散随机方程收敛框架，本文明确指出既有框架要求的全局一致控制与当前格式可用的停时局部化控制不匹配。因此作者发展了以连续方程型插值为紧性对象的路径：连续插值具有预期时间正则性，分段常值插值保留离散弱形式所需的结构，再证明两者在相关拓扑中渐近等价。这不是简单更换紧性定理，而是针对局部控制重组紧性论证。其代价是需要处理多种插值、停时、联合紧性和 Jakubowski 表示，证明结构较为专门。

与欧氏空间中的 Kraichnan 协方差分析相比，作者将环面上的显式三角 Fourier 模态直接代入离散二次型，结合 Riemann 和一致逼近、连续缺陷积分正性及格点分析建立 (4.58)。作者称该环面有限维强制性估计据其所知是新的，并将其视作对 [5] 中协方差方法的环面对应物，而非同一估计的直接推论。尺度条件 $K\geq\sqrt{2}\gamma N$（$\gamma>1+\sqrt{2}$）明确要求噪声截断频率高于空间解析频率；这是捕捉高频噪声正则化所付出的分辨率条件。现有摘录没有展示该条件是否最优，也没有给出减弱该条件的结果，故不能断言其尖锐性。

最后，本文的结论受弱解理论本身的限制。作者明确说明由于这一正则性层级下缺乏完整唯一性理论，收敛结果采取子列形式；因此不能推出整个离散族收敛到唯一极限，也不能从本文声称路径唯一性或数值解对初值的稳定连续依赖。摘要和摘录强调的是存在弱鞅解及其由方案逼近的构造性意义。关于噪声正则化的物理解释、向其他几何或噪声模型推广的可能性，现有摘录未提供可据以断言的定量结果；作者未显式展开的比较应保留为未确认事项。

2. 二维 Navier-Stokes 方程中的激励与可辨识性

Excitation and Identifiability in the 2D Navier-Stokes equations

作者: Vincent R. Martinez, Sarah Strikwerda, Xiang Wan

类别: math.OC, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31230v1

摘要翻译: 参数估计的收敛性经典地要求“持续激励”，即对依赖于轨线的信号施加非退化条件。据我们所知，本文首次为非线性偏微分方程建立了此类激励理论。针对不可压流体二维 Navier–Stokes 方程中如何根据谱观测识别未知黏度的问题，我们提出的激励条件可计算、可先验验证，并且在缩放意义下是尖锐的。我们采用数据同化方法来处理未知初始状态，并通过最小化流体速度的低模观测与同化这些观测的 nudging 型滤波器的低模投影之间的损失，筛选候选黏度。我们框架的主要新意在于区分出一个新的泛函 $[?25l\mathsf{W}[?25h$，它表示滤波器所对应的灵敏度变量对涡量平方积分所做的功；我们的整个分析都围绕这一泛函展开。我们证明，$[?25l\mathsf{W}[?25h$ 与相应 Gauss–Newton 迭代的可观测性 Gramian 显式可比，并进而在观测损失的每个临界点建立如下二分结论：要么 $[?25l\mathsf{W}[?25h$ 超过某个阈值，此时候选黏度满足一个显式误差估计，该估计与 $[?25l\mathsf{W}[?25h$ 和观测密度 N 成反比，而与初始条件之间的误差成正比；要么 $[?25l\mathsf{W}[?25h$ 低于该阈值，此时观测对参数更新在定量意义下不敏感，且用户能够检测到这种不敏感性。值得注意的是，我们的方法基于能量，因此有望推广到许多其他非线性耗散系统。

Comment: 49 pages

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

1. 研究对象和背景

本文研究二维不可压缩流体中未知运动黏度的参数反演, 特别关注有限空间分辨率的速度观测能否保证黏度可辨识, 以及如何在观测轨迹上验证这一点. 背景模型是在周期区域 $\mathbb{T}^2=[0,2\pi]^2$ 上的均值为零、散度为零的速度场 $u$, 满足

$$ \partial_t u+\nu Au+B(u,u)=f,\qquad u(0)=u_0, $$

其中 $\nu>0$ 是未知运动黏度, $A=P_\sigma(-\Delta)$ 是 Stokes 算子, $P_\sigma$ 是 Leray 投影, $B(a,b)=P_\sigma(a\cdot\nabla)b$, 而 $f$ 为给定的无散外力. 观测为有限个 Fourier 模态在时间区间 $[0,\tau]$ 上的完整时间序列

$$ \mathcal{O}_{N,\tau}=\{P_Nu(t):0\leq t\leq\tau\}, $$

其中 $P_N$ 投影到波数 $|k|\leq N$ 的模态. 因而问题不是从完整流场恢复参数, 而是在初始状态 $u_0$ 未知、空间观测被截断的情况下, 利用单条观测轨迹估计一个标量参数 $\nu$.

未知初值会与未知参数相互混淆: 候选黏度不正确造成的模型误差, 可能在有限时间内被初始状态误差部分补偿; 反过来, 即使真实黏度已知, 若状态未正确初始化, 直接比较未同化模型轨迹与观测也未必得到小残差. 作者采用连续数据同化中的 nudging 机制, 对每个候选黏度 $\hat\nu$ 构造滤波状态 $v(t;\hat\nu)$:

$$ \partial_t v+\hat\nu Av+B(v,v)=f-\mu P_Nv+\mu P_Nu, $$
$$ v(0)=P_Nu(0)+Q_Nv_0,\qquad Q_N=I-P_N. $$

这里 $\mu>0$ 是用户选择的 nudging 增益, $v_0$ 指定未观测高模态的初值. 右端反馈只使用已观测的低模态, 旨在使滤波状态追踪参考流. 参数估计再通过观测损失

$$ J(\hat\nu;\tau)=\frac12\int_0^\tau\|P_Nv(t;\hat\nu)-P_Nu(t)\|_{L^2}^2\,dt $$

在可行区间 $\mathcal{A}=[\nu,(1+\rho)\nu]$ 上进行. 论文分析的不只是全局极小点, 还包括满足一阶约束最优性条件的候选集合 $\mathcal{C}_\tau$, 因此结论可用于判断优化过程可能停留的临界候选参数.

研究动机来自参数估计中的 persistency of excitation 条件. 在有限维系统中, 这类条件要求轨迹含有足够的信息, 使参数变化能在观测中产生可检测的响应; 缺少激发时, 参数到输出的映射可能局部或定量地近乎不变, 即便优化算法本身运行正常也无法可靠恢复参数. 对流体方程而言, 这一问题同时牵涉非线性演化、无限维状态、部分观测及初值误差. 既有二维 Navier–Stokes nudging 理论主要回答在观测分辨率和反馈强度足够时, 滤波状态是否同步或重构参考状态; 状态同步并不自动给出黏度参数的可辨识性, 更不提供可由能量量直接检验的参数激发阈值.

本文的出发点是把滤波状态对候选黏度的敏感度 $\partial_{\hat\nu}v$ 与涡量/涡度能量联系起来. 作者引入

$$ \mathsf{W}(\hat\nu;\tau)=\int_0^\tau(\Delta v(t;\hat\nu),\partial_{\hat\nu}v(t;\hat\nu))_{L^2}\,dt, $$

并注意到分部积分给出

$$ \mathsf{W}(\hat\nu;\tau)=-\frac{d}{d\hat\nu}\mathsf{E}(\hat\nu;\tau),\qquad \mathsf{E}(\hat\nu;\tau)=\frac12\int_0^\tau\|\nabla v(t;\hat\nu)\|_{L^2}^2\,dt. $$

因此 $\mathsf{W}$ 是滤波状态总 enstrophy 对候选黏度的负导数, 也可解释为敏感度变量对 enstrophy 所做的功. 论文试图填补的关键空白, 是把抽象的观测可辨识性与一个能由滤波轨迹和敏感度计算的能量量联系起来, 给出“激发充分时参数误差可控, 激发不足时观测对参数更新不敏感”的定量二分. 作者将结果定位为非线性 PDE 中 persistency-of-excitation 型理论的首次发展; 这一“首次”是作者在论文中的表述, 不应理解为对所有 PDE 参数识别工作的穷尽性文献判断.

2. 主要定理或主要结果

主要结论集中在 Theorem 3.11 及其前置的 Proposition 3.3. 论文考虑二维周期 Navier–Stokes 方程、有限 Fourier 模态观测及 nudging 滤波, 以参考黏度 $\nu$ 为真值, 候选值 $\hat\nu$ 取自 $\mathcal{A}=[\nu,(1+\rho)\nu]$. 所有估计均要求处于作者称为 detectability 的参数范围. 按文中给出的尺度关系, 该范围为

$$ \mu\gtrsim \nu N^2,\qquad N\gtrsim \frac{U}{\nu}, $$

其中 $U$ 是刻画参考流大小的量; 提供的摘录没有列出其精确定义及上述比较常数, 因而这里仅保留论文明确给出的尺度形式.

首先, 作者证明新的 work functional 与 Gauss–Newton 迭代的 observability Gramian 可比较. 记

$$ G(\hat\nu;\tau)=\int_0^\tau\|P_N\partial_{\hat\nu}v(t;\hat\nu)\|_{L^2}^2\,dt\geq 0. $$

敏感度方程的能量平衡将 $G$、$\mathsf{W}$、非线性项、阻尼项和终端敏感度联系起来; 在 detectability 条件下得到

$$ \frac{\mu}{2}G(\hat\nu;\tau) \leq \mathsf{W}(\hat\nu;\tau) \leq \mu G(\hat\nu;\tau)+O\!\left(\frac{V^2\tau}{\mu}\right), $$

其中 $V$ 表示滤波流的大小, 文中说明其可由与 $\mu$ 无关的先验界控制. 因而在该参数范围内, work functional 不仅非负, 而且在控制误差项意义下反映 Gramian 的非退化程度. 这一比较把参数到观测的微分响应, 即 $G$, 转换为基于滤波轨迹 enstrophy 的能量量 $\mathsf{W}$.

其次, Theorem 3.11 给出对每个一阶候选极小点 $\hat\nu\in\mathcal{C}_\tau$ 的激发二分. 令

$$ \eta_0=\|\nabla(v_0-u_0)\|_{L^2}. $$

存在无量纲阈值 $\omega_*>0$, 其显式依赖于 $\eta_0$, 可行区间相对宽度 $\rho$ 以及参考速度场大小, 但不依赖观测分辨率 $N$ 和具体候选值 $\hat\nu$. 若

$$ \mathsf{W}(\hat\nu;\tau)\geq\omega_*\tau, $$

则黏度误差与观测损失满足论文所示的估计

$$ (\hat\nu-\nu)^2\mathsf{W}(\hat\nu;\tau) \leq O\!\left(\frac{\eta_0^2}{N}\right), \qquad J(\hat\nu;\tau) \leq O\!\left(\frac{\eta_0^2}{N^2\mu}\right). $$

这表明在给定初值误差下, 激发越强, 对 $|\hat\nu-\nu|$ 的控制越好; 增大 $N$ 也改善误差界. 这里的 $O(\cdot)$ 记号按原文保留, 摘录未提供隐含常数的完整表达式, 不宜据此补写精确常数. 另一方面, 若

$$ \mathsf{W}(\hat\nu;\tau)<\omega_*\tau, $$

则有

$$ \frac1\tau\int_0^\tau\|\partial_{\hat\nu}v(t;\hat\nu)\|_{L^2}^2\,dt \leq O\!\left(\frac{\omega_*}{\mu}\right). $$

因此, 未达到阈值并不被解释成另一个参数恢复定理, 而是给出可检测的低敏感度诊断: 平均而言, 参数扰动对滤波状态的影响较小, 观测对参数更新不够敏感.

这些估计对整个候选集合 $\mathcal{C}_\tau$ 成立, 而非只对某个已知的局部极小点成立. 若定义激发候选集

$$ \mathcal{V}_\tau=\{\hat\nu\in\mathcal{A}:\mathsf{W}(\hat\nu;\tau)\geq\omega_*\tau\}, $$

则激发分支上的候选临界点满足上述误差估计. 在 $N\to\infty$ 的极限下, 论文据此推出: 当激发候选集非空时, $\mathcal{V}_\tau\cap\mathcal{C}_\tau$ 最终收缩到真值 $\nu$; 并且 $\nu\in\mathcal{M}_\tau$, 其中 $\mathcal{M}_\tau$ 是观测损失在可行区间上的极小点集合. 作者还强调, 在其所述的最小激发情形, 这一关于真值为极小点的结论不需要先取无限长观测窗 $\tau\to\infty$.

把 $\mathsf{P}(\hat\nu;\tau)=\mathsf{W}(\hat\nu;\tau)/\tau$ 作为平均 enstrophy 敏感度, 激发条件也可写作 $\mathsf{P}\geq\omega_*$. 文中据此区分有限数据下的实用判据与极限意义下的可辨识性判据: 对给定候选值检查 $\mathsf{W}\geq\omega_*\tau$ 或 $\mathsf{P}\geq\omega_*$, 可作为可计算的激发诊断; 若相应条件对 $N\to\infty$ 或 $\tau\to\infty$ 的候选集合一致成立, 则得到作者所称的 bona fide identifiability criterion. 此外, 只要存在一个候选临界点达到阈值, 即 $\mathcal{V}_\tau\cap\mathcal{C}_\tau\neq\varnothing$, 该候选点便可使用激发分支的参数误差界; 作者将这类存在性解释为 locatability criterion.

最后, Section 4 讨论阈值尺度的尖锐性. 在小振幅外力情形, 参考方程有唯一全局吸引稳态, 并且 nudging 方程及敏感度方程也有相应的唯一全局吸引稳态. 作者对稳态平均功率给出下界, 其量级与 $\omega_*$ 相匹配. 对 Kolmogorov shear forcing, 可显式构造参考稳态, 并证明前述下界在量级上尖锐至显式前因子; 还构造了低于激发阈值的情形, 其中敏感度变量的时间平均平方范数严格小于 $\omega_*\mu^{-1}$ 的相应尺度, 与二分结论中的低敏感度分支一致. 摘录没有呈现这些定理的完整假设和常数, 因而无法在此补列其精确参数条件.

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

证明的核心障碍是区分“滤波状态能跟踪观测”与“观测足以识别参数”. 前者由低模态 nudging 和 Foias–Prodi 型 detectability 控制; 后者要求候选黏度变化能在低模态观测中产生非退化响应. 直接使用 Gauss–Newton Gramian

$$ G(\hat\nu;\tau)=\int_0^\tau\|P_N\partial_{\hat\nu}v\|_{L^2}^2\,dt $$

固然自然, 但它涉及敏感度方程解的观测部分, 不容易直接由轨迹能量作先验检验. 本文的关键设计不是仅仅估计这个 Gramian, 而是找到与之相连、具有能量解释的 $\mathsf{W}$, 再围绕 $\mathsf{W}$ 建立全篇参数估计.

具体地, 对 nudging 方程关于 $\hat\nu$ 求导. 若记 $z=\partial_{\hat\nu}v$, 则其线性化方程包含黏度扰动产生的强迫 $-Av$, 线性化对流项 $B(z,v)+B(v,z)$, 黏性耗散 $\hat\nu Az$ 及观测反馈 $\mu P_Nz$. 初值对 $\hat\nu$ 不变时, $z(0)=0$. 对该方程与 $z$ 作 $L^2$ 内积, 可以得到敏感度能量平衡. 摘录明确列出的恒等式为

$$ \frac12|z(\tau)|^2+\hat\nu\int_0^\tau\|z(t)\|_{L^2}^2\,dt +\mu G(\hat\nu;\tau) =\int_0^\tau(B(z,v),z)\,dt+\mathsf{W}(\hat\nu;\tau). $$

这里 nudging 项恰好产生 $\mu G$, 参数微分的扩散项经过分部积分则形成 $\mathsf{W}=\int(\Delta v,z)$. 非线性项 $\int(B(z,v),z)$ 是能量平衡中需要控制的主要项; 作者在二维情形利用先验界、二维非线性估计以及 detectability 参数条件, 把它吸收或估计到可控误差中, 从而得到 $G$ 与 $\mathsf{W}$ 的比较. Proposition 3.3 给出的结果是

$$ \frac{\mu}{2}G\leq\mathsf{W}\leq\mu G+O(V^2\tau/\mu). $$

因此 $\mathsf{W}$ 的下界不仅给出灵敏度能量, 也表明滤波状态总 enstrophy 随候选参数变化的单调性强度; 当 $G\geq0$ 时, 该比较特别推出 detectability 范围内 $\mathsf{W}\geq0$.

第二个技术环节是把参数差 $\hat\nu-\nu$ 与滤波模型误差、观测残差联系起来. 对参考解与 nudged 解作差, 误差方程同时含有初值差、黏度差造成的项以及由反馈产生的低模态耗散. 作者先建立 model error 的精细表示和残差估计, 再将其与临界点条件结合. 其逻辑是: 一方面, 一阶最优性条件限制损失对参数方向的变化; 另一方面, 非零 $\mathsf{W}$ 通过敏感度能量恒等式保证参数变化确实在系统响应中留下可量化信号. 两者合用, 才能将黏度误差乘上激发量控制, 而不仅得到一个依赖不可观测状态的粗略界. 完整的 model-error 表示、残差引理及其精确函数空间假设在提供的 PDF 摘录中未出现, 因而这里不臆造其逐项公式.

第三个环节是阈值构造及其对候选点的一致性. 作者将初始状态不匹配量

$$ \eta_0=\|\nabla(v_0-u_0)\|_{L^2} $$

纳入估计, 并将可行参数区间宽度 $\rho$、参考流大小和该初值误差组合成阈值 $\omega_*$. 重要之处在于 $\omega_*$ 不依赖 $N$ 或具体候选 $\hat\nu$, 因此可在所有 $\hat\nu\in\mathcal{C}_\tau$ 上统一使用. 这使结论成为对整个候选临界集合的 dichotomy, 而非依赖已知真值邻域的局部逆函数定理. 当 $\mathsf{W}$ 超过阈值时, 误差界中 $\mathsf{W}$ 位于分母所体现的有效控制作用; 当低于阈值时, 则改为给出 $z$ 的时间平均小量估计, 将失败的识别条件转译成用户能够发现的敏感度不足.

锐性分析也有实质性的技术作用, 并非只提供示例. 在小振幅外力条件下, 作者先证明参考动力系统以及同化状态、敏感度系统存在唯一全局吸引稳态, 从而可以用稳态功率刻画长时间平均激发. 随后对线性化 Navier–Stokes 算子加 nudging 项建立新的 coercive 型估计, 用它推出稳态功率下界, 并显示下界与最小阈值具有相同尺度. 在 Kolmogorov 流这一可显式求解的剪切流设置中, 精确稳态结构使得这些尺度关系可以进一步检验, 并可展示低于阈值时敏感度确实偏小. 作者称相关 coercive 估计可能有助于其他 nudging 同化研究; 这是作者对潜在用途的判断, 而非本文已经证明的推广定理.

整体上, 证明框架依赖二维耗散方程的先验估计、敏感度方程的能量恒等式、低模态 nudging 的 detectability 尺度、模型误差与残差控制, 以及临界点的一阶条件. 其创新在于把这些要素围绕 enstrophy 灵敏度 $\mathsf{W}$ 组织成一条可量化的论证链: 参数微分产生敏感度方程, 敏感度能量平衡把 $\mathsf{W}$ 连到观测 Gramian, 临界点分析把该量转成参数误差界, 而低 $\mathsf{W}$ 分支则获得可检测的参数不敏感性结论. 论文声称方法以能量估计为基础, 有望迁移到其他非线性耗散系统; 但摘录中没有给出针对其他方程的完整验证.

4. 与之前工作的比较

与连续数据同化方面的工作相比, 本文把 Azouani–Olson–Titi 的 nudging 状态重构框架及相关分析从“未知初值下重构状态”推进到“未知初值与未知参数同时存在时识别黏度”. 既有工作 [7,17] 研究 nudging 或相关连续同化滤波的状态重构; [11] 等工作处理噪声, 另有研究扩展到不同强迫、观测方式或模型误差. 这些结果提供了观测模态数与 nudging 强度足够时状态误差衰减的基础, 但状态同步本身并不保证参数到观测映射的可逆性. 本文在相同的二维 Navier–Stokes 观测背景下新增参数敏感度分析及激发阈值, 将 detectability 与 identifiability 区分开来: $\mu,N$ 的尺度条件确保滤波器能利用观测, 而 $\mathsf{W}$ 是否超过 $\omega_*\tau$ 决定是否能对黏度误差作定量控制.

在参数迭代与 Gauss–Newton 分析方面, 作者承接了 [22] 中基于同化的黏度估计思路及其数值观察, 也联系 [72] 对相关 Newton 型迭代的收敛分析、[74,76] 对损失函数与 Gauss–Newton 结构的解释. 标准 Gauss–Newton 分母是观测 Gramian $G$, 其非退化是经典局部收敛和辨识分析的重要条件; 相关理论通常要求 Gramian 非退化, 但在 PDE 中这一条件直接验证并不容易. 本文的推进是证明在 detectability 尺度下 $\mathsf{W}$ 与 $G$ 可比较, 进而用滤波 enstrophy 对参数的导数作为能量化替代量. 这不是说 Gramian 不再重要, 而是提供了一个与 Gramian 相联系、较便于先验分析及从滤波轨迹评估的量. 与 [2] 对 Lorenz 63 有限维系统的同类优化与可辨识性研究相比, 本文处理的是无限维、非线性的二维 PDE, 并须同时控制状态模型误差和空间截断; 其代价是目前只处理一个未知标量参数.

与既有二维 Navier–Stokes 可辨识性结果 [15] 相比, 本文强调有限时间窗上的判据. 摘录指出 [15, Theorem 5.3] 的判据与到全局吸引子的距离有关, 且要求全局吸引子不包含零状态; 该文的框架建立在全局吸引子上, 自然使用双无限时间观测 $\{P_Nu(t):t\in\mathbb{R}\}$. 本文则通过同化状态和敏感度, 在有限区间 $[0,\tau]$ 上定义 $\mathsf{W}$ 与损失, 并给出基于有限窗轨迹的定量二分. 这使结果更适于有限数据诊断, 但仍有 detectability 条件、初值误差和激发阈值等限制; 不能据此断言任意有限时间观测均可辨识黏度.

在控制理论语境中, nudging 系统可视为 Luenberger observer, 观测算子为有限秩有界投影 $P_N$, 增益为 $\mu P_N$. 经典 detectability 说明反馈能否稳定状态误差; 文中联系 [7] 的 Foias–Prodi 型机制, 将足够大的 $\mu$ 和足够精细的 $N$ 解释为 PDE 情形的 detectability 条件. 传统有限维非线性 observer 的状态与参数估计研究, 如 [1,4,5,61,84], 在作者所述框架中仍主要面对有限维系统, 并且有些方法依赖状态到参数关系的线性化或秩条件. 本文面对的映射 $\hat\nu\mapsto v(\cdot;\hat\nu)$ 是非线性的, 因而需要从敏感度能量和非线性 PDE 估计中推导激发阈值. 另一方面, 本文观测算子 $P_N$ 有界且有限秩, 作者明确指出这避免了边界观测、无界观测算子理论中常见的精细 trace 正则性困难; 因而本文并未解决边界或局部观测下的同类参数反演问题.

与最优控制和强约束 4DVAR 相比, 常见做法是同时优化初值和参数, 并在损失中加入初值背景项或 Tikhonov 正则化, 以控制未观测高模态. 论文给出的对照泛函含有初值变量及惩罚项 $\lambda\|\hat u_0-P_Nu(0)\|^2/2$. 这种方法需要处理对无限维初值的导数, 通常需要伴随方程或额外正则化. 本文固定滤波器初值的观测低模态部分, 让 nudging 反馈处理未知初值带来的状态差异, 从而将优化变量缩减为一维黏度参数, 并避免在该损失中显式优化完整初值. 这带来计算和分析上的简化, 但没有消除对初值错配的依赖: $\eta_0$ 仍显式进入阈值和误差界, 且滤波初始高模态仍由用户选择的 $v_0$ 决定.

与逆问题理论相比, 作者把定量稳定性与单射性解释为条件良定性: PDE 正问题保证参数到观测映射有定义, 激发阈值提供条件性的唯一性, 参数误差估计则给出逆映射稳定性. 其机制与唯一延拓在稳定化问题中的作用有结构类比, 但不是通常基于局部空间观测的唯一延拓定理. 对 Fourier 投影观测, $P_N\hat v\equiv0$ 并不一般推出 $\hat v\equiv0$; 作者在 Kolmogorov 流讨论中强调这种逐点模态观测的唯一延拓式命题可能失败. 本文的贡献因此是以激发量给出定量、轨迹依赖的替代判据, 而不是证明谱观测具有一般唯一延拓性质.

相对于稳态黏度反演研究 [91], 本文对象为时间依赖轨迹及有限时间观测, 而 [91] 的结果在作者描述中建立于小振幅外力、稳态 Navier–Stokes 逆问题, 并分析适定性、稳定性及离散化下优化方案的收敛. 论文将 [91] 视为小外力 regime 中的数值伴随方向; 本文则在小外力情形下把阈值锐性作为更一般时间依赖理论的一个检验, 并进一步研究 Kolmogorov 流. 但本文的 PDE 层面定量结论并不等同于对任意数值离散方案的收敛证明; 摘录仅提到数值误差可相对直接纳入分析, 具体内容指向 Remark 3.15, 其细节从提供内容中无法确认.

适用范围和代价方面, 本文的明确理论范围是二维周期、均值为零的不可压 Navier–Stokes 方程, 单一未知常黏度, 有限 Fourier 模态观测, 且要求 nudging 参数与分辨率满足 detectability 尺度. 阈值依赖初值错配、可行参数区间和参考流大小; 观测分辨率虽不进入 $\omega_*$, 但进入误差估计, 而 $\mu$ 同时影响滤波与敏感度界. 论文没有在所给内容中证明多参数情形、三维 Navier–Stokes、含噪观测下的同等激发定理, 也没有给出对一般观测算子都成立的结论. 作者明确把多参数扩展列为未来工作. 对这些方向的可推广性, 只能说能量方法预计可能适用, 不能视为已经建立的结果.

3. 临界 Fourier-Herz 空间中不可压 Navier-Stokes 方程解的渐近剖面

Asymptotic profiles for solutions to the incompressible Navier-Stokes equations in the critical Fourier-Herz spaces

作者: Ikki Fukuda, Ryosuke Nakasato

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30902v1

摘要翻译: 我们研究全空间 $mathbb{R}^d$ ($d\ge 2$) 上不可压 Navier-Stokes 方程初值问题解的大时间渐近行为。当初始速度属于 $L^1(\mathbb{R}^d)$ 时，众所周知，由于散度自由相容条件，其空间积分为零。考虑到这一性质，Carpio (1996) 和 Fujigaki-Miyakawa (2001) 推导出了强解的高阶渐近公式。本文研究初始速度属于 Fourier-Herz 空间 $\widehat{L}^1(\mathbb{R}^d)$ 时解的大时间行为；与 $L^1(\mathbb{R}^d)$ 相比，该空间提供了更宽泛的框架，并在低频区域施加了额外条件。特别地，在二维情形下，我们证明会出现一种不同于 Carpio 和 Fujigaki-Miyakawa 所得到的渐近剖面。

Comment: 20 pages

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

1. 研究对象和背景

论文研究全空间 $\mathbb{R}^d$, $d\ge 2$ 上不可压 Navier–Stokes 初值问题的长时间渐近行为. 未知量为速度 $u(x,t)$ 与压力 $p(x,t)$, 方程为

$$ \begin{cases} \partial_t u-\mu\Delta u+\operatorname{div}(u\otimes u)+\nabla p=0,\\ \operatorname{div}u=0,\\ u|_{t=0}=u_0, \end{cases} $$

其中 $\mu>0$ 为黏性系数. 对速度作 Helmholtz 投影后, 解满足 mild 公式

$$ u(t)=e^{\mu t\Delta}u_0-\int_0^t e^{\mu(t-\tau)\Delta}P_\sigma\operatorname{div}(u\otimes u)(\tau)\,d\tau. $$

核心问题是刻画 $t\to\infty$ 时解的首项及次阶项, 并判断非线性项如何改变线性热扩散给出的渐近剖面. 经典研究常用初值的空间衰减, 例如 $u_0\in L^1$ 或带权 $\langle x\rangle u_0\in L^1$, 以定义初值的空间矩并从低频 Taylor 展开导出渐近公式. 但不可压约束对这些矩有特殊影响: 若 $u_0\in L^1$ 且 $\operatorname{div}u_0=0$, 则 $\widehat u_0$ 在原点连续, 且 $\xi\cdot\widehat u_0(\xi)=0$, 因而 $\int_{\mathbb{R}^d}u_0(x)\,dx=\widehat u_0(0)=0$. 这使得通常的热核质量项消失, 需要进一步分析更高阶矩和非线性贡献.

本文转而使用傅里叶–Herz 空间. 记 $\dot\Delta_j$ 为 Littlewood–Paley 频率块, 则

$$ \|f\|_{\dot B^s_{p,\sigma}(\widehat L)} =\left\|\left\{2^{sj}\|\dot\Delta_j f\|_{\widehat L^p}\right\}_{j\in\mathbb Z}\right\|_{\ell^\sigma}, \qquad \|f\|_{\widehat L^p}=\|\widehat f\|_{L^{p'}}. $$

尺度变换 $u_\lambda(t,x)=\lambda u(\lambda^2t,\lambda x)$ 下, $\dot B^{-1}_{\infty,2}(\widehat L)$ 是本文使用的临界数据空间. 既有结果已在该空间的小初值条件下建立全局 mild 解、解析性和时间衰减; 本文利用其中的全局控制, 研究更精细的渐近展开. 作者还要求 $u_0\in\dot B^0_{1,\infty}(\widehat L)$, 这是用于获得所需衰减估计的附加条件.

本文要处理的关键低频信息不是空间积分, 而是指定 $\widehat u_0$ 在 $|\xi|\to0$ 时相对于 Helmholtz 投影的行为. 具体地, 假设存在非零向量 $M_0\in\mathbb R^d$, 使

$$ \lim_{|\xi|\to0}|\widehat u_0(\xi)-P(\xi)M_0|=0, \qquad P(\xi)=I-\frac{\xi\otimes\xi}{|\xi|^2}. $$

更强的假设要求某个 $\alpha>0$ 满足

$$ |\widehat u_0(\xi)-P(\xi)M_0|=O(|\xi|^\alpha),\qquad |\xi|\to0. $$

由于 $P(\xi)$ 随方向变化, $P(\xi)M_0$ 一般在原点没有方向无关的连续延拓. 因此, 该条件表达的是广义低频矩, 而非 $L^1$ 情形中的普通积分矩. 对应的线性剖面是 Oseen 核作用于 $M_0$:

$$ S_\mu(t)M_0=\mathcal F^{-1}\!\left[e^{-\mu t|\xi|^2}P(\xi)M_0\right]. $$

这项研究的重要性在于, 它把渐近展开从依赖物理空间衰减和可积矩的框架, 推进到允许较弱空间衰减、直接以低频 Fourier 行为描述初值的框架, 并揭示二维中二次非线性产生的临界对数修正. 本文的结论仍建立在小临界范数、额外的 $\dot B^0_{1,\infty}(\widehat L)$ 条件及相应低频假设之上, 并非对任意临界空间初值都成立.

2. 主要定理或主要结果

以下结论针对小初值全局 mild 解. 定理 2.1–2.2 假设 $u_0\in\dot B^{-1}_{\infty,2}(\widehat L)\cap\dot B^0_{1,\infty}(\widehat L)$, $\operatorname{div}u_0=0$, 且 $\|u_0\|_{\dot B^{-1}_{\infty,2}(\widehat L)}$ 足够小. 还假设上述弱低频条件 (A), 并对更强结论使用带速率的条件 (B). 参数范围 (C) 为: $p=1$ 时 $s\ge0$; $1<p<\infty$ 时 $s>-\min\{d/p',1+d/p\}$; $p=\infty$ 时 $s\ge-1$.

定理 2.1 给出首阶剖面. 在条件 (A) 和 (C) 下, 对 $t\to\infty$,

$$ \lim_{t\to\infty}t^{\frac d2(1-1/p)+\frac s2} \big\||\nabla|^s(u(t)-S_\mu(t)M_0)\big\|_{\widehat L^p}=0. $$

因此, 主导项是 $S_\mu(t)M_0$, 而不是经典 $L^1$ 展开中的热核质量项. 若再假设 (B), 即低频误差为 $O(|\xi|^\alpha)$, 定理 2.2 将收敛改进为定量余项: 对 $t\ge1$, 当 $d\ge3$ 时

$$ \big\||\nabla|^s(u(t)-S_\mu(t)M_0)\big\|_{\widehat L^p} \lesssim t^{-\frac d2(1-1/p)-\frac s2-\frac{\min\{\alpha,1\}}2}; $$

当 $d=2$ 时

$$ \big\||\nabla|^s(u(t)-S_\mu(t)M_0)\big\|_{\widehat L^p} \lesssim t^{-\frac d2(1-1/p)-\frac s2} \left(t^{-\alpha/2}+t^{-1/2}\log(1+t)\right). $$

定理 2.4 给出 $d\ge3$ 的二阶展开. 在定理 2.2 的假设之外, 要求 $\alpha>1$ 及 $1<p\le\infty$. 则

$$ \lim_{t\to\infty}t^{\frac d2(1-1/p)+\frac{s+1}{2}} \left\||\nabla|^s\left(u(t)-S_\mu(t)M_0 +\partial_kS_\mu(t)\int_0^\infty\!\int_{\mathbb R^d}(u_k u)(x,\tau)\,dx\,d\tau\right)\right\|_{\widehat L^p}=0. $$

换言之, 在此符号约定下, 次阶项为 $-\partial_kS_\mu(t)$ 乘以非线性时间积分. 它比首项多衰减 $t^{-1/2}$; 积分振幅编码了整个演化过程的非线性作用.

定理 2.6 处理 $d=2$. 在定理 2.2 的假设之外, 要求 $\alpha\ge1$; 此时对 $1\le p\le\infty$,

$$ \big\||\nabla|^s\big(u(t)-S_\mu(t)M_0+(\log t)\partial_kS_\mu(t)A_k[M_0]\big)\big\|_{\widehat L^p} \lesssim t^{-\frac d2(1-1/p)-\frac{s+1}{2}},\qquad t\ge2, $$

其中

$$ A_k[M_0]=\frac{\pi}{16\mu}\left(2M_0^kM_0+|M_0|^2e_k\right). $$

故二维次阶剖面为 $-(\log t)\partial_kS_\mu(t)A_k[M_0]$, 且其系数由 $M_0$ 显式确定. 推论 2.7 进一步指出, 在定理 2.6 的条件下,

$$ \big\||\nabla|^s(u(t)-S_\mu(t)M_0)\big\|_{\widehat L^p} =(C_*+o(1))t^{-\frac d2(1-1/p)-\frac{s+1}{2}}\log t, $$

其中 $C_*=\big\||\nabla|^s\partial_kS_\mu(1)A_k[M_0]\big\|_{\widehat L^p}\ne0$. 因而该对数因子是首阶剖面逼近误差中的真实渐近量级, 而非证明估计中的损失.

上述 $s,p$ 还须满足各定理引用的条件 (C). 对更强的线性估计, 论文使用 $p>1$ 时 $s>-d/p'$, $p=1$ 时 $s\ge0$; 这些条件来自低频积分的可积性. 定理 2.4、2.6 的结论不能与其较强的低频假设及维数、$p$ 范围分离理解.

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

证明的主线是把解分解为线性热演化与 Duhamel 非线性项, 分别识别其渐近行为. 论文第 4 节先直接在 Fourier 变量中处理线性项. 由于 $S_\mu(t)M_0$ 的 Fourier 变换是 $e^{-\mu t|\xi|^2}P(\xi)M_0$, 线性余项的核心表达式为

$$ |\xi|^s e^{-\mu t|\xi|^2} \big(\widehat u_0(\xi)-P(\xi)M_0\big). $$

证明将频率区域分为 $|\xi|\le\delta$ 与 $|\xi|>\delta$. 低频区使用 (A) 的任意小误差或 (B) 的 $O(|\xi|^\alpha)$ 控制, 再通过高斯积分获得 $t^{-d/(2p')-(s+\alpha)/2}$ 型衰减; 高频区由 $e^{-\mu t|\xi|^2}$ 得到指数衰减. $p=1$ 时 Fourier–Lebesgue 范数对应频率上的 $L^\infty$ 范数, 因而改用对 $|\xi|^{s+\alpha}e^{-\mu t|\xi|^2}$ 直接取上确界. 这一处理既说明首项由低频极限决定, 也把低频逼近阶 $\alpha$ 转换为时间衰减阶.

非线性部分记为

$$ N[u](t)=\int_0^t e^{\mu(t-\tau)\Delta}P_\sigma\operatorname{div}(u\otimes u)(\tau)\,d\tau. $$

对首阶结果, 论文的关键任务是证明该项相对于 $S_\mu(t)M_0$ 更快衰减. 定理 2.1、2.2 所需的误差范围由命题 5.1 的估计控制:

$$ \||\nabla|^sN[u](t)\|_{\widehat L^p} \lesssim \begin{cases} t^{-\frac d2(1-1/p)-\frac s2-\frac12},&d\ge3,\\ t^{-\frac d2(1-1/p)-\frac s2-\frac12}\log(1+t),&d=2, \end{cases} \qquad t\ge1. $$

该估计使用临界空间小数据所给出的全局 mild 解控制及已有的 $L^p$–$L^1$ 型时间衰减, 再结合 Fourier–Herz 空间的乘积估计处理 $u\otimes u$. 论文明确引用的工具包括: 正则性指标 $s>0$ 时将乘积控制为 $\|f\|_{\widehat L^\infty}\|g\|_{\dot B^s_{p,\sigma}(\widehat L)}$ 及对称项的双线性估计, 以及在相应 $s,p$ 范围内以 $\dot B^{d/p}_{p,\infty}(\widehat L)$ 控制乘积的估计. 这些工具的作用是把二次项转化为可由已知时间衰减控制的因子, 而非仅仅形式上将非线性视为小量.

二阶展开需要更细地提取 Duhamel 项的主导部分. 在 $d\ge3$ 时, 非线性产生的有效振幅由

$$ \int_0^\infty\!\int_{\mathbb R^d}(u_k u)(x,\tau)\,dx\,d\tau $$

给出. 证明必须控制有限时间积分与时间无穷积分之间的误差, 并在 Duhamel 核中分离出带一个空间导数的 Oseen 剖面. 条件 $\alpha>1$ 确保初值的低频余项不会与该次阶项同阶; 二维定理则要求 $\alpha\ge1$.

二维的困难在于非线性时间积分处于临界衰减边界, 使得通常在高维可积的时间贡献累积为对数. 作者不仅估计出 $\log t$ 上界, 还计算了相应的首个对数系数 $A_k[M_0]$. 为处理投影热核在二维的空间尾部, 论文专门证明 Oseen 核加权估计

$$ \int_{\mathbb R^2}|x|\,|S_\mu(x,t)|^2\,dx\lesssim t^{-1/2}. $$

证明利用缩放将问题归结到 $t=1$, 再以二维径向势 $\Phi_\mu=(-\Delta)^{-1}G_\mu$ 表示 Oseen 核, 得到远场 $|S_\mu(y,1)|\lesssim |y|^{-2}$、原点附近有界的估计. 该加权控制服务于二维非线性项的精细渐近计算, 也是显式对数系数推导中的专门工具.

因此, 论文的实质性技术推进不是单纯套用衰减估计, 而是把低频 Fourier 渐近、临界小数据全局控制与 Duhamel 项的时间累积结构衔接起来: 低频分析识别线性首项, 乘积估计控制非线性余项, 再通过维数相关的时间可积性区分 $d\ge3$ 的有限次阶振幅与 $d=2$ 的对数修正. 从提供内容中无法确认第 5 节每一步时间区间分割与所有中间引理的完整细节, 因此不进一步声称某种未在摘录中明确展示的分解公式.

4. 与之前工作的比较

与经典空间可积初值理论相比, 最根本的差别是渐近数据的刻画方式不同. Carpio 的展开及 Fujigaki–Miyakawa 的后续结果以 $L^1$ 或带权 $L^1$ 等空间衰减条件为基础, 其系数由初值空间矩及非线性时间积分决定. 对 $u_0\in L^1$ 的散度自由场, $\int u_0=0$, 因而热核质量项 $G_\mu(t)\int u_0$ 消失; Fujigaki–Miyakawa 的高阶展开能在这种消失后继续提供非平凡描述. 本文不以普通空间积分矩作为首项系数, 而以低频条件 $\widehat u_0(\xi)\sim P(\xi)M_0$ 定义广义低频量, 从而得到 $S_\mu(t)M_0$ 作为非零首项. 这避免了直接要求 $u_0$ 有相应空间矩, 但并不意味着完全没有数据限制: 仍须有 $u_0\in\dot B^0_{1,\infty}(\widehat L)$、小的临界范数以及指定的低频渐近行为.

对一阶渐近, 本文证明的不仅是弱条件下的归一化误差趋于零, 还在更强低频速率假设下给出定量改进. $d\ge3$ 时, 线性低频误差带来 $t^{-\alpha/2}$ 级别, 而非线性可达到的次阶衰减封顶于 $t^{-1/2}$, 因而余项指数出现 $\min\{\alpha,1\}$. $d=2$ 时非线性贡献再出现 $\log(1+t)$. 这些结论覆盖 $\widehat L^p$ 范数及导数阶数 $s$ 的明确范围, 与只在某个物理空间范数陈述的结果有所不同.

高维二阶展开与 Fujigaki–Miyakawa 型公式有共同之处: 非线性积分 $\int_0^\infty\int (u_k u)$ 进入次阶系数, 并由 Oseen 核的导数承载. 但两者的剖面结构并不相同. 经典 $L^1$ 展开在质量项消失后, 初值一阶矩项和非线性贡献可能同阶; 本文则先有 $S_\mu(t)M_0$ 首项, 次阶导数 Oseen 项比首项额外衰减 $t^{-1/2}$. 因此本文呈现的是层级分明的首项—次项结构, 而非将同一衰减阶的多个贡献合并为首阶剖面. 作者还指出 Fujigaki–Miyakawa 型结果在 Fourier–Herz 空间中的相关工作见文献 [29]; 本文在该背景上进一步处理低频广义矩和二阶渐近.

二维结果相较 Carpio 的公式有更明确的差异. Carpio 曾得到带 $\log t$ 的形式展开, 但在其所述散度自由且 $L^1$ 的条件下, $M=\int u_0=0$, 相关的首项及对数项均可能恒等消失, 因而该形式不能提供对应阶的非平凡剖面. 本文的二维对数项则为

$$ -(\log t)\partial_kS_\mu(t)A_k[M_0],\qquad A_k[M_0]=\frac{\pi}{16\mu}(2M_0^kM_0+|M_0|^2e_k), $$

在本文非零 $M_0$ 的设定下给出显式、非平凡的次阶剖面. 推论 2.7 进一步证明逼近首项的误差确实具有 $t^{-\frac d2(1-1/p)-\frac{s+1}{2}}\log t$ 的渐近量级, 说明对数不是估计技术造成的损失. 作者将其解释为二维二次非线性时间可积性的临界效应; 这属于由定理支持的渐近机制说明, 与关于其他抛物方程中类似现象的类比应区分看待.

适用范围和代价方面, 本文结果适用于全空间 $d\ge2$ 的小临界数据全局 mild 解, 不是大数据 Leray–Hopf 弱解的渐近理论. 首阶结果要求 $\dot B^{-1}_{\infty,2}(\widehat L)\cap\dot B^0_{1,\infty}(\widehat L)$, 二阶结果还要求更强的低频逼近阶: 高维 $\alpha>1$, 二维 $\alpha\ge1$. 因而从经典空间矩条件转向低频条件, 是改变并拓宽刻画框架, 不是无条件放宽所有假设. $p,s$ 的范围也受低频可积性、乘积估计和导数损失约束; 高维二阶结论仅陈述于 $1<p\le\infty$, 而二维结论允许 $1\le p\le\infty$. 论文未在所给摘录中定量比较新旧假设之间的严格包含关系, 也未证明本文条件覆盖所有不具空间衰减的初值; 这些更强的包含性判断从提供内容中无法确认.

4. 基于涡量的平面定常 Navier–Stokes 方程的尖锐基本速度估计

Sharp Basic Velocity Estimates for the Plane Steady Navier–Stokes Equations through Vorticity

作者: Mikhail Korobkov, Zhixiang Zhuang

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30844v1

摘要翻译: 我们考察平面区域中一般定常 Navier-Stokes 解的一些新速度估计。具体而言，我们得到了一些最优估计，用环域内涡量的 $L^2$ 范数来控制两个同心圆上速度平均值之差。通过一些出乎意料的反例，我们证明了这些估计是尖锐的。

Comment: 无

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

1. 研究对象和背景

本文研究二维平面区域中的稳态、无外力 Navier–Stokes 方程。区域 $\Omega\subset\mathbb R^2$ 可以无界, 未知量是速度场 $w=(w_1,w_2)$ 与压力 $p$, 方程为

$$ -\Delta w+(w\cdot\nabla)w+\nabla p=0,\qquad \nabla\cdot w=0\quad\text{in }\Omega. $$

作者主要考虑具有有限 Dirichlet 积分的 $D$-解, 即

$$ \int_\Omega |\nabla w|^2\,dx<\infty. $$

在二维, 这一条件本身不能推出速度有界, 甚至不能推出其圆周平均有界。与三维中由 Sobolev 嵌入获得的控制不同, 平面上有限 Dirichlet 能量允许函数以对数尺度增长。此困难与二维障碍物绕流问题及 Stokes 悖论密切相关: 对外区域中的流动, 如何在一般情形下控制速度, 尤其是如何理解无穷远处的行为, 仍是重要问题。本文不直接解决障碍物绕流的存在性或渐近问题, 而是研究一类局部但基础的速度估计, 为理解平面稳态流的速度控制提供更精细的定量信息。

设 $S_r=\{z:|z|=r\}$, 并用 $w(r)$ 表示速度在圆 $S_r$ 上的平均值。考虑同心圆之间的环域

$$ \Omega_{r_1,r_2}=\{z\in\mathbb R^2:r_1<|z|<r_2\},\qquad 0\le r_1<r_2, $$

以及其中的 Dirichlet 能量和涡量能量

$$ D=\int_{\Omega_{r_1,r_2}}|\nabla w|^2,\qquad D_\omega=\int_{\Omega_{r_1,r_2}}|\omega|^2, \qquad \omega=\partial_y w_1-\partial_x w_2. $$

核心问题是估计 $|w(r_2)-w(r_1)|$。对一般 $W^{1,2}$ 函数, 径向积分与 Cauchy–Schwarz 不等式给出

$$ |f(r_2)-f(r_1)|\le \frac{1}{\sqrt{2\pi}} \left(\ln\frac{r_2}{r_1}\right)^{1/2} \left(\int_{\Omega_{r_1,r_2}}|\nabla f|^2\right)^{1/2}. $$

其中的对数因子随半径比增大而增长。Navier–Stokes 方程的结构能显著改善这一一般函数估计。Gilbarg–Weinberger 对压力圆周平均证明了由能量控制的估计; 后续工作又得到稳态速度的基本估计。在环域情形, 该估计具有形式

$$ |w(r_1)-w(r_2)|\le C\sqrt{\ln(2+\mu)}\,D^{1/2}, \qquad \mu=\frac{1}{r_1m},\quad m=\max_{r\in[r_1,r_2]}|w(r)|. $$

而当解定义在整个圆盘内时, 已有结果可去掉这一对数因子, 得到 $|w(r_1)-w(r_2)|\le C D^{1/2}$。这些结果控制的是完整速度梯度能量, 但对散度为零的二维速度场, 圆周平均差还满足涡量表示公式

$$ w(r_2)-w(r_1)=\frac{1}{2\pi}\int_{\Omega_{r_1,r_2}} \omega(z)\frac{z^\perp}{|z|^2}\,dx\,dy, \qquad z^\perp=(-y,x). $$

因此, 圆周平均差直接由涡量决定; 若环域内涡量为零, 该差就为零。仅以 $D$ 控制的既有估计没有显露这一结构, 在涡量相对较小时可能不够精细。本文的目标是把已有基本速度估计改写为显式依赖涡量能量的界, 并检验能否完全去掉完整梯度能量的依赖。作者证明了带有额外对数损失的正向估计, 同时构造反例否定了仅以涡量 $L^2$ 范数控制圆周平均差的自然猜想。换言之, 研究所补的空白既包括对原有基本估计的涡量敏感化, 也包括厘清这种敏感化能够达到的边界。

2. 主要定理或主要结果

记 $D=\int_{\Omega_{r_1,r_2}}|\nabla w|^2$, $D_\omega=\int_{\Omega_{r_1,r_2}}|\omega|^2$, 其中 $w$ 是相应区域上的稳态 Navier–Stokes 解。文中常数 $C$ 为普适常数。

圆盘内的估计 (定理 4). 若 $w$ 是圆盘 $B_R$ 内的解, 则对任意 $0\le r_1<r_2\le R$, 有

$$ |w(r_1)-w(r_2)|\le C\left[\ln\left(2+\frac{D}{D_\omega}\right)\right]^{1/2}D_\omega^{1/2}. $$

当 $D_\omega=0$ 时, 圆周平均差由涡量表示公式直接为零; 上式中的比值可按正能量情形理解。该结果将差值控制在涡量能量尺度上, 但保留了 $D/D_\omega$ 的对数损失。

环域内的估计 (定理 5). 若 $w$ 是环域 $\Omega_{r_1,r_2}$ 内的解, 并记

$$ m=\max_{r\in[r_1,r_2]}|w(r)|,\qquad \mu=\frac{1}{r_1m}, $$

则文中给出的估计为

$$ |w(r_1)-w(r_2)|\le C \left[\ln(2+\mu)\right]^{1/2} \left[\ln\left(2+\frac{D}{D_\omega}\right)\right]^{1/2} D_\omega^{1/2}. $$

相较于圆盘情形, 环域估计多出反映原有环域基本速度估计的因子 $\sqrt{\ln(2+\mu)}$。定理陈述中参数 $\mu$ 使用 $r_1$, 因而在退化端点 $r_1=0$ 的解释需依照原定理的适用约定; 所提供摘录没有进一步说明这一端点约定。

插值型推论 (推论 6). 对任意 $\delta\in(0,1)$, 已有的圆盘和环域基本速度估计仍成立, 若将右端的 $D$ 替换为

$$ D_\omega^\delta D^{1-\delta}. $$

对应常数只依赖于 $\delta$, 以及环域估计中原有的普适常数。此结论给出介于完整梯度能量与涡量能量之间的连续插值控制。

仅由涡量控制的猜想不成立 (定理 8). 作者否定了如下圆盘估计的普适版本:

$$ |w(1)-w(0)|\le C\|\omega\|_{L^2(B_1)}. $$

具体地, 存在一列在 $B_1$ 内满足稳态 Navier–Stokes 方程的解 $w_k$, 其涡量为 $\omega_k$, 且

$$ \frac{|w_k(1)-w_k(0)|}{\|\omega_k\|_{L^2(B_1)}} =\frac{1}{2\pi\|\omega_k\|_{L^2(B_1)}} \left|\int_{B_1}\omega_k(z)\frac{z^\perp}{|z|^2}\,dx\,dy\right| \longrightarrow\infty. $$

因此, 定理 4 中的圆盘圆周平均差不可能在一般解类上仅由涡量的 $L^2$ 范数以线性尺度控制。该反例不否认涡量表示公式, 而是说明其积分核 $|z|^{-1}$ 与允许涡量分布之间可能产生任意大的相对比值。

用于构造反例的存在性结果 (定理 9). 取固定调和函数 $\psi_A\in C^2(\overline{B_1})$, 令 $\widehat w_A=\nabla^\perp\psi_A$, 并给定 $\alpha\in\mathbb R$ 及 $\omega_1\in H^1(B_1)$。若

$$ \|\omega_1\|_{L^2(B_1)}\le C_*,\qquad \widehat w_A\cdot\nabla\omega_1\equiv0\quad\text{in }B_1, $$

则剖面系统 (1.11) 至少存在一个 $\omega\in H^1(B_1)$ 解, 并满足

$$ \|\nabla\omega\|_{L^2(B_1)} \le C\left( \|\omega_1\|_{L^4(B_1)}\|\omega_1\|_{L^2(B_1)} +\|\nabla\omega_1\|_{L^2(B_1)} \right). $$

$C_*,C$ 是普适常数, 特别是不依赖于 $\alpha$、$\psi_A$ 或 $\omega_1$。这一关于 $\alpha$ 一致的先验估计使作者能够令 $\alpha\to\infty$, 进而实现定理 8 的构造。上述结果是论文明确陈述的主要定量结论; 所提供摘录未包含定理 5 与定理 9 的完整证明细节。

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

正向估计的关键不是重新从稳态方程直接推导基本速度估计, 而是把既有基本估计、涡量积分恒等式及解析部分的空间衰减组合起来。困难在于, 涡量恒等式给出的直接估计会产生半径比的对数, 而完整能量估计虽然可借助 Navier–Stokes 结构改善, 却没有显式体现涡量的大小。证明需要在这两种控制之间选择适当的中间半径, 让局部涡量估计与解析部分的衰减恰好平衡。

作者先利用 Cauchy–Pompeiu 表示, 将速度写成两部分 $w=w_A+w_B$: $w_A$ 是复解析部分, $w_B$ 是由涡量生成的部分, 并选择常数使 $w_A$ 在每个同心圆上的平均值为零。相应地, 圆周平均差由涡量决定, 而 $w_A$ 的作用体现在对完整速度梯度能量的分解控制中。对涡量生成部分, 使用平面 Biot–Savart/Cauchy–Pompeiu 表示的能量等距关系

$$ \|\nabla w_B\|_{L^2(\mathbb R^2)}=\|\omega\|_{L^2(\Omega_{r_1,r_2})}. $$

在解析部分上, 作者先证明最大值控制: 当 $0<4r_1<r_2$ 时, 零圆周平均条件、在合适内外圆上的能量选圆以及调和函数最大值原理给出

$$ \sup_{\Omega_{2r_1,r_2/2}}|w_A|\le 4\sqrt{D_A}, \qquad D_A=\int_{\Omega_{r_1,r_2}}|\nabla w_A|^2. $$

再用 Cauchy 积分公式估计 $\nabla w_A$, 得到中间尺度上的分段衰减

$$ |\nabla w_A(z)|\le \begin{cases} C\sqrt{D_A}\,r_1|z|^{-2},&4r_1<|z|<\sqrt{r_1r_2},\\ C\sqrt{D_A}\,r_2^{-1},&\sqrt{r_1r_2}<|z|<r_2/4. \end{cases} $$

这是把解析性转化为可用于径向能量估计的核心步骤: 靠近内边界的衰减是 $|z|^{-2}$ 型, 靠近外侧则由 $r_2^{-1}$ 控制。它使作者可以在不同尺度截断解析部分的贡献, 而不必让整个环域的对数损失直接进入最终界。

定理 4 的证明先处理易于控制的情形。若 $D\le100D_\omega$, 直接使用已知圆盘基本估计 $|w(r_1)-w(r_2)|\le C\sqrt D$。若 $D$ 远大于 $D_\omega$, 且半径比不足以产生目标对数, 则直接对涡量恒等式应用 Cauchy–Schwarz。剩余情形中, 通过边界小区间的涡量估计将问题缩减到较内部的圆周平均差, 并以几何平均半径 $r_0=\sqrt{r_1r_2}$ 分成内外两段。内侧选取

$$ r_*\asymp r_1\sqrt{D/D_\omega}, $$

外侧选取

$$ r_{**}\asymp r_2\sqrt{D_\omega/D}. $$

在这些尺度上, 解析部分能量的贡献被压到 $\sqrt{D_\omega}$ 量级; 剩余的涡量贡献通过径向 Cauchy–Schwarz 产生 $\sqrt{\ln(D/D_\omega)}$。这解释了最终对数损失的来源, 也说明证明并未得到纯粹的 $D_\omega^{1/2}$ 界。

定理 5 除上述两种能量机制外, 还需处理环域基本估计中的速度参数 $\mu$。作者先将问题约化到圆周平均速度在所考虑区间内不发生大幅降至极小值的情形: 若 $|w(r)|$ 从最大值下降到其五分之一, 则选取达到该阈值的子环域, 在该子环域应用估计, 并利用子环域参数不大于原参数的单调性恢复结论。所提供文本在定理 5 证明的这一约化处截断, 因而其后如何逐步结合解析部分衰减和双对数尺度, 无法从摘录中完整确认; 不宜据此补写未展示的论证。

反例的技术障碍更本质: 需要构造满足无外力稳态 Navier–Stokes 方程的精确解, 不能任意指定一个具有所需积分性质的涡量。作者通过涡量剖面系统解决这一问题。令

$$ w=\alpha\widehat w_A+w_B, \qquad \widehat w_A=\nabla^\perp\psi_A, $$

其中 $\widehat w_A$ 是无旋、无散的调和速度, $w_B$ 由涡量经 Biot–Savart 型积分给出。复变量表示下有

$$ \partial_{\bar z}w=\partial_{\bar z}w_B=\frac{i\omega}{2}. $$

选取边界涡量 $\omega_1$ 沿 $\psi_A$ 的等值线为常数, 即 $\widehat w_A\cdot\nabla\omega_1=0$, 可使放大项 $\alpha\widehat w_A\cdot\nabla\omega$ 与参考剖面相容。对所得椭圆输运方程, 作者采用 Leray–Schauder 不动点定理, 沿经典稳态方程存在性论证建立解, 并得到不依赖于 $\alpha$ 的 $H^1$ 先验估计。这种对大参数一致的界是构造反例的关键创新点: 没有它, 将调和背景流速度放大至无穷大时, 无法保证剖面解保持受控。

具体构造时取 $\psi_A=xy$, 选取支撑于 $B_1\cap\{xy>t_0\}$ 的光滑初始涡量 $\omega_1$, 使

$$ \left|\int_{B_1}\omega_1(z)\frac{z^\perp}{|z|^2}\,dx\,dy\right| \gg \|\omega_1\|_{L^2(B_1)}. $$

将其乘以足够小的常数, 以满足定理 9 的小 $L^2$ 假设, 同时保持积分与范数之比的巨大性。随后令 $\alpha_k\to\infty$, 由一致先验估计取得相应解 $\omega_k$, 并证明 $\omega_k\to\omega_1$ 于 $L^2(B_1)$。于是涡量积分所产生的大比值传递到真实 Navier–Stokes 解, 得到定理 8。故本文的反例不是单纯的函数空间反例, 而是通过一个保持参数一致性的存在性定理, 将所需涡量分布嵌入真实方程解中。

4. 与之前工作的比较

本文首先以两类既有估计为基准。第一类是一般 Sobolev 函数的径向估计, 其右端含 $\sqrt{\ln(r_2/r_1)}$; 该对数来自单纯的径向 Cauchy–Schwarz, 并不能利用方程结构消去。第二类是 Gilbarg–Weinberger 关于压力圆周平均的结果, 其以环域 Dirichlet 能量控制压力平均差。本文讨论的是速度圆周平均, 并利用二维不可压缩场的涡量结构改进已有速度基本估计, 研究对象和控制机制均不同于压力估计。

针对稳态 Navier–Stokes 速度, 文中引用此前的基本估计 (定理 2), 其环域版本为

$$ |w(r_1)-w(r_2)|\le C\sqrt{\ln(2+\mu)}\,\sqrt D, \qquad \mu=(r_1m)^{-1},\quad m=\max_{r\in[r_1,r_2]}|w(r)|. $$

作者指出该估计在外区域半径比增大时, 不会像一般 Sobolev 估计那样仅因 $r_2/r_1$ 增大就产生发散系数; 它利用了方程结构, 且对于小来流速度的障碍物绕流, 文中引用先前结果说明其对数因子在一般意义下不能删去。圆盘中的既有结果 (定理 3) 则有更强结论 $|w(r_1)-w(r_2)|\le C\sqrt D$, 不含该对数因子。本文定理 4、5 的推进不是取代这些估计, 而是在其基础上把完整梯度能量 $D$ 转换成涡量能量 $D_\omega$, 并分别保留圆盘与环域原有的尺度结构。

这一推进有明确代价。定理 4 并未得到猜想中的纯粹 $C\sqrt{D_\omega}$, 而是多出

$$ \sqrt{\ln(2+D/D_\omega)}. $$

环域定理 5 还保留原有的 $\sqrt{\ln(2+\mu)}$, 因而两种损失分别反映涡量相对于完整梯度能量的比例以及环域速度尺度。推论 6 进一步给出 $D_\omega^\delta D^{1-\delta}$ 型控制, 但常数依赖于 $\delta$; 当 $\delta$ 趋近端点时, 文中结论不保证常数仍一致。因而本文建立的是涡量主导、但仍依赖全梯度能量的近似最优控制, 不是只用涡量能量的闭合估计。

定理 8 对“只用涡量”的可能改进给出强否定答案。尽管圆周平均差在恒等式上完全由涡量决定, 该线性泛函在一般解类中并不由 $\|\omega\|_2$ 一致控制: 在圆盘上甚至存在真实稳态解序列, 使平均差与涡量 $L^2$ 范数之比趋于无穷。这一结果比指出某个抽象积分算子在函数空间上无界更强, 因为反例满足原非线性方程。另一方面, 反例只展示不存在普适线性估计, 不直接给出该比值的具体增长速率, 也不否定带有 $D/D_\omega$ 对数损失的定理 4。文中将其表述为定理 4 尺度不可进一步简化的重要证据, 但没有证明定理 4 的全部常数或对数增长在所有参数区间上都是最优的。

反例依赖定理 9 的特殊构造, 其假设包括初始涡量的小 $L^2$ 范数以及沿调和流线不变。它并非对任意给定涡量数据的普遍实现定理, 也没有声称适用于任意区域或边界条件。构造发生在 $B_1$, 并利用可任意放大的调和背景速度参数 $\alpha$; 关键代价是采用专门设计的边界剖面及相容条件。定理 9 的先验界对 $\alpha$ 一致, 但所提供内容没有给出解的唯一性, 因而不能将该构造理解为对任意初值都有稳定或唯一响应。

文中明确提及的相关工作包括 Gilbarg–Weinberger 的压力估计及其二维流动分析, 此前关于环域速度基本估计的工作, 圆盘内无对数因子的估计, 以及小来流障碍物绕流中对数损失不可删的结果; 作者还引用相关综述和经典稳态方程存在性理论。所提供 PDF 摘录没有列出参考文献表的完整书目信息, 因此除文中所述作者、编号与结论外, 无法可靠补充各编号文献的完整题名或对所有先前结果作逐条比较。总体而言, 本文的主要推进是将既有基本速度估计改造成涡量敏感的带对数估计, 并以满足方程的反例证明纯涡量线性控制失败; 其适用范围仍是二维无外力稳态解及文中具体圆盘、环域设定, 对非定常问题、带外力问题及更一般边界条件的推广, 从提供内容中无法确认。

5. Boussinesq 方程的 Korteweg-de Vries 近似的长时间有效性

Long-time validity of the Korteweg-de Vries approximation for the Boussinesq equation

作者: Younghun Hong, Junyeong Jang

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30786v1

摘要翻译: 我们研究良性 Boussinesq 方程的 Korteweg-de Vries (KdV) 近似。在 Hong 和 Yang \cite{HY2024} 建立的低正则局部时间论证基础上，我们证明该近似在对数长时间区间上仍然有效。关键在于一个正则性保持论证，它利用 Boussinesq 方程的缩放后守恒律对局部近似进行迭代。

Comment: 24 pages

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

1. 研究对象和背景

本文讨论良性Boussinesq方程在长波、小振幅尺度下由两支反向传播的KdV波描述的严格性，以及这种近似能维持多长时间。原方程为

$$ \partial_t^2 U=\partial_x^2U-\partial_x^4U-\partial_x^2(U^2),\qquad (t,x)\in\mathbb R\times\mathbb R. $$

在长波尺度下，令 $U(t,x)=\epsilon^2u_\epsilon(\epsilon^3t,\epsilon x)$，其中 $0<\epsilon\ll1$。相应的缩放方程是

$$ \epsilon^4\partial_t^2u_\epsilon=\partial_x^2(1-\epsilon^2\partial_x^2)u_\epsilon-\epsilon^2\partial_x^2(u_\epsilon^2). $$

这一极限并非单向传播：线性Boussinesq色散关系包含向左和向右传播的两支。将解分解为两种移动坐标中的剖面后，形式上有

$$ u_\epsilon(t,x)=u_\epsilon^+(t,x-t/\epsilon^2)+u_\epsilon^-(t,x+t/\epsilon^2), $$

而两支剖面在 $\epsilon\to0$ 时分别趋向KdV解。本文使用的KdV方程为

$$ 2\partial_tw^\pm\mp\partial_x^3w^\pm\mp2w^\pm\partial_xw^\pm=0. $$

相应地，缩放后的耦合剖面系统可写为

$$ \partial_tu_\epsilon^\pm=\pm\frac{\partial_x^3}{\langle\epsilon\partial_x\rangle+1}u_\epsilon^\pm\pm\frac12\frac{\partial_x}{\langle\epsilon\partial_x\rangle}(u_\epsilon^2)\left(t,x\pm\frac{t}{\epsilon^2}\right), $$

其中 $u_\epsilon$ 是两支剖面在相应移动坐标中的合成，$\langle\epsilon\partial_x\rangle$ 的傅里叶符号为 $\langle\epsilon\xi\rangle=(1+\epsilon^2\xi^2)^{1/2}$。问题的核心是：在初值正则性较低时，如何以定量误差证明这个形式极限，并且不只在固定的局部时间内有效，而能随 $\epsilon\to0$ 延长有效时间。

这一问题处在调制近似与低正则色散方程理论的交界处。KdV方程是许多长波模型的有效方程，严格控制原模型解与KdV解的差异，既能说明渐近模型在何种数据类别中可信，也能为从复杂系统提取主导动力学提供误差依据。过去的近似论证常要求较高的Sobolev正则性，有时还要求加权空间条件，并且结论多限于局部时间。Hong–Yang此前针对Boussinesq–KdV问题建立了仅略高于 $L^2$ 的局部时间论证，并去除了加权范数假设；但局部结果能否反复应用，取决于控制解的守恒量是否能控制局部理论所需的范数。

这里恰好存在一个尺度上的障碍。原方程守恒能量

$$ E_0[U]=\frac12\int_{\mathbb R}\left(U^2+(\partial_xU)^2+(\partial_x^{-1}\partial_tU)^2\right)\,dx-\frac13\int_{\mathbb R}U^3\,dx. $$

缩放后对应的守恒量为

$$ E_\epsilon[u_\epsilon]=\frac12\int_{\mathbb R}\left(u_\epsilon^2+\epsilon^2(\partial_xu_\epsilon)^2+\epsilon^4(\partial_x^{-1}\partial_tu_\epsilon)^2\right)\,dx-\frac{\epsilon^2}{3}\int_{\mathbb R}u_\epsilon^3\,dx. $$

当 $\epsilon$ 很小时，导数项带有小系数，因此该能量主要提供一致的低阶、特别是 $L^2$ 级控制，并不能直接给出一致的 $H^1$ 控制。若局部近似的存在时间依赖高阶 $H^s$ 范数，那么仅凭守恒量就无法保证每次重启局部论证时都有统一时间步长。本文试图补上的空白，是在低正则局部近似框架中建立高阶正则性的持续性，使局部寿命只依赖守恒量控制的低范数，而高阶范数只影响局部误差大小；在此基础上利用缩放守恒律迭代局部近似，得到对数长时间有效性。论文还将该机制作为更一般的分析思路提出，但所给文本中的严格定量结论针对Boussinesq–KdV模型。

2. 主要定理或主要结果

主要结论是定理1.1。设 $0<s\le5$，并假设初始剖面满足一致的加权Sobolev有界性

$$ \sup_{0<\epsilon\le1}\sum_{\pm}\|\langle\epsilon\partial_x\rangle u_{\epsilon,0}^\pm\|_{H^s_x}<\infty, \qquad \sum_{\pm}\|w_0^\pm\|_{H^s_x}<\infty. $$

令 $u_\epsilon^\pm\in C(\mathbb R;H^s_x)$ 为耦合Boussinesq剖面系统(1.5)的解，令 $w^\pm\in C(\mathbb R;H^s_x)$ 为KdV方程(1.6)的解，初值分别为 $u_{\epsilon,0}^\pm$ 和 $w_0^\pm$。则存在 $\epsilon_0\in(0,1]$ 及常数 $K_0,K_s>0$，使得对 $0<\epsilon\le\epsilon_0$ 和所有 $t\in\mathbb R$，有

$$ \|u_\epsilon^\pm(t)-w^\pm(t)\|_{L^2_x} \le e^{K_0|t|}\left(K_0\|u_{\epsilon,0}^\pm-w_0^\pm\|_{L^2_x}+K_s\epsilon^{\min\{2s/5,\,1/2\}}\right). $$

常数的依赖具有迭代论证所需的区分：$K_0$ 与 $\epsilon,t$ 无关，且其指数增长率只由低阶界

$$ \sup_{0<\epsilon\le1}\sum_{\pm}\|\langle\epsilon\partial_x\rangle u_{\epsilon,0}^\pm\|_{L^2_x}, \qquad \sum_{\pm}\|w_0^\pm\|_{L^2_x} $$

决定；$K_s$ 可进一步依赖假设中的高阶 $H^s$ 界。因而高阶正则性影响误差常数，但不决定每次局部迭代可采用的时间步长。误差的 $\epsilon$ 次幂随 $s$ 增加，在 $s=5/4$ 处由 $2s/5$ 达到 $1/2$，之后饱和于 $1/2$；文中给出的结论不宣称更高正则性会继续改善这一幂次。

定理直接给出对数时间有效性。若初始误差为零，或更一般地满足

$$ \|u_{\epsilon,0}^\pm-w_0^\pm\|_{L^2_x}\lesssim\epsilon^{\min\{2s/5,\,1/2\}}, $$

则对任意

$$ 0<c<\frac{1}{K_0}\min\left\{\frac{2s}{5},\frac12\right\}, $$

有

$$ \lim_{\epsilon\to0}\sup_{|t|\le c\log(1/\epsilon)} \|u_\epsilon^\pm(t)-w^\pm(t)\|_{L^2_x}=0. $$

这里的时间是缩放时间；结论是误差在该对数增长区间内趋于零，而不是在任意更长时间上趋于零。更一般地，误差估计本身给出允许时间尺度约为

$$ |t|\ll \frac1{K_0}\log\frac1{K_0\|u_{\epsilon,0}^\pm-w_0^\pm\|_{L^2_x}+K_s\epsilon^{\min\{2s/5,1/2\}}}. $$

作者还将结论换回原始变量。若 $U_\epsilon$ 是原Boussinesq方程的解，初值由两支剖面按论文所列公式重构，且 $W_\epsilon^\pm$ 是以相同剖面初值为初值的KdV解，则对 $T>0$ 有

$$ \sup_{|t|\le T/\epsilon^3} \left\|U_\epsilon(t,x)-\epsilon^2\sum_{\pm}W_\epsilon^\pm\bigl(\epsilon^3t,\epsilon(x\mp t)\bigr)\right\|_{L^2_x} \le K_s\epsilon^{3/2+\min\{2s/5,1/2\}}e^{K_0T}. $$

该式中的 $\epsilon^{3/2}$ 来自空间变量缩放下的 $L^2$ 范数因子，$T=\epsilon^3|t|$ 是缩放时间。因而对数缩放时间 $T\lesssim\log(1/\epsilon)$ 对应原变量中 $|t|\lesssim\epsilon^{-3}\log(1/\epsilon)$ 的时间尺度；在上述初值匹配条件下，指数因子与误差幂次共同决定具体可取的对数常数。定理适用于 $0<s\le5$，但所给摘录未列出各中间局部估计的完整形式，因而无法据此补写其更精细的依赖或端点之外的结论。

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

证明需要克服的主要障碍不是单纯得到局部误差，而是使局部近似可以用一个由守恒量控制的统一时间步长反复重启。缩放能量只给低阶控制，若沿用高正则局部理论而使寿命依赖 $H^s$ 范数，高阶范数在迭代中可能增长，局部时间步长也随之缩短，无法直接获得所需的长时间结论。作者的核心处理是把固定点构造改造成低范数度量下的高正则有界集：高阶正则性在解的构造过程中保留，但不进入收缩度量和局部寿命的决定因素。

第一步是将二阶方程改写成适合色散分析的两支一阶系统。作者对线性部分对角化，分离向左、向右传播的模式，并在各自移动坐标中消除主导平移。剖面系统的线性色散相位为

$$ s_\epsilon(\xi)=\frac{\xi}{\epsilon^2}(\langle\epsilon\xi\rangle-1)=\frac{\xi^3}{\langle\epsilon\xi\rangle+1}, $$

当 $\epsilon\to0$ 时趋于Airy相位 $s(\xi)=\xi^3/2$。相应线性群记作 $S_\epsilon^\pm(t)$，KdV线性群记作 $S^\pm(t)$。通过Duhamel公式，作者使用耦合积分方程

$$ u_\epsilon^\pm(t)=S_\epsilon^\pm(t)u_{\epsilon,0}^\pm\pm\frac12\int_0^t S_\epsilon^\pm(t-t_1)\frac{\partial_x}{\langle\epsilon\partial_x\rangle} \left(u_\epsilon^\pm(t_1)+e^{\pm 2t_1\partial_x/\epsilon^2}u_\epsilon^\mp(t_1)\right)^2dt_1. $$

这种积分形式避免直接在微分方程层面对不同线性算子作高阶导数比较，也显式呈现了两支间带快速平移的耦合项；它是低正则估计和极限比较的起点。

第二步是在与色散相位匹配的Bourgain空间中建立局部理论。空间 $X^{s,b}_{\epsilon,\pm}$ 的时空傅里叶权重为 $\langle\xi\rangle^s\langle\tau\pm s_\epsilon(\xi)\rangle^b$，KdV空间则使用Airy相位。证明利用 $b>1/2$ 时的时间连续性嵌入、线性流估计、时间截断估计及Duhamel非齐次项估计，将二次非线性映射控制在相应时空范数中。由于两类解对应不同色散相位，不能不加处理地在同一个Bourgain范数中比较；论文指出在近似估计中采用频率局部化的辅助方程处理这一问题，并沿用Hong–Yang的傅里叶分析框架。不过，所给摘录没有展示该辅助方程的具体定义、频率阈值或相位误差估计，故这些细节从提供内容中无法确认。

第三步是持续正则性的关键构造。作者把固定点空间取为同时满足高阶与低阶界的集合，例如

$$ X=\{u\in X^{s,b}_{\epsilon,\pm}:\|u\|_{X^{s,b}_{\epsilon,\pm}}\le M_s,\ \|u\|_{X^{0,b}_{\epsilon,\pm}}\le M_0\}, $$

但以较弱的距离 $d(u,v)=\|u-v\|_{X^{0,b}_{\epsilon,\pm}}$ 赋予度量。作者证明这种集合在弱度量下完备，于是可在低正则范数中实施压缩映射，同时把高阶范数作为不变集约束传播。完备性并非形式细节：仅有高范数有界的集合相对于较弱范数的极限，必须仍属于高正则集合，才能保证固定点既存在于低阶空间又保有所需的 $H^s$ 正则性。论文借助弱收敛与下半连续性所给出的泛函分析判据建立这一点。

为使非线性估计不把高阶导数分散到多个因子上，作者进一步精化双线性估计，将导数集中在单一因子。所用频率基本不等式是

$$ \left\langle\sum_{j=1}^m\xi_j\right\rangle^s\le C_{s,m}\sum_{\ell=1}^m\langle\xi_\ell\rangle^s,\qquad s\ge0. $$

在多重卷积估计中，这一分配方式使高阶权重可由一个输入承担，其余因子只需低阶控制。由此，局部存在时间可以仅依赖 $L^2$ 级界，而高阶 $H^s$ 量在同一时间区间内持续传播；论文进一步证明高阶范数增长至多为时间的指数增长。摘录将这一结论标为命题3.7，但未提供命题的具体不等式常数与证明细节。

最后，作者把精化后的局部近似按统一时间步长迭代。缩放守恒律提供所需的低阶一致控制，故每次可使用相同量级的局部寿命；持续正则性允许每次重启仍处于局部理论的高正则类中。局部误差经反复估计后累积成定理中的指数因子 $e^{K_0|t|}$，而局部渐近误差保留为 $\epsilon^{\min\{2s/5,1/2\}}$。所以指数误差估计与对数时间结论不是额外的独立机制，而是低范数控制局部寿命、高范数只进入误差、再结合守恒律迭代的结果。

4. 与之前工作的比较

与Schneider关于Boussinesq–KdV近似的早期严格论证相比，本文所处的研究脉络强调降低正则性并延长有效时间。论文明确指出Schneider的结果建立在高正则性假设下；Hong–Yang随后利用傅里叶分析和色散空间技术，将局部时间论证推进到仅略高于 $L^2$ 的正则性，并移除了加权范数假设。本文以Hong–Yang的低正则局部结果为基础，主要新增内容不是重新建立一套不同的渐近模型，而是通过持续正则性改造局部理论，使其可在守恒律仅提供低阶控制的情况下迭代。因此，正则性层面的进展应归于前一局部理论，本文的主要推进在时间区间与迭代机制。

与Hong–Yang的局部有效性结论相比，本文给出对任意缩放时间的指数型误差上界，并由此推出在初值误差与渐近误差同阶时，误差趋零可至少延续到 $c\log(1/\epsilon)$ 的缩放时间。其具体上界为

$$ \|u_\epsilon^\pm(t)-w^\pm(t)\|_{L^2_x} \lesssim e^{K_0|t|}\left(\|u_{\epsilon,0}^\pm-w_0^\pm\|_{L^2_x}+\epsilon^{\min\{2s/5,1/2\}}\right), $$

其中常数写法按定理可吸收固定因子。与局部结果相比，这是有效时间从固定的局部区间推进到随小参数增长的对数尺度；但它并不意味着误差在所有长时间都小。指数因子限制了可用时间：当 $|t|$ 达到超过相应误差幂次所允许的对数阈值时，该估计本身不再保证误差趋零。

与BBM近似的情形相比，差异在于守恒量的强弱与局部理论所需正则性的匹配。论文指出，BBM方程的守恒量能给出一致 $H^1$ 控制，而将BBM近似的严格论证正则性要求降至 $H^1$ 后，守恒量足以直接支持局部迭代，已有工作[15]因而得到指数时间误差估计。Boussinesq缩放能量则为

$$ E_\epsilon[u_\epsilon]=\frac12\int\left(u_\epsilon^2+\epsilon^2|\partial_xu_\epsilon|^2+\epsilon^4|\partial_x^{-1}\partial_tu_\epsilon|^2\right)dx-\frac{\epsilon^2}{3}\int u_\epsilon^3dx, $$

导数控制带有 $\epsilon$ 小因子，无法直接给出一致 $H^1$ 界。因此本文不能简单照搬BBM的直接迭代，而必须让局部寿命只依赖低阶守恒控制，同时另行传播高阶正则性。这说明方法的针对性在于处理低阶守恒量与高阶近似误差之间的错配。

从更广泛的文献背景看，作者将方法置于Kenig–Ponce–Vega傅里叶分析、Bourgain空间以及低正则调制近似的发展脉络中，并提及Hong–Kwak–Yang关于Fermi–Pasta–Ulam系统导出KdV的工作，以及Schneider–Wayne的早期高正则结果。论文也指出，长波极限下守恒律退化为低阶控制的困难并非Boussinesq独有，Fermi–Pasta–Ulam问题中亦有类似现象。本文的贡献因此包含一个可迁移的证明组织方式：把低正则固定点度量与高正则有界集结合，再利用守恒量反复启动局部近似。作者称相关持续正则性策略也用于其他问题[13,14]，但从提供内容无法确认这些工作的具体假设和结论，故不宜据此作更细的技术比较。

结论的适用范围和代价也应明确。本文处理一维实值良性Boussinesq方程、其两支耦合剖面系统与KdV极限，数据范围为 $0<s\le5$，误差在 $L^2_x$ 中估计，且初始剖面需要一致的 $\langle\epsilon\partial_x\rangle H^s$ 控制。虽然相对于高正则理论及加权数据要求有明显放宽，但并非无正则性假设；结论也没有给出更高维、其他非线性结构或不同方程的一般定理。误差率在 $s\ge5/4$ 后饱和于 $\epsilon^{1/2}$，显示提高正则性并未在该结果中带来无限改善。指数时间误差的增长率依赖低阶范数，而误差前的常数仍可依赖高阶界；因此这并不是无条件的统一长时间稳定性结论。

就所给PDF摘录而言，作者明确比较了Schneider的高正则局部论证、Hong–Yang的低正则局部论证及BBM的守恒律迭代优势，但没有在摘录中给出这些前人定理的完整数值误差率、精确正则阈值或逐项假设表，也没有展开与所有调制近似文献的系统比较。这些未提供的细节从提供内容中无法确认。

6. $C^1$ Yudovich 类中定常与匀速旋转平面欧拉流的刚性

Rigidity of stationary and uniformly rotating planar Euler flows in the $C^1$ Yudovich class

作者: Boquan Fan, Yuchen Wang, Chunjing Xie, Weicheng Zhan

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30664v1

摘要翻译: 我们证明了二维不可压缩欧拉方程定常解与匀速旋转解的一个刚性定理，其涡量满足

Comment: 无

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

1. 研究对象和背景

本文研究平面不可压缩欧拉方程的定常涡量与匀速旋转涡量剖面。以涡量 $\omega(t,x)$ 为未知量，速度由 Biot-Savart 定律从涡量恢复，方程为

$$ \partial_t\omega+v\cdot\nabla\omega=0,\qquad v=K*\omega,\qquad K(x)=\frac{1}{2\pi}\frac{(-x_2,x_1)}{|x|^2},\quad x\in\mathbb R^2. $$

匀速旋转解是指存在角速度 $\Omega\in\mathbb R$，使得

$$ \omega(t,x)=\omega_0(R_{-\Omega t}x), $$

其中 $R_\theta$ 表示平面旋转。该形式代入欧拉方程后，旋转剖面 $\omega_0$ 必须满足

$$ \bigl(K*\omega_0+\Omega x^\perp\bigr)\cdot\nabla\omega_0=0, $$

本文采用 $x^\perp=(x_2,-x_1)$ 的约定。定常情形对应 $\Omega=0$，方程退化为 $(K*\omega_0)\cdot\nabla\omega_0=0$。研究的关键问题是：在什么条件下，上述输运平衡迫使涡量剖面具有径向对称性？径向涡量显然满足任意角速度下的剖面方程，因此刚性定理要排除的是满足一定结构条件的非径向解。

本文考虑的函数类为

$$ \omega_0\in C^1(\mathbb R^2)\cap L^1(\mathbb R^2)\cap L^\infty(\mathbb R^2). $$

这是 Yudovich 理论中自然的涡量可解性框架之一：$L^1\cap L^\infty$ 控制总涡量与涡量幅度，而额外的 $C^1$ 正则性使作者能够使用经典流函数、正则等值线及其几何性质。本文不是要证明时间演化的适定性，而是分类满足稳态或旋转剖面方程的经典剖面。特别重要的是，$L^1$ 条件并不要求紧支撑，因而允许涡量在空间无穷远处具有非紧支撑尾部。

若 $\psi$ 是标准流函数，满足 $-\Delta\psi=\omega_0$，则相对流函数可取为

$$ \Psi=\psi+\frac{\Omega}{2}|x|^2, \qquad -\Delta\Psi=\omega_0-2\Omega. $$

剖面输运方程意味着涡量沿相对流函数的等值线保持不变。右侧的有效源 $\omega_0-2\Omega$ 若具有固定符号，便为椭圆方程的水平集几何提供了单侧结构。本文的角速度条件正是

$$ \Omega\leq\frac12\inf_{\mathbb R^2}\omega_0 \quad\text{或}\quad \Omega\geq\frac12\sup_{\mathbb R^2}\omega_0, $$

即 $\omega_0-2\Omega$ 不变号。固定符号的有效源使适当符号归一化后的相对流函数成为超调和函数，并且在无穷远趋于 $-\infty$；因此其超水平集有界，正则水平集可按有限个 Jordan 域及其边界曲线分析。

已有刚性理论的重要一类对象是涡旋斑块和紧支撑光滑涡量。Fraenkel 对定常单连通涡旋斑块证明径向对称；后续工作对旋转涡旋斑块的角速度区间及凸性等条件作了推进。Gómez-Serrano、Park、Shi 和 Yao 对旋转涡旋斑块在相应角速度范围内证明了径向刚性，并得到阈值尖锐性；Fan、Wang 和 Zhan 则对紧支撑 $C^2$ 旋转涡量证明了与本文相似的角速度阈值型结果，允许涡量变号。此类结果依赖紧支撑假设的原因之一，是可以直接利用支撑边界或有限区域的几何信息。对于非紧支撑尾部，单靠边界论证并不能直接完成分类。

本文试图填补的空白，正是把角速度阈值下的径向刚性从紧支撑情形推进到 $C^1\cap L^1\cap L^\infty$ 的 Yudovich 型全空间框架，同时涵盖定常情形。难点在于输运方程只保证涡量在每条连通的正则流函数等值曲线上为常数；同一水平值的不同连通分支上，这些常数起初可以不同。因此，不能预先把问题化为单一的半线性方程 $-\Delta\Phi=f(\Phi)$。论文的策略是先保留各个水平集分支的信息，利用分支上的 Bernoulli 关系和全局无穷远分析证明分支最终不能分裂，再导出半线性关系并调用径向对称性定理。

2. 主要定理或主要结果

主要结论是定理 1.1。设 $\Omega\in\mathbb R$，且

$$ \omega_0\in C^1(\mathbb R^2)\cap L^1(\mathbb R^2)\cap L^\infty(\mathbb R^2) $$

满足旋转剖面方程

$$ \bigl(K*\omega_0+\Omega x^\perp\bigr)\cdot\nabla\omega_0=0 \quad\text{in }\mathbb R^2, $$

以及单侧角速度条件

$$ \Omega\leq\frac12\inf_{\mathbb R^2}\omega_0 \quad\text{或}\quad \Omega\geq\frac12\sup_{\mathbb R^2}\omega_0. $$

则存在 $x_0\in\mathbb R^2$ 和函数 $W:[0,\infty)\to\mathbb R$，使得

$$ \omega_0(x)=W(|x-x_0|)\qquad\text{对所有 }x\in\mathbb R^2. $$

若 $\Omega\neq0$，对称中心必须是原点，即可取 $x_0=0$。定理不要求 $\omega_0$ 紧支撑，也不要求其本身非负；真正的符号假设施加在有效源 $\omega_0-2\Omega$ 上。特别地，当 $\Omega=0$ 且定常涡量为单号时，条件自动成立，因此每个 $C^1\cap L^1\cap L^\infty$ 类单号定常涡量都关于某一点径向对称。定理本身陈述的是径向对称，摘要及证明纲要还说明所得流函数的相关超水平集具有径向单调结构；不应将此误写为涡量 $W$ 必然单调，因为所给定理没有作出这一结论。

证明中的核心中间结论可概括如下。对非零剖面选取符号 $\sigma\in\{-1,1\}$，定义

$$ \Phi=\sigma\Psi, \qquad \rho=\sigma(\omega_0-2\Omega)\geq0. $$

通过对流函数作适当归一化，得到

$$ -\Delta\Phi=\rho,\qquad \nabla^\perp\Phi\cdot\nabla\rho=0, \qquad \Phi(x)\to-\infty\quad (|x|\to\infty), $$

其中 $\Phi\in C^2(\mathbb R^2)$，$\rho\in C^1(\mathbb R^2)\cap L^\infty(\mathbb R^2)$。作者对每个正则值 $t<\max\Phi$ 考察超水平集 $E_\Phi(t)=\{\Phi>t\}$ 的连通分支 $D$。每个分支的边界是 Jordan 曲线，且 $\rho$ 在该边界上为常数 $c_D(t)$。因为 $\rho\nabla\Phi$ 无旋，存在 Bernoulli 函数 $H$ 满足

$$ \nabla H=\rho\nabla\Phi. $$

相应的分支公式为

$$ I(t):=\sum_{D\in\mathcal C_\Phi(t)}\int_D\bigl(H(x)-h_D(t)\bigr)\,dx =\int_t^{\max\Phi}A(s)\,ds, $$

其中 $h_D(t)=H|_{\partial D}$，而正则值处

$$ A(t)=\sum_{D\in\mathcal C_\Phi(t)}c_D(t)|D|. $$

另定义分支质量平方和

$$ S(t)=\sum_{D\in\mathcal C_\Phi(t)}\left(\int_D\rho(x)\,dx\right)^2. $$

利用散度定理、余面积公式、Cauchy-Schwarz 不等式及平面等周不等式，作者得到测度分解

$$ -dS=8\pi A(t)\,dt+\mu, $$

其中 $\mu$ 是非负 Radon 测度。它既记录正则水平集上的等周不等式等几何估计的亏损，也记录临界水平处分支分裂造成的质量平方损失。对正则值有

$$ \delta(t):=S(t)-8\pi I(t)=\mu((t,\max\Phi))\geq0. $$

因此，只要证明存在正则值序列 $t_k\to-\infty$ 使 $\delta(t_k)\to0$，单调性与非负性便推出 $\mu=0$。这排除了正则超水平集的分裂，并使作者能够得到几乎处处的半线性关系 $\rho=f(\Phi)$，其中 $f$ 是有界非负 Borel 函数。再结合等周不等式取等，适用 P.-L. Lions 的径向对称性定理，推出径向对称。

对于 $\Omega=0$，令 $\Gamma=\int_{\mathbb R^2}\rho\,dx$。在 Bernoulli 函数于无穷远归一化后，质量平方分解给出

$$ \Gamma^2=8\pi\int_{\mathbb R^2}H\,dx+\mu(( -\infty,\max\Phi)), $$

而全局 Pohozaev 恒等式给出 $\Gamma^2=8\pi\int_{\mathbb R^2}H\,dx$，故 $\mu=0$。在 $\Omega\neq0$ 时，归一化后的有效源含正的常数背景，因而其全空间总质量为无穷，不能照搬有限总质量论证。作者改为在低水平集上估计 $\delta(t)$：旋转二次项使足够低的等值线在大圆环上成为径向图，边界 Pohozaev 恒等式把亏损控制在流函数的角向能量内；沿适当圆环选取水平值并利用角向能量衰减，得到 $\delta(t_k)\to0$，继而仍有 $\mu=0$。此外，渐近关系 $\nabla\Phi(x)=-\frac b2x+o(1)$（$b=2|\Omega|$）把对称中心固定在原点。

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

证明的首要障碍是不能一开始把输运约束改写成全局单值的半线性椭圆方程。由 $\nabla^\perp\Phi\cdot\nabla\rho=0$ 只能推出 $\rho$ 在每条连通的正则水平曲线上为常数；同一个水平值的不同分支，涡量常数可能不同。临界点处水平集还会发生拓扑改变。因此，直接令 $\rho=f(\Phi)$ 会预先假定需要证明的连通性与全局一致性。作者的关键组织方式是先构造一个按连通分支记账的几何框架，再证明分支分裂的总缺陷为零。

首先，作者把旋转剖面方程转成具有固定符号源的椭圆—输运系统。设相对流函数 $\Psi=\psi+\frac{\Omega}{2}|x|^2$，则 $-\Delta\Psi=\omega_0-2\Omega$。依据角速度条件选取 $\sigma$，令 $\Phi=\sigma\Psi$、$\rho=\sigma(\omega_0-2\Omega)$，从而 $-\Delta\Phi=\rho\geq0$。同时，原输运关系给出 $\nabla^\perp\Phi\cdot\nabla\rho=0$。对流函数的对数势作适当归一化后，$\Phi(x)\to-\infty$，于是所有超水平集有界。固定符号的超调和结构还排除了正则超水平分支中的孔洞：正则等值集由有限条 Jordan 曲线组成，每个超水平分支由一条曲线围成。这些事实使得后续的质量和等周估计可以逐分支进行。

其次，利用输运方程导出 Bernoulli 结构。因为

$$ \operatorname{curl}(\rho\nabla\Phi)=\nabla^\perp\Phi\cdot\nabla\rho=0, $$

平面全空间上存在 $H$ 使 $\nabla H=\rho\nabla\Phi$。在每条连通的正则等值曲线上，$H$ 与 $\rho$ 都是常数。作者构造分支上的一维表示，并将其积分为 $I(t)=\int_t^{\max\Phi}A(s)\,ds$。重要的技术点是跨越临界值时不能简单假设分支延拓唯一；作者需要沿水平集分支追踪 Bernoulli 常数，并证明分支合并或分裂时积分公式仍成立。PDF 摘录说明论文第 2 节详细处理了该问题，但未提供该节全部引理证明，因此具体的分支追踪技术细节从提供内容中无法进一步确认。

第三，作者引入分支质量平方函数 $S(t)$。对每个分支 $D$，由 $-\Delta\Phi=\rho$ 和散度定理，分支内的源质量可转化为边界上的法向导数积分；对水平集使用余面积公式，再结合 Cauchy-Schwarz 与平面等周不等式，得到正则水平区间上的质量平方微分不等式。它的关键常数为 $8\pi$，并最终形成

$$ -dS=8\pi A(t)\,dt+\mu. $$

这里不把不等式损失丢弃，而是将其编码为非负测度 $\mu$。这一步是论证的核心创新之一：$\mu$ 同时检测两类不同机制，一类是正则水平集的几何不圆或相关不等式不取等，另一类是临界水平处正质量分支发生分裂所造成的质量平方下降。若分裂为质量 $m_1,m_2$ 的两个分支，$S$ 至少损失 $2m_1m_2$，因此在 $\mu$ 中产生正原子。于是 $\mu=0$ 不仅意味着几何估计处处取等的倾向，也排除了具有正质量的拓扑分支分裂。

接着，作者将局部的几何缺陷与无穷远分析连接起来。由于 $S(t)\to0$ 当 $t\uparrow\max\Phi$，可将累计缺陷写为 $\delta(t)=S(t)-8\pi I(t)=\mu((t,\max\Phi))$。因此证明目标化为在低水平找到 $\delta$ 趋于零的序列。定常情形中 $\rho\in L^1$，设总质量为 $\Gamma$，在 $t\to-\infty$ 时 $S(t)\to\Gamma^2$。作者通过减去径向对数项 $-\frac{\Gamma}{2\pi}\log|x|$，并在一列扩张圆环上证明余项梯度能量趋于零，建立全局 Pohozaev 恒等式。将该恒等式与质量平方分解比较，直接迫使 $\mu$ 总质量为零。

非零角速度下的困难实质不同：归一化源可写为 $\rho=b+\sigma\omega_0$，其中 $b=2|\Omega|>0$，它有常数背景，故 $\int\rho=\infty$，定常情形的有限总质量极限不可用。作者利用流函数中的二次项，证明足够低的等值线在大圆环上可以表示为径向图；再用边界 Pohozaev 恒等式将 $\delta(t)$ 控制在这些等值线上的角向能量。减去径向对数部分不改变角向导数，结合余项在适当圆环上的梯度能量衰减，并对穿过圆环的水平值取平均，选出 $t_k\to-\infty$ 使 $\delta(t_k)\to0$。这是处理无限背景质量而仍能消除缺陷测度的主要新步骤。摘要与证明纲要给出了此框架，但 PDF 摘录没有展示相关边界估计的完整常数和逐步推导，故这些更细的定量版本从提供内容中无法确认。

最后，$\mu=0$ 导致正则超水平集连通，并由等周等号条件推出适当水平集为圆盘。连通性使 $\rho$ 可以写成 $\rho=f(\Phi)$，其中 $f$ 为有界非负 Borel 函数，于是问题降为带有有界 Borel 非线性的半线性椭圆方程。作者调用 P.-L. Lions 的径向对称性定理，并利用嵌套超水平集确保各水平的圆盘具有共同中心；可能存在的外部区域满足 $\rho=0$，通过调和延拓处理。对 $\Omega\neq0$，流函数的二次远场项进一步使中心只能是预先给定的旋转中心原点；对 $\Omega=0$，平移不变性使中心可以自由选取。

4. 与之前工作的比较

与涡旋斑块刚性结果相比，本文最直接的推进是从具有清晰支撑边界的斑块问题转向全空间涡量剖面。Fraenkel 的结果针对定常、单连通涡旋斑块；Hmidi 对凸的单连通旋转斑块建立相应刚性；Gómez-Serrano、Park、Shi 和 Yao 后来在所述角速度范围内证明平面旋转涡旋斑块的径向刚性，不再需要凸性假设，并说明斑块问题中的阈值具有尖锐性。斑块模型的未知量主要由区域边界决定，而本文处理的是一般经典函数型涡量，涡量可在整个平面上非零，也可有非紧支撑尾部。因而本文的结论不依赖有限条支撑边界曲线或紧支撑边界条件，而依赖 $L^1$、$L^\infty$、$C^1$ 正则性及有效源的单侧符号条件。

与最接近的紧支撑光滑涡量结果相比，Fan、Wang 和 Zhan 的工作已对紧支撑 $C^2$ 旋转涡量在 $\Omega\leq\frac12\inf\omega_0$ 或 $\Omega\geq\frac12\sup\omega_0$ 时证明径向对称，并允许涡量变号。本文保留基本相同的角速度阈值结构，但将紧支撑与 $C^2$ 假设推进为 $C^1\cap L^1\cap L^\infty$ 的全空间框架。这个推进不是简单把紧支撑估计取极限：非紧支撑时没有可直接使用的外边界，定常情形要处理对数远场，旋转情形还出现常数背景源导致无限总质量。因此作者分别建立圆环上的渐近能量论证和低水平边界 Pohozaev 估计，以证明缺陷测度消失。代价是仍需 $C^1$ 正则性、可积有界性和固定符号有效源；论文并未声称适用于任意有界弱解、仅有限能量解或任意变号旋转剖面。

就参数范围而言，本文的条件可用椭圆方程解释为 $\omega_0-2\Omega$ 固定符号，而不只是单纯的符号涡量条件。对 $\Omega=0$，它覆盖单号定常涡量，并给出对称中心未指定的结论。对 $\Omega\neq0$，不论取哪一侧阈值，只要条件成立，对称中心都被固定在原点。这一中心定位与定常情况的平移自由度不同，反映了旋转剖面方程中指定的旋转轴。文中提及的紧支撑非负光滑旋转涡量结果曾覆盖 $\Omega\leq0$ 的情形；本文在其函数型全空间设定中给出更一般的阈值表述，但适用前提仍是有效源单号，不能据此推出阈值互补区域内也有刚性。

相较于 Ruiz 对紧支撑定常速度场、De Regibus、Esposito 和 Ruiz 对有限动能定常流的工作，以及 Hamel 和 Nadirashvili 基于流线和停滞集结构的 Liouville 型定理，本文采用的是针对流函数超水平集的分支质量方法。既有研究有的要求支撑呈环状或边界非退化，有的依赖停滞集连通性、有限动能或无穷远附加条件；本文的主假设则是涡量属于 $C^1\cap L^1\cap L^\infty$ 且有效椭圆源固定符号。需要注意，这并不意味着本文无条件分类所有定常流：作者在引言中明确指出，光滑变号定常涡量已有非径向紧支撑例子；其他文献也表明一般平滑稳态未必存在全局单值的涡量—流函数关系。因此本文通过角速度阈值排除了这类灵活性，而不是消除了所有非径向稳态。

与移动平面法或直接利用支撑边界的策略相比，本文先证明水平集连通，再把涡量表示为流函数的函数。这个次序尤其适合全空间非紧支撑问题，也避开了预先假设全局单值关系的风险。其代价是证明依赖二维等周不等式、Bernoulli 分支公式、Pohozaev 恒等式以及无穷远渐近分析；目前给出的结论是定性刚性，没有给出“接近圆形则接近径向”的定量稳定估计，也没有给出刚性阈值之外的分类或分岔结果。论文提及互补角速度范围内已有非径向旋转和周期解构造，故不能把本文结论外推到所有 $\Omega$。

最后，摘要和引言强调本结果对非紧支撑尾部的覆盖以及 Yudovich 框架的自然性，但没有声称其弱正则性达到 Yudovich 的最低可解性假设 $L^1\cap L^\infty$：刚性定理还要求 $C^1$。从提供内容中无法确认作者是否证明了角速度阈值在一般非紧支撑函数型涡量类别中的尖锐性；文中明确提到的尖锐性是既有涡旋斑块结果，而非本文对全部函数型剖面的尖锐性结论。

7. 具有非线性速度对齐的动力学群集模型的等温流体力学极限

Isothermal hydrodynamic limit for kinetic flocking models with nonlinear velocity alignment

作者: Roman Shvydkoy, Changhui Tan

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30457v1

摘要翻译: 我们研究具有非线性速度对齐和强局部 Fokker–Planck 强迫项的动力学群集模型的流体力学极限。当 Knudsen 数趋于零时，动力学密度收敛到一个局部 Maxwellian，其中心由满足带等温压力的可压缩 Euler 系统的宏观量给出，压力为 $p = \s ρ$；与此同时，宏观对齐项由微观速度律与标准高斯函数的卷积得到，从而产生一种不同的非线性。这种差异已在 \cite{black2025hydrodynamic} 中得到观察。

Comment: 无

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

1. 研究对象和背景

本文研究带局部 Fokker–Planck 松弛和非线性速度对齐的 Vlasov–Fokker–Planck 型动力学方程，在 Knudsen 数 $\varepsilon\to0$ 时的水动力极限。空间区域取单位环面 $\Omega=\mathbb T^d$，速度变量属于 $\mathbb R^d$，动力学密度 $f_\varepsilon(t,x,v)\ge0$ 满足

$$ \partial_t f_\varepsilon+v\cdot\nabla_x f_\varepsilon+\nabla_v\cdot\big(F(f_\varepsilon)f_\varepsilon\big) =\frac1\varepsilon\nabla_v\cdot\big((v-u_\varepsilon)f_\varepsilon+\sigma\nabla_v f_\varepsilon\big), $$

其中 $\sigma>0$ 固定，局部密度和平均速度由

$$ \rho_\varepsilon=\int_{\mathbb R^d}f_\varepsilon\,dv, \qquad \rho_\varepsilon u_\varepsilon=\int_{\mathbb R^d}v f_\varepsilon\,dv $$

定义。对齐力为

$$ F(f)(t,x,v)=\int_{\Omega\times\mathbb R^d}\varphi(x-x')\Phi(v'-v)f(t,x',v')\,dx'\,dv', $$

其中 $\varphi$ 是非负、偶对称的通信核，$\Phi$ 是奇的单调速度对齐律。代表性例子是 $p$ 阶对齐律 $\Phi(z)=|z|^{p-2}z$，$p\ge2$；$p=2$ 对应线性 Cucker–Smale 型相互作用。右端的局部 Fokker–Planck 项以速率 $1/\varepsilon$ 将速度分布拉向以局部平均速度为中心、温度由 $\sigma$ 固定的高斯平衡态。与无速度扩散的压力无关模型不同，这里极限平衡不是单速度的 Dirac 质量，而是局部 Maxwellian：

$$ \mu_\varepsilon(t,x,v)=\rho_\varepsilon(t,x)(2\pi\sigma)^{-d/2} \exp\!\left(-\frac{|v-u_\varepsilon(t,x)|^2}{2\sigma}\right). $$

这一区别决定了宏观极限保留等温压力。若 $f_\varepsilon$ 接近局部 Maxwellian，速度协方差趋向 $\sigma\rho_\varepsilon I_d$，因而宏观连续方程和动量方程的预期形式为

$$ \partial_t\rho+\nabla_x\cdot(\rho u)=0, $$
$$ \partial_t(\rho u)+\nabla_x\cdot(\rho u\otimes u)+\sigma\nabla_x\rho =\rho A_\Psi[\rho,u], $$

其中

$$ A_\Psi[\rho,u](x)=\int_\Omega\varphi(x-y)\Psi(u(y)-u(x))\rho(y)\,dy. $$

关键之处在于，非线性对齐律在 Maxwellian 平衡下并不直接保持为原来的 $\Phi$。两个局部高斯速度涨落之差仍是高斯分布，因此有效宏观对齐律为

$$ \Psi(z)=(4\pi\sigma)^{-d/2}\int_{\mathbb R^d}\Phi(z-\zeta) \exp\!\left(-\frac{|\zeta|^2}{4\sigma}\right)d\zeta. $$

例如 $d=1$、$\sigma=1$ 且 $\Phi(z)=z^3$ 时，$\Psi(z)=z^3+6z$。因此问题不仅是证明速度分布集中到 Maxwellian，还要识别热涨落对非线性力矩闭合的影响，并控制实际力矩与 Maxwellian 闭合之间的误差。

这一问题位于动力学集群模型与可压缩 Euler 型宏观模型之间。在线性对齐下，力矩可由密度和动量精确闭合；非线性对齐下，第一速度矩一般不能仅由 $\rho_\varepsilon,u_\varepsilon$ 表示。本文将该误差写成 Maxwellian discrepancy，并同时处理动量方程中的 Reynolds 应力缺陷。已有工作中，线性等温极限已有相应理论；非线性压力无关极限则可利用单速度平衡使热速度涨落消失。固定正温度的等温情形兼具非线性闭合缺陷与非消失的热涨落，本文的目标是在这一设定下，对满足熵条件的弱动理学解相对于给定光滑宏观解严格证明收敛，并给出相对熵定量速率。本文假设宏观极限解已存在且足够光滑，并不在所给结果中证明该解的存在时间或其形成机制。

2. 主要定理或主要结果

主要结论是定理 2.2。假设 (H1)–(H3) 成立：通信核满足 $\varphi\in W^{1,\infty}(\mathbb T^d)$、$\varphi\ge0$、$\varphi(x)=\varphi(-x)$；速度律 $\Phi$ 局部 $C^1$、奇且单调，满足多项式增长

$$ |\Phi(z)|\le C(1+|z|^{p-1}),\qquad |\nabla\Phi(z)|\le C(1+|z|^{p-2}),\qquad p\ge2, $$

以及对任意 $B>0$ 存在 $c_B>0,C_B\ge0$ 使得

$$ (\Phi(z)-\Phi(w))\cdot(z-w)\ge c_B|z-w|^p-C_B, \qquad |w|\le B, $$

并满足文中的散度控制条件 $|\nabla\cdot\Phi(z)|\le\delta z\cdot\Phi(z)+C_\delta$。假设极限系统存在 $[0,T]$ 上的光滑正密度解 $(\rho,u)$，且 $\rho,u,\nabla_xu,\nabla_x\log\rho$ 有界、总质量为 $1$。初始动理学数据质量为 $1$、自由能一致有界，并满足 well-prepared 条件

$$ \mathcal H(f_\varepsilon^0\mid\mu^0)\le C_0\varepsilon. $$

对每个 $0<\varepsilon\le1/2$，若 $f_\varepsilon$ 是满足定义 2.1 弱形式及熵不等式的 admissible weak solution，则存在与 $\varepsilon$ 无关的 $C_T$，使得

$$ \sup_{0\le t\le T}\mathcal H(f_\varepsilon\mid\mu)(t) +\frac{1}{4\varepsilon}\int_0^T I(f_\varepsilon(t))\,dt \le C_T\varepsilon\bigl(1+|\log\varepsilon|^{p-2}\bigr). $$

这里 $\mu=\rho(2\pi\sigma)^{-d/2}\exp(-|v-u|^2/(2\sigma))$，相对熵按文中约定带温度因子 $\sigma$；$I(f_\varepsilon)$ 是相对于局部 Maxwellian $\mu_\varepsilon$ 的调制 Fisher 信息：

$$ I(f_\varepsilon)=\int_{\mathbb T^d\times\mathbb R^d} \frac{|\sigma\nabla_v f_\varepsilon+(v-u_\varepsilon)f_\varepsilon|^2}{f_\varepsilon}\,dx\,dv. $$

因此结论还给出 $\int_0^T I(f_\varepsilon)\,dt\le C_T\varepsilon^2(1+|\log\varepsilon|^{p-2})$。由 Csiszár–Kullback 不等式，定理推出

$$ \sup_{0\le t\le T}\|f_\varepsilon(t)-\mu(t)\|_{L^1(\mathbb T^d\times\mathbb R^d)}^2 \le C_T\varepsilon\bigl(1+|\log\varepsilon|^{p-2}\bigr), $$

以及宏观量的 $L^1$ 收敛率

$$ \sup_{0\le t\le T} \left(\|\rho_\varepsilon-\rho\|_{L^1(\mathbb T^d)} +\|\rho_\varepsilon u_\varepsilon-\rho u\|_{L^1(\mathbb T^d)}\right) \le C_T\sqrt{\varepsilon}\bigl(1+|\log\varepsilon|^{(p-2)/2}\bigr). $$

方程的宏观力矩闭合关系也被明确识别。若定义

$$ R_\varepsilon=\int_{\mathbb R^d}(v-u_\varepsilon)\otimes(v-u_\varepsilon)f_\varepsilon\,dv -\sigma\rho_\varepsilon I_d, $$

及

$$ G_\varepsilon=B_\Phi[f_\varepsilon]-B_\Phi[\mu_\varepsilon], $$

则动量方程可写为

$$ \partial_t(\rho_\varepsilon u_\varepsilon) +\nabla_x\cdot(\rho_\varepsilon u_\varepsilon\otimes u_\varepsilon) +\sigma\nabla_x\rho_\varepsilon =\rho_\varepsilon A_\Psi[\rho_\varepsilon,u_\varepsilon] +G_\varepsilon-\nabla_x\cdot R_\varepsilon. $$

定理的相对熵估计控制了这些缺陷对宏观极限的影响。需注意，$f_\varepsilon$ 的弱解存在性不是定理 2.2 的结论；作者明确将完整非线性速度律类下此类 admissible weak solution 的构造留待后续工作。定理也以光滑正密度宏观解作为假设，不声称对一般弱宏观解成立。

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

证明的主要障碍不是单纯证明 Fokker–Planck 松弛使 $f_\varepsilon$ 接近 $\mu_\varepsilon$，而是控制非线性对齐力矩在速度平均后产生的闭合误差。对线性 $\Phi$，由于 $f_\varepsilon$ 与 $\mu_\varepsilon$ 在每个空间点具有相同的质量和一阶矩，力矩差 $G_\varepsilon$ 恒为零；对非线性 $\Phi$，速度涨落参与力矩，因而不能仅凭低阶矩相同消去误差。固定 $\sigma>0$ 时，平衡速度涨落本身也不趋于零，极限力律必须先经过高斯卷积形成 $\Psi$。

证明首先从弱解的自由能不等式建立统一先验估计。作者使用增长条件 $|\nabla\cdot\Phi(z)|\le\delta z\cdot\Phi(z)+C_\delta$，将熵产生项 $Q_\Phi$ 吸收到非负对齐耗散 $D_\Phi$ 中，得到自由能、二阶速度矩的一致控制，以及一般初值下 $\int_0^T I(f_\varepsilon)\,dt\lesssim\varepsilon$。well-prepared 初值下的更强 Fisher 信息估计则由主定理进一步给出。对弱解进行零阶和一阶速度矩测试，得到连续方程和含有 $R_\varepsilon,G_\varepsilon$ 的动量平衡式。应力项利用恒等式

$$ R_\varepsilon=\int_{\mathbb R^d} \bigl[2\sigma\nabla_v\sqrt{f_\varepsilon}+(v-u_\varepsilon)\sqrt{f_\varepsilon}\bigr] \otimes (v-u_\varepsilon)\sqrt{f_\varepsilon}\,dv $$

以及 Cauchy–Schwarz 估计为 $\|R_\varepsilon\|_{L^1_x}\lesssim I(f_\varepsilon)^{1/2}$；所需的内部二阶矩由自由能统一控制。

核心框架是相对于光滑极限解构造的等温相对熵。作者将其拆为动理学自由能与宏观修正项，后者由 $\rho_\varepsilon,u_\varepsilon$ 和极限量 $\rho,u$ 组成。将宏观修正项求导，再与动理学自由能不等式合并，压力项相消，从而避免直接对弱变量 $u_\varepsilon$ 求导。所得不等式中，惯性误差由宏观相对熵控制，应力缺陷由 Fisher 信息平方根控制；对齐部分则拆成宏观对齐耗散和非线性余项。宏观对齐耗散的非负性依赖 $\varphi$ 的对称性以及 $\Psi$ 的奇性、单调性。

真正困难的是把非线性余项与宏观对齐耗散放在一起估计，而不是孤立地对 $G_\varepsilon$ 作粗糙上界。作者先按动理学平均相对速度与极限相对速度之差划分空间点对。对两者差异较大的点对，局部 coercivity 使合并后的贡献具有非正符号，因而可丢弃；剩余点对上，动理学平均相对速度有界。随后再按速度相对变量划分有界区域与远场区域：远场动理学项保留有利符号，可舍去，只留下 Maxwellian 的高斯尾部；有界区域内，利用 transport–Fisher 型估计及分部积分恒等式，将非线性 Maxwellian discrepancy 归约到局部 Maxwellian Fisher 信息的控制。该安排利用了耗散结构，避免要求额外的一致高阶速度矩，这是论文特别强调的性质。

速度截断半径选在 $\sqrt{|\log\varepsilon|}$ 的尺度。这样高斯尾部随截断半径指数衰减，而有界区域内由 $p$ 阶增长导致的常数损失表现为 $|\log\varepsilon|^{p-2}$，最终得到 $\varepsilon(1+|\log\varepsilon|^{p-2})$ 的相对熵速率。这里的对数损失是非线性多项式增长估计与高斯尾部平衡的结果；$p=2$ 时该因子退化为常数。提供的文本在第 3 页概述了这一核心估计策略，但未给出第 5 节逐项估算的完整推导，因此不能据此补写更精细的中间不等式或各项常数依赖。

此外，对弱解的论证不是从光滑解的形式计算直接跳步：文章将测试函数的矩计算解释为通过速度截断逼近，在分布意义下成立；相对熵推导则以 admissible weak solution 所满足的熵不等式为出发点。论文说明完整非线性模型的弱解构造暂不处理，故这一证明是对满足定义 2.1 的解类作条件性极限论证，而不是同时建立该解类的存在性。

4. 与之前工作的比较

与线性等温极限相比，本文的直接参照是 Karper、Mellet 和 Trivisa 的工作 [12]。在线性速度律 $\Phi(z)=Lz$ 下，局部质量和动量已足以确定对齐力矩，非线性 Maxwellian discrepancy 消失；本文保留同一类强局部 Fokker–Planck 松弛和等温压力 $p(\rho)=\sigma\rho$，但允许满足多项式增长、单调性与 coercivity 条件的非线性速度律，并必须识别高斯平均后的有效宏观律 $\Psi$。作者指出 $p=2$ 时恢复已知线性结果。定理给出的统一表达式在 $p=2$ 时为相对熵 $O(\varepsilon)$；当 $p>2$ 时得到 $O(\varepsilon(1+|\log\varepsilon|^{p-2}))$，相应宏观 $L^1$ 误差为 $O(\sqrt\varepsilon(1+|\log\varepsilon|^{(p-2)/2}))$。非线性推广的代价是出现随非线性阶数增长的对数损失。

与 Black 和 Tan [2] 的非线性压力无关水动力极限相比，差别既在平衡态，也在极限闭合律。压力无关问题的强松弛产生单速度平衡 $f=\rho\delta_{v=u}$；速度涨落在极限中消失，因此形式上宏观非线性仍是原始 $\Phi$。本文研究固定 $\sigma>0$ 的等温模型，平衡为 Maxwellian，热涨落持续存在，故极限律变为 $\Psi=\Phi*\bar\mu_{2\sigma}$，且必须控制 $G_\varepsilon=B_\Phi[f_\varepsilon]-B_\Phi[\mu_\varepsilon]$。作者指出这一宏观非线性差异已在 [2] 中观察到；本文的推进是对 admissible weak kinetic solutions 和非真空光滑宏观解严格证明该等温闭合及定量相对熵收敛。形式上令 $\sigma\to0$ 时高斯卷积趋于 Dirac 质量、$\Psi\to\Phi$，可联系压力无关图景，但这只是模型间的形式关系，不应与本文定理的固定 $\sigma>0$ 结论混同。

模型背景方面，论文将线性速度对齐与 Cucker–Smale/Vlasov-alignment 文献联系，并提及 Figalli 和 Kang [7] 对线性压力无关极限的严格推导，以及 Black 和 Tan [2] 对非线性压力无关情形的拓展。本文相对于这些工作所处理的新组合，是局部热扩散导致的等温压力与非线性对齐闭合误差同时出现。其可处理范围由假设 (H1) 刻画：$\Phi$ 奇、单调，具有 $p$ 阶 coercivity 和相应多项式增长，并满足散度控制；这包含原型 $|z|^{p-2}z$。因此结果不只针对一个显式幂次律，但也不是任意非线性速度律的普遍定理。通信核取有界、对称且非负，宏观解要求光滑、正密度及有界导数，初值必须以 $O(\varepsilon)$ 相对熵 well-prepared；这些条件限制了结论的适用范围。

时间尺度上，结论是在预先给定的光滑宏观解存在区间 $[0,T]$ 上一致成立，常数 $C_T$ 可依赖于该区间上的宏观正则性等量，定理没有给出随 $\varepsilon$ 延长的时间区间，也不覆盖激波形成后的弱宏观解。空间上，正式定理在环面上陈述。作者提到全空间推广还需控制动理学熵负部并通过空间截断论证，例如借助类似 [12] 的约束条件；因此不能把本文环面结果直接宣称为全空间定理。

结论的逻辑也有明确代价：本文假设 $f_\varepsilon$ 已属于满足弱方程和熵不等式的 admissible 解类，并未在本文中为完整非线性速度律类别构造这些弱解；该构造被留作未来工作。与 [12] 中线性 Vlasov–Fokker–Planck 弱解理论及 [20] 对更一般通信方案的解构造、重整化和能量等式结果相关的背景，不能据此推断本文已完成非线性等温模型的存在性理论。本文的贡献主要是“给定 admissible 弱解后的定量极限证明”，而非动力学弱解存在性、宏观解存在性或一般弱解稳定性。

论文摘要和所提供正文还提及关于非线性平均场极限的 [16]、压力无关或 Euler-alignment 相关文献及其他模型背景，但提供内容未给出这些工作的完整定理假设与速率，故无法可靠逐项比较其误差阶或正则性条件。对这些未展开的文献比较，从提供内容中无法确认；作者也未在所给摘录中系统比较本文速率与所有相关工作的最优性。

8. 三维定常欧拉流附近的 Floquet 不稳定性

Floquet instability near a stationary Euler flow in 3D

作者: Mimi Dai

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30401v1

摘要翻译: 我们研究三维（3D）不可压缩欧拉方程在 Gavrilov 型光滑紧支撑定常流附近的线性化。Baldi 的保体积作用–角变量坐标图使得可以约化不变环面上的高频振幅方程。利用 Baldi 的坐标图，我们确定了由紧支撑截断保证存在的一个临界环面。我们证明，对于法向分量足够大的固定轴对称余向量，该环面上的约化无迹周期余循环是双曲的，并具有正的 Floquet 指数。随后，通过局部化的几何光学构造，我们将这种 Floquet 不稳定性提升到线性化欧拉方程，从而得到空间 $H^s$（$scgeq 0$）中的指数本质范数增长。

Comment: 无

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

1. 研究对象和背景

论文研究三维不可压缩欧拉方程在一类光滑紧支撑定常流附近的线性稳定性。方程为

$$ \partial_t u+(u\cdot\nabla)u+\nabla p=0,\qquad \nabla\cdot u=0,\qquad x\in\mathbb R^3. $$

所考察的定常解记为 $(U,P)$, 满足

$$ (U\cdot\nabla)U+\nabla P=0,\qquad \nabla\cdot U=0, $$

而关键未知量是其上的速度扰动 $v$ 及线性化压力 $q$. 线性化系统为

$$ \partial_t v+\nabla_Uv+\nabla_vU+\nabla q=0,\qquad \nabla\cdot v=0. $$

这里 $\nabla$ 是欧氏协变导数. 论文关注的不是一般定常流的抽象稳定性判据, 而是 Gavrilov 构造的轴对称、有非零旋流、支撑在细环形壳层内的特定紧支撑流, 并试图证明其线性化演化在 $H^s$ 中具有指数增长的本质范数.

Gavrilov 的构造从圆周附近的局部解析定常流出发. 对适当的压力变量 $P_0$, 可取光滑截断函数 $\chi$, 并令 $U=\chi(P_0)U_0$, $P=W(P_0)$, 其中 $W'=\chi^2$. 因为 $U_0\cdot\nabla P_0=0$, 截断沿不变压力面进行, 从而仍得到定常欧拉解. 选择支撑于正压力值的窄区间的 $\chi$, 可去除局部公式在核心圆周及外部区域的不光滑部分, 得到 $C^\infty_0$ 定常流. Constantin、La、Vicol 将这一思路置于 Grad–Shafranov 结构下系统化: 轴对称流函数满足相应的 Grad–Shafranov 方程, 并要求速度模长为流函数的函数, 以保证压力沿流线不变, 从而允许上述截断. 因而本文所研究的对象具有明确的几何结构, 但其紧支撑性也带来截断边界与作用区间的问题.

Baldi 对 Gavrilov 流核心附近的粒子运动给出了保体积的作用–角坐标. 在坐标 $y=(\sigma,\beta,I)$ 中, $I$ 为作用变量, 两个角变量分别以只依赖 $I$ 的频率匀速运动:

$$ \dot\sigma=\Omega_1(I),\qquad \dot\beta=\Omega_2(I),\qquad \dot I=0, $$

且坐标变换 $\Phi$ 满足 $\det D\Phi=1$. 不截断压力具有形式 $P_0\circ\Phi=K(I)$, 而局部化后的压力为 $P\circ\Phi=W(K(I))$. 对小作用量, Baldi 的频率关系为

$$ \Omega_1(I)=\chi(K(I))K'(I),\qquad \Omega_2(I)=\sqrt I\,R(I)\Omega_1(I), $$

其中 $K(I)=I+O(I^3)$, $R(I)=1+\frac74 I+O(I^2)$. 因此在非零截断作用区内, 流粒子沿不变环面运动, 而频率由截断函数及作用量共同决定.

研究这类线性化问题的重要背景是短波几何光学. 对高频扰动作 $v_\varepsilon\sim a(x,t)e^{iS(x,t)/\varepsilon}$, 可导出沿流粒子轨道演化的相位协向量与振幅系统, 即 bicharacteristic–amplitude 系统. 它把无限维线性化方程的高频主导行为编码为有限维线性 cocycle. Friedlander–Vishik 及 Lifschitz–Hameiri 的工作建立了这一短波分析框架, 并给出利用粒子流拉伸或振幅动力学判断流体不稳定性的途径. 但对本文的紧支撑 Gavrilov 流, 关键问题不是仅仅套用一般判据, 而是能否从截断结构本身找出一个明确的不稳定轨道, 并证明其振幅 cocycle 确实有正 Floquet 指数, 再将该有限维增长定量传递回精确的线性化欧拉演化。

本文所补的空白正是这一具体构造与传递链条. 截断使子午角频率在活动壳层内为正、在截断的平坦边界处为零, 因而产生内部临界环面. 该环面上的粒子仍在运动, 并非停滞集; 但作者选择轴对称高频协向量, 使其投影扭转在临界环面上消失, 从而得到需要细致分析的周期振幅系统. 作者据摘要及定理说明, 对固定轴对称协向量 $(n,0,\ell)$, 在法向分量 $\ell$ 足够大时, 临界环面上的约化无迹周期 cocycle 双曲, 并具有正 Floquet 指数. 随后以局部化几何光学将该指数增长提升为线性化欧拉演化在 $H^s$ 中的指数本质范数增长. 由此得到的机制不同于粒子轨道本身具有正 Lyapunov 指数, 也不依赖存在孤立的不稳定特征值。

2. 主要定理或主要结果

主要结论是定理 1.1. 固定截断参数 $\chi_*>0$, 对所有充分小的 $\iota$, 取论文构造的局部化 Gavrilov 定常流 $U_\iota$, 并记其线性化欧拉演化算子为 $G_\iota(t)$. 存在一个临界环面 $I_c$, 一个非零整数 $n$, 一个充分大的 $\ell>0$, 以及与轴对称协向量 $(n,0,\ell)$ 相联系的常数 $\sigma_\ell>0$, 使得对每个 $s\geq0$,

$$ \liminf_{t\to\infty}\frac1t\log\left\|G_\iota(t)\big|_{H^s_{\mathrm{axi}}}\right\|_{\mathrm{ess},H^s}\geq\sigma_\ell>0. $$

这里 $H^s_{\mathrm{axi}}$ 是 $H^s(\mathbb R^3;\mathbb R^3)$ 中轴对称、无散度向量场构成的闭子空间, $\|\cdot\|_{\mathrm{ess},H^s}$ 表示相对于 $H^s$ 的本质算子范数. 因而结论是本质范数至少以指数率增长, 而不只是某个有限维特征模态或单个测试函数的增长. 摘录未给出 $\sigma_\ell$ 关于 $\ell$, $\iota$ 或 $\chi_*$ 的显式公式, 也没有给出定量下界的具体常数, 因此不能据此补写其精确依赖关系。

定理 1.1 还推出时间依赖测试数据的版本: 对每个 $0<\delta<\sigma_\ell$ 及所有充分大的观测时刻 $T$, 存在实值、光滑、轴对称、无散度的 $v_{0,T}$, 满足 $\|v_{0,T}\|_{H^s}=1$, 且

$$ \|G_\iota(T)v_{0,T}\|_{H^s}\geq e^{(\sigma_\ell-\delta)T}. $$

此处初值可以依赖于观测时刻 $T$, 因而这条结论本身不是对某个预先固定初值的逐时刻指数下界。

针对固定初值, 推论 1.2 使用一致有界原理与 Baire 范畴论给出泛性结论. 对每个 $s\geq0$, 存在实 Banach 空间 $H^s_{\mathrm{div}}(\mathbb R^3;\mathbb R^3)$ 中的剩余集 $\mathcal R_s$, 使得每个非零 $v_0\in\mathcal R_s$ 都满足

$$ \limsup_{t\to\infty}\frac1t\log\bigl(1+\|G_\iota(t)v_0\|_{H^s}\bigr)\geq\sigma_\ell. $$

轴对称实无散度子空间内也有相同的剩余集结论. 因此泛性固定初值的 $H^s$ 轨道沿某些趋于无穷的时刻至少达到该指数增长率; 尤其 $\sup_{t\geq0}\|G_\iota(t)v_0\|_{H^s}=\infty$. 这里是 limsup 结论, 不能加强表述为所有大时间均有相同下界。

摘要和定理还说明有限维来源: 在临界环面上, 对固定轴对称协向量 $(n,0,\ell)$ 且法向分量 $\ell$ 充分大, 约化的无迹周期 cocycle 为双曲系统, Floquet 指数为正. 论文的第 4 节标题进一步表明作者分析了法向极限 cocycle、其恒等 monodromy、首个法向修正及其平均生成元, 然后证明临界环面上的 Floquet 双曲性. 但提供摘录没有列出该平均生成元矩阵、Floquet 指数的展开式、余项阶或 $\ell$ 的显式阈值, 故这些数值型结论从提供内容中无法确认。严格证明的结论是线性化演化的本质范数指数增长及上述泛性固定数据结论; 非线性不稳定性只被列为 companion paper 的研究方向, 并非本文在所给内容中证明的结果。

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

证明的关键障碍是从一条特殊的高频射线得到可验证的振幅增长, 并确保这种有限维增长不是坐标退化或截断边界造成的伪效应. 作用–角坐标在核心 $I=0$ 退化, 而紧支撑截断在活动区外平坦消失; 所以分析必须限制在严格包含于活动区的闭环形带上, 同时又要利用截断造成的频率变化来找到临界环面. Corollary 2.2 正是提供这一框架: 对 $[I_-,I_+]\Subset G_\chi$, 坐标映射是解析、保体积的微分同胚, 且背景速度的坐标分量为 $(\Omega_1(I),\Omega_2(I),0)$. 这使后续计算可以避开核心坐标奇性和截断的平坦边界。

第一步是把线性化欧拉算子准确写入保体积作用–角坐标. 由于 $\det D\Phi=1$, 体积元为 $d\sigma\,d\beta\,dI$, 向量场 $v=v^ae_a$ 的散度简化为 $\partial_av^a$. 但协变导数并不因此变成平直导数: 速度变形项仍含有由坐标度量产生的 Christoffel 符号. 具体地,

$$ (\nabla_Uv+\nabla_vU)^a= \Omega\cdot\nabla_\theta v^a+v^I\partial_IU^a+2\Gamma^a_{bc}U^bv^c, \qquad \Omega=(\Omega_1,\Omega_2). $$

作者因而先展开坐标映射 $\Phi$、度量 $g_{ab}$、逆度量和 Christoffel 符号的小作用量渐近式. 摘录给出 $r=\sqrt{2I}$ 下 $\rho,\zeta,\gamma$ 的 $r$ 展开, 以及对应的几何系数. 这些计算的作用不是单纯整理坐标, 而是为振幅方程中的连接项提供足够精确的首个非平凡贡献; 保体积恒等式 $\det D\Phi=1$ 同时约束展开系数, 用来核对迹与散度结构。

第二步是提取轴对称角模并构造高频 bicharacteristic–amplitude 系统. 论文第 3 节引入 angular modes 与 localized frequency profile, 第 4.1 节沿作用–角射线写出位置、相位协向量和速度极化振幅的耦合系统, 第 4.2 节再将约化动力学压到二维平面. 由摘要可确认, 作者固定协向量的角向分量为零, 选择 $(n,0,\ell)$, 并要求法向分量 $\ell$ 足够大. 这样得到的重点是临界环面上的 $2\times2$ 周期无迹振幅 cocycle. 由于摘录未提供该 cocycle 的完整系数公式, 不能进一步声称具体的矩阵元素或某个未展示的能量恒等式。

真正的 Floquet 分析采用逐阶法向渐近. 第 4.3 节处理法向极限 cocycle, 并得到其 monodromy 为恒等矩阵; 这意味着仅看极限系统并不会给出双曲增长, 必须辨认使 monodromy 从单位矩阵偏离的首个法向修正. 第 4.4 节计算这一首修正的平均生成元, 第 4.5 节据此找出截断产生的临界环面并证明 Floquet 双曲性. 从摘要可确定这一分析的核心结论是: 对固定轴对称协向量、充分大的法向频率, 临界环面上的约化 trace-free periodic cocycle 为双曲并具有正 Floquet 指数. 这一步的实质是把一个极限上退化为恒等 monodromy 的问题, 转化为考察首个有效平均修正的特征方向与符号; 但平均矩阵的具体形式及误差控制阶数在提供内容中没有呈现, 不宜补造。

临界环面的产生也具有结构性. 截断使 $\Omega_1(I)$ 在活动壳层内非零, 而在平坦边界处趋于零. 因此频率剖面须在内部出现临界作用量. 作者利用该临界点使选定轴对称模式的投影扭转消失, 同时粒子仍以非零速度沿环面运动. 这避开了把不稳定性归因于停滞点或粒子流指数拉伸的解释, 并使周期振幅系统能在固定环面上作 Floquet 分析。

最后一步是从短波振幅增长传递到精确线性化欧拉演化. 作者在第 5 节对 Floquet 不稳定射线作局部化几何光学构造, 目标是把振幅指数增长实现为 $H^s$ 中的高频扰动, 并据此获得本质算子范数的指数下界. 这一“本质”结论排除了仅由有限维紧扰动或孤立模态解释增长的可能. 摘要称该构造将 Floquet 不稳定性 lift 到线性化方程; 引言则把它置于 Vishik 的算子转移理论及 Shvydkoy 的伪微分框架背景下. 提供内容未包含第 5 节的误差估计、波包尺度选择、截止函数参数或近似解残差阶, 因而这些技术细节无法确认。得到算子增长后, 再以一致有界原理与 Baire 范畴论将随时间变化的测试向量结论转成剩余集中的固定数据 limsup 增长, 这是定理到推论的函数分析环节, 并不声称每个初值均不稳定。

4. 与之前工作的比较

与 Gavrilov 的工作相比, 本文不再以构造紧支撑定常解本身为目标, 而是以其局部化流为背景研究线性化稳定性. Gavrilov 的关键贡献是通过沿压力不变面截断局部定常流得到光滑紧支撑解; 本文利用的正是这种截断对作用–角频率的影响, 进一步识别由紧支撑局部化所强制产生的内部临界环面. Constantin、La、Vicol 将相同思路推广为 Grad–Shafranov 框架下的局部化构造, 并讨论定常解的正则性与能量耗散等问题; 本文的问题则转向特定平滑紧支撑解上的高频线性动力学, 没有据此重新改进定常解存在或正则性门槛。Peralta-Salas 与 Slobodeanu 关于局部可截断稳态流的轴对称性及凸环面几何结果, 属于稳态结构分类背景, 并非本文 Floquet 不稳定性证明的直接替代结果。

与 Baldi 的工作相比, 本文把 Baldi 提供的保体积作用–角坐标作为分析基础, 而不是重新建立粒子动力学的可积性. Baldi 描述粒子轨道、坐标映射及频率 $\Omega_1,\Omega_2$; Dai 则在这些坐标中计算线性化欧拉振幅方程的几何连接项, 并研究特定协向量上的周期 cocycle. 关键新增内容是把频率截断 profile 的临界点与振幅 Floquet 双曲性联系起来, 再证明相应高频增长能在精确线性化演化中留下本质范数增长. 两篇工作的定量关系, 例如 Baldi 展开系数在最终 Floquet 指数中的显式贡献, 从提供内容中无法确认。

与 Friedlander–Vishik 及 Lifschitz–Hameiri 的短波理论相比, 本文沿用高频几何光学把线性化 PDE 化为 bicharacteristic–amplitude 系统这一基本路线, 但推进点是对一个明确的紧支撑 Gavrilov 族实际构造不稳定射线, 而不是假设或抽象要求存在具有正增长率的射线. 经典轴对称可积流判据常以流线曲率、测地挠率、螺度等量构造充分不稳定条件; 本文强调的是截断生成的临界环面与足够大的法向协向量所带来的 Floquet 双曲性. 其背景粒子流在该可积环面上的经典 Lyapunov 指数为零, 但受散度约束的 Euler 振幅 cocycle 仍有正 Floquet 指数. 因此该机制不要求粒子轨道本身具有指数分离, 也不依赖孤立不稳定特征值。Friedlander–Vishik 关于零粒子 Lyapunov 指数时仍可能出现流体振幅增长的既有观察, 已在相关文献中出现; 本文的贡献是对局部化 Gavrilov 族给出定量的临界环面与 Floquet 双曲性论证, 而不是首次提出这类概念可能性。

在从有限维 cocycle 到 PDE 增长的层面, Vishik 建立了振幅 cocycle 与演化算子本质谱半径之间的算子理论联系, Shvydkoy–Vishik 讨论三维 Euler cocycle 的 Lyapunov–Oseledets 指数与演化群谱的联系, Shvydkoy 则发展平流型方程的伪微分传递框架及其后续谱理论. 本文不只是报告一个振幅指数, 而是把该机制用于局部化构造, 给出每个 $s\geq0$ 的 $H^s$ 本质范数指数增长下界, 且可先限制到实轴对称无散度子空间. 与仅有谱不稳定或有限时刻模态增长相比, 本质范数结论表达了高频方向上的非紧增长; 但提供摘录没有展示作者对上述一般传递定理所作的具体改进程度, 也没有足够信息比较其常数或正则性代价。

结果范围方面, 本文结论覆盖所有 $s\geq0$, 但依赖特定的局部化 Gavrilov 流、固定截断参数及充分小的 $\iota$, 并需选取非零角向整数 $n$ 与充分大的法向分量 $\ell$. 证明因此是构造性的但具有明显参数约束, 不是对所有三维定常欧拉流的普遍不稳定性定理. 对单个观测时刻可以选取依赖 $T$ 的测试波包并获得 $e^{(\sigma_\ell-\delta)T}$ 下界; 对固定初值, 所给结论只是在剩余集上成立的 limsup 指数增长, 不能据此推出每个非零初值都增长, 也不能推出固定初值对所有充分大时间均满足逐时刻指数下界. 摘录没有给出 $\sigma_\ell$ 的显式大小、临界作用量位置、可允许 $\ell$ 的定量阈值或几何光学误差阶, 故这些比较维度无法确认。最后, 引言提及非线性分析由 companion paper 处理; 本文给出的严格结论是线性不稳定及其固定数据推论, 不应与非线性不稳定性证明混同。

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

1. 具有长程相互作用的一维 Hartree 方程的渐近行为

Asymptotic Behavior of One-Dimensional Hartree Equations with the Long-Range Interaction

作者: Changhun Yang, Chanjin You

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31515v1

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

摘要翻译

中文标题: 具有长程相互作用的一维 Hartree 方程的渐近行为

摘要翻译: 我们研究具有软库仑相互作用的一维 Hartree 方程的小解的长时间行为。Hartree 方程是描述多体量子系统演化的有效平均场模型，而软库仑势则提供了库仑相互作用的一维正则化近似：它消除了原点处的奇异性，同时保留了其长程特性。

由于相互作用核在空间中的衰减缓慢，该方程呈现长程非线性效应。我们证明，对于充分小且指数局部化的初始数据，全局存在性和修正散射成立。一个显著特点是，与众所周知的三维库仑 Hartree 方程所产生的对数相位修正不同，一维软库仑相互作用产生的主导相位修正具有 $(\log t)^2$ 阶，随后还有一个低阶的对数修正。

证明在解析函数空间中进行。初始数据的指数局部化转化为变换后剖面的解析性，并采用一种基于具有时间依赖解析半径的生成函数的方法，以补偿非线性相位方程所产生的导数损失。

2. 无散场的海森堡不确定性原理及二阶 Caffarelli–Kohn–Nirenberg 不等式的最优稳定性层级

Optimal stability hierarchies of the Heisenberg Uncertainty Principle for solenoidal fields and of the second order Caffarelli–Kohn–Nirenberg inequalities

作者: Anh Do, Tuan Duong, Nguyen Lam, Guozhen Lu, Van Hoang Nguyen

类别: math.AP, math.CA

发布日期: 2026-09-25

arXiv ID: 2609.31579v1

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

摘要翻译

中文标题: 无散场的海森堡不确定性原理及二阶 Caffarelli–Kohn–Nirenberg 不等式的最优稳定性层级

摘要翻译: 2018 年，V. Maz’ya 提出了 75 个分析与偏微分方程中的公开问题。其中一个问题源于流体力学中的相关问题，涉及无散向量场的海森堡不确定性原理的最佳常数。Cazacu、Flynn 和 Lam 在二维情形回答了这一问题，随后 Hamamoto 通过建立具有显式常数的尖锐不等式，将结果推广至所有维数 $N\ge3$。然而，难度更大的稳定性问题在所有维数下仍未解决。本文的第一个主要结果是建立一个尖锐的稳定性层级：我们证明其亏损项控制到极值函数集合的距离，其中 $N\ge4$ 时尖锐稳定性常数为 $\frac12(N-\sqrt{N^2-4N+12})$，$N=3$ 时为 $1$；此外，我们还给出一系列余项，它们度量到由显式给出的更大极向场与环向场族的距离。第二个主要结果是，在直线 $1+a+b=0$ 上，针对权重 $|x|^{-2a}$ 和 $|x|^{-2b}$ 建立二阶 $L^2$-Caffarelli–Kohn–Nirenberg 不等式，并得到其稳定性常数 $\frac12\min\{4(1+a),\sqrt{(N+2a)^2+4N-4}-(N+2a)\}$。当第二个数较小时，该 CKN 不等式是尖锐的。我们发展了一种新方法：将与球面模态相联系的四阶一维问题通过实数阶 Fourier–Hankel 变换转化为一阶不等式，而这些不等式的尖锐稳定性正是加权 Gaussian Poincaré 不等式的尖锐稳定性。逆变换则给出以 Kummer 函数表示的极值函数。作为这一新方法的附带结果，我们还以显著更简洁的方式证明了 Hamamoto 的一维不等式及其关于无散场的尖锐不确定性原理，并给出了尖锐余项和等号成立情形。

3. 测度值映射的半线性 Neumann 边值问题

Semilinear Neumann boundary value problems for measure-valued maps

作者: Aleksei Kroshnin, Hugo Lavenant, Dmitry Vorotnikov

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

发布日期: 2026-09-25

arXiv ID: 2609.31353v1

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

摘要翻译

中文标题: 测度值映射的半线性 Neumann 边值问题

摘要翻译: 我们研究具有非齐次 Neumann 边界条件的非线性 Poisson 系统的 Wasserstein 提升。我们利用涉及测度值迹的 $p$ 阶矩的新估计，证明极小化子的存在性。我们还推导出最优性条件；这些条件可视为无压 Euler 方程的矩阵值推广。在适当的凸性假设下，以及当源区域为一维时，确定性极小化子存在；而我们证明，在更高维情形会出现真正的测度值极小化子。特别地，我们构造了若干例子，其中提升后的变分问题与经典变分问题具有不同的最小值，从而排除了提升问题的解为确定性的可能性，尽管并不存在任何先验机制来强制产生测度值行为。据我们所知，这一现象在非线性椭圆理论中尚属首次。

4. 黎曼流形上 $u_{t}=Δu^{m}+au\log u+bu$ 的 Hamilton 梯度估计与 Liouville 定理

Hamilton’s gradient estimates and Liouville theorems for $u_{t}=Δu^{m}+au\log u+bu$ on Riemannian manifolds

作者: Jun Sun, Jiaming Yang

类别: math.AP, math.DG

发布日期: 2026-09-25

arXiv ID: 2609.31393v1

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

摘要翻译

中文标题: 黎曼流形上 $u_{t}=Δu^{m}+au\log u+bu$ 的 Hamilton 梯度估计与 Liouville 定理

摘要翻译: 本文应用 Nash-Moser 迭代和 Saloff-Coste 的 Sobolev 不等式，为完备黎曼流形 $(M,g)$ 上一类如下形式的非线性抛物方程的正解推导梯度估计： \begin{equation*} u_{t}=Δu^{m}+a(x, t)u\log u+b(x, t)u, \quad 1-\frac{2}{n}<m<1 \quad \text{or}\quad 1<m<1+\frac{1}{1+\sqrt{2n}}, \end{equation*} 其中流形的维数为 $n$，且 $a(x, t)$、$b(x, t)$ 是 $C^{1}$ 函数。当 $a=b=0$ 时，若 $m>1$，该方程为多孔介质方程（PME）；若 $m<1$，则为快速扩散方程（FDE）。作为推论，对于相应的常系数方程，当 \begin{equation*} 1-\frac{1}{\sqrt{n-1}}<m<1\quad \text{and}\quad 1<m<1+\frac{1}{\sqrt{n-1}}. \end{equation*} 时，我们得到 Liouville 型定理。

5. 任意指数下黎曼流形上多孔介质方程与快速扩散方程的 Hamilton 型梯度估计

Hamilton-Type Gradient Estimates for Porous Medium and Fast Diffusion Equations on Riemannian Manifolds for all exponents

作者: Jun Sun, Jiaming Yang

类别: math.DG, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31373v1

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

摘要翻译

中文标题: 任意指数下黎曼流形上多孔介质方程与快速扩散方程的 Hamilton 型梯度估计

摘要翻译: 我们建立了具有下有界 Ricci 曲率的完备黎曼流形上方程 $u_t=Δu^m$ 的正 $C^{2,1}$ 解的局部 Hamilton 型梯度估计。对于每个固定的 $m>1$ 以及每个固定的 $0<m<1$，存在 $β=β(m,n)>0$ 和 $C=C(m,n)>0$，使得在 $B_{2R}(x_0)\times(t_0-T,t_0]$ 中满足 $0<u\leqslant A$ 的解满足以下局部梯度估计：\begin{equation*} \sup_{B_R(x_0)\times(t_0-T/2,t_0]}|\nabla u^β| \leqslant C A^β\left(\frac1R+\sqrt{k}+\frac{A^{(1-m)/2}}{\sqrt T}\right), \end{equation*} 其中 $\mathrm{Ric}_M\geqslant-k$。证明中采用一种内禀的定量两分法，在每个指定的正值点周围找到一个解的上、下界之比受控的柱体。随后将该柱体上的局部梯度估计与停止论证相结合。由此可得，在具有非负 Ricci 曲率的连通完备流形上，每个一致有界的正古代解都是常数。

6. 准线性薛定谔方程的数学研究

A mathematical study of a quasi-linear Schr{ö}dinger equation

作者: Meriem Bahhi, Jonas Lampart, Simona Rota Nodari

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31139v1

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

摘要翻译

中文标题: 准线性薛定谔方程的数学研究

摘要翻译: 我们研究一种准线性薛定谔型方程，该方程源于对原子核内粒子动力学的描述。在模型参数满足适当假设的条件下，我们证明了正径向解的存在性、唯一性和非退化性。此外，我们还分析了这些解对参数的依赖性。

7. 关于二维中的 Shnirelman 不等式：不可逆行为与 Euler-Lagrange 变分表述

On Shnirelman’s inequality in 2D: Irreversible behavior and Euler-Lagrange variational formulation

作者: Stefan Schiffer, Martina Zizza

类别: math.AP, math.OC

发布日期: 2026-09-25

arXiv ID: 2609.31294v1

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

摘要翻译

中文标题: 关于二维中的 Shnirelman 不等式：不可逆行为与 Euler-Lagrange 变分表述

摘要翻译: 本文证明了二维正方形 $[0,1]^2$ 上 Shnirelman 不等式的一个版本。更确切地说，给定一个目标流体构型 $f$（一个保体积微分同胚），我们构造了一个散度为零的速度场 $v\in L^q_tL^p_x$，其流将 $f$ 与恒等映射连接起来，并且其 $L^q_{t} L^p_x$-范数实际上可以由 $f$ 与 $id$ 的 $L^p$-差控制。然而，向量场的更高正则性可能会受到影响。与高维情形的主要区别在于，由于维数 $ν=2$ 中的拓扑障碍，我们必须允许不同的流体轨迹在空间中“相交”。我们还观察到，这一选择会导致欧拉描述下的动力学中出现不可逆行为。

8. 基于 Cosserat 梁理论的蛇类运动仿生控制模型

Bio-inspired control model for snake locomotion based on Cosserat beam theory

作者: Eliot Thys, Frédéric Boyer, Yacine Chitour, Petri Kokkonen, Vincent Lebastard, Swann Marx, Clément Moreau

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31036v1

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

摘要翻译

中文标题: 基于 Cosserat 梁理论的蛇类运动仿生控制模型

摘要翻译: 本文从数学角度重新审视文献 \cite{boyer2010poincare} 和 \cite{boyer2011macrocontinuous} 中给出的流体环境下蛇形身体运动建模。利用 Cosserat 杆理论，我们根据约束和驱动的类型推导出若干组偏微分方程系统。

9. 三维凸域中散焦立方薛定谔方程的大数据全局适定性

Large-Data Global Well-Posedness for the Defocusing Cubic Schrödinger Equation in 3D Convex Domains

作者: Len Meas

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30920v1

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

摘要翻译

中文标题: 三维凸域中散焦立方薛定谔方程的大数据全局适定性

摘要翻译: 我们证明了三维 Friedlander 模型区域上的能量次临界散焦立方非线性薛定谔方程（NLS），在齐次 Dirichlet 边界条件下，对于强制能量空间 $H_0^1(Ω)$ 中任意大的初始数据，具有全局适定性。该构型的相空间几何具有严格凸边界，并存在非空掠射集；广义耳语廊焦散会将高频波包限制在边界层内。这些集中现象使尖锐线性 Strichartz 估计中出现内在的导数损失，从而构成一项重要的微局部障碍，使得无法通过经典微扰框架直接在能量层面闭合非线性 Duhamel 迭代。 为弥合守恒律与线性理论之间的正则性缺口，我们构造了一族适配边界的连续–离散 Bourgain–Strichartz 限制空间 $X^{s,b}$，该空间基于模型 Friedlander 算子的 Dirichlet 实现的谱分解。法向变量通过离散 Airy 谱模态加以分解，而切向变量则进行连续傅里叶分析。在这一函数空间框架内，我们实施精细的高–低频分解，在次能量空间 $H_0^s(Ω)$ 中为高频余项建立一个受扰动的非线性方程，其中 $ \frac{1}{2}<s<1. $ 初始高频数据满足如下定量次能量衰减估计：$ \|P_{>λ}u_0\|_{H_0^s(Ω)} \label{eq:abs_decay} \lesssim λ^{s-1}\|u_0\|_{H_0^1(Ω)}, $，这提供了一个抵消导数损失的小参数。

10. 关于 Tzitzéica 方程柯西问题的研究：孤子分解猜想与渐近分析

On the Cauchy problem for the Tzitzéica equation: Soliton resolution conjecture and asymptotic analysis

作者: Shou-Fu Tian, Jia-Fu Tong

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

发布日期: 2026-09-25

arXiv ID: 2609.30858v1

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

摘要翻译

中文标题: 关于 Tzitzéica 方程柯西问题的研究：孤子分解猜想与渐近分析

摘要翻译: 我们研究 Tzitzéica 方程的柯西问题。该方程是一个重要的可积模型，源自仿射微分几何，并刻画真仿射球面。近期，Huang、Wang 和 Zhu (Math. Ann. \textbf{395}, 18 (2026)) 报道了纯连续谱情形下 Tzitzéica 方程的长时间渐近结果。受他们工作的启发，我们对 Tzitzéica 方程的孤子分解猜想与渐近分析进行了深入研究。与以往工作相比，通过发展 $bar{\partial}$-非线性最速下降法，我们揭示了 Tzitzéica 方程的解依赖于 $ξ:=x/t$ 呈现出三个不同的渐近区域。对于不存在驻相点的区域 $ξ\in(-\infty,-1]\cup[1,\infty)$，我们严格证明，Tzitzéica 方程的解以代数速度衰减至零，误差为 $O(t^{-1})$。在区域 $ξ\in(-1,1)$ 中，相位函数 $θ(z)$ 有两个驻相点。相应的渐近近似可由一个 $N$-孤子解以及孤子解与色散项之间的相互作用项刻画，其剩余误差阶为 $O(t^{-3/4})$。在过渡区域 $\{(x,t):1-\varepsilon<|ξ|<1\}$ 中，相应的渐近近似可由 $N$-孤子解连同孤子解与色散项之间的相互作用项刻画，其剩余误差阶为 $O(t^{-3/4})$。我们的结果对 Tzitzéica 方程解在不同区域的长时间渐近行为给出了严格分析，并证实了 Tzitzéica 方程的孤子分解猜想。

11. 部分谐势约束下非线性薛定谔方程的逆散射问题

An inverse scattering problem for the nonlinear Schrödinger equation under partial harmonic confinement

作者: Pranav Kumar, Li Li

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30806v1

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

摘要翻译

中文标题: 部分谐势约束下非线性薛定谔方程的逆散射问题

摘要翻译: 我们研究一类非线性薛定谔方程的逆问题：该方程在某些空间方向上具有谐势，在其余方向上则自由传播。非线性项具有未知幂次和实值空间依赖系数。对于加权能量空间中的小入射数据，我们构造波算子，并证明当系数依赖于受约束变量或自由变量时，该波算子都能确定幂次和系数。证明利用波算子的首个非线性项以及集中型乘积数据。我们说明，系数依赖于自由变量时，幂次需要满足更强的下界。我们还讨论出射波算子、小数据散射算子，以及依赖于所有空间变量的系数。

12. 波速与初始源同时反演的利普希茨稳定性

Lipschitz stability in the simultaneous recovery of wave speed and initial source

作者: Shiqi Ma, Qi Zhai

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30767v1

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

摘要翻译

中文标题: 波速与初始源同时反演的利普希茨稳定性

摘要翻译: 我们考虑一个双曲方程的反问题。目标是根据一组狄利克雷–诺伊曼柯西数据，同时反演空间变系数（其平方根表示波速）和初始位移。通过引入适当的变换并对所得方程关于时间求导，我们推导出一个关于未知系数差与初始位移差的耦合系统。结合加权能量估计与 Carleman 估计，我们在适当的 Sobolev 范数下，为这两个未知量建立了 Lipschitz 稳定性估计。证明依赖于先验正则性假设、初始数据的非退化条件，以及控制系数差与初始位移差之间相互作用的结构性不等式。我们的结果将仅针对系数的稳定性结果推广到了系数与初始位移同时识别的问题。

13. 双调和非线性薛定谔方程的长时间约化

Long-time reduction for the biharmonic nonlinear Schrödinger equation

作者: Younghun Hong, Junyeong Jang

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30681v1

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

摘要翻译

中文标题: 双调和非线性薛定谔方程的长时间约化

摘要翻译: 本文研究高阶色散消失极限下，一维双调和非线性薛定谔方程向三次非线性薛定谔方程的约化。在中间正则性范围 $0<s<2$ 内，能量守恒律无法控制 $H^s$–范数；我们证明了 $H^s$–解的长时间 $L^2$–收敛，并给出了随时间呈指数增长的逼近误差界。论证主要依赖于正则性保持性。该结果提供了一个简单例子，说明在极限问题中，无需依赖高阶守恒律也可运用正则性保持性。

14. 消失黏性 Hamilton-Jacobi 方程周期均匀化中的收敛速率

Convergence rates in the periodic homogenization of vanishing-viscous Hamilton-Jacobi equations

作者: Kai Qin

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30693v1

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

摘要翻译

中文标题: 消失黏性 Hamilton-Jacobi 方程周期均匀化中的收敛速率

摘要翻译: 我们研究动量 Hessian 在每一点均正定的 Hamilton-Jacobi 方程中，周期均匀化与消失黏性在同一尺度发生时的收敛速率。对于有界 Lipschitz 初值和一类光滑 Hamiltonian，我们在固定时间区间上建立了最优收敛速率 $O(\varepsilon|\log\varepsilon|)$。二阶项发挥着关键作用：椭圆型胞元问题为有效 Hamiltonian $\bar{H}$ 及其 Legendre 对偶 $\bar{L}$ 提供了高阶导数界。证明完全基于 PDE，采用校正子展开和黏性比较；其核心思路则来自 Hopf–Lax 公式的控制论解释。对于固定目标点 $(x,t)$，一个极小化起点选定特征速度及其对偶动量；它们引导我们构造光滑的比较剖面，而这些剖面的梯度为一阶校正子的动量变量提供取值。因此，即使在有效解不光滑的点，我们也无需对其求导。

15. 楔形体绕流稳态超声速势流中斜音速激波的稳定性

Stability of oblique sonic shocks in steady supersonic potential flow past a wedge

作者: Geng Lai

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30666v1

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

摘要翻译

中文标题: 楔形体绕流稳态超声速势流中斜音速激波的稳定性

摘要翻译: 当均匀稳态超声速来流冲击直楔形体时，若楔角小于脱体角，流场中将形成两类满足熵条件的稳态斜激波：下游流为超声速或亚声速的弱激波，以及下游流为亚声速的强激波。关于下游流为亚声速或超声速的斜激波稳定性，已有许多研究结果。然而，由于激波后方的流动对扰动非常敏感，人们对下游流为音速的斜激波稳定性所知甚少。本文在楔形体为直线且来流非均匀的假设下，研究斜音速激波的稳定性。我们将该问题化约为一个退化双曲型自由边界问题。利用经典的阿斯科利-阿尔泽拉定理和对角线论证，构造出一个使下游流呈音速-超声速状态的自由边界问题的局部 Lipschitz 连续解。该自由边界问题的主要困难在于音速线上出现的双曲退化。通过采用加权特征分解方法，我们建立了下游流在音速线附近的一个精细估计。

16. 费米-帕斯塔-乌拉姆-秦古系统的长时间 Korteweg-de Vries 近似

Long-time Korteweg-de Vries approximation for the Fermi-Pasta-Ulam-Tsingou system

作者: Younghun Hong, Junyeong Jang

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30668v1

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

摘要翻译

中文标题: 费米-帕斯塔-乌拉姆-秦古系统的长时间 Korteweg-de Vries 近似

摘要翻译: 我们证明，当晶格间距 $h$ 趋于零时，在长度为 $$ 的时间区间上，无限费米–帕斯塔–乌拉姆–秦古（FPUT）系统的一般解可以在 $L^2$ 中由两个反向传播的 Korteweg–de Vries（KdV）波近似，其中时间区间的长度为 $$。这解决了第一作者及其合作者提出的一个开放问题~\cite{HKY2021}。我们的证明结合了 FPUT 守恒律、$L^2$ 中关于 $h$ 一致的局部适定性理论以及 Sobolev 正则性的保持性。我们还引入一个频率局部化的辅助方程，以克服比较与 FPUT 流和 KdV 流相关的 Fourier 限制范数时遇到的困难。

17. H 型群上 Folland-Stein 不等式的极值函数与稳定性

Extremals and stability for the Folland-Stein inequality on H-type groups

作者: Zhipeng Yang

类别: math.AP, math.CA

发布日期: 2026-09-24

arXiv ID: 2609.30531v1

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

摘要翻译

中文标题: H 型群上 Folland-Stein 不等式的极值函数与稳定性

摘要翻译: 我们分类了 H 型群上尖锐 $L^2$ Folland-Stein 不等式的所有实值极值函数，并证明了一个全局 Bianchi-Egnell 稳定性估计。标准泡的线性化核恰好由其平移导数和伸缩导数组成。我们还得到了对称模态之上谱隙的显式界，并刻画了局部稳定性渐近展开中的最优系数。上述结果适用于任意正交 Clifford 模。

18. 具有对数律的伪相对论薛定谔方程基态：定性性质与非相对论极限

Ground states for pseudo-relativistic Schrödinger equations with logarithm-law: qualitative properties and nonrelativistic limit

作者: Pingkang Guan, Yuanhong Wei

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30516v1

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

摘要翻译

中文标题: 具有对数律的伪相对论薛定谔方程基态：定性性质与非相对论极限

摘要翻译: 本文研究具有对数非线性的伪相对论薛定谔方程 \begin{equation*} \left(\sqrt{-c^{2}Δ+m^{2}c^{4}}-mc^{2}\right)u+μu=u\log|u| \quad\text{in }\mathbb{R}^{N}. \end{equation*} 通过幂律近似，我们建立了作用量基态的存在性及其定性性质。随后证明，当 $c\to\infty$ 时，这些态在 $H^{1}(\mathbb{R}^{N})$ 中收敛到 Gausson。对于充分大的 $c$，我们进一步证明，作用量基态除平移和符号改变外唯一，且具有非退化性。最后，我们证明作用量基态与能量基态之间存在精确对应关系，从而将上述结果推广到任意给定质量下的能量基态。

19. 薛定谔演化的相位恢复

Phase retrieval for Schrödinger evolutions

作者: Ben Pineau, João P. G. Ramos, Mitchell A. Taylor

类别: math.AP, math.CA

发布日期: 2026-09-24

arXiv ID: 2609.30511v1

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

摘要翻译

中文标题: 薛定谔演化的相位恢复

摘要翻译: 考虑一维薛定谔演化

$$ i\partial_tw+\partial_x^2w=Vw,\hspace{5mm}(t,x)\in [0,T]\times\mathbb{R}, $$

其中 $V$ 是取实值的时变势，属于 $L^1_tL^\infty_x+L^2_{t,x}$。我们证明，上述方程的一个 $L^2$ 解由其质量密度 $|w|^2$ 唯一确定，但相差一个全局模为一的相位因子。对于一大类势，这回答了 Jaming 最近提出的问题，并推广恢复了他关于自由薛定谔演化的早期结果。我们的证明简单而稳健。我们结合上述薛定谔演化的基本局部平滑性质、Strichartz 估计，以及 Ionescu 和 Kenig 唯一延拓定理的逐条带版本，证明只要 $u$ 和 $v$ 是具有势 $V$ 的薛定谔方程的 $L^2$ 解且具有相同的质量密度，交互函数

$$ F(t,x,y)=u(t,x)v(t,y)-u(t,y)v(t,x) $$

就在 $[0,T]\times\mathbb{R}^2$ 上恒等于零。该方法可直接推广至若干重要的非线性模型，例如三次 NLS 方程；取 $V=|u|^2$ 时，该方程在 $L^2(\mathbb{R})$ 中整体适定。有趣的是，我们还证明，上述相位恢复性质在二维及更高维中失效，即使对于具有 Schwartz 初值的自由方程也如此。

20. Verigin问题的一个相场模型：弱解与锐界面极限

A phase-field model for the Verigin problem: weak solutions and sharp-interface limit

作者: Anna Kubin, Alice Marveggio

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30412v1

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

摘要翻译

中文标题: Verigin问题的一个相场模型：弱解与锐界面极限

摘要翻译: 我们提出一种新的相场模型，用以近似含相变的Verigin问题。对于固定的界面厚度 $$，我们通过一种时间离散变分格式构造弱解；该格式将针对密度的Wasserstein最小化运动与针对相场的Schätzle选择过程相结合，从而得到强紧性和几乎极小性质。随后，我们研究锐界面极限 $$。对于准备良好的初始数据，我们建立密度和相场的强紧性，并证明时间积分的弥散界面能收敛到周长。我们将极限演化识别为含相变的Verigin问题的一个弱解，并证明其满足最优能量耗散不等式。

21. 三维波导上共振双曲薛定谔流的尖锐临界 Strichartz 估计

Sharp critical Strichartz estimates for a resonant hyperbolic Schrödinger flow on the 3D waveguide

作者: Xi Cen, Zhezhi Zhang

类别: math.AP, math-ph

发布日期: 2026-09-24

arXiv ID: 2609.30437v1

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

摘要翻译

中文标题: 三维波导上共振双曲薛定谔流的尖锐临界 Strichartz 估计

摘要翻译: 我们研究波导 $$ 上共振双曲薛定谔流的临界 Strichartz 估计。对于频率局限在至多 $N$ 的初始数据，我们证明单位时间 $L^2\to L^{10/3}$ 算子范数的尖锐阶为 $N^{1/5}$。我们进一步刻画其对观测时长的依赖，得到短时间和长时间区域中的尖锐估计，并刻画临界估计保持无损的频率相关时间区间。在中间区域 $N^{-1}<T<1$ 中，关于 $T$ 的精确依赖仍未解决。上界通过将椭圆柱面 $$ 上的全局端点估计提升到 Hilbert 值数据，同时保留其时间可和性，并结合周期 Bernstein 不等式与插值获得。下界则利用周期频率格点中的一个精确整数零方向，以及一个经过适当缩放的欧几里得波包。最后，我们通过 Bloch–Floquet 转移论证，从紧致环面理论恢复已知的椭圆型单位时间估计，并在频率相关的短时间区间上得到一个无损的椭圆型估计。

不相关论文（基本信息）

1. 到达时间方程的最优正则性

Optimal regularity for the arrival time equation

作者: Yiqi Huang, Jingze Zhu

类别: math.DG, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31610v1

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

摘要: For a convex mean curvature flow in $\mathbb{R}^{n+1}$, the corresponding arrival time function solves a degenerate elliptic equation that becomes singular at extinction. We prove that it is $C^{24,α}$ for any $0<α<1$ in the plane, with a logarithmic modulus for its twenty-fourth derivatives. We prove that this bound is sharp combining regularity obstruction in our companion paper. We obtain asymptotic expansions of the rescaled flow to any prescribed order, including derivative estimates for the remainder. Combined with the work of Sesum for $n\ge 2$, this also gives the sharp $C^{2,2/n}$ regularity for convex arrival time functions in higher dimensions.

2. 到达时间方程的一般正则性障碍

General regularity obstructions for the arrival time equation

作者: Yiqi Huang, Jingze Zhu

类别: math.DG, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31599v1

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

摘要: We study the regularity of the arrival time function associated with mean curvature flow, formulated as a degenerate elliptic equation encoding the singular structure of the flow. We identify a systematic mechanism obstructing higher regularity, arising from higher-order asymptotic expansions near spherical singularities. We construct uncountably many low-regularity arrival time functions in all dimensions. In particular, we settle the question of smoothness in the planar case by showing that even for convex curve shortening flow the arrival time need not be $C^{25}$. Our companion paper proves that every such arrival time is $C^{24,α}$ for any $0<α<1$. Hence the obstruction at order $25$ gives the optimal regularity threshold. Our approach introduces new analytic ingredients, including a precise correspondence between elliptic asymptotic expansions and parabolic long-time asymptotics of the associated rescaled mean curvature flow, a complexification of the arrival time equation, and a method for prescribing higher-order asymptotic expansions.

3. 基于等参超曲面的伯努利自由边界问题奇异性模型

Singularity models for the Bernoulli free boundary problem from isoparametric hypersurfaces

作者: Benjy Firester, Raphael Tsiamis, Zihui Zhao

类别: math.AP, math.DG

发布日期: 2026-09-25

arXiv ID: 2609.31617v1

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

摘要: We develop a general construction of homogeneous solutions to the Bernoulli free boundary problem, as well as general extremal domains on the sphere, from isoparametric foliations of the sphere. Our construction produces rich families of infinitely many new examples with sophisticated topologies connected to minimal surfaces of the sphere by smooth families of interpolating capillary surfaces, which include novel singularity models in low dimensions and recover most known homogeneous one-phase solutions. We also introduce a geometric representation for every such homogeneous map and establish a geometric rigidity theorem: for every prescribed isoparametric foliation, and in particular for cohomogeneity-one subgroups of $O(n)$, its radial geometry uniquely determines the corresponding homogeneous solution. Finally, we study the density of the constructed solutions and conjecture that the lowest density among non-flat cones in a given dimension $n \geq 5$ is attained by an $O(k) \times O(n-k)$-invariant cone.

4. 标量 Hencky 塑性中的周期均匀化

Periodic homogenization in scalar Hencky plasticity

作者: Gilles A. Francfort, Alessandro Giacomini

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31404v1

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

摘要: This paper investigates the periodic homogenization of a two-phase elasto-plastic material in the framework of $Γ$-convergence. It establishes the link between the homogenized functional and that which was previously obtained in [12] in the context of two-scale convergence. Although there is no gap between the two formulations and that involving two-scale has the hallmark of a (two-scale) elasto-plastic problem, we are unable to positively ascertain the elasto-plastic character of the homogenized energy except in a one-dimensional setting.

5. 自由边界极小曲面与常平均曲率曲面的拓扑与狄利克雷谱

Topology and Dirichlet spectrum of free boundary minimal and CMC surfaces

作者: Alcides de Carvalho, Roney Santos

类别: math.DG, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31304v1

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

摘要: We find bounds for the genus and for the number of boundary components of a surface that is either index one as a free boundary minimal surface or stable as a free boundary constant mean curvature surface in a Riemannian three-manifold with convex boundary and non-negative Ricci curvature, provided the first Dirichlet eigenvalue of its Jacobi operator is non-negative. As an application, we establish strong topological restrictions on such surfaces in compact strictly convex domains of three-dimensional normal homogeneous spaces.

6. 具有双相 Orlicz 增长的泛函的局部极小化子的部分正则性

Partial regularity for local minimizers of functionals with double-phase Orlicz growth

作者: Wenrui Chang, Shenzhou Zheng

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31242v1

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

摘要: In this paper, we consider the local minimizers to a class of non-autonomous functional with double-phase Orlicz growth: \begin{equation*} u \in W^{1,1}(Ω; \mathbb{R}^N) \mapsto \int_Ω\Big( {\varphi_1}\big(|Du|{\mathbb{A}^u}\big) + a(x),{\varphi_2} \big(|Du|{\mathbb{A}^u}\big) \Big), dx, \end{equation*} where $|Du|_{\mathbb{A}^u} := \big\langle \mathbb{A}(x,u)\,Du, Du \big\rangle^{\frac{1}{2}}$ with $\mathbb{A}(x,u) = \big\{A^{αβ}_{ij}(x, u)\big\}_{i,j = 1\cdots N}^{α,β= 1\cdots n}$ as the uniformly elliptic bounded symmetric tensor field, and $\varphi_k(\cdot)$ for $k=1,2$ are two different $N$-functions. We prove the partial regularity of their local minimizers if the continuity modulus of $a(\cdot)$ and $\varphi_k(\cdot)$ are fit for the gap conditions, and $\mathbb{A}(x,u)$ meets an optimal regularity in $(x,u)$. This is an extension from the single Orlicz growth functional to the double-phase one based on an essential improvement of several important inequalities.

7. 边界诱导重置下的二维目标搜索

2D target search with boundary-induced resetting

作者: Kevin Chen, Paul C. Bressloff

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

发布日期: 2026-09-25

arXiv ID: 2609.31085v1

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

摘要: Most studies of stochastic search processes with resetting focus on spontaneous resetting events that occur at a random sequence of times typically generated by a Poisson process. One of the major consequences of spontaneous resetting is that it yields a finite mean first passage time (MFPT) in unbounded domains, resulting in an optimal resetting rate that minimises the MFPT. An optimal resetting rate can also occur in bounded domains provided that the reset point is not too far from the target. Recently, there has been growing interest in event-based rather than spontaneous resetting, where resetting is triggered when a specific threshold is reached. Identifying the threshold with a physical non-target boundary then leads to so-called boundary-induced resetting. In this paper, we explore the effects of boundary-induced resetting on single-particle diffusive search in a two-dimensional (2D) bounded domain $Ω$ containing one or more small interior targets. The interior target boundaries are totally absorbing, whereas the exterior boundary $\partial Ω$ is sticky. That is, whenever the particle reaches a point on $\partial Ω$ it remains attached for a random waiting time $τ$ after which it immediately resets to a fixed interior point $\x_0\in Ω$. Using renewal theory, matched asymptotic analysis and Green’s functions we calculate the unconditional MFPT for absorption by any of the targets and compare this with the corresponding MFPT of a search process on the same domain with a reflecting surface $\partial Ω$ and no resetting. We show that boundary-induced resetting only reduces the MFPT relative to the no resetting case if the searcher’s reset point $\x_0$ is sufficiently close to one of the targets.

8. 位移单调线性-二次 Nash 系统的长时间与大群体极限

Long-time and large-population limits of displacement monotone linear-quadratic Nash systems

作者: Davide Francesco Redaelli

类别: math.OC, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31091v1

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

摘要: For linear-quadratic $N$-dimensional Nash systems with mean-field-like scaling and strongly displacement monotone data, we prove uniform-in-time a priori derivative estimates exhibiting the expected scaling as $N \to \infty$ in order to pass the system to the limit. As a consequence, in a heterogeneous (i.e., non-symmetric) mean-field setting, we obtain both uniform-in-time convergence of the Nash system on any arbitrarily long horizon $[0,T]$ (as $N \to \infty$), and convergence of both the $N$-dimensional Nash system and the corresponding master equation to their respective ergodic counterparts (as $T \to \infty$).

9. 半球中的非齐次曲率流

Non-homogeneous curvature flows in a hemisphere

作者: Hongyi Sheng, Weimin Sheng, Jiazhuo Yang

类别: math.DG, math.AP

发布日期: 2026-09-25

arXiv ID: 2609.31023v1

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

摘要: Let S^{n+1}{+} be the open hemisphere of the unit sphere S^{n+1} centred at o. We study the non-homogeneous curvature flow X_t=-f(r) sigma_k^{alpha} {nu} of smooth, closed, strictly convex hypersurfaces enclosing o, where r is the geodesic distance to o. We consider both the supercritical regime beta>1+k {alpha} and the critical regime beta=1+k{alpha}, where beta is the growth order of the profile at the origin, f(r) almost equals r^{beta} as r descends to 0. Under the structural condition that f^{1/(1+k{alpha})} is convex, we prove long-time existence and preservation of strict convexity. The normalized radial function converges smoothly and exponentially to a constant: to 1 and to R{infty}>0, resp. in different two cases. Thus the normalized radial graphs become round, while the original hypersurfaces contract to o.

10. 非线性 Stein 定理的统一方法

A unified approach to nonlinear Stein theorems

作者: Arka Mallick, Swarnendu Sil

类别: math.AP, math.FA

发布日期: 2026-09-25

arXiv ID: 2609.30949v1

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

摘要: We study regularity for the variable exponent quasilinear $\mathfrak{p}\left(x\right)$-Laplace type equation on domains in Heisenberg groups and Euclidean spaces. We establish borderline continuity estimates for the appropriate first order derivatives of solutions. More precisely, we show that for any weak solution $u \in \mathbb{E}W^{1, \mathfrak{p}(\cdot)}\left(Ω\right)$ of \begin{align*} \operatorname{div}{\mathbb{E}} \left( \mathfrak{a}(x) \lvert \nabla{\mathbb{E}} u \rvert^{\mathfrak{p}\left(x\right)-2} \nabla_{\mathbb{E}} u \right) = f \qquad \text{ in } Ω, \end{align*} where $Ω\subset \mathbb{E}^{n}$, $\mathfrak{a}$ is a uniformly positive bounded scalar function and the exponent function $\mathfrak{p}$ is uniformly bounded away from $1$ and $\infty,$ $\nabla_{\mathbb{E}}u$ is continuous in $Ω$ as soon as $f \in L^{\left( Q_{\mathbb{E}}, 1\right)}\left(Ω\right)$ and $\mathfrak{a}, \mathfrak{p}$ satisfies some conditions regarding the summability of their mean-oscillations. Here $\mathbb{E}^{n}$ is either the Heisenberg group $\mathbb{H}_{n}$ or the Euclidean space $\mathbb{R}^{n}$ and $\operatorname{div}_{\mathbb{E}}$, $\nabla_{\mathbb{E}}$, $Q_{\mathbb{E}}$ stands for the corresponding divergence, gradient and homogeneous dimension, respectively. We treat the elliptic and subelliptic cases in a unified manner using Euclidean techniques and our conditions on $\mathfrak{a}$ and $\mathfrak{p}$ are new and weaker than all the known sufficient conditions even in the Euclidean case. However, all the known sufficient conditions imply our conditions, achieving yet another unification.

11. 时间分数阶演化方程阶数反演的稳定性

Stability in recovering the order of a time-fractional evolution equation

作者: Ravshan Ashurov, Masahiro Yamamoto

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30799v1

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

摘要: We consider the recovery of the fractional order in a time-fractional evolution equation with a source of the form $p(t)f$, where the temporal factor $p$ is also unknown but belongs to a bounded admissible class and has a prescribed nonzero value at $t=0$. From one scalar observation we prove a uniform small-time expansion and obtain a Hölder estimate for the fractional order. More precisely, if a fixed derivative order $γ$ lies below every admissible fractional order, then the stability exponent is $α_+/(2α_+-γ)$; for the undifferentiated observation this gives the exponent $1/2$. We also show that differentiating the data to an order larger than all admissible fractional orders is supercritical: the corresponding derivative is unbounded near $t=0$ whenever the two fractional orders are different. A uniform estimate for the two-parameter Mittag–Leffler function, including the diagonal case, is proved in the appendix.

12. 初始值未知时确定时间源因子和分数阶次的唯一性

Uniqueness in determining the temporal source factor and the fractional order with an unknown initial value

作者: Ravshan Ashurov, Masahiro Yamamoto

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30794v1

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

摘要: We consider inverse problems for a time-fractional diffusion equation with a second-order elliptic operator and the homogeneous Dirichlet boundary condition in a bounded domain of $\mathbb R^N$. The initial value is unknown and is not assumed to be zero. The main result concerns the simultaneous determination of the fractional order and the time-dependent source factor from one weighted spatial observation. For irrational orders in $(0,1)$ and bounded source factors satisfying an exponential integrability condition, we prove that a nonzero observation given for all positive times uniquely determines both unknowns. We also describe the information about the initial value determined by this observation. A counterexample shows that uniqueness of the source factor may fail at the classical order $α=1$, even when the observation is nonzero. In the second part, we establish two uniqueness results for the fractional order when the elliptic operators, initial values, and source terms in the compared equations may be different. The first result is proved by the Laplace transform, whereas the second uses short-time observations and the asymptotic behavior of the Mittag–Leffler functions.

13. 几乎周期介质中 Fisher-KPP 方程过渡前沿平均速度的一个上界及前沿宽度无界性

An Upper Bound for the Mean Speed of Transition Fronts and Unbounded Front Widths for Fisher KPP Equations in Almost Periodic Media

作者: Xing Liang, Linfeng Xu, Qi Zhou, Tao Zhou

类别: math.AP

发布日期: 2026-09-25

arXiv ID: 2609.30760v1

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

摘要: In this paper, we investigate transition fronts and spreading solutions of Fisher–KPP equations in one-dimensional almost periodic media, both on the real line and on the lattice. Let $λ_1$ be the supremum of the spectrum of the linearized operator acting on $L^2(\R)$ or $\ell^2(\Z)$, respectively, and let $L(λ_1)$ denote the spatial Lyapunov exponent at the spectral parameter $λ_1$. We prove that, if $L(λ_1)>0$, the global mean speed of any transition front is at most $λ_1/L(λ_1)$. Moreover, if a solution of the Cauchy problem spreads faster than this bound, its transition width is unbounded along a sequence of times. This occurs for a class of initial data with slowly decaying exponential tails. This contrasts with periodic media, where pulsating fronts exist at every speed above the minimal speed. As applications, we consider almost Mathieu coefficients and a continuous quasiperiodic coefficient with two frequencies. Together with the results of Nadin–Rossi and Liang–Wang–Zhou–Zhou \citep{LWZZ24}, our results provide an almost complete picture of the admissible speeds of generalized transition fronts, while leaving the critical cases open.

14. 线性抛物型方程的定量衰减

Quantitative Decay for Linear Parabolic Equations

作者: Igor Kukavica, Qi Xu

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30630v1

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

摘要: We establish sharp temporal lower bounds for the full spatial $L^2$-norm of energy solutions to $\partial_tu-Δu=vu$ on $\mathbb{R}^n$, with $n\geq3$. For global solutions, nonvanishing at a single time implies the lower bound $ce^{-Ct}$ when $v$ is small in the scale-invariant space $L^\infty_tL^{n/2}_x$, and $ce^{-Ct^2}$ throughout the subcritical range $v\in L^\infty_tL^p_x$, with $p>n/2$, without smallness. The linear time exponent is optimal already for the free heat equation, while we show that the quadratic exponent is optimal for bounded complex potentials on $\mathbb{R}^3$. No additional regularity of the potential is required. The proof uses a convexified mixed Carleman estimate whose norms are designed together with the final time truncation. Integration over spatial centers removes the heat-kernel weight and yields quantitative comparisons of the full solution norms. To prove quadratic sharpness in the Euclidean energy class, we construct a time-dependent function $g$ that localizes Meshkov’s periodic parabolic example while keeping the induced potential bounded across the zeros of the periodic profile.

15. 流行率与抗原演化速率在生态-演化 SIR 模型中共同决定

Prevalence and the rate of antigenic evolution are jointly determined in an eco-evolutionary SIR model

作者: Daniel Gomez, Joshua B. Plotkin

类别: q-bio.PE, math.AP, math.DS, nlin.PS

发布日期: 2026-09-24

arXiv ID: 2609.30515v1

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

摘要: Influenza A persists by continually changing the antigens that its host population recognizes. How fast it evolves and how many hosts it infects are not independent, because the immunity that selects new variants is produced by the infections the pathogen has already caused. We analyze a model of Andreasen et.\ al.\ that captures this eco-evolutionary feedback, and we extend it to the case where infection may kill the host. In the model infected and recovered hosts are distributions over a one-dimensional antigenic space coupled by cross-immunity. The endemic state is a traveling pulse: the infected distribution keeps a fixed shape and advances through antigenic space at constant speed. In the limit of weak antigenic mutation we construct the pulse asymptotically and find that: (i) the speed of antigenic evolution decreases as prevalence rises, (ii) less standing antigenic variation accompanies either faster or slower antigenic evolution, depending on whether pathogens differ in immune evasion or host lifespan (iii) the total number of infected hosts is largest at an intermediate virulence when the host birth rate is sufficiently low.

16. 忆阻型 FitzHugh–Nagumo 神经网络的连续数据同化

Continuous Data Assimilation for Memristive FitzHugh–Nagumo Neural Networks

作者: Jing Tian

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30480v1

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

摘要: In this study, we investigate continuous data assimilation for a class of partly diffusive FitzHugh–Nagumo neural networks with memristive feedback and linear synaptic coupling. The nudging algorithm implemented here only adds a feedback term to the membrane-potential equations, while the recovery and memductance variables remain unobserved. We obtain weighted energy estimates that give sufficient conditions for exponential convergence of the data-assimilated solution to the reference solution under full or coarse membrane potential observations. Moreover, we prove energy-class uniqueness in dimensions one and two. In dimension three, the uniqueness criterion is an $L^4$ membrane-potential initial data. We also establish an observational-noise estimate for locally square-integrable forcing, yielding an asymptotic error bound when the noise energy is eventually bounded.

17. 非退化梯度 Young 测度与非线性膜的降维

Nondegenerate gradient Young measures and dimension reduction for nonlinear membranes

作者: Gokul G. Nair, David Padilla-Garza, Marco Picchi Scardaoni

类别: math.AP

发布日期: 2026-09-24

arXiv ID: 2609.30411v1

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

摘要: We consider the rigorous derivation of membrane theories from three-dimensional nonlinear elasticity. A salient feature of our framework is the inclusion of an orientation-preservation constraint together with energy growth that penalizes vanishing local volume. We embed the 3D variational problem into a space of parametrized measures, and obtain a membrane $Γ$-limit defined on a class of nondegenerate gradient Young measures. This formulation offers a twofold advantage: first, it uniquely identifies a non-relaxed membrane energy density capable of capturing the fine oscillations (microstructure) of minimizing sequences; second, the resulting membrane density preserves the unbounded energy growth near degenerate configurations. We further show that the classical relaxed membrane energy of Le Dret and Raoult is recovered as the minimum of our functional over all measures with a prescribed barycenter. The main difficulty lies in generating such measures by maps whose gradients have rank two almost everywhere, which we address through piecewise isometric (origami-like) constructions.