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

1. 欧拉方程中从熵波到涡片

From entropy waves to vortex sheets for the Euler equations

作者: Yuri Trakhinin, Wei Wang, Zhifei Zhang, Wenbin Zhao

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32893v1

摘要翻译: 我们构造了全可压缩欧拉方程的分片光滑解，这些解从熵波演化为涡片，同时保持严格为正的 Taylor 跳跃。对于每个整数 $κ\geq1$ ，我们得到满足压力和速度的初始物质相容性直至 $$κ$$ 阶的解，其切向速度跳跃具有一个非零的 $t^{κ+1}$ 阶首项。这些解关于切向变量是周期的，并且可在具有不可穿透壁面的薄层区域中光滑实现。我们还构造了局部非周期切向区域上的解，这些解满足任意阶的初始物质相容性，但其切向速度跳跃在每个充分小的正时间均非零，并且在 $$t=0$$ 处是平坦的。有限阶构造采用速度势的法向 Cauchy-Kowalevski 问题，而全阶构造则通过齐次化约化为一个光滑双曲 Cauchy 问题。

Comment: 14 pages. All comments are welcome

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

1. 研究对象和背景

本文研究三维全可压缩 Euler 方程中接触间断的生成机制, 特别考察初始时刻具有连续压力和连续全速度的熵波, 是否会在正时间演化为带有切向速度跳跃的涡片. 方程组为

\begin{cases} \partial_t\rho+\nabla\cdot(\rho u)=0,\\ \partial_t(\rho u)+\nabla\cdot(\rho u\otimes u)+\nabla p=0,\\ \partial_t E+\nabla\cdot((E+p)u)=0, \end{cases} \qquad p=A\rho^\gamma e^S,\quad E=\frac{p}{\gamma-1}+\frac12\rho|u|^2,

其中 $\rho,u,p,E,S$ 分别为密度、速度、压力、总能量密度和熵, $A>0,\gamma>1$ . 在每个相内, 对正密度和正压力, 可等价使用

D_t\rho=-\rho\nabla\cdot u,\qquad D_tu=-\rho^{-1}\nabla p,\qquad D_tp=-\gamma p\nabla\cdot u, \qquad D_t=\partial_t+u\cdot\nabla.

自由界面写为 $\Gamma_f=\{y=f(t,x)\}$ , 两侧区域分别记为 $\Omega_\pm(t)$ , 法向量未归一化形式为 $N=(-\nabla_x f,1)$ . 接触间断的相对法向质量通量为零, 因而 Rankine–Hugoniot 条件化为

[p]=0,\qquad u^+\cdot N=u^-\cdot N=f_t.

本文所说的熵波是 $$[u]=0$$ 而 $[\rho]\ne0$ 的接触间断; 当界面上的切向速度跳跃 $[u_\tau]$ 非零时, 则出现涡片. 因此问题并非一般意义上的已存在涡片稳定性, 而是从初始熵波状态出发, 研究界面两侧速度何时、以何种时间阶失去连续性. 论文关注的初始兼容条件是压力和全速度沿两侧物质导数的迹匹配:

[(D_t^j p)_0]=0,\qquad [(D_t^j u)_0]=0,\qquad 0\le j\le\kappa.

这些条件由相内 Euler 方程递归确定, 是若要经典地延续压力及全速度连续性所必需的初始条件. 核心问题是: 任意给定有限阶, 甚至在光滑范畴中所有阶的初始物质兼容性, 加上严格正的 Taylor 跳跃, 是否足以迫使接触间断继续是熵波? Taylor 跳跃取为

T=n\cdot[\nabla p]=\frac{[\nabla_Np]}{|N|},\qquad \nabla_Np=N\cdot\nabla p.

它的正性对应自由边界问题中的 Taylor 符号条件, 通常与界面演化的良定性结构相关. 作者构造的解表明, 正 Taylor 跳跃并不阻止熵波转化为涡片.

这一机制与文献 [13] 中关于熵波的分析直接相关. 若全速度连续性在一段时间内保持, 沿共同物质界面轨迹对动量关系作切向投影, 会得到

[\rho^{-1}]\nabla_{\Gamma_f}q=0,

其中 $$q$$ 为共同压力迹. 因 $[\rho]\ne0$ , 维持速度连续要求 $$q$$ 在每个连通界面分支上切向常数. 若 $$q$$ 有切向变化, 两侧压力加速度因 $\rho^{-1}$ 不同而不同, 因而可产生切向滑移. 本文把这一机制转化为精确的局部构造, 并进一步研究压力迹在初始时刻具有高阶消失或无限阶平坦时的影响.

具体地, 有限阶构造在切向周期变量 $x\in\mathbb T^{d-1}$ 上进行, 并先在界面邻域构造, 再延拓到带不可穿透壁面的条形区域 $\mathbb T^{d-1}\times(-L,L)$ . 其界面为加速平面 $$f(t)=bt^2/2$$ , 压力迹为 $q=P_*+\varepsilon t^\kappa a(x)$ , 其中 $$a$$ 为非恒定解析函数. 无限阶构造则在局部非周期切向区间 $I\subset(0,\infty)$ 上进行, 采用曲界面 $$f(t,x)=bt^2/(2x)$$ 和平坦扰动压力迹 $P_*+\varepsilon e^{-x^2/t^2}$ ( $$t>0$$ ), 并令其在 $$t=0$$ 等于 $$P_*$$ . 两种构造分别回答有限阶兼容和光滑无限阶兼容下的生成问题; 后者说明, 即使所有初始物质导数匹配, 光滑但非解析的时间平坦扰动仍可能在每个充分小的正时间产生非零滑移。

2. 主要定理或主要结果

主要结论为定理 2.2 和定理 2.3.

定理 2.2 处理有限阶兼容. 给定整数 $\kappa\ge1$ , $a\in C^\omega(\mathbb T^{d-1};\mathbb R)$ 且 $\nabla_xa\not\equiv0$ , 以及常数 $$b,P_*>0$$ , $$r_->r_+>0$$ , 存在 $\varepsilon_0,T,\delta>0$ , 使每个 $|\varepsilon|\le\varepsilon_0$ 对应一个定义在界面邻域

0\le t\le T,\qquad 0<\pm(y-f(t))<\delta,\qquad x\in\mathbb T^{d-1},\qquad f(t)=\frac12bt^2

上的接触解. 两相场在界面处光滑, 且光滑依赖于 $\varepsilon$ . 界面共同压力迹为 $q(t,x)=P_*+\varepsilon t^\kappa a(x)$ , 法向速度为 $$bt=f_t$$ , 并且 $\rho^\pm=R_\pm(p^\pm)$ , 其中

R_\pm(p)=r_\pm\left(\frac p{P_*}\right)^{1/\gamma}.

相内熵为常数, 两侧熵不同. 初始数据满足 $$[(D_t^jp)_0]=[(D_t^ju)_0]=0$$ 对 $0\le j\le\kappa$ . 但是当 $\varepsilon\ne0$ 时, 第一个非零切向速度物质导数差出现在 $\kappa+1$ 阶:

[(D_t^{\kappa+1}U)_0] =-\varepsilon\kappa!\left(\frac1{r_+}-\frac1{r_-}\right)\nabla_xa.

Taylor 跳跃在界面上满足

T=b\big(R_-(q)-R_+(q)\big)\ge\theta>0,

其中 $\theta$ 可取与 $\varepsilon$ 无关. 对每个固定非负整数 $$s$$ , 切向速度跳跃具有展开

[U](t)=-\frac{\varepsilon t^{\kappa+1}}{\kappa+1} \left(\frac1{r_+}-\frac1{r_-}\right)\nabla_xa+R_\varepsilon(t), \qquad \|R_\varepsilon(t)\|_{H^s}\le C_s|\varepsilon|t^{\kappa+2}.

因此缩短 $$T$$ 后存在 $$c_s>0$$ , 使

\|[U](t)\|_{H^s(\mathbb T^{d-1})}\ge c_s|\varepsilon|t^{\kappa+1},\qquad 0<t\le T.

初始体场还满足, 对每个固定整数 $m\ge0$ ,

\sum_\pm\|U^\pm_{0,\varepsilon}-(\bar\rho^\pm,0,\bar p^\pm)\|_{H^m(\Omega^\pm_{0,\delta})}\le C_m|\varepsilon|,

其中参考熵波剖面为 $\bar p^\pm(\xi)=P_*(1-\beta_\pm\xi)^{\gamma/(\gamma-1)}$ , $\bar\rho^\pm(\xi)=r_\pm(1-\beta_\pm\xi)^{1/(\gamma-1)}$ , $\xi=y-bt^2/2$ . 对任意预先给定的条形宽度 $$L>0$$ , 可把局部构造平滑延拓为 $\mathbb T^{d-1}\times(-L,L)$ 上的分片 $C^\infty$ 解, 在 $y=\pm L$ 满足 $$u_3=0$$ , 且上述滑移估计仍成立. 寿命和常数可以依赖 $\kappa$ . $\varepsilon=0$ 时得到持续为熵波的参考解; $\varepsilon\ne0$ 时, 任意指定有限 Sobolev 范数下初始扰动可为 $O(|\varepsilon|)$ , 但正时间立即产生非零切向滑移.

定理 2.3 给出光滑无限阶兼容的局部构造. 固定 $\gamma>1$ , $$b,P_*>0$$ , $$r_->r_+>0$$ , $\varepsilon>0$ 和紧区间 $I\subset(0,\infty)$ , 存在 $T,\delta>0$ 及局部的分片 $C^\infty$ 解. 解在每一相内等熵且无旋, 界面为

f(t,x)=\frac{bt^2}{2x},\qquad x\in I,

初始迹为 $u^\pm_0=0$ , $p^\pm_0=P_*$ , $\rho^\pm_0=r_\pm$ , 并且所有阶的初始物质兼容条件均成立. Taylor 跳跃在 $0\le t\le T$ 上一致正; 初始时

T(0,x)=\frac{b(r_--r_+)}x.

以切向正向为 $$x$$ 增大的单位切向量时, 滑移满足

[u_\tau](t,x,f(t,x)) =\varepsilon\left(\frac1{r_+}-\frac1{r_-}\right) \frac{t}{x}e^{-x^2/t^2}(1+o(1)),\qquad t\downarrow0,

且 $$o(1)$$ 对 $x\in I$ 一致. 该跳跃对每个充分小的 $$t>0$$ 非零, 但在 $$t=0$$ 平坦到任意阶. 共同压力迹为 $P_*+\varepsilon e^{-x^2/t^2}$ ( $$t>0$$ ), 在 $$t=0$$ 取 $$P_*$$ . 取 $\varepsilon=0$ 得到同一界面上的参考熵波, 其压力迹恒为 $$P_*$$ 并持续保持速度连续。定理 2.3 是局部非周期构造; 摘录未声称其可延拓到带壁面的整个条形区域。

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

证明的共同策略是先指定界面运动和共同压力迹, 再由运动学条件及势函数关系构造两侧界面速度数据, 最后分别实现两相的体解. 这样可直接控制希望生成滑移的界面历史, 避免先任意选取体初值、再反向追踪自由边界的困难. 两个定理的体解实现方式不同, 且论文明确指出两种证明彼此独立.

有限阶构造中, 取加速坐标 $\xi=y-f(t)$ , 并使用速度势表示

u^\pm=bt\,e_y+\nabla_{x,\xi}\psi^\pm,\qquad p^\pm=H_\pm^{-1}\left(-\psi_t^\pm-\frac12|\nabla\psi^\pm|^2-b\xi\right),\qquad \rho^\pm=R_\pm(p^\pm),

其中 $H_\pm'(p)=R_\pm(p)^{-1}$ , 且 $H_\pm(P_*)=0$ . 先在界面上解 Hamilton–Jacobi 方程

\partial_tA^\pm+\frac12|\nabla_xA^\pm|^2+H_\pm(q)=0, \qquad A^\pm(0,x)=0,

并令 $\psi^\pm(t,x,0)=A^\pm(t,x)$ , $\psi^\pm_\xi(t,x,0)=0$ . 第一项数据给出切向速度, 第二项数据配合 $$bt$$ 给出共同法向速度; Hamilton–Jacobi 方程保证两侧压力迹都等于预设的 $$q$$ . 由于 $$a$$ 解析, 这些界面数据可由解析 Cauchy 定理构造.

关键障碍是要从界面迹恢复满足完整 Euler 方程的体场, 并确认法向方向确实是可用的 Cauchy 方向. 由势函数 Bernoulli 关系和压力演化方程得到拟线性二阶势方程

\psi_{tt}+2\nabla\psi\cdot\nabla\psi_t +\nabla\psi\cdot(D^2\psi)\nabla\psi-c^2\Delta\psi+b\psi_\xi=0.

令 $w=\psi_\xi$ , $z=\nabla_x\psi$ , 可将其整理为以 $\xi$ 为演化变量的正规形

(c^2-w^2)\psi_{\xi\xi} =\psi_{tt}+2w\psi_{\xi t}+2z\cdot\nabla_x\psi_t +2wz\cdot\nabla_x\psi_\xi +z\cdot(D_x^2\psi)z-c^2\Delta_x\psi+bw.

界面数据上 $$w=0$$ , $c^2=\gamma q/R_\pm(q)>0$ , 故系数 $$c^2-w^2$$ 非零. 这正是可在法向变量应用二阶 Cauchy–Kowalevski 定理的条件. 解析唯一性使局部构造可在环面坐标图间拼接, 紧致性给出共同管状邻域. 随后通过 Bernoulli 关系验证动量方程, 势方程给出压力方程, $R_\pm'(p)=R_\pm(p)/(\gamma p)$ 给出质量方程; 因而所得并非仅满足某个等熵子系统的形式解, 而是完整可压缩 Euler 解. 法向压力导数满足 $R_\pm(q)^{-1}\partial_\xi p^\pm=-b$ , 从而获得正 Taylor 跳跃.

滑移的产生由界面 Hamilton–Jacobi 方程的切向梯度直接读出:

\partial_tv^\pm+(v^\pm\cdot\nabla_x)v^\pm =-\frac{\varepsilon t^\kappa}{R_\pm(q)}\nabla_xa.

递归展开给出

v^\pm(t,x)=-\frac{\varepsilon t^{\kappa+1}}{(\kappa+1)r_\pm}\nabla_xa(x) +\varepsilon t^{\kappa+2}B^\pm(t,x,\varepsilon),

其中余项解析且在固定有限阶切向导数下有界. 两侧压力加速度的密度因子不同, 因而首个速度差恰在 $t^{\kappa+1}$ 阶出现. 兼容性验证使用界面迹上的物质导数 $L_\pm=\partial_t+v^\pm\cdot\nabla_x$ . 因 $v^\pm=O(t^{\kappa+1})$ , 对 $j\le\kappa+1$ 有 $(L_\pm^jG)(0)=\partial_t^jG(0)$ ; 结合压力迹和速度迹的 Taylor 展开即可逐阶验证条件, 并得到首个不匹配阶. 对任意固定 $$s$$ , 展开余项为 $O(|\varepsilon|t^{\kappa+2})$ 后, 用反三角不等式推出滑移的 $$H^s$$ 下界.

定理 2.2 的壁面实现还需要处理局部构造与全区域的差异. 作者把管状邻域中的初始场用截断函数延拓到全条形区域, 保持正密度和正压力, 再对两相分别使用对称双曲 Euler 系统的局部适定性. 通过有限传播速度和局部唯一性, 截断影响在短时间内到不了较小界面邻域及壁面附近, 因而局部精确接触构造保留, 同时壁面处速度法向分量为零. 对初值关于 $\varepsilon$ 的 $O(|\varepsilon|)$ 控制来自联合解析依赖及截断操作.

无限阶构造不能直接照搬有限阶的解析法向 Cauchy–Kowalevski 方法: $e^{-x^2/t^2}$ 在 $$t=0$$ 光滑平坦, 却非解析, 正是要保留的关键时间行为. 作者改用二维相似变量 $$s=t/x$$ , $$z=y/x$$ 和齐次势函数

\phi^\pm(t,x,y)=xF^\pm(s,z).

势方程因此严格化为

AF_{ss}+2BF_{sz}+CF_{zz}=0,

其中系数 $$A,B,C$$ 由速度和声速构成. 取界面 $$z=g(s)=bs^2/2$$ 后, 压力迹是 $Q(s)=P_*+\varepsilon\chi(s)$ , $\chi(s)=e^{-1/s^2}$ 对 $$s>0$$ 、 $\chi(s)=0$ 对 $s\le0$ . 通过界面势迹 $a^\pm(s)=F^\pm(s,g(s))$ 和接触条件确定 Cauchy 数据. 核心不是把平坦数据解析化, 而是利用约化方程在参考状态附近、以横向坐标为演化变量仍为双曲型这一性质, 在光滑范畴中解 Cauchy 问题. 从两侧边界势之差的方程提取切向滑移, 得到 $t e^{-x^2/t^2}/x$ 量级. 因此此法既实现非解析平坦压力迹, 又保留精确界面条件和正 Taylor 跳跃。

4. 与之前工作的比较

本文与既有工作的差别首先在研究问题上. 既有涡片稳定性文献通常以已经存在的非零切向速度跳跃为背景态, 研究其扰动是否稳定. 文中提到 Coulombel–Secchi 对二维等熵涡片建立线性估计和局部非线性存在性结果 [3,4]; Morando–Trebeschi 对非等熵流建立相应线性估计 [10], Morando、Trebeschi 和 Wang 证明适当超声速假设下的非线性稳定性 [11]. 这些结果处理的是涡片的延续或稳定, 且非线性论证需要 Nash–Moser 迭代以弥补线性估计中的导数损失. 本文则从初始切向速度连续、密度有跳跃的熵波出发, 构造速度连续性实际破坏的解; 它不以涡片稳定性为结论, 也不宣称解决一般扰动下的自由边界适定性.

其次, 本文把此前关于熵波的结构分析推进到显式的生成构造. 作者在 [13] 中研究 Euler 坐标下的熵波, 通过切向界面导数的演化方程把 Taylor 符号条件与双曲性联系起来, 并得到无导数损失的先验估计; 那些估计针对压力和全速度在界面两侧保持连续的解. [13] 还讨论了保持连续性所需的兼容条件及其失效可能性. 本文不推翻熵波持续解的先验理论, 而是构造另一类在初始时刻满足给定兼容条件、随后不再保持速度连续的解, 从而指出不能仅凭初始兼容性和正 Taylor 跳跃推出熵波性质持续. 文中给出的机制与关系 $[\rho^{-1}]\nabla_{\Gamma_f}q=0$ 一致: 维持速度连续会约束共同压力迹切向常数, 而非恒定压力迹使两侧切向压力加速度不同.

在兼容阶数方面, 有限阶定理对任意预先指定的 $\kappa\ge1$ 都能令压力与全速度的初始物质导数匹配到 $\kappa$ 阶, 同时把首个切向速度缺陷精确推迟到 $\kappa+1$ 阶, 且滑移主项为 $t^{\kappa+1}$ . 这比只展示某个低阶不兼容例子更精细, 因为它给出兼容阶数与涡片生成时间阶之间的明确对应, 并提供任意固定 $$H^s$$ 范数中的余项控制和下界. 无限阶构造进一步表明, 光滑范畴的所有初始物质 jets 相同仍不保证时间延续; 滑移可为非零但对 $$t=0$$ 平坦. 这一结论依赖光滑非解析函数, 不应误读为解析范畴中也存在相同的无限阶反例; 摘录没有证明解析初值的相应结论.

在 Taylor 条件和扰动大小方面, 两族构造都保持严格正的 Taylor 跳跃. 有限阶情形中该下界可选为与 $\varepsilon$ 无关, 而初始数据在每个固定有限阶 Sobolev 范数下与参考熵波相差 $O(|\varepsilon|)$ . 因此破坏熵波性质并不要求先失去 Taylor 符号条件, 也不要求在指定有限阶范数中远离参考波. 但这不是对所有范数一致的小扰动结论: 常数可依赖所考察的 Sobolev 阶数, 寿命和常数可依赖 $\kappa$ , 不能据此断言存在统一控制所有阶范数的单一小扰动尺度.

在空间范围和边界设置方面, 定理 2.2 的有限阶构造从周期切向变量出发, 并进一步通过截断、局部双曲演化和有限传播速度实现到任意预给定宽度的有壁条形区域, 壁面满足不可穿透条件. 定理 2.3 的无限阶构造只在局部非周期切向区间 $I\subset(0,\infty)$ 和界面邻域内成立, 并采用 $$f=bt^2/(2x)$$ 的曲界面; 摘录没有把它延拓到周期域或带固定壁面的完整条形区域. 因而无限阶结果在兼容性阶数上更强, 但全局空间实现范围比有限阶结果窄.

方法上的代价也不同. 有限阶结果要求切向振幅 $$a$$ 为解析函数, 借助以法向变量为 Cauchy 变量的 Cauchy–Kowalevski 理论构造体势; 它能获得周期性及有壁条形区域的实现, 但解析性是这一构造框架的重要条件. 无限阶结果允许所需的平坦非解析时间迹, 依赖齐次约化后横向演化的双曲性, 代价是局部化、特定相似结构和曲界面, 不能直接视为一般自由界面的统一构造定理.

此外, 接触波的耗散稳定性是另一条相关但不同的研究线. 文中列举 Huang、Matsumura 和 Xin 关于一维可压缩 Navier–Stokes 弱黏性接触波在零超额质量条件下渐近稳定的结果 [6], Huang、Xin 和 Yang 对一般小扰动的 Navier–Stokes 与 Boltzmann 情形的工作 [7], 以及 Kang、Vasseur 和 Wang 对三维平面熵波在 Navier–Stokes–Fourier 零耗散极限类中的稳定性与唯一性结果 [8]. 这些工作包含黏性或零耗散极限框架, 与本文的无黏全 Euler、自由界面涡片生成构造并非同一问题; 本文不从耗散极限推出涡片, 也不提供黏性稳定性结论. 文中对这些前人工作只作问题背景区分, 未给出可直接比较的统一参数阈值、最优时间尺度或定量常数, 因而此类更细的性能比较从提供内容中无法确认。

2. gSQG 补丁方程的自相似解

Self-similar solutions for the gSQG patch equations

作者: Javier Gomez-Serrano, Marcus Pasquariello

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32857v1

摘要翻译: 无黏性广义表面准地转（gSQG）方程是一族由标量 $$α$$ 参数化的主动标量方程，其中包括二维 Euler 方程（ $$α= 0)$$ ）和 SQG 方程 $$(α= 1)$$ 。我们证明，对每个 $α\in (\frac{7}{4}, 2)$ 和每个 $$T > 0$$ ，都存在一个自相似补丁解，其边界在 $0 \leq t < T$ 内是一条光滑的无界图形，并在时间 $$T$$ 形成角点奇异性。

Comment: 77 pages

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

1. 研究对象和背景

本文研究无黏广义表面准地转方程（gSQG）的斑块解。其主动标量模型为

\partial_t\theta+u\cdot\nabla\theta=0,\qquad u=-\nabla^\perp(-\Delta)^{-1+\alpha/2}\theta,\qquad (x,t)\in\mathbb R^2\times\mathbb R_+,

其中 $\alpha\in(0,2)$ 。 $\alpha=0$ 对应二维 Euler 方程， $\alpha=1$ 对应 SQG 方程；在 $\alpha\in(1,2)$ 时，速度场相对于标量具有更强的奇异性。斑块设定取 $\theta(\cdot,t)=\mathbf 1_{\Omega(t)}$ ，于是标量的输运等价于区域边界的演化。论文考虑无界下方区域

\Omega(t)=\{(x_1,x_2)\in\mathbb R^2:x_2<F(x_1,t)\},

其边界由图函数 $$F$$ 描述。相应的前沿方程为

\partial_tF(x,t)=c_\alpha\int_{\mathbb R}\frac{\partial_xF(x,t)-\partial_xF(x',t)}{\big((x-x')^2+(F(x,t)-F(x',t))^2\big)^{\alpha/2}}\,dx',

其中 $c_\alpha>0$ 。当 $\alpha\in(1,2)$ 且 $\partial_xF$ 有界并满足 Lipschitz 正则性时，积分绝对收敛，不要求 $$F$$ 在无穷远衰减。平坦界面附近的线性化（按文中规范化）是色散方程

\partial_tF=|D|^{\alpha-1}\partial_xF,

其色散关系为 $\xi|\xi|^{\alpha-1}$ 。因此，问题并非仅是几何曲率失控：界面的非局部色散与非线性共同决定了奇异尺度及其形成方式。

研究动机来自斑块边界的有限时奇异性问题。平滑解和斑块解在 $\alpha\in(1,2)$ 的局部适定性已有相当成熟的理论，但在全平面、无外力的情形下，平滑初始界面是否会在有限时间形成奇异性仍是开放问题。已有数值工作提出了相似塌缩、细丝不稳定性及曲率爆破等情景；另一方面，关于 splash 型自接触的若干结果排除了界面保持正则而仅发生触碰的情形。这使得局部、自相似的边界奇异性成为值得检验的机制。本文的目标是给出严格构造：从随时间演化且在每个 $$t<T$$ 都光滑的单个无界斑块界面出发，在有限时刻 $$T$$ 形成角点，而不是仅以数值图像或形式相似解描述可能的爆破。

前沿方程具有平移与缩放不变性；特别是

F(x,t)\mapsto c^{-1}F(cx,c^\alpha t),\qquad c>0,

而一次齐次角点 $x_-:=\max\{-x,0\}$ 与该尺度相容。论文所构造的角点斜率为 $\sqrt\epsilon$ ，其中 $\epsilon$ 是小参数，故其本质是平坦界面的小斜率扰动。构造所采用的相似尺度为 $L(t)=(\alpha(T-t))^{1/\alpha}$ ，对应界面幅度与空间尺度同阶。本文试图弥补的空白由此很明确：在 $\alpha\in(1,2)$ 的强奇异速度范围内，严格构造一个全平面、单斑块、无外力的有限时角点形成解，并给出其剖面的远场结构和曲率爆破率。需注意，结论是关于特定的自相似斑块解，并不证明一般光滑初值的有限时奇异性，也不解决所有 $\alpha$ 的问题。

2. 主要定理或主要结果

主要结论是定理 1.1。若 $\alpha\in(7/4,2)$ ，则存在 $\epsilon_0>0$ ，使得对每个 $0<\epsilon<\epsilon_0$ 及每个 $$T>0$$ ，前沿方程存在经典解 $F:\mathbb R\times(-\infty,T)\to\mathbb R$ ，且具有自相似形式

F(x,t)=\sqrt\epsilon\,L(t)f_\epsilon\!\left(\frac{x}{L(t)}\right),\qquad L(t)=\big(\alpha(T-t)\big)^{1/\alpha}.

剖面 $f_\epsilon\in C^\infty(\mathbb R)$ 。在负半轴，其远场接近斜率为 $$-1$$ 的直线，且对 $0\le j\le2$ 有

|\partial^j(f_\epsilon(y)+y)|\lesssim\langle y\rangle^{-\alpha+1-j},\qquad y\le0.

在正半轴，剖面趋于零，并满足

|f_\epsilon(y)|\lesssim\langle y\rangle^{-\alpha+1},\qquad |\partial^{j+1}f_\epsilon(y)|\lesssim\langle y\rangle^{-\frac{\alpha}{2(\alpha-1)}+\frac{j}{\alpha-1}},\quad 0\le j\le1, \qquad y\ge0.

这里 $\langle y\rangle$ 表示通常的权重 $\sqrt{1+y^2}$ ；隐含常数的具体依赖关系在所给摘录中未完整列出。界面趋于角点的定量估计是

\|F(\cdot,t)-\sqrt\epsilon\,x_-\|_{L^\infty(\mathbb R)}\lesssim\sqrt\epsilon\,(T-t)^{1/\alpha}.

此外，曲率在角点附近至少以 $(T-t)^{-1/\alpha}$ 的速率发散：对每个 $\delta>0$ ，存在 $c=c(\epsilon)>0$ 和 $t_\delta<T$ ，使得

\|\partial_x^2F(\cdot,t)\|_{L^\infty(-\delta,\delta)}\ge c(T-t)^{-1/\alpha},\qquad t_\delta<t<T.

因此，严格证明的是 $$t<T$$ 时的光滑自相似演化、 $L^\infty$ 意义下趋向小斜率角点，以及局部曲率的下界爆破；摘录并未给出曲率的匹配上界，故不能据此声称曲率恰好等于该阶。定理还允许任意指定 $$T>0$$ 。由方程的反射时间对称性 $F(x,t)\mapsto F(-x,-t)$ ，作者指出该结果也可读作角点的瞬时平滑：在反向时间表述中，初始角点 $\sqrt\epsilon x_+$ 对应于之后的光滑界面。这是定理的对称推论，不是另一套独立构造。

参数限制 $\alpha>7/4$ 与所选加权空间中的指数条件有关。作者要求存在 $\gamma$ 满足

\frac{\alpha+2}{2(\alpha-1)}<\gamma\le2\alpha-1,

此区间非空恰对应 $\alpha>7/4$ 。论文明确不主张 $$7/4$$ 是实际最优阈值；该限制是当前近似解修正及函数空间闭合所产生的技术条件。

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

证明采用两阶段扰动构造：先造残差足够小的近似自相似剖面，再对其线性化并解修正方程。相似变量取 $\tau=\log(1/L(t)$ 、 $$y=x/L(t)$$ ，并写 $$F=L f$$ ；再把剖面按小斜率因子归一化后，定常剖面满足文中方程 (2.2)。其核心线性算子为

L_0(f)=|D|^{\alpha-1}f'+f-yf'.

困难在于原前沿算子是拟线性的、导数阶数较高，直接围绕平坦界面做固定点会面对不利的导数映射性质。作者利用非线性核的结构，将算子分解为

J_\epsilon(f)=A_\epsilon(f')|D|^{\alpha-1}f'+\widetilde J_\epsilon(f),\qquad A_\epsilon(q)=(1+\epsilon q^2)^{-\alpha/2},

其中 $\widetilde J_\epsilon$ 是低一阶项：在 $f\in W^{2,\infty}_{\rm loc}$ 时，核中的差分抵消使其局部奇性比原主项弱一阶。对定常方程除以 $A_\epsilon(f')$ ，遂得到以 $$L_0$$ 为主部、非线性项导数阶数较低的半线性形式。这一改写是固定点分析能够闭合的关键，而不是简单地把高阶项称为小扰动。

近似解从 $$L_0(f_0)=0$$ 的显式傅里叶表示出发：

f_0(y)=-\frac y2+\frac1\pi\int_0^\infty\frac{1-\cos(y\xi-\Phi_0(\xi))}{\xi^2}\,d\xi, \qquad \Phi_0(\xi)=\frac{\xi|\xi|^{\alpha-1}}{\alpha}.

其负端呈慢变尾部，正端则由驻相产生 Airy 型快速振荡，振幅约为 $y^{-\alpha/(2(\alpha-1))}$ ，相位约为 $\Psi_0(y)=\frac{\alpha-1}{\alpha}y^{\alpha/(\alpha-1)}$ 。近似剖面分成线性种子、慢变修正和振荡修正。慢变部分以 $O(\epsilon)$ 的 $|y|^{-\alpha+1}$ 尾部抵消非线性算子作用于增长尾部产生的慢衰减响应。振荡部分则不能只保留线性相位：非线性作用会产生与基频共振的三次项，以及增长尾部与正向振荡尾部的耦合项。作者分别加入三次谐波校正和相位修正；其中对数相位项抵消共振三次响应，幂次相位项抵消尾部耦合。再引入复参数 $\mu$ 调整正频振荡分量，以便控制线性逆算子无法直接容纳的两个实方向的残差。

修正方程在非对称加权空间中处理。 $X_\gamma$ 控制修正量及其两阶导数， $Y_\gamma$ 控制残差及其一阶导数；正负无穷端采用不同导数权重，以分别适配正端快速振荡和负端慢变行为。一个核心障碍是 $L_0^{-1}$ 在完整 $Y_\gamma$ 上不具备所需衰减。其傅里叶表示含有

\widehat{L_0^{-1}G}(\xi)=\frac{e^{-i\Phi_0(\xi)}}{\xi^2}\int_0^\xi\eta e^{i\Phi_0(\eta)}\widehat G(\eta)\,d\eta.

高频端的积分常数一般不消失，会留下仅按 $y^{-(\alpha+2)/(2(\alpha-1))}$ 衰减的驻相尾部；临界权重甚至可能多出对数因子。因此作者把允许的强迫限制到满足共振消去条件的子空间

\widetilde Y_\gamma=\{G\in Y_\gamma:\ell(G)=0\},\qquad \ell(G)=\int_0^\infty\eta e^{i\Phi_0(\eta)}\widehat G(\eta)\,d\eta,

并在该子空间建立逆估计。为使非线性残差也满足此相容条件，使用 Lyapunov–Schmidt 约化：把强迫投影为共维二的部分与其补空间，先在投影后的子空间用压缩映射求修正，再利用复参数 $\mu$ 解有限维的两个实相容方程。因而 $\mu$ 不只是记号上的自由度，而是消除逆算子共振障碍所需的参数。

线性逆估计也不是单一的粗糙核界。作者将其分解为低频、高频驻相和高频非驻相三部分。低频处符号在零点不光滑，需单独作加权 $L^\infty$ 乘子估计；高频驻相项显式保留振荡核及其权重；远离驻点的高频项则视为具有光滑符号的伪微分算子并估计。非线性估计中，局部项由加权导数控制，非局部项通过近远场分割，并利用核衰减和差分抵消控制。另一个不能忽略的结构是增长尾部与远端修正之间的非局部耦合：尽管其系数为 $O(\epsilon)$ ，该项在给定权重下仍不一定映到 $Y_\gamma$ ，故被并入线性算子；通过可逆共轭变换 $L_\epsilon=V_\epsilon L_0U_\epsilon$ 将其化归自由算子。摘录将近似解、线性逆和非线性估计分别组织在定理 3.1、3.2、3.4 及相应章节中；未提供这些定理的全部技术假设与常数细节。

4. 与之前工作的比较

相对于已有局部理论，本文的贡献不是重新证明 $\alpha\in(1,2)$ 的局部适定性。作者在引言中指出，该范围内光滑解的 Sobolev 局部适定性已有 Chae、Constantin、Córdoba、Gancedo、Wu 等人的结果，临界指数附近亦有后续发展；斑块局部存在理论也已有多项工作。本文处理的是不同问题：构造特定自相似轨道，使一个单斑块图形界面在有限时间形成角点，并给出其剖面尾部及曲率爆破下界。结果不能直接推出一般初值的爆破，也不与常规局部适定性矛盾，因为每个 $$t<T$$ 的解仍光滑，失去光滑性发生在极限时刻。

相对于小扰动全局结果，Córdoba、Gómez-Serrano、Ionescu 的工作表明，在 $\alpha\in(1,2)$ 中，平坦界面的充分小、光滑且衰减扰动可全局演化并发生色散及修正散射。本文的剖面在小斜率意义下也接近平坦界面，但其远场不是通常的快速衰减扰动：一侧趋近线性斜坡，另一侧带有特定幂次尾部和高频振荡。因此它并不落在上述小而衰减数据结论直接排除的情形中。本文揭示了色散并非只能导致扰动散开；在特定的非衰减、非对称自相似背景及精细相位调制下，色散结构也参与奇异剖面的构造。

相对于此前的奇异性证据，Córdoba、Fontelos、Mancho、Rodrigo 关于双斑块的数值研究观察到间距塌缩与边界曲率增长，Scott、Dritschel 等工作报告了 SQG 斑块细丝不稳定性和尺度不变奇异情景。这些是数值机制或候选情景；本文给出全平面前沿方程中单个无界斑块的严格自相似构造，但参数范围是 $\alpha\in(7/4,2)$ ，几何类型是图形角点形成，不能视为对所有数值情景的证明。与有边界问题中的 Kiselev、Ryzhik、Yao、Zlatoš 等人关于小 $\alpha$ 斑块有限时奇异性，以及 Zlatoš 关于半平面光滑数据爆破的结果相比，本文不依赖物理边界；但它研究的是斑块边界角点，且只覆盖较靠近 $$2$$ 的 $\alpha$ 范围。Jeon、Zlatoš 对全平面相互接触斑块的结果涉及多个斑块接触机制，也与本文单一图形界面的自相似角点不同。

与现有自相似解及奇异边界研究相比，Elgindi、Jeong 关于二维 Euler 角点涡斑的理论、García、Gómez-Serrano 关于 $\alpha\in(0,1)$ 的自相似螺旋解，以及相关齐次稳态研究，提供了相邻背景；本文则把严格自相似构造放到强奇异速度的 $\alpha\in(1,2)$ 区间，并具体处理前沿方程中的振荡尾部和非局部耦合。文中还提及强迫 gSQY 方程的有限时爆破构造，但那些结果允许外力；本文的构造是无外力的，却是特殊斑块自相似解，并非一般有限能量光滑数据的爆破结论。本文解的无界图形、非衰减斜坡尾部以及斑块结构均是适用范围的一部分；摘录未给出有限能量性质，不能将其归入有限能量解的结论。

最重要的限制是阈值 $\alpha>7/4$ 。作者明确说明该阈值来自权重指数必须满足 $\frac{\alpha+2}{2(\alpha-1)}<\gamma\le2\alpha-1$ ，而非已证明的方程内在临界值。他们指出，增加近似解校正项可能把上界推进至 $\gamma\le\alpha+1$ ，从而有望覆盖 $\alpha>(1+\sqrt{33})/4\approx1.686$ ；继续降低阈值还可能需要强迫满足额外矩条件（例如 $\int G=0$ ）并进一步改善残差衰减。这些是作者讨论的可能改进，不是本文已证明的结果。因此本文的推进在严格性和机制解释上显著，但代价是较窄的参数区间、专门设计的渐近剖面、有限维相容性调节及强烈依赖权重的技术框架。

3. 可压缩 Navier-Stokes-Poisson 系统在 Navier 滑移边界条件下的一致正则性

Uniform Regularity for the Compressible Navier-Stokes-Poisson System under Navier-Slip Boundary Conditions

作者: Hairong Liu, Jinlong Wu

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32454v1

摘要翻译: 本文证明了在三维有界光滑区域内，满足一般 Navier 滑移边界条件的可压缩 Navier-Stokes-Poisson 系统存在局部时间强解。我们证明该解在共法 Sobolev 空间中一致有界，并得到与黏性系数无关的一致正则性估计。作为推论，通过强紧性论证，我们得到相应可压缩 Euler-Poisson 系统的无黏极限。

Comment: 无

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

1. 研究对象和背景

本文研究三维有界光滑区域 $\Omega\subset\mathbb{R}^3$ 内的可压缩 Navier–Stokes–Poisson 系统。未知量为密度 $\rho^\varepsilon$ 、速度 $u^\varepsilon$ 和静电势 $\varphi^\varepsilon$ ； $\varepsilon\in(0,1]$ 是黏性参数，压力满足 $p^\varepsilon=(\rho^\varepsilon)^\gamma$ ，其中 $\gamma>1$ 。具体方程为

\begin{cases} \rho_t^\varepsilon+\operatorname{div}(\rho^\varepsilon u^\varepsilon)=0,\\ \rho^\varepsilon u_t^\varepsilon+\rho^\varepsilon u^\varepsilon\cdot\nabla u^\varepsilon+\nabla p^\varepsilon =\mu\varepsilon\Delta u^\varepsilon+(\lambda+\mu)\varepsilon\nabla\operatorname{div}u^\varepsilon+\rho^\varepsilon\nabla\varphi^\varepsilon,\\ \Delta\varphi^\varepsilon=\rho^\varepsilon-1, \end{cases} \qquad x\in\Omega,

其中 $\mu>0$ 、 $2\mu+3\lambda>0$ ，Poisson 方程中的常数背景电荷密度归一化为 $$1$$ 。速度满足一般 Navier–滑移边界条件，电势采用齐次 Neumann 条件：

u^\varepsilon\cdot n=0,\qquad (S u^\varepsilon\cdot n)_\tau=-(A u^\varepsilon)_\tau,\qquad \partial_n\varphi^\varepsilon=0\quad\text{于 }\partial\Omega,

这里 $S u=\frac12(\nabla u+(\nabla u)^T)$ ， $$n$$ 为外法向， $$A$$ 是边界上的光滑对称矩阵， $\tau$ 表示切向分量。对光滑解，速度条件也可写为涡量形式 $u^\varepsilon\cdot n=0$ 、 $n\times\omega^\varepsilon=[Bu^\varepsilon]_\tau$ ，其中 $\omega^\varepsilon=\nabla\times u^\varepsilon$ 、 $$B=2(A-S(n))$$ 。初始电势由 $\Delta\varphi_0=\rho_0-1$ 、 $\partial_n\varphi_0=0$ 确定，因此还需满足 Neumann 可解性条件 $\int_\Omega(\rho_0-1)\,dx=0$ 。

该模型描述带电可压缩流体与其自洽静电场的耦合，在等离子体、半导体及电离气体等模型中出现。相较于通常的可压缩 Navier–Stokes 系统，Poisson 方程使电场由密度非局部决定，而力项 $\rho\nabla\varphi$ 又反馈进入动量方程；高阶能量估计因而必须同时处理流体未知量、电势的椭圆正则性及边界条件。研究的问题是：在黏性趋于零时，能否在与 $\varepsilon$ 无关的时间区间上控制解的正则性，并据此严格导出对应的可压缩 Euler–Poisson 极限。

形式上令 $\varepsilon\to0$ ，黏性项消失，极限系统为

\begin{cases} \rho_t+\operatorname{div}(\rho u)=0,\\ \rho u_t+\rho u\cdot\nabla u+\nabla p=\rho\nabla\varphi,\\ \Delta\varphi=\rho-1,\qquad p=\rho^\gamma, \end{cases} \qquad u\cdot n=0,\quad\partial_n\varphi=0\quad\text{于 }\partial\Omega.

有边界时，黏性与无黏方程的边界机制不同，常规各向同性 Sobolev 范数一般不能直接给出适用于零黏性极限的一致估计。Navier–滑移条件与 Euler 的不可穿透条件相容，因而比无滑移边界更适于研究该极限；但曲边界仍会造成法向正则性控制的困难。论文沿用 conormal Sobolev 框架，以切于边界的向量场导数刻画可一致控制的方向，再利用方程和边界涡量结构恢复法向控制。本文试图填补的具体空白，是把可压缩 Navier–Stokes 系统的滑移边界一致正则性方法扩展到含非局部 Poisson 耦合的三维有界光滑区域，并在无阻尼、齐次 Neumann 电势边界条件下证明局部一致正则性及强收敛的 Euler–Poisson 极限。齐次 Neumann 条件使电势能量与分部积分相容，但并不消除高阶估计中电场与密度、速度的耦合困难。

2. 主要定理或主要结果

主要结论是定理 1.1 的一致正则性和定理 1.2 的无黏极限。设 $m\ge6$ ， $\Omega$ 为 $C^{m+2}$ 区域，边界矩阵 $A\in C^{m+1}(\partial\Omega)$ 。初值满足密度正下界与上界

0<\frac{1}{\widetilde C_0}\le \rho_0(x)\le \widetilde C_0<\infty,

相容条件及对所有 $0<\varepsilon\le1$ 一致的初始范数界 $\sup_{0<\varepsilon\le1}\|(p_0,u_0,\nabla\varphi_0)\|_{X_m^\varepsilon}\le\widetilde C_0'$ 。其中 $p_0=\rho_0^\gamma$ ； $X_m^\varepsilon$ 是论文定义的 conormal 能量范数，包含 $\|(p,u,\nabla\varphi)\|_{H^m}$ 、 $\|\nabla u\|_{H^{m-1}}$ 、压力的时间导数与 $\Delta p$ 控制、 $\|(\nabla u,\nabla\varphi)\|_{H^{1,\infty}}$ ，以及带权高阶项 $\varepsilon\|\nabla\partial_t^{m-1}p\|_{L^2}^2$ 和 $\varepsilon\|\Delta p\|_{H^2}^2$ 。文中 $$H^m$$ 表示包含时间导数 $\partial_t$ 与空间 conormal 导数的时空 conormal 范数，纯空间 conormal 范数另行记为 $\|\cdot\|_m$ 。

定理 1.1 断言存在 $$T_0>0$$ 和常数 $\widetilde C_1>0$ ，仅依赖于初值界、 $$m$$ 、区域和边界矩阵，而与 $\varepsilon$ 无关，使得每个 $0<\varepsilon\le1$ 的初边值问题在 $$[0,T_0]$$ 上存在唯一强解，且

\frac{1}{2\widetilde C_0}\le\rho^\varepsilon(t,x)\le2\widetilde C_0.

更具体地，解满足下列对 $\varepsilon$ 一致的完整估计：

\begin{aligned} \sup_{0\le t\le T_0}\Big(&\|(p^\varepsilon,u^\varepsilon,\nabla\varphi^\varepsilon)(t)\|_{H^m}^2 +\|\nabla u^\varepsilon(t)\|_{H^{m-1}}^2 +\sum_{k=0}^{m-2}\|\partial_t^k\nabla p^\varepsilon(t)\|_{m-1-k}^2 +\|\Delta p^\varepsilon(t)\|_{H^1}^2\\ &+\| (\nabla u^\varepsilon,\nabla\varphi^\varepsilon)(t)\|_{H^{1,\infty}}^2 +\varepsilon\|\nabla\partial_t^{m-1}p^\varepsilon(t)\|_{L^2}^2 +\varepsilon\|\Delta p^\varepsilon(t)\|_{H^2}^2\Big)\\ &+\int_0^{T_0}\|\nabla\partial_t^{m-1}p^\varepsilon\|_{L^2}^2\,dt +\int_0^{T_0}\|\Delta p^\varepsilon\|_{H^2}^2\,dt +\varepsilon\int_0^{T_0}\|\nabla u^\varepsilon\|_{H^m}^2\,dt\\ &+\varepsilon\sum_{k=0}^{m-2}\int_0^{T_0}\|\nabla^2\partial_t^k u^\varepsilon\|_{m-1-k}^2\,dt +\varepsilon^2\int_0^{T_0}\|\nabla^2\partial_t^{m-1}u^\varepsilon\|_{L^2}^2\,dt \le\widetilde C_1. \end{aligned}

定理 1.2 在相同假设下给出同一时间区间 $$[0,T_0]$$ 上的 Euler–Poisson 强解，满足 $\sup_{0\le t\le T_0}\|(p,u,\nabla\varphi)(t)\|_{X_m^0}<\infty$ 。以相同初值为初始数据，黏性解趋于该极限解，且收敛在时间上一致、在空间的 $$L^2$$ 与 $L^\infty$ 范数下成立：

\lim_{\varepsilon\to0}\sup_{0\le t\le T_0} \left( \|(p^\varepsilon,u^\varepsilon,\nabla\varphi^\varepsilon)-(p,u,\nabla\varphi)\|_{L^2(\Omega)} +\|(p^\varepsilon,u^\varepsilon,\nabla\varphi^\varepsilon)-(p,u,\nabla\varphi)\|_{L^\infty(\Omega)} \right)=0.

定理没有给出关于 $\varepsilon$ 的定量收敛速率；其结论是通过一致估计和强紧性得到的强收敛。文中给出的结果为局部时间结论，未声称无阻尼系统的全局存在或长期衰减。

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

证明的核心难点是，在边界附近 conormal 导数不控制一般法向导数，而动量方程中的 Poisson 力又要求同时控制密度产生的电场。作者采用分层先验估计，而不是直接对完整系统作一个单一的高阶能量估计。以下技术细节依据论文引言及先验估计部分的证明路线说明；所给 PDF 摘录只展示了基本能量估计的开头，后续各引理的逐项计算并未完整提供。

首先，建立适用于边界的 conormal 框架。局部边界坐标下使用切向向量场及退化法向场 $Z_3=\phi(z)\partial_z$ ，其中 $\phi(0)=0$ ，并令 $Z_0=\partial_t$ 。这使 conormal 向量场在边界处不产生不可控的法向边界项，同时允许以时空 conormal 范数组织高阶估计。证明依赖 conormal 乘积估计、各向异性 Sobolev 嵌入和迹估计；这些工具把高阶 $$L^2$$ 控制与低阶 $L^\infty$ 控制结合起来，处理非线性乘积及边界迹。对向量场，还使用 div–curl 型椭圆估计，以涡量、散度、低阶范数和边界迹恢复 conormal 正则性。

其次，Poisson 耦合必须纳入能量而不能造成导数损失。对 $\Delta\varphi=\rho-1$ 和齐次 Neumann 条件，作者使用椭圆估计控制电场的 conormal 正则性，例如将 $\nabla^2\varphi$ 的 conormal 范数归结为密度源项及低阶电场范数。低阶基本能量中，动量方程与连续性方程、Poisson 方程配合，使静电能 $\frac12\|\nabla\varphi\|_{L^2}^2$ 与流体能量共同出现。摘录给出的基本估计形式为

\sup_{0\le\tau\le t}\left\{ \int_\Omega\left(\frac12\rho|u|^2+\frac{1}{\gamma-1}\rho^\gamma\right)d x +\frac12\|\nabla\varphi(\tau)\|_{L^2}^2\right\} +c\varepsilon\int_0^t\|\nabla u\|_{L^2}^2d\tau \le E(0)+C\int_0^t\|u\|_{L^2}^2d\tau.

这里 Neumann 条件允许在椭圆方程相关项中分部积分，而电场能被直接并入能量；高阶层面则需估计 conormal 微分作用于 $\rho\nabla\varphi$ 后产生的交换子和乘积项。作者指出，这些项由 conormal 乘积估计与 Neumann Poisson 椭圆估计控制，额外的 $\|\nabla\varphi\|_{H^m}^2$ 正是相较于纯可压缩 Navier–Stokes 情形新增的能量组成。

再次，压缩性使 $\operatorname{div}u$ 不为零，因此不能只靠涡量估计恢复速度的法向正则性。证明路线还要控制 $\operatorname{div}u$ 、压力梯度和纯时间导数，尤其借助压力方程获得不带黏性退化的纯时间导数估计。随后通过速度的涡量方程以及 Navier–滑移条件的涡量形式处理法向速度导数：引入满足齐次 Dirichlet 边界条件的边界涡量量，估计其中来自电场的源项，并使用先前得到的 conormal 流体估计和 Poisson 椭圆估计闭合。这一步利用的是边界条件的结构，而非把法向导数简单地视为可由 conormal 范数控制。

最后，作者通过各向异性嵌入、边界涡量估计和电势椭圆估计导出 $L^\infty$ 型控制，并单独建立 $\Delta p$ 的一致估计以控制压力的 Lipschitz 正则性。上述估计合并后，以近似构造和连续性论证得到定理 1.1；再由一致界、紧性及极限方程识别得到定理 1.2。摘录并未给出紧性论证使用的具体紧嵌入定理、收敛子列的全部拓扑细节或逐项误差估计，因此这些细节从提供内容中无法确认。特别地，本文证明的关键价值并非仅是使用能量估计或紧性，而是把电场纳入 conormal 高阶能量，并将压力散度控制、边界涡量法向估计及 Poisson 椭圆正则性依次衔接，使常数和存在时间不依赖于黏性参数。

4. 与之前工作的比较

论文把自身定位在滑移边界下零黏性极限的 conormal 正则性框架内。无边界情形中，当 Euler 解保持光滑时，黏性趋零通常可在标准 Sobolev 框架下处理；物理边界会改变问题结构。无滑移条件与无黏 Euler 的不可穿透条件不相容，通常伴随边界层，相关极限除解析初值、初始涡量远离边界或特定稳态等特殊情形外仍很困难。Navier–滑移条件与 Euler 边界条件相容，因而已有工作能够在有界区域建立一致正则性。

在不可压缩流方面，Masmoudi 与 Rousset [22] 在一般区域发展 conormal Sobolev 估计，并通过紧性论证零黏性极限。可压缩情形方面，Wang、Xin 与 Yong [34] 将该框架用于一般 Navier–滑移边界下的可压缩 Navier–Stokes 方程，Wang [33] 又将相关分析用于全可压缩 Navier–Stokes 系统。本文明确以这些工作为主要分析背景；相较于它们，新增困难不是单纯更换边界条件，而是要控制 Poisson 方程所定义的非局部电场，并估计动量方程中 $\rho\nabla\varphi$ 对高阶 conormal 能量、法向估计和 $L^\infty$ 控制的影响。本文为能量空间增加 $\|\nabla\varphi\|_{H^m}$ ，并以 Neumann 椭圆估计处理该耦合。

与 Navier–Stokes–Poisson 的其他奇异极限研究相比，论文介绍了全空间或周期区域中向不可压缩 Euler、Euler–Poisson 等系统的极限工作 [31, 6, 11, 18]，以及半空间问题 [12]。边界电势条件会显著改变边界层机制：Ju 与 Xu [12] 研究半空间中的准中性及零黏性极限；Ju、Luo 与 Xu [13] 处理速度 Navier–滑移条件与电势 Dirichlet 条件下的黏性等奇异极限，论文指出其主要困难包括由指定边界电势产生的强密度边界层，以及它与较弱速度边界层的相互作用。本文采用齐次 Neumann 电势条件，作者强调这种设定下电静能量的分部积分结构更自然，且不存在由指定边界电势引发的同类强边界层机制。不过这并不意味着边界与 Poisson 耦合完全无碍：高阶电场仍需与密度和速度一致控制。

与 Gao、Peng、Wei、Yao [7] 的工作相比，本文处理的是不含阻尼项的系统，并考虑一般三维有界光滑区域、一般 Navier–滑移边界以及齐次 Neumann 电势条件。[7] 研究带阻尼项的滑移边界模型，获得全局时间一致界和带收敛速率的无黏极限；论文指出阻尼的稳定作用支持其全局估计。本文没有阻尼，因此结论限于与黏性无关的局部时间 $$[0,T_0]$$ 。本文的定理只证明 $$L^2$$ 与 $L^\infty$ 意义下的强收敛，并未给出关于 $\varepsilon$ 的收敛速率。这是可直接辨认的结论差异：本文覆盖无阻尼的局部适定性与一般有界曲域，但没有声称带阻尼工作所获得的全局性和定量收敛率。

假设与结论的边界也需准确区分。本文要求初始密度有正的上下界、Poisson Neumann 相容条件成立，并要求初始数据在所定义的高阶 conormal 范数中对 $\varepsilon$ 一致有界；结论依赖 $m\ge6$ 、 $\Omega\in C^{m+2}$ 和 $A\in C^{m+1}(\partial\Omega)$ 。它没有据此证明标准各向同性 Sobolev 范数下所有法向导数均一致有界，也没有处理无滑移条件下的极限问题。作者还在备注中说明，可将 conormal 框架调整到电势 Dirichlet 条件 $\varphi|_{\partial\Omega}=0$ 的绝缘边界设置，但此时 $\partial_n\varphi$ 通常非零，低阶散度估计会出现边界项 $\int_{\partial\Omega}\rho\,\partial_n\varphi\,\operatorname{div}u\,dS$ 。作者建议先取得高阶纯时间导数控制，再回头控制该迹项；这属于备注所述的可调整分析路线，并非本文定理 1.1–1.2 已陈述的边界条件结论。总体而言，本文的推进是在现有可压缩滑移边界 conormal 理论上纳入 Poisson 耦合，代价是高阶初值与边界光滑性要求以及仅有局部时间、无定量收敛率；更细的定量比较或其他文献假设的完整对照，作者未在所提供摘录中显式展开。

4. 广义相对论中具有自由边界的正压液态恒星的线性稳定性与转折点原理

Linear Stability and a Turning Point Principle for Barotropic Liquid Stars with a Free Boundary in General Relativity

作者: Zeming Hao, Jianchao Xu

类别: math.AP, gr-qc

发布日期: 2026-10-23

arXiv ID: 2610.32138v1

摘要翻译: 我们研究描述被真空包围的液态恒星的 Einstein-Euler 系统中球对称、渐近平坦稳态的线性稳定性。状态方程具有正的表面密度、正的因果声速以及渐近线性的高密度极限。我们构造了平衡态族，并推导出带 Robin 边界条件的线性化方程。线性化算子具有离散实谱，其负 Morse 指数等于增长模态的数目。我们证明，对于中心红移较大的稳态，线性化算子存在增长模态；而对于中心红移较小的稳态，不存在此类增长模态。平衡态族关于参数的导数给出线性化方程的一个正则内部零解，其边界缺陷是总质量关于参数的导数的非零倍数。由此得到转折点原理：当且仅当质量达到临界点时，核是一维的；穿越特征值的符号由质量—半径曲线的弯曲方向决定。

Comment: 29 pages

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

1. 研究对象和背景

本文研究广义相对论中的球对称液态星: 流体占据一个随时间运动的有界区域, 外部是真空, 而流体侧的表面密度为正, 因而密度在自由界面处发生跃变. 研究对象是渐近平坦时空中的 Einstein-Euler 系统. 理想流体的能量动量张量为

T_{\alpha\beta}=(\rho+p)u_\alpha u_\beta+p g_{\alpha\beta},

并满足

\bar R_{\alpha\beta}-\frac12g_{\alpha\beta}\bar R=8\pi T_{\alpha\beta},\qquad \nabla_\alpha T^{\alpha\beta}=0.

这里 $$g$$ 是时空度量, $$u$$ 是单位类时四速度, $\rho$ 和 $$p$$ 分别是质量能量密度与压强. 状态方程取 $p=P(\rho)$ , 流体区域记为 $$B(t)$$ , 真空边界满足压强为零且流体速度与边界相切:

p=P(\rho),\qquad p|_{\partial B(t)}=0,\qquad u|_{\partial B(t)}\in T(t,\partial B(t)).

文章采用的液体型状态方程满足 $\rho_0>0$ , $P\in C^2([\rho_0,\infty))$ , $P(\rho_0)=0$ , $0<P'(\rho)\leq 1$ . 因此表面流体侧有 $\rho=\rho_0$ , 而真空侧 $\rho=p=0$ . 另外, 存在 $c_*^2\in(0,1]$ 使得高密度时状态方程渐近线性: $|P(\rho)-c_*^2\rho|$ 有界, $P'(\rho)$ 以幂次速度趋于 $$c_*^2$$ , 且 $P''(\rho)\to0$ . 这些条件兼顾正的表面密度、正且不超光速的声速以及高密度渐近行为. 硬相流体 $P(\rho)=\rho-\rho_0$ 是其中的典型模型, 但本文的论证覆盖更一般的压强依赖型状态方程.

在球对称 Schwarzschild 坐标下, 度量写为

g=-e^{2\mu(t,r)}dt^2+e^{2\lambda(t,r)}dr^2+r^2d\omega^2,

其中 $d\omega^2$ 是单位球面度量. 静态平衡时流体速度为零, Einstein-Euler 方程化为 Tolman-Oppenheimer-Volkoff 型系统:

e^{-2\lambda}(2r\lambda'-1)+1=8\pi r^2\rho,

e^{-2\lambda}(2r\mu'+1)-1=8\pi r^2p, \qquad (\rho+p)\mu'+p'=0.

质量函数 $$m$$ 满足

e^{-2\lambda}=1-\frac{2m(r)}r,\qquad m(r)=4\pi\int_0^r s^2\rho(s)\,ds.

静态星体在首个压强为零的位置 $R_\kappa$ 截止, 总质量为 $M_\kappa=m_\kappa(R_\kappa)$ ; 外部由真空 Schwarzschild 区域描述. 通过相对论焓

Q(\rho)=\int_{\rho_0}^{\rho}\frac{P'(s)}{s+P(s)}\,ds

可将静态 Euler 方程积分, 并以中心红移 $\kappa=Q(\rho_\kappa(0))$ 标记一参数平衡态族.

核心问题是这些平衡态在径向扰动下是否具有指数增长模态, 以及谱稳定性如何随平衡态参数变化. 对随流体运动的径向坐标作扰动参数化 $\eta(y,t)=y(1+\zeta(y,t))$ 后, 线性化系统具有带自由表面 Robin 边界条件的形式

y e^{2\lambda_\kappa-2\mu_\kappa}\ddot\zeta+L_\kappa\zeta=0.

因此问题转化为带权广义特征值问题. 其困难不只在于引力与流体的耦合, 还在于液体表面的密度跳跃会给线性化边界贡献表面项, 改变算子的边界条件与谱结构. 本文试图在这一自由边界液体模型中建立可用的线性谱理论, 证明小中心红移与大中心红移区域的稳定性差异, 并把径向不稳定模态数的变化与总质量曲线的转折联系起来. 这里的稳定性结论是线性谱意义下的结论, 不能直接等同于一般扰动下的非线性长期稳定性.

2. 主要定理或主要结果

平衡态构造. 在文中给定的状态方程假设下, 对每个中心红移 $\kappa>0$ 构造唯一球对称静态星体. 星体半径 $R_\kappa$ 有限, 流体侧表面密度为 $\rho_0$ , 真空侧密度为零, 总质量为 $M_\kappa=m_\kappa(R_\kappa)$ . 高密度分析中作者也使用中心压强参数化, 并明确指出它与中心红移参数化是单调对应的重新参数化, 不能把两种记号下的参数定义混同.
线性化谱与稳定性. 对每个平衡态, 线性化径向位移满足

y e^{2\lambda_\kappa-2\mu_\kappa}\ddot\zeta+L_\kappa\zeta=0,

并满足自由表面 Robin 条件. 相应线性化算子具有离散实谱; 其负 Morse 指数等于指数增长模态的数目. 定理 1.1 证明: 当中心红移 $\kappa$ 充分小时, 不存在指数增长模态; 当 $\kappa$ 充分大时, 至少存在一个指数增长模态. 这里的小参数区间和大参数阈值依赖于所给状态方程, 摘录未给出显式阈值或增长率. “小红移稳定”指无增长模态的线性谱稳定性, 不是非线性稳定性结论.

转折点处的核刻画. 令 $$y$$ 为平衡态的径向变量, 在固定 $$y$$ 下对平衡态族求参数导数, 定义边际位移

v_\kappa(y)=-\frac{\partial_\kappa m_\kappa(y)}{4\pi y^3(\rho_\kappa(y)+p_\kappa(y))}.

定理 1.2(i) 给出其弱形式边界缺陷恒等式: 对任意 $\chi\in H^1$ ,

\langle L_\kappa v_\kappa,\chi\rangle_{\Lambda_\kappa} =-b_\kappa\,\partial_\kappa M_\kappa\,\chi(R_\kappa), \qquad b_\kappa=\frac{\Lambda_\kappa(R_\kappa)M_\kappa}{4\pi\rho_0R_\kappa^3(R_\kappa-2M_\kappa)}>0.

因此

\ker L_\kappa\ne\{0\}\quad\Longleftrightarrow\quad \partial_\kappa M_\kappa=0.

在质量临界点, 核恰为一维且由 $v_\kappa$ 张成; 非临界点核为零. 该结论中的 $\langle\cdot,\cdot\rangle_{\Lambda_\kappa}$ 是文中线性化算子的带权双线性形式.

不稳定模态数的转折规律. 若 $\kappa_0$ 是非退化质量临界参数, 即

\partial_\kappa M_\kappa\big|_{\kappa_0}=0, \qquad \partial_\kappa^2 M_\kappa\big|_{\kappa_0}\ne0,

则 $\partial_\kappa R_\kappa|_{\kappa_0}\ne0$ . 不稳定模态数 $n_u(\kappa)$ 只能在质量函数 $M_\kappa$ 的极值处改变, 且每经过一个这样的非退化临界点只改变一个模态. 令

\mathcal P(\kappa)=(\partial_\kappa M_\kappa)(\partial_\kappa R_\kappa).

跨越临界点时, 若 $\mathcal P$ 的符号由负变正, $$n_u$$ 增加 $$1$$ ; 若由正变负, $$n_u$$ 减少 $$1$$ . 等价地, 穿越零点的特征值符号满足

\operatorname{sgn}\nu_*(\kappa)=-\operatorname{sgn}\big[(\partial_\kappa M_\kappa)(\partial_\kappa R_\kappa)\big].

大中心红移的定量渐近信息. 在高密度极限, 对合适的缩放环带 $[r_\kappa^1,r_\kappa^2]$ , 其中 $r_\kappa^j=\kappa^{\alpha_j}e^{-2\kappa}$ 且 $0<\alpha_1<\alpha_2<1/4$ , 作者证明平衡态系数接近由渐近线性状态方程确定的常数态. 具体地, 对某个 $\delta>0$ ,

r^2\rho_\kappa\to B_*,\quad r^2p_\kappa\to A_*,\quad \frac{m_\kappa}{r}\to U_*,\quad r\mu_\kappa'\to\beta_*,

并且上述量以及 $e^{2\lambda_\kappa}$ 、相关导数和 $P'(\rho_\kappa)$ 的误差均为 $O(\kappa^{-\delta})$ . 常数由 $$c_*^2$$ 决定, 例如 $$D_*=1+6c_*^2+c_*^4$$ , $B_*=c_*^2/(2\pi D_*)$ , $A_*=(c_*^2)^2/(2\pi D_*)$ , $$U_*=2c_*^2/D_*$$ . 这一渐近结论用于构造负方向, 并非另一个独立的稳定性阈值公式.

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

证明的主要技术难点在于液体自由表面处的密度跃变. 气态星体的边界密度趋于零, 而这里流体侧密度等于 $\rho_0>0$ . 线性化后这一跃变进入表面项, 产生 Robin 边界条件; 因而不能直接照搬气态模型中带有零均值约束的算子框架. 作者的处理重点是保留该边界项, 建立适配它的加权线性化算子与变分形式, 并将负特征值与增长模态数对应起来. 摘录明确给出了算子具有离散实谱以及负 Morse 指数计数的结论, 但未展示完整函数空间定义、闭性证明和各项系数的详细推导, 因此这些技术细节从提供内容中无法进一步确认.

小红移稳定性通过控制线性化二次型完成. 文章将最小广义特征值表述为带权变分问题: 在 $\|\chi\|_{W_\kappa}=1$ 的归一化下考察 $\langle L_\kappa\chi,\chi\rangle_{\Lambda_\kappa}$ . 对充分小的中心红移, 作者证明该二次型对所有可容许扰动非负, 从而没有负特征值, 也就没有指数增长模态. 其关键不是单纯引用一般能量估计, 而是依托小红移平衡态的系数控制来获得对整个容许函数类的符号结论. 摘录未给出该控制的逐项估计, 故不宜补写其具体常数或误差阶.

大红移不稳定性的障碍是中心附近出现奇异高密度尺度, 直接以原变量估计背景解不足以判断二次型符号. 作者先利用渐近线性状态方程 $P(\rho)\sim c_*^2\rho$ , 比较真实平衡态与线性方程状态 $p=c_*^2\rho$ 对应的解. 在缩放变量下, 真实压强与比较压强之差由状态方程的有界余项控制; 文中给出相应比较估计, 再以 Gronwall 不等式传递到平衡态系数. 对比较模型, 令 $w_1=r^2\rho$ , $$w_2=m/r$$ , $\tau=\log r$ , 静态方程化为二维自治系统. 该系统的正平衡点为

Z_*=(B_*,U_*)=\left(\frac{c_*^2}{2\pi D_*},\frac{2c_*^2}{D_*}\right).

利用轨道有界性、Buchdahl 不等式、Dulac 判据和 Poincare-Bendixson 理论, 作者证明物理轨道趋于该平衡点, 并由线性化特征值得到任意 $0<\gamma<\eta_*$ 下的幂次衰减. 进一步在缩放环带上证明真实星体具有相同的主导渐近形状, 以此选取局部化试验函数, 使线性化二次型为负. 负方向推出负特征值, 因而出现增长模态. 提供摘录给出了这一证明路线与系数渐近, 但没有给出试验函数的完整表达式及其二次型逐项计算.

转折点原理的核心技巧是对整族平衡态求参数导数, 而不是在边界处寻找一个预设的零模. 固定径向变量 $$y$$ 对质量函数求导, 得到

v_\kappa(y)=-\frac{\partial_\kappa m_\kappa(y)}{4\pi y^3(\rho_\kappa(y)+p_\kappa(y))}.

对静态方程微分可知 $v_\kappa$ 在星体内部正则并满足零频线性化方程. 关键区别是: 它一般不满足自由表面的 Robin 条件. 将该条件写为 $D_\kappa f=0$ , 其中

D_\kappa f=yf'+(3-y\mu_\kappa')f,

对表面恒等式求导后, $v_\kappa$ 的边界缺陷正比于 $\partial_\kappa M_\kappa$ . 这就把“内部零解”与“真正算子核”精确地区分开来: 只有质量曲线临界时边界缺陷消失. 结合正则内部解的唯一性, 得到核维数刻画.

最后, 通过线性化算子的对称性和特征值随 $\kappa$ 的连续变化, 作者把核出现转化为唯一特征值穿越零的问题. 在非退化质量临界点附近, 质量导数的符号变化与半径导数共同决定穿越方向, 给出 $\operatorname{sgn}\nu_*=-\operatorname{sgn}[(\partial_\kappa M_\kappa)(\partial_\kappa R_\kappa)]$ . 这一步把谱变化与质量-半径曲线的弯折方向连接起来. 摘录未提供特征值导数的完整计算, 因而具体的边界配对推导无法从提供内容中确认.

4. 与之前工作的比较

本文与气态相对论星体的谱分析之间, 最明确的比较对象是 Hadzic、Lin 和 Rein 的工作 [6,7]. 据本文引言, [6] 建立了 Einstein-Euler 系统线性化动力学的谱分析并证明中心红移充分大时不稳定; [7] 建立了稳定性与质量-半径曲线之间的转折点原理. 两者研究的是气态模型, 密度在边界处连续趋零. 本文研究的液态模型在表面有 $\rho_0>0$ 的密度跃变, 因此线性化算子多出表面项并满足 Robin 边界条件. 作者指出, 气态框架中的零均值约束以及相应需要扣除的绕数指标在这里不适用; 取而代之的是构造新的正则内部零解, 再用它的边界缺陷识别质量临界点. 因此推进不只是把同一结论换到另一种状态方程, 还包括处理不同的自由边界谱结构.

就液态星的稳定性而言, 直接前身是 Hao 和 Miao [3] 对硬相模型的研究. 本文引言称 [3] 已证明硬相液态星在小中心红移下线性稳定、大中心红移下线性不稳定, 并在大密度区域使用局部化试验函数进行变分分析. 本文把稳定性结论推广到满足一般条件的 barotropic 状态方程, 其中高密度主导斜率为 $$c_*^2$$ , 同时要求正表面密度、正且因果的声速以及指定的高密度逼近性质. 大红移分析继承了“在高密度区域寻找负方向”的思路, 但需要额外证明一般状态方程的解能够在缩放环带上逼近渐近线性比较模型. 此外, 本文在该一般模型上建立了转折点原理, 而引言将硬相工作的重点描述为两端红移区域的稳定性判断及硬相模型的相关非线性不稳定性结果.

背景理论方面, 文中提及 Fournodavlos 和 Schlue [5] 构造硬相静态星体族, 并证明小星体在固定粒子数变分下局部最小化质量能量泛函、线性化能量一致有界; Miao 和 Shahshahani [8] 以及 Hao、Huo 和 Miao [9] 建立相关自由边界模型的适定性. 这些结果为动态自由边界模型提供基础, 但本文处理的是特定平衡态附近的径向线性谱问题, 其结论不能与自由边界非线性初值问题的长时间适定性或非线性稳定性混为一谈. 对硬相情形, 引言称 [3] 还结合 [8] 的非线性估计得到大密度区域的非线性不稳定性; 本文摘要和所给定理则明确陈述线性增长模态及谱转折规律, 并未声称一般状态方程下相同的非线性结论.

转折点思想本身可追溯至 Zel’dovich、Wheeler 以及 Harrison、Thorne、Wakano 和 Wheeler 的相对论恒星研究, 其基本图像是稳定模态数只在总质量沿平衡态族达到极值时改变, 而质量-半径曲线的弯折方向决定增减. 相关 Newtonian 情形的工作在文中列为 [10,11,12,17]. 本文的贡献是为具有密度跃变和 Robin 自由表面条件的液态模型给出与该图像相符的谱论证明, 尤其把“质量临界”严格落实为核非平凡, 并把穿越方向写成 $\partial_\kappa M_\kappa$ 与 $\partial_\kappa R_\kappa$ 的乘积符号. 这比只用质量曲线作经验判据更具体, 但定理所陈述的方向变化结论针对非退化质量临界点; 摘录没有给出退化临界点处模态数变化的分类.

适用范围与代价也需区分. 本文覆盖的状态方程比单一硬相律更一般, 但仍依赖严格的单调性、因果声速界、正表面密度和特定高密度渐近条件; 不能据此推断任意状态方程或非球对称扰动都适用. 分析限定于球对称径向谱稳定性, 不直接处理非径向模态. 大红移结论给出“充分大”时存在增长模态, 小红移结论给出“充分小”时无增长模态, 但提供内容没有显式给出阈值, 也没有说明中间参数区间的全部谱结构. 质量转折结论刻画模态数发生变化的必要位置与非退化穿越方向, 本身不提供每个平衡态的绝对不稳定模态数; 还需结合某个参数区间的基准模态数才能确定整条平衡态族上的计数. 文中对这些限制之外更广泛的适用性未作明确证明, 不应从所给摘录外推.

5. 具有超线性 BD 黏性和远场真空的三维大振幅球对称退化可压缩 Navier-Stokes 方程的整体正则性

Global regularity for 3D large spherically symmetric degenerate compressible Navier-Stokes equations with superlinear BD viscosity and far-field vacuum

作者: Xiangdi Huang, Bingshu Li

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32062v1

摘要翻译: 超线性密度依赖黏性给可压缩 Navier-Stokes 方程的整体正则性带来了根本障碍，尤其是在存在远场真空时：黏性扩散在无穷远处消失，因此，仅排除有界区域上的真空不足以在全局控制流动。对于 $α\in (1,11.7)$ 和 $γ\in [2α-1,3α-1)$ ，我们证明了任意大振幅、强、有限质量初值所对应解的整体正则性，其中初值在每个有限点处密度为正，并在无穷远处处于真空状态。首先，我们引入关于速度、有效速度和质量的新缩放变量，随后利用中心附近的若干新型量，并采用精细的最大值原理，证明速度在中心处具有 Lipschitz 连续性，从而消除速度梯度的奇异性。另一个关键机制是对物质体积畸变进行定量控制，且这种控制在真空尾部仍然有效。对物理速度梯度和有效速度梯度的保守估计给出了在无穷远处一致趋于零的变差尾部。这些估计确定了一条累积 Jacobian 变差较小的物质轨迹，并导出外部体积 Jacobian 的双侧界，从而使沿粒子路径演化的密度与其初始值可比。在由初始密度确定的空间尺度上，所得退化方程具有一致的抛物型估计，无需全局密度下界即可恢复对物理形变的控制。结合紧集上的正性估计和单侧 Glaeser 不等式，这闭合了一个齐次延拓论证，而无需对密度幂本身施加 $$L^2$$ -可积性。解保持质量守恒和远场真空，并且在每个紧时空柱上均保持光滑且严格为正。

Comment: 66 pages

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

1. 研究对象和背景

本文研究三维全空间中的等熵可压缩 Navier–Stokes 方程. 未知量为密度 $\rho(t,x)$ 和速度 $$u(t,x)$$ , 压力取 $P(\rho)=\rho^\gamma$ . 黏性系数依赖密度, 并满足 Bresch–Desjardins (BD) 结构关系 $\mu(\rho)=\rho^\alpha$ , $\lambda(\rho)=(\alpha-1)\rho^\alpha$ . 论文采用如下对称梯度应力约定:

\begin{cases} \rho_t+\operatorname{div}(\rho u)=0,\\ (\rho u)_t+\operatorname{div}(\rho u\otimes u)+\nabla \rho^\gamma =\operatorname{div}\left[\dfrac{\rho^\alpha}{2}\bigl(\nabla u+(\nabla u)^T\bigr)+(\alpha-1)\rho^\alpha(\operatorname{div}u)I\right]. \end{cases}

初值满足 $(\rho,u)|_{t=0}=(\rho_0,u_0)$ , 并要求远场状态为真空:

(\rho,u)(t,x)\longrightarrow (0,0)\qquad (|x|\to\infty).

研究限于球对称解 $\rho(t,x)=\rho(t,r)$ , $$u(t,x)=u(t,r)x/r$$ , $$r=|x|$$ . 初始密度在每个有限位置严格为正, 但具有有限总质量且趋于零于无穷远. 因而这里的真空是远场真空, 不是允许初始时有限位置出现真空. 球心处还需满足与光滑笛卡尔向量场相容的奇偶条件: $\rho$ 关于径向变量为偶函数, $$u$$ 为奇函数, 特别是 $$u(t,0)=0$$ .

问题的重要性在于, 对 $\alpha>1$ , 单位质量上的有效黏性扩散系数为

\frac{\mu(\rho)}{\rho}=\rho^{\alpha-1}\longrightarrow 0\qquad (\rho\to0).

所以即使能在每个有界区域内证明密度有正下界, 也不能据此获得无穷远区域的统一抛物性或速度梯度控制. 远场密度趋零时, 方程在整个空间上退化; 全局正则性不能仅靠局部排除真空来完成. 另一方面, 球对称方程包含 $$u/r$$ 等几何项, 球心处径向速度梯度的控制也不能直接从远离原点的估计推出. 论文要同时处理球心的几何奇性和远场退化, 并在大数据、有限质量及无穷远真空的条件下阻止有限时间内的奇性和有限位置真空形成.

在径向变量中, 令 $$d=u_r+2u/r$$ , 则连续性方程为 $\rho_t+(\rho u)_r+2\rho u/r=0$ . 论文还引入有效速度

v=u+\frac{\alpha}{\alpha-1}\nabla \rho^{\alpha-1},\qquad v=u+\alpha\rho^{\alpha-2}\rho_r

(第二式为径向表达), 从而

\partial_r\rho^{\alpha-1}=\frac{\alpha-1}{\alpha}(v-u).

BD 结构使有效速度具有额外能量估计, 但在 $\alpha>1$ 时, 控制 $\rho^{\alpha-1}$ 的空间变化并不自动排除有限位置密度降为零, 因为该幂在 $\rho\downarrow0$ 时仍可保持有限. 因此需要把有效速度关系与流映射的演化结合起来控制密度.

已有理论提供了若干背景, 但与本文设定并不等价. 常黏性小扰动理论和允许真空的大数据弱解理论, 本身不保证强解的点态密度界及导数正则性. BD 熵结构已用于退化系统的弱解存在, 但弱解框架也不直接给出本文所需的强正则性延拓控制. 对径向强解, 既有工作主要覆盖有界区域、常黏性或线性及次线性密度黏性; 全空间真空问题的相关整体结果则采用小初始密度并配合膨胀速度, 或要求远场有正密度背景. 因此, 本文试图填补的具体空白是: 对超线性 BD 黏性 $\alpha>1$ , 在全空间远场真空、有限质量且初值可任意大的径向情形下, 构造并延拓整体强解, 同时不假设全局正密度下界或预设密度衰减率.

2. 主要定理或主要结果

核心结论是定理 1.2. 假设参数满足

1<\alpha<11.7,\qquad 2\alpha-1\leq\gamma<3\alpha-1.

初值为球对称强数据, 在每个有限球上具有与光滑笛卡尔中心延拓相容的正则性, 并满足

0<\rho_0\in L^1\cap L^\infty,\qquad \rho_0(x)\to0\ (|x|\to\infty),\qquad u_0\in H^3,

\nabla\log\rho_0\in L^\infty,\qquad \nabla^2\log\rho_0,\ \nabla^3\log\rho_0\in L^2,

以及 $\nabla\rho_0^{\gamma-1}\in H^2$ . 则方程存在唯一的整体球对称解; 对任意有限 $$T$$ , 解属于文中定义的齐次解类 $$X_T$$ , 且在每个有限时空柱体上为强解, 包括球心处. 唯一性是在具有相同原始初值、且属于该强解类的径向解中成立.

定理给出的定量点态结论是: 对任意有限 $$T$$ 和 $$R$$ , 存在依赖于数据、 $$T$$ 、 $$R$$ 的常数 $C_T<\infty$ 和 $c_{T,R}>0$ , 使得

\sup_{0\leq t\leq T,\ x\in\mathbb R^3}\rho(t,x)\leq C_T, \qquad \inf_{0\leq t\leq T,\ |x|\leq R}\rho(t,x)\geq c_{T,R}.

因此没有有限时间内有限位置的真空形成, 但定理不声称存在全空间统一的正密度下界. 同时

\rho\in C([0,T];L^1),\qquad \int_{\mathbb R^3}\rho(t,x)\,dx=\int_{\mathbb R^3}\rho_0(x)\,dx,

并且 $\rho,u,v\to0$ 在 $t\in[0,T]$ 上对远场一致成立. 数据大小不受小量条件限制.

证明中还得到下列用于延拓的关键估计. 令 $$r(t,a)$$ 为初始半径 $$a$$ 出发的粒子轨迹, 满足 $\partial_t r=u(t,r)$ , $$r(0,a)=a$$ , 并令

J(t,a)=\frac{r(t,a)^2\partial_a r(t,a)}{a^2}.

连续性方程给出 $\rho(t,r(t,a))J(t,a)=\rho_0(a)$ , 且 $\partial_tJ=J(u_r+2u/r)(t,r(t,a))$ . 作者证明在充分大的初始半径上, $$J$$ 有与半径无关的双边界 $C_T^{-1}\leq J(t,a)\leq C_T$ ; 因此沿相应粒子轨迹, 当前密度与初始密度可比较. 另外, 在有限延拓区间上得到外部变差尾部消失:

\lim_{R\to\infty}\sup_{0\leq t<T}\int_R^\infty\bigl(|u_r(t,r)|+|v_r(t,r)|\bigr)\,dr=0.

这一估计用于选取累计雅可比变化较小的外部物质轨迹, 而不是直接预设远场密度的衰减速度.

作者使用的齐次解类不要求 $\rho^{(\alpha-1)/2}$ 或 $\rho^{(\gamma-1)/2}$ 属于 $L^2(\mathbb R^3)$ , 也不要求它们的一阶导数具有相应的全局平方可积性. 文中能量量 $$E$$ 只纳入 $$u$$ 的 $$H^3$$ 范数、 $\nabla\rho^{\gamma-1}$ 及两种密度幂的二、三阶导数等项, 不含两种密度幂本身的 $$L^2$$ 范数或其一阶导数的 $$L^2$$ 范数. 这使整体结论适用于仅有有限质量、但密度幂未必平方可积的情形.

摘要和定理陈述支持上述整体存在、唯一性、质量守恒、局部正性、远场衰减及外部控制结论. PDF 摘录只提供到第 10 页, 未包含后续各节的完整引理和证明细节; 因此不能从所给内容进一步确认外部抛物估计及 Glaeser 不等式的全部具体常数、范数指数或各中间命题编号.

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

证明的主要困难不是获得普通能量界, 而是把这些界转化为整个无界区域上的速度变形控制. 在远场, $\rho^{\alpha-1}$ 使单位质量扩散退化; 球心处则有 $$u/r$$ 、 $$u_r-u/r$$ 等径向几何结构. 作者将证明分成球心和有界区域的正性控制、外部流映射控制、密度依赖尺度上的抛物正则性三部分, 最后通过延拓准则闭合.

首先利用普通能量和 BD 有效能量. 普通能量控制 $\int(\rho|u|^2/2+\rho^\gamma/(\gamma-1))$ 以及黏性耗散; 有效速度能量进一步给出密度梯度型耗散. 摘录中的两条恒等式为

\int_{\mathbb R^3}\left(\frac{\rho|u|^2}{2}+\frac{\rho^\gamma}{\gamma-1}\right)(t)dx +\int_0^t\!\int_{\mathbb R^3}\rho^\alpha\left(|\nabla u|^2+(\alpha-1)(\operatorname{div}u)^2\right)dxds=E_{\rm ord}(0),

\int_{\mathbb R^3}\left(\frac{\rho|v|^2}{2}+\frac{\rho^\gamma}{\gamma-1}\right)(t)dx +\alpha\gamma\int_0^t\!\int_{\mathbb R^3}\rho^{\alpha+\gamma-3}|\nabla\rho|^2dxds=E_{\rm BD}(0).

这些估计提供 $\rho^\gamma$ 、有效速度及密度梯度耗散的控制, 但不直接给出 $\rho$ 的全局上界或正下界. 为得到密度上界, 作者借用并改造 Huang–Meng–Zhang 的径向矩估计, 在全空间外部处理压力项, 同时控制物理速度和有效速度的高阶径向矩, 再借助 $$v-u$$ 与 $\partial_r\rho^{\alpha-1}$ 的关系回收密度上界. 所给摘录明确显示, 选择矩指数时需满足

4p^2(\alpha-1)^2<(2p-1)\alpha(4\alpha-2),\qquad p>q>\frac32.

在 $$p=3/2$$ 处相应的余量为 $-\alpha^2+14\alpha-9$ , 这为参数上界 $\alpha<11.7$ 提供了文中展示的矩估计来源之一. 摘录未给出该节各矩估计的完整推导, 不宜将这一条件解释成唯一的参数限制来源.

其次, 有界区域上的正性与球心正则性需要避免径向奇性. 文中以质量坐标及物理、有效速度的耦合比较估计控制球内密度, 并对倒数密度作物质导数估计. 球心附近作者引入速度、有效速度和质量的重标度变量, 并构造适配球心的量, 再通过精细最大值原理证明速度在中心具有 Lipschitz 连续性. 其作用是控制速度梯度的球心极限, 而不是把径向坐标中的 $$1/r$$ 项当作普通有界系数处理. 摘录未列出这些新变量的具体定义及最大值函数, 因而其精确形式从提供内容中无法确认.

第三, 远场的核心创新是对物质体积畸变作定量控制. 定义流映射的体积雅可比 $$J$$ , 则 $\rho J=\rho_0$ 沿轨迹成立. 但 $$J$$ 的方程含有物理速度散度, 在证明外部梯度控制之前不能先验积分获得双边界. 作者因此推导 $$u_r$$ 与 $$v_r$$ 的保守型方程及变差估计, 证明远场的径向变差尾部对时间一致趋零. 由此选择一条累计雅可比变动较小的外部物质轨迹, 再结合雅可比方程和归一化有效应变, 先在大初始半径区域建立 $$J$$ 的双边控制. 这一步把“无法先控制远场速度梯度”的循环依赖改造成先控制外部变差尾部、再确定密度与初值可比、最后恢复抛物控制的顺序.

第四, 有了沿物质轨迹的密度可比性后, 作者按初始密度决定空间尺度. 在大初始半径 $$a_c$$ 附近取

\ell=\rho_0(a_c)^{(\alpha-1)/2},\qquad a=a_c+\ell y,

并保持时间变量不变. 在这一尺度上, 关于 $J^{-\alpha}$ 的方程具有一致有界且远离零的扩散系数, 从而可使用一维内部抛物 $$L^p$$ 估计及 Schauder 估计. 这些估计的常数不依赖中心位置 $$y_0$$ 或向外趋于无穷的终止时刻; 这使估计可以在密度越来越小、尺度 $\ell$ 趋于零的远场保持统一. 抛物估计先控制 $J^{-\alpha}$ 的时间导数和空间二阶导数, 再结合连续性方程把控制转回物理应变 $$u_r+2u/r$$ , 由此恢复无须全局密度下界的变形控制.

最后, 局部紧致区域上的正性估计与外部估计拼接, 得到 $\int_0^T\|\nabla u(t)\|_\infty^2dt<\infty$ . 再将其与 $\rho^{\alpha-1}$ 的单侧 Glaeser 型不等式结合, 控制延拓准则所需的密度幂梯度, 完成齐次强解类中的延拓. 关键之处是能量与延拓空间均不强加密度幂本身的 $$L^2$$ 可积性. 摘录给出了上述证明框架及抛物估计类型, 但未展示外部变差方程、Glaeser 不等式及最终范数估计的完整计算, 对这些中间环节的更细公式从提供内容中无法确认.

4. 与之前工作的比较

论文在引言中把结果放在常黏性整体理论、BD 退化系统弱解理论、全空间真空强解及径向强解几条线上比较. 与常黏性情形相比, Merle–Raphaël–Rodnianski–Szeftel 构造了特定数据下的球心内爆, 说明强解未必自动整体正则. 本文并未否定该类机制, 而是在不同的密度依赖黏性结构 $\mu=\rho^\alpha$ , $\lambda=(\alpha-1)\rho^\alpha$ 和指定参数范围内证明整体正则; 作者也提及 Chen–Liu–Zhu 在某些次线性参数区域构造非线性黏性内爆. 因而结论应理解为特定 BD 超线性模型下的整体结果, 不是对一般可压缩黏性流无奇性问题的普遍回答.

与 $\alpha\leq1$ 的相关径向强解工作相比, 差异在于黏性是否超线性以及远场扩散的性质. Huang–Meng–Zhang 在三维有界球上处理任意大初值的次线性区间并建立强解不爆破结果; Zhang 及 Lei 的结果也涉及次线性指数范围, 其中 Lei 还考虑全空间非真空远场密度. Chen–Zhang–Zhu 对全空间无穷远真空的线性黏性情形证明了大径向数据整体正则, 三维压力范围为 $1<\gamma<3$ ; 其后续工作处理非线性次线性黏性. 这些结果均没有覆盖本文 $\alpha>1$ 的远场退化机制. 在超线性情形, $\mu(\rho)/\rho=\rho^{\alpha-1}$ 随真空趋于零, 因而有界区域的正性不能替代外部控制; 本文的实质推进是引入流映射体积畸变及按初始密度选定的外部抛物尺度, 而非简单延用线性或次线性的局部估计.

与已有全空间远场真空整体强解结果相比, 本文不要求初始密度足够小并配合膨胀型初始速度. 引言列出的 Xin–Zhu 结果允许速度振幅较大, 但要求小初始密度和膨胀初始速度; Huang–Li–Zhang 的整体强解结果则要求初始密度与足够大的正常数可比. 本文初始密度具有有限质量并趋于零, 可在每个有限位置为正但没有全局正下界, 且数据可以任意大. 代价是必须限制在球对称类, 并采用明确的初始 Sobolev 与对数导数条件及参数区间. 因此它拓展的是对称性受限但远场退化更强、初值不要求小或膨胀的整体强解理论, 而不是一般多维非对称大数据理论.

与弱解存在理论相比, Lions、Feireisl–Novotný–Petzeltová 以及后续 BD 结构相关工作处理了较宽泛的大数据和真空情形, 但弱解存在不等于点态正则性、有限位置严格正性或强解唯一性. Bresch–Vasseur–Yu 等对物理对称应力和非线性黏性律的弱解结果涵盖了广泛的黏性幂次, 其结论类型与本文不同. 本文获得每个有限时空柱体上的强解正则性、有限位置密度正下界、质量守恒和强解类内唯一性, 但只对径向数据成立, 且需初始密度在有限点严格为正. 它不覆盖初始时有限位置已有真空的强解问题.

参数和函数空间方面, 本文要求 $1<\alpha<11.7$ 及 $2\alpha-1\leq\gamma<3\alpha-1$ . 文中展示的径向矩相干性条件解释了 $\alpha$ 上界的一个来源; 压力指数范围则是定理的明确限制. 结果不声称这些区间最优, 所给摘录也没有给出端点不可达或反例, 故不能据此断言范围尖锐. 与常见有限能量空间不同, 作者特意采用齐次密度幂解类, 不要求 $\rho^{(\alpha-1)/2}$ 、 $\rho^{(\gamma-1)/2}$ 本身属于 $$L^2$$ , 这避免把有限质量不蕴含的平方可积性当作隐藏假设. 另一方面, 初始对数导数条件 $\nabla\log\rho_0\in L^\infty$ 、 $\nabla^2\log\rho_0,\nabla^3\log\rho_0\in L^2$ 仍是明确代价; 作者特别说明这些全局对数范数的传播不是定理结论, 延拓依靠的是其他范数控制.

总的说, 文章相对既有工作的推进集中在机制而非单纯扩大指数区间: 将球心 Lipschitz 控制、外部速度变差尾部、物质体积雅可比双边界、密度决定的内禀抛物尺度以及单侧 Glaeser 控制连成闭合的延拓链条. 其局限也清楚: 结论依赖球对称性, 不给出全空间正密度下界, 不处理初始有限位置真空, 且参数范围及初值正则性有要求. PDF 摘录给出了引言中的文献比较, 但未提供每篇前作全部假设与定理的完整陈述; 除文中明确写出的范围外, 更细的技术条件差异从提供内容中无法确认.

6. 具有局部化扰动的 Stokes 与 Navier-Stokes 系统的快速稳定化

Rapid stabilization for Stokes and Navier-Stokes systems with localized disturbance

作者: Christian Calle, Patricio Guzmán, Hugo Parada, Jefferson Prada

类别: math.OC, math.AP

发布日期: 2026-10-22

arXiv ID: 2610.31790v1

摘要翻译: 本文研究存在未知空间局部化扰动时 Stokes 系统以及二维和三维 Navier-Stokes 系统的快速稳定化。我们提出一种多值反馈控制策略，将 Lyapunov 频率方法与多值符号算子相结合。该方法将基于 Lyapunov 的稳定性分析与谱不等式联系起来，而引入符号算子则提供了对扰动的鲁棒性。闭环系统由一个微分包含所支配；我们借助极大单调算子理论证明了该系统的适定性。对于任意大的扰动界，该方法保证了 Stokes 系统的全局快速指数稳定，以及二维和三维 Navier-Stokes 系统的局部快速指数稳定。在三维情形下，对于在 $$L^2$$ -norm 中充分小、且在 $$H^1$$ -norm 中具有预先给定的任意大界的初值，我们获得了全局强解和快速指数稳定性。

Comment: 无

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

1. 研究对象和背景

本文研究有界光滑连通区域 $\Omega\subset\mathbb{R}^n$ 中不可压缩流体的内部反馈稳定化. 控制只作用于非空开子集 $\omega\subset\Omega$ , 同一作用通道上还存在未知扰动. 速度为 $$y$$ , 压力为 $$p$$ , 控制为 $$u$$ , 扰动为 $$d$$ , 初值为 $$y_0$$ . 文中考虑的 Stokes 系统和 Navier–Stokes 系统分别为

\begin{cases} \partial_t y-\Delta y+\nabla p=\mathbf{1}_\omega(u+d), & \Omega\times(0,\infty),\\ \nabla\cdot y=0, & \Omega\times(0,\infty),\\ y=0, & \partial\Omega\times(0,\infty),\\ y(0)=y_0, & \Omega, \end{cases}

\begin{cases} \partial_t y-\Delta y+(y\cdot\nabla)y+\nabla p=\mathbf{1}_\omega(u+d), & \Omega\times(0,\infty),\\ \nabla\cdot y=0, & \Omega\times(0,\infty),\\ y=0, & \partial\Omega\times(0,\infty),\\ y(0)=y_0, & \Omega. \end{cases}

其中 $\mathbf{1}_\omega$ 表示将控制域上的函数延拓为零后嵌入整个区域. 扰动与控制通道匹配, 但其具体取值未知; 控制器只使用已知上界 $$D>0$$ . 假设

d\in L^1_{\mathrm{loc}}(0,\infty;L^2(\omega)),\qquad \|d(\cdot,t)\|_{L^2(\omega)}\le D\quad\text{对几乎处处的 }t>0.

因此本文并不要求扰动随时间衰减, 也不要求扰动振幅很小. 目标是在仅知道 $$D$$ 而不知道 $$d$$ 的情况下, 构造依赖状态的反馈, 使速度以预先指定的速率指数衰减. “快速”在此指衰减参数 $\lambda>0$ 可任意预先指定, 而不是只证明某个由系统耗散决定的固定衰减率.

分析使用标准的无散空间 $H=\overline{\{v\in C_0^\infty(\Omega)^n:\nabla\cdot v=0\}}^{L^2}$ 与 $V=\overline{\{v\in C_0^\infty(\Omega)^n:\nabla\cdot v=0\}}^{H_0^1}$ . Leray 投影记作 $$P$$ , Stokes 算子为 $S=-P\Delta$ , 其定义域为 $D(S)=H^2(\Omega)^n\cap V$ . Stokes 算子具有紧预解式, 因而可取 $$H$$ 的正交特征基 $\{e_j\}$ , 满足 $Se_j=\tau_j e_j$ , 且 $0<\tau_1\le\tau_2\le\cdots$ . 对给定 $\lambda$ , 取低频空间 $H_N=\operatorname{span}\{e_j:\tau_j\le\lambda\}$ , 并记相应正交投影为 $$P_N$$ .

问题的重要性在于, 低频不稳定或衰减缓慢的分量需要通过局部控制迅速压制, 但局部控制区域并不覆盖整个流体域; 同时未知扰动可持续作用, 不能简单地作为已知外力抵消. 对 Navier–Stokes 系统还需处理无散约束、压力以及二次对流项. 三维情形尤其受到强解全局适定性问题的限制. 既有频率 Lyapunov 方法能借助局部观测的谱不等式构造任意指定衰减率的有限维反馈, 但如何将其与对未知且不小扰动的鲁棒性结合, 并保持非线性闭环系统适定, 是本文试图推进的核心空白. 摘录明确给出的贡献是把频率 Lyapunov 设计与多值符号反馈结合, 并据此处理 Stokes 及二维、三维 Navier–Stokes 系统; 对证明细节超出所提供 PDF 摘录的部分, 下文不作推断.

2. 主要定理或主要结果

记 $C_1\ge1$ 为谱不等式常数, 其只依赖于 $(\Omega,\omega)$ . 对每个 $\lambda>0$ , 文中给出反馈

G_\lambda(y)=-\gamma_\lambda P_Ny-D\,\operatorname{sign}_{L^2(\Omega)}(\mathbf{1}_\omega C_\lambda y),

其中

\gamma_\lambda=\lambda C_1e^{C_1\sqrt{\lambda}},\qquad \mu_\lambda=\left(\frac{\gamma_\lambda}{\lambda}\right)^2=C_1^2e^{2C_1\sqrt{\lambda}},\qquad C_\lambda y=\mu_\lambda P_Ny+P_N^\perp y.

符号算子是多值映射: 对 $f\ne0$ , $\operatorname{sign}_{L^2}(f)=\{f/\|f\|_{L^2}\}$ ; 对 $$f=0$$ , $\operatorname{sign}_{L^2}(0)=\{\eta:\|\eta\|_{L^2}\le1\}$ . 这里的反馈写法与论文中的算子记号一致, 其中作为权重的 $C_\lambda$ 由低频加权定义, 而稳定性估计中的常数也记作 $C_\lambda$ ; 摘录将二者同记, 具体语境可区分.

定理 1.1 (Stokes). 对任意 $n\ge2$ , $\lambda>0$ , 任意满足上述假设的扰动和任意 $y_0\in H$ , 闭环 Stokes 系统存在唯一广义解

y\in C([0,\infty);H)\cap L^2_{\mathrm{loc}}(0,\infty;V),

并满足

\|y(t)\|_H\le C_\lambda e^{-\lambda t/4}\|y_0\|_H,\qquad t\ge0,

其中稳定性常数可取 $C_\lambda=\sqrt{\mu_\lambda}=C_1e^{C_1\sqrt{\lambda}}$ . 因而结论对任意大小的扰动上界 $$D$$ 成立, 且初值不需要小性条件.

定理 1.2 (二维 Navier–Stokes). 对 $$n=2$$ 和任意 $\lambda>0$ , 存在 $C_\lambda\ge1$ 及 $\rho_\lambda>0$ . 对每个满足扰动假设的 $$d$$ 以及满足 $\|y_0\|_H\le\rho_\lambda$ 的初值, 闭环系统存在唯一广义解

y\in C([0,\infty);H)\cap L^2_{\mathrm{loc}}(0,\infty;V),

并且同样有 $\|y(t)\|_H\le C_\lambda e^{-\lambda t/4}\|y_0\|_H$ . 这是局部于初值的快速稳定结论; 摘录中的定理陈述没有说明 $\rho_\lambda$ 对 $$D$$ 的依赖, 因而不能据此补充其依赖关系.

定理 1.3 (三维 Navier–Stokes). 对 $$n=3$$ , 任取 $\lambda>0$ , $$D>0$$ 和 $$R>0$$ , 存在 $C_\lambda\ge1$ 及 $\rho_{\lambda,D,R}>0$ . 若 $y_0\in V$ 且

\|y_0\|_V\le\frac{R}{2},\qquad \|y_0\|_H\le\rho_{\lambda,D,R},

则闭环系统有唯一全局强解, 对每个 $$T>0$$ ,

y\in H^1(0,T;H)\cap L^2(0,T;D(S))\cap C([0,T];V),

并满足

\|y(t)\|_H\le C_\lambda e^{-\lambda t/4}\|y_0\|_H,\qquad \sup_{t\ge0}\|y(t)\|_V\le\frac{R}{\sqrt{2}}.

该结论允许预先指定任意大的 $$H^1$$ 初值上界 $$R$$ , 但要求 $$L^2$$ 初值足够小; 小性阈值依赖 $\lambda,D,R$ . 这不是对任意三维初值的全局强解定理. 三个定理中的指数衰减率均为 $\lambda/4$ , 而常数 $C_\lambda=C_1e^{C_1\sqrt{\lambda}}$ 随指定速率增长. 摘录未给出 $\rho_{\lambda,D,R}$ 的显式公式或其定量尺度, 因而不能进一步报告该阈值的误差阶或增长率.

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

证明设计的主线是将两类功能不同的反馈叠加: 低频线性反馈产生快速衰减, 多值符号反馈消除未知扰动在 Lyapunov 能量中的正贡献. 核心困难不是单纯控制线性半群, 而是控制只能局部作用, 同时扰动未知且闭环控制含有不连续的多值算子; 在 Navier–Stokes 情形还需证明非线性解不会离开适用的能量区域.

首先, 作者以 Stokes 特征分解选取所有满足 $\tau_j\le\lambda$ 的低频模态. 低频部分通过谱不等式获得局部可控性:

\|P_Ny\|_{L^2(\omega)}^2\ge C_1^{-1}e^{-C_1\sqrt{\lambda}}\|P_Ny\|_H^2.

反馈增益取 $\gamma_\lambda=\lambda C_1e^{C_1\sqrt{\lambda}}$ , 恰使低频阻尼与上述谱不等式结合后给出量级为 $\lambda$ 的耗散. 高低频在局部观测区域中并不正交, 所以不能只对低频方程做估计; 证明还估计高频方程中由 $-\gamma_\lambda P_Ny$ 引起的交叉项, 并用高频特征值满足 $\tau_j>\lambda$ 所给出的 $\|P_N^\perp y\|_V^2\ge\lambda\|P_N^\perp y\|_H^2$ 吸收该项.

为组织这些估计, 论文使用加权 Lyapunov 泛函

V(y)=\mu_\lambda\|P_Ny\|_H^2+\|P_N^\perp y\|_H^2=\|y\|_\mu^2,

其中低频权重 $\mu_\lambda=(\gamma_\lambda/\lambda)^2$ . 该权重使低频阻尼、高频耗散以及低高频耦合项能够按同一尺度比较. 由谱不等式和 Young 不等式, 对强解可导出文中式 (3.4) 的能量不等式

\frac{d}{dt}V(y(t))\le-\frac12\|y\|_V^2-\frac{\lambda}{2}V(y(t)) +2(\mathbf{1}_\omega C_\lambda y,\,\widehat u+d)_{L^2(\Omega)}.

取 $\widehat u=-D\operatorname{sign}_{L^2}(\mathbf{1}_\omega C_\lambda y)$ . 对任意 $\eta\in\operatorname{sign}_{L^2}(\mathbf{1}_\omega C_\lambda y)$ , 有

(\eta,\mathbf{1}_\omega C_\lambda y)_{L^2}=\|\mathbf{1}_\omega C_\lambda y\|_{L^2},\qquad \|\eta\|_{L^2}\le1.

另一方面, $\|d(t)\|_{L^2(\omega)}\le D$ . 因而符号反馈贡献至多为 $-D\|\mathbf{1}_\omega C_\lambda y\|_{L^2}$ , 扰动贡献至多为 $D\|\mathbf{1}_\omega C_\lambda y\|_{L^2}$ , 两者在该能量配对中抵消, 得到

\frac{d}{dt}V(y(t))\le-\frac12\|y\|_V^2-\frac{\lambda}{2}V(y(t)).

再利用 $\|y\|_H^2\le V(y)\le\mu_\lambda\|y\|_H^2$ , 得到衰减率 $\lambda/4$ 及常数 $\sqrt{\mu_\lambda}$ . 这里“抵消”是 Lyapunov 配对中的鲁棒抵消, 不是声称反馈逐点重构或精确辨识未知函数 $$d$$ .

第二个障碍是多值反馈使闭环方程成为微分包含, 因此经典光滑反馈方程的适定性理论不能直接照搬. 文中把方程投影到无散空间, 并用包含 $\operatorname{sign}_{L^2}$ 的算子表达闭环系统; 摘录明确说明 Stokes 闭环系统的适定性通过极大单调算子理论建立, 并通过对正则初值和正则扰动的逼近, 将强解上的先验衰减估计传递到广义解. 多值符号图的单调性及强收敛—弱收敛闭性是所列预备性质, 与该框架相配合. 但所提供内容只到第 10 页, 未包含第 4 节的算子定义、极大单调性验证及连续依赖证明细节; 对这些证明步骤的具体形式, 从提供内容中无法确认.

第三个障碍来自对流项 $$B(y,y)$$ . 该项在未加权的 $$H$$ 能量测试中具有经典消去性质 $$b(y,y,y)=0$$ , 但加权测试 $C_\lambda y$ 一般不再完全消去. 作者给出的关键估计为

|\langle B(y),C_\lambda y\rangle_{V',V}\le a_\lambda\|y\|_H^3,

其中 $a_\lambda=(\mu_\lambda-1)M_\lambda$ , 而 $M_\lambda$ 控制有限维低频部分的 $\|\nabla P_Ny\|_{L^\infty}$ . 该估计利用低频子空间有限维、特征函数光滑以及对流项的反对称性, 把加权造成的误差压到仅依赖 $$H$$ 范数的三次项. 二维中, 通过小 $$H$$ 初值和相应先验估计控制该项. 摘录未提供第 5 节的完整不等式推导, 因此不能更精确地给出二维闭合条件.

三维证明还需同时控制较弱的 $$H$$ 范数和较强的 $$V$$ 范数. 作者引入径向截断 $R_\delta y$ 及截断非线性 $B_\delta(y)=B(R_\delta y,R_\delta y)$ , 先构造截断系统的强解, 再证明解不会达到截断阈值, 从而去掉截断. 摘录明确指出, 在三维中结合 $$H$$ 衰减与 $$V$$ 有界性, 并在高阶能量估计中保留完整 Stokes 耗散; 这是在允许任意给定的 $$V$$ 上界及扰动界时, 仍从足够小的 $$H$$ 初值推出全局强解和快速稳定的关键策略. 具体高阶微分不等式、截断阈值的选择及小性常数的表达式未在提供页中出现, 不宜补写.

总体上, 本文可确认的技术创新不只是“使用能量估计”, 而是把谱不等式决定的低频快速阻尼、针对有界匹配扰动的多值符号项、以及为 Navier–Stokes 非线性和三维强解所需的截断与双范数控制组合为一个闭环论证. 其代价也反映在反馈高增益、衰减常数对 $\lambda$ 的指数型依赖, 以及非线性情形的初值限制上.

4. 与之前工作的比较

本文的直接线性稳定化背景是 Xiang 的频率 Lyapunov 方法 [52, 53]. 该方法将 Lebeau–Robbiano 型谱不等式与作用于有限维低频空间的反馈结合, 以获得任意指定的指数衰减率; 本文沿用这一低频设计, 在 Stokes 算子特征模态上构造 $-\gamma_\lambda P_Ny$ , 并用加权 Lyapunov 泛函处理局部控制导致的低高频耦合. 因此对 Stokes 部分的推进, 主要不是另造谱稳定化机制, 而是在该机制中加入未知局部匹配扰动, 同时保留可指定的快速衰减率.

鲁棒性设计与文中明确提到的多维热方程工作 [22] 最为接近. 作者说明本文受 [22, 52] 启发, 将频率 Lyapunov 方法与多值符号算子结合, 从热方程推进到具有向量值无散约束、压力和 Stokes 算子的流体方程, 并进一步处理 Navier–Stokes 对流非线性. 与只讨论已知控制设计或一般扰动衰减的路线相比, 本文针对的是控制通道上的未知局部扰动, 只使用其先验界 $$D$$ , 不要求扰动小或趋于零. 但依据摘录, 不应把能量配对抵消夸大为对扰动函数的逐点精确辨识; 严格陈述是闭环解满足与无扰动型相同的快速衰减估计.

与滑模控制、主动扰动抑制及内模原理等工作相比, 本文所述问题设定各有侧重. 文中将滑模控制 [7, 29, 32] 作为匹配扰动鲁棒性背景, 将主动扰动抑制 [56] 用于控制通道扰动下的多维热方程, 将内模原理 [55] 用于有限维外系统生成扰动的边界输出调节. 本文不假设扰动由已知形式的有限维外系统生成, 而只设定其 $L^2(\omega)$ 范数有界; 相应地, 反馈采用多值符号项并形成微分包含. 摘录没有提供这些方法之间的定量控制代价比较, 也未声称在所有扰动类别上优于这些既有方法.

Navier–Stokes 稳定化已有较长研究脉络, 文中列举了内部和边界反馈以及随机系统方面的工作, 包括 [8, 9, 52, 38]、[13, 4]、[16, 14, 20, 3] 等. 其中若干研究关注局部稳定、相对于参考轨道的稳定、Riccati 或 Lyapunov 型反馈、临界函数空间中的适定性, 或边界控制. 作者特别指出, 既有三维结果在初值正则性和局部化控制方面各有不同限制; 本文目标则是在同一框架中处理未知局部匹配扰动与任意指定的指数衰减率. 文中还提及 $H^\infty$ 流体控制 [15] 的扰动衰减观点, 并将本文目标表述为对未知匹配扰动的拒斥加指定衰减率. 两类目标的性能指标不同, 摘录没有给出可直接比较的增益或代价界.

适用范围方面, Stokes 结果对任意维数 $n\ge2$ 、任意 $$H$$ 初值成立; 二维 Navier–Stokes 结果仅对 $$H$$ 中足够小的初值给出全局广义解与快速稳定; 三维结果要求初值属于 $$V$$ , $$V$$ 范数不超过给定的 $$R/2$$ , 且 $$H$$ 范数进一步小于依赖 $\lambda,D,R$ 的阈值, 从而得到全局强解及一致 $$V$$ 有界性. 所以三维结果不是消除 Navier–Stokes 全局正则性难题, 而是在特定小 $$H$$ 初值、受控 $$V$$ 初值范围内建立全局强解. 允许 $$R$$ 和 $$D$$ 任意给定, 但需要相应地缩小 $$H$$ 初值阈值; 这是相较于线性 Stokes 情形的实质代价.

快速衰减的定量代价也需一并考虑: 目标速率 $\lambda/4$ 可任意指定, 但稳定性常数为 $C_\lambda=C_1e^{C_1\sqrt{\lambda}}$ , 而低频反馈增益为 $\gamma_\lambda=\lambda C_1e^{C_1\sqrt{\lambda}}$ . 因而更快的衰减对应更大的反馈增益和更大的前置常数; 这并非对控制成本无代价的任意快稳定. 谱不等式常数只依赖 $(\Omega,\omega)$ , 但摘要和所给定理未给出控制能量或饱和约束分析.

中文标题: 拟线性 Hardy-Hénon 方程的尖峰自相似解及爆破后的延拓

摘要翻译: 我们考虑指数和维数满足

N\geq3, \quad σ\in(-2,\infty), \quad m>1, \quad p>p_S(σ)=\frac{m(N+2σ+2)}{N-2}.

的拟线性 Hardy-Hénon 方程

u_t=Δu^m+|x|^σu^p, \quad (x,t)\in\mathbb{R}^N\times(0,T),

。我们建立了指数的锐利范围 $$(m,p,σ)$$ ，从而保证存在具有后向自相似形式

arXiv ID: 2610.32707v1

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

摘要: We study the following prescribed scalar curvature problem

-Δu=K(y)u^{2^*-1},\ \ u>0 \,\ \ \text{in}\ \ \mathbb{R}^N, \tag*{$(*)$}

where $2^*=\frac{2N}{N-2}$ . For $N\geq4$ , we construct a positive $$K(y)$$ with maximum set $\{|y|=1\}$ and critical order of flatness $$m=N-2$$ at $$|y|=1$$ , such that problem $$(*)$$ admits only finitely many positive solutions in $D^{1,2}(\mathbb{R}^N)$ . To this end, we prove compactness of the solution set of $$(*)$$ by blow-up analysis and then apply Sard’s theorem to choose $$K$$ for which each solution is nondegenerate and hence isolated, so that the solution set of $$(*)$$ is finite. Our results give a negative answer to Wei–Yan’s conjecture in [J. Wei and S. Yan, JFA, 2010].

5. 确定性初始条件下的二维窄逃逸问题：平均逃逸时间与逃逸位置概率

The Two-Dimensional Narrow Escape Problem with Deterministic Initial Condition: Mean Escape Time and Escape location probabilities

作者: Loïs Delande, Tony Lelièvre, Boris Nectoux, Thomas Normand

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

发布日期: 2026-10-23

arXiv ID: 2610.32829v1

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

摘要: This paper presents a rigorous asymptotic analysis of the two-dimensional narrow escape problem for a Brownian particle initiated at a deterministic position x. While we address the mean escape time, the main novelty lies in our derivation of the asymptotic escape distribution. We derive explicit asymptotic formulas that characterize the dependence of the escape behavior on the relative configuration of the starting point and the window centers. These results remain valid as the starting position approaches the boundary windows and include refined error estimates. The formulas are obtained with two independent methods, each of them having their own interest. A central feature of this work is the rigorous treatment of the underlying mixed boundary value problems in the absence of a contact angle at the escape window boundaries. By leveraging recent results on the maximal regularity of the gradient of the solutions of such mixed boundary value problems, we ensure the analytical validity of our approach despite the inherent loss of regularity of these solutions at the escape window boundaries. Finally, we apply these general results to provide a full sharp asymptotic equivalents for regimes where window sizes scale at specific relative rates.

6. 界面问题中渗透率的作用：不可渗透屏障与自由扩散

The effect of the permeability in interface problems: impermeable barriers and free diffusion

作者: Pablo Álvarez-Caudevilla, Cristina Brändle, Mónica Molina-Becerra, Antonio Suárez

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32607v1

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

摘要: Focusing on reaction-diffusion models with an interface region dividing the domain, we study the stability of systems as the permeability parameter approaches zero or infinity under Kedem-Katchalsky boundary conditions imposed on the membrane or interface region. These boundary conditions establish that the flux across the interface manifold is proportional to the density jump across that region, with the constant of proportionality defined as the permeability parameter. We analyze both limiting behaviors: when permeability vanishes and when infinite flux is permitted. This analysis is carried out for both the well-known symmetric case, where permeability is identical in both directions across the membrane, and the asymmetric case, where it depends on the direction of the flux. For the asymmetric case, we obtain novel and surprising results that reveal a richer variety of behaviors in this class of interface problems and that differ from the scalar models.

7. 热含量系数的逆问题与连续性

Inverse problems and continuity of heat content coefficients

作者: Ryoki Endo, Katie Gittins, Jacob Walker

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32613v1

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

摘要: We consider the small-time asymptotic expansions for the Dirichlet and Open heat contents of convex planar domains. We prove that the small-time asymptotic expansion for the Open heat content can determine a triangle among triangles, a parallelogram among parallelograms, and an acute trapezium among acute trapezia of the same height. In addition, for a sequence of convex polygons converging to a convex set with smooth boundary in the sense of Hausdorff, we show that the coefficients of order $$t$$ in the small-time asymptotic expansions of the Dirichlet and the Open heat contents converge to those of the limiting set. While, for a sequence of convex sets with smooth boundary converging to a convex polygon in the sense of Hausdorff, we show that the coefficients of order $$t$$ do not converge to that of the limiting polygon.

8. 从动理学激活动力学到 Holling 型反应扩散捕食者-猎物系统

From kinetic activation dynamics to Holling type reaction-diffusion predator-prey systems

作者: Paolo Buttà

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32603v1

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

摘要: We consider a kinetic model of predator-prey dynamics in which predation is characterized by a sequence of behavioral transitions. Prey and predators are described by distribution functions in phase space, with the predator population further structured into internal states corresponding to inactive, active, and saturated individuals. We show that different macroscopic functional responses emerge from distinct mechanisms governing predator activation. More precisely, activation driven uniquely by pairwise encounters leads to a saturating Holling type II response, whereas activation additionally modulated by local prey density produces a generalized sigmoidal Holling type III response. This is achieved by means of a diffusive scaling and a formal hydrodynamic limit, from which we derive a closed reaction-diffusion system of Rosenzweig-MacArthur type for the prey and total predator spatial densities.

9. 取值于网络的能量与几何拟凸性

Energies with values in networks and geometric quasiconvexity

作者: Andrea Braides, Jan Kristensen, Sergio Scalabrino

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32444v1

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

摘要: We study the asymptotic behaviour, as $\varepsilon\to 0^+$ , of the $$p$$ -Dirichlet energy of maps $u\colonΩ\subset\mathbb{R}^2 \to \mathbb{R}^2$ constrained to take values in the rescaled periodic network $\varepsilon\mathcal Z$ , where $\mathcal Z$ is a standard square network in $\mathbb{R}^2$ . The oscillations are carried by the constraint in the target, not by the energy density, and their effect depends on $$p$$ . For $$1<p<2$$ the constraint disappears in the limit: the $$Γ$$ -limit is an integral functional on the whole of $\rm{W}^{1,p}(Ω, \mathbb{R}^2 )$ , whose density is given by an asymptotic cell formula on tilted cells. This density is equal to $$(|a_1|+|a_2|)^p$$ on rank-one matrices $a\otimesν$ , $$|ν|=1$$ , and exceeds its convex envelope by a quantity comparable to $|\det A|\,|A|^{p-2}$ , with constants of order $(2-p)^{-1}$ . Consequently, although the original energy density is convex, the limit density is quasiconvex but neither convex nor polyconvex: quasiconvexity is generated by the geometry of the constraint. For $p\ge2$ the constraint survives in the limit as $\det\nabla u=0$ , and the $$Γ$$ -limit is the degenerate energy $\int_Ω(|\nabla u_1|+|\nabla u_2|)^p\, \rm{d} x$ , that is, the $$p$$ -energy of $$u$$ considered as a map into the normed plane $(\mathbb{R}^2,|\cdot|_1)$ , the metric limit of the networks.

10. 弹性八字结并非平面

The elastic figure-eight knot is not planar

作者: Tatsuya Miura

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

发布日期: 2026-10-23

arXiv ID: 2610.32320v1

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

摘要: We prove that every elastic knot in any non-torus twist knot class is not planar. In particular, this disproves a numerically motivated planarity conjecture for elastic figure-eight knots.

11. 随机抛物型方程的几何性质与有界域上的 Prékopa-Leindler 不等式

Geometry property of stochastic parabolic equations and Prékopa-Leindler inequality on bounded domains

作者: Guangying Lv

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32288v1

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

摘要: We study the preservation of spatial concavity and log-concavity for stochastic parabolic equations posed on bounded convex domains. For the linear noise, by using Minkowski convolutions, we establish the convexity preservation of solution almost surely. For the nonlinear noise, by using Itô formula, the impact of noise is given. Moreover, we establish the relationship between the weighted (stochastic) Prékopa-Leindler inequality and properties of (stochastic) parabolic equation on bounded domain, which is new even for deterministic case.

12. 圆球面上共形爱因斯坦–狄拉克系统的耦合扰动

Coupling perturbations of the conformal Einstein-Dirac system on round spheres

作者: Ke Xu, Tian Xu, Minbo Yang

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32172v1

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

摘要: We study the conformal Einstein-Dirac system on the round sphere $$S^m$$ , $m\geq3$ , with coupling coefficient $1+\eps k$ , where $k\in C^2(S^m,\R)$ and $\eps$ is small. Morse and degree conditions yield solutions with positive scalar component, and a generic nondegeneracy condition gives multiplicity modulo spinorial symmetries. For a suitable fixed coefficient, we obtain a noncompact solution sequence with infinitely many distinct energy values. A Kazdan-Warner identity gives nonexistence for a smooth Morse potential and its small $$C^2$$ -perturbations, even with positive coupling. The existence proofs use a local reduction in the interior of the positive cone, accommodating the negative powers occurring when $m\geq5$ .

13. 环域上多重情形设定下 $$p(x)$$ -Laplacian 抛物型方程的适定性

Well-posedness for $$p(x)$$ -Laplacian parabolic equations under a multiple-regime setting on an annulus

作者: Yitian Wang, Runzhang Xu, Chao Yang

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32083v1

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

摘要: This paper proposes an energy-variational framework for the well-posedness of the initial boundary value problem of the parabolic equation with variable exponents on an annulus, under the multiple-regime setting consisting of the subcritical, critical and supercritical regimes of the variable exponents. The multiple-regime setting leads to the appearance of the critical and supercritical regions of positive measure, which cause the failure of the classical Sobolev embedding and the lack of compactness. To address the challenges posed by multiple regimes to well-posedness within the energy framework, we recover the compact Sobolev embedding for a class of radially symmetric functions by the Strauss inequality and obtain the local existence of the solution via the combination of the subdifferential technique and the Principle of Symmetric Criticality for any given radially symmetric initial data. Additionally, the uniqueness result is also given with the recovered Sobolev embedding. Moreover, examining the impact on well-posedness of solutions caused by the variable exponents crossing three sensitive regions, we establish a variational structure of the potential well theory for the evolution equation under the multiple-regime setting on the annulus. For subcritical and critical initial energy levels, we consider the non-positive and positive initial energy cases separately. We obtain the finite-time blow-up of the solution with non-positive initial energy and the threshold condition for global existence and finite-time blow-up of the solution with positive initial energy. For supercritical initial energy level, several sufficient criteria for the global existence and finite-time blow-up of the solution are given. Finally, the above results, including global existence and finite-time blow-up, are generalized to a more general locally symmetric domain containing an annulus.

14. 从单个终态全局恢复标量扩散系数

Global recovery of a scalar diffusivity from a single final state

作者: Ching-Lung Lin, Yi-Hsuan Lin

类别: math.AP

发布日期: 2026-10-23

arXiv ID: 2610.32063v1

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

摘要: We prove global uniqueness and Lipschitz stability for a scalar diffusivity in the heat equation from a single final state at any prescribed positive time. On each prescribed uniformly elliptic $C^{1,α}$ -bounded class with known coefficient boundary values, we choose a constant initial datum and a spatially constant exponential boundary input that work for the entire class. The result holds on smooth bounded connected domains in every dimension $n\ge2$ and gives an affirmative answer to \cite[Problem~8.2]{Alessandrini} within these classes. The stability estimate controls the $$L^2$$ diffusivity difference by the $$H^2$$ final-state difference, and this data norm is optimal within the Sobolev scale for the chosen experiment.

15. 满足容量密度条件的区域上的加权与非局部泊松方程

Weighted and nonlocal Poisson equations on domains with capacity density conditions

作者: Takanobu Hara

类别: math.AP, math.FA

发布日期: 2026-10-23

arXiv ID: 2610.32070v1

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

摘要: We establish necessary and sufficient conditions for global Hölder solvability of weighted Poisson equations and related fractional problems. On domains satisfying a capacity density condition, a nonnegative measure with bounded support depth admits a bounded globally Hölder continuous zero-boundary solution if and only if it satisfies an interior Morrey-type growth condition. Conversely, uniform boundary decay estimates for solutions with localized nonnegative sources characterize the capacity density condition. Analogous characterizations are established for the restricted and spectral fractional Laplacians. Our approach is based on weighted strong barriers and also yields two-weight Hardy inequalities, compact trace embeddings, and global Hölder solvability for nonlinear Dirichlet-to-Neumann operators.

16. 利比里亚社区卫生项目中药品最后一公里配送的动态库存配给模型

A Dynamic Inventory Rationing Model for Last-Mile Delivery of Medication in Liberia’s Community Health Program

作者: Baris Ata, Yufan Chen, Robert Montgomery

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

发布日期: 2026-10-23

arXiv ID: 2610.32051v1

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

摘要: We study the dynamic allocation of scarce medications in Liberia’s community health program. A central warehouse supplies a health facility and the Community Health Assistants (CHAs) it supervises. When a CHA runs out of medication, patients are referred to the facility, incurring travel costs if they complete the referral and remain untreated otherwise. Allocating inventory to CHAs allows more patients to receive treatment without traveling, while retaining inventory centrally preserves flexibility to respond to uncertain demand. The allocation problem must balance these considerations while accounting for fairness across locations. We formulate this problem as a finite-horizon stochastic control problem whose objective combines the costs of referrals and untreated patients with a demand-weighted fairness utility. Because standard dynamic programming methods are computationally intractable in this setting, we develop a continuous-time approximation that involves controlling the drift rate of a reflected Brownian motion. We then adapt the deep splitting method, using deep neural networks to approximate the value function and its gradient and compute a dynamic allocation policy. We calibrate the model using data on malaria cases among children in the Fangonda region of Lofa County, Liberia, and evaluate the proposed policy in the original discrete-time model. In the baseline setting, our policy reduces operating costs by about 6% relative to the strongest benchmark, the fairness-only policy, while retaining approximately 98% of its fairness utility. Additional experiments show that higher referral costs shift inventory toward CHAs and reduce local unmet demand, while greater demand variability increases our policy’s performance over the fairness-only benchmark. Our model and computational approach provide a framework for studying similar last-mile allocation problems.

17. 三维中一个非拟凸的框架无关秩一凸函数

A frame-indifferent rank-one convex function in three dimensions which is not quasiconvex

作者: Sergio Conti, Georg Dolzmann, Patrick Dondl

类别: math.AP

发布日期: 2026-10-22

arXiv ID: 2610.31886v1

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

摘要: We construct a function $W:\mathbb{R}^{3\times3}\to\mathbb{R}$ which is frame-indifferent, $$W(QF)=W(F)$$ for all $Q\in\mathrm{SO}(3)$ , rank-one convex, but not quasiconvex. In a neighbourhood of $\mathrm{SO}(3)$ the function is an explicit polynomial of degree four in the strain $F^TF-\mathrm{Id}$ ; the failure of quasiconvexity occurs at $F=\mathrm{Id}$ through a third-order mechanism of Šverák type, with test functions whose gradients are arbitrarily close to $\mathrm{Id}$ . The polynomial is obtained from a geometrically linear example, a rank-one convex, non-quasiconvex quartic polynomial of the symmetric part of the gradient, by evaluating it at the strain $\frac12(F^TF-\mathrm{Id})$ , adding a uniaxial prestress and a cubic correction, and enlarging the quartic term. The local function is extended to all of $\mathbb{R}^{3\times3}$ by adding a frame-indifferent polyconvex function which vanishes near $\mathrm{SO}(3)$ . This work was carried out with the assistance of the AI model Claude Opus 5.5; the proofs of the main results have been formally verified in Lean 4.

18. 强耦合抛物系统在可测集上的可观测性

Observability from measurable sets for strongly coupled parabolic systems

作者: Xinyi Chen, Kévin Le Balc’h

类别: math.OC, math.AP

发布日期: 2026-10-22

arXiv ID: 2610.31876v1

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

摘要: We characterize observability of strongly coupled parabolic systems from arbitrary space-time measurable sets of positive measure. For constant coercive diffusion matrices and arbitrary constant zeroth-order coupling and observation matrices, observability holds if and only if the generalized Kalman rank condition is satisfied at every Dirichlet eigenvalue. The diffusion matrix need not be diagonalizable, and the number of components is arbitrary. The observation is measured in $$L^1$$ , yielding null controllability with bounded controls by duality. The proof combines a Remez-type inequality for exponential polynomials with a quantitative analysis of the mode-dependent Kalman condition. The resulting finite-dimensional observability estimate is uniform over the spatial modes up to a polynomial frequency loss. A spectral inequality and the Lebeau–Robbiano strategy then give observability of the full system.

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

1. 欧拉方程中从熵波到涡片#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. gSQG 补丁方程的自相似解#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 可压缩 Navier-Stokes-Poisson 系统在 Navier 滑移边界条件下的一致正则性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 广义相对论中具有自由边界的正压液态恒星的线性稳定性与转折点原理#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 具有超线性 BD 黏性和远场真空的三维大振幅球对称退化可压缩 Navier-Stokes 方程的整体正则性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 具有局部化扰动的 Stokes 与 Navier-Stokes 系统的快速稳定化#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 含奇异外势的三维 Hartree 方程的修正散射#

摘要翻译#

2. A2 条件与 Helson–Szego 条件等价性的直接证明#

摘要翻译#

3. 由薄桁架穿孔区域中的 Stokes 方程均匀化#

摘要翻译#

4. 二维辐射输运方程的角点正则性与渐近行为#

摘要翻译#

5. p < q 时加权 Triebel–Lizorkin Hardy 不等式的锐性转变#

摘要翻译#

6. $\mathbb{S}^{3}$ 上三次共形波动方程的新时间周期解：I 球对称性#

摘要翻译#

7. 慢扩散型双重非线性方程内在抛物型 Harnack 不等式的刻画#

摘要翻译#

8. Riesz 变换的有界性蕴含图的抛物一致可整性#

摘要翻译#

9. 双曲空间上的双线性 Strichartz 估计与三次 NLS#

摘要翻译#

10. 拟线性 Hardy-Hénon 方程的尖峰自相似解及爆破后的延拓#

摘要翻译#

11. 具有调和势和组合幂型非线性的半线性薛定谔方程#

摘要翻译#

12. 全空间中截断软势玻尔兹曼方程的最优时间衰减#

摘要翻译#

13. 非线性项支撑于超曲面上的三次非线性薛定谔方程的适定性#

摘要翻译#

14. 带小孔流形上的 Gagliardo–Nirenberg 极值函数#

摘要翻译#

15. 集体动力学的最优图子函数控制#

摘要翻译#

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

1. 周期项链图上的平稳双稳态前沿：存在性、稳定性与多重性#

2. 热点：从光滑到非光滑的单连通平面区域#

3. 正整数格上的 Hardy 不等式：维数转变#

4. 临界平坦阶下指定标量曲率问题的有限多个解#

5. 确定性初始条件下的二维窄逃逸问题：平均逃逸时间与逃逸位置概率#

6. 界面问题中渗透率的作用：不可渗透屏障与自由扩散#

7. 热含量系数的逆问题与连续性#

8. 从动理学激活动力学到 Holling 型反应扩散捕食者-猎物系统#

9. 取值于网络的能量与几何拟凸性#

10. 弹性八字结并非平面#

11. 随机抛物型方程的几何性质与有界域上的 Prékopa-Leindler 不等式#

12. 圆球面上共形爱因斯坦–狄拉克系统的耦合扰动#

13. 环域上多重情形设定下 $p(x)$-Laplacian 抛物型方程的适定性#

14. 从单个终态全局恢复标量扩散系数#

15. 满足容量密度条件的区域上的加权与非局部泊松方程#

16. 利比里亚社区卫生项目中药品最后一公里配送的动态库存配给模型#

17. 三维中一个非拟凸的框架无关秩一凸函数#

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

1. 欧拉方程中从熵波到涡片

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. gSQG 补丁方程的自相似解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 可压缩 Navier-Stokes-Poisson 系统在 Navier 滑移边界条件下的一致正则性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 广义相对论中具有自由边界的正压液态恒星的线性稳定性与转折点原理

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 具有超线性 BD 黏性和远场真空的三维大振幅球对称退化可压缩 Navier-Stokes 方程的整体正则性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 具有局部化扰动的 Stokes 与 Navier-Stokes 系统的快速稳定化

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 含奇异外势的三维 Hartree 方程的修正散射

摘要翻译

2. A2 条件与 Helson–Szego 条件等价性的直接证明

摘要翻译

3. 由薄桁架穿孔区域中的 Stokes 方程均匀化

摘要翻译

4. 二维辐射输运方程的角点正则性与渐近行为

摘要翻译

5. p < q 时加权 Triebel–Lizorkin Hardy 不等式的锐性转变

摘要翻译

6. $\mathbb{S}^{3}$ 上三次共形波动方程的新时间周期解：I 球对称性

摘要翻译

7. 慢扩散型双重非线性方程内在抛物型 Harnack 不等式的刻画

摘要翻译

8. Riesz 变换的有界性蕴含图的抛物一致可整性

摘要翻译

9. 双曲空间上的双线性 Strichartz 估计与三次 NLS

摘要翻译

10. 拟线性 Hardy-Hénon 方程的尖峰自相似解及爆破后的延拓

摘要翻译

11. 具有调和势和组合幂型非线性的半线性薛定谔方程

摘要翻译

12. 全空间中截断软势玻尔兹曼方程的最优时间衰减

摘要翻译

13. 非线性项支撑于超曲面上的三次非线性薛定谔方程的适定性

摘要翻译

14. 带小孔流形上的 Gagliardo–Nirenberg 极值函数

摘要翻译

15. 集体动力学的最优图子函数控制

摘要翻译

不相关论文（基本信息）

1. 周期项链图上的平稳双稳态前沿：存在性、稳定性与多重性

2. 热点：从光滑到非光滑的单连通平面区域

3. 正整数格上的 Hardy 不等式：维数转变

4. 临界平坦阶下指定标量曲率问题的有限多个解

5. 确定性初始条件下的二维窄逃逸问题：平均逃逸时间与逃逸位置概率

6. 界面问题中渗透率的作用：不可渗透屏障与自由扩散

7. 热含量系数的逆问题与连续性

8. 从动理学激活动力学到 Holling 型反应扩散捕食者-猎物系统

9. 取值于网络的能量与几何拟凸性

10. 弹性八字结并非平面

11. 随机抛物型方程的几何性质与有界域上的 Prékopa-Leindler 不等式

12. 圆球面上共形爱因斯坦–狄拉克系统的耦合扰动

13. 环域上多重情形设定下 $$p(x)$$ -Laplacian 抛物型方程的适定性

14. 从单个终态全局恢复标量扩散系数

15. 满足容量密度条件的区域上的加权与非局部泊松方程

16. 利比里亚社区卫生项目中药品最后一公里配送的动态库存配给模型

17. 三维中一个非拟凸的框架无关秩一凸函数