Home » 论文总结

2026年02月17日论文分析

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

February 17, 2026 · 58 min · 12091 words · P0lyno3ial

Table of Contents

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

生成时间: 2026年02月17日 04:36:00

重点关注论文数量: 8

了解领域论文数量: 4

不相关论文数量: 34

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

1. 可压缩流体-颗粒流动的小德博拉数极限

The small Deborah number limit for the compressible fluid-particle flows

作者: Zhendong Fang, Kunlun Qi, Huanyao Wen

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14412v1

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

详细分析

中文标题: 可压缩流体-颗粒流动的小德博拉数极限

研究对象和背景

本文研究的是可压缩流体-颗粒两相耦合系统在小德博拉数（Deborah number）极限下的宏观渐近行为，即当颗粒弛豫时间相对于流体动力学时间尺度趋于零时，该微观-宏观耦合系统的水动力学极限。其物理背景涵盖喷雾燃烧（如柴油发动机）、气溶胶沉降、生物医学输运、矿物加工等典型多相流动现象。在此类模型中，一类相被视为连续介质（流体），另一类则由大量离散微粒构成（颗粒相），二者通过动量交换相互作用。

论文所研究的数学模型是Vlasov-Fokker-Planck方程与可压缩Navier-Stokes方程的强耦合系统（VFP-CNS），其无量纲化形式为：

\begin{cases} \partial_t f^\varepsilon + \dfrac{1}{\varepsilon} v \cdot \nabla_x f^\varepsilon = \dfrac{1}{\varepsilon^2} \operatorname{div}_v \left[ \nabla_v f^\varepsilon + (v - \varepsilon u^\varepsilon) f^\varepsilon \right], \\ \partial_t (\rho^\varepsilon u^\varepsilon) + \operatorname{div}_x (\rho^\varepsilon u^\varepsilon \otimes u^\varepsilon) + \nabla_x P(\rho^\varepsilon) + L u^\varepsilon = \dfrac{1}{\varepsilon} \displaystyle\int_{\mathbb{R}^3} (v - \varepsilon u^\varepsilon) f^\varepsilon \, dv, \\ \partial_t \rho^\varepsilon + \operatorname{div}_x (\rho^\varepsilon u^\varepsilon) = 0, \end{cases}

其中 $f^\varepsilon = f^\varepsilon(t,x,v) \ge 0$ 是颗粒在位置 $x \in \mathbb{R}^3$ 、速度 $v \in \mathbb{R}^3$ 处的概率密度分布函数； $\rho^\varepsilon = \rho^\varepsilon(t,x) > 0$ 和 $u^\varepsilon = u^\varepsilon(t,x)$ 分别为流体密度与速度场；压力函数取绝热形式 $P(\rho) = A \rho^\gamma$ ，其中 $$A > 0$$ , $\gamma > 1$ ；算子 $L u = -\mu \Delta_x u - (\mu + \lambda) \nabla_x \operatorname{div}_x u$ 是标准Lamé算子，对应牛顿流体本构关系，且满足 $\mu > 0$ , $\mu + \lambda > 0$ 。该系统定义在全空间 $\mathbb{R}^3$ 上，且要求解在无穷远处趋于平衡态：

(f^\varepsilon, u^\varepsilon, \rho^\varepsilon) \to (M(v), 0, 1), \quad \text{as } |x| \to +\infty,

其中 $M(v) = (2\pi)^{-3/2} e^{-|v|^2/2}$ 是标准归一化Maxwellian分布。初始数据被假定为围绕稳态的小扰动：

(f^\varepsilon, u^\varepsilon, \rho^\varepsilon)|_{t=0} = (f^\varepsilon_{\mathrm{in}}, u^\varepsilon_{\mathrm{in}}, \rho^\varepsilon_{\mathrm{in}}).

特别地，作者采用如下尺度参数设定：

\chi = 1,\quad \mathrm{Ma} = \varepsilon,\quad \mathrm{Re} = 1,\quad \mathrm{De} = \varepsilon^2,\quad P = \varepsilon^2,

这对应于轻颗粒（light particles）情形，即颗粒质量密度远小于流体密度： $\rho_P / \rho_F = O(\varepsilon^2)$ 。此尺度源自Goudon-Jabin-Vasseur的经典工作，并适配于“气泡区（bubbling regime）”建模。

该模型的关键物理机制在于：颗粒受布朗运动（Fokker-Planck扩散项）、外部流体拖曳力（ $(v - \varepsilon u^\varepsilon)f^\varepsilon$ ）及自身惯性共同作用；而流体则通过积分型动量源项 $\int (v - \varepsilon u^\varepsilon) f^\varepsilon \, dv$ 感知颗粒集体运动。当德博拉数 $\mathrm{De} \to 0$ ，意味着颗粒响应极快——其速度迅速弛豫至局部流体速度附近，从而在极限中退化为一种“被动示踪粒子”，其分布由流场决定。

因此，本文旨在严格建立如下奇异极限过程：当 $\varepsilon \to 0$ 时，上述VFP-CNS系统收敛至一个宏观双连续介质模型——即Navier-Stokes-Smoluchowski系统（NSS）：

\begin{cases} \partial_t n_0 + \operatorname{div}_x (u_0 n_0) = \Delta_x n_0, \\ \partial_t (\rho_0 u_0) + \operatorname{div}_x (\rho_0 u_0 \otimes u_0) + \nabla_x P(\rho_0) + L u_0 + \nabla_x n_0 = 0, \\ \partial_t \rho_0 + \operatorname{div}_x (\rho_0 u_0) = 0, \end{cases}

其中 $$n_0 = n_0(t,x)$$ 表示颗粒数密度（或质量浓度），它满足带对流-扩散结构的Smoluchowski方程；而流体部分则保留完整的可压缩Navier-Stokes结构，并新增一个由颗粒浓度梯度驱动的外力项 $\nabla_x n_0$ ，体现了颗粒-流体间动量反馈的宏观体现。

值得注意的是，该极限系统定义在相同全空间 $\mathbb{R}^3$ 上，边界条件为自然无穷远衰减条件 $(n_0,u_0,\rho_0)\to(1,0,1)$ 。函数空间设定为高阶Sobolev空间，例如 $$H^6_x$$ 对于宏观变量、 $H^4_{x,v}$ 对于微观余项，以支撑能量估计所需的正则性。

主要定理或主要结果

论文的核心结论由以下两个命题组成：

Proposition 1.1（NSS系统的全局适定性）：设初始数据 $(m^{\mathrm{in}}_0, u^{\mathrm{in}}_0, h^{\mathrm{in}}_0) := (n^{\mathrm{in}}_0 - 1, u^{\mathrm{in}}_0, \rho^{\mathrm{in}}_0 - 1)$ 满足
(i) $1 + \inf_{x\in\mathbb{R}^3} h^{\mathrm{in}}_0(x) > 0$ （保证密度正性）；
(ii) $(m^{\mathrm{in}}_0, u^{\mathrm{in}}_0, h^{\mathrm{in}}_0) \in H^6_x \times H^6_x \times H^6_x$ ，
且其范数足够小： $\| (m^{\mathrm{in}}_0, u^{\mathrm{in}}_0, h^{\mathrm{in}}_0) \|_{H^6_x}^2 \le \delta_0$ 。
则NSS系统存在唯一全局经典解 $$(m_0,u_0,h_0)$$ ，满足正则性：

(m_0,u_0) \in C^0([0,\infty); H^6_x) \cap C^1([0,\infty); H^4_x),\quad h_0 \in C^0([0,\infty); H^6_x) \cap C^1([0,\infty); H^5_x),

并满足一致能量耗散不等式：

\frac{1}{2}\frac{d}{dt} E_{\mathrm{ma}}(t) + C_0 D_{\mathrm{ma}}(t) \le 0,

其中 $E_{\mathrm{ma}}(t), D_{\mathrm{ma}}(t)$ 是定义于(3.5)-(3.6)的宏观能量与耗散泛函。

Theorem 1.1（VFP-CNS到NSS的小De极限）：令 $$(m_0,u_0,h_0)$$ 为Proposition 1.1所给NSS解。若VFP-CNS初始数据形如：

\begin{aligned} f^\varepsilon_{\mathrm{in}} &= (1+m^{\mathrm{in}}_0) M + \varepsilon \big[ (v\cdot u^{\mathrm{in}}_0)(1+m^{\mathrm{in}}_0) - v\cdot\nabla_x m^{\mathrm{in}}_0 \big] M + \varepsilon g^{\mathrm{in}} \sqrt{M}, \\ u^\varepsilon_{\mathrm{in}} &= u^{\mathrm{in}}_0 + \varepsilon u^{\mathrm{in}}, \\ \rho^\varepsilon_{\mathrm{in}} &= 1 + h^{\mathrm{in}}_0 + \varepsilon \rho^{\mathrm{in}}, \end{aligned}

且满足小初值条件：

\|g^{\mathrm{in}}\|_{H^4_{x,v}}^2 + \|(u^{\mathrm{in}},\rho^{\mathrm{in}})\|_{H^4_x}^2 + \|m^{\mathrm{in}}_0\|_{H^6_x}^2 + \|(u^{\mathrm{in}}_0,h^{\mathrm{in}}_0)\|_{H^5_x}^2 \le \delta,

则存在 $\varepsilon_0 > 0$ ，使得对任意 $\varepsilon \in (0,\varepsilon_0]$ ，VFP-CNS系统存在唯一全局解 $(f^\varepsilon,u^\varepsilon,\rho^\varepsilon)$ ，可表示为：

\begin{aligned} f^\varepsilon &= (1+m_0) M + \varepsilon \big[ (v\cdot u_0)(1+m_0) - v\cdot\nabla_x m_0 \big] M + \varepsilon g \sqrt{M}, \\ u^\varepsilon &= u_0 + \varepsilon u, \\ \rho^\varepsilon &= 1 + h_0 + \varepsilon \rho, \end{aligned}

其中余项 $(g,u,\rho)$ 属于如下函数类：

\begin{aligned} g &\in C^0([0,\infty); H^4_{x,v}) \cap L^2([0,\infty); H^5_{x,v}), \\ u &\in C^0([0,\infty); H^4_x) \cap L^2([0,\infty); H^5_x), \\ \rho &\in C^0([0,\infty); H^4_x) \cap L^2([0,\infty); H^4_x), \end{aligned}

并满足统一能量估计：

\sup_{t\ge0} E(t) + C_0 \int_0^\infty D(s)\,ds \lesssim E(0),

其中总能量泛函 $$E(t)$$ 与耗散泛函 $$D(t)$$ 定义为：

E(t) = \|g\|_{H^4_{x,v}}^2 + \|(u,\rho)\|_{H^4_x}^2 + \|m_0\|_{H^6_x}^2 + \|(u_0,h_0)\|_{H^5_x}^2,

\begin{aligned} D(t) =& \varepsilon^2 \left( \sum_{|\alpha|+|\beta|\le4} \|\partial_x^\alpha \nabla_v^\beta (I-P)g\|_\nu^2 + \|b - \varepsilon u\|_{H^4_x}^2 \right) \\ &+ \|(\nabla_x u, \operatorname{div}_x u)\|_{H^4_x}^2 + \varepsilon \|(\nabla_x b, \operatorname{div}_x b)\|_{H^3_x}^2 \\ &+ \|(\partial_t u, \nabla_x \rho, \nabla_x a)\|_{H^3_x}^2 + \|(\nabla_x m_0, \nabla_x u_0)\|_{H^5_x}^2 + \|(\partial_t \nabla_x m_0, \partial_t u_0, \nabla_x h_0, \operatorname{div}_x u_0)\|_{H^4_x}^2. \end{aligned}

作为直接推论，利用Sobolev嵌入 $H^2_{x,v} \hookrightarrow L^\infty_{x,v}$ ，得到逐点收敛速率：

Corollary 1.1：在上述条件下，

|f^\varepsilon(t,x,v) - (1+m_0(t,x)) M(v)| + |u^\varepsilon(t,x) - u_0(t,x)| + |\rho^\varepsilon(t,x) - (1+h_0(t,x))| \lesssim \varepsilon,

对所有 $(t,x,v) \in \mathbb{R}_+ \times \mathbb{R}^3 \times \mathbb{R}^3$ 成立。

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

a) 核心分析学工具

本文证明框架建立在Hilbert展开法与宏观-微观分解（macro-micro decomposition） 的双重基石之上，辅以精细的能量估计技术。

首先，区别于以往依赖相对熵方法（relative entropy method） 的工作（如Mellet-Vasseur [32]），作者采用Caflish型Hilbert展开策略：将原始解按小参数 $\varepsilon$ 展开为

f^\varepsilon = f_0 + \varepsilon f_1 + \varepsilon^2 f_2 + \cdots,\quad u^\varepsilon = u_0 + \varepsilon u_1 + \cdots,\quad \rho^\varepsilon = \rho_0 + \varepsilon \rho_1 + \cdots.

但关键创新在于，展开中心并非自由态，而是以NSS解 $(n_0,u_0,\rho_0)$ 构造的精确近似解。具体而言，作者提出如下新型展开形式：

\sqrt{f^\varepsilon} = \sqrt{M}(1 + m_0) + \varepsilon \sqrt{M} \big[ (v\cdot u_0)(1 + m_0) - v\cdot \nabla_x m_0 \big] + \varepsilon g,

即引入加权平方根变换以匹配Fokker-Planck算子的对称结构。这一构造使余项 $$g$$ 自动满足零均值条件，并极大简化线性化算子的谱结构。

其次，整个分析依赖于L^2_{x,v}-based加权Sobolev能量估计。作者定义了五类能量泛函：宏观能量 $E_{\mathrm{ma}}$ 、四类微观能量 $E_{\mathrm{mi},K,i}$ （ $$i=1,2,3,4$$ ）以及Fisher型交叉项 $E_{\mathrm{mi},F}$ ，并通过精心选取权重系数 $\lambda_i$ 构造总能量 $E(t) = \sum_{i=1}^4 \lambda_i E_{\mathrm{mi},K,i} + \lambda_5 E_{\mathrm{mi},F} + \lambda_6 E_{\mathrm{ma}}$ 。该泛函设计确保各阶导数项能被相应耗散项控制。

第三，证明中大量使用投影算子理论。记 $P = P_0 \oplus P_1$ 为 $$L^2_v$$ 上向 $\operatorname{Span}\{\sqrt{M}, v\sqrt{M}\}$ 的正交投影，其中 $P_0 g = a \sqrt{M},\; P_1 g = (v\cdot b)\sqrt{M}$ 。由此可将 $$g$$ 分解为宏观部分 $$(a,b)$$ 与微观部分 $$(I-P)g$$ 。线性化Fokker-Planck算子 $$L$$ 具有如下基本性质：

\langle L g, g \rangle_{x,v} \ge c_0 \| (I-P)g \|_\nu^2 + \|b\|_{L^2_x}^2,

该不等式是获得速度方向耗散的关键。此外，利用 $$(I-P)$$ 的自伴性与幂等性（见(1.9)-(1.10)），可在高阶导数估计中保持算子结构不变。

最后，局部适定性由压缩映射原理保证：对固定 $\varepsilon > 0$ ，将余项系统(3.2)视为关于 $(g,u,\rho)$ 的拟线性演化方程组，在适当加权Sobolev空间中构造迭代序列，并利用先验估计验证其收缩性。

b) 技术性难题与应对技巧

本文面临三大技术挑战，作者分别提出了精巧对策：

挑战一：可压缩性引入的强耦合奇异性。与不可压情形不同，连续性方程 $\partial_t \rho + \operatorname{div}(\rho u) = 0$ 与动量方程深度耦合，导致余项系统中出现高阶奇异项，如 $\frac{1}{\varepsilon^2} \nabla_x P(1+h+\varepsilon\rho) - \nabla_x P(1+h)$ 。传统Hilbert展开无法抵消此类 $\varepsilon^{-2}$ 项。

对策：重构展开形式 + 高阶Taylor展开补偿。作者观察到，若仅用标准展开，压力差将产生 $\varepsilon^{-1} \nabla_x (P'(\rho_0)\rho)$ 项，破坏能量闭合。为此，他们显式计算：

\nabla_x P(1+h+\varepsilon\rho) - \nabla_x P(1+h) = \varepsilon P'(1+h) \nabla_x \rho + \varepsilon^2 \int_0^1 P''(1+h+s\varepsilon\rho) \rho \nabla_x \rho \, ds,

并将首项纳入余项方程的线性主部，次项作为高阶小量处理。同时，新展开中 $$f_0,f_1$$ 的构造已内嵌了对 $P(\rho)$ 的一阶响应，使余项 $$g$$ 的演化方程中不再含 $\varepsilon^{-2}$ 压力源。

挑战二：速度空间奇异性与各向异性耗散。Fokker-Planck算子 $$L$$ 在 $$v$$ 方向仅提供 $\|\nabla_v (I-P)g\|_{L^2}$ 耗散，但宏观变量 $$a,b$$ 的演化需从 $$(I-P)g$$ 中提取信息。如何将 $$v$$ -导数耗散转化为 $$x$$ -导数控制？

对策：混合导数能量泛函 + Fisher信息交叉项。作者定义了含 $\partial_x^\alpha \partial_v^\beta (I-P)g$ 的耗散项，并引入Fisher型泛函：

E_{\mathrm{mi},F}(t) = \|(a,b)\|_{H^3_x}^2 + 2 \sum_{|\alpha|=0}^3 \sum_{i,j=1}^3 \langle \partial_x^\alpha (\partial_{x_i} b_j + \partial_{x_j} b_i), [\partial_x^\alpha (I-P)g](v_i v_j - \delta_{ij}) \sqrt{M} \rangle_{x,v},

该交叉项巧妙桥接了 $$v$$ -动量矩与 $$x$$ -梯度，使 $\|\nabla_x b\|$ 可由 $\|\nabla_v (I-P)g\|$ 控制，从而实现各向异性耗散的“传导”。

挑战三：非线性项的高阶匹配与余项闭合。余项方程(3.2)含大量非线性耦合项，如 $R_3 = -[\rho \partial_t u_0 + \cdots + \varepsilon^2 \rho u \cdot \nabla_x u]$ ，其结构复杂，易破坏能量不等式。

对策：分层能量估计 + 权重动态调整。作者将能量泛函分为四层： $E_{\mathrm{mi},K,1}$ 控制 $$g$$ 的 $$L^2$$ 范数与 $$u$$ 的加权 $$L^2$$ 范数； $E_{\mathrm{mi},K,2}$ 控制 $$u$$ 的时间导数； $E_{\mathrm{mi},K,3}$ 控制 $\rho$ 的时间与空间导数； $E_{\mathrm{mi},K,4}$ 控制 $$(I-P)g$$ 的高阶 $$v$$ -导数。每层估计中，通过选择足够大的权重 $\lambda_i$ ，使上一层的非线性项被下一层的耗散吸收。例如， $$R_3$$ 中的 $\varepsilon^2 \rho u \cdot \nabla_x u$ 项，其大小为 $O(\varepsilon^2 \|g\|_{H^4} \|u\|_{H^5}^2)$ ，可通过 $\lambda_2$ 与 $\lambda_4$ 的比例关系压制。

与之前工作的比较

a) 突破性改进与优势体现

本文相较既有文献实现了三重实质性突破：

第一，收敛模式的根本强化：从弱收敛到逐点收敛。此前所有rigorous hydrodynamic limit结果（如Mellet-Vasseur [32], Ballew [2], Choi-Jung [13]）均基于相对熵方法，其本质是利用解与近似解之间的模化能量（modulated energy）构造Lyapunov泛函。该方法虽能证明 $$L^1$$ 或 $$L^2$$ 意义下的收敛，但受限于熵不等式的积分特性，无法获得逐点（pointwise）或一致（uniform）收敛性。而本文通过Hilbert展开+能量估计，不仅证明了解的存在性与唯一性，更导出显式收敛速率 $\varepsilon$ ，并借助Sobolev嵌入 $H^2_{x,v} \hookrightarrow L^\infty_{x,v}$ 得到：

\sup_{t,x,v} \big| f^\varepsilon(t,x,v) - n_0(t,x) M(v) \big| \lesssim \varepsilon.

这是首个对可压缩VFP-CNS系统实现逐点收敛的rigorous结果，显著增强了极限模型的物理可解释性——例如，在数值模拟中可直接用 $$n_0 M(v)$$ 替代 $f^\varepsilon$ 而无需担心局部振荡。

第二，初始数据正则性要求的实质性降低。Mellet-Vasseur [31] 对VFP-CNS全局弱解的要求是 $L^1 \cap L\log L$ 初值，而相对熵方法通常需 $$H^1$$ 或更高；本文虽仍需 $$H^6_x$$ 宏观初值，但对微观初值 $g^{\mathrm{in}}$ 仅需 $H^4_{x,v}$ ，且该空间允许非紧支集、非小振幅的扰动（只要范数小）。更重要的是，本文的初值相容性条件（即(1.16)式）显式构造了与NSS解的最优匹配，避免了传统方法中因初值不匹配导致的初始层（initial layer）问题。

第三，模型适用范围的拓展与物理机制的深化。以往工作多聚焦于不可压流体（如Goudon-Jabin-Vasseur [22]）或无压力流体（如Carrillo-Goudon [11]），而本文首次在完全可压缩框架下处理小De极限，完整保留了压力项 $P(\rho) = A\rho^\gamma$ 的非线性效应，并揭示了其如何通过 $\nabla_x P(\rho)$ 项影响颗粒-流体动量交换。此外，作者明确区分了“气泡区（bubbling）”与“流动区（flowing）”，指出当前工作为后续攻克更具挑战性的流动区极限（需处理激波、真空等奇性）提供了关键技术模板，如Hilbert展开的重构思想、混合导数泛函的设计等。

综上，本文不仅填补了可压缩流体-颗粒系统小De极限理论的空白，更通过方法论创新，为多尺度偏微分方程的奇异极限研究树立了新的技术标杆。

2. 带加性噪声的随机Navier–Stokes方程的高阶离散格式

A Higher Order Discretization for the Stochastic Navier–Stokes equations with additive Noise

作者: L. Banas, D. Breit, A. Chaudhary, A. Prohl

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

发布日期: 2026-02-15

arXiv ID: 2602.14268v1

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

详细分析

中文标题: 带加性噪声的随机Navier–Stokes方程的高阶离散格式

研究对象和背景

本文研究的是二维环面 $\mathbb{T}^2 \subset \mathbb{R}^2$ 上带加性噪声的随机Navier–Stokes方程（stochastic Navier–Stokes equations, SNSE）, 其原始形式为如下耦合系统：

\begin{cases} du(t) + \big[(u(t)\cdot\nabla)u(t) - \nu\Delta u(t) + \nabla p(t)\big]\,dt = \Phi\,dW(t), & \text{in } Q_T = (0,T) \times \mathbb{T}^2, \\ \operatorname{div} u = 0, & \text{in } Q_T, \\ u(0) = u_0, & \text{in } \mathbb{T}^2, \end{cases}

其中：

$u(t,x) \in \mathbb{R}^2$ 是速度场, $p(t,x) \in \mathbb{R}$ 是压力函数；
$\nu > 0$ 是运动粘性系数；
$$ W(t) $$ 是定义在完备概率空间 $(\Omega,\mathcal{F},(\mathcal{F}_t)_{t\geq 0},\mathbb{P})$ 上的柱状Wiener过程, 即
$W(t) = \sum_{k\geq 1} W_k(t) e_k,$ 其中 $\{e_k\}_{k\geq 1}$ 是可分Hilbert空间 $$ U $$ 的标准正交基， $\{W_k\}_{k\geq 1}$ 是相互独立的标准布朗运动；
加性噪声由线性算子 $\Phi \in \mathcal{L}_2(U; H^2(\mathbb{T}^2)^2)$ 驱动，且满足无散度条件（solenoidal）：
$\Phi W(t,x) := \sum_{k\geq 1} W_k(t) \phi_k(x), \quad \phi_k \in H^2(\mathbb{T}^2)^2,\ \operatorname{div}\phi_k = 0.$

该模型描述了在随机外力扰动下不可压缩流体的演化行为。与乘性噪声不同，加性噪声不依赖于解本身，因此在数学结构上更“友好”，但仍因Wiener过程的几乎必然仅属于 $C^{1/2^-}$ 类（Hölder连续但不可微） 而带来严重的时间正则性障碍。

论文特别强调其分析框架建立在周期边界条件（即定义在 $\mathbb{T}^2$ 上）与函数空间的强正则性假设之上，这是为了充分利用确定性Navier–Stokes理论中的经典正则性结果（如Ladyzhenskaya不等式、Poincaré不等式、Sobolev嵌入等）。关键函数空间包括：

$H := \{ v \in L^2(\mathbb{T}^2)^2 : \operatorname{div} v = 0 \}$ ：零散度平方可积向量场空间；
$V := H \cap H^1(\mathbb{T}^2)^2$ ：零散度H¹向量场空间，其范数为 $\|v\|_V = \|\nabla v\|_{L^2}$ ；
$Q := \{ q \in L^2(\mathbb{T}^2) : \int_{\mathbb{T}^2} q\,dx = 0 \}$ ：零均值压力空间；
$\mathcal{L}_2(U;X)$ ：从 $$ U $$ 到 Hilbert 空间 $$ X $$ 的Hilbert–Schmidt算子空间。

值得注意的是，论文明确指出其理论分析严格限定在二维环面 $\mathbb{T}^2$ ，以规避复杂边界带来的正则性损失，并便于使用Fourier方法与谱投影技术（如Helmholtz–Leray投影 $P_{\mathrm{HL}}$ ）。尽管数值实现可推广至有界域（如单位正方形 $$ (0,1)^2 $$ ）并采用Dirichlet边界条件（见第5节），但收敛性证明所依赖的先验估计（如引理2.2与2.3）本质上依赖于周期性与谱结构。

此外，作者强调：由于加性噪声的特殊结构，可通过随机变换（stochastic transformation） 将原SPDE转化为一个含随机系数的确定性偏微分方程（random PDE），这是本工作的核心思想起点。具体地，令

y(t) := u(t) - \Phi W(t),

则 $$ y $$ 满足如下随机PDE系统（见(1.3)或(2.7)）：

\begin{cases} \partial_t y + \big[(y+\Phi W)\cdot\nabla\big](y+\Phi W) - \nu\Delta(y+\Phi W) + \nabla p = 0, & \text{in } Q_T, \\ \operatorname{div} y = 0, & \text{in } Q_T, \\ y(0) = u_0, & \text{in } \mathbb{T}^2. \end{cases}

该变换的关键优势在于：虽然 $$ u $$ 本身仅具有 $C^{1/2^-}$ 时间正则性，但 $$ y $$ 在适当条件下可达到一阶时间可微性（即 $\partial_t y \in L^2(0,T;H)$ 或 $$ L^2(0,T;V) $$ ），从而为构造高阶时间离散格式提供了理论基础。这也解释了为何传统Euler型方法对SNSE仅能达到 $O(\sqrt{\tau})$ 强收敛率，而本文目标是突破至 $O(\tau^{1.5})$ 。

主要定理或主要结果

本文的核心理论成果是以下强收敛性定理，它确立了所提出的Modified Crank–Nicolson格式（Algorithm 1）对速度与压力的高阶逼近精度。

Theorem 3.1（速度场强收敛率，原文Section 3.3）
设 $$ u $$ 是SNSE (1.1) 在Definition 2.1意义下的唯一弱路径解， $\{u_n\}_{n=0}^N$ 是Algorithm 1生成的离散速度序列。假设初始数据与噪声满足：
$u_0 \in L^4(\Omega; H^4(\mathbb{T}^2)^2) \cap L^{36}(\Omega; V), \quad \Phi \in \mathcal{L}_2(U; H^4(\mathbb{T}^2)^2).$

则存在常数 $$ C > 0 $$ ，使得对任意足够小的时间步长 $\tau$ ，有
$\mathbb{E}\Bigg[ \max_{0 \le n \le N} \big\| u(t_n) - u_n \big\|_{L^2(\mathbb{T}^2)^2}^2 \Bigg]^{1/2} \le C \tau^{3/2}.$

即，速度迭代 $$ u_n $$ 以强收敛阶 $\frac{3}{2}$ 逼近真解 $$ u(t_n) $$ 。

Theorem 3.2（压力场强收敛率，原文Section 3.4）
在相同假设下，设 $\{p_n\}_{n=1}^N$ 为Algorithm 1输出的压力迭代，定义局部时间平均压力
$\bar{p}_n := \frac{1}{\tau} \int_{t_{n-1}}^{t_n} p(s)\,ds.$

则存在常数 $$ C > 0 $$ ，使得
$\mathbb{E}\Bigg[ \max_{1 \le n \le N} \big\| \bar{p}_n - p_n \big\|_{Q}^2 \Bigg]^{1/2} \le C \tau^{3/2}.$

即，压力迭代 $$ p_n $$ 以强收敛阶 $\frac{3}{2}$ 逼近真解压力的局部时间平均值 $\bar{p}_n$ 。

这两个定理共同构成了本文最核心的理论贡献，首次将非线性随机Navier–Stokes方程的时间离散格式的强收敛阶从经典的 $O(\tau)$ （隐式Euler）或 $O(\sqrt{\tau})$ （显式Euler）提升至 $O(\tau^{3/2})$ ，且适用于真实物理相关的加性噪声设定。

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

a) 核心分析学工具

本文的误差分析建立在一套精密的随机解析框架之上，融合了确定性PDE理论、随机分析与数值分析三重工具：

Helmholtz–Leray 投影与变分结构：
利用正交投影 $P_{\mathrm{HL}}: L^2(\mathbb{T}^2)^2 \to H$ 消除压力项，将动量方程写为纯速度演化形式：
$\partial_t y + P_{\mathrm{HL}}\big[(y+\Phi W)\cdot\nabla\big](y+\Phi W) - \nu \Delta y = f_{\mathrm{noise}},$
其中 $f_{\mathrm{noise}} = \nu \Delta(\Phi W) - P_{\mathrm{HL}}[L_W(y)]$ 含有随机源项。该结构保证了能量恒等式的可行性。
改进的时间正则性估计（Lemma 2.3）：
关键引理证明了变换后变量 $$ y $$ 具有比 $$ u $$ 更优的时变性质。例如，在足够正则的 $\Phi$ 和 $$ u_0 $$ 下，
$\partial_t y \in L^2(0,T; V), \quad \text{a.s.},$
且满足先验界：
$\int_0^T \|\partial_t \nabla y\|_{L^2}^2\,dt \le C\left( \|\Phi W\|_{L^\infty_t H^3_x}^2 + \|u\|_{L^\infty_t H^3_x}^2 + \cdots \right).$
此估计是构造高阶格式的先决条件——若 $$ y $$ 不具 $$ C^1 $$ 时间正则性，则Crank–Nicolson无法获得超收敛。
多尺度网格与高阶数值积分：
为匹配 $$ y $$ 的高正则性与 $\Phi W$ 的低正则性，算法采用异构离散策略：
- 对光滑部分（如 $\nu \Delta y$ ）使用Crank–Nicolson格式（二阶精度）；
- 对粗糙部分（如 $\Delta(\Phi W)$ ）使用精细网格上的梯形/矩形积分（微网格步长 $\tau^2$ ），使其截断误差达 $O(\tau^3)$ ，从而不破坏整体 $O(\tau^{3/2})$ 阶。
截断与扰动论证（truncation and perturbation argument）：
由于非线性项导致离散格式缺乏直接的能量守恒，作者引入高概率集
$\Omega_R := \left\{ \omega \in \Omega : \sup_{t \in [0,T]} \|y(t,\omega) + \Phi W(t,\omega)\|_{H^2} \le R \right\},$
并在其上进行一致误差估计，再利用切比雪夫不等式控制补集贡献。这是处理非线性SPDE高阶方法缺乏离散能量泛函这一根本困难的标准技术。

b) 处理非线性与奇异性技巧

非线性项的线性化处理：
在离散格式(3.3)中，对对流项 $C(a,b,c) = \int (a\cdot\nabla)b \cdot c$ ，取
$a = y^\star_{n+1/2} + \Phi I^n_W, \quad b = y_{n+1/2} + \Phi I^n_W, \quad c = \varphi,$
其中 $y^\star_{n+1/2}$ 是BDF2外推， $y_{n+1/2}$ 是中点平均。此设计使格式关于 $y_{n+1}$ 线性，避免了全隐式牛顿迭代，同时保留了非线性结构的几何特性（如能量守恒近似）。
噪声项的三重积分近似（ $I^n_{W^2}$ ）：
为精确捕捉 $\int (y+\Phi W) \otimes (y+\Phi W)$ 中的交叉项 $\Phi W \otimes \Phi W$ ，定义矩阵值三重积分
$Q^n_{W^2} := \frac{1}{\tau^3} \int_{t_n}^{t_{n+1}} \int_{t_n}^{t_{n+1}} \int_{t_n}^{t_{n+1}} \big[\Phi W(t)-\Phi W(s)\big] \otimes \big[\Phi W(t)-\Phi W(r)\big]\, ds\,dt\,dr,$
并在 $\tau^2$ -细网上用矩形法则逼近。其均方误差满足
$\mathbb{E}\big[ \|Q^n_{W^2} - I^n_{W^2}\|_{L^2}^2 \big] \le C \tau^3,$
这是达成 $\tau^{3/2}$ 整体阶的关键一环。
压力重构与 inf-sup 稳定性：
为处理压力项，算法采用混合有限元离散（如Taylor–Hood元），并借助Babuska–Brezzi理论确保离散inf-sup条件成立。压力收敛性并非指向瞬时值 $$ p(t_n) $$ （因其缺乏时间正则性），而是指向局部时间平均 $\bar{p}_n$ ，这与 $\partial_t p$ 不存在的事实完全兼容。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重实质性突破：

收敛阶的首次提升：
此前，对SNSE的强收敛分析几乎全部止步于 $O(\tau)$ （如[5]中的隐式Euler）或 $O(\sqrt{\tau})$ （显式Euler）。本文是首个为非线性SNSE建立 $O(\tau^{3/2})$ 强收敛率的完整理论框架。该结果直接回应了[12]中线性热方程高阶格式能否推广至Navier–Stokes的开放问题。
正则性假设的精巧平衡：
虽然要求 $u_0 \in H^4$ 、 $\Phi \in H^4$ ，看似苛刻，但作者明确指出：这些是为补偿高阶方法缺失离散能量泛函而引入的技术性假设（见Section 3引言）。相比而言，低阶方法（如隐式Euler）可在 $u_0 \in V$ 、 $\Phi \in \mathcal{L}_2(U;H)$ 下工作。本文的改进在于：以适度增强的数据正则性为代价，换取计算效率的显著提升（允许更大 $\tau$ ，见Figure 1与Section 5讨论）。
方法普适性的拓展：
文献[12]仅处理线性方程；本文首次将随机变换+多尺度积分+非线性线性化策略系统应用于非线性、非自治、强耦合的流体力学SPDE。此外，Section 3.5与4指出：该框架可自然延拓至一般噪声（非solenoidal） 及Stokes极限（忽略对流项），显示其方法论具有宽广适应性。

综上，本文不仅是一项技术性突破，更标志着随机流体数值分析从“可用”迈向“高效可靠” 的重要一步。

3. 具有诺伊曼边界条件的二维外区域上二次拟线性波动方程的长时间光滑解

Long time smooth solutions of 2-D quadratic quasilinear wave equations in exterior domains with Neumann boundary conditions

作者: Fei Hou, Huicheng Yin, Meng Yuan

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14019v1

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

详细分析

中文标题: 具有诺伊曼边界条件的二维外区域上二次拟线性波动方程的长时间光滑解

研究对象和背景

本文研究的是定义在 $\mathbb{R}^2$ 中一个紧致障碍物 $O \subset \mathbb{R}^2$ 的外部区域 $K = \mathbb{R}^2 \setminus O$ 上、满足齐次诺伊曼边界条件的二维二次拟线性波动方程的初边值问题（IBVP）。该问题具有深刻的物理动机，其解描述了在不可渗透刚性障碍物存在下、低振幅扰动在二维空间中传播的长期行为。

设 $t \geq 0$ , $x = (x_1, x_2) \in K$ , 记 $\partial = (\partial_t, \partial_{x_1}, \partial_{x_2}) = (\partial_0, \partial_1, \partial_2)$ , 波算子为 $\Box = \partial_t^2 - \Delta$ , 其中 $\Delta = \partial_{x_1}^2 + \partial_{x_2}^2$ 。障碍物 $$O$$ 是凸的、光滑的紧集，且包含原点；其边界 $\partial K = \partial O$ 是光滑闭曲线；单位外法向量记为 $$n = (0, n_1(x), n_2(x))$$ 。初始数据 $(u_0, u_1) \in C^\infty(K)$ 满足 $\operatorname{supp}(u_0, u_1) \subset \{x \in K : |x| \le M_0\}$ ，即具有紧支集。

所研究的主方程为：

\begin{cases} \Box u = Q(\partial u, \partial^2 u), & (t,x) \in [0,+\infty)\times K, \\ \displaystyle \frac{\partial u}{\partial n}(t,x) = 0, & (t,x) \in [0,+\infty)\times \partial K, \\ (u, \partial_t u)(0,x) = (\varepsilon u_0, \varepsilon u_1)(x), & x \in K. \end{cases}

其中非线性项 $Q(\partial u, \partial^2 u)$ 是关于一阶与二阶导数的二次齐次多项式，形如：

Q(\partial u, \partial^2 u) = \sum_{\alpha,\beta,\mu=0}^{2} Q_{\alpha\beta\mu}\, \partial_\alpha \partial_\beta u \, \partial_\mu u + \sum_{\alpha,\beta,\mu,\nu=0}^{2} Q_{\alpha\beta\mu\nu}\, \partial_\alpha \partial_\beta u \, \partial_\mu u \, \partial_\nu u,

系数满足对称性 $Q_{\alpha\beta\mu} = Q_{\beta\alpha\mu}$ , $Q_{\alpha\beta\mu\nu} = Q_{\beta\alpha\mu\nu} = Q_{\alpha\beta\nu\mu}$ 。为保证诺伊曼边界条件与非线性结构相容，作者引入可容许性条件（admissible condition）：对任意在边界上满足 $\partial_n w = \partial_n \psi = 0$ 的 $$C^1$$ 函数 $w,\psi$ ，要求

\sum_{\alpha,\beta,\mu=0}^{2} Q_{\alpha\beta\mu}\, n_\alpha \, \partial_\beta w \, \partial_\mu \psi = 0, \quad \sum_{\alpha,\beta,\mu,\nu=0}^{2} Q_{\alpha\beta\mu\nu}\, n_\alpha \, \partial_\beta w \, \partial_\mu \psi \, \partial_\nu \psi = 0,

在 $[0,+\infty)\times \partial K$ 上恒成立。该条件确保非线性项在边界处不破坏诺伊曼约束，是构造高阶能量估计的前提。

此外，为使解在时间方向保持光滑，需满足相容性条件（compatibility conditions）：令 $J_k u = \{ \partial^\alpha u : |\alpha| \le k \}$ ，并由方程与初始数据递推定义 $\partial_t^k u(0,x) = G_k(J_k u_0, J_{k-1} u_1)$ ，则要求对所有 $0 \le k \le 2N$ ，均有 $\partial_n G_k \equiv 0$ 在 $\partial K$ 上成立。

函数空间设定为加权 Sobolev 类：初始数据属于高阶 Sobolev 空间 $H^{2N+1}(K) \times H^{2N}(K)$ ，解 $$u$$ 落在时空函数类 $\bigcap_{j=0}^{2N+1} C^j([0,T], H^{2N+1-j}(K))$ ，其中 $T = \varepsilon^{\delta-2}$ 或 $T = +\infty$ ，取决于非线性结构是否满足更强的零条件（null condition）。

特别地，在 Theorem 1.2 中，非线性项进一步被限制为含 $$Q_0$$ 型零形式：

Q(\partial u, \partial^2 u) = 2 \sum_{\alpha=0}^{2} C_\alpha Q_0(\partial_\alpha u, u) + \sum_{\alpha,\beta,\mu,\nu=0}^{2} Q_{\alpha\beta\mu\nu}\, \partial_\alpha \partial_\beta u \, \partial_\mu u \, \partial_\nu u,

其中 $Q_0(f,h) = \partial_t f \, \partial_t h - \sum_{i=1}^{2} \partial_i f \, \partial_i h$ ，并附加第二零条件：

\sum_{\alpha,\beta,\mu,\nu=0}^{2} Q_{\alpha\beta\mu\nu}\, \xi_\alpha \xi_\beta \xi_\mu \xi_\nu \equiv 0, \quad \forall (\xi_0,\xi_1,\xi_2) \in \{\pm1\} \times S^1.

该条件对应于“强衰减”非线性，是获得全局解的关键代数结构。

主要定理或主要结果

Theorem 1.1（长时间存在性）
假设障碍物 $$O$$ 是凸的，且非线性项 $$Q$$ 满足可容许性条件 (1.3) 与相容性条件 (1.4)。则对任意固定 $\delta \in (0,1)$ ，存在 $\varepsilon_0 > 0$ ，使得当 $0 < \varepsilon \le \varepsilon_0$ 且初始数据满足

\|u_0\|_{H^{2N+1}(K)} + \|u_1\|_{H^{2N}(K)} \le 1, \quad N \ge 19,

时，初边值问题 (1.1) 存在唯一的光滑解

u \in \bigcap_{j=0}^{2N+1} C^j\big([0,\varepsilon^{\delta-2}],\, H^{2N+1-j}(K)\big).

Theorem 1.2（全局存在性）
在 Theorem 1.1 假设基础上，若 $$Q$$ 具有 $$Q_0$$ -型零形式结构 (1.6)，并满足第二零条件 (1.7)，则存在 $\varepsilon_0 > 0$ ，使得当 $0 < \varepsilon \le \varepsilon_0$ 且

\|u_0\|_{H^{2N+1}(K)} + \|u_1\|_{H^{2N}(K)} \le 1, \quad N \ge 40,

时，问题 (1.1) 存在全局光滑解

u \in \bigcap_{j=0}^{2N+1} C^j\big([0,\infty),\, H^{2N+1-j}(K)\big),

且满足如下时空点态衰减估计：

\sum_{|a|\le N+1} |\partial Z^a u| \le C\varepsilon \langle x\rangle^{-1/2}\langle t-|x|\rangle^{-1}\ln^2(2+t),

\sum_{|a|\le N+1}\sum_{i=1}^{2} |\overline{\partial}_i Z^a u| \le C\varepsilon \langle x\rangle^{-1/2}\langle t+|x|\rangle^{-0.999},

\sum_{|a|\le N+1} |Z^a u| \le C\varepsilon \langle t+|x|\rangle^{-1/2}\langle t-|x|\rangle^{-0.499},

|\partial_t \partial_t u| \le C\varepsilon \langle x\rangle^{-1/2}\langle t-|x|\rangle^{-1},

其中 $Z = \{\partial, \Omega\}$ ， $\Omega = x_1 \partial_2 - x_2 \partial_1$ 为角动量向量场， $\langle x\rangle = \sqrt{1+|x|^2}$ ， $\overline{\partial}_i = \frac{x_i}{|x|}\partial_t + \partial_i$ 为好导数（tangent to outgoing light cone），以及如下局部能量衰减估计：

\sum_{|a|\le 2N-29} \|\partial^a u\|_{L^2(K_R)} \le C_R \varepsilon (1+t)^{-1}\ln(2+t),

\sum_{|a|\le 2N-30} \|\partial^a \partial_t u\|_{L^2(K_R)} \le C_R \varepsilon (1+t)^{-1},

其中 $K_R = K \cap \{ |x| \le R \}$ ， $$R > 1$$ 固定。

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

a) 核心分析学工具

本文采用基于向量场方法的能量估计框架，结合加权时空衰减估计与精细的边界处理技术。其核心工具包括：

向量场方法（Vector Field Method）：引入完整向量场族 $Z = \{\partial_t, \partial_{x_1}, \partial_{x_2}, \Omega\}$ ，其中 $\Omega$ 刻画旋转对称性；通过交换 $$Z^a$$ 与 $\Box$ ，建立高阶导数的控制。由于外区域缺乏平移不变性，需借助截断函数与局部化技巧处理边界附近区域。
加权能量泛函（Weighted Energy Estimates）：定义带权重 $\langle t-r\rangle^p \langle t+r\rangle^q$ 的能量，例如
$E_{N}(t) = \sum_{|a|\le N} \int_K \left[ (\partial_t Z^a u)^2 + |\nabla Z^a u|^2 \right] \langle x\rangle^{2\sigma} dx,$
其中 $\sigma > 0$ 用于增强近边界区域的控制能力。该加权设计直接应对诺伊曼边界下缺乏 Poincaré 不等式的困难。
Good Derivatives 与 Null Structure 的耦合分析：定义“好导数” $\overline{\partial}_i = \frac{x_i}{|x|}\partial_t + \partial_i$ ，它们沿光锥 $$|x|=t$$ 切向，具有比普通导数更优的衰减性质。论文证明：在零条件下， $$Q$$ 中含 $\overline{\partial}$ 的项主导衰减，而坏导数 $\partial_t - \frac{x_i}{|x|}\partial_i$ （径向导数）仅以对数修正形式出现，从而控制 $\ln^2(2+t)$ 因子。
时空点态估计（Pointwise Spacetime Decay）：基于二维外区域上Neumann 波方程的新型 Green 函数估计与辐射场表示，建立如下关键引理：
$|\partial Z^a u(t,x)| \lesssim \frac{1}{\langle x\rangle^{1/2}} \left( \frac{1}{\langle t-|x|\rangle} + \frac{\ln(2+t)}{\langle t+|x|\rangle} \right) \cdot \text{(weighted energy)},$
此估计比经典自由波估计 $O(\langle t-|x|\rangle^{-1})$ 多出 $\ln^2$ 修正，源于 Neumann 边界反射导致的慢速能量滞留。
Gronwall-type Bootstrap Argument：设定先验假设（bootstrap assumption）
$\sum_{|a|\le N+1} |\partial Z^a u| \le A \varepsilon \langle x\rangle^{-1/2}\langle t-|x|\rangle^{-1}\ln^2(2+t),$
并利用上述能量与点态估计闭环验证其自洽性，最终确定 $$A$$ 可取为绝对常数。

b) 关键技术性突破

“好未知量”（Good Unknowns）的系统构造：针对 Neumann 边界下无法直接控制 $\partial_n u$ 的困难，作者引入一系列新变量：
$w = \partial_t u - \frac{x_i}{|x|}\partial_i u, \quad v = \partial_i u - \frac{x_i}{|x|}\partial_r u, \quad \text{etc.},$
这些变量自动满足齐次 Dirichlet 型边界条件或具有更高阶衰减，从而绕过法向导数奇异性。
Neumann 问题的新型局部能量衰减引理：证明了对任意 $$R>1$$ ，
$\int_0^t \int_{K_R} |\partial u(s,x)|^2 dx ds \lesssim \varepsilon^2 \ln(2+t),$
该对数增长是二维 Neumann 外问题的本质特征，区别于 Dirichlet 情形下的指数衰减或三维情形的代数衰减。
非线性项的精细分解与零条件利用：将 $$Q$$ 分解为“坏-坏”、“坏-好”、“好-好”三类组合，并证明：仅“坏-坏”项贡献主导对数因子，其余项可被吸收进能量闭包；第二零条件则彻底消除四次坏导数组合，保障全局存在。
截断与 Duhamel 原理的协同使用：将解写为 $u = u_{\text{lin}} + \int_0^t \mathcal{K}(t-s) Q(\cdots) ds$ ，其中 $\mathcal{K}$ 为 Neumann 波传播子；利用外区域上的局部光滑效应（Local Smoothing Effect） 控制积分核的 $L^2 \to L^\infty$ 映射范数。

与之前工作的比较

a) 突破性改进

填补二维 Neumann 外问题的长期存在性空白：此前，三维情形（ $$n=3$$ ）的全局/长时间存在性已由 Metcalfe–Sogge、Wang–Yu 等人解决；二维 Dirichlet 情形由 authors 自身在 [17] 中完成；但二维 Neumann 情形自 1990 年代以来始终未解，本文首次给出 $\varepsilon^{\delta-2}$ 时间存在性（Thm 1.1）及全局存在性（Thm 1.2），属领域内里程碑式突破。
放宽对初始数据正则性与支集的要求：相比早期工作（如 [29,39] 对线性 Neumann 波的 $$H^3$$ 要求），本文仅需 $H^{2N+1} \times H^{2N}$ ，且允许任意紧支集（不限于星形区域），关键在于向量场方法与加权能量对障碍物几何的鲁棒性。
揭示二维 Neumann 边界的本质对数障碍：以往猜测二维 Neumann 解寿命仅为 $O(\varepsilon^{-2})$ ，本文不仅证实此下界，更精确刻画其为 $\varepsilon^{\delta-2}$ （任意 $\delta<1$ ），并指出 $\ln^2(2+t)$ 修正不可去除——这是由二维 Green 函数在边界附近的对数奇异性决定的，属维度特有现象。
应用广度显著拓展：结果直接适用于：
- 二维等熵无旋可压 Euler 方程（polytropic 与 Chaplygin 气体），
- 二维相对论膜方程（relativistic membrane equation），
- 二维经典膜方程（membrane equation），所有这些模型在物理中描述弹性膜振动、宇宙弦动力学等，此前均无二维外区域 Neumann 全局理论。

综上，本文不仅解决了长期悬而未决的基础 PDE 问题，更发展了一套适用于低维、非紧、带反射边界的拟线性波动系统的通用分析范式，为后续研究高维 Neumann 问题、非均匀介质中的波传播等开辟了新路径。

4. 一维磁流体动力学模型的稳定性与不稳定性

Stability and instability of a one-dimensional MHD model

作者: Nicola De Nitti, Jie Guo, Quansen Jiu

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.13995v1

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

详细分析

中文标题: 一维磁流体动力学模型的稳定性与不稳定性

研究对象和背景

本文研究的是一类一维磁流体动力学（MHD）模型，其原始系统由 Dai 等人于 2023 年提出，记为 (1.1)–(1.2)，定义在一维环面 $\mathbb{T} := [-\pi,\pi]$ 上，具有周期边界条件。该模型旨在捕捉三维理想不可压 MHD 系统中关键的涡量拉伸（vortex stretching）与磁场耦合效应，同时保持数学可处理性。

系统具体形式如下：

\begin{cases} \partial_t \omega^+ + a u^- \partial_\theta \omega^+ = (\omega^+ H \omega^- + p \omega^- H \omega^+) - q (\omega^+ H \omega^-), \\ \partial_t \omega^- + a u^+ \partial_\theta \omega^- = (\omega^- H \omega^+ + p \omega^+ H \omega^-) - q (\omega^- H \omega^+), \\ \partial_\theta u^\pm = H \omega^\pm, \quad u^\pm(t,0) = 0, \end{cases}

其中 $a,p,q \in \mathbb{R}$ 是参数； $$H$$ 是定义在 $\mathbb{T}$ 上的 Hilbert 变换，即

H f(\theta) := \mathrm{p.v.}\, \frac{1}{2\pi} \int_{-\pi}^{\pi} \cot\left( \frac{\theta - \vartheta}{2} \right) f(\vartheta)\, d\vartheta;

而 $u^\pm$ 通过 Biot–Savart 型关系由 $\omega^\pm$ 确定： $\partial_\theta u^\pm = H \omega^\pm$ ，且满足零均值归一化条件 $u^\pm(t,0)=0$ 。

特别地，论文聚焦于参数取值

a = 1,\quad p + q = 1,\quad p,q \geq 0,

这一设定使得当 $\omega^+ \equiv \omega^-$ （即无磁场情形， $b \equiv 0$ ）时，系统退化为著名的 De Gregorio 模型（DG 模型）：

\partial_t \omega + u \partial_\theta \omega = \omega H \omega,\quad \partial_\theta u = H \omega.

进一步，当 $\omega^\pm$ 接近第一激发态（first excited state）

\Omega^\pm(\theta) := -\sin(2\theta),\quad U^\pm(\theta) := \frac{1}{2}\sin(2\theta),

时，作者引入扰动变量：

\omega^\pm(t,\theta) = -\sin(2\theta) + \eta^+(t,\theta) \pm \eta^-(t,\theta),\quad u^\pm(t,\theta) = \frac{1}{2}\sin(2\theta) + v^+(t,\theta) \pm v^-(t,\theta),

并导出扰动满足的演化方程组（2.1）–（2.3），其线性化主部为：

\begin{cases} \partial_t \eta^+ = \frac{1}{2} \big\{ \eta^+, \sin(2\theta) \big\} + \big\{ v^+, \sin(2\theta) \big\}, \\ \partial_t \eta^- = \frac{1}{2} \big\{ \eta^-, \sin(2\theta) \big\} - \big\{ v^-, \sin(2\theta) \big\} - 2q \big( \sin(2\theta) H \eta^- + \cos(2\theta) \eta^- \big), \end{cases}

其中 $\{f,g\} := f \partial_\theta g - g \partial_\theta f$ 为 Lie 括号。

值得注意的是，该问题的研究区域是紧致流形 $\mathbb{T}$ ，而非全空间 $\mathbb{R}$ ；初始数据被限定在一类特殊的加权 Sobolev 空间 $$H_2$$ 中：

H_2 := \left\{ \eta \in H^1(\mathbb{T}) : \eta \text{ 是奇函数},\ \int_{\mathbb{T}} \frac{|\partial_\theta \eta|^2}{\sin^2 \theta}\, d\theta < \infty \right\},

其内积定义为

\langle \xi, \eta \rangle_{H_2} := \int_{-\pi}^{\pi} \rho(\theta)\, \partial_\theta \xi\, \partial_\theta \eta\, d\theta,\quad \rho(\theta) := \frac{1}{4\pi \sin^2 \theta}.

该权重 $\rho(\theta)$ 在 $\theta = 0,\pm\pi$ 处具有二次奇异性，反映了扰动在临界点附近的强局部行为，也使得标准 Sobolev 嵌入失效，必须构造适配的正交基 $\{e_{2,k}\}_{k \ge 1}$ ，其中

e_{2,k}(\theta) := \frac{\sin((k+2)\theta)}{k+2} - \frac{\sin(k\theta)}{k}.

此空间选择并非任意，而是源于对 DG 模型激发态稳定性分析的自然推广（见 [LLR20, GJ25]），它能精确刻画解在 $\theta=0$ 附近因 Hilbert 变换与 $\sin(2\theta)$ 结构引发的非均匀衰减特性。

综上，本工作的核心对象是一个带奇异权重、定义在环面上的一维双分量非线性非局部 PDE 系统，其物理动机是理解三维 MHD 中磁场如何调制涡量动力学，数学动机则是探索非线性竞争机制（输运 vs 涡量拉伸）在激发态邻域的精细平衡结构。

主要定理或主要结果

论文建立了关于该一维 MHD 模型在第一激发态 $\Omega^\pm = -\sin(2\theta)$ 附近的一系列严格稳定性/不稳定性结论，分为线性与非线性两个层面。

线性化系统的全局适定性与不稳定性：

Theorem 3.1（线性化问题的整体存在唯一性）：设 $\rho^{1/2} \partial_\theta \eta_0^\pm \in H^m(\mathbb{T})$ ， $$m > 3$$ 为整数，则对任意 $$T > 0$$ ，线性化系统 (2.7) 存在唯一经典奇解 $\eta^\pm$ ，满足
$\rho^{1/2} \partial_\theta \eta^\pm \in C([0,T]; H^m(\mathbb{T})) \cap C^2([0,T]; H^{m-2}(\mathbb{T})).$
Theorem 3.2（更一般初值下的线性适定性）：若 $\eta_0^\pm \in H^m$ ， $$m > 3$$ ，则同样存在唯一经典解 $\eta^\pm \in C([0,T]; H^m) \cap C^2([0,T]; H^{m-2})$ 。
Theorem 3.3（线性不稳定性）：设初值满足 $\|\eta_0^+\|_{H_2} \neq 0$ 且 $\langle L^+ \eta_0^+, \eta_0^+ \rangle_{H_2} \ge 0$ ，其中 $$L^+$$ 是主线性算子（见 (2.5)），则存在绝对常数 $\lambda_1 > 0$ ，使得对应解满足指数下界
$\|\eta^+(t)\|_{H_2} + \|\eta^-(t)\|_{H_2} > E_1^{1/2}(t) > 0,\quad \forall t > 0,$
其中
$E_1(t) = \|\eta_0^+\|_{H_2}^2 + \frac{1}{\sqrt{\lambda_1}} \langle L^+ \eta_0^+, \eta_0^+ \rangle_{H_2} e^{2\sqrt{\lambda_1}t} + \left( \|\eta_0^+\|_{H_2}^2 - \frac{1}{\sqrt{\lambda_1}} \langle L^+ \eta_0^+, \eta_0^+ \rangle_{H_2} \right) e^{-2\sqrt{\lambda_1}t}.$

非线性系统的局部适定性与不稳定性/稳定性：

Theorem 3.6（非线性问题的局部适定性）：设 $$m > 3$$ ，初值满足 $I_m(0) := \|\rho^{1/2} \partial_\theta \eta_0^+\|_{H^m}^2 + \|\rho^{1/2} \partial_\theta \eta_0^-\|_{H^m}^2 = \delta^2$ ，则存在 $$T_0 > 0$$ 和唯一经典奇解 $\eta^\pm$ ，满足
$\rho^{1/2} \partial_\theta \eta^\pm \in C([0,T_0]; H^m) \cap C^2([0,T_0]; H^{m-2}),$
且解的尺度化变量 $u^\pm := \rho^{1/2} \partial_\theta \eta^\pm$ 满足一致小能量控制： $\sup_{t \in [0,T_0]} \|u^\pm(t)\|_{H^m} \le C(T_0)\delta$ 。
Theorem 3.7（非线性不稳定性）：对任意 Lipschitz 函数 $F(y) \le K y$ ，存在初值 $\eta_0^\pm$ 满足 $I_m(0) < \delta$ 且 $0 \le \langle L^+ u_0^+, u_0^+ \rangle_{L^2} \le \lambda_1 \|u_0^+\|_{L^2}^2$ ，使得对应解在时间 $t_K \in (0,T_0)$ 达到
$\|u^+(t_K)\|_{L^2} > F(I_m(0)).$
这表明系统在 Lipschitz 拓扑下不稳定。
Theorem 3.10（非线性稳定性）：若初值形如 $\eta_0^+(\theta) = \sum_{k \ge 1} \eta_{0,2k}^+ e_{2,2k}(\theta)$ （即仅含偶数下标基函数分量），且 $\eta_0^- \in H_2$ ，则存在 $\varepsilon > 0$ ，当 $I_2(0) \le \varepsilon$ 时，非线性系统 (2.7) 全局适定，且满足指数衰减
$I_2(t) := \|\eta^+(t)\|_{H_2}^2 + \|\eta^-(t)\|_{H_2}^2 \lesssim e^{-t/2} I_2(0),\quad \forall t \ge 0.$

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

a) 核心分析学工具

论文综合运用了多种现代 PDE 分析技术，构成一个严密的逻辑链条：

Galerkin 逼近法：为证明线性化系统的整体适定性（Theorem 3.1），作者在加权空间 $$H_2$$ 的正交基 $\{e_{2,k}\}$ 上构造有限维逼近解 $\eta_n^\pm = \sum_{k=1}^n \eta_{2,k}^\pm(t) e_{2,k}(\theta)$ ，将无限维 PDE 转化为常微分方程组。通过能量估计获得 $\eta_{2,k}^\pm$ 及其导数的一致先验界（见 (4.18)–(4.20)），再利用紧性（Ascoli–Arzelà）取极限，最终恢复原方程解。此方法的关键在于基函数 $e_{2,k}$ 的精心设计，使其能自动吸收权重 $\rho$ 引发的奇异性。
谱分析与二次型正定性论证：线性不稳定性（Theorem 3.3）的证明不依赖于传统特征值计算（因 $$L^+$$ 非自伴且定义域奇异），而是采用傅里叶系数的二阶 ODE 方法。将 $\eta^\pm$ 展开为 $\sum_k \eta_{2,k}^\pm(t) e_{2,k}$ 后，计算 $\frac{d}{dt}\|\eta^\pm\|_{H_2}^2 = \langle L^\pm \eta^\pm, \eta^\pm \rangle_{H_2}$ ，发现其可化为仅含平方项的级数：
$\frac{1}{2}\frac{d}{dt}\|\eta^+\|_{H_2}^2 = \sum_{k \ge 1} (d_{2,k}^+ - d_{2,k+2}^+) |\eta_{2,k}^+(t)|^2,$
其中系数 $d_{2,k}^+ - d_{2,k+2}^+ = -\frac{1}{2} - \frac{2k^2 - 4k - 8}{(k+2)^2 k^2}$ 。关键观察是：当 $$k=1$$ 时该系数为正（ $\frac{11}{18} > 0$ ），而 $k \ge 2$ 时为负。这暗示 $\eta_{2,1}^+$ 分量主导不稳定增长。作者进而推导出 $\|\eta^+\|_{H_2}^2$ 满足二阶 ODE：
$\frac{d^2}{dt^2}\|\eta^+\|_{H_2}^2 = \sum_{k \ge 1} (d_{2,k}^+ - d_{2,k+2}^+) \frac{d}{dt}|\eta_{2,k}^+(t)|^2,$
并通过对该 ODE 进行截断与正定二次型分解（见 (3.4)），证明存在 $\lambda_1,\lambda_2 > 0$ 使得
$4\lambda_1 \|\eta^+\|_{H_2}^2 < \frac{d^2}{dt^2}\|\eta^+\|_{H_2}^2 < 4\lambda_2 \|\eta^+\|_{H_2}^2.$
最后应用比较定理（Lemma A.6）得出指数下界。
连续性论证（continuity argument）与小初值能量估计：非线性稳定性（Theorem 3.10）的证明依赖于将非线性项视为扰动。作者首先证明线性部分对特定初值（仅含 $e_{2,2k}$ 分量）具有指数衰减，再通过精细估计非线性项 $$N_1,N_2$$ （见 (2.2)–(2.3)）在 $$H_2$$ 范数下的高阶小量性质（如 $\|N_1\|_{H_2} \lesssim \|\eta^+\|_{H_2}^2 + \|\eta^-\|_{H_2}^2$ ），最后结合连续性论证：假设解在时间区间 $$[0,T]$$ 内保持小能量，则可延拓至更大时间，从而推出全局存在性。

b) 处理技术性难题的技巧

奇异权重的处理：空间 $$H_2$$ 的权重 $\rho(\theta) = (4\pi \sin^2 \theta)^{-1}$ 在 $\theta = 0,\pm\pi$ 处发散，导致标准乘子估计失效。作者通过选用基函数 $e_{2,k}$ ，使 $\partial_\theta e_{2,k}$ 在奇点处具有可控的 $\sim |\theta|^{-1}$ 行为，从而保证 $\langle L^\pm \eta, \eta \rangle_{H_2}$ 的积分收敛，并能分离出主导不稳定模态。
非线性项的强奇异性控制：非线性项 $N_1 = \{\eta^+,v^+\} - \{\eta^-,v^-\}$ 包含 $\partial_\theta$ 导数与 Hilbert 变换的复合，其导数在 $\theta=0$ 附近有更强奇性。作者利用 $v^\pm = H^{-1} \partial_\theta \eta^\pm$ 的积分表示及 $\sin(2\theta)$ 的零点结构，将 $$N_1$$ 控制在 $$H_2$$ 中，关键引理（如 Lemma 5.3）表明：只要 $\eta^\pm$ 属于 $$H_2$$ ，则 $$N_1$$ 的加权 $$L^2$$ 范数可被 $\|\eta^\pm\|_{H_2}^2$ 控制。
Lipschitz 不稳定性框架的构建：为证明 Theorem 3.7，作者不直接估计非线性解，而是构造一个缩放序列 $\eta^\pm_\varepsilon(t,\theta) = \varepsilon \eta^\pm(t/\varepsilon,\theta)$ ，并证明当 $\varepsilon \to 0$ 时，该序列在适当拓扑下收敛到线性化解。这一“线性化极限”过程严格验证了非线性不稳定性由线性机制主导，是处理弱非线性扰动的经典策略。

与之前工作的比较

a) 突破与改进

本文是对 Guo–Jiu 2025 年关于 De Gregorio 模型（DG）激发态稳定性的开创性工作 [GJ25] 的实质性推广与深化，其突破体现在三个维度：

模型维度的提升：[GJ25] 仅研究单分量 DG 模型 $\partial_t \omega + u \partial_\theta \omega = \omega H \omega$ ，而本文首次将激发态稳定性分析拓展至双分量、强耦合的 MHD 模型。这不仅是形式上的推广，更揭示了磁场如何通过参数 $$q$$ 调制不稳定性阈值——例如，当 $$q > 0$$ 时， $$L^-$$ 算子的耗散增强（因 $$-2q Q$$ 项），但 $$L^+$$ 的不稳定性机制保持不变，从而实现了“磁场抑制 $\eta^-$ 不稳定但无法阻止 $\eta^+$ 不稳定”的新物理图景。
初值空间的优化与松弛：[GJ25] 要求初值在 $$H_2$$ 中且满足严格的符号条件（如 $\eta_0^+$ 与 $\sin\theta$ 同号），而本文的 Theorem 3.3 仅需 $\langle L^+ \eta_0^+, \eta_0^+ \rangle_{H_2} \ge 0$ ，这是一个纯谱条件，覆盖了更广泛的振荡初值（如 Remark 3.5 中的 $e_{2,1} + e_{2,k}$ 组合）。此外，Theorem 3.10 的稳定性结论针对的是子空间 $\mathrm{span}\{e_{2,2k}\}$ ，比 [GJ25] 的全局稳定性假设更弱，体现了对解结构的更精细刻画。
技术工具的创新性整合：[GJ25] 的核心是分析 $\frac{d^2}{dt^2}\|\eta\|_{H_2}^2$ 的符号，但未处理双分量耦合。本文不仅将此方法成功移植到耦合系统，还发展了跨分量能量耦合估计：在证明 Theorem 3.10 时，作者发现 $\eta^+$ 的偶数模态衰减可“拖拽” $\eta^-$ 一同衰减，这依赖于对非线性项 $$N_2$$ 中 $\eta^- H \eta^+$ 项的精细 cancellation 估计，是此前工作中未曾涉及的新技术难点。

综上，本文不仅将 DG 模型的理论推向更真实的物理模型（MHD），更在数学上建立了一套适用于奇异加权、非局部、多分量竞争系统的稳定性分析范式，为未来研究更高维激发态（ $\kappa \ge 3$ ）及粘性 MHD 模型奠定了坚实基础。

5. 四维空间中轴对称无旋欧拉方程解的整体正则性

Global regularity for axisymmetric, swirl-free solutions of the Euler equation in four dimensions

作者: Evan Miller

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.13963v1

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

详细分析

中文标题: 四维空间中轴对称无旋欧拉方程解的整体正则性

研究对象和背景

本文研究的是四维欧几里得空间 $\mathbb{R}^4$ 中的不可压理想流体运动，其数学模型为不可压欧拉方程（incompressible Euler equation）。该方程描述无粘性、不可压缩流体的速度场 $u:\mathbb{R}^4\times[0,T)\to\mathbb{R}^4$ 与压力 $p:\mathbb{R}^4\times[0,T)\to\mathbb{R}$ 的耦合演化，基本形式如下：

\partial_t u + (u\cdot\nabla)u + \nabla p = 0,

\nabla\cdot u = 0.

其中， $\partial_t$ 表示时间偏导， $(u\cdot\nabla)u$ 是非线性对流项， $\nabla p$ 为压力梯度， $\nabla\cdot u = 0$ 是不可压条件（即质量守恒）。该系统定义在全空间 $\mathbb{R}^4$ 上，无边界、无初边值约束，仅要求初始速度场 $u^0 = u(\cdot,0)$ 具有足够正则性与衰减性。

论文聚焦于一类具有高度对称性的解：轴对称、无旋（swirl-free）解。所谓“轴对称”，指解关于 $\mathbb{R}^4$ 中某固定轴（取为 $$x_4$$ -轴）旋转不变；“无旋”指速度场在环绕该轴方向（即“swirl”分量）为零。在四维情形下，采用广义柱坐标 $$(r,z)$$ ，其中

r = \sqrt{x_1^2 + x_2^2 + x_3^2},\quad z = x_4,

并定义单位向量场：

e_r = \frac{1}{r}(x_1,x_2,x_3,0),\quad e_z = (0,0,0,1).

则轴对称无旋速度场形如：

u(x) = u_r(r,z)\,e_r + u_z(r,z)\,e_z.

这一结构自动满足对称性保持性：若初值 $$u^0$$ 具有该形式，则欧拉方程的解在存在时间内亦保持此形式。值得注意的是，该对称性将原四维向量场问题降维为二维 $$(r,z)$$ 平面上的标量函数系统，极大简化分析，同时仍保留核心非线性结构（特别是涡量拉伸效应）。

在该对称框架下，散度自由条件化为：

\partial_r u_r + \partial_z u_z + \frac{2}{r}u_r = 0,

其中系数 $$2 = d-2$$ 来自四维空间中球面测度的 Jacobian $r^{d-2} = r^2$ 。相应地，标量涡量（scalar vorticity） 定义为：

\omega(r,z) = \partial_r u_z - \partial_z u_r,

它刻画了速度梯度的反对称部分，并满足如下关键演化方程：

\partial_t \omega + (u\cdot\nabla)\omega - 2\frac{u_r}{r}\omega = 0.

此处 $2\frac{u_r}{r}\omega$ 即为涡量拉伸项（vortex stretching term），是三维及以上欧拉方程中导致奇点可能形成的核心机制。当 $$d=4$$ 时，该系数恰为 $$d-2 = 2$$ ，凸显四维的临界性。

为消除拉伸项，引入加权涡量 $\frac{\omega}{r^2}$ ，可得纯输运方程：

\left(\partial_t + u\cdot\nabla\right)\left(\frac{\omega}{r^2}\right) = 0.

这表明 $\frac{\omega}{r^2}$ 沿流线守恒，因此其任意范数（如 $L^{p,q}$ ）在时间上不变——这是后续全局控制的关键。

函数空间设定极为关键：初值 $$u^0$$ 属于Sobolev 空间 $H^s_{\text{df}}(\mathbb{R}^4)$ （ $$s>4$$ ），即具有 $$s$$ 阶导数平方可积且散度为零的向量场；对应涡量 $\omega^0$ 则自然属于 $$L^2$$ 型空间。但传统方法（如 $L^\infty$ 控制）在此失效：即使 $$u^0$$ 是 Schwartz 函数， $\frac{\omega^0}{r^2}$ 也可能无界（因 $$r=0$$ 处的奇性）。例如文中反例：

u(x) = e^{-r^2-z^2}\big(r(1-2z^2)e_r - z(3-2r^2)e_z\big),

计算得 $\omega(r,z) = 4rz(4-r^2-z^2)e^{-r^2-z^2}$ ，故 $\frac{\omega}{r^2} \sim \frac{z}{r}$ 当 $r\to 0$ ，显然不属于 $L^\infty$ 。因此，必须寻求比 $L^\infty$ 更弱、却仍能控制 $$u_r/r$$ 的范数。

论文选择 Lorentz 空间 $L^{2,1}(\mathbb{R}^4)$ （以测度 $$r^2drdz$$ 计，等价于 $\mathbb{R}^4$ 上的标准 Lebesgue 测度），其准范数为：

\|f\|_{L^{2,1}} = \int_0^\infty \mu\big(\{|f|>\tau\}\big)^{1/2}\,d\tau.

该空间严格介于 $$L^2$$ 与 $L^{2,\infty}$ （弱 $$L^2$$ ）之间，具备良好插值性质与对偶性（ $(L^{2,1})^* = L^{2,\infty}$ ），且对初值正则性要求自然：只要 $u^0\in H^s$ , $$s>4$$ ，并具多项式衰减（如 $|\omega^0|\lesssim (1+r^2+z^2)^{-2}$ ），即可保证 $\frac{\omega^0}{r^2}\in L^{2,1}$ 。这使得结果具有广泛适用性，覆盖绝大多数物理意义明确的光滑初值。

主要定理或主要结果

论文的核心成果是以下全局正则性定理：

Theorem 1.1.
设初值 $u^0 \in H^s_{\text{df}}(\mathbb{R}^4)$ ， $$s > 4$$ ，且为轴对称、无旋；进一步假设加权涡量满足

\frac{\omega^0}{r^2} \in L^{2,1}(\mathbb{R}^4).

则欧拉方程存在唯一的全局光滑解

u \in C\big([0,+\infty); H^s(\mathbb{R}^4)\big),

且对任意 $$t>0$$ ，涡量满足指数型控制：

\|\omega(\cdot,t)\|_{L^\infty} \leq \left\|\frac{\omega^0}{r^2}\right\|_{L^\infty} \exp\!\left(C\left\|\frac{\omega^0}{r^2}\right\|_{L^{2,1}} t\right),

其中 $$C>0$$ 为绝对常数。此外，该假设对所有满足 $|\omega^0(r,z)| \lesssim (1+r^2+z^2)^{-2}$ 的初值均成立。

该定理的证明依赖两个支柱性引理：

涡量拉伸控制引理（Theorem 2.5）：存在常数 $$C>0$$ ，使得
$\left\|\frac{u_r}{r}\right\|_{L^\infty(\mathbb{R}^4)} \leq C \left\|\frac{\omega}{r^2}\right\|_{L^{2,1}(\mathbb{R}^4)}.$
加权涡量守恒命题（Proposition 3.1）：
$\left\|\frac{\omega(\cdot,t)}{r^2}\right\|_{L^{2,1}} = \left\|\frac{\omega^0}{r^2}\right\|_{L^{2,1}},\quad \forall\, t\geq 0.$

二者结合，立即给出 $\|\omega(\cdot,t)\|_{L^\infty} \leq \|\omega^0/r^2\|_{L^\infty} \exp(C\|\omega^0/r^2\|_{L^{2,1}} t)$ ，再由 Beale–Kato–Majda 准则（若解在有限时间爆破，则 $\int_0^{T_{\max}} \|\omega(\cdot,t)\|_{L^\infty}\,dt = \infty$ ）可知解永不失控，从而全局存在。

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

a) 核心分析学工具

论文综合运用了现代调和分析与偏微分方程理论中的若干关键工具：

Lorentz 空间理论：作为比 Lebesgue 空间更精细的尺度， $L^{p,q}$ 空间提供了对函数奇性与衰减的双重刻画。本文特别依赖 $L^{2,1}$ 与其对偶 $L^{2,\infty}$ 的广义 Hölder 不等式（Proposition 1.8）：
$\left|\int f g \,d\mu\right| \leq C \|f\|_{L^{2,1}} \|g\|_{L^{2,\infty}}.$

此不等式将速度分量 $$u_r/r$$ 的 $L^\infty$ 估计转化为对核函数 $g_{r,z}$ 的 $L^{2,\infty}$ 范数与加权涡量 $\omega/r^2$ 的 $L^{2,1}$ 范数的乘积控制，是连接几何结构与分析估计的桥梁。
积分核估计与尺度不变性分析：速度通过 Biot–Savart 律由涡量重构。在四维轴对称无旋情形下，径向速度 $$u_r$$ 的积分表示含奇异核 $g_{r,z}(r',z')$ （见 Proposition 2.3）。证明其 $L^{2,\infty}(r^2drdz)$ 一致有界（Lemma 2.4）是技术核心。作者利用缩放不变性（ $g_{a,b}(r,z) = a^{-2}g(r/a,z/a)$ ）将一般情况归约为单位核 $g_{1,0}$ 的估计，并通过分区域（ $\tau\geq 1$ 与 $\tau<1$ ）精细控制水平集测度 $\mu\{g>\tau\}$ ，最终得到 $\|g\|_{L^{2,\infty}} \leq 8\pi$ 。
输运方程守恒律：加权涡量 $\omega/r^2$ 满足纯输运方程 $\partial_t(\omega/r^2) + u\cdot\nabla(\omega/r^2)=0$ ，且速度场 $$u$$ 散度为零，故流映射 $X(\cdot,t)$ 是体积保持微分同胚。由此直接推出 $L^{2,1}$ 范数守恒（Proposition 3.1），无需额外能量估计。
Beale–Kato–Majda（BKM）准则：该经典爆破判据指出，强解 $u\in C([0,T);H^s)$ 在 $T<\infty$ 处爆破当且仅当 $\int_0^T \|\nabla u(\cdot,t)\|_{L^\infty}\,dt = \infty$ 。在轴对称无旋情形下， $\|\nabla u\|_{L^\infty} \lesssim \|\omega\|_{L^\infty}$ ，故只需控制 $\|\omega\|_{L^\infty}$ 。

b) 处理非线性与奇性的主要技巧

权重消除拉伸项：传统方法试图直接控制 $$u_r/r$$ ，但其表达式含 $$1/r$$ 奇性。本文洞察到 $\omega/r^2$ 的输运结构，将问题转化为对 $$u_r/r$$ 关于 $\omega/r^2$ 的单线性估计，绕过直接处理 $$1/r$$ 奇异性的困难。
核函数的分段精细估计：为证 $\|g_{r,z}\|_{L^{2,\infty}}$ 一致有界，作者将水平集 $\{g>\tau\}$ 分为两区域：当 $\tau\geq 1$ 时，利用分母中 $$(r-1)^2+z^2$$ 项主导，得到测度 $\lesssim \tau^{-2}$ ；当 $\tau<1$ 时，利用 $$1+r^2+z^2$$ 项主导，结合球坐标积分，仍得 $\mu\{g>\tau\}\lesssim \tau^{-2}$ 。这种按尺度划分、匹配主导项的策略是处理各向异性奇异积分的标准范式。
从局部适定性到全局延拓：借助 Kato 的局部适定性定理（Theorem 1.11），解在 $$H^s$$ 中短期存在；BKM 准则与前述指数控制共同保证 $$H^s$$ 范数在任意有限时间区间上有界，从而可无限延拓。

与之前工作的比较

a) 突破与改进

本文的根本性突破在于将全局正则性所需的初值条件，从先前工作强制要求的 $L^\infty$ 有界性，大幅削弱为 $L^{2,1}$ 可积性。

先前工作的局限：Lim–Jeong（2021）、Choi–Jeong–Lim（2019）及作者与 Tsai（2022）等均要求 $\frac{\omega^0}{r^2}\in L^\infty$ 。该条件极强：如前所述，即使 $$u^0$$ 是 Schwartz 函数， $\omega^0/r^2$ 仍可在 $$r=0$$ 处无界。这意味着这些结果不适用于 generic 初值，仅覆盖紧支撑或特殊对称性的子集。
本文的实质性放宽： $L^{2,1}$ 是比 $L^\infty$ 弱得多的条件。事实上，由 Sobolev 嵌入 $H^2(\mathbb{R}^4)\hookrightarrow L^{4}(\mathbb{R}^4)$ 及插值， $u^0\in H^s$ , $$s>4$$ 自动蕴含 $\omega^0\in H^{s-1}\subset L^p$ 对某些 $$p>2$$ ，进而 $\omega^0/r^2\in L^{2,1}$ 只需温和衰减（如 $(1+r^2+z^2)^{-2}$ ）。这使定理适用于几乎所有物理上合理的光滑初值。
维度特异性意义：在 $$d=3$$ 时，Yudovich 定理已用 $\omega^0/r\in L^\infty$ 解决全局正则性；在 $d\geq 5$ ，Elgindi 等数值与理论工作暗示爆破可能性增大。本文首次在 $$d=4$$ 这一临界维度确立 generic 全局正则性，填补了高维欧拉方程理论的关键空白。
方法论启示：本文示范了如何用 Lorentz 空间替代 $L^\infty$ 来驾驭涡量拉伸项，为其他高维流体力学问题（如 Navier–Stokes、Boussinesq）提供了新思路。其核估计技术亦可推广至更一般的各向异性算子。

6. 二维与三维欧拉方程的精确非定常解

Exact non-stationary solutions of the Euler equations in two and three dimensions

作者: Patrick Heslin, Stephen C. Preston

类别: math.AP, math.DG

发布日期: 2026-02-14

arXiv ID: 2602.13929v1

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

详细分析

中文标题: 二维与三维欧拉方程的精确非定常解

研究对象和背景

本文研究的是不可压缩理想流体在紧致黎曼流形 $$(M,g)$$ （可能带边界）上的运动，其数学模型为经典的欧拉方程组。该方程组首次由 Euler 于 1757 年提出，描述无粘性、不可压缩流体的速度场 $$u(t,x)$$ 与压力场 $$p(t,x)$$ 的耦合演化：

\partial_t u + \nabla_u u = -\nabla p,

\operatorname{div} u = 0,

u \parallel \partial M.

其中 $\nabla_u u$ 表示协变导数沿 $$u$$ 方向的作用， $\operatorname{div} u = 0$ 是不可压缩性约束，而 $u \parallel \partial M$ 表示速度场在边界上切向（即无穿透条件）。这些方程定义在时间区间 $t \in \mathbb{R}$ 与空间区域 $$M$$ 上，且要求解具有全局存在性、光滑性与非定常性——这正是本文的核心目标。

论文特别关注两类典型空间：二维曲面（如平面环面 $\mathbb{T}^2$ 、单位圆盘 $$D^2$$ 、标准球面 $$S^2$$ 、双曲圆盘 $$H^2$$ ）与三维流形（如标准三维球面 $$S^3$$ 、圆丛、环面丛等）。所有分析均在光滑函数空间中进行：记 $X_\mu(M)$ 为 $$M$$ 上所有光滑、散度为零、且在边界上切向的向量场构成的空间；进一步将其正交分解为：

X_\mu(M) = X_{\mu,\text{ex}}(M) \oplus H(M),

其中 $$H(M)$$ 是调和场子空间（即同时满足 $\operatorname{div} u = 0$ 与 $\operatorname{curl} u = 0$ 的场），而 $X_{\mu,\text{ex}}(M)$ 是其 $$L^2$$ 正交补空间，对应于“恰当”部分（exact part），即可以表示为某势函数的旋度或斜梯度。该分解源于 Hodge 分解理论，是整个构造的基础。

在二维情形下，每个 $v \in X_{\mu,\text{ex}}(M)$ 可唯一写成斜梯度形式：

v = \nabla^\perp \psi = (\star d\psi)^\sharp,

其中 $\psi$ 称为流函数（stream function）， $\star$ 是 Hodge 星算子， $\sharp$ 是音乐同构。此时惯性算子定义为 Hodge Laplacian：

A = \Delta = (d\delta + \delta d)(\cdot^\flat)^\sharp.

而在三维情形下，惯性算子取为 旋度算子：

A = \operatorname{curl} = (\star d(\cdot^\flat))^\sharp.

这一统一记号 $$A$$ 是 Arnold 几何框架的关键桥梁：它将动能内积 $\langle u, v \rangle_{L^2} = \int_M g(u,v)\, \mu$ 与右不变 $$L^2$$ 度量联系起来，并使欧拉方程重写为简洁的抽象形式：

\partial_t (Au) + [u, Au] = 0, \quad \operatorname{div} u = 0, \quad u \parallel \partial M.

此处 $[u, v] = \nabla_u v - \nabla_v u$ 是向量场的李括号。该形式凸显了欧拉方程作为体积保持微分同胚群 $\operatorname{Diff}_\mu(M)$ 上测地线方程的本质——这是 Arnold (1966) 开创性的几何观点，后由 Ebin–Marsden 在 Sobolev 框架中严格实现。

值得注意的是，本文不依赖于小初值、线性化或对称简化假设；相反，它致力于构造显式、全局、光滑、真正非定常（即非稳态）的解。传统上，已知的显式解极为稀少：例如平面上的 Kelvin 波、球面上的 Rossby–Haurwitz 波、旋转剪切流等，但它们或是定常的、或是仅存在于特殊对称背景下、或是依赖奇异涡量分布。本文突破在于：在更广泛的几何背景下（包括弯曲空间与高维流形），系统性地构造出一大类新解，且其存在性由流形自身的对称结构（Killing 场及其谱性质）所决定。

此外，所有构造均要求解属于 $C^\infty$ 类，并在整个时间轴 $\mathbb{R}$ 上定义（global-in-time），而非仅局部存在。这种强正则性与全局性，使得所得解不仅具理论价值，也为后续稳定性分析、数值验证及物理建模提供了坚实基础。

主要定理或主要结果

本文最核心的成果是以下定理，它给出了构造显式非定常欧拉解的充分代数条件：

Theorem 2.2.
设 $$(M,g)$$ 是一个（可能带边界的）黎曼流形， $u_0 \in X_{\mu,\text{ex}}(M)$ 是欧拉方程的一个定常解，即满足 $$[u_0, A u_0] = 0$$ 。若存在复向量场 $z = v + i w \in \mathbb{C} \otimes X_{\mu,\text{ex}}(M)$ （其中 $v,w \in X_{\mu,\text{ex}}(M)$ ）以及实常数 $\lambda, \alpha, \zeta$ ，使得

K_{u_0} z = -i \lambda z, \quad A z = \alpha z, \quad [u_0, z] = i \zeta z,

其中 $K_{u_0} : v \mapsto A^{-1}[v, A u_0]$ 是与 $$u_0$$ 关联的余伴随算子（coadjoint operator），则曲线

U(t) = u_0 + e^{i(\lambda - \zeta)t} z

是欧拉方程的复值解；其实部与虚部

U_R(t) = u_0 + \cos\big((\lambda - \zeta)t\big) v - \sin\big((\lambda - \zeta)t\big) w,

U_I(t) = \sin\big((\lambda - \zeta)t\big) v + \cos\big((\lambda - \zeta)t\big) w

分别是欧拉方程本身及其在 $$U_R$$ 处线性化的实值解。进一步地， $$U_R$$ 是定常解当且仅当 $\lambda = \zeta$ 。

该定理的威力在于将偏微分方程的求解问题转化为有限维李代数上的同时本征值问题。其关键前提是存在一个 Killing 场 $$X$$ 满足 $$AX$$ 仍是 Killing 场。对此，作者给出如下结构性结论：

Lemma 2.7.
若 $X \in X_{\mu,\text{ex}}(M)$ 是 Killing 场，且 $$AX$$ 也是 Killing 场，则三个算子

$K_X : z \mapsto A^{-1}[z, A X]$ ,
$A : z \mapsto A z$ ,
$\operatorname{ad}_X : z \mapsto -[X, z]$

在 $\mathbb{C} \otimes X_{\mu,\text{ex}}(M)$ 中存在一个离散的同时本征基（discrete simultaneous eigenbasis）。

由此导出分类结果：

二维完全分类（Proposition 3.1 + Sections 3.2–3.5）：若二维流形 $$M$$ 存在 Killing 场 $$X$$ 使得 $\Delta X$ 仍为 Killing 场，则 $$M$$ 必为常曲率曲面（constant curvature surface），即局部等距于 $\mathbb{R}^2$ 、 $$S^2$$ 或 $$H^2$$ 。所有经典与新解（Kelvin 波、Rossby–Haurwitz 波、双曲 Kelvin 波）均在此框架下统一生成。
三维部分分类（Section 4）：若三维流形 $$M$$ 存在 Killing 场 $$X$$ 使得 $\operatorname{curl} X$ 仍为 Killing 场，则必有 $g(X, \operatorname{curl} X)$ 为常数（Lemma 4.1）；进一步，若 $\operatorname{curl} X = f X$ ，则 $$|X|_g$$ 与 $$f$$ 均为常数。典型例子包括 $$S^3$$ 、某些圆丛与环面丛。

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

a) 核心分析学工具

本文并非采用传统 PDE 分析中的能量估计、紧性论证或不动点原理，而是立足于无穷维李群几何与谱理论的交叉。其核心工具包括：

Arnold 的几何力学框架：将 $\operatorname{Diff}_\mu(M)$ 视为无限维李群，赋予右不变 $$L^2$$ 度量，欧拉方程即为其测地线方程。此视角将动力学问题转化为群作用与共轭类的研究。
Hodge–de Rham 理论：用于建立 $X_\mu(M)$ 的正交分解，定义惯性算子 $$A$$ ，并保证其自伴性与紧逆性（Lemma 2.1）。特别地，在二维中 $A = \Delta$ 具有离散谱 $\{\alpha_k\}$ ，对应本征函数构成 $$L^2$$ 标准正交基；在三维中 $A = \operatorname{curl}$ 虽非紧，但限制在 $X_{\mu,\text{ex}}$ 上仍具有良好谱性质。
同时对角化技术（simultaneous diagonalization）：关键引理 2.7 的证明依赖于三个算子两两交换且均为（斜）自伴的事实。具体而言：
- 由 $$[X, AX] = 0$$ 得 $[\operatorname{ad}_X, A] = 0$ ；
- 由 $K_X = A^{-1} \operatorname{ad}_{AX}$ 及 $$[X, AX] = 0$$ 得 $[K_X, \operatorname{ad}_X] = 0$ ；
- Corollary 2.6 给出 $$K_X$$ 是 $$L^2$$ 斜自伴的， $\operatorname{ad}_X$ 在每个 $$A$$ -本征子空间 $$E_k$$ 上亦斜自伴（通过分部积分与边界条件验证）。

从而可在每个有限维 $$E_k$$ 上施行标准线性代数——寻找共同本征向量，最终拼成全局同时本征基。

Bochner 技术与 Weitzenböck 公式：用于二维分类（Proposition 3.1）。Weitzenböck 公式 $\Delta(v^\flat) = \nabla^*\nabla v^\flat + \operatorname{Ric}(v^\flat)$ 结合 Killing 条件 $\nabla^2 X = -R(\cdot,X)\cdot$ ，推出 $\Delta X = 2 \operatorname{Ric}(X)$ 。再结合 $\Delta X$ 仍为 Killing 的假设，强制曲率张量作用后保持对称性，最终导出 $\kappa$ 为常数。

b) 技术性难题的处理技巧

非线性项的精确封闭（exact closure of nonlinearity）：欧拉方程的非线性项 $$[u,Au]$$ 通常导致无限维耦合。本文巧妙选取扰动方向 $$z$$ 为 $$A$$ 、 $K_{u_0}$ 与 $\operatorname{ad}_{u_0}$ 的共同本征场，使得非线性项在展开中自动落回同一一维复子空间：计算显示 $[u_0, z] \propto z$ ， $[z, A u_0] \propto z$ ， $$[z, A z] = 0$$ ，从而整个方程在 $\operatorname{span}_\mathbb{C}\{z\}$ 上封闭，无需截断或近似。
边界条件的自然兼容：所有构造均基于 $X_{\mu,\text{ex}}(M)$ 中的场，其定义已隐含切向边界条件；而 Killing 场自动满足 $g(\nabla_v X, \nu) = 0$ （ $\nu$ 为单位外法向），保证了 $\Delta X$ 或 $\operatorname{curl} X$ 仍切向，从而维持解的空间适定性。
复化解到实解的提取机制：利用复指数振荡因子 $e^{i(\lambda-\zeta)t}$ 与实部/虚部分离，既保留非定常性（只要 $\lambda \neq \zeta$ ），又确保实解满足原始实方程。该技巧避免了复化带来的物理意义模糊，且与线性化分析无缝衔接。
广义科里奥利力解释：将 $K_{u_0}$ 视为由背景流 $$u_0$$ 诱导的“广义科里奥利力算子”，其谱 $\lambda$ 决定了扰动的旋转频率。这一物理图像将 Rossby–Haurwitz 波等经典现象纳入统一几何框架，超越了传统旋转坐标系的局限。

与之前工作的比较

a) 突破与改进

本文相较以往工作实现了三重根本性突破：

第一，维度与几何背景的根本拓展。
此前绝大多数显式解局限于 $\mathbb{R}^2$ 或 $$S^2$$ 上的高度对称情形（如轴对称、平移不变），且难以推广至三维或负曲率空间。本文首次系统构建了：

在双曲空间 $$H^2$$ 上的 Kelvin 波（Example 3.5），其频移项含双曲特征值 $\alpha_{n,m} = -(1/4 + \beta_{n,m}^2)$ ，揭示负曲率导致的反向色散；
在三维球面 $$S^3$$ 上的高维 Rossby–Haurwitz 类解，具有更高多重性与丰富的时间周期结构；
在非平凡纤维丛（如 circle bundles over surfaces）上的解，体现拓扑非平凡性对流体动力学的深刻影响。

此类构造此前未见报道，标志着从“平坦/球面对称”到“任意常曲率/纤维化几何”的范式跃迁。

第二，正则性与全局性的双重保障。
既有显式解常依赖弱解、分布涡量或无穷能量（如点涡、偶极子），而本文所有解均为 $C^\infty$ 光滑、有限动能、全局定义。这得益于：

基于紧致流形上的离散谱理论，避免了非紧区域的连续谱困难；
利用 Killing 场的刚性约束，规避了奇异性产生机制（如涡拉伸爆炸）；
不依赖小参数展开，故解的有效性不随振幅衰减。

第三，判据的普适性与可计算性。
以往判别非定常性的方法多依赖数值模拟或特定 ansatz，而本文给出纯代数判据： $\lambda \neq \zeta$ 即保证非定常。该判据可通过计算 Killing 场的 Lie 导数与 Laplacian/curl 的交换关系直接验证，已在多个例子中显式实现（如 $$S^2$$ 上 $\lambda = 2n/(m(m+1)) \neq n = \zeta$ 当 $$m > 1$$ ）。这为自动化搜索新解提供了算法基础。

综上，本文不仅大幅扩充了已知精确解的版图，更建立了连接微分几何、李群表示、流体力学与谱理论的全新方法论，为未来探索湍流前兆、模式选择与几何稳定性开辟了道路。

7. 黏性-色散激波的收缩：零黏度-毛细作用极限

Contraction of viscous-dispersive shocks: Zero viscosity-capillarity limits

作者: Namhyun Eun, Moon-Jin Kang, Jeongho Kim

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13788v1

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

详细分析

中文标题: 黏性-色散激波的收缩：零黏度-毛细作用极限

研究对象和背景

本文研究的核心对象是一维Navier–Stokes–Korteweg (NSK) 系统，该系统描述了具有黏性（viscosity）与毛细作用（capillarity） 的可压缩流体动力学行为，尤其适用于量子流体、两相界面扩散模型及超流体等物理场景。系统在拉格朗日质量坐标 $$(t,x)$$ 下写为：

\begin{cases} v_t - u_x = 0, \\ u_t + p(v)_x = \big( \mu(v) u_x \big)_x + \left[ \kappa(v) \left( -\dfrac{v_{xx}}{v} + \dfrac{5(v_x)^2}{2v^6} \right) - \dfrac{\kappa'(v)(v_x)^2}{2v^5} \right]_x. \end{cases} \tag{1.1}

其中：

$$v(t,x) > 0$$ 是比容（specific volume），即密度 $\rho = v^{-1}$ 的倒数；
$u(t,x) \in \mathbb{R}$ 是流体速度；
$p(v) = v^{-\gamma}$ 是压强函数，对应多方气体（polytropic gas）， $\gamma \ge 1$ ；
$\mu(v) = b v^{-\alpha}$ 是黏性系数，允许在真空附近（ $v \to \infty$ ）退化（degenerate）， $\alpha \ge 0$ ；
$\kappa(v) = c \mu(v)^2 v^3 = c b^2 v^{3 - 2\alpha}$ 是毛细系数，由物理一致性要求与黏性耦合， $$c > 0$$ 。

该系统源于 Dunn–Serrin 的热力学推导，并可由非线性薛定谔方程经 Madelung 变换导出，因此在量子流体力学中具有根本意义；例如当 $\kappa(v) = v$ （即 $\alpha = 1$ , $c = b^{-2}$ ）时，它退化为 Euler–Korteweg 系统，对应于 Bose–Einstein 凝聚体或超流 He-II 的宏观描述。

研究区域为全空间 $\mathbb{R}$ ，时间 $t \ge 0$ ，初值满足大初值扰动（large initial perturbation） 条件，但不作小参数假设。解所处的函数空间为：

X_T := \Big\{ (v,u) \,\Big|\, v - \tilde{v} \in C([0,T]; H^2(\mathbb{R})) \cap L^2(0,T; H^3(\mathbb{R})),\; 0 < v^{-1} \in L^\infty((0,T)\times\mathbb{R}), \\ u - \tilde{u} \in C([0,T]; H^1(\mathbb{R})) \cap L^2(0,T; H^2(\mathbb{R})) \Big\},

其中 $(\tilde{v}, \tilde{u})$ 是连接两个常状态 $(v_\pm, u_\pm)$ 的黏性-色散激波（viscous-dispersive shock），即行波解 $(v,u)(t,x) = (\tilde{v}, \tilde{u})(x - \sigma t)$ ，满足 ODE 系统：

\begin{cases} -\sigma \tilde{v}' - \tilde{u}' = 0, \\ -\sigma \tilde{u}' + p(\tilde{v})' = \big( \mu(\tilde{v}) \tilde{u}' \big)' + \left[ \kappa(\tilde{v}) \left( -\dfrac{\tilde{v}''}{\tilde{v}} + \dfrac{5(\tilde{v}')^2}{2\tilde{v}^6} \right) - \dfrac{\kappa'(\tilde{v})(\tilde{v}')^2}{2\tilde{v}^5} \right]'. \end{cases} \tag{1.4}

边界条件为 $\lim_{\xi \to \pm\infty} (\tilde{v}, \tilde{u})(\xi) = (v_\pm, u_\pm)$ ，且 $(v_\pm, u_\pm)$ 满足 Rankine–Hugoniot 条件 与 Lax 熵条件：

\begin{cases} -\sigma (v_+ - v_-) - (u_+ - u_-) = 0, \\ -\sigma (u_+ - u_-) + p(v_+) - p(v_-) = 0, \end{cases} \quad\text{且}\quad \text{either } v_- > v_+,\, u_- > u_+ \text{ or } v_- < v_+,\, u_- > u_+.

该激波是非单调、振荡型（dispersive） 的，区别于纯黏性 Navier–Stokes 激波——后者无高频振荡结构；而 NSK 激波因毛细项引入二阶色散效应，导致剖面呈现“震荡尾迹（oscillatory tail）”，是激波结构理论中黏性与色散竞争平衡的典型范例。

特别地，作者考虑零黏度-毛细极限（zero viscosity-capillarity limit），即对缩放参数 $\nu > 0$ ，考察如下尺度化系统：

\begin{cases} v^\nu_t - u^\nu_x = 0, \\ u^\nu_t + p(v^\nu)_x = \nu \big( \mu(v^\nu) u^\nu_x \big)_x + \nu^2 \left[ \kappa(v^\nu) \left( -\dfrac{v^\nu_{xx}}{v^\nu} + \dfrac{5(v^\nu_x)^2}{2(v^\nu)^6} \right) - \dfrac{\kappa'(v^\nu)(v^\nu_x)^2}{2(v^\nu)^5} \right]_x. \end{cases} \tag{1.6}

目标是证明：当 $\nu \to 0$ 时，解 $(v^\nu, u^\nu)$ 弱收敛至Riemann 熵解：

(\bar{v}, \bar{u})(x - \sigma t) = \begin{cases} (v_-, u_-), & x - \sigma t < 0, \\ (v_+, u_+), & x - \sigma t > 0, \end{cases} \tag{1.7}

该解是等熵 Euler 系统

\begin{cases} v_t - u_x = 0, \\ u_t + p(v)_x = 0, \end{cases} \tag{1.8}

的激波型弱解。本文首次在无小初值、无 BV 正则性假设、仅依赖相对熵控制的前提下，建立了该极限的唯一性与稳定性（up to shifts），填补了该方向长期存在的空白。

主要定理或主要结果

论文确立了三个相互支撑的核心定理，构成完整逻辑链：

定理 1.1（收缩性定理）

设 $\gamma, \alpha$ 满足

1 \le \gamma \le \frac{5}{4}, \quad 0 \le \alpha \le \gamma \le 1 + \frac{\alpha}{3},

且 $c \in [0, 9/100]$ （即 $\tau_1 = \frac{1 + \sqrt{1 - 4c}}{2} \ge 9/10$ ）。给定左常状态 $(v_-, u_-) \in \mathbb{R}_+ \times \mathbb{R}$ ，存在 $\varepsilon_0, \delta_0 > 0$ ，使得对任意 $\varepsilon \in (0,\varepsilon_0)$ 、 $\lambda \in (\delta_0^{-1}\varepsilon, \delta_0)$ 、右状态 $$(v_+, u_+)$$ 满足 $|p(v_+) - p(v_-)| = \varepsilon$ 与 Rankine–Hugoniot 条件，令 $(\tilde{v}, \tilde{u})$ 为对应黏性-色散激波， $(v,u) \in X_T$ 为其大初值扰动解，且初始相对熵有限：

\int_\mathbb{R} E\big((v_0,u_0) \,\big|\, (\tilde{v}, \tilde{u})\big)\,dx < \infty,

其中

E\big((v,u) \,\big|\, (\tilde{v}, \tilde{u})\big) := \frac{1}{2}\left(u - \frac{\tau_1 \gamma}{v^{\alpha+1}} v_x - \tilde{u} + \frac{\tau_1 \gamma}{\tilde{v}^{\alpha+1}} \tilde{v}' \right)^2 + Q(v \,|\, \tilde{v}),

$$Q(v|w) = Q(v) - Q(w) - Q'(w)(v-w)$$ ， $Q(v) = \begin{cases} \frac{v^{-\gamma+1}}{\gamma-1}, & \gamma > 1, \\ -\log v, & \gamma = 1. \end{cases}$

则存在：

光滑单调权函数 $a:\mathbb{R} \to \mathbb{R}_+$ ，满足 $\lim_{x \to \pm\infty} a(x) = 1 + a_\pm$ ， $|a_+ - a_-| = \lambda$ ；
动态位移 $X \in W^{1,1}(0,T)$ ， $$X(0)=0$$ ；
常数 $$C > 0$$ （独立于 $\varepsilon,\lambda,\delta_0,T$ ），

使得对所有 $t \in [0,T]$ ，

\int_\mathbb{R} a(\xi) E\big((v,u)(t,\xi+X(t)) \,\big|\, (\tilde{v}, \tilde{u})(\xi)\big)\,d\xi + \text{dissipative terms} \le \int_\mathbb{R} a(\xi) E\big((v_0,u_0)(\xi) \,\big|\, (\tilde{v}, \tilde{u})(\xi)\big)\,d\xi.

该不等式中耗散项包含加权 $$L^2$$ 范数控制 $Q(v|\tilde{v})$ 、 $p(v)-p(\tilde{v})$ 与有效速度差 $h - \tilde{h}$ 的梯度，详见 (1.12)。

定理 1.2（全局存在性）

在相同参数约束下，对任意满足 $v_0 - \tilde{v} \in H^2$ , $u_0 - \tilde{u} \in H^1$ , $v_0^{-1} \in L^\infty$ 的初值，NSK 系统存在唯一全局强解 $(v,u) \in X_T$ 对所有 $$T>0$$ 成立。

定理 1.3（零极限的唯一性与稳定性）

设 $\{(v^\nu, u^\nu)\}_{\nu > 0}$ 是 (1.6) 的光滑初值解序列，初值按 (1.15) “良好准备（well-prepared）”，则存在子列弱-* 收敛至 Radon 测度 $(v_\infty, u_\infty) \in L^\infty(0,T; L^\infty + \mathcal{M})$ ，且存在位移 $X_\infty \in \mathrm{BV}(0,T)$ ，使得对 a.e. $$t$$ ，

\int_\mathbb{R} dQ(v_\infty | \bar{v}(\cdot - X_\infty(t))) + \frac{1}{2} \int_\mathbb{R} |u_\infty(t,x) - \bar{u}(x - X_\infty(t))|^2 \,dx \le C E_0.

进一步， $|X_\infty(t) - \sigma t| \le C \sqrt{E_0}/|v_+ - v_-|$ ，表明 Riemann 激波在零黏度-毛细极限下唯一且稳定（模平移）。

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

a) 核心分析工具

相对熵方法（Relative Entropy Method）
这是处理守恒律大初值稳定性问题的基石工具。本文采用的相对熵基于 Bresch–Desjardins 型熵，定义为：
$\eta\big((v,h) \,\big|\, (\tilde{v}, \tilde{h})\big) := \frac{1}{2}(h - \tilde{h})^2 + Q(v \,|\, \tilde{v}),$
其中 $h = u - \tau_1 \mu(v) v_x / v = u - \tau_1 \gamma v^{-\alpha-1} v_x$ 是有效速度（effective velocity）。该构造巧妙地将 NSK 系统转化为两个退化抛物型方程：
$\begin{cases} v_t - h_x = -\tau_1 \big( v^{\gamma - \alpha} p(v)_x \big)_x, \\ h_t + p(v)_x = \tau_2 \gamma \big( v^{-\alpha - 1} h_x \big)_x, \end{cases} \tag{2.4}$
其中 $\tau_2 = 1 - \tau_1$ 。该变换消除了原始系统中高阶非线性耦合项（如 $\kappa(v) v_{xxx}$ ），使能量估计可行。
$$a$$ -收缩法（ $$a$$ -contraction with shifts）
该方法由 Vasseur 等人发展，用于处理激波稳定性中非均匀权重与动态位移。其核心是构造时空依赖的权函数 $$a(x)$$ 和位移 $$X(t)$$ ，使得泛函
$\mathcal{F}(t) := \int_\mathbb{R} a(x) \, \eta\big((v,h)(t,x+X(t)) \,\big|\, (\tilde{v}, \tilde{h})(x)\big)\,dx$
满足 $\frac{d}{dt}\mathcal{F}(t) \le 0$ 。关键在于选择 $$a(x)$$ 满足 $a' \sim \mathrm{sgn}(x)$ ，以补偿激波穿越时熵通量的跳跃；位移 $$X(t)$$ 则通过 ODE 控制： $\dot{X}(t) \sim \int a' \cdot \text{flux difference}$ 。
紧性原理与弱收敛分析
为证零极限，作者利用收缩性导出一致估计： $\{v^\nu\}$ 在 $L^\infty_t(L^1_x)$ 中有界， $\{h^\nu\}$ 在 $$L^2_t(H^1_x)$$ 中有界。结合 Banach–Alaoglu 定理，得弱-* 收敛 $v^\nu \overset{*}{\rightharpoonup} v_\infty \in L^\infty_t(\mathcal{M}_x)$ 。广义相对熵 (1.14) 的下半连续性保证极限满足稳定性不等式。

b) 技术难点与处理技巧

退化黏性与毛细耦合的处理：因 $\mu(v) = b v^{-\alpha}$ 在 $v \to \infty$ （真空）退化，标准 Sobolev 嵌入失效。作者利用 $v^{-1} \in L^\infty$ 的先验控制，将所有估计建立在 $$v$$ 的正下界上，并借助 $$Q(v|w)$$ 的局部二次性（Lemma 3.2）绕过 $$v$$ 的 $L^\infty$ 控制。
色散项的振荡破坏能量凸性：毛细项 $\kappa(v) v_{xxx}$ 不满足最大值原理，传统能量法失效。通过有效速度 $$h$$ 变换，将色散效应“吸收”进新黏性项 $\tau_2 \gamma (v^{-\alpha-1} h_x)_x$ ，从而恢复抛物结构。
相对熵的奇异部分处理：极限中可能出现空化（cavitation），即 $v_\infty$ 含 Dirac 质点。作者用 Radon–Nikodym 分解 $$dv = v_a dx + dv_s$$ ，并定义广义相对熵测度 $dQ(v|\bar{v}) = Q(v_a|\bar{v})dx + |Q'(V)| dv_s$ ，其中 $$V$$ 取极值以控制奇异部增长。
位移函数的构造：不同于固定位移， $$X(t)$$ 需满足 $\dot{X}(t) \lesssim \varepsilon^{-2} (f(t)+1)$ ，其中 $f \in L^1$ 由初始熵控制。该设计确保 $$X$$ 有界变差（BV），从而可提取收敛子列。

与之前工作的比较

a) 突破与改进

首次实现“大初值 + 零黏度-毛细极限”的唯一性
此前工作（如 Charve–Haspot [10], Germain–LeFloch [28]）仅证明解收敛至某个熵解，但无法排除多个不同熵解的可能性；而本文通过收缩性直接推出极限必为 Riemann 激波 $\bar{v}$ （模平移），解决了唯一性鸿沟。
摆脱 BV 或小初值假设
经典激波稳定性（如 Glimm 方案、Bressan 的 BV 理论）要求初值具小总变差；而本文仅需 $$L^2$$ -型相对熵有限，适用于任意大的 $H^2 \times H^1$ 扰动，极大拓展适用范围。
统一处理黏性与色散的标度
以往文献（如 [2, 49]）多假设 $\mu \sim \nu$ , $\kappa \sim \nu^2$ ，但未证明该标度下极限的稳定性；本文不仅确认该标度最优，更证明收缩性不依赖 $\nu$ ，从而自然导出极限稳定性。
物理模型的严格验证
本文为量子流体中“经典激波由微观色散-黏性机制选择”的猜想提供了首个严格数学证据：NSK 解族在 $\nu \to 0$ 下唯一收敛至 Lax 激波，证实该激波是物理上可实现的、稳定的宏观结构。

综上，该工作将相对熵方法、 $$a$$ -收缩法与退化抛物理论深度融合，为高阶耗散-色散系统的渐近分析树立了新范式。

8. Shinbrot型判据与可压缩Navier-Stokes方程的能量守恒

Shinbrot Type Criteria for Energy Conservation of the Compressible Navier-Stokes Equations

作者: Ruxuan Chen, Qi Zhang, Zhikang Zhang, Xiongbo Zheng

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13696v1

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

详细分析

中文标题: Shinbrot型判据与可压缩Navier-Stokes方程的能量守恒

研究对象和背景

本文研究的是可压缩Navier–Stokes方程组（Compressible Navier–Stokes equations）在有界区域 $\Omega \subset \mathbb{R}^N$ （ $$N=2,3$$ ）或环面 $\mathbb{T}^N$ 上的弱解的能量守恒问题。该方程组描述了粘性、可压缩流体的运动，是流体力学中最基本且最具挑战性的偏微分方程模型之一。

论文在第2页明确给出了所研究的控制方程（即系统(1)）：

\partial_t(\rho u) + \mathrm{div}(\rho u \otimes u) - \mu \Delta u - (\mu + \lambda) \nabla \mathrm{div}\, u + \nabla P = 0,

\partial_t \rho + \mathrm{div}(\rho u) = 0,

其中：

$\rho = \rho(t,x) \geq 0$ 是密度函数；
$u = u(t,x) \in \mathbb{R}^N$ 是速度场；
$P = \rho^\gamma$ 是压力项， $\gamma > 1$ 为绝热指数；
$\mu > 0$ , $2\mu + N\lambda \geq 0$ 是常数粘性系数（满足物理意义的椭圆性条件）；
初始数据为 $\rho|_{t=0} = \rho_0(x),\; (\rho u)|_{t=0} = \rho_0 u_0(x)$ ，且定义 $$u_0 = 0$$ 在 $\{x \mid \rho_0(x) = 0\}$ 上（即真空处速度被自然归零）。

该系统定义在时空域 $(0,T) \times \Omega$ 上，边界条件依情形而定：

若 $\Omega$ 是具有 $$C^1$$ 边界的有界区域，则采用无滑移边界条件 $u|_{\partial \Omega} = 0$ ；
若 $\Omega = \mathbb{T}^N$ （即周期性情形），则无需边界处理，所有函数视为周期延拓。

论文强调，尽管该系统存在全局弱解（参见文献[5,6,21,26,31]），但其能量等式一般不成立——即经典能量守恒律：

\int_\Omega \left( \frac{1}{2} \rho |u|^2 + \frac{\rho^\gamma}{\gamma - 1} \right) dx \Big|_{t=0}^{t} = \int_0^t \int_\Omega \left( \mu |\nabla u|^2 + (\mu + \lambda) |\mathrm{div}\, u|^2 \right) dx dt

并非对任意弱解自动满足，而仅在附加正则性条件下才成立。这种现象被称为能量耗散异常（energy dissipation anomaly），与Onsager猜想密切相关：在湍流中，即使粘性趋于零（Euler极限），若解正则性不足，仍可能产生非零能量耗散。

在不可压情形中，Lions（1969）、Serrin（1962）、Shinbrot（1974）等人建立了经典的能量守恒判据。例如，Shinbrot证明：若速度 $u \in L^p_t L^q_x$ 满足

\frac{2}{p} + \frac{2}{q} \leq 1,\quad q \geq 4,

或（当 $3 \leq q < 4$ 时）

\frac{1}{p} + \frac{3}{q} \leq 1,

则不可压Navier–Stokes弱解满足能量等式。这些条件显著弱于Serrin型判据（如 $\frac{2}{p} + \frac{N}{q} \leq 1$ ），因其摆脱了维度依赖，更具普适性。

然而，在可压情形下，由于密度 $\rho$ 可变、真空可能出现、非线性耦合更强（如 $\rho u \otimes u$ 、 $\nabla P = \nabla(\rho^\gamma)$ ），能量守恒分析面临根本困难。特别是时间导数项 $\partial_t(\rho u)$ 的处理需时空联合磨光（spatio-temporal mollification），引发新的时间方向上的交换子误差（temporal commutator error）。此前工作（如Yu [33], Chen–Zhang–Zheng [9], Ye–Zhu [32]）均依赖对 $\partial_t \rho$ 的强可积性估计（例如要求 $$u$$ 满足Serrin型条件以推出 $\partial_t \rho \in L^p_t L^q_x$ ），从而使用经典Lions交换子估计。但此类假设强于Shinbrot型，形成“正则性鸿沟”。

此外，论文特别关注真空允许情形（即 $\rho$ 可为零），此时 $\sqrt{\rho}$ 的正则性成为关键——Bresch–Desjardins熵结构暗示 $\nabla \sqrt{\rho} \in L^\infty_t L^{3/2}_x$ 是自然且较弱的假设（见(10)），远弱于传统要求的 $\sqrt{\rho} \in L^\infty_t H^1_x$ 。

综上，本文的研究背景是：在更弱的Shinbrot型速度正则性假设下，克服可压系统特有的时空耦合与真空奇异性，建立能量守恒的充分条件，从而弥合可压与不可压、Shinbrot与Serrin两类判据之间的理论断层。

主要定理或主要结果

论文的核心结果是 Theorem 1.1，其完整陈述如下：

Theorem 1.1 设 $(\rho,u)$ 是方程组 (1)–(2) 在 Definition 2.2 意义下的弱解。假设密度满足
$0 \leq \rho(t,x) \leq \bar{\rho} < \infty,\quad \nabla \sqrt{\rho} \in L^\infty\big(0,T; L^{3/2}(\Omega)\big),$

初始速度满足
$u_0 \in L^{\frac{6N}{6-N}}(\Omega),\quad N=2,3,$

且速度 $$u$$ 满足 Shinbrot 型可积性条件：
$u \in L^p\big(0,T; L^q(\Omega)\big),\quad \text{其中 } \begin{cases} \dfrac{1}{p} + \dfrac{3}{q} \leq 1, & \text{若 } 3 \leq q < 4, \\ \dfrac{2}{p} + \dfrac{2}{q} \leq 1, & \text{若 } 4 \leq q < \infty. \end{cases}$

则对任意 $t \in [0,T]$ ，能量等式 (3) 成立：
$\int_\Omega \left( \frac{1}{2} \rho |u|^2 + \frac{\rho^\gamma}{\gamma - 1} \right) dx \Big|_{0}^{t} = \int_0^t \int_\Omega \left( \mu |\nabla u|^2 + (\mu + \lambda) |\mathrm{div}\, u|^2 \right) dx dt.$

该定理适用于常粘性（constant viscosity）情形，并在 Remark 1.4 中指出方法可推广至退化粘性（degenerate viscosity）情形（见 Section 5 的 Theorem 5.3）。值得注意的是，条件中未要求 $\partial_t \rho$ 具有任何 $$L^p_t L^q_x$$ 可积性，亦未要求 $\sqrt{\rho} \in H^1$ ，仅需 $\nabla \sqrt{\rho} \in L^\infty_t L^{3/2}_x$ ，这比 Yu [33] 和 Chen et al. [9] 所需的 $L^\infty_t H^1_x$ 显著更弱。

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

a) 核心分析学工具

本文构建了一个两阶段能量守恒框架：

局部能量等式（local energy equality）：在测试函数 $\varphi \in C_c^\infty((0,T)\times\Omega)$ 下，证明分布意义下的等式
$$ \frac{1}{2} \int_0^T \int_\Omega \partial_t \varphi , \rho |u|^2 , dx dt
- \int_0^T \int_\Omega \varphi \Big( \rho^\gamma \mathrm{div}, u + \mu |\nabla u|^2 + (\mu + \lambda) (\mathrm{div}, u)^2 \Big) dx dt
- \int_0^T \int_\Omega \nabla \varphi \cdot \Big( \frac{1}{2} \rho |u|^2 u + \rho^\gamma u + \mu u \cdot \nabla u + (\mu + \lambda) u , \mathrm{div}, u \Big) dx dt = 0. $$
  此步骤的关键在于选择特殊测试函数 $(\varphi u^\varepsilon)^\varepsilon$ （先空间磨光再时间磨光），从而分离时空误差。
全局化（globalization）：通过截断函数 $\varphi = \psi_\tau \phi_\delta$ （ $\psi_\tau$ 时间截断， $\phi_\delta$ 空间边界截断）并取极限 $\tau,\delta \to 0$ ，将局部等式提升至全域。

支撑该框架的三大工具为：

时空联合磨光算子：定义 $f^\varepsilon(t,x) = (\eta_{x,\varepsilon} \ast_x f)^\varepsilon_t = \eta_{t,\varepsilon} \ast_t (\eta_{x,\varepsilon} \ast_x f)$ ，区别于纯空间磨光，能同时控制 $\partial_t(\rho u)$ 与非线性对流项。
紧性引理（Aubin–Lions型）：Lemma 2.4 给出嵌入 $W^{s_0,r_0}_t X \cap W^{s_1,r_1}_t Y \hookrightarrow C_t B$ 的紧性，用于处理 $\rho u$ 在负索伯列夫空间中的收敛性。
分布值时间导数的对偶估计：利用 Definition 2.1 中的负索伯列夫空间 $W^{-1,q_1}(\Omega)$ 范数，将 $\partial_t \rho$ 视为有界线性泛函，避免对其直接赋范。

b) 关键技术：弱型交换子估计

本文最核心的创新是 Lemma 2.5 —— 一个弱型时间交换子估计（weak-type temporal commutator estimate）：

Lemma 2.5 设 $f \in L^{p_1}_t L^{q_1}_x$ , $\partial_t f \in L^{p_1}_t W^{-1,q_1}_x$ , $g \in L^{p_2}_t W^{1,q_2}_x$ , $\varphi \in L^{\bar{p}}_t W^{1,\bar{q}}_x$ ，且满足
$\frac{1}{\bar{p}} + \frac{1}{p_1} + \frac{1}{p_2} = 1,\quad \frac{1}{\bar{q}} + \frac{1}{q_1} + \frac{1}{q_2} = 1.$

则当 $p_2, q_2, \bar{p}, \bar{q} < \infty$ 时，有
$\int_0^T \int_\Omega \varphi \Big[ \partial_t(f g)^\varepsilon - \partial_t(f g^\varepsilon) \Big] dx dt \to 0,\quad \text{as } \varepsilon \to 0.$

该引理的精妙之处在于：

它不预设 $\partial_t f \in L^{p_1}_t L^{q_1}_x$ ，而仅需 $\partial_t f \in L^{p_1}_t W^{-1,q_1}_x$ ，即允许 $\partial_t f$ 仅为分布；
收敛性由三因子乘积范数控制： $\|\partial_t f\|_{L^{p_1}_t W^{-1,q_1}_x} \|g\|_{L^{p_2}_t W^{1,q_2}_x} \|\varphi\|_{L^{\bar{p}}_t W^{1,\bar{q}}_x}$ ；
证明中采用时间平均技巧：将 $\partial_t(f g)^\varepsilon - \partial_t(f g^\varepsilon)$ 表示为 $\int_{t-\varepsilon}^{t+\varepsilon} \partial_t \eta_t\big(\frac{t-s}{\varepsilon}\big) \langle \partial_t f(\tau), g(s) \varphi(t) \rangle d\tau ds dt$ ，再利用对偶配对与Hölder不等式完成估计。

此引理直接应用于主项 $I_1 = \int \varphi u^\varepsilon \big[ \partial_t(\rho u)^\varepsilon - \partial_t(\rho u^\varepsilon) \big]$ ，其中 $f = \rho u$ , $$g = u$$ , $\varphi = \varphi u^\varepsilon$ 。由质量方程 $\partial_t \rho = -\mathrm{div}(\rho u)$ ，立即得 $\partial_t(\rho u) \in L^p_t W^{-1,q}_x$ （因 $\rho u \in L^p_t L^q_x$ ），从而绕过对 $\partial_t \rho$ 的强估计，实现Shinbrot型条件下的误差消没。

此外，论文还运用了：

Bresch–Desjardins 熵估计：隐含 $\nabla \sqrt{\rho} \in L^\infty_t L^{3/2}_x$ ，用于控制压力项 $\nabla P = \gamma \rho^{\gamma-1} \nabla \rho = 2\gamma \rho^{\gamma-1/2} \nabla \sqrt{\rho}$ ；
边界截断技术（Lemma 2.3）：利用 $\mathrm{dist}(x,\partial\Omega)$ 加权范数控制近壁区解的爆破行为；
插值不等式：如 $L^p_t L^q_x \cap L^2_t H^1_x \hookrightarrow L^4_t L^4_x$ ，用于处理临界情形。

与之前工作的比较

a) 突破性改进

本文宣称的突破集中于三点：

第一，正则性假设的实质性放宽。
先前工作均要求速度满足Serrin型条件或其端点：

Yu [33] 要求 $u \in L^p_t L^q_x$ 满足 $\frac{1}{p} + \frac{5}{q} \leq \frac{5}{12},\, q \geq 6$ （即 $$(p,q) = (12/5,6)$$ ）；
Chen et al. [9] 要求 $p \geq 4,\, q \geq 6$ （即 $$(p,q) = (4,6)$$ ）；
Ye–Zhu [32] 达到 $$(p,q) = (4,4)$$ ，但额外要求 $\nabla \sqrt{\rho} \in L^4_t L^4_x$ 。

而本文的 Shinbrot 型条件覆盖了更大范围的 $$(1/p,1/q)$$ 平面（见 Fig. 1 阴影区）：

当 $$q = 3$$ 时，允许 $p = \infty$ （即 $u \in L^\infty_t L^3_x$ ）；
当 $$q = 4$$ 时，允许 $$p = 2$$ （即 $u \in L^2_t L^4_x$ ），而此前工作最低需 $$p = 4$$ ；
特别地， $$(p,q) = (3,9/2)$$ 、 $$(5/2,5)$$ 等新点均被首次覆盖。

第二，对 $\partial_t \rho$ 正则性的完全解除。
所有先前工作均通过 Serrin 条件推出 $\mathrm{div}(\rho u) \in L^p_t L^q_x$ ，再结合 $\rho \in L^\infty$ 得 $\partial_t \rho \in L^p_t L^q_x$ ，从而应用 Lions 交换子估计：

\|\partial_t(f g)^\varepsilon - \partial_t(f g^\varepsilon)\|_{L^r_{t,x}} \lesssim \|\partial_t f\|_{L^{p_1}_t L^{q_1}_x} \|g\|_{L^{p_2}_t L^{q_2}_x}.

但本文 Lemma 2.5 仅需 $\partial_t f \in L^{p_1}_t W^{-1,q_1}_x$ ，而由质量方程 $\partial_t \rho = -\mathrm{div}(\rho u)$ 及 $\rho u \in L^p_t L^q_x$ ，自动有 $\partial_t \rho \in L^p_t W^{-1,q}_x$ ，无需任何额外假设。这是方法论层面的根本革新。

第三，对真空与边界的统一处理。
[9] 专门处理有界区域，但依赖强正则性；本文在相同边界条件下，用 $\phi_\delta$ 截断与 Lemma 2.3 加权估计，将真空（ $\rho=0$ ）与边界效应同时纳入弱型框架，且未牺牲判据强度。

综上，本文不仅将可压Navier–Stokes的能量守恒判据推进至Shinbrot型，更通过弱型交换子估计这一新工具，为未来研究（如可压MHD、量子Navier–Stokes）提供了可迁移的技术范式。

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

1. 从Cucker–Smale模型大尺度极限导出的动理学方程研究

On the kinetic equation arising from the large-scale limit of the Cucker-Smale model

作者: Ruicheng Cheng, Seung-Yeal Ha, Jaemoon Lee, Zhenfu Wang

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14685v1

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

摘要翻译

中文标题: 从Cucker–Smale模型大尺度极限导出的动理学方程研究

摘要翻译:
本文提出一种大尺度缩放视角，用于从相互作用粒子系统中推导介观尺度动力学，并将其应用于Cucker–Smale集群模型。与经典平均场极限（导致具有空间非局部（卷积型）对齐力的Vlasov型Cucker–Smale方程）不同，本文所采用的缩放方式导出一类新型动理学方程：其对齐场在空间变量上为局部的，仅在速度变量上保持非局部性。针对空间均匀情形，我们获得了该方程的显式解，并定量刻画了集群收敛速率；针对空间非均匀情形，我们在Sobolev空间 $W^{1,\infty}$ 及Hölder型空间 $C_b^{1,\alpha}$ 中建立了局部适定性理论，特别凸显了因缺乏卷积结构而引发的额外技术困难。进一步地，当相互作用强度足够小时，我们证明了全局适定性，并确立了解在正向时间方向上的 $$L^1$$ 渐近完备性。最后，我们研究了单速（mono-kinetic）解，构造了若干有限时间爆破（blow-up）情形。

2. 紧致李群上伪微分算子的核估计与弱 (1,1)-有界性

Kernel estimates and weak (1,1)-boundedness of pseudo-differential operators on compact Lie groups

作者: Duván Cardona, Rafik Yeghoyan, Michael Ruzhansky

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14638v1

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

摘要翻译

中文标题: 紧致李群上伪微分算子的核估计与弱 (1,1)-有界性

摘要翻译:
设 $$G$$ 为一紧致李群， $\widehat{G}$ 为其酉对偶。本文建立了定义在非交换相空间 $G \times \widehat{G}$ 上、属于全局 Hörmander 类 $S^{-n(1-\rho)/2}_{\rho,\delta}(G)$ 的伪微分算子的弱 $$(1,1)$$ 连续性（即从 $$L^1(G)$$ 到弱 $L^1(G) = L^{1,\infty}(G)$ 的有界性），其中阶数为 $-n(1-\rho)/2$ （ $n = \dim G$ ）。我们的方法核心在于推导此类算子核函数的精细估计。进一步，我们利用所得核估计，为该类算子的 $$H^1(G)$$ – $$L^1(G)$$ 连续性给出一种新的证明；该证明适用于完整的参数范围： $0 \leq \delta \leq \rho \leq 1$ ，且 $\rho \neq 0$ ， $\delta \neq 1$ 。算子所满足的条件由非交换相空间 $G \times \widehat{G}$ 上的 Hörmander 符号类 $S^m_{\rho,\delta}(G) := S^m_{\rho,\delta}(G \times \widehat{G})$ 刻画——这是欧氏空间中经典 $(\rho,\delta)$ -类在紧致李群上的自然推广。本文所有结果均在上述完整参数范围内予以陈述。
作为该有界性结果的应用，我们针对三维球面 $SU(2) \cong \mathbb{S}^3$ 上的次拉普拉斯算子 $\mathcal{L}_{\mathrm{sub}} = X^2 + Y^2$ 及热型算子 $$T = Z - X^2 - Y^2$$ ，给出了端点型先验 $$L^1$$ 估计；此类估计无法借助经典的 Hörmander 伪微分算子理论获得。更精确地，考虑如下次椭圆初值问题：
\begin{equation}\label{IVP:abstract} \begin{cases} Tu = f, & \ u,f \in \mathscr{D}’(SU(2)) := (C^\infty(SU(2)))’, & \end{cases} \end{equation}
我们证明：若 $f \in W^{1,-\frac{1}{4}}(SU(2))$ （即 $$f$$ 属于指数为 $$-1/4$$ 的 Sobolev 空间 $W^{1,p}$ 的对偶意义下的广义函数），则解 $$u$$ 必属于弱 $$L^1$$ 空间 $L^{1,\infty}(SU(2))$ 。

3. Bethe-Zel’dovich-Thompson型流体中二维定常超音速流内的斜激波相互作用

Oblique wave interactions in 2D steady supersonic flows of Bethe-Zel’dovich-Thompson fluids

作者: Geng Lai

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14144v1

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

摘要翻译

中文标题: Bethe-Zel’dovich-Thompson型流体中二维定常超音速流内的斜激波相互作用

摘要翻译: 本文研究二维半无限扩张通道内的定常超音速流动。假设流体在通道壁面上满足无滑移边界条件，且扩张通道入口处的流场状态给定。当工质为多方理想气体时，该问题可归结为若干稀疏简单波的相互作用，进而可通过特征线法证明通道内全局经典解的存在性。本文则考虑工质为非凸的Bethe-Zel’dovich-Thompson（BZT）型流体。此类流体的热力学与动力学行为可能显著区别于多方理想气体；例如，物理上容许的稀疏激波（rarefaction shocks）可能出现。依据来流状态及扩张通道张角的不同，共存在十三种互不相同的斜波相互作用构型，其中包括由激波与中心简单波复合而成的斜复合波。本文系统研究了上述各类斜波相互作用，借助特征分解法与霍多格拉夫（hodograph）变换方法，在扩张通道内构造出全局的、分片光滑的超音速解，并进一步获得了这些解的精细结构（包括各波系的几何构型、强度分布及相互作用机制等），而不仅限于解的存在性证明。本文所得结果与所采用的方法亦可推广应用于具非凸状态方程气体的若干二维黎曼问题。

4. 具有有限亏格代数几何背景的散焦非线性薛定谔方程的Painlevé XXXIV型渐近行为

Painlevé XXXIV asymptotics for the defocusing nonlinear Schrödinger equation with a finite-genus algebro-geometric background

作者: Engui Fan, Gaozhan Li, Yiling Yang, Lun Zhang

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

发布日期: 2026-02-14

arXiv ID: 2602.13589v1

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

摘要翻译

中文标题: 具有有限亏格代数几何背景的散焦非线性薛定谔方程的Painlevé XXXIV型渐近行为

摘要翻译: 本文研究具有有限亏格代数几何背景的散焦非线性薛定谔方程（NLS）的柯西问题。我们在四个不同的时空区域中导出了该初值问题解的长时间渐近展开。结果表明：渐近展开的首项（至多相差一个常数因子）即为背景解，但其参数发生了一个平移；而次主导项在不同区域中则以不同的速率衰减。特别值得注意的是，在两个过渡区域中，该次主导项的量级为 $\mathcal{O}(t^{-1/3})$ ，且其系数显式地依赖于Painlevé XXXIV型超越函数的一个积分。我们的结论是通过对相关黎曼–希尔伯特问题实施非线性最速下降法分析而严格建立的。

不相关论文（基本信息）

1. 平均场Kuramoto-Daido模型的相变与线性稳定性分析

Phase transitions and linear stability for the mean-field Kuramoto-Daido model

作者: Kyunghoo Mun, Matthew Rosenzweig

类别: math.AP, math-ph

发布日期: 2026-02-16

arXiv ID: 2602.14954v1

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

摘要: We consider the mean-field noisy Kuramoto-Daido model, which is a McKean-Vlasov equation on the circle with bimodal interaction $W(θ)=\cosθ+m\cos2θ$ for $m\ge 0$ and interaction strength $$K$$ , generalizing the celebrated noisy Kuramoto model corresponding to $$m=0$$ . Our first contribution is to characterize the phase transition threshold $K_{c}$ by comparing it to the linear stability threshold $K_\# = \min (1, m^{-1})$ of the uniform distribution. When $m \leq 1/2,$ $K_{c}=1$ , coinciding with that of the Kuramoto model. On the other hand, for $m \geq 2$ , we show $K_c= m^{-1}$ . We also classify the regimes in which the phase transition is continuous or discontinuous. Our second contribution is to analyze the linear stability of a global minimizer $$q$$ (the ordered phase'') of the mean-field free energy in the supercritical regime $K>1$. This stationary solution of the Kuramoto-Daido equation is unique up to translation invariance and distinct from the uniform distribution (the disordered phase’’). Our approach extends the Dirichlet form method of Bertini et al. from the unimodal to bimodal setting. In particular, for $m \leq 1.590 \times 10^{-4}$ and $$K>1$$ , we show an explicit lower bound on the spectral gap of the linearized McKean-Vlasov operator at $$q$$ . To our knowledge, this is the first rigorous stability analysis for this class of models with bimodal interactions.

2. 关于Hessian商方程 $$σ_2/σ_1$$ 的刘维尔定理

On Liouville’s theorem for the Hessian quotient equation $$σ_2/σ_1$$

作者: Siyuan Lu, Marcin Sroka

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14946v1

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

摘要: We prove Liouville’s theorem for semi-convex entire solutions to Hessian quotient equation $$σ_2/σ_1=1$$ in $\mathbb{R}^n$ . The proof is based on the observation that after rewriting the quotient operator as the $$σ_2$$ operator, acting on a new function, one can refer to the recent result of Shankar and Yuan on Liouville’s theorem for $$σ_2$$ equation.

3. 大初值问题下热-压电耦合系统的解的存在性及相应反问题的正向算子分析

Existence of large-data solutions to a thermo-piezoelectric system and forward operator analysis for associated inverse problems

作者: Torben J. Fricke, Raphael Kuess, Felix Meyer

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14848v1

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

摘要: We consider an inverse problem governed by the initial-boundary value problem for the thermoviscoelastic Kelvin-Voigt system \begin{align*}\left{ \begin{array}{l} ρ(z,t) u_{tt}- \left(Γ(Θ) u_{zt} +p(z,t) u_z + βΘ_z\right)z=0\ b(z,t) Θ_t-\left(k(z,t)Θ_z\right)z - Γ(Θ) u{zt}^2+βΘu{zt}=0, \end{array} \right. \end{align*} in an open bounded interval $Ω\subset\mathbb{R}$ , for the evolution of the displacement variable $$u$$ , and the temperature $Θ\geq 0$ . Assuming the material coefficients $$ρ$$ , $$Γ$$ , $$p$$ , $$b$$ and $$k$$ are strictly positive and bounded and $β\in\mathbb{R}$ a global-in-time existence result is established for weak solutions. The present manuscript demonstrates that this can be achieved under energy- and entropy-minimal assumptions, in the sense that global weak solutions are shown to exist for any initial data

u_0\in W^{1,2}(Ω),\quad u_{0t}\in L^2(Ω)\quad\text{and}\quad 0\leΘ_0\in L^2(Ω).

The qualitative analysis of the evolution problem then allows to model and analyze the structural properties of the corresponding forward operator that naturally arises in inverse parameter identification settings. Therein, two modeling approaches of the observation operator as approximations of the electrical surface charge are presented and results on their well-definedness and boundedness are established. With the results on well-definedness and boundedness of the model operator, established in this paper as well, results on well-definedness, boundedness and continuous Fréchet differentiability of the forward operator are presented.

4. 从区间上识别Bergman空间函数

Identifying Bergman space functions from intervals

作者: Andreas Hartmann, Marcu-Antone Orsoni

类别: math.CV, math.AP, math.FA, math.OC

发布日期: 2026-02-16

arXiv ID: 2602.14801v1

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

摘要: We characterize functions of a Bergman space on a square by their values and derivatives on the diagonals. This problem is connected with the reachable space of the one-dimensional heat equation on a finite interval with boundary $$L^2$$ -controls.

5. 关于混合退化/奇异情形下双非线性抛物方程的变号解的Hölder连续性

On the Hölder continuity of signed solutions to doubly nonlinear parabolic equations in the mixed degenerate/singular cases

作者: Igor I. Skrypnik

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14800v1

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

摘要: We prove the Hölder continuity of sign-changing solutions to the equation of the type

\frac{\partial}{\partial t}\big(|u|^{q-1} u\big)- div\Big(|D u|^{p-2}\,D u\Big)=0,

where numbers $$p$$ , $$q$$ satisfy the conditions

0<q<p-1\quad \text{and}\quad p<2,

q>p-1\quad\text{and}\quad p>2.

Our proof uses new versions of the integral Harnack type inequalities for sign-changing solutions.

6. 极端Reissner-Nordström时空上标量波动方程散射解的匹配条件

Matching conditions for scattering solutions of scalar wave equations on extremal Reissner-Nordström spacetimes

作者: Yannis Angelopoulos, Istvan Kadar

类别: math.AP, gr-qc

发布日期: 2026-02-16

arXiv ID: 2602.14712v1

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

摘要: We study scattering solutions $$φ$$ of the linear wave equation on extremal Reissner-Nordström spacetimes, satisfying the following properties: i) $$φ$$ attains a prescribed radiation field $ψ_{\mathcal{I}}$ through future null infinity, which decays at an inverse polynomial rate; ii) $$φ$$ is regular in the exterior region up to and including the future event horizon, i.e. $φ\in C^N$ , where $N\gg1$ is independent of the decay rate of $ψ_{\mathcal{I}}$ . We prove that such solutions exist for arbitrary $$N$$ , and that they are not unique. The proof consists of: 1) finding an approximate solution $φ_{\mathrm{app}}$ with fast decaying error; 2) the use of backwards energy estimates in order to correct $φ_{\mathrm{app}}$ to an exact solution. Extremality is used only in the second step. The methods of the linear case described above are then used to show the same results for semilinear equations where the nonlinearity satisfies the null condition, as well as to geometries describing the hyperbolic orbit of multiple extremal Reissner-Nordström black holes.

7. 非线性变分传输问题中梯度的Morrey估计

Morrey estimates for the gradient in non-linear variational transmission problems

作者: Luca Esposito, Lorenzo Lamberti

类别: math.AP, math.OC

发布日期: 2026-02-16

arXiv ID: 2602.14658v1

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

摘要: We study a class of variational transmission problems driven by nonlinear energies with discontinuous coefficients across a prescribed interface. The model setting consists of integral functionals of the form

\mathcal{F}(u;E)=\int_Ωσ_E(x)\,F(\nabla u)\,dx,

where the coefficient $$σ_E$$ takes two constant values on complementary regions separated by a $$C^1$$ hypersurface, and the integrand $$F$$ satisfies standard $$p$$ -growth and monotonicity conditions with $$p>2$$ . In this nonlinear variational framework, we establish local Morrey-space regularity for the gradient of local minimizers, proving that $\nabla u\in L^{2,λ}_{\mathrm{loc}}(Ω)$ for every $0\leqλ<n$ , provided $2<p<\frac{2n}{n-2}$ . The proof is based on quantitative decay estimates for the energy near the interface, first obtained in a flat configuration and then extended to the general case by a suitable approximation argument.

8. Weingarten曲面在无穷远点的渐近行为

Asymptotic behavior at infinity of Weingarten surfaces

作者: Aires E. M. Barbieri, José A. Gálvez, Yuanyuan Lian, Kai Zhang

类别: math.DG, math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14651v1

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

摘要: We derive the asymptotic expansion at infinity for embedded ends of uniformly elliptic Weingarten surfaces with finite total curvature in $\mathbb{R}^3$ , and we establish a maximum principle at infinity. Furthermore, we solve the Dirichlet problem for the uniformly elliptic Weingarten equation in dimension two on strictly convex bounded domains.

9. 基于热等离子体纳米颗粒作执行器的热方程反馈镇定与跟踪控制

Feedback Stabilization and Tracking for Heat Equations Using Thermo-Plasmonic Nanoparticles as Actuators

作者: Arpan Mukherjee, Sérgio S. Rodrigues, Mourad Sini

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14581v1

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

摘要: We propose a feedback strategy to track prescribed heat profiles using plasmonic nanoparticles as actuators. Starting from a thermo–plasmonic Maxwell–heat model, we use a time-domain discrete effective description in which the generated heat is approximated by a superposition of heat kernels centered at particle locations with amplitudes governed by a coupled Volterra system. We recast this dynamics as a heat equation on a bounded domain with finitely many point actuators and design a tracking feedback based on pointwise evaluations of $\mathcal A^{-1}y$ , where $\mathcal A=I-A_0$ and $$A_0$$ is the Neumann diffusion operator. Working in the natural $$V'$$ setting with $V=D(\mathcal A)$ , we prove exponential stabilization of the tracking error via distribution-actuator theory. For non-equilibrium reference profiles, we add a constant feedforward term and a low-mode fixed-point pre-compensation on $$X_N$$ , ensuring exact steady matching on $$X_N$$ and an explicit bound on the residual tail mismatch.

10. 正则性沿不稳定流形的传播

Propagation of regularity along unstable manifolds

作者: Thibault Lefeuvre, Rafael Potrie

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

发布日期: 2026-02-16

arXiv ID: 2602.14558v1

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

摘要: Let $\varphi_t : M \to M$ be a flow on a smooth closed connected manifold $$M$$ that preserves and expands a foliation $$F$$ . We establish a theorem of propagation of regularity along the leaves of $$F$$ for sections of vector bundles satisfying a transport equation involving the generator of a cocycle over $\varphi_t$ . As a consequence, we prove a regularity result for Pollicott-Ruelle resonant states: if such state is smooth in restriction to a piece of an unstable leaf, then it is in fact smooth over the entire manifold. We also announce further applications related to joint integrability of extreme bundles of partially hyperbolic diffeomorphisms. The proofs rely on a leafwise semiclassical pseudodifferential calculus adapted to a foliated space, which may be of independent interest.

11. 关于一类非线性薛定谔型方程的第一特征值

On the first eigenvalue of a nonlinear Schrödinger type equation

作者: Ardra A

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14545v1

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

摘要: We consider an eigenvalue problem for the generalized nonlinear Schrödinger type operator with the Robin boundary condition as given below. \begin{equation*} \label{ab-Robin p-Laplace evp with potential term_intro} \left{ \begin{split} -Δ_p u+V(x)|u|^{p-2}u&=λ|u|^{p-2}u\quad &&\mathrm{in} ~~Ω, |\nabla u|^{p-2}\frac{\partial u}{\partialη}+β|u|^{p-2}u&=0\quad &&\mathrm{on}~~\partialΩ, \end{split} \right. \end{equation*} where $Δ_p u := \operatorname{div}(|\nabla u|^{p-2}\nabla u)$ is the $$p$$ -Laplace operator, $$Ω$$ is a bounded domain in $\mathbb{R}^n$ with smooth boundary, $V \in C^1(\mathbb{R}^n),$ $$ η$$ denotes the outward unit normal, and $$ β$$ is a positive real constant. We study the properties of its first eigenvalue with respect to the potential $$V$$ , the boundary parameter $$β$$ as well as the domain. First, we establish some properties of the smallest eigenvalue $$λ_1(V)$$ with respect to the potential. We then prove the differentiability of $$λ_1(V)$$ with respect to the Robin boundary parameter $$β$$ and give an explicit formula for this derivative, which is then used to investigate some monotonicity properties of $$λ_1(V).$$ We also obtain a shape derivative formula for the smallest eigenvalue. Using these derivatives, we also study domain monotonicity properties of the first eigenvalue.

12. Ricci曲率有界条件下的盖尔范德反问题

Gel’fand’s inverse problem under Ricci curvature bounds

作者: Shouhei Honda, Jinpeng Lu

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

发布日期: 2026-02-16

arXiv ID: 2602.14527v1

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

摘要: The classical Gel’fand’s inverse problem asks whether a Riemannian manifold is uniquely determined by the knowledge of the heat kernel on any open subset of the manifold. We study this inverse problem in the non-smooth setting in the framework of ${\rm RCD}(K,N)$ spaces, namely, metric-measure spaces with synthetic Riemannian Ricci curvature bounded below by $$K$$ and dimension bounded above by $$N$$ . We establish the unique solvability of Gel’fand’s inverse problem for the class of compact ${\rm RCD}(K,N)$ spaces whose regular set admits $$C^1$$ -Riemannian structure. As an application, we obtain the stability of Gel’fand’s inverse problem in the class of closed Riemannian manifolds with bounded Ricci curvature, diameter and volume bounded from below. We note that the results are new even for Einstein orbifolds and (weighted) Riemannian manifolds with non-smooth boundary.

13. 具有阿哈罗诺夫–玻姆势的磁诺伊曼问题：特征值的边界渐近行为与分裂现象

Magnetic Neumann problems with Aharonov-Bohm potentials: boundary asymptotics of eigenvalues and splitting phenomena

作者: Veronica Felli, Prasun Roychowdhury, Giovanni Siclari

类别: math.AP

发布日期: 2026-02-16

arXiv ID: 2602.14522v1

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

摘要: We study a planar magnetic Schrödinger operator with an Aharonov-Bohm vector potential, under Neumann boundary conditions. Through a gauge transformation, the corresponding eigenvalue problem can be formulated in terms of the Laplacian on a fractured domain, where the fracture lies along the segment connecting the pole to its projection on the boundary. As the pole approaches the boundary, we prove that the eigenvalues converge to those of the Neumann Laplacian and the variation exhibits a logarithmic vanishing rate. In the case of multiple eigenvalues, when the pole approaches a fixed point of the boundary, we observe a splitting phenomenon, with the largest branch separating from the others.

14. 矩阵加权的庞加莱型不等式及其在齐型空间上对数Hajłasz–Besov空间中的应用

Matrix-Weighted Poincaré-Type Inequalities with Applications to Logarithmic Hajłasz–Besov Spaces on Spaces of Homogeneous Type

作者: Ziwei Li, Dachun Yang, Wen Yuan

类别: math.FA, math.AP, math.CA

发布日期: 2026-02-16

arXiv ID: 2602.14361v1

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

摘要: Let ${\mathcal {X}}$ be a space of homogeneous type. In this article, based on the reducing operators of matrix $$A_p$$ -weights, the authors introduce the vector-valued Hajłasz gradient sequences and establish some related matrix-weighted Poincaré-type inequalities on ${\mathcal {X}}$ . As an application, the authors introduce the matrix-weighted logarithmic Besov spaces on ${\mathcal {X}}$ and establish their pointwise characterization via Hajłasz gradient sequences. The novelty of this article lies in that, by means of both the $$A_p$$ dimension and its properties of matrix $$A_p$$ -weights and the wavelet reproducing formula with exponential decay of P. Auscher and T. Hytönen, all the main results get rid of the dependence on the reverse doubling conditions of both weights and ${\mathcal {X}}$ under consideration and these results are also completely new even for unweighted logarithmic Besov spaces on ${\mathcal {X}}$ .

15. 由闭合型微机电系统模型导出的分数阶椭圆方程解的存在性与非存在性

Existence and nonexistence of solutions for fractional elliptic equations arising from closed MEMS model

作者: Huyuan Chen, Jialei Jiang, Jun Wang

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14346v1

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

摘要: The objective of our paper is to investigate fractional elliptic equations of the form $(-Δ)^s u=\frac{λ}{(a-u)^2}$ within a bounded domain $$Ω$$ , subject to zero Dirichlet boundary conditions. Here, $s\in(0,1)$ , $$λ>0$$ , and the function $$a$$ vanishes at the boundary while satisfying additional conditions. This problem originates from Micro-Electromechanical Systems (MEMS) devices, particularly when the elastic membrane makes contact with the ground plate at the boundary. We establish both existence and nonexistence results, illustrating how the boundary decay of the membrane influences the solutions and pull-in voltage.

16. 非紧双曲曲面上薛定谔方程的可观测性与半经典控制

Observability and Semiclassical Control for Schrödinger Equations on Non-compact Hyperbolic Surfaces

作者: Xin Fu, Yulin Gong, Yunlei Wang

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

发布日期: 2026-02-15

arXiv ID: 2602.14316v1

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

摘要: We study the observability of the Schrödinger equation on $$X$$ , a non-compact covering space of a compact hyperbolic surface $$M$$ . Using a generalized Bloch theory, functions on $$X$$ are identified as sections of flat Hilbert bundles over $$M$$ . We develop a semiclassical analysis framework for such bundles and generalize the result of semiclassical control estimates in [Dyatlov and Jin, Acta Math., 220 (2018), pp. 297-339] to all flat Hilbert bundles over $$M$$ , with uniform constants with respect to the choice of bundle. Furthermore, when the Riemannian cover $X \to M$ is a normal cover with a virtually Abelian deck transformation group $$Γ$$ , we combine the uniform semiclassical control estimates on flat Hilbert bundles with the generalized Bloch theory to derive observability from any $$Γ$$ -periodic open subsets of $$X$$ . We also discuss applications of the uniform semiclassical control estimates in spectral geometry.

17. 子解梯度的局部正则性

Partial regularity of the gradient for subsolutions

作者: Aram Hakobyan, Michael Poghosyan, Henrik Shahgholian

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14305v1

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

摘要: We prove that the gradient of any bounded subharmonic function is upper semi-continuous, provided that its super-level sets can be touched from the exterior by uniform $C^{1,\text{Dini}}$ domains at every point. This idea extends to a class of general operators, as well as to the boundary behaviour of the gradient of solutions of the Dirichlet problem in a domain whose boundary satisfy this geometric condition.

18. 吸引-排斥趋化模型中的非局部物流项与非线性产生项：整体适定性分析

Nonlocal logistics and nonlinear productions in an attraction-repulsion chemotaxis model: analysis of the global well-posedness

作者: Rafael Díaz Fuentes, María Victoria Redondo Neble, Giuseppe Viglialoro

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14227v1

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

摘要: This paper investigates a {three-component} chemotaxis system involving both attraction and repulsion effects, as well as a nonlocal logistic-type source term. Mathematically, if $$u=u(x,t)$$ , $$v = v(x,t)$$ and $$w = w(x,t)$$ denote the cell distribution, and the attractive and the repulsive chemical signals, the model is then described by \begin{equation*} \begin{cases} u_t = Δu - χ\nabla \cdot (u \nabla v) + ξ\nabla \cdot (u \nabla w) + a u^α- b u^α\int_Ωu^β, & x \in Ω, \ t > 0, τv_t = Δv - v + f(u), & x \in Ω, \ t > 0, τw_t = Δw - w + g(u), & x \in Ω, \ t > 0. \end{cases} \end{equation*} Here, $Ω\subset \mathbb{R}^n$ ( $n \geq 1$ ) is a bounded smooth domain, $τ\in\{0,1\}$ , $$a,b,α,β,χ,ξ>0$$ , the production functions $$f(u)$$ and $$g(u)$$ are assumed to satisfy algebraic growth conditions of order $\ell$ and $$ρ$$ , generalizing prototypes of the form $u^\ell$ and $$u^ρ$$ , $\ell,ρ>0$ . The work is devoted to proving the global existence and boundedness of classical solutions under a suitable balance between the signal production exponents $\ell, ρ$ and the nonlocal damping exponents $$α, β$$ , for regular enough initial data and zero-flux boundary restrictions. In this regard, two main theorems are established for the cases where the chemical signals satisfy either elliptic ( $$τ=0$$ ) or parabolic ( $$τ=1$$ ) partial differential equations, highlighting how sufficiently strong nonlocal damping prevents the formation of singularities in time. We extend the results obtained in [Chiyo et al., Appl. Math. Optim. 89:9 (2024)], where the fully parabolic ( $$τ=1$$ ) and only attraction version is studied. In our context, we establish well-posedness of the system and the long-time behavior of solutions.

19. 双曲空间与黎曼流形上的对数索博列夫不等式、庞加莱不等式及贝克纳不等式

Logarithmic Sobolev, Poincaré and Beckner Inequalities on Hyperbolic Spaces and Riemannian Manifolds

作者: Anh Xuan Do, Debdip Ganguly, Nguyen Lam, Guozhen Lu

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

发布日期: 2026-02-15

arXiv ID: 2602.14202v1

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

摘要: We investigate several functional and geometric inequalities on the hyperbolic space $\mathbb{H}^N$ , with a primary emphasis on logarithmic Sobolev inequalities, Poincaré inequalities, and Beckner-type inequalities, all studied within the framework of the AB program. The main analytical tool employed throughout this paper is symmetrization. More precisely, our approach relies on an improved version of the Pólya-Szegö inequality on the hyperbolic space, obtained through a careful comparison of the gradient norms of rearranged functions in the hyperbolic and Euclidean settings. For Beckner-type inequalities, we adopt a semigroup approach based on sharp estimates for the heat semigroup, leading to refined interpolation inequalities between Poincaré and logarithmic Sobolev inequalities. Finally, we extend our results beyond hyperbolic space to a class of Riemannian model manifolds $\mathbb{M}^N$ satisfying the centered isoperimetric inequality. This shows that the inequalities and methods developed in this work are robust and rely mainly on geometric and isoperimetric properties, rather than on the specific structure of hyperbolic space itself.

20. 关于共形 $$Q$$ -曲率方程的刘维尔型定理

Liouville theorems for conformal $$Q$$ -curvature equations

作者: Meiqing Xu, Hui Yang

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14072v1

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

摘要: In this paper, we study the non-existence of positive solutions for the following conformal $$Q$$ -curvature equation \begin{equation*} (-Δ)^σu = K(x) u^{\frac{n+2σ}{n-2σ}} \quad \text{in } \mathbb{R}^n, \end{equation*} where $σ\in (0, n/2)$ is a real number. When $$σ=1$$ , this equation reduces to the well-known scalar curvature equation arising from the prescribed scalar curvature problem. For general $σ\in (0, n/2)$ , it appears in the study of prescribing $$Q$$ -curvature. We establish Liouville theorems under various assumptions on the $$Q$$ -curvature $$K(x)$$ by developing a unified approach applicable to all $σ\in (0, n/2)$ . Our method successfully addresses the challenges posed by the absence of ODE tools in the fractional regime and the lack of a classification of Delaunay-type singular solutions for the general fractional Yamabe equation.

21. 具有退化扩散的离散非线性碎裂方程解的存在性

Existence for the Discrete Nonlinear Fragmentation Equation with Degenerate Diffusion

作者: Saumyajit Das, Ram Gopal Jaiswal

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14070v1

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

摘要: A mathematical model for the discrete nonlinear fragmentation (collision-induced breakage) equation with diffusion is studied. The existence of global weak solutions is established in arbitrary spatial dimensions without assuming a strictly positive lower bound on the diffusion coefficients, extending previous results that were restricted to one-dimensional domains and relied on uniformly positive diffusion. The analysis is carried out under boundedness assumptions on the collision and breakage kernels. The proof is based on the construction of a suitable regularized system, combined with weak $$L^2$$ a priori estimates and compactness arguments in $$L^1$$ , which allow the passage to the limit in the nonlinear fragmentation operator.

22. 三维空间中Hessian商方程的内点Hessian估计

Interior Hessian estimates for Hessian quotient equations in dimension three

作者: Heming Jiao, Zhenan Sui

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.14064v1

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

摘要: In this paper, we establish the interior Hessian estimates for $$2$$ -convex solutions to $\frac{σ_2}{σ_1} (D^2 u) = ψ(x,u)$ in dimension three. In higher dimensions ( $n \geq 4$ ), we prove the interior Hessian estimates for semi-convex solutions. We provide a new method to prove the doubling inequality for smooth solutions in dimensions three and four. In higher dimensions ( $n\geq 5$ ) the doubling inequality is proved under an additional dynamic semi-convexity condition which is the same to that in \cite{SY2025}. The method also applies to the equation $σ_2 (D^2 u) = ψ(x, u, \nabla u)$ .

23. 分数阶合作系统的谱理论与传染病模型中的阈值动力学

Spectral Theory of Fractional Cooperative Systems and Threshold Dynamics in Epidemic Models

作者: Cong-Bang Trang, Hoang-Hung Vo

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.13986v1

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

摘要: Spectral analysis has long been recognized as a fundamental tool for studying the existence, uniqueness, and qualitative behavior of solutions to semilinear elliptic and parabolic equations, as well as their long-time dynamics. In modern mathematics, fractional Laplacians are widely used to model nonlocal or long-range diffusion processes arising in biology, including anomalous movement, long-distance dispersal, and Levy-flight migration of organisms, cells, and epidemics. In this paper, we employ the spectral fractional Laplacian introduced by Caffarelli and Stinga (2016) to develop the eigentheory for a cooperative system describing an infectious epidemic process and to analyze its long-term behavior. Using Fredholm theory and related analytical techniques, building in part on ideas of Lam and Lou (2016), we establish a sharp criterion ensuring the existence and simplicity of the principal eigenvalue, together with variational characterizations and consequences for the validity of maximum principles. We further derive the asymptotic behavior of the principal eigenvalue with respect to diffusion coefficients, fractional orders, and domain scaling, complementing recent developments by Zhao and Ruan (2023) and Feng, Li, Ruan, and Xin (2024). As an application of this spectral framework, we prove the existence, uniqueness, and threshold-type long-time dynamics of solutions to an endemic reaction-diffusion system with fractional diffusion, providing a perspective that differs from earlier approaches such as Hsu and Yang (2013). Our results contribute to the growing interaction between spectral theory and nonlocal analysis, in line with recent advances in the area.

24. 扰动晶格图上的约束变分问题

Constrained variational problems on perturbed lattice graphs

作者: Weiqi Guan

类别: math.AP

发布日期: 2026-02-15

arXiv ID: 2602.13978v1

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

摘要: In this paper, we solve some constrained variational problems on perturbed lattice graphs $$G$$ . The first problem addresses the existence of ground state normalized solutions to Schrödinger equations \begin{equation*} \left{ \begin{aligned} &-Δ_{G} u+λu=\vert u\vert^{p-2}u,x\in G &\Vert u\Vert_{l^2(G)}^2=a. \end{aligned} \right. \end{equation*} We prove that if the graph is obtained by deleting finite edges in lattice graphs while maintaining connectivity, then there exists a threshold $α_G\in[0,\infty)$ such that there do not exist ground state normalized solution if $$0<a<α_G$$ , and there exists a ground state normalized solution if $$a>α_G.$$ If the graph is obtained by adding finite edges $E^{'}$ to lattice graphs, we prove that there exist $E^{'}$ and $$a_1$$ such that for all $$a>a_1,$$ there do not exist ground state normalized solutions. The second problem concerns the existence of an extremal function for the Sobolev inequality. If the graph $$G$$ is obtained by deleting finite edges in lattice graphs while maintaining connectivity, for the Sobolev super-critical regime, we prove that there exists an extremal function. for the Sobolev critical regime, we prove that there exists $$G$$ such that extremal can be attained. If the graph is obtained by adding finite edges $E^{'}$ to lattice graphs, we prove that there exists $E^{'}$ such that there does not exist an extremal function.

25. 复Sobolev空间中函数的精细性质

Fine properties of functions in complex Sobolev spaces

作者: Ngoc Cuong Nguyen

类别: math.CV, math.AP

发布日期: 2026-02-15

arXiv ID: 2602.13968v1

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

摘要: We study comprehensively local properties of functions in complex Sobolev spaces on a bounded open subset of $\mathbb{C}^n$ . The main tool is the corresponding functional capacity for the space which is inspired by the global one due to Vigny (2007). An inequality between this capacity and the Bedford-Taylor capacity for plurisubharmonic functions is proved, which is sharp as far as the exponents are concerned. Moreover, it is shown that the functional capacity is a Choquet capacity. The Alexander-Taylor type inequality for the capacity is also proved. This allows us to strengthen the results in the works of Dinh, Marinescu and Vu (2023), Vigny and Vu (2024). Lastly, the Moser-Trudinger type inequality in this space is characterized by the volume-capacity inequality.

26. 涉及非局部算子的Lane–Emden型不等式的刘维尔型定理

Liouville-type theorems for Lane–Emden inequalities involving nonlocal operators

作者: T. Kim, T. Lee

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13822v1

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

摘要: We establish a Liouville-type theorem for nonnegative weak supersolutions to $\mathcal{L}_K u = u^q$ in $\mathbb{R}^n$ , where $\mathcal{L}_K$ is a translation-invariant integro-differential operator of order $$2s$$ with $s \in (0,1)$ . The kernel $$K$$ is assumed to be even and satisfy uniform ellipticity bounds. We prove that the only nonnegative supersolution is the trivial one $u \equiv 0$ in the range $1 < q \le \frac{n}{n-2s}$ for $$n > 2s$$ (and for all $$q > 1$$ when $n \le 2s$ ). Our proof is elementary and relies on a test function method combined with a dyadic decomposition of the nonlocal tail. Notably, our argument does not rely on the maximum principle or the fundamental solution.

27. 强竞争非线性薛定谔系统的整体解

Entire solutions to a strongly competitive nonlinear Schrödinger system

作者: Pierpaolo Esposito, Pablo Salgado, Angela Pistoia, Giusi Vaira

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13753v1

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

摘要: We build infinitely-many non-radial positive solutions to the Schrödinger system \begin{equation*} \left{\begin{aligned} &-Δu_1+u_1=u_1^{{\mathfrak p} }-Λu_1^{a_1} u_2^{a_2}\ \hbox{in}\ \mathbb R^N\ &-Δu_2+u_2=u_2^{{\mathfrak p} }-Λu_1^{b_1}u_2^{b_2} \ \hbox{in}\ \mathbb R^N\ \end{aligned}\right. \end{equation*} with sub-critical $\mathfrak p$ -growth as $Λ\to +\infty$ . The profile of each component is the sum of several copies of the positive solution to $-ΔU+U=U^{{\mathfrak p} }$ in $\mathbb R^N$ , centered at suitable {\em peaks} whose mutual distances diverge as $$Λ$$ increases. More precisely, given two concentric regular polygons with $$k$$ sides and very large radii, the peaks of the first component are arranged along the edges of the {\em outer} polygon, alternated with those of the second component, and along the $$k$$ rays joining the vertices of the two polygons. To the best of our knowledge, this provides the first example of non-radial positive solutions for strongly competitive Schrödinger systems in the whole space.

28. 基于预解算子方法的定常福克–普朗克方程中椭圆正则性研究

Resolvent approaches to elliptic regularity in stationary Fokker-Planck equations

作者: Haesung Lee

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13727v1

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

摘要: This paper investigates the local regularity of solutions to stationary Fokker-Planck equations on an open set $U \subset \mathbb{R}^d$ with $d \geq 2$ . A central objective is to relax the classical assumptions on the coefficients by focusing on the case where the drift vector field $\mathbf{G}$ is only assumed to be locally square-integrable, i.e. $\mathbf{G} \in L^2_{loc}(U, \mathbb{R}^d)$ , the diffusion matrix $A = (a_{ij})_{1 \leq i,j \leq d}$ is assumed to be locally uniformly strictly elliptic and bounded, with coefficients satisfying $\frac{a_{ij}+a_{ji}}{2} \in VMO_{loc}(U)$ for all $1 \leq i,j \leq d$ and ${\rm div}A \in L^2_{loc}(U, \mathbb{R}^d)$ . Our main result shows that any locally bounded function $h \in L^\infty_{loc}(U)$ satisfying the stationary Fokker-Planck equation $L^*(h\, dx) = 0$ must in fact belong to the local Sobolev space $H^{1,2}_{loc}(U)$ . The proof is based on the construction of a sub-Markovian resolvent associated with the principal elliptic operator $L^A := \operatorname{trace}(A \nabla^2)$ , combined with delicate energy-type inequalities. In particular, we show that the density $$h$$ can be realized as the weak $H^{1,2}$ -limit of images of resolvent operators.

29. 各向异性Calderón问题在Painlevé-Liouville型黎曼流形上的稳定性

Stability in the anisotropic Calderón problem for Painlevé-Liouville Riemannian manifolds

作者: Thierry Daudé, Niky Kamran, François Nicoleau

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

发布日期: 2026-02-14

arXiv ID: 2602.13703v1

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

摘要: We study the question of stability of the global and partial anisotropic Calderón inverse problems for the class of Painlevé-Liouville Riemannian manifolds, that is compact $$n$$ -dimensional manifolds with boundary $$(M,g)$$ , where $M=[0,1]\times K\,$ , $$K$$ is any smooth closed connected orientable manifold of dimension $$n-1$$ endowed with a Riemannian metric $$g_K$$ , and $g=α^4 g_{0}$ is any conformal deformation of the product metric $g_{0}=dx^2+g_{K}$ on $$M$$ which is compatible with the Painlevé block-separability of the Laplace-Beltrami operator $Δ_{g_0}$ . Given a pair of Painlevé-Liouville Riemannian manifolds $$(M,g)$$ and $(M,\tilde{g})$ satisfying some technical hypothesis, denoting the corresponding Dirichlet-to-Neumann maps by $Λ_{g}$ and $Λ_{\tilde{g}}$ , and assuming that $\lVert Λ_{g}-Λ_{\tilde{g}}\rVert_{\mathcal{B}(H^{1/2}(\partial M), H^{-1/2}(\partial M))}\ = ε$ , we show a logarithmic stability result for the global anisotropic Calderón problem which says that there exists constants $$C$$ and $$0<θ<1$$ such that $\| α- \tildeα \|_{C^{0,r}(M)} \leq C \left( \ln \frac{1}ε \right)^{-θ}$ for some $$0<r<1$$ . Similar results are obtained for the partial anisotropic Calderón problem, corresponding to the case where the data are measured on only one connected component of the boundary.

30. 基于边界控制法的时间域—频域波边界数据转换公式

A Formula for Time-to-Frequency Wave Boundary Data Conversion by the Boundary Control Method

作者: Yang Yang

类别: math.AP

发布日期: 2026-02-14

arXiv ID: 2602.13612v1

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

摘要: Given the wave equation on a compact Riemannian manifold with boundary, we derive an explicit reconstruction procedure to represent the frequency-domain Neumann-to-Dirichlet map in terms of the time-domain Neumann-to-Dirichlet map at any non-eigenfrequency. If the wave equation is exactly controllable, we derive an explicit formula to compute the former from the latter. The derivation is based on the boundary control method and requires only knowledge on the boundary of the manifold. The formula is stable when the level of regularization is fixed. The numerical feasibility is validated using one-dimensional examples in both Euclidean and non-Euclidean geometries.

31. 具有局部奇异结构阻尼的两个耦合欧拉-伯努利方程的正则性与稳定性

Regularity and stability of two coupled Euler-Bernoulli equations with a localized singular structural damping

作者: K. Ammari, F. Hassine, L. Tebou

类别: math.AP

发布日期: 2026-02-13

arXiv ID: 2602.13525v1

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

摘要: This paper is concerned with the study of regularity and stability properties of two Euler-Bernoulli beam equations with localized singular damping. Under suitable regularity assumptions on the damping coefficient, we establish Gevrey regularity for the semigroup generated by the associated operator. Furthermore, for a broader class of damping mechanisms, including less regular damping, we derive uniform stability result. These findings provide a detailed description of the long-term behavior of the corresponding dynamical systems.

32. 球对称黎曼流形上的径向索博列夫嵌入

Radial Sobolev embeddings on spherically symmetric Riemannian manifolds

作者: João Marcos do Ó, Guozhen Lu, Raoní Ponciano

类别: math.AP

发布日期: 2026-02-13

arXiv ID: 2602.13501v1

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

摘要: We study Sobolev spaces of radial functions on spherically symmetric Riemannian manifolds. Using geodesic polar coordinates, we give a sharp one-dimensional reduction: a radial function belongs to the Sobolev space on the manifold if and only if its radial representation lies in an associated weighted Sobolev space on an interval, with weights determined explicitly by the metric. This characterization allows us to prove optimal Sobolev-type embeddings for radial functions into weighted Lebesgue spaces on both bounded and unbounded spherically symmetric manifolds. As further consequences, we establish new radial lemmas and decay estimates that capture the precise behaviour of radial Sobolev functions near the origin and at infinity. Our results unify and extend the classical radial embeddings in Euclidean and hyperbolic spaces.

33. 全非线性非局部问题的强极大值原理

Strong maximum principle for fully nonlinear nonlocal problems

作者: Juan Pablo Cabeza, Gabrielle Nornberg, Disson dos Prazeres

类别: math.AP

发布日期: 2026-02-13

arXiv ID: 2602.13425v1

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

摘要: In this paper, we study solvability and qualitative properties of nonnegative solutions for a sublinear nonlocal problem with fully nonlinear structure in the form

\mathcal{M}^{\pm}[u]+a(x)u^{q}(x)=0 \; \text{ in }Ω,\qquad u\geq 0 \; \text{ in }Ω.

Here $Ω\subset \mathbb{R}^n$ is a bounded $C^{1,1}$ convex domain, $\mathcal{M}^{ \pm}$ stands for nonlocal Pucci extremal operators defined in a class $\mathcal{L}_*$ of homogeneous kernels, $q\in(0,1)$ , and $$a$$ is a possibly sign-changing weight. We introduce a new nonlocal hypothesis on the negative part of the solution outside the domain, which together with the negative part of the potential, influences the formation of dead cores and cannot be removed. Our approach relies on uniform bounds from below of the maximum of nontrivial solutions through Liouville theorems, and on a Hopf lemma for viscosity solutions driven by fully nonlinear operators, which we also prove.

34. 基于神经网络从数据中学习偏微分方程的功能性组成部分

Learning functional components of PDEs from data using neural networks

作者: Torkel E. Loman, Yurij Salmaniw, Antonio Leon Villares, Jose A. Carrillo, Ruth E. Baker

类别: cs.LG, math.AP

发布日期: 2026-02-13

arXiv ID: 2602.13174v1

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

摘要: Partial differential equations often contain unknown functions that are difficult or impossible to measure directly, hampering our ability to derive predictions from the model. Workflows for recovering scalar PDE parameters from data are well studied: here we show how similar workflows can be used to recover functions from data. Specifically, we embed neural networks into the PDE and show how, as they are trained on data, they can approximate unknown functions with arbitrary accuracy. Using nonlocal aggregation-diffusion equations as a case study, we recover interaction kernels and external potentials from steady state data. Specifically, we investigate how a wide range of factors, such as the number of available solutions, their properties, sampling density, and measurement noise, affect our ability to successfully recover functions. Our approach is advantageous because it can utilise standard parameter-fitting workflows, and in that the trained PDE can be treated as a normal PDE for purposes such as generating system predictions.

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

1. 可压缩流体-颗粒流动的小德博拉数极限#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 技术性难题与应对技巧#

与之前工作的比较#

a) 突破性改进与优势体现#

2. 带加性噪声的随机Navier–Stokes方程的高阶离散格式#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 处理非线性与奇异性技巧#

与之前工作的比较#

a) 突破与改进#

3. 具有诺伊曼边界条件的二维外区域上二次拟线性波动方程的长时间光滑解#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 关键技术性突破#

与之前工作的比较#

a) 突破性改进#

4. 一维磁流体动力学模型的稳定性与不稳定性#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 处理技术性难题的技巧#

与之前工作的比较#

a) 突破与改进#

5. 四维空间中轴对称无旋欧拉方程解的整体正则性#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 处理非线性与奇性的主要技巧#

与之前工作的比较#

a) 突破与改进#

6. 二维与三维欧拉方程的精确非定常解#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 技术性难题的处理技巧#

与之前工作的比较#

a) 突破与改进#

7. 黏性-色散激波的收缩：零黏度-毛细作用极限#

详细分析#

研究对象和背景#

主要定理或主要结果#

定理 1.1（收缩性定理）#

定理 1.2（全局存在性）#

定理 1.3（零极限的唯一性与稳定性）#

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

a) 核心分析工具#

b) 技术难点与处理技巧#

与之前工作的比较#

a) 突破与改进#

8. Shinbrot型判据与可压缩Navier-Stokes方程的能量守恒#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 关键技术：弱型交换子估计#

与之前工作的比较#

a) 突破性改进#

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

1. 从Cucker–Smale模型大尺度极限导出的动理学方程研究#

摘要翻译#

2. 紧致李群上伪微分算子的核估计与弱 (1,1)-有界性#

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

1. 可压缩流体-颗粒流动的小德博拉数极限

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 技术性难题与应对技巧

与之前工作的比较

a) 突破性改进与优势体现

2. 带加性噪声的随机Navier–Stokes方程的高阶离散格式

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 处理非线性与奇异性技巧

与之前工作的比较

a) 突破与改进

3. 具有诺伊曼边界条件的二维外区域上二次拟线性波动方程的长时间光滑解

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 关键技术性突破

与之前工作的比较

a) 突破性改进

4. 一维磁流体动力学模型的稳定性与不稳定性

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 处理技术性难题的技巧

与之前工作的比较

a) 突破与改进

5. 四维空间中轴对称无旋欧拉方程解的整体正则性

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 处理非线性与奇性的主要技巧

与之前工作的比较

a) 突破与改进

6. 二维与三维欧拉方程的精确非定常解

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 技术性难题的处理技巧

与之前工作的比较

a) 突破与改进

7. 黏性-色散激波的收缩：零黏度-毛细作用极限

详细分析

研究对象和背景

主要定理或主要结果

定理 1.1（收缩性定理）

定理 1.2（全局存在性）

定理 1.3（零极限的唯一性与稳定性）

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

a) 核心分析工具

b) 技术难点与处理技巧

与之前工作的比较

a) 突破与改进

8. Shinbrot型判据与可压缩Navier-Stokes方程的能量守恒

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 关键技术：弱型交换子估计

与之前工作的比较

a) 突破性改进

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

1. 从Cucker–Smale模型大尺度极限导出的动理学方程研究

摘要翻译

2. 紧致李群上伪微分算子的核估计与弱 (1,1)-有界性