生成时间: 2026年10月22日 08:02:11

运行状态: 完成

计划处理论文数量: 48

已完成论文数量: 48

未完成论文数量: 0

重点关注论文数量: 7

了解领域论文数量: 16

不相关论文数量: 25

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

1. 二维非齐次不可压 Euler 方程的适定性理论及一类全局时间解

A well-posedness theory for the 2D inhomogeneous incompressible Euler equations and a class of global-in-time solutions

作者: Ibrahim Suleiman

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.30096v1

摘要翻译: 我们研究二维非齐次不可压 Euler 方程的 Cauchy 问题。我们证明，在临界正则性下，对于一类 $$m$$ 重对称函数，且允许速度场在空间无穷远处至多线性增长，该 Cauchy 问题局部适定。这是 P.L. Lions 在齐次、但不要求对称的情形下提出的问题。一个关键要素是变系数椭圆理论，该理论通过拟共形映射以及某些尺度不变 Hölder 空间中的 Schauder 估计建立。局部适定性使我们能够导出满足 $$0$$ -齐次涡量和密度的解在 $\mathbb{S}^1$ 上所满足的一维系统。对于该系统，我们可以证明一类远离平衡态的初值具有整体适定性。据我们所知，此前尚无任何非微扰的时变解类被证明对二维非齐次 Euler 方程具有整体适定性。最后，对于某些带符号的全局解，我们证明在无限时间极限下，密度均匀化与平均弛豫之间存在二分性。

Comment: 47 pages, 6 figures

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

1. 研究对象和背景

本文研究二维全平面 $\mathbb{R}^2$ 上密度可变、速度不可压的理想流体。未知量为密度 $\rho(t,x)>0$ 、速度 $u(t,x)\in\mathbb{R}^2$ 和压力 $$p(t,x)$$ ，控制方程为

\begin{cases} \partial_t\rho+u\cdot\nabla\rho=0,\\ \partial_t(\rho u)+u\cdot\nabla(\rho u)+\nabla p=0,\\ \nabla\cdot u=0. \end{cases}

给定散度为零的初始速度 $$u_0$$ 和有正下界、有限上界的初始密度 $0<c_0\leq \rho_0(x)\leq c_1$ ，问题是建立 Cauchy 问题的存在性、唯一性与连续依赖理论，并进一步确定某些非平凡解是否能对所有时间存在。与常密度 Euler 方程不同，变密度情形中通常不能仅靠涡量的输运获得闭合估计；密度梯度与压力项会进入涡量演化，因此二维常密度问题的全局理论并不能直接移植。

适合本文分析的涡量变量不是速度涡量 $\omega=\nabla^\perp\cdot u$ ，而是动量涡量 $\Omega=\nabla^\perp\cdot(\rho u)$ 。令 $u=\nabla^\perp\psi$ ，则 $\Omega=\nabla\cdot(\rho\nabla\psi)$ ，形式上可将原方程改写为

\begin{cases} \partial_t\rho+u\cdot\nabla\rho=0,\\ \partial_t\Omega+u\cdot\nabla\Omega=\dfrac12\nabla|u|^2\cdot\nabla^\perp\rho,\\ \nabla\cdot(\rho\nabla\psi)=\Omega,\qquad u=\nabla^\perp\psi. \end{cases}

因此，核心反演问题是由 $\Omega$ 和 $\rho$ 恢复流函数 $\psi$ ，也就是处理系数随空间变化、且系数还随时间演化的椭圆方程。全平面上允许速度在无穷远处线性增长，源项也不要求衰减，这使通常依赖衰减条件的 Biot-Savart 积分公式不能直接适用。本文将离散旋转对称性作为补偿条件：若 $$O_m$$ 表示转角 $2\pi/m$ 的旋转，则标量函数满足 $$f(O_mx)=f(x)$$ ，向量场满足 $$v(O_mx)=O_mv(x)$$ 。作者在尺度不变的 Hölder 空间 $\mathring C^{k,\gamma}(\mathbb{R}^2)$ 中工作；这些空间刻画角向正则性，允许函数不随半径衰减，且与线性增长的速度尺度相容。

研究的另一条主线是径向齐次解。这里的“径向齐次”指密度与动量涡量为零次齐次，而流函数为二次齐次：

\rho(t,r,\theta)=P(t,\theta),\qquad \Omega(t,r,\theta)=g(t,\theta),\qquad \psi(t,r,\theta)=r^2G(t,\theta).

代入后，二维演化约化为圆周 $\mathbb{S}^1$ 上的非局部一维系统：

\begin{cases} \partial_tg+2G\partial_\theta g=-(4G^2+(\partial_\theta G)^2)\partial_\theta P,\\ \partial_tP+2G\partial_\theta P=0,\\ \partial_\theta(P\partial_\theta G)+4PG=g. \end{cases}

常密度时，该约化退化为已有的径向齐次二维 Euler 方程所对应的一维模型；变密度时新增的非局部非线性源项正是全局估计的主要困难。既有工作已在多种正则性设定下建立非齐次 Euler 的局部理论，也有延拓判据和长时间存在结果，但一般大数据情形的全局适定性仍是开放问题。常密度二维 Euler 的全局理论则有 Wolibner、Yudovich 等经典结果，其机制依赖涡量输运及相应的速度反演。本文试图补足两处空白：一是在允许线性增长、临界尺度正则性下建立变系数椭圆反演与局部适定性；二是在一个不要求小扰动、但具有径向齐次和特定对称、单调结构的数据类中证明整体时间适定性。

2. 主要定理或主要结果

变系数椭圆反演（定理 1.4）。设 $0<\gamma<1$ ， $\rho\in\mathring C^{1,\gamma}_m(\mathbb{R}^2)$ ， $0<\inf\rho$ ， $\|\rho\|_{\mathring C^{1,\gamma}}\leq M$ ，且 $f\in\mathring C^{0,\gamma}_m(\mathbb{R}^2)$ 。定义密度对比度

\alpha(\rho)=\frac{\sup_{\mathbb{R}^2}\rho}{\inf_{\mathbb{R}^2}\rho}.

若整数对称阶满足 $m>2\alpha(\rho)$ ，则方程 $\nabla\cdot(\rho\nabla\psi)=f$ 有唯一的 $$m$$ 重对称解，满足二次、一次增长界

|\psi(x)|\leq C\|f\|_{\mathring C^{0,\gamma}}|x|^2, \qquad |\nabla\psi(x)|\leq C\|f\|_{\mathring C^{0,\gamma}}|x|,

且 $\nabla^2\psi\in\mathring C^{0,\gamma}$ ，并有

\|\nabla^2\psi\|_{\mathring C^{0,\gamma}} \leq C\|f\|_{\mathring C^{0,\gamma}}.

常数依赖于 $\gamma,m,\inf\rho,\sup\rho,M$ 。这里的对称阶限制依赖密度对比度；当 $\rho$ 为常数时条件退化为 $$m>2$$ ，与常系数结果相符。

局部适定性（定理 1.5）。设 $\Omega_0\in\mathring C^{1,\gamma}_m$ ， $\rho_0\in\mathring C^{2,\gamma}_m$ ，且 $0<\inf\rho_0$ 、 $m>2\alpha(\rho_0)$ 。则存在

T=T\bigl(\gamma,m,\inf\rho_0,\sup\rho_0, \|\rho_0\|_{\mathring C^{2,\gamma}}, \|\Omega_0\|_{\mathring C^{1,\gamma}}\bigr)>0,

使系统 $$(1.12)$$ 有唯一解；对每个 $0<\beta<\gamma$ ，

\Omega\in L^\infty(0,T;\mathring C^{1,\gamma}_m) \cap C([0,T];\mathring C^{1,\beta}_m),

\rho\in L^\infty(0,T;\mathring C^{2,\gamma}_m) \cap C([0,T];\mathring C^{2,\beta}_m).

这是尺度不变 Hölder 框架中的局部理论；原文说明所列初值正则性可能并非最优。

径向齐次约化及整体适定性（定理 1.12）。若初值 $g_0\in C^{1,\gamma}(\mathbb{S}^1)$ 、 $P_0\in C^{2,\gamma}(\mathbb{S}^1)$ 且 $\inf P_0>0$ ，令 $\alpha=\sup P_0/\inf P_0$ 。作者要求： $$(g_0,P_0)$$ 为 $$m$$ 重对称， $m>4\sqrt{\alpha}$ ； $$g_0$$ 关于 $$0$$ 为奇对称， $$P_0$$ 关于 $$0$$ 为偶对称；在基本扇区 $[0,\pi/m]$ 上 $$P'_0$$ 固定符号；并且 $$P'_0$$ 的支集与扇区端点保持正距离，即存在 $\varepsilon>0$ 使

\operatorname{supp}(P'_0)\subset [-\pi/m+\varepsilon,\,\pi/m-\varepsilon].

在这些假设下，一维系统有唯一整体解；任意 $$T>0$$ 、 $0<\beta<\gamma$ 时， $$g$$ 与 $$P$$ 分别属于 $L^\infty(0,T;C^{1,\gamma}_m)\cap C([0,T];C^{1,\beta}_m)$ 和 $L^\infty(0,T;C^{2,\gamma}_m)\cap C([0,T];C^{2,\beta}_m)$ 。由于 $$P$$ 沿特征输运，密度对比度保持不变。单调性假设还给出文中明确写出的下界 $g(t,\theta)\geq-\|g_0\|_{L^\infty}$ ，适用于基本扇区及作者采用的符号约定。

回到二维方程（定理 1.13）。若初始 $\Omega_0,\rho_0$ 为上述径向齐次数据，具有定理 1.12 的结构，且 $m>\max\{2\alpha(\rho_0),4\sqrt{\alpha(\rho_0)}\}$ ，则 $$(1.12)$$ 有唯一整体解，任意有限 $$T$$ 上保持定理 1.5 所列的 $\mathring C^{1,\gamma}$ 、 $\mathring C^{2,\gamma}$ 有界正则性及较低指数 Hölder 连续性。借助动量涡量与速度的椭圆反演，这对应原始非齐次 Euler 方程的一类整体时间解。摘要还称，对某些带符号的整体解证明了“密度均匀化”与“平均弛豫”之间的二分；但所给 PDF 摘录未包含相应定理、带符号条件及收敛量的定义，故无法从提供内容中确认该二分的精确数学表述与速率。

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

证明首先要处理一个不同于标准 Biot-Savart 反演的椭圆问题。方程 $\nabla\cdot(\rho\nabla\psi)=f$ 没有显式平移不变核， $\rho$ 还可不衰减，因而不能直接把常系数奇异积分估计套用到变系数算子上。作者把问题拆成唯一性与先验估计两部分，并利用二维特有的准共形结构解决前者。

唯一性论证的关键不是简单的能量积分：在无穷远增长至多为 $$|x|^2$$ 的类中，齐次方程未必能仅凭有限能量条件排除非零解。作者把标量导电系数方程 $\nabla\cdot(\rho\nabla\psi)=0$ 与复平面上的准共形映射联系起来。其核心结论是，给定 $$m$$ 重对称的正系数 $\rho$ ，相应解在拓扑变换下对应于一个调和函数 $\phi$ ，而增长阶由密度对比度控制，可至多为 $2\alpha(\rho)$ 。对称性使调和函数的角向频率只能取与 $$m$$ 重旋转相容的频率；当 $m>2\alpha(\rho)$ 时，增长界与允许的角向频率不相容，Liouville 型结论遂迫使 $\phi$ 为常数，进而 $\psi=0$ （在相应归一化下）。因此，对称阶条件既用于消除齐次解，也明确反映系数对比度的代价。常系数时该条件回到 $$m>2$$ 。

存在性和二阶 Schauder 估计则处理尺度不变空间中的变系数算子。作者在 $\mathring C^{0,\gamma}$ 框架发展先验估计，再用紧性—反证论证排除估计常数失控的可能；该论证依赖前述唯一性，以及常系数对照问题的估计。随后通过连续性方法从常系数情形延拓至一般正系数，得到 $\|\nabla^2\psi\|_{\mathring C^{0,\gamma}}$ 对源项的控制。其作用不仅是椭圆方程本身的可解性，也为输运方程中的速度导数提供闭合估计。PDF 摘录未给出 Schauder 先验估计的逐项证明细节，故无法进一步确认其局部冻结系数或覆盖论证的具体实现。

有了椭圆反演后，作者对 $\rho$ 与 $\Omega$ 的输运型系统建立局部理论。证明框架是对线性输运方程在 $\mathring C^{k,\gamma}$ 空间中的估计与椭圆估计相结合，再通过迭代构造解。其结构性难点在于速度 $u=\nabla^\perp\psi$ 由当前密度决定，故输运速度与密度、动量涡量之间存在非线性反馈；同时 $\Omega$ 方程右端含 $\nabla|u|^2\cdot\nabla^\perp\rho$ 。椭圆估计负责控制速度的空间正则性，输运估计负责传播密度和动量涡量正则性，而较低指数 $\beta<\gamma$ 的时间连续性是定理所给出的结论形式。径向齐次性再由方程的不变性与局部解唯一性传播，从而把二维解严格联系到圆周上的约化系统。

全局部分的技术核心是控制一维系统的延拓判据，而不是把局部存在直接宣称为整体存在。椭圆边值问题 $\partial_\theta(P\partial_\theta G)+4PG=g$ 先通过 Green 函数构造并估计，以从 $$(g,P)$$ 控制 $$G$$ 。随后分析特征线与非局部源项。密度满足输运方程，故其正性和对比度保持；对称性使特征线受限于基本扇区。对 $$P'_0$$ 的单调性，作者利用输运结构和比较原理得到 $g\geq-\|g_0\|_{L^\infty}$ 。然而 $$g$$ 方程中的 $$-(4G^2+(G')^2)P'$$ 会造成额外增长，单凭密度输运并不能控制所有导数。假设 $$P'_0$$ 的支集避开扇区端点，使密度在边界邻域为常数；结合特征线被对称性限制在扇区内，这为延拓判据所需的边界控制提供条件。作者明确指出这一支集假设是当前证明的技术条件，并猜测它可能并非必要。所给摘录只说明 Green 函数估计和延拓判据的章节安排，没有呈现判据的精确公式及其逐步推导，因此不能据此补写未显示的定量增长估计。

4. 与之前工作的比较

首先，本文的局部理论与 Lions 提出的问题相关，但不能说解决了其原始问题。Lions 所问的是常密度情形下，不要求离散旋转对称、涡量有界且速度至多线性增长时的适定性。本文处理的是变密度方程，并允许线性增长速度，但依赖 $$m$$ 重对称性和尺度不变 Hölder 正则性；椭圆反演还要求 $m>2\alpha(\rho)$ 。因此本文在方程系数上更一般，却在初值几何结构上更受限。作者将此称为在对称类中的临界正则性局部理论，而不是对 Lions 无对称问题的完整回答。

与 Elgindi 和 Jeong [13] 的常密度离散对称理论相比，本文沿用 $\mathring C^{k,\gamma}$ 尺度不变空间以及对称性对无穷远反演的作用，但核心椭圆算子从 $\Delta$ 变成 $\nabla\cdot(\rho\nabla\cdot)$ 。常密度时可利用 Biot-Savart 核的对称化以及其角向抵消；变密度后没有可直接使用的显式核。本文以二维准共形映射处理唯一性，并发展变系数 Schauder 估计处理存在性。代价是对称阶不仅须大于 $$2$$ ，还须随密度对比度增大至 $m>2\alpha(\rho)$ 。因此推广的主要增量是变系数椭圆反演，而不是取消对称性限制。

在常密度二维 Euler 方面，Wolibner 与 Yudovich 的全局理论依赖 $\partial_t\omega+u\cdot\nabla\omega=0$ ，以及衰减条件下的 Biot-Savart 反演，或相应弱解框架。本文处理的是不要求远场衰减、允许速度线性增长的尺度不变情形，并使用动量涡量 $\Omega=\nabla^\perp\cdot(\rho u)$ 。这使研究对象更符合全平面流动可有非衰减远场结构的动机，但不能据此推得无对称或一般变密度数据的全局理论。本文的整体结果只适用于满足径向齐次、离散对称、奇偶性、单调性及端点支集条件的特殊数据族。

既有非齐次 Euler 局部适定性结果见文中引述的 [4]、[5]、[8]、[11]、[12]、[17]、[22]，较长时间存在工作包括 [3]、[12]；Fanelli [15] 使用过动量涡量变量来建立爆破判据。本文相较这些结果的具体函数空间、初值假设和时间范围逐项比较，所给摘录没有提供足够信息，不能概括为统一的假设放宽。可以确认的区别是：本文在特定对称尺度不变框架中建立变系数椭圆反演与局部理论，并借径向齐次约化得到该数据类的整体存在，而非仅给出一般延拓判据。

全局结果与扰动型工作 [9]、[26] 及行波工作 [25] 也有不同。[9] 在周期通道 $\mathbb{T}\times\mathbb{R}$ 上，对常密度 Couette 流的小 Gevrey 扰动证明整体适定；[26] 在有限通道上研究带密度剖面的单调剪切流附近 Gevrey 扰动；[25] 在有限通道的单调剪切流附近构造 Sobolev 正则行波。本文所称的推进在于给出非小数据的时间依赖整体解类，而不是只分析某个稳定剪切流附近的扰动或构造行波。相应代价是解类高度结构化：二维数据必须径向齐次，且约化初值还满足 $m>4\sqrt{\alpha}$ 、单调性与端点附近密度导数为零等条件。作者也指出数值观察显示去掉端点支集条件后可能出现集中增长，并猜测该条件或许非必要；这仍是猜想，不是定理。

文中还把结果与常密度径向齐次一维系统的工作 [14] 联系起来：该工作研究 $m\geq4$ 的常密度约化系统长时间行为，本文加入可变密度后，新的非局部源项使估计更困难，并得到具有密度对比度依赖的对称阶条件。关于摘要所称的带符号整体解之“密度均匀化或平均弛豫”二分，提供的 PDF 内容尚未达到相应章节，无法核验具体假设、极限对象、收敛拓扑或速率；作者对前人是否已有同类结果的判断也主要是摘要中的文献认知陈述，不能由当前摘录独立验证。

2. 具有支撑延伸至边界的涡量扰动的库埃特流的非线性无黏阻尼

Nonlinear Inviscid Damping of Couette Flow with Vorticity Perturbations Supported up to the Boundary

作者: Siming He, Sameer Iyer, Weiren Zhao

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29975v1

摘要翻译: 我们证明了有限通道中近库埃特流的非线性无黏阻尼，其中涡量扰动在边界处无穷阶消失。为研究边界效应，我们将涡量分解为支撑向边界扩展的二进层，并采用归纳式牛顿迭代方案追踪各层之间的非线性交互。最后，我们建立了精确的椭圆估计，以量化不同层所产生的速度。

Comment: 72 pages

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

1. 研究对象和背景

本文研究二维不可压 Euler 方程中有限通道内 Couette 流 $$(y,0)$$ 的非线性长期稳定性。空间区域为周期通道 $\Omega=\mathbb T_{2\pi}\times[-1,1]$ ，速度扰动记为 $$U=(U^x,U^y)$$ ，涡量扰动记为 $\omega=\partial_xU^y-\partial_yU^x$ 。作者考察的方程是

\begin{cases} \partial_t\omega+y\partial_x\omega+U\cdot\nabla\omega=0,\\ U=\nabla^\perp\Delta_D^{-1}\omega,\\ U^y(t,x,\pm1)=0,\\ \omega(0,x,y)=\omega_{\rm in}(x,y), \end{cases}

其中 $\nabla^\perp=(-\partial_y,\partial_x)$ ， $\Delta_D^{-1}$ 是通道上的 Dirichlet Laplace 逆算子。边界条件表示流体不穿透墙面；文中还采用零通量、零调和流规范，使流函数在两侧边界均取零，从而确定所用的 Biot–Savart 算子。

Couette 流的线性化具有显式混合机制。对非零水平 Fourier 模 $$k$$ ，线性涡量满足

\omega_k(t,y)=e^{-ikty}\omega_{{\rm in},k}(y).

涡量的 $$L^2$$ 大小不会因该输运而衰减，但其在垂直方向形成越来越细的尺度；通过椭圆反演恢复速度时，细尺度混合转化为速度衰减。线性问题的典型各向异性估计为

\|P_{\ne0}U^x(t)\|_{L^2}\lesssim \langle t\rangle^{-1}\|P_{\ne0}\omega_{\rm in}\|_{L_x^2H_y^1},\qquad \|U^y(t)\|_{L^2}\lesssim \langle t\rangle^{-2}\|P_{\ne0}\omega_{\rm in}\|_{L_x^2H_y^2}.

这种在无黏性条件下由相位混合和椭圆平滑产生的速度衰减称为无黏阻尼。非线性问题的困难在于不同时间、频率之间的共振相互作用可能部分逆转混合，形成 echo 链并导致正则性损失。因此，在无界通道上建立非线性阻尼的工作使用了精细的时间频率权重和较强的 Gevrey 正则性。已有有限通道结果则通常要求扰动与边界保持固定正距离，使得边界附近不必参与非线性估计。

这一支撑条件有实质性的边界原因。对一般剪切流，边界处的剪切曲率与涡量迹可能引发随时间增长的边界效应。文中回顾的 Zillinger 结果表明，在某些单调剪切背景下，若边界 $$b$$ 满足 $u''(b)\omega_{{\rm in},k}(b)\ne0$ ，则沿背景剪切输运的涡量模不可能具有一致有界的 $$H^2$$ 范数；其边界法向导数会出现对数增长。精确 Couette 剪切满足 $$u''=0$$ ，但非线性演化的平均流会改变有效背景剪切，故不能仅凭线性模型排除边界困难。

本文的问题是：能否将有限通道既有阻尼定理中的“扰动紧支撑于远离边界的区域”放宽为“扰动支撑可延伸至边界”，同时仍获得全局非线性无黏阻尼？作者并不只假设定性边界平坦，而要求涡量在边界无限阶消失，并给出定量加权控制。令 $$d(y)=1-|y|$$ ，假设初值属于相应 Gevrey 空间且满足

\|\omega_{\rm in}\|_{G^s_{3\lambda_0;2\sigma+8}} +\|e^{c_0d(y)^{-\alpha_*}}\omega_{\rm in}\|_{L^2(\Omega)}\le\varepsilon, \qquad \frac12<s<1,\quad \alpha_*>\frac{s}{1-s}.

该指数权重意味着距边界约 $2^{-n}$ 的薄层内，初始涡量的 $$L^2$$ 质量至多为 $\varepsilon e^{-c2^{n\alpha_*}}$ 。本文试图弥补的空白，正是以这一可量化的无限阶边界消失代替固定正距离支撑，并在存在边界效应的有限通道中证明非线性阻尼。论文还强调，单纯的定性平坦性不足以表达其估计所需的逐尺度小量。

2. 主要定理或主要结果

主要结论为定理 1.1。设 $$1/2<s<1$$ 、 $\alpha_*>s/(1-s)$ ，且实值初始涡量边界平坦，即 $\partial_y^j\omega_{\rm in}(x,\pm1)=0$ 对每个 $j\ge0$ 成立。存在 $\sigma_0=\sigma_0(s,\alpha_*)$ ，使得对所有 $\sigma\ge\sigma_0$ ，固定 $\lambda_0,c_0>0$ 后，存在 $\varepsilon_0>0$ 、 $\lambda_\infty\in(0,\lambda_0/2)$ 及 $c_\infty\in(0,c_0)$ ；若 $0<\varepsilon\le\varepsilon_0$ 且

\|\omega_{\rm in}\|_{G^s_{3\lambda_0;2\sigma+8}} +\|e^{c_0d(y)^{-\alpha_*}}\omega_{\rm in}\|_{L^2(\Omega)}\le\varepsilon,

则方程有全局解，并满足下列结论。定理中的常数可依赖于 $s,\alpha_*,\sigma,\lambda_0,c_0$ ，但不依赖于 $$t$$ 和 $\varepsilon$ 。

第一，边界加权涡量范数一致受控：

\sup_{t\ge0}\|e^{c_\infty d(y)^{-\alpha_*}}\omega(t)\|_{L^2(\Omega)}\le C\varepsilon.

因此，初始数据所具有的定量边界消失在非线性演化下仍可一致传播，但指数常数由 $$c_0$$ 降为较小的 $c_\infty$ 。

第二，涡量在随剪切移动的坐标中保持 Gevrey 控制并散射。令

\Phi(t,y)=\int_0^t\langle U^x\rangle(\tau,y)\,d\tau,

其中 $\langle U^x\rangle=P_0U^x$ 是水平平均。存在 $\omega_\infty\in G^s_{\lambda_\infty;\sigma-6}$ ，使得

\sup_{t\ge1}\|\omega(t,x+ty+\Phi(t,y),y)\|_{G^s_{\lambda_\infty;\sigma}}\le C\varepsilon,

且

\lim_{t\to\infty} \|\omega(t,x+ty+\Phi(t,y),y)-\omega_\infty(x,y)\|_{G^s_{\lambda_\infty;\sigma-6}}=0.

这里散射是在由 Couette 平移和演化平均流共同确定的坐标中成立；结论容许 Gevrey 指数参数从初始条件中的 $3\lambda_0$ 降至 $\lambda_\infty<\lambda_0/2$ ，并在收敛范数中损失六阶多项式权重。

第三，非零水平模的水平速度与垂直速度分别满足

\|P_{\ne0}U^x(t)\|_{L^2(\Omega)}\le C\varepsilon\langle t\rangle^{-1},\qquad \|U^y(t)\|_{L^2(\Omega)}\le C\varepsilon\langle t\rangle^{-2}.

作为推论，非零水平模涡量满足 $P_{\ne0}\omega(t)\rightharpoonup0$ 于 $L^2(\Omega)$ 的弱收敛。定理没有声称涡量强收敛至零；其散射结论是在移动坐标下收敛至一般非零的极限剖面。

参数条件方面，作者明确指出 $\alpha_*>s/(1-s)$ 是当前证明使用的充分条件，并非已知最优条件。文中讨论了通过更仔细计算将其改进至 $\alpha_*\ge s/(1-s)$ 的可能性，也提到使用此前提出的 Fourier 乘子后或可覆盖 $$s=1/2$$ ；这些属于作者提出的改进方向，并非定理 1.1 已证明的范围。对有限阶边界消失是否会导致非线性不稳定，文中也未给出结论。

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

证明的核心困难是边界附近的局部化会使椭圆反演估计随薄层变窄而恶化，而非线性 echo 控制又要求在 Gevrey 尺度上估计频率相互作用。作者的策略不是把全通道问题直接视为一个无边界问题，而是按距边界的尺度分解涡量，分别追踪每层的演化和其产生的速度，再利用初始边界权重提供的极小量补偿局部椭圆估计中的损失。

首先，作者构造沿边界递进的 dyadic 层分解。令 $$d(v)$$ 为变换后区间端点的距离， $d_m=2^{-m}$ ，以 Gevrey- $\rho$ 截断函数构造壳层分割 $\{\Psi_m\}$ ，其支集满足 $d(v)\simeq2^{-m}$ 。对初始数据，这一分解配合边界加权假设给出第 $$n$$ 层增量的估计

\|f_{\rm in}^{(n)}\|_{G^s_{\lambda_0;\sigma+2}} +\|h_{\rm in}^{(n)}\|_{G^s_{\lambda_0;\sigma+2}} +\|g_{\rm in}^{(n)}\|_{G^s_{\lambda_0;\sigma+2}} \le C\varepsilon e^{-c_{\rm in}2^{n\alpha_*}}.

其中 $$f,h,g$$ 是移动坐标中的涡量剖面、坐标变形和平均流变量。有限截断层数 $$N$$ 下先建立估计，再令 $N\to\infty$ ；估计常数须与 $$N$$ 无关。边界层虽无限多，但其初始增量具有双指数型可求和衰减。

其次，对层增量采用 Newton 型迭代，而不是对每一层单独求解原方程后再粗略相加。一般模型中，若 $F^{[n]}=\sum_{m\le n}F^{(m)}$ ，增量方程安排为

(\partial_t+y\partial_x-E(F^{[n]})\cdot\nabla)F^{(n)} =E(F^{(n)})\cdot\nabla F^{[n-1]}.

由于 $$E$$ 对 $$F$$ 线性，左右两项重组后使每个部分和 $F^{[n]}$ 精确满足相应的截断非线性方程。这样既保留每一增量的局部支撑信息，又把增量之间的非线性耦合组织成可迭代控制的形式；右端涉及较早的部分和，估计时利用稍大的 Gevrey 半径吸收导数损失。对 Euler 方程，还须同时处理依赖解的椭圆算子，因此作者在椭圆方程层面取精确差值，而不把不同迭代层的椭圆系数误当成相同。

第三，证明在长时间区域采用随平均剪切移动的 $$(z,v)$$ 坐标。其定义为

z=x-tv(t,y),\qquad v(t,y)=y+\frac1t\int_0^t\langle U^x\rangle(\tau,y)\,d\tau.

若 $f,\varphi,g,h$ 分别表示新坐标中的涡量、流函数、平均速度变量及坐标变形，则剖面方程包含

\partial_t f=(1+h)(\partial_v\varphi_{\ne0}\partial_zf-\partial_z\varphi_{\ne0}\partial_vf)-g\partial_vf,

以及平均变量方程和带时变系数的椭圆方程

\Delta_t\varphi=f,\qquad \Delta_t=\partial_z^2+(1+h)^2(\partial_v-t\partial_z)^2 +(1+h)(\partial_vh)(\partial_v-t\partial_z).

这组坐标吸收主导的 Couette 混合与平均剪切漂移，使非线性相互作用可以用 profile Fourier 变量分析；同时，椭圆算子本身会随 $$h$$ 改变，故其差分估计是证明中的一项实质工作。文中还利用边界平均速度的守恒，使变换后的区间端点固定。

第四，边界局部椭圆估计需要另一种 Gevrey 指数 $\rho$ ，满足 $s<\rho<1$ 。尺度为 $2^{-m}$ 的 Gevrey- $\rho$ 截断核提供频率差的衰减 $e^{-\lambda_m^\Delta|\eta-\xi|^\rho}$ ，其中 $\lambda_m^\Delta=c_\Delta2^{-m\rho}$ 。将 Gevrey- $$s$$ 权重从输出频率移至输入频率时会产生至多 $e^{C|\eta-\xi|^s}$ 的代价。优化两种指数后出现局部化损失，量级为

e^{C2^{m\alpha^\dagger}},\qquad \alpha^\dagger=\frac{\rho s}{\rho-s}.

作者选择 $\rho$ 使 $\alpha^\dagger<\alpha_*$ 。于是第 $$n$$ 层初始质量 $e^{-c2^{n\alpha_*}}$ 足以压过所有 $m\le n$ 的椭圆损失 $e^{C2^{m\alpha^\dagger}}$ ，并留下对 $$n$$ 可求和的衰减。这个指数竞争解释了定理对边界消失阶的要求，是层分解方案得以闭合的关键，而非仅是技术性的光滑性假设。

第五，控制各层的支撑还依赖垂直速度的时间衰减。作者在证明框架中使用类似

|U^y(t,x,y)|\lesssim\frac{\varepsilon(1-|y|)}{(1+t)^2}

的边界相对估计。沿特征线积分 $\partial_t\log(1-|y(t)|)$ 后，距墙距离至多损失一个与时间无关的因子，因此每个有限层的支撑不会在有限时间内撞到边界。这里的相对距离因子十分重要：仅有 $$U^y$$ 的绝对衰减不足以控制趋近边界的任意薄层。论文指出，Euler 情形还需通过精确椭圆估计验证不同层产生的速度具有所需控制。

最后，全文将这些局部椭圆估计与 reaction/transport multiplier 估计、能量恒等式结合，以统一吸收层间非线性相互作用和坐标变形项。目录列出了相关估计与闭合步骤，但所提供摘录未包含第 3 至第 7 节的具体乘子定义、完整引理陈述及逐项能量不等式。因此，不能仅依据摘录进一步确认每一种 echo 共振项的精确乘子结构或各误差项的具体时间可积指数；可以确认的是，证明以 Gevrey 时间频率控制、层依赖椭圆估计和 Newton 增量能量闭合为主要框架。

4. 与之前工作的比较

本文相对于既有结果最直接的推进，是将有限通道 Couette 流的非线性无黏阻尼从“初始扰动紧支撑于远离墙面的固定区域”推进到“支撑可以延伸至墙面，但涡量必须在墙面无限阶消失并满足定量指数权重”。Ionescu–Jia 将 Bedrossian–Masmoudi 的 Couette 非线性阻尼框架扩展到有限通道时要求扰动与边界保持固定正距离。对稳定单调剪切流，Ionescu–Jia 与 Masmoudi–Zhao 的结果同样采用相关量支撑远离边界的条件；变密度有限通道结果也要求剪切、背景密度和扰动满足类似紧支撑条件。本文的层分解正面处理这些结果通过支撑假设排除掉的边界薄层。

与无界通道上的 Bedrossian–Masmoudi 工作相比，本文面对的新增困难不是单纯改变空间区间，而是 Dirichlet 椭圆反演与边界薄层尺度之间的相互作用。无界域中建立非线性阻尼的理论通过适当 Gevrey 正则性与不平衡 Fourier 乘子控制 echo 链；本文仍处于 Gevrey 框架，但要额外为不同边界距离尺度建立局部化椭圆估计。层越靠近墙面，局部截断的 Gevrey 半径越小，椭圆核的频率衰减越弱，因而会产生随层数增长的损失。作者以 $e^{c_0d^{-\alpha_*}}$ 形式的边界小量支付这一损失。这一做法换取了允许支撑接触边界，但代价是显著强于一般 Sobolev 小数据或仅有定性边界平坦性的初值条件。

本文也与 Zillinger 的边界驱动障碍形成互补。Zillinger 指出，对一般平滑剪切流，若 $u''(b)\omega_{{\rm in},k}(b)\ne0$ ，沿背景剪切输运的涡量在边界附近会出现导数增长，因而不能期待相应的统一 $$H^2$$ 控制；相关非线性障碍还需要附加的平均剪切非退化性及边界非线性项衰减条件。本文研究的是精确 Couette 参考流附近，并要求初始涡量在墙面无限阶消失，因此直接排除了上述边界涡量迹非零的情形。它证明在这一强边界消失条件下，非线性边界效应仍可被控制；但没有消除 Zillinger 障碍所针对的一般非零边界迹情形，也没有证明有限阶消失足够。

在参数范围上，定理要求 $$1/2<s<1$$ 及 $\alpha_*>s/(1-s)$ ，并假设较高的 Gevrey 正则性和足够大的多项式权重指数 $\sigma$ 。这是为边界尺度损失与 Gevrey- $$s$$ 权重之间的竞争而付出的代价。文中明确将 $\alpha_*=s/(1-s)$ 及 $$s=1/2$$ 列为可能改进方向，而不是现有结论。相较于无界域或远离边界的有限域结果，本文扩大的是空间支撑的允许范围，不是降低正则性门槛；对于有限阶边界消失，作者指出其是否引起不稳定仍为开放问题，并给出线性层面的提示，但未将其证明为非线性不稳定定理。

衰减结论方面，本文获得 $\|P_{\ne0}U^x(t)\|_{L^2}\lesssim\varepsilon\langle t\rangle^{-1}$ 与 $\|U^y(t)\|_{L^2}\lesssim\varepsilon\langle t\rangle^{-2}$ ，并在随剪切和平流平均量移动的坐标中证明涡量散射。摘要和定理并未声称这些速度衰减率相较既有 Couette 结果达到更优时间阶；主要新意是这些结论在边界可被数据触及的条件下仍成立。至于同既有结果相比最优常数、完整正则性阈值或各项误差估计的精细强弱，提供的摘录没有给出足够信息作进一步判定。

作者还在备注中提出，希望将结论推广到满足 Ionescu–Jia 与 Masmoudi–Zhao 假设的一类稳定单调剪切流，但指出仍需建立相应的 wave-operator 估计和非线性闭合；该推广并未在本文证明。特别是既有工作中关于 $$u''$$ 支撑远离边界的条件，作者认为只有在明确给出定量边界假设并证明相应算子估计后才能放宽。因此，本文结论目前应理解为精确 Couette 背景下的边界层突破，而非一般单调剪切流、非单调剪切流或变密度模型的边界版本。

3. PDE理论引导的神经算子与二维各向异性 Boussinesq 系统的长期弛豫

A PDE-theory-guided neural operator and the long-time relaxation of the two-dimensional anisotropic Boussinesq system

作者: Akram Moustafa, Jiahong Wu, Xiaoming Zheng

类别: math.AP, math.NA

发布日期: 2026-10-21

arXiv ID: 2610.29938v1

摘要翻译: 我们基于对二维周期盒上具有各向异性耗散的 Boussinesq 系统长期行为的严格分析，设计了一种神经算子。该算子将振幅为 $\varepsilon$ 的静水平衡小扰动映射为晚期时刻的涡量。我们证明，流动会弛豫至分层剪切流，其中水平平均涡量与其初始值 $\barω_0$ 的差保持在 $O(\varepsilon^2)$ 以内，而相对于该平均值的振荡量级为 $\varepsilon$ 。因此，网络的输出由 $\barω_0$ 与这两部分组成，并分别按其大小进行缩放；这两部分的大小分别增长为 $\varepsilon^2$ 和 $\varepsilon$ 。若不采用这种引导，训练要么无法收敛，要么无法捕捉 $\varepsilon^2$ 的增长。在这种引导下，同一个网络能够在 $2\times10^3$ 次测试运行中逐样本复现该量；这些测试运行的振幅跨度为一个因子 $$23$$ ，网络的中位误差为 $$2.2$$ 百分比。由于该网络速度快且关于初始数据可微，我们还用它搜索初始数据空间；对每个候选初始数据均由求解器进行检验后，搜索找到了使该量达到测试运行中最大值 $$4.4$$ 倍的初始数据，并将定理中常数的下界提高了相应倍数。

Comment: 无

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

1. 研究对象和背景

本文研究二维周期区域上的部分耗散 Boussinesq 系统在静水压平衡附近的小扰动演化，以及如何将其严格长时间分析用于构造神经算子。模型的速度方程仅有水平方向黏性，温度方程仅有竖直方向扩散；各向异性使两类耗散都存在未直接受耗散控制的模态。令静水压平衡为 $(U_{he},P_{he},\Theta_{he})=(0,\tfrac12 x_2^2,x_2)$ ，并以 $(u,p,\theta)$ 表示相对于该平衡的扰动，则在 $\mathbb T^2=[0,2\pi]^2$ 上考虑

\begin{aligned} \partial_t u+u\cdot\nabla u&=-\nabla p+\nu\partial_{11}u+\theta e_2,\\ \partial_t\theta+u\cdot\nabla\theta+u_2&=\eta\partial_{22}\theta,\\ \nabla\cdot u&=0. \end{aligned}

其中 $$u=(u_1,u_2)$$ 是速度扰动， $\theta$ 是温度扰动， $\nu,\eta>0$ 分别为水平黏性和竖直热扩散系数。以涡量 $\omega=\partial_1u_2-\partial_2u_1$ 和流函数 $\psi$ 改写，得到

\partial_t\omega+u\cdot\nabla\omega=\nu\partial_{11}\omega+\partial_1\theta, \qquad -\Delta\psi=\omega, \qquad u=(\partial_2\psi,-\partial_1\psi),

并仍有 $\partial_t\theta+u\cdot\nabla\theta+u_2=\eta\partial_{22}\theta$ 。研究的关键未知量不仅是瞬时解，也包括其长时间极限，尤其是涡量最终留下的分层剪切结构，以及该极限如何依赖初始速度和温度。

小数据假设为散度自由初速度及零均值初始扰动满足

\int_{\mathbb T^2}\theta_0\,dx=\int_{\mathbb T^2}u_{2,0}\,dx=\int_{\mathbb T^2}u_{1,0}\,dx=0, \qquad \|u_0\|_{H^2}+\|\theta_0\|_{H^2}\leq\varepsilon.

零均值条件并非纯技术约定：全局均值中的 $$u_2$$ 与 $\theta$ 满足一个振荡耦合关系，不会一般地衰减；要求这两个均值初始为零，才排除该有限维不衰减分量。 $$u_1$$ 的全局均值则由流动守恒。本文使用沿 $$x_1$$ 的平均与振荡分解

\bar g(x_2,t)=\frac1{2\pi}\int_0^{2\pi}g(x_1,x_2,t)\,dx_1, \qquad \widetilde g=g-\bar g.

水平黏性作用于 $$x_1$$ 依赖的振荡模态，却湮灭于水平平均；竖直热扩散则无法直接消除与 $$x_2$$ 无关的温度模态。因此，常规依赖全梯度耗散和 Poincaré 不等式的衰减论证不能直接闭合。相较于全耗散情形中扰动整体趋于零，这种各向异性系统可能保留仅依赖 $$x_2$$ 的涡量，即分层剪切流。

这一问题建立在部分耗散 Boussinesq 系统小数据全局正则性与稳定性理论之上。文中指出，条带区域上的相关稳定性结果见文献 [1]，其他部分耗散配置见 [14]，无热扩散的浮力驱动情形见 [30,6]；但这些稳定性结论本身并未回答解趋向何种长时间状态，也未量化极限由初值决定的部分。周期盒上的一致 $$H^2$$ 稳定性近期由 [7] 建立，而条带上的能量方法不能直接移植到周期盒，须使用适合该区域的稳定性机制。本文在此基础上分析长时间极限，并将所得的幅度分层结构作为神经算子的输出设计依据。

应用方面，神经算子可快速近似由初值到晚时刻解的映射，但直接预测完整涡量会遇到尺度严重不平衡：初始水平平均 $\bar\omega_0$ 通常远大于其长时间变化，而平均变化仅为二阶小量；涡量振荡则是一阶量。网络若将完整末态作为单一目标，容易把需要检验的微小平均变化淹没在整体预测误差中。此外，温度的晚时刻解对初值具有强烈的模态依赖和病态敏感性。本文试图弥补的空白，正是把定理所揭示的近似守恒量、不同分量的幅度阶以及温度的线性衰减特征，具体转化为网络输出参数化，并进一步用快速代理模型探索分析估计常数的下界。

2. 主要定理或主要结果

严格分析的核心结果是定理 3.1。假设初值满足上述零均值、散度自由和 $$H^2$$ 小数据条件，且 $\varepsilon\leq\varepsilon_*(\nu,\eta)$ ，则依据附录 A 中引用自 [7] 的全局稳定性结果，系统有唯一全局解，并且

\sup_{t\geq0}\bigl(\|u(t)\|_{H^2}+\|\theta(t)\|_{H^2}\bigr)\leq C_*\varepsilon, \qquad \int_0^\infty\left(\nu\|\partial_1u\|_{H^2}^2+\eta\|\partial_2\theta\|_{H^2}^2\right)dt\leq C_\nu\varepsilon^2.

这些估计提供长时间分析所需的全局控制。零均值条件沿流保持。

定理 3.1(a) 断言，竖直速度和温度在任意低于 $$H^2$$ 的 Sobolev 正则度中趋于零：对每个 $$s<2$$ ，

u u u_2(t)\longrightarrow0,\qquad \theta(t)\longrightarrow0 \quad\text{于 }H^s(\mathbb T^2),\qquad t\to\infty.

此处公式中的首项为 $$u_2(t)$$ ，即结论是 $u_2(t)\to0$ ；上述排版若将其读作其他量均属笔误。原定理没有给出该非线性衰减的统一时间速率。

定理 3.1(b) 给出涡量的分层极限。存在 $\bar\omega_\infty=\bar\omega_\infty(x_2)\in H^1(\mathbb T^2)$ ，使得对每个 $$s<1$$ ，

\omega(\cdot,t)\longrightarrow\bar\omega_\infty \quad\text{于 }H^s(\mathbb T^2),\qquad t\to\infty.

极限由初始水平平均涡量加上沿整个演化轨道累积的非线性修正给出：

\bar\omega_\infty=\bar\omega_0+\int_0^\infty\partial_{22}\bigl(\widetilde u_1\widetilde u_2\bigr)(\cdot,t)\,dt,

其中积分在 $H^{-2}$ 中收敛。并且存在仅依赖于 $\nu$ 的常数 $$C$$ ，使

\sup_{t\geq0}\|\bar\omega(\cdot,t)-\bar\omega_0\|_{H^{-2}}\leq C\varepsilon^2, \qquad \|\bar\omega_\infty-\bar\omega_0\|_{H^{-2}}\leq C\varepsilon^2.

关键演化恒等式为

\partial_t\bar\omega=\partial_{22}\bigl(\widetilde u_1\widetilde u_2\bigr).

因此，平均涡量并非严格守恒；其变化由振荡速度产生的二次 Reynolds 应力驱动。文中说明振荡部分具有 $O(\varepsilon)$ 的大小并在时间上平方可积，从而二次应力产生 $O(\varepsilon^2)$ 的平均漂移。对于固定 $$u_0$$ ，不同可 admissible 温度初值产生的极限之差至多为 $2C\varepsilon^2$ （在 $H^{-2}$ 中）；这表明温度初值对极限剖面的影响只在二阶量级上受控。定理没有给出 $\bar\omega_\infty$ 关于初值的显式闭式公式，也没有给出非线性收敛速率。

命题 3.2 分析线性化后的模态衰减。对 $k=(k_1,k_2)\neq0$ ，傅里叶系数满足

\frac{d\widehat u_2}{dt}=-\nu k_1^2\widehat u_2+\frac{k_1^2}{|k|^2}\widehat\theta, \qquad \frac{d\widehat\theta}{dt}=-\widehat u_2-\eta k_2^2\widehat\theta.

若 $k_1k_2\neq0$ ，则

|\widehat u_2(t)|^2+\frac{k_1^2}{|k|^2}|\widehat\theta(t)|^2 \leq e^{-2\min\{\nu,\eta\}t} \left(|\widehat u_2(0)|^2+\frac{k_1^2}{|k|^2}|\widehat\theta(0)|^2\right).

若 $k_1=0,k_2\neq0$ ，则 $\widehat u_2\equiv0$ 且 $\widehat\theta(t)=e^{-\eta k_2^2t}\widehat\theta(0)$ 。特殊模态 $$k_2=0$$ 则呈现慢衰减：当 $\nu k_1^2>2$ 时，温度系数是两个指数模态之和，其中慢特征值为

\lambda_+=\frac{-2}{\nu k_1^2+\sqrt{\nu^2k_1^4-4}}, \qquad \frac1{\nu k_1^2}\leq|\lambda_+|\leq\frac2{\nu k_1^2}<1.

故 $|k_1|\to\infty$ 时 $\lambda_+\to0$ ，线性温度衰减不存在对所有模态统一适用的指数速率。该结论解释了为何晚时刻温度虽整体趋零，却可能包含对高水平频率初值极敏感的成分。

神经算子部分是数值实验结果，而非上述解析定理的推论。作者用 Fourier 神经算子学习从 $(\omega_0,\theta_0)$ 到 $$T=100$$ 时涡量的映射，取 $\nu=\eta=1$ ，使用 $128\times128$ 伪谱计算及 $$10^4$$ 组训练样本。网络不直接预测完整末态，而以已知的 $\bar\omega_0$ 为基底，分别学习并按其预期幅度缩放水平平均变化 $\bar\omega(T)-\bar\omega_0$ 与振荡部分，二者分别对应 $\varepsilon^2$ 和 $\varepsilon$ 尺度。摘要报告，该结构在振幅跨度约 $$23$$ 倍的 $2\times10^3$ 组测试运行中，对定理所控制的量取得 $2.2\%$ 的中位误差；未采用该分析指导的训练不收敛或不能捕捉 $\varepsilon^2$ 增长。作者还使用网络对初值作快速、可微的搜索，并以原数值求解器逐一复核候选，找到使该量达到测试样本最大值 $$4.4$$ 倍的初值，从而将定理常数的数值下界提高同一倍数。这是对常数下界的计算性改进，不是对定理上界或解析常数的证明改进。

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

证明的结构主线是先从周期盒上的小数据全局稳定性获得统一控制，再利用水平平均/振荡分解识别耗散可达的分量与可能残留的分量。提供的文本说明，稳定性输入是附录 A 中定理 A.1（引自 [7]）：全局 $$H^2$$ 有界和耗散积分有限。本文据此研究长时间行为。由于附录 A、B 的具体推导并未包含在所提供的 PDF 摘录中，关于稳定性估计如何闭合的逐步能量恒等式、定理 3.1 的完整技术引理及其精细证明，从提供内容中无法确认；不能据此臆测附录中的具体加权能量或交换子估计。

长时间平均的处理则在正文中明确给出。对速度第一分量取 $$x_1$$ 平均，压力梯度项和 $\partial_{11}u_1$ 项消失；利用不可压条件将平流项写成散度，并利用全局均值条件推出 $\bar u_2=0$ 。于是

\partial_t\bar u_1=-\partial_2\overline{u_1u_2} =-\partial_2(\widetilde u_1\widetilde u_2), \qquad \bar\omega=-\partial_2\bar u_1,

从而得到 $\partial_t\bar\omega=\partial_{22}(\widetilde u_1\widetilde u_2)$ 。这一步把长时间极限问题转化为控制振荡速度所产生的二次 Reynolds 应力。平均部分不受水平黏性直接耗散，因而不能期望其像振荡部分一样直接衰减；可证明的结论是其总漂移在 $H^{-2}$ 中为二阶小量。双重导数与负阶 Sobolev 范数相配，使应力积分能够在弱拓扑中控制平均涡量变化。应力是振荡速度的乘积，故其二次结构解释了 $O(\varepsilon^2)$ 而非 $O(\varepsilon)$ 的平均漂移。关于乘积估计和时间积分的具体常数推导，须参见未提供的附录 B，无法从摘录进一步核实。

另一个障碍是证明振荡部分及温度趋于零，但各向异性耗散不能单独控制所有频率。水平平均/振荡分解与竖直方向的对应投影，分别隔离 $\partial_1$ 黏性和 $\partial_2$ 热扩散能够控制的模态；正文列出水平与竖直 Poincaré 不等式，并利用不可压条件将 $$u_2$$ 归入水平振荡。在线性化问题中，作者对每个 Fourier 模态写出 $(\widehat u_2,\widehat\theta)$ 的耦合常微分系统。对 $k_1k_2\neq0$ ，适当加权的二次量

|\widehat u_2|^2+\frac{k_1^2}{|k|^2}|\widehat\theta|^2

具有至少 $e^{-2\min\{\nu,\eta\}t}$ 的衰减；而 $$k_2=0$$ 时，特征值显式展示了 $|k_1|^{-2}$ 级的慢衰减。这种逐模态谱分析既补足了长时间行为的解释，也直接揭示晚时刻温度预测任务的病态性。非线性结论只有趋零而无统一衰减率，与该无统一线性速率现象相容；作者没有把线性速率误当作非线性定理。

神经算子的设计创新不是简单把方程残差加入损失，而是把定理给出的“哪些量已知、哪些量小、各小量如何随 $\varepsilon$ 缩放”写进输出表示。由于 $\bar\omega_0$ 是已知输入，目标改写为初始平均加两种不同尺度的修正：平均修正除以 $\varepsilon^2$ ，振荡修正除以 $\varepsilon$ 。这样避免网络在预测完整涡量时被较大的背景平均掩盖二阶小量，也减少训练集中不同振幅目标量级差别过大的问题。作者选择只输出涡量，不预测晚时刻温度：命题 3.2 表明 $$k_2=0$$ 的温度高水平频率慢模态衰减率约为 $1/(\nu k_1^2)$ ，导致从初始数据到固定晚时刻温度的映射严重病态；数值实验还显示该温度主要集中在这些模态。训练该目标未收敛的观察是实验事实，不应与理论不可能性结论混为一谈。

实现上，网络为四层 Fourier 神经算子，隐藏层在 Fourier 空间保留低频模态并学习复权重；训练数据由 $128\times128$ 伪谱离散和二阶半隐式 BDF 时间推进产生，扩散隐式处理、非线性与耦合项显式外推，积分至 $$T=100$$ 。训练样本在十个振幅上生成，并以相同初始涡量配不同温度的方式形成样本组；测试涵盖训练振幅及其间的插值振幅。作者还比较无理论结构的网络，并让代理模型对初值可微，以便进行快速搜索；每个候选最终由原求解器复核。这种“网络提出候选、可信求解器确认”的流程避免把代理模型预测误差直接当成新解析结论。具体优化器、完整损失定义、网络超参数和搜索算法细节在提供的摘录中没有全部出现，相关细节从提供内容中无法确认。

4. 与之前工作的比较

与部分耗散 Boussinesq 稳定性文献相比，本文并非首次证明小扰动全局存在。作者明确将周期盒上的一致 $$H^2$$ 稳定性归于 [7]，并将其作为长时间分析的前提；[1] 处理条带 $\mathbb T\times\mathbb R$ 上水平速度耗散、竖直温度扩散的相关系统，[14] 研究相关部分耗散配置，[30,6] 研究无热扩散的浮力驱动情形。本文推进之处在于从“解始终存在且有界”进一步识别其周期盒长时间状态： $$u_2$$ 和 $\theta$ 趋零，涡量趋向仅依赖 $$x_2$$ 的分层极限，并以 $H^{-2}$ 定量控制平均涡量相对初始值的 $O(\varepsilon^2)$ 漂移。

与 [30]（及 [6]）的比较尤其清楚。文中所述的 [30] 研究二维周期盒上速度具有全各向同性黏性、温度无扩散的系统，结果是速度衰减而温度趋向显式的 $$x_2$$ 函数。本文把速度耗散改为仅水平方向，并加入竖直热扩散，残留结构的角色随之反转：温度趋零，涡量趋向 $$x_2$$ 分层剖面。这里的极限剖面不是已知的初值显式函数，而由

\bar\omega_\infty-\bar\omega_0 =\int_0^\infty\partial_{22}(\widetilde u_1\widetilde u_2)\,dt

隐式依赖整个非线性轨道；本文提供的是其二阶小量控制，而非闭式表达式。与两种场均有全 Laplace 耗散的标准情形相比，后者可以借助耦合能量恒等式及零均值 Poincaré 不等式得到整体指数衰减；本系统的耗散核非平凡，因而出现持久的分层涡量结构，也不能沿用那种直接的整体衰减结论。

在假设和适用范围上，定理针对 $\mathbb T^2$ 、 $$H^2$$ 小初值、指定的零均值条件及本文给定的各向异性配置。小量阈值依赖 $\nu,\eta$ ，而定理 3.1(b) 中漂移常数据文中陈述仅依赖 $\nu$ 。对温度初值的比较也有明确但有限的含义：固定 $$u_0$$ ，不同 admissible $\theta_0$ 所产生的极限差在 $H^{-2}$ 中至多为 $2C\varepsilon^2$ ；这并不表示极限与温度无关，也不提供其一阶或更高阶精确依赖公式。结论中的收敛拓扑是 $u_2,\theta$ 的 $$H^s$$ （ $$s<2$$ ）及涡量的 $$H^s$$ （ $$s<1$$ ）；摘录未给出更强拓扑、定量非线性收敛速率或非小扰动结果。因此不能将线性化中部分模态的指数率提升为整个非线性解的统一衰减率。

与数值计算工作相比，文中提到 [2] 对全耗散 Boussinesq 系统发展并分析二阶数值格式，其长时间模拟显示扰动指数趋零。本文的数值求解器服务于数据生成与定理量检验，研究重点不是提出新的数值收敛定理；它考察的是部分耗散下理论预测的分层极限和小量尺度。数值设置使用 $\nu=\eta=1$ 、有限分辨率和有限终止时间 $$T=100$$ ，所以实验结果不能替代解析的 $t\to\infty$ 结论，也不能据此断言离散误差已对所有振幅和频率得到理论控制。

与已有神经算子工作的差别在于理论信息进入模型的方式。文中回顾 Fourier 神经算子、DeepONet、流体预测应用、物理信息网络、守恒律结构化网络、渐近保持模型和边界层修正网络等工作，并指出许多方法利用方程、形式极限或守恒结构；本文突出的是将严格长时间分析所给出的近似守恒量与不同分量的振幅阶直接用于输出分解和归一化，而不只是将 PDE 残差加入损失。实验表明，不采用此分解的训练未能收敛或漏掉 $\varepsilon^2$ 尺度；采用后，在摘要报告的 $2\times10^3$ 测试运行及约 $$23$$ 倍振幅跨度上，对所考察定理量的中位误差为 $2.2\%$ 。这是所报告数据分布、离散求解器和训练设置下的经验精度，不是对任意初值集合的算子逼近定理。

最后，本文把代理模型用于常数下界的数值搜索，与流体稳定性问题中有限时间增长优化的既有思路相呼应；文中点名 [25,10]，并提及机器学习辅助构造数学猜想或反例及流体奇异结构候选的工作 [5,35,28,38,37]。其可信性保障在于每个搜索候选都重新交由数值求解器检查，而非接受网络输出作为证据。搜索找到的目标量是测试集最大值的 $$4.4$$ 倍，因此把定理常数的数值下界提高同一倍数；它没有确定最佳常数，也没有改善解析上界。相应代价与局限包括对训练分布、有限网格与有限时间的依赖，以及候选验证仍需调用求解器。文中未在所提供摘录中给出该搜索的完整算法、误差认证或离散误差严谨界，因此这些方面从提供内容中无法确认；作者亦未显式建立该神经算子在连续函数空间上的泛化保证。

4. 黏性 Gröbli 三角形的长时间动力学

Long-time dynamics of viscous Gröbli triangles

作者: Michele Dolce, Antonio Tirotta

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.30038v1

摘要翻译: 1877 年，Gröbli 发现，平面上的三个点涡可以作自相似运动，形成一个旋转的三角形；该三角形要么无限扩张，要么在有限时间内塌缩。我们研究这些构型的黏性对应情形，即以排列成 Gröbli 三角形的三个点涡为初值的二维 Navier–Stokes 演化。在扩张情形中，黏性使涡核扩散的速率与三角形增长的速率相同，因此相对于自身大小而言，各涡始终彼此充分分离。基于这一观察，我们证明，在高雷诺数下，对于随雷诺数某个幂次呈指数增长的时间，解始终接近三个 Lamb–Oseen 涡；它们的中心位于无黏三角形各顶点附近。在塌缩情形中，同样的描述一直成立到接近无黏碰撞的时刻。作为扩张情形下长时间结果的应用，我们在黏性消失极限下实现了涡旋爆发：这是点涡的一种非唯一性机制，其中单个涡分裂为一个扩张的 Gröbli 三角形。

Comment: 68 pages, 4 figures, SageMath script in the supplementary material

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

1. 研究对象和背景

论文研究二维不可压 Navier–Stokes 方程中由三个点涡初值产生的涡团演化, 特别是其无黏极限对应的格罗布利三角形自相似扩张或塌缩。涡量形式的方程为

\partial_t\omega+\{\psi,\omega\}=\nu\Delta\omega,\qquad \Delta\psi=\omega,\qquad x\in\mathbb R^2,

其中 $u=\nabla^\perp\psi$ , $\nu>0$ 为运动黏度, $\{f,g\}=\nabla^\perp f\cdot\nabla g$ 。初始涡量是有限个原子测度, 本文取三个点涡:

\omega_{\rm in}=\sum_{i=1}^3\Gamma_i\delta_{z_i^{\rm in}},\qquad \Gamma_i\ne0.

对于这种初值, 黏性方程有全局解; 当黏性很小时, 其涡团中心通常由 Helmholtz–Kirchhoff 点涡系统近似:

z_i'(t)=\sum_{j\ne i}\frac{\Gamma_j}{2\pi}\frac{(z_i(t)-z_j(t))^\perp}{|z_i(t)-z_j(t)|^2},\qquad z_i(0)=z_i^{\rm in}.

点涡模型是二维流体中描述集中涡量相互作用的基本有限维模型, 既保留了流体的若干关键结构, 也产生碰撞、相对平衡与自相似运动等现象。格罗布利在 1877 年发现, 三个点涡在特定环量和初始几何条件下可形成始终相似的三角形, 同时旋转并按平方根时间尺度扩张, 或在有限时间内塌缩。本文将三个点涡的黏性演化与这一有限维轨道相比较, 关心的不是单纯验证某个固定时间区间上的无黏近似, 而是确定该近似在扩张或接近碰撞时能维持多久, 并据此研究黏性趋零时点涡非唯一延拓机制的实现。

将总涡量中心置于原点, 并用 $$d$$ 表示初始三角形的最大边长。格罗布利轨道可写为

z_i(t)=d_{\rm HK}(t)e^{i\theta(t)}n_i,\qquad \max_{i\ne j}|n_i-n_j|=1,

其中 $$n_i$$ 固定, 而

d_{\rm HK}(t)=d\sqrt{1+\operatorname{sign}(\alpha)\frac{t}{T_{\rm adv}}},\qquad \theta(t)=\frac{\beta}{\alpha}\log\left(1+\operatorname{sign}(\alpha)\frac{t}{T_{\rm adv}}\right),\qquad T_{\rm adv}=\frac{d^2}{|\alpha|\Gamma}.

这里 $\Gamma=\max_i|\Gamma_i|$ , $\alpha,\beta$ 由归一化环量和三角形构型决定, $\alpha>0$ 对应扩张, $\alpha<0$ 对应塌缩。格罗布利自相似解要求满足点涡系统的相应守恒量条件, 包括

\sum_{1\le i<j\le3}\Gamma_i\Gamma_j=0,\qquad \sum_{1\le i<j\le3}\Gamma_i\Gamma_j|z_i^{\rm in}-z_j^{\rm in}|^2=0,

并需取相容的初始几何构型。黏性单涡的基本解是 Lamb–Oseen 涡:

\omega(x,t)=\frac{\Gamma_i}{\nu t}G\left(\frac{x-z_i}{\sqrt{\nu t}}\right),\qquad G(\xi)=\frac1{4\pi}e^{-|\xi|^2/4}.

因此自然的近似是三个此类涡核的叠加, 涡核中心跟随点涡轨道。问题的难点在于扩散尺度 $\sqrt{\nu t}$ 随时间增长, 而点涡间距也随格罗布利轨道变化。令 $\delta=\nu/\Gamma$ , 本文使用随轨道尺度变化的有效纵横比

\epsilon(t)=\frac{\sqrt{\nu t}}{d_{\rm HK}(t)},\qquad \epsilon(t)^2=\frac{\delta}{|\alpha|}\frac{t/T_{\rm adv}}{1+\operatorname{sign}(\alpha)t/T_{\rm adv}}.

扩张时 $\epsilon(t)$ 始终为 $O(\sqrt\delta)$ , 即使 $$t$$ 很大, 涡核相对于三角形尺度仍小; 塌缩时 $\epsilon(t)$ 则在 $t\to T_{\rm adv}$ 时增大。既有黏性点涡逼近理论通常在任意固定的无碰撞时间区间成立, 或在特殊平移、旋转构型中推进到接近扩散时间的尺度。扩张格罗布利轨道不属于刚性相对平衡, 但其几何尺度与扩散同步增长, 为突破固定扩散时间障碍提供了结构性可能。本文的另一目标是利用这一长时间逼近, 从一族初始间距趋零的黏性三涡构造出点涡意义下的 vortex burst: 单个原子涡量在极限中分裂成扩张的三角形分支。

2. 主要定理或主要结果

主要长时间结论是定理 1.2 的简化表述; 更强版本为定理 5.1。设三个初始点涡满足格罗布利自相似构型条件, 其无黏轨道为扩张型或塌缩型。令 $\delta=\nu/\Gamma$ , $$d$$ 为初始三角形尺度, $\epsilon(t)=\sqrt{\nu t}/d_{\rm HK}(t)$ 。存在仅依赖归一化环量 $\{\gamma_i\}$ 与归一化构型 $\{n_i\}$ 、而与 $$d$$ 无关的充分小常数 $\delta_0,p_0,p'_0$ 及常数 $$C_0$$ , 使得当 $\delta\le\delta_0$ 且固定 $p\in(0,p_0)$ , 或固定 $p'\in(0,p'_0)$ 时, 可构造调制中心 $$Z_i(t)$$ , 满足

\frac1\Gamma\left\|\omega(t)-\frac1{\nu t}\sum_{i=1}^3\Gamma_iG\left(\frac{\cdot-Z_i(t)}{\sqrt{\nu t}}\right)\right\|_{L^1(\mathbb R^2)} \le C_0\bigl(\epsilon(t)^2+\delta\sqrt{\epsilon(t)}\bigr),

且

\max_i\frac{|Z_i(t)-z_i(t)|}{d_{\rm HK}(t)}\le C_0\sqrt\delta\,\epsilon(t)^2.

扩张情形上述估计对

0<t<T_{\rm adv}\exp(\delta^{-p})

成立; 塌缩情形对

0<t<(1-\delta^{p'})T_{\rm adv}

成立。常数对初始尺度 $$d$$ 一致, 是后续爆发极限的重要条件。该结论表明扩张轨道上的黏性解在拉伸指数级长时间内仍由三个 Lamb–Oseen 涡近似, 而塌缩情形可逼近无黏碰撞时刻, 但仍与碰撞时刻保持相对间距 $\delta^{p'}T_{\rm adv}$ 。作者强调塌缩结论不是本文最强的已知版本: 文中提及 Nguyen 的独立结果可对任意固定 $p'\in(0,1)$ 得到相应逼近。

第二个主要结果是定理 1.8 的黏性趋零爆发极限。设 $d_\nu\to0$ , 并满足

\frac{d_\nu^2}{\Gamma\exp((\Gamma/\nu)^p)}\longrightarrow+\infty \qquad(\nu\to0),

其中 $$p>0$$ 取在上述扩张定理允许的范围内。取适当的初始旋转相位 $\beta_\nu$ , 使其抵消 $T_{\rm adv,\nu}=d_\nu^2/(\alpha\Gamma)$ 趋零时轨道相位中的对数奇异项, 并令

\omega_\nu^{\rm in}=\sum_{i=1}^3\Gamma_i\delta_{d_\nu e^{i\beta_\nu}n_i}.

则 $\omega_\nu^{\rm in}\rightharpoonup^*\Gamma_{\rm tot}\delta_0$ , 其中 $\Gamma_{\rm tot}=\sum_i\Gamma_i$ , 而对每个固定 $$T>0$$ ,

\omega_\nu\rightharpoonup^*\omega_0^{\rm b}\quad\text{于 }C([0,T];\mathcal M(\mathbb R^2)),

其中极限分支为

\omega_0^{\rm b}(t)=\sum_{i=1}^3\Gamma_i\delta_{z_{0,i}(t)},\qquad z_{0,i}(t)=\sqrt{\alpha\Gamma t}\exp\left(i\beta^*+i\frac{\beta}{\alpha}\log\frac{t}{t_*}\right)n_i\quad(t>0),

并约定 $z_{0,i}(0)=0$ 。通过选择允许的扩张构型及极限相位, 可从同一个初始极限原子 $\Gamma_{\rm tot}\delta_0$ 选出不同的扩张三涡分支。该结论是严格的黏性趋零收敛结论; 文中关于之后黏性解最终弛豫为单个 Lamb–Oseen 涡的讨论属于已有长时间行为与可能机制, 不是定理 1.2 或定理 1.8 所证明的内容。

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

证明以每个涡核为中心的自相似变量为基础, 并同时调制涡核形状与涡心位置。写成

\omega(x,t)=\sum_{i=1}^3\frac{\Gamma_i}{\nu t} \bigl(\Omega_{{\rm app},i}+\delta w_i\bigr) \left(\frac{x-Z_i(t)}{\sqrt{\nu t}},t\right),\qquad Z_i(t)=z_i(t)+\rho_i(t).

其中 $$z_i(t)$$ 是精确的点涡轨道, $\rho_i$ 是调制位移, $$w_i$$ 是相对于近似涡核的余项。通过施加零一阶矩条件消除平移模态, 使涡心调制可以从形状误差中分离出来。近似剖面按有效纵横比展开为

\Omega_{{\rm app},i}=G+\sum_{n=2}^4\epsilon(t)^n\Omega_{n,i}.

其关键不是只取高斯主项, 而是把二阶椭圆形变纳入参考剖面。四阶近似的代入残差满足文中给出的估计

\delta^{-1}\mathrm{NS}(\Omega_{\rm app})=O\bigl(\delta^{-2}\epsilon(t)^5+\epsilon(t)\bigr).

扩张时 $\epsilon^2\lesssim\delta$ , 所以该误差在相关长时间范围内仍为 $O(\sqrt\delta)$ ; 塌缩时误差随 $\epsilon$ 放大, 需限制 $$p'$$ 足够小。与固定尺度问题相比, 扩张几何使有效纵横比不随时间无界增长, 因而固定四阶近似即可控制所需时间尺度, 无须随接近扩散时间而不断提高展开阶数。

更深的线性障碍在于 Lamb–Oseen 涡单独提供的稳定性不足。线性化算子具有示意展开

\delta^{-1}L_{\Omega_{\rm app}}=\delta^{-1}L_G+\delta^{-1}\epsilon(t)^2L_2+\cdots.

在目标长时间尺度上, $\delta^{-1}\epsilon^2$ 未必小, 因此不能把二阶形变项当作微扰直接吸收。作者把 $G+\epsilon^2\Omega_2$ 作为主要参考状态, 围绕该状态构造受 Arnold 型泛函框架启发的加权能量估计: 通过近似建立修正涡量与总流函数之间的函数关系, 使线性化无黏部分近似反自伴, 同时保留黏性耗散的强制性。这样, 强的二阶效应进入主算子, 剩余线性误差从 $\delta^{-1}\epsilon^2$ 降至例如 $\delta^{-1}\epsilon^3$ , 在扩张情形可视为小量。相对平衡构型可利用对称性建立所需关系; 格罗布利轨道缺乏这些刚性对称性, 本文的工作之一是在一般动态构型中显式识别到二阶所需的结构与抵消。余项方程中的非线性项则结合该加权估计通过 bootstrap 控制。

另一个核心障碍是调制中心并非被动误差。移动坐标系产生 $T_{\rho_i}(\Omega_{\rm app})$ 等项, 其大小不能仅凭三角形边长接近来控制。证明还需同时掌握

\frac{|\rho_i(t)|}{d_{\rm HK}(t)},\qquad \frac{|t\partial_t\rho_i(t)|}{d_{\rm HK}(t)},

第二项用于排除几何形状看似稳定、但中心快速振荡的情况。作者因此分析格罗布利点涡轨道的线性化, 在适配扩张或塌缩的自相似变量中将时变问题化为一个有限维常系数谱问题, 归结为研究 $3\times3$ 矩阵。所得线性传播子满足

\|U_L(t,s)\|\lesssim \left(1+\left|\log\frac{1+\varsigma t/T_{\rm adv}}{1+\varsigma s/T_{\rm adv}}\right|\right) \left(\frac{d_{\rm HK}(t)}{d_{\rm HK}(s)}\right)^{\varsigma},\qquad \varsigma=\operatorname{sign}(\alpha).

把这一传播子界与涡核余项的 bootstrap 估计结合, 可推出中心及其时间导数的相对误差界。扩张情况下出现 $\bigl(\log(e+t/T_{\rm adv})\bigr)^2$ 型损失, 这是将误差沿时间积分并迭代控制所付出的代价, 也解释了为何证明只能达到 $T_{\rm adv}\exp(\delta^{-p})$ 的拉伸指数尺度, 而不能据此直接推出更长时间有效性。塌缩时则出现 $(1-t/T_{\rm adv})^{-1}$ 型损失, 反映碰撞前线性化放大。

最后, 爆发极限的证明不是另建一套涡核分析。点涡轨道有显式公式, 经相位校正后在固定 $$t>0$$ 收敛到目标自相似分支; 长时间 Lamb–Oseen 近似和中心调制界给出黏性涡量与离散点涡分布的控制, 再结合 Kantorovich–Rubinstein 型距离估计及集中紧性框架, 得到测度空间中对 $t\in[0,T]$ 一致的弱星收敛。这里初始间距条件保证固定物理时间相对于 $T_{\rm adv,\nu}$ 已足够大, 同时仍落在可用的长时间逼近范围内。

4. 与之前工作的比较

在一般黏性点涡逼近方面, Marchioro 的早期工作与 Gallay 的结果奠定了从 Dirac 涡量出发的描述框架; Gallay 证明大 Reynolds 数下解可由沿点涡轨道运动的 Lamb–Oseen 涡叠加近似, 并覆盖无黏轨道存在且无碰撞的任意固定时间区间。Ceci 与 Seis 的 Wasserstein 集中估计适用于更一般的集中初值, 时间可随 Reynolds 数对数增长。本文不改变初值为原子、且构型严格满足格罗布利条件这一较强前提, 而是利用该特定轨道的尺度结构, 在扩张情形把有效时间推进到 $T_{\rm adv}\exp(\delta^{-p})$ 。因此其推进重点是特定构型上的时间尺度, 不是一般初值理论的统一延拓。

与长时间 Lamb–Oseen 逼近的特殊构型结果相比, 文中提到的平移涡偶、总环量非零的一般涡对及刚性旋转涡晶体工作分别由 [13]、[43]、[12] 研究, 相关逼近可推进到任意接近扩散时间的尺度。那些构型的运动具有平移或相对平衡等刚性, 稳定性分析可利用相应对称结构; 格罗布利三角形持续改变尺度并旋转, 不属于相对平衡。本文的关键增益来自 $d_{\rm HK}(t)\sim\sqrt t$ 与黏性核宽 $\sqrt{\nu t}$ 同步增长, 使 $\epsilon(t)$ 在扩张时保持 $O(\sqrt\delta)$ 。代价是结论依赖严格的自相似三涡初值、环量与形状条件, 且扩张时间指数 $$p$$ 仅被证明存在于某个小区间, 文中未给出可用的显式下界。

与一般相对平衡长时间分析的比较还揭示了方法差异。既有研究常以相对于固定最小涡间距定义纵横比, 该量按 $\sqrt{\nu t}/d$ 增长, 因而 $t\sim T_{\rm diff}=d^2/\nu$ 是重要障碍; 要接近该时间尺度, 相关工作需构造更高阶近似。本文改以变化中的格罗布利边长定义有效纵横比, 在扩张时避免这一增长, 所以四阶近似可覆盖更长区间。与此同时, 动态线性化造成的调制中心控制与累积对数损失成为新瓶颈, 使作者得到的是拉伸指数时间而非任意时间的近似。作者还指出扩张涡团最终会弛豫为总环量对应的单个 Lamb–Oseen 涡; Gallay–Rodrigues 的定量弛豫上界为 Reynolds 数的双指数尺度。本文识别的是弛豫发生前的长暂态, 不证明其持续时间最优, 也不证明扩张结构最终如何改变形状或何时合并。

塌缩结果的范围应与扩张结果及后续工作区分。本文证明到 $(1-\delta^{p'})T_{\rm adv}$ , 其中 $$p'$$ 只保证为充分小的正数。作者明确提到 Nguyen [36] 独立得到任意固定 $p'\in(0,1)$ 的塌缩逼近, 因而本文塌缩定理不是当时最强范围。Nguyen 的论证以 $\epsilon(t)^2$ 归一化余项, 并引入适合碰撞的对数时间 $s=-\log(1-t/T_{\rm adv})$ ; 本文则以 $\delta$ 归一化, 主要为扩张问题设计, 塌缩结论由相同框架直接调整而来。 $$p'=1$$ 附近涡间距达到扩散长度量级, 作者将其视为扰动 Lamb–Oseen 涡方法的自然障碍, 并未声称解决碰撞后的黏性动力学。文中还讨论 Gallay–Šverák 关于某些无黏塌缩构型受初值扰动后可能转为近碰撞再反弹扩张的结果; 本文从精确点涡初值出发, 并未处理这类不良备初值下的黏性稳定性。

在无黏爆发机制方面, Grotto 与 Pappalettera [25] 已在广义点涡框架下利用涡爆发构造 Euler 弱解的非唯一性: 单个点涡可在碰撞时刻后接续为扩张格罗布利三角形。本文的贡献不是重新证明该点涡模型中的非唯一性, 而是给出 Navier–Stokes 黏性解到这类分支的 vanishing viscosity 实现。与 Ceci–Seis 的对数时间集中估计所对应的初始间距限制相比, 本文要求

d_\nu\to0,\qquad \frac{d_\nu^2}{\Gamma\exp((\Gamma/\nu)^p)}\to+\infty,

允许 $d_\nu$ 按 Reynolds 数倒数的任意多项式幂趋零, 从而可以在固定正时间把初始收缩的三个涡推进到与黏度无关的分离尺度。该结果还与 [28] 关于更大空间尺度自相似观察下总环量集中于原点的结论形成对照: 本文在平方根尺度 $t^{1/2}$ 上得到非平凡的三个原子分支, 即该临界尺度可保留自相似结构。与 [2] 中通过受迫 Euler 系统和临界扰动选择非径向极限的机制不同, 本文利用的是 Dirac 初值可同时对应静止单涡与扩张格罗布利分支这一点涡非唯一性结构。局限在于结果针对精心选择的相位与间距族, 并不能推出任意逼近同一原子测度的 Navier–Stokes 初值都会选择该分支; 关于其他数据的选择规律, 以及扩张三角形更长时间后的黏性演化, 从提供内容中无法确认。

5. 涡丝的 Arnold 型稳定性分析

An Arnold-Type Stability Analysis of Vortex Filaments

作者: Masashi Aiki, Mitsuo Higaki

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29353v1

摘要翻译: 我们证明，在局部诱导方程（LIE）下，开角为 $θ\le π$ 的圆弧以及闭合圆具有非线性轨道稳定性；对于圆弧，还施加了给定的端点约束。我们的 Arnold 型变分方法将两个守恒量——弯曲能和流体冲量的一个分量——组合成一个泛函，并证明对于所有容许曲线，该泛函与其到相应轨道的 $$H^2$$ 距离的平方可比。对于角度为 $$θ>π$$ 的圆弧，我们利用逆广义 Hasimoto 变换证明其非线性轨道不稳定性。

Comment: 31 pages

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

1. 研究对象和背景

本文研究三维不可压缩理想流体中细涡丝的局部诱导运动模型, 即局部诱导方程（Localized Induction Equation, LIE）。曲线 $x(s,t)\in\mathbb{R}^3$ 以弧长 $$s$$ 参数化, 其演化由切向量与曲率向量的叉积给出：

x_t=x_s\times x_{ss},\qquad |x_s|=1.

这里 $$x_s$$ 是单位切向量, $x_{ss}$ 描述曲率向量。LIE 是细涡丝自诱导运动的经典模型；论文将稳定性对象限定为圆弧和闭合圆, 它们都是形状保持、以速度 $R^{-1}e_3$ 沿 $$e_3$$ 方向平移的精确解。给定半径 $$R>0$$ , 基准曲线为

x_R(s,t)=x_R^0(s)+\frac{t}{R}e_3,\qquad x_R^0(s)=R\cos\!\left(\frac{s}{R}\right)e_1+R\sin\!\left(\frac{s}{R}\right)e_2.

对圆弧, 固定开角 $0<\theta<2\pi$ , 区间长度为 $L=\theta R$ , 空间区间为 $$I=(0,L)$$ , 并施加切向端点条件

x_s(0,t)=e_2,\qquad x_s(L,t)=b_\theta:=-\sin\theta\,e_1+\cos\theta\,e_2.

初始端点还须满足 $x_0(0)\cdot e_2=0$ 与 $x_0(L)\cdot b_\theta=0$ , 使端点分别处于基准圆弧对应的固定平面内。对闭圆, 取周期域 $\mathbb T=\mathbb R/(2\pi R\mathbb Z)$ , 不施加端点条件。问题是：当初始曲线满足弧长参数约束及相应边界或周期条件时, 它与基准圆弧或圆的轨道在 $$H^2$$ 意义下接近, 是否能在任意长时间内保持接近？此外, 对哪些圆弧开角会发生非线性轨道不稳定？

这类问题的重要性在于, LIE 是涡丝动力学的几何模型, 而稳定性刻画了理想化细涡丝形状对扰动的长期响应。论文也把问题放在更广泛的涡动力学背景中：欧拉方程、Navier–Stokes 方程和 Gross–Pitaevskii 方程均有涡丝或涡环相关研究, 但这些模型与 LIE 的方程结构及可用守恒量不同, 不能直接替代本模型中的稳定性分析。LIE 的曲率与挠率演化又可通过 Hasimoto 变换联系到聚焦三次非线性薛定谔方程（NLS）；有限长度涡丝的端点条件对应于 NLS 的齐次 Neumann 边界条件。这种联系为理解开角阈值附近的调制不稳定性提供背景, 但本文的稳定性结论本身仍是对 LIE 曲线及其轨道距离作出的陈述。

作者采用的变分对象由弯曲能与流体冲量的 $$e_3$$ 分量组成：

B[x]=\frac12\int_D|x_{ss}|^2\,ds,\qquad P_3[x]=\frac12\int_D e_3\cdot(x\times x_s)\,ds,

并取

\mathcal L[x]=B[x]-\frac1{R^3}P_3[x].

在所考虑的圆弧边界条件下, $$B$$ 与 $$P_3$$ 均沿 LIE 守恒；其中 $$P_3$$ 的守恒依赖端点约束。闭圆情形则在周期域上使用同一组合。系数 $R^{-3}$ 的选取使基准圆弧或圆成为相应 admissible 类中的临界构形, 且 $\mathcal L[x_R^0]=0$ 。研究的关键不只是证明某个局部二次变分正定, 而是要在满足单位切速和边界或周期约束的整个 admissible 类上, 将该守恒泛函与到对称轨道的平方 $$H^2$$ 距离定量比较。

已有结果提供了直接背景。文献 [3] 对圆弧给出过非线性稳定估计： $0<\theta<\pi$ 时可在 $$H^2$$ 中模去沿 $$e_3$$ 平移控制扰动；但当 $\pi\leq\theta<2\pi$ 时, 其类似结论要求额外的中点反射对称性。文献 [4] 已证明闭圆在非对称扰动下的轨道稳定性, 但其陈述使用初始流体冲量相对基准值及高 Fourier 模态的定量条件, 证明中分别以相对能量控制 $|n|\geq2$ 模态、以向量流体冲量控制 $n=\pm1$ 模态。本文试图推进的空白, 是以同一个守恒变分泛函统一处理稳定圆弧与闭圆, 证明其对轨道距离具有全局定量控制, 并在圆弧问题中识别开角 $\theta=\pi$ 的稳定性阈值及阈值以上的不稳定性。

2. 主要定理或主要结果

主要结论可分为圆弧稳定、圆弧不稳定和闭圆稳定三部分。

圆弧稳定（定理 A 第 1 项；详细结果见定理 2.5、2.9）。令 $0<\theta\leq\pi$ , $I=(0,\theta R)$ 。初值 $x_0\in H^4(I)^3$ 满足 $|x_{0,s}|=1$ , 端点位置条件

x_0(0)\cdot e_2=0,\qquad x_0(L)\cdot b_\theta=0,

以及相容条件

x_{0,s}(0)=e_2,\quad x_{0,s}(L)=b_\theta,\quad (x_{0,s}\times x_{0,sss})(0)=(x_{0,s}\times x_{0,sss})(L)=0.

相应 LIE 初边值问题存在唯一全局解。定义稳定轨道

\mathcal A_\theta=\begin{cases} x_R^0+Re_3,&0<\theta<\pi,\\ x_R^0+Re_1+Re_3,&\theta=\pi, \end{cases} \qquad d(x,\mathcal A_\theta)=\inf_{y\in\mathcal A_\theta}\|x-y\|_{H^2(I)}.

则存在仅依赖于 $R,\theta$ 的常数 $$c,C>0$$ , 使所有 $t\geq0$ 满足

c\,d(x_0,\mathcal A_\theta)\leq d(x(t),\mathcal A_\theta)\leq C\,d(x_0,\mathcal A_\theta).

因此稳定性是对整个 admissible 初值类的轨道估计, 不要求初始扰动小；半圆 $\theta=\pi$ 也稳定, 且不需要文献 [3] 中的中点反射对称性。

圆弧不稳定（定理 A 第 2 项；详细结果见定理 2.15）。令 $\pi<\theta<2\pi$ 。存在仅依赖于 $R,\theta$ 的常数 $C_0,\eta>0$ , 使得对每个 $0<\varepsilon<1/2$ , 都可构造光滑初值 $x_0^{(\varepsilon)}$ , 满足单位弧长条件、上述端点位置与相容条件, 且

\|x_0^{(\varepsilon)}-x_R^0\|_{H^4(I)}\leq C_0\varepsilon.

相应全局解在某个时间 $T_\varepsilon$ 达到与稳定轨道保持正距离的状态：

d(x^{(\varepsilon)}(T_\varepsilon),\mathcal A_\theta)\geq\eta, \qquad T_\varepsilon\text{ 的量级为 }\log(1/\varepsilon).

这给出非线性轨道不稳定, 并表明阈值 $\theta=\pi$ 将文中证明的稳定与不稳定角度区间分开。摘要及概述将此机制联系到区间上聚焦三次 NLS 的调制不稳定性, 具体构造使用 Akhmediev breather 与逆广义 Hasimoto 变换。

闭圆稳定（定理 B；详细结果见定理 3.5）。令 $L=2\pi R$ , $\mathbb T=\mathbb R/(L\mathbb Z)$ , 初值 $x_0\in H^4(\mathbb T)^3$ 满足 $|x_{0,s}|=1$ , 并令 $Q_z^\alpha$ 表示绕 $$e_3$$ 轴转角 $\alpha$ 的旋转。定义

\mathcal C=\{Q_z^\alpha x_R^0+\tau:\alpha\in\mathbb R,\ \tau\in\mathbb R^3\},\qquad d(x,\mathcal C)=\inf_{y\in\mathcal C}\|x-y\|_{H^2(\mathbb T)}.

对应唯一全局解满足

c\,d(x_0,\mathcal C)\leq d(x(t),\mathcal C)\leq C\,d(x_0,\mathcal C),\qquad t\geq0,

其中 $$c,C>0$$ 仅依赖于 $$R$$ 。轨道同时模去空间平移及绕对称轴旋转；对于基准圆, 后者等价于弧长参数的相位平移。与文献 [4] 的表述相比, 本定理不要求另行假设初始流体冲量满足特定定量条件, 而直接以初始轨道距离控制全时间轨道距离。

变分定量比较。对 $0<\theta\leq\pi$ 的 admissible 圆弧以及 admissible 闭圆, 文中证明存在依赖固定几何的正数 $$c,C$$ , 使

c\,d(x,\mathcal S)^2\leq\mathcal L[x]\leq C\,d(x,\mathcal S)^2,

其中 $\mathcal S$ 分别为 $\mathcal A_\theta$ 或 $\mathcal C$ , 且 $\mathcal L$ 沿 LIE 守恒。该比较不要求 $$x$$ 靠近基准曲线。结合守恒性即得上下轨道距离估计。论文还指出, 对任一非零初始扰动, 解不会渐近趋于相应基准圆弧或圆。另有一个关于保持弧长的受迫方程的先验轨道距离估计作为闭圆证明的推论；提供的摘录未给出该估计的具体公式与假设, 因而无法进一步列明。

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

证明的主线是构造同时利用守恒量与几何约束的 Arnold 型 Lyapunov 泛函, 再将其与商去对称性后的 $$H^2$$ 距离建立定量等价。这里的困难不只是计算基准圆的二阶变分：需要处理单位弧长约束, 圆弧的端点位置与切向条件, 以及闭圆的平移、旋转和参数相位对称性；还要使所得控制对整个 admissible 类成立, 而不是只在基准解的小邻域中成立。

第一项关键结构是弯曲能与流体冲量的组合。作者取

\mathcal L[x]=\frac12\int_D|x_{ss}|^2\,ds-\frac1{2R^3}\int_D e_3\cdot(x\times x_s)\,ds.

弯曲能 $$B$$ 衡量曲率的 $$L^2$$ 大小, $$P_3$$ 是流体冲量沿圆的平移轴方向的分量。系数 $R^{-3}$ 被选为使参考构形满足变分驻定。圆弧情形, 对端点相容的变分 $\phi$ 有 $\phi_s(0)=\phi_s(L)=0$ 及相应端点法向约束；通过分部积分和圆的微分恒等式可验证 $D\mathcal L(x_R^0)[\phi]=0$ 。由于 $\mathcal L$ 是二次泛函且 $\mathcal L[x_R^0]=0$ , 对 $x=x_R^0+\phi$ 得到精确关系 $\mathcal L[x]=\mathcal L[\phi]$ , 没有 Taylor 展开的高阶余项。这一精确二次结构是将静态强制性转化为动态稳定性的关键。

第二项是守恒性与边界项的匹配。将 LIE 对 $$s$$ 求导得 $x_{st}=x_s\times x_{sss}$ , 切向边界条件保证 $x_{st}$ 在两端为零, 从而分部积分证明 $$B$$ 守恒。对 $$P_3$$ 的时间导数, 需要控制边界项 $e_3\cdot(x\times x_t)$ ；端点位置约束与固定端点切向条件共同给出相应边界消失条件, 再利用 $$|x_s|=1$$ 所导出的 $(x_t\times x_s)=x_{ss}$ 及端点切向量差的结构得到 $$P_3$$ 守恒。由此 $\mathcal L$ 守恒。特别是, 圆弧问题中的冲量守恒并非脱离边界条件的形式恒等式, 证明明确使用了端点约束。

第三项也是稳定证明的定量核心, 是在 admissible 曲线上证明 $\mathcal L$ 与轨道距离平方上下比较。下界给出定量刚性：若 $\mathcal L[x]$ 小, 则 $$x$$ 必须在 $$H^2$$ 中接近其全部极小构形；上界则防止该泛函相对于轨道距离过度增长。论文摘要和引言明确说明这种比较对所有 admissible 曲线成立, 不需要小扰动假设。对闭圆, 比较还须正确消除平移与绕 $$e_3$$ 轴旋转（等价于参数相位）的退化方向；对半圆, 需选取相应的轨道代表以处理额外的几何退化。给定摘录没有展示各角度区间的具体不等式推导、辅助引理、傅里叶分解或常数的显式表达, 因此不能据此确认作者在逐项估计中使用的精确技术细节。

不稳定证明走另一条路线：稳定侧的 $\mathcal L$ 强制性在 $\theta>\pi$ 时失效, 单靠该泛函无法给出距离控制。作者利用 Hasimoto 型对应, 将 LIE 曲线问题与聚焦三次 NLS 的有限区间问题联系起来, 选取显式 Akhmediev breather, 再通过逆广义 Hasimoto 变换构造满足 LIE 端点条件的涡丝解。构造的初值与基准弧在 $$H^4$$ 中相差 $O(\varepsilon)$ , 但在 $O(\log(1/\varepsilon))$ 时间后与稳定轨道相距至少固定常数 $\eta$ 。因此该方法提供的是显式小扰动在对数时间内放大的非线性不稳定机制, 而不是仅从线性化谱得到的稳定性判断。摘录未提供 breather 参数、逆变换公式及保持端点条件的逐步验证, 故这些细节无法确认。

最后, 闭圆证明还导出保持弧长的受迫演化的先验轨道距离估计。其思路是沿用上述泛函与轨道距离的定量比较, 但受迫情形下守恒律一般不再可直接使用；具体如何估计强迫项及其时间积分, 在提供的摘录中未展开。总体而言, 真正的证明框架不是泛称的能量法, 而是“边界兼容的守恒组合 + 参考构形的精确驻定二次结构 + 全局定量刚性 + 守恒推出全时间轨道控制”；不稳定侧则以 NLS breather 的显式增长和逆几何变换补足变分框架失去强制性的区域。

4. 与之前工作的比较

与圆弧稳定性结果 [3] 相比, 本文最直接的推进是覆盖了半圆情形 $\theta=\pi$ 而不要求中点反射对称性。文献 [3] 已在 $0<\theta<\pi$ 给出模去沿 $$e_3$$ 平移的全时间 $$H^2$$ 稳定估计, 且不要求 admissible 初值小；对 $\pi\leq\theta<2\pi$ 的类似估计则依赖额外的中点反射对称性。本文通过共同的变分泛函在 $0<\theta\leq\pi$ 证明稳定, 并对 $\theta>\pi$ 给出不稳定构造, 因而在其所研究的角度范围内明确呈现 $\theta=\pi$ 阈值。需要区分的是, 对 $\theta>\pi$ 的不稳定性证明不等同于此前在对称子类上的稳定估计；二者讨论的扰动类别和结论不同。

就估计形式而言, 作者不仅给出从初始距离到后续距离的上界, 还证明上下双边控制

c\,d(x_0,\mathcal S)\leq d(x(t),\mathcal S)\leq C\,d(x_0,\mathcal S).

这来自对守恒泛函与平方轨道距离的双边比较, 并使静态结论具有定量刚性含义： $\mathcal L[x]^{1/2}$ 控制 $d(x,\mathcal S)$ 。其适用类不局限于足够小的邻域。论文还据此指出非零初始扰动不会渐近趋于基准构形。提供的摘录没有给出 [3] 的精确常数、比较不等式或其下界性质, 因此无法进一步逐项比较两篇文章常数的依赖或定量强弱。

对闭圆, 文献 [4] 已经证明非对称扰动下的 $$H^2$$ 轨道稳定, 模去空间平移及绕对称轴旋转。该工作通过相对能量控制 Fourier 模态 $|n|\geq2$ , 并以向量流体冲量控制 $n=\pm1$ 模态；其定理表述包含初始流体冲量相对参考值和高模态的定量条件。本文并非首次证明闭圆的非对称轨道稳定, 而是改变了证明组织方式并加强了结论的直接表述：同一标量组合 $B-R^{-3}P_3$ 对闭圆整个 admissible 类具有与轨道距离平方相当的控制, 稳定估计可直接用初始轨道距离表示, 不需单列文献 [4] 中的初始冲量条件。本文仍依赖弧长约束和固定半径参考几何, 其结果并不意味着任意曲线或任意 LIE 解都具有此类强制性。

相较于通常只在平衡态附近分析二次变分的局部 Arnold 型方法, 本文强调的是全局 admissible 类上的比较不等式。由于泛函是二次的, 且一阶变分在适配边界条件的扰动方向上严格为零, 作者得到精确展开而非带有未控制余项的局部 Taylor 估计。这是其能够不附加小扰动条件证明稳定的结构原因。不过, “全局”指满足单位弧长、指定端点或周期条件的曲线类, 并不表示离开这些约束后仍成立。比较常数依赖固定几何参数；摘录未提供其关于 $R,\theta$ 的显式退化率, 因而不能判断当 $\theta$ 趋近阈值或其他几何极限时估计是否一致。

不稳定性方面, 文献 [14] 研究闭合涡丝的线性稳定性, 结合 Hasimoto 映射与 NLS 的平方本征函数；本文讨论的是非线性、全时间轨道稳定性以及开角超过 $\pi$ 时的非线性轨道不稳定, 两者问题层次不同。文献 [6] 的 Akhmediev breather 提供聚焦 NLS 的显式调制不稳定图像, 文献 [7] 证明了 Akhmediev breather 的不稳定性, 而本文利用逆广义 Hasimoto 变换把相关机制落实到满足端点条件的涡丝。文献 [37] 已有无边界涡丝的相关 breather 构造；本文的特色是将构造用于有限长度圆弧并满足这里的端点约束。具体而言, 初始 $$H^4$$ 误差为 $O(\varepsilon)$ , 在 $O(\log(1/\varepsilon))$ 时间内产生固定量级的轨道偏离。摘录未列出 [37] 的精确假设和时间尺度, 故不能进一步断言本文在所有定量指标上优于该工作。

论文的背景也与更广泛的变分稳定性文献相连。Arnold 的能量-Casimir思想通过守恒结构和约束下的能量强制性研究流体稳定性；流体冲量也曾用于平移涡结构的变分研究, 如 Friedman–Turkington 对涡环的构造及 Choi 对 Hill 球形涡的稳定性分析。Gallay–Šverák 关于平面涡量和 Oseen 涡的工作强调以二次型强制性量化接近极值态的距离。本文与这些工作相通之处是把守恒变分结构转化为定量轨道控制, 不同之处在于对象是由曲线弯曲能与涡丝流体冲量组成的 LIE 泛函, 且必须处理弧长约束及有限区间端点条件。摘录并未提供这些文献与本文之间更具体的定理级假设比较, 不宜推断其方法可直接推广到欧拉、黏性方程或其他边界几何。

总的说, 本文的主要增量不是首次建立圆弧或闭圆稳定, 而是将稳定圆弧和闭圆纳入一个以 $\mathcal L$ 为核心的变分框架, 获得对 admissible 类的全局定量刚性与双边轨道距离估计；圆弧侧去除了半圆稳定所需的反射对称性, 闭圆侧免除了单独的初始冲量定量条件, 并补充了 $\theta>\pi$ 时具对数增长时间尺度的非线性不稳定结果。代价和局限是结论依赖特定的 LIE、固定半径几何、单位弧长约束及所指定的边界条件；不稳定构造依赖 Hasimoto 对应和显式 breather, 而受迫问题的具体估计以及强制性证明中的细部技术, 从提供内容中无法确认。

6. 黏性与热传导理想气体径向对称定常态的无黏极限问题

Inviscid Limit Problem for Radially Symmetric Stationary States of Viscous and Heat-Conductive Ideal Gas

作者: Itsuko Hashimoto

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29074v1

摘要翻译: 我们研究球外域 $\mathbb{R}^n$ ( $n \geq 3$ ) 中黏性热传导多方理想气体入流与出流问题的径向对称定常解的零黏性极限，并施加给定的边界条件和远场条件。假设远场速度为零，且所考虑的解位于远场状态的一个充分小邻域内。在零黏性极限下，动量方程中的黏性项消失，而能量方程中的热传导项仍然保留。对于出流问题，我们证明黏性定常解一致收敛于相应的热传导无黏系统的定常解，收敛速率为 $$O(μ)$$ 。对于入流问题，内边界附近会出现边界层。我们证明，黏性定常解可由极限定常解与边界层剖面的叠加一致逼近。这里，边界层仅出现在密度和速度中，而温度不存在边界层。此外，对于任意 $α\in (0,1)$ ，逼近误差为 $$O(μ^α)$$ 阶。

Comment: 无

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

1. 研究对象和背景

本文研究外区域中的可压缩理想气体径向对称定常流，当剪切黏性与体黏性共同趋于零时，黏性定常解如何逼近相应的无黏极限解。空间区域为 $\Omega=\{x\in\mathbb{R}^n:|x|>r_0\}$ , $n\geq 3$ . 状态变量是密度 $\rho>0$ , 径向速度 $$u$$ , 温度 $\theta>0$ , 状态方程为 $p=R\rho\theta$ , 内能为 $e=c_V\theta$ . 模型保留热传导系数 $\kappa>0$ , 令径向黏性系数 $\mu=2\nu+\lambda>0$ , 并假设 $\nu>0$ , $2\nu+n\lambda\geq0$ .

径向对称设定为 $U(t,x)=x|x|^{-1}u(t,r)$ , $$r=|x|$$ . 定常方程可写为

\begin{aligned} &(r^{n-1}\rho u)_r=0,\\ &(r^{n-1}\rho u^2)_r+r^{n-1}p_r =\mu r^{n-1}\left(\frac{(r^{n-1}u)_r}{r^{n-1}}\right)_r,\\ &\{r^{n-1}\rho u(c_V\theta+\tfrac12u^2)\}_r+(r^{n-1}pu)_r\\ &\qquad=\kappa(r^{n-1}\theta_r)_r+\mu r^{n-1}\left(\frac{(r^{n-1}u)_r}{r^{n-1}}\right)_r u+r^{n-1}\Psi, \end{aligned}

其中 $\Psi=2\nu u_r^2+2\nu(n-1)(u/r)^2+\lambda(u_r+(n-1)u/r)^2$ 是非负黏性耗散项. 远场状态规定为 $(\rho,u,\theta)\to(\rho_+,0,\theta_+)$ 当 $r\to\infty$ . 远场速度取零并非额外的独立限制: 由定常质量守恒， $r^{n-1}\rho u$ 为常数，且 $\rho\to\rho_+>0$ , 因而在 $n\geq3$ 时必有 $u\to0$ .

内边界 $$r=r_0$$ 的条件依流向而异. 流入情形 $$u_b>0$$ 时，黏性问题给定 $(\rho,u,\theta)(r_0)=(\rho_b,u_b,\theta_b)$ ; 流出情形 $$u_b<0$$ 时只给定 $(u,\theta)(r_0)=(u_b,\theta_b)$ , 密度由解确定. $$u_b=0$$ 是不可渗透壁面情形，文中给出其显式解，但主要极限分析针对非零流速的流入与流出.

研究的关键是 $\mu\to0$ 而 $\kappa$ 固定. 极限中动量方程的黏性项消失，但能量方程仍含热传导，因此极限不是绝热欧拉系统，而是质量方程、无黏动量方程与导热能量方程耦合的系统：

\begin{aligned} &(r^{n-1}\rho_Eu_E)_r=0,\\ &(r^{n-1}\rho_Eu_E^2)_r+r^{n-1}(p_E)_r=0,\\ &\{r^{n-1}\rho_Eu_E(c_V\theta_E+\tfrac12u_E^2)\}_r+(r^{n-1}p_Eu_E)_r =\kappa(r^{n-1}(\theta_E)_r)_r, \end{aligned}

其中 $p_E=R\rho_E\theta_E$ , 远场状态仍为 $(\rho_+,0,\theta_+)$ , 且温度边界值为 $\theta_E(r_0)=\theta_b$ . 亚声速范围内，极限系统的质量-动量部分在内边界只能施加一个边界条件，而黏性流入问题却指定了密度和速度两个量. 这类边界条件数目的失配是流入边界层产生的机制. 流出问题则可直接以 $$u_E(r_0)=u_b$$ 选定极限解，无须边界层修正.

既有基础包括作者此前建立的小扰动径向定常黏性解存在唯一性及其空间衰减估计 [4]. 本文的目标是在该理论上刻画黏性定常态的零黏性极限，区分流入、流出两种边界机制，并具体确定流入边界层中哪些物理量发生修正、边界层的空间尺度以及一致逼近误差. 该问题与一般非定常无黏极限文献所讨论的情形不同: 这里聚焦外区域、径向对称、定常小扰动解，而且热传导在极限中保留.

2. 主要定理或主要结果

以下结果均针对 $n\geq3$ , $$u_+=0$$ , $\kappa>0$ 固定及充分小的小扰动数据. 定理 1.2 与定理 1.3 中的小量阈值和常数可取为与 $\mu\in(0,1]$ 无关；流入边界层定理中的常数还可依赖于所选的 $\alpha$ .

流出情形, 定理 1.2. 若 $u_b\leq0$ 且 $|u_b|+|\theta_b-\theta_+|\leq\varepsilon_0$ , 则极限导热无黏系统存在唯一光滑定常解 $(\rho_E,u_E,\theta_E)$ , 满足 $$u_E(r_0)=u_b$$ , $\theta_E(r_0)=\theta_b$ , 以及

|\rho_E(r)-\rho_+|+|\theta_E(r)-\theta_+|\leq C r^{-(n-2)}(|u_b|+|\theta_b-\theta_+|),

并且

C^{-1}r^{-(n-1)}|u_b|\leq |u_E(r)|\leq Cr^{-(n-1)}|u_b|,\qquad r\geq r_0.

黏性定常解 $(\rho,u,\theta)$ 与该极限解在整个外区域上一致收敛，且有一阶黏性误差估计

\sup_{r\geq r_0}|(\rho,u,\theta)(r)-(\rho_E,u_E,\theta_E)(r)|\leq C\mu.

因此流出情形不存在需要加入的边界层项.

流入情形, 定理 1.3. 若 $$u_b>0$$ 且 $|u_b|+|\rho_b-\rho_+|+|\theta_b-\theta_+|\leq\varepsilon_0$ , 极限解由温度边界条件和质量通量匹配条件选定:

\theta_E(r_0)=\theta_b,\qquad \rho_E(r_0)u_E(r_0)=\rho_bu_b.

该解唯一、光滑，并满足

|\rho_E(r)-\rho_+|+|\theta_E(r)-\theta_+| \leq Cr^{-(n-2)}(|u_b|+|\theta_b-\theta_+|),

以及

C^{-1}r^{-(n-1)}|u_b|\leq|u_E(r)|\leq Cr^{-(n-1)}|u_b|.

令 $\epsilon=\rho_bu_b$ , $y=(r-r_0)/\mu$ . 由边界层方程

(\rho_Bu_B)_y=0,\qquad (\rho_Bu_B^2)_y+(p_B)_y=(u_B)_{yy},\qquad \kappa(\theta_B)_{yy}=0,

并施加远离边界时匹配 $(\rho_E(r_0),u_E(r_0),\theta_b)$ 、在 $$y=0$$ 匹配流入边界数据 $(\rho_b,u_b,\theta_b)$ , 得到唯一光滑边界层解. 其温度严格为常数 $\theta_B(y)\equiv\theta_b$ , 而密度、速度指数趋近极限解的边界值：

\begin{aligned} |\rho_B(y)-\rho_E(r_0)|+|u_B(y)-u_E(r_0)| &\leq C\exp\left(-\frac{\nu_0}{\epsilon}y\right)\\ &\quad\times(|u_b|+|\rho_b-\rho_+|+|\theta_b-\theta_+|),\qquad y\geq0. \end{aligned}

由 $y=(r-r_0)/\mu$ , 衰减长度为 $O(\mu)$ 量级. 黏性解由极限解与边界层修正叠加逼近. 对任意 $\alpha\in(0,1)$ ,

\sup_{r\geq r_0}\left|\rho(r)-\rho_E(r)-\rho_B\left(\frac{r-r_0}{\mu}\right)+\rho_E(r_0)\right|\leq C\mu^\alpha,

\sup_{r\geq r_0}\left|u(r)-u_E(r)-u_B\left(\frac{r-r_0}{\mu}\right)+u_E(r_0)\right|\leq C\mu^\alpha,

并且

\sup_{r\geq r_0}|\theta(r)-\theta_E(r)|\leq C\mu^\alpha.

叠加式中减去 $\rho_E(r_0)$ 与 $$u_E(r_0)$$ 是为避免把极限解的边界值重复计算. 文中给出的误差是每个 $\alpha<1$ 均成立的 $O(\mu^\alpha)$ , 并未据此宣称 $O(\mu)$ 或 $O(\mu|\log\mu|)$ .

已有定常解背景. 定理 1.1 回顾 [4] 的黏性定常解结果: 对流入数据的小量条件为 $|u_b|+|\rho_b-\rho_+|+|\theta_b-\theta_+|\leq\varepsilon_0$ ; 对流出数据则为 $|u_b|+|\theta_b-\theta_+|\leq\varepsilon_0$ . 两种情形下密度、温度偏离远场的量均为 $O(r^{-(n-2)})$ , 速度为 $O(|u_b|r^{-(n-1)})$ . 本文的极限定理在此基础上进一步给出极限解、边界层和关于 $\mu$ 的一致误差率.

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

证明的结构性困难首先来自方程退化方式的不均匀: 动量方程失去二阶黏性项，但能量方程中的 $\kappa$ -热传导仍保留. 因此不能把三条方程作为同一类黏性扩散方程直接取极限. 同时，定常质量守恒把速度与密度耦合为常数通量关系，而流入、流出边界所给的条件数不同，决定了是否必须引入边界层.

作者先利用质量守恒

r^{n-1}\rho u=\epsilon_0

消去速度，并改用比容 $v=1/\rho$ 及远场扰动 $\eta=v-v_+$ , $\chi=\theta-\theta_+$ . 例如归一化 $$r_0=1$$ 后，速度可写为 $u=\epsilon_0v/r^{n-1}$ . 这一替换将动量方程改写成关于 $\eta$ 的方程，并利用远场条件积分，得到带有从 $$r$$ 到无穷远积分项的表达式. 能量方程也经积分改写为温度的一阶方程，其中积分常数由 $\chi(1)=\chi_b$ 确定. 这种积分形式的作用是把无穷远边界条件和径向衰减直接嵌入方程，避免在无界区间上另行处理不兼容的两端数据. 前一工作 [4] 已用近似序列和压缩映射建立黏性问题的存在唯一性，并给出对 $\mu\in(0,1]$ 一致的衰减估计; 本文的分析以这些统一估计为基础.

流出情形的关键是构造与黏性解相匹配的极限定常解. 在质量守恒下，指定 $$u_E(r_0)=u_b$$ 可确定通量；动量方程的无黏关系再把比容扰动与温度扰动联系起来，能量方程则保留热扩散. 文中把这些关系改写为积分方程，并证明小扰动解存在唯一，同时得到 $r^{-(n-2)}$ 的状态量衰减和 $r^{-(n-1)}$ 的速度衰减. 之后将黏性解与极限解作差，目标是控制整个 $r\geq r_0$ 上的差，并获得 $C\mu$ 的一致界. 提供的 PDF 摘录尚未包含第 3 节完整证明，因此无法从所给内容确认该差值估计中逐项使用的能量恒等式、权函数或闭合不等式细节; 不应据此补写未显示的技术步骤.

流入情形的核心处理是正确选择极限边界数据. 黏性解的通量在边界由 $r_0^{n-1}\rho_bu_b$ 固定，而无黏质量-动量子系统只能接受一个内边界条件. 作者因此对极限解施加通量匹配 $\rho_E(r_0)u_E(r_0)=\rho_bu_b$ , 另以 $\theta_E(r_0)=\theta_b$ 处理热方程. 密度与速度的剩余边界失配由内边界层承担，而非强行令极限解同时满足过多的边界值.

边界层采用拉伸坐标 $y=(r-r_0)/\mu$ . 在此尺度下保留动量扩散的主导项，形式上令 $\mu\to0$ 得到

(\rho_Bu_B)_y=0,\quad(\rho_Bu_B^2)_y+(p_B)_y=(u_B)_{yy},\quad\kappa(\theta_B)_{yy}=0.

匹配条件分别规定 $$y=0$$ 处的黏性边界数据和 $y\to\infty$ 时的极限解边界值. 温度方程的有界匹配解为 $\theta_B\equiv\theta_b$ , 因而温度没有独立的边界层; 密度和速度则通过质量通量与无黏动量关系形成非平凡剖面. 定理 1.3 给出该剖面的指数衰减率 $\exp(-\nu_0y/\epsilon)$ , 也即物理变量中的 $\exp[-\nu_0(r-r_0)/(\epsilon\mu)]$ .

最后，作者比较黏性解与“极限解加边界层减去匹配边界值”的复合近似. 由于边界层在 $O(\mu)$ 区域内变化快，直接用普通的全域光滑误差展开并不能自动给出一致逼近; 必须同时利用外区极限解的衰减、边界层的指数局域性和黏性解关于 $\mu$ 一致的小扰动估计. 定理 1.3 证明的结论是任意 $\alpha\in(0,1)$ 下的 $O(\mu^\alpha)$ 一致误差. 不过所给摘录止于第 3 节开端，未包含第 4 节及完整证明，故无法确认从复合近似残差到 $O(\mu^\alpha)$ 的具体引理、插值或分区估计，也无法判断该次优幂率是否由特定误差项所限. 文中明确可见的实质分析框架是质量通量降维、积分型扰动方程、拉伸变量边界层和复合近似比较; 更细的证明机制从提供内容中无法确认.

4. 与之前工作的比较

本文在引言中把相关既有工作分成一般可压缩 Navier-Stokes 无黏极限研究与作者此前针对外区域径向定常流的工作. 作者列举了文献 [1,6,5,8,9,10] 作为一般无黏极限研究的代表，但提供的摘录没有给出这些文献的具体定理、假设及结论，因而无法可靠地逐篇比较其正则性门槛、维数、时间区间或收敛拓扑. 本文明确强调的差别是研究对象为流入、流出边界条件下的外区域径向对称定常解；作者称此前除其文献 [3] 外没有相应的此类结果.

最直接的前序工作是作者的 [3], 该文处理相应的径向定常外区域问题，但气体为 barotropic 黏性气体. 当前论文将分析扩展到具有独立温度变量的多方理想气体，并加入热传导及黏性耗散. 这一扩展不仅增加一个未知量: 当 $\mu\to0$ 时，动量方程的黏性项消失，能量方程中的热传导项却不消失，故极限方程仍是导热系统. 相比只需处理质量与动量关系的 barotropic 情形，本文必须同时构造温度极限剖面并协调其边界值与质量-动量部分的条件数.

边界条件方面，论文的推进在于分别识别两类流向的极限结构. 流出时内边界只规定 $u_b,\theta_b$ , 极限问题可直接取 $$u_E(r_0)=u_b$$ , $\theta_E(r_0)=\theta_b$ , 并证明无边界层的一致 $O(\mu)$ 收敛. 流入时黏性问题还规定 $\rho_b$ , 而极限质量-动量系统只能接受一个边界条件. 作者以通量匹配 $\rho_E(r_0)u_E(r_0)=\rho_bu_b$ 选取外层极限解，余下的密度、速度边界值差由厚度 $O(\mu)$ 的边界层吸收. 与此同时，温度边界数据已由导热极限方程处理，边界层温度恒为 $\theta_b$ . 这种“只有密度和速度出现边界层”的结论是本文相对于单纯形式取极限更具体的结构信息.

收敛率与适用范围也有明确区分. 流出问题获得整个外区域上的 $O(\mu)$ 一致误差，而流入复合近似仅证明对每个 $\alpha\in(0,1)$ 有 $O(\mu^\alpha)$ 误差. 因而不能把两种流向概括成相同的一阶收敛; 流入结果对边界层作了显式修正，但其已证明的余项幂率较弱. 论文的结论建立在 $n\geq3$ , 远场速度为零、边界与远场状态充分接近、径向对称和定常假设之上，不能据此推断一般多维非对称流、非小扰动解、非定常演化或声速/超声速边界流也满足相同结论.

关于定常解存在性，定理 1.1 引用 [4] 的结果: 小扰动黏性定常解唯一存在，并有密度、温度 $r^{-(n-2)}$ 衰减和速度 $r^{-(n-1)}$ 衰减; 对 $\mu\in(0,1]$ 相关小量阈值及常数可统一选择. 本文不是重新建立完整的黏性定常解理论，而是在这一统一基础上构造无黏导热极限解、流入边界层并证明极限误差. 作者还提到 [7] 已证明流出定常解在小扰动下的渐近稳定性，但该稳定性问题与本文定常零黏性极限的论证目标不同.

本文明确指出将 $\mu\to0$ 时固定 $\kappa>0$ . 同时令热传导系数趋零的联合极限可能改变能量方程结构，并可能出现不同的热边界层; 这被作者列为有待研究的问题，而非本文已证明的结论. 此外，提供的文本没有列出 [3]、[4]、[7] 和一般文献 [1,5,6,8,9,10] 的完整书目信息或逐项结果，故无法进一步核验前人工作的精确参数范围、速率及假设差异; 对这些细节应标记为从提供内容中无法确认.

7. 轴对称 Navier–Stokes 方程的一个对数型正则性判据

A logarithmic regularity criterion for the axisymmetric Navier-Stokes equations

作者: Stan Palasek

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29012v1

摘要翻译: 对于 $\mathbb R^3$ 中 Navier–Stokes 方程的轴对称解，环量 $$Γ=ru^θ$$ 满足最大值原理，由此得到的关于 $$Γ$$ 的先验界具有尺度不变性。尽管如此，全局正则性问题仍完全悬而未决。以环量表述的最尖锐已知正则性判据略微次临界，即在对称轴附近满足 $|Γ|\lesssim(\log r^{-1})^{-k}$ ；Lei 和 Zhang 证明了 $$k=2$$ 条件下的正则性，随后 Wei 将指数进一步降低至 $$k=3/2$$ 。本文将这一判据进一步加强为 $$k=5/4$$ 。论证沿用对 $$ω^r/r$$ 和 $$ω^θ/r$$ 进行能量估计的常规路线，但现在对 $$ω^r/r$$ 方程中的涡量拉伸项给出了更尖锐的界，其中损失仅与内层涡核和远场之间的尺度数目的平方根成正比。

Comment: 16 pages, 1 figure

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

1. 研究对象和背景

本文研究三维不可压缩、轴对称纳维–斯托克斯方程的条件正则性问题, 重点是旋流在对称轴附近的衰减能否控制解的全局正则性. 速度场写为柱坐标形式 $u=u^r e_r+u^\theta e_\theta+u^z e_z$ , 方程取单位黏性系数:

\partial_t u-\Delta u+u\cdot\nabla u+\nabla p=0,\qquad \nabla\cdot u=0\quad\text{于 }\mathbb R^3.

令 $$r$$ 为到对称轴的距离, 并定义环量 $\Gamma=ru^\theta$ . 对应涡量分量记为 $\omega^r,\omega^\theta,\omega^z$ . 轴对称结构使 $\Gamma$ 满足最大值原理, 即 $\|\Gamma(t)\|_{L^\infty}\leq\|\Gamma_0\|_{L^\infty}$ . 该控制在纳维–斯托克斯缩放 $u_\lambda(x,t)=\lambda u(\lambda x,\lambda^2t)$ 下具有临界性, 因而是轴对称问题中一个重要的先验结构. 但仅有这一临界界并未解决全局正则性问题, 该问题仍然开放.

为表达旋流如何影响涡量, 引入无旋流速度 $$b=u^r e_r+u^z e_z$$ , 以及重标度涡量 $\Omega=\omega^\theta/r$ . 由流函数 $\psi$ 恢复 $$b$$ 的关系为

-\left(\Delta-\frac{1}{r^2}\right)\psi=\omega^\theta=r\Omega,\qquad u^r=-\partial_z\psi,\qquad u^z=\frac{1}{r}\partial_r(r\psi).

环量和 $\Omega$ 满足耦合的漂移扩散系统

\partial_t\Gamma+b\cdot\nabla\Gamma-\left(\Delta-\frac{2}{r}\partial_r\right)\Gamma=0,

\partial_t\Omega+b\cdot\nabla\Omega-\left(\Delta+\frac{2}{r}\partial_r\right)\Omega =\frac{\partial_z(\Gamma^2)}{r^4}.

当 $u^\theta=0$ 时, 方程退化为 $\omega^\theta/r$ 的漂移扩散问题, 全局正则性已有经典结果. 有旋流时, $\omega^\theta$ 方程出现由 $u^\theta$ 产生的源项, 使这一机制不再直接闭合. 因此, 一个自然的研究方向是加强环量在轴附近的控制, 并寻求仅依赖旋流衰减的正则性判据.

已有工作表明, 若 $|\Gamma|\lesssim r^\alpha$ 且 $\alpha>0$ , 可由 Prodi–Serrin–Ladyzhenskaya 型条件推出正则性. 这促使研究者考察比临界有界性略强、但只含对数衰减的条件

|\Gamma(t,r,z)|\lesssim\left(\log\frac{1}{r}\right)^{-k}\qquad (r\to0).

该条件的 $$k=0$$ 对应环量有界这一临界先验估计; $$k>0$$ 则给出对数意义上的次临界改善. Q. S. Zhang 证明任意 $$k>2$$ 时正则, Lei–Zhang 将指数推进至 $$k=2$$ , Wei 随后推进至 $$k=3/2$$ . 本文试图进一步降低这一指数, 核心空白在于改进涡量方程中拉伸项的估计, 以降低能量论证对对数尺度数的损失. 研究对象主要是从 $$H^2$$ 轴对称无散度初值出发的强解, 并以其最大存在时间是否有限为判据; 文中还讨论一个局部适定弱解的推广版本。

2. 主要定理或主要结果

主要结论是定理 1.1. 存在绝对常数 $$c_*>0$$ , 使得若 $u_0\in H^2(\mathbb R^3)$ 为轴对称无散度初值, $\Gamma_0\in L^\infty$ , 且 $u\in C([0,T_*);H^2)$ 是相应强解, $$T_*$$ 为最大存在时间, 则对任意 $\delta_0\in(0,1/2)$ , 若

|\Gamma(t,r,z)|\leq c_*\left(\log\frac{1}{r}\right)^{-5/4} \quad\text{对所有 }0<r\leq\delta_0,\ z\in\mathbb R,\ t\in(0,T_*),

便有 $T_*=\infty$ . 该结论将单纯环量对数衰减判据中的指数从先前的 $$3/2$$ 降至 $$5/4$$ .

定理 1.2 给出用于推出定理 1.1 的定量版本. 令

J=\frac{\omega^r}{r}=-\frac{\partial_z u^\theta}{r},\qquad \Omega=\frac{\omega^\theta}{r},

并以 $dx=r\,dr\,dz$ 定义轴对称积分范数 $\|g\|_2^2=\int |g|^2r\,dr\,dz$ . 记

M=\|\Gamma_0\|_{L^\infty},\qquad M_1=(1+M)\|u_0\|_2,\qquad M_0=(\|J_0\|_2+\|\Omega_0\|_2)M_1^3.

存在绝对常数 $$c_0,C_0>0$$ , 若 $r_0,\varepsilon,S\geq0$ 满足

\varepsilon^4(1+S)^5\leq c_0, \qquad S\geq\log\left(C_0\varepsilon^{-5}\max\{M_0^{1/4},r_0^{-1/2}M_1\}\right),

且在整个存在时间内

|\Gamma(t,r,z)|\leq\varepsilon \quad\text{于 }\{r\leq r_0\}\times(0,T_*),

则 $T_*=\infty$ . 这里 $$S$$ 表示证明可容忍的、从内侧涡核到外场之间的对数尺度数. 对给定小 $\varepsilon$ , 条件 $\varepsilon^4(1+S)^5\leq c_0$ 允许的最大尺度约为 $S_*(\varepsilon)=c_0^{1/5}\varepsilon^{-4/5}-1$ . 若记 $X=\max\{M_0^{1/4},r_0^{-1/2}M_1\}$ , 两个约束在 $$X$$ 大时对应 $\varepsilon\lesssim(\log X)^{-5/4}$ . 令 $r_0\to0$ 并取 $\varepsilon=c_*(\log r_0^{-1})^{-5/4}$ , 即可得到定理 1.1.

作者还给出初值层面的推论: 若 $\|\Gamma_0\|_{L^\infty}\leq\varepsilon\leq c_0^{1/4}$ , 且 $S_*(\varepsilon)\geq\log(C_0\varepsilon^{-5}M_0^{1/4})$ , 则由环量最大值原理可直接应用定量判据, 得到全局正则性. 文中指出, 对大 $$M_0$$ , 这改善了 Wei 的相应初值小性条件.

另有两项扩展性结论. 其一, 定理 1.1 的点态对数衰减可放宽为平均型条件. 令 $\Gamma^\sharp(r)=\sup_{z,t}|\Gamma(t,r,z)|$ , 若

A:=\limsup_{\ell\to\infty}\ell^{3/4}\Gamma^\sharp(e^{-\ell})<\infty,

且

\limsup_{\ell\to\infty}\ell^{3/2}\int_\ell^{2\ell}\Gamma^\sharp(e^{-\ell'})^2\,d\ell' \leq c\min\{1,A^{-2}\},

则作者说明同样可推出正则性. 这允许环量在稀疏的一列环带上达到约 $(\log r^{-1})^{-3/4}$ 的大小, 只要指定对数窗口上的平方平均衰减足够快. 其二, 对 Lei–Ren 的局部定量判据, 轴对称适定弱解在 $$Q(1)$$ 中满足

\|\Gamma\|_{L^\infty(Q(1))} \leq C^{-1}\big((\min\{H,G\}+1)\log(G+2)\big)^{-5/4}

时, 可得其在 $$Q(1/2)$$ 中正则, 其中 $$G,H$$ 是该文所用的局部临界量。此处是对既有局部定理指数的改进, 并非定理 1.1 的全局强解结论。

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

证明沿用 $J=\omega^r/r$ 与 $\Omega=\omega^\theta/r$ 的耦合能量路线, 但关键在于分别处理两个非线性相互作用, 并改善拉伸项的尺度损失. 这两个量满足

\partial_t\Omega+b\cdot\nabla\Omega-\left(\Delta+\frac{2}{r}\partial_r\right)\Omega =-2\frac{u^\theta}{r}J,

\partial_tJ+b\cdot\nabla J-\left(\Delta+\frac{2}{r}\partial_r\right)J =(\omega^r\partial_r+\omega^z\partial_z)\frac{u^r}{r}.

相应能量恒等式提供 $\|\nabla J\|_2^2$ 与 $\|\nabla\Omega\|_2^2$ 耗散, 但还包含

T_1=-4\int\frac{u^\theta}{r}J\Omega\,dx, \qquad T_2=\int J(\omega^r\partial_r+\omega^z\partial_z)\frac{u^r}{r}\,dx.

作者取加权能量 $A(t)=\|J\|_2^2+\mu\|\Omega\|_2^2$ , 通过选择权重 $\mu$ 同时吸收 $$T_1,T_2$$ . 困难在于 $|u^\theta|/r=|\Gamma|/r^2$ 在轴附近具有强奇异权重, 而 $J,\Omega$ 一般不在轴上消失, 因此不能直接使用带 $r^{-2}$ 权重的标准 Hardy 不等式.

控制 $$T_1$$ 时, 作者采用 Wei 的截断 Hardy 型不等式. 在 $r\leq r_1$ 假设 $|\Gamma|\leq\varepsilon$ , 并取 $r_1\leq\varepsilon K/\lambda(t)$ , 其中

\lambda(t)=\big(\|J(t)\|_2\|u^\theta(t)/r\|_2\big)^{1/2} \leq\big(\|J(t)\|_2\|\nabla u(t)\|_2\big)^{1/2}.

截断不等式给出, 对光滑轴对称 $h\in H^1$ ,

\int\frac{|u^\theta|}{r}|h|^2\,dx \leq L\|\partial_rh\|_2^2 +O\left(\frac{L+M}{r_1^2}\int_{r\geq r_1/2}|h|^2\,dx\right), \qquad L=\varepsilon(1+\log K)^2.

其核心不是普通的点态 Hardy 控制, 而是先在 $r\leq r_1$ 使用环量小性, 再在外侧使用最大值原理给出的 $|\Gamma|\leq M$ , 并利用截断函数将外场误差转化为能量可控项. 内涡核半径由 $\lambda(t)$ 决定; 引理 2.3 给出的

\int_0^\rho\sup_z|u^\theta(t,r,z)|\,dr\leq\rho\lambda(t)

保证了在轴附近积分不会因仅有 $|\Gamma|\leq\varepsilon$ 而发散. 这一估计改进了 Wei 使用的速度尺度控制, 因为它将对 $$z$$ 的上确界移入了径向积分所需的估计中。

改进的关键在 $$T_2$$ . 由分部积分及涡量关系, 将其重写为

T_2=\int u^\theta\left(\partial_zJ\,\partial_r\frac{u^r}{r} -\partial_rJ\,\partial_z\frac{u^r}{r}\right)dx.

两个乘积不能用同一种方式估计. 对含 $\partial_r(u^r/r)$ 的一项, 写 $u^\theta=\Gamma/r$ , 并利用 $r^{-1}\partial_r(u^r/r)$ 是 $$u^r/r$$ 的 Hessian 分量. 引理 2.2 给出

\|\nabla^2(u^r/r)\|_2\leq\|\partial_z\Omega\|_2,

故该项可直接由 $\varepsilon\|\partial_zJ\|_2\|\partial_z\Omega\|_2$ 控制, 不需 Hardy 不等式. 对另一项, 作者不把旋流权重放在轴向导数上后再做径向 Hardy 估计, 而是先沿 $$z$$ 方向积分. 引理 3.2 的迹估计为

\sup_{r>0}\int_{\mathbb R}|\partial_zg(r,z)|^2\,dz \lesssim\|\Delta g\|_2^2,

并用于 $$g=u^r/r$$ . 它将每个固定半径圆柱上的轴向 $$L^2$$ 范数控制在体积 $$H^2$$ 信息之下. 于是环量只需出现在径向量 $\int_0^{r_1}\sup_z\Gamma^2\,dr/r$ 中; 涡核处用引理 2.3 细致处理, 涡核外则由 $|\Gamma|\leq\varepsilon$ 得到该量为 $O(\varepsilon^2S)$ , 其中 $S\simeq\log K$ .

因此, 拉伸项满足的系数尺度为 $\Sigma\sim\varepsilon\sqrt{S}$ , 而不是 Wei 方法中的 $\Sigma\sim\varepsilon S$ . 与 $$T_1$$ 的 Hardy 系数 $L\sim\varepsilon S^2$ 合并, 吸收两项要求

\Sigma^2\lesssim\mu\lesssim L^{-2}, \qquad\text{从而}\qquad \Sigma L\lesssim1.

新估计使这一条件成为 $\varepsilon^2S^5\lesssim1$ ? 更准确地代入 $\Sigma\sim\varepsilon\sqrt S$ 与 $L\sim\varepsilon S^2$ , 得 $\Sigma L\sim\varepsilon^2S^{5/2}\lesssim1$ , 即

\varepsilon^4S^5\lesssim1.

这正是定理 1.2 中的尺度约束, 并经由所需尺度 $S\sim\log r_0^{-1}$ 转换为对数指数 $$5/4$$ . 能量论证还需控制截断区域之外的误差; 作者以 $|\Gamma|\leq M$ 及速度能量估计控制远场, 这产生定理中 $$S$$ 对 $$M_0,r_0,M_1$$ 的下界. 最后, 对照文中采用的延拓判据, 只需在最大存在时间之前一致控制 $\|\Omega(t)\|_2$ , 即可排除有限时刻失去正则性。

4. 与之前工作的比较

就环量对数衰减判据而言, 本文的直接比较对象是 Q. S. Zhang、Lei–Zhang 和 Wei 的结果. Q. S. Zhang 在 $|\Gamma|\lesssim(\log r^{-1})^{-k}$ 下证明任意 $$k>2$$ 时正则; Lei–Zhang 达到 $$k=2$$ , 其结论针对 $H^{1/2}$ 初值并与旋流的形式有界性条件相关; Wei 对任意 $u_0\in H^2$ 且 $\Gamma_0\in L^\infty$ 将指数降至 $$k=3/2$$ . 本文保持 $$H^2$$ 强解框架及环量轴附近的点态小性思路, 将指数进一步降至 $$k=5/4$$ . 因而推进主要来自拉伸项估计的改进, 并非取消轴对称性或建立无条件全局正则性。

相较 Wei 的具体技术, $$T_1$$ 仍使用截断 Hardy 不等式, 其损失仍为 $L\sim\varepsilon S^2$ ; 改动发生在 $$T_2$$ . Wei 使用同一类加权 Hardy 控制拉伸项, 得到 $\Sigma\sim\varepsilon S$ , 从而要求 $\varepsilon^2S^3\lesssim1$ , 对应 $$k=3/2$$ . 本文通过将 $$T_2$$ 改写为两个导数乘积, 一项用 $$u^r/r$$ 的 Hessian 控制, 另一项用轴向迹不等式控制, 将损失降为 $\Sigma\sim\varepsilon\sqrt S$ . 与未改进的估计相比, 允许的尺度窗口由 $$S$$ 的 $\varepsilon^{-2/3}$ 型量级扩展为 $\varepsilon^{-4/5}$ 型量级, 因而得到 $$5/4$$ . 这说明改进来自对涡量拉伸结构的非对称处理, 而不是单纯增强已有能量不等式。

定量结论也有区别. 定理 1.2 不只给出轴附近的环量阈值, 还明确要求 $$S$$ 足以覆盖由初值和截断半径决定的远场尺度:

S\geq\log\left(C_0\varepsilon^{-5}\max\{M_0^{1/4},r_0^{-1/2}M_1\}\right), \qquad \varepsilon^4(1+S)^5\leq c_0.

因此, 结论仍带有数据依赖的尺度代价, 并没有证明任意有界环量都足以正则. 作者还指出初值判据比 Wei 的一个对应小性阈值更好, 但阈值仍含初值量 $$M_0$$ ; 仅由临界范数构造类似判据是否可能, 文中明确称为开放问题, 并提及已有相关进展。不能据此将本文解读为已得到只依赖临界数据的全局定理。

本文与其他轴对称正则性准则的关系也有边界. Chen–Fang–Zhang 的结果由 $r^du^\theta\in L_t^pL_x^q$ 及 $2/p+3/q\leq1-d$ , $0\leq d<1$ 的 Prodi–Serrin–Ladyzhenskaya 型条件推出正则; 幂次衰减 $|\Gamma|\leq Cr^\alpha$ 是其中一种充分条件. 本文则关注更接近临界的对数尺度, 不要求环量有正幂次衰减, 代价是必须满足明确的对数小性和强解初值假设. 与 Lei–Ren 的局部定量判据相比, 本文方法可将局部结论中的指数从 $$3/2$$ 改进至 $$5/4$$ , 但其局部阈值依赖 $$G,H$$ 及对数因子, 并非一个与局部临界量无关的普适小常数结论。

摘要和引言也列举了其他类型的结果, 例如 Koch–Nadirashvili–Seregin–Šverák 及 Chen–Strain–Tsai–Yau 对 $|u|\leq C/r$ 类条件的正则性, Pan 对无旋流速度的近轴增长控制, Li–Pan 同时约束 $$b$$ 与 $\Gamma$ 的判据, 以及 Zhang 关于单侧径向速度条件的结果. 这些准则控制的对象、结构和适用设定与本文不同; 作者未将它们逐一作同一量纲或假设层面的定量比较, 因而不能从所提供文本断言本文在所有轴对称判据意义下都更强。

中文标题: 上部快速扩散范围内非局部抛物方程的 Harnack 型不等式

摘要翻译: 我们研究一类形如

(u^q)_t+\mathcal Lu=0,

的双重非线性非局部抛物方程的非负弱解，其中 $\mathcal L$ 是具有对称可测核的非局部 $$p$$ -Laplace 型算子。我们考虑快速扩散范围 $$1<p<q+1$$ ，并在上部快速扩散区域 $\widetilde p_c<p<q+1$ 中建立 Harnack 型估计，其中 $\widetilde p_c$ 是由时间非线性与非局部扩散之间的相互作用确定的临界指数。特别地，这些估计适用于非局部抛物 $$p$$ -Laplace 方程和非局部多孔介质方程。

9. 半线性正则 Euler-Poisson-Darboux-Tricomi 方程的整体存在性，II

Global existence for the semilinear regular Euler-Poisson-Darboux-Tricomi equation, II

作者: Yuequn Li, Fei Guo

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29623v1

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

摘要翻译

中文标题: 半线性正则 Euler-Poisson-Darboux-Tricomi 方程的整体存在性，II

摘要翻译: 在我们此前的工作 [25] 中，我们证明了在阻尼项占主导、强于质量项的情形下，具有幂次非线性 $$|u|^p$$ 的半线性 Euler-Poisson-Darboux-Tricomi 方程的解整体存在。与此相反，本文研究相反的情形，即质量项占主导、强于阻尼项时该方程整体解的存在性。我们证明，若

arXiv ID: 2610.29852v1

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

摘要: We give a finite negative-shift criterion for preserving a prescribed exterior noncollapsing bound under fully nonlinear curvature flows whose hypersurfaces need not be convex. The criterion uses the optimized two-point expression of Andrews and Langford. We verify it by a uniform global two-matrix estimate for the explicit speed

F_{\varepsilon}(A)= tr A- \varepsilon\frac{\norm{A^{\circ}}^{2}}{ tr A}

on the open cone

Γ=\{A\in\Sym(n): tr A>0, \norm{A^{\circ}}< tr A\}.

Here $A^\circ$ denotes the trace-free part of $$A$$ . This cone contains the positive cone and also matrices with a negative eigenvalue, while the displayed speed is concave and non-convex. For an explicit range of parameters depending on the initial exterior bound, the cone and that bound are preserved throughout the smooth embedded existence interval. Together with the known interior estimate for concave speeds, this gives two-sided tangent-ball bounds for this nonlinear family, including closed nonconvex initial hypersurfaces.

8. 振动薄板的瑞利猜想：更简洁的平面证明与非线性推广

Rayleigh’s conjecture for vibrating plates: simpler planar proof and nonlinear extensions

作者: Dorin Bucur, Richard S. Laugesen

类别: math.SP, math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29840v1

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

摘要: We give a simpler proof of Nadirashvili’s theorem (Rayleigh’s conjecture) that among vibrating clamped plates of prescribed area in the plane, the disk has the lowest fundamental tone. The proof uses a one-shot Talenti type comparison for sign-changing solutions of the Poisson equation with mean-zero data, established in a companion paper. The comparison considers just two rearranged sources, rather than optimizing over all admissible sources as in the traditional two-ball technique. This ``paired comparison’’ simplification yields immediate nonlinear extensions: the disk minimizes the clamped quotient $\|Δu\|_q/\|u\|_p$ among bounded open sets of prescribed area, for every $1\le p,q\le\infty$ . Equality is characterized for $$q>1$$ , while exceptional behavior at $$q=1$$ is identified. In higher dimensions, we prove ball minimality for $p=\infty$ and conjecture minimality for all $$p,q$$ within the compact embedding range. We also establish the corresponding interval inequality in one dimension.

9. 具有零均值数据的泊松方程的椭圆重排

Elliptic rearrangement for Poisson equations with mean-zero data

作者: Dorin Bucur, Richard S. Laugesen

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29783v1

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

摘要: The $$L^p$$ norm of a sign-changing Dirichlet solution of Poisson’s equation is shown to be bounded by the maximum of the $$L^p$$ norms of the ball solutions obtained by signed rearrangement of the source and its negative, assuming the source term has mean value zero. This paired comparison on bounded open sets extends Talenti’s elliptic rearrangement principle to solutions that change sign. In dimensions one and two, the paired comparison holds for all $1\le p\le\infty$ , and it holds in every dimension for $p=\infty$ . The inequalities are sharp and equality cases are characterized. In dimensions $n \geq 3$ , paired comparison can fail for every finite $$p$$ although it remains valid for $p\ge(3n-1)/4$ in the special case of bang-bang sources with mean zero. The mean-zero source assumption is crucial when $p<\infty$ . Without it, paired comparison can fail even for bang-bang sources. For sources of arbitrary mean in one dimension, however, the $L^\infty$ paired comparison still holds and so does a comparison for $$L^p$$ norms of the derivative when $1\le p\le3$ .

10. 三维库仑晶格能量的体心立方最优性

BCC optimality for the three-dimensional Coulomb lattice energy

作者: C. B. Muratov, M. Novaga

类别: math-ph, math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29732v1

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

摘要: We prove that the body-centered cubic lattice uniquely minimizes the renormalized Coulomb energy among three-dimensional Bravais lattices of fixed covolume, up to orthogonal transformations. Using Ewald summation, we write the energy, to within an additive constant, as a sum of two positive, rapidly convergent series in real and reciprocal space. Their positivity provides finite lower bounds for all competitors. A passage to the dual lattice reduces the problem to a face-centered cubic comparison. A short basis vector estimate, a certified local Hessian, and interval bounds on a compact five-dimensional domain complete the computer-assisted proof.

11. 任意空间维数下变系数时间分数阶扩散方程中反应系数的恢复

Recovering the reaction coefficient in a time-fractional diffusion equation with variable coefficients in any spatial dimension

作者: D. K. Durdiev

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29763v1

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

摘要: We explore the inverse problem of determining the reaction coefficient $$q(x',t)$$ in the time-fractional diffusion equation ${}^CD_t u-Lu+q(x',t)u=f$ in $\mathbb{R}^n\times(0,T]$ , $$x=(x',x_n)$$ , $$0<α<1$$ , from the trace of the solution on the hyperplane $\{x_n=0\}$ . Here ${}^CD_t$ is the Gerasimov–Caputo derivative and $$L$$ is a uniformly elliptic second-order operator with variable, $$x_n$$ -independent coefficients and without zero-order term. For $$L=Δ$$ this problem was studied in dimensions $$n=1,2,3$$ by methods based on pointwise bounds of the kernels, which are not available for $n\ge4$ . Our approach relies on the Eidelman–Kochubei fundamental solution: although its singularity on the diagonal depends on the dimension, we show that it satisfies integral bounds of the same form for all $n\ge1$ . Using these bounds, we prove existence, uniqueness and stability of regular solutions of the direct problem with a Hölder continuous reaction coefficient. Without assuming additional smoothness of the solution, we show that the inverse problem is equivalent to a closed system of integral equations for $$q$$ and the derivatives $\partial_{x_n}u$ , $\partial^2_{x_n}u$ . For the inverse problem we establish a conditional Lipschitz stability estimate on the whole time interval and global uniqueness in every dimension, and prove local existence of a solution. In the case $a_{in}\equiv0$ , $b_n\equiv0$ we obtain a stability estimate with explicit constants, which reduces to the known estimates for the Laplacian. An explicit two-dimensional example with a mixed coefficient and a normal drift shows that the hypotheses of the main theorems are compatible and that the corresponding terms in the trace condition cannot be omitted.

12. Minkowski 正则化下总等周剖面的单调性与演化

Monotonicity and evolution of total isoperimetric profiles under Minkowski regularization

作者: Zakaria Fattah, Ilias Ftouhi, Giorgio Saracco

类别: math.AP, math.DG

发布日期: 2026-10-21

arXiv ID: 2610.29704v1

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

摘要: We study the evolution of the total isoperimetric profile of a planar convex body under volume-preserving Minkowski regularization. Several scalar geometric quantities are known to be monotone along this flow; we prove that the total isoperimetric profile is monotone as well, and characterize the equality cases. We also establish an explicit reparametrization law, by which, for nonnegative times, the profile of the regularized body is determined by that of the initial one, and the order of profiles of bodies of equal area is preserved. As consequences, we obtain the uniform convergence of the profile to the one of the ball of equal area; its joint $\mathrm C^1$ regularity in the volume and the regularization parameter; the fact that it solves a nonlocal quasilinear Hamilton–Jacobi equation; and sharp inequalities for the profile.

13. 多重调和狄利克雷问题的尖锐布雷齐斯–梅尔不等式：极值测度与临界渐近

Sharp Brezis–Merle Inequalities for Polyharmonic Dirichlet Problems: Extremal Measures and Critical Asymptotics

作者: João Marcos do Ó, Federica Sani, Victor Rafael Silva

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29626v1

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

摘要: We establish the Brezis–Merle inequality for very weak solutions of polyharmonic Dirichlet problems with integrable data on smooth bounded domains in conformal dimensions. The inequality is sharp both with respect to the admissible range of the exponent and to the quantity appearing on the right-hand side, which is characterized as the maximum of an associated one-pole Green functional. The supremum defining this quantity is not attained by any integrable datum, but remains unchanged when the admissible class is extended to Radon probability measures; in this enlarged class, the extremals are Dirac masses, and their concentration points are characterized in terms of the associated Green function. On the unit ball, we prove uniqueness of the extremal measure and obtain an explicit optimal bound extending the sharp result for the Laplacian. Finally, we determine the critical asymptotics of the sharp constant on the right-hand side and show that the possible limiting concentration points are selected by the maximum points of the polyharmonic Robin function.

14. $$L^p$$ -不完备黎曼流形上薛定谔算子的 Liouville 性质与实现

$$L^p$$ -Liouville properties and realizations of Schrödinger operators on incomplete Riemannian manifolds

作者: Andrea Bisterzo, Salvatore Lincastri, Giona Veronelli

类别: math.AP, math.DG

发布日期: 2026-10-21

arXiv ID: 2610.29514v1

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

摘要: The aim of the paper is to highlight and standardize in a self-contained framework the connections between some properties associated with the Laplace-Beltrami operator on possibly incomplete Riemannian manifolds. We present both new and well-known characterizations of the positivity-preserving property, essential self-adjointness, and Liouville-type results for $$L^2$$ -functions, using the latter to establish a clear hierarchy among these properties. Furthermore, we prove the stability of the essential self-adjointness of the Laplace-Beltrami operator on incomplete Riemannian manifolds under the removal of closed (possibly non-compact) subsets from a complete manifold, subject to suitable conditions on their (generalized) Hausdorff measure. With appropriate modifications, this framework is adapted throughout the paper to the general $$L^p$$ -case for $p\in (1,+\infty)$ .

15. $\mathbb{R}^2$ 上 $$L^p$$ -扰动的 Moser–Trudinger 泛函的阈值达到与尖锐渐近行为

Threshold Attainment and Sharp Asymptotics for $$L^p$$ -Perturbed Moser–Trudinger Functionals on $\mathbb{R}^2$

作者: Kaiwen Guo, Yanjun Liu

类别: math.AP, math.FA

发布日期: 2026-10-21

arXiv ID: 2610.29402v1

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

摘要: In this paper, we study the critical Moser–Trudinger supremum

S(λ,p) = \sup_{\substack{ u\in H^1(\mathbb{R}^2)\\ \int_{\mathbb{R}^2} (|\nabla u|^2+u^2)\,dx\le 1 }} \int_{\mathbb{R}^2} \left(e^{4πu^2}-1-λ|u|^p\right)\,dx, \qquad p\ge 2.

Building on the classification away from the thresholds by Chen, Jiang, Lu and Zhu (2026), we resolve the endpoint cases: both thresholds for $$p=2$$ and the finite threshold for $2<p\le4$ are attained. We also prove compactness of centred radial extremals and one-sided value expansions governed by the selected endpoint $$L^p$$ mass. For $$p>4$$ , we quantify how extremals spread out as $λ\to\infty$ by establishing the two-term expansion

S(λ,p)=4π+\mathcal C_pλ^{-2/(p-4)} +D_pλ^{-4/(p-4)}+o\bigl(λ^{-4/(p-4)}\bigr),

Here $\mathcal C_p$ is the maximum of a limiting variational problem at unit $$L^2$$ mass. An explicit correction involving the sixth moment and the full $$H^1$$ normalization determines $$D_p$$ and selects the limiting maximizers. Rescaled centred radial extremals converge strongly to these profiles in $H^1(\mathbb{R}^2)$ along subsequences. For every fixed $$λ>λ_*(4)$$ , we further establish

\lim_{p\downarrow4}(p-4)\log\bigl(S(λ,p)-4π\bigr) =2\log\frac{λ_{\mathrm{loc}}}λ<0,

where $λ_{\mathrm{loc}}$ is the local critical penalty determined by the sharp Gagliardo–Nirenberg constant. The gradient energy has the same logarithmic rate, and the centred radial profiles, normalized in mass and gradient energy, converge along subsequences to sharp Gagliardo–Nirenberg optimizers. An exact transfer identity extends these endpoint and asymptotic results to every constant planar norm.

16. 具有极强退化的热核的上界与下界

Upper and lower bounds for heat kernels with very strong degeneracy

作者: Hector Bouton

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29325v1

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

摘要: We consider the following degenerate parabolic equation on a bounded domain of $\mathbb{R}^d$ containing $$0$$ : $\partial_t u - \nabla\cdot(\varphi(t, x) \nabla u)= 0$ , where $\varphi$ is a rough weight behaving like $$|x|^θ$$ near $$0$$ for some $θ\geq 2$ , and $$u$$ satisfies Neumann boundary condition. We prove that this equation admits a heat kernel with upper and lower Gaussian-type bounds.

17. 自由边界触及退化边界时的尖锐渐近：永续 Heston 看跌期权

Sharp asymptotics for a free boundary touching a degenerate boundary: the perpetual Heston put

作者: Takwon Kim, Junkee Jeon, Jehan Oh

类别: math.AP, math.PR

发布日期: 2026-10-21

arXiv ID: 2610.29253v1

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

摘要: We study the exercise boundary of the perpetual American put in the Heston model near zero variance. There the pricing operator degenerates, and the exercise boundary, a free boundary in the stock–variance half-plane, touches the boundary of the half-plane. No Feller condition is imposed. The boundary is analytic at positive variance and extends as a regular $$C^1$$ curve tangent to the zero-variance axis. Its approach to the axis follows one of three laws, according to whether the Feller ratio $$2κθ/σ^2$$ is larger than, equal to, or smaller than one: a logarithmic correction to linear order, an exponential correction in the square root of the logarithm, or a fractional power. In the first two cases the leading coefficients are explicit and the tangent is not Hölder continuous at the axis. In the third case the exponent is selected by a possibly nonsymmetric stable operator on the axis, and the amplitude is finite and positive.

18. 球面值分数阶调和映射热流的部分正则性

Partial regularity of heat flow of fractional harmonic maps to spheres

作者: Hengrong Du, Tao Huang, Armin Schikorra, Changyou Wang

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29191v1

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

摘要: In this paper, we introduce the notation of stationary heat flow of $$s$$ -harmonic maps from $\R^m$ into the unit sphere $\mathbb{S}^{\ell-1}\hookrightarrow \mathbb{R}^\ell$ , for $$0<s<1$$ , framed through the parabolic version of Caffarelli-Silvestre’s extension. We prove that any stationary heat flow of $$s$$ -harmonic maps is smooth away from a closed set of vanishing $$(m+2-2s)$$ -dimensional parabolic Hausdorff measure. This extends to the fractional settings the partial regularity theorem on the heat flow of harmonic maps into spheres by Chen-Li-Lin and Feldman.

19. 某些拟线性椭圆不等式的孤立奇点 II

Isolated singularities of certain quasi-linear elliptic inequalities II

作者: Shiguang Ma, Shengyang Zang

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29113v1

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

摘要: In this paper, we study the asymptotic behavior of isolated singularities for quasilinear elliptic equations of $$p$$ -Laplacian type. We focus on the $$p$$ -Laplace equation with the critical Serrin exponent $$n(p-1)/(n-p)$$ and the $$n$$ -Laplace equation with critical exponential growth. We prove that any positive solution either admits a $C^{1,α}$ extension to the origin or exhibits a precise singular asymptotic behavior. These results provide a complement to our previous study on $$p$$ -Laplacian type inequalities.

20. $$k$$ -Hessian 方程的孤立奇点

Isolated Singularities of k-Hessian Equations

作者: Shiguang Ma, Shengyang Zang

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29079v1

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

摘要: We classify isolated singularities of nonpositive, $$k$$ -admissible classical solutions to $$k$$ -Hessian equations without radial symmetry assumptions. For $1\le k<n/2$ , we consider $$F_k[u]=(-u)^p$$ with $1<p\le nk/(n-2k)$ . Below the Serrin exponent, we obtain fundamental-solution asymptotics and characterize removability by the vanishing of the leading coefficient. At the Serrin exponent, we determine the exact logarithmic correction and its universal coefficient for every nonremovable singularity. In the critical dimension $$k=n/2$$ , we establish logarithmic asymptotics and a removability criterion for $F_k[u]=f(x,u,\nabla u)$ with continuous sources satisfying $0\le f(x,s,ξ)\le C|ξ|^\ell e^{-αs}$ , where $$α>0$$ and $0\le\ell<n$ . The proofs combine Hessian-measure compactness, Wolff-potential estimates, a global Bôcher theorem, and a bounded-variation analysis of spherical extrema.

21. SIS流行病模型的时间周期解

Time-periodic solutions of an SIS epidemic model

作者: Fang Li, Pengfei Song, Kexin Zhang, Maolin Zhou

类别: math.AP

发布日期: 2026-10-21

arXiv ID: 2610.29062v1

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

摘要: We prove the existence of strictly positive, nonstationary time-periodic solutions for the spatially heterogeneous SIS epidemic model with standard incidence. On a bounded interval with Neumann boundary conditions, we construct smooth, time-independent transmission and recovery rates satisfying $$β(x)>γ(x)>0$$ , together with diffusion coefficients $$d_I>d_S>0$$ , such that periodic solutions coexist with the unique endemic equilibrium for every prescribed positive total population. Thus, pointwise dominance of transmission over recovery does not ensure convergence to equilibrium. Our approach combines an inverse construction of an unstable endemic equilibrium with a coefficient path that preserves both the equilibrium and the pointwise inequalities. A change in the number of unstable eigenvalues along this path yields small-amplitude periodic solutions through a center-manifold reduction and a generalized Hopf bifurcation argument. The proof uses an AI-assisted construction of an unstable endemic equilibrium under the pointwise condition $$β>γ>0$$ .

22. 用于能量核最大均值差异的跳跃修正 Wasserstein 梯度流

A Jump-Modified Wasserstein Gradient Flow for Maximum Mean Discrepancy with Energy Kernels

作者: Ansh Desai, Jason Liu, Srujan Roplekar, Matthew Rosenzweig

类别: math.AP, cs.LG, math.PR, stat.ML

发布日期: 2026-10-21

arXiv ID: 2610.29007v1

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

摘要: We study a jump-modified Wasserstein gradient flow of the squared maximum mean discrepancy generated by the energy kernels $$K(z)=-|z|^q$$ , $$0<q<2$$ , in every dimension $d\ge2$ as well as the one-dimensional range $1\le q<2$ . Under subcritical density and finite-moment assumptions, we prove global well-posedness. The added jumps transfer mass toward smaller values of the first variation and dissipate the same MMD energy as the transport term. Under our hypotheses on the target and jump rates, every Borel probability measure with finite MMD energy and zero jump dissipation is equal to the target. Combining this rigidity with compactness and a second-moment identity yields convergence to the target from finite-second-moment data, without additional uniform-in-time density or moment assumptions. We also construct a global collision-free particle process whose energy is nonincreasing along each path. A quantitative mean-field estimate applies to general initial empirical approximations, including deterministic configurations; propagation of chaos for exchangeable data and long-time target approximation follow as consequences.

23. 强梯度吸收与初始衰减防止爆破

Strong gradient absorption and initial decay prevent blow-up

作者: Philippe Souplet, Berikbol T. Torebek

类别: math.AP

发布日期: 2026-10-20

arXiv ID: 2610.28913v1

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

摘要: Consider positive classical solutions of the Cauchy problem for $u_t-Δu=u^p-μ|\nabla u|^q$ , with sufficiently smooth initial data $$u_0$$ , in the strong absorption regime, namely $q\ge p>1$ and $$μ>0$$ . It has remained a long standing open problem whether finite time blow-up can occur if $$u_0$$ decays at space infinity. We here solve this problem, showing that the answer in negative. To this end, we establish a novel, universal Bernstein-type gradient estimate, of independent interest, valid only in the range $q\ge p$ . Moreover, the decay assumption on $$u_0$$ turns out to be essentially optimal. Indeed we show that the solution blows up in finite time whenever $$u_0$$ has a positive limit at infinity, and we determine the blow-up behavior of $$u$$ at space infinity under an additional assumption.

24. 欧几里得半空间中 Morse 指数为一的自由边界极小超曲面

Free boundary minimal hypersurfaces of Morse index one in Eucildean half-space

作者: Junzhang Li

类别: math.DG, math.AP

发布日期: 2026-10-20

arXiv ID: 2610.28824v1

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

摘要: We prove that a complete, connected, two-sided, embedded free boundary minimal hypersurface of finite total curvature and Morse index one in a half-space of RN, with nonempty boundary is a higher-dimensional half-catenoid. We also prove that a complete, connected, two-sided, embedded free boundary minimal hypersurface of finite total curvature and Morse index one in a quarter-space is the higher-dimensional quarter-catenoid.

25. Hong–Krahn–Szegő 不等式的尖锐定量稳定性

Sharp quantitative stability of the Hong–Krahn–Szegő inequality

作者: Yuxuan Zhou, Wenming Zou

类别: math.AP

发布日期: 2026-10-20

arXiv ID: 2610.28757v1

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

摘要: We establish sharp quantitative stability estimates for the Hong–Krahn–Szegő inequality for arbitrary open sets of positive finite measure in $\mathbb R^d$ , $d\ge2$ . In dimensions $d\ge3$ , the optimal stability exponent is $$(d+2)/2$$ , while in dimension two the sharp stability modulus is $t^2/\log(e/t)$ . These results improve the quantitative stability estimate of Brasco and Pratelli~\cite[Theorem 3.5]{BP} and resolve Open Problem 7.56 of Brasco and De Philippis~\cite[Open Problem 7.56]{BDP}.

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

1. 二维非齐次不可压 Euler 方程的适定性理论及一类全局时间解#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 具有支撑延伸至边界的涡量扰动的库埃特流的非线性无黏阻尼#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. PDE理论引导的神经算子与二维各向异性 Boussinesq 系统的长期弛豫#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 黏性 Gröbli 三角形的长时间动力学#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 涡丝的 Arnold 型稳定性分析#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 黏性与热传导理想气体径向对称定常态的无黏极限问题#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 轴对称 Navier–Stokes 方程的一个对数型正则性判据#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 通过纳什嵌入研究紧黎曼流形上的平方梯度非线性饱和空间#

摘要翻译#

2. 紧致正规齐性空间上双线性薛定谔方程的小时间近似可控性#

摘要翻译#

3. 带加性噪声的周期随机 Korteweg-de Vries 方程无条件适定性的正规形方法#

摘要翻译#

4. 零矩约束下 CR 球面上的锐利 Moser-Onofri 不等式及其在 CR Moser-Onofri 不等式稳定性中的应用#

摘要翻译#

5. 黏性-色散热传导可压缩流体一般黎曼解的时间渐近稳定性#

摘要翻译#

6. 用于对流-扩散问题的超收敛准间断 Galerkin 方法#

摘要翻译#

7. 高秩谱投影算子估计与 Stein–Tomas 限制定理#

摘要翻译#

8. 上部快速扩散范围内非局部抛物方程的 Harnack 型不等式#

摘要翻译#

9. 半线性正则 Euler-Poisson-Darboux-Tricomi 方程的整体存在性，II#

摘要翻译#

10. 哈代算子的最大速度估计#

摘要翻译#

11. 辐射气体模型系统的局部适定性、爆发现象与不适定性#

摘要翻译#

12. Schwarzschild 黑洞背景下非线性 Dirac 系统的旋量适配几何方法#

摘要翻译#

13. 圆周上 Hamilton 偏微分方程的有效稳定性与尖锐对数时间尺度#

摘要翻译#

14. 关于具有应力扩散的二维不可压缩流体的有限维吸引子#

摘要翻译#

15. 非线性薛定谔方程时空锥中的散射#

摘要翻译#

16. 紧支撑聚集-扩散方程的尖锐定量收敛#

摘要翻译#

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

1. 关于流与时空轨迹的严格逼近#

2. 库仑 Landau 方程的边界驱动非平衡态#

3. 关于一类全非线性奇异退化椭圆方程的 ABP 估计及其应用#

4. Hartree 型能量下球的刚性及基于连续性论证的二阶形状导数#

5. $S^2\times S^2$ 的连通和上爱因斯坦度量的一种变分构造#

6. 通过 Stone-Weierstrass 型定理实现紧流形上的饱和，以及受迫非线性偏微分方程#

7. 非凸超曲面全非线性曲率流的有限平移逆凹性与外部非塌缩#

8. 振动薄板的瑞利猜想：更简洁的平面证明与非线性推广#

9. 具有零均值数据的泊松方程的椭圆重排#

10. 三维库仑晶格能量的体心立方最优性#

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

1. 二维非齐次不可压 Euler 方程的适定性理论及一类全局时间解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 具有支撑延伸至边界的涡量扰动的库埃特流的非线性无黏阻尼

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. PDE理论引导的神经算子与二维各向异性 Boussinesq 系统的长期弛豫

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 黏性 Gröbli 三角形的长时间动力学

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 涡丝的 Arnold 型稳定性分析

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 黏性与热传导理想气体径向对称定常态的无黏极限问题

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 轴对称 Navier–Stokes 方程的一个对数型正则性判据

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 通过纳什嵌入研究紧黎曼流形上的平方梯度非线性饱和空间

摘要翻译

2. 紧致正规齐性空间上双线性薛定谔方程的小时间近似可控性

摘要翻译

3. 带加性噪声的周期随机 Korteweg-de Vries 方程无条件适定性的正规形方法

摘要翻译

4. 零矩约束下 CR 球面上的锐利 Moser-Onofri 不等式及其在 CR Moser-Onofri 不等式稳定性中的应用

摘要翻译

5. 黏性-色散热传导可压缩流体一般黎曼解的时间渐近稳定性

摘要翻译

6. 用于对流-扩散问题的超收敛准间断 Galerkin 方法

摘要翻译

7. 高秩谱投影算子估计与 Stein–Tomas 限制定理

摘要翻译

8. 上部快速扩散范围内非局部抛物方程的 Harnack 型不等式

摘要翻译

9. 半线性正则 Euler-Poisson-Darboux-Tricomi 方程的整体存在性，II

摘要翻译

10. 哈代算子的最大速度估计

摘要翻译

11. 辐射气体模型系统的局部适定性、爆发现象与不适定性

摘要翻译

12. Schwarzschild 黑洞背景下非线性 Dirac 系统的旋量适配几何方法

摘要翻译

13. 圆周上 Hamilton 偏微分方程的有效稳定性与尖锐对数时间尺度

摘要翻译

14. 关于具有应力扩散的二维不可压缩流体的有限维吸引子

摘要翻译

15. 非线性薛定谔方程时空锥中的散射

摘要翻译

16. 紧支撑聚集-扩散方程的尖锐定量收敛

摘要翻译

不相关论文（基本信息）

1. 关于流与时空轨迹的严格逼近

2. 库仑 Landau 方程的边界驱动非平衡态

3. 关于一类全非线性奇异退化椭圆方程的 ABP 估计及其应用

4. Hartree 型能量下球的刚性及基于连续性论证的二阶形状导数

5. $S^2\times S^2$ 的连通和上爱因斯坦度量的一种变分构造

6. 通过 Stone-Weierstrass 型定理实现紧流形上的饱和，以及受迫非线性偏微分方程

7. 非凸超曲面全非线性曲率流的有限平移逆凹性与外部非塌缩

8. 振动薄板的瑞利猜想：更简洁的平面证明与非线性推广

9. 具有零均值数据的泊松方程的椭圆重排

10. 三维库仑晶格能量的体心立方最优性