生成时间: 2025年12月19日 03:59:49

重点关注论文数量: 5

了解领域论文数量: 2

不相关论文数量: 24

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

1. Uniform vanishing damping limit for the 2D inviscid Oldroyd-B model with fractional stress tensor diffusion

分数阶应力张量扩散下二维无粘Oldroyd-B模型的一致衰减阻尼极限

作者: Chen Liang, Zhaonan Luo, Zhaoyang Yin

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16436v1

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

详细分析

中文标题: 分数阶应力张量扩散下二维无粘Oldroyd-B模型的一致衰减阻尼极限

研究对象和背景

本文研究的是带有分数阶应力张量扩散的二维无粘Oldroyd-B模型，其控制方程如下：

$$ \begin{cases} \partial_t u + u \cdot \nabla u + \nabla P = \mathrm{div}\, \tau, & \mathrm{div}\, u = 0, \\ \partial_t \tau + u \cdot \nabla \tau + a\tau + Q(\nabla u, \tau) + (-\Delta)^\beta \tau = D(u), \\ u|_{t=0} = u_0, \quad \tau|_{t=0} = \tau_0. \end{cases} \tag{1.2} $$

其中：

  • $ u(t,x) $ 是速度场，
  • $ \tau(t,x) $ 是对称的约束应力张量，
  • $ P $ 是压力，
  • $ a \in [0,1] $ 是阻尼系数（damping coefficient），
  • $ \beta \in [\frac{1}{2},1) $ 是分数阶扩散的指数，$ (-\Delta)^\beta $ 表示分数拉普拉斯算子，定义为傅里叶乘子：
    $$ \mathcal{F}\left[(-\Delta)^\beta \tau\right](\xi) = |\xi|^{2\beta} \hat{\tau}(\xi), $$
  • $ D(u) = \frac{\nabla u + (\nabla u)^T}{2} $ 是变形张量，
  • $ \Omega(u) = \frac{\nabla u - (\nabla u)^T}{2} $ 是涡量张量，
  • $ Q(\nabla u, \tau) = \tau \Omega(u) - \Omega(u) \tau - b(D(u)\tau + \tau D(u)) $，其中 $ b \in [-1,1] $。当 $ b=0 $ 时称为共旋（co-rotational）情形。

该系统描述了不可压缩、无粘性（$\nu = 0$）、具有分数阶扩散机制的黏弹性流体动力学行为。与经典的Navier-Stokes或Oldroyd-B模型不同，此处没有速度的粘性项（即无 $\Delta u$），但引入了对应力张量 $\tau$ 的分数阶耗散 $(-\Delta)^\beta \tau$，这在物理上可解释为分子链内部摩擦导致的非局部能量耗散。

数学挑战与研究动机

传统Oldroyd-B模型中，若缺乏速度粘性（$\nu = 0$）和足够强的应力耗散（如 $\beta < 1/2$），通常难以建立全局适定性。然而，本文发现：即使在无粘情况下，只要 $ \beta \in [\frac{1}{2},1) $，结合适当的阻尼项 $a\tau$，仍可在临界Sobolev空间中获得全局正则解。

特别地，作者关注的是阻尼趋于零时解的行为一致性——即所谓的“一致衰减阻尼极限”问题。具体而言，考虑当 $ a \to 0^+ $ 时，带阻尼解 $ (u_a, \tau_a) $ 是否一致收敛到无阻尼解 $ (u_0, \tau_0) $，并给出收敛速率。

这一问题不仅涉及渐近分析，还要求对所有 $ a \in [0,1] $ 建立统一的时间衰减速率估计，从而保证极限过程的一致性。

空间区域与函数空间设定

研究是在全空间 $ \mathbb{R}^2 $ 上进行的，初始数据位于Sobolev空间 $ H^s $ 中，并附加Besov型低频条件 $ \dot{B}^{-1}_{2,\infty} $，以控制低频部分的能量分布，这对时间衰减估计至关重要。

尤其值得注意的是，对于临界正则性的情形，作者选取 $ s = 2\beta $ 或 $ s > 1 $（当 $ \beta = 1 $ 时）。这种选择使得能量估计恰好处于尺度不变的空间中，体现了“临界性”。

此外，在证明时间衰减时，使用了改进的Fourier分裂方法（improved Fourier splitting method），该方法通过将频率空间分为低频与高频区域，分别处理能量演化，是获得最优衰减速率的关键工具。

另一个重要技术点是构造了包含交叉项 $ \langle \tau, -\Lambda^{-2+2\beta} \nabla u \rangle $ 的修正能量泛函，用于生成速度场的额外耗散项，这是克服无粘条件下缺乏速度耗散的核心技巧。

综上所述，本工作的数学框架具有以下特征：

  • 无粘性（$\nu = 0$）；
  • 分数阶应力扩散（$\mu = 0$, 但有 $(-\Delta)^\beta \tau$）；
  • 弱阻尼（$a \in [0,1]$）；
  • 二维空间（$d=2$）；
  • 临界正则性假设（$s = 2\beta$）；
  • 一致估计目标：关于参数 $a$ 的均匀性。

这些问题组合在一起构成了一个高度非线性的、非局部的、且需保持参数一致性的偏微分方程系统分析难题。

主要定理或主要结果

定理 1.1（临界全局正则性）

设 $ d = 2 $, $ a \in [0,1] $, $ \frac{1}{2} \leq \beta < 1 $, $ s \geq 2\beta $ 或 $ \beta = 1 $, $ s > 1 $。令初值 $ (u_0, \tau_0) \in H^s $ 且 $ \tau_0 $ 对称。存在小常数 $ \delta > 0 $，使得当

$$ \|(u_0, \tau_0)\|_{H^s} \leq \delta $$

时，系统 (1.2) 存在唯一全局强解 $ (u_a, \tau_a) \in C([0,\infty); H^s) $。进一步地，对任意 $ t > 0 $，有：

$$ \frac{d}{dt} \left( \|(u_a, \tau_a)\|_{H^s}^2 + 2k \langle -\nabla u_a, \tau_a \rangle_{H^{s-\beta}} \right) + \frac{k}{2} \|\nabla u_a\|_{H^{s-\beta}}^2 + \|\Lambda^\beta \tau_a\|_{H^s}^2 \leq 0, $$

其中 $ k > 0 $ 是充分小的常数。

注记：此结果揭示了分数阶扩散 $ (-\Delta)^\beta $ 在 $ \beta \in [\frac{1}{2},1) $ 范围内足以补偿无粘性带来的困难，实现全局适定性；但当 $ \beta = 1 $ 时，该方法失效，需更高正则性。

定理 1.4（无阻尼情形下的最优时间衰减速率）

设 $ d = 2 $, $ a = 0 $，且初值满足定理1.1的条件。若额外假设 $ (u_0, \tau_0) \in \dot{B}^{-1}_{2,\infty} $，则存在 $ C > 0 $，使得对任意 $ t > 0 $ 和 $ s_1 \in [0,s] $，有：

$$ \|\Lambda^{s_1}(u_0, \tau_0)(t)\|_{L^2} \leq C(1+t)^{-\frac{s_1}{2} + \frac{1}{2\beta}}. $$

若进一步假设 $ (u_0, \tau_0) \in H^{s+\beta} $ 且 $ \left| \int_{\mathbb{R}^2} (u_0, \tau_0) dx \right| > 0 $，则存在 $ C_\beta > 0 $，使得：

$$ \|\Lambda^{s_1}(u_0, \tau_0)(t)\|_{L^2} \geq C_\beta (1+t)^{-\frac{s_1}{2} + \frac{1}{2\beta}}, $$

即上下界匹配，衰减速率是最优的。

定理 1.5（一致时间衰减与可积性）

设 $ d = 2 $, $ a \in [0,1] $, $ \frac{1}{2} \leq \beta < 1 $, $ s \geq 1 + \beta $，初值满足定理1.1条件且 $ (u_0, \tau_0) \in \dot{B}^{-1}_{2,\infty} $。则存在常数 $ C > 0 $，使得：

$$ \|(u_a, \tau_a)\|_{L^2} + (1+t)^{1/2} \|\nabla(u_a, \tau_a)\|_{L^2} \leq C(1+t)^{-1/2}. $$

若 $ s > 1 + \beta $，还有：

$$ \int_0^\infty \|\nabla(u_a, \tau_a)\|_{L^\infty} dt \leq C. $$

此估计对所有 $ a \in [0,1] $ 成立，具有一致性，是后续极限分析的基础。

定理 1.6（一致阻尼消失极限）

在相同条件下，若 $ \frac{1}{2} < \beta < 1 $，则对任意 $ \alpha \in [0,1] $，有：

$$ \|\nabla^\alpha(u_a - u_0, \tau_a - \tau_0)\|_{L^\infty([0,\infty); L^2)} \leq C a^{\frac{\alpha}{\beta(\alpha + 3) - 1}}. $$

当 $ \beta = \frac{1}{2} $ 时，对任意 $ \varepsilon \in (0,1) $，有：

$$ \|\nabla^\alpha(u_a - u_0, \tau_a - \tau_0)\|_{L^\infty([0,\infty); L^2)} \leq C a^{1 - \varepsilon}. $$

这表明当阻尼系数趋于零时，解一致收敛于无阻尼解，且给出了明确的收敛速率。

定理 1.7（迹函数 tr$\tau_a$ 的锐利一致衰减）

令 $ \tilde{\tau}_a = \mathrm{tr}\,\tau_a $，则有：

$$ \|\tilde{\tau}_a(t)\|_{L^2} \leq C(1+t)^{-\frac{1}{2\beta}}. $$

若 $ \frac{1}{2} < \beta < 1 $ 且 $ s > 1 + \beta $，则：

$$ \|\tilde{\tau}_a - \tilde{\tau}_0\|_{L^\infty([0,\infty); L^2)} \leq C a^{\frac{1}{2\beta}}. $$

并在命题5.10中说明该速率是最优的（sharp），即存在初值使下界成立。

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

a) 核心分析工具

  1. 能量估计方法（Energy Estimates）

    全局适定性的基础是构造合适的能量泛函并建立先验估计。由于系统无速度粘性，标准动能耗散缺失，因此必须借助交叉耦合项 $ \langle -\nabla u, \tau \rangle $ 构造修正能量：

    $$ \mathcal{E}(t) = \|(u, \tau)\|_{H^s}^2 + 2k \langle -\nabla u, \tau \rangle_{H^{s-\beta}}, $$

    其中 $ k > 0 $ 小。通过对时间求导并利用分部积分、交换子估计等手段，可从 $ \|\Lambda^\beta \tau\|^2 $ 项中提取出 $ \|\nabla u\|^2 $ 的耗散，形成闭合的能量不等式。

  2. 改进的Fourier分裂方法（Improved Fourier Splitting Method）

    用于证明时间衰减速率。基本思想是将频率空间划分为低频集 $ B(t) = \{ |\xi| \leq g(t) \} $ 与高频补集，利用低频部分由初值主导、高频受耗散控制的特性，递推获得逐层衰减。

    特别地，作者引入了一个新的加权能量估计来替代复杂的谱分析，简化了下界的证明。

  3. Littlewood-Paley 分解与 Besov 空间理论

    文中多次使用齐次Besov空间 $ \dot{B}^s_{p,q} $ 来刻画初值的低频性质。例如，$ (u_0, \tau_0) \in \dot{B}^{-1}_{2,\infty} $ 意味着其低频部分满足某种均匀有界性，这对启动Fourier分裂至关重要。

    同时，利用该分解处理非线性项中的乘积与交换子。

  4. 交换子与乘积估计（Commutator and Product Estimates）

    非线性项如 $ [\Lambda^s, u \cdot \nabla]u $ 的控制依赖于经典引理：

    $$ \|[\Lambda^s, f]g\|_{L^p} \lesssim \|\Lambda^s f\|_{L^{p_1}} \|g\|_{L^{p_2}} + \|\nabla f\|_{L^{p_3}} \|\Lambda^{s-1} g\|_{L^{p_4}}, $$

    以及乘积估计：

    $$ \|\Lambda^s(fg)\|_{L^p} \lesssim \|f\|_{L^{p_1}} \|\Lambda^s g\|_{L^{p_2}} + \|g\|_{L^{p_3}} \|\Lambda^s f\|_{L^{p_4}}, $$

    其中 $ \frac{1}{p} = \frac{1}{p_i} + \frac{1}{p_j} $。这些工具允许将高阶导数作用于非线性项并有效控制。

  5. Gagliardo-Nirenberg 不等式

    在二维情形下，用于插值不同范数：

    $$ \|\Lambda^s f\|_{L^p} \lesssim \|\Lambda^{s_1} f\|_{L^2}^{1-\theta} \|\Lambda^{s_2} f\|_{L^2}^{\theta}, $$

    其中 $ s + \frac{2}{p} - 1 = (1-\theta)s_1 + \theta s_2 $，$ 0 < \theta < 1 $。这对连接 $ L^2 $、$ L^\infty $ 和 Sobolev 范数非常关键。

  6. 谱分析方法（Spectral Analysis）

    用于改进 $ \mathrm{tr}\,\tau_a $ 的衰减速率。注意到 $ \tilde{\tau}_a = \mathrm{tr}\,\tau_a $ 满足近似分数热方程：

    $$ \partial_t \tilde{\tau}_a + (-\Delta)^\beta \tilde{\tau}_a + a \tilde{\tau}_a = Q_1(u_a, \nabla \tau_a) + Q_2(\nabla u_a, \tau_a), $$

    其中右边为二次项。利用谱投影将解分解为低频（扩散主导）与高频（阻尼主导），可精细分析其长期行为。

b) 关键技术难点与突破

  1. 无粘性下的速度耗散生成

    由于 $ \nu = 0 $，无法直接从 $ \Delta u $ 获取速度耗散。作者通过构造交叉项 $ \langle \tau, -\Lambda^{-2+2\beta} \nabla u \rangle $，将其加入能量泛函，从而在时间导数中诱导出 $ \|\Lambda^\beta u\|^2 $ 类型的耗散项。这是整个能量结构得以闭合的关键。

  2. 阻尼参数 $ a $ 的一致性处理

    当 $ a > 0 $ 时，虽然有额外阻尼项 $ a\tau $，但在估计交叉项时会出现形如 $ a \langle \tau_a, \Lambda^{-2+2\beta} \nabla u_a \rangle $ 的项，而该项不能被现有耗散控制（因 $ \|\Lambda^{-2+2\beta} \nabla u_a\|_{L^2} $ 不受控）。为此，作者放弃追求“最优”衰减速率，转而采用迭代方式建立一致但稍弱的衰减估计，最终仍能满足积分可积性需求。

  3. 临界正则性边界效应

    当 $ \beta = 1 $ 时，原有能量估计失效，因为需要 $ \|u\|_{L^\infty} $ 控制，而这在 $ H^1 $ 中不封闭。这说明 $ \beta = 1 $ 是一个临界阈值，也凸显了分数阶指数的选择敏感性。

  4. trace项的锐利估计

    利用 $ \mathrm{div}\,D u = 0 $，可知 $ \mathrm{tr}(D u) = \mathrm{div}\,u = 0 $，从而 $ \mathrm{tr}(\partial_t \tau + u\cdot\nabla\tau + (-\Delta)^\beta \tau + a\tau) = 0 $ 加上线性化后的结构，使得 $ \tilde{\tau}_a $ 接近分数热方程解，便于应用谱方法获得更优衰减。

与之前工作的比较

放宽的假设条件

  1. 降低初始数据正则性要求

    相较于早期工作如 Elgindi–Rousset [10] 要求 $ H^s $ ($ s > 2 $) 或 Besov空间中的严格小性条件，本文在临界Sobolev空间 $ H^{2\beta} $ 中建立全局适定性，显著降低了正则性门槛。特别是当 $ \beta = \frac{1}{2} $ 时，仅需 $ H^1 $ 初值即可。

  2. 推广维度至二维临界情形

    多数已有结果集中在三维或高正则性情形。本文专注于二维，并系统研究 $ \beta \in [\frac{1}{2},1) $ 区间内的过渡行为，填补了该参数范围的研究空白。

  3. 允许 $ a = 0 $ 的极限情形

    多数含阻尼模型要求 $ a > 0 $ 才能获得全局解。本文首次在统一框架下处理 $ a \in [0,1] $，包括 $ a = 0 $ 的无阻尼极限，为渐近分析提供基础。

结果强度提升

  1. 建立了一致性估计而非逐点收敛

    以往文献多证明对每个固定 $ t $，$ (u_a, \tau_a)(t) \to (u_0, \tau_0)(t) $，但未必一致。本文首次给出在时间区间 $ [0,\infty) $ 上一致成立的收敛速率，更强且更具应用价值。

  2. 获得 sharp 衰减速率

    不仅证明了 $ \|\tilde{\tau}_a - \tilde{\tau}_0\|_{L^\infty_t L^2_x} \lesssim a^{1/(2\beta)} $，还构造反例说明该速率不可改进，达到“最优”，这在同类工作中罕见。

  3. 发展了无需复杂谱分析的新方法

    以往下界估计常依赖精细谱投影。本文提出一种新型加权能量估计法，避免了繁复计算，技术更简洁，可能适用于其他系统。

与其他模型对比

  • 与 Wu–Zhao [29] 相比：他们研究的是 Besov 框架下的全局适定性，本文推进到 Sobolev 临界空间，并加入阻尼极限分析。
  • 与 Chen–Luo 等 [6] 相比：他们研究整数阶扩散（$ \beta = 1 $）情形，得到 $ a^{d/4} $ 的消失速率。本文将其推广至分数阶 $ \beta \in [\frac{1}{2},1) $，并获得依赖于 $ \beta $ 的新速率 $ a^{\alpha / (\beta(\alpha+3)-1)} $，更具一般性。

总之，本文在模型广度、正则性深度、估计精度与技术新颖性四个方面均实现了显著超越。

2. Low-Mach-number limit for multiphase flows

多相流动的低马赫数极限

作者: Cassandre Lebot

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16286v1

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

详细分析

中文标题: 多相流动的低马赫数极限

研究对象和背景

本文旨在系统研究多相流动（multiphase flows）在低马赫数极限（low-Mach-number limit）下的渐近行为。马赫数 $M$ 是一个无量纲参数，定义为流体速度与声速之比：$M = |u|/c$。当 $M \to 0$ 时，流体运动远慢于声波传播，此时可期望压缩性效应减弱，系统趋于不可压流动。这一极限过程在理论分析与数值模拟中均有重要意义：一方面有助于理解不可压方程的来源；另一方面，低马赫数模型能显著简化数值算法并允许更大的时间步长。

论文首先回顾了单相流体中的低马赫数极限理论，随后重点聚焦于两相流系统，考虑了多种建模方式，包括具有代数压力闭合（algebraic closure）或偏微分方程压力闭合（PDE closure）的情形，并区分了单速度场与双速度场模型、等熵与非等熵情形。

单相流模型（等熵情况）

在等熵假设下，压缩性Navier-Stokes方程组描述如下：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) + \nabla p(\rho) - \mu \Delta u - (\lambda + \mu)\nabla \mathrm{div}\, u = 0, \end{cases} $$

其中 $\rho \geq 0$ 为密度，$u$ 为速度场，压力满足状态方程 $p(\rho) = \rho^\gamma$，$\gamma > 1$，粘性系数 $\mu > 0$，$\lambda + 2\mu > 0$。该系统在区域 $\Omega \subset \mathbb{R}^d$（$d \geq 2$）上定义，通常考虑周期边界条件或全空间情形。

引入小参数 $\varepsilon > 0$ 表示马赫数，通过尺度变换：

$$ u(t,x) = \varepsilon u_\varepsilon(\varepsilon^2 t, \varepsilon x), \quad \rho(t,x) = \rho_\varepsilon(\varepsilon^2 t, \varepsilon x), $$

得到无量纲化系统：

$$ \begin{cases} \partial_t \rho_\varepsilon + \mathrm{div}(\rho_\varepsilon u_\varepsilon) = 0, \\ \partial_t (\rho_\varepsilon u_\varepsilon) + \mathrm{div}(\rho_\varepsilon u_\varepsilon \otimes u_\varepsilon) + \frac{1}{\varepsilon^2} \nabla p(\rho_\varepsilon) - \mu_\varepsilon \Delta u_\varepsilon - (\lambda_\varepsilon + \mu_\varepsilon)\nabla \mathrm{div}\, u_\varepsilon = 0. \end{cases} $$

形式上令 $\varepsilon \to 0$，由动量方程中 $1/\varepsilon^2$ 项主导得 $\nabla p(\rho_\varepsilon) \to 0$，从而极限密度 $\rho_0$ 满足 $p(\rho_0) = C(t)$，即常值压力。结合质量守恒方程可推出 $\mathrm{div}\, u_0 = 0$，最终极限系统为不可压Navier-Stokes方程：

$$ \begin{cases} \mathrm{div}\, u_0 = 0, \\ \rho_0 (\partial_t u_0 + (u_0 \cdot \nabla) u_0) - \mu \Delta u_0 + \nabla \pi = 0, \end{cases} $$

其中 $\pi$ 为压力极限项，来源于 $(\rho_\varepsilon^\gamma - \bar{\rho}^\gamma)/\varepsilon^2$ 的收敛。

若初始数据满足 $\mathrm{div}\, u_\varepsilon = O(\varepsilon)$ 且 $\rho_\varepsilon \to \bar{\rho}$（常数），称为“适定初值”（well-prepared data）；否则为“不适定初值”，可能激发高速传播的声波。本文主要关注适定初值情形。

非等熵单相流

对于非等熵Euler系统：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) + \nabla p(\rho, S) = 0, \\ \partial_t S + u \cdot \nabla S = 0, \end{cases} $$

其中 $S$ 为熵，压力 $p = p(\rho,S)$，例如取 $p(\rho,S) = \rho^\gamma e^S$。类似地进行尺度变换并引入变量 $q_\varepsilon$ 使得 $p = \bar{p} e^{\varepsilon q_\varepsilon}$，可得：

$$ \begin{cases} a_\varepsilon (\partial_t q_\varepsilon + u_\varepsilon \cdot \nabla q_\varepsilon) + \frac{1}{\varepsilon} \mathrm{div}\, u_\varepsilon = 0, \\ r_\varepsilon (\partial_t u_\varepsilon + u_\varepsilon \cdot \nabla u_\varepsilon) + \frac{1}{\varepsilon} \nabla q_\varepsilon = 0, \\ \partial_t S_\varepsilon + u_\varepsilon \cdot \nabla S_\varepsilon = 0, \end{cases} $$

其中 $a_\varepsilon, r_\varepsilon$ 依赖于 $\rho_\varepsilon, S_\varepsilon$。形式极限 $\varepsilon \to 0$ 给出：

$$ \mathrm{div}\, u_0 = 0, \quad \nabla q_0 = 0 \Rightarrow q_0 = \text{const}, \quad p_0 = \text{const}, $$

进而得到不可压Euler系统耦合输运方程：

$$ \begin{cases} r_0 (\partial_t u_0 + u_0 \cdot \nabla u_0) + \nabla \Pi = 0, \\ \partial_t S_0 + u_0 \cdot \nabla S_0 = 0, \\ \mathrm{div}\, u_0 = 0. \end{cases} $$

两相流模型（单速度场）

考虑两相 $+$ 和 $-$，其密度 $\rho_i$、体积分数 $\alpha_i$、压力 $p_i(\rho_i)$，共用速度场 $u$。系统分为两类：

(1) 压力代数闭合系统

$$ \begin{cases} \alpha_+ + \alpha_- = 1, \\ \partial_t (\alpha_+ \rho_+) + \mathrm{div}(\alpha_+ \rho_+ u) = 0, \\ \partial_t (\alpha_- \rho_-) + \mathrm{div}(\alpha_- \rho_- u) = 0, \\ \partial_t ((\alpha_+ \rho_+ + \alpha_- \rho_-) u) + \mathrm{div}((\alpha_+ \rho_+ + \alpha_- \rho_-) u \otimes u) \\ \quad + \nabla(\alpha_+ p_+(\rho_+) + \alpha_- p_-(\rho_-)) - 2\,\mathrm{div}\left( (\alpha_+ \mu_+ + \alpha_- \mu_-) D(u) \right) \\ \quad - \nabla\left( (\alpha_+ \lambda_+ + \alpha_- \lambda_-) \mathrm{div}\, u \right) = (\alpha_+ \rho_+ + \alpha_- \rho_-) f, \\ p_+(\rho_+) = p_-(\rho_-), \end{cases} $$

其中 $D(u) = \frac{\nabla u + (\nabla u)^T}{2}$ 为应变率张量，最后的压力平衡条件是典型的界面平衡假设。

(2) 压力PDE闭合系统

与上述不同的是，取消压力相等假设，代之以演化方程：

$$ \partial_t \alpha_+ + u \cdot \nabla \alpha_+ = \frac{\alpha_+ \alpha_-}{\mu} (p_+(\rho_+) - p_-(\rho_-)), $$

其余方程相同。此模型允许压力不等，通过该项驱动相间转移。

对应的能量估计分别为：

  • 代数闭合：

    $$ \frac{d}{dt} \int_\Omega \left[ \frac{1}{2} (\alpha_+ \rho_+ + \alpha_- \rho_-)|u|^2 + \frac{\alpha_+ \rho_+^{\gamma_+}}{\gamma_+ - 1} + \frac{\alpha_- \rho_-^{\gamma_-}}{\gamma_- - 1} \right] dx + \text{viscous terms} \leq \int (\cdots) f \cdot u\,dx $$

  • PDE闭合：额外多出一项耗散项： $$

  • \int_0^t \int_\Omega \frac{\alpha_+ \alpha_-}{\mu} (p_+(\rho_+) - p_-(\rho_-))^2 dx dt, $$ 体现相间压力差引起的能量耗散。

区域 $\Omega$ 可为有界域（具周期边界）或全空间 $\mathbb{R}^d$，函数空间一般取 Sobolev 空间如 $H^s$ 或能量类空间。由于涉及弱解与极限过程，常需使用紧性方法，在 $L^p$ 弱拓扑或测度值解框架下工作。

主要定理或主要结果

本文未明确陈述严格定理（Theorem 或 Proposition 形式），而是以形式推导（formal derivation）为主，总结已有结果并提出开放问题。因此，其“主要结果”体现在对各类两相系统在低马赫数极限下的形式渐近分析及其物理与数学结构的揭示。

核心结论如下：

  1. 对于单相流：在适定初值条件下，压缩性Navier-Stokes/Euler方程在 $\varepsilon \to 0$ 极限下形式收敛至不可压Navier-Stokes/Euler系统，压力趋于常数，速度场满足不可压条件。

  2. 对于两相单速度系统（代数闭合）：在低马赫数极限下，质量守恒方程保持不变，动量方程中压力梯度项导致极限压力为常数，速度场趋于不可压，即 $\mathrm{div}\, u_0 = 0$，极限系统为不可压两相流模型，压力均匀。

  3. 对于两相单速度系统（PDE闭合）：极限过程中，除了获得 $\mathrm{div}\, u_0 = 0$ 外，压力差驱动项 $\frac{\alpha_+ \alpha_-}{\mu}(p_+ - p_-)$ 在 $\varepsilon \to 0$ 下仍保留，但受控于新的平衡机制。形式极限中，该项可能趋于零或产生约束，取决于尺度。

  4. 对于双速度两相系统：当两相具有不同速度时，低马赫数极限更为复杂。论文指出，若两相速度均以相同尺度缩放，则极限中可能出现两个不可压速度场，或通过松弛项强制速度一致。是否收敛至共同不可压流依赖于相对速度的初始尺度。

  5. 适定初值的重要性：文中强调，若初始数据不满足 $\mathrm{div}\, u_\varepsilon = O(\varepsilon)$ 或密度接近常值，则会产生高频声波，需用滤波方法或几何光学处理，超出本文范围。

综上，虽无严格定理，但作者系统梳理了从压缩性到不可压两相系统的形式极限路径，为后续严格数学分析提供了蓝图。

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

a) 核心分析学工具

  1. 形式渐近分析（Formal Asymptotics）
    本文最核心的方法是形式极限分析，即不证明解的存在性或收敛性，而直接对控制方程进行渐近展开，观察 $\varepsilon \to 0$ 时主导项的行为。这是研究奇异极限的经典起点，尤其适用于多尺度PDE系统。

  2. 能量估计与能量泛函
    能量不等式是贯穿全文的核心工具。例如，单相Navier-Stokes的能量不等式：

    $$ \frac{d}{dt} \int_\Omega \left( \frac{1}{2} \rho_\varepsilon |u_\varepsilon|^2 + \frac{1}{\varepsilon^2(\gamma-1)} \rho_\varepsilon^\gamma \right) dx + \int_0^t \int_\Omega \left( \mu_\varepsilon |\nabla u_\varepsilon|^2 + \lambda_\varepsilon (\mathrm{div}\, u_\varepsilon)^2 \right) dx dt \leq C $$

    提供了关于 $u_\varepsilon$ 和 $\rho_\varepsilon$ 的先验估计。在两相系统中，作者详细计算了能量耗散结构（见公式(3.8)-(3.9)），揭示体积分数演化如何影响总能量变化。这种显式能量恒等式推导是判断模型物理合理性的关键。

  3. 紧性原理（Compactness Arguments）
    尽管本文未给出证明，但引用了大量已有工作（如[Lions and Masmoudi, 1998]）使用弱紧性方法：利用能量估计获得解序列的有界性，再通过 Aubin-Lions 引理、Rellich-Kondrachov 定理等提取强收敛子列。特别地，在周期域中，Fourier 分析与均匀估计结合可用于处理声波振荡。

  4. 滤波方法（Filtering Techniques）
    对于不适定初值，声波以 $O(1/\varepsilon)$ 速度传播，形成快速振荡。[Schochet, 1986] 等人发展的滤波技术用于分离慢变（不可压）模式与快变（声波）模式。本文虽未展开，但提及此法为处理一般初值的关键。

  5. 尺度变换与无量纲化
    所有分析基于尺度变换：将原始变量按特征尺度重标，使方程中出现小参数 $\varepsilon$。这一步骤将物理问题转化为数学上的奇异扰动问题，是连接可压与不可压模型的桥梁。

  6. 变分结构与耗散机制识别
    在PDE闭合模型中，新增项 $\frac{\alpha_+\alpha_-}{\mu}(p_+ - p_-)^2$ 显式出现在能量不等式中，表明系统具有更强的内在耗散。识别此类结构有助于后续构造弱解或证明稳定性。

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

  1. 压力梯度项的奇异性
    动量方程中 $\frac{1}{\varepsilon^2} \nabla p(\rho_\varepsilon)$ 是最大阶项。处理技巧是将其视为约束力：在极限中迫使 $\nabla p(\rho_0) = 0$，从而限制密度行为。这要求对压力状态方程有良好控制（如凸性）。

  2. 非线性项的收敛问题
    如何证明 $\rho_\varepsilon u_\varepsilon \otimes u_\varepsilon \to \rho_0 u_0 \otimes u_0$？常用技巧是利用弱-强收敛乘积引理：若 $\rho_\varepsilon \rightharpoonup \rho_0$ 弱收敛，$u_\varepsilon \to u_0$ 强收敛，则乘积收敛。为此需通过能量估计获得 $u_\varepsilon$ 的更高正则性。

  3. 体积分数与密度耦合的复杂性
    在两相系统中，$\alpha_i$ 与 $\rho_i$ 相互耦合，导致质量方程非标准。例如在PDE闭合中，$\partial_t \alpha_+$ 含压力差，需精细估计该项在极限中的行为。作者通过积分恒等式（如(3.8)-(3.9)）逐项分析，确保能量结构封闭。

  4. 双速度系统的松弛机制
    当两相速度不同时，系统可能包含速度松弛项（如 $-\frac{1}{\tau}(u_+ - u_-)$）。在低马赫数极限中，若松弛时间 $\tau = O(\varepsilon)$，则极限中两速度被迫相等；若 $\tau$ 固定，则可能出现两个独立不可压场。此尺度匹配技巧是处理多速度系统的核心。

  5. 边界条件的影响
    在有界域中，若边界条件为 $u \cdot n = 0$，声波可能被吸收或反射。[Desjardins et al., 1999] 利用几何假设（如凸域）证明声波能量衰减，从而获得强收敛。这是克服边界振荡的关键技术。

与之前工作的比较

a) 突破与改进

  1. 首次系统综述两相流低马赫数极限
    此前仅有少数文献（如[Varsakelis and Papalexandris, 2011]、[Narbona-Reina et al., 2024]）零星研究此问题。本文是首个对该领域进行全面分类与形式推导的工作，涵盖代数/PDE闭合、单/双速度、等熵/非等熵等多种模型，填补了综述空白。

  2. 推广建模范式
    相较于传统Baer–Nunziato模型，本文讨论的平均化两相模型（如[Bresch et al., 2022]）更具物理基础，源于均匀化技术而非经验假设。这标志着从“现象建模”向“推导建模”的转变。

  3. 强调能量结构的一致性
    作者不仅写出方程，更详细推导了各类系统的能量不等式（如(3.2)、(3.11)），并指出其与单相系统的继承关系（见Remark 3.3）。这种对热力学一致性的关注，超越了仅关注形式极限的早期工作。

  4. 识别PDE闭合的额外耗散机制
    在压力PDE闭合模型中，发现新增的 $(p_+ - p_-)^2$ 耗散项（见(3.11)加粗项），这在代数闭合中不存在。这一观察提示PDE闭合模型可能具有更好的数学适定性，为后续存在性证明提供线索。

  5. 统一处理框架
    本文将单相、两相、等熵、非等熵置于同一分析框架下，使用相同的尺度变换与极限逻辑，体现了方法论的普适性。相比以往针对特定系统的孤立研究，更具系统性。

改进体现于假设放宽与结果强化

  • 对初值假设的潜在放宽：虽然本文限于适定初值，但通过引用[Ukai, 1986]、[Grenier, 1997]等工作，暗示未来可推广至一般初值，只需处理声波层。

  • 维度推广：现有严格结果多限于一维或周期域，本文讨论高维情形，推动向 $\mathbb{R}^d$（$d=2,3$）发展。

  • 结果强度：形式推导为后续强收敛或误差估计奠定基础。例如，若能量估计一致，则可用弱-强收敛结合获得几乎处处收敛。

总之，本文虽以综述与形式分析为主，但通过清晰分类、能量结构揭示与开放问题提出，为该领域的严格数学理论发展提供了重要方向。

3. Well-Posedness for Low Regularity Solutions to the g-SQG Equation with Regular Level Sets

关于具有正则水平集的g-SQG方程低正则性解的适定性研究

作者: Junekey Jeon, Andrej Zlatos

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16128v1

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

详细分析

中文标题: 关于具有正则水平集的g-SQG方程低正则性解的适定性研究

研究对象和背景

本文研究的是平面上的广义表面准地转方程（generalized Surface Quasi-Geostrophic equation, 简称 g-SQG），其数学形式为如下主动标量偏微分方程：

$$ \partial_t \theta + u(\theta) \cdot \nabla \theta = 0, $$

其中未知函数 $\theta = \theta(t,x): I \times \mathbb{R}^2 \to \mathbb{R}$ 是一个在时间区间 $I$ 上定义的空间标量场，而速度场 $u(\theta)$ 由下式给出：

$$ u(\theta) := -\nabla^\perp (-\Delta)^{-1+\alpha} \theta, \quad \text{其中 } \alpha \in \left(0, \frac{1}{2}\right), $$

这里 $\nabla^\perp := (-\partial_{x_2}, \partial_{x_1})$，且记 $(x_1, x_2)^\perp = (-x_2, x_1)$。通过积分核表示，该速度场也可写成卷积形式：

$$ u(\theta_t; x) = c_\alpha \int_{\mathbb{R}^2} \frac{(x - y)^\perp}{|x - y|^{2 + 2\alpha}} \theta_t(y) \, dy = \int_{\mathbb{R}^2} \nabla^\perp K(x - y) \theta_t(y) \, dy, $$

其中 $K(x) = -\frac{c_\alpha}{2\alpha |x|^{2\alpha}}$，$c_\alpha > 0$ 是仅依赖于 $\alpha$ 的常数，$\theta_t := \theta(t, \cdot)$。

当参数 $\alpha = 0$ 时，此方程退化为二维不可压缩欧拉方程中的涡度方程；当 $\alpha = \frac{1}{2}$ 时，则对应经典的SQG方程。因此，g-SQG可以被视为连接这两个经典模型的一族中间情形，在大气与海洋动力学中有重要物理意义。

该方程的一个关键特征是它是一个输运型PDE：解沿由自身诱导的速度场 $u(\theta)$ 的流线被传输。这意味着从Lagrangian视角出发理解其解结构尤为自然。具体而言，若初始数据为 $\theta_0$，则形式上解可表示为 $\theta_t = \theta_0 \circ \Phi_t^{-1}$，其中 $\Phi_t$ 是向量场 $u(\theta_t)$ 所生成的流映射，满足：

$$ \partial_t \Phi_t = u(\theta_t) \circ \Phi_t, \quad \Phi_0 = \mathrm{Id}. $$

由于算子 $(-\Delta)^{-1+\alpha}$ 在 $\alpha > 0$ 时具有奇异性，即使 $\theta$ 属于 $L^\infty \cap L^1$，对应的 $u(\theta)$ 通常仅为 $(1 - 2\alpha)$-Hölder 连续，不足以保证ODE流的存在唯一性。因此，传统的Cauchy-Lipschitz理论无法直接应用，必须对解空间施加额外正则性条件以确保适定性。

以往研究表明，当 $\theta_0 \in C^\beta(\mathbb{R}^2)$ 且 $\beta > 2\alpha$ 时，$u(\theta)$ 成为Lipschitz连续，从而局部适定性可在 $L^1 \cap C^\beta$ 中建立 [4]。然而，这类结果要求较高的整体函数正则性，排除了许多几何上有趣但不光滑的情形，例如具有分片常数结构的“patch”解。

事实上，“patch问题”——即初始数据为有限个互不相交、边界足够光滑的有界区域上的特征函数之和——也已被广泛研究。在这种情况下，只要所有初始patch边界 $\partial \Theta_n^0$ 是彼此分离的简单闭曲线并属于适当Sobolev类（如 $H^2$ 或更光滑），就可以证明局部适定性 [7,8]。这是因为此时动态完全由速度场在边界上的限制决定，而由于 $\theta$ 在切向近似恒定，$u(\theta)$ 沿边界表现出更高阶的正则性。

受此启发，本文提出一个核心问题：是否存在一类低正则性的连续解，其本身可能仅具 $C^{2\alpha}$ 正则性（即处于经典理论失效边缘），但由于其水平集边界保持高正则性（特别是 $H^2$ 类曲线，意味着曲率平方可积），仍能保证方程的局部甚至长时间适定性？

为此，作者引入了一种新的分析框架，基于“广义层蛋糕分解”（generalized layer cake representation）。对于任意 $\theta \in L^1 \cap L^\infty$，标准层蛋糕公式为：

$$ \theta(x) = \int_0^{\sup \theta} \mathbf{1}_{\{\theta > \lambda\}}(x) \, d\lambda - \int_{\inf \theta}^0 \mathbf{1}_{\{\theta < \lambda\}}(x) \, d\lambda. $$

而广义版本允许非嵌套的集合族：设 $(L, \mu)$ 为带符号测度的可测空间，$\Theta \subseteq \mathbb{R}^2 \times L$ 为乘积σ-代数中的集合，若对每个 $x \in \mathbb{R}^2$，

$$ \theta(x) = \int_L \mathbf{1}_{\Theta_\lambda}(x) \, d\mu(\lambda), \quad \text{其中 } \Theta_\lambda := \{x \in \mathbb{R}^2 : (x,\lambda) \in \Theta\}, $$

则称 $(\Theta, \mu)$ 是 $\theta$ 的一个广义层蛋糕表示。这种构造特别适合处理具有多个连通分支的水平集，尤其是当这些边界均为 $H^2$ 曲线时。

本文所考虑的函数空间本质上包含那些模为某个满足

$$ \int_0^1 \frac{\rho(s)}{s^{1+2\alpha}} \, ds < \infty $$

的连续递增函数 $\rho: [0,\infty) \to [0,\infty)$（满足 $\rho(0)=0$）的函数。典型例子包括 $\rho(s) = s^{2\alpha} |\log s|^{-p}$ 对 $p > 1$，这比 $C^{2\alpha}$ 稍强一点。值得注意的是，尽管此类函数整体上可能接近临界正则性，但若其水平集边界具备 $H^2$ 正则性，并满足某些密度控制条件（如 $R_\mu(\Theta), Q(\Theta) < \infty$，见后文定义），就能获得更强的动力学控制。

此外，文中采用的解概念是Lagrangian解，即存在连续流映射 $\Phi \in C(I; C(\mathbb{R}^2; \mathbb{R}^2))$，使得 $\theta_t = \Phi_{t*} \theta_0 := \theta_0 \circ \Phi_t^{-1}$，且 $\Phi_t$ 是保测同胚（因 $\nabla \cdot u = 0$）。这一设定天然兼容层蛋糕表示下的变换性质：$(\Phi_* \Theta)_\lambda = \Phi(\Theta_\lambda)$。

综上所述，本文的研究对象是一类介于经典高正则解与patch解之间的过渡情形：既不要求整体函数高度光滑，也不允许间断跳跃，而是强调几何结构性质——即所有水平集边界保持 $H^2$ 正则性。目标是在尽可能弱的整体函数正则性假设下，建立此类解的局部乃至部分全局适定性，并揭示奇点形成机制的本质必须涉及水平集几何正则性的丧失，而非简单的“碰撞”或“堆积”。

主要定理或主要结果

论文的核心成果是以下两个定理：

定理 1.3（一般低正则性解的局部适定性）

设初值 $\theta_0 \in L^1(\mathbb{R}^2) \cap L^\infty(\mathbb{R}^2)$ 存在一个广义层蛋糕表示 $(\Theta, \mu)$，使得量

$$ L_\mu(\Theta) := \sup_{x \in \mathbb{R}^2} \int_L \frac{d|\mu|(\lambda)}{d(x, \partial \Theta_\lambda)^{2\alpha}} < \infty. $$

则存在包含0的开时间区间 $I$ 和唯一的Lagrangian解 $\theta: I \to L^1 \cap L^\infty$ 到 g-SQG 方程，对应流映射 $\Phi \in C(I; C(\mathbb{R}^2; \mathbb{R}^2))$，满足对任意紧子区间 $J \subset I$，有 $\sup_{t \in J} L_\mu(\Phi_{t*} \Theta) < \infty$。令 $I$ 为极大这样的区间，则解唯一且独立于 $(\Theta, \mu)$ 的选择；进一步，对任意紧 $J \subset I$，有 $\sup_{t \in J} \max\{ \|\Phi_t\|_{\dot{C}^{0,1}}, \|(\Phi_t)^{-1}\|_{\dot{C}^{0,1}} \} < \infty$；而对于任一端点 $T$，要么 $|T| = \infty$，要么 $\lim_{t \to T} L_\mu(\Phi_{t*} \Theta) = \infty$。

注：条件 $L_\mu(\Theta) < \infty$ 蕴含 $\theta_0 \in C^{2\alpha}(\mathbb{R}^2)$，且足以推出 $u(\theta_0)$ 为Lipschitz，从而保障流的存在性。

定理 1.5（具有 $H^2$ 水平集的低正则性解的适定性）

设 $\theta_0$ 有一个广义层蛋糕表示 $(\Theta, \mu)$，其中每个 $\Theta_\lambda$ 是以 $z_\lambda \in \mathrm{PSC}(\mathbb{R}^2)$（正定向简单闭曲线等价类）为边界的有界开集，且满足：

  • $L_\mu(\Theta) < \infty$,
  • $R_\mu(\Theta) := \sup_{\lambda} \ell(z_\lambda)^{1/2} \int_L \frac{d|\mu|(\lambda')}{\ell(z_{\lambda'})^{1/2} \Delta(z_\lambda, z_{\lambda'})^{2\alpha}} < \infty$,
  • $Q(\Theta) := \sup_{\lambda} \ell(z_\lambda) \|z_\lambda\|_{\dot{H}^2}^2 < \infty$,

其中 $\ell(\gamma)$ 表示曲线长度，$\|\gamma\|_{\dot{H}^2}$ 为其 $L^2$ 意义下的曲率范数，$\Delta(\gamma_1, \gamma_2) = d(\mathrm{im}(\gamma_1), \mathrm{im}(\gamma_2))$ 为像集之间的距离。

令 $\theta$ 为定理1.3给出的解。则：

(i) 对所有 $t \in I$，$\Phi_{t*} \Theta$ 仍由简单闭曲线组成，且对任意紧 $J \subset I$，$\sup_{t \in J} R_\mu(\Phi_{t*} \Theta) < \infty$。

(ii) 存在包含0的子区间 $I' \subset I$，使得对任意紧 $J \subset I'$，$\sup_{t \in J} Q(\Phi_{t*} \Theta) < \infty$。令 $I'$ 为极大这样的区间，则其端点 $T$ 满足：若 $|T| < \infty$，则必有 $\lim_{t \to T} L_\mu(\Phi_{t*} \Theta) = \infty$ 或 $\lim_{t \to T} Q(\Phi_{t*} \Theta) = \infty$。

(iii) 若进一步假设 $\alpha \leq \frac{1}{6}$，$|\mu|(L) < \infty$，$\int_L \ell(\partial \Theta_\lambda) d|\mu|(\lambda) < \infty$，以及 $\inf_{\lambda \in L} |\Theta_\lambda| > 0$，则对上述 $I'$ 的任意有限端点 $T$，必有 $\limsup_{t \to T} Q(\Phi_{t*} \Theta) = \infty$。

结论(iii)表明：只要初始水平集具有严格正面积且其他积分条件满足，解不能仅因水平集“靠近”或“堆积”而破裂；任何有限时间爆破都必须伴随着某个水平集失去 $H^2$ 正则性（即曲率 $L^2$ 范数趋于无穷）。

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

a) 核心分析工具

本文主要使用Lagrangian方法结合逼近方案来构建解。基本思路是通过光滑化核 $K_\varepsilon(x) = \chi(\varepsilon^{-1}|x|)K(x)$ 构造正则化速度场 $u_\varepsilon(\theta)$，求解相应的ODE系统 $\partial_t \Phi_t^\varepsilon = u_\varepsilon(\theta_t^\varepsilon) \circ \Phi_t^\varepsilon$，然后取极限 $\varepsilon \to 0$。为了控制收敛性，作者依赖以下几个关键估计：

  • Lipschitz性估计：引理2.1证明了 $ \|D(u_\varepsilon(\theta))\|_{L^\infty} \leq C_\alpha L_\mu(\Theta) $，这是整个适定性理论的基础。
  • 稳定性估计：引理2.2给出了不同配置下速度差的控制：
    $$ \|u(\Phi_{1*}\theta) - u(\Phi_{2*}\theta)\|_{L^\infty} \leq C_\alpha (L_{\mu}(\Phi_{1*}\Theta) + L_{\mu}(\Phi_{2*}\Theta)) \|\Phi_1 - \Phi_2\|_{L^\infty}, $$
    这使得可以通过压缩映射原理在适当的函数空间中建立局部解的存在唯一性。
  • 误差估计：引理2.3给出 $ \|u(\theta) - u_\varepsilon(\theta)\|_{L^\infty} \leq C_\alpha \|\theta\|_{L^\infty} \varepsilon^{1 - 2\alpha} $，用于验证逼近解收敛到真实解。

b) 关键技术与难点处理

最大挑战在于如何在整体函数正则性极低的情况下，维持水平集几何结构的稳定性。作者的关键洞见是将适定性问题转化为对三个量 $L_\mu, R_\mu, Q$ 的一致控制：

  • $L_\mu(\Theta)$ 控制整体Hölder正则性；
  • $R_\mu(\Theta)$ 控制不同尺度曲线间的“拥挤程度”，防止小尺度曲线密集靠近大尺度曲线导致奇性积累；
  • $Q(\Theta)$ 直接衡量所有水平集的 $H^2$ 正则性强度。

通过追踪这些量在流作用下的演化（利用 $\Phi_{t*}\Theta$ 的变换性质），作者能够证明：只要 $Q(\Phi_{t*}\Theta)$ 保持有界，水平集就不会坍缩或产生角点，从而避免非紧性或拓扑变化带来的困难。

特别地，在证明 $Q(\Phi_{t*}\Theta)$ 的局部有界性时，需精细估计曲率沿流的传播行为。虽然文中未显式写出曲率方程，但从 $Q$ 的定义可知，其控制等价于控制所有 $\|\kappa_\lambda(t)\|_{L^2(ds)}^2$ 的一致上界，其中 $\kappa_\lambda(t)$ 是 $\Phi_t(\partial \Theta_\lambda)$ 的曲率。这一控制依赖于速度场在曲线上的正则性提升，而这正是 $H^2$ 假设所提供的几何优势。

此外，作者巧妙利用了广义层蛋糕表示的灵活性：允许非嵌套集合族，使他们可以处理复杂拓扑结构（如嵌套patch或多重连通水平集），而不局限于单调层蛋糕。

与之前工作的比较

a) 突破与改进

本文的主要贡献体现在以下几方面：

  1. 显著降低整体正则性要求：相比传统理论要求 $\theta_0 \in C^\beta$ 且 $\beta > 2\alpha$，本文允许 $\theta_0$ 仅处于 $C^{2\alpha}$ 边缘（甚至略弱），仅需满足积分型模条件。这是对经典局部适定性结果 [4] 的边际改进（见定理1.3注释）。

  2. 首次建立具有 $H^2$ 水平集的连续解的适定性理论：以往关于patch解的工作 [7,8,11] 集中于分片常数函数，而本文处理的是取值连续、但水平集边界规则的情形。这填补了高正则连续解与patch解之间的空白。

  3. 揭示奇点形成的几何本质：定理1.5(iii) 表明，当 $\alpha \leq 1/6$ 且初始数据满足一定积分条件时，解只能在水平集失去 $H^2$ 正则性时停止存在，不可能仅仅因为水平集相互靠近或堆积而破裂。这排除了所谓“pile-up”机制作为有限时间爆破的原因，深化了对g-SQG动力学的理解。

  4. 推广到任意几何构型：不同于以往研究水平集接近速率的工作 [6] 多限于特殊对称设置（如图形状水平集），本文适用于任意 $H^2$ 简单闭曲线组成的水平集族，无论其相对位置或拓扑如何。

总之，本文通过引入几何导向的函数空间（基于层蛋糕与 $H^2$ 曲线），发展出一套新的适定性框架，不仅拓宽了解的存在范围，而且明确了真正导致解破裂的内在机制——即几何正则性的崩溃，而非单纯的函数值聚集。

4. Global weak solutions of 3D compressible magnetohydrodynamic equations subject to large external potential forces with discontinuous initial data and vacuum

具有间断初值和真空的三维可压缩磁流体动力学方程在大外部势力作用下的全局弱解

作者: Geyuan Chen, Xin Zhong

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16121v1

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

详细分析

中文标题: 具有间断初值和真空的三维可压缩磁流体动力学方程在大外部势力作用下的全局弱解

研究对象和背景

本文研究的是三维有界区域中带有大外部势力、不连续初始数据以及允许真空状态存在的可压缩磁流体动力学（Magnetohydrodynamics, MHD）方程组的全局弱解存在性问题。所考虑的系统是一个耦合了流体力学与电磁场演化的偏微分方程组，描述导电流体（如等离子体）在磁场影响下的运动规律。

该系统的数学形式如下，在一个光滑、单连通的有界区域 $\Omega \subset \mathbb{R}^3$ 中：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho u) = 0, \\ (\rho u)_t + \mathrm{div}(\rho u \otimes u) - \mu \Delta u - (\mu + \lambda)\nabla \mathrm{div}\,u + \nabla P = (\nabla \times H) \times H + \rho f, \\ H_t - \nabla \times (u \times H) + \nu \nabla \times (\nabla \times H) = 0, \\ \mathrm{div}\,H = 0, \end{cases} \tag{1.1} $$

其中未知函数分别为：密度 $\rho = \rho(x,t)$、速度向量场 $u = (u_1,u_2,u_3)(x,t)$、压力 $P = A\rho^\gamma$（$A>0$, $\gamma > 1$），以及磁场向量场 $H = (H_1,H_2,H_3)(x,t)$。粘性系数满足物理条件：

$$ \mu > 0, \quad 2\mu + 3\lambda \geq 0, $$

而 $\nu > 0$ 是电阻率系数。外力项假设为保守力，即 $f = \nabla \phi(x)$，其中 $\phi(x)$ 是给定的空间依赖势函数。

初始条件为：

$$ (\rho, \rho u, H)(x,0) = (\rho_0, \rho_0 u_0, H_0)(x), \quad x \in \Omega, \tag{1.2} $$

并配备Navier-slip 边界条件：

$$ u \cdot n = 0, \quad \mathrm{curl}\,u \times n = 0, \quad H \cdot n = 0, \quad \mathrm{curl}\,H \times n = 0, \quad x \in \partial\Omega, \; t > 0, \tag{1.3} $$

其中 $n$ 是边界 $\partial\Omega$ 上的单位外法向量。

值得注意的是，这些边界条件比常见的 Dirichlet 条件更复杂，它们保留了切向滑移的可能性，同时保证了某些高阶算子的良好定义性。特别地，$(\mathrm{curl}\,u) \times n = 0$ 和 $(\mathrm{curl}\,H) \times n = 0$ 对分析 vorticity 和 magnetic current 的行为至关重要。

本文允许初始密度 $\rho_0$ 出现真空状态（即 $\rho_0 = 0$ 在某些点或集合上成立）且具有大的振荡（large oscillations），这极大地增加了系统在数学上的奇异性与退化性。此外，外力 $f = \nabla\phi$ 被允许是“大的”，这意味着不能通过小参数摄动方法处理，必须发展新的估计技术来克服其对解结构的强烈扰动。

为了构建平衡态，作者引入稳态解 $(\rho_s, 0, 0)$ 满足静止方程：

$$ \nabla P(\rho_s) = \rho_s \nabla \phi, \quad \int_\Omega \rho_s dx = \int_\Omega \rho_0 dx. \tag{1.9} $$

在适当正则性条件下（$\phi \in H^2$ 且满足某种平均不等式），这样的 $\rho_s$ 存在唯一，并满足 $0 < \underline{\rho} \leq \rho_s \leq \overline{\rho}$，从而作为参考密度用于构造能量泛函。

系统的总初始能量定义为：

$$ C_0 \triangleq \int_\Omega \left( \frac{1}{2} \rho_0 |u_0|^2 + G(\rho_0, \rho_s) + \frac{1}{2} |H_0|^2 \right) dx, \tag{1.10} $$

其中 $G(\rho, \rho_s)$ 是相对于稳态的势能密度：

$$ G(\rho, \rho_s) \triangleq \int_{\rho_s}^\rho \frac{P(\xi) - P(\rho_s)}{\xi^2} d\xi, $$

并满足控制关系：

$$ C^{-1}(\rho - \rho_s)^2 \leq G(\rho, \rho_s) \leq C(\rho - \rho_s)^2. $$

本文的核心挑战在于：即使在没有真空的情况下，可压缩 MHD 方程本身就因强非线性、速度-磁场耦合及约束条件 $\mathrm{div}\,H=0$ 而极为困难；一旦允许真空出现，密度可能趋于零，导致质量方程退化、动量方程中对流项失去控制，传统 Sobolev 嵌入不再有效。再加上物理边界的出现和大外力的影响，使得以往在全空间或周期情形下发展的方法难以直接应用。

因此，本文的研究对象不仅是经典的 MHD 系统，而且是在真实物理边界环境下、面对高度非正则初始数据（含真空、间断）、受强外场驱动的一类极端情形，目标是建立此类条件下全局弱解的存在性。

主要定理或主要结果

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

Theorem 1.1. 设 $\Omega \subset \mathbb{R}^3$ 是一个单连通的有界区域，其光滑边界 $\partial\Omega$ 具有有限多个二维连通分支。设势函数 $\phi$ 满足条件：

$$ \phi \in H^2, \quad \int_\Omega \left( \phi - \inf_\Omega \phi \right) dx < \frac{A\gamma}{\gamma - 1} \int_\Omega \rho_0 dx, \tag{1.6} $$

此时存在唯一的稳态密度 $\rho_s$ 满足 (1.9)，且 $0 < \underline{\rho} \leq \rho_s \leq \overline{\rho}$。

对于任意给定的正常数 $M_1, M_2$ 及 $\hat{\rho} \geq \overline{\rho} + 1$，假设初始数据 $(\rho_0, u_0, H_0)$ 满足：

$$ 0 \leq \inf_\Omega \rho_0 \leq \sup_\Omega \rho_0 \leq \hat{\rho}, \quad (u_0, H_0) \in H^1_\omega, \quad \mathrm{div}\,H_0 = 0, $$

$$ \|\nabla u_0\|_{L^2}^2 \leq M_1, \quad \|\nabla H_0\|_{L^2}^2 \leq M_2, \tag{1.16} $$

其中 $H^1_\omega \triangleq \{ v \in H^1 \mid (v\cdot n)|_{\partial\Omega} = (\mathrm{curl}\,v \times n)|_{\partial\Omega} = 0 \}$。

则存在仅依赖于 $\mu, \lambda, \nu, \gamma, A, \hat{\rho}, \Omega, \|\phi\|_{H^2}, \inf_\Omega \phi, M_1, M_2$ 的正常数 $\varepsilon$，使得当初始能量满足：

$$ C_0 \leq \varepsilon, \tag{1.17} $$

时，初边值问题 (1.1)–(1.3) 存在一个全局弱解 $(\rho, u, H)$，在 Definition 1.1 意义下成立，并满足：

$$ 0 \leq \rho(x,t) \leq 2\hat{\rho}, \quad \forall x \in \Omega, \; t \geq 0. \tag{1.18} $$

进一步地，对任意 $r \in [1,\infty), p \in [1,6]$，存在常数 $C, \eta > 0$，使得当 $r \in (2,\infty), p \in (2,6]$ 时有：

$$ \lim_{t \to \infty} \int_\Omega \left( |\rho - \rho_s|^r + \sqrt{\rho}|u|^4 + |H|^p \right) dx = 0. \tag{1.19} $$

此定理表明：只要初始能量足够小，即便初始密度含有真空、允许大幅震荡、且存在大外部势力，仍可在三维有界域中构造出满足 Navier-slip 边界条件的全局弱解，并揭示了解随时间趋于稳态的趋势。

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

a) 核心分析工具

论文采用了基于先验估计+紧性方法的经典框架，结合现代流体力学中处理低正则性问题的技术路线。主要工具包括：

  • 有效粘性通量（Effective Viscous Flux）：定义为

    $$ F \triangleq (2\mu + \lambda)\mathrm{div}\,u - (P - P(\rho_s)) - \frac{1}{2}|H|^2. $$

    这一组合在压缩流中起着关键作用，它将压力梯度与速度散度联系起来，使动量方程可重写为：

    $$ \rho \dot{u} - (\rho - \rho_s)\nabla\phi = \nabla F - \mu \mathrm{curl}\,\mathrm{curl}\,u + H \cdot \nabla H. $$

    利用这一表达式，结合椭圆理论，可以推导出 $\nabla F$ 的 $L^p$ 估计，进而控制 $\nabla u$ 的高阶范数。

  • Bogovskii 算子：用于处理不可压约束下的压力估计。由于边界非平坦，标准傅里叶方法失效，Bogovskii 提供了一个从 $L^p$ 到 $W^{1,p}_0$ 的线性算子 $B[f]$，满足 $\mathrm{div}\,B[f] = f$，从而可用于测试函数构造以获得 $\|P\|_{L^2}$ 控制。

  • Lamé 系统正则性理论：针对形如 $-\mu\Delta u - (\mu+\lambda)\nabla\mathrm{div}\,u = f$ 的弹性系统，在 Navier-slip 边界下，利用引理 2.3 得到 $W^{k+2,q}$ 正则性估计，这对闭合先验估计至关重要。

  • Gagliardo–Nirenberg 不等式与向量场分解：借助引理 2.1 和 2.4–2.5，在单连通区域内利用 $\|\nabla v\|_{L^q} \leq C(\|\mathrm{div}\,v\|_{L^q} + \|\mathrm{curl}\,v\|_{L^q})$ 的估计技巧，避免使用迹定理带来的边界项困难。

  • Zlotnik 不等式（Lemma 2.8）：这是获得密度一致上界的利器。若某函数 $y(t)$ 满足 $y'(t) = g(y) + b'(t)$，且 $g(\infty) = -\infty$，$b(t_2)-b(t_1) \leq N_0 + N_1(t_2-t_1)$，则 $y(t)$ 有界。本文将其应用于密度演化，结合 $F$ 的结构性质，最终证明 $\rho \leq 2\hat{\rho}$。

b) 关键技术难点与应对策略

  1. 大外部势力的处理：传统小外力可通过摄动处理，但本文中 $\nabla\phi$ 可任意大。解决方式是将外力效应吸收进稳态密度 $\rho_s$，并通过定义相对能量 $G(\rho, \rho_s)$ 将系统“中心化”于平衡态附近。再利用 $F$ 的结构消去 $\nabla\phi$ 项，转化为对 $(\rho - \rho_s)$ 的控制。

  2. 边界积分项的消除：Navier-slip 条件下会出现难以估计的边界积分项，如 $\int_{\partial\Omega} F (u_t \cdot n)\,dS$。作者利用 $u\cdot n = 0$ 推出 $u\cdot\nabla u\cdot n = -u\cdot\nabla n\cdot u$，从而得到 $(\dot{u} + (u\cdot\nabla n)\times u_\perp)\cdot n = 0$，规避了一阶导数出现在边界上的问题。

  3. 真空导致的退化性：当 $\rho \to 0$ 时，$\rho u_t$ 项消失，常规能量方法失效。作者采用加权时间函数 $\sigma(t) = \min\{1,t\}$ 构造加权能量泛函 $A_1(T), A_2(T)$，并结合 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 的控制来维持动量变化的有效性。

  4. 非线性项的精细估计：特别是磁场相关项如 $H\cdot\nabla H$、$(\nabla\times H)\times H$，需利用 Hölder 不等式与 Sobolev 嵌入进行协调估计。例如，

    $$ \|H\cdot\nabla H\|_{L^p} \leq \|H\|_{L^\infty} \|\nabla H\|_{L^p}, $$

    再结合 $\|H\|_{L^\infty} \leq C(\|H\|_{L^q} + \|\nabla H\|_{L^r})$ 完成控制。

  5. 逼近解的构造与极限过程：通过正则化初始密度（$\rho_0^\delta \geq \delta > 0$）构造局部经典解，利用统一先验估计延拓至全局，最后令 $\delta \to 0$ 取弱极限。紧性来自各量在适当空间中的有界性，结合 Aubin-Lions 引理和弱收敛方法完成极限过渡。

综上，本文通过巧妙构造有效粘性通量、合理使用加权能量泛函、严格控制边界效应与非线性耦合项，成功突破了真空、大外力、物理边界三重障碍。

与之前工作的比较

a) 突破与改进之处

本文相较于已有工作有多方面显著推进：

  1. 从全空间推广到有界域：此前关于中间型弱解的工作（如 Suen–Hoff [27], Chen–Suen [20]）均限于 $\mathbb{R}^3$ 或周期情形。本文首次将这类结果扩展到具有真实物理边界的有界区域，必须处理大量由边界引起的表面积分和技术障碍，属于非平凡推广。

  2. 允许大外部势力：多数文献要么忽略外力，要么要求其充分小。而本文明确允许 $\nabla\phi$ 为“大”的，仅需满足积分型条件 (1.6)。这种设定更符合地球物理与天体物理中旋转主导流动的实际场景（如大气环流），具有更强的应用价值。

  3. 放宽初始数据条件：相比早期要求初始密度远离真空（如 Matsumura–Nishida 类解），本文完全允许 $\rho_0 = 0$ 出现，并接受大的密度振荡。这使得模型适用于激波、空腔形成等极端现象。

  4. 弱解类别的强化：不同于 Lions–Feireisl 类的有限能量弱解（正则性较低），也区别于小扰动下的经典解，本文所得弱解属于“中间型”——虽不要求密度光滑，但仍能定义粒子轨迹（在非真空区），且具备更好的衰减性质（见 (1.19)）。这填补了强解与弱解之间的空白。

  5. 边界条件的新处理方式：采用 Navier-slip 而非 Dirichlet 条件，更适合滑移界面建模。同时，文中发展了一套完整的边界项估计技术，为后续研究类似边界问题提供了范式。

总之，本文不仅在理论上统一并推广了多个先前结果，还在技术层面实现了对多重困难的同时克服，标志着可压缩 MHD 在复杂环境下的适定性理论取得重要进展。

5. Low Regularity Well-Posedness of Cauchy Problem for Two-Dimensional Relativistic Euler Equation

二维相对论欧拉方程柯西问题的低正则性适定性

作者: Huali Zhang

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16090v1

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

详细分析

中文标题: 二维相对论欧拉方程柯西问题的低正则性适定性

研究对象和背景

本文研究的是二维空间中相对论欧拉方程（relativistic Euler equations）的柯西问题在低正则性初始数据下的适定性。该系统描述了在狭义相对论框架下，可压缩流体（如高能密度物质、天体物理中的等离子体）的动力学演化。

系统的原始形式基于能量-动量张量守恒律。设 $%$ 表示能量密度，$u^\alpha$ 表示四维速度，满足归一化条件 $u_\alpha u^\alpha = -1$（其中指标升降由平直Minkowski度规 $m_{\alpha\beta} = \mathrm{diag}(-1,1,1)$ 完成）。令 $Q^{\alpha\beta} = (p + %)u^\alpha u^\beta + p m^{\alpha\beta}$ 为应力-能量张量，其守恒给出运动方程：

$$ \partial_\alpha Q^{\alpha\beta} = 0. $$

将此方程投影到与 $u^\alpha$ 平行和正交的方向上，得到如下系统（见原文(1.2)）：

$$ \begin{cases} u^\kappa \partial_\kappa % + (p+%) \partial_\kappa u^\kappa = 0, \\ (p+%) u^\kappa \partial_\kappa u^\alpha + (m^{\alpha\kappa} + u^\alpha u^\kappa) \partial_\kappa p = 0. \end{cases} $$

文中考虑的状态方程为 $p = p(%) = %^A$，其中 $A \geq 1$ 是绝热常数。当 $A > 1$ 时称为多变流体（polytropic fluid），当 $A=1$ 时称为刚硬流体（stiff fluid），此时声速等于光速，具有特殊的几何结构。

为了研究低正则性解，作者引入了一组“好变量”（good variables），包括：

  • 对数焓（logarithmic enthalpy）：$h = \log H$，其中 $H = (p+%)/q$ 是每粒子焓，$q$ 是数密度。
  • 重缩放速度（rescaled velocity）：$\mathbf{v} = e^h \mathbf{u}$，这使得速度场在低正则性分析中表现出更好的代数结构。
  • 声速：$c_s = \sqrt{dp/d%}$，依赖于状态方程。

通过这些变换，原系统被重写为一个耦合的波-输运系统。关键的新变量是声学度规 $g_{\alpha\beta}$ 和涡量 $\mathbf{w}_\alpha$。声学度规定义为：

$$ g_{\alpha\beta} = \Theta \left( c_s^2 m_{\alpha\beta} + e^{-2h}(c_s^2 - 1) v_\alpha v_\beta \right), $$

其中 $\Theta = \left[ c_s^2 - e^{-2h}(c_s^2 - 1)(v^0)^2 \right]^{-1}$。该度规决定了扰动传播的因果结构（即声锥）。涡量定义为：

$$ w_\alpha = \epsilon_{\alpha\beta\gamma} \partial^\beta v^\gamma, $$

其中 $\epsilon_{\alpha\beta\gamma}$ 是Levi-Civita符号。这是二维情形下自然的涡量定义，不同于三维的旋度。

最终，系统被转化为如下形式（见原文(1.19)）：

$$ \begin{cases} \square_g h = D, \\ \square_g v_\alpha = -c_s^2 \Theta \epsilon_{\alpha\beta\gamma} \partial^\beta w^\gamma + Q_\alpha, \\ v^\kappa \partial_\kappa w_\alpha = w^\kappa \partial_\alpha v_\kappa - w_\alpha \partial_\kappa v^\kappa, \end{cases} $$

其中 $\square_g = g^{\alpha\beta} \partial^2_{\alpha\beta}$ 是关于声学度规 $g$ 的达朗贝尔算子，$D$ 和 $Q_\alpha$ 是关于 $(h, \mathbf{v})$ 及其导数的二次项（见原文(1.20)-(1.21)），右侧的输运方程描述了涡量沿流线的演化。

本文的研究区域是全平面 $\mathbb{R}^2$，无边界，因此是全空间上的柯西问题。初始数据在 $t=0$ 给出，属于特定的Sobolev空间和Besov空间。特别地，作者关注的是低于经典局部适定性理论要求的低正则性（low regularity）情形。

关键的函数空间包括：

  • $H^s(\mathbb{R}^2)$：标准Sobolev空间，范数为 $\|f\|_{H^s} = \|f\|_{L^2} + \|f\|_{\dot{H}^s}$，其中齐次范数 $\|f\|_{\dot{H}^s}^2 = \sum_{j \in \mathbb{Z}} 2^{2js} \|P_j f\|_{L^2}^2$，$P_j$ 是Littlewood-Paley投影。
  • $L^8(\mathbb{R}^2)$：用于控制 $\nabla \mathbf{w}_0$，这是一个非标准但技术上必要的条件，源于交换子估计（commutator estimate）的需求。
  • $\dot{B}^1_{2,1}(\mathbb{R}^2)$：齐次Besov空间，其范数为 $\|f\|_{\dot{B}^1_{2,1}} = \sum_{j \in \mathbb{Z}} 2^j \|P_j f\|_{L^2}$，出现在小初值全局适定性的结果中。

一个重要的假设是初始数据满足有界性和严格双曲性条件（见原文(1.23)）：

$$ |h_0, \mathbf{v}_0| \leq C_0, \quad 0 < c_0 < c_s < 1, $$

这保证了声学度规在整个演化过程中保持时间类，并且系统是严格双曲的。

主要定理或主要结果

本文建立了四个主要定理，分别对应不同正则性和物理情形下的适定性结果。

定理 2.1（一般多变流体，高正则性涡量）
设 $7/4 < s \leq 15/8$，$s_0 > 3/2$。给定初值 $(h_0, \mathbf{v}_0, \mathbf{w}_0)$ 满足(1.23)且

$$ \|(h_0, \mathbf{v}_0)\|_{H^s} + \|\mathbf{w}_0\|_{H^{s_0}} + \|\nabla \mathbf{w}_0\|_{L^8} \leq M_0, $$

则存在 $T>0$ 和 $M_1>0$，使得柯西问题(1.19)在 $[0,T]$ 上存在唯一解 $(h,\mathbf{v}) \in C([0,T];H^s) \cap C^1([0,T];H^{s-1})$，$\mathbf{w} \in C([0,T];H^{s_0}) \cap C^1([0,T];H^{s_0-1})$。解满足能量估计、Strichartz估计以及线性化方程的一致估计。

定理 2.2（一般多变流体，临界正则性涡量）
设 $7/4 < s \leq 15/8$。给定初值 $(h_0, \mathbf{v}_0, \mathbf{w}_0)$ 满足(1.23)且

$$ \|(h_0, \mathbf{v}_0)\|_{H^s} + \|\mathbf{w}_0\|_{H^{3/2}} + \|\nabla \mathbf{w}_0\|_{L^8} \leq M_0, $$

则存在 $T^*>0$ 和 $M_2>0$，使得问题(1.19)在 $[0,T^*]$ 上存在唯一解，且解满足类似定理2.1的能量和Strichartz估计。此定理的关键在于将涡量正则性从 $H^{3/2+}$ 放宽至 $H^{3/2}$。

定理 2.3（刚硬流体，一般涡量）
设 $p(%)=%$（即 $A=1$），初值 $(h_0, \mathbf{v}_0, \mathbf{w}_0) \in H^{7/4+\epsilon}(\mathbb{R}^2) \times H^{7/4+\epsilon}(\mathbb{R}^2) \times H^{1+\epsilon'}(\mathbb{R}^2)$，其中 $\epsilon' < \epsilon \leq 1/8$。则柯西问题(2.7)局部适定，且 $(dh, d\mathbf{v}) \in L^4 L^\infty \cap L^4 \dot{B}^{\epsilon'}_{\infty,2}$。

定理 2.4（刚硬流体，无旋流动）
设 $p(%)=%$ 且初值涡量为零（$\mathbf{w}_0 \equiv 0$）。若 $(h_0, \mathbf{v}_0) \in H^{1+\epsilon}(\mathbb{R}^2)$，则问题(2.8)局部适定。进一步，若 $(h_0, \mathbf{v}_0) \in \dot{B}^1_{2,1}(\mathbb{R}^2)$ 且足够小，则解全局存在。

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

a) 核心分析工具

本文的核心方法是将Smith-Tataru关于拟线性波动方程的低正则性理论推广到相对论欧拉方程这一更复杂的耦合系统。主要工具包括：

  1. Strichartz估计：这是证明低正则性适定性的基石。对于线性波动方程 $\square_g f = 0$，当系数（即度规 $g$）本身不够光滑时，经典的Strichartz估计不再成立。Smith-Tataru [61] 的突破在于，即使 $g \in H^{7/4+}$（在二维），仍可通过波包分解（wave packet decomposition）和特征超曲面（characteristic hypersurfaces）分析，建立弱形式的Strichartz估计，例如 $\|f\|_{L^4_t L^\infty_x} \lesssim \|f[0]\|_{H^r}$ 对某些 $r > 7/4$ 成立。本文中，作者将此方法应用于声学度规 $g$ 依赖于未知函数 $h$ 和 $\mathbf{v}$ 的情形。

  2. 能量估计：用于控制解在 $H^s$ 空间中的大小。由于系统是非线性的，需要设计适当的能量泛函并进行迭代。能量估计通常结合Gronwall不等式来获得解的时间存在区间。

  3. 半经典分析（semiclassical analysis）：这是处理定理2.2中涡量正则性降至 $H^{3/2}$ 的关键技术。受Bahouri-Chemin [9] 和 Ai-Ifrim-Tataru [2] 启发，作者利用频率截断构造近似解序列。在每个频率尺度上，Strichartz估计可能带有依赖于频率的常数，导致求和时出现导数损失（loss of derivatives）。通过精细分析这种损失，可以证明近似解在较低正则性空间中收敛，从而得到原方程的解。

  4. 椭圆估计（elliptic estimates）：在二维情形下，作者引入了新的椭圆估计来处理涡量 $\mathbf{w}$ 与速度 $\mathbf{v}$ 之间的关系。特别是，通过定义算子 $P$（见原文(1.17)），将速度分解为 $v_\alpha = v_{+\alpha} + v_{-\alpha}$，其中 $v_-$ 由涡量决定，$v_+$ 则更规则。这种分解允许将 $v_+$ 的正则性与 $h$ 和 $w$ 联系起来。

  5. Bootstrap论证（bootstrap argument）：这是一种先验假设再证明自洽的方法。假设解在某个时间段内满足一定的正则性估计，然后利用方程和各种不等式证明这些估计实际上可以被改进或保持，从而延长解的存在时间。

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

  1. 克服二维涡量结构的困难：在三维中，标准的 $u$-涡量公式能自然产生有效的输运结构。但在二维，直接推广会失效。作者借鉴[17,19]的思想，通过引入重缩放速度 $\mathbf{v}$ 和对数焓 $h$，揭示了一个隐藏的波-输运结构，并定义了适合二维的涡量 $\mathbf{w}_\alpha = \epsilon_{\alpha\beta\gamma} \partial^\beta v^\gamma$，从而解决了这一结构性难题。

  2. 处理涡量正则性临界下降：从 $H^{3/2+}$ 到 $H^{3/2}$ 是一个临界下降。直接应用Smith-Tataru方法会失败，因为无法获得速度的Strichartz估计。作者的技巧是使用频率截断：将初值 $\mathbf{w}_0 \in H^{3/2}$ 截断为光滑序列 $\mathbf{w}_0^{(N)} \in H^{3/2+}$。对每个 $N$，应用定理2.1得到解 $(h^{(N)}, \mathbf{v}^{(N)})$ 在短时间 $T_N$ 内存在。关键是分析 $T_N$ 如何依赖于 $N$。利用半经典分析，发现Strichartz估计在每个频段上有可控的损失，使得总和仍收敛。这允许将解延拓到一个与 $N$ 无关的统一时间 $T^*$，并证明序列在适当空间中收敛。

  3. 利用特殊状态方程的简化：当 $p(%)=%$ 时，声速 $c_s \equiv 1$，声学度规退化为平直Minkowski度规 $g=m$。这使得波动算子 $\square_g$ 成为标准的 $\square$，其Strichartz估计是已知的（无需依赖系数的正则性）。这极大地简化了分析，允许将对数焓 $h$ 和速度 $\mathbf{v}$ 的正则性降至 $H^{7/4+}$，同时将涡量正则性进一步降至 $H^{1+}$。

  4. 无旋情形的“波映照”结构：当 $\mathbf{w} \equiv 0$ 时，系统(2.8)中出现了类似于“波映照”（wave map）的零条件（null condition）的非线性项（如 $-2\partial_\beta h \partial^\beta h$）。这类非线性在能量估计中可以被更好地控制，从而允许更低的局部正则性要求 $H^{1+}$，甚至在小初值下达到Besov空间 $\dot{B}^1_{2,1}$ 中的全局适定性。

与之前工作的比较

a) 突破与改进

本文的主要突破体现在以下几个方面：

  1. 首次建立二维相对论欧拉方程的低正则性适定性：尽管三维情形已有工作（如Disconzi-Speck [19], Zhang [78]），但二维情形因其独特的涡量动力学和Strichartz估计的不同而长期未被解决。本文填补了这一空白，是该领域的开创性工作。

  2. 显著放宽了对初始涡量的正则性要求：与之前的许多工作相比，本文对涡量的要求更低。例如，在Newtonian可压Euler方程中，Zhang [76,77] 要求涡量在 $H^{3/2}$ 以上。本文的定理2.2成功地将要求降至临界空间 $H^{3/2}$，这是一个实质性的改进。这得益于半经典分析和频率求和技巧的有效应用。

  3. 推广了Smith-Tataru方法的应用范围：Smith-Tataru [61] 的理论原本针对单个拟线性波动方程。本文将其成功应用于一个包含波方程和输运方程的耦合系统，并处理了度规依赖于多个未知函数（$h$ 和 $\mathbf{v}$）的复杂情况，展示了该方法的强大普适性。

  4. 揭示了特殊物理状态的优越数学结构：当 $A=1$（刚硬流体）时，本文指出声学度规变为平直，从而将问题降为半线性系统。这不仅使结果更强（更低的正则性要求），也凸显了物理模型选择对数学分析难度的深刻影响。此前的工作并未在相对论背景下充分探讨这一特例的数学优势。

  5. 实现了与经典结果的最优匹配：在刚硬流体情形，定理2.3中 $h$ 和 $\mathbf{v}$ 的正则性 $H^{7/4+}$ 与Smith-Tataru [61] 关于拟线性波方程的结果一致，而涡量的 $H^{1+}$ 正则性与Bourgain-Li [12] 关于不可压Euler方程的结果相呼应。这表明本文的结果在各自框架下是最优的（sharp）。

综上所述，本文通过精巧的变量变换、对先进分析工具的娴熟运用以及对物理模型的深刻理解，在二维相对论流体力学的低正则性理论中取得了重要进展，其结果在正则性要求上比前人工作更为宽松，且在特殊情形下达到了最优。

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

1. Global existence and stability of near-affine solutions of compressible elastodynamics

可压缩弹性动力学近仿射解的整体存在性与稳定性

作者: Xianpeng Hu, Yuanzhi Tu, Changyou Wang, Huanyao Wen

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16505v1

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

摘要翻译

中文标题: 可压缩弹性动力学近仿射解的整体存在性与稳定性

摘要翻译: 本文证明了对于 $ \mathbb{R}^d $（其中 $ d=2,3 $）上可压缩非线性弹性动力学方程组的柯西问题，当初始扰动为足够小的 $ H^3 $-扰动且接近一个仿射解时，存在唯一的整体强解。此外，我们还建立了该解的渐近行为。

2. Homogenization of a micropolar fluid past a porous media with non-zero spin boundary condition

具有非零自旋边界条件的多孔介质中微极流体的均匀化

作者: Francisco J. Suárez-Grau

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16353v1

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

摘要翻译

中文标题: 具有非零自旋边界条件的多孔介质中微极流体的均匀化

摘要翻译: 本文研究了一种微极流体在由尺寸为 $\varepsilon$ 的周期性分布障碍物穿孔的介质中的流动问题。考虑了非齐次的微旋转边界条件：假设障碍物边界的微旋转与速度场在边界上的旋转率成正比。文中分析了解的存在性与唯一性。此外，当 $\varepsilon$ 趋于零时，通过取极限过程，推导出了多孔介质理论中经典的微极Darcy定律的一个类比形式。

不相关论文（基本信息）

1. Stability under lamination and polycrystalline effective conductivity

层状结构下的稳定性与多晶有效电导率

作者: Nathan Albin, Vincenzo Nesi, Mariapia Palombaro

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16787v1

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

摘要: We prove the stability under lamination of a set of real, symmetric 3$\times$3 matrices that can be viewed as a subset of the effective conductivities of a polycrystal. Constructed in a companion paper, such set in combination with several previous constructions provides the best inner bound known so far on the $G$-closure of a three dimensional polycrystal.

2. Unconditional uniqueness of Hardy–Hénon parabolic equations on Herz spaces

Herz空间上Hardy–Hénon抛物方程解的无条件唯一性

作者: Naoya Hatano, Masahiro Ikeda

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16711v1

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

摘要: In this paper, we introduce the unconditional uniqueness of solutions in Herz spaces for the Hardy–Hénon parabolic equation, which is a semilinear heat equation with a power-type weight in the nonlinear term $|x|^γ|u|^{α-1}u$. It is expected that the power-type weight in the nonlinear term can be effectively handled within Herz spaces. In fact, our result in Herz spaces $\dot{K}^s_{q,r}({\mathbb R}^n)$ relaxes the endpoint case $q=α$ and the large interpolation exponent case $r\ge q$ compared to previous results.

3. The capillary Christoffel-Minkowski problem

毛细现象中的Christoffel-Minkowski问题

作者: Xinqun Mei, Guofang Wang, Liangjun Weng

类别: math.AP, math.DG

发布日期: 2025-12-18

arXiv ID: 2512.16655v1

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

摘要: In this article, we introduce a $k$-th capillary area measure for capillary convex bodies in the Euclidean half-space, which serves as a boundary counterpart to the classical concept of area measure (see, e.g., \cite[Chapter 8]{Sch}). We then propose a Christoffel-Minkowski problem for capillary convex bodies, to find a capillary convex body in the Euclidean half-space with a prescribed $k$-th capillary area measure. This problem is equivalent to solving a Hessian-type equation with a Robin boundary value condition. We then establish the existence and uniqueness of a smooth solution under a natural sufficient condition.

4. Existence and stability of discretely self-similar blowup for a wave maps type equation

一类波映照型方程离散自相似爆破解的存在性与稳定性

作者: Irfan Glogić, David Hilditch, David Wallauch

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

发布日期: 2025-12-18

arXiv ID: 2512.16623v1

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

摘要: We study finite-time blowup for a nonlinear wave equation for maps from the Minkowski space $\mathbb{R}^{1+d}$ into the 1-sphere $\mathbb{S}^1$, whose nonlinearity exhibits a null-form structure. We construct, for every dimension $d \geq 1$, a countable family of discretely self-similar blowup solutions, which are even for $d=1$ and radial for $d \geq 2$. The main contribution of the paper is a detailed nonlinear stability analysis of this family of solutions. For $d \geq 2$, we consider radial data, while in $d=1$ we allow for general perturbations. After linearizing around the self-similar profiles in similarity variables, we construct resolvents of the resulting highly non-self-adjoint operators through Liouville-Green transformations and precise Volterra-type asymptotics. The construction itself, which occupies most of the paper, is technically challenging, as it is performed in arbitrary dimensions and for a countable family of operators in each. Combined with a detailed spectral analysis of the linearized operators, this yields sharp semigroup bounds and allows us to establish nonlinear stability of all discretely self-similar profiles in all dimensions, with precise co-dimension determined by the unstable spectrum. To our knowledge, this is the first result on the existence and stability of discretely self-similar blowup for a geometric wave equation.

5. Manifold submetries from compact homogeneous spaces

紧致齐性空间上的流形次测度

作者: Samuel Lin, Ricardo A. E. Mendes, Marco Radeschi

类别: math.DG, math.AC, math.AP, math.SP

发布日期: 2025-12-18

arXiv ID: 2512.16606v1

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

摘要: We show that singular Riemannian foliations, or, more generally, manifold submetries, defined on a compact normal homogeneous space, have algebraic nature. Moreover, in this case there exists a one-to-one correspondence between algebras of algebraic functions preserved by the Laplace–Beltrami operator, and manifold submetries. A key intermediate result is that, for any manifold submetry on a compact normal homogeneous space, the vector field given by the mean curvature of the fibers is basic, in the sense that it is related to a vector field in the base.

6. Regularity for fully nonlinear elliptic equations in generalized Orlicz spaces

广义Orlicz空间中完全非线性椭圆型方程的正则性

作者: Sun-Sig Byun, Jeongmin Han, Mikyoung Lee

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16600v1

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

摘要: In this paper, we establish an optimal global Calderón-Zygmund type estimate for the viscosity solution to the Dirichlet boundary problem of fully nonlinear elliptic equations with possibly nonconvex nonlinearities. We prove that the Hessian of the solution is as integrable as the nonhomogeneous term in the setting of a given generalized Orlicz space even when the nonlinearity is asymptotically convex with respect to the Hessian of the solution.

7. Constructing steering-type solutions for higher order Cauchy-Riemann equations in $\mathbb{R}^{m+1}$

构造 $\mathbb{R}^{m+1}$ 中高阶柯西-黎曼方程的导引型解

作者: Daniel Alfonso Santiesteban, Dixan Peña Peña, Ricardo Abreu Blaya

类别: math.AP, math-ph

发布日期: 2025-12-18

arXiv ID: 2512.16594v1

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

摘要: The multidimensional Cauchy-Riemann operator provides a framework for studying higher order partial differential equations in $\mathbb{R}^{m+1}$, whose solutions include polymonogenic and polyharmonic functions, among others. In this work, we aim to explicitly construct solutions to such systems, generated from families of complex valued functions which are closed under conjugation and under the action of the complex Cauchy-Riemann operator. Moreover, we prove that precisely some of these solutions also satisfy homogeneous linear differential equations involving the so-called hypercomplex derivative.

8. Liouville-type Theorems for Stable Solutions of the Hénon-Lane-Emden System

关于Hénon-Lane-Emden系统稳定解的刘维尔型定理

作者: Long-Han Huang, Wenming Zou

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16566v1

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

摘要: We investigate the Hénon-Lane-Emden system defined by $- Δu=|x|^a |v|^{p-1}v$ and $- Δv=|x|^b |u|^{q-1}u$ in $\mathbb{R}^N \!\setminus\! \{0\}$. We begin by establishing a general Liouville-type theorem for the subcritical case. Then we prove that the Hénon-Lane-Emden conjecture is valid for solutions stable outside a compact set, provided that $0 < \min\,\{p, q\} < 1$, or $0 \leq a - b \leq (N-2)(p - q)$, or $N \leq \frac{2(p+q+2)}{pq-1} + 10$. Additional Liouville-type theorems for the subcritical case are also obtained. Furthermore, we address the supercritical case. To our knowledge, these results constitute the first Liouville-type theorems for this class of solutions in the Hénon-Lane-Emden system. As a by-product, several existing results in the literature are refined.

9. Ground states for the Hartree energy functional in the critical case

临界情形下Hartree能量泛函的基态解

作者: Tommaso Pistillo

类别: math-ph, math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16513v1

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

摘要: We consider the problem of finding a minimizer $u$ in $ H^1(\mathbb{R}^3)$ for the Hartree energy functional with convolution potential $w$ in $L^\infty(\mathbb{R}^3)+L^{3/2,\infty}(\mathbb{R}^3)$ with $L^\infty$ part vanishing at infinity. This class includes sums of potentials of the kind $-\frac{1}{|x|^α}$, $0<α\le2$, together with the case $w$ in $L^{3/2}(\mathbb{R}^3)$. We prove the existence of such groundstates for a wide range of $L^2$ masses. We also establish basic properties of the groundstates, i.e.~positivity and regularity. Lastly, we exploit the estimates we derived for the stationary problem to prove global well-posedness of the associated evolution problem and orbital stability of the set of ground states.

10. On the distribution kernels of Toeplitz operators on CR manifolds

关于CR流形上Toeplitz算子的分布核

作者: Chin-Yu Hsiao, Ood Shabtai

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

发布日期: 2025-12-18

arXiv ID: 2512.16506v1

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

摘要: We study the distribution kernel of a Toeplitz operator associated with a classical pseudodifferential operator on a compact, embeddable, strictly pseudoconvex CR manifold. The main result consists of a formula for the values at the diagonal of the second coefficient in the expansion of the symbol of the kernel. We also establish asymptotic expansions for Toeplitz operators on the positive part of a compact not necessary strictly pseudoconvex CR orbifold under certain natural assumptions.

11. Normalized solutions for a class of fractional Choquard equations with mixed nonlinearities

一类具混合非线性的分数阶Choquard方程的归一化解

作者: Shaoxiong Chen, Zhipeng Yang, Xi Zhang

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16438v1

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

摘要: In this paper we study the following fractional Choquard equation with mixed nonlinearities:

$$ \left\{ \begin{array}{l} (-Δ)^s u = λu + α\left( I_μ* |u|^q \right) |u|^{q-2} u + \left( I_μ* |u|^p \right) |u|^{p-2} u, \quad x \in \mathbb{R}^N, \\[4pt] \displaystyle \int_{\mathbb{R}^N} |u|^2 \,\mathrm{d}x = c^2 > 0. \end{array} \right. $$

Here $N > 2s$, $s \in (0,1)$, $μ\in (0, N)$, and the exponents satisfy

$$ \frac{2N - μ}{N} < q < p < \frac{2N - μ}{N - 2s}, $$

while $α> 0$ is a sufficiently small parameter, $λ\in \mathbb{R}$ is the Lagrange multiplier associated with the mass constraint, and $I_μ$ denotes the Riesz potential. We establish existence and multiplicity results for normalized solutions and, in addition, prove the existence of ground state normalized solutions for $α$ in a suitable range.

12. Quantitative stratification and optimal regularity for harmonic almost complex structures

调和几乎复结构的定量分层与最优正则性

作者: Chang-Yu Guo, Ming-Lun Liu, Chang-Lin Xiang

类别: math.AP, math.DG

发布日期: 2025-12-18

arXiv ID: 2512.16341v1

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

摘要: In a recent interesting work [15], W.Y. He established the important partial regularity theory and the almost optimal higher regularity theory for energy minimizing harmonic almost complex structures. Based on a new observation on the structure of equations, we give an easier new proof of the partial regularity theorem, and adapting the powerful quantitative stratification method of Naber-Valtorta [22], we further prove the rectifiability of singular stratum of energy minimizing harmonic almost complex structures. Based on this, we establish an optimal regularity theory, which improves the corresponding result of He.

13. Nekhoroshev type stability for Ultra-differential Hamiltonian in $L^2$ space

$L^2$空间中超微分Hamilton系统的Nekhoroshev型稳定性

作者: Bingqi Yu, Li Yong

类别: math.AP, math.DS

发布日期: 2025-12-18

arXiv ID: 2512.16332v1

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

摘要: This paper combines the decay of high modes with the smallness introduced by high orders, leading to a normal form lemma for infinite-dimensional Hamiltonian systems under ultra-differentiable regularity. We prove the sub-exponential stability time of a wide class of Hamiltonian PDEs, including the Schrödinger equation with convolution potentials, fractional-order Schrödinger equations, and beam equations with metrics. When the conditions are equivalent to previous ones, the stability time we obtain reaches Bourgain’s predicted optimal bound. Furthermore, we approach earlier results under lower conditions. These results are discussed within a general framework we propose, which applies to the ultra-differential class.

14. Bifurcating domains for an overdetermined eigenvalue problem in cylinders

柱体中一类超定特征值问题的分岔区域

作者: Yuanyuan Lian, Filomena Pacella, Pieralberto Sicbaldi

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16319v1

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

摘要: We study an overdetermined eigenvalue problem for domains $Ω$ contained in the half-cylinder $Σ=ω\times (0, +\infty)$, based on a bounded regular domain $ω\subset \mathbb{R}^{N-1}$. It is easy to see that in any bounded cylinder $Ω_{t}=ω\times (0, t)$, $t > 0$, the eigenvalue problem admits a one-dimensional positive eigenfunction which satisfies the overdetermined boundary conditions. The aim of the paper is to construct other domains $Ω\subset Σ$ for which there exists a positive eigenfunction that is a solution of the overdetermined problem. This is achieved by showing that branches of such domains bifurcate from the ``trivial’’ domains $Ω_{t_j}$ at the values $t_{j} = \fracπ{2\sqrt{σ_j}}$ where $σ_j$ ($j\geq 1$) is a simple Neumann eigenvalue of the Laplace operator on $ω\subset \mathbb{R}^{N-1}$. The solutions can be reflected with respect to $ω$ to generate nontrivial solutions in a cylinder.

15. A unified proof of sharp bounds for the Jacobi heat kernel with trace and estimates of multiplicative constants

关于带迹的雅可比热核尖锐界的一个统一证明及乘性常数估计

作者: Adam Nowak, Paweł Plewa, Tomasz Z. Szarek

类别: math.AP, math.CA

发布日期: 2025-12-18

arXiv ID: 2512.16306v1

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

摘要: We give a unified and optimized proof of the sharp bounds for the Jacobi heat kernel, which were obtained gradually in several papers in recent years. We lay particular emphasis on tracing and estimating all constants appearing throughout the entire reasoning. This allows us to quantitatively control the multiplicative constants in the Jacobi heat kernel bounds in terms of the parameters involved. Consequently, analogous control extends to a number of interrelated heat kernels. In particular, we obtain quantitative control in terms of the associated dimension for the spherical heat kernel and for all other heat kernels on compact rank one symmetric spaces.

16. Nekhoroshev type stability for non-local semilinear Schrödinger equations

非局部半线性薛定谔方程的Nekhoroshev型稳定性

作者: Bingqi Yu, Li Yong

类别: math.AP, math.DS

发布日期: 2025-12-18

arXiv ID: 2512.16299v1

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

摘要: This paper investigates Nekhoroshev-type stability for solutions of ultra-differentiable regularity in Schrödinger equations with non-local nonlinear terms, employing the method of rational normal forms. We establish the first rigorous results for logarithmic ultra-differentiable regularity in infinite-dimensional Hamiltonian systems without external parameters. Under Gevrey class regularity assumptions, we achieve the stability times matching Bourgain’s conjectured optimal stability time in \cite{B04}. Furthermore, we introduce a novel global vector field norm adapted to the rational normal form framework. This norm eliminate the need for degree tracking during the iteration process, thereby enabling a unified treatment of nonlinear terms.

17. On the existence of full dimensional KAM tori for 1D periodic nonlinear Schrödinger equation

关于一维周期非线性薛定谔方程的全维KAM环面的存在性

作者: Yuan Wu

类别: math.AP, math.DS

发布日期: 2025-12-18

arXiv ID: 2512.16283v1

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

摘要: In this paper, we will prove the existence of full dimensional tori for 1-dimensional nonlinear Schrödinger equation \begin{eqnarray}\label{maineq0} \mathbf{i}u_{t}-u_{xx}+V*u+εf(x)|u|^{4}u=0,\ x\in\mathbb{T}=\mathbb{R}/2π\mathbb{Z}, \end{eqnarray} with boundary conditions, where $V*$ is the Fourier multiplier, and $f(x)$ is Gevrey smooth. Here the radius of the invariant tori satisfies a slower decay, i.e.

$$ I_n\sim e^{-2\ln^σ|n|}, \mbox{as}\ n\rightarrow\infty, $$

for any $ σ> 2, $ which extends results of Bourgain \cite{BJFA2005} and Cong \cite{cong2024} to the case that the nonlinear perturbation depends explicitly on the space variable $x$.

18. Numerical reconstruction of Schrödinger equations with quadratic nonlinearities

具有二次非线性项的薛定谔方程的数值重构

作者: Khaoula El Maddah, Matti Lassas, Teemu Tyni

类别: math.NA, math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16269v1

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

摘要: We introduce a numerical framework for reconstructing the potential in two dimensional semilinear elliptic PDEs with power type nonlinearities from the nonlinear Dirichlet to Neumann map. By applying higher order linearization method, we compute the Fourier data of the unknown potential and then invert it to recover $q$. Numerical experiments show accurate reconstructions for both smooth and discontinuous test cases.

19. A Poisson Formula for the Wave Propagator on Schwarzschild-de Sitter Backgrounds

黑洞-德西特背景上波传播子的泊松公式

作者: Izak Oltman, Ben Pineau

类别: math.AP

发布日期: 2025-12-18

arXiv ID: 2512.16054v1

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

摘要: This paper proposes a Poisson formula for the wave propagator of the Schwarzschild–de Sitter (SdS) metric. That is done by proving a Poisson formula relating wave propagators and scattering resonances for a class of non-compactly supported potentials on the real line. That class includes the Regge–Wheeler potentials obtained from separation of variables for SdS. The novelty lies in allowing non-compact supports – all exact Poisson formulae of Lax–Phillips, Melrose, and other authors required compactness of the support of the perturbation.

20. Information theory and discriminative sampling for model discovery

信息论与用于模型发现的判别性采样

作者: Yuxuan Bao, J. Nathan Kutz

类别: cs.IT, cs.LG, math.AP, math.DS

发布日期: 2025-12-17

arXiv ID: 2512.16000v1

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

摘要: Fisher information and Shannon entropy are fundamental tools for understanding and analyzing dynamical systems from complementary perspectives. They can characterize unknown parameters by quantifying the information contained in variables, or measure how different initial trajectories or temporal segments of a trajectory contribute to learning or inferring system dynamics. In this work, we leverage the Fisher Information Matrix (FIM) within the data-driven framework of {\em sparse identification of nonlinear dynamics} (SINDy). We visualize information patterns in chaotic and non-chaotic systems for both single trajectories and multiple initial conditions, demonstrating how information-based analysis can improve sampling efficiency and enhance model performance by prioritizing more informative data. The benefits of statistical bagging are further elucidated through spectral analysis of the FIM. We also illustrate how Fisher information and entropy metrics can promote data efficiency in three scenarios: when only a single trajectory is available, when a tunable control parameter exists, and when multiple trajectories can be freely initialized. As data-driven model discovery continues to gain prominence, principled sampling strategies guided by quantifiable information metrics offer a powerful approach for improving learning efficiency and reducing data requirements.

21. A Fourier analysis for $(θ,T)$-periodic functions and applications

关于 $(θ,T)$-周期函数的傅里叶分析及其应用

作者: André Pedroso Kowacs, Marielle Aparecida Silva

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15974v1

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

摘要: We develop a Fourier analysis for a generalization of the class of periodic functions, often referred to as $(θ, T)$-periodic functions, and prove several properties and inequalities related to the Fourier transform, including a type of Poincaré inequality, which extend the periodic case. As an application, we employ this analysis to show that a continuous linear operator acting on smooth $(θ, T)$-periodic functions is globally hypoelliptic/solvable if and only if the corresponding operator which acts on periodic functions is globally hypoelliptic/solvable, and characterize the global hypoellipticity/solvability of a class of first order differential operators acting on the set of smooth $(θ, T)$-periodic functions.

22. Fundamental Properties and Embedding Results in a Novel $(Φ_x, ψ)$-Fractional Musielak Space with an Application to Nonlocal BVP

一类新型 $(Φ_x, ψ)$-分数阶Musielak空间的基本性质与嵌入结果及其在非局部边值问题中的应用

作者: Ayoub Kasmi, El Houssine Azroul, Mohammed Shimi

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15972v1

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

摘要: In this paper, we introduce and study a novel class of generalized $(Φ_x,ψ)$-fractional Musielak spaces $\mathcal{K}_{Φ_x}^{α, β, ψ}$, which extends classical fractional spaces and offers the flexibility to model heterogeneous and nonlinear phenomena with memory and nonlocal effects. A detailed and rigorous analysis of their functional structure is carried out. Several new properties and embedding results are established, highlighting the originality of the proposed framework and its relevance to nonlocal BVPs. To illustrate the significance of this functional setting, we prove the existence of nontrivial solutions to a nonlinear fractional differential problem under an Ambrosetti–Rabinowitz type condition, using the mountain pass theorem. Our results provide new perspectives for the analysis of nonlocal and nonhomogeneous equations in variable-exponent and Musielak-Orlicz settings.

23. Fokas-type closed-form solution formulae for Sobolev-type equations with time-dependent coefficients

具有时变系数的Sobolev型方程的Fokas型闭式解公式

作者: Andreas Chatziafratis

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

发布日期: 2025-12-17

arXiv ID: 2512.15937v1

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

摘要: We analytically derive novel explicit integral representations for the solution of nonhomogeneous initial-boundary-value problems for a large category of evolution partial differential equations of Sobolev-Galpern type with generic temporally variable coefficients, satisfying suitable mild conditions, and with arbitrary data in classical function spaces. This work is based on the careful implementation of the pioneering Fokas unified transform methodology alongside its recently-proposed extension for solving a class of linear evolution equations with dispersion relation of specific polynomial type and time-dependent coefficients. We herein effectively extend those techniques to a special collection of evolution equations with time-dependent coefficients and mixed spatiotemporal derivatives, which induce rational dispersion relations. The new approach is exhibited in detail through illustrative generation of closed-form solutions for a multitude of such equations (such as Milne-Taylor-Barenblatt-Coleman-Ting-Chen-type, Benjamin-Bona-Mahony-type, as well as numerous higher-order variants) posed on the half-line. Challenging technical difficulties of complex-analytic and of algebraic flavour naturally emerge due the presence of mixed-derivative terms, and these are appropriately resolved in each case. The new formulas are of utility in subsequent investigation of qualitative properties and analysis of nonlinear counterparts too. Further extensions, generalizations, rigorous aspects and implementations to other types of problems as well will soon be reported in forthcoming publications.

24. Variational solutions of the Dirichlet problem, Lebesgue’s cusp and non-local properties

Dirichlet问题的变分解、Lebesgue尖点与非局部性质

作者: Wolfgang Arendt, Daniel Daners, Manfred Sauter

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15682v1

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

摘要: A recent result from [AtES24] allows one to define variational solutions of the Dirichlet problem for general continuous boundary data. We establish basic properties of this notion of solution and show that it coincides with the Perron solution. Variational solutions can elegantly be characterised in terms of the given boundary function when the variational solution has finite energy. However, it is impossible to decide in terms of the regularity of the given boundary function when a classical solution exists. We demonstrate this by analysing Lebesgue’s cusp, and more precisely Lebesgue’s domain which is associated with the potential of a thin rod with mass density going to zero at one end. We also show that the non-continuity of the Perron solution at a singular point is a generic and non-local property.