生成时间: 2026年09月11日 07:33:02

运行状态: 完成

计划处理论文数量: 43

已完成论文数量: 43

未完成论文数量: 0

重点关注论文数量: 10

了解领域论文数量: 14

不相关论文数量: 19

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

1. 可压缩等熵 Navier-Stokes 方程的点态衰减极限

Limiting Pointwise Decay for the compressible isentropic Navier-Stokes equations

作者: Rui Li, Xu Song, Guanghui Wang, Lingda Xu

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11707v1

摘要翻译: 我们研究一维可压缩等熵 Navier-Stokes 方程中常数态的小局部扰动的长时间点态行为。在减去两个 Burgers 扩散波、所有高阶扩散波的收敛和以及跨族粘性修正之后，我们对精确的物理余项证明了一个保锥的点态估计，特别是

Comment: 无

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

1. 研究对象和背景

本文研究一维可压缩等熵 Navier-Stokes 方程（也称为粘性 p-系统）在常状态附近的小幅局部扰动解的长时间渐近行为。具体而言，考虑如下方程组：

\begin{cases} v_t - u_x = 0, \\ u_t + p(v)_x = \nu \left(\frac{u_x}{v}\right)_x, \end{cases}

其中 $$v$$ 为比容， $$u$$ 为速度， $$x$$ 为拉格朗日质量坐标。假设压强函数 $$p$$ 在 $$v=1$$ 附近满足 $$p'(1) < 0$$ 且 $p''(1) \neq 0$ ，即具有真实的非线性特征。该模型是部分耗散双曲-抛物守恒律的原型，仅速度方程直接具有抛物性，而比容的耗散通过双曲耦合间接传递。

该问题的核心难点在于，线性化系统具有两个声速，导致局部扰动分裂为沿相反特征方向传播的两个扩散波包。由于 $p''(1) \neq 0$ ，每个真实非线性族的二次自相互作用在扩散尺度上存活，将线性的热波转变为被对流输运的 Burgers 扩散波。这些 Burgers 波不仅描述了低频部分的非线性有效动力学，而且其质量守恒性决定了波的振幅。然而，全局 $$L^p$$ 估计往往只能捕捉主导的质量携带波，却掩盖了波之间的代数衰减尾部信号。因此，必须区分接收特征锥、相反锥、两锥之间区域以及外部区域的空间几何结构，才能准确描述系统的长时间行为。

已有文献中，Chern 和 Liu 证明了向热波和 Burgers 扩散波叠加的收敛性；Liu 和 Zeng 建立了相应的 $$L^1$$ 渐近理论；Liu 和 Yu 发展了一维气体动力学的详细 Green 函数理论。Van Baalen, Popović 和 Wayne 以及 Koike 进一步发展了高阶渐近修正理论，构建了任意高阶的扩散波层级。Koike 的结果表明，随着保留的高阶波数量 $$n$$ 的增加，幂律指数 $eta_n = 2 - 2^{-(n+1)}$ 趋近于极限值 2，对应全局 $L^\\infty$ 估计中的 $(1+t)^{-1}$ 衰减率。然而，Koike 的估计是针对固定有限阶展开的，常数依赖于阶数，无法直接取极限 $n \to \infty$ 。Koike 在其论文的备注中提出了一个关于极限逐点权重的猜想，但并未给出证明。本文旨在填补这一空白，即在减去整个收敛的高阶波层级和跨族粘性修正后，直接估计物理余项的逐点行为，从而验证或推进这一极限衰减率的猜想。

2. 主要定理或主要结果

本文的主要结果是定理 2.2，它给出了物理余项 $$R_i(x,t)$$ 的锥保持逐点估计。具体地，设 $N \ge 5$ ，初始数据满足小性条件 $E_N \le \epsilon_0$ ，则对于 $$i=1,2$$ ，有如下估计：

|R_i(x, t)| \le C E_N \log(2 + t) \Psi_i(x, t),

其中 $$E_N$$ 度量初始数据的大小， $ig|\Psi_i(x,t)$ 是一个精细定义的锥分辨权重函数。特别地，由此可得全局 $L^\infty$ 范数的衰减估计：

\|R_i(t)\|_{L^\infty} \le C E_N (1+t)^{-1} \log(2+t).

此外，还得到了导数的逐点估计：

\|R_{i,x}(t)\|_{L^\infty} \le C E_N (1+t)^{-3/2} \log(2+t),

以及余项在 Sobolev 空间 $$H^N$$ 中的均匀有界性：

\sup_{t \ge 0} \|(R_1, R_2)(t)\|_{H^N} \le C E_N.

这里的权重函数 $\Psi_i(x,t)$ 保留了两个对流空间变量 $x - \lambda_i(1+t)$ 和 $x - \lambda_{3-i}(1+t)$ ，分别对应两个特征族的传播方向。这种“锥保持”性质意味着估计能够区分接收锥、相反锥、两锥之间区域和外部区域，从而精确捕捉不同区域的信号强度。例如，在接收锥 $|x - \lambda_i(1+t)| = O(\sqrt{1+t})$ 附近，权重表现为高斯型的快速衰减；而在相反锥或两锥之间，权重则表现为代数衰减。这种几何结构的保留是全局 $L^\infty$ 估计所不具备的。值得注意的是，虽然结果达到了极限代数衰减率 $(1+t)^{-1}$ ，但存在一个对数损失 $\log(2+t)$ 。这与 Koike 提出的猜想相比稍弱，因为猜想可能期望更优的对数权重或无对数损失，但本文的结果首次严格证明了在减去整个收敛层级后余项的极限衰减行为。

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

证明的核心障碍在于如何处理同族的一阶反馈项。在引入空间原函数 $W_i(x,t) = \int_{-\infty}^x R_i(y,t) dy$ 后，其演化方程中包含关键的反馈项 $\theta_i W_{i,x}$ ，其中 $\theta_i$ 是第 $$i$$ 个族的 Burgers 扩散波。由于 $\theta_i$ 以非可积速率 $O((1+t)^{-1/2})$ 衰减，直接扰动处理该项会导致对数发散，从而破坏极限估计。作者的真正创新点是应用一族 Cole-Hopf 变换来消除这一临界反馈。具体地，定义 Cole-Hopf 权重 $\phi_i$ 并令 $Z_i = \phi_i W_i$ 。利用 Cole-Hopf 恒等式，可以精确消除 $\theta_i W_{i,x}$ 项，而不是对其进行估计。这种变换不仅是表示 Burgers 轮廓的工具，更是一种族级共轭，消除了耦合余项方程中的主要变系数障碍。

然而，仅靠 Cole-Hopf 消去并不足以达到极限估计。Duhamel 表示法中 $$Z_i$$ 的主部高斯输出系数在时间上仅是边界可积的。因此，作者引入了第二个结构性步骤：时间依赖分解

Z_i(x, t) = a_i(t) \Gamma_i(x, 1+t) + Z^{rem}_i(x, t),

其中 $\Gamma_i$ 是第 $$i$$ 个族的对流热核。系数 $$a_i(t)$$ 并非简单的空间质量，而是直接从 Duhamel 源中提取所有主部同族高斯系数以及由普通非线性项、异常反馈和显式轮廓源生成的锥-时间局部化相反族高斯系数。通过这一定义， $Z^{rem}_i$ 中不再包含潜在的对数主部高斯项。闭合非线性 bootstrap 后，模态系数满足 $|a_i(t)| \le C E_N \log(2+t)$ 和 $|a'_i(t)| \le C E_N (1+t)^{-1}$ ，而余项 $Z^{rem}_i$ 满足无对数损失的估计 $\|\partial_x^k Z^{rem}_i(t)\|_{L^\infty} \le C E_N (1+t)^{-(k+1)/2}$ 。对数损失被限制在提取的高斯模态中，而通过恢复公式 $R_i = \phi_i^{-1} (Z_{i,x} + \frac{\theta_i}{2d} Z_i)$ ，导数结构将这些界转化为物理余项的极限衰减。

在线性传播方面，作者构造了一个适应于两个特征族的近似 Green 函数。通过将族级 Cole-Hopf 权重插入常状态 Green 矩阵，得到共轭系统的近似 Green 函数。该构造同时适应两个传播方向和部分耗散粘性耦合，其与共轭算子的交换子具有端点论证所需的散度和零质量结构。最后，结合 Kawashima 补偿子和高斯加权能量估计，闭合了均匀高阶界。

4. 与之前工作的比较

与既有工作相比，本文的主要推进在于从“固定有限阶展开”到“整个收敛层级的极限”的转变。Koike [6] 建立了递归高阶扩散波层级，并证明了在每个固定有限阶修正后的逐点余项估计。这些估计捕捉了主导 Burgers 波指数小的区域内的幂律渐近，但其常数依赖于保留的阶数 $$n$$ ，因此不能直接取 $n \to \infty$ 的极限。Van Baalen, Popović 和 Wayne [16] 在 $$L^2$$ 框架下识别了对数加权的无限阶极限，但并未提供相应的逐点估计。Koike 在其论文的备注 2.5 中提出了一个关于极限逐点权重的猜想，但未予证明。本文首次严格证明了在减去整个收敛层级和跨族粘性修正后，物理余项的锥保持逐点估计，达到了极限代数衰减率 $(1+t)^{-1}$ ，尽管存在一个对数损失 $\log(2+t)$ 。

与 Liu 和 Yu [10] 的 Green 函数理论相比，本文利用了他们的正则-奇异 Green 函数估计作为分析基础，但进一步引入了 Cole-Hopf 共轭和高斯模态提取技术来处理非线性反馈。Liu 和 Yu 的工作侧重于线性传播和奇异性传播，而本文的重点是非线性余项的精细渐近结构。与 Deng 和 Yu [3] 的正则-奇异 Green 函数构造相比，本文的 Green 函数是专门为共轭系统设计的，以适应 Cole-Hopf 变换后的方程结构。

与 Chern 和 Liu [1, 2] 的早期工作相比，本文不仅关注主导波的收敛性，还深入研究了高阶修正和跨族交互作用对衰减率的影响。Chern 和 Liu 的修正揭示了累积跨族交互作用可以改变甚至决定最优衰减率，但本文通过精确的余项估计，量化了这一效应在极限情况下的表现。与 Xin 和 Zeng [18] 的一般扰动框架相比，本文针对一维可压缩流的具体结构进行了优化，利用了 Burgers 波和 Cole-Hopf 变换的特殊性质，从而获得了更精细的逐点估计。

总体而言，本文的贡献在于提供了一种新的分析框架，结合了 Cole-Hopf 共轭、高斯模态提取和自适应 Green 函数，从而能够处理极限情形下的临界反馈问题。这一方法不仅适用于当前的 p-系统，也可能推广到其他部分耗散双曲-抛物系统。然而，代价是证明的复杂性显著增加，且结果中存在一个对数损失，这可能源于模态系数的边界可积性。未来工作可能致力于消除这一对数损失，或者将结果推广到更高维度或其他非线性结构。

2. 二维不可压欧拉方程柯西问题已报道多解的数值研究

Numerical Study of Reported Multiple Solutions to a Cauchy Problem for the 2d Incompressible Euler Equations

作者: Allen Tesdall, Richard Sanders

类别: math.AP, math.NA

发布日期: 2026-09-10

arXiv ID: 2609.11694v1

摘要翻译: Bressan 和 Shen (2021) 与 Bressan 和 Murray (2020) 共同提出了一个位于 $\mathbb{R}^2$ 上的柯西问题的简单例子，该例子由两个关键参数刻画，旨在产生不可压欧拉方程的非唯一弱解。他们建议，对其示例初值施加某种极限过程并应用于两个特定的近似，可得到两个不同的解。在此，我们提供了一个经过精心设计的数值程序所提供的证据，以支持这一引人入胜的论断是否成立，或在何时成立。本研究中遇到的一些数值挑战包括：(i) 所提出的初值在原点的邻域内无界；(ii) 所提出的问题定义在整个 $\mathbb{R}^2$ 上；(iii) 真实解中存在的若干基本对称性必须被捕捉到；(iv) 所提出的示例依赖于两个可能密切相关的关键参数。我们证明，为了捕捉精确自由空间问题的核心对称性，在数学上需要一个极大的计算区域。相应地，我们实施了一种嵌套网格策略，专门针对其在当前问题中的特殊适用性和有效性进行设计。我们的策略还允许我们在解精度最为关键的区域实现极高的网格细化。在我们的嵌套网格框架内，我们设计了一个高阶、稳健且快速的底层离散化方案。我们考虑了广泛的可允许参数对集合，但在其中仅发现一对可能满足所建议的论断。

Comment: 37 pages, 20 figures

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

1. 研究对象和背景

本文的研究对象是二维不可压欧拉方程的柯西问题。该方程组描述了无粘、不可压缩流体在 $\mathbb{R}^2$ 中的运动，其数学形式为：

u_t + (u \cdot \nabla)u + \nabla p = 0, \quad \text{div } u = 0.

其中 $$u(t, x)$$ 为流体速度场， $$p(t, x)$$ 为压力标量函数， $$x=(x,y)$$ 。

为了更清晰地分析涡旋结构，通常将上述方程改写为流函数 $\psi$ 和涡度 $\omega$ 的形式：

\omega_t + u \cdot \nabla \omega = 0, \quad \nabla^2 \psi = \omega,

其中速度与流函数的关系定义为 $u = \nabla^\perp \psi \equiv (-\psi_y, \psi_x)$ ，且涡度 $\omega = \text{curl } u = \nabla^2 \psi$ 。初始条件由给定的涡度分布 $\omega_0(x) = \omega(0, x)$ 确定。

Bressan 和 Shen (2021) 以及 Bressan 和 Murray (2020) 提出了一个旨在展示非唯一弱解的简单示例。他们考虑的初始涡度数据在极坐标 $(r, \theta)$ 下具有特定的对称性和奇异性，形式为：

\omega_0(x) = r^{-\alpha}\phi(\theta),

其中参数 $\alpha \in (0, 2)$ 描述了径向依赖性，角向部分 $\phi(\theta)$ 满足以下性质：

\phi(\theta + \pi) = \phi(\theta), \quad \phi(\theta) = 0 \text{ if } \theta \in [\frac{\pi}{4}, \pi].

值得注意的是，在 $\phi \neq 0$ 的射线上，当 $|x| \to 0$ 时，初始涡度 $\omega_0(x) \to \infty$ 。这种在原点附近的无界性是问题的核心特征之一。

为了处理这一奇异性，作者引入了两个有界的正则化初始数据族 $\omega_{0,1,\epsilon}$ 和 $\omega_{0,2,\epsilon}$ ：

\omega_{0,1,\epsilon}(x) = \begin{cases} \omega_0(x) & \text{if } |x| > \epsilon \\ \epsilon^{-\alpha} & \text{if } |x| \le \epsilon \end{cases}, \quad \omega_{0,2,\epsilon}(x) = \begin{cases} \omega_0(x) & \text{if } |x| > \epsilon \\ 0 & \text{if } |x| \le \epsilon. \end{cases}

这两个序列在 $L^p_{loc}$ 范数下（对于任意 $1 \le p < 2/\alpha$ ）收敛于原始奇异数据 $\omega_0$ 。区别在于， $\omega_{0,1,\epsilon}$ 在一个小圆盘内保持恒定高涡度，而 $\omega_{0,2,\epsilon}$ 在同一区域内涡度为零。

Bressan 等人推测，随着 $\epsilon \to 0$ ，分别由这两个正则化数据产生的解会收敛到两个不同的自相似极限解：一个包含单个螺旋中心（位于原点），另一个包含两个关于原点对称的螺旋中心。如果这一推测成立，它将证明该柯西问题存在“不可治愈”的病态性（即非唯一性）。然而，先前的分析依赖于复杂的域分解方法和小范围数值模拟，缺乏对大范围自由空间问题中对称性和尺度行为的严格验证。本文旨在通过高精度的数值实验来检验这一非唯一性声明的有效性，填补从有限域近似到无限域极限过程中的理论空白。

2. 主要定理或主要结果

本文的主要结论基于对一系列允许参数对的广泛数值探索。尽管考虑了多种可能的参数组合，但研究发现只有一组参数可能满足 Bressan 等人提出的非唯一性声明。

具体而言，作者定义了如下两个问题：

问题 1：系统 (2) 配合初始数据 $\omega_{0,1,\epsilon}$ 。
问题 2：系统 (2) 配合初始数据 $\omega_{0,2,\epsilon}$ 。

主要结果包括：

对称性验证：作者证明了要捕捉精确自由空间问题的核心对称性，必须使用极其大的计算域。在较小的域上，螺旋中心的时间轨迹无法满足关键的对称性约束。通过使用嵌套网格策略和大尺度域，作者展示了计算解确实高精度地满足了这些对称性。
自相似性检查：根据缩放对称性，若存在自相似极限 $\Omega(y)$ ，则对于 $$t > 0$$ ，应有 $t \omega_{\epsilon}(t, x) \sim \Omega(x/t^{1/\alpha})$ 。特别地，如果极限解包含非零位置的螺旋中心 $y^* \neq 0$ ，则峰值涡度位置 $x_{\epsilon}(t)$ 应近似满足：
$x_{\epsilon}(t) \sim t^{1/\alpha} y^*.$
这意味着螺旋中心应沿从原点发出的直线运动，且具有可观测的速度。作者通过追踪数值解中峰值涡度的位置随时间和 $\epsilon$ 的变化，检验了这一行为是否成立。
参数敏感性：在测试了广泛的参数范围后，只有特定的参数配置显示出与理论预测一致的行为。这表明非唯一性并非普遍现象，而是依赖于初始数据的具体细节（如 $\alpha$ 的值和 $\phi(\theta)$ 的形状）。
数值证据的限制：虽然数值结果支持在某些特定条件下存在两个不同解的可能性，但作者强调这是“证据”而非严格证明。由于计算资源的限制，无法完全消除边界效应或达到真正的 $\epsilon \to 0$ 极限，因此结论带有条件性。

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

本研究的核心挑战在于如何在有限的计算资源下准确模拟定义在整个 $\mathbb{R}^2$ 上的奇异初值问题，同时保持解的关键对称性和尺度特性。作者开发了一套专门针对此问题的数值框架，主要包括以下几个方面：

嵌套网格策略 (Nested-Grid Strategy)：为了平衡大尺度域的需求和局部高精度解析的需要，作者设计了一种嵌套网格方法。该方法允许在靠近原点（解最复杂、精度要求最高）的区域进行极端网格细化，同时在远离原点的区域使用较粗的网格以覆盖更大的物理域。这种策略不仅提高了计算效率，还确保了远场边界效应对核心区域的影响最小化。
高阶离散方案：
- 空间离散：采用有限体积法结合 minmod 斜率限制器来实现二阶精度的空间离散。涡度守恒律被离散为： $\frac{\partial \omega}{\partial t} + \frac{\partial}{\partial x}(u \omega) + \frac{\partial}{\partial y}(v \omega) = 0.$ 数值通量使用迎风格式构造，并结合 minmod 限制器防止振荡： $s_{L|R}^{i,j} = \omega_{i,j} \mp \frac{1}{2} \text{minmod}(\omega_{i+1,j} - \omega_{i,j}, \omega_{i,j} - \omega_{i-1,j}).$
- 时间推进：使用 Heun 的两步法进行时间积分，这是一种显式二阶 Runge-Kutta 方法： $\tilde{\omega} = \omega^n - \Delta t F(u^n, \omega^n),$ $\omega^{n+1} = \omega^n - \frac{1}{2} \Delta t (F(u^n, \omega^n) + F(\tilde{u}, \tilde{\omega})).$ 每个时间步需要两次拉普拉斯算子的逆运算，通过多重网格 V-cycle 算法高效求解泊松方程 $\nabla^2 \psi = \omega$ 。
对称性监控：作者利用缩放对称性公式 $\rho x_{\epsilon}(t) = x_{\epsilon_0}(\rho^\alpha t)$ 来验证数值解的正确性。通过比较不同 $\epsilon$ 值下螺旋中心轨迹的缩放一致性，可以检测出数值误差或边界效应导致的偏差。这种方法提供了一种内置的诊断工具，用于评估计算结果的可靠性。
大规模并行计算：为了处理巨大的计算域和高阶精度带来的计算量，作者优化了代码的执行速度，使得能够长时间积分并处理精细网格。这对于捕捉缓慢演化的自相似结构至关重要。

4. 与之前工作的比较

与 Bressan 和 Shen (2021) 以及 Bressan 和 Murray (2020) 的先前工作相比，本文在以下几个方面做出了显著推进：

计算域的扩展：先前的研究主要在较小的有限域上进行数值模拟，例如 $[-0.2, 0.2] \times [-0.06, 0.06]$ ，并使用硬壁边界条件 ( $\psi=0$ )。这种设置可能会引入显著的边界效应，扭曲真实的自由空间解。本文明确指出，要捕捉精确自由空间问题的核心对称性，必须使用极其大的计算域。作者通过对比不同域大小下的结果，证明了小域模拟可能导致错误的对称性表现，从而质疑先前某些结论的普适性。
对称性的严格验证：在先前的工作中，对称性主要是作为理论假设或初步观察被提及。本文则将对称性验证作为数值方法的核心组成部分。通过系统地追踪螺旋中心的时间轨迹，并利用缩放不变性公式进行交叉验证，作者提供了更严格的证据来确认或否定自相似极限的存在。这种方法比单纯的涡度轮廓图更具说服力。
参数空间的全面探索：先前的研究主要关注特定的参数配置，以展示非唯一性的可能性。本文则对更广泛的参数对进行了系统性研究。结果显示，非唯一性仅在极少数参数组合下可能出现，这与先前暗示的更广泛适用性形成对比。这一发现强调了初始数据细节对解的唯一性的重要性。
数值方法的改进：本文开发的嵌套网格策略和高阶离散方案比先前使用的均匀网格有限差分法更加高效和准确。特别是，嵌套网格允许在关键区域获得极高的分辨率，同时保持全局计算的可行性。这使得作者能够更可靠地解析靠近原点的奇异行为和远场的衰减特性。
局限性与代价：尽管本文在数值精度和域大小上有所提升，但仍受限于计算资源，无法真正达到 $\epsilon \to 0$ 的极限。此外，虽然使用了大域，但边界条件仍需在无穷远处适当处理，这可能引入微小的误差。作者诚实地指出，他们的结果是“证据”而非“证明”，并呼吁未来工作结合更多理论分析来巩固这些数值发现。

3. 时空膨胀下 Favre 滤波可压缩 Navier–Stokes 系统的微分同变不变性

Diffeomorphism Invariance of the Favre-Filtered Compressible Navier–Stokes System under Spacetime Dilation

作者: Yuanya Li

类别: math.AP, math-ph

发布日期: 2026-09-10

arXiv ID: 2609.11566v1

摘要翻译: 我们证明，包含亚网格尺度闭合项的 Favre 滤波可压缩 Navier–Stokes 方程组，在任意参数 $$L>0$$ 下的坐标 $$(Lt,Lx_1,Lx_2,Lx_3)$$ 中写出的系统，可由标准坐标 $$(t,x_1,x_2,x_3)$$ 中的系统通过底层时空流形上的一个光滑微分同胚的作用而得到。证明始于对无穷维流形及光滑场构型 Fréchet 流形的自洽处理。随后，我们定义单参数膨胀群 $\{φ_L\}_{L>0}$ ，验证每个 $$φ_L$$ 均为微分同胚，并计算其切空间与余切空间作用。被诱导到构型丛截面上的拉回作用被证明与全导数算子交换，这意味着 $$φ_L$$ 的一阶射影延拓将系统的方程子流形映射到自身。逐方程验证——连续性、动量与能量——确认每一个残差均乘以公共的非零因子 $L^{-1}$ （或者等价地，在变换后的坐标中恒为零），从而解集得以保持。我们得出结论：由 $$L>0$$ 参数化的系统族满足微分同胚变换表达式，且对于所有 $$L$$ ，解空间是典范同构的。

Comment: 14 pages

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

1. 研究对象和背景

本文的研究对象是带有亚网格尺度（subgrid-scale, SGS）封闭项的Favre滤波可压缩Navier-Stokes方程组。该方程组通常用于大涡模拟（Large Eddy Simulation, LES），旨在通过滤除小尺度湍流运动来降低计算成本，同时保留大尺度流动的动力学特征。在标准坐标 $$(t, x_1, x_2, x_3)$$ 下，设 $$ar{ ho}$$ 为滤波密度， $ar{u}_i$ 为滤波速度， $$ar{p}$$ 为滤波压力， $$ar{E}$$ 为滤波总能量， $ar{ au}_{ij}$ 为滤波粘性应力张量， $ar{q}_i$ 为滤波热通量。此外，引入SGS模型变量： $au^{sgs}_{ij}$ 为SGS粘性应力， $H^{sgs}_i$ 为SGS热通量， $ar{ au}^{sgs}_{ij}$ 和 $ar{q}^{sgs}_i$ 可能代表更高级别的修正项或特定模型的闭合变量（文中记号略有重叠，统一视为场分量），以及源项 $\dot{Q}_eta$ 。

方程组由连续性方程、动量方程和能量方程组成。以连续性方程为例，其形式为：

\frac{\partial \bar{\rho}}{\partial t} + \frac{\partial (\bar{\rho} \bar{u}_i)}{\partial x_i} = 0

动量方程包含惯性项、压力梯度、粘性力及SGS应力散度：

\frac{\partial (\bar{\rho} \bar{u}_i)}{\partial t} + \frac{\partial}{\partial x_j} (\bar{\rho} \bar{u}_i \bar{u}_j + \bar{p}\delta_{ij} - \bar{\tau}_{ij}) + \frac{\partial \tau^{sgs}_{ij}}{\partial x_j} = 0

能量方程则涉及总能量的输运、做功及SGS热通量和源项：

\frac{\partial (\bar{\rho} \bar{E})}{\partial t} + \frac{\partial}{\partial x_i} ((\bar{\rho} \bar{E} + \bar{p})\bar{u}_i + \bar{q}_i - \bar{u}_j \bar{\tau}_{ij}) + \frac{\partial H^{sgs}_i}{\partial x_i} + \frac{\partial \sigma^{sgs}_{ij}}{\partial x_i} + \dot{Q}_eta = 0

在物理建模和数值分析中，理解方程在不同尺度变换下的行为至关重要。特别是当研究多尺度湍流或进行非均匀网格计算时，确定方程组是否在某种几何变换（如缩放/膨胀）下保持形式不变（即协变性或不变性）是一个基础且关键的问题。如果方程组在时空膨胀 $(t, x) \mapsto (Lt, Lx)$ 下具有微分同胚不变性，意味着解的空间在所有尺度 $$L > 0$$ 下是同构的，这为跨尺度分析提供了坚实的理论基础。

然而，传统的处理往往侧重于具体的代数变换验证，缺乏从无穷维流形和射影丛（jet bundle）视角的系统性几何框架描述。特别是对于包含复杂非线性项和SGS封闭项的系统，如何严格证明其在一般微分同胚作用下的映射性质，以及如何处理源项（source term）在坐标变换中的权重问题，是理论上的难点。本文试图填补这一空白，通过建立基于Fréchet流形的无限维微分几何框架，严格证明该Favre滤波系统在一参数膨胀群作用下的微分同胚不变性，并澄清源项变换的几何意义。

2. 主要定理或主要结果

本文的主要结果可以概括为以下几个定理和命题：

膨胀映射的微分同胚性质 (Theorem 2.2): 对于任意参数 $$L > 0$$ ，定义时空膨胀映射 $\phi_L: M \to M$ 为 $\phi_L(t, x_1, x_2, x_3) = (Lt, Lx_1, Lx_2, Lx_3)$ 。作者证明了 $\phi_L$ 是时空流形 $$M$$ 上的一个光滑微分同胚（ $C^\infty$ diffeomorphism）。其逆映射为 $\phi_L^{-1} = \phi_{1/L}$ 。雅可比矩阵为 $D\phi_L = L I_4$ ，行列式为 $$L^4 > 0$$ ，表明该变换保持定向。
一参数群结构 (Theorem 2.5): 集合 $\{\phi_L\}_{L>0}$ 构成 $$M$$ 上的光滑一参数微分同胚群。其无穷小生成元为向量场 $X_{dil} = t \frac{\partial}{\partial t} + x_i \frac{\partial}{\partial x_i}$ 。
截面拉回与导数的交换律 (Lemma 3.4 & Proposition 3.7): 定义了配置丛上的拉回作用 $\Phi_L(u) = \phi_L^* u$ ，其中 $(\phi_L^* u)(t,x) = u(Lt, Lx)$ 。核心引理指出，拉回运算与偏导数算子满足链式法则关系：
$\frac{\partial}{\partial y} (\phi_L^* f) = L \phi_L^* \left( \frac{\partial f}{\partial Y} \right)$
其中 $$y$$ 是原坐标， $$Y$$ 是对应的膨胀后坐标。进一步，作者将 $\phi_L$ 延拓到一阶射影丛 $$J^1E$$ 上，得到延拓映射 $j^1 \phi_L$ ，其在坐标 $$(p, u, u_t, u_x)$$ 下的作用为 $(\phi_L(p), u, L^{-1}u_t, L^{-1}u_x)$ 。
残差映射的等变性 (Theorem 4.7): 这是文章的核心结论。设 $P: C \to C^\infty(M, \mathbb{R}^5)$ 为整个方程组的残差映射（包括连续性、动量、能量方程）。作者证明了对于任意 $$L > 0$$ 和 $u \in C$ ，残差映射满足以下等变关系：
$$ P(\Phi_L(u)) = L \phi_L^* (P(u)) $$
具体而言：
- 连续性方程: $R_0[\Phi_L(u)](t,x) = L R_0[u](Lt, Lx)$ 。
- 动量方程: $R_i[\Phi_L(u)](t,x) = L R_i[u](Lt, Lx)$ 。
- 能量方程: 这里需要仔细处理源项。若将源项 $\dot{Q}_eta$ 视为场的一部分，其变换需符合密度权重。作者指出，严格的等变性要求源项作为密度变换，或者在被动坐标观点下，方程形式保持不变。最终结论是：解集被保持，即 $P(\Phi_L(u)) = 0 \iff P(u) = 0$ 。
解空间的典范同构: 基于上述等变性，作者得出结论：对于所有 $$L > 0$$ ，参数化方程族的解空间是典范同构的（canonically isomorphic）。这意味着在膨胀坐标 $$(T, X) = (Lt, Lx)$$ 下写出的方程组，与标准坐标下的方程组描述的是同一类几何对象，只是坐标表示不同。

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

本文的证明方法建立在严谨的无穷维微分几何和分析框架之上，主要技术障碍和创新点如下：

Fréchet流形框架的建立: 作者首先自洽地介绍了Fréchet空间和Fréchet流形的定义（第1节）。由于PDE的解空间通常是无穷维的，标准的Banach流形工具不够用。作者采用Bastiani-Kriegl-Michor定义的Fréchet光滑映射概念，将场配置空间 $C = C^\infty(M, V)$ 视为Fréchet流形。这一步骤为后续讨论无穷维空间上的微分同胚和切丛提供了严格的拓扑基础。特别地，作者引用了Nash-Moser隐函数定理（Proposition 1.5(v)），虽然在本例中由于线性性质未直接使用，但确立了局部子流形结构的合法性。
射影丛（Jet Bundle）语言的应用: 传统PDE不变性分析往往直接对代数表达式进行代换。本文的创新在于引入了射影丛 $$J^1E$$ 的语言（Definition 3.6）。通过将PDE系统定义为射影丛中的子流形 $E \subset J^1E$ ，微分同胚的作用可以通过射影延拓 $j^1 \phi_L$ 来描述。这种方法不仅使得“方程映射到自身”这一概念几何化，而且清晰地分离了场值本身（ $$u$$ ）及其导数（ $$u_t, u_x$$ ）的变换规则。例如，Proposition 3.7明确指出，在射影坐标下，导数部分乘以 $L^{-1}$ ，这与主动拉回视角下的 $$L$$ 因子形成对偶，避免了符号混淆。
主动与被动观点的统一: 作者在Remark 3.2中明确区分了“主动”（Active）和“被动”（Passive）两种观点。主动观点是将场拉回到原坐标系产生新场；被动观点是在新坐标系中重写同一几何场。证明过程中，作者交替使用这两种观点以确保一致性。特别是在处理能量方程的源项时（Section 4.3），这种区分至关重要。作者指出，源项 $\dot{Q}_eta$ 作为一个标量场，在被动坐标变换下不需要额外的Jacobian因子，但在主动拉回视图中，它是场元组的一部分，随整体变换。这种细致的几何解释解决了以往文献中关于源项权重的歧义。
逐项残差验证与等变性推导: 尽管有几何框架，作者仍然进行了详细的逐项代数验证（Section 4）。利用链式法则引理（Lemma 3.4），作者展示了每个偏导数项在拉回后都会提取出一个公因子 $$L$$ 。例如，对于连续性方程， $\frac{\partial}{\partial t}(\phi_L^* \bar{\rho}) = L \frac{\partial \bar{\rho}}{\partial T}|_{(Lt,Lx)}$ 。通过累加这些项，证明了残差映射的整体等变性 $P(\Phi_L(u)) = L \phi_L^* P(u)$ 。这一推导过程严密地连接了微分几何结构与具体的PDE代数形式。
源项处理的特殊技巧: 能量方程中包含非导数源项 $\dot{Q}_eta$ ，这使得简单的线性等变性受到挑战。作者在Definition 4.8中引入了“增强场元组”的概念，建议将源项视为具有适当密度的场，或者在被动坐标下将其视为不变的标量场。通过这种方式，作者证明了尽管源项不产生 $$L$$ 因子，但方程的零解集（solution set）依然被保持，从而确保了物理意义的不变性。

4. 与之前工作的比较

从提供的内容来看，本文并未显式列出大量特定的前人工作名称进行比较，但其方法论和结论隐含了对现有PDE不变性分析方法的推进。以下是基于文本内容的对比分析：

从代数验证到几何框架的推进: 许多关于Navier-Stokes方程或其滤波版本对称性的研究通常局限于具体的坐标变换代数操作，即直接代入 $$t'=Lt, x'=Lx$$ 并整理系数。本文的不同之处在于，它构建了一个完整的无穷维微分几何背景（Fréchet流形、射影丛）。这种处理方式不仅证明了当前的膨胀不变性，还为未来研究其他更复杂的微分同胚变换（如旋转、平移、甚至非线性的共形变换）提供了通用的理论模板。相比之下，纯代数方法难以推广到非线性的坐标变换，因为链式法则会变得极其复杂，而射影丛语言可以系统地处理高阶导数和联络项。
对SGS封闭项系统的严格处理: 标准的不可压或可压缩Navier-Stokes方程的对称性研究较为成熟。然而，Favre滤波系统引入了额外的SGS变量（如 $\tau^{sgs}_{ij}, H^{sgs}_i$ 等），这些变量通常依赖于滤波器的选择和本地流动状态。本文的一个贡献在于将这些SGS项视为配置丛纤维的一部分，并证明它们在膨胀变换下与其他基本场（密度、速度等）具有相同的变换规律（即作为普通标量或张量场的拉回）。这澄清了SGS模型在尺度变换下的行为，暗示如果SGS模型本身是基于局部梯度构建的，它们也应遵循类似的协变律。这一点在数值模拟的多尺度分析中具有重要的指导意义，而在以往的LES理论综述中，这种几何视角的阐述较少见。
源项变换权的澄清: 在处理包含源项的能量方程时，许多工程数学文献可能会忽略源项在坐标变换中的权重变化，或者错误地赋予其导数项相同的权重。本文通过区分主动/被动观点，并引入密度权重的概念（Definition 4.8），精确地指出了源项 $\dot{Q}_\beta$ 应作为标量场处理（在被动观点下不变），而不是像体积力那样变换。这一细节对于确保能量守恒定律在变换后的坐标系中形式正确至关重要。虽然这在经典力学中是已知的，但在结合SGS封闭项的现代湍流模型中，明确这一几何事实有助于避免数值实现中的潜在错误。
局限性与适用范围: 需要注意的是，本文主要关注的是形式不变性（Form Invariance）或协变性（Covariance），即方程的形式在变换下保持不变（除了可能的整体缩放因子），而不是指解的具体数值不变。此外，证明依赖于场的光滑性（ $C^\infty$ ），这对于湍流这种本质上是奇异的物理现象来说是一个理想化的假设。在实际应用中，解可能仅在分布意义下存在，此时Fréchet流形框架可能需要调整为Sobolev空间或其他更弱的拓扑结构。不过，对于理论分析和光滑初始值的短期演化，这一结果是严格且有效的。

综上所述，本文的主要推进在于将Favre滤波可压缩Navier-Stokes系统的膨胀不变性置于一个严格的无穷维微分几何框架内，通过射影丛语言统一了场及其导数的变换规则，并澄清了源项的几何权重，从而为多尺度湍流模拟的理论基础提供了更清晰的几何解释。

4. 具有水平粘性的二维 Navier–Stokes 方程的各向异性 Onsager 判据

An Anisotropic Onsager Criterion for the Two-Dimensional Navier–Stokes Equations with Horizontal Viscosity

作者: Siyu Liang

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11488v1

摘要翻译: 我们考虑在 $\mathbb T^2$ 或 $\mathbb R^2$ 上，具有水平粘性的二维 Navier–Stokes 方程的分布弱解 $u\in L^\infty_tL^2_x$ ，且未假设 $\partial_1u\in L^2_tL^2_x$ 。我们证明 Lions 可积性条件 $u\in L^4_tL^4_x$ 连同在非耗散方向上的 Onsager 临界正则性（即 Besov 正则性仅施加于竖直方向的 $u\in L^3_tB^{1/3,v}_{3,\infty}$ ）蕴含 $\partial_1u\in L^2_tL^2_x$ ，并给出关于水平耗散的定量界。因此，耗散正则性是方程的结果，而不必作为解定义的一部分。此外，若 $u\in L^3_tB^{1/3,v}_{3,c_0}$ ，则 $$u$$ 属于各向同性临界 Onsager 空间 $L^3_tB^{1/3}_{3,c_0}$ ，满足能量等式，并且在 $$L^2$$ 中连续。这里 $$c_0$$ 表示相应的二进 Besov 序列在高频处趋于零。该证明将水平能量通量吸收进耗散项与针对竖直能量通量的一维交换子估计相结合。

Comment: 无

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

1. 研究对象和背景

本文研究二维不可压缩 Navier-Stokes 方程在仅具有水平方向粘性的情形下的能量守恒与正则性. 具体考虑的方程为

\partial_t u + \text{div}(u \otimes u) + \nabla p - \nu \partial_1^2 u = 0, \nabla \cdot u = 0,

定义在环面 $\mathbb T^2$ 或全平面 $\mathbb R^2$ 上. 这种各向异性耗散模型出现在地球物理流体动力学中, 其中水平与垂直方向的有效粘性存在巨大差异. 该方程处于完全粘性 Navier-Stokes 方程与无粘 Euler 方程之间: 耗散项 $\nu \partial_1^2 u$ 仅控制水平方向的导数 $\partial_1 u$ , 而不控制垂直方向的导数 $\partial_2 u$ . 因此, 经典的 Ladyzhenskaya 估计不再适用, 且抛物型光滑化仅在水平频率方向起作用.

对于完全粘性 Navier-Stokes 方程, Lions [21] 证明了若解满足 $u \in L^4(0, T; L^4)$ , 则能量等式成立; Shinbrot [24] 将此条件推广至 $u \in L^p(0, T; L^q)$ 且 $2/p + 2/q \le 1, q \ge 4$ . 对于 Euler 方程, Onsager 猜想指出 $$1/3$$ 是临界正则性指数: Constantin, E 和 Titi [9] 证明了若 $u \in L^3(0, T; B^{\alpha}_{3,\infty})$ 且 $\alpha > 1/3$ , 则能量守恒; Cheskidov, Constantin, Friedlander 和 Shvydkoy [6] 进一步将其推广至临界空间 $L^3(0, T; B^{1/3}_{3,c_0})$ , 其中 $$c_0$$ 表示对应的 dyadic Besov 序列在高频处趋于零. 低于临界指数时, 存在不守恒能量的弱解.

在二维完全粘性情形下, 每个 Leray-Hopf 解均满足能量等式. 然而, 对于仅含水平粘性的方程, 既有的论证均不适用. Demmel 和 Wiedemann [10] 最近在 $\mathbb R^2$ 上证明了: 若假设解已满足自然能量级别 $u \in L^\infty(0, T; L^2)$ 且 $\partial_1 u \in L^2(0, T; L^2)$ , 则能量等式成立. 但该文未处理周期域 $\mathbb T^2$ , 且其结论依赖于预先假设 $\partial_1 u \in L^2(0, T; L^2)$ .

本文的核心问题在于: 从更弱的分布解 $u \in L^\infty(0, T; L^2)$ 出发, 在不预先假设 $\partial_1 u \in L^2(0, T; L^2)$ 的情况下, 确定何种额外正则性能保证达到这一能量级别. 作者试图填补的空白是: 结合 Lions 可积性条件与非耗散方向（即垂直方向）上的 Onsager 临界正则性, 能否推导出耗散方向（即水平方向）的正则性, 从而使得耗散正则性成为方程本身的推论, 而非解定义的一部分.

2. 主要定理或主要结果

设 $$T > 0$$ , $\nu > 0$ , $u_0 \in L^2_\sigma(\Omega)$ , 其中 $\Omega$ 为 $\mathbb T^2$ 或 $\mathbb R^2$ . 令 $$u$$ 为方程 (2.2) 在 $$[0, T]$$ 上以 $$u_0$$ 为初值的弱解（见 Definition 2.6）, 并假设

u \in L^4(0, T; L^4(\Omega)), \tag{2.8}

u \in L^3(0, T; B^{1/3, v}_{3,\infty}(\Omega)). \tag{2.9}

这里 $B^{s, v}_{p, r}$ 表示仅在垂直变量 $$x_2$$ 上进行 Littlewood-Paley 分解的各向异性 Besov 空间.

主要结果 (Theorem 2.8) 分为两部分:

(i) 在上述假设下, 有 $\partial_1 u \in L^2(0, T; L^2(\Omega))$ , 且存在普适常数 $$C$$ （独立于 $\nu, T, u$ ）使得

\nu \int_0^T \|\partial_1 u(t)\|^2_{L^2} dt \le \|u_0\|^2_{L^2} + \frac{C}{\nu} \|u\|^4_{L^4(0,T; L^4)} + C \int_0^T \|u(t)\|^3_{B^{1/3, v}_{3,\infty}} dt. \tag{2.10}

这表明水平耗散量可由初值能量、Lions 范数以及垂直方向 Besov 范数定量控制.

(ii) 若进一步假设

u \in L^3(0, T; B^{1/3, v}_{3, c_0}(\Omega)), \tag{2.11}

则:

$$u$$ 属于各向同性临界 Onsager 空间 $L^3(0, T; B^{1/3}_{3, c_0}(\Omega))$ ;
能量等式 $\frac{1}{2}\|u(t)\|^2_{L^2} + \nu \int_0^t \|\partial_1 u(s)\|^2_{L^2} ds = \frac{1}{2}\|u_0\|^2_{L^2} \tag{2.13}$ 对所有 $t \in [0, T]$ 成立;
$u \in C([0, T]; L^2(\Omega))$ .

值得注意的是, 上述结论均未假设各向同性临界正则性或 $$L^2$$ 连续性, 而是由方程结构与给定的各向异性条件推导得出. 此外, 通过插值, 假设 (2.8) 可由 $u \in L^p(0, T; L^q(\Omega))$ 且 $2/p + 2/q \le 1, q \ge 4$ 推出, 因此涵盖了 Shinbrot 型条件. 特别地, 作者在联合工作 [20] 中构造的全局确定性解（初值 $u_0 \in \tilde H^{0,1}(\Omega)$ ）在每个有限时间区间上均满足 (2.8) 和 (2.11).

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

证明的核心技术障碍在于: 耗散项仅作用于水平方向, 无法直接控制垂直方向的非线性项. 传统方法依赖 Ladyzhenskaya 不等式或各向同性抛物型估计, 在此处均失效. 作者的真正创新点在于将水平能量通量吸收进耗散项, 并利用一维交换子估计处理垂直能量通量.

证明框架始于对场 $$u$$ 在水平频率 $$N$$ 和垂直频率 $$Q$$ 处截断后的能量平衡. 具体而言, 使用一维 Littlewood-Paley 分解定义水平截断算子 $$S_N^h$$ 和垂直截断算子 $$S_Q^v$$ , 并考虑双重截断 $$S_N^h S_Q^v u$$ . 能量平衡中的关键项包括水平通量和垂直通量.

首先, 水平通量的处理是关键一步. 由于耗散项 $\nu \partial_1^2 u$ 的存在, 水平方向的能量输运可以被耗散所吸收. 交换子误差由 $\|u\|_{L^4}^2$ 控制, 这得益于 Lions 可积性条件 $u \in L^4_t L^4_x$ . 这一步骤确保了水平方向的能量不会发散, 并为后续取极限 $N \to \infty$ 提供了基础.

其次, 垂直通量的处理更为微妙. 由于没有垂直耗散, 必须利用垂直方向的正则性来控制此项. 作者引入了一维交换子估计, 其形式类似于 Constantin-E-Titi [9] 和 Cheskidov et al. [6] 中的技术, 但适配于各向异性设定. 具体地, 定义垂直增量 $\delta_z^v f(x_1, x_2) := f(x_1, x_2 - z) - f(x_1, x_2)$ , 并建立如下估计:

\|\delta_z^v S_J^v g\|_{L^3} \lesssim 2^{2J/3} \|g\|_{B^{1/3, v}_{3,\infty}} |z|, \tag{3.19}

\|\delta_z^v (g - S_J^v g)\|_{L^3} \lesssim \left(\sup_{j>J} \epsilon_j(g)\right) |z|^{1/3}, \tag{3.20}

其中 $\epsilon_j(g) := 2^{j/3} \|\Delta_j^v g\|_{L^3}$ . 这些估计表明, 当 $$|z|$$ 较小时, 垂直增量的 $$L^3$$ 范数受控, 且控制常数依赖于 $B^{1/3, v}_{3,\infty}$ 范数. 指数 $$1/3$$ 的出现正是 Onsager 临界指数的体现.

通过让 $N \to \infty$ 然后 $Q \to \infty$ , 作者得到水平耗散估计 (2.10). 随后, 能量等式直接从垂直截断后的平衡推导得出, 无需使用各向同性空间. 最后, 在 $$c_0$$ 细化条件下, 作者证明了 $$u$$ 属于各向同性临界 Onsager 空间 $L^3(0, T; B^{1/3}_{3, c_0})$ , 并由此得出 $$L^2$$ 连续性和能量等式.

关键技术工具包括:

一维 Littlewood-Paley 分解: 分别在水平和垂直变量上应用, 以分离不同方向的频率行为.
Bernstein 不等式: 用于控制单变量 Fourier 支集受限函数的 $$L^p$$ 范数.
交换子估计: 特别是 Lemma 3.8 中的一维交换子估计, 是处理垂直通量的核心.
核函数估计: Lemma 3.7 中对截断核 $$K_M$$ 的矩估计, 确保了一维估计在各向异性设定下的有效性.
插值与密度论证: 用于从截断解过渡到完整解, 并确保弱收敛性.

4. 与之前工作的比较

与既有 PDE/分析理论文献相比, 本文的主要推进体现在以下几个方面:

弱化先验假设: Demmel 和 Wiedemann [10] 的结果要求预先假设 $\partial_1 u \in L^2(0, T; L^2)$ , 即解已处于自然能量级别. 本文则从更弱的分布解 $u \in L^\infty(0, T; L^2)$ 出发, 不假设任何耗散方向的正则性, 而是通过 Lions 可积性条件和垂直方向 Onsager 临界正则性推导出 $\partial_1 u \in L^2(0, T; L^2)$ . 这意味着耗散正则性是方程结构的后果, 而非解定义的一部分, 从而显著放宽了对解的先验要求.
统一处理两个域: Demmel 和 Wiedemann [10] 仅处理了 $\mathbb R^2$ 的情形, 而本文同时处理了 $\mathbb T^2$ 和 $\mathbb R^2$ . 尽管在技术上两者相似, 但周期域涉及 Fourier 级数而非 Fourier 变换, 需要额外的周期性 Bernstein 不等式和 Poisson 求和公式. 作者通过统一的核函数表示 (Definition 3.1) 和估计 (Lemma 3.7) 实现了这一统一.
各向异性正则性的精确刻画: 既有工作通常考虑各向同性正则性, 如 $u \in L^3(0, T; B^{1/3}_{3, c_0})$ . 本文则引入了仅在垂直方向测量的各向异性 Besov 空间 $B^{1/3, v}_{3, \infty}$ 和 $B^{1/3, v}_{3, c_0}$ , 并证明这些较弱条件足以推导出各向同性临界正则性. 这揭示了方程结构如何补偿非耗散方向的正则性不足.
与全局适定性结果的联系: 作者在联合工作 [20] 中建立了初值 $u_0 \in \tilde H^{0,1}(\Omega)$ 的全局适定性, 这需要初值具有一个完整的垂直 $$L^2$$ 导数. 本文则表明, 对于已有的弱解, 只需其在垂直方向具有 $$1/3$$ 阶的时间积分 Besov 正则性即可保证能量等式. 这提供了一个独立的正则性判据, 不依赖于初值的 Sobolev 正则性.
代价与局限: 本文的结论依赖于 Lions 可积性条件 $u \in L^4(0, T; L^4)$ 和垂直方向 Onsager 临界正则性. 虽然这些条件在物理上是合理的, 但在数学上仍属较强假设. 此外, 能量等式的证明直接依赖于水平耗散估计, 而未使用各向同性 Onsager 空间, 这表明各向异性结构在能量守恒中的作用比先前认为的更为根本. 未来工作可能探索是否能进一步弱化 Lions 条件或降低垂直方向的正则性要求.

5. 一种双极 Euler-Poisson 系统的全局存在性与时间衰减：含一种无压且无阻尼流体

Global existence and time decay for a bipolar Euler-Poisson system with one pressureless and undamped fluid

作者: Young-Pil Choi, Houzhi Tang, Weiyuan Zou

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.10959v1

摘要翻译: 我们研究三维双极 Euler–Poisson 系统的 Cauchy 问题，其中一种流体是无压且无阻尼的，而另一种流体则受动量弛豫作用。对于常数平衡态的充分小的光滑扰动，在关于无压流体的初始速度无旋性的假设下，我们证明了光滑解的全局存在性与唯一性，并给出了代数时间衰减估计。主要难点在于无压流体的速度仅通过 Poisson 耦合间接耗散，并且该机制在高频处强烈退化，导致正则性损失结构。我们通过结合精细的 Green 函数估计、低-中-高频分解以及适应于非对称正则性层级的高阶非线性能量估计来克服这一困难。该结果建立了在这种非对称情形下的全局小数据理论，在此情形中同一流体同时缺乏压力与阻尼。

Comment: 31 pages

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

1. 研究对象和背景

本文研究三维双极 Euler-Poisson 系统的柯西问题。该模型用于描述半导体器件或等离子体中两种带电粒子的流体动力学行为，其中两个组分通过自洽电场相互作用。具体而言，系统由以下方程组成：

\begin{cases} \partial_t \rho_1 + \text{div}(\rho_1 u_1) = 0, \\ \partial_t (\rho_1 u_1) + \text{div}(\rho_1 u_1 \otimes u_1) = \rho_1 \nabla \Phi, \\ \partial_t \rho_2 + \text{div}(\rho_2 u_2) = 0, \\ \partial_t (\rho_2 u_2) + \text{div}(\rho_2 u_2 \otimes u_2) + \nabla P(\rho_2) = -\rho_2 \nabla \Phi - \rho_2 u_2, \\ \Delta \Phi = \rho_1 - \rho_2. \end{cases}

这里 $\rho_i$ 和 $$u_i$$ ( $$i=1,2$$ ) 分别表示密度和速度， $\Phi$ 是电势。压力定律为 $P(\rho_2) = \rho_2^\gamma$ ，其中 $\gamma > 1$ 。

该系统的一个显著特征是高度的不对称性：第一个流体组分既无压强项（pressureless）也无动量松弛项（undamped），而第二个流体组分同时具有压强和阻尼。这种设定源于物理上一种组分温度极低或弛豫率极小的极限情形。在标准的对称双极模型中，压强和阻尼共同为两个组分提供稳定性机制；然而在此处，由于第一种流体缺乏这两种机制，其耗散只能通过 Poisson 耦合间接从第二种流体传递过来。这种间接耗散机制在高频率下表现出强烈的退化性，导致了所谓的“正则性损失结构”（regularity-loss structure）。这意味着，通常依赖于均匀耗散的估计方法在此失效，因为高频模态的衰减率与波数 $|\xi|$ 相关，表现为 $e^{-\eta t |\xi|^{-4}}$ 的形式，而非均匀的指数衰减。因此，要在 $\mathbb{R}^3$ 上建立全局光滑解的存在性并证明其长时间代数衰减，需要克服这一非对称耗散带来的核心分析障碍。

2. 主要定理或主要结果

本文的主要结果概括于定理 1.1 中。假设初始速度 $u_{10}$ 是无旋的（即 $\text{curl } u_{10} = 0$ ），并且初始数据 $(a_{10}, a_{20}, u_{10}, u_{20})$ 在适当的 Sobolev 空间范数下足够小。具体地，定义小量参数：

\varepsilon_N := \|a_{20}\|_{\dot{H}^{-1}} + \left\|a_{10} - \frac{1}{c}a_{20}\right\|_{\dot{H}^{-1}} + \|a_{10}\|_{H^{N-1}} + \|(u_{10}, a_{20}, u_{20})\|_{H^N},

其中 $c = \sqrt{\gamma}$ ， $N \geq 10$ 。若 $\varepsilon_N \leq \varepsilon_0$ （某个依赖于 $\gamma$ 和正性常数的小量），则柯西问题存在唯一的全局经典解。解的正则性满足：

a_1 \in C([0, \infty); H^{N-1}) \cap C^1([0, \infty); H^{N-2}), (u_1, a_2, u_2) \in C([0, \infty); H^N) \cap C^1([0, \infty); H^{N-1}).

此外，解满足如下的代数时间衰减估计：

\|(a_1, a_2)(t)\|_{L^2} \leq C \varepsilon_N (1+t)^{-1/2},

\|\nabla(a_1, a_2)(t)\|_{H^1} \leq C \varepsilon_N (1+t)^{-1},

\|(\text{div } u_1, \text{div } u_2)(t)\|_{H^2} \leq C \varepsilon_N (1+t)^{-1},

\|(a_1, u_1, a_2, u_2)(t)\|_{W^{1,\infty}} \leq C \varepsilon_N (1+t)^{-5/4}.

这些结果表明，尽管缺乏直接的耗散项，但通过 Poisson 耦合的间接作用，系统仍然展现出良好的长时间衰减性质，且衰减速率与初始数据的 Sobolev 范数线性相关。

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

为了克服非对称耗散导致的正则性损失，作者采用了结合精细 Green 函数估计、低-中-高频分解以及高阶非线性能量估计的策略。

首先，作者对线性化系统进行频谱分析。通过引入纵向变量 $q_i = \text{div } u_i$ ，分析了不同频率区域的本征值分布。在低频区，存在一个扩散模态，其余模态稳定；在中频区，谱与虚轴均匀分离；而在高频区，两个本征值的实部以 $-|\xi|^{-4}$ 的速率趋近于零。这种高频退化意味着直接的能量估计无法提供 $$u_1$$ 的时间可积耗散。为了弥补这一缺陷，作者利用 Green 函数的精细渐近展开，证明了线性解的时间衰减是以牺牲空间导数为代价获得的，即获得 $(1+t)^{-l/4}$ 的时间衰减需要额外的 $$l$$ 个高频导数。

其次，为了闭合高阶能量估计，作者构造了一个适应于非对称正则性层级的能量泛函。关键变量是电荷变量的线性部分 $b = a_1 - \frac{1}{c}a_2$ 。通过将 $$b$$ 置于比 $$(u_1, a_2, u_2)$$ 低一阶的 Sobolev 空间中，作者利用了方程结构中最高阶导数的精确抵消。具体来说，在计算高阶能量时，Poisson 力项在 $$b, u_1, u_2$$ 的方程中产生的最高阶项相互抵消，从而避免了额外增加 $$a_1$$ 的正则性要求。

最后，作者利用了无旋条件 $\text{curl } u_1 = 0$ 。由于 Poisson 力是梯度场，它只作用于 $$u_1$$ 的纵向分量，而不影响横向分量（涡度）。通过分析涡度方程，作者证明了如果初始涡度为零，则涡度在整个演化过程中保持为零。这使得可以通过 Hodge 分解将 $$u_1$$ 重构为其散度的函数，从而将针对散度的 Green 函数估计推广到完整的速度场。结合低-中-高频分解，作者在 Duhamel 积分中分别处理不同频率区域的非线性项，最终通过 Bootstrap 论证闭合了全局存在性和衰减估计。

4. 与之前工作的比较

与既有文献相比，本文的研究对象更为退化。在标准的双极 Euler-Poisson 系统中（如文献 [1], [15], [23]），两个流体组分均受到压强和阻尼的作用，这提供了强大的稳定机制，使得全局存在性和衰减估计相对容易建立。相比之下，本文中的第一个流体既无压强也无阻尼，这种双重缺失显著增加了问题的难度。

最接近的相关工作是文献 [25]，该文研究了部分阻尼的两流体 Euler-Poisson 系统在 $\mathbb{T}^3$ 上的全局存在性。然而，[25] 中的模型保留了两个组分的压强，仅在一个速度上施加阻尼。本文的系统更为退化，因为无阻尼组分也是无压的。此外，[25] 是在紧致域 $\mathbb{T}^3$ 上研究的，主要关注全局正则性；而本文是在非紧域 $\mathbb{R}^3$ 上研究的，必须同时处理全局正则性和长时间衰减，这需要更精细的频率分析和 Green 函数估计来捕捉耗散在不同频率下的传递机制。

对于单流体 Euler-Poisson 系统，已有工作（如文献 [8]）利用 Poisson 耦合产生的色散 Klein-Gordon 结构建立了无旋小数据的全局理论。但对于双流体系统，尤其是当其中一个流体完全缺乏耗散机制时，现有的理论尚不完善。本文填补了这一空白，首次建立了在这种高度非对称、双重退化情形下的全局小数据理论。与文献 [16] 中关于部分阻尼 Euler-Maxwell 系统的结果相比，本文不仅处理了阻尼的非对称性，还引入了无压特性，导致高频动态和正则性结构的根本变化。

6. Sobolev 空间中具有退化临界点的二维非定常 Prandtl 系统适定性

Well-posedness of the two-dimensional unsteady Prandtl system in Sobolev space with degenerate critical points

作者: Shi-Yong Zhu, Ya-Guang Wang

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.10936v1

摘要翻译: 本文致力于经典 Prandtl 方程在有限阶 Sobolev 空间中的适定性问题。对于具有退化临界点和一般流出条件的初始数据，我们通过引入一种新的迭代格式和线性消去法，得到了 Prandtl 方程在 Sobolev 空间中关于时间的局部存在性与唯一性。这一结果表明，Oleinik 单调性条件并非 Prandtl 方程在 Sobolev 空间中适定的必要条件，并为证明二维非定常边界层中零剪应力不一定导致边界层分离提供了证据。

Comment: 无

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

1. 研究对象和背景

边界层理论是现代流体力学的基石之一，其核心数学问题在于证明 Navier-Stokes 方程在低雷诺数极限下的渐近展开的有效性。这一有效性依赖于边界层内主导方程——Prandtl 方程的适定性。本文关注的是二维非稳态 Prandtl 方程在周期性水平方向上的半无限域 $\Omega = \{(x, y) | x \in \mathbb{T}, y \in [0, +\infty)\}$ 上的经典解存在唯一性问题。该方程组由动量方程和不可压缩条件组成：

\begin{cases} \partial_t u + u \partial_x u + v \partial_y u = \partial_{yy} u - \partial_x P, \\ \partial_x u + \partial_y v = 0, \\ u|_{y=0} = v|_{y=0} = 0, \\ \lim_{y \to +\infty} u = U(t, x), \\ u|_{t=0} = u_{in}(x, y). \end{cases}

其中 $$(u, v)$$ 是边界层内的速度场， $$U(t, x)$$ 和 $$P(t, x)$$ 分别是外部 Euler 流动在边界处的切向速度和压力迹，它们通过 Bernoulli 定律 $\partial_t U + U \partial_x U = -\partial_x P$ 相关联。

Prandtl 方程的主要分析困难源于其退化抛物性质和非局部性。由于不可压缩条件 $v = -\int_0^y \partial_x u \, d\eta$ ，非线性项 $v \partial_y u$ 涉及对 $$u$$ 的水平导数的积分，导致在能量估计中出现水平导数的损失（loss of tangential derivatives）。这种损失使得在无结构假设的情况下，Prandtl 方程在有限阶 Sobolev 空间中可能是不适定的。例如，文献 [4] 证明了当背景流速具有非退化临界点时，线性化问题在 Sobolev 空间中强不适定。

历史上，Oleinik 通过引入单调性条件 $\partial_y u > 0$ ，利用 Crocco 变换将方程转化为单个退化抛物方程，并结合极大值原理，建立了 Hölder 空间中的局部适定性。近期工作 [1] 和 [30] 进一步在 Sobolev 空间中重建了这一结果，关键在于引入了“好未知量”（good unknown） $W = \partial_y u \partial_y (\frac{u-U}{\partial_y u})$ ，利用其对偶结构抵消了导数损失。然而，这些方法严格依赖于 Oleinik 单调性条件。当初始数据 $u_{in}$ 在法向变量 $$y$$ 上具有退化临界点（即 $\partial_y u_{in}(y^*) = \partial_{yy} u_{in}(y^*) = 0$ ）时，单调性被破坏，上述基于 $$W$$ 的方法失效。

尽管存在不适定性的反例，但文献 [4] 的不适定性结论主要针对非退化临界点。对于退化临界点的情况，适定性是否成立尚属开放问题。物理上，临界点的存在意味着壁面剪切应力为零，传统观点认为这可能导致边界层分离。然而，实验和数值计算表明，二维非稳态边界层中零剪切应力并不必然导致分离。因此，从数学上证明在退化临界点条件下 Prandtl 方程的适定性，不仅填补了理论空白，也为理解非稳态边界层的稳定性提供了重要依据。本文旨在解决这一空白，即在初始数据具有退化临界点且外流一般的情况下，建立 Prandtl 方程在 Sobolev 空间中的局部存在唯一性。

2. 主要定理或主要结果

本文的主要结果是定理 2.1，它证明了在特定结构假设下，二维非稳态 Prandtl 方程在有限阶 Sobolev 空间中的局部存在唯一性。

主要假设：

初始数据的退化临界点结构 (Assumption 2.1)：设存在 $$y^* > 0$$ 及区间 $[\check{y}, \hat{y}]$ 包含 $$y^*$$ ，使得初始速度 $u_{in}(x, y)$ 满足：
$\partial_y u_{in}(x, y^*) = \partial_{yy} u_{in}(x, y^*) = 0, \quad \forall x \in \mathbb{T},$
$\partial_{yyy} u_{in} \ge C > 0 \quad \text{在 } [\check{y}, \hat{y}] \text{ 上},$
$\partial_y u_{in} \ge c_{in} \varpi(y) e^{-\delta \langle y \rangle}, \quad \forall x \in \mathbb{T}, y \ge 0,$
其中 $\varpi(y)$ 是一个非负光滑函数，在 $$y^*$$ 附近表现为 $$(y-y^*)^2$$ ，在其他区域为 1。这表明 $u_{in}$ 在 $$y^*$$ 处有二阶零点，但在其他区域保持正性且具有指数衰减。
初始数据的高阶正则性与权重控制 (Assumption 2.2)：定义修正后的初始扰动 $\tilde{u}_{in} = u_{in} - u_0(0, x, y)$ ，其中 $$u_0$$ 是一个精心构造的背景函数。要求 $\tilde{u}_{in}$ 在加权 Sobolev 空间 $\hat{H}^s_\psi$ 中具有有限范数，并且其高阶水平导数受到特定的奇性权重 $\omega_a^{in}$ 和 $P_{in}$ 的控制。具体地，对于某个奇整数 $s \ge 11$ ，存在常数 $C_{in}$ 使得：
$\|\tilde{u}_{in}\|_{\hat{H}^s_{\psi_{in}}} + \|\omega_{\epsilon_0}^{in} P_{in} \partial_y \partial_x^s \tilde{u}_{in}\|_{L^2_{\hat{\psi}_{in}}} + \dots \le C_{in},$
这里权重函数 $P_{in}$ 和 $\omega_a^{in}$ 专门设计用于处理 $$y^*$$ 附近的奇性。
外流条件：外流速度 $$U(t, x)$$ 是正的，且属于适当的光滑函数空间 $C([0, T]; H^{s+2}(\mathbb{T})) \cap C^1([0, T]; H^{s+1}(\mathbb{T})) \cap \bigcap_{i=2}^{(s-3)/2} C^i([0, T]; H^{s-2i}(\mathbb{T}))$ 。

主要结论： 在上述假设下，存在 $t^* \in (0, T]$ ，使得 Prandtl 方程在 $$[0, t^*]$$ 上存在唯一的经典解 $$u$$ 。令 $$w = u - u_0$$ ，则 $$w$$ 满足如下正则性：

w \in L^\infty((0, t^*), \hat{H}^s_\psi(\Omega)) \cap L^2((0, t^*), H^s_{\hat{\psi}} \cap H^{s+1}_{\check{\psi}}(\Omega)),

其中权重函数定义为 $\psi(t, y) = e^{(\frac{3}{4} - \Lambda t)\delta \langle y \rangle}$ ， $\hat{\psi} = \langle y \rangle^{-1/2} \psi$ ， $\check{\psi} = \langle y \rangle^{-1} \psi$ ， $\langle y \rangle = 1+y$ 。这里的权重 $\psi$ 随时间演化，旨在闭合能量估计并控制非局部项。

此外，文章指出该结果可以直接推广到有限个退化临界点的情形。这一结果明确表明，Oleinik 单调性条件并非 Prandtl 方程在 Sobolev 空间中适定的必要条件。

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

本文的核心技术挑战在于如何在缺乏全局单调性的情况下克服非局部项导致的水平导数损失。作者提出了一种新颖的迭代方案结合线性抵消机制，具体分为以下几个关键步骤：

1. 修正增量迭代方案 (Corrected Increment Iteration Scheme)： 直接对原方程进行能量估计是不可行的，因为单调性缺失。作者采用一种增量迭代方法，将解 $$u$$ 分解为一系列修正增量的和： $u = \sum_{n=0}^\infty u_n$ 。其中 $$u_0$$ 是一个预先构造的背景函数，它在临界点 $$y^*$$ 附近添加了“补丁”以恢复某种形式的单调性，并在远处匹配外流。对于 $n \ge 1$ ， $$u_n$$ 满足一系列半线性修正增量问题 (AEn)。这些问题的设计目的是使其结构接近线性且无退化，从而便于建立先验估计。例如，第一阶修正增量 $$u_1$$ 满足的问题 (AE1) 中，非线性项被移至源项，同时引入了人工粘性项 $\epsilon_0^2 \vartheta \partial_{xx} u_1$ 和一个截断的非局部项 $h \partial_y^{-1}[\partial_x u_1] \partial_y u_1$ 。截断函数 $$h(y)$$ 仅在远离边界的区域取值为 1，从而避免了在边界附近因非局部性导致的衰减损失。

2. 线性抵消与奇异权重的引入： 在修正增量问题中，虽然通过添加“补丁” $\varsigma$ 恢复了局部的类 Oleinik 结构，但当补丁参数趋于零时，标准的“好未知量”无法提供一致估计。作者的关键创新在于观察到，在退化临界点附近，随着时间演化，单调性会立即恢复（即 $\partial_y u$ 的行为类似于 $$C((y-y^*)^2 + t)$$ ）。利用这一观察，作者引入了一组新的奇异权重函数 $$P$$ 和 $\omega_a$ ，它们精确地捕捉了 $$y^*$$ 附近的奇性行为。定义新未知量 $$U = P w$$ ，其中 $$w$$ 是与最高阶水平导数相关的量。通过证明 $$P^2 w$$ 可以一致地控制 $P^2 \partial_x^s (u-u_0)$ ，作者建立了关于 $$U$$ 的一致加权能量估计。这一步骤克服了 $$P$$ 在初始时刻的强奇异性，使得估计在补丁参数和人工粘性参数趋于零时保持一致。

3. 一致先验估计与后验验证： 为了证明迭代序列的收敛性，作者首先在第 4 节建立了修正增量问题 (AEn) 的一致先验估计（Uniform a priori estimates）。这包括低阶和高阶水平导数的估计，以及 $L^\infty_t$ 范数的控制。这些估计不依赖于迭代步数 $$n$$ ，仅依赖于初始数据和外流的界。随后，在第 5 节，作者对这些先验估计进行了后验验证（A posteriori estimates），通过归纳法证明如果第 $$n$$ 步的近似解满足某些小性条件，那么第 $$n+1$$ 步的增量 $u_{n+1}$ 也将满足相应的小性条件。这确保了迭代过程可以逐步推进。

4. 收敛性论证： 最后，作者证明了修正增量序列 $$u_n$$ 在适当的 Sobolev 空间中构成柯西列。关键在于证明 $\|u_n\|_{L^\infty_t \hat{H}^s_\psi}$ 以几何级数速率衰减（如公式 3.8 所示）。由于 $u_n|_{t=0}=0$ ( $n \ge 2$ ) 且源项 $F_{n-1}$ 中除了一个线性项外都是非线性的，通过选取足够小的时间 $$t^*$$ 和足够小的补丁参数，可以实现这种快速衰减。最终，取极限得到原非线性问题的经典解。

4. 与之前工作的比较

本文工作与既有文献相比，主要在以下几个方面取得了显著推进：

1. 突破 Oleinik 单调性条件的限制： 经典的 Prandtl 方程适定性理论，无论是 Oleinik 的原始工作 [31, 32] 还是近期的 Sobolev 空间工作 [1, 30]，都严格依赖于单调性条件 $\partial_y u > 0$ 。这些工作通过引入“好未知量”来抵消导数损失，但该技巧在单调性失效时不再适用。本文首次证明了在初始数据具有退化临界点（即违反单调性）的情况下，Prandtl 方程在 Sobolev 空间中仍然是局部适定的。这一结果直接回应了文献 [4] 等指出的不适定性风险，明确了不适定性并非在所有非单调情形下都会发生，而是取决于临界点的退化性质。

2. 与 Gevrey 类工作的对比： 由于在 Sobolev 空间中缺乏单调性会导致不适定或需要极强的正则性，许多研究者转向 Gevrey 类函数空间来建立适定性 [2, 6, 25, 27]。Gevrey 类的解析性余量可以吸收导数损失。相比之下，本文工作在更一般的 Sobolev 空间中，无需假设解析性或 Gevrey 正则性，仅要求有限阶 Sobolev 正则性（ $s \ge 11$ ）。这使得结果在物理应用中更具普适性，因为实际流体通常不具备解析性。

3. 与非稳态边界层分离理论的关联： 在物理意义上，临界点对应于零剪切应力。传统观点认为零剪切应力是导致边界层分离的充分条件。然而，本文的数学结果表明，在二维非稳态情形下，即使存在零剪切应力（退化临界点），方程依然适定，至少在短时间内不会发生数学意义上的“破裂”。这与稳态边界层不同，后者中零剪切应力往往标志着分离的开始。这一发现与实验观察和数值模拟结果一致，为非稳态边界层分离的判据提供了新的数学视角。值得注意的是，本文结果并不与文献 [10] 中的非线性不稳定结果冲突；后者指出在某些特殊背景下渐近展开可能失效，但这并不否定 Prandtl 方程本身的适定性。

4. 技术方法的创新与代价： 与 [1] 和 [30] 使用固定的“好未知量”不同，本文引入了依赖于时间和空间的动态权重函数 $\psi(t, y)$ 以及针对临界点奇性的特殊权重 $$P$$ 和 $\omega_a$ 。这种复杂的权重结构是处理退化奇性的必要代价。此外，迭代方案中引入的人工粘性和“补丁”虽然在极限过程中消失，但在中间步骤增加了分析的复杂度。尽管如此，这种方法成功地将适定性范围从单调类扩展到了具有退化临界点的非单调类，代表了 Prandtl 方程分析理论的重要进展。

7. 二维趋化性Navier-Stokes方程的对数敏感性边界层. II. 粘性消失极限

Boundary layer of 2D Chemotaxis Navier-Stokes equations with logarithmic Sensitivity. II. viscous vanishing limit

作者: Hui Wang, Lingling Zhao

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11027v1

摘要翻译: 这是一项关于二维半空间中耦合趋化性-Navier-Stokes系统边界层解的两部分工作的第二部分。在本文中，我们研究二维半平面内，奇异趋化性-流体方程的边界层解关于化学扩散-粘性参数 $\varepsilon$ 的收敛性，且满足滑移边界条件。更具体地，我们证明当 $\varepsilon>0$ （粘性对流系数）时，边界层随着 $\varepsilon\rightarrow0$ 收敛于外层（ $\varepsilon=0$ 时的解）与内层的叠加。外层和内层剖面如第一部分\cite{WWZ}中所述被显式导出。此外，耦合趋化性-Navier-Stokes系统在法向Sobolev空间中的适定性结果将在附录中给出。它们回答了两部分工作中第一部分提到的问题。本研究有助于理解实验中观察到的好氧细菌向水-气界面的趋化运动，并丰富趋化性流体模型中边界层的理论结果。

Comment: 无

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

1. 研究对象和背景

本文研究的是二维半平面 $\Omega = \mathbb{R}^2_+ = \{(x, y) \in \mathbb{R}^2 | y > 0\}$ 上耦合的趋化-纳维-斯托克斯 (Chemotaxis-Navier-Stokes) 系统在滑移边界条件下的粘性消失极限问题。该模型旨在描述好氧细菌在流体中向水-气界面聚集时的生物运动与流体力学相互作用。具体考虑如下系统：

\begin{cases} \partial_t n - D\Delta n + u \cdot \nabla n + \chi \nabla \cdot \left(n \frac{\nabla c}{c}\right) = 0, & \Omega \times (0, T), \\ \partial_t c - \varepsilon \Delta c + u \cdot \nabla c + nc = 0, & \Omega \times (0, T), \\ \partial_t u - \varepsilon \Delta u + (u \cdot \nabla)u + \nabla p = n\nabla \phi, & \Omega \times (0, T), \\ \nabla \cdot u = 0, & \Omega \times (0, T). \end{cases}

其中 $$n$$ 为细胞密度， $$c$$ 为氧气浓度， $$u$$ 为流体速度， $$p$$ 为压力。 $$D>0$$ 和 $\varepsilon \ge 0$ 分别为细胞和氧气的扩散系数， $\chi > 0$ 为趋化系数。为简化起见，文中假设 $D=\chi=1$ 且 $\nabla \phi = e_2 = (0, 1)$ 。

该系统的核心难点在于趋化项中的对数奇异性（logarithmic sensitivity），即当 $c \to 0$ 时 $\nabla c / c$ 发散。为了处理这一奇异性并建立统一的分析框架，作者采用了 Cole-Hopf 型变换 $v_\varepsilon = -\nabla \ln c = -\nabla c / c$ 。这一变换将原奇异系统转化为一个非奇异的守恒律系统：

\begin{cases} \partial_t n_\varepsilon - \Delta n_\varepsilon + u_\varepsilon \cdot \nabla n_\varepsilon - \nabla \cdot (n_\varepsilon v_\varepsilon) = 0, \\ \partial_t v_\varepsilon - \varepsilon \Delta v_\varepsilon + \nabla(u_\varepsilon \cdot v_\varepsilon) + \nabla(\varepsilon|v_\varepsilon|^2 - n_\varepsilon) = 0, \\ \partial_t u_\varepsilon - \varepsilon \Delta u_\varepsilon + (u_\varepsilon \cdot \nabla)u_\varepsilon + \nabla p_\varepsilon = n_\varepsilon e_2, \\ \nabla \cdot u_\varepsilon = 0, \end{cases}

并施加相应的滑移边界条件。对于 $\varepsilon > 0$ ，边界条件为 $(\partial_y n_\varepsilon + n_\varepsilon v_{\varepsilon, 2})|_{y=0} = 0$ , $(v_{\varepsilon, 2}, u_{\varepsilon, 2})|_{y=0} = 0$ , $(\partial_y v_{\varepsilon, 1}, \partial_y u_{\varepsilon, 1})|_{y=0} = 0$ ；而对于 $\varepsilon = 0$ 的非粘性极限情况，边界条件变为 $(\partial_y n_0 + n_0 v_{0, 2})|_{y=0} = 0$ , $u_{0, 2}|_{y=0} = 0$ 。

物理背景上，边界层的产生源于好氧细菌在水面附近聚集引起的粘性效应，这种聚集进一步调节了流速和氧气溶解动力学。因此，理解 $\varepsilon \to 0$ 时解的行为至关重要。已有文献分别研究了 Keller-Segel 系统和 Navier-Stokes 系统的边界层问题，但针对耦合系统 (1.3) 在滑移边界条件下的粘性消失极限研究尚属空白。特别是，由于 $n_\varepsilon$ 的边界条件在 $\varepsilon > 0$ 和 $\varepsilon = 0$ 时保持一致，故 $n_\varepsilon$ 分量不产生内部边界层，而 $v_\varepsilon$ 和 $u_\varepsilon$ 则存在显著的边界层结构。本文旨在严格证明当 $\varepsilon \to 0$ 时，粘性解收敛于外层解（ $\varepsilon=0$ 的解）与内层边界层剖面的叠加，并给出明确的收敛速率。

2. 主要定理或主要结果

本文的主要结果概括为以下两个定理：

定理 1.1 (变换后系统的收敛性): 假设初始数据满足 $(n^{in}, v^{in}, u^{in}) \in H^{14}_{xy} \times H^{14}_{xy} \times H^{15}_{xy}$ 。设 $$(n_0, v_0, u_0)$$ 为第一部分工作中得到的 $\varepsilon=0$ 系统的解，且 $0 < T \le T_{max}$ 。则存在 $\varepsilon_0 > 0$ ，使得对任意 $0 < \varepsilon \le \varepsilon_0$ ，问题 (1.3)-(1.4) $$_1$$ 存在唯一解 $(n_\varepsilon, v_\varepsilon, u_\varepsilon)$ ，且在 $$[0, T]$$ 上满足如下 $L^\infty_{xy}$ 收敛估计：

\sup_{t \in [0, T]} \| n_\varepsilon(t, x, y) - n_0(t, x, y) \|_{L^\infty_{xy}} \le C \varepsilon^{1/4},

\sup_{t \in [0, T]} \| v_\varepsilon(t, x, y) - v_0(t, x, y) - (0, v^{B,0}_2)(t, x, \frac{y}{\sqrt{\varepsilon}}) \|_{L^\infty_{xy}} \le C \varepsilon^{1/4},

\sup_{t \in [0, T]} \| u_\varepsilon(t, x, y) - u_0(t, x, y) \|_{L^\infty_{xy}} \le C \varepsilon^{1/2},

其中 $$C$$ 为与 $\varepsilon$ 无关的常数， $(v^{B,0}_1, v^{B,0}_2)$ 为内层边界层剖面，且此处 $v^{B,0}_1 \equiv 0$ 。

定理 1.2 (原始系统的收敛性): 假设初始数据满足 $(n^{in}, \ln c^{in}, u^{in}) \in H^{14}_{xy} \times H^{15}_{xy} \times H^{15}_{xy}$ ，且 $n^{in} \ge 0, c^{in} > 0$ 。则原始趋化-纳维-斯托克斯系统 (1.10) 在 $$[0, T]$$ 上存在唯一解 $(n_\varepsilon, c_\varepsilon, u_\varepsilon)$ ，且满足：

\sup_{t \in [0, T]} \| c_\varepsilon(t, x, y) - c_0(t, x, y) \|_{L^2_{xy}} \le C \varepsilon^{1/2},

\sup_{t \in [0, T]} \| c_\varepsilon(t, x, y) - c_0(t, x, y) \|_{L^\infty_{xy}} \le C \varepsilon^{1/4}.

此外， $n_\varepsilon$ 和 $u_\varepsilon$ 的收敛性与定理 1.1 一致。特别地，边界层主要出现在 $c_\varepsilon$ 的导数中，即：

\sup_{t \in [0, T]} \| \nabla c_\varepsilon(t, x, y) - \nabla c_0(t, x, y) - (0, c_0(t, x, y)v^{B,0}_2)(t, x, \frac{y}{\sqrt{\varepsilon}}) \|_{L^\infty_{xy}} \le C \varepsilon^{1/4}.

这表明化学扩散率 $\varepsilon$ 在系统中扮演重要角色，不可忽略。

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

证明的核心障碍在于处理粘性耗散项与非粘性非线性项之间的不平衡，特别是在缺乏一般 Sobolev 空间均匀高阶正则性的情况下。作者的真正创新点在于引入共法线 Sobolev 空间 (Conormal Sobolev Spaces) 来建立粘性解的弱正则性估计，从而克服直接能量估计失效的问题。

Cole-Hopf 变换与非奇异化: 首先通过 $v_\varepsilon = -\nabla \ln c$ 消除对数奇异性，将问题转化为非奇异系统。这一步不仅简化了方程结构，还揭示了 $v_\varepsilon$ 必须是无旋的 ( $\nabla \times v_\varepsilon = 0$ )，这一内在约束在边界条件的匹配中起到了关键作用。
边界层展开与余项方程推导: 作者根据 Prandtl 边界层理论，将解展开为外层解、内层边界层剖面以及高阶余项之和。具体形式为：
$n_\varepsilon = n_0 + \sqrt{\varepsilon}n^{B,1} + \varepsilon n^{B,2} + \varepsilon^{1/2}H_\varepsilon,$
$v_\varepsilon = v_0 + (0, v^{B,0}_2) + \sqrt{\varepsilon}v^{B,1} + \varepsilon(v^{B,2}_1, 0) + \varepsilon^{1/2}V_\varepsilon,$
$u_\varepsilon = u_0 + \sqrt{\varepsilon}(u^{B,1}_1, 0) + \varepsilon u^{B,2} + \varepsilon^{3/2}(0, u^{B,3}_2) + \varepsilon^{1/2}U_\varepsilon.$
通过代入原方程，推导出余项 $(H_\varepsilon, V_\varepsilon, U_\varepsilon)$ 满足的线性化抛物型系统 (2.3)。这些余项方程包含复杂的非线性耦合项，如 $\nabla(U_\varepsilon \cdot V_\theta)$ 和 $\varepsilon^{1/2}\nabla(|V_\varepsilon|^2)$ 。
共法线 Sobolev 空间的构造与应用: 由于在一般 Sobolev 空间中无法获得粘性解关于 $\varepsilon$ 的一致高阶正则性（受限于 $\nabla(u_\varepsilon \cdot v_\varepsilon)$ 项），作者定义了基于加权导数 $\partial^\alpha = \partial_x^{\alpha_1} \psi(y)^{\alpha_2} \partial_y^{\alpha_2}$ 的共法线 Sobolev 空间 $\tilde{H}^m(\mathbb{R}^2_+)$ ，其中 $\psi(y)$ 是一个在边界附近线性增长、远处饱和的光滑函数。在此空间中，作者建立了粘性解的局部适定性和一致能量估计（定理 2.1）。这一框架允许作者在控制边界附近的高阶法向导数时，利用权重函数抵消部分奇异性，从而获得足以进行余项估计的正则性。
涡度估计与 Leray 投影: 在处理余项方程中的压力项时，作者利用了不可压缩条件 $\nabla \cdot U_\varepsilon = 0$ ，引入涡度 $\omega_\varepsilon = \partial_x u_{\varepsilon, 2} - \partial_y u_{\varepsilon, 1}$ 并通过最大原理控制其 $L^\infty$ 范数。同时，使用 Leray 正交投影处理压力项，避免了直接求解泊松方程带来的正则性损失。
各向异性 Sobolev 嵌入与截止函数: 在估计余项的法向导数时，作者使用了各向异性 Sobolev 嵌入不等式，并结合截止函数 $\phi$ 来处理边界项的高阶正则性要求。这些技术共同确保了余项在 $L^2_{xy}$ 和 $L^\infty_{xy}$ 范数下的收敛性估计得以成立。

4. 与之前工作的比较

本文工作相较于既有文献的推进主要体现在以下几个方面：

从独立系统到耦合系统: 既往工作主要集中在独立的 Keller-Segel 系统或 Navier-Stokes 系统的边界层问题上。例如，Li-Zhao [25] 验证了无流体影响下 Keller-Segel 系统的边界层行为；Hou-Wang-Zhao [14] 和 Hou-Liu-Wang-Wang [10] 在一维区间上证明了 Keller-Segel 系统的边界层收敛性及最优速率 $O(\varepsilon^{1/2})$ ；Hou-Wang [13] 在二维半平面上推导了 Keller-Segel 系统的外层和内层结构。然而，这些工作均未考虑流体速度 $$u$$ 对趋化过程的反馈作用。本文首次针对耦合的趋化-纳维-斯托克斯系统 (1.3) 在二维半平面上的边界层问题进行了严格分析，填补了这一空白。
边界条件的差异: 与 [11, 12, 23] 等研究不同，本文考虑的是滑移边界条件 (Slip Boundary Conditions)，而非 Dirichlet 或 Neumann 条件。滑移边界条件在物理上更符合自由表面或水-气界面的实际情况，但其数学处理更为复杂，因为边界条件本身依赖于未知量 $v_\varepsilon$ 和 $u_\varepsilon$ 的法向导数。本文成功处理了这种非齐次、非局部的边界条件对边界层结构的影响。
对数敏感性的处理: 与传统的 Keller-Segel 系统相比，本文处理的系统包含对数敏感性项 $\nabla \cdot (n \nabla c / c)$ 。这一项在 $$c=0$$ 处具有奇异性，使得标准的能量方法难以直接应用。作者通过 Cole-Hopf 变换将其转化为非奇异系统，并在共法线 Sobolev 空间中建立了统一的估计框架。这一方法不同于以往仅针对光滑敏感性或线性敏感性的研究，为处理奇异趋化项提供了新的分析工具。
收敛速率与正则性代价: 在 Navier-Stokes 系统的边界层研究中，Iftimie-Sueur [16] 和 Wang-Xin-Zang [43] 分别获得了 $$L^2$$ 和 $L^\infty$ 收敛性，但通常假设固定的滑移长度或特定的单调性条件。本文在更一般的滑移条件下，获得了 $n_\varepsilon$ 和 $v_\varepsilon$ 的 $L^\infty$ 收敛速率 $O(\varepsilon^{1/4})$ 以及 $u_\varepsilon$ 的 $O(\varepsilon^{1/2})$ 。虽然 $n_\varepsilon$ 和 $v_\varepsilon$ 的 $L^\infty$ 速率低于 $$L^2$$ 速率 $O(\varepsilon^{1/2})$ ，但这主要是由于在共法线空间中只能获得较弱的正则性所致。这一结果揭示了在处理强非线性耦合和奇异项时，高阶正则性的丧失是不可避免的代价。
与第一部分工作的联系: 本文是两部分工作的第二部分。第一部分 [41] 主要关注 $\varepsilon=0$ 时系统的适定性以及边界层剖面的显式推导。本文则在此基础上，严格证明了 $\varepsilon > 0$ 时解向 $\varepsilon=0$ 解及边界层剖面的收敛性，并回答了第一部分中提出的关于耦合系统适定性的开放问题。附录中提供的共法线 Sobolev 空间中的适定性结果，进一步完善了该理论体系。

8. $$ rac{1}{2}$$ 阶 Euler 方程在 $$ rac{3}{4}$$ 空间上的稳定奇点

Stable Singularity of the Euler Equations on $\mathbb{R}^3$

作者: Adarsh Ganeshram, Valentin Duruisseaux, Anima Anandkumar

类别: math.AP

发布日期: 2026-09-09

arXiv ID: 2609.10867v1

摘要翻译: 我们提供了关于无界域上三维 Euler 方程存在稳定的有限时间奇点的证据。使用带有自相似假设的物理信息神经网络（PINN），我们在临界爆破速率 $$0.5$$ 处找到了 Euler 系统的一个近似奇异剖面，并使用样条表示对其进行了验证。与所获得剖面相关的输运场表明，可以在整个域内建立线性阻尼，这为候选剖面的整体稳定性提供了强有力的证据。我们还建立了证明该近似自相似剖面非线性稳定性的完整框架，将分析简化为一组庞大但有限的显式估计和可计算常数。

Comment: 107 pages

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

1. 研究对象和背景

本文研究的核心问题是三维不可压缩欧拉方程在无限域 $\mathbb{R}^3$ 上是否允许从光滑初始数据出发，在有限时间内形成奇点（即解失去正则性）。该问题与著名的 Navier-Stokes 方程的千禧年问题密切相关。欧拉方程描述了无粘流体运动，其形式为：

\partial_t \tilde{u} + (\tilde{u} \cdot \nabla) \tilde{u} = -\nabla \tilde{p}, \quad \nabla \cdot \tilde{u} = 0,

其中 $\tilde{u}$ 是速度场， $\tilde{p}$ 是压力。通过取旋度，可以得到涡量方程：

\partial_t \tilde{\omega} + \tilde{u} \cdot \nabla \tilde{\omega} = \tilde{\omega} \cdot \nabla \tilde{u}, \quad \tilde{\omega} := \nabla \times \tilde{u}.

尽管在有界域或特定非光滑初始条件下已证明存在有限时间爆破，但在无限域 $\mathbb{R}^3$ 上，对于光滑初始数据的有限时间奇点形成仍是开放问题。一些非存在性结果（如 Chae [11, 12]）表明，如果涡量满足超代数衰减条件，则不存在有限时间爆破。然而，这些假设在实际物理场景中可能过于严格。

为了探索可能的奇点结构，研究者通常采用自相似假设（self-similar ansatz）。具体而言，假设解在接近爆破时间 $$T$$ 时具有如下形式：

u_\bullet(\tilde{r}, \tilde{z}, t) = (T-t)^{c_u} U\left(\frac{\tilde{r}}{(T-t)^\lambda}, \frac{\tilde{z}-\tilde{z}_c(t)}{(T-t)^\lambda}\right),

\omega_\bullet(\tilde{r}, \tilde{z}, t) = (T-t)^{c_\omega} \Omega\left(\frac{\tilde{r}}{(T-t)^\lambda}, \frac{\tilde{z}-\tilde{z}_c(t)}{(T-t)^\lambda}\right),

其中 $\lambda$ 是空间缩放指数， $c_u, c_\omega$ 是振幅指数， $\tilde{z}_c(t)$ 是奇点中心的轴向位置。根据 Constantin et al. [17] 的理论，临界缩放指数 $\lambda = 0.5$ 是数学上最有可能出现稳定自相似奇点的“区分”区域。

本文旨在通过数值方法发现一个在 $\lambda = 0.5$ 下的近似自相似奇点剖面，并建立其线性稳定性框架，从而为证明非线性稳定性提供基础。这一工作填补了现有文献中缺乏对无限域上欧拉方程稳定自相似奇点的定量证据的空白。

2. 主要定理或主要结果

近似奇点剖面的发现: 作者使用物理信息神经网络（PINN）结合行波自相似假设，在轴对称设定下发现了三维欧拉方程的一个高精度近似自相似奇点剖面 $(U, \Omega, \Psi)$ 。该剖面满足以下稳态剖面方程：
$U + (\lambda r + U_r) \partial_r U + (C + \lambda z + U_z) \partial_z U = 2U \partial_z \Psi,$
$(1+\lambda)\Omega + (\lambda r + U_r) \partial_r \Omega + (C + \lambda z + U_z) \partial_z \Omega = 2U \partial_z U,$
$-E \Psi = \Omega,$
其中 $E = \partial_r^2 + \frac{3}{r}\partial_r + \partial_z^2$ 是椭圆算子， $\lambda = 0.5$ 是临界缩放指数， $$C$$ 是行波速度。通过优化，作者将 PDE 残差的均方误差（MSE）降低至 $10^{-10}$ 量级，并通过样条函数表示进一步认证了剖面的准确性，使得 $$L^2$$ 和 $L^\infty$ 范数下的残差分别控制在 $O(10^{-5})$ 和 $O(10^{-3})$ 。
线性阻尼的建立: 作者证明了在该近似剖面附近，线性化系统在整个计算域内具有线性阻尼特性。具体而言，通过优化加权能量泛函中的权重函数，作者在大部分区域获得了低阶阻尼，而在两个关键的低输运区域（轴上固定点和离轴固定点附近）分别通过奇异权重和高阶导数的有利符号实现了局部阻尼。这一结果表明，候选剖面在微扰下是稳定的，为后续的非线性稳定性分析提供了关键的量化依据。
非线性稳定性框架的建立: 作者构建了一个完整的非线性稳定性证明框架，将稳定性问题归结为一组有限但显式的估计和可计算常数的认证。这些估计包括低阶和高阶能量估计，以及线性项和非线性项的控制。通过区间算术（interval arithmetic）和 Arb 库，作者对这些常数进行了严格认证，确保所有估计在数值上是可靠的。虽然完整的非线性稳定性证明仍在进行中，但该框架已经为最终证明奠定了坚实基础。

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

物理信息神经网络（PINN）与优化策略: 作者采用了先进的 PINN 方法来发现近似奇点剖面。为了克服传统 PINN 训练中的困难（如收敛到平凡解、损失景观复杂等），作者引入了一系列创新技术：
- 自适应损失权重调度：动态调整不同损失项的相对重要性，以平衡竞争目标。
- FP64 精度计算：使用 64 位浮点精度减少舍入误差，提高梯度质量和收敛性。
- PINN Boosting：通过逐步添加新网络来纠正前一阶段的残差误差，将一个全局优化问题分解为一系列更简单的修正问题。
- 曲率感知优化器：使用 SS-eSOAP 和 SS-Broyden 优化器，利用二阶曲率信息改善条件数和收敛速度，将 MSE 从 $10^{-4}$ 降低至 $10^{-10}$ 。
- Sinh 变换与自适应配点：通过 sinh 变换高效采样无限域，并根据残差大小动态重新采样配点，集中计算资源于最难解析的区域。
样条函数表示与区间算术认证: 为了将数值结果转化为严格的数学证明，作者将 PINN 输出的剖面转换为分段多项式样条函数。样条函数的局部低阶结构支持精确的微分和残差估计，避免了直接认证神经网络输出的复杂性。作者使用 Arb 库进行任意精度区间算术，严格控制了所有数值误差，确保了剖面依赖常数的可靠性。
加权能量方法与参数优化: 在稳定性分析中，作者采用了加权 $$L^2$$ 和 Sobolev 范数来定义能量泛函。与传统方法不同，作者没有手动选择权重函数，而是将其参数化为高斯基函数、奇异分量、远场神经网络和全局修正的组合，并通过数值优化确定最优参数。这种自动化方法显著简化了稳定性证明中的试错过程，并确保阻尼机制在所有关键区域都得到充分控制。
动态重标度与调制参数: 作者引入了动态重标度框架，将固定的自相似参数替换为随时间变化的调制参数（如缩放速率、平移速率和振幅归一化）。这使得可以在归一化变量中研究长时演化，避免了物理变量中振幅发散和长度尺度坍缩带来的数值困难。通过线性化动态重标度方程，作者能够系统地分析扰动行为，并将稳定性问题转化为对线性化和非线性项的精细估计。

4. 与之前工作的比较

与 Constantin et al. [17] 的比较: Constantin et al. 的理论工作指出，临界缩放指数 $\lambda = 0.5$ 是欧拉方程自相似爆破的最可能区域，并排除了 $\lambda < 0.5$ 和 $\lambda > 0.5$ 的稳定自相似奇点。本文通过数值发现了一个在 $\lambda = 0.5$ 下的高精度近似剖面，验证了这一理论预测。此外，Constantin et al. 强调了局部外向性质（local outgoing property）的重要性，而本文通过优化权重函数，成功地在轴上和离轴固定点附近建立了所需的阻尼机制，弥补了理论分析与数值实现之间的差距。
与 Hou et al. [20] 的比较: Hou et al. 研究了带有对流参数的欧拉方程族，并在弱对流情形（ $\epsilon \leq 0.3$ ）下观察到了行波自相似奇点，但在强对流情形（ $\epsilon = 0.4$ ）下未能观察到类似行为。本文通过 PINN 方法在整个对流强度范围 $\epsilon \in [0, 1]$ 内找到了高精度近似剖面，表明先前动态演化中消失的场景并不意味奇点剖面的不存在，而是由于前向动力学未能达到该状态。这一结果扩展了已知奇点结构的适用范围，并为理解强对流下的奇点形成提供了新的视角。
与 Chen et al. [3, 10, 18] 的比较: Chen et al. 在有界域或特定假设下建立了欧拉方程爆破的稳定性理论。本文借鉴了其稳定性框架，但针对无限域 $\mathbb{R}^3$ 和无边界条件进行了重大改进。特别是，作者通过数值优化权重函数和其他可调参数，避免了传统方法中繁琐的手动试错过程，提高了证明的可操作性和可扩展性。此外，本文首次将 PINN 与区间算术结合，实现了从数值发现到严格认证的无缝衔接，为计算机辅助证明在 PDE 中的应用树立了新范例。
与非存在性结果的比较: Chae [11, 12] 的非存在性定理要求涡量具有超代数衰减，而本文发现的剖面仅表现出多项式衰减，因此不违反这些定理。这表明，有限时间奇点可能在更广泛的初值空间中存在，只要它们不满足过强的衰减条件。这一发现挑战了以往认为无限域上不可能存在有限时间爆破的观点，并为未来放宽非存在性结果的假设提供了方向。

9. 三维球对称可压缩 Navier-Stokes 方程中任意大初值强解全局存在与爆破的尖锐符号判据

Sharp sign criterion for global existence and blow-up of strong solutions with arbitrarily large initial data in the 3D spherically symmetric compressible Navier-Stokes equations

作者: Xiangdi Huang, Bingshu Li

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11304v1

摘要翻译: 我们建立了三维球对称可压缩 Navier-Stokes 方程在实心球域上强解全局存在与有限时间爆破的尖锐判据，该判据仅依赖于有效速度的初始符号。粘性系数假设满足 Bresch-Desjardins 结构 $$μ=ρ^α$$ , $$λ=(α-1)ρ^α$$ ，其中有效速度由 $v=u+αρ^{α-2}ρ_r$ 给出。此前，高维情形下强解的全局存在性结果仅限于 $α\le 1$ 的情形，且端点 $$α=1$$ 对应粘性 Saint-Venant（浅水）系统。本文将该全局存在性理论拓展至超临界区域 $$α>1$$ 。具体而言，只要初始有效速度在球面上非负，我们便证明了任意大初值下强解的全局存在性。主要难点在于推导密度的均匀下界，这需要处理球心处的奇异估计。为此，我们在球心附近利用若干新颖的量，并采用改进的最大值原理，首次建立了球心处速度的 Lipschitz 连续性。随后，我们对适当构造的函数推导了 Dini-Gronwall 型不等式，从而闭合了密度下界的估计。相反地，我们可以构造一族初值，使得初始有效速度在球面上为负，从而对应的强解在有限时间内发生爆破，且在爆破时刻真空恰好出现在球面上。我们的结果建立了一个清晰的两分法：全局正则性与奇异性形成之间的界限完全由有效速度的初始符号决定。

Comment: 50 pages

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

1. 研究对象和背景

本文研究的是在三维实心球 $B_R \subset \mathbb{R}^3$ 上, 具有依赖于密度的粘滞系数的等熵可压缩 Navier-Stokes 方程。具体而言, 考虑如下方程组:

\begin{cases} \partial_t \rho + \text{div}(\rho u) = 0, \\ \partial_t(\rho u) + \text{div}(\rho u \otimes u) + \nabla P(\rho) - \text{div}(\mu(\rho) Du) - \nabla(\lambda(\rho) \text{div } u) = 0, \end{cases}

其中 $Du = \frac{1}{2}(\nabla u + (\nabla u)^T)$ , 状态方程为 $P(\rho) = \rho^\gamma$ , 粘滞系数满足 Bresch-Desjardins (BD) 结构:

\mu(\rho) = \rho^\alpha, \quad \lambda(\rho) = (\alpha - 1)\rho^\alpha.

边界条件为无滑移条件 $$u(t, x) = 0$$ on $\partial B_R$ . 我们关注球对称解, 即 $\rho(t, x) = \rho(t, r)$ , $u(t, x) = u(t, r)\frac{x}{r}$ , 其中 $r = |x| \in [0, R]$ . 此时方程化为径向形式:

\begin{cases} \rho_t + (\rho u)_r + \frac{2}{r}\rho u = 0, \\ (\rho u)_t + (\rho u^2)_r + \frac{2}{r}\rho u^2 + (\rho^\gamma)_r - \left[\frac{\alpha}{r^2}\rho^\alpha(r^2 u)_r\right]_r + \frac{2}{r}(\rho^\alpha)_r u = 0. \end{cases}

真空形成是可压缩粘性流动经典可解性的主要障碍之一. 当粘滞系数随密度消失时, $\rho \to 0$ 的过程同时导致抛物型平滑作用的退化以及对速度控制能力的减弱. 核心问题在于: 在超临界区域 $\alpha > 1$ , 初始光滑且严格正的密度是否能保持远离零的下界, 或者是否会在球内某处形成真空?

经典的常数粘滞理论提供了广泛背景. Matsumura 和 Nishida 建立了小扰动下的全局光滑解 [21], 而 Lions [20] 以及 Feireisl-Novotný-Petzeltová [10] 发展了大初值及可能含真空的全局有限能量弱解理论. 对于依赖密度的粘滞系数, 关键的结构进展是 Bresch 和 Desjardins [1] 发现的额外熵, 并与 Lin [2] 共同发展. 相关的修正速度将物理速度与密度梯度耦合, 从而获得了泛型粘滞律下不可得的估计.

这种 BD 结构支撑了大量弱解理论. 例如, Mellet 和 Vasseur 获得了稳定性估计 [22], Guo-Jiu-Xin 构造了全局球对称弱解 [11], Li-Xin [19] 和 Vasseur-Yu [24] 处理了多维退化系统, 而 Bresch-Vasseur-Yu 后来涵盖了物理对称应力及广泛的非线性粘滞律 [3]. 特别是, 后一理论包含了 $\alpha = 1$ 两侧的幂律. 因此, 全局弱解的存在性并非本文所解决的问题: 弱解可能包含真空, 并不能决定一个初始为正的经典分支是否保持与零均匀分离.

一维结果揭示了有效变量在强解理论中的互补作用. Haspot 在亚临界范围获得了全局强解 [14], 而 Constantin-Drivas-Nguyen-Pasqualotto [8] 以及 Burtea-Haspot [4] 在关于某个有效量的一边或 Oleinik 型假设下, 将结果扩展到 $\alpha > 1$ . 然而, 一维问题既没有球心处的坐标奇异性, 也没有此处存在的有限半径边界约束.

球对称性将系统简化为一个径向变量, 但并未将其转化为非奇异的一维问题. 权重 $$r^2$$ , 商 $$u/r$$ , 以及差 $$u_r - u/r$$ 使得原点成为一个真正的分析难点. 早期的球对称真空结果通常将异常行为定位在此中心. Xin-Yuan 证明在其非相互作用假设下, 真空不会在原点之外形成, 并给出了一个关于速度正则性的独立判据以排除真空形成 [25]; Guo-Wang-Wang 获得了大初值的球对称理论及另一种粘滞律, 其中正性在 $$r=0$$ 附近受到保护 [12]. 相比之下, 本文构造的密度损失发生在外边界, 且仅发生于此.

最近的相关结果是 Zhang 在同一 BD 问题、固定球域及无滑移条件下的大初值经典理论. Zhang 证明了任意严格正初值在以下范围内的全局径向强解:

N=2: \frac{4}{5} \le \alpha < 1, \gamma > 1; \quad N=3: \frac{7}{8} \le \alpha < 1, 1 < \gamma < \frac{9\alpha-6}{3}.

Huang-Meng-Zhang [16] 将此扩大至

N=2: \frac{1}{2} < \alpha < 1, \gamma > 1;

以及, 利用其论文中定义的 $n_3(\alpha)$ , 至

N=3: 0.686 < \alpha < 1, 1 < \gamma < \frac{6\alpha-3+3-5\alpha}{2n_3(\alpha)}.

他们还处理了 $\alpha=1$ 时的二维端点情形, 并在部分 $\alpha \ge 1$ 范围内构造了允许真空的全局径向弱解. Lei 随后获得了经典解, 并证明了在以下范围内真空不会在有限时间内形成:

N=2: 0.54369 < \alpha < 1, \gamma > 1; \quad N=3: 0.67661 < \alpha < 1, 1 < \gamma < 6\alpha-3.

Chen-Zhang-Zhu 证明了远场含真空的大径向正则解, 范围为 $1/2 < \alpha < 1$ [6], 而 X. Zhang 研究了三维等温外问题, 远场含真空 [29]. 这些域和远场状态不同于此处考虑的固定边界正问题, 但这些结果清晰地展示了现有无限制径向经典理论中的阈值.

此外, 还存在其他结构下 $\alpha > 1$ 的多维强解结果. Yu [27] 和 Huang-Li-Zhang [15] 处理了围绕足够大的正密度水平的三维柯西问题; Xin-Zhu [26] 研究了无穷远处含真空的一类特殊膨胀速度; Du-Guo-Xu [9] 证明了有界域上带有滑移边界条件、大外部势及足够大的密度下界的结果. 移动界面和自由边界的构造构成了另一个邻近理论: Guo-Xu-Zhang [13] 在移动无应力界面上拼接内部正解和外部解, 而 Chen-Zhang-Zhu [7] 处理了基于距离的更广泛的真空类, 包括物理真空自由边界.

最后, Merle-Raphaël-Rodnianski-Szeftel [23] 的坍缩构造及其由 Chen-Liu-Zhu [5] 给出的退化粘滞对应物涉及原点处密度的爆破, 而非边界上真空的形成.

接下来解释为何 $\alpha = 1$ 是排除真空的结构阈值. 定义 BD 势 $\Phi_\alpha$ 为

\Phi'_\alpha(\rho) = \alpha \rho^{\alpha-2}, \quad \Phi_\alpha(\rho) = \begin{cases} \frac{\alpha}{\alpha-1}(\rho^{\alpha-1}-1), & \alpha \neq 1, \\ \log \rho, & \alpha = 1. \end{cases}

径向有效速度满足

v - u = \partial_r \Phi_\alpha(\rho).

因此, 对物理速度和有效速度的适当估计控制了变换后密度变量的空间变化.

对于 $\alpha \le 1$ , 有

\Phi_\alpha(\rho) \to -\infty \quad \text{当 } \rho \downarrow 0.

在此变换变量中, 真空因此无限远离任何正密度值. 由于正质量守恒保证密度在球内某处至少达到其正空间平均值, 对该奇异变量空间变化的足够强控制可以防止密度在有界域上达到零. 在现有的 $\alpha < 1$ 径向论证中, 这一原则通过对密度的邻近奇异负幂次进行估计来实现; 在二维端点 $\alpha = 1$ 处, 则通过对 $\partial_r \log \rho = v-u$ 的控制来实现; 参见 [16, Sections 3.3–3.4 and 4.2]. 而对于 $\alpha > 1$ , $\Phi_\alpha$ 在 $\rho = 0$ 处保持有限, 因此用于下界的此机制消失.

同一阈值也反映在单位质量的粘性扩散尺度上:

\frac{\mu(\rho)}{\rho} = \rho^{\alpha-1}.

它在 $\alpha < 1$ 时靠近真空更强, 在 $\alpha = 1$ 时为一阶, 而在 $\alpha > 1$ 时消失. 因此, 在超临界区域, 奇异密度壁垒恰好在粘性正则化最弱的地方丢失.

边界迹解释了质量坐标中相关的是哪种剖面. 对于每个正的径向经典解, 下面导出的有效方程为

Dt v + F_1(\rho)(v-u) = 0, \quad F_1(\rho) = \frac{\gamma}{\alpha} \rho^{\gamma-\alpha}.

由于 $$u(t,R)=0$$ , 固定边界是一条物质轨迹, 且

\partial_t v(t,R) + F_1(\rho(t,R))v(t,R) = 0, \quad v(t,R) = v_0(R) \exp\left(-\int_0^t F_1(\rho(s,R)) ds\right).

因此, $$v(t,R)$$ 的符号在整个正经典寿命期间保持不变.

为了理解相应的密度几何, 引入

\delta(t,r) := \int_r^R \rho(t,s)s^2 ds.

因此 $\delta$ 是 $$r$$ 和边界之间包含的质量; 质量坐标的系统表述将在下文发展. 由于 $\delta_r = -\rho r^2$ , 有效速度的恒等式给出

(\rho^\alpha)_\delta = -\frac{v-u}{r^2}.

在无滑移边界处, 这变为

(\rho^\alpha)_\delta(t,0) = -\frac{v(t,R)}{R^2}.

因此, 有限的非零负边界迹 $$v(t,R)<0$$ 选择了有限的正向内质量斜率:

\rho^\alpha(t,\delta) = \rho(t,R)^\alpha + A(t)\delta + o(\delta), \quad A(t) = -\frac{v(t,R)}{R^2} > 0.

如果边界密度趋于零, 而该斜率保持非退化, 则候选终端切线为

\rho^\alpha \sim A\delta.

剩下的问题是, 这样的质量剖面能否占据一个具有有限物理厚度的衣领. 一个质量切片 $d\delta$ 具有径向厚度

d(R-r) = \frac{d\delta}{\rho r^2}.

对于 $\rho^\alpha \sim A\delta$ , 有 $\rho \sim (A\delta)^{1/\alpha}$ , 因此, 在有限边界附近 $r \sim R$ ,

R-r \sim \frac{1}{R^2 A^{1/\alpha}} \int_0^\delta s^{-1/\alpha} ds.

当 $\alpha > 1$ 时,

\int_0^\delta s^{-1/\alpha} ds = \frac{\alpha}{\alpha-1} \delta^{(\alpha-1)/\alpha} < \infty.

等价地, 在物理空间中对应的行为是

\rho(r) \sim C(R-r)^{1/(\alpha-1)}.

对于 $\alpha \le 1$ , 同样的积分发散, 因此仿射 $\rho^\alpha$ 剖面无法将有界半径边界上的正内部密度连接到真空. 这正是后续章节中构造所利用的有限容量机制.

我们的第一个定理使用了有利的符号. 对于 $\alpha > 1$ 和 $\gamma > \alpha$ , 如果一个正的径向经典解在有限时间区间上有界的密度上限, 且 $v_0(R) \ge 0$ , 那么密度在整个球上在该区间内远离零. 证明结合了边界条件、与中心兼容的辅助量、倒数密度变量的改进最大原理以及 Dini-Gronwall 估计. 在下面指定的指数范围内, 可以在每个有限时间区间上获得密度的上界. 然后, 条件性下界为该范围内具有非负边界迹的每个相容数据提供全局正径向经典解.

我们的构造定理实现了相反的符号. 在相同的指数范围 (7) 中, 我们构造光滑相容径向数据, 其密度严格为正且 $v_{0,k}(R) < 0$ , 其最大正经典寿命满足 $T_k^* \le C k^{(\alpha-1)/\alpha}$ . 密度在终端时间一致收敛到一个非平凡剖面, 其零集恰好为 $\{R\}$ . 因此, 崩溃是固定边界上的有限时间空化.

一起, 这些结果将初始边界有效速度识别为超临界径向 BD 问题的边界序参量. 它们展示了由同一迹选择的两条严格相反的边界密度机制, 其符号被保留. 注意, 负边界迹的结果是一个显式构造, 而不是对所有具有 $$v_0(R) < 0$$ 的数据的分类.

2. 主要定理或主要结果

本文的主要结果包括两个核心定理和一个推论, 分别对应于非负边界迹的全局存在性和负边界迹的有限时间爆破.

定理 2.1 (非负边界迹的密度下界) 设 $\alpha > 1$ 且 $\gamma > \alpha$ . 设 $(\rho, u)$ 是方程 (1) 在 $[0, T) \times B_R$ 上的最大正径向经典解, 其中 $T \in (0, \infty)$ 是其作为正径向经典解的最大存在时间. 假设

\sup_{0 \le t < T} \|\rho(t)\|_{L^\infty(B_R)} =: \rho_T < \infty,

并且初始有效速度具有非负边界迹:

v_0(R) \ge 0.

那么存在常数 $$c_T > 0$$ , 仅依赖于 $R, \alpha, \gamma, T, \rho_T$ 和初始数据, 使得

\rho(t, r) \ge c_T \quad (0 \le t < T, 0 \le r \le R).

备注 2.2: 定理 2.1 的结论是一个当 $$t$$ 趋近于最大存在时间 $$T$$ 时仍然均匀的下界. 这样的界不能仅仅从经典正性得出: $$[0, T)$$ 上的正性只保证在每个紧时间区间 $[0, \tau]$ ( $\tau < T$ ) 上有正的下界.

推论 2.3 (非负边界迹的全局解) 假设 (7) 成立, 即

1 < \alpha < 11.7, \quad 2\alpha - 1 \le \gamma < 3\alpha - 1.

设 $(\rho_0, u_0)$ 是满足命题 3.3 假设的相容径向初始数据, $\rho_0$ 严格为正, 并且假设 $v_0(R) \ge 0$ . 那么相应的最大正径向经典解是全局的:

T^* = \infty.

此外, 对于每个有限的 $$T_0 > 0$$ , 存在常数 $0 < c_{T_0} \le C_{T_0} < \infty$ 使得

c_{T_0} \le \rho(t, r) \le C_{T_0} \quad (0 \le t \le T_0, 0 \le r \le R).

备注 2.4: 该定理的意义在于确立了以下事实: 对于超临界情况 $\alpha > 1$ , 如果初始有效速度在边界上非负, 则三维可压缩 Navier-Stokes 方程存在唯一的全局强解. 这是多维设置中任意大初值强解全局存在的首个结果. 与 Haspot 的一维结果 [4] 相比, 高维系统引入了额外的几何奇异项, 带来了实质性困难. 主要困难在于建立密度的均匀下界, 这需要处理球心处速度的奇异估计——一种在一维设置中不出现的奇异性. 为了克服这一点, 我们利用了中心附近的几个新量, 并采用改进的最大原理首次建立了速度在中心的 Lipschitz 连续性, 从而解决了速度梯度的奇异性. 随后, 我们为适当构造的函数推导了 Dini-Gronwall 型不等式, 从而封闭了密度下界估计. 这一论证关键依赖于初始有效速度在边界上的非负性.

定理 2.5 (内部密度保持与真空局部化) 设 $\alpha > 1$ , $\gamma > \alpha$ , 且 $0 < T < \infty$ . 设 $(\rho, u)$ 是 $[0, T) \times B_R$ 上的正径向经典解, 满足密度上界 (8). 设 $$m$$ 和 $$M$$ 分别是径向质量坐标和守恒总质量, 定义如 (13). 不对 $$v_0(R)$$ 施加任何符号条件. 那么存在有限常数 $$C_X$$ , 仅依赖于 $R, \alpha, \gamma, T, \rho_T$ 和初始数据, 使得

\frac{(M-m)^{1/\alpha}}{\rho(t,m)} \le C_X \quad (0 \le t < T, 0 \le m \le M).

因此, 对于每个 $\varepsilon \in (0, R)$ , 存在 $c_{T,\varepsilon} > 0$ , 具有相同的常数依赖性加上 $\varepsilon$ , 使得

\rho(t, r) \ge c_{T,\varepsilon} \quad (0 \le t < T, 0 \le r \le R-\varepsilon).

特别地, 如果 $t_j \uparrow T$ 且 $\rho(t_j, r_j) \to 0$ , 则 $r_j \to R$ . 此外, 如果 $\rho(t, \cdot) \to \rho_*$ 在 $$[0, R]$$ 上一致收敛当 $t \uparrow T$ , 且 $\rho_*$ 有零点, 则

\{r \in [0, R] : \rho_*(r) = 0\} = \{R\}.

备注 2.6: 在定理 2.5 的假设下, 密度的均匀上界阻止了球内部（包括中心）的有限时间真空形成. 事实上, 无论初始边界有效速度 $$v_0(R)$$ 的符号如何, 密度在每个固定的较小球 $B_{R-\varepsilon}$ 上始终远离零. 因此, 任何沿其密度当 $t \uparrow T$ 时趋于零的点序列必须接近外边界 $\partial B_R$ . 如果存在均匀的终端密度剖面并包含真空, 其真空集恰好为 $\partial B_R$ . 该定理因此确定了有限时间真空形成的唯一可能位置.

定理 2.7 (负边界迹数据的有限时间边界真空) 假设 (7). 那么存在 $$k_0, C > 0$$ 使得, 对于每个 $$0 < k < k_0$$ , 可以选择光滑相容径向初始数据, 其密度严格为正且

v_{0,k}(R) < 0.

相应的最大正经典解存在于 $$[0, T_k^*)$$ , 其中

0 < T_k^* \le C k^{(\alpha-1)/\alpha} < \infty.

此外, 存在非负、非平凡的终端剖面 $\rho_{k,*}$ 使得

\rho_k(t) \longrightarrow \rho_{k,*} \quad \text{uniformly on } B_R \quad \text{as } t \uparrow T_k^*,

并且

\{r \in [0, R] : \rho_{k,*}(r) = 0\} = \{R\}.

注意, 我们不在 $$k$$ 上取极限, $$k$$ 只是一个指示我们有一族初始数据的索引.

备注 2.8 (符号机制的范围): 定理 2.1 的条件是在所考虑的有限时间区间上密度的上界, 其常数 $$c_T$$ 不一定在 $T \to \infty$ 时均匀. 推论 2.3 成为无条件的是因为该上界在 (7) 中可用. 相反, 定理 2.7 构造了一个具有负边界迹的特定族; 它不断言每个具有 $$v_0(R) < 0$$ 的数据都会形成真空.

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

本文的证明架构分为两部分, 分别处理“好”边界情况（非负边界迹）和“坏”边界情况（负边界迹）, 两者使用相同的边界迹但结构不同.

(i) 好边界情况：非负边界迹下的密度下界

主要技术障碍在于建立密度的均匀下界, 特别是在球心 $$r=0$$ 附近. 球对称性引入了权重 $$r^2$$ , 商 $$u/r$$ , 以及差 $$u_r - u/r$$ , 使得原点成为一个真正的分析难点. 在超临界区域 $\alpha > 1$ , BD 势 $\Phi_\alpha(\rho)$ 在 $\rho=0$ 处保持有限, 因此无法像 $\alpha \le 1$ 那样通过变换变量的奇异壁垒来直接获得密度下界. 此外, 单位质量的粘性扩散尺度 $\mu(\rho)/\rho = \rho^{\alpha-1}$ 在 $\alpha > 1$ 时随 $\rho \to 0$ 而消失, 意味着粘性正则化在最需要它的地方最弱.

作者的真正创新点在于引入了一组新的辅助量来处理中心奇异性, 并利用改进的最大原理和 Dini-Gronwall 估计来封闭密度下界.

横向应变率函数的引入与估计: 由于光滑径向对称性强制 $$u(t,0)=v(t,0)=0$$ , 中心值本身不包含关于速度第一非平凡行为的信息. 其径向展开为
$u(t,r) = u_r(t,0)r + O(r^3), \quad v(t,r) = v_r(t,0)r + O(r^3).$
因此, 相关的中心量是斜率 $$u_r(t,0)$$ 和 $$v_r(t,0)$$ . 作者将这些中心斜率扩展到整个球, 将 $$u/r$$ 和 $$v/r$$ 视为新的未知数. 定义横向应变率函数:
$U = \frac{u}{r}, \quad V = \frac{v}{r}.$
它们在 $$r=0$$ 处的值通过光滑径向极限理解. 这些商在 $$r=0$$ 处光滑, 并保留了在 $$u$$ 和 $$v$$ 的中心恒为零的值中丢失的信息. 还定义散度缩写:
$d := \text{div } u = u_r + 2\frac{u}{r}, \quad \omega := \text{div } v = v_r + 2\frac{v}{r}.$
这些横向应变率函数满足如下方程:
$Dt U = \alpha \rho^{\alpha-1} \left(U_{rr} + \frac{4}{r}U_r\right) + (V-U)[\alpha r (U)_r + (3\alpha-2)U - F_1(\rho)] - U^2,$
$Dt V = -UV - F_1(\rho)(V-U).$
命题 5.1 对这些横向应变率函数进行了估计, 控制了有效压力.
有效压力的控制: 命题 5.2 利用横向应变率函数的估计来控制有效压力. 有效压力与密度和速度相关, 其控制是下一步的关键.
倒数密度的最大原理与 Dini-Gronwall 估计: 命题 5.3 通过 Dini-Gronwall 论证封闭了倒数密度的最大不等式. 具体而言, 作者构造了适当的函数, 并推导了 Dini-Gronwall 型不等式, 从而封闭了密度下界估计. 这一论证关键依赖于初始有效速度在边界上的非负性. 边界条件 $$u(t,R)=0$$ 使得 $$v(t,R)$$ 的符号在整个正经典寿命期间保持不变, 从而为非负边界迹的情况提供了稳定的边界行为.
速度在中心的 Lipschitz 连续性: 作者首次建立了速度在中心的 Lipschitz 连续性, 从而解决了速度梯度的奇异性. 这是通过改进的最大原理实现的, 利用了中心附近的新量. 这一结果对于处理球心处的奇异估计至关重要.

(ii) 坏边界情况：负边界迹下的有限时间边界真空

主要技术障碍在于构造一族具有负边界迹的初始数据, 并证明其在有限时间内在边界上形成真空. 作者利用了“有限容量”机制.

密度上界与密度损失判据: 命题 4.1 提供了密度的上界和密度损失的判据. 这对于控制解的行为至关重要.
相容数据的构造与负通量平台: 引理 6.3 构造了相容数据, 创建了一个负通量平台. 这意味着在边界附近, 有效速度 $$v$$ 保持负值, 导致质量向内流动.
观察衣领上的缩放最大原理与通量持久性: 命题 6.4 在一个观察衣领上保持了足够的通量. 作者使用缩放最大原理和论证来确保负通量在一段时间内持续存在.
有限容量导致的矛盾: 命题 6.6 利用了有限容量机制. 当 $\alpha > 1$ 时, 积分
$\int_0^\delta s^{-1/\alpha} ds = \frac{\alpha}{\alpha-1} \delta^{(\alpha-1)/\alpha} < \infty$
是有限的. 这意味着, 如果边界密度趋于零, 而质量斜率保持非退化, 则对应的物理厚度是有限的. 因此, 负通量会导致质量在有限时间内完全流出边界附近的衣领, 从而导致真空形成.
终端迹的构造与局部化: 命题 6.7 构造并局部化了终端密度迹. 作者证明了密度在终端时间一致收敛到一个非平凡剖面, 其零集恰好为 $\{R\}$ , 即真空仅在边界上形成.

总结来说, 本文的核心工具包括:

横向应变率函数 $$U=u/r$$ 和 $$V=v/r$$ : 用于处理中心奇异性, 保留速度信息.
改进的最大原理: 用于建立速度在中心的 Lipschitz 连续性和倒数密度的估计.
Dini-Gronwall 估计: 用于封闭密度下界估计.
有限容量机制: 用于证明负边界迹下的有限时间真空形成.
质量坐标: 用于系统地描述解的行为, 特别是在边界附近.

4. 与之前工作的比较

本文与既有 PDE/分析理论文献的比较可以从以下几个维度展开:

参数范围的推进: 此前的全局存在结果在高维设置中仅限于 $\alpha \le 1$ 的情况, 其中端点 $\alpha=1$ 对应粘性 Saint-Venant (浅水) 系统. 具体来说, Zhang [28] 证明了 $$N=2$$ 时 $4/5 \le \alpha < 1$ 和 $$N=3$$ 时 $7/8 \le \alpha < 1$ 的全局径向强解. Huang-Meng-Zhang [16] 将此扩大到 $$N=2$$ 时 $1/2 < \alpha < 1$ 和 $$N=3$$ 时 $0.686 < \alpha < 1$ . Lei [17] 进一步缩小了下界, 证明了 $$N=2$$ 时 $0.54369 < \alpha < 1$ 和 $$N=3$$ 时 $0.67661 < \alpha < 1$ 的真空不形成. 本文首次将全局存在理论扩展到超临界区域 $\alpha > 1$ , 具体范围为 $1 < \alpha < 11.7$ 且 $2\alpha-1 \le \gamma < 3\alpha-1$ . 这是多维设置中 $\alpha > 1$ 的首个全局存在结果.
初始数据的大小: 本文的结果适用于任意大的初始数据, 只要初始有效速度在边界上非负. 与此对比, 一些先前工作可能需要小扰动假设或对初始数据施加额外限制. 例如, Matsumura-Nishida [21] 的小扰动理论不适用于大初值. 虽然一些弱解理论允许大初值, 但它们不保证经典解的正性. 本文的强解全局存在性是大初值的重要进展.
真空形成的位置与机制: 早期球对称真空结果通常将异常行为定位在球心 $$r=0$$ . 例如, Xin-Yuan [25] 证明真空不会在原点之外形成, Guo-Wang-Wang [12] 的理论中 positivity 在 $$r=0$$ 附近受到保护. 本文的显著不同之处在于, 构造的密度损失发生在外边界 $$r=R$$ , 且仅发生于此. 定理 2.5 明确证明了在任何情况下, 有限时间真空形成只能发生在边界, 而不能发生在内部或中心. 这与先前工作的焦点形成了鲜明对比.
符号判据的尖锐性: 本文建立了一个清晰的二分法: 全局正则性与奇异性形成, 完全由初始有效速度的符号决定. 如果 $v_0(R) \ge 0$ , 则解全局存在; 如果 $$v_0(R) < 0$$ , 则可以构造在有限时间内在边界上形成真空的解. 这种基于单一边界量的尖锐判据在先前工作中尚未出现. 例如, Haspot [14] 和 Constantin-Drivas-Nguyen-Pasqualotto [8] 的一维结果依赖于关于有效量的一边或 Oleinik 型假设, 而非简单的符号条件. 本文的符号判据更为简洁和直接.
技术方法的创新: 与先前工作相比, 本文的技术方法有几个关键创新:
- 横向应变率函数: 引入 $$U=u/r$$ 和 $$V=v/r$$ 来处理中心奇异性, 这在先前的高维球对称理论中未被充分利用.
- 速度在中心的 Lipschitz 连续性: 首次建立了这一结果, 解决了速度梯度的奇异性. 先前工作可能依赖于较弱的正则性或回避了这一问题.
- Dini-Gronwall 估计: 用于封闭密度下界估计, 结合改进的最大原理. 这种方法在处理超临界区域的粘性退化问题时提供了新的工具.
- 有限容量机制: 用于证明负边界迹下的有限时间真空形成. 这一机制利用了 $\alpha > 1$ 时积分的有限性, 是超临界区域特有的.
与弱解理论的区分: 本文强调, 全局弱解的存在性并非所解决的问题. 弱解可能包含真空, 并不能决定一个初始为正的经典分支是否保持与零均匀分离. 因此, 本文的强解全局存在性是对弱解理论的重要补充, 提供了关于经典解行为的更多信息.
与自由边界问题的区别: 本文处理的是固定边界上的无滑移条件, 而自由边界问题（如 Guo-Xu-Zhang [13] 和 Chen-Zhang-Zhu [7]）涉及移动界面. 本文的真空形成发生在固定边界上, 是一种“空化”现象, 而非自由边界的运动. 这使得问题在数学结构上有所不同.
局限性:
- 定理 2.1 是条件性的, 依赖于密度的上界. 只有在指数范围 (7) 内, 上界才可用, 从而得到无条件的全局存在性.
- 定理 2.7 是一个显式构造, 不断言所有具有 $$v_0(R) < 0$$ 的数据都会形成真空. 因此, 符号判据在负方向上是充分的, 但不是必要的.
- 参数范围 $1 < \alpha < 11.7$ 和 $2\alpha-1 \le \gamma < 3\alpha-1$ 仍有改进空间. 未来的工作可能会扩大这一范围.

总结来说, 本文在参数范围、初始数据大小、真空形成位置、符号判据尖锐性以及技术方法等方面均对既有文献做出了重要推进. 它不仅填补了 $\alpha > 1$ 区域的多维强解全局存在空白, 还揭示了一个简洁的物理图像: 初始边界有效速度的符号决定了流体的长期行为.

10. $\mathbb{R}^3$ 上 Euler 方程奇异性稳定性框架

Stability Framework for the Singularity of the Euler Equations on $\mathbb{R}^3$

作者: Valentin Duruisseaux, Adarsh Ganeshram, Robert J. George, Anima Anandkumar

类别: math.AP, math.NA

发布日期: 2026-09-09

arXiv ID: 2609.10860v1

摘要翻译: 在一项伴随的数值研究中，我们在无界区域 $\mathbb{R}^3$ 上找到了 Euler 方程的高精度奇异剖面。本手稿通过详细建立证明近似自相似剖面（非线性）稳定性的框架来补充该研究，将分析化简为大量但有限的显式估计和可计算常数。在严格认证论证中出现的估计与常数的条件下，该框架完成了稳定性证明，并且关键在于允许将稳定的重标度剖面重构为原始变量下的一个容许解，该解在有限物理时间内变得奇异。随着稳定性和重构机制的建立，剩余工作主要是定量的：认证显式常数，并在必要时锐化选定的解析估计以完成稳定性证明。

Comment: 112 pages

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

1. 研究对象和背景

在偏微分方程分析中, 一个核心问题是: 从光滑初值出发的解是否能在有限时间内失去正则性, 即发生奇异性形成 (singularity formation) 或爆破 (blowup). 特别是, 三维不可压缩 Navier-Stokes 方程是否存在有限时间奇点是流体力学中最著名的未解决问题之一, 也是千禧年问题之一. 与 Navier-Stokes 方程密切相关的是三维不可压缩 Euler 方程的奇异性问题.

本文研究的对象是定义在无界域 $\mathbb{R}^3$ 上的三维轴对称带旋涡 (axisymmetric with swirl) 不可压缩 Euler 方程:

\partial_t \tilde{u} + (\tilde{u} \cdot \nabla) \tilde{u} = -\nabla \tilde{p}, \quad \nabla \cdot \tilde{u} = 0,

其中 $\tilde{u}(x,t) \in \mathbb{R}^3$ 是速度场, $\tilde{p}(x,t) \in \mathbb{R}$ 是压力标量. 通过取旋度, 可以得到涡量形式:

\partial_t \tilde{\omega} + \tilde{u} \cdot \nabla \tilde{\omega} = \tilde{\omega} \cdot \nabla \tilde{u}, \quad \tilde{\omega} := \nabla \times \tilde{u}.

在柱坐标 $(\hat{r}, \theta, \hat{z})$ 下, 利用轴对称性 ( $\partial_\theta = 0$ ), 并引入重标度的旋涡和流函数变量:

u^\bullet := \frac{u_\theta}{\hat{r}}, \quad \omega^\bullet := \frac{\omega_\theta}{\hat{r}}, \quad \psi^\bullet := \frac{\psi_\theta}{\hat{r}},

可以消除柱坐标中的 $1/\hat{r}$ 奇异性. 于是得到关于 $(u^\bullet, \omega^\bullet, \psi^\bullet)$ 的简化轴对称系统:

\partial_t u^\bullet + u_r \partial_{\hat{r}} u^\bullet + u_z \partial_{\hat{z}} u^\bullet = 2u^\bullet \partial_{\hat{z}} \psi^\bullet,

\partial_t \omega^\bullet + u_r \partial_{\hat{r}} \omega^\bullet + u_z \partial_{\hat{z}} \omega^\bullet = 2u^\bullet \partial_{\hat{z}} u^\bullet,

-\left[ \partial_{\hat{r}}^2 + \frac{3}{\hat{r}} \partial_{\hat{r}} + \partial_{\hat{z}}^2 \right] \psi^\bullet = \omega^\bullet.

为了研究可能的有限时间奇点, 本文采用行波自相似 ansatz (traveling self-similar ansatz):

u^\bullet(\hat{r}, \hat{z}, t) = (T-t)^{c_u} u\left( \frac{\hat{r}}{(T-t)^\lambda}, \frac{\hat{z}-\hat{z}_c(t)}{(T-t)^\lambda} \right),

\omega^\bullet(\hat{r}, \hat{z}, t) = (T-t)^{c_\omega} \omega\left( \frac{\hat{r}}{(T-t)^\lambda}, \frac{\hat{z}-\hat{z}_c(t)}{(T-t)^\lambda} \right),

\psi^\bullet(\hat{r}, \hat{z}, t) = (T-t)^{c_\psi} \psi\left( \frac{\hat{r}}{(T-t)^\lambda}, \frac{\hat{z}-\hat{z}_c(t)}{(T-t)^\lambda} \right),

其中 $$T$$ 是爆破时间, $\lambda$ 是空间爆破率, $\hat{z}_c(t)$ 是行波中心, $c_u, c_\omega, c_\psi$ 是振幅指数. 通过平衡各项的量级, 确定这些指数为:

c_u = -1, \quad c_\omega = -1-\lambda, \quad c_\psi = -1+\lambda.

将上述 ansatz 代入原方程, 并设定 $\lambda = 1/2$ , 可以得到关于重标度变量 $(U, \Omega, \Psi)$ 的稳态剖面方程:

U + (\lambda r + U_r) \partial_r U + (C + \lambda z + U_z) \partial_z U = 2U \partial_z \Psi,

(1+\lambda) \Omega + (\lambda r + U_r) \partial_r \Omega + (C + \lambda z + U_z) \partial_z \Omega = 2U \partial_z U,

-E \Psi = \Omega,

其中 $E := \partial_r^2 + \frac{3}{r}\partial_r + \partial_z^2$ 是椭圆算子, $$C$$ 是行波速度参数.

已有研究表明, 在有边界的情况下 (Chen and Hou, 2025, 2026) 或在非光滑初值条件下 (Elgindi, 2021; Shkoller, 2026; Chen, 2026a,b), 可以证明 Euler 方程存在有限时间爆破. 然而, 对于无界域 $\mathbb{R}^3$ 上的光滑初值, 发现近似奇异剖面的难度极大. 一种常见的策略是通过数值方法构造近似解, 然后利用稳定性理论证明存在一个与该近似解足够接近的奇异解. 自相似 ansatz 为此提供了一条可行的路径.

在伴随的数值研究中 (Ganeshram et al., 2026), 作者使用物理信息神经网络 (PINNs) 获得了高精度的自相似剖面, 并将其转换为分段多项式样条表示, 以便进行严格的认证. 样条表示允许精确的区间内微分和评估, 并使用任意精度区间算术库 Arb 严格控制数值误差. 这使得剖面残差可以被严格界定为 $$L^2$$ 范数下的 $O(10^{-5})$ 和 $L^\infty$ 范数下的 $O(10^{-3})$ .

然而, 建立稳定性是至关重要且复杂的. 当满足全局流出性质 (global outgoing property) 时, 稳定性分析会简化, 因为这种膨胀效应会将自相似流动推离奇异核心, 并提供有利的高阶阻尼结构. 但对于带有自相似 ansatz 的 Euler 系统, Constantin et al. (2026) 提供了强有力的理论证据表明, 子午线流动的非平凡不动点很可能存在, 并且"众所周知会造成巨大的麻烦". 这促使作者用更弱的局部条件替代严格的全局流出条件, 该局部条件仅要求自相似流动远离这些不动点, 防止扰动在这些点集中, 从而为稳定性分析提供有利的结构. 然而, 在这种较弱流出性质下建立稳定性需要比先前工作精细得多的分析.

2. 主要定理或主要结果

本文建立了近似爆破剖面非线性稳定性的证明框架, 将该论证归约为大量但有限的显式、可严格验证的估计集合. 具体而言, 主要结果包括:

动态重标度方程的建立: 基于行波自相似 ansatz, 推导了适用于 Euler 方程的动态重标度方程. 这些方程作为围绕给定近似剖面进行稳定性分析的起点.
线性化分析: 在样条近似剖面周围对适当的动态重标度方程进行线性化, 并在合适的加权 $$L^2$$ 和 Sobolev 范数下表述稳定性定理. 引入了调制方程 (modulation equations) 来处理尺度一致性变量, 并定义了闭合的线性化变量.
加权范数和能量估计: 定义了低阶和高阶加权能量泛函, 并导出了互补的低阶和高阶能量估计, 分别针对线性项和非线性项. 这些估计具有显式的可计算界.
稳定性定理: 证明了单半径稳定性定理 (Single-Radius Stability Theorems) 和双半径稳定性定理 (Two-Radii Stability Theorem). 这些定理在临界标度 $\lambda = 1/2$ 下成立, 控制了调制参数在其自相似值附近的稳定性.
计算机辅助证明策略: 提出了一种计算机辅助策略来证明稳定性, 将证明归约为对大量显式估计的严格认证. 这些估计可以通过区间算术或 Arb 计算进行严格认证.
重构机制: 一旦稳定性和重构机制建立, 剩余的工作主要是定量的: 认证显式常数, 并在必要时锐化选定的解析估计以关闭稳定性证明. 在临界标度 $\lambda = 1/2$ 下, 控制调制参数在其自相似值附近将导致候选剖面的稳定性, 并在重构后得到一个在原变量中有限时间变得奇异的容许解.
LeanPDE 形式化: 并行开发了一个 Lean 形式化, 重点在于符号推导和分析步骤的证明. 使用 Mathlib 和 TorchLean, 该形式化将 Lean 中的符号验证与在 Julia 中使用 Arb 区间算术进行的严格数值计算连接起来, 形成了机器检查 PDE 分析的综合框架的一部分.

关键假设包括: 使用 $\lambda = 1/2$ 的临界标度, 以及较弱的局部流出条件而非严格的全局流出条件. 代表性公式包括动态重标度方程、线性化方程、加权能量泛函的定义, 以及稳定性定理中的显式常数界.

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

本文的核心技术障碍在于如何在没有严格全局流出条件的情况下, 证明近似自相似剖面的非线性稳定性. 作者的真正创新点在于发展了一套完整的分析框架, 将稳定性证明归约为大量显式、可计算的估计, 并利用计算机辅助方法进行严格认证.

1. 动态重标度和线性化框架:

作者首先建立了适用于 Euler 方程的动态重标度方程. 这是稳定性分析的基础, 因为它将原始的时间依赖问题转化为围绕固定剖面的扰动问题.
在线性化过程中, 作者引入了尺度一致性变量 (scale-consistent variables) 和调制方程 (modulation equations). 这些工具处理了由于自相似标度引起的变量之间的耦合, 确保了线性化系统的闭合性.
关键的技术细节在于如何选择适当的变量组合, 使得线性化后的系统具有良好的结构, 便于后续的能量估计.

2. 加权范数和能量方法:

作者定义了低阶和高阶加权能量泛函. 低阶能量捕捉主要的动力学行为, 而高阶能量提供额外的正则性和阻尼.
互补的低阶和高阶能量估计是关键创新. 单独的低阶能量不足以关闭稳定性论证, 需要高阶能量的贡献来补偿某些不利项.
权重函数的选择受到严格约束, 以确保能量估计的有效性. 这些约束在附录 B.5 中详细讨论.

3. 流函数椭圆估计:

由于 Euler 方程中包含椭圆型的流函数关系, 对流函数的估计至关重要. 作者发展了专门的流函数椭圆估计 (streamfunction elliptic estimates), 包括认证设置、奇异积分约定、约化核和二阶导数算子的完整处理.
附录 B.6 详细展示了这些估计的推导, 包括点常数、加权算子常数和块界的计算. 这些估计将源项的导数转化为能量范数, 是关闭非线性估计的关键.

4. 计算机辅助证明策略:

作者提出了一个四层的认证管道: 从 PINN 数值发现到样条基表示, 再到区间证书, 最后是计算机辅助稳定性证明.
样条表示的选择是关键的, 因为它允许精确的区间内微分和评估, 可以使用标准严格工具 (如区间算术) 进行评估和界定.
附录 C-F 详细展示了低阶线性估计、高阶线性估计、低阶非线性估计和高阶非线性估计的具体实现, 包括矩阵构造、尾部闭合、局部矩阵证书和直接项的处理.
附录 G 提供了稳定性定理的完整证明, 展示了如何将各个估计组合起来以关闭整个论证.

5. Lean 形式化:

作者并行开发了 Lean 形式化, 使用 Mathlib 和 TorchLean 进行符号验证. 这与 Julia 中的 Arb 区间算术计算相结合, 形成了一个综合的机器检查框架.
这种形式化不仅提高了证明的可信度, 还为未来的 PDE 分析提供了一个可复用的基础设施.

6. 局部流出条件:

作者用较弱的局部流出条件替代了严格的全局流出条件. 这一改变使得稳定性分析更加复杂, 但更具一般性. 局部条件仅要求自相似流动远离子午线流动的非平凡不动点, 防止扰动在这些点集中.
这一创新使得框架能够应用于更广泛的场景, 尽管需要更精细的分析来处理由此带来的额外困难.

4. 与之前工作的比较

本文与既有工作的比较可以从多个维度展开:

1. 与有边界 Euler 方程爆破结果的比较:

Chen and Hou (2025, 2026) 在有边界情况下证明了 Euler 方程的有限时间爆破. 这些结果依赖于固体壁面施加的约束来产生快速的涡量增长 (Luo and Hou, 2014a,b).
本文的研究对象是无界域 $\mathbb{R}^3$ 上的 Euler 方程, 这是一个更具挑战性的问题, 因为没有边界提供的约束机制. 因此, 发现近似奇异剖面的难度极大, 需要新的方法和工具.

2. 与非光滑初值爆破结果的比较:

Elgindi (2021), Shkoller (2026), Chen (2026a,b) 在非光滑初值条件下证明了 Euler 方程的有限时间爆破.
本文关注的是光滑初值的情况, 这是更接近实际物理场景的问题. 非光滑初值的爆破机制可能不直接适用于光滑初值, 因此需要独立的发展.

3. 与 Navier-Stokes 方程的比较:

Navier-Stokes 方程的有限时间爆破是一个更为困难的问题, 因为粘性会耗散与小尺度结构相关的任何潜在奇异性 (Tao, 2007, 2015).
本文专注于 Euler 方程, 避免了粘性带来的复杂性, 但仍面临非线性项的挑战. Euler 方程的奇异性研究可以为理解 Navier-Stokes 方程提供重要的见解.

4. 与先前稳定性工作的比较:

Chen et al. (2021) 在先前的工作中建立了满足全局流出性质的自相似流的稳定性. 这种全局流出性质提供了有利的膨胀效应, 将自相似流动推离奇异核心, 并产生有利的高阶阻尼.
本文的创新在于用较弱的局部流出条件替代了严格的全局流出条件. 这一改变使得框架更具一般性, 但需要更精细的分析. Constantin et al. (2026) 的理论证据表明, 子午线流动的非平凡不动点很可能存在, 这使得全局流出条件难以满足. 因此, 局部条件是一个必要的妥协.

5. 与数值方法的比较:

在伴随的数值研究 (Ganeshram et al., 2026) 中, 作者使用 PINNs 获得高精度剖面, 并将其转换为样条表示.
与传统的谱方法或全局多项式表示相比, 样条表示在局部准确性与严格认证之间取得了更好的平衡. 样条的局部性质允许仅在剖面快速变化的区域细化网格, 而不必不必要地增加其他地方近似的复杂性.
这种表示方式使得残差、范数和导数的严格认证成为可能, 这对于稳定性论证至关重要.

6. 与计算机辅助证明的比较:

本文提出的计算机辅助证明策略结合了符号验证 (Lean) 和严格数值计算 (Arb 区间算术), 形成了一个综合的机器检查框架.
这种方法论的创新在于它不仅提高了证明的可信度, 还为未来的 PDE 分析提供了一个可复用的基础设施. 通过将证明归约为大量显式、可计算的估计, 作者使得整个论证过程更加透明和可验证.

7. 局限性和代价:

本文的框架仍然依赖于对大量显式常数的严格认证, 这是一个计算密集型的任务. 虽然模块化结构允许在不改变基本稳定性机制的情况下进行调整, 但完成所有认证可能需要大量的计算资源.
局部流出条件的引入使得分析更加复杂, 可能需要更多的技巧来处理由此带来的额外困难.
目前的结果是条件性的: 只有在严格认证所有估计和常数之后, 稳定性证明才能完全关闭. 因此, 本文提供了一个强大的框架, 但最终结论的实现仍有待完成.

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

1. 随机风驱动边界条件下原始方程的全局适定性

Global well-posedness of the primitive equations with stochastic wind-driven boundary conditions

作者: Tarek Zoechling

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11734v1

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

摘要翻译

中文标题: 随机风驱动边界条件下原始方程的全局适定性

摘要翻译: 考虑受有限个布朗运动驱动的 Neumann 风应力作用的三维原始方程。对于任意光滑空间风廓线，在不假设初始数据具有垂直导数的前提下，证明了全局逐路径存在性与唯一性，从而将 Binz, Hieber, Hussein, and Saal~\cite{BHHS-24} 针对此类数据与风廓线的局部理论进行了推广。该证明结合了局部最大 $\rL^2$ 正则性构造与边界随机卷积的加权估计。

2. 奇异非局部输运的定量扩散极限

Quantitative Diffusive Limits for Singular Nonlocal Transport

作者: Andrea Agazzi, Giuseppe Bruno, Federico Pasqualotto, Philippe Rigollet

类别: math.AP, stat.ML

发布日期: 2026-09-10

arXiv ID: 2609.11837v1

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

摘要翻译

中文标题: 奇异非局部输运的定量扩散极限

摘要翻译: 我们研究闭连通黎曼流形上的非局部连续性方程

\partial_t μ_b = \operatorname{div}\!\left( μ_b \nabla \log \bigl((I-b^2Δ)^{-1}μ_b\bigr) \right)

对于光滑严格正的初始数据，我们证明当 $b \to 0$ 时，其全局解以最优且关于时间一致的速率收敛到热流 $$μ(t)$$ ：

\sup_{t\ge0}\|μ_b(t)-μ(t)\|_{L^1}\le Cb^2.

关键估计是 $$b$$ -加权的高阶解析能量的一致耗散，该估计在缺乏 Wasserstein 梯度流结构的情况下仍能得到指数弛豫。在圆上，我们还分析了相应的确定性 $$N$$ -粒子动力学。弱–强调制能量论证给出

中文标题: 全次临界长程情形下非线性 Schrödinger 方程的修正波算子

摘要翻译: 我们针对非线性 Schrödinger 方程 $i\partial_tu+\frac12Δu=|u|^pu$ 在全次临界长程情形 $$0<p<2/d$$ 下构造了修正波算子，其中末态数据 $$W(x)$$ 为小、非零且解析，其对数梯度有界。此前的结果仅对某些特定类别的 Cauchy 数据建立了大时间渐近行为；此外，当 $$p<1/d$$ 时，精确的渐近展开式仍未被确定。当 $$1/d<p<2/d$$ 时，我们的结果给出了如下近似：

\frac{1}{(it)^{\frac{d}{2}}} e^{\frac{i|x|^2}{2t}} W\left(\frac{x}{t}\right) \exp\left[ -i\frac{t^{1-\frac{dp}{2}}-1}{1-\frac{dp}{2}} \left|W\left(\frac{x}{t}\right)\right|^p \right]

该波算子通过半径递减的解析空间中的迭代过程构造而成。当 $p\le 1/d$ 时，我们通过 Fuchsian 方程与输运方程耦合的有限截断来构造剖面（profile）。该剖面仍遗留一个长程三角耦合，其终端积分不保持所需的快速衰减类。该构造给出了 $2\le q\le\infty$ 时的定量 $$L^q$$ 渐近行为，并在预设的解析渐近类中证明了唯一性。核心的新要素包括一种消除该长程耦合的非线性末态正规形（normal form），以及在特定解析空间中的混合迭代方法。

13. 半空间上的分数阶 Hardy–Maz’ya 不等式

Fractional Hardy–Maz’ya inequality on a half-space

作者: Michał Kijaczko, Antoni Szczukiewicz

类别: math.AP, math.FA

发布日期: 2026-09-09

arXiv ID: 2609.10832v1

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

摘要翻译

中文标题: 半空间上的分数阶 Hardy–Maz’ya 不等式

摘要翻译: 本文的主要目的是在半空间上给出著名的 Maz’ya 不等式的分数阶对应形式，即

\int_{\mathbb{R}^{d}_{+}}\int_{\mathbb{R}^{d}_{+}}\frac{|u(x)-u(y)|^p}{|x-y|^{d+sp}}dy\,dx\ge\mathcal{D}_{d,s,p}\int_{\mathbb{R}^{d}_{+}}\frac{|u(x)|^p}{x_d^{sp}}dx+C_{d,s,p,τ}\int_{\mathbb{R}^{d}_{+}}\frac{|u(x)|^p}{x_{d}^{sp-τ}\left(x_{d-1}^2+x_d^2\right)^{τ/2}}dx,

其中 $\mathcal{D}_{d,s,p}$ 表示半空间 $\mathbb{R}^{d}_{+}$ 上分数阶 Hardy 不等式的最佳常数。我们还在 Sobolev–Bregman 形式的设定下得到了类似的结果。

14. 带非线性阻尼的半线性波动方程的最优寿命估计

Sharp Lifespan Estimates for a Semilinear Wave Equation with Nonlinear Damping

作者: Firas Kaabi

类别: math.AP

发布日期: 2026-09-09

arXiv ID: 2609.10819v1

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

摘要翻译

arXiv ID: 2609.11494v1

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

摘要: In this paper, we study the coupled singular weighted elliptic system

\left\{ \begin{aligned} -\operatorname{div}(|x|^{-2a}\nabla u)+λ_1 \frac{u}{|x|^{2a}}&=βp\frac{|v|^q|u|^{p-2}u}{|x|^{b(p+q)}}+\frac{|u|^{2^{\sharp}-2}u}{|x|^{b2^{\sharp}}},\quad\text{in}~\mathbb{R}^N, -\operatorname{div}(|x|^{-2a}\nabla v)+λ_2\frac{v}{|x|^{2a}}&=βq\frac{|u|^p|v|^{q-2}v}{|x|^{b(p+q)}}+\frac{|v|^{2^\sharp-2}v}{|x|^{b{2^\sharp}}},\quad\text{in}~\mathbb{R}^N, \displaystyle\int_{\mathbb{R}^N}\frac{|u|^2}{|x|^{2a}} dx=ρ_1^2,\quad &\int_{\mathbb{R}^N}\frac{|v|^2}{|x|^{2a}}dx=ρ_2^2. \end{aligned} \right.

where $N\ge3$ , $β\in\mathbb{R}$ , $$ρ_1,ρ_2>0$$ , $\max\{0,\tfrac{N-4}{2}\}\le a<\tfrac{N-2}{2}$ , $$a<b<a+1$$ , $2^{\sharp}:=\frac{2N}{N-2(1+a-b)}$ , $$p,q>1$$ and $2<p+q<2^{\sharp}$ . Here $2^\sharp$ is the critical exponent of the Caffarelli–Kohn–Nirenberg inequality, and $λ_1,λ_2\in\mathbb{R}$ are unknown, since the two weighted masses are prescribed. We show that the ground state level is not attained for $β\le0$ . For $$β>0$$ we obtain positive normalized ground states, with positive Lagrange multipliers, in the mass subcritical, mass critical, and mass supercritical regimes, for explicit ranges of $$β$$ . In the mass supercritical regime, a threshold $β_0\ge0$ appears: a ground state exists for every $$β>β_0$$ and does not exist for $$0<β<β_0$$ . We give conditions on $$p,q$$ under which $$β_0=0$$ , and we prove that $$β_0>0$$ when $p,q\ge2$ . Finally, we describe the ground states as $β\to0^+$ and as $β\to+\infty$ : after a dilation, they converge to the ground states of the limit system without critical terms, or one component vanishes and the other concentrates on an extremal of the Caffarelli–Kohn–Nirenberg inequality.

10. 调和函数节点体积的 sharp 下界

A sharp lower bound for the nodal volume of harmonic functions

作者: Zhehui Wang

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11444v1

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

摘要: Let $$u$$ be a non-zero real-valued harmonic function in $B_4\subset\mathbb{R}^n$ with $n\geq3$ and $$u(0)=0$$ . We prove that

\mathcal{H}^{n-1}\bigl(\{u(x)=0\}\cap B_2\bigr)\ge C_n\mathcal{N},

where $$C_n$$ is a positive constant depending only on $$n$$ , and $\mathcal{N}$ is the doubling index defined by

\mathcal{N}=\log_2\frac{\sup_{B_1}|u|}{\sup_{B_{\frac{1}{2}}}|u|}.

The linear dependence on $\mathcal{N}$ is optimal. This estimate confirms a folklore conjecture on the nodal volume of harmonic functions.

11. 最优几何不等式与全非线性共形流

Optimal geometric inequalities and fully nonlinear conformal flows

作者: Yuxin Ge, Guofang Wang, Wei Wei

类别: math.DG, math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11421v1

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

摘要: We establish sharp Sobolev-type geometric inequalities on $\mathbb{S}^n$ involving the total $$σ_k$$ -curvatures $\int_{\mathbb{S}^n}σ_k(g)\,dv_g$ . These results extend the optimal inequalities of Guan–Wang~\cite{GWDuke} from the cone $\mathcal{C}_k$ to the strictly larger cone $\mathcal{C}_{k-1}$ , thereby enlarging the range of admissible conformal metrics. Our approach is variational and is implemented through a fully nonlinear conformal flow. Working in $\mathcal{C}_{k-1}$ introduces substantial analytic difficulties; in particular, one must obtain $$C^2$$ a priori estimates while simultaneously verifying that the flow remains parabolic. We resolve these issues via a carefully designed test function and by applying the maximum principle to the maximal eigenvalue of the Hessian matrix. As applications, we solve two open problems in dimensions 3 and 4. Finally, we give examples to show that these inequalities cannot be extended to $\mathcal{C}_{k-2}$ .

12. 临界指数下非局部能量均匀化中凸奇异项的涌现

Emergence of a convex strange term via homogenization of non-local energies at the critical exponent

作者: Giuseppe Cosma Brusca, Giuliana Fusco

类别: math.AP, math.OC

发布日期: 2026-09-10

arXiv ID: 2609.11371v1

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

摘要: We derive the $$Γ$$ -limit of convolution-type and discrete energies subject to Dirichlet boundary conditions on periodically perforated domains at the critical exponent. We assume that the length-scale of the non-local interactions is much smaller than the side-length of the cubic perforations and prove that a separation of scales occurs. Exploiting the analogies between the variational frameworks of our interest, we employ a unified argument to overcome the technical difficulties that arise from the scaling invariance of the energies. Our multiscale analysis yields a novel observation that is not related to the non-local nature of the functionals, but rather to the analysis at the critical exponent: we prove that the energy density of the $\textit{strange term}$ is convex, even in the vector-valued setting.

13. Helmholtz 方程的笛卡尔 PML 截断在高频率下具有指数精度

Cartesian PML truncation for the Helmholtz equation is exponentially accurate at high frequency

作者: Martin Averseng, Jeffrey Galkowski, Euan A. Spence

类别: math.AP, math.NA

发布日期: 2026-09-10

arXiv ID: 2609.11343v1

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

摘要: For a wide class of Helmholtz scattering problems (namely, those fitting in the black-box framework of Sjöstrand–Zworski) we prove that the error incurred by approximating the radiation condition by a Cartesian perfectly-matched layer (PML) decreases exponentially in the width of the PML, the strength of the PML scaling, and the frequency.

14. 由梯度值确定的自由边界的正则性. 第一部分：三阶渐近展开

Regularity of a free boundary determined by the value of the gradient. Part 1: Third order asymptotics

作者: Frida Fejne

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11324v1

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

摘要: This is the first part of a project concerning the regularity of the free boundary of a function, where the free boundary depends on the gradient of the function. We study the minimizer of the expression \begin{equation*} J(u) : = \int_{B_1} F(\nabla u) \ dx, \end{equation*} where $$F(x)$$ is a uniformly convex function whose second derivatives might jump at $$|x|=1$$ . This results in an Euler-Lagrange equation that varies over the free boundary, which we define as

Γ= \partial \{x \in B_1: |\nabla u| >1 \} \cap \partial \{x \in B_1: |\nabla u| <1 \}.

We consider two-phase flat points for which the second derivatives of $$F$$ jump over $$Γ$$ . In this paper, we show that under some regularity and non-degeneracy assumptions, a minimizer $$u$$ can be expressed as

u(x)= a + ν\cdot x+ δp(x) + δεq(x),

where $a \in \mathbb{R}$ , $ν\in \mathbb{R}^n$ , $$0 < δ, ε<< 1$$ . Here $$p$$ is a $$C^1$$ function, which consists of one polynomial in the upper half ball and another polynomial in the lower half ball. The function $$q$$ is a rest term with bounded $$L^2$$ norm. Furthermore, assuming that we have a sequence, $$u^j$$ , of minimizers, we show that $$q^j$$ converges weakly to a function $$q^0$$ that satisfies a certain PDE. In addition, we show that $$q^0$$ is $$C^3$$ in $$B_r^+$$ and $$B_r^-$$ , respectively, for $$0<r<1$$ . This is the main result of this paper which is intended to be used to show regularity of the free boundary.

15. 关于涉及非局部无穷拉普拉斯算子的演化方程

On an evolutionary equation involving the nonlocal infinity Laplacian

作者: Frida Fejne

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11311v1

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

摘要: We study the evolutionary equation

\frac{\partial u}{\partial t} = (\mathcal{L}_{\infty})^αu \quad \text{in} \quad Ω,

where $(\mathcal{L}_{\infty})^αu$ denotes a nonlocal infinity Laplacian acting on the function $$u$$ , $0<α\leq 1$ and $$Ω$$ is a bounded open set in $\mathbb{R}^{n+1}$ . We prove existence and uniqueness using Perron’s method for the Dirichlet problem when $$Ω$$ is a cylinder and $$0<α<1$$ .

16. 带有低阶 Hardy-Littlewood-Sobolev 项和 $$L^2$$ -临界项的分数阶 Choquard 方程的锐质量阈值

A Sharp Mass Threshold for Fractional Choquard Equations with Lower Hardy-Littlewood-Sobolev and $$L^2$$ -Critical Terms

作者: Yongpeng Chen, Zhipeng Yang

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11284v1

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

摘要: For $\frac12\le s<1$ , we study a fractional Choquard equation with prescribed mass and nonlinearities at the lower Hardy-Littlewood-Sobolev and $$L^2$$ -critical exponents. The sharp Hardy-Littlewood-Sobolev and Choquard Gagliardo-Nirenberg inequalities determine an explicit critical mass $$a_*$$ . For $0<a\le a_*$ , we compute the exact infimum of the constrained energy and prove that it is not attained and that no normalized solution exists. For $$a>a_*$$ , the energy is unbounded from below on the mass sphere, while the Pohozaev set is nonempty.

17. 自由边界问题中的黏性与变分方法：献给 Sandro Salsa 的纪念文集

Viscosity and variational approaches in free boundary problems: A Volume in Honor of Sandro Salsa

作者: Daniela De Silva, Fausto Ferrari

类别: math.AP

发布日期: 2026-09-10

arXiv ID: 2609.11205v1

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

摘要: Free boundary problems (FBPs) are differential equations where the domain depends on the solution, with applications in diverse fields such as flame propagation, financial mathematics, and tumor growth. Classical examples include the Bernoulli problem, the obstacle problem, and the Stefan problem. The field has advanced significantly through variational and geometric approaches, notably by Alt, Caffarelli, and collaborators, who developed monotonicity formulas and viscosity solutions to address free boundary regularity. Extensions have included nonlinear operators, inhomogeneous problems, and lower-dimensional cases like the thin Bernoulli problem. Evolution problems, such as the Hele-Shaw flow, also connect to FBPs. To celebrate these developments and honor Sandro Salsa’s pioneering contributions, a special session on FBPs was held at the AMS-UMI joint meeting in Palermo (2024). This is the editorial of a volume dedicated to Sandro Salsa on the occasion of his 75th birthday.

18. 具有快速增长非线性的 $$p$$ -Laplacian Schrödinger–Maxwell 系统

A $$p$$ -Laplacian Schrödinger–Maxwell system with rapidly growing nonlinearities

作者: Genival da Silva

类别: math.AP

发布日期: 2026-09-09

arXiv ID: 2609.10878v1

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

摘要: Let $Ω\subset\R^N$ , $N\ge2$ , be a bounded domain and let $1<p<\infty$ . Inspired by the exponential case treated in \cite{BO2026}, we study the quasilinear Schrödinger–Maxwell system

\begin{cases} -Δ_p u+ψG'(u)=f &\text{in }Ω,\\ -Δ_pψ=G(u) &\text{in }Ω,\\ u=ψ=0 &\text{on }\partialΩ, \end{cases}

where $G\in C^1(\R)$ is even, convex, nonnegative and nontrivial, with $$G(0)=0$$ , and satisfies

G(t)\lesssim 1+|G'(t)| \qquad\text{for large }|t|.

Under the assumption

|f|\ln(1+|f|)\in L^1(Ω),

we prove the existence of a finite-energy weak solution and show that the solution is a saddle point of the associated functional.

19. 具有内部点不连续有界解的奇异漂移热方程

Heat equation with singular drift which has bounded solution discontinuous at an interior point

作者: N. V. Krylov

类别: math.AP, math.PR

发布日期: 2026-09-09

arXiv ID: 2609.10828v1

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

摘要: The goal of this paper is to present a one dimensional heat equation with a singular drift for which solutions of the Cauchy problem with smooth initial data are bounded and are discontinuous at an interior point (common to all solutions).

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

1. 可压缩等熵 Navier-Stokes 方程的点态衰减极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 二维不可压欧拉方程柯西问题已报道多解的数值研究#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 时空膨胀下 Favre 滤波可压缩 Navier–Stokes 系统的微分同变不变性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 具有水平粘性的二维 Navier–Stokes 方程的各向异性 Onsager 判据#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 一种双极 Euler-Poisson 系统的全局存在性与时间衰减：含一种无压且无阻尼流体#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. Sobolev 空间中具有退化临界点的二维非定常 Prandtl 系统适定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 二维趋化性Navier-Stokes方程的对数敏感性边界层. II. 粘性消失极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. $ rac{1}{2}$ 阶 Euler 方程在 $ rac{3}{4}$ 空间上的稳定奇点#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

9. 三维球对称可压缩 Navier-Stokes 方程中任意大初值强解全局存在与爆破的尖锐符号判据#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

10. $\mathbb{R}^3$ 上 Euler 方程奇异性稳定性框架#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 随机风驱动边界条件下原始方程的全局适定性#

摘要翻译#

2. 奇异非局部输运的定量扩散极限#

摘要翻译#

3. 重新审视 Schrödinger-Klein-Gordon 系统#

摘要翻译#

4. 高斯函数并不总是最大化Schrödinger方程的混合范数Strichartz不等式#

摘要翻译#

5. 关于最大 Riesz 变换的无维 $L^p$ 估计的更短证明#

摘要翻译#

6. 随机多孔介质方程中支撑传播的抛物-双曲分解#

摘要翻译#

7. 多物种 Wasserstein 梯度流的控制策略#

摘要翻译#

8. 非自治系数部分耗散双曲系统的长时间动力学#

摘要翻译#

9. 带反应项的连续性方程的拉格朗日叠加原理#

摘要翻译#

10. 某些常系数一维线性双曲平衡律的最小零控制时间及相关核方程的性质#

摘要翻译#

11. Neumann-Neumann 波形松弛法求解含时滞双曲型偏微分方程#

摘要翻译#

12. 全次临界长程情形下非线性 Schrödinger 方程的修正波算子#

摘要翻译#

13. 半空间上的分数阶 Hardy–Maz’ya 不等式#

摘要翻译#

14. 带非线性阻尼的半线性波动方程的最优寿命估计#

摘要翻译#

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

1. 可压缩等熵 Navier-Stokes 方程的点态衰减极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 二维不可压欧拉方程柯西问题已报道多解的数值研究

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 时空膨胀下 Favre 滤波可压缩 Navier–Stokes 系统的微分同变不变性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 具有水平粘性的二维 Navier–Stokes 方程的各向异性 Onsager 判据

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 一种双极 Euler-Poisson 系统的全局存在性与时间衰减：含一种无压且无阻尼流体

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. Sobolev 空间中具有退化临界点的二维非定常 Prandtl 系统适定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 二维趋化性Navier-Stokes方程的对数敏感性边界层. II. 粘性消失极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. $$ rac{1}{2}$$ 阶 Euler 方程在 $$ rac{3}{4}$$ 空间上的稳定奇点

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

9. 三维球对称可压缩 Navier-Stokes 方程中任意大初值强解全局存在与爆破的尖锐符号判据

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

10. $\mathbb{R}^3$ 上 Euler 方程奇异性稳定性框架

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 随机风驱动边界条件下原始方程的全局适定性

摘要翻译

2. 奇异非局部输运的定量扩散极限

摘要翻译

3. 重新审视 Schrödinger-Klein-Gordon 系统

摘要翻译

4. 高斯函数并不总是最大化Schrödinger方程的混合范数Strichartz不等式

摘要翻译

5. 关于最大 Riesz 变换的无维 $$L^p$$ 估计的更短证明

摘要翻译

6. 随机多孔介质方程中支撑传播的抛物-双曲分解

摘要翻译

7. 多物种 Wasserstein 梯度流的控制策略

摘要翻译

8. 非自治系数部分耗散双曲系统的长时间动力学

摘要翻译

9. 带反应项的连续性方程的拉格朗日叠加原理

摘要翻译

10. 某些常系数一维线性双曲平衡律的最小零控制时间及相关核方程的性质

摘要翻译

11. Neumann-Neumann 波形松弛法求解含时滞双曲型偏微分方程

摘要翻译

12. 全次临界长程情形下非线性 Schrödinger 方程的修正波算子

摘要翻译

13. 半空间上的分数阶 Hardy–Maz’ya 不等式

摘要翻译

14. 带非线性阻尼的半线性波动方程的最优寿命估计

摘要翻译