Home » 论文总结

2026年02月10日论文分析

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

February 10, 2026 · 68 min · 14262 words · P0lyno3ial

Table of Contents

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

生成时间: 2026年02月10日 05:00:22

重点关注论文数量: 12

了解领域论文数量: 4

不相关论文数量: 30

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

1. 小 $BV \cap L^1$ 初值下动态燃烧中一维可压 Navier–Stokes 方程组的整体适定性

Global well-posedness for one-dimensional compressible Navier–Stokes system in dynamic combustion with small $BV\cap L^1$ initial data

作者: Siran Li, Haitao Wang, Jianing Yang

类别: math.AP, math-ph

发布日期: 2026-02-09

arXiv ID: 2602.08867v1

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

详细分析

中文标题: 小 $BV \cap L^1$ 初值下动态燃烧中一维可压 Navier–Stokes 方程组的整体适定性

1. 研究对象和背景

本文研究的是一类描述动态燃烧过程中一维可压缩反应气体混合物演化的偏微分方程组——即 Lagrangian 坐标下的 1D 可压 Navier–Stokes–Combustion 系统。该系统耦合了流体动力学（质量、动量守恒）与热力学/化学反应过程（能量守恒、组分输运），其物理变量包括：比容 $v = 1/\rho$ （密度 $\rho$ 的倒数）、速度 $$u$$ 、温度 $\theta$ 和反应物质量分数 $$z$$ 。

系统在 Lagrangian 坐标 $(x,t) \in \mathbb{R} \times [0,\infty)$ 下写为：

\begin{cases} v_t - u_x = 0, \\ u_t + p_x = \left(\mu \dfrac{u_x}{v}\right)_x, \\ \theta_t + p u_x - \mu \dfrac{u_x^2}{v} = \left(\nu \dfrac{\theta_x}{v}\right)_x + q K \phi(\theta) z, \\ z_t + K \phi(\theta) z = \left(D \dfrac{z_x}{v^2}\right)_x. \end{cases} \tag{0.4}

其中，压力 $$p$$ 和内能 $$e$$ 满足理想气体状态方程：

p = \dfrac{a \theta}{v}, \quad e = c_v \theta,

且绝热指数 $\gamma = c_p / c_v > 1$ ；函数 $\phi(\theta) \ge 0$ 是 Lipschitz 连续的反应速率函数，刻画 Arrhenius 型化学反应强度；常数 $\mu > 0$ , $\nu > 0$ , $$D > 0$$ , $$K > 0$$ , $q \in \mathbb{R}$ 分别表示粘性系数、热传导系数、扩散系数、反应速率系数、反应热差。

该模型源于燃烧理论中对预混火焰传播、爆轰波结构及不稳定性分析的建模需求，具有典型的多物理场强耦合、非线性、部分耗散（仅 $v,u,\theta$ 含二阶耗散项，而 $$z$$ 仅具扩散项） 特征。特别地， $$z$$ -方程不含时间导数的显式耗散机制（如 $\theta_{xx}$ ），导致其正则性更弱，是技术难点之一。

论文考虑的是全实轴 $\mathbb{R}$ 上的 Cauchy 问题，无边界条件；初始数据 $(v_0,u_0,\theta_0,z_0)$ 被假定为围绕恒定平衡态 $(1,0,1,0)^\top$ 的小扰动，且属于交空间：

(v_0 - 1,\, u_0,\, \theta_0 - 1,\, z_0) \in (L^1 \cap BV)(\mathbb{R})^4.

这里 $BV(\mathbb{R})$ 表示有界变差函数空间，其范数定义为：

\|w\|_{BV} = \|w\|_{L^\infty} + \mathrm{Tot.Var.}(w;\mathbb{R}),

其中 $\mathrm{Tot.Var.}(w;\mathbb{R}) = \sup \sum_i |w(x_{i+1}) - w(x_i)|$ 对所有分划取上确界。该空间允许初值含激波型间断（jump discontinuities），远弱于 Sobolev 空间 $$H^1$$ 或 $W^{1,1}$ ，但强于纯 $$L^1$$ 。因此，本文工作处于弱解理论与双曲-抛物耦合系统的低正则性分析前沿。

值得注意的是，解的空间正则性框架（Definition 0.1）并非传统 Sobolev 类型，而是基于 $$BV$$ 时间演化性质与加权 $L^\infty$ 估计的混合范数：

$v(\cdot,t)-1 \in C([0,t_\sharp); L^1 \cap BV)$ ，
$u,\theta,z$ 及其空间导数满足 $\sqrt{t}\, u_x \in L^\infty$ , $t\, u_t \in L^\infty$ 等带时间权重的尖锐估计，反映热核主导的扩散衰减行为。

此设定直接呼应 Liu–Yu（2022）与 Wang–Yu–Zhang（2022）对等熵/全 Navier–Stokes 系统的 $$BV$$ 理论，但首次将之拓展至含化学反应源项与组分输运的四维系统，显著增加了线性化算子的结构复杂性与 Green 函数分析难度。

2. 主要定理或主要结果

论文建立了如下三个核心定理，构成完整的局部存在性、唯一性、稳定性与整体适定性理论：

定理 0.2（局部存在性与正则性）

存在普适常数 $\delta > 0$ ，若初值满足小性条件：

\|v_0 - 1\|_{L^1 \cap BV} + \|u_0\|_{L^1 \cap BV} + \|\theta_0 - 1\|_{L^1 \cap BV} + \|z_0\|_{L^1 \cap BV} < \delta,

则存在 $t_\sharp > 0$ ，使得系统 (0.4) 在 $\mathbb{R} \times [0,t_\sharp)$ 上存在弱解 $(v,u,\theta,z)$ ，满足 Definition 0.1 中的所有正则性要求，并有如下统一控制：

\max\Big\{ \|u(\cdot,t)\|_{L^1 \cap L^\infty},\; \sqrt{t}\,\|u_x(\cdot,t)\|_{L^\infty},\; t\,\|u_t(\cdot,t)\|_{L^\infty},\; \dots \Big\} \le 2C\delta,

对 $v-1,\theta-1,z$ 类似。此外，通量函数：

\frac{\mu u_x}{v} - p,\quad \frac{\nu \theta_x}{c_v v} - \int_{-\infty}^x \left( \frac{p u_y}{v} - \frac{\mu u_y^2}{c_v v} - \frac{q K \phi(\theta) z}{c_v v} \right) dy,\quad \frac{D z_x}{v^2}

关于 $$x$$ 是全局 Lipschitz 连续的；且 $$v$$ 关于时间满足 Hölder 连续性：

\|v(\cdot,t) - v(\cdot,s)\|_{BV} \le 2C\delta \frac{(t-s)|\log(t-s)|}{\sqrt{t}},\quad \|v(\cdot,t) - v(\cdot,s)\|_{L^1} \le 2C\delta (t-s).

定理 0.3（稳定性与唯一性）

存在更小的 $\delta > 0$ ，对任意两组初值 $(v_0^\epsilon,u_0^\epsilon,\theta_0^\epsilon,z_0^\epsilon)$ 与 $(v_0^\iota,u_0^\iota,\theta_0^\iota,z_0^\iota)$ 满足上述小性条件，则对应弱解 $(v^\epsilon,u^\epsilon,\theta^\epsilon,z^\epsilon)$ 与 $(v^\iota,u^\iota,\theta^\iota,z^\iota)$ 在公共存在时间 $$[0,t_*)$$ 内满足稳定性估计：

\sum_{\alpha \in \{v,u,\theta,z\}} \|\alpha^\epsilon(\cdot,t) - \alpha^\iota(\cdot,t)\|_{L^1} \le C_2 \sum_{\beta} \|\beta_0^\epsilon - \beta_0^\iota\|_{L^1 + L^\infty + BV},

从而保证解的唯一性。

定理 0.4（整体存在性与大时间行为）

存在 $\delta > 0$ ，若初值满足小性条件，则局部解可延拓为全局弱解 $(v,u,\theta,z)$ 定义于 $\mathbb{R} \times [0,\infty)$ ，并具有最优衰减速率：

\sqrt{t+1}\,\Big( \|v(\cdot,t)-1\|_{L^\infty} + \|u(\cdot,t)\|_{L^\infty} + \|\theta(\cdot,t)-1\|_{L^\infty} + \|z(\cdot,t)\|_{L^\infty} \Big) \le C_3 \delta,

\sqrt{t+1}\,\Big( \|u_x(\cdot,t)\|_{L^\infty} + \|\theta_x(\cdot,t)\|_{L^\infty} + \|z_x(\cdot,t)\|_{L^\infty} \Big) \le C_3 \delta.

该 $O(t^{-1/2})$ 衰减率与热方程基本解的衰减一致，被证明是最优的。

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

a) 核心分析工具

论文采用 “迭代构造 + Green 函数精细分析” 双轨策略，深度融合现代 PDE 分析中的多项关键技术：

Banach 不动点定理：用于局部解构造。作者设计一个精心加权的 Banach 空间 $\mathcal{X}_{t_\sharp}$ ，其范数包含 $$BV$$ 、 $$L^1$$ 、加权 $L^\infty$ 及其时间导数项，确保非线性映射 $T: (v,u,\theta,z) \mapsto (\tilde{v},\tilde{u},\tilde{\theta},\tilde{z})$ 是压缩映射（Lemma 2.2, 2.6）。关键在于对 $$z$$ -方程的处理：利用反应项 $K \phi(\theta) z$ 的 Lipschitz 性与扩散项的 $$L^1$$ 散度结构，建立 $\|z\|_{L^1 \cap BV}$ 的闭合先验估计。
BV 函数的分解理论：依据 Radon–Nikodým 定理，将 $$v$$ 分解为绝对连续部分 $$v_a$$ 、跳跃部分 $$v_j$$ 与奇异连续部分 $$v_c$$ （见 (1.2)–(1.3)）。由于 $$v_j$$ 由 Heaviside 函数叠加而成，其演化受特征线与粘性正则化共同支配，需单独估计（Theorem 0.2 中的 jump 控制）。
Green 函数方法：这是本文最核心的工具。线性化系统 (3.2) 是一个 $4 \times 4$ 的部分退化双曲-抛物系统：
$$ \partial_t U + A \partial_x U = \partial_x (B \partial_x U) + F, $$
其中 $U = (v,u,\theta,z)^\top$ ，矩阵 $$A$$ 有三个双曲模与一个慢速 $$z$$ -模 $\lambda_4 = 0$ 。作者构造其 Green 矩阵 $$G(x,t;y,0)$$ ，并分频分析：
- 高频分析 ( $\eta \to \infty$ )：通过 Fourier 变换与谱分析，证明 $$G$$ 的奇异性被 $e^{-c \eta^2 t}$ 控制，对应热扩散主导；
- 低频分析 ( $\eta \to 0$ )：利用谱投影与 Jordan 块分解，提取 $$z$$ -模的慢速传播特性，得到 $G_{4j}$ 的 $O(t^{-1/2})$ 衰减（Lemma 4.1）。
加权能量估计与最大值原理：对 $$z$$ 应用比较原理，结合 $0 \le z_0 \le 1$ 与反应项非负性，得 $0 \le z(x,t) \le 1$ ；对 $v,\theta$ 利用最小模原理控制正性（隐含于 (0.6)）。

b) 技术性难题的克服

$$z$$ -方程的弱耗散性： $$z_t$$ 项无二阶导数，导致标准能量法失效。作者引入修正通量 $\mathcal{F}_z = D z_x / v^2$ ，并证明其 Lipschitz 连续性（Theorem 0.2），从而将 $$z$$ 的 $$BV$$ 控制转化为对 $\mathcal{F}_z$ 的控制。
Green 矩阵的精细估计：由于 $$z$$ 引入新模态， $$G$$ 不再是块对角形式。作者发展了 “模态分离 + 逐元素估计” 技术：对 $G_{4j}$ （ $$j=1,2,3,4$$ ）分别估计，尤其对 $G_{44}$ （ $$z$$ 自身响应）使用低频展开：
$G_{44}(x,t;y,0) = \frac{1}{\sqrt{4\pi D t}} e^{-\frac{(x-y)^2}{4 D t}} + \text{error},$
并严格控制误差项的 $$L^1$$ 与 $$BV$$ 范数（Lemma 3.5, 4.1）。
非线性项的 $$BV$$ 乘积估计：对形如 $p_x = (a \theta / v)_x = a (\theta_x / v - \theta v_x / v^2)$ 的项，需估计 $$v_x / v^2$$ 的 $$BV$$ 范数。利用 $v \approx 1$ 及 $v \in BV$ ，通过链式法则与 Young 不等式得：
$\left\| \frac{v_x}{v^2} \right\|_{L^1} \le C \|v_x\|_{L^1},\quad \mathrm{Tot.Var.}\left( \frac{v_x}{v^2} \right) \le C \|v\|_{BV}.$

4. 与之前工作的比较

a) 突破与改进

本文是对 Liu–Yu（2022）与 Wang–Yu–Zhang（2022）工作的实质性推广，主要突破体现在三方面：

模型维度升级：前作处理的是 2×2 等熵 NS（0.7） 与 3×3 NS–Fourier（0.9） 系统，而本文处理 4×4 NS–Combustion 系统，首次将 $$BV$$ 理论应用于含化学反应源项与组分输运的模型。这不仅增加未知函数个数，更使线性化算子从严格抛物（NS–Fourier）退化为部分双曲-抛物，Green 矩阵出现零速模 $\lambda_4 = 0$ ，其长时间行为需全新分析。
初值正则性要求的保持与适配：尽管系统更复杂，本文仍维持与前作相同的初值空间 $L^1 \cap BV$ ，未额外要求更高正则性（如 $$H^1$$ 或 $W^{1,p}$ ）。这得益于对 $$z$$ 的特殊处理：不依赖其高阶导数，而仅用 $$L^1$$ 与 $$BV$$ 控制其扩散与反应耦合效应。
结果强度的增强：除整体存在性外，本文首次给出 $$z$$ 的最优 $O(t^{-1/2})$ 一致衰减（Theorem 0.4），并证明其通量 $\mathcal{F}_z$ 的 Lipschitz 连续性（Theorem 0.2），这是前作未涉及的新结论。此外，对解的时间 Hölder 连续性（(0.14)）的刻画，比单纯 $L^\infty_t L^\infty_x$ 更精细，揭示了 $$BV$$ 解的时间振荡本质。

综上，本文将 $$BV$$ 方法成功植入反应流体力学领域，为理解燃烧过程中激波、火焰锋面与扩散层的相互作用提供了严格的数学基础。

2. 抽象积分微分方程及其应用

Abstract integrodifferential equations and applications

作者: Bruno de Andrade, Marcos Gabriel de Santana

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08691v1

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

详细分析

中文标题: 抽象积分微分方程及其应用

研究对象和背景

本文研究的核心对象是一类抽象积分微分方程（abstract integrodifferential equation），其一般形式为：

\begin{cases} u'(t) = \displaystyle\int_0^t g(t-s) A u(s)\,ds + f(t,u(t)), & t \geq 0, \\ u(0) = u_0, \end{cases} \tag{1.3}

其中 $A : D(A) \subset X_0 \to X_0$ 是定义在复Banach空间 $$ X_0 $$ 上的闭、无界、稠定线性算子，且 $$ -A $$ 是扇形算子（sectorial operator）；核函数 $g \in L^1_{\mathrm{loc}}([0,\infty);\mathbb{C})$ 是非恒零、Laplace 可积函数；非线性项 $$ f $$ 定义在由 $$ -A $$ 生成的插值尺度 $\{X_\alpha\}_{\alpha \geq 0}$ 上，满足局部 Lipschitz 或更弱的正则性条件。

该方程并非孤立抽象构造，而是源于两类重要物理模型：具遗传粘性的 Navier-Stokes 方程与具记忆效应的反应-扩散方程。

第一类动机是三维有界光滑区域 $\Omega \subset \mathbb{R}^3$ 上的 Navier-Stokes 方程，其粘性项不再采用经典 Newtonian 形式 $-\nu \Delta u$ ，而是引入材料记忆核 $$ g $$ ，即：

\begin{cases} u_t(x,t) = \displaystyle\int_0^t g(t-s) \Delta u(x,s)\,ds - u \cdot \nabla u - \nabla p + h(x,t), & (x,t) \in \Omega \times (0,\infty), \\ \mathrm{div}\, u = 0, & \text{in } \Omega \times (0,\infty), \\ u(x,t) = 0, & (x,t) \in \partial\Omega \times (0,\infty), \\ u(x,0) = u_0(x), & x \in \Omega. \end{cases} \tag{1.1}

此处 $u : \Omega \times (0,\infty) \to \mathbb{R}^3$ 表示流体速度场， $$ p $$ 为压力， $$ h $$ 为外力，而 $$ g $$ 是刻画非牛顿流体或湍流统计粘弹性行为的材料函数。特别地，当取幂律型核 $(1 * g)(t) = t^\alpha / \Gamma(\alpha+1)$ （ $\alpha \in [0,1)$ ）时，对应于分数阶粘性模型，其物理意义在于描述具有时间依赖松弛机制的流变响应。

第二类动机来自Gurtin-Pipkin 型热传导模型，用于描述均匀各向同性刚性热导体中具记忆效应的热流。它将 Fourier 定律 $q = -c \nabla u$ 替换为：

q(x,t) = -\int_0^\infty g(t-s) \nabla u(x,s)\,ds,

从而导出如下演化方程：

u_t(x,t) = \int_0^t g(t-s) \Delta u(x,s)\,ds + f(t,u(x,t)), \quad (x,t) \in \Omega \times (0,\infty). \tag{1.2}

该方程中非线性项 $$ f $$ 可为超线性型（如 $|u|^{p-1}u$ ， $$ p > 1 $$ ）或梯度型（如 $|\nabla u|^q$ ， $$ q > 1 $$ ），分别刻画强非平衡热源或非局部能量耗散机制。

论文特别强调：由于核积分项破坏了半群性质（即 $S(t+s) \neq S(t)S(s)$ ），传统的 $$ C_0 $$ -半群理论无法直接适用；必须借助 Prüss 提出的解族（resolvent family）理论，即一族满足特定卷积型再生公式的算子族 $\{S(t)\}_{t \geq 0}$ ，使得 mild 解可表示为：

u(t) = S(t)u_0 + \int_0^t S(t-s) f(s,u(s))\,ds.

在函数空间选取上，作者采用以 Stokes 算子 $A = -P\Delta$ （ $$ P $$ 为 Leray 投影）或 Laplace 算子 $-\Delta$ 的分数幂空间 $X_\alpha = D((-A)^\alpha)$ 构成的插值尺度 $\{X_\alpha\}_{\alpha \geq 0}$ 。例如，在 Navier-Stokes 情形中，初始数据取于 $X_1 = D(A) \hookrightarrow H^2(\Omega)^3 \cap H^1_0(\Omega)^3$ ；在反应-扩散情形中，则分别取于 Lebesgue 空间 $L^q(\Omega)$ 或 Besov 空间 $\dot{B}^s_{p,r}(\mathbb{R}^N)$ ，这体现了对初始正则性要求的精细刻画——既非统一要求高阶 Sobolev 正则性，亦非仅限于 $$ L^2 $$ 框架，而是依据非线性结构与记忆核衰减特性进行适配。

此外，边界条件均为齐次 Dirichlet 型（如 $u|_{\partial\Omega} = 0$ ），区域 $\Omega$ 可为有界光滑域或全空间 $\mathbb{R}^N$ ，后者需借助齐次空间与 Littlewood-Paley 分解技术处理。

主要定理或主要结果

本文核心结论集中于三类定理，分别对应解族正则性、局部适定性与临界情形下的 ε-正则化解。

Theorem 2.2（解族的插值尺度正则性估计）：
设 $$ -A $$ 是角为 $\psi_0 \in [0,\pi/2)$ 的扇形算子，核 $$ g $$ 满足假设 (B1)–(B3)，其中关键参数 $\zeta_g > 1$ 满足

\frac{1}{\hat{g}(\lambda)} = O(|\lambda|^{\zeta_g - 1}), \quad |\lambda| \to \infty,\ \lambda \in \Sigma(\omega,\eta+\pi/2).

则对任意 $0 \leq \theta \leq \gamma \leq 1$ 及 $\tau_0 > 0$ ，存在常数 $$ M > 0 $$ ，使得

\|S(t)\|_{\mathcal{B}(X_\gamma,X_{1+\theta})} \leq M t^{-\zeta_g(1+\theta-\gamma)}, \quad t \in (0,\tau_0].

该估计揭示了记忆核主导的正则化速率：指数 $\zeta_g$ 越大，衰减越快，正则性提升越显著。

Theorem 3.1（超临界情形局部适定性）：
设 $\gamma > 1 - \frac{1}{\zeta_g}$ ， $u_0 \in X_1$ ，且 $f(t,\cdot): X_1 \to X_\gamma$ 局部 Lipschitz 连续，满足

\|f(t,x) - f(t,y)\|_{X_\gamma} \leq c(r) \|x-y\|_{X_1}, \quad \forall x,y \in X_1,\ \|x\|_{X_1},\|y\|_{X_1} \leq r.

则存在 $\tau > 0$ ，唯一 mild 解 $u \in C([0,\tau];X_1)$ ，且满足连续依赖性与最大存在时间 $\tau_{\max}$ 的爆破替代：若 $\tau_{\max} < \infty$ ，则

\limsup_{t \to \tau_{\max}^-} \|u(t)\|_{X_{1+\theta}} = +\infty, \quad \forall \theta \in \left(0,\gamma - 1 + \frac{1}{\zeta_g}\right).

进一步， $u \in C((0,\tau];X_{1+\theta})$ ，且对任意紧集 $J \subset X_1$ ，有

\lim_{t \to 0^+} t^{\zeta_g \theta} \sup_{u_0 \in J} \|u(t;u_0)\|_{X_{1+\theta}} = 0.

Theorem 3.4 & 3.5（临界情形 ε-正则化解）：
当 $\gamma \leq 1 - \frac{1}{\zeta_g}$ 时，引入 ε-正则映射概念：存在 $\varepsilon \in (0,1)$ ， $\gamma(\varepsilon) \in (0,1)$ ，使得 $f(t,\cdot): X_{1+\varepsilon} \to X_{\gamma(\varepsilon)}$ 局部 Lipschitz。则对任意 $u_0 \in X_1$ ，存在 $\tau > 0$ 及唯一 ε-正则 mild 解 $$ u $$ ，满足

u \in C((0,\tau];X_{1+\varepsilon}), \quad \lim_{t \to 0^+} t^{\zeta_g \varepsilon} \|u(t)\|_{X_{1+\varepsilon}} = 0.

且该解可延拓至最大时间，并在 $\tau_{\max} < \infty$ 时发生 $X_{1+\varepsilon}$ -范数爆破。

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

a) 核心分析学工具

本文构建适定性理论的基石是不动点方法与插值空间理论的深度结合。具体而言：

Banach 不动点定理被应用于映射
$(T\psi)(t) = S(t)u_0 + \int_0^t S(t-s) f(s,\psi(s))\,ds$

在加权连续空间 $\mathcal{X}_\theta = \{ \psi \in C((0,\tau];X_{1+\theta}) : \sup_{t \in (0,\tau]} t^{\zeta_g \theta} \|\psi(t)\|_{X_{1+\theta}} < \infty \}$ 上。权重 $t^{\zeta_g \theta}$ 的引入，正是为了吸收解族 $$ S(t) $$ 的奇异性 $t^{-\zeta_g(1+\theta-\gamma)}$ ，使非线性积分项收敛。
抽象插值理论（real/complex interpolation functors）是贯穿全文的框架性工具。作者采用 $X_\alpha = F_\alpha(X_k,X_{k+1})$ （ $k \in \mathbb{N}_0$ ）定义插值尺度，并反复运用重迭代定理（reiteration theorem）与插值不等式（interpolation inequality）：
$\|T\|_{\mathcal{B}(X_\gamma,X_{1+\theta})} \leq C \|T\|_{\mathcal{B}(X_0,X_1)}^{1+\theta-\gamma} \|T\|_{\mathcal{B}(X_0)}^{\gamma-\theta},$

将解族在任意插值空间上的范数控制归结为其在端点空间 $\mathcal{B}(X_0,X_1)$ 与 $\mathcal{B}(X_0)$ 上的估计，后者由 Theorem 2.1 提供。
Laplace 变换与复分析技术用于建立解族的解析延拓。通过验证 $\lambda/\hat{g}(\lambda)$ 落入 $$ -A $$ 的预解集 $\rho(-A)$ ，并利用扇形算子的预解估计
$\|(\lambda + A)^{-1}\|_{\mathcal{B}(X_0)} \leq \frac{C}{|\lambda|}, \quad \lambda \in \mathbb{C}\setminus \Sigma[0,\psi_1],$

结合 Cauchy 积分公式，构造出解析解族 $$ S(z) $$ ，进而导出其在扇形内的增长控制。

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

非线性项的奇异性处理：当 $\gamma \leq 1 - \frac{1}{\zeta_g}$ ，标准 Lipschitz 条件导致卷积积分发散，因
$\int_0^t (t-s)^{-\zeta_g(1-\gamma)} ds = \infty.$

文章创新性地引入 Arrieta-Carvalho 型 ε-正则映射，即允许非线性作用于更高正则空间 $X_{1+\varepsilon}$ ，同时要求解在 $t \to 0^+$ 时满足加权初值条件 $t^{\zeta_g \varepsilon} \|u(t)\|_{X_{1+\varepsilon}} \to 0$ 。该条件本质上是对解在原点附近“瞬时正则化强度”的量化约束，确保非线性积分收敛。
记忆核的精细刻画：针对物理中常见的幂律核 $g(t) = k t^{\alpha-1} e^{c t}$ ，作者在附录中严格验证其 Laplace 变换 $\hat{g}(\lambda)$ 满足 (B3) 中 $\zeta_g = 1+\alpha$ 的渐近行为，从而将抽象理论与具体模型无缝衔接。此分析依赖于 Gamma 函数的渐近展开与复平面上分支点的精确处理。
一致紧性论证：为证明解关于初值的连续依赖性及 $t \to 0^+$ 的一致小量估计（Proposition 2.2），作者利用 $X_1 \hookrightarrow X_{1+\beta}$ 的紧嵌入（当 $\Omega$ 有界时），结合 $$ S(t) $$ 的衰减估计，对紧集 $J \subset X_1$ 进行有限覆盖，再逐片控制，避免了对初值整体正则性的过度依赖。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重实质性推进：

第一，对初始数据正则性要求的系统性降低与精细化分类。
经典工作如 Fujita-Kato（Navier-Stokes）或 Brezis-Cazenave（热方程）均要求初始数据位于 $$ H^1 $$ 或 $$ L^p $$ 的临界空间（如 $\dot{H}^{1/2}$ 或 $$ L^N $$ ）。而本文通过插值尺度 $X_\alpha$ 与参数 $\zeta_g$ ，将适定性阈值明确划分为：

超临界区 $\gamma > 1 - \frac{1}{\zeta_g}$ ：只需 $u_0 \in X_1$ （即 $H^2 \cap H^1_0$ 类），解自动获得 $X_{1+\theta}$ 正则性；
临界区 $\gamma = 1 - \frac{1}{\zeta_g}$ ：引入 ε-正则性，允许 $u_0 \in X_1$ ，但解在 $$ t>0 $$ 即属 $X_{1+\varepsilon}$ ；
亚临界区 $\gamma < 1 - \frac{1}{\zeta_g}$ ：传统方法失效，本文首次给出可操作的 ε-正则化解理论。
此三分法远超 Barbu-Sritharan（仅 $$ L^2 $$ 局部解）与 de Andrade et al.（仅幂律核 $\alpha \in [0,1)$ ）的局限，构成首个适用于一般 Laplace 可积核的统一框架。

第二，解的正则性提升幅度由材料函数定量决定。
Theorem 2.2 显示，正则化速率 $t^{-\zeta_g(1+\theta-\gamma)}$ 直接由核的 Laplace 渐近阶 $\zeta_g$ 控制。例如，当 $g(t) \sim t^{\alpha-1}$ ，则 $\zeta_g = 1+\alpha$ ，故 $\alpha \to 0$ （接近 Newtonian 粘性）时，正则化最弱； $\alpha \to 1^-$ （强记忆）时，正则化最强。这一材料-正则性对应原理此前仅在特定模型（如 [17]）中零星出现，本文将其提升为普适规律。

第三，应用广度与空间框架的拓展。
先前关于记忆型 PDE 的工作多局限于 $$ L^2 $$ 或 $$ L^p $$ 框架（如 [14,36]），而本文首次将 Besov 空间（用于梯度型非线性）与 Lebesgue 空间（用于超线性）统一纳入插值尺度分析，并通过重迭代定理建立 $X_\alpha$ 与 $\dot{B}^s_{p,r}$ 、 $H^{s,p}$ 的等价嵌入，从而将 Navier-Stokes 与反应-扩散问题置于同一抽象平台，实现跨模型的正则性理论整合。

3. 由分数布朗运动驱动的快慢Navier-Stokes系统的平均化动力学与Wong-Zakai逼近

Averaging Dynamics and Wong-Zakai approximations for a Fast-Slow Navier-Stokes System Driven by fractional Brownian Motion

作者: Eliseo Luongo, Francesco Triggiano

类别: math.PR, math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08680v1

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

详细分析

中文标题: 由分数布朗运动驱动的快慢Navier-Stokes系统的平均化动力学与Wong-Zakai逼近

1. 研究对象和背景

本文研究一类耦合的快慢Navier-Stokes系统，其核心特征在于：快变量受加性分数布朗运动（fractional Brownian motion, fBm）驱动，而慢变量通过非线性耦合受快变量影响。该系统建模了湍流中多尺度相互作用的随机本质，尤其关注非马尔可夫、具有长程记忆效应的噪声对宏观流体行为的平均化效应。

系统定义在 $$d$$ -维环面 $\mathbb{T}^d$ （ $$d=2$$ 或 $$3$$ ）上，所有向量场均满足零均值与不可压缩条件（即 $\nabla \cdot u = 0$ , $\nabla \cdot v = 0$ ）。空间函数属于Sobolev型Hilbert空间 $H = \{ u \in L^2(\mathbb{T}^d;\mathbb{R}^d) : \nabla \cdot u = 0,\, \int_{\mathbb{T}^d} u = 0 \}$ ，并进一步引入Leray投影算子 $\Pi$ 和Stokes算子 $A = \Pi \Delta$ 。噪声为发散自由、空间相关的 $$H$$ -值分数Wiener过程 $W^H = \sum_{i} \sigma_i e_i W^{H,i}_t$ ，其中 $\{e_i\}$ 是 $H^\xi$ 中正交基（ $\xi > d/2 + 2$ ）， $\{\sigma_i\}$ 满足迹类条件，且协方差算子 $Q = \sum_i \sigma_i^2 e_i \otimes e_i$ 与阻尼算子 $$C$$ 可交换； $$C$$ 是定义在 $$L^2$$ 上的线性闭算子，满足谱条件 $\sup_i \lambda_i = \lambda_1 < 0$ 及负定性 $-\langle w, C w \rangle \gtrsim \|w\|_{H^\gamma}^2$ 。

系统原始形式（见第1页公式(1.5)）为：

\begin{cases} du^\epsilon = \nu \Delta u^\epsilon \, dt - (u^\epsilon + v^\epsilon) \cdot \nabla u^\epsilon \, dt + \nabla p^\epsilon \, dt, \\ dv^\epsilon = \big( \nu \Delta v^\epsilon + \epsilon^{-1} C v^\epsilon - (u^\epsilon + v^\epsilon) \cdot \nabla v^\epsilon + \nabla q^\epsilon \big) \, dt + \epsilon^{-\alpha} dW^H_t, \\ \nabla \cdot u^\epsilon = 0,\quad \nabla \cdot v^\epsilon = 0. \end{cases}

此处 $\epsilon \ll 1$ 是尺度分离参数，控制快变量 $v^\epsilon$ 的时间尺度远快于慢变量 $u^\epsilon$ ； $\alpha \in \{1,\,1+H\}$ 是噪声强度标度，其选择取决于 Hurst 参数 $H \in \left(\frac{1}{3},1\right)$ ，体现两种物理机制的竞争：Wong–Zakai近似（对应 $\alpha = 1$ ）与多尺度分析（对应 $\alpha = 1 + H$ ）。特别地，当 $H = \frac{1}{2}$ 时二者重合，退化为经典布朗情形。

值得注意的是，该系统不采用Stratonovich或Itô积分解释，而是在粗糙路径框架下构造解。原因在于：当 $H \neq \frac{1}{2}$ 时，fBm 不是半鞅，Itô微积分失效；且其样本轨道 Hölder 正则性仅为 $C^{H-\delta}$ ，对 $H \in \left(\frac{1}{3},\frac{1}{2}\right)$ ，甚至无法定义 Young 积分，必须诉诸 $$p$$ -粗糙路径理论（ $$p > 1/H$$ ）。因此，解被理解为无界粗糙驱动器（unbounded rough driver, URD）意义下的弱解（Definition 2.3），其核心是将非线性项 $b(v^\epsilon, u^\epsilon) = \Pi(v^\epsilon \cdot \nabla u^\epsilon)$ 重写为关于辅助过程 $y^\epsilon_t = \int_0^t \epsilon^{-\alpha + H} w^\epsilon_s \, ds$ 的粗糙积分。

为解析快慢解耦，作者采用关键分解：

v^\epsilon = \epsilon^{-\alpha + H} w^\epsilon + r^\epsilon,

其中 $w^\epsilon$ 是分数Ornstein-Uhlenbeck（fOU）过程，满足线性SDE：

dw^\epsilon_t = \epsilon^{-1} C w^\epsilon_t \, dt + \epsilon^{-H} dW^H_t,

而余项 $r^\epsilon$ 承载高阶非线性效应。该分解将原系统转化为等价三元组 $(u^\epsilon, w^\epsilon, r^\epsilon)$ 动力学（见第4页公式(1.8)），使得后续对 $w^\epsilon$ 的遍历性分析与对 $r^\epsilon$ 的能量估计成为可能。

整个分析立足于概率泛函空间：对每个 $\epsilon \in (0,1)$ ，存在一个概率空间 $(\Omega^\epsilon, \mathcal{F}^\epsilon, \mathbb{P}^\epsilon)$ 支撑 $$Q$$ -分数Wiener过程 $W^{\epsilon,H}$ ，并在其上构造概率弱解（probabilistically weak solution）。在二维情形下，因解唯一性成立，可提升为概率强解（probabilistically strong solution）。

2. 主要定理或主要结果

论文的核心结论由以下定理与命题给出：

Theorem 1.1（主定理）
设 $H \in \left(\frac{1}{3},1\right)$ ，初始数据 $u^\epsilon_0, v^\epsilon_0$ 满足：

$\nabla \cdot u^\epsilon_0 = \nabla \cdot v^\epsilon_0 = 0$ ，
$\int_{\mathbb{T}^d} u^\epsilon_0 = \int_{\mathbb{T}^d} v^\epsilon_0 = 0$ ，
$\sup_{\epsilon \in (0,1)} \left( \|u^\epsilon_0\|_{L^2_x} + \sqrt{\epsilon}\, \|v^\epsilon_0\|_{L^2_x} \right) < +\infty$ ，
且算子 $$Q, C$$ 满足第2.2节假设（可交换、负定、迹类等）。则存在概率弱粗糙路径解 $(u^\epsilon, (-C)^{-1} dW^H, \epsilon)$ 满足系统(1.8)，且当 $\epsilon \to 0$ 时，依分布收敛于如下极限方程的解：
若 $\alpha = 1$ 且 $H > \frac{1}{2}$ ，则 $u^\epsilon \xrightarrow{d} u$ ，其中 $$u$$ 是如下含输运噪声的Young型PDE的弱解：
$\begin{cases} du = (\nu \Delta u - u \cdot \nabla u) \, dt + d\nabla p - d\big[(-C)^{-1} W^H\big] \cdot \nabla u, \\ \nabla \cdot u = 0. \end{cases}$
若 $\alpha = 1 + H$ 且 $H < \frac{1}{2}$ ，则 $u^\epsilon \xrightarrow{d} u$ ，其中 $$u$$ 是如下确定性PDE附加Itô–Stokes漂移项的弱解：
$\begin{cases} \partial_t u = \nu \Delta u - u \cdot \nabla u - r \cdot \nabla u + \nabla p, \\ \nabla \cdot u = 0, \end{cases}$
其中漂移场 $r \in H$ 显式依赖于 $$C, Q, H$$ （见Lemma 3.4），且由二次泛函 $\int_0^T b(w^\epsilon_s, w^\epsilon_s) \, ds$ 的遍历极限给出。

Remark 1.2（二维强化）
当 $$d = 2$$ 时，因预极限与极限方程均具唯一性，上述依分布收敛可加强为依概率收敛，且解可在任意支撑 $$W^H$$ 的概率空间上构造（概率强解）。

Proposition 1.4（特殊噪声下的Wong–Zakai扩展）
若噪声具有分离结构 $W^H_t(x) = f(x) W^{H,1}_t$ ，其中 $$f$$ 满足 $\nabla \cdot f = 0$ 、 $\int f = 0$ 且 $f \cdot \nabla f = 0$ （如 $f(x) = a \sin(l \cdot x)$ ， $a \perp l$ ），则对任意 $H \in \left(\frac{1}{3}, \frac{1}{2}\right)$ 与 $\alpha = 1$ ，Theorem 1.1 中 $H > \frac{1}{2}$ 的Wong–Zakai型收敛仍成立。

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

a) 核心分析学工具

粗糙路径理论（Rough Path Theory）：
鉴于 $H \in \left(\frac{1}{3},1\right)$ ，fBm 的正则性不足，传统随机微积分失效。作者采用 $$p$$ -粗糙路径框架（ $p \in (1/H, 3)$ ），将 $y^\epsilon_t = \int_0^t \epsilon^{-\alpha + H} w^\epsilon_s \, ds$ 提升为其典范 $$p$$ -粗糙路径 lift $(Y^{1,\epsilon}, Y^{2,\epsilon})$ ，满足 Chen 关系 $\delta Y^{2,\epsilon}_{s,u,t} = Y^{1,\epsilon}_{s,u} \otimes Y^{1,\epsilon}_{u,t}$ 。此 lift 构成无界粗糙驱动器（URD） $(A^{1,\epsilon}, A^{2,\epsilon})$ ，其中
$A^{1,\epsilon}_{s,t} \phi = b(Y^{1,\epsilon}_{s,t}, \phi), \quad A^{2,\epsilon}_{s,t} \phi = b^{(2)}(Y^{2,\epsilon}_{s,t}, \phi),$

用于定义 Davie 型粗糙积分 $\int_s^t b(dy^\epsilon_r, u^\epsilon_r) = \sum_{k=1}^2 A^{k,\epsilon}_{s,t} u^\epsilon_s + u^{\natural,\epsilon}_{s,t}$ 。该框架统一处理 $H < \frac{1}{2}$ （需 rough integral）与 $H > \frac{1}{2}$ （可降为 Young integral）情形。
Gaussian Rough Path 理论：
对 $H < \frac{1}{2}$ ，fOU 过程 $w^\epsilon$ 的协方差结构复杂，无法直接应用 Markov 遍历定理。作者调用 Gaussian rough path 理论（[31,30]），利用 fBm 的自相似性与分数协方差核 $R_H(s,t) = \frac{1}{2}(|s|^{2H} + |t|^{2H} - |t-s|^{2H})$ ，结合 Malliavin 微积分技巧，严格构造 $Y^{1,\epsilon}, Y^{2,\epsilon}$ 的收敛性。
Galerkin 逼近与随机紧性（Stochastic Compactness）：
对固定 $\epsilon$ ，通过截断 Fourier 模式构造 Galerkin 近似解 $(u^{\epsilon,N}, w^{\epsilon,N}, r^{\epsilon,N})$ ，建立一致先验估计（见第5节）：
$\sup_{\epsilon} \mathbb{E}\left[ \|u^\epsilon\|_{L^\infty_t L^2_x}^2 + \|u^\epsilon\|_{L^2_t H^1_x}^2 \right] < \infty, \quad \sup_{\epsilon} \mathbb{E}\left[ \|w^\epsilon\|_{L^2_t H^\gamma_x}^2 \right] < \infty.$

利用 Aldous 准则与 Jakubowski–Skorokhod 定理，在广义 Skorokhod 空间中提取子列，获得概率弱解的存在性。
遍历性与二次泛函收敛（Ergodicity of Quadratic Functionals）：
关键突破在于处理 $r^\epsilon$ 的极限。作者证明：当 $H < \frac{1}{2}$ ， $\int_0^T b(w^\epsilon_s, w^\epsilon_s) \, ds$ 的渐近行为由 fOU 的分数遍历定理控制。借助 $$C$$ 与 $$Q$$ 的可交换性，将无限维问题约化为一维 fOU 的二次变差估计（Lemma 3.1）：
$\lim_{\epsilon \to 0} \mathbb{E}\left[ \left\| \int_0^T b(w^\epsilon_s, w^\epsilon_s) \, ds - T r \right\|_{H^{-1}}^2 \right] = 0,$

其中 $r = \sum_i \sigma_i^2 \lambda_i^{-1} b(e_i, e_i)$ 即 Itô–Stokes 漂移。

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

非线性项 $b(v^\epsilon, u^\epsilon)$ 的奇异性控制：
直接估计 $b(v^\epsilon, u^\epsilon)$ 在 $L^1_t H^{-1}_x$ 中发散。作者通过分解 $v^\epsilon = \epsilon^{-\alpha + H} w^\epsilon + r^\epsilon$ ，将该项拆为 $b(\epsilon^{-\alpha + H} w^\epsilon, u^\epsilon) + b(r^\epsilon, u^\epsilon)$ 。前者由粗糙积分理论控制；后者通过 $r^\epsilon$ 的高阶能量估计（ $\|r^\epsilon\|_{L^2_t H^\gamma} \lesssim \epsilon^{1/2}$ ）吸收。
尺度参数 $\alpha$ 的双重角色处理：
$\alpha$ 同时影响噪声强度与 fOU 的扩散速率。作者通过精细比较 $\epsilon^{-\alpha + H}$ 与 $\epsilon^{-1}$ 的相对大小，发现：当 $\alpha = 1$ ， $w^\epsilon$ 快速弛豫至“分数白噪声”；当 $\alpha = 1 + H$ ， $w^\epsilon$ 缓慢演化，其二次效应累积为确定性漂移。这一 dichotomy 由 $$H$$ 与 $\frac{1}{2}$ 的大小关系决定。
Hurst 参数临界值 $H = \frac{1}{2}$ 的统一处理：
为避免对 $H < \frac{1}{2}$ 与 $H > \frac{1}{2}$ 分别构造不同积分理论，作者全程采用 Davie 的 URD 公式，其本质是将粗糙积分视为有界线性泛函作用于测试函数，从而绕过路径逐点定义的困难。

4. 与之前工作的比较

a) 突破与改进

首次建立 fBm 驱动快慢NS系统的平均化理论：
此前工作（如 [13,15,23,24]）均假设驱动噪声为标准布朗运动（ $H = \frac{1}{2}$ ），依赖 Itô 公式与 Markov 性质。本文首次将分析拓展至全范围 $H \in \left(\frac{1}{3},1\right)$ ，揭示 Hurst 参数如何根本性改变平均化极限的结构： $H > \frac{1}{2}$ 产生输运噪声， $H < \frac{1}{2}$ 产生 Itô–Stokes 漂移， $H = \frac{1}{2}$ 二者共存。
放宽对噪声结构的各向同性假设：
文献 [23,24] 要求协方差算子 $$Q$$ 各向同性以消去漂移项。本文证明：漂移项的出现与否取决于 $$H$$ ，而非 $$Q$$ 的对称性；即使 $$Q$$ 各向异性，只要 $H < \frac{1}{2}$ ，漂移 $$r$$ 依然非零且显式可算。
解决非马尔可夫框架下的 Wong–Zakai 逼近问题：
经典 Wong–Zakai 定理（[29]）仅适用于半鞅噪声。本文证明：对 fBm，当 $\alpha = 1$ ，fOU 过程 $w^\epsilon$ 的典范 lift 依粗糙路径度量收敛于 $$W^H$$ 的 lift，从而确立分数 Wong–Zakai 近似的有效性，这是对随机分析基础理论的重要拓展。
克服多尺度分析中非线性耦合的障碍：
相比 [42] 中简化的快变量线性方程 $dv^\epsilon = \epsilon^{-1} C v^\epsilon dt + \epsilon^{-1} dW^H_t$ ，本文保留完整的非线性项 $-(u^\epsilon + v^\epsilon)\cdot \nabla v^\epsilon$ ，技术难度剧增。作者通过 Proposition 1.4 给出在特定物理可实现噪声（如单模 $f(x) W^{H,1}_t$ ）下的完整结果， bridging 理论严谨性与建模现实性。

综上，本文不仅填补了分数噪声下流体平均化理论的空白，更发展了一套融合粗糙路径、Gaussian 分析与流体力学的通用框架，为研究具有记忆效应的复杂系统提供了新范式。

4. 几何Navier-Stokes方程的反问题

Inverse problem for the geometric Navier-Stokes equations

作者: Yavar Kian, Lauri Oksanen, Ziyao Zhao

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08644v1

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

详细分析

中文标题: 几何Navier-Stokes方程的反问题

研究对象和背景

本文研究的是定义在紧致带边黎曼流形 $$(M,g)$$ 上的几何Navier–Stokes方程（geometric Navier–Stokes equations），其核心是将经典欧氏空间中的不可压缩流体模型推广至具有非平凡曲率的弯曲时空结构中。该方程系统不仅保留了物理上关键的不可压缩性约束（即速度场的余散度为零），还通过Hodge拉普拉斯算子 $\Delta_H$ 和Levi-Civita联络诱导的非线性对流项 $$R(u)$$ ，自然地融合了流形的内蕴几何结构。

具体而言，论文在第3页明确给出了所研究的初边值问题：

\begin{cases} \partial_t u - \Delta_H u + R(u) + dp = f, & \text{in } [0,T] \times M, \\ d^* u = 0, & \text{in } [0,T] \times M, \\ u|_{[0,T]\times \partial M} = 0, & \\ u|_{\{0\}\times M} = 0,\quad p|_{\{0\}\times M} = 0. & \end{cases}

其中：

$$u$$ 是定义在 $[0,T] \times M$ 上的1-形式（对应于切向量场的速度场）；
$$p$$ 是标量压力函数；
$$f$$ 是给定的局部源项，支撑集包含于 $(0,T] \times \omega$ ， $\omega \subset M$ 是一个非空开子集；
$\Delta_H = -dd^* - d^*d$ 是作用于微分形式上的Hodge拉普拉斯算子；
$$d^*$$ 是余微分算子，由Hodge星算子 $\ast$ 定义为
$d^* = (-1)^{n(k+1)+1} \ast d \ast,$ 其中 $n = \dim M$ , $$k=1$$ （因 $$u$$ 是1-形式）；
$R(u) = (\nabla_{u^*} u^*)^*$ ，其中 $\nabla$ 是 $$(M,g)$$ 上的Levi-Civita联络， $$u^*$$ 是与向量场 $$u$$ 对应的1-形式（通过度量 $$g$$ 的升降指标实现），因此 $$R(u)$$ 是一种几何化的对流项，体现曲率对流体运动的影响；
边界条件 $u|_{\partial M}=0$ 是无滑移边界条件（no-slip condition），在微分形式语言中等价于切向分量消失；
初始条件 $u(0,\cdot)=0$ 、 $p(0,\cdot)=0$ 表明系统从静止状态启动。

该方程组定义在Hölder空间 $C^{1,\alpha}([0,T]; C^{2,\alpha}\Omega^1(M))$ 中，而非通常的Sobolev框架，这为后续非线性分析与线性化技术提供了更精细的正则性控制。特别地，作者引入了如下函数空间：

$W_0 := \{ u \in C^\infty\Omega^1(M) \mid d^* u = 0 \}$ ，
$V := \{ v \in H^1\Omega^1(M) \mid d^* v = 0 \}$ ，
$$H$$ ： $$W_0$$ 在 $L^2\Omega^1(M)$ 中的闭包。

这些空间构成经典的de Rham复形中满足散度自由（divergence-free）约束的子空间，是处理不可压缩流体问题的自然框架。值得注意的是，由于 $$d^* u = 0$$ ，压力项 $$dp$$ 不再是独立变量，而是由正交分解唯一确定（见Lemma 2.2）：对任意 $f \in H^{-1}\Omega^1(M)$ ，若 $\langle f,v\rangle = 0$ 对所有 $v \in W_0$ 成立，则存在 $p \in L^2(M)$ 使得 $$f = dp$$ 。

此外，区域 $$M$$ 是紧致、光滑、连通、带边黎曼流形，其边界 $\partial M$ 非空；观测仅在时间终态 $$t=T$$ 、且仅限于空间局部区域 $\omega \subset M$ 内进行，即所谓“固定时刻局部观测”（fixed-time local observation）。这种设定显著弱于传统反问题中要求的全边界测量（如Dirichlet-to-Neumann映射），也不同于封闭流形上的源到解映射，因而更具挑战性与现实意义——例如在地球物理中，仅能获取地表某小片区域在某一时刻的风速数据，却需推断整个大气层的几何结构（如地形曲率、有效重力场）。

综上，本研究对象是一个高度几何化、强非线性、带约束的抛物型系统，其解的存在性、唯一性依赖于流形的整体几何性质（如曲率、拓扑、谱特性），而反问题的目标正是从极有限的局部终值数据中唯一恢复流形的等距类（isometry class），即其完整的黎曼度量结构。

主要定理或主要结果

论文的核心成果是以下唯一性定理，出现在第3页的Theorem 1.1：

Theorem 1.1. 设 $$(M_1,g_1)$$ 与 $$(M_2,g_2)$$ 是两个 $$n$$ 维光滑连通带边黎曼流形， $\omega_j \subset M_j$ 是非空连通开集（ $$j=1,2$$ ）。记 $$N_j$$ 为对应于 $$(M_j,g_j)$$ 与 $\omega_j$ 的局部源到终值映射（local source-to-final-value map）：
$N_j : B_\delta \to L^2\Omega^1(\omega_j), \quad N_j(f) = u_f(T,\cdot)|_{\omega_j},$
其中 $$u_f$$ 是方程 (2) 在源 $f \in B_\delta = \{ g \in C^\alpha(M) : \operatorname{supp}(g) \subset (0,T]\times \omega_j,\ \|g\|_{C^\alpha} \le \delta \}$ 下的唯一解。若存在微分同胚 $\Phi : \omega_1 \to \omega_2$ 满足
$N_1 \circ \Phi^* = \Phi^* \circ N_2,$
则 $$(M_1,g_1)$$ 与 $$(M_2,g_2)$$ 等距同构（isometric）。

该定理的假设极为精炼：仅需两个流形上定义的局部观测数据算子在某个小邻域间通过拉回映射 $\Phi^*$ 相容，即可推出全局几何等价。结论的强度在于它不依赖于任何先验的流形拓扑或曲率信息，亦无需知道维数 $$n$$ （只要两者相同），且适用于任意 $n \ge 2$ 。

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

a) 核心分析学工具

论文的证明骨架建立在三大支柱之上：非线性系统的线性化理论、椭圆/抛物系统的谱分析、以及边界控制（Boundary Control, BC）方法的几何化推广。

首先，在非线性处理上，作者采用高阶Fréchet导数线性化策略。设 $$u_f$$ 是方程 (2) 关于源 $$f$$ 的解，则考虑其在 $$f=0$$ 处的一阶变分 $v = \partial_f u_f|_{f=0} \cdot h$ ，可导出线性化Stokes系统：

\begin{cases} \partial_t v - \Delta_H v + R'(0)[v] + dq = h, \\ d^* v = 0, \\ v|_{\partial M} = 0,\quad v(0,\cdot)=0. \end{cases}

此处 $$R'(0)[v]$$ 是 $$R(u)$$ 在 $$u=0$$ 处的导数，因 $$R(u)$$ 含联络系数，故仍具几何本质。Proposition 5.2 严格证明了：原始非线性映射 $$N$$ 的一阶线性化，恰好等于该线性Stokes系统的源到解映射 $$P$$ 在终值 $$t=T$$ 上的限制。这一结论是连接非线性反问题与线性反问题的桥梁。

其次，在谱分析层面，作者深入研究了定常Stokes算子 $$L$$ （定义于第8页）：

L : L^2\Omega^1(M) \to V,\quad \langle d(Lf), dw\rangle = \langle f, w\rangle,\ \forall w \in V.

利用Lax-Milgram定理及紧嵌入 $H^1 \hookrightarrow L^2$ ，证明 $$L$$ 是紧自伴正定算子，从而拥有离散谱 $\{1/\lambda_j\}_{j=1}^\infty$ ，对应标准正交特征1-形式 $\{\phi_{jk}\}$ 满足：

-\Delta_H \phi_{jk} + dp_{jk} = \lambda_j \phi_{jk},\quad d^*\phi_{jk}=0.

关键洞见在于：非定常Stokes系统的源到解算子 $$P$$ ，可通过 $$L$$ 的谱投影精确表示（见Section 5）。具体地，对任意 $h \in C^\alpha$ ，有

P(h)(T,\cdot) = \sum_{j=1}^\infty e^{-\lambda_j T} \sum_{k=1}^{m_j} \langle h, \phi_{jk}\rangle \phi_{jk},

即 $$P$$ 本质上是 $$L$$ 的热半群 $e^{-tL}$ 在 $$t=T$$ 的作用。这一步将抛物问题转化为椭圆算子的谱数据，为下一步引入双曲类比铺平道路。

最后，也是最富创新性的部分，是将上述谱表示解析延拓至虚时间，构造一个辅助双曲Stokes系统（auxiliary hyperbolic Stokes system）：

\partial_{tt} w - \Delta_H w + R'(0)[w] + dr = 0,

其对应的波传播算子 $\mathcal{U}(t) = \cos(t\sqrt{L})$ 具有有限传播速度（finite speed of propagation），这是BC方法赖以生存的物理基础。作者在Section 6中严格建立了该双曲系统的条件可控性（conditional controllability）：对任意 $$T_0 > 0$$ ，存在时间 $$T < T_0$$ （取决于流形直径与曲率下界），使得由 $\omega$ 出发的控制造波可在 $$T$$ 时刻覆盖 $$M$$ 的任意点。此性质虽弱于全局可控性，但已足以支撑BC框架。

b) 技术性难题的克服

最大难点源于散度自由约束 $$d^* u = 0$$ 。它导致：

解空间 $$V$$ 非全空间，传统BC方法中基于 $$H^1$$ 能量的控制论工具失效；
压力项 $$dp$$ 无法被消去，必须作为拉格朗日乘子显式处理；
特征函数系 $\{\phi_{jk}\}$ 仅在 $$V$$ 中正交，而非在 $$L^2$$ 中，使投影算子缺乏标准正交基的简洁性。

作者通过三项关键技术突破：

几何正交分解引理（Lemma 2.2）：将任意 $f \in H^{-1}\Omega^1$ 分解为 $f = f_{\text{div}} + dp$ ，其中 $f_{\text{div}} \in H$ （散度自由部分）， $$dp$$ 为梯度场。这保证了 $$L$$ 作用于 $f_{\text{div}}$ 时压力自动被吸收。
特征函数局部线性无关性（Lemma 3.3）：证明 $\{\phi_{jk}|_\omega\}$ 在任意开集 $\omega$ 上线性无关，从而确保局部观测能激发全部模态，避免“盲区”。
BC方法的向量丛适配：将经典BC中的“边界响应算子”替换为局部源激励下的内部响应算子 $\mathcal{R}_\omega$ ，并定义其“可控子空间”为 $\mathcal{K}_\omega(T) := \overline{\operatorname{span}}\left\{ \mathcal{U}(t)h : h \in C^\infty_0(\omega),\ t \in [0,T] \right\} \subset H.$ 利用流形的距离函数 $\rho(x,y) = \inf\{\text{length of curves from }x\text{ to }y\}$ ，证明 $\mathcal{K}_\omega(T)$ 恰好等于所有满足 $\operatorname{supp}(v) \subset \{x \in M : \rho(x,\omega) \le T\}$ 的 $v \in H$ 。由此，通过比较两个流形的 $\mathcal{K}_{\omega_j}(T)$ ，即可重建其距离函数，进而由Myers–Steenrod定理得到等距。

与之前工作的比较

a) 突破与改进

本文实现了三重根本性突破：

第一，首次将反问题拓展至非欧几里得几何框架。 此前所有关于Navier–Stokes反问题的研究（如[37, 57, 63]）均局限于欧氏区域 $\Omega \subset \mathbb{R}^n$ ，且多针对定常情形或仅重构粘性系数 $\mu$ 。本文则首次处理动态、非线性、弯曲空间中的完整度量重构，将反问题从“参数识别”升维至“几何识别”。

第二，观测方式取得质的简化。 经典结果（如[2, 49]）依赖全边界Dirichlet-to-Neumann映射，即需在全部 $\partial M$ 上施加并测量数据；而本文仅需内部局部区域 $\omega$ 在单一时刻 $$T$$ 的速度场。这种“稀疏、瞬时、内部”观测更贴近实际（如气象卫星仅扫描局部云图），技术难度却指数级上升——因丧失了时间演化与边界反射提供的丰富信息。

第三，方法论上开创BC方法在向量值系统中的应用。 BC方法此前仅成功用于标量波动方程、Maxwell系统（电磁场）、Dirac系统（旋量场），但从未触及Stokes或Navier–Stokes这类带约束的流体系统。本文不仅首次应用，更发展出一套约束兼容的BC理论，包括：约束下的可控子空间刻画、散度自由波的传播锥定义、以及基于Hodge分解的响应算子构造。该框架可推广至其他受约束的物理系统（如弹性力学中的位移-应变相容性条件）。

综上，本文不仅解决了长期悬而未决的几何Navier–Stokes反问题，更开辟了一条将现代几何分析、谱理论与控制论深度融合的新路径。

5. 雷诺数对薄多孔介质中非牛顿流体流动的影响

Influence of the Reynolds number on non-Newtonian flow in thin porous media

作者: Maria Anguiano, Matthieu Bonnivard, Francisco J. Suarez-Grau

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08555v1

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

详细分析

中文标题: 雷诺数对薄多孔介质中非牛顿流体流动的影响

研究对象和背景

本文研究的是广义牛顿流体（generalized Newtonian fluid）在三维薄多孔介质 $\Omega_\varepsilon \subset \mathbb{R}^3$ 中的定常流动，其本构关系由 Carreau 黏度律 描述，控制方程为非线性定常 Navier–Stokes 系统。该系统具有强非线性、低维几何约束（厚度 $\varepsilon \ll 1$ ）与周期微结构（固相圆柱直径亦为 $\varepsilon$ ）三重耦合特征，属于典型的多尺度奇摄动问题。

具体地，设 $\omega \subset \mathbb{R}^2$ 是光滑有界连通区域，定义薄域：

Q_\varepsilon = \omega \times (0,\varepsilon),\quad \Omega_\varepsilon = Q_\varepsilon \setminus S_\varepsilon,

其中 $S_\varepsilon = \bigcup_{k' \in K_\varepsilon} T_{k',\varepsilon}$ 是周期分布的固相圆柱集合，每个 $T_{k',\varepsilon} = (\varepsilon k' + \varepsilon T') \times (0,\varepsilon)$ ， $T' \Subset Y' = (-1/2,1/2)^2$ 为单位胞腔中的固相部分。因此 $\Omega_\varepsilon$ 是一个厚度为 $\varepsilon$ 、孔隙率固定、微结构尺度亦为 $\varepsilon$ 的薄多孔介质。

流体速度场 $u_\varepsilon = (u'_\varepsilon, u_{\varepsilon,3})$ 与压力 $p_\varepsilon$ 满足如下定常 Navier–Stokes 系统（见原文 (2.1) 式）：

\begin{cases} -\mu\varepsilon^\gamma \operatorname{div}\big( \eta_r(D[u_\varepsilon]) D[u_\varepsilon] \big) + (u_\varepsilon \cdot \nabla) u_\varepsilon + \nabla p_\varepsilon = f & \text{in } \Omega_\varepsilon, \\ \operatorname{div} u_\varepsilon = 0 & \text{in } \Omega_\varepsilon, \\ u_\varepsilon = 0 & \text{on } \partial Q_\varepsilon \cup \partial S_\varepsilon. \end{cases}

其中外力 $$f(x) = (f'(x'), 0)$$ ，且 $f' \in L^\infty(\omega)^2$ ；黏度函数 $\eta_r$ 遵循 Carreau 定律（见原文 (1.1) 式）：

\eta_r(D[u]) = (\eta_0 - \eta_\infty)\big(1 + \lambda |D[u]|^2\big)^{\frac{r}{2}-1} + \eta_\infty,\quad 1 < r < +\infty,

其中 $|D[u]|^2 = \operatorname{Tr}(D[u] D[u]^t)$ ， $D[u] = \frac{1}{2}(\nabla u + \nabla u^t)$ 为对称形变张量；参数满足 $\eta_0 > \eta_\infty > 0$ , $\lambda > 0$ 。该黏度模型涵盖两大类非牛顿流体：

当 $$1 < r < 2$$ 时， $\eta_r$ 随剪切率增大而减小，称为 伪塑性流体（shear-thinning）；
当 $$r > 2$$ 时， $\eta_r$ 随剪切率增大而增大，称为 胀塑性流体（shear-thickening）；
当 $$r = 2$$ 时， $\eta_r \equiv \eta_0$ ，退化为牛顿流体。

关键物理参数是雷诺数 $\mathrm{Re}_\varepsilon \sim \mu^{-1}\varepsilon^{-\gamma}$ ，其尺度由指数 $\gamma \in \mathbb{R}$ 控制。本文核心问题是：当 $\varepsilon \to 0$ 时，惯性项 $(u_\varepsilon \cdot \nabla) u_\varepsilon$ 在何种 $\gamma$ 下对宏观均质化模型产生实质性影响？

为处理薄域与周期结构的双重奇异性，作者引入拉伸变量 $y_3 = x_3/\varepsilon$ ，定义重缩放函数：

\tilde{u}_\varepsilon(x',y_3) = u_\varepsilon(x',\varepsilon y_3),\quad \tilde{p}_\varepsilon(x',y_3) = p_\varepsilon(x',\varepsilon y_3),\quad \text{defined on } \widetilde{\Omega}_\varepsilon = \omega_\varepsilon \times (0,1),

并相应定义缩放微分算子 $D_\varepsilon$ , $\operatorname{div}_\varepsilon$ , $\nabla_\varepsilon$ （见原文 (2.2)–(2.3)）。重缩放后系统变为：

\begin{cases} -\mu\varepsilon^\gamma \operatorname{div}_\varepsilon\big( \eta_r(D_\varepsilon[\tilde{u}_\varepsilon]) D_\varepsilon[\tilde{u}_\varepsilon] \big) + (\tilde{u}_\varepsilon \cdot \nabla_\varepsilon) \tilde{u}_\varepsilon + \nabla_\varepsilon \tilde{p}_\varepsilon = f & \text{in } \widetilde{\Omega}_\varepsilon, \\ \operatorname{div}_\varepsilon \tilde{u}_\varepsilon = 0 & \text{in } \widetilde{\Omega}_\varepsilon, \\ \tilde{u}_\varepsilon = 0 & \text{on } \partial \Omega \cup \partial T_\varepsilon. \end{cases}

此形式凸显了惯性项的放大效应：由于 $\nabla_\varepsilon$ 含 $\varepsilon^{-1}\partial_{y_3}$ ，惯性项实际贡献阶为 $\varepsilon^{-1}$ ，需与 $\varepsilon^\gamma$ 平衡。

函数空间设定亦具特殊性：对 $1 < r \leq 2$ ，解属于 $(u_\varepsilon,p_\varepsilon) \in H^1_0(\Omega_\varepsilon)^3 \times L^2_0(\Omega_\varepsilon)$ ；对 $$r > 2$$ ，则属于 $(u_\varepsilon,p_\varepsilon) \in W^{1,r}_0(\Omega_\varepsilon)^3 \times L^{r'}_0(\Omega_\varepsilon)$ ，其中 $$1/r + 1/r' = 1$$ 。这种依赖 $$r$$ 的可积性选择源于 Carreau 黏度在高剪切与低剪切下的渐近行为差异，是分析非线性本构的关键技术前提。

主要定理或主要结果

本文的核心结论由两个主定理给出，分别对应伪塑性/牛顿流体与胀塑性流体，并统一指出临界雷诺数为 $\mathrm{Re}_c \sim 1/\varepsilon$ ，即 $\gamma = 1$ 。

定理 2.2（伪塑性或牛顿流体）：设 $1 < r \leq 2$ ， $\gamma \leq 1$ ，并令

C(r) = \begin{cases} \frac{3}{2}, & \gamma \leq \frac{3}{4}, \\ \frac{3}{2\gamma}, & \frac{3}{4} < \gamma \leq 1. \end{cases}

则存在 $\tilde{u} \in H^1(0,1;L^2(\omega)^3)$ （满足 $\tilde{u}_3 = 0$ 于 $\omega \times \{0,1\}$ ）、 $\tilde{P} \in L^{C(r)}(\omega)$ ，使得重缩放解的延拓 $(\tilde{u}_\varepsilon,\tilde{P}_\varepsilon)$ 满足：

\varepsilon^{\gamma-2}\tilde{u}_\varepsilon \rightharpoonup \tilde{u} \text{ weakly in } H^1(0,1;L^2(\omega)^3),\quad \tilde{P}_\varepsilon \to \tilde{P} \text{ strongly in } L^{C(r)}(\Omega).

定义二维滤速 $V'(x') = \int_0^1 \tilde{u}'(x',y_3)\,dy_3$ ，则：

若 $\gamma < 1$ ， $(V',\tilde{P})$ 是如下线性二维达西律的唯一解： $\begin{cases} V'(x') = \dfrac{1}{\mu \eta_0} A \big(f'(x') - \nabla_{x'} \tilde{P}(x')\big) & \text{in } \omega, \\ \operatorname{div}_{x'} V'(x') = 0 & \text{in } \omega, \\ V'(x') \cdot n = 0 & \text{on } \partial \omega, \end{cases}$ 其中渗透率张量 $A = (A_{ij}) \in \mathbb{R}^{2\times 2}$ 由局部 Stokes 问题 (2.7) 定义。
若 $\gamma = 1$ ， $(V',\tilde{P})$ 满足非线性 Carreau 型二维达西律： $\begin{cases} V'(x') = U\big(f'(x') - \nabla_{x'} \tilde{P}(x')\big) & \text{in } \omega, \\ \operatorname{div}_{x'} V'(x') = 0 & \text{in } \omega, \\ V'(x') \cdot n = 0 & \text{on } \partial \omega, \end{cases}$ 其中 $U(\xi') = \int_{Y_f} w'_{\xi'}(y)\,dy$ ， $w'_{\xi'}$ 由非线性局部问题 (2.10) 给出。

定理 2.5（胀塑性流体）：设 $$r > 2$$ ， $\gamma \leq 1$ ，并令

C(r) = \begin{cases} \frac{3r}{6-r}, & 2 < r < \frac{9}{4}, \\ r', & r \geq \frac{9}{4}. \end{cases}

则结论类似：当 $\gamma < 1$ 时得线性达西律；当 $\gamma = 1$ 时得非线性 Carreau 型达西律，但此时解空间为 $(V',\tilde{P}) \in L^r(\omega)^2 \times (L^{r'}(\omega) \cap W^{1,r'}(\omega))$ ，且局部问题 (2.10) 的解属于 $W^{1,r}_\#(Y_f)^3 \times L^{r'}_\#(Y_f)$ 。

综上，所有情形下，临界雷诺数均为 $\mathrm{Re}_c \sim \varepsilon^{-1}$ ；当 $\gamma \leq 1$ （即 $\mathrm{Re}_\varepsilon \lesssim \varepsilon^{-1}$ ），惯性项在极限中消失，有效模型为达西型；当 $\gamma > 1$ ，惯性项将出现在均质化方程中（文中未完全推导，仅指出预期存在）。

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

a) 核心分析学工具

本文证明严格建立在均匀化理论（homogenization theory）与渐近分析（asymptotic analysis）框架之上，结合多重尺度技术与泛函分析估计。核心工具包括：

Poincaré 与 Korn 不等式（引理 3.1）：在薄多孔域 $\widetilde{\Omega}_\varepsilon$ 上，对任意 $v \in W^{1,q}_0(\widetilde{\Omega}_\varepsilon)^3$ ，有
$\|v\|_{L^q} \leq C\varepsilon \|D_\varepsilon v\|_{L^q},\quad \|D_\varepsilon v\|_{L^q} \leq C \|D_\varepsilon[v]\|_{L^q},$
此不等式是获得速度先验估计 $\|\tilde{u}_\varepsilon\|_{L^2} \lesssim \varepsilon^{2-\gamma}$ 的基石。
De Rham 表示定理：用于从弱形式中重构压力延拓 $P_\varepsilon$ ，关键在于证明 $\langle \nabla p_\varepsilon, R_\varepsilon v \rangle$ 定义了一个有界线性泛函于 $W^{1,q}_0(Q_\varepsilon)^3$ ，从而存在 $P_\varepsilon \in L^{C(r)}(Q_\varepsilon)$ 满足 $\nabla P_\varepsilon = F_\varepsilon$ 。
紧性嵌入与弱收敛理论：通过缩放后的先验估计，得到 $\varepsilon^{\gamma-2}\tilde{u}_\varepsilon$ 在 $H^1(0,1;L^2(\omega)^3)$ 中弱收敛， $\tilde{P}_\varepsilon$ 在 $L^{C(r)}(\Omega)$ 中强收敛，进而对测试函数取极限。
Unfolding 方法（虽未显式命名，但思想贯穿）：通过变量拉伸与周期延拓，将振荡系数问题转化为固定域上的收敛问题，是处理周期微结构的标准技术。

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

最大技术难点在于惯性项的精细估计，因其含双线性结构 $(\tilde{u}_\varepsilon \cdot \nabla_\varepsilon)\tilde{u}_\varepsilon$ ，且 $\nabla_\varepsilon$ 含 $\varepsilon^{-1}$ 导数。原文 (3.14) 将其分解为三项：

\int_{\widetilde{\Omega}_\varepsilon} (\tilde{u}_\varepsilon \cdot \nabla_\varepsilon)\tilde{u}_\varepsilon \cdot \widetilde{R}_\varepsilon v \,dx'dy_3 = -\int \tilde{u}_\varepsilon \otimes \tilde{u}_\varepsilon : D_{x'} \widetilde{R}_\varepsilon v + \frac{1}{\varepsilon}\int \partial_{y_3}\tilde{u}_{\varepsilon,3}\, \tilde{u}_\varepsilon \cdot \widetilde{R}_\varepsilon v + \frac{1}{\varepsilon}\int \tilde{u}_{\varepsilon,3}\, \partial_{y_3}\tilde{u}_\varepsilon \cdot \widetilde{R}_\varepsilon v.

为控制这些项，作者采用插值技巧与分段优化可积性：

对 $1<r\leq 2$ ，利用 $$L^2$$ 速度估计与 $$L^q$$ 测试函数，通过 Hölder 不等式及 $L^2 \hookrightarrow L^{C(r)}$ 嵌入，将惯性项控制为 $C \varepsilon^\alpha \|v\|_{W^{1,q}}$ ，并证明当 $\gamma \leq 1$ 时 $\alpha \geq 0$ ，故可忽略；
对 $$r>2$$ ，因速度本身属 $$L^r$$ ，作者创新性地构造 $$C(r)$$ 使 $$L^r$$ 与 $L^{C(r)}$ 对偶匹配，并对 $$2<r<9/4$$ 使用 $\frac{3r}{6-r}$ 这一“超临界”指数，确保惯性项仍可被吸收。

此外，限制算子 $R_\varepsilon$ （引理 3.4）是连接薄层 $Q_\varepsilon$ 与穿孔域 $\Omega_\varepsilon$ 的桥梁，其构造保证了散度保持与稳定性（估计 (3.5)），是压力延拓可行性的关键。

与之前工作的比较

a) 突破与改进

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

第一，几何模型的统一与强化。此前工作如 Boughanim–Tapiéro [22] 研究无孔薄板（thin slab），得临界雷诺数 $\varepsilon^{-3/2}$ （伪塑性）或 $\varepsilon^{-1-2/r}$ （胀塑性）；Bourgeat–Mikelić [23] 研究常规三维周期多孔介质，得 $\varepsilon^{-3/2}$ 或 $\varepsilon^{-1-r/4}$ 。本文首次处理兼具薄性与周期穿孔的双重结构，发现临界雷诺数统一降为 $\varepsilon^{-1}$ ，表明厚度效应主导于孔隙率效应，这是几何建模上的本质推进。

第二，非线性处理的普适性提升。[22,23] 均要求 $$r$$ 落于特定区间以闭合估计，而本文通过自适应选取压力空间 $L^{C(r)}$ （见 (2.4)(2.11)），对所有 $$r>1$$ 实现完整分析：对 $1<r\leq 2$ ，用 $L^{3/(2\gamma)}$ 改进标准 $$L^2$$ ；对 $$2<r<9/4$$ ，用 $\frac{3r}{6-r}>r'$ 突破传统对偶限制。这使得结论不再依赖 $$r$$ 的“友好”取值，显著拓宽适用范围。

第三，物理洞见的深化。本文明确指出：当 $\mathrm{Re}_\varepsilon \lesssim \varepsilon^{-1}$ ，流动完全由黏性主导，宏观表现为达西律；一旦超越此阈值，惯性不可忽略，需更高阶模型（如 Brinkman 或 Forchheimer）。这一临界尺度的精确刻画，为工业中聚合物注塑、泥浆过滤等过程的模型简化提供了严格的数学判据。

6. 受点源驱动的黏性Burgers方程：弱极限的显式表达式

Viscous Burgers equation driven by point source: a formula for the weak limit

作者: Smritikana Pal, Manas R. Sahoo

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08496v1

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

详细分析

中文标题: 受点源驱动的黏性Burgers方程：弱极限的显式表达式

研究对象和背景

本文研究的核心对象是带狄拉克点源项的黏性Burgers方程，其形式为：

\begin{cases} u_t + \left(\frac{u^2}{2}\right)_x = \epsilon u_{xx} + \delta(x)\delta(t), & x \in \mathbb{R},\, t > 0, \\ u(x,0) = u_0(x), & x \in \mathbb{R}, \end{cases}

其中 $\epsilon > 0$ 是黏性系数， $\delta(x)\delta(t)$ 表示位于原点 $$(0,0)$$ 的时空点源（在原文中简记为 $\delta$ ，上下文明确其为 $\delta(x)\delta(t)$ 或更常取为 $\delta(x)$ 作用于空间变量、时间上恒定；第1页(1.2)式与第2页(1.6)式表明实际考虑的是 $\delta(x)$ ，即仅空间点源，且不显含时间脉冲），而 $u_0 \in L^\infty(\mathbb{R})$ 是给定的有界初值函数。该方程属于非齐次守恒律范畴，可写成一般形式：

u_t + f(u)_x = \epsilon u_{xx} + \alpha(x,t)\delta,

其中通量函数 $f(u) = \frac{u^2}{2}$ 是严格凸的， $\alpha(x,t) \equiv 1$ （或局部为常数），而 $\delta$ 代表分布意义下的点源——这使得方程在 $$x=0$$ 处出现本质奇异性，解在该点附近无法具有经典正则性，甚至不能属于 $H^1_{\mathrm{loc}}$ 类。

值得注意的是，该模型并非抽象构造，而是具有明确物理动机：论文引言指出，此类方程可用于描述异质介质中的流体运动，例如石油工业中岩石孔隙结构突变区域的渗流问题；此时点源可解释为注入井或抽采井的瞬时/稳态流率贡献。从数学结构看，它等价于具有不连续通量的Hamilton-Jacobi方程的梯度形式，即若令 $$v_x = u$$ ，则 $$v$$ 满足：

v_t + \frac{1}{2}(v_x)^2 = \epsilon v_{xx} + H(x),

其中 $$H(x)$$ 为Heaviside函数（因 $\delta = H_x$ 分布意义下成立），从而将问题嵌入更广阔的不连续哈密顿系统框架中。

研究区域为全空间 $\mathbb{R}_x \times (0,\infty)_t$ ，无边界条件，但因点源存在，需在分布意义下理解解：即要求 $u_\epsilon \in C((0,\infty); H^1(\mathbb{R}) \cap L^\infty(\mathbb{R})) \cap C^\infty(\mathbb{R}\setminus\{0\} \times (0,\infty))$ ，满足如下弱公式（见原文(2.3)式）：

\int_0^\infty \int_{\mathbb{R}} \left[ u_\epsilon \phi_t + \frac{u_\epsilon^2}{2} \phi_x - \epsilon u_\epsilon \phi_x \right] dx dt + \int_0^\infty \phi(0,t)\,dt - \int_{\mathbb{R}} u_0(x)\phi(x,0)\,dx = 0,

对任意 $\phi \in C_c^\infty(\mathbb{R} \times [0,\infty))$ 成立。此处关键在于：测试函数 $\phi$ 在 $$x=0$$ 处无需特殊处理，但源项 $\delta(x)$ 直接贡献了 $\int_0^\infty \phi(0,t)\,dt$ 项，体现了点源的“边界型”影响——尽管定义域无边界，奇点本身扮演了虚拟边界角色。

与标准黏性Burgers方程（ $\epsilon u_{xx}$ 项存在但无 $\delta$ ）不同，该模型的零黏性极限 $\epsilon \to 0^+$ 不再收敛到经典的Lax-Oleinik熵解。原因在于：点源破坏了通量的全局光滑性与守恒结构，导致极限必须同时平衡输运效应、扩散效应与源强集中效应。正如摘要所强调：“弱极限关联于一个包含三类泛函的变分问题”，这在无源情形下是不存在的——标准Lax-Oleinik公式仅涉及单一分段极小化（对应特征线到达时间），而此处需在右半平面、左半平面及跨源路径三类候选中择优，反映出点源诱导的解结构分叉。

此外，初值 $$u_0$$ 仅被假设为有界可测函数（见(1.2)式脚注），未要求连续性、单调性或BV正则性，这显著弱于传统熵解理论所需的初始数据条件（如 $u_0 \in BV$ ）。因此，该工作本质上是在最粗略的函数空间中构建并刻画奇异扰动下的渐近行为，凸显其分析难度与普适性价值。

主要定理或主要结果

本文的核心成果由两个主定理构成，分别给出黏性解的显式表达式与零黏性弱极限的变分刻画。

定理 2.1（黏性正则化解的显式公式）
设 $\epsilon > 0$ ，初值 $u_0 \in L^\infty(\mathbb{R})$ 。则问题(1.2)的解 $u_\epsilon(x,t)$ 可表示为：

u_\epsilon(x,t) = -2\epsilon \frac{\partial_x \theta_\epsilon(x,t)}{\theta_\epsilon(x,t)},

其中辅助函数 $\theta_\epsilon$ 分区域定义为：

\theta_\epsilon(x,t) = \begin{cases} R_\epsilon(x,t), & x > 0, \\ L_\epsilon(x,t), & x < 0, \end{cases}

且 $R_\epsilon$ , $L_\epsilon$ 由以下公式给出：

R_\epsilon(x,t) = e^{-\frac{t}{2\epsilon}} \left[ \frac{1}{2\sqrt{\pi \epsilon t}} \int_0^\infty \left( e^{-\frac{(x-\xi)^2}{4\epsilon t}} - e^{-\frac{(x+\xi)^2}{4\epsilon t}} \right) \theta_0^\epsilon(\xi)\,d\xi + \int_0^t \left( e^{\frac{\tau}{2\epsilon}} g(\tau) \right)' \operatorname{erfc}\left( \frac{x}{2\sqrt{\epsilon(t-\tau)}} \right) d\tau + g(0)\operatorname{erfc}\left( \frac{x}{2\sqrt{\epsilon t}} \right) \right],

L_\epsilon(x,t) = \frac{1}{2\sqrt{\pi \epsilon t}} \int_0^\infty \left( e^{-\frac{(x-\xi)^2}{4\epsilon t}} - e^{-\frac{(x+\xi)^2}{4\epsilon t}} \right) \theta_0^\epsilon(-\xi)\,d\xi + \int_0^t g'(\tau) \operatorname{erfc}\left( \frac{-x}{2\sqrt{\epsilon(t-\tau)}} \right) d\tau + g(0)\operatorname{erfc}\left( \frac{-x}{2\sqrt{\epsilon t}} \right),

其中 $\theta_0^\epsilon(x) = \exp\left( -\frac{1}{2\epsilon} \int_0^x u_0(y)\,dy \right)$ ，而边界函数 $$g(t)$$ 由下述Volterra型积分方程确定：

g(t) = g(0) e^{-\frac{t}{4\epsilon}} I_0\left( -\frac{t}{4\epsilon} \right) + \frac{\epsilon^3}{2\sqrt{\pi}} \int_0^t \frac{1 - e^{-\tau}}{\tau^{3/2}} f(t-\tau)\,d\tau,

且 $$f(t)$$ 显式依赖于 $\theta_0^\epsilon$ 与 $$g(0)$$ 。该定理断言 $u_\epsilon$ 属于 $C((0,\infty); H^1(\mathbb{R}) \cap L^\infty(\mathbb{R}))$ ，并在 $\mathbb{R}\setminus\{0\}$ 上无穷次可微。

定理 3.1（零黏性弱极限的变分公式）
当 $\epsilon \to 0^+$ 时， $u_\epsilon$ 在分布意义下收敛于某函数 $$u$$ ，其在 $$x>0$$ 和 $$x<0$$ 区域分别由以下三重极小化问题给出：

对 a.e. $(x,t) \in (0,\infty) \times (0,\infty)$ ，

u(x,t) = \min\left\{ U_1^R(x,t),\, U_2^R(x,t),\, U_3^R(x,t) \right\} + t,

其中

U_1^R(x,t) = \min_{\substack{0 \leq \tau < t \\ 0 \leq u < \tau \\ \xi \geq 0}} \left\{ \frac{x^2}{2(t-\tau)} + \frac{\xi^2}{2(\tau - u)} - \int_0^\xi u_0(z)\,dz \right\},

U_2^R(x,t) = \min_{\substack{0 \leq \tau < t \\ \xi \geq 0}} \left\{ \frac{x^2}{2(t-\tau)} + \frac{\xi^2}{2\tau} + \int_0^{-\xi} u_0(z)\,dz \right\},

U_3^R(x,t) = \min_{\xi \geq 0} \left\{ \frac{(x-\xi)^2}{2t} + \int_0^\xi u_0(z)\,dz \right\}.

对 a.e. $(x,t) \in (-\infty,0) \times (0,\infty)$ ，

u(x,t) = \min\left\{ U_1^L(x,t),\, U_2^L(x,t),\, U_3^L(x,t) \right\},

其中

U_1^L(x,t) = \min_{\substack{0 \leq \tau < t \\ 0 \leq u < \tau \\ \xi \geq 0}} \left\{ \frac{x^2}{2(t-\tau)} + \frac{\xi^2}{2(\tau - u)} + \int_0^\xi u_0(z)\,dz \right\},

U_2^L(x,t) = \min_{\substack{0 \leq \tau < t \\ \xi \geq 0}} \left\{ \frac{x^2}{2(t-\tau)} + \frac{\xi^2}{2\tau} + \int_0^{-\xi} u_0(z)\,dz \right\},

U_3^L(x,t) = \min_{\xi \geq 0} \left\{ \frac{(x+\xi)^2}{2t} + \int_0^{-\xi} u_0(z)\,dz \right\}.

该极限 $$u$$ 是无黏性方程 $u_t + (u^2/2)_x = \delta(x)$ 的弱解，且与Adimurthi等人通过其他方法得到的解一致。

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

a) 核心分析学工具

本文的证明骨架建立在Hopf-Cole变换、热核表示、变分原理与分布收敛理论四大支柱之上。

首先，Hopf-Cole变换是处理黏性Burgers方程的基石工具。作者令：

\theta_\epsilon(x,t) = \exp\left( -\frac{1}{2\epsilon} \int_0^x u_\epsilon(y,t)\,dy \right),

将拟线性抛物型PDE $u_t + uu_x = \epsilon u_{xx} + \delta(x)$ 转化为线性扩散-吸收方程：

\theta_{\epsilon,t} = \epsilon \theta_{\epsilon,xx} - \frac{1}{2\epsilon} H(x) \theta_\epsilon,

其中 $$H(x)$$ 为Heaviside函数。此变换不仅线性化方程，更将点源项转化为空间依赖的衰减项，使问题可分解为左右半平面的子问题。

其次，在每个半平面内，作者采用热核方法求解。由于方程在 $$x>0$$ 和 $$x<0$$ 上退化为标准热方程（除边界耦合外），解可表为初值热核卷积与边界热核响应的叠加。关键技巧在于：将非齐次边界条件 $R_\epsilon(0,t) = L_\epsilon(0,t) = g(t)$ 的影响，通过误差函数 $\operatorname{erfc}$ 显式编码。例如，对右半平面，解表示为：

R_\epsilon(x,t) = \int_0^\infty \Phi_\epsilon(x-\xi,t)\theta_0^\epsilon(\xi)\,d\xi + \int_0^t \partial_t \left[ e^{\frac{\tau}{2\epsilon}} g(\tau) \right] \operatorname{erfc}\left( \frac{x}{2\sqrt{\epsilon(t-\tau)}} \right) d\tau + \cdots,

其中 $\Phi_\epsilon(z,t) = \frac{1}{\sqrt{4\pi \epsilon t}} e^{-z^2/(4\epsilon t)}$ 为热核。 $\operatorname{erfc}$ 的引入巧妙地将边界记忆效应（即历史 $g(\tau)$ 对当前状态的影响）转化为可积分形式，避免了传统格林函数法的奇异性困境。

第三，变分原理是通往零黏性极限的核心桥梁。作者观察到：当 $\epsilon \to 0$ ，Hopf-Cole变换中的对数项 $-2\epsilon \log \theta_\epsilon$ 的渐近行为由大偏差原理（Laplace方法）主导。具体而言，对右半平面， $-2\epsilon \log R_\epsilon(x,t)$ 的极限等价于以下变分问题的极小值：

\lim_{\epsilon \to 0} \left( -2\epsilon \log R_\epsilon(x,t) \right) = \inf_{\gamma} \left\{ \frac{1}{2} \int_0^t \dot{\gamma}(s)^2 ds + \int_0^t \mathbf{1}_{\{\gamma(s)=0\}} ds + \int_0^{\gamma(0)} u_0(z)\,dz \right\},

其中 $\gamma$ 为连接初值位置与 $$(x,t)$$ 的路径。此泛函的三项分别对应：动能（输运代价）、在源点驻留时间（源强贡献）、初值累积效应。通过路径分类（是否接触源点、何时接触、从哪侧接近），自然导出三类候选泛函 $$U_1^R, U_2^R, U_3^R$$ ，体现了解的多尺度竞争机制。

最后，分布收敛理论确保极限过程的合法性。作者利用 $u_\epsilon$ 的一致有界性（由 $L^\infty$ 初值与最大值原理保证）及 $H^1_{\mathrm{loc}}$ 估计，结合紧嵌入定理（ $H^1 \hookrightarrow L^2_{\mathrm{loc}}$ 紧），证得存在子列弱收敛；再借助变分公式的唯一性，确认整个序列收敛。所有估计均在不依赖 $$u_0$$ 高阶正则性的前提下完成。

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

首要难题是点源导致的解不连续性与边界层冲突。作者采用区域分解+匹配法：将 $\mathbb{R}$ 拆分为 $$x>0$$ 与 $$x<0$$ ，在每区内求解齐次热方程，再通过未知函数 $$g(t)$$ 匹配解在 $$x=0$$ 处的值。此 $$g(t)$$ 并非自由指定，而是由通量连续性条件 $R_{\epsilon,x}(0,t) = L_{\epsilon,x}(0,t)$ 确定——这源于原始方程在分布意义下对 $$u_x$$ 的要求。计算显示，该条件导出关于 $$g(t)$$ 的奇异积分方程，其核含 $(t-\tau)^{-3/2}$ 奇性。作者的突破在于：将该方程改写为卷积型Volterra方程，并运用Laplace变换将其转化为代数方程，最终通过逆Laplace变换与修正Bessel函数 $$I_0$$ 给出显式解。此技巧绕开了直接求解奇异积分方程的困难。

其次，零黏性极限中三类泛函的竞争分析极具挑战。作者并未直接证明逐点收敛，而是先建立 $-2\epsilon \log \theta_\epsilon$ 的上下界，分别由三类泛函控制，再利用Gamma-收敛思想证明下界可达性。例如，为证 $$U_3^R$$ 的极小化器可达，作者构造一条从初值点 $\xi$ 出发、以恒定速度 $\frac{x-\xi}{t}$ 运动至 $$(x,t)$$ 的直线路径，并验证其代价趋近 $$U_3^R$$ ；类似构造覆盖所有三类情形，确保极限恰为三者之最小。

与之前工作的比较

a) 突破与改进

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

第一，放宽了对初值正则性的要求。经典Lax-Oleinik理论（如Lax[10]、Oleinik[12,13]）要求 $u_0 \in BV(\mathbb{R})$ 或至少为Lipschitz连续，以保证特征线不交叉、熵解唯一。而本文仅假设 $u_0 \in L^\infty(\mathbb{R})$ ，甚至允许其为任意有界可测函数（如Dirichlet函数）。这意味着解的存在性与极限公式对初值的“病态性”完全鲁棒，极大拓展了理论适用范围。此前，Chung等[5]虽也处理点源，但依赖“一阶接触”技巧，隐含对初值光滑性的需求；而本文通过Hopf-Cole变换与热核表示，将正则性负担转移至线性方程，成功规避了非线性初值敏感性。

第二，首次给出带点源Burgers方程零黏性极限的完整显式变分公式**。先前工作如Karlsen-Risebro-Towers[8]仅证明极限存在并满足某种熵条件，但未给出构造性公式；Ostrov[14]聚焦于通量正则化极限，与本文的“固定奇异通量+黏性消失”路径不同。本文的三类泛函 $U_i^{R/L}$ 不仅提供计算工具，更揭示了点源如何重塑解的几何结构： $$U_3$$ 对应经典特征线（避开源点）， $$U_1$$ 对应“源点反射路径”（在源点停留后折返）， $$U_2$$ 对应“源点瞬时穿越路径”。这种基于路径类型的分类极小化，是无源情形所不具备的全新数学现象。

第三，解决了与Adimurthi等[1,4]解的一致性问题。Adimurthi等人通过动力系统方法与界面熵条件构造解，但未给出显式公式。本文不仅证明其解与vanishing viscosity limit重合，更通过变分公式为其提供了直观的物理诠释：三类泛函分别对应流体粒子的三种最优运动策略——直行、反弹、穿越源点。这一一致性确认了不同数学路径（PDE正则化 vs 动力系统）对同一物理问题的内在统一性，增强了理论的坚实性。

综上，本文通过精巧融合经典工具（Hopf-Cole, 热核）与现代分析（变分法, Gamma-收敛），在最弱初值假设下，为带点源的Burgers方程建立了首个完整的显式渐近理论，填补了该领域长期存在的空白。

7. 关于网络上压缩因子非常数的混合气体稳态解存在性的研究

On the Existence of Steady States for Blended Gas Flow with Non-Constant Compressibility Factor on Networks

作者: Simone Göttlich, Michael Schuster, Alena Ulke

类别: math.NA, math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08481v1

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

详细分析

中文标题: 关于网络上压缩因子非常数的混合气体稳态解存在性的研究

研究对象和背景

本文研究的是氢气–天然气混合气体在管道网络（gas pipeline networks）中传输的稳态（steady-state）流动问题。其物理建模基于等温欧拉方程组（isothermal Euler equations），该系统描述了可压缩、粘性、单相、二元混合气体在一维管道中的质量与动量守恒。与传统纯天然气模型不同，本文的关键创新在于：压缩因子 $$ Z $$ 不再是常数，而是依赖于混合物的组成（以质量分数 $\eta$ 表征）和压力 $$ p $$ 的函数 $Z(\eta,p)$ 。这一非线性、非局部、非理想的状态方程深刻改变了压力–密度关系，使得经典分析工具失效。

论文明确给出的基础控制方程为混合气体的质量守恒与动量平衡系统（见原文(1.2a)–(1.2c)及推广至一般二元混合物的(2.2a)–(2.2c)）。设 $\rho_{\text{NG}}$ , $\rho_{\text{H}_2}$ 分别为天然气与氢气的密度， $q_{\text{NG}}$ , $q_{\text{H}_2}$ 为其质量通量，则总密度与总通量定义为

\rho = \rho_{\text{NG}} + \rho_{\text{H}_2}, \quad q = q_{\text{NG}} + q_{\text{H}_2}.

在**稳态（ $\partial_t = 0$ ）**假设下，时间导数项消失，系统退化为常微分方程（ODE）系统。对任一管道边 $e \in \mathcal{E}_{\text{pipe}}$ ，长度为 $$ L_e > 0 $$ ，其上的稳态流满足：

\frac{d}{dx} q_e(x) = 0,

\frac{d}{dx} \big( \eta_e \rho_e(x) \big) = 0,

\frac{d}{dx} p_e(x) + \frac{\lambda_{\text{fr}}}{2D_e} \frac{|q_e| q_e}{\rho_e(x)} + \frac{q_e^2}{\rho_e^2(x)} \frac{d}{dx} \rho_e(x) = 0.

其中 $\eta_e \in [0,1]$ 是氢气的质量分数， $$ D_e > 0 $$ 为管径， $\lambda_{\text{fr}} > 0$ 为达西–魏斯巴赫摩擦系数。关键物理关系由广义状态方程给出：

p_e(x) = \frac{RT}{M(\eta_e)} Z\big(\eta_e, p_e(x)\big) \rho_e(x),

其中 $$ R $$ 为通用气体常数， $$ T > 0 $$ 为恒定温度， $M(\eta_e)$ 是混合物摩尔质量（由式(2.5)确定），而 $Z: [0,1] \times [0,\infty) \to (0,\infty)$ 满足如下基本正则性假设（原文Assumption 1 & 2）：

$Z(\eta,p)$ 关于 $$ p $$ 可微且可积；
$Z(\eta,p)$ 关于 $\eta$ 可微；
$Z(\eta,p) > 0$ 对所有 $(\eta,p) \in [0,1]\times[0,\infty)$ 成立。

该状态方程隐含地定义了 $\rho_e$ 关于 $$ p_e $$ 的非线性映射，从而将动量方程完全转化为仅含 $$ p_e $$ 的一阶ODE。经代数整理与链式法则应用（见原文(2.9)），最终得到压力损失方程：

\frac{d}{dx} \left[ \int_{p_0}^{p_e(x)} \frac{p'}{Z(\eta_e,p')} dp' + \frac{RT}{M(\eta_e)} q_e^2 \Big( \ln\big|Z(\eta_e,p_e(x))\big| - \ln p_e(x) \Big) \right] = -\frac{\lambda_{\text{fr}} RT}{2D_e M(\eta_e)} |q_e| q_e.

为简化记号，定义辅助函数：

F(\eta,q,p) := \int_{p_0}^{p} \frac{p'}{Z(\eta,p')} dp' + \frac{RT}{M(\eta)} q^2 \Big( \ln\big|Z(\eta,p)\big| - \ln p \Big).

则上述方程等价于：

\frac{d}{dx} F(\eta_e, q_e, p_e(x)) = -\frac{\lambda_{\text{fr}} RT}{2D_e M(\eta_e)} |q_e| q_e.

积分后即得隐式压力剖面公式（原文(2.11)）：

F(\eta_e, q_e, p_e(x)) = F(\eta_e, q_e, p_e(0)) - \frac{\lambda_{\text{fr}} RT}{2D_e M(\eta_e)} |q_e| q_e x.

该公式是全文分析的核心——它将压力沿管道的演化完全编码于函数 $$ F $$ 的反函数中。

研究区域为有向连通图 $\mathcal{G} = (\mathcal{V}, \mathcal{E})$ ，其中节点集 $\mathcal{V}$ 表示管网接点（供气/用气点），边集 $\mathcal{E} = \mathcal{E}_{\text{pipe}} \cup \mathcal{E}_{\text{comp}}$ 分为管道边与压缩机边。每条管道边 $$ e $$ 定义在空间区间 $$ [0,L_e] $$ 上；压缩机边 $e \in \mathcal{E}_{\text{comp}}$ 则被建模为零维节点间跳跃算子，满足压力提升关系 $p_{h(e)} = \gamma_e p_{f(e)}$ ，其中 $\gamma_e \ge 1$ 为压缩比。

边界条件体现为节点耦合条件，包含三类物理约束：

质量守恒：对每个节点 $v \in \mathcal{V}$ ，有
$\sum_{e \in \mathcal{E}} a(v,e) q_e = b_v,$

其中 $$ a(v,e) $$ 是节点–边关联矩阵（见(2.1)）， $b_v \in \mathbb{R}$ 为节点负载（负值表示供气，正值表示耗气）。
完美混合与组分守恒：节点处流入气体瞬时均匀混合，故流出各边的氢气质量分数相同，记为 $\eta_v$ ，并由加权平均给出：
$\eta_v = \frac{ \sum_{e \in \mathcal{E}(v)} \big(a(v,e) q_e\big)^+ \eta_e + \big(b_v^-\big) \zeta_v }{ \sum_{e \in \mathcal{E}(v)} \big(a(v,e) q_e\big)^+ + b_v^- },$

其中 $\zeta_v \in [0,1]$ 是供气节点 $v \in \mathcal{V}_{<0}$ 的输入氢气分数， $(\cdot)^+ = \max\{\cdot,0\}$ ， $(\cdot)^- = \max\{-\cdot,0\}$ 。
压力连续性：所有与节点 $$ v $$ 相连的边 $$ e $$ 在端点处压力相等，即
$p_e(x_e) = p_v, \quad \text{其中 } x_e = 0 \text{ if } v = f(e),\; x_e = L_e \text{ if } v = h(e).$

特别地，由于 $\eta_e$ 由上游节点决定（见(2.13)），即 $\eta_e = \eta_{f(e)}$ 若 $q_e \ge 0$ ，否则 $\eta_e = \eta_{h(e)}$ ，这导致 $\eta_e$ 关于流量方向具有符号依赖的不连续性，构成分析的主要难点。

函数空间方面，稳态解被理解为：对每条边 $$ e $$ ，压力 $p_e \in C^0([0,L_e])$ ，流量 $q_e \in \mathbb{R}$ ，节点压力 $p_v \in \mathbb{R}_+$ ，质量分数 $\eta_v \in [0,1]$ ，且整体满足上述耦合条件与物理约束。进一步，为保证适定性，作者引入亚音速流域（subsonic regime）作为自然解空间限制：要求马赫数 $\mathrm{Ma} < 1$ ，即
$\frac{q_e^2}{\rho_e^2} < \sigma_e^2 = \frac{RT}{M(\eta_e)} \left[ Z(\eta_e,p_e) - p_e \partial_p Z(\eta_e,p_e) \right],$
等价于
$p_e^2 - \frac{RT}{M(\eta_e)} q_e^2 \left[ Z(\eta_e,p_e) - p_e \partial_p Z(\eta_e,p_e) \right] > 0.$
该不等式定义了关键的亚音速区域 $\mathcal{D}_{\text{subsonic}} \subset [0,1] \times \mathbb{R} \times [0,\infty)$ ，并在后续证明中用于确保 $$ F $$ 的严格单调性（Lemma 2.2）。

主要定理或主要结果

本文的核心理论成果是以下存在性定理，其结论覆盖含压缩机的树状网络与单环网络，且适用于一大类非恒定压缩因子模型：

Theorem 3.1（原文Section 3 主要结论）
设 $\mathcal{G} = (\mathcal{V}, \mathcal{E})$ 为一个有向连通图，其中 $\mathcal{E} = \mathcal{E}_{\text{pipe}} \cup \mathcal{E}_{\text{comp}}$ ，且满足：
(i) 压缩因子 $Z(\eta,p)$ 满足 Assumption 1 与 Assumption 2，且 $Z(\eta,p) > 0$ ；
(ii) 所有管道参数 $$ D_e > 0 $$ , $\lambda_{\text{fr}} > 0$ , $$ L_e > 0 $$ 给定；
(iii) 所有压缩机比 $\gamma_e \ge 1$ 给定；
(iv) 节点负载 $b_v \in \mathbb{R}$ 满足 $\sum_{v \in \mathcal{V}} b_v = 0$ （全局质量守恒）；
(v) 输入氢气分数 $\zeta_v \in [0,1]$ 给定；
(vi) 网络结构为树状（acyclic）或恰含一个独立回路（single cycle）。

则存在一组 $(p_v)_{v \in \mathcal{V}} \in \mathbb{R}_+^{|\mathcal{V}|}$ , $(q_e)_{e \in \mathcal{E}} \in \mathbb{R}^{|\mathcal{E}|}$ , $(\eta_v)_{v \in \mathcal{V}} \in [0,1]^{|\mathcal{V}|}$ ，使得：

对每条管道边 $$ e $$ ，压力剖面 $$ p_e(x) $$ 由 Lemma 2.4 给出，即
$p_e(x) = g\big( \eta_e, q_e, F(\eta_e, q_e, p_{f(e)}) - \tfrac{\lambda_{\text{fr}} RT}{2D_e M(\eta_e)} |q_e| q_e x \big);$
对每条压缩机边 $$ e $$ ，满足 $p_{h(e)} = \gamma_e p_{f(e)}$ ；
所有节点耦合条件（质量守恒、完美混合、压力连续）均被满足；
且所有解位于亚音速区域 $\mathcal{D}_{\text{subsonic}}$ 内。

该定理未断言唯一性或稳定性，但其存在性结论突破了此前工作（如[35]）仅适用于恒定 $$ Z $$ 且无压缩机的局限。

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

a) 核心分析学工具

本文的证明框架建立在**隐函数定理（Implicit Function Theorem）与连续性论证（continuity argument）之上，辅以拓扑度理论（degree theory）与不动点方法（fixed-point method）**的思想内核。

首要工具是隐函数定理的应用（Corollary 2.3）。由 Lemma 2.2 已证：在亚音速区域内，函数

F(\eta,q,p) = \int_{p_0}^{p} \frac{p'}{Z(\eta,p')} dp' + \frac{RT}{M(\eta)} q^2 \big( \ln Z(\eta,p) - \ln p \big)

关于 $$ p $$ 严格单调递增且可微，故其关于 $$ p $$ 的偏导 $\partial_p F > 0$ 。由此，对任意固定 $(\eta,q)$ ，映射 $p \mapsto F(\eta,q,p)$ 是局部双射，其局部逆 $g(\eta,q,\cdot)$ 连续。此性质使压力剖面 $$ p_e(x) $$ 成为 $p_{f(e)}$ 与 $$ q_e $$ 的连续函数，为后续构造整体解提供基础。

其次，核心存在性证明采用参数延拓法（parameter continuation）与连续性论证。作者将整个网络的求解问题视为一个高维非线性方程组 $\Phi(\mathbf{u}) = 0$ ，其中未知向量 $\mathbf{u} = (p_v, q_e, \eta_v)$ 包含所有节点压力、边流量与节点组分。关键观察是：当压缩因子退化为常数（即理想气体情形）时， $$ F $$ 具有显式表达（如 Remark 2.1 所示），此时 $\Phi$ 是连续且良定义的，且已知存在解（[35]）。随后，作者构造一族插值压缩因子

Z_\theta(\eta,p) = \theta Z(\eta,p) + (1-\theta) Z_{\text{const}}, \quad \theta \in [0,1],

并考虑对应族方程 $\Phi_\theta(\mathbf{u}) = 0$ 。若能证明解集 $\{ \mathbf{u}_\theta : \Phi_\theta(\mathbf{u}_\theta) = 0 \}$ 在 $\theta \in [0,1]$ 上构成一条紧致、连通、无断裂的连续路径，则 $\theta = 1$ 时解必存在。为此，需验证两点：
(i) $\Phi_\theta$ 的雅可比矩阵在解附近非奇异（由隐函数定理与 $\partial_p F > 0$ 保证）；
(ii) 解在 $\theta \in [0,1]$ 上一致有界（即先验估计）。后者通过亚音速约束与物理参数的正性（ $D_e, \gamma_e, b_v$ 有界）实现：因 $$ p_v > 0 $$ 、 $$ q_e $$ 受质量守恒线性约束、 $\eta_v \in [0,1]$ ，故 $\mathbf{u}$ 落在某个有界闭集中。

对于含单环网络，作者进一步引入环流变量（circulation variable） $\psi$ ，将环上压力匹配条件转化为关于 $\psi$ 的标量方程，并利用中间值定理（Intermediate Value Theorem）证明其零点存在——此即典型的度为1的映射思想。

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

本文面临三大技术挑战，其应对策略极具启发性：

挑战一：压缩因子 $Z(\eta,p)$ 非显式可积，导致 $$ F $$ 无闭式解。
传统方法（如[35]）依赖 $$ F $$ 的显式反函数（如 $F(p) \propto p^2$ 或 $\ln p$ ）来解析追踪压力变化。本文彻底放弃显式计算，转而构建抽象函数 $$ F $$ 的结构性质：通过亚音速条件强制 $\partial_p F > 0$ ，从而仅凭单调性与连续性就足以定义反函数 $$ g $$ 及其连续依赖性。这是“几何代替代数”的典范。

挑战二：质量分数 $\eta_e$ 关于流量 $$ q_e $$ 的符号依赖不连续性。
由于 $\eta_e = \eta_{f(e)}$ 当 $q_e \ge 0$ ，而 $\eta_e = \eta_{h(e)}$ 当 $$ q_e < 0 $$ （式(2.13)）， $\eta_e$ 在 $$ q_e = 0 $$ 处发生跳跃，导致 $F(\eta_e,q_e,p)$ 在 $$ q_e = 0 $$ 处不可微，破坏标准不动点迭代的压缩性。作者的精妙技巧是：将流量方向视为预设参数，即先固定一个符号模式 $\boldsymbol{\sigma} = (\operatorname{sgn} q_e)_{e \in \mathcal{E}}$ ，在每个符号象限内 $\eta_e$ 成为固定节点变量的线性组合，从而使 $\Phi$ 在该开集内光滑；再通过有限枚举所有可能的 $\boldsymbol{\sigma}$ （因网络有限，符号组合数有限），并证明至少存在一个 $\boldsymbol{\sigma}$ 使得对应方程有解。此即“分片光滑+穷举论证”。

挑战三：压缩机引入的非线性跳跃边界条件。
压缩机边不满足微分方程，而施加离散跳跃 $p_{h(e)} = \gamma_e p_{f(e)}$ ，破坏了纯管道网络的单调压力衰减结构。作者将其嵌入节点压力连续性条件：对压缩机边 $$ e $$ ，直接设定 $p_{h(e)}$ 与 $p_{f(e)}$ 的比例关系，将其视为对节点压力向量 $$ (p_v) $$ 的线性约束。结合质量守恒的线性方程，整个系统仍保持足够自由度，且 $\gamma_e \ge 1$ 的保序性确保压力始终为正，维持亚音速区域的封闭性。

与之前工作的比较

a) 突破与改进

本文相较既有文献（特别是[35]、[7]、[4]）实现了三重实质性突破：

第一，模型普适性的根本拓展。
先前工作[35]仅处理恒定压缩因子 $Z \equiv \text{const}$ ，其本质是理想气体律，导致压力–密度呈线性关系， $$ F(p) $$ 可显式积分。本文首次在任意满足可微性与正性假设的 $Z(\eta,p)$ 下建立存在性，涵盖线性模型（如AGA公式）、二次模型（Papay公式）乃至更一般的非多项式模型（如对应真实实验数据的拟合函数）。这意味着结论不再依赖于工程近似，而扎根于热力学基本原理。

第二，网络拓扑复杂性的显著提升。
[35] 仅针对无压缩机的纯输运网络（transport networks），其图结构为森林（forest），无反馈回路。本文将结论推广至含可控元件（compressor stations）的网络，并处理单环结构。压缩机的引入意味着压力可升可降，打破了传统“压力沿流向单调递减”的单调性，迫使作者发展新的全局耦合分析技术，而非依赖路径积分的逐段估计。

第三，数学工具的范式转换。
以往研究（如[7],[4]）侧重于数值验证，通过仿真展示不同 $$ Z $$ 模型对压力分布的影响（如图1），但未提供理论保障。本文则完成从“现象观察”到“严格证明”的跃迁，其核心贡献在于：将物理约束（亚音速性）转化为数学条件（ $\partial_p F > 0$ ），进而驱动隐函数与连续性论证。这种“物理引导数学”的范式，为未来研究高压、超临界、多相混合流奠定了方法论基础。

综上，本文不仅填补了混合气体稳态理论的空白，更提供了一套可推广的分析框架，适用于任何由状态方程驱动的网络化传输系统。

8. 相对熵与玻尔兹曼方程的弱可压缩纳维–斯托克斯动力学

Relative entropy and slightly compressible Navier-Stokes dynamics of the Boltzmann equation

作者: Yuhan Chen, Ning Jiang

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.07957v1

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

详细分析

中文标题: 相对熵与玻尔兹曼方程的弱可压缩纳维–斯托克斯动力学

研究对象和背景

本文研究的核心对象是玻尔兹曼方程（Boltzmann equation）在三维环面 $\mathbb{T}^3$ 上的小马赫数（small Mach number）渐近极限，其目标是建立该动力学系统向弱可压缩纳维–斯托克斯–傅里叶系统（slightly compressible Navier-Stokes-Fourier system）的定量收敛性，并进一步通过低马赫数极限过渡到不可压纳维–斯托克斯–傅里叶系统（incompressible Navier-Stokes-Fourier, INSF）。整个分析框架立足于相对熵方法（relative entropy method），而非传统的Chapman-Enskog或Hilbert展开。

论文考虑的尺度化玻尔兹曼方程（rescaled Boltzmann equation）为：

\begin{cases} \displaystyle \frac{1}{\tau_\varepsilon} \partial_t f_\varepsilon + v \cdot \nabla_x f_\varepsilon = C(f_\varepsilon, f_\varepsilon), & x \in \mathbb{T}^3,\ v \in \mathbb{R}^3, \\ f_\varepsilon(0,x,v) = f^{\text{in}}_\varepsilon(x,v) \ge 0, \end{cases} \tag{BE}

其中 $\tau_\varepsilon = \varepsilon$ 是马赫数（Mach number）， $\varepsilon > 0$ 是克努森数（Knudsen number），即分子平均自由程与宏观特征长度之比。碰撞算子 $$C(f,g)$$ 定义为：

C(f,g)(v) = \iint_{\mathbb{R}^3 \times \mathbb{S}^2} B(v-v_1,\sigma) \big[ g_1 f_1' + f_1 g_1' - g f - f g \big] \, d\sigma \, dv_1,

其中 $f = f(v),\ f_1 = f(v_1),\ f_1' = f(v_1'),\ f' = f(v')$ ，且 $(v,v_1) \mapsto (v',v_1')$ 表示弹性碰撞的前向散射变换：

v' = \frac{v+v_1}{2} + \frac{|v-v_1|}{2}\sigma,\quad v_1' = \frac{v+v_1}{2} - \frac{|v-v_1|}{2}\sigma,\quad \sigma \in \mathbb{S}^2.

散射截面 $B(v-v_1,\sigma)$ 仅依赖于相对速度 $$v-v_1$$ 及其与单位向量 $\sigma$ 的夹角余弦 $\frac{v-v_1}{|v-v_1|}\cdot \sigma$ ，属标准硬球或硬势型假设（cutoff hard potentials）。

为研究流体极限，作者引入局部麦克斯韦分布（local Maxwellian）：

M(\rho,u,\theta)(v) = \frac{\rho}{(2\pi \theta)^{3/2}} \exp\left(-\frac{|v-u|^2}{2\theta}\right),

其中 $(\rho,u,\theta)$ 是宏观变量：密度、宏观速度、温度。当它们为时空函数 $$(t,x)$$ 时， $$M$$ 称为局部麦克斯韦分布；若为常数（如 $$M(1,0,1)$$ ），则称全局麦克斯韦分布。

关键设定是：令 $f_\varepsilon = M_\varepsilon + \varepsilon M_\varepsilon g_\varepsilon$ ，其中 $M_\varepsilon = M(\rho_\varepsilon,u_\varepsilon,\theta_\varepsilon)$ ，且宏观状态满足弱可压缩纳维–斯托克斯–傅里叶系统（记为 $\mathrm{(CNS)}_\varepsilon$ ）：

\begin{cases} \varepsilon \partial_t \rho_\varepsilon + \nabla_x \cdot (\rho_\varepsilon u_\varepsilon) = 0, \\ \varepsilon \partial_t u_\varepsilon + u_\varepsilon \cdot \nabla_x u_\varepsilon + \dfrac{1}{\rho_\varepsilon} \nabla_x (\rho_\varepsilon \theta_\varepsilon) = \varepsilon \nabla_x \cdot \big[\mu_\varepsilon(\rho_\varepsilon,\theta_\varepsilon)\, \sigma(u_\varepsilon)\big], \\ \varepsilon \partial_t \theta_\varepsilon + u_\varepsilon \cdot \nabla_x \theta_\varepsilon + \dfrac{2}{3} \theta_\varepsilon \nabla_x \cdot u_\varepsilon = \dfrac{\varepsilon}{\rho_\varepsilon} \Big[ \mu_\varepsilon(\rho_\varepsilon,\theta_\varepsilon)\, \sigma(u_\varepsilon):\sigma(u_\varepsilon) + \nabla_x \cdot \big(\kappa_\varepsilon(\rho_\varepsilon,\theta_\varepsilon) \nabla_x \theta_\varepsilon \big) \Big], \end{cases} \tag{CNS$_\varepsilon$}

其中应力张量 $\sigma(u) = \nabla_x u + (\nabla_x u)^\top - \frac{2}{3}(\nabla_x \cdot u) I$ ，粘性系数 $\mu_\varepsilon$ 与热导率 $\kappa_\varepsilon$ 由线性化玻尔兹曼算子 $$L_M$$ 的谱结构定义：

\mu(\rho,\theta) = \frac{1}{10} \int_{\mathbb{R}^3} A(V) : \widehat{A}(V)\, M\, dv,\quad \kappa(\rho,\theta) = \frac{2}{15} \int_{\mathbb{R}^3} B(V) \cdot \widehat{B}(V)\, M\, dv,

其中 $V = \frac{v-u}{\sqrt{\theta}}$ ， $A(V) = V \otimes V - \frac{1}{3}|V|^2 I$ ， $B(V) = V\big(\frac{|V|^2}{2} - \frac{5}{2}\big)$ ，而 $\widehat{A}, \widehat{B}$ 是 $$L_M$$ 在正交补空间 $\mathcal{N}^\perp$ 中的预解元解。

特别地，作者设宏观初值为小扰动形式：

(\rho_\varepsilon, u_\varepsilon, \theta_\varepsilon)(0,x) = (1+\varepsilon \tilde{\rho}^{\text{in}},\ \varepsilon \tilde{u}^{\text{in}},\ 1+\varepsilon \tilde{\theta}^{\text{in}}),

从而确保系统处于小马赫数 regime（ $Ma = \varepsilon$ ），对应Navier-Stokes 时间尺度（而非Euler尺度 $\tau_\varepsilon = 1$ ）。

研究区域为三维周期域 $\mathbb{T}^3$ ，即无边界效应，避免了边界层等复杂现象；初始数据要求非负、可积、具有有限质量、动量、能量及熵，且满足相对熵初值趋于零的条件：

\frac{1}{\varepsilon^2} H(f^{\text{in}}_\varepsilon \mid M^{\text{in}}_\varepsilon) \to 0 \quad \text{as } \varepsilon \to 0.

该条件等价于初始分布 $f^{\text{in}}_\varepsilon$ 在相对熵意义下以 $O(\varepsilon^2)$ 精度逼近局部麦克斯韦分布 $M^{\text{in}}_\varepsilon$ ，属于适度准备初值（moderately prepared data），介于完全准备（well-prepared）与未准备（ill-prepared）之间。

主要定理或主要结果

论文的核心结论由如下主定理给出：

Theorem 1.1.
设 $f_\varepsilon = M_\varepsilon + \varepsilon M_\varepsilon g_\varepsilon$ 是方程 (BE) 的一族解，满足时间尺度 $\tau_\varepsilon = \varepsilon$ 及假设 (2.1)–(2.4)（即局部守恒律与相对熵不等式）；设 $(\rho_\varepsilon, u_\varepsilon, \theta_\varepsilon) = (1+\varepsilon \tilde{\rho}_\varepsilon, \varepsilon \tilde{u}_\varepsilon, 1+\varepsilon \tilde{\theta}_\varepsilon)$ 是 $\mathrm{(CNS)}_\varepsilon$ 的一族解。则对任意 $t \ge 0$ ，有如下一致先验估计：

\frac{1}{\varepsilon^2} H(f_\varepsilon \mid M_\varepsilon)(t) \le \frac{1}{\varepsilon^2} H(f^{\text{in}}_\varepsilon \mid M^{\text{in}}_\varepsilon) + O(\varepsilon), \tag{1.4}

\int_0^t \int_{\mathbb{T}^3} \left[ \frac{1}{\varepsilon^4} D(f_\varepsilon) - \frac{1}{2} \mu_\varepsilon \sigma(u^b_\varepsilon) : \sigma(u^b_\varepsilon) - \frac{5}{2} \kappa_\varepsilon |\nabla_x \theta^b_\varepsilon|^2 \right] dx\,ds \lesssim \frac{1}{\varepsilon^2} H(f^{\text{in}}_\varepsilon \mid M^{\text{in}}_\varepsilon) + O(\varepsilon), \tag{1.5}

\int_0^t \int_{\mathbb{T}^3} \left[ \frac{1}{2} \mu_\varepsilon \sigma(\tilde{u}_\varepsilon - u^b_\varepsilon) : \sigma(\tilde{u}_\varepsilon - u^b_\varepsilon) + \frac{5}{2} \kappa_\varepsilon |\nabla_x \tilde{\theta}_\varepsilon - \nabla_x \theta^b_\varepsilon|^2 \right] dx\,ds \lesssim \frac{1}{\varepsilon^2} H(f^{\text{in}}_\varepsilon \mid M^{\text{in}}_\varepsilon) + O(\varepsilon). \tag{1.6}

其中 $M_\varepsilon = M(\rho_\varepsilon,u_\varepsilon,\theta_\varepsilon)$ ， $(\rho^b_\varepsilon, u^b_\varepsilon, \theta^b_\varepsilon)$ 是 $g_\varepsilon$ 的矩：

\rho^b_\varepsilon = \langle g_\varepsilon \rangle,\quad u^b_\varepsilon = \langle v g_\varepsilon \rangle,\quad \theta^b_\varepsilon = \left\langle \frac{|v|^2 - 3}{3} g_\varepsilon \right\rangle,

且 $$D(f)$$ 是玻尔兹曼熵耗散泛函：

D(f) = \frac{1}{4} \iiint_{\mathbb{R}^3 \times \mathbb{R}^3 \times \mathbb{S}^2} (f'_1 f' - f_1 f) \log \frac{f'_1 f'}{f_1 f} \, d\sigma\, dv_1\, dv.

作为推论，在初值相对熵满足 $\frac{1}{\varepsilon^2} H(f^{\text{in}}_\varepsilon \mid M^{\text{in}}_\varepsilon) \to 0$ 时，得到如下强收敛结果：

相对熵趋于零： $\displaystyle \lim_{\varepsilon \to 0} \frac{1}{\varepsilon^2} H(f_\varepsilon \mid M_\varepsilon)(t) = 0$ ；
耗散项集中于流体部分： $\displaystyle \lim_{\varepsilon \to 0} \int_0^t \int_{\mathbb{T}^3} \left[ \frac{1}{\varepsilon^4} D(f_\varepsilon) - \frac{1}{2} \mu \sigma(u^b_\varepsilon) : \sigma(u^b_\varepsilon) - \frac{5}{2} \kappa |\nabla_x \theta^b_\varepsilon|^2 \right] dx\,ds = 0$ ；
宏观变量差趋于零： $\displaystyle \lim_{\varepsilon \to 0} \int_0^t \int_{\mathbb{T}^3} \left[ \mu |\nabla_x (\tilde{u}_\varepsilon - u^b_\varepsilon)|^2 + \kappa |\nabla_x (\tilde{\theta}_\varepsilon - \theta^b_\varepsilon)|^2 \right] dx\,ds = 0$ 。

这表明：在相对熵范数下， $f_\varepsilon$ 以 $O(\varepsilon)$ 速率逼近 $M_\varepsilon$ ，且其宏观矩 $(\rho^b_\varepsilon, u^b_\varepsilon, \theta^b_\varepsilon)$ 以强 $$L^2$$ 拓扑收敛于 $(\tilde{\rho}_\varepsilon, \tilde{u}_\varepsilon, \tilde{\theta}_\varepsilon)$ ，从而确立了玻尔兹曼方程到弱可压缩NSF系统的定量、稳定、可控制的渐近一致性。

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

a) 核心分析学工具

本文采用的主导工具是相对熵方法（relative entropy method），其本质是一种Lyapunov型能量估计，但作用于概率密度之间的信息距离，而非传统 $$L^2$$ 或 Sobolev 范数。定义如下：

H(f \mid g) = \iint_{\mathbb{T}^3 \times \mathbb{R}^3} \left( f \log \frac{f}{g} - f + g \right) dx\,dv.

该泛函天然具备非负性（Gibbs不等式）、凸性与耗散结构：对玻尔兹曼方程， $H(f \mid M)$ 随时间单调递减，且其耗散率 $-\frac{d}{dt} H(f \mid M)$ 控制着 $$D(f)$$ —— 即碰撞导致的微观无序增长。

论文的关键创新在于：将“靶函数” $$g$$ 设为依赖于 $\varepsilon$ 的局部麦克斯韦分布 $M_\varepsilon$ ，而非固定全局麦克斯韦分布（如 $$M(1,0,1)$$ ）。由此导出的调制相对熵（modulated relative entropy）：

\frac{1}{\varepsilon^2} H(f_\varepsilon \mid M_\varepsilon)

成为衡量 $f_\varepsilon$ 与动态演化流体解之间偏差的合适度量。其时间导数经细致计算后分解为三类项：

输运项（transport terms）：源于 $M_\varepsilon$ 的时空调变，需用宏观方程 $\mathrm{(CNS)}_\varepsilon$ 抵消；
耗散项（dissipation terms）：含 $D(f_\varepsilon)$ 及粘性/热传导耗散 $\mu \sigma(u^b):\sigma(u^b),\ \kappa |\nabla \theta^b|^2$ ；
余项（remainder terms）：高阶小量 $O(\varepsilon)$ ，需借助线性化算子 $$L_M$$ 的谱分解与投影技术控制。

为此，作者频繁使用Hilbert空间投影：将 $g_\varepsilon$ 分解为 $P_L g_\varepsilon + P_K g_\varepsilon$ ，其中 $$P_L$$ 投影至碰撞不变子空间 $\mathcal{N} = \operatorname{span}\{1,v,|v|^2\}$ ， $$P_K$$ 投影至正交补 $\mathcal{N}^\perp$ 。利用 $$L_M$$ 在 $\mathcal{N}^\perp$ 上的正定性（即存在 $\lambda > 0$ 使得 $\langle L_M h, h \rangle \ge \lambda \|h\|_{L^2(M)}^2$ ），可控制 $P_K g_\varepsilon$ 的大小。

另一核心工具是Grönwall不等式：由微分不等式

\frac{d}{dt} \left( \frac{1}{\varepsilon^2} H(f_\varepsilon \mid M_\varepsilon) \right) \le C \left( \frac{1}{\varepsilon^2} H(f_\varepsilon \mid M_\varepsilon) \right) + O(\varepsilon)

直接推出指数型控制，再结合初值小性，获得一致有界性。

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

非线性碰撞项的处理：
玻尔兹曼碰撞算子 $$C(f,f)$$ 具有双线性、非局部、非交换特性。本文避开直接估计 $$C$$ ，转而利用相对熵耗散恒等式：
$-\frac{d}{dt} H(f \mid M) = D(f) - \iint C(f,f) \log \frac{M}{f} \,dx\,dv,$
并将右侧第二项通过Taylor展开与矩匹配转化为宏观变量的偏微分表达式。关键一步是引理 2.1：精确展开 $v \cdot \nabla_x \log M_\varepsilon$ ，揭示其与NSF方程中对流项、压力梯度、粘性应力及热通量的对应关系：
$v \cdot \nabla_x \log M(\rho,u,\theta) = \underbrace{\cdots}_{\text{Euler项}} + \underbrace{A(V) : \nabla_x u + B(V) \cdot \frac{\nabla_x \theta}{\sqrt{\theta}}}_{\text{NS修正项}}.$
多尺度耦合的分离：
由于 $M_\varepsilon$ 本身含 $\varepsilon$ -参数，其对时间/空间的导数产生 $\varepsilon$ -阶项。作者通过精细的渐近展开（见第9页公式 (2.15)–(2.18)），将 $f_\varepsilon = M_\varepsilon(1+\varepsilon g_\varepsilon)$ 代入各项，逐阶提取 $O(1), O(\varepsilon), O(\varepsilon^2)$ 项，并利用矩恒等式（如 $\langle A(v) \rangle = 0$ , $\langle B(v) \rangle = 0$ ）消去虚假发散项。
余项的闭环控制：
所有 $O(\varepsilon)$ 余项最终归结为形如 $\varepsilon \langle \widehat{A}(v), Q(P_K g_\varepsilon, P_L g_\varepsilon) \rangle$ 的内积。借助线性化碰撞算子的再生性（即 $$Q(P_K g, P_L g)$$ 仍落于 $\mathcal{N}^\perp$ ），再用 Cauchy-Schwarz 与谱 Gap 不等式，将其控制为 $\varepsilon \|P_K g_\varepsilon\|_{L^2(M)}^2 \lesssim \varepsilon D(f_\varepsilon)$ ，从而被主耗散项吸收。
声波效应的规避：
对于环面上的未准备初值，经典方法会遭遇声波振荡破坏强收敛。本文通过两步分解策略绕过此障碍：先证 $f_\varepsilon \to M_\varepsilon$ （压缩NSF），再证 $M_\varepsilon \to M_{\text{incomp}}$ （低马赫极限）。相对熵框架天然兼容此分解，因 $H(f_\varepsilon \mid M_\varepsilon)$ 不敏感于快变声模。

与之前工作的比较

a) 突破与改进

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

第一，首次建立玻尔兹曼方程到可压缩NSF系统的定量相对熵收敛。
此前工作（如 Golse & Saint-Raymond [21,22]、Levermore & Masmoudi [35]）聚焦于紧性论证（compactness argument），即从Boltzmann解族中抽取子列，证明其弱- $\star$ 收敛至某个满足NSF弱解的极限。此类方法无法给出收敛速率，亦难以刻画误差结构。而本文通过相对熵方法，不仅证明收敛，更获得显式余项 $O(\varepsilon)$ 与耗散控制（式 (1.4)–(1.6)），实现首次定量刻画。

第二，摆脱对Hilbert/Chapman-Enskog展开的依赖。
传统方法（如 Guo [25–27]、Jiang & Xiong [30]）需构造具特定形式的近似解（如 $f_\varepsilon \approx M + \varepsilon M g^{(1)} + \varepsilon^2 M g^{(2)}$ ），并验证其满足方程至某阶精度。这要求初值严格匹配展开式，适用范围受限。本文则对任意满足守恒律与初值熵小性的经典解均适用，无需预先假设解的结构，显著拓宽了适用场景。

第三，统一处理“可压缩”与“不可压”极限的桥梁角色。
以往研究常割裂两类极限：BGL纲领处理不可压极限，Caflisch-Nishida 处理可压Euler极限。本文提出的 $\mathrm{(CNS)}_\varepsilon$ 系统作为中介模型，其自身即为Chapman-Enskog二阶近似，且当 $\varepsilon \to 0$ 时自动退化为INSF（通过低马赫极限）。图1.1所示的“三角关系”首次在相对熵框架下被严格量化，为理解不同尺度间的连续过渡提供了新范式。

综上，本文不仅是技术路线的更新，更是物理图像与数学语言的双重深化：它确认了相对熵不仅是热力学概念，更是连接微观动力学与宏观流体的普适性度量工具，为未来研究量子玻尔兹曼方程、活性物质输运等广义动力学系统的渐近理论开辟了新路径。

9. 具有中性嵌入特征值的剪切流的不稳定性

Instability of shear flows with neutral embedded eigenvalues

作者: Hui Li, Siqi Ren, Yuxi Wang, Guoqing Zhang

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.07807v1

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

详细分析

中文标题: 具有中性嵌入特征值的剪切流的不稳定性

研究对象和背景

本文研究的是二维不可压欧拉方程在环带与实轴乘积区域 $\mathbb{T} \times \mathbb{R}$ 上的线性稳定性问题，聚焦于一类单调剪切流 $$V_s = (b(y), 0)$$ 的扰动演化行为。该背景流是稳态解，其线性化系统导出经典的 Rayleigh 方程，其本质是研究算子谱结构对扰动增长性的决定作用。

原始控制方程为二维不可压欧拉方程：

\begin{cases} \partial_t V + V \cdot \nabla V + \nabla P = 0, \\ \nabla \cdot V = 0, \\ V(0,x,y) = V_{\text{in}}(x,y), \end{cases}

其中速度场 $$V = (V_1,V_2)$$ ，压力为 $$P$$ ，涡量定义为 $\Omega = \partial_x V_2 - \partial_y V_1$ 。在涡量形式下，方程简化为输运型方程：

\partial_t \Omega + V \cdot \nabla \Omega = 0.

对剪切流 $$V_s = (b(y),0)$$ 进行线性化，设扰动涡量为 $\omega(t,x,y)$ ，对应扰动流函数 $\psi$ 满足 $\Delta \psi = \omega$ ，则线性化方程为：

\partial_t \omega + b(y)\partial_x \omega - b''(y)\partial_x \psi = 0,

即

\partial_t \omega + ik R_{b,k}\, \omega = 0,\quad \omega(0,y) = \omega_{\text{in}}(y),

其中 $k \in \mathbb{Z}\setminus\{0\}$ 是流向波数， $R_{b,k}$ 是Rayleigh 算子：

R_{b,k} = b(y)\,\mathrm{Id} - b''(y)(\partial_y^2 - k^2)^{-1}.

该算子作用于 $L^2(\mathbb{R})$ （或加权空间），其谱结构决定了线性稳定性：若存在特征值 $$c$$ 满足 $\operatorname{Im} c > 0$ ，则系统存在指数增长模态，称为谱不稳定；若谱仅含实轴上的连续谱，则传统上认为系统线性稳定（如 Couette 流 $$b(y)=y$$ ）。然而，本文关注一个长期被忽视但物理意义重大的中间情形：Rayleigh 算子在实轴上具有嵌入特征值（embedded eigenvalue），即存在 $c_* \in \mathbb{R} \cap \operatorname{Ran}(b)$ ，使得 $(R_{b,k} - c_*)\phi = 0$ 在 $H^2(\mathbb{R})$ 中有非零解 $\phi$ 。此时系统既无指数增长，也不满足经典线性稳定结论——这种状态被称为中性态（neutral state）。

论文假设剪切流 $$b(y)$$ 满足严格单调性、高阶正则性及可积性：

单调性：存在 $c_m \in (0,1)$ ，使得 $0 < c_m \leq b'(y) \leq c_m^{-1}$ 对所有 $y \in \mathbb{R}$ 成立；
正则性与衰减： $b''(y) \in H^4(\mathbb{R})$ ；
嵌入特征值存在性：对某固定波数 $$k=1$$ ， $R_{b,1}$ 恰有一个实嵌入特征值 $$c_* = 0$$ ，且无其他点谱。

边界条件体现为解在无穷远处趋于零： $\phi(y) \to 0$ 当 $y \to \pm \infty$ ，这等价于要求 $\phi \in H^2(\mathbb{R})$ 。值得注意的是，由于 Rayleigh 算子是非局部的（含 $(\partial_y^2 - k^2)^{-1}$ ），其定义域需仔细处理；文中隐含使用 Sobolev 空间 $H^2(\mathbb{R})$ 或加权空间以确保紧性与自伴性缺失下的谱分析可行性。

典型例子包括双曲正切流 $b(y) = \tanh y$ ，其在 $$k=1$$ 模式下具有唯一嵌入特征值 $$c_* = 0$$ 。该流虽单调、光滑、有界，却因嵌入特征值导致谱密度函数出现奇异性，破坏标准衰减估计。

综上，本文的研究对象是一个定义在 $\mathbb{T} \times \mathbb{R}$ 上、关于单调剪切流的线性化欧拉系统，其核心数学对象是 Rayleigh 算子 $R_{b,k}$ 的谱理论，特别关注当其具有实轴上嵌入特征值时所引发的非模态不稳定性（non-modal instability），该现象源于算子的非正规性（non-normality），而非传统谱不稳定性。

主要定理或主要结果

论文建立两个关键定理，刻画中性剪切流的两类不稳定性：

定理 1.2（一般中性态的强不稳定性）
设背景剪切流 $$b(y)$$ 满足假设 1.1（单调、 $$H^4$$ 正则、存在单个嵌入特征值 $$c_*$$ ，几何重数为 1）。则对任意 $$M > 0$$ ，存在单位初值 $\omega_{\text{in}} \in L^\infty(\mathbb{R}) \cap L^2(\mathbb{R})$ （即 $\|\omega_{\text{in}}\|_{L^\infty} + \|\omega_{\text{in}}\|_{L^2} = 1$ ）及时间 $$T > 0$$ ，使得线性化方程 $\partial_t \omega + ik R_{b,k}\, \omega = 0$ 的解 $\omega(t,y)$ 满足：

\|\omega(T,\cdot)\|_{L^\infty(\mathbb{R})} \geq M,\quad \|\omega(T,\cdot)\|_{L^2(\mathbb{R})} \geq M.

该结论表明：即使无指数增长，系统仍可呈现任意大尺度的代数增长（transient growth），属于经典意义上的Lyapunov 不稳定性。

定理 1.3（多重嵌入特征值诱导的线性增长）
在定理 1.2 假设基础上，若嵌入特征值 $$c_*$$ 的代数重数大于 1（即为多重特征值），则存在初值 $\omega_{\text{in}} \in L^\infty(\mathbb{R}) \cap L^2(\mathbb{R})$ 及常数 $$C > 0$$ ，使得解满足：

\|\omega(t,\cdot)\|_{L^\infty(\mathbb{R})} \geq C t,\quad \|\omega(t,\cdot)\|_{L^2(\mathbb{R})} \geq C t,\quad \forall t \geq 0.

该结论揭示了一种更强的不稳定性机制：线性时间增长，其根源在于算子存在广义特征函数（associated function） $\eta$ 满足 $(R_{b,1} - c_*)\eta = \omega_*$ ，从而生成类似 Jordan 块的演化行为。

此外，论文还给出粘性情形的推论：

推论 1.6（小粘性下的不稳定性保持）
在相同假设下，对任意 $$M > 0$$ ，存在初值 $\omega_{\text{in}}$ （ $\|\omega_{\text{in}}\|_{L^2}=1$ ）、时间 $$T>0$$ 及粘性系数 $\nu_0 > 0$ ，使得对所有 $0 \leq \nu \leq \nu_0$ ，Orr-Sommerfeld 方程 $\partial_t \omega + O_{b,k,\nu}\,\omega = 0$ 的解满足：

\|\omega(T,\cdot)\|_{L^2(\mathbb{R})} \geq M.

这说明嵌入特征值导致的不稳定性在小粘性扰动下仍鲁棒存在，从而否定了“粘性可抑制此类不稳定性”的直观猜测。

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

a) 核心分析学工具

本文的证明框架深度融合了谱理论、奇异积分估计、渐近分析与对偶性论证，其技术内核可概括为以下四类工具：

Rayleigh 方程基本解构造与 Wronskian 分析
关键是构造 Rayleigh 方程 $(b(y)-c)(\partial_y^2 - 1)\phi - b''(y)\phi = 0$ 的两族指数衰减解 $\phi^\pm(y,c)$ ，分别满足 $\phi^\pm(y,c) \to 0$ 当 $y \to \pm \infty$ 。由此定义 Wronskian：
$W(c) = \int_{-\infty}^{+\infty} \frac{1}{\phi(y',c)^2}\, dy'.$
引理 2.2 指出： $$c$$ 是特征值当且仅当 $$W(c)=0$$ 。但由于 $$W(c)$$ 在实轴上非解析，作者引入修正 Wronskian $\widetilde{W}(c)$ （见式 (3.8)），其在复平面上解析，从而可应用复分析工具（如 Cauchy 积分公式、留数定理）精确计算极限 $\lim_{c_i \to 0^\pm} \frac{\widetilde{W}(c)}{c}$ ，这是提取投影系数 $P(\omega_{\text{in}})$ 的关键。
投影分解与谱表示公式
解 $\omega(t,y)$ 被分解为嵌入特征模部分 $\omega_1$ 和连续谱部分 $\omega_2$ ：
$\omega(t,y) = e^{-i c_* t} P(\omega_{\text{in}})\, \omega_*(y) + e^{-i R_{b,1} t} (\omega_{\text{in}} - P(\omega_{\text{in}})\, \omega_*(y)).$
其中投影系数形式为奇异主值积分：
$P(\omega_{\text{in}}) \sim \mathrm{P.V.} \int_{\mathbb{R}} \frac{\omega_{\text{in}}(z)}{z}\, dz.$
通过选取支撑在 $$[1/Z,1]$$ 上的非负初值 $\omega_{\text{in}}$ ，可使该奇异积分任意放大，而 $\|\omega_{\text{in}}\|_{L^2}$ 保持为 1，实现“小初值→大输出”。
对偶性论证与无粘阻尼估计
为控制连续谱部分 $\omega_2$ ，作者采用 Wei–Zhang–Zhao 发展的对偶方法：将流函数 $\Psi = \Delta^{-1} \omega_2$ 分解，并利用剪切流的混合效应证明：
$\|\Psi_2(t,\cdot)\|_{L^2} \lesssim \frac{1}{t}.$
结合特征模部分 $\Psi_1$ 的静态性（因 $$c_*=0$$ ），得出 $\|\Psi(t,\cdot)\|_{L^2} \approx |P(\omega_{\text{in}})|$ ，从而 $\|\omega(t,\cdot)\|_{L^2} \gtrsim |P(\omega_{\text{in}})|$ ，完成定理 1.2 的证明。
广义特征函数构造与 Jordan 型演化
当代数重数 $$>1$$ 时，作者证明存在关联函数 $\eta(y)$ 满足 $(R_{b,1} - c_*)\eta = \omega_*$ 。其构造依赖于 $\phi^\pm$ 关于参数 $$c$$ 的导数在 $$c=c_*$$ 处的匹配条件（见第 7 页）。由此，初值取 $\omega_{\text{in}} = \eta$ 时，解显式为：
$\omega(t,y) = t e^{-i c_* t} \omega_*(y) + i e^{-i c_* t} \eta(y),$
直接给出线性增长。此步骤本质上是将无限维算子的Jordan 结构（缺陷性）与有限维矩阵的广义特征向量理论类比。

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

克服 Wronskian 非解析性：通过引入“好导数” $\partial_G = \partial_y + \frac{\partial_{c_r}}{b'(y)}$ （式 (2.11)–(2.12)），将 $$1/(b(y)-c_r)^2$$ 的奇异性转化为可控的 $|y-y_c|^{-2}$ 型，再结合 $\phi_1$ 的精细估计（引理 2.5），确保 $\widetilde{W}(c)$ 的解析延拓可行。
处理奇异积分发散：利用初值 $\omega_{\text{in}}$ 的符号一致性与支撑集中性，将主值积分转化为带符号振荡可忽略的正部主导积分，实现 $|P(\omega_{\text{in}})| \gg \|\omega_{\text{in}}\|_{L^2}$ 。
统一处理连续谱演化：避免直接估计半群 $e^{-i R_{b,1} t}$ ，转而通过对偶变量 $\Psi$ 利用频散混合效应，将问题降为一维振荡积分估计，大幅降低技术难度。

与之前工作的比较

a) 突破与改进

本文的核心突破在于首次系统揭示并严格证明了中性剪切流固有的、由非正规性驱动的不稳定性机制，填补了从“谱稳定”到“谱不稳定”之间的理论空白。

假设条件的实质性放宽：以往关于剪切流稳定性的经典工作（如 Rayleigh 准则、Howard 半圆定理、Fjørtoft 改进）均聚焦于排除不稳定谱，而对嵌入特征值情形几乎未触及。本文不仅接受其存在，更将其作为研究起点，且仅需 $b'' \in H^4$ （远弱于早期工作要求的 $C^\infty$ 或解析性），极大拓展了适用范围。
结果强度的根本性提升：
- 对比 Couette 流（ $$b(y)=y$$ ）：已知其无嵌入特征值，满足无粘阻尼 $\|\omega(t)\|_{L^2} \lesssim \|\omega_{\text{in}}\|_{L^2}$ 。本文则证明：一旦引入嵌入特征值（如 $\tanh y$ ）， $$L^2$$ 范数可无界增长，颠覆了“单调流必稳定”的朴素认知。
- 对比 谱不稳定流（如 Poiseuille 流）：其不稳定性由 $\operatorname{Im} c > 0$ 导致指数增长，而本文的不稳定性是纯实谱下的代数/线性增长，物理机制完全不同（非模态 vs 模态）。
- 对比 Lin (2003)、Zillinger (2016) 等关于嵌入特征值的工作：他们主要研究其存在性或对邻近谱的影响，而本文首次定量刻画其动力学后果，并区分简单与多重情形，建立完整的不稳定性分类。
新现象的发现与命名：本文明确将此类不稳定性归因为 Rayleigh 算子的非正规性，并类比数值线性代数中的“伪谱（pseudospectrum）”概念，指出即使谱在实轴上，其伪谱可大幅向右延伸，导致巨大瞬态增长。这一观点为后续研究提供了全新范式。

综上，本文不仅是技术上的精进，更是概念上的革新：它确立了中性态剪切流是一种独立的、内在不稳定的流态，其不稳定性强度介于谱稳定与谱不稳定之间，且具有深刻的物理与数学内涵。

10. 不可压缩非牛顿流体适定弱解的一个加权正则性判据

A Weighted Regularity Criterion for Suitable Weak Solutions of Incompressible Non-Newtonian Fluids

作者: Jae-Myoung Kim

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07567v1

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

详细分析

中文标题: 不可压缩非牛顿流体适定弱解的一个加权正则性判据

研究对象和背景

本文研究的是三维欧氏空间 $\mathbb{R}^3$ 中不可压缩剪切变稀型非牛顿流体的运动行为，其数学模型由如下耦合偏微分方程组描述：

\begin{cases} \partial_t u - \nabla \cdot S(Du) + (u \cdot \nabla) u + \nabla \pi = 0, & \text{in } Q_T := \mathbb{R}^3 \times (0,T), \\ \operatorname{div} u = 0, & \text{in } Q_T, \\ u(x,0) = u_0(x), & x \in \mathbb{R}^3. \end{cases} \tag{1.1}

其中 $u : \mathbb{R}^3 \times (0,T) \to \mathbb{R}^3$ 表示速度场, $\pi : \mathbb{R}^3 \times (0,T) \to \mathbb{R}$ 是压力函数, 而 $S = (S_{ij})_{i,j=1}^3$ 是依赖于变形张量 $Du = \frac{1}{2}(\nabla u + \nabla u^\top)$ 的非线性应力张量。该模型属于经典的 “幂律流体”（power-law fluids） 类别，其本构关系取为：

S(Du) = (\mu_0 + \mu_1 |Du|^{q-2}) Du, \quad 1 < q < +\infty, \tag{1.3}

其中 $\mu_0 > 0$ , $\mu_1 > 0$ 为粘性系数常数。当 $$q = 2$$ 时退化为经典牛顿流体（即不可压缩 Navier–Stokes 方程），而 $q \ne 2$ 则刻画了非线性剪切依赖粘性：若 $$1 < q < 2$$ ，称作剪切变稀（shear-thinning）；若 $$q > 2$$ ，则为剪切增稠（shear-thickening）。本文聚焦于 $$q > 2$$ 的情形。

与经典 Navier–Stokes 方程不同，该系统中扩散项 $\nabla \cdot S(Du)$ 具有 $$q$$ -拉普拉斯结构，导致能量估计中自然出现 $$L^q$$ 型梯度范数，而非 $$L^2$$ 。这使得传统 Sobolev 嵌入、插值不等式及正则性理论需重新校准。特别地，由于 $$q$$ -增长非线性项在高频或奇点附近表现出更强的奇异行为，对解的局部结构分析提出了更高要求。

研究区域为全空间 $\mathbb{R}^3$ ，无边界条件；初始数据 $u_0 \in L^2(\mathbb{R}^3)$ 满足零散度约束 $\operatorname{div} u_0 = 0$ ，属于典型的适定弱解框架（suitable weak solution），即满足广义能量不等式（见 Definition 2.1）。该框架由 Caffarelli–Kohn–Nirenberg 对 Navier–Stokes 方程开创性工作所确立，并被推广至非牛顿情形（如文献 [3], [9]）。

值得注意的是，本文采用的测试函数为带反距离权重 $|x - x_0|^{-1}$ 的形式，该权重在 $$x = x_0$$ 处具有临界奇异性，不属于标准 Sobolev 空间 $W^{1,2}$ ，因此不能直接代入弱公式。作者通过引入径向截断函数 $\eta_\varepsilon$ 构造正则化测试函数 $\phi_\varepsilon(x,t) = \eta_\varepsilon(x) |x - x_0|^{-1} u(x,t)$ ，从而在严格意义下建立能量型不等式（见第4页脚注1）。这种处理方式体现了对加权空间中分布解理论的精细把握。

此外，文中频繁使用混合范数空间 $L^{r,s}(Q_T)$ ，定义为：

\|f\|_{L^{r,s}(Q_T)} := \left( \int_0^T \left( \int_{\mathbb{R}^3} |f(x,\tau)|^r dx \right)^{s/r} d\tau \right)^{1/s},

用于刻画时空方向上不同可积性的平衡，是连接局部正则性与全局加权控制的关键桥梁。

综上，本文的研究对象是一个具有强非线性耗散机制、定义在无界区域、初始数据仅具 $$L^2$$ 正则性、且要求满足加权能量不等式的不可压缩非牛顿流体系统；其核心挑战在于克服 $|x-x_0|^{-1}$ 权重带来的奇异性，同时协调 $$q$$ -增长结构与经典 Caffarelli–Kohn–Nirenberg 不等式的兼容性。

主要定理或主要结果

论文的核心结论以如下定理形式呈现：

Theorem 1.1
设 $$q > 2$$ 。假设初始速度 $u_0 \in H^1(\mathbb{R}^3)$ 且满足 $\operatorname{div} u_0 = 0$ ，并满足

\sup_{x_0 \in \mathbb{R}^3} \| |x - x_0|^{-1/2} u_0 \|_{L^2(\mathbb{R}^3)} < \infty.

若弱解 $$u$$ 满足以下两个加权可积性条件：

\sup_{x_0 \in \mathbb{R}^3} \int_0^T \| |x - x_0|^\beta \nabla u(\cdot,\tau) \|_{L^{r,s}(\mathbb{R}^3)} \, d\tau < \infty,

\sup_{x_0 \in \mathbb{R}^3} \int_0^T \| |x - x_0|^{-4/3} |\nabla u(\cdot,\tau)|^{q-1} \|_{L^{3/2}(\mathbb{R}^3)} \, d\tau < \infty,

其中参数满足：

\frac{3}{r} + \frac{5q - 6}{2s} = \frac{5q - 8}{2} + 1 - \beta, \quad 1 < s < \infty,

\frac{3}{5q - 8} + 1 - \beta < r < \infty, \quad -1 \le \beta < \frac{5q - 6}{2}.

则该弱解 $$u$$ 在时间区间 $$(0,T]$$ 上是半正则的（semi-regular），即其奇点集的抛物型 Hausdorff 测度为零，且对任意 $(x,t) \in \mathbb{R}^3 \times (0,T]$ ，存在邻域内 $u \in C^\infty$ 。

该定理本质上给出了一个加权梯度判据：只要速度梯度 $\nabla u$ 在某种带幂权 $|x-x_0|^\beta$ 的混合范数下可控，即可排除奇异点的形成。这是对经典 Serrin 型判据（如 $\|u\|_{L^\alpha_t L^\beta_x}$ 控制）在非牛顿背景下的重要加权推广。

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

a) 核心分析学工具

论文构建证明体系所依赖的三大支柱工具如下：

Caffarelli–Kohn–Nirenberg (CKN) 加权插值不等式
这是全文技术中枢。引理 2.1 给出如下关键估计：对任意光滑紧支函数 $$u$$ ，
$\| |x|^{-\gamma} u \|_{L^p} \le C \| |x|^{-\alpha} \nabla u \|_{L^2}^\theta \| |x|^{-\beta} u \|_{L^q}^{1-\theta},$
其中参数需满足相容性条件（2.2）–（2.3）。该不等式将加权 $$L^p$$ 范数控制为加权梯度 $$L^2$$ 范数与加权函数 $$L^q$$ 范数的插值，为将高阶非线性项（如 $$|u|^3$$ , $|\nabla u|^{q-1} u$ ）降阶至可处理形式提供了刚性框架。尤其在估计 $$I_3$$ （含压力项）时，它使 $\||x-x_0|^{-2/3} u\|_{L^3}$ 可进一步分解为 $\||x-x_0|^{-1/q} \nabla u\|_{L^q}$ 和 $\||x-x_0|^{-1/2} u\|_{L^2}$ 的乘积，从而嵌入最终的能量递推结构。
Stein 奇异积分估计
用于处理压力项 $\Phi = \pi$ 的加权 $L^\gamma$ 估计。由动量方程取散度得 Poisson 型压力方程：
$-\Delta \Phi = \operatorname{div}\operatorname{div}(u \otimes u + |Du|^{q-2} Du).$
Stein 在 [11] 中证明：对 $1 < \gamma < \infty$ 且 $-\frac{3}{\gamma} < \delta < \frac{3(\gamma-1)}{\gamma}$ ，有
$\| |x - x_0|^\delta \Phi \|_{L^\gamma} \lesssim \| |x - x_0|^\delta |u|^2 \|_{L^\gamma} + \| |x - x_0|^\delta |Du|^{q-1} \|_{L^\gamma}.$
此估计允许将压力的加权范数完全由速度及其梯度控制，避免了对压力单独做先验估计的困难。
$\varepsilon$ -正则性理论与广义能量不等式
基于定义 2.1 中的“适定弱解”概念，作者调用已知的 $\varepsilon$ -正则性定理（Proposition 2.1）：若对某点 $$(x,t)$$ ，
$\limsup_{r \to 0} r^{2p - 5} \iint_{Q_r(x,t)} |\nabla u|^p \, dx dt \le \varepsilon,$
则 $$(x,t)$$ 是正则点。本文通过加权能量估计（2.10）导出
$\sup_{x_0} \int_0^T \int_{\mathbb{R}^3} |x - x_0|^{-1} |\nabla u|^q \, dx dt < \infty,$
再结合 $B_r(x) \subset \mathbb{R}^3$ 的包含关系，推出
$\frac{1}{r^{5 - 2q}} \int_{t - r^2}^t \int_{B_r(x)} |\nabla u|^q \, dy d\tau \to 0 \quad \text{as } r \to 0,$
从而满足 $\varepsilon$ -条件，完成正则性提升。

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

奇异性测试函数的正则化：因 $|x - x_0|^{-1}$ 不属于 $W^{1,2}_{\mathrm{loc}}$ ，不能直接作为测试函数。作者构造截断函数 $\eta_\varepsilon$ ，令 $\phi_\varepsilon = \eta_\varepsilon |x - x_0|^{-1} u$ ，利用 $|\nabla \eta_\varepsilon| \lesssim \varepsilon^{-1}$ 及支撑集远离 $$x_0$$ 的性质，在极限 $\varepsilon \to 0$ 下恢复原始加权能量恒等式（2.4）——此为全文逻辑自洽性的基石。
非线性项的分层控制策略：对形如 $|Du|^{q-1} : \nabla(|x - x_0|^{-1}) u$ 的项，作者将其拆分为
$|Du|^{q-1} |x - x_0|^{-q/q_1} \cdot |x - x_0|^{-1 - 1/q_1} |u|,$
再用 Hölder 不等式配对，并借助 Young 不等式吸收高阶项（如 $\frac{1}{16}\int |\nabla u|^q |x - x_0|^{-1}$ ），保留低阶加权 $$L^2$$ 项（如 $\||x - x_0|^{-1} u\|_{L^2}^q$ ），最终实现闭环估计。
混合范数与参数自由度的协同设计：通过设定 $L^{r,s}$ 中的 $$r,s$$ 满足特定线性关系（见 Theorem 1.1 条件），使得插值后各项在 Gronwall 不等式中能统一衰减，避免因参数失配导致的指数爆炸。例如，选择 $\kappa = -\frac{1}{6+(3-q)s-2q}$ 保证了 $\||x-x_0|^\kappa u\|_{L^{\delta/(1-\theta_1)}}$ 可控于初始加权 $$L^2$$ 范数。

与之前工作的比较

a) 突破与改进

本文相较已有文献实现了三重实质性推进：

从牛顿到非牛顿的加权判据迁移：此前 Z. Yong 在 [12] 中建立了三维 Navier–Stokes 方程的加权正则性判据，但仅适用于 $$q = 2$$ 的线性扩散情形。本文首次将 CKN 不等式与非牛顿 $$q$$ -结构深度融合，给出适用于任意 $$q > 2$$ 的统一加权判据，填补了该方向的理论空白。
放宽对初始数据的正则性要求：经典结果（如 [9]）通常要求 $u_0 \in W^{2,2}$ 或更高；而本文仅需 $u_0 \in H^1$ 且满足加权 $$L^2$$ 条件 $\sup_{x_0} \||x-x_0|^{-1/2} u_0\|_{L^2} < \infty$ ，该条件比 $$L^2$$ 更弱（因权重衰减），却足以启动加权能量迭代。
提升判据的“局域适应性”：以往 Serrin 型判据（如 [1] 中针对 Ostwald–de Waele 流体的结果）依赖全局 $L^\alpha_t L^\beta_x$ 范数，无法识别解在某点附近的局部奇性。本文判据以 $\sup_{x_0 \in \mathbb{R}^3}$ 形式出现，天然具备点态敏感性，可精确刻画每个潜在奇点 $$x_0$$ 处的梯度爆破速率，从而实现真正意义上的局部正则性刻画。

综上，该工作不仅拓展了 CKN 不等式在非线性 PDE 中的应用疆界，更确立了一类新型加权正则性理论范式，为后续研究复杂流变流体（如 Bingham、Herschel–Bulkley 模型）提供了可复用的技术蓝图。

11. 关于二维定常紧支集欧拉流中常值涡度的研究

On two-dimensional steady compactly supported Euler flows with constant vorticity

作者: Changfeng Gui, Jun Wang, Wen Yang, Yong Zhang

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07407v1

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

详细分析

中文标题: 关于二维定常紧支集欧拉流中常值涡度的研究

研究对象和背景

本文研究的是二维定常、不可压缩、具有常值涡度（constant vorticity）的欧拉流，其数学模型由稳态欧拉方程在具有自由边界的紧支集区域上给出。核心物理图像是：流体被限制在一个环形（annular）或两相分层（two-phase）有界区域内，所有流线闭合，且涡度 $\gamma$ 为常数——这对应于物理中“均匀旋转”或“刚性旋转”类流动，如涡旋核、旋转液柱等。

论文聚焦三类具有不同边界结构与耦合机制的自由边界椭圆型问题，均以流函数 $\psi$ 为基本未知量。在二维不可压缩流中，速度场可表示为 $(u,v) = (\psi_y,-\psi_x)$ ，满足连续性 $\nabla \cdot (u,v) = 0$ ；而涡度 $\omega = \partial_x v - \partial_y u = -\Delta \psi$ ，故常值涡度条件即转化为 Poisson 方程：

\Delta \psi = -\gamma \quad \text{in the fluid domain}.

但需注意：原文中约定符号习惯为 $\Delta \psi = \gamma$ （见(1.2)式），即定义涡度为 $\omega = \Delta \psi$ ，因此我们统一采用该约定。以下三类问题构成全文研究骨架：

第一类：部分超定问题（Partially overdetermined problem）

对应系统 (1.2)，定义在环形区域 $\Omega \setminus D_0$ 上，其中 $D_0 = B_\lambda$ 是半径为 $\lambda \in (0,1)$ 的同心圆盘， $\Omega$ 是外边界待定的光滑单连通区域（含 $$D_0$$ ）。流函数 $\psi$ 满足：

\begin{cases} \Delta \psi = \gamma & \text{in } \Omega \setminus D_0, \\ \psi = 0 & \text{on } \partial \Omega, \\ \psi = 1 & \text{on } \partial D_0, \\ |\nabla \psi|^2 = Q > 0 & \text{on } \partial \Omega. \end{cases}

该系统称为Bernoulli 自由边界问题。前三个条件是标准 Dirichlet 边值，而第四个是非线性 Neumann 型条件（因 $|\nabla \psi|^2$ 依赖于梯度模长平方），使问题成为部分超定：在 $\partial \Omega$ 上同时指定 $\psi$ 和 $|\nabla \psi|$ 的大小（但未指定方向），自由度超出一般椭圆边值问题，因而解的存在性高度依赖区域几何。

第二类：两相超定问题（Two-phase overdetermined problem）

对应系统 (1.3)，定义在包含内核 $D \subset \Omega$ 的单连通区域 $\Omega$ 上， $$D$$ 亦为待定区域（初始取为 $B_\lambda$ ）。流场在内核 $$D$$ 与外部 $\Omega \setminus D$ 中具有不同常值涡度 $\gamma_1 \ne \gamma_2$ ，流函数分别为 $\psi_1, \psi_2$ ，并满足跨界面的连续性与通量匹配条件：

\begin{cases} \Delta \psi_1 = \gamma_1 & \text{in } D, \\ \Delta \psi_2 = \gamma_2 & \text{in } \Omega \setminus D, \\ \psi_1 = \psi_2 & \text{on } \partial D, \\ \frac{1}{\gamma_1} \partial_\nu \psi_1 = \frac{1}{\gamma_2} \partial_\nu \psi_2 & \text{on } \partial D, \\ \psi_2 = 0 & \text{on } \partial \Omega, \\ |\nabla \psi_2|^2 = Q > 0 & \text{on } \partial \Omega. \end{cases}

该模型模拟具有不同旋转特性的双层流体结构（如旋转核+剪切外壳），其超定性体现在：对同一边界 $\partial \Omega$ 同时施加 Dirichlet 和非线性 Neumann 条件；对内界面 $\partial D$ 则施加两个独立的匹配条件（位势连续 + 法向通量比例连续），属典型的两相超定椭圆问题。

第三类：完全超定问题（Truly overdetermined problem）

对应系统 (1.4)，同样定义在环形区域 $\Omega \setminus D$ 上，但内外边界均为自由边界，且各自配备独立的 Neumann 条件：

\begin{cases} \Delta \psi = \gamma & \text{in } \Omega \setminus D, \\ \psi = 0 & \text{on } \partial \Omega, \\ \partial_\nu \psi = q_{\mathrm{out}} & \text{on } \partial \Omega, \\ \psi = 1 & \text{on } \partial D, \\ \partial_\nu \psi = q_{\mathrm{in}} & \text{on } \partial D. \end{cases}

此处 $\nu$ 表示外法向（故在 $\partial D$ 上 $\nu$ 指向 $$D$$ 外部，即朝向 $\Omega \setminus D$ 内部），因此该系统在每条自由边界上均同时指定 Dirichlet 与 Neumann 数据，属严格意义上的完全超定问题，经典理论中仅在球对称情形下存在解。

所有三类问题均设定在 $\mathbb{R}^2$ 中，区域为有界、单连通（或环状）、 $C^{2,\alpha}$ 光滑；解空间取为 Hölder 空间 $C^{2,\alpha}(\overline{\Omega \setminus D})$ ；自由边界通过法向扰动参数化：设基准域为同心圆环 $B_1 \setminus B_\lambda$ ，则扰动后外边界为

\partial \Omega_\eta = \{ x + \eta(x)\nu_0(x) \mid x \in \partial B_1 \}, \quad \eta \in C^{2,\alpha}(\partial B_1),

内边界为

\partial D_\xi = \{ x + \xi(x)\nu_0(x) \mid x \in \partial B_\lambda \}, \quad \xi \in C^{2,\alpha}(\partial B_\lambda).

这种形状微分（shape derivative）框架将自由边界问题转化为关于 $(\eta,\xi)$ 的非线性泛函方程，是全文分析的几何基础。

主要定理或主要结果

本文三大主定理分别对应三类问题，核心是建立局部分支曲线（local bifurcation curves），刻画当控制参数（常值涡度 $\gamma$ 或 $\gamma_2$ ）变化时，非平凡（nontrivial）解如何从径向对称基准解（standard annulus / concentric disks）中产生。

定理 1.1（部分超定问题的分支）

对任意 $\lambda \in (0,1)$ ，存在 $$s_0 > 0$$ 及一条解曲线 $(\gamma(s), \eta(s))$ ， $s \in (-s_0,s_0)$ ，满足：
(1) 当 $$s = 0$$ 时， $(\gamma(0), \eta(0)) = (\gamma_*, 0)$ ，其中
$\gamma_* = \frac{4}{\lambda^2 - 2\lambda^2 \ln \lambda - 1} < -4,$

且此时 $\Omega_{\eta(0)} \setminus B_\lambda = B_1 \setminus B_\lambda$ 为标准同心环；
(2) 当 $s \ne 0$ 充分小时， $\eta(s) = s \alpha \cos\theta + o(s)$ ， $\alpha \ne 0$ ，对应非环形区域
$\Omega_{\eta(s)} \setminus B_\lambda = \{ x \in \mathbb{R}^2 \mid \lambda < |x| < 1 + \eta(s) \}.$

此即从径向解出发的非平凡分支。

定理 1.3（两相超定问题的分支）

设 $D = B_\lambda$ ，对任意 $\lambda \in (0,1)$ 及 $\gamma_1 \in \mathbb{R}\setminus\{0\}$ ，存在 $$s_0 > 0$$ 及解曲线 $(\gamma_2(s), \eta(s))$ ，满足：
(1) $(\gamma_2(0), \eta(0)) = (\gamma_1, 0)$ ，基准解为同心圆盘 $B_\lambda \cup (B_1 \setminus B_\lambda) = B_1$ ；
(2) $s \ne 0$ 时， $\eta(s) = s \alpha \cos\theta + o(s)$ ，区域 $B_\lambda \cup (\Omega_{\eta(s)} \setminus B_\lambda)$ 非同心；
(3) 若额外要求 $\partial_{\nu\nu} \psi_2 = m$ （曲率相关条件）在 $\partial \Omega_\eta$ 上成立，则解存在当且仅当区域为同心圆盘——此为刚性结论（rigidity）。

定理 1.5（完全超定问题的分支）

对任意 $\lambda \in (0,1)$ ，存在 $$s_0 > 0$$ 及解曲线 $(\gamma(s), (\eta(s), \xi(s)))$ ，满足：
(1) $(\gamma(0), \eta(0), \xi(0)) = (\tilde{\gamma}, 0, 0)$ ，其中 $\tilde{\gamma}$ 取下列两值之一：
$\tilde{\gamma} = \frac{4}{\lambda^2 - 2\lambda^2 \ln \lambda - 1} \quad \text{或} \quad \tilde{\gamma} = \frac{4 - 4\lambda^{-2}}{2\lambda^2 \ln \lambda + \lambda - 2 + \lambda^2 - 2 - 2 \ln \lambda};$
(2) $s \ne 0$ 时， $(\eta(s), \xi(s)) = s(\alpha \cos\theta, \beta \cos\theta) + o(s)$ ， $\alpha \ne \beta$ ，区域非同心环。

此外，三类问题均配有对应的刚性定理（Remarks 1.2(iii), 1.4(ii), 1.6(ii)）及稳定性定理（Corollaries 1.1–1.3），表明：当 $\gamma$ 属于某正则区间（如 $\gamma \ge 0$ 且小），则唯一解必为径向对称；且在 Neumann 数据受小扰动 $\rho(\theta)$ 时，解仍存在并保持接近基准区域，且扰动响应呈线性主导：

\eta(\rho) = \sum_{k \ge 1} \frac{\tau_k}{\sigma_k} \cos(k\theta) + o(\|\rho\|_Y).

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

a) 核心分析学工具

本文构建解的核心策略是隐函数定理（Implicit Function Theorem）与 Crandall–Rabinowitz 分支理论（CR theory） 的协同运用。具体流程如下：

抽象化建模：将自由边界扰动 $\eta,\xi$ 视为未知函数，将原 PDE 系统转化为 Banach 空间间的非线性算子方程
$\mathcal{F}(\gamma; \eta, \xi) = 0,$

其中 $\mathcal{F}: \mathbb{R} \times X \to Y$ ， $X = C^{2,\alpha}(\partial B_1) \times C^{2,\alpha}(\partial B_\lambda)$ ， $Y = C^{0,\alpha}(\partial B_1) \times C^{0,\alpha}(\partial B_\lambda)$ 等。
线性化与 Fredholm 性质：计算 Fréchet 导数 $D_{(\eta,\xi)} \mathcal{F}(\gamma_*; 0,0)$ 。关键观察是：该线性算子在极坐标下可分离变量，其在 Fourier 模式 $\cos(k\theta), \sin(k\theta)$ 上的作用为乘以标量系数 $\sigma_k(\gamma)$ （见 (3.12), (4.13), (5.17)）。例如，对问题 (1.2)，有
$\sigma_k(\gamma) = \frac{k \left[ (1-\lambda^{2k}) \gamma - 4k \right]}{(1-\lambda^{2k}) \ln \lambda}.$

当 $\sigma_k(\gamma_*) = 0$ 且 $\partial_\gamma \sigma_k(\gamma_*) \ne 0$ 时，零空间维数为 1，满足 CR 分支定理条件。
形状导数（Shape Derivative）技术：区别于传统 Hanzawa 变换，作者直接对区域扰动求导，导出边界上的线性化算子。例如，对 Dirichlet 数据扰动，有
$\left.\frac{d}{ds}\psi_{\eta(s)}\right|_{s=0} = w \quad \text{满足} \quad \Delta w = 0,\; w|_{\partial B_1} = -\partial_\nu \psi_* \cdot \eta,$

再结合 Bernoulli 条件的线性化，最终得到封闭的边界积分方程。此法避免了坐标变换带来的高阶导数损失，是全文技术亮点。
谱分析与特征值计算：所有 $\sigma_k(\gamma)$ 的零点即为分支点。作者显式求解方程 $\sigma_k(\gamma) = 0$ ，证明对每个 $k \in \mathbb{N}_+$ ，均存在实数解 $\gamma_k^*$ ，且当 $$k=1$$ 时 $\gamma_1^* < 0$ ，当 $k \ge 2$ 时 $\gamma_k^* > 0$ （见 Remark 1.2(i)），揭示负涡度驱动低频失稳，正涡度激发高频模态的物理机制。

b) 技术性难题处理技巧

非线性边界条件的处理：对 $|\nabla \psi|^2 = Q$ 类条件，作者采用切向微分+法向展开：将 $\psi$ 在边界附近作 Taylor 展开 $\psi(x + t\nu) = \psi(x) + t \partial_\nu \psi(x) + \frac{t^2}{2} \partial_{\nu\nu}\psi(x) + \cdots$ ，再代入 $|\nabla \psi|^2 = (\partial_\nu \psi)^2 + |\nabla_\tau \psi|^2$ ，利用切向导数在边界上的消失性（因 $\psi$ 为常数），将非线性项降为 $\partial_\nu \psi$ 的二次式，从而线性化为 $\partial_\nu \psi \cdot \delta(\partial_\nu \psi)$ 。
两相界面匹配的奇异性控制：在问题 (1.3) 中， $\frac{1}{\gamma_i} \partial_\nu \psi_i$ 的连续性导致跨界面导数跳跃。作者引入加权 Hölder 空间 $C^{2,\alpha}_{\mathrm{even}}$ （偶函数空间），并利用 Green 函数对调和函数的表示，将界面条件转化为关于 $\psi_1,\psi_2$ 边界值的奇异积分方程，再借助 Schur 引理 控制 Calderón–Zygmund 奇异积分算子的范数。
刚性证明中的最大值原理与 Reichel 定理：对非负解，由最大值原理得 $0 < \psi < 1$ ，再引用 Reichel 关于超定椭圆问题的对称化定理（Theorem 6.3 in Appendix），断言若在边界上同时满足 Dirichlet 与 Neumann 条件，则区域必为球（圆）——这是连接 PDE 解性质与几何刚性的关键桥梁。

与之前工作的比较

a) 突破与改进

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

第一，研究对象的根本转向：从开流域（strip）到闭流域（compact support）。以往二维定常水波理论（如 Constantin–Strauss [8], Wahlén–Weber [55]）几乎全部聚焦于无限深或有限深的带状区域 $\mathbb{R} \times (0,d)$ ，其自由表面为图 $y = \eta(x)$ 。而本文首次系统研究紧支集、闭流线、自由内外边界的环形/盘形区域，填补了“旋转液滴”“涡旋环”等物理模型的严格数学空白。此类区域缺乏平移不变性，传统 Fourier 方法失效，迫使作者发展基于形状导数与球谐展开的新框架。

第二，超定问题的构造范式革新。此前关于椭圆超定问题（Serrin 问题、Schiffer 问题）的工作（如 [47,50,37]）多假设区域凸性或使用移动平面法，解的存在性依赖强几何先验。本文则反其道而行之：不预设对称性，而将对称性作为分支结果；通过精确计算线性化算子的谱，主动构造非对称解，并给出其存在的充要参数条件（ $\gamma = \gamma_k^*$ ）。这标志着超定问题研究从“唯一性论证”迈向“存在性构造”。

第三，参数选择与验证简化。多数分支工作将几何参数（如内半径 $\lambda$ ）作为分支变量（如 [1,22]），导致线性化算子含 $\lambda$ -导数，计算繁复。本文创新地选取物理参数——常值涡度 $\gamma$ 为分支变量，使得线性化算子 $\sigma_k(\gamma)$ 成为 $\gamma$ 的仿射函数，CR 条件 $\partial_\gamma \sigma_k \ne 0$ 自动满足，极大简化了验证流程。正如 Remark 1.6(iii) 所强调：“This choice significantly simplifies the subsequent verification”。

综上，本文不仅拓展了欧拉流与超定问题的适用范围，更提供了一套可推广至三维、其他涡度分布乃至可压缩流的自由边界分析新范式。

12. 不可压磁流体力学的一个全局稳定性结果

A global stability result for incompressible magnetohydrodynamics

作者: Livio Pizzocchero, Emanuele Tassi

类别: math.AP, math-ph

发布日期: 2026-02-06

arXiv ID: 2602.06933v1

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

详细分析

中文标题: 不可压磁流体力学的一个全局稳定性结果

研究对象和背景

本文研究的是齐次、不可压、无界（周期）区域上的磁流体动力学（Magnetohydrodynamics, MHD）方程组，其数学模型描述了导电流体在磁场作用下的运动行为。具体而言，该系统由速度场 $$u(x,t)$$ 与磁场 $$b(x,t)$$ 构成，二者均为定义在 $$d$$ 维环面 $\mathbb{T}^d := (\mathbb{R}/2\pi\mathbb{Z})^d$ （其中 $d \in \{2,3,\dots\}$ ）上的向量场，满足如下偏微分方程组：

\begin{cases} \displaystyle \frac{\partial u}{\partial t} = \nu \Delta u - (u \cdot \nabla) u + (b \cdot \nabla) b - \nabla \left(p + \frac{|b|^2}{2}\right), \\ \displaystyle \frac{\partial b}{\partial t} = \eta \Delta b - (u \cdot \nabla) b + (b \cdot \nabla) u, \\ \nabla \cdot u = 0, \quad \nabla \cdot b = 0. \end{cases} \tag{1.2}

这里 $\nu > 0$ 和 $\eta > 0$ 分别为正粘性系数（viscosity）和正电阻率（resistivity）， $$p = p(x,t)$$ 是压力项， $\Delta$ 是拉普拉斯算子， $\nabla$ 是梯度算子， $(u \cdot \nabla) u := \sum_{s=1}^d u_s \partial_s u$ 表示对流导数。该系统是经典 Navier–Stokes（NS）方程的自然推广：当 $b \equiv 0$ 时，第一式退化为 NS 方程。

为消除压力项并突出结构特性，作者采用Leray 投影算子 $$L$$ ，将方程投影到散度为零的向量场上，从而得到等价的无压演化系统：

\begin{cases} \displaystyle \frac{\partial u}{\partial t} = \nu \Delta u - L((u \cdot \nabla) u) + L((b \cdot \nabla) b), \\ \displaystyle \frac{\partial b}{\partial t} = \eta \Delta b - L((u \cdot \nabla) b) + L((b \cdot \nabla) u). \end{cases} \tag{1.4}

进一步地，通过引入复合变量 $v = (u,b) \in H^p_{\Sigma_0} \times H^p_{\Sigma_0} =: H^p_{\Sigma_0}$ （见下文），该系统可重写为统一形式：

\frac{d v}{d t} = A v + \mathcal{P}(v,v), \tag{2.29.a}

其中线性算子 $$A$$ 定义为 $A(u,b) := (\nu \Delta u, \eta \Delta b)$ ，非线性双线性映射 $\mathcal{P}$ 由

\mathcal{P}((u,b),(w,g)) := \big(P(u,w) - P(b,g),\; P(u,g) - P(b,w)\big) \tag{2.24}

给出，而 $P(\cdot,\cdot) := -L((\cdot \cdot \nabla)\cdot)$ 是 NS 框架下的基本双线性映射。

函数空间的选择极为关键：作者工作于光滑、零均值、散度自由的周期向量场空间，即 Fréchet 空间

H^\infty_{\Sigma_0} := \bigcap_{p \in \mathbb{R}} H^p_{\Sigma_0},

其中每个 Sobolev 型空间 $H^p_{\Sigma_0}$ 定义为

H^p_{\Sigma_0} := \left\{ \ell \in \mathcal{D}'(\mathbb{T}^d,\mathbb{R}^d) \,\middle|\, \nabla \cdot \ell = 0,\; \langle \ell \rangle = 0,\; \sum_{k \in \mathbb{Z}^d \setminus \{0\}} |k|^{2p} |\hat{\ell}_k|^2 < +\infty \right\}, \tag{2.2}

其范数为 $\|\ell\|_p := \left( \sum_{k \neq 0} |k|^{2p} |\hat{\ell}_k|^2 \right)^{1/2}$ ，对应 Fourier 系数 $\hat{\ell}_k$ 。特别地， $H^p_{\Sigma_0}$ 中的零均值与散度自由条件确保了 Leray 投影 $$L$$ 的良定义性与连续性；而 Fréchet 结构则允许对任意高阶 Sobolev 范数进行统一控制。

研究区域为环面 $\mathbb{T}^d$ ，而非欧氏空间 $\mathbb{R}^d$ 或有界域，这带来若干优势：避免边界效应、保证 Fourier 展开绝对收敛、便于构造 Beltrami 类初值，并天然满足零均值要求（因周期性）。此外，所有解均被假定为 $C^\infty$ 光滑且时间全局存在（即定义于 $[0,+\infty)$ ），符合物理中“强解”或“经典解”的标准。

综上，本文的研究对象是一个具有正耗散（ $\nu,\eta>0$ ）、周期边界、零均值、散度自由约束的非线性耦合抛物型系统，其核心问题在于：给定一个已知的全局衰减解 $$v(t)$$ （即 $\|v(t)\|_p \to 0$ 当 $t \to +\infty$ ），是否所有足够接近它的初始数据 $$u_0$$ 都能生成另一个全局衰减解 $$u(t)$$ ？若成立，则称该解具有全局稳定性（global stability）。这一问题本质上属于非线性动力系统的鲁棒性分析，旨在刻画解对初值扰动的长期容忍能力。

主要定理或主要结果

本文的核心成果是以下全局稳定性定理：

Theorem 5.1（全局稳定性结果）
设 $v_0 \in E^\infty_{\Sigma_0}$ （即 $$v_0$$ 生成一个全局衰减 MHD 解 $v \in C^\infty([0,+\infty), H^\infty_{\Sigma_0})$ ），记 $$v(t)$$ 为其对应的解。对任意实数 $p \in \mathbb{R}$ ，令 $J_p \geq 0$ 满足
$\int_0^{+\infty} \|v(t)\|_p \, dt \leq J_p. \tag{5.1}$

再任取 $$n > d/2 + 1$$ ，并定义
$\rho_n := \frac{\mu}{\widehat{G}_n} \exp\big( -\widehat{G}_n J_n - \widehat{K}_n J_{n+1} \big), \tag{5.2}$

其中 $\mu := \min(\nu,\eta)$ ， $\widehat{G}_n, \widehat{K}_n$ 是非线性映射 $\mathcal{P}$ 在 $H^n_{\Sigma_0}$ 上的定量常数（见 (2.25)(2.26)）。则有：
(i) 若 $u_0 \in H^\infty_{\Sigma_0}$ 满足 $\|u_0 - v_0\|_n < \rho_n$ ，则 $$u_0$$ 亦属于 $E^\infty_{\Sigma_0}$ ，即它生成一个全局衰减解 $$u(t)$$ ；
(ii) 进一步，设 $\delta_n := \|u_0 - v_0\|_n$ ，则对所有 $t \geq 0$ ，有如下显式估计：
$\|u(t) - v(t)\|_n \leq \frac{e^{\widehat{G}_n J_n + \widehat{K}_n J_{n+1}}}{1 - \delta_n / \rho_n} \, \delta_n \, e^{-\mu t}; \tag{5.4}$

对任意 $$p > n$$ ，还有
$\|u(t) - v(t)\|_p \leq \exp\!\left( \widehat{G}_p J_p + \widehat{K}_p J_{p+1} + \frac{\widehat{G}_{pn} \delta_n / \rho_n}{1 - \delta_n / \rho_n} \right) \delta_p \, e^{-\mu t},$

其中 $\delta_p := \|u_0 - v_0\|_p$ ；
(iii) 若更强地假设 $\delta_n \leq \rho_n / 2$ ，则上述估计可简化为线性依赖形式：
$\|u(t) - v(t)\|_n \leq 2 e^{\widehat{G}_n J_n + \widehat{K}_n J_{n+1}} \delta_n e^{-\mu t}, \quad \|u(t) - v(t)\|_p \leq e^{\widehat{G}_p J_p + \widehat{K}_p J_{p+1} + \widehat{G}_{pn}} \delta_p e^{-\mu t}. \tag{5.5}$

该定理首次在任意维数 $d \geq 2$ 、任意高阶 Sobolev 范数 $p \in \mathbb{R}$ 下，给出了完全显式、定量、指数衰减型的全局稳定性结论，且不依赖于初值的小性假设（仅需相对于某个固定阶 $$n$$ 的距离小于一个显式阈值 $\rho_n$ ）。

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

a) 核心分析学工具

本文证明的核心框架建立在作者前期发展的后验近似解理论（a posteriori theory of approximate solutions）之上，该理论源自对 Navier–Stokes 方程近似解误差控制的经典方法（如 [30, 32]），并被系统拓展至 MHD 情形 [31]。其本质是一种基于 Riccati 型微分不等式的迭代控制策略，而非传统不动点或紧性论证。

设 $$u_a(t)$$ 是 MHD Cauchy 问题的一个近似解（未必满足初始条件），其误差被量化为：

微分误差： $e(u_a)(t) := \frac{d u_a}{dt} - A u_a - \mathcal{P}(u_a,u_a)$ ，
初值误差： $$u_0 - u_a(0)$$ 。

\frac{d^+ R_n}{dt} \geq -\mu R_n + (\widehat{G}_n D_n + \widehat{K}_n D_{n+1}) R_n + \widehat{G}_n R_n^2 + \epsilon_n, \quad R_n(0) \geq \delta_n. \tag{3.3}

该不等式的关键在于：其右端同时含线性项 $-\mu R_n$ （来自耗散）、双线性项 $(\cdots) R_n$ （来自非线性对误差的放大）、二次项 $\widehat{G}_n R_n^2$ （来自非线性自耦合）及外源项 $\epsilon_n$ 。求解此 Riccati 方程（或其等式情形）即可获得误差上界 $$R_n(t)$$ ，进而推出 $$u(t)$$ 的存在时间下界及距离估计。

对于更高阶范数 $$p > n$$ ，控制不等式退化为线性微分不等式：

\frac{d^+ R_p}{dt} \geq -\mu R_p + (\widehat{G}_p D_p + \widehat{K}_p D_{p+1} + \widehat{G}_{pn} R_n) R_p + \epsilon_p, \tag{3.5}

其解由显式公式给出：

R_p(t) = e^{-\mu t + A_p(t)} \left( \delta_p + \int_0^t e^{\mu s - A_p(s)} \epsilon_p(s) \, ds \right), \quad A_p(t) := \int_0^t \big( \widehat{G}_p D_p + \widehat{K}_p D_{p+1} + \widehat{G}_{pn} R_n \big) \, ds. \tag{3.7}

因此，整个证明逻辑链条为：将精确解 $$v(t)$$ 视作近似解 $$u_a(t)$$ ，则其微分误差为零，初值误差即为 $\delta_n = \|u_0 - v_0\|_n$ ，增长估计 $D_p(t) = \|v(t)\|_p$ 由衰减性提供 → 代入控制不等式 → 解出 $$R_n(t), R_p(t)$$ → 得到全局存在性与定量估计。

该方法完全规避了传统能量法中难以处理的非线性高阶耦合，也无需紧性或弱收敛——它是一种纯 ODE 驱动的、确定性的、可计算的误差传播分析。

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

MHD 系统的非线性项 $\mathcal{P}(v,v)$ 同时包含速度-速度、磁场-磁场、速度-磁场三类耦合，比 NS 更复杂。作者通过以下技巧实现可控性：

结构类比与常数量化：利用 MHD 与 NS 在 Leray 投影下的深刻相似性，将 $\mathcal{P}$ 的估计归结为 NS 基本映射 $$P$$ 的估计。文献 [27–29, 31] 提供了对 $$P$$ 的完全显式 Sobolev 常数上界：
$\|P(\ell,w)\|_p \leq \frac{1}{2} K_{pn} \big( \|\ell\|_p \|w\|_{n+1} + \|\ell\|_n \|w\|_{p+1} \big), \quad |\langle P(\ell,w), w \rangle_p| \leq \frac{1}{2} G_{pn} \big( \|\ell\|_p \|w\|_n + \|\ell\|_n \|w\|_p \big) \|w\|_p.$
并证明 $\widehat{K}_{pn} = \sqrt{2} K_{pn},\; \widehat{G}_{pn} = \sqrt{2} G_{pn}$ （见 (2.27)），从而将 MHD 常数完全由 NS 常数量化。
零均值与散度自由的 Fourier 表征：在环面上， $H^p_{\Sigma_0}$ 的范数直接由 Fourier 系数加权和定义（见 (2.2)）。这使得所有算子（ $\Delta$ , $$L$$ , $\mathcal{P}$ ）的估计均可转化为对 $$|k|^p$$ 权重的代数操作，避免了实空间中乘积估计（如 Moser 型不等式）的模糊常数，确保了全显式性。
衰减初值的构造技术：为验证定理适用性，作者引入广义 Beltrami 初值对 $$(u_0,b_0)$$ ，其 Fourier 支集位于锥面 $k \cdot u_0(k) = k \cdot b_0(k) = 0$ 上，且满足 $k \times u_0(k) = \alpha(k) b_0(k),\; k \times b_0(k) = \beta(k) u_0(k)$ 。此类初值虽可任意大（ $\|u_0\|_n$ 无界），但其非线性相互作用被几何约束极大削弱，从而仍产生全局衰减解。这是突破“小初值”限制的关键构造。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了多维度实质性突破：

维度普适性：此前关于 MHD 全局稳定性的结果（如 [42] 在 $\mathbb{R}^3$ ，[6] 在 $\mathbb{T}^3$ ）均局限于三维。本文将结论推广至任意 $d \geq 2$ 维环面，且所有估计常数均明确依赖于 $$d$$ （通过 $K_{pn}, G_{pn}$ 的 $$d$$ -依赖性体现），这是首次实现的维度无关性。
正则性自由度：文献 [42, 6] 均使用整数阶 Sobolev 范数（如 $$H^1, H^r$$ ），而本文工作于全实数阶 $H^p_{\Sigma_0}$ ，并给出对任意 $p \in \mathbb{R}$ 的一致指数衰减估计（见 (5.4)(5.5)）。这意味着不仅解本身，其任意阶导数（甚至分数阶导数）均以相同速率衰减，远超现有工作中仅控制有限阶导数的结论。
完全显式与定量：以往结果（如 [33, 42]）通常只断言存在某个 $\delta > 0$ ，但未给出其表达式；[6] 虽给出指数衰减形式 $e^{-\theta \mu t}$ ，但 $\theta \in (0,1)$ 为抽象常数。本文的 $\rho_n$ （见 (5.2)）及所有估计常数均完全由物理参数 $\nu,\eta$ 、Sobolev 常数 $\widehat{G}_n,\widehat{K}_n$ 和参考解 $$v$$ 的积分范数 $$J_p$$ 显式构成，具备可计算性与物理可解释性。
初值大小的实质性放宽：虽然定理仍要求 $\|u_0 - v_0\|_n < \rho_n$ ，但 $\rho_n$ 的构造（含 $\exp(-\widehat{G}_n J_n)$ ）表明：若参考解 $$v(t)$$ 衰减越快（即 $$J_n$$ 越小），则允许的扰动范围 $\rho_n$ 反而越大。更关键的是，作者在 Section 6 构造的广义 Beltrami 初值可任意大（ $\|v_0\|_n$ 任意大），却仍满足 $v_0 \in E^\infty_{\Sigma_0}$ ，从而使得该稳定性结果适用于大初值背景下的扰动分析，突破了“小初值”范式的长期束缚。
方法论革新：相比 [6] 使用的传统能量法与 Gronwall 不等式，本文采用的后验近似解框架将稳定性问题彻底转化为 ODE 控制问题，不仅简化了证明，还自然导出了高阶范数估计，且无需对解的长时间行为做额外先验假设（如 [42] 中的 $\int_0^\infty \|\nabla v(t)\|_{L^2}^4 dt < \infty$ ），仅需最自然的衰减条件 $\|v(t)\|_n \to 0$ （等价于指数衰减，见 Proposition 4.2）。

综上，本文不仅在技术细节上全面超越前人，更在哲学层面将 MHD 全局稳定性从一个“存在性断言”提升为一个“可计算、可验证、可设计”的工程化理论。

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

1. 用于模拟人体皮肤体温调节的演化血管中Stokes-输运耦合系统的推导与分析

Derivation and analysis of a Stokes-transport system in evolving vessels modeling thermoregulation in human skin

作者: Kilian Hacker, Maria Neuss-Radu

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08788v1

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

摘要翻译

中文标题: 用于模拟人体皮肤体温调节的演化血管中Stokes-输运耦合系统的推导与分析

摘要翻译: 本文研究一个定义在演化区域上的Stokes流动与对流-扩散输运过程的耦合系统，其边界条件允许流体流入与流出。区域的演化由输运过程本身所驱动，从而构成一个完全耦合的非线性问题。本研究旨在建模人体皮肤中血液流动的热调控机制。为此，该模型纳入了温度依赖性的生化物质生成过程、继而引发的血管舒张与收缩效应，以及由此导致的对流换热系数变化。我们采用不动点方法证明了该耦合系统弱解的存在性与唯一性，该方法可有效处理其中由物理机制所诱发的强非线性耦合。

2. 全直线上的阻尼BBM方程在低正则性空间中的强迫振荡

Forced oscillation of a damped BBM equation posed on whole line in low regularity spaces

作者: Chun Ho Lau, Taige Wang

类别: math.AP, math.DS

发布日期: 2026-02-09

arXiv ID: 2602.08327v1

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

摘要翻译

中文标题: 全直线上的阻尼BBM方程在低正则性空间中的强迫振荡

摘要翻译: 本文在低正则性Sobolev空间 $H^\ell$ （其中 $\ell \in [0,1)$ ）中，建立了带强迫项与阻尼项的一维Benjamin–Bona–Mahony（BBM）方程Cauchy问题的时周期解的存在性与稳定性结果。我们借助I-能量方法所导出的先验估计，在 $H^\ell$ 空间中为线性化问题获得所需估计；进而将对流项视为线性问题的扰动项，从而求解原Cauchy问题。

3. 基于分支理论方法研究非局部Kuramoto–Sivashinsky方程

A bifurcation theory approach to the nonlocal Kuramoto-Sivashinsky equation

作者: Pablo Cubillos, Rafael Granero-Belinchón y Juan Carlos Sampedro

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.08107v1

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

摘要翻译

中文标题: 基于分支理论方法研究非局部Kuramoto–Sivashinsky方程

摘要翻译:
我们在一维环面 $\mathbb{T}$ 上研究非局部Kuramoto–Sivashinsky方程：

u_t + u u_x = \Lambda^{r} u - \varepsilon \Lambda^{s} u, \qquad x \in \mathbb{T},

其中 $\varepsilon > 0$ ， $$s > 1$$ ， $r \in [-1, s)$ 。首先，我们证明了该方程在初值属于 $H^{3}(\mathbb{T})$ 时的局部适定性与整体适定性。随后，我们分析其稳态问题，证明平凡解分支（即零解分支）在临界参数值 $\varepsilon_k = k^{r-s}$ （ $k \in \mathbb{N}$ ）处发生分支。利用Crandall–Rabinowitz分支定理，我们得到从每个分岔点 $(\varepsilon_k, 0)$ 出发的光滑局部非平凡平衡解曲线，并精确计算了分支方向。为研究这些分支的整体延拓性，我们导出一致先验估计，并应用基于Fitzpatrick–Pejsachowicz–Rabier度理论（适用于指标为零的Fredholm映射）的整体替代定理。特别地，针对由首个临界点 $\varepsilon_1 = 1$ 分岔出的连通分支，我们证明其在 $\varepsilon$ -轴上的投影包含区间 $(2^{r-s}, 1)$ ，从而在此参数范围内确立了非平凡稳态解的存在性。最后，我们辅以数值延拓结果，直观呈现分支图结构及相应解的剖面形态。

4. 由 $|D|^{\frac{1}{2}}\xi$ 驱动的广义分数阶奇异Burgers方程的适定性

Well-posedness of Generalized Fractional Singular Burgers equation driven by $|D|^{\frac{1}{2}}ξ$

作者: Shuolin Zhang, Zhaonan Luo, Zhaoyang Yin

类别: math.AP, math.PR

发布日期: 2026-02-07

arXiv ID: 2602.07492v1

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

摘要翻译

中文标题: 由 $|D|^{\frac{1}{2}}\xi$ 驱动的广义分数阶奇异Burgers方程的适定性

摘要翻译: 本文研究由算子 $|D|^{\frac{1}{2}}\xi$ 驱动的分数阶奇异Burgers方程的广义解。我们构建了一个刻画该广义解所满足方程的理论框架，并将此类方程命名为广义分数阶奇异Burgers方程（Generalized Fractional Singular Burgers equation, GFSB）；在此框架下，我们证明了GFSB的局部适定性。最后，我们证明：当参数 $\gamma > \frac{3}{2}$ 时，GFSB的解即为原分数阶奇异Burgers方程的广义解。

不相关论文（基本信息）

1. $$(s,p)$$ -泊松型方程解梯度的精确可积性

Sharp gradient integrability for $$(s,p)$$ -Poisson type equations

作者: Verena Bögelein, Frank Duzaar, Naian Liao, Kristian Moring

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08944v1

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

摘要: We prove local $W^{1,q}$ -regularity for weak solutions to fractional $$p$$ -Laplacian type equations with right-hand side $f\in L^r_{\mathrm{loc}}(Ω)$ . Assuming $$p>1$$ , $s\in(0,1)$ , and $$sp'>1$$ , solutions belong to $W^{1,q}_{\mathrm{loc}}(Ω)$ for the optimal exponent $$q=q(n,p,s,r)$$ . We obtain quantitative local gradient estimates involving nonlocal tail terms. The optimality of $$q$$ is confirmed by a counterexample.

2. 半球面上扩张型调和映射流的存在性

Existence of expanding harmonic map flows to hemispheres

作者: Xuanyu Li

类别: math.AP, math.DG

发布日期: 2026-02-09

arXiv ID: 2602.08932v1

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

摘要: We show the existence of non-trivial self-expanding harmonic map flows starting from non-energy-minimizing 0-homogeneous maps to a regular ball or a closed hemisphere. In particular, given a non-minimizing but stationary 0-homogeneous harmonic map $$u_0$$ to a closed hemisphere, we construct infinitely many different weak solutions to harmonic map flow starting from $$u_0$$ , all of which satisfy the parabolic monotonicity formula. This answers a question of Struwe.

3. 齐次拉梅系统的刚性

Rigidity of homogeneous Lamé systems

作者: Joonas Ilmavirta, Teemu Saksala, Lili Yan

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08860v1

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

摘要: In this short paper, we show that any Lamé system whose Dirichlet-to-Neumann map for the elastic wave equation agrees with the one arising from the homogeneous Lamé system must actually be homogeneous. We do not need to impose any assumptions for the Lamé coefficients that we aim to recover. We use the fact that the homogeneous system gives rise to a geometry that is both simple and admits a strictly convex foliation.

4. 带边界项的加权Hardy-Sobolev型不等式

Weighted Hardy-Sobolev type inequalities with boundary terms

作者: João Marcos do Ò, Marcelo Furtado, Everaldo Medeiros, Jesse Ratzkin

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08702v1

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

摘要: In this paper we establish a new class of weighted Hardy-Sobolev type inequalities under mild monotonicity assumptions on the weight function. As a consequence, we derive the corresponding weighted Sobolev and trace-type inequalities. These results play an important role in the analysis of elliptic problems with Neumann or Robin boundary conditions in unbounded domains.

5. 构造性条件归一化流

Constructive conditional normalizing flows

作者: Borjan Geshkovski, Domènec Ruiz-Balet

类别: math.OC, cs.LG, math.AP, math.PR

发布日期: 2026-02-09

arXiv ID: 2602.08606v1

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

摘要: Motivated by applications in conditional sampling, given a probability measure $$μ$$ and a diffeomorphism $$φ$$ , we consider the problem of simultaneously approximating $$φ$$ and the pushforward $φ_{\#}μ$ by means of the flow of a continuity equation whose velocity field is a perceptron neural network with piecewise constant weights. We provide an explicit construction based on a polar-like decomposition of the Lagrange interpolant of $$φ$$ . The latter involves a compressible component, given by the gradient of a particular convex function, which can be realized exactly, and an incompressible component, which – after approximating via permutations – can be implemented through shear flows intrinsic to the continuity equation. For more regular maps $$φ$$ – such as the Knöthe-Rosenblatt rearrangement – we provide an alternative, probabilistic construction inspired by the Maurey empirical method, in which the number of discontinuities in the weights doesn’t scale inversely with the ambient dimension.

6. 能量临界Hartree方程双泡解的构造

Construction of two-bubble solutions for the energy-critical Hartree equation

作者: Jacek Jendrej, Xuemei Li, Guixiang Xu

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08490v1

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

摘要: We construct a pure two-bubble solution for the focusing, energy-critical Hartree equation in space dimension $N \geq 7$ . The constructed solution is spherically symmetric, global in (at least) the negative time direction and asymptotically behaves as a superposition of two ground states (or bubbles) both centered at the origin, with the ratio of their length scales converging to $$0$$ and the phases of the two bubbles form the right angle. The main arguments are the modulation analysis, the bootstrap argument and the topological argument. The main novelty with respect to existing constructions of pure two-bubble solutions is the nonlocal interaction, which is more complex to analyze.

7. 一类Plateau逼近问题的极小值解的存在性与正则性

Existence and Regularity of Minimizers for a Plateau Approximation Problem

作者: Eve Machefert

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08476v1

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

摘要: In this paper, we study the functional introduced by the author in collaboration with Bonnivard, Bretin, and Lemenant, which is designed to approximate Plateau’s problem. We establish the existence of a minimizer and prove its H{ö}lder regularity. Our results may be viewed as a generalization to higher-dimensional surfaces of the one-dimensional work of Bonnivard, Lemenant, and Millot on the approximation of the Steiner problem.

8. 一维准周期介质中频率域内的波传播——极限吸收原理

Wave propagation in the frequency regime in one-dimensional quasiperiodic media -Limiting absorption principle

作者: Pierre Amenoagbadji, Sonia Fliss, Patrick Joly

类别: math.AP, math-ph

发布日期: 2026-02-09

arXiv ID: 2602.08442v1

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

摘要: We study the one-dimensional Helmholtz equation with (possibly perturbed) quasiperiodic coefficients. Quasiperiodic functions are the restriction of higher dimensional periodic functions along a certain (irrational) direction. In classical settings, for real-valued frequencies, this equation is generally not well-posed: existence of solutions in L 2 is not guaranteed and uniqueness in L $\infty$ may fail. This is a well-known difficulty of Helmholtz equations, but it has never been addressed in the quasiperiodic case. We tackle this issue by using the limiting absorption principle, which consists in adding some imaginary part (also called absorption) to the frequency in order to make the equation well-posed in L 2 , and then defining the physically relevant solution by making the absorption tend to zero. In previous work, we introduced a definition of the solution of the equation with absorption based on Dirichlet-to-Neumann (DtN) boundary conditions. This approach offers two key advantages: it facilitates the limiting process and has a direct numerical counterpart. In this work, we first explain why the DtN boundary conditions have to be replaced by Robin-to-Robin boundary conditions to make the absorption go to zero. We then prove, under technical assumptions on the frequency, that the limiting absorption principle holds and we propose a numerical method to compute the physical solution.

9. Heisenberg群上混合局部与非局部退化椭圆方程的 $C^{1,\alpha}$ 正则性

$C^{1,α}$ -regularity for Mixed Local and Nonlocal Degenerate Elliptic Equations in the Heisenberg Group

作者: Junli Zhang

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08398v1

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

摘要: The regularity theory for equations combining both local and nonlocal operators in sub-Riemannian geometries is a huge challenge. In this paper, we investigate the $C^{1,α}$ -regularity of weak solutions to mixed local and nonlocal degenerate elliptic equations on the Heisenberg group. We first derive a sophisticated iteration scheme of Morrey-type by leveraging horizontal difference quotient combined with the Nirenberg difference quotient and fractional Sobolev-type inequality on the Heisenberg group. Then, the Hölder continuity of the weak solutions is established by applying the local boundedness, the iteration scheme of Morrey-type, an iterative method and the Morrey inequality. Finally, we use the Hölder continuity in conjunction with Theorem 1.2 from Mukherjee and Zhong[18] to prove the $C^{1,α}$ -regularity of weak solutions.

10. Heisenberg群上一类非局部Sobolev不等式的余项估计与尖锐定量稳定性

Remainder terms and sharp quantitative stability for a nonlocal Sobolev inequality on the Heisenberg group

作者: Wenjing Chen, Zexi Wang

类别: math.AP

发布日期: 2026-02-09

arXiv ID: 2602.08375v1

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

摘要: In this paper, we study the following nonlocal Sobolev inequality on the Heisenberg group \begin{equation}\label{eq:HLS} S_{HL}(Q,μ) \left(\int_{\mathbb{H}^{n}}\int_{\mathbb{H}^{n}}\frac{|u(ξ)|^{Q^{\ast}μ}|u(η)|^{Q^{\ast}μ}}{|η^{-1}ξ|^μ}{d}ξ{d}η\right)^{\frac{1}{Q^{\ast}μ}}\leq \int{\mathbb{H}^{n}}|\nabla{\mathbb{H}}u|^{2}dξ,\quad \forall , u\in S^{1,2}(\mathbb{H}^{n}), \end{equation} where $$Q=2n+2$$ is the homogeneous dimension of the Heisenberg group $\mathbb{H}^{n}$ , $n\geq1$ , $μ\in(0,Q)$ , $Q^{\ast}_μ=\frac{2Q-μ}{Q-2}$ is the upper critical exponent in the sense of the Hardy-Littlewood-Sobolev inequality and the Folland-Stein-Sobolev inequality on the Heisenberg group, $S_{HL}(Q,μ)$ is the sharp constant of \eqref{eq:HLS}, and $S^{1,2}(\mathbb{H}^{n})$ is the Folland-Stein-Sobolev space. %of the nonlocal-Sobolev inequality. It is well-known that, up to a translation and suitable scaling, \begin{equation}\label{eq:abs} -Δ{\mathbb{H}} u=\left(\int_{\mathbb{H}^{n}}\frac{|u(η)| ^{Q^{\ast}_μ}}{|η^{-1}ξ|^μ}{d}η\right)|u|^{Q_μ^*-2}u,~~u\in S^{1,2}(\mathbb{H}^{n}) \end{equation} is the Euler-Lagrange equation corresponding to the associated minimization problem. On the one hand, we show the existence of a gradient-type remainder term for inequality \eqref{eq:HLS} when $Q\geq4$ , $μ\in (0,4]$ , and as a corollary, derive the existence of a remainder term in the weak $L^{\frac{Q}{Q-2}}$ -norm on bounded domains. On the other hand, we establish the quantitative stability of critical points for equation \eqref{eq:abs} in the multi-bubble case when $$Q=4$$ and $μ\in (2,4)$ .

11. $$L^p$$ Dirichlet 可解性的稳定性：关于区域的小双 Lipschitz 变换

Stability of $$L^p$$ Dirichlet solvability under small bi-Lipschitz transformations of domains

作者: Joseph Feneuil, Linhan Li, Jinping Zhuge

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.08115v1

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

摘要: We show that small bi-Lipschitz deformations of a Lipschitz domain (with possibly large Lipschitz constant) preserve the solvability of the Dirichlet problem for the Laplacian with boundary data in $$L^p$$ , for the same value of $$p>1$$ . As a consequence, for all $p\in(1,\infty)$ , we obtain the solvability of the $$L^p$$ Dirichlet problem for small Lipschitz perturbations of convex domains, thereby unifying two fundamentally different settings in which such results were previously known: convex and $$C^1$$ domains. The key ingredient and novelty of our approach is a construction of a change of variables based on a non-constant basis derived from the Green function, which encodes the geometry of the base domain.

12. 关于数学生物学中某类含两个核函数模型的适定性问题

On the well-posedness of a certain model with two kernels appearing in the mathematical biology

作者: Messoud Efendiev, Vitali Vougalter

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.08102v1

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

摘要: The work is devoted to establishing the global well-posedness in $W^{(1,2),2}(R\times R^{+})$ of the integro-differential problem involving the two nonlocal terms describing the diffusion and the production in the biological system in the presence of the transport term. Such model is relevant to the cell population dynamics in the Mathematical Biology. The proof is based on a fixed point technique.

13. 一维拟线性克莱因-戈登流解寿命的改进上界

Enhanced lifespan bounds for 1D quasilinear Klein-Gordon flows

作者: Hongjing Huang, Mihaela Ifrim, Daniel Tataru

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.08055v1

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

摘要: In this article we consider one-dimensional scalar quasilinear Klein–Gordon equations with general nonlinearities, on both $\mathbb{R}$ and $\mathbb{T}$ . By employing a refined modified-energy framework of Ifrim and Tataru, we investigate long time lifespan bounds for small data solutions. Our main result asserts that solutions with small initial data of size $$ε$$ persist on the improved cubic timescale $|t| \lesssim ε^{-2}$ and satisfy sharp cubic energy estimates throughout this interval. We also establish difference bounds on the same time scale. In the case of $\mathbb{R}$ , we are further able to use dispersion in order to extend the lifespan to $ε^{-4}$ . This generalizes earlier results obtained by Delort in the semilinear case.

14. 奇异Grushin方程的可观测性性质

Observability properties of the singular Grushin equation

作者: Roman Vanlaere

类别: math.AP, math.OC

发布日期: 2026-02-08

arXiv ID: 2602.08044v1

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

摘要: We study the observability properties of the Grushin equation with an inverse square potential, whose singularity occurs at the boundary of two-dimensional rectangular domains or in the interior of the domain in higher dimensions. In some specific configurations of the observation set, we obtain the exact minimal time of observability. The analysis we present relies on recent Carleman estimates obtained by K. Beauchard, J. Dardé, and S. Ervedoza. As a byproduct of these results, we observe, for the heat equation associated to the Laplace-Beltrami operator on almost-Riemannian manifolds, a dependence of the minimal time of observability on the dimension of the singularity.

15. 具有毛细型边界条件的分数阶平均曲率流古典解的存在性

Existence of the classical solution to the fractional mean curvature flow with capillary-type boundary conditions

作者: Linlin Fan, Peibiao Zhao

类别: math.DG, math.AP

发布日期: 2026-02-08

arXiv ID: 2602.07989v1

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

摘要: Wang, Weng and Xia[Math. Ann. 388 (2024), no. 2] studied a mean curvature type flow for the smooth, embedded capillary hypersurfaces with a constant contact angle $θ\in(0,π)$ and confirmed the existence of solutions by the standard PDE theory. In the present paper, we study a fractional mean curvature flow for $C^{1,1}$ -regular hypersurfaces with a capillary-type boundary condition and obtain the short time existence by the fixed point argument.

16. 具有非正则数据的输运-扩散方程及其在二阶Hamilton-Jacobi偏微分方程稳定性估计中的应用

Transport-diffusion equations with irregular data and applications to stability estimates for second-order Hamilton-Jacobi PDEs

作者: Gianmarco Giovannardi, Alessandro Goffi

类别: math.AP

发布日期: 2026-02-08

arXiv ID: 2602.07969v1

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

摘要: This paper studies quantitative uniqueness properties in $$L^p$$ spaces for Fokker-Planck and transport-diffusion equations under two new assumptions on their velocity field $$b=b(x,t)$$ . We first prove $$L^p$$ -stability estimates for advection-diffusion PDEs when $\mathrm{div}(b)\in L^r_t(L^q_x)$ with $r\in[2,\infty]$ and $q\in[n/2,\infty)$ satisfying the compatibility condition $n/(2q)+1/r\leq 1$ . We then prove a stability result in $L^\infty$ for solutions of viscous transport equations when $\mathrm{div}(b(t))$ fails to be integrable in time. We apply these properties to obtain new continuous dependence estimates for viscous Hamilton-Jacobi equations via integral methods. One of the main novelties in this latter setting is that the constants of the estimates are all explicit with respect to the data of the problem. These imply new uniqueness properties for diffusive Hamilton-Jacobi equations without relying on the theory of viscosity solutions.

17. 紧致曲面上的格林–瓦瑟斯坦不等式

Green–Wasserstein Inequality on Compact Surfaces

作者: Maja Gwozdz

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

发布日期: 2026-02-08

arXiv ID: 2602.07843v1

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

摘要: Let $$(M,g)$$ be a compact connected two-dimensional Riemannian manifold without boundary. In this note, we answer a question posed by Steinerberger: can one remove the $\sqrt{\log n}$ factor in the two-dimensional Green–Wasserstein inequality while keeping the unrenormalized off-diagonal Green term? We show that this is impossible on any compact connected surface: there is no inequality of the same form that holds uniformly over point sets with an $O(n^{-1/2})$ remainder for all $$n$$ . We argue by contradiction and combine a second-moment estimate for the random Green energy of i.i.d. samples with the semi-discrete random matching asymptotics of Ambrosio–Glaudo.

18. 利用随机几何图求解边界数据的凸包络

Finding the convex envelope of a boundary datum using random geometric graphs

作者: Aurelia Deshayes, Nicolás Frevenza, Alfredo Miranda, Julio D. Rossi

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07696v1

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

摘要: In this paper we approximate the convex envelope of a boundary datum inside a bounded domain in the Euclidean space. We work with a random graph that is obtained as random points with uniform distribution that are connected by proximity ( $x\sim y$ when $$|x-y|<r$$ ). On the graph we solve an equation (that approximate the first eigenvalue of the Hessian of a smooth function) with an exterior datum. Under appropriate assumptions on $$r$$ we show that the unique solution to the equation in the graph converges to the convex envelope of the boundary datum as the number of points goes to infinity.

19. 积分型Harnack估计与奇异分数阶扩散方程解的熄灭速率

Integral Harnack estimates and the rate of extinction of singular fractional diffusion

作者: Filippo M. Cassanello, Simone Ciani, Antonio Iannizzotto

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07647v1

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

摘要: We prove several integral Harnack-type inequalities for local weak solutions of parabolic equations with measurable and bounded coefficients, describing singular s-fractional p-Laplacian diffusion. Then we apply the aforementioned estimates to evaluate the decay rate of the local mass and supremum of the solutions as they approach a possible extinction time. Yet we show consistency of our general decay estimates by studying the extinction phenomenon for weak solutions of the Cauchy-Dirichlet problem, by means of an approximation procedure that carefully avoids the use of an integrable time derivative.

20. 非柱状区域中耦合退化方程组的纳什-斯塔克尔伯格可控性

Nash-Stackelberg controllability for coupled systems of degenerate equations in non-cylindrical domains

作者: Alfredo S. Gamboa, Juan Limaco, Luis P. Yapu

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07582v1

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

摘要: In this paper we investigate the Hierarchical null controllability of a coupled degenerate semilinear parabolic equation in domains which are moving in time. We show the local null controllability of the semilinear system using Liusternik’s inverse function theorem. Nevertheless, the main difficulty is to adapt a Carleman estimate for the controllability of the linearized otimality system, using a Carleman inequality for degenerate non-autonomous equation obtanied by the authors previously.

21. Dirichlet 问题作为 Poisson 问题的边界情形：一个精确的逼近结果

The Dirichlet problem as the boundary of the Poisson problem: A sharp approximation result

作者: Mihalis Mourgoglou, Bruno Poggi

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07560v1

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

摘要: On a bounded domain $Ω\subset\mathbb R^{n+1}$ , $n\geq2$ , satisfying the corkscrew condition and with Ahlfors regular boundary, we characterize the dual space to the space ${\bf N}_{2,p}$ of functions $$u$$ whose Kenig-Pipher modified non-tangential maximal operator $\mathcal N_2(u)$ lies in $L^p(\partialΩ)$ , $p\in(1,\infty)$ . We find that

({\bf N}_{2,p})^*={\bf C}_{2,p'}\oplus L^{p'}(\partialΩ),\qquad\text{and that}\qquad L^{p'}(\partialΩ)=\partial^{\operatorname{weak}-*}{\bf C}_{2,p'}\,/\,{\bf C}_{2,p'},

where ${\bf C}_{2,p'}$ is a certain $L^{p'}$ -Carleson space and $$p'$$ is the Hölder conjugate of $$p$$ . This answers a question considered by Hytönen and Rosén. Inspired by this result and the recently understood characterizations of the $$L^p$$ -solvability of the Dirichlet problem in terms of the Poisson problem by Mourgoglou, Poggi, and Tolsa, we show a novel approximation result: for an arbitrary elliptic operator $L=-\operatorname{div} A\nabla$ with a not necessarily symmetric matrix $$A$$ of real bounded measurable coefficients, the solution space to the Dirichlet problem with data in $L^p(\partialΩ)$

\left\{\begin{aligned}-\operatorname{div} A\nabla u&=0,\quad&\text{in }&Ω,\\u&=g,\quad&\text{on }&\partialΩ,\end{aligned}\right.

lies on the weak- $$*$$ boundary in ${\bf N}_{2,p}$ of the solution space to the Poisson problem

\left\{\begin{aligned}-\operatorname{div} A\nabla w&=-\operatorname{div} F,\qquad&\text{in }&Ω,\\ w&=0,\qquad&\text{on }&\partialΩ,\end{aligned}\right.

with $F\in{\bf C}_{2,p}$ , provided that the Dirichlet problem for $$L$$ with data in $L^p(\partialΩ)$ is solvable in $$Ω$$ . This approximation result is sharp and new even for the Laplacian and on the unit ball.

22. 双曲空间中周长约束下Robin拉普拉斯算子的精确估计

Sharp estimates for the Robin Laplacian under a perimeter constraint in hyperbolic space

作者: Daguang Chen, Shan Li

类别: math.DG, math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07510v1

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

摘要: In this paper, we establish a lower bound, in terms of the isoperimetric deficit, for the first eigenvalue of the Robin Laplacian with negative boundary parameter on horospherically convex bounded domains in hyperbolic space, which implies that the geodesic ball maximizes this eigenvalue among all such domains. Furthermore, we derive upper bounds for the first eigenvalue of the Robin Laplacian with positive boundary parameter on horospherically convex bounded domains in hyperbolic space.

23. 具有空间增长非线性的聚焦非均匀薛定谔方程的归一化驻波解

Normalized Standing Waves for the Focusing Inhomogeneous Schrödinger Equation with Spatially Growing Nonlinearity

作者: Mohamed Majdoub, Tarek Saanouni

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07505v1

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

摘要: We study the focusing inhomogeneous nonlinear Schrödinger equation

i\partial_t u + Δu = -|x|^b |u|^{p-1}u ,\quad (t,x)\in (0,\infty)\times\mathbb{R}^N,

with $$b>0$$ and $$p>1$$ . Due to the spatial growth of the nonlinearity, standard compactness arguments do not apply and new difficulties arise. We first characterize ground state standing waves via a variational approach on the Nehari manifold and we establish some sharp stability and instability properties. In the $$L^2$$ -subcritical regime, we prove the existence of normalized ground states by solving a constrained energy minimization problem in the radial energy space, and we show that the resulting set of minimizers is orbitally stable under the flow. In contrast, in the $$L^2$$ -critical and supercritical regimes, ground state standing waves are shown to be strongly unstable by finite-time blow-up. Our results extend classical stability and instability theory for nonlinear Schrödinger equations to the case of spatially growing inhomogeneous nonlinearities.

24. 一般椭圆型与抛物型边值问题在非正则区域上的先验估计

A priori estimates for general elliptic and parabolic boundary value problems over irregular domains

作者: Maria R. Lancia, Alejandro Vélez-Santiago

类别: math.AP

发布日期: 2026-02-07

arXiv ID: 2602.07485v1

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

摘要: We investigate the realization of a myriad of general local and nonlocal inhomogeneous elliptic and parabolic boundary value problems over classes of irregular regions. We present a unified approach in which either local or nonlocal Neumann, Robin, and Wentzell boundary value problems are treated simultaneously. We establish solvability and global regularity results for both the stationary and time-dependent heat equations governed by general differential operators with unbounded measurable coefficients and various boundary conditions at once, first on a general framework, and then by presenting concrete important examples of irregular domains, Wentzell-type boundary conditions, and nonlocal maps. As a consequence, we develop a priori estimates for multiple differential equations under various situations, which are tied to a large number of applications performed over real world regions, such heat transfer, electrical conductivity, stable-like processes (probability theory), diffusion of medical sprays in the bronchial trees, and oceanography (among many others).

25. 实Bers嵌入在实直线上的研究：Fisher-Rao线性化、Schwarzian曲率与散射坐标

Real Bers embedding on the line: Fisher-Rao linearization, Schwarzian curvature, and scattering coordinates

作者: Hy Lam

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

发布日期: 2026-02-07

arXiv ID: 2602.07373v1

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

摘要: We develop a real-analytic counterpart of the Bers embedding for the Fréchet Lie group $\Diff^{-\infty}(\R)$ of decay-controlled diffeomorphisms of the line, and establish its connection to $$L^p$$ Fisher-Rao geometry on densities. For $p\in[1,\infty)$ , the $$p$$ -root map $\varphi\mapsto p(\varphi'^{1/p}-1)$ isometrically linearizes the homogeneous $\dot W^{1,p}$ Finsler metric on $\Diff^{-\infty}(\R)$ , yielding explicit geodesics and a canonical flat connection whose Eulerian geodesic equation is the generalized Hunter-Saxton equation; for $p=\infty$ , logarithmic coordinates $\varphi\mapsto\log\varphi'$ provide a global isometry and the Schwarzian derivative emerges as the projective curvature. We construct a real Bers map $β^{-\infty}\colon\Diff^{-\infty}(\R)/\Aff(\R)\to W^{\infty,1}(\R)$ via this Schwarzian, prove it is a Fréchet-smooth injective immersion whose linearization admits a tame right inverse given by an explicit Volterra operator, and characterize its image through Sturm-Liouville spectral theory. We introduce an $$L^p$$ -Schwarzian family $$S_p$$ that interpolates between affine and projective cocycles, establish full asymptotic expansions as $p\to\infty$ in Fréchet and Orlicz-Sobolev scales, and extend the Bers embedding to Orlicz diffeomorphism groups. Through the Jacobian correspondence, these structures transfer to a manifold of densities asymptotic to Lebesgue measure, where the nonlinear Eulerian transport reduces to a pointwise Riccati law and the Schwarzian becomes the score curvature governing Fisher information. The compact-manifold $$L^p$$ Fisher-Rao linearization of Bauer, Bruveris, Harms, and Michor is recalled as a guiding framework.

26. 二阶退化双曲型方程的非齐次边值问题及其应用

Non-homogeneous boundary value problems for second-order degenerate hyperbolic equations and their application

作者: Donghui Yang, Jie Zhong

类别: math.AP

发布日期: 2026-02-06

arXiv ID: 2602.07271v1

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

摘要: We study second-order hyperbolic equations with degenerate elliptic operators and non-homogeneous Dirichlet boundary inputs. We establish existence and regularity of weak solutions in weighted Sobolev spaces under mild assumptions on the degenerate weight. A Dirichlet map is constructed for the degenerate elliptic operator, leading to a solution theory that extends classical approaches to the degenerate setting. In particular, we derive energy estimates and well-posedness for boundary inputs of low regularity (in appropriate trace spaces), even though the classical Dirichlet-to-Neumann framework is not directly applicable in the degenerate setting. As an application, we prove an approximate controllability criterion, which generalizes the Hilbert Uniqueness Method to degenerate wave equations. Our framework accommodates higher-dimensional degenerate waves, non-homogeneous boundary conditions, and weighted functional analysis. We also illustrate how our criterion connects to higher-dimensional Grushin equations and waves with single-point degeneracy, and we highlight the remaining unique continuation/observability issue as an open problem.

27. 耦合热弹性系统中模态逼近的稳定性与收敛性：理论与数值模拟

Stability and Convergence of Modal Approximations in Coupled Thermoelastic Systems: Theory and Simulation

作者: I. Essadeq, S. Nafiri, S. Benjelloun, A. E. Fettouh

类别: math.NA, math.AP

发布日期: 2026-02-06

arXiv ID: 2602.07224v1

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

摘要: In this work, we review and analyze both the theoretical and numerical aspects of strongly and weakly coupled thermoelastic systems. By employing spectral analysis techniques and establishing uniform resolvent estimates, we derive uniform polynomial decay rates for the associated semigroups under a suitable class of boundary conditions. Particular attention is paid to the role of modal approximations in energy analysis. The theoretical results are complemented by numerical experiments that illustrate how the regularity of initial data, smooth versus nonsmooth, affects the observed decay rates, providing deeper insight into the interplay between spectral structure and energy dissipation.

28. 一类Brezis–Peletier型结果：关于分数阶Robin函数

A Brezis and Peletier type result for the fractional Robin function

作者: Sidy M. Djitte, Franck Sueur

类别: math.AP

发布日期: 2026-02-06

arXiv ID: 2602.07221v1

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

摘要: This paper is devoted to the Laplacian operator of fractional order $s\in (0,1)$ in several dimensions. We consider the equation $$(-Δ)^su=f(x,u)$$ in $$Ω$$ , $$u=0$$ in $$Ω^c$$ and establish a representation formula for partial derivatives of solutions in terms of the normal derivative $$u/δ^s$$ . As a consequence, we prove that solutions to the overdetermined problem $$(-Δ)^su=f(x,u)$$ in $$Ω$$ , $$u=0$$ in $$Ω^c$$ , and $$u/δ^s=0$$ on $\partialΩ$ are globally Lipschitz continuous provided that $$2s>1$$ . We also prove a Pohozaev-type identity for the Green function and, in particular, obtain a formula for the gradient of the Robin function, which extends to the fractional setting some results obtained by Brézis and Peletier in \cite{Bresiz} in the classical case of the Laplacian. Finally, an application to the nondegeneracy of critical points of the fractional Robin function in symmetric domains is discussed.

29. 分数阶拉普拉斯算子的逐点阿达马变分公式

Pointwise Hadamard variational formula for the fractional Laplacian

作者: Sidy M. Djitte, Franck Sueur

类别: math.AP

发布日期: 2026-02-06

arXiv ID: 2602.07214v1

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

摘要: We establish pointwise formulas for the shape derivative of solutions to the Dirichlet problem associated with the fractional Laplacian. Specifically, we consider the equation $$(-Δ)^s u = h$$ in $$Ω$$ and $$u=0$$ in $$Ω^c$$ , where the right-hand side $$h$$ is either a Dirac delta distribution or a Lipschitz function. In both cases, we prove that the corresponding solution is shape differentiable in every direction and we derive a formula for the pointwise value of its shape derivative. These formulas involve integral on the domain’s boundary and fractional Neumann’s traces. This extends to the case of the fractional Laplacian the well-known Hadamard variational formula for the standard Laplacian. Our argument is in the spirit of \cite{Ushikoshi, Kozono-Ushikoshi} and is based on PDEs techniques.

30. 薛定谔本征函数的半经典局域化

Semiclassical localization of Schrödinger’s eigenfunctions

作者: Sébastien Campagne

类别: math.AP

发布日期: 2026-02-06

arXiv ID: 2602.07128v1

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

摘要: This article addresses the microlocalization of eigenfunctions for the semiclassical Schrödinger operator $$-h^2Δ+V$$ on closed Riemann surfaces with real bounded potentials. Our primary aim is to establish quantitative bounds on the spatial concentration of these eigenfunctions, extending classical results, typically restricted to smooth potentials, to the more general case where the potential is merely bounded. Our main result provides an explicit exponential bound for the $$L^2$$ -norm of eigenfunctions on the entire surface in terms of their $$L^2$$ -norm on an arbitrary open subset with an exponential weight of $Ch^{-1}\log(h)^2$ . This bound improves upon previous estimates for non-smooth potentials that was an exponential weight of $Ch^{-4/3}$ . Our proof is based on a recent approach of the Landis conjecture develop by Logunov, Malinnikova, Nadirashvili and Nazarov (2025).

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

1. 小 $BV \cap L^1$ 初值下动态燃烧中一维可压 Navier–Stokes 方程组的整体适定性#

详细分析#

1. 研究对象和背景#

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

定理 0.2（局部存在性与正则性）#

定理 0.3（稳定性与唯一性）#

定理 0.4（整体存在性与大时间行为）#

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

a) 核心分析工具#

b) 技术性难题的克服#

4. 与之前工作的比较#

a) 突破与改进#

2. 抽象积分微分方程及其应用#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

3. 由分数布朗运动驱动的快慢Navier-Stokes系统的平均化动力学与Wong-Zakai逼近#

详细分析#

1. 研究对象和背景#

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

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

a) 核心分析学工具#

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

4. 与之前工作的比较#

a) 突破与改进#

4. 几何Navier-Stokes方程的反问题#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

b) 技术性难题的克服#

与之前工作的比较#

a) 突破与改进#

5. 雷诺数对薄多孔介质中非牛顿流体流动的影响#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

6. 受点源驱动的黏性Burgers方程：弱极限的显式表达式#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

7. 关于网络上压缩因子非常数的混合气体稳态解存在性的研究#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

8. 相对熵与玻尔兹曼方程的弱可压缩纳维–斯托克斯动力学#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

9. 具有中性嵌入特征值的剪切流的不稳定性#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

1. 小 $BV \cap L^1$ 初值下动态燃烧中一维可压 Navier–Stokes 方程组的整体适定性

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

定理 0.2（局部存在性与正则性）

定理 0.3（稳定性与唯一性）

定理 0.4（整体存在性与大时间行为）

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

a) 核心分析工具

b) 技术性难题的克服

4. 与之前工作的比较

a) 突破与改进

2. 抽象积分微分方程及其应用

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

3. 由分数布朗运动驱动的快慢Navier-Stokes系统的平均化动力学与Wong-Zakai逼近

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

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

a) 核心分析学工具

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

4. 与之前工作的比较

a) 突破与改进

4. 几何Navier-Stokes方程的反问题

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

b) 技术性难题的克服

与之前工作的比较

a) 突破与改进

5. 雷诺数对薄多孔介质中非牛顿流体流动的影响

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

6. 受点源驱动的黏性Burgers方程：弱极限的显式表达式

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

7. 关于网络上压缩因子非常数的混合气体稳态解存在性的研究

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

8. 相对熵与玻尔兹曼方程的弱可压缩纳维–斯托克斯动力学

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

9. 具有中性嵌入特征值的剪切流的不稳定性

详细分析

研究对象和背景

主要定理或主要结果