Home » 论文总结

2026年07月14日论文分析

共有 10 篇重点关注论文, 22 篇论文可以了解, 22 篇不相关论文

July 14, 2026 · 56 min · 11762 words · P0lyno3ial

Table of Contents

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

生成时间: 2026年07月14日 05:21:45

运行状态: 完成

计划处理论文数量: 54

已完成论文数量: 54

未完成论文数量: 0

重点关注论文数量: 10

了解领域论文数量: 22

不相关论文数量: 22

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

1. 关于三维Peskin问题在临界正则性下的全局渐近稳定性

On the global asymptotic stability for the 3D Peskin Problem at critical regularity

作者: Eduardo García-Juárez, Susanna V. Haziot, Po-Chun Kuo, Yoichiro Mori, Han Zhou

类别: math.AP, math-ph

发布日期: 2026-07-13

arXiv ID: 2607.11731v1

摘要翻译: 我们证明了三维Peskin问题的全局适定性与渐近稳定性，该问题建模了一个浸没于不可压Stokes流体中的闭弹性膜。我们处理初始数据属于最优正则性空间 $W^{1,\infty}(\mathbb{S}^2)$ ，该空间允许包含无穷多个角点。这些初始构型通过流的抛物正则化效应被瞬时去奇异性，在所有 $$t > 0$$ 时刻变为光滑。随后我们建立了解在 $$C^1$$ 拓扑下以指数速率收敛至一个平移并缩放后的共形球面。

Comment: 2 figures

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

1. 研究对象和背景

本文研究三维Peskin问题，即一个二维闭合弹性膜 $X = (X^1, X^2, X^3) : [0,T] \times \mathbb{S}^2 \to \mathbb{R}^3$ 浸没于不可压Stokes流体中的自由边界流固耦合模型。该模型源于Peskin浸没边界法，其数学核心是边界积分形式的非局部非线性演化方程：

(\partial_t X^i)(t, x) = \int_{\mathbb{S}^2} G_{ij}\big(X(t,x) - X(t,p)\big) \, \Delta_{\mathbb{S}^2} X^j(t,p) \, d\mu(p), \quad i=1,2,3,\; x \in \mathbb{S}^2,

其中 $\Delta_{\mathbb{S}^2}$ 是球面上的Laplace-Beltrami算子， $d\mu$ 是单位球面诱导测度， $G_{ij}(y) = \frac{1}{8\pi} \left( \frac{\delta_{ij}}{|y|} + \frac{y_i y_j}{|y|^3} \right)$ 是三维Stokeslet核。物理上，该方程描述Hookean型弹性膜在Stokes流中受自身曲率驱动的运动：膜张力由Dirichlet能量 $\int_{\mathbb{S}^2} |\nabla_{\mathbb{S}^2} X|^2 d\mu$ 生成，而流体响应通过Stokeslet以非局部方式反作用于界面。

关键未知量是映射 $$X$$ 本身——它同时编码界面几何（像集）与参数化（作为 $\mathbb{S}^2$ 到 $\mathbb{R}^3$ 的嵌入）。这导致一个根本性区别：与仅依赖形状的流体界面问题（如Muskat或Stokes液滴）不同，Peskin问题的稳态不仅包含欧氏对称性（平移、旋转、缩放），更因二维Dirichlet能量在共形变换下的不变性，产生额外的有限维稳态流形。具体而言， $\mathbb{S}^2$ 上保向共形微分同胚群 $\mathrm{SO}^+(3,1)$ 作用于恒等映射 $$X_0(p)=p$$ ，生成一族稳态 $S_\Lambda(p)$ ；再叠加物理平移 $b \in \mathbb{R}^3$ 与缩放 $$r > 0$$ ，得到完整的10维稳态流形：

\mathcal{M}_\mathrm{steady} = \big\{ S_{r,b,\Lambda}(p) := b + r S_\Lambda(p) \,\big|\, r > 0,\, b \in \mathbb{R}^3,\, \Lambda \in \mathrm{SO}^+(3,1) \big\}.

该流形的切空间由10个线性无关的零模构成：3个平移生成元 $\{\partial_{b_k}\}$ 、3个旋转生成元 $\{V_{ij}\}$ （见(2.7)）、1个缩放生成元 $\partial_r$ ，以及3个特殊共形生成元（对应 $\mathrm{SO}^+(3,1)$ 的非紧子群）。这些零模正是线性化算子 $\mathcal{N}_1$ 的核（Lemma 3.5），其存在使得标准能量方法失效——扰动沿 $\mathcal{M}_\mathrm{steady}$ 方向不衰减，必须发展新的正交分解框架。

已有理论主要集中于二维情形：Lin–Tong [21] 与 Mori–Rodenberg–Spirn [25] 建立了次临界正则性（如 $C^{1+\alpha}$ ）下的局部适定性与小初值收敛；Chen–Nguyen [7] 将初值空间推进至临界Besov空间 $\dot{B}^1_{\infty,\infty}$ ；García-Juárez–Haziot [13] 在允许角点的临界空间中实现全局适定性。但所有这些工作均未解决三维情形的核心障碍：(i) 几何复杂性导致共形稳态流形维度跃升至10维（二维仅为4维），(ii) 非线性项 $G_{ij}(X(x)-X(p)) \Delta_{\mathbb{S}^2} X^j(p)$ 在 $W^{1,\infty}$ 初值下高度奇异——当 $$X$$ 仅 Lipschitz 连续时， $$X(x)-X(p)$$ 可退化为零，使Stokeslet核 $|X(x)-X(p)|^{-3}$ 失去可积性；(iii) 标准Sobolev或Besov空间无法控制该非线性在球面上的谱行为，尤其当高频分量 $$P'_j Y$$ 与低频 $$P'_k Y$$ 相互作用时，经典乘积估计 $\|fg\|_{L^\infty} \le \|f\|_{L^\infty}\|g\|_{L^\infty}$ 完全失效。因此，本文试图填补的关键空白是：在最优初值正则性 $W^{1,\infty}(\mathbb{S}^2)$ （允许无限多角点）下，建立三维Peskin问题的全局适定性，并证明解指数收敛至10维共形稳态流形——这要求同时克服几何、分析与代数三重障碍：构造适配球面共形结构的函数空间、控制Stokeslet非线性的临界奇异性、并实现稳态流形与其正交补空间的精确解耦。

2. 主要定理或主要结果

本文在初值满足 $\|Y_0\|_{L^\infty} + \|\nabla Y_0\|_{L^\infty} \le \varepsilon \le \varepsilon_0$ 的条件下，证明以下核心结论：

(i) 【局部适定性】存在唯一解 $$X = X_0 + Y$$ 定义在 $$[0,1]$$ 上，且满足 $\|Y\|_{Z[0,1]} \lesssim \varepsilon$ ，其中函数空间 $$Z[0,1]$$ 由范数

\|f\|_{Z[0,1]} := \sup_{t \in [0,1]} \Big\{ \|f(t)\|_{L^\infty} + \|\nabla f(t)\|_{L^\infty} + \sup_{j \in \mathbb{Z}_+} 2^j (1 + 2^j t)^{\delta_0} \|P'_j f(t)\|_{L^\infty} \Big\}, \quad \delta_0 = \frac{1}{8},

定义， $$P'_j$$ 为基于球谐函数的谱Littlewood-Paley投影（见(2.35)）。该范数的关键创新在于时间权重 $(1 + 2^j t)^{\delta_0}$ ：它精确捕获了高频分量 $$P'_j Y$$ 的抛物正则化速率 $\sim 2^{-j}$ ，从而在 $W^{1,\infty}$ 初值下实现即时光滑化（ $$t>0$$ 时 $Y(t) \in C^\infty$ ）。

(ii) 【全局适定性与渐近稳定性】解可延拓为全局解 $X \in C((0,\infty); C^1(\mathbb{S}^2))$ ，且存在最终状态参数 $(r_\infty, b_\infty, \Lambda_\infty) \in (0,\infty) \times \mathbb{R}^3 \times \mathrm{SO}^+(3,1)$ 满足 $|(r_\infty,b_\infty,\Lambda_\infty) - (1,0,I_4)| \lesssim \varepsilon$ ，使得

\|X(t) - S_{r_\infty,b_\infty,\Lambda_\infty}\|_{C^1(\mathbb{S}^2)} \le C_\lambda \varepsilon e^{-a_2 / r_\infty \, t}, \quad \forall t \ge 1,

其中 $$a_2 = 8/35$$ ， $S_{r,b,\Lambda}(\alpha) = b + r S_\Lambda(\alpha)$ 。该指数衰减率由线性化算子在稳态流形正交补上的最优谱间隙决定（Proposition 5.11），且最终半径 $r_\infty$ 由体积守恒律 $\int_{\mathbb{S}^2} \det(\partial_\theta X \wedge \partial_\phi X) \, d\theta d\phi = \frac{4}{3}\pi r_\infty^3$ 唯一确定。

(iii) 【稳态分类】在单位球邻域内，共形球是唯一的稳态；结合Changyou Wang的通信结果，三维Peskin问题的所有稳态恰为 $\mathcal{M}_\mathrm{steady}$ 中的元素，此结论源于二维Dirichlet能量的共形不变性与Hopf定理。

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

证明的核心障碍有三：(1) 非线性项 $N^i(X) = \int_{\mathbb{S}^2} G_{ij}(X(x)-X(p)) \Delta_{\mathbb{S}^2} X^j(p) \, d\mu(p)$ 在 $W^{1,\infty}$ 初值下具有临界奇异性——当 $$X$$ 仅为Lipschitz时， $|X(x)-X(p)| \sim |x-p|$ ，导致被积函数含 $|x-p|^{-3}$ 因子，在 $\mathbb{S}^2$ 上不可积；(2) 线性化算子 $\mathcal{N}_1$ 具有10维零空间，标准能量估计无法控制沿稳态流形的漂移；(3) 球面非交换几何结构使经典欧氏Littlewood-Paley理论失效，需重建适配共形对称性的谱分析框架。

作者的真正创新点在于构建了一个三位一体的分析体系：

第一，提出谱Littlewood-Paley框架 $$Z[0,T]$$ 。该空间范数中的时间权重 $(1 + 2^j t)^{\delta_0}$ 并非随意选取，而是精确匹配Stokes流的抛物正则化效应：对频率 $\sim 2^j$ 的分量，热核衰减尺度为 $t \sim 2^{-2j}$ ，故权重 $2^j (1 + 2^j t)^{1/8}$ 在 $t \ll 2^{-j}$ 时主导 $L^\infty$ 范数，在 $t \gg 2^{-j}$ 时主导导数范数，从而统一控制瞬态奇点与长期正则性。关键引理4.1利用球面矢量场 $V_{12}, V_{23}, V_{31}$ （见(2.7)）的散度为零性及Laplace-Beltrami算子的二次表示 $\Delta_{\mathbb{S}^2} = V_{12}^2 + V_{23}^2 + V_{31}^2$ （见(2.10)），通过反复分部积分将非线性估计转化为对 $$P'_j Y$$ 的 $L^\infty$ 控制，避免了对 $\nabla Y$ 的直接估计，突破了 $W^{1,\infty}$ 的临界壁垒。

第二，实现稳态流形的精确结构解耦。作者将解 $$X(t)$$ 投影到10维稳态流形 $\mathcal{M}_\mathrm{steady}$ 及其正交补 $\mathcal{M}_\perp$ ：设 $p(t) = (r(t), b(t), \Lambda(t))$ 为流形上最近点参数，则扰动 $Y_\perp(t) = X(t) - S_{p(t)}$ 满足正交约束 $\langle Y_\perp(t), \psi_k(S_{p(t)}) \rangle_{L^2} = 0$ （ $\psi_k$ 为 $\mathcal{M}_\mathrm{steady}$ 的切向基）。通过计算 $\frac{d}{dt} p(t)$ 的显式公式（利用 $\mathcal{N}_1$ 的零空间结构），证明 $$p(t)$$ 收敛至极限 $p_\infty$ ，且 $Y_\perp(t)$ 满足线性化方程 $\partial_t Y_\perp = \mathcal{N}_1(Y_\perp) + \mathcal{R}(Y_\perp)$ ，其中余项 $\mathcal{R}$ 被证明是高阶小量。这一解耦严格分离了中性方向（由ODE控制）与耗散方向（由谱间隙控制）。

第三，建立球面谱 multiplier 理论。在Section 6中，作者证明：对任意满足 $|m(k)| \lesssim (1+k)^{-\sigma}$ 的谱乘子 $$m$$ ，算子 $T_m f = \sum_{k\ge 0} m(k) f_k$ （ $f_k \in H_k$ ）在 $L^p(\mathbb{S}^2)$ 上有界，且其核 $K_m(x,y) = \sum_{k\ge 0} m(k) Z_k(x,y)$ 满足 $|K_m(x,y)| \lesssim d(x,y)^{-2+\sigma}$ （ $$d$$ 为球面距离）。该理论为线性化算子 $\mathcal{N}_1$ 的谱分解（Proposition 3.4）及非线性估计提供了严格基础，其证明依赖于zonal spherical harmonics $$Z_k$$ 的精确渐近（(2.29)-(2.31)）与Legendre多项式 $$P_k$$ 的振荡性质。

4. 与之前工作的比较

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

【正则性维度】此前所有三维结果均限于次临界空间：García-Juárez–Po-Chun–Mori–Strain [15] 仅在 Hölder 空间 $C^{1+\alpha}$ 中建立局部适定性，其证明依赖于 $$X$$ 的 $C^{1+\alpha}$ 一致性控制以保证 $|X(x)-X(p)| \gtrsim |x-p|^{1+\alpha}$ ，从而避免Stokeslet奇点。本文将初值空间降至 $W^{1,\infty}$ ——这是最优临界空间，因 $W^{1,\infty}$ 恰好是 Lipschitz 映射空间，而角点（即不可微点）是弹性膜物理建模的自然特征。代价是必须放弃经典乘积估计，转而依赖谱Littlewood-Paley框架与球面矢量场分部积分技术。

【稳态结构处理】二维工作 [21,25] 仅需处理4维稳态流形（3维欧氏对称+1维缩放），其正交补空间的谱间隙可直接由Fourier级数计算。本文面对10维共形流形，首次实现其精确解耦：以往方法（如模态稳定性分析）仅能证明收敛至某稳态，但无法确定最终参数 $(r_\infty,b_\infty,\Lambda_\infty)$ ；本文通过构造流形上测地投影并求解参数ODE，不仅证明收敛，且给出 $r_\infty$ 由体积守恒决定、 $b_\infty$ 由质心守恒决定、 $\Lambda_\infty$ 由共形不变量决定的完整图景。特别地，指数率 $e^{-a_2 / r_\infty t}$ 中的 $$a_2 = 8/35$$ 是线性化算子在 $\mathcal{M}_\perp$ 上的最优谱半界，其计算依赖于球谐函数 $$H_k$$ 上 $\mathcal{N}_1$ 的精确本征值（Proposition 3.4指出 $\mathcal{N}_1$ 在 $$H_k$$ 上作用为乘子 $$-k(k+1)/(2k+1)$$ ，故最小非零本征值在 $$k=2$$ 时取得 $$-6/5$$ ，经缩放得 $$a_2$$ ）。

【适用范围与局限】本文结果要求初值 $$C^1$$ -接近单位球，而 [13] 在二维中允许初值远离圆但保持弦弧条件。三维中弦弧条件难以定义，故本文采用更严格的邻域假设。此外，指数衰减仅对 $t \ge 1$ 成立，初始阶段 $t \in (0,1)$ 的行为由抛物正则化主导，此时解从角点初值瞬时光滑化，但衰减速率尚未显现。数值验证（Section 7）确认了理论预测的 $e^{-ct}$ 衰减，但未检验大初值情形——作者明确指出“不排除远离单位球存在其他稳态”，此为当前理论的固有局限。相较二维Navier-Stokes情形 [36]，本文完全忽略流体惯性（即使用Stokes而非Navier-Stokes方程），故不适用于高雷诺数流动；但这也使得线性化算子成为纯耗散算子，为指数稳定性提供了必要条件。

2. 空间异质环境下混合反应-扩散传染病模型中的波前传播. 第 II 部分：脉动行波

Front propagation in hybrid reaction-diffusion epidemic models with spatial heterogeneity. Part II: Pulsating traveling waves

作者: Quentin Griette, Hiroshi Matano

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11869v1

摘要翻译: 我们考虑一个一维空间中的双物种反应-扩散系统，该系统源自空间周期环境中具有两种病原体（野生型与突变型）的流行病学模型。该系统具有混合性质：部分合作、部分竞争，但又不完全属于其中任一类型；因此比较原理不成立。在前期工作中，我们研究了该系统初值问题解的传播性质，并证明了左右两个方向的波前传播速度 $$ c^*_R$$ 与 $$ c^*_L$$ 可通过某些主特征值表征；同时研究了空间周期 $L \to 0$ 时的均匀化极限，并讨论了波前后方解的长时间行为。本文中，我们证明：对任意 $c\geq c^*_R$ （相应地， $c\geq c^*_L$ ），存在向右方向（相应地，向左方向）传播、速度为 $$c$$ 的脉动行波，其中 $$c^*_R$$ 与 $$c^*_L$$ 即前述左右方向的传播速度。我们还证明：任意行波的前缘具有形式线性分析所预期的指数衰减率，从而将 Hamel 2008 的部分结果推广至方程组情形。最后，我们给出一个例子，其中 $c^*_R \neq c^*_L$ ；该例子通过考虑一个多尺度奇异极限问题构造而成。这一结果凸显了我们的系统与标量 KPP 型方程之间的显著差异。

Comment: arXiv admin note: substantial text overlap with arXiv:2408.07501

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

1. 研究对象和背景

本文研究如下一维空间周期异质环境中的两物种反应-扩散系统：

\begin{cases} \partial_t u = \big(\sigma(x)\, \partial_x u\big)_x + \big(r_u(x) - \kappa_u(x)(u+v)\big)u + \mu_v(x)v - \mu_u(x)u, \\ \partial_t v = \big(\sigma(x)\, \partial_x v\big)_x + \big(r_v(x) - \kappa_v(x)(u+v)\big)v + \mu_u(x)u - \mu_v(x)v, \end{cases} \quad t>0,\ x\in\mathbb{R},

其中 $u(t,x),v(t,x)\ge0$ 分别表示野生型与突变型病原体的空间密度；所有系数 $\sigma(x)>0$ , $r_u(x),r_v(x)\in\mathbb{R}$ , $\kappa_u(x),\kappa_v(x)>0$ , $\mu_u(x),\mu_v(x)>0$ 均为 $$L$$ -周期连续函数，且 $\sigma\in C^1(\mathbb{R})$ 。该系统源于SIS型传染病模型（1.2），经总宿主数 $$N=S+I_1+I_2$$ 恒定假设约化而来，故其动力学完全由（1.1）刻画。

该系统具有本质的“混合性”（hybrid nature）：当 $$(u,v)$$ 较小时，交叉项 $\mu_v v - \mu_u u$ 与 $\mu_u u - \mu_v v$ 主导，系统呈合作性（cooperative），即 $\partial_u f_v \ge 0$ , $\partial_v f_u \ge 0$ ；而当 $$(u,v)$$ 较大时，竞争项 $-\kappa_u(u+v)u$ 与 $-\kappa_v(u+v)v$ 主导，系统转为竞争性（competitive）。因此，整个系统既非整体合作亦非整体竞争，标准比较原理（comparison principle）全局失效——这是分析的核心障碍。相比之下，经典KPP型标量方程或纯合作/纯竞争系统均满足比较原理，其行波存在性、单调性、唯一性等可借助上下解法、单调迭代或单调动力学方法处理；而本系统的非单调耦合结构导致传统工具失灵。

研究动机来自抗药性病原体在空间异质环境（如抗生素梯度、禽舍笼格排布、农田农药斑块）中的传播建模。空间异质性不仅改变局部增殖率 $$r_*(x)$$ 和竞争强度 $\kappa_*(x)$ ，更通过影响突变速率 $\mu_*(x)$ 调控基因漂变方向，从而可能诱发不对称传播现象。已有工作[15]已建立该系统的传播速度理论：定义线性化系统

\begin{cases} \partial_t u = \big(\sigma(x)\, \partial_x u\big)_x + r_u(x)u + \mu_v(x)v - \mu_u(x)u, \\ \partial_t v = \big(\sigma(x)\, \partial_x v\big)_x + r_v(x)v + \mu_u(x)u - \mu_v(x)v, \end{cases}

其 $\lambda$ -周期主特征值 $k(\lambda)$ （见定义2）为严格凸解析函数，右、左传播速度由Gärtner-Freidlin型公式给出：

c_R^* = \min_{\lambda>0} \frac{k(\lambda)}{\lambda}, \quad c_L^* = \min_{\lambda<0} \frac{k(\lambda)}{-\lambda}.

当系数满足对称条件（如 $\mu_u=\mu_v$ 或所有系数偶函数）时， $k(\lambda)=k(-\lambda)$ ，故 $$c_R^*=c_L^*$$ ；但一般情形下 $k(\lambda)$ 非偶， $$c_R^*$$ 与 $$c_L^*$$ 可能不等——这与标量KPP方程根本不同（后者因Fredholm替代性必有 $$c_R^*=c_L^*$$ ）。然而，[15]仅处理初值问题的渐近传播行为，未触及行波解本身的存在性、精细渐近结构及后端收敛性。本文旨在填补这一关键空白：（i）证明脉动行波（pulsating traveling waves）在任意 $c\ge c_R^*$ （右行）或 $c\ge c_L^*$ （左行）下的存在性；（ii）确立行波前沿（leading edge）的指数衰减率与线性预测完全一致，即 $u(t,x)\sim \alpha \varphi_{\lambda_c}(x)e^{-\lambda_c(x-ct)}$ 等；（iii）揭示 $c_R^*\neq c_L^*$ 的构造性例子及其多尺度奇异性机制；（iv）在齐次或快振荡系数下，证明行波后端收敛于唯一正稳态。这些结果将混合系统的行波理论从初值问题推进至解的内在几何结构层面，为理解空间异质性如何塑造病原体传播的定向性与不对称性提供严格数学基础。

2. 主要定理或主要结果

行波存在性（定理2.13）：设系数满足假设1且周期主特征值 $\lambda_1^{\text{per}}>0$ ，则右行脉动行波 $\big(u(t,x),v(t,x)\big)$ 存在当且仅当波速 $c\ge c_R^*$ ；左行脉动行波存在当且仅当 $c\ge c_L^*$ 。此处 $$c_R^*,c_L^*$$ 由 $k(\lambda)$ 定义为 $c_R^*=\min_{\lambda>0}k(\lambda)/\lambda$ , $c_L^*=\min_{\lambda<0}k(\lambda)/(-\lambda)$ 。特别地，当 $$c_R^*=0$$ （定理2.15），存在右定常波（stationary front）；若 $$c_R^*>0$$ ，则 $$c_R^*$$ 是最小可行波速，且临界波 $$c=c_R^*$$ 的存在性不依赖于超临界波列的极限过程 $c\to c_R^*$ ，而是通过直接构造专属上下解实现。
前沿渐近衰减（定理2.17）：对任意右行波，　　(i) 若 $$c>c_R^*$$ （超临界），令 $\lambda_c>0$ 为方程 $k(\lambda)=\lambda c$ 的最小正根， $(\varphi_{\lambda_c},\psi_{\lambda_c})$ 为其对应的 $\lambda_c$ -周期主特征向量，则存在 $\alpha>0$ 使得
$\lim_{x\to+\infty}\sup_{t\in[0,L/c]}\max\left\{\left|u(t,x)e^{\lambda_c(x-ct)}-\alpha\varphi_{\lambda_c}(x)\right|,\ \left|v(t,x)e^{\lambda_c(x-ct)}-\alpha\psi_{\lambda_c}(x)\right|\right\}=0.$
　　(ii) 若 $$c=c_R^*$$ （临界），令 $\lambda^*>0$ 为 $k(\lambda)=\lambda c_R^*$ 的唯一正根，则存在 $\alpha>0$ 使得
$\lim_{x\to+\infty}\sup_{t\in[0,L/c_R^*]}\max\left\{\left|\frac{u(t,x)e^{\lambda^*(x-c_R^*t)}}{x-c_R^*t}-\alpha\varphi_{\lambda^*}(x)\right|,\ \left|\frac{v(t,x)e^{\lambda^*(x-c_R^*t)}}{x-c_R^*t}-\alpha\psi_{\lambda^*}(x)\right|\right\}=0.$
即临界波呈准指数衰减（quasi-exponential decay），含代数修正因子 $$1/(x-c_R^*t)$$ ，符合线性色散关系预测。
后端收敛性（定理2.20）：当所有系数为常数时，设矩阵 $A=\begin{pmatrix}r_u-\mu_u & \mu_v\\ \mu_u & r_v-\mu_v\end{pmatrix}$ 的主特征值 $\lambda_A>0$ ，则系统存在唯一正稳态 $$(u^*,v^*)$$ ，且任意右行波 $\big(U(x-ct),V(x-ct)\big)$ 满足
$\lim_{\xi\to-\infty}\big(U(\xi),V(\xi)\big)=(u^*,v^*).$
对快振荡系数 $r_*^\varepsilon(x)=r_*^1(x/\varepsilon)$ 等，其同质化极限（2.30）的主特征值 $\bar{\lambda}_A>0$ 时，行波后端亦收敛于该极限系统的唯一正稳态。
不对称传播速度（定理2.22）：存在参数配置使 $c_R^*\neq c_L^*$ 。具体通过多尺度奇异极限构造：设 $r_u(x)=r_0+\varepsilon r_1(x)$ , $r_v(x)=r_0-\varepsilon r_1(x)$ , 其余系数适当选取，令 $\varepsilon\to0$ 并缩放变量，证明 $k(\lambda)$ 在 $\lambda>0$ 与 $\lambda<0$ 区域的极小值点不对称，从而 $c_R^*\neq c_L^*$ 。此例凸显系统与标量方程的本质差异。

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

本文证明的核心障碍在于系统缺乏全局比较原理，导致标准上下解法、单调迭代法及相平面分析均不可直接应用。作者发展了一套针对混合系统的精细化上下解构造框架，其创新点体现在三个层面：

第一，超临界波（ $$c>c_R^*$$ ）的上下解设计。上解取为线性化系统（2.4）的 $\lambda_c$ -指数模态解： $\overline{U}(t,x)=A\varphi_{\lambda_c}(x)e^{-\lambda_c(x-ct)}$ , $\overline{V}(t,x)=A\psi_{\lambda_c}(x)e^{-\lambda_c(x-ct)}$ ，其满足 $\mathcal{F}[\overline{U},\overline{V}]\le0$ （ $\mathcal{F}$ 为系统右端算子），因非线性项在前沿可被线性主导。下解构造更为精巧：借鉴Hamel[17]对标量方程的技巧，但避免简单取 $\max\{0,\text{sub-sol}\}$ （此操作在系统中会破坏耦合结构）。作者构造形如 $\underline{U}=a\varphi_{\lambda_c}e^{-\lambda_c(x-ct)}-b\varphi_{\tilde{\lambda}}e^{-\tilde{\lambda}(x-ct)}$ , $\underline{V}=a\psi_{\lambda_c}e^{-\lambda_c(x-ct)}-b\psi_{\tilde{\lambda}}e^{-\tilde{\lambda}(x-ct)}$ 的差分形式，其中 $\tilde{\lambda}>\lambda_c$ 确保第二项衰减更快，且系数 $$a,b$$ 的选择使 $\mathcal{F}[\underline{U},\underline{V}]\ge0$ 在 $$x$$ 充分大时成立。此构造利用了 $\lambda$ -周期主特征向量的正性及 $k(\lambda)$ 的凸性，是系统特有技术。

第二，临界波（ $$c=c_R^*$$ ）的专属上下解。此时 $\lambda^*$ 满足 $k'(\lambda^*)=c_R^*$ （由极小值条件），线性化算子在 $\lambda^*$ 处退化，标准指数上解失效。作者构造上解为 $\overline{U}=A\varphi_{\lambda^*}(x)(x-c_R^*t)^{-1}e^{-\lambda^*(x-c_R^*t)}$ , $\overline{V}=A\psi_{\lambda^*}(x)(x-c_R^*t)^{-1}e^{-\lambda^*(x-c_R^*t)}$ ，引入代数衰减因子以匹配临界色散关系；下解则采用 $\underline{U}=a\varphi_{\lambda^*}e^{-\lambda^*(x-c_R^*t)}-b\varphi_{\lambda^*}e^{-\lambda^*(x-c_R^*t)}(x-c_R^*t)^{-\delta}$ （ $\delta>0$ 小），通过精细估计非线性剩余项证明其为有效下解。此方法绕开了常规的 $c\to c_R^*$ 极限论证，是本文最原创的技术贡献。

第三，行波前沿渐近的刚性证明（定理2.17）。作者在移动坐标系 $\xi=x-ct$ 下将行波方程重写为椭圆型系统，并引入加权空间 $X_\lambda=\{w\in C^0_{\text{per}}([0,L];\mathbb{R}^2):\|w\|_\lambda=\sup_{\xi\ge0}|w(\xi)|e^{\lambda\xi}<\infty\}$ 。对超临界情形，定义算子 $\mathcal{T}:X_{\lambda_c}\to X_{\lambda_c}$ 为线性化算子的逆，证明其为压缩映射，从而解的渐近展开唯一确定；对临界情形，需转入加权空间 $Y_{\lambda^*}=\{w:\|w\|_{Y}=\sup_{\xi\ge1}|w(\xi)|\xi e^{\lambda^*\xi}<\infty\}$ ，并利用 $k''(\lambda^*)>0$ （凸性）控制高阶项。整个框架依赖于 $\lambda$ -周期主特征问题的谱理论（Krein-Rutman定理）、 $k(\lambda)$ 的严格凸性（命题2.3(iii)）及线性化算子在加权空间中的Fredholm性质，构成一套完整的系统行波渐近分析工具箱。

4. 与之前工作的比较

本文与既有文献的关键差异体现在四个维度：

与标量周期反应-扩散方程的比较：对单一方程 $u_t=(\sigma(x)u_x)_x+f(x,u)$ ，Weinberger[32]、Berestycki-Hamel-Roques[5]等已建立完备的脉动行波理论，其最小波速 $$c^*$$ 由 $\min_{\lambda}k(\lambda)/\lambda$ 给出，且左右速度恒等（因Fredholm替代性保证 $k(\lambda)$ 为偶函数）。本文首次在非标量框架下证明 $$c_R^*$$ 与 $$c_L^*$$ 可严格不等（定理2.22），并通过多尺度奇异极限显式构造反例，揭示了系统耦合结构对传播各向异性的根本影响——这是标量理论无法企及的深度。
与纯合作/纯竞争系统的比较：Wang[28]、Girardin[12]等对齐次系数的混合系统建立了行波存在性，但其证明严重依赖比较原理在局部区域的有效性。本文将此推广至空间周期情形，核心突破在于：当比较原理全局失效时，通过 $\lambda$ -周期主特征向量的正性与 $k(\lambda)$ 的凸性，构造出适配系统混合结构的上下解（见方法部分），而非强行恢复比较性。代价是证明复杂度剧增，且无法直接继承单调性结论（故定理2.19需额外小参数条件 $\kappa_{\max}/\kappa_{\min}\le\mu_{\min}/r_{\max}$ 才能证唯一性）。
与前作[15]的比较：[15]解决了初值问题的传播速度刻画（定理2.7），但未触及行波解本身。本文（定理2.13）证明传播速度 $$c_R^*,c_L^*$$ 恰为行波存在的最小波速，实现了“宏观传播”与“微观行波”两个层面的严格对应。此外，[15]仅知解在后端有界，本文（定理2.20）在齐次或快振荡条件下证明其收敛于唯一正稳态，补全了行波的整体结构图景。
与Hamel[17]的比较：Hamel对标量KPP方程证明了前沿指数衰减，其下解构造基于单个指数函数的差分。本文将此思想推广至系统，但面临两大挑战：(i) 差分下解需同时控制两个分量，且当一者为负时另一者仍需为正，故需双参数精细调节；(ii) 临界情形的代数修正因子 $$1/(x-ct)$$ 在系统中需与周期特征向量耦合，构造难度远超标量情形。本文的临界波直接存在性证明（非极限法）亦是对[17]技术的重要超越。综上，本文在假设上未强于[15]（仅要求 $\lambda_1^{\text{per}}>0$ ），但在结论深度（行波存在性、前沿刚性衰减、后端收敛、不对称速度）与技术原创性（混合系统上下解框架）上实现了系统性突破。

3. 非齐次抛物型方程的有限时间爆破

Finite time blow-up for an inhomogeneous parabolic equation

作者: Kaiqiang Zhang

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11165v1

摘要翻译: 我们考虑非齐次非线性热方程

Comment: 无

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

1. 研究对象和背景

本文研究三维空间中带空间非齐次源项的非线性热方程

\partial_t u - \Delta u = |u|^{p-1}u + f(x), \quad x \in \mathbb{R}^3, \; p > 5,

其中 $f \in L^\infty(\mathbb{R}^3) \cap C^{0,1}(\mathbb{R}^3)$ 是一个有界且 Lipschitz 连续的外力项。该方程是经典 Fujita 型非线性热方程 $\partial_t u - \Delta u = |u|^{p-1}u$ 的扰动版本，其关键特征在于：当 $f \not\equiv 0$ 时，方程完全丧失尺度不变性——即不存在 $\lambda > 0$ 使得若 $$ u(t,x) $$ 是解，则 $u_\lambda(t,x) = \lambda^{2/(p-1)} u(\lambda^2 t, \lambda x)$ 仍是解。这一结构性缺失直接破坏了自相似变量变换 $s = -\log(T-t), \; y = x/\sqrt{T-t}$ 后所得系统的时间自治性：在齐次情形下，该变换将原方程化为关于 $v(s,y) = (T-t)^{1/(p-1)} u(t,x)$ 的自治椭圆型动力系统

\partial_s v = \Delta v - \frac{1}{2} y \cdot \nabla v - \frac{1}{p-1} v + |v|^{p-1} v,

其平衡点即为自相似剖面 $\Phi_n$ ，满足

-\Delta \Phi_n + \Lambda \Phi_n - \Phi_n^p = 0, \quad \Lambda \Phi_n = \frac{2}{p-1} \Phi_n + y \cdot \nabla \Phi_n.

而在非齐次情形下，经相同变换后得到

\partial_s v = \Delta v - \frac{1}{2} y \cdot \nabla v - \frac{1}{p-1} v + |v|^{p-1} v + e^{-\frac{p}{p-1}s} f\big(e^{-s/2} y\big),

右端显含 $$ s $$ -依赖项 $e^{-ps/(p-1)} f(e^{-s/2} y)$ ，导致系统非自治。该扰动虽随 $s \to +\infty$ （即 $t \to T^-$ ）指数衰减，但其空间振荡结构与 $\Phi_n$ 的渐近行为 $\Phi_n(y) \sim |y|^{-2/(p-1)}$ 相互耦合，使得标准谱分析与模态控制失效。特别地， $$ f $$ 的 Lipschitz 连续性仅保证 $\|f(x_1)-f(x_2)\| \leq L|x_1-x_2|$ ，无法推出更高正则性或紧支性，故不能通过局部化或截断消去其影响；而 $L^\infty$ 有界性又不足以压制非线性项 $|v|^{p-1}v$ 在大 $$ |y| $$ 区域的主导作用。因此，核心困难在于：如何在缺乏尺度对称性的前提下，仍能构造出以 $\Phi_n$ 为极限剖面的爆破解，并刻画其稳定性结构。已有工作如 Bandle–Levine–Zhang [1] 仅给出爆破判据（如 $\int_{\mathbb{R}^3} f > 0$ 且 $$ p < 3 $$ 时全局爆破），而 [27,42–45] 仅估计寿命 $T \lesssim \|u_0\|_{L^\infty}^{-(p-1)}$ ，均未触及剖面收敛性、爆破点定位及初始数据流形的几何结构。本文旨在填补这一空白：建立首个严格证明的、具有精确自相似剖面收敛性的非齐次爆破解存在性与有限余维稳定性理论。

2. 主要定理或主要结果

定理 1.1（主结果）：设 $$ d = 3 $$ ， $$ p > 5 $$ ， $f \in L^\infty(\mathbb{R}^3) \cap C^{0,1}(\mathbb{R}^3)$ 。则对任意大整数 $n \gg 1$ ，任意固定 $x_0 \in \mathbb{R}^3$ ，及充分小 $\lambda_0 > 0$ ，存在一个余维数为 $$ n $$ 的 Lipschitz 流形 $M_{n,\lambda_0,x_0} \subset L^\infty(\mathbb{R}^3)$ ，其元素均为非径向初值，使得对任意 $u_0 \in M_{n,\lambda_0,x_0}$ ，对应解 $$ u(t,x) $$ 在有限时间 $$ T $$ 爆破，且满足

\frac{1}{8\lambda_0^2} \leq T \leq 2\lambda_0^2.

更进一步，存在调制参数 $x(t) \in \mathbb{R}^3$ 及余项 $\bar{u}(t,y)$ ，使得

u(t,x) = \frac{1}{[2(T-t)]^{1/(p-1)}} \left[ \Phi_n(y) + \bar{u}(t,y) \right], \quad y = \frac{x - x(t)}{\sqrt{2(T-t)}},

其中余项满足

\lim_{t \to T^-} \|\bar{u}(t)\|_{L^\infty(\mathbb{R}^3)} = 0,

且调制中心收敛于爆破点：

x(t) \to x(T), \quad t \to T^-.

此外，爆破时间映射 $T: M_{n,\lambda_0,x_0} \to \mathbb{R}$ 关于 $L^\infty$ 拓扑是 Lipschitz 连续的：

|T(u_0) - T(u_0')| \lesssim \|u_0 - u_0'\|_{L^\infty(\mathbb{R}^3)}.

流形构造细节：令 $L_n = -\Delta + \Lambda - p \Phi_n^{p-1}$ 为围绕 $\Phi_n$ 的线性化算子，其在加权空间 $L^2_\rho$ （权函数 $\rho(y) = e^{-|y|^2/2}$ ）中具有谱分解：存在 $$ n+1 $$ 个负特征值 $-\mu_{n+1,n} < \cdots < -\mu_{2,n} < -\mu_{1,n} = -2$ ，对应球对称特征函数 $\psi_{j,n}$ ；另有三个零特征值对应平移模态 $\partial_k \Phi_n$ （ $$ k=1,2,3 $$ ）。定义正交补空间

V_n^\perp = \left\{ \bar{u} \in L^\infty(\mathbb{R}^3): \, (\bar{u}, \psi_{j,n})_\rho = 0, \; j = 2,\dots,n+1 \right\},

则流形表示为

M_{n,\lambda_0,x_0} = \left\{ \lambda_0^{-2/(p-1)} \left( \Phi_n + \bar{u}_0 + \sum_{j=2}^{n+1} a_j(\bar{u}_0) \psi_{j,n} \right)\!\left( \frac{x-x_0}{\lambda_0} \right): \, \bar{u}_0 \in V_n^\perp, \; \|\bar{u}_0\|_{L^\infty} \ll 1 \right\},

其中 $a_j: V_n^\perp \to \mathbb{R}$ 是 Lipschitz 函数。值得注意的是， $\psi_{1,n} = \Lambda \Phi_n / \|\Lambda \Phi_n\|_\rho$ 对应尺度变化模态， $\partial_k \Phi_n$ 对应位置移动模态，二者由调制参数 $\lambda(t), x(t)$ 动态控制，故不需施加初始正交条件；唯余 $$ n $$ 个不稳定模态 $\psi_{2,n},\dots,\psi_{n+1,n}$ 需通过余维 $$ n $$ 条件消除。

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

本文的核心技术障碍在于：非齐次项 $$ f(x) $$ 彻底破坏尺度不变性，导致自相似变量下的动力系统非自治，从而无法直接套用 Collot–Raphaël–Szeftel [6] 中基于自治流的有限余维稳定性框架。作者的真正创新点在于发展了一套“非自治扰动下的模态分离与误差吸收”机制，其关键步骤如下：

第一，引入动态调制与广义自相似坐标系。定义 $s = -\log(T-t)$ ， $y = (x - x(t))/\sqrt{2(T-t)}$ ，并设

v(s,y) = [2(T-t)]^{1/(p-1)} u(t,x) = \Phi_n(y) + \varepsilon(s,y),

其中 $\varepsilon$ 满足正交条件 $(\varepsilon, \psi_{j,n})_\rho = 0$ （ $j=1,\dots,n+1$ ）及 $(\varepsilon, \partial_k \Phi_n)_\rho = 0$ （ $$ k=1,2,3 $$ ）。此正交性由引理 3.1（几何分解）保证，其证明依赖隐函数定理作用于有限维投影映射 $$ G $$ ，关键在于利用 $$ L_n $$ 的自伴性及特征函数正交性 $(\psi_j, \psi_k)_\rho = \delta_{jk}$ ，确保 Jacobian 矩阵可逆。

第二，建立非自治项的精细衰减估计。将方程重写为

\partial_s \varepsilon = L_n \varepsilon + \mathcal{F}(s,y) + \mathcal{N}(\varepsilon),

其中非自治扰动项为

\mathcal{F}(s,y) = e^{-\frac{p}{p-1}s} f\big(x(s) + e^{-s/2} y\big) - \text{投影修正项},

非线性项 $\mathcal{N}(\varepsilon)$ 满足 $|\mathcal{N}(\varepsilon)| \lesssim |\varepsilon|^2 + |\varepsilon|^p$ 。利用 $f \in C^{0,1}$ ，可证

\|\mathcal{F}(s)\|_{L^2_\rho} \lesssim e^{-\frac{p}{p-1}s} \left(1 + e^{-s/2} \|y\|_{L^2_\rho}\right) \lesssim e^{-\frac{p}{p-1}s},

且该衰减速率 $e^{-ps/(p-1)}$ 远快于线性化算子最慢衰减模态 $e^{-2s}$ （因 $$ p>5 $$ 故 $$ p/(p-1) > 2 $$ ），从而 $\mathcal{F}$ 在整个 Bootstrap 区间内为低阶扰动。

第三，设计分层 Bootstrap 闭合方案。设定三组先验估计：(i) $L^\infty$ 控制 $\|\varepsilon(s)\|_{L^\infty} \leq e^{-\delta s}$ ；(ii) 能量控制 $\|\varepsilon(s)\|_{H^1_\rho} \leq e^{-\gamma s}$ （ $\gamma < 2$ ）；(iii) 调制参数控制 $|\partial_s \lambda| + |\partial_s x| \leq e^{-\kappa s}$ 。利用命题 2.2 的谱隙估计

(L_n \varepsilon, \varepsilon)_\rho \geq c_n \|\varepsilon\|_{H^1_\rho}^2 - \frac{1}{c_n} \sum_{j=1}^{n+1} (\varepsilon, \psi_{j,n})_\rho^2 - \frac{1}{c_n} \sum_{k=1}^3 (\varepsilon, \partial_k \Phi_n)_\rho^2,

结合正交条件消去不稳定模态投影项，再吸收 $\mathcal{F}$ 和 $\mathcal{N}$ 的贡献，最终闭合所有估计。特别地， $\mathcal{F}$ 的快速衰减允许选择足够大的初始 $s_0 = -\log(T)$ （即足够小的 $\lambda_0$ ），使扰动在整个演化过程中始终可控，且无需对 $\|f\|_{L^\infty}$ 施加小性假设。

4. 与之前工作的比较

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

对称性假设的彻底放宽：此前所有自相似爆破解构造（如 [6,7,15,24]）均严格依赖方程的尺度不变性，以此保障自相似变量下系统的自治性，进而将问题归约为围绕稳态 $\Phi_n$ 的非线性稳定分析。本文首次在完全丧失尺度不变性的框架下实现相同目标，其代价是引入非自治扰动 $\mathcal{F}(s,y)$ 并发展相应的误差吸收技术。相较而言，[40,47] 中处理的类似问题仅给出形式推导或数值验证，而本文提供严格解析证明。
稳定性结构的重构：在齐次情形 [6] 中，余维数 $$ n $$ 源于消除 $$ n $$ 个球对称不稳定模态 $\psi_{2,n},\dots,\psi_{n+1,n}$ ；而平移与尺度模态 $\partial_k \Phi_n$ 、 $\Lambda \Phi_n$ 由对称性自然生成，无需初始约束。本文中，由于 $$ f(x) $$ 破坏平移与尺度对称性，这些模态不再对应守恒律，但作者通过动态调制参数 $x(t), \lambda(t)$ 将其纳入控制，而非剔除——这本质上将原本由对称性承担的‘自由度’转嫁为动力系统的反馈调节，从而维持了相同的余维数 $$ n $$ ，但稳定性流形的几何结构更复杂（依赖 $$ f $$ 的具体形式）。
适用范围与精度提升：以往工作 [1,27,42–45] 仅给出爆破存在性及寿命上界 $T \lesssim \|u_0\|_{L^\infty}^{-(p-1)}$ ，而本文首次确立剖面收敛性 $\|\bar{u}(t)\|_{L^\infty} \to 0$ 及爆破点定位 $x(t) \to x(T)$ ，且给出爆破时间 $$ T $$ 关于初值的 Lipschitz 连续性，这是对爆破动力学的定量化飞跃。此外，要求 $$ p>5 $$ 与 [6] 一致，确保 $\Phi_n$ 存在且谱隙 $$ c_n > 0 $$ ，但本文未扩展至临界指数 $$ p = 5 $$ 或亚临界情形，此为当前方法的局限。
技术代价与开放问题：为应对非自治性，本文必须假设 $f \in C^{0,1}$ （Lipschitz 连续），而 [1] 中仅需 $f \in L^\infty$ 即可得爆破判据；若 $$ f $$ 仅为 $L^\infty$ ，则 $\mathcal{F}(s,y)$ 的振荡可能破坏 $H^1_\rho$ 估计的闭合。此外，流形 $M_{n,\lambda_0,x_0}$ 的正则性仅证明为 Lipschitz，而 [28] 中阈值流形可达 $$ C^1 $$ 甚至更高，本文未解决 $$ a_j $$ 的高阶可微性，此为遗留问题。从提供内容中无法确认是否可将结果推广至 $d \neq 3$ 或 $p \leq 5$ ，作者亦未显式展开与 [14,25] 等其他非齐次模型的比较。

4. 非均匀不可压缩热传导流体在穿孔区域上的均匀化

Homogenization of a non-homogeneous incompressible heat-conducting fluid in perforated domains

作者: Yong Lu, Jiaojiao Pan, Luqi Wang

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11036v1

摘要翻译: 本文研究了三维非均匀不可压缩热传导流体在具有亚临界尺寸孔洞的穿孔区域上的均匀化问题，其中粘性系数与热传导系数均被假定为依赖于温度。孔洞直径为 $\varepsilon^α \ (α>3)$ 量级，其中 $\varepsilon > 0$ 是度量孔洞间相互距离的小参数。我们证明：当 $\varepsilon\to 0$ 时，速度、密度与温度的极限行为由无孔洞的均匀区域上原始系统所支配。

Comment: Homogenization; perforated domains; non-homogeneous incompressible heat- conducting fluid

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

1. 研究对象和背景

本文研究三维非均匀不可压缩导热流体在穿孔区域 $\Omega_\varepsilon = \Omega \setminus \bigcup_{k \in K_\varepsilon} T_{\varepsilon,k}$ 上的均匀化极限行为，其中 $\Omega \subset \mathbb{R}^3$ 是 $C^{2,\beta}$ 有界区域，穿孔 $T_{\varepsilon,k} = x_{\varepsilon,k} + a_\varepsilon T_0$ 是以 $a_\varepsilon = \varepsilon^\alpha$ （ $\alpha > 3$ ）为直径、中心间距为 $O(\varepsilon)$ 的周期性小障碍； $K_\varepsilon = \{ \varepsilon k \in \mathbb{Z}^3 \mid Q_{k,\varepsilon} \subset \Omega \}$ ， $Q_{k,\varepsilon} = \varepsilon (-\frac{1}{2}, \frac{1}{2})^3 + \varepsilon k$ 。该模型由如下非齐次 Navier-Stokes-Boussinesq 型系统刻画：

$$ \begin{cases} \partial_t \varrho_\varepsilon + \operatorname{div}x (\varrho\varepsilon u_\varepsilon) = 0, \ \operatorname{div}x u\varepsilon = 0, \ \partial_t (\varrho_\varepsilon u_\varepsilon) + \operatorname{div}x (\varrho\varepsilon u_\varepsilon \otimes u_\varepsilon) + \nabla_x P_\varepsilon = \operatorname{div}x S(\Theta\varepsilon, \nabla_x u_\varepsilon) - \Theta_\varepsilon \nabla_x F, \ S(\Theta_\varepsilon, \nabla_x u_\varepsilon) = \mu(\Theta_\varepsilon) \big( \nabla_x u_\varepsilon + \nabla_x^\top u_\varepsilon \big), \

\operatorname{div}x (\kappa\varepsilon \nabla_x \Theta_\varepsilon) = \nabla_x F \cdot u_\varepsilon, \end{cases} \quad \text{in } (0,T) \times \Omega_\varepsilon, $$

其中 $\varrho_\varepsilon(t,x) > 0$ 为质量密度， $u_\varepsilon(t,x) \in \mathbb{R}^3$ 为速度场， $\Theta_\varepsilon(t,x)$ 为温度；黏性系数 $\mu : \mathbb{R} \to [\underline{\mu}, \overline{\mu}]$ 是 Lipschitz 连续且正的函数，依赖于温度；热传导系数分段定义为 $\kappa_\varepsilon(x) = \kappa_f > 0$ （流体区 $\Omega_\varepsilon$ ）与 $\kappa_s > 0$ （固体障碍区 $\Omega \setminus \Omega_\varepsilon$ ），并设 $\kappa_f \neq \kappa_s$ ；外力项 $$F$$ 表示重力势，满足 $\nabla_x F \in L^\infty(\Omega; \mathbb{R}^3)$ ；边界条件为无滑移 $u_\varepsilon|_{\partial \Omega_\varepsilon} = 0$ 与齐次 Dirichlet 温度 $\Theta_\varepsilon|_{\partial \Omega} = 0$ ；初始数据 $(\varrho_0^\varepsilon, u_0^\varepsilon)$ 满足扩展后收敛性： $\widetilde{\varrho}_0^\varepsilon \to \varrho_0$ 强 $L^\infty(\Omega)$ ， $\widetilde{u}_0^\varepsilon \to u_0$ 强 $L^2(\Omega; \mathbb{R}^3)$ ，且 $0 < \underline{\varrho} \leq \widetilde{\varrho}_0^\varepsilon \leq \overline{\varrho}$ 。该系统可视为高 Péclet 数下经典热方程 $\partial_t(\varrho_\varepsilon \Theta_\varepsilon) + \operatorname{div}_x(\varrho_\varepsilon \Theta_\varepsilon u_\varepsilon) - \operatorname{div}_x(\kappa \nabla_x \Theta_\varepsilon) = \nabla_x F \cdot u_\varepsilon$ 的准静态近似（即忽略时间导数项 $\partial_t(\varrho_\varepsilon \Theta_\varepsilon)$ ）。问题核心在于：当 $\varepsilon \to 0$ 时，穿孔尺度 $a_\varepsilon = \varepsilon^\alpha$ （ $\alpha > 3$ ）属于亚临界（subcritical）情形——对应于经典均匀化参数 $\sigma_\varepsilon = \varepsilon^{\frac{3 - \alpha}{2}} \to +\infty$ ，此时穿孔体积占比 $|\Omega \setminus \Omega_\varepsilon| \sim \varepsilon^{3\alpha - 3} \to 0$ 极快，但其几何分布仍显著干扰弱解紧性。已有工作表明：在临界情形（ $\alpha = 3$ ，即 $\sigma_\varepsilon \to \sigma_* \in (0,\infty)$ ）下，Feireisl–Lu–Sun [12] 已证得极限系统为含摩擦项的 Brinkman 型方程；而在超临界（ $\alpha < 3$ ）下则退化为 Darcy 定律。然而，对亚临界非均匀导热流体，尚无关于温度依赖黏性与分段热导率耦合效应下的严格均匀化结果——尤其缺乏对非线性项 $\varrho_\varepsilon u_\varepsilon \otimes u_\varepsilon$ 、变系数应力张量 $S(\Theta_\varepsilon, \nabla_x u_\varepsilon)$ 及热源项 $\nabla_x F \cdot u_\varepsilon$ 在穿孔趋零极限下的一致强收敛控制。本文旨在填补这一空白，确立在 $\alpha > 3$ 下极限系统完全恢复原始方程（即无有效阻力项），并给出关键的强收敛速率估计。

2. 主要定理或主要结果

在假设 $\alpha > 3$ 、初始数据满足 (1.12)–(1.13)、黏性函数 $\mu \in \operatorname{Lip}(\mathbb{R})$ 且 $0 < \underline{\mu} \leq \mu(\cdot) \leq \overline{\mu} < \infty$ 、热导率 $\kappa_f, \kappa_s > 0$ 的条件下，设 $(\varrho_\varepsilon, u_\varepsilon, \Theta_\varepsilon)$ 为系统 (1.7)–(1.11) 的有限能量弱解，则存在子列使得零延拓速度 $\widetilde{u}_\varepsilon$ 、正延拓密度 $\widetilde{\varrho}_\varepsilon$ 与温度 $\Theta_\varepsilon$ 满足：

\widetilde{\varrho}_\varepsilon \to \varrho \quad \text{强 } C([0,T]; L^p(\Omega)) \quad \forall 1 \leq p < \infty,

\widetilde{u}_\varepsilon \to u \quad \text{强 } L^2((0,T) \times \Omega; \mathbb{R}^3) \quad \text{且弱 } L^2(0,T; W^{1,2}_0(\Omega; \mathbb{R}^3)),

\Theta_\varepsilon \to \Theta \quad \text{强 } L^2(0,T; W^{1,2}_0(\Omega)).

极限三元组 $(\varrho, u, \Theta)$ 是完整系统 (1.22)–(1.24) 的弱解，其动量方程右端无额外摩擦项：

\partial_t (\varrho u) + \operatorname{div}_x (\varrho u \otimes u) + \nabla_x P = \operatorname{div}_x \big[ \mu(\Theta)(\nabla_x u + \nabla_x^\top u) \big] + \Theta \nabla_x F.

更精确地，定理 1.2 给出余项估计：对任意 $\varphi \in C_c^\infty([0,T) \times \Omega; \mathbb{R}^3)$ 且 $\operatorname{div}_x \varphi = 0$ ，延拓动量方程满足

\int_0^T \int_\Omega \big[ \widetilde{\varrho}_\varepsilon \widetilde{u}_\varepsilon \cdot \partial_t \varphi + (\widetilde{\varrho}_\varepsilon \widetilde{u}_\varepsilon \otimes \widetilde{u}_\varepsilon) : \nabla_x \varphi \big] \,dx dt - \int_0^T \int_\Omega S(\Theta_\varepsilon, \nabla_x \widetilde{u}_\varepsilon) : \nabla_x \varphi \,dx dt = \int_0^T \int_\Omega \Theta_\varepsilon \nabla_x F \cdot \varphi \,dx dt - \int_\Omega \widetilde{\varrho}_0^\varepsilon \widetilde{u}_0^\varepsilon \cdot \varphi(0,\cdot) \,dx + \langle G_\varepsilon, \varphi \rangle,

其中余项满足 $|\langle G_\varepsilon, \varphi \rangle| \leq C \varepsilon^\sigma \big( \|\partial_t \varphi\|_{L^{4/3}(0,T; L^2(\Omega))} + \|\nabla_x \varphi\|_{L^4(0,T; L^{r_1}(\Omega))} + \|\varphi(0,\cdot)\|_{L^{r_2}(\Omega)} \big)$ ，其中 $\sigma = \frac{(3-q)\alpha - 3}{q} > 0$ ， $q \in (2,3)$ 接近 2， $r_1 \in (2,3)$ 满足 $$1/r_1 = 5/6 - 1/q$$ ， $$r_2 > 1$$ 满足 $$1/2 = 1/q + 1/r_2$$ 。该 $\varepsilon^\sigma$ 阶余项直接源于穿孔几何引起的 Bogovskii 算子误差（引理 3.1 中 $\|B_\varepsilon(f)\|_{W^{1,q}_0} \lesssim (1 + \varepsilon^{\frac{(3-q)\alpha - 3}{q}}) \|f\|_{L^q}$ ），是亚临界尺度下获得强收敛的关键量化指标。

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

证明的核心障碍在于：穿孔导致标准 Sobolev 嵌入与紧性失效，且非线性项 $\varrho_\varepsilon u_\varepsilon \otimes u_\varepsilon$ 、变系数应力 $S(\Theta_\varepsilon, \nabla_x u_\varepsilon)$ 及热源 $\nabla_x F \cdot u_\varepsilon$ 在 $\Omega_\varepsilon$ 上缺乏一致可积性。作者的真正创新点在于构建了一个适配亚临界穿孔尺度的“双层测试函数校正框架”，其技术枢纽包含三重精密分析：

第一，密度强收敛的 DiPerna-Lions 紧性升级。不同于经典 $L^\infty$ 密度下仅得弱- $\star$ 收敛，作者利用 (1.12) 中初始密度的常数延拓性质及 renormalized 方程 (1.20)，选取截断函数 $b(\widetilde{\varrho}_\varepsilon) = [\widetilde{\varrho}_\varepsilon - \overline{\varrho}]_+$ 与 $-[\widetilde{\varrho}_\varepsilon - \underline{\varrho}]_-$ ，结合能量不等式 (1.19) 推出一致上下界 $0 < \underline{\varrho} \leq \widetilde{\varrho}_\varepsilon \leq \overline{\varrho}$ ，从而激活 Lions [21] 中的紧性定理 2.4：当 $\partial_t \widetilde{\varrho}_\varepsilon + \operatorname{div}_x(\widetilde{\varrho}_\varepsilon \widetilde{u}_\varepsilon) = 0$ 且 $\widetilde{u}_\varepsilon$ 一致有界于 $L^2(0,T; W^{1,2}_0)$ 时， $\widetilde{\varrho}_\varepsilon$ 强收敛于 $$C([0,T]; L^p)$$ ，此为后续所有非线性项处理的基础。

第二，速度强收敛的 Bogovskii–Cut-off 耦合估计。由于直接延拓 $\widetilde{u}_\varepsilon$ 不满足 $\operatorname{div}_x \widetilde{u}_\varepsilon = 0$ ，作者引入穿孔依赖的截断函数族 $g_\varepsilon$ （满足 $g_\varepsilon = 0$ 于穿孔、 $g_\varepsilon = 1$ 于远离穿孔区域、 $|\nabla_x g_\varepsilon| \leq C \varepsilon^{-\alpha}$ ），将测试函数 $\varphi$ 分解为 $g_\varepsilon \varphi + (1-g_\varepsilon)\varphi$ 。关键突破在于： $g_\varepsilon \varphi$ 不满足散度自由，故不能代入原弱公式；作者调用穿孔域 Bogovskii 算子 $B_\varepsilon$ （引理 3.1），构造 $\varphi_2 = B_\varepsilon(\varphi \cdot \nabla_x g_\varepsilon)$ ，令 $\varphi_1 = g_\varepsilon \varphi - \varphi_2$ ，则 $\operatorname{div}_x \varphi_1 = 0$ 且 $\varphi_1$ 为合法测试函数。此构造将余项 $G_\varepsilon$ 精确分解为十项 $$I_1$$ – $I_{10}$ ，每一项均通过插值不等式、Sobolev 嵌入 $W^{1,q} \hookrightarrow L^{q^*}$ （ $$1/q^* = 1/q - 1/3$$ ）及 Bogovskii 算子范数估计 $\|B_\varepsilon(f)\|_{W^{1,q}_0} \leq C(1 + \varepsilon^{\frac{(3-q)\alpha - 3}{q}}) \|f\|_{L^q}$ 控制，最终导出 $\varepsilon^\sigma$ 阶余项（ $\sigma > 0$ ）。

第三，温度强收敛的椭圆正则性提升。热方程 $-\operatorname{div}_x(\kappa_\varepsilon \nabla_x \Theta_\varepsilon) = \nabla_x F \cdot u_\varepsilon$ 在穿孔上为变系数椭圆边值问题。作者利用 $\kappa_\varepsilon$ 的分段常数结构及 $u_\varepsilon$ 的 $$L^2$$ 一致有界性，结合 Poincaré 不等式 (2.6) 得 $\|\Theta_\varepsilon\|_{L^2(0,T; W^{1,2}_0)} \leq C$ ；再借助 Aubin-Lions 引理及 $\nabla_x F \in L^\infty$ ，从弱收敛 $\Theta_\varepsilon \rightharpoonup \Theta$ 推出强收敛 $\Theta_\varepsilon \to \Theta$ in $L^2(0,T; W^{1,2}_0)$ ，此为处理 $\mu(\Theta_\varepsilon)$ 非线性依赖的关键。

4. 与之前工作的比较

本文与既有工作的本质差异体现在四个维度：

穿孔尺度与极限类型：相比 Allaire [1,2] 对稳态 Stokes/Navier-Stokes 的经典三分法（临界 $\sigma_\varepsilon \to \sigma_*$ 导出 Brinkman，超临界 $\sigma_\varepsilon \to 0$ 导出 Darcy），本文首次将亚临界 $\sigma_\varepsilon \to +\infty$ 结果推广至非均匀导热流体。特别地，Feireisl–Lu–Sun [12] 处理了临界情形（ $\alpha = 3$ ），其极限含 Brinkman 摩擦项 $\sigma_*^2 u$ ；而本文 $\alpha > 3$ 使 $\sigma_\varepsilon$ 增长更快，余项 $G_\varepsilon$ 的 $\varepsilon^\sigma$ 阶衰减（ $\sigma > 0$ ）彻底消除有效阻力，极限严格回归原始系统——这是对 [12] 的互补与完成。
物理模型复杂性：此前亚临界结果（如 [6,9,11,27,33]）仅针对不可压缩 Navier-Stokes 或均匀密度流体，未涉及温度依赖黏性 $\mu(\Theta)$ 与分段热导率 $\kappa_\varepsilon$ 的耦合。本文必须同时控制 $S(\Theta_\varepsilon, \nabla_x u_\varepsilon)$ 的变系数振荡与热方程右端 $\nabla_x F \cdot u_\varepsilon$ 的弱收敛—强收敛混合项，这要求比 [27] 中纯流体情形更精细的 Bogovskii 算子估计（引理 3.1 中 $\varepsilon$ 幂次依赖 $\alpha$ 与 $$q$$ ）。
收敛强度与速率：Mikelić [30] 与 H"ofer–Nečasová–Oschmann [18] 在超临界下仅得弱收敛或定量 Darcy 收敛；本文则建立 $u_\varepsilon$ 的 $$L^2$$ 强收敛与 $\Theta_\varepsilon$ 的 $W^{1,2}$ 强收敛，并给出显式余项阶 $\varepsilon^\sigma$ （ $\sigma = \frac{(3-q)\alpha - 3}{q} > 0$ ），这是亚临界下首个具速率的导热流体均匀化结果。代价是要求 $\alpha > 3$ （而非 $\alpha > 2$ ），因 $\sigma > 0$ 需 $(3-q)\alpha > 3$ ，而 $$q > 2$$ 保证 Sobolev 嵌入有效。
技术工具拓展：Allaire [1] 使用振荡试函数法，但难以处理非线性变系数；本文的 Bogovskii–cut-off 框架（第 3.1 节）是对 [9,27] 中穿孔 Bogovskii 算子的实质性升级：前者仅用于稳态 Stokes，本文将其嵌入演化系统的动量方程推导，并与 renormalized 密度理论耦合，形成适用于非均匀导热流体的统一分析范式。从提供内容中无法确认是否对二维情形或更低正则性初始数据进行推广。

5. $L^2(\mathbb{R}^2)$ 适定性与密度块问题的对数Lipschitz正则性

$L^2(\mathbb{R}^2)$ Well-Posedness and Logarithmic Lipschitz Regularity for the Density Patch Problem

作者: Alessandro Violini

类别: math.AP

发布日期: 2026-07-12

arXiv ID: 2607.10676v1

摘要翻译: 我们研究二维非均匀不可压Navier–Stokes方程组在真空情形下的密度块问题，初始数据由一个Lipschitz密度块和一个属于 $L^2(\mathbb{R}^2)$ 的无散度速度场构成。

Comment: 无

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

1. 研究对象和背景

本文研究二维非均匀不可压Navier–Stokes方程在真空情形下的密度块（density patch）问题，即初始密度为指示函数 $\rho_0(x) = \mathbf{1}_D(x)$ 、初速度 $u_0 \in L^2_\sigma(\mathbb{R}^2)$ 的Cauchy问题。所考虑的方程组为

\begin{cases} \partial_t(\rho u) + \operatorname{div}(\rho u \otimes u) + \nabla P = \Delta u, \\ \partial_t \rho + \operatorname{div}(\rho u) = 0, \\ \operatorname{div} u = 0, \end{cases} \tag{1.1}

其中 $\rho = \rho(t,x) \ge 0$ 为密度， $u = u(t,x) \in \mathbb{R}^2$ 为速度场， $$P = P(t,x)$$ 为压力，粘性系数已归一化为1。初始数据满足

\rho(0,x) = \mathbf{1}_D(x), \quad u(0,x) = u_0(x) \in L^2_\sigma(\mathbb{R}^2), \tag{1.2}

其中 $D \subset \mathbb{R}^2$ 是有界Lipschitz区域， $L^2_\sigma(\mathbb{R}^2)$ 表示弱散度自由的 $$L^2$$ 向量场空间。该系统在物理上刻画一个由真空包围的均质不可压流体块的演化；其核心几何问题是：若初始界面 $\partial D$ 具有Lipschitz正则性，该正则性是否在时间演化中被保持？Lions在[20]中明确提出此问题，并证明存在全局Leray–Hopf弱解 $(\rho,u)$ 满足能量不等式

\frac{1}{2} \int_{\mathbb{R}^2} \rho |u(t)|^2 \,dx + \int_0^t \int_{\mathbb{R}^2} |\nabla u|^2 \,dx ds \le \frac{1}{2} \int_{\mathbb{R}^2} \rho_0 |u_0|^2 \,dx, \tag{1.3}

且密度始终保持块结构 $\rho(t) = \mathbf{1}_{D(t)}$ ，体积守恒 $$|D(t)| = |D|$$ 。然而，Lions未解决唯一性与界面正则性的定量传播问题。在严格正密度情形 $0 < c \le \rho_0(x) \le C$ 下，Hao–Shao–Wei–Zhang [16] 在 $u_0 \in L^2_\sigma$ 水平建立了全局适定性；但该方法严重依赖密度下界以控制动量方程的退化性。而在真空情形（即 $\rho_0$ 可取零值），方程退化为仅在支撑集 $\operatorname{supp}\rho$ 上定义的抛物型系统，导致经典Sobolev嵌入（如Ladyzhenskaya不等式 $\|v\|_{L^4} \lesssim \|v\|_{L^2}^{1/2} \|\nabla v\|_{L^2}^{1/2}$ ）无法直接应用于差值估计，因为加权范数 $\|\sqrt{\rho} \delta u\|_{L^2}$ 仅控制 $\delta u$ 在 $\operatorname{supp}\rho$ 上的行为。此前工作如Danchin–Mucha [11] 要求 $u_0 \in H^1$ ，[22] 推广至 $\mathbb{R}^2$ 但需额外几何假设；[23] 在临界Besov空间 $u_0 \in \dot{B}^0_{2,1}$ 下建立适定性，其关键在于利用 $\int_0^\infty \|\nabla u(t)\|_{L^\infty} dt < \infty$ 推出流映射一致Lipschitz连续，从而保持界面Lipschitz性。但该积分条件在 $$L^2$$ 初值下失效：热方程反例表明 $u_0 \in L^2$ 时 $\nabla u$ 未必属于 $L^1_{\mathrm{loc}}([0,\infty); L^\infty)$ 。因此，本文试图填补的核心空白是：在自然能量空间 $u_0 \in L^2_\sigma$ 下，建立密度块问题的全局唯一性（从而全局适定性），并发展适用于退化系统的弱正则性理论——特别是超越Lipschitz的对数Lipschitz估计——以定量刻画界面几何演化的稳定性。

2. 主要定理或主要结果

本文在两个层面取得严格证明的结果：

唯一性与全局适定性：设初始数据满足 $\rho_0(x) = \mathbf{1}_D(x)$ ， $u_0 \in L^2_\sigma(\mathbb{R}^2)$ ，其中 $D \subset \mathbb{R}^2$ 为有界Lipschitz区域。定义时间加权能量泛函

A_0(t,u) := \|\sqrt{\rho} u(t)\|_{L^2}^2 + \int_0^t \|\nabla u(s)\|_{L^2}^2 ds,

A_1(t,u) := t \|\nabla u(t)\|_{L^2}^2 + \int_0^t s \big( \|\nabla^2 u(s)\|_{L^2}^2 + \|\nabla P(s)\|_{L^2}^2 + \|\sqrt{\rho} \dot{u}(s)\|_{L^2}^2 \big) ds,

A_2(t,u) := t^2 \big( \|\sqrt{\rho} \dot{u}(t)\|_{L^2}^2 + \|\nabla^2 u(t)\|_{L^2}^2 \big) + \int_0^t s^2 \|\nabla \dot{u}(s)\|_{L^2}^2 ds,

A_3(t,u) := t^3 \|\nabla \dot{u}(t)\|_{L^2}^2 + \int_0^t s^3 \|\nabla^2 \dot{u}(s)\|_{L^2}^2 ds,

其中 $\dot{u} = \partial_t u + (u \cdot \nabla) u$ 为物质导数。称Leray–Hopf解 $(\rho,u)$ 为“即时强解”（immediately strong solution），若 $\sup_{i=0,1,2,3} \sup_{t>0} A_i(t,u) < \infty$ 。在此框架下，定理A断言：对任意满足(1.6)的初值，至多存在一个即时强解。结合[23, Theorem 1.2]的存在性结果，即得全局适定性：存在唯一即时强解 $(\rho,u)$ ，满足 $\rho u \in C([0,\infty); L^2)$ ， $\rho \in C([0,\infty); L^p_{\mathrm{loc}})$ （ $\forall p < \infty$ ），且能量恒等式(1.3)取等号。

对数Lipschitz正则性：定理B证明，在以下任一情形下，速度场 $$u$$ 满足对数Lipschitz估计：(i) $\rho$ 严格正（即 $0 < c \le \rho_0 \le C$ ）且 $(\rho,u)$ 为唯一Leray–Hopf解；(ii) $\rho_0 = \mathbf{1}_D$ 且 $(\rho,u)$ 为唯一即时强解。具体地，对任意 $\eta \in (0,1/2)$ ，存在 $\gamma \in L^1_{\mathrm{loc}}([0,\infty))$ ，使得对几乎处处 $$t > 0$$ 及所有 $x,y \in \mathbb{R}^2$ 满足 $$0 < |x-y| < 1$$ ，有

|u(t,x) - u(t,y)| \le \gamma(t) |x-y| \big(-\log|x-y|\big)^{1-\eta}. \tag{1.10}

作为推论，流映射 $X(t,\cdot)$ 属于 $L^\infty_{\mathrm{loc}}([0,\infty); C^{1-\varepsilon}(\mathbb{R}^2))$ 对任意 $\varepsilon \in (0,1)$ 成立，从而保证块边界Hausdorff维数守恒： $\dim_H(\partial D(t)) = 1$ ， $\forall t \ge 0$ 。

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

本文的证明包含两大核心技术障碍与对应创新：

障碍一：真空下差值估计的闭合问题。标准唯一性论证需控制能量差 $\frac{1}{2}\|\sqrt{\rho_1} \delta u(t)\|_{L^2}^2 + \int_0^t \|\nabla \delta u\|_{L^2}^2 ds$ ，其右端含两项： $I(t) = \left| \int_0^t \int_{\mathbb{R}^2} \rho_1 \delta u \otimes \delta u : \nabla u_2 \,dx ds \right|$ 和 $J(t) = \left| \int_0^t \int_{\mathbb{R}^2} \delta \rho \, \dot{u}_2 \cdot \delta u \,dx ds \right|$ 。在真空情形，经典Ladyzhenskaya不等式 $\|\sqrt{\rho_1} \delta u\|_{L^4} \lesssim \|\sqrt{\rho_1} \delta u\|_{L^2}^{1/2} \|\nabla \delta u\|_{L^2}^{1/2}$ 失效，因 $\|\sqrt{\rho_1} \delta u\|_{L^2}$ 不控制 $\delta u$ 在真空区的行为。作者的创新在于构造局部化Gagliardo–Nirenberg型估计：利用流映射的抛物位移性质 $|X_i(t,x) - x| \le C \sqrt{t}$ （ $$i=1,2$$ ），证明引理2.1：对任意 $f \in \dot{H}^1 \cap L^2(D)$ ，有

\|f\|_{L^4(X(t,D))} \le C \big( t^{1/4} \|\nabla f\|_{L^2} + \|f\|_{L^2(D)}^{1/2} \|\nabla f\|_{L^2}^{1/2} + \|f\|_{L^2(D)} \big).

此估计将 $$L^4$$ 范数限制在演化块 $$X_i(t,D)$$ 上，且不依赖 $\nabla u_i \in L^1_t L^\infty_x$ 。应用此式于 $f = \delta u$ ，可控制 $$I(t)$$ 。

障碍二： $\delta \rho$ 项的精细处理。 $$J(t)$$ 涉及密度差 $\delta \rho = \rho_1 - \rho_2$ ，其演化由 $\partial_t \delta \rho + \operatorname{div}(\delta \rho u_2) = -\operatorname{div}(\rho_1 \delta u)$ 决定，故需估计 $\delta u$ 在复合流 $Y_\tau(s,x) = X_2(s; \tau, X_1(\tau,x))$ 作用下的集合上的范数。关键难点在于 $$X_2$$ 非Lipschitz，传统链式法则失效。作者的突破性工具是加权流梯度控制：通过Proposition 2.4证明，对即时强解，有 $\|\nabla X_2(t,s,\cdot)\|_{L^\infty} \le w_M(t,s) := \exp\big(M \sqrt{\log(t/s)}\big)$ ，其中 $$M$$ 依赖于 $$C_u$$ 。此权重源于 $\int_0^\infty s \|\nabla u_2(s)\|_{L^\infty}^2 ds < \infty$ 的隐含信息，并经Cauchy–Schwarz导出。进而，Corollary 2.5给出 $\|\nabla \phi_i(s)\|_{L^p} \le w_{M_2}(t,s) \|\nabla g_i\|_{L^p}$ （ $\phi_i$ 为后向输运解），使 $$J(t)$$ 可分解为 $\int_0^t \int_{\mathbb{R}^2} \rho_1 \delta u \cdot \nabla(\phi_1 \phi_2)$ 并利用前述局部估计闭合Gronwall不等式。

障碍三： $$L^2$$ 初值下的对数正则性。Chemin–Lerner [5] 的Fourier方法在退化密度下失效。作者采用原子分解技术：利用插值 $L^2 = [\dot{H}^{-s}, \dot{H}^{s}]_{1/2,2}$ ( $s \in (0,1]$ )，将 $u_0 = \sum_j u_{j,0}$ 分解为满足 $c_j = 2^{-j/2} \|u_{j,0}\|_{\dot{H}^s} + 2^{j/2} \|u_{j,0}\|_{\dot{H}^{-s}}$ 且 $\|c_j\|_{\ell^2} < \infty$ 的散度自由原子。对每个原子 $u_{j,0}$ ，线性化系统(1.19)的解 $$u_j$$ 满足抛物衰减(1.20)： $\|\nabla u_j(t)\|_{L^\infty} \le C \min\{ t^{s/2-1} \|u_{j,0}\|_{\dot{H}^s}, t^{-s/2-1} \|u_{j,0}\|_{\dot{H}^{-s}} \}$ 。通过选取截断尺度 $N \approx -\log_2|x-y|$ 并估计 $A(t) = \sum_j (1+|j|)^{\eta-1} \|\nabla u_j(t)\|_{L^\infty}$ 与 $B(t) = \sum_j 2^{|j|/(sp)} (1+|j|)^{\eta-1} \|\nabla u_j(t)\|_{L^p}$ 的 $L^1_{\mathrm{loc}}$ 性，最终导出(1.10)。

4. 与之前工作的比较

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

适定性框架的尺度推进：[16] 在 $u_0 \in H^1$ 下建立适定性，[22] 将空间放宽至 $\mathbb{R}^2$ 但保留 $$H^1$$ 正则性；[23] 将初值空间推进至临界Besov空间 $\dot{B}^0_{2,1}$ ，其核心是利用 $\int_0^\infty \|\nabla u\|_{L^\infty} dt < \infty$ 导出流映射一致Lipschitz性。本文将初值空间彻底降至自然能量空间 $L^2_\sigma$ ，这是该问题的最小正则性要求。代价是放弃Lipschitz正则性，转而建立更弱的对数Lipschitz估计(1.10)，这反映了 $$L^2$$ 数据固有的正则性极限。

唯一性证明的技术革新：此前工作如[23]依赖 $\nabla u \in L^1_t L^\infty_x$ 控制差值，而本文在 $\nabla u \notin L^1_t L^\infty_x$ 前提下，通过抛物位移估计(1.7)与加权流梯度 $$w_M(t,s)$$ 构造局部化Gagliardo–Nirenberg不等式，首次实现 $$L^2$$ 水平的唯一性闭合。该方法不依赖Fourier分析，具有更强的模型无关性。

几何正则性的量化精度：Danchin–Mucha [11] 在 $u_0 \in H^1$ 下证明 $X \in C^{1,\gamma}$ ，[15] 推广至 $$H^s$$ ( $s > \gamma$ )，[23] 达到 $\dot{B}^\gamma_{2,1}$ 。本文在 $$L^2$$ 下获得 $X \in L^\infty_t C^{1-\varepsilon}_x$ ，虽弱于Hölder连续，但足以保证Hausdorff维数守恒(1.12)。值得注意的是，当 $\eta \to 0^+$ 时，(1.13)退化为 $|X(t,x)-X(t,y)| \le |x-y|^{\alpha(t)}$ ，导致 $\dim_H(\partial D(t)) \le e^{\|\gamma\|_{L^1(0,t)}}$ ，即维数可能指数增长；而 $\eta > 0$ 的严格条件正是本文确保维数恒为1的关键。

方法论的普适性拓展：Chemin–Lerner [5] 的频率分解法受限于常系数算子与正密度假设；本文的原子分解法基于插值空间 $[\dot{H}^{-s}, \dot{H}^{s}]_{1/2,2}$ ，且衰减估计(1.20)已在[9,23]中验证于退化情形，表明该框架可推广至更广泛的非均匀流体模型。从提供内容中无法确认该方法在三维情形的可行性，亦未显式展开与Vishik [19] 原子分解工作的比较细节。

6. 演化边界时滞方程

Evolutionary boundary delay equations

作者: Bernhard Aigner

类别: math.AP

发布日期: 2026-07-11

arXiv ID: 2607.10307v1

摘要翻译: 近期关于具有状态依赖非齐次项的演化型偏微分方程的适定性结果被推广至一大类更广的问题，该类问题纳入了非自治材料行为。这一推广使得构建一个能够容纳与演化型偏微分方程相关的边界时滞问题的框架成为可能。所得结果允许对具有状态依赖时滞的 Dirichlet 型与 Neumann 型边界条件进行特设处理。通过引入扩展状态空间，可容纳在本构律中含非自治时滞、在强迫项中含状态依赖时滞的复杂边界条件。该方法为若干涉及状态依赖时滞的经典边界条件提供了首个系统性的适定性分析。该理论的可行性与普适性通过其在具 Dirichlet、Neumann、Robin、Wentzell-Robin 与 Leontovich 边界条件的抛物型和双曲型偏微分方程中的应用得以展示。

Comment: 无

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

1. 研究对象和背景

本文研究的是一类具有非自治性与状态依赖时滞的演化型边值问题，其核心对象是定义在指数加权空间 $L^2_\rho(\mathbb{R}; H)$ 上的抽象演化方程

\left( \partial_t M + N + A \right) u = F(\cdot, u(\cdot)), \tag{1}

其中 $H = H_0 \times H_1$ 是 Hilbert 乘积空间； $M, N \in \mathcal{L}(L^2_\rho(\mathbb{R}; H))$ 是非自治有界线性算子（满足演化算子结构，见定义 2.4）； $A: H \supseteq \mathrm{dom}(A) \to H$ 是纯空间算子，分两类：(i) 最大耗散扩张的斜对称算子 $\mathring{A} = \begin{pmatrix} 0 & \mathring{D} \\ \mathring{G} & 0 \end{pmatrix}$ ，其中 $\mathring{D}: H_1 \supseteq \mathrm{dom}(\mathring{D}) \to H_0$ , $\mathring{G}: H_0 \supseteq \mathrm{dom}(\mathring{G}) \to H_1$ 满足 $\mathring{D} \subseteq -\mathring{G}^*$ , $\mathring{G} \subseteq -\mathring{D}^*$ ，该设定对应 Dirichlet/Neumann 型边界条件的抽象边界数据空间框架；(ii) 斜自伴算子 $A = \begin{pmatrix} 0 & -C^* \\ C & 0 \end{pmatrix}$ ，其中 $C: H_0 \supseteq \mathrm{dom}(C) \to H_1$ 是闭稠定线性算子，生成抽象 div-grad 系统，适用于 Robin、Wentzell–Robin、Leontovich 等复杂边界条件。关键未知量为状态函数 $u: \mathbb{R} \to H$ ，其历史映射定义为 $u(t): [-h, 0] \to H$ , $s \mapsto u(t+s)$ ，即前向演化依赖于长度为 $$h>0$$ 的有限回溯区间上的状态轨迹。右端项 $$F(t, u(t))$$ 含状态依赖时滞：典型形式包括 $F(t, u(t)) = I_\rho G(\cdot, u(\cdot))(t)$ （卷积型）或 $F(t, u(t)) = G(t, I_\rho u(t))$ （预解型），其中 $I_\rho \varphi(t) = \int_{-\infty}^t e^{-\rho(t-s)} \varphi(s) \, ds$ 是加权 Volterra 积分算子， $G: [0,\infty) \times L^2(-h,0; H) \to H$ 满足几乎一致 Lipschitz 连续性（定义 2.10）与正则性保持性（定义 2.11）。边界时滞体现为：Dirichlet 条件 $u|_{\partial\Omega}(t) = g(t - \tau(u(t)))$ 、Neumann 条件 $\partial_\nu u(t) = h(t - \sigma(u(t)))$ 或更一般的 Leontovich 条件 $\alpha(t) u|_{\partial\Omega}(t) + \beta(t) \partial_\nu u(t) = f(t - \delta(u(t), t))$ ，其中延迟函数 $\tau, \sigma, \delta$ 显含当前状态 $$u(t)$$ 或时间 $$t$$ ，导致非线性、非局部、非自治耦合。此类问题在热传导（含时变热导率 $$a(t,x)$$ ）、电磁波传播（含时变介电常数 $\epsilon(t,x)$ 、磁导率 $\mu(t,x)$ 、电导率 $\sigma(t,x)$ ）及粘弹性材料建模中广泛出现，但经典理论（如半群方法、能量估计）难以处理状态依赖延迟与非自治材料律的联合效应。已有工作如 [3] 仅覆盖自治情形（ $$M,N$$ 与 $$t$$ 无关），且未系统处理边界条件中的时滞；而传统泛函微分方程理论（如 Hale & Lunel）局限于 ODE 或简单 PDE，缺乏对抽象演化框架下最大耗散算子 $$A$$ 与非自治 $$M,N$$ 的统一处理能力。本文旨在填补这一空白：建立首个系统性框架，证明在非自治演化方程 (1) 上叠加状态依赖边界时滞后，仍能保证局部存在唯一性，并覆盖 Dirichlet、Neumann、Robin、Wentzell–Robin、Leontovich 等五类经典边界条件。

2. 主要定理或主要结果

本文在非自治演化方程框架下，建立了含状态依赖边界时滞的边值问题的局部存在唯一性理论。主要结果分为两部分：首先，对一般非自治演化方程 (1)，作者推广了 [3] 中的自治时滞结果，得到两个关键定理：定理 2.14（卷积型右端项）：设 $$G$$ 满足假设 2.13(A)， $\Phi \in W^{1,\infty}(-h,0; H)$ 与延拓 $$Z$$ 满足假设 2.13(B)， $$A, M, M', N$$ 满足定理 2.6 及定理 2.9（ $$k=1,2$$ ）的条件，且 $f \in H^2_{\tilde{\rho}}(0,\infty; H)$ 满足 $$f(0)=0$$ 与 $G(0,\Phi) = -f'(0)$ ，则方程 $v = S_\rho \left[ f + I_\rho G(\cdot, \pi_\alpha((v+Z)(\cdot))) \right]$ 在 $H^2_\rho(0,T; H)$ 中存在唯一局部解 $$u = v + Z$$ ，其中 $$T > 0$$ 依赖于 $\|\Phi'\|_\infty$ 、Lipschitz 常数 $L_\alpha$ 及 coercivity 常数 $$c$$ ，具体尺度关系为 $\rho > \max\left\{ \frac{L_\alpha^2}{2c^2}, \tilde{\rho} \right\}$ 。定理 2.15（预解型右端项）：在相同假设下，若 $f \in H^1_{\tilde{\rho}}(0,\infty; H)$ 且 $f(0) = -G(0,\Phi)$ ，则方程 $v = S_\rho \left[ f + G(\cdot, I_\rho v + Z)(\cdot) \right]$ 在 $H^1_\rho(0,T; H)$ 中存在唯一局部解， $$T$$ 满足 $\|v' + Z'\|_{L^\infty(0,T)} < \alpha$ 。其次，针对边界条件，作者构建了统一框架：对类型 (i) 算子 $$A$$ ，将非齐次 Dirichlet/Neumann 边界条件 $u|_{\partial\Omega}(t) = g(t - \tau(u(t)))$ 归化为方程 (1) 中 $$F$$ 的特定构造，证明其解属于 $H^1_\rho(0,T; H_0) \times H^1_\rho(0,T; H_1)$ ；对类型 (ii) 算子 $$A$$ ，引入扩展状态空间 $\mathcal{H} = H_0 \times H_1 \times L^2(-h,0; H_0) \times L^2(-h,0; H_1)$ ，将 Robin 条件 $\alpha(t) u|_{\partial\Omega} + \beta(t) \partial_\nu u = f(t - \delta(u(t),t))$ 编码为扩展算子 $\mathcal{A}: \mathcal{H} \supseteq \mathrm{dom}(\mathcal{A}) \to \mathcal{H}$ 的最大耗散性，并证明 $\partial_t \mathcal{M} + \mathcal{N} + \mathcal{A}$ 在 $L^2_\rho(\mathbb{R}; \mathcal{H})$ 上可逆。应用层面，定理 2.15 直接推出：对非自治热方程 $\partial_t \theta - \nabla \cdot (a(t,x) \nabla \theta) = Q(t,x)$ ，若 $$a$$ 满足例 2.16 的 (i)-(iv) 条件（含 $\|a_t\|_{L^\infty} < \infty$ ），且边界条件为 $\theta|_{\partial\Omega}(t) = G(t, I_\rho \theta(t))$ ，则存在 $$T>0$$ 使得解 $\theta \in H^1(0,T; L^2(\Omega))$ ；对非自治 Maxwell 方程 $\partial_t (\epsilon E, \mu H) + (\sigma E, 0) + (0,-\nabla\times H, \nabla\times E, 0) = (j_0, 0)$ ，若 $\epsilon, \mu, \sigma$ 满足例 2.17 的 (i)-(v) 条件（含 $\|\epsilon_{tt}\|, \|\sigma_t\| < \infty$ ），且边界条件为 $E_\parallel|_{\partial\Omega}(t) = G_1(t, I_\rho E(t), I_\rho H(t))$ ， $\nu \times H|_{\partial\Omega}(t) = G_2(t, I_\rho E(t), I_\rho H(t))$ ，则解 $(E,H) \in H^1(0,T; L^2(\Omega)^3 \times L^2(\Omega)^3)$ 存在唯一。所有结果均给出显式 $$T$$ 的下界估计： $T \gtrsim \min\left\{ \left( \frac{c}{L_\alpha} \right)^2, \frac{1}{\|\Phi'\|_\infty + \|Z'\|_{L^\infty}} \right\}$ 。

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

证明的核心障碍在于：非自治性破坏了经典 Picard 定理所需的 $$C_0$$ -半群结构，而状态依赖时滞导致右端项 $$F$$ 非 Lipschitz（因 $u(t - \tau(u(t)))$ 对 $$u$$ 的 Fréchet 导数涉及 $\tau'$ 与链式法则，产生无界因子），且边界时滞使问题丧失标准抛物/双曲正则性。作者的真正创新点在于三重技术整合：第一，采用 Picard 演化方程理论（[29]）的非自治版本，将问题置于指数加权空间 $L^2_\rho(\mathbb{R}; H)$ 与 $H^k_\rho(\mathbb{R}; H)$ ，利用时间导数 $\partial_{t,\rho} = I_\rho^{-1}$ 的代数性质（定义 2.2）规避半群构造，转而依赖算子 $\partial_t M + N + A$ 的 coercivity 估计（定理 2.6, 2.9 的条件 (ii)）。第二，针对状态依赖时滞，设计“截断-投影”不动点框架：定义凸闭集 $V_\alpha = \{ \phi \in H^1(-h,0; H) \mid \|\phi'\|_{L^\infty} \leq \alpha \}$ （定义 2.10），并利用其上度量投影 $\pi_\alpha$ （引理 [3, lem. 3.3]）构造修正映射 $\Gamma_{\rho,\alpha}(v) = S_\rho f + S_\rho F(\cdot, \pi_\alpha((v+Z)(\cdot)))$ （公式 (4)），将非 Lipschitz 问题转化为 $V_\alpha$ 上的 Lipschitz 问题，再通过连续性论证恢复原始解（定理 2.14 证明末段）。第三，为处理边界时滞，提出“扩展状态空间”方法：对 div-grad 系统，将延迟历史 $u_t(\theta) = u(t+\theta)$ （ $\theta \in [-h,0]$ ）作为独立变量嵌入新 Hilbert 空间 $\mathcal{H}$ ，并将边界条件编码为 $\mathcal{H}$ 上的耗散算子 $\mathcal{A}$ ，从而将时滞边界问题转化为无时滞的抽象演化问题 $\partial_t \mathcal{M} + \mathcal{N} + \mathcal{A}$ ，其 coercivity 由原算子 $$A$$ 的耗散性与延迟算子 $S_\tau$ 的有界性（引理 2.19）共同保障。关键估计包括：对非自治热方程，利用 $\partial_t(a^{-1}\phi) = (a^{-1})_t \phi + a^{-1} \phi'$ 及 $\|(a^{-1})_t\| \leq \|a^{-1}\|^2 \|a_t\|$ 得到 $\mathrm{Re}\langle a^{-1}\phi, \phi \rangle_{H^1_\rho} \geq \left( \frac{c}{\|a\|^2} - \frac{\|a_t\|}{\rho c^2} \right) \|\phi\|^2_{H^1_\rho}$ （第 9 页），确保 $\rho$ 足够大时 coercivity 成立；对 Maxwell 方程，类似估计 $\mathrm{Re}\langle (\partial_t \epsilon + \sigma)\phi, \phi \rangle_{H^1_\rho} \geq \left( d - \frac{\|\epsilon_{tt}\| + \|\sigma'\|}{\rho} \right) \|\phi\|^2_{H^1_\rho}$ （第 10 页）。整个框架的枢纽是解算子 $S_\rho = (\partial_t M + N + A)^{-1}$ 的因果性（theorem 2.6）与正则性提升性质（theorem 2.9），它们保证 $\Gamma_{\rho,\alpha}$ 将 $H^k_\rho$ 映入自身，使压缩映射原理适用。

4. 与之前工作的比较

本文相较既有工作实现了三重实质性推进。首先，对比自治时滞理论 [3]：[3] 仅处理 $$M,N$$ 与 $$t$$ 无关的情形，其 coercivity 条件 (ii) 简化为 $\mathrm{Re}\langle (\partial_t M + N)\phi, \phi \rangle_{L^2_\rho} \geq c \|\phi\|^2_{L^2_\rho}$ ，而本文定理 2.6, 2.9 要求非自治版本 $\mathrm{Re}\langle (\partial_{t,\rho} M^\rho + N^\rho)\phi, \phi \rangle_{H^k_\rho} \geq c \|\phi\|^2_{H^k_\rho}$ ，需额外验证 $$M', N$$ 的演化算子性质（definition 2.4）及 $M^\rho$ 在 $H^k_\rho$ 上的有界性（theorem 2.9 条件 (i)），这直接允许 $a(t,x), \epsilon(t,x)$ 等系数含时变导数（例 2.16, 2.17）。其次，对比传统泛函微分方程理论：Hale & Lunel [10] 的 ODE 框架要求 $\tau$ 为 Lipschitz 函数且 $\tau' < 1$ （assumption 2.18），但无法处理 PDE 的无穷维算子 $$A$$ ；而本文将 $\tau$ 纳入 $$G$$ 的参数，通过 $\pi_\alpha$ 投影控制 $\|u'\|_{L^\infty}$ ，将 ODE 的 $\tau' < 1$ 条件升维为 PDE 解的 $W^{1,\infty}$ 正则性要求，适用范围从有限维推广至 $L^2(\Omega)$ 等无限维空间。第三，对比边界条件研究：[5] 的边界三元组（boundary triples）理论仅适用于自治 $$A$$ 的缺陷指数相等情况，而本文类型 (ii) 的扩展状态空间方法绕过缺陷指数限制，直接构造耗散算子 $\mathcal{A}$ ，首次系统处理 Wentzell–Robin（含二阶边界导数）与 Leontovich（含阻抗边界）等强耦合时滞条件。代价方面：本文要求初始历史 $\Phi \in W^{1,\infty}(-h,0; H)$ （而非 $$H^1$$ ），且解正则性为 $$H^1$$ 或 $$H^2$$ （非更高阶），这是为控制 $\|u'\|_{L^\infty}$ 所必需；同时，coercivity 估计中的 $\rho$ 下界依赖于 $\|a_t\|, \|\epsilon_{tt}\|$ 等时变系数范数，当系数剧烈振荡时 $$T$$ 可能极小。局限在于：未处理分布时滞（distributed delay）或无限时滞（infinite delay），且对超临界非线性（如 $\tau(u) = |u|^{p}$ , $$p>1$$ ）未给出适定性结果——作者未显式展开比较，从提供内容中无法确认是否涵盖此类情形。

7. 一个与平衡律耦合的输运方程的存在性与唯一性

Existence and Uniqueness of a Transport Equation coupled with a Balance law

作者: Shuvam Banerjee, Manas R. Sahoo

类别: math.AP

发布日期: 2026-07-11

arXiv ID: 2607.10158v1

摘要翻译: 本文研究了一个具有不连续系数的输运方程的存在性与唯一性，该输运方程源于含时变源项的平衡律 $\partial_t u + \mathrm{div}(f(u)) = S(t,x,u)$ 的解。我们还给出一个简化证明，利用变分公式为该平衡律提供一个显式解公式。该系统可视为描述宇宙大尺度结构形成的经典一维模型的推广形式。

Comment: 22 pages

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

1. 研究对象和背景

本文研究如下一维非齐次双曲耦合系统：

\begin{cases} u_t + f(u)_x = \alpha(t) u, & x \in \mathbb{R},\ t > 0, \\ \rho_t + (f'(u) \rho)_x = 0, & x \in \mathbb{R},\ t > 0, \end{cases} \tag{1.3}

其中 $f \in C^2(\mathbb{R}; \mathbb{R})$ 严格凸且具有超线性增长： $\lim_{|u| \to \infty} f(u)/|u| = +\infty$ ； $\alpha \in L^\infty([0,\infty); \mathbb{R})$ ；初值 $u_0, \rho_0$ 为有界可测函数。该系统是零压气体动力学模型（Zeldovich 模型）的广义形式，其物理背景源于宇宙大尺度结构形成的流体力学建模：当压力可忽略时，物质粒子在引力作用下沿特征线自由运动，速度场 $$u$$ 满足含时变源项的守恒律，而密度 $\rho$ 沿 $$u$$ -依赖的特征方向满足输运方程。特别地，当 $\alpha(t) \equiv 0$ 且 $$f(u) = u^2/2$$ 时，(1.3) 退化为经典一维零压方程组

\begin{cases} u_t + (u^2/2)_x = 0, \\ \rho_t + (u \rho)_x = 0. \end{cases}

该模型在数学上极具挑战性：第一方程是非线性双曲守恒律带线性源项，其解通常在有限时间内形成激波并出现 Dirac 测度型奇性；第二方程则为系数 $$f'(u)$$ 不连续（因 $$u$$ 含激波）的输运方程，其解 $\rho$ 一般仅为 Radon 测度或 BV 函数的导数（即 $\rho = R_x$ ，其中 $R(x,t) = \int_0^x \rho_0(y)\,dy$ ）。已有工作如 Adimurthi–Ghoshal–Veerasamy [1] 基于 Hamilton–Jacobi 方程 $U_t + f(U_x) = \alpha(t) U$ 的 Lax–Oleinik 型变分公式，通过动态规划原理构造弱解，但要求 $$U$$ 具有 Lipschitz 连续性以保证弱导数存在；LeFloch [11] 则借助 Volpert 乘积处理不连续系数下的弱解定义，但未给出显式构造。本文试图填补的关键空白在于：(i) 对含时变源项 $\alpha(t)u$ 的广义平衡律，建立不依赖动态规划原理、直接从 Legendre 变换出发的变分公式；(ii) 在 $$u$$ 仅属 $L^\infty_{\mathrm{loc}} \cap \mathrm{BV}_{\mathrm{loc}}$ 的极端正则性下，为输运方程 $$R_t + f'(u) R_x = 0$$ （等价于 $\rho_t + (f'(u)\rho)_x = 0$ ）建立基于 Lebesgue–Stieltjes 积分的弱解理论，并证明其唯一性；(iii) 统一处理 $$u$$ 和 $\rho$ 的耦合奇性传播机制——即激波如何驱动密度集中（如 delta-shock），这在宇宙学中对应物质坍缩成星系团的数学刻画。

2. 主要定理或主要结果

本文在以下假设下建立了完整解理论： $f \in C^2(\mathbb{R})$ 严格凸、超线性； $\alpha \in L^\infty([0,\infty))$ ； $u_0, \rho_0 \in L^\infty(\mathbb{R})$ 。主要结果如下：

变分公式与 $$u$$ 的存在性（Theorem 3.1, Theorem 4.5）：定义 $\beta(t) := \int_0^t \alpha(s)\,ds$ ，并引入 $$h$$ -函数 $h : \mathbb{R} \times (0,\infty) \to \mathbb{R}$ ，其由 $$h$$ -曲线 $X'(\theta) = f'(h(x-y,t) e^{\beta(\theta)})$ ， $$X(0)=y$$ ， $$X(t)=x$$ 隐式定义。则对任意 $(x,t) \in \mathbb{R} \times (0,\infty)$ ，泛函
$F(x,y,t) := \int_0^t f^*\big(f'(h(x-y,t) e^{\beta(\theta)})\big) e^{-\beta(\theta)}\,d\theta + \int_0^y u_0(z)\,dz$
在 $y \in \mathbb{R}$ 上几乎处处达到唯一最小值 $$y^*(x,t) = y_*(x,t)$$ （Proposition 4.4(4)）。由此定义
$u(x,t) := h\big(x - y^*(x,t), t\big) e^{\beta(t)},$
则 $u = \nu$ 是方程 $u_t + f(u)_x = \alpha(t)u$ 的熵解，满足弱形式
$\int_0^\infty \int_{\mathbb{R}} \big[u \phi_t + f(u) \phi_x\big]\,dx dt + \int_{\mathbb{R}} u_0(x) \phi(x,0)\,dx = -\int_0^\infty \int_{\mathbb{R}} \alpha(t) u(x,t) \phi(x,t)\,dx dt,$
对所有 $\phi \in C_c^\infty(\mathbb{R} \times [0,\infty))$ 成立。
$\rho$ 的存在唯一性（Section 5）：设 $R(x,t) := \int_0^x \rho_0(z)\,dz$ ，则 $$R$$ 满足输运方程 $$R_t + f'(u) R_x = 0$$ 的弱解定义为：对任意 $\phi \in C_c^\infty(\mathbb{R} \times (0,\infty))$ ，
$\int_0^\infty \int_{\mathbb{R}} R \phi_t\,dx dt - \int_0^\infty \int_{\mathbb{R}} f'(u(x,t)) \phi(x,t)\,dR(x) dt = 0,$
其中 $$dR$$ 为 Lebesgue–Stieltjes 测度；且 $\|R(\cdot,t) - R_0\|_{L^\infty} \to 0$ 当 $t \to 0^+$ 。在此框架下，作者证明：若 $(u,R) \in L^\infty_{\mathrm{loc}}((0,\infty); \mathrm{BV}_{\mathrm{loc}}(\mathbb{R}; \mathbb{R}^2))$ 且映射 $t \mapsto \int_{\mathbb{R}} \phi(x) R(x,t)\,dx$ 连续，则弱解 $$R$$ 唯一（Theorem 5.1）。进而 $\rho = R_x$ （分布意义下）是原系统 (1.3) 的唯一解，且满足 $\rho(x,t) = \rho_0(y^*(x,t))$ （Theorem 3.1）。

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

本文的核心技术障碍在于：(i) 传统 Lax–Oleinik 公式依赖于 Hamilton–Jacobi 方程 $U_t + f(U_x) = \alpha(t)U$ 的 Lipschitz 连续解 $$U$$ ，但 $\alpha(t) \neq 0$ 时， $$U$$ 的弱导数 $$U_x = u$$ 的存在性无法由标准动态规划原理直接保证；(ii) 输运方程 $$R_t + f'(u) R_x = 0$$ 的系数 $$f'(u)$$ 在激波处不连续，导致经典特征线方法失效，且 Volpert 乘积难以与变分结构兼容。

作者的真正创新点在于三重构造：

$$h$$ -曲线与 $$h$$ -函数的几何重构：摒弃传统特征线 $$X'(s) = f'(u(X(s),s))$$ ，转而定义 $$h$$ -曲线为 $X'(\theta) = f'(h(x-y,t) e^{\beta(\theta)})$ ，其参数化显式嵌入了源项累积效应 $\beta(t)$ 。关键引理 Proposition 2.4 证明 $$h(x,t)$$ 关于 $$x$$ 严格递增且满足超线性增长 $\lim_{|x|\to\infty} \frac{1}{x} \int_0^x h(y,t)\,dy = \infty$ ，这确保了变分泛函 $$F(x,y,t)$$ 的强制性（coercivity），从而保证最小值存在。
免动态规划的变分潜力重构：Proposition 4.1 给出潜力 $$F$$ 的等价表示
$F(x,y,t) = F(x,0,t) + \int_0^y \big[u_0(z) - h(x-z,t)\big]\,dz,$
其核心在于利用 Legendre 变换恒等式 $$u f'(u) - f^*(f'(u)) = f(u)$$ 消去积分中的 $$f^*$$ 项，将问题转化为比较 $$u_0$$ 与 $$h$$ -函数的“斜率差”。这一重构使最小值点 $$y^*(x,t)$$ 的单调性（Proposition 4.4）和几乎处处唯一性（Lemma 4.3 的非交性）得以直接从 $$h$$ 的单调性导出，完全规避了动态规划原理。
Lebesgue–Stieltjes 弱解框架：针对输运方程，作者未采用 Sobolev 空间或测度解的传统定义，而是将 $$R$$ 视为 $\mathrm{BV}_{\mathrm{loc}}$ 函数，其弱解定义为
$\int_0^\infty \int_{\mathbb{R}} R \phi_t\,dx dt = \int_0^\infty \int_{\mathbb{R}} f'(u) \phi\,dR dt,$
其中右侧为 Lebesgue–Stieltjes 积分。证明存在性时，作者构造中间潜力 $G(x,y,t) = (x-y)h(x-y,t) - \int_0^t f(h(x-y,t) e^{\beta(s)}) e^{-\beta(s)}\,ds$ ，并证明 $\partial_t G = -f(h e^{\beta(t)}) e^{-\beta(t)}$ （式 (4.12)），从而获得 $$F$$ 关于 $$t$$ 的 Lipschitz 连续性（Step 2），这是验证弱解满足方程的前提。唯一性证明（Section 5）则依赖于 $$y^*, y_*$$ 的单调性及 $$h$$ -曲线的非交性，通过反证法排除两个不同解对应的 $$h$$ -曲线在内部相交的可能性。

4. 与之前工作的比较

本文与既有工作的实质性推进体现在四个维度：

源项处理的突破：Adimurthi et al. [1] 仅处理 $\alpha(t) \equiv 0$ 情形，其变分公式依赖于 $$U$$ 的 Lipschitz 连续性，而该性质在 $\alpha \neq 0$ 时需额外验证。本文通过引入 $\beta(t)$ 和 $$h$$ -函数，将源项完全吸收到特征线参数化中，使得 $u(x,t) = h(x-y^*,t) e^{\beta(t)}$ 的构造天然兼容 $\alpha \in L^\infty$ ，无需任何小初值或衰减假设。代价是 $$h$$ -函数的定义更复杂，但换来的是更普适的框架。
弱解定义的革新：LeFloch [11] 使用 Volpert 乘积定义 $f'(u) \rho$ ，要求 $$u$$ 为 BV 函数且 $\rho$ 为测度，但未给出显式公式。本文采用 Lebesgue–Stieltjes 积分定义 $$R$$ 的弱解，将 $\rho$ 的奇性完全编码于 $$R$$ 的跳跃集上，并通过 $$y^*(x,t)$$ 的显式表达 $\rho(x,t) = \rho_0(y^*(x,t))$ 实现解的完全构造。此方法避免了 Volpert 乘积的技术繁琐性，但要求 $R \in \mathrm{BV}_{\mathrm{loc}}$ ，限制了初始密度 $\rho_0$ 的正则性（仅需有界可测）。
时间尺度与适用范围：[1] 的解仅在短时间存在，因动态规划原理对 $$t$$ 有隐式限制；本文定理 4.5 和 5.1 的结论对所有 $$t > 0$$ 成立，得益于 $$h$$ -函数的全局定义和 $$F$$ 的强制性。然而，本文未讨论解的大时间行为（如自相似渐近），而 [8–10] 对 $\alpha=0$ 情形有深入研究。
代价与局限：本文的 $$h$$ -函数依赖于 $$f$$ 的严格凸性与超线性性，无法推广至非凸通量（如交通流模型）；且唯一性证明要求解属于 $L^\infty_{\mathrm{loc}} \cap \mathrm{BV}_{\mathrm{loc}}$ ，而 [1] 在更弱的 $L^\infty$ 框架下亦得唯一性。此外，作者未显式展开与 [2,3]（多维不连续系数输运方程）的比较，从提供内容中无法确认本文方法在高维的可推广性。

8. 具间断通量的粘性守恒律的存在性、唯一性与大时间行为

Existence, uniqueness and Large time behavior of a viscous conservation law with discontinuous flux

作者: Smritikana Pal, Manas Ranjan Sahoo

类别: math.AP

发布日期: 2026-07-11

arXiv ID: 2607.10219v1

摘要翻译: 本文研究一类具间断通量的粘性守恒律。我们定义了一种弱解概念，通过一个显式公式证明了该弱解的存在性，并证明了该弱解的唯一性。此外，我们得到了解的大时间行为。解的渐近极限是一个稳态解。此处，由对流项主导的大时间渐近行为或收敛于一个常数稳态解，或收敛于一个非常数稳态解，这取决于通量函数在阈值处（即两个通量函数的交点）的行为。这一现象在单通量情形下不会发生。

Comment: 26 pages, 12 figures

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

1. 研究对象和背景

本文研究如下一维具不连续通量的黏性守恒律初值问题：

\partial_t u + \partial_x \left[ (1 - H(x)) \frac{u^2}{2} + H(x) \frac{(u + \alpha)^2}{2} \right] - \partial_{xx} u = 0, \quad x \in \mathbb{R},\ t > 0,

u(x,0) = u_0(x), \quad x \in \mathbb{R},

其中 $$H(x)$$ 是 Heaviside 函数， $\alpha \neq 0$ 为实参数， $u_0 \in L^1_{\mathrm{loc}}(\mathbb{R})$ 。该方程是经典黏性 Burgers 方程在异质介质中的自然推广，其通量函数在 $$x=0$$ 处发生跳跃：左侧通量为 $$f(u) = u^2/2$$ ，右侧通量为 $g(u) = (u+\alpha)^2/2$ ，二者在唯一交点 $u = -\alpha/2$ 处相等——即满足所谓的“交汇条件”（meeting condition）。这一结构源于两相渗流模型（如带毛细扩散的 Buckley–Leverett 型方程），其中 $$x=0$$ 表征不同岩层界面， $\alpha$ 刻画渗透率或孔隙率的突变幅度，而 $-u_{xx}$ 项代表由毛细力、分子扩散或机械弥散引起的物理黏性效应。

从PDE理论角度看，该问题处于三个经典方向的交汇处：(i) 非线性抛物型方程的适定性；(ii) 不连续系数/非齐次通量下的弱解理论；(iii) 大时间渐近分析。已有工作对三类情形有较深入研究：(a) 光滑通量的黏性守恒律（如Hopf对Burgers方程的开创性工作[1]及后续[2–8]）；(b) 无黏性（ $\varepsilon=0$ ）不连续通量守恒律（见[12,13]，其核心挑战在于熵条件的选择与激波/接触间断的构造）；(c) 特殊源项（如 $\delta$ -函数源）的黏性Burgers方程（如Chung等[17]处理的 $u_t + u u_x = u_{xx} + \delta$ ，对应 $$f=g=u^2/2$$ 但含奇异源）。然而，兼具通量不连续性与黏性项的一般情形，尤其当左右通量在某状态值交汇时的大时间行为，尚未建立系统理论。关键空白在于：(1) 如何定义既兼容通量跳跃又体现黏性正则性的弱解概念；(2) 在交汇点 $u=-\alpha/2$ 处，通量导数 $$f'(u)=u$$ 与 $g'(u)=u+\alpha$ 符号相反，导致特征线在界面处呈现“入流-出流混合”结构，这使标准能量估计与最大值原理失效；(3) 大时间极限是否收敛至常数稳态，抑或出现非平凡稳态（如线性或分段线性分布），其判据完全依赖于 $\alpha$ 的符号与初始数据的整体积分性质，而既有文献（如[14–16]对Buckley–Leverett方程的研究）多假设通量单调或忽略黏性对稳态结构的决定性影响。本文旨在填补此空白，通过构造显式解公式，严格刻画交汇型不连续通量下黏性守恒律的全局适定性与稳态选择机制。

2. 主要定理或主要结果

本文在初始数据满足 $u_0 \in L^1((−∞, 0))$ 、 $u_0 + \alpha \in L^1((0, ∞))$ 及辅助函数 $\theta(\xi) := \phi_0(−\xi) + \phi_0(\xi)$ （其中 $\phi_0$ 由初值通过 $\phi_0(x) = \exp\left(−\frac{1}{2} \int_0^x u_0(y)\,dy\right)$ for $$x<0$$ 及 $\phi_0(x) = \exp\left(−\frac{1}{2} \int_0^x (u_0(y)+\alpha)\,dy\right)$ for $$x>0$$ 定义）满足 $\theta \in C^2(\mathbb{R}) \cap W^{2,\infty}(\mathbb{R})$ 的假设下，证明以下核心结论：

(1) 存在性与唯一性（Theorem 1.1 & Theorem 2.1）：存在唯一的弱解 $u \in C([0,\infty); H^1(\mathbb{R}) \cap L^\infty(\mathbb{R}))$ 满足积分恒等式

\int_0^\infty \int_{-\infty}^\infty \left[ u \zeta_t + \left( (1-H(x)) \frac{u^2}{2} + H(x) \frac{(u+\alpha)^2}{2} - u_x \right) \zeta_x \right] dx dt + \int_{-\infty}^\infty u_0(x) \zeta(x,0) dx = 0,

对任意 $\zeta \in C_c^\infty(\mathbb{R} \times [0,\infty))$ 成立。该解由显式公式给出：

u(x,t) = \begin{cases} -2 \dfrac{\phi_x(x,t)}{\phi(x,t)}, & x < 0,\\ -2 \dfrac{\psi_x(x,t)}{\psi(x,t)} - \alpha, & x > 0, \end{cases}

其中 $\phi, \psi$ 由热方程初边值问题（Lemma 2.1）及中间函数 $$G(t)$$ （满足Caputo半阶微分方程 $\mathcal{D}_t^{1/2} G(t) = \frac{\alpha}{4} G(t) + \frac{f(t)}{4}$ ，其解含Mittag-Leffler函数 $E_{1/2}, E_{1/2,1/2}$ ）显式表示（见(1.3)及(2.4)）。

(2) 大时间渐近行为：

当 $\alpha < 0$ 时（Theorem 1.2），解收敛至非恒定稳态： $\lim_{t \to \infty} u(x,t) = \begin{cases} \dfrac{-2\alpha \phi_0(-\infty)}{2(\phi_0(-\infty)+\phi_0(\infty)) + \alpha \phi_0(-\infty) x}, & x < 0,\\ \dfrac{\alpha (2\phi_0(-\infty) - \alpha \phi_0(\infty) x)}{\alpha \phi_0(\infty) x - 2(\phi_0(-\infty)+\phi_0(\infty))}, & x > 0. \end{cases}$ 该稳态在 $$x=0$$ 处连续（因黏性正则化），且依赖于初值在 $\pm\infty$ 处的指数积分极限 $\phi_0(\pm\infty)$ 。
当 $\alpha > 0$ 时（Theorem 1.3），解收敛至常数稳态： $\lim_{t \to \infty} u(x,t) = -\frac{\alpha}{2}, \quad \forall x \in \mathbb{R}.$ 此极限恰为左右通量 $$f(u), g(u)$$ 的唯一交点，且收敛速率受 $$G(t)$$ 的衰减控制（文中隐含 $G(t) \sim t^{-1/2} e^{-\alpha^2 t /16}$ ）。

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

本文的核心技术障碍在于：通量在 $$x=0$$ 处不连续，导致标准抛物型理论（如Schauder估计、Nash–Aronson热核估计）无法直接应用；同时，交汇条件 $f(-\alpha/2)=g(-\alpha/2)$ 虽保证通量连续，但其导数 $f'(-\alpha/2) = -\alpha/2$ 与 $g'(-\alpha/2) = \alpha/2$ 符号相反，使界面处特征线呈双向传播，破坏了单向能量耗散结构。作者的真正创新点在于构建了一套基于Hopf–Cole变换、热方程初边值分解与分数阶微分方程耦合的分析框架：

(1) 区域分解与Hopf–Cole线性化：将全空间问题分解为左半平面 ( $$x<0$$ ) 与右半平面 ( $$x>0$$ ) 上的两个热方程初边值问题（Lemma 2.1），其边界值 $$G(t)$$ 未知。通过定义 $U = -2\log\phi$ , $V = -2\log\psi - \alpha x$ ，将原非线性问题转化为 $$U,V$$ 满足的Burgers型方程（Problem 1 & 2），从而将通量不连续性“吸收”进边界条件，而非方程系数中。

(2) 界面匹配与Caputo半阶微分方程：为保证解在 $$x=0$$ 处的弱连续性（即 $$u_x$$ 的匹配），作者引入关键正则性条件 $G(t) = \frac{2}{\alpha} (\phi_x(0+,t) - \psi_x(0-,t))$ （Lemma 2.2）。将此条件代入热方程解的显式表达式后，经细致计算（步骤2），发现 $$G(t)$$ 满足Caputo型半阶微分方程 $\mathcal{D}_t^{1/2} G(t) = \frac{\alpha}{4} G(t) + \frac{f(t)}{4}$ ，其中 $$f(t)$$ 由初值 $\phi_0$ 决定。此方程的解由Mittag-Leffler函数 $E_{1/2}(z) = e^{z^2} \operatorname{erfc}(-z)$ 和 $E_{1/2,1/2}(z) = \frac{1}{\sqrt{\pi}} + z e^{z^2} \operatorname{erfc}(-z)$ 显式给出（(2.11)），其渐近行为 $G(t) \sim t^{-1/2} e^{-\alpha^2 t /16}$ 直接主导大时间尺度下的稳态选择。

(3) 稳态选择的符号判据：作者通过分析 $$G(t)$$ 的符号与衰减速率，揭示了 $\alpha$ 符号的物理意义：当 $\alpha > 0$ 时，界面处 $$f'$$ 为负、 $$g'$$ 为正，特征线“向外”发散，黏性 $-u_{xx}$ 主导并强制解均匀化至交汇点 $-\alpha/2$ ；当 $\alpha < 0$ 时，特征线“向内”汇聚，初值信息被“囚禁”在界面附近，导致非恒定稳态（Theorem 1.2中分式线性形式）的出现。此判据的严格证明依赖于对 $$G(t)$$ 渐近展开（Sections 3 & 4）与显式解公式中各项（如 $\phi_x/\phi$ ）在 $t\to\infty$ 时的逐点极限分析，其核心是利用 $\operatorname{erfc}(z)$ 的渐近展开 $\operatorname{erfc}(z) \sim \frac{e^{-z^2}}{\sqrt{\pi} z} (1 - \frac{1}{2z^2} + \cdots)$ 控制所有积分项。

4. 与之前工作的比较

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

(1) 与光滑通量文献（如Hopf[1], Liu[5], Huang et al.[11]）的比较：前述工作均假设通量 $$F(u)$$ 光滑且通常凸，其大时间渐近由初值的总质量与通量导数 $$F'(u)$$ 的全局单调性决定（如收敛至常数或稀疏波）。本文则首次系统处理通量本身在空间上不连续的情形，且当左右通量交汇时， $$F'(u)$$ 在界面处变号，导致稳态不再由全局凸性控制，而由局部符号 $\alpha$ 决定——这是光滑通量框架下不可能出现的现象（摘要明确指出“This phenomenon does not occur in the single-flux case”）。

(2) 与无黏性不连续通量文献（如Adimurthi et al.[12], Towers[13]）的比较：这些工作聚焦于熵解的存在唯一性及激波结构，其核心是设计合适的熵条件（如vanishing viscosity, adapted entropy）。本文则保留黏性项，并证明其不仅提供正则性，更从根本上改变大时间极限的定性行为：无黏性情形下，交汇型通量可能导致接触间断或非唯一熵解；而本文中黏性项强制解在 $$x=0$$ 处连续，并唯一确定稳态类型（常数或非常数），代价是需处理分数阶微分方程这一新工具。

(3) 与特殊黏性模型（如Chung et al.[17]）的比较：[17]研究 $u_t + u u_x = u_{xx} + \delta$ ，其本质是通量相同（ $$f=g=u^2/2$$ ）但含奇异源，故稳态由 $\delta$ 平衡，与交汇结构无关。本文的 $\alpha$ 参数直接调控通量交汇点的几何与动力学性质，使问题具有新的参数敏感性。此外，[17]未涉及 $\alpha$ 符号导致的稳态分支现象。

(4) 与Buckley–Leverett型研究（如[14–16]）的比较：这些工作虽涉及不连续通量，但多忽略黏性或仅作形式渐近，缺乏对黏性如何精确塑造稳态的严格分析。本文通过显式解与Caputo方程，首次给出了稳态的解析表达式（Theorem 1.2, 1.3），并量化了收敛速率（隐含于 $$G(t)$$ 的指数衰减）。其代价是假设初值满足较强的可积性与正则性（ $\theta \in C^2 \cap W^{2,\infty}$ ），而[14–16]常允许更广义的初值。从提供内容中无法确认本文结果是否可推广至更高维或更一般的交汇通量（如 $$f,g$$ 为多项式且有多重交点），作者亦未显式展开此类比较。

9. 带加性噪声的三维 Navier–Stokes 方程的误差分析

Error analysis for 3D Navier–Stokes equations with additive noise

作者: Soumya Ranjan Behera, Dominic Breit, Abhishek Chaudhary

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

发布日期: 2026-07-10

arXiv ID: 2607.09926v1

摘要翻译: 我们考虑带加性随机强迫项的三维随机 Navier–Stokes 方程。我们研究基于半隐式 Euler 格式以及 Crank-Nicolson 格式的时域离散化方法。我们证明：在局部时间范围内，前者以速率 1 收敛，后者以速率 3/2 收敛。这些理论结果得到了大量数值模拟的验证。据我们所知，这是首次对三维随机 Navier–Stokes 方程开展数值模拟。

Comment: 无

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

1. 研究对象和背景

本文研究三维周期域 $\mathbb{T}^3$ 上带加性噪声的随机Navier-Stokes方程（stochastic Navier-Stokes equations, sNSE）：

\begin{cases} du = \nu \Delta u\, dt - (u \cdot \nabla) u\, dt - \nabla p\, dt + \Phi\, dW & \text{in } O_T = (0,T) \times \mathbb{T}^3, \\ \mathrm{div}\, u = 0 & \text{in } O_T, \\ u(0) = u_0 & \text{in } \mathbb{T}^3, \end{cases} \quad \mathbb{P}\text{-a.s.}

其中 $\nu > 0$ 为运动粘性系数， $u : \Omega \times [0,T] \times \mathbb{T}^3 \to \mathbb{R}^3$ 为速度场， $$p$$ 为压力， $$W$$ 是定义在可分Hilbert空间 $$U$$ 上的柱状Wiener过程， $\Phi \in \mathcal{L}_2(U; H^2 \cap H)$ 是Hilbert-Schmidt型扩散算子， $u_0 \in L^2(\Omega; H^2(\mathbb{T}^3) \cap H)$ 是 $\mathcal{F}_0$ -可测初值。该系统在物理上建模湍流中不可预知的大尺度外源扰动，如热驱动、边界粗糙度或环境随机激励；数学上，加性噪声不依赖解本身，避免了乘性噪声引起的非线性耦合与局部Lipschitz障碍，但其路径仅具 $C^{\alpha}$ 正则性（对任意 $\alpha < 1/2$ ），严重制约时间离散精度。

关键难点在于三维sNSE缺乏全局强解理论：与确定性情形一致，仅知存在局部强路径解 $(u, \tau)$ ，其中停时 $\tau = \sup\{ t > 0 : \|u(t)\|_{H^2} < \infty \}$ 满足 $\mathbb{P}(\tau > 0) = 1$ ，且 $\mathbb{P}(\tau < \infty) > 0$ （可能爆破）。因此所有收敛性分析必须限制在停时 $\tau_R := \inf\{ t \geq 0 : \|u(t)\|_{H^2} \geq R \} \wedge T$ 上， $R \gg 1$ 。已有工作集中于二维情形：[7,13] 在乘性噪声下证得Euler格式在概率意义下收敛阶 $\alpha < 1/2$ ；[9] 利用Yosida正则化与小步长假设 $\tau \leq \tau_0(u_0, \Phi)$ ，在加性噪声下对二维sNSE证明Euler格式强收敛阶 $$1-$$ （即任意 $\varepsilon > 0$ 下 $O(\tau^{1-\varepsilon})$ ）；[1] 进一步通过随机PDE变换 $y(t) = u(t) - \Phi W(t)$ 将原问题转化为含确定性右端项的随机偏微分方程（random PDE）：

\partial_t y = \nu \Delta y - P_{\mathrm{HL}}[(y \cdot \nabla) y] + \nu \Delta[\Phi W] - P_{\mathrm{HL}}[L^W(y)], \quad \mathrm{div}\, y = 0,

其中 $P_{\mathrm{HL}}$ 为Helmholtz-Leray投影， $L^W(y) = (\nabla[\Phi W]) y + (\nabla y)[\Phi W] + (\nabla[\Phi W])[\Phi W]$ 。该变换使 $$y$$ 具有更高时间正则性：若 $\Phi \in \mathcal{L}_2(U; H^3)$ ，则 $\partial_t y \in L^2(0,T; V)$ 且 $\partial_t y \in C^{1/2}([0,T]; L^{r/2}(\Omega); V)$ （见引理5(c)），从而允许Crank-Nicolson格式达到 $$3/2-$$ 阶。然而，三维情形长期空白：(i) 所有已知强收敛结果（如[8]中 $$1/2-$$ 阶Euler估计）均未突破二维框架；(ii) 缺乏适用于局部强解的高阶时间离散误差分析框架；(iii) 无任何数值实验验证——因Monte Carlo模拟需大量样本与高分辨率三维网格，计算成本极高。本文旨在填补此三重空白，核心挑战在于：如何在停时 $\tau_R$ 上建立不依赖小步长假设的误差界，并处理三维非线性对流项 $(u \cdot \nabla) u$ 与随机项 $\Phi dW$ 的强耦合导致的高阶能量估计失效。

2. 主要定理或主要结果

本文在周期边界条件下，对三维加性噪声sNSE建立了两个时间离散格式的强收敛误差估计，所有结论均在停时 $\tau_R$ 上成立， $R \gg 1$ 。定理12：设 $u_0 \in L^r(\Omega; H^2 \cap H), \Phi \in \mathcal{L}_2(U; H^2 \cap H)$ ，对半隐式Euler格式 $$u^n$$ （定义见第3节），存在常数 $C = C(r,T,\Phi,u_0)$ 使得对任意 $\varepsilon > 0$ ，

\mathbb{E}\left[ \sup_{0 \leq n \leq M} \|u(t_n \wedge \tau_R) - u^n\|_{L^2_x}^r + \sum_{n=1}^M \tau \|\nabla(u(t_n \wedge \tau_R) - u^n)\|_{L^2_x}^2 \right]^{1/r} \leq C R^{6} \tau^{1-\varepsilon},

即强收敛阶为 $$1-$$ ，且无需传统小步长条件 $\tau \leq \tau_0$ 。定理14：设 $u_0 \in L^r(\Omega; H^4), \Phi \in \mathcal{L}_2(U; H^4)$ ，对Crank-Nicolson格式 $$u^n$$ （定义见第4节），存在常数 $C = C(r,T,\Phi,u_0)$ 使得对任意 $\varepsilon > 0$ ，

\mathbb{E}\left[ \sup_{0 \leq n \leq M} \|u(t_n \wedge \tau_R) - u^n\|_{L^2_x}^r + \sum_{n=1}^M \tau \|\nabla(u(t_n \wedge \tau_R) - u^n)\|_{L^2_x}^2 \right]^{1/r} \leq C R^{12} \tau^{3/2-\varepsilon},

即强收敛阶为 $$3/2-$$ 。定理18进一步给出压力近似 $$p^n$$ 的误差估计：在相同假设下，

\mathbb{E}\left[ \sum_{n=1}^M \tau \|p(t_n \wedge \tau_R) - p^n\|_{L^2_x}^2 \right]^{1/2} \leq C R^{6} \tau^{1-\varepsilon}.

所有常数 $$C$$ 显式依赖于 $T, r, \Phi, u_0$ ，但独立于 $$R$$ 和 $\tau$ 。数值验证（第5节）确认：在固定空间步长 $$h = 1/16$$ 下，IE格式动能误差随 $\tau$ 减小呈线性衰减，CN格式呈 $\tau^{3/2}$ 衰减，与理论一致。

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

本文的核心创新在于构建了一个适配三维局部强解的、基于随机PDE变换的路径依赖误差分析框架，克服了传统方法在停时上的技术断裂。第一，作者摒弃了[9]中依赖小步长 $\tau \leq \tau_0$ 的迭代能量估计，转而采用[1]启发的变换策略但进行根本性重构：定义 $y(t \wedge \tau_R) = u(t \wedge \tau_R) - \Phi W(t \wedge \tau_R)$ ，并严格证明（引理5）在 $\Phi \in \mathcal{L}_2(U; H^3)$ 下， $\partial_t y \in L^2(0,T; V)$ 且满足 $\mathbb{E}[\sup_{t \leq \tau_R} \|\partial_t y(t)\|_{V}^r] \leq C R^{6r}$ （推论6），这为高阶时间离散提供了关键正则性输入。第二，针对Euler格式（定理12），作者将误差 $e^n = u(t_n \wedge \tau_R) - u^n$ 分解为 $e^n = (y(t_n \wedge \tau_R) - y^n) + \Phi(W(t_n \wedge \tau_R) - W^n)$ ，其中 $$y^n$$ 是 $$y$$ 的Euler近似。关键突破在于：利用 $$y$$ 的 $C^{1,\alpha}$ 正则性（ $\alpha < 1/2$ ）和引理20的路径依赖Brownian积分误差界 $\|W(t_n) - W^n\|_{L^2_x} \leq C \tau^{1/2} \log^{1/2}(1/\tau)$ ，结合对流项线性化技巧 $$C(y,y,e^n) - C(y^n,y^n,e^n) = C(e^n,y,e^n) + C(y^n,e^n,e^n)$$ ，并通过Gagliardo-Nirenberg插值 $\|e^n\|_{L^4} \leq c \|e^n\|_{L^2}^{1/2} \|\nabla e^n\|_{L^2}^{1/2}$ 控制非线性项，最终导出 $\ell^\infty(L^2) \cap \ell^2(H^1)$ 型误差递推不等式，经Grönwall型引理（在停时上）得到 $O(\tau^{1-\varepsilon})$ 阶。第三，对Crank-Nicolson格式（定理14），障碍在于 $\partial_t y$ 仅 $C^{1/2}$ 连续，无法直接应用经典CN稳定性分析。作者引入创新性工具：引理21建立的路径依赖Brownian二次变差误差界 $\|\int_{t_{n-1}}^{t_n} (W(s) - W(t_{n-1})) ds\|_{L^2_x} \leq C \tau^{3/2} \log^{1/2}(1/\tau)$ ，该界比标准 $O(\tau^{3/2})$ 更优，补偿了三维非线性带来的额外损失；同时，将CN格式应用于变换后方程(9)的光滑部分 $\nu \Delta y$ 和 $P_{\mathrm{HL}}[(y \cdot \nabla) y]$ ，而将粗糙项 $\nu \Delta[\Phi W]$ 和 $P_{\mathrm{HL}}[L^W(y)]$ 以一阶Euler方式离散，形成混合离散策略。误差方程中出现的 $\partial_{tt} y$ 项由引理5(c)的 $C^{1/2}$ 估计控制，最终通过精细的能量配对与停时截断，获得 $O(\tau^{3/2-\varepsilon})$ 阶。第四，压力估计（定理18）依赖于离散inf-sup稳定性与速度误差的 $$H^1$$ 控制，结合连续压力的正则性 $\|p\|_{L^2(0,\tau_R; L^2)} \leq C R^3$ （由连续方程推导），实现 $O(\tau^{1-\varepsilon})$ 收敛。

4. 与之前工作的比较

本文与既有文献的比较需从三个维度展开：理论深度、技术适用性与数值实践。首先，相较于二维sNSE工作：[7,13] 仅获概率收敛阶 $$1/2-$$ ，且要求乘性噪声；[9] 在加性噪声下证得Euler阶 $$1-$$ ，但依赖小步长 $\tau \leq \tau_0(u_0,\Phi)$ ，该条件在实际计算中难以验证且限制自适应步长策略；本文定理12彻底消除此假设，通过停时上的统一能量估计实现更稳健的理论保障。其次，对比[1]的二维Crank-Nicolson结果：[1]在 $L^2(\Omega; L^2(0,T; H))$ 框架下证明 $$3/2-$$ 阶，其分析基于全局解存在性及均方意义下的误差控制；本文则首次将该高阶框架推广至三维，但面临根本性差异——三维解仅局部存在，故无法使用全局期望，必须在停时 $\tau_R$ 上进行路径依赖分析。为此，作者发展了新的技术：将[1]的均方估计替换为停时上的 $L^r(\Omega)$ 估计（ $r \geq 2$ ），并引入引理20–21的路径依赖Brownian积分界，这些界在 $\tau \to 0$ 时以 $\log^{1/2}(1/\tau)$ 因子弱于确定性情形，但足以支撑 $$3/2-$$ 阶。代价是常数显式依赖于 $$R$$ （如 $C R^{12}$ ），反映三维爆破风险对误差放大的量化影响，这是二维工作中不存在的固有局限。第三，与三维弱解数值工作[11]对比：[11]仅证明近似解在子列意义下弱收敛于martingale解，属概率测度层面的收敛，无法提供路径wise误差界；本文则首次建立强路径收敛，直接量化单条样本轨道的逼近精度，为不确定性量化（UQ）提供理论基础。最后，数值层面，本文是首篇报道三维sNSE Monte Carlo模拟的研究（第5节），其创新在于：(i) 采用MPI embarrassingly parallel架构，每个样本独立分配至MPI进程，规避了三维模拟的通信瓶颈；(ii) 使用Taylor-Hood $$P_2/P_1$$ 元在单位立方体上实现Dirichlet边界，突破论文理论中周期边界假设，验证算法的工程鲁棒性；(iii) 提供首套三维随机lid-driven cavity的系综平均可视化（vorticity contours, streamlines, iso-surfaces），揭示噪声强度 $\mu$ 对涡结构的调制机制——这在[5,6,25]等湍流建模文献中长期缺失。然而，本文未处理乘性噪声情形，亦未给出全局收敛性（因 $\mathbb{P}(\tau < \infty) > 0$ ），且空间离散误差未被分析，这些是未来工作的自然延伸。

10. 三维带双调和阻尼的罚 Navier–Stokes 系统的整体存在性与最优衰减

Global existence and optimal decay for a three-dimensional penalized Navier–Stokes system with biharmonic damping

作者: Kabiru Michael Adeyemo, Mohamed Majdoub, Subha Pal

类别: math.AP, math-ph

发布日期: 2026-07-12

arXiv ID: 2607.10657v1

摘要翻译: 我们研究一个三维抛物型系统，该系统作为不可压 Navier–Stokes 方程的超粘性与罚逼近而出现。该模型融合了三种互补的耗散机制：经典粘性扩散、双调和（超粘性）正则化以及散度罚项。此外，在非线性对流项中引入了一类 Temam 型校正项，以补偿罚过程所引发的弱可压效应。

Comment: 无

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

1. 研究对象和背景

本文研究如下三维抛物型系统：

\partial_t u - \nu \Delta u + N(u) - \frac{1}{\varepsilon} \nabla \mathrm{div}\, u + \beta \Delta^2 u = 0, \quad (t,x) \in \mathbb{R}_+ \times \mathbb{R}^3,

其中速度场 $u = u(t,x) \in \mathbb{R}^3$ ， $\nu > 0$ 为运动粘性系数， $\beta > 0$ 为超粘性（hyperviscosity）系数， $\varepsilon > 0$ 为罚参数，非线性项取 Temam 型斜对称形式

N(u) := (u \cdot \nabla) u + \frac{1}{2} u \, \mathrm{div}\, u.

该模型是不可压 Navier–Stokes 方程（NS）

\begin{cases} \partial_t u - \nu \Delta u + (u \cdot \nabla) u + \nabla p = 0, \\ \mathrm{div}\, u = 0 \end{cases}

的一种超粘性-罚逼近：它通过三项协同耗散机制替代原始 NS 的压力约束与弱非线性控制——经典粘性 $-\nu \Delta u$ 、双调和阻尼 $\beta \Delta^2 u$ （即四阶耗散）、以及发散罚项 $-\varepsilon^{-1} \nabla \mathrm{div}\, u$ 。Temam 校正项 $\frac{1}{2} u \, \mathrm{div}\, u$ 的引入至关重要：它使 $$N(u)$$ 满足 $$L^2$$ 意义下的无源性（cancellation），即对任意光滑向量场 $$u$$ ，有

\int_{\mathbb{R}^3} N(u) \cdot u \, dx = 0,

从而导出无条件能量恒等式

\frac{1}{2} \frac{d}{dt} \|u(t)\|_{L^2}^2 + \nu \|\nabla u(t)\|_{L^2}^2 + \beta \|\Delta u(t)\|_{L^2}^2 + \frac{1}{\varepsilon} \|\mathrm{div}\, u(t)\|_{L^2}^2 = 0.

此恒等式是整个分析的基石：它不依赖于 $\mathrm{div}\, u = 0$ ，且所有耗散项均具正号，尤其罚项 $\varepsilon^{-1} \|\mathrm{div}\, u\|_{L^2}^2$ 强制 $\mathrm{div}\, u \to 0$ 在 $L^2_{t,x}$ 意义下以速率 $O(\sqrt{\varepsilon})$ ；而双调和项 $\beta \|\Delta u\|_{L^2}^2$ 提供 $$H^2$$ 级先验控制，远超 Leray–Hopf 解的 $$H^1$$ 控制。

该模型的重要性在于其处于三个经典研究方向的交汇点：(i) 超粘性正则化（Lions [12] 证明 $(-\Delta)^\gamma$ 当 $\gamma \ge 5/4$ 可恢复三维整体正则性）；(ii) 罚方法（Temam [23] 用 $-\varepsilon^{-1} \nabla \mathrm{div}\, u$ 替代不可压约束，消除压力并简化数值实现）；(iii) 斜对称非线性（Temam 校正确保 $$L^2$$ 能量封闭，避免压缩性导致的能量注入）。然而，既有工作多聚焦单一机制：如 Lions 型超粘性未处理罚结构；Temam 罚法通常仅配经典粘性（ $-\nu \Delta u$ ），缺乏高阶耗散；而吸收型阻尼（如 $\alpha |u|^{r-1}u$ ）属零阶非线性修正，与双调和耗散的频域机制本质不同。本文填补的关键空白在于：在同时包含双调和耗散、罚项与斜对称非线性的耦合框架下，建立关于罚参数 $\varepsilon$ 一致的全局适定性与最优大时间衰减理论。这一一致性至关重要——它保证当 $\varepsilon \to 0$ 时，解族 $\{u^\varepsilon\}$ 在适当拓扑下具有紧性，为严格证明 $\varepsilon \to 0$ 极限收敛至不可压超粘性 NS 解（含压力恢复）提供坚实基础；而现有文献中，罚法与超粘性的联合分析几乎空白，尤其缺乏 $\varepsilon$ -一致估计。

2. 主要定理或主要结果

本文确立三项核心结果，全部具有关于罚参数 $\varepsilon > 0$ 的一致性：

整体弱解存在性（Theorem 1.3）：对任意初值 $u_0 \in L^2(\mathbb{R}^3; \mathbb{R}^3)$ ，存在全局弱解 $$u$$ 满足定义 1.1，即对任意 $$T > 0$$ ，
$u \in L^\infty(0,T; L^2) \cap L^2(0,T; H^2) \cap L^4(0,T; L^4), \quad \partial_t u \in L^{4/3}(0,T; H^{-2}), \quad \mathrm{div}\, u \in L^2(0,T; L^2),$
且满足弱形式 (1.6) 与能量不等式 (1.7)。关键在于所有先验估计（包括 $\|u\|_{L^\infty_t L^2_x}$ 、 $\|\nabla u\|_{L^2_{t,x}}$ 、 $\|\Delta u\|_{L^2_{t,x}}$ 、 $\varepsilon^{-1/2} \|\mathrm{div}\, u\|_{L^2_{t,x}}$ ）均独立于 $\varepsilon$ 。
一致小初值整体强解（Theorem 1.4）：存在仅依赖于 $\nu, \beta$ 的阈值 $\delta_0 > 0$ （与 $\varepsilon$ 无关），使得当 $u_0 \in H^2(\mathbb{R}^3; \mathbb{R}^3)$ 且 $\|u_0\|_{H^2} \le \delta_0$ 时，存在唯一全局强解
$u \in C([0,\infty); H^2) \cap L^2_{\mathrm{loc}}([0,\infty); H^4),$
并满足一致 $$H^2$$ 能量估计：
$\sup_{t \ge 0} \|u(t)\|_{H^2}^2 + c_0 \int_0^\infty \left( \|\nabla u(t)\|_{H^2}^2 + \|\Delta u(t)\|_{H^2}^2 + \frac{1}{\varepsilon} \|\mathrm{div}\, u(t)\|_{H^2}^2 \right) dt \le \|u_0\|_{H^2}^2,$
其中 $c_0 = c_0(\nu, \beta) > 0$ 亦与 $\varepsilon$ 无关。该结果首次在罚框架下实现 $$H^2$$ 级整体正则性，且小初值条件不随 $\varepsilon \to 0$ 恶化。
最优大时间衰减（Theorem 1.6）：存在 $\delta_* \in (0, \delta_0]$ （仅依赖 $\nu, \beta$ ），当 $u_0 \in L^1 \cap H^2$ 且 $\|u_0\|_{L^1} + \|u_0\|_{H^2} \le \delta_*$ 时，对应强解满足对所有 $t \ge 0$ ，
$\|u(t)\|_{L^2} \le C (1+t)^{-3/4}, \quad \|\nabla^k u(t)\|_{L^2} \le C (1+t)^{-3/4 - k/2}, \quad k = 1,2,$
其中常数 $$C > 0$$ 依赖于 $\nu, \beta$ 和 $\|u_0\|_{L^1 \cap H^2}$ ，但完全独立于 $\varepsilon$ 。这些速率与三维热方程 $e^{t\nu \Delta} u_0$ 对 $L^1 \cap L^2$ 初值的最优衰减率 $\|\nabla^k e^{t\nu \Delta} u_0\|_{L^2} \lesssim t^{-3/4 - k/2} \|u_0\|_{L^1}$ 完全一致，且被证明是 sharp 的（Remark 1.7）：低频主导衰减（ $|\xi| \lesssim t^{-1/2}$ ），此时双调和项 $\beta |\xi|^4$ 相对于 $\nu |\xi|^2$ 可忽略，故衰减速率由热流决定；而罚项虽影响高频，但不加速大时间行为。

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

证明的核心障碍源于三项耦合机制的相互作用：(i) 罚项 $-\varepsilon^{-1} \nabla \mathrm{div}\, u$ 在 $\varepsilon \to 0$ 时成为强奇异扰动，其弱形式中的 $\varepsilon^{-1} \|\mathrm{div}\, u\|_{L^2}^2$ 项虽提供控制，但直接估计会引入 $\varepsilon^{-1}$ 因子，破坏一致性；(ii) 双调和项 $\beta \Delta^2 u$ 虽增强耗散，但其高阶导数在非线性估计中产生 $\|\nabla u\|_{L^\infty}$ 类临界项，需精细处理；(iii) 斜对称非线性 $$N(u)$$ 的分解 $N(u) = \mathrm{div}(u \otimes u) - \frac{1}{2} u \, \mathrm{div}\, u$ （引理 2.3）揭示了其傅里叶表现： $\widehat{\mathrm{div}(u \otimes u)}(\xi) \sim i \xi \cdot \widehat{u \otimes u}(\xi)$ 含 $|\xi|$ 因子，而 $\widehat{u \, \mathrm{div}\, u}(\xi)$ 则需利用罚耗散控制。

作者的真正创新点在于构建一套完全 $\varepsilon$ -一致的分析框架：

弱解存在性（Theorem 1.3）：摒弃标准 Leray 投影（因其会湮灭罚项），采用 Friedrichs 型 Fourier 截断近似 $$J_n$$ （Step 1）。关键技巧在于：利用 $$J_n$$ 的正交性及 Bernstein 不等式，证明截断非线性 $$J_n N(u_n)$$ 在有限维空间 $$X_n$$ 上 Lipschitz 连续，从而获得局部解；继而通过能量恒等式 (1.5) 直接导出 $\varepsilon$ -一致先验估计（Step 2），无需对罚项做任何上界控制；最后应用 Aubin–Lions–Simon 引理（Step 3），其紧性依赖于 $\partial_t u_n \in L^{4/3}(0,T; H^{-2})$ 的一致有界性——这由方程本身及 $\varepsilon$ -一致的 $$L^2_t H^2_x$$ 估计保证，而非依赖 $\varepsilon^{-1}$ 。
强解存在唯一性（Theorem 1.4）：采用半群方法与能量法结合。首先，线性算子 $L_\varepsilon = \nu \Delta - \beta \Delta^2 + \varepsilon^{-1} \nabla \mathrm{div}$ 的傅里叶乘子 $M_\varepsilon(\xi) = (\nu |\xi|^2 + \beta |\xi|^4) I + \varepsilon^{-1} \xi \otimes \xi$ 满足统一的下界 $\langle M_\varepsilon(\xi) z, z \rangle \ge \nu |\xi|^2 |z|^2$ （引理 2.4），故其生成的半群 $S_\varepsilon(t)$ 满足 $\varepsilon$ -一致的热型估计 $\|\nabla^m S_\varepsilon(t) f\|_{L^2} \le C_m t^{-m/2} \|f\|_{L^2}$ （公式 2.6）。其次，在 $$H^2$$ 能量估计中，将非线性项 $\int \nabla^2 N(u) \cdot \nabla^2 u$ 分解并利用 $N(u) = \mathrm{div}(u \otimes u) - \frac{1}{2} u \, \mathrm{div}\, u$ ，其中 $\mathrm{div}(u \otimes u)$ 项通过 Sobolev 嵌入 $H^2 \hookrightarrow L^\infty$ 控制为 $\|u\|_{H^2}^2 \|\nabla u\|_{L^2}$ ，而 $u \, \mathrm{div}\, u$ 项则被罚耗散 $\varepsilon^{-1} \|\mathrm{div}\, u\|_{H^2}^2$ 吸收。最终，小初值条件 $\|u_0\|_{H^2} \le \delta_0$ 使非线性贡献被 $\nu \|\nabla u\|_{H^2}^2 + \beta \|\Delta u\|_{H^2}^2$ 主导，且 $\delta_0$ 的选取仅需 $\nu, \beta$ ，与 $\varepsilon$ 无关（Remark 1.5）。
最优衰减（Theorem 1.6）：采用 Schonbek 的 Fourier 分裂法，但需克服非散度自由的挑战。核心是将频率空间分裂为低频区 $|\xi| \le g(t)$ 与高频区 $|\xi| > g(t)$ ，其中 $g(t) = (1+t)^{-1/2}$ 。低频部分 $\|P_{\mathrm{low}} u(t)\|_{L^2}$ 通过 Duhamel 公式及 $$N(u)$$ 的分解估计： $\widehat{\mathrm{div}(u \otimes u)}$ 的 $|\xi|$ 因子使其在低频可被 $\|u\|_{L^2}^2$ 控制，而 $\widehat{u \, \mathrm{div}\, u}$ 则利用 $\varepsilon^{-1} \|\mathrm{div}\, u\|_{L^2}^2$ 的耗散与 $\|u\|_{L^2}$ 的衰减耦合。高频部分通过 $$H^2$$ 能量估计与 Gagliardo–Nirenberg 插值 $\|\nabla u\|_{L^\infty} \lesssim \|u\|_{H^2}^{1/2} \|\Delta u\|_{H^2}^{1/2}$ 控制。全程所有常数均追溯至 $\nu, \beta$ 及 $$L^1$$ 初值范数，彻底规避 $\varepsilon$ 依赖。

4. 与之前工作的比较

本文与既有工作的实质性推进体现在四个维度：

与经典罚方法（Temam [23], Shen [19]）比较：传统罚法仅含 $-\nu \Delta u - \varepsilon^{-1} \nabla \mathrm{div}\, u$ ，其弱解理论要求 $\varepsilon$ 固定，且 $\varepsilon \to 0$ 极限的收敛性证明常需额外假设（如 $\varepsilon$ -一致 $$H^1$$ 估计）。本文首次在罚框架中引入双调和项 $\beta \Delta^2 u$ ，其提供的 $$H^2$$ 正则性不仅强化了先验估计，更关键的是使 $\varepsilon$ -一致 $$H^2$$ 能量估计（Theorem 1.4）与最优衰减（Theorem 1.6）成为可能——经典罚法因缺乏高阶耗散，无法获得 $\|\Delta u\|_{L^2_{t,x}}$ 一致控制，故无法建立 $\varepsilon$ -一致强解理论。
与超粘性 NS（Lions [12], Tao [22]）比较：Lions 证明 $\partial_t u + (-\Delta)^{5/4} u + (u \cdot \nabla) u + \nabla p = 0$ 具有任意初值的整体正则性，但其模型仍含压力 $$p$$ 与不可压约束 $\mathrm{div}\, u = 0$ ，分析中压力项带来技术复杂性。本文模型 $\partial_t u + L_\varepsilon u + N(u) = 0$ 彻底消除了压力，且通过 Temam 校正维持 $$L^2$$ 能量封闭。更重要的是，Lions 型结果不涉及罚参数，故其估计天然不具 $\varepsilon$ -一致性；而本文所有结果均显式构造为 $\varepsilon$ -一致，为 $\varepsilon \to 0$ 极限分析铺平道路。
与阻尼型 NS（Cai–Jiu [4], Zhou [30], Xu–Zhou [27]）比较：阻尼项 $\alpha |u|^{r-1}u$ 属零阶非线性修正，其耗散机制与频域无关；而双调和项 $\beta \Delta^2 u$ 是线性、四阶、频域选择性耗散（主要作用于高频）。本文证明：尽管机制迥异，但双调和项与罚项的组合能产生与最优阻尼（ $$r=3$$ ）相当的 $\varepsilon$ -一致小初值整体正则性（Theorem 1.4），且进一步给出最优大时间衰减（Theorem 1.6）——而多数阻尼型工作（如 [4, 27]）未系统研究衰减率或其 $\varepsilon$ -一致性。
与 Fourier 分裂法（Schonbek [17], Wiegner [26]）比较：Schonbek 对经典 NS 的衰减分析依赖于 $\mathrm{div}\, u = 0$ 及压力消失，其低频估计基于 $\widehat{u \cdot \nabla u}(\xi) \sim i \xi \cdot \widehat{u \otimes u}(\xi)$ 的 $|\xi|$ 因子。本文必须处理 $\mathrm{div}\, u \neq 0$ ，故引入 $$N(u)$$ 的关键分解 $N(u) = \mathrm{div}(u \otimes u) - \frac{1}{2} u \, \mathrm{div}\, u$ （引理 2.3）：前者保留 $|\xi|$ 因子，后者则通过罚耗散 $\varepsilon^{-1} \|\mathrm{div}\, u\|_{L^2}^2$ 控制。此分解是本文独有的技术，确保 Fourier 分裂法在非散度自由框架下仍能导出与热方程相同的最优速率（Remark 1.7）。此外，所有常数均被显式追踪为 $\varepsilon$ -独立，而 Schonbek 原始工作未考虑参数一致性问题。

中文标题: 具有反平方渐近性的径向薛定谔算子分数幂的尖锐分段幂次Lorentz估计

摘要翻译: 设 $$H=-Δ+V(|x|)$$ 是定义在 $\mathbb{R}^d$ （ $d\ge 2$ ）上的非负径向薛定谔算子，其正调和函数满足 $U(r)\simeq r^{-σ_0}$ （当 $0<r\le 1$ ）与 $U(r)\simeq r^{-σ_\infty}$ （当 $r\ge 1$ ），其中 $-d/2<σ_0,σ_\infty<d/2$ 。在假设 Ishige、Kabeya 与 Ouhabaz 所建立的双向基态热核估计成立的前提下，我们确定了 $H^{-s/2}$ 的核 $$K_s^H(x,y)$$ 具有干净双向估计

K_s^H(x,y)\simeq |x-y|^{s-d}U(|x|)U(|y|)/[U(|x|+|x-y|)U(|y|+|x-y|)]

的最大开区间，即 $0<s<\min\{d,d-2σ_0,d-2σ_\infty\}$ 。在此区间内，我们对分段幂次估计

\|w_{-β_0,-β_\infty}H^{-s/2}f\|_{L^{q,v}}\lesssim \|w_{α_0,α_\infty}f\|_{L^{p,u}}

给出了关于 $1<p,q<\infty$ 与 $1\le u,v\le\infty$ 的完备充要分类。该结果涵盖带符号的基态指数、整个 $$q<p$$ 范围、所有单侧权等式、两种尺度等式以及同时发生的端点角点情形。尺度等式由 $u\le v$ 控制（包括 $$q<p$$ 情形）；输入或输出幂次端点分别要求 $$u=1$$ 或 $v=\infty$ ；而在同侧幂次/尺度角点处，唯一容许的对是 $(u,v)=(1,\infty)$ 。证明综合运用了干净核分析、局部 Lorentz–Hardy–Littlewood–Sobolev 估计、秩一端点论证、几何环形序列空间、三角矩阵定理以及九块分解法。

9. 具有Hardy奇点的临界Schrödinger–Poisson系统的归一化态：多重性与半经典集中

Normalized states for a critical Schrödinger–Poisson system with Hardy singularity: multiplicity and semiclassical concentration

作者: Khaled Khachnaoui

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11277v1

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

摘要翻译

中文标题: 具有Hardy奇点的临界Schrödinger–Poisson系统的归一化态：多重性与半经典集中

摘要翻译: 本文研究了在 $\mathbb{R}^3$ 中带有奇异Hardy势的临界Schrödinger–Poisson系统之归一化解的存在性、多重性与半经典集中现象。更精确地，我们考虑如下系统：

\begin{cases} -\varepsilon^2\Delta u+ \left(V(x)-\dfrac{κ\varepsilon^2}{|x|^2}\right)u -φ|u|^3u =λu+μ|u|^{q-2}u+|u|^4u, & \text{in } \mathbb{R}^3, \\[1mm] -\varepsilon^2\Delta φ=|u|^5, & \text{in } \mathbb{R}^3, \end{cases}

并施加给定的质量约束

中文标题: 通过多尺度磨光度量实现的奇异平均场极限

中文标题: Fabes-Kenig-Serapioni 奇异-退化型抛物方程的加权 $W^{1,p}$ -估计

摘要翻译: 我们研究一类散度型二阶抛物方程的 Dirichlet 边值问题，其系数矩阵具有由 Muckenhoupt 权类刻画的奇异与退化行为。该框架是 Fabes、Kenig 与 Serapioni 所开创的奇异-退化椭圆方程的抛物类比。在系数的部分加权平均振荡满足小性假设的前提下，我们在适当定义的加权 Sobolev 空间中建立了弱解的存在性、唯一性以及局部内点与边界正则性估计。证明依赖于系数冻结技巧以及 Caffarelli 与 Peral 引入的水平集方法。此外，我们构建了所需的加权 Sobolev 空间框架及相关加权不等式。最后，利用紧性论证，我们证明：这些方程的解在加权 Sobolev 范数意义下，局部地与其冻结系数对应解保持接近。

不相关论文（基本信息）

1. 面向高维偏微分方程的Inf-Sup神经网络

Inf-Sup Neural Networks for High Dimensional PDEs

作者: Ziren Chen, Hailiang Liu

类别: math.NA, math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11718v1

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

摘要: Solving partial differential equations (PDEs) in high dimensions remains challenging due to the curse of dimensionality. We propose a neural-network-based framework that reformulates PDEs as inf–sup optimization problems through the introduction of a Lagrange multiplier. The primal solution and the associated Lagrange multiplier are parameterized by two networks and are computed via an iterative saddle-point optimization procedure. We prove the theoretical equivalence between the proposed optimization formulation and the original PDE problem, and we derive rigorous error estimates that quantify the total approximation error in terms of the network approximation error, statistical (sampling) error, and optimization error. Numerical experiments demonstrate the accuracy, stability, and efficiency of the proposed method for solving high-dimensional PDEs.

2. Dirichlet 扰动下的 Serrin 问题：几何紧性与尖锐平面稳定性

Serrin’s Problem under Dirichlet Perturbations: Geometric Compactness and Sharp Planar Stability

作者: Qinfeng Li, Weihong Xie, Hang Yang

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11812v1

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

摘要: In earlier work, we posed a stability question for Serrin’s overdetermined problem under Dirichlet perturbations and proved that the answer is negative in dimensions $n\ge3$ . Here we resolve the question in the planar convex class and obtain a sharp quantitative theory without any a priori geometric nondegeneracy. Let $$u_Ω$$ solve

-Δu_Ω=1\ \text{in }Ω,\qquad \partial_νu_Ω=-\frac{|Ω|}{P(Ω)}\ \text{on }\partialΩ, \qquad \int_{\partialΩ}u_Ω\,dσ=0,

and set $O(Ω):=\text{osc}_{\partial Ω}u_Ω$ . We construct fixed-area annuli with $O(Ω_k)\to0$ that remain far from every disk, showing that convexity is essential in dimension two. By contrast, if $Ω_k\subset\mathbb R^2$ are convex, $$|Ω_k|=π$$ , and $O(Ω_k)\to0$ , then, up to translations, $$Ω_k$$ converges in Hausdorff distance to the unit disk. Moreover,

R_Ω-r_Ω+\inf_{z\in\mathbb R^2}d_H(Ω,B_1(z)) \le C\,O(Ω)

for all planar convex $$Ω$$ with $$|Ω|=π$$ and sufficiently small $$O(Ω)$$ , and the linear order is optimal. The proof combines a new mechanism excluding long-thin degeneration, the rough-domain Serrin rigidity theorem of Figalli–Zhang, new tangential-gradient and linear boundary-growth estimates, a boundary $$P$$ -function estimate, and the reverse-Serrin identity of Magnanini–Molinarolo–Poggesi. We also study the weaker deficit

A(Ω):=\frac1{P(Ω)}\int_{\partialΩ}u_Ω,dσ-\min_{\partialΩ}u_Ω.

In the planar convex class, $A(Ω_k)\to0$ still forces convergence to a disk, and

R_Ω-r_Ω+\inf_z d_H(Ω,B_1(z)) \le C A(Ω)^{2/3}

for $$|Ω|=π$$ and sufficiently small $$A(Ω)$$ .

3. 面向噪声型积分-发放神经元群体平均场极限的随机 Fokker–Planck 方程

A stochastic Fokker–Planck equation for the mean-field limit of a population of noisy integrate-and-fire neurons

作者: Ben Hambly, Aldaïr Petronilia, Christoph Reisinger, Andreas Søjmark

类别: math.AP, math.PR

发布日期: 2026-07-13

arXiv ID: 2607.11685v1

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

摘要: We study a densely connected system of excitatory integrate-and-fire neurons which are subject to common noise. The system incorporates a gradual transmission of action potentials and captures the re- and hyperpolarization phases through a random refractory period followed by a reset to a randomized level below the rest potential. As the number of neurons tends to infinity, we show that there is weak convergence to a unique membrane potential density governed by a stochastic Fokker–Planck equation with a well-defined spike transmission rate. The latter is driven by the mean cumulative spike count, which is non-differentiable but shown to satisfy a generalized flux condition. We obtain the uniqueness of the Fokker–Planck equation from energy estimates in the dual of the first Sobolev space. Finally, we give a conditional McKean–Vlasov representation of the membrane potential density as the law of a representative neuron given the common noise.

4. 具有次Lipschitz Hessian的Bellman方程

Bellman Equations with Sub-Lipschitz Hessians

作者: Xavier Fernández-Real

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11618v1

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

摘要: We construct homogeneous solutions with non-Lipschitz Hessian for finite, constant-coefficient Bellman equations. First, for every $σ\in(0,1)$ , we find two uniformly elliptic matrices $A_1,A_2\in\mathcal{S}^4$ and a nonzero $$(2+σ)$$ -homogeneous solution $$u$$ of

\max\bigl\{{\rm tr}\,(A_1D^2u),{\rm tr}\,(A_2D^2u)\bigr\}=0 \qquad\text{in }\mathbb{R}^4.

Second, in $\mathbb{R}^2$ we construct three matrices satisfying ${\rm Id}_2\leq A_j\leq3{\rm Id}_2$ for which the corresponding Bellman equation admits a homogeneous solution with a non-Lipschitz Hessian. In particular, solutions to convex fully nonlinear uniformly elliptic equations are not in $C^{2,1}$ , and not even in $C^{2, 1-\varepsilon}$ for $\varepsilon > 0$ small.

5. 矩结构分块编码的周期性有限差分算子

Moment-Structured Block Encodings of Periodic Finite-Difference Operators

作者: Jishnu Mahmud, Rebekah Herrman

类别: quant-ph, math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11596v1

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

摘要: Block encoding is the standard technique for accessing matrix data in quantum linear-algebra algorithms. Its implementation directly affects its subnormalization, which in turn controls the algorithm’s success probability, simulation time, and downstream costs. Explicit construction of block encodings with provably optimal subnormalization exists for only a handful of operators, with bespoke calculations used in its design. In this work, we develop a framework to block encode translation-invariant finite-difference operators on a periodic grid. These operators are the finite-difference discretizations of the constant-coefficient partial differential equations that sit at the core of scientific computing. We show that the moment order of these stencils can be used to simultaneously determine the continuum operator approximated, the vanishing order of the Fourier symbol, and the cost of the block encoding. From there, we derive a closed-form optimality criterion as a function of the stencil coefficients, which certifies whether the construction attains the optimal subnormalization for an entire operator family, uniformly in grid size, and quantifies the gap when it does not. The framework subsumes optimal constructions for the Laplacian operator in the literature and can be used to certify new instances at higher even orders, including the biharmonic operator. Furthermore, we derive success-probability floors parameterized by spectral properties of the operator’s symbol and find explicit constants for the block encoding of the advection-diffusion family for which no prior explicit spatial block encoding exists.

6. 带势的局部-非局部算子首加权特征值的优化

Optimization for the first weighted eigenvalue of local-nonlocal operators with a potential

作者: R. Lakshmi, Sekhar Ghosh, Ratan Kr. Giri

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11527v1

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

摘要: In this work, we study the optimization problem for the first eigenvalue of mixed local nonlocal operators plus a potential $$V$$ . We begin by investigating the existence and fundamental properties of the eigenvalues, with special emphasis on the first eigenvalue. Finally, we discuss the dependence of the first eigenvalue on the potential function and establish the existence of optimal potentials within certain admissible classes. In particular, we show the existence of a unique maximizer and a minimizer of the first eigenvalue on any bounded, closed, and convex subset of $$L^q(Ω)$$ . Moreover, these results enable us to characterize the maximizers and minimizers in the closed unit ball of $$L^q(Ω)$$ , as well as in the class of rearrangements of any $V\in L^q(Ω)$ .

7. 弱耦合椭圆系统中的图相容幂凹性

Graph-Compatible Power Concavity in Weakly Coupled Elliptic Systems

作者: Jiahuan Li, Hanpeng Zhou

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11513v1

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

摘要: We introduce graph-compatible powers determined by the coupling graph for multicomponent weakly coupled elliptic systems and establish concavity of the associated transformed components through a concave-envelope method and a weighted viscosity comparison principle. We further establish a componentwise constant-rank theorem for the transformed Hessians, which yields strict power concavity under additional assumptions.

8. Lieb-Oxford 界的非紧支集极小化子的不存在性

Absence of non-compactly supported minimisers for the Lieb-Oxford bound

作者: Simone Di Marino, Rodrigue Lelotte

类别: math.AP, math-ph

发布日期: 2026-07-13

arXiv ID: 2607.11440v1

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

摘要: We prove the minimisers of the Lieb-Oxford bound (if any) for a fixed (and finite) number of particles $N \geq 1$ are necessarily compactly supported, extending a result proved by E. H. Lieb and S. Oxford in the one-particle case.

9. 一种考虑三相的、具有Herschel-Bulkley流变律的瞬态深度平均熔岩流模型

A transient depth-averaged lava flow model with a Herschel-Bulkley rheology accounting for three phases

作者: Julie Binard, Alain Burgisser, Enrique D. Fernández-Nieto, Gladys Narbona-Reina

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11383v1

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

摘要: This study presents a three-phase suspension lava flow model with a Herschel-Bulkley rheology. The suspension contains crystals and gas bubbles, and two closures for the evolution of the crystal volume fraction are considered. The first closure minimizes the complexity of the system by treating crystal fraction as a transported quantity subject to relaxation towards an equilibrium state. This closure avoids a parametrization of the many heat transfer mechanisms a lava flow is subjected to. The other closure is on the lava temperature considering four heat transfer mechanisms and a prescribed temperature–crystallinity relationship. We deduce from this system and solve numerically a one-dimensional depth-averaged model. A comparison with a pre-existing model based on real lava flow data suggests that the prediction of flow parameters done with the multi-parametric evolution of temperature yields more accurate results than those obtained with the other closure. The transient nature of our model correctly predicts that confined lava traveling down an irregular steep slope yields a series of cascading fill-then-breakout lumps that causes the overall flow to be pulsatory. These pulses dominate the local dynamics and preclude a strict steady state to be reached. Theoretically, taking gas bubbles into account is best done with a general rheological relationship valid at any capillary number. In the conditions explored herein, bubbles modulate viscosity within a factor 2 with a shear thinning behavior, decelerating slow flows and accelerating fast flows. When a simplified rheology treating bubbles as hard spheres was used, the only dynamic parameter affected was bulk viscosity.

10. 年龄结构化资源模型中的竞争排斥与共存

Competitive exclusion and coexistence in a model with age-structured resource

作者: Simon Girel, Quentin Richard

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11247v1

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

摘要: In this paper, we introduce and analyse a model in which a shared, age-structured resource is exploited by two competing populations through age-dependent interaction kernels. It provides a general framework for studying resource-mediated competition in epidemiology, immunology and ecology. We use integrated semigroups theory to establish the well-posedness of the model and its fundamental properties, including positivity and the existence of a global attractor. Then, by means of Lyapunov functionals, persistence results and generalised LaSalle principle, we fully characterize the asymptotic dynamics in terms of the competitors’s basic reproduction numbers and invasive fitnesses. Our results demonstrate that in addition to the classical competitive exclusion, the system can exhibit globally asymptotically stable coexistence, an outcome allowed by a mechanism of age-mediated niche partitioning.

11. 从类空初始数据形成陷俘曲面

Formation of Trapped Surfaces from Spacelike Initial Data

作者: Xuantao Chen, Sergiu Klainerman

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

发布日期: 2026-07-13

arXiv ID: 2607.11236v1

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

摘要: We make use of the free data formalism developed in \cite{CK25,CK26} to provide a direct construction of short-pulse type Cauchy data. The construction of the spacelike short-pulse follows from the local existence result established in \cite{CK26}. The forward integration construction in \cite{CK26} allows us to show that such data can be extended to a set of asymptotically flat Cauchy data. This greatly extends the result of Li–Yu \cite{LiYu}.

12. Brezis-Nirenberg 方程的精细爆破分析及其应用： $N\geq4$ 时的单气泡情形

A refined blow-up analysis of the Brezis-Nirenberg equation and its application: The one-bubble case for $N\geq4$

作者: Rui He, Xiangqing Liu, Juncheng Wei, Yuanze Wu

类别: math.AP

发布日期: 2026-07-13

arXiv ID: 2607.11132v1

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

摘要: In this paper, we consider the famous Brezis-Nirenberg equation \begin{eqnarray*} \left{ \aligned &-Δu=λu+|u|^{\frac{4}{N-2}}u,\quad&\mbox{in},, Ω,\ &u=0,\quad&\mbox{on},, \partialΩ, \endaligned \right. \end{eqnarray*} where $N\geq3$ is the dimension, $Ω\subset\mathbb{R}^N$ is a bounded domain with smooth boundary $\partialΩ$ and $$λ>0$$ is a parameter. By developing a refined blow-up analysis based on the inverse reduction argument developed in \cite{WW2019,WW2019-2}, we classify, for the fist time, the Struwe decomposition of the Brezis-Nirenberg equation in the one-bubble case as the parameter $$λ$$ varies for $N\geq4$ . As applications, we prove that the $$4d$$ Brezis-Nirenberg equation has a nontrivial solution (least energy solution) for $λ\inσ(-Δ)$ in general bounded domains, where $$σ(-Δ)$$ is the spectrum of $$-Δ$$ in $$H^1_0(Ω)$$ . Our result completes the existence theory of the Brezis-Nirenberg equation for $N\geq4$ in \cite{AP2025,CFP1985,CFS,CSS1986,CW2005,CSZ2012,SWW2009,TYZ2022} since 1984.

13. 通过重耦合影子流实现的均场控制中的最优收敛速率

The optimal rate of convergence in mean field control via recoupled shadow flows

作者: Sebastian Munoz

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

发布日期: 2026-07-13

arXiv ID: 2607.11062v1

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

摘要: We prove that the value functions of the $$N$$ -particle stochastic optimal control problem converge uniformly to the value function of the corresponding mean field control problem at the empirical-measure rate ( $N^{-1/d}$ for $d\geq3$ and $N^{-1/2}\sqrt{\log N}$ for $$d=2$$ ), for mean field costs that are merely Lipschitz continuous in the $$1$$ -Wasserstein distance. For $d\geq2$ the rate is optimal in this class, which covers problems whose mean field optimizers are neither unique nor stable. This proves the rate conjectured by Daudin, Delarue and Jackson, and removes the semiconcavity hypothesis made there. The proof is control-theoretic: from each realization of an $$N$$ -particle control we build a pathwise Fokker–Planck flow – a shadow flow – which, repeatedly recoupled to the particles by optimal transport, shadows the empirical measure at the empirical-measure rate. The same rate holds with additive common noise, uniformly in its intensity. Finally, in dimension one we discover that the empirical-measure benchmark is not optimal: cooperating particles beat it, and the optimal polynomial exponent is $$4/7$$ , strictly between the accuracy of independent samples and that of quantization by freely placed points. The proofs combine an anticipating correction of the shadow flow, a Gibbs law implementing the cooperation, and a Schrödinger ground-state estimate for the optimality.

14. 透镜空间上热流的最终 Morse 极小性

Eventual Morse Minimality of Heat Flow on Lens Spaces

作者: Carlos A. Cadavid, Dairo J. Hernández, Juan Diego Vélez

类别: math.DG, math.AP

发布日期: 2026-07-12

arXiv ID: 2607.10837v1

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

摘要: We prove that heat flow generically produces Morse-minimal functions on round lens spaces. The result gives a rigorous explanation of a phenomenon first suggested by computational experiments: after high-frequency modes have been suppressed by diffusion, generic heat evolutions on lens spaces tend to settle into Morse functions with the smallest possible number of critical points. Precisely, for every lens space $$L(p,q)$$ with $p\geq2$ , $1\leq q\leq p/2$ , and $$(p,q)=1$$ , there is an open dense set of initial data $f\in L^2(L(p,q),\mathbb R)$ such that the heat evolution $e^{tΔ}f$ , where $$Δ$$ is taken in the non-positive sign convention, is, for all sufficiently large $$t$$ , a Morse function with exactly four critical points, of indices $$0,1,2,3$$ . The proof analyzes the asymptotic spectral expansion of the heat flow. In the most delicate case $$1<q<p/2$$ , the leading term is Morse–Bott with two critical circles, and the higher heat modes break these circles through a resonant Fourier mechanism. An arithmetic estimate shows that the fundamental reduced frequency dominates all higher multiples, and generic nonvanishing of the corresponding resonant coefficients gives the minimal critical-point count.

15. Cahn-Hilliard 时间步的四场辅助重构：分析与协调有限元离散化

A four-field auxiliary reformulation of a Cahn-Hilliard time step: Analysis and conforming finite element discretization

作者: Marvin Fritz

类别: math.NA, math.AP

发布日期: 2026-07-12

arXiv ID: 2607.10724v1

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

摘要: We propose a four-field auxiliary reformulation of a time-discrete Cahn-Hilliard step. The construction is motivated by the scalar trace structure underlying two-dimensional Rafetseder-Zulehner decompositions and represents the scalar quantity $$-Δc$$ in the form $2p+div\,u$ through an auxiliary scalar field $$p$$ and an auxiliary vector field $$u$$ . The resulting auxiliary problem is a mixed second-order Stokes/elasticity-type system, while the evolution equation for the phase field retains its standard mass-conserving gradient-flow structure. We derive a continuous four-field formulation that is equivalent to the classical convex-splitting mixed time step. We also state a conforming finite element discretization and prove one-step spatial estimates for the phase-field variables together with a stable auxiliary block estimate containing an explicit weak-Laplacian recovery defect. The numerical experiments verify the expected phase-field convergence in a manufactured setting, mass conservation, energy decay, and projected consistency of the auxiliary block.

16. 通过分支理论研究 McKean-Vlasov 相变的连续性与不连续性

Continuity and Discontinuity of McKean-Vlasov Phase Transitions via Bifurcation Theory

作者: Junlang Hu, Zhenxin Liu

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

发布日期: 2026-07-12

arXiv ID: 2607.10723v1

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

摘要: It is well known that the McKean-Vlasov stochastic differential equation with a symmetric double-well potential exhibits a continuous phase transition. In contrast, for an asymmetric double-well potential, the system undergoes a discontinuous phase transition, in which an invariant measure abruptly appears in the shallower well and subsequently splits into two as the temperature decreases. In this work, we systematically investigate the types of phase transitions driven by non-convex confining potentials and cooperative interactions in $\mathbb{R}^n$ . Our main results consist of a pitchfork bifurcation theorem characterizing continuous phase transitions and a saddle-node bifurcation theorem governing discontinuous phase transitions. In addition, despite the dimension-dependent nature of phase transitions, we are able to extend Dawson’s criterion for phase transition points – originally formulated for quadratic interactions in $$1$$ -dimensional space – to $$n$$ -dimensional space. This extension directly relates the critical parameter and bifurcation direction to the eigenvalue and eigenvector of the critical covariance matrix, respectively. Finally, we apply our theoretical results to the aforementioned double-well models, to a class of four-well models in two-dimensional Euclidean space that exhibit multiple phase transitions, and to an attractive Gaussian interaction model.

17. 描述非线性连续体振动的偏微分方程的谱子流形约化

Spectral submanifold reduction for PDEs describing nonlinear continuum vibrations

作者: Gergely Buza, George Haller

类别: math.DS, math.AP

发布日期: 2026-07-12

arXiv ID: 2607.10675v1

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

摘要: We prove the existence of spectral submanifolds in nonlinear partial differential equations (PDEs) describing forced-damped continuum vibrations. Our results cover structural vibrations subject to generalized structural damping. Based on these results, backbone curves and forced response curves can be rigorously extracted from PDEs, without a priori discretization. On two examples involving an elastic beam and a thin plate, we show that backbone and forced response curves can even be calculated by hand directly from the PDEs.

18. 平面的各向异性等周双铺砌

Anisotropic isoperimetric double tilings of the plane

作者: Francesco Nobili, Matteo Novaga, Emanuele Paolini

类别: math.AP, math.MG

发布日期: 2026-07-12

arXiv ID: 2607.10636v1

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

摘要: We study periodic partitions of the plane into two distinct cells minimizing the anisotropic $\ell_1$ -perimeter. For a rectangular lattice $$G$$ , we compute explicitly the $(G,\ell_1)$ -isoperimetric profile and classify all minimizers. When one cell has small area, the optimal tiling is generated by a square and a chipped rectangle, while in the remaining regime the tiling is generated by two adjacent rectangles. Further minimizing over all possible planar lattices, we show that the $\ell_1$ -isoperimetric profile is attained by the Pythagorean double tiling of two axis-aligned squares sharing a vertex. This configuration is unique unless the two cells are assigned the same area. Finally, we prove that the limiting (non periodic) partitions obtained by sending one volume to infinity are locally $\ell_1$ -isoperimetric.

19. 非各向同性分数阶拉普拉斯方程、亥姆霍兹方程与尤卡瓦方程的Walk-on-Cubes蒙特卡洛模拟

Walk-on-Cubes Monte Carlo Simulation for nonisotropic fractional Laplace, Helmholtz, and Yukawa equations

作者: Antti Rasila, Tommi Sottinen, Yaotong Yuan

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

发布日期: 2026-07-12

arXiv ID: 2607.10528v1

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

摘要: We study the nonisotropic fractional analogs of Laplace, Helmholtz and Yukawa equations. We provide a Duffin correspondence for the Yukawa equation and a Feynman–Kac reconstruction for the Helmholtz equation. The foundation of our analysis is the fact that the nonisotropic fractional Laplace equation is related to a symmetric $$α$$ -stable Lévy process with independent identically distributed components. By using this relation we provide a Walk-on-Cubes algorithm that simulates the solutions of Helmholtz and Yukawa equations.

20. 四维Anderson模型：临界随机偏微分方程的一个案例研究

The four-dimensional Anderson model: a case study for critical SPDEs

作者: Yu Deng, Hao Shen

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

发布日期: 2026-07-11

arXiv ID: 2607.10105v1

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

摘要: We study the weakly coupled elliptic Anderson model with spatial white noise on the four-dimensional torus, which provides a basic example of a critical SPDE requiring renormalization at arbitrarily high orders. With coupling $λ|\log\varepsilon|^{-\frac12}$ where $$λ>0$$ is sufficiently small, we prove that the Green’s function of the corresponding random Schrödinger operator, suitably centered and rescaled, converges to a centered Gaussian random field with explicit covariance. The main difficulty is that, for such critical models, one must expand up to order $|\log\varepsilon|$ , while the perturbative expansion contains factorially many pairings and a growing number of renormalization terms. To overcome this, we construct a truncated renormalized parametrix and prove sharp high-order bounds for its remainder. A central ingredient is a multiscale analysis based on a new version of Hepp trees, combined with new estimates for summations over permutations. These estimates reveal a precise balance between logarithmic losses from scale summation and factorial gains from the structure of primitive pairings. The methods developed here are intended as a first step toward a general theory for critical SPDEs with weak couplings.

21. 等周最优控制中的乘子敏感性

Multiplier Sensitivity in Isoperimetric Optimal Control

作者: Ye Liang

类别: math.AP

发布日期: 2026-07-10

arXiv ID: 2607.09944v1

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

摘要: We study finite-horizon optimal control problems with scalar isoperimetric constraints from a function-space duality perspective. Controls are treated as elements of weakly compact subsets of $L^\infty(0,T;\mathbb R^m)$ , while the state equation induces a control-to-state map into $W^{1,\infty}(0,T;\mathbb R^n)$ . For linear dynamics, concave payoff, and an affine isoperimetric functional, we prove that the constrained value function has an interval domain, is concave, and admits a Fenchel–Moreau dual representation. This yields a superdifferential formula identifying the negative of the dual multiplier with the sensitivity of the value function with respect to the constraint level. A constraint qualification is then used to identify the dual multiplier with the normal Pontryagin multiplier of the augmented isoperimetric system. We also treat linear-quadratic problems with a single quadratic equality constraint by reducing them to quadratic forms on a Hilbert space. The resulting analysis separates the validity of the envelope formula from the regularity needed for Riccati synthesis, showing that sensitivity may persist even when the modified control-weight operator becomes singular.

22. 具有最大值型收益的反射最优停时：测度值停时收益与被杀死的预解表示

Reflected Optimal Stopping with a Max-Type Payoff: Measure-Valued Stopping Gains and Killed Resolvent Representation

作者: Ye Liang

类别: math.AP, math.PR

发布日期: 2026-07-10

arXiv ID: 2607.09987v1

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

摘要: We study an infinite-horizon optimal stopping problem for a two-dimensional normally reflected diffusion in the quadrant with payoff $G(x_1,x_2)=x_1\vee αx_2$ . The problem combines three features that complicate the usual free-boundary analysis: reflection on the coordinate axes, a genuinely two-dimensional stopping region, and a nonsmooth max-type reward. We formulate the associated reflected obstacle problem, prove a verification theorem under explicit Itô–Krylov–Tanaka admissibility and measure-superharmonicity assumptions, and derive a conditional epigraph structure for the stopping set. The main technical point is that the stopping-gain object $Γ=c+rG-\mathcal LG$ is a signed measure rather than a function. Its diagonal component is $Γ^Δ(dx) = -\frac{n^\top a(x)n}{2\sqrt{1+α^2}}σ_Δ(dx)$ , $$n=(1,-α)$$ , which shows that pointwise stopping-gain sign conditions must be interpreted with care. We also prove that the correct potential representation is the killed-resolvent formula $V(x)=G(x)-R_r^{\mathcal C}Γ(x)$ , rather than the unrestricted reflected resolvent. A constant-coefficient reflected Brownian example illustrates the diagonal singular term explicitly.

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

1. 关于三维Peskin问题在临界正则性下的全局渐近稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 空间异质环境下混合反应-扩散传染病模型中的波前传播. 第 II 部分：脉动行波#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 非齐次抛物型方程的有限时间爆破#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 非均匀不可压缩热传导流体在穿孔区域上的均匀化#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. $L^2(\mathbb{R}^2)$ 适定性与密度块问题的对数Lipschitz正则性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. 演化边界时滞方程#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 一个与平衡律耦合的输运方程的存在性与唯一性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. 具间断通量的粘性守恒律的存在性、唯一性与大时间行为#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

9. 带加性噪声的三维 Navier–Stokes 方程的误差分析#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

10. 三维带双调和阻尼的罚 Navier–Stokes 系统的整体存在性与最优衰减#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 大势能下波动方程的一致可控性#

摘要翻译#

2. 弱相互作用聚焦Gibbs测度在谐势下的相变#

摘要翻译#

3. 关于具有缓慢衰减数据的非局域非线性薛定谔方程的逆散射问题#

摘要翻译#

4. 从扩散到传输 via EDP-收敛：一种典型的多尺度极限#

摘要翻译#

5. 粗糙区域上散度测度张量场的 Gauss-Green 公式#

摘要翻译#

6. 分数阶对数 $p$-拉普拉斯算子的正则性#

摘要翻译#

7. Nakao问题在二维和三维空间中的小初值全局时间存在性#

摘要翻译#

8. 具有反平方渐近性的径向薛定谔算子分数幂的尖锐分段幂次Lorentz估计#

摘要翻译#

9. 具有Hardy奇点的临界Schrödinger–Poisson系统的归一化态：多重性与半经典集中#

摘要翻译#

10. 抽象波动方程的弱可观测性刻画#

摘要翻译#

11. 等离子体-壁面动力学：适定性#

摘要翻译#

12. 关于 Calogero-Moser 导数非线性 Schrödinger 方程的直接散射理论#

摘要翻译#

13. 有界域上带 Neumann 边界条件的一类梯度型 Schrödinger 系统的归一化解#

摘要翻译#

14. 基于激波追踪格式的 $p$-系统弱解的存在性与稳定性估计#

摘要翻译#

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

1. 关于三维Peskin问题在临界正则性下的全局渐近稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 空间异质环境下混合反应-扩散传染病模型中的波前传播. 第 II 部分：脉动行波

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 非齐次抛物型方程的有限时间爆破

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 非均匀不可压缩热传导流体在穿孔区域上的均匀化

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. $L^2(\mathbb{R}^2)$ 适定性与密度块问题的对数Lipschitz正则性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 演化边界时滞方程

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 一个与平衡律耦合的输运方程的存在性与唯一性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 具间断通量的粘性守恒律的存在性、唯一性与大时间行为

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

9. 带加性噪声的三维 Navier–Stokes 方程的误差分析

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

10. 三维带双调和阻尼的罚 Navier–Stokes 系统的整体存在性与最优衰减

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 大势能下波动方程的一致可控性

摘要翻译

2. 弱相互作用聚焦Gibbs测度在谐势下的相变

摘要翻译

3. 关于具有缓慢衰减数据的非局域非线性薛定谔方程的逆散射问题

摘要翻译

4. 从扩散到传输 via EDP-收敛：一种典型的多尺度极限

摘要翻译

5. 粗糙区域上散度测度张量场的 Gauss-Green 公式

摘要翻译

6. 分数阶对数 $$p$$ -拉普拉斯算子的正则性

摘要翻译

7. Nakao问题在二维和三维空间中的小初值全局时间存在性

摘要翻译

8. 具有反平方渐近性的径向薛定谔算子分数幂的尖锐分段幂次Lorentz估计

摘要翻译

9. 具有Hardy奇点的临界Schrödinger–Poisson系统的归一化态：多重性与半经典集中

摘要翻译

10. 抽象波动方程的弱可观测性刻画

摘要翻译

11. 等离子体-壁面动力学：适定性

摘要翻译

12. 关于 Calogero-Moser 导数非线性 Schrödinger 方程的直接散射理论

摘要翻译

13. 有界域上带 Neumann 边界条件的一类梯度型 Schrödinger 系统的归一化解

摘要翻译

14. 基于激波追踪格式的 $$p$$ -系统弱解的存在性与稳定性估计

摘要翻译