生成时间: 2025年11月10日 04:03:02

重点关注论文数量: 8

了解领域论文数量: 1

不相关论文数量: 29

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

1. Classical Euler flows generate the strong Guderley imploding shock wave

经典欧拉流产生强古德雷利内爆激波

作者: Giorgio Cialdea, Steve Shkoller, Vlad Vicol

类别: math.AP, physics.flu-dyn, 35L67, 76L05

发布日期: 2025-10-22

arXiv ID: 2510.19688v2

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

详细分析

中文标题: 经典欧拉流产生强古德雷利内爆激波

研究对象和背景

本文研究的是可压缩欧拉方程（compressible Euler equations）在径向对称情形下的解，特别是从光滑、无激波的初始数据出发，如何演化出一个著名的自相似内爆激波——古德雷利激波（Guderley imploding shock wave）。这一现象在高能量密度物理、惯性约束聚变以及天体物理中的超新星爆发等场景中具有重要应用。

所研究的系统是 $d = 2$ 或 $3$ 维空间上的理想气体模型下的欧拉方程，其守恒律形式如下：

其中 $\rho > 0$ 是密度，$u : \mathbb{R}^d \times \mathbb{R} \to \mathbb{R}^d$ 是速度场，$E$ 是总能量，$p$ 是压力。系统的封闭通过理想气体状态方程实现：

其中 $\gamma > 1$ 是绝热指数（adiabatic exponent），文中考虑 $\gamma \in (1,3]$ 的范围。

为了便于分析非等熵流动，作者引入了伪熵平方根变量 $b = e^{S/2}$，其中 $S$ 是比熵。由此定义重缩放声速 $\sigma = \rho^\alpha b$，其中 $\alpha = \frac{\gamma - 1}{2}$。利用这些变量，原欧拉方程可转化为对称形式：

在径向对称假设下，令 $r = |x|$，并设 $u = u(r,t)e_r$，$\rho = \rho(r,t)$，$b = b(r,t)$，$\sigma = \sigma(r,t)$，则上述系统简化为一维形式：

进一步引入黎曼变量：$z = u - \sigma$（快声学变量），$w = u + \sigma$（慢声学变量），并结合微分黎曼变量（Differentiated Riemann Variables, DRV）：

得到一组关于梯度演化的偏微分方程组（见公式 (1.9a)-(1.9c)），这对方程奇异性分析至关重要。

本文关注的核心现象是古德雷利自相似内爆激波。该解描述了一个强激波从远处向原点对称收缩，并在有限时间 $t=0$ 处发生坍缩的过程。其自相似结构由以下变量定义：

其中 $\lambda > 1$ 是必须通过求解非线性特征值问题确定的自相似指数。原始变量具有如下自相似试探解：

函数 $U(\xi), C(\xi), R(\xi)$ 满足一组常微分方程系统（ODEs），并在 $\xi = 1$ 处满足 Rankine-Hugoniot 跳跃条件，在某个临界点 $\xi_s > 1$ 处满足正则性条件（分子与分母同时为零），从而唯一确定 $\lambda$ 和解轮廓。

尽管 Guderley 解早在1942年就被提出，并广泛用于数值模拟和物理建模中，但长期以来缺乏严格的 PDE 理论支持。直到最近 Jang-Liu-Schrecker [18] 才首次严格证明了这类自相似激波解的存在性和解析性。

本文的研究目标正是填补另一个长期空白：是否存在一组经典、光滑、无激波的初始数据，使得对应的唯一弱熵解能够自然演化出这个强自相似激波？即验证 Guderley 解是否真的属于某个“吸引子”的吸引域。

为此，作者构造的解定义在时空区域 $\mathbb{R}^d \times (T_{\mathrm{in}}, T_{\mathrm{fin}}]$ 上，其中 $T_{\mathrm{in}} < T_* < T_{\mathrm{fin}} < 0$。初始时刻 $t = T_{\mathrm{in}} < 0$ 给定的数据位于函数空间 $C^{1,1/3}(\mathbb{R}_+)$ 中，意味着一阶导数 Hölder 连续，指数为 $1/3$。特别地，存在一个“静止核”区域 $r < (-T_{\mathrm{fin}})^{1/\lambda}$，其中流体处于静止状态 $(u,\rho,c)=(0,1,0)$。边界条件隐含在自相似匹配中：在 $t = T_{\mathrm{fin}}$ 时，解与 Guderley 剖面完全吻合，包括激波位置、强度及所有阶的导数。

主要定理或主要结果

文章证明的主要定理是 Theorem 1.3，其陈述如下：

定理 1.3（主定理）
固定空间维度 $d \in \{2,3\}$，任意给定最终匹配时间 $T_{\mathrm{fin}} < 0$，则存在初始时间 $T_{\mathrm{in}} < T_*$ 和预激波形成时间 $T_* \in (T_{\mathrm{in}}, T_{\mathrm{fin}})$，以及对应的时间 $t = T_{\mathrm{in}}$ 上的径向初始数据集 $(w_{\mathrm{in}}, z_{\mathrm{in}}, b_{\mathrm{in}})$，满足以下性质：

• 正则性：数据属于空间 $C^{1,1/3}(\mathbb{R}_+)$。
• 静止核：在核心区 $r < (-T_{\mathrm{fin}})^{1/\lambda}$ 内，$(u_{\mathrm{in}}, \rho_{\mathrm{in}}, b_{\mathrm{in}}) \equiv (0,1,0)$。

该初始数据生成唯一的径向欧拉解，其演化分为五个阶段：

1. 光滑演化：当 $t \in [T_{\mathrm{in}}, T_*)$ 时，解 $(w,z,b)$ 光滑演化，保持在 $C^{1,1/3}$ 类中。
2. 预激波形成：在精确时刻 $t = T_*$，于某空间点 $r_* > 0$ 处形成预激波奇点。主导黎曼变量 $z$ 出现 $C^{1/3}$ 尖点，局部渐近行为为：
   $$ z(r,T_*) = z(r_*,T_*) + a(r - r_*)^{1/3} + b(r - r_*)^{2/3} + O(|r - r_*|), $$
   而 $w$ 和 $b$ 保持更高正则性。
3. 激波发展：当 $t \in (T_*, T_{\mathrm{fin}}]$ 时，从预激波点 $(r_*, T_*)$ 发展出不连续激波，沿轨迹 $\{r = s(t)\}$ 传播，且 $s(T_*) = r_*$。流动成为唯一的“正规激波解”（regular shock solution），在激波两侧无限光滑。
4. 与古德雷利解匹配：在时间 $t = T_{\mathrm{fin}}$，解 $(\rho,u,E)$（满足 Lax 熵条件）和激波前沿 $s(T_{\mathrm{fin}})$ 完全匹配经典的古德雷利自相似内爆激波解。此后所有时间 $t \in [T_{\mathrm{fin}}, 0)$，演化与古德雷利解一致。
5. 最终内爆：在 $t = 0$ 时，演化以原点 $r = 0$ 处的速度 $u \to -\infty$ 和能量 $E \to +\infty$ 结束，而密度 $\rho$ 保持有界。

作为推论，结合已知的反射古德雷利解全局存在性 [17]，得到：

推论 1.4（全局时间欧拉演化）
存在初始时间 $T_{\mathrm{in}} < 0$ 和属于 $C^{1,1/3}(\mathbb{R}_+)$ 的径向初始数据 $(u_{\mathrm{in}}, \rho_{\mathrm{in}}, E_{\mathrm{in}})$，使得对应的唯一解 $(u,\rho,E)$ 在区间 $[T_{\mathrm{in}}, \infty)$ 上全局存在，并经历完整物理序列：光滑演化 → 激波形成 → 强古德雷利内爆（$t=0$）→ 反射为向外传播的自相似爆炸激波（$t>0$）。

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

a) 核心分析学工具

本文的证明策略极为精细，融合了几何 PDE 方法、动力系统思想与迭代构造技术。其核心工具包括：

• 反向时间激波发展理论（Backward-in-time shock development）：这是论文最核心的技术创新之一。传统上，激波一旦形成便不可逆；但作者巧妙地将问题“倒放”，从目标状态——即 $t = T_{\mathrm{fin}}$ 处的强古德雷利激波——开始，反向构造一个弱激波，再继续反向追踪至预激波奇点 $T_*$。这种“逆向工程”允许他们控制激波强度的增长路径。

• 不动点定理与迭代方案：整个构造依赖于在适当函数空间中建立收缩映射。例如，在附录 A 和 B 中，作者定义了两个函数空间 $\mathcal{D}$ 和 $\mathcal{B}$，分别用于处理不同阶段的动力学。通过精心设计的迭代格式，证明了解的存在唯一性。

• 特征线方法与渐近分析：由于系统具有双曲结构，作者深入使用特征线追踪变量演化。特别是在预激波附近，主导变量 $z$ 沿 $\lambda_1 = u - \alpha \sigma$ 特征线演化，其梯度 $\mathring{z}$ 满足输运型方程。通过对 $\mathring{z}$ 的精细估计，揭示出 $C^{1/3}$ 奇性的必然出现。

• 高阶匹配与渐近展开：在 $t = T_{\mathrm{fin}}$ 处实现“所有阶匹配”要求极高精度。作者不仅匹配函数值和一阶导数，还确保所有高阶导数连续衔接。这需要对 Guderley 剖面进行泰勒展开，并反向构造匹配的初值。

• 正规激波解框架（Regular shock solution framework）：基于 Christodoulou-Miao [6] 等人的工作，作者采用“正规激波解”概念，保证激波前后流场足够光滑，跳跃满足 Rankine-Hugoniot 条件和 Lax 熵条件。

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

• 从强激波反推至预激波的奇性控制：直接从强激波反向熄灭会导致加速度发散（见第4页：“would force the shock accelerations to become singular”）。为此，作者引入两段式激波曲线：在 $(T_*, T_\circ]$ 使用高度正则的曲线 $r = c(t)$ 实现 $C^\infty$ 过渡到 $C^{1/3}$ 预激波；在 $(T_\circ, T_{\mathrm{fin}}]$ 使用 $r = s(t)$ 连接弱激波与强 Guderley 激波。这种分段构造避免了物理上不合理的不对称奇性。

• 非线性项的处理：DRV 方程 (1.9) 中包含大量非线性耦合项，如 $(w^2 - z^2)/r$、$\rho^\alpha \partial_r b$ 等。作者通过加权能量估计和 Sobolev 嵌入控制这些项的增长，尤其是在 $r \to 0$ 或接近奇点时。

• 静止核的兼容性：核心区 $(u,\rho,c)=(0,0,1)$ 是退化的（声速为零），不能直接嵌入标准双曲理论。作者通过构造外部扰动逐渐衰减至核区的方式来实现平滑连接。

• 唯一性论证：在附录 C 中，作者证明从一个 $C^{1/3}$ 预激波出发，只能发展出唯一的弱激波。这依赖于对特征速度差异的精细控制和 Gronwall 不等式。

与之前工作的比较

a) 突破与改进

本文相对于先前工作的突破体现在多个层面：

• 首次实现从光滑数据到强自相似激波的严格构造：此前虽有 [6,9,10,26,47] 研究从光滑数据形成激波，但仅限于弱激波 regime，即激波强度随 $(t - T_*)^{1/2}$ 缓慢增长。本文成功跨越至强非线性 regime，使激波强度可控地增强直至完全匹配 Guderley 剖面。

• 验证了 Guderley 解的“吸引子”角色：长久以来，Guderley 解被认为是一个“普适吸引子”，但缺乏数学基础。本文首次严格证明其吸引域非空，且包含一大类 $C^{1,1/3}$ 初始数据，从而为中间渐近理论提供了坚实支撑。

• 克服线性不稳定性悖论：Guderley 解本身对非径向扰动是线性不稳定的 [4]，这似乎与其作为吸引子的角色矛盾。本文表明：尽管最终会偏离，但在坍缩前的中间阶段，径向对称解确实会被吸引至 Guderley 剖面，体现了“中间渐近”思想的本质。

• 构造了完整的物理演化路径：这是首个已知的欧拉解，完整经历了 (i) 光滑初值 → (ii) 梯度灾难（预激波）→ (iii) 弱激波 → (iv) 强激波 → (v) 自相似匹配 → (vi) 内爆与反射 的全过程，极大丰富了多维守恒律解的图景。

• 降低对初始数据的要求：相比要求更高正则性的早期工作，本文允许初始数据仅为 $C^{1,1/3}$，更贴近物理实际。

综上所述，本论文不仅解决了数学物理中的一个长期开放问题，而且发展了一套强大的反向构造方法，为未来研究其他自相似波现象（如爆轰波、接触间断演化等）提供了范式。

2. Well-posedness of initial boundary value problems for 2D compressible MHD equations in domains with corners

角点区域中二维可压缩磁流体动力学方程初边值问题的适定性研究

作者: Wen Guo, Ya-Guang Wang

类别: math.AP, math-ph, math.MP, 35L04, 76W05, 35A01

发布日期: 2025-11-07

arXiv ID: 2511.05343v1

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

详细分析

中文标题: 角点区域中二维可压缩磁流体动力学方程初边值问题的适定性研究

研究对象和背景

本文研究的是定义在具有角点（corners）的有界区域上的二维可压缩理想磁流体动力学（compressible ideal magnetohydrodynamic, MHD）方程组的初边值问题（initial boundary value problem, IBVP）的适定性（well-posedness）。这类问题在数学理论和物理应用中均具有重要意义。从物理角度看，该模型描述了等离子体在刚性、完全导电壁内的运动行为，而实际中的聚变装置横截面或金属容器常呈现多边形结构，因此自然引入了带角点的几何区域。从数学角度看，由于角点的存在破坏了边界光滑性，传统的基于光滑向量场（如切向与法向导数）的分析方法不再适用，从而带来了深刻的解析障碍。

所研究的控制方程为如下二维可压缩理想MHD系统：

其中 $ u = (u_1,u_2)^t $ 是速度场，$ b = (b_1,b_2)^t $ 是磁场，$ p $ 是压力，$ s $ 是熵，密度 $ R(p,s) $ 由物态方程给出，且假设 $ Q(p,s) = \partial_p R > 0 $，保证系统的双曲性。此外，要求初始数据满足 $ \nabla \cdot b_0 = 0 $ 和边界条件 $ b \cdot \nu = 0 $，则在整个演化过程中这些约束保持成立，避免了系统的超定性。

初始条件为：

边界条件为不可穿透与完全导电壁条件：

其中 $ \nu $ 是边界上光滑部分的外法向量，$ \{q_n\}_{n=1}^N $ 是边界上的角点。

研究区域 $ \Omega \subset \mathbb{R}^2 $ 是一个有界、单连通的“曲线多边形”（curvilinear polygon），其边界由有限条光滑曲线段连接而成，在连接处形成角点，每个角点 $ q_n $ 处的角度记为 $ \omega_n \in (0,\pi) $，即均为凸角。这种非光滑几何结构导致两个根本困难：

1. 缺乏全局光滑的切向向量场：在角点附近，无法定义在整个边界上非零且光滑的切向向量场，这使得传统用于估计切向导数的各向异性Sobolev空间（如Chen [5]中引入的）无法直接使用。
2. 法向导数耦合严重：在角点邻域内，不同方向的法向导数相互耦合，导致不能像光滑边界情形那样通过方程逐阶估计法向导数。

为克服上述困难，作者引入了一类全新的适应于角点几何的各向异性Sobolev空间 $ H^m_*(\Omega) $。该空间的核心思想是利用构造性的加权切向导数来替代标准的空间导数。具体地，通过单位分解技术构造两个在整个区域 $ \Omega $ 上光滑定义的向量场 $ w_1(x), w_2(x) $，它们在远离角点的边界部分与真实切向一致，在角点附近则分别沿两个边界的法向方向衰减（例如用距离函数乘以旋转后的梯度），从而保证其与局部法向正交。

由此定义空间：

其范数为：

这一构造的关键优势在于：在远离角点处，$ H^m_* $ 退化为经典的各向异性Sobolev空间；在角点邻域（如第一象限 $ x_1>0,x_2>0 $），其范数等价于加权形式 $ \|x_1^{\alpha_1}x_2^{\alpha_2}\partial_1^{\beta_1}\partial_2^{\beta_2}f\|_{L^2} $，权重体现了对边界奇性的补偿。

此外，文章强调了解的相容性条件（compatibility conditions）的重要性：初始数据必须满足直到 $ m-1 $ 阶的边界相容性，以确保时间导数在边界上有良好定义。这也是为何定理中要求初始数据属于 $ H^m_* $ 而非普通 $ H^m $ 的原因——后者不足以控制高阶迹。

综上所述，本工作的核心挑战是在非光滑区域中处理强耦合的双曲型系统，其难点不仅来自MHD本身的复杂结构（速度与磁场的强非线性耦合），更源于几何奇异性对函数空间和能量估计框架的根本限制。

主要定理或主要结果

本文的主要结果由以下定理给出：

定理 1.1（主定理）
设整数 $ m \ge 8 $，区域 $ \Omega $ 满足Notation 1.1中的条件，即为具有有限个凸角的曲线多边形。假设每个角点 $ q_n $ 处的角度 $ \omega_n $ 满足：

其中 $ [\cdot] $ 表示取整函数。若初始数据 $ (u_0,b_0,p_0,s_0) \in H^m_*(\Omega) $ 满足问题 (1.1)-(1.3) 在 $ \partial\Omega \setminus \{q_1,\dots,q_N\} $ 上直到 $ m-1 $ 阶的相容性条件，则存在 $ T > 0 $，使得初边值问题存在唯一解：

该定理断言了解的局部存在性、唯一性与正则性，且解的空间正则性随时间导数阶数递减，符合双曲系统的典型特征。值得注意的是，角度条件 $ \omega_n < \pi/[m] $ 是为了防止解在角点附近出现过强的奇性而导致正则性损失；排除某些特定角度则是出于技术性证明需要（见附录A中Proposition 3.9的证明）。

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

a) 核心分析工具

1. 新型各向异性Sobolev空间 $ H^m_* $：这是全文的基础框架。该空间允许在角点附近进行有效的能量估计，其嵌入性质（如 $ H^4_* \hookrightarrow L^\infty $）为非线性项的估计提供了支撑。

2. 对偶方法（Lax–Phillips型）与弱解构造：针对线性化问题，作者首先在高阶切向空间中通过对偶论证建立弱解的存在性。这种方法绕开了经典Galerkin逼近在非光滑区域中的困难。

3. 正则化与迹保持平滑过程：为验证弱解实为强解，作者设计了一系列保持边界迹不变的平滑算子，逐步提升解的正则性，并结合紧性论证得到“弱等于强”的唯一性结果。

4. Div-Curl结构与Helmholtz型分解：这是处理法向导数估计的核心工具。利用MHD方程中的结构，将梯度 $ \nabla u $、$ \nabla b $ 分解为散度与旋度部分：

   $$ u = \nabla \phi + \nabla^\perp \psi, \quad b = \nabla \eta + \nabla^\perp \zeta, $$

   结合边界条件 $ u\cdot\nu = b\cdot\nu = 0 $ 可知 $ \phi,\eta $ 满足Neumann问题，$ \psi,\zeta $ 满足Dirichlet问题。由此可通过估计 $ \nabla \cdot u $、$ \nabla \times u $、$ \nabla \cdot b $、$ \nabla \times b $ 来重构法向导数。

5. 能量估计中的“借乘子”技巧（borrowing multiplier）：在估计高阶切向导数时，会出现难以控制的法向导数项。作者巧妙利用动量方程中压力梯度与速度的关系，将此类项转化为关于 $ \nabla p $ 的内积，并通过选取适当的乘子将其“转移”到压力项上，从而闭合能量不等式（见Lemma 2.10）。

6. Picard迭代与收敛性分析：对于非线性问题，采用Picard迭代格式，依赖于线性化问题的适定性及非线性项在 $ H^m_* $ 中的良好乘积性质（由Corollary 2.1保证），最终在 $ m \ge 8 $ 时实现收敛。

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

• 法向导数估计的障碍：由于角点处多个法向方向共存，传统通过方程代入估计 $ \partial_\nu u $ 的方法失效。作者转而利用MHD特有的结构：从连续性方程得 $ \nabla \cdot u = -(\partial_t + u\cdot\nabla)\log R $，仅含切向导数；从动量方程投影得 $ \nabla p $ 的表达式；再结合 $ \nabla \times u $ 的演化方程，最终通过Helmholtz分解反演出 $ \nabla u $ 的完整信息。

• 高阶切向导数的能量闭合问题：当对解作用高阶切向导数时，会产生涉及法向导数的交换子项。这些项无法直接估计，但恰好出现在能量泛函中与 $ \nabla p $ 的配对。通过引入辅助乘子并利用方程结构，可将这部分“坏项”吸收进已控项中。

• 角点附近的奇性控制：通过构造带有距离权重的切向导数 $ \partial_{w_i} = w_i \cdot \nabla $，自动引入了对边界奇性的抑制机制。例如在角点邻域 $ w_i \sim x_i \nabla $，故 $ \partial_{w_i} $ 相当于加权导数，有效缓解了解在角点处可能的奇性增长。

• 弱解到强解的升级：标准正则化方法会破坏边界条件。为此，作者发展了一套保迹平滑技术，即先在内部用磨光子卷积，再通过求解椭圆边值问题修正边界值，确保每一步近似解仍满足原始边界条件，从而保证极限解的有效性。

与之前工作的比较

a) 突破与改进

本文相较于已有工作实现了多项关键突破：

1. 首次在角点区域建立双曲型系统的适定性：此前关于双曲系统在非光滑区域的研究极为稀少。Osher [18,19]虽在70年代尝试用Kreiss对称化方法处理线性双曲组，但其条件难以验证且存在隐含的正则性损失；Huang–Temam [9,10]仅在矩形区域对浅水波方程获得 $ L^2 $ 适定性，未触及高阶正则性。本文首次在一般角点区域中建立了完整的高阶能量估计与强解存在性，填补了该领域的空白。

2. 推广了各向异性Sobolev空间的应用范围：Chen [5]等人发展的 $ H^m_* $ 空间原本仅适用于光滑边界。本文通过构造几何适应的切向场 $ w_1,w_2 $，将其成功推广至角点情形，赋予其新的几何意义，成为处理此类问题的标准工具之一。

3. 解决了MHD在完全导电边界下的正则性难题：Yanagisawa–Matsumura [28]首次指出，在 $ b\cdot\nu = 0 $ 条件下，标准Sobolev空间中线性化MHD可能是不适定的。他们改用各向异性空间解决了光滑区域的问题。本文在此基础上进一步克服了几何非光滑性带来的额外困难，将这一成果扩展至更具现实意义的多边形区域。

4. 比Godin [7]的结果更强且方法更普适：Godin曾研究二维可压缩Euler方程在角点区域的适定性，但他依赖于方程的特殊结构（如curl-vorticity关系），且要求角度足够小。本文处理的是更强耦合的MHD系统，且通过Div-Curl-Helmholtz框架提供了一种更系统、可推广的方法。此外，本文明确给出了角度阈值 $ \pi/[m] $，而不仅是“充分小”。

5. 对初始数据的要求更清晰合理：相比早期工作中模糊的“光滑性”假设，本文明确指出需初始数据属于 $ H^m_* $ 并满足高阶相容性条件，这既反映了物理合理性（边界行为匹配），也揭示了数学本质（时间导数的迹存在性）。

6. 无需结构简化或耗散机制：Benoit [1,2]等人的工作依赖于特殊的边界结构或添加人工耗散项。本文完全在理想（无粘、无电阻）MHD框架下完成证明，更具物理真实性。

综上，本文不仅在结果上实现了从光滑到非光滑区域的重大跨越，更在方法论上发展出一套适用于非光滑域上双曲系统的通用分析框架，有望推广至弹性力学、相对论流体力学等领域。

3. Extreme internal waves: gravity currents and overturning fronts

极端内波：重力流与翻转锋面

作者: Robin Ming Chen, Samuel Walsh, Miles H. Wheeler

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.05329v1

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

详细分析

中文标题: 极端内波：重力流与翻转锋面

研究对象和背景

本文研究的是二维不可压缩、无旋双层流体系统中的水动力激波（hydrodynamic bores），即一种在稳定分层流体中传播的前沿型行波解。这类解出现在两层密度不同的理想流体之间，界面为自由边界，上下由刚性水平壁限制，忽略表面张力，仅考虑恒定重力作用。

物理上，这种模型广泛用于描述海洋或大气中由于密度跃层（如温跃层或盐跃层）引发的内波现象。当较重流体侵入较轻流体时，可形成所谓的重力流（gravity current），例如火山喷发后的云团、河口羽流或油污扩散等自然过程。数学上，该问题归结为求解双层自由边界欧拉方程组的稳态解。

设上层流体密度为 $\rho_1 > 0$，下层为 $\rho_2 > \rho_1$，分别占据区域 $D_1$ 和 $D_2$，其交界面记为 $S = \partial D_1 \cap \partial D_2$。引入移动坐标系，使得波以恒定速度前进，故解不依赖于时间。令上游极限 $x \to -\infty$ 处的速度场趋于水平常数 $( -c, 0 )$，其中 $c>0$ 是波速；上游厚度分别为 $h_1, h_2 > 0$，总深度归一化为 $h_1 + h_2 = 1$。

通过量纲分析，定义无量纲弗劳德数（Froude number）：

其中 $g > 0$ 为重力加速度。进一步假设每层流体均为不可压且无旋，因此可用伪流函数 $\psi$ 描述速度场：$\mathbf{v}_i = (\partial_y \psi, -\partial_x \psi)$ 在第 $i$ 层成立，且满足拉普拉斯方程。

整个系统的控制方程如下：

最后一个方程是动态边界条件，来源于压力连续性和伯努利定律。值得注意的是，对于存在激波解的情形，弗劳德数 $F$ 必须取特定值：

这是通过对共轭流分析得出的唯一允许非平凡激波存在的参数值。

文中还特别考虑了Boussinesq极限，即 $\rho_1 \approx \rho_2$ 的情形。在此极限下，密度差极小，可近似设 $\rho_1 = \rho_2 =: \rho$，同时将动态条件简化为：

此类解被称为“光滑水动力激波”，并分为两类：抬升激波（elevation bore），其界面在上游达到最小高度；以及下沉激波（depression bore），其上游高度最大。作者此前的工作 [11] 已通过全局分支理论构造出两条大振幅解曲线：$C_{\text{elev}}$ 和 $C_{\text{depr}}$。

研究的关键在于理解这些解族在振幅趋于极端时的行为。具体而言，是否存在翻转（overturning）——即界面出现垂直切线？或者是否趋向于某种奇异构型，如接触刚性壁形成重力流？

为此，文章引入了一种更广义的解概念：域变分解（domain variational solution）。这类解不要求经典正则性，而是作为某个能量泛函在内部变分下的临界点来定义，允许界面出现尖点或接触边界的情况。这使得可以在极限过程中处理非光滑几何结构。

此外，问题的几何特性极为重要。例如，在可能的双重停滞点（double stagnation point）处，上下两层流体的速度同时趋于零，导致动态条件退化，标准的自由边界正则性理论失效。此时需要发展新的几何分析工具，特别是关于速度场在停滞点附近的衰减速率估计。

空间正则性方面，经典解要求 $\psi \in C^{2+\alpha}_b(D_1) \cap C^{2+\alpha}_b(D_2) \cap C_b(D)$，$\eta \in C^{2+\alpha}(\mathbb{R})$，其中 $\eta(x)$ 表示界面 $S$ 的函数形式。但在极限情况下，仅能期望 $\eta$ 为Lipschitz连续，甚至可能出现非图表示的过顶结构。

综上所述，本文的研究对象是一个高度非线性的自由边界椭圆系统，具有物理背景明确的双层分层结构，关注其大振幅极限行为，尤其是界面是否翻转或趋近于重力流配置。这一问题长期困扰学界，因传统方法难以刻画解分支的最终形态。

主要定理或主要结果

本文证明了三个核心定理，揭示了水动力激波在大振幅极限下的确切行为。

定理 1.1（翻转激波）：对任意 $0 < \rho_2 \leq \rho_1$，存在一条全局连续曲线

由满足正则性条件 (1.4) 的经典解构成。若 $\rho_2 < \rho_1$，则沿此曲线极限 $s \to -\infty$ 时，界面必然发生翻转：

这意味着抬升激波族在极限下界面斜率无界，出现垂直切线。

定理 1.2（重力流或翻转）：对同一参数范围，存在另一条全局曲线

对应下沉激波。若 $\rho_2 < \rho_1$，则当 $s \to \infty$ 时，要么发生翻转：

要么存在有界序列 $\{x_n\} \subset \mathbb{R}$ 和 $s_n \nearrow \infty$，使得

其中 $(\psi, S)$ 是重力流问题 (1.5) 的非平凡域变分解，且 $S$ 为 Lipschitz 函数图像，与上壁相切接触。

定理 1.3（Boussinesq极限下的接触角）：在 $\rho_1 = \rho_2$ 的 Boussinesq 极限下，

• 沿 $C_{\text{elev}}$，若不翻转，则极限为重力流，界面与下壁成60°角；
• 沿 $C_{\text{depr}}$，若不翻转，则极限为重力流，界面与上壁成60°角。

这首次严格证明了冯·卡门（von Kármán）关于重力流接触角为 $60^\circ$ 的著名猜想。

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

a) 核心分析工具

论文采用的主要数学框架是几何分析与自由边界正则性理论的结合。核心工具包括：

首先，域变分结构被用来定义弱解。作者借鉴 Va˘rva˘ruca 与 Weiss [54] 的思想，将流函数 $\psi$ 视为某能量泛函的临界点，允许在停滞点附近进行内变分。这种方法绕开了经典解必须满足的高阶正则性要求，使极限过程得以闭合。

其次，单调性公式（monotonicity formula）起着决定性作用。受 Alt–Caffarelli–Friedman 方法启发，作者构建适用于双相问题的频率型公式，用以控制不同尺度下动能的分布。尽管经典的乘积型单调性不适用，但通过对称化和平均化技巧，仍可获得有效信息。

最关键的技术是吹胀分析（blowup argument）。考虑可能的双重停滞点（即上下流体速度均趋于零的点），通过缩放变换 $\psi_r(x,y) = \psi(rx, ry)/r^{3/2}$ 分析其渐近行为。保持动态条件不变的自然缩放指数为 $3/2$。若原问题不存在翻转，则可通过紧性提取收敛子列至某个齐次极限解。

然而，直接应用 [54] 中的方法面临根本障碍：在单相水中，Bernoulli 定律给出 $|\nabla \psi|^2 = O(y)$ 的逐点衰减估计（称为 Bernstein 型不等式），但在双层情形中，动态条件仅控制速度跳跃，而非各相速度大小。为此，作者发展了一个全新的平均衰减估计：

该估计对两类解族均成立，是全文最关键的创新之一。它通过利用解的单调性结构（如水平方向单调、垂直方向符号固定）和椭圆估计得到，而非依赖复分析或共形映射。

接着使用Oddson 比较原理建立梯度的逐点下界：在锥形区域 $D_\varepsilon^\pm$ 内，

此下界与上述上界矛盾，从而排除双重停滞的可能性，迫使翻转发生。

最后，对于重力流极限，作者引入一个辅助延拓问题：将接近上壁的界面“打开”，将流场延拓为均匀层流。这样可在远离真实界面的区域应用标准椭圆理论，统一控制 $\psi$ 的 Lipschitz 范数，避免因几何塌缩导致估计失效。

b) 技术难点与解决技巧

主要技术挑战在于处理退化自由边界和潜在停滞点。

第一类困难来自动态条件在停滞点处的奇异性。当 $|\nabla \psi_1|, |\nabla \psi_2| \to 0$ 时，方程失去椭圆性，标准 Schauder 估计失效。作者通过引入变分结构，将问题转化为几何测度论框架下的极小曲面类问题，从而启用 Blow-up 和单调性方法。

第二类困难是双相系统缺乏有效的逐点能量估计。为克服这一点，作者创造性地使用积分平均法代替点态分析。他们证明：只要界面是单调 Lipschitz 图像，即使没有显式表达式，也能通过对 Green 恒等式和 Hardy 不等式的精细估计，得出动能的平均衰减速率。

第三，为了处理界面趋近刚性壁的情形，常规方法会因域退化而崩溃。作者设计的延拓技巧巧妙规避了这一问题：将上层流体外推为上游层流状态，并证明这种人工构造不会破坏关键估计。这体现了对物理结构的深刻洞察——远场行为主导局部正则性。

此外，作者充分利用了激波解的定性性质：如水平单调性、垂直对称性破缺、上下游不对称等。这些先验性质不仅简化了分析，也为吹胀极限提供了额外约束，确保极限对象具有足够结构以便分类。

综上，本文融合了现代偏微分方程中的多个高级工具：变分法、几何测度论、椭圆正则性理论、渐近分析，并加以创新改造，成功解决了长期悬而未决的翻转问题。

与之前工作的比较

a) 突破与改进

本文相较以往研究实现了多项根本性突破。

首先，在翻转现象的严格证明方面，此前虽有大量数值证据（如 Dias 与 Vanden-Broeck [17]）显示内波可能翻转，但始终缺乏数学证明。尽管 Constantin–Va˘rva˘ruca–Strauss [12] 等人构造了含涡度的周期波族并允许过顶结构，但他们使用的是直接构造法，而非从平滑解分支延续而来。本文则是首次在全局分支框架下，对O(1) 重力强度的无涡双层系统，证明了解族必然经历翻转。这填补了从光滑解到奇异构型之间的理论空白。

其次，针对重力流极限中的接触角问题，冯·卡门早在1940年就基于复势展开提出 $60^\circ$ 猜想，后经 Benjamin [4] 等人推广，但始终停留在形式推导层面。本文首次给出了严格的数学证明，且适用于广泛的解类，不限于小扰动或特殊对称情况。更重要的是，证明不依赖复分析，因而更具普适性，可推广至含涡情形。

再者，本文发展的平均衰减估计

是对自由边界正则性理论的重要补充。以往工作（如 [54]）依赖 $|\nabla u|^2 = O(\text{dist})$ 的点态估计，而这在双相问题中不可得。本文通过放弃点态控制、改用积分平均，既保留了足够信息用于吹胀分析，又适应了物理系统的本质限制。这一思路有望应用于其他多相流或复合材料界面问题。

最后，本文所用方法大幅降低了对解正则性的先验假设。不同于早期研究需假设解解析或高度光滑，本文仅需 Lipschitz 连续性即可启动分析。这使得结论适用于更广泛的极限过程，增强了结果的鲁棒性。

相比之下，Dávila 等人 [15] 使用胶合方法构造小重力下的过顶孤立波，虽具开创性，但局限于弱非线性 regime，无法回答有限重力下是否存在翻转。而本文恰恰在 $F = O(1)$ 条件下完成证明，更具物理现实意义。

总之，本文不仅解决了两个经典猜想（翻转与 $60^\circ$ 接触角），更提供了一套可用于分析大振幅波极限行为的新工具链，标志着自由边界水波理论的重大进展。

4. Existence and weak-strong uniqueness of measure solutions to Euler-alignment/Aw–Rascle–Zhang model of collective behaviour

欧拉对齐/Aw–Rascle–Zhang集体行为模型测度解的存在性与弱-强唯一性

作者: Jakub Woźnicki, Ewelina Zatorska

类别: math.AP, 35Q35, 35Q70, 35A02, 35B35

发布日期: 2025-11-07

arXiv ID: 2511.05326v1

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

详细分析

中文标题: 欧拉对齐/Aw–Rascle–Zhang集体行为模型测度解的存在性与弱-强唯一性

研究对象和背景

本文研究的是多维 欧拉对齐系统（Euler-alignment system）及其形式等价的非局部 Aw–Rascle–Zhang（ARZ）交通流模型。这两个系统均用于描述具有集体行为（collective behaviour）的宏观动力学，例如动物群集、人群运动或交通流中的自组织现象。

所研究的核心方程是如下压力为零的可压缩欧拉型系统：

其中：

• $\rho(t,x)$ 表示在时间 $t$ 和位置 $x \in \mathbb{R}^d$ 处的密度；
• $u(t,x)$ 是相应的速度场；
• $K : \mathbb{R}^d \to \mathbb{R}$ 是一个给定的通信势函数，其Hessian矩阵 $D^2K = \left( \partial^2_{x_i x_j} K \right)_{i,j=1}^n$ 构成所谓的矩阵值通信核（matrix-valued communication kernel），该结构由 Shu 与 Tadmor [51] 引入，用以建模个体之间的预判动力学（anticipation dynamics）。

这种 Hessian 型核的关键特征在于它能诱导出类似于“弹性力”的相互作用，从而促进一致运动（flocking behavior），尤其是在吸引力主导的情形下。

值得注意的是，上述系统没有经典的压力项（即无 $p(\rho)$ 项），因此被称为压力为零系统（pressureless system）。这使得标准的流体力学分析工具（如基于熵不等式或 Sobolev 嵌入的能量方法）难以直接应用，因为缺乏对密度的紧性控制。

一个重要的观察是：通过引入一个新的变量——期望速度场 $w$，定义为：

可以将原欧拉对齐系统重写为如下形式的非局部 ARZ 系统：

这一形式与经典的 Aw–Rascle–Zhang 交通模型高度相似。在一维情形中，传统 ARZ 模型包含守恒律形式的质量与广义动量方程，且 $w - u = p(\rho)$ 被解释为一种“压力”或“驾驶偏好偏移”，通常取为密度的单调增函数（如 $p(\rho) = \rho^\gamma$）。但在本文中，这个偏移是由非局部梯度卷积 $(\nabla K)\star\rho$ 给出的，体现了长程交互与前瞻性调整机制。

本文的研究区域为全空间 $\mathbb{R}^d$，未加边界条件（因为空间无界），初始数据满足一定的可积性和矩条件。具体而言，假设初始密度 $\rho_0 \in \mathcal{P}(\mathbb{R}^d)$（概率测度空间），初始速度 $u_0 \in L^2(\mathbb{R}^d; d\rho_0)$，并满足二阶矩有限：

此外，要求通信核 $K$ 满足以下部分或全部条件：

1. $D^2K \in C_b(\mathbb{R}^d)$：Hessian 有界连续；
2. $D^2K$ 为偶函数且半正定；
3. $D^2K$ 是 Lipschitz 连续的。

这些假设确保了非局部对齐力的良好定义性以及能量估计的有效性。

由于系统缺乏压力项，传统的弱解框架面临严重挑战，尤其是质量集中（concentration effects）可能导致动量方程中出现奇异项。为此，作者采用测度解（measure solution）的概念，在 Radon 测度空间中构建解。此时，密度 $\rho_t$ 是一族随时间演化的概率测度，而动量方程中的非线性项可能产生一个额外的缺陷测度 $\mu_t$，用来捕捉速度场的集中现象。

文中还引入了 flat metric $d_f$ 来度量测度之间的距离，并利用 Wasserstein-2 距离 $W_2$ 分析解的稳定性。函数空间方面，主要涉及：

• 测度空间 $M(\mathbb{R}^d)$、正测度空间 $M_+(\mathbb{R}^d)$、概率测度空间 $\mathcal{P}(\mathbb{R}^d)$；
• 向量/矩阵值测度空间 $M(\mathbb{R}^d; \mathbb{R}^{d\times d}_{\rm sym})$；
• 时间连续路径空间 $C([0,T]; (\mathcal{P}(\mathbb{R}^d), d_f))$ 或 $C([0,T]; (\mathcal{P}_2(\mathbb{R}^d), W_2))$；
• $L^p$ 空间关于测度 $d\rho_t$ 的加权版本。

综上所述，本文研究的是一个在全空间 $\mathbb{R}^d$ 上定义的、无压力、含非局部 Hessian 核对齐力的双曲型系统，目标是在较弱的正则性条件下建立全局测度解的存在性，并证明其与经典解的一致性。

主要定理或主要结果

定理 2.6（测度解的存在性）

设初始数据 $(\rho_0, u_0)$ 满足 Assumption 2.1，且核 $K$ 满足 Assumption 2.2(i)，即 $D^2K \in C_b(\mathbb{R}^d)$。则存在一组三元组 $(\rho, u, \mu)$，构成系统 (1.1) 的测度解，满足 Definition 2.4 中的所有条件，包括：

• 连续性方程 (2.1) 成立；
• 动量方程 (2.2) 在缺陷测度 $\mu$ 下成立；
• 能量不等式 (2.3) 成立；若 $D^2K$ 半正定，则更强的能量等式 (2.4) 成立。

进一步地，此解可通过一个退化粘性 Navier–Stokes 近似的零粘性极限得到。

同时，令 $w$ 由关系式：

几乎处处定义。则存在矩阵值测度 $\nu$，使得四元组 $(\rho, w, u, \nu)$ 构成非局部 ARZ 系统 (1.3)-(1.4) 的测度解（按 Definition 2.5）。

定理 2.7（弱-强唯一性）

设核 $K$ 满足 Assumptions 2.2(i)-(iii)，即 $D^2K \in C_b \cap {\rm Lip}$ 且为偶函数、半正定。假设存在一组强解 $(r,v) \in C([0,T]; \mathcal{P}(\mathbb{R}^d)) \times L^\infty(0,T; W^{1,\infty}(\mathbb{R}^d))$，满足 $r > 0$，并求解：

又设一列初始数据 $(\rho^n_0, u^n_0)$ 满足：

• $\|\rho^n_0 - r_0\|_{\rm TV} \to 0$，
• $\||x|^2(\rho^n_0 - r_0)\|_{\rm TV} \to 0$，
• $\displaystyle\int |u^n_0 - v_0|^2 \, d\rho^n_0 \to 0$。

令 $(\rho_n, u_n, \mu_n)$ 为对应初值的测度解序列（由 Theorem 2.6 得到）。则存在常数 $C > 0$（依赖于 $T, \|D^2K\|_{W^{1,\infty}}, \|v\|_{W^{1,\infty}}$），使得对几乎所有 $t \in (0,T)$，

特别地，当初始误差趋于零时，测度解收敛到强解：

推论 2.8 与 2.9（形式等价性的严格验证）

若存在经典强解，则任意从相同初始数据出发的测度解必然与其一致（密度相等、速度 a.e. 相等、缺陷测度为零）。由此可得：欧拉对齐系统与非局部 ARZ 系统不仅形式等价，而且在共同的粘性逼近下生成相同的极限解，并在强解存在时完全重合。

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

a) 核心分析工具

1. 退化粘性逼近法（Degenerate pressureless Navier–Stokes approximation）
   作者并未使用粒子逼近或动理学层次过渡的方法，而是构造了一个带有密度依赖粘性的近似系统：

   $$ \begin{cases} \partial_t \rho^N + \mathrm{div}_x (\rho^N u^N) = 0, \\ \partial_t (\rho^N u^N) + \mathrm{div}_x (\rho^N u^N \otimes u^N) - \frac{1}{N} \mathrm{div}_x (\rho^N D u^N) + \rho^N u^N (D^2K \star \rho^N) - \rho^N D^2K \star (\rho^N u^N) = 0. \end{cases} $$

   其中粘性项为 $\frac{1}{N} \mathrm{div}_x (\rho^N D u^N)$，属于Bresch–Desjardins 类型粘性。这类粘性具有特殊结构，允许推导出BD 熵估计，从而获得 $\sqrt{\rho^N} \nabla u^N \in L^2$ 的一致估计，进而提升密度的正则性（如 $\nabla \sqrt{\rho^N} \in L^\infty L^2$）。

2. Bresch–Desjardins 熵估计
   利用该估计可得：

   $$ \sup_N \|\nabla \sqrt{\rho^N}\|_{L^\infty(0,T; L^2)} < \infty, $$

   结合 Sobolev 嵌入，得到 $\rho^N$ 在 $L^p_{\rm loc}$ 中的紧性（$p < \infty$），这对处理非线性项至关重要。

3. Mellet–Vasseur 估计
   提供对速度场的 $L^\infty L^2(d\rho^N)$ 控制，结合 BD 估计，可推出 $u^N$ 在加权空间中的有界性，从而保证 $\rho^N u^N$ 的紧性。

4. Flat metric 与 Wasserstein 距离下的紧性理论
   使用 flat metric $d_f$ 和 Wasserstein 距离 $W_2$ 来控制测度序列的弱收敛。特别是借助 Proposition 3.7，若测度族tight 且总变差一致有界，则 flat 收敛等价于弱收敛。这为通过近似解提取极限提供了拓扑框架。

5. 相对熵方法（Relative entropy method）
   用于证明弱-强唯一性。不同于传统依赖凸熵势的方法（如 Dafermos 对标量守恒律的工作），本文采用基于 Wasserstein 距离平方 的泛函作为“相对熵”，构造 Lyapunov 型量：

   $$ \mathcal{E}(t) := \int |u_n - v|^2 \, d\rho_n + W_2^2(\rho_n, r) + \mathrm{tr}(\mu_n(\mathbb{R}^d)), $$

   并证明其满足微分不等式 $\frac{d}{dt} \mathcal{E}(t) \leq C \mathcal{E}(t)$，从而由 Gronwall 不等式得出唯一性。

6. 缺陷测度控制技术
   在动量方程中，对流项 $\rho^N u^N \otimes u^N$ 可能在极限下产生集中效应。作者通过能量不等式与近似系统的耗散结构，控制该集中表现为一个正半定矩阵测度 $\mu$，并在弱-强唯一性中将其纳入估计。

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

1. 无压力导致的紧性缺失
   传统 Euler 方程中压力提供对密度的控制，但此处无压力。解决办法是利用密度依赖粘性带来的 BD 熵，获得 $\nabla \sqrt{\rho}$ 的一致估计，从而避免使用压力相关的 Sobolev 方法。

2. 非局部对齐项的收敛问题
   非局部项 $\int D^2K(x-y)(u(x)-u(y)) \rho(y) dy$ 含双重积分，需处理乘积测度 $\rho \otimes \rho$ 的弱收敛。作者利用 $u \in L^\infty(L^2(d\rho))$ 与 tightness，结合 $D^2K$ 的光滑性，通过截断和逼近论证完成极限过程。

3. 缺陷测度的识别与消除
   当 $u^N$ 出现震荡或集中时，$\rho^N u^N \otimes u^N$ 的弱极限小于 $\rho u \otimes u$，差值即为缺陷测度 $\mu$。作者通过能量不等式保留该项，并在弱-强唯一性中证明：若强解存在，则 $\mu \equiv 0$。

4. 测度解中速度的定义方式
   由于 $\rho$ 可能奇异（如 Dirac 浓缩），速度 $u$ 仅在 $d\rho_t$ 意义下定义。作者始终在 $L^2(d\rho_t)$ 框架下操作，避免逐点定义问题。

5. ARZ 与欧拉对齐之间的桥梁构建
   关键在于验证变换 $w = u + \nabla K \star \rho$ 在测度解意义下仍保持方程等价性。作者通过对偶测试函数方法，在极限中恢复动量方程 $\partial_t(\rho w) + \mathrm{div}_x(\rho w \otimes u) = 0$，并引入新缺陷测度 $\nu$ 来吸收残余项。

与之前工作的比较

a) 突破与改进之处

1. 首次在全空间建立全局测度解
   此前多数关于非局部欧拉系统的工作（如 [12], [14]）集中在有界或周期域，而本文在 $\mathbb{R}^d$ 上工作，更贴近实际应用场景（如开放空间中的群体运动），也带来更大的技术挑战（如无穷远处的行为控制）。

2. 无需小初值假设
   许多已知结果（如 [11], [20], [52]）要求初始数据足够小才能保证整体适定性。本文通过粘性逼近与能量估计，完全去除小初值限制，仅需初始二阶矩与能量有限。

3. 统一处理欧拉对齐与非局部 ARZ
   尽管 [18] 观察到两者的形式等价性，但缺乏严格的数学联系。本文首次证明二者可通过同一粘性逼近获得，并在弱-强唯一性下完全一致，实现了两种模型的动态等价性验证。

4. 推广至矩阵值通信核
   以往多数对齐模型使用标量权重 $\Psi(|x-y|)$，而本文考虑 $D^2K$ 作为矩阵核，能描述更复杂的交互方向性（如预判方向、弹性响应），更具物理真实性。

5. 弱-强唯一性适用于压力为零系统
   传统相对熵方法依赖压力势的凸性，而本文创新性地使用 Wasserstein 距离 + 加权速度差 构造相对熵，成功应用于无压力情形，拓展了该方法的适用范围。

6. 避免凸积分方法的病态性
   [16] 使用凸积分构造了无限多个弱解，说明局部 ARZ 系统不适定。本文通过引入粘性逼近和缺陷测度，选择物理上有意义的解（来自粘性极限），并通过弱-强唯一性排除其他非物理解。

综上，本文在解的存在性范围、模型等价性的严格性、方法的普适性等方面均超越了已有文献，为集体行为模型的宏观分析提供了坚实基础。

5. On global regular axially-symmetric solutions to the Navier-Stokes equations in a cylinder

关于柱体内Navier-Stokes方程的整体正则轴对称解的研究

作者: Wiesław J. Grygierzec, Wojciech M. Zajączkowski

类别: math.AP, 35A01, 35B01, 35B65, 35Q30, 76D03, 76D05

发布日期: 2025-11-07

arXiv ID: 2511.05098v1

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

详细分析

中文标题: 关于柱体内Navier-Stokes方程的整体正则轴对称解的研究

研究对象和背景

本文研究的是三维不可压缩 Navier-Stokes方程 在有限圆柱区域 $\Omega \subset \mathbb{R}^3$ 上的 轴对称整体正则解 的存在性与先验估计问题。所考虑的方程组为：

其中 $v = v(x,t)$ 是速度场，$p = p(x,t)$ 是压力函数，$f$ 是外力项，$\nu > 0$ 是粘性系数。区域 $\Omega$ 定义为有限圆柱体：

其边界分为两部分：侧壁 $S_1$ 和上下底面 $S_2$，即：

由于问题具有 轴对称性（axially symmetric），作者采用柱坐标 $(r,\varphi,z)$ 进行分析，其中：

在该坐标系下，速度向量可分解为：

其中 $v_r, v_\varphi, v_z$ 分别表示径向、周向（旋流）和轴向分量。

对应的 涡量（vorticity）定义为：

其各分量由速度梯度导出，并满足特定表达式（见原文(1.5)）。

本文的关键在于所设定的 非标准边界条件，这在经典无滑移边界条件下难以获得全局估计。具体地，边界条件如下：

• 在侧壁 $S_1$ 上：

  $$ v_r = v_\varphi = \omega_\varphi = 0; $$

• 在上下底面 $S_2$ 上：

  $$ v_z = \omega_\varphi = \partial_z v_\varphi = 0. $$

这些边界条件不同于传统的 无滑移条件（$v=0$ on $\partial\Omega$），而属于一种“混合型”或“自然型”边界条件，最早由 Ladyzhenskaya 提出，在使用 涡量-流函数方法（vorticity-stream function formulation）时显得更为合适。例如，$\omega_\varphi|_{S}$ 的出现是推导过程中自然产生的量，因此设定其为零有助于控制高阶导数。

此外，初始条件为：

其中 $v_0$ 是给定的无散度向量场，且满足上述边界条件。

为了简化分析并聚焦于关键结构，作者引入了所谓的 旋流变量（swirl variable）：

这一变换能够有效消除柱坐标中因 $r=0$ 轴上的奇异性带来的技术困难。类似地，他们还引入了归一化后的涡量分量：

并证明这两个量满足一个耦合的抛物型偏微分系统（见原文(1.11)-(1.12)），这是整个估计过程的核心。

值得注意的是，由于轴对称解可能在对称轴 $r=0$ 处产生奇异性，作者引用了 Liu-Wang 的 渐近展开理论（见(1.21)），假设解在 $r=0$ 附近具有光滑展开形式，如：

以及相应的流函数 $\psi$ 和修正流函数 $\psi_1 = \psi / r$ 的展开。这类展开保证了即使在 $r=0$ 处，物理量仍保持有界性和正则性，从而避免伪奇点。

函数空间方面，文章主要工作在 能量类空间 中进行，特别是定义了一个关键的能量范数：

其具体形式为：

即时间上 $L^\infty$、空间上 $L^2$ 的范数加上梯度的空间时间 $L^2$ 范数。这种范数常用于控制解的正则性。

目标是建立如下类型的 全局先验估计：

其中 $\phi$ 是仅依赖于初值和外力的数据（记作“data”）的一个递增正函数。这个估计一旦成立，结合局部解的存在性和连续性论证，即可推出 整体正则解的存在性。

特别指出的是，作者明确说明：“我们无法对无滑移边界条件推导任何全局类型估计”，突显了当前边界条件选择的技术必要性。

综上所述，本文的研究对象是一个在有限圆柱域中带有特殊混合边界条件的轴对称 Navier-Stokes 方程系统，通过引入归一化涡量变量和流函数方法，旨在克服 $r=0$ 奇异性与非线性耦合带来的挑战，最终实现整体正则解的构造。

主要定理或主要结果

论文的主要结果体现在 Theorem 1.2 中，尽管它以一种间接方式呈现，但实质上构成了全文的核心结论。该定理断言：

Theorem 1.2. 设 Notation 1.1 中定义的所有量均为有限。若存在常数 $c_0 > 0$，使得以下两个条件之一成立：

1. （大旋流情形）$\lambda(\sigma') \geq c_0$，其中 $\sigma' > 6$，
2. （小旋流情形）$\lambda(\sigma') \leq c_0$，

其中 $\lambda(s)$ 定义为：

则存在仅依赖于数据和 $c_0$ 的递增正函数 $\phi_1, \phi_2$，使得能量泛函

满足估计：

此结果意味着无论旋流 $v_\varphi$ 相对于 $L^s$ 范数是“大”还是“小”，都可以获得统一的先验界。结合局部解延拓思想，这足以推出 整体正则解的存在性。

更进一步，文中多个引理构成支撑该定理的关键步骤：

• Lemma 3.2 给出当 $\lambda(\sigma') \geq c_0 > 0$ 时，$X(t)$ 可被 $\phi_1(1/c_0, \mathrm{data})$ 控制；
• Lemma 4.1 和 Remark 4.2 表明当 $\lambda(\sigma') \leq c_0$ 时，$X(t)$ 可被 $\phi_2(c_0, \mathrm{data})$ 控制；
• 因此，不论哪种情况发生，总能获得一致界。

最终结论虽未直接写成“存在唯一的整体正则解”，但从上下文及引用文献 [Z1][Z2][OZ] 可知，此类估计一旦建立，配合局部存在性与唯一性理论，即可推出整体适定性。

因此，本文的主定理实质为：

在所设边界条件下，对于足够正则的初值 $v_0$ 和外力 $f$，三维轴对称 Navier-Stokes 方程在有限圆柱域中存在唯一的整体正则解。

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

a) 核心分析学工具

本文采用了多种现代偏微分方程分析技术，围绕 能量估计 展开，构建了一套系统的先验估计框架。主要工具包括：

1. 能量不等式（Energy Inequalities）
   最基础的是经典的动能估计（Lemma 2.4）：

   $$ |v(t)|_{2,\Omega}^2 + \nu \int_{\Omega_t} \left( |\nabla v|^2 + \frac{v_r^2 + v_\varphi^2}{r^2} \right) dx dt \leq C D_1^2, $$

   其中 $D_1$ 包含初值和外力的 $L^2$ 范数。该项提供了速度的基本能量控制。

2. 最大值原理（Maximum Principle）
   对旋流变量 $u = r v_\varphi$ 所满足的输运-扩散方程（见(1.9)）应用最大值原理，得到：

   $$ |u(t)|_{\infty,\Omega} \leq D_2, $$

   即旋流在 $L^\infty$ 意义上有界。这是控制非线性项的关键一步。

3. 椭圆正则性估计（Elliptic Regularity Estimates）
   利用修正流函数 $\psi_1 = \psi / r$ 满足的方程：

   $$ -\Delta \psi_1 - \frac{2}{r} \partial_r \psi_1 = \Gamma, $$

   并结合边界条件 $\psi_1|_S = 0$，应用 H² 和 H³ 椭圆估计（Lemmas 2.7–2.8），得到：

   $$ \|\psi_1\|_{H^2(\Omega)} \lesssim |\Gamma|_{2,\Omega},\quad \|\psi_1\|_{H^3(\Omega_t)} \lesssim |\Gamma_{,z}|_{2,\Omega_t}. $$

   这些估计将流函数的高阶正则性转化为对 $\Gamma = \omega_\varphi / r$ 的控制。

4. Hardy 不等式与插值不等式
   文中多次使用 Hardy 型不等式（Lemma 2.1）来处理 $1/r$ 奇性，例如：

   $$ \left\| \frac{f}{r} \right\|_{L^p(\Omega)} \leq C \|\nabla f\|_{L^p(\Omega)}, \quad \text{当 } f|_{r=R} = 0. $$

   此外，Sobolev 插值（Lemma 2.2）和 Hardy 插值（Lemma 2.3）被用来连接不同 $L^p$ 范数，特别是在处理非线性项时至关重要。

5. Gronwall-type 推理与连续性方法
   所有估计最终通过 Gronwall 引理整合，结合局部解的存在性，利用连续性延拓论证得出整体解的存在性。

b) 技术难点与应对技巧

1. $r=0$ 处的奇异性处理
   柱坐标下 $1/r$ 项在 $r=0$ 发散，是主要障碍。作者通过两种方式克服：

   • 引入归一化变量 $\Phi = \omega_r / r$, $\Gamma = \omega_\varphi / r$，使原奇异项变为有界；
   • 假设解在轴附近具有 Liu-Wang 渐近展开（1.21），确保所有物理量在 $r=0$ 解析光滑。

2. 非线性项的控制
   非线性项如 $v \cdot \nabla \Phi$ 在估计中需精细处理。作者通过将速度用 $\psi_1$ 表示（见(1.19)），并将梯度项拆解为低阶与高阶部分，再结合 Sobolev 嵌入与插值，实现控制。

3. 旋流主导情形的分离讨论
   关键创新在于引入 $\lambda(s)$ 函数，将问题分为“大旋流”与“小旋流”两类。在每种情形下采用不同的估计策略：

   • 当旋流较大时，利用其主导效应压制其他不稳定项；
   • 当旋流较小时，则退化为接近二维情形，可用更强的嵌入不等式。

4. 边界条件的设计优势
   特殊边界条件（如 $\omega_\varphi = 0$ on $S$）使得在边界上许多高阶导数自动消失，极大简化了积分边界项的处理，否则将在分部积分中引发额外困难。

综上，本文通过变量替换、函数空间选取、分类讨论与经典不等式组合，成功绕开了传统方法难以突破的瓶颈。

与之前工作的比较

本文相较于先前研究做出了显著改进，主要体现在 放宽假设条件 与 强化结论适用范围 两个方面。

首先，对比作者 Zajączkowski 本人早期工作 [Z1][Z2]，本文去除了两个关键限制：

1. 取消了 Serrin 型附加条件
   在 [Z1][Z2] 中，要求：

   $$ \psi_1|_{r=0} = 0 \quad \text{和} \quad \frac{|v_\varphi|_{s,\infty,\Omega_t}}{|v_\varphi|_{\infty,\Omega_t}} \geq c_0 > 0,\ \forall s > 0, $$

   即强制旋流在轴上为零且在所有 $L^s$ 空间中均匀有界远离零。这些条件缺乏物理动机，仅为数学便利而设。本文通过引入 $\lambda(s)$ 的分类讨论，完全摆脱了这些人为约束。

2. 推广边界条件适用性
   [Z1] 使用周期性边界条件，[Z2] 虽然也考虑柱域，但仍需额外控制条件。本文直接在物理意义更合理的混合边界条件下完成证明，增强了模型的实际可解释性。

其次，相比 Ozanski-Zajączkowski [OZ] 的工作，后者虽已去除 $\psi_1|_{r=0}=0$ 条件，但仍保留旋流非退化条件（即 $\lambda(s) \geq c_0$）。本文进一步将其推广至允许旋流任意小的情形，实现了真正意义上的 无条件整体正则性（在给定边界条件下）。

此外，与 Constantin-Foias-Zhou [CFZ] 等经典工作相比，本文不仅继承了他们关于 $(\omega_r/r, \omega_\varphi/r)$ 控制速度的思想，还通过更精细的能量层级分析，避免了对初始数据过高正则性的依赖。

最后，本文明确指出：在标准无滑移边界条件下目前尚无法获得类似全局估计，这反衬出所选边界条件的优越性——它既保留物理合理性，又具备数学可控性。

综上，本文的突破在于：首次在无需额外非退化或消失性假设的前提下，证明了有限圆柱域中轴对称 Navier-Stokes 方程的整体正则解存在性，标志着该方向研究的重要进展。

6. Waterborne epidemics via a new coupled SIR–Pathogen–Navier-Stokes system: Mathematical modeling, nonlinear analysis and numerical simulation

基于新型耦合SIR–病原体–纳维–斯托克斯系统的水媒传染病传播：数学建模、非线性分析与数值模拟

作者: Mohamed Mehdaoui, Yassine Ouzrour

类别: math.AP, 35K57, 35Q92, 49K20, 76Zxx, 92D25, 92D30

发布日期: 2025-11-06

arXiv ID: 2511.04841v1

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

详细分析

中文标题: 基于新型耦合SIR–病原体–纳维–斯托克斯系统的水媒传染病传播：数学建模、非线性分析与数值模拟

研究对象和背景

本文研究一类由流体介质驱动的水媒传染病（waterborne epidemics）的空间传播机制，提出并分析一个全新的耦合偏微分方程系统——SIR–Pathogen–Navier–Stokes（简称 SIRPNS）系统。该模型将经典的流行病学SIR模型与描述流体运动的不可压缩纳维–斯托克斯方程（Navier–Stokes equations）相结合，并引入环境中的病原体浓度作为中间变量，从而构建了一个包含宿主动力学、病原体传输与流体动力学三者相互作用的完整框架。

传统的传染病模型如Kermack-McKendrick提出的SIR模型仅考虑个体之间的直接接触传播，难以刻画通过水或空气等流体媒介进行间接传播的情形。对于诸如霍乱、痢疾等依赖水源传播的疾病，环境中存活的病原体可随水流扩散，导致感染在空间上呈现非均匀分布。已有研究通常将病原体视为被动标量（passive scalar），仅受扩散或衰减影响，而忽略了其对流体物理性质（如粘度、密度）可能产生的反馈效应。本文的核心创新在于同时考虑两个关键机制：

1. 流体输运效应：宿主密度和病原体浓度均受到流体速度场 $U(x,t)$ 的对流作用；
2. 双向耦合机制：流体的有效粘度 $\nu(C)$ 是病原体浓度 $C$ 的函数，即病原体的积累会改变流体的力学行为，形成从流行病学到流体力学的反向反馈。

设 $\Omega \subset \mathbb{R}^d$（$d=2$ 或 $3$）为有界光滑区域，代表所研究的地理或水域空间；时间区间为 $(0,T)$，其中 $T>0$。定义如下变量：

• $S(x,t)$: 易感者密度
• $I(x,t)$: 感染者密度
• $R(x,t)$: 康复者密度
• $C(x,t)$: 环境中病原体浓度
• $U(x,t)$: 流体速度场
• $p(x,t)$: 流体压力

总人口密度为 $N(x,t) = S(x,t) + I(x,t) + R(x,t)$。

根据文中假设（A₁–A₆），完整的SIRPNS系统如下所示：

边界条件为：

初始条件为：

值得注意的是，该系统的非线性来源于多个方面：

• Navier-Stokes 方程本身的非线性对流项 $(U \cdot \nabla)U$；
• 病原体浓度依赖的粘度函数 $\nu(C)$ 出现在应力张量中，导致动量方程具有非标准非线性；
• 流行病学部分的非线性感染率 $\beta(C)\frac{SI}{N}$，体现剂量响应关系；
• 所有物种均受速度场 $U$ 的对流作用，形成强耦合结构。

此外，函数空间的选择也体现了问题的复杂性。作者在分析中使用了以下空间：

• $H = \overline{W}^{L^2(\Omega)}$, 即 $L^2$ 中散度为零且满足无滑移边界的向量场闭包；
• $V = \overline{W}^{H^1_0(\Omega)}$, 对应能量空间；
• 各种 $L^p(0,T;X)$ 和 $H^1(\Omega)$、$W^{m,p}(\Omega)$ 等Sobolev空间。

特别地，弱解的正则性要求为：

• $U \in L^\infty(0,T;V) \cap L^2(0,T;H)$,
• $C,S,I,R \in L^\infty(0,T;L^2(\Omega)) \cap L^2(0,T;H^1(\Omega))$,
• 时间导数属于对偶空间，如 $\partial_t U \in L^2(0,T;H')$。

这些空间设定符合Navier-Stokes方程的经典框架，但由于粘度 $\nu(C)$ 的依赖性和多场强耦合，传统方法面临挑战。

本模型的应用背景涵盖自然灾害后的疫情爆发、城市供水系统污染、河流湖泊中的病原体扩散等问题。通过引入流体动力学，能够更真实地模拟病原体在复杂地形或流动环境中的传播路径，进而为公共卫生干预提供理论支持。

主要定理或主要结果

定理 3.1（弱解的存在性）

假设条件：

1. （A1）粘度函数 $\nu(C)$ 和传染率函数 $\beta(C)$ 连续，且存在正常数 $\nu_1,\nu_2,\beta_1,\beta_2 > 0$，使得
   $$ 0 < \nu_1 \leq \nu(s) \leq \nu_2, \quad 0 < \beta_1 \leq \beta(s) \leq \beta_2, \quad \forall s \in \mathbb{R}. $$
2. （A2）初始值满足：$U_0 \in L^2(\Omega)$, $(C_0,S_0,I_0,R_0) \in (L^2(\Omega))^4$，且几乎处处非负。
3. （A3）外力项 $f \in L^2(0,T;L^2(\Omega))$。

结论： 在上述条件下，SIRPNS系统（2.4）存在一个满足定义3.1的全局弱解 $(U,C,S,I,R)$，其正则性如前所述，且保持生物合理性（即各密度非负）。

定理 3.2（二维情形下的唯一性）

在额外假设下（未完全写出但隐含于分析中），当空间维度 $d=2$ 时，该弱解是唯一的。这得益于二维Navier-Stokes方程的良好性态以及Gagliardo-Nirenberg不等式的更强控制能力。

此外，论文还进行了数值模拟，验证了以下现象：

• 病原体浓度升高可显著改变流场结构（因粘度变化）；
• 流体对流能加速或引导疫情的空间蔓延；
• 存在“粘度-传播正反馈”机制：高病原体浓度 → 高粘度 → 局部滞留增强 → 更多感染 → 更高浓度。

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

a) 核心分析学工具

1. Faedo-Galerkin 方法
   这是处理非线性偏微分方程（尤其是Navier-Stokes类方程）存在性的经典方法。作者构造有限维逼近解序列 $(U_n, C_n, S_n, I_n, R_n)$，将其投影到由正交基张成的有限维子空间中，转化为常微分方程组求解局部解。具体而言，选取 $V \times H^1_0(\Omega)^4$ 的正交基 $\{(\xi_k, \ell_k)\}$，令近似解为：

   $$ U_n(t) = \sum_{k=1}^n u_k(t)\xi_k(x),\quad C_n(t) = \sum_{k=1}^n c_k(t)\ell_k(x),\quad \dots $$

   代入原方程后得到ODE系统 $\frac{d}{dt}U_n = F_n(U_n)$，利用Carathéodory定理保证局部可解性。

2. 紧性原理（Aubin-Lions引理）
   为了通过极限过程 $n \to \infty$ 获得真解，必须证明逼近序列在适当拓扑下相对紧。作者多次应用Aubin-Lions紧嵌入引理（Lemma 3.1），其形式为： 若 $X_0 \Subset X \hookrightarrow X_1$，且函数族 $G$ 在 $L^p(0,T;X_0)$ 有界，时间导数 $\partial_t G$ 在 $L^1(0,T;X_1)$ 有界，则 $G$ 在 $L^p(0,T;X)$ 相对紧。

   特别地，用于控制 $U_n$ 的时间导数 $\partial_t U_n \in L^2(0,T;H')$，结合空间正则性 $U_n \in L^2(0,T;H)$，可得 $U_n$ 在 $L^2(Q_T)$ 中强收敛。

3. 能量估计与先验估计
   对每个方程进行测试函数法（test function method），建立统一的能量框架。例如，对Navier-Stokes部分取 $Z = U_n$，利用不可压性和边界条件消去非线性项 $(U_n \cdot \nabla)U_n \cdot U_n$，得到：

   $$ \frac{1}{2}\frac{d}{dt}\|U_n\|^2_{L^2} + \|\sqrt{\nu(C_n)}\nabla U_n\|^2_{L^2} \leq \|f\|_{L^2}\|U_n\|_{L^2}, $$

   再结合 $\nu(C_n) \geq \nu_1 > 0$，可得 $U_n$ 在 $L^\infty(0,T;L^2) \cap L^2(0,T;H^1)$ 一致有界。

   类似地，对 $C_n, S_n, I_n, R_n$ 使用各自的测试函数（如 $\psi_C = C_n$），并利用最大值原理思想（见Lemma 3.3），通过取负部 $u_-$ 作为测试函数，证明非负性保持。

4. Gagliardo-Nirenberg 不等式与Young不等式
   文中引用了如下嵌入不等式：

   $$ \|Z(t)\|_{L^4(\Omega)} \leq c \|Z(t)\|^\zeta_{H^1(\Omega)} \|Z(t)\|^{1-\zeta}_{L^2(\Omega)},\quad \zeta = d/4, $$

   在估计非线性项（如 $(U \cdot \nabla)U$ 或 $\beta(C)SI/N$）时至关重要。配合Young不等式：

   $$ ab \leq \delta a^p + c(\delta)b^q, $$

   可吸收高阶项至左侧，完成闭合估计。

5. 单调算子与弱收敛技术
   在极限过程中，需处理非线性项如 $\nu(C_n)\nabla U_n$。由于 $\nu(C_n)$ 弱收敛、$\nabla U_n$ 弱收敛，乘积未必收敛。为此，作者借助紧性得出 $C_n \to C$ 几乎处处，故 $\nu(C_n) \to \nu(C)$ 几乎处处，再结合有界性，利用Lebesgue控制收敛定理获得强收敛，从而保证非线性项的极限正确传递。

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

1. 粘度非线性项 $\mathrm{div}(\nu(C)\nabla U)$ 的处理
   此项打破了Navier-Stokes方程的标准结构，使得经典变分形式不再适用。作者的关键技巧是：利用 $C_n$ 的紧性（来自Aubin-Lions）获得 $C_n \to C$ 强收敛于 $L^2(Q_T)$，从而 $\nu(C_n) \to \nu(C)$ 几乎处处且一致有界，因此 $\nu(C_n)\nabla U_n \rightharpoonup \nu(C)\nabla U$ 在 $L^2$ 中成立。

2. 保持生物合理性的非负性证明
   为避免近似解出现负值（违反生物学意义），作者采用“截断+测试负部”的技巧。例如，在方程中将 $\beta(C)\frac{SI}{N}$ 替换为 $\beta(C)\frac{S_+ I_+}{C_+ + S_+ + I_+}$，然后选择 $\psi_S = -S_{n,-}$ 作为测试函数，推导出：

   $$ \frac{1}{2}\frac{d}{dt}\int_\Omega |S_{n,-}|^2 dx + D_S \int_\Omega |\nabla S_{n,-}|^2 dx + \eta \int_\Omega |S_{n,-}|^2 dx = -\int_\Omega \Lambda S_{n,-} dx - \cdots \leq 0, $$

   初始非负即推出 $S_{n,-}(t) = 0$ 对所有 $t$ 成立，即 $S_n \geq 0$。

3. 对流项 $(U \cdot \nabla)U$ 的控制
   在能量估计中，该项看似会产生增长项，但实际上由于不可压性 $\nabla \cdot U = 0$，有：

   $$ \int_\Omega (U \cdot \nabla)U \cdot U dx = -\frac{1}{2}\int_\Omega |U|^2 (\nabla \cdot U) dx = 0, $$

   故可被完美消除，这是Navier-Stokes分析中的经典技巧。

4. 分数型非线性项 $\frac{SI}{N}$ 的处理
   当 $N \to 0$ 时该项可能奇异，但作者假设总人口有正下界或通过正则化处理。在实际估计中，利用 $\beta(C) \leq \beta_2$ 和Hölder不等式控制：

   $$ \left|\int \beta(C)\frac{SI}{N} \phi dx\right| \leq \beta_2 \|\phi\|_{L^\infty} \|S\|_{L^2} \|I\|_{L^2}, $$

   并结合Sobolev嵌入控制测试函数。

与之前工作的比较

a) 突破与改进

本文相较于以往研究实现了多项实质性突破：

1. 首次实现SIR模型与Navier-Stokes方程的全耦合建模
   已有文献（如 [6–9, 11, 13, 16–18, 20, 21, 24–26, 28–32]）大多将病原体视为被动扩散-反应方程中的浓度场，忽略其对流效应；即使考虑流体输运，也常假定速度场预先给定（如恒定流或Stokes流）。本文则将速度场 $U$ 本身作为未知函数，通过Navier-Stokes方程动态演化，真正实现宿主-病原体-流体三重耦合。

2. 引入病原体浓度依赖的粘度函数 $\nu(C)$，建立双向反馈机制
   这是本文最突出的创新点。传统模型认为流体参数恒定，而现实中微生物增殖、有机物释放会显著增加液体粘稠度（如黏液环境、生物膜形成）。本文明确提出 $\nu = \nu(C)$，使流行病过程反过来影响流体行为，形成闭环反馈。这种非标准非线性粘度项此前在流行病模型中从未出现。

3. 数学分析上的严格性超越同类模型
   多数环境流行病模型只做数值模拟或线性稳定性分析，缺乏全局存在性证明。本文采用Faedo-Galerkin+紧性方法，严格证明了全局弱解的存在性，并在二维情形下获得唯一性，达到了Navier-Stokes类方程的现代PDE分析水平。

4. 放宽对病原体输运机制的简化假设
   以往工作常假设病原体仅靠扩散传播，或置于预设流场中被动输运。本文明确指出：“当病原体积累改变流体特性时，现有模型无法捕捉此类反馈”。通过耦合NS方程，自然包含了对流主导、涡旋捕获、边界层效应等复杂输运机制。

5. 更贴近现实的生物机制建模

   • 引入剂量响应函数 $\beta(C)$，而非常数传染率；
   • 考虑感染者排菌率 $\alpha I$；
   • 包含出生/死亡项 $\Lambda, \eta$，允许人口动态变化；
   • 所有宿主类别均可移动（扩散+对流），反映人类行为。

6. 维度适应性更强
   虽然唯一性仅证得二维情形（受限于NS方程本身），但存在性适用于 $d=2,3$，更具普适性。相比之下，许多反应扩散流行病模型仅在低维或周期域中可解。

综上所述，本文不仅在建模层面提出了前所未有的耦合结构，在数学理论上也达到了较高严谨程度，填补了“流体力学-传染病学”交叉领域的空白，为后续研究提供了坚实的理论基础和数值范式。

7. Attractors and their dimensions for the 3D Fractional Navier–Stokes–Voigt Equations

三维分数阶Navier–Stokes–Voigt方程的吸引子及其维数

作者: Alexei Ilyin, Varga Kalantarov, Sergey Zelik

类别: math.AP, 35Q30, 35B41, 37L30

发布日期: 2025-11-06

arXiv ID: 2511.04783v1

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

详细分析

中文标题: 三维分数阶Navier–Stokes–Voigt方程的吸引子及其维数

研究对象和背景

本文研究的是三维分数阶Navier–Stokes–Voigt方程（fractional Navier–Stokes–Voigt equations, fNSV），其数学形式如下：

其中 $ u = u(x,t) = (u_1,u_2,u_3) $ 是速度场，$ p = p(t,x) $ 是压力函数，$ \nu > 0 $ 为运动粘性系数，$ h \in H^{-1}_\sigma(\Omega) $ 是外力项，$ \alpha > 0 $ 是正则化参数，$ s \in [\frac{1}{2},1] $ 表示分数阶Stokes算子 $ A^s $ 的阶数。非线性项定义为：

该系统在有界光滑区域 $ \Omega \subset \mathbb{R}^3 $ 上考虑，边界条件为齐次Dirichlet边界条件。

这个系统的显著特征在于时间导数前引入了带有分数幂的正则化项 $ \alpha^s A^s \partial_t u $，这使得它成为经典Navier–Stokes方程的一种自然推广与Voigt型正则化。当 $ s = 1 $ 时，方程退化为经典的Navier–Stokes–Voigt（NSV）系统，最早由O. Ladyzhenskaya于1966年提出作为对三维Navier–Stokes方程的数值稳定化手段。后来发现这类模型也适用于描述水溶性聚合物溶液的流动行为，并被广泛用于流体力学中的数值模拟中，尤其是在小 $ \alpha $ 情况下能有效抑制高频振荡。

从PDE结构上看，fNSV系统通过添加高阶时间导数项增强了耗散机制，从而避免了经典Navier–Stokes方程中可能出现的解的爆破问题。这种正则化方式不仅提升了方程的适定性，还使其长期动力学行为更易于分析。事实上，在 $ s \geq \frac{1}{2} $ 的条件下，可以证明全局弱解的存在唯一性和强连续性。

在函数空间方面，作者使用了一系列基于Stokes算子 $ A $ 的Hilbert空间族 $ H^s_A $，这些空间由 $ A $ 的谱分解构造而来。设 $ \{e_n\}_{n=1}^\infty $ 为 $ A $ 在空间 $ H $（即 $ L^2 $ 中的无散向量场且满足零边界条件的空间）上的标准正交特征基，对应特征值为 $ \{\lambda_n\} $，则任意函数 $ u \in H $ 可表示为：

并由此定义Sobolev型空间：

特别地，$ V_s := H^s_A $ 被用作相空间，而 $ A^s $ 定义为：

值得注意的是，这里的分数幂是谱意义下的（spectral fractional Stokes operator），而非区域外或非局部意义上的分数拉普拉斯。这一选择保证了 $ A^s $ 与 $ A $ 共享相同的特征基，极大简化了后续关于吸引子维度估计的技术处理。

此外，文中强调了负索博列夫空间 $ H^{-s}_\sigma(\Omega) $ 的微妙性质：当 $ s > \frac{1}{2} $ 时，$ H^{-s}_\sigma $ 并不能嵌入到全空间分布 $ D'(\Omega) $ 中，这意味着传统的分布理论框架需谨慎使用。为此，作者采用Leray–Helmholtz投影 $ \Pi $ 将 $ L^2 $ 空间正交分解为无旋与梯度部分，并借助此投影定义负阶空间中的“投影”作用，确保变分形式的合理性。

综上所述，本文的研究对象是一个在物理上有明确背景、数学上具有良好适定性的偏微分方程系统，其核心挑战在于刻画其长期动力学行为——特别是全局吸引子的分形维度上界，而这正是本文的主要贡献所在。

主要定理或主要结果

本文的核心成果体现在以下定理中，给出了全局吸引子分形维度（fractal dimension）的精确上界估计，依赖于关键物理参数如Grashof数 $ G $ 和正则化参数相关的 $ G_1 $。

定理 1.1（主要结果）

设 $ \Omega \subset \mathbb{R}^3 $ 为光滑有界区域，$ h \in L^2(\Omega) $。令 $ \mathcal{A} $ 为问题 (1.1) 的全局吸引子。则：

• 当 $ s \in [\frac{3}{4},1] $ 时，有

  $$ \dim_f(\mathcal{A}) \leq C \left( \frac{1 + G_1^s}{G_1^s} \right)^{3/2} G^{3/2}, $$

• 当 $ s \in [\frac{1}{2}, \frac{3}{4}] $ 时，有

  $$ \dim_f(\mathcal{A}) \leq C \left( \frac{1 + G_1^s}{G_1^s} \right)^{3/2} G^{6(1-s)}, $$

• 特别地，当 $ s = 1 $ 且 $ G_1^3 G^4 \ll 1 $ 时，

  $$ \dim_f(\mathcal{A}) \leq C (G G_1^3)^{6/13}. $$

其中常数 $ C $ 为绝对常数，$ G $ 为三维Grashof数：

而 $ G_1 $ 是一个新引入的无量纲参数：

这里 $ \lambda_1 $ 是Stokes算子的第一特征值。

上述估计在 $ \alpha \lambda_1 $ 充分小时成立。作者指出，即使在经典情形 $ s=1 $ 下，所得结果也显著优于此前文献中的估计，例如比[26]中的

更为精细。

此外，文中说明目前尚无吸引子的下界估计，因此无法判断这些上界是否最优。但在极限 $ G_1 \to \infty $ 情形下，推测下界应为 $ O(G) $，而当前上界为 $ O(G^{3/2}) $，仍存在差距。

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

本文在估计吸引子维度时采用了现代无穷维动力系统理论的标准框架：首先建立系统的耗散性与紧性，从而得到全局吸引子的存在性；然后利用Lyapunov指数求和法结合迹公式来控制吸引子的维度。具体而言，维度估计的关键步骤依赖于对线性化系统的横截扩张体积的增长率进行控制，最终归结为对某些谱不等式的应用。

a) 核心分析工具

1. Lieb–Thirring不等式

   这是本文最重要的工具之一。经典的Lieb–Thirring不等式提供了关于Schrodinger算子负特征值之和的上界，而在本文中被用于估计向量场梯度的积分与特征值分布之间的关系。具体形式如下：

   $$ \sum_{j=1}^N \int_\Omega |\nabla \phi_j(x)|^2 dx \geq C |\Omega|^{-2/3} N^{5/3}, $$

   其中 $ \{\phi_j\} $ 是某个子空间的标准正交基。该不等式反映了特征函数的“空间震荡频率”随模态数量增加的趋势，是连接几何与谱信息的桥梁。

2. Cwikel–Lieb–Rosenblum（CLR）不等式

   用于估计负特征值个数的CLR不等式在本文中被推广至分数阶设置。对于算子 $ -\Delta + V(x) $，其负特征值个数满足：

   $$ N(-\Delta + V) \leq C \int_\Omega (V_-(x))^{3/2} dx, $$

   其中 $ V_- = \max\{0, -V\} $。作者将其与分数阶算子结合，发展出适用于 $ A^s $ 的版本，用于控制扰动后的谱分布。

3. Berezin–Li–Yau型不等式

   用于估计特征值累加和的下界，形式为：

   $$ \sum_{j=1}^N \lambda_j \geq C |\Omega|^{-2/3} N^{5/3}. $$

   这类不等式帮助作者建立了从谱数据到空间维度的定量联系。

4. 分数阶Sobolev嵌入与插值不等式

   利用 $ V_s \hookrightarrow L^{6/(3-2s)} $ 和 $ V_1 \hookrightarrow L^6 $ 等嵌入关系，控制非线性项 $ (u,\nabla)u $ 的可积性。例如，在能量估计中需要验证：

   $$ (u,\nabla)u \in L^2((0,T); V_{s - 3/2}), $$

   这依赖于 $ s \geq \frac{1}{2} $ 才能保证 $ V_{s - 3/2} \subset V_{-1} $。

5. 迹公式与Lyapunov指数求和

   吸引子维度的经典估计方法基于线性化流的迹：

   $$ \dim_f(\mathcal{A}) \leq 1 + \sup_{u_0 \in \mathcal{A}} \sup_{t > 0} \frac{1}{t} \int_0^t \sum_{j=1}^N \mu_j(u(s)) ds, $$

   其中 $ \mu_j $ 是线性化算子的第 $ j $ 个Lyapunov指数。通过对Jacobi方程的分析，将此和转化为对算子迹的估计，进而应用上述谱不等式。

b) 技术难点与处理技巧

1. 处理双参数依赖性

   与经典Navier–Stokes方程仅依赖Grashof数 $ G $ 不同，NSV系统引入了第二个无量纲参数 $ G_1 = \alpha \lambda_1 $，导致维度估计必须同时考虑 $ G $ 与 $ G_1 $ 的相对大小。作者通过细致划分参数区间（如 $ G_1 \ll 1 $ 或 $ G_1 \gg 1 $）分别优化估计，实现了对不同物理场景的统一覆盖。

2. 分数阶算子带来的技术复杂性

   分数阶 $ A^s $ 导致能量空间 $ V_s $ 的正则性介于整数阶之间，使得传统的整数阶估计不再适用。作者通过引入广义Gagliardo–Nirenberg不等式和谱投影分解，将高阶导数控制转化为特征值求和问题，再利用前述谱不等式完成估计。

3. 非线性项的精细估计

   非线性项 $ ((v,\nabla)u_2,v) $ 在唯一性证明中出现，需用Hölder不等式与Sobolev嵌入估计：

   $$ |((v,\nabla)u_2,v)| \leq \|v\|_{L^3} \|\nabla u_2\|_{L^2} \|v\|_{L^6} \leq C \|\nabla u_2\|_{L^2} \|v\|_{V_{1/2}}^2. $$

   再利用 $ V_{1/2} \subset V_s $（因 $ s \geq 1/2 $），将其吸收进主耗散项。

4. 能量恒等式的严格推导

   文中通过测试函数取为解本身，并利用 $ (1+\alpha^s A^s)^{1/2} u \in L^2(V_{1-s}) $ 和分布导数的延拓技术，严格证明了能量恒等式：

   $$ \frac{d}{dt} \left( \|u(t)\|^2 + \alpha^s \|u(t)\|^2_{V_s} \right) + 2\nu \|\nabla u(t)\|^2 = 2\langle h, u(t) \rangle, $$

   这是建立耗散性与吸引子存在的基础。

与之前工作的比较

本文相较于已有工作做出了若干重要改进：

1. 推广至分数阶情形

   此前关于Navier–Stokes–Voigt系统的吸引子维度研究主要集中于 $ s=1 $ 的整数阶情况（如[26],[31]）。本文首次系统研究了 $ s \in [\frac{1}{2},1] $ 的分数阶版本，填补了该领域的空白。

2. 显著改进经典情形下的估计精度

   即使在 $ s=1 $ 的经典设定下，本文给出的估计也远优于先前结果。例如，相比[26]中依赖于 $ (\alpha \lambda_1)^{-3} $ 的立方衰减，本文通过更精细的谱分析得到了更优的 $ (G G_1^3)^{6/13} $ 型估计，尤其在 $ \alpha \lambda_1 \to 0 $ 极限下表现更好。

3. 引入先进的谱不等式工具

   以往研究多依赖于Weyl渐近律或粗糙的特征值计数，而本文系统运用了Lieb–Thirring、CLR、Berezin–Li–Yau等现代谱不等式组合，实现了从几何到谱再到动力学维度的精确传递。

4. 统一处理多种参数极限

   旧有估计往往只适用于特定参数范围（如 $ \alpha \to 0 $ 或 $ \alpha \to \infty $），而本文通过分段估计策略，能够统一处理 $ G_1 \to 0 $、$ G_1 \to \infty $ 等多种极限情形，增强了结果的普适性。

5. 为未来下界估计提供方向

   虽然目前缺乏吸引子维度的下界，但作者指出可通过构造类似“Vishik涡旋”的特殊解来逼近下界，提示未来可能实现上下界的匹配。

综上，本文不仅在技术深度上超越前人，也在物理适用性和数学严谨性之间取得了良好平衡，代表了该领域当前的前沿水平。

8. Finite time blow-up for a multi-dimensional model of the Kiselev-Sarsam equation

一个关于Kiselev-Silvestre方程的多维模型的有限时间爆破

作者: Wanwan Zhang

类别: math.AP

发布日期: 2025-11-06

arXiv ID: 2511.04660v1

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

详细分析

中文标题: 一个关于Kiselev-Silvestre方程的多维模型的有限时间爆破

研究对象和背景

本文研究了一类多维非局部活性标量方程，其形式为：

其中 $n \geq 2$ 是空间维度，常数 $g > 0$ 表示重力加速度因子，而 $\mathcal{R}_a$ 是一个非局部奇异积分算子，定义如下：

其中 $a > 0$ 是一个固定参数，$\Gamma(\cdot)$ 是标准的 Gamma 函数，P.V. 表示柯西主值积分。该算子 $\mathcal{R}_a$ 可以看作是经典 Riesz 变换与一种“正则化”版本之间的差。

这个模型被称为 Kiselev-Sarsam 方程的多维推广。原始的一维 Kiselev-Sarsam 方程最早由 Kiselev 和 Sarsam 在文献 [14] 中提出，作为二维不可压缩多孔介质（IPM）方程的一个简化一维类比模型。IPM 方程描述的是在达西定律和重力作用下，由不可压缩流体输运的密度场 $\rho(x,t)$ 的演化，其形式为：

其中 $P$ 是压力项。由于 IPM 方程本身的数学分析极为困难，特别是关于解是否会在有限时间内形成奇性的开放问题，因此构造并研究其简化模型成为理解此类非线性非局部现象的重要途径。

当空间维度 $n=1$ 时，上述多维模型退化为著名的 Kiselev-Sarsam 方程：

其中 $H_a$ 是一个非局部算子，称为带参数的希尔伯特型变换，其核函数为：

该算子 $H_a$ 也可以被视作经典的 Hilbert 变换 $H$ 与一个光滑核（类似于共轭 Poisson 核）之差。特别地，当 $a \to 0$ 时，$H_a \to 0$；当 $a \to \infty$ 时，$H_a \to H$。因此，它提供了一个从平凡算子到典型奇异积分算子的连续插值。

与此密切相关的还有著名的 Córdoba-Córdoba-Fontelos (CCF) 方程：

该方程也源于对欧拉方程或 SQG 方程的建模。CCF 方程已被证明对于某些光滑偶函数初值会发生梯度爆破（即 $\|\nabla \theta(t)\|_{L^\infty} \to \infty$ 在有限时间内发生），这主要得益于作者在 [4] 中发现的一个关键不等式：

适用于满足一定对称性和单调性条件的函数 $f$。

类似地，高维情形下的 CCF 推广为：

其中 $R\theta = (R^{(1)}\theta, \dots, R^{(n)}\theta)$ 是 $\theta$ 的 Riesz 变换。Balodis 和 Córdoba 在 [2] 中证明了该方程的局部适定性以及对径向初值的有限时间爆破。

本文提出的多维 Kiselev-Sarsam 模型 $\partial_t \rho + g \mathcal{R}_a \rho \cdot \nabla \rho = 0$ 正是在这一系列工作基础上的自然推广：它将一维的 $H_a$ 算子推广到了高维的 $\mathcal{R}_a$ 算子。值得注意的是，$\mathcal{R}_a$ 同样具有插值性质：当 $a \to 0$ 时，$\mathcal{R}_a \to 0$；当 $a \to \infty$ 时，$\mathcal{R}_a \to R$，即趋向于经典的 Riesz 变换。因此，该模型统一地连接了零算子、正则化奇异积分与完全奇异积分三类情形。

文章所考虑的空间区域是整个欧几里得空间 $\mathbb{R}^n$，没有边界条件限制。函数空间方面，初始数据 $\rho_0$ 被假设属于 Sobolev 空间 $H^s(\mathbb{R}^n)$，其中 $s > \frac{n}{2} + 1$，这是保证局部适定性所需的基本正则性条件。此外，在研究爆破问题时，特别关注一类具有径向对称性且单调不减的光滑紧支撑初值，并进一步要求其满足一个加权 $L^2$ 型条件（见后文 Theorem 1.2）。这些结构假设对于建立关键的非线性项下界估计至关重要。

主要定理或主要结果

本文建立了两个核心结果：局部适定性与有限时间梯度爆破。

定理 1.1（局部适定性）
设 $n \geq 2$，$a,g > 0$。对于任意初值 $\rho_0 \in H^s(\mathbb{R}^n)$，其中 $s > \frac{n}{2} + 1$，存在依赖于 $\|\rho_0\|_{H^s}$ 的时间 $T = T(\|\rho_0\|_{H^s}) > 0$，使得方程 (1.1) 在区间 $[0,T)$ 上存在唯一解 $\rho \in C([0,T); H^s(\mathbb{R}^n)) \cap \mathrm{Lip}((0,T); H^{s-1}(\mathbb{R}^n))$。

该定理表明，在足够高的 Sobolev 正则性下，该非局部方程存在唯一的局部经典解。

定理 1.2（有限时间梯度爆破）
设 $n \geq 2$，$a,g > 0$，$\delta \in (0,1)$。假设初值 $\rho_0: \mathbb{R}^n \to \mathbb{R}$ 是光滑的、径向对称的、单调不减的且具有紧支撑。此外，假设其满足以下加权积分条件：

其中 $C_\delta > 0$ 是仅依赖于 $\delta$ 的常数。则对应于该初值的解 $\rho$ 将在有限时间内发生梯度爆破，即存在有限时间 $T^* < \infty$，使得

该定理是本文的核心成果，首次证明了多维 Kiselev-Sarsam 模型中存在导致有限时间奇性形成的初值。

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

a) 核心分析学工具

论文在证明过程中综合运用了多种现代偏微分方程分析技术。

首先，在证明局部适定性（Theorem 1.1）时，采用了经典的 能量方法 结合 Sobolev 空间理论。具体而言，通过建立一系列先验估计来控制解的 Sobolev 范数增长。关键步骤包括：

• Grönwall 不等式：用于从微分不等式推导出解范数的指数型增长上界。
• Kato-Ponce 交换子估计：这是处理非线性项 $\mathcal{R}_a \rho \cdot \nabla \rho$ 高阶导数时的核心工具。引理 3.3 给出了如下形式的估计：
  $$ \|\Lambda^s(fg) - f \Lambda^s g\|_{L^p} \leq C \left( \|\Lambda^s f\|_{L^{p_1}} \|g\|_{L^{p_2}} + \|\nabla f\|_{L^{p_3}} \|\Lambda^{s-1} g\|_{L^{p_4}} \right), $$
  其中 $\Lambda^s = (-\Delta)^{s/2}$ 是分数阶拉普拉斯算子。此估计允许作者将高阶导数作用于乘积项时产生的“误差”（即交换子）用低阶范数控制，从而封闭能量估计。
• BMO 空间与对数不等式：为了建立 Beale-Kato-Majda (BKM) 型爆破准则（Proposition 4.2），作者引用了一个重要的极限 Sobolev 不等式（引理 4.1）：
  $$ \|f\|_{L^\infty} \leq C (1 + \|f\|_{\mathrm{BMO}})(1 + \log^+ \|f\|_{W^{s,p}}) $$
  该不等式将 $L^\infty$ 范数与 BMO 范数及高阶 Sobolev 范数联系起来，结合 $\mathcal{R}_a$ 将 $L^\infty$ 映射到 BMO 的有界性（$\|\mathcal{R}_a f\|_{\mathrm{BMO}} \leq C \|f\|_{L^\infty}$），从而得到梯度范数与 Sobolev 范数之间的对数型控制关系。

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

在证明有限时间爆破时，论文借鉴了 CCF 类方程的研究思想，但面临高维和非完全奇异算子的新挑战。其核心技术路径如下：

1. 保持对称性：利用方程的结构和解的唯一性，证明若初值是径向且单调不减的，则解在其生命周期内始终保持这一性质（引理 4.4）。这一点至关重要，因为它允许后续估计在径向坐标下进行，极大简化了计算。

2. 最大值原理：证明 $L^\infty$ 范数守恒（引理 4.3），即 $\|\rho(t)\|_{L^\infty} = \|\rho_0\|_{L^\infty}$。这为后续加权估计提供了稳定的上界。

3. 构造加权能量泛函：受一维情形中加权 Hilbert 变换不等式的启发，作者考虑形如 $\int \frac{(\rho(x)-\rho(0))^2}{|x|^{n+\delta}} dx$ 的泛函。通过对该泛函求导并分析其演化，目标是证明其导数有一个正的下界。

4. 非线性项的下界估计：这是最核心的技术难点。作者需要估计如下形式的积分：

   $$ -\int_{\mathbb{R}^n} \frac{(\mathcal{R}_a \rho \cdot \nabla \rho)(x)}{|x|^{n+\delta}} (\rho(x) - \rho(0)) \, dx. $$

   利用 $\mathcal{R}_a$ 的积分核表达式和解的径向单调性，通过细致的积分估计（包括主值积分的处理和核函数的点态符号分析），证明该非线性项可以被一个正的加权 $L^2$ 范数所控制，即存在常数 $C > 0$ 使得

   $$ -\int \cdots \geq C \int \frac{(\rho(x) - \rho(0))^2}{|x|^{n+1+\delta}} dx. $$

5. 微分不等式与爆破论证：结合上述下界估计和最大值原理，最终得到关于加权能量 $E(t)$ 的微分不等式：

   $$ \frac{d}{dt} E(t) \geq c E(t)^2, $$

   其中 $c > 0$。一旦初始能量 $E(0)$ 足够大（即满足定理中的条件），该不等式立即推出 $E(t)$ 在有限时间内趋于无穷，从而迫使 $\|\nabla \rho(t)\|_{L^\infty}$ 爆破。

与之前工作的比较

a) 突破与改进

本文的主要贡献在于将一维的 Kiselev-Sarsam 模型成功推广至高维，并首次证明了其存在有限时间爆破现象。相比已有工作，其突破体现在以下几个方面：

1. 维度的实质性推广：此前关于 Kiselev-Sarsam 方程的研究（如 [14], [19]）均局限于一维情形。本文构建了合理的高维算子 $\mathcal{R}_a$，使其既保持原算子的傅里叶乘子结构（$\widehat{\mathcal{R}_a f}(\xi) = -i \frac{\xi}{|\xi|}(1-e^{-2\pi a |\xi|}) \hat{f}(\xi)$），又具备良好的奇异积分性质。这种构造本身就是一个重要进展。

2. 填补研究空白：尽管高维 CCF 方程（对应 $a \to \infty$ 极限）已被 Balodis-Córdoba 和 Dong 等人研究过，但带有正则化参数 $a>0$ 的中间情形——即本文的 $\mathcal{R}_a$ 模型——在高维的爆破问题尚属空白。本文首次系统研究了这类“部分奇异”非局部输运方程的奇性形成机制。

3. 方法的适应性与创新：虽然证明思路借鉴了 CCF 方程中利用加权不等式的方法，但由于 $\mathcal{R}_a$ 不再是齐次的 Calderón-Zygmund 算子，传统的点态不等式（如 Li-Rodrigo 对 CCF 的证明）不再直接适用。作者必须重新分析 $\mathcal{R}_a$ 的核函数行为，并结合 $a>0$ 带来的额外衰减特性，发展出适用于此类正则化奇异积分的新估计技术。

4. 物理意义的延伸：由于原始 Kiselev-Sarsam 模型是对二维 IPM 的一维模拟，本文的高维模型可能为理解更高维多孔介质流动或其他非局部输运过程提供新的数学视角，尽管其直接物理对应尚待探索。

综上所述，本文不仅在数学上完成了从一维到高维的关键跨越，而且在非线性非局部 PDE 的奇性形成理论中，为一类新的插值型模型建立了完整的局部理论与爆破结果，显著拓展了该领域的研究范围。

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

1. Asymptotic Stability of Solutions to the Forced Higher-order Degenerate Parabolic Equations

强迫高阶退化抛物方程解的渐近稳定性

作者: Jinhong Zhao, Bin Guo

类别: math.AP, cs.NA, math-ph, math.MP, math.NA, 76A05, 76A20, 35B40, 35Q35, 35K35

发布日期: 2025-11-07

arXiv ID: 2511.04947v1

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

摘要翻译

中文标题: 强迫高阶退化抛物方程解的渐近稳定性

摘要翻译: 本文研究一类四阶拟线性退化抛物型方程，在含时与不含时非齐次外力作用下的解的行为，该方程用于模拟“完全润湿”情形下非牛顿流体在固体表面的薄膜流动。利用高阶抛物型方程的正则性理论和能量方法，我们建立了正弱解的整体存在性，并刻画了解的长时间行为。

具体而言，对于受含时外力 $f(t,x)$ 作用的幂律型薄膜问题，我们证明了弱解收敛到 $\bar{u}_0 + \frac{1}{|\Omega|}\int_{0}^t \int_{\Omega} f(s,x) \, {\rm d}x \, {\rm d}s$，并给出了收敛速率，其中 $\bar{u}_0$ 表示初始数据的空间平均值。与文献 \cite{JJCLKN}（Jansen 等，2023）中齐次情形的结果相比，该结果明确揭示了非齐次外力对解的收敛速率的影响。对于不依赖时间的外力 $f(x)$，我们证明了弱解与线性函数 $\bar{u}_0 + \frac{t}{|\Omega|}\int_\Omega f(x)\, {\rm d}x$ 之差是一致有界的。当外力为常数 $f_0$ 时，我们进一步表明：在剪切增稠情形下，弱解会在有限时间内精确等于 $\bar{u}_0 + t f_0$；而在剪切减薄和牛顿流体情形下，弱解分别以多项式速率和指数速率趋近于 $\bar{u}_0 + t f_0$。随后，针对Ellis律描述的薄膜问题，我们发现其解的行为类似于牛顿流体的解。

最后，我们通过数值模拟验证了主要的理论分析结果。

不相关论文（基本信息）

1. The higher-order fractional Schrödinger equation with nonlinear local perturbations: Uniqueness

具有非线性局部扰动的高阶分数阶薛定谔方程：唯一性

作者: Giovanni Covi, Ru-Yu Lai, Lili Yan

类别: math.AP, 35R30, 35Q55

发布日期: 2025-11-07

arXiv ID: 2511.05384v1

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

摘要: We study the higher-order fractional Schr"odinger equation with local nonlinear perturbations and investigate both the forward and inverse problems. We establish both the Sobolev $H^s$ and H"older $C^s$ estimates for the well-posedness of the nonlinear problem, based on the corresponding estimates derived for the linear fractional Schr"odinger equation. For the inverse problem, we show that the local nonlinear perturbations can be uniquely determined from the Dirichlet-to-Neumann map, by using the higher-order linearization and the unique continuation property of the fractional Laplace operator.

2. Balanced Dynamics in Strongly Coupled Networks

强耦合网络中的平衡动力学

作者: Cristobal Quininao, Jonathan Touboul

类别: math.AP, nlin.AO

发布日期: 2025-01-20

arXiv ID: 2501.11769v2

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

摘要: Many mathematical models of interacting agents assume that individual interactions scale down in proportion to the network size, ensuring that the combined input received from the network does not diverge. In theoretical neuroscience, Sompolinsky and Van Vreeswijk proposed in 1996 that, should these scalings be violated (and under appropriate conditions), the system may not diverge but rather approach a balanced state where the inputs to each neuron compensate each other (in neuroscience, where inhibitory currents compensate the excitatory ones). We come back to this observation and formulate here a mathematical conjecture for the occurrence of such behaviors in general stochastic systems of interacting agents. From a mathematical viewpoint, this conjecture can be viewed as a double-limit problem in the space of probability measures, which we discuss in detail, as it provides several possible mathematical avenues for proving this result. We provide some numerical and theoretical explorations of the conjecture in classical models of neuronal networks. Moreover, we provide a complete proof of an asymptotic result consistent with one of the double-limit problems in a one-dimensional model with separable coupling inspired by models of chemically-coupled neurons. This proof relies on asymptotic methods, and particularly desingularization techniques used in some PDEs, that we apply here to the mean-field limit of the network as the coupling is made to diverge. From the applications viewpoint, this theory provides an alternative, minimalistic explanation for the widely observed balance of excitation and inhibition in the cerebral cortex not requiring the assumption of the existence of complex regulatory mechanisms.

3. A time-fractional Fisher-KPP equation for tumor growth: Analysis and numerical simulation

用于肿瘤生长的分数阶Fisher-KPP方程：分析与数值模拟

作者: Marvin Fritz, Nikos I. Kavallaris

类别: math.AP, cs.NA, math.NA, 35A01, 35R11, 65M12, 65M60, 92C50

发布日期: 2025-11-07

arXiv ID: 2511.05312v1

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

摘要: We study a time-fractional Fisher-KPP equation involving a Riemann-Liouville fractional derivative acting on the diffusion term, as derived by Angstmann and Henry (Entropy, 22:1035, 2020). The model captures memory effects in diffusive population dynamics and serves as a framework for tumor growth modeling. We first establish local well-posedness of weak solutions. The analysis combines a Galerkin approximation with a refined a priori estimate based on a Bihari-Henry-Gronwall inequality, addressing the nonlinear coupling between the fractional diffusion and the reaction term. For small initial data, we further prove global well-posedness and asymptotic stability. A numerical method based on a nonuniform convolution quadrature scheme is then proposed and validated. Simulations demonstrate distinct dynamical behaviors compared to conventional formulations, emphasizing the physical consistency of the present model in describing tumor progression.

4. The weighted isoperimetric inequality and Sobolev inequality outside convex sets

凸集外部的加权等周不等式与索博列夫不等式

作者: Lu Chen, Jiali Lan

类别: math.AP

发布日期: 2025-10-02

arXiv ID: 2510.01647v2

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

摘要: In this paper, we establish a weighted capillary isoperimetric inequality outside convex sets using the $\lambda_w$-ABP method. The weight function $w$ is assumed to be positive, even, and homogeneous of degree $\alpha$, such that $w^{1/\alpha}$ is concave on $\R^n$. Based on the weighted isoperimetric inequality, we develop a technique of capillary Schwarz symmetrization outside convex sets, and establish a weighted P'{o}lya-Szeg"{o} principle and a sharp weighted capillary Sobolev inequality outside convex domain. Our result can be seen as an extension of the weighted Sobolev inequality in the half-space established by Ciraolo-Figalli-Roncoroni in \cite{CFR}.

5. On partial diffusion and mixing without hypoellipticity

论部分扩散与非 hypoellipticity 条件下的混合

作者: Xu’an Dou, Delphine Salort, Didier Smets

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.05280v1

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

摘要: A simple Markov process is considered involving a diffusion in one direction and a transport in a transverse direction. Quantitative mixing rate estimates are obtained with limited assumptions about the transport field, which might be highly irregular and/or highly degenerate, in particular quite far from satisfying an hypoellipticity type assumption.

6. Regularized Reconstruction of Scalar Parameters in Subdiffusion with Memory via a Nonlocal Observation

基于非局部观测的带记忆性次扩散中标量参数的正则化重构

作者: Andrii Hulianytskyi, Sergei Pereverzyev, Sergii Siryk, Nataliya Vasylyeva

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

发布日期: 2025-11-07

arXiv ID: 2511.05277v1

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

摘要: In the paper, we propose an analytical and numerical approach to identify scalar parameters (coefficients, orders of fractional derivatives) in the multi-term fractional differential operator in time, $\mathbf{D}_t$. To this end, we analyze inverse problems with an additional nonlocal observation related to a linear subdiffusion equation $\mathbf{D}_{t}u-\mathcal{L}_{1}u-\mathcal{K}*\mathcal{L}_{2}u=g(x,t),$ where $\mathcal{L}_{i}$ are the second order elliptic operators with time-dependent coefficients, $\mathcal{K}$ is a summable memory kernel, and $g$ is an external force. Under certain assumptions on the given data in the model, we derive explicit formulas for unknown parameters. Moreover, we discuss the issues concerning to the uniqueness and the stability in these inverse problems. At last, by employing the Tikhonov regularization scheme with the quasi-optimality approach, we give a computational algorithm to recover the scalar parameters from a noisy discrete measurement and demonstrate the effectiveness (in practice) of the proposed technique via several numerical tests.

7. The Gross-Pitaewskii equation with time and space dependent coefficients

具有时空依赖系数的Gross-Pitaevskii方程

作者: Federico Lai

类别: math.AP

发布日期: 2025-09-07

arXiv ID: 2509.06001v3

链接: http://arxiv.org/abs/2509.06001v3

摘要: We study the existence of weak solutions of a generalized Gross-Pitaewskii equation, with time and space dependent coefficients that could blow up or vanish asymptotically in time, with initial data not necessarily segregated. We also study the asymptotic behavior in time of these solutions, including cases in which the segregation phenomenon appears.

8. Sharp conditions for the BBM formula and asymptotics of heat content-type energies

BBM公式与热含量型能量渐近性的精确条件

作者: Luca Gennaioli, Giorgio Stefani

类别: math.AP, math.FA, Primary 46E35. Secondary 26A33, 26D10, 35K08

发布日期: 2025-02-20

arXiv ID: 2502.14655v2

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

摘要: Given $p\in[1,\infty)$, we provide sufficient and necessary conditions on the non-negative measurable kernels $(\rho_t)_{t\in(0,1)}$ ensuring convergence of the associated Bourgain-Brezis-Mironescu (BBM) energies $(\mathscr{F}_{t,p})_{t\in(0,1)}$ to a variant of the $p$-Dirichlet energy on $\mathbb R^N$ as $t\to0^+$ both in the pointwise and in the $\Gamma$-sense. We also devise sufficient conditions on $(\rho_t)_{t\in(0,1)}$ yielding local compactness in $L^p(\mathbb R^N)$ of sequences with bounded BBM energy. Moreover, we give sufficient conditions on $(\rho_t)_{t\in(0,1)}$ implying pointwise and $\Gamma$-convergence and compactness of $(\mathscr{F}_{t,p})_{t\in(0,1)}$ when the limit $p$-energy is of non-local type. Finally, we apply our results to provide asymptotic formulas in the pointwise and $\Gamma$-sense for heat content-type energies both in the local and non-local settings.

9. Classification of fractional, singular Yamabe metrics on a twice punctured sphere I

两倍穿孔球面上分数阶奇异Yamabe度量的分类（I）

作者: João Henrique Andrade, Azahara DelaTorre, João Marcos do Ò, Jesse Ratzkin, Juncheng Wei

类别: math.DG, math.AP, 35J60, 35R11, 53C21, 53A30, 35B33, 35B40

发布日期: 2025-11-07

arXiv ID: 2511.05225v1

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

摘要: The Delaunay metrics form a family of conformally flat, constant fractional Q-curvature metrics on a twice-punctured sphere. They are all (after a M"obius transformation) rotationally symmetric and periodic, and admit several elegant variational descriptions. We prove that, when s is close to but less than 1, any complete, conformally flat constant Q-curvature metric on a twice-punctured sphere is a Delaunay metric. Along the way, we prove a sharp a priori bound for the conformal factor of these metrics, which may be of independent interest.

10. The Kato problem for weighted elliptic and parabolic operators of higher order

高阶加权椭圆与抛物型算子的Kato问题

作者: Guoming Zhang

类别: math.AP

发布日期: 2025-03-03

arXiv ID: 2503.01636v4

链接: http://arxiv.org/abs/2503.01636v4

摘要: We solve the Kato square root problem for parabolic operators of arbitrary order $2m$ whose coefficients are allowed to depend on both space and time in a merely measurable way and possess boundedness and ellipticity controlled by a Muckenhoupt $A_{2}-$weight. Notably, the proof applies to the weighted Kato problem within an elliptic framework.

11. Quenching phenomena in a system of non-local diffusion equations

非局部扩散方程组中的淬灭现象研究

作者: José M. Arrieta, Raúl Ferreira, Sergio Junquera

类别: math.AP

发布日期: 2025-03-31

arXiv ID: 2503.23994v2

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

摘要: In this paper we study the quenching phenomena occurring in a non-local diffusion system of two equations with intertwined singular absorption terms of the type $u^{-p}$. We prove that there exists a range of multiplicative parameters for which every solution presents quenching, while outside this range there are both global and quenching solutions. We also characterize in terms of the exponents of the absorption terms when the quenching is simultaneous or non-simultaneous and obtain the quenching rates.

12. Variational models of robust optimal transport

鲁棒最优传输的变分模型

作者: Luigi De Masi, Andrea Marchese, Annalisa Massaccesi

类别: math.AP, math.OC, 49Q10, 49Q20, 49Q22, 90B06

发布日期: 2025-11-07

arXiv ID: 2511.05146v1

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

摘要: This paper introduces two variational formulations for a model of robust optimal transport, that is, the problem of designing optimal transport networks that are resilient to potential damages, balancing construction costs against the benefit of maintaining partial functionality when parts of the network are damaged. We propose a Eulerian formulation, where the network is modeled by a rectifiable measure and recovery plans are represented by 1-dimensional normal currents. This framework allows for changes in the direction of the transportation in response to damages but restricts damages to be characteristic functions of closed sets. We also propose a Lagrangian formulation, where the network is a traffic plan (that is, a measure on the space of Lipschitz curves) and recovery plans are sub-traffic plans. This approach prescribes the network’s orientation but allows for a wider class of damages. We prove existence of minimizers in both settings. The two models are compared through examples that illustrate their main differences: the Eulerian formulation’s necessity for an unoriented network to achieve existence, the Lagrangian formulation’s ability to handle general damages and its requirement for a positive distance between the supports of the source and target measures.

13. Null controllability of damped nonlinear wave equation

阻尼非线性波动方程的零可控性

作者: Yan Cui, Peng Lu, Yi Zhou

类别: math.AP, math.OC

发布日期: 2022-11-28

arXiv ID: 2211.15449v4

链接: http://arxiv.org/abs/2211.15449v4

摘要: In this paper, we investigate the null controllability of nonlinear wave systems. Initially, we employ a combination of the Galerkin method and a fixed point theorem to establish the null controllability for semi-linear wave equations with nonlinear functions that are dependent on velocities, under the geometric control condition. Subsequently, utilizing a novel iterative method, we demonstrate the null controllability for a class of quasi-linear wave systems in a constructive manner. Lastly, we present a control result for a class of fully nonlinear wave systems, serving as an application.

14. Sharp remainder terms and stability of weighted Hardy-Poincaré and Heisenberg-Pauli-Weyl inequalities related to the Baouendi-Grushin operator

与Baouendi-Grushin算子相关的加权Hardy-Poincaré不等式和Heisenberg-Pauli-Weyl不等式的精确余项及其稳定性

作者: Yerkin Shaimerdenov, Nurgissa Yessirkegenov, Amir Zhangirbayev

类别: math.AP, 26D10, 35J70

发布日期: 2025-08-22

arXiv ID: 2508.16380v3

链接: http://arxiv.org/abs/2508.16380v3

摘要: In this paper, we obtain sharp remainder terms for the Hardy-Poincar'e inequalities with general non-radial weights in the setting of Baouendi-Grushin vector fields (see Theorem 2.5). It is worth emphasizing that all of our results are new both in the Baouendi-Grushin and standard Euclidean settings. The method employed allows us to not only unify, but also improve the results of Kombe and Yener [KY18] for any $1<p<\infty$ while holding true for complex-valued functions and providing explicit constants (Corollary 2.7). As a result, we are able to obtain sharp remainder terms to many known weighted Hardy-type inequalities (see Section 3.1). Aside from weighted Hardy-type inequalities, we also recover a sharp remainder formula for the $L^{p}$-Poincar'e inequality (Corollary 3.5). In the special case of radial weights, we are naturally able to introduce the notion of Baouendi-Grushin $p$-Bessel pairs (see Definition 2.9). Finally, we apply the technique to establish the sharp remainder term of the Heisenberg-Pauli-Weyl inequality in $L^{p}$ for $1<p<\infty$ (Corollary 3.13), which includes the sharp constant. This makes it possible to obtain the $L^{p}$-analogue for $2\leq p < n$ (Theorem 3.17) of a stability result by Cazacu, Flynn, Lam and Lu [CFLL24].

15. Inverse problem of determining a time-dependent coefficient in the time-fractional subdiffusion equation

确定时间分数阶次扩散方程中时变系数的反问题

作者: Ravshan Ashurov, Elbek Husanov

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.05011v1

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

摘要: This paper explores the forward and inverse problems for a fractional subdiffusion equation characterized by time-dependent diffusion and reaction coefficients. Initially, the forward problem is examined, and its unique solvability is established. Subsequently, the inverse problem of identifying an unknown time-dependent reaction coefficient is addressed, with rigorous proofs of the existence and uniqueness of its solution. Both problems’ existence and uniqueness are demonstrated using Banach’s contraction mapping theorem. Notably, this work is the first to investigate direct and inverse problems for such equations with time-dependent coefficients.

16. Uniqueness for phaseless inverse elastic scattering problem for periodic structures

周期结构弹性散射相位缺失反问题的唯一性

作者: Youzi He, Wei Wu, Hongyi Dang

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.04999v1

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

摘要: This paper establishes uniqueness results of inverse elastic scattering problem with phaseless near-field data in periodic structures in $\mathbb{R}^2$ and periodic/biperiodic structures in $\mathbb{R}^3$. We use a superposition of two point sources in each periodic unit with different positions as the incident field, and measures the phaseless near-field data on a line parallel to $x_1$-axis in 2D, or on a plane parallel to $(x_1,x_2)$-plane in 3D. We first calculate the explicit formula of quasi-periodic/biperiodic Green’s functions of Lam'{e} system in $\mathbb{R}^3$. Then, to establish the uniqueness results, the reciprocity relations for point sources, scattered fields, and total fields are derived. Finally, with the help of Rayleigh’s expansion, the uniqueness results are proved. The quasi-periodic/biperiodic Green’s functions of Lam'{e} system in $\mathbb{R}^3$, the reciprocity relations, and Rayleigh’s expansion in $\mathbb{R}^3$ are novel results as important by-products in the proof process.

17. Jump problem for generalized Lamé-Navier systems in $\mathbb{R}^m$

$\mathbb{R}^m$中广义Lamé-Navier系统的跳跃问题

作者: Daniel Alfonso Santiesteban, Ricardo Abreu Blaya, Daniel Alpay

类别: math.AP, 30G35

发布日期: 2025-11-07

arXiv ID: 2511.04959v1

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

摘要: This paper is devoted to study a fundamental system of equations in Linear Elasticity Theory: the famous Lam'e-Navier system. The Clifford algebra language allows us to rewrite this system in terms of the Euclidean Dirac operator, which at the same time suggests a very natural generalization involving the so-called structural sets. Our interest lies mainly in the jump problem for these elastic systems. A generalized Teodorescu transform, to be introduced here, provides the means for obtaining the explicit solution of the jump problem for a very wide classes of regions, including those with a fractal boundary.

18. Positive solutions of elliptic systems with superlinear nonlinearities on the boundary

具有超线性非线性的椭圆系统在边界上的正解

作者: Shalmali Bandyopadhyay, Maya Chhetri, Briceyda Delgado, Nsoki Mavinga, Rosa Pardo

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.04943v1

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

摘要: We consider elliptic systems with superlinear and subcritical boundary conditions and a bifurcation parameter as a multiplicative factor. By combining the rescaling method with degree theory and elliptic regularity theory, we prove the existence of a connected branch of positive weak solutions that bifurcates from infinity as the parameter approaches zero. Furthermore, under additional conditions on the nonlinearities near zero, we obtain a global connected branch of positive solutions bifurcating from zero, which possesses a unique bifurcation point from infinity when the parameter is zero. Finally, we analyze the behavior of this branch and discuss the number of positive weak solutions with respect to the parameter using bifurcation theory, degree theory, and sub- and super-solution methods.

19. Representation formula, regularity, and decay of solutions for sub-diffusion equations

次扩散方程解的表示公式、正则性与衰减性

作者: Sandro Coriasco, Giovanni Girardi, Stevan Pilipović

类别: math.AP

发布日期: 2025-11-07

arXiv ID: 2511.04885v1

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

摘要: We study regularity and decay properties for the solutions of the Cauchy problem for time-fractional partial differential equations, with tempered initial data, belonging to suitable (weighted) Sobolev spaces, associated with a differential operator on space variables with polynomially bounded coefficients. We obtain a representation formula for the solution, modulo time-regular functions, smooth and rapidly decreasing with respect to the space variables. By means of the representation formula, the (decay and smoothness) singularities of the solution of the homogeneous Cauchy problem can be controlled, in terms of (global) wavefront sets of the initial data.

20. Integrated Local Energy Decay for Waves with Time-Dependent Damping

带有时变阻尼的波方程的局部能量综合衰减

作者: Perry Kleinhenz, Michael McNulty

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

发布日期: 2025-09-30

arXiv ID: 2510.00179v2

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

摘要: We prove integrated local energy decay for solutions of the damped wave equation with time-dependent damping satisfying an appropriate generalization of the geometric control condition on asymptotically flat, stationary space-times. We first obtain a high frequency estimate, which we prove via a positive commutator estimate using an escape function explicitly constructed in terms of the damping around individual space-time trajectories. We combine the high frequency estimate with low and medium frequency results for the undamped problem, then we handle the damping term as a perturbation to obtain local energy decay.

21. A nonlinear homogenization-based perspective on the soft modes and effective energies of some conformal metamaterials

基于非线性均匀化的共形超材料软模与有效能量研究视角

作者: Xuenan Li, Robert V. Kohn

类别: math.AP, math-ph, math.MP, 49Nxx

发布日期: 2025-09-21

arXiv ID: 2509.16907v2

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

摘要: There is a growing mechanics literature concerning the macroscopic properties of mechanism-based mechanical metamaterials. This amounts mathematically to a homogenization problem involving nonlinear elasticity. A key goal is to identify the “soft modes” of the metamaterial. We achieve this goal using methods from homogenization for some specific 2D examples – including discrete models of the Rotating Squares metamaterial and the Kagome metamaterial – whose soft modes are compressive conformal maps. The innovation behind this achievement is a new technique for bounding the effective energy from below, which takes advantage of the metamaterial’s structure and symmetry.

22. Schrodinger-Poisson-Slater equations with nonlinearity subscaled near zero

非线性在零附近次标度化的薛定谔-泊松-斯莱特方程

作者: Shibo Liu, Kanishka Perera

类别: math.AP, 35J91, 35J20, 58E05

发布日期: 2025-07-24

arXiv ID: 2507.18568v3

链接: http://arxiv.org/abs/2507.18568v3

摘要: We study the following zero-mass Schr{"o}dinger-Poisson-Slater equation

with nonlinearity subscaled near zero in the sense that $f (| x |, t) \approx a | t |^{p - 2} t$ as $| t | \rightarrow 0$ for some $p\in\big(\frac{18}{7},3\big)$. A nonzero solution is obtained via Morse theory when the nonlinearity is asymptotically scaled at infinity. For this purpose we prove an abstract result on the critical groups at infinity for functionals satisfying the geometric assumptions of the scaled saddle point theorem of Mercuri & Perera [arXiv:2411.15887]. For the case that $f (| x |, \cdot)$ is odd, a sequence of solutions are obtained via a version of Clark’s theorem due to Kajikiya [J.
Funct.\ Anal.\ 225 (2005) 352–370].

23. Dynamics of the Energy-Critical Nonlinear Schrödinger System in ${\mathbb R}^{4}$

四维欧几里得空间中能量临界非线性薛定谔系统的动力学

作者: Alex H. Ardila

类别: math.AP

发布日期: 2025-11-06

arXiv ID: 2511.04839v1

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

摘要: In this paper, we investigate the dynamics of radial solutions at threshold energy for a 3-component Schr"{o}dinger system with cubic nonlinearity in four dimensions. The main difference from the cases previously addressed in the literature is that, in our system, the kernel of the imaginary part $L_I$ of the linearized operator $-i{\mathcal L}=L_{R}+iL_{I}$ has dimension 2. To overcome this difficulty, we carry out a detailed study of the coercivity properties of these operators. We also introduce a new modulation parameter associated with the additional eigenfunction in the kernel of the operator $L_{I}$, which enables us to perform the modulation analysis and establish the uniqueness of exponentially decaying solutions to the linearized equation.

24. Criticality, splitting theorems under spectral Ricci bounds and the topology of stable minimal hypersurfaces

谱Ricci曲率条件下的临界性、分裂定理及稳定极小超曲面的拓扑

作者: Giovanni Catino, Luciano Mari, Paolo Mastrolia, Alberto Roncoroni

类别: math.DG, math.AP

发布日期: 2024-12-17

arXiv ID: 2412.12631v4

链接: http://arxiv.org/abs/2412.12631v4

摘要: In this paper we prove general criticality criteria for operators $\Delta + V$ on manifolds with more than one end, where $V$ bounds the Ricci curvature, and a related spectral splitting theorem extending Cheeger-Gromoll’s one. Our results give new insight on Li-Wang’s theory of manifolds with a weighted Poincar'e inequality. We apply them to study stable and $\delta$-stable minimal hypersurfaces in manifolds with non-negative bi-Ricci or sectional curvature, in ambient dimension up to $5$ and $6$, respectively. In the special case where the ambient space is $\mathbb{R}^4$, we prove that a $1/3$-stable minimal hypersurface must either have one end or be a catenoid, and that proper, $\delta$-stable minimal hypersurfaces with $\delta > 1/3$ must be hyperplanes.

25. Fractional Schrödinger-Poisson-Slater equations in Coulomb-Sobolev spaces

库仑-索博列夫空间中的分数阶薛定谔-泊松-斯莱特方程

作者: Elisandra Gloss, Carlo Mercuri, Kanishka Perera, Bruno Ribeiro

类别: math.AP, math.FA, 35R11, 35J60, 35A15, 35B33, 35J20

发布日期: 2025-11-06

arXiv ID: 2511.04829v1

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

摘要: We prove existence and multiplicity results for the fractional Schroedinger–Poisson–Slater equation $(-\Delta)^s u + (I_\alpha * u^2)u = f(|x|,u)$ in $\mathbb{R}^N$, where $0<s<1$ and $\alpha \in (1,N)$. We seek solutions in a fractional Coulomb-Sobolev space and employ new tools in critical point theory that link the behavior of $f$ at zero and at infinity to the scaling properties of the left-hand side. For several regimes of $f$, we establish compactness for an associated action functional and obtain multiple solutions as critical points, with the number governed by the interaction of $f$ with a sequence of eigenvalues $\{\lambda_k\}$ defined via the $\mathbb{Z}_2$ cohomological index of Fadell and Rabinowitz (rather than the classical Krasnosel’skii genus). In this fractional setting we also prove new regularity results and necessary conditions for the existence of solutions.

26. The strongly nonlocal Allen-Cahn problem

强非局部Allen-Cahn问题

作者: Erisa Hasani, Stefania Patrizi

类别: math.AP

发布日期: 2025-11-06

arXiv ID: 2511.04818v1

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

摘要: We study the sharp interface limit of the fractional Allen-Cahn equation

where $\varepsilon >0$, $\mathcal{I}^s_n=-c_{n,s}(-\Delta )^s$ is the fractional Laplacian of order $2s\in(0,1)$ in $\mathbb{R}^n$, and $W$ is a smooth double-well potential with minima at 0 and 1. We focus on the singular regime $s\in(0,\frac{1}{2})$, corresponding to strongly nonlocal diffusion. For suitably prepared initial data, we prove that the solution $ u^\varepsilon $ converges, as $\varepsilon\to0$, to the minima of $W$ with the interface evolving by fractional mean curvature flow. This establishes the first rigorous convergence result in this regime, complementing and completing previous work for $s\geq \frac{1}{2}$.

27. Self-consistent equations and quantum diffusion for the Anderson model

安德森模型中的自洽方程与量子扩散

作者: Adam Black, Reuben Drogin, Felipe Hernández

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

发布日期: 2025-06-06

arXiv ID: 2506.06468v2

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

摘要: We consider the Anderson tight-binding model on $\mathbb{Z}^d$, $d\geq 2$, with Gaussian noise and at low disorder $\lambda>0$. We derive a diffusive scaling limit for the entries of the resolvent $R(z)$ at imaginary part $\operatorname*{Im} z\sim\lambda^{2+\kappa_d}$, $\kappa_d>0$, with high probability. As consequences, we establish quantum diffusion (in a time-averaged sense) for the Schr"{o}dinger propagator at the longest timescale known to date and improve the best available lower bounds on the localization length of eigenfunctions. Our results for $d=2$ are the first quantum diffusion results for the Anderson model on $\mathbb{Z}^2$. The proof avoids the use of diagrammatic expansions and instead proceeds by analyzing certain self-consistent equations for $R(z)$. This is facilitated by new estimates for $\|R(z)\|_{\ell^p\rightarrow \ell^q}$ that control the recollisions.

28. A priori estimates and $η-$compactness for anisotropic Ginzburg-Landau minimizers with tangential anchoring

各向异性Ginzburg-Landau极小化子在切向锚定条件下的先验估计与$η-$紧性

作者: Lia Bronsard, Andrew Colinet, Dominik Stantejsky, Lee van Brussel

类别: math.AP, math.CA, 49K20, 35B45, 35B38, 35E20, 49S05

发布日期: 2025-11-06

arXiv ID: 2511.04672v1

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

摘要: We consider minimizers $u_\varepsilon$ of the Ginzburg-Landau energy with quadratic divergence or curl penalization on a simply-connected two-dimensional domain $\Omega$. On the boundary, strong tangential anchoring is imposed. We prove a priori estimates for $u_\varepsilon$ in $L^\infty$ uniform in $\varepsilon$ and that the Lipschitz constant of $u_\varepsilon$ blows up like $\varepsilon^{-1}$. We then deduce compactness for a subsequence that converges to an $\mathbb{S}^1-$valued map with either one interior point defect or two boundary half-defects. We conclude our study with a proof that no boundary vortices can occur in the divergence penalized case.

29. Generic regularity of intermediate complex structure limits

中间复结构极限的通有正则性

作者: Yang Li, Valentino Tosatti

类别: math.DG, math.AP, 53C55, 14J32, 32Q25, 35J96

发布日期: 2025-11-06

arXiv ID: 2511.04651v1

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

摘要: We study certain polarized degenerations of Calabi-Yau manifolds near an intermediate complex structure limit, and improve the potential $C^0$-convergence to a metric convergence result on the generic region for the corresponding collapsing Ricci-flat K"ahler metrics.