生成时间: 2026年02月06日 04:29:50

重点关注论文数量: 4

了解领域论文数量: 0

不相关论文数量: 19

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

1. 可压缩磁流体动力学方程在临界Besov空间中的大体积黏度极限

Large bulk viscosity limit for compressible MHD equations in critical Besov spaces

作者: Gennaro Ciampa, Donatella Donatelli, Giada Pellecchia

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05878v1

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

详细分析

中文标题: 可压缩磁流体动力学方程在临界Besov空间中的大体积黏度极限

1. 研究对象和背景

本文研究的是 二维与三维全空间 $\mathbb{R}^d$（$d=2,3$）中可压缩磁流体动力学（Compressible Magnetohydrodynamics, C-MHD）方程组 在 大体积黏度（large bulk viscosity）极限下的渐近行为。该系统描述了导电流体（如太阳风、电离层等离子体）在磁场作用下的耦合运动，其物理建模融合了流体力学与电磁学的基本原理。

论文在引言第1节明确给出了所研究的Cauchy问题——即如下可压缩MHD方程组 $(\text{C-MHD})$：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho v) = 0, \\ \rho v_t + \rho v \cdot \nabla v - \mu \Delta v - (\lambda + \mu)\nabla \mathrm{div}\,v + \nabla P = b \cdot \nabla b - \frac{1}{2} \nabla(|b|^2), \\ b_t + b \,\mathrm{div}\,v + v \cdot \nabla b - b \cdot \nabla v = \eta \Delta b, \\ \mathrm{div}\,b = 0, \\ \rho|_{t=0} = \rho_0,\quad v|_{t=0} = v_0,\quad b|_{t=0} = b_0. \end{cases} $$

其中未知函数为：密度 $\rho = \rho(t,x) > 0$、速度场 $v = v(t,x) \in \mathbb{R}^d$、磁场 $b = b(t,x) \in \mathbb{R}^d$；压力满足绝热律 $P = P(\rho) = A \rho^\gamma$，$\gamma > 1$；黏性系数 $\mu > 0$（剪切黏度），$\lambda$（体积黏度），并定义总黏度 $\nu := \lambda + 2\mu > 0$；$\eta > 0$ 是磁扩散率。初始数据 $(\rho_0,v_0,b_0)$ 满足 $b_0$ 为 散度自由向量场，且 $\rho_0$ 有正下界与上界：$c_0 \le \rho_0 \le c_0^{-1}$。

该系统具有典型的双曲–抛物混合结构：连续性方程为双曲型，动量方程含二阶黏性项（强抛物性），磁场演化方程为对流–扩散型，而约束条件 $\mathrm{div}\,b = 0$ 是由麦克斯韦方程在低频近似下导出的“冻结场”条件。

研究区域为整个欧氏空间 $\mathbb{R}^d$，无边界，因此无需考虑边值条件；但函数空间的选择极为关键：作者采用 临界齐次Besov空间 $ \dot{B}^{s}_{p,q}(\mathbb{R}^d) $，这是由系统的自然缩放不变性所决定的。令密度扰动 $a := \rho - 1$，则系统在变换

$$ a_\kappa(t,x) = a(\kappa^2 t, \kappa x),\quad v_\kappa(t,x) = \kappa v(\kappa^2 t, \kappa x),\quad b_\kappa(t,x) = \kappa b(\kappa^2 t, \kappa x) $$


下保持形式不变。由此导出临界正则性要求：

$$ a \in \dot{B}^{d/p}_{p,1}(\mathbb{R}^d),\quad v,b \in \dot{B}^{d/p-1}_{p,1}(\mathbb{R}^d). $$

特别地，在 $L^2$ 基础上取 $p=2$，得到最常用临界框架：

$$ a \in \dot{B}^{d/2}_{2,1}(\mathbb{R}^d),\quad v,b \in \dot{B}^{d/2-1}_{2,1}(\mathbb{R}^d). $$

这一空间选择具有深刻意义：它恰好处于强解存在性的“门槛”位置——低于此正则性，经典能量法失效；高于此，则可通过Sobolev嵌入获得更高正则性，但会丢失尺度不变性。因此，在临界Besov空间中建立全局适定性，是检验理论深度的核心指标。

本文聚焦于 大体积黏度极限 $\lambda \to \infty$ 下的行为分析。注意：体积黏度 $\lambda$ 不出现在方程的缩放变换中，故不影响临界空间定义，但显著增强对压缩波（声波）的耗散能力。直观上，当 $\lambda \gg 1$ 时，密度扰动被强烈抑制，系统应趋近于不可压缩MHD（I-MHD）：

$$ \begin{cases} V_t + V \cdot \nabla V - \mu \Delta V + \nabla \Pi = B \cdot \nabla B, \\ B_t + V \cdot \nabla B - B \cdot \nabla V = \eta \Delta B, \\ \mathrm{div}\,V = 0,\quad \mathrm{div}\,B = 0, \\ V|_{t=0} = \mathbb{P} v_0,\quad B|_{t=0} = \mathbb{P} b_0, \end{cases} $$

其中 $\mathbb{P}$ 是Leray投影算子，将向量场正交分解为散度自由部分（用 $V,B$ 表示）与梯度部分（用 $Qv,Qb$ 表示）。值得注意的是，由于初始磁场 $b_0$ 已满足 $\mathrm{div}\,b_0 = 0$，故 $B_0 = b_0$，但记号上仍区分大小写以强调不可压缩框架。

因此，本研究的本质是：在临界正则性假设下，证明当 $\lambda \to \infty$ 时，C-MHD 解一致收敛到 I-MHD 解，并给出显式的收敛速率估计；进一步，利用该稳定性结果，构造具有磁重联（magnetic reconnection） 现象的全局光滑解——这是首次在可压缩MHD框架中实现的严格解析构造。

2. 主要定理或主要结果

论文建立了三个核心定理，构成逻辑递进链：从 $L^2$-型临界空间（Theorem A）推广至一般 $L^p$-型临界空间（Theorem B），最终应用于磁重联构造（Theorem C）。

Theorem A（$L^2$-临界情形）

设 $d \ge 2$，初始数据满足：

$$ v_0,b_0 \in \dot{B}^{d/2-1}_{2,1}(\mathbb{R}^d),\quad a_0 := \rho_0 - 1 \in \dot{B}^{d/2}_{2,1}(\mathbb{R}^d),\quad c_0 \le \rho_0 \le c_0^{-1}, $$

且对应的不可压缩初值 $V_0 = \mathbb{P} v_0$, $B_0 = \mathbb{P} b_0$ 生成 I-MHD 的唯一全局强解 $(V,B)$，定义常数： $$ M := |V|{L^\infty(\mathbb{R}+;\dot{B}^{d/2-1}{2,1})} + |V_t,\mu\nabla^2 V|{L^1(\mathbb{R}+;\dot{B}^{d/2-1}{2,1})}

  • |B|{L^\infty(\mathbb{R}+;\dot{B}^{d/2-1}{2,1})} + |B_t,\eta\nabla^2 B|{L^1(\mathbb{R}+;\dot{B}^{d/2-1}{2,1})}. $$

若 $\nu \ge \mu$ 且存在大常数 $C$ 使得

$$ C e^{C(M+M^2)}\left( \|Q v_0\|_{\dot{B}^{d/2-1}_{2,1}} + \|a_0\|_{\dot{B}^{d/2-1}_{2,1}} + \nu \|a_0\|_{\dot{B}^{d/2}_{2,1}} + M^2 + \mu^2 \right) \le \mu \nu, $$

则 C-MHD 存在唯一全局强解 $(\rho,v,b)$，满足：

  • $v,b \in C_b(\mathbb{R}_+;\dot{B}^{d/2-1}_{2,1}) \cap L^1(\mathbb{R}_+;\dot{B}^{d/2+1}_{2,1})$，
  • $a \in C(\mathbb{R}_+;\dot{B}^{d/2-1}_{2,1} \cap \dot{B}^{d/2}_{2,1}) \cap L^2(\mathbb{R}_+;\dot{B}^{d/2}_{2,1})$，

且当 $a_0 = 0$（即 $\rho_0 \equiv 1$）时，有定量收敛率： $$ \sqrt{\frac{\nu}{\mu}} \Big( |\rho-1|{L^\infty(\mathbb{R}+;\dot{B}^{d/2}_{2,1})}

  • |\nabla^2 Q v|{L^1(\mathbb{R}+;\dot{B}^{d/2-1}_{2,1})}
  • |P v - V|{L^\infty(\mathbb{R}+;\dot{B}^{d/2-1}_{2,1})}
  • |b - B|{L^\infty(\mathbb{R}+;\dot{B}^{d/2-1}_{2,1})} \Big)
    $$ $$
  • \sqrt{\frac{\mu}{\nu}} \Big( |P v_t - V_t,\mu\nabla^2(P v - V)|{L^1(\mathbb{R}+;\dot{B}^{d/2-1}_{2,1})}
  • |b_t - B_t,\eta\nabla^2(b - B)|{L^1(\mathbb{R}+;\dot{B}^{d/2-1}_{2,1})} \Big) \le C. $$

Theorem B（$L^p$-临界推广）

将空间推广至 $ \dot{B}^{d/p-1}_{p,1} $，允许 $2 \le p < 4$（$d=2$）或 $2 \le p \le \min\{4,2d/(d-2)\}$（$d \ge 3$），通过高低频分解处理不同尺度的正则性不匹配问题，获得类似收敛率，但需分别控制低频 $a_\ell$ 与高频 $a_h$ 分量。

Theorem C（磁重联应用）

给定任意 $M,T>0$, $0<c<1$，存在光滑初值 $(\rho_0,u_0,b_0)$ 满足：

$$ \|u_0\|_{\dot{B}^{1}_{2,1}} \simeq \|b_0\|_{\dot{B}^{1}_{2,1}} \simeq M,\quad a_0 \in \dot{B}^{3/2}_{2,1}(\mathbb{R}^3),\quad c \le \rho_0 \le c^{-1}, $$

使得 C-MHD 存在唯一全局光滑解 $(\rho,u,b)$，且磁力线在 $t=0$ 与 $t=T$ 之间发生拓扑等价性破坏——即不存在 $\mathbb{R}^3$ 上的同胚映射将 $b(0,\cdot)$ 的积分曲线族映为 $b(T,\cdot)$ 的积分曲线族。

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

a) 核心分析学工具

本文证明策略继承并发展了 [10] 对可压缩Navier–Stokes方程的大黏度极限分析，核心是 参考流方法（reference flow method） 与 Leray投影分解 相结合的能量估计框架。

  • Leray投影分解：将速度场 $v = Pv + Qv$ 分解为散度自由部分 $Pv$（对应不可压缩模态）与梯度部分 $Qv = -\nabla(-\Delta)^{-1}\mathrm{div}\,v$（对应声波模态）。该分解天然适配临界Besov空间，因 $P,Q$ 均为零阶Fourier乘子，连续作用于 $\dot{B}^{s}_{p,1}$（$s \le d/2$）。

  • 高斯–格林型能量泛函设计：针对 $Pv - V$ 和 $b - B$ 构造差分能量，利用 $P$ 的正交性消去非物理压力项；对 $Qv$ 和 $a$ 构造“声波能量”，重点控制 $\|\nabla^2 Qv\|_{L^1_t \dot{B}^{d/2-1}_{2,1}}$ 与 $\|a\|_{L^\infty_t \dot{B}^{d/2}_{2,1}}$，后者依赖于体积黏度 $\lambda$ 提供的强阻尼。

  • 最大正则性理论（Maximal Regularity）：引用 Theorem 2.14（Stokes系统在Besov空间的 $L^1_t$-最大正则性），将 $Qv$ 满足的椭圆–抛物耦合方程转化为：

    $$ Qv_t - \nu \Delta Qv = \mathcal{F}(v,b,a), $$

    并直接获得 $\|Qv\|_{L^\infty_t \dot{B}^{d/2-1}_{2,1}} + \|\nabla^2 Qv\|_{L^1_t \dot{B}^{d/2-1}_{2,1}} \lesssim \|\mathcal{F}\|_{L^1_t \dot{B}^{d/2-1}_{2,1}}$，这是闭合先验估计的关键。

  • 实插值与Bernstein不等式：利用 Proposition 2.15 的插值不等式

    $$ \|f\|_{\dot{B}^{1/2}_{2,1}} \le 4 \|f\|_{L^2}^{1/2} \|f\|_{\dot{H}^1}^{1/2}, $$

    将 $L^2$ 能量与 $\dot{H}^1$ 耗散耦合；同时，Bernstein不等式 (2.7)(2.8) 用于精细控制各频率块的导数增益与衰减。

b) 非线性项处理技巧

  • Bony分解与乘积估计：对所有双线性项（如 $v \cdot \nabla v$, $b \cdot \nabla b$, $v \cdot \nabla b$）应用 Bony 分解 $gh = T_g h + T_h g + R(g,h)$，再调用 Lemma 2.7–2.9 的连续性估计。例如，对 $T_v \nabla v$，利用 $\|T_v \nabla v\|_{\dot{B}^{d/2-1}_{2,1}} \lesssim \|v\|_{\dot{B}^{d/2}_{2,1}} \|\nabla v\|_{\dot{B}^{d/2-1}_{2,1}}$，再由插值控制。

  • Commutator估计：对形如 $[v \cdot \nabla, \dot{\Delta}_j] f$ 的交换子，使用 Lemma 2.10：

    $$ \sum_j 2^{j(d/2-1)} \|[v \cdot \nabla, \dot{\Delta}_j] f\|_{L^2} \lesssim \|\nabla v\|_{L^\infty} \|f\|_{\dot{B}^{d/2-1}_{2,1}}, $$

    再结合 $\dot{B}^{d/2}_{2,1} \hookrightarrow L^\infty$（由 (2.10)），将 $L^\infty$ 范数控制为 Besov 范数。

  • 高低频分解技术（Theorem B）：为突破 $L^2$ 框架限制，引入频率截断 $f_\ell = \sum_{2^k \nu \le 1} \dot{\Delta}_k f$, $f_h = \sum_{2^k \nu > 1} \dot{\Delta}_k f$（见 (2.17)），分别用 $\dot{B}^{d/2-1}_{2,1}$ 控制低频 $a_\ell$ 的慢变性，用 $\dot{B}^{d/p}_{p,1}$ 控制高频 $a_h$ 的强耗散，避免 $p \ne 2$ 时的嵌入失效。

  • Composition估计：对压力项 $P(\rho) = A \rho^\gamma$，写作 $P = A(1+a)^\gamma$，利用 Lemma 2.13：

    $$ \|P(a) - A\|_{\dot{B}^{d/2}_{2,1}} \lesssim \|a\|_{\dot{B}^{d/2}_{2,1}},\quad \text{当 } a \text{ 有界}, $$

    从而将非线性压力转化为关于 $a$ 的线性主导项。

4. 与之前工作的比较

a) 突破与改进

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

第一，正则性假设的临界性突破。 此前关于可压缩MHD全局解的工作（如 [41,61]）均要求 $H^3$ 或更高 Sobolev 正则性，隐含 $L^2$-基底的强衰减；而 [20] 虽进入临界Besov框架，但仅处理小初值。本文首次在 任意大的临界Besov初值 下建立全局适定性，且不依赖小参数，仅需体积黏度 $\lambda$ 足够大——这是对“临界正则性 vs 大数据”矛盾的根本性调和。

第二，极限过程的定量性与强收敛性。 文献 [58] 使用紧性方法证明弱解的大体积黏度极限，但仅得弱收敛与存在性，无速率、无唯一性、无正则性传递。本文在强解框架下，不仅证明 $L^\infty_t \dot{B}^{d/2-1}_{2,1} \cap L^1_t \dot{B}^{d/2+1}_{2,1}$ 中的强收敛，更给出 显式代数速率 $\mathcal{O}(\sqrt{\mu/\nu})$ 与 $\mathcal{O}(\sqrt{\nu/\mu})$，揭示了体积黏度与剪切黏度的竞争平衡机制。

第三，物理现象构造的首创性。 磁重联的严格数学构造此前仅见于不可压缩MHD（[5,7,45]），依赖于 $2D$ 涡量输运或 $3D$ 精巧对称性。本文首次将重联解稳定地提升至可压缩框架，证明即使存在显著密度变化（$a \ne 0$），只要体积黏度主导，重联拓扑仍能鲁棒保持。这为理解真实天体物理中（如日冕物质抛射）的可压缩重联提供了首个严格模型。

综上，本文将可压缩MHD的数学理论推进至临界正则性、定量极限、物理可实现性三位一体的新高度。

2. 半空间中具有局部边界数据的Stokes方程组的流场反向现象

Flow reversal of the Stokes system with localized boundary data in the half space

作者: Tongkeun Chang, Kyungkeun Kang, Chanhong Min

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05858v1

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

详细分析

中文标题: 半空间中具有局部边界数据的Stokes方程组的流场反向现象

研究对象和背景

本文研究的是非稳态Stokes方程组（unsteady Stokes system）在上半空间 $\mathbb{R}^n_+ := \mathbb{R}^{n-1} \times \mathbb{R}_+$（其中 $n \geq 2$）中的解行为，特别关注当边界数据具有时空紧支集（space-time localized）且仅作用于平面边界 $\{x_n = 0\}$ 时所诱导的速度场符号变化现象，即所谓的 flow reversal（流场反向）。该现象表现为：在某个固定时刻或时间区间内，速度分量沿法向变量 $x_n$ 发生至少一次变号；更进一步地，在靠近边界与远离边界的区域，不同分量呈现相反的符号关系——这是边界层分离（boundary layer separation）的重要前兆。

论文考虑如下初边值问题：

$$ \begin{cases} \partial_t w - \Delta w + \nabla p = 0, & \text{div } w = 0, \quad \text{in } \mathbb{R}^n_+ \times (0,\infty), \\ w(x,t=0) = 0, & \\ w(x',x_n=0,t) = g(x',t), & \text{on } \partial\mathbb{R}^n_+ \times (0,\infty), \end{cases} $$

其中 $w : \mathbb{R}^n_+ \times (0,\infty) \to \mathbb{R}^n$ 是流体速度场，$p : \mathbb{R}^n_+ \times (0,\infty) \to \mathbb{R}$ 是压力，$g : \mathbb{R}^{n-1} \times (0,\infty) \to \mathbb{R}^n$ 是给定的非零、光滑（仅关于空间变量）、时空紧支的边界数据。具体地，作者选取单一分量非零的形式：

$$ g(y',s) = \psi(y') \phi(s) e_n, $$

其中 $\psi \in C^\infty_c(B'_1)$ 非负、$\phi \in C^\infty([0,1)) \cap L^1(\mathbb{R})$ 满足支撑在 $(0,1)$ 内，并在 $(\frac{1}{2},1)$ 上取幂律形式 $\phi(s) = (1-s)^a$，参数 $a > -1$ 控制奇异性强度；记 $M = \int_0^{1/2} \phi(s)\,ds$。

值得注意的是，该问题设定具有几个关键特殊性：

  • 零初始数据：$w(\cdot,0) \equiv 0$，排除了由初始扰动引发的动态效应，使现象纯粹源于边界驱动；
  • 非齐次但局部化的边界条件：$g$ 不恒为零，但其支撑集在空间与时间上均受限，从而模拟“瞬时、局域”的注入或抽吸；
  • 无滑移边界条件仅施加于平面边界：即不涉及曲面或复杂几何，便于利用对称性和Green函数结构；
  • 函数空间未预设正则性：解的正则性需从表示公式中逐点推导，尤其关注近边界区域（如 $x_n \ll 1$）的奇异行为；
  • 维度 $n \geq 2$，涵盖物理上重要的二维与三维情形。

这种设定区别于经典热方程或Laplace方程的局部正则性理论——由于Stokes系统的非局部性（nonlocal effect），即使在远离 $ \operatorname{supp} g $ 的区域，解也可能表现出法向导数爆炸（$\|\partial_{x_n} w\|_{L^\infty} = \infty$）或速度本身无界（$\|w\|_{L^\infty} = \infty$），正如文献[5]所示。因此，传统Caccioppoli型不等式在边界附近失效，亟需发展新的点态估计工具。

此外，该工作紧密联系流体力学中长期关注的核心问题：边界层分离机制。实验与数值研究表明，分离点常伴随速度切向分量变号（即所谓 separation point），而理论建模（如Prandtl方程、triple-deck模型）亦指出，逆压梯度（adverse pressure gradient，即 $\partial_x p_E > 0$）是触发分离的关键诱因。本文则试图在线性化模型（Stokes系统）中，通过构造显式反例，揭示即使在无非线性项的情形下，“流场反向”仍可由纯边界驱动自发产生，从而为理解Navier–Stokes中分离的线性起源提供严格依据。

主要定理或主要结果

论文建立了三个核心定理，分别刻画了分离点的存在性、零点渐近位置及内部反转点的存在性：

定理 1.6（分离点存在性）

设 $-1 < a \leq -\frac{1}{2}$，$w$ 是上述Stokes系统以 $g = g_n e_n$ 为边界数据的解。则存在充分大的常数 $N > 1$，使得对任意 $|x'| > N$，点 $(x',0,1)$ 是 $w_i$（$1 \le i \le n-1$）的一个分离点，即满足定义1.1中关于切向速度在边界附近变号的四条几何条件。

定理 1.8（零点渐近行为）

记 $x^*_{kn}(x',t)$ 为 $w_k(x',x_n,t)=0$ 关于 $x_n$ 的解（即第 $k$ 个分量的零点）。则存在 $N>1$ 和 $\delta = \delta(|x'|)>0$，使得当 $|x'|>N$ 且 $0<|t-1|<\delta$ 时，对 $i \ne n$ 有：

$$ x^*_{ni} \approx \begin{cases} (a+1)^{\frac{1}{2}} |x’|^{1-2a} |t-1|^{2a} \left( \left|a+\frac{1}{2}\right| |x’|^2 \gtrsim 1 \right) \

  • \left|a+\frac{1}{2}\right|^{\frac{1}{2}} |x’|^{1-2a} |t-1|^{2a} \left( \left|a+\frac{1}{2}\right| |x’|^2 \lesssim 1 \right), & -1<a<-\frac{1}{2} \ |x’|^{-1} \ln \frac{1}{|t-1|}, & a = -\frac{1}{2} \ |x’|^{\frac{1}{2}} (t-1)^{\frac{1}{2}}, & a > -\frac{1}{2},, t>1 \end{cases} $$

而法向分量 $x^*_{nn}$ 具有更复杂的指数-高斯混合渐近式，含 $e^{-|x'|^2/(t-1)}$ 项。

定理 1.10（反转点存在性）

对同一解 $w$，存在 $c>0$ 与 $0<\delta \ll 1$，使得当 $|x'| \ge c \sqrt{1+t}$ 时：

  • 若 $-1<a<0$，则对每个 $t>1$，$w_i(x',x_n,t)$（$i \ne n$）至少有一个反转点，$w_n(x',x_n,t)$ 至少有两个；
  • 若 $a>0$，则对每个 $t>1-\delta$，上述结论同样成立。

这些定理共同表明：仅靠局部、单向的边界注入，即可在Stokes流中激发出多尺度的符号振荡结构——切向分量至少一变号，法向分量至少两变号，且二者在近边界与远边界处符号相反，构成典型的“回流”（recirculation）模式。

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

a) 核心分析学工具

本文证明高度依赖于Green函数表示法与精细点态估计，而非能量方法或抽象泛函框架。其主干工具包括：

  • Stokes半空间Green张量显式公式：基于Golovkin与Solonnikov的经典工作，解可表为：

    $$ w_i(x,t) = \sum_{j=1}^n \int_0^t \int_{\mathbb{R}^{n-1}} K_{ij}(x-y',t-s) g_j(y',s)\,dy'ds, $$

    其中 $K_{ij}$ 由热核 $\Gamma$、Newtonian核 $N$ 及其导数组合而成（见公式(2.2)）。此表示允许将解分解为若干可独立控制的项：$w^{(L)}$（主奇异积分项）、$w^B$（边界修正项）、$w^N$（静态势项）、$w^G$（热核梯度项）。

  • 带权积分估计引理：如引理2.1给出关键积分

    $$ I_{m,k}(x,t) = \int_0^{x_n} z_n^k e^{-\frac{(x_n - z_n)^2}{4t}} \left( |x'|^2 z_n + z_n^2 \right)^{-m/2} dz_n $$

    的渐近等价式，其核心是分区域控制（$z_n \lesssim t$ 与 $z_n \gtrsim t$）并利用缩放与Laplace方法。

  • 卷积型核估计引理（引理2.2）：将 $\Gamma'$ 与 $\partial_{z_i}\partial_{z_j} N$ 的卷积分解为主项 $\partial_{x_i}\partial_{x_j} N(x)$ 加误差项 $J_\ell$，后者被证明具有更高阶衰减（如 $|J_1| \lesssim t x_n / |x|^{n+2}|x'|$）。这使得对 $L_{ij}$ 的控制转化为对 $N$ 导数的控制。

  • 奇异积分重写与主值处理（命题2.4）：通过引入截断球 $B_\epsilon(0)$ 并分析边界积分，将 $L_{ij}$ 精确分解为 $L_{ij} = \widetilde{L}_{ij} + \text{correction terms}$，其中 $\widetilde{L}_{ij}$ 是主奇异部分，其余为可显式计算的低阶项（如 $\partial_{x_n} \Gamma$ 或 $B_{in}$）。

  • 零点存在性判据（定义2.9与引理2.7）：借助中间值定理，将反转点存在性归结为证明某函数在区间两端异号；而零点位置估计则转化为求解形如 $\theta^a \ln \theta = -M$ 的超越方程，引理2.7为此类方程提供了精确的渐近根界。

b) 技术性难题处理技巧

  • 奇异性分离与主导项提取：针对不同 $a$ 值（$a < -1/2$, $a = -1/2$, $a > -1/2$），分别识别积分中起主导作用的区域（如 $s \approx 1$ 或 $s \approx 0$），并将被积函数按主导幂次展开，舍弃高阶小量。

  • 空间-时间耦合尺度分析：引入变量替换 $s = 1 - \sigma (t-1)$ 或 $z_n = \lambda \sqrt{t}$，使积分变量与参数 $x',t$ 形成统一尺度，从而显露真正的渐近主导项。

  • 符号控制的几何构造：为验证分离点定义中的不等式链，作者并非直接估计 $w_i$，而是构造辅助函数（如 $f_i(t) = \sup\{x_n > 0 : w_i(x',x_n,t) > 0\}$），并利用单调性与比较原理证明其收敛性与变号性。

  • 高斯衰减与代数衰减的竞争平衡：在 $x^*_{nn}$ 的估计中，$e^{-|x'|^2/(t-1)}$ 项与幂律项相互竞争，作者通过分区域讨论（$x_n \lesssim |x'|$ vs $x_n \gtrsim |x'|$）实现精确匹配。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重实质性推进：

  • 首次在纯线性系统中严格证实“流场反向”的边界驱动起源：以往关于分离点的理论（如Prandtl方程[29,9]、triple-deck[12,14]）均建立在非线性近似或无穷Reynolds数极限之上；而本文在完全线性、有限Reynolds数（即Stokes）框架下，仅凭局部边界数据即构造出符号变号解，表明分离的“种子”可独立于非线性对流项而存在。

  • 突破正则性假设限制：先前关于Stokes边界正则性的研究（如[22,3,6]）多要求边界数据足够光滑（$C^\infty$ in $t$）或力项满足特定衰减。本文仅假设 $g$ 关于时间连续但不可微（$\phi \in C^\infty([0,1))$ 但 $\phi^{(k)}$ 在 $s=1$ 处可能发散），且通过幂律奇异性参数 $a$ 系统刻画了正则性临界阈值：当 $a \le -1/2$ 时出现分离点，$a > 0$ 时需调整时间窗口，揭示了奇异性强度与动力学响应之间的定量对应。

  • 提供首个半空间Stokes解的全局点态渐近图谱：不同于仅关注 $L^p$ 或 $H^s$ 范数的文献，本文对每个速度分量 $w_i$ 给出了在 $(x',x_n,t)$ 全参数空间中的逐点等价式（$\approx$），覆盖近边界（$x_n \to 0$）、远场（$|x'| \to \infty$）、临界时间（$t \to 1$）等所有关键区域，为后续数值验证与物理建模提供了精确基准。

综上，本文不仅填补了线性化边界层理论的关键空白，更以严谨分析揭示了流体不稳定性的一种普适机制：局部边界扰动可通过非局部流体响应，在远离扰动源的空间中诱发全局拓扑重构——这一思想对微流控、生物流动及湍流初生研究均具深远启示。

3. 具有非Lipschitz泊松跳跃噪声的二维Boussinesq方程的弱平均与强平均原理

Weak and strong averaging principle for 2D Boussinesq equations with non-Lipschitz Poisson jump noise

作者: Yangyang Shi, Dong Su, Hui Liu

类别: math.AP, math.PR

发布日期: 2026-02-05

arXiv ID: 2602.05696v1

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

详细分析

中文标题: 具有非Lipschitz泊松跳跃噪声的二维Boussinesq方程的弱平均与强平均原理

1. 研究对象和背景

本文研究的是定义在二维环面 $\mathbb{T}^2 = [0,2\pi] \times [0,2\pi]$ 上的 二维Boussinesq方程组，该方程组描述了不可压缩粘性流体中速度场与温度场之间的耦合演化，是流体力学中刻画热对流现象的核心模型之一。其确定性形式在引言第1节中给出：

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u = \Delta u - \nabla p + \theta e_2, \\ \partial_t \theta + (u \cdot \nabla) \theta = \Delta \theta, \\ \nabla \cdot u = 0, \\ u(0,x,y) = u_0(x,y),\quad \theta(0,x,y) = \theta_0(x,y), \end{cases} $$

其中 $u(t,x,y) = (u_1,u_2)$ 是速度向量场, $p = p(t,x,y)$ 是标量压力, $\theta(t,x,y)$ 是温度（或密度扰动）标量场, $e_2 = (0,1)$ 表示重力方向（竖直方向），$\Delta$ 是拉普拉斯算子，$\nabla$ 是梯度算子。所有函数均满足周期边界条件，即定义在 $\mathbb{T}^2$ 上，且空间变量具有 $2\pi$-周期性。

为便于分析，作者对第一式取旋度（curl），引入涡度变量 $j := \nabla^\perp \cdot u = \partial_x u_2 - \partial_y u_1$，从而将原始速度-压力系统转化为关于 $(j,\theta)$ 的等价系统。这一转化消除了压力项，并简化了非线性结构；同时，由于 $\nabla \cdot u = 0$，速度场可通过 Biot–Savart 律由涡度重构：$u = \nabla^\perp (-\Delta)^{-1} j$。

本文的关键创新在于：在上述系统中引入了多尺度随机扰动，具体为两类非Lipschitz泊松跳跃噪声，驱动一个快慢耦合的随机偏微分方程（SPDE）系统。该系统以小参数 $\varepsilon > 0$ 刻画时间尺度分离，形式如下（见原文公式(1.2)）：

$$ \begin{cases} \partial_t j^\varepsilon + (u^\varepsilon \cdot \nabla) j^\varepsilon = \Delta j^\varepsilon + f(j^\varepsilon,\theta^\varepsilon) + \partial_x \theta^\varepsilon + \displaystyle\int_{|z|<1} \sigma_1(j^\varepsilon,\theta^\varepsilon,z)\,\tilde{\eta}_1(dt,dz), \\ \partial_t \theta^\varepsilon + \dfrac{1}{\varepsilon}(u^\varepsilon \cdot \nabla)\theta^\varepsilon = \dfrac{1}{\varepsilon}\Delta \theta^\varepsilon + \displaystyle\int_{|z|<1} \sigma_2(j^\varepsilon,\theta^\varepsilon,z)\,\tilde{\eta}_2^\varepsilon(dt,dz), \\ j^\varepsilon(0) = j_0,\quad \theta^\varepsilon(0) = \theta_0. \end{cases} $$

此处，$\tilde{\eta}_1$ 和 $\tilde{\eta}_2^\varepsilon$ 是定义在可测空间 $(H,\mathcal{B}(H))$ 上的补偿泊松随机测度（compensated Poisson random measures），其强度测度分别为 $\nu_1,\nu_2$，满足 $\int_{|z|<1}|z|^2\,\nu_i(dz) < \infty$（即二阶矩有限）。特别地，$\tilde{\eta}_2^\varepsilon(dt,dz) = \tilde{\eta}_2\left(\frac{dt}{\varepsilon},dz\right)$，体现温度变量 $\theta^\varepsilon$ 的快速演化特性，而涡度 $j^\varepsilon$ 相对缓慢——这是典型的双尺度系统（two-scale system）结构。

函数 $f$, $\sigma_1$, $\sigma_2$ 均不满足经典 Lipschitz 条件，而是满足非Lipschitz型增长与连续性条件（见 Assumptions 1–3），例如对 $f$ 要求：

$$ \|f(u_1,v_1) - f(u_2,v_2)\|_H^{2p} \le \kappa\big( \|u_1-u_2\|_H^{2p} + \|v_1-v_2\|_H^{2p} \big), $$

其中 $\kappa:\mathbb{R}_+\to\mathbb{R}_+$ 是凹、连续、非减、满足 Osgood 条件 $\int_{0^+} \frac{du}{\kappa(u)} = \infty$ 的模函数。此类条件显著弱于 Lipschitz 连续性，允许更强的奇异性（如平方根型、对数型非线性），常见于建模真实物理噪声（如脉冲型热源、间歇性湍流激发）。

函数空间设定方面，作者采用标准 Sobolev-Hilbert 框架：设 $H = L^2_0(\mathbb{T}^2) = \{ f \in L^2(\mathbb{T}^2):\int_{\mathbb{T}^2} f\,dxdy = 0 \}$，配备内积 $(\cdot,\cdot)$；对 $\alpha > 0$，记 $H^\alpha = D((-\Delta)^{\alpha/2})$，范数为 $\|u\|_{H^\alpha} = \|u\|_H + \|(-\Delta)^{\alpha/2}u\|_H$。特别地，$H^1$ 对应能量空间，$H^{-1}$ 为其对偶。所有解均取值于 $H \times H$，初始数据 $(j_0,\theta_0) \in H \times H$（部分估计中需 $j_0 \in H^1$）。

综上，本文研究对象是一个具有非Lipschitz跳跃噪声、快慢时间尺度分离、周期边界、涡度-温度耦合结构的随机 Boussinesq 系统。其物理背景涵盖大气热对流、地幔对流、工程冷却系统等存在突发性热扰动与多时间尺度相互作用的场景。

2. 主要定理或主要结果

本文核心成果围绕平均原理（averaging principle）展开，分为弱平均（convergence in probability）与强平均（$2p$-th mean convergence）两类。主要定理如下：

Theorem 3.2（适定性与先验估计）
假设 Assumptions 1–4 成立，且 $(j_0,\theta_0) \in H \times H$。则对任意 $T > 0$ 及某 $p \ge 1$，方程组 (1.2) 存在唯一 $H \times H$-值右连左极（càdlàg）适应解 $(j^\varepsilon,\theta^\varepsilon)$，并满足如下一致先验估计：

$$ \mathbb{E}\left[ \sup_{t \in [0,T]} \|j^\varepsilon\|_H^{2p} \right] + \mathbb{E}\left[ \int_0^T \|j^\varepsilon\|_H^{2(p-1)} \|\nabla j^\varepsilon\|_H^2 \,ds \right] \le C_T \left(1 + \|j_0\|_H^{2p} + \|\theta_0\|_H^{2p} \right), $$
$$ \sup_{t \in [0,T]} \mathbb{E}\left[ \|\theta^\varepsilon\|_H^{2p} \right] + \mathbb{E}\left[ \int_0^T \|\theta^\varepsilon\|_H^{2(p-1)} \|\nabla \theta^\varepsilon\|_H^2 \,ds \right] \le C_T \left(1 + \|j_0\|_H^{2p} + \|\theta_0\|_H^{2p} \right). $$

此外，对停时 $\tau_M := \inf\{ t \ge 0 : \|j^\varepsilon\|_{H^1} > M \}$，有正则性估计：

$$ \mathbb{E}\left[ \sup_{t \in [0,\tau_M \wedge T]} \|j^\varepsilon\|_{H^1}^2 \right] + \mathbb{E}\left[ \int_0^{\tau_M \wedge T} \|\Delta j^\varepsilon\|_H^2 \,ds \right] \le C_{T,M} \left(1 + \|j_0\|_{H^1}^2 + \|\theta_0\|_H^2 \right). $$

Theorem 4.1（温度变量的遍历性）
在固定“冻结”涡度 $j$ 下，快变温度方程的“冻结方程”（frozen equation）具有唯一不变概率测度 $\mu_j$，且对任意初值 $\theta_0 \in H$，其解 $\theta_t^{j}$ 满足：

$$ \left| \mathbb{E}\big[\varphi(\theta_t^{j})\big] - \int_H \varphi(y)\,\mu_j(dy) \right| \le C_j \, e^{-\gamma_j t} \|\varphi\|_{\mathrm{Lip}}, \quad \forall \varphi \in \mathrm{Lip}(H), $$

即指数混合（exponential mixing），其中 $\gamma_j > 0$ 依赖于 $j$。

Theorem 5.1（平均原理）
设 $\bar{j}$ 是如下平均方程的解：

$$ \partial_t \bar{j} + (\bar{u} \cdot \nabla) \bar{j} = \Delta \bar{j} + \bar{f}(\bar{j}) + \partial_x \bar{\theta} + \int_{|z|<1} \bar{\sigma}_1(\bar{j},z)\,\tilde{\eta}_1(dt,dz), $$

其中 $\bar{f}(j) = \int_H f(j,\theta)\,\mu_j(d\theta),\; \bar{\sigma}_1(j,z) = \int_H \sigma_1(j,\theta,z)\,\mu_j(d\theta)$，且 $\bar{\theta} = \int_H \theta\,\mu_j(d\theta)$ 为温度不变分布的均值。则：

  • 弱平均：当 $\varepsilon \to 0$ 时，$j^\varepsilon \to \bar{j}$ 在概率意义下收敛，即对任意 $\delta > 0$，

    $$ \lim_{\varepsilon \to 0} \mathbb{P}\left( \sup_{t \in [0,T]} \|j^\varepsilon(t) - \bar{j}(t)\|_H > \delta \right) = 0. $$

  • 强平均：若进一步假设 $\sigma_2$ 关于 $\theta$ 是 Lipschitz 的（即 $L_{\sigma_2} < \infty$），则存在常数 $C > 0$ 使得：

    $$ \mathbb{E}\left[ \sup_{t \in [0,T]} \|j^\varepsilon(t) - \bar{j}(t)\|_H^{2p} \right] \le C \varepsilon^{1/2}. $$

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

a) 核心分析学工具

本文证明框架建立在随机分析、泛函分析与偏微分方程先验估计的深度交叉之上，关键工具包括：

  • Itô 公式（含 Lévy 噪声）：全文能量估计（如 Theorem 3.2 中的 $2p$-阶矩估计）均基于广义 Itô 公式。对 $F(x) = \|x\|_H^{2p}$ 应用 Itô–Lévy 公式，产生漂移项、扩散项及跳项三部分。跳项处理依赖于补偿测度性质与 Lévy 测度的二阶矩有限性，如：

    $$ \mathbb{E}\left[ \int_0^t \int_{|z|<1} \Big( \|x+z\|_H^{2p} - \|x\|_H^{2p} - 2p\|x\|_H^{2(p-1)}(x,z)_H \Big)\,\nu(dz)ds \right]. $$

    此处利用了 Assumption 2 中的非Lipschitz增长控制及引理 (2.1) 的 Taylor 型余项估计。

  • Burkholder–Davis–Gundy (BDG) 不等式：用于控制随机积分项的高阶矩。例如，在估计 $j^\varepsilon$ 的 $2p$-阶矩时，对鞅项

    $$ M_t = \int_0^t \int_{|z|<1} \|j^\varepsilon\|_H^{2(p-1)} (j^\varepsilon,\sigma_1(j^\varepsilon,\theta^\varepsilon,z))\,\tilde{\eta}_1(ds,dz) $$

    应用 BDG 得到：

    $$ \mathbb{E}\left[ \sup_{s \le t} |M_s| \right] \le C_p \mathbb{E}\left[ \left( \int_0^t \int_{|z|<1} \|j^\varepsilon\|_H^{4p-2} \|\sigma_1\|_H^2 \,\nu_1(dz)ds \right)^{1/2} \right], $$

    再结合 Young 不等式与吸收技巧完成闭环。

  • Gronwall 不等式（随机版本）：几乎所有先验估计最终归结为形如

    $$ \frac{d}{dt} \mathbb{E}[Y_t] \le -\lambda_\varepsilon \mathbb{E}[Y_t] + C \mathbb{E}[Z_t], \quad \lambda_\varepsilon \sim \varepsilon^{-1}, $$

    的微分不等式。通过 Gronwall 引理得到指数衰减控制，如 (3.4) 式所示：

    $$ \mathbb{E}[\|\theta^\varepsilon(t)\|_H^{2p}] \le C e^{-\lambda_{p,\gamma} t/\varepsilon} + \frac{C}{\varepsilon} \int_0^t e^{-\lambda_{p,\gamma}(t-s)/\varepsilon} \mathbb{E}[\|j^\varepsilon(s)\|_H^{2p}]\,ds. $$

  • 紧性论证（Aldous–Kurtz 准则）：为证平均原理，需证 $\{j^\varepsilon\}_{\varepsilon>0}$ 在 $D([0,T];H)$ 中紧。Lemma 3.1 给出模连续性估计：

    $$ \mathbb{E}\left[ \|j^\varepsilon(t) - j^\varepsilon(t-\delta)\|_H^2 \right] \le C_{T,M} \delta^{1/2}, $$

    结合有界性 $\sup_\varepsilon \mathbb{E}[\sup_{t} \|j^\varepsilon(t)\|_H^{2p}] < \infty$，由 Aldous 准则得紧性，进而存在收敛子列。

  • Krylov–Bogoliubov 方法与 Doeblin 条件：用于证明快变温度方程的遍历性（Theorem 4.1）。通过构造 Lyapunov 函数 $V(\theta) = \|\theta\|_H^2$ 并验证耗散性（Assumption 4 保证 $\lambda_p > 0$），结合 Doeblin 小集条件，推出唯一不变测度 $\mu_j$ 的存在性与指数收敛性。

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

  • 非Lipschitz 非线性项的处理：传统平均原理要求系数满足 Lipschitz 条件以控制误差项 $\mathbb{E}\|f(j^\varepsilon,\theta^\varepsilon) - \bar{f}(j^\varepsilon)\|$。本文采用 Osgood 型模连续性 + 遍历性平均的稳定性：利用 $\mu_j$ 对 $j$ 的连续依赖性（由 Assumption 1 的凹模 $\kappa$ 控制），结合 $\theta^\varepsilon$ 的快速混合，证明

    $$ \left\| \mathbb{E}[f(j,\theta^\varepsilon_t)] - \bar{f}(j) \right\|_H \le C \kappa\big( \mathbb{E}\|j^\varepsilon - j\|_H^{2p} \big)^{1/2p} + C e^{-\gamma t/\varepsilon}. $$

    此处 $\kappa$ 的 Osgood 性质确保误差可被吸收。

  • 对流项 $(u^\varepsilon \cdot \nabla) j^\varepsilon$ 的控制：该非线性项破坏能量守恒，需借助 Ladyzhenskaya 不等式 与插值：

    $$ |((u^\varepsilon \cdot \nabla) j^\varepsilon,\Delta j^\varepsilon)| \le C \|j^\varepsilon\|_{H^1} \|j^\varepsilon\|_{H^2}^2, $$

    再通过停时 $\tau_M$ 截断，将高阶范数控制在 $M$ 内，实现局部一致估计。

  • 快慢变量耦合的解耦技巧：为分离尺度，作者引入修正函数（corrector）$\chi(j,\theta) = \theta - \bar{\theta}(j)$，并证明 $\chi$ 关于 $\theta$ 指数衰减。利用此性质，将 $j^\varepsilon$ 的误差分解为慢变部分与快变振荡部分，后者由遍历性平均消除。

4. 与之前工作的比较

a) 突破与改进

本文相较已有文献实现了三重实质性突破：

  • 噪声类型的拓展：此前关于 Boussinesq 方程平均原理的工作（如 [36–39]）均假设高斯白噪声或Lipschitz 连续的 Lévy 噪声。本文首次处理非Lipschitz泊松跳跃噪声，允许 $\sigma_i$ 含有 $\sqrt{\|\theta\|_H}$、$\log(1+\|\theta\|_H)$ 等奇异依赖，更贴合实际物理中的脉冲型热扰动（如火山喷发、核反应瞬态加热）。技术上，这迫使作者放弃传统的 Picard 迭代与 Lipschitz 压缩映射，转而采用 Osgood 型单调性+截断+紧性 的组合策略。

  • 非线性结构的放宽：既有工作（如 Du et al. [40]）要求快变方程为抽象反应扩散型 $\partial_t \theta = \frac{1}{\varepsilon} A \theta + \frac{1}{\varepsilon} g(j,\theta) + \text{noise}$，其中 $A$ 为强椭圆算子。本文中快变温度方程含对流项 $\frac{1}{\varepsilon}(u^\varepsilon \cdot \nabla)\theta^\varepsilon$，破坏了标准半群生成性，使冻结方程的遍历性证明异常困难。作者通过发展 依赖于涡度 $j$ 的加权 Lyapunov 函数 与 $j$-参数化 Doeblin 小集，克服了该障碍。

  • 收敛速率的量化：多数平均原理仅给出弱收敛（如 [23,33]），而本文在附加 Lipschitz 条件下获得显式强收敛速率 $\mathbb{E}\sup_t \|j^\varepsilon - \bar{j}\|_H^{2p} = O(\varepsilon^{1/2})$，优于经典 Khasminskii 方法的 $O(\varepsilon^{1/3})$。这得益于对跳项误差的精细估计：利用 Lévy 测度的二阶矩与 $\sigma_1$ 的非Lipschitz 模，将跳误差控制为 $\varepsilon^{1/2} \|\sigma_1\|_{\mathrm{Lip}}$ 量级。

综上，本文将随机 Boussinesq 方程的平均原理推进至更真实、更奇异、更复杂的噪声与非线性框架，为多尺度流体模型的降维与不确定性量化提供了新范式。

4. 三维外区域中具退化黏性及远场真空的可压磁流体力学方程初边值问题的整体强解

Strong solutions to the initial-boundary-value problem of compressible MHD equations with degenerate viscosities and far field vacuum in 3D exterior domains

作者: Jiaxu Li, Boqiang Lü, Bing Yuan

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05264v1

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

详细分析

中文标题: 三维外区域中具退化黏性及远场真空的可压磁流体力学方程初边值问题的整体强解

1. 研究对象和背景

本文研究的是三维外区域 $\Omega = \mathbb{R}^3 \setminus \overline{D}$ 中，具退化黏性系数与远场真空的可压磁流体力学（MHD）方程组的初边值问题（IBVP）。该系统描述了在电磁场作用下、具有可变密度的导电流体（如等离子体）的动力学行为。其物理机制比单纯流体更复杂，因速度场 $u$ 与磁场 $H$ 相互耦合、非线性反馈强烈，且在密度趋于零（即真空）时，黏性系数退化，导致控制方程出现严重奇异性。

系统的基本控制方程为：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho u) = 0, \\ (\rho u)_t + \mathrm{div}(\rho u \otimes u) + \nabla P = \mathrm{div}\, T + H \cdot \nabla H - \nabla \left( \frac{|H|^2}{2} \right), \\ H_t - \eta \Delta H = \mathrm{curl}(u \times H), \\ \mathrm{div}\, H = 0, \end{cases} \tag{1.1} $$

其中：

  • $\rho = \rho(x,t) \ge 0$ 是密度；
  • $u = (u_1,u_2,u_3)(x,t)$ 是速度场；
  • $P = a \rho^\gamma$（$a>0,\gamma >1$）是绝热压力；
  • $H = (H_1,H_2,H_3)(x,t)$ 是磁场；
  • $\eta > 0$ 是磁扩散系数；
  • $T$ 是黏性应力张量，定义为
    $$ T = 2\mu(\rho) D(u) + \lambda(\rho)\, \mathrm{div}\, u\, I, \tag{1.2} $$
    其中 $D(u) = \frac{1}{2}\big(\nabla u + (\nabla u)^T\big)$ 是形变率张量，$I$ 是 $3\times3$ 单位阵；
  • 黏性系数取为密度的幂函数：
    $$ \mu(\rho) = \mu \rho^\delta,\quad \lambda(\rho) = \lambda \rho^\delta,\quad 0 < \delta < 1, \tag{1.3} $$
    且满足物理相容性条件 $\mu > 0$, $2\mu + 3\lambda \ge 0$。这种选择导致当 $\rho \to 0$（即远场真空）时，$\mu(\rho),\lambda(\rho) \to 0$，故称退化黏性。

研究区域 $\Omega$ 是 $\mathbb{R}^3$ 中去掉一个光滑有界单连通区域 $D$ 后的外区域，即 $\Omega = \mathbb{R}^3 \setminus \overline{D}$。这是典型的无界但非全空间的几何设置，带来边界效应与无穷远衰减双重挑战。

边界条件采用Navier-滑移型速度边界条件与理想导体型磁场边界条件：

  • 速度边界条件（Navier-slip）：

    $$ u \cdot n = 0,\quad \mathrm{curl}\, u \times n = -A(x) u,\quad \text{on } \partial\Omega, \tag{1.6} $$

    其中 $n$ 是指向 $\Omega$ 外侧的单位法向量，$A$ 是定义在 $\partial\Omega$ 上的对称矩阵，刻画边界摩擦效应。

  • 磁场边界条件（perfect conducting）：

    $$ H \cdot n = 0,\quad \mathrm{curl}\, H \times n = 0,\quad \text{on } \partial\Omega. \tag{1.7} $$

  • 远场行为（far-field vacuum）：

    $$ (\rho, u, H) \to (0, 0, 0),\quad \text{as } |x| \to \infty. \tag{1.8} $$

    这意味着初始密度在无穷远处趋于零，即存在真空状态，而非传统“非真空小扰动”框架。

函数空间设定亦具特殊性：解需满足加权 Sobolev 正则性，例如要求 $\rho^{\frac{1-\delta}{2}} u \in L^2$, $\nabla \rho^{\delta - \frac{1}{2}} \in D^{1} = \{ f \in L^1_{\mathrm{loc}} : \|\nabla f\|_{L^2} < \infty \}$, 以及关键新变量 $J \triangleq \rho^{-1+\frac{2\delta}{3}} H \in H^1$。这些空间选择直接对应于退化黏性带来的奇异性控制需求——不能用标准 $H^1$ 或 $W^{1,p}$ 控制 $\nabla u$，而必须引入密度权重以平衡 $\rho^\delta$ 的消失。

特别地，文中定义的密度加权平均磁场

$$ J \triangleq \rho^{-1+\frac{2\delta}{3}} H $$


是全文技术核心之一。它将磁场与密度退化耦合显式分离，使 $J$ 满足一个标准抛物型方程（见(1.20)），从而规避了原始方程中 $\rho^{-1+\frac{2\delta}{3}} \nabla H$ 类奇异性项。

综上，本工作聚焦于一个兼具三重困难的数学模型：
（i）几何困难：外区域 $\Omega$ —— 非紧致、含边界、需处理无穷远衰减；
（ii）奇异性困难：$\delta \in (0,1)$ 导致 $\mu(\rho) = \mu \rho^\delta$ 在 $\rho=0$ 处不可微，且 $\nabla \rho^{\delta-1}$ 在真空附近爆炸；
（iii）耦合困难：MHD 系统中 $u$ 与 $H$ 通过 Lorentz 力 $H\cdot\nabla H$ 与感应方程 $\partial_t H = \mathrm{curl}(u\times H) + \eta \Delta H$ 强耦合，远超 Navier–Stokes 框架。

2. 主要定理或主要结果

本文的核心结论是以下局部强解存在唯一性定理：

Theorem 1.1.
设参数 $(\gamma,\delta)$ 满足

$$ \gamma > 1,\quad 0 < \delta < 1. \tag{1.11} $$


若初始数据 $(\rho_0, u_0, H_0)$ 满足如下正则性与相容性条件：

  • 正则性条件：

    $$ \begin{cases} \rho_0^{\frac{1-\delta}{2}} u_0 \in L^2,\quad u_0 \in D^1 \cap D^2,\\ H_0 \in H^2,\quad \mathrm{div}\, H_0 = 0,\\ J_0 \triangleq H_0 / \rho_0^{1+\frac{2\delta}{3}} \in H^1,\\ \rho_0^{\gamma-1+\frac{2\delta}{3}} \in D^1 \cap D^2,\quad \nabla \rho_0^{\delta-\frac{1}{2}} \in D^1,\\ \rho_0^{\delta-\frac{1}{2}} \in L^6(\Omega_0), \end{cases} \tag{1.12} $$

    其中 $\Omega_0 = \Omega \cap B_{2R_0}$ 是截断区域；

  • 相容性条件：存在 $g_0 \in L^2$，使得

    $$ L(u_0) = \rho_0^{\frac{1-\delta}{2}} g_0, \tag{1.13} $$

    其中 $L(u) = \mu \Delta u + (\mu+\lambda)\nabla \mathrm{div}\, u$ 是 Lamé 算子。

则存在 $T_0 > 0$，使得初边值问题 (1.1)–(1.8) 在 $[0,T_0]\times \Omega$ 上存在唯一强解 $(\rho,u,H)$，满足：

$$ \begin{aligned} &\rho^{\gamma-1+\frac{2\delta}{3}} \in L^\infty([0,T_0]; W^{1,6}\cap D^1\cap D^2),\\ &\nabla \rho^{\delta-\frac{1}{2}} \in L^\infty([0,T_0]; D^1),\quad \rho^{\delta-\frac{1}{2}} \in L^\infty([0,T_0]; L^6(\Omega_0)),\\ &u \in L^\infty([0,T_0]; D^1\cap D^2)\cap L^2([0,T_0]; D^3),\\ &H \in L^\infty([0,T_0]; H^2)\cap L^2([0,T_0]; D^3),\\ &J \in L^\infty([0,T_0]; H^1)\cap L^2([0,T_0]; D^2),\\ &u_t \in L^2([0,T_0]; D^1),\quad H_t \in L^2([0,T_0]; H^1),\quad J_t \in L^2([0,T_0]; L^2),\\ &\rho^{\frac{1-\delta}{2}} u_t,\, H_t \in L^\infty([0,T_0]; L^2),\\ &\rho^{\delta-\frac{1}{2}} L u,\; \rho^{\delta-\frac{1}{2}} \nabla \mathrm{div}\, u \in L^2([0,T_0]; H^1). \end{aligned} \tag{1.14} $$

该定理不仅给出解的存在唯一性，更揭示了一个关键物理现象：由 $J_0 \in H^1$ 可推出 $H_0$ 比 $\rho_0^{\frac{1+\delta}{2}}$ 衰减更快（Remark 1.3）；而定理保证 $J \in L^\infty([0,T_0]; H^1)$，说明该磁场衰减速率优势在整个演化过程中保持不变。这表明磁场在抑制退化黏性奇异性方面具有主动稳定作用，是区别于纯 Navier–Stokes 系统的本质特征。

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

a) 核心分析学工具

本文采用经典能量方法结合加权估计与先验估计闭包策略，辅以现代椭圆/抛物型理论工具。具体包括：

  • 加权 Sobolev 空间理论：针对退化黏性，引入密度权重空间 $D^{k,p}_\omega = \{f : \|\omega f\|_{W^{k,p}} < \infty\}$，如 $\rho^{\delta-\frac{1}{2}}$ 加权控制 $\nabla u$；

  • Lamé 方程正则性引理（Lemma 2.3）：对速度满足的椭圆系统 $L u = f$，在 Navier-slip 边界下建立高阶估计

    $$ \|\nabla^{k+2} u\|_{L^p} \le C \big( \|f\|_{W^{k,p}} + \|\nabla u\|_{L^2} \big), $$

    是控制 $u$ 高阶导数的基础；

  • 带曲率边界的 Poisson 方程估计（Lemma 2.4）：对 $J$ 满足的 $\Delta J = F$，在边界含 $\psi = \nabla \log \rho$ 的非齐次条件 $ \mathrm{curl}\, J \times n = -\frac{1+\delta}{2}(\psi \cdot n) J$ 下，建立

    $$ \|\nabla^2 J\|_{L^2} \le C \big( \|F\|_{L^2} + (1+\|\psi\|_{D^1}^2) \|J\|_{H^1} \big), $$

    该估计将 $\psi$ 的低阶范数控制高阶 $J$，是克服边界奇异性关键；

  • Gagliardo–Nirenberg 不等式（Lemma 2.1）与 div-curl 分解（Lemma 2.2）：用于将 $L^p$ 范数控制为散度与旋度组合，特别适用于 $u\cdot n = 0$ 的滑移边界；

  • 截断区域逼近法（annular domain $\Omega_R$）：先在有界区域 $\Omega_R = \Omega \cap B_R$ 上构造光滑近似解，再令 $R \to \infty$ 取极限，结合一致先验估计实现紧性收敛；

  • 时间加权能量泛函：定义

    $$ E(t) = 1 + \|\rho^{\frac{1-\delta}{2}} u\|_{L^2} + \|\nabla u\|_{L^2} + \|H\|_{H^1} + \|J\|_{H^1} + \|\rho^{\frac{1-\delta}{2}} u_t\|_{L^2} + \cdots, $$

    并证明 $\sup_{[0,T]} E(t) \le M$ 对某 $T>0$ 成立，从而闭合估计。

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

  • 退化黏性奇异性控制：原始动量方程含 $\mathrm{div}(\rho^\delta \nabla u)$，当 $\rho \to 0$ 时系数消失。作者乘以 $\rho^{-\delta}$ 得到

    $$ \rho^{1-\delta}(u_t + u\cdot \nabla u) - L u = -\rho^{1-\frac{2\delta}{3}} \nabla \rho^{\gamma-1+\frac{2\delta}{3}} - \cdots, \tag{1.17} $$

    此时 $L u$ 项保留，而右端奇异性被密度幂次吸收。

  • 磁场相关奇异性消除：项 $\rho^{-1+\frac{2\delta}{3}} H \cdot \nabla H$ 在 $\rho \to 0$ 时发散。引入 $J = \rho^{-1+\frac{2\delta}{3}} H$ 后，利用恒等式

    $$ \rho^{-1+\frac{2\delta}{3}} H \cdot \nabla H = J \cdot \nabla H,\quad \rho^{-1+\frac{2\delta}{3}} \nabla H \cdot H = \nabla H \cdot J, $$

    并推导出 $J$ 满足的抛物方程 (1.20)，其源项仅含 $J,\psi,\nabla J$，无 $\rho^{-\alpha}$ 奇异。

  • 边界项处理：$J$ 的边界条件含 $\psi\cdot n$，导致 $\int_{\partial\Omega} (\mathrm{curl}\, J \times n)\cdot J_t \, dS$ 出现 $\psi J_t$ 乘积。作者将其还原为 $H_t$ 表达式：

    $$ J_t = \rho^{-1+\frac{2\delta}{3}} H_t - H \, u \cdot \nabla \rho^{-1+\frac{2\delta}{3}} + \cdots, $$

    利用 $\rho$ 在边界有正下界（Lemma 2.5）及 $\nabla H_t \in L^2$，再应用迹定理控制。

  • 唯一性证明中的负幂困难：比较两个解时需估计 $\bar{g} = \rho_1^{-\delta} - \rho_2^{-\delta}$，其本身不具正则性。作者利用初始时刻 $\bar{g}(0)=0$ 及方程结构，将 $\bar{g}$ 的演化嵌入 $H$ 差量估计中，绕过直接控制 $\bar{g}$。

4. 与之前工作的比较

a) 突破与改进

本文在三个维度上显著推进了现有理论：

  • 首次建立 $\delta \in (0,1)$ 的 MHD 外区域强解理论：此前，Liu–Luo–Zhong [28] 仅处理 $\delta = 1$ 情形；Li–Lü–Yuan [24] 解决了 Navier–Stokes 对应问题，但未涉及磁场耦合。本文是首个将退化黏性、外区域、真空、磁场四者统一处理的工作。

  • 大幅放宽对初始磁场衰减的要求：以往结果（如 [19,32]）要求 $H_0 \in H^2$ 且无加权，隐含 $H_0$ 衰减慢于 $\rho_0$；而本文通过 $J_0 \in H^1$ 要求，等价于 $H_0 = O(\rho_0^{\frac{1+\delta}{2}})$，即磁场比密度更高阶衰减，这更符合物理中磁场在稀薄等离子体中快速弥散的图像。

  • 揭示磁场的“奇异性压制”机制：Theorem 1.1 的结论 $J \in L^\infty(H^1)$ 表明，初始的 $H/\rho^{\frac{1+\delta}{2}}$ 正则性在演化中不变，而 Navier–Stokes 中并无类似量。这说明磁场通过洛伦兹力提供额外阻尼，抑制 $\nabla \rho^{\delta-1}$ 类奇异性，是纯流体理论所不具备的新自由度。

综上，本文不仅填补了 $\delta \in (0,1)$ MHD 外区域理论的空白，更通过构造 $J$ 变量，首次量化了磁场在控制退化黏性奇异性中的主导作用，为后续研究真空 MHD 的全局适定性提供了关键范式。

不相关论文（基本信息）

1. 具有时滞的线性波动方程解的显式表示

Explicit representation of solutions to a linear wave equation with time delay

作者: Javad A. Asadzade, Jasarat J. Gasimov, Nazim I. Mahmudov, Ismail T. Huseynov

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05906v1

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

摘要: This paper develops an explicit spectral representation for solutions of a one-dimensional linear wave equation with a constant time delay. The model is considered on a bounded interval with non-homogeneous Dirichlet boundary data and a prescribed history function. To accommodate the loss of global smoothness in time caused by delay terms, solutions are understood in a \textit{stepwise classical sense}, allowing jump discontinuities in the second time derivative at multiples of the delay while maintaining continuity of the solution and its first time derivative. By combining separation of variables with Sturm-Liouville expansions, the delayed PDE is reduced to a family of scalar second-order delay differential equations. Using delay-dependent fundamental solutions, we derive closed-form representation formulas for the modal dynamics and reconstruct the PDE solution as a Fourier series. Convergence conditions guaranteeing uniform convergence and admissibility of termwise differentiation in space are established. A numerical example demonstrates the practical computation of truncated series solutions and their visualization.

2. 基于微局部Kakeya–Nikodym估计的算术本征函数的 $L^q$-范数界

$L^q$-norm bounds for arithmetic eigenfunctions via microlocal Kakeya-Nikodym estimate

作者: Jiaqi Hou, Xiaoqi Huang

类别: math.NT, math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05697v1

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

摘要: Let $X$ be a compact arithmetic hyperbolic surface, and let $ψ$ be an $L^2$-normalized Hecke-Maass form on $X$ with sufficiently large spectral parameter $λ$. We give a new proof to obtain some power saving for the global $L^6$-norm $\|ψ\|_{L^6(X)}\lesssim_ελ^{\frac{1}{6}-\frac{1}{144}+ε}$ over the local bound $\|ψ\|_{L^6(X)}\lesssimλ^{\frac16}$ of Sogge. Using the $L^\infty$-norm bound $\|ψ\|_{L^\infty(X)}\lesssim_ελ^{\frac{5}{12}+ε}$ of Iwaniec and Sarnak and harmonic analysis tools, we reduce the $L^6$-norm problem to a microlocal $L^6$ Kakeya-Nikodym estimate for $ψ$. Finally, we establish an improved microlocal $L^6$ Kakeya-Nikodym estimate via arithmetic amplification developed by Iwaniec and Sarnak.

3. 基于Hörmander向量场构造的高阶齐次 hypoelliptic 算子的基本解

Fundamental solution for higher order homogeneous hypoelliptic operators structured on Hörmander vector fields

作者: Stefano Biagi, Marco Bramanti

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05647v1

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

摘要: We introduce and study a new class of higher order differential operators defined on $\mathbb{R}^{n}$, which are built with Hörmander vector fields, homogeneous w.r.t. a family of dilations (but not left invariant w.r.t. any structure of Lie group) and have a structure such that a suitably lifted version of the operator is hypoelliptic. We call these operators ‘‘generalized Rockland operators’’. We prove that these operators are themselves hypoelliptic and, under a natural condition on the homogeneity degree, possess a global fundamental solution $Γ\left( x,y\right) $ which is jointly homogeneous in $\left( x,y\right) $ and satisfies sharp pointwise estimates. Our theory can be applied also to some higher order heat-type operators and their fundamental solutions.

4. 具有弱温度依赖参数的热黏弹性演化问题的大时间存在性

Large time existence in a thermoviscoelastic evolution problem with mildly temperature-dependent parameters

作者: Felix Meyer

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05640v1

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

摘要: We consider \begin{align*} \label{HS} \left{ \begin{array}{l} u_{tt} = (γ(Θ) u_{xt})x + a (γ(Θ) u_x)x +(f(Θ))x, \[1mm] Θ_t = DΘ{xx} + Γ(Θ) u{xt}^2 + F(Θ) u{xt}, \end{array}\right. \qquad \qquad (\star) \end{align*} under Neumann boundary conditions for $u$ and Dirichlet boundary conditions for $Θ$ in a bounded interval $Ω\subset\mathbb{R}$. \abs This model is a generalization of the classical system for the description of strain and temperature evolution in a thermo-viscoelastic material following a Kelvin-Voigt material law, in which $γ\equiv Γ$ and $f\equiv F$. Different variations of this model have already been analyzed in the past and the present study draws upon a known result concerning the existence of classical solutions, which are local in time, for suitably smooth initial data, arbitrary $a>0$, $D>0$ and $γ,f\in C^2([0,\infty))$ as well as $Γ,F\in C^1([0,\infty))$ with $γ>0,Γ\ge0$ and $F(0)=0$. Our work focuses on proving that existence times for classical solutions can be arbitrarily large, assuming sublinear temperature dependencies of $γ$ and $f$, and further $|F(s)|\le C_F(1+s)^α$ for some $C_F>0$ and $α\in(0,1)$. In particular, for any given $T_\star$, initial mass $M$ and $0<\underlineγ<\overlineγ$, there exists a constant $δ_\star(M,T_\star,a,D, Ω, \underlineγ, \overlineγ,C_F,α)>0$, such that if

$$ \underlineγ\leγ\le \overlineγ\quad\mbox{ and }\quad 0\le Γ\le \overlineγ\quad \mbox{ as well as } \quad\|γ'\|_{L^\infty([0,\infty))}\le δ_\star \quad \mbox{ and }\quad \|f'\|_{L^\infty([0,\infty))}\le δ_\star $$

hold, the maximal existence time of the classical solution to $(\star)$ surpasses $T_\star$.

5. 基于相容性映射不动点的非局部Fokker–Planck方程数值定态解

Numerical stationary states for nonlocal Fokker-Planck equations via fixed points of consistency maps

作者: José A. Carrillo, Yurij Salmaniw, Antonio León Villares

类别: math.NA, math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05632v1

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

摘要: We propose a fixed-point-based numerical framework for computing stationary states of nonlocal Fokker-Planck-type equations. Instead of discretising the differential operators directly, we reformulate the stationary problem as a nonlinear fixed-point map built from the original PDE and its nonlocal interaction terms, and solve the resulting finite-dimensional problem with a matrix-free Newton-Krylov method. We compare implementations using the analytic Frechet derivative of this map with a simple central-difference approximation. Because the method does not rely on time evolution, it is agnostic to dynamical stability and can detect both stable and unstable stationary states. Its accuracy is determined mainly by the numerical treatment of convolutions and quadrature, rather than by differentiation stencils. We apply the approach to three model problems with linear diffusion, use existing analytical results to verify the outputs, and reproduce known bifurcation diagrams, as well as new bifurcation behaviour not previously observed in this kind of problem.

6. 具有温度依赖参数的一维热粘弹性模型的整体光滑解

Global smooth solutions in a one-dimensional thermoviscoelastic model with temperature-dependent paramaters

作者: Felix Meyer

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05621v1

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

摘要: This manuscript is concerned with the system \begin{align*} \left{ \begin{array}{l} u_{tt} = (γ(Θ) u_{xt})x + (a(x,t) u_x)x +(f(Θ))x, \[1mm] Θ_t = DΘ{xx} + γ(Θ) u{xt}^2 + f(Θ) u{xt}, \end{array} \right. \end{align*} which is used to describe thermoviscoelastic developments in one-dimensional Kelvin-Voigt materials. \abs It is assumed that $a,γ$ and $f$ are sufficiently smooth functions that satisfy

$$ c_γ<γ(ζ)<C_γ, \quad γ''(ζ) \le 0,\quad f(0)=0, \quad |f'(ζ)|\le C_f \quad \mbox{ and } |f(ζ)|\le C_f(1+ζ)^α\quad \mbox{ for all }ζ\ge 0 $$

and some positive constants $c_γ,C_γ,C_f>0$ and $α\in (0,5/6)$. Under these conditions, this study then establishes a result on the existence of global classical solutions for sufficiently smooth but arbitrarily large initial data.

7. 基于磨光核构造的Carnot环面上线性抛物型方程的正则性结果

Regularity results for linear parabolic equations on Carnot tori via mollifier kernel construction

作者: Yiming Jiang, Yawei Wei, Yiyun Yang

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05498v1

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

摘要: This paper first proves the existence, uniqueness and regularity of the solution to a class of linear backward parabolic equations on Carnot tori, namely the periodic linear parabolic equation on Carnot groups. Such groups are non-commutative and typical examples of sub-Riemannian manifolds. Moreover, we apply the results for this equation to its dual equation (i.e., the Fokker-Planck-Kolmogorov equation in the general form), and derive the existence, uniqueness and regularity of its weak solution. To obtain the regularity results for solutions to the linear parabolic equation and its dual equation, firstly, we construct several families of mollifiers adapted respectively to the Hörmander vector fields generating Carnot groups, Carnot tori and dual spaces of non-isotropic Hölder spaces; secondly, we use the theory of singular integral operators to establish stronger a priori regularity for the solutions.

8. 完全匹配层方法在热弹性波散射问题中的收敛性分析

Convergence of the PML method for thermoelastic wave scattering problems

作者: Qianyuan Yin, Changkun Wei, Bo Zhang

类别: math.AP, math.NA

发布日期: 2026-02-05

arXiv ID: 2602.05497v1

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

摘要: This paper is concerned with the thermoelastic obstacle scattering problem in three dimensions. A uniaxial perfectly matched layer (PML) method is firstly introduced to truncate the unbounded scattering problem, leading to a truncated PML problem in a bounded domain. Under certain constraints on model parameters, the well-posedness for the truncated PML problem is then proved except possibly for a discrete set of frequencies, based on the analytic Fredholm theory. Moreover, the exponential convergence of the uniaxial PML method is established in terms of the thickness and absorbing parameters of PML layer. The proof is based on the PML extension technique and the exponential decay properties of the modified fundamental solution. As far as we know, this is the first convergence result of the PML method for the time-harmonic thermoelastic scattering problem.

9. 关于电阻猜想

On the Resistance Conjecture

作者: Sylvester Eriksson-Bique

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

发布日期: 2026-02-05

arXiv ID: 2602.05477v1

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

摘要: We give an affirmative answer to the resistance conjecture on characterization of parabolic Harnack inequalities in terms of volume doubling, upper capacity bounds and a Poincaré inequalities. The key step is to show that these three assumptions imply the so called cutoff Sobolev inequality, an important inequality in the study of anomalous diffusions, Dirichlet forms and re-scaled energies in fractals. This implication is shown in the general setting of $p$-Dirichlet Spaces introduced by the author and Murugan, and thus a unified treatment becomes possible to proving Harnack inequalities and stability phenomena in both analysis on metric spaces and fractals and for graphs and manifolds for all exponents $p\in (1,\infty)$. As an application, we also show that a Dirichlet space satisfying volume doubling, Poincaré and upper capacity bounds has finite martingale dimension and admits a type of differential structure similar to the work of Cheeger. In the course of the proof, we establish methods of extension and characterizations of Sobolev functions by Poincaré-inequalities, and extend the methods of Jones and Koskela to the general setting of $p$-Dirichlet spaces.

10. Lyapunov区域之间调和拟共形映射的Lipschitz正则性

Lipschitz regularity of harmonic quasiconformal maps between Lyapunov domains in $\mathbb{R}^n$

作者: Anton Gjokaj, David Kalaj

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05436v1

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

摘要: We prove that every sense-preserving harmonic $K$–quasiconformal homeomorphism $f\colon D\toΩ$ between Lyapunov domains (equivalently, bounded $C^{1,α}$ domains) in $\mathbb{R}^n$, $α\in(0,1]$, is globally Lipschitz on $\overline D$. The argument is based on a boundary iteration scheme: an initial Hölder modulus for the boundary trace (coming from quasiconformality) is improved via the $C^{1,α}$ graph representation of $\partialΩ$, yielding higher Hölder regularity for the normal component. This boundary gain is converted into a near-boundary gradient bound for harmonic functions through a basepoint boundary Hölder-to-gradient estimate obtained by flattening the boundary and using local harmonic-measure bounds. Quasiconformality then propagates the resulting control from one component to the full differential, and iteration gives boundedness of $|Df|$ up to the boundary. Along the way we briefly survey several standard tools from the theory of quasiconformal harmonic mappings (QCH), including boundary Hölder continuity, distortion of derivatives, and boundary-to-interior propagation principles that enter the iteration.

11. 具有新型自由边界条件的非局部传染病模型的动力学，第一部分：传播-消亡二分性

Dynamics of a nonlocal epidemic model with a new free boundary condition, part 1: Spreading-vanishing dichotomy

作者: Yao Chen, Yihong Du, Wan-Tong Li, Rong Wang

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05331v1

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

摘要: This paper investigates the long-time dynamics of a nonlocal epidemic model with free boundaries, where a pathogen with density $u(t,x)$ and the infected humans with density $v(t,x)$ evolve according to a reaction-diffusion system with nonlocal diffusion over a one dimensional interval $[g(t), h(t)]$, which represents the epidemic region expanding through its boundaries $x=g(t)$ and $x=h(t)$, known as free boundaries. Such a model with free boundary conditions based on those of Cao et al. \cite{fb27} was considered by several works. Inspired by recent works of Feng et al. \cite{fb20} and Long et al. \cite{fb5}, we propose a new free boundary condition, where the expansion rate of the epidemic region, determined by $h'(t)$ and $g'(t)$, is proportional to a linear combination of the outward flux of the pathogen $u$ through the range boundary (as in \cite{fb27}) and the weighted total population of infected individuals $v$ within the region (as in \cite{fb5}). We prove that the system under this new free boundary condition is well-posed, and its long-time dynamical behavior is characterized by a spreading-vanishing dichotomy. Moreover, we obtain sharp criteria for this dichotomy, including a sharp threshold in terms of the initial data $(u_0,v_0)$; and by studying a related eigenvalue problem, we also find a sharp threshold in terms of the diffusion rate, which complements related results in Nguyen and Vo \cite{fb7}. This is Part $1$ of a two part series. In Part $2$, we will determine the spreading speed of the model when spreading occurs, and for some typical classes of kernel functions, we will obtain the precise rates of accelerated spreading.

12. 时间周期介质中若干单稳自由边界问题的精确传播轮廓

Precise propagation profile for some monostable free boundary problems in time-periodic media

作者: Yihong Du, Zhuo Ma, Zhi-Cheng Wang

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05328v1

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

摘要: We consider reaction-diffusion equations of the form \begin{equation*} u_t - d u_{xx} = f(t,u), \quad t>0,\ \ x \in [g(t), h(t)], \end{equation*} where $f(t,u)$ is periodic in $t$ and monostable in $u$, and the interval $[g(t), h(t)]$ represents the one dimensional population range of a species with density $u(t,x)$ at time $t$ and spatial location $x$. The free boundaries $x=g(t)$ and $x=h(t)$ evolve subject to a ``preferred population density" condition at the habitat edges. Analogous to the traveling wave solutions in the corresponding Cauchy problem, semi-wave solutions play a fundamental role in understanding the propagation phenomena governed by the free boundary problem here. But in contrast to the Cauchy problem, where the KPP condition plays a subtle role in the precise approximation of its solution (with compactly supported initial function) by the traveling wave solution with minimal speed, here we prove the existence and uniqueness of a semi-wave in a general monostable setting, and obtain a precise description of the convergence of the solution toward the semi-wave as time goes to infinity, where the KPP condition plays no special role. Previously, such a sharp result was proved for a free boundary model only when $f$ is autonomous ($f=f(u)$, see \cite{D} or \cite{DL15} for a related free boundary model), or a less precise result was obtained in the time-periodic case under an extra strong KPP condition on $f$ (see \cite{MDW}, or \cite{DGP} for a related free boundary model). This work appears to be the first to prove the sharp convergence result for a general monostable free boundary problem in a heterogeneous environment, and we believe the methods developed here should have applications to related free boundary problems in heterogeneous media with nonlinearities more general than those of KPP type.

13. 随机时间分数阶慢扩散与快扩散方程的强局部非确定性

Strong local nondeterminism for stochastic time-fractional slow and fast diffusion equations

作者: Le Chen, Cheuk Yin Lee, Panqiu Xia

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

发布日期: 2026-02-05

arXiv ID: 2602.05317v1

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

摘要: We study a class of stochastic time-fractional equations on $\mathbb{R}^d$ driven by a centered Gaussian noise, involving a Caputo time derivative of order $β>0$, a fractional (power) Laplacian of order $α>0$, and a Riemann-Liouville time integral of order $γ\ge0$ acting on the noise. The noise is fractional in time (index $H$) and Riesz-type in space (index $\ell$). We derive sharp Dalang-type necessary and sufficient conditions for the existence of a random field solution across almost full parameter range $(α,β,γ;H,\ell)$. Under the Dalang-type conditions, we prove sharp variance bounds for temporal and spatial increments, as well as strong local nondeterminism in time in several regimes (two-sided version for $β=1$ and for parts of the case $β=2$; one-sided version for $0<β<2$) and strong local nondeterminism in space for the whole range of parameters. As applications, we derive exact uniform and local moduli of continuity, Chung-type laws of the iterated logarithm, and quantitative bounds on small ball probabilities. Along the way, we obtain sharp asymptotics for the fundamental solution kernels at $0$ and $\infty$, which may be of independent interest.

14. 一类非线性耦合薛定谔系统的无穷多个新解

Infinitely many new solutions for a nonlinear coupled Schrödinger system

作者: Qingfang Wang, Mingxue Zhai

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05283v1

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

摘要: We revisit the following nonlinear Schrödinger system \begin{align*}\begin{cases} -ε^{2}Δu +P(x) u= μ_1 u^3 +βuv^2, &\text{in};\mathbb {R}^3,\ -ε^{2}Δv+Q(x) v= μ_2 v^3 +βu^2v, &\text{in};\mathbb{ R}^3, \end{cases} \end{align*} where $ε$ is a positive parameter, $P(x),\,Q(x)$ are the potential functions, $μ_1>0$, $μ_2>0$ and $β\in\mathbb R$ is a coupling constant. Employing the finite dimensional reduction method, we prove that there are new kind of synchronized and segregated solutions, which concentrate both in a bounded domain and near infinity, and present a special structure. Moreover, by applying the local Pohozaev identities and some point-wise estimates of the errors, we prove that the new kind of synchronized solutions are non-degenerate, which is of great interest independently. One of the main difficulties of Schrödinger system come from the interspecies interaction between the components, which never appear in the study of single equation. Secondly, prior to the construction of new solutions, we shall verify the non-degeneracy of the solutions established in [Peng-Pi, Discrete Contin. Dyn. Syst., 2016] for the Schrödinger systems.

15. 高维相互作用粒子采样的Radon–Wasserstein梯度流

Radon–Wasserstein Gradient Flows for Interacting-Particle Sampling in High Dimensions

作者: Elias Hess-Childs, Dejan Slepčev, Lantian Xu

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

发布日期: 2026-02-05

arXiv ID: 2602.05227v1

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

摘要: Gradient flows of the Kullback–Leibler (KL) divergence, such as the Fokker–Planck equation and Stein Variational Gradient Descent, evolve a distribution toward a target density known only up to a normalizing constant. We introduce new gradient flows of the KL divergence with a remarkable combination of properties: they admit accurate interacting-particle approximations in high dimensions, and the per-step cost scales linearly in both the number of particles and the dimension. These gradient flows are based on new transportation-based Riemannian geometries on the space of probability measures: the Radon–Wasserstein geometry and the related Regularized Radon–Wasserstein (RRW) geometry. We define these geometries using the Radon transform so that the gradient-flow velocities depend only on one-dimensional projections. This yields interacting-particle-based algorithms whose per-step cost follows from efficient Fast Fourier Transform-based evaluation of the required 1D convolutions. We additionally provide numerical experiments that study the performance of the proposed algorithms and compare convergence behavior and quantization. Finally, we prove some theoretical results including well-posedness of the flows and long-time convergence guarantees for the RRW flow.

16. 高阶Hardy-Sobolev-Maz’ya不等式极值函数的存在性与对称性

Existence and symmetry of extremals for the high order Hardy-Sobolev-Maz’ya inequalities

作者: Guozhen Lu, Chunxia Tao

类别: math.AP

发布日期: 2026-02-05

arXiv ID: 2602.05203v1

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

摘要: In this article, we establish the existence of an extremal function for the k-th order critical Hardy-Sobolev-Maz’ya (HSM) inequalities on the upper half space $\mathbb{R}^{n+1}_{+}$ when $k\ge 2$ and $n\geq 2k+2$:

$$ \int_{\mathbb{R}^{n}_{+}}|\nabla^{k}u|^2dx-\prod_{i=1}^{k}\frac{\left(2i-1\right)^2}{4}\int_{\mathbb{R}^{n}_{+}}\frac{u^2}{x_1^{2k}}dx\geq C_{n,k,\frac{2n}{n-2k}} \left(\int_{\mathbb{R}^{n}_{+}}|u|^{\frac{2n}{n-2k}}dx\right)^{\frac{n-2k}{n}}. $$

The analysis of this extremal problem is challenging due to the presence of the higher order derivatives, the lack of translation invariance, the inapplicability of rearrangement techniques on the upper half-space, and the presence of a Hardy singularity along the boundary. To overcome these difficulties, instead of directly considering the HSM inequality on the upper half space, we establish the existence of an extremal for its equivalent version: Poincaré-Sobolev inequality on the hyperbolic space. We develop a novel duality theory of the minimizing sequences, the concentration-compactness principle for radial functions in the hyperbolic setting, which combines with the Helgason-Fourier analysis and the Riesz rearrangement inequality on the hyperbolic space, to resolve the lack of compactness issue. As an application, we also obtain the existence of positive symmetric solutions for the high order Brezis-Nirenberg equation on the entire hyperbolic space associated with the GJMS operators $P_k$ (i.e., when $k\ge 2$):

$$ P_{k}\left(f\right)-αf=|f|^{p-2}f $$

at the critical situation $α=\prod\limits_{i=1}^{k}\frac{\left(2i-1\right)^2}{4}$ when either $2k+2\leq n$ and $p=\frac{2n}{n-2k}$ or $2k<n$ and $2<p<\frac{2n}{n-2k}$.

17. 随机初值下聚焦非线性薛定谔方程最大振幅解的Painlevé普适类

Painlevé Universality classes for the maximal amplitude solution of the Focusing Nonlinear Schrödinger Equation with randomness

作者: Aikaterini Gkogkou, Guido Mazzuca, Kenneth D. T-R McLaughlin

类别: math.AP, math-ph, nlin.SI

发布日期: 2026-02-04

arXiv ID: 2602.05101v1

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

摘要: We establish universality for extremal solutions of the focusing nonlinear Schrödinger equation. Extremal solutions are $N$-soliton solutions that achieve the theoretical maximal amplitude and diverge as $N \to \infty$. We consider extremal solutions with the discrete eigenvalues randomly drawn from sub-exponential distributions, and identify two distinct universality classes, determined by the macroscopic structure of the spectrum: the Painlevé–III rogue-wave solution, where the eigenvalues take the form $λ_j = v_j + i μ_j$, and the Painlevé–V rogue wave solution, where $λ_j = -ζ\, j + v_j + i μ_j$, with $0 < ζ< 1$. (In both cases, $μ_{j}$ and $v_{j}$ are subexponential random variables.) Universality can then be summarized as follows: independently of the specific distribution of the eigenvalues, the rescaled solutions converge locally to a deterministic profile governed by the Painlevé-III equation in the first regime, and the Painlevé-V equation in the second. These results demonstrate that the formation of Painlevé-type rogue waves is a universal phenomenon robust to randomness.

18. 洛伦兹力方程的 $S^1$-指标理论

$\bf{S^1}$-index theory for the Lorentz force equation

作者: Cristian Bereanu, Alexandru Pîrvuceanu

类别: math.AP, math.CA

发布日期: 2026-02-04

arXiv ID: 2602.05015v1

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

摘要: In this paper we prove that the $S^1$-invariance of the Poincaré action functional associated to the Lorentz force equation gives the existence of multiple critical points which are periodic solutions with a fixed period. To do this, we prove an abstract multiplicity result which is based upon the Lusternik-Schnirelman method with the $S^1$-index. The corresponding result in the context of the Fadell-Rabinowitz index is proved in Ekeland and Lasry (Ann. Math., 112 (1980)). The main feature of our abstract result is that it allows us to consider nonsmooth functionals satisfying only a weak compactness condition well adapted to the Poincaré functional.

19. 半线性波动方程中势能项与阻尼系数的重构

Reconstruction of potential and damping coefficients in a semi-linear wave equation

作者: Rahul Bhardwaj, Mandeep Kumar, Manmohan Vashisth

类别: math.AP

发布日期: 2026-02-04

arXiv ID: 2602.04822v1

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

摘要: In this article, we investigate an inverse problem for a semi-linear wave equation posed in $\mathbb{R}^{n+1}$, with $n \geq 2$. Our primary objective is to reconstruct the damping coefficient, the linear potential, and the nonlinear potential from the associated Dirichlet-to-Neumann map. The analysis is based on a higher-order linearization method. As a key step, we establish the existence of suitable asymptotic solutions, crucial for reconstructing the nonlinear potential. In addition, we also provide a detailed study of the corresponding forward problem.