Home » 论文总结

2026年01月28日论文分析

共有 9 篇重点关注论文, 5 篇论文可以了解, 38 篇不相关论文

January 28, 2026 · 61 min · 12893 words · P0lyno3ial

Table of Contents

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

生成时间: 2026年01月28日 13:43:15

重点关注论文数量: 9

了解领域论文数量: 5

不相关论文数量: 38

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

1. 一类双参数双曲松弛系统向不可压Navier-Stokes方程的收敛性

Convergence of a two-parameter hyperbolic relaxation system toward the incompressible Navier-Stokes equations

作者: Qian Huang, Christian Rohde, Ruixi Zhang

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19846v1

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

详细分析

中文标题: 一类双参数双曲松弛系统向不可压Navier-Stokes方程的收敛性

研究对象和背景

本文研究的核心对象是一类双参数双曲松弛系统，其目标是作为不可压Navier-Stokes方程（incompressible Navier-Stokes equations, INS）的渐近逼近模型。该系统融合了两种经典数值/建模思想：人工可压缩性方法（artificial compressibility, AC）与一阶应力松弛机制（first-order relaxation of the stress-like tensor），从而在保持有限传播速度与因果性的同时，规避传统INS中因强约束 $\nabla \cdot u = 0$ 所引发的全局压力Poisson求解困难。

设空间域为 $$d$$ -维环面 $\mathbb{T}^d = \mathbb{R}^d / \mathbb{Z}^d$ ，其中 $$d = 2$$ 或 $$3$$ ；时间区间为 $$[0,T)$$ 。记未知函数 $(p^{\epsilon,\delta}, u^{\epsilon,\delta}, U^{\epsilon,\delta})$ 分别为标量压力、向量值速度场及矩阵值类应力张量，满足如下初值问题：

\begin{cases} \epsilon \partial_t p^{\epsilon,\delta} + \nabla \cdot u^{\epsilon,\delta} = 0, & \text{(AC型连续性修正)} \\ \partial_t u^{\epsilon,\delta} + \nabla \cdot U^{\epsilon,\delta} + \nabla p^{\epsilon,\delta} = 0, & \text{(动量守恒)} \\ \delta \partial_t U^{\epsilon,\delta} + \nabla u^{\epsilon,\delta} = u^{\epsilon,\delta} \otimes u^{\epsilon,\delta} - U^{\epsilon,\delta}, & \text{(应力松弛动力学)} \\ u^{\epsilon,\delta}(0,\cdot) = u_0^{\epsilon,\delta},\quad p^{\epsilon,\delta}(0,\cdot) = p_0^{\epsilon,\delta},\quad U^{\epsilon,\delta}(0,\cdot) = U_0^{\epsilon,\delta}. \end{cases} \tag{1.2}

其中 $\epsilon > 0$ 和 $\delta > 0$ 是两个独立的小参数，分别控制人工可压缩性强度与应力松弛速率。当 $\epsilon,\delta \to 0$ 时，形式上系统 (1.2) 趋于经典的不可压Navier-Stokes方程：

\begin{cases} \nabla \cdot u = 0, \\ \partial_t u + \nabla \cdot (u \otimes u) = -\nabla p + \Delta u, \\ u(0,\cdot) = u_0. \end{cases} \tag{1.1}

注意：论文中粘性系数被归一化为 $$1$$ ，且对压力施加零均值条件 $\int_{\mathbb{T}^d} p(t,x)\,dx = 0$ 以保证唯一性。该设定下，(1.1) 在二维情形 $$d=2$$ 具有全局光滑解；而在三维情形 $$d=3$$ ，仅对小初值（如 $u_0 \in H^{1/2}(\mathbb{T}^3)$ 或 $L^3(\mathbb{T}^3)$ ）能保证全局存在性，一般初值仅知局部存在 Leray–Hopf 弱解或短时光滑解。

本文所考虑的函数空间为周期 Sobolev 空间 $H^s(\mathbb{T}^d)$ ，并采用标准记号：

$\|u\|_{C^0 H^s} := \sup_{t \in [0,T]} \|u(t,\cdot)\|_{H^s}$ ，
$\|u\|_{C^0 H^s \cap C^1 H^{s-1}} := \sup_{t \in [0,T]} \left( \|u(t,\cdot)\|_{H^s}^2 + \|\partial_t u(t,\cdot)\|_{H^{s-1}}^2 \right)^{1/2}$ 。

特别地，在三维分析中频繁使用 Hodge 分解恒等式：

\|\nabla v\|_{L^2(\mathbb{T}^3)}^2 = \|\operatorname{div} v\|_{L^2(\mathbb{T}^3)}^2 + \|\operatorname{curl} v\|_{L^2(\mathbb{T}^3)}^2, \tag{2.26}

以及 Sobolev 嵌入 $H^1(\mathbb{T}^3) \hookrightarrow L^6(\mathbb{T}^3)$ ，即

\|\phi\|_{L^6(\mathbb{T}^3)} \le C \|\nabla \phi\|_{L^2(\mathbb{T}^3)}, \quad \text{for } \phi = \omega^{\epsilon,\delta} - \omega \text{ or } \omega. \tag{2.27}

此外，为统一处理误差项，引入尺度化误差变量：

v = u^{\epsilon,\delta} - u,\quad q = \sqrt{\epsilon}(p^{\epsilon,\delta} - p),\quad V = \sqrt{\delta}(U^{\epsilon,\delta} - U), \tag{2.24}

它们满足一个对称双曲系统：

\begin{cases} \partial_t q + \frac{1}{\sqrt{\epsilon}} \nabla \cdot v = -\frac{1}{\sqrt{\epsilon}} \partial_t p, \\ \partial_t v + \frac{1}{\sqrt{\delta}} \nabla \cdot V + \frac{1}{\sqrt{\epsilon}} \nabla q = 0, \\ \partial_t V + \frac{1}{\sqrt{\delta}} \nabla v = 2 u \odot v + \sqrt{\delta}\, v \otimes v - \frac{1}{\delta} V - \sqrt{\delta}\, \partial_t U, \end{cases} \tag{2.25}

其中 $u \odot v = \frac{1}{2}(u \otimes v + v \otimes u)$ 。该系统的对称双曲结构是进行能量估计的基础。

主要定理或主要结果

本文共建立三类核心收敛结果，按初始扰动幅度分为两个主范式：

【小扰动情形】 $\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = o(1)$

定理 2.1（ $$H^1$$ -速度收敛， $$d=3$$ ）
假设存在常数 $a \in [0,2)$ 和 $$C_0 > 0$$ ，使得初值满足：

\|u_0^{\epsilon,\delta} - u_0\|_{H^1}^2 + (\epsilon+\delta)^a \|\nabla^2(u_0^{\epsilon,\delta} - u_0)\|_{L^2}^2 \le C_0 (\epsilon+\delta),

\|(p_0^{\epsilon,\delta}, U_0^{\epsilon,\delta})\|_{H^1}^2 + (\epsilon+\delta)^a \|\nabla^2(p_0^{\epsilon,\delta}, U_0^{\epsilon,\delta})\|_{L^2}^2 \le C_0.

则对任意 $$T > 0$$ ，只要 $u^{\epsilon,\delta}, u$ 在 $$[0,T]$$ 上存在，就有：

\sup_{t \in [0,T]} \|u^{\epsilon,\delta}(t,\cdot) - u(t,\cdot)\|_{H^1(\mathbb{T}^3)}^2 \le C_T (\epsilon+\delta).

推论 2.3（ $$L^2$$ -速度收敛， $$d=3$$ ）
若将上述 $$H^1$$ 条件替换为 $$L^2$$ 初值控制（ $a \in [0,4)$ ），则仍有：

\sup_{t \in [0,T]} \|u^{\epsilon,\delta}(t,\cdot) - u(t,\cdot)\|_{L^2(\mathbb{T}^3)}^2 \le C_T (\epsilon+\delta).

定理 2.5（压力 $$H^1$$ -收敛， $$d=3$$ ）
若进一步假设 $\delta \le C_0 \sqrt{\epsilon}$ 及更强初值正则性：

\|u_0^{\epsilon,\delta} - u_0\|_{H^1} + \sqrt{\epsilon}\,\|p_0^{\epsilon,\delta} - p_0\|_{H^1} + \sqrt{\delta}\,\|U_0^{\epsilon,\delta} - U_0\|_{H^1} \le C_0(\epsilon+\delta),

\|\nabla \operatorname{div} u_0^{\epsilon,\delta}\|_{L^2} + \delta\big( \|\nabla^2 p_0^{\epsilon,\delta}\|_{L^2} + \|\nabla \operatorname{div} U_0^{\epsilon,\delta}\|_{L^2} \big) \le C_0(\epsilon+\delta),

\|\nabla^2(u_0^{\epsilon,\delta}, \sqrt{\epsilon}p_0^{\epsilon,\delta}, \sqrt{\delta}U_0^{\epsilon,\delta})\|_{L^2} \le C_0 (\epsilon+\delta)^{-a}, \quad a \in [0,1),

则：

\sup_{t \in [0,T]} \Big( \|u^{\epsilon,\delta}(t,\cdot) - u(t,\cdot)\|_{H^1} + \sqrt{\epsilon}\,\|p^{\epsilon,\delta}(t,\cdot) - p(t,\cdot)\|_{H^1} \Big) \le C_T (\epsilon+\delta).

特别地，若 $\delta = o(\sqrt{\epsilon})$ ，则 $\|p^{\epsilon,\delta} - p\|_{H^1} \to 0$ 。

【大扰动情形】 $\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = \mathcal{O}(1)$

定理 2.6 与 2.7（ $$d=3$$ ，两类不同尺度结构）
二者均要求 $\epsilon \lesssim \delta$ （Thm 2.6 更强： $\epsilon \lesssim \delta^2$ ），并分别对初值的 $$L^2$$ 、旋度、散度、压力与应力张量施加精细的 $\delta$ -依赖控制。例如 Thm 2.6 给出：

\|u^{\epsilon,\delta}(t,\cdot) - u(t,\cdot)\|_{L^2}^2 \le C \delta^{3/2},\quad \|\operatorname{curl}(u^{\epsilon,\delta} - u)(t,\cdot)\|_{L^2}^2 + \delta \|\operatorname{curl}(u^{\epsilon,\delta} - u)(t,\cdot)\|_{H^1}^2 \le C \delta^{1/2},

\|\operatorname{div} u^{\epsilon,\delta}(t,\cdot)\|_{L^2}^2 + \delta \|\operatorname{div} u^{\epsilon,\delta}(t,\cdot)\|_{H^1}^2 \le C \delta^{-1/2} \epsilon.

定理 2.8（ $$d=2$$ ，唯一结构）
在二维情形下，仅需单一尺度： $\epsilon \lesssim \delta$ ，并得：

\|u^{\epsilon,\delta}(t,\cdot) - u(t,\cdot)\|_{L^2(\mathbb{T}^2)}^2 \le C \delta.

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

a) 核心分析学工具

对称双曲系统的经典能量估计：
对误差变量 $$W = (q,v,V)$$ 构造能量泛函
$E(t) = \|W(t,\cdot)\|_{H^1}^2 + (\epsilon+\delta)^a \|\nabla^2 W(t,\cdot)\|_{L^2}^2,$
利用系统 (2.25) 的对称性，导出微分不等式：
$\frac{d}{dt}E(t) \le C\Big( \sqrt{\epsilon+\delta}\,E(t) + E(t) + (\epsilon+\delta)^{-a/2}E(t)^2 \Big).$
此类 Gronwall 型不等式结合初始能量小性与参数依赖的 blow-up 准则（引理 2.9），可同时控制解的存在时间下界与误差大小。
中间仿射系统（intermediate affine system）：
为突破压力缺乏直接能量控制的瓶颈，作者构造辅助系统 (3.7)，其右端含线性项 $u \otimes u$ （而非非线性 $u^{\epsilon,\delta} \otimes u^{\epsilon,\delta}$ ）。通过分别估计 $\|p^{\epsilon,\delta} - p'\|_{H^1}$ 与 $\|p' - p\|_{H^1}$ ，再三角不等式合成，成功获得压力的定量收敛率。此技术将非线性耦合解耦为“松弛主导”与“极限主导”两部分。
调制能量（modulated energy）结构：
针对 $\mathcal{O}(1)$ 初值情形，常规 $$H^s$$ 能量估计失效（存在时间随初值反比衰减）。作者引入源自 [4] 的调制能量：
$\mathcal{E}_{\mathrm{mod}}(t) = \frac{1}{2}\|v\|_{L^2}^2 + \frac{\delta}{2}\|\operatorname{curl} v\|_{L^2}^2 + \cdots,$
其耗散项含 $\|\nabla \operatorname{curl} u^{\epsilon,\delta}\|_{L^2}^2$ ，不受人工可压缩性破坏，成为高阶正则性控制的关键。
Bootstrap 论证与 Hodge 分解：
在 $\mathcal{O}(1)$ 情形中，作者设计精巧的 bootstrap 假设：假设 $\|\operatorname{div} u^{\epsilon,\delta}\|_{L^2} \lesssim \delta^{-1/2}\epsilon$ ，再利用 Hodge 恒等式 (2.26) 将 $\|\nabla v\|_{L^2}$ 控制为旋度与散度之和，最终闭环验证假设。该论证依赖于 $\epsilon = o(\delta)$ 的尺度分离。

b) 关键技术难点突破

压力估计的“双重小参数障碍”：传统 AC 方法（单参数 $\epsilon$ ）可得 $p^{\epsilon} \rightharpoonup p$ 弱收敛，但无速率；而单参数松弛模型（如 (1.3)）未定义压力。本文首次在双参数框架下，通过中间系统桥接两者，实现 $\sqrt{\epsilon}\|p^{\epsilon,\delta} - p\|_{H^1} = \mathcal{O}(\epsilon+\delta)$ 。
大初值下的全局适定性恢复：利用调制能量中 $\delta \|\operatorname{curl} v\|_{L^2}^2$ 的强耗散，克服 $U^{\epsilon,\delta}$ -阻尼不足以压制 $v \otimes v$ 非线性的缺陷，实现 $t \in [0,T_0]$ 上一致控制，并借助 $\epsilon = o(\delta)$ 抑制可压缩性扰动。
三维旋度向量化处理：区别于二维标量涡量，三维 $\operatorname{curl} u$ 是向量场，作者通过 $\|\nabla \operatorname{curl} u\|_{L^2} \sim \|\Delta u\|_{L^2}$ 及嵌入 $H^1 \hookrightarrow L^6$ ，将高阶非线性项控制在 $$L^2$$ 框架内。

与之前工作的比较

a) 突破与改进

首次建立双参数同时极限的严格理论：
此前工作多限于单参数极限：AC 方法（ $\epsilon \to 0$ 固定 $\delta$ ）或松弛方法（ $\delta \to 0$ 固定 $\epsilon$ ）。本文首次处理 $(\epsilon,\delta) \to (0,0)$ 的联合极限，揭示参数间必须满足 $\delta = o(\sqrt{\epsilon})$ （压力收敛）或 $\epsilon = o(\delta)$ （大初值全局恢复）等隐含相容性条件，填补了双参数奇摄动理论的空白。
压力收敛的首个定量结果：
文献 [10,19,31] 对 AC 方法仅给出弱*收敛 $p^\epsilon \stackrel{*}{\rightharpoonup} p$ ；文献 [4,14,27] 对松弛模型甚至未定义压力。本文定理 2.5 是首个在三维中给出 $\|p^{\epsilon,\delta} - p\|_{H^1} = \mathcal{O}(\sqrt{\epsilon+\delta})$ 的定量估计，且明确刻画了 $\epsilon$ - $\delta$ 耦合机制。
大初值鲁棒性突破：
以往能量法要求初值 $\mathcal{O}((\epsilon+\delta)^\alpha)$ （ $\alpha>0$ ），本文定理 2.6–2.8 将初值扰动放宽至 $\mathcal{O}(1)$ ，并仍保证全局时间存在性（ $T^{\epsilon,\delta} \ge T_0$ 与 $\epsilon,\delta$ 无关），这在数值启动（initial inconsistency）场景中具有直接应用价值。
维度普适性与结构识别：
发现三维存在两类收敛速率结构（Thm 2.6 vs. 2.7），源于旋度 $$H^1$$ 范数与 $$L^2$$ 范数的不同标度行为；而二维仅一种结构（Thm 2.8）。这揭示了高维流体力学中几何自由度增加导致的渐近复杂性跃升，为后续多尺度建模提供理论判据。

综上，本文不仅推进了松弛逼近理论的深度，更构建了一套可迁移的双参数奇摄动分析范式，涵盖能量估计、中间系统构造、调制能量与 bootstrap 闭环等完整工具链，为可压缩/不可压过渡、湍流统计模拟等前沿方向奠定基础。

2. 基于凸积分框架的低马赫数极限普适性研究

Universality in the Low Mach number limit via a convex integration framework

作者: Robin Ming Chen, Alexis Vasseur, Dehua Wang, Cheng Yu

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19744v1

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

详细分析

中文标题: 基于凸积分框架的低马赫数极限普适性研究

研究对象和背景

本文研究的核心对象是可压缩欧拉方程组在低马赫数极限下的奇异极限行为，其数学本质是考察当声速趋于无穷（即 Mach 数 $\delta \to 0$ ）时，可压缩流体动力学系统如何趋近于不可压缩欧拉系统。该问题属于偏微分方程中的奇异摄动理论与弱解理论交叉领域，具有深刻的物理动机与数学挑战。

论文从经典流体力学出发，明确给出了两类基本方程组：
首先是 $$n=2$$ 或 $$3$$ 维环面 $\mathbb{T}^n$ 上的不可压缩欧拉方程组（Incompressible Euler equations）：

\begin{cases} \partial_t u + u \cdot \nabla u + \nabla \pi = 0, \\ \operatorname{div} u = 0, \end{cases} \tag{1.2}

其中 $u = u(x,t) \in \mathbb{R}^n$ 是速度场， $\pi = \pi(x,t)$ 是压力项，定义在时空区域 $(x,t) \in \mathbb{T}^n \times [0,T]$ 上。该系统描述理想、无粘、密度恒定的流体运动；其解空间设定为 $u \in L^\infty(0,T;L^2(\mathbb{T}^n))$ ，且满足分布意义下的散度自由条件及动量守恒（见 Definition 1.1）。特别地，文中未施加任何光滑性假设——允许 $$u$$ 是任意 $$L^2$$ 意义下的弱解，甚至可能包含激波、涡旋等非正则结构。

其次，对应地引入经尺度变换后的可压缩欧拉方程组（Rescaled compressible Euler equations），其形式为：

\begin{cases} \partial_t \rho_\delta + \operatorname{div}(\rho_\delta v_\delta) = 0, \\ \partial_t (\rho_\delta v_\delta) + \operatorname{div}(\rho_\delta v_\delta \otimes v_\delta) + \frac{1}{\delta^2} \nabla P(\rho_\delta) = 0, \end{cases} \tag{1.3}

其中状态方程取为 等熵压强律 $P(\rho_\delta) = a \rho_\delta^\gamma$ ，参数 $\gamma > 1$ ， $$a > 0$$ ；而尺度参数 $\delta > 0$ 表征 Mach 数大小。通过变量替换 $\rho_\delta(x,t) = \rho(x,\delta t),\; v_\delta(x,t) = \delta v(x,\delta t)$ ，原物理可压缩系统被重标度为上述形式，使得 $\delta \to 0$ 对应于“慢时间+小速度”极限，从而凸显不可压缩性主导效应。

值得注意的是，本文所考虑的区域为周期性环面 $\mathbb{T}^n$ ，即无边界情形（free boundary or periodic boundary condition），这避免了边界层效应带来的技术复杂性，同时便于应用傅里叶分析与调和分析工具。函数空间方面，对可压缩系统要求初始数据 $(\rho_0,m_0) \in C_0(\mathbb{T}^n) \times L^2(\mathbb{T}^n)$ ，并构造出满足 Definition 1.2 的弱解 $(\rho_\delta, \rho_\delta v_\delta)$ ，其中关键能量型范数包括：

$\rho_\delta \in L^\infty(0,T;L^\gamma(\mathbb{T}^n))$ ，
$\sqrt{\rho_\delta} v_\delta \in L^\infty(0,T;L^2(\mathbb{T}^n))$ ，
动量 $\rho_\delta v_\delta$ 在测试函数空间中满足分布意义的守恒律。

这种设定完全脱离传统强解框架（如 Klainerman–Majda 的光滑初值理论）或相对熵方法所需的高正则性与适定性前提，转而立足于现代弱解存在性理论，尤其是由 De Lellis–Székelyhidi 开创、并在 [5] 中推广至一般 Reynolds 应力张量的凸积分（convex integration）方法。

凸积分本质上是一种构造性反例技术，它揭示了在 $L^\infty$ 或 $$L^2$$ 类弱解中，Euler 方程存在病态非唯一性：对稠密初值集，存在无限多个满足能量不等式的弱解。这一现象颠覆了经典适定性观念，也促使作者提出一个全新视角的问题：若不可压缩极限常被视为“平均行为”或“宏观描述”，那么是否每个不可压缩弱解 $$u$$ 都能作为某个可压缩弱解族 $\{(\rho_\delta, v_\delta)\}_\delta$ 的极限？换言之，不可压缩系统是否是一个普适吸引子（universal attractor）？

该问题直指当前理论的结构性缺失：以往工作（如 Feireisl–Klingenberg–Markfelder 的 DMV 解理论）仅证明“某类可压缩解收敛到某个不可压缩解”，但无法回答“是否所有不可压缩解都可被逼近”。本文首次将凸积分从“制造病态解”的工具，升华为“实现任意目标极限”的通用逼近引擎，从而将低 Mach 数极限的研究提升至密度型泛函分析高度。

主要定理或主要结果

本文的核心成果由以下定理完整刻画：

Theorem 1.1. 设 $$n = 2$$ 或 $$3$$ ， $\gamma > 1$ ，给定任意 $$T > 0$$ 及任意满足 Definition 1.1 的不可压缩 Euler 方程弱解 $u \in L^\infty(0,T;L^2(\mathbb{T}^n))$ ，其初值为 $u_0 \in L^2(\mathbb{T}^n)$ 。则对任意 $$T_1 < T$$ ，存在 $\delta_0 > 0$ ，使得对任意 $0 < \delta < \delta_0$ ，存在一族初值 $(\rho_0^\delta,m_0^\delta) \in C_0(\mathbb{T}^n) \times L^2(\mathbb{T}^n)$ ，满足如下性质：

对每个固定 $\delta > 0$ ，存在无穷多个弱解 $(\rho_\delta, \rho_\delta v_\delta)$ 满足 Definition 1.2，即求解 (1.3)，且以 $(\rho_0^\delta,m_0^\delta)$ 为初值；
下列强收敛性成立： $\rho_\delta \to 1 \quad \text{in } C^0([0,T_1] \times \mathbb{T}^n),$ $\rho_\delta v_\delta \to u \quad \text{in } L^2(0,T_1;L^2(\mathbb{T}^n)),$ $m_0^\delta \rightharpoonup u_0 \quad \text{weakly in } L^2(\mathbb{T}^n).$

该定理断言：任意 $$L^2$$ 意义下的不可压缩弱解 $$u$$ ，均可被一族由凸积分构造的可压缩弱解强逼近。这是首个在无能量不等式约束、无初值正则性假设、任意 $\gamma > 1$ 条件下建立的普适性密度结果，彻底解决了“不可压缩极限是否普适”的根本问题。

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

a) 核心分析学工具

本文的证明骨架建立在三大支柱之上：

凸积分框架（Convex Integration Framework）：继承并深化了 De Lellis–Székelyhidi 及作者前期工作 [5] 的思想，将 PDE 求解转化为在某个凸集上寻找扰动向量场的问题。具体而言，先构造一个“松弛解”（subsolution） $(\rho_0,V_0,R_0)$ 满足线性化系统（3.1）与 Reynolds 应力误差项 $$R_0$$ ，再通过迭代添加高频振荡扰动 $$V_k$$ ，逐步减小 $$R_k$$ 直至消失。此过程依赖于Baire 范畴定理：在适当拓扑（弱- $\ast$ 拓扑）下，满足饱和条件（3.3）的解构成剩余集（residual set），从而保证无穷多解的存在性。
Reynolds 应力控制技术：区别于传统凸积分仅关注 $$L^1$$ 或 $$C^0$$ 范数衰减，本文引入关键创新—— $$L^2$$ 型约束条件（1.4）：
$\left| \int_{[0,T]\times \mathbb{T}^n} \frac{V \cdot \delta V}{\rho_\delta} \,dxdt \right| \leq \frac{1}{8} \int_{[0,T]\times \mathbb{T}^n} \frac{|V|^2}{\rho_\delta} \,dxdt.$
此不等式强制扰动动能受控于 Reynolds 张量的迹，导出统一估计：
$\int_{[0,T]\times \mathbb{T}^n} \frac{|V|^2}{\rho_\delta} \,dxdt \leq C \int_{[0,T]\times \mathbb{T}^n} \operatorname{tr} R_\delta \,dxdt.$
结合 subsolution 构造中 $\operatorname{tr} R_\delta = O(\delta)$ ，立即得到扰动能量一致趋于零，保障强收敛。
Subsolution 提升技术（Lifting Scheme）：为连接不可压缩与可压缩系统，需将 $(u_\varepsilon,\pi_\varepsilon,R_\varepsilon)$ 映射为 $(\rho_\delta,v_\delta,\tilde{R}_\delta)$ 。关键技术在于设计密度修正项 $\varrho_{\delta,K^*}$ （见 (2.3a)），使其满足质量守恒与平均密度为 1 的约束（2.3b），并利用 Bourgain–Brezis 方程求解器处理 $\partial_t \rho_\delta + \operatorname{div}(\rho_\delta v_\delta)=0$ ，确保 $\rho_\delta > 0$ 且保持正定性。

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

Mach 数耦合难题：传统凸积分中 $\delta$ 视为固定参数，而此处需让整个构造一致依赖于 $\delta$ 并在 $\delta \to 0$ 时保持稳定性。作者通过精细选取 mollification 参数 $\alpha_\varepsilon$ 与修正幅值 $$K^*$$ ，使所有误差项（如 $\|\rho_\delta - 1\|_{C^0},\; \|v_\delta - u_\varepsilon\|_{L^2}$ ）均以 $O(\varepsilon + \delta)$ 控制，最终令 $\varepsilon = \delta$ 实现统一速率。
压力项奇异性处理：方程 (1.3) 含 $\delta^{-2}\nabla \rho_\delta^\gamma$ ，当 $\rho_\delta \to 1$ 时产生大梯度。作者巧妙将其拆分为 $\delta^{-2}\nabla(\rho_\delta^\gamma - 1) = \nabla \pi_\varepsilon + \text{error}$ ，并将主部吸收进不可压缩压力项，余项由 $\tilde{R}_\delta$ 承载，从而规避直接估计高阶导数。
弱- $\ast$ 收敛与非线性项匹配：在 Baire 类论证中，需验证扰动序列在弱- $\ast$ 下收敛时，非线性项 $V \otimes V$ 仍能逼近目标 $$R_0$$ 。为此，作者发展了改进的 $$L^1$$ -coercivity 引理（Lemma 3.1），借助局部立方体离散化与最小特征值控制，证明扰动 $$V_i$$ 的 $$L^1$$ 范数受 $\operatorname{tr} M_i$ 控制，从而闭环迭代。

与之前工作的比较

a) 突破与改进

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

从“选择性收敛”到“普适逼近”：此前所有结果（如 Ebin [10], Klainerman–Majda [15], Feireisl et al. [12]）均属“正向极限”——给定可压缩解族，证其收敛到某个不可压缩解。本文首次实现“逆向构造”：对任一不可压缩弱解 $$u$$ ，显式构造可压缩解族逼近之。这将低 Mach 数极限由“筛选机制”升格为“全息映射”。
彻底摆脱能量不等式束缚：[5] 中要求 $\gamma \leq 1 + \frac{2}{n}$ 以保证能量不等式成立；而本文因不施加该约束，将适用范围拓展至全部 $\gamma > 1$ （Remark 1.1），覆盖从单原子气体 ( $\gamma = 5/3$ ) 到软物质 ( $\gamma \to 1^+$ ) 的广泛物理模型。
初值正则性归零：以往工作（包括 DMV 解）至少要求有限能量初值，而本文初值 $(\rho_0^\delta,m_0^\delta)$ 可取 $C_0 \times L^2$ ，甚至可逼近 Dirac 型奇异性，真正实现“任意弱解均可被实现”。

综上，本文不仅提供了一个强大新工具，更重塑了我们对奇异极限本质的理解：它不再是某种“理想化约化”，而是可压缩系统在弱拓扑下对不可压缩系统的满射式稠密嵌入。

3. 关于排斥性无压欧拉-泊松系统的渐近行为

On the asymptotic behavior of the Repulsive Pressureless Euler-Poisson System

作者: Nicholas Biglin, Joseph Crachiola, Jack Curtis, Thomas Kunz, Omkar Maralappanavar, Adrian Tudorascu

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19733v1

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

详细分析

中文标题: 关于排斥性无压欧拉-泊松系统的渐近行为

1. 研究对象和背景

本文研究的是一维排斥性无压欧拉-泊松系统（Repulsive Pressureless Euler-Poisson System），记为 repulsive PEP。该系统建模了定义在实轴 $\mathbb{R}$ 上的质量分布的演化动力学，其中每个质量元通过一个长程排斥相互作用势施加向外的力。与经典欧拉方程不同，该系统无压力项（即 $p \equiv 0$ ），因此动量输运完全由粒子间非局部排斥力驱动；同时，它又区别于无相互作用的压力less Euler 系统（即自由粒子流），因为其动量方程右侧包含一个非零的 Poisson 型源项。

该系统的弱形式（distributional formulation） 在 Introduction 第 1.1 节及第 2.2 节中明确给出，对应如下偏微分方程组：

\begin{cases} \partial_t \rho + \partial_y (\rho v) = 0, \\ \partial_t (\rho v) + \partial_y (\rho v^2) = \dfrac{1}{2} \big( \mathrm{sgn} * \rho \big) \rho, \end{cases} \tag{1.1}

其中：

$\rho = \rho(t,y) \ge 0$ 是质量密度，取值为 $\mathbb{R}_+$ 上的 Borel 概率测度（即 $\rho_t \in \mathcal{P}_2(\mathbb{R})$ ，具有有限二阶矩）；
$$v = v(t,y)$$ 是速度场，定义在 $\rho_t$ -a.e. 处，且满足 $v_t \in L^2(\rho_t)$ ；
$\mathrm{sgn}(y) = \begin{cases} -1, & y < 0 \\ 0, & y = 0 \\ 1, & y > 0 \end{cases}$ ，卷积 $(\mathrm{sgn} * \rho)(y) := \int_{\mathbb{R}} \mathrm{sgn}(y - z)\, \rho(dz)$ 表示净排斥加速度：对位于 $$y$$ 的质量元，其加速度等于右侧总质量减去左侧总质量（因 $\int \mathrm{sgn}(y-z)\,\rho(dz) = \rho((y,\infty)) - \rho((-\infty,y))$ ）。

该模型源自宇宙学中大尺度结构形成的简化描述，但与 Zeldovich 提出的吸引性版本（右侧为 $-\frac{1}{2}(\mathrm{sgn}*\rho)\rho$ ）本质不同：此处力是排斥的，导致粒子倾向于彼此远离，而非聚团。

论文特别聚焦于粘性粒子解（sticky particle solutions, SPS）——即离散情形下，当有限个粒子发生碰撞时，它们立即粘合为一个新粒子，新质量为原质量之和，新速度由动量守恒确定：

m_{\text{new}} = \sum_{i\in I} m_i,\quad v_{\text{new}} = \frac{1}{m_{\text{new}}} \sum_{i\in I} m_i v_i.

该规则自然诱导出经验测度解 $(\rho_t,v_t)$ ，满足 (1.1) 在分布意义下成立。值得注意的是，SPS 并非对所有初始数据都存在或唯一；事实上，第 2.3 节通过 Example 2.3 和 2.4 明确指出：排斥性 PEP 系统在粘性框架下是不适定的（ill-posed）：存在初始数据序列 $\{(\rho_0^n,v_0^n)\}$ 收敛于 $(\rho_0,v_0)$ ，但其对应的 SPS 解序列 $\{(\rho_t^n,v_t^n)\}$ 在 Wasserstein 距离 $$W_p$$ 下不收敛，甚至发散（如 $W_2(\rho_t^n,\rho_t) \sim t^2$ ）。这迫使作者引入更广义的解类：广义粘性粒子解（generalized sticky particle solutions），定义为所有离散 SPS 解在 $$W_2$$ 拓扑下的闭包 $\overline{\mathcal{S}}$ 。

函数空间设定极为关键：解始终取值于概率测度空间 $\mathcal{P}_2(\mathbb{R})$ ，而非传统 Sobolev 或 Hölder 空间。这是由于无压极限下，解天然呈现集中奇异性（Dirac delta 簇），无法用常规函数描述。文中反复使用 Wasserstein-2 距离：

W_2(\mu,\nu) = \left( \inf_{\gamma \in \Pi(\mu,\nu)} \int_{\mathbb{R}^2} |x-y|^2 \, d\gamma(x,y) \right)^{1/2},

其中 $\Pi(\mu,\nu)$ 是以 $\mu,\nu$ 为边缘的耦合测度集。等价地，利用累积分布函数 $M_\mu(y)=\mu((-\infty,y])$ 及其广义逆 $N_\mu(m)=\inf\{y:M_\mu(y)>m\}$ ，有：

W_2^2(\mu,\nu) = \int_0^1 |N_\mu(m) - N_\nu(m)|^2 \, dm.

这一表示将测度演化转化为单位区间 $$(0,1)$$ 上的函数演化，是分析的核心几何工具。

此外，系统具有Galilean 不变性（Proposition 2.2, Item 2）：若定义初始质心 $Y_0 = \int y\,\rho_0(dy)$ 与总动量 $m_0 = \int v_0(y)\,\rho_0(dy)$ ，则平移后的 $(\bar{\rho}_t,\bar{v}_t)$ 仍是解，且满足 $\int y\,\bar{\rho}_t(dy) = \int v_t(y)\,\bar{\rho}_t(dy) = 0$ 。因此，不失一般性，全文常设 $$Y_0 = m_0 = 0$$ ，将分析归约至零质心、零总动量的中心化情形。

边界条件隐含于测度的紧支性或衰减性假设中（如 Example 2.1 中初始支集为 $\{-1,1\}$ ），而全局存在性无需显式边界条件，因 $\mathbb{R}$ 无边界。

2. 主要定理或主要结果

本文的核心成果围绕能量演化、有限时间坍缩（finite-time collapse）与完美平衡态（perfect states）展开，具体如下：

Theorem 4.5（必要坍缩条件）：设初始为 $$n$$ 粒子离散分布 $\rho_0 = \sum_{i=1}^n m_i \delta_{y_i(0)}$ ， $$v_0(y_i(0)) = v_i(0)$$ ，且满足 $$Y_0 = m_0 = 0$$ 。若存在有限时间 $T < \infty$ 使得 $\rho_T = \delta_{y_*}$ （单点坍缩），则必有：
$\sum_{i=1}^n m_i y_i^2(0) < \frac{1}{2} \left( \sum_{i=1}^n m_i |v_i(0)| \right)^2.$
Theorem 4.6（充分坍缩条件，离散情形）：同上设定，若存在常数 $$C > 0$$ 使得对所有 $$i$$ ，
$|v_i(0)| \le C - |y_i(0)|,$
且 $C^2 > 2 \sum_i m_i y_i^2(0)$ ，则坍缩在有限时间内发生。
Theorem 4.7（充分坍缩条件，一般情形）：对一般 $\rho_0 \in \mathcal{P}_2(\mathbb{R})$ ， $v_0 \in L^2(\rho_0)$ ，若
$\int y^2 \, \rho_0(dy) < \frac{1}{2} \left( \int |v_0(y)| \, \rho_0(dy) \right)^2,$
且 $$v_0$$ 满足某种 Lipschitz 型控制，则坍缩发生。
Theorem 4.1（二次包络）：若解 $\rho_t$ 在 $$[0,T]$$ 内坍缩，则对所有 $t \in [0,T]$ ，其支撑必须包含于区间 $$[-a(t), a(t)]$$ ，其中 $$a(t) = A + Bt + Ct^2$$ ，且系数 $$A,B,C$$ 由初始二阶矩与动能决定。此为坍缩的必要空间约束。
Existence and Uniqueness of Perfect States（“完美态”）：存在唯一的解 $(\rho_t,v_t)$ ，其Hamiltonian（总能量）
$\mathcal{H}(t) = \frac{1}{2} \int v^2(t,y) \, \rho_t(dy) + \frac{1}{4} \iint |y-z| \, \rho_t(dy)\rho_t(dz)$
恒为常数（即 $\mathcal{H}'(t) \equiv 0$ ），且 $\rho_t \to \delta_0$ 弱收敛。此类解被称为 perfect，代表能量最保守的坍缩路径（最慢收敛）。
Proposition 2.2（运动学律）：在无碰撞区间 $$(0,t_0)$$ 内， $$n$$ 粒子解的轨迹满足：
$y_i(t) = y_i(0) + t v_i(0) + \frac{t^2}{4} \left( \sum_{j: y_j(0)<y_i(0)} m_j - \sum_{j: y_j(0)>y_i(0)} m_j \right).$
Ill-posedness Result（隐含于 Example 2.3–2.4）：粘性解类 $\mathcal{S}$ 在 $$W_2$$ 拓扑下不闭，且不稳定；其闭包 $\overline{\mathcal{S}}$ 中的解可能不满足粘性条件，故需广义化。

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

a) 核心分析学工具

论文深度融合了最优传输理论、测度微分方程与动力系统视角。

Wasserstein-2 几何：全篇以 $$W_2$$ 距离为度量基础，利用其 Hilbert 结构（通过 $N_\mu$ 映射到 $$L^2(0,1)$$ ）将测度演化线性化。例如，Proposition 2.2 的证明依赖于 $$W_2$$ 梯度流解释：方程 (1.1) 可视为 $\mathcal{P}_2(\mathbb{R})$ 上某泛函的梯度流，但因排斥势非凸，标准 McCann 理论不适用，故转向构造性方法。
粒子追踪与碰撞图（collision graph）：对离散解，将粒子轨迹 $$y_i(t)$$ 视为抛物线（由 Proposition 2.2 给出），碰撞时刻即为 $$y_i(t) = y_j(t)$$ 的最小正根。通过分析抛物线族的相对位置（如 Example 2.1 的图 2.1），判断“交叉后是否再分离”，从而判定粘性是否被破坏。此属初等动力系统分析，但需精细代数运算。
能量估计与 Hamiltonian 守恒：定义总能量 $\mathcal{H}(t)$ 后，直接计算其导数：
$\frac{d}{dt}\mathcal{H}(t) = \int v \, \partial_t v \, \rho_t(dy) + \frac{1}{2} \iint \mathrm{sgn}(y-z) (v(y)-v(z)) \, \rho_t(dy)\rho_t(dz).$
利用动量方程 (1.1)₂ 代入，可证：仅当 $$v$$ 为 $$y$$ 的线性函数且 $\rho_t$ 对称时， $\mathcal{H}'(t) = 0$ 。此即 perfect state 的特征。
投影公式失效的严格证伪：针对文献中吸引性系统的成功投影法（ $N_t = \mathrm{proj}_{\mathcal{K}}(N_0 + t v_0(N_0) + \frac{t^2}{2}(0.5-\mathrm{Id}))$ ），作者构造反例（Example 2.1）并证明：若强行将同一公式用于排斥情形，所得 $$N_t$$ 在交叉后“退投影”，导致粒子分离，违反粘性。此涉及凸包分析（convex envelope $F^{**}$ ）与单调重排理论。

b) 技术性难题处理技巧

奇异性处理：Dirac 测度的乘积 $\rho \cdot \rho$ 在分布意义下无定义，但方程 (1.1)₂ 中的 $(\mathrm{sgn}*\rho)\rho$ 通过配对测试函数定义为：
$\int \varphi(y) (\mathrm{sgn}*\rho)(y) \, \rho(dy) = \iint \varphi(y) \mathrm{sgn}(y-z) \, \rho(dy)\rho(dz),$
此二重积分对 $\varphi \in C_c^1$ 恒有意义，因 $\mathrm{sgn}$ 是有界函数。
粘性条件的弱化：为处理 ill-posedness，作者不追求单个解的稳定性，而研究解集的闭包性质。利用 $$W_2$$ 的完备性，任何 Cauchy 序列 $\{(\rho_t^n,v_t^n)\} \subset \mathcal{S}$ 收敛于某 $(\rho_t,v_t) \in \overline{\mathcal{S}}$ ，后者虽未必粘性，但保留了能量与动量的基本守恒律。
坍缩判据的优化：Theorem 4.5 的必要条件源于对 $\int y^2 \rho_t(dy)$ 的二阶导数估计：
$\frac{d^2}{dt^2} \int y^2 \rho_t(dy) = 2 \int v^2 \rho_t(dy) - \frac{1}{2} \iint |y-z| \rho_t(dy)\rho_t(dz) \le 2 \mathcal{H}(t),$
若初始二阶矩过小而动能过大，则曲率负，强制坍缩。此为比较原理的测度版本。

4. 与之前工作的比较

a) 突破与改进

本文最关键的突破在于首次系统建立排斥性无压欧拉-泊松系统的渐近理论，填补了吸引性与无相互作用情形之间的巨大空白。

与吸引性 PEP 的对比：在吸引情形（如 [1], [9]），投影公式成立，解由 $N_t = \mathrm{proj}_{\mathcal{K}}(N_0 + t v_0(N_0) - \frac{t^2}{2}(0.5-\mathrm{Id}))$ 唯一确定，且所有解均收敛于单点，速率受初始数据控制。本文证明：排斥情形无全局投影表示，且收敛速率可任意慢（呼应 [14]），但更深刻的是揭示了排斥与粘性的根本冲突——前者驱散，后者要求聚合，导致解空间结构剧变。
与压力less Euler 的对比：后者无长程力，粒子匀速运动直至碰撞（[8], [10]）。本文显示：即使加入排斥力，只要足够弱，仍可坍缩；但机制完全不同——非靠初始速度对齐，而靠排斥力诱导的相对加速度差异（见 Proposition 2.2 中的 $\frac{t^2}{4}(\cdots)$ 项）。
对初始数据正则性的放松：先前工作（如 [4], [2]）多要求 $\rho_0$ 为 Dirac 组合或绝对连续。本文直接在 $\mathcal{P}_2(\mathbb{R})$ 中工作，允许任意概率测度，且通过离散逼近（Example 2.4）证明：连续初值可被离散粘性解逼近，但极限未必粘性，从而厘清了逼近策略的适用边界。
结果强度的提升：此前关于排斥系统的结论多限于局部存在性或数值观察。本文给出首个严格的坍缩充要条件（Theorems 4.5–4.7）、精确的二次包络（Theorem 4.1）及完美态的存在唯一性，构成完整的渐近分类框架。

综上，本文不仅解决了特定方程的数学问题，更深化了人们对非局部相互作用、奇异性传播与解稳定性之间张力的理解，为相关模型（如带排斥核的 Keller-Segel、集体动力学）提供了范式性工具。

4. Vlasov-Fokker-Planck-Navier-Stokes 系统的光子-粒子定量极限

Quantitative light-particle limit for the Vlasov-Fokker-Planck-Navier-Stokes system

作者: Young-Pil Choi, Jinwook Jung

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19110v1

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

详细分析

中文标题: Vlasov-Fokker-Planck-Navier-Stokes 系统的光子-粒子定量极限

研究对象和背景

本文研究的是 Vlasov–Fokker–Planck–Navier–Stokes (VFP-NS) 耦合系统 在 light-particle（轻粒子）渐近极限 下的定量流体力学极限行为。该系统建模了大量微小、弥散的带电或中性粒子（如气溶胶、喷雾液滴、光子类准粒子等）在粘性不可压缩流体中的集体运动，广泛应用于化学工程、环境科学（如废水处理）、生物医学（如药物输送）以及机械系统（如燃烧室喷雾动力学）等领域。

系统的核心在于两类物理对象的强耦合：

粒子相 由其在相空间 $(x,\xi)\in\Omega\times\mathbb{R}^2$ 上的分布函数 $f^\varepsilon=f^\varepsilon(t,x,\xi)\ge0$ 描述，满足 Vlasov–Fokker–Planck 方程；
流体相 由速度场 $v^\varepsilon=v^\varepsilon(t,x)$ 和压力 $P^\varepsilon=P^\varepsilon(t,x)$ 描述，满足不可压缩 Navier–Stokes 方程。

其中参数 $\varepsilon>0$ 表征粒子质量相对于流体特征质量的尺度，即“轻粒子”假设：粒子惯性极小，其动量弛豫时间尺度为 $O(\varepsilon)$ ，远快于流体演化时间尺度。因此，当 $\varepsilon\to0$ 时，粒子分布迅速趋向局部平衡态，而流体则趋近于一个解耦的、无外力驱动的 Navier–Stokes 流，同时粒子密度 $\rho^\varepsilon=\int_{\mathbb{R}^2} f^\varepsilon\,d\xi$ 满足一个对流-扩散方程。

论文中明确给出的原始系统（公式 (1.1)）为：

\begin{cases} \displaystyle \partial_t f^\varepsilon + \frac{1}{\varepsilon}\xi\cdot\nabla_x f^\varepsilon = \frac{1}{\varepsilon^2}\nabla_\xi\cdot\big(\nabla_\xi f^\varepsilon + (\xi - \varepsilon v^\varepsilon)f^\varepsilon\big), \\[1em] \displaystyle \partial_t v^\varepsilon + v^\varepsilon\cdot\nabla_x v^\varepsilon + \nabla_x P^\varepsilon - \Delta_x v^\varepsilon = \frac{1}{\varepsilon}\int_{\mathbb{R}^2} (\xi - \varepsilon v^\varepsilon) f^\varepsilon\,d\xi, \\[1em] \nabla_x\cdot v^\varepsilon = 0, \end{cases}

定义在时空域 $(t,x,\xi)\in [0,T]\times\Omega\times\mathbb{R}^2$ 上，其中 $\Omega = \mathbb{T}^2$ （二维环面）或 $\Omega = \mathbb{R}^2$ （全空间）。该设定具有显著的几何与分析双重普适性：既涵盖周期边界下天然的紧性结构，又处理非紧全空间中必须克服的空间衰减与无界性挑战。

值得注意的是，该系统不施加任何外部势场或 confinement 势，粒子仅受流体拖曳力 $F_d = \xi - \varepsilon v^\varepsilon$ 驱动，属于典型的无外场、弱耦合、慢流体-快粒子机制。这种设定使得传统基于势能控制或谱间隙的方法失效，亟需发展新的稳定性框架。

解的空间正则性要求严格遵循 weak entropy solution（弱熵解） 框架（Definition 1.1），其核心是能量-熵联合控制：

$f^\varepsilon \in L^\infty(0,T; L^1(\Omega\times\mathbb{R}^2,(1+|x|+|\xi|^2)dxd\xi)\cap L\log L(\Omega\times\mathbb{R}^2))$ ，即要求一阶空间矩 $\int\!\!\int |x|f^\varepsilon\,dxd\xi$ 和 LlogL 熵 $\int\!\!\int f^\varepsilon\log f^\varepsilon\,dxd\xi$ 一致有界；
$v^\varepsilon \in L^\infty(0,T; L^2(\Omega)) \cap L^2(0,T; H^1(\Omega))$ ，体现流体动能与耗散的平衡；
满足分布意义下的方程及总自由能不等式（Lemma 2.2）： $\mathcal{F}[f^\varepsilon(t),v^\varepsilon(t)] + \frac{1}{\varepsilon^2}\int_0^t\!\!\int\!\!\int |\nabla_\xi f^\varepsilon + (\xi - \varepsilon v^\varepsilon)f^\varepsilon|^2\,dxd\xi ds + \int_0^t\!\!\int |\nabla_x v^\varepsilon|^2\,dx ds \le \mathcal{F}[f^\varepsilon_0,v^\varepsilon_0],$ 其中自由能泛函定义为： $\mathcal{F}[f,v] := \iint_{\Omega\times\mathbb{R}^2}\left(\frac{|\xi|^2}{2} + \log f\right)f\,dxd\xi + \frac{1}{2}\int_\Omega |v|^2\,dx.$

该能量结构不仅保障解的存在性（见 [9,30]），更构成后续相对熵方法的基石。特别地，由于 $\Omega = \mathbb{R}^2$ 时缺乏紧嵌入 $H^1 \hookrightarrow L^\infty$ ，作者引入了关键的 weighted Trudinger–Moser 不等式（Lemma 2.1），以控制形如 $\int \rho^\varepsilon |v^\varepsilon|^2 dx$ 的混合项——这在二维情形下无法由标准 Sobolev 嵌入获得，是全空间分析得以成立的技术前提。

综上，本研究的对象是一个定义在双空间域（位置 $x\in\Omega$ 、速度 $\xi\in\mathbb{R}^2$ ）上的高度非线性、强耦合、奇摄动偏微分方程组；其背景是真实工业与自然场景中普遍存在的稀疏粒子-连续流体相互作用；其数学挑战集中于：跨尺度耦合的刚性、相空间高维性、全空间非紧性、以及缺乏外势导致的缺乏先验衰减。

主要定理或主要结果

本文的核心成果由 Theorem 1.1 给出，它建立了 light-particle 极限下从 VFP-NS 系统到宏观耦合系统的首个定量收敛理论。该定理包含两种初始数据准备情形：well-prepared（良备） 与 scaled well-prepared（尺度良备），分别对应不同强度的初始误差控制。

设 $$T>0$$ ， $(\rho,v)$ 是如下极限系统（公式 (1.3)）的唯一经典解：

\begin{cases} \partial_t \rho + \nabla_x\cdot(\rho v) = \Delta_x \rho, \\ \partial_t v + v\cdot\nabla_x v + \nabla_x P - \Delta_x v = 0, \\ \nabla_x\cdot v = 0, \end{cases}

并满足正则性条件：

\rho \in L^\infty(0,T; \mathcal{P}(\Omega)\cap L^\infty(\Omega)), \quad v,\nabla_x\log\rho \in L^\infty(0,T; W^{1,\infty}(\Omega)) \cap W^{1,\infty}(0,T; L^\infty(\Omega)).

令 $\{(f^\varepsilon,v^\varepsilon)\}_{\varepsilon>0}$ 为系统 (1.1) 在 $$[0,T]$$ 上的一族弱熵解，初值 $\{(f^\varepsilon_0,v^\varepsilon_0)\}_{\varepsilon>0}$ 满足一致界：

\sup_{\varepsilon>0}\Big( \| |x| f^\varepsilon_0 \|_{L^1} + \|f^\varepsilon_0\|_{L^1\cap L\log L} + \|v^\varepsilon_0\|_{L^2} \Big) < \infty.

（I）良备初值情形（Well-prepared initial data）

若初值满足：

\mathcal{H}[f^\varepsilon_0 \,|\, M_{\rho,0}](0) + \|v^\varepsilon_0 - v_0\|_{L^2}^2 \to 0 \quad \text{as } \varepsilon \to 0,

其中 $M_{\rho,0}(\xi) := \frac{\rho(x)}{2\pi}\exp\big(-\frac{|\xi|^2}{2}\big)$ 是以极限密度 $\rho$ 和零流速为参数的局部麦克斯韦分布，则有如下强收敛与定量估计：

$f^\varepsilon \to M_{\rho,0}$ in $L^\infty(0,T; L^1(\Omega\times\mathbb{R}^2))$ ,
$\int_{\mathbb{R}^2} f^\varepsilon\,d\xi \to \rho$ in $L^\infty(0,T; L^1(\Omega))$ ,
$\frac{1}{\varepsilon}\int_{\mathbb{R}^2} \xi f^\varepsilon\,d\xi \to \rho(v - \nabla_x\log\rho)$ in $L^1((0,T)\times\Omega)$ ,
$v^\varepsilon \to v$ in $L^\infty(0,T; L^2(\Omega)) \cap L^2(0,T; H^1(\Omega))$ .

且存在常数 $$C>0$$ （仅依赖于初值、 $$T$$ 和 $(\rho,v)$ ）使得： $$ |f^\varepsilon - M_{\rho,0}|{L^\infty(0,T;L^1)}^2 + |\rho^\varepsilon - \rho|{L^\infty(0,T;L^1)}^2 + |v^\varepsilon - v|_{L^\infty(0,T;L^2)\cap L^2(0,T;H^1)}^2 \

\left|\frac{1}{\varepsilon}\int_{\mathbb{R}^2}\xi f^\varepsilon,d\xi - \rho(v - \nabla_x\log\rho)\right|{L^1((0,T)\times\Omega)}^2 \ \le C\Big( \mathcal{H}f^\varepsilon_0 ,|, M_{\rho,0} + |v^\varepsilon_0 - v_0|{L^2}^2 \Big) + C\varepsilon^2. $$

（II）尺度良备初值情形（Scaled well-prepared initial data）

若初值满足更弱的条件：

d_{\mathrm{BL}}(\rho^\varepsilon_0,\rho_0) \to 0, \quad \varepsilon^2\Big( \mathcal{H}[f^\varepsilon_0 \,|\, M_{\rho,0}](0) + \|v^\varepsilon_0 - v_0\|_{L^2}^2 \Big) \to 0,

其中 $d_{\mathrm{BL}}$ 是有界 Lipschitz 距离，则得到弱收敛与最优弱拓扑定量率：

$f^\varepsilon \rightharpoonup M_{\rho,0}$ in $L^2(0,T;\mathcal{M}(\Omega\times\mathbb{R}^2)) \cap L^1((0,T)\times\Omega\times\mathbb{R}^2)$ ,
$\int_{\mathbb{R}^2} f^\varepsilon\,d\xi \overset{*}{\rightharpoonup} \rho$ in $L^\infty(0,T;\mathcal{M}(\Omega)) \cap L^1((0,T)\times\Omega)$ ,
$\frac{1}{\varepsilon}\int_{\mathbb{R}^2} \xi f^\varepsilon\,d\xi \rightharpoonup \rho(v - \nabla_x\log\rho)$ in $L^1((0,T)\times\Omega)$ ,
$v^\varepsilon \to v$ in $L^2_{\mathrm{loc}}((0,T)\times\Omega)$ （局部强收敛）。

并有如下最优弱收敛速率：

\int_0^T d_{\mathrm{BL}}^2(f^\varepsilon(t), M_{\rho,0}(t))\,dt \le C d_{\mathrm{BL}}^2(\rho^\varepsilon_0,\rho_0) + C\varepsilon^2\Big( \mathcal{H}[f^\varepsilon_0 \,|\, M_{\rho,0}](0) + \|v^\varepsilon_0 - v_0\|_{L^2}^2 \Big) + C\varepsilon^2,

\sup_{0\le t\le T} d_{\mathrm{BL}}^2(\rho^\varepsilon(t),\rho(t)) \le C d_{\mathrm{BL}}^2(\rho^\varepsilon_0,\rho_0) + C\varepsilon^2\Big( \mathcal{H}[f^\varepsilon_0 \,|\, M_{\rho,0}](0) + \|v^\varepsilon_0 - v_0\|_{L^2}^2 \Big) + C\varepsilon^4.

该定理首次在 light-particle 极限中实现了显式收敛速率、统一适用环面与全空间、以及最优弱拓扑精度（ $\varepsilon^2$ 对密度、 $\varepsilon^2$ 对分布函数），填补了此前仅有定性紧性结果（Goudon–Jabin–Vasseur, 2004）的空白。

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

a) 核心分析学工具

本文的证明骨架建立在 modulated relative entropy method（调制相对熵方法） 之上，这是近年来处理 kinetic-fluid 渐近极限最有力的定量工具之一，但此前仅成功应用于 fine-particle（超声速/双曲）标度，而从未用于 light-particle（抛物/扩散）标度。其根本创新在于将经典相对熵 $\mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}]$ 进行动态调制，即选取时刻依赖的参考麦克斯韦分布 $M_{\rho^\varepsilon,u^\varepsilon}$ ，其中宏观场 $(\rho^\varepsilon,u^\varepsilon)$ 由 $f^\varepsilon$ 自身通过矩运算定义：

\rho^\varepsilon(t,x) := \int_{\mathbb{R}^2} f^\varepsilon(t,x,\xi)\,d\xi, \quad u^\varepsilon(t,x) := \frac{1}{\rho^\varepsilon}\int_{\mathbb{R}^2} \xi f^\varepsilon(t,x,\xi)\,d\xi.

由此构造的 modulated entropy：

\mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}] = \iint f^\varepsilon \log\frac{f^\varepsilon}{M_{\rho,u}}\,dxd\xi - \iint (f^\varepsilon - M_{\rho,u})\,dxd\xi,

其时间导数可被精细分解为三部分：
(i) 熵耗散项（positive definite）： $\frac{1}{\varepsilon^2}\iint |\nabla_\xi f^\varepsilon + (\xi - \varepsilon v^\varepsilon)f^\varepsilon|^2\,dxd\xi$ ，主导 $\varepsilon\to0$ 的快速弛豫；
(ii) 流体-动能耦合误差项：含 $\int \rho^\varepsilon |v^\varepsilon - u^\varepsilon|^2\,dx$ 及 $\int \rho^\varepsilon |u^\varepsilon - v|^2\,dx$ ，需用 weighted Trudinger–Moser 控制；
(iii) 极限系统剩余项：涉及 $\partial_t\rho + \nabla_x\cdot(\rho v) - \Delta_x\rho$ 等，因 $(\rho,v)$ 满足精确方程而消失。

关键引理（Section 3.2）导出如下 relative entropy inequality：

\frac{d}{dt}\mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}] + \frac{1}{\varepsilon^2}\iint |\nabla_\xi f^\varepsilon + (\xi - \varepsilon v^\varepsilon)f^\varepsilon|^2\,dxd\xi \le C\Big( \mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}] + \|v^\varepsilon - v\|_{L^2}^2 \Big) + \mathcal{O}(\varepsilon^2).

该不等式结合 Grönwall 引理，直接给出 $\mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}](t) \le C\varepsilon^2$ ，再经 Csiszár–Kullback–Pinsker 不等式：

\|f^\varepsilon - M_{\rho,u}\|_{L^1}^2 \le 2\mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}],

即得 $f^\varepsilon$ 的 $$L^1$$ 定量收敛。

对流体部分，作者构造 modulated energy：

\mathcal{E}^\varepsilon(t) := \frac{1}{2}\|v^\varepsilon - v\|_{L^2}^2 + \mathcal{H}[f^\varepsilon \,|\, M_{\rho,u}],

并证明其满足闭合微分不等式：

\frac{d}{dt}\mathcal{E}^\varepsilon(t) + \|\nabla_x(v^\varepsilon - v)\|_{L^2}^2 \le C\mathcal{E}^\varepsilon(t) + \mathcal{O}(\varepsilon^2),

从而获得 $v^\varepsilon$ 的 $L^\infty_t L^2_x \cap L^2_t H^1_x$ 定量收敛。

此外，为处理全空间情形，作者系统使用了 Dunford–Pettis 定理（保证 $f^\varepsilon$ 、 $\rho^\varepsilon$ 、 $\frac{1}{\varepsilon}m^\varepsilon$ 的 $$L^1$$ 一致可积性）与 Aubin–Lions 紧性引理（对 $v^\varepsilon$ 提供 $L^2_{\mathrm{loc}}$ 强紧性），这是弱收敛结论的基础。

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

首要技术难点是 二维全空间中混合项 $\int \rho^\varepsilon |v^\varepsilon|^2 dx$ 的一致控制。标准 Sobolev 嵌入 $H^1(\mathbb{R}^2) \not\hookrightarrow L^\infty$ 失效，而 $\rho^\varepsilon \in L\log L$ 仅提供次指数可积性。作者引入的 weighted Trudinger–Moser 不等式（Lemma 2.1）巧妙绕过此障碍：

\int_0^T\!\!\int_\Omega \rho^\varepsilon |v^\varepsilon|^2\,dxdt \le C\left(1 + \sup_{t}\int_\Omega \rho^\varepsilon|\log\rho^\varepsilon|\,dx\right)\|v^\varepsilon\|_{L^2(0,T;H^1)}^2.

该不等式利用了 $\rho^\varepsilon$ 的熵有界性与 $v^\varepsilon$ 的耗散控制，是闭合所有先验估计的关键。

其次，为提取 $\frac{1}{\varepsilon}\int \xi f^\varepsilon\,d\xi$ 的极限，作者采用 moment expansion + Hilbert expansion 启发式校正（Appendix A）：形式展开 $f^\varepsilon = M_{\rho,0} + \varepsilon f_1 + \mathcal{O}(\varepsilon^2)$ ，代入动量方程可得首阶校正 $$f_1$$ 满足 Poisson 方程，其积分恰好给出 $\rho(v - \nabla_x\log\rho)$ 。这解释了为何 $\varepsilon^2$ 是最优速率——首阶项不消失，但二阶项才贡献收敛误差。

最后，为在 BL 距离下获得 $\varepsilon^2$ 速率，作者发现：密度 $\rho^\varepsilon$ 的首阶 Hilbert 校正恒为零（Remark 1.3），故误差始于 $\varepsilon^2$ 阶，这与 BL 距离对低频模式的敏感性完美匹配，从而实现拓扑最优性。

与之前工作的比较

a) 突破与改进

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

第一，首次建立 light-particle 极限的定量理论。 此前所有工作（Goudon–Jabin–Vasseur [2004], Choi–Kwon [2018], Goudon–Jung [2021]）均仅用 compactness 方法 得到子列收敛，无法给出任何速率。本文通过适配 light-particle 标度的 relative entropy 结构，首次导出显式 $\mathcal{O}(\varepsilon^2)$ 速率，将“存在性”提升至“可计算性”。

第二，彻底摆脱对周期域的依赖。 所有先前定量结果（如 [31,42]）均局限于 $\mathbb{T}^d$ 或有界区域，因其依赖 Poincaré 不等式与紧嵌入。本文通过 weighted Trudinger–Moser 不等式 + 一阶矩控制，首次在 $\mathbb{R}^2$ 上完成全套估计，解决了“无 confinement 下如何控制空间扩散”的开放问题。

第三，实现弱拓扑下的最优速率。 现有文献中，即使对 simpler VFP 方程，BL 距离的定量估计也仅达 $\mathcal{O}(\varepsilon)$ （如 [11]）。本文通过揭示密度首阶校正的消失性，将 BL 速率提升至 $\mathcal{O}(\varepsilon^2)$ ，且证明该速率不可改进（Remark 1.3），达到拓扑与尺度双重最优。

综上，本文不仅填补了 light-particle 极限定量理论的空白，更将 kinetic-fluid 渐近分析推进至统一、最优、全空间的新阶段。

5. 液晶与拓扑涡度：适定弱解的光滑性

Liquid crystals and topological vorticity: smoothness of mild solutions

作者: Fanghua Lin, Yannick Sire, Yantao Wu, Yifu Zhou

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18726v1

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

详细分析

中文标题: 液晶与拓扑涡度：适定弱解的光滑性

研究对象和背景

本文研究一类耦合非局部抛物-输运-几何演化系统，其物理背景是向列型液晶流（nematic liquid crystal flow）在二维平面上的涡度形式建模。核心动机在于：传统液晶模型（如Lin-Liu系统）虽具清晰物理意义，但其速度场 $$u$$ 为向量场，导致分析困难；而通过引入拓扑涡度 $\theta$ （即标量涡度），可将流体动力学部分降维为一个标量拟地转方程（SQG-type equation），同时保留与取向场 $$d$$ 的强几何耦合。该框架不仅揭示了液晶流动中拓扑缺陷（vortices）与取向奇点之间的内在联系，也为处理Navier-Stokes类方程的正则性问题提供新视角。

论文聚焦于如下关键系统（见原文第2页公式(2)及第4页具体化版本(3)）：

\begin{cases} \partial_t \theta + u \cdot \nabla \theta + (-\Delta)^a \theta = (-\Delta)^{-1/2} \nabla^\perp \nabla \cdot \Xi, \\ u = R^\perp \theta = \nabla^\perp (-\Delta)^{-1/2} \theta, \quad \nabla \cdot u = 0, \\ \partial_t d + u \cdot \nabla d = \Delta d + |\nabla d|^2 d, \end{cases}

其中：

$\theta : \mathbb{R}^2 \times (0,T) \to \mathbb{R}$ 是标量涡度（topological vorticity），替代传统速度场 $$u$$ ；
$d : \mathbb{R}^2 \times (0,T) \to \mathbb{S}^2$ 是单位取向矢量场，满足球面约束 $|d| \equiv 1$ ；
$R^\perp = \nabla^\perp (-\Delta)^{-1/2}$ 是二维旋转Riesz变换，将标量涡度映射为无散速度场；
$(-\Delta)^a$ （ $a \in (\frac{1}{2},1)$ ）为分数阶拉普拉斯算子，定义为傅里叶乘子 $|\xi|^{2a}$ 或超奇异积分： $(-\Delta)^a \theta(x) = c(2,a) \, \mathrm{P.V.} \int_{\mathbb{R}^2} \frac{\theta(x)-\theta(y)}{|x-y|^{2+2a}} \, dy;$
应力张量 $\Xi$ 定义为： $\Xi := \nabla d \odot \nabla d - \frac{1}{2} |\nabla d|^2 I,$ 其中 $(\nabla d \odot \nabla d)_{ij} = \partial_i d \cdot \partial_j d$ ，体现取向场的能量密度与各向异性应力；
右端项 $(-\Delta)^{-1/2} \nabla^\perp \nabla \cdot \Xi$ 是奇异性高达一阶导数的非局部源项，本质是 $\Xi$ 的“涡度型”散度，刻画液晶中拓扑缺陷驱动的宏观涡旋生成机制。

该系统定义在全空间 $\mathbb{R}^2$ 上，无边界条件，初始数据为：

\theta_0 \in L^{\frac{2}{2a-1}}(\mathbb{R}^2), \quad \nabla d_0 \in L^2(\mathbb{R}^2).

值得注意的是，文中明确指出不考虑 Dirichlet 或 Neumann 边界（因分数阶算子在边界处引发非平凡奇性），且所有分析均基于整体空间 $\mathbb{R}^2$ —— 这保证了 Riesz 变换的 $$L^p$$ 有界性（ $1<p<\infty$ ）及分数热核的良好衰减性质。

从物理角度看，该模型试图捕捉两个核心现象：

液晶中拓扑涡度的守恒与演化： $\theta$ 不再是经典流体力学中的 $\omega = \nabla^\perp \cdot u$ ，而是被赋予新的几何含义——它与取向场 $$d$$ 的拓扑荷（topological charge） 密切相关， $\nabla^\perp \nabla \cdot \Xi$ 正是描述该荷输运的通量；
分数耗散对奇点形成的抑制效应：当 $a > \frac{1}{2}$ （亚临界情形）， $(-\Delta)^a$ 提供比经典 Laplacian 更强的高频阻尼，有望阻止涡度集中与爆破；而 $a = \frac{1}{2}$ （临界 SQG）则处于正则性阈值，需精细估计。

此外，作者强调该系统与多个重要模型存在深刻联系：

与表面准地转方程（SQG） 直接对应（当 $$d$$ 恒定时退化为纯 SQG）；
与调和映射热流（heat flow of harmonic maps） 耦合，后者本身在二维已知存在有限时间奇点；
与反铁磁体（anti-ferromagnet）模型 类比，其“拓扑涡度”由 $\Omega = m \cdot (\partial_1 m \wedge \partial_2 m)$ 定义，满足类似守恒律 $\partial_t \Omega = -\nabla^\perp \nabla \cdot \widetilde{\Xi}$ ；
与Harmonic-Ricci 流 形式相似，暗示几何分析工具（如共形变换、曲率控制）可能迁移应用。

综上，本研究对象是一个兼具非局部性、几何约束、强非线性与高阶奇异性的耦合系统，其数学挑战在于：如何在分数耗散不足以完全压制非线性输运的前提下，仍能建立全局光滑解。

主要定理或主要结果

论文的核心成果是建立系统(3)的全局光滑性（global smoothness），即在适当初值假设下，任意有限时间内解自动属于 $C^\infty$ 类。该结论以三个递进式定理呈现，构成完整正则性链：

定理 1.1（主定理）
设 $(\theta,d)$ 是系统(3)的一个mild 解（即半群意义下的积分解），且满足以下条件：

初值 $\theta_0 \in L^{\frac{2}{2a-1}}(\mathbb{R}^2)$ ， $\nabla d_0 \in L^2(\mathbb{R}^2)$ ；
$\theta \in L^2 \cap L^{p_0}((0,T)\times\mathbb{R}^2)$ ，其中 $p_0 > \frac{2a+2}{2a-1}$ ；
$|\nabla d| \in L^2 \cap L^{q_0}((0,T)\times\mathbb{R}^2)$ ，其中 $$q_0 > 4$$ ；
则该解是经典的：
$\theta,\, u,\, d \in C^\infty((0,T)\times\mathbb{R}^2).$

该定理的“近最优性”体现在：当 $$a=1$$ （即恢复经典 Laplacian 耗散），条件退化为 $$p_0 > 4$$ ， $$q_0 > 4$$ ，与 Lin-Liu 原始系统中已知的 $$L^4$$ 正则性准则一致（见[LLW10]）。因此，本文将正则性阈值连续延拓至分数阶亚临界区域。

为证明定理1.1，作者分解为三步，每步本身具有独立价值：

定理 1.2（高阶可积性）
在同一定理1.1的假设下，可推出：

\theta,\, u,\, \nabla d \in L^p((0,T)\times\mathbb{R}^2), \quad \forall p \in [2,\infty).

即解的可积性可任意提升至 $$L^p$$ ，为后续 Sobolev 嵌入铺路。

定理 1.3（Hölder 连续性）
在同一定理1.1的假设下，存在 $\alpha > 0$ ，使得：

\theta,\, u \in C^\alpha_{\mathrm{loc}}((0,T)\times\mathbb{R}^2).

此零阶正则性是启动标准Bootstrap 机制的关键跳板。

最终，结合 Schauder 估计、分数阶 De Giorgi-Nash-Moser 理论及调和映射热流的先验估计，即可迭代获得任意阶导数的 Hölder 连续性，从而证得 $C^\infty$ 光滑性。

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

a) 核心分析学工具

分数阶热半群理论与核估计：
定义 $P_{a,t} f(x) = \int_{\mathbb{R}^2} p_a(t,x-y) f(y) \, dy$ ，其中 $$p_a$$ 是 $\partial_t + (-\Delta)^a$ 的基本解。作者引用 [JS20] 给出关键点态估计：
$|\nabla^k p_a|(t,x) \lesssim \delta((0,0),(t,x))^{-2-k}, \quad |\partial_t p_a|(t,x) \lesssim \delta((0,0),(t,x))^{-2},$

其中 $\delta((s,y),(t,x)) = \max\{ |t-s|^{1/(2a)}, |x-y| \}$ 是抛物型距离。此估计支撑了所有卷积型先验估计的收敛性。
Morrey 空间与抛物位势理论（见附录）：
为控制非局部项 $(-\Delta)^{-1/2}\nabla^\perp\nabla\cdot\Xi$ ，作者系统运用抛物 Morrey 空间 $M^{p,\lambda}_{\mathrm{par}}$ 的嵌入性质：若 $f \in M^{p,\lambda}_{\mathrm{par}}$ ，则 $(-\Delta)^{-\beta} f$ 属于更高阶 Morrey 空间。这绕开了 $$L^p$$ 框架下分数积分算子的端点失效问题。
广义Stroock-Varopoulos不等式与分数能量估计：
对 $\theta$ 方程测试 $\theta$ 得能量恒等式：
$\frac{1}{2}\frac{d}{dt}\int |\theta|^2 + \int |(-\Delta)^{a/2}\theta|^2 = \int F \theta,$

其中右端需用分数 Gagliardo-Nirenberg 插值：
$\left|\int F \theta\right| \leq \|F\|_{L^r} \|\theta\|_{L^{r'}} \lesssim \|\nabla d\|_{L^{2r}}^2 \|\theta\|_{L^{r'}},$

再结合 $\|\theta\|_{L^{r'}}$ 的抛物嵌入控制，实现闭环。
调和映射热流的尖锐先验估计：
对 $$d$$ 方程，利用其作为球面值映射的特性，推导出：
$\frac{1}{2}\frac{d}{dt}\int |\nabla d|^2 + \int |\Delta d + |\nabla d|^2 d|^2 = -\int (u\cdot\nabla d)\Delta d.$

右端通过分部积分与 Riesz 变换 $$L^p$$ 有界性转化为 $\|u\|_{L^p} \|\nabla d\|_{L^q} \|\nabla^2 d\|_{L^r}$ ，再用插值吸收。

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

奇异性匹配与平衡律设计：
非局部源项 $F = (-\Delta)^{-1/2}\nabla^\perp\nabla\cdot\Xi$ 的奇性阶为 $$1$$ （即比 $\theta$ 方程主算子 $(-\Delta)^a$ 低 $$2a-1 < 1$$ 阶），直接估计发散。作者通过能量平衡律（10） 强制要求：
$\int F \theta + \int (u\cdot\nabla d)\Delta d = 0,$

将 $$F$$ 的奇异性“转移”至 $u\cdot\nabla d$ 项，并利用 $u = R^\perp \theta$ 的非局部性，将二者统一为 $\Xi$ 的共轭配对，从而在能量层面自动抵消。
尺度临界性分析与参数三分类：
明确区分临界（ $a=\alpha=1/2$ ）、亚临界（ $$a>1/2$$ ）、超临界（ $$a<1/2$$ ）三类。本文仅处理亚临界情形，此时 $(-\Delta)^a$ 的耗散优势足以压制 $$F$$ 的奇性，而临界情形需 Caffarelli-Vasseur 方法，超临界则尚无通用工具。
mild 解框架规避弱解紧性障碍：
不采用传统的 Galerkin 逼近与紧性论证（易受分数阶算子非局部性干扰），而直接在抛物型 $$L^p$$ -空间中构造压缩映射：定义解算子
$\mathcal{T}(\theta,d)(t) = P_{a,t}\theta_0 + \int_0^t P_{a,t-s}\big[ (-\Delta)^{-1/2}\nabla^\perp\nabla\cdot\Xi - u\cdot\nabla\theta \big] ds,$

并利用热核估计证明 $\mathcal{T}$ 在合适加权 $$L^p$$ 空间中为压缩映射，从而获得局部 mild 解存在唯一性。

与之前工作的比较

a) 突破与改进

本文相较 Lin-Liu 原始工作 [Lin89, LLW10] 及 SQG 相关文献 [CMT94, CV10]，实现了三重实质性推进：

正则性假设的维度解耦与分数化：
经典结果要求 $u \in L^4_{t,x}$ （即速度场四次可积），而本文将该条件转化为标量涡度 $\theta$ 的 $L^{p_0}_{t,x}$ 条件，且 $$p_0$$ 可随 $$a$$ 连续变化。当 $a \to 1^-$ ， $p_0 \to 4^+$ ，完美衔接；当 $$a = 0.6$$ ，仅需 $$p_0 > 3.5$$ ，显著弱于 $$L^4$$ 。这是首次在液晶模型中实现耗散阶数与正则性阈值的定量联动。
几何耦合奇性的系统性处理：
以往工作（如 [HL01, KP96]）多单独处理 $$d$$ -方程或 $\theta$ -方程，而本文首创将 $\nabla^\perp\nabla\cdot\Xi$ 视为拓扑涡度的内生源项，并设计匹配的分数阶 Biot-Savart 关系 $u = \nabla^\perp(-\Delta)^{-1+\alpha}\theta$ ，使整个系统的能量结构保持闭合。这种“几何-分析双适配”设计，为同类耦合系统（如磁流体、活性物质）提供了范式。
工具层面的跨领域融合：
首次将抛物 Morrey 空间理论（源自偏微分方程正则性研究）与液晶拓扑缺陷分析（源自材料科学）深度结合。例如，利用 Morrey 嵌入 $M^{p,\lambda} \hookrightarrow C^\alpha$ 直接导出定理1.3，避免了传统 De Giorgi 迭代中繁琐的缩放论证，极大提升了技术透明度。

综上，本文不仅解决了特定模型的正则性问题，更构建了一套面向非局部-几何耦合系统的通用分析框架，其思想与技术有望推广至三维液晶、量子液晶及生物膜动力学等前沿方向。

6. 关于各向异性Boussinesq方程在半空间 $\mathbb{R}_+^3$ 中临近静水平衡态解的衰减性

On decay of solutions to the anisotropic Boussinesq equations near the hydrostatic balance in half space $\mathbb{R}_+^3$

作者: Wangrong Yang, Aibin Zang

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18481v1

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

详细分析

中文标题: 关于各向异性Boussinesq方程在半空间 $\mathbb{R}_+^3$ 中临近静水平衡态解的衰减性

1. 研究对象和背景

本文研究的是三维不可压缩各向异性Boussinesq方程组在上半空间 $\mathbb{R}_+^3 = \{(x_1,x_2,x_3) \in \mathbb{R}^3 \mid x_3 > 0\}$ 上的初边值问题，特别关注其在静水平衡态（hydrostatic balance）附近的小扰动解的全局稳定性与大时间衰减行为。该模型是地球物理流体力学中描述大气与海洋运动的核心方程之一，广泛用于模拟热对流、地转平衡、Ekman层动力学等现象。

论文从标准 Boussinesq 方程出发，其原始形式（见第1页公式 (1.1)）为：

\begin{cases} \partial_t u + u \cdot \nabla u + \nabla P = \nu \Delta u + \Theta e_3, \\ \partial_t \Theta + u \cdot \nabla \Theta = \kappa \Delta \Theta, \\ \nabla \cdot u = 0, \end{cases} \quad x \in \mathbb{R}^3,\; t > 0,

其中 $$u = (u_1,u_2,u_3)$$ 是速度场， $$P$$ 是压力， $\Theta$ 是温度（或密度扰动）， $$e_3 = (0,0,1)$$ 表示重力方向， $\nu \ge 0$ 和 $\kappa \ge 0$ 分别为粘性系数与热扩散系数。当 $\nu = \kappa = 0$ 时系统退化为 Euler-Boussinesq；当二者均为正时，经典理论可建立局部适定性；而当仅部分方向存在耗散时，则称为各向异性耗散（anisotropic dissipation）——这正是本文的关键设定。

具体地，作者考虑的是仅具水平耗散（horizontal dissipation）的情形：即 Laplace 算子被替换为水平拉普拉斯 $\Delta_h = \partial_{x_1}^2 + \partial_{x_2}^2$ ，而垂直方向 $$x_3$$ 上无粘性/热扩散机制。此时方程组变为（见第2页公式 (1.6)）：

\begin{cases} \partial_t u + u \cdot \nabla u + \nabla p = \Delta_h u + \theta e_3, \\ \partial_t \theta + u \cdot \nabla \theta = \Delta_h \theta, \\ \nabla \cdot u = 0, \end{cases} \quad x \in \mathbb{R}_+^3,\; t > 0,

其中扰动变量定义为 $\theta = \Theta - \Theta^{(0)}$ , $p = P - P^{(0)}$ ，对应于静水平衡态

\Theta^{(0)}(x) = x_3,\quad P^{(0)}(x) = \frac{1}{2}x_3^2,\quad u^{(0)} = 0.

该平衡态满足静力学关系 $\partial_{x_3} P^{(0)} = \Theta^{(0)}$ ，并自然满足边界条件。

研究区域 $\mathbb{R}_+^3$ 是一个典型的无界但具有边界的非紧区域，其边界 $\{x_3 = 0\}$ 引入了关键的几何约束。为此，作者施加如下Navier 类滑移边界条件（见 (1.3)–(1.4)）：

u_3 = 0,\quad \partial_{x_3} u_1 = \partial_{x_3} u_2 = 0,\quad \text{on } x_3 = 0,

\theta = 0,\quad \text{on } x_3 = 0.

这些条件既非完全无滑移（Dirichlet），亦非自由滑移（Neumann），而是介于两者之间的一种物理合理的中间情形，反映了流体在固-气/固-水界面处的部分滑移效应。特别值得注意的是，该边界条件导致法向速度分量 $$u_3$$ 为零，但切向导数 $\partial_{x_3} u_h$ 消失，从而保证了水平梯度算子 $\nabla_h$ 在边界上的良好行为。

函数空间方面，初始数据取在 Sobolev 空间 $H^3(\mathbb{R}_+^3)$ 中，并要求满足相容性条件： $\nabla \cdot u_0 = 0$ ， $\partial_{x_3}^2 \theta_0|_{x_3=0} = 0$ 。解的存在性与唯一性被证明在 $L^\infty(0,\infty; H^3(\mathbb{R}_+^3))$ 中成立，且能量控制依赖于初始小性假设 $\|u_0\|_{H^3} + \|\theta_0\|_{H^3} \le \varepsilon$ 。这种高阶正则性要求对于处理非线性项及边界耦合至关重要。

此外，由于缺乏垂直耗散，系统表现出强烈的各向异性结构：水平方向具备强阻尼，而垂直方向仅靠对流与压力梯度维持平衡。这一特性使得经典的谱方法或均匀衰减估计失效，必须发展新的频率分解技术——如后文所述，作者采用水平 Fourier 变换 + 垂直 Fourier 余弦/正弦变换联合策略，以精确刻画不同方向的传播与耗散机制。

综上，本研究的对象是一个兼具物理真实性、数学挑战性与结构特殊性的偏微分方程系统：它处于地球物理建模前沿，又处于现代非线性 PDE 分析的难点交汇处——非紧区域、弱耗散、边界效应、各向异性、强非线性四者交织，构成极具代表性的分析平台。

2. 主要定理或主要结果

本文建立了两个核心定理，分别刻画了解的全局稳定性与最优代数衰减率。

定理 1.1（全局稳定性）

设初始数据 $(u_0,\theta_0) \in H^3(\mathbb{R}_+^3)$ 满足边界条件 (1.3)–(1.4)，且 $\nabla \cdot u_0 = 0$ ， $\partial_{x_3}^2 \theta_0|_{x_3=0} = 0$ 。则存在 $\varepsilon > 0$ ，使得若
$\|u_0\|_{H^3} + \|\theta_0\|_{H^3} \le \varepsilon,$

则初边值问题 (1.6) 存在唯一全局解 $(u,\theta) \in L^\infty(0,\infty; H^3(\mathbb{R}_+^3))$ ，并满足如下一致能量估计：
$$ \sup_{t \ge 0}\left( |u(t)|{H^3}^2 + |\theta(t)|{H^3}^2 \right)
\int_0^\infty \left( |\nabla_h u(\tau)|{H^3}^2 + |\nabla_h \theta(\tau)|{H^3}^2 \right) d\tau \le C \varepsilon^2. $$

该定理表明，在静水平衡附近的小扰动将永远保持在 $$H^3$$ 小邻域内，且水平梯度的能量持续耗散，体现了系统的Lyapunov 稳定性。

定理 1.2（衰减性）

设 $9/10 < \sigma < 1$ ，初始数据 $(u_0,\theta_0)$ 还满足低频控制条件：
$\|\Lambda_h^{-\sigma} u_0\|_{L^2} + \|\Lambda_h^{-\sigma} \theta_0\|_{L^2} \le \varepsilon, \quad \|\partial_{x_3} \Lambda_h^{-\sigma} u_0\|_{L^2} + \|\partial_{x_3} \Lambda_h^{-\sigma} \theta_0\|_{L^2} \le \varepsilon,$

其中 $\Lambda_h^{-\sigma}$ 是水平 Fourier 乘子算子，符号为 $|\xi_h|^{-\sigma}$ 。则对任意 $1-\sigma \le \delta < \sigma - \frac{1}{2}$ ，全局解满足如下衰减估计：
$\|u(t)\|_{L^2} + \|\theta(t)\|_{L^2} \le C \varepsilon (1+t)^{-\frac{\sigma}{2}+\delta},$

$\|\partial_{x_3} u(t)\|_{L^2} + \|\partial_{x_3} \theta(t)\|_{L^2} \le C \varepsilon (1+t)^{-\frac{\sigma}{2}+3\delta},$

$\|\nabla_h u(t)\|_{L^2} + \|\nabla_h \theta(t)\|_{L^2} \le C \varepsilon (1+t)^{-\frac{\sigma}{2}+1+\delta}.$

该结果首次在半空间+Navier 边界+水平耗散框架下，获得了扰动解及其各向异性导数的最优代数衰减率，揭示了水平耗散主导、垂直方向受边界反射调制的深层动力学机制。

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

a) 核心分析学工具

本文综合运用了多种现代非线性 PDE 分析技术，形成一套严密递进的技术链：

Bootstrapping 能量方法：这是证明全局稳定性的主干工具。作者构造带时间权重的能量泛函
$$ E(t) = \sup_{0 \le \tau \le t} \left( |u(\tau)|{H^3}^2 + |\theta(\tau)|{H^3}^2 \right)
- \int_0^t \left( |\nabla_h u|{H^3}^2 + |\nabla_h \theta|{H^3}^2 \right) d\tau, $$
  并通过反复应用 Sobolev 嵌入、插值不等式与乘积估计（引理 2.1–2.2），推导出闭合型不等式 $E(t) \le E(0) + C E^{3/2}(t)$ ，再借助标准 bootstrapping 论证完成控制。
各向异性插值与奇异积分估计：由于系统缺失垂直耗散，传统 $$H^m$$ 能量无法直接控制 $\partial_{x_3} u$ 。作者引入关键引理 2.2，给出如下一维 Gagliardo–Nirenberg 型不等式：
$\|f\|_{L^\infty(\mathbb{R}_+)} \le C \|f\|_{L^2}^{1/2} \|f'\|_{L^2}^{1/2},$

并推广至三维三重乘积估计：
$\left| \int_{\mathbb{R}_+^3} f g h \,dx \right| \le C \|f\|_{L^2}^{1/2} \|\partial_{x_1} f\|_{L^2}^{1/2} \|g\|_{L^2}^{1/2} \|\partial_{x_2} g\|_{L^2}^{1/2} \|h\|_{L^2}^{1/2} \|\partial_{x_3} h\|_{L^2}^{1/2}.$

此类估计精准捕获了各向异性导数间的补偿效应，是控制非线性对流项 $\int u \cdot \nabla u \cdot \Delta u$ 的基石。
Fourier 变换联合技术：为获得衰减率，作者对水平变量 $$x_h = (x_1,x_2)$$ 施行标准 Fourier 变换，对垂直变量 $$x_3$$ 则分别使用 Fourier sine 变换（对应 Dirichlet 条件）与 cosine 变换（对应 Neumann/滑移条件）。由此将原 PDE 转化为含参数 $\xi_h$ 的常微分方程组，进而导出线性化系统的 Green 函数结构。该策略绕开了半空间上 Green 函数难以显式表达的困难，实现了频域层面的“解耦”。
Duhamel 原理与核估计：非线性解由 Duhamel 公式表示为
$(u,\theta)(t) = e^{t\mathcal{L}}(u_0,\theta_0) + \int_0^t e^{(t-s)\mathcal{L}} \mathcal{N}(u,\theta)(s)\,ds,$

其中 $\mathcal{L}$ 是线性算子， $\mathcal{N}$ 是非线性项。作者证明该公式可分解为 七个退化、各向异性的卷积核，每个核在频域中呈现 $|\xi_h|^{-2}$ 或 $|\xi_h|^{-1}$ 奇异性。通过对这些核施行精细的 $$L^1$$ – $$L^2$$ 型估计（利用 Marcinkiewicz 插值与 Mikhlin 乘子定理），最终导出时间衰减幂律。

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

边界项处理：在对高阶导数做能量估计时，边界积分 $\int_{x_3=0} \partial_{x_3}^3 p\, \partial_{x_3}^3 u_3 \,d\sigma$ 显著阻碍估计闭合。作者通过引理 2.4 证明：在 Navier 条件与相容性 $\partial_{x_3}^2 \theta_0|_{x_3=0}=0$ 下，有 $\partial_{x_3}^3 p|_{x_3=0} = 0$ ，从而消去该危险项。此结论依赖于对边界上水平散度与垂直导数关系的精细追踪。
低频控制设计：为获得衰减，必须控制解的低频部分。作者引入广义 Riesz 势 $\Lambda_h^{-\sigma}$ ，并要求初始数据属于负阶 Sobolev 空间 $W^{-\sigma,2}_h$ 。这一设定巧妙规避了因缺失垂直耗散而导致的 $$L^1$$ – $L^\infty$ 估计失效问题，使低频能量可随时间增长而缓慢释放。
非线性项的各向异性分解：对形如 $\partial^\alpha (u \cdot \nabla u)$ 的高阶导数项，作者按导数落在 $$u$$ 或 $\nabla u$ 上进行分类，并利用 $\nabla_h$ 的强耗散优势，将高阶垂直导数 $\partial_{x_3}^k$ 与低阶水平导数 $\nabla_h^l$ 配对，形成“水平主导、垂直辅助”的估计范式。

4. 与之前工作的比较

a) 突破与改进

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

首次在半空间中建立水平耗散 Boussinesq 的全局稳定性：此前工作（如 [10], [15], [20]）主要集中于全空间 $\mathbb{R}^3$ 或周期盒 $\mathbb{R}^2 \times \mathbb{T}$ ，而半空间因边界反射、法向约束与各向异性耦合，长期被视为技术禁区。本文通过 Navier 边界条件与水平–垂直频域分离，成功跨越此障碍。
显著放宽对初始数据的正则性要求：经典结果（如 [2], [6]）通常要求 $$H^s$$ 初值满足 $$s > 5/2$$ 或更高；而本文仅需 $$H^3$$ ，且允许 $\partial_{x_3}^2 \theta_0$ 在边界消失（而非强制为零），更符合物理实际。
获得迄今最优的各向异性衰减率：相比 Wu–Zhang [22] 在 $\mathbb{R}^2 \times \mathbb{T}$ 中得到的 $(1+t)^{-1/2}$ 衰减，本文在更难的半空间中仍达到 $(1+t)^{-\sigma/2 + \delta}$ （ $\sigma \to 1^-$ 时趋近 $(1+t)^{-1/2 +}$ ），且区分了水平梯度与垂直导数的不同衰减速率，揭示了耗散各向异性诱导的解结构分层。

综上，本文不仅填补了半空间各向异性 Boussinesq 理论的空白，更提供了一套可推广至其他退化耗散系统的通用分析范式。

7. Navier-Stokes–Cahn–Hilliard 系统的一类松弛近似方法的合理性证明

Justification of a Relaxation Approximation for the Navier-Stokes-Cahn-Hilliard System

作者: Jan Giesselmann, Jens Keim, Fabio Leotta, Christian Rohde

类别: math.NA, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18463v1

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

详细分析

中文标题: Navier-Stokes–Cahn–Hilliard 系统的一类松弛近似方法的合理性证明

1. 研究对象和背景

本文研究的是不可压缩两相流在扩散界面（diffuse-interface）建模框架下的核心数学模型——Navier–Stokes–Cahn–Hilliard (NSCH) 系统，及其一类低阶松弛近似（relaxation approximation） 的严格数学合理性（justification）。该系统描述两种不可压缩、互不混溶（immiscible）且密度匹配（matched densities） 的粘性流体在恒温条件下的演化动力学，特别适用于涉及拓扑变化（如液滴合并、破裂、Ostwald熟化）的界面流动问题。

论文在第2.1节明确定义了经典 NSCH 系统：设 $d \in \{2,3\}$ ， $\Omega \subset \mathbb{R}^d$ 是具有光滑边界的有界开集（亦可取周期域 $\Omega = \mathbb{T}^d$ ），时间区间为 $$ (0,T) $$ 。NSCH 系统由如下耦合偏微分方程组构成：

\begin{cases} \nabla \cdot u = 0, \\ u_t + (u \cdot \nabla) u + \nabla p = -c \nabla \mu(c) + \nu \Delta u, \\ c_t + \nabla \cdot (c u) = \nabla \cdot \left( \nabla \mu(c) \right), \end{cases} \quad \text{in } \Omega_T := (0,T) \times \Omega.

其中未知函数为：速度场 $u = u(t,x) \in \mathbb{R}^d$ 、压力扰动 $p = p(t,x) \in \mathbb{R}$ 、相场变量 $c = c(t,x) \in \mathbb{R}$ （作为相指示器， $c \approx \pm 1$ 分别对应两相主体， $c \in (-1,1)$ 过渡层宽度由小参数 $\gamma > 0$ 控制）。化学势 $\mu(c)$ 定义为：

\mu(c) := W'(c) - \gamma \Delta c,

而自由能密度 $W : \mathbb{R} \to \mathbb{R}$ 取标准双阱多项式：

W(c) = \frac{1}{4}(c^2 - 1)^2, \quad c \in \mathbb{R}.

该选择保证了热力学一致性：系统耗散总自由能（动能 + Van der Waals 型界面能）。初始条件为：

u(0,\cdot) = u_0, \quad c(0,\cdot) = c_0 \quad \text{in } \Omega.

边界条件采用自然物理设定（当 $\partial \Omega \neq \emptyset$ ）：

u = 0, \quad \nabla c \cdot n = 0, \quad \nabla \mu(c) \cdot n = 0 \quad \text{on } (0,T) \times \partial \Omega.

其中 $$ n $$ 为外法向量；Neumann 条件 $\nabla c \cdot n = 0$ 意味着接触角固定为 $90^\circ$ 。在分析部分（Section 3），作者聚焦于周期边界条件 $\Omega = \mathbb{T}^d$ ，以规避边界处理的技术复杂性，但指出结果可推广至物理边界情形（Remark 3.15）。

NSCH 系统虽被广泛用于多相流模拟，其数值实现却面临严峻挑战：一方面，不可压约束 $\nabla \cdot u = 0$ 引入非局部性；另一方面，Cahn–Hilliard 方程含四阶空间导数 $\Delta^2 c$ （因 $\nabla \cdot \nabla \mu(c) = \Delta W'(c) - \gamma \Delta^2 c$ ），导致刚性与高计算成本，尤其在高速流动（high-velocity flow） 区域。

为此，作者考虑一种全低阶松弛近似系统（Section 2.2），其核心思想是融合三类经典近似：
(i) 人工可压缩性（artificial compressibility）：松弛不可压约束为一阶双曲型演化方程；
(ii) 摩擦型近似（friction-type approximation）：将 Cahn–Hilliard 流量 $j = \nabla \mu(c)$ 视为受阻尼控制的通量变量；
(iii) 筛选泊松松弛（screened Poisson relaxation）：用二阶椭圆算子替代原三阶/四阶结构。

该近似系统引入四维小参数向量 $\varepsilon = (\alpha, \beta, \delta, \tau) \in (0,1)^4$ ，对应不同物理尺度的松弛强度。其方程组为：

\begin{cases} \displaystyle \frac{1}{\alpha} p_t^\varepsilon + \nabla \cdot u^\varepsilon = 0, \\ \displaystyle u_t^\varepsilon + (u^\varepsilon \cdot \nabla) u^\varepsilon + \frac{1}{2} (\nabla \cdot u^\varepsilon) u^\varepsilon + \nabla p^\varepsilon = -c^\varepsilon \nabla \left( W'(c^\varepsilon) + \frac{1}{\beta}(c^\varepsilon - \omega^\varepsilon) \right) + \nu \Delta u^\varepsilon + (\lambda + \nu) \nabla (\nabla \cdot u^\varepsilon), \\ \displaystyle c_t^\varepsilon + \nabla \cdot (c^\varepsilon u^\varepsilon) + \nabla \cdot v^\varepsilon = \tau \Delta c^\varepsilon, \\ \displaystyle \frac{1}{\delta} v_t^\varepsilon + \nabla \left( W'(c^\varepsilon) + \frac{1}{\beta}(c^\varepsilon - \omega^\varepsilon) \right) = -\frac{1}{\delta} v^\varepsilon - \frac{1}{\delta \beta} \nabla \omega^\varepsilon + \nabla (\nabla \cdot v^\varepsilon), \\ \displaystyle -\gamma \Delta \omega^\varepsilon + \frac{1}{\beta} \omega^\varepsilon = \frac{1}{\beta} c^\varepsilon. \end{cases}

未知函数新增：Cahn–Hilliard 通量 $v^\varepsilon$ 和松弛变量 $\omega^\varepsilon$ 。关键改进在于作者引入三项额外耗散项： $(\lambda + \nu) \nabla (\nabla \cdot u^\varepsilon)$ 、 $\tau \Delta c^\varepsilon$ 、 $\nabla (\nabla \cdot v^\varepsilon)$ ，它们在 $|\varepsilon| \to 0$ 时 vanish，但对能量估计至关重要。

函数空间设定上，NSCH 经典解要求 $u \in C^1((0,T); C(\Omega)^d) \cap C^0([0,T]; C^2(\Omega)^d)$ 、 $c \in C^1((0,T); C(\Omega)) \cap C^0([0,T]; C^4(\Omega))$ ；而松弛系统弱解（Definition 3.2）仅需 $u \in L^2(0,T; H^1(\Omega)^d)$ 、 $c \in L^2(0,T; H^1(\Omega)) \cap L^\infty(0,T; L^4(\Omega))$ 、 $v \in L^2(0,T; H_{\mathrm{div}}(\Omega)^d)$ 、 $\omega \in H^1((0,T) \times \Omega)$ ，显著降低正则性要求，体现“低阶”本质。

2. 主要定理或主要结果

本文的核心理论成果是 Theorem 3.13，它建立了松弛近似系统解向 NSCH 系统解的强收敛性，并给出定量收敛速率。

Theorem 3.13 (Convergence of the relaxation approximation)
Let $\Omega = \mathbb{T}^d$ , $$ T > 0 $$ , and assume that the NSCH system (2.1)–(2.4) admits a classical solution $$ (p,u,c) $$ satisfying
$p \in C^0(\overline{\Omega_T}),\; u \in C^1((0,T); C(\Omega)^d) \cap C^0([0,T]; C^2(\Omega)^d),\; c \in C^1((0,T); C(\Omega)) \cap C^0([0,T]; C^4(\Omega)).$

Let $U^\varepsilon = (p^\varepsilon, u^\varepsilon, c^\varepsilon, v^\varepsilon, \omega^\varepsilon)$ be a weak solution to the relaxation system (2.7)–(2.9) in the sense of Definition 3.2, with initial data satisfying
$\|p^\varepsilon(0,\cdot) - p(0,\cdot)\|_{L^2}^2 + \|u^\varepsilon(0,\cdot) - u(0,\cdot)\|_{L^2}^2 + \|c^\varepsilon(0,\cdot) - c(0,\cdot)\|_{L^2}^2 + \|v^\varepsilon(0,\cdot)\|_{L^2}^2 + \|\omega^\varepsilon(0,\cdot) - c(0,\cdot)\|_{H^1}^2 \leq C_0 |\varepsilon|^2$

for some constant $$ C_0 > 0 $$ .
Then, there exists a constant $$ C > 0 $$ , independent of $\varepsilon$ , such that for all $t \in [0,T]$ ,
$\|u^\varepsilon(t,\cdot) - u(t,\cdot)\|_{L^2}^2 + \|c^\varepsilon(t,\cdot) - c(t,\cdot)\|_{L^2}^2 \leq C |\varepsilon|^2.$

In particular, as $|\varepsilon| \to 0$ , $u^\varepsilon \to u$ and $c^\varepsilon \to c$ strongly in $L^\infty(0,T; L^2(\Omega))$ .

该定理的关键在于：

无需粘性假设：结论对 $\nu = 0$ （即 Euler–Cahn–Hilliard 极限）同样成立；
定量误差估计：误差以 $|\varepsilon|^2$ 阶收敛，明确依赖所有四个松弛参数；
弱–强唯一性框架：利用相对熵方法，将松弛解与给定的强 NSCH 解比较，避免对松弛系统自身解的存在性/唯一性作全局假设。

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

a) 核心分析学工具

本文证明的基石是 相对熵方法（relative entropy method），这是处理奇异性极限（singular limits）、弱–强唯一性（weak–strong uniqueness）及模型约化合理性（model reduction justification）的强有力工具。其核心思想是构造一个非负的、关于两个解差的泛函（即相对熵），并证明其时间演化满足 Gronwall 不等式，从而控制误差增长。

作者定义了两类关键能量泛函：

原始相对能量（full relative energy） $\eta_\varepsilon(U,U^\star)$ （公式 (3.8)），包含所有变量差的二次型及 $$ W(c) $$ 的非凸部分；
约化相对能量（reduced relative energy） $\widetilde{\eta}_\varepsilon$ （公式 (3.9)），通过添加 $\frac{1}{2}(c-c^\star)^2$ 并改写 $$ W $$ 为四次主导项（ $\frac{1}{4}c^4$ ），确保其全局非负性，为 Gronwall 引理铺平道路。

证明中至关重要的技术是 Gronwall-type argument on $\widetilde{\eta}_\varepsilon$ ，这要求精确控制 $\frac{d}{dt} \widetilde{\eta}_\varepsilon$ 。由于 $\widetilde{\eta}_\varepsilon = \eta_\varepsilon + \frac{1}{2}(c-c^\star)^2$ ，作者分别估计二者的时间导数（Lemma 3.4 与 Proposition 3.7）。

Lemma 3.4 给出了 $\eta_\varepsilon$ 的积分不等式，其右端包含大量非线性项（如 $(c-c^\star)^3$ , $(c-c^\star)^4$ , $\nabla(c-c^\star) \cdot \nabla c^\star$ 等）。处理这些项依赖于：
- Taylor 展开与恒等式技巧：例如，利用 $c^4 - \bar{c}^4 = (c^3 + c^2 \bar{c} + c \bar{c}^2 + \bar{c}^3)(c-\bar{c})$ （公式 (3.11)）分解高阶差；
- 分部积分与散度形式提取：将许多高阶项重写为 $\nabla \cdot (\cdots)$ ，使其在积分后转化为边界项（周期域下自动消失）或可控余项；
- Young 不等式与 Sobolev 嵌入：对立方项 $(c-c^\star)^3$ ，利用 $|ab| \leq \frac{1}{2}a^2 + \frac{1}{2}b^2$ 及 $\|c-c^\star\|_{L^6} \lesssim \|c-c^\star\|_{H^1}$ 控制（Remark 3.5）。
Proposition 3.7 单独处理 $(c-c^\star)^2$ 项，因其在弱解框架下不直接由能量不等式 (3.4) 控制。作者巧妙地将 $c^\varepsilon$ 的方程与 $c^\star$ 相减，再乘以 $(c^\varepsilon - c^\star)$ 并积分，得到其演化方程，其中关键项 $-\tau |\nabla(c^\varepsilon - c^\star)|^2$ 提供所需耗散。

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

首要难题是 $$ W(c) $$ 的非凸性。标准相对熵要求凸熵，但 $W(c) = \frac{1}{4}(c^2-1)^2$ 在 $$ c=0 $$ 处凹，导致 $\eta_\varepsilon$ 不定号，无法直接 Gronwall。作者的突破性技巧是：

构造“约化”相对能量 $\widetilde{\eta}_\varepsilon$ ，用凸的 $\frac{1}{4}c^4$ 替代 $$ W(c) $$ ，并添加 $\frac{1}{2}(c-c^\star)^2$ （Proposition 3.3 证明其非负性）；
利用恒等式 $\widetilde{\eta}_\varepsilon = \eta_\varepsilon + \frac{1}{2}(c-c^\star)^2$ ，将问题拆解为两个可估计的部分。

第二大难题是 弱解框架下高阶非线性项的处理。例如， $(c^\varepsilon)^3$ 或 $\nabla \omega^\varepsilon \otimes \nabla \omega^\varepsilon$ 出现在估计中。作者采用：

椭圆耦合技巧：利用 $\omega^\varepsilon$ 满足的 Poisson 型方程 $-\gamma \Delta \omega^\varepsilon + \frac{1}{\beta} \omega^\varepsilon = \frac{1}{\beta} c^\varepsilon$ ，将 $\omega^\varepsilon$ 的高阶导数用 $c^\varepsilon$ 表达（公式 (3.32)）；
向量恒等式 (3.33)：对形如 $\varphi \Phi \cdot \nabla \Delta \varphi$ 的项，展开为散度项加低阶项，避免直接估计四阶导数。

第三大难题是 初始数据的相容性。松弛系统初始值 $v^\varepsilon_0$ 、 $\omega^\varepsilon_0$ 需与 NSCH 解 $$ (u,c) $$ 匹配。作者在 Proposition 3.12 中给出了非平凡的初始层构造：令 $v^\varepsilon_0 = \nabla \mu(c_0)$ ， $\omega^\varepsilon_0$ 为 $-\gamma \Delta \omega + \frac{1}{\beta}\omega = \frac{1}{\beta}c_0$ 的解，确保初始误差为 $O(|\varepsilon|)$ ，而非 $$ O(1) $$ （Remark 3.11 强调此构造对获得最优 $|\varepsilon|^2$ 速率的必要性）。

4. 与之前工作的比较

a) 突破与改进

本文的创新性体现在首个对全低阶、双曲–椭圆混合型松弛模型（含独立四参数）向 NSCH 系统的严格奇异性极限分析。此前工作存在以下局限：

[30]（原始模型提出者）：仅进行形式渐近分析，未提供收敛性证明；其模型缺少本文引入的三项耗散（ $\lambda$ , $\tau$ , $\nabla(\nabla \cdot v)$ ），导致能量估计中关键负项缺失，无法闭合 Gronwall 不等式。
[13, 18]（其他松弛方案）：或采用两阶段近似（先摩擦后 Poisson），或仅松弛 Cahn–Hilliard 部分，未同时处理 NSCH 的全部刚性源项（不可压 + 四阶）。本文的单一方程组框架更简洁，且证明覆盖整个系统耦合效应。
相对熵应用的拓展：以往在 Cahn–Hilliard 相关极限中（如 [10, 11, 20]），相对熵仅用于凸自由能情形，或需强解假设。本文首次在非凸双阱 $$ W $$ 下，通过构造 $\widetilde{\eta}_\varepsilon$ ，成功将相对熵方法推广至 NSCH 松弛极限，且不依赖粘性 $\nu > 0$ （Theorem 3.13 对 $\nu = 0$ 成立），而多数现有结果（如 [23] 的适定性）要求 $\nu > 0$ 。
解的正则性要求大幅降低：先前的模型合理性证明（如人工可压缩性对 NS 的极限 [40]）常要求 NS 解属 $$ C^3 $$ 或更高。本文仅需 NSCH 解满足 $c \in C^0_tC^4_x$ 、 $u \in C^0_tC^2_x$ ，这是物理上可实现的正则性（如数值验证中的 Ostwald 熟化解）。
数值验证超越理论假设：Section 4 的数值实验不仅确认了 $O(|\varepsilon|^2)$ 收敛率，更在理论未覆盖的高速、大变形流动（如气泡合并、 $|u| \sim \mathcal{O}(1)$ ）中展示了松弛模型的鲁棒性，表明其实用价值远超当前理论保证范围。

综上，本文不仅填补了 NSCH 松弛近似的严格数学基础空白，更通过精巧的能量构造与估计技巧，为多相流模型降阶与高效数值算法的设计提供了坚实的理论支撑与新范式。

8. 一种用于四维血流磁共振成像衍生流动问题的稳定化有限元方法

A stabilized finite element method for a flow problem arising from 4D flow magnetic resonance imaging

作者: Gabriel Barrenechea, Cristian Cárcamo, Abner Poza

类别: math.NA, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18454v1

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

详细分析

中文标题: 一种用于四维血流磁共振成像衍生流动问题的稳定化有限元方法

1. 研究对象和背景

本文研究的核心问题是：如何从含噪声的4D Flow MRI观测速度场 $\mathbf{u}_m$ 出发，稳健地估计其与真实物理流场之间的观测误差（即“噪声”） $\mathbf{w}$ ，并同步重构不可直接测量的压力场 $$ p $$ 。该问题源于心血管疾病的无创诊断需求——临床亟需在不进行侵入性导管测压的前提下，从非侵入式、但固有含噪的MRI数据中提取可靠的血流动力学生物标志物（如压力梯度、壁面剪应力、涡量等）。因此，研究对象并非原始Navier–Stokes方程本身，而是其一个结构化的扰动模型，该模型将观测误差 $\mathbf{w}$ 视为待求未知量，并将已知的MRI重建速度 $\mathbf{u}_m$ 视为输入数据。

论文从经典的不可压缩Navier–Stokes方程出发：

\begin{cases} \rho \mathbf{u}_t - \mu \Delta \mathbf{u} + \rho (\nabla \mathbf{u}) \mathbf{u} + \nabla p = \mathbf{0} & \text{in } \Omega \times (0,T], \\ \nabla \cdot \mathbf{u} = 0 & \text{in } \Omega \times (0,T], \\ \mathbf{u} = \mathbf{0} & \text{on } \partial\Omega \times (0,T], \\ \mathbf{u}(0) = \mathbf{u}_0 & \text{in } \Omega, \end{cases}

其中 $\Omega \subset \mathbb{R}^d$ （ $$ d=2 $$ 或 $$ 3 $$ ）是一个有界、多面体、Lipschitz连续边界的开区域，代表血管腔； $\rho > 0$ 为血液密度， $\mu > 0$ 为运动粘度； $\mathbf{u}$ 为真实（未知）速度场， $$ p $$ 为对应压力。关键建模步骤在于对速度场进行可解-不可解分解：设 $\mathbf{u} = \mathbf{u}_m + \mathbf{w}$ ，其中 $\mathbf{u}_m$ 是由MRI数据插值得到的已知函数（视为给定数据），而 $\mathbf{w}$ 是待求的观测误差/建模残差。这一分解基于三条合理假设：

(H1) $\mathbf{u} = \mathbf{u}_m + \mathbf{w}$ 在 $\Omega$ 内成立；
(H2) $\nabla \cdot \mathbf{w} = -\nabla \cdot \mathbf{u}_m =: g$ ，即 $\mathbf{w}$ 的散度由 $\mathbf{u}_m$ 的非零散度（因MRI重建不满足精确不可压性）决定；
(H3) $\mathbf{w} = \mathbf{0}$ 在 $\partial\Omega$ 上，即误差在血管壁处消失（符合物理直觉）。

将上述分解代入Navier–Stokes方程，并忽略高阶非线性项 $(\nabla \mathbf{w}) \mathbf{w}$ （或通过线性化处理），即可导出以 $(\mathbf{w}, p)$ 为未知量的修正Oseen型方程组（即文中的模型问题 $(\mathbf{P})$ ）：

\begin{cases} \sigma \mathbf{w} - \mu \Delta \mathbf{w} + \rho (\nabla \mathbf{u}_m) \mathbf{w} + \rho (\nabla \mathbf{w})(\mathbf{a} + \mathbf{u}_m) + \nabla p = \mathbf{f} & \text{in } \Omega, \\ \nabla \cdot \mathbf{w} = g & \text{in } \Omega, \\ \mathbf{w} = \mathbf{0} & \text{on } \partial\Omega. \end{cases}

此处 $\sigma = \rho / \tau > 0$ 是时间离散化引入的正则化参数（ $\tau$ 为固定时间步长）， $\mathbf{a} \in W^{1,\infty}(\Omega)^d$ 是一个给定的无散度对流场（例如，在Newton迭代中可取为前一步的 $\mathbf{w}$ 近似）， $\mathbf{f} \in L^2(\Omega)^d$ 和 $g \in L^2(\Omega)$ 是由 $\mathbf{u}_m$ 及其时间导数生成的源项。该系统本质上是一个带强对流项和源项的奇异摄动Oseen问题，其数学挑战在于：一方面， $\mathbf{u}_m$ 通常仅具低正则性（如分片常数或分片线性），导致系数 $\nabla \mathbf{u}_m$ 在 $L^\infty$ 意义下可能无界；另一方面，临床实践要求采用等阶有限元空间（如 $\mathbb{P}_k-\mathbb{P}_k$ ）同时逼近 $\mathbf{w}$ 和 $$ p $$ ，而这直接违反了标准的Babuška–Brezzi inf-sup稳定性条件，使得经典Galerkin方法不稳定。

为此，作者定义了适配的函数空间：设 $H_0^1(\Omega)^d$ 为零边值Sobolev空间， $L_0^2(\Omega) := \{ q \in L^2(\Omega) : \int_\Omega q = 0 \}$ 为零均值压力空间；其变分形式在 $H \times Q$ 上表述，其中 $H := H_0^1(\Omega)^d$ ， $Q := L_0^2(\Omega)$ 。该设定确保了弱解的适定性框架，且与物理约束（不可压性、无滑移边界）严格一致。

2. 主要定理或主要结果

Theorem 1（连续问题的适定性）
假设 $\mathbf{u}_m \in W^{1,\infty}(\Omega)^d$ 且满足

\sigma > 4\rho \|\nabla \mathbf{u}_m\|_{\infty,\Omega},

则模型问题 $(\mathbf{P})$ 的变分形式 (3.3)–(3.4) 在 $H \times Q$ 中存在唯一解 $(\mathbf{w}, p)$ 。

Lemma 2（离散格式的稳定性）
设 $\sigma > 4\rho \|\nabla \mathbf{u}_m\|_{\infty,\Omega}$ 成立，且稳定化参数 $\delta$ 满足

\delta < \min\left\{ \frac{1}{4}, \frac{1}{2 C_{\text{inv}}^2} \right\},

其中 $C_{\text{inv}}$ 是逆不等式常数。则稳定化双线性型 $B_{\text{stab}}(\cdot,\cdot)$ 满足强制性估计：

B_{\text{stab}}(\mathbf{w}_h,p_h;\mathbf{w}_h,p_h) \ge \frac{1}{4} |||(\mathbf{w}_h,p_h)|||^2, \quad \forall (\mathbf{w}_h,p_h) \in H_h \times Q_h,

其中离散范数定义为

|||(\mathbf{v}_h,q_h)|||^2 := \sigma \|\mathbf{v}_h\|_{0,\Omega}^2 + \mu |\mathbf{v}_h|_{1,\Omega}^2 + \lambda \|\nabla \cdot \mathbf{v}_h\|_{0,\Omega}^2 + \sum_{T \in \mathcal{T}_h} \tau_T \|L_h(\mathbf{v}_h,q_h)\|_{0,T}^2.

Lemma 3（一致性误差估计）
设 $(\mathbf{w},p)$ 与 $(\mathbf{w}_h,p_h)$ 分别为连续问题 $(\mathbf{P})$ 和离散问题 (4.9) 的解，且 $\mathbf{a}_h$ 对 $\mathbf{a}$ 的逼近满足 (5.24)。则对任意 $(\mathbf{v}_h,q_h) \in H_h \times Q_h$ ，有

B_{\text{stab}}(\mathbf{w}-\mathbf{w}_h,p-p_h;\mathbf{v}_h,q_h) \le C \sigma^{-1/2} \rho h^k \|\mathbf{a}\|_{k+1,\Omega} \|\mathbf{w}\|_{2,\Omega} \, |||(\mathbf{v}_h,q_h)|||.

Corollary（最优误差估计，隐含于分析中）
在上述假设下，存在常数 $$ C > 0 $$ ，使得

\sigma^{1/2} \|\mathbf{w} - \mathbf{w}_h\|_{0,\Omega} + \mu^{1/2} |\mathbf{w} - \mathbf{w}_h|_{1,\Omega} + \|p - p_h\|_{0,\Omega} \le C h^k \left( \|\mathbf{w}\|_{k+1,\Omega} + \|p\|_{k,\Omega} \right).

这表明所提方法对速度和压力均达到最优收敛阶 $\mathcal{O}(h^k)$ 。

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

a) 核心分析学工具

本文的理论分析建立在变分方法与Babuška–Brezzi理论的坚实基础上。首先，通过引入加权内积与修正的双线性型 $\tilde{a}(\cdot,\cdot)$ ，作者证明了其在 $$ H $$ 上的椭圆性（coercivity），这是利用Lax–Milgram定理保证解存在唯一的前提。关键技巧在于巧妙运用向量恒等式引理（Lemma 1）：

((\nabla \mathbf{w}) \boldsymbol{\alpha}, \mathbf{v}) = -((\nabla \mathbf{v}) \boldsymbol{\alpha}, \mathbf{w}) - ((\nabla \cdot \boldsymbol{\alpha}) \mathbf{w}, \mathbf{v}),

该恒等式将非对称对流项 $((\nabla \mathbf{w}) \mathbf{a}, \mathbf{v})$ 转化为对称部分与一阶散度项之和，从而允许控制其大小。结合 $\mathbf{a}$ 的无散度性（ $\nabla \cdot \mathbf{a} = 0$ ）及小参数 $\sigma$ 的主导作用，最终导出 $\tilde{a}(\mathbf{v},\mathbf{v}) \gtrsim \sigma \|\mathbf{v}\|^2 + \mu |\mathbf{v}|^2$ 。

其次，针对压力-速度耦合项 $b(\mathbf{v},q) = -(q, \nabla \cdot \mathbf{v})$ ，作者直接援引经典结论：该双线性型在 $H \times Q$ 上满足inf-sup条件（Ladyzhenskaya–Babuška–Brezzi condition），即存在 $\beta > 0$ 使得

\sup_{\mathbf{v} \in H} \frac{b(\mathbf{v},q)}{\|\mathbf{v}\|_{1,\Omega}} \ge \beta \|q\|_{0,\Omega}, \quad \forall q \in Q.

这保证了压力变量的稳定性，是应用Babuška–Brezzi理论的关键。

对于离散分析，作者构建了一个定制化的加权范数 $|||\cdot|||$ ，它不仅包含标准的 $$ H^1 $$ -速度与 $$ L^2 $$ -压力项，还显式嵌入了稳定化项的局部 $$ L^2 $$ -贡献 $\sum_T \tau_T \|L_h(\mathbf{v}_h,q_h)\|_{0,T}^2$ 。该范数的设计使得稳定化项 $S_{\text{press}}$ 能被精确吸收进主能量估计，而无需依赖传统 $$ L^2 $$ -投影或Riesz表示定理。

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

首要技术难点是等阶元的inf-sup失效。作者采用Galerkin Least-Squares (GLS)-型稳定化，但进行了重要改造：其压力稳定项 $S_{\text{press}}$ 并非作用于残差 $\nabla \cdot \mathbf{w}_h - g$ ，而是作用于整个微分算子残差 $\sigma \mathbf{w}_h - \mu \Delta \mathbf{w}_h + L_h(\mathbf{w}_h,p_h) - \mathbf{f}$ 。这种设计使稳定项能同时控制对流、扩散与压力梯度，显著增强鲁棒性。

其次，针对非一致性（inconsistency） ——由于用离散场 $\mathbf{a}_h$ 代替连续场 $\mathbf{a}$ ，导致 $B_{\text{stab}}$ 不再满足 $B_{\text{stab}}(\mathbf{w},p; \mathbf{v}_h,q_h) = F_{\text{stab}}(\mathbf{v}_h,q_h)$ ——作者发展了一套精细的一致性误差分解。通过反复应用恒等式 (3.7) 与逆不等式 (5.18)，将所有 $\mathbf{a} - \mathbf{a}_h$ 相关项统一上界为 $\mathcal{O}(h^k)$ 阶，从而证明误差受制于最佳逼近误差，而非被稳定化项污染。

最后，为处理低正则性数据 $\mathbf{u}_m$ （如分片常数），作者在稳定参数 $\tau_T$ 的设计中引入尺度律

\tau_T = \frac{\delta h_T^2}{\sigma h_T^2 + \mu},

该形式自动在对流主导区（ $\sigma h_T^2 \gg \mu$ ）退化为 $\mathcal{O}(h_T^2/\sigma)$ ，在扩散主导区（ $\mu \gg \sigma h_T^2$ ）退化为 $\mathcal{O}(h_T^2/\mu)$ ，实现了对多尺度物理的自适应响应。

4. 与之前工作的比较

a) 突破与改进

本文相较于既有工作（如 [21], [5], [18]）实现了三重实质性推进：

第一，问题建模的临床针对性更强。 以往MRI驱动的流体力学反演多聚焦于直接压力泊松求解（如基于 $\nabla p = -\rho \mathbf{u}_t - \rho (\mathbf{u} \cdot \nabla) \mathbf{u} + \mu \Delta \mathbf{u}$ ），但该途径严重依赖 $\mathbf{u}_m$ 的高阶导数，对噪声极度敏感。本文则构建了一个物理一致的、以误差 $\mathbf{w}$ 为基本未知量的正向模型，天然内嵌质量守恒修正（通过 $g = -\nabla \cdot \mathbf{u}_m$ ），从根本上缓解了由重建伪影导致的系统性偏差。

第二，稳定化机制更具普适性与鲁棒性。 文献 [5] 提出的GLS格式要求 $\mathbf{a}$ 精确已知且光滑，而本文明确处理了 $\mathbf{a}_h$ 的离散逼近误差，并通过引理3给出了其对整体误差的定量影响，这是此前工作未覆盖的。此外，所设计的 $\tau_T$ 尺度律比经典SUPG或CIP方法更契合4D flow MRI中常见的低雷诺数、瞬态脉动流特性。

第三，理论保证更完整。 [21] 虽提出类似等阶元方案，但未给出最优阶误差估计；本文不仅证明了离散解的存在唯一性（Lemma 2），更通过一致性误差控制（Lemma 3）与Strang引理框架，首次为该类MRI衍生Oseen问题建立了完整的先验最优收敛理论。这为算法在临床数据上的可信度提供了坚实的数学背书。

综上，本文不是对现有数值方法的简单移植，而是一次面向医学影像物理特性的、从建模、离散到分析的全链条创新，显著提升了4D flow MRI后处理的精度、稳定性与理论严谨性。

9. 主动湍流中Tonner-Tu-Swift-Hohenberg方程的整体吸引子及其性质

The global attractor of the Toner-Tu-Swift-Hohenberg equations of active turbulence and its properties

作者: Daniel W. Boutros, Kolluru Venkata Kiran, John D. Gibbon, Rahul Pandit

类别: physics.flu-dyn, cond-mat.soft, math.AP, math.DS, nlin.CD

发布日期: 2026-01-26

arXiv ID: 2601.18233v1

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

详细分析

中文标题: 主动湍流中Tonner-Tu-Swift-Hohenberg方程的整体吸引子及其性质

研究对象和背景

本文研究的核心对象是 Toner-Tu-Swift-Hohenberg (TTSH) 方程组, 一类用于建模主动湍流（active turbulence）——特别是悬浮液中细菌群集（bacterial swarming）——的非线性偏微分方程系统。该方程融合了不可压 Navier-Stokes、Toner-Tu 和 Swift-Hohenberg 方程的关键特征，并引入了两个物理上至关重要的机制：线性驱动（linear driving）与负拉普拉斯扩散叠加超耗散（negative Laplacian diffusion combined with hyper-dissipation）。这种结构使其能自然刻画主动物质中由能量注入引发的自发涡旋生成与时空混沌行为。

在周期域 $\mathbb{T}^d$ （ $$d = 2,3$$ ）上，TTSH 方程组写作：

\begin{cases} \partial_t u + (u \cdot \nabla) u + \Gamma_0 \Delta u + \Gamma_2 \Delta^2 u + \nabla p = -\left( \alpha + \beta |u|^2 \right) u, \\ \nabla \cdot u = 0, \end{cases}

其中 $u: \mathbb{T}^d \times [0,T] \to \mathbb{R}^d$ 是粗粒化速度场， $p: \mathbb{T}^d \times [0,T] \to \mathbb{R}$ 是压力， $\Gamma_0 > 0$ 、 $\Gamma_2 > 0$ 、 $\beta > 0$ 是正参数，而 $\alpha \in \mathbb{R}$ 可正可负（实验中常取 $\alpha < 0$ ，对应能量注入主导）。注意：原文中将材料导数简化为 $\partial_t + u \cdot \nabla$ （即设 Toner-Tu 参数 $\lambda = 1$ ），这是对群体运动中“推-拉”效应的简化处理。

该系统的数学区域与函数空间设定具有关键意义。全文工作于单位周期立方体 $\mathbb{T}^d = [0,1]^d$ 上，采用无散度约束（ $\nabla \cdot u = 0$ ），并始终在 Sobolev 空间框架下分析解的适定性与长期行为。特别地，初始数据取自 $L^2(\mathbb{T}^d)$ ，而解的空间正则性被提升至 $H^2(\mathbb{T}^d)$ （甚至更高），这直接源于 $\Delta^2 u$ 项提供的四阶超耗散（hyper-dissipation），其强度远超标准粘性 $\Delta u$ 。这种高阶耗散不仅确保解的光滑性，更从根本上控制了解的长时间演化，使得动力系统具有有限维整体吸引子（finite-dimensional global attractor）成为可能。

TTSH 方程的物理动机源于对细菌悬浮液中观测到的无序涡旋态（disordered vortex states）的理论刻画。线性稳定性分析揭示了一个核心尺度：Swift-Hohenberg 尺度

\eta_{\text{sh}} := \sqrt{\frac{\Gamma_2}{\Gamma_0}},

它由反扩散项 $+\Gamma_0 \Delta u$ 与超耗散项 $+\Gamma_2 \Delta^2 u$ 的竞争决定。具体而言，围绕零解 $$u=0$$ 的线性化给出增长率为

\Sigma(k) = -\alpha + \Gamma_0 |k|^2 - \Gamma_2 |k|^4,

其最大正值出现在波数 $|k_c| \approx \sqrt{\Gamma_0 / \Gamma_2} = \eta_{\text{sh}}^{-1}$ 处。因此， $\eta_{\text{sh}}$ 成为系统中最不稳定的涡旋特征长度，亦即数值与实验中反复观测到的涡旋直径尺度。这一现象在细菌湍流中已被广泛验证，但此前缺乏严格的数学基础。本文正是要为该尺度的主导性（predominance）提供rigorous theoretical foundation。

此外，TTSH 方程还具备丰富的非线性结构：非线性项 $-(\alpha + \beta|u|^2)u$ 包含线性驱动（ $-\alpha u$ ）与Landau 型饱和非线性（ $-\beta|u|^2 u$ ），前者在 $\alpha < 0$ 时持续注入能量，后者则抑制无限增长；而对流项 $(u \cdot \nabla)u$ 与高阶耗散共同作用，形成类似经典湍流的能量级串（energy cascade），但在二维情形下表现为逆向级串（inverse cascade），即能量向大尺度输运，导致长程速度关联——这与实验中观察到的 $O(10\,\eta_{\text{sh}})$ 尺度的大涡结构一致。

综上，TTSH 方程是一个兼具物理真实性（描述主动湍流）、数学挑战性（高阶非线性、强耦合、无散度约束）与结构丰富性（线性不稳定性、多尺度竞争、能量注入/耗散平衡）的理想模型。本文将其置于严格的无穷维动力系统框架下，聚焦于其长期渐近行为的本质维度，即整体吸引子的存在性、有限维性及其Lyapunov 维数（ $d_L(\mathcal{A})$ ）的定量估计。

主要定理或主要结果

本文的核心数学成果由以下三个定理构成，分别确立了解的全局适定性、整体吸引子的存在性，以及其维数的显式上界估计。

Theorem 1（全局适定性）.
设 $$d = 2$$ 或 $$3$$ ，初始速度场 $u_0 \in L^2(\mathbb{T}^d)$ 。则 TTSH 方程组存在唯一的全局弱解

u \in C([0,T]; L^2(\mathbb{T}^d)) \cap L^2((0,T); H^2(\mathbb{T}^d)) \cap L^4((0,T); L^4(\mathbb{T}^d)),

对任意 $$T > 0$$ 成立。进一步，若 $u_0 \in H^2(\mathbb{T}^d)$ ，则解为强解且属于 $C([0,T]; H^2(\mathbb{T}^d))$ 。

Theorem 2（整体吸引子存在性）.
TTSH 动力系统在相空间 $H = \{ u \in L^2(\mathbb{T}^d) : \nabla \cdot u = 0 \}$ 上存在一个紧致、连通、有限维的整体吸引子 $\mathcal{A} \subset H$ 。即：
(i) $\mathcal{A}$ 是不变集（ $S(t)\mathcal{A} = \mathcal{A},\ \forall t \geq 0$ ）；
(ii) 对任意有界集 $B \subset H$ ，有 $\mathrm{dist}(S(t)B, \mathcal{A}) \to 0$ 当 $t \to \infty$ ；
(iii) $\dim_H(\mathcal{A}) = \dim_F(\mathcal{A}) = d_L(\mathcal{A}) < \infty$ ，其中 $\dim_H$ 为 Hausdorff 维数， $\dim_F$ 为分形维数。

Theorem 3（吸引子维数估计）.
在参数满足 $\Gamma_2^{-1}, \Gamma_0 \gg |\alpha|, \beta$ 的物理相关区域内，整体吸引子 $\mathcal{A}$ 的 Lyapunov 维数 $d_L(\mathcal{A})$ 满足如下显式上界：

二维情形（ $$d = 2$$ ）：
$d_L(\mathcal{A}) \le c \max\left\{ \left( \frac{\Gamma_0}{\Gamma_2} - \alpha \Gamma_2^{-1} \right)^{1/2};\ \left| \alpha - \frac{1}{2}\frac{\Gamma_0}{\Gamma_2} \right|^{3/8} \Gamma_2^{-5/8} \beta^{-1/4} \right\},$
三维情形（ $$d = 3$$ ）：
$d_L(\mathcal{A}) \le c \max\left\{ \left( \frac{\Gamma_0}{\Gamma_2} - \alpha \Gamma_2^{-1} \right)^{3/4};\ 2^{9/4} \Gamma_2^{-1/2} \beta^{-1/2} \left| \alpha - \frac{1}{2}\frac{\Gamma_0}{\Gamma_2} \right|^{1/2} \right\},$
其中 $$c > 0$$ 为绝对常数。特别地，当 $\Gamma_0/\Gamma_2 \gg |\alpha|$ 时，主导项为 $(\Gamma_0/\Gamma_2)^{1/2}$ （2D）或 $(\Gamma_0/\Gamma_2)^{3/4}$ （3D），从而导出分辨率尺度 $\eta_{\text{res}} \sim \eta_{\text{sh}} = \sqrt{\Gamma_2/\Gamma_0}$ ，严格证实了 Swift-Hohenberg 尺度在主动湍流中的根本性作用。

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

a) 核心分析学工具

本文的证明体系建立在现代无穷维动力系统与非线性 PDE 理论的坚实基础上，综合运用了以下几类核心工具：

1. Galerkin 截断与紧性论证：
为构造弱解，作者采用经典的 Galerkin 方法，将方程投影到 Stokes 算子前 $$N$$ 个特征模上，得到有限维 ODE 系统（见公式 (2.1)）。其解的存在唯一性由 ODE 理论保证。随后，通过一系列先验能量估计（见 (2.5), (2.12), (2.14)），获得关于 $$u_N$$ 在 $L^\infty_t L^2_x$ 、 $$L^2_t H^2_x$$ 、 $$L^4_t L^4_x$$ 等范数的一致（关于 $$N$$ ）有界性。这些估计最终允许应用 Banach-Alaoglu 定理（弱*紧性）与 Aubin-Lions 引理（强紧性），从 Galerkin 序列 $\{u_N\}$ 中抽取子列，收敛到某个 $$u$$ ，并验证其满足 TTSH 方程的弱形式。

2. 能量方法与插值不等式：
所有先验估计均基于对不同阶数的 $$L^2$$ -内积（如 $\langle \cdot, u_N \rangle$ , $\langle \cdot, -\Delta u_N \rangle$ , $\langle \cdot, \Delta^2 u_N \rangle$ ）进行操作，从而导出关于 $\|u_N\|_{L^2}$ 、 $\|\nabla u_N\|_{L^2}$ 、 $\|\Delta u_N\|_{L^2}$ 的微分不等式。在此过程中，大量使用了 Sobolev 插值不等式，例如：

\|\nabla u\|_{L^2} \le \|u\|_{L^2}^{1/2} \|\Delta u\|_{L^2}^{1/2}, \quad \|u\|_{L^6} \lesssim \|u\|_{H^1}, \quad \|u\|_{L^\infty} \lesssim \|u\|_{H^2} \quad (\text{for } d=2,3),

以及 Young 不等式 $ab \le \varepsilon a^p + C_\varepsilon b^q$ 来吸收坏项（如 (2.2), (2.4), (2.17)）。

3. Lieb-Thirring 不等式与迹公式：
为估计整体吸引子的维数，作者采用了 Constantin-Foias 理论的经典路径：利用线性化算子的迹公式（trace formula）与Lieb-Thirring 不等式来控制其特征值衰减。具体而言，对 TTSH 线性化算子 $$L(t)$$ ，其第 $$m$$ 个特征值 $\lambda_m(t)$ 满足

\sum_{j=1}^m \lambda_j(t) \le \operatorname{Tr} L(t) = \int_{\mathbb{T}^d} \operatorname{div} F'(u(t,x)) \, dx,

其中 $F(u) = -(\alpha + \beta|u|^2)u + \Gamma_0 \Delta u + \Gamma_2 \Delta^2 u$ 。通过精细计算 $\operatorname{Tr} L(t)$ 并应用 Lieb-Thirring 型不等式（如 $\sum_j |\mu_j|^{p} \le C \|\nabla^k u\|_{L^q}^r$ ），最终导出维数上界（见 §3.1–3.3）。

4. 吸引子理论与 Lyapunov 指数谱：
整体吸引子 $\mathcal{A}$ 的存在性依赖于系统具有吸收集（absorbing ball）与渐近紧性（asymptotic compactness）。前者由能量估计 (2.6)–(2.8) 给出：存在 $$k_1$$ 使得对充分大的 $$t$$ ， $\|u(t)\|_{L^2} \le k_1$ ；后者则由高阶估计 (2.13)–(2.14) 保证。Lyapunov 维数 $d_L(\mathcal{A})$ 定义为最大整数 $$m$$ 加上比值 $$q$$ ，使得 $\sum_{j=1}^m \ell_j + q \ell_{m+1} = 0$ ，其中 $\ell_j$ 是 Lyapunov 指数谱（按降序排列）。数值模拟（§4）通过伪谱法计算该谱，验证了理论界的合理性。

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

1. 负扩散项的处理：
TTSH 方程中 $+\Gamma_0 \Delta u$ 项是反耗散（anti-diffusive），会破坏能量单调递减性。作者通过引入重定义参数 $\alpha_R = \alpha - \frac{1}{2}\Gamma_0^2/\Gamma_2$ （见 (2.3)），将反耗散效应部分吸收进有效线性驱动系数，从而在能量估计中构造出正定的耗散项（如 (2.5) 中的 $\frac{1}{2}\Gamma_2 \|\Delta u_N\|_{L^2}^2$ ），这是克服负扩散的关键代数技巧。

2. 高阶非线性与对流项的耦合控制：
三维情形下，对流项 $(u \cdot \nabla)u$ 与 $\Delta^2 u$ 的耦合导致边界项难以消失。作者利用Sobolev 嵌入与插值（如 $\|u\|_{L^3} \lesssim \|u\|_{L^2}^{1/2} \|\nabla u\|_{L^2}^{1/2}$ ）及Gagliardo-Nirenberg 不等式，将非线性项控制为 $\|\nabla u\|_{L^2}^a \|\Delta u\|_{L^2}^b$ 形式，再用 Young 不等式配平（见 (2.12), (2.17)），确保高阶耗散项 $\Gamma_2 \|\Delta^2 u\|_{L^2}^2$ 能主导所有坏项。

3. 无散度约束的隐式利用：
在二维情形下，作者明确利用了 $\nabla \cdot u = 0$ 导致的对流项正交性： $\int_{\mathbb{T}^2} (u \cdot \nabla u) \cdot \Delta u \, dx = 0$ （见 (2.11)），极大简化了梯度估计。此恒等式是二维 Navier-Stokes 理论的核心，在 TTSH 中同样发挥着“免费”控制作用。

4. 数值与解析的交叉验证：
为增强结论的可信度，作者不仅给出解析维数上界，还通过伪谱法（pseudospectral method）在 $$d=2$$ 进行直接数值模拟，计算 Lyapunov 谱，并确认其总和（即 Lyapunov 维数）落在理论区间内。这种“解析—数值”双轨验证，是当代复杂系统研究的典范。

与之前工作的比较

a) 突破与改进

本文相较于已有文献，实现了三方面的实质性突破：

第一，首次为 TTSH 方程建立严格的整体吸引子理论。
此前，Wensink 等人 [66] 提出了 TTSH 模型并进行了线性稳定性与数值研究；[67] 证明了其在全空间 $\mathbb{R}^3$ 上的全局适定性，但未涉及周期域上的吸收集构造与吸引子存在性。本文首次在物理上更相关的有界周期域 $\mathbb{T}^d$ 上，完整建立了从适定性到吸引子存在的逻辑链条，并给出了可计算的维数上界，填补了该模型数学理论的关键空白。

第二，显著放宽了对初始数据正则性的要求。
以往关于主动流体方程（如 ITT 方程 [3, 28]）的结果，往往需要 $$H^1$$ 或 $$H^2$$ 初值以保证强解存在。而本文 Theorem 1 仅假设 $u_0 \in L^2$ ，便获得了 $C([0,T]; L^2) \cap L^2((0,T); H^2)$ 的弱解，这得益于 $\Delta^2 u$ 项提供的强大正则化效应。该结果将 TTSH 的适定性门槛降至与经典 Navier-Stokes 相当的水平，极大地拓展了其适用范围。

第三，维数估计首次与物理尺度精确对接。
现有文献（如 KS 方程 [35] 或 Ginzburg-Landau 方程 [4, 23]）虽有吸引子维数估计，但其上界常含抽象常数或与物理参数关系模糊。本文的 Theorem 3 则首次导出显式、无量纲、与 Swift-Hohenberg 尺度 $\eta_{\text{sh}}$ 直接关联的维数公式。例如，2D 维数上界中的 $(\Gamma_0/\Gamma_2)^{1/2}$ 正比于 $\eta_{\text{sh}}^{-1}$ ，从而由 Landau-Lifshitz 关系 $N \sim (L/\eta_{\text{res}})^d$ 直接得出 $\eta_{\text{res}} \sim \eta_{\text{sh}}$ 。这不仅是数学估计的优化，更是对“为何实验中涡旋尺寸锁定在 $\eta_{\text{sh}}$ ”这一核心物理问题的第一个严格解答。

综上，本文通过融合现代 PDE 分析、动力系统理论与高性能数值模拟，不仅完善了 TTSH 方程的数学基础，更架起了严格数学与复杂物理现象之间一座坚实的桥梁。

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

1. 具有密度依赖黏性应力张量的可压缩非牛顿型Korteweg系统的耗散解

Dissipative Solutions to a Compressible Non-Newtonian Korteweg System with Density-Dependent Viscous Stress Tensor

作者: Didier Bresch, Christophe Lacave, Maja Szlenk

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19442v1

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

摘要翻译

中文标题: 具有密度依赖黏性应力张量的可压缩非牛顿型Korteweg系统的耗散解

摘要翻译: 本文的主要目标是证明：当考虑毛细效应时，在周期区域 $\T^d$ （其中 $$d=2,3$$ ）上，描述具有密度依赖黏性的可压缩黏塑性流动的系统存在耗散解。我们推导并计算了相对熵不等式，并由此证明了该系统耗散解的存在性以及弱解—强解唯一性（weak–strong uniqueness）。本结果将近期关于欧拉–Korteweg系统与纳维–斯托克斯–Korteweg系统之间联系的研究进展——该进展针对黏性依赖于密度的牛顿型流动（参见 D. Bresch, M. Gisclon, I. Lacroix-Violet, Arch. Rational Mech. Anal. (2019)）——推广至非牛顿型流动情形。

2. 无穷维空间中的随机持续性

Stochastic Persistence in Infinite Dimensions

作者: Juraj Foldes, Declan Stacy

类别: math.PR, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19145v1

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

摘要翻译

中文标题: 无穷维空间中的随机持续性

摘要翻译:
受无穷维生态学与生物学模型（如反应-扩散型随机偏微分方程（SPDEs）和随机泛函微分方程）的启发，我们建立了一套关于随机持续性（即种群共存）的一般性判据，该判据以平均Lyapunov函数为特征——此前此类结果仅在有限维情形下得到确立。为将本成果应用于SPDEs，我们系统分析了其射影过程，并综合运用了温和解方法（特别是随机卷积技术）与变分方法（特别是Lyapunov函数构造）。此外，我们的分析还依赖于反应-扩散型SPDEs若干非平凡的适定性与非负性结果；我们以极大一般性陈述并严格证明了这些结果，显著拓展了文献中已知的相应结论。最后，我们给出了若干应用实例，包括生态学模型（Lotka–Volterra系统）、传染病模型（SIR模型）以及湍流模型。特别值得注意的是，正如常微分随机系统（SDEs）情形所示，在Lotka–Volterra模型中，种群共存性同样由入侵率（invasion rates）所决定。

3. 重力作用下斯托克斯流中界面的不稳定动力学

Unstable Interface Dynamics for Gravity Stokes Flow

作者: Francisco Gancedo, Rafael Granero-Belinchón, Zhongtian Hu, Elena Salguero, Yao Yao

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18460v1

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

摘要翻译

中文标题: 重力作用下斯托克斯流中界面的不稳定动力学

摘要翻译: 本文研究了由两种密度不同的不可压流体所构成的界面在重力作用下遵循经典斯托克斯定律演化时所呈现的若干不稳定行为。在不稳定构型（即高密度流体位于低密度流体上方）下，我们证明了界面长度或曲率在无穷时间内持续增长。我们通过数值模拟验证了上述理论结果，模拟结果与所预测的增长现象高度吻合。而在稳定构型（即高密度流体位于低密度流体下方）下，我们进一步指出：某些初始界面构型将在有限时间内演化进入不稳定区域。

4. 具有二次反应动力学的耦合达西–福希海默–布林克曼系统的整体存在性与有限时间爆破

Global Existence and Finite-Time Blow-Up for a Coupled Darcy-Forchheimer-Brinkman System with Quadratic Reaction Dynamics

作者: Sahil Kundu, Manmohan Vashisth, Manoranjan Mishra

类别: math.AP

发布日期: 2026-01-25

arXiv ID: 2601.17984v1

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

摘要翻译

中文标题: 具有二次反应动力学的耦合达西–福希海默–布林克曼系统的整体存在性与有限时间爆破

摘要翻译:
本文研究一类非线性耦合系统，该系统将达西–福希海默–布林克曼（Darcy–Forchheimer–Brinkman）方程与一个对流–扩散–反应方程相耦合，其物理背景为多孔介质中的反应性物质输运。该模型包含非线性黏度耦合项、福希海默惯性阻力项、对流输运机制以及一个二次型反应项。我们首先针对任意初值证明了局部时间弱解的存在性。在物理上自然的初值条件 $0 \leq c_0 \leq 1$ 下，我们建立了浓度变量的最大值原理，从而在二维与三维空间中获得了弱解的整体存在性与唯一性。对于正则性更高的初值，我们进一步证明了强解的存在性、唯一性及其关于初值的连续依赖性。在此正则性框架下，浓度在所有 $$L^p$$ -范数（ $1 \leq p \leq \infty$ ）中均以一致的指数衰减速率趋近于零。相反，若初值满足 $$c_0 > 1$$ ，我们则证明了解的有限时间爆破现象，并给出了爆破时间的一个显式上界估计。最后，本文基于有限元方法开展数值模拟，直观地验证了理论所预言的指数衰减行为与有限时间爆破现象。

5. 具有指数型非线性的耦合达西–对流–扩散系统的整体适定性与数值逼近

Global Well-Posedness and Numerical Approximation of a Coupled Darcy-Convection-Diffusion System with Exponential Nonlinearity

作者: Sahil Kundu, Amiya K. Pani, Manoranjan Mishra

类别: math.AP, math-ph

发布日期: 2026-01-25

arXiv ID: 2601.17968v1

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

摘要翻译

中文标题: 具有指数型非线性的耦合达西–对流–扩散系统的整体适定性与数值逼近

摘要翻译:
本文研究多孔介质中由密度差驱动的流动，重点关注黏度对比度、密度对比度及线性吸附效应的作用机制。在此设定下，上层流体比下层流体更重且更具黏性；在重力作用下，该系统失稳并产生指状（finger-like）不稳定性结构。该现象在数学上通过将达西定律与一个对流–扩散–反应方程耦合来刻画。模型中的非线性主要源于黏度对浓度的依赖关系以及对流输运项。利用Galerkin逼近法与截断技术，我们在二维和三维情形下均证明了弱解唯一存在性。此外，借助极值原理可证得浓度变量的非负性。进一步，我们分析了解的长时间行为，并证明：当 $t \to \infty$ 时，浓度在 $$ L^p $$ 范数（对所有 $1 \le p \le \infty$ ）下以指数速率收敛于零。为补充理论分析，我们基于压力公式开展数值模拟；通过追踪总动能与混合度量随时间的演化，分别讨论了系统的不稳定性特征与混合效率。本研究揭示：尽管增大密度对比度会增强总动能，但其边际增幅随密度对比度的持续增加而逐渐减弱；类似地，吸附虽能抑制混合，但其抑制效能亦随吸附强度的进一步提升趋于饱和。这些非线性敏感性特征既被我们的理论估计所预示，亦得到数值模拟结果的明确验证。

不相关论文（基本信息）

1. 关于非局部薄膜方程的α-熵解：存在性与有限传播速度

On $$α$$ -entropy solutions of a nonlocal thin film equation: existence and finite speed of propagation

作者: Antonio Segatti, Roman Taranets

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19808v1

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

摘要: We consider an initial-boundary value problem for a class of nonlocal thin film equations governed by the spectral fractional Laplacian with homogeneous Neumann boundary conditions. We were the first to establish an $$α$$ -entropy estimate for nonlocal thin film equations, which yields essential a priori bounds for the regularity and long-time behavior of weak solutions. By developing a localized version of this estimate, we prove finite speed of propagation, showing that the support of nonnegative solutions remains compact for positive times. Furthermore, we find a sufficient condition for a waiting time phenomenon, whereby the solution remains identically zero in a region for a nontrivial time interval. These results highlight new features in the interaction between nonlocal effects and classical thin film dynamics.

2. 物理学中出现的椭圆方程的稳定解与有限莫尔斯指标解的先验估计

A priori estimates of stable and finite Morse index solutions to elliptic equations that arise in Physics

作者: J. Silverio Martínez-Baena

类别: math.AP, math-ph

发布日期: 2026-01-27

arXiv ID: 2601.19801v1

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

摘要: This thesis studies qualitative properties of solutions to nonlinear elliptic equations of Poisson type with Dirichlet boundary conditions that arise from some physical phenomena, with a particular focus on regularity, stability, and multiplicity of solutions. Building on the modern framework of solution stability and Morse index theory, the work investigates how these notions influence regularity in nonlinear elliptic problems. A central contribution is the construction of a counterexample showing that bounded radial Morse index does not prevent singular behavior of solutions in dimensions three through nine, challenging a natural extension of the Brezis-Vázquez regularity conjecture. In addition, optimal regularity results are established for radial solutions of a non-autonomous Hardy-Hénon equation, identifying the precise range of dimensions for which regularity holds. The thesis also addresses existence and multiplicity results for elliptic equations involving nonlinearities with spatially vanishing coefficients. Under suitable assumptions, the existence of multiple distinct solutions is proved using variational and topological methods. Finally, the thesis outlines several directions for future research, including extensions of stability-based regularity techniques to non-autonomous problems and potential applications of these techniques to field theories arising in theoretical physics.

3. 基于非局部梯度的一般非局部弹性理论及其与Eringen模型的联系

A general theory of nonlocal elasticity based on nonlocal gradients and connections with Eringen’s model

作者: J. C. Bellido, G. García-Sáez

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19797v1

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

摘要: We develop a general theory of nonlocal linear elasticity based on nonlocal gradients with general radial kernels. Starting from a nonlocal hyperelastic energy functional, we perform a formal linearization around the identity deformation to obtain a system of nonlocal linear elasticity equations. We establish the existence and uniqueness of weak solutions for both Dirichlet and Neumann boundary conditions, proving a general Korn-type inequality for nonlocal gradients. We show that this framework encompasses Eringen’s nonlocal elasticity model as a particular case, establishing an explicit connection between the two formulations. Finally, we prove localization results demonstrating that solutions to the nonlocal problems converge to their classical local counterparts in two different regimes: as the interaction horizon vanishes and, in the fractional case, as the fractional parameter approaches one. These results provide a comprehensive and unified mathematical foundation for nonlocal elasticity theories.

4. 分数阶方程定性性质研究方法

Methods in studying qualitative properties of fractional equations

作者: Wenxiong Chen, Yahong Guo, Congming Li

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19783v1

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

摘要: In this paper, we systematically review a series of effective methods for studying the qualitative properties of solutions to fractional equations. Beginning with the pioneering extension method and the method of moving planes in integral forms, we introduce a variety of direct methods, including the direct method of moving planes, the method of moving spheres, blow-up and rescaling techniques, the sliding method, regularity lifting, and approaches for interior and boundary regularity estimates. To elucidate the core ideas behind these methods, we employ simple examples that demonstrate how they can be applied to investigate qualitative properties of solutions. We also provide a comparative discussion of their respective strengths and limitations. It is our hope that this paper will serve as a useful handbook for researchers engaged in the study of fractional equations.

5. 非柱状区域上分数阶半线性波动方程的研究

On fractional semilinear wave equations in non-cylindrical domains

作者: Mauro Bonafini, Van Phu Cuong Le, Riccardo Molinarolo

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18355v2

链接: http://arxiv.org/abs/2601.18355v2

摘要: In this paper, we investigate a class of semilinear wave equations in non-cylindrical time-dependent domains, subject to exterior homogeneous Dirichlet conditions. Under mild regularity and monotonicity assumptions on the evolving spatial domains, we establish existence of weak solutions by two different methods: a constructive time-discretization scheme and a penalty approach. The analysis applies to nonlocal fractional Laplacians and potentials with Lipschitz continuous gradient, and to vector-valued maps.

6. 具有吸引性汤川势的赝相对论费米系统

Pseudo-relativistic fermionic systems with attractive Yukawa potential

作者: Bin Chen, Yujin Guo, Phan Thành Nam, Dong Hao Ou Yang

类别: math-ph, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19751v1

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

摘要: We study the Hartree-Fock and Hartree-Fock-Bogoliubov theories for a large fermionic system with the pseudo-relativistic kinetic energy and an attractive Yukawa interaction potential. We prove that the system is stable if and only if the total mass does not excess a critical value, and investigate the existence and properties of ground states in both sub-critical and critical mass regimes.

7. 基于格林函数的封闭框架：用于泊松方程及其广义上下解

A Green’s Function-Based Enclosure Framework for Poisson’s Equation and Generalized Sub- and Super-Solutions

作者: Kazuaki Tanaka, Ryoga Iwanami, Kaname Matsue, Hiroyuki Ochiai

类别: math.NA, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19682v1

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

摘要: This paper presents a novel framework for enclosing solutions of Poisson’s equation based on generalized sub- and super-solutions constructed using fundamental solutions. The conventional definition of sub- and super-solutions based on variational inequalities often fails for natural function classes such as piecewise linear functions and encounters theoretical difficulties in non-convex polygonal domains, where H^2 regularity is lost because of corner singularities. To overcome these limitations, we introduce the concept of ``Green-representable solutions’’ utilizing test functions constructed from fundamental solutions. This framework enables a new formulation of sub- and super-solutions that permits rigorous pointwise evaluation. For one-dimensional problems, we derive explicit constructions of the test functions. For two-dimensional polygonal domains, we employ the Method of Fundamental Solutions to generate test functions. The approach is validated through numerical experiments in both settings, including non-convex polygons. The results demonstrate that the proposed method yields strict and accurate pointwise enclosures of the true solution, even for problems with discontinuous source terms or geometric singularities.

8. 伪富克斯群的定量傅里叶不确定性原理与谱隙

Quantitative FUP and spectral gap for quasi-Fuchsian group

作者: Long Jin, An Zhang, Hong Zhang

类别: math.CA, math.AP, math.GT, math.SP

发布日期: 2026-01-27

arXiv ID: 2601.19663v1

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

摘要: We derive an explicit formula for the exponent $$β$$ in the higher-dimensional fractal uncertainty principle (FUP) established by Cohen 2023, quantifying its dependence on the porosity parameter $$ν$$ of the Fourier support. This quantitative version of FUP yields an explicit essential spectral gap for convex co-compact hyperbolic 3-manifolds arising from quasi-Fuchsian groups, thereby refining the result of Tao 2025. Our result extends the earlier work of Jin-Zhang 2020 to higher dimensions.

9. 非迹类噪声的精确界及其在随机偏微分方程中的应用

Sharp bounds for non-trace class noise and applications to SPDEs

作者: Antonio Agresti, Fabian Germ, Mark Veraar

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

发布日期: 2026-01-27

arXiv ID: 2601.19639v1

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

摘要: In the study of stochastic PDEs with colored, non-trace class space-time noise, one frequently encounters Gaussian series of the form

\sum_{n\geq 1} γ_n μ_n f_n,

where $(γ_n)_{n}$ is a sequence of standard independent Gaussian variables, $$g$$ is an $L^η(\mathcal{O})$ function, $(μ_n)_{n}$ is a sequence of scalars, and $$(f_n)_n$$ is an orthonormal system in $L^2(\mathcal{O})$ where $\mathcal{O} \subseteq \mathbb{R}^d$ is an open set. In this manuscript, we establish necessary and sufficient conditions for the above sum to converge in Bessel potential spaces $H^{-s,q}(\mathcal{O})$ . The latter can be interpreted as a Sobolev embedding for Gaussian series. Our main theorem is formulated using weighted sequence spaces that encode the $L^\infty$ -growth of the orthonormal system $(f_n)_{n}$ , a feature that is crucial for obtaining sharp estimates. We apply our results to the stochastic heat equation with additive non-trace class noise. In this case, our conditions capture the scaling relationship between the heat operator and the coloring of the noise.

10. 大N极限下二维O(N)线性σ模型的质量生成

Mass generation for the two dimensional O(N) Linear Sigma Model in the large N limit

作者: Matías G. Delgadino, Scott A. Smith

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

发布日期: 2026-01-27

arXiv ID: 2601.19630v1

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

摘要: This work studies the $$O(N)$$ Linear Sigma Model on $\mathbb{R}^{2}$ under a scaling dictated by the formal $$1/N$$ expansion. We show that in the large $$N$$ limit, correlations decay exponentially fast, where the acquired mass decays exponentially in the inverse temperature. In fact, each marginal converges to a massive Gaussian Free Field (GFF) on $\mathbb{R}^{2}$ , quantified in the $$2$$ -Wasserstein distance with a weighted $H^{1}(\mathbb{R}^{2})$ cost function. In contrast to prior work on the torus via parabolic stochastic quantization, our results hold without restrictions on the coupling constants, allowing us to also obtain a massive GFF in a suitable double scaling limit. Our proof combines the Feyel/Üstünel extension of Talagrand’s inequality with some classical tools in Euclidean Quantum Field Theory.

11. 自由间断泛函在有界变形函数空间中具有线性增长时的Γ-收敛性与均匀化

$$Γ$$ -convergence and homogenisation for free discontinuity functionals with linear growth in the space of functions with bounded deformation

作者: Gianni Dal Maso, Davide Donati

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19591v1

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

摘要: We study the $$Γ$$ -convergence of sequences of free discontinuity functionals with linear growth defined in the space ${\rm BD}$ of functions with bounded deformation. We prove a compactness result with respect to $$Γ$$ -convergence and outline the main properties of the $$Γ$$ -limits, which lead to an integral representation result. The corresponding integrands are obtained by taking limits of suitable minimisation problems on small cubes. These results are then used to study the deterministic and stochastic homogenisation problem for a large class of free discontinuity functionals defined in ${\rm BD}$ .

12. 一类双线性波动方程的近似可控性及最短控制时间

Approximate controllability of a bilinear wave equation and minimum time

作者: Karine Beauchard, Thomas Perrin, Eugenio Pozzoli

类别: math.OC, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19544v1

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

摘要: We study the global approximate controllability (GAC) of a Klein-Gordon wave equation, posed on the torus $\mathbb{T}^d$ of arbitrary dimension $d\in \mathbb{N}^*$ , with bilinear control potentials supported on the first $$(2d+1)$$ -Fourier modes. Let $Z(W_0)\subset \mathbb{T}^d$ be the set of essential zeroes of the initial state $W_0\in H^1\times L^2(\mathbb{T}^d)$ , and $r(W_0)\geq 0$ be the maximum radius of a ball of $\mathbb{T}^d$ contained in $$Z(W_0)$$ . Due to finite speed of propagation, the minimum control time starting from $$W_0$$ is necessarily larger than or equal to $$r(W_0)$$ . We prove the following three facts. In low dimensions $d \in \{1,2\}$ : the minimum time for GAC from $W_0 \neq 0$ is equal to $$r(W_0)$$ . In any dimensions $d\geq 3$ : the minimum time for GAC from $$W_0$$ is zero if $$Z(W_0)$$ has zero Lebesgue measure; and the GAC in sufficiently large time from all $W_0\neq 0$ . The proof strategy consists in combining Lie bracket techniques \emph{à la Agrachev-Sarychev} with the propagation of well-prepared positive states.

13. 无对称性条件下的波动映射爆破稳定性

Blowup stability of wave maps without symmetry

作者: Roland Donninger, Frederick Moscatelli

类别: math.AP, math-ph

发布日期: 2026-01-27

arXiv ID: 2601.19516v1

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

摘要: We study wave maps from $$(1+d)$$ -dimensional Minkowski space into the $$d$$ -sphere without any symmetry assumptions. There exists an explicit self-similar blowup solution and we prove that this solution is asymptotically stable under small perturbations of the initial data. The proof is fully rigorous and requires no numerical input whatsoever.

14. 高维无对称性自相似波映射的模稳定性

Mode stability of self-similar wave maps without symmetry in higher dimensions

作者: Roland Donninger, Frederick Moscatelli

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

发布日期: 2026-01-27

arXiv ID: 2601.19515v1

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

摘要: We consider wave maps from $$(1+d)$$ -dimensional Minkowski space into the $$d$$ -sphere. For every $d \geq 3$ , there exists an explicit self-similar solution that exhibits finite time blowup. This solution is corotational and its mode stability in the class of corotational functions is known. Recently, Weissenbacher, Koch, and the first author proved mode stability without symmetry assumptions in $$d =3$$ . In this paper we extend this result to all $d \geq 4$ . On a technical level, this is the first successful implementation of the quasi-solution method where two additional parameters are present.

15. 关于若干非局部非线性扩散问题

On some nonlocal, nonlinear diffusion problems

作者: M. M. Chipot, A. Luthra, S. A. Sauter

类别: math.AP, math.NA

发布日期: 2026-01-27

arXiv ID: 2601.19478v1

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

摘要: This note is devoted to some nonlocal, nonlinear elliptic problems with an emphasis on the computation of the solution of such problems, reducing it in particular to a fixed point argument in R. Errors estimates and numerical experiments are provided.

16. 一般紧致障碍物散射问题的积分方程方法：波数显式估计

Integral equation methods for scattering by general compact obstacles: wavenumber-explicit estimates

作者: Simon N. Chandler-Wilde, Siavash Sadeghi

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19456v1

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

摘要: There has been significant recent interest in understanding the dependence on the wavenumber, $$k$$ , of boundary integral operators (BIOs), supported on some set $Γ\subset \mathbb{R}^n$ , that arise in the solution of BVPs for the Helmholtz equation, $$Δu + k^2 u=0$$ . Recently, for the Dirichlet BVP with data $$g$$ , Caetano et al (2025) have proposed an integral equation (IE) $$A_kφ=g$$ that applies for arbitrary compact $$Γ$$ . This formulation is a generalisation of standard first kind IEs, where the BIO is $$S_k$$ , the single-layer BIO on a surface $$Γ$$ , that apply when $$Γ$$ is the boundary of a Lipschitz domain or a screen. In this paper we study the dependence of $$A_k$$ on $$k$$ , showing that, for $k\geq k_0>0$ , $\|A_k\|\leq ck$ while $\|A_k^{-1}\| \leq c'k$ if $$Γ$$ is star-shaped, where $$c, c'>0$$ depend only on $$k_0$$ and $$Γ$$ . Amongst other bounds we also show that: (i) on the one hand, given any mildly increasing unbounded positive sequence $$(k_m)$$ and any unbounded sequence $$(a_m)$$ , there exists $$Γ$$ , with connected complement, such that $\|A_{k_m}^{-1}\|\geq a_m$ for every $$m$$ ; (ii) on the other hand, for every $Γ\subset \mathbb{R}^n$ and $k_0,\varepsilon, δ>0$ , there exists $$c>0$$ and $E\subset [k_0,\infty)$ , with Lebesgue measure $m(E)\leq \varepsilon$ , such that $\|A_{k}^{-1}\|\leq c k^{2n+2+δ}$ on $[k_0,\infty)\setminus E$ , i.e., the growth of $\|A_{k}^{-1}\|$ is at worst polynomial in $$k$$ if one avoids a set $$E$$ of arbitrarily small measure. As a corollary of these results we obtain the first $$k$$ -explicit bounds on $\|S_k^{-1}\|$ and the condition number of $$S_k$$ for the case that $$Γ$$ is the boundary of a Lipschitz domain, or a screen not contained in a hyperplane, and analogous estimates for the case that $$Γ$$ is a $$d$$ -set (and so of Hausdorff dimension $$d$$ ), for non-integer values of $$d$$ .

17. 多发性硬化症约束型聚集-扩散-反应模型弱解的存在性

Existence of Weak Solutions to a Constrained Aggregation-Diffusion-Reaction Model for Multiple Sclerosis

作者: S. Fagioli, M. Kamath Katapady

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19427v1

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

摘要: We establish an existence result for weak solutions to an aggregation-diffusion-reaction equation with a constraint, arising in the modelling of multiple sclerosis. The model is derived from a general chemotaxis-type framework and describes the time evolution of the density of activated macrophages, which is subject to attraction by oligodendrocytes. The latter are governed by a constraint equation. The proof relies on a variational splitting scheme that isolates the transport (aggregation-diffusion) and reaction contributions. The structure of the constraint makes it possible to recover the oligodendrocyte density as the limit of a sequence of characteristic functions.

18. 一个精确的单项式Caffarelli–Kohn–Nirenberg不等式

A sharp monomial Caffarelli-Kohn-Nirenberg inequality

作者: Francesco Pagliarin

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19359v1

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

摘要: We consider a monomial Caffarelli-Kohn-Nirenberg inequality, find the optimal constant and classify the optimizers under an integrated curvature dimension condition. We take advantage of the $$Γ$$ -calculus to exploit geometrical techniques to tackle the problem and regularity results to justify some integration by parts. A symmetry-breaking result is also provided.

19. 关于带乘性噪声的非线性薛定谔方程的爆破问题

On blow up NLS with a multiplicative noise

作者: Chenjie Fan, Junzhe Wang

类别: math.AP, math.PR

发布日期: 2026-01-27

arXiv ID: 2601.19330v1

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

摘要: It is of significant interest to understand whether a noise will speed up or prevent blow up. Under certain nondegenerate conditions, \cite{dD2005Blowup} proved a multiplicative noise will speed up blow up of NLS, in the sense that, blow up can happen in any short time with positive probability. We prove that such probability is indeed quite small, and provide a large deviation type upper bound.

20. 关于通过无散度变换研究线性椭圆方程适定性的若干注记

Remarks on well-posedness for linear elliptic equations via divergence-free transformation

作者: Haesung Lee

类别: math.AP, math.FA

发布日期: 2026-01-27

arXiv ID: 2601.19317v1

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

摘要: This paper investigates the well-posedness of linear elliptic equations, focusing on the divergence-free transformation introduced in the author’s recent work [J. Math. Anal. Appl. 548 (2025), 129425]. By comparing this approach with classical bilinear form methods, we demonstrate that while standard techniques encounter limitations in handling zero-order coefficients $c \in L^1(U)$ , the divergence-free transformation successfully establishes well-posedness in this setting. Furthermore, utilizing the Riesz-Thorin interpolation theorem between the cases $c \in L^1(U)$ and $c \in L^{\frac{2d}{d+2}}(U)$ , we establish the existence and uniqueness of weak solutions under the assumption $c \in L^s(U)$ for $s \in [1, \frac{2d}{d+2}]$ .

21. 用于积分-发放模型的、在无穷远处具有超线性漂移项的福克-普朗克方程

A Fokker-Planck equation with superlinear drift at infinity for Integrate-and-Fire model

作者: Benoît Perthame, Clément Rieutord, Delphine Salort

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19282v1

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

摘要: The Integrate-and-Fire model is a Fokker-Planck equation arising in neuroscience. It describes the evolution of the probability density of the neuronal membrane potential and fitting has shown that the inclusion of a em superlinear drift provides the most realistic description. To make sense of this, we propose to set the equation on the full line, the neural activity being described by the flux at infinity. This framework serves as a model extension of the classical Noisy Integrate-and-Fire model, with a fixed firing potential. We first establish the well-posedness of the solution, establish the boundary condition at infinity which is the major difficulty. Then, state rigorously the entropy dissipation property. Finally, using Doeblin’s method, we prove the exponential convergence of the solution toward the unique stationary state in full generality.

22. 平面上的变系数椭圆结构：输运动力学、刚性理论与函数论

Variable Elliptic Structures on the Plane: Transport Dynamics, Rigidity, and Function Theory

作者: Daniel Alayón-Solarz

类别: math.CV, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19274v1

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

摘要: We study variable elliptic structures in the plane defined by a smoothly varying quadratic relation i^2 + beta(x,y) i + alpha(x,y) = 0, and the associated first order operator dbar = 1/2 (dx + i dy). Differentiating the structure relation yields explicit expressions for the derivatives of i(x,y) in terms of the coefficient functions alpha and beta, leading to a universal transport system governing their admissible variations. In the elliptic regime this system reduces to a forced complex Burgers equation for a scalar spectral parameter encoding the structure coefficients. We identify a rigidity condition under which the transport becomes conservative, and show that in this regime the generalized Cauchy Riemann operator satisfies a Leibniz rule and admits a factorization of the associated second order operator into first order components. As a consequence, classical tools of planar complex analysis, including Cauchy Pompeiu type formulas, integral representations, and elliptic second order operators, reappear in a variable coefficient setting with explicit structure. The theory is developed at the level of direct computation, emphasizing transparency of the integrability mechanism and the interplay between transport dynamics, rigidity, and function theory.

23. 与晶界运动相关的Łojasiewicz–Simon不等式及其应用

The Łojasiewicz-Simon inequality related to grain boundary motion and its applications

作者: Masashi Mizuno, Ayumi Sakiyama, Keisuke Takasao

类别: math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19226v1

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

摘要: In this paper, we study the Łojasiewicz-Simon gradient inequality for the mathematical model of grain boundary motion. We first derive a curve shortening equation with time-dependent mobility, which guarantees the energy dissipation law for the grain boundary energy, including the difference between orientations of the constituent grains as a state variable. Next, we discuss the Łojasiewicz-Simon gradient inequality for the grain boundary energy. Finally, we give applications of the inequality to the energy.

24. 关于一类刻画布朗运动的多维非线性时间分数阶福克-普朗克方程

On a Class of Multi-Dimensional Non-linear Time-Fractional Fokker-Planck Equations Capturing Brownian Motion

作者: Neetu Garg, Varsha R

类别: math.NA, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19211v1

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

摘要: The time-fractional Fokker-Planck equation is a key model for characterizing anomalous diffusion, stochastic transport, and non-equilibrium statistical mechanics with applications in finance, chaotic dynamics, optical physics, and biological systems. In this work, we develop a semi-analytical solution for the multi-dimensional time-fractional Fokker-Planck equation employing the Laplace residual power series method. This method blends the Laplace transform and the traditional residual power series method, guaranteeing efficient solutions incorporating the memory and nonlocal effects. To validate the accuracy and effectiveness of the approach, we address several examples, including non-linear problems in multi-dimensions, and analyze the evolution of errors. The numerical simulations are compared with existing methods to confirm the adopted method’s strength. The smooth and stable error evolution promises that the suggested method is a powerful tool for analyzing time-fractional Fokker-Planck equations.

25. 条带型算子与抽象柯西问题

Strip-type operators and abstract Cauchy problems

作者: Nikolaos Roidos

类别: math.FA, math.AP

发布日期: 2026-01-27

arXiv ID: 2601.19075v1

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

摘要: We consider the non-homogeneous abstract linear Schrödinger and wave equations with zero initial conditions, defined by operators of strip-type and parabola-type in Banach spaces, respectively, and establish the well-posedness of classical solutions in appropriate vector-valued Sobolev-Slobodetskii spaces. We obtain analogous results for two extensions of these equations by replacing the previously mentioned boundedness properties of the associated operators with $$R$$ -boundedness. As an application, we consider an abstract semilinear wave equation and establish the existence and uniqueness of classical solutions to this problem for short times.

26. 基于截断Lamb核的弹性半空间上有限圆形薄板的动力响应分析

Dynamic Response of a Finite Circular Plate on an Elastic Half-Space Using the Truncated Lamb Kernel

作者: Greyson Meares, Sage Meiling, Charis Tsikkou

类别: math.AP, math.SP

发布日期: 2026-01-26

arXiv ID: 2601.19031v1

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

摘要: We develop an exact operator formulation for the dynamic interaction between a finite circular elastic plate and an elastic half-space. Classical analyses, beginning with Lamb’s representation of the half-space response, typically assume an infinite plate and rely on diagonalization of the soil operator via the continuous Hankel transform. For a plate of finite radius $$R$$ , however, both traction and displacement are supported only on $0 \le r \le R$ , leading to the spatially truncated Lamb operator

\mathscr{M}(ω) = χ_{[0,R]} \, T(ω)\, χ_{[0,R]},

where $$T(ω)$$ is the Hankel multiplier involving the Rayleigh denominator $$Ω(ξ,ω)$$ . Truncation destroys the diagonal structure of $$T(ω)$$ and introduces real-axis singularities associated with the Rayleigh pole, in addition to square-root branch points at $$ξ= k_T$$ and $$ξ= k_L$$ . We represent the action of $\mathscr{M}(ω)$ on a finite-disk Bessel basis $\{ φ_n(r) = A_{1,n} J_0(λ_n r) + A_{2,n} I_0(λ_n r)\},$ which satisfies the free-edge boundary conditions of the plate, and derive explicit expressions for the resulting matrix elements. These involve integrals of the Lamb kernel evaluated as Cauchy principal values, with residue contributions corresponding to radiation damping in the half-space. The resulting operator matrix is dense but spectrally convergent. Its inversion yields a complete frequency-domain solution for finite-radius plates. The analysis reproduces Chen et al.’s finite-radius experiments for small $$R$$ , approaches the infinite-radius limit as $R \to \infty$ , and quantifies finite-radius corrections. To our knowledge, this is the first exact operator-level treatment of finite-radius plate-half-space interaction that retains the full nonlocal Lamb kernel.

27. Kähler–Einstein 流形乘积上的 ℂℙ¹-丛上的 Ricci 流

Ricci Flow on CP1-bundles over a Product of Kähler-Einstein Manifolds

作者: Frederick Tsz-Ho Fong, Hung Tran

类别: math.DG, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18931v1

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

摘要: In this paper, we study the Ricci flow on CP1-bundles over a product of Kähler-Einstein manifolds whose initial metric is constructed by the ansatz used in works by M. Wang et. al. We prove that the ansatz is preserved along the Ricci flow. Furthermore, in the Kähler case, we proved that Type I finite-time singularity must occur under such an ansatz.

28. 面积最小的毛细管锥面

Area-minimizing capillary cones

作者: Benjy Firester, Raphael Tsiamis, Yipeng Wang

类别: math.DG, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18794v1

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

摘要: We construct non-flat minimal capillary cones with bi-orthogonal symmetry groups for any dimension and contact angle. These cones interpolate between rescalings of a singular solution to the one-phase problem and the free-boundary cone obtained by halving a Lawson cone along a hyperplane of symmetry. The existence and uniqueness of such cones is proved by solving a nonlinear free boundary equation parametrized by the contact angle and obtaining monotonicity properties for the solutions. The constructed cones are minimizing in ambient dimension $$8$$ or higher, for appropriate contact angles, demonstrating that the regularity theory for minimizing capillary hypersurfaces can have singularities in codimension $$7$$ and completing the capillary regularity theory for contact angles near $$π/2$$ . We further develop the connection between capillary hypersurfaces and solutions of the one-phase problem, consequently producing new examples of singular minimizing free boundaries for the Alt-Caffarelli functional.

29. 受趋化性启发的空气传播传染病传输偏微分方程模型：基于流行病学动机的数学与数值分析

Chemotaxis-inspired PDE models of airborne infectious disease transmission: epidemiologically-motivated mathematical and numerical analyses

作者: Alex Viguerie, Malú Grave, Alvaro L. G. A. Coutinho, Alessandro Veneziani, Thomas J. R. Hughes

类别: q-bio.PE, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18703v1

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

摘要: Partial differential equation (PDE) models for infectious diseases, while less common than their ordinary differential equation (ODE) counterparts, have found successful applications for many years. Such models are typically of reaction-diffusion type, and model spatial propagation as a diffusive process. However, given the complex nature of human mobility, such models are limited in their ability to describe airborne infectious diseases in human populations. Recent work has advocated for the inclusion of an additional chemotaxis-type term as an alternative; spatial propagation of infection fronts is assumed additionally to flow from low-to-high concentrations of susceptible populations. The present work extends the study of such models by providing an epidemiologically interpretable analysis, directly connecting model behavior to information readily available to policymakers. In particular, we derive a spatially-aware basic reproduction number, which accounts for spatial heterogeneity in population density. Furthermore, we discuss several important aspects concerning the numerical solution of the model, including the introduction of a stabilization scheme. Finally, we perform a series of simulation studies in the Italian region of Lombardy (severely affected by the COVID-19 outbreak in 2020) and in the US state of Georgia, in which we demonstrate the model’s potential to better capture important spatiotemporal dynamics observed in real-world data compared to pure reaction-diffusion models.

30. 无穷图上的半线性扩散方程：耗散情形与Lipschitz情形

Semilinear Diffusion Equations on Infinite Graphs: The Dissipative and Lipschitz Cases

作者: Elvise Berchio, Davide Bianchi, Alberto G. Setti, Maria Vallarino

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18549v1

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

摘要: We study a class of semilinear diffusion equations on infinite, connected, weighted graphs, focusing on two types of nonlinearities: monotone decreasing and Lipschitz continuous. Under minimal structural assumptions on the graph, we establish existence, uniqueness, and regularity of mild solutions for initial data in $\ell^p$ spaces, with $1\leq p<\infty$ . Our approach relies on time discretization via an implicit Euler scheme and an exhaustion technique using Dirichlet subgraphs. As a by-product, we obtain existence and uniqueness results for a related time-independent equation. Finite-time extinction and positivity for solutions under a specific forcing term are also proved.

31. 一个“刷子问题”：涉及细长区域与图上偏微分方程的均匀化理论

A brush problem. Homogenization involving thin domains and PDEs in graphs

作者: José M. Arrieta, Joaquín Domínguez-de-Tena

类别: math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18430v1

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

摘要: This work analyses the homogenization of a linear elliptic equation with Neumann boundary conditions in a comb/brush domain, composed of a fixed base and a family of thin teeth. The teeth are defined as rescalings of order less than or equal to $\varepsilon$ of a model tooth of arbitrary shape. Periodicity in their distribution is not assumed; instead, the existence of an asymptotic limit density $$θ$$ , which may vanish in certain regions, is assumed. The convergence analysis is performed using an adaptation of the unfolding operator method to a non-periodic framework. Finally, it is shown that, under certain conditions on the geometry of the teeth, the resulting limit problem can be interpreted as a differential equation on a graph.

32. 随机Ginzburg-Landau近似的高维空间有效性——基于Wiener代数的方法

Validity of the stochastic Ginzburg-Landau approximation in higher space dimensions -A Wiener algebra approach-

作者: Anna Logioti, Guido Schneider

类别: math.PR, math.AP

发布日期: 2026-01-26

arXiv ID: 2601.18406v1

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

摘要: We consider an anisotropic $$d$$ -dimensional Swift-Hohenberg model $\mathcal{O}(\varepsilon^2)$ -close to the first instability, where $0 < \varepsilon \ll 1$ is a small perturbation parameter. This model for pattern formation is perturbed with additive noise in time and space. By a multiple scaling ansatz we derive a stochastic $$ d $$ -dimensional Ginzburg-Landau equation for the approximate description of the bifurcating solutions. We prove the validity of the approximation by this amplitude equation on its natural time scale in case of $$ d $$ -dimensional periodic domains of length $\mathcal{O}(1/\varepsilon)$ for the Swift-Hohenberg model under suitable conditions on the additive noise. In detail, we prove the validity of this approximation for noise whose set of Fourier coefficients with respect to $$ x $$ is in $\ell^1$ for fixed $t \geq 0$ . Moreover, we improve existing approximation results in the sense that the stable part of the noise can be larger.

33. 自由边界Willmore问题的稳定性

Stability of the free boundary Willmore problem

作者: Anna Dall’Acqua, Fabian Rupp, Reiner Schätzle, Manuel Schlierf

类别: math.AP, math.DG

发布日期: 2026-01-26

arXiv ID: 2601.18388v1

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

摘要: We study the Willmore problem with free boundary by means of a new Łojasiewicz-Simon gradient inequality for functionals on infinite dimensional manifolds. In contrast to previous works, we do not rely on a gradient-like representation of the Fréchet derivative, but merely on an inequality. For the free boundary Willmore flow, we prove that solutions starting sufficiently close to a local minimizer exist for all times and converge. In the static setting, we prove quantitative stability of free boundary Willmore immersions and a local rigidity result in a neighborhood of free boundary minimal surfaces.

34. 调和映射到CAT(0)-空间的热流

Heat flow of harmonic maps into CAT( $$0$$ )-spaces

作者: Fang-Hua Lin, Antonio Segatti, Yannick Sire, Changyou Wang

类别: math.AP, math.DG

发布日期: 2026-01-26

arXiv ID: 2601.18046v1

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

摘要: We introduce a new approach to prove the global existence and uniqueness of suitable weak solutions of the heat flow of harmonic mappings into CAT(0) metric spaces. Our method allows also to prove Lipschitz continuity in spatial variables for such solutions into a wide class of CAT $$(0)$$ spaces, answering a long-standing open problem in the field. Our approach is based on an elliptic regularization of the gradient flow of the Dirichlet energy and even in the case of smooth Riemannian targets provides a novel viewpoint, together with a new Dynamical Variational Principle and a new proof of the celebrated Eells-Sampson theorem. The spatial Lipschitz regularity for such weak solutions is achieved by fully exploiting the variational structure of the problem at the regularized level and introducing a parabolic frequency function of Almgren-Poon type. Our contribution is the first instance of the use of monotonicity methods for parabolic deformations of maps into singular targets.

35. 非线性非局部算子及其相关椭圆型与抛物型问题定性性质的统一视角

A unified view of nonlinear, nonlocal operators and qualitative properties of associated elliptic and parabolic problems

作者: Ralph Chill, Mahamadi Warma

类别: math.AP

发布日期: 2026-01-25

arXiv ID: 2601.18028v1

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

摘要: We put together a general framework to deal with elliptic and parabolic equations associated with (nonlinear) nonlocal (fractional order) operators. Many well-known nonlocal operators enter into our framework, and in addition one may introduce many other, new nonlocal operators that have not yet been considered in the literature. We use the abstract theory of (nonlinear) semigroups generated by subgradients of proper, lower semicontinuous and convex functionals on Hilbert spaces to build a rigorous and applicable framework that works for many classical elliptic operators but also nonlocal or sometimes fractional operators. After recalling the notion of a nonlinear semigroup generated by subgradients and $$j$$ -subgradients of the associated energy functions, we introduce a general class of (nonlinear) nonlocal elliptic type operators and define rigorously subgradients and $$j$$ -subgradients of such functionals that generate (nonlinear) submarkovian semigroups and hence, the abstract Cauchy problem associated with these subgradients and/or $$j$$ -subgradients is wellposed. The existence and the qualitative properties of solutions to these Cauchy problems and the corresponding semigroups are investigated. More precisely, we show some comparison and maximum principles, submarkovian, domination, ultracontractivity properties, and some Hölder type estimates for these semigroups of operators. These results are usually useful in several branches of pure and applied partial differential equations. We finish the paper by giving several examples of nonlocal operators in Euclidean spaces, graphs, metric random walk spaces, fractional Brownian motions, and Lévy flights, that fit in our general framework.

36. 基于预解式方法的阿坦纳纳–巴莱安努演化方程研究：拉普拉斯符号、适定解与正则性

Resolvent Approach to Atangana–Baleanu Evolution Equations: Laplace Symbols, Mild Solutions, and Regularity

作者: Mohamed Wakrim

类别: math.AP

发布日期: 2026-01-25

arXiv ID: 2601.17992v1

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

摘要: Fractional evolution equations with memory terms are widely used to model anomalous diffusion, viscoelastic response, and hereditary dynamics in physics, biology, and engineering. Among the recently introduced operators, the Atangana–Baleanu (AB) derivatives have attracted considerable attention due to their non-singular Mittag–Leffler kernels. However, their analytic treatment remains limited, as the AB kernel does not fall within the classical Volterra or Bernstein-function frameworks. This paper develops a unified resolvent approach for AB-type evolution equations in Banach spaces. Using a Laplace-domain formulation inspired by Hille–Phillips theory, we introduce a fractional resolvent associated with the AB kernel and establish optimal bounds on sectorial contours. Under the natural condition $$β<1+α$$ , we construct an AB–Mittag–Leffler resolvent family and obtain a complete representation of mild solutions to the AB Cauchy problem. Sharp stability and regularity estimates of Mittag–Leffler type are derived, including fractional-domain bounds. Numerical illustrations confirm the predicted decay, and connections with non-autonomous operators, maximal $$L^p$$ -regularity, and weighted AB kernels are outlined. The results place AB-type equations within a functional-analytic framework comparable to the classical theory for Caputo and Volterra models.

37. Landau-Lifshitz-Gilbert方程在螺旋态邻域中局域扰动的时间衰减估计

Time decay estimates for localized perturbations around a helical state for the Landau-Lifshitz-Gilbert equation

作者: Ikkei Shimizu

类别: math.AP

发布日期: 2026-01-25

arXiv ID: 2601.17941v1

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

摘要: We study the dynamics of the Landau–Lifshitz–Gilbert equation with the Dzyaloshinskii–Moriya interaction. The equation admits a family of exact stationary solutions, referred to as helical states, which are periodic in one spatial variable and constant in the others. We investigate the dynamical stability of a helical state with respect to perturbations belonging to suitable Lebesgue and Sobolev spaces. Under a smallness assumption on the initial perturbation, we prove global existence and time decay estimates for solutions, demonstrating that the above helical state is stable. The analysis of the relevant linear operator is carried out via the Bloch–Fourier-wave decomposition, where the eigenvalue problem for the reduced operator is characterized by certain Mathieu equations.

38. 基于广义极化张量消零结构的隐身层合材料设计

Cloaking laminate design based on GPT-vanishing structures

作者: Eleanor Gemida, Mikyoung Lim

类别: math.AP

发布日期: 2026-01-25

arXiv ID: 2601.17932v1

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

摘要: We propose a near-cloaking design which is a lamination of a finite number of layers of isotropic materials. The proposed design is an approximation of a cloaking material obtained by pushing forward a multi-coated structure for which the coating cancels the generalized polarization tensors (GPTs) up to several leading orders. The enhanced cloaking effect achieved by the GPT-vanishing structure permits a coarser microscale requirement and reduces the contrast in the constituent isotropic materials, thereby improving constructibility of cloaking laminates compared with designs based on a non-coated structure.

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

1. 一类双参数双曲松弛系统向不可压Navier-Stokes方程的收敛性#

详细分析#

研究对象和背景#

主要定理或主要结果#

【小扰动情形】$\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = o(1)$#

【大扰动情形】$\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = \mathcal{O}(1)$#

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

a) 核心分析学工具#

b) 关键技术难点突破#

与之前工作的比较#

a) 突破与改进#

2. 基于凸积分框架的低马赫数极限普适性研究#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

3. 关于排斥性无压欧拉-泊松系统的渐近行为#

详细分析#

1. 研究对象和背景#

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

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

a) 核心分析学工具#

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

4. 与之前工作的比较#

a) 突破与改进#

4. Vlasov-Fokker-Planck-Navier-Stokes 系统的光子-粒子定量极限#

详细分析#

研究对象和背景#

主要定理或主要结果#

（I）良备初值情形（Well-prepared initial data）#

（II）尺度良备初值情形（Scaled well-prepared initial data）#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

5. 液晶与拓扑涡度：适定弱解的光滑性#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

a) 核心分析学工具#

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

与之前工作的比较#

a) 突破与改进#

6. 关于各向异性Boussinesq方程在半空间 $\mathbb{R}_+^3$ 中临近静水平衡态解的衰减性#

详细分析#

1. 研究对象和背景#

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

定理 1.1（全局稳定性）#

定理 1.2（衰减性）#

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

a) 核心分析学工具#

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

4. 与之前工作的比较#

a) 突破与改进#

7. Navier-Stokes–Cahn–Hilliard 系统的一类松弛近似方法的合理性证明#

详细分析#

1. 研究对象和背景#

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

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

a) 核心分析学工具#

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

4. 与之前工作的比较#

a) 突破与改进#

8. 一种用于四维血流磁共振成像衍生流动问题的稳定化有限元方法#

详细分析#

1. 研究对象和背景#

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

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

a) 核心分析学工具#

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

4. 与之前工作的比较#

a) 突破与改进#

9. 主动湍流中Tonner-Tu-Swift-Hohenberg方程的整体吸引子及其性质#

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

1. 一类双参数双曲松弛系统向不可压Navier-Stokes方程的收敛性

详细分析

研究对象和背景

主要定理或主要结果

【小扰动情形】 $\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = o(1)$

【大扰动情形】 $\|u_0^{\epsilon,\delta} - u_0\|_{L^\infty} = \mathcal{O}(1)$

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

a) 核心分析学工具

b) 关键技术难点突破

与之前工作的比较

a) 突破与改进

2. 基于凸积分框架的低马赫数极限普适性研究

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

3. 关于排斥性无压欧拉-泊松系统的渐近行为

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

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

a) 核心分析学工具

b) 技术性难题处理技巧

4. 与之前工作的比较

a) 突破与改进

4. Vlasov-Fokker-Planck-Navier-Stokes 系统的光子-粒子定量极限

详细分析

研究对象和背景

主要定理或主要结果

（I）良备初值情形（Well-prepared initial data）

（II）尺度良备初值情形（Scaled well-prepared initial data）

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

5. 液晶与拓扑涡度：适定弱解的光滑性

详细分析

研究对象和背景

主要定理或主要结果

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

a) 核心分析学工具

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

与之前工作的比较

a) 突破与改进

6. 关于各向异性Boussinesq方程在半空间 $\mathbb{R}_+^3$ 中临近静水平衡态解的衰减性

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

定理 1.1（全局稳定性）

定理 1.2（衰减性）

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

a) 核心分析学工具

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

4. 与之前工作的比较

a) 突破与改进

7. Navier-Stokes–Cahn–Hilliard 系统的一类松弛近似方法的合理性证明

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

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

a) 核心分析学工具

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

4. 与之前工作的比较

a) 突破与改进

8. 一种用于四维血流磁共振成像衍生流动问题的稳定化有限元方法

详细分析

1. 研究对象和背景

2. 主要定理或主要结果

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

a) 核心分析学工具

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

4. 与之前工作的比较

a) 突破与改进

9. 主动湍流中Tonner-Tu-Swift-Hohenberg方程的整体吸引子及其性质