生成时间: 2026年05月25日 07:08:36

运行状态: 完成

计划处理论文数量: 22

已完成论文数量: 22

未完成论文数量: 0

重点关注论文数量: 5

了解领域论文数量: 11

不相关论文数量: 6

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

1. 利用 Lambert W 函数求解 Navier-Stokes 方程的一个简单情形

Solution of a Simple Case of the Navier-Stokes Equations via Employing the Lambert W Function

作者: József Vass

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23305v1

摘要翻译: 本文旨在引入一个由两个变量构成的奇异函数，该函数可通过 Lambert W 函数 $$W$$ 表达；此函数可被推广，以在特定区域上满足无粘运动的 Euler 方程 $\partial_t u + (u \cdot \nabla) u = -\nabla p$ , 且压力 $$p$$ 与空间变量无关。

Comment: This paper was presented at the Scientific Student Conference of Eötvös Loránd University on Nov. 26, 2003 with the same title. (Contains 6 pages with 0 figures.)

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

1. 研究对象和背景

本文研究的是不可压Navier-Stokes方程在零粘性极限下的一个特殊可解情形——即Euler方程的显式解析解构造问题。具体而言，作者考虑如下 $$ n $$ 维欧拉方程（无粘、无外力、压力仅依赖于时间）：

\partial_t u + (u \cdot \nabla) u + \nabla p = 0, \quad \text{div}\, u = 0,

但明确指出其构造的解 $u = (u_1,\dots,u_n)$ 不满足不可压条件 $\operatorname{div} u = 0$ ，而仅满足无压（或压力空间均匀）的输运型方程：

\partial_t u_k + \sum_{i=1}^n u_i \partial_{x_i} u_k = 0, \quad k = 1,\dots,n.

该方程等价于向量场 $$ u $$ 沿其自身流的李导数为零，即 $\mathcal{L}_u u = 0$ ，亦即 $$ u $$ 是其自身流的守恒量；它本质上是特征线方程 $\dot{X}(t) = u(t,X(t))$ 的伴随输运方程，对应于速度场在拉格朗日坐标下恒定的情形。因此，该模型并非标准Euler方程，而是其一个退化子系统：压力梯度项被完全消去，且不可压约束被主动放弃。作者坦承：‘I am not entirely convinced that this requirement [incompressibility] is really a basic characteristic of fluids in nature’，暗示其物理动机在于质疑不可压假设的普适性，并探索可压但压力空间齐次的动力学结构。

研究对象的核心是一个二元超越函数 $\Omega(x,y)$ ，定义为：

\Omega(x, y) := \frac{y}{x} - W\left( -\frac{1}{x} e^{y/x} \right),

其中 $$ W $$ 为Lambert $$ W $$ 函数的主支，满足 $w = W(z) \iff z = w e^w$ ，定义域为 $z \geq -e^{-1}$ 。该函数并非初等函数的复合，而是超越方程 $e^\omega = x \omega - y$ 的隐式解——这一点由引理中恒等式 $e^{\Omega(x,y)} = x \Omega(x,y) - y$ 严格确立。因此， $\Omega$ 本质上是标量非线性双曲型PDE的隐式特征解：它编码了直线族 $\omega \mapsto x\omega - y$ 与指数曲线 $\omega \mapsto e^\omega$ 之交点的光滑分支。其定义域由Lambert $$ W $$ 的分支点决定：

\operatorname{Dom}(\Omega) = \left\{ (x,y) \in \mathbb{R}^2 : x < 0 \; \text{or} \; \left( x > 0 \; \text{and} \; y \leq x \ln\frac{x}{e} \right) \right\},

该区域在 $$ (x,y) $$ 平面上呈分片结构，包含左半平面全集及右半平面中位于抛物线型边界 $y = x \ln(x/e)$ 下方的开区域，边界上 $$ W $$ 取值 $$ -1 $$ ，对应分支点。此定义域限制直接导致后续速度场 $u_k(t,x) = \Omega(t, x_k)$ 的时空定义域为各坐标方向定义域的交：

\operatorname{Dom}(u) = \bigcap_{k=1}^n \left\{ (t,x) \in \mathbb{R}^{n+1} : t < 0 \; \text{or} \; \left( t > 0 \; \text{and} \; x_k \leq t \ln\frac{t}{e} \right) \right\}.

该区域在时空 $$ (t,x) $$ 中非凸、非柱状，且随维度 $$ n $$ 增长而急剧收缩，例如当 $n \geq 2$ 时，对任意固定 $$ t>0 $$ ，允许的 $$ x $$ 必须满足所有分量 $x_k \leq t \ln(t/e)$ ，故为一无界但高度各向异性的楔形区域。问题的重要性在于：尽管千年难题级别的全局正则性仍遥不可及，但构造任何非平凡、非稳态、非线性且具有明确定义域的解析解，均有助于检验数值格式、理解奇性形成机制、并为能量估计提供反例或基准。本文填补的空白恰在于：首次将Lambert $$ W $$ 函数系统性地嵌入流体力学PDE的解结构中，揭示了一类由超越函数主导的、具有精确隐式特征的输运动力学，其存在性不依赖于小初值、尺度不变性或对称约化，而源于代数-超越方程的可解性本身。

2. 主要定理或主要结果

本文证明了两个核心定理，全部建立在函数 $\Omega(x,y)$ 的微分恒等式基础之上。

定理1（Euler型输运解）：设 $u_k(t,x) := \Omega(t, x_k)$ ， $k = 1,\dots,n$ ，定义域为 $\operatorname{Dom}(u) = \bigcap_{k=1}^n \operatorname{Dom}(u_k)$ ，其中 $\operatorname{Dom}(u_k) = \{ (t,x) \in \mathbb{R}^{n+1} : t < 0 \; \text{or} \; (t > 0 \; \text{and} \; x_k \leq t \ln(t/e)) \}$ 。则在 $\operatorname{int}\, \operatorname{Dom}(u)$ 内，向量场 $u = (u_1,\dots,u_n)$ 满足如下 $$ n $$ 个耦合的一阶拟线性偏微分方程：

\frac{\partial u_k}{\partial t} + \sum_{i=1}^n u_i \frac{\partial u_k}{\partial x_i} = 0, \quad k = 1,\dots,n.

该结论的证明直接引用引理中的关键恒等式 $\partial_1 \Omega(x,y) + \Omega(x,y) \, \partial_2 \Omega(x,y) = 0$ ，并注意到 $\partial_{x_i} u_k = \delta_{ik} \, \partial_2 \Omega(t,x_k)$ ，从而 $\sum_i u_i \partial_{x_i} u_k = u_k \partial_{x_k} u_k = \Omega(t,x_k) \, \partial_2 \Omega(t,x_k)$ ，与时间导数 $\partial_t u_k = \partial_1 \Omega(t,x_k)$ 相加即得零。

定理2（质量守恒律的适配密度）：定义标量函数

\rho(t,x) := \prod_{k=1}^n \frac{1}{e^{\Omega(t,x_k)} - t}, \quad \operatorname{Dom}(\rho) := \operatorname{int}\, \operatorname{Dom}(u).

则在 $\operatorname{int}\, \operatorname{Dom}(u)$ 内， $(\rho, u)$ 满足连续性方程（质量守恒）：

\frac{\partial \rho}{\partial t} + \operatorname{div}(\rho u) = 0.

证明中关键使用了 $\Omega$ 的恒等式 $e^{\Omega} = t \Omega - x_k$ （因 $u_k = \Omega(t,x_k)$ ），从而 $e^{\Omega(t,x_k)} - t = t \Omega(t,x_k) - x_k - t$ ，并结合定理1的零输运性质，最终导出 $\partial_t \rho + \langle u, \nabla \rho \rangle = -\rho \, \operatorname{div} u$ ，再经直接计算验证 $\operatorname{div}(\rho u) = \rho \, \operatorname{div} u + \langle u, \nabla \rho \rangle$ ，故成立。值得注意的是，此处 $\rho$ 并非给定的物态方程输出，而是为适配 $$ u $$ 而逆向构造的、保证质量守恒的密度；其正则性依赖于 $e^{\Omega} \neq t$ ，即排除 $\Omega = \ln t$ 的奇性集，该集恰好落在定义域边界 $x_k = t \ln(t/e)$ 上。

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

本文的核心障碍并非传统PDE中的非线性破坏、耗散缺失或奇性传播，而在于：如何从一个超越函数出发，系统性地导出其满足的偏微分关系，并赋予其物理语境。作者的真正创新点有三：

第一，Lambert W函数的PDE适配性重构。作者并未将 $$ W $$ 视为黑箱特殊函数，而是将其嵌入一个具体的超越方程 $e^{\omega} = x \omega - y$ ，并证明其解 $\omega = \Omega(x,y)$ 天然满足一阶线性PDE $\partial_x \Omega + \Omega \, \partial_y \Omega = 0$ 。这一发现的关键在于对 $$ W $$ 函数导数的精细处理：利用 $$ W'(z) = W(z)/(z(1+W(z))) $$ （文中虽未显式写出，但证明中 $\partial_1 \Omega$ 与 $\partial_2 \Omega$ 的推导隐含此式），并结合链式法则与 $\Omega$ 的定义，将复杂的超越导数化归为代数运算。这构成了全文的基石引理，其证明过程本身即是一种方法论创新——将特殊函数的分析性质翻译为PDE的对称性。

第二，时空变量的非对称分离技巧。标准分离变量法要求解形如 $$ u(t,x) = T(t) X(x) $$ ，但此处作者采用 $u_k(t,x) = \Omega(t, x_k)$ ，即对每个分量 $$ k $$ ，仅将时间 $$ t $$ 与第 $$ k $$ 个空间变量 $$ x_k $$ 配对，其余变量 $$ x_j $$ （ $j \neq k$ ）完全不出现。这种‘单变量耦合’策略绕开了多维耦合的复杂性，使每个 $$ u_k $$ 仅依赖于二维子空间 $$ (t,x_k) $$ ，从而可直接复用 $\Omega$ 的二维PDE性质。然而，代价是定义域变为各二维区域的交集，导致高维时空域严重受限，且解不具备旋转或平移对称性。

第三，可压流框架下的逆向密度构造。面对 $$ u $$ 明显不满足 $\operatorname{div} u = 0$ 的事实，作者未强行施加不可压约束，而是转向可压流框架，将 $$ u $$ 视作给定速度场，反求使其满足连续性方程的密度 $\rho$ 。其构造 $\rho = \prod_k (e^{\Omega(t,x_k)} - t)^{-1}$ 绝非随意猜测：它源于对 $\partial_t \rho + \operatorname{div}(\rho u) = 0$ 进行形式分离，假设 $\rho$ 为乘积形式 $\rho = \prod_k \rho_k(x_k)$ ，并利用 $$ u_k $$ 仅依赖 $$ (t,x_k) $$ 的特性，将PDE分解为 $$ n $$ 个独立的常微分方程。每个因子 $\rho_k$ 需满足 $\partial_t \rho_k + u_k \partial_{x_k} \rho_k = -\rho_k \partial_{x_k} u_k$ ，而由 $\Omega$ 的恒等式 $\partial_{x_k} u_k = \partial_2 \Omega$ 及 $e^{u_k} = t u_k - x_k$ ，可解出 $\rho_k \propto (e^{u_k} - t)^{-1}$ 。这一构造体现了对输运算子本质的深刻把握：它将PDE求解问题转化为对已知向量场的几何测度（密度）的适配问题。

4. 与之前工作的比较

本文与既有文献的差异需从数学、物理和方法论三个层面审视。首先，与经典Euler/NS解析解工作的比较：标准解析解（如Taylor-Green涡、Kida涡、Shear flow）均依赖于强对称性（周期性、自相似性、线性剪切）或小参数展开（如WKB近似），其解空间通常属于Sobolev类或解析函数类，定义域为全时空或柱状区域。而本文解定义在高度各向异性的楔形区域 $\operatorname{Dom}(u)$ 内，且其正则性由Lambert $$ W $$ 的分支点决定——在边界 $x_k = t \ln(t/e)$ 上， $\partial_2 \Omega$ 发散，解出现梯度爆炸（ $\partial_{x_k} u_k \to \infty$ ），这与经典解的全局光滑性形成鲜明对比。文中未与任何具体前人工作（如Bardos & Titi 2007关于弱解唯一性，或Constantin 2008关于临界正则性）进行技术性比较，但从结果看，该解无法用于检验这些理论，因其不满足不可压性，且定义域不满足标准能量估计所需的紧支或周期性条件。

其次，与Lambert W应用文献的比较：参考文献[1–3]集中于渐近分析、算法复杂度或物理模型的代数求解（如燃烧速率、酶动力学），均将 $$ W $$ 作为求解超越代数方程的工具，而非构造PDE解的基石。本文首次将 $$ W $$ 提升为PDE解的核心构件，其解 $\Omega$ 不是某个ODE的积分常数，而是PDE本身的原始解函数，其微分性质（如引理中的PDE）是内禀的，而非派生的。这标志着 $$ W $$ 从“计算辅助函数”到“分析基本函数”的范式跃迁，与文中宣称“should have its proper place among the basic elementary functions”形成呼应。

最后，与可压Euler方程工作的比较：本文定理2构造的 $(\rho, u)$ 对，虽满足连续性方程，但未验证其是否满足完整的可压Euler动量方程 $\partial_t (\rho u) + \operatorname{div}(\rho u \otimes u) + \nabla p = 0$ 。由于 $$ u $$ 不满足 $\operatorname{div} u = 0$ ，若要闭合系统，需额外指定状态方程 $p = p(\rho)$ 并验证动量平衡。作者未进行此项工作，故该解目前仅构成可压流的一个‘半解’（half-solution）：具备质量守恒，但动量守恒未验证。相较之下，Serre（2009）或Chen & Perepelitsa（2012）对可压Euler的研究均基于严格的弱解框架与熵条件，而本文解缺乏此类结构，亦未讨论激波、接触间断等典型可压现象。因此，其主要贡献在于揭示了一种新的、由超越函数驱动的解结构存在性，而非提供一个物理上自洽的完整流体模型。从提供内容中无法确认该解是否能推广至粘性情形（即完整Navier-Stokes），亦未讨论其与Clay Millennium Prize Problem中正则性问题的潜在关联。

2. Surface Quasi-Geostrophic 方程中移动刚性区域的局部适定性

Local well-posedness for a moving rigid region in Surface Quasi-Geostrophic equations

作者: Ludovic Godard-Cadillac, Arnab Roy

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23289v1

摘要翻译: 我们引入并分析一类在存在移动刚性障碍物情形下的 Surface Quasi-Geostrophic (SQG) 方程。该模型的动机一方面来自活动标量方程中奇异结构的涡-波型渐近理论，另一方面来自呈现类刚性相干区域的地球物理现象，例如气旋眼或长寿命的大气尘埃云。我们考虑由给定刚性运动生成的时变外部区域中的临界 SQG 方程，并通过一种适配不可渗透性约束的非局部椭圆公式重构速度场。活动标量被假定在刚性区域内及其边界的某个邻域内保持常值，从而产生一种与输运动力学相容的平台结构。对单个移动障碍物，我们在 Sobolev 空间 $$H^k$$ （ $k\geq 4$ ）中建立了经典解的局部适定性，并同时证明了解的唯一性、局部稳定性以及一个爆破判据。分析依赖于在适配坐标系下的重构，将问题化归为固定区域上的问题，结合分数阶椭圆算子的积分表示、正则化程序以及非线性不动点论证。一个核心难点源于 $s=\frac{1}{2}$ 情形下 SQG Biot-Savart 核的临界奇异性：此时在移动边界附近重构速度需要借助交换子估计。我们进一步证明了平台性质的传播，并导出了先验估计，以同时控制标量梯度支集的大小与解的 Sobolev 范数。据我们所知，本工作提供了首套关于带移动刚性障碍物的 SQG 方程的适定性理论，亦是朝向从 SQG 流中收缩的刚体严格推导点涡型动力学迈出的第一步。

Comment: 无

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

1. 研究对象和背景

本文研究的是带移动刚性障碍物的临界表面准地转（Surface Quasi-Geostrophic, SQG）方程，其核心数学模型为在时变外区域 $F(t) := \mathbb{R}^2 \setminus S(t)$ 上定义的非线性输运方程：

\partial_t \omega(t,x) + v(t,x) \cdot \nabla \omega(t,x) = 0, \quad (t,x) \in [0,T) \times F(t),

其中主动标量 $\omega : [0,T) \times \mathbb{R}^2 \to \mathbb{R}$ 在刚性区域 $S(t) \subset \mathbb{R}^2$ 内恒为常数（设为零），即

\omega(t,x) = 0, \quad \forall x \in S(t),

且在 $$S(t)$$ 的某个邻域内亦为常数——该几何假设称为“平台性质”（plateau property），形式化为

R(t) := \inf\left\{ |x - y| : x \in S(t),\, y \in F(t) \cap \operatorname{supp}(\nabla \omega(t,\cdot)) \right\} > 0.

速度场 $$v$$ 由分数阶 Biot–Savart 律重构，对应临界指数 $s = \frac{1}{2}$ ：

v = -\nabla^\perp \Lambda^{-1} \omega, \quad \text{其中 } \Lambda := (-\Delta)^{1/2},

但需满足运动边界上的不可渗透条件：

v(t,x) \cdot n(x) = \dot{\theta}(t)\, (x - h(t))^\perp \cdot n(x) + \dot{h}(t) \cdot n(x), \quad \forall x \in \partial S(t),

其中 $h(t) \in \mathbb{R}^2$ 和 $\theta(t) \in \mathbb{R}$ 分别为刚体中心位置与朝向角， $$n(x)$$ 是 $\partial S(t)$ 上向外单位法向量。为同时满足输运结构与边界约束，作者引入分解型流函数 $\psi = \psi_2 - \psi_1$ ，其中 $\psi_1$ 满足分数椭圆方程

\Lambda \psi_1(t,x) = \omega(t,x) + \Lambda \psi_2(t,x), \quad x \in F(t),

且 $\psi_1|_{S(t)} = 0$ ；而修正项 $\psi_2$ 显式构造为

\psi_2(t,x) = \chi\big(\operatorname{dist}(x,S(t))\big) \left[ \frac{\dot{\theta}(t)}{2} |x - h(t)|^2 - x \cdot \dot{h}^\perp(t) \right],

其中 $\chi \in C^\infty_c([0,\infty);[0,1])$ 是截断函数，满足 $\chi \equiv 1$ 在 $$[0,a/2]$$ 、 $\chi \equiv 0$ 在 $[a,\infty)$ ， $$a>0$$ 为固定宽度。由此导出的速度为 $v = \nabla^\perp \psi_2 - \nabla^\perp \psi_1$ 。该构造确保 $$v$$ 在 $\partial S(t)$ 上自动满足不可渗透条件，且在 $$F(t)$$ 内恢复标准 SQG 关系 $v = -\nabla^\perp \Lambda^{-1} \omega$ 。问题本质在于：在 $$s=1/2$$ 临界情形下，Biot–Savart 核 $K(x) = -c\, x^\perp / |x|^3$ 非局部可积，导致 $\nabla^\perp \Lambda^{-1}$ 在边界附近缺乏经典正则性；而刚体运动使定义域 $$F(t)$$ 时变，破坏平移不变性，阻碍傅里叶方法直接应用；加之不可渗透约束是非局部、非齐次且随时间演化的，使得传统椭圆理论（如 $$L^p$$ 理论或 Schauder 估计）无法直接移植。已有工作仅处理静止边界情形（如 [12,35,47]），而 Euler 方程中移动障碍物的适定性理论（[25,28,29]）依赖于 Laplace 算子的局部 Green 函数 $G_1(x) = \frac{1}{2\pi} \log \frac{1}{|x|}$ 及其光滑核 $K_1(x) = \frac{1}{2\pi} \frac{x^\perp}{|x|^2}$ ，其可积性规避了临界奇异性；SQG 的分数 Green 函数 $G_{1/2}(x) = c/|x|$ 导致核 $$K(x)$$ 具有 $|x|^{-2}$ 奇异性，恰处于 Calderón–Zygmund 理论失效边缘。本文试图填补的关键空白是：建立首个针对临界 SQG 方程、含预设运动轨迹的刚性障碍物的局部经典解理论，为未来严格推导 SQG 点涡动力学（作为收缩刚体极限）提供基础框架——该目标此前因缺乏运动障碍物下的分数椭圆边值问题适定性工具而长期悬置。

2. 主要定理或主要结果

在单个刚性障碍物 $S(t) = R_{\theta(t)}(S_0 - h(0)) + h(t)$ 、临界指数 $s = \frac{1}{2}$ 、平台性质初始成立及 $\omega_0 \in H^k(\mathbb{R}^2)$ （ $k \geq 4$ ）的假设下，本文证明了如下主定理（Theorem 1.1）：(i) 局部存在唯一经典解 $\omega \in C^1([0,T) \times \mathbb{R}^2)$ ，且 $\omega \in X_T^k$ ，其中 Banach 空间范数定义为 $\|\omega\|_{X_T^k} := \sup_{t \in [0,T)} \|\omega(t,\cdot)\|_{H^k(F(t))}$ ；(ii) 爆破判据：若最大存在时间 $T^* < \infty$ ，则必有 $\int_0^{T^*} \|\omega(t,\cdot)\|_{H^2(F(t))} \, dt = +\infty$ ；(iii) 局部稳定性：对另一初值 $\omega_0^\dagger$ 对应的解 $\omega^\dagger$ ，存在依赖于 $\omega_0$ 、 $$T$$ 及障碍物运动 $h,\theta$ 的常数 $$C > 0$$ ，使得当 $\|\omega_0 - \omega_0^\dagger\|_{L^2(F(0))}$ 足够小时，

\|\omega(t,\cdot) - \omega^\dagger(t,\cdot)\|_{L^2(F(t))} \leq C \, \|\omega_0 - \omega_0^\dagger\|_{L^2(F(0))}, \quad \forall t \in [0,T).

此外，论文还证明了平台性质的传播性：若 $$R(0) > 0$$ ，则存在 $$T > 0$$ 使得 $R(t) \geq R(0)/2 > 0$ 对所有 $t \in [0,T)$ 成立；并建立了先验估计： $\|\omega(t,\cdot)\|_{H^k(F(t))} \leq C(T, \omega_0, h, \theta) \exp\big(C(T, \omega_0, h, \theta) \int_0^t \|\nabla v(s,\cdot)\|_{L^\infty(F(s))} ds\big)$ ，其中关键控制量 $\|\nabla v\|_{L^\infty}$ 通过 $\|\omega\|_{H^2}$ 与平台宽度 $$R(t)$$ 的下界联合控制。

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

证明的核心障碍有三重：一是临界分数核 $K(x) = -c x^\perp / |x|^3$ 在 $\partial S(t)$ 附近导致 $\nabla^\perp \Lambda^{-1} \omega$ 缺乏 $L^\infty$ 控制，传统奇异积分估计失效；二是时变区域 $$F(t)$$ 破坏傅里叶分析与标准能量估计框架；三是不可渗透条件耦合了刚体运动学与流体动力学，使速度重构成为非局部、非自治的边值问题。作者的真正创新点在于提出一套协同技术体系：首先，通过拉格朗日坐标变换 $y = R_{-\theta(t)}(x - h(t))$ 将时变外区域 $$F(t)$$ 映射至固定参考域 $\widetilde{F} = \mathbb{R}^2 \setminus \widetilde{S}$ ，将原问题转化为固定域上的非自治输运方程，代价是引入额外的旋转与平移对流项（见(11)式）。其次，针对分数椭圆算子 $\Lambda$ 在刚体邻域的奇异性，作者设计了基于平台性质的“几何分离”策略：利用 $$R(t) > 0$$ 将 $\operatorname{supp}(\nabla \omega(t,\cdot))$ 与 $\partial S(t)$ 严格分隔，从而在 $\nabla^\perp \Lambda^{-1} \omega$ 的积分表示

\nabla^\perp \Lambda^{-1} \omega(x) = \text{p.v.} \int_{F(t)} K(x-y) \omega(y) \, dy

中，对 $$x$$ 接近 $\partial S(t)$ 的情形，积分核 $$K(x-y)$$ 在 $y \in \operatorname{supp}(\nabla \omega)$ 上远离奇点，得以应用加权 Hardy–Littlewood 极大函数控制。第三，为处理临界情形下 $\nabla^\perp \Lambda^{-1}$ 的弱连续性，作者发展了一套精细的交换子估计：对截断函数 $\chi$ 与 $\Lambda^{-1}$ 的复合，证明了

\| [\chi, \Lambda^{-1}] \nabla \omega \|_{L^2} \leq C \|\nabla \chi\|_{L^\infty} \|\omega\|_{L^2},

该估计依赖于 $\chi$ 的支集宽度 $$a$$ 与平台宽度 $$R(t)$$ 的比值，并成为控制 $\nabla v$ 的 $$L^2$$ -模的关键环节。最后，迭代构造采用正则化-紧性-不动点三重框架：对 $\varepsilon > 0$ ，用磨光算子 $J_\varepsilon$ 正则化 $\omega$ ，求解正则化输运方程，利用平台性质获得一致 $$H^k$$ 先验界，再通过 Aubin–Lions 引理提取强收敛子列，最终在 $C([0,T]; H^{k-1})$ 中验证不动点映射 $\mathcal{T}: \omega \mapsto \widetilde{\omega}$ 的压缩性，其中 $\widetilde{\omega}$ 由正则化速度场输运生成。整个框架的枢纽在于平台性质的传播证明：通过特征线追踪与距离函数演化方程 $\frac{d}{dt} \operatorname{dist}(X(t), S(t)) = -v(X(t)) \cdot n_{S(t)}(X(t)) + O(|\dot{h}| + |\dot{\theta}|)$ ，结合不可渗透条件与平台初始宽度，导出 $$R(t)$$ 的微分不等式 $\frac{d}{dt} R(t) \geq -C \|v\|_{W^{1,\infty}} R(t)$ ，从而保证 $$R(t)$$ 在短时内保持正性。

4. 与之前工作的比较

本文与既有工作的根本推进体现在三个维度：首先是问题设定维度。此前关于 SQG 的边界问题仅限于静止刚性边界（如 [12,35,47]），其分析依赖于 $\Lambda$ 在固定外区域上的谱理论与 $$L^p$$ 有界性；而本文首次处理预设运动轨迹的刚体，彻底摆脱了静态几何的对称性假设。相较 Euler 方程中移动障碍物工作（[25,28,29]），本文的突破在于克服了分数阶算子的临界奇异性：Euler 的 $K_1(x) \sim |x|^{-1}$ 属于弱 $$L^2$$ ，允许直接使用 Calderón–Zygmund 不等式；SQG 的 $K(x) \sim |x|^{-2}$ 则要求更精细的交换子估计与几何分离，这是本文方法论的核心差异。其次是技术工具维度。文献 [16,2] 建立了分数椭圆算子在静止 Lipschitz 区域上的 $C^{\alpha}$ 估计，但未处理时变边界与不可渗透约束；本文通过自适应坐标变换与平台性质，将临界正则性问题转化为固定域上的可控扰动，此策略在 SQG 文献中属首次。再次是物理建模维度。本文明确区分了“障碍物运动给定”与“障碍物运动由流体反作用力决定”两类模型：前者（本文）为点涡极限提供可控的中间步骤，后者（如 Euler 中的 (6)–(7) 式）在 SQG 中因非局域性而尚未建模；作者指出，若强行引入类似 (6)–(7) 的动力学，需额外处理压力 $$p$$ 的非局域表达式，这会引发更高阶的奇异性，目前尚无数学工具支撑。本文的代价在于：解的存在时间 $$T$$ 依赖于初始平台宽度 $$R(0)$$ 与障碍物运动 $C^\infty$ 范数，且 $$H^k$$ 正则性要求 $k \geq 4$ （高于 Euler 中 $k \geq 2$ 的典型要求），反映了临界分数阶对初始数据光滑性的更高需求。此外，爆破判据仅控制 $$H^2$$ 范数而非更低阶范数，暗示 $$H^2$$ 可能是临界正则性阈值——这与 SQG 在全空间中的已知结果（[11]）一致，但本文将其推广至运动障碍物背景。从适用范围看，本文结果严格限定于单障碍物与 $$s=1/2$$ ，而作者明确指出：多障碍物情形需处理障碍物间相互作用的长程耦合，超临界 $$s > 1/2$$ 可沿用本文框架，但亚临界 $$s < 1/2$$ 需引入低频-高频分解（类 [11] 方法），这些延伸均未在本文中展开。

3. 单速度Baer-Nunziato模型的压力与温度松弛极限

Pressure and temperature relaxation limit for a one-velocity Baer-Nunziato model

作者: Cosmin Burtea, Timothée Crin-Barat, Pierre Gonin–Joubert

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23693v1

摘要翻译: 处于力学与热力学非平衡态的两相流动力学由一个具刚性相互作用项的一阶拟线性偏微分方程组描述，该系统为部分耗散型，且其相互作用项对应于快速松弛尺度。本文从数学角度分析单速度Baer-Nunziato模型所导出的压力与温度松弛奇异极限问题。该问题涉及两个小参数。我们为此系统提出一种一致对称化方法，从而严格证明强松弛极限的存在性，并为经典解建立收敛速率估计。

Comment: 无

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

1. 研究对象和背景

本文研究的是一类描述非机械平衡与非热平衡两相流的双曲-耗散耦合系统，其核心数学对象是Kapila等人在[25]中导出的单速度、双压、双温Baer-Nunziato（BN）模型。该模型由如下一阶拟线性偏微分方程组给出：

\begin{cases} \partial_t \alpha_\pm + u \cdot \nabla \alpha_\pm = \pm \dfrac{\alpha_+ \alpha_-}{\mu \varepsilon} (p_+ - p_-), \\ \partial_t (\alpha_\pm \rho_\pm) + \mathrm{div}(\alpha_\pm \rho_\pm u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) + \nabla(\alpha_+ p_+ + \alpha_- p_-) = 0, \\ \partial_t (\alpha_\pm \rho_\pm e_\pm) + \mathrm{div}(\alpha_\pm \rho_\pm e_\pm u) + \alpha_\pm p_\pm \, \mathrm{div}\, u = \mp \dfrac{\alpha_+ \alpha_-}{\mu \varepsilon} p_- (p_+ - p_-) \mp \dfrac{\theta_+ - \theta_-}{\mu}. \end{cases} \tag{BN}

其中状态变量为： $\alpha_\pm \in (0,1)$ 为体积分数且满足 $\alpha_+ + \alpha_- = 1$ ； $\rho_\pm > 0$ 为相密度； $u \in \mathbb{R}^d$ 为公共速度； $p_\pm > 0$ 为相压力； $e_\pm > 0$ 为内能； $\theta_\pm > 0$ 为温度。物理闭包取理想气体物态方程 $p_\pm = R_\pm \rho_\pm \theta_\pm$ 及线性内能关系 $e_\pm = c_{v\pm} \theta_\pm$ ，其中绝热常数满足 $\gamma_\pm = 1 + R_\pm / c_{v\pm}$ 且 $\gamma_+ \neq \gamma_-$ （见(1.1)式）。参数 $\mu, \varepsilon > 0$ 表征松弛时间尺度： $\mu \varepsilon$ 对应压力弛豫时间， $\mu$ 对应温度弛豫时间；文献[25]指出前者远快于后者，故 $\varepsilon \ll 1$ 固定，而 $\mu \to 0$ 是本文所考察的奇异极限过程。

该模型的数学挑战根植于其结构性缺陷：(i) 方程含非守恒乘积项（如 $u \cdot \nabla \alpha_\pm$ ），源于界面交换；(ii) 自然熵泛函不严格凸（参见[16, 7, 10]），导致经典Lax-Friedrichs对称化定理失效；(iii) 系统同时含两个小参数 $\mu, \varepsilon$ ，构成双重奇异极限问题，其刚性源项 $\mathcal{O}(1/(\mu \varepsilon))$ 和 $\mathcal{O}(1/\mu)$ 在常规能量估计中引发灾难性增长；(iv) 极限系统(K)为退化守恒律：

\begin{cases} \partial_t (\alpha_\pm \rho_\pm) + \mathrm{div}(\alpha_\pm \rho_\pm u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) + \nabla p = 0, \\ \partial_t (\alpha_+ \rho_+ e_+ + \alpha_- \rho_- e_-) + \mathrm{div}((\alpha_+ \rho_+ e_+ + \alpha_- \rho_- e_-)u) + p \, \mathrm{div}\, u = 0, \\ p = R_+ \rho_+ \theta_+ = R_- \rho_- \theta_-, \quad \theta_+ = \theta_-, \quad e_\pm = c_{v\pm} \theta_\pm, \quad \rho = \alpha_+ \rho_+ + \alpha_- \rho_-. \end{cases} \tag{K}

(K)中压力与温度已统一，但其变量间存在强代数约束（如 $p = R_+ \rho_+ \theta_+$ 与 $p = R_- \rho_- \theta_-$ 联立隐含 $\rho_+, \rho_-, \theta_+$ 非独立），致使(K)本身并非标准双曲系统，而需通过变量约化（如引入质量分数 $y_\pm = \alpha_\pm \rho_\pm / \rho$ 与平均量 $\langle \cdot \rangle = y_+ (\cdot)_+ + y_- (\cdot)_-$ ）重构为六维系统。已有工作如Yong[33,34]为单参数松弛极限建立了对称化框架，要求系统具备严格凸熵及特定零阶项结构，但(BN)因熵退化与非守恒性，无法直接满足其假设。Crin-Barat–Danchin[12]、Danchin[13]等虽处理了多类松弛系统，但均未覆盖(BN)这种双参数、非守恒、退化熵的复合结构。因此，本文旨在填补这一关键空白：构建一个对 $\mu, \varepsilon$ 一致有效的对称化框架，从而严格证明当 $\mu \to 0$ 时，(BN)的经典解强收敛至(K)的解，并给出定量收敛速率。

2. 主要定理或主要结果

本文在 Sobolev 空间 $H^s(\mathbb{R}^d)$ （ $$s > d/2 + 1$$ ）框架下，对初值 $U^\mu(0) = U_0$ 满足 $\|U_0\|_{H^s} \leq M$ 且位于相空间开集 $\mathcal{U}_\eta$ 内（即 $|\delta p| < \eta$ , $|\delta \theta| < \eta$ ，见(2.46)式）的情形，证明了以下核心结论：

定理 1（一致局部适定性与对称化）：存在 $$T > 0$$ 与 $\eta > 0$ ，使得对任意 $\mu \in (0, \mu_0]$ , $\varepsilon \in (0, \varepsilon_0]$ ，系统(BN)在时间区间 $$[0, T]$$ 上存在唯一经典解 $U^\mu \in C([0,T]; H^s) \cap C^1([0,T]; H^{s-1})$ ，且其 $$H^s$$ 范数关于 $\mu, \varepsilon$ 一致有界： $\sup_{t \in [0,T]} \|U^\mu(t)\|_{H^s} \leq C(M)$ 。该结论依赖于作者构造的显式对称化变换 $\Phi: (\alpha_+, \rho_\pm, \theta_\pm) \mapsto (U_1, U_2, w, y_+, \rho)$ ，其中新变量定义为：

w := \frac{\langle R \rangle}{\langle c_v \rangle} \langle e \rangle, \quad U_1 := f \, \delta p, \quad U_2 := \frac{p}{\rho} - w - U_1,

其中 $f = \dfrac{\gamma_+ - \gamma_-}{\rho} \left( \dfrac{\gamma_+}{\alpha_+} + \dfrac{\gamma_-}{\alpha_-} \right)$ ， $\delta p = p_+ - p_-$ ， $\delta \theta = \theta_+ - \theta_-$ ， $\langle \cdot \rangle$ 为质量平均。在此坐标下，(BN)被重写为对称形式 $\partial_t U + \sum_{j=1}^d H_j(U) \partial_{x_j} U + D_{\varepsilon,\mu} U = R_{\varepsilon,\mu}$ ，其中 $D_{\varepsilon,\mu} = \mathrm{diag}\left( \dfrac{\langle R/\gamma \rangle}{\langle R \rangle \gamma_+ \gamma_- \rho w}{\mu \varepsilon}, \dfrac{\langle c_v \rangle \langle R \rangle}{\langle R/\gamma \rangle} \dfrac{1}{\rho c_{v+} c_{v-} y_+ y_- \gamma_+ \gamma_-} \dfrac{1}{\mu}, 0, 0, 0, 0 \right)$ ，且对称矩阵 $S = \mathrm{diag}(a_1^{-1}, a_2^{-1}, a_3^{-1}, 1, p/\rho^2, \rho)$ 满足 $$S H_j$$ 对称（见(2.23)式）。

定理 2（强松弛极限与收敛速率）：设 $U^\mu$ 为(BN)的解， $$U^0$$ 为极限系统(K)在相同初值 $U^0(0) = \Phi(U_0)$ 下的解，则存在常数 $$C > 0$$ 使得

\sup_{t \in [0,T]} \|U^\mu(t) - U^0(t)\|_{H^{s-1}} \leq C \left( \|U_0^\mu - U_0^0\|_{H^{s-1}} + \sqrt{\mu} \right),

且若初值更光滑（ $U_0 \in H^{s+1}$ ），则可提升至

\sup_{t \in [0,T]} \|U^\mu(t) - U^0(t)\|_{H^{s-2}} \leq C \left( \|U_0^\mu - U_0^0\|_{H^{s-2}} + \mu \right).

该速率 $\mathcal{O}(\sqrt{\mu})$ （ $H^{s-1}$ 范数）与 $\mathcal{O}(\mu)$ （ $H^{s-2}$ 范数）是本文的核心定量结果，优于Yong[33]中仅在首阶修正项获得的衰减率。

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

本文的核心技术障碍在于：(BN)系统既非守恒型，其自然熵又不严格凸，且含双重小参数刚性源项，导致标准对称化方法（如Lax-Friedrichs或Yong条件）完全失效。作者的真正创新点在于提出了一种“松弛自适应对称化”（relaxation-adapted symmetrization）框架，其关键步骤与工具如下：

第一步：物理驱动的变量分解与快慢分离。作者观察到质量分数 $y_\pm$ 满足无源输运方程 $\partial_t y_\pm + u \cdot \nabla y_\pm = 0$ ，而平均内能 $\langle e \rangle$ 的演化（见(2.24)式）不含奇异项： $\partial_t \langle e \rangle + u \cdot \nabla \langle e \rangle + (p/\rho) \, \mathrm{div}\, u = 0$ 。这提示将 $w = \langle R \rangle / \langle c_v \rangle \langle e \rangle$ （其物理量纲与 $p/\rho$ 相同）作为慢变量。剩余部分 $p/\rho - w$ 则承载了非平衡效应，作者将其精确展开为 $p/\rho - w = (\gamma_+ - \gamma_-) c_{v+} c_{v-} / \langle c_v \rangle \, y_+ y_- \, \delta \theta$ （见(2.31)式），表明其大小与温度差 $\delta \theta$ 同阶。此分解天然地将系统解耦为慢变量 $(w, y_+, \rho, u)$ 与快变量 $$(U_1, U_2)$$ ，其中 $$U_1, U_2$$ 被刻意设计为 vanish as $\delta \theta \to 0$ ，从而在极限 $\mu \to 0$ 下自动趋于零。

第二步：奇异项的代数抵消机制。这是全文最精巧的技术。作者发现，若直接对 $p/\rho - w$ 或 $\delta p$ 建立方程，其右端会出现无法控制的 $\mathcal{O}(1/(\mu \varepsilon))$ 项（见(2.32)(2.33)式）。为消除此障碍，作者构造了两个辅助变量： $U_1 = f \, \delta p$ 与 $U_2 = p/\rho - w - U_1$ ，其中系数 $$f$$ 的选取（见(2.36)式）并非任意，而是通过符号计算（SageMath）严格求解出的代数条件——它必须使得 $$U_2$$ 的演化方程中 $\delta p/(\mu \varepsilon)$ 项被完全抵消。这一设计将原本分散在多个方程中的刚性耗散集中到 $$U_1$$ 的方程中（其耗散率 $\sim 1/(\mu \varepsilon)$ ），而 $$U_2$$ 则获得一个更温和的 $\mathcal{O}(1/\mu)$ 耗散率（见(2.37)式第二行）。这种“耗散分配”是实现一致估计的关键。

第三步：全局微分同胚与一致对称化。为确保变量变换 $\Phi$ 可逆且保持正则性，作者证明了其Jacobi行列式在 $\delta p = \delta \theta = 0$ 处非零（见(2.47)式），并利用引理2.1（基于隐函数定理的局部扰动引理）得出：只要初值满足 $|\delta p| < \eta$ , $|\delta \theta| < \eta$ ，则 $\Phi$ 在邻域 $\mathcal{U}_\eta$ 上是 $C^\infty$ 微分同胚。这一全局可逆性保证了新坐标系下的能量估计可无损地回传至原变量，从而得到关于 $\mu, \varepsilon$ 一致的 $$H^s$$ 先验界。整个证明框架不依赖于任何紧性论证或弱解技术，而是建立在经典解的 $$L^2$$ -能量估计之上，其核心引理是针对对称化系统 $S \partial_t U + S \sum_j H_j \partial_{x_j} U + S D_{\varepsilon,\mu} U = S R_{\varepsilon,\mu}$ 所做的加权能量估计，其中对角对称矩阵 $$S$$ 的选取（见(2.23)式）确保了耗散项 $\int (S D_{\varepsilon,\mu} U) \cdot U$ 产生正定贡献 $\gtrsim \mu^{-1} \|U_1\|_{L^2}^2 + \mu^{-1} \|U_2\|_{L^2}^2$ ，而源项 $S R_{\varepsilon,\mu}$ 则被 $\|U_1\|_{L^2} + \|U_2\|_{L^2}$ 控制，从而通过Grönwall不等式封闭估计。

4. 与之前工作的比较

本文与既有工作的推进主要体现在四个维度：

1. 与Yong[33,34]工作的比较：Yong为单参数松弛极限（如 $\partial_t V + A(V) \partial_x V + Q(V)/\mu = 0$ ）建立了系统性理论，其核心是要求系统存在严格凸熵及满足特定代数条件的对称化矩阵。然而，(BN)模型因其非守恒结构与退化熵，无法满足Yong的假设。本文通过构造新的变量 $(U_1, U_2, w, y_+, \rho, u)$ ，本质上绕过了对原始熵的要求，转而利用物理驱动的耗散结构（ $\delta p$ 与 $\delta \theta$ 的显式耗散，见(2.14)式）来构建对称化，从而将Yong的框架推广到了非守恒、退化熵的双参数情形。代价是：本文的对称化是高度问题特异的，其变量变换 $\Phi$ 依赖于 $\gamma_\pm, R_\pm, c_{v\pm}$ 等具体物性参数，不具备Yong理论的普适性。

2. 与Crin-Barat–Danchin[12]及Danchin[13]工作的比较：这些工作处理了如多孔介质流、非平衡气体等松弛系统，但其模型通常具有守恒形式或至少具备一个严格凸熵。例如，[12]中处理的系统可写为 $\partial_t W + \mathrm{div} F(W) + G(W)/\mu = 0$ ，其守恒结构允许使用标准熵方法。而(BN)的非守恒性（如 $u \cdot \nabla \alpha_\pm$ 项）迫使本文发展出一套全新的代数技巧来处理非守恒乘积，其核心是将非守恒项吸收进新变量的对流部分（见(2.15)式），而非像[12]那样直接估计。适用范围上，本文覆盖了更复杂的两相流模型，但时间尺度受限于局部适定性（ $$T$$ 依赖于初值大小），而[13]在某些情形下获得了全局解。

3. 与Kawashima–Xu[32]及Beauchard–Zuazua[3]工作的比较：这些工作侧重于全局存在性，常依赖于精细的频谱分析或Lyapunov函数。本文则聚焦于强极限的定量刻画，其收敛速率 $\mathcal{O}(\sqrt{\mu})$ 在 $H^{s-1}$ 中是全新的。Yong[33]的速率仅在首阶修正项 $U^\mu = U^0 + \mu U^1 + \cdots$ 中出现，而本文直接给出了原始解的差 $U^\mu - U^0$ 的速率，且该速率源于快变量 $$U_1, U_2$$ 的本征耗散（ $\|U_i(t)\|_{L^2} \lesssim e^{-c t / \mu} \|U_i(0)\|_{L^2} + \sqrt{\mu}$ ），这比[32]中依赖于色散效应的衰减更为直接。

4. 与Baer–Nunziato模型自身文献的比较：此前关于(BN)的工作，如[17, 9, 28]，主要集中于单速度情形的局部适定性或Riemann问题，但均未处理双参数松弛极限。本文首次将压力与温度松弛统一纳入一个数学框架，并揭示了二者的时间尺度分离（ $\mu \varepsilon$ vs $\mu$ ）是构造一致对称化的必要条件（见Remark 2.2）。从提供内容中无法确认本文结果是否可推广至两速度BN模型，因作者明确指出单速度情形的保守性（可恢复守恒形式）是其方法的关键前提。

4. 一类具有超二次非线性的非局部问题的山路算法

A Mountain-Pass Algorithm for Nonlocal Problems with Super-quadratic Nonlinearities

作者: Loic Cappanera, Gabriela Jaramillo, Joshua M. Siktar

类别: math.AP

发布日期: 2026-05-21

arXiv ID: 2605.23106v1

摘要翻译: 本文考虑一类非线性方程 $-\mathcal{L} u(x) = f(x, u(x))$ ，其中非线性项 $$f$$ 具有超二次增长性，且非局部算子 $\mathcal{L}$ 由一类特殊的径向对称 $$L^1$$ 卷积核生成，该类核具有有限的二阶矩。对该算子所作的假设较为宽松，允许涵盖生物与物理应用中广泛使用的多种核函数，包括代数衰减型核与变号核。

Comment: 34 pages

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

1. 研究对象和背景

本文研究定义在有界光滑区域 $\Omega \subset \mathbb{R}^n$ 上的非局部非线性方程

-\mathcal{L} u(x) = f(x, u(x)), \quad x \in \Omega,

其中非局部算子 $\mathcal{L}$ 由径向对称、 $$L^1$$ 可积且具有有限二阶矩的卷积核 $\gamma(|x-y|)$ 生成：

\mathcal{L}u(x) = \int_{\mathbb{R}^n} \big(u(y) - u(x)\big) \gamma(|x - y|) \, dy.

该算子不具奇异性（即核函数无原点奇点），区别于分数阶拉普拉斯 $(-\Delta)^s$ ；其傅里叶符号 $\widehat{\mathcal{L}}(\xi) = \widehat{\gamma}(\xi) - 1$ 满足 Assumption 2.1 与 2.3：存在 $\xi_0 > 0$ ，使得 $\widehat{\mathcal{L}}(\xi)$ 在复平面条带 $|\operatorname{Re}(\xi)| < \xi_0$ 上解析且一致有界，并在 $|\xi| \to 0$ 时具有渐近展开

\widehat{\mathcal{L}}(\xi) = -|\xi|^2 + O(|\xi|^4).

这一结构保证 $\mathcal{L}$ 在低频端等价于经典拉普拉斯 $-\Delta$ ，但允许高频行为自由变化，从而涵盖代数衰减核（如 $\gamma(z) \sim |z|^{-n-\alpha},\, \alpha \in (0,2)$ ）与变号核（如神经场模型中的兴奋-抑制耦合）。关键的是，作者将一般核 $\gamma$ 分解为扩散主部 $\gamma_e$ 与小扰动部 $\gamma_s$ ：

\gamma = \gamma_e + \gamma_s, \quad \|\gamma_s\|_{L^1} - \Gamma_s < \beta,

其中 $\Gamma_s = \int \gamma_s$ ，而 $\beta > 0$ 是主部双线性型 $B_e[u,u] = -\langle \mathcal{L}_e u, u \rangle_{L^2(\Omega)}$ 的最优 coercivity 常数。该分解使算子可写为

\mathcal{L}u = \Gamma_e (\tilde{\gamma}_e * u - u) + \gamma_s * u - \Gamma_s u,

从而将非局部问题嵌入到 $H^1(\Omega)$ 框架中——这是克服 $L^2(\Omega)$ 上缺乏紧嵌入的根本策略。边界条件采用两类齐次非局部形式：(i) 非局部 Dirichlet 条件 $$u = 0$$ on $\Omega^c = \mathbb{R}^n \setminus \Omega$ ，对应函数空间 $X_D^1 = \{ u \in H^1(\Omega) : u|_{\Omega^c} = 0 \}$ ；(ii) 非局部 Neumann 条件 $\mathcal{L}u = 0$ on $\Omega^c$ ，对应空间 $X_N^1 = \{ u \in H^1(\mathbb{R}^n) : \mathcal{L}u = 0 \text{ on } \Omega^c \}$ ，其范数含非局部半范 $|u|_{B_N}^2 = \frac{1}{2} \iint_{\mathbb{R}^n \times \mathbb{R}^n} (u(y)-u(x))^2 \gamma(|x-y|) \, dy dx$ 。非线性项 $$f(x,u)$$ 被假设为超二次的：存在 $$p > 2$$ 及 $$C_f > 0$$ ，使得对所有 $x \in \Omega$ , $u \in \mathbb{R}$ ，有

|f(x,u)| \leq C_f \big( |u| + |u|^{p-1} \big), \quad \text{且} \quad 0 < \theta F(x,u) \leq u f(x,u) \quad \text{for } |u| \geq R_0,

其中 $F(x,u) = \int_0^u f(x,s)\,ds$ ， $\theta > 2$ 。该 Ambrosetti-Rabinowitz 条件是构造山路几何的关键，它确保能量泛函 $J(u) = \frac{1}{2} B_{D/N}[u,u] - \int_\Omega F(x,u)\,dx$ 满足“山谷-山脊”结构： $$J(0)=0$$ ， $J(u) \geq \alpha > 0$ 当 $\|u\|_{H^1} = \rho$ 小时，而存在 $e \in X_{D/N}^1$ 满足 $\|e\|_{H^1} > \rho$ 且 $$J(e) < 0$$ 。研究动机源于生物与物理建模：种群扩散模型（ $\gamma$ 控制种子传播范围）、神经场模型（变号核刻画兴奋/抑制突触）、相变模型（非局部 Landau 泛函 $\mathcal{E}[u] = \frac{1}{2} \iint (u(y)-u(x))^2 \gamma(|x-y|) dy dx - \int W(u) dx$ ）均导出此类稳态方程。已有工作多聚焦于分数阶情形（如 $\mathcal{L} = (-\Delta)^s$ ）或紧支核，而本文填补了对一般 $$L^1$$ 、有限二阶矩、可能变号核的非局部算子，在超二次非线性下系统应用极小极大理论的空白——尤其在于如何在缺乏标准 Sobolev 紧嵌入（因 $\mathcal{L}$ 不产生分数阶正则性）的前提下，通过重构双线性型至 $$H^1$$ 框架并利用核的扩散主导结构来恢复 $$(PS)_c$$ 序列的紧性。

2. 主要定理或主要结果

本文在 Assumptions 2.1, 2.3, 2.11 及 $$f$$ 满足超二次增长与 Ambrosetti-Rabinowitz 条件（即存在 $\theta > 2$ , $$R_0 > 0$$ 使得 $0 < \theta F(x,u) \leq u f(x,u)$ for $|u| \geq R_0$ ）下，证明了以下核心定理：

Theorem 3.1（山路解存在性）：设 $\Omega \subset \mathbb{R}^n$ 为有界光滑区域， $\gamma$ 满足 Assumption 2.11，则对非局部 Dirichlet 或 Neumann 边界条件，能量泛函

J(u) = \frac{1}{2} B_{D/N}[u,u] - \int_\Omega F(x,u)\,dx, \quad u \in X_{D/N}^1,

满足山路几何结构与 Palais-Smale 条件 $$(PS)_c$$ （即任意满足 $J(u_k) \to c$ 且 $J'(u_k) \to 0$ in $(X_{D/N}^1)'$ 的序列 $\{u_k\}$ 在 $X_{D/N}^1$ 中存在强收敛子列）。因此，存在非平凡临界点 $u_0 \in X_{D/N}^1 \setminus \{0\}$ 满足 $$J'(u_0) = 0$$ ，即 $-\mathcal{L} u_0 = f(x, u_0)$ 在 $\Omega$ 中弱成立。

Theorem 4.1（改进存在性与数值基础）：若进一步假设 $$f(x,u)$$ 关于 $$u$$ 为奇函数（ $$f(x,-u) = -f(x,u)$$ ）且满足更强的单调性条件 $\partial_u f(x,u) > 0$ for $$u > 0$$ ，则存在一个径向路径 $\gamma(t) = t v$ （ $v \in X_{D/N}^1$ 固定）使得 $$J$$ 沿此路径有唯一最大值点 $$t_v > 0$$ ，且该点显式可表为

t_v = \left( \frac{B_{D/N}[v,v]}{\int_\Omega f(x,v) v \, dx} \right)^{\frac{1}{p-2}} \quad \text{（当 } f(x,u) \sim |u|^{p-2}u \text{ 时）}.

此性质构成第5节数值算法的基础：梯度下降方向被约束在单位球面 $S = \{ u \in X_{D/N}^1 : \|u\|_{H^1} = 1 \}$ 上，并通过求解非局部线性方程 $-\mathcal{L} w = f(x,u)$ 更新法向分量，避免了标准牛顿法在非局部框架下的发散问题。误差分析表明，该算法在 $$H^1$$ 范数下达到 $$O(h^2)$$ 收敛阶（ $$h$$ 为空间离散步长），且迭代次数与问题维数 $$n$$ 呈亚线性关系。

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

本文的核心技术障碍在于：非局部算子 $\mathcal{L}$ 由 $$L^1$$ 核生成，其自然作用空间为 $L^2(\Omega)$ ，但 $L^2(\Omega)$ 到 $L^p(\Omega)$ （ $$p>2$$ ）的嵌入不紧，导致标准变分法所需的 Palais-Smale $$(PS)_c$$ 条件失效——这是应用山路定理的最大障碍。作者的真正创新点并非简单套用已知定理，而是通过三重结构性重构，将问题强行纳入 $H^1(\Omega)$ 框架并恢复紧性：

第一，双线性型域提升：放弃 $$L^2$$ 自然域，明确定义双线性型 $B_{D/N}[u,\phi] = \frac{1}{2} \iint_{\mathbb{R}^n \times \mathbb{R}^n} (u(y)-u(x))\gamma(|x-y|)(\phi(y)-\phi(x)) \, dy dx$ 于 $$H^1$$ 型空间 $X_{D/N}^1$ 上。关键引理（Propositions 2.5, 2.9, 2.12）证明了该双线性型在 $$H^1$$ 范数下是强制的（coercive）与有界的：存在 $\beta_2 > 0$ , $$C_2 > 0$$ 使得

\beta_2 \|u\|_{H^1(\Omega)}^2 \leq B_{D/N}[u,u] \leq C_2 \|u\|_{H^1(\Omega)}^2, \quad \forall u \in X_{D/N}^1.

此强制性直接来自 Assumption 2.3 的傅里叶符号展开 $\widehat{\mathcal{L}}(\xi) = -|\xi|^2(1 + L_1(\xi))$ ，其中 $L_1(\xi) = O(|\xi|^2)$ ，故低频主导项 $-|\xi|^2$ 提供了 $\|\nabla u\|_{L^2}^2$ 等价控制；高频部分由 $$L_1$$ 的解析性保证其作为小扰动不破坏强制性。

第二，核分解与扰动控制：为处理一般 $\gamma$ （含变号、代数衰减），作者引入关键分解 $\gamma = \gamma_e + \gamma_s$ （Assumption 2.11）。主部 $\gamma_e$ 生成满足 Assumptions 2.1 & 2.3 的扩散算子 $\mathcal{L}_e$ ，其双线性型 $$B_e$$ 具有已知强制常数 $\beta$ ；扰动部 $\gamma_s$ 被严格量化为“小”：要求 $\|\gamma_s\|_{L^1} - \Gamma_s < \beta$ 。由此导出

B_{D/N}[u,u] = B_e[u,u] + \underbrace{\int_{\Omega} (\gamma_s * u) u \, dx - \Gamma_s \|u\|_{L^2}^2}_{\text{扰动项}} \geq (\beta - (\|\gamma_s\|_{L^1} - \Gamma_s)) \|u\|_{L^2}^2,

从而将强制性从 $$B_e$$ 传递至 $B_{D/N}$ 。此估计是全文紧性论证的基石，它不依赖于核的正定性，允许 $\gamma_s < 0$ （如抑制性神经连接）。

第三， $$(PS)_c$$ 序列的紧性恢复：给定 $$(PS)_c$$ 序列 $\{u_k\} \subset X_{D/N}^1$ ，由强制性得 $\|u_k\|_{H^1}$ 有界，故存在弱收敛子列 $u_k \rightharpoonup u$ in $$H^1$$ 。难点在于证强收敛。作者利用超二次条件与强制性：由 $J'(u_k) \to 0$ 得

B_{D/N}[u_k, u_k - u] = \int_\Omega f(x,u_k)(u_k - u) \, dx + o(1) \|u_k - u\|_{H^1}.

左边由强制性及 $u_k \rightharpoonup u$ 得 $\liminf B_{D/N}[u_k, u_k - u] \geq \beta_2 \liminf \|u_k - u\|_{H^1}^2$ ；右边利用 $$f$$ 的超二次增长与 Lions 引理（或通过 $H^1 \hookrightarrow L^p$ 紧嵌入，因 $\Omega$ 有界）证得 $\int f(x,u_k)(u_k - u) \to 0$ 。联立即得 $\|u_k - u\|_{H^1} \to 0$ 。此论证完全规避了分数阶空间的复杂性，是方法论上的实质性推进。

4. 与之前工作的比较

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

1. 算子类别的根本拓展：先前关于非局部山路定理的工作（如 [5, 6, 17, 24, 28, 42]）几乎全部针对分数阶拉普拉斯 $(-\Delta)^s$ （ $s \in (0,1)$ ）或其推广（如 $\mathcal{L}u = \int (u(y)-u(x)) |x-y|^{-n-2s} dy$ ）。这些算子具有齐次性、尺度不变性及明确的 Sobolev 嵌入 $W^{s,2} \hookrightarrow L^p$ （ $p \leq 2n/(n-2s)$ ），其 $$(PS)_c$$ 条件可直接由分数阶 Rellich-Kondrachov 紧嵌入保证。而本文处理的是更一般的 $$L^1$$ 卷积算子，其核无需齐次性、无需奇异、可代数衰减或变号——这涵盖了大量实际模型（如神经场、生态扩散），但丧失了分数阶的统一正则性理论。作者通过 $$H^1$$ 框架重构而非 $W^{s,2}$ ，是对此类算子变分分析的首次系统性建立。

2. 边界条件的非局部化处理：经典工作（如 [18, 19, 21]）使用局部 Dirichlet 条件 $u|_{\partial \Omega}=0$ ，其自然空间为 $H_0^1(\Omega)$ ，紧嵌入自动成立。本文的非局部 Dirichlet 条件 $$u=0$$ on $\Omega^c$ 对应空间 $X_D^1 = H^1(\Omega)$ （非 $$H_0^1$$ ），其范数不含边界项；非局部 Neumann 条件 $\mathcal{L}u = 0$ on $\Omega^c$ 更无局部类比，需全新定义空间 $$X_N^1$$ 及其非局部半范 $|\cdot|_{B_N}$ 。Proposition 2.9 与 Lemma 2.10 证明的 Poincaré 不等式 $\|u\|_{L^2}^2 \leq \beta B_N[u,u]$ 是对此类边界条件的首个强制性结果，代价是要求主部核 $\gamma_e > 0$ （Neumann 情形），而 Dirichlet 情形无此限制。

3. 假设的放宽与代价：相比 [1, 2, 18, 19] 中对非线性 $$f$$ 的严格凸性或 $$C^2$$ 光滑性要求，本文仅需连续性与超二次增长（AR 条件），更贴近物理模型（如 $$f(u) = u^3 - u$$ ）。但代价是：(i) 必须假设核具有有限二阶矩（排除 $\alpha \geq 2$ 的超慢衰减核）；(ii) 扰动部 $\gamma_s$ 的“小性”需显式量化（ $\|\gamma_s\|_{L^1} - \Gamma_s < \beta$ ），而 $\beta$ 依赖于主部核的二阶矩 $\alpha = \int |z|^2 \gamma_e(z) dz$ （Remark 2.13 给出 $\beta \sim \alpha / d(\Omega)^2$ ），这意味着在大区域 $\Omega$ 上，允许的扰动幅度更小；(iii) 数值算法（Theorem 4.1）要求 $$f$$ 为奇函数，此限制未在存在性定理中出现，是为获得显式径向最大值点而引入的额外代价。

4. 与分数阶工作的互补性：本文未与 [17, 24, 28] 等分数阶工作进行直接比较，但从方法论看，其优势在于适用性广（涵盖非奇异核），劣势在于无法处理 $$s < 1/2$$ 的强非局部情形（因 $$H^1$$ 框架无法控制过弱的正则性）。作者明确指出（第2页）：“In contrast to other works that model the operator $\mathcal{L}$ as a fractional Laplacian […] our analytical results are based on integral operators […] generated by non-singular, radially symmetric kernels”，这一定位清晰划定了本文的贡献边界：它不是对分数阶理论的改进，而是为另一大类重要非局部算子建立了独立的变分分析范式。

5. 半空间中可压 Navier–Stokes 方程的整体一致正则性与垂直粘性消失极限

Global uniform regularity and vanishing vertical viscosity limit for the compressible Navier–Stokes equations in the half-space

作者: Jincheng Gao, Lianyun Peng, Jiahong Wu, Zheng-an Yao

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23233v1

摘要翻译: 在大地物理流体（如大尺度海洋动力学）中，垂直粘性通常远小于水平粘性。这种各向异性自然引出如下问题：当垂直粘性系数 $\varepsilon \to 0$ 时，全各向异性可压 Navier–Stokes 方程的解是否收敛于一个仅具水平耗散的极限系统之解，且该极限能否在整体时间区间上被严格验证？此前工作仅在局部时间或不可压情形下回答了该问题。本文解决了上半空间中带 Navier 滑移边界条件的三维可压 Navier–Stokes 方程的这一问题。本文对常数平衡态的小扰动建立了两个主要结果：第一，证明存在唯一整体时间解，其余法 Sobolev 范数在所有 $t \ge 0$ 及所有 $\varepsilon \in (0,1)$ 下保持一致有界；第二，严格验证了整体垂直粘性消失极限。更精确地，我们证明该解收敛于水平耗散可压 Navier–Stokes 系统的一个整体解。这首次为无限时间区间上的可压流体各向异性粘性极限提供了严格的数学验证。

Comment: 无

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

1. 研究对象和背景

本文研究三维上半空间 $\mathbb{R}^3_+ := \{x = (x_1, x_2, x_3) : x_3 > 0\}$ 中带Navier滑移边界条件的各向异性可压Navier-Stokes方程组：

\begin{cases} \partial_t \rho^\varepsilon + \operatorname{div}(\rho^\varepsilon u^\varepsilon) = 0, \\ \rho^\varepsilon \partial_t u^\varepsilon + \rho^\varepsilon u^\varepsilon \cdot \nabla u^\varepsilon - \Delta_h u^\varepsilon - \varepsilon \partial_3^2 u^\varepsilon - \nabla \operatorname{div} u^\varepsilon + \nabla P(\rho^\varepsilon) = 0, \end{cases} \tag{1.1}

其中 $\Delta_h := \partial_1^2 + \partial_2^2$ 为水平Laplace算子， $\varepsilon \in (0,1)$ 是趋于零的垂直粘性系数。该模型源于大尺度海洋动力学等地球物理流体，其典型参数尺度显示：垂直粘性系数约为 $$1$$ – $10^3\, \text{cm}^2/\text{s}$ ，而水平粘性达 $$10^5$$ – $10^8\, \text{cm}^2/\text{s}$ ，即存在强各向异性。因此，自然提出核心问题：当 $\varepsilon \to 0$ 时，解 $(\rho^\varepsilon, u^\varepsilon)$ 是否收敛于水平耗散极限系统

\begin{cases} \partial_t \rho^0 + \operatorname{div}(\rho^0 u^0) = 0, \\ \rho^0 \partial_t u^0 + \rho^0 u^0 \cdot \nabla u^0 - \Delta_h u^0 - \nabla \operatorname{div} u^0 + \nabla P(\rho^0) = 0, \end{cases} \tag{1.5'}

（对应文中方程(1.5)在扰动变量下的形式）？且该收敛能否在全局时间区间 $[0,\infty)$ 上严格建立？

关键难点在于：(i) 边界条件失配——原系统(1.1)满足Navier滑移条件 $u^\varepsilon_3|_{x_3=0} = 0$ , $\partial_3 u^\varepsilon_h|_{x_3=0} = 0$ ，而极限系统(1.5’)仅要求 $u^0_3|_{x_3=0} = 0$ ，缺失对 $\partial_3 u^0_h$ 的约束，导致标准Sobolev范数无法控制垂直导数的一致性；(ii) 密度方程无垂直耗散项， $\partial_3^2 \rho^\varepsilon$ 仅以 $\varepsilon \|\partial_3^2 \rho^\varepsilon\|_{L^2}$ 形式出现在能量泛函中（见Remark 1.4），远弱于各向同性情形下 $\|\partial_3^2 \rho^\varepsilon\|_{L^2}$ 的控制，致使 $\|\partial_3 \rho^\varepsilon\|_{L^\infty}$ 必须作为独立项纳入能量；(iii) 可压性引入密度-速度强耦合，压力项 $\nabla P(\rho^\varepsilon)$ 将速度耗散部分转移至密度，但该机制在垂直方向严重退化。已有工作仅覆盖局部时间或不可压情形：Wang-Xin-Yong [46] 在有界域中对各向同性粘性极限建立了局部时间一致估计；Iftimie-Planas [22] 和 Gao-Peng-Wu-Yao [13] 分别在不可压半空间中解决了滑移边界下的局部/全局各向异性极限问题。然而，可压情形下的全局时间各向异性极限仍是空白。本文旨在填补这一根本性空白，其核心挑战是：在 $\varepsilon \to 0$ 过程中，同时克服边界层效应、垂直耗散缺失、非线性耦合放大三重障碍，建立不依赖于 $\varepsilon$ 的全局一致先验估计，并由此导出强紧性以完成极限过程。

2. 主要定理或主要结果

本文在小初值扰动框架下建立两项严格结果：

定理1.1（全局一致正则性）：设整数 $m \ge 5$ ，常数 $\zeta \in (0, \frac{1}{34}]$ ，初值 $(\varrho_0, u_0)$ 满足小性条件 $E(0) \le \delta_0$ （ $\delta_0 > 0$ 充分小）。则系统(1.3)–(1.4)存在唯一全局解 $(\varrho^\varepsilon, u^\varepsilon)$ ，且对所有 $t \ge 0$ 与 $\varepsilon \in (0,1)$ 有统一估计

E(t) + \int_0^t D(\tau) \, d\tau \le C \delta_0, \tag{1.8}

其中能量泛函 $$E(t)$$ 与耗散泛函 $$D(t)$$ 定义为

\begin{aligned} E(t) &:= \| (\varrho^\varepsilon, u^\varepsilon)(t) \|_{H^m_{\mathrm{co}}}^2 + \| (\partial_3 \varrho^\varepsilon, \partial_3 u^\varepsilon)(t) \|_{H^{m-1}_{\mathrm{co}}}^2 \\ &\quad + \| w^\varepsilon_h(t) \|_{L^\infty}^3 + \| \partial_3 \varrho^\varepsilon(t) \|_{L^\infty}^3 \\ &\quad + \| \Lambda_h^{-(1-\zeta)} (\varrho^\varepsilon, u^\varepsilon)(t) \|_{L^2}^2 + \| \Lambda_h^{-(1-\zeta)} \partial_3 (\varrho^\varepsilon, u^\varepsilon)(t) \|_{L^2}^2 \\ &\quad + \| \partial_t (u^\varepsilon, w^\varepsilon_h, \partial_3 u^\varepsilon_3)(t) \|_{H^2_{\tan}}^2 + \| Z_3 \partial_t (w^\varepsilon_h, \partial_3 u^\varepsilon_3)(t) \|_{L^2}^2 + \varepsilon \| \partial_3^2 \varrho^\varepsilon(t) \|_{H^2_{\mathrm{co}}}^2, \end{aligned}

\begin{aligned} D(t) &:= \| \nabla_h u^\varepsilon(t) \|_{H^m_{\mathrm{co}}}^2 + \| \operatorname{div} u^\varepsilon(t) \|_{H^m_{\mathrm{co}}}^2 + \varepsilon \| \partial_3 u^\varepsilon(t) \|_{H^m_{\mathrm{co}}}^2 \\ &\quad + \| \partial_3 \nabla_h u^\varepsilon(t) \|_{H^{m-1}_{\mathrm{co}}}^2 + \| \partial_3 \operatorname{div} u^\varepsilon(t) \|_{H^{m-1}_{\mathrm{co}}}^2 + \varepsilon \| \partial_3^2 u^\varepsilon(t) \|_{H^{m-1}_{\mathrm{co}}}^2 \\ &\quad + \| \partial_t (\nabla_h u^\varepsilon, \operatorname{div} u^\varepsilon, \nabla_h w^\varepsilon_h, \nabla \partial_3 u^\varepsilon_3)(t) \|_{H^2_{\tan}}^2 + \| Z_3 \partial_t (\nabla_h w^\varepsilon_h, \nabla \partial_3 u^\varepsilon_3)(t) \|_{L^2}^2 \\ &\quad + \varepsilon \big( \| \partial_3 \partial_t u^\varepsilon(t) \|_{H^2_{\tan}}^2 + \| \partial_3 \partial_t w^\varepsilon_h(t) \|_{H^2_{\tan}}^2 + \| \partial_3 Z_3 \partial_t w^\varepsilon_h(t) \|_{L^2}^2 \big) \\ &\quad + \| \nabla \varrho^\varepsilon(t) \|_{H^{m-1}_{\mathrm{co}}}^2 + \varepsilon \| \partial_3^2 \varrho^\varepsilon(t) \|_{H^2_{\mathrm{co}}}^2. \end{aligned}

此外，解还满足两层时间衰减估计：

(1 + t)^{1-\zeta} E_{\tan}(t) + (1 + t)^{2-2\zeta} \overline{E}_{\tan}(t) + \int_0^t (1 + \tau)^{1-2\zeta} D_{\tan}(\tau) \, d\tau + \int_0^t (1 + \tau)^{2-2\zeta} \overline{D}_{\tan}(\tau) \, d\tau \le C \delta_0, \tag{1.9}

其中 $E_{\tan}(t) := \| (\varrho^\varepsilon, u^\varepsilon)(t) \|_{H^m_{\tan}}^2 + \| \partial_3 (\varrho^\varepsilon, u^\varepsilon)(t) \|_{H^{m-1}_{\tan}}^2$ ， $\overline{E}_{\tan}(t) := \| \nabla_h (\varrho^\varepsilon, u^\varepsilon)(t) \|_{H^{m-2}_{\tan}}^2 + \| \nabla_h \partial_3 (\varrho^\varepsilon, u^\varepsilon)(t) \|_{H^{m-3}_{\tan}}^2$ ，且 $$C > 0$$ 独立于 $$t$$ 与 $\varepsilon$ 。

定理1.2（全局垂直粘性消失极限）：在定理1.1条件下，解序列 $(\varrho^\varepsilon, u^\varepsilon)$ 存在子列 $(\varrho^{\varepsilon_n}, u^{\varepsilon_n})$ （ $\varepsilon_n \to 0$ ）使得对任意 $$T > 0$$ ，

(\varrho^{\varepsilon_n}, u^{\varepsilon_n}) \to (\varrho^0, u^0) \quad \text{在 } C\big([0,T]; H^{m-1}_{\mathrm{co},\mathrm{loc}}\big) \text{ 中},

其中极限 $(\varrho^0, u^0)$ 是水平耗散系统(1.5)–(1.6)的全局解。注：由于极限系统(1.5)的全局唯一性尚未解决（Remark 1.5），此处仅能保证子列收敛。

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

证明的核心障碍在于：(i) 边界层导致垂直导数无法用标准Sobolev范数控制；(ii) 密度方程缺乏垂直耗散，迫使 $\|\partial_3 \varrho^\varepsilon\|_{L^\infty}$ 成为关键控制量，但其最大值原理需 $$L^1$$ -in-time 可积性；(iii) 各向异性耗散结构使传统能量方法失效。作者的真正创新点是构建了一套自洽的“两层衰减驱动的 $L^\infty$ 闭合”分析框架，具体体现为三个关键技术突破：

第一，新型混合能量泛函 $$E(t)$$ 的构造：为克服边界失配，采用余法向Sobolev空间 $H^m_{\mathrm{co}}$ ，并显式包含 $\|w^\varepsilon_h\|_{L^\infty}^3$ 与 $\|\partial_3 \varrho^\varepsilon\|_{L^\infty}^3$ 项（而非常规的平方项）。此立方幂次源于最大值原理应用中线性项 $\sim (1+\varrho^\varepsilon)\partial_t u^\varepsilon_3$ 、 $\Delta_h u^\varepsilon_3$ 、 $\varepsilon \partial_3 \nabla_h \cdot u^\varepsilon_h$ 的吸收需求（Remark 1.3, (2.5)式）；若仅为平方，则无法压制线性贡献。同时，泛函中嵌入负阶水平Fourier乘子范数 $\|\Lambda_h^{-(1-\zeta)} (\varrho^\varepsilon, u^\varepsilon)\|_{L^2}^2$ ，为后续衰减估计提供插值支点。

第二，“两层时间衰减”机制的设计与实现：定义切向能量 $E_{\tan}(t)$ 与高阶切向能量 $\overline{E}_{\tan}(t)$ ，通过插值不等式 $\|f\|_{L^2} \lesssim \|\Lambda_h^{-(1-\zeta)} f\|_{L^2}^{\theta} \|\nabla_h f\|_{L^2}^{1-\theta}$ （ $\theta = 1-\zeta$ ）将负阶范数与耗散 $D_{\tan}(t)$ 关联，导出非线性微分不等式

\frac{d}{dt} \widetilde{E}_{\tan}(t) + \kappa C_0^{-\frac{1}{1-\zeta}} \widetilde{E}_{\tan}(t)^{1 + \frac{1}{1-\zeta}} \le 0,

从而得到首层衰减 $\widetilde{E}_{\tan}(t) \lesssim C_0 (1+t)^{-(1-\zeta)}$ （(2.7)式）。再以此为基础，在更高阶能量上重复论证，获得第二层更快衰减 $\overline{E}_{\tan}(t) \lesssim C_0 (1+t)^{-(2-2\zeta)}$ （(2.9)式）。该两层结构的关键在于：首层衰减权重 $(1+\tau)^{1-2\zeta}$ 与第二层权重 $(1+\tau)^{2-2\zeta}$ 恰好保证了 $$L^1$$ -in-time 积分 $\int_0^t \|\partial_3 \varrho^\varepsilon\|_{L^\infty} \|\partial_3^2 u^\varepsilon_h\|_{L^\infty} \, d\tau$ 与 $\int_0^t \|\nabla_h \varrho^\varepsilon\|_{L^\infty} \|\partial_3 \operatorname{div} u^\varepsilon\|_{L^\infty} \, d\tau$ 的收敛性（(2.2)–(2.4)式），这是闭合 $L^\infty$ 估计的充要条件。

第三，垂直导数的“方程提取”策略：为避免直接估计 $\|\partial_3^2 u^\varepsilon_h\|_{L^\infty}$ （会引入 $\varepsilon^{-1}$ 奇异性），作者从水平动量方程中直接提取 $\varepsilon \|\partial_3^3 u^\varepsilon_h\|_{H^2_{\tan}}$ 的控制：

\varepsilon \|\partial_3^3 u^\varepsilon_h\|_{H^2_{\tan}} \lesssim \left\| \partial_3 \left( (1+\varrho^\varepsilon) \partial_t u^\varepsilon_h + (1+\varrho^\varepsilon) u^\varepsilon \cdot \nabla u^\varepsilon_h - \Delta_h u^\varepsilon_h - \nabla_h \operatorname{div} u^\varepsilon + (1+\varrho^\varepsilon)^2 \nabla_h \varrho^\varepsilon \right) \right\|_{H^2_{\tan}},

右侧由 $E(t)^{1/2} + D(t)^{1/2}$ 控制。类似地， $\partial_3 \operatorname{div} u^\varepsilon$ 由垂直动量方程(2.3)显式解出。这些“方程提取”步骤将高阶垂直导数的估计转化为对已知能量项的控制，是绕过边界层奇异性、实现一致性的核心技术。

4. 与之前工作的比较

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

时间尺度与适用范围：此前所有可压情形的结果均为局部时间。Wang-Xin-Yong [46] 在有界域中对各向同性粘性极限 $\varepsilon (\Delta_h + \partial_3^2 + \nabla \operatorname{div})$ 建立了局部一致估计，但其系统(1.10)含完整拉普拉斯项，而本文系统(1.1)的耗散仅为 $\Delta_h + \varepsilon \partial_3^2 + \nabla \operatorname{div}$ ，后者在垂直方向仅具 $\varepsilon$ -阶耗散（Remark 1.4）。Masmoudi-Sun-Wang-Zhang [34] 考虑了Dirichlet边界下不可压与垂直粘性消失的联合极限，但同样限于局部时间。本文首次将可压各向异性极限推进至全局时间 $[0,\infty)$ ，且适用于无限域 $\mathbb{R}^3_+$ ，这是质的飞跃。

假设条件与代价：相比不可压类比工作 [13]，本文必须处理密度扰动 $\varrho^\varepsilon$ 这一新增未知量。其无垂直耗散，且与速度通过压力强耦合，导致能量泛函中必须显式包含 $\|\partial_3 \varrho^\varepsilon\|_{L^\infty}^3$ 项（而 [13] 无需此类项）。此外，[13] 的衰减率可达 $(1+t)^{-1}$ ，而本文因密度耦合的弱耗散，衰减速率受限于 $\zeta \le \frac{1}{34}$ ，最优为 $(1+t)^{-1+\frac{1}{34}}$ （Remark 1.1）。此速率限制是处理可压性的必然代价，亦反映在插值不等式(3.22)的精细计算中。

中文标题: 含分数阶拉普拉斯算子的非线性Hartree方程解的渐近行为

摘要翻译: 本文研究如下非局部问题\begin{equation*}\left\lbrace \begin{aligned} &A_{s} u=(|x|^{-(n-2s)}\ast u^{2_{s}^{\sharp}-1-ε})u^{2_{s}^{\sharp}-2-ε} \quad\quad\hspace{3.5mm} \mbox{in}\hspace{2mm}Ω,\ &u>0\quad\quad \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\hspace{2mm}\mbox{in}\hspace{2mm}Ω,\ &u=0\quad \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\hspace{2mm}\mbox{on}\hspace{2mm}\mathbb{R}^n\setminusΩ, \end{aligned} \right.\end{equation*}其中 $$Ω$$ 是 $\mathbb{R}^n$ 中光滑有界区域， $$0<s<1$$ ， $n\in(2s,\min\{6s,n+2s\})$ ， $$ε>0$$ 充分小， $2_{s}^{\sharp}-1=(n+2s)/(n-2s)$ ，且 $A_{s}$ 表示定义在 $$Ω$$ 上、具外部零Dirichlet边界条件的分数阶拉普拉斯算子 $(-Δ)^{s}$ 。上述问题可约化为次临界分数阶系统

中文标题: 非线性Hartree型方程正驻波的对称性与分类

摘要翻译: 本文研究强耦合系统

\begin{cases} \displaystyle - \Delta u + \tau u = \frac{2 p}{p + q} \left( I_\alpha\ast |v|^q \right) |u|^{p - 2} u &\text{in} ~ \mathbb{R}^N, \\ \\ \displaystyle - \Delta v + \eta v = \frac{2 q}{p + q} \left( I_\alpha\ast |u|^p \right) |v|^{q - 2} v &\text{in} ~ \mathbb{R}^N, \end{cases}

的正解的一些定性性质，其中 $\tau, \eta > 0$ , $N \in \mathbb{N}$ , $0 < \alpha < N$ ,

\max \left\{1, \frac{2 \alpha}{N}\right\} < p, q < 2^* \quad \text{and} \quad \frac{2 (N + \alpha)}{N} < p + q < 2_\alpha^*,

此处 $I_\alpha$ 表示Riesz位势，

2^* := \begin{cases} \infty, &\text{if} ~ N \in \{1, 2\}, \\ \frac{2 N}{N - 2}, &\text{if} ~ N \geq 3, \end{cases} \quad \text{and} \quad 2_\alpha^* := \begin{cases} \infty, &\text{if} ~ N \in \{1, 2\}, \\ \frac{2 (N + \alpha)}{N - 2}, &\text{if} ~ N \geq 3. \end{cases}

更精确地，利用移动平面法，我们证明：当 $p, q \geq 2$ 时，该系统的正解是径向对称且严格径向递减的；并在 $$p = q$$ 且 $\tau = \eta$ 的情形下，得到了正基态解的分类结果。

10. ZK方程在半带域上的低正则性适定性

Low-regularity well-posedness for the ZK equation on a half-strip

作者: E Avelino, G Doronin

类别: math.AP

发布日期: 2026-05-21

arXiv ID: 2605.23092v1

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

摘要翻译

中文标题: ZK方程在半带域上的低正则性适定性

摘要翻译: 本文研究定义在半带域上、带有线性输运项且具非齐次边界条件的Zakharov–Kuznetsov方程。借助适配ZK色散结构的Bourgain型空间、各向异性光滑性估计以及边界迹估计，我们在 $$L^2$$ 中建立了该方程的局部适定性。

11. 关于一维三次非线性薛定谔方程的终态问题

On the final-state problem for the 1D cubic NLS

作者: Gong Chen, Yongyu Qiang

类别: math.AP

发布日期: 2026-05-21

arXiv ID: 2605.23063v1

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

摘要翻译

中文标题: 关于一维三次非线性薛定谔方程的终态问题

摘要翻译: 我们考虑一维三次非线性薛定谔方程

\ii\partial_tu+\frac12\partial_{xx}u=\la|u|^2u,\,λ=\pm 1

，并针对小渐近数据求解其终态问题（即修正波算子问题）。更精确地，给定一个小的 $$W(ξ)$$ ，我们构造一个解 $$u$$ ，使得 \begin{equation*} u\rightarrow (2π)^{-1/2}(\ii t)^{-1/2}e^{\ii x^2/(2t)}, W!\Big(\frac{x}{t}\Big)\exp(-\ii\la|W(x/t)|^2\log t). \end{equation*} 关键在于，我们设计了一个压缩映射，从而可在 Kato–Pusateri \cite{KP} 的精神下对 $$w$$ 进行分析，其中强迫项仅依赖于终态数据 $$W$$ 。该方案易于推广至其他终态问题，前提是相应的正向问题具有完备理论。

不相关论文（基本信息）

1. 可逆化学反应扩散系统存在光滑解的简单证明

Simple proofs for the existence of smooth solutions to a reaction-diffusion system modeling reversible chemistry

作者: Hector Bouton, Laurent Desvillettes, Helge Dietert

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23709v1

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

摘要: We present in this work a very short proof for the existence, uniqueness and smoothness in dimensions $d\leq 3$ of the system of reaction diffusion $\partial\_t a\_i - d\_i Δa\_i = (-1)^i (a\_1 a\_3 - a\_2 a\_4)$ , where $a\_i \geq 0$ model the concentrations of chemical species undergoing a chemical reaction and diffusing (each with its diffusion rate $d\_i > 0$ ) in a bounded container.

2. 拟线性椭圆方程在奇点集附近的逐点估计

Pointwise Estimates Near Singular Sets for Quasilinear Elliptic Equations

作者: Juan Pablo Alcon Apaza

类别: math.AP, math.FA

发布日期: 2026-05-22

arXiv ID: 2605.23835v1

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

摘要: In this work, we study the removability of boundary singular sets for certain classes of quasilinear elliptic equations in domains $$Ω$$ of an $$n$$ -dimensional Finsler manifold ( $\mathcal{M}, F, \vartheta$ ). We work with Lipschitz functions $$ρ_1$$ and $$ρ_2$$ satisfying distance-type properties; in particular, $F(\cdot, \boldsymbol{\nabla} ρ_1) \leq 1$ and $F(\cdot, \boldsymbol{\nabla} ρ_2) \leq 1$ a.e. in $\mathcal{M}$ . The singular set is defined by $Γ=ρ_1^{-1}(\{0\})$ . The model problem is $-Δ_{p(x)} u+|u|^{q-1} u=0$ in domains of $\mathbb{R}^n \cong \mathbb{R}^d \times \mathbb{R}^{n-d} \cong ρ_1^{-1}(\{0\}) \times ρ_2^{-1}(\{0\})$ , where $ρ_1(x)=|(x_{d+1}, \ldots, x_n)|$ and $ρ_2(x)=|(x_1, \ldots, x_d)|$ . The main tool in our analysis is the estimate

|u(x)| \leq \mathbf{C} ρ_1(x)^{-τ}

near $$Γ$$ for weak solutions $u \in W_{loc}^{1, p(x)}(\barΩ \backslash(Γ\cup Σ) ; \vartheta) \cap L_{loc}^{\infty}(\barΩ \backslash(Γ\cup Σ))$ , where the constants $\mathbf{C}>0$ and $$τ>0$$ converge to positive values as $p^{+} \rightarrow 1$ . This estimate is a key ingredient in proving that the singularity at $$Γ$$ is removable. Moreover, in a bounded domain $$Ω$$ , using this estimate and assuming that, for every variable exponent satisfying $1<p^{-} \leq p^{+}<\min \{2, q+1\}$ , there exists a weak solution $u_p \in W_{loc}^{1, p(x)}(Ω; \vartheta) \cap L_{loc}^{\infty}(Ω)$ of

-\operatorname{div}\left(|\boldsymbol{\nabla} u_p|_F^{p-2} \boldsymbol{\nabla} u_p\right)+|u_p|^{q-1} u_p=0 \quad \text { in } Ω,

we prove that, for every $U \Subset Ω$ , there exists a subsequence $\{u_{p_m}\}$ , with $p_m^{+} \rightarrow 1$ , that converges to a solution $u \in B V(U ; \vartheta) \cap L^{q+1}(U ; \vartheta)$ of

-Δ_1 u+|u|^{q-1} u=0 \quad \text { in } U .

3. Brenier 映射的整体估计

Global estimates on the Brenier map

作者: Andrea Bidoia

类别: math.AP, math.PR

发布日期: 2026-05-22

arXiv ID: 2605.23731v1

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

摘要: Caffarelli’s contraction theorem and the analogous Laplacian result in [arXiv:2411.12109, arXiv:2501.11382] are two examples of how log-Hessian bounds on probability densities yield estimates on the derivative of the corresponding Brenier map with optimal dimensional dependence. The main goal of this paper is to extend such phenomenon to a broader class of convex estimates such as norms.

4. Z峰重味物理项目中的超精细粒度量热器性能

Ultra-Granular Calorimeter Performances for the Heavy Flavor Physics Program at the Z Peak

作者: J. -C Brient

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23667v1

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

摘要: Various decays of the B mesons are here used to establish the performances of an ultra-granular electromagnetic calorimeter for heavy flavour physics at an electron positron accelerator running at the Z peak. The silicon-tungsten electromagnetic calorimeter of the ILD concept is used for this purpose, enhanced by a timing capability. When possible, a $$π$$ 0 mass fit of $$γ$$ $$γ$$ system is performed to improve the $$π$$ 0 energy resolution. It is also shown that in the presence in the final state of a photon without a $$π$$ 0 , the identification of genuine photon(s) versus fake photon(s) coming from K L ’s, neutron’s, debris of hadronic shower or other high energy $$π$$ 0 , is essential. It allows for a good precision and a good signal over background ratio for this kind of physics. The possible impact for the tau physics is discussed.

5. 具有不同 p-增长非齐次算子的非线性透射特征值问题

Nonlinear Transmission Eigenvalue Problems with Nonhomogeneous Operators of Different p-Growth

作者: Luminita Barbu, Raluca-Gabriela Turtoi

类别: math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23266v1

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

摘要: Let $Ω\subset \mathbb{R}^N$ , $N \ge 2$ , be a bounded domain with Lipschitz boundary, divided by a Lipschitz hypersurface $$Σ$$ into two open, disjoint Lipschitz subdomains $$Ω_1$$ and $$Ω_2$$ . We study a nonlinear transmission eigenvalue problem driven by nonhomogeneous operators with $$p_i$$ - growth in each subdomain $$Ω_i$$ , $$i=1,2$$ , and subject to continuity and flux transmission conditions across the interface $$Σ$$ . The real parameter $$λ$$ appears both in the equations and in the nonlinear boundary conditions. Using variational methods, we prove the existence of an unbounded sequence of eigenvalues. Under additional assumptions, we establish that the set of eigenvalues coincides with the entire interval $(0,\infty)$ . As a particular case, we obtain the corresponding eigenvalue results for the associated single-domain problem.

6. 二维声学超材料中亚波长共振的 Fourier–Galerkin 方法

Fourier–Galerkin Methods for Subwavelength Resonances in 2D Acoustic Metamaterials

作者: Jinghao Cao

类别: math.NA, math.AP

发布日期: 2026-05-22

arXiv ID: 2605.23251v1

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

摘要: We present a Fourier–Galerkin framework for the analysis and computation of subwavelength resonances in two-dimensional scattering problems in finite domains. Starting from the boundary integral formulation, we project the operator onto Fourier modes and derive an explicit finite-dimensional effective matrix whose singularity characterizes the resonant frequencies. In the subwavelength regime, we obtain asymptotic expansions of this matrix in terms of $$ω$$ and the material contrast, identifying the leading-order operators and their kernel structure. This reduction transforms the resonance problem into a low-dimensional nonlinear eigenvalue problem, avoiding large-scale discretizations and global root-search procedures. The entries of the effective matrix are explicitly computable and admit fast evaluation using FFT-based quadrature. The resulting approach provides an efficient and robust computational framework for resonances in general smooth geometries.

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

1. 利用 Lambert W 函数求解 Navier-Stokes 方程的一个简单情形#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. Surface Quasi-Geostrophic 方程中移动刚性区域的局部适定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 单速度Baer-Nunziato模型的压力与温度松弛极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 一类具有超二次非线性的非局部问题的山路算法#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 半空间中可压 Navier–Stokes 方程的整体一致正则性与垂直粘性消失极限#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. KdV 和 KP 方程在负阶 Sobolev 空间中的范数膨胀#

摘要翻译#

2. 具有二维节点集乘积结构的椭圆方程解之比的先验正则性估计#

摘要翻译#

3. 带磁漂移的双线性薛定谔方程的小时间近似可控性#

摘要翻译#

4. 含分数阶拉普拉斯算子的非线性Hartree方程解的渐近行为#

摘要翻译#

5. Hamilton-Jacobi 方程在 Neumann 边界条件下具有分数模的最优半凹性#

摘要翻译#

6. 具有非线性边界相互作用的薛定谔方程的 Strichartz 估计#

摘要翻译#

7. 关于旋度-测度场的切向迹#

摘要翻译#

8. 广义Ostrovsky方程的空间衰减与非线性光滑性#

摘要翻译#

9. 非线性Hartree型方程正驻波的对称性与分类#

摘要翻译#

10. ZK方程在半带域上的低正则性适定性#

摘要翻译#

11. 关于一维三次非线性薛定谔方程的终态问题#

摘要翻译#

不相关论文（基本信息）#

1. 可逆化学反应扩散系统存在光滑解的简单证明#

2. 拟线性椭圆方程在奇点集附近的逐点估计#

3. Brenier 映射的整体估计#

4. Z峰重味物理项目中的超精细粒度量热器性能#

5. 具有不同 p-增长非齐次算子的非线性透射特征值问题#

6. 二维声学超材料中亚波长共振的 Fourier–Galerkin 方法#

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

1. 利用 Lambert W 函数求解 Navier-Stokes 方程的一个简单情形

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. Surface Quasi-Geostrophic 方程中移动刚性区域的局部适定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 单速度Baer-Nunziato模型的压力与温度松弛极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 一类具有超二次非线性的非局部问题的山路算法

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 半空间中可压 Navier–Stokes 方程的整体一致正则性与垂直粘性消失极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. KdV 和 KP 方程在负阶 Sobolev 空间中的范数膨胀

摘要翻译

2. 具有二维节点集乘积结构的椭圆方程解之比的先验正则性估计

摘要翻译

3. 带磁漂移的双线性薛定谔方程的小时间近似可控性

摘要翻译

4. 含分数阶拉普拉斯算子的非线性Hartree方程解的渐近行为

摘要翻译

5. Hamilton-Jacobi 方程在 Neumann 边界条件下具有分数模的最优半凹性

摘要翻译

6. 具有非线性边界相互作用的薛定谔方程的 Strichartz 估计

摘要翻译

7. 关于旋度-测度场的切向迹

摘要翻译

8. 广义Ostrovsky方程的空间衰减与非线性光滑性

摘要翻译

9. 非线性Hartree型方程正驻波的对称性与分类

摘要翻译

10. ZK方程在半带域上的低正则性适定性

摘要翻译

11. 关于一维三次非线性薛定谔方程的终态问题

摘要翻译

不相关论文（基本信息）

1. 可逆化学反应扩散系统存在光滑解的简单证明

2. 拟线性椭圆方程在奇点集附近的逐点估计

3. Brenier 映射的整体估计

4. Z峰重味物理项目中的超精细粒度量热器性能

5. 具有不同 p-增长非齐次算子的非线性透射特征值问题

6. 二维声学超材料中亚波长共振的 Fourier–Galerkin 方法