生成时间: 2026年02月12日 04:39:38

重点关注论文数量: 4

了解领域论文数量: 1

不相关论文数量: 18

重点关注论文(完整分析)

1. 多维可压缩Navier-Stokes–Korteweg方程组柯西问题在大初值下的整体强解

Global strong solutions with large initial data for the Cauchy problem of the multi-dimensional compressible Navier-Stokes-Korteweg system

作者: Xiangdi Huang, Yongteng Gu, Muxi Lei

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10700v1

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

详细分析

中文标题: 多维可压缩Navier-Stokes–Korteweg方程组柯西问题在大初值下的整体强解


研究对象和背景

本文研究的是多维可压缩Navier-Stokes–Korteweg(NSK)系统,即带内聚力(capillarity)效应的可压缩粘性流体模型。该系统最初由Dunn与Serrin于1985年基于热力学扩展与梯度理论严格推导而出,亦被称为量子Navier-Stokes方程——因其在半经典极限下可导出量子流体力学中的Bohm势项,广泛应用于超流、冷原子气体、半导体载流子输运及相变界面动力学等物理场景。

系统定义在全空间 $\mathbb{R}^N$($N=2$ 或 $3$)上,描述密度 $\rho = \rho(x,t) > 0$ 与速度 $u = u(x,t) \in \mathbb{R}^N$ 的演化,其控制方程为:

$$ \begin{cases} \partial_t \rho + \nabla \cdot (\rho u) = 0, \\ \partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla P(\rho) = \nabla \cdot S + \nabla \cdot K. \end{cases} \tag{1.1} $$

其中压力 $P(\rho)$ 满足等熵状态方程 $P(\rho) = \rho^\gamma$,$\gamma \ge 1$ 为绝热指数;粘性应力张量 $S$ 由Newton本构律给出:

$$ S = 2\mu(\rho) D(u) + \lambda(\rho) (\nabla \cdot u) I, \quad D(u) = \frac{1}{2}(\nabla u + \nabla u^\top), \tag{1.2} $$

且粘性系数满足物理约束 $\mu(\rho) > 0$、$2\mu(\rho) + N \lambda(\rho) \ge 0$;为利用BD熵结构并保证系统稳定性,进一步假设其满足BD关系:

$$ \lambda(\rho) = 2\big(\rho \mu'(\rho) - \mu(\rho)\big). $$

而Korteweg张量 $K$ 刻画了密度梯度诱导的内聚力效应,其散度形式为:

$$ \nabla \cdot K = \nabla \left( \rho \kappa(\rho) \Delta \rho + \frac{\kappa(\rho) + \rho \kappa'(\rho)}{2} |\nabla \rho|^2 \right) - \nabla \cdot \big( \kappa(\rho) \nabla \rho \otimes \nabla \rho \big), \tag{1.3} $$

其中 $\kappa(\rho) > 0$ 是毛细系数。本文聚焦于具有物理意义的浅水型粘性+反平方毛细情形:

$$ \mu(\rho) = \rho,\quad \lambda(\rho) = 0,\quad \kappa(\rho) = \frac{1}{\rho}, $$

此时Korteweg项简化为 $\nabla \cdot K = \nabla \cdot (\rho \nabla \nabla \log \rho)$,整个系统重写为:

$$ \begin{cases} \partial_t \rho + \nabla \cdot (\rho u) = 0, \\ \partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla (\rho^\gamma) = \nabla \cdot (2\rho D(u)) + \nabla \cdot (\rho \nabla \nabla \log \rho). \end{cases} \tag{1.5} $$

该系统配备柯西初值条件:

$$ \rho(x,0) = \rho_0(x), \quad u(x,0) = u_0(x), \quad x \in \mathbb{R}^N, \tag{1.6} $$

并要求非真空远场条件:

$$ \rho(x,t) \to \bar{\rho} > 0, \quad u(x,t) \to 0, \quad \text{as } |x| \to \infty. \tag{1.7} $$

值得注意的是,本文不施加任何对称性(如径向、球对称)、周期性或小初值假设;初始数据允许任意大的 $L^\infty$-范数,仅需满足:

$$ 0 < \underline{\varrho} \le \rho_0(x) \le \overline{\varrho} < \infty, \quad (\rho_0 - \bar{\rho}) \in H^3(\mathbb{R}^N), \quad u_0 \in H^2(\mathbb{R}^N). \tag{1.9} $$

函数空间设定极具挑战性:解需落在非齐次Sobolev-Besov混合空间中,例如 $\rho - \bar{\rho} \in C([0,T]; H^3) \cap L^2(0,T; H^4)$,$u \in C([0,T]; H^2) \cap L^2(0,T; H^3)$。关键难点在于:全空间无紧性,远场非零导致标准能量法失效;密度可能出现振荡甚至退化(尽管初值远离真空,但演化中仍可能趋近零);Korteweg项含三阶导数且具强非线性耦合,破坏抛物正则性。

因此,本文的研究区域是无界全空间 $\mathbb{R}^N$,边值条件为自然远场条件(non-vacuum far-field),而非Dirichlet/Neumann边界;函数空间选择强调临界正则性与时空可积性平衡,依赖Littlewood-Paley分解与Chemin-Lerner型范数(如 $L^q_T B^s_{p,r}$),以精细刻画高频耗散与低频传播的协同效应。


主要定理或主要结果

本文的核心成果是以下整体强解存在唯一性定理:

Theorem 1.1.
设 $N = 2$ 或 $N = 3$,绝热指数 $\gamma$ 满足:

$$ \begin{cases} \gamma \in [1,\infty), & \text{若 } N = 2, \\ \gamma \in [1,8/3), & \text{若 } N = 3. \end{cases} \tag{1.8} $$

令 $\bar{\rho} > 0$ 为密度远场常数,初始数据 $(\rho_0,u_0)$ 满足:

$$ 0 < \underline{\varrho} \le \rho_0(x) \le \overline{\varrho} < \infty, \quad (\rho_0 - \bar{\rho}) \in H^3(\mathbb{R}^N), \quad u_0 \in H^2(\mathbb{R}^N). $$

则柯西问题 (1.5)–(1.7) 存在唯一全局强解 $(\rho,u)$,满足对任意 $0 < T < \infty$ 及 $(x,t) \in \mathbb{R}^N \times [0,T]$:

  • 密度有界性:存在依赖于 $T$ 与初值的常数 $C(T) > 0$,使得

    $$ (C(T))^{-1} \le \rho(x,t) \le C(T); \tag{1.10} $$
  • 正则性:

    $$ \begin{aligned} &(\rho - \bar{\rho}) \in C([0,T]; H^3) \cap L^2(0,T; H^4), \\ &\rho_t \in C([0,T]; H^1) \cap L^2(0,T; H^2), \\ &u \in C([0,T]; H^2) \cap L^2(0,T; H^3), \\ &u_t \in L^\infty(0,T; L^2) \cap L^2(0,T; H^1). \end{aligned} \tag{1.11} $$

该定理首次在全空间、无对称性、任意大的非真空初值下,建立了NSK系统的全局强解,突破了此前所有工作的框架限制。


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

a) 核心分析学工具

本文构建了一个多层次、递进式的分析框架,融合现代调和分析与非线性偏微分方程前沿技术:

  • Littlewood-Paley理论与Besov空间分析:
    引入非齐次Littlewood-Paley分解 $\Delta_j$ 与部分和 $S_j$,定义Besov范数:

    $$ \|u\|_{B^s_{p,r}} := \begin{cases} \left( \sum_{j \ge -1} 2^{jsr} \|\Delta_j u\|_{L^p}^r \right)^{1/r}, & r < \infty, \\ \sup_{j \ge -1} 2^{js} \|\Delta_j u\|_{L^p}, & r = \infty. \end{cases} \tag{2.4} $$

    并采用Chemin-Lerner型时空范数:

    $$ \|u\|_{L^q_T B^s_{p,r}} := \left\| \left( 2^{js} \|\Delta_j u\|_{L^q_T L^p} \right)_{j \ge -1} \right\|_{\ell^r}. \tag{2.5} $$

    此工具使作者能分离不同频率尺度的行为,对高阶导数项(如 $\nabla \cdot (\rho \nabla \nabla \log \rho)$)进行精细估计。

  • 抛物型极大正则性(Maximal Regularity):
    针对密度方程的抛物结构(经有效速度变换后),使用如下热方程极大正则性引理:
    若 $u_0 \in B^s_{p,r}$,$f \in L^{q_2}_T B^{s-2+2/q_2}_{p,r}$,则解 $u$ 满足

    $$ \|u\|_{L^{q_1}_T B^{s+2/q_1}_{p,r}} \le C \left( \|u_0\|_{B^s_{p,r}} + \|f\|_{L^{q_2}_T B^{s-2+2/q_2}_{p,r}} \right). \tag{2.19} $$

    该估计是建立密度上界的关键,将 $\rho$ 控制在 $L^\infty_T L^\infty_x$ 中。

  • BD熵与基本能量恒等式:
    构造BD型能量泛函:

    $$ E_0 = \int_{\mathbb{R}^3} \left( \frac{1}{2} \rho_0 |u_0|^2 + \Pi(\rho_0) + 2 |\nabla \sqrt{\rho_0}|^2 \right) dx, \quad \Pi(\rho) = \frac{1}{\gamma} \rho^\gamma - \frac{1}{\gamma} \bar{\rho}^\gamma - \bar{\rho}^{\gamma-1}(\rho - \bar{\rho}), \tag{3.1}–(3.2) $$

    并证明其耗散性:

    $$ \frac{d}{dt} \int \left( \frac{1}{2} \rho |u|^2 + \Pi(\rho) + 2 |\nabla \sqrt{\rho}|^2 \right) dx + 2 \int \rho |D(u)|^2 dx = 0. \tag{3.10} $$

    结合 $\Pi(\rho)$ 的等价性引理(Lemma 3.1–3.2),导出 $\rho - \bar{\rho} \in L^\infty_T L^2_x \cap L^2_T H^1_x$ 等先验控制。

b) 关键技术性技巧

  • 有效速度变换(Effective Velocity):
    定义 $v = u + \nabla \log \rho$,将原系统转化为“准线性抛物-输运”耦合系统:

    $$ \begin{cases} \partial_t \rho + \nabla \cdot (\rho v) - \Delta \rho = 0, \\ \rho \partial_t v + \rho u \cdot \nabla v + \nabla P(\rho) = \nabla \cdot (\rho \nabla v). \end{cases} \tag{1.13} $$

    此变换消除了Korteweg项的三阶导数,将高阶效应“吸收”进新变量 $v$ 的扩散项中,是处理强色散非线性的枢纽。

  • 截断De Giorgi迭代法(Truncated De Giorgi Iteration):
    为克服全空间远场非零导致的经典De Giorgi方法失效,作者构造自适应截断序列:

    $$ \rho^{-1}_{(k_n)} := \max\{ \rho^{-1} - k_n,\, 0 \}, \quad k_n = M(1 - 2^{-n}) + 2 \|\rho_0^{-1}\|_{L^\infty}, $$

    并定义能量泛函:

    $$ U_{n,T} := \|\rho^{-1}_{(k_n)}\|_{L^\infty_T L^2}^2 + \|\nabla \rho^{-1}_{(k_n)}\|_{L^2_T L^2}^2. \tag{1.15} $$

    利用迭代引理(Lemma 2.5)证明 $U_{n,T} \to 0$,从而获得一致正的密度下界 $\rho \ge \underline{\rho}(T) > 0$。

  • 精化的Moser迭代与最优插值不等式:
    为控制 $v$ 的 $L^\infty$ 范数,作者使用最优插值不等式(Lemma 2.3):

    $$ \|u\|_{L^\infty_T B^{s}_{p,1}} \le \frac{C}{s_2 - s_1} \left( \frac{1}{\theta} \|u\|_{L^\infty_T B^{s_1}_{p,\infty}}^\theta \|u\|_{L^\infty_T B^{s_2}_{p,\infty}}^{1-\theta} \right), \tag{2.13} $$

    结合Moser迭代,建立关键不等式:

    $$ \|v\|_{L^\infty} \lesssim \sqrt{ \log\left( e^2 + \|\rho^{-1}\|_{L^\infty} \right) }, $$

    该不等式是连接密度上下界与解正则性的桥梁。


与之前工作的比较

a) 突破性改进

本文实现了三重根本性突破:

  • 首例全空间大初值全局强解:
    此前所有关于NSK系统的全局存在性结果均受限于:① 一维情形(Charve-Haspot);② 周期域(Huang [14]);③ 小初值(Danchin-Desjardins);④ 特殊对称性(Gu-Huang 径向解)。本文是首个在 $\mathbb{R}^2$ 与 $\mathbb{R}^3$ 上,对任意大的非真空初值,建立强解全局存在性的结果,且无需任何几何假设。

  • 彻底解决密度正性问题:
    Haspot [10] 与 Yu-Wu [22] 虽获密度上界,但因无法控制 $\|\rho^{1/(p+2)} v\|_{L^\infty_T L^{p+2}}$ 的 $p$-独立性,未能获得下界;本文通过截断De Giorgi迭代,将下界估计从“依赖 $p$ 的Moser型”升级为“统一、自适应、全时空一致”的正性控制,解决了困扰该领域十余年的核心障碍。

  • 维度与参数范围的实质性拓展:
    对 $N = 3$,此前周期域结果要求 $\gamma < 3$(Li-Xin),而本文将阈值提升至 $\gamma < 8/3$;对临界情形 $\gamma = 1$(等温流),本文首次将其纳入统一框架,指出方法可直接适用(Remark 1.1),而此前工作均回避此情形。

综上,本文不仅提供了技术范式的革新(截断De Giorgi + Chemin-Lerner + 有效速度),更在物理建模意义上确立了NSK系统在真实无界环境中的数学适定性,为后续研究量子流体湍流、界面不稳定性等前沿问题奠定了严格基础。


2. 动态滑移边界条件下非自治Navier-Stokes流的稳定性分析

Stabilization of nonautonomous Navier-Stokes flows under dynamic slip boundary conditions

作者: Buddhika Priyasad, Sérgio S. Rodrigues

类别: math.AP, math.OC

发布日期: 2026-02-11

arXiv ID: 2602.10678v1

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

详细分析

中文标题: 动态滑移边界条件下非自治Navier-Stokes流的稳定性分析


研究对象和背景

本文研究的是非自治(time-dependent)不可压缩Navier–Stokes方程组在动态滑移边界条件(dynamic slip boundary conditions)下的指数稳定化问题。其核心目标是:对任意给定的、可能随时间演化的参考轨迹 $ \tilde{y}(x,t) $,设计一个有限维、空间局域化、显式构造的反馈控制律,使得受控系统解 $ y(x,t) $ 以指数速率收敛于该轨迹,即满足

$$ \left\| (y(\cdot,t),\, y(\cdot,t)|_{\partial\Omega}) - (\tilde{y}(\cdot,t),\, \tilde{y}(\cdot,t)|_{\partial\Omega}) \right\|_X \le C e^{-\mu(t-s)} \left\| (y(\cdot,s),\, y(\cdot,s)|_{\partial\Omega}) - (\tilde{y}(\cdot,s),\, \tilde{y}(\cdot,s)|_{\partial\Omega}) \right\|_X, $$


对所有 $ t \ge s \ge 0 $ 成立,其中 $ X = L^2(\Omega,\mathbb{R}^d) \times L^2(\partial\Omega,T\partial\Omega) $,$ d \in \{2,3\} $。

所研究的原始系统为定义在光滑有界区域 $ \Omega \subset \mathbb{R}^d $ 上的如下方程组:

$$ \begin{cases} \partial_t \tilde{y} = \nu \Delta \tilde{y} - \langle \tilde{y} \cdot \nabla \rangle \tilde{y} + \nabla \tilde{p} + f, & \mathrm{div}\,\tilde{y} = 0 \quad \text{in } \Omega \times \mathbb{R}_+, \\ \beta \partial_t \tilde{y}|_{\partial\Omega} = - D \tilde{y}|_{\partial\Omega} - \nu L \tilde{y} + g, & (\tilde{y} \cdot n)|_{\partial\Omega} = 0 \quad \text{on } \partial\Omega \times \mathbb{R}_+, \\ \tilde{y}(\cdot,0) = \tilde{y}_0, \quad \tilde{y}(\cdot,0)|_{\partial\Omega} = \tilde{y}_{0\partial}, \end{cases} \tag{1.1} $$

其中 $ \tilde{y}(x,t) \in \mathbb{R}^d $ 是流体速度场,$ \tilde{p}(x,t) \in \mathbb{R} $ 是压力;$ \nu > 0 $ 为运动粘性系数;$ \beta > 0 $ 是边界动力学的时间尺度参数;$ n(x) $ 是 $ \partial\Omega $ 上单位外法向量;$ T\partial\Omega $ 表示切丛;$ D: V \to V' $ 是定义在切向函数空间 $ V \coloneqq L^2(\partial\Omega,T\partial\Omega) $ 上的对称正定算子,满足

$$ \alpha \coloneqq \min_{\xi \in V} \frac{\langle D\xi,\xi\rangle_{V',V}}{|\xi|_{L^2(\partial\Omega,T\partial\Omega)}^2} > 0. $$


关键结构在于动态边界条件:第一式中 $ \beta \partial_t \tilde{y}|_{\partial\Omega} $ 项表明边界值本身满足一个演化方程,而非静态约束(如经典 Navier 条件 $ \tilde{y}_\tau = -\kappa L\tilde{y} $)。算子 $ L $ 具有广义 Green 公式

$$ (\Delta w,z)_{L^2(\Omega,\mathbb{R}^d)} = -(Lw,Lz)_{L^2(\Omega,\mathbb{R}^{d_1})} + (Lw,z)_{L^2(\partial\Omega,T\partial\Omega)}, $$


对所有散度为零且切向于边界的光滑向量场 $ w,z $ 成立。文中明确指出,该框架统一涵盖多种物理意义明确的边界类型:

  • 动态 Navier 条件:$ Lw = [(\nabla w)n]_\tau $,即切向应力;
  • 动态涡量型(vorticity-type)条件:当 $ d=3 $ 时 $ Lw = -n \times \mathrm{curl}\,w $;当 $ d=2 $ 时 $ Lw = (\mathrm{curl}\,w) n^\perp $;
  • 动态 Neumann 条件:$ Lw = [(\nabla w)n]_\tau $,但此时 $ L $ 取为全梯度投影。

此外,区域 $ \Omega $ 允许为多连通域(multi-connected domain),即其边界 $ \partial\Omega $ 由有限个互不相交的闭曲面组成(如环形域、带孔圆盘等),这显著增加了谱分析与函数空间构造的复杂性。

函数空间设置亦具特殊性:作者引入扩展 Hilbert 空间

$$ \mathcal{H} \coloneqq L^2(\Omega,\mathbb{R}^d) \times L^2(\partial\Omega,T\partial\Omega), $$


并定义能量空间

$$ \mathcal{V} \coloneqq \left\{ (h,\varphi) \in H^1(\Omega,\mathbb{R}^d) \times L^2(\partial\Omega,T\partial\Omega) \mid h|_{\partial\Omega} = \varphi \right\}, $$


其中 $ H^1 $-范数由 $ \|h\|_{H^1}^2 + \|Lh\|_{L^2}^2 $ 等价刻画(Assumption 3.2)。这一设置将内部状态与边界状态耦合为单一变量 $ Y = (y,\varphi) \in \mathcal{V} $,从而将原系统重写为抽象进化方程,为后续算子理论分析奠定基础。


主要定理或主要结果

论文的核心结论由 Theorem 2.2 给出(后推广为更一般的 Theorem 6.1):

Theorem 2.2(指数稳定化主定理)
设目标轨迹 $ \tilde{y} \in L^\infty(\mathbb{R}_+, V) $ 满足 Assumption 4.1(即 $ \|\nabla \tilde{y}\|_{L^\infty(\mathbb{R}_+;L^2(\Omega))} < \infty $),则存在 $ M_* \in \mathbb{N} $,使得对任意 $ M \ge M_* $,存在 $ \lambda_*(M) \ge 0 $,使得对所有 $ \lambda \ge \lambda_*(M) $,闭环系统

$$ \begin{cases} \partial_t y = \nu \Delta y - \langle y \cdot \nabla \rangle y + \nabla p + f + U_\diamond K_{\lambda,M}(y - \tilde{y}), & \mathrm{div}\,y = 0, \\ \beta \partial_t y|_{\partial\Omega} = -D y|_{\partial\Omega} - \nu L y + g, & (y \cdot n)|_{\partial\Omega} = 0, \\ y(\cdot,0) = y_0, \quad y(\cdot,0)|_{\partial\Omega} = y_{0\partial}, \end{cases} $$


存在唯一弱解 $ y $,并满足指数稳定估计

$$ \left\| (y(\cdot,t),\, y(\cdot,t)|_{\partial\Omega}) - (\tilde{y}(\cdot,t),\, \tilde{y}(\cdot,t)|_{\partial\Omega}) \right\|_{\mathcal{H}} \le e^{-\mu(t-s)} \left\| (y(\cdot,s),\, y(\cdot,s)|_{\partial\Omega}) - (\tilde{y}(\cdot,s),\, \tilde{y}(\cdot,s)|_{\partial\Omega}) \right\|_{\mathcal{H}}, $$


对所有 $ t \ge s \ge 0 $ 成立,其中衰减率 $ \mu = \beta^{-1}\alpha $,常数 $ C = 1 $。

该定理的关键假设包括:目标轨迹有界且梯度 $ L^2 $-有界;控制作用通过 $ m = d \cdot \#P $ 个空间局域化内部执行器实现(每个支撑集为小区域 $ \omega_i \subset \Omega $ 上的指示函数乘以标准基向量);反馈律为显式构造的斜投影型(oblique projection)线性反馈 $ K_{\lambda,M} $,不依赖于系统谱信息。


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

a) 核心分析学工具

论文采用一套融合算子半群理论、变分框架与Galerkin逼近的综合方法。其技术骨架可分解为以下四层:

第一层:抽象演化系统的构造
通过 Helmholtz 投影 $ \Pi $ 将原系统投影至散度为零切向子空间 $ \mathcal{H} $,并定义扩展 Stokes 算子 $ A: \mathcal{D}(A) \subset \mathcal{H} \to \mathcal{H} $,其双线性型为

$$ a\left( \begin{bmatrix} w \\ \varphi \end{bmatrix}, \begin{bmatrix} z \\ \psi \end{bmatrix} \right) \coloneqq \nu (Lw,Lz)_{L^2(\Omega,\mathbb{R}^{d_1})} + (w,z)_{L^2(\Omega,\mathbb{R}^d)} + (D\varphi,\psi)_{L^2(\partial\Omega,T\partial\Omega)}. $$


由 Assumption 3.3(紧嵌入 $ \mathcal{V} \hookrightarrow \mathcal{H} $)及 Lax-Milgram 引理,$ A $ 是自伴、正定、具有紧逆的无界算子,故存在正交完备特征基 $ \{e_j\} $ 与递增特征值 $ \{\gamma_j\} $。此为后续 Galerkin 截断与能量估计提供谱离散基础。

第二层:Galerkin 逼近与弱解存在性
对自由系统(无控)与受控系统,均采用标准 Galerkin 方法:取 $ Y_N(t) = \sum_{j=1}^N c_j^N(t) e_j \in \mathrm{span}\{e_1,\dots,e_N\} $,代入弱形式后获得常微分方程组。利用 Aubin–Lions 紧性引理(因 $ \mathcal{V} \hookrightarrow \mathcal{H} \hookrightarrow \mathcal{V}' $ 三重嵌入)从 $ Y_N $ 序列中抽取强收敛子列,结合 $ B(y,y) $ 的双线性连续性($ |B(y,z)|_{\mathcal{V}'} \le C \|y\|_{\mathcal{V}} \|z\|_{\mathcal{V}} $)完成极限过程,证明弱解存在唯一性(Section 4)。

第三层:Lyapunov 能量泛函与单调性增强
核心技巧在于构造合适的能量泛函。令误差 $ Y = (y,\varphi) = (y-\tilde{y},\, y|_{\partial\Omega}-\tilde{y}|_{\partial\Omega}) $,则其满足

$$ \frac{d}{dt} \begin{bmatrix} y \\ \beta \varphi \end{bmatrix} + A Y + A_{\mathrm{rc}} Y + N(Y) = F_{\lambda,M} Y. $$


关键一步是证明:对足够大的执行器数量 $ M $ 与增益 $ \lambda $,闭环算子 $ A + F_{\lambda,M} $ 在 $ \mathcal{V} $ 上满足强制性(coercivity),即存在 $ \mu > 0 $ 使得

$$ \left( (A + F_{\lambda,M})Y, Y \right)_{\mathcal{H}} \ge \mu \|Y\|_{\mathcal{H}}^2. $$


此性质由 Section 5 中的“局部化算子单调性引理”保证:通过精细设计执行器支撑集 $ \omega_i $ 的几何分布(满足 Assumption 5.1 的 Sobolev 扩展性质),可使 $ \Pi U_m $ 在 $ \mathcal{H} $ 中稠密逼近,从而斜投影 $ P_{\Pi U_m} $ 控制住 $ A $ 的低频部分。

第四层:指数衰减的 Gronwall 型估计
在获得强制性后,对 $ \|Y(t)\|_{\mathcal{H}}^2 $ 求导,利用非线性项 $ N(Y) $ 的局部 Lipschitz 性($ \|N(Y_1)-N(Y_2)\|_{\mathcal{V}'} \le C_R \|Y_1-Y_2\|_{\mathcal{H}} $,其中 $ R $ 为球半径)及 Young 不等式,导出微分不等式

$$ \frac{d}{dt} \|Y\|_{\mathcal{H}}^2 + 2\mu \|Y\|_{\mathcal{H}}^2 \le C \|Y\|_{\mathcal{H}}^3. $$


再结合局部存在性与先验有界性(由目标轨迹有界性及初始误差控制),应用改进的 Gronwall 不等式(含立方项)得到全局指数衰减(Section 6)。

b) 非线性与奇异性处理技巧

  • 非线性对流项 $ B(y,y) $ 的处理:采用经典 $ H^1 $-$ L^2 $-配对估计,利用 $ \|B(y,z)\|_{\mathcal{V}'} \le C \|y\|_{\mathcal{V}} \|z\|_{\mathcal{H}} $($ d=2,3 $),避免高阶导数爆炸。
  • 动态边界耦合的奇异性:通过扩展空间 $ \mathcal{H} $ 将 $ y $ 与 $ y|_{\partial\Omega} $ 视为整体变量,使 $ \beta \partial_t y|_{\partial\Omega} $ 项自然融入时间导数,消除了传统方法中需单独处理边界时间导数的困难。
  • 无谱假设的突破:摒弃传统基于线性化算子谱分解的控制器设计(要求稳态解、简单特征值等),转而利用斜投影 $ \Pi U_m $ 的几何逼近能力,仅需执行器数量 $ M $ 足够大即可,适用于任意时间依赖轨迹。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重实质性突破:

第一,边界条件类型的本质拓展
此前关于 Navier–Stokes 稳定化的工作(如 [56]、[7–9])几乎全部局限于静态边界条件($ \beta = 0 $),例如经典 Dirichlet、Navier 或 curl 型条件。而本文首次系统处理动态滑移边界($ \beta > 0 $),其数学难点在于:边界值自身满足 ODE,导致系统成为微分-代数混合系统(DAE),传统半群理论失效。作者通过构造扩展空间 $ \mathcal{H} $ 与扩展算子 $ A $,将 DAE 转化为标准抽象 Cauchy 问题,为动态边界情形建立全新分析范式。

第二,目标轨迹的普适性提升
经典结果(如 [10,14,35])多针对稳态解(steady state) 或周期解,依赖线性化算子的谱隙(spectral gap)性质。本文则处理任意 $ L^\infty(\mathbb{R}_+;V) $ 类时间依赖轨迹,完全摆脱谱分析束缚。这源于其反馈律不依赖于任何线性化算子的特征结构,而仅需目标轨迹的梯度有界性(Assumption 4.1),极大拓宽了适用场景(如追踪湍流瞬态、工程启停过程等)。

第三,控制器构造的实用性强化
既有有限维控制器(如 [20,24])常要求执行器覆盖整个区域或满足特定谱条件(如 $ \Phi_i $ 为特征函数)。本文控制器仅需 $ m = d \cdot \#P $ 个空间局域化、形状规则(矩形/三角形支撑)的执行器,且增益 $ \lambda $ 可显式选取($ \lambda \ge \lambda_*(M) $),无需在线辨识或迭代优化。更重要的是,其稳定性率 $ \mu = \beta^{-1}\alpha $ 显式依赖于边界动力学参数,揭示了“越快的边界响应($ \beta \to 0 $),越快的系统收敛”这一物理直观,且与静态极限 $ \beta = 0 $ 下任意高速率可实现的结果一致([56]),验证了理论自洽性。

综上,本文不仅填补了动态边界 Navier–Stokes 稳定化研究的空白,更提供了一套可迁移至其他非自治偏微分方程(如 Boussinesq、MHD)的通用框架。


3. 旋转三维Boussinesq方程在剪切流(Couette流)与层结温度背景下的稳定性阈值

Stability threshold for 3D Boussinesq equations with rotation near the Couette flow and stratified temperature

作者: Wenting Huang, Zekai Luo, Ying Sun, Xiaojing Xu

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10591v1

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

详细分析

中文标题: 旋转三维Boussinesq方程在剪切流(Couette流)与层结温度背景下的稳定性阈值


研究对象和背景

本文研究的是三维旋转Boussinesq方程组在高雷诺数 $\mathbf{Re} \gg 1$ 下围绕Couette流 $(y,0,0)^\top$ 与垂直层结温度场 $\Theta_s(z) = 1 + \alpha^2 z$ 的非线性稳定性阈值问题。该系统建模了地球大气与海洋中同时受科里奥利力(旋转效应) 和 浮力(层结热效应) 驱动的不可压流体运动,是地转流体力学的核心模型之一。

论文所研究的原始方程组为如下定义在区域

$$ \Omega = \mathbb{T} \times \mathbb{R} \times \mathbb{T} $$


上的三维Boussinesq方程(含旋转):

$$ \begin{cases} \partial_t \mathbf{V} + \mathbf{V} \cdot \nabla \mathbf{V} - \nu \Delta \mathbf{V} + \nabla P + \beta \mathbf{e}_3 \times \mathbf{V} = \Theta \, \mathbf{e}_3, \\ \partial_t \Theta + \mathbf{V} \cdot \nabla \Theta - \mu \Delta \Theta = 0, \\ \nabla \cdot \mathbf{V} = 0, \end{cases} $$

其中:

  • $\mathbf{V} = (V_1,V_2,V_3)^\top$ 是速度场;
  • $P$ 是压力;
  • $\Theta$ 是温度(或位温扰动);
  • $\mathbf{e}_3 = (0,0,1)^\top$ 是竖直单位向量;
  • $\beta \in \mathbb{R}\setminus\{0\}$ 表征旋转强度(科里奥利参数),对应角速度为 $\beta/2$;
  • $\nu = \mathbf{Re}^{-1} > 0$ 是运动粘性系数;
  • $\mu$ 是热扩散率;为简化,取 $\mu = \nu$, $g = 1$;
  • $\alpha > 0$ 是Brunt–Väisälä频率,刻画层结稳定性:$\alpha^2 = -\frac{1}{\rho_0}\frac{d\rho}{dz} > 0$,对应稳定层结(暖流体在上、冷流体在下),重力提供恢复力,支持内重力波传播。

稳态解(基流)为:

$$ \mathbf{V}_s(y) = (y,0,0)^\top, \quad \Theta_s(z) = 1 + \alpha^2 z, $$

其对应压力为

$$ P_s(y,z) = z + \frac{1}{2}(\alpha^2 z^2 - \beta y^2). $$

为研究扰动稳定性,引入扰动变量:

$$ \mathbf{u} = \mathbf{V} - \mathbf{V}_s, \quad \theta = \frac{1}{\alpha}(\Theta - \Theta_s), \quad p = P - P_s. $$

代入原系统并利用不可压条件 $\nabla \cdot \mathbf{u} = 0$,可得扰动方程(即论文(1.5)式):

$$ \begin{cases} \partial_t \mathbf{u} + y \partial_{x_1} \mathbf{u} - \nu \Delta \mathbf{u} + \nabla p_L + \begin{pmatrix} 0 \\ (1-\beta) u_2 \\ \beta u_1 \end{pmatrix} = -\mathbf{u} \cdot \nabla \mathbf{u} - \nabla p_{NL}, \\ \partial_t \theta + y \partial_{x_1} \theta - \nu \Delta \theta - \alpha u_3 = -\mathbf{u} \cdot \nabla \theta, \\ \nabla \cdot \mathbf{u} = 0, \\ (\mathbf{u},\theta)(0,\mathbf{x}) = (\mathbf{u}_{\mathrm{in}}, \theta_{\mathrm{in}})(\mathbf{x}). \end{cases} $$

此处线性压力 $p_L$ 由椭圆投影确定:

$$ \Delta p_L = -2 \partial_{x_1} u_2 + \beta (\partial_{x_1} u_2 - \partial_y u_1) - \alpha \partial_z \theta, $$

而非线性压力 $p_{NL}$ 满足 $\Delta p_{NL} = -\partial_i u_j \partial_j u_i$。

值得注意的是,该问题定义在无限长通道 $\mathbb{T}_{x_1} \times \mathbb{R}_y \times \mathbb{T}_z$ 上,即:

  • $x_1$ 和 $z$ 方向为周期边界(代表水平方向闭合或大尺度平均);
  • $y$ 方向为无界实轴(模拟剪切流的横向展布,避免边界层干扰);
  • 初始数据属于混合范数空间 $H^{N+2} \cap W^{N+3,1}$(Sobolev + Wiener空间),体现对低正则性但高衰减性初始扰动的容许性。

关键物理参数是Bradshaw–Richardson数:

$$ B_\beta := \beta(\beta - 1), $$

它统一刻画旋转与层结耦合的稳定性本质:当 $B_\beta > 0$(即 $\beta < 0$ 或 $\beta > 1$)时,系统处于动力学稳定区;而 $B_\beta = 0$($\beta = 0$ 或 $1$)对应临界情形,易出现抬升效应(lift-up) 或退化色散。

本文聚焦于 $B_\beta > 0$ 且 $B_\beta \ne \alpha^2$ 的情形——前者保证线性稳定性与色散波存在,后者避免色散关系 $h(\alpha,\beta,\eta,l) = \sqrt{\alpha^2 \eta^2 + B_\beta l^2}$ 在 $\eta = 0$ 处退化为零(即避免“静止波”奇点)。


主要定理或主要结果

论文建立了两个核心定理,分别刻画线性与非线性稳定性。

定理 1.1(线性稳定性)

假设 $\nu > 0$, $\alpha > 0$, $\beta \in \mathbb{R}$ 满足 $B_\beta > 0$,且 $\mathbf{u}_{\mathrm{in}}$ 是散度自由光滑向量场。令 $\mathbf{u} = \mathbf{u}_\ne + \tilde{\mathbf{u}}_0 + \mathbf{u}_0$、$\theta = \theta_\ne + \tilde{\theta}_0 + \theta_0$ 分别为扰动的速度与温度按 Fourier 模式分解(非零 $x_1$-频、单零模、双零模)。则线性化系统满足以下估计:

(1) 非零模的增强耗散与无粘阻尼:对任意 $s \ge 0$,

$$ \left\| \mathbf{U}_\ne^1(t) \right\|_{H^s} \lesssim \max\left\{ |\beta - 1|, \left|\frac{1}{\beta - 1}\right| \right\} e^{-\frac{1}{6}\nu^{1/3}t} \left( \|\mathbf{u}_{\mathrm{in}}\|_{H^{s+3}} + \|\theta_{\mathrm{in}}\|_{H^{s+2}} \right), $$

$$ \langle t \rangle \left\| \mathbf{U}_\ne^2(t) \right\|_{H^s} \lesssim \max\left\{ |\beta - 1|, \left|\frac{1}{\beta - 1}\right| \right\} e^{-\frac{1}{6}\nu^{1/3}t} \left( \|\mathbf{u}_{\mathrm{in}}\|_{H^{s+5}} + \|\theta_{\mathrm{in}}\|_{H^{s+4}} \right), $$

$$ \left\| (\mathbf{U}_\ne^3, \Theta_\ne)(t) \right\|_{H^s} \lesssim e^{-\frac{1}{6}\nu^{1/3}t} \left( \|\mathbf{u}_{\mathrm{in}}^3\|_{H^{s+3}} + \left|\frac{1}{\beta - 1}\right| \|\omega_{\mathrm{in}}^3\|_{H^{s+2}} + \|\theta_{\mathrm{in}}\|_{H^{s+2}} \right). $$

(2) 单零模($k = 0$)的精确色散解构:对 $l \ne 0$,Fourier 变换后有显式表达式,特征根为

$$ \lambda_1 = -\nu(\eta^2 + l^2), \quad \lambda_\pm = -\nu(\eta^2 + l^2) \pm i \frac{h(\alpha,\beta,\eta,l)}{\sqrt{\eta^2 + l^2}}, \quad h = \sqrt{\alpha^2 \eta^2 + B_\beta l^2}. $$

特别地,当 $l = 0$(双零模),解退化为纯扩散+振荡耦合:

$$ u_1^0(t,y) = e^{\nu \partial_y^2 t} u_{\mathrm{in}}^1, \quad u_2^0 = 0, $$

$$ \begin{pmatrix} u_3^0 \\ \theta^0 \end{pmatrix}(t,y) = e^{\nu \partial_y^2 t} \begin{pmatrix} \cos(\alpha t) & -\sin(\alpha t) \\ \sin(\alpha t) & \cos(\alpha t) \end{pmatrix} \begin{pmatrix} u_{\mathrm{in}}^3 \\ \theta_{\mathrm{in}} \end{pmatrix}. $$

(3) 单零模的色散估计:

$$ \left\| (\tilde{u}_2^0, \tilde{u}_3^0)(t) \right\|_{L^\infty(\mathbb{R}_y \times \mathbb{T}_z)} \lesssim e^{-\nu t} t^{-1/3} \left\| \widehat{\mathbf{u}}_{\mathrm{in}}^0 \right\|_{W^{3,1}(\mathbb{R}_\eta \times \mathbb{Z}_l)}. $$

定理 1.3(非线性稳定性)

假设 $\nu \in (0,1)$, $\beta \in \mathbb{R}$ 满足 $B_\beta > \frac{1}{4}$ 且 $B_\beta \ne \alpha^2$;对任意 $s \ge 7$、$N \ge s + 4$,存在 $\delta = \delta(N,s,\alpha,\beta) > 0$,使得若初始扰动满足零均值条件:

$$ \int_{\mathbb{T}^2} u_{\mathrm{in}}^3 \, dx_1 dz = \int_{\mathbb{T}^2} \theta_{\mathrm{in}} \, dx_1 dz = 0, $$

且

$$ \|\mathbf{u}_{\mathrm{in}}\|_{H^{N+2} \cap W^{N+3,1}} + \|\theta_{\mathrm{in}}\|_{H^{N+1} \cap W^{N+3,1}} = \varepsilon < \delta \, \nu^{14/15}, $$

则非线性系统存在唯一全局解 $(\mathbf{u},\theta)$,并满足:

  • 全局能量控制:

    $$ \|(\mathbf{U},\Theta)\|_{L^\infty_t H^s} + \nu^{1/6} \|(\mathbf{U}_\ne,\Theta_\ne)\|_{L^2_t H^s} + \nu^{1/2} \|\nabla_L(\mathbf{U},\Theta)\|_{L^2_t H^s} \lesssim \varepsilon; $$
  • 无粘阻尼与增强耗散:

    $$ \left\| (\mathbf{U}_\ne^1,\mathbf{U}_\ne^3,\Theta_\ne)(t) \right\|_{L^2} + \langle t \rangle \left\| \mathbf{U}_\ne^2(t) \right\|_{L^2} \lesssim b_{\alpha,\beta} \, e^{-\lambda \nu^{1/3} t} \varepsilon, $$

    其中 $b_{\alpha,\beta} = \frac{1}{\alpha} \sqrt{ \frac{B_\beta + 1}{2} \max\left\{ \frac{1}{B_\beta - 1}, \frac{B_\beta - 1}{1} \right\} }$, $\lambda = \frac{1}{162} \sqrt{ \frac{B_\beta - 1}{B_\beta + 1} }$;

  • 非线性色散:

    $$ \|(\tilde{u}_2^0,\tilde{u}_3^0)(t)\|_{W^{2,\infty}} \lesssim (qt)^{-1/3} e^{-\nu t} \varepsilon_0 + q^{-1/3} \nu^{-2/3} \varepsilon^2, $$

    其中 $q = \frac{|B_\beta - \alpha^2|}{\alpha}$,$\varepsilon_0$ 为单零模初值大小。

该定理断言:只要初始扰动大小低于 $\mathbf{Re}^{-14/15}$ 阈值,则系统保持非线性稳定,不发生转捩。


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

a) 核心分析学工具

论文融合了现代流体力学稳定性理论的多项前沿技术,形成一套自洽的“色散–耗散–非线性校正”三重框架。

1. 坐标系变换与拉格朗日描述
采用经典移动坐标系 $(\tilde{x}_1, y, z) = (x_1 - t y, y, z)$,将剪切输运项 $\mathbf{V}_s \cdot \nabla$ 转化为常系数导数 $\partial_t - t \partial_y$。此变换使 Fourier 符号变为时变函数:

$$ p_L(t,k,\eta,l) = k^2 + (\eta - k t)^2 + l^2, \quad \partial_t p_L = -2k(\eta - k t), $$

从而自然引出时间依赖的频率调制,为后续构造乘子提供基础。

2. Fourier 乘子法与“好未知量”构造
为克服线性耦合项(如 $\beta u_1$, $\alpha \theta$, $\partial_z \theta$)导致的能量估计闭合困难,作者构造三类“好未知量”:

$$ Q_\ne \coloneqq (-|\nabla_L|^{3/2} |\nabla_{L,h}|^{-1}) u_\ne^3, \quad K_\ne \coloneqq (-i \sqrt{\beta - 1} \, |\nabla_L|^{1/2} |\nabla_{L,h}|^{-1}) \omega_\ne^3, \quad H_\ne \coloneqq (-|\nabla_L|^{1/2}) \theta_\ne, $$

其中 $\nabla_L = (\partial_x, \partial_y - t \partial_x, \partial_z)$, $\nabla_{L,h} = (\partial_x, \partial_y - t \partial_x)$。这些未知量通过精心选取的 Fourier 乘子 $C_Q, C_K, C_H$ 实现:

  • 抵消配对耦合项(如 $\beta u_1$ 与 $\partial_z \theta$);
  • 抑制由压力引起的线性拉伸项 $\partial_x \partial_y \Delta^{-1}$ 的增长;
  • 提取最优耗散率 $\nu^{1/3}$(源于 $\nu \Delta$ 与色散时间尺度 $t \sim \nu^{-1/3}$ 的共振)。

3. 增强耗散与无粘阻尼的联合估计
借助加权能量泛函:

$$ \mathcal{E}_s(t) = \|Q_\ne\|_{H^s}^2 + \|K_\ne\|_{H^s}^2 + \|H_\ne\|_{H^s}^2 + \text{交叉项}, $$

并设计时变乘子 $m(t,k,\eta,l)$ 满足

$$ \partial_t m + \nu p_L m \gtrsim \nu^{1/3} m, $$

从而获得指数衰减 $e^{-c \nu^{1/3} t}$。该策略继承并推广了 Bedrossian–Masmoudi–Vicol 对 Couette 流的“乘子能量法”。

4. Strichartz 型色散估计
对单零模,利用色散算子 $A = B_\beta R_l^3 + \nu(\partial_y^2 - l^2)$(其中 $R_l^3 = -i l (l^2 - \partial_y^2)^{-1/2}$)的半群性质,建立如下 Strichartz 型估计:

$$ \left\| e^{t A} h_l \right\|_{L^2_t L^\infty_y} + \nu^{1/2} |l| \left\| e^{t A} h_l \right\|_{L^1_t L^\infty_y} \lesssim \nu^{-1/3} |l|^{-1/6} \|h_l\|_{L^2_y}. $$

该估计是获得 $t^{-1/3}$ 色散衰减的关键。

5. 准线性化与正则性提升技巧
针对 $u_1^0$ 与 $\theta^0$ 缺乏色散的问题,采用准线性化:将非线性项 $\mathbf{u} \cdot \nabla u_1$、$\mathbf{u} \cdot \nabla \theta$ 视为外力,利用已证的 $u_2^0, u_3^0$ 的色散性控制其 $L^1_t$ 范数,再通过 Duhamel 公式与 Sobolev 嵌入获得 $H^s$ 正则性提升。

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

  • 退化色散的消除:在无旋情形($\beta = 0$)下,垂直层结导致色散关系在 $\eta = 0$ 处奇异(分母为零);本文证明:旋转引入 $B_\beta l^2$ 项,使 $h(\eta,l) > 0$ 对所有 $(\eta,l) \ne (0,0)$ 成立,彻底消除奇点。

  • 浮力诱导的竖直正则性损失:$\alpha u_3$ 项在 $z$ 方向不平滑,作者通过将 $\theta$ 与 $u_3$ 耦合进同一未知量 $H_\ne$,并利用 $\nabla_L^{1/2}$ 权重平衡竖直导数损失。

  • 非线性耦合的高频控制:对形如 $u_j \partial_j u_i$ 的项,采用频谱分割:低频部分用 $L^\infty$ 色散控制,高频部分用 $H^N$ 能量控制,并借助 $W^{N+3,1}$ 初值保证 $L^1$ 可积性,规避 Besov 嵌入缺陷。


与之前工作的比较

a) 突破与改进

本文在多个维度上实现了对既有文献的实质性超越:

1. 首次建立旋转+垂直层结下的最优稳定性阈值
此前工作多限于单一效应:

  • Bedrossian–Germain–Masmoudi(2015)对无旋 Navier–Stokes 给出 $\gamma \le 3/2$;
  • Coti Zelati–Del Zotto(2020)对水平层结 Boussinesq 得到 $\gamma \le 11/12$;
  • Li–Sun–Wang–Wei–Zhang(2023)对旋转 Navier–Stokes 达到 $\gamma \le 5/6$。
    本文首次在三维旋转+垂直层结设定下获得 $\gamma = \frac{14}{15} = 0.933\ldots$,逼近理论极限 $1$,且明确指出该阈值由单零模色散与非线性相互作用共同决定。

2. 放宽对初始数据的正则性要求
以往工作(如 [12], [45])依赖 Gevrey 或高阶 Sobolev 范数(如 $H^{100}$);本文仅需 $H^{N+2} \cap W^{N+3,1}$,其中 $N \ge 11$,显著降低对初值光滑性的苛刻要求,更贴近物理实际(如观测误差、数值离散)。

3. 彻底解决“垂直层结奇点”难题
Oxley–Kerswell(2022)实验发现旋转可抑制层结不稳定性,但数学上未解决 $\eta = 0$ 奇点。本文通过旋转–层结协同色散机制 $h = \sqrt{\alpha^2 \eta^2 + B_\beta l^2}$,首次给出严格证明:旋转完全修复垂直层结的色散退化,为地转湍流建模提供理论基石。

4. 统一处理线性耦合与非线性反馈
相比 [16,17] 中需额外假设 $\int u_{\mathrm{in}}^3 = \int \theta_{\mathrm{in}} = 0$,本文在定理 1.3 中仍保留该条件以优化阈值,但在 Remark 1.4 指出:即使取消该条件,仍可得 $\gamma \le 1$,表明结果具有鲁棒性。

综上,本文不仅推进了 Boussinesq 系统稳定性理论的前沿,更构建了一套可推广至更复杂地转模型(如分层旋转 MHD、湿对流)的通用分析范式。


4. 旋转球面上点涡附近的约束性结果

Confinement results near point vortices on the rotating sphere

作者: Martin Donati, Emeric Roulley

类别: math.AP, physics.ao-ph

发布日期: 2026-02-10

arXiv ID: 2602.10061v1

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

详细分析

中文标题: 旋转球面上点涡附近的约束性结果


研究对象和背景

本文研究的是旋转单位球面 $ S^2 \subset \mathbb{R}^3 $ 上的 二维不可压欧拉方程(Euler equation),其物理背景是大尺度大气与海洋动力学中高环量涡旋的演化行为。由于地球自转效应显著,模型需引入恒定角速度 $ \gamma \in \mathbb{R} $ 的均匀旋转,从而构成一个典型的旋转流体系统。

论文考虑的是以绝对涡度 $ \zeta $ 为基本变量的活动标量形式。设 $ x = (x_1,x_2,x_3) \in S^2 $,记 $ |x|_{\mathbb{R}^3} $ 为欧氏范数,$ d\sigma(x) $ 为球面上的标准面积测度,则所研究的方程为:

$$ \partial_t \zeta(t,x) + u(t,x) \cdot \nabla \zeta(t,x) = 0, $$

其中速度场 $ u $ 由Biot–Savart 律给出:

$$ u(t,x) = \frac{1}{2\pi} \int_{S^2} \frac{x \wedge y}{|x - y|_{\mathbb{R}^3}^2} \, \zeta(t,y)\, d\sigma(y) + \gamma \, e_3 \wedge x. $$

该公式体现了两个关键物理成分:第一项为由绝对涡度诱导的无旋速度场(对应于球面上格林函数 $ G(x,y) = \frac{1}{2\pi}\log|x-y|_{\mathbb{R}^3} $ 的梯度旋度),第二项为科里奥利力产生的背景旋转流($ e_3 = (0,0,1) $ 是垂直轴方向单位向量)。值得注意的是,此 Biot–Savart 核在 $ x=y $ 处具有对数奇异性,但因球面紧致性,其积分仍良定义于 $ L^p $ 空间中。

系统还需满足高斯约束条件(Gauss constraint):

$$ \int_{S^2} \zeta(t,x)\, d\sigma(x) = 0, $$

这是由球面紧致性和不可压性 $ \nabla \cdot u = 0 $ 所强制要求的守恒律;它等价于总环量为零,在点涡近似下体现为强度和为零:

$$ \sum_{i=1}^{N} \Gamma_i = 0. $$

此外,边界条件自动满足——因为 $ S^2 $ 是无边闭流形,无需额外指定边界行为;而函数空间选取为 $ \zeta \in L^1 \cap L^\infty(S^2) $,初始数据被假定为在若干点附近高度集中,即满足如下结构化假设(Hypothesis 1.2):

  • 存在常数 $ M > 0 $、$ \varepsilon_0 > 0 $ 及 $ \eta \geq 2 $,使得对所有 $ \varepsilon \in (0,\varepsilon_0) $,
    • $ \|\zeta_0^\varepsilon\|_{L^\infty(S^2)} \leq M \varepsilon^{-\eta} $,
    • $ \zeta_0^\varepsilon = \sum_{i=1}^{N} \zeta_0^{\varepsilon,i} $,其中 $ \Gamma_i = \int_{S^2} \zeta_0^{\varepsilon,i}(x)\,d\sigma(x) $ 满足 $ \sum_i \Gamma_i = 0 $,
    • $ \operatorname{supp}(\zeta_0^{\varepsilon,i}) \subset B(x_i^0,\varepsilon) \subset \mathbb{R}^3 $,即每个分量支撑在半径为 $ \varepsilon $ 的球邻域内。

这种“尖锐集中”初值模拟了现实中强涡旋结构(如极地涡旋、热带气旋),其数学本质是考察点涡动力学作为真实欧拉解的渐近近似是否有效,以及该近似的时间尺度有多长。

研究区域 $ S^2 $ 具有常正曲率,这使其几何性质既不同于平面 $ \mathbb{R}^2 $(零曲率),也区别于一般曲面(变曲率)。特别地,文献指出:在非常曲率情形下,点涡运动的一阶项会包含曲率导数,导致动力学更复杂;而在球面上,常曲率保证了Biot–Savart核仅依赖于弦距 $ |x-y|_{\mathbb{R}^3} $,从而使得分析可沿用平面上发展成熟的涡度约束理论框架,同时必须克服球面拓扑与非欧几里得距离带来的新困难。

综上,本文的研究对象是一个定义在紧致黎曼流形 $ S^2 $ 上、受刚性旋转驱动、具有对数奇异核的非线性输运型偏微分方程;其核心问题是:当初始涡度以 $ \varepsilon $-尺度集中在 $ N $ 个点附近时,解是否长时间保持局部化?能否被点涡轨迹所“捕获”?这一问题直接关联到天气预报模型中涡旋追踪的可靠性与气候系统中涡旋寿命的理论估计。


主要定理或主要结果

本文建立了四个核心定理,分别刻画点涡系统的全局存在性概率、一般情形下的对数约束时间尺度、该尺度的最优性,以及特殊构型下的幂律增强约束。

定理 1.2(点涡碰撞的不可能性)
设 $ N \in \mathbb{N} \setminus \{0,1\} $。则对于几乎所有的初始位置 $ (x_1^0,\dots,x_N^0) \in (S^2)^N $(关于乘积归一化面积测度 $ P_N $),点涡系统

$$ \frac{d}{dt}x_i(t) = \sum_{\substack{j=1 \\ j \ne i}}^{N} \frac{2\Gamma_j}{\pi} \frac{x_i(t) \wedge x_j(t)}{|x_i(t) - x_j(t)|_{\mathbb{R}^3}^2} + \gamma \, e_3 \wedge x_i(t) $$

具有全局解,即不会发生有限时间碰撞。

定理 1.3(对数尺度涡度约束)
设 $ (x_1^0,\dots,x_N^0) \in (S^2)^N $ 两两不同,$ (\Gamma_1,\dots,\Gamma_N) \in \mathbb{R}^N \setminus \{0\} $ 满足 $ \sum_i \Gamma_i = 0 $,且对应点涡解 $ (x_i(t))_{1\le i \le N} $ 全局存在并满足最小距离控制:

$$ \inf_{t \ge 0} \, \min_{i \ne j} |x_i(t) - x_j(t)|_{\mathbb{R}^3} \ge d_0 > 0. $$

则对任意 $ \beta < \frac{1}{2} $,存在 $ \varepsilon_0 = \varepsilon_0(\beta,d_0) > 0 $ 和 $ \alpha = \alpha(\beta,d_0) > 0 $,使得对所有 $ \varepsilon \in (0,\varepsilon_0) $,满足 Hypothesis 1.2 的初值 $ \zeta_0^\varepsilon $ 对应的欧拉解 $ \zeta^\varepsilon(t,x) $ 满足退出时间估计:

$$ \tau_{\varepsilon,\beta} \triangleq \sup \left\{ t \ge 0 : \forall s \in [0,t],\; \operatorname{supp}(\zeta^\varepsilon(s,\cdot)) \subset \bigcup_{i=1}^{N} B(x_i(s),\varepsilon^\beta) \right\} \ge \alpha |\ln \varepsilon|. $$

定理 1.4(对数约束的最优性)
存在某组 $ (x_1^0,\dots,x_N^0) $ 与 $ (\Gamma_1,\dots,\Gamma_N) $,满足 Gauss 条件与最小距离条件,使得对某个 $ \beta_0 < \frac{1}{2} $、$ \eta \ge 2 $,及任意 $ \beta \in (\beta_0,1) $,存在 $ \alpha_0 > 0 $,使得对充分小的 $ \varepsilon > 0 $,存在初值 $ \zeta_0^\varepsilon $ 满足 Hypothesis 1.2,但其退出时间满足:

$$ \tau_{\varepsilon,\beta} \le \alpha_0 |\ln \varepsilon|. $$

定理 1.5(特殊构型下的幂律约束)
存在某组 $ (x_1^0,\dots,x_N^0) $ 与 $ (\Gamma_1,\dots,\Gamma_N) $,满足 Gauss 条件与最小距离条件,使得对任意 $ \beta < \frac{1}{2} $ 与 $ \alpha < \min\{\beta, 2 - 4\beta\} $,存在 $ \varepsilon_0 > 0 $,使得对所有 $ \varepsilon \in (0,\varepsilon_0) $,对应解满足:

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

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

a) 核心分析学工具

本文综合运用了多种现代偏微分方程与动力系统分析技术:

  • Hamiltonian 结构与能量方法:点涡系统 (1.7) 被证明是 Hamilton 系统,其哈密顿量为

    $$ H(x_1,\dots,x_N) = \sum_{i \ne j} \frac{\Gamma_i \Gamma_j}{4\pi} \log|x_i - x_j|_{\mathbb{R}^3} + \gamma e_3 \cdot \sum_i \Gamma_i x_i. $$


    这一结构不仅保障了相空间体积守恒(Liouville 定理),也为构造 Lyapunov 型泛函、控制相对距离提供了基础。

  • 正则化与逼近论证(regularization and approximation):为处理碰撞奇点,作者引入光滑截断对数函数 $ \ln_\varepsilon(r) $,构造正则化哈密顿量 $ H_\varepsilon $,从而获得全局光滑流 $ S_t^\varepsilon $。随后通过单调收敛定理与测度连续性将正则化解极限过渡至原系统,这是证明“几乎必然无碰撞”的关键步骤。

  • 矩估计(moment estimates)与支撑增长控制:定义第 $ p $ 阶矩

    $$ M_p^\varepsilon(t) \triangleq \sum_{i=1}^N \int_{S^2} |x - x_i(t)|_{\mathbb{R}^3}^p \, \zeta^{\varepsilon,i}(t,x)\, d\sigma(x), $$


    并利用输运方程性质与 Biot–Savart 核的奇异性阶数($ \sim |x-y|^{-1} $),推导出

    $$ \frac{d}{dt} M_p^\varepsilon(t) \lesssim M_p^\varepsilon(t)^{1 - \frac{1}{p}} \cdot \left( \sum_{j \ne i} \frac{|\Gamma_j|}{|x_i - x_j|} \right). $$


    结合最小距离假设,得到 $ M_p^\varepsilon(t) \lesssim \varepsilon^{p\beta} e^{Ct} $,进而控制支撑直径增长。

  • 概率测度论与随机过程工具:将初始位置视为随机变量,赋予乘积归一化面积测度 $ P_N $。利用Markov 不等式、停时(stopping time)理论及自然滤波 $ \mathcal{F}_t^\varepsilon $,将碰撞事件测度上界转化为对正则化泛函 $ \phi_\varepsilon(X) \triangleq \sum_{i \ne j} \exp(-\eta \ln_\varepsilon|x_i - x_j|) $ 的期望估计,最终导出 $ P_N(C_N) = 0 $。

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

  • 对数奇异性处理:Biot–Savart 核在 $ x=y $ 处为 $ \sim |x-y|^{-1} $,远强于平面上的 $ \sim |x-y|^{-2} $(后者对应 $ \nabla^\perp \log|x-y| $)。作者巧妙利用球面嵌入 $ \mathbb{R}^3 $ 的几何特性,将 $ |x-y|_{\mathbb{R}^3} $ 视为弦距而非测地距离,从而避免引入更复杂的球面格林函数展开;并通过精细的积分分割(near/far field decomposition)分离奇异与正则部分。

  • 紧致流形上的全局分析:区别于平面情形需处理无穷远行为,球面紧致性自动排除了无限分离或逃逸,简化了能量控制;但同时也丧失平移不变性,迫使所有估计必须依赖于一致几何常数(如曲率、注入半径)。为此,全文大量使用球面三角不等式与 $ |x-y|_{\mathbb{R}^3} \approx \operatorname{dist}_{S^2}(x,y) $ 的双侧等价性(当距离小时)。

  • 多尺度匹配与退出时间精确定义:退出时间 $ \tau_{\varepsilon,\beta} $ 的设定兼顾了物理意义($ \varepsilon^\beta $ 控制“模糊半径”)与数学可行性($ \beta < 1/2 $ 确保矩估计闭合)。证明中采用迭代比较法:先证低阶矩控制,再通过 Gronwall 引理提升至支撑控制,并反复利用点涡解的稳定性假设来压制扰动放大。

  • 线性稳定性与非线性约束的桥接:在 Theorem 1.5 的证明中,作者将点涡构型的线性稳定性(通过 Hessian 矩阵 $ D^2 H $ 的特征值分析)与非线性约束时间联系起来。具体而言,若 $ D^2 H $ 正定,则系统呈现“超稳定”行为,允许构造指数型 Lyapunov 函数,从而将对数时间提升至幂律时间。


与之前工作的比较

a) 突破与改进

本文首次将涡度约束理论系统性地拓展至旋转球面这一具有重要地球物理意义的几何背景下,填补了该领域长期存在的空白。此前所有严格结果均局限于平面 $ \mathbb{R}^2 $、环面 $ \mathbb{T}^2 $ 或有界区域,而球面因其非平凡拓扑与曲率,传统方法面临根本性挑战。

  • 几何推广的实质性突破:以往工作(如 [9,21,22])严重依赖平面上的傅里叶分析、平移对称性及 $ \mathbb{R}^2 $ 中的位势理论。本文完全摒弃这些工具,转而构建基于球面调和分析与嵌入空间几何的新框架。例如,Biot–Savart 律的显式表达(含 $ x \wedge y $ 项)及格林函数 $ \log|x-y|_{\mathbb{R}^3} $ 的选择,均为球面特有,无法简单类比平面情形。

  • 放宽对初始数据的对称性要求:先前关于幂律约束的工作(如 [9,57])往往要求初值接近径向对称或位于特定对称边界(如圆盘中心)。本文 Theorem 1.5 的条件仅依赖于点涡构型本身的线性稳定性,而此类构型(如 Thomson 多边形、极地涡晶)天然存在于球面,无需人为施加初值对称性。

  • 去除对涡强度的额外限制:在平面与环面情形中,“无碰撞”结论常需假设 $ \sum \Gamma_i^2 \neq (\sum \Gamma_i)^2 $ 等非退化条件(见 [26,54])。本文得益于球面紧致性,完全取消该假设,证明对任意满足 Gauss 条件的强度组合,碰撞均为零测事件。

  • 统一性与教学完整性:作者明确指出,本文重写了此前散见于十余篇论文中的证明,提供了一个自洽、统一、 pedagogically transparent 的版本。例如,矩估计、支撑增长、对数约束的三步逻辑链被清晰拆解并强化,使后续研究者得以在此坚实基础上继续推进。

综上,本文不仅是已有理论的“球面类比”,更是通过深刻挖掘几何结构、重构分析工具、打通线性/非线性桥梁,实现了从技术适配到原理创新的跃迁。


了解领域论文(摘要翻译)

1. 多孔介质中的一类流体-固体相互作用问题

A fluid-solid interaction problem in porous media

作者: Diego Alonso-Orán, Rafael Granero-Belinchón

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

发布日期: 2026-02-11

arXiv ID: 2602.10822v1

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

摘要翻译

中文标题: 多孔介质中的一类流体-固体相互作用问题

摘要翻译: 本文针对描述弹性膜下方达西流动的弹性Muskat自由边界问题,推导了相应的渐近界面模型。在界面陡度较小的弱非线性情形下,我们获得了刻画自由边界动力学至二次阶精度的非局部演化方程。在长波薄液膜极限下,我们将运动学条件改写为通量形式,对移动区域进行坐标平坦化处理,并由此导出一类润滑型方程。此外,我们在适当的Wiener函数空间中建立了上述模型的适定性理论。


不相关论文(基本信息)

1. 度量测度空间上局部与非局部 $p$-能量的截断型索博列夫不等式

Cutoff Sobolev inequalities for local and non-local $p$-energies on metric measure spaces

作者: Meng Yang

类别: math.AP, math.FA, math.MG, math.PR

发布日期: 2026-02-11

arXiv ID: 2602.10990v1

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

摘要: For $p>1$, we study subordination phenomena for local and non-local regular $p$-energies on metric measure spaces. Under suitable geometric assumptions, we show that if a local regular $p$-energy satisfies a Poincaré inequality together with a cutoff Sobolev inequality with scaling function $Ψ$, then all associated stable-like non-local $p$-energies with scaling functions strictly below $Ψ$ are regular and satisfy the corresponding non-local cutoff Sobolev inequalities. Moreover, if a stable-like non-local regular $p$-energy with scaling function $Ψ$ satisfies the corresponding non-local cutoff Sobolev inequality, then the same conclusion holds for all associated stable-like non-local $p$-energies with scaling functions below $Ψ$. These results provide a non-linear extension of the classical subordination principle beyond the Dirichlet form framework.


2. 基于模型参考自适应的数据同化方法:面向线性与非线性动力系统

Data assimilation via model reference adaptation for linear and nonlinear dynamical systems

作者: Benedikt Kaltenbach, Christian Aarset, Tram Thi Ngoc Nguyen

类别: math.OC, math.AP, math.DS

发布日期: 2026-02-11

arXiv ID: 2602.10920v1

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

摘要: We address data assimilation for linear and nonlinear dynamical systems via the so-called \emph{model reference adaptive system}. Continuing our theoretical developments in \cite{Tram_Kaltenbacher_2021}, we deliver the first practical implementation of this approach for online parameter identification with time series data. Our semi-implicit scheme couples a modified state equation with a parameter evolution law that is driven by model-data residuals. We demonstrate four benchmark problems of increasing complexity: the Darcy flow, the Fisher-KPP equation, a nonlinear potential equation and finally, an Allen-Cahn type equation. Across all cases, explicit model reference adaptive system construction, verified assumptions and numerically stable reconstructions underline our proposed method as a reliable, versatile tool for data assimilation and real-time inversion.


3. 离散定量等容不等式:涨落估计

Discrete Quantitative Isocapacitary Inequality: Fluctuation Estimates

作者: Marco Cicalese, Leonaard Kreutz, Imteyaz Mansoor

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10918v1

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

摘要: The classical isocapacitary inequality states that, among all sets of fixed volume, the ball uniquely minimizes the capacity. While this result holds in the continuum, it fails in the discrete setting, where the isocapacitary problem may admit multiple minimizers. In this paper we establish quantitative fluctuation estimates for the discrete isocapacitary problem on subsets of $\mathbb{Z}^d$ as their cardinality diverges. Our approach relies on a careful extension of the associated variational problem from the discrete to the continuum setting, combined with sharp (continuum) quantitative isocapacitary inequalities.


4. ε-调和映射的广义能量恒等式与颈部长度

The generalised energy identity and length of necks for $\varepsilon$-harmonic maps

作者: Andrew M. Roberts

类别: math.DG, math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10914v1

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

摘要: In this paper we find analogues for $\varepsilon$-harmonic maps to the generalised energy identity and the existence of geodesic necks result discovered by Yuxiang Li and Youde Wang for $α$-harmonic maps. In particular there exist specific quantities depending only on $\varepsilon$ and the bubbling radius which entirely determine if the full energy identity holds and if a neck forms. In the case these fail we can calculate the energy lost and the length of the geodesic neck based on only these quantities and the biharmonic energy of the bubble.


5. 低能量 $\varepsilon$-调和映射到标准球面

Low energy $\varepsilon$-harmonic maps into the round sphere

作者: Andrew M. Roberts

类别: math.DG, math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10913v1

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

摘要: In this paper we classify the low energy $\varepsilon$-harmonic maps from the surfaces of constant curvature with positive genus into the round sphere. We find that all such maps with degree $\pm1$ are all quantitively close to a bubble configuration with bubbles forming at special points on the domain with bubbling radius proportional to $\varepsilon^{1/4}$.


6. 散度型椭圆方程的内点 $BMO$ 正则性

Interior $BMO$ regularity for elliptic equations in divergence form

作者: Yuanyuan Lian

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10772v1

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

摘要: In this note, we establish the interior $BMO$ regularity of weak solutions to uniformly elliptic equations in divergence form. Moreover, the assumptions on the coefficients are nearly optimal.


7. 含时依赖磁场的薛定谔方程解的波前集

Wave front set of solutions to Schrödinger equations with time-dependent magnetic fields

作者: Fumihito Abe, Ryo Muramatsu

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10741v1

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

摘要: In this paper, we determine the wave front set of solutions to the Schrödinger equation with time-dependent magnetic fields. We considered time-dependent and `not so small’ magnetic fields through the method using the wave packet transform established by K. Kato, M. Kobayashi and S. Ito. Furthermore, we checked that the fundamental solution of the Schrödinger equation in a spatially decaying magnetic field has no singularities as a consequence of our result.


8. 高斯分布间Sinkhorn散度的Wasserstein梯度流

The Wasserstein gradient flow of the Sinkhorn divergence between Gaussian distributions

作者: Mathis Hardion, Théo Lacombe

类别: math.AP, math.OC

发布日期: 2026-02-11

arXiv ID: 2602.10726v1

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

摘要: We study the Wasserstein gradient flow of the Sinkhorn divergence when both the source and the target are Gaussian distributions. We prove the existence of a flow that stays in the class of Gaussian distributions, and is unique in the larger class of measures with strongly-concave and smooth log-densities. We prove that the flow globally converges toward the target measure when the source’s covariance matrix is not singular, and provide counter-examples to global convergence when it is, giving a first answer to an open question raised in [Carlier et al. 2024, \S4.2]. When the covariance matrix of the source distribution commutes with that of the target, we derive more quantitative results that showcase exponential convergence toward the target when the source and the target share their support, but dropping to linear rates (O(t^{-1})) if the target is concentrated on a strict subspace of the source’s support.


9. 各向同性广义伯格斯模型松弛张量的显式表达式及其谱逆变换

Explicit form of relaxation tensor for isotropic extended Burgers model and its spectral inversion

作者: Youjun Deng, Ching-Lung Lin, Gen Nakamura

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10636v1

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

摘要: Concerning the anelastic nature of Earth, the quasi-static extended Burgers model (abbreviated by q-EBM), an integro-differential system, is used to study the free oscillation of Earth (abbreviated by FOE). In this paper, we first provide a general method to obtain an explicit form of the relaxation tensor for inhomogeneous isotropic q-EBM. Then, we apply it to compute the eigenvalues of the free oscillation of Earth, assuming that Earth is a unit ball modeled as a homogeneous and isotropic q-EBM. So far, an analytical and systematic way to compute the eigenvalues of the FOE has been missing when modeling Earth as a q-EBM. In particular, we compute some clusters of eigenvalues (abbreviated by C-ev’s). To be more precise, integrating by parts with respect to time of the q-EBM under the assumption that the initial strain is zero, the q-EBM becomes the sum of two terms. The first term, called the instantaneous term, doesn’t have any integration with respect to time, but the second term, called the memory term, has such an integration. Then, consider the eigenvalues of the instantaneous part of the q-EBM. The eigenfunctions of C-ev’s share the same eigenfunctions of the instantaneous part. However, the C-ev’s may be shifted from the eigenvalues of the instantaneous part. Further, we analyze the structure of C-ev’s and provide an inversion formula identifying the q-EBM from the C-ev’s.


10. 具有退化强制项的障碍问题中自由边界的正则性

Free boundary regularity in obstacle problems with a degenerate forcing term

作者: Yong Liu

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10617v1

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

摘要: In this paper, we consider the properties of a special free boundary point in the following obstacle problem: The Laplacian of u equals f(x) multiplied by the characteristic function of the set where u is positive within the two-dimensional unit ball, where $f(x)=|x|$ is a degenerate forcing term. The key challenge stems from the degeneracy of $f(x)$, which leads to a more complex structure of the free boundary compared to the classical setting. To analyze it, we introduce the epiperimetric inequality developed by Weiss (Invent Math 138:23-50, 1999). Although this powerful tool was firstly introduced for the classical obstacle problem characterized by $f(x)>0$ in B_1, it also proves effective in our degenerate setting. This allows us to first obtain the decay rate of the Weiss energy for all blow-ups at the origin, which in turn implies the uniqueness of the blow-up profiles. With this uniqueness established, we then prove a very weak directional monotonicity properties satisfied by the solutions. This finally yields the regularity of the free boundary at the origin if the origin is a regular point.


11. Heisenberg群上非局部抛物型p-Laplace方程弱解的尾部最优Hölder估计

Optimal in tail Hölder estimates for weak solutions of the nonlocal parabolic p-Laplace equations on the Heisenberg group

作者: Debraj Kar

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10612v1

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

摘要: We prove the Hölder continuity for weak solutions to parabolic p-Laplace equations on the Heisenberg group. We deduce this result while considering an optimal tail condition.


12. 具有非强制性势的非线性狄拉克方程的定态周期解

Stationary periodic solutions to Nonlinear Dirac equations with non-coercive potentials

作者: Fuping Zhang, Ruijun Wu

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10509v1

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

摘要: We obtain periodic solutions for nonlinear Dirac equations with a nonlinear term that is not necessarily coercive. This amounts to study the equation on a three-dimensional torus. The Palais-Smale condition is enhanced by involving a coercive perturbation. Uniform estimates for the critical levels as well as the Sobolev norms for the perturbed solutions are obtained, making it possible to pass to a limit which gives a nontrivial solution.


13. $C^{1,1}$ 开集上的分数阶 Hardy 不等式

Fractional Hardy inequalities on $C^{1,1}$ open sets

作者: Abdelrazek Dieb, Remi Yvant Temgoua

类别: math.AP

发布日期: 2026-02-11

arXiv ID: 2602.10463v1

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

摘要: Let $Ω$ be a bounded open set of class $C^{1,1}$ in $\mathbb{R}^N$ and $s\in(\frac{1}{2}, 1)$. We study a family of fractional Hardy-type inequalities \begin{equation} \frac{c_{N,s}}{2}\displaystyle\iint_{Ω\timesΩ}\frac{(u(x)-u(y))^2}{|x-y|^{N+2s}}\ dxdy-\displaystyleλ\int_Ωu^2\ dx\geq C\displaystyle\int_Ω\frac{u^2}{δ^{2s}}\ dx,~~~\quad\forallλ\in\mathbb{R},~~~~~~~(0.1) \end{equation} with $u\in C_c^\infty(Ω)$ and $C=C(Ω,s,N,λ)>0$. We show that the best constant in $(0.1)$ is achieved if and only if $λ>λ^*(s,Ω)$, for some $λ^*(s,Ω)\in\mathbb{R}$. As a by-product, we derive in particular that the best constant in Hardy inequality $μ_{N,s}(Ω)$ is achieved if and only if $μ_{N,s}(Ω)<\mathfrak{h}_{N,s}$, with $\mathfrak{h}_{N,s}$ being the best constant for the fractional Hardy inequality in the half space. Moreover, if $Ω$ is a convex open set, we obtain a lower bound for $λ^*(s,Ω)$ in terms of the volume of $Ω$. Specifically, we prove that $λ^*(s,Ω)\geq a(N,s)|Ω|^{-\frac{2s}{N}}$ with an explicit constant $a(N,s)>0$. For general bounded $C^{1,1}$ open sets, we prove instead that $λ^*(s,Ω)\geq0$ when $s$ is close to $\frac{1}{2}$. The aforementioned result is proved after showing that $μ_{N,s}(Ω)=\mathfrak{h}_{N,s}$ for $s$ close to $\frac{1}{2}$. In particular, we deduce that, whenever $s$ is sufficiently close to $\frac{1}{2}$, the Hardy constant $μ_{N,s}(Ω)$ is never achieved, hence, behaves differently from that in the local case. This result is completely new in the fractional setting, and was known only for convex open sets for the full range $s\in(\frac{1}{2}, 1)$.


14. 时间尺度上的非线性椭圆型狄利克雷边值问题

Nonlinear elliptic Dirichlet boundary value problems on time scales

作者: Shalmali Bandyopadhyay, F. Ayça Çetinkaya, Tom Cuchta

类别: math.AP, math.SP

发布日期: 2026-02-10

arXiv ID: 2602.10335v1

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

摘要: We establish existence and uniqueness results for nonlinear elliptic Dirichlet boundary value problems on n-dimensional time scale domains. Time scales provide a unified framework that encompasses continuous, discrete, and hybrid settings. Under a Lipschitz condition on the nonlinearity bounded by the first eigenvalue, we prove existence and uniqueness using the contraction mapping theorem. Under a weaker one-sided growth condition, we establish existence using the Leray–Schauder fixed point theorem. To apply these functional analytic methods, we reformulate the problem as an operator equation, which requires developing the spectral theory for the Dirichlet Laplacian with mixed nabla-delta derivatives. We establish self-adjointness, positivity, and completeness of eigenfunctions, and the product eigenfunctions form a complete orthonormal basis in the n-dimensional setting.


15. 具有超高斯解的新对数幂型非线性薛定谔方程

New logarithmic power nonlinear Schrödinger equations with super-Gaussons

作者: Hadi Susanto

类别: nlin.PS, math.AP

发布日期: 2026-02-10

arXiv ID: 2602.10263v1

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

摘要: We introduce a new class of nonlinear Schrödinger equations with a logarithmic-power nonlinearity that admits exact localized solutions of super-Gaussian form. The resulting stationary states possess flat-top profiles with sharp edges and are referred to as super-Gaussons, in analogy with the Gaussian Gaussons of the classical logarithmic NLS (log-NLS). The model, which we call the logarithmic-power NLS (logp-NLS), is parameterized by an exponent $p\geq1$ that controls the degree of flatness of the soliton core and the sharpness of its decay. Mathematically, $p$ interpolates between the standard log-NLS ($p=1$) and increasingly flat-top profiles as $p$ increases, while physically it governs the stiffness of an underlying logarithmic-power compressibility law. The proposed equation is constructed so as to admit super-Gaussian stationary states and can be interpreted a posteriori within a generalized pressure-law framework, thereby extending the log-NLS. We investigate the dynamics of super-Gaussons in one spatial dimension through numerical simulations for various values of $p$, demonstrating how this parameter regulates both the internal structure of the soliton and its collision dynamics. The logp-NLS thus generalizes the standard log-NLS by admitting a broader family of localized states with distinctive structural and dynamical properties, suggesting its relevance for flat-top solitons in nonlinear optics, Bose-Einstein condensates, and related nonlinear media.


16. 离散波导中横向缺陷引起的波散射

Wave scattering by a transversal defect in a discrete waveguide

作者: Elena Medvedeva, Raphael Assier, Anastasia Kisil

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

发布日期: 2026-02-10

arXiv ID: 2602.10227v1

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

摘要: We study wave scattering by a finite transversal strip in a discrete square-lattice waveguide with Dirichlet boundary conditions imposed on the strip and the waveguide walls. The setting is motivated as a discrete analogue of the classical continuous waveguide problem with a screen. The corresponding Wiener–Hopf formulation leads to an equation with a $4 \times 4$ matrix kernel, which reduces to a $2 \times 2$ matrix kernel under some symmetry assumptions. The factorisation prospects of this kernel are discussed, but this route is not followed. Instead, an exact analytical solution is obtained using the pole removal technique. This contrasts with the continuous case, where only approximate solutions are currently available. The reflection and transmission coefficients resulting from an incident duct mode are computed with an accuracy up to $10^{-13}$, showing consistency with theoretical predictions from continuous waveguide theory. In particular, full reflection and zero transmission are recovered as the frequency approaches the cut-off value for the incident mode. Finally, the solution is validated against a numerical computation of the diffraction problem via the Boundary Algebraic Equations method with a tailored lattice Green’s function.


17. 无 $Δ_2$-条件的Orlicz-Sobolev空间中广义多调和算子的特征值问题

An eigenvalue problem for a generalized polyharmonic operator in Orlicz-Sobolev spaces without the $Δ_2$-condition

作者: Ignacio Ceresa Dussel, Julián Fernández Bonder, Pablo Ochoa

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.10077v1

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

摘要: In this paper, we consider a generalized polyharmonic eigenvalue problem of the form $A(u)= λh(u)$ in a bounded smooth domain with Dirichlet boundary conditions in the setting of higher-order Orlicz-Sobolev spaces. Here, $A$ is a very general operator depending on $u$ and arbitrary higher-order derivatives of $u$, whose growth is governed by an Orlicz function, and $h$ is a lower order term. Combining the theories of pseudomonotone operators with complementary systems, we prove that this eigenvalue problem has an infinite number of eigenfunctions and that the corresponding sequence of eigenvalues tends to infinity. We point out that the $Δ_2$-condition is not assumed for the involved Orlicz functions. Finally, we prove a first regularity result for eigenfunctions by following a De Giorgi’s iteration scheme.


18. 一类不满足细致平衡条件的凝并-碎裂方程解向平衡态的收敛性

Convergence to equilibrium for a class of coagulation-fragmentation equations without detailed balance

作者: Apratim Bhattacharya, Sebastian Throm

类别: math.AP

发布日期: 2026-02-10

arXiv ID: 2602.10059v1

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

摘要: We prove convergence to equilibrium for a class of coagulation-fragmentation equations that do not satisfy a detailed balance condition. More precisely, we consider perturbations of constant rate kernels. Our result provides in particular explicit convergence rates. Uniqueness of the stationary states is proven as well for the considered class of kernels.