生成时间: 2025年10月17日 05:38:30

重点关注论文数量: 4

了解领域论文数量: 3

不相关论文数量: 25

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

1. Global stability for compressible isentropic Navier-Stokes equations in 3D bounded domains with Navier-slip boundary conditions

具有Navier滑移边界条件的三维有界区域中可压缩等熵Navier-Stokes方程的整体稳定性

作者: Yang Liu, Guochun Wu, Xin Zhong

类别: math.AP

发布日期: 2025-04-23

arXiv ID: 2504.17136v4

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

详细分析

中文标题: 具有Navier滑移边界条件的三维有界区域中可压缩等熵Navier-Stokes方程的整体稳定性

1. 研究对象和背景

本文研究的是三维有界光滑区域 $\Omega \subset \mathbb{R}^3$ 中，具有 Navier型滑移边界条件 的可压缩等熵Navier–Stokes方程（Compressible Isentropic Navier–Stokes Equations）的整体稳定性问题。这类方程描述了粘性、可压缩流体在受限空间中的运动规律，在空气动力学、气象学、工程流体力学等领域具有重要应用。

所考虑的控制方程如下：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho u) = 0, \\ (\rho u)_t + \mathrm{div}(\rho u \otimes u) + \nabla P = \mu \Delta u + (\mu + \lambda) \nabla \mathrm{div}\, u, \end{cases} \quad \text{in } \Omega \times (0,\infty), $$

其中：

  • $\rho = \rho(x,t) \geq 0$ 表示流体的密度；
  • $u = (u_1,u_2,u_3)$ 是速度向量场；
  • 压力 $P = a \rho^\gamma$，$\gamma > 1$，$a>0$，文中不失一般性设 $a=1$，即 $P = \rho^\gamma$；
  • $\mu > 0$ 和 $\lambda$ 分别为剪切粘性和体积粘性系数，满足物理条件：$2\mu + 3\lambda \geq 0$。

该系统定义在一个简单连通、光滑、有界的三维区域 $\Omega$ 上，并配备如下初始数据：

$$ \rho(x,0) = \rho_0(x), \quad \rho u(x,0) = \rho_0 u_0(x), \quad x \in \Omega. $$

特别地，作者考虑的是 Navier滑移边界条件：

$$ u \cdot n = 0, \quad \mathrm{curl}\, u \times n = 0, \quad \text{on } \partial\Omega \times (0,\infty), $$

其中 $n$ 是 $\partial\Omega$ 上的单位外法向量。这一边界条件比经典的 无滑移边界条件 $u = 0$ 更加适用于粗糙或带微结构的边界情形，允许流体沿边界的切向发生滑动，而法向仍保持无穿透（$u \cdot n = 0$）。值得注意的是，此处的滑移条件通过旋度项 $\mathrm{curl}\, u \times n = 0$ 来体现，这与常见的 Navier 摩擦型边界 $(D(u)n + \vartheta u)_{\rm tan} = 0$ 在特定几何下等价，但数学处理上更具挑战性。

从函数空间角度看，本文关注的是所谓的“强解”（strong solution），其正则性高于弱解，但仍允许包含真空状态（即 $\rho = 0$ 的点存在）。根据文中的定义（Definition 1.1），强解 $(\rho, u)$ 需满足一定的Sobolev空间正则性，例如：

  • $\rho \in C([0,T]; W^{1,q})$, $\rho_t \in C([0,T]; L^q)$ 对某个 $q \in (3,6]$，
  • $u \in C([0,T]; H^2) \cap L^2(0,T; W^{2,q})$,
  • $\sqrt{\rho} u_t \in L^\infty(0,T; L^2)$, $u_t \in L^2(0,T; H^1)$.

这种设定使得解足够光滑以进行高阶能量估计，同时又能容纳初始可能含真空的情形——这是当前研究的一个核心难点。

此外，文中引入了一个关键的能量泛函：

$$ C_0 \triangleq \int_\Omega \left( \frac{1}{2} \rho_0 |u_0|^2 + G(\rho_0) \right) dx, \quad \text{其中 } G(\rho) \triangleq \int_{\bar{\rho}}^\rho \frac{P(s) - P(\bar{\rho})}{s^2} ds, $$

且 $\bar{\rho}_0 \triangleq \frac{1}{|\Omega|} \int_\Omega \rho_0 dx$ 为平均密度，保持守恒（由质量守恒可得 $\int_\Omega \rho(x,t)dx = \int_\Omega \rho_0 dx$ 恒成立）。

本研究的核心问题是：当初始扰动不必小（即所谓“大初值”情形），且允许出现真空时，是否存在全局强解？这些解是否会在长时间趋于平衡态？收敛速率如何？

以往关于可压缩Navier-Stokes方程的研究多集中于两类情形：

  1. Cauchy问题（全空间 $\mathbb{R}^3$ 或环面 $\mathbb{T}^3$）：利用傅里叶分析、色散估计等工具获得衰减率；
  2. 无滑移边界条件下的有界域：此时速度在边界消失，便于使用Poincaré不等式等经典工具。

然而，对于一般有界区域中带有滑移边界条件的大初值、含真空解的整体稳定性，此前尚无结果。本文正是填补了这一空白，成为首个在该框架下建立指数稳定性的研究成果。

尤其值得注意的是，Navier滑移边界条件破坏了速度场在边界上的消失性质，导致传统的能量方法难以直接应用，必须发展新的技巧来控制边界项和非线性耦合效应。例如，在估计时间导数 $u_t$ 时，会出现边界积分项如 $\int_{\partial\Omega} F (u_t \cdot n) dS$，这类项无法直接消去，需借助滑移条件推导出新的恒等式（如 $ (u_t + (u\cdot\nabla n)\times u_\perp )\cdot n = 0 $）来进行精细估计。

综上所述，本文的研究对象是一个在三维有界光滑区域中，具有Navier滑移边界条件、允许真空出现的可压缩等熵Navier-Stokes系统，目标是证明其大初值强解的整体指数稳定性，并揭示真空结构的持久性。这一工作不仅拓展了已有理论的应用范围，也推动了对复杂边界条件下流体长期行为的理解。

2. 主要定理或主要结果

本文的主要成果由两个定理构成：

定理 1.1（整体指数稳定性）

假设初始数据 $(\rho_0 \geq 0, u_0)$ 满足：

$$ K \triangleq \|\rho_0 - \bar{\rho}_0\|_{L^2} + \|\sqrt{\rho_0} u_0\|_{L^2} + \|\nabla u_0\|_{H^1} < \infty. $$

设 $(\rho, u)$ 是问题 (1.1)–(1.3) 的一个全局强解，且满足密度的一致上界：

$$ \sup_{t \geq 0} \|\rho(\cdot, t)\|_{L^\infty} \leq \hat{\rho}, \quad \text{对某个正常数 } \hat{\rho}. $$

则存在仅依赖于 $\hat{\rho}$ 和 $K$ 的正常数 $C_1, \eta_1$，使得对所有 $t \geq 0$，

$$ \|(\rho - \bar{\rho}_0, \sqrt{\rho} u, \nabla u, \sqrt{\rho} \dot{u})(\cdot, t)\|_{L^2} \leq C_1 e^{-\eta_1 t}. $$

进一步，若初始密度远离零，即存在 $\rho_* > 0$ 使得 $\inf_{x \in \Omega} \rho_0(x) \geq \rho_*$，则存在依赖于 $\rho_*, \hat{\rho}, K$ 的常数 $C_2, \eta_2 > 0$，使得

$$ \|\rho(\cdot, t) - \bar{\rho}_0\|_{L^\infty} \leq C_2 e^{-\eta_2 t}. $$

此定理表明：只要密度被一致上有界（允许真空），强解将指数收敛到平衡态 $(\bar{\rho}_0, 0)$；若初始密度远离真空，则密度本身也在 $L^\infty$ 范数下指数衰减。

定理 1.2（真空不可消失）

在定理 1.1 所有条件下，若初始时刻存在真空，即 $\inf_{x \in \Omega} \rho_0(x) = 0$，则对任意 $t \geq 0$，

$$ \inf_{x \in \Omega} \rho(x,t) = 0. $$

这意味着：一旦初始分布中某点为真空（即使只是一个点），那么在整个演化过程中，真空状态永不消失。结合Gagliardo–Nirenberg不等式还可推出梯度 $\|\nabla \rho(\cdot,t)\|_{L^r}$ ($r>3$) 将随时间呈指数增长，暗示可能出现剧烈的空间振荡或奇异性传播。

这两个定理共同刻画了解的长期行为：既表现出宏观上的指数稳定性，又保留了微观上的结构刚性（真空持续存在）。

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

a) 核心分析工具

  1. 能量方法与Lyapunov泛函构造

    文章采用改进的能量估计框架，构建合适的Lyapunov型能量泛函来捕捉系统的耗散机制。不同于传统仅用基本能量 $E(t)$，作者引入修正能量：

    $$ M_1(t) = D \int_\Omega \left( \frac{1}{2} \rho |u|^2 + G(\rho) \right) dx - \int_\Omega \rho u \cdot B[\rho - \bar{\rho}_0] dx, $$

    其中 $B$ 是 Bogovskii算子（见下文），用于处理压力项与密度扰动之间的耦合。该泛函等价于 $\|(\rho - \bar{\rho}_0, \sqrt{\rho} u)\|_{L^2}^2$，并通过组合基本能量不等式与动量方程测试得到：

    $$ \frac{d}{dt} M_1(t) + C M_1(t) + \|\mathrm{div}\, u\|_{L^2}^2 + \|\mathrm{curl}\, u\|_{L^2}^2 \leq 0, $$

    从而得出低阶变量的指数衰减。

  2. Bogovskii算子

    在处理有界区域的压力估计时，Bogovskii算子 $B: L^p(\Omega) \to W^{1,p}_0(\Omega)$ 是关键技术之一。它满足 $\mathrm{div}\, B[f] = f - \bar{f}$，且具有连续性 $\|\nabla B[f]\|_{L^p} \leq C \|f\|_{L^p}$。作者利用它将 $\|P - \bar{P}\|_{L^2}$ 控制在 $\|\rho - \bar{\rho}\|_{L^2}$ 和 $\|\rho u\|_{L^4}$ 之下，进而实现非线性项的闭合估计。

  3. Div-Curl引理与椭圆正则性

    利用滑移条件 $u \cdot n = 0$ 和 $\mathrm{curl}\, u \times n = 0$，结合 div-curl系统估计（Lemma 2.5），可得：

    $$ \|\nabla u\|_{L^q} \leq C \left( \|\mathrm{div}\, u\|_{L^q} + \|\mathrm{curl}\, u\|_{L^q} \right). $$

    这一不等式绕过了速度在边界不为零带来的困难，是连接低阶与高阶范数的关键桥梁。

  4. 有效粘性通量 $F$ 与涡量 $\omega$ 的估计

    定义：

    $$ F \triangleq (2\mu + \lambda)\mathrm{div}\, u - (P - \bar{P}), \quad \omega \triangleq \nabla \times u, $$

    它们满足椭圆边值问题：

    $$ \Delta F = \mathrm{div}(\rho \dot{u}), \quad \frac{\partial F}{\partial n} = \rho \dot{u} \cdot n \quad \text{on } \partial\Omega. $$

    利用这一结构及 $L^p$ 椭圆理论，可导出：

    $$ \|\nabla F\|_{L^p} \leq C \|\rho \dot{u}\|_{L^p}, \quad \|\nabla \omega\|_{L^p} \leq C (\|\rho \dot{u}\|_{L^p} + \|\omega\|_{L^p}), $$

    从而将高阶导数控制在物质导数 $\dot{u}$ 上。

  5. Gagliardo–Nirenberg不等式与Sobolev嵌入

    多处使用插值不等式，如：

    $$ \|f\|_{L^p} \leq C \|f\|_{L^2}^{\theta} \|\nabla f\|_{L^2}^{1-\theta}, \quad p \in [2,6], $$

    特别是在边界条件下 $C=0$ 若 $f \cdot n = 0$ 或 $f \times n = 0$，这对消除常数项至关重要。

b) 技术性突破与难点克服

  1. 边界项处理：滑移诱导的恒等式

    关键技巧来自观察：由 $u \cdot n = 0$ 在 $\partial\Omega$ 上，可推出：

    $$ u \cdot \nabla u \cdot n = - u \cdot (\nabla n) \cdot u, $$

    进一步得到：

    $$ (\dot{u} + (u \cdot \nabla n) \times u_\perp ) \cdot n = 0 \quad \text{on } \partial\Omega, $$

    其中 $u_\perp = -u \times n$。这一恒等式使作者能对 $\dot{u}$ 应用Poincaré不等式，避免边界积分发散。

  2. 真空下的物质导数估计

    在估计 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 时，需处理形如 $\int_{\partial\Omega} F (\dot{u} \cdot n) dS$ 的项。借助上述边界恒等式，将其转化为内积形式并加以控制，最终建立形如：

    $$ \frac{d}{dt} \|\sqrt{\rho} \dot{u}\|_{L^2}^2 + \mu \|\nabla \dot{u}\|_{L^2}^2 \leq C \|\nabla u\|_{L^2}^2, $$

    再结合已知的 $\|\nabla u\|_{L^2}$ 小性（对大 $t$ 成立），构造出新的Lyapunov不等式，完成指数衰减证明。

  3. 密度下界保持

    为证 $L^\infty$ 收敛，需先证密度下界一致正。作者修改[43]的方法，利用 $G(\rho)$ 的凸性和已有的上界 $\hat{\rho}$，结合Gronwall不等式证明：若初始 $\rho_0 \geq \rho_* > 0$，则 $\rho(x,t) \geq c > 0$ 对所有 $t$ 成立。

  4. 真空不可消失性的反证法

    假设某时刻 $\inf \rho(x,t) > 0$，结合连续性与紧性论证，可推出 $\rho(x,t)$ 在某个邻域内严格正，再利用方程的结构（特别是沿粒子轨迹的ODE分析）导出矛盾，从而证明真空一旦出现便永久存在。

4. 与之前工作的比较

a) 突破与改进

本文相较于先前工作有多项显著推进：

  1. 首次在一般有界域中处理大初值+滑移边界+真空三重挑战的整体稳定性
    此前结果如[8]虽在Navier边界下建立全局存在性，但要求初值能量小；而[42,43]在环面或全空间中研究大解稳定性，依赖傅里叶分析。本文摆脱了小初值限制，并将结论推广至更现实的有界区域+滑移边界。

  2. 放宽对密度的先验假设
    Padula[36]要求密度上下均有正下界，后被[37]削弱。本文仅需上界一致有界即可获得 $L^2$ 指数收敛，大大放松了条件。

  3. 获得更强的收敛模态
    相比[8]中仅得 $L^r$ ($r < \infty$) 衰减，本文在额外正则性下获得 $L^\infty$ 指数收敛，且速率明确。

  4. 揭示真空结构的动力学刚性
    以往研究多避开真空或假设远离零，本文正面回答：真空一旦出现，永不消失。这是对解结构本质的新认知，具有深刻的物理意义——信息不会因扩散而完全抹除。

  5. 发展适用于滑移边界的新型能量估计技术
    提出利用 $u \cdot \nabla n \cdot u$ 结构处理边界项，为后续研究此类边界提供了范式。

综上，本文不仅在数学上实现了多重要素的统一处理，更在物理图像上深化了我们对可压缩流体长期演化的理解。

2. Global stability for the compressible isentropic magnetohydrodynamic equations in 3D bounded domains with Navier-slip boundary conditions

具有Navier滑移边界条件的三维有界区域中可压缩等熵磁流体动力学方程的整体稳定性

作者: Yang Liu, Guochun Wu, Xin Zhong

类别: math.AP

发布日期: 2025-05-26

arXiv ID: 2505.19749v3

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

详细分析

中文标题: 具有Navier滑移边界条件的三维有界区域中可压缩等熵磁流体动力学方程的整体稳定性

1. 研究对象和背景

本文研究的是在三维（3D）有界光滑区域 $\Omega \subset \mathbb{R}^3$ 中，带有 Navier型滑移边界条件 的 可压缩等熵磁流体动力学（Magnetohydrodynamics, MHD）方程组的整体稳定性问题。这类系统描述了在磁场作用下导电流体的运动行为，广泛应用于天体物理、等离子体物理以及工程中的高温流体模拟。

所考虑的数学模型由以下偏微分方程组构成：

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

其中未知函数分别为：密度 $\rho = \rho(x,t) \geq 0$，速度场 $u = (u_1,u_2,u_3)(x,t)$，磁感应场 $H = (H_1,H_2,H_3)(x,t)$，压力 $P = a\rho^\gamma$（$a > 0$, $\gamma > 1$），$\mu > 0$ 为剪切粘性系数，$\lambda$ 为体积粘性系数，满足物理约束：

$$ 2\mu + 3\lambda \geq 0. $$

而 $\nu > 0$ 是电阻率系数，在文中不失一般性地设为 $\nu = 1$。

该系统结合了 可压缩等熵Navier-Stokes方程 和 Maxwell电磁方程 的简化形式。电场 $E$ 可通过公式 $E = \nu \nabla \times H - u \times H$ 表示，因此未显式出现在系统中。

初始条件设定如下：

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

边界条件采用 Navier滑移类型，具体为：

$$ 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边界条件（即 $u|_{\partial\Omega} = 0$）相比，Navier滑移条件更贴近实际物理场景，如壁面粗糙或低摩擦界面的情况。然而这也带来了分析上的挑战，因为能量耗散结构不再直接控制速度梯度的 $L^2$ 范数，必须借助 广义Poincaré不等式 和 div-curl估计 来恢复正则性。

此外，文章所处的空间框架是定义在有界区域 $\Omega$ 上的强解空间。作者引入了如下Sobolev空间：

  • $L^p = L^p(\Omega)$,
  • $W^{k,p} = W^{k,p}(\Omega)$,
  • $H^k = W^{k,2}$, 并特别定义了一个满足边界条件的速度函数空间：
    $$ \widetilde{H}^1 \triangleq \left\{ v \in H^1(\Omega) \mid (v \cdot n)|_{\partial\Omega} = (\mathrm{curl}\, v \times n)|_{\partial\Omega} = 0 \right\}. $$

值得注意的是，本文关注的是 大初值解的整体稳定性，而非小扰动下的局部存在性或渐近行为。所谓“整体稳定性”，指的是当时间趋于无穷时，解以指数速率收敛到某个平衡态，即使初始数据并非接近稳态的小扰动。

平衡态在此为常数状态 $(\bar{\rho}_0, 0, 0)$，其中 $\bar{\rho}_0 = \frac{1}{|\Omega|} \int_\Omega \rho_0 dx$ 是平均密度。目标是证明：只要密度在时间上一致有上界（即 $\sup_{t \geq 0} \|\rho(\cdot,t)\|_{L^\infty} \leq \hat{\rho} < \infty$），则解将 在 $L^2$ 意义下指数衰减 到该平衡态。

这一类结果属于“有界性蕴含衰减”的研究范式——即若某关键物理量（这里是密度）保持有界，则系统自动趋于稳定。这在处理具有真空或大振荡初值的问题中尤为重要。

此外，论文还进一步假设初始密度远离真空（即 $\inf_{x \in \Omega} \rho_0(x) \geq \rho_* > 0$），从而得到更强的结果：密度和磁场在 $L^\infty$ 范数下也呈指数衰减。

综上所述，本文的研究对象是一个在 三维有界光滑区域 上的 可压缩等熵MHD系统，配备 Navier滑移边界条件，研究其 大初值强解的长时间行为与整体稳定性。核心难点在于克服边界的非消失性、非线性耦合项的影响，以及在没有小初值假设的情况下建立有效的能量耗散机制。

2. 主要定理或主要结果

论文的核心结论是以下定理：

定理 1.1（整体稳定性结果）
设初始数据 $(\rho_0 \geq 0, u_0, H_0)$ 满足

$$ K \triangleq \|\rho_0 - \bar{\rho}_0\|_{L^2} + \|\sqrt{\rho_0} u_0\|_{L^2} + \|\nabla u_0\|_{H^1} + \|\nabla H_0\|_{H^1} < \infty, \tag{1.9} $$

且 $(\rho, u, H)$ 是问题 (1.1)–(1.3) 的全局强解，并满足密度的一致上界条件：

$$ \sup_{t \geq 0} \|\rho(\cdot,t)\|_{L^\infty} \leq \hat{\rho}, \tag{1.10} $$

其中 $\hat{\rho} > 0$ 为常数。则存在仅依赖于 $\hat{\rho}$ 和 $K$ 的正常数 $C_1, \eta_1 > 0$，使得对所有 $t \geq 0$，

$$ \|(\rho - \bar{\rho}_0, \sqrt{\rho}u, H, \nabla u, \nabla H, \sqrt{\rho}\dot{u}, H_t)(\cdot,t)\|_{L^2} \leq C_1 e^{-\eta_1 t}. \tag{1.11} $$

进一步地，若初始密度远离真空，即存在 $\rho_* > 0$ 使得

$$ \inf_{x \in \Omega} \rho_0(x) \geq \rho_*, $$

则存在仅依赖于 $\rho_*, \hat{\rho}, K$ 的常数 $C_2, \eta_2 > 0$，使得

$$ \|(\rho - \bar{\rho}_0, H)(\cdot,t)\|_{L^\infty} \leq C_2 e^{-\eta_2 t}. \tag{1.12} $$

此外，作者指出此类强解的存在性已在先前工作 [3] 中被证明，只要初始能量足够小即可保证密度上界成立，从而满足条件 (1.10)。

作为推论，文章还得出关于真空持续性的结论：

推论 1.1
若初始密度包含真空点（即 $\inf_{x \in \Omega} \rho_0(x) = 0$），则对任意 $t \geq 0$，仍有

$$ \inf_{x \in \Omega} \rho(x,t) = 0. $$

这表明真空不会随时间消失，除非初始时刻就完全远离真空。

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

a) 核心分析学工具

本文采用了多种现代偏微分方程分析技术，主要包括：

  • 能量方法与Lyapunov泛函构造：这是全文最核心的技术路线。作者构建了一系列复合能量泛函（如 $M_1(t)$），将基本能量、有效粘性通量与Bogovskii算子相结合，形成一个等价于目标范数平方的正定泛函，并证明其满足形如 $\frac{d}{dt} M(t) + c M(t) \leq 0$ 的微分不等式，从而导出指数衰减。

  • Bogovskii算子的应用：由于密度偏离平均值 $\rho - \bar{\rho}_0$ 的积分可能非零，标准Poincaré不等式无法使用。为此，作者利用 Bogovskii算子 $\mathcal{B}$ 将压力项 $(P - \bar{P})(\rho - \bar{\rho}_0)$ 控制住。该算子能构造向量场使其散度等于给定零均值函数，且保持边界为零，是处理非均匀密度扰动的关键工具（见引理2.2）。

  • Gagliardo-Nirenberg不等式与Sobolev嵌入：用于控制非线性项，尤其是高阶导数与低阶项之间的插值估计。例如文中多次用到：

    $$ \|f\|_{L^p} \lesssim \|\nabla f\|_{L^2}^{\theta} \|f\|_{L^2}^{1-\theta},\quad p \in [2,6]. $$

  • div-curl系统估计（引理2.4）：对于向量场 $u$，若已知 $\mathrm{div}\, u$ 和 $\mathrm{curl}\, u$，可在有界区域上恢复 $\nabla u$ 的 $L^p$ 范数。这是处理Navier边界条件下速度正则性的关键技术之一。

  • 广义Poincaré不等式（引理2.3）：适用于满足 $(u \cdot n)|_{\partial\Omega} = 0$ 或 $(u \times n)|_{\partial\Omega} = 0$ 的向量场，允许在不假设零边界值的前提下控制 $L^p$ 范数由梯度主导。

  • Gronwall不等式与闭合论证：整个证明采用逐步提升策略：先证低阶项（密度、速度、磁场）指数衰减，再利用此信息改进高阶导数（梯度、时间导数）的能量估计，最终实现全范数控制。

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

  1. 处理边界项的新方法：
    在估计 $\int F_t (u_t \cdot n) dS$ 类似边界积分时，传统方法失效。作者受文献[2]启发，利用恒等式：

    $$ u \cdot \nabla u \cdot n = -u \cdot \nabla n \cdot u,\quad \text{on } \partial\Omega, $$

    推出修正后的物质导数满足 $(\dot{u} + (u \cdot \nabla n) \times u_\perp) \cdot n = 0$，从而将边界项转化为可控形式，关键是控制 $\int_\partial\Omega (u \cdot \nabla n \cdot u) F dS$，并通过迹定理与Sobolev嵌入完成估计。

  2. 有效粘性通量 $F$ 与涡量 $\omega$ 的精细估计：
    定义

    $$ F \triangleq (2\mu + \lambda)\mathrm{div}\, u - (P - \bar{P}) - \frac{1}{2}|H|^2,\quad \omega \triangleq \nabla \times u, $$

    这两个量在控制 $\nabla u$ 中起决定性作用。通过椭圆估计（引理2.4），可得：

    $$ \|\nabla u\|_{L^p} \lesssim \|\rho \dot{u}\|_{L^2} + \|H \cdot \nabla H\|_{L^2} + \cdots, $$

    使得即使在滑移边界下也能获得速度梯度的高阶控制。

  3. 材料导数 $\dot{u} = u_t + u \cdot \nabla u$ 的运用：
    引入 $\dot{u}$ 可自然吸收对流项，使动量方程变为：

    $$ \rho \dot{u} = \cdots, $$

    更便于进行高阶能量估计。尤其在估计 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 的时间衰减时，这种方法避免了直接处理 $u_t$ 与非线性项的复杂耦合。

  4. 时间反演思想与小时间后启动Lyapunov不等式：
    由于初期解可能较大，作者观察到当 $t$ 足够大时 $\|\nabla u\|_{L^2}$ 已充分小，此时可通过连续性闭合能量不等式，构造出有效的Lyapunov泛函（如(3.28)式），再结合Gronwall引理推出指数衰减。

  5. 密度下界估计：
    为获得 $L^\infty$ 收敛性，需确保密度远离真空。作者修改了[17,25]的方法，利用阻尼机制（damping mechanism）证明：若初始密度有正下界，则在整个演化过程中仍保持正下界，进而可用最大值原理或比较定理控制 $\|\rho - \bar{\rho}_0\|_{L^\infty}$。

4. 与之前工作的比较

本文相较于已有研究成果有多项显著突破：

a) 改进之处

  1. 从环面到一般有界域的推广：
    Zhu-Zi [28] 的工作仅限于三维环面 $\mathbb{T}^3$，依赖周期性结构和傅里叶分析工具。本文首次在 一般光滑有界区域 中建立同类结果，更具普适性。

  2. 去除粘性系数限制 $2\mu > \lambda$：
    许多早期工作（包括[28]）要求 $2\mu > \lambda$ 以保证某些能量结构正定。本文通过更精细的估计完全摆脱了这一限制，仅需物理合理的 $2\mu + 3\lambda \geq 0$。

  3. 无需 $L^6$-范数对磁场的统一控制：
    [28] 要求 $\|H\|_{L^6}$ 一致有界作为额外前提，而本文仅需密度有上界即可推出所有量的指数衰减，大大弱化了假设。

  4. 取消零总动量假设：
    Torus情形常需 $\int_\Omega \rho u dx = 0$ 来消除刚体运动模态。本文在有界域中通过几何与边界条件自然抑制此类模式，无需额外限制。

  5. 首次针对大解在一般有界域的整体稳定性结果：
    文章宣称这是 第一个在三维一般有界区域中关于可压缩MHD大强解的整体稳定性结果，填补了领域空白。

b) 结果强度对比

  • 相比于 Chen-Huang-Peng-Shi [3] 在相同边界条件下的结果，本文不仅得到了相同阶的 $L^2$ 指数衰减，还进一步获得了 $\|\sqrt{\rho} \dot{u}\|_{L^2}$ 和 $H_t$ 的衰减，说明解的正则性更高。
  • 相比 Wu-Zhang-Zou [23] 的弱解衰减结果，本文处理的是 强解，且衰减速率为 指数级 而非代数衰减。
  • 对比 Li-Xu-Zhang [15] 的小能量解存在性，本文不要求初值接近平衡态，而是接受 大初值，只要密度有上界即可。

总之，本文在 适用范围更广（一般有界域）、假设更弱（无 $L^6(H)$ 控制、无 $2\mu > \lambda$）、结论更强（指数衰减至 $L^\infty$）等方面全面超越前人工作，代表了当前可压缩MHD稳定性理论的重要进展。

3. Quantitative stability of a class of explicit steady Euler flows in a disk

一类显式定常欧拉流在圆盘内的定量稳定性

作者: Fatao Wang, Guodong Wang

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14394v1

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

详细分析

中文标题: 一类显式定常欧拉流在圆盘内的定量稳定性

1. 研究对象和背景

本文研究的是二维不可压缩理想流体在单位圆盘 $D \subset \mathbb{R}^2$ 上的 欧拉方程（Euler equations）的解的 轨道稳定性，特别是针对一类具有显式结构的定常解。具体而言，作者关注的是该类定常解在 $L^2$ 范数下的定量稳定性估计，即当初始扰动足够小时，解在演化过程中与原始定常解的旋转轨道之间的距离可以被精确控制。

所考虑的流体模型为密度恒为1的理想流体，在区域 $D = \{x \in \mathbb{R}^2 : |x| < 1\}$ 内运动，并满足边界上的不可穿透条件（impermeability condition），即速度场在边界 $\partial D$ 上的法向分量为零。这一物理条件自然地对应于速度场由流函数通过旋度关系生成。

论文采用的是涡度形式的二维欧拉方程：

$$ \begin{cases} \partial_t \omega + v \cdot \nabla \omega = 0, & t \in \mathbb{R},\ x = (x_1,x_2) \in D, \\ v = \nabla^\perp G\omega, \\ \omega(0,\cdot) = \omega_0, \end{cases} $$

其中：

  • $\omega(t,x)$ 是流体的涡度（vorticity），定义为 $\omega = \nabla \times v = \partial_{x_2}v_1 - \partial_{x_1}v_2$；
  • $v(t,x)$ 是速度场；
  • $\nabla^\perp := (\partial_{x_2}, -\partial_{x_1})$ 表示逆时针旋转90度的梯度算子；
  • $G$ 是负拉普拉斯算子 $-\Delta$ 在 $D$ 上带有零狄利克雷边界条件的逆算子，即对于任意函数 $\omega$，$G\omega$ 满足：
    $$ \begin{cases} -\Delta (G\omega) = \omega & \text{in } D, \\ G\omega|_{\partial D} = 0. \end{cases} $$

这个系统描述了理想流体中涡度沿粒子轨迹守恒的动力学行为。由于 $v = \nabla^\perp G\omega$，速度场自动满足不可压缩性 $\nabla \cdot v = 0$ 和边界不可穿透条件。

一个关键的事实是：该欧拉系统拥有多个守恒量，这些量在光滑解的演化下保持不变。文中明确指出并利用了以下三个核心守恒律：

  1. 动能（Kinetic Energy）：
    $$ E(\omega) := \frac{1}{2} \int_D \omega G\omega \, dx; $$
  2. 涡量分布测度守恒（Level Set Measure Conservation）： 对任意实数 $s$，集合 $\{x \in D \mid \omega(t,x) > s\}$ 的二维勒贝格测度不随时间变化；
  3. 流体冲量矩（Moment of Fluid Impulse）：
    $$ I(\omega) := \int_D |x|^2 \omega(x)\, dx. $$

由第(2)条可推出一个重要推论：涡度的各阶 $L^p$ 范数均守恒，尤其是enstrophy（涡量平方积分）：

$$ J(\omega) := \int_D \omega^2 dx $$

在整个演化过程中保持不变。

本文的研究重点是一类特殊的显式定常解。所谓“定常”是指 $\partial_t \omega = 0$，即涡度不随时间变化。这类解构成了动力系统的平衡态或周期轨道的基础。

设 $J_n(s)$ 为第一类 $n$ 阶贝塞尔函数，$j_{n,k}$ 为其第 $k$ 个正零点，记 $j := j_{1,1} \approx 3.831706$。定义如下三维线性空间：

$$ V := \mathrm{span}\left\{ J_0(jr),\ J_1(jr)\cos\theta,\ J_1(jr)\sin\theta \right\}, $$

其中 $(r,\theta)$ 是极坐标。根据文献[11]，该空间中的任意元素 $\bar{\omega} \in V$ 都是欧拉方程的一个定常解。

特别地，作者考虑的是 $V$ 中的一类旋转轨道（rotational orbit）。固定实数 $A,B \in \mathbb{R}$ 且 $B \geq 0$，定义轨道：

$$ O := \left\{ A J_0(jr) + B J_1(jr)\cos(\theta + \alpha) \mid \alpha \in \mathbb{R} \right\}. $$

这表示将某个特定方向的模态绕原点旋转所得的所有函数构成的集合。显然，这样的轨道在 $L^2(D)$ 中是紧致的。

文章的目标是对该轨道的**$L^2$-轨道稳定性给出定量估计**。给定一个光滑解 $\omega(t,x)$，令 $\omega_t := \omega(t,\cdot)$，并用 $\mathrm{dist}_2(v,K)$ 表示函数 $v$ 到集合 $K$ 在 $L^2(D)$ 中的距离：

$$ \mathrm{dist}_2(v,K) := \inf_{w \in K} \|v - w\|_{L^2(D)}. $$

已有的定性结果（Theorem A）表明：对任意 $1 < p < \infty$ 和 $\epsilon > 0$，存在 $\delta > 0$，使得只要初值满足 $\mathrm{dist}_p(\omega_0, O) < \delta$，就有 $\mathrm{dist}_p(\omega_t, O) < \epsilon$ 对所有时间成立。这是经典的轨道稳定性结论，基于紧性论证和守恒量的组合。

本文的主要贡献在于将此结果提升至定量水平，特别是在 $p=2$ 的情形下，明确揭示了 $\delta$ 与 $\epsilon$ 之间的依赖关系，从而提供更强、更具实用性的稳定性控制。

此外，问题所在的函数空间也值得注意。作者主要工作在 $L^2(D)$ 空间及其子空间上，尤其是均值为零的子空间：

$$ \mathring{L}^2(D) := \left\{ v \in L^2(D) \mid \int_D v\,dx = 0 \right\}. $$

这是因为 $J_0(jr)$ 在单位圆盘上的积分非零，但整体分析中经常需要处理均值部分与零均值部分的分离。事实上，$V$ 并不完全包含在 $\mathring{L}^2(D)$ 中，因为 $J_0(jr)$ 本身有非零平均。

区域 $D$ 是单位圆盘，具有高度的径向对称性，这也导致特征函数可以用贝塞尔函数展开。这种几何结构使得许多谱分析工具（如Poincaré型不等式）可以显式计算或估计。

综上所述，本文的研究对象是一个定义在单位圆盘上的二维欧拉方程，其动力学受多个守恒量约束，目标是分析一类由贝塞尔函数构成的显式定常解在其旋转轨道附近的长期稳定性行为，并通过精细的能量-Casimir方法获得定量的 $L^2$ 轨道稳定性估计。

2. 主要定理或主要结果

本文的核心成果是如下定理，它给出了关于旋转轨道 $O$ 的定量 $L^2$-轨道稳定性：

定理 1.1（Quantitative $L^2$-Orbital Stability）
设 $O$ 为由公式 (1.5) 定义的旋转轨道：

$$ O := \left\{ A J_0(jr) + B J_1(jr)\cos(\theta + \alpha) \mid \alpha \in \mathbb{R} \right\}, $$

其中 $A,B \in \mathbb{R}$ 且 $B \geq 0$。则存在一个通用常数 $C > 0$（与 $A,B,\omega$ 无关），使得对任意 $\varepsilon > 0$ 和任意光滑解 $\omega(t,x)$ 满足 $\mathrm{dist}_2(\omega_0, O) \leq \varepsilon$，有：

(i) 若 $B \ne 0$，则对所有 $t \in \mathbb{R}$，

$$ \mathrm{dist}_2(\omega_t, O) \leq C B^{-1} \left( A^2 + B^2 \right)^{1/2} \varepsilon + C B^{-1} \varepsilon^2; $$

(ii) 若 $B = 0$，则对所有 $t \in \mathbb{R}$，

$$ \mathrm{dist}_2(\omega_t, O) \leq C |A|^{1/2} \varepsilon^{1/2} + C \varepsilon. $$

该定理的关键意义在于：

  • 当 $B \ne 0$ 时，扰动增长被控制在 $O(\varepsilon)$ 阶，尽管系数中含有 $B^{-1}$，暗示当 $B \to 0$ 时稳定性退化；
  • 当 $B = 0$ 时，即定常解完全径向对称（仅含 $J_0(jr)$ 项），稳定性降为 $O(\varepsilon^{1/2})$，说明此时稳定性显著减弱。

此外，文中还证明了一个辅助性的空间 $V$ 的定量稳定性命题，作为主要定理的基石：

命题 3.1（Quantitative Stability of $V$ under Mean-Zero Perturbations）
设 $\omega(t,x)$ 是欧拉方程的光滑解，且 $\omega_t \in \mathring{L}^2(D)$（即均值为零）。则对所有 $t \in \mathbb{R}$，

$$ \mathrm{dist}_2(\omega_t, V) \leq \sqrt{ \frac{\Lambda_2}{\Lambda_2 - \Lambda_1} } \mathrm{dist}_2(\omega_0, V), $$

其中 $\Lambda_1 = j^2$，$\Lambda_2 > \Lambda_1$ 分别是某个带均值约束的拉普拉斯特征值问题的第一和第二特征值。

此命题表明：若初始涡度接近空间 $V$，则其演化始终不会远离 $V$，且偏差可通过初始偏差线性控制。

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

a) 核心分析学工具

本文采用了经典的能量-Casimir方法（Energy-Casimir method），这是由Vladimir Arnold发展起来的一种用于研究流体力学中定常解稳定性的变分框架。其基本思想是构造一个结合动能与其他守恒量（如enstrophy）的复合泛函，并利用其极值性质来控制扰动的增长。

具体来说，作者引入了一个修正的能量泛函：

$$ EC(v) := 2E(v) - \frac{1}{\Lambda_1} J(v) = \int_D v Gv\,dx - \frac{1}{\Lambda_1} \int_D v^2 dx, $$

其中 $\Lambda_1 = j^2$ 是关键的谱参数。由于 $E$ 和 $J$ 均为守恒量，故 $EC(\omega_t)$ 在时间演化下保持不变。

另一个核心技术是谱投影与Poincaré型不等式。作者考虑如下带边界和均值约束的特征值问题：

$$ \begin{cases} -\Delta u = \Lambda u & \text{in } D, \\ u|_{\partial D} = c \quad (\text{自由常数}), \\ \int_D u\,dx = 0. \end{cases} $$

对应的算子 $T: \mathring{L}^2(D) \to \mathring{L}^2(D)$ 定义为 $Tv = Gv - \frac{1}{|D|}\int_D Gv\,dx$，其特征值满足 $0 < \Lambda_1 < \Lambda_2 < \cdots \to \infty$。引理2.2断言 $\Lambda_1 = j^2$，且相应的特征空间正是 $V$。

由此得到两个关键的Poincaré型不等式：

$$ \int_D v Gv\,dx \leq \frac{1}{\Lambda_1} \int_D v^2 dx, \quad \forall v \in \mathring{L}^2(D), $$

等号成立当且仅当 $v \in V$；

$$ \int_D v Gv\,dx \leq \frac{1}{\Lambda_2} \int_D v^2 dx, \quad \forall v \in V^\perp \cap \mathring{L}^2(D). $$

这些不等式允许将 $Gv$ 的“能量”与 $L^2$ 范数联系起来，从而实现对扰动的精细估计。

此外，作者使用了正交分解技术：将当前涡度 $\omega_t$ 投影到 $V$ 上得 $\hat{\omega}_t = P_1 \omega_t$，剩余部分 $g_t = \omega_t - \hat{\omega}_t$ 属于 $V^\perp$。通过守恒量 $EC$ 的不变性，结合上述不等式，可得：

$$ \|g_t\|_{L^2}^2 \leq \frac{\Lambda_2}{\Lambda_2 - \Lambda_1} \|g_0\|_{L^2}^2, $$

从而完成对 $\mathrm{dist}_2(\omega_t, V)$ 的定量估计。

b) 解决技术难题的关键技巧

本文最精巧的技术在于如何从 $\mathrm{dist}_2(\omega_t, V)$ 进一步推导出 $\mathrm{dist}_2(\omega_t, O)$，即从整个空间 $V$ 收敛到具体的轨道 $O$。

为此，作者采用两步策略：

  1. 将 $\omega_t$ 分解为 $\hat{\omega}_t \in V$ 与余项；
  2. 控制 $\hat{\omega}_t$ 到轨道 $O$ 的距离。

关键难点在于：虽然 $V$ 是有限维的，但轨道 $O$ 只是其中的低维子集（一维流形）。为了控制 $\hat{\omega}_t$ 接近 $O$，必须利用守恒量来限制其参数漂移。

作者巧妙地使用了两个守恒量：

  • 冲量矩 $I(\omega)$：用于控制 $J_0(jr)$ 项的系数 $A$；
  • enstrophy $J(\omega)$：用于控制 $J_1(jr)\cos\theta$ 项的系数 $B$。

具体而言，设 $\tilde{\omega}_t \in O$ 为 $\omega_t$ 在轨道上的最近点，则通过三角不等式和投影估计，有：

$$ \|\omega_t - \tilde{\omega}_t\|_{L^2} \lesssim \varepsilon + \|\hat{\omega}_t - \tilde{\omega}_t\|_{L^2}. $$

而 $\hat{\omega}_t$ 和 $\tilde{\omega}_t$ 同属 $V$，可分别写成：

$$ \hat{\omega}_t = A' J_0 + B' J_1 \cos(\theta+\alpha'),\quad \tilde{\omega}_t = A J_0 + B J_1 \cos(\theta+\alpha). $$

利用 $I(\omega_t) = I(\omega_0)$ 和 $I$ 在旋转下的不变性，结合Hölder不等式，证明了：

$$ |I(\tilde{\omega}_t) - I(\hat{\omega}_t)| \leq C \varepsilon \Rightarrow |A - A'| \leq C \varepsilon. $$

类似地，利用 $J(\omega_t) = J(\omega_0)$ 得到：

$$ |B^2 - B'^2| \leq C(|A|^2 + B^2)^{1/2} \varepsilon + C \varepsilon^2. $$

当 $B \ne 0$ 时，利用 $B + B' \geq B$ 可得 $|B - B'| \leq C B^{-1} (\cdots)$；当 $B = 0$ 时，则只能得到 $|B'| \leq C \varepsilon^{1/2}$，从而出现半次幂衰减。

这一技巧的核心在于：用全局守恒量反推局部参数的稳定性，避免了复杂的动态分析。

4. 与之前工作的比较

a) 突破与改进

本文相对于已有工作（特别是文献[11]中的Theorem A）的主要突破在于：

  • 从定性到定量：Theorem A仅断言存在性意义上的稳定性（“对任意 $\epsilon$ 存在 $\delta$”），而本文给出了 $\delta$ 与 $\epsilon$ 的显式依赖关系，实现了真正的可计算稳定性。
  • 方法革新：Theorem A依赖于紧性论证（compactness argument），本质上是非构造性的；而本文采用定量分析+守恒量控制参数的方法，更具解析性和透明度。
  • 揭示稳定性退化机制：本文首次明确指出，随着定常解的径向对称性增强（即 $B \to 0$），稳定性指数从 $O(\varepsilon)$ 退化为 $O(\varepsilon^{1/2})$，说明更对称的流动反而更不稳定。

这一点在Remark 4.5中被强调，并与Elgindi[7]在环面上的结果相呼应：越对称的流（如纯剪切流）越容易失稳。

此外，本文的方法相较于高阶Casimir泛函（如四次型）更为简洁——仅使用线性的冲量矩 $I(\omega)$ 即可完成参数控制，大大简化了计算复杂度。

b) 局限与未来方向

尽管取得了重要进展，本文仍有局限：

  • 结果仅适用于 $L^2$ 范数，难以推广到其他 $L^p$ 空间；
  • 所考虑的定常解属于特定低维空间 $V$，尚不能处理更高模态或更复杂结构的解；
  • 对 $B=0$ 情形无法改进至 $O(\varepsilon)$，作者在Remark 4.5中指出这是本质困难，即使引入其他 $C^1$ Casimir 泛函也无法改善。

然而，这种方法已被成功应用于球面S²上的Rossby-Haurwitz波（见Delgadino-Melzi[6]），显示出良好的可移植性。未来可能的方向包括：

  • 将此方法推广至任意形状的二维环面（flat torus）上的第一特征模态；
  • 研究非线性不稳定性阈值；
  • 探索在数值模拟中验证此类定量估计的可能性。

总之，本文通过精巧地结合经典守恒律与谱分析，为流体力学中定常解的稳定性理论提供了新的定量视角，标志着从“是否存在稳定”迈向“稳定到什么程度”的重要一步。

4. Stability and large-time behavior for the N-Dimensional Euler-FENE dumbbell model near an equilibrium

N维欧拉-FENE哑铃模型在平衡态附近的稳定性及大时间行为

作者: Zheng-an Yao, Ruijia Yu

类别: math.AP, Primary: 35Q30, Secondary: 76A10, 35B40, 35Q30

发布日期: 2025-07-05

arXiv ID: 2507.03962v2

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

详细分析

中文标题: N维欧拉-FENE哑铃模型在平衡态附近的稳定性及大时间行为

1. 研究对象和背景

本文研究的是 N维欧拉-FENE哑铃模型（Euler-FENE dumbbell model） 在其稳态解 $(0, \psi_\infty)$ 附近的小扰动的全局存在性、稳定性以及大时间衰减性质。该模型描述了一类典型的流体-聚合物耦合系统，其中聚合物分子被理想化为由两个“珠子”通过弹性弹簧连接而成的“哑铃”结构，其内部构型由向量 $R \in B(0,1) \subset \mathbb{R}^N$ 表示。

系统的宏观动力学由流体速度场 $u(x,t)$ 描述，而微观构型则通过分布函数 $\Psi(x,R,t)$ 来刻画，后者满足一个Fokker-Planck型方程。整个系统是一个微宏观耦合系统，结合了不可压缩流体的动力学与聚合物构型演化的统计规律。

原始的 FENE 模型包含粘性项（即 Navier-Stokes 型速度方程），但在本工作中，作者考虑的是无速度耗散的情形，即将粘性系数设为零，从而得到所谓的 Euler-FENE 模型：

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u = -\nabla p + \mathrm{div}\, \tau, \\ \mathrm{div}\, u = 0, \\ \partial_t \Psi + (u \cdot \nabla_x)\Psi = \nu \Delta_R \Psi + \mathrm{div}_R\left[ -(\nabla u) R \Psi + \nabla_R U \Psi \right], \\ \tau(x,t) = \int_{B} R \otimes (\nabla_R U) \Psi(x,R,t)\, dR, \\ u(x,0) = u_0(x),\quad \Psi(x,R,0) = \Psi_0(x,R), \end{cases} \tag{1} $$

其中：

  • $x \in \mathbb{R}^N$, $N \geq 2$,
  • $R \in B := B(0,1)$ 是单位球内的聚合物伸展向量，
  • 势能函数取为 $U(R) = -k \log(1 - |R|^2)$，对应经典的 FENE 势，保证了聚合物有限可伸展性；
  • $\nu > 0$ 是中心质量扩散率（center-of-mass diffusion rate）；
  • $\beta = 1$ 已归一化；
  • $\tau$ 是聚合物对流体施加的额外应力张量；
  • 边界条件为：在 $\partial B(0,1)$ 上，$( -(\nabla u)R \Psi + \nabla_R \Psi + \nabla_R U \Psi ) \cdot n = 0$，以确保概率守恒；
  • 平衡态分布为：
    $$ \psi_\infty(R) = \frac{e^{-U(R)}}{\int_B e^{-U(R)}\,dR} = \frac{(1 - |R|^2)^k}{\int_B (1 - |R|^2)^k\,dR}, $$
    它满足 $L \psi_\infty = 0$，其中 $L \psi := \mathrm{div}_R(\nabla_R \psi + \nabla_R U \psi)$ 是 Fokker-Planck 算子。

本文关注系统 (1) 在稳态 $(u, \Psi) = (0, \psi_\infty)$ 附近的扰动行为。令：

$$ \psi = \Psi - \psi_\infty, $$

由于 $\nabla_x \psi_\infty = 0$，代入后可得扰动系统：

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u = -\nabla p + \mathrm{div}\, \tau, \quad \mathrm{div}\, u = 0, \\ \partial_t \psi + (u \cdot \nabla_x) \psi = \nu \Delta_R \psi + \mathrm{div}_R\left[ -(\nabla u) R (\psi + \psi_\infty) \right] + L \psi, \\ \tau = \int_B R \otimes (\nabla_R U) \psi\, dR, \\ u(x,0) = u_0(x),\quad \psi(x,R,0) = \psi_0(x,R), \\ (\nabla_R \psi + \nabla_R U \psi) \cdot n = 0 \quad \text{on } \partial B. \end{cases} \tag{2} $$

值得注意的是，此系统中没有流体粘性项（$\mu = 0$），因此速度方程本质上是非线性欧拉型方程，这使得分析极为困难。事实上，若不考虑与聚合物的耦合，单纯的 $N$ 维欧拉方程在 $u=0$ 附近并不稳定——例如，Kiselev 和 Šverák [8] 曾证明二维情形下涡量梯度可双指数增长。

然而，在本文中，正是由于聚合物分布 $\psi$ 与流体速度 $u$ 的强耦合作用，系统展现出一种新的波结构（wave structure），这种结构来源于平衡态 $\psi_\infty$ 的对称性和非线性耦合机制，它提供了隐含的耗散效应，从而能够控制速度的增长并实现全局稳定性。

此外，文中所用的函数空间具有特殊加权结构：

  • 在 $R$ 变量上使用加权 $L^2$ 空间：$L^2(dR / \psi_\infty)$，定义范数为
    $$ \|f\|_{L^2}^2 = \iint_{\mathbb{R}^N \times B} |f(x,R)|^2 \frac{dR dx}{\psi_\infty(R)}, $$
  • 引入加权 $\dot{H}^1$ 范数：
    $$ \|f\|_{\dot{H}^1}^2 = \iint_{\mathbb{R}^N \times B} \left| \nabla_R \left( \frac{f}{\psi_\infty} \right) \right|^2 \psi_\infty(R)\, dR dx, $$
    这种空间结构源于熵方法和 Hardy 型不等式，对于处理 $U(R)$ 在 $|R|\to 1$ 处的奇异性至关重要。

特别地，当 $k > 1$ 时，关键常数如

$$ c_1(k) = -2k \int_B \frac{R_1^3 \partial_{R_1} \psi_\infty}{1 - |R|^2} dR, \quad c_2(k) = -2k \int_B \frac{R_1^2 R_2 \partial_{R_2} \psi_\infty}{1 - |R|^2} dR $$

均为正且有限，并满足 $c_1(k) = 3c_2(k) > 0$，这一性质将在揭示波结构中起决定性作用。

综上所述，本文的研究对象是一个高维、无粘性、强非线性、微宏观耦合、带奇异势能和有限域约束的偏微分方程组，其数学挑战极高，尤其体现在缺乏标准耗散机制的情况下如何建立全局稳定性与衰减估计。

2. 主要定理或主要结果

Theorem 1.1（小初值全局存在性与稳定性）

假设 $k > 1$，$s > N/2 + 1$，存在一个小常数 $\varepsilon > 0$，使得对任意初值 $u_0 \in H^s(\mathbb{R}^N)$，$\mathrm{div}\, u_0 = 0$，$\psi_0 \in H^s(\mathbb{R}^N; L^2)$，若满足

$$ \|u_0\|_s + \|\psi_0\|_{s,L^2} \leq \varepsilon, $$

则系统 (2) 存在唯一的全局解 $(u,\psi)$。进一步地，对所有 $t > 0$，有如下能量估计成立：

$$ \|u(t)\|_s^2 + \|\psi(t)\|_{s,L^2}^2 + \int_0^t \left( \|\nabla u(s)\|_{s-1}^2 + \nu \|\nabla \psi(s)\|_{s,L^2}^2 + \|\psi(s)\|_{s,\dot{H}^1}^2 \right) ds \leq C \varepsilon^2. $$

注：这里的能量泛函定义为

$$ E_1(t) := \sup_{0 \leq s \leq t} \left( \|u(s)\|_s^2 + \|\psi(s)\|_{s,L^2}^2 \right) + \int_0^t \left( \nu \|\nabla \psi(s)\|_{s,L^2}^2 + \|\psi(s)\|_{s,\dot{H}^1}^2 \right) ds, $$

和

$$ E_2(t) := \int_0^t \|\nabla u(s)\|_{s-1}^2 ds, $$

分别衡量解的大小和速度梯度的时间可积性。

Theorem 1.2（大时间衰减估计）

在 Theorem 1.1 的条件下，若进一步要求初值满足

$$ \|u_0\|_{L^1} + \|\psi_0\|_{L^1(L^2)} \leq \varepsilon, $$

则解满足以下衰减估计：

$$ \|u(t)\|_{1} \leq C \varepsilon (1+t)^{-N/4}, \tag{4} $$

$$ \|\psi(t)\|_{1,L^2} \leq C \varepsilon (1+t)^{-N/4 + 1/2}. $$

注：此处 $\|\cdot\|_1$ 指 $H^1(\mathbb{R}^N)$ 范数，$\|\cdot\|_{1,L^2}$ 指 $H^1_x(L^2_R)$ 范数。

Remark 关键观察

  • 尽管系统中没有显式的速度耗散项（即无 $\Delta u$），但通过波结构诱导的隐式耗散机制，速度仍表现出与热方程相同的衰减速率 $(1+t)^{-N/4}$。
  • 该速率与具有完整粘性的 FENE 模型一致，说明聚合物耦合本身可以替代粘性提供有效的耗散。
  • 这是首次对无粘 Euler-FENE 模型建立严格的全局稳定性与最优衰减结果。

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

a) 核心分析工具

（1）Sobolev 空间框架与高阶能量估计

作者采用标准的 $H^s$ 高阶能量方法进行估计。但由于缺少速度耗散，传统能量估计无法闭合。为此，引入两个能量泛函：

  • $E_1(t)$：控制解本身的 $H^s$ 范数及其关于 $\psi$ 的耗散；
  • $E_2(t)$：专门用于控制 $\nabla u$ 的时间积分，这是稳定性的关键。

目标是证明：

$$ E_1(t) \lesssim E_1(0) + C E_1(t)^{3/2} + C E_2(t)^{3/2}, $$

$$ E_2(t) \lesssim C E_1(0) + C E_1(t) + C E_1(t)^{3/2} + C E_2(t)^{3/2}, $$

然后通过连续性论证完成闭合。

（2）波结构分析（Wave Structure Analysis）

这是本文最核心的技术创新。通过对扰动系统进行细致分析，发现速度场 $u$ 满足某种阻尼波动型方程。具体来说，定义变形率张量 $Du = \frac{1}{2}(\nabla u + (\nabla u)^T)$，从第二方程积分 $R_j R_k$ 后可得：

$$ C(k)[Du]_{j,k} = \frac{d}{dt} \int \psi R_j R_k dR + \cdots, $$

其中左边出现了 $Du$，右边包含了 $\psi$ 的时间导数。利用此关系，可以将 $\|\nabla u\|$ 的时间积分转化为对 $\psi$ 及其导数的估计。

更进一步，应用 Leray 投影 $P = I - \nabla \Delta^{-1} \mathrm{div}$ 消去压力，并计算 $\partial_t P \mathrm{div}\, \tau$，最终导出如下二阶演化系统：

$$ \partial_t^2 u - \nu \partial_t \Delta u - c_2(k) \Delta u + P \mathrm{div} \int_B L\psi \, R \otimes \nabla_R U \, dR = N_3, $$

$$ \partial_t^2 \Upsilon - \nu \partial_t \Delta \Upsilon - c_2(k) \Delta \Upsilon - \partial_t P \mathrm{div} \int_B L\psi \, R \otimes \nabla_R U \, dR = N_4, $$

其中 $\Upsilon = P \mathrm{div}\, \tau$，$N_3, N_4$ 为非线性项。

这个形式揭示了 $u$ 所满足的阻尼波动方程结构，其中 $-\nu \partial_t \Delta u$ 提供了类似热传导的耗散，$-c_2(k)\Delta u$ 类似于恢复力项。尽管原系统是欧拉型，但由于耦合，速度获得了“记忆效应”和有效耗散。

（3）Fourier 分裂法（Fourier Splitting Method）

为了获得大时间衰减率，作者结合波结构与 Fourier 分析技巧。将频率空间分为低频区 $|\xi| \leq g(t)$ 与高频区 $|\xi| > g(t)$，其中 $g(t)$ 是适当选取的时间依赖函数（如 $g(t)^2 \sim (1+t)^{-1}$）。

利用低频部分的能量估计和 Sobolev 嵌入，再结合波结构提供的额外正则性，最终导出：

$$ \frac{d}{dt} \left( \|\hat{u}(\xi)\|^2 + \cdots \right) + c |\xi|^2 \|\hat{u}(\xi)\|^2 \lesssim \cdots, $$

在低频区域积分即可得到 $(1+t)^{-N/4}$ 衰减。

（4）Hardy 型不等式与 Poincaré 不等式

处理 $U(R)$ 在边界 $|R|=1$ 处的奇异性是关键技术难点。作者引用 Masmoudi [14] 提出的加权 Hardy 不等式：

Lemma 2.2（Hardy 型不等式）
若 $\psi \in L^2 \cap \dot{H}^1$ 且 $\int_B \psi\, dR = 0$，则

$$ \left( \int_B \frac{|\psi|}{1 - |R|} dR \right)^2 \lesssim \int_B \left| \nabla_R \left( \frac{\psi}{\psi_\infty} \right) \right|^2 \psi_\infty\, dR. $$

若 $k > 1$，还有

$$ \int_B \frac{|\psi|^2}{(1 - |R|^2)^2} dR \lesssim \|\psi\|_{\dot{H}^1}^2. $$

这些不等式允许控制涉及 $\nabla_R U = \frac{2k R}{1 - |R|^2}$ 的奇异项，是闭合高阶估计的基础。

b) 技术难点与应对策略

（1）无速度耗散导致的标准能量框架失效

标准欧拉方程不具备耗散，难以控制 $\|\nabla u\|$ 的增长。本文通过挖掘耦合产生的波结构，发现 $\|\nabla u\|$ 的时间积分可通过 $\psi$ 的耗散项控制。例如，在估计 $E_2(t)$ 时，利用：

$$ \|Du\|_{s-1}^2 \lesssim \sum_{j,k} \left| \int \psi R_j R_k dR \right| + \cdots, $$

进而将其与 $\|\psi\|_{s,L^2}, \|\nabla \psi\|_{s,L^2}$ 关联起来。

（2）非线性项的精细分解

对于非线性项如 $\nabla^s(u \cdot \nabla u)$，使用经典Coifman-Meyer型估计：

$$ \|\nabla^s(fg) - f \nabla^s g\|_{L^2} \lesssim \|\nabla f\|_{L^\infty} \|\nabla^{s-1} g\|_{L^2} + \|g\|_{L^\infty} \|\nabla^s f\|_{L^2}, $$

结合 Sobolev 嵌入 $H^s \hookrightarrow L^\infty$（当 $s > N/2$）来控制。

（3）构造两层估计以处理时间导数项

在衰减估计中出现形如

$$ \int \partial_t(\psi R \otimes R) : Du\, dR dx $$

的交叉项，含有时间导数。作者创新性地引入 $\dot{H}^1$ 层估计，并构造复合能量：

$$ \mathcal{E}(t) = \|u\|_1^2 + \|\psi\|_{1,L^2}^2 + \eta(t) \int (\psi R \otimes R) : Du\, dR dx, $$

其中 $\eta(t)$ 为小权重函数，从而吸收交叉项的影响。

4. 与之前工作的比较

a) 突破与改进

（1）首次在无粘情况下建立 Euler-FENE 模型的全局稳定性

以往关于 FENE 模型的研究大多假设有完整的速度耗散（即 Navier-Stokes-FENE）。例如：

  • Masmoudi [14,15] 研究了有粘情形下的局部/全局适定性；
  • Luo-Yin [11,12] 和 Luo-Zhang [13] 建立了有粘情况下的 $L^2$ 衰减。

本文的突破在于完全去除粘性项（$\mu = 0$），仅保留聚合物扩散 $\nu > 0$，却依然得到了全局解和最优衰减，表明聚合物耦合本身就足以稳定流体运动。

（2）揭示并利用“波结构”作为新型稳定机制

此前工作未曾系统分析过此类耦合系统中的“波结构”。本文首次明确指出，平衡态 $\psi_\infty$ 的各向同性与非线性耦合共同激发了一个类波动的传播机制，它在数学上表现为速度满足带有记忆和阻尼的二阶方程。

这一发现不仅解决了稳定性问题，也为理解其他无耗散耦合系统（如 Vlasov-Euler、kinetic-fluid 系统）提供了新思路。

（3）**推广到任意维度 $N\geq 2$ **

多数早期结果集中在 $N=2$ 或 $N=3$，而本文适用于所有 $N \geq 2$，且衰减速率 $(1+t)^{-N/4}$ 符合标准色散估计规律，体现出方法的普适性。

（4）最优衰减速率匹配粘性模型

令人惊讶的是，即使没有粘性，速度衰减速率仍达到 $(1+t)^{-N/4}$，与热方程相同，且与有粘 FENE 模型一致。这说明聚合物引起的隐式耗散在长时间尺度下等效于粘性耗散，是对“粘性非必需”的深刻验证。

（5）降低对初始数据的要求

相比某些需要高正则性或小 $L^1$ 初值的工作，本文只需 $H^s$ 小初值（$s > N/2 + 1$）即可建立全局存在性；而在衰减部分才额外要求 $L^1$ 性质，结构清晰合理。

综上所述，本文在无速度耗散的极端条件下，通过深入挖掘耦合系统内在的波结构，成功建立了 Euler-FENE 模型的全局稳定性与最优衰减，实现了从“不稳定欧拉方程”到“稳定耦合系统”的跨越，是流体-聚合物耦合理论的一项重要进展。

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

1. Unique continuation and stabilization for nonlinear Schrödinger equations under the Geometric Control Condition

在几何控制条件下非线性薛定谔方程的唯一延拓性与稳定性

作者: Cristóbal Loyola

类别: math.AP, 35A20, 35B60, 93B05, 93B07, 35Q55

发布日期: 2025-10-16

arXiv ID: 2510.14632v1

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

摘要翻译

中文标题: 在几何控制条件下非线性薛定谔方程的唯一延拓性与稳定性

摘要翻译: 本文证明了在几何控制条件成立的区域 $\omega$ 上，具有解析非线性的半线性薛定谔方程的解在有限时间内整体解析性的传播。我们的方法改进了Laurent与作者近期提出的一种技术，该技术结合了控制论方法与Galerkin逼近，以从可观测性成立的区域出发实现解析性的时域传播。作为主要推论，我们得到了当解在区域 $\omega$ 上消失时，定义在二维和三维紧流形上的次临界半线性薛定谔方程的唯一延拓性。此外，仅需在观测区域上满足几何控制条件，即可实现半全局控制与稳定化。特别地，这一结果肯定地回答了Dehman、Gérard和Lebeau在2006年提出的关于非线性情形下的一个公开问题。

2. Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry

基于岸线数据重构浮标处的非线性波及其在分段斜坡地形下海啸反演问题中的应用

作者: Oleksandr Bobrovnikov, Madison Jones, Shriya Prasanna, Josiah Smith, Alexei Rybkin, Efim Pelinovsky

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14177v1

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

摘要翻译

中文标题: 基于岸线数据重构浮标处的非线性波及其在分段斜坡地形下海啸反演问题中的应用

摘要翻译: 本文讨论如下一类反问题：给定海啸波的爬升（run-up）数据，我们能否恢复其初始形态？我们在非线性浅水方程的框架下研究该问题，该模型被广泛用于海啸传播与淹没过程的研究。此前研究表明，在无限倾斜海底地形的情形下，可以从岸线观测数据中恢复初始水位位移和速度 \cite{Rybkin23,Rybkin24,Rybkin25}。本文考虑有限倾斜的海底地形情形，证明了可以从岸线数据中恢复位于虚拟浮标处的边界条件（即水位位移和速度）。进一步地，针对更为复杂的分段斜坡地形，我们探讨将浅水方程与Boussinesq方程进行拼接的方法，以在包含色散效应的模型中恢复初始条件。

3. Exponential and algebraic decay in Euler–alignment system with nonlocal interaction forces

具有非局部相互作用力的欧拉—对齐系统中的指数与代数衰减

作者: José A. Carrillo, Young-Pil Choi, Dowan Koo, Oliver Tse

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.14123v1

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

摘要翻译

中文标题: 具有非局部相互作用力的欧拉—对齐系统中的指数与代数衰减

摘要翻译: 本文研究无压强欧拉系统在非局部速度对齐及相互作用力下的长时间行为，旨在刻画经典解在一般相互作用势 $W$ 和通信权重下的渐近收敛性。我们在三种情形下建立了定量收敛结果。在一维情形中，若势函数为 $(\lambda,\Lambda)$-凸（即满足一致上下二次界），则有界通信权重导致指数衰减，而弱奇异通信权重导致精确的代数衰减速率。对于库仑—二次型势函数 $W(x) = -|x| + \frac{1}{2}|x|^2$，我们证明了在有界通信权重下具有指数收敛性，而在奇异通信权重下具有代数上界估计。在多维情形中，若势函数为一致 $(\lambda,\Lambda)$-凸，则有界权重仍导致指数衰减，奇异权重则产生改进的代数衰减。在所有情况下，密度（至多相差一个平移）收敛于相互作用能量的极小化子，而速度趋于一个均匀常数。一个统一的特征是：收敛速率仅依赖于通信权重的局部性质——有界核导致指数收敛，而弱奇异核导致代数收敛速率。因此，我们的结果为具有一般相互作用势的欧拉—对齐动力学的渐近行为提供了全面的描述。

不相关论文（基本信息）

1. Vortex lines interaction in the three-dimensional magnetic Ginzburg–Landau model

三维磁性Ginzburg-Landau模型中涡旋线的相互作用

作者: Carlos Román, Etienne Sandier, Sylvia Serfaty

类别: math.AP, 35Q56, 82D55, 35J50, 49K10

发布日期: 2025-10-16

arXiv ID: 2510.14910v1

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

摘要: We complete our study of the three dimensional Ginzburg–Landau functional with magnetic field, in the asymptotic regime of a small inverse Ginzburg–Landau parameter $\varepsilon$, and near the first critical field $H_{c_1}$ for which the first vortex filaments appear in energy minimizers. Under a nondegeneracy condition, we show a next order asymptotic expansion of $H_{c_1}$ as $\varepsilon \to 0$, and exhibit a sequence of transitions, with vortex lines appearing one by one as the intensity of the applied magnetic field is increased: passing $H_{c_1}$ there is one vortex, then increasing $H_{c_1}$ by an increment of order $\log |\log\varepsilon|$ a second vortex line appears, etc. These vortex lines accumulate near a special curve $\Gamma_0$, solution to an isoflux problem. We derive a next order energy that the vortex lines must minimize in the asymptotic limit, after a suitable horizontal blow-up around $\Gamma_0$. This energy is the sum of terms where penalizations of the length of the lines, logarithmic repulsion between the lines and magnetic confinement near $\Gamma_0$ compete. This elucidates the shape of vortex lines in superconductors.

2. Adapting Noise to Data: Generative Flows from 1D Processes

适应数据的噪声：源自一维过程的生成流

作者: Jannis Chemseddine, Gregor Kornhardt, Richard Duong, Gabriele Steidl

类别: stat.ML, cs.LG, math.AP

发布日期: 2025-10-14

arXiv ID: 2510.12636v2

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

摘要: We introduce a general framework for constructing generative models using one-dimensional noising processes. Beyond diffusion processes, we outline examples that demonstrate the flexibility of our approach. Motivated by this, we propose a novel framework in which the 1D processes themselves are learnable, achieved by parameterizing the noise distribution through quantile functions that adapt to the data. Our construction integrates seamlessly with standard objectives, including Flow Matching and consistency models. Learning quantile-based noise naturally captures heavy tails and compact supports when present. Numerical experiments highlight both the flexibility and the effectiveness of our method.

3. Stable Type I blow-up for the one-dimensional wave equation with time-derivative nonlinearity

一类具时间导数非线性的二维波动方程的稳定I型爆破

作者: Oliver Gough

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14815v1

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

摘要: We study finite-time blow-up for the one-dimensional nonlinear wave equation with a quadratic time-derivative nonlinearity,

$$ u_{tt}-u_{xx}=(u_t)^2,\qquad (x,t)\in\mathbb R\times[0,T). $$

Building on the work of Ghoul, Liu, and Masmoudi \cite{ghoul2025blow} on the spatial-derivative analogue, we establish the non-existence of smooth, exact self-similar blow-up profiles. Instead we construct an explicit family of \emph{generalised self-similar} solutions, bifurcating from the ODE blow-up, that are smooth within the past light cone and exhibit type-I blow-up at a prescribed point $(x_0,T)$. We further prove asymptotic stability of these profiles under small perturbations in the energy topology. In particular, these profiles verify that the spatially homogeneous ODE blow-up is not asymptotically stable.

4. The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system

化学趋向性与警戒趋向性对捕食者-猎物系统全局存在性与稳定性的同时影响

作者: Gnanasekaran Shanmugasundaram, Jitraj Saha, Rafael Díaz Fuentes

类别: math.AP, 35A01, 35A09, 35B40, 92C17, 37N25, 65N06

发布日期: 2025-10-16

arXiv ID: 2510.14728v1

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

摘要: This study examines a fully parabolic predator-prey chemo-alarm-taxis system under homogeneous Neumann boundary conditions in a bounded domain $\Omega \subset \mathbb{R}^n$ with a smooth boundary $\partial\Omega$. Under specific parameter conditions, it is shown that the system admits a unique, globally bounded classical solution. The convergence of the solution is established through the construction of an appropriate Lyapunov functional. In addition, numerical simulations are presented to validate the asymptotic behaviour of the solution. The results highlight the significant role of chemotaxis and alarm-taxis coefficients in determining the existence and stability of predator-prey models, as discussed in the literature.

5. Stochastic-Periodic Homogenization of a Class of Minimization Problems in Orlicz-Sobolev Spaces

Orlicz-Sobolev空间中一类极小化问题的随机周期均匀化

作者: Dongho Joseph, Fotso Tachago Joel, Tchinda Takougoum Franck

类别: math.AP, math.PR, 35B27, 35B40, 49J55, 37A05, 46E30

发布日期: 2023-11-13

arXiv ID: 2311.10103v4

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

摘要: We develop the stochastic two-scale convergence method in the framework of Orlicz-Sobolev spaces, in order to deal with the homogenization of coupled stochastic-periodic problems in such spaces. One fundamental in this topic is the extension of compactness results for this method to the Orlicz setting. For the application, we show that the sequence of minimizers of a class of highly oscillatory minimizations problems involving integral functionals with convex and nonstandard growth integrands, converges to the minimizer of a homogenized problem.

6. Generalized comparison principle for contact Hamilton-Jacobi equations

接触哈密顿-雅可比方程的广义比较原理

作者: Gengyu Liu, Jianlu Zhang

类别: math.DS, math.AP, 35B40, 35B51, 35F21, 37J50, 37J55, 49L25

发布日期: 2025-09-22

arXiv ID: 2509.17310v2

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

摘要: In this paper, we discuss all the possible pairs $(u,c)\in C(M,\mathbb R)\times\mathbb R$ solving (in the sense of viscosity) the contact Hamilton-Jacobi equation

$$ H (x, d_xu, u) = c,\quad x\in M $$

of which $M$ is a closed manifold and the continuous Hamiltonian $H: (x,p,u)\in T^*M\times\mathbb R\rightarrow\mathbb R$ is convex, coercive in $p$ but merely non-decreasing in $u$. Firstly, we propose a comparison principle for solutions by using the dynamical information of Mather measures. We then describe the structure of $\mathfrak C$ containing all the $c\in\mathbb R$ makes previous equation solvable. We also propose examples to verify the optimality of our approach.

7. Travelling waves for Maxwell’s equations in nonlinear and symmetric media

非线性对称介质中麦克斯韦方程组的行波解

作者: Jarosław Mederski, Jacopo Schino

类别: math.AP

发布日期: 2024-06-03

arXiv ID: 2406.01433v3

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

摘要: We look for travelling wave fields

$$ E(x,y,z,t)= U(x,y) \cos(kz+\omega t)+ \widetilde U(x,y)\sin(kz+\omega t),\quad (x,y,z)\in\mathbb{R}^3,\, t\in\mathbb{R}, $$

satisfying Maxwell’s equations in a nonlinear and cylindrically symmetric medium. We obtain a sequence of solutions with diverging energy that is different from that obtained by McLeod, Stuart, and Troy. In addition, we consider a more general nonlinearity, controlled by an \textit{N}-function.

8. Maxwell’s equations with mixed impedance boundary conditions

具混合阻抗边界条件的麦克斯韦方程组

作者: Ben Schweizer, David Wiedemann

类别: math.AP, 35Q61, 78A25, 35A01

发布日期: 2025-10-16

arXiv ID: 2510.14600v1

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

摘要: We study the time-harmonic Maxwell equations on bounded Lipschitz domains with an impedance boundary condition. The impedance coefficient can be matrix valued such that, in particular, a polarization dependent impedance is modeled. We derive a Fredholm alternative for this system. As a consequence, we obtain the existence of weak solutions for arbitrary sources when the frequency is not a resonance frequency. Our analysis covers the case of singular impedance coefficients.

9. Nonlocal critical problems with mixed boundary conditions and nearly resonant perturbations

具有混合边界条件和近共振扰动的非局部临界问题

作者: Eduardo Colorado, Giovanni Monica Bisci, Alejandro Ortega, Luca Vilasi

类别: math.AP, Primary: 35R11, 35A15, 35S15, 49J35, Secondary: 35J61, 35B33, 58E05

发布日期: 2025-09-29

arXiv ID: 2509.25581v2

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

摘要: We consider the following nonlocal critical problem with mixed Dirichlet-Neumann boundary conditions, \begin{equation} \left{ \begin{array}{ll} (-\Delta)^su=\lambda u+|u|^{2_s^*-2}u &\text{in}\ \Omega,\ \mkern+38.5mu u=0& \text{on}\ \Sigma_{\mathcal{D}},\ \mkern+24mu \displaystyle \frac{\partial u}{\partial \nu}=0 &\text{on}\ \Sigma_{\mathcal{N}}, \end{array} \right. \end{equation} where $(-\Delta)^s$, $s\in (1/2,1)$, is the spectral fractional Laplacian operator, $\Omega\subset\mathbb{R}^N$, $N>2s$, is a smooth bounded domain, $2_s^*=\frac{2N}{N-2s}$ denotes the critical fractional Sobolev exponent, $\lambda>0$ is a real parameter, $\nu$ is the outwards normal to $\partial\Omega$, $\Sigma_{\mathcal{D}}$, $\Sigma_{\mathcal{N}}$ are smooth $(N-1)$–dimensional submanifolds of $\partial\Omega$ such that $\Sigma_{\mathcal{D}}\cup\Sigma_{\mathcal{N}}=\partial\Omega$, $\Sigma_{\mathcal{D}}\cap\Sigma_{\mathcal{N}}=\emptyset$ and $\Sigma_{\mathcal{D}}\cap\overline{\Sigma}_{\mathcal{N}}=\Gamma$ is a smooth $(N-2)$–dimensional submanifold of $\partial\Omega$. By employing a $\nabla$-theorem we prove the existence of multiple solutions when the parameter $\lambda$ is in a left neighborhood of a given eigenvalue of $(-\Delta)^s$.

10. Families of surfaces with constant ratio of principal curvatures and Plateau’s problem

主曲率比为常数的曲面族与普拉托问题

作者: Mikhail Skopenkov, Khusrav Yorov

类别: math.DG, math.AP, 53A10, 53C42, 35J66, 35J96

发布日期: 2025-10-16

arXiv ID: 2510.14527v1

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

摘要: This work is on surfaces with a constant ratio of principal curvatures. These CRPC surfaces generalize minimal surfaces but are much more challenging to construct. We propose a construction of a family of such surfaces containing a given minimal surface without flat points. This leads to a partial solution of Plateau’s problem for CRPC surfaces. We obtain analogous results in isotropic geometry. This work illustrates a general approach to solving Euclidean problems by starting with their isotropic analogs. Besides, we apply the method of successive approximations and analytic majorization.

11. An $L^\infty$-variational problem involving the Fractional Laplacian

一个涉及分数阶拉普拉斯算子的 $L^\infty$-变分问题

作者: Simone Carano, Roger Moser

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14476v1

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

摘要: For $s\in(0,1)$ and an open bounded set $\Omega\subset\mathbb R^n$, we prove existence and uniqueness of absolute minimisers of the supremal functional

$$ E_\infty(u)=\|(-\Delta)^s u\|_{L^\infty(\mathbb R^n)}, $$

where $(-\Delta)^s$ is the Fractional Laplacian of order $s$ and $u$ has prescribed Dirichlet data in the complement of $\Omega$. We further show that the minimiser $u_\infty$ satisfies the (fractional) PDE

$$ (-\Delta)^s u_\infty=E_\infty(u_\infty)\,\mathrm{sgn}f_\infty \qquad\mbox{in }\Omega, $$

for some analytic function $f_\infty\in L^1(\Omega)$ obtained as the restriction of an $s$-harmonic measure $\mu$ in $\Omega$.

12. Small-time approximate controllability of the logarithmic Schr'‘dinger equation

对数薛定谔方程的小时间近似能控性研究

作者: Karine Beauchard, Rémi Carles, Eugenio Pozzoli

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

发布日期: 2025-10-16

arXiv ID: 2510.14461v1

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

摘要: We consider Schr{"o}dinger equations with logarithmic nonlinearity and bilinear controls, posed on $\mathbb{T}^d$ or $\mathbb{R}^d$. We prove their small-time global $L^2$-approximate controllability. The proof consists in extending to this nonlinear framework the approach introduced by the first and third authors in \cite{beauchard-pozzoli2} to control the linear equation: it combines the small-time controllability of phases and gradient flows. Due to the nonlinearity, the required estimates are more difficult to establish than in the linear case. The proof here is inspired by WKB analysis. This is the first result of (small-time) global approximate controllability, for nonlinear Schr{"o}dinger equations, with bilinear controls.

13. On some Elliptic and Parabolic Problems Involving the Anisotropic $p(u)$-Laplacian

关于涉及各向异性 $p(u)$-Laplace 算子的一些椭圆型与抛物型问题

作者: Kaushik Bal, Shilpa Gupta

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14432v1

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

摘要: We investigate a class of elliptic and parabolic partial differential equations driven by p(u) laplacian. This dependence necessitates the use of variable exponent Sobolev spaces specifically tailored to the anisotropic framework. For the elliptic case, we establish the existence of a weak solution by employing the theory of pseudomonotone operators in conjunction with suitable approximation techniques. In the parabolic setting, the existence of a weak solution is obtained via a time discretization scheme and Schauder fixed-point theorem, supported by a priori estimates and compactness arguments.

14. Maximal regularity of Dirichlet problem for the Laplacian in Lipschitz domains

拉普拉斯算子在利普希茨域中狄利克雷问题的最大正则性

作者: Chérif Amrouche, Mohand Moussaoui

类别: math.AP

发布日期: 2025-09-10

arXiv ID: 2509.08543v2

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

摘要: The focus of this work is on the homogeneous and non-homogeneous Dirichlet problem for the Laplacian in bounded Lipschitz domains (BLD). Although it has been extensively studied by many authors, we would like to return to a number of fundamental questions and known results, such as the traces and the maximal regularity of solutions. First, to treat non-homogeneous boundary conditions, we rigorously define the notion of traces for non regular functions. This approach replaces the non-tangential trace notion that has dominated the literature since the 1980s. We identify a functional space E = {v\in H^{1/2}(\Omega);\nabla v\in [H^1/2(\Omega)]’} for which the trace operator is continuous from $E$ into $L^2(\Gamma)$. Second, we address the regularity of solutions to the Laplace equation with homogeneous Dirichlet conditions. Using specific equivalent norms in fractional Sobolev spaces and Grisvard’s results for polygons and polyhedral domains, we prove that maximal regularity $H^{3/2}$ holds in any BLD $\Omega$, for all right-hand sides in the dual of $H^{1/2}_{00}(\Omega)$. This conclusion contradicts the prevailing claims in the literature since the 1990s. Third, we describe some criteria which establish new uniqueness results for harmonic functions in Lipschitz domains. In particular, we show that if $u\in H^{1/2}(\Omega)$ or $u\in W^{1, 2N/(N+1)}(\Omega)$, is harmonic in $\Omega$ and vanishes on $\Gamma$, then $u= 0$. These criteria play a central role in deriving regularity properties. Finally, we revisit the classical Area Integral Estimate. Using Grisvard’s work and an explicit function given by Necas, we show that this inequality cannot hold in its stated form. Since this estimate has been widely used to argue that $H^{3/2}$-regularity is unattainable for data in the dual of $H^{1/2}_{00}(\Omega)$, our counterexample provides a decisive clarification.

15. Parabolic PDEs on a fixed domain with evolving subdomains: function spaces and well-posedness

固定区域上具有演化子区域的抛物型偏微分方程：函数空间与适定性

作者: Van Chien Le, Karel Van Bockstal

类别: math.AP, 35A15, 35K15, 35R05

发布日期: 2025-10-16

arXiv ID: 2510.14373v1

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

摘要: This paper develops the necessary ingredients for the variational approach of initial boundary-value problems of parabolic partial differential equations on a fixed spatial domain containing evolving subdomains. In particular, we introduce function spaces for the variational solution that extend standard Sobolev-Bochner spaces to account for a coefficient associated with the time derivative that may be discontinuous across the evolving interface. We further show the density of smooth functions in these spaces by extending the mollification technique and the Reynolds transport theorem, and establish the corresponding “embedding” theory and an integration by parts formula. Finally, we prove the well-posedness of the space-time variational formulation in the natural setting using the Banach-Necas-Babuska theorem.

16. Viscosity solutions posed on star-shaped network with Kirchhoff’s boundary condition: Well-posedness

具有基尔霍夫边界条件的星形网络上的粘性解：适定性

作者: Isaac Ohavi

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14364v1

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

摘要: The aim of this work is to establish the well-posedness of fully nonlinear partial differential equations (PDE) posed on a star-shaped network, having nonlinear Kirchhoff’s boundary condition at the vertex, and possibly degenerate. We obtain a comparison theorem, for discontinuous viscosity solutions, following the recent ideas obtained by Ohavi for second order problems, building test functions at the vertex solutions of Eikonal equations with well-designed coefficients. Another strong result obtained in this contribution is to show that any generalized Kirchhoff’s viscosity solution introduced by Lions-Souganidis, is indeed a Kirchhoff’s viscosity solution. In other terms, the values of the Hamiltonians are not required at the vertex in the analysis of these types of PDE systems.

17. Sobolev regularity for the perturbed fractional 1-Laplace equations in the subquadratic case

次二次情形下扰动分数阶1-拉普拉斯方程的Sobolev正则性

作者: Dingding Li, Chao Zhang

类别: math.AP

发布日期: 2025-10-16

arXiv ID: 2510.14346v1

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

摘要: This work investigates the Sobolev regularity of solutions to perturbed fractional 1-Laplace equations. Under the assumption that weak solutions are locally bounded, we establish that the regularity properties are analogous to those observed in the superquadratic case. By introducing the threshold $\frac{p-1}{p}$, we divide the range of the parameter $s_p$ into two distinct scenarios. Specifically, for any $s_p\in \left(0, \frac{p-1}{p}\right]$ and $q\ge p$, we demonstrate that the solutions possess $W_{\rm loc}^{\gamma, q}$-regularity for all $\gamma\in \left(0, \frac{s_p p}{p-1}\right)$ and the $W_{\rm loc}^{1, q}$-regularity for any $s_p\in \left(\frac{p-1}{p}, 1\right)$ and $q\ge p$, respectively. Our analysis relies on the nonlocal finite-difference quotient method combined with a Moser-type iteration scheme, which provides a systematic approach to the regularity theory for such nonlocal and singular problems.

18. Propagation speed of traveling waves for diffusive Lotka-Volterra system with strong competition

强竞争下扩散型Lotka-Volterra系统行波解的传播速度

作者: Ken-Ichi Nakamura, Toshiko Ogiwara

类别: math.AP, q-bio.PE

发布日期: 2025-10-16

arXiv ID: 2510.14311v1

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

摘要: We study the propagation speed of bistable traveling waves in the classical two-component diffusive Lotka-Volterra system under strong competition. From an ecological perspective, the sign of the propagation speed determines the long-term outcome of competition between two species and thus plays a central role in predicting the success or failure of invasion of an alien species into habitats occupied by a native species. Using comparison arguments, we establish sufficient conditions determining the sign of the propagation speed, which refine previously known results. In particular, we show that in the symmetric case, where the two species differ only in their diffusion rates, the faster diffuser prevails over a substantially broader parameter range than previously established. Moreover, we demonstrate that when the interspecific competition coefficients differ significantly, the outcome of competition cannot be reversed by adjusting diffusion or growth rates. These findings provide a rigorous theoretical framework for analyzing invasion dynamics, offering sharper mathematical criteria for invasion success or failure.

19. Spectral comparison and splitting theorems for the infinity-Bakry-Emery Ricci curvature

关于无穷Bakry-Emery Ricci曲率的谱比较与分裂定理

作者: Jia-Yong Wu

类别: math.DG, math.AP

发布日期: 2025-09-27

arXiv ID: 2509.23182v2

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

摘要: In this paper, we prove the diameter comparison, the global weighted volume comparison and the splitting theorem in weighted manifolds when the infinity-Bakry-Emery Ricci curvature has a lower bound in the spectrum sense. Our results extend Antonelli-Xu’s spectral Bonnet-Myers and Bishop-Gromov theorems, and Antonelli-Pozzetta-Xu’s spectral splitting theorem to weighted manifolds. Our results are also some supplements of Chu-Hao’s spectral diameter and global volume comparisons, and Yeung’s spectral splitting theorem in weighted manifolds.

20. Rigidity of Asymptotically Hyperboloidal Initial Data Sets with Vanishing Mass

质量为零的渐近双曲型初始数据集的刚性

作者: Sven Hirsch, Hyun Chul Jang, Yiyue Zhang

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

发布日期: 2024-11-11

arXiv ID: 2411.07357v2

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

摘要: In Special Relativity, massless objects are characterized as either vacuum states or as radiation propagating at the speed of light. This distinction extends to General Relativity for asymptotically flat initial data sets (IDS) $(M^n, g, k)$, where vacuum is represented by slices of Minkowski space, and radiation is modeled by slices of $pp$-wave spacetimes. In contrast, we demonstrate that asymptotically hyperboloidal IDS with zero mass must embed isometrically into Minkowski space, with no possible IDS configurations modeling radiation in this setting. Our result holds under the most general assumptions. The proof relies on precise decay estimates for spinors on level sets of spacetime harmonic functions and works in all dimensions.

21. Multiscaling asymptotic behavior of solutions to random high-order heat equations

随机高阶热方程解的多尺度渐近行为研究

作者: Maha Mosaad A Alghamdi, Nikolai Leonenko, Andriy Olenko

类别: math.PR, math.AP, 60F05, 60H15, 60G60

发布日期: 2025-10-15

arXiv ID: 2510.14153v1

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

摘要: This paper studies high-order partial differential equations with random initial conditions that have both long-memory and cyclic behavior. The cases of random initial conditions with the spectral singularities, both at zero (representing classical long-range dependence) and at non-zero frequencies (representing cyclic long-range dependence), are investigated. Using spectral methods and scaling techniques, it is proved that, after proper rescaling and normalization, the solutions converge to Gaussian random fields. For each type of equation, spectral representations and covariance functions of limit fields are given. For odd-order equations, we apply the kernel averaging of solutions to obtain nonexplosive and nondegenerate limits. It is shown that the different limit fields are determined by the even or odd orders of the equations and by the presence or absence of a spectral singularity at zero. Several numeric examples illustrate the obtained theoretical results.

22. Time-harmonic scattering of plane waves from an infinite periodically inhomogeneous medium

平面波在无限周期非均匀介质中的时谐散射

作者: Guanghui Hu, Andreas Rathsfeld, Jiayi Zhang, Ruming Zhang

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.14070v1

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

摘要: We propose a new radiation condition for an infinite inhomogeneous two-dimensional medium which is periodic in the vertical direction and remains invariant in the horizontal direction. The classical Rayleigh-expansion radiation condition does not apply to our case, because this would require the medium to be inhomogeneous in a half plane. We utilize the Floquet theory to derive upward/downward wave modes and define radiation conditions by expansions w.r.t. these modes. The downward radiation conditions leads to a downward Dirichlet-to-Neumann map which can be used to truncate the infinite inhomogeneous domain in the vertical direction. So we prove mapping properties of the upward/downward Dirichlet-to-Neumann maps based on the asymptotic behavior of high-order wave modes. Finally, we verify the strong ellipticity of the sesquilinear form corresponding to the new scattering problem and show the unique solvability for all wavenumbers with the exception of a countable set of numbers bounded below by a small positive constant.

23. The Fujita exponent for a heat equation with mixed local and nonlocal nonlinearities on the Heisenberg group

Heisenberg群上带有局部与非局部混合非线性项的热方程的Fujita指数

作者: Ahmad Z. Fino, Mokhtar Kirane

类别: math.AP, 35K55, 35B44, 35A01, 26A33

发布日期: 2025-07-22

arXiv ID: 2507.16411v2

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

摘要: This article deals with the problems of local and global solvability for a semilinear heat equation on the Heisenberg group involving a mixed local and nonlocal nonlinearity. The characteristic features of such equations, arising from the interplay between the geometric structure of the Heisenberg group and the combined nonlinearity, are analyzed in detail. The need to distinguish between subcritical and supercritical regimes is identified and justified through rigorous analysis. On the basis of the study, the author suggests precise conditions under which local-in-time mild solutions exist uniquely for regular, nonnegative initial data. It is proved that global existence holds under appropriate growth restrictions on the nonlinear terms. To complement these results, it is shown, by employing the capacity method, that solutions cannot exist globally in time when the nonlinearity exceeds a critical threshold. As a result, the Fujita exponent is formulated and identified as the dividing line between global existence and finite-time blow-up. In addition, lifespan estimates were obtained in the supercritical regime, providing insight into how the size of the initial data influences the time of blow-up.

24. Diameter bounds in 3d Type I Ricci flows

三维I型Ricci流中的直径界

作者: Panagiotis Gianniotis

类别: math.DG, math.AP

发布日期: 2025-10-15

arXiv ID: 2510.14019v1

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

摘要: We prove that a three dimensional compact Ricci flow that encounters a Type I singularity has uniformly bounded diameter up to the singular time, thus giving an affirmative answer - for Type I singularities - to a conjecture of Perelman. To achieve this, we introduce a concept of a neck-region for a Ricci flow, analogous to the neck-regions introduced by Jiang-Naber and Cheeger-Jiang-Naber, in the study of Ricci limit spaces. We then prove that the associated packing measure is, in a certain sense, Ahlfors regular, a result that holds in any dimension.

25. Law of Large Numbers and Central Limit Theorem for random sets of solitons of the focusing nonlinear Schrödinger equation

聚焦非线性薛定谔方程的随机孤子集合的大数定律与中心极限定理

作者: Manuela Girotti, Tamara Grava, Ken D. T-R McLaughlin, Joseph Najnudel

类别: math-ph, math.AP, math.MP, math.PR, nlin.PS, nlin.SI

发布日期: 2024-11-26

arXiv ID: 2411.17036v5

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

摘要: We study a random configuration of $N$ soliton solutions $\psi_N(x,t;\boldsymbol{\lambda})$ of the cubic focusing Nonlinear Schr"odinger (fNLS) equation in one space dimension. The $N$ soliton solutions are parametrized by $2N$ complex numbers $(\boldsymbol{\lambda}, \boldsymbol{c})$ where $\boldsymbol{\lambda}\in\mathbb{C}_+^N$ are the eigenvalues of the Zakharov-Shabat linear operator, and $ \boldsymbol{c}\in\mathbb{C}^N\backslash \{0\}$ are the norming constants of the corresponding eigenfunctions. The randomness is obtained by choosing the complex eigenvalues to be i.i.d. random variables sampled from a probability distribution with compact support in the complex plane. The corresponding norming constants are interpolated by a smooth function of the eigenvalues. Then we consider the expectation of the random measure associated to this random spectral data. Such expectation uniquely identifies, via the Zakharov-Shabat inverse spectral problem, a solution $\psi_\infty(x,t)$ of the fNLS equation. This solution can be interpreted as a soliton gas solution. We prove a Law of Large Numbers and a Central Limit Theorem for the differences $\psi_N(x,t;\boldsymbol{\lambda})-\psi_\infty(x,t)$ and $|\psi_N(x,t;\boldsymbol{\lambda})|^2-|\psi_\infty(x,t)|^2$ when $(x,t)$ are in a compact set of $\mathbb R\times\mathbb R^+$; we additionally compute the correlation functions.