生成时间: 2026年02月11日 04:54:25

重点关注论文数量: 3

了解领域论文数量: 4

不相关论文数量: 16

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

1. 强迫Navier-Stokes方程的定量估计及其应用

Quantitative estimates for the forced Navier-Stokes equations and applications

作者: Tobias Barker, Henry Popkin

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09951v1

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

详细分析

中文标题: 强迫Navier-Stokes方程的定量估计及其应用

研究对象和背景

本文研究的核心对象是三维不可压 Navier-Stokes 方程（NSE）及其耦合系统——Boussinesq 方程（Bous.），二者均定义在全空间 $\mathbb{R}^3 \times (0,\infty)$ 上。这些方程描述了粘性、不可压缩流体的动力学行为，是流体力学中最基本且最具挑战性的非线性偏微分方程系统之一。

论文在引言第1节明确给出了所研究的方程组：

$$ \text{(NSE)}\quad \begin{cases} \partial_t v - \Delta v + (v \cdot \nabla) v + \nabla p = 0, \\ \operatorname{div}(v) = 0, \end{cases} $$

其中 $v = v(x,t) \in \mathbb{R}^3$ 是速度场，$p = p(x,t) \in \mathbb{R}$ 是压力场，$\Delta$ 是拉普拉斯算子，$(v \cdot \nabla)v$ 表示非线性对流项。该方程具有著名的尺度不变性：若 $(v,p)$ 是解，则对任意 $\lambda > 0$，函数

$$ v_\lambda(x,t) := \lambda v(\lambda x, \lambda^2 t), \quad p_\lambda(x,t) := \lambda^2 p(\lambda x, \lambda^2 t) $$


也是解，对应初始数据 $v_{0,\lambda}(x) = \lambda v_0(\lambda x)$ 和强迫项 $F_\lambda(x,t) = \lambda^3 F(\lambda x, \lambda^2 t)$。这一性质决定了所谓“临界空间”——即范数在该缩放下保持不变的空间，如 $L^3(\mathbb{R}^3)$、$\dot{H}^{\frac{1}{2}}(\mathbb{R}^3)$、$BMO^{-1}(\mathbb{R}^3)$ 等——在正则性理论中起决定性作用。

与之密切相关的是 Boussinesq 方程组，它刻画了受温度（或密度）扰动影响的热对流现象：

$$ \text{(Bous.)}\quad \begin{cases} \partial_t v - \Delta v + v \cdot \nabla v + \nabla p = \theta e_3, \\ \partial_t \theta - \Delta \theta + v \cdot \nabla \theta = 0, \\ \operatorname{div}(v) = 0, \end{cases} $$

其中 $\theta = \theta(x,t) \in \mathbb{R}$ 表示温度（或密度）场，$e_3 = (0,0,1)^T$ 是标准基向量。注意到第一式右侧出现的 $\theta e_3$ 可视为由温度驱动的外力项，因此 Boussinesq 系统可被看作一类典型的强迫 Navier-Stokes 方程，其强迫项 $F = \theta e_3$ 具有特定结构与正则性。

论文特别关注的是Leray–Hopf 弱解与适合弱解（suitable weak solution）两类解的概念。根据 Definition 1.3.1 与 1.3.2，解需满足分布意义下的方程、能量不等式（global/local energy inequality），以及适当的可积性条件：

  • $v \in L^\infty_{\mathrm{loc}}((0,T); L^2(\mathbb{R}^3)) \cap L^2_{\mathrm{loc}}((0,T); \dot{H}^1(\mathbb{R}^3))$，
  • $p \in L^{3/2}_{\mathrm{loc}}(\mathbb{R}^3 \times (0,T))$，
  • 强迫项 $F$ 属于 $L^2_t H^1_x \cap L^6_{t,x}$ 类空间（见第3页与第7页），这是由 Bogovskii 截断方法引入的局部化强迫所自然导出的正则性。

本文聚焦于奇点形成机制（singularity formation）与爆破速率估计（blow-up rate）。经典结果 ESS03 指出：若 Leray–Hopf 解在有限时间 $T^* > 0$ 处首次失去正则性（即存在 $(x_*, T^*)$ 使得 $v \notin L^\infty(Q_r(x_*,T^*))$ 对所有小 $r>0$ 成立），则必有

$$ \limsup_{t \uparrow T^*} \|v(\cdot,t)\|_{L^3(\mathbb{R}^3)} = \infty. $$


但该结论仅为定性断言；Tao 在 [Tao21] 中首次实现了定量强化，证明存在绝对常数 $c > 0$，使得

$$ \limsup_{t \uparrow T^*} \frac{\|v(\cdot,t)\|_{L^3(\mathbb{R}^3)}}{\left( \log\log\log \frac{1}{T^* - t} \right)^c} = \infty. $$


该结果依赖于Carleman 不等式替代传统的反向唯一性论证，并建立关键的定量正则性估计：若 $\|v\|_{L^\infty_t L^3_x} \le M$，则

$$ \|v(\cdot,t)\|_{L^\infty_x} \le t^{-1/2} \exp\exp\exp(M^c). $$

本文将上述框架推广至含强迫项的情形，尤其针对由局部化过程产生的低正则强迫 $F$。这不仅要求重新构建适用于强迫系统的 Carleman 型估计，还需处理因强迫放大效应带来的大尺度发散问题。此外，区域设定始终为 $\mathbb{R}^3$（无边界），初值属于 $L^2 \cap L^\infty$ 或更一般的 Sobolev 类，解空间强调局部可积性与能量控制，而非全局光滑性。

主要定理或主要结果

本文确立了两个核心定理，分别对应 Navier-Stokes 与 Boussinesq 系统的局部化超临界爆破准则与定量 $L^3$ 爆破速率。

Theorem A（第3页，公式 (6)）：
存在普适常数 $\theta \in (0,1)$，使得如下成立：设 $v: \mathbb{R}^3 \times (0,\infty) \to \mathbb{R}^3$ 是 $\mathbb{R}^3 \times (0,\infty)$ 上的适合 Leray–Hopf 弱解，且在 $T^* > 0$ 处首次爆破，$(x_*, T^*)$ 为其一奇异点。则对任意 $\delta > 0$，有

$$ \limsup_{t \uparrow T^*} \int_{B_\delta(x_*)} \frac{|v(x,t)|^3}{\left( \log\log\log\left( e^{e^{3e}} + |v(x,t)| \right) \right)^\theta} \, dx = \infty. $$


该定理是 [BP21b] 中全局 Orlicz 型爆破准则的局部化版本，将积分区域从全空间 $\mathbb{R}^3$ 缩小至任意小邻域 $B_\delta(x_*)$，从而更贴近数值模拟与物理观测的实际需求。

Theorem B（第4页，公式 (11)）：
设 $(v,\theta,p)$ 是 Boussinesq 方程在 $\mathbb{R}^3 \times (0,T^*)$ 上属于能量类的光滑解，初值满足 $v_0 \in L^2 \cap L^\infty$, $\theta_0 \in L^2 \cap L^\infty$，且 $v$ 在 $T^*$ 处爆破。则存在普适常数（隐含于指数中），使得

$$ \limsup_{t \uparrow T^*} \frac{\|v(\cdot,t)\|_{L^3(\mathbb{R}^3)}}{\left( \log\log\log \frac{1}{(T^* - t)^{1/4}} \right)^{1/614}} = \infty. $$


此即对 Boussinesq 系统中速度场临界范数的首个定量爆破下界，其形式与 Tao 的 NSE 结果平行，但指数经精细调整以适应耦合结构。

此外，文中还建立了支撑上述定理的关键定量正则性估计（Section 3），即针对强迫 NSE 的如下结论（隐含于 Proposition 2.6.1–2.6.2 及后续推导中）：
若 $v$ 满足 $\|v\|_{L^\infty_t L^3_x(Q_1)} \le M$ 且强迫项 $F$ 满足 $\|F\|_{L^6_{t,x}(Q_1)} \le \varepsilon$（$\varepsilon$ 充分小），则在适当选取的“正则环形区域”上，可控制 $v$ 的高阶导数与涡量 $\omega = \nabla \times v$，并导出类 Gaussian 下界（见公式 (15)）。

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

a) 核心分析学工具

本文的技术骨架建立在三大支柱之上：Carleman 不等式、Caccioppoli 型局部能量估计 与 Serrin 型迭代引导（bootstrapping）。

  • Carleman 不等式 是全文最核心的工具，用于替代传统反向唯一性与唯一延拓定理的定量版本。具体而言，作者使用两类 Carleman 权重函数：

    • 第一类对应唯一延拓性（Proposition 2.6.1）：对满足 $|\partial_t \omega - \Delta \omega| \le C^{-1}|\omega| + C^{-1/2}|\nabla \omega| + |\nabla \times F|$ 的涡量 $\omega$，有
      $$ \int_{-T_1}^0 \int_{B_R \setminus B_{R/2}} |\omega|^2 \, dx ds \gtrsim T_1^{-1/2} e^{-C M T_1 R^2}. $$
      此即公式 (15) 的严格表述，它提供了关键的 Gaussian 下界。
    • 第二类对应反向唯一性（Proposition 2.6.2）：在环形区域上控制解在 $t=0$ 时刻的大小，结合第一类不等式完成闭环。

  • Caccioppoli 不等式（Proposition 2.1.4）用于处理低正则强迫下的梯度控制。由于 $F$ 仅属 $L^6_{t,x}$，无法直接控制 $\nabla^2 v$，故需借助如下形式的局部估计：

    $$ \int_I \int_{B_{(2)}} |\nabla W|^2 \, dx ds \lesssim \frac{1}{\operatorname{dist}(\partial B_{(1)}, \partial B_{(2)})^2} \int_I \int_{B_{(1)}} |W|^2 \, dx ds + \|G\|_{L^2}^2 + \|H\|_{L^2} \|W\|_{L^2}, $$


    其中 $W$ 为 Stokes 系统解，$G,H$ 含 $F$ 及其导数。该不等式避免了对高阶导数的直接依赖，是克服强迫低正则性的关键技术。

  • Serrin 型 bootstrapping（Lemma 2.1.1–2.1.2）提供从 $L^3$ 控制到 $L^\infty$ 的逐级提升。例如，由 $\|v\|_{L^\infty_t L^3_x} \le M$ 与 $\|F\|_{L^6_{t,x}} \le \varepsilon$，先得 $\|\nabla v\|_{L^6_{t,x}} \lesssim \varepsilon + M^2$，再代入涡量方程 $\partial_t \omega - \Delta \omega = \operatorname{div}(\omega \otimes v - v \otimes \omega) + \nabla \times F$，利用 Calderón–Zygmund 理论获得 $\|\omega\|_{C^{0,1/6}_p}$ 控制，最终导出 $\|\nabla v\|_{L^\infty}$ 估计。

b) 技术性难题的解决技巧

首要障碍是：Carleman 不等式中强迫项在大尺度 $R \gg 1$ 下被指数放大，导致标准推导失效。作者提出三重策略予以克服：

  1. 尺度截断与“正则环形区域”构造（Proposition 2.5.1）：
    利用强迫项的次临界性（subcriticality），即 $\|F\|_{L^6_{t,x}}$ 在缩放下按 $\lambda^{1/2}$ 衰减，可选取尺度 $R \in [2M T_1^{1/2}, M T_1^{1/2} \exp(M^{607})]$，使得 Carleman 不等式适用的环形区域

    $$ \{ 4\hat{R} < |x| < M^{1200} \hat{R} \} \times (-T_1/2, 0) $$


    存在且可控。此即公式 (17)，它显式限定了 $R$ 的上界，从而将强迫项放大因子 $T e^{C M' T R^2}$ 控制在可接受范围。

  2. 强迫项的小性补偿：
    通过缩放论证（rescaling argument），可假设 $\|F\|_{L^6_{t,x}} \le \varepsilon$ 极小。结合前述 $R$ 的上界，迫使放大后的强迫项

    $$ T e^{C M' T R^2} \| \nabla \times F \|_{L^2(B_R \times (-T,0))} $$


    仍远小于主项 $T_1^{-1/2} e^{-C M T_1 R^2}$，从而可忽略。

  3. Caccioppoli 与低正则估计耦合：
    针对 $F \in L^6_{t,x}$ 导致的 $\nabla^2 v$ 缺失，作者在 Carleman 推导中不直接估计 $\nabla^2 v$，而是用 Caccioppoli 不等式将 $\|\nabla v\|_{L^2}$ 控制在 $\|v\|_{L^2} + \|F\|_{L^2}$ 量级，并利用插值与 Hölder 不等式将高阶项降为可处理形式，例如

    $$ \int |\nabla v| |\nabla \phi| |v| \, dx \lesssim \|\nabla v\|_{L^2} \|\nabla \phi\|_{L^6} \|v\|_{L^3}. $$

与之前工作的比较

a) 突破与改进

本文在多个维度实现显著推进：

  • 首次建立强迫 NSE 的定量正则性估计：此前 [Tao21]、[BP21a]、[Bar23] 均假设无强迫或强迫紧支且远离奇点区域；而本文处理的是由局部化自然诱导的、全局定义但低正则的强迫 $F \in L^2_t H^1_x \cap L^6_{t,x}$，且必须在整个正则区域内一致控制。这是对定量方法适用范围的根本性拓展。

  • Orlicz 准则的局部化：[BP21b] 证明了全局积分型爆破

    $$ \int_{\mathbb{R}^3} \frac{|v|^3}{\left( \log\log\log(e^{e^{3e}} + |v|) \right)^\theta} \, dx \to \infty, $$


    但该量在数值计算中难以捕捉局域奇点。本文 Theorem A 将其收缩至任意小球 $B_\delta(x_*)$，使准则具备可验证性与物理可解释性，呼应 Hou 等关于数值奇点检测的讨论 [Hou23]。

  • Boussinesq 系统的定量爆破：虽已有大量关于 Boussinesq 正则性准则的工作（如 [CV07], [ZL08]），但均属定性或亚临界范畴。Theorem B 首次给出 $L^3$ 范数的定量发散速率，且其形式与 NSE 完全平行，揭示了热对流与纯流体运动在奇点结构上的深层共性。

  • 技术层面的精细化：相较于 [Bar23] 中“强迫远离关注区域”的简化设定，本文直面强迫与奇点共存的困难情形，通过尺度依赖的环形区域设计、强迫次临界性的主动利用及Caccioppoli 与 Carleman 的协同嵌套，为未来处理更复杂耦合系统（如 MHD、nematic liquid crystals）提供了可复用的方法论模板。

2. 具有Navier滑移-摩擦边界条件的外部区域中Oseen流的大时间衰减

Large time decay of the Oseen flow in exterior domains subject to the Navier slip-with-friction boundary condition

作者: Toshiaki Hishida

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09804v1

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

详细分析

中文标题: 具有Navier滑移-摩擦边界条件的外部区域中Oseen流的大时间衰减

研究对象和背景

本文研究的是三维无界外部区域 $\Omega \subset \mathbb{R}^3$ 中，满足 Navier滑移-摩擦边界条件 的 Oseen方程 的大时间行为。该问题源于不可压缩黏性流体（Navier-Stokes流）在障碍物外部的运动建模，其中障碍物占据有界闭集 $K = \mathbb{R}^3 \setminus \Omega$，$\partial\Omega$ 是光滑边界（属于 $C^{2,1}$ 类），且 $\Omega$ 为连通无界区域。

物理背景是：当流体在无穷远处以恒定向量 $u_\infty \in \mathbb{R}^3$ 定常流出时，若将参考系平移至随 $u_\infty$ 运动，则原Navier-Stokes系统可线性化为Oseen系统，其未知量为扰动速度场 $u(x,t)$ 和扰动压力 $p(x,t)$，满足：

$$ \begin{cases} \partial_t u = \Delta u + \eta \cdot \nabla u - \nabla p, & \text{div } u = 0 \quad \text{in } \Omega \times (0,\infty), \\ u(\cdot,0) = f, & \\ \nu \cdot u|_{\partial\Omega} = 0, & \\ [(2Du)\nu]_\tau + \alpha u_\tau|_{\partial\Omega} = 0, & \\ \lim_{|x|\to\infty} u(x,t) = 0, \end{cases} \tag{1.4--1.5} $$

其中：

  • $\eta = u_\infty$ 是常向量，代表远场恒定流速；
  • $Du = \frac{1}{2}(\nabla u + (\nabla u)^\top)$ 是对称梯度（变形张量）；
  • $\nu(x)$ 是 $\partial\Omega$ 上的外法向单位向量；
  • $W_\tau = W - (\nu \cdot W)\nu$ 表示任意向量场 $W$ 的切向分量；
  • $\alpha = \alpha(x) \geq 0$ 是定义在 $\partial\Omega$ 上的连续函数，称为摩擦系数，控制滑移程度；
  • 边界条件 $(1.5)_2$ 即 Navier slip-with-friction condition，等价于经典应力条件：
    $$ [T(u,p)\nu]_\tau + \alpha u_\tau = 0, \quad \text{其中 } T(u,p) = 2Du - pI, $$
    但值得注意的是，该条件不显含压力 $p$，仅通过切向应力与切向速度耦合。

该边界条件具有重要物理意义：当 $\alpha \equiv 0$ 时为完全滑移（full-slip）；当 $\alpha \to \infty$ 时形式上恢复无滑移（no-slip） 条件 $u|_{\partial\Omega}=0$；而中间值 $\alpha > 0$ 则对应实际中粗糙表面或微结构诱导的有限滑移长度 $1/\alpha$。与传统无滑移条件相比，Navier条件更适于描述刚体碰撞流体系统中的有限时间接触行为——此类现象在无滑移假设下被证明不可能发生（见文献[17]），而在Navier条件下已被证实可行（见[15]）。

函数空间设定亦具特色：解取值于solenoidal Lebesgue空间

$$ L^q_\sigma(\Omega) := \overline{C^\infty_{0,\sigma}(\Omega)}^{\|\cdot\|_{L^q(\Omega)}}, $$


其中 $C^\infty_{0,\sigma}(\Omega)$ 表示支集紧、散度为零的光滑向量场全体；该空间可由经典的Helmholtz分解刻画：

$$ L^q(\Omega) = L^q_\sigma(\Omega) \oplus \{ \nabla p \in L^q(\Omega); \, p \in L^q_{\mathrm{loc}}(\Omega) \}, $$


并存在Fujita–Kato投影 $P_q : L^q(\Omega) \to L^q_\sigma(\Omega)$。相应地，Stokes算子 $A$ 与 Oseen算子 $L$ 定义为：

$$ D(A) = \{ u \in L^q_\sigma(\Omega) \cap W^{2,q}(\Omega); \; [(2Du)\nu]_\tau + \alpha u_\tau = 0 \text{ on } \partial\Omega \}, \quad Au = -P_q \Delta u, $$

$$ D(L) = D(A), \quad Lu = -P_q[\Delta u + \eta \cdot \nabla u]. $$

关键难点在于：外部区域 $\Omega$ 的非紧性导致谱理论复杂化；边界条件含位置依赖的摩擦系数 $\alpha(x)$，破坏均匀性；且需处理各向异性衰减结构（如wake效应）；此外，远场非零（即 $\eta \neq 0$）使得基本解不再具有卷积平移不变性，迫使必须精细分析复参数 $\lambda$ 邻域（尤其是虚轴附近）的预解式正则性。

主要定理或主要结果

论文确立了两个核心定理，分别关于预解集结构与大时间衰减估计：

Theorem 2.1（预解集包含性）
设 $\eta \in \mathbb{R}^3 \setminus \{0\}$，$\alpha \in C(\partial\Omega)$ 非负，且满足

$$ \alpha(x) \geq \frac{\eta \cdot \nu(x)}{2}, \quad \forall x \in \partial\Omega. \tag{2.5} $$


则对任意 $q \in (1,\infty)$，有

$$ \mathbb{C} \setminus \{0\} \subset \rho(-L), $$


即除原点外，整个复平面属于 $-L$ 的预解集；进一步，若还满足几何加强条件

$$ \alpha(x) + \kappa_0(x) \geq \frac{\eta \cdot \nu(x)}{2}, \quad \forall x \in \partial\Omega, \tag{2.7} $$


其中 $\kappa_0(x) \leq 0$ 是边界 $\partial\Omega$ 在 $x$ 处沿 $\nu(x)$ 方向的主曲率最小值（Weingarten映射 $-\nabla N$ 的最小特征值），则

$$ \mathbb{C} \setminus S_\eta \subset \rho(-L), \quad \text{其中 } S_\eta := \left\{ \lambda \in \mathbb{C};\; |\eta|^2 \operatorname{Re}\lambda + (\operatorname{Im}\lambda)^2 \leq 0 \right\}. \tag{2.8} $$

Theorem 2.2（$L^q$–$L^r$ 大时间衰减）
在相同假设 (2.5) 下，对任意 $1 < q \leq r \leq \infty$（当 $j=0$）或 $1 < q \leq r < \infty$（当 $j=1$），存在常数 $C = C(m,\alpha,q,r,\Omega) > 0$，使得对所有 $t > 0$、$|\eta| \leq m$ 及 $f \in L^q_\sigma(\Omega)$，成立

$$ \| \nabla^j e^{-tL} f \|_r \leq C t^{-\frac{j}{2} - \frac{3}{2}\left( \frac{1}{q} - \frac{1}{r} \right)} (1 + t)^{\max\left\{1 - \frac{3}{2r}, 0\right\}} \|f\|_q, \quad j = 0,1. \tag{2.9} $$


特别地，当 $\eta = 0$（Stokes情形）时，该估计退化为已知最优结果；当 $\alpha \to \infty$ 时，形式上恢复无滑移情形的经典衰减率。

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

a) 核心分析学工具

  1. 预解式分析与谱理论：全文围绕算子 $-L$ 的预解式 $(\lambda + L)^{-1}$ 展开。关键在于证明其在 $\lambda = 0$ 附近（特别是虚轴 $\operatorname{Re}\lambda = 0$）的一致有界性与Hölder连续性。为此，作者采用局部能量估计+截断技巧（cut-off procedure）：将全空间分解为内域 $B_R$ 与外域 $\Omega \setminus B_R$，分别构造内解与外解，并通过全局拼接获得整体估计。

  2. Weingarten映射与曲率补偿技术：为处理边界项中出现的形变能 $\|Du\|_2^2$ 与曲率耦合项，在内域估计中引入引理3.1，建立切向导数与曲率算子的关系：

    $$ [(\nabla u)^\top \nu]_\tau = -(\nabla N) u \quad \text{on } \partial\Omega, $$

    进而导出关键不等式：

    $$ 2\|Du\|_{2,D}^2 \geq \|\nabla u\|_{2,D}^2 + \int_{\partial D} \kappa_0 |u|^2 d\sigma. $$

    此式将几何信息 $\kappa_0(x)$ 显式嵌入能量估计，成为验证唯一性的枢纽。

  3. Fourier乘子理论与全空间预解式正则性：利用全空间Oseen预解式的显式Fourier表示：

    $$ u(\cdot,\lambda) = \mathcal{F}^{-1} \left[ \frac{|\xi|^2 I - \xi \otimes \xi}{(\lambda + |\xi|^2 - i \eta \cdot \xi)|\xi|^2} \widehat{f}(\xi) \right], $$

    结合Mikhlin-Hörmander乘子定理，得到对 $\lambda \in \mathbb{C} \setminus S_\eta$ 的 $W^{2,q}$ 估计（式(3.7)–(3.9)），并进一步建立沿虚轴的一致有界性（式(3.10)）与1/2-Hölder连续性（式(3.12)）：

    $$ \sup_{\tau \in \mathbb{R}} \|E_\eta(i\tau)\|_{\mathcal{L}(L^q_{[R]}, W^{2,q}(B_3))} \leq C, $$

    $$ \|E_\eta(i(\tau+h)) - E_\eta(i\tau)\| \leq C |h|^{1/2}. $$

  4. 解析半群理论与扰动论证：基于Stokes算子 $-A$ 生成解析半群的已知结论（Shibata–Shimada, 2003），将 $-L = -A - P_q(\eta \cdot \nabla \cdot)$ 视为小摄动，利用Neumann级数论证得到预解式在扇形区域 $\Sigma_{\varepsilon,\lambda_0}$ 内的存在性与估计（式(2.4)）。

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

  • 非齐次边界与无穷远条件的齐次化：通过变量替换 $u \mapsto u - u_\infty$ 将原始Navier-Stokes问题转化为以 $0$ 为远场的Oseen问题，同时将非齐次边界条件吸收进初始数据与源项处理中。

  • 压力项的消去与投影技巧：借助Fujita–Kato投影 $P_q$，将压力梯度项 $\nabla p$ 自动剔除出方程，使分析聚焦于 $L^q_\sigma$ 空间上的纯速度演化，避免处理压力的奇异性。

  • 局部能量衰减的构造：核心创新在于将“大时间衰减”问题转化为“预解式在 $\lambda = 0$ 附近的行为”问题。通过证明 $(\lambda + L)^{-1}$ 在 $\{\operatorname{Re}\lambda \geq 0\} \setminus \{0\}$ 上的一致有界性，并结合复积分表示（如Bromwich积分）与Phragmén–Lindelöf型极大值原理，最终导出指数型衰减的积分核控制，再经插值获得 $L^q$–$L^r$ 型代数衰减。

  • 对摩擦系数 $\alpha(x)$ 的非均匀处理：未要求 $\alpha$ 为常数，而是允许其为边界上连续函数，并通过逐点不等式 (2.5)、(2.7) 将其与局部几何量（法向分量、曲率）耦合，从而实现对任意形状障碍物的普适性。

与之前工作的比较

a) 突破与改进

本文最根本的推进在于首次在Navier边界条件下建立了Oseen半群的大时间 $L^q$–$L^r$ 衰减估计，填补了此前文献的空白。以往工作可分为三类：

  • 无滑移情形（$\alpha \to \infty$）：Kobayashi–Shibata (2005) 已完整解决，证明 $\mathbb{C} \setminus S_\eta \subset \rho(-L)$ 并导出 (2.9) 型估计。但其方法严重依赖边界刚性约束 $u|_{\partial\Omega}=0$，无法推广至滑移情形。

  • 完全滑移情形（$\alpha \equiv 0$）：Shimada–Yamaguchi (2009) 仅对 Stokes 算子（$\eta = 0$）获得衰减，且明确指出：一旦 $\eta \neq 0$，Oseen 算子在完全滑移下不具备大时间衰减性——因缺乏足够耗散机制抑制切向运动。本文通过引入正摩擦 $\alpha(x) > 0$，成功激活耗散通道，使 $\eta \neq 0$ 时衰减成为可能。

  • 一般 $\alpha > 0$ 的先验研究：虽有大量关于Navier边界的Stokes/N-S适定性工作（如[1,6,11,16,28–29,36,38–39]），但全部假设远场为零（$u_\infty = 0$），即未触及Oseen框架。本文首次将 $\alpha(x)$ 与 $\eta$、$\nu(x)$、$\kappa_0(x)$ 三者置于统一不等式 (2.5)/(2.7) 下，揭示了滑移强度、来流方向、边界凸性之间的精妙平衡关系：只有当局部摩擦足以抵消来流对边界的“剥离效应”（由 $\eta \cdot \nu/2$ 刻画）时，系统才保持稳定并呈现代数衰减。

具体而言，改进体现在：

  • 假设条件显著放宽：相比无滑移情形需 $\eta \neq 0$ 即可，本文要求 $\alpha(x) \geq \eta \cdot \nu(x)/2$，允许 $\alpha$ 在背风面（$\eta \cdot \nu < 0$）自然为零，更符合物理直觉；
  • 几何适应性增强：引入曲率修正项 $\kappa_0(x)$，使结果适用于任意 $C^{2,1}$ 边界，而非仅限于凸域或小扰动球体；
  • 结果强度提升：所得衰减率 (2.9) 与无滑移情形完全一致，表明Navier条件在适当摩擦下可完全复现无滑移的长期动力学特性，这是此前任何工作均未证实的。

综上，本文不仅完成了从Stokes到Oseen、从无滑移到Navier的双重跨越，更提炼出支配外部流动稳定性的统一几何-力学判据，为流体力学边界条件设计提供了严格的数学基准。

3. 从斜向波强迫到条带强化：一种基于扰动的频率响应框架

From oblique-wave forcing to streak reinforcement: A perturbation-based frequency-response framework

作者: Dušan Božić, Anubhav Dwivedi, Mihailo R. Jovanović

类别: physics.flu-dyn, math.AP, math.DS, math.OC

发布日期: 2026-02-09

arXiv ID: 2602.09137v1

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

详细分析

中文标题: 从斜向波强迫到条带强化：一种基于扰动的频率响应框架

研究对象和背景

本文研究的是不可压缩牛顿流体在壁面约束剪切流中的转捩动力学，核心系统为Navier–Stokes (NS) 方程组线性化于层流基流 $ \bar{u}(y) = (U(y),0,0) $ 附近的扰动演化。具体考虑的是压力驱动的平行平板通道流（Poiseuille 流），其无量纲化的层流速度剖面为

$$ U(y) = 1 - y^2, $$


定义在对称区间 $ y \in [-1,1] $ 上，对应上下壁面位于 $ y = \pm 1 $。所有物理量已按通道半宽 $ \delta $、最大层流速度 $ \bar{U} $、惯性时间尺度 $ \delta/\bar{U} $ 和压力尺度 $ \rho \bar{U}^2 $ 进行无量纲化，从而引入雷诺数

$$ \mathrm{Re} = \frac{\bar{U} \delta}{\nu}, $$


其中 $ \nu $ 为运动粘性系数。

设速度与压力扰动分别为 $ u(x,y,z,t) = (u,v,w) $ 和 $ p(x,y,z,t) $，则其满足如下非定常不可压缩 Navier–Stokes 方程（见原文式 (1)）：

$$ \partial_t u = \Delta u - (\bar{u} \cdot \nabla) u - (u \cdot \nabla) \bar{u} - (u \cdot \nabla) u - \nabla p + \frac{1}{\mathrm{Re}} d, $$


$$ \nabla \cdot u = 0. $$


此处 $ d = (d_u,d_v,d_w) $ 是外加扰动源项；$ \Delta = \nabla \cdot \nabla $ 是拉普拉斯算子；$ \nabla $ 表示梯度；$ \partial_t $ 为时间偏导。该方程组定义在三维空间域 $ \Omega = \mathbb{R}_x \times [-1,1]_y \times \mathbb{R}_z $ 上，并配备无滑移与无穿透边界条件：

$$ u(x,y=\pm 1,z,t) = 0. $$

为便于输入–输出分析，作者采用标准的壁面法向速度–涡量表示 $ q = (v,\eta) $，其中 $ \eta = \partial_z u - \partial_x w $，将 NS 方程转化为状态空间形式（原文式 (4a)）：

$$ \partial_t q^{(n)} = A q^{(n)} + B d^{(n)}, $$


$$ u^{(n)} = C q^{(n)}. $$


这里 $ A $ 是 Orr–Sommerfeld/Squire 算子，即 NS 方程关于层流基流 $ \bar{u} $ 的线性化算子；$ B $ 和 $ C $ 分别刻画输入作用方式与状态到速度场的映射。边界条件为

$$ v^{(n)}(x,y=\pm 1,z,t) = \partial_y v^{(n)}(x,y=\pm 1,z,t) = \eta^{(n)}(x,y=\pm 1,z,t) = 0. $$

进一步地，在流向恒定（streamwise-constant）假设下（即 $ \partial_x (\cdot) = 0 $），系统可降维至 $ (y,z) $ 平面，引入流函数 $ \psi $ 满足 $ (v,w) = (\partial_z \psi, -\partial_y \psi) $，并导出如下耦合演化方程（原文式 (5)）：

$$ \partial_t \psi^{(n)} = \Delta^{-1} \Delta^2 \psi^{(n)} + \frac{1}{\mathrm{Re}} d_\psi^{(n)}, $$


$$ \partial_t u^{(n)} = -U'(y) \partial_z \psi^{(n)} + \Delta u^{(n)} + \frac{1}{\mathrm{Re}} d_u^{(n)}. $$


其中 $ \Delta = \partial_{yy} + \partial_{zz} $ 是二维拉普拉斯，$ \Delta^2 = \partial_{yyyy} + 2 \partial_{yyzz} + \partial_{zzzz} $ 是双调和算子；强迫项 $ d_\psi^{(n)} $ 由式 (6) 给出：

$$ d_\psi^{(n)} = \Delta^{-1} \left( \partial_z d_v^{(n)} - \partial_y d_w^{(n)} \right). $$


该模型专用于刻画流向条带（streamwise streaks） 的生成与演化，是全文分析的核心简化框架。

研究背景聚焦于亚临界转捩（subcritical transition） —— 即在经典模态不稳定性（如 Tollmien–Schlichting 波）尚未出现的雷诺数范围内，流动仍可通过非模态放大机制自发失稳。实验与直接数值模拟（DNS）反复证实：斜向波（oblique waves） —— 具有非零流向与展向波数 $ (k_x,k_z) \neq (0,k_z) $ 的三维非定常扰动 —— 是触发此类转捩最有效的路径之一。它们通过二次非线性相互作用激发稳态条带，后者再经二级不稳定性（secondary instability） 破碎为湍流。然而，传统线性理论（如 resolvent 分析）仅能描述无限小扰动的瞬态增长，无法定量刻画有限幅值下的条带形成、相位调控及转捩阈值；而完全非线性模拟计算成本高昂且物理解释性弱。因此，本文旨在构建一个兼具机理透明性与计算高效性的统一框架，弥合线性非模态分析与经典模态不稳定性理论之间的鸿沟。

主要定理或主要结果

本文未以传统“定理–证明”形式陈述结论，而是通过系统性扰动分析建立若干严格成立的结构化命题，其核心成果可归纳为以下五点：

  1. 二阶条带生成定理：对于幅值为 $ \epsilon $ 的斜向波强迫 $ d^{(1)} $，其二次非线性耦合 $ d^{(2)} = - (u^{(1)} \cdot \nabla) u^{(1)} $ 在 $ O(\epsilon^2) $ 阶产生稳态流向条带 $ u_s^{(2)} $，其空间结构精确由流向恒定 resolvent 算子 $ H_{(0,k_z)}(\omega=0) $ 的第二输出奇异函数 $ \vartheta_{(0,k_z),2}(y) $ 所刻画。即存在常数 $ \alpha $ 使得

    $$ u_s^{(2)}(y,z) = \alpha \, \vartheta_{(0,2k_z),2}(y) \cos(2k_z z). $$

  2. 高阶条带强化/衰减判据：在 $ O(\epsilon^n) $ 阶（$ n \geq 3 $），斜向波与已生成条带的非线性耦合产生新条带分量 $ u_s^{(n)} $，其空间结构仍由同一 $ \vartheta_{(0,2k_z),2}(y) $ 主导，但相位取决于低阶扰动间的相对相位：若相位对齐，则发生条带强化（streak reinforcement）；若反相，则导致衰减（attenuation）。数学上体现为叠加系数的符号一致性。

  3. 弱非线性失效阈值定理：存在一临界强迫幅值 $ \epsilon_{\mathrm{cr}} $，当 $ \epsilon > \epsilon_{\mathrm{cr}} $ 时，扰动级数发散，标志弱非线性近似失效。该阈值与 DNS 中首次观测到持续非定常性的临界点精确吻合。

  4. 非模态–模态统一性定理：上述 $ \epsilon_{\mathrm{cr}} $ 同时对应于条带畸变基流的二级不稳定性起始点。即，由扰动分析自洽导出的条带修正 $ \bar{u} + \epsilon^2 u_s^{(2)} + \cdots $，其线性稳定性分析所给出的最不稳定模态增长率在 $ \epsilon = \epsilon_{\mathrm{cr}} $ 处由负变正。这建立了 $ \textbf{resolvent-driven non-modal amplification} $ 与 $ \textbf{classical modal instability} $ 的直接定量联系，无需预设条带幅值。

  5. Shanks 变换加速收敛定理：对原始发散的扰动级数应用 Shanks 变换后，所得加速序列在 $ \epsilon \lesssim \epsilon_{\mathrm{cr}} $ 范围内保持良好收敛性，且能定量复现 DNS 中条带能量演化与平均流修正，验证了框架的预测能力。

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

a) 核心分析学工具

本文的数学骨架是多尺度摄动理论（multi-scale perturbation theory） 与频响算子谱分析（frequency-response operator spectral analysis） 的深度融合。其严谨性根植于以下几类工具：

  • Hilbert 空间上的线性算子理论：将 NS 方程扰动系统建模为作用于 $ L^2([-1,1]) $ 空间的无穷维线性算子 $ A $、$ B $、$ C $ 的组合。所有 resolvent 算子 $ H_k(\omega) = C(i\omega I - A_k)^{-1}B_k $ 均被严格定义在适当的 Sobolev 空间（如 $ H^2_0 $）上，确保其紧性与 Fredholm 性质，从而保证奇异值分解（SVD）的存在性与完备性。

  • Fourier–Laplace 变换与谱参数化：利用流向 $ x $ 与展向 $ z $ 的平移不变性，对扰动进行 Fourier 变换，将偏微分方程组解耦为一族参数化常微分方程（ODE）系统，每个系统由波数 $ k = (k_x,k_z) $ 和频率 $ \omega $ 标记。此步骤将问题转化为对一族 $ y $-方向边值问题的分析，极大降低了维度。

  • 奇异值分解（SVD）与最优模态识别：对每个 $ H_k(\omega) $ 实施 SVD，得到输入/输出正交基 $ \{ \delta_{k,j} \}, \{ \vartheta_{k,j} \} $ 及奇异值 $ \sigma_{k,j}(\omega) $。主奇异值 $ \sigma_{k,1}(\omega) $ 刻画最大可能增益，而高阶模态（如 $ j=2 $）则编码特定物理机制（如 lift-up）的主导结构。

  • 解析延拓与 Galilean 变换技巧：为处理非定常斜向波，引入移动坐标系 $ \bar{x} = x - ct $，将时间依赖项 $ e^{i(k_x x + \omega t)} $ 转化为稳态形式 $ e^{ik_x \bar{x}} $，此时频响算子变为 $ G_k(c) = H_k(-c k_x) $。该变换将非定常响应问题转化为稳态边值问题，使数值求解可行，并保留全部物理信息。

b) 技术性难题的解决技巧

  • 非线性项的系统截断与层次化建模：将非线性对流项 $ -(u \cdot \nabla) u $ 按 $ \epsilon $ 幂次展开，得到递推关系式 (3b)：

    $$ d^{(n)} = - \sum_{r=1}^{n-1} (u^{(r)} \cdot \nabla) u^{(n-r)}. $$


    此式将高阶动力学完全由低阶响应驱动，形成严格耦合的线性系统链（见图 2）。该结构避免了传统弱非线性理论中需预设基流修正的循环论证，实现了“从层流基流出发，自洽生成条带”的逻辑闭环。

  • 流向恒定近似的精准定位与验证：虽假设 $ \partial_x u = 0 $ 以简化条带动力学，但作者通过完整三维扰动展开（式 (18)）明确指出：二阶条带 $ u_s^{(2)} $ 的主导分量确为 $ (k_x,k_z) = (0,2k_z) $，且其能量远超其他 $ O(\epsilon^2) $ 分量（如 $ (2k_x,0) $ 的二维模态）。这从数学上证成该简化的合理性。

  • 发散级数的数值解析延拓：针对 $ \epsilon \to \epsilon_{\mathrm{cr}} $ 时摄动级数失效的问题，采用 Shanks 变换 这一非线性序列加速技术。对部分和 $ S_n = \sum_{j=1}^n \epsilon^j u^{(j)} $ 构造变换

    $$ S_n^{\mathrm{Shanks}} = \frac{S_{n-1} S_{n+1} - S_n^2}{S_{n-1} - 2S_n + S_{n+1}}, $$


    有效压制高频振荡，显著拓宽了有效预报范围，使理论与 DNS 在临界区域实现定量一致。

  • 高精度数值实现：采用 ultraspherical spectral method 与 spectral integration scheme 离散 $ y $-方向，结合 Chebyshev 配置点（71 点），确保 Orr–Sommerfeld/Squire 算子的高保真离散，使 $ H_\infty $ 范数与奇异函数的计算误差控制在机器精度量级。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重范式跃迁：

第一，超越线性 resolvent 分析的定性局限。传统 resolvent 理论（如 McKeon & Sharma, JFM 2010）虽能识别 $ \vartheta_{(0,k_z),1} $ 作为最优条带结构，却无法解释其起源（为何是第二模态而非第一？）、幅值标度律（为何是 $ \epsilon^2 $？）及相位敏感性（为何某些斜向波组合强化条带而另一些削弱？）。本文通过二阶摄动，首次严格证明：斜向波的二次耦合天然激发 $ \vartheta_{(0,2k_z),2} $，其物理根源在于 lift-up 机制中涡量–速度的奇偶性匹配，从而将“最优结构”从经验观察升华为可推导的必然结果。

第二，革新二级稳定性分析的本构基础。经典二级稳定性（Herbert, AGARD 1988；Andersson et al., JFM 2001）需人为指定条带幅值，构造畸变基流 $ \bar{u} + U_s(y)\cos(k_z z) $，再对其线性化。此做法引入主观性，且无法回答“条带幅值多大时二级不稳定性启动”。本文则完全摒弃预设，从第一性原理出发，由 $ \epsilon $ 展开自洽导出 $ U_s(y) \sim \epsilon^2 \vartheta_{(0,2k_z),2}(y) $，并将二级不稳定性阈值 $ \epsilon_{\mathrm{cr}} $ 定义为扰动级数的收敛半径，实现了从“假设条带”到“生成并预测条带”的根本转变。

第三，建立非模态与模态理论的定量桥梁。以往工作（如 Schmid & Henningson, Stability and Transition in Shear Flows）将二者视为割裂范式：前者描述瞬态，后者描述渐近。本文则证明：非模态放大是模态不稳定的前置引擎——斜向波通过非线性生成条带，条带自身又通过修改平均剪切率诱发模态不稳定性；且两者的临界参数（$ \epsilon_{\mathrm{cr}} $）严格重合。这终结了长期存在的“非模态 vs 模态”二元争论，提出统一的“非线性频响框架”。

综上，本文不仅在技术细节上（如 Shanks 加速、ultraspherical 离散）优于前人，更在哲学层面完成了从“现象分类”到“机制溯源”、从“模型拟合”到“原理推演”的深刻跃进。

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

1. $N$维重力驱动的Muskat问题

The $N$-dimensional gravity driven Muskat problem

作者: Bogdan-Vasile Matioc, Georg Prokert

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09786v1

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

摘要翻译

中文标题: $N$维重力驱动的Muskat问题

摘要翻译: 本文研究Muskat问题，该问题描述了在占据全空间 $\mathbb{R}^{N+1}$（其中 $N \geq 2$）的均匀多孔介质中，由重力驱动的两种互不混溶、不可压缩流体的运动。两流体之间的界面表示为定义在 $\mathbb{R}^N$ 上的一个函数的图像。我们将该问题重新表述为关于该函数的非线性、非局部演化方程，其中出现的奇异积分源于速度场与压力场的位势表示。借助调和分析的相关结果，我们证明：在Rayleigh–Taylor条件所确定的开集上，该演化具有抛物型性质。进一步，我们运用此类问题的抽象理论，证明Muskat问题在此开集上定义了一个半流（semiflow），且该半流作用于所有次临界Sobolev空间 $H^s(\mathbb{R}^N)$（其中 $s > s_c$），此处临界指数为 $s_c = 1 + N/2$。此外，我们还获得了直至 $\mathrm{C}^\infty$ 类的抛物正则化效应（即解在任意正时间后自动属于 $\mathrm{C}^\infty$）。

2. 具有时滞的一侧型与抛物型BLO空间及其在Muckenhoupt $A_1$ 权和双重非线性抛物方程中的应用

One-Sided and Parabolic BLO Spaces with Time Lag and Their Applications to Muckenhoupt $A_1$ Weights and Doubly Nonlinear Parabolic Equations

作者: Weiyi Kong, Dachun Yang, Wen Yuan

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

发布日期: 2026-02-10

arXiv ID: 2602.09741v1

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

摘要翻译

中文标题: 具有时滞的一侧型与抛物型BLO空间及其在Muckenhoupt $A_1$ 权和双重非线性抛物方程中的应用

摘要翻译: 本文首先引入一侧型BLO空间 $\mathrm{BLO}^+(\mathbb{R})$，并分别借助一侧型Muckenhoupt权类 $A_1^+(\mathbb{R})$ 及一侧型John–Nirenberg不等式对该空间进行刻画。以此为基础，我们建立了 $\mathrm{BLO}^+(\mathbb{R})$ 函数的Coifman–Rochberg型分解，并证明该空间的定义不依赖于所涉两个区间之间的距离；这一性质进一步导出了 $\mathrm{BLO}^+(\mathbb{R})$ 关于一侧型BMO空间 $\mathrm{BMO}^+(\mathbb{R})$ 的刻画（即Bennett型引理）。作为应用，我们证明了任一 $\mathrm{BMO}^+(\mathbb{R})$ 函数均可分解为两个 $\mathrm{BLO}^+(\mathbb{R})$ 函数之和，并给出了 $\mathrm{BLO}^+(\mathbb{R})$ 中函数到 $L^\infty(\mathbb{R})$ 的距离的显式描述。最后，作为高维类比，我们引入具时滞的抛物型BLO空间 $\mathrm{PBLO}_\gamma^-(\mathbb{R}^{n+1})$，并将上述所有一维结果系统地推广至该空间；进一步地，在应用层面，我们不仅建立了 $\mathrm{PBLO}_\gamma^-(\mathbb{R}^{n+1})$ 与双重非线性抛物方程解之间的联系，还就“抛物距离函数的负对数属于 $\mathrm{PBLO}_\gamma^-(\mathbb{R}^{n+1})$”这一性质，以集合的弱多孔性（weak porosity）为条件，给出了一个必要性判据。

3. 静电效应对黏性可压缩流体临界正则性及长时间行为的影响

Electrostatic effects on critical regularity and long-time behavior of viscous compressible fluids

作者: Ling-yun Shou, Zihao Song

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09332v1

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

摘要翻译

中文标题: 静电效应对黏性可压缩流体临界正则性及长时间行为的影响

摘要翻译:
本文研究 $\mathbb{R}^d$（$d \geq 2$）中可压缩Navier–Stokes–Poisson方程组，该方程组是描述等熵可压缩流体与自洽静电势耦合的经典模型。我们指出：由于静电耦合所蕴含的底层Klein–Gordon结构，其对解的长时间动力学行为具有显著影响。作为首要结果，我们在全频段 $L^p$ 型临界Besov空间中证明了该柯西问题在平衡态邻域内初值下的整体适定性——且无需借助双曲对称化技巧。相较于若干里程碑式工作（如Charve与Danchin, Arch. Rational Mech. Anal., 198 (2010), 233–271；Chen、Miao与Zhang, Commun. Pure Appl. Math., 63 (2010), 1173–1224；Haspot, Arch. Rational Mech. Anal., 202 (2011), 427–460）中对无Poisson项情形的研究，我们消除了低频段额外的 $L^2$ 可积性假设，并将容许的 $p$ 值范围拓展至最优区间 $1 \leq p < 2d$。据我们所知，这是可压缩流体力学中首个允许初始速度场在所有频率尺度上均呈现高度振荡性的整体适定性结果。

进一步地，受Poisson耦合驱动，密度与速度在低频段展现出本质不同的渐近行为。受此特征启发，我们提出一类通用的 $L^p$ 型低频假设，并在此框架下建立了全局解趋近平衡态的最优收敛速率。对于一大类指标参数，该假设所导出的衰减速率快于经典 $L^1$ 框架下所得结果。为此，我们发展了一种时间加权能量方法——该方法本身具有独立意义，且使我们得以在不附加初始数据小性条件的前提下，精确捕获最大可能的衰减估计。

4. 二维Vlasov–Poisson方程的离子–电子态

Ions-electrons-states for the two-component Vlasov-Poisson equation

作者: Emeric Roulley

类别: math.AP, physics.plasm-ph

发布日期: 2026-02-10

arXiv ID: 2602.09293v1

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

摘要翻译

中文标题: 二维Vlasov–Poisson方程的离子–电子态

摘要翻译:
本文建立了关于一维双组分Vlasov–Poisson方程行波周期解的局部与全局分歧结果。此类解在相空间中表现为条带状的离子与电子分布区域，它们以相干方式传播，并从空间均匀、但速度依赖的平衡层中产生。分歧图的拓扑结构取决于底层平衡态的几何性质以及所选取的傅里叶模态：在不同情形下，分歧图呈现两条或四条解分支。在所有情形中，该分歧均为叉形（pitchfork）分歧；而在具有对称性的构型中，平衡点附近的局部结构具有双曲几何特征。我们进一步证明：这些局部构造的解分支均可延拓为全局分支。本工作推广了作者此前关于纯电子情形的研究——在该研究中，离子被建模为静止的中和背景。而当前允许离子与电子两类粒子均作动力学演化，导致分析维度升高、结构更为复杂。最后，通过一个仿射变量变换，我们揭示了该模型与一维双组分Euler–Poisson系统之间的内在联系；这一联系进而使得我们亦能为此Euler–Poisson模型构造出任意小振幅乃至大振幅的行波周期解。

不相关论文（基本信息）

1. 一种力学-化学耦合的图式形成模型的稳定性与分岔分析

Stability and bifurcation analysis in a mechanochemical model of pattern formation

作者: Szymon Cygan, Anna Marciniak-Czochra, Finn Münnich, Dietmar Oelz

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09998v1

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

摘要: We analyze the stability and bifurcation structure of steady states in a mechanochemical model of pattern formation in regenerating tissue spheroids. The model couples morphogen dynamics with tissue mechanics via a positive feedback loop: mechanical stretching enhances morphogen production, while morphogen concentration modulates tissue elasticity. Global strain conservation implements a nonlocal inhibitory effect, realizing a mechanochemical variant of the local activation–long-range inhibition mechanism. For exponential elasticity-morphogen coupling, the system admits a variational formulation. We prove existence of nonconstant steady states for small diffusion and uniqueness of the homogeneous state for large diffusion. Linear stability analysis shows that only unimodal patterns are stable, while multimodal solutions are unstable. Bifurcation analysis reveals subcritical and supercritical pitchforks, with fold bifurcations generating bistable regimes. Our results demonstrate that mechanochemical feedback provides a robust mechanism for single-peaked pattern formation without requiring a second diffusible inhibitor.

2. 动态规划原理的粘性解框架及其应用

A Viscosity Framework for Dynamic Programming Principles and Applications

作者: Félix del Teso, Julio D. Rossi, Jorge Ruiz-Cases

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09946v1

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

摘要: In this work we introduce a viscosity-based notion of solution for general approximation schemes associated with partial differential equations, such as dynamic programming principles~(DPPs). A key feature of our approach is that it bypasses any measurability requirement on solutions of the DPP, an assumption that is often difficult to verify and may even fail in relevant examples. We establish a comparison principle between classical strict supersolutions and viscosity subsolutions of the DPP, which yields stability results under minimal and natural hypotheses. As a consequence, we prove existence of viscosity solutions of the DPP and their convergence to viscosity solutions of a PDE that is consistent with the underlying approximation scheme. Moreover, we show that solutions of the limiting PDE admit an asymptotic expansion encoded by the approximation operator. Finally, we demonstrate that a broad class of local, nonlocal, and nonlinear partial differential equations fits into our framework, recovering known examples in the literature and completing gaps in the existing literature.

3. 混合区域中双非线性方程的正则性

Regularity for Doubly Nonlinear Equations in the Mixed Regime

作者: Simone Ciani, Eurica Henriques, Mariia Savchenko, Igor I. Skrypnik, Yevgeniia Yevgenieva

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09906v1

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

摘要: We study the local Hölder continuity of nonnegative solutions to doubly nonlinear equations by introducing a new technique that allows us to treat the cases where the equation is both singular and degenerate, up to specific Barenblatt numbers. Our argument relies on a new integral $L^1$-$L^1$ Harnack estimate, of independent interest.

4. 基于Kuttler–Sigillito型的微分形式几何特征值估计

Geometric eigenvalue estimates of Kuttler-Sigillito type on differential forms

作者: Rodolphe Abou Assali

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

发布日期: 2026-02-10

arXiv ID: 2602.09876v1

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

摘要: We introduce a new biharmonic Steklov problem on differential forms with Dirichlet-type boundary conditions and show that it is elliptic. We prove the existence of a discrete spectrum for this problem and give variational characterizations for eigenvalues associated to it. We establish eigenvalue estimates known as Kuttler-Sigillito inequalities, that connect the eigenvalues of different problems on differential forms with curvature quantities on the manifold.

5. 多组分动力学模型：广义热力学动力学演化（GENERIC）框架与费舍尔信息

Multi-species kinetic models: GENERIC formulation and Fisher information

作者: Manh Hong Duong, Zihui He

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09875v1

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

摘要: In this paper, we study the GENERIC structures of multi-species spatially inhomogeneous Boltzmann and Landau equations with Bose-Einstein, Maxwell-Boltzmann, and Fermi-Dirac statistics. In addition, under suitable assumptions on the collision kernels, we show that the Fisher information for the multi-species spatially homogeneous Boltzmann equation is non-increasing in time.

6. 一种通过测度空间中的凸化与对偶性求解双重非线性方程的直接方法

A direct method for doubly nonlinear equations via convexification in spaces of measures and duality

作者: Alessandro Pinzi, Filippo Riva, Giuseppe Savaré

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09808v1

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

摘要: Existence of solutions to doubly nonlinear equations in reflexive Banach spaces is established by resorting to a global-in-time variational approach inspired by De Giorgi’s principle, which characterizes the associated flows as null-minimizers of a suitable energy-dissipation functional defined on trajectories. In contrast to the celebrated minimizing movements scheme, the proposed strategy does not rely on any time-discretization or iterative constructions. Instead, it provides a direct method based on the relaxation of the problem in spaces of measures, constrained by the continuity equation: in this procedure, no gap is introduced due to the Ambrosio’s superposition principle. Within this weak convex framework, the validity of the null-minimization property is recovered through two further steps. First, a careful application of the Von Neumann minimax theorem yields an identification of the dual problem as a supremum over the set of smooth and bounded cylinder functions, solving an Hamilton-Jacobi-type inequality. Secondly, a suitable “backward boundedness” property of solutions to such Hamilton-Jacobi system gives a proper bound of the dual problem, ensuring that the minimum value of the original functional is actually zero. The proposed strategy naturally extends to non-autonomous equations, encompassing time- and space-dependent dissipation potentials and time-dependent potential energies.

7. 单稳态反应-扩散方程中解向推动型行波解的收敛性及水平集的行为

Convergence to pushed fronts and the behavior of level sets in monostable reaction-diffusion equations

作者: Ryo Kiyono

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09806v1

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

摘要: We study the behavior of solutions of a monostable reaction-diffusion equation $u_t=Δ_x u +u_{yy} +f(u)$ ($x \in \mathbb{R}^{n-1}$, $y \in \mathbb{R}$, $t>0$), with the unstable equilibrium point $0$ and the stable equilibrium point $1$. Under the condition that the corresponding one-dimensional equation has a pushed front $Φ_{c^*}(z)$ with $Φ_{c^*}(-\infty)=1$, $Φ_{c^*}(\infty)=0$, we show that the solution $u(x,y,t)$ approaches $Φ_{c^*}(y-γ(x,t))$ for some $γ(x,t)$ as $t \to \infty$, if initially $u(x,y,0)$ decays sufficiently fast as $y \to \infty$ and is bounded below by some positive constant near $y=-\infty$. It is also shown that $γ(x,t)$ is approximated by the mean curvature flow with a drift term.

8. 亚极端克尔时空上玻色型特科尔斯基因子场的昂鲁态

The Unruh state for bosonic Teukolsky fields on subextreme Kerr spacetimes

作者: Dietrich Häfner, Christiane K. M. Klein

类别: math-ph, gr-qc, math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09796v1

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

摘要: We perform the quantization of Teukolsky scalars of spin $0$, $\pm 1$, and $\pm 2$ within the algebraic approach to quantum field theory. We first discuss the classical phase space, from which we subsequently construct the algebra. This sheds light on which fields are conjugates of each other. Further, we construct the Unruh state for this theory on Kerr and show that it is Hadamard on the black hole exterior and the interior up to the inner horizon. This shows not only that Hadamard states exist for this theory, but also extends the existence and Hadamard property of the Unruh state to (bosonic) Teukolsky fields on Kerr, where such a result was previously missing.

9. 具有拉皮尼-帕普拉尔表面锚定效应的相变

Phase Transition With Rapini-Papoular Surface Anchoring

作者: Shun Li, Yong Yu

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09778v1

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

摘要: We analyze the dynamical (in)stability of nematic liquid crystals in the presence of external magnetic fields and Rapini-Papoular surface potential. The P-HAN transition is investigated using a simplified 3D Ericksen-Leslie system. We find the thickness threshold of the P-HAN transition. If the thickness of the nematic layer exceeds this threshold, there is a global-in-time suitable weak solution converging exponentially to a nontrivial equilibrium state as time tends to infinity. If the thickness is no more than the threshold, the global-in-time suitable weak solution has a trivial long-time asymptotic limit. Our results rigorously justify the P-HAN transition discussed in the physics literature.

10. 非单调零折扣问题中基于静态类别引导的基本解选择

Static class-guided selection of elementary solutions in non-monotone vanishing discount problems

作者: Panrui Ni, Jun Yan, Maxime Zavidovique

类别: math.AP, math.DS

发布日期: 2026-02-10

arXiv ID: 2602.09697v1

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

摘要: We study a generalized vanishing discount problem for Hamilton–Jacobi equations, removing the standard monotonicity assumption, either in a global sense or when integrated against all Mather measures. Specifically, we consider

$$ λa(x)u(x)+H(x,Du(x))-Aλ=c_0, $$

with a suitably chosen constant $A>0$. By appropriately changing the signs of the function $a(x)$ on different static classes associated with $H$, we show that the maximal viscosity solution converges uniformly as $λ\to 0^+$ and that all elementary solutions of the stationary equation

$$ H(x,Du(x))=c_0 $$

can be selected as limits. This provides the first result for selecting multiple viscosity solutions in vanishing discount problems beyond the usual monotonicity and integral assumptions, as long as $a(x)$ is positive on one static class. Our results highlight the crucial role of static classes in controlling the asymptotic behavior of viscosity solutions. Previously, under usual monotonicity assumptions, only a single solution could be selected (as discussed in \cite{GL}), whereas our approach allows controlled selection of multiple solutions via static class-guided discount coefficients.

11. 下有界Ricci曲率的非紧黎曼流形上的精确Sobolev不等式与Moser–Trudinger不等式

Sharp Sobolev and Moser-Trudinger inequalities on noncompact Riemannian manifolds with Ricci curvature bounded below

作者: Carlo Morpurgo, Liuyu Qin

类别: math.AP, math.DG

发布日期: 2026-02-10

arXiv ID: 2602.09677v1

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

摘要: We establish Sobolev and Moser-Trudinger inequalities with best constants on noncompact Riemannan manifolds with Ricci curvature bounded below, and positive injectivity radius.

12. 某类分数阶调和函数再生核希尔伯特空间上的 Lions 型公式

A Lions’ type formula for some reproducing kernel Hilbert spaces of fractional harmonic functions

作者: Sidy M. Djitte, Franck Sueur

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09669v1

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

摘要: In \cite{Lions}, J. L. Lions considered a reproducing kernel Hilbert space (RKHS) of harmonic functions on a regular domain with Sobolev traces and obtained a formula that expresses the kernel of this space as an integral on the boundary of some derivatives of the Green function associated with the Laplace operator and the homogeneous Dirichlet boundary condition. This result was simplified and extended later by Englis, Lukkassen, Peetre, and Persson in \cite{ELPL} to more general elliptic systems of even orders. In particular, they emphasized that the resemblance between Lions’ type formula and the Hadamard variational formula only appears when the operator is of order $2$. In this paper, we investigate some RKHS of $a$-harmonic functions, where $a$ in $(0,1)$ refers to a fractional exponent of the Laplace operator. For such fractional order pseudo-differential operators, the local nonhomogeneous Dirichlet problem can be addressed by means of some $a$-transmission Sobolev spaces, which were introduced by Hörmander in the sixties and recently developed by Grubb in a series of papers. We deduce from these works a fractional Poisson formula, which is applied to obtain a Lions’ type formula. We observe, in particular, that despite the order of the operator not being $2$, this formula resembles the Hadamard variational formula that we prove in the companion paper \cite{sidy-franck_1}. As a complementary remark, we observe that for a family of RKHS associated with the steady Stokes system, a second order system, there is also a Lions’ type formula for their two-point kernels, which turn out not to be similar to the corresponding Hadamard variation formula.

13. 时间周期性盖尔范德–希洛夫空间中的微分复形

Differential Complexes in Time-Periodic Gelfand-Shilov Spaces

作者: Fernando de Ávila Silva, Marco Cappiello, Alexandre Kirilov, Pedro Meyer Tokoro

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09646v1

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

摘要: We study the global solvability of a class of differential complexes on the product manifold $\mathbb{T}^m \times \mathbb{R}^n$ associated with systems of evolution operators of the form $L_r = \partial_{t_r} + ia_r(t)P(x,D_x), r=1,\ldots,m,$ where the coefficients $a_r$ are real-valued Gevrey functions on the torus and $P(x,D_x)$ is a globally elliptic normal differential operator on $\mathbb{R}^n$. Within the framework of time-periodic Gelfand–Shilov spaces, we introduce a natural differential complex generated by these operators and investigate its solvability in both functional and ultradistributional settings. We provide a complete characterization of global solvability in terms of a Diophantine condition involving the constant part of the associated $1$-form and the spectrum of $P$. We also analyze global hypoellipticity of the complex. These results extend previous works on scalar operators and constant coefficient systems to the setting of differential complexes with time-dependent real coefficients.

14. 随机演化介质中非局部方程的均匀化：扩散近似

Homogenization of nonlocal equations in randomly evolving media. Diffusion approximation

作者: Marina Kleptsyna, Andrey Piatnitski, Alexandre Popier

类别: math.AP, math.PR

发布日期: 2026-02-10

arXiv ID: 2602.09584v1

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

摘要: The paper deals with homogenization and higher order approximations of solutions to nonlocal evolution equations of convolution type whose coefficients are periodic in the spatial variables and random stationary in time. We assume that the convolution kernel has finite moments up to order three. Under proper mixing assumptions, we study the limit behavior of the normalized difference between solutions of the original and the homogenized problems and show that this difference converges to the solution of a linear stochastic partial differential equation.

15. 时周期合作系统的主谱理论与时间非局部扩散下的渐近分析

Principal spectral theory and asymptotic analysis for time-periodic cooperative systems with temporally nonlocal dispersal

作者: Hao Wu, Wan-Tong Li, Jian-Wen Sun, Hoang-Hung Vo

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09436v1

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

摘要: This paper investigates the principal spectral theory and the asymptotic behavior of the principal spectrum point for a class of time-periodic cooperative systems with nonlocal dispersal operators, incorporating both coupled and uncoupled nonlocal terms. By applying the theory of resolvent positive operators and their perturbations, we first establish criteria for the existence of the principal eigenvalue. We then construct sequences of smooth upper and lower approximating matrix-valued functions, each of whose corresponding operators satisfies the principal eigenvalue existence condition. This approximation framework allows the principal spectrum point to effectively substitute for the principal eigenvalue in characterizing the global dynamics of the nonlinear system. Moreover, it facilitates the study of the asymptotic behavior of the principal spectrum point with respect to parameters under fairly general assumptions. Subsequently, for systems with both coupled and uncoupled nonlocal terms, we analyze the asymptotic behavior of the principal spectrum point in terms of the dispersal rate, dispersal range, and frequency. Finally, we illustrate the applicability of our theoretical results through a Zika virus model and a stem cell model.

16. Baouendi 算子的强唯一延拓性结果

A strong unique continuation result for the Baouendi operator

作者: Agnid Banerjee, Nicola Garofalo

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.09322v1

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

摘要: We establish a strong unique continuation property for the subelliptic Baouendi operator under the presence of zero-order perturbations satisfying an almost Hardy-type growth condition. In particular, the admissible class includes both $L^\infty_{\mathrm{loc}}$ and singular potentials. We prove that any solution vanishing to infinite order at a point of the degeneracy manifold of the operator must be identically zero. The result holds extends to variable-coefficient operators with intrinsic Lipschitz regularity. A notable feature of the proof is that it relies exclusively on $L^2$ Carleman estimates combined with the classical Hardy inequality.