生成时间: 2025年10月31日 03:44:06

重点关注论文数量: 3

了解领域论文数量: 4

不相关论文数量: 21

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

1. On action rate admissibility criteria

论作用率可接受性准则

作者: Heiko Gimperlein, Michael Grinfeld, Robin J. Knops, Marshall Slemrod

类别: math.AP

发布日期: 2025-03-05

arXiv ID: 2503.03491v4

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

详细分析

中文标题: 论作用率可接受性准则

1. 研究对象和背景

本文研究的是二维等熵可压缩欧拉流体方程的黎曼初值问题，其核心目标是提出一种新的解的选择准则（admissibility criterion），以从众多非唯一弱解中识别出物理上合理的经典激波解——即所谓的“2-激波解”（2-shock solution）。这一问题源于近年来在流体力学偏微分方程理论中出现的一个深刻现象：利用凸积分方法（convex integration）可以构造出大量满足能量不等式但高度震荡、非光滑的“野解”（wild solutions），这些解与经典激波解共享相同的初始数据和守恒律，却表现出完全不同的动力学行为。

该系统描述的是无粘、可压缩气体的运动。考虑二维空间中的质量密度 $\rho(t,x) > 0$ 和速度场 $v(t,x) \in \mathbb{R}^2$，其中 $(t,x) \in [0,\infty) \times \mathbb{R}^2$。系统的控制方程为等熵可压缩欧拉方程组：

$$ \partial_t \rho + \mathrm{div}_x (\rho v) = 0, $$
$$ \partial_t (\rho v) + \mathrm{div}_x (\rho v \otimes v) + \nabla_x p(\rho) = 0. $$

这里，压力函数取为幂律形式 $p(\rho) = \rho^\gamma$，其中绝热指数 $\gamma$ 满足 $1 < \gamma \leq 3$。对应的比内能 $\varepsilon(\rho)$ 定义为 $\varepsilon'(\rho) = p(\rho)/\rho^2$，从而总能量密度为 $\rho \varepsilon(\rho) + \frac{1}{2}\rho |v|^2$。

初始条件为经典的黎曼型初值，即在 $t=0$ 时，状态在 $x_1=0$ 处发生跳跃：

$$ (\rho, v)(0, x) = \begin{cases} (\rho_-, v_-), & x_1 < 0, \\ (\rho_+, v_+), & x_1 > 0, \end{cases} $$

且假设 $v_\pm = (0, v_\pm^{(2)})$，即速度仅在 $x_2$ 方向有分量，这保证了问题具有某种对称性。这种设定下，存在一个自相似的经典解——2-激波解，它由两个向外传播的激波构成，在中间区域形成一个恒定的中间态 $(\rho_m, v_m)$。

然而，Chiodaroli 与 Kreml 的开创性工作表明，当满足一定条件（如文中所提的不等式 (3.6)）时，除了这个经典解外，还存在无穷多个满足能量耗散不等式的弱解（weak solutions），这些解通过凸积分技术构造而来，被称为“野解”。更令人震惊的是，某些野解的能量耗散率甚至低于经典2-激波解，因此无法被传统的Dafermos熵速率准则（entropy rate admissibility criterion）排除。

此外，Markfelder 与 Pellhammer 进一步构造了所谓的“混合解”（hybrid solutions），它们先经历一段类似野解的动力学过程（包含高频振荡），然后切换回经典波结构。由于这类解在长时间内的总作用量可能小于2-激波解，导致先前提出的最小作用原理选择准则（Least Action Admissibility Principle, LAAP）也失效。

因此，当前的研究背景是一个典型的非唯一性危机：对于同一个物理初始状态，数学上允许多个看似合理但行为迥异的解共存。这就迫切需要一个新的、更具区分力的可接受性准则来筛选出物理上真实发生的解。

论文特别关注的函数类是所有满足相同初始数据的熵弱解（entropic weak solutions），即同时满足动量守恒和能量不等式 (3.4) 的分布意义下的解。研究区域为全平面 $\mathbb{R}^2$，时间区间为 $[0,T]$ 或局部邻域 $[t_0,t]$，未施加额外边界条件，因问题本身具有平移不变性和自相似结构。

关键难点在于：传统基于全局能量或作用量极小化的准则（如LAAP）依赖于终时间 $t_1$ 的选择，而混合解正是利用这一点，在某个时间段后“回归正统”，使得其总作用量最终低于经典解。因此，作者提出应转向局部时间尺度上的比较，考察解在初始时刻附近的作用量增长速率。

2. 主要定理或主要结果

本文的核心成果并非以标准“定理”形式陈述，而是通过一系列定义和命题构建了一个新的可接受性框架，并应用于具体问题得出明确结论。其主要结果可归纳如下：

主要结果一（新准则的定义与关系）

引入三个递进的可接受性准则：

  • LAAP₀（局部最小作用原则）：设 $S$ 是某一演化方程初值问题的解集。称解 $u \in S$ 是 LAAP₀-可接受的，如果对任意其他解 $v \in S$, $v \ne u$，存在时间 $t_1(v) > t_0$，使得在区间 $(t_0, t_1)$ 内总有：

    $$ A(u)(t_0, t) \leq A(v)(t_0, t) $$

    若严格不等式在某子区间成立，则称 $u$ 严格偏好于 $v$。

  • ARAC（作用率可接受性准则）：若作用量 $t \mapsto A(u)(t_0, t)$ 在 $t = t_0^+$ 可导，则称 $u$ 满足 ARAC，若对所有 $v \ne u$ 有：

    $$ \left.\frac{d}{dt}A(u)(t_0, t)\right|_{t=t_0^+} \leq \left.\frac{d}{dt}A(v)(t_0, t)\right|_{t=t_0^+} $$

  • sARAC（严格作用率可接受性准则）：存在最小整数 $k \in \mathbb{N}$，使得前 $k-1$ 阶导数相等，但第 $k$ 阶导数满足：

    $$ \left.\frac{d^k}{dt^k}A(u)(t_0, t)\right|_{t=t_0^+} < \left.\frac{d^k}{dt^k}A(v)(t_0, t)\right|_{t=t_0^+} $$

并证明了以下逻辑蕴含关系：

Proposition 2.1:

  1. 若 $u$ 是 LAAP₀-可接受的，则 $u$ 满足 ARAC；
  2. 若 $u$ 满足 sARAC，则 $u$ 是严格 LAAP₀-可接受的。

主要结果二（应用结论）

在二维等熵可压缩欧拉方程的特定黎曼问题中，考虑包含以下三类解的集合 $S$：

  • 经典的2-激波解；
  • Chiodaroli-Kreml 构造的凸积分野解；
  • Markfelder-Pellhammer 构造的混合解；

则对于适当选取的参数范围（如 $\rho_1$ 接近 $\rho_m$ 且 $\varepsilon_2 > 0$ 小量），2-激波解是唯一满足 sARAC（从而也是 LAAP₀-可接受）的解。

换句话说，尽管这些野解和混合解在全局作用量上可能更小，但在初始时刻 $t=0^+$ 的作用量增长率上，2-激波解的增长最慢，因而被新准则选中。

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

a) 核心分析学工具

本文的方法论融合了变分法思想、分布理论与渐近分析，其核心工具围绕“作用量”的定义及其时间导数展开。

首先，作用量（action）被定义为拉格朗日量在时空区域上的积分。对于本文考虑的流体系统，拉格朗日量为：

$$ \mathcal{L} = \frac{1}{2}\rho |v|^2 - \rho \varepsilon(\rho) $$

即动能减去势能（内部能量）。在区域 $\Lambda = \Omega \times [t_0, t]$ 上的作用量为：

$$ A(\rho, v)(t_0, t) = \int_{t_0}^{t} \int_{\Omega} \mathcal{L} \, dx d\tau $$

这是经典力学中最小作用原理的直接推广。虽然欧拉方程本身不是某个作用量泛函的临界点（因其非保守），但作者仍借用此框架作为选择准则的基础。

其次，为了处理弱解之间的比较，必须使用分布意义下的解概念。文中所有解（包括野解和混合解）都被理解为满足守恒律和能量不等式 (3.4) 的分布解。特别是凸积分构造的解，本质上是通过极限过程得到的，依赖于弱收敛和振荡补偿效应。

第三，新准则的关键在于右导数的存在性与比较。为此，作者假设作用量函数 $t \mapsto A(u)(t_0, t)$ 在 $t = t_0^+$ 足够光滑（至少一阶可导），从而可以比较不同解的作用量增长率。这属于局部时间渐近分析的一种形式。

最后，文中隐含使用了参数化技巧与扰动分析。通过对扇形子解（fan sub-solution）进行参数化（$\rho_1, \varepsilon_2$），并将野解视为对经典解的小扰动，作者能够计算两者之间拉格朗日量的差异：

$$ L_{\mathrm{diff}}(\rho_1) := \mathcal{L}_{\mathrm{2-shock}} - \mathcal{L}_{\mathrm{wild}} $$

并在局部区域内积分以估计作用量差的增长率。

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

本文面临的主要技术挑战是如何在高度非光滑、非唯一解的集合中建立有效的排序机制。传统全局准则（如LAAP）失败的根本原因在于非紧性与长时行为主导短时判断。

为此，作者采用的关键技巧是：

  1. 从全局到局部的视角转换：放弃对终时间 $t_1$ 固定的作用量极小化，转而考察 $t \to t_0^+$ 时的瞬时增长率。这种“初瞬间作用率最小化”的思想类似于 Dafermos 的熵速率准则，但将“熵”替换为“作用量”。

  2. 排除混合解的构造机制：Markfelder-Pellhammer 的混合解之所以能击败LAAP，是因为它们在后期切换为低能量结构，从而降低总作用量。但由于这种切换发生在有限时间之后，其作用量在 $t=0^+$ 附近的导数仍由前期的高频振荡部分决定，而这一部分通常具有较高的动能（因速度震荡），导致初始作用率较大。因此，LAAP₀ 自动排除了这类解。

  3. 利用子解参数化进行显式估计：通过将野解构造中的常数 $C$ 表达为 $\rho_1$ 和 $\varepsilon_2$ 的函数（见公式 (3.19)），并结合已知关系 $\beta^2(\rho_1) + \varepsilon_1(\rho_1) + \varepsilon_2 = C$，作者可以在 $\rho_1$ 接近 $\rho_m$ 时进行泰勒展开，精确估算作用量差异的一阶项。

  4. 对称性与区域分解：由于问题在 $x_2$ 方向上均匀，作者可在有限宽度 $[-L_3, L_3] \times [\ell_1, \ell_2]$ 上计算作用量，避免发散。同时，将空间划分为三个扇区 $P_-, P_1, P_+$，分别处理经典解与野解的行为差异。

  5. 利用凸积分的周期性结构：文中指出（Remark 3.1(2)），corrugation 解可取为 $x_1$-方向周期函数。这意味着在平均意义上，其附加速度场的平方贡献一个正的常数项，增加动能，进而提高作用量增长率。

综上，作者通过将变分思想本地化、结合参数扰动分析、利用构造细节进行显式估计，成功绕开了凸积分带来的非唯一性障碍。

4. 与之前工作的比较

a) 相对于传统熵准则的改进

传统的熵条件（如Lax条件、Liu条件）主要用于一维标量或双曲守恒律，而在高维或多维系统中难以直接应用。Dafermos 提出的熵速率准则（entropy rate admissibility criterion）主张选择使熵产生率最大的解，其物理动机是自然界倾向于最快达到平衡的状态。

然而，Chiodaroli-Kreml 的研究表明，某些野解的熵耗散率反而低于2-激波解，因此不被该准则排除。这暴露了熵速率准则在高维可压缩流中的局限性：它不能有效抑制高频震荡解的出现。

本文提出的 LAAP₀ / ARAC 准则 则提供了一种替代路径：不是看谁“更快耗散能量”，而是看谁“最少增加作用量”。由于作用量包含动能项 $\frac{1}{2}\rho |v|^2$，而野解往往伴随速度的剧烈震荡，即使平均动量为零，其速度平方的平均值也会显著增大，导致作用量快速增长。因此，该准则天然地惩罚震荡行为，从而偏好光滑或分片常数的激波结构。

b) 相对于原始LAAP的突破

在作者之前的论文 [9] 中提出的 LAAP 原则要求在整个时间区间 $[t_0, t_1]$ 上作用量最小。虽然在某些情况下能选出2-激波解，但它严重依赖于终时间 $t_1$ 的选择。正如 Markfelder-Pellhammer 所示，只要让混合解在足够大的 $t_1$ 后趋于稳定低能态，就能使其总作用量低于经典解。

本文提出的 LAAP₀ 是对LAAP的根本性改进：

  • 从全局变为局部：不再依赖终时间，只关心初始时刻的动态响应；
  • 更强的稳定性：不受后期行为影响，避免被“延迟优势”欺骗；
  • 更具物理直觉：类似于“惯性最小”或“启动成本最低”的原则，符合经典解作为“最自然演化路径”的直觉。

更重要的是，LAAP₀ 成功将 Markfelder-Pellhammer 的混合解排除在外，因为它们的“启动阶段”本质上仍是野解，具有高作用率。

c) 对凸积分构造的回应

凸积分技术的强大之处在于它揭示了弱解集合的病态丰富性：即使满足能量不等式，解也不唯一。这动摇了经典PDE理论中“熵解唯一”的信念。

本文并未否定这些解的数学存在性，而是提出：存在性不等于物理可实现性。通过引入基于作用率的新准则，作者实际上在说：“虽然你们都是弱解，但只有那个启动最平稳的才是自然界会选择的。”

这种方法的优势在于：

  • 不需要提高正则性假设（如要求解属于 $BV$ 或 Sobolev 空间）；
  • 不改变原方程或引入人工粘性；
  • 保留了所有熵弱解作为候选，仅通过一个额外的变分性判据进行筛选。

总之，本文的突破体现在范式转变：从“寻找最强正则性条件以保证唯一性”转向“在非唯一解集中建立优先序”。这是一种更为实用和物理导向的研究策略，尤其适用于那些已被证明存在多重弱解的重要流体力学模型。

2. Long-time Confinement near Special Vortex Crystals

特殊涡旋晶体附近的长期约束

作者: Martin Donati

类别: math.AP, 35Q31, 76B47, 34D20

发布日期: 2024-10-18

arXiv ID: 2410.14221v3

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

详细分析

中文标题: 特殊涡旋晶体附近的长期约束

1. 研究对象和背景

本文研究的是二维不可压缩、无粘性流体中涡量集中于特殊涡旋晶体结构附近的长时间动力学行为。所涉及的基本方程是二维欧拉方程（incompressible Euler equations）在涡量形式下的表达式，其数学模型如下：

$$ \begin{cases} \partial_t \omega(x,t) + u(x,t) \cdot \nabla \omega(x,t) = 0 & \text{on } \mathbb{R}^2 \times \mathbb{R}_+, \\ u(x,t) = (K * \omega(t))(x), & \end{cases} \tag{Eu} $$

其中涡量场 $\omega : \mathbb{R}^2 \times \mathbb{R}_+ \to \mathbb{R}$ 满足初始条件 $\omega(0) = \omega_0$，速度场 $u$ 通过Biot-Savart律由涡量卷积得到，核函数为：

$$ K(x) = \frac{x^\perp}{2\pi |x|^2}, \quad x^\perp := (-x_2, x_1). $$

该系统描述了理想流体中涡量沿速度场的输运过程：$\omega$ 在拉格朗日轨迹上保持守恒。

论文关注一类特殊的解，其初始涡量集中在若干小区域（称为“涡旋斑块”或“vortex blobs”）内，且这些斑块的位置接近某个特定的点涡构型（point-vortex configuration）——即所谓的涡旋晶体（vortex crystals）。所谓“晶体”，是指这些点涡以刚体方式旋转运动，彼此间距离不变，形成稳定的几何构型。

具体地，考虑一族初始数据 $(\omega_\varepsilon)_\varepsilon$，满足以下假设（文中记为 Hypothesis 2.1）：

  • $\omega_\varepsilon \in L^1 \cap L^\infty(\mathbb{R}^2)$；
  • 总涡量分解为 $N$ 个非负部分：$\omega_\varepsilon = \sum_{i=1}^N \omega_{\varepsilon,i}$，每个分量满足：
    $$ \int_{\mathbb{R}^2} \omega_{\varepsilon,i}(x)\,dx = \gamma_i, $$
    即第 $i$ 个斑块具有固定强度 $\gamma_i \in \mathbb{R}\setminus\{0\}$；
  • 支集包含在半径为 $\varepsilon$ 的圆盘内：$\mathrm{supp}\,\omega_{\varepsilon,i} \subset D(z_i, \varepsilon)$，其中 $D(z,r) := \{x \in \mathbb{R}^2 : |x - z| < r\}$；
  • 涡量大小受控：$|\omega_\varepsilon| \leq M \varepsilon^{-\eta}$，$\eta \geq 2$。

这表示当 $\varepsilon \to 0$ 时，每个斑块趋于一个狄拉克测度：$\omega_\varepsilon \rightharpoonup \sum_{i=1}^N \gamma_i \delta_{z_i}$，弱收敛于测度意义下。

此时，在经典结果（如[12]）支持下，对应的解 $\omega_\varepsilon(t)$ 的演化在短时间尺度上逼近**点涡系统（Point-Vortex System, PVS）**的动力学：

$$ \frac{d}{dt} z_i(t) = \sum_{j \ne i} \gamma_j K(z_i(t) - z_j(t)) = \sum_{j \ne i} \gamma_j \frac{(z_i(t) - z_j(t))^\perp}{2\pi |z_i(t) - z_j(t)|^2}. \tag{PVS} $$

这是一个哈密顿系统，描述 $N$ 个理想点涡之间的相互作用。然而，该系统的解并不总是全局存在的，可能因多个涡旋发生**碰撞（collapse）**而爆破。但若总环量满足适当非共振条件（如所有子集环量和非零），则几乎处处初值不会导致碰撞。

本文特别关注一类能产生相对平衡（relative equilibrium）的构型，即涡旋晶体：系统整体作刚体旋转，形状不变。设初始位置为：

$$ z_k^* = \zeta^k, \quad k = 1,\dots,N-1; \quad z_N^* = 0, \quad \text{其中 } \zeta = e^{2\pi i / (N-1)}, $$

即前 $N-1$ 个涡位于单位圆上的正 $(N-1)$ 边形顶点，最后一个位于中心。外部涡强度均取为 1，中心涡强度 $\gamma_N$ 待定。

这类构型被称为“带中心涡的正多边形涡旋阵列”。文献 [4] 已证明，当 $\gamma_N$ 落在某个区间内时，该构型具有线性和非线性稳定性。但本文不取整个稳定区间，而是选择一个特定值：

$$ \gamma_N = \frac{(N-2)(N-6)}{12}, \tag{3.3} $$

使得每个涡所受来自其他涡的应变率（strain）在一级近似下为零。这种构型表现出更强的稳定性，从而允许对连续介质解进行更精细的长时间控制。

值得注意的是，当 $N=6$ 时，$\gamma_6 = 0$，意味着中心无涡，对应一个正五边形涡旋结构，这也是唯一一个无需中心涡即可实现强稳定性的例子。

目标问题是：若初始涡量集中于这样的稳定涡旋晶体附近（支集在 $\varepsilon$ 尺度），那么其演化过程中，各斑块是否能在很长一段时间内仍保持在 $\varepsilon^\beta$（$\beta < 1/2$）的小邻域内？换言之，是否存在长时间的空间约束（confinement）现象？

为此定义约束时间：

$$ \tau_{\varepsilon,\beta} := \sup \left\{ t \geq 0 : \forall s \in [0,t],\ \mathrm{supp}\,\omega_\varepsilon(s) \subset \bigcup_{i=1}^N D(z_i(s), \varepsilon^\beta) \right\}, $$

即解的支集始终被限制在随点涡轨迹移动的 $\varepsilon^\beta$ 圆盘内的最长时间。

以往结果（如[3]）表明一般情况下 $\tau_{\varepsilon,\beta} \gtrsim |\ln \varepsilon|$，这是最优阶；但在某些特殊情形（如单涡、自相似扩张构型）可提升至 $\varepsilon^{-\alpha}$ 阶。本文的核心贡献正是将此类代数阶长时间约束推广到一类新的稳定涡旋晶体构型，而不要求系统扩张或孤立。

函数空间方面，研究的是 $L^1 \cap L^\infty(\mathbb{R}^2)$ 中的整体弱解，这类解的存在唯一性由Yudovich理论保证（见[18]）。由于涡量仅被输运，故 $L^p$ 范数守恒，特别是支集始终保持有界。

边界条件方面，问题设定在全平面 $\mathbb{R}^2$ 上，无物理边界，属于无界区域问题。这也带来了长程相互作用的技术挑战，尤其是远场影响不能忽略。

综上所述，本文的研究对象是一个在全平面上演化的二维欧拉方程，初始涡量高度集中于接近某种特殊稳定涡旋晶体的 $N$ 个小区间内，目标是证明这些斑块在极长时间尺度（$\varepsilon^{-\alpha}$）内不会显著扩散或偏离其理想轨迹。

2. 主要定理或主要结果

本文的主要结果是以下两个定理的组合，分别从抽象框架与具体构造角度给出。

Theorem 2.6（主定理）

对于任意 $N \geq 4$，存在唯一的中心涡强度 $\gamma_N \in \mathbb{R}$（由公式 (3.3) 给出），使得：对任何满足 Hypothesis 2.1 的初始族 $(\omega_\varepsilon)_\varepsilon$，其对应的欧拉方程解 $\omega_\varepsilon(t)$ 满足如下长期约束估计：

对任意 $\beta < 1/2$，任意 $\alpha < \min(\beta/2, 2 - 4\beta)$，存在 $\varepsilon_0 > 0$，使得对所有 $\varepsilon \in (0, \varepsilon_0)$，

$$ \tau_{\varepsilon,\beta} \geq \varepsilon^{-\alpha}. $$

这意味着：只要初始斑块足够小（$\varepsilon \ll 1$），它们将在至少 $\varepsilon^{-\alpha}$ 的时间内保持在宽度为 $\varepsilon^\beta$ 的轨道管内。

特别地，当 $N=6$ 时，$\gamma_6 = 0$，对应一个正五边形涡旋结构（无中心涡），也满足此结论。

Theorem 3.4（一般性约束定理）

设点涡构型 $Z^* = (z_1^*, \dots, z_N^*)$ 满足以下两个关键假设：

  • Hypothesis 3.1（零一阶应变）：

    $$ \sum_{j \ne i} \gamma_j J_K(z_i^* - z_j^*) = 0, \quad \forall i = 1,\dots,N, $$

    其中 $J_K(x)$ 是 Biot-Savart 核 $K(x)$ 的雅可比矩阵。这一条件确保每个涡在其位置处感受到的流场梯度（即局部应变）为零，从而避免斑块因剪切而快速变形。

  • Hypothesis 3.3（线性稳定性）： 构型 $Z^*$ 是旋转坐标系下的平衡态，且其稳定性矩阵 $S = \Gamma^{-1} J D^2 H(Z^*) + \nu J$ 在不变子空间 $V^\perp$ 上可对角化为纯虚数块：

    $$ S|_{V^\perp} \sim \begin{pmatrix} 0 & b \\ -b & 0 \end{pmatrix}, $$

    表示扰动不会指数增长。

则对任何满足 Hypothesis 2.1 的初始数据族，其解满足 $\tau_{\varepsilon,\beta} \geq \varepsilon^{-\alpha}$，其中 $\alpha < \min(\beta/2, 2 - 4\beta)$。

Theorem 3.5（构造性验证）

对上述正多边形构型（含中心涡），当 $\gamma_N = \frac{(N-2)(N-6)}{12}$ 时，它同时满足 Hypotheses 3.1 和 3.3。

因此，结合 Theorem 3.4 与 Theorem 3.5 即得 Theorem 2.6。

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

a) 核心分析工具

（1）点涡系统的哈密顿结构与稳定性理论

论文充分利用了点涡动力学的辛结构（symplectic structure）。将 $Z = (z_1, \dots, z_N) \in (\mathbb{R}^2)^N$ 视为 $2N$ 维向量，引入质量矩阵 $\Gamma = \mathrm{diag}(\gamma_1, \gamma_1, \dots, \gamma_N, \gamma_N)$ 和标准辛矩阵 $J$，则 (PVS) 可写为：

$$ \Gamma \frac{d}{dt} Z = J \nabla H(Z), \quad H(Z) = \frac{1}{2} \sum_{i \ne j} \gamma_i \gamma_j \ln |z_i - z_j|. $$

在旋转坐标系 $Z(t) = e^{\nu t J} \tilde{Z}(t)$ 下，方程变为：

$$ \Gamma \frac{d}{dt} \tilde{Z} = J (\nabla H(\tilde{Z}) + \nu \Gamma \tilde{Z}), $$

相对平衡对应 $\nabla H(Z^*) + \nu \Gamma Z^* = 0$。围绕平衡态线性化后得到稳定性矩阵 $S = \Gamma^{-1} J D^2 H(Z^*) + \nu J$。

作者引用文献 [13] 的结果，将相空间分解为不变子空间 $V = \mathrm{span}\{Z^*, JZ^*\}$（对应平移与缩放模态）和正交补 $V^\perp$（物理扰动方向）。要求 $S|_{V^\perp}$ 只有纯虚特征值，即线性稳定性，防止扰动放大。

（2）中心涡量轨迹的误差估计（Proposition 4.1）

这是全文最关键的中间结果。考虑一个受摄动的点涡系统：

$$ \left\| \Gamma \frac{d}{dt} Z_\varepsilon - J \nabla H(Z_\varepsilon) - \nu J \Gamma Z_\varepsilon \right\| \leq C \varepsilon, $$

初始偏差 $|Z_\varepsilon(0) - Z^*| \leq \varepsilon$。命题断言：只要 $\alpha < \beta/2$，就有

$$ |Z_\varepsilon(t) - Z^*| \leq C \varepsilon^{2\beta - 2\alpha}, \quad \forall t \leq \varepsilon^{-\alpha}. $$

证明依赖于将扰动分解到各个不变子空间 $W_i \subset V^\perp$，利用 $S|_{W_i}$ 的反对称性（即能量守恒型结构）来控制增长率。例如，在 $W_i$ 分量 $Y_i$ 上：

$$ \frac{d}{dt} |Y_i|^2 = 2 Y_i \cdot \frac{dY_i}{dt} = 2 Y_i \cdot (S Y_i + O(|Z_\varepsilon - Z^*|^2) + O(\varepsilon)). $$

由于 $Y_i \cdot S Y_i = 0$，主导项消失，剩下高阶小量，从而导出 $dY_i/dt = O(\varepsilon^{2\beta})$，积分得 $|Y_i(t)| = O(\varepsilon^{2\beta - \alpha})$。

（3）外场 Lipschitz 估计（Lemma 5.1）

为了控制单个斑块的演化，需估计其它斑块产生的速度场 $F_i^\varepsilon(x,t)$ 在局部区域的光滑性。关键引理指出：当 $x,y \in D(z_i^*(t), \varepsilon^\beta)$ 时，

$$ |F_i^\varepsilon(x,t) - F_i^\varepsilon(y,t)| \leq C \varepsilon^\beta |x - y|. $$

这依赖于 Hypothesis 3.1 提供的“零应变”性质，使速度梯度在中心点为零，再结合支集小与强度固定的事实，得以获得比通常更好的Lipschitz常数。

（4）矩方法与能量型估计

引入每个斑块的中心涡量 $b_i^\varepsilon(t)$ 和归一化二阶矩：

$$ I_i^\varepsilon(t) = \frac{1}{\gamma_i} \int |x - b_i^\varepsilon(t)|^2 \omega_{\varepsilon,i}(x,t)\,dx, \quad I^\varepsilon(t) = \max_i I_i^\varepsilon(t). $$

利用输运方程可推导 $I_i^\varepsilon$ 的微分不等式：

$$ \frac{d}{dt} I_i^\varepsilon \lesssim \|DF_i^\varepsilon\|_{L^\infty} I_i^\varepsilon + o(\varepsilon^2), $$

再结合 $F_i^\varepsilon$ 的低应变特性，最终得出 $I^\varepsilon(t) \leq C \varepsilon^2$，即斑块尺寸平方仍为 $\varepsilon^2$ 阶。

b) 关键技术难点与处理技巧

（1）克服一般情形下的对数时间瓶颈

已知一般点涡构型只能保证 $\tau_{\varepsilon,\beta} \gtrsim |\ln \varepsilon|$（Buttà–Marchioro [3]），原因是即使无碰撞，斑块间相互作用会引起累积漂移。本文通过选取零一阶应变构型，使得每个斑块所在位置的速度场梯度为零，从而抑制了剪切效应，使斑块自身形变更慢。

（2）处理非线性耦合与误差传播

实际斑块的中心 $b_i^\varepsilon(t)$ 并不严格遵循点涡方程，而是满足带误差项的动力学：

$$ \frac{d}{dt} b_i^\varepsilon = F_i^\varepsilon(b_i^\varepsilon,t) + R_i^\varepsilon(t), \quad |R_i^\varepsilon| = O(\varepsilon^{2\beta}). $$

作者巧妙地将此视为对理想点涡系统的摄动，并应用 Proposition 4.1 来控制 $B^\varepsilon(t) = (b_1^\varepsilon, \dots, b_N^\varepsilon)$ 偏离 $Z^*(t)$ 的程度，从而闭环论证支集不逃逸。

（3）精确计算特殊构型的稳定性参数

第6节通过显式计算验证了当 $\gamma_N = \frac{(N-2)(N-6)}{12}$ 时，既满足零应变条件（H3.1），又落入[4]给出的线性稳定区间（2.5）。这一选择极为精细，例如当 $N=5$ 时 $\gamma_5 = -1$，$N=7$ 时 $\gamma_7 = 5/12 > 0$，体现了对称性与稳定性之间的微妙平衡。

4. 与之前工作的比较

a) 突破与改进之处

（1）摆脱对“快速分离”或“单涡”的依赖

此前工作如 Buttà–Marchioro [3] 中，获得 $\varepsilon^{-\alpha}$ 时间尺度的关键前提是点涡之间迅速远离（如自相似扩张构型），从而使相互作用减弱。类似地，[7] 的结果局限于单一涡斑（$N=1$）在特殊域中的情况。

本文的重大突破在于：即使所有点涡保持有限距离、共同旋转（无扩张），也能实现代数级长时间约束。这极大地扩展了适用范围，贴近真实湍流中观察到的稳定涡旋团。

（2）提出并利用“零一阶应变”新条件

以往研究中，“稳定性”多指线性谱稳定性或能量凸性。本文首次强调“局部流场梯度为零”的重要性，即 Hypothesis 3.1：

$$ \sum_{j \ne i} \gamma_j J_K(z_i^* - z_j^*) = 0. $$

这一条件确保每个涡斑处于“驻点”附近，不受剪切力作用，从而抑制内部变形。这是比单纯稳定性更强的几何条件，也是实现长期约束的核心机制。

（3）改进时间尺度依赖关系

本文所得时间下界为 $\varepsilon^{-\alpha}$，其中 $\alpha < \min(\beta/2, 2 - 4\beta)$。相比[3][7]中的类似估计，此处允许 $\beta < 1/2$（传统极限），且 $\alpha$ 可达约 $0.4$（当 $\beta \approx 0.4$ 时），优于仅得 $\alpha < \beta$ 的早期结果。

（4）发现正五边形的特殊地位

一个深刻观察是：当 $N=6$ 时，所需中心强度 $\gamma_6 = 0$，即正五边形本身已是最佳稳定构型。这说明在 $N=5$ 时，无需额外中心涡即可实现最强约束，暗示五重对称在二维涡旋组织中有独特优势，可能与自然界中常见五边形云系结构有关。

b) 对假设条件的放宽与强化

  • 未放宽初始数据正则性：仍要求 $\omega_\varepsilon \in L^\infty$ 且支集紧致，延续Yudovich框架。
  • 未降低维度或改变区域：仍在 $\mathbb{R}^2$ 上研究，未限于有界域。
  • 反而加强了构型要求：不是对任意稳定构型成立，而是仅对满足精确零应变条件的特定 $\gamma_N$ 成立。这是一种“牺牲普适性换取精度”的策略。

但正因如此，论文揭示了一类最优稳定涡旋晶体的存在，并为其动力学持久性提供了严格的数学解释，推动了从“存在性”到“鲁棒性”的理解深化。

3. Global weak solutions and incompressible limit of two-dimensional isentropic compressible magnetohydrodynamic equations with ripped density and large initial data

二维等熵可压缩磁流体动力学方程在撕裂密度和大初值下的全局弱解及不可压极限

作者: Shuai Wang, Guochun Wu, Xin Zhong

类别: math.AP

发布日期: 2025-10-01

arXiv ID: 2510.00812v2

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

详细分析

中文标题: 二维等熵可压缩磁流体动力学方程在撕裂密度和大初值下的全局弱解及不可压极限

研究对象和背景

本文研究的是二维全平面 $\mathbb{R}^2$ 上的等熵可压缩磁流体动力学（Magnetohydrodynamic, MHD）方程组，其数学形式如下：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho u) = 0, \\ (\rho u)_t + \mathrm{div}(\rho u \otimes u) + \nabla P = \mu \Delta u + (\mu + \lambda) \nabla \mathrm{div}\, u + B \cdot \nabla B - \frac{1}{2} \nabla |B|^2, \\ B_t + u \cdot \nabla B - B \cdot \nabla u + B \mathrm{div}\, u = \nu \Delta B, \\ \mathrm{div}\, B = 0, \end{cases} \quad (x,t) \in \mathbb{R}^2 \times (0,+\infty), \tag{1.1} $$

其中未知函数分别为：密度 $\rho = \rho(x,t) \geq 0$、速度场 $u = (u_1,u_2)(x,t)$、磁场 $B = (B_1,B_2)(x,t)$。压力由状态方程给出：$P = a \rho^\gamma$，其中常数 $a > 0$，绝热指数 $\gamma > 1$。粘性系数满足物理条件：$\mu > 0$，$\mu + \lambda \geq 0$；电阻率 $\nu > 0$。

该系统描述了导电流体在电磁场作用下的运动规律，耦合了可压缩 Navier-Stokes 方程与 Maxwell 方程的简化模型，在天体物理、地球物理和等离子体物理中有广泛应用。

考虑以下 Cauchy 问题 的初始数据：

$$ (\rho, u, B)|_{t=0} = (\rho_0, u_0, B_0)(x), \quad x \in \mathbb{R}^2, \tag{1.2} $$

并假设远场行为为：

$$ (\rho_0, u_0, B_0)(x) \to (\tilde{\rho}, 0, 0), \quad \text{当 } |x| \to \infty, \tag{1.3} $$

其中参考密度 $\tilde{\rho} > 0$ 是一个正常数。

值得注意的是，本文允许初始密度具有“撕裂（ripped）结构”，即允许存在非平凡的真空区域（interior vacuum），且对密度没有额外正则性要求（如 Lipschitz 或 BV 正则）。这使得分析更具挑战性，因为方程在真空附近退化或奇异。

此外，初始能量可以是任意大的，不依赖于小扰动框架。这是对传统小能量理论的重大突破。

文中定义了总能量泛函：

$$ \int_{\mathbb{R}^2} \left( \frac{1}{2} \rho |u|^2 + \frac{1}{2} |B|^2 + G(\rho) \right) dx + \int_0^T \int_{\mathbb{R}^2} \left[ \mu |\nabla u|^2 + (\mu + \lambda)(\mathrm{div}\, u)^2 + \nu |\nabla B|^2 \right] dxd\tau \leq C_0, $$

其中势能项 $G(\rho)$ 定义为：

$$ G(\rho) \triangleq \int_{\tilde{\rho}}^\rho \frac{P(v) - P(\tilde{\rho})}{v^2} dv. $$

从能量估计可见，当体积粘性系数 $\lambda \to \infty$ 时，$\mathrm{div}\, u$ 被强烈压制，从而启发我们研究所谓的不可压极限：即当 $\lambda \to \infty$ 时，原可压缩 MHD 解趋于某个不可压 MHD 系统的解。

目标极限系统是非齐次不可压磁流体力学方程：

$$ \begin{cases} \varrho_t + v \cdot \nabla \varrho = 0, \\ \varrho v_t + \varrho v \cdot \nabla v + \nabla \Pi = \mu \Delta v + b \cdot \nabla b - \frac{1}{2} \nabla |b|^2, \\ b_t + v \cdot \nabla b = b \cdot \nabla v + \nu \Delta b, \\ \mathrm{div}\, v = \mathrm{div}\, b = 0, \end{cases} \tag{1.5} $$

初值为 $(\varrho, v, b)|_{t=0} = (\rho_0, v_0, B_0)$，其中 $v_0 = \mathcal{P}u_0$ 是 $u_0$ 的 Leray-Helmholtz 投影（无散部分）。

研究区域为全平面 $\mathbb{R}^2$，边界条件为空间无穷远处趋于常态（far-field condition），而非周期或有界域中的 Dirichlet/Neumann 条件。这一设定更接近真实物理场景，但技术上更为复杂，尤其在使用紧性方法和插值不等式时需特别处理。

函数空间方面，作者关注的是满足特定正则性的弱解，属于 Hoff 类型的中间正则性解（intermediate weak solutions），其特征在于：

  • 密度可有跳跃或真空；
  • 动量 $\rho u$ 和磁场 $B$ 具有一定 Sobolev 正则性；
  • 使用有效粘性通量 $F$、涡量 $\omega$ 等结构变量进行精细估计。

这类解比 Lions-Feireisl 型有限能量弱解更强（路径可定义），但弱于经典强解（允许密度不连续），适用于处理大初值和真空情形。

主要定理或主要结果

定理 1.1（全局弱解的存在性）

设初始数据 $(\rho_0, u_0, B_0)$ 满足：

  1. 有界密度：$0 \leq \inf \rho_0 \leq \sup \rho_0 \leq \hat{\rho}$，
  2. 初始能量控制：
    $$ C_0 + \mu \| \nabla u_0 \|^2_{L^2} + (\mu + \lambda) \| \mathrm{div}\, u_0 \|^2_{L^2} + \nu \| \nabla B_0 \|^2_{L^2} \leq M, $$
    其中 $C_0$ 为初始总能量，$M > 0$ 不一定小；
  3. $\mathrm{div}\, B_0 = 0$ in $\mathcal{D}'(\mathbb{R}^2)$。

则存在仅依赖于 $\tilde{\rho}, \hat{\rho}, a, \gamma, \nu, \mu$ 的正常数 $D$，使得当体积粘性系数 $\lambda$ 满足：

$$ \lambda \geq \exp\left\{ (2+M)\exp\left[ D(1+C_0)^6 \right] \right\}, \tag{1.11} $$

时，Cauchy 问题 (1.1)–(1.3) 存在一个全局弱解 $(\rho, u, B)$，满足 Definition 1.1，并具备如下正则性：

$$ \begin{cases} 0 \leq \rho(x,t) \leq 2\hat{\rho}, \quad \text{a.e. on } \mathbb{R}^2 \times [0,\infty), \\ (\rho - \tilde{\rho}, \sqrt{\rho} u, B) \in C([0,\infty); L^2(\mathbb{R}^2)), \\ (u, B) \in L^\infty(0,\infty; H^1(\mathbb{R}^2)), \\ (\nabla^2 P, \nabla F, \sqrt{\rho} \dot{u}, B_t, \nabla^2 B) \in L^2(\mathbb{R}^2 \times (0,\infty)), \\ \sigma^{1/2} (\sqrt{\rho} \dot{u}, B_t, \nabla^2 B) \in L^\infty(0,\infty; L^2(\mathbb{R}^2)), \\ \sigma^{1/2} (\nabla \dot{u}, \nabla B_t) \in L^2(\mathbb{R}^2 \times (0,\infty)), \end{cases} \tag{1.12} $$

其中 $\sigma = \min\{1,t\}$，$\dot{f} = f_t + u \cdot \nabla f$ 表示物质导数。

关键点：即使初始能量 $C_0$ 和 $M$ 可以很大，只要 $\lambda$ 足够大，仍可构造全局弱解。且允许真空存在。

定理 1.2（不可压极限）

令 $\{(\rho_\lambda, u_\lambda, B_\lambda)\}$ 为定理 1.1 所得的一族解。则存在子列 $\{\lambda_k\} \to \infty$，使得当 $k \to \infty$ 时，

  • $\rho_{\lambda_k} \to \varrho$ 在 $L^2(\mathbb{R}^2)$ 中关于时间一致收敛；
  • $u_{\lambda_k} \to v$, $B_{\lambda_k} \to b$ 在 $\mathbb{R}^2 \times (0,\infty)$ 的任何紧集上一致收敛；

且极限 $(\varrho, v, b)$ 是非齐次不可压 MHD 方程 (1.5) 的全局弱解（Definition 1.2），满足：

$$ \begin{aligned} & \varrho \in L^\infty(\mathbb{R}^2 \times (0,\infty)), \quad (\sqrt{\varrho} v, b) \in L^\infty([0,\infty); L^2(\mathbb{R}^2)), \\ & (\nabla v, \nabla b) \in L^2(\mathbb{R}^2 \times (0,T)) \quad \forall T>0, \\ & (\varrho - \tilde{\rho}) \in C([0,\infty); L^2(\mathbb{R}^2)), \quad v \in L^2(\mathbb{R}^2 \times (0,T)). \end{aligned} \tag{1.14}–(1.15) $$

意义：首次建立了二维可压缩 MHD 在大体积粘性下的不可压极限，推广了 Navier-Stokes 情形的结果。

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

a) 核心分析工具

  1. 有效粘性通量（Effective Viscous Flux）：
    定义 $F \triangleq (2\mu + \lambda)\mathrm{div}\, u - (P(\rho) - P(\tilde{\rho})) - \frac{1}{2}|B|^2$。它是连接速度散度与压力的关键桥梁。利用椭圆估计：

    $$ \Delta F = \mathrm{div}(\rho \dot{u} - B \cdot \nabla B), $$

    可将 $\nabla u$ 分解为 $F$ 和涡量 $\omega = \mathrm{curl}\, u$ 的组合，从而获得高阶可积性估计（见 Lemma 2.5）。

  2. Desjardins 型对数插值不等式：
    关键引理（Lemma 2.3）指出：若 $0 \leq \rho \leq \hat{\rho}$，$u \in H^1$，则

    $$ \| \sqrt{\rho} u \|_{L^4}^2 \leq C(\hat{\rho})(1 + \| \sqrt{\rho} u \|_{L^2}) \| u \|_{H^1} \sqrt{ \ln\left(2 + \| u \|_{H^1}^2 \right) }. $$

    这个不等式允许我们在没有真空附近良好控制的情况下，通过 $\| u \|_{H^1}$ 控制加权 $L^4$ 范数，对于克服真空带来的奇异性至关重要。

  3. Gagliardo-Nirenberg 不等式：用于嵌入估计，例如：

    $$ \| f \|_{L^p} \leq C \| f \|_{L^2}^{2/p} \| \nabla f \|_{L^2}^{1 - 2/p}, \quad p \in [2,\infty). $$

  4. Leray-Helmholtz 投影：将向量场分解为无散（$\mathcal{P}$）和梯度（$\mathcal{Q}$）部分，便于分离声波模态与不可压模态，在极限过程中提取收敛子列。

  5. 紧性论证（Compactness Arguments）：基于能量估计和 Aubin-Lions 引理等，从近似解序列中提取强收敛子列，最终通过弱下半连续性和 mollifier 收敛性证明极限满足弱解定义。

  6. Lagrangian 坐标技巧：用于控制密度上界。通过追踪粒子轨迹，结合 ODE 分析，避免直接使用最大值原理。

b) 技术难点与应对策略

  1. 密度上界的统一估计：
    最大挑战之一是证明 $\|\rho\|_{L^\infty}$ 不爆破。作者引入先验假设：

    $$ \sup_{R^2 \times [0,T]} \rho \leq 2\hat{\rho}, \quad \frac{1}{(2\mu + \lambda)^2} \int_0^T \| P(\rho) - P(\tilde{\rho}) \|_{L^4}^4 dt \leq 2, $$

    并试图闭合此假设。关键是通过估计 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 来控制 $\|F\|_{L^4}$，再借助 Desjardins 不等式得到 $\|\nabla u\|_{L^4}$，从而用 Gronwall 方法完成闭合。

  2. 材料导数 $\dot{u}$ 的估计困难：
    因 $\lambda$ 出现在 $\mathrm{div}\, u$ 项中，难以直接控制 $\|\dot{u}\|_{L^\infty(0,T;L^2)}$。作者借鉴 [13] 的思想，转而估计：

    $$ \frac{d}{dt} \int (\mu + \lambda)(\mathrm{div}\, u)^2 dx, $$

    而非 $\int (\mu + \lambda)(\mathrm{div}\, \dot{u})^2 dx$，避免出现高阶时间导数。同时利用：

    $$ \mathrm{div}\, \dot{u} = \mathrm{div}\, u_t + u \cdot \nabla \mathrm{div}\, u - u_i \partial_j u_j, $$

    将非线性项转化为可控形式。

  3. 磁场耦合项的处理：
    如 $B \cdot \nabla B$、$B_t$、$u \cdot \nabla B$ 等非线性项带来额外复杂性。作者通过：

    • 从 $B$ 方程导出 $\|B_t\|_{L^2}$ 和 $\|\nabla^2 B\|_{L^2}$ 的估计；
    • 利用 $\|B\|_{L^4} \leq C \|B\|_{L^2}^{1/2} \|\nabla B\|_{L^2}^{1/2}$ 控制交叉项；
    • 在估计 $\rho \dot{u}$ 时，将 $B \cdot \nabla B$ 视为外力项参与能量估计。

  4. 真空与大初值共存的挑战：
    当初值含真空且能量大时，传统的 Sobolev 嵌入失效。解决方式是：

    • 使用加权范数 $\|\sqrt{\rho} u\|$ 替代 $\|u\|$；
    • 结合 Desjardins 不等式提升可积性；
    • 利用远场非真空条件（$\rho \to \tilde{\rho} > 0$）激活对数不等式。

  5. 不可压极限的紧性获取：
    通过均匀估计（uniform in $\lambda$）获得：

    • $\{u_\lambda\}$ 在 $L^2(0,T;H^1)$ 中有界；
    • $\{\partial_t (\rho_\lambda u_\lambda)\}$ 在负阶 Sobolev 空间有界； 再结合 Arzelà-Ascoli 和 Aubin-Lions 引理，得到强收敛性，进而验证极限满足弱解方程。

与之前工作的比较

本文相较于已有工作有多重显著改进：

  1. 放宽初始能量限制：
    相比于 Suen-Hoff [26]、Li-Yu [16] 等要求“小能量”的结果，本文允许任意大的初始能量，仅需 $\lambda$ 足够大即可保证全局存在。这是本质性突破。

  2. 允许内部真空且无需零动量假设：
    Danchin-Mucha [2] 在环面 $T^2$ 上研究 Navier-Stokes 时虽允许真空，但要求初始总动量为零。本文在 $\mathbb{R}^2$ 上取消了这一限制，适用范围更广。

  3. 推广至 MHD 系统：
    文献 [1,2] 仅针对 Navier-Stokes 方程。本文成功将结果推广到强耦合的 MHD 系统，需处理 $u \cdot \nabla B$、$B \cdot \nabla u$、$B \cdot \nabla B$ 等新非线性项，技术难度更高。

  4. 首次建立二维可压缩 MHD 的不可压极限：
    尽管 Lions-Feireisl 理论在三维有不可压极限，但在二维大初值情形下尚属空白。本文首次证明当 $\lambda \to \infty$ 时，解收敛到非齐次不可压 MHD 解，具有开创意义。

  5. 摆脱有界域限制：
    [2] 的方法依赖于环面上的 Poincaré 不等式和紧嵌入，无法直接推广到 $\mathbb{R}^2$。本文通过远场非真空条件和对数不等式克服此障碍。

  6. 更高的解正则性：
    得到的弱解属于 Hoff 类，比 Lions 型弱解具有更强的正则性（如 $\nabla^2 B \in L^2$），可用于进一步研究唯一性或稳定性。

综上所述，本文在大初值、真空、全平面、强耦合 MHD 等多重困难下，利用有效粘性通量 + 对数插值 + Lagrangian 技巧，实现了全局存在性和不可压极限的双重突破，推动了可压缩流体力学理论的发展。

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

1. A Coupled Generalized Korteweg-de Vries System Driven by White Noise

受白噪声驱动的耦合广义Korteweg-de Vries系统

作者: Aissa Boukarou

类别: math.AP

发布日期: 2025-10-27

arXiv ID: 2510.23843v2

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

摘要翻译

中文标题: 受白噪声驱动的耦合广义Korteweg-de Vries系统

摘要翻译: 本文研究了受白噪声驱动的耦合广义Korteweg-de Vries系统的Cauchy问题。我们证明了当初始数据属于 $ H^{s} \times H^{s} $ 且 $ s>1/2 $ 时，该问题的局部适定性。本文所使用的关键工具包括Bourgain空间中的多线性估计、It^o公式以及不动点论证。我们的结果改进了Gomes和Pastor \cite{gomes2021solitary} 所得到的局部适定性结论。

2. Improved Gevrey Class Regularity of the Kadomtsev Petviashvili Equation

五阶Kadomtsev-Petviashvili方程解的Gevrey类正则性改进

作者: Aissa Boukarou, Lamia Seghour

类别: math.AP

发布日期: 2025-10-30

arXiv ID: 2510.26669v1

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

摘要翻译

中文标题: 五阶Kadomtsev-Petviashvili方程解的Gevrey类正则性改进

摘要翻译: 本文改进并推广了Boukarou等人 \cite{boukarou1} 关于一类五阶Kadomtsev-Petviashvili（KP）型方程解的Gevrey正则性所得到的结果。我们建立了该方程在 $2+1$ 维情形下解关于时间变量的Gevrey正则性，通过一种新的分析方法获得了更精确的结果。假设初始数据在空间变量中具有阶数 $\sigma \geq 1$ 的Gevrey正则性，我们证明了相应的解在时间变量中具有阶数为 $5\sigma$ 的Gevrey正则性。此外，我们还证明了函数 $u(x, y, t)$ 作为 $t$ 的函数，不属于任何满足 $1 \leq z < 5\sigma$ 的Gevrey类 $G^z$。证明过程同时处理了变量 $x$、$y$ 和 $t$，采用了控制级数（majorant series）方法，并精确追踪了高阶色散项 $\partial_x^5 u$ 与低阶项 $\alpha \partial_x^3 u$、$\partial_x^{-1} \partial_y^2 u$ 以及 $u \partial_x u$ 的共同影响。

3. Time-periodic boundary effects on the shocks for scalar conservation laws

时间周期性边界条件对标量守恒律中激波的影响

作者: Yuan Yuan

类别: math.AP

发布日期: 2025-10-30

arXiv ID: 2510.26153v1

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

摘要翻译

中文标题: 时间周期性边界条件对标量守恒律中激波的影响

摘要翻译: 本文研究了定义在半直线 $(-\infty,0)$ 上的标量守恒律方程在激波速度 $s<0$ 情形下，无粘与有粘激波在时间周期性边界条件下的渐近稳定性问题，此类边界条件源于流体力学中的经典活塞问题。尽管该问题具有重要意义，但时间周期性边界条件如何影响黎曼解的长时间行为至今尚不清楚。本工作通过严格的数学分析填补了这一空白，证明了在无粘和有粘两种情形下，解的渐近状态不仅由平移后的背景（有粘）激波所决定，还与由时间周期性边界条件诱导出的时间周期性边界解相耦合。我们的分析表明，这种耦合效应表现为一种传播的“边界波”，并对激波动力学产生显著影响。

4. A linear, unconditionally stable, second order decoupled method for the Ericksen-Leslie model with SAV approach

基于SAV方法的Ericksen-Leslie模型的一种线性、无条件稳定、二阶解耦格式

作者: Ruonan Cao, Nianyu Yi

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

发布日期: 2025-03-25

arXiv ID: 2503.19424v3

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

摘要翻译

中文标题: 基于SAV方法的Ericksen-Leslie模型的一种线性、无条件稳定、二阶解耦格式

摘要翻译: 本文提出了一种求解Ericksen-Leslie模型的二阶精度、线性、完全解耦且无条件能量稳定的数值格式。该方法结合了压力修正法与标量辅助变量（SAV）技术。我们严格证明了所提出格式的无条件能量稳定性。此外，通过多个数值实验验证了该格式的收敛阶数、稳定性以及计算效率。

不相关论文（基本信息）

1. Quantitative Lorentzian isoperimetric inequalities

定量洛伦茨等周不等式

作者: Christian Lange, Jonas W. Peteranderl

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

发布日期: 2025-10-30

arXiv ID: 2510.26755v1

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

摘要: We establish optimal stability estimates in terms of the Fraenkel asymmetry with universal dimensional constants for a Lorentzian isoperimetric inequality due to Bahn and Ehrlich and, as a consequence, for a special version of a Lorentzian isoperimetric inequality due to Cavalletti and Mondino. For the Bahn–Ehrlich inequality the Fraenkel asymmetry enters the stability result quadratically like in the Euclidean case while for the Cavalletti–Mondino inequality the Fraenkel asymmetry enters linearly. As it turns out, refining the latter inequality through an additional geometric term allows us to recover the more common quadratic stability behavior. Along the way, we provide simple self-contained proofs for the above isoperimetric-type inequalities.

2. Degenerate singular Kirchhoff problems in Musielak-Orlicz spaces

Musielak-Orlicz空间中退化的奇异Kirchhoff问题

作者: Umberto Guarnotta, Patrick Winkert

类别: math.AP

发布日期: 2024-11-15

arXiv ID: 2411.09970v2

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

摘要: In this paper we study quasilinear elliptic Kirchhoff equations driven by a non-homogeneous operator with unbalanced growth and right-hand sides that consist of sub-linear, possibly singular, and super-linear reaction terms. Under very general assumptions we prove the existence of at least two solutions for such problems by using the fibering method along with an appropriate splitting of the associated Nehari manifold. In contrast to other works our treatment is very general, with much easier and shorter proofs as it was done in the literature before. Furthermore, the results presented in this paper cover a large class of second-order differential operators like the $p$-Laplacian, the $(p,q)$-Laplacian, the double phase operator, and the logarithmic double phase operator.

3. Homogenization of elasto-plastic plate equations with vanishing hardening

硬化效应消失的弹塑性板方程的均匀化

作者: Marin Bužančić, Igor Velčić, Josip Žubrinić

类别: math.AP, 74C05, 74G65, 74K20, 49J45, 74Q09, 35B27

发布日期: 2025-06-09

arXiv ID: 2506.08215v2

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

摘要: We study the asymptotic behavior of thin heterogeneous elastoplastic plates in the framework of linearized elastoplasticity, focusing on the regime where the plate thickness vanishes much faster than the characteristic scale of the material’s periodic microstructure. In contrast to earlier analyzes that required restrictive geometric assumptions on admissible yield surfaces, our approach accommodates general relations between phases without imposing any specific ordering. The analysis proceeds in two main steps. First, we rigorously derive a heterogeneous plate model with both isotropic and kinematic hardening through a dimension reduction procedure based on evolutionary $\Gamma$-convergence. This result extends existing plate models for homogeneous materials to the heterogeneous setting and allows for general forms of hardening and dissipation potentials. In the second step, we perform two-scale homogenization while simultaneously letting the hardening tend to zero. This process yields an effective elasto-perfectly plastic plate model and, crucially, provides a characterization of the dissipation potential at the interfaces between different phases. The resulting dissipation functional takes the form of a non-local inf-convolution of the traces of plastic strains on both sides of the interface, reflecting the Kirchhoff-Love structure of admissible displacements.

4. A Convolutional Hierarchical Deep-learning Neural Network (C-HiDeNN) Framework for Non-linear Finite Element Analysis

用于非线性有限元分析的卷积分层深度学习神经网络（C-HiDeNN）框架

作者: Yingjian Liu, Monish Yadav Pabbala, Jiachen Guo, Chanwook Park, Gino Domel, Wing Kam Liu, Dong Qian

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

发布日期: 2025-09-02

arXiv ID: 2509.02435v2

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

摘要: We present a framework for the Convolutional Hierarchical Deep-learning Neural Network (C-HiDeNN) tailored for nonlinear finite element analysis. Building upon the structured foundation of HiDeNN, C-HiDeNN introduces a convolution operator to enhance numerical approximation. A distinctive feature of C-HiDeNN is its higher-order accurate approximation achieved through an expanded set of parameters, such as the polynomial order ‘p,’ dilation parameter ‘a,’ patch size ’s,’ and nodal position ‘X’. These parameters function as the functional equivalents of weights and biases within each C-HiDeNN patch. In addition, C-HiDeNN can be selectively applied to regions requiring high resolution to adaptively improve local prediction accuracy. To demonstrate the effectiveness of this framework, we provide numerical examples in the context of nonlinear finite element analysis. The results show that our approach achieves significantly higher accuracy than conventional Finite Element Method (FEM) while substantially reducing computational costs.

5. On a parabolic curvature lower bound generalizing Ricci flows

关于推广里奇流的抛物曲率下界的研究

作者: Marco Flaim, Erik Hupp

类别: math.DG, math.AP

发布日期: 2025-07-07

arXiv ID: 2507.05032v2

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

摘要: Optimal transport plays a major role in the study of manifolds with Ricci curvature bounded below. Some results in this setting have been extended to super Ricci flows, revealing a unified approach to analysis on Ricci nonnegative manifolds and Ricci flows. However we observe that the monotonicity of Perelman’s functionals ($\mathcal{F}$, $\mathcal{W}$, reduced volume), which hold true for Ricci flows and Ricci nonnegative manifolds, cannot be strictly generalized to super Ricci flows. In 2010 Buzano introduced a condition which still generalizes Ricci flows and Ricci nonnegative manifolds, and on which Perelman’s monotonicities do hold. We provide characterizations of this condition using optimal transport and understand it heuristically as Ricci nonnegativity of the space-time. This interpretation is consistent with its equivalence to Ricci nonnegativity on Perelman’s infinite dimensional manifold. More precisely, we prove that for smooth evolutions of Riemannian manifolds, this condition is equivalent to a Bochner inequality (resembling Perelman’s Harnack inequality but for the forward heat flow), a gradient estimate for the heat flow, a Wasserstein contraction along the adjoint heat flow, the convexity of a modified entropy along Wasserstein geodesics, and an Evolutionary Variational Inequality (EVI). The optimal transport statements use Perelman’s $L$ distance as cost, as first studied on Ricci flows by Topping and by Lott. We also consider dimensionally improved and weighted versions of these conditions. The dimensional Bochner inequality and all gradient estimates for the forward heat equation, along with the EVIs, appear to be new even for general Ricci flows, and are related to the Hamiltonian perspective on the $L$ distance. Most of our proofs do not use tensor calculus or Jacobi fields, suggesting the possibility of future extensions to more singular settings.

6. On a semi-discrete model of Maxwell’s equations in three and two dimensions

关于三维与二维情况下麦克斯韦方程组的半离散模型

作者: Volodymyr Sushch

类别: math-ph, cs.NA, math.AP, math.MP, math.NA, 39A12, 39A70, 35Q61

发布日期: 2025-10-30

arXiv ID: 2510.26427v1

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

摘要: In this paper, we develop a geometric, structure-preserving semi-discrete formulation of Maxwell’s equations in both three- and two-dimensional settings within the framework of discrete exterior calculus. This approach preserves the intrinsic geometric and topological structures of the continuous theory while providing a consistent spatial discretization. We analyze the essential properties of the proposed semi-discrete model and compare them with those of the classical Maxwell’s equations. As a special case, the model is illustrated on a combinatorial two-dimensional torus, where the semi-discrete Maxwell’s equations take the form of a system of first-order linear ordinary differential equations. An explicit expression for the general solution of this system is also derived.

7. Tangential approach in the Dirichlet problem for elliptic equations

椭圆型方程狄利克雷问题中的切向逼近方法

作者: Jonathan Bennett, Arnaud Dumont, Andrew J. Morris

类别: math.AP, 35J25 (Primary) 35J25, 35B30, 42B37 (Secondary)

发布日期: 2025-10-30

arXiv ID: 2510.26400v1

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

摘要: It is well-known that solvability of the $\mathrm{L}^{p}$-Dirichlet problem for elliptic equations $Lu:=-\mathrm{div}(A\nabla u)=0$ with real-valued, bounded and measurable coefficients $A$ on Lipschitz domains $\Omega\subset\mathbb{R}^{1+n}$ is characterised by a quantitative absolute continuity of the associated $L$-harmonic measure. We prove that this local $A_{\infty}$ property is sufficient to guarantee that the nontangential convergence afforded to $\mathrm{L}^{p}$ boundary data actually improves to a certain \emph{tangential} convergence when the data has additional (Sobolev) regularity. Moreover, we obtain sharp estimates on the Hausdorff dimension of the set on which such convergence can fail. This extends results obtained by Dorronsoro, Nagel, Rudin, Shapiro and Stein for classical harmonic functions in the upper half-space.

8. Asymptotic meshes from $r$-variational adaptation methods for static problems in one dimension

一维静态问题中基于$r$-变分自适应方法的渐近网格

作者: Darith Hun, Nicolas Moës, Heiner Olbermann

类别: math.NA, cs.NA, math.AP, 49J45, 49Q10, 65N30, 65N50, 74S05

发布日期: 2025-10-30

arXiv ID: 2510.26375v1

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

摘要: We consider the minimization of integral functionals in one dimension and their approximation by $r$-adaptive finite elements. Including the grid of the FEM approximation as a variable in the minimization, we are able to show that the optimal grid configurations have a well-defined limit when the number of nodes in the grid is being sent to infinity. This is done by showing that the suitably renormalized energy functionals possess a limit in the sense of $\Gamma$-convergence. We provide numerical examples showing the closeness of the optimal asymptotic mesh obtained as a minimizer of the $\Gamma$-limit to the optimal finite meshes.

9. Simulation of the magnetic Ginzburg-Landau equation via vortex tracking

通过涡旋追踪模拟磁性Ginzburg-Landau方程

作者: Thiago Carvalho Corso, Gaspard Kemlin, Christof Melcher, Benjamin Stamm

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

发布日期: 2025-10-30

arXiv ID: 2510.26334v1

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

摘要: This paper deals with the numerical simulation of the 2D magnetic time-dependent Ginzburg-Landau (TDGL) equations in the regime of small but finite (inverse) Ginzburg-Landau parameter $\epsilon$ and constant (order $1$ in $\epsilon$) applied magnetic field. In this regime, a well-known feature of the TDGL equation is the appearance of quantized vortices with core size of order $\epsilon$. Moreover, in the singular limit $\epsilon \searrow 0$, these vortices evolve according to an explicit ODE system. In this work, we first introduce a new numerical method for the numerical integration of this limiting ODE system, which requires to solve a linear second order PDE at each time step. We also provide a rigorous theoretical justification for this method that applies to a general class of 2D domains. We then develop and analyze a numerical strategy based on the finite-dimensional ODE system to efficiently simulate the infinite-dimensional TDGL equations in the presence of a constant external magnetic field and for small, but finite, $\epsilon$. This method allows us to avoid resolving the $\epsilon$-scale when solving the TDGL equations, where small values of $\epsilon$ typically require very fine meshes and time steps. We provide numerical examples on a few test cases and justify the accuracy of the method with numerical investigations.

10. Numerical simulations of oscillations for axisymmetric solar backgrounds with differential rotation and gravity

具有较差自转和引力的轴对称太阳背景下振荡的数值模拟

作者: Ha Pham, Florian Faucher, Damien Fournier, Hélène Barucq, Laurent Gizon

类别: astro-ph.SR, astro-ph.IM, math.AP

发布日期: 2025-08-26

arXiv ID: 2508.19386v5

链接: http://arxiv.org/abs/2508.19386v5

摘要: Local helioseismology comprises of imaging and inversion techniques employed to reconstruct the dynamic and interior of the Sun from correlations of oscillations observed on the surface, all of which require modeling solar oscillations and computing Green’s kernels. In this context, we implement and investigate the robustness of the Hybridizable Discontinuous Galerkin (HDG) method in solving the equation modeling stellar oscillations for realistic solar backgrounds containing gravity and differential rotation. While a common choice for modeling stellar oscillations is the Galbrun’s equation, our working equations are derived from an equivalent variant, involving less regularity in its coefficients, working with Lagrangian displacement and pressure perturbation as unknowns. Under differential rotation and axisymmetric assumption, the system is solved in azimuthal decomposition with the HDG method. Compared to no-gravity approximations, the mathematical nature of the wave operator is now linked to the profile of the solar buoyancy frequency N which encodes gravity, and leads to distinction into regions of elliptic or hyperbolic behavior of the wave operator at zero attenuation. While small attenuation is systematically included to guarantee theoretical well-posedness, the above phenomenon affects the numerical solutions in terms of amplitude and oscillation pattern, and requires a judicious choice of stabilization. We investigate the stabilization of the HDG discretization scheme, and demonstrate its importance to ensure the accuracy of numerical results, which is shown to depend on frequencies relative to N, and on the position of the Dirac source. As validations, the numerical power spectra reproduce accurately the observed effects of the solar rotation on acoustic waves.

11. Complete spectrum of the Robin eigenvalue problem on the ball

球域上Robin特征值问题的完备谱

作者: Cancan Chen, Guowei Dai, Yingxin Sun

类别: math.AP, math.SP, 34L20, 34L10, 35J05, 35P10, 35P20

发布日期: 2025-10-30

arXiv ID: 2510.26331v1

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

摘要: We investigate the following Robin eigenvalue problem \begin{equation*} \left{ \begin{array}{ll} -\Delta u=\mu u,, &\text{in},, B,\ \partial_\texttt{n} u+\alpha u=0 &\text{on},, \partial B \end{array} \right. \end{equation*} on the unit ball of $\mathbb{R}^N$. We obtain the complete spectral structure of this problem. In particular, for $\alpha>0$, we find that the first eigenvalue is $k_{\nu,1}^2$ and the second eigenvalue is exactly $k_{\nu+1,1}^2$, where $k_{\nu+l,m}$ is the $m$th positive zero of $kJ_{\nu+l+1}(k)-(\alpha+l) J_{\nu+l}(k)$. Moreover, when $\alpha\in(-1,0)$, the first eigenvalue is $-\widehat{k}_{\nu,1}^2$ where $\widehat{k}_{\nu,1}$ denotes the unique zero of $\alpha I_{\nu}(k)+kI_{\nu+1}(k)$, and the second eigenvalue is exactly $k_{\nu+1,1}^2$. Furthermore, for $\alpha=-1$, the first eigenvalue is $-\widehat{k}_{\nu,1}^2$ and the second eigenvalue is exactly $0$. Our conclusions indicate the ratio $\mu_2/\mu_1$ may be positive, negative or zero according to the suitable ranges of the parameter $\alpha$.

12. Coupling local and nonlocal total variation flow for image despeckling

耦合局部与非局部全变分流的图像去斑方法

作者: Yi Ran, Zhichang Guo, Kehan Shi, Qirui Zhou, Jingfeng Shao, Martin Burger, Boying Wu

类别: math.AP

发布日期: 2025-10-30

arXiv ID: 2510.26296v1

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

摘要: Nonlocal equations effectively preserve textures but exhibit weak regularization effects in image denoising, whereas local equations offer strong denoising capabilities yet fail to protect textures. To integrate the advantages of both approaches, this paper investigates a coupled local-nonlocal total variation flow for image despeckling. We establish the existence and uniqueness of the weak solution for the proposed equation. Several properties, including the equivalent forms of the weak solution and its asymptotic behavior, are derived. Furthermore, we demonstrate that the weak solutions of the proposed equation converge to the weak solution of the classical total variation flow under kernel rescaling. The importance of coupling is highlighted through comparisons with local and nonlocal models for image despeckling.

13. Sharp embeddings and existence results for Logarithmic $p$-Laplacian equations with critical growth

具有临界增长的对数型 $p$-Laplace 方程的精确嵌入与存在性结果

作者: Rakesh Arora, Jacques Giacomoni, Hichem Hajaiej, Arshi Vaishnavi

类别: math.AP

发布日期: 2025-10-30

arXiv ID: 2510.26286v1

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

摘要: In this paper, we derive a new $p$-Logarithmic Sobolev inequality and optimal continuous and compact embeddings into Orlicz-type spaces of the function space associated with the logarithmic $p$-Laplacian. As an application of these results, we study a class of Dirichlet boundary value problems involving the logarithmic $p$-Laplacian and critical growth nonlinearities perturbed with superlinear-subcritical growth terms. By employing the method of the Nehari manifold, we prove the existence of a nontrivial weak solution. Lastly, we conduct an asymptotic analysis of a weighted nonlocal, nonlinear problem governed by the fractional $p$-Laplacian with superlinear or sublinear type non-linearity, demonstrating the convergence of least energy solutions to a non-trivial, non-negative least energy solution of a Brezis-Nirenberg type or logistic-type problem, respectively, involving the logarithmic $p$-Laplacian as the fractional parameter $s \to 0^+$. The findings in this work serve as a nonlinear analogue of the results reported in \cite{Angeles-Saldana, Arora-Giacomoni-Vaishnavi, Santamaria-Saldana}, thereby extending their scope to a broader variational framework.

14. Quantization property of n-Laplacian mean field equation and sharp Moser-Onofri inequality

n-拉普拉斯平均场方程的量子化性质与精确的Moser-Onofri不等式

作者: Lu Chen, Guozhen Lu, Bohan Wang

类别: math.AP

发布日期: 2024-06-02

arXiv ID: 2406.00743v3

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

摘要: In this paper, we are concerned with the following $n$-Laplacian mean field equation

$$ \left\{ {\begin{array}{*{20}{c}} { - \Delta_n u = \lambda e^u} & {\rm in} \ \ \Omega, \\ {\ \ \ \ u = 0} &\ {\rm on}\ \partial \Omega, \end{array}} \right. $$

[] where $\Omega$ is a smooth bounded domain of $\mathbb{R}^n \ (n\geq 2)$ and $- \Delta_n u =-{\rm div}(|\nabla u|^{n-2}\nabla u)$. We first establish the quantization property of solutions to the above $n$-Laplacian mean field equation. As an application, combining the Pohozaev identity and the capacity estimate, we obtain the sharp constant $C(n)$ of the Moser-Onofri inequality in the $n$-dimensional unit ball $B^n:=B^n(0,1)$,

$$ \mathop {\inf }\limits_{u \in W_0^{1,n}(B^n)}\frac{1}{ n C_n}\int_{B^n} | \nabla u|^n dx- \ln \int_{B^n} {e^u} dx\geq C(n), $$

which extends the result of Caglioti-Lions-Marchioro-Pulvirenti in \cite{Caglioti} to the case of $n$-dimensional ball. Here $C_n=(\frac{n^2}{n-1})^{n-1} \omega_{n-1}$ and $\omega_{n-1}$ is the surface measure of $B^n$. For the Moser-Onofri inequality in a general bounded domain of $\mathbb{R}^n$, we apply the technique of $n$-harmonic transplantation to give the optimal concentration level of the Moser-Onofri inequality and obtain the criterion for the existence and non-existence of extremals for the Moser-Onofri inequality.

15. A parallel solver for random input problems via Karhunen-Loève expansion and diagonalized coarse grid correction

基于Karhunen-Loève展开与对角化粗网格校正的随机输入问题并行求解器

作者: Dou Dai, Qiuqi Li, Huailing Song

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

发布日期: 2025-10-30

arXiv ID: 2510.26180v1

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

摘要: This paper is dedicated to enhancing the computational efficiency of traditional parallel-in-time methods for solving stochastic initial-value problems. The standard parareal algorithm often suffers from slow convergence when applied to problems with stochastic inputs, primarily due to the poor quality of the initial guess. To address this issue, we propose a hybrid parallel algorithm, termed KLE-CGC, which integrates the Karhunen-Lo`{e}ve (KL) expansion with the coarse grid correction (CGC). The method first employs the KL expansion to achieve a low-dimensional parameterization of high-dimensional stochastic parameter fields. Subsequently, a generalized Polynomial Chaos (gPC) spectral surrogate model is constructed to enable rapid prediction of the solution field. Utilizing this prediction as the initial value significantly improves the initial accuracy for the parareal iterations. A rigorous convergence analysis is provided, establishing that the proposed framework retains the same theoretical convergence rate as the standard parareal algorithm. Numerical experiments demonstrate that KLE-CGC maintains the same convergence order as the original algorithm while substantially reducing the number of iterations and improving parallel scalability.

16. A one-dimensional Stefan problem for the heat equation with a nonlinear boundary condition

具有非线性边界条件的热传导方程的一维Stefan问题

作者: Kensho Araya, Kazuhiro Ishige

类别: math.AP

发布日期: 2025-10-30

arXiv ID: 2510.26088v1

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

摘要: We study the one-dimensional one-phase Stefan problem for the heat equation with a nonlinear boundary condition. We show that all solutions fall into one of three distinct types: global-in-time solutions with exponential decay, global-in-time solutions with non-exponential decay, and finite-time blow-up solutions. The classification depends on the size of the initial function. Furthermore, we describe the behavior of solutions at the blow-up time.

17. Bochner-Riesz means on a conical singular manifold

锥形奇性流形上的Bochner-Riesz平均

作者: Qiuye Jia, Junyong Zhang, Jiqiang Zheng

类别: math.AP, math.SP

发布日期: 2025-10-30

arXiv ID: 2510.26059v1

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

摘要: We prove a sharp $L^p$-boundedness criterion for Bochner-Riesz multipliers on flat cones $X = (0,\infty) \times \mathbb{S}_\sigma^1$. The operator $S_\lambda^\delta(\Delta_X)$ is bounded on $L^p(X)$ for $1 \leq p \leq \infty$, $p \neq 2$, if and only if $\delta > \delta_c(p,2) = \max\left{ 0, 2\left| 1/2

  • 1/p \right| - 1/2 \right}$. This result is also applicable to the infinite sector domain with Dirichlet or Neumann boundary, resolving the critical exponent problem in this wedge setting.

18. Non-Markovian heat flows on directed hypergraphs

有向超图上的非马尔可夫热流

作者: Delio Mugnolo

类别: math.DS, math.AP, math.CO, 47D06, 05C50, 05C65, 39A12

发布日期: 2025-10-20

arXiv ID: 2510.17497v2

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

摘要: We introduce a semigroup framework for Laplacians on directed hypergraphs, extending the classical heat flow models on graphs and establishing hypergraphs as prototypical models for non-Markovian diffusion. We apply spectral surgery methods to derive eigenvalue bounds, thus describing large-time behaviour of the heat flow. Unlike on standard graphs, heat flows on directed hypergraphs may lose positivity and/or $\infty$-contractivity, yet can recover them eventually or asymptotically under specific combinatorial configurations: examples based on duals of oriented graph and realisations of the Fano plane illustrate these phenomena. Our approach combines combinatorial, order-theoretic and linear-algebraic methods.

19. A game-theoretic approach to the parabolic normalized p-Laplacian obstacle problem

抛物型归一化p-拉普拉斯障碍问题的博弈论方法

作者: Hamid El Bahja

类别: math.PR, math.AP, 5B65, 47G20, 31B05

发布日期: 2025-10-29

arXiv ID: 2510.25999v1

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

摘要: This paper establishes a probabilistic representation for the solution of the parabolic obstacle problem associated with the normalized $p$-Laplacian. We introduce a zero-sum stochastic tug-of-war game with noise in a space-time cylinder, where one player has the option to stop the game at any time to collect a payoff given by an obstacle function. We prove that the value functions of this game exist, satisfy a dynamic programming principle, and converge uniformly to the unique viscosity solution of the continuous obstacle problem as the step size $\varepsilon$ tends to zero.

20. Gradient Flow Sampler-based Distributionally Robust Optimization

基于梯度流采样器的分布鲁棒优化

作者: Zusen Xu, Jia-Jie Zhu

类别: math.OC, math.AP, stat.ML

发布日期: 2025-10-29

arXiv ID: 2510.25956v1

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

摘要: We propose a mathematically principled PDE gradient flow framework for distributionally robust optimization (DRO). Exploiting the recent advances in the intersection of Markov Chain Monte Carlo sampling and gradient flow theory, we show that our theoretical framework can be implemented as practical algorithms for sampling from worst-case distributions and, consequently, DRO. While numerous previous works have proposed various reformulation techniques and iterative algorithms, we contribute a sound gradient flow view of the distributional optimization that can be used to construct new algorithms. As an example of applications, we solve a class of Wasserstein and Sinkhorn DRO problems using the recently-discovered Wasserstein Fisher-Rao and Stein variational gradient flows. Notably, we also show some simple reductions of our framework recover exactly previously proposed popular DRO methods, and provide new insights into their theoretical limit and optimization dynamics. Numerical studies based on stochastic gradient descent provide empirical backing for our theoretical findings.

21. Solutions to Second-Order Nonlocal Evolution Equations Governed by Non-Autonomous Forms

由非自治形式支配的二阶非局部发展方程的解

作者: Sajid Ullah, Vittorio Colao

类别: math.AP

发布日期: 2025-10-29

arXiv ID: 2510.25881v1

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

摘要: Our main contributions include proving sufficient conditions for the existence of solution to a second order problem with nonzero nonlocal initial conditions, and providing a comprehensive analysis using fundamental solutions and fixed-point techniques. The theoretical results are illustrated through applications to partial differential equations, including vibrating viscoelastic membranes with time-dependent material properties and nonlocal memory effects.