生成时间: 2026年07月07日 06:36:26

运行状态: 完成

计划处理论文数量: 55

已完成论文数量: 55

未完成论文数量: 0

重点关注论文数量: 12

了解领域论文数量: 26

不相关论文数量: 17

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

1. 可能退化粘性Hamilton–Jacobi方程的vanishing viscosity逼近的收敛速率

Convergence Rates for Vanishing Viscosity Approximations of Possibly Degenerate Viscous Hamilton–Jacobi Equations

作者: Junkai Qi

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.05186v1

摘要翻译: 我们在环面 $\mathbb{T}^d$ 上研究可能退化的粘性Hamilton–Jacobi方程的vanishing viscosity逼近的定量收敛速率。极限方程含有空间依赖的扩散系数 $a(x) \geq 0$，该系数允许在某些点处消失。在Hamiltonian的标准结构假设下，我们首先证明了逐点收敛速率为 $O(\epsilon |\log \epsilon|)$。接着我们证明：当误差与一光滑概率密度函数作配对测试时，对数损失可被消除，从而得到平均意义下的 $O(\epsilon)$ 收敛速率。证明基于非线性伴随方法、加权Hessian估计以及伴随密度的熵估计。

Comment: 17 pages

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

1. 研究对象和背景

本文研究定义在平坦环面 $\mathbb{T}^n$ 上的退化粘性Hamilton–Jacobi方程的消失粘性逼近的定量收敛速率问题。核心研究对象是如下退化抛物型初值问题：

其中扩散系数 $a \in C^2(\mathbb{T}^n)$ 满足 $a(x) \geq 0$，且允许在非空子集上严格为零（即退化），$H \in C^2(\mathbb{T}^n \times \mathbb{R}^n)$ 是满足均匀凸性条件 $D_{pp}^2 H(x,p) \geq \theta I_n$（$\theta > 0$）及标准增长条件 $|D_x H(x,p)| \leq C(1+|p|^2)$ 的Hamiltonian，$g \in C^2(\mathbb{T}^n)$。该方程描述了在空间非均匀、甚至局部无扩散机制下的最优控制或动力学平均场演化；当 $a \equiv 0$ 时退化为一阶Hamilton–Jacobi方程，而当 $a \geq a_0 > 0$ 时则为非退化粘性方程。其解 $u$ 被理解为粘性解（viscosity solution），这是由Crandall–Lions理论保证的唯一稳定解。

对应的消失粘性正则化族为：

由于 $a(x)+\varepsilon > 0$，对每个 $\varepsilon > 0$，(1.2) 是一致椭圆抛物方程，存在经典 $C^{2,1}$ 解 $u^\varepsilon$；由粘性解稳定性知 $u^\varepsilon \to u$ 一致收敛。但经典理论仅给出定性收敛，无法回答关键定量问题：收敛速度多快？这一问题在PDE正则性理论、最优控制数值逼近、以及与随机微分方程（如小噪声扩散过程）的联系中均具根本意义。例如，在随机最优控制中，$u^\varepsilon$ 对应带白噪声扰动的值函数，$\varepsilon \to 0$ 对应噪声强度趋于零的渐近，收敛速率直接决定近似策略的误差界。

已有工作主要分三类：(i) 非退化情形（$a \geq a_0 > 0$）下，利用标准能量估计或Bernstein方法可得 $O(\sqrt{\varepsilon})$ 速率，这是由线性化算子 $\partial_t - b^\varepsilon \cdot D + (a+\varepsilon)\Delta$ 的最小特征值 $\sim \varepsilon$ 所决定的尺度；(ii) 完全退化情形（$a \equiv 0$）下，Cirant–Goffi [14] 利用非线性伴随方法与熵估计获得 $O(\varepsilon |\log \varepsilon|)$ 点态速率，其关键在于伴随密度 $\sigma^\varepsilon$ 在终点 $t=\tau$ 处的Dirac质量导致的奇异性，表现为 $\|D\sigma^\varepsilon(\cdot,t)\|_{L^1} \sim \varepsilon^{-1}|\log \varepsilon|^{1/2}$；(iii) 周期均匀化等变体中，Qian–Sprekeler–Tran–Yu [16] 得到 $O(\sqrt{\varepsilon})$，Liu–Tran–Yu [17] 在二次Hamiltonian下改进至 $O(\varepsilon |\log \varepsilon|)$。然而，对空间依赖且可能退化的扩散 $a(x) \geq 0$ 这一更贴近物理模型（如各向异性介质中的输运）的情形，既往工作存在显著空白：Bernstein方法失效——因对(1.2)求导产生不可控项 $D a \cdot D^2 u^\varepsilon$；标准伴随估计仅给出 $O(\sqrt{\varepsilon})$；而 $a(x)$ 的零点集结构（如是否为孤立点、低维流形）如何影响速率亦未被系统刻画。本文旨在填补这一空白，建立首个适用于一般 $a \in C^2$, $a \geq 0$ 的定量收敛理论，并揭示退化几何（通过 $|Da|^2 \leq C a$）与对数速率之间的内在联系。

2. 主要定理或主要结果

本文在假设 (1.3) 与 (1.4) 下，证明了两个主定理：

定理 1.1（点态收敛速率）：设 $u^\varepsilon$ 为 (1.2) 的解，$u$ 为 (1.1) 的粘性解，则存在常数 $C > 0$（仅依赖于 $n, T, \theta, C$ 及 $\|g\|_{C^2}, \|a\|_{C^2}, \|H\|_{C^2}$，与 $\varepsilon$ 无关），使得

该速率在 $\varepsilon \to 0^+$ 时为最优可能的点态速率，且对任意退化结构 $a(x) \geq 0$ 成立。

定理 1.2（加权平均收敛速率）：设 $r \in C^\infty(\mathbb{T}^n)$ 满足 $r \geq 0$ 且 $\int_{\mathbb{T}^n} r(x)\,dx = 1$，则存在常数 $C > 0$（同上，且与 $r$ 无关），使得

特别地，若 $r$ 有界远离零（如 $r \geq c > 0$），则右侧简化为 $O(\varepsilon)$；若 $r$ 光滑且紧支，则 $\int r|\log r|$ 有限，故仍为 $O(\varepsilon)$。该结果表明：对光滑概率密度 $r$ 的测试，对数因子完全消失，获得真正的一阶速率。证明中关键参数范围是 $0 < \varepsilon < 1$，且时间区间 $[0,T]$ 固定；所有常数不依赖于 $\varepsilon$，但隐含依赖于 $T$（通过 $\|Du^\varepsilon\|_{L^\infty}$ 的Bernstein界 (1.5)）。文中未讨论 $T \to \infty$ 的大时间行为，亦未给出 $\varepsilon$-独立的显式常数。

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

本文的核心技术障碍源于退化扩散项 $a(x)\Delta u$：当对 (1.2) 求导以分析 $\partial_\varepsilon u^\varepsilon$ 时，线性化源项为 $\Delta u^\varepsilon$，而伴随方法将其转化为对 $\Delta u^\varepsilon$ 的加权积分估计；但标准加权估计仅能给出 $O(\sqrt{\varepsilon})$，因 $\int_0^T \int_{\mathbb{T}^n} |D^2 u^\varepsilon|^2 \sigma^\varepsilon \,dx dt \lesssim 1$ 与 $\int_0^T \int_{\mathbb{T}^n} \frac{\sigma^\varepsilon}{a+\varepsilon} \,dx dt \lesssim \varepsilon^{-1}$ 的乘积导致 $\|\partial_\varepsilon u^\varepsilon\|_{L^\infty} \lesssim \varepsilon^{-1/2}$。为突破此瓶颈，作者构建了三层递进的技术框架：

第一层：退化几何控制引理（Lemma 2.1）。利用 $a \in C^2$, $a \geq 0$ 的结构，证明 $|Da(x)|^2 \leq C a(x)$ 对所有 $x$ 成立。此不等式是全文基石——它将 $Da$ 项“吸收”进 $a$ 中，使后续估计中出现的 $\frac{|Da|^2}{a+\varepsilon}$ 有界，从而避免了 $Da$ 的奇异性。证明基于 $a^{1/2} \in \mathrm{Lip}$ 及零点处 $Da=0$ 的事实，是纯分析的精细构造，而非泛泛的Sobolev嵌入。

第二层：加权二阶导数估计（Proposition 2.5）。为获得 $O(\varepsilon |\log \varepsilon|)$，需证明 $\int_{t_1}^{t_2} \int_{\mathbb{T}^n} \chi_\tau(t) |D^2 u^\varepsilon|^2 \sigma^\varepsilon \,dx dt \lesssim 1 + |\log \varepsilon|$，其中 $\chi_\tau(t) = t(\tau-t)$ 是端点衰减权重，$[t_1,t_2] \subset [\varepsilon, \tau-\varepsilon]$ 是精心选取的时间区间。关键创新在于对方程 $\phi^\varepsilon = \Delta u^\varepsilon$ 应用加权能量估计：

其中 $L^\varepsilon f = f_t + b^\varepsilon \cdot Df - (a+\varepsilon)\Delta f$。难点在于 $2 \chi_\tau Da \cdot D\phi^\varepsilon$ 项：积分分部后生成 $\int \chi_\tau \phi^\varepsilon Da \cdot D\sigma^\varepsilon \,dx dt$，这要求控制 $D\sigma^\varepsilon$。作者引入Fisher信息 $\int A^\varepsilon \frac{|D\sigma^\varepsilon|^2}{\sigma^\varepsilon} \,dx$ 并结合熵估计 $E(t) = \int \sigma^\varepsilon \log \sigma^\varepsilon \,dx$，利用 $E'(t)$ 的结构（见 (2.23)）与 $|Da|^2 \leq C a$，最终导出 (2.16)：

此即连接 $D\sigma^\varepsilon$ 与 $D^2 u^\varepsilon$ 的关键不等式，其代价是引入对数因子。

第三层：伴随密度的时空奇异性刻画（Lemma 2.4）。为闭合上述估计，必须量化 $\sigma^\varepsilon$ 的奇异性。作者将时间反转后的方程 (2.17) 视为具有椭圆常数 $\lambda_0 = \varepsilon$ 的抛物方程，引用文献 [18] 的热核上界，得到 $\log \sigma^\varepsilon(x,t) \lesssim 1 + |\log \varepsilon| + |\log(\tau-t)|$，进而推出熵 $\int \sigma^\varepsilon |\log \sigma^\varepsilon| \,dx \lesssim 1 + |\log \varepsilon|$ 及 $L^1$ 梯度界 $\|D\sigma^\varepsilon(\cdot,t_i)\|_{L^1} \lesssim \varepsilon^{-1}(1 + |\log \varepsilon|)^{1/2}$（见 (2.15)）。这些估计精准捕获了Dirac初始数据在退化扩散下的传播特性，是去除对数因子的前提。

4. 与之前工作的比较

本文与既有工作的比较需从四个维度展开：

1. 与非退化情形（$a \geq a_0 > 0$）的比较：经典工作（如[1–3]）在 $a \geq a_0$ 下仅得 $O(\sqrt{\varepsilon})$，因其伴随密度 $\sigma^\varepsilon$ 的奇异性被 $a_0$ 控制，$\int_0^T \int \frac{\sigma^\varepsilon}{a+\varepsilon} \,dx dt \lesssim a_0^{-1}$。本文在 $a \geq 0$ 下将速率提升至 $O(\varepsilon |\log \varepsilon|)$，代价是引入对数因子及更精细的熵/Fisher信息分析。若强行套用非退化方法，会得到 $\|\partial_\varepsilon u^\varepsilon\|_{L^\infty} \lesssim \varepsilon^{-1/2}$，远逊于本文结果。

2. 与完全退化情形（$a \equiv 0$）的比较：Cirant–Goffi [14] 在 $a \equiv 0$ 下证明 $O(\varepsilon |\log \varepsilon|)$，但其证明依赖于 $a \equiv 0$ 时方程的特殊结构（如 $b^\varepsilon$ 的特定形式），且未处理 $Da$ 项。本文将 [14] 的框架推广至 $a(x) \geq 0$，核心创新是 Lemma 2.1（$|Da|^2 \leq C a$）及由此衍生的加权估计 (2.16)，这在 $a \equiv 0$ 时自动成立但无意义。因此，本文不仅是推广，更是对退化几何与速率关系的首次系统揭示。

3. 与周期均匀化工作的比较：Qian–Sprekeler–Tran–Yu [16] 在均匀化背景下得 $O(\sqrt{\varepsilon})$，其困难在于振荡系数；Liu–Tran–Yu [17] 在二次Hamiltonian下得 $O(\varepsilon |\log \varepsilon|)$，但要求 $a$ 为常数且 $H$ 二次。本文设定更一般（$a(x)$ 任意 $C^2$ 非负，$H$ 仅需均匀凸），且不依赖周期结构，适用范围更广。代价是未利用均匀化特有的细尺度振荡补偿效应。

4. 与加权平均估计的比较：Han–Tu [4] 对状态约束方程得 $O(\sqrt{\varepsilon})$，Camilli–Goffi–Mendico [11] 用伴随法得 $L^p$ 速率，但均未触及 $O(\varepsilon)$。本文定理 1.2 的 $O(\varepsilon)$ 是首个针对退化粘性HJ方程的线性速率，其本质是用光滑 $r$ 替换Dirac质量 $\delta_{x_0}$，消除了伴随密度在 $t=T$ 的奇异性，从而规避了对数损失。文中明确指出此改进“removes the endpoint singularity”，与 [14] 的 $O(\varepsilon |\log \varepsilon|)$ 形成互补：前者是点态最优，后者是平均最优。从提供内容中无法确认本文结果是否可推广至 $L^p$ 或 Sobolev 范数下的收敛速率。

2. 关于零质量初值条件下对流-扩散方程整体解的最优衰减率

On the optimal decay rate of global solutions to the convection-diffusion equation with zero-mass initial data

作者: Yi C. Huang, Ryunosuke Kusaba

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.05158v1

摘要翻译: 我们研究零质量初值条件下对流-扩散方程整体解的最优衰减率。我们在零质量条件下建立了整体解的改进衰减估计。此外，我们推导出整体解的自相似渐近剖面。该结果给出了达到改进衰减率的充要条件，并表明该衰减率是最优的。

Comment: 21 pages

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

1. 研究对象和背景

本文研究如下 Cauchy 问题所描述的对流-扩散方程的全局解的长时间渐近行为：

其中 $u: [0,\infty) \times \mathbb{R}^n \to \mathbb{R}$ 为未知标量函数，$a = (a_1,\dots,a_n) \in \mathbb{R}^n \setminus \{0\}$ 为固定非零向量，非线性项 $f:\mathbb{R}\to\mathbb{R}$ 是次数 $p > 1$ 的（正）齐次函数，取形如 $f(\xi) = \xi^p$、$|\xi|^p$ 或 $|\xi|^{p-1}\xi$；当 $f(\xi)=\xi^p$ 时可设 $p$ 为整数。该方程建模了具有质量守恒律的流体动力学过程：

其中线性扩散项 $-\Delta u$ 与非线性对流项 $a \cdot \nabla f(u)$ 分别刻画沿方向 $a$ 的扩散与非线性输运。本文聚焦于零质量初值条件 $M_0(u_0) = 0$ 下的衰减行为——这一条件在不可压 Navier–Stokes 方程中天然满足（因 $\nabla \cdot u_0 = 0$ 蕴含 $M_0(u_0) = 0$），故其渐近结构对理解高维非线性抛物系统中对流与扩散的竞争机制至关重要。

已有理论将 $p$ 相对于临界指数 $1 + 1/n$ 进行分类：当 $p < 1 + 1/n$ 时，解以更快速率衰减，渐近剖面 $A_{1,\text{sub}}$ 满足退化椭圆型方程 $\partial_t A - \Delta_\perp^a A = a \cdot \nabla f(A)$；当 $p = 1 + 1/n$ 时，渐近剖面 $A_{1,\text{cri}}$ 满足原方程自身；当 $p > 1 + 1/n$ 时，渐近剖面退化为热核 $A_{1,\text{sup}}(t) = M_0(u_0) G_t$，即线性主导。对应衰减率为：

且已知这些速率在 $M_0(u_0) \neq 0$ 时最优。然而，当 $M_0(u_0) = 0$ 时，所有前述一阶渐近剖面 $A_{1,*}$ 恒为零，导致原有衰减估计失效，亟需建立新的尺度。此前工作（如 [11]）仅在 $p = 1 + 1/n$ 且附加高斯权 $e^{|x|^2/8} u_0 \in L^2$ 条件下得到 $t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2} + \varepsilon}$ 型估计，未达精确阶且适用范围窄。本文旨在填补这一空白：在零质量条件下，确定全局解的最优代数衰减率，并构造其自相似渐近剖面，从而完整刻画 $M_0(u_0) = 0$ 时对流-扩散方程的渐近主导结构。

2. 主要定理或主要结果

本文在零质量初值 $M_0(u_0) = 0$ 下建立了三个核心定理：

定理 2.1（上界估计）：设 $p > 1 + \frac{1}{n+1}$，$u_0 \in L^1_1(\mathbb{R}^n) \cap L^\infty(\mathbb{R}^n)$ 满足 $M_0(u_0) = 0$，$u \in X$ 为 (P) 的全局解。则： (1) 若 $p > 1 + \frac{1}{n}$，则对任意 $q \in [1,\infty]$，存在常数 $C_0 > 0$ 使得

(2) 若 $p \le 1 + \frac{1}{n}$，则存在 $\varepsilon_0 > 0$ 与 $C_1 > 0$（独立于 $\|u_0\|_{L^1_1 \cap L^\infty}$），使得当 $\|u_0\|_{L^1_1 \cap L^\infty} < \varepsilon_0$ 时，对任意 $q \in [1,\infty]$ 有

定理 2.3（自相似渐近剖面）：在同设定下，若 $u$ 满足先验控制 $\sup_{q \in [1,\infty]} \sup_{t>0} (1+t)^{\frac{n}{2}(1 - \frac{1}{q}) + \frac{1}{2}} \|u(t)\|_{L^q} < \infty$，则对任意 $q \in [1,\infty]$，有

其中渐近剖面为

且 $A(t)$ 具有自相似结构 $A(t) = t^{-1/2} \delta_t A(1)$，其中 $\delta_t \phi(x) = t^{-n/2} \phi(t^{-1/2}x)$。

定理 2.4（最优性与充要条件）：在定理 2.3 条件下，对任意 $q \in [1,\infty]$，有

且 $\|A(1)\|_{L^q} > 0$ 当且仅当存在 $j \in \{1,\dots,n\}$ 使得

特别地，该结果证实了命题 1.7 中估计 (1.8) 在 $M_0(u_0)=0$ 时的最优性，因此时 $A_{2,\text{sub}}(t) \equiv 0$。

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

证明的核心障碍在于：零质量条件 $M_0(u_0)=0$ 导致标准热半群 $e^{t\Delta}u_0$ 的主项消失，迫使必须从 Duhamel 公式

中提取更高阶贡献。作者的关键创新点在于构建并严格控制一个二阶自相似剖面 $A(t)$，其本质是将非线性项 $f(u)$ 视为源项后，在热流作用下的首阶响应：

1. 加权 $L^1$ 初值的精细热估计：利用 $u_0 \in L^1_1$ 及 $M_0(u_0)=0$，通过分部积分与热核梯度估计，获得

   $$ \|e^{t\Delta}u_0\|_{L^q} \lesssim t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}}, \quad \|\nabla e^{t\Delta}u_0\|_{L^q} \lesssim t^{-\frac{n}{2}(1 - \frac{1}{q}) - 1}. $$

   此为定理 2.1 的基础，且比文献 [11] 的 $\varepsilon$-损失更优。

2. 非线性项的渐近积分重构：定义 $\psi_0 = \int_0^\infty f(u(\tau)) \, d\tau$，并证明其收敛性。关键引理（见附录）表明：在 $p > 1 + \frac{1}{n+1}$ 下，由定理 2.1 的先验衰减可推出 $\|f(u(t))\|_{L^1} \lesssim t^{-\frac{n(p-1)}{2} - \frac{p}{2}}$，其积分在 $t\to\infty$ 时绝对收敛，从而 $\psi_0 \in L^1$。这使 $A(t)$ 中的 $M_0(\psi_0)$ 有定义。

3. Duhamel 项的自相似分解：将非线性积分拆为

   $$ \int_0^t a \cdot \nabla e^{(t-\tau)\Delta} f(u(\tau)) \, d\tau = a \cdot \nabla G_t \int_0^t M_0(f(u(\tau))) \, d\tau + R(t), $$

   其中余项 $R(t)$ 包含 $f(u(\tau))$ 的中心矩修正。利用 $M_0(u_0)=0$ 及 $f$ 的齐次性，证明 $\int_0^t M_0(f(u(\tau))) \, d\tau = M_0(\psi_0) + o(1)$，进而

   $$ a \cdot \nabla G_t \int_0^t M_0(f(u(\tau))) \, d\tau = a M_0(\psi_0) \cdot \nabla G_t + o(t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}}). $$

   结合 $e^{t\Delta}u_0$ 的一阶展开 $-M_1(u_0) \cdot \nabla G_t + o(t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}})$，即得 $A(t)$ 的显式形式。

4. 唯一性与最优性论证：利用 $A(t)$ 的自相似性 $A(t) = t^{-1/2} \delta_t A(1)$，证明任何具有相同衰减率的自相似剖面必与 $A(t)$ 重合；再通过直接计算 $\|A(t)\|_{L^q} = t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}} \|A(1)\|_{L^q}$，结合定理 2.3 的收敛性，立即导出定理 2.4 的极限等式及非零条件。

4. 与之前工作的比较

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

1. 衰减速率的精确化与范围扩展：

• 文献 [11] 在 $p = 1 + \frac{1}{n}$ 且 $e^{|x|^2/8}u_0 \in L^2$ 下仅得 $t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2} + \varepsilon}$，本文在更弱的 $L^1_1 \cap L^\infty$ 初值下取得无 $\varepsilon$ 损失的精确阶 $t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}}$，并将适用 $p$ 范围从单点 $p=1+1/n$ 扩展至 $p > 1 + \frac{1}{n+1}$（注意 $1 + \frac{1}{n+1} < 1 + \frac{1}{n}$）。
• 对 $p \le 1 + \frac{1}{n}$ 的情形，本文首次给出小初值下的统一衰减率，而此前工作（如 [9,10,3]）在 $M_0(u_0)=0$ 时未处理此范围。

2. 渐近剖面的结构性突破：

• 此前所有渐近剖面（$A_{1,*}, A_{2,*}$）均依赖 $M_0(u_0) \neq 0$，在零质量时恒为零。本文提出的 $A(t) = -(M_1(u_0) - a M_0(\psi_0)) \cdot \nabla G_t$ 是首个非平凡的、由初值一阶矩 $M_1(u_0)$ 与非线性累积量 $M_0(\psi_0)$ 共同驱动的剖面，揭示了零质量下对流效应如何通过非线性积分反馈至主导项。
• 该剖面与 Navier–Stokes 方程中的 Oseen 流剖面 $-M_1(u_0) \cdot \nabla G_t$ 形式一致，但多出修正项 $a M_0(\psi_0) \cdot \nabla G_t$，定量刻画了非线性对流对线性扩散主导结构的扰动，为类比不可压流提供了新依据。

3. 最优性的完备刻画：

• 文献 [20,27] 在 $M_0(u_0) \neq 0$ 下确立了各阶渐近展开的最优性，但对 $M_0(u_0)=0$ 的情形，命题 1.7 中的估计 (1.8) 的最优性长期悬而未决。本文定理 2.4 直接证明：当 $M_0(u_0)=0$ 时，$A_{2,\text{sub}}(t) \equiv 0$，故 (1.8) 的衰减率 $t^{-\frac{n}{2}(1 - \frac{1}{q}) - \frac{1}{2}}$ 实际对应于 $\|u(t)\|_{L^q}$ 的最优速率，从而彻底解决了该开放问题。
• 同时，本文给出了 attainability 的充要条件（式 (2.6)），指出即使 $M_1(u_0) = 0$，只要 $M_0(\psi_0) \neq 0$ 且 $a_j \neq 0$，仍可达到最优衰减，这比单纯要求 $M_1(u_0) \neq 0$ 更精细，体现了非线性反馈的核心作用。

代价方面，本文要求初值属于加权空间 $L^1_1$（而非仅 $L^1$），这是提取 $M_1(u_0)$ 所必需；且对 $p \le 1 + \frac{1}{n}$ 需小初值假设，这与该参数区间的强非线性不稳定性一致。从提供内容中无法确认本文方法能否推广至 $p \le 1 + \frac{1}{n+1}$ 的情形。

3. 不可压旋转Navier–Stokes方程整体解的Sharp衰减估计

Sharp decay estimates for global solutions to the incompressible rotating Navier–Stokes equations

作者: Mikihiro Fujii, Yang Li, Jiang Xu

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.05237v1

摘要翻译: 本文研究三维不可压旋转Navier–Stokes方程，并建立其整体解的sharp $L^p$ 衰减估计。我们指出：对 $2<p<\infty$，最优 $L^p$ 衰减率严格快于通过在 $L^2$ unitary identity 与 $L^\infty$ dispersive estimates 之间插值得到的已有结果中的衰减率，尽管端点情形已被证明是sharp的。此外，通过对某一特定初值数据的下界估计，我们也证明了衰减率的最优性。其内在机制在于由Coriolis力所引发的振荡积分的各向异性退化。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压旋转Navier-Stokes方程的初值问题：

其中 $u(t,x) : [0,\infty) \times \mathbb{R}^3 \to \mathbb{R}^3$ 为未知速度场，$P(t,x)$ 为压力，$u_0 \in L^1(\mathbb{R}^3) \cap \dot{H}^{1/2}(\mathbb{R}^3)$ 满足 $\mathrm{div}\, u_0 = 0$，$\Omega \in \mathbb{R}$ 为绕竖直方向 $e_3 = (0,0,1)$ 的旋转角速度。该模型刻画了快速旋转流体（如大气与海洋动力学）中科里奥利力主导的耗散-色散平衡机制，其线性化算子为 Stokes–Coriolis 算子 $-\Delta + \Omega \mathbb{P} e_3 \times \mathbb{P}$，其中 $\mathbb{P} = I - \nabla (-\Delta)^{-1} \mathrm{div}$ 为 Helmholtz 投影。由 Hieber–Shibata [17] 给出的显式半群表示为

其中 Fourier 乘子 $P_\pm(\xi) = \frac{1}{2}(I \pm i R(\xi))$，$R(\xi) = |\xi|^{-1} \begin{pmatrix} 0 & -\xi_3 & \xi_2 \\ \xi_3 & 0 & -\xi_1 \\ -\xi_2 & \xi_1 & 0 \end{pmatrix}$。关键物理与数学特征在于色散演化群 $\{ e^{i\tau D_3 / |D|} \}_{\tau \in \mathbb{R}}$ 所诱导的振荡积分结构，其相位函数 $\Psi(\xi) = \xi_3 / |\xi|$ 具有各向异性退化性：其 Hessian 矩阵的行列式满足

该退化非均匀——在水平方向 $|\xi_h|$ 小或垂直方向 $|\xi_3|$ 小的区域，退化程度不同，导致经典 Littman 引理所给出的 $L^\infty$ 衰减率 $(1+|\tau|)^{-1}$ 仅反映最坏点处的退化阶，而未计入退化集的测度收缩效应。此前工作（如 [10, 22, 23]）通过 $L^2$ 守恒性 $\|e^{i\tau D_3/|D|} \phi\|_{L^2} = \|\phi\|_{L^2}$ 与 $L^\infty$ 色散估计 $\|e^{i\tau D_3/|D|} \phi\|_{L^\infty} \leq C (1+|\tau|)^{-1}$ 插值得到插值型衰减率

该速率在端点 $p=2$（能量守恒）与 $p=\infty$（Koh–Lee–Takada [23] 证明其尖锐性）处确为最优，但对 $2 < p < \infty$ 是否最优长期悬而未决。已有非线性结果（如 [10, 29]）均以此插值速率作为色散贡献，导出解的 $L^p$ 衰减形如 $t^{-\frac{3}{2}(1-1/p)-m/2} (1+|\Omega|t)^{-(1-2/p)}$，其物理含义是将整个频率环 $\{1/4 \leq |\xi| \leq 4\}$ 视为具有统一退化阶的“最差情形”，忽略了退化集 $\{ |\det \nabla^2 \Psi| \simeq 2^k \}$ 的体积随 $k$ 指数衰减这一关键几何事实。本文旨在填补这一理论空白：严格证明对 $2<p<\infty$，真实衰减速率严格快于插值速率，其根源在于相位函数的各向异性退化结构——退化越强的频率区域体积越小，该体积补偿效应在 $L^p$ 插值中产生额外衰减因子。

2. 主要定理或主要结果

本文核心结论分为线性与非线性两部分，全部建立在初始数据 $u_0 \in L^1(\mathbb{R}^3) \cap \dot{H}^{1/2}(\mathbb{R}^3)$、$\mathrm{div}\, u_0 = 0$ 及大旋转参数 $|\Omega| \geq \Omega_0(u_0)$ 假设下。

线性部分（Theorem 3.2）： 对任意 $m \in \mathbb{N} \cup \{0\}$，Stokes–Coriolis 半群 $T_\Omega(t) u_0$ 满足：

1. 一般衰减： 若 $u_0 \in L^1$，则
   $$ \|\nabla^m T_\Omega(t) u_0\|_{L^p} = o\left( t^{-\frac{3}{2}(1-1/p)-m/2} D_p(|\Omega| t) \right), \quad t \to \infty, $$
   其中 $D_p(\tau)$ 定义为
   $$ D_p(\tau) := \begin{cases} (1 + |\tau|)^{-\frac{3}{2}(1-2/p)}, & 2 \leq p < 4, \\ (1 + |\tau|)^{-3/4} (\log(e + |\tau|))^{1/4}, & p = 4, \\ (1 + |\tau|)^{-(1-1/p)}, & 4 < p \leq \infty. \end{cases} $$
2. 增强衰减： 若进一步有 $|x| u_0(x) \in L^1$，则存在常数 $C = C(p,m)$ 使得
   $$ \|\nabla^m T_\Omega(t) u_0\|_{L^p} \leq C t^{-\frac{3}{2}(1-1/p)-m/2 + 1/2} D_p(|\Omega| t) \| |x| u_0 \|_{L^1}, \quad \forall t > 0. $$
   特别地，对 $p=4$，$L^{4,\infty}$ 范数可消除对数修正：
   $$ \|\nabla^m T_\Omega(t) u_0\|_{L^{4,\infty}} \leq C t^{-\frac{3}{2}(1-1/4)-m/2 + 1/2} (1 + |\Omega| t)^{-3/4} \| |x| u_0 \|_{L^1}. $$

非线性部分（Theorems 1.1 & 1.2）： 方程 (1.1) 存在唯一温和解 $u \in C([0,\infty); \dot{H}^{1/2}) \cap C((0,\infty); L^\infty)$，且满足：

• 一般衰减（Theorem 1.1）： 对 $2 \leq p \leq \infty$，
  $$ \|\nabla^m u(t)\|_{L^p} = \begin{cases} o\left( t^{-\frac{3}{2}(1-1/p)-m/2} (1 + |\Omega|t)^{-\frac{3}{2}(1-2/p)} \right), & 2 \leq p < 4, \\ o\left( t^{-\frac{3}{2}(1-1/4)-m/2} (1 + |\Omega|t)^{-3/4} (\log(e + |\Omega|t))^{1/4} \right), & p = 4, \\ o\left( t^{-\frac{3}{2}(1-1/p)-m/2} (1 + |\Omega|t)^{-(1-1/p)} \right), & 4 < p \leq \infty. \end{cases} $$
• 增强衰减（Theorem 1.2）： 若 $|x| u_0 \in L^1$，则
  $$ \|\nabla^m u(t)\|_{L^p} \leq \begin{cases} C t^{-\frac{3}{2}(1-1/p)-m/2 + 1/2} (1 + |\Omega|t)^{-\frac{3}{2}(1-2/p)}, & 2 \leq p < 4, \\ C t^{-\frac{3}{2}(1-1/4)-m/2 + 1/2} (1 + |\Omega|t)^{-3/4} (\log(e + |\Omega|t))^{1/4}, & p = 4, \\ C t^{-\frac{3}{2}(1-1/p)-m/2 + 1/2} (1 + |\Omega|t)^{-(1-1/p)}, & 4 < p \leq \infty, \end{cases} $$
  且对 $p=4$ 有 $\|\nabla^m u(t)\|_{L^{4,\infty}} \leq C t^{-\frac{3}{2}(1-1/4)-m/2 + 1/2} (1 + |\Omega|t)^{-3/4}$。

最优性验证（Theorem 1.4）： 对 $\Omega = 1$，存在紧支光滑初值 $u_0(x) = \varepsilon e^{-|x|^2/4} (x_2, -x_1, 0)$（$\varepsilon > 0$ 充分小），使得解满足下界估计：

以及 $\|u(t)\|_{L^{4,\infty}} \geq c t^{-\frac{3}{2}(1-1/4)-1/2} (1 + t)^{-3/4}$，对所有 $t \geq T(p)$ 成立。这严格证实了前述上界估计的尖锐性。

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

本文的核心技术障碍在于：经典插值方法无法捕获振荡积分 $e^{i\tau D_3/|D|} \phi$ 在 $L^p$ ($2<p<\infty$) 中的真实衰减行为，因其隐含假设整个支撑集具有相同退化阶，而实际相位 $\Psi(\xi) = \xi_3/|\xi|$ 的 Hessian 行列式 $\det \nabla^2 \Psi = |\xi_h|^2 \xi_3 / |\xi|^9$ 在频率空间中呈现各向异性退化——退化最强处（如 $\xi_h \approx 0$ 或 $\xi_3 \approx 0$）的集合测度极小。作者的真正创新点在于构建了一套分层退化分析框架，将频率支撑集 $\{1/4 \leq |\xi| \leq 4\}$ 按退化程度 $|\det \nabla^2 \Psi| \simeq 2^k$ 进行 dyadic 分解，并分别估计每一层的贡献，最终通过求和揭示体积补偿效应。

具体而言，证明框架包含三个关键步骤：

第一步：各向异性频率分解（Section 3.1, Lemma 3.4）。 引入截断函数 $\rho_h(\xi) = \varrho(|\xi_h|/|\xi|)$, $\rho_v(\xi) = \varrho(\xi_3/|\xi|)$, $\rho_0 = 1 - \rho_h - \rho_v$，将振荡积分分解为三部分：

• 非退化部分 $I_0$：在 $\mathrm{supp}\, a_0 = \mathrm{supp}\, \hat{\phi} \rho_0$ 上，$|\det \nabla^2 \Psi| \gtrsim 1$，直接应用标准平稳相位引理（Lemma 2.4）得 $\|I_0(\tau)\|_{L^\infty} \lesssim |\tau|^{-3/2}$；
• 水平退化部分 $I^h_k$：对应 $|\xi_h| \simeq 2^{k/2}$，通过变量替换 $(\eta_h, \eta_3) = (2^{-k/2} \xi_h, \xi_3)$ 将相位缩放为非退化形式 $\Psi^h_{\pm,k}(\eta)$，并验证其 Hessian 行列式一致有下界，从而应用 Lemma 2.4 得 $\|I^h_k(\tau)\|_{L^\infty} \lesssim \min\{2^k, 2^{-k/2} |\tau|^{-3/2}\}$；
• 竖直退化部分 $I^v_k$：对应 $|\xi_3| \simeq 2^k$，类似地通过 $(\eta_h, \eta_3) = (\xi_h, 2^{-k} \xi_3)$ 缩放，得 $\|I^v_k(\tau)\|_{L^\infty} \lesssim \min\{2^k, 2^{-k/2} |\tau|^{-3/2}\}$。 此分解精准刻画了退化集的几何结构：$|\{ |\det \nabla^2 \Psi| \simeq 2^k \}| \simeq 2^k$，即退化越强（$k \to -\infty$），集合体积越小。

第二步：分层 $L^p$ 插值与求和（Lemma 3.5）。 对每一层 $I^\#_k$，利用 $L^2$-守恒 $\|I^\#_k\|_{L^2} \lesssim 2^{k/2}$ 与上述 $L^\infty$ 估计插值得到

关键洞察在于，对 $2 < p < 4$，指数 $2/p - 1/2 > 0$，求和 $\sum_{k \leq 1} \|I^\#_k\|_{L^p}$ 由非退化部分（$k$ 大）主导，结果为 $|\tau|^{-\frac{3}{2}(1-2/p)}$；对 $4 < p \leq \infty$，指数 $2/p - 1/2 < 0$，求和集中在过渡尺度 $2^k \simeq |\tau|^{-1}$，结果为 $|\tau|^{-(1-1/p)}$；对 $p=4$，两指数均为零，出现对数发散，需引入弱 $L^{4,\infty}$ 范数或 Lorentz 插值（Remark 1.3）来控制。

第三步：分布函数估计与 $L^{4,\infty}$ 优化（Proof of Lemma 3.5, Step 2）。 为消除 $p=4$ 的对数项，作者采用分布函数技巧：对给定 $\lambda > 0$，定义集合 $E(\tau, \lambda) = \{ k \leq 1 : \min\{2^k, 2^{-k/2}|\tau|^{-3/2}\} > \delta \lambda \}$，则 $\{ |I^\#(\tau,x)| > \lambda \}$ 被 $\sum_{k \in E} I^\#_k$ 的水平集控制。利用 $|E(\tau, \lambda)| \lesssim \log(e + |\tau|)/\log(1/\delta)$ 及每项 $L^\infty$ 界，最终得到 $\|I^\#(\tau)\|_{L^{4,\infty}} \lesssim |\tau|^{-3/4}$。此技巧将对数发散转化为弱范数下的干净衰减，是处理临界情形的关键工具。

综上，本文的分析框架超越了传统插值，是一种融合各向异性频率分解、分层平稳相位估计、体积补偿求和及分布函数控制的全新方法，首次系统揭示了旋转流体中色散效应的精细 $L^p$ 结构。

4. 与之前工作的比较

本文与既有工作的根本区别在于对色散效应衰减速率的定量刻画，其推进体现在假设、速率、适用范围与技术深度四个维度：

1. 衰减速率的实质性推进： 此前所有工作（[10, 22, 23, 29]）均基于插值估计 $\|e^{i\tau D_3/|D|} \phi\|_{L^p} \lesssim (1+|\tau|)^{-(1-2/p)}$，由此导出非线性解的 $L^p$ 衰减为 $t^{-\frac{3}{2}(1-1/p)-m/2} (1+|\Omega|t)^{-(1-2/p)}$。本文证明该速率对 $2<p<\infty$ 严格非最优：当 $2 \leq p < 4$ 时，真实速率快 $|\Omega|t$ 的 $\frac{1}{2}(1-2/p)$ 阶（例如 $p=3$ 时，新速率 $(1+|\Omega|t)^{-1/4}$ vs 旧速率 $(1+|\Omega|t)^{-1/3}$）；当 $4 < p \leq \infty$ 时，新速率 $(1+|\Omega|t)^{-(1-1/p)}$ 仍快于旧速率 $(1+|\Omega|t)^{-(1-2/p)}$（例如 $p=6$ 时，新速率 $(1+|\Omega|t)^{-5/6}$ vs 旧速率 $(1+|\Omega|t)^{-2/3}$）。这一差距源于对退化集测度的精确利用，是质的飞跃。

2. 假设条件的等价性与普适性： 本文的初值正则性假设 $u_0 \in L^1 \cap \dot{H}^{1/2}$ 与 [10, 29] 一致，未作任何加强；大旋转参数 $|\Omega| \geq \Omega_0(u_0)$ 的要求亦为标准（见 [7, 9]）。因此，速率提升无需额外代价，是纯粹的分析深化。值得注意的是，本文首次建立了 $L^p$ 衰减的下界（Theorem 1.4），而此前文献（包括 [23]）仅对端点 $p=\infty$ 给出下界，对 $2<p<\infty$ 的下界完全空白。这彻底解决了该问题的最优性争议。

3. 适用范围的扩展与临界处理： 本文不仅给出 $L^p$ 估计，更通过 $L^{4,\infty}$ 范数和 Lorentz 插值（Remark 1.3）精细刻画了 $p=4$ 这一临界情形，指出对数修正 $\log^{1/4}(e+|\Omega|t)$ 是 $L^4$ 框架下的固有现象，但可在弱 $L^{4,\infty}$ 框架中完全消除。这种对临界空间的区分处理，在 [10, 29] 中完全缺失，后者仅给出统一的 $L^p$ 估计。

4. 技术局限与未覆盖情形： 本文结果依赖于大 $\Omega$ 假设，对小 $\Omega$ 或 $\Omega=0$（即经典 Navier-Stokes）不适用，这与 [7, 9] 的全局适定性框架一致。此外，证明基于初值 $u_0 \in L^1$，未涉及更广的 $\dot{B}^{-1}_{\infty,\infty}$ 类初值（如 [12]），也未讨论分数阶 Laplacian 或 MHD 等推广模型（虽在引言提及 [1, 21]，但未将其衰减结果改进）。从提供内容中无法确认本文方法能否直接推广至这些情形。最后，本文聚焦于 $L^p$ 衰减，未涉及解的渐近形态（如 Fujigaki–Miyakawa 展开 [11] 的改进）或长时间稳定性，这些是未来工作方向。

4. 无磁粘性和角粘性的三维不可压磁微极流体在条形区域中竖直带电稳态磁场的稳定性

Stability of vertically charged steady magnetic field in 3D incompressible magneto-micropolar fluids without magnetic and angular viscosity in a strip domain

作者: Zefu Feng, Guangyi Hong, Jiahong Wu, Kun Zhao

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.05166v1

摘要翻译: 本文旨在研究无磁粘性和角粘性的三维不可压磁微极方程在条形区域上的正则性与稳定性问题。磁微极系统建模了存在磁场时的导电流体微极流体。由于缺乏磁扩散和角耗散，甚至连小初值的整体适定性结果都难以证明，更不用说一般大初值的整体正则性了。本文构造了一个稳态构型，围绕该构型可证明任意扰动均具有整体正则性与稳定性。更精确地说，任一靠近垂直于水平边界的稳态磁场的小扰动，均诱导出唯一的全局经典解。此外，该解在时间趋于无穷时以几乎指数速率收敛至该稳态。这些结果似乎是关于本文所考虑的磁微极方程的首批严格整体结果。

Comment: 无

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

1. 研究对象和背景

本文研究的是定义在条形区域 $\Omega = \mathbb{R}^2 \times (0,1)$ 上的三维不可压磁微极流体方程组，其核心特征是磁扩散系数 $\nu = 0$ 与角粘性系数 $\kappa_1 = \kappa_2 = \kappa_3 = 0$。该模型由如下初边值问题（IBVP）描述：

其中 $\tilde{u} = (\tilde{u}_1, \tilde{u}_2, \tilde{u}_3)$ 为流体速度，$\tilde{p}$ 为压力，$\tilde{b} = (\tilde{b}_1, \tilde{b}_2, \tilde{b}_3)$ 为磁场，$\tilde{w} = (\tilde{w}_1, \tilde{w}_2, \tilde{w}_3)$ 为微旋转角速度；$\mu > 0$ 为动力粘性系数，$\zeta > 0$ 为涡旋粘性系数。该系统是经典不可压磁流体动力学（MHD）方程与微极流体角动量方程的耦合：当 $\zeta = 0$ 且 $\tilde{w} \equiv 0$ 时，退化为 MHD；当 $\tilde{b} \equiv 0$ 时，则退化为标准微极流体方程。

问题的关键难点在于方程 (1.1b) 缺乏 Laplace 型磁扩散项 $\nu \Delta \tilde{b}$，而 (1.1c) 缺乏角动量耗散项 $\kappa_1 \Delta \tilde{w} + (\kappa_2 + \kappa_3) \nabla (\nabla \cdot \tilde{w})$。这导致两个根本性障碍：(i) 磁场 $\tilde{b}$ 的 Sobolev 范数（如 $\|\nabla \tilde{b}\|_{L^2}$）无法通过能量方法进行先验控制，其演化由输运-拉伸方程 $\partial_t \tilde{b} = \nabla \tilde{u}^T \tilde{b}$ 驱动，本质上是双曲型；(ii) 速度 $\tilde{u}$ 与角速度 $\tilde{w}$ 通过反对称算子 $\nabla \times$ 强耦合（见 (1.1a) 和 (1.1c)），即 $\zeta \nabla \times \tilde{w}$ 作为 $\tilde{u}$ 的源项，而 $\zeta \nabla \times \tilde{u}$ 又驱动 $\tilde{w}$ 的演化。这种结构在常规欧拉坐标系下会引发导数损失（loss of derivatives）：对 $\tilde{w}$ 直接做高阶能量估计需控制 $\|\nabla \times \tilde{u}\|_{H^k}$，但 $\nabla \times \tilde{u}$ 的正则性又依赖于 $\tilde{u}$ 在更高阶范数下的控制，形成闭合困难。

在数学上，这一系统的全局适定性远比全耗散情形（$\nu, \kappa_i > 0$）困难。已有工作表明：对于全耗散二维/三维磁微极系统，Rojas-Medar-Boldrini [27]、Xiang-Yang [35] 等已建立解的存在唯一性；Yamazaki [37] 处理了二维零角粘性情形；Lin-Xiang [21] 则在二维全耗散框架下研究了围绕平衡态 $(0, e_1, 0)$ 的稳定性。然而，对于 (1.1) 这类双重退化（零磁扩散+零角粘性）的三维系统，即使小初值，传统能量方法亦告失败——因为 $\|\nabla \tilde{b}\|_{L^2}$ 的增长会破坏小初值框架，而 $\tilde{w}$ 的导数损失又使 $\tilde{u}$ 的耗散无法有效传递。因此，一个自然的问题是：是否存在某种非平凡的稳态背景，使得围绕其的小扰动能获得全局存在性与稳定性？

本文选取的稳态背景是垂直于水平边界 $x_3 = 0, 1$ 的恒定磁场，即 $\tilde{b}_{\text{eq}} = e_3 = (0,0,1)$。这是一个 (1.1) 的平凡稳态解（因 $\nabla \tilde{b}_{\text{eq}} = 0$，故 $\tilde{b}_{\text{eq}} \cdot \nabla \tilde{u} = \tilde{u} \cdot \nabla \tilde{b}_{\text{eq}} = 0$）。物理上，该设定对应于条形域中施加垂直方向的外磁场，是实验室等离子体与微流控器件中的典型构型。本文试图填补的空白正是：在三维、无磁扩散、无角粘性、强 $\tilde{u}$-$\tilde{w}$-$\tilde{b}$ 耦合的极端条件下，证明围绕该物理相关稳态的小扰动具有全局经典解、指数型衰减及结构稳定性——这是该模型首个严格成立的全局结果。

2. 主要定理或主要结果

定理 2.1 给出了本文的核心结论：设初边值问题 (1.7)（即 Lagrangian 坐标系下的简化系统）的初始数据满足相容性条件，且存在整数 $N \geq 4$ 和充分小常数 $\varepsilon_0 > 0$，使得初始能量满足 $E_{2N}(0) + F_{2N}(0) \leq \varepsilon_0$，其中能量泛函 $E_{2N}(t)$ 与 $F_{2N}(t)$ 定义为

则存在唯一的全局经典解 $(\eta, u, p, w)$ 满足

对任意 $t > 0$ 成立，其中 $D_{2N}(t)$ 为耗散泛函，$J_{2N}(t) := \|u\|_{4N+1}^2 + \|\nabla p\|_{4N-1}^2 + \|w\|_{4N}^2$，$\theta > 0$ 为任意小正数。由此可直接推出关键衰减估计（推论 2.1）：

该衰减速率随 $N$ 增大而加快；当 $N$ 充分大时，$(1+t)^{4-2N}$ 的衰减等价于“几乎指数衰减”（almost exponential decay）。进一步，由关系式 $b = e_3 + \partial_3 \eta$（见 (1.6)）及上述估计，可得磁场扰动 $b - e_3 = \partial_3 \eta$ 满足

即 $b(t) \to e_3$ 在 $H^{2(N+1)}$ 范数下以代数速率收敛，且当 $N$ 大时亦呈几乎指数形式。值得注意的是，该结论要求背景磁场严格垂直于边界（$b_0 = e_3$），但作者在 Remark 2.2 中指出，结果可能对小角度倾斜的情形仍成立，尽管技术上更复杂。

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

本文的证明框架建立在 Lagrangian 坐标系与两层能量法（two-tier energy method）之上，其核心创新点在于通过 Helmholtz 投影解耦压力并重构角速度的正则性估计，从而克服双重退化带来的导数损失与非耗散性障碍。

第一，Lagrangian 坐标系的引入（见 (1.3)-(1.4)）将非线性对流项 $\tilde{u} \cdot \nabla \tilde{u}$ 吸收到时间导数中，消除了欧拉坐标下能量估计中难以处理的 $\tilde{u} \cdot \nabla$ 项。然而，这引入了变形梯度 $A = [(I_3 + \nabla \eta)^{-1}]^T$，导致 Laplacian 项变为 $\Delta_A u = \operatorname{div}_A \nabla_A u$，其展开包含 $\nabla^2 \eta$ 的高阶项。由于 $\eta$ 满足无粘性方程 $\partial_t \eta = u$，其本身不具耗散，故 $\|\nabla^2 \eta\|$ 无法被直接控制。本文的关键洞察（Remark 2.4）在于：稳态背景 $e_3$ 所诱导的项 $\partial_3^2 \eta$（见 (1.7b)）天然提供了对 $\eta$ 的垂直二阶耗散。该结构仅在磁场非零时出现，是稳定性的物理来源。

第二，Helmholtz 投影 $P$ 是处理压力与角速度耦合的核心工具。作者并未直接对 (1.7b) 应用 Stokes 算子（因其含未知 $w$ 项），而是先作用 $P$ 消去压力梯度 $\nabla p$，得到

此式将问题转化为对 $P(\nabla_A \times w)$ 的估计。但 $P$ 是正交投影，$\|\nabla_A \times w\|_{H^k}$ 无法由 $\|P(\nabla_A \times w)\|_{H^k}$ 直接控制，因余项 $\nabla \phi$（满足 $\Delta \phi = 0$, $\partial_3 \phi|_{x_3=d} = (\partial_1 w_2 - \partial_2 w_1)|_{x_3=d}$）的水平梯度 $\nabla_h \phi$ 需要单独估计。为此，作者发展了引理 3.2 所述的各向异性 Helmholtz 投影估计：对 $W = \nabla \times w = PW + \nabla \phi$，有

其中 $\|\cdot\|_{k,\lambda}$ 为各向异性范数（见 (2.1)）。该估计的证明（第 9–10 页）利用了条形域上 Poisson 方程 $\Delta \phi = 0$ 的 Fourier 分析：对低频 $|\xi| \leq 1$，$\|\nabla_h \phi\|$ 被边界值 $f = \partial_1 w_2 - \partial_2 w_1$ 控制，并通过迹不等式 $\|w(x_h, d)\|_{L^2_{x_h}} \lesssim \|w\|_0 + \|w\|_0^{1/2} \|\partial_3 w\|_0^{1/2}$ 化为 $\|w\|_0 \|w\|_1$；对高频 $|\xi| \gg 1$，则利用 $\partial_3 \hat{\phi}$ 与 $|\xi| \hat{\phi}$ 的精确关系（见 (3.15)）及 $\hat{f}$ 的边界积分表示，导出 $\||\xi| \hat{\phi}\|^2 \lesssim \|\partial_3 \hat{\phi}\|^2 + \sigma \|\hat{W}\|^2$。此估计是本文最关键技术突破，它允许作者将 $\|w\|_{H^m}$ 的控制最终归结为 $\|P(\nabla \times w)\|_{H^m}$ 与 $\|\partial_3 \phi\|_{H^m}$ 的控制。

第三，$\partial_3 \phi$ 的估计源自角动量方程 (1.7c) 的第三分量。作者通过对 (1.7c) 取 $\nabla \times$ 并分析其第三分量，提取出 $\partial_3 \phi$ 的演化方程，从而将其 $H^m$ 范数由 $\|P(\nabla \times w)\|$ 和 $\|u\|$ 控制（见 Lemma 6.2）。最终，整个能量闭合依赖于对 $P \partial_3^2 \eta$ 的估计（Lemma 6.1），该估计利用了 (1.7b) 中 $\partial_3^2 \eta$ 项与 $P \Delta_A u$ 的耗散配对，以及前述各向异性估计对 $w$ 的闭环控制。

4. 与之前工作的比较

本文与既有工作的比较可从三个维度展开：

1. 与全耗散磁微极系统比较：Rojas-Medar-Boldrini [27] 和 Xiang-Yang [35] 对 $\nu, \kappa_i > 0$ 的三维系统建立了全局适定性，其证明依赖于 $\nu \Delta \tilde{b}$ 和 $\kappa_i \Delta \tilde{w}$ 提供的完整耗散。本文则完全移除了这两项，将问题提升至双曲-椭圆耦合的临界情形，技术上需放弃标准能量法，转而依赖背景磁场诱导的结构耗散与 Lagrangian 框架。

2. 与非电阻 MHD 系统比较：Lin-Xu-Zhang [20] 和 Xu-Zhang [36] 对 $\nu = 0, \zeta = 0, \tilde{w} \equiv 0$ 的三维非电阻 MHD，在类似条形域中证明了围绕 $e_3$ 的稳定性，其核心是发现 $\partial_3^2 \eta$ 项与 $\Delta u$ 的配对耗散。本文的推进在于：增加了角速度 $\tilde{w}$ 及其与 $\tilde{u}$ 的 $\nabla \times$ 耦合，这引入了新的导数损失障碍。作者通过 Helmholtz 投影与各向异性估计（引理 3.2）解决了此问题，而此前 MHD 工作中无需处理 $\tilde{w}$。

3. 与压缩性对应模型比较：作者明确引用了文献 [11]（压缩性磁微极系统），指出其使用 Lagrangian 坐标与两层能量法。但关键差异在于：压缩性模型中压力由状态方程确定，而不可压模型的压力 $p$ 由椭圆方程 $\Delta p = \operatorname{div} \operatorname{div} (\tilde{u} \otimes \tilde{u} - \tilde{b} \otimes \tilde{b} - \zeta \tilde{w} \times)$ 决定，显式耦合 $\tilde{u}$ 与 $\tilde{w}$。因此，[11] 中用于估计 $\tilde{w}$ 的 div-curl 方法在此失效。本文的代价是引入了更复杂的各向异性 Helmholtz 分析（引理 3.2），其证明长度与技术深度远超 [11]。此外，[11] 未涉及条形域的边界效应，而本文的 $\partial_3 \phi$ 边界条件（见 (3.9)）是各向异性估计的关键输入，这在全空间或周期域中不存在。

综上，本文的推进是本质性的：在保持三维、不可压、条形域物理设定的前提下，首次将稳定性结果从全耗散、非电阻 MHD 推广至双重退化的磁微极系统。其代价是技术框架的高度特化——Lagrangian 坐标、Helmholtz 投影、各向异性 Fourier 分析——这些工具目前尚难推广至一般有界域或非垂直背景磁场，正如 Remark 2.2 所承认的局限性。

5. Navier–Stokes 方程在 $\dot{B}_{\infty,σ}^{-1}(\mathbb{R}^3)$ 中具有大垂直速度的全局解

Global solutions to the Navier–Stokes equations with large vertical velocities in $\dot{B}_{\infty,σ}^{-1}(\mathbb{R}^3)$

作者: Mikihiro Fujii

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04918v1

摘要翻译: 本文研究三维不可压 Navier–Stokes 方程的 Cauchy 问题，证明了在如下框架下存在唯一全局解：速度场的水平分量在某些临界 Besov 空间（包括经典的 Fujita–Kato 类）中充分小，而垂直分量在较宽泛的函数空间 $\dot{B}_{\infty,σ}^{-1}(\mathbb{R}^3)$（其中 $1 \leq σ< \infty$）中可以任意大；已知 Navier–Stokes 方程在此空间中是不适定的。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压Navier-Stokes方程的Cauchy问题：

其中 $u = (u_1,u_2,u_3): [0,\infty) \times \mathbb{R}^3 \to \mathbb{R}^3$ 为未知速度场，$P: (0,\infty) \times \mathbb{R}^3 \to \mathbb{R}$ 为压力，$a = (a_1,a_2,a_3)$ 为满足 $\mathrm{div}\, a = 0$ 的给定初值。该方程具有尺度不变性：若 $(u,P)$ 是解，则对任意 $\lambda > 0$，$u_\lambda(t,x) := \lambda u(\lambda^2 t, \lambda x)$、$P_\lambda(t,x) := \lambda^2 P(\lambda^2 t, \lambda x)$ 亦为解。因此，一个Banach空间 $X \subset \mathcal{S}'(\mathbb{R}^3)$ 被称为尺度临界空间，当且仅当 $\|u_\lambda(0,\cdot)\|_X = \|u(0,\cdot)\|_X$ 对所有 $\lambda > 0$ 成立。经典例子包括 $\dot{H}^{1/2}(\mathbb{R}^3)$、$L^3(\mathbb{R}^3)$ 及 Besov 空间 $\dot{B}_{p,q}^{3/p-1}(\mathbb{R}^3)$（$1 \leq p < \infty$, $1 \leq q \leq \infty$），以及端点空间 $\dot{B}_{\infty,q}^{-1}(\mathbb{R}^3)$ 和 $\mathrm{BMO}^{-1}(\mathbb{R}^3)$。

Fujita–Kato 理论确立了在临界空间中构造小初值全局解与大初值局部解的基本范式。Koch–Tataru 证明了在 $\mathrm{BMO}^{-1}$ 中小初值的全局适定性，这是目前已知最大的保证适定性的临界空间。然而，关键障碍在于端点空间 $\dot{B}_{\infty,q}^{-1}(\mathbb{R}^3)$ —— 已有工作 [2,16,17] 严格证明：对任意 $1 \leq q \leq \infty$，Navier-Stokes 方程在 $\dot{B}_{\infty,q}^{-1}(\mathbb{R}^3)$ 中是不适定的（ill-posed），即解映射 $a \mapsto u$ 在该空间中不连续；更甚者，存在初值序列 $a_n \to 0$ 于 $\dot{B}_{\infty,q}^{-1}$，但对应解在任意小时间区间上无界。这一现象源于双线性项 $\int_0^t e^{(t-\tau)\Delta} \mathbb{P}(u \cdot \nabla v)(\tau)\,d\tau$ 在 $E_{\infty,q}(0,T)$ 范数下缺乏一致有界性，即不存在常数 $C$ 使得

因此，在 $\dot{B}_{\infty,q}^{-1}$ 中构造非平凡解（尤其全局解）被视为根本性困难。现有突破均需额外结构：Chemin–Gallagher–Paicu [4] 利用Gevrey正则性假设，在 $\dot{B}_{\infty,\infty}^{-1}$ 中构造全局解，但其初值需基于某 $L^2$-Besov空间上的小Gevrey向量场；Iwabuchi–Nakamura [9] 则要求 $\|a\|_{\dot{B}_{\infty,\infty}^{-1}}$ 小于 $\log^{-1}(e + \|a\|_{\dot{B}_{p,\infty}^{3/p-1}})$（$3 < p < \infty$），本质仍是端点范数的小性控制。

本文旨在突破这一范式：不强加Gevrey条件，不依赖端点范数的小性，而允许垂直分量 $a_3$ 在整个 $\dot{B}_{\infty,\sigma}^{-1}(\mathbb{R}^3)$（$1 \leq \sigma < \infty$）中任意大，同时仅要求水平分量 $a_h = (a_1,a_2)$ 在较弱临界空间 $\dot{B}_{p,q}^{3/p-1}(\mathbb{R}^3)$（$3/2 < p < 3$）中足够小。这触及了“各向异性容错”的极限——因 $\dot{B}_{\infty,\sigma}^{-1}$ 中函数可极度振荡且非 $L^p$ 可积（如 $a_3(x) = A \cos(k_1 x_1 + k_2 x_2) \in \dot{B}_{\infty,\sigma}^{-1} \setminus \dot{B}_{p,q}^{3/p-1}$），传统能量或傅里叶方法完全失效。核心科学空白在于：如何在端点空间不适定性未被消除的前提下，利用不可压约束 $\partial_{x_3} u_3 = -\nabla_h \cdot u_h$ 所诱导的水平-垂直耦合结构，将非线性效应“降阶”并实现可控估计。

2. 主要定理或主要结果

定理 1.1（全局适定性）：设参数 $p,q,r,\theta,\sigma$ 满足

则存在依赖于这些参数的正常数 $\eta = \eta(p,q,r,\theta,\sigma)$ 与 $C = C(p,q,r,\theta,\sigma)$，使得对任意满足 $\mathrm{div}\, a = 0$ 的初值 $a = (a_h,a_3)$，若

则方程 (1.1) 存在唯一全局解 $u = (u_h,u_3)$，其中

且满足先验估计

此处 $E_{p,q}(0,\infty)$ 为 Chemin–Lerner 空间，其范数定义为 $\|u\|_{E_{p,q}(T_1,T_2)} := \|u\|_{\widetilde{L}^\infty(T_1,T_2; \dot{B}_{p,q}^{3/p-1})} + \|u\|_{\widetilde{L}^1(T_1,T_2; \dot{B}_{p,q}^{3/p+1})}$；类似定义 $E_{\infty,\sigma}$。值得注意的是，(1.3) 中的指数型小性条件是本质的：当 $\|a_3\|_{\dot{B}_{\infty,\sigma}^{-1}}$ 增大时，$\|a_h\|_{\dot{B}_{p,q}^{3/p-1}}$ 需以 $\exp(-C \|a_3\|^{r/\theta})$ 速率衰减。参数范围 $3/2 < p < 3$ 是技术关键，源于引理 2.3 中对 $u_3 \partial_{x_3} u_h$ 项的估计所需 $3/p - 2 + 2/\rho' < 0$，即 $\rho < 2p/3$；而 $r > 2$ 保证了时间积分的 Hölder 共轭 $r'$ 有限。此外，$\sigma < \infty$ 是必要限制（见 Remark 1.2(1)），因 $\sigma = \infty$ 时 $\dot{B}_{\infty,\infty}^{-1}$ 的不适定性更强，现有框架无法处理。

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

本文的核心障碍在于：标准双线性估计在 $\dot{B}_{\infty,\sigma}^{-1}$ 中彻底失效，而 $a_3$ 的任意大小排除了任何端点范数小性假设。作者的真正创新点在于三重技术组合：

1. 各向异性方程重构与线性化策略：利用不可压条件 $\partial_{x_3} u_3 = -\nabla_h \cdot u_h$，将原方程 (1.1) 重写为

   $$ \partial_t u_h - \Delta u_h + (u_h \cdot \nabla_h) u_h + u_3 \partial_{x_3} u_h + \nabla_h P = 0,\tag{1.4} $$
   $$ \partial_t u_3 - \Delta u_3 + (u_h \cdot \nabla_h) u_3 - u_3 \nabla_h \cdot u_h + \partial_{x_3} P = 0.\tag{1.5} $$

   关键观察是：(1.5) 中唯一含 $u_3$ 二次项的非线性项为 $-u_3 \nabla_h \cdot u_h$，而 $\nabla_h \cdot u_h = -\partial_{x_3} u_3$，故该项实为 $u_3 \partial_{x_3} u_3$ —— 但此形式在 (1.5) 中并未显式出现；相反，通过代入 $\nabla_h \cdot u_h = -\partial_{x_3} u_3$，(1.5) 中实际出现的 $u_3$-二次项仅为 $-u_3 \nabla_h \cdot u_h = u_3 \partial_{x_3} u_3$，然而作者刻意保留 $-u_3 \nabla_h \cdot u_h$ 形式，并将其视为关于 $u_h$ 的线性项（因 $u_h$ 是主控变量）。这使得 $u_3$ 的演化方程在形式上近似线性，从而可期望获得形如

   $$ \left\|\int_0^t e^{(t-\tau)\Delta} (u_3 \nabla_h \cdot u_h)(\tau)\,d\tau\right\|_{E_{\infty,\sigma}(0,T)} \leq C \|u_h\|_{E_{p,q}(0,T)} \|u_3\|_{E_{\infty,\sigma}(0,T)}\tag{1.6} $$

   的估计 —— 这正是引理 2.2 第三项所证，且其成立依赖于 $3/2 < p < 3$（确保 $\dot{B}_{p,q}^{3/p-1}$ 嵌入 $L^\infty$ 的弱形式）。

2. 时间区间分解法（Time Decomposition Lemma 2.1）：为克服 $\|u_3\|_{E_{\infty,\sigma}}$ 可能很大的问题，作者不采用传统的加权能量法（如 [13,18,19] 中的 $t^\alpha$ 权），而是引入关键引理 2.1：对任意 $f \in \widetilde{L}^r(0,T; \dot{B}_{p,\sigma}^s)$ 和 $\varepsilon > 0$，可将 $[0,T]$ 分解为 $N+1$ 个子区间 $[T_n,T_{n+1}]$，使得 $\|f\|_{\widetilde{L}^r(T_n,T_{n+1}; \dot{B}_{p,\sigma}^s)} \leq \varepsilon$。此处 $N \approx \|f\|_{\widetilde{L}^r}^r / \varepsilon^r$。在迭代构造中，作者在每个子区间上应用 Banach 不动点定理，要求 $\|u_3\|_{\widetilde{L}^r(T_n,T_{n+1}; \dot{B}_{\infty,\sigma}^{-1+2/r})}$ 足够小（由引理 2.1 保证），从而控制非线性项 $u_3 \nabla_h \cdot u_h$。该方法将全局存在性归结为证明 $N$ 有限，而这由 (1.3) 中的指数型小性条件保障：因 $\|u_h\|_{E_{p,q}} \lesssim \eta$，且 $\|u_3\|_{E_{\infty,\sigma}}$ 的增长被 $\eta^2$ 控制，故总“坏区间”数 $N$ 有界。

3. 精细化的仿积（para-product）分析与混合范数估计：为验证 (1.6) 及其他关键估计，作者发展了一套针对混合 Besov 空间的仿积工具（引理 2.2–2.4）。例如，在估计 $\int_0^t e^{(t-\tau)\Delta} (u_h \cdot \nabla_h) u_3\,d\tau$ 时，将乘积分解为 $T_{u_h} \nabla_h u_3 + R(u_h, \nabla_h u_3) + T_{\nabla_h u_3} u_h$，并分别利用引理 A.1（热半群映射性质）、A.2（仿积估计）及 Besov 空间嵌入 $\dot{B}_{p,q}^{3/p-1+2/r} \hookrightarrow \dot{B}_{\infty,\infty}^{-1+2/r}$（因 $p < 3$）来控制各项。特别地，对最坏项 $u_3 \partial_{x_3} u_h$（引理 2.3），作者采用双重 Hölder 指数 $\rho$ 与 $r$：因 $3/p - 2 + 2/\rho' < 0$，$\dot{B}_{p,q}^{3/p-2+2/\rho'}$ 可控于 $L^\infty$，从而 $\|T_{\partial_{x_3} u_h} u_3\|_{\widetilde{L}^1 \dot{B}_{p,q}^{3/p-1}} \lesssim \|u_3\|_{\widetilde{L}^\rho \dot{B}_{\infty,\sigma}^{-1+2/\rho}} \|u_h\|_{\widetilde{L}^{\rho'} \dot{B}_{p,q}^{3/p-1+2/\rho'}}$；其余项则用 $r$-范数控制。这种双重尺度处理是突破 $p=2$ 限制的关键。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开：

1. 与端点空间不适定性工作 [2,16,17] 的对比：这些工作构造了初值 $a_{Q,r}(x) = Q r^{-1/2} \sum_{k=1}^r (0, \cos(h_k x_1), \cos(h_k x_1 - x_2))^\top$，其第二分量 $a_2$ 的高频振荡驱动了解的爆炸。本文的初值示例 $a(x) = (0,0,A \cos(k_1 x_1 + k_2 x_2)) + \delta \psi$ 却将大振荡完全置于 $a_3$，而 $a_h = \delta \psi_h$ 保持小且光滑。Remark 1.2(4) 明确指出：$a_{Q,r}$ 的 $a_2$ 是不适定性的“本质角色”，而本文方案规避了该机制，表明不适定性并非均匀分布于所有分量，而是对水平分量更敏感。代价是：本文无法处理 $\sigma = \infty$（因 [2,16,17] 的反例恰在 $\dot{B}_{\infty,\infty}^{-1}$），且必须牺牲 $a_h$ 的任意性。

2. 与 Chemin–Gallagher–Paicu [4] 的对比：[4] 在 $\dot{B}_{\infty,\infty}^{-1}$ 中构造全局解，但要求初值形如 $a = a^{\mathrm{Gevrey}} + a^{\mathrm{small}}$，其中 $a^{\mathrm{Gevrey}}$ 属于某 $L^2$-Besov 空间上的 Gevrey 类（即傅里叶系数衰减如 $e^{-c|\xi|^{1/s}}$）。本文完全摒弃 Gevrey 假设，仅需 $a_h \in \dot{B}_{p,q}^{3/p-1}$（$p<3$），适用范围更广（例如包含各向异性振荡 $\cos(k_1 x_1 + k_2 x_2)$）。但代价是：[4] 允许 $\sigma = \infty$，而本文限定 $\sigma < \infty$；且 [4] 的 $a_h$ 可大，本文则必须小。

3. 与 Iwabuchi–Nakamura [9] 的对比：[9] 要求 $\|a\|_{\dot{B}_{\infty,\infty}^{-1}} \leq C_p / \log(e + \|a\|_{\dot{B}_{p,\infty}^{3/p-1}})$，本质是 $\dot{B}_{\infty,\infty}^{-1}$ 小性与 $\dot{B}_{p,\infty}^{3/p-1}$ 大性之间的对数补偿。本文则反转逻辑：$a_3$ 可任意大于 $\dot{B}_{\infty,\sigma}^{-1}$，而 $a_h$ 的小性需以指数速率补偿 $\|a_3\|$（见 (1.3)）。这使本文能处理 $a_3 \notin \dot{B}_{p,q}^{3/p-1}$ 的情形（如示例），而 [9] 要求 $a \in \dot{B}_{p,\infty}^{3/p-1}$。代价是：[9] 适用于 $p>3$，本文限于 $p<3$；且 [9] 允许 $\sigma = \infty$，本文不允许。

4. 与各向异性 Navier-Stokes 工作 [13,18,19] 的对比：这些工作针对水平粘性方程（即 $\Delta_h u$ 代替 $\Delta u$），其方法依赖于时间加权能量估计，但仅适用于 $p=2, r=4$ 且 $u_3$ 与 $u_h$ 使用相同空间。本文通过时间分解法与混合范数估计，将框架推广至 $3/2 < p < 3$ 及任意 $r>2$，并允许 $u_3$ 使用 $\dot{B}_{\infty,\sigma}^{-1}$（比 $\dot{B}_{2,1}^{0}$ 更粗糙）。代价是：本文的解映射仅在 $a_h \times a_3$ 空间上 Lipschitz，而 [13,18,19] 的加权法可能给出更强连续性。

6. 具接触间断的定常Euler流在具有一般上游数据的无穷长喷管中的存在性

Steady Euler flows with contact discontinuities in infinitely long nozzles with general upstream data

作者: Jun Chen, Xuemei Deng, Xiaoguang You

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04804v1

摘要翻译: 我们研究二维无穷长喷管中的定常可压Euler流，其中在无穷远处的分片光滑上游数据允许存在特征间断。除亚音速性条件外，未对数据施加任何额外约束。我们建立了与一条光滑接触间断曲线相关联的亚音速弱解的存在性与唯一性。原问题被重述为一个具不连续系数的散度型椭圆方程，使得接触间断条件天然地保留在该椭圆问题的解中。我们进一步研究下游渐近行为，并证明流场的收敛速率与喷管壁的收敛速率一致。

Comment: 无

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

1. 研究对象和背景

本文研究二维无穷长喷管 $\mathcal{N}_\zeta = \{(x_1,x_2) \in \mathbb{R}^2 \mid -1 + \zeta W_-(x_1) < x_2 < 1 + \zeta W_+(x_1)\}$ 中的定常可压欧拉流，其中 $W_\pm \in C^{2,\alpha}(\mathbb{R})$ 满足 $\|W_\pm\|_{2,\alpha;\mathbb{R}} \leq C$ 且 $W_\pm(x_1) \to 0$ 当 $x_1 \to -\infty$。喷管上下壁 $\Gamma_\zeta^\pm = \{(x_1,x_2) \mid x_2 = h_\zeta^\pm(x_1) = \pm 1 + \zeta W_\pm(x_1)\}$ 上施加滑移边界条件 $\mathbf{u} \cdot \mathbf{n}_\pm = 0$。所考察的控制方程为二维定常可压欧拉系统：

其中 $\mathbf{u} = (u_1,u_2)$ 为速度场，$\rho > 0$ 为密度，$p > 0$ 为压力，$E = |\mathbf{u}|^2/2 + p/((\gamma-1)\rho)$ 为单位质量总焓，$\gamma > 1$ 为绝热指数。对光滑解，伯努利函数 $B = |\mathbf{u}|^2/2 + \gamma p / ((\gamma-1)\rho)$ 与熵函数 $S = \gamma p / ((\gamma-1)\rho^\gamma)$ 沿流线守恒：$\operatorname{div}_x(\rho \mathbf{u} B) = 0$，$\operatorname{div}_x(\rho \mathbf{u} S) = 0$。声速 $c = \sqrt{\gamma p / \rho}$，马赫数 $M = |\mathbf{u}|/c$；本文限定于亚音速流，即 $\sup_{\mathcal{N}_\zeta^\pm} M < 1$。

关键建模特征在于上游无穷远 $x_1 \to -\infty$ 处给定分片光滑初值 $U_\ell = (u_\ell, 0, p_0, \rho_\ell)$，其在 $x_2 = x_2^* \in (-1,1)$ 处具有特征型间断（即接触间断），满足 $u_\ell, \rho_\ell \in C^{1,\alpha}([-1,x_2^*]) \cap C^{1,\alpha}([x_2^*,1])$，且 $\inf_{[-1,x_2^*) \cup (x_2^*,1]} u_\ell > 0$，$\inf \rho_\ell > 0$，$\sup u_\ell / c_\ell < 1$，其中 $c_\ell = \sqrt{\gamma p_0 / \rho_\ell}$。该间断导致物理量在 $x_2 = x_2^*$ 处跳跃，但要求压力 $p_0$ 为全局正常数以保证动量与能量方程相容。由此定义下、上质量通量 $m_- = \int_{-1}^{x_2^*} \rho_\ell(s) u_\ell(s)\,ds > 0$，$m_+ = \int_{x_2^*}^1 \rho_\ell(s) u_\ell(s)\,ds > 0$。

未知自由界面 $\Gamma_\zeta^* = \{(x_1,x_2) \mid x_2 = \omega^*(x_1)\}$ 将喷管分为两子区域 $\mathcal{N}_\zeta^\pm$，并满足接触间断条件：

其中 $[f] = f|_{\mathcal{N}_\zeta^+} - f|_{\mathcal{N}_\zeta^-}$ 为跨 $\Gamma_\zeta^*$ 的跳跃。由于解在 $\Gamma_\zeta^*$ 上不连续，需寻求弱解：对任意 $v \in C_0^\infty(\mathcal{N}_\zeta)$，有 $\int_{\mathcal{N}_\zeta} \mathbf{F} \cdot \nabla v \,dx = 0$，其中 $\mathbf{F}$ 为对应守恒律的通量向量。

该问题的重要性在于：它代表了真实气流中普遍存在的接触间断（如激波后熵层分离、不同气团交汇）在非均匀几何约束下的基本适定性问题。传统方法常依赖上游数据的强结构假设（如单调性、凸性、小扰动），而工程实际中上游往往含复杂间断与大涡量。因此，核心空白在于：能否在仅保留物理本质条件（亚音速性、正质量通量）的前提下，建立具一般分片 $C^{1,\alpha}$ 上游数据的亚音速弱解的存在唯一性，并刻画接触间断曲线 $\omega^*$ 的最优正则性（超越 Lipschitz）及下游渐近行为。

2. 主要定理或主要结果

定理1.2（存在唯一性）：存在 $\zeta_+ \in \mathbb{R}_+ \cup \{+\infty\}$ 与 $\zeta_- \in \mathbb{R}_- \cup \{-\infty\}$，使得对任意 $\zeta \in (\zeta_-, \zeta_+)$，问题1.1存在唯一解 $(U, \omega^*)$，其中 $U \in C^{1,\alpha}(\mathcal{N}_\zeta^+) \cap C^{1,\alpha}(\mathcal{N}_\zeta^-)$，$\omega^* \in C^{2,\alpha}(\mathbb{R})$，且满足 $\sup_{\mathcal{N}_\zeta^\pm} M < 1$ 与 $\inf_{\mathcal{N}_\zeta^\pm} (\rho u_1) > 0$。若某端点有限（如 $\zeta_+ < +\infty$），则对任一序列 $\zeta_n \to \zeta_+$，必有 $\sup_{\mathcal{N}_{\zeta_n}^\pm} M(\zeta_n) \to 1$ 或 $\inf_{\mathcal{N}_{\zeta_n}^\pm} (\rho(\zeta_n) u_1(\zeta_n)) \to 0$。

定理1.3（正则性提升）：设整数 $k \geq 2$，若 $U_\ell \in C^{k,\alpha}([-1,x_2^*]) \cap C^{k,\alpha}([x_2^*,1])$ 且 $W_\pm \in C^{k+1,\alpha}(\mathbb{R})$，则对任意 $\zeta \in (\zeta_-, \zeta_+)$，解 $U(\zeta)$ 在 $\mathcal{N}_\zeta^\pm$ 上属于 $C^{k,\alpha}$ 类，且 $\omega^* \in C^{k+1,\alpha}(\mathbb{R})$。

下游渐近性：若喷管壁 $W_\pm(x_1)$ 在 $x_1 \to +\infty$ 时存在极限，则解 $U$ 存在下游渐近状态；更精确地，若 $W_\pm(x_1) = O((1+x_1)^{-\beta})$（$\beta > 0$），则通过构造奇异屏障函数并应用极值原理，证明扰动 $\varphi - \varphi_r$（其中 $\varphi$ 为逆流函数，$\varphi_r$ 为参考流）满足 $|\varphi(x_1,x_2) - \varphi_r(x_1,x_2)| = O((1+x_1)^{-\beta})$，即下游收敛速率与壁面衰减速率完全匹配。

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

核心障碍在于：（1）接触间断导致系数在自由界面处不连续，传统椭圆理论失效；（2）上游数据无单调/凸性约束，无法用经典比较原理控制解的梯度；（3）自由边界与间断耦合，需同时处理几何非线性与间断传播。

作者的根本创新在于将 Euler–Lagrange 变换与逆流函数 $\varphi(y_1,y_2) = x_2$ 结合，将原自由边界问题转化为固定带域 $\mathcal{N} = \mathbb{R} \times (-m_-, m_+)$ 上的拟线性椭圆方程。变换 $y = T(x)$ 定义为 $y_1 = x_1$, $y_2 = \int_{h_\zeta^-(x_1)}^{x_2} \rho u_1(x_1,s)\,ds - m_-$，将 $\Gamma_\zeta^*$ 映为固定直线 $\Gamma^* = \mathbb{R} \times \{0\}$，并将欧拉系统（1.1）化为 Lagrangian 坐标下的守恒律（2.2）。关键推导是：由 $\varphi_{y_2} = 1/(\rho u_1) > 0$ 与 $\varphi_{y_1} = u_2/u_1$，结合伯努利 $B$ 与熵 $S$ 仅依赖于 $y_2$ 的事实，导出密度 $\rho$ 满足隐式方程 $X(\nabla \varphi) = Q(S,B,\rho)$，其中 $X(\xi) = \xi_2/(1 + \frac{1}{2}\xi_2^2)$，$Q(\eta,\rho) = \eta_2 \rho^2 - \eta_1 \rho^{\gamma+1}$。此方程在亚音速区域 $\mathcal{D} = \{(S,B,\nabla \varphi) \mid X(\nabla \varphi) < Q_{\max}(S,B)\}$ 内定义光滑函数 $\rho = \rho(S,B,\nabla \varphi)$，从而所有流变量均可由 $\varphi$ 表示：$u_1 = 1/(\rho \varphi_{y_2})$, $u_2 = \varphi_{y_1}/(\rho \varphi_{y_2})$, $p = \frac{\gamma-1}{\gamma} S \rho^\gamma$。

由此得到 $\varphi$ 满足的散度型椭圆方程：

其中系数 $a_{ij} = a_{ij}(S,B,\nabla \varphi)$ 在 $\Gamma^*$ 上不连续，但满足一致椭圆性：存在 $\lambda > 0$ 使得 $\sum_{i,j} a_{ij} \xi_i \xi_j \geq \lambda |\xi|^2$ 对所有 $\xi \in \mathbb{R}^2$ 成立。作者利用 De Giorgi–Nash 理论处理间断系数，证明 $\varphi \in C^{1,\alpha}(\mathcal{N}^\pm) \cap C(\mathcal{N})$；再通过 Schauder 估计与系数的 $C^{1,\alpha}$ 依赖性，建立 $\varphi \in C^{2,\alpha}(\mathcal{N}^\pm)$，从而 $\omega^* \in C^{2,\alpha}(\mathbb{R})$——此正则性高于文献[9]的 Lipschitz 性。

为处理非结构上游数据，作者发现：$\varphi$ 的高阶 Hölder 范数本质上由亚音速性 $\sup M < 1$ 和无滞止性 $\inf (\rho u_1) > 0$ 控制，而非依赖于 $U_\ell$ 的单调性。这一观察能使先验估计独立于上游的精细结构，仅需 $C^{1,\alpha}$ 分片光滑性。最后，为获得尖锐下游衰减，作者构造奇异屏障函数 $\Psi(y_1,y_2) = C (1+y_1)^{-\beta} \chi(y_2)$（$\chi$ 为截断函数），利用极值原理比较 $\varphi - \varphi_r$ 与 $\Psi$，严格证明 $|\varphi - \varphi_r| \leq C (1+y_1)^{-\beta}$。

4. 与之前工作的比较

与 Chen–Feldman [7,8] 的位势流结果相比，本文处理全欧拉系统，允许非零涡量与熵变，且无需小 Bernoulli 数假设。与 Xie–Xin [17] 相比，本文取消了上游 Bernoulli 函数“足够小”的限制；与 Chen–Deng–Xiang [6] 相比，本文不要求熵与 Bernoulli 函数的 $C^{1,1}$ 振幅“足够小”，亦不预设质量通量在特定区间内。Du–Xie–Xin [13] 依赖上游水平速度 $u_\ell$ 的二阶导数符号（凸性条件），而本文彻底移除该结构性假设；Chen [4] 与 Duan–Luo [14] 虽推广至非等熵与轴对称情形，但仍继承凸性要求。最接近的工作是 Chen et al. [9]，其允许大涡量与间断，但强制上游数据满足单调性及壁面单调性，且仅获 $\omega^*$ 的 Lipschitz 正则性；本文则证明 $\omega^* \in C^{2,\alpha}$，显著提升正则性。

在方法论上，[9] 使用 Euler–Lagrange 变换但未显式构建 $\rho$ 的隐式方程，其系数正则性分析较粗；本文通过 $X(\nabla \varphi) = Q(S,B,\rho)$ 的精细分析，确立系数对 $\nabla \varphi$ 的 $C^{1,\alpha}$ 依赖性，为 Schauder 估计提供坚实基础。下游渐近方面，[6][9] 仅给出存在性，未量化收敛速率；本文首次证明衰减速率与壁面完全匹配，且通过奇异屏障函数实现尖锐估计，代价是需 $W_\pm \in C^{2,\alpha}$ 以支撑 $\beta$ 阶衰减。局限在于：仍要求上游 $U_\ell$ 分片 $C^{1,\alpha}$，尚未处理 $L^\infty$ 初值或激波型间断；且最大存在区间 $(\zeta_-,\zeta_+)$ 的有限性判据（1.14）未给出显式阈值，仅定性描述临界态。

7. 可压缩流体-结构相互作用热力学的适定性与正则性

Well-posedness and regularity in thermodynamics of compressible fluid-structure interactions

作者: D. Breit, P. M. Ngougoue Ngougoue

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04752v1

摘要翻译: 我们考虑一类具有一般粘性、可压缩且导热的流体与位于流体区域边界的弹性壳体之间的相互作用。该流体由可压缩 Navier-Stokes-Fourier 方程组描述，而壳体则遵从一个粘弹性梁方程演化。二者通过运动学边界条件及力的平衡关系耦合。我们的第一个结果是：该非线性偏微分方程耦合系统在光滑函数空间中具有局部时间适定性。进一步，我们证明了一个条件正则性判据，其精神类似于不可压 Euler 方程的 Beale-Kato-Majda 条件，并辅以密度与温度的有界性；该判据建立在该耦合系统的弱-强唯一性结果之上。

Comment: 无

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

1. 研究对象和背景

本文研究三维空间中粘性、可压缩、热传导流体与位于其边界上的弹性壳体的耦合动力学系统。该系统由两个核心子系统构成：(1) 流体部分满足完整的可压缩 Navier-Stokes-Fourier 方程组，定义在随时间演化的区域 $\Omega_{\eta(t)} \subset \mathbb{R}^3$ 上；(2) 壳体部分由定义在二维参数域 $\omega \subset \mathbb{R}^2$ 上的粘弹性梁方程描述。二者通过运动学边界条件与力平衡条件强耦合。

流体状态变量为速度场 $v : I \times \Omega_{\eta} \to \mathbb{R}^3$、密度 $\rho : I \times \Omega_{\eta} \to \mathbb{R}^+$ 和绝对温度 $\vartheta : I \times \Omega_{\eta} \to \mathbb{R}^+$，满足如下方程组（即原文(1.2)）：

其中压力采用理想气体状态方程 $p(\rho, \vartheta) = \rho \vartheta$；比热容 $c_v = 1/(\gamma - 1) > 0$，$\gamma > 1$；热导率 $\kappa > 0$；粘性应力张量为

且满足物理约束 $\mu > 0$, $\lambda + \frac{2}{3}\mu \ge 0$。

壳体位移函数 $\eta : I \times \omega \to \mathbb{R}$ 满足粘弹性梁方程（即原文(1.1)）：

其中 $\varphi_\eta(y) = \varphi(y) + \eta(t,y) n(y)$ 是边界参数化映射，$n$ 与 $n_\eta$ 分别为参考边界 $\partial \Omega$ 与变形边界 $\partial \Omega_\eta$ 的单位外法向量，$\tau = S(\nabla v) - p(\rho, \vartheta) I_{3\times3}$ 为柯西应力张量。初始条件为 $\eta(0,\cdot) = \eta_0$, $\partial_t \eta(0,\cdot) = \eta^*$。

耦合机制体现于两个关键条件：(i) 运动学边界条件 $v \circ \varphi_\eta = (\partial_t \eta) n$，确保流体速度在边界处等于壳体法向速度；(ii) 热流边界条件 $\nabla \vartheta \cdot n_\eta \circ \varphi_\eta^{-1} = 0$，即绝热边界（无热通量穿过壳体）。

该模型在水弹性力学、航空航天工程及生物力学中具有根本重要性，其数学挑战源于三重非线性：几何非线性（区域 $\Omega_\eta$ 随 $\eta$ 动态变形）、物理非线性（Navier-Stokes-Fourier 中的对流项、非线性压力与粘性耗散）以及耦合非线性（边界条件将流体应力与壳体动力学显式关联）。已有工作存在显著空白：对于不可压流体，弱解存在性已在一般几何下建立（[16,27,28]）；对于等熵可压流体，弱解存在性亦有结果（[4]），且局部强解与条件正则性准则已在平板情形（flat plate）和二维壳体中取得（[21,29,25]）；但对于热传导可压流体-壳体系统，其局部强解的存在性仍是开放问题——此前仅见于平板情形（[23]），而对一般曲面壳体（elastic shell）尚未解决；更关键的是，尚无任何关于该系统的条件正则性准则（conditional regularity criterion）或弱-强唯一性（weak-strong uniqueness）结果。本文旨在填补这三项根本性空白：(1) 在光滑函数空间中建立局部强解的存在唯一性；(2) 给出一个融合 Beale-Kato-Majda 型条件与热力学量有界性的 blow-up 判据；(3) 构造适配热力学结构的相对熵泛函，证明弱-强唯一性，从而为条件正则性提供逻辑基础。

2. 主要定理或主要结果

本文在严格数学框架下建立了三个层次的核心结果：

1. 局部强解存在唯一性（Theorem 2.1）：设初始数据 $(\rho_0, v_0, \vartheta_0, \eta_0, \eta^*)$ 满足

且满足相容性条件 $\|\eta_0\|_{L^\infty(\omega)} < L$（$L$ 为参考边界的管状邻域半径）及运动学相容性 $v_0 \circ \varphi_{\eta_0} = \eta^* n$ on $\omega$，则存在 $T^* > 0$，使得系统 (1.1)–(1.4) 在 $I^* = (0, T^*]$ 上存在唯一强解 $(\rho, v, \vartheta, \eta)$，其正则性为：

2. 加速度估计与高阶先验估计（Theorems 2.2 & 2.3）：设强解满足以下四类假设： (A1) Serrin型动量控制：$\rho^{1/2} v \in L^s(I^*; L^r(\Omega_\eta))$，其中 $2/s + 3/r \le 1$, $r \in (3, \infty]$, $s \in [2, \infty)$； (A2) 可压性与热效应控制：(a) $\|\operatorname{div} v\|_{L^1(I^*; L^\infty(\Omega_\eta))} < \infty$，(b) $\|\vartheta\|_{L^\infty(I^*; L^\infty(\Omega_\eta))} < \infty$； (A3) 结构几何正则性：$\eta \in L^\infty(I^*; C^1(\omega))$； (A4) Beale-Kato-Majda型控制：$\int_{I^*} \left( \|\nabla \rho\|^2_{L^\infty(\Omega_\eta)} + \|\nabla v\|^2_{L^\infty(\Omega_\eta)} + \|\nabla \vartheta\|^2_{L^\infty(\Omega_\eta)} \right) dt < \infty$。

在此假设下，Theorem 2.2 给出加速度能量估计： $$ E_{\text{acc}} := \sup_{I^} \int_\omega \left( |\partial_t \nabla_y \eta|^2 + |\nabla_y \Delta_y \eta|^2 \right) dy + \frac{\mu}{2} \sup_{I^} \int_{\Omega_\eta} |\nabla v|^2 dx + \frac{\lambda + \mu}{2} \sup_{I^*} \int_{\Omega_\eta} |\operatorname{div} v|^2 dx \

• \int_{I^} \int_\omega \left( |\partial_t^2 \eta|^2 + |\partial_t \Delta_y \eta|^2 \right) dy dt + \int_{I^} \int_{\Omega_\eta} \left( |\nabla^2 v|^2 + \rho |\partial_t v|^2 + |\nabla p|^2 \right) dx dt \ \lesssim \int_{I^} \int_{\Omega_\eta} |\partial_t \vartheta|^2 dx dt + \int_\omega \left( |\nabla_y \eta^|^2 + |\nabla_y \Delta_y \eta_0|^2 \right) dy \
• \int_{\Omega_{\eta_0}} \left( |\nabla v_0|^2 + |\operatorname{div} v_0|^2 + \vartheta_0 \rho_0 \operatorname{div} v_0 \right) dx. $$

Theorem 2.3 在 (A1)–(A4) 及相容性条件 $(\partial_t v) \circ \varphi_\eta|_{t=0} = (\partial_t^2 \eta) n|_{t=0}$ 下，给出高阶先验估计 $E_{\text{high}} \lesssim \|\eta_0\|^2_{W^{5,2}(\omega)} + \|\eta^*\|^2_{W^{3,2}(\omega)} + \|v_0\|^2_{W^{3,2}(\Omega_{\eta_0})} + \|\rho_0\|^2_{W^{3,2}(\Omega_{\eta_0})} + \|\vartheta_0\|^2_{W^{3,2}(\Omega_{\eta_0})}$，其中 $E_{\text{high}}$ 包含 $\partial_t^2 \nabla_y \eta$, $\partial_t \nabla_y \Delta_y \eta$, $\partial_t \nabla v$, $\partial_t \operatorname{div} v$, $\partial_t \nabla^2 v$, $\partial_t^2 v$, $\partial_t^3 \eta$, $\partial_t^2 \Delta_y \eta$, $\partial_t \Delta_y^2 \eta$, $\Delta_y^3 \eta$, $\nabla^4 v$, $\nabla^4 \vartheta$, $\nabla^3 p$, $\nabla^3 v$, $\partial_t \nabla^2 \vartheta$, $\partial_t \nabla^2 \rho$ 等所有四阶空间导数与二阶时间导数的 $L^2$ 范数。

3. 弱-强唯一性与条件正则性（Theorem 2.4 & Corollary 2.1）：设 $(\rho_1, v_1, \vartheta_1, \eta_1)$ 为有限能量弱解，$(\rho_2, v_2, \vartheta_2, \eta_2)$ 为强解，二者初值相同，且 $\eta_1 \in L^\infty(I; C^1(\omega))$，则二者恒等。由此推出 Corollary 2.1：若弱解满足 (A1)–(A4) 及相容性条件，则其必为强解，且在 $L^\infty(I; C^1(\omega))$ 类变形中唯一。

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

本文的证明框架围绕三大技术障碍展开，其创新点在于对热传导可压流体-壳体耦合系统特有的非线性与几何复杂性的精细处理：

1. 固定域变换与Hanzawa映射的精细分析：由于区域 $\Omega_\eta$ 依赖于未知函数 $\eta$，标准PDE理论无法直接应用。作者采用Hanzawa变换 $\Psi_\eta$ 将动态区域映射回固定参考域 $\Omega$（公式(2.5)）。关键创新在于Proposition 2.1中对 $\Psi_\eta$ 及其逆映射的Sobolev范数估计：$\|\Psi_\eta\|_{W^{k,p}(\Omega)} \lesssim 1 + \|\eta\|_{W^{k,p}(\omega)}$，且差估计 $\|\Psi_\eta - \Psi_\zeta\|_{W^{k,p}(\Omega)} \lesssim \|\eta - \zeta\|_{W^{k,p}(\omega)}$。此估计依赖于参考边界 $\partial \Omega$ 的曲率与管状邻域半径 $L$，为后续线性化与不动点论证提供了严格的函数空间基础。特别地，该映射要求 $\|\eta\|_{W^{1,\infty}(\omega)} < L$ 以保证映射的双射性与正则性，这是区别于固定域问题的关键几何约束。

2. 内能方程的转化与最大正则性估计：在固定域上，内能方程（即(1.2)第三式）经变换后变为一个系数依赖于 $\eta$ 的抛物型方程，其系数仅具有限正则性（$\eta \in W^{5,2}(\omega)$ 导致系数在 $W^{2,2}$ 空间）。作者在Section 3.3与4.2中发展了一套针对此类低正则系数抛物方程的最大正则性理论，其核心是估计(4.25)：利用 $\eta$ 的临界正则性（来自加速度估计）来控制系数的 $L^\infty$ 模，并结合Agmon-Douglis-Nirenberg理论，得到 $\vartheta$ 的 $W^{2,2}_t W^{4,2}_x$ 正则性。这一分析比等熵情形（[25]）更为精细，因为热传导项 $\kappa \Delta \vartheta$ 与粘性耗散项 $S(\nabla v):\nabla v$ 的耦合引入了额外的非线性源项，必须通过(A4)中的 $\|\nabla \vartheta\|_{L^\infty}$ 控制。

3. 热力学相对熵泛函的构造与弱-强唯一性证明：为克服温度自由度带来的困难，作者摒弃了等熵情形中使用的压力势 $P(\rho) = a\rho^\gamma$，而引入了基于吉布斯关系的“弹道自由能”（ballistic free energy）：

其中 $e = c_v \vartheta_1$, $s = c_v \log \vartheta_1 - \log \rho_1 + s_0$。该泛函是热力学系统中自然的能量比较工具，其二阶展开定义了相对熵 $E_{\text{rel}}$（公式(2.19)）。证明的关键步骤是推导余项 $R$ 的估计（公式(2.20)–(2.21)），其中作者将温度相关项的处理完全参照Feireisl-Jin-Novotný [10]，而流体-结构耦合项则沿用[25]的策略，并最终通过Grönwall不等式完成唯一性论证。这一构造是本文最本质的创新，它首次将热力学框架完整嵌入流体-结构相互作用的相对熵方法中。

4. 与之前工作的比较

本文在多个维度上显著推进了流体-结构相互作用（FSI）理论，其贡献需与三类既有工作细致比较：

1. 与不可压FSI工作的比较：对于不可压流体（$\rho = \text{const}$），弱解存在性已建立（[16,27,28]），且Ladyzhenskaya-Prodi-Serrin型正则性准则已被证明（[3]）。然而，这些结果严重依赖于 $\operatorname{div} v = 0$ 所带来的椭圆性提升与压力作为拉格朗日乘子的特殊结构。本文处理的可压情形彻底失去了这一优势，密度演化方程 $\partial_t \rho + \operatorname{div}(\rho v) = 0$ 引入了强烈的对流-扩散耦合，且压力 $p = \rho \vartheta$ 成为真正的状态变量。因此，本文的局部强解存在性（Theorem 2.1）是首个适用于三维一般曲面壳体的可压热传导FSI结果，而此前[3]的Serrin准则无法直接推广，因其未考虑温度与密度的联合演化。

2. 与等熵可压FSI工作的比较：等熵情形（$p = a \rho^\gamma$）的弱解存在性见[4]，局部强解见[29]，条件正则性见[25]。本文的突破在于：(i) 物理模型升级：将压力律从 $p = a \rho^\gamma$ 升级为 $p = \rho \vartheta$，使系统成为封闭的热力学系统，内能方程成为独立方程；(ii) 正则性准则扩展：[25]的BKM型条件仅含 $\|\nabla \rho\|_{L^\infty}$ 与 $\|\nabla v\|_{L^\infty}$，而本文新增 $\|\nabla \vartheta\|_{L^\infty}$（Assumption A4），并证明其必要性——因 $\dot{p} + p \operatorname{div} v = \rho \dot{\vartheta}$（公式(2.15)），压力演化不再由连续性方程单独决定，必须控制温度的时间导数，而这又依赖于内能方程中的 $S(\nabla v):\nabla v$ 项，最终导出A4；(iii) 唯一性框架重构：[25]使用等熵压力势 $P(\rho)$，而本文构造的 $H_{\vartheta_2}(\rho, \vartheta)$ 是热力学一致的泛函，其成功依赖于吉布斯关系的精确实现，这是对[10]工作的非平凡推广。

3. 与热传导可压FSI的先前尝试比较：弱解存在性已在[5]中建立，但其压力律为更复杂的 $p \sim \vartheta \rho + \rho^\gamma + \vartheta^4$，且未涉及正则性问题。本文的局部强解（Theorem 2.1）是首个针对理想气体律 $p = \rho \vartheta$ 的结果，其代价是初始数据需处于更高正则性空间（$W^{3,2}$ for $\rho,v,\vartheta$ and $W^{5,2}$ for $\eta$），这反映了热传导项带来的额外正则性要求。此外，Corollary 2.1所给出的条件正则性是一个‘全系统’判据，它同时要求密度、速度、温度的梯度 $L^\infty$ 可积，这比Nash猜想（[12]中仅需 $\|\rho\|_{L^\infty}, \|v\|_{L^\infty}, \|\vartheta\|_{L^\infty}$ 有界）更强，但也更切合实际——因为FSI中边界运动会放大内部奇点。最后，本文明确指出，将正则性门槛降至更低Sobolev空间（如[6,38]对固定域所做的）仍是开放问题，这凸显了动态边界带来的根本性困难。

8. 一种针对带随机扰动数据的Oberbeck-Boussinesq系统变体的连续数据同化方法

A continuous data assimilation method for a variant of Oberbeck-Boussinesq system with randomly perturbed data

作者: Eduard Feireisl, Madalina Petcu

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04458v1

摘要翻译: 我们证明了连续数据同化方法在维数 $d=2,3$ 下对Oberbeck-Boussinesq系统的收敛性。我们的工作假设是参考解有界，而同步化解以弱意义满足方程。主要工具是针对随机问题的相对能量不等式。

Comment: 无

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

1. 研究对象和背景

本文研究的是Oberbeck-Boussinesq（OB）系统的变体——修正Oberbeck-Boussinesq（MOB）系统，其物理背景源于弱分层Navier-Stokes-Fourier系统的奇异极限[3]，广泛应用于气象学与地球物理流体力学中温度驱动对流建模。该系统描述速度场 $ u = u(t,x) \in \mathbb{R}^d $ 与温度扰动 $ \theta = \theta(t,x) \in \mathbb{R} $ 的耦合演化，其原始形式（见(1.1)）为：

定义在有界区域 $ Q \subset \mathbb{R}^d $（$ d = 2,3 $），边界条件为：

其中 $ G = G(x) $ 是重力势，满足 $ \int_Q G \,dx = 0 $ 与 $ \Delta_x G = 0 $（见(2.1)(2.2)）；$ \theta_B $ 是边界的温度偏差分布；非局部边界条件反映原系统缩放前Dirichlet条件的奇异极限结构[3,7]。为处理该非齐次边界，作者引入变换 $ \Theta = \theta + \alpha \fint_Q \theta \,dx - \vartheta_B $，其中 $ \vartheta_B $ 是 $ \theta_B $ 的调和延拓（$ \Delta_x \vartheta_B = 0 $ in $ Q $, $ \vartheta_B|_{\partial Q} = \theta_B $），从而将系统等价重构为具有齐次Dirichlet边界条件的MOB系统（见(2.4)–(2.6)）：

其中算子 $ M $ 定义为 $ M(\Theta) = \Theta - \frac{\alpha}{\alpha+1} \fint_Q \Theta \,dx $，是 $ L^2(Q) $ 上有界、自伴、正定且具连续逆的线性算子（见(2.6)(2.7)）。该MOB系统已被证明存在全局弱解：对任意初值 $ u_0 \in L^2(Q;\mathbb{R}^d) $（满足 $ \mathrm{div}_x u_0 = 0 $）、$ \Theta_0 \in C(Q) $，存在解 $ (u,\Theta) $ 满足

当 $ d=2 $ 时，若数据与区域足够正则，则存在唯一强解；而 $ d=3 $ 时，即使对确定性MOB系统，其强解的存在唯一性仍是未解决问题（类比Navier-Stokes千禧难题）。本文的研究对象正是该MOB系统的连续数据同化问题：给定一个“观测解” $ (u,\Theta) $（满足全局有界性假设 $ \sup_{t\in(T_-,T_+)} (\|u(t,\cdot)\|_{L^\infty} + \|\Theta(t,\cdot)\|_{L^\infty}) \leq E $，见(3.1)），目标是构造一个“同步解” $ (\tilde{u},\tilde{\Theta}) $，使其在预测时段 $ (T,T_+) $ 内以高精度逼近 $ (u,\Theta) $，仅依赖于其在采样时段 $ (0,T) $ 内的随机扰动观测值 $ (\bar{I}_\delta(u), \bar{I}_\delta(\Theta)) = (I_\delta(u), I_\delta(\Theta)) + (\xi_1, \xi_2) $。此处 $ I_\delta $ 是正交投影型插值算子（见(3.3)），误差 $ (\xi_1,\xi_2) $ 分解为确定性部分 $ (e_{D,u}, e_{D,\Theta}) $ 与随机部分——由两个独立柱状Wiener过程 $ W_u = \sum_k v_k \chi_{u,k} $、$ W_\Theta = \sum_k h_k \chi_{\Theta,k} $ 驱动，其扩散系数 $ v_k \in L^\infty(Q;\mathbb{R}^d), h_k \in L^\infty(Q) $ 满足 $ \sum_k (\|v_k\|_{L^\infty}^2 + \|h_k\|_{L^\infty}^2) \leq E_S $（见(3.6)(3.7)）。同步系统（见(3.8)）即为此带随机强迫与nudging项的MOB系统：

由于 $ d=3 $ 时同步系统的强解存在性未知，作者采用弱martingale解框架（Definition 3.1），其核心是满足动量与热方程的积分恒等式（(3.12)(3.13)）及动能/热能不等式（(3.14)(3.15)）。本文试图填补的关键空白在于：在物理最相关但数学最困难的三维情形下，为带有随机观测误差的OB类系统建立严格、定量的连续数据同化收敛理论；此前工作（如[2,4]）集中于确定性或白噪声强度极小的情形，且多限于二维或线性化模型，缺乏对随机扰动、非线性耦合与三维不适定性三重挑战的统一处理。

2. 主要定理或主要结果

定理4.1（连续数据同化）设 $ Q \subset \mathbb{R}^d $（$ d = 2,3 $）为 $ C^{2,\alpha} $ 类有界区域，$ (u,\Theta) $ 为MOB系统（2.4）–（2.5）在时间区间 $ (T_-, T_+) $ 上满足有界性假设（3.1）的解，$ (\tilde{u}, \tilde{\Theta}) $ 为同步系统（3.8）–（3.10）在 $ (0,T_+) $ 上按定义3.1给出的弱martingale解。则对任意 $ \gamma > 0 $，存在仅依赖于 $ \mu, \kappa $ 及插值算子 $ I_\delta $ 具体形式的阈值 $ \Lambda_0 > 0 $、$ \delta_0 > 0 $，使得对所有 $ \tau \in (T, T_+) $、任意 $ \Lambda \geq \Lambda_0 $、$ 0 < \delta \leq \delta_0 $，有

$$ \mathbb{E}\left[ | (u - \tilde{u})(\tau, \cdot) |{L^2(Q;\mathbb{R}^d)}^2 + | (\Theta - \tilde{\Theta})(\tau, \cdot) |{L^2(Q)}^2 \right] \leq C \exp(K(\tau - T_+)) \Bigg[ \exp(-\gamma T) \max\left{ E, |u_0|_{L^2(Q;\mathbb{R}^3)}, |\tilde{\Theta}0|{L^2(Q)} \right}

• \exp(K(\tau - T_+)) \Lambda^2 \int_0^T \int_Q \mathbb{E}\left[ |e_{D,u}|^2 + |e_{D,\Theta}|^2 \right] dx dt
• \exp(K(\tau - T_+)) \Lambda^2 \left( \frac{\sigma_u^2}{2} \sum_{k=1}^\infty \int_0^T \int_Q |v_k|^2 dx dt + \frac{\sigma_\Theta^2}{2} \sum_{k=1}^\infty \int_0^T \int_Q h_k M^{-1}(h_k) dx dt \right) \Bigg]. $$

其中常数 $ K $ 仅依赖于 $ E, \mu, \kappa $；符号 $ \lesssim $ 表示不等式右侧乘以某个仅依赖于前述参数的绝对常数。该估计明确揭示了三重误差源：（i）指数衰减项 $ \exp(-\gamma T) $，控制初始误差与观测有界性贡献，其衰减速率 $ \gamma $ 可任意大，但需增大 $ \Lambda_0 $；（ii）确定性观测误差项，其尺度为 $ \Lambda^2 \int_0^T \mathbb{E}[|e_{D,u}|^2 + |e_{D,\Theta}|^2] $；（iii）随机误差项，其尺度为 $ \Lambda^2 (\sigma_u^2 + \sigma_\Theta^2) \sum_k \int_0^T (\|v_k\|_{L^2}^2 + \|h_k M^{-1}(h_k)\|_{L^1}) dt $，其中 $ M^{-1}(h_k) = h_k + \alpha \fint_Q h_k dx $（见(2.7)）体现MOB系统非局部算子对随机扰动的调制效应。结论成立的关键前提是观测解的 $ L^\infty $-有界性（3.1），此假设在 $ d=2 $ 时由吸引子性质自动满足，在 $ d=3 $ 时则为必要工作假设（参见Fefferman[8]）。该结果首次在三维、随机、非线性、非自治背景下，为OB系统建立了具有显式误差结构的同化收敛性。

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

证明的核心障碍在于：同步系统（3.8）是高度非线性的随机偏微分方程组，其弱martingale解缺乏经典正则性，而观测解 $ (u,\Theta) $ 仅具 $ L^\infty $-有界性，无法直接代入标准能量估计框架。作者的根本创新在于将相对能量不等式（relative energy inequality） 这一源自流体力学变分稳定性的工具，系统性地拓展至随机PDE语境，并针对MOB系统的特殊结构（非局部算子 $ M $、重力耦合项 $ -\Theta \nabla_x G $、非齐次边界处理）进行精密适配。具体而言：

第一，构建MOB专属的相对能量泛函（见(5.1)）：

此泛函不仅包含标准动能差 $ \frac{1}{2}|\tilde{u}-w|^2 $，更关键的是其温度部分 $ (M(\tilde{\Theta}) - M(T))(\tilde{\Theta} - T) $ 精确捕获了MOB系统中非局部算子 $ M $ 的双线性结构，确保其正定性（因 $ M $ 正定）且能与热方程（3.13）中的 $ M(\tilde{\Theta}) $ 项匹配。

第二，将观测解 $ (u,\Theta) $ 作为相对能量中的“参考场” $ (w,T) $。尽管 $ (u,\Theta) $ 仅为有界弱解，但其正则性（由（3.2）保证的 $ W^{2-2/q,q} $ 估计）足以使其满足相对能量推导所需的测试函数条件（5.1），这是证明可行性的基石。

第三，对相对能量泛函应用Itô公式，并严格利用同步系统的弱形式（3.12）–（3.15）与观测解满足的MOB方程（2.4）进行逐项抵消。此过程产生大量交叉项，其关键处理在于：（i）利用 $ M $ 的自伴性与正定性控制 $ \tilde{\Theta} $ 与 $ \Theta $ 的耦合；（ii）通过Gronwall引理与Poincaré不等式将nudging项 $ -\Lambda \int (I_\delta(\tilde{u}) - I_\delta(u)) \cdot (\tilde{u} - u) $ 和 $ -\Lambda \int (I_\delta(\tilde{\Theta}) - I_\delta(\Theta)) (\tilde{\Theta} - \Theta) $ 主导为负定耗散项；（iii）将随机积分项 $ \Lambda \sum_k \int \cdots d\chi_{u,k} $、$ \Lambda \sum_k \int \cdots d\chi_{\Theta,k} $ 通过鞅性质取期望后消失，而其二次变分（即Itô校正项）精确生成误差估计中的 $ \Lambda^2 \sigma_u^2 \sum_k \int |v_k|^2 $ 与 $ \Lambda^2 \sigma_\Theta^2 \sum_k \int h_k M^{-1}(h_k) $ 项，这正是随机误差的定量来源。

第四，引入参数 $ \gamma $ 的精细平衡：在Gronwall处理中，通过添加人工耗散项 $ \gamma \mathcal{E} $ 并选择足够大的 $ \Lambda $ 使总耗散主导所有非负项（包括由 $ \nabla_x G $、$ \vartheta_B $ 等产生的非线性增长项），最终导出以 $ \exp(-\gamma T) $ 为特征的指数收敛率。整个框架完全规避了对同步解正则性的先验要求，仅依赖于弱解的积分形式与相对能量的凸性结构，是处理随机非线性PDE同化问题的范式级技术。

4. 与之前工作的比较

本文与既有工作的本质区别体现在三个维度：

（1）问题设定与适用范围：Azouani-Olson-Titi [2] 开创了确定性连续同化理论，但仅适用于线性或二维Navier-Stokes方程，其nudging项为 $ -\Lambda (\tilde{u} - u) $，要求观测解为强解且系统满足谱间隙条件；Bessaih-Olson-Titi [4] 将其推广至随机情形，但模型为简化版（如忽略重力耦合 $ -\theta \nabla_x G $）且仅处理加性白噪声，未涉及插值算子 $ I_\delta $ 引入的离散化误差。本文首次将同化理论置于完整的三维Oberbeck-Boussinesq框架下，显式包含重力驱动项、非局部边界条件及MOB系统特有的算子 $ M $，其结论直接适用于气象学中真实的对流建模场景。

（2）假设强度与数学代价：[2,4] 要求观测解属于特定Sobolev空间（如 $ H^1 $ 或更高），而本文仅需 $ L^\infty $-有界性（3.1），这在三维中是远弱于经典正则性假设的物理合理条件（如大气流速的实测上界）。代价是收敛速率 $ \exp(-\gamma T) $ 的指数 $ \gamma $ 与nudging参数 $ \Lambda $ 耦合：增大 $ \gamma $ 以加速收敛必须同步增大 $ \Lambda $，导致误差估计中 $ \Lambda^2 $ 倍放大的确定性与随机误差项（见(4.1)）。这揭示了同化精度与鲁棒性间的内在权衡，是此前工作未量化揭示的深层结构。

（3）技术工具与理论深度：[2,4] 依赖经典能量估计与Gronwall不等式，其随机分析限于简单Itô积分。本文的核心工具——为随机PDE定制的相对能量不等式——由Feireisl团队在[5,9]中发展，本文是其首次在数据同化领域的实质性应用。该工具的威力在于：它不要求解的高阶正则性，仅需弱解的积分形式与参考解的有界性，从而绕过了三维OB系统强解存在性这一悬而未决的难题。此外，对非局部算子 $ M $ 的相对能量项设计、对柱状Wiener过程二次变分的精确计算（含 $ M^{-1}(h_k) $ 项），均体现了比前人更精细的分析层次。从提供内容中无法确认本文是否与近期基于机器学习的同化方法（如[12]）进行比较，作者亦未显式展开与[10]（关于MOB吸引子）的关联讨论，其比较焦点明确锁定在经典同化文献[2,4]上。

9. Gevrey 不稳定性在无粘流入-流出问题中

Gevrey instability in the inviscid inflow-outflow problem

作者: Igor Kukavica, Wojciech Ożański, Qi Xu

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04604v1

摘要翻译: 我们在周期通道 $\mathbb{T}\times (0,1)$ 上考虑二维不可压 Euler 方程，边界条件为流入-流出速度条件 $u=(0,1)$，即在 $\mathbb{T} \times \{0,1 \}$ 上成立。我们还在 $\mathbb{T}\times \{ 0 \}$ 上施加入流涡度边界条件 $ω=η$，其中 $η$ 是给定的。我们证明：只要 $η$ 属于 Gevrey 类，则该问题在任意 Gevrey 指数 $s>1$ 对应的 Gevrey 空间中整体适定。这表明：流入-流出速度边界条件局部唯一确定解当且仅当解被置于解析类中。特别地，一旦离开解析类，通过适当选取 $η$，已在任意 Gevrey 类中出现解的非唯一性。因此，若指定解析的流入-流出速度，则恰好存在唯一一个解析解及一族连续依赖于 $s$ 的 $s$-Gevrey 解（对每个 $s>1$）。此外，该结果蕴含：若 $η$ 是解析的，则唯一整体解可在所有 $t>0$ 时刻失去空间解析性，但仍对所有 $s$ 保持 $s$-Gevrey 正则性。

Comment: 25 pages

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

1. 研究对象和背景

本文研究二维不可压Euler方程在周期通道 $\Omega = \mathbb{T} \times (0,1)$ 上的流入-流出边值问题：

其中速度场 $u = (u_1, u_2)$ 满足刚性流入-流出边界条件

即法向速度恒为1；同时在下边界 $S := \mathbb{T} \times \{0\}$（流入端）施加入流涡量边界条件

注意：上边界 $\widetilde{S} = \mathbb{T} \times \{1\}$（流出端）不施加涡量条件，仅由 $u_2=1$ 约束。该设定源于物理与数学双重动机：一方面，经典无粘流入-流出问题（如Zaj\c{a}czkowski [Z1–Z4]、Antontsev-Kazhikov-Monakhov [AKM1, AKM2]、Petcu [P]、Gie-Kelliher-Mazzucato [GKM1, GKM2]等系列工作）长期面临根本性适定性障碍——因Euler方程在边界上导出的椭圆压力系统与Gromeka-Lamb形式共同给出关系式（文中(1.4)）

而自由度为四，故需额外两个边界约束才能封闭系统。在三维情形中，这体现为入流部分要求 $\omega$ 属于 $\mathrm{curl}$ 的像且 $\mathrm{div}\, \omega$ 守恒（见[GKM2, (1.6)]）。现有两类主流处理方式：(i) 在入流边界完整指定 $u$（如[GKM2, Theorem 1.2]在 $C^{N+1,\alpha}$ 中建立局部适定性，但需高阶相容性条件）；(ii) 指定入流涡量（如[Z1]在3D Sobolev空间、[GKM2, Theorem 1.3]在3D $C^{k,\alpha}$ 中建立局部适定性）。本文聚焦二维简化模型，其核心结构由涡量输运方程 $\partial_t \omega + u \cdot \nabla \omega = 0$ 与Biot-Savart重构 $u = u[\omega]$ 所决定：给定 $\omega$，速度通过求解如下耦合椭圆系统唯一确定

该重构隐含关键事实：$u$ 的正则性由 $\omega$ 控制，且 $\omega$ 在流出边界 $\widetilde{S}$ 的行为不贡献能量增长（见(1.13)–(1.14)），故仅需在 $S$ 施加 $\eta$ 即可驱动系统演化。本文试图填补的关键空白是：在非解析正则性框架下，流入-流出速度边界条件是否足以唯一确定解？ 已知[KOS]证明：若 $u$ 初始为解析函数（$s=1$ Gevrey类），则仅凭 $u_2=1$ 即可获得局部唯一解析解（Theorem 1.1），无需任何相容性条件；但该结果对解析性具有本质依赖。本文首次严格证明：一旦离开解析类（即进入任意 $s>1$ Gevrey类），仅指定 $u_2=1$ 将导致无限多解——其自由度恰好由入流涡量 $\eta$ 的Gevrey正则性所参数化。因此，问题本质从‘速度边界条件是否充分’转变为‘何种正则性层级下边界信息才足以封闭系统’。

2. 主要定理或主要结果

定理1.2（Gevrey空间中的整体适定性）：设 $s > 1$，若入流涡量 $\eta$ 满足Gevrey型估计

且初始涡量 $\omega_0$ 满足 $\|D^\alpha \omega_0\|_{L^2(\Omega)} \leq N_0(0) N(0)^{|\alpha|} |\alpha|!^s$（$\forall \alpha \in \mathbb{N}_0^3$），则存在唯一全局解 $u \in C([0,\infty); H^k(\Omega))$（$\forall k \geq 0$），其涡量 $\omega$ 满足Gevrey正则性估计

其中 $L_0(t), L(t) \geq 1$ 仅依赖于 $N_0, N$ 及 $s$。特别地，$L(t)$ 具体形式为 $L(t) = K(t) \tau(t)^{-1}$，其中Gevery半径 $\tau(t)$ 满足指数衰减律

且 $K(t) = J_2(t)$ 是由 $\eta$ 的前二阶Gevrey常数控制的连续增函数。推论1.3（Gevrey空间中的非唯一性）：对任意 $s > 1$，若仅给定流入-流出速度 $u_2 = 1$（即问题(1.1)–(1.3)），则对任一初始数据 $u_0 \in G^s$，均存在无穷多个解 $u(t) \in G^s$（$\forall t > 0$）；其构造方式为：任取满足上述Gevrey条件的 $\eta$，由定理1.2生成对应解。推论1.4（解析初值产生非解析解）：若 $\omega_0 = 0$ 且 $\eta(t)$ 对每个 $t \in (0,T]$ 是空间解析的（即满足 $s=1$ 的Gevrey估计），则由定理1.2生成的唯一解 $u(t)$ 在任意 $t > 0$ 均不解析（即其Gevrey半径 $\tau(t) \to 0$），尽管它对所有 $s > 1$ 仍保持 $G^s$ 正则性。

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

证明的核心障碍在于：Gevrey框架下，标准能量估计无法控制高阶导数的阶乘爆炸增长，且边界涡量输入 $\eta$ 在 $S$ 上的Gevrey正则性需通过非局部Biot-Savart重构精确传递至内部 $\omega$ 和 $u$。作者的真正创新点在于构建了一套双尺度Gevrey系数递归控制体系，其关键工具与技术包括：（1）边界-内部Gevrey系数分离：定义边界系数 $a_\alpha(t) := \frac{\epsilon^\alpha}{|\alpha|!^s \tau(t)^{|\alpha|-2}_+} \|D^\alpha \omega\|_{L^2(S)}$ 与内部系数 $b_\alpha(t) := \frac{\epsilon^\alpha}{|\alpha|!^s \tau(t)^{|\alpha|-2}_+} \|D^\alpha \omega\|_{L^2(\Omega)}$（见(2.1)–(2.2)），其中 $\epsilon = (1, \epsilon_2(t), 1)$ 且 $\epsilon_2(t) = \frac{1}{4K(t)}$ 动态衰减，$\tau(t)$ 为时变Gevrey半径。此分离将边界输入的正则性与内部传播的耗散效应解耦。（2）速度的阶数削减估计：证明关键不等式（(2.4)–(2.5)）

该估计揭示了 $u$ 的 $|\alpha|$ 阶导数由 $\omega$ 的至多 $|\alpha|-1$ 阶导数控制，是连接涡量Gevrey范数与速度非线性项的关键桥梁。（3）非线性项的Gevrey压缩引理：对非线性项 $NL(\alpha) = -\frac{\epsilon^{2\alpha}}{|\alpha|!^{2s} \tau(t)^{2|\alpha|-4}} \int_\Omega \left(D^\alpha(u \cdot \nabla \omega) - u \cdot \nabla D^\alpha \omega\right) D^\alpha \omega$，作者利用组合恒等式（引理3.1）及 $s>1$ 的凸性，证明（(2.16)）

其中 $b_{m+1}(t) = \max_{|\beta| \leq m+1} b_\beta(t)$。该估计的关键在于：当 $s>1$ 时，$\binom{n}{k}^{1-s}$ 的衰减足以压制 $k$ 阶导数与 $n-k$ 阶导数的乘积带来的阶乘增长，从而实现非线性项对高阶Gevrey系数的“压缩”效应。（4）动态Gevrey半径设计：选择 $\tau(t)$ 满足 $\dot{\tau}/\tau = -8\tilde{C} K(t)^6$，使得在Gevrey能量不等式（(2.15)）中，耗散项 $|\alpha| \frac{\dot{\tau}}{\tau} b_\alpha^2$ 主导非线性项 $NL(\alpha)$ 与边界项 $a_\alpha^2$，从而保证 $b_\alpha(t)$ 的一致有界性（见(2.18)–(2.19)）。这一设计将Gevrey正则性的维持转化为对 $K(t)$ 积分增长的精细控制。

4. 与之前工作的比较

本文与既有工作的差异体现在三个根本维度：（1）正则性阈值的精确刻画：[KOS]证明解析类（$s=1$）下 $u_2=1$ 足以保证局部唯一性，但未探讨 $s>1$ 的失效机制；本文首次严格证明 $s=1$ 是临界阈值——对任意 $s>1$，$u_2=1$ 均不足以唯一确定解，且非唯一性由 $\eta$ 的Gevrey参数连续谱实现。这比[GKM2]在Hölder空间的局部适定性更精细，因其揭示了正则性损失的定量结构。（2）边界条件策略的根本转向：此前工作（如[Z1]、[GKM2, Thm 1.3]）将入流涡量 $\eta$ 视为辅助工具以恢复适定性，其目标仍是证明存在唯一解；本文则将 $\eta$ 提升为自由参数，证明其Gevrey正则性直接参数化解空间，从而确立‘Gevrey不稳定性’的本质——适定性丧失并非因技术缺陷，而是Gevrey类内在的非解析性所致。（3）时间尺度与适用范围的拓展：[KOS]的解析解仅在有限时间 $T_0 \sim \tau_0/M$ 内存在；本文定理1.2给出全局Gevrey解，其Gevrey半径 $\tau(t)$ 虽指数衰减（$\tau(t) \sim e^{-C\int_0^t K^6}$），但对任意 $t>0$ 保持正，故解在所有时刻均属 $G^s$（$s>1$）。代价在于：解析性永久丧失（推论1.4），即 $\tau(t) \to 0$ 作为 $t \to \infty$ 的必然结果，这与[KOS]中解析解可能长期保持解析形成鲜明对比。局限性在于：结果局限于二维周期通道及恒定流入-流出速度 $u_2=1$；对一般 $u_2 = g(x_1)$ 或非周期域，$u[\omega]$ 的重构及Gevrey估计需全新分析框架。此外，作者未显式展开与Sobolev框架（如[Z1], [GKM2]）的比较，但从方法看，Gevrey估计蕴含所有Sobolev范数有界，故本文结果强于纯Sobolev适定性，但弱于解析适定性。

10. 基于梯度流的欧拉型非线性弹性与扩散的离散化

Discretisation of Eulerian nonlinear elasticity and diffusion using gradient flows

作者: Andrea Zafferi, Dirk Peschka

类别: math.NA, math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03829v1

摘要翻译: 在本研究中，我们提出了一种通用的能量型建模方法，用于具有拉格朗日与欧拉两种构型下扩散输运特性的粘性多孔弹性材料。我们的研究通过在欧拉构型中引入参考映射（reference map）概念，导出了更精细的弱形式。我们提出并实现了一种新型的结构保持型离散策略，采用混合有限元方法。本文通过对拉格朗日与欧拉格式进行对比分析，展示了所提方法在空间与时间上的数值收敛性，从而证明了该方法的鲁棒性与实用性。此外，在欧拉型多相流（特别是准静态欧拉–欧拉型）背景下，我们的研究证明了在多孔粘弹性介质中存在孤波型流体波（solitary fluid waves）。这一能量型方法为耦合多孔弹性、流动与扩散的热力学建模及其相应离散格式提供了理论基础。

Comment: 无

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

1. 研究对象和背景

本文研究的核心对象是粘弹性多孔介质中耦合非线性弹性变形与扩散输运的动力学模型，其数学本质是一类在欧拉坐标系下定义的、具有明确热力学结构的非线性偏微分方程组。该模型以能量变分框架为出发点，将系统状态变量取为欧拉位移场 $u(t,x) \in H^1_0(\Omega; \mathbb{R}^d)$ 和（隐含的）参考映射 $\alpha(t,x)$，其中 $\Omega \subset \mathbb{R}^d$（文中重点处理 $d=2$）为固定空间域，对应于变形后构型；而参考构型 $\bar{\Omega}$ 通过可逆流图 $\bar{\chi}(t,\bar{x})$ 与之关联，满足 $\Omega = \bar{\chi}(t,\bar{\Omega})$。关键物理量包括欧拉变形梯度 $F(t,x) = (\nabla \alpha)^{-1} = (\mathrm{Id} - \nabla u)^{-1}$，其雅可比行列式 $J(t,x) = \det F(t,x)$ 表征局部体积变化；以及压缩性新胡克型自由能密度

其中 $C = F^\top F$ 为右柯西-格林张量，$\mu > 0$、$\lambda > 0$ 分别为剪切模量与体积模量。总自由能泛函为

这里 $f(x)$ 代表外体力密度（如重力）。系统演化由梯度流结构 $(Q, H, R)$ 控制：状态空间 $Q = \{ u \in H^1_0(\Omega; \mathbb{R}^d) \}$，能量泛函 $H$ 如上，耗散势 $R$ 取 Darcy 型形式

其中 $\kappa(u) > 0$ 为渗透率。由此导出的弱形式演化方程为

其强形式为

其中 $\sigma = J^{-1} \partial_F \bar{H}_{\text{elastic}} F^\top$ 是柯西应力张量。该模型的重要性在于：（1）它统一处理大变形（finite strain）下的粘弹性响应与孔隙流体/溶质扩散，超越了小应变 Biot 理论的适用范围；（2）其欧拉表述避免了 Lagrangian 方法中网格畸变导致的数值不稳定性，对地质学中的孔隙波（porosity waves）和生物软组织建模尤为关键；（3）它天然嵌入 GENERIC 等热力学框架，保证了能量耗散律 $\frac{d}{dt} H(u(t)) = - \| \dot{u} \|_{R}^2 \leq 0$ 的离散保持。然而，现有文献中，欧拉框架下的非线性弹性-扩散耦合缺乏严格、一致的变分结构推导，尤其在参考映射的引入、能量与耗散在坐标变换下的结构保持、以及混合有限元离散的鞍点结构构造等方面存在理论空白。本文旨在填补这一空白，建立从 Lagrangian 到 Eulerian 梯度流的系统性、结构保持的变分映射，并据此构建鲁棒的数值方案。

2. 主要定理或主要结果

本文在以下核心假设下证明了若干关键结论：(i) 域 $\Omega \subset \mathbb{R}^2$ 为有界 Lipschitz 区域；(ii) 材料参数 $\mu, \lambda > 0$，渗透率 $\kappa(u) \equiv 1$（或满足 Lipschitz 连续性）；(iii) 初始位移 $u_0 = 0$，边界条件为齐次 Dirichlet 条件 $u|_{\partial \Omega} = 0$；(iv) 时间离散采用一阶向后欧拉格式，空间离散采用 $P_k$ ($k=1,2$) 有限元。主要结果如下：

定理 3.1（能量最小化解的一致性）：设 $\bar{u} \in \bar{Q} = H^1_0(\bar{\Omega}; \mathbb{R}^2)$ 和 $u \in Q = H^1_0(\Omega; \mathbb{R}^2)$ 分别为 Lagrangian 和 Eulerian 自由能泛函 $\bar{H}$ 和 $H$ 的极小值点，则存在 $L^2$ 投影 $\operatorname{proj}_{L^2}(\bar{u} \circ \bar{\chi}^{-1})$，使得误差估计成立

其中 $h$ 为网格尺寸，$C$ 为与 $h$ 无关的常数。图3的数值实验验证了该理论收敛阶。

定理 3.2（梯度流解的结构保持离散一致性）：对于时间步长 $\tau$ 和空间步长 $h$，Lagrangian 离散格式 (24) 与 Eulerian 离散格式 (25) 的解序列 $\{\bar{u}_h^n\}$ 和 $\{u_h^n\}$ 满足离散能量耗散律

且二者在 Bochner 范数 $L^2(0,T; L^2(\Omega))$ 下的差值满足

定理 4.1（孔隙波存在性）：在准静态 Euler-Euler 型多相流背景下，针对粘弹-多孔介质，本文通过分析该梯度流系统的行波解，证明了孤立孔隙波（solitary porosity wave）的存在性。具体而言，存在非平凡解 $u(z), \phi(z)$（其中 $z = x - ct$ 为行波坐标），满足耦合方程

并具有渐近行为 $\lim_{|z| \to \infty} u(z) = 0$, $\lim_{|z| \to \infty} \phi(z) = \phi_\infty$，其波形由能量泛函 $H$ 的非凸性及耗散项共同决定。

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

本文的核心技术障碍在于：如何在非线性、大变形的欧拉框架下，严格地、结构保持地将 Lagrangian 梯度流 $(\bar{Q}, \bar{H}, \bar{R})$ 映射为 Eulerian 梯度流 $(Q, H, R)$，并确保该映射能自然导出稳定的混合有限元离散格式。作者的真正创新点并非提出新的 PDE 模型，而是构建了一套系统性的、基于变分原理的坐标变换理论，其关键工具与分析框架如下：

1. 结构保持的变分映射框架（Section 2.2）：这是全文最核心的创新。作者没有采用传统的协变导数或 Lie 导数等张量分析方法，而是发展了一套抽象的“闭包条件”（closure condition）体系。对于状态变量映射 $q = T(\bar{q})$，要求能量、耗散势、双线性形式均满足特定的相容性条件。例如，对于耗散势，要求 $R(T(\bar{q}), L\bar{v}) = \bar{R}(\bar{q}, \bar{v})$，其中 $L = D_{\bar{q}}T$ 是 Fréchet 导数（即变形梯度的欧拉表示 $F = (\nabla \alpha)^{-1}$）。这一条件直接导出了欧拉耗散势 $R_{\text{Darcy}}(u,w) = \frac{1}{2} \int_\Omega \kappa(u) |F w|^2 dx$，其被积函数中显式出现的 $F$ 因子是传统方法无法自然产生的，却是保证能量耗散律 $\frac{d}{dt}H = -\|\dot{u}\|_R^2$ 成立的关键。该框架还区分了两种映射方法：Method 1 适用于已知所有双线性形式的情形；Method 2（文中最终采用）则允许在未知欧拉双线性形式 $b$ 的情况下，通过指定辅助空间 $W$ 上的映射 $L_W$ 来反推 $b$，这为处理含对流导数的复杂问题提供了灵活性。

2. 参考映射 $\alpha$ 的双重角色：$\alpha$ 不仅是坐标变换的数学工具，更是欧拉框架下定义变形梯度 $F = (\nabla \alpha)^{-1}$ 和雅可比 $J = \det F$ 的物理基础。作者通过恒等式 $\partial_t \alpha = -(\nabla \alpha) v$ 和 $\partial_t u = -\partial_t \alpha$，将欧拉速度场 $v$ 与位移场 $u$ 的时间演化紧密耦合，从而绕开了在欧拉坐标下直接定义物质导数 $D/Dt$ 的困难。这一技巧使得整个梯度流结构能在纯欧拉变量下自洽表述。

3. 鞍点问题的结构化离散：作者将梯度流方程 (4) $\dot{q} = -K(q)DH(q)$ 重构为弱形式的鞍点问题 (6)-(7)，其中辅助变量 $\eta \in W$ 对应广义力，双线性形式 $k(\eta,w)$ 和 $b(v,\eta)$ 构成一个对称正定系统。这种结构天然适配混合有限元法，且能严格保持离散能量耗散律 (8)。数值实现中，作者利用 FEniCS 框架自动微分功能，直接对泛函 $H$ 和 $R$ 进行变分，避免了传统方法中手动推导复杂应力表达式的错误风险。

4. 与之前工作的比较

本文与既有工作的比较需从多个维度展开：

与经典 Lagrangian 方法（如 [17,18,19]）的比较：Lagrangian 方法在数学上成熟，能自然处理大变形，但其数值实现严重依赖于网格随材料点移动，易在大变形下产生网格畸变和单元翻转，导致计算崩溃。本文的 Eulerian 方法将计算域固定在空间 $\Omega$ 上，完全规避了此问题，显著提升了数值鲁棒性，代价是引入了参考映射 $\alpha$ 这一额外的、非局部的变量，增加了问题的耦合复杂度。

与传统 Eulerian 张量方法（如 [20,15,21,22]）的比较：后者通常基于左/右 Cauchy-Green 张量 $B = FF^\top$ 或 $C = F^\top F$ 的协变 Lie 导数来构建本构方程，虽能保持客观性，但其变分结构不清晰，难以保证能量耗散律的离散保持。本文通过结构保持映射，将 Lagrangian 的清晰能量结构完整地移植到 Eulerian 框架，使 $F$ 和 $J$ 的定义直接源于 $\alpha$，而非作为独立张量场求解，从而在保证物理正确性的同时，获得了严格的热力学一致性。

与基于 GENERIC 的离散工作（如 [39,40,41,42]）的比较：GENERIC 框架强调可逆与不可逆过程的分离，但其离散化常面临算子 $K(q)$ 的构造难题。本文的梯度流框架 $(Q,H,R)$ 是 GENERIC 在纯耗散情形下的特例，其优势在于 $R$ 的凸性保证了 $K(q)$ 的良定性，且作者提出的闭包条件为 $K(q)$ 在坐标变换下的映射提供了普适准则，这是 GENERIC 文献中未系统解决的。

与 Biot 理论及小应变模型（如 [31,32,33]）的比较：Biot 理论是线性化的基石，但其失效于 $|\nabla u| \gtrsim 0.3$ 的有限应变 regime。本文的压缩性新胡克模型能精确刻画 $|\nabla \bar{u}| \approx 1.3$ 的强变形（见图2），将适用范围从 infinitesimal strain 推进到 finite strain，代价是方程的非线性程度剧增，对数值方法的收敛性和稳定性提出了更高要求。

与单纯数值比较工作（如 [43,44]）的比较：文献 [43,44] 也探讨了坐标变换，但其焦点在于 Hamiltonian 系统，且未深入到耗散结构 $R$ 的映射细节。本文首次系统地、定量地给出了耗散势 $R$ 在 Lagrangian-Eulerian 变换下的闭包条件 (10)，并将其应用于 Darcy 型耗散，这是对前人工作的实质性推进。从提供内容中无法确认本文方法与更前沿的任意拉格朗日-欧拉（ALE）或浸入边界法（IBM）在相同问题上的精度与效率对比。

11. 关于 Navier-Stokes 方程的极弱适定解与 Leray 解的一个注记

A remark on very weak suitable solutions and Leray solutions of the Navier-Stokes equations

作者: Oscar Jarrín

类别: math.AP

发布日期: 2026-07-03

arXiv ID: 2607.03602v1

摘要翻译: 我们引入 Navier-Stokes 方程的极弱适定解（very weak suitable solutions）的概念。此处，速度场与压力满足使局部能量不等式在分布意义下成立的最小条件。一个广为人知但至今仍具相关性的问题是：寻找充分条件，以保证极弱解实际上就是 Leray 解。我们在局部 Morrey 空间的通用框架内，利用其局部能量不等式，建立了此类充分条件。局部 Morrey 空间提供了一个通用框架，其中包含了 Navier-Stokes 方程理论分析中其他有用的功能性设定，例如 Lebesgue 空间、Lorentz 空间、齐次 Morrey 空间以及抛物型 Morrey 空间。作为副产品，我们还导出了一些用于研究所得 Leray 解之唯一性与正则性的充分条件。

Comment: 17 pages

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

1. 研究对象和背景

本文研究三维全空间 $\mathbb{R}^3$ 上不可压Navier-Stokes方程的解的正则性提升问题，核心方程为

其中 $\vec{u}: [0,T] \times \mathbb{R}^3 \to \mathbb{R}^3$ 为速度场，$P: [0,T] \times \mathbb{R}^3 \to \mathbb{R}$ 为压力，所有物理常数归一化为1。研究对象是两类解之间的关系：一类是仅满足最弱分布意义下可定义性的“极弱解”（very weak solution），另一类是经典Leray-Hopf弱解（Leray solution）。前者由Definition 1.1定义：要求 $\vec{u} \in L^2_{\mathrm{loc}}([0,T] \times \mathbb{R}^3)$、初始值 $\vec{u}_0 \in \mathcal{D}'(\mathbb{R}^3)$，且对任意 $\vec{\varphi} \in \mathcal{D}([0,T] \times \mathbb{R}^3)$ 满足 $\mathrm{div}\, \vec{\varphi} = 0$，有

此处压力 $P$ 由 $-\partial_t \vec{u} + \Delta \vec{u} - \mathrm{div}(\vec{u} \otimes \vec{u})$ 的无旋性隐式确定（[21, Lemma 6.3]）。该定义允许 $\vec{u}_0$ 为任意分布，故解的存在性无需能量控制，但缺乏先验估计。而Leray解（Definition 1.2）要求更强的函数空间正则性：$\vec{u} \in L^\infty([0,T]; L^2(\mathbb{R}^3)) \cap L^2([0,T]; \dot{H}^1(\mathbb{R}^3))$，并满足全局能量不等式

该不等式蕴含 $\vec{u}$ 在 $L^2$ 意义下连续于 $t=0$，且整体受控于初始能量。一个长期悬而未决的基础问题是：在何种附加条件下，一个极弱解自动成为Leray解？这本质上是正则性提升（regularity upgrading）问题——从仅能定义方程的分布解，跃升至具有能量守恒结构和Sobolev可微性的经典弱解。已有工作多依赖Ladyzhenskaya-Prodi-Serrin型条件 $\vec{u} \in L^p_t L^q_x$ 满足 $2/p + 3/q = 1$，例如[12]处理亚临界情形 $2/p + 3/q < 1$；[8]在临界情形下额外要求解满足mild公式；[16]引入加权连续性 $t^{1/p} \vec{u} \in C([0,T]; L^q)$；[14]允许初值仅属 $L^2$ 但需在 $[\delta,T]$ 上满足 $L^p_t L^q_x$；[7]甚至处理超临界 $L^4_t L^4_x$（但需强收敛初值）。然而，这些条件均局限于Lebesgue或Lorentz空间框架，未能统一涵盖具有不同衰减行为的解（如各向异性衰减、局部奇点强度等）。本文试图填补的关键空白在于：建立一个统一的、更广谱的函数空间框架，使得在该框架下，结合Navier-Stokes方程固有的局部能量结构，即可系统地实现从极弱解到Leray解的跃迁。为此，作者引入“极弱适定解”（very weak suitable solution）概念（Definition 1.3）：要求 $\vec{u} \in L^3_{\mathrm{loc}}$、$P \in L^{3/2}_{\mathrm{loc}}$，且由局部能量等式导出的分布

是一个非负局部有限测度。这一条件直接继承自Caffarelli-Kohn-Nirenberg理论中对Leray解的“适定性”（suitability）要求，但不要求 $\vec{u}$ 属于能量空间，从而将适定性结构推广至极弱解范畴。问题转化为：在何种泛函空间假设下，该非负测度 $\mu \geq 0$ 足以强制 $\vec{u}$ 满足全局能量不等式（3）及相应Sobolev正则性？本文选择局部Morrey空间 $M^{p,\gamma}_{t,x}$ 作为载体，其范数定义为

该空间刻画函数在球 $|x| < R$ 上的平均 $L^p$ 范数随 $R$ 的衰减速率，参数 $\gamma$ 控制渐近行为。它天然包含Lebesgue空间（$\gamma = 3$）、加权Lebesgue空间（$L^p_\gamma$）、Lorentz空间、齐次Morrey空间 $\dot{M}^{p,q}$ 及抛物Morrey空间 $M^{p,q}_{t,x}$，构成一个高度统一的分析平台。因此，本文的研究背景是：在Navier-Stokes理论中，如何利用局部能量不等式的结构性约束（即 $\mu \geq 0$），在一个比传统 $L^p_t L^q_x$ 更广泛、更具几何适应性的函数空间（局部Morrey空间）中，建立极弱解向Leray解的正则性跃迁准则。

2. 主要定理或主要结果

本文的核心结果是Theorem 1.1：设 $(\vec{u}, P)$ 是Navier-Stokes方程(1)的极弱适定解（Definition 1.3），若满足以下条件，则 $(\vec{u}, P)$ 是Leray解（Definition 1.2）。

(1) 初始值 $\vec{u}_0 \in L^2(\mathbb{R}^3)$；且对任意紧集 $\Omega \subset \mathbb{R}^3$，映射 $t \mapsto \vec{u}(t, \cdot)$ 在区间 $(0,T]$ 上关于 $L^2(\Omega)$ 弱拓扑连续，在 $t=0$ 处强连续。

(2) 速度场 $\vec{u}$ 属于局部Morrey空间 $M^{p,\gamma}_{t,x}([0,T] \times \mathbb{R}^3)$，其中参数满足 $0 < \gamma < 3 \leq p < \infty$ 且

特别地，当 $p=3$ 时，条件简化为 $\gamma < 1$；当 $p \to 9/2^-$ 时，$\gamma \to 0^+$。结合 $0<\gamma<3$ 和 $p \geq 3$，参数范围被严格限定为 $3 \leq p < 9/2$。

该定理的证明实质上建立了全局能量不等式(3)：通过构造一个时空截断测试函数 $\Phi_{t_0,t_1,\varepsilon,R}(t,x) = \alpha_{t_0,t_1,\varepsilon}(t) \varphi_R(x)$，其中 $\alpha$ 在时间上光滑截断于 $[t_0,t_1]$，$\varphi_R$ 在空间上光滑截断于 $|x| < R$，并将 $\Phi$ 代入局部能量等式(6)。利用 $\mu \geq 0$ 的非负性，得到关键的局部能量估计：

令 $R \to +\infty$，第二项因 $R^{-1/3} \to 0$ 消失；第三项的指数 $3(\frac{\gamma}{p} - \frac{3}{p} + \frac{2}{3}) < 0$ 由条件(8)保证，故亦趋于0，从而得到 $\|\vec{u}(t, \cdot)\|_{L^2}^2 + 2 \int_0^t \|\nabla \otimes \vec{u}(s, \cdot)\|_{L^2}^2 \, ds \leq \|\vec{u}_0\|_{L^2}^2$，即(3)成立。进而由标准论证（如[21, Theorem 12.2]）可将解延拓至全局时间。此外，Proposition 1.1给出了压力项的刻画：若极弱解 $(\vec{u}, P)$ 满足 $\vec{u} \in M^{p,\gamma}_{t,x}([0,T] \times \mathbb{R}^3)$，其中 $0 < \gamma < 3/2$ 且 $2 < p < \infty$，则压力 $P$ 必然由Riesz变换表示为

Corollary 1.1给出Liouville型唯一性：若 $\vec{u}_0 = 0$，则 $\vec{u} \equiv 0$。Corollary 1.2和1.3则提供了进一步的唯一性与正则性：若在更小的时间区间 $[0,T_1]$ 上，$\vec{u}$ 还属于另一局部Morrey空间 $M^{q,\delta}_{t,x}$，其中 $q \geq \frac{3p}{p-2}$ 且 $\delta > 0$，则 $\vec{u}$ 在 $[0,T_1]$ 上唯一，且 $\vec{u} \in C^\infty((0,T_1] \times \mathbb{R}^3)$。

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

本文证明的核心障碍在于：极弱适定解仅知 $\vec{u} \in L^3_{\mathrm{loc}}$，远不足以直接控制非线性项 $\mathrm{div}(\vec{u} \otimes \vec{u})$ 在能量估计中的贡献；同时，局部能量不等式(6)中的测度 $\mu$ 仅提供非负性信息，而非具体的大小控制。传统方法（如在 $L^p_t L^q_x$ 中使用插值或Gagliardo-Nirenberg不等式）在此失效，因为局部Morrey空间不具备尺度不变性，且其范数涉及空间平均而非逐点控制。作者的真正创新点在于构建了一套基于局部Morrey空间特性的、精细的截断-极限分析框架，其关键工具与技术分述如下。

第一，时空双截断测试函数的构造与优化。证明中使用的测试函数 $\Phi_{t_0,t_1,\varepsilon,R}(t,x) = \alpha_{t_0,t_1,\varepsilon}(t) \varphi_R(x)$ 具有双重目的：时间截断 $\alpha_{t_0,t_1,\varepsilon}$ 确保能量估计在任意子区间 $[t_0,t_1]$ 上成立，并在 $\varepsilon \to 0$ 时恢复指示函数 $\mathbf{1}_{[t_0,t_1]}$；空间截断 $\varphi_R$ 将估计局限在球 $|x| < R$ 内，其梯度与Laplacian满足 $\|\nabla \varphi_R\|_{L^\infty} \lesssim R^{-1}$、$\|\Delta \varphi_R\|_{L^\infty} \lesssim R^{-2}$。这一构造将非线性项的贡献精确地分离为三部分：来自 $\partial_t |\vec{u}|^2$ 的边界项（由时间截断处理）、来自 $|\nabla \otimes \vec{u}|^2$ 的耗散项（保留为正）、以及来自 $\mathrm{div}(|\vec{u}|^2 \vec{u} + 2P \vec{u})$ 的通量项（需用Morrey范数控制）。该函数的设计直接服务于后续的尺度分析。

第二，局部Morrey空间上的关键嵌入与算子有界性。证明中反复依赖Lemma 2.1和Lemma 2.2所建立的深刻性质。Lemma 2.1表明 $M^{p,\gamma}_{t,x}$ 可被加权Lebesgue空间夹逼：$L^p([0,T]; L^p_\gamma(\mathbb{R}^3)) \subset M^{p,\gamma}_{t,x} \subset L^p([0,T]; L^p_\delta(\mathbb{R}^3))$ 对 $0 < \gamma < \delta$ 成立。这使得作者能将 $\vec{u}$ 的Morrey范数控制转化为对加权空间范数的控制，为后续应用Hardy-Littlewood极大算子 $M$ 和Riesz变换 $R_i$ 的有界性铺平道路。Lemma 2.2的核心结论是：当 $0 < \gamma < 3$ 且 $1 < p < \infty$ 时，$M$ 和 $R_i$ 在 $M^{p,\gamma}_{t,x}$ 上连续，即 $\|R_i \phi\|_{M^{p,\gamma}_{t,x}} \lesssim \|\phi\|_{M^{p,\gamma}_{t,x}}$。这一有界性并非平凡，其证明基于实插值理论：$M^{p,\gamma}_{t,x} = [L^p_t L^p_x, L^p_t(L^p_\delta)_x]_{\gamma/\delta, \infty}$，而 $R_i$ 在端点空间 $L^p_t L^p_x$ 和 $L^p_t(L^p_\delta)_x$ 上均有界（由经典理论），故在插值空间上亦有界。此性质是Proposition 1.1得以成立的基石。

第三，非线性项的Morrey尺度精细估计。在局部能量估计中，最关键的非线性项是 $\int_0^t \int_{|x|<R} \mathrm{div}(|\vec{u}|^2 \vec{u}) \Phi \, dx ds$。作者并未简单粗暴地使用Hölder不等式，而是利用测试函数的梯度估计 $|\nabla \Phi| \lesssim R^{-1}$ 和 $\varphi_R$ 的支撑性质，将其上界化为

随后，利用Morrey范数定义(7)，对环形区域 $R/2 < |x| < R$ 上的积分进行控制：

对于一般 $p \geq 3$，则需用Hölder不等式将 $|\vec{u}|^3$ 拆为 $|\vec{u}|^{p} \cdot |\vec{u}|^{3-p}$，再结合Morrey嵌入 $M^{p,\gamma}_{t,x} \subset L^p_t L^p_x$ 和加权空间嵌入，最终导出估计中出现的 $R^{3(\gamma/p - 3/p + 2/3)}$ 项。这一指数的符号完全由条件(8)决定，是整个论证的枢纽。

第四，压力项的Riesz表示的Morrey框架重建。Proposition 1.1的证明是独立的技术亮点。其核心障碍是：在 $M^{p,\gamma}_{t,x}$ 中，$\vec{u} \otimes \vec{u}$ 仅属 $L^{p/2}_t L^{p/2}_{\mathrm{loc}}$，而Riesz变换 $R_i R_j$ 作用其上需更强的可积性。作者的创新在于引入插值思想：先由 $\vec{u} \in M^{p,\gamma}_{t,x}$ 推出 $\vec{u} \in L^p_t L^p_\delta(\mathbb{R}^3)$（Lemma 2.1），再得 $\vec{u} \otimes \vec{u} \in L^{p/2}_t L^{p/2}_\delta(\mathbb{R}^3)$，最后利用 $R_i R_j$ 在加权空间 $L^{p/2}_\delta$ 上的有界性（Lemma 2.2）得到 $Q = \sum R_i R_j(u_i u_j) \in L^{p/2}_t L^{p/2}_\delta(\mathbb{R}^3)$。进而证明 $\nabla(P-Q)$ 是调和多项式，再由其属于 $L^p_\delta + L^{p/2}_\delta$ 而必为零，从而确立(14)。这一论证将经典压力表示理论成功移植到全新的Morrey框架。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开：函数空间框架、假设强度、时间尺度与适用范围、以及代价与局限。

首先，在函数空间框架上，本文实现了根本性拓展。此前所有相关工作（[12, 8, 16, 14, 7, 15]）均局限于标准Lebesgue空间 $L^p_t L^q_x$ 或其变体（如Lorentz空间 $L^{q,r}_x$）。例如，[14]要求 $\vec{u} \in L^p_t L^q_x$ 满足 $2/p + 3/q = 1$；[7]处理 $L^4_t L^4_x$；[15]处理 $L^4_t L^4_x$ 并获得能量等式。而本文的 $M^{p,\gamma}_{t,x}$ 是一个统一的母空间，其嵌入关系（9）–（11）明确显示：$L^p_t L^q_x \subset L^p_t L^{q,r}_x \subset M^{p,\gamma}_{t,x}$、$L^p_t \dot{M}^{p,q}_x \subset M^{p,\gamma}_{t,x}$、$M^{p,q}_{t,x} \subset M^{p,\gamma}_{t,x}$（经Lemma 1.1）。这意味着，任何满足前述经典条件的解，只要其参数 $p,q$ 符合(8)，自动属于 $M^{p,\gamma}_{t,x}$，从而被本文结论所覆盖。更重要的是，$M^{p,\gamma}_{t,x}$ 能容纳 $L^p_t L^p_\gamma(\mathbb{R}^3)$（加权空间），这允许解在无穷远处具有幂律衰减 $|x|^{-\gamma}$，而 $L^p_t L^q_x$ 无法刻画此类各向异性行为。因此，本文框架显著拓宽了可处理的解的类别。

其次，在假设强度上，本文对初值和解的正则性要求更为精炼。以往工作大多要求初值 $\vec{u}_0$ 属于更强的空间：[12]要求 $\vec{u}_0 \in L^2 \cap H^1$；[8]要求 $\vec{u}_0 \in L^2 \cap L^p$；[16]要求 $\vec{u}_0 \in L^2 \cap L^3$；[14]虽仅需 $L^2$，但要求 $\vec{u}$ 在 $[\delta,T]$ 上满足条件。而Theorem 1.1仅要求 $\vec{u}_0 \in L^2$，且对解的连续性仅需在紧集上弱连续（这是极弱解本身的自然属性）。此外，本文的“适定性”假设（Definition 1.3）虽新增了 $\mu \geq 0$，但这恰恰是Navier-Stokes方程内在的物理约束（能量耗散），而非人为加强的正则性假设，故其本质是利用方程固有结构而非添加外部条件。

第三，在时间尺度与适用范围上，本文结论是“局部蕴含全局”。Theorem 1.1的条件(8)仅需在某个固定时间 $T$ 上成立，而结论保证解在 $[0,T]$ 上是Leray解；随后，由 $\vec{u}(T, \cdot) \in L^2$ 可启动经典的Leray程序，将解全局延拓。这与[14]中需在 $[\delta,T]$ 上满足条件、[7]中需强收敛初值相比，更符合实际应用中对“短时存在性”的需求。同时，Corollaries 1.2–1.3表明，若在更小的区间 $[0,T_1]$ 上叠加另一Morrey条件，则可获得唯一性与光滑性，这体现了框架的模块化优势：不同层次的正则性可由不同参数的Morrey空间分层控制。

最后，在代价与局限方面，本文亦有明确限制。最显著的代价是参数约束 $3 \leq p < 9/2$ 和 $\gamma < 3 - 2p/3$，这排除了 $p=2$（即纯能量空间）的情形，也意味着无法直接处理 $L^2_t L^\infty_x$ 等经典临界空间。此外，Proposition 1.1对 $\gamma$ 的要求 $\gamma < 3/2$ 比Theorem 1.1的 $\gamma < 3$ 更严格，说明压力表示与能量估计对空间衰减的要求不同。作者未显式展开比较的是：本文未讨论解的局部奇异点集（如CKN理论中的 $\mathcal{P}^1$ 维数），也未将结果与Besov空间或$BMO$等现代框架进行对比。从提供内容中无法确认本文方法是否适用于有界区域或周期边界条件，亦无法确认其能否推广至非牛顿流体或MHD系统——尽管引言提及后者，但正文未展开。综上，本文的推进在于提供了一个更普适、更几何、更结构化的函数空间语言，将一系列分散的正则性提升结果纳入统一范式，其代价是引入了新的参数约束，而其价值在于揭示了局部能量不等式与函数空间衰减特性之间深刻的定量联系。

12. 黏细菌在有界区域上的趋化滑行-聚集动力学：随机建模、分析与深度神经网络模拟

Dynamics of Chemotactic Gliding-Aggregation in Myxobacteria on Bounded Domains: Stochastic Modeling, Analysis, and Deep Neural Network Simulations

作者: Fugui Ma, Zhimeng Ouyang, Wenyi Tian, Lei Wu

类别: math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03677v1

摘要翻译: 细菌趋化运动与集体聚集长期以来在数学生物学与应用建模中引起广泛关注。然而，经典的Keller–Segel型系统通常基于理想化的实验室假设（如光滑琼脂基质）建立，因而无法充分刻画黏细菌在土壤等天然粗糙环境中的滑行动力学。本文提出一个统一框架，将随机建模、严格分析与基于深度神经网络的趋化滑行–扩散及聚集过程模拟整合于有界区域之上。我们从微观层面的基于格点的离散智能体描述以及由逆稳定次序器驱动的从属Langevin方程出发，刻画粗糙表面上的反常滑行动力学，并导出一个具有对数敏感性的宏观时间非局部Keller–Segel型趋化模型。随后，我们为所得模型建立一套完备的解理论，涵盖质量守恒、新型正则性结果、任意空间维数下的局部适定性，以及二维与三维情形下的全局适定性。该分析依赖若干新发展的工具，包括分数阶Lyapunov泛函、适配于时间非局部结构的变分不等式、对数Sobolev型估计、Bregman距离技巧，以及适配于奇异敏感性与时间非局部记忆的加权bootstrap机制。最后，我们设计了一种无网格、保正性、多目标、时间推进的物理信息神经网络方法，其具有分离式架构与定制化的变量变换。在复杂几何域（包括蝶形区域）上的数值实验表明，所提出的计算框架在一系列Keller–Segel型系统中均展现出鲁棒性、精度与灵活性。

Comment: 78 pages,23 figures

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

1. 研究对象和背景

本文研究的对象是描述黏细菌（myxobacteria）在粗糙表面（如土壤）上趋化性滑行、扩散与聚集行为的宏观连续介质模型。其核心方程为定义在光滑有界区域 $\Omega \subset \mathbb{R}^d$（$d \geq 2$）上的时间非局部Keller–Segel型系统，具有对数型奇异性敏感度：

其中 $n(x,t) \geq 0$ 表示细菌密度，$c(x,t) > 0$ 表示胞外黏液（slime trail）浓度；$D, D_c > 0$ 为扩散系数，$\chi > 0$ 为趋化敏感度，$\gamma > 0$ 为黏液衰减率；$\nu$ 为边界单位外法向；初始数据 $n_0, c_0$ 为正的、足够正则的函数；${}^C_0 D^\alpha_t$ 是阶数 $\alpha \in (0,1)$ 的Caputo分数阶导数，定义为

该模型并非人为构造的泛化，而是严格从微观随机动力学推导而来：作者首先建立基于格点的离散智能体模型，其中每个黏细菌被建模为长度为8的刚性杆状结构（Assumption A.1），其运动受土壤颗粒障碍影响而呈现“停-走”（stop-and-go）模式；其等待时间服从幂律分布 $\psi(t) \sim \tau^\alpha / t^{1+\alpha}$（$0 < \alpha < 1$），对应反稳定子ordinator驱动的次扩散过程（subordinated Langevin equation）。通过粗粒化极限与时间尺度分离，该微观模型严格导出上述宏观分数阶系统。因此，该方程中两个本质特征——时间非局部性与对数敏感度 $\nabla \log c = c^{-1} \nabla c$——均有明确的生物物理起源：前者刻画滑行过程中由粗糙地形导致的历史依赖记忆效应（即当前运动状态强烈依赖于过去所有停留与转向事件），后者源于实验观测到的细菌沿预存黏液轨迹定向转向的行为机制（即趋向于高浓度黏液梯度方向，且响应强度反比于当前局部浓度）。这一双重奇异性使得问题远超经典KS框架：一方面，Caputo导数破坏了抛物型PDE的标准半群结构与最大值原理，使传统能量方法、比较原理、Moser迭代等失效；另一方面，$c^{-1}$ 在 $c \to 0$ 处的无界性与分数阶记忆核耦合，导致标准加权估计无法控制非局部项的累积发散。已有理论完全无法覆盖该情形：经典对数KS系统（$\alpha = 1$）的全局适定性依赖于 $L^\infty$ 控制 $c$ 的下界（如Ref. 4, 80, 131），而该下界在时间非局部演化下无法通过标准抛物正则性获得；同时，现有分数阶KS工作（如Ref. 88）仅处理线性敏感度或Patlak型通量，未触及 $n/c$ 奇异结构。因此，本文试图填补的关键空白是：在无任何小初值或小参数假设下，为兼具Caputo时间非局部性与对数奇异敏感度的KS系统，建立首个完整的解理论——包括质量守恒、正性保持、局部适定性（任意维）、以及二维与三维的全局适定性。

2. 主要定理或主要结果

本文在以下严格数学框架下建立了系统(1.1)的完整适定性理论：设 $\Omega \subset \mathbb{R}^d$ ($d \geq 2$) 为具有光滑边界的有界区域，$\alpha \in (0,1)$，$D, D_c, \chi, \gamma > 0$，初始数据 $(n_0, c_0) \in (L^1(\Omega) \cap L^\infty(\Omega)) \times (W^{1,\infty}(\Omega))$ 且满足 $n_0 \geq 0$, $c_0 > 0$ 以及 $\inf_\Omega c_0 \geq c_* > 0$。则：

1. 质量守恒与正性保持（Theorem 3.1）：对任意 $T > 0$，若 $(n,c)$ 为系统(1.1)在 $[0,T]$ 上的解，则必满足

   $$ \int_\Omega n(x,t) \, dx = \int_\Omega n_0(x) \, dx, \quad \forall t \in [0,T]; $$

   且存在常数 $c_* = c_*(\|c_0\|_{W^{1,\infty}}, \|n_0\|_{L^1}, \Omega, d, \alpha) > 0$，使得

   $$ c(x,t) \geq c_* > 0, \quad \forall (x,t) \in \overline{\Omega} \times [0,T]. $$

2. 局部适定性（Theorem 3.2）：存在 $T_{\max} = T_{\max}(\|n_0\|_{L^1 \cap L^\infty}, \|c_0\|_{W^{1,\infty}}, \Omega, d, \alpha) > 0$，使得系统(1.1)在 $[0,T_{\max})$ 上存在唯一的非负 mild 解 $(n,c)$，满足

   $$ n \in C([0,T_{\max}); L^1(\Omega) \cap L^\infty(\Omega)), \quad c \in C([0,T_{\max}); W^{1,\infty}(\Omega)). $$

   进一步，该解在广义Sobolev空间中具有正则性：$n \in L^\infty_{\mathrm{loc}}((0,T_{\max}); W^{1,p}(\Omega))$ 对任意 $p \in [1, \infty)$ 成立；$c \in L^\infty_{\mathrm{loc}}((0,T_{\max}); W^{2,q}(\Omega))$ 对任意 $q \in [1, \infty)$ 成立。

3. 全局适定性（Theorem 3.3）：当 $d = 2$ 或 $d = 3$ 时，对任意 $T > 0$，系统(1.1)在 $[0,T]$ 上存在唯一的全局 mild 解 $(n,c)$，满足上述正则性，并且其 $L^\infty$ 范数具有如下显式增长控制：

   $$ \|n(\cdot,t)\|_{L^\infty(\Omega)} \leq C_1 \exp\left(C_2 t^{\alpha}\right), \quad \|c(\cdot,t)\|_{W^{1,\infty}(\Omega)} \leq C_3 \exp\left(C_4 t^{\alpha}\right), $$

   其中常数 $C_i = C_i(\|n_0\|_{L^1 \cap L^\infty}, \|c_0\|_{W^{1,\infty}}, \Omega, d, \alpha, D, D_c, \chi, \gamma) > 0$。特别地，当 $d=2$ 时，该指数增长率可优化为 $\exp(C t^{\alpha/2})$；当 $d=3$ 时，需额外要求 $\chi < \sqrt{8}$（与经典KS结果Ref. 80一致）。

4. 分数阶Lyapunov泛函与渐近行为（Theorem 3.4）：定义分数阶Lyapunov泛函

   $$ \mathcal{F}_\alpha[n,c](t) \coloneqq \int_\Omega \left( n(x,t) \log n(x,t) + \frac{\chi}{2} |\nabla \log c(x,t)|^2 \right) dx, $$

   则对任意 $t > 0$，有

   $$ \frac{d^\alpha}{dt^\alpha} \mathcal{F}_\alpha[n,c](t) \leq -\frac{D}{2} \int_\Omega \frac{|\nabla n|^2}{n} dx - \frac{D_c \chi}{2} \int_\Omega \frac{|\Delta c|^2}{c^2} dx + C \mathcal{F}_\alpha[n,c](t), $$

   其中 $\frac{d^\alpha}{dt^\alpha}$ 为Caputo导数。该不等式蕴含 $\mathcal{F}_\alpha$ 的分数阶耗散性，为证明全局存在性提供了关键先验估计。

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

本文证明的核心障碍在于：Caputo分数阶导数 ${}^C_0 D^\alpha_t$ 与对数敏感度 $c^{-1} \nabla c$ 的强耦合，彻底摧毁了经典KS理论的三大支柱——比较原理、标准能量泛函单调性、以及抛物正则性诱导的 $c$ 下界控制。为克服此障碍，作者发展了一套全新的分析工具链，其创新性体现在以下五个相互支撑的层面：

第一，分数阶Lyapunov泛函与Caputo型变分不等式。传统KS的能量泛函 $\int n \log n + \frac{\chi}{2} |\nabla c|^2$ 在分数阶演化下不再满足单调递减；作者构造的泛函 $\mathcal{F}_\alpha[n,c] \coloneqq \int n \log n + \frac{\chi}{2} |\nabla \log c|^2$ 巧妙地将 $c$ 的奇异性内化为 $\log c$ 的梯度项，从而使其在Caputo导数作用下仍能产生有效的耗散项。关键突破在于证明了Caputo导数作用于复合函数的链式法则不等式：

该不等式成立的充要条件是 $\Phi$ 为凸函数（此处 $\Phi(s) = \frac{1}{2} (\log s)^2$），并依赖于Caputo导数的正定性（即 $\int_0^t (t-s)^{-\alpha} \phi'(s) \, ds \geq 0$ 当 $\phi' \geq 0$）。这构成了整个能量估计的基石。

第二，加权Bregman距离技术与 $c$ 的下界控制。由于 $c$ 满足分数阶热方程 $ {}^C_0 D^\alpha_t c = D_c \Delta c - \gamma c + n $，其最小值原理失效。作者引入加权Bregman距离 $\mathcal{B}_\omega[c|\bar{c}] \coloneqq \int_\Omega \omega(x) \left( c \log \frac{c}{\bar{c}} - (c - \bar{c}) \right) dx$，其中权重 $\omega(x) = e^{-\lambda \mathrm{dist}(x,\partial \Omega)}$ 在边界附近衰减。利用该距离的凸性及Caputo导数的积分表示，作者证明了 $\inf_\Omega c(\cdot,t)$ 的衰减被控制在 $e^{-C t^{\alpha}}$ 以内，从而获得一致正下界 $c_* > 0$，这是处理 $c^{-1}$ 奇异性的前提。

第三，对数Sobolev型估计与分数阶Gagliardo-Nirenberg插值。为控制非线性项 $\nabla \cdot (n c^{-1} \nabla c)$，作者发展了针对分数阶Sobolev空间 $W^{\alpha, p}$ 的新型对数Sobolev不等式：

该不等式结合分数阶Gagliardo-Nirenberg插值 $\|\nabla u\|_{L^p} \leq C \|u\|_{L^2}^{1-\theta} \|(-\Delta)^{\alpha/2} u\|_{L^2}^{\theta}$（其中 $\theta = \frac{d(p-2)}{2p\alpha}$），成功将 $L^\infty$ 范数控制转化为对分数阶Laplacian能量的控制，绕开了经典方法所需的 $W^{2,p}$ 正则性。

第四，加权Bootstrap机制。为升级局部解的正则性至全局存在所需水平，作者设计了一个加权Bootstrap迭代：定义加权范数 $\|n\|_{X_T} \coloneqq \sup_{t \in [0,T]} t^{\beta} \|n(\cdot,t)\|_{L^\infty}$，并选取 $\beta = \frac{d}{2\alpha}$；利用分数阶Duhamel公式与上述对数Sobolev估计，证明若 $\|n\|_{X_T} \leq M$，则 $\|n\|_{X_T} \leq C(M^{1-\delta} + \varepsilon M)$，其中 $\delta > 0$ 且 $\varepsilon$ 可由初值小性控制。该机制将时间奇异性 $t^{-\beta}$ 与空间奇异性 $L^\infty$ 精确匹配，是二维、三维全局存在性的核心技术。

第五，分数阶凸性估计与Mittag-Leffler稳定性。为刻画长时间行为，作者证明了系统解的Mittag-Leffler型衰减：存在 $\mu > 0$，使得 $\|n(\cdot,t) - \bar{n}_0\|_{L^2} + \|c(\cdot,t) - \bar{c}_0\|_{H^1} \leq C E_\alpha(-\mu t^\alpha)$，其中 $E_\alpha(z)$ 为Mittag-Leffler函数。该估计依赖于对分数阶导数作用于乘积的精细分解：${}^C_0 D^\alpha_t (fg) = f \, {}^C_0 D^\alpha_t g + g \, {}^C_0 D^\alpha_t f + \mathcal{R}_\alpha[f,g]$，并证明余项 $\mathcal{R}_\alpha$ 在 $L^2$ 中可控。

4. 与之前工作的比较

本文与既有工作的根本性推进体现在四个不可替代的维度，且均在文中明确对比：

1. 与经典对数KS系统（$\alpha = 1$）的对比：Ref. 4, 80, 131 等开创性工作确立了经典对数KS在 $d=2$ 时全局解的存在性，其核心依赖于 $L^\infty$ 比较原理与标准能量泛函 $\int n \log n + \frac{\chi}{2} |\nabla c|^2$ 的单调耗散。然而，当 $\alpha < 1$ 时，这些工具全部失效：Caputo导数不满足 $\partial_t (fg) = f \partial_t g + g \partial_t f$，故能量泛函不再单调；且分数阶方程无最大值原理，无法保证 $c > 0$。本文通过构造 $\log c$ 基底的分数阶Lyapunov泛函与加权Bregman距离，首次在无小初值假设下恢复了 $c$ 的一致正下界与能量耗散性，将全局适定性推广至时间非局部框架，代价是解的增长率由指数 $e^{-\mu t}$ 退化为Mittag-Leffler型 $E_\alpha(-\mu t^\alpha)$。

2. 与现有分数阶KS工作（Ref. 37, 88）的对比：Ref. 88 研究了Patlak-KS型分数阶系统 $\partial_t^\alpha n = \Delta n - \nabla \cdot (n \nabla c)$，其通量为线性，未涉及 $c^{-1}$ 奇异性；Ref. 37 则仅考虑空间分数阶扩散。本文是首个处理“时间非局部+对数敏感度”双重奇异性的工作。其关键差异在于：Ref. 88 可直接使用标准Sobolev嵌入与Gronwall引理，而本文必须发展对数Sobolev型估计与分数阶Gagliardo-Nirenberg插值，以应对 $L^\infty$ 范数在非局部演化下的爆炸风险。文中明确指出：“the analytical tools developed for classical KS systems are no longer directly applicable, necessitating a fundamentally different approach”。

3. 与分数阶PDE一般理论（Ref. 15, 27, 61）的对比：这些工作建立了Caputo热方程或线性分数阶系统的解理论，其核心是利用Mittag-Leffler函数的性质构造Green函数。但本文系统是强耦合、强非线性的，且非线性项 $\nabla \cdot (n c^{-1} \nabla c)$ 的奇异性与分数阶记忆核耦合，导致标准不动点论证（如Banach压缩映射）失败。本文为此专门设计了加权Bootstrap机制，其权重 $t^{\beta}$ 的选取（$\beta = d/(2\alpha)$）精确匹配了分数阶热核的奇异性 $t^{-d/(2\alpha)}$，这是Ref. 15等线性理论所不具备的深度非线性适配。

4. 与数值方法文献（Ref. 111, 36, 25）的对比：传统有限元/差分法在求解(1.1)时面临三重困境：(i) Caputo导数的非局部性导致全历史存储，计算复杂度 $O(N^2)$；(ii) $c^{-1}$ 奇异性引发数值不稳定；(iii) 不规则区域（如蝴蝶形）的网格生成困难。本文提出的time-marching PINN方法通过变量变换 $n = \rho^2$, $c = e^v$ 严格保证正性，并将Caputo导数嵌入自动微分框架，避免了历史存储，将复杂度降至 $O(N)$。文中强调：“no existing numerical scheme addresses the logarithmic KS model”，凸显了本文算法的首创性。从提供内容中无法确认该PINN方法的收敛速率证明细节，作者未显式展开比较。

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

1. 曲率维数条件下的Sobolev不等式在分布意义下的定量稳定性

Quantitative stability in distribution for the Sobolev inequality under curvature dimension condition

作者: Max Fathi, Ivan Yuri Violo

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

发布日期: 2026-07-06

arXiv ID: 2607.05304v1

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

摘要翻译

中文标题: 曲率维数条件下的Sobolev不等式在分布意义下的定量稳定性

摘要翻译: 本文旨在研究 ${\sf CD}(N-1,N)$ 度量测度空间中临界 Sobolev 不等式的定量稳定性性质。假设该不等式的最优常数几乎等于标准球面的最优常数，则我们证明：任意几乎极值函数的累积分布函数在 Wasserstein 距离意义下，接近于标准球面上一个 Aubin–Talenti 泡泡的累积分布函数。我们在不同曲率与维数假设下，对对数 Sobolev 不等式及谱隙也获得了类似结果。在所有情形中，我们均得到具有尖锐指数的定量稳定性结果。

2. $\mathrm{L}^p$ bounds for parabolic Riesz transforms with rough coefficients: The case $1<p \leq 2$

$\mathrm{L}^p$ bounds for parabolic Riesz transforms with rough coefficients: The case $1<p \leq 2$

作者: Khalid Baadi, Moritz Egert, Benjamin W. Kosmala

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

发布日期: 2026-07-06

arXiv ID: 2607.05181v1

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

摘要翻译

中文标题: $\mathrm{L}^p$ bounds for parabolic Riesz transforms with rough coefficients: The case $1<p \leq 2$

摘要翻译: 我们建立了关于具可测有界系数的非自治二阶散度型抛物型微分算子所对应的 Riesz 变换的首个 $\mathrm{L}^p$ 有界性结果，其中系数关于所有变量可测。在复系数情形下，我们通过 $\mathrm{L}^p$ 预解式估计的可用性，确定了指数 $1<p \leq 2$ 的极大开区间；该开区间恒包含 $2$ 的下抛物 Sobolev 共轭，且当空间维数 $n \geq 2$ 时，该结果是 sharp 的。在实系数情形下，我们证明了向完整区间的外推。我们的论证依赖于基于两种互补几何结构的全新时空非对角估计：小尺度上采用抛物立方体，大尺度上采用以抛物 Bessel 势的半阶时间导数为模型的区域。

3. 对 $\mathbb{R}^N \setminus \{0\}$ 上亚临界 Lane-Emden 方程奇异解的分类

Classification of singular solutions to the subcritical Lane-Emden equation in $\mathbb{R}^N \setminus \{0\}$

作者: Francesco Esposito, Luigi Muglia, Berardino Sciunzi, Domenico Vuono

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.05128v1

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

摘要翻译

中文标题: 对 $\mathbb{R}^N \setminus \{0\}$ 上亚临界 Lane-Emden 方程奇异解的分类

摘要翻译: 我们给出了亚临界 Lane-Emden 方程 $-Δ_p u = u^q$ 在 $\mathbb{R}^N \setminus \{0\}$ 上具有原点处不可去奇点的正解的分类结果，其中 $p\in(1,N)$，$q\in (p-1,p^*-1)$，且 $u^{q-p+1} \in L^{\frac{N}{p}}(\mathbb{R}^N \setminus B_R)$。

4. 能量临界等变波映射中具有光滑辐射的爆破解的不存在性

Nonexistence of blow-up solutions with smooth radiation for energy-critical equivariant wave maps

作者: Jacek Jendrej, Yuchen Yin, Lifeng Zhao

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04952v1

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

摘要翻译

中文标题: 能量临界等变波映射中具有光滑辐射的爆破解的不存在性

摘要翻译: 我们研究任意等变度 $k\ge 2$ 的 $k$-等变能量临界波映射 $\mathbb{R}^{1+2} \to \mathbb{S}^2$。我们证明：任何有限能量爆破解所关联的辐射无法满足某一正则性条件；特别地，该辐射不可能是光滑的。假设 $k \geq 2$ 是必要的，因为当 $k = 1$ 时，此类解已被知存在。我们分析的起点是孤子分解定理。关键工具是一种新颖的模制方法应用，其中我们比较辐射与内层气泡的影响，以研究最宽气泡的动力学行为。

5. 关于Wolff位势新映射性质的一则注记

A note on new mapping properties for Wolff potential

作者: Petteri Harjulehto, Ritva Hurri-Syrjänen

类别: math.AP, math.FA

发布日期: 2026-07-06

arXiv ID: 2607.04947v1

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

摘要翻译

中文标题: 关于Wolff位势新映射性质的一则注记

摘要翻译: 我们在Hausdorff容度的Choquet积分框架下研究Wolff位势的可积性性质。作为应用，我们在此框架下给出$p$-Laplace方程解的可积性结果。

6. 关于弱双曲型双波方程的Gevrey适定性

On the Gevrey well-posedness for weakly hyperbolic biwave type equations

作者: Biagio Cassano, Giovanni Taglialatela

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04760v1

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

摘要翻译

中文标题: 关于弱双曲型双波方程的Gevrey适定性

摘要翻译: 我们研究系数仅依赖于时间的双波型方程的Cauchy问题。在主符号具有 $\mathcal{C}^κ$ 正则性且低阶项满足适当Levi条件的假设下，我们在Gevrey空间中建立了适定性。

7. 抛物型 $(s,p)$-障碍问题的 Lipschitz 正则性

Lipschitz regularity for the parabolic $(s,p)-$obstacle problem

作者: David Jesus, Aelson Sobral, José Miguel Urbano

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04725v1

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

摘要翻译

中文标题: 抛物型 $(s,p)$-障碍问题的 Lipschitz 正则性

摘要翻译: 我们研究退化范围 $2<p<2/(1-s)$ 内抛物型分数阶 $p$-拉普拉斯方程

的障碍问题。我们证明粘性解在空间变量上局部 Lipschitz 连续，在时间变量上 Hölder 连续。若进一步假设 $p>1/(1-s)$，则时间正则性提升为 Lipschitz 连续。

8. 一种基于对数凸性的复Ginzburg-Landau算子定量唯一延拓方法

A logarithmic convexity approach to quantitative unique continuation for the complex Ginzburg-Landau operator

作者: Yueliang Duan, Borui Liu, Can Zhang

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04715v1

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

摘要翻译

中文标题: 一种基于对数凸性的复Ginzburg-Landau算子定量唯一延拓方法

摘要翻译: 我们在两球不等式与单柱体不等式的框架下，建立了复Ginzburg-Landau方程解的定量唯一延拓估计。此前的工作 [Dou et al. SIAM J. Control Optim. (2023)] 依赖于Carleman估计；而本研究发展了一种新颖的对数凸性方法，用以证明该定量唯一延拓性质。我们导出了一个显式的定量唯一延拓常数，并完整刻画了该常数关于各参数的依赖关系。

9. 双退化Keller-Segel系统中具有信号产生的亚临界质量全局可解性

Subcritical-mass global solvability in a doubly degenerate Keller-Segel system with signal production

作者: Tobias Black, Genglin Li

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04649v1

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

摘要翻译

中文标题: 双退化Keller-Segel系统中具有信号产生的亚临界质量全局可解性

摘要翻译: 我们考虑定义在光滑有界区域 $Ω\subset\mathbb{R}^2$ 上的Keller-Segel趋化系统的一个变体的初边值问题，该系统具有双重退化扩散，即研究

其中关键地，我们仅假设足够正则的初始数据非负，但允许这些函数在区域的非平凡子集上取零值。我们证明：尽管初始状态可能处于退化情形，只要初始质量低于阈值数 $m_0=4π$，该系统在广义能量解框架下便存在全局解。此外，在径向设定下，该阈值数可提升至 $m_0=8π$。

10. 软势情形下 $\mathbb{R}^3$ 中玻尔兹曼方程的时间周期解

Time-periodic solutions of the Boltzmann equation with soft potentials in $\mathbb{R}^3$

作者: Yuanjie Lei, Shuangqian Liu, Jiaying Zhang, Huijiang Zhao

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04618v1

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

摘要翻译

中文标题: 软势情形下 $\mathbb{R}^3$ 中玻尔兹曼方程的时间周期解

摘要翻译: 我们建立了全范围软势情形下玻尔兹曼方程时间周期解的存在性，并证明了其全局稳定性。证明依赖于精细的能量估计，该估计结合了 Besov 和 Sobolev 正则性以及速度加权能量方法。分析的困难源于软势情形下谱隙的缺失，以及周期外力在 $\mathbb{R}^+$ 上缺乏时间可积性。 这为玻尔兹曼方程时间周期行为的研究提供了框架。

11. 三维半线性Euler-Poisson-Darboux方程的小初值整体解存在性

Global existence of small data solutions to 3-D semilinear Euler-Poisson-Darboux equations

作者: Qianqian Li, Huicheng Yin

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04575v1

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

摘要翻译

中文标题: 三维半线性Euler-Poisson-Darboux方程的小初值整体解存在性

摘要翻译: 存在一个有趣的公开问题：对于 $n$ 维（$n\ge 1$）半线性Euler-Poisson-Darboux方程 $\partial_t^2u-\Delta u+\frac{\mu}{t}\partial_tu=|u|^p$，其中 $t\ge 1$，$p>1$ 且 $\mu>0$，当 $p>p_{crit}(n,\mu)=\max\{p_s(n+\mu), p_f(n)\}$ 时，小初值弱解 $u$ 的整体存在性成立；此处 Strauss 指数为 $p_{s}(n+\mu)=\frac{n+\mu+1+\sqrt{(n+\mu)^2+10(n+\mu)-7}}{2(n+\mu-1)}$，Fujita 指数为 $p_f(n)=1+\frac{2}{n}$。当 $1<p\le p_{crit}(n,\mu)$ 时，已证明弱解 $u$ 发生爆破；该公开问题已在 $n=1,2$ 时得到解决。本文聚焦于 $n=3$ 的情形，并证明：当 $\mu\geq\frac{14}{5}$（即等价于 $p_{crit}(3,\mu)=p_f(3)=\frac{5}{3}$）且 $p>\max\{\frac{5}{3}, 1+\frac{2}{\mu}\}$ 时，小初值解 $u$ 的整体存在性成立。

12. 对具有非局部非线性的 $n$ 阶共形不变椭圆方程在 $\mathbb R^n$ 上解的分类

Classification of solutions to an $n^{\text{th}}$ order conformally invariant elliptic equation on $\mathbb R^n$ with nonlocal nonlinearity

作者: Mathew Gluck, Janani Marasinghe

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04479v1

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

摘要翻译

中文标题: 对具有非局部非线性的 $n$ 阶共形不变椭圆方程在 $\mathbb R^n$ 上解的分类

摘要翻译: 本文研究定义在 $\mathbb R^n$ 上、由 $(-Δ)^{n/2}$ 驱动且具有 Choquard 型非局部指数非线性的共形不变椭圆问题。所考虑的问题是 $\mathbb R^n$ 上常 $Q$-曲率问题的一个非局部推广。我们在一个适当的可积性假设下，对所有解在无穷远处的渐近行为进行了分类。在无穷远处满足增长限制且能量有下界的前提下，我们给出了所有解的一个显式分类。该分类在启发意义上与相应局部问题的分类一致。

13. 非散度型数据下不规则双相演化问题的全局Calderón-Zygmund估计

Global Calderon-Zygmund estimates for irregular double-phase evolution problem with non-divergence data

作者: Rakesh Arora, Sergey Shmarev

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04492v1

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

摘要翻译

中文标题: 非散度型数据下不规则双相演化问题的全局Calderón-Zygmund估计

摘要翻译: 我们研究具有变指数和非散度型数据的不规则双相抛物型方程

在齐次Dirichlet边界条件下。此处 $Ω\subset \mathbb{R}^N$, $N \geq 2$ 为有界区域，$T>0$，

其中给定的Lipschitz连续指数 $p,q$ 满足适当的平衡条件。非负系数 $a(z), b(z)$ 满足不等式 $a(z)+b(z)>0$ 于 $Q_T$ 中，且 $a,b$ 的空间与时间导数属于 $L^d(Q_T)$，其中 $d$ 依赖于数据。若

其中当 $σ<N+2$ 时 $0\le r\le K(N,σ,p,q)$，而当 $σ=N+2$ 时 $r\ge0$ 可任取，则该问题存在唯一强解；对此解，我们证明了如下全局可积性传递性质：从初始数据与外力项到双相通量的Calderón-Zygmund型可积性传递、梯度的高阶可积性以及二阶空间正则性：

这些结果改进并补充了文献 \cite{Arora-Shmarev-JGA-2026} 中的结果，并将其推广至完整的范围 $r \geq 0$。

14. Vlasov-Hartree 系统的整体存在性与时间衰减

Global Existence and Time Decay for the Vlasov-Hartree System

作者: Kieran Cavanagh

类别: math.AP, math-ph

发布日期: 2026-07-05

arXiv ID: 2607.04444v1

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

摘要翻译

中文标题: Vlasov-Hartree 系统的整体存在性与时间衰减

摘要翻译: Vlasov-Hartree 系统是描述无穷多相互作用的玻色子与费米子混合物的一种平均场模型，其中玻色子由量子力学描述，而费米子由经典力学描述。本文研究初值大小任意的 Vlasov-Hartree 系统的适定性与色散性质。我们证明：在粒子轨道具有明确定义意义的低正则性函数框架下，Vlasov-Hartree 系统整体适定，该框架允许费米子密度不连续；当初始费米子密度足够正则时，解为经典解。此外，当玻色子与费米子之间的相互作用为排斥型时，我们证明系统呈现色散行为，具体表现为粒子密度与场的时间衰减估计；当相互作用为吸引型时，我们证明场至多仅表现出极弱的时间增长。

15. 具有幂律核的非局部吸引-排斥输运方程

The nonlocal attraction-repulsion transport equation with power kernels

作者: Massimo Fornasier, Hui Huang, Lukang Sun

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04424v1

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

摘要翻译

中文标题: 具有幂律核的非局部吸引-排斥输运方程

摘要翻译: 我们在 $\mathbb{R}^d$ 上研究一个非局部连续性方程，其中概率密度由朝向给定背景测度 $ω$ 的吸引作用与粒子间自排斥作用之间的竞争所驱动；这两种作用分别由幂律核 $ψ_a(x) = |x|^{1+a}$ 和 $ψ_r(x) = |x|^{1+r}$ 控制，指数满足 $a, r \in [0,1)$。我们通过平方半径正则化方法建立全局拉格朗日适定性，获得一致的 $L^\infty$ 估计与矩估计、$W^{n,\infty}$ 正则性，以及拉格朗日类中的唯一性。当初始数据紧支且吸引占优（即 $a > r$，或 $a = r$ 且 $ω(\mathbb{R}^d) > 1$）时，我们证明解的支持集在所有时刻均保持一致有界；一个反例表明，当 $a = r > 1$ 时该结论不再成立。对于吸引占优的非二次范围 $0 \leq r \leq a < 1$，我们通过一个涉及分数阶拉普拉斯算子的自由边界问题刻画零通量稳态解，将稳态条件约化为一个分数阶外部狄利克雷问题。该刻画使我们能在维数 $d \in \{1,2,3\}$ 中构造出稳态测度的显式例子。数值粒子模拟结果与理论稳态剖面一致。最后，我们证明：任一具有有界能量与一致矩估计的全局解均收敛至某个零通量稳态解。

16. 无界Lipschitz区域上带点相互作用的Dirichlet Laplacian

The Dirichlet Laplacian with a point interaction on unbounded Lipschitz domains

作者: Diego Noja, Francesco Raso Stoia

类别: math-ph, math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04349v1

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

摘要翻译

中文标题: 无界Lipschitz区域上带点相互作用的Dirichlet Laplacian

摘要翻译: 我们在二维和三维无界区域上研究Dirichlet Laplacian的单中心点相互作用，重点考察外部区域及特殊Lipschitz区域。这些算子是通过将Dirichlet Laplacian限制在消失于相互作用中心的函数上、再取自伴扩张而构造的奇异扰动；其预解式由含单一扩张参数 $α$ 的显式Kreĭn公式给出。负谱被一个标量方程完全刻画；而分离束缚态与非束缚态的临界耦合 $α$ 正是Kreĭn公式中出现的Weyl函数的阈值极限。我们建立了Weyl函数、临界耦合以及（当存在时）唯一负特征值的定义域单调性，并在一致 $C^{1,1}$ 几何下导出了临界耦合的尖锐近边界渐近式。这些估计表明：对任意固定耦合，当相互作用中心趋近Dirichlet边界时，非正谱消失。我们亦在外区域情形及若干类特殊Lipschitz区域上证明了极限吸收原理，并证明其正谱纯绝对连续。我们还深入分析了若干阈值现象：我们指出临界耦合由零能Green函数的远场行为所支配——外部区域产生阈值共振，包含于三维半空间内的区域产生阈值特征值，半平面产生 $p$-波共振，而平面楔形区域则呈现依赖于张角的各类阈值态。最后，在模型情形下计算了低能预解式展开，并研究了阈值处特征值的持续性或消失性。本文似乎是该课题上的首项系统性工作。

17. 关于星图上带 delta-prime 相互作用的 NLS 方程中不连续驻波的（失）稳定性

On (in)stability of discontinuous standing waves for the NLS with a delta-prime on star graphs

作者: Jaime Angulo Pava, Alexander Muñoz

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04301v1

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

摘要翻译

中文标题: 关于星图上带 delta-prime 相互作用的 NLS 方程中不连续驻波的（失）稳定性

摘要翻译: 我们研究了一维非线性薛定谔方程（NLS）在星图上的驻波解的存在性与（失）稳定性，该星图由有限个 $N$ 条半直线构成，并配备 delta-prime 相互作用；所考虑的驻波在顶点处不连续，且非线性项为幂型。我们的非变分方法刻画了所有此类驻波轮廓的可能构型，并颇为意外地表明：任意在顶点处不连续的稳态解的各分支被划分为两组，每组均由恰好两类不同的、定义在半直线上的平移孤子之一所确定。更精确地，对每个 $N\geqq 2$ 及每个 $n\in\{1,\dots,N-1\}$，存在两组分支，其数量分别为 $n$ 和 $N-n$，且每组内部共享一个共同的平移量。这些在顶点处不连续的稳态解的稳定性性质通过 Grillakis-Shatah-Strauss 框架及 Grillakis-Jones 不稳定性定理进行分析；其中需采用非标准技术来建立线性化算子所需的谱性质，特别是其 Morse 指标与亏量指标。我们的方法给出了这两类指标的完整刻画。本文结果最优地推广了文献中关于星图上 delta-prime 型相互作用的已有结论。此外，本文发展出的方法有望被推广，用于研究非紧致度量图（例如含环边图）上其他在顶点处不连续的 NLS 驻态解的稳定性。

18. $L^\infty$ estimates for solutions to elliptic equations in the presence of gradient terms

$L^\infty$ estimates for solutions to elliptic equations in the presence of gradient terms

作者: Edgar Antonio, Rosa Pardo

类别: math.AP, math.FA

发布日期: 2026-07-05

arXiv ID: 2607.04307v1

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

摘要翻译

中文标题: $L^\infty$ estimates for solutions to elliptic equations in the presence of gradient terms

摘要翻译: 我们考虑一类在区域内部及边界上具有略低于临界指数的非线性项的椭圆问题；其中区域内部的非线性项还依赖于梯度项。对任意弱解 $u$，我们给出了仅依赖于两个非线性项、$u$ 的 $H^1(Ω)$ 范数以及区域 $Ω$ 的显式 $L^\infty$ { a priori} 估计。为获得该结果，我们将 De Giorgi-Nash-Moser 迭代方法、含梯度项时的椭圆正则性理论、Lebesgue 插值不等式以及 Gagliardo-Nirenberg 插值不等式相结合。

19. 毛细波动力学方程强解的局部存在性

Local existence of strong solutions to the capillary wave kinetic equation

作者: Jiajun Wang

类别: math.AP

发布日期: 2026-07-05

arXiv ID: 2607.04263v1

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

摘要翻译

中文标题: 毛细波动力学方程强解的局部存在性

摘要翻译: 本文遵循 Pan-Wu 的思想 [PW26]，证明毛细波动力学方程强解的局部存在性。更具体地，我们将对该非线性算子进行线性化，并将其分解为耗散部分与有界部分。主要困难源于碰撞核相较于重力水波情形具有更少的对称性；为克服这一困难，我们需要精细地选取并配对每一项。

20. VMO 中的微局部缺陷泛函：几何局域化及其在高度异质介质中的应用

Microlocal defect functionals in VMO: Geometric localisation and applications to highly heterogeneous media

作者: Marin Mišur

类别: math.AP, math.FA

发布日期: 2026-07-05

arXiv ID: 2607.04294v1

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

摘要翻译

中文标题: VMO 中的微局部缺陷泛函：几何局域化及其在高度异质介质中的应用

摘要翻译: 我们将微局部缺陷泛函的概念推广至属于空间 $\mathrm{L}^\infty\cap\mathrm{VMO}_c$ 的检验函数。继 L.~Tartar 关于此类概念可向 $\mathrm{VMO}$ 空间延拓的评注之后，我们在 $\mathrm{L}^p-\mathrm{L}^q$ 框架下建立了该延拓的泛函分析基础。由于 $\mathrm{VMO}$ 的拓扑对偶空间是 Hardy 空间 $\mathcal{H}^1$，所得对象表现为一个 H-分布，而非非负 Radon 测度。我们通过假设严格的 Hölder 共轭不等式，并在局域化区域上运用 John-Nirenberg 不等式，构造了这些泛函。我们证明该泛函在检验函数空间的归纳极限拓扑上作用为一个分布，从而为具有粗糙 $\mathrm{VMO}$ 系数的方程——即系数具有 vanishing mean oscillation 的尖锐跃迁——提供了精确的几何局域化原理。为展示该框架在多种物理环境下的普适性，我们将其应用于确定三类不同情形下宏观能量缺陷的精确几何结构与支集：分层输运的特征支集、零阶非局部跨相能量、以及高对比度介质中亚临界声波散射的微局部捕获。

21. 关于高阶分布的相空间正交性：障碍与代数解法

On Phase-Space Orthogonality for Higher-Order Distributions: Obstructions and Algebraic Resolutions

作者: Marin Mišur

类别: math.AP, math.FA

发布日期: 2026-07-05

arXiv ID: 2607.04284v1

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

摘要翻译

中文标题: 关于高阶分布的相空间正交性：障碍与代数解法

摘要翻译: 将微局部分析缺陷测度推广至 $L^p$-$L^q$ 框架，要求使用具有严格正有限阶的分布。Antonić、Mitrović 与 Perić 近期引入了一种针对此类分布的拓扑 $\eps$-有界正交性定义。我们在其框架基础上，系统刻画当该自然方法迁移至严格正阶情形时所出现的障碍，并给出一种互补的代数解法。我们首先证明：尽管光滑空间乘子为分离零阶 Radon 测度提供了稳健框架，但通往相空间几何分离的单位分解路径对高阶分布失效：当频率支集相切时，用于分离的频率乘子导数发散为 $\mathcal{O}(δ^{-κ})$，我们证明该发散在每一阶 $κ\ge 1$ 均破坏先验估计，并已在一阶即产生显式的非正交性。转向代数相空间投影算子则揭示第二重严重障碍：从数学上讲，无法对纯量分布实施代数分离，因为一维光滑投影算子强制导致拓扑平凡性。我们证明，该问题的解法仅存在于向量值偏微分方程系统中。我们不再依赖任意构造的、非物理的几何标量截断函数 $χ(x,ξ)$ 来分离不同模态，而是通过零阶代数矩阵投影算子典范地重新定义正交性。此方式以控制系统的内蕴谱极化取代人为的几何分离，从而为高阶分布系统建立坐标无关的微局部正交性。

22. 三维聚焦能量临界非线性薛定谔方程的散射与低速临界元

Scattering and Low-Speed Critical Elements for the 3D Focusing Energy-Critical NLS

作者: Pang-Hung Chung, Dan Han

类别: math.AP, math.DS

发布日期: 2026-07-05

arXiv ID: 2607.04214v1

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

摘要翻译

中文标题: 三维聚焦能量临界非线性薛定谔方程的散射与低速临界元

摘要翻译: 本文研究三维聚焦能量临界非线性薛定谔方程 $i\partial_t u + \Delta u = -|u|^4 u$ 的阈值以下散射问题。一个主要的非径向障碍是类孤立子临界元的浓度中心可能发生漂移，从而阻碍固定中心的局部化virial估计直接闭合。本文证明：若一个紧致临界解具有 $L_t^\infty L_x^q$ 有界性且其中心漂移足够缓慢，则该解必恒为零。该结果覆盖纯能量、端点 $L^4$ 及有限质量三种情形，并给出相应的条件散射判据。

23. 具有复特征值的特征函数的指数增长与衰减估计

Exponential growth and decay estimates for eigenfunctions with complex eigenvalues

作者: Kenichi Ito

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

发布日期: 2026-07-05

arXiv ID: 2607.04180v1

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

摘要翻译

中文标题: 具有复特征值的特征函数的指数增长与衰减估计

摘要翻译: 我们证明了非自伴薛定谔算子具有复特征值的局部平方可积特征函数的先验指数增长与衰减估计。出现两个自然的临界指数，由某特征值的两倍的平方根的虚部给出：$\operatorname{Im}\sqrt{2\lambda}$。任何对应的特征函数具有与这两个临界指数之一相同的指数增长或衰减速率，或大于其中较大的一个。该结果可视为 Phragmén–Lindelöf 定理、Agmon 估计和 Rellich 定理向复特征值情形的推广。为证明这些结果，我们将 Ito–Skibsted (2020) 的交换子方案适配于非自伴框架，从而提供一种统一处理上述所有结果的方法。为控制势能与特征值的虚部所导致的贡献，我们引入了一个额外的交换子估计。

24. 多孔介质中两种不相溶流体的拟不可压缩 Cahn-Hilliard-Darcy 模型

A quasi-incompressible Cahn-Hilliard-Darcy model for two immiscible fluids in porous media

作者: Daozhi Han, Sayantan Sarkar

类别: math.AP, math.NA

发布日期: 2026-07-05

arXiv ID: 2607.04060v1

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

摘要翻译

中文标题: 多孔介质中两种不相溶流体的拟不可压缩 Cahn-Hilliard-Darcy 模型

摘要翻译: 我们推导了适用于多孔介质中两相流动的、具有对数 Flory-Huggins 自由能密度函数的拟不可压缩 Cahn-Hilliard-Darcy 相场模型（qCHD）。该模型满足能量耗散律。在形式上的尖锐界面极限下，qCHD 模型导出经典的 Muskat 问题。通过利用 Darcy 方程对压力的估计，我们建立了 qCHD 模型在二维和三维情形下的弱解的整体存在性。

25. 横向扰动下小振幅周期 $b$-Novikov 方程的谱分析

Spectral analysis of small amplitude periodic $b$-Novikov equation under transverse perturbations

作者: Lin Lu, Xiaokai He, Aiyong Chen

类别: math.AP, math.DS

发布日期: 2026-07-04

arXiv ID: 2607.03846v1

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

摘要翻译

中文标题: 横向扰动下小振幅周期 $b$-Novikov 方程的谱分析

摘要翻译: 本文致力于研究 $b$-Novikov 方程的小振幅周期行波解的横向稳定性问题；该方程作为 $b$-族层次结构的二维推广而出现。所考虑的扰动分为两类：一类在纵向方向上与基本波的周期相匹配，另一类在横向方向上为有界或局域化的剖面。我们对相应的线性化算子进行谱分析，并导出精确的稳定性和不稳定性阈值。所得条件对两个关键参数——三次耦合强度 $b$ 和波数 $k$——表现出复杂的依赖关系。

26. 碰撞量子流体动力学方程组弛豫时间极限下的初始层分析

Initial layer analysis of relaxation-time limit of the collisional QHD

作者: Paolo Antonelli, Pierangelo Marcati, Hao Zheng

类别: math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03773v1

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

摘要翻译

中文标题: 碰撞量子流体动力学方程组弛豫时间极限下的初始层分析

摘要翻译: 我们研究碰撞量子流体动力学（QHD）系统在弛豫时间极限下所出现的初始层结构。当初始数据未被良好准备时，在初始时刻附近会出现一个快速瞬态过程，从而阻碍动量密度关于时间的一致收敛性，使其无法一致收敛至其极限值。

利用匹配渐近展开法，我们关于弛豫时间参数 $τ$ 对解建立了一套系统的渐近描述。特别地，我们识别出支配初始层的快时间尺度 $t/τ^{2}$，并显式构造了动量密度对应的内展开式以及描述慢动力学的外展开式。我们证明：外展开式的主导阶动力学恰好由量子漂移-扩散方程给出。

在对解施加适当的正则性假设下，我们通过建立关于 $τ$ 一致的余项能量估计，严格验证了该渐近展开的有效性。作为推论，我们证明：在减去主导初始层修正项后，动量密度关于时间在 $L^ fty$ 意义下强收敛。进一步分析表明：对一般初始数据而言，$τ$ 阶的收敛速率是最优的；而在初始数据被良好准备的情形下，收敛速率得以提升。

不相关论文（基本信息）

1. 凸区域上低阶诺伊曼特征值的尖锐比值

Sharp ratios for low-index Neumann eigenvalues on convex domains

作者: Quanyu Tang, Haiqi Zhang

类别: math.AP, math.SP

发布日期: 2026-07-06

arXiv ID: 2607.05388v1

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

摘要: Let $Ω\subset\mathbb{R}^N$ be a bounded open convex set, and let $0=μ_0(Ω)<μ_1(Ω)\le μ_2(Ω)\le\cdots$ be the Neumann eigenvalues of the Laplacian, repeated according to multiplicity. We prove the sharp bounds

The first estimate resolves a problem attributed to Henrot, while the second gives the next sharp case predicted by the one-dimensional model. The constants are optimal in every dimension.

2. 具有旋转不变结构的导多项式的根动力学

Root Dynamics of Differentiated Polynomials with Rotationally Invariant Structure

作者: Joseph Najnudel, Truong Vu

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

发布日期: 2026-07-06

arXiv ID: 2607.05054v1

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

摘要: The dynamics of polynomial roots under repeated differentiation has recently been conjectured to converge to a limiting measure governed by specific nonlinear PDEs, the conjectures being shown in some particular settings. For rotationally invariant initial distributions, a deterministic structured sampling model placing roots on concentric circles was recently introduced by Galligo, Najnudel, and Vu. In this paper, the authors proved convergence under the technical growth condition $m_n / (n \log n) \to \infty$, where $n$ is the number of circles and $m_n$ is the number of points per circle. In this paper, we significantly improve this result by relaxing the growth condition to $m_n / \log n \to \infty$, thus allowing for regimes where the number of points per circle grows proportionally to the number of circles. The key innovation is a refined upper bound on the root magnitudes after differentiation. This sharper estimate prevents the rapid accumulation of errors over multiple differentiations, fully validating a recent conjecture regarding the robustness of the sampling scheme.

3. 一维高对比度谐振器链的长程相互作用与安德森局域化

Long-range interactions and Anderson localisation for one-dimensional high-contrast resonator chain

作者: Habib Ammari, Silvo Barandun, Jiayu Qiu, Alexander Uhlmann

类别: math-ph, math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04971v1

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

摘要: Spectral and transport properties of high-contrast resonator systems can be described in the subwavelength regime in terms of the so-called capacitance operator. In this paper, we consider an infinitely periodic chain of high-contrast resonators in three dimensions. The first result is a precise estimate of the off-diagonal decay rate of the capacitance operator $C$. Importantly, we demonstrate that the decay rate is long-range and critical: as $|n-m|\to\infty$, \begin{equation*} C(n,m)\sim \frac{1}{|n-m|\log^2|n-m|}, \end{equation*} which is $\ell^1$ summable but slower than quadratic. This borderline decay of the off-diagonal entries makes the present proof of Anderson localisation with arbitrary disorder, which is observed numerically in this paper, out of reach; we hope that this physical example of classical wave systems with critical long-range interactions provides new insight in the field of Anderson localisation. As the second main result, based on the off-diagonal decay estimate, we prove a strong convergence of the finite capacitance operator, which corresponds to a truncated chain, to the capacitance operator as the size of the truncated chain grows to infinity. Using this strong convergence, we improve the results of [Ammari et al., SIAM J. Math. Anal., 2023 and Bull. London Math Soc., 2025] by presenting a rigorous estimate of the convergence rate of the spectrum.

4. 抛物型方程中的阈值、碎片化与对称化

Thresholds, fragmentation and symmetrization in parabolic equations

作者: Matthieu Alfaro, Yihong Du, François Hamel, Lionel Roques

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04807v1

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

摘要: This paper is mainly concerned with the large-time dynamics of bounded nonnegative solutions of reaction-diffusion equations on the real line with nonlinearities mainly of the bistable-type. We first consider initial data of the type $α\mathbf{1}_I$, namely scalar multiples of indicator functions of intervals $I$. For each amplitude $α$, the existence of a threshold length $L^*(α)$ separating the extinction and the persistence of the solutions at large time is known. We here address the question as to whether the limit of the threshold sizes $L^*(α)$ as $α\to+\infty$ is positive or zero. We provide sufficient conditions under which this limit is positive, and others under which it is zero. Secondly, when the limit is positive, we show that some fragmented initial data, which are equally distributed as $α\mathbf{1}_I$, give rise to solutions persisting at large time, whereas the solutions emanating from $α\mathbf{1}_I$ go to extinction at large time. This result shows an unexpected favourable effect of the fragmentation of the initial datum on the large-time dynamics. Last, we consider the mass concentration principle for parabolic equations, which states that nonnegative solutions can be controlled from above in an integral sense by the solutions emanating from the Schwarz symmetrically decreasing rearrangements of the initial data. Since the pioneering results of [Bandle, 1976] and [Alvino, Trombetti, Lions, 1990], this principle is known to hold in bounded domains with Dirichlet boundary conditions under different types of assumptions of the coefficients of the equation. We show that this mass concentration principle is not valid in general, even for simple equations of the type $\partial_t u = \partial_{xx} u + f(u)$.

5. 具有通常增长性的正交各向异性泛函的Lipschitz正则性

Lipschitz regularity for orthotropic functionals with general growth

作者: Mikyoung Lee, Jihoon Ok, Bianca Stroffolini

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04773v1

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

摘要: We study the local Lipschitz regularity of local minimizers for a class of degenerate orthotropic functionals with $\varphi$-growth, where $\varphi$ is a general N-function. Unlike standard isotropic functionals, the ellipticity of the associated Euler-Lagrange equation degenerates separately in each coordinate direction, presenting significant anisotropic difficulties. Furthermore, the general N-function setting lacks the algebraic scale invariance available in the classical orthotropic $p$-Laplacian case. Despite these structural difficulties, we prove that local minimizers are locally Lipschitz continuous. Our approach relies on a regularized approximation scheme, mixed-direction Caccioppoli inequalities, and a carefully designed Moser-type iteration that incorporates an interpolation argument to bridge the gaps between consecutive integrability exponents.

6. 具有 $L^1$ 数据的复 Hessian 方程的稳定性与强收敛性

Stability and strong convergence for complex Hessian equations with $L^1$ data

作者: Truong Dinh Dat

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

发布日期: 2026-07-06

arXiv ID: 2607.04704v1

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

摘要: We study complex $m$-Hessian equations on bounded hyperconvex domains with right-hand side in $L^1(Ω)$. The main contribution of this paper is a strong stability result for weak solutions in the Hessian energy sense. More precisely, if $f_j \to f$ in $L^1(Ω)$ and $u_j, u$ are the corresponding solutions, then

This provides convergence in the natural energy topology associated to the Hessian operator, which is significantly stronger than convergence in capacity. For completeness, we also recall the existence of solutions and stability in capacity, which follow from known results in the literature.

7. 仿射截面层析成像在具有受限大边界数据的 $k$-Hessian 方程反源问题中的应用

Affine section tomography for inverse source problems in $k$-Hessian equations with restricted large boundary data

作者: Yi-Hsuan Lin

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04662v1

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

摘要: We consider an inverse source problem for the $k$-Hessian equation \begin{equation*} σ_k(D^2u)=f(x) \end{equation*} on the $k$-admissible branch in a smooth uniformly convex domain in $\mathbb R^n$, where $2\le k\le n$. We prove that the nonlinear Dirichlet-to-Neumann map determines the positive smooth source from its values on the restricted large-data rays $tφ_E|_{\partialΩ}$, where $E\in Gr(k-1,n)$, $φ_E(x)=|P_Ex|^2/2$, and $t$ is sufficiently large. The Hessian of each boundary profile has exactly $k-1$ large directions and is flat on $V=E^\perp$. Thus, the leading profile lies on a rank $k-1$ face of the $k$-Hessian structure, and the first source-dependent correction is governed by the missing directions. More precisely, this correction solves fiberwise Poisson equations on the affine sections $Ω\cap(y+V)$ of dimension $q=n-k+1$. We prove that these sectionwise solutions patch smoothly through glancing points where the sections collapse, and we obtain boundary normal derivative asymptotics by local barriers. The boundary flux of the correction gives the section integrals $\int_{Ω\cap(y+V)}f\,d\mathcal H^q$. Varying $E$ yields the affine $q$-plane Radon transform of the zero extension of $f$. We give an explicit reconstruction formula through the Fourier slice identity, and the injectivity of the affine Radon transform gives uniqueness. The endpoint $k=n$ recovers the Monge–Ampère chord/X-ray geometry, while the range $n\ge3$ and $2\le k<n$ gives inverse source results for genuinely non-determinant Hessian equations.

8. 在预设高斯曲率方程的反问题中的规范刚性

Gauge rigidity in an inverse problem for the prescribed Gaussian curvature equation

作者: Yi-Hsuan Lin

类别: math.AP

发布日期: 2026-07-06

arXiv ID: 2607.04657v1

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

摘要: We prove a uniqueness result for an inverse boundary problem associated with the prescribed Gaussian curvature equation

for graphs over a planar domain. The comparison is made on a common open class of smooth boundary values for which both admissible Dirichlet problems are well posed. We show that, if the corresponding nonlinear Dirichlet-to-Neumann maps agree on this class and the two prescribed curvatures have the same first boundary jet, then the curvatures agree in the whole domain. The main difficulty comes from a gauge obstruction already present at the first linearization. In logarithmic form, the linearized equation is a two-dimensional non-divergence form elliptic equation with drift. Its boundary data determine the coefficients only up to a boundary-fixing change of variables and a scalar gauge factor. For the prescribed Gaussian curvature equation this gauge is not an artifact of the method: the gradient dependence leaves a residual gauge which cannot be removed by the first variation. The proof uses the nonlinear structure to remove this remaining gauge. We derive a second-linearized interaction identity in the first-linearized gauge. In this identity, the covariant Hessian of the gauge map appears in the leading part. Testing the identity with two complementary CGO families gives a pair of residual equations for the gauge variables. Combined with the drift and conductivity identities from the first linearization, these equations form a closed system for the residual gauge. A boundary unique continuation argument for this system gives the trivial gauge and hence determines the prescribed curvature.

9. 希尔伯特空间中的非交换各向异性扩散。II. 对数梯度的全局闭包、下界及纳米系统应用

Noncommutative Anisotropic Diffusion in Hilbert Space. II. Global Closure of the Logarithmic Gradient, Lower Bounds, and Nanosystem Applications

作者: E. Yu. Shchetinin, A. A. Shevchuk, S. I. Salpagarov

类别: math.AP, math.PR

发布日期: 2026-07-05

arXiv ID: 2607.04131v1

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

摘要: This second part of the series develops the statistical and applied layer of the theory built in Part I [1]. Unlike current Hilbert diffusion models [2-5], we focus not only on well-posedness of the infinite-dimensional generative dynamics, but on the noncommutative $A$-geometry, explicit entropy constants, and checkable lower bounds. The analytic estimate of Part I reduces stability of the backward evolution to control of a validation error $\mathcal{E}_{val}$. We prove three results. First, we construct a cylindrically weak denoising label for the infinite-dimensional score field, consistent with the exact logarithmic gradient. Second, the local parametric closure is replaced by a global nonparametric score closure, its complexity controlled via the Dudley entropy integral and uniform empirical bounds in $L^2(ν_*;A)$. Third, for the trace-smoothed nonlinear class we construct minimax bounds by the Le Cam-Assouad method, showing the root statistical rate cannot be improved without additional quadratic structure. The final section applies the theory to anisotropic diffusion in a nanosystem model and checks the constants $c_A,C_A,C_{LSI}^A$ independently of any smallness condition, establishing explicit accuracy orders: the parametric trace-smoothed minimax score risk is $p/M$, while the uniform error of the risk functional is $\sqrt{p/M}$, showing the statistical plateau of order $p/M$ is not a proof artifact. The applied layer closes with an independent analytic benchmark for the isotropic case, a comparison with classical cell homogenization, and an approximation theorem for smooth logarithmic gradients by $A$-adapted spectral networks.

10. $L^p$-极值Teichmüller映射在黎曼曲面之间是微分同胚

$L^p$-Extremal Teichmüller mappings between Riemann surfaces are diffeomorphisms

作者: Gaven Martin, Cong Yao

类别: math.CV, math.AP

发布日期: 2026-07-04

arXiv ID: 2607.04051v1

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

摘要: We consider minimisers in the homotopy class of a homeomorphism $f_0:R\to S$ between analytically finite Riemann surfaces with minimal $L^p$- conformal energy

The problem was first raised by Ahlfors in his celebrated proof of Teichmüller’s theorem-the case $p=\infty$, but the existence, topological regularity and analytic regularity of these $L^p$ minimisers remained unknown for all $1<p<\infty$. Ahlfors established weak existence for $p\geq 2$. Here we prove that for all p, $1\leq p<\infty$, such minimisers exist, are unique and are diffeomorphisms. They are quasiconformal but not diffeomorphic at $p=\infty$.

11. De Giorgi 测度解与平均曲率流：一种极小化运动方法

De Giorgi varifold solutions to Mean Curvature Flow: a minimizing movements approach

作者: Tim Laux, Andrea Poiatti

类别: math.AP, math.DG

发布日期: 2026-07-04

arXiv ID: 2607.03930v1

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

摘要: We propose an alternative existence proof of global weak solutions to mean curvature flow and volume preserving mean curvature flow. We prove for the first time for a minimizing movements scheme the unconditional convergence towards a varifold solution, here a De Giorgi solution. The argument is purely variational and does not rely on comparison principles. The key novelty is an alternative proxy for the completely degenerate $L^2$ distance that is more robust than the one of Almgren-Taylor-Wang and Luckhaus-Sturzenhecker.

12. 趋化驱动的肿瘤–免疫空间模式与稳定性的计算肿瘤学

Computational Oncology of Chemotaxis-Driven Tumour–Immune Spatial Patterning and Stability

作者: Zonghao Liu, Jiguang Yu, Lei Su, Louis Shuo Wang, Yang Du, Jingfeng Liu

类别: math.AP, math.NA

发布日期: 2026-07-04

arXiv ID: 2607.03813v1

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

摘要: Spatial tumour–immune heterogeneity is a key feature of solid-tumour progression, immune infiltration, and immune exclusion. We develop a computational oncology model in which tumour cells, immune effector cells, and a chemokine signal interact through a reaction–diffusion–chemotaxis system on a bounded tissue domain with no-flux boundaries. Chemokine is produced by tumour cells and tumour–immune contact, recruits immune cells, and guides chemotactic migration. After nondimensionalization, we establish positivity, a tumour-density bound, and immune/chemokine mass estimates. We identify the tumour-free equilibrium, derive the immune-control threshold $σ_0>δ$, and reduce coexistence to a scalar equation. Linear stability analysis about coexistence yields a mode-wise dispersion relation in which chemotaxis appears as a wavenumber amplified coupling, producing finite-wavelength instability above a critical sensitivity. A conservative finite-volume scheme with upwind chemotactic flux verifies the thresholds, dominant unstable modes, sensitivity maps, positivity, convergence, and residual consistency.

13. 具有环境反馈与固定招募通量的尺度结构种群的最优收获

Optimal Harvesting of Size-Structured Populations with Environmental Feedback and Fixed Recruitment Flux

作者: Louis Shuo Wang, Jiguang Yu

类别: math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03790v1

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

摘要: We study a nonlinear size-structured transport model with distributed harvesting and prescribed recruitment flux, where environmental feedback is determined by a scalar population functional. After establishing global well-posedness on $L^1$, we reduce stationary equilibria to a scalar closure equation. This reduction reveals that loss of equilibrium uniqueness occurs through a generic fold, mathematically characterizing critical population transitions. On uniformly nonresonant equilibrium branches, we prove the existence of optimal stationary harvesting policies via the direct method of the calculus of variations. We then derive a boundary-corrected adjoint equation and establish an identity equating equilibrium sensitivity with the adjoint loop gain. This relation yields explicit criteria for the persistence and creation of optimal harvesting thresholds. Collectively, these results provide a unified analytical framework connecting environmental feedback, equilibrium structure, and optimal harvesting.

14. Riesz空间中强单调算子方程的解

Solutions for Strongly Monotone Operator Equations in Riesz Spaces

作者: Xian Xu, Baoxia Qin

类别: math.FA, math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03743v1

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

摘要: This paper is devoted to the study of solutions for a class of operator equations governed by strongly monotone operators on real Riesz spaces continuously embedded into Banach lattices. By exploiting the intrinsic lattice structure of Banach lattices, we establish refined growth assumptions on nonlinear terms to guarantee that suitable neighbourhoods of positive and negative cones are invariant under the descending flow. Combining the descending flow invariant set technique with the theory of strongly monotone operators, we derive abstract existence theorems: the operator equation possesses at least one positive solution, one negative solution and one sign-changing solution. These abstract results are further applied to $(p,q)$-Laplacian boundary value problems, yielding corresponding multiplicity conclusions on positive, negative and sign-changing solutions.

15. Riesz 空间中的抽象不定问题及其应用

Abstract Indefinite Problems in Riesz Spaces with Its Applications

作者: Xian Xu, Baoxia Qin

类别: math.FA, math.AP

发布日期: 2026-07-04

arXiv ID: 2607.03742v1

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

摘要: This paper investigates the existence of critical points for functionals defined on a Hilbert space $X$ which is continuously embedded into a Banach lattice $E$. A lattice decomposition of $E$ is constructed, which possesses both order disjointness and inner-product orthogonality. Accordingly, a corresponding decomposition of the Hilbert space $X$ is obtained. Under this decomposition, the associated functional satisfies the energy collapse condition and order-preserving property on certain subspaces, while exhibiting coerciveness on others. By combining the descending flow invariant set method with Morse theory, we establish the existence of multiple critical points for the abstract indefinite problems. Finally, applications to elliptic boundary value problems are provided.

16. LRX-PINN：一种具有集成柯西激活函数的层解析XNet物理信息神经网络，用于对流主导问题

LRX-PINN: A Layer-Resolving XNet Physics-Informed Neural Network with Integrated Cauchy Activations for Convection-Dominated Problems

作者: Zihao Guo, Xin Li, Zhihong Xia

类别: math.AP, cs.LG

发布日期: 2026-07-04

arXiv ID: 2607.03682v1

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

摘要: Convection-dominated convection-diffusion problems often develop thin layers, where the solution has sharp transition profiles and its derivatives are highly localized. This creates a structural mismatch for standard physics-informed neural networks (PINNs), whose trial spaces are not designed to match the value–derivative structure of such layers. We propose a Layer-Resolving XNet Physics-Informed Neural Network (LRX-PINN) based on integrated Cauchy activations. The proposed basis is transition-type at the solution level, while its derivative recovers a localized Cauchy kernel. We show that this structure matches the scaling of convection-dominated layers, inherits the Cauchy approximation mechanism at the derivative-profile level, and identifies $d/\|w\|$ as the effective physical width of a ridge neuron. For analytic layer profiles, this yields derivative-stable exponential approximation in the stretched coordinate and a layer-scaled estimate for the strong residual of the singularly perturbed operator. Numerical experiments on several convection-dominated benchmarks show that LRX-PINN achieves higher accuracy than PIKAN and Fourier-feature PINNs while using less than $30\%$ of their trainable parameters. On more challenging benchmarks, embedding the proposed representation into hp-VPINN-based frameworks further improves the best results obtained by existing hp-VPINN-based baselines without changing their original loss functionals or stabilization strategies. These results show that neural representations aligned with layer structure provide a compact and effective approach for convection-dominated problems.

17. Hölder映射在Pfaff约束下

Hölder maps under Pfaffian constraints

作者: Armin Schikorra

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

发布日期: 2026-07-04

arXiv ID: 2607.03667v1

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

摘要: Given a one form $λ$ in $\mathbb{R}^N$ and $f: \mathbb{S}^{n} \to \mathbb{R}^N$ with $f^\ast λ= 0$ we discuss the maximal Hölder regularity of extensions $F: \mathbb{B}^{n+1} \to \mathbb{R}^N$ such that $F^\astλ= 0$ in distributional sense. Our analysis applies to the Heisenberg groups $\mathbb{H}_n$. It implies in particular that for all $n \geq 1$ any smooth horizontal map $f: \mathbb{S}^{n} \to \mathbb{H}_n$ can be extended to a $C^α$-map $F: \mathbb{B}^{n+1} \to \mathbb{H}_n$ for some $α> 1/2$. Moreover, if $n \geq 3$ we find $C^α$-embeddings from $\mathbb{B}^{n+1}$ into $\mathbb{H}_n$ for some $α> \frac{1}{2}$.