生成时间: 2025年11月14日 18:44:20

重点关注论文数量: 4

了解领域论文数量: 2

不相关论文数量: 21

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

1. Instantaneous Type~I blow-up and non-uniqueness of smooth solutions of the Navier–Stokes equations

作者: Alexey Cheskidov, Mimi Dai, Stan Palasek

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09556v1

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

详细分析

中文标题: Navier–Stokes方程光滑解的瞬时I型爆破与非唯一性

研究对象和背景

本文研究的是不可压缩Navier–Stokes方程（Navier–Stokes equations, NSE）在环面 $\mathbb{T}^d$ 上的初值问题，其中 $d \geq 2$。该方程组描述了粘性流体的速度场 $u(x,t)$ 和压力场 $p(x,t)$ 的演化过程，其数学形式如下：

$$ \begin{cases} \partial_t u - \nu \Delta u + \mathrm{div}(u \otimes u) + \nabla p = 0, \\ \mathrm{div}\, u = 0, \end{cases} \quad \text{on } \mathbb{T}^d \times [0,T], $$

初始条件为 $u(0) = u_0$，且 $u_0$ 是一个给定的光滑、无散度（divergence-free）向量场。参数 $\nu > 0$ 表示流体的粘性系数，在不失一般性的前提下，作者通常将其归一化为 $\nu = 1$。

这里的空间域是 $d$ 维环面 $\mathbb{T}^d = (\mathbb{R}/2\pi\mathbb{Z})^d$，即周期边界条件下的有界区域。这种设定广泛用于数学流体力学中，以避免无穷远边界带来的技术复杂性，同时保留非线性相互作用的核心特征。

所谓“经典解”（classical solution），是指在时空 $(x,t) \in \mathbb{T}^d \times [0,T]$ 上无限次可微（smooth）的解。Leray 在1934年的开创性工作 [Ler34] 中证明：对于任意光滑初始数据，总存在唯一的局部经典解；而对于有限能量初始数据，则存在全局弱解（称为Leray–Hopf解），但其正则性和唯一性在三维情形下仍是著名的千禧年难题之一。

长期以来，人们关注的一个核心问题是：是否可能从光滑初始数据出发，在有限时间内发生爆破（blow-up），即某个范数趋于无穷。例如，若

$$ \lim_{t \to T_*^-} \|u(t)\|_{L^\infty} = \infty, $$

则称解在时间 $T_*$ 发生有限时间爆破。当 $d=2$ 时，已知不会发生此类爆破；而当 $d \geq 3$ 时，这一问题至今未解。

本文引入并严格构造了一种全新的爆破现象——瞬时I型爆破（instantaneous Type I blow-up）。所谓“瞬时”，指的是解在某个时刻 $T_*$ 从右侧突然爆破，而在此前和此后的每一时刻都保持空间上的无限光滑性。更精确地说，解满足以下性质：

  • 在 $\mathbb{T}^d \times ([0,T] \setminus \{T_*\})$ 上是经典解；
  • 在 $t = T_*$ 时仍属于 $C^\infty(\mathbb{T}^d)$，即无空间奇异性；
  • 但从右侧逼近 $T_*$ 时，$L^\infty$ 范数发散：
    $$ \lim_{t \to T_*^+} \|u(t)\|_{L^\infty} = \infty. $$

这与传统意义上的有限时间爆破形成鲜明对比：传统的爆破是从左侧逼近 $T_*$ 时范数趋于无穷，意味着经典解无法延拓过 $T_*$；而这里的爆破是“向后看”是光滑的，“向前看”却立即变得无界。

进一步地，这种瞬时爆破并非孤立现象，而是引发了解的非唯一性。具体而言，由于解在 $t = T_*$ 仍然是光滑的，可以以此为初值应用标准局部适定性理论，得到一条始终有界的经典解 $U(t)$，它与原解 $u(t)$ 在 $[0,T_*]$ 上一致，但在 $(T_*,T]$ 上不同。因此，同一个光滑初始数据 $u_0$ 对应多个不同的弱解，从而违反了Cauchy问题的唯一性。

尤为关键的是，这些新构造的解不仅在 $L^\infty$ 意义下爆破速率达到临界尺度（Type I），而且满足一系列近乎最优的正则性条件，处于现有唯一性准则的临界边缘。例如：

  • 满足对数弱化的Ladyzhenskaya–Prodi–Serrin条件：
    $$ \int_0^T \frac{\|u(t)\|_{L^\infty}^2}{1 + (\log\log(e + \|u(t)\|_{L^\infty}))^c} dt < \infty; $$
  • 满足Beale–Kato–Majda型条件（带对数修正）；
  • 属于 $L^2([0,T]; L^p)$ 对所有 $p < \infty$ 成立；
  • 当 $d \geq 3$ 时，还满足 $u \in L^\infty_t \dot{W}^{-1,\infty}_x$，这是一个关键控制空间。

此外，解的时间连续性被控制在 $BMO^{-1}$ 的弱-*拓扑中，说明其虽爆破但仍具一定“物理合理性”。

论文强调，这种瞬时爆破机制本质上是非线性的，不能由线性热方程产生。即使热方程解可能出现类似 $1/\sqrt{t - T_*}$ 的增长，但若初始数据光滑，则不可能出现从右侧爆破的现象。真正的驱动力来自于非局部频率间的非线性相互作用，特别是逆向能量级联（inverse energy cascade）——即能量从极高波数（无限高频）自发流向低频模式的过程。

值得注意的是，这是首个明确构造出由经典Navier–Stokes流实现的完整逆向能量级联的工作，且该级联跨越无限多个频率尺度，起始于“无穷大波数”，终止于宏观尺度。这一机制打破了传统关于能量只能从前向后传递的认知，并成为实现非唯一性的核心引擎。

综上所述，本文的研究对象不仅是Navier–Stokes方程本身，更是围绕其解的长时间行为、正则性边界与决定论失效的可能性展开的深刻探索。通过构造一种前所未见的“瞬时爆破+非唯一性”现象，作者揭示了即使在解处处光滑的前提下，经典流体力学也可能丧失预测能力，挑战了我们对偏微分方程确定性本质的理解。

主要定理或主要结果

定理 1.1（瞬时I型爆破的存在性）

设 $d \geq 2$，$u_0 \in C^\infty(\mathbb{T}^d)$ 且 $\mathrm{div}\, u_0 = 0$。则存在 $T > 0$，使得对任意 $T_* \in [0,T)$，存在一个弱解 $u(t)$ 定义在 $\mathbb{T}^d \times [0,T]$ 上，满足：

  1. 对 $t \in [0,T_*]$，$u$ 是方程的经典解，且 $u(0) = u_0$；
  2. 对 $t \in (T_*, T]$，$u$ 是经典解，且存在常数 $c > 0$ 及趋于 $T_*^+$ 的时间序列 $t_n \to T_*^+$，使得
    $$ \|u(t_n)\|_{L^\infty} \geq \frac{c}{\sqrt{t_n - T_*}}, \quad \|\nabla \wedge u(t_n)\|_{L^\infty} \geq \frac{c}{t_n - T_*}; $$
  3. 解 $u(t)$ 关于时间弱-*连续，取值于 $BMO^{-1}$；
  4. 爆破为Type I型：存在常数 $C > 0$，使得对所有 $t \in (T_*, T]$，
    $$ \|u(t)\|_{L^\infty} \leq \frac{C}{\sqrt{t - T_*}}, \quad \|\nabla u(t)\|_{L^\infty} \leq \frac{C}{t - T_*}; $$
    特别地，当 $d \geq 3$ 时，还有 $u \in L^\infty_t \dot{W}^{-1,\infty}_x$；
  5. 解满足对数弱化的正则性条件：
    $$ \int_0^T \frac{\|u(t)\|_{L^\infty}^2}{1 + (\log\log(e + \|u(t)\|_{L^\infty}))^c} dt < \infty, $$
    $$ \int_0^T \frac{\|\nabla u(t)\|_{L^\infty}}{1 + (\log\log(e + \|\nabla u(t)\|_{L^\infty}))^c} dt < \infty, $$
    并且 $u \in L^2([0,T]; L^p)$ 对所有 $p < \infty$ 成立。

定理 1.6（光滑解的非唯一性）

设 $U$ 是任意定义在 $\mathbb{T}^d \times [0,T]$ 上的经典解（$d \geq 2$, $T \in (0,1]$）。则对任意 $T_* \in [0,T)$，存在一族弱解 $(u^{(\sigma)})_{\sigma \in [0,1]}$，定义在 $\mathbb{R}^d \times [0,T]$ 上，满足：

  • $u^{(0)} \equiv U$；
  • 对所有 $\sigma \in (0,1]$，有
    $$ u^{(\sigma)}\big|_{[0,T_*]} \equiv U, \quad u^{(\sigma)}\big|_{(T_*,T]} \not\equiv U; $$
  • 每个 $u^{(\sigma)}(t) \in C^\infty(\mathbb{R}^d)$ 对所有 $t \in [0,T]$ 成立；
  • 正则性包括：
    $$ u^{(\sigma)} \in L^{2,\infty}_t L^\infty_x \cap L^2_{t,\mathrm{loc}} L^p_x \quad \forall p < \infty, $$
    且当 $d \geq 3$ 时，$u^{(\sigma)} \in L^\infty_t \dot{W}^{-1,\infty}_x$；
  • 存在一个趋于零的离散序列 $\sigma_n \to 0^+$，使得 $u^{(\sigma_n)}$ 是 $2\pi$-周期的，并在弱-* $BMO^{-1}(\mathbb{T}^d)$ 中收敛到 $U$。

这意味着：从同一光滑初始数据出发，可以构造出一个连续统的彼此不同的光滑弱解，它们在 $T_*$ 之前完全重合，之后分道扬镳，从而彻底破坏了唯一性。

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

a) 核心分析学工具

本文的构造建立在高度精细的频率空间分析与渐近展开方法之上，结合了多种现代PDE技巧。

首先，作者采用仿射框架（affine framework）构建解的主部（principal part），使用一系列集中在特定频率尺度 $N_k$ 上的“建筑块”（building blocks）。这些块函数具有良好的几何结构（如螺旋状或剪切流形态），并通过精心设计的相位函数（velocity potential）来控制其传播与叠加方式。

为了处理非线性项 $u \cdot \nabla u$ 所产生的误差，作者发展了一套逐阶误差估计体系。关键在于利用振荡积分和双线性估计来控制高频交互项。特别地，他们运用了类似于para-differential calculus的思想，将速度场分解为低频与高频部分，分离出主导项与可积误差。

其次，整个构造依赖于一个不动点论证（fixed point argument），用于修正主解使其精确满足Navier–Stokes方程。为此，作者考虑线性化算子（Stokes算子）生成的半群，并建立相应的半群估计（semigroup estimates）在 $BMO^{-1}$ 和 $\dot{W}^{-1,\infty}$ 等临界空间中的表现。这类估计确保了扰动项可以在足够小的邻域内通过压缩映射原理求解。

另一个核心技术是频率级联控制。作者设计了一个递归结构，使能量能够从极高波数 $N_k \to \infty$ 逐步向下传递至低频。这一过程依赖于非线性项中三波共振机制的有效激发，尤其是在反向级联方向上保持能量守恒的同时实现放大效应。为此，他们引入了自相似缩放律与多尺度迭代方案，其中频率尺度呈双重指数增长（$N_k \sim \exp(\exp(k))$），以保证各级之间的去耦合。

此外，作者充分利用了弱紧性原理（weak compactness）来提取极限解。由于构造过程中涉及无限多个频率模态的叠加，直接验证方程满足较为困难。因此，他们先构造近似解序列，再利用 $L^\infty_t BMO^{-1}_x$ 中的弱-*紧性取极限，并借助弱解的稳定性（stability of weak solutions under weak convergence）完成最终解的验证。

最后，为了刻画解的时间不连续性而又不失物理意义，作者采用了Koch–Tataru路径空间 $X_{[T_*,T]}$，即满足某种积分型控制的空间，例如：

$$ \sup_{t \in (T_*,T]} \sqrt{t - T_*} \|u(t)\|_{L^\infty} < \infty. $$

这类空间天然适配Type I爆破行为，并允许使用标准的局部适定性理论进行拼接。

b) 技术性难题的克服

最根本的技术挑战是如何在保持解处处空间光滑的前提下，实现从右侧的$L^\infty$爆破。这要求所有奇异性仅出现在时间变量中，而非空间导数爆炸。

为此，作者提出了一种创新机制：在$t = T_*$时刻，通过非线性相互作用，瞬间从无限高波数注入能量，启动一个逆向能量级联。这个过程不是来自初始数据，而是源于方程自身的非线性反馈机制。具体来说，他们在$t = T_*$附近设计一组高频振荡模态，其振幅按$\sim 1/\sqrt{t - T_*}$增长，波数趋于无穷，从而在$L^\infty$中产生爆破，但由于快速振荡，其空间导数可通过相消干涉保持可控。

为控制非线性项带来的误差，作者采用几何相位调制技术，使得交叉项在积分后高度震荡而衰减。同时，利用对称性破缺（symmetry breaking）机制，让不同频率组件之间产生非平凡耦合，驱动能量向下转移。

此外，为了避免粘性耗散过早抑制高频活动，构造必须在极短时间内完成能量传递，使得非线性效应主导于扩散效应。这通过选择适当的频率增长率与时间尺度匹配来实现。

最终，通过将主解与一个由Stokes半群控制的小修正项相结合，作者完成了精确解的构造，并证明其满足所需的所有正则性与爆破性质。

与之前工作的比较

a) 突破与改进

本文的贡献是革命性的，体现在以下几个方面：

第一，首次构造出从光滑初值出发的瞬时I型爆破。以往关于Navier–Stokes爆破的研究几乎全部集中在“有限时间爆破”（finite-time blow-up），即从左侧逼近$T_*$时范数发散。而本文展示的是“从右侧爆破”，且解在爆破时刻仍为$C^\infty$，这是前所未有的现象。相比之下，Jia–Šverák提出的自相似解虽然也表现出Type I行为，但其初值本身就是奇异的（$(-1)$-齐次），不属于光滑范畴。

第二，实现了光滑解的非唯一性。此前的非唯一性构造主要基于凸积分（convex integration）方法，如Buckmaster–Vicol的工作[BV19]，但所得解仅为$H^\epsilon$正则，远低于经典解水平。本文的解在整个时空上都是空间光滑的，仅在时间上失去连续性，因此突破了“非唯一性必伴随低正则性”的传统认知。

第三，揭示了逆向能量级联可在经典NSE中实现。尽管在湍流理论中逆级联广为人知（如2D Euler），但在粘性Navier–Stokes中一直缺乏严格构造。本文首次证明存在一个经典流动，能将能量从无限高频传至低频，且全过程符合方程动力学。这与Palasek在Obukhov模型中的工作[Pal24]形成呼应，但后者仅为简化模型，本文则是真实PDE的结果。

第四，构造位于唯一性准则的尖锐边界。例如，解满足对数弱化的$L^2_t L^\infty_x$条件，接近Koch–Tataru空间的极限。这意味着现有的正则性标准（如Ladyzhenskaya–Prodi–Serrin或Beale–Kato–Majda）即使稍作放松，也无法保证唯一性，从而说明这些准则已是最优。

第五，适用范围广泛：结果对所有维度$d \geq 2$成立，不限于三维；初始数据任意光滑，无需特殊对称性或小性假设。

综上，本文不仅提供了新的反例，更提出了一个新的物理图像：决定论的崩溃可以发生在完全光滑的背景下，仅由非线性引起的瞬时能量注入触发。这一发现或将重塑我们对流体方程适定性理论的基本理解。

2. An existence theory for solitary waves on a ferrofluid jet

一种关于铁磁流体射流上孤立波存在性的理论

作者: Mark D. Groves, Dag Nilsson, Leon Schütz

类别: math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08799v1

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

详细分析

中文标题: 一种关于铁磁流体射流上孤立波存在性的理论

研究对象和背景

本文研究的是轴对称孤立波（axisymmetric solitary waves）在铁磁流体射流表面上的存在性问题。具体地，考虑一个无限长、不可压缩、无粘、无旋的铁磁流体柱状射流，其包围一根通有恒定电流的金属导线。该系统处于静磁场中，磁场由导线中的电流产生，方向为方位角方向（azimuthal），即围绕 $ z $-轴旋转对称。

流体区域定义为：

$$ S_1 = \{ (r, \theta, z) \in \mathbb{R}^3 \mid 0 < r < R + \eta(\theta, z, t) \}, $$

其中 $ r $ 是径向坐标，$ \theta $ 是方位角，$ z $ 是轴向坐标，$ R > 0 $ 是未扰动时射流的半径，而 $ \eta(\theta, z, t) $ 表示自由表面相对于圆柱面 $ r = R $ 的扰动。外部区域 $ S_2 = \{ r > R + \eta \} $ 被假定为真空。

由于物理系统的旋转对称性，作者专注于轴对称波，即所有变量不依赖于 $ \theta $，因此 $ \eta = \eta(z, t) $，速度势 $ \phi = \phi(r, z, t) $，磁势也具有相应对称性。

铁磁流体满足非线性磁化规律，其磁感应强度与磁场的关系为：

$$ \mathbf{B}_1 = \mu_0 (\mathbf{H}_1 + \mathbf{M}_1(\mathbf{H}_1)), \quad \mathbf{B}_2 = \mu_0 \mathbf{H}_2, $$

其中 $ \mu_0 $ 是真空磁导率，磁化强度 $ \mathbf{M}_1 $ 满足：

$$ \mathbf{M}_1(\mathbf{H}_1) = m_1(|\mathbf{H}_1|) \frac{\mathbf{H}_1}{|\mathbf{H}_1|}, $$

且 $ m_1(s) \geq 0 $ 是给定的非负函数。

引入磁势 $ \psi_1, \psi_2 $ 和速度势 $ \phi $，使得：

$$ \mathbf{H}_1 = -\nabla \psi_1, \quad \mathbf{H}_2 = -\nabla \psi_2, \quad \mathbf{v} = \nabla \phi. $$

控制方程如下：

在流体内部 $ 0 < r < R + \eta(z,t) $：

$$ \nabla \cdot (\mu(|\nabla \psi_1|) \nabla \psi_1) = 0, $$

$$ \Delta \phi = 0; $$

在外部真空区域 $ r > R + \eta(z,t) $：

$$ \Delta \psi_2 = 0. $$

边界条件位于自由表面 $ r = R + \eta(z,t) $ 上：

  • 连续性条件（磁场与磁感应场连续）：

    $$ \psi_2 - \psi_1 = 0, $$

    $$ \psi_{2n} - \mu(|\nabla \psi_1|) \psi_{1n} = 0; $$

  • 流体动力学边界条件（运动学）：

    $$ -\eta_t + \phi_r - \phi_z \eta_z = 0; $$

  • 动力学平衡条件（力平衡）：

    $$ \phi_t + \frac{1}{2} |\nabla \phi|^2 - \mu_0 \nu(|\nabla \psi_1|) + 2\sigma \kappa - \frac{1}{2} \mu_0 (\mu(|\nabla \psi_1|) - 1)^2 = c_0, $$

    其中 $ \nu(s) $ 是磁能密度函数，$ \sigma $ 是表面张力系数，$ \kappa $ 是平均曲率，$ c_0 $ 是积分常数。

通过量纲归一化，并考虑行波解形式：$ \phi(r,z,t) = \phi(r, z - ct) $，$ \eta(z,t) = \eta(z - ct) $，系统被简化为仅依赖于单变量 $ z $ 的稳态问题。

最终得到的核心方程是通过Zakharov-Craig-Sulem型非局部公式推导出的一个关于自由表面扰动 $ \eta(z) $ 的非局部微分方程：

$$ K(\eta) - c^2 L(\eta) = 0, $$

其中 $ K(\eta) $ 和 $ L(\eta) $ 是由 Dirichlet-Neumann 映射构造的非线性泛函，具体形式为：

$$ K(\eta) = -\gamma \left( \nu\left( \frac{1}{1+\eta} \right) - \nu(1) \right) + \frac{\eta}{(1+\eta)\sqrt{1+\eta_z^2}} - \frac{\eta_{zz}}{(1+\eta_z^2)^{3/2}} - 1, $$

$$ L(\eta) = -\frac{1}{2}(K(\eta)\eta + \tfrac{1}{2}K(\eta)\eta^2)^2 + \frac{1}{2(1+\eta_z^2)}(\eta_z - \eta_z K(\eta)\eta - \tfrac{1}{2}\eta_z K(\eta)\eta^2)^2 + K(\eta)\eta + \tfrac{1}{2}K(\eta)\eta^2, $$

而算子 $ K(\eta) \xi = - (G(\eta)^{-1} \xi)_z $，其中 $ G(\eta) $ 是Dirichlet-Neumann算子，定义如下：

设 $ \phi $ 是以下边值问题的解：

$$ \Delta \phi = 0 \quad \text{in } 0 < r < 1 + \eta(z), $$

$$ \phi|_{r=1+\eta} = \Phi, $$

则

$$ G(\eta)\Phi := (1+\eta)(1+\eta_z^2)^{1/2} \frac{\partial \phi}{\partial n}\bigg|_{r=1+\eta} = (1+\eta)(\phi_r - \eta_z \phi_z)|_{r=1+\eta}. $$

此非局部结构是本文分析的关键难点之一。此外，问题定义在无限长柱面上，自变量 $ z \in \mathbb{R} $，目标是寻找满足 $ \eta(z) \to 0 $ 当 $ |z| \to \infty $ 的非平凡孤立波解。

值得注意的是，文中使用的函数空间并非标准 Sobolev 空间，而是特别构造的空间以适应 Fourier 分析和小参数展开的需求。特别是定义了空间：

$$ Z = \left\{ \eta \in \mathcal{S}'(\mathbb{R}) \,\middle|\, \|\eta\|_Z := \|\hat{\eta}_1\|_{L^1(\mathbb{R})} + \|\eta_2\|_{L^2(\mathbb{R})} < \infty \right\}, $$

其中 $ \eta_1 = \chi(D)\eta $，$ \eta_2 = (1 - \chi(D))\eta $，$ \chi $ 是某个特征函数，用于将 $ \eta $ 的 Fourier 变换分解为集中在低频或共振频率附近的部分和其余部分。这种分解对于后续的Lyapunov-Schmidt约化至关重要。

另一个关键点是使用了径向 Sobolev 空间来处理轴对称边界值问题。作者引入了基于 Bessel 算子的函数类 $ H^{(q)}_m((0,1)\times\mathbb{R};\mathbb{C}) $，并利用 Wirtinger 型复微分算子 $ \partial_\tau, \partial_{\bar{\tau}} $ 来刻画模态函数的行为。这使得他们能够精确估计 Green 函数并建立解算子的解析性。

总之，这是一个涉及非线性偏微分方程组、非局部算子、自由边界问题以及特殊几何对称性（轴对称）的复杂系统。其数学挑战在于如何将高维 PDE 系统降维为单一非局部 ODE，并在此基础上严格证明孤立波的存在性。

主要定理或主要结果

定理 1.1 （强表面张力情形下的 KdV 型孤立波）

假设 $ 1 < \gamma < 9 $，且波速满足 $ c^2 = c_0^2(1 - \varepsilon^2) $，其中 $ c_0^2 = \frac{1}{2}(\gamma - 1) $。则对所有足够小的 $ \varepsilon > 0 $，存在一个对称的 Korteweg-de Vries 孤立波解 $ \eta(z) $ 满足：

$$ \eta(z) = \varepsilon^2 \zeta_{\mathrm{KdV}}(\varepsilon z) + o(\varepsilon^2), \quad \text{一致于 } z \in \mathbb{R}, $$

其中

$$ \zeta_{\mathrm{KdV}}(Z) = -\frac{3}{2d_0} \mathrm{sech}^2\left( \sqrt{\frac{c_0^2}{9 - \gamma}} Z \right), $$

且

$$ d_0 = \frac{1}{2c_0^2} \left( \gamma - \frac{3}{2} \gamma \nu''(1) - \frac{3}{2} \right). $$

注：当 $ d_0 > 0 $ 时为凹陷波（depression wave），$ d_0 < 0 $ 时为隆起波（elevation wave）。

定理 1.2 （弱表面张力情形下的 NLS 型孤立波）

假设 $ \omega > 0 $，

$$ \gamma = 1 - \omega^2 + \frac{2\omega f(\omega)}{f'(\omega)}, \quad c^2 = \frac{2\omega}{f'(\omega)} = c_0^2(1 - \varepsilon^2), $$

其中 $ f(k) = \frac{|k| I_0(|k|)}{I_1(|k|)} $，$ I_\nu $ 是第 $ \nu $ 阶修正贝塞尔函数。则对所有足够小的 $ \varepsilon > 0 $，存在两个对称的非线性薛定谔方程型孤立波解满足：

$$ \eta(z) = \pm \varepsilon \zeta_{\mathrm{NLS}}(\varepsilon z) \cos(\omega z) + o(\varepsilon), \quad \text{一致于 } z \in \mathbb{R}, $$

其中

$$ \zeta_{\mathrm{NLS}}(Z) = \left( \frac{2a_2}{a_3} \right)^{1/2} \mathrm{sech}\left( \left( \frac{a_2}{a_1} \right)^{1/2} Z \right), $$

且 $ a_1, a_2, a_3 > 0 $ 是依赖于 $ \omega $ 的正常数。

这两个定理共同构成了论文的主要成果：在不同参数范围内，严格证明了铁磁流体射流上存在两类不同的小振幅轴对称孤立波——一类类似于经典水波中的 KdV 波（长波极限），另一类则是调制周期波包型的 NLS 波（短波共振情形）。这些解是非平凡、光滑、衰减于无穷远处的。

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

本文采用了一套结合非局部公式化、谱分解、Lyapunov-Schmidt 约化与隐函数定理的现代泛函分析方法，系统地将原始复杂的自由边界问题转化为可处理的有限维约化方程，并最终借助已知模型方程的非退化解构造近似解。

a) 核心分析工具

  1. Zakharov-Craig-Sulem 非局部公式化

    作者借鉴经典水波理论中的 ZCS 方法，将原三维势流问题通过引入 Dirichlet-Neumann 算子 $ G(\eta) $，降维为仅关于自由表面高度 $ \eta(z) $ 的非局部方程。这种方法避免了直接求解整个域内的 Laplace 方程，而是将其“编码”进一个边界算子中。这是处理自由边界问题的标准利器。

  2. Fourier 分析与色散关系研究

    在弱非线性近似下，线性化方程给出色散关系：

    $$ c^2 = \frac{\gamma - 1 + k^2}{f(k)}, \quad f(k) = \frac{|k| I_0(|k|)}{I_1(|k|)}. $$

    通过对 $ f(k) $ 的单调性分析（见附录 A），发现当 $ 1 < \gamma \leq 9 $ 时最小相速度出现在 $ k = 0 $；而当 $ \gamma > 9 $ 时最小值出现在某个 $ k = \omega > 0 $。这一分岔现象决定了后续应使用何种渐近模型（KdV 或 NLS）。

  3. Lyapunov-Schmidt 分解

    这是全文最核心的技术之一。令 $ g(k) = \gamma - 1 + k^2 - c_0^2 f(k) $，注意到 $ g(\pm \omega) = g'(\pm \omega) = 0 $，$ g''(\omega) > 0 $，说明线性算子在 $ k = \pm \omega $ 处有双重零模。于是将解 $ \eta $ 分解为：

    $$ \eta = \eta_1 + \eta_2, \quad \hat{\eta}_1 = \chi(D) \hat{\eta}, \quad \hat{\eta}_2 = (1 - \chi(D)) \hat{\eta}, $$

    其中 $ \chi $ 是支撑在 $ (-\omega - \delta, -\omega + \delta) \cup (\omega - \delta, \omega + \delta) $ 上的特征函数（若 $ \omega = 0 $ 则为 $ (-\delta, \delta) $）。这样 $ \eta_1 $ 包含“慢变调制”成分，$ \eta_2 $ 代表高频修正。

    将原方程 $ K(\eta) - c^2 L(\eta) = 0 $ 投影到两个子空间：

    $$ \chi(D)[K(\eta_1 + \eta_2) - c^2 L(\eta_1 + \eta_2)] = 0, $$

    $$ (1 - \chi(D))[K(\eta_1 + \eta_2) - c^2 L(\eta_1 + \eta_2)] = 0. $$

    利用压缩映射原理，可在 $ \varepsilon $ 小时从第二个方程解出 $ \eta_2 = \eta_2(\eta_1) $，代入第一个方程得约化方程：

    $$ \chi(D)\left[ K(\eta_1 + \eta_2(\eta_1)) - c^2 L(\eta_1 + \eta_2(\eta_1)) \right] = 0. $$

  4. 尺度变换与极限模型逼近

    引入小参数 $ \varepsilon $，进行伸缩变换：

    • 对于 $ 1 < \gamma < 9 $：令 $ \eta_1(z) = \varepsilon^2 \zeta(\varepsilon z) $
    • 对于 $ \gamma > 9 $：令 $ \eta_1(z) = \frac{1}{2} \varepsilon [\zeta(\varepsilon z)e^{i\omega z} + \overline{\zeta(\varepsilon z)} e^{-i\omega z}] $

    经过细致的 Taylor 展开与符号计算，约化方程变为：

    $$ \varepsilon^{-2} g(\varepsilon D)\zeta + 2c_0^2 \zeta + 2c_0^2 d_0 \chi_0(\varepsilon D)\zeta^2 + O(\varepsilon^{1/2}) = 0 \quad (\text{KdV 类}), $$

    $$ \varepsilon^{-2} g(\omega + \varepsilon D)\zeta + a_2 \zeta - a_3 \chi_0(\varepsilon D)(|\zeta|^2 \zeta) + O(\varepsilon^{1/2}) = 0 \quad (\text{NLS 类}), $$

    其中 $ \chi_0 $ 支撑在 $ (-\delta, \delta) $。这些被称为“全色散版本”的 KdV/NLS 方程，保留了原始系统的完整线性结构。

  5. 隐函数定理的应用

    已知标准 KdV 和 NLS 方程的孤立波解 $ \zeta_{\mathrm{KdV}}, \zeta_{\mathrm{NLS}} $ 是非退化的（即其线性化算子的核仅为平移生成的一维空间）。结合反射对称性（$ \eta(z) \mapsto \eta(-z) $），限制在偶函数空间后，可消除奇函数方向的零模。于是可用 Banach 空间中的隐函数定理，将真实解视为这些显式解的小扰动，从而完成存在性证明。

b) 技术难点与解决技巧

  1. Dirichlet-Neumann 算子的解析性

    传统方法难以直接处理非平直边界上的 DN 算子正则性。作者创新地引入“压平变换”（flattening transformation）：

    $$ r' = \frac{r}{1 + \eta(z)}, \quad u(r', z) = \phi(r, z), $$

    将变动区域映射到固定带状域 $ (0,1) \times \mathbb{R} $，并将原问题转化为带有非线性项的椭圆边值问题：

    $$ D_1 D_0 u + u_{zz} = D_1 F_1(\eta, u) + \partial_z F_2(\eta, u), $$

    $$ D_0 u|_{r=1} = F_1(\eta, u) + \xi_z. $$

    此处 $ D_j $ 为 Bessel 算子，$ F_1, F_2 $ 由 $ \eta $ 和 $ u $ 构成的非线性项。

  2. Green 函数估计与新型函数空间设计

    解算子表达为积分形式：

    $$ u = \mathcal{F}^{-1} \left[ ik \int_0^1 G(r, \tilde{r}) \hat{F}_2 - D_{\tilde{r}} G(r, \tilde{r}) \hat{F}_1 \, d\tilde{r} - ik G(r,1) \hat{\xi} \right], $$

    其中 Green 函数 $ G $ 由修正贝塞尔函数 $ I_\nu, K_\nu $ 构造。通过对 $ I_\nu, K_\nu $ 的精细渐近估计，在新定义的径向 Sobolev 空间 $ H^{(q)}_m $ 中建立了统一估计：

    $$ \|S(F_1, F_2, \xi)\|_\star \lesssim \|F_1\|_{H^1_{(1)}} + \|F_2\|_{H^1_{(0)}} + \|\xi\|_{H^{3/2}}, $$

    并进一步证明 $ F_1, F_2 $ 关于 $ (\eta, u) $ 是解析映射（Lemma 2.2），从而保证 $ K(\eta) $ 的解析性。

  3. 非局部项的控制与余项估计

    在约化过程中，必须精确估计高阶非线性项和余项。作者利用 $ \eta \in Z $ 空间蕴含 $ \eta \in W^{1,\infty} $ 的事实（因 $ \hat{\eta}_1 \in L^1 \Rightarrow \eta_1 \in C_b^1 $），确保所有非线性组合均有界。同时利用 $ \varepsilon $ 小量控制误差项至 $ O(\varepsilon^{1/2}) $，足以支持隐函数定理应用。

与之前工作的比较

本文相较于以往研究做出了若干重要突破：

  1. 首次严格存在性证明

    尽管此前已有大量关于铁磁流体射流波的研究，如 Bashtovoi et al. [2] 使用经验模型方程，Bourdin et al. [5] 开展实验观测，Blyth & Parau [3], Guyenne & Parau [13] 数值模拟，Groves & Nilsson [12] 使用空间动力系统方法，但本文是首个基于第一性原理（从完整控制方程出发）、通过严格泛函分析手段证明孤立波存在的工作。

  2. 放宽模型假设：一般磁化律

    早期理论常假设磁化强度为线性或 Langevin 形式。本文允许一般的非线性磁化律 $ m_1(s) $，仅要求其光滑性，大大增强了结果的普适性。这体现在泛函 $ \nu(s) $ 的任意性上。

  3. 克服非局部性障碍

    相比于 Blyth & Parau [4] 提出的另一种非局部公式，本文采用更接近水波理论的 ZCS 框架，并成功证明了相应 Dirichlet-Neumann 算子的解析性，这是构造摄动解的基础。此前尚无对此类轴对称磁流体问题中 DN 算子正则性的系统研究。

  4. 统一处理两种波型

    许多文献只关注某一种波（如驻波或行波）。本文在一个框架下统一处理了 $ \gamma < 9 $ 和 $ \gamma > 9 $ 两种情形，分别导出 KdV 和 NLS 型孤立波，揭示了参数变化导致的分岔机制。

  5. 改进函数空间技术

    作者发展了适用于轴对称问题的新型 Sobolev 空间理论（基于 Bessel 算子模态分解），解决了传统 Cartesian 坐标系下难以处理柱面对称性的难题。这种空间设计使得 Green 函数估计和非线性项控制成为可能。

  6. 获得更强的渐近精度

    与形式渐近法不同，本文不仅给出存在性，还提供了一致渐近展开式（如 $ \eta = \varepsilon^2 \zeta_{\mathrm{KdV}}(\varepsilon z) + o(\varepsilon^2) $），且误差可控。这意味着所构造的解真正逼近理想模型方程的解。

综上所述，本文在数学严谨性、模型广度、技术深度方面均显著超越了前人工作，为铁磁流体力学中的非线性波现象提供了坚实的理论基础。

3. Well-posedness for a diffuse interface model of non-Newtonian two-phase flows

非牛顿两相流 diffuse interface 模型的适定性研究

作者: Fang Li, Duan Xingyu, Guo Zhenhua

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.08876v1

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

详细分析

中文标题: 非牛顿两相流 diffuse interface 模型的适定性研究

研究对象和背景

本文研究的是在三维有界区域中，描述两种部分可混溶、非均匀、不可压缩且具有非牛顿特性的粘性流体混合演化的 Navier-Stokes-Cahn-Hilliard (NSCH) 系统。这类系统广泛用于模拟宏观上不相容但存在微观扩散效应的两相流体动力学行为，尤其是在涉及液滴形成、合并与破裂等拓扑结构变化的过程中，传统的 sharp interface 模型难以处理界面奇异性，而 diffuse interface（或称相场）方法通过引入一个连续的浓度差函数 $\phi$ 来光滑地刻画界面过渡层，从而避免了显式追踪界面的困难。

所考虑的数学模型为如下耦合偏微分方程组：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0 & \text{in } \Omega \times (0,T), \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) - \mathrm{div}( \nu(\phi)(1+|Du|^2)^{\frac{p-2}{2}} Du ) + \nabla P = -\mathrm{div}(\nabla\phi \otimes \nabla\phi) & \text{in } \Omega \times (0,T), \\ \partial_t (\rho \phi) + \mathrm{div}(\rho u \phi) = \Delta \mu & \text{in } \Omega \times (0,T), \\ \rho \mu = -\Delta \phi + \rho \Psi'(\phi) & \text{in } \Omega \times (0,T), \end{cases} $$

其中：

  • $\rho(x,t)$ 表示混合物的总密度；
  • $u(x,t)$ 是质量平均速度；
  • $P(x,t)$ 是压力；
  • $\phi(x,t)$ 是两种流体的浓度差（order parameter），取值通常在 $[-1,1]$ 范围内表示不同相；
  • $\mu(x,t)$ 是化学势；
  • $Du = \frac{1}{2}(\nabla u + (\nabla u)^t)$ 是应变率张量；
  • $\nu(\phi)$ 是依赖于浓度的粘度系数，满足 $0 < \nu_* \leq \nu(s) \leq \nu^* < \infty$，且 $\nu \in W^{1,\infty}(\mathbb{R})$；
  • 自由能势 $\Psi(s)$ 取为 Landau 势：$\Psi(s) = \frac{1}{4}(s^2 - 1)^2$，这是典型的双阱势，对应于两个稳定相 $s = \pm 1$。

该系统的物理基础来源于热力学第二定律和能量耗散机制。整个系统的自由能通常定义为：

$$ E(\phi, u) = \int_\Omega \left( \frac{1}{2} \rho |u|^2 + \frac{1}{2} |\nabla \phi|^2 + \rho \Psi(\phi) \right) dx, $$

其时间导数沿解满足耗散关系，体现了系统向平衡态演化的过程。

边界条件为 Neumann 型：

$$ u = 0, \quad \partial_n \phi = 0, \quad \partial_n \mu = 0 \quad \text{on } \partial\Omega \times (0,T), $$

初始条件为：

$$ \rho|_{t=0} = \rho_0, \quad u|_{t=0} = u_0, \quad \phi|_{t=0} = \phi_0. $$

值得注意的是，本文的关键特征之一是允许初始密度 $\rho_0$ 包含零点（即允许真空出现），这在分析上带来了显著挑战，因为动量方程中的惯性项 $\partial_t(\rho u)$ 和对流项 $\mathrm{div}(\rho u \otimes u)$ 在 $\rho = 0$ 处失去控制，常规的能量估计不再直接适用。此外，粘度 $\nu$ 显式依赖于浓度 $\phi$，增强了系统的非线性耦合程度。

研究区域 $\Omega \subset \mathbb{R}^3$ 是一个具有 $C^2$ 边界的有界域；而在强解的存在性部分，则设定为三维环面 $\mathbb{T}^3$，便于使用周期性结构进行谱分析。函数空间方面，文中大量使用 Sobolev 空间 $W^{s,q}(\Omega)$、Lebesgue 空间 $L^q(\Omega)$，以及针对不可压流体构造的 Solenoidal 函数空间 $H_\sigma$ 和 $V_\sigma$，分别为 $C_{0,\sigma}^\infty(\Omega)$ 在 $L^2$ 和 $H^1$ 范数下的闭包。

特别地，由于密度可能退化（即 $\rho$ 可以为零），作者采用加权空间如 $\sqrt{\rho} u \in L^\infty(0,T;L^2(\Omega))$ 来控制速度场，并借助 Gagliardo-Nirenberg 不等式、椭圆正则性理论 和 单调性方法 来克服退化带来的技术障碍。

综上所述，本工作的核心在于建立一个包含密度退化、浓度相关粘度及 Landau 势的非牛顿两相流 diffuse interface 模型的全局弱解存在性和局部强解存在性，填补了此前文献中关于此类复杂情形下适定性分析的空白。

主要定理或主要结果

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

设 $\Omega \subset \mathbb{R}^3$ 为具有 $C^2$ 边界的有界区域，$T > 0$，且 $p \in (5/2, +\infty)$。假设初始数据 $(\rho_0, u_0, \phi_0, \mu_0)$ 满足：

$$ 0 \leq \rho_0 \in L^\infty(\Omega), \quad \phi_0 \in H^1(\Omega), \quad \sqrt{\rho_0} u_0 \in L^2(\Omega), \quad \nabla \mu_0 \in L^2(\Omega). $$

则初边值问题 (1.1)-(1.4) 存在一个全局弱解 $(\rho, u, \phi, \mu)$，定义在 $[0,T]$ 上，满足以下正则性：

$$ \rho \in C([0,T]; L^r(\Omega)) \cap L^\infty(\Omega \times (0,T)), \quad \partial_t \rho \in L^\infty(0,T; H^{-1}(\Omega)), $$

$$ \sqrt{\rho} u \in L^\infty(0,T; L^2(\Omega)), \quad u \in L^p(0,T; W^{1,p}(\Omega)), $$

$$ \sqrt{\rho} \phi \in L^\infty(0,T; L^2(\Omega)), \quad \phi \in L^\infty(0,T; H^2(\Omega)) \cap C([0,T]; H^1(\Omega)), $$

$$ \sqrt{\rho} \mu \in L^\infty(0,T; L^2(\Omega)), \quad \mu \in L^\infty(0,T; H^1(\Omega)) \cap L^2(0,T; H^2(\Omega)). $$

并且以下积分恒等式成立：

(i) 对所有 $\psi \in C_c^\infty(\Omega \times (0,T))$，

$$ \int_0^T \int_\Omega \rho \partial_t \psi \,dxdt + \int_0^T \int_\Omega \rho u \cdot \nabla \psi \,dxdt = 0; $$

(ii) 对所有 $w \in C^1([0,T); W^{1,p}_0(\Omega))$，

$$ -\int_0^T \int_\Omega \rho u \cdot \partial_t w \,dxdt - \int_0^T \int_\Omega \rho u \otimes u : \nabla w \,dxdt + \int_0^T \int_\Omega \nu(\phi)(1+|Du|^2)^{\frac{p-2}{2}} Du : Dw \,dxdt $$

$$ = \int_\Omega \rho_0 u_0 \cdot w(0)\,dx + \int_0^T \int_\Omega \nabla\phi \otimes \nabla\phi : \nabla w \,dxdt; $$

(iii) 对所有 $w \in C^1([0,T); H^1(\Omega))$，

$$ -\int_0^T \int_\Omega \rho \phi \partial_t w \,dxdt - \int_0^T \int_\Omega \rho u \phi \cdot \nabla w \,dxdt + \int_0^T \int_\Omega \nabla \mu \cdot \nabla w \,dxdt = \int_\Omega \rho_0 \phi_0 w(0)\,dx; $$

(iv) 几乎处处在 $\Omega \times (0,T)$ 上有：

$$ \rho \mu = -\Delta \phi + \rho \Psi'(\phi); $$

此外，$\partial_n \phi = 0$ 几乎处处在 $\partial\Omega \times (0,T)$ 上成立。

定理 1.2（局部强解的存在性）

设 $\Omega = \mathbb{T}^3$，$T > 0$，且 $p \in [5/2, 3)$。假设初始数据 $(\rho_0, u_0, \phi_0, \mu_0)$ 满足：

$$ \rho_0 \in L^\infty(\Omega), \quad 0 \leq \rho_0 \leq \rho_*, \quad u_0 \in W^{1,p}(\Omega), \quad \nabla \mu_0 \in L^2(\Omega), \quad \phi_0 \in H^1(\Omega). $$

则存在 $T_0 > 0$（仅依赖于初始数据），使得问题 (1.1)-(1.4) 在 $[0,T_0]$ 上存在一个强解 $(\rho, u, \phi, \mu)$，满足：

$$ \rho \in C([0,T_0]; L^r(\Omega)) \cap L^\infty(\Omega \times (0,T_0)), \quad \forall r \in [1,\infty), $$

$$ u \in C([0,T_0]; L^q(\Omega)) \cap L^\infty(0,T_0; W^{1,p}(\Omega)) \cap L^2(0,T_0; H^2(\Omega)) \cap H^1(0,T_0; H^1(\Omega)), \quad \forall q \in [1,6), $$

$$ \phi \in C([0,T_0]; W^{1,q}(\Omega)) \cap L^\infty(0,T_0; H^2(\Omega)) \cap L^2(0,T_0; W^{2,\infty}(\Omega)), $$

$$ \mu \in L^\infty(0,T_0; H^1(\Omega)) \cap L^2(0,T_0; H^2(\Omega)). $$

此外，压力 $P$ 满足 $P \in L^2(0,T_0; W^{1, \frac{3p}{2p-4}}(\Omega))$。

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

a) 核心分析工具

本文的主要证明策略基于 半Galerkin逼近方案（semi-Galerkin scheme） 与 单调性方法（monotonicity method） 的结合。具体而言，作者首先对原始系统进行正则化处理：将初始密度 $\rho_0$ 正则化为严格正的光滑函数 $\rho_{0\delta}$，并对速度和相场变量分别在有限维子空间中进行投影逼近。这种构造保留了动量方程的 Galerkin 投影特性，同时保持 Cahn-Hilliard 方程的结构不变，因此称为“半”Galerkin 方法。

在此框架下，通过经典的常微分方程理论可以证明每个近似系统存在唯一解。接下来的关键步骤是对这些近似解建立统一的能量估计，以支持紧性论证和极限过程。

主导的能量恒等式来自系统的时间演化律：

$$ \frac{d}{dt} \int_\Omega \left( \frac{1}{2} \rho_m |u_m|^2 + \frac{1}{2} |\nabla \phi_m|^2 + \rho_m \Psi(\phi_m) \right) dx + \int_\Omega \nu(\phi_m)(1+|Du_m|^2)^{\frac{p-2}{2}} |Du_m|^2 dx + \int_\Omega |\nabla \mu_m|^2 dx = 0. $$

这一恒等式提供了多个关键先验估计：

  • $\sqrt{\rho_m} u_m$ 在 $L^\infty(L^2)$ 中有界；
  • $u_m$ 在 $L^p(W^{1,p})$ 中有界；
  • $\phi_m$ 在 $L^\infty(H^1)$ 中有界；
  • $\mu_m$ 在 $L^2(H^1)$ 中有界。

为了处理密度可能为零的情况，作者引入了 加权范数控制，例如利用 $\sqrt{\rho} u$ 和 $\sqrt{\rho} \phi$ 来替代传统意义上的 $u$ 或 $\phi$ 的整体控制。进一步，借助 Gagliardo-Nirenberg 插值不等式 和 Poincaré 型不等式（见 Lemma 2.3），即使在密度退化的情况下也能获得速度场的 $L^q$ 控制。

另一个关键技术是 椭圆正则性引理（Lemma 2.1），它用于从化学势方程 $\rho \mu = -\Delta \phi + \rho \Psi'(\phi)$ 推出 $\phi$ 的高阶正则性。由于 $\rho$ 可能为零，不能直接应用标准椭圆估计，但通过对 $\mu$ 的控制并结合权重，仍可得到 $\phi \in L^\infty(H^2)$。

最后，在极限过程中，非线性项（如 $\nu(\phi)(1+|Du|^2)^{(p-2)/2} Du$）的收敛性依赖于 单调性方法（monotonicity method）。该方法的核心思想是：若算子 $A(u)$ 是单调的，且近似解 $u_m \rightharpoonup u$ 弱收敛，同时 $\langle A(u_m), u_m \rangle \to \langle A(u), u \rangle$，则可推出 $A(u_m) \rightharpoonup A(u)$。文中引用了 Browder-Minty 类型的结果（见 Lemma 2.5），确保粘性应力张量的弱极限正确传递。

b) 技术难点与解决技巧

  1. 密度退化问题（Vacuum case）
    当 $\rho_0$ 允许为零时，动量方程中 $\partial_t(\rho u)$ 和 $\mathrm{div}(\rho u \otimes u)$ 的定义变得敏感。作者通过引入 $\sqrt{\rho} u$ 这一加权量，在能量空间中绕过真空点的影响。同时，在测试函数选取时谨慎处理，避免除以 $\rho$。

  2. 浓度依赖粘度 $\nu(\phi)$
    粘度不再是常数或仅依赖于空间，而是随 $\phi$ 变化，导致动量方程高度非线性。但由于假设 $\nu \in W^{1,\infty}$ 且有正下界，可通过一致估计控制其影响。

  3. Landau 势的增长性
    $\Psi'(\phi) = \phi^3 - \phi$ 属于三次非线性，虽比对数势更易处理，但仍需精细估计。文中通过 Sobolev 嵌入 $H^1 \hookrightarrow L^6$ 和插值技巧控制高次项。

  4. 非牛顿流体的 $p$-Laplacian 结构
    应力项 $(1+|Du|^2)^{(p-2)/2} Du$ 在 $p \ne 2$ 时不具备线性结构。当 $p > 2$ 为剪切增稠，$p < 2$ 为剪切稀化。文中要求 $p > 5/2$ 以保证足够的可积性，使 $Du \in L^p$ 能嵌入到合适的空间。

  5. 极限过程中的紧性缺失
    由于缺乏时间导数的一致估计（尤其对 $u$），无法直接使用 Aubin-Lions 引理。为此，作者转而利用 弱收敛 + 单调性 的组合来恢复非线性项的极限。

与之前工作的比较

a) 突破与改进

本文相较于已有工作有多项实质性推进：

  1. 首次允许初始密度含真空
    此前 Giorgini-Temam [10] 等人在非均匀 NSCH 模型中均假设 $\rho_0 \geq \rho_- > 0$，以避免退化。本文完全去除此限制，适用于更广泛的物理场景（如空腔填充、气液界面）。

  2. 推广至非牛顿流体类型
    经典 NSCH 多限于牛顿流体（$p=2$）。本文考虑 power-law 型非牛顿流体，涵盖剪切稀化（$p<2$）与增稠（$p>2$）行为，更具实际意义（如聚合物溶液、泥浆流动）。

  3. 结合 Landau 势与浓度依赖粘度
    同时处理非常规粘度结构与多项式势能的情形此前鲜见成果。多数研究要么用 Flory-Huggins 对数势，要么假设常粘度。

  4. 在三维情况下建立全局弱解
    尽管 Kim 等人 [20] 研究过非牛顿 NSCH，但要求 $p \geq 3d/(d+2)$ 且未考虑密度变化。本文在 $p > 5/2$ 下实现三维全局弱解，条件更宽松。

  5. 局部强解在周期域上的构造
    利用 $\mathbb{T}^3$ 的谱性质，通过更高阶估计获得强解，为后续研究爆破准则或长时间行为奠定基础。

综上，本文不仅拓展了非牛顿两相流模型的数学适定性理论，也为含真空、变粘度、非线性流变行为的实际问题提供了坚实的分析框架。

4. On the Singular Limit in Hibler’s Sea Ice Model

关于Hibler海冰模型中的奇异极限

作者: Robert Denk, Franz Gmeineder, Matthias Hieber

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09327v1

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

详细分析

中文标题: 关于Hibler海冰模型中的奇异极限

研究对象和背景

本文研究的是Hibler海冰模型（Hibler’s sea ice model）中动量平衡方程的奇异极限（singular limit），即在不进行应力正则化的前提下，建立该模型的能量驱动解的存在性理论。该模型是描述极地海冰动力学的核心数学工具之一，广泛应用于气候模拟与地球物理流体力学。

论文聚焦于以下动量平衡方程：

$$ m\dot{u} = \mathrm{div}(\sigma) - mg\nabla H + \tau_{\mathrm{atm}} + \tau_{\mathrm{ocean}}(u) \quad \text{in } \Omega_T = (0,T)\times\Omega, $$

其中：

  • $ u: \Omega_T \to \mathbb{R}^2 $ 是水平方向上的海冰速度场；
  • $ m > 0 $ 是单位面积上的冰质量（文中假设为常数）；
  • $ \dot{u} $ 表示物质导数（忽略对流项 $ u\cdot\nabla u $ 和科里奥利力后简化为时间导数）；
  • $ \sigma $ 是粘塑性应力张量（viscoplastic stress），由Hibler提出，具有高度非线性和退化结构；
  • $ g $ 是重力加速度，$ H $ 是海表动态高度；
  • $ \tau_{\mathrm{atm}} $ 和 $ \tau_{\mathrm{ocean}}(u) $ 分别表示大气风应力和海洋拖曳力。

海洋力被建模为带截断的非线性函数：

$$ \tau_{\mathrm{ocean}}(u) = \rho_{\mathrm{ocean}} C_{\mathrm{ocean}} \eta(|U_{\mathrm{ocean}} - u|) R_{\mathrm{ocean}}(U_{\mathrm{ocean}} - u), $$

其中 $ \eta(s) $ 在 $ s=0 $ 附近线性增长，在无穷远处以 $ s^{-\gamma} $ 截断（$ \gamma \in (0,1) $），以控制其在无穷远处的增长行为，便于能量估计。

关键挑战来自应力项 $ \sigma $ 的形式。根据Hibler原始模型，定义Hibler变形张量 $ Tu $ 为：

$$ Tu := \varepsilon_D u + \frac{1}{\sqrt{2}e}\,\mathrm{tr}(\varepsilon u)\, \mathbf{1}_{2\times2}, $$

其中 $ \varepsilon u = \frac{1}{2}(\nabla u + \nabla u^\top) $ 是对称梯度，$ \varepsilon_D u $ 是其无迹部分，$ e \in (0,1) $ 是屈服椭圆的轴比。由此定义的粘性系数为：

$$ \zeta(\varepsilon,P) = \frac{P}{2\Delta(\varepsilon)}, \quad \Delta(\varepsilon) = |Tu|. $$

最终得到的应力表达式为：

$$ \sigma = \zeta(Tu)\, Tu = P \frac{Tu}{2|Tu|}. $$

注意到当 $ |Tu| \to 0 $ 时，$ \zeta \to \infty $，这导致应力在小变形区域趋于无穷大——这是典型的塑性行为。特别地，这种结构意味着应力仅依赖于变形的方向而非大小，类似于1-Laplace型结构或总变差流（total variation flow）。因此，该模型天然属于非光滑、不可微的能量泛函驱动系统。

然而，传统的分析方法通常采用正则化处理，例如将 $ |Tu| $ 替换为 $ \sqrt{\varepsilon + |Tu|^2} $，从而获得可微能量并使用极大正则性等技术（如 [BDHDH22] 中所做）。但此类正则化会抹除裂缝形成、剪切带等真实塑性现象。

本文的核心目标正是克服这一障碍：在不引入任何正则化的前提下，处理真正的奇异Hibler应力 $ \sigma = P \frac{Tu}{2|Tu|} $，即所谓“奇异极限”。

为此，作者必须跳出Sobolev空间框架。原因在于：若要求速度场 $ u \in W^{1,1}(\Omega;\mathbb{R}^2) $，则其梯度绝对连续，无法刻画裂纹或集中变形。而若尝试使用全变差空间 $ BV(\Omega;\mathbb{R}^2) $，又面临一个根本性障碍——Ornstein非不等式（Ornstein’s Non-Inequality）：不存在Korn型不等式在 $ L^1 $ 框架下成立，即无法从对称梯度 $ Eu $ 的Radon测度性质推出全梯度 $ \nabla u $ 的测度性。

因此，自然的选择是进入有界变形函数空间（Functions of Bounded Deformation, BD）：

$$ BD(\Omega) := \left\{ u \in L^1(\Omega;\mathbb{R}^2) : Eu \in \mathcal{M}_{\mathrm{fin}}(\Omega;\mathbb{R}^{2\times2}_{\mathrm{sym}}) \right\}, $$

其中 $ Eu = \varepsilon u $ 是对称分布梯度，取值为有限对称矩阵值Radon测度。注意 $ BV \subsetneq BD $，且 $ BD $ 不保证全梯度为测度。

在此空间中，原能量泛函

$$ F[u] = \int_\Omega |Tu|\,dx $$

不再适用于包含奇异部分的速度场。为此，作者引入松弛能量（relaxed energy）：

$$ F^*[u] = \int_\Omega F\left( \frac{dTu}{d\mathcal{L}^2} \right) d\mathcal{L}^2 + \int_\Omega F^\infty\left( \frac{dTs_u}{d|Ts_u|} \right) d|Ts_u| + \int_{\partial\Omega} F^\infty(-\mathrm{tr}_\Omega(u)\otimes T\nu_{\partial\Omega})\,d\mathcal{H}^1, $$

其中：

  • $ Tu = T[Eu] $ 是通过线性映射作用于 $ Eu $ 得到的Hibler变形测度；
  • $ Ts_u $ 是其奇异部分；
  • $ F^\infty(z) = \lim_{t\searrow 0} tF(tz) $ 是** recession function **（余函数），用于捕捉集中效应；
  • 边界项涉及迹算子 $ \mathrm{tr}_\Omega(u) $ 和边界法向 $ \nu_{\partial\Omega} $。

这个松弛能量允许能量沿低维集合（如裂纹或边界）集中，并保持下半连续性，是构造变分解的关键。

此外，由于考虑Dirichlet边界条件 $ u=0 $ on $ \partial\Omega $，但在能量方法中需允许比较函数具有非零边界值，这就引出了一个深刻的技术难题：如何用体积分逼近边界能量项？作者发展了一种新的非线性迹表达式的约化方案（reduction scheme for nonlinear trace expressions），实现了所谓的“体近似结果”（bulk approximation result），这是证明存在性的核心一步。

综上所述，本文的研究对象是一个基于真实Hibler应力的、定义在 $ BD(\Omega) $ 空间上的非线性演化系统，其主要特征包括：

  • 奇异粘塑性应力导致能量泛函不可微；
  • 解的空间正则性低于Sobolev类，允许出现裂纹和集中变形；
  • 能量松弛引入了体、奇异和边界三部分贡献；
  • 动力学表现为能量驱动的梯度流型演化。

主要定理或主要结果

本文的主要存在性结果体现在两个定理中：

定理 5.2 （能量驱动解的存在性）

设 $ \Omega \subset \mathbb{R}^2 $ 为有界Lipschitz域，$ T>0 $，质量 $ m>0 $ 和压力 $ P>0 $ 为常数。假设外力项满足适当可积性条件（如 $ f \in L^2(0,T;L^2(\Omega)) $），海洋力满足上述截断增长条件。则对于任意初始数据 $ u_0 \in BD(\Omega) $，存在一个能量驱动解 $ u: [0,T] \to BD(\Omega) $，使得：

  1. $ u \in L^\infty(0,T;L^2(\Omega)) \cap C_w([0,T];BD(\Omega)) $，
  2. $ u(0,\cdot) = u_0 $，
  3. $ u $ 满足松弛能量泛函 $ F^* $ 所诱导的梯度流不等式（见定义5.1），
  4. $ u $ 是动量平衡方程的弱解，即对所有测试函数 $ \varphi \in C_c^\infty((0,T)\times\Omega;\mathbb{R}^2) $，
$$ \int_0^T \int_\Omega m \partial_t u \cdot \varphi \,dxdt = \int_0^T \langle T^*(\sigma), \varphi \rangle dt + \cdots $$

成立（具体形式略去外力项）。

此定理首次在未正则化的Hibler应力下建立了全局时间存在的能量驱动解。

定理 5.19 （作为正则化问题极限的解）

设 $ u_\varepsilon $ 是正则化问题（将 $ |Tu| $ 替换为 $ \sqrt{\varepsilon + |Tu|^2} $）的解序列，则当 $ \varepsilon \to 0 $ 时，$ u_\varepsilon $ 在适当拓扑下收敛到某个 $ u $，该 $ u $ 是上述能量驱动解，满足定理5.2的所有性质。

这表明奇异极限下的解可以被物理上合理的正则化解逼近，增强了模型的稳健性。

此外，在第6节中，作者提出了更广义的“弱变分解”（Weak-Variational Solutions, WVS）概念，适用于质量非常数的情形，为未来工作奠定基础。

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

a) 核心分析工具

  1. 变分方法与能量驱动演化
    作者采用能量驱动解（energy-driven solution）框架，源自De Giorgi的最小作用原理和Ambrosio-Gigli-Savaré的度量梯度流理论。这类解不要求方程逐点成立，而是通过能量耗散不等式来刻画演化过程：

    $$ \mathcal{E}(u(t)) + \int_0^t \Phi(u'(s))\,ds \leq \mathcal{E}(v) + \int_0^t \Psi(u(s),v)\,ds, $$

    其中 $ \mathcal{E} $ 是能量泛函，$ \Phi $ 是耗散机制，$ \Psi $ 是功率项。这种方法特别适合处理非凸、非光滑能量系统。

  2. 函数空间理论：BD空间与弱-收敛*
    使用 $ BD(\Omega) $ 空间是本文的根本选择。其上的弱*-收敛（$ u_j \rightharpoonup^* u $ in $ BD $）指：

    • $ u_j \to u $ in $ L^1(\Omega) $，
    • $ Eu_j \rightharpoonup^* Eu $ in $ \mathcal{M}_{\mathrm{fin}} $（窄收敛）。

    该收敛方式下，能量泛函 $ F^* $ 具有下半连续性，这对取极限至关重要。

  3. Reshetnyak下半连续定理
    松弛能量 $ F^* $ 的构造依赖于经典结果：若 $ \mu_j \rightharpoonup^* \mu $ 为向量测度，则

    $$ \int_\Omega f^\infty\left( \frac{d\mu}{d|\mu|} \right) d|\mu| \leq \liminf_{j\to\infty} \int_\Omega f\left( \frac{d\mu_j}{d\mathcal{L}^n} \right) d\mathcal{L}^n, $$

    当 $ |\mu_j|(\Omega) \to \infty $ 时有效。这保证了奇异部分的能量不会在极限中消失。

  4. 三重逼近策略
    存在性证明采用三层逼近：

    • 初始值光滑化：$ u_0^{(k)} \in C^\infty \cap BD $，
    • 粘性正则化：加入 $ \delta \Delta u $ 项使问题强椭圆，
    • 应力正则化：替换 $ |Tu| \mapsto \sqrt{\varepsilon + |Tu|^2} $。

    然后依次令 $ k\to\infty $, $ \delta\to0 $, $ \varepsilon\to0 $，每步利用紧性提取收敛子列。

  5. Poincaré-Sobolev嵌入与刚体运动商空间
    利用 $ BD(\Omega) \hookrightarrow L^{2/(2-1)} = L^2(\Omega) $ 的紧嵌入（Lemma 3.3），确保了解的时间连续性与空间可积性。

b) 技术难点与创新技巧

  1. 边界项的体近似定理（Theorem 4.4）
    最关键技术突破是如下结果：对于 $ u \in BD(\Omega) $，其边界能量项

    $$ \int_{\partial\Omega} F^\infty(-\mathrm{tr}(u)\otimes T\nu)\,d\mathcal{H}^1 $$

    可以被一系列内部函数 $ u_j \in LD(\Omega) $ 的体积分逼近：

    $$ \lim_{j\to\infty} \int_\Omega F(Tu_j)\,dx = F^*[u]. $$

    该结果打破了Ornstein非不等式对边界浓度的限制，是首个将“非线性迹表达式”完全还原为体积分的工作。

  2. 非线性迹表达式的约化方案
    作者构造了一个精细的延拓算子与局部坐标变换，将边界附近的 $ Tu $ 展开为切向与法向分量的组合，再通过渐近分析分离出主导项。这一技巧独立于具体物理模型，具有普遍意义。

  3. Hibler-Anzellotti配对（Section 6）
    为了定义弱解中的 $ \langle T^*(\sigma), \varphi \rangle $，作者推广了Anzellotti关于$ L^\infty \cdot \mathrm{div}BV $ 配对的思想，建立 $ T^*(\sigma) \in L^2(0,T;W^{-1,\infty}) $ 与 $ BD $ 函数之间的对偶配对。

  4. 避免最大值原理与Moser迭代
    由于缺乏全梯度控制，传统抛物型高阶正则性工具失效。作者转而依赖能量估计与紧性方法，接受较低正则性但更贴近物理现实。

与之前工作的比较

a) 突破与改进

  1. 摆脱正则化假设
    此前多数工作（如 [BDHDH22], [LTT22]）均采用 $ \sigma_\varepsilon = \frac{Tu}{\sqrt{\varepsilon + |Tu|^2}} $ 类型的正则化，以便应用极大正则性或线性化技术。这本质上将问题限制在 $ W^{1,p} $ 空间内，排除了裂纹形成的可能。本文首次在 $ \varepsilon = 0 $ 的奇异极限下建立严格数学理论，真正实现了对塑性行为的建模。

  2. 提升解的物理真实性
    以往基于Sobolev空间的解自动具有连续梯度，不能反映海冰中常见的断裂、剪切带等现象。而本文解属于 $ BD $ 空间，允许 $ Eu $ 包含奇异测度（如沿曲线的集中），能自然描述这些地质特征。

  3. 拓宽函数空间框架
    将解空间从 $ W^{1,p} $ 推广至 $ BD $，不仅是正则性的降低，更是建模范式的转变。它标志着海冰模型分析从“光滑流体”视角转向“脆性材料”视角。

  4. 解决Ornstein障碍下的边界问题
    在 $ BV $ 或 $ BD $ 框架下，由于Ornstein非不等式，一般无法控制全梯度，使得边界迹难以处理。本文通过发展新的体近似技术，成功绕过此障碍，解决了长期悬而未决的边界能量集中问题。

  5. 统一多种塑性模型
    文中提出的框架不仅适用于标准Hibler模型，还可容纳Mohr-Coulomb型修正（图1所示），即小变形弹性、大变形塑性的混合行为。这为构建更真实的海冰本构关系提供了数学平台。

  6. 为数值模拟提供理论基础
    能量驱动解可通过正则化序列逼近，这为设计稳定数值格式（如有限元结合正则化）提供了理论支持。同时，“体近似”思想也可转化为数值积分策略。

总之，本文完成了从“理想化可微模型”到“真实奇异模型”的跨越，推动了粘塑性介质数学理论的发展，尤其在地球物理应用中具有里程碑意义。

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

1. Singular, finite-time $L^2$ attractors for odd, smooth solutions of Burgers equation on the torus

环面上Burgers方程奇光滑解的奇异有限时间 $L^2$ 吸引子

作者: Evan Miller

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09528v1

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

摘要翻译

中文标题: 环面上Burgers方程奇光滑解的奇异有限时间 $L^2$ 吸引子

摘要翻译: 本文证明了某个特定函数 $F$ 的正数倍——该函数在原点处具有跳跃间断，呈奇异性——对于一维无粘Burgers方程的一般奇光滑解而言，在 $L^2$ 空间中构成有限时间全局吸引子。此外，导出该结果的恒等式还为超临界范围 $0 < \alpha < \frac{1}{2}$ 内分形Burgers方程的有限时间爆破现象提供了另一种证明途径。该证明基于解与全局吸引子 $F$ 之间内积所定义的Lyapunov泛函的下界估计。我们还将证明，这一性质对一类更广泛的奇函数也成立，这类函数在区间 $(0, \pi)$ 上严格递增。

2. Comments on the paper `Modelling and nonclassical symmetry analysis of a complex porous media flow in a dilating channel'

关于论文《扩张通道中复杂多孔介质流动的建模与非经典对称性分析》的评述

作者: Roman Cherniha

类别: physics.flu-dyn, math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08656v1

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

摘要翻译

中文标题: 关于论文《扩张通道中复杂多孔介质流动的建模与非经典对称性分析》的评述

摘要翻译: 本文针对近期发表的论文《扩张通道中复杂多孔介质流动的建模与非经典对称性分析》（Physica D 481 (2025) 134834）进行评述，该文研究了一个描述通过多孔介质的非定常二维粘性不可压缩流体流动的模型。该研究的主要理论成果在于求解一个通过简化原模型导出的四阶偏微分方程的李对称性和非经典对称性。本文指出，该研究所得出的主要理论结果是不完整且具有误导性的。

不相关论文（基本信息）

1. Fast integral methods for the Neumann Green’s function: applications to capture and signaling problems in two dimensions

Neumann格林函数的快速积分方法：在二维捕获与信号问题中的应用

作者: Sanchita Chakraborty, Jeremy Hoskins, Alan E. Lindsay

类别: math.NA, math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08458v2

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

摘要: We present a high order numerical method for the solution of the Neumann Green’s function in two dimensions. For a general closed planar curve, our computational method resolves both the interior and exterior Green’s functions with the source placed either in the bulk or on the surface – yielding four distinct functions. Our method exactly represents the singular nature of the Green’s function by decomposing the singular and regular components. In the case of the interior function, we exactly prescribe an integral constraint which is necessary to obtain a unique solution given the arbitrary constant solution associated with Neumann boundary conditions. Our implementation is based on a fast integral method for the regular part of the Green’s function which allows for a rapid and high order discretization for general domains. We demonstrate the accuracy of our method for simple geometries such as disks and ellipses where closed form solutions are available. To exhibit the usefulness of these new routines, we demonstrate several applications to open problems in the capture of Brownian particles, specifically, how the small traps or boundary windows should be configured to maximize the capture rate of Brownian particles.

2. Microlocal analysis of the non-relativistic limit of the Klein–Gordon equation: Asymptotics

克莱因-戈登方程非相对论极限的微局部分析：渐近性

作者: Andrew Hassell, Qiuye Jia, Ethan Sussman, Andras Vasy

类别: math.AP, math-ph

发布日期: 2025-11-11

arXiv ID: 2511.08724v1

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

摘要: This is the less technical half of a two-part work in which we introduce a robust microlocal framework for analyzing the non-relativistic limit of relativistic wave equations with time-dependent coefficients, focusing on the Klein–Gordon equation. Two asymptotic regimes in phase space are relevant to the non-relativistic limit: one corresponding to what physicists call ``natural’’ units, in which the PDE is approximable by the free Klein–Gordon equation, and a low-frequency regime in which the equation is approximable by the usual Schrödinger equation. As shown in the companion paper, combining the analyses in the two regimes gives global estimates which are uniform as the speed of light goes to infinity. In this paper, we derive asymptotics from those estimates. Our framework differs from those in previous works in that ours is based on spacetime phase-space analysis.

3. Nonlinear Dirac equations on noncompact quantum graphs with potentials: Multiplicity and Concentration

具有势函数的非紧量子图上的非线性Dirac方程：多重性与集中现象

作者: Guangze Gu, Ziwei Li, Michael Ruzhansky, Zhipeng Yang

类别: math.AP, math-ph

发布日期: 2025-11-12

arXiv ID: 2511.09285v1

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

摘要: In this paper, we study the existence and multiplicity of solutions to the following class of nonlinear Dirac equations (NLDE) on noncompact quantum graphs:

$$ -i\,\varepsilon c\,σ_1\,\partial_x u + m c^2 σ_3 u + V(x)\,u = f(|u|)\,u, \quad x\in \mathcal{G}, \tag{P} $$

where $V:\mathcal{G}\to\mathbb{R}$ and $f:\mathbb{R}\to\mathbb{R}$ are continuous, $\varepsilon>0$ is a semiclassical parameter, $m>0$ denotes the mass, and $c>0$ the speed of light. Here $σ_1,σ_3$ are Pauli matrices, and $\mathcal{G}$ is a noncompact quantum graph. We prove that when $\varepsilon$ is sufficiently small, the number of solutions to $(P)$ is at least the number of global minima of $V$. Moreover, these solutions exhibit semiclassical concentration: as $\varepsilon\to0$, their concentration points approach the set of global minima of $V$.

4. Dynamical Formation of Apparent Horizons due to Boundary Effect in Vacuum Einstein Gravity

真空爱因斯坦引力中由于边界效应导致视界的动力学形成

作者: Puskar Mondal, Shing-Tung Yau

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

发布日期: 2025-11-12

arXiv ID: 2511.09508v1

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

摘要: We prove that an apparent horizon can form as a result of Einsteinian evolution in pure vacuum spacetime starting from regular initial data free of apparent horizons due to pure boundary effects. We adapt a Cauchy-double-null framework and use the boundary generalized mean curvature condition for the existence of an interior apparent horizon imposed by the author S-T Yau in \cite{yau}. In particular, we prove that the condition of \cite{yau} can be met dynamically starting from a configuration that does not verify the same through a focusing mechanism. This is the first part of a two-part sequence, and in the sequel, we will focus on explicitly constructing the Cauchy data.

5. Propagation of chaos for the Landau equation via microcanonical binary collisions

通过微正则二体碰撞研究朗道方程的混沌传播

作者: Kai Du

类别: math.AP, math.PR

发布日期: 2025-11-12

arXiv ID: 2511.09035v1

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

摘要: We develop a fully constructive, conservative, and collision-level realization of Kac’s program for the spatially homogeneous Landau equation across the full interaction range, including the Coulomb case. Our model is the microcanonical binary-collision (MBC) process: a reversible pure-jump $N$-particle Markov process that is Landau-native, realizing the grazing-collision mechanism via small conservative rotations of relative velocities. The analysis hinges on two critical structural pillars: a Fisher-information dissipation mechanism that extends the Guillén–Silvestre paradigm (Acta Math. 234:315-375, 2025) to a genuinely conservative particle system, yielding robust control of singular configurations, and a quantitative self-averaging principle that enforces a coherent deterministic emergence of the Landau flow from the microscopic dynamics. We prove propagation of chaos in the joint mean-field and grazing-collision limit ($N\to\infty, h\to 0$), identifying any limit point with the unique global solution to the Landau equation. Furthermore, we rigorously derive the Landau master equation as the grazing-collision limit of the MBC process. To the best of our knowledge, this provides the first fully conservative, Landau-native binary-collision model rigorously shown to produce the Landau equation over the entire interaction range.

6. Boundedness for the wave equation on $C^1$ stationary axisymmetric perturbations of Kerr

关于Kerr度规的$C^1$类定常轴对称扰动上波动方程解的有界性研究

作者: Yakov Shlapentokh-Rothman, Mihai Tohaneanu

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

发布日期: 2025-11-11

arXiv ID: 2511.08751v1

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

摘要: On the full range of sub-extremal Kerr exterior spacetimes we give a new proof of energy boundedness for high-frequency projections of solutions to the wave equation onto trapped frequencies. A key feature of the new estimate is that it circumvents the use of an integrated local energy decay (ILED) statement. As an illustration of the robustness of the estimate, we use it to establish energy boundedness for solutions to the wave equation on stationary and axisymmetric metrics which are merely $C^1$ close to a sub-extremal Kerr spacetime. We show explicitly that such perturbed metrics may possess stably trapped null geodesics, and thus one does not expect ILED statements to hold.

7. Singular solutions and bifurcation diagram of semilinear elliptic equations with general nonlinearity in two dimensions

二维具有一般非线性的半线性椭圆方程的奇异解与分歧图

作者: Hiroaki Kikuchi, Kenta Kumagai

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.08961v1

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

摘要: In this paper, we investigate semilinear elliptic equations with general exponential-type nonlinearities in two dimensions. For such nonlinearities, we establish two main results. The first is the construction of a singular solution. Recently, Fujishima, Ioku, Ruf, and Terraneo [10] proved the existence of singular solutions under certain assumptions for nonlinearities. We succeed in relaxing these conditions by providing the precise asymptotic form of a singular solution. Our second result concerns the bifurcation diagram of regular solutions. While the bifurcation structure has been extensively studied in three or higher dimensions, comparatively little was known in two dimensions until recently. In [18], the second author proved that the bifurcation curve possesses infinitely many turning points for supercritical analytic nonlinearities. In the present work, we refine this analysis by showing the bifurcation curve oscillates infinitely many times around some point, without assuming analyticity of the nonlinearities. The novelty of our approach lies in the introduction of a generalized Emden-type transformation.

8. A note on the a.e. second-order differentiability of rank-one convex functions

关于单凸函数几乎处处二阶可微性的一个注记

作者: Jonas Hirsch

类别: math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08397v2

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

摘要: In the Euclidean setting, the well-known Alexandrov theorem states that convex functions are twice differentiable almost everywhere. In this note, we extend this theorem to rank-one convex functions. Our approach is novel in that it draws more from viscosity techniques developed in the context of fully nonlinear elliptic equations. As a byproduct, the original Alexandrov theorem can essentially be reduced to the a.e. differentiability of one-dimensional monotone functions, as presented in the appendix.

9. A few techniques to achieve invisibility in waveguides

实现波导中隐身的几种技术

作者: Lucas Chesnel

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09172v1

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

摘要: The aim of this lecture is to consider a concrete problem, namely the identification of situations of invisibility in waveguides, to present techniques and tools that may be useful in various fields of applied mathematics. To be more specific, we will be interested in the propagation of acoustic waves in guides which are unbounded in one direction. In general, the diffraction of an incident field in such a structure in presence of an obstacle generates a reflection and a transmission characterized by some scattering coefficients. Our goal will be to play with the geometry, the frequency and/or the index material to control these scattering coefficients. We will explain how to:

  • develop a continuation method based on the use of shape derivatives to construct invisible defects;
  • exploit complex resonances located closed to the real axis to hid obstacles;
  • construct a non self-adjoint operator whose eigenvalues coincide with frequencies such that there are incident fields whose energy is completely transmitted. Our approaches will mainly rely on techniques of asymptotic analysis as well as spectral theory for self-adjoint and non self-adjoint operators. Most of the results will be illustrated by numerical experiments.

10. Rademacher’s Theorem for Calderon-Zygmund-type Spaces

Calderon-Zygmund型空间上的Rademacher定理

作者: Thomas Lamby

类别: math.FA, math.AP, math.CA

发布日期: 2025-11-12

arXiv ID: 2511.09159v1

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

摘要: Rademacher’s Theorem is a classical result stating that a Lipschitz function on $\mathbb{R}^d$ possesses a total differential almost everywhere. This implies that if $f$ is a function defined on $\mathbb{R}^d$ belonging to the Calderón–Zygmund space $T^\infty_1(x)$ for every $x \in \mathbb{R}^d$, then $f \in t^\infty_1(x)$ for almost every $x \in \mathbb{R}^d$. The main purpose of this paper is to extend this property to a broader functional framework. More precisely, we replace the Lipschitz condition by the assumption that $f$ belongs to a general weighted Calderón–Zygmund space $T^p_φ(x)$ for $x \in E$, where $E$ is a measurable subset of $\mathbb{R}^d$ and $φ$ is a weight with fractional indices. We then show that, under suitable assumptions, the function belongs to $t^p_φ(x)$ for almost every $x \in E$. The result can be obtained either by imposing additional hypotheses or by working within a suitably adapted $t^p_φ$ space. Whenever relevant, we also provide counterexamples illustrating the sharpness and the limitations of the statement.

11. On a partial data inverse problem for the semi-linear wave equation

关于半线性波动方程的一类偏数据反问题

作者: Boya Liu, Weinan Wang

类别: math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08794v1

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

摘要: We show that a partial Dirichlet-to-Neumann map, where the measurement set is arbitrarily small, uniquely determines the time-dependent nonlinearity of order three or higher in a semi-linear wave equation up to natural obstructions on a Lorentzian manifold with boundary. In particular, we do not impose any geometric or size restrictions on the measurement set. The proof relies on the technique of higher order linearization combined with the construction of Gaussian beams with reflections on the boundary.

12. Doubling variables and uniqueness of probability solutions to degenerate stationary Kolmogorov equations

变量加倍法与退化定常Kolmogorov方程概率解的唯一性

作者: V. I. Bogachev, S. V. Shaposhnikov, D. V. Shatilovich

类别: math.AP, math.PR

发布日期: 2025-11-11

arXiv ID: 2511.08781v1

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

摘要: We obtain sufficient conditions for the uniqueness of a probability solution to the stationary Kolmogorov equation with a degenerate diffusion matrix. We employ the method of doubling variables known in stochastic analysis directly to the Kolmogorov equation.

13. On Fractional Anisotropic Musielak-Sobolev Spaces with Applications to Nonlocal Eigenvalue Problems

关于各向异性Musielak-Sobolev分数阶空间及其在非局部特征值问题中的应用

作者: Mohammed Srati

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09256v1

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

摘要: In this paper, we introduce and study a new class of fractional modular function spaces, called \emph{Fractional Anisotropic Musielak–Sobolev Spaces}, which generalize both the fractional Anisotropic Orlicz–Sobolev spaces and the Anisotropic fractional Sobolev spaces with variable exponent. These spaces are designed to handle anisotropic and heterogeneous behaviors that naturally arise in nonlocal and nonlinear models. We develop their fundamental properties and embedding results, establishing a solid variational framework. As an application, we investigate a class of nonlocal anisotropic eigenvalue problems involving variable growth and direction-dependent fractional integro-differential operators. We prove the existence of eigenvalues by means of critical point theory and modular analysis. Our results extend and unify several existing models in the theory of nonlocal partial differential equations.

14. On the Dirichlet problem for the degenerate $k$-Hessian equation

关于退化 $k$-Hessian 方程的狄利克雷问题

作者: Yasheng Lyu

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09205v1

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

摘要: This paper investigates the existence of a global $C^{1,1}$ solution to the Dirichlet problem for the $k$-Hessian equation with a nonnegative right-hand side $f$, focusing on the required conditions for $f$. A well-known counterexample demonstrates the sharpness of the condition $f^{1/(k-1)}\in C^{1,1}(\overline{Ω_{0}})$ with $f\geq0$ in a domain $Ω_{0}\SupsetΩ$. We establish the existence for all $2\leq k\leq n-1$ under the condition $f^{1/(k-1)}\in C^{1,1}(\overline{Ω_{0}})$ with $f\geq0$ in $Ω_{0}$, and under an a priori assumption that $Δu\geq1$. Surprisingly, higher smoothness of $f$ permits a better exponent to be achieved. The condition $f^{3/(2k-2)}\in C^{2,1}(\overline{Ω_{0}})$ with $f\geq0$ in $Ω_{0}$ is also sharp, as demonstrated by the same counterexample. For the Monge-Ampère equation (the case $k=n$), we establish the existence under the optimal condition $f^{3/(2n-2)}\in C^{2,1}(\overline{Ω_{0}})$ with $f\geq0$ in $Ω_{0}$. For any $5\leq k\leq n-1$, we prove the existence under the condition $f^{3/(2k-2)}\in C^{2,1}(\overline{Ω_{0}})$ with $f\geq0$ in $Ω_{0}$, and under the a priori assumption $Δu\geq1$. Moreover, we obtain the existence for all $2\leq k\leq n-1$ under the condition $f^{3/(2k)}\in C^{2,1}(\overline{Ω_{0}})$ with $f\geq0$ in $Ω_{0}$.

15. Existence and uniqueness of solutions to SDEs with state-dependent variable exponent

具有状态依赖变指数的随机微分方程解的存在性与唯一性

作者: Mustafa Avci

类别: math.PR, math.AP

发布日期: 2025-11-12

arXiv ID: 2511.08882v1

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

摘要: We study a time-inhomogeneous, nonlinear stochastic differential equation whose drift and diffusion terms involve state-dependent variable exponents. This novel structure introduces flexibility, allowing the system to adapt dynamically to its current state and better capture real-world phenomena compared to classical models such as GBM and CEV. While this flexibility enhances modeling accuracy, it also creates analytical challenges. Using a fixed-point approach, we establish existence and uniqueness under suitable conditions. The model generalizes GBM and CEV while exhibiting distinct behavior for small states.

16. On the existence of non-negative weak solutions for $1D$ fourth order equations of gradient flow type

关于一类一维梯度流型四阶方程非负弱解的存在性研究

作者: Stefanos Georgiadis, Stefano Spirito

类别: math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08776v1

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

摘要: In this paper, we consider a family of one-dimensional fourth order evolution equations arising as gradient flows of the Korteweg energy, i.e. the $L^2$-norm of the first derivative of some power of the density. This family of equations generalizes the Quantum-Drift-Diffusion equation and the Thin-Film equation. We prove the global-in-time existence of {\em non-negative} weak solutions without requiring any upper bound on the exponent of the power of the density in the energy.

17. Ginzburg-Landau minimizers with high topological degrees in an annulus

环形区域中具有高拓扑度的金兹堡-朗道极小化子

作者: Amandine Aftalion, Rémy Rodiac

类别: math.AP, math-ph

发布日期: 2025-11-11

arXiv ID: 2511.08744v1

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

摘要: Motivated by recent experiments on fermionic rings, we study the asymptotic behaviour of minimizers of the Ginzburg-Landau (GL) energy in an annulus with a Dirichlet data which depends on the GL parameter on the outer boundary. We show that there is a critical degree of order $|\ln \varepsilon|$ under which the ground state displays a giant vortex and above which minimizers exhibit a combination of a giant vortex and vortices which tend to the outer boundary as the GL parameter tends to zero. Our analysis relies on the construction of suitable upper and lower bounds, on the extension to a slightly bigger annulus and on the minimization of the mean-field energy appearing in the lower bound. In order to be able to derive the minimum of this energy we use the symmetry of the domain and criticality with respect to inner variations.

18. Inverse problems for time-fractional Schrödinger equations

时间分数阶薛定谔方程的逆问题

作者: S. E. Chorfi, F. Et-tahri, L. Maniar, M. Yamamoto

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

发布日期: 2025-11-11

arXiv ID: 2511.08701v1

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

摘要: We study some inverse problems for time-fractional Schrödinger equations involving the Caputo derivative of fractional order $α\in (0,1)$. We prove refined uniqueness results from sets of positive Lebesgue measure for various problems by weakening the regularity of initial data.

19. Inverse Source Problem for a Nonlinear Parabolic Equation via Paralinearization

通过拟线性化方法求解非线性抛物型方程的反源问题

作者: Hu Xirui

类别: math.AP

发布日期: 2025-11-11

arXiv ID: 2511.08460v2

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

摘要: This thesis investigates an inverse source problem for a semilinear parabolic equation that includes both a quadratic reaction and a convection-type nonlinearity depending on the unknown function and its gradient. Based on partial boundary observations of the time derivative and its spatial gradient on a nonempty portion of the boundary, together with interior data at a fixed time, the goal is to recover an unknown spatial source term. On a short time interval, the nonlinear part is treated by Bony paralinearization technique, which transforms the original nonlinear problem into a linear parabolic equation with small, bounded coefficients and a negligible remainder term. By differentiating the equation with respect to time, a linearized model for the time derivative of the solution is obtained, which is suitable for the application of Carleman estimates. Using these tools, the work establishes conditional stability of Holder-logarithmic type for the recovery of the source term under the same observation geometry as in the corresponding linear case. The analysis combines techniques from Littlewood-Paley theory, Besov space estimates, and Carleman inequalities, showing that the stability structure of the linear inverse source problem extends naturally to the nonlinear setting under smallness and short-time assumptions.

20. The exterior Steklov problem for Euclidean domains

欧几里得区域上的外部Steklov问题

作者: Lukas Bundrock, Alexandre Girouard, Denis S. Grebenkov, Michael Levitin, Iosif Polterovich

类别: math.SP, math.AP

发布日期: 2025-11-12

arXiv ID: 2511.09490v1

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

摘要: We investigate the Steklov eigenvalue problem in an exterior Euclidean domain. First, we present several formulations of this problem and establish the equivalences between them. Next, we examine various properties of the exterior Steklov eigenvalues and eigenfunctions. One of our main findings is an Escobar-type lower bound for the first exterior Steklov eigenvalue on convex domains in dimensions three and higher. This bound is expressed in terms of the principal curvatures of the boundary and is sharp, with equality attained for a ball. Moreover, it implies the existence of a sequence of convex domains with fixed volume and the first exterior Steklov eigenvalues tending to infinity. This contrasts with the interior case, as well as with the two-dimensional exterior case, for which we show that an analogue of the Weinstock isoperimetric inequality holds.

21. Singularly Weighted X-ray Tensor Tomography

奇异加权X射线张量断层成像

作者: Jonathan Kay, François Monard

类别: math.AP

发布日期: 2025-11-12

arXiv ID: 2511.08871v1

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

摘要: If $d$ is a boundary defining function for the Euclidean unit disk and $I$ denotes the geodesic X-ray transform, for $γ\in (-1,1)$, we study the singularly-weighted X-ray transforms $I_m d^γ$ acting on symmetric $m$-tensors. For any $m$, we provide a sharp range decomposition and characterization in terms of a distinguished Hilbert basis of the data space, that comes from earlier studies of the Singular Value Decomposition for the case $m=0$. Since for $m\ge 1$, the transform considered has an infinite-dimensional kernel, we fully characterize this kernel, and propose a representative for an $m$-tensor to be reconstructed modulo kernel, along with efficient procedures to do so. This representative is based on a new generalization of the potential/conformal/transverse-tracefree decomposition of tensor fields in the context of singularly weighted $L^2$-topologies.