生成时间: 2025年12月18日 03:58:54

重点关注论文数量: 7

了解领域论文数量: 4

不相关论文数量: 24

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

1. Estimates for the 2D Navier-Stokes equations: the effects of forcing

关于二维Navier-Stokes方程的估计：外力作用的影响

作者: Ritwik Mukherjee, John D. Gibbon, Dario Vincenzi

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

发布日期: 2025-12-17

arXiv ID: 2512.15188v1

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

详细分析

中文标题: 关于二维Navier-Stokes方程的估计：外力作用的影响

研究对象和背景

本文研究的是二维不可压缩Navier–Stokes方程在周期性边界条件下受外部强迫项影响的动力学行为，重点关注能量耗散率、涡度耗散率以及全局吸引子维数的上界估计，并探讨这些量如何依赖于雷诺数（Reynolds number, $Re$）以及强迫项的正则性。所考虑的方程组如下：

$$ \partial_t u + u \cdot \nabla u + \nabla p = \nu \Delta u + f, \quad \nabla \cdot u = 0, $$

定义在二维周期域 $\Omega = [0, \ell]^2$ 上，其中 $u(x,t)$ 是速度场，$p(x,t)$ 是压力，$\nu > 0$ 是运动粘性系数，$f(x)$ 是一个时间无关、无散度且均值为零的外部强迫力。

由于研究的是二维情形，系统具有两个二次守恒量（对强解而言）：动能和涡度（enstrophy）。这导致了与三维湍流截然不同的动力学现象——即 Kraichnan–Leith–Batchelor (KLB) 双级联理论 所预测的行为：存在从中间尺度向大尺度的 逆向能量级联（inverse energy cascade），同时伴随着从小尺度到更小尺度的 直接涡度级联（direct enstrophy cascade）。这一理论框架是理解二维湍流统计性质的核心。

为了分析系统的长期行为，作者引入了标量涡度 $\omega = e_3 \cdot (\nabla \times u)$，从而将原方程转化为涡度形式：

$$ \partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega + \phi, \quad \text{其中 } \phi = e_3 \cdot (\nabla \times f). $$

该系统的数学分析通常使用广义Grashof数 $Gr$ 作为控制参数，其定义基于强迫项的大小。例如，当 $f \in L^2(\Omega)$ 时，

$$ Gr = \frac{\ell^3 \|f\|_{L^2}}{\nu^2}. $$

然而，在实验测量和统计湍流理论中，普遍采用的是雷诺数 $Re$，它由流动自身的响应（如平均速度或能量耗散率）决定：

$$ Re = \frac{U \ell}{\nu}, \quad \text{其中 } U^2 = \langle \|u\|_{L^2}^2 \rangle. $$

因此，要实现数学结果与物理观测之间的有效比较，必须将传统的以 $Gr$ 表达的估计转换为以 $Re$ 为变量的形式。Doering 和 Foias 在1983年首次完成了这种转换，但前提是强迫项满足较高的正则性条件，例如 $f \in \dot{H}^2(\Omega)$，即其二阶导数平方可积。

本文的关键动机在于：已有数值模拟表明，即使在二维情况下，强迫项的正则性程度（regularity）也会显著影响湍流的现象学特征。例如，Fractal Grid 实验中使用的宽频带强迫会破坏经典级联所需的尺度分离。因此，作者将现有估计推广至一类更广泛的强迫函数空间——Sobolev空间 $\dot{H}^s(\Omega)$，其中 $-1 \leq s \leq 2$。

这里 $\dot{H}^s(\Omega)$ 是齐次Sobolev空间，定义为所有均值为零的函数构成的空间，其范数为：

$$ \|g\|_{\dot{H}^s}^2 = \sum_k (\ell k)^{2s} |\hat{g}_k|^2, $$

其中 $\hat{g}_k$ 是 $g$ 的傅里叶系数，$k = |k|$ 是波数模长。这个定义中的 $\ell$ 因子确保了范数具有正确的量纲。参数 $s$ 刻画了函数的光滑性：$s$ 越大，函数越光滑；$s = 0$ 对应 $L^2$ 范数，$s = 2$ 对应二阶导数平方可积的条件，而 $s = -1$ 是目前已知保证弱解存在的最小正则性要求。

具体地，强迫项被设定为 $f(x) = F \Phi(\ell^{-1}x)$，其中 $F$ 是幅值，$\Phi$ 是无量纲、零均值、无散度且满足 $\|\Phi\|_{\dot{H}^s} = 1$ 的“形状”函数。

文中关注的主要物理量包括：

  • 时间平均的能量耗散率：$\displaystyle \varepsilon = \nu \ell^{-2} \left\langle \|\nabla u\|_{L^2}^2 \right\rangle = \nu \left\langle \|\omega\|_{L^2}^2 \right\rangle$；
  • 时间平均的涡度耗散率：$\displaystyle \chi = \nu \ell^{-2} \left\langle \|\nabla \omega\|_{L^2}^2 \right\rangle = \nu \ell^{-4} \left\langle \|u\|_{\dot{H}^2}^2 \right\rangle$；
  • 能量耗散长度尺度 $\eta_\varepsilon = (\nu^3 / \varepsilon)^{1/4}$ 和涡度耗散长度尺度 $\eta_\chi = (\nu^3 / \chi)^{1/6}$；
  • 全局吸引子的分形维数 $d_f(\mathcal{A})$，用于量化系统的“自由度”数量。

特别值得注意的是，只有当强迫项足够规则（如 $f \in L^2$）时，二维Navier–Stokes系统才被证明存在紧致的全局吸引子 $\mathcal{A}$。对于更粗糙的强迫（如 $s < 0$），虽然弱解存在，但是否仍存在吸引子尚不明确。

综上所述，本文的研究对象是一个在周期区域上受 $\dot{H}^s$ 类强迫驱动的二维Navier–Stokes系统，目标是建立关于 $\varepsilon$、$\chi$ 和 $d_f(\mathcal{A})$ 的 $Re$-依赖型上界估计，并揭示随着强迫正则性参数 $s$ 的变化，系统呈现出三种不同的动力学机制。

主要定理或主要结果

本文并未以传统方式陈述一个单一的“定理”，而是通过一系列推导得出针对不同 $s$ 区间的渐近估计。其核心结论可总结如下：

主要结果：设 $f \in \dot{H}^s(\Omega)$，$-1 \leq s \leq 2$，且 $\|\Phi\|_{\dot{H}^s} = 1$。则对于高雷诺数极限 $Re \gg 1$，以下估计成立（忽略对数因子和形状相关的常数）：

  1. 当 $0 \leq s < 2$ 时：

    $$ \frac{\chi \ell^3}{U^3} \lesssim Re^{\frac{2}{2+s}}, \quad \frac{\varepsilon \ell}{U^3} \lesssim Re^{-\frac{s}{2+s}}; $$

  2. 当 $-1 \leq s < 0$ 时：

    $$ \frac{\varepsilon \ell}{U^3} \lesssim Re^{s+2}. $$

此外，吸引子维数的估计为：

$$ d_f(\mathcal{A}) \lesssim Re^{\frac{2(4+s)}{3(2+s)}} (1+\ln Re)^{1/3}, \quad \text{for } 0 \leq s < 2. $$

这些结果揭示了三个不同的动力学区域：

  • $s \geq 2$：经典双级联区域，$\chi$ 收敛，$\varepsilon \to 0$；
  • $0 < s < 2$：中间区域，$\chi$ 随 $Re$ 增长，$\varepsilon$ 减小但慢于幂律；
  • $-1 \leq s < 0$：高耗散区域，$\varepsilon$ 可随 $Re$ 增长，类似于三维非平方可积强迫的情形。

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

a) 核心分析工具：

本文综合运用了多种现代PDE分析技术，主要包括：

  • 能量方法与Gronwall引理：通过对涡度方程 $(\partial_t \omega + u\cdot\nabla\omega = \nu\Delta\omega + \phi)$ 乘以 $\omega$ 并积分，得到涡度平衡方程：

    $$ \frac{1}{2}\frac{d}{dt}\|\omega\|_{L^2}^2 + \nu \|\omega\|_{\dot{H}^1}^2 = \langle \omega, \phi \rangle, $$

    再利用Gronwall引理证明时间平均下时间导数项消失，从而导出 $\chi = \langle \omega, \phi \rangle$。

  • Parseval恒等式与傅里叶分析：利用傅里叶展开将内积 $\langle \omega, \phi \rangle$ 转化为波数空间求和：

    $$ \left| \int_\Omega \omega \phi \,dx \right| \leq \sum_k k^2 |\hat{u}_k| |\hat{f}_k|, $$

    进而结合Sobolev范数表达式进行估计。

  • 插值不等式（Interpolation Inequality）：这是处理不同Sobolev范数之间关系的关键。例如，对于 $0 \leq s \leq 2$，有：

    $$ \|u\|_{\dot{H}^{2-s}} \leq \|u\|_{\dot{H}^2}^{1 - s/2} \|u\|_{L^2}^{s/2}, $$

    该不等式通过Hölder不等式在傅里叶系数上直接证明。

  • Cauchy–Schwarz 不等式与时均处理：在估计时间平均量时，反复使用：

    $$ \langle AB \rangle \leq \langle A^2 \rangle^{1/2} \langle B^2 \rangle^{1/2}, $$

    以及Jensen不等式 $\langle X^\alpha \rangle \geq \langle X \rangle^\alpha$（当 $\alpha \in (0,1)$）来控制非线性组合的时间平均。

b) 关键技术难点与处理技巧：

  1. 从 $\dot{H}^s$ 强迫到 $Re$-依赖估计的转换：
    核心挑战是如何用 $Re$ 控制强迫幅值 $F$。为此，作者采用Doering–Foias策略：选取一个测试函数 $\psi$（如 $\Phi$ 的Galerkin截断），将Navier–Stokes方程与 $\psi$ 做内积并时间平均，得到：

    $$ F \left| \int \Phi \cdot \psi \,dx \right| \leq \|\nabla \psi\|_{L^\infty} U^2 + \nu \|\Delta \psi\|_{L^2} U, $$

    从而导出 $F \lesssim \ell^{-1} U^2 + \nu \ell^{-2} U$，即 $F \lesssim \nu \ell^{-2} Re (1 + Re^{-1})$。这一估计使得最终结果仅依赖于 $Re$ 和强迫形状。

  2. 适应不同正则性区间的分段处理：
    对于 $0 \leq s < 2$，利用插值连接 $\|u\|_{\dot{H}^{2-s}}$ 与 $\|u\|_{\dot{H}^2}$ 和 $\|u\|_{L^2}$，从而将 $\chi$ 与 $U$ 和 $\chi$ 自身联系起来，形成闭合不等式并求解。而对于 $-1 \leq s < 0$，无法获得 $\chi$ 的界，故转而直接分析速度方程，沿用Cheskidov等人在三维中的方法。

  3. 吸引子维数的估计：
    基于已有的几何不等式 $d_f(\mathcal{A}) \lesssim (\ell / \eta_\chi)^2 (1 + \ln(\ell / \eta_\chi))^{1/3}$，代入新得到的 $\eta_\chi \sim Re^{-3(4+s)/(6(2+s))}$ 即得结果。

与之前工作的比较

a) 突破与改进：

本文最显著的贡献是将Reynolds数依赖型估计从 $f \in \dot{H}^2$ 推广到 $f \in \dot{H}^s$，$-1 \leq s < 2$，填补了二维情形下强迫正则性影响研究的空白。

  • 相较于Alexakis–Doering [34] 和 Gibbon–Pavliotis [35] 的工作（要求 $s \geq 2$），本文极大地放宽了对强迫项正则性的假设。他们只覆盖了最光滑的一端，而本文系统研究了从 $s=-1$（最粗糙）到 $s=2$ 的整个谱系。

  • 与三维中Cheskidov–Doering–Petrov [16] 的结果类比，本文首次在二维中揭示了类似“相变”的结构：当 $s < 0$ 时，$\varepsilon$ 可随 $Re$ 增长，打破了传统认为二维湍流中 $\varepsilon$ 应趋于零的预期，显示出强迫粗糙性带来的本质改变。

  • 此外，本文统一了多个先前孤立的结果：当 $s=2$ 时恢复[34,35]的经典估计；当 $s=0$ 时给出新的中间行为；当 $s<0$ 时继承三维思想但应用于二维模型。

因此，本文不仅推广了已有理论，更重要的是揭示了强迫正则性 $s$ 作为一个连续参数，如何调控湍流耗散机制和系统复杂性，为连接数学分析与物理模拟提供了更精细的桥梁。

2. Continuous data assimilation for 2D stochastic Navier-Stokes equations

二维随机Navier-Stokes方程的连续数据同化

作者: Hakima Bessaih, Benedetta Ferrario, Oussama Landoulsi, Margherita Zanella

类别: math.PR, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15184v1

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

详细分析

中文标题: 二维随机Navier-Stokes方程的连续数据同化

研究对象和背景

本文研究的是二维不可压缩随机Navier-Stokes方程（2D stochastic incompressible Navier-Stokes equations）在连续数据同化（Continuous Data Assimilation, CDA）框架下的长期行为，特别是通过nudging算法（又称反馈控制或松弛法）将一个近似解 $U(t)$ 同步到真实解 $u(t)$ 的收敛性问题。该工作将经典的确定性CDA理论推广至包含加性噪声与乘性噪声的随机情形，建立了期望意义下乃至路径意义上的一致收敛结果。

所考虑的物理系统是定义在一个有界光滑区域 $D \subset \mathbb{R}^2$ 上的速度场 $u(t,x)$ 和压力场 $p(t,x)$，满足如下随机偏微分方程组：

对于乘性噪声的情形：

$$ \begin{cases} du(t) - \nu \Delta u(t)\,dt + (u(t) \cdot \nabla) u(t)\,dt + \nabla p\,dt = f\,dt + G(u(t))\,dW(t), \\ \mathrm{div}\, u(t) = 0, \\ u(0) = u_0, \end{cases} $$

而对于加性噪声的情形则为：

$$ \begin{cases} du(t) - \nu \Delta u(t)\,dt + (u(t) \cdot \nabla) u(t)\,dt + \nabla p\,dt = f\,dt + G\,dW(t), \\ \mathrm{div}\, u(t) = 0, \\ u(0) = u_0. \end{cases} $$

其中，$\nu > 0$ 是运动粘性系数，$f$ 是外力项（假设属于 $V^*$），$u_0 \in H$ 是初始速度，且边界条件为无滑移Dirichlet边界条件，即 $u|_{\partial D} = 0$。随机扰动由 $W(t)$ 表示，它是一个取值于某个可分Hilbert空间 $U$ 的柱Wiener过程（cylindrical Wiener process）。算子 $G$ 描述了噪声的结构：当 $G$ 依赖于解 $u$ 时称为乘性噪声；否则为加性噪声。

为了进行数据同化，作者引入了一个插值算子 $R_h: [H^1(D)]^2 \to [L^2(D)]^2$，用于模拟对真实状态 $u(t)$ 在粗空间尺度 $h > 0$ 下的部分观测。该算子满足如下逼近性质：

$$ \|\phi - R_h(\phi)\|_{[L^2(D)]^2} \leq c_0 h \|\phi\|_{[H^1(D)]^2}, \quad \forall \phi \in [H^1(D)]^2, $$

其中 $c_0 > 0$ 是常数。典型的例子包括基于有限体积元的空间平均或低频Fourier模态投影。

在此基础上，定义所谓的nudging方程（也称数据同化方程），其形式为：

$$ dU(t) + [\nu A U(t) + B(U(t), U(t))] dt = f\,dt + G(U(t))\,dW(t) - \mu \Pi R_h(U(t) - u(t))\,dt, $$

其中 $A = -\Pi \Delta$ 是Stokes算子，$\Pi$ 是Leray-Helmholtz投影，$B(u,v)$ 是由非线性对流项诱导的双线性算子，$\mu > 0$ 是nudging参数，控制反馈强度。

目标是分析误差 $r(t) := U(t) - u(t)$ 是否随时间趋于零。由于系统具有强非线性和随机性，直接求解困难，而CDA提供了一种仅利用稀疏观测信息重构全局动力学的方法。

函数空间方面，文章采用标准的Solenoidal空间设置：

  • $H := \overline{\{v \in [C^\infty_0(D)]^2 : \mathrm{div}\,v = 0\}}^{L^2}$，赋予范数 $\| \cdot \|_H$ 和内积 $(\cdot,\cdot)_H$；
  • $V := \overline{\{v \in [C^\infty_0(D)]^2 : \mathrm{div}\,v = 0\}}^{H^1}$，赋予范数 $\|v\|_V := \|\nabla v\|_{[L^2]^2}$。

由此形成变分三重嵌入：

$$ V \subset H \equiv H^* \subset V^*, $$

所有嵌入均为紧且稠密。Stokes算子 $A$ 是从 $V$ 到 $V^*$ 的同构，并具有离散谱 $\{\lambda_n\}_{n\in\mathbb{N}}$，满足 $0 < \lambda_1 \leq \lambda_2 \leq \cdots \to \infty$，对应的标准正交特征基为 $\{e_n\}$。Poincaré不等式成立：

$$ \lambda_1 \|u\|_H^2 \leq \|u\|_V^2, \quad \forall u \in V. $$

此外，非线性项通过三线性形式 $b(u,v,w)$ 定义：

$$ b(u,v,w) = \int_D (u \cdot \nabla)v \cdot w\,dx, $$

并诱导出双线性映射 $B: V \times V \to V^*$，满足关键的反对称性：

$$ \langle B(u,v), w \rangle_{V^*,V} = - \langle B(u,w), v \rangle_{V^*,V}, \quad \text{且} \quad \langle B(u,v), v \rangle = 0. $$

这一性质在能量估计中至关重要，使得非线性项在自身上的贡献消失。

噪声结构方面，文章对乘性噪声提出了三种不同增长条件的假设：

  1. 有界噪声（Bounded noise）：存在 $K_B > 0$ 使得 $\|G(u)\|_{L_{HS}(U,H)}^2 \leq K_B$；
  2. 次线性增长噪声（Sublinear growth）：存在 $K_S, \gamma \in (0,1)$ 使得 $\|G(u)\|_{L_{HS}(U,H)}^2 \leq K_S + K_S' \|u\|_{H^\gamma}^{2\gamma}$；
  3. 线性增长噪声（Linear growth）：存在 $K_L$ 使得 $\|G(u)\|_{L_{HS}(U,H)}^2 \leq K_L + K_L' \|u\|_H^2$。

这些假设反映了现实中不同强度的随机扰动建模需求，例如模型截断误差、测量噪声或环境波动等。值得注意的是，这三种情况之间存在蕴含关系：(L) ⇒ (S) ⇒ (B)，但对应的收敛结果却越来越弱，体现出噪声增长对同步能力的负面影响。

综上所述，本文的研究对象是一类带有随机强迫的二维Navier-Stokes系统，其数学挑战在于结合无限维随机分析、非线性PDE理论与控制论思想，研究如何在仅有部分空间观测的情况下，通过设计合适的反馈机制实现对真实轨道的有效逼近。这种设定在气象预报、海洋建模等实际应用中具有重要意义，因为真实系统总是受到未解析尺度和观测误差的影响，必须借助随机模型来更真实地刻画不确定性。

主要定理或主要结果

本文的核心成果体现在以下几个定理中，分别针对乘性噪声和加性噪声两种情形给出了数据同化系统的收敛性保证。

乘性噪声情形

定理 4.2 设 Assumption 2.1 成立，且插值算子 $R_h$ 满足 (2.5)。令 $u$ 和 $U$ 分别为方程 (4.1) 和 (4.2) 的解。

(i) 若满足有界噪声条件 (B)，且参数 $\mu, h$ 满足：

$$ \frac{C_B}{\nu^2} + \frac{L}{\nu} < \mu \leq \frac{\nu}{c_0^2 h^2}, $$

其中

$$ C_B := 2\left(K_B + \frac{1}{\nu} \|f\|_{V^*}^2\right), $$

则有：

$$ \mathbb{E} \|u(t) - U(t)\|_H^2 \to 0 \quad \text{指数快速}, \quad t \to +\infty. $$

(ii) 若满足次线性增长条件 (S)，且 $\mu, h$ 满足：

$$ \frac{C_S}{\nu^2} + \frac{L}{\nu} < \mu \leq \frac{\nu}{c_0^2 h^2}, $$

其中

$$ C_S := 2\left(1 + K_S + (1-\gamma)\left(\frac{K_S'}{\lambda_1^\gamma \nu^{1-\gamma}}\right)^{\frac{1}{1-\gamma}} + \frac{1}{\nu} \|f\|_{V^*}^2\right), $$

则对任意 $p > 0$，

$$ \mathbb{E} \|u(t) - U(t)\|_H^2 \to 0 \quad \text{以速率 } O(t^{-p}), \quad t \to +\infty. $$

(iii) 若满足线性增长条件 (L)，且进一步要求：

$$ \nu > \frac{6K_L}{\lambda_1}, $$

同时 $\mu, h$ 满足：

$$ \frac{C_L}{\nu \left(\nu - \frac{K_L}{\lambda_1}\right)} + \frac{L}{\nu} < \mu \leq \frac{\nu}{c_0^2 h^2}, $$

其中

$$ C_L := 2\left(1 + K_L + \frac{1}{\nu} \|f\|_{V^*}^2\right), $$

则对任意 $p \in \left(0, \frac{\nu \lambda_1}{4K_L} - \frac{1}{2}\right)$，

$$ \mathbb{E} \|u(t) - U(t)\|_H^2 \to 0 \quad \text{以速率 } O(t^{-p}), \quad t \to +\infty. $$

加性噪声情形

定理 5.2 在加性噪声下，若 $G \in L_{HS}(U,H)$，且 $\mu, h$ 满足：

$$ \mu > \frac{C_A}{\nu^2}, \quad \mu \leq \frac{\nu}{c_0^2 h^2}, $$

其中 $C_A = 2\left(\|G\|_{L_{HS}}^2 + \frac{1}{\nu} \|f\|_{V^*}^2\right)$，则

$$ \mathbb{E} \|u(t) - U(t)\|_H^2 \to 0 \quad \text{指数快速}. $$

定理 5.3 在相同条件下，还成立几乎处处路径收敛：

$$ \|u(t) - U(t)\|_H^2 \to 0 \quad \text{指数快速}, \quad \mathbb{P}\text{-a.s. as } t \to +\infty. $$

这些结果表明：只要 nudging 参数 $\mu$ 足够大（但不能过大以免破坏稳定性）、观测分辨率 $h$ 足够细，就能实现同步。然而，噪声类型显著影响收敛模式——乘性噪声只能在期望中收敛，且收敛速度随噪声增长而退化；而加性噪声允许更强的路径一致收敛。

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

本文证明的核心在于建立关于误差 $r(t) = U(t) - u(t)$ 的Foias-Prodi型估计，即控制 $\mathbb{E} \|r(t)\|_H^2$ 随时间衰减的行为。为此，作者综合运用了无穷维随机分析、先验能量估计、停时技巧以及Grönwall型不等式的推广。

首先，在乘性噪声情形下，误差过程 $r(t)$ 满足以下随机微分方程：

$$ dr(t) + [\nu A r(t) + B(U,U) - B(u,u)] dt = -\mu \Pi R_h r(t) dt + [G(U) - G(u)] dW(t). $$

对其应用Itô公式于泛函 $\|r(t)\|_H^2$，得到：

$$ d\|r(t)\|^2_H + 2\nu \|r(t)\|^2_V dt = -2\mu (r(t), \Pi R_h r(t))_H dt - 2\langle B(r,u), r \rangle dt + \|G(U) - G(u)\|_{L_{HS}}^2 dt + 2(G(U)-G(u))dW(t), r(t). $$

接下来的关键步骤是对各项进行逐项估计。利用插值算子的逼近性质和Young不等式，可以证明：

$$ -2\mu (r, \Pi R_h r)_H \leq -\mu \|r\|_H^2 + \mu c_0^2 h^2 \|r\|_V^2. $$

再结合Poincaré不等式 $\|r\|_V^2 \geq \lambda_1 \|r\|_H^2$，可得：

$$ -2\mu (r, \Pi R_h r)_H \leq -\mu (1 - c_0^2 h^2 \lambda_1) \|r\|_H^2 + \mu c_0^2 h^2 (\|r\|_V^2 - \lambda_1 \|r\|_H^2). $$

因此，若选择 $h$ 足够小使得 $c_0^2 h^2 \lambda_1 < 1$，则负反馈项主导。

对于非线性项 $\langle B(r,u), r \rangle$，利用经典的Navier-Stokes估计：

$$ |\langle B(r,u), r \rangle| \leq C \|r\|_H \|r\|_V \|u\|_V \leq \nu \|r\|_V^2 + \frac{C}{\nu} \|r\|_H^2 \|u\|_V^2, $$

从而将其吸收进耗散项并产生一个依赖于 $\|u(s)\|_V^2$ 的积分因子。

噪声项的处理是区分乘性和加性情形的关键。在乘性噪声下，扩散项 $[G(U)-G(u)]dW(t)$ 不为零，其二次变差贡献为 $\|G(U)-G(u)\|_{L_{HS}}^2 dt$。利用Lipschitz假设 (2.4)，有：

$$ \|G(U)-G(u)\|_{L_{HS}}^2 \leq L \|r\|_H^2. $$

而在期望中方差项消失，故只需控制漂移部分。

最终，取期望后得到：

$$ \frac{d}{dt} \mathbb{E} \|r(t)\|_H^2 + (\mu - \tilde{C}) \mathbb{E} \|r(t)\|_H^2 \leq \mathbb{E} \left[ \|r(t)\|_H^2 \left( \frac{C}{\nu} \|u(t)\|_V^2 - \kappa \right) \right], $$

其中 $\kappa = \mu - L/\nu - C_B/\nu^2$。此时需要对 $\|u(t)\|_V^2$ 做一致先验估计（见附录A），确保其时间平均有界。这是通过构造适当的Lyapunov泛函并使用停时 $\tau_M = \inf\{t: \|u(t)\|_V^2 > M\}$ 来实现的。一旦获得 $\sup_t \mathbb{E} \|u(t)\|_V^2 < \infty$，即可应用Grönwall引理得出指数或多项式衰减。

特别地，在线性增长噪声下，需更强的粘性条件 $\nu > 6K_L / \lambda_1$ 以保证耗散压倒噪声增长，否则无法闭合估计。

在加性噪声情形下，最大优势在于：由于 $G$ 不依赖于解，误差方程中的噪声项完全抵消，即：

$$ d(u - U) + \cdots = -\mu \Pi R_h (u - U) dt, $$

不含随机积分项。因此可以直接对 $\|r(t)\|_H^2$ 进行路径wise估计，无需取期望。结合一致有界的先验估计（附录B），即可直接推出路径指数收敛。

总之，核心技术工具包括：

  • Itô公式在Hilbert空间的应用；
  • Foias-Prodi型能量估计；
  • Grönwall不等式的随机版本；
  • 停时方法处理爆炸风险；
  • 插值算子与Poincaré不等式的协同使用；
  • 对非线性项的经典控制技术（如Agmon不等式在高阶估计中隐含使用）。

与之前工作的比较

本文相对于已有文献的主要突破体现在以下几个方面：

首先，它是首个系统地将AOT型连续数据同化理论推广至一般乘性噪声随机Navier-Stokes方程的工作。此前的结果多集中于确定性系统（如[2]）或仅考虑加性噪声（如[5],[7]）。本文不仅涵盖加性情形作为特例，更重要的是处理了更为现实但也更难的乘性噪声，揭示了噪声结构对同步机制的根本影响。

其次，本文根据噪声的增长性建立了分层收敛理论：有界噪声 ⇒ 指数收敛；次线性增长 ⇒ 任意阶多项式收敛；线性增长 ⇒ 受限范围内的多项式收敛。这种精细分类在以往工作中未曾出现，体现了对噪声效应的深刻理解。尤其是在线性增长情形下发现必须满足 $\nu > C K_L / \lambda_1$ 才能收敛，说明粘性耗散必须足够强以对抗噪声的不稳定性，这是全新的物理洞察。

第三，相比数值研究（如[24],[10]）仅提供仿真支持，本文提供了严格的数学证明，特别是在路径收敛方面，明确区分了“期望收敛”与“几乎必然收敛”的适用范围，澄清了概念混淆。

第四，技术上，作者成功将用于研究不变测度唯一性和遍历性的工具（如[18],[23]中的方法）迁移到数据同化问题中，展示了不同领域间的深刻联系。例如，对解的一致二阶矩估计正是遍历理论中的经典技术。

最后，本文的结果在极限情况下恢复了确定性理论（令噪声强度→0），验证了其一致性。同时，所得条件中的常数显式依赖于 $\nu, \lambda_1, \|f\|, K_\cdot$ 等物理参数，具有明确的应用指导意义。

综上，本文不仅拓展了CDA的适用范围至随机系统，而且深化了对“随机强迫、粘性耗散与反馈控制三者平衡”的理解，为构建鲁棒的数据驱动预测模型提供了坚实的理论基础。

3. Nonlinear asymptotic stability and optimal decay rate around the three-dimensional Oseen vortex filament

三维Oseen涡丝周围的非线性渐近稳定性及最优衰减率

作者: Te Li, Ping Zhang, Yibin Zhang

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15040v1

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

详细分析

中文标题: 三维Oseen涡丝周围的非线性渐近稳定性及最优衰减率

研究对象和背景

本文研究的是在高雷诺数（high-Reynolds-number）条件下，三维不可压缩Navier–Stokes方程在Oseen涡丝（Oseen vortex filament）附近的长期动力学行为。该问题的核心在于理解强剪切流场中非轴对称扰动的演化机制，尤其是由涡旋拉伸（vortex stretching）、剪切（shearing）与混合（mixing）共同作用所引发的增强耗散（enhanced dissipation）现象。

考虑如下以涡度形式表示的三维不可压缩Navier–Stokes方程：

$$ \begin{cases} \partial_t \eta + v \cdot \nabla \eta - \eta \cdot \nabla v = \nu \Delta \eta, & (t,x,z) \in \mathbb{R}_+ \times \mathbb{R}^2 \times \mathbb{T}, \\ v = \nabla \times (-\Delta)^{-1} \eta, \end{cases} $$

其中 $\eta$ 表示涡度场，$v$ 表示速度场，$\nu > 0$ 是粘性系数。初始数据为集中于 $x=0$ 的涡丝：

$$ \eta(t=0) = \alpha \delta_{x=0}, $$

这里 $\delta_{x=0}$ 是支撑在直线 $\{x=0\} \subset \mathbb{R}^2 \times \mathbb{T}$ 上的发散自由向量值测度，$\alpha \in \mathbb{R}$ 是环量强度。系统定义在柱面区域 $\mathbb{R}^2 \times \mathbb{T}$ 上，其中 $\mathbb{T} = \mathbb{R}/(2\pi\mathbb{Z})$ 是单位圆周，表明在 $z$ 方向具有周期性边界条件。

高雷诺数情形对应于：

$$ \mathrm{Re} = |\alpha|/\nu \gg 1. $$

作者不失一般性地通过尺度变换将问题归一化为 $\nu = 1$、$\alpha \gg 1$，从而聚焦于大环量极限下的动力学。

该系统的自相似解——即3D Oseen涡丝——具有如下形式：

$$ \eta_G(t,x) = \alpha \begin{pmatrix} 0 \\ 0 \\ G(x/\sqrt{t}) \end{pmatrix}, \quad G(\xi) = (4\pi)^{-1} e^{-|\xi|^2/4}. $$

其对应的自相似速度场为：

$$ v_G(t,x) = \frac{\alpha}{\sqrt{t}} \begin{pmatrix} v_G(x/\sqrt{t}) \\ 0 \end{pmatrix}, \quad v_G(\xi) = \frac{\xi^\perp}{2\pi |\xi|^2} \left(1 - e^{-|\xi|^2/4}\right). $$

为了分析长期行为，引入自相似坐标：

$$ \tau = \log t, \quad \xi = x / \sqrt{t}, $$

并定义重标度后的涡度与速度场：

$$ H(\tau,\xi,z) = e^\tau \eta(e^\tau, e^{\tau/2} \xi, z), \quad V(\tau,\xi,z) = e^{\tau/2} v(e^\tau, e^{\tau/2} \xi, z). $$

在此坐标下，原方程转化为：

$$ \begin{cases} \partial_\tau H + V \cdot \nabla H - H \cdot \nabla V = (\mathcal{L} + e^\tau \partial_z^2) H, \\ V = \nabla \times (-\Delta)^{-1} H, \quad \mathrm{div}\, H = 0, \end{cases} $$

其中算子定义为：

$$ \mathcal{L} = \Delta_\xi + \frac{\xi}{2} \cdot \nabla_\xi + 1, \quad \nabla = (\partial_{\xi_1}, \partial_{\xi_2}, e^{\tau/2} \partial_z), \quad \Delta = \Delta_\xi + e^\tau \partial_z^2. $$

在该框架下，Oseen涡丝退化为平衡态 $(\alpha \Omega_G, \alpha U_G)$，其中：

$$ \Omega_G(\xi) = (0,0,G(\xi))^\top, \quad U_G(\xi) = (v_G(\xi), 0)^\top. $$

令 $H = \alpha \Omega_G + \Omega$，则扰动 $\Omega$ 满足非线性演化方程：

$$ \partial_\tau \Omega - (\mathcal{L} + e^\tau \partial_z^2)\Omega + \alpha \left( U_G \cdot \nabla \Omega + U \cdot \nabla \Omega_G - \Omega \cdot \nabla U_G - \Omega_G \cdot \nabla U \right) = -U \cdot \nabla \Omega + \Omega \cdot \nabla U, $$

其中 $U = \nabla \times (-\Delta)^{-1} \Omega$。

论文的关键挑战在于处理三维几何中的各向异性：横向（$\xi$）与轴向（$z$）的尺度不同，且存在涡旋拉伸项 $\eta \cdot \nabla v$，这在二维情况下不存在。此外，由于初始数据是奇异测度，需在适当的函数空间中建立适定性。

为此，作者采用各向异性的自相似函数空间。定义加权空间：

  • $L^2(m)$：配备范数 $\|w\|^2 = \int (1+|\xi|^2)^m |w(\xi)|^2 d\xi$；
  • $L^p(\infty)$：配备范数 $\|w\| = \|G^{-1/2} w\|_{L^p(d\xi)}$；
  • $Y = L^2(\infty)$，并定义子空间：
    $$ Y_0 = \{ w \in Y : \int w d\xi = 0 \}, \quad Y_1 = \{ w \in Y_0 : \int \xi w d\xi = 0 \}, \quad Y_2 = \{ w \in Y_1 : \int |\xi|^2 w d\xi = 0 \}. $$

进一步，考虑关于 $z$ 变量的Wiener代数结构：

$$ B_z(X) := \left\{ f : \|f\|_{B_z(X)} = \sum_{\zeta \in \mathbb{Z}} \| \hat{f}(\cdot, \zeta) \|_X < \infty \right\}, $$

其中 $\hat{f}$ 是关于 $z$ 的傅里叶变换。最终，扰动分量 $\Omega_z$ 属于空间 $Z \subset B_z(Y)$，满足零质量与零动量条件：

$$ \int_{\mathbb{R}^2 \times \mathbb{T}} \Omega_z d\xi dz = 0, \quad \int_{\mathbb{R}^2 \times \mathbb{T}} \xi \Omega_z d\xi dz = 0. $$

这些条件可通过调整环量 $\alpha$ 和平移不变性实现，因此不失一般性。

整个分析建立在自相似坐标系与各向异性函数空间的结合之上，使得能够精确刻画剪切诱导的混合效应及其对粘性耗散的放大作用。

主要定理或主要结果

定理 1.1（非线性渐近稳定性）

设初始扰动 $\Omega(0) \in B_z(Y)$ 满足：

$$ \mathrm{div}\, \Omega(0) = 0, \quad \int \Omega_z(0) d\xi dz = 0, \quad \int \xi \Omega_z(0) d\xi dz = 0. $$

则存在仅依赖于 $\alpha$ 的常数 $\kappa, \alpha_0 > 0$，以及 $\varepsilon(\alpha), C_\alpha > 0$，使得当 $\|\Omega(0)\|_{B_z(Y)} \leq \varepsilon(\alpha)$ 时，系统 (1.10) 存在唯一解 $\Omega(\tau) \in C([0,\infty); B_z(Y))$，且对所有 $\tau \geq 0$ 满足：

$$ (1+\tau)^{-1} e^{\tau/2} \|M_\tau \Omega_0(\tau)\|_{B_z(Y)} + e^{\kappa |\alpha|^{1/2} \tau} \|M_\tau \Omega_{\neq}(\tau)\|_{B_z(Y)} \leq C_\alpha \|\Omega(0)\|_{B_z(Y)}, $$

其中 $M_\tau = \exp\left( (e^{\tau/2} - 1)|D_z| \right)$ 是关于 $z$ 变量的傅里叶乘子，$\Omega_0$ 与 $\Omega_{\neq}$ 分别表示轴对称与非轴对称部分。

此定理表明：非轴对称扰动以指数速率 $e^{-\kappa |\alpha|^{1/2} \tau}$ 衰减，远快于热扩散的代数速率 $e^{-c\tau}$，体现了增强耗散机制。

定理 1.2（物理坐标的推论）

在原始变量下，若 $e^{|x|^2/8}(\eta - \alpha \eta_G)(1) \in B_z(L^2_x)$ 且满足相应的守恒律，则存在唯一解满足：

$$ (1+\log t)^{-1} t^{3/2 - 1/p} \|(\eta - \alpha \eta_G)_0(t)\|_{B_z(L^p_x)} + t^{\kappa |\alpha|^{1/2} + 1 - 1/p} \|(\eta - \alpha \eta_G)_{\neq}(t)\|_{B_z(L^p_x)} \leq C_\alpha \|e^{|x|^2/8}(\eta - \alpha \eta_G)(1)\|_{B_z(L^2_x)}. $$

特别地，非轴对称部分的衰减速率为 $t^{-\kappa |\alpha|^{1/2}}$，这是最优衰减率。

定理 1.4（谱下界最优估计）

定义线性化Oseen算子在正交补空间上的限制：

$$ \mathcal{L}_\perp - \alpha \Lambda_\perp = (\mathcal{L} - \alpha \Lambda)|_{(\mathrm{Ker}\,\Lambda)^\perp}, $$

其谱下界为：

$$ \Sigma(\alpha) = \inf \left\{ \mathrm{Re}\, z : z \in \sigma(-\mathcal{L}_\perp + \alpha \Lambda_\perp) \right\}. $$

则存在与 $\alpha$ 无关的常数 $C_1, C_2 > 0$，使得当 $|\alpha| \to \infty$ 时，

$$ C_1 |\alpha|^{1/2} \leq \Sigma(\alpha) \leq C_2 |\alpha|^{1/2}. $$

这一结果完全解决了 Gallay–Maekawa 在 [12] 中提出的猜想的第一部分，即 $\Sigma(\alpha) = O(|\alpha|^{1/2})$。

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

a) 核心分析工具

  1. 自相似坐标变换：将时间渐进行为转化为稳定态附近的演化问题。通过 $\tau = \log t$, $\xi = x/\sqrt{t}$ 将扩散主导的长时行为转化为谱分析问题。

  2. 各向异性函数空间 $B_z(Y)$：结合横向加权空间 $Y = L^2(\infty)$ 与轴向Wiener代数结构，有效控制高频振荡与局部奇异性。

  3. 谱理论与非自伴算子分析：研究非局部算子 $\mathcal{L} - \alpha \Lambda$ 的谱结构。其中 $\mathcal{L}$ 自伴，$\Lambda$ 斜自伴，其组合代表从自伴到非自伴的过渡。

  4. 复化谐振子逼近法：受 [9,21] 启发，构造一族复化的谐振子算子来逼近 $\mathcal{L}_k - \alpha \Lambda_k$（第 $k$ 个傅里叶模），通过复形变方法（complex scaling）识别达到最优谱界的本征值。

  5. Fourier分解与模态分析：将扰动按角度 $\theta$ 和轴向 $z$ 傅里叶展开：

    $$ \Omega(\xi,z) = \sum_{k \in \mathbb{Z}} \Omega_k(r) e^{ikz + in\theta}, $$

    并逐模分析算子 $(\mathcal{L} - \alpha \Lambda - ik e^\tau \partial_z)$ 的谱性质。

  6. 相对紧摄动理论：利用 $\alpha \Lambda$ 相对于 $\mathcal{L}$ 的相对紧性，证明谱界估计可传递至非线性系统。

b) 技术性技巧

  1. 显式构造极值谱点：作者不仅证明谱下界 $\Sigma(\alpha) \gtrsim |\alpha|^{1/2}$，更显式构造一个谱点其真实部恰好位于 $-O(|\alpha|^{1/2})$，从而确立最优性。

  2. 处理非局部项 $\Lambda$：$\Lambda \omega = v_G \cdot \nabla \omega + u \cdot \nabla G$ 包含非局部速度 $u = \nabla^\perp \Delta^{-1} \omega$，难以直接估计。作者通过将 $\nabla^\perp \Delta^{-1}$ 视为奇异积分算子，并在 $Y$ 空间中建立其有界性来克服。

  3. 各向异性能量估计：设计依赖于 $\tau$ 的加权能量：

    $$ \mathcal{E}(\tau) = \|\Omega_0(\tau)\|_{B_z(Y)}^2 + e^{2\kappa |\alpha|^{1/2} \tau} \|\Omega_{\neq}(\tau)\|_{B_z(Y)}^2, $$

    利用线性算子的谱间隙控制非线性项。

  4. 非线性项的小初值控制：当 $\|\Omega(0)\|_{B_z(Y)} \leq \varepsilon(\alpha)$ 时，非线性项 $U \cdot \nabla \Omega$ 可被吸收进线性衰减项，保证全局存在性。

  5. 模态分离与渐近匹配：对每个 $k \neq 0$ 模，分析其主导行为，发现最大增长率出现在 $k \sim |\alpha|^{1/2}$ 附近，从而导出整体最优衰减率。

与之前工作的比较

a) 突破与改进

  1. 首次证明三维Oseen涡丝的非线性稳定性具最优衰减率：
    虽然 [2] 建立了3D Oseen涡丝的全局稳定性，但未给出定量衰减率。本文首次证明非轴对称扰动以 $t^{-\kappa |\alpha|^{1/2}}$ 衰减，这是目前最精确的定量描述。

  2. 解决Gallay–Maekawa猜想：
    在 [12] 中提出 $\Sigma(\alpha) = O(|\alpha|^{1/2})$、$\Psi(\alpha) = O(|\alpha|^{1/3})$ 的猜想。本文结合 [25] 的结果，完全证实第一部分，即谱下界确为 $|\alpha|^{1/2}$ 阶，而非更低阶如 $|\alpha|^{1/3}$。

  3. 超越二维类比：
    在二维Lamb–Oseen涡中，Gallay [11] 利用伪谱界得到 $e^{-c|\alpha|^{1/3}\tau}$ 衰减。他猜测可能改进至 $e^{-\kappa|\alpha|^{1/2}\tau}$。本文不仅实现此改进，还证明其最优性，表明三维剪切混合机制比二维更强。

  4. 直接处理非局部算子的技术突破：
    多数谱分析针对微分算子，而 $\Lambda$ 含非局部项 $\nabla^\perp \Delta^{-1}$。本文发展了一套适用于此类非局部非自伴算子的复形变方法，为未来研究提供新工具。

  5. 放宽对初始数据的要求：
    相较于经典Sobolev空间，本文在 $B_z(Y)$ 中工作，允许初始扰动在横向仅具指数权重控制，无需高阶光滑性，更贴近物理实际。

综上所述，本文在三维涡丝稳定性问题上实现了从定性存在性到定量最优衰减的跨越，为高雷诺数流体中相干结构的持久性提供了严格的数学基础。

4. On global classical and weak solutions with arbitrary large initial data to the multi-dimensional viscous Saint-Venant system and compressible Navier-Stokes equations subject to the BD entropy condition under spherical symmetry

关于多维粘性Saint-Venant系统和可压缩Navier-Stokes方程在球对称下满足BD熵条件且具任意大的初始数据时的经典整体解与弱整体解的研究

作者: Xiangdi Huanga, Weili Meng, Xueyao Zhang

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15029v1

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

详细分析

中文标题: 关于多维粘性Saint-Venant系统和可压缩Navier-Stokes方程在球对称下满足BD熵条件且具任意大的初始数据时的经典整体解与弱整体解的研究

研究对象和背景

本文研究的是多维（二维和三维）可压缩等熵Navier-Stokes方程，其粘性系数依赖于密度，并满足所谓的 Bresch–Desjardins (BD) 熵条件。特别地，作者关注该类系统在球对称设定下的全局适定性问题，即是否存在定义在整个时间区间 $[0,\infty)$ 上的经典解或弱解，即使初始数据具有任意大的能量（不接近平衡态），且允许出现真空状态。

所研究的主控方程是如下形式的可压缩等熵Navier-Stokes系统：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) + \nabla P(\rho) - \mathrm{div}\left( \mu(\rho) D u \right) - \nabla\left( \lambda(\rho) \mathrm{div}\, u \right) = 0, \end{cases} $$

其中：

  • $\rho(x,t) \geq 0$ 是流体的密度；
  • $u(x,t) \in \mathbb{R}^N$ 是速度场；
  • $P(\rho) = \rho^\gamma$ ($\gamma > 1$) 是压力函数；
  • $D u = \frac{1}{2}(\nabla u + \nabla u^T)$ 是变形张量；
  • $\mu(\rho), \lambda(\rho)$ 分别为剪切粘性系数和体积粘性系数，满足物理条件：
    $$ \mu(\rho) \geq 0, \quad \mu(\rho) + N \lambda(\rho) \geq 0. $$

在本工作中，粘性系数被具体取为幂律形式并满足 BD熵关系：

$$ \mu(\rho) = \rho^\alpha, \quad \lambda(\rho) = (\alpha - 1)\rho^\alpha, $$

这恰好对应于经典的 BD熵结构：

$$ \lambda(\rho) = \rho \mu'(\rho) - \mu(\rho), $$

这一结构使得系统拥有一种特殊的数学熵——BD熵估计，可用于控制密度梯度的行为，在处理真空问题中起着关键作用。

值得注意的是，当 $\alpha = 1$ 时，上述模型退化为：

$$ \mu(\rho) = \rho, \quad \lambda(\rho) = 0, $$

此时系统被称为 粘性Saint-Venant系统，用于描述浅水波运动（shallow water equations）。因此，$\alpha=1$ 被称为 端点情形（endpoint case），它在物理上非常重要但在数学分析中更具挑战性。

研究区域 $\Omega$ 取为 $\mathbb{R}^N$ 中以原点为中心的球体 $B_R$，即：

$$ \Omega = B_R := \{ x \in \mathbb{R}^N : |x| < R \}, \quad N = 2 \text{ 或 } 3. $$

考虑径向对称（radially symmetric）解，即设：

$$ \rho(x,t) = \rho(r,t), \quad u(x,t) = u(r,t) \frac{x}{r}, \quad r = |x|. $$

在此对称假设下，原偏微分方程组可简化为关于 $(r,t)$ 的一维双曲-抛物耦合系统：

$$ \begin{cases} \rho_t + (\rho u)_r + \dfrac{N{-}1}{r} \rho u = 0, \\[1em] (\rho u)_t + (\rho u^2)_r + (\rho^\gamma)_r + \dfrac{N{-}1}{r} \rho u^2 - \dfrac{\alpha}{r^{N{-}1}} \left( r^{N{-}1} \rho^\alpha u_r \right)_r + \dfrac{N{-}1}{r} (\rho^\alpha)_r u = 0. \end{cases} $$

对应的初边值条件为：

$$ \rho(r,0) = \rho_0(r), \quad \rho u(r,0) = m_0(r), $$

以及边界条件：

$$ \rho u(0,t) = \rho u(R,t) = 0, $$

表示在中心与外边界处无质量通量。

由于系统包含奇性项 $\frac{N{-}1}{r}$，尤其在 $r=0$ 处可能出现奇异性，因此如何控制解在原点附近的正则性和避免真空形成成为核心难点之一。此外，若初始密度远离真空（即 $\rho_0 \geq \underline{\rho} > 0$），则希望证明解始终保持正密度；而若允许初始真空，则需研究真空是否会消失、何时消失等问题。

文中还指出，这类带有密度依赖粘性的系统长期以来面临严重困难：例如，对于二维浅水方程（$\alpha=1$, $\gamma=2$），即使在光滑、远离真空的大初值下，整体光滑解的存在性在过去一直是一个公开难题。本文的重要贡献正是在径向对称框架下解决了这一长期悬而未决的问题。

主要定理或主要结果

定理 2.1（远离真空的经典整体解，$N{-}1/N < \alpha < 1$）

设 $N = 2$ 或 $3$，参数 $(\alpha, \gamma)$ 满足：

  • 当 $N=2$：$0.5 < \alpha < 1$，$\gamma > 1$；
  • 当 $N=3$：$0.686 < \alpha < 1$，且 $1 < \gamma < 6\alpha - 3 + \frac{3-5\alpha}{2n_3(\alpha)}$，

其中 $n_3(\alpha)$ 是由 $\alpha_{3,\pm}(n)=\alpha$ 定义的函数（见后文定义）。假设初始数据 $(\rho_0, u_0)$ 为径向对称，满足：

$$ 0 < \underline{\rho} \leq \rho_0(r) \leq \overline{\rho} < \infty, \quad (\rho_0, u_0) \in H^3(\Omega), \quad u_0|_{\partial\Omega} = 0. $$

则初边值问题 (1)–(3) 存在唯一全局径向对称经典解 $(\rho,u)$，满足对任意 $T>0$，

$$ \rho^{-1}(x,t) + \rho(x,t) \leq C(T), $$

且解属于适当的空间（见 Definition 2.6），包括 $\rho \in C([0,T];H^3)$，$u \in C([0,T];H^3 \cap H_0^1)$ 等。

注记：此结果改进了 Zhang [65] 的工作，将二维中 $\alpha$ 的范围从 $[0.8,1)$ 扩展到最优物理范围 $(0.5,1)$，因为 $\mu + 2\lambda = \alpha \rho^\alpha \geq 0$ 要求 $\alpha \geq 0.5$。

定理 2.3（端点情形 $\alpha=1$ 的经典整体解，含浅水方程）

在二维情形下，取 $\alpha = 1$（即粘性Saint-Venant系统），$\gamma > 1$，初始数据同上且远离真空。则存在唯一的全局经典解 $(\rho,u)$，满足对所有 $(x,t) \in \Omega \times [0,\infty)$，

$$ \rho(x,t) \geq c > 0, $$

即密度有一致正下界，从而避免真空产生。

意义：这是首次在二维粘性浅水方程中，对任意大的光滑初值建立整体经典解的存在性，解决了该领域一个长期开放问题。

定理 2.4 与 2.5（允许真空的全局弱解，$\alpha \geq 1$）

设 $N=2$ 或 $3$，$\alpha \geq 1$，$\gamma > 1$，初始数据满足 $\rho_0 \in L^\infty$, $m_0 = \rho_0 u_0$，能量有限。则存在一个全局弱解 $(\rho,u)$，满足：

  • $\rho \in L^\infty(0,\infty; L^\infty(\Omega))$，
  • $u \in L^2(0,\infty; H^1(\Omega))$，
  • 解在分布意义下满足方程，
  • 并且存在新的时间无关 $L^4$ 估计：
    $$ \| \sqrt{\rho} u \|_{L^\infty_t L^4_x} + \| \nabla \rho^{\alpha/2} \|_{L^\infty_t L^2_x} \leq C. $$

此类弱解比以往工作（如 Guo-Jiu-Xin, Vasseur-Yu, Bresch-Vasseur-Yu）中的解具有更高的正则性。

定理 2.6（真空在有限时间内消失）

进一步地，对于上述 $\alpha \geq 1$ 构造的弱解，若初始总能量有限，则存在某个有限时间 $T_* > 0$，使得当 $t > T_*$ 时，$\rho(x,t) > 0$ 几乎处处成立。即：真空状态会在有限时间内完全消失。

这是首次将一维中已知的“真空消失”现象推广到高维球对称情形。

定理 2.7（仅可能在原点处出现真空的弱解，$\alpha=1$）

当 $\alpha = 1$ 时，若初始密度在原点附近有一定奇性但整体有界，可通过构造环形逼近区域上的解并令内半径趋于零，得到一类全局弱解，其性质如下：

  • 任何可能的真空只能出现在原点附近；
  • 在非真空区域上，解实际上是强解；
  • 此类解适用于浅水方程模型。

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

a) 核心分析学工具

  1. BD熵估计（Bresch–Desjardins Entropy）
    本文最关键的数学工具是利用 BD 结构导出的新型能量型估计。通过引入所谓的“有效速度”（effective velocity）变量：

    $$ w = u + \nabla (\phi(\rho)), \quad \text{其中 } \phi'(\rho) = \frac{2\alpha}{\alpha} \rho^{\alpha-1}, $$

    即 $w = u + 2\nabla(\rho^{\alpha-1})$（当 $\alpha \ne 1$），可以构造一个新的泛函：

    $$ \int \frac{1}{2} \rho |w|^2 dx + \int \frac{\rho^\gamma}{\gamma-1} dx, $$

    其时间导数可被控制，从而获得对 $\nabla (\rho^{\alpha/2})$ 和 $\mathrm{div}\,u$ 的先验估计。这种估计独立于时间，是控制高维临界嵌入的关键。

  2. 临界 Sobolev 嵌入与二维特殊性
    在二维情形下，有紧嵌入 $H^1 \hookrightarrow L^p$ 对所有 $p < \infty$ 成立，但在 $p=\infty$ 不成立。然而，结合 BD 熵提供的 $\|\nabla \rho^{\alpha/2}\|_{L^2}$ 控制，以及径向对称性带来的额外衰减（如 $|u(r)| \lesssim \|\nabla u\|_{L^2} \sqrt{|\log r|}$），可在 $N=2$ 下克服临界性障碍。这是为何主要结果仅在二维能覆盖端点 $\alpha=1$ 的根本原因。

  3. 参数延拓法（Parameter-continuation method）
    针对 $\alpha < 1$ 时在 $r=0$ 处可能出现的奇性，作者设计了一种新颖的连续性方法：固定其他参数，将 $\alpha$ 视为变量，从小扰动开始建立局部解，再通过先验估计将其延拓至整个区间 $(N{-}1)/N < \alpha < 1$。这种方法避免了直接处理奇点。

  4. 加权估计与 $r$-权重函数技巧
    为了控制中心附近的密度行为，作者引入了一个精心设计的 $r$-加权函数，例如考虑量：

    $$ \eta(r) = r^\delta, \quad \delta > 0, $$

    并估计 $\eta \rho$ 或 $\eta \nabla \rho$ 的范数，从而抑制原点处可能的集中或塌陷。这对处理 $\alpha < 1$ 时的低密度区域至关重要。

  5. 紧性方法与稳定性理论
    对于弱解的构造，采用标准的逼近-紧性框架：先添加人工粘性或正则化项构造近似解序列，再利用 Aubin–Lions 引理、De la Vallée Poussin 定理等获得紧性，最后取极限恢复原方程。关键在于获得逼近解的一致估计，特别是：

    • 时间无关的 $L^\infty$ 密度上界；
    • 空间梯度的均匀控制；
    • 时间导数的弱紧性。

  6. 新建立的时间无关 $L^4$ 估计
    本文一大突破是建立了：

    $$ \| \sqrt{\rho} u \|_{L^\infty_t L^4_x} + \| \nabla \rho^{\alpha/2} \|_{L^\infty_t L^2_x} \leq C(E_0), $$

    这个估计不随时间增长，显著强于以往文献中的结果，使得能够证明真空在有限时间内消失。

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

  1. 端点情形 $\alpha=1$ 的正密度下界估计
    当 $\alpha=1$，传统方法失效，因 $\rho^{\alpha-1} = 1$ 不再提供信息。作者转而估计所谓“有效速度”的 $L^\infty$ 范数：

    $$ v = u + \nabla \rho, $$

    利用其满足的演化方程和最大值原理，结合 Gronwall 不等式，推导出：

    $$ \|v\|_{L^\infty} \leq C \Rightarrow \inf \rho \geq c > 0. $$

    这是处理浅水方程的核心突破。

  2. 真空消失的定量分析
    通过构造 Lyapunov 型泛函：

    $$ \mathcal{V}(t) = \int \rho |\log \rho|^2 dx, $$

    并估计其耗散率，发现当 $\alpha \geq 1$ 时，粘性项会产生足够强的正则化效应，迫使 $\rho$ 在有限时间内远离零。这是首次在高维给出真空消失的定量时间估计。

  3. 原点奇性的规避策略（Theorem 2.7）
    对于 $\alpha=1$ 且允许在原点出现真空的情形，作者不在全球区域内直接求解，而是先在环形区域 $A_\iota = \{ \iota < r < R \}$ 上构造解 $(\rho_\iota, u_\iota)$，然后让 $\iota \to 0$。由于去掉原点，避免了奇性，同时利用一致估计保证极限存在。最终极限解在原点外光滑，仅可能在 $r=0$ 处退化。

  4. 割断函数（cut-off functions）的设计
    在 Theorem 2.7 的证明中，为获得逼近解的统一估计，作者引入一系列与 $\iota$ 无关的截断函数 $\chi_\iota(r)$，使其在 $r \sim \iota$ 区域平滑过渡，从而分离奇性区与正则区，实现估计的独立性。

与之前工作的比较

a) 突破与改进之处

  1. 首次解决二维粘性浅水方程的整体适定性问题
    尽管一维浅水方程已有丰富结果（Haspot, Cao-Li-Zhu 等），但二维情形由于缺乏有效的最大值原理和更强的非线性耦合，长期未能取得进展。本文首次在径向对称大初值下建立整体经典解，填补了该领域的空白。

  2. 显著放宽了对 $\alpha$ 的限制条件
    以往工作（如 Zhang [65]）要求 $\alpha \geq 0.8$ 才能控制密度下界，而本文通过更精细的能量结构分析，将下限降至物理极限 $0.5$，达到最优范围。这是因为当 $\alpha = 0.5$ 时，$\mu + 2\lambda = 0$，导致耗散机制崩溃，故无法进一步降低。

  3. 将真空消失现象从一维推广到高维
    Li-Li-Xin [41] 在一维中发现弱解会在有限时间后变为强解，即真空消失。本文首次在二维和三维中构造出具有相同性质的弱解，并给出严格证明，标志着该现象的高维可推广性。

  4. 构造出更高正则性的弱解
    相较于 Lions、Feireisl 类型的 renormalized 弱解（正则性较低），或 Vasseur-Yu 中的弱解（仅在 $L^2_t H^1_x$ 控制速度），本文构造的弱解具有：

    • 密度一致有界（$L^\infty$）；
    • 密度梯度受控（$\nabla \rho^{\alpha/2} \in L^\infty_t L^2_x$）；
    • 新的 $L^4$ 型估计； 因此正则性更强，更接近物理直观。

  5. 统一处理端点与非端点情形
    多数前期工作要么排除 $\alpha=1$，要么单独处理。本文通过区分 $\alpha<1$ 与 $\alpha=1$ 的不同估计策略，实现了对整个区间 $((N{-}1)/N, \infty)$ 的系统性覆盖，体现了方法的普适性。

  6. 无需小初值或渐近平衡假设
    经典的 Matsumura-Nishida 理论要求初值接近常状态，而本文结果适用于任意大的初值，只要远离真空或能量有限，大大扩展了适用范围。

综上所述，本文不仅解决了多个长期存在的公开问题，还在技术层面提出了多项创新工具，推动了密度依赖粘性可压缩流体理论的发展，尤其是在球对称框架下的全局适定性研究方面达到了新的高度。

5. Global well-posedness of the three-dimensional non-isentropic compressible magnetohydrodynamic equations under a scaling-invariant smallness condition

在尺度不变小性条件下三维非等熵可压缩磁流体动力学方程的整体适定性

作者: Lin Xu, Xin Zhong

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15021v1

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

详细分析

中文标题: 在尺度不变小性条件下三维非等熵可压缩磁流体动力学方程的整体适定性

研究对象和背景

本文研究的是三维空间 $\mathbb{R}^3$ 中带有远场真空条件的 非等熵可压缩磁流体动力学（MHD）方程 的柯西问题。这类方程描述了在磁场作用下，具有黏性、热传导性质的导电流体（如等离子体、液态金属或电解质）的运动规律。系统耦合了可压缩流体力学与电磁场理论，表现出高度非线性和强耦合特性。

该系统的控制方程如下：

$$ \begin{cases} \rho_t + \mathrm{div}(\rho u) = 0, \\ (\rho u)_t + \mathrm{div}(\rho u \otimes u) - \mu \Delta u - (\lambda + \mu)\nabla \mathrm{div}\, u + \nabla P = (\mathrm{curl}\, b) \times b, \\ c_v \rho (\theta_t + u \cdot \nabla \theta) + P \mathrm{div}\, u - \kappa \Delta \theta = 2\mu |D(u)|^2 + \lambda (\mathrm{div}\, u)^2 + \nu |\mathrm{curl}\, b|^2, \\ b_t - \mathrm{curl}(u \times b) = \nu \Delta b, \\ \mathrm{div}\, b = 0, \end{cases} \tag{1.1} $$

其中未知函数分别为：密度 $\rho \geq 0$、速度 $u \in \mathbb{R}^3$、绝对温度 $\theta \geq 0$ 和磁场 $b \in \mathbb{R}^3$。压力 $P$ 由理想气体状态方程给出：

$$ P = R \rho \theta, $$

其中 $R > 0$ 是气体常数。变形张量定义为：

$$ D(u) = \frac{1}{2} \left( \nabla u + (\nabla u)^\top \right), $$

$(\nabla u)^\top$ 表示梯度矩阵的转置。正常数 $c_v$、$\kappa$、$\nu$ 分别表示定容比热、热传导系数和磁扩散率。黏性系数 $\mu$ 和 $\lambda$ 满足物理约束：

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

考虑如下初始条件：

$$ (\rho, u, \theta, b)(x, 0) = (\rho_0, u_0, \theta_0, b_0)(x), \quad x \in \mathbb{R}^3, \tag{1.2} $$

以及远场行为（即真空条件）：

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

值得注意的是，本文所研究的情形允许初始密度在无穷远处趋于零（即“远场真空”），这使得解可能在全局范围内不远离平衡态，增加了分析难度。

作者关注的核心问题是：在何种条件下，上述系统存在唯一的整体强解？所谓强解，是指满足一定正则性的解，在几乎处处意义下方程成立，并且其时间导数和空间导数属于适当的 $L^p$ 空间。

根据文中的定义（Definition 1.1），强解 $(\rho, u, \theta, b)$ 在 $[0,T]$ 上需满足以下正则性要求（对某个 $q \in (3,6]$）：

  • $\rho \in C([0,T]; W^{1,q} \cap H^1),\quad \rho_t \in C([0,T]; L^2 \cap L^q)$；
  • $(u,\theta,b) \in C([0,T]; D^2 \cap D^1) \cap L^2([0,T]; D^{2,q})$；
  • $(u_t, \theta_t, b_t) \in L^2([0,T]; D^1),\quad (\sqrt{\rho}u_t, \sqrt{\rho}\theta_t, b_t) \in L^\infty([0,T]; L^2)$。

此外，还需满足兼容性条件（Compatibility Condition）以保证局部解的存在性：

$$ \begin{cases} -\mu \Delta u_0 - (\lambda + \mu)\nabla \mathrm{div}\, u_0 + \nabla(R \rho_0 \theta_0) - \mathrm{curl}\, b_0 \times b_0 = \sqrt{\rho_0} g_1, \\ -\kappa \Delta \theta_0 - 2\mu |D(u_0)|^2 - \lambda (\mathrm{div}\, u_0)^2 - \nu |\mathrm{curl}\, b_0|^2 = \sqrt{\rho_0} g_2, \end{cases} \tag{1.7} $$

其中 $g_1, g_2 \in L^2(\mathbb{R}^3)$。这一条件通常用于处理初值与算子之间的奇异性，确保解在 $t=0$ 处具有良好正则性。

已有大量文献研究此类系统。例如：

  • Pu–Guo [20] 在初值接近常数平衡态时建立了光滑解的整体存在性；
  • Liu–Zhong [19] 在总能量较小的情况下获得了整体强解及其代数衰减速率；
  • 文献 [16] 要求 $3\mu > \lambda$ 并对 $\|\rho_0\|_{L^\infty} + \|b_0\|_{L^3}$ 小加以限制；
  • Liu–Zhong [17] 提出了一个尺度不变量（scaling-invariant quantity），但仍然依赖于人工假设 $3\mu > \lambda$。

而本文的关键突破在于：完全移除了对黏性系数的人工限制 $3\mu > \lambda$，从而更符合物理实际（许多流体的体积黏度远大于剪切黏度）。同时，作者构造了一个新的尺度不变初值量，在保持不变性的同时适用于更广泛的黏性参数范围。

主要定理或主要结果

文章的主要结论是以下定理：

定理 1.1（整体适定性）
设 $q \in (3,6]$，初始数据 $(\rho_0 \geq 0, u_0, \theta_0 \geq 0, b_0)$ 满足正则性条件：

$$ \rho_0 \in H^1 \cap W^{1,q},\quad (\sqrt{\rho_0} u_0, \sqrt{\rho_0} \theta_0) \in L^2,\quad (u_0, \theta_0) \in D^1 \cap D^2,\quad b_0 \in H^2,\quad \mathrm{div}\, b_0 = 0, \tag{1.6} $$

以及兼容性条件 (1.7)，其中 $g_1, g_2 \in L^2(\mathbb{R}^3)$。

令 $\bar\rho \triangleq \|\rho_0\|_{L^\infty}$，并定义三个辅助量：

$$ S'_0 \triangleq \|\rho_0\|_{L^3} + (\bar\rho^2 + \bar\rho)\left( \|\sqrt{\rho_0} u_0\|_{L^2}^2 + \|b_0\|_{L^2}^2 \right), $$
$$ S''_0 \triangleq \|\nabla u_0\|_{L^2}^2 + (\bar\rho + 1)\|\sqrt{\rho_0} \theta_0\|_{L^2}^2 + \|\nabla b_0\|_{L^2}^2 + \|b_0\|_{L^4}^4, $$

再定义整体尺度不变量：

$$ N_0 \triangleq \left(1 + \bar\rho + \frac{1}{\bar\rho}\right) S'_0 S''_0. $$

则存在仅依赖于物理参数 $R, \mu, \lambda, \nu, \kappa, c_v$ 的正常数 $\varepsilon_0 > 0$，使得若

$$ N_0 \leq \varepsilon_0, $$

则柯西问题 (1.1)–(1.3) 存在唯一的一个全局强解 $(\rho, u, \theta, b)$。

特别地，此结果不要求 $3\mu > \lambda$。

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

a) 核心分析工具

  1. Bootstrap论证法（Bootstrapping Argument）
    作者采用经典的抽象Bootstrap原理（见Lemma 2.2），通过设定一个先验假设（hypothesis）$H(t)$，然后证明它蕴含更强的结论 $C(t)$，再利用连续性闭合论证区间。具体来说，文中设定了如下先验估计：

    $$ \rho \leq 2\bar\rho, \quad S_t + (\bar\rho^2 + \bar\rho) \sup_{s \in [0,t]} \|\rho(s)\|_{L^3} \int_0^t \|\nabla \theta\|_{L^2}^2 ds \leq 2\varepsilon, $$

    其中 $S_t$ 是随时间演化的能量型泛函。通过一系列估计验证在此假设下可推出更优估计（如密度上界降至 $3\bar\rho/2$，积分项减半），从而完成bootstrap闭合。

  2. 精细的能量估计（Delicate Energy Estimates）
    与传统基于总能量的方法不同，作者分别对动量方程乘以 $u$、温度方程乘以 $\theta$、磁场方程乘以 $b$ 进行测试，得到多个耦合的能量不等式。例如基本能量估计（Lemma 3.1）：

    $$ \|\sqrt{\rho} u(t)\|_{L^2}^2 + \|b(t)\|_{L^2}^2 + \int_0^T \left( \|\nabla u\|_{L^2}^2 + \|\nabla b\|_{L^2}^2 \right) dt \leq C\left( \|\sqrt{\rho_0} u_0\|_{L^2}^2 + \|b_0\|_{L^2}^2 \right) + C \sup_t \|\rho(t)\|_{L^3}^2 \int_0^T \|\nabla \theta\|_{L^2}^2 dt. $$

    此处出现了关键结构：$\|\rho\|_{L^3}$ 控制耗散项增长，成为后续估计的基础。

  3. Gagliardo–Nirenberg不等式
    多次使用该不等式进行中间插值，例如：

    $$ \|f\|_{L^6} \leq C \|\nabla f\|_{L^2}, \quad \|f\|_{L^\infty} \leq C \|f\|_{L^6}^{1/2} \|\nabla f\|_{L^6}^{1/2} \leq C \|\nabla f\|_{L^2}^{1/2} \|\nabla^2 f\|_{L^2}^{1/2}. $$

  4. 有效粘性通量（Effective Viscous Flux）技术
    定义：

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

    利用椭圆正则性（Lemma 2.4）将速度梯度控制在高阶范数中：

    $$ \|\nabla u\|_{L^6} \leq C\left( \|\rho\|_{L^\infty}^{1/2} \|\sqrt{\rho} \dot{u}\|_{L^2} + \|\rho\|_{L^\infty} \|\nabla \theta\|_{L^2} + \||b||\nabla b|\|_{L^2} \right). $$

  5. 特征线方法（Characteristic Method）
    用于控制密度 $L^3$ 范数的增长。通过对质量方程沿流线分析，结合积分形式的演化估计（如 Lemma 3.2），最终获得 $\|\rho(t)\|_{L^3}$ 的一致有界性。

b) 技术难点与创新技巧

  1. 克服 $3\mu > \lambda$ 条件的缺失
    原先工作 [17] 依赖此条件获取速度的 $L^4$ 加权估计。本文通过引入新结构：

    $$ (\bar\rho^2 + \bar\rho) \sup_{t} \|\rho(t)\|_{L^3} \int_0^T \|\nabla \theta\|_{L^2}^2 dt, $$

    并将其纳入bootstrap框架，成功避免对该黏性条件的需求。

  2. 磁场项的复杂耦合处理
    方程中出现 $(\mathrm{curl}\, b) \times b$ 项，导致动量方程右侧产生非局部、非线性力。作者通过细致分解（如 Lemma 3.4 中的 $K_2(t)$）分离出所有涉及 $b, \nabla b, \nabla^2 b$ 的交叉项，并借助Sobolev嵌入和Young不等式进行吸收。

  3. 温度方程作为主控方程
    以往多用总能量估计，但此处改用温度方程直接估计 $\|\nabla \theta\|_{L^2}$，从而绕开对 $E = |u|^2/2 + c_v \theta$ 的依赖，增强了灵活性。

  4. 构造完整的尺度不变量
    关键观察是：在缩放变换

    $$ \rho_\tau(x,t) = \rho(\tau x, \tau^2 t),\quad u_\tau = \tau u(\tau x, \tau^2 t),\quad \theta_\tau = \tau^2 \theta(\tau x, \tau^2 t),\quad b_\tau = \tau b(\tau x, \tau^2 t) $$

    下，文中定义的 $N_0$ 保持不变。这是唯一非平凡的尺度不变变换，说明 $N_0$ 是自然的临界量。

  5. 避免对初始数据施加额外小性
    相较于之前要求 $\|\rho_0\|_{L^1}, \|b_0\|_{L^2}$ 小，本文的小性只体现在组合量 $N_0$ 上，更具结构性和不变性。

与之前工作的比较

a) 改进之处

  1. 彻底移除人工黏性条件 $3\mu > \lambda$
    这是最显著的改进。原条件 $3\mu > \lambda$ 即 $\lambda < 3\mu$，等价于体积黏度 $\zeta = 2\mu/3 + \lambda$ 满足 $\zeta < 3\mu$，但在物理上许多材料（如水蒸气、聚合物熔体）的 $\zeta \gg \mu$。因此，移除此限制使结果更具普适性和物理合理性。

  2. 提出新的尺度不变初值量
    文献 [17] 使用的量为：

    $$ \hat{\rho} \left( \|\rho_0\|_{L^3} + \hat{\rho}^2 (\|\sqrt{\rho_0} u_0\|_{L^2}^2 + \|b_0\|_{L^2}^2) \right) \left( \|\nabla u_0\|_{L^2}^2 + \hat{\rho} (\|\sqrt{\rho_0} E_0\|_{L^2}^2 + \|\nabla b_0\|_{L^2}^2) \right), $$

    而本文构造的新量 $N_0$ 不仅保持尺度不变性，还通过因子 $(1 + \bar\rho + 1/\bar\rho)$ 实现对低密度与高密度情形的统一控制，更加稳健。

  3. 放宽对初始能量的要求
    不同于 [19] 要求“总能量小”，本文仅需某个特定组合量 $N_0$ 小，允许某些单项（如动能或磁场能）较大，只要整体平衡即可。

  4. 更强的解类：强解而非弱解
    结果给出的是强解，具有更高的正则性（如 $\nabla u \in L^\infty(L^2), \sqrt{\rho} u_t \in L^\infty(L^2)$），优于仅存在的弱解。

  5. 无需额外边界条件
    在全空间 $\mathbb{R}^3$ 中处理真空问题，无需像有界域那样施加Navier-slip或Dirichlet边界条件，更具挑战性。

综上所述，本文不仅在数学上推进了非等熵可压缩MHD系统的理论，也使其更贴近真实物理场景，为后续研究提供了新的分析范式。

6. On self-similar singular solutions to a vorticity stretching equation

关于涡量拉伸方程的自相似奇异解

作者: Dapeng Du, Jingyu Li, Xinyue Shi

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15017v1

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

详细分析

中文标题: 关于涡量拉伸方程的自相似奇异解

研究对象和背景

本文研究的是一个简化的二维非局部偏微分方程模型，用于模拟三维不可压缩欧拉方程中的涡量拉伸（vorticity stretching）机制。该机制是三维流体中可能产生有限时间奇点（blow-up）的关键物理过程之一，而在二维情形下由于涡量不被拉伸，通常不会发生奇点形成。因此，构造并分析能够体现涡量拉伸效应的低维模型，成为理解三维欧拉方程是否存在奇异解的重要途径。

论文所考虑的核心方程是一个非线性非局部演化方程：

$$ \omega_t = Z_{11} \omega \cdot \omega, $$

其中 $\omega = \omega(x,t)$ 是定义在 $\mathbb{R}^2$ 上的标量函数，代表某种简化后的涡量；而 $Z_{11}$ 是一个Calderón-Zygmund算子，具体形式为：

$$ Z_{11} = \partial_{11} \Delta^{-1}, $$

即先对函数作用拉普拉斯逆算子 $\Delta^{-1}$（相当于卷积以牛顿位势），再对其结果求关于第一个空间变量的二阶偏导数 $\partial_{11} = \frac{\partial^2}{\partial x_1^2}$。这个算子具有零阶齐次性，在傅里叶侧对应乘子：

$$ \mathcal{F}(Z_{11} f)(\lambda) = \frac{\lambda_1^2}{|\lambda|^2} \hat{f}(\lambda), $$

这表明 $Z_{11}$ 是一个典型的奇异积分算子，其符号属于 $L^\infty$ 但不光滑，在某些方向上退化——特别是当 $\lambda_1 = 0$ 时，符号为零，导致算子在垂直于 $x_1$-轴的方向上“失去强度”。这种退化性（degeneracy）是本文分析中的主要技术障碍。

该方程最早由 Constantin 和 Sun 在文献 [3] 中提出作为二维非局部模型，形式与本文一致：

$$ \omega_t = Z_{11} \omega \cdot \omega. $$

他们证明了对于一大类负初值存在全局正则解。然而近年来的研究，如 Bianchini 与 Elgindi [1] 的工作表明，该方程在适当初值下会发生有限时间爆破。本文在此基础上进一步深入，旨在构建一类特殊的自相似奇异解（self-similar singular solutions），并给出新的爆破准则。

所谓自相似解，是指其空间结构随时间缩放保持不变的解。这类解常出现在非线性演化方程的奇点形成过程中，被认为是奇点形成的“普适模式”。本文关注如下形式的自相似解：

$$ \omega(x,t) = \frac{1}{T - t} Q(x), $$

其中 $T > 0$ 是爆破时间，$Q(x)$ 是固定的剖面函数（profile function）。将此形式代入原方程后可得关于 $Q$ 的静态方程：

$$ Z_{11} Q \cdot Q = Q. $$

目标是证明此类解的存在性，并揭示其结构特征。

研究区域为全空间 $\mathbb{R}^2$，未施加任何边界条件，因此属于无界域上的整体问题。函数空间方面，作者主要工作在 $L^2(\mathbb{R}^2)$ 空间及其子集上。特别地，定理 1.1 中要求剖面函数 $Q \in L^2(\mathbb{R}^2)$，且满足支撑集限制：$Q|_{\mathbb{R}^2 \setminus A} = 0$，即 $Q$ 仅在某个有界可测集 $A \subset \mathbb{R}^2$ 上非零。此外还要求 $Z_{11} Q|_A = 1$，这意味着在集合 $A$ 内部，算子 $Z_{11}$ 作用于 $Q$ 后恒等于 1。

值得注意的是，尽管 $Q$ 属于 $L^2$，但由于 $Z_{11}$ 是零阶算子，$Z_{11} Q$ 也应在 $L^2$ 中，但此处要求其在 $A$ 上为常数 1，若 $A$ 有正测度，则意味着 $Z_{11} Q$ 在局部为光滑函数。结合 $Q$ 具有紧支撑这一事实，这就涉及到了函数与其像之间的局部-整体关系，从而引出了所谓的谱不确定性原理（spectral uncertainty principle）。

该原理本质上断言：一个函数及其傅里叶变换不能同时在相空间中高度局域化。更精确地说，如果一个函数在物理空间中有紧支撑，其傅里叶变换不可能也在某个紧集上有紧支撑，除非它是零函数。本文使用的是 Logvinenko-Sereda 定理的一个推论，它给出了定量版本的这一思想。

综上所述，本文的研究对象是一个具有退化奇异积分算子的非线性演化方程，目标是在全平面 $\mathbb{R}^2$ 上构造具有紧支撑剖面的自相似奇异解，并建立新的有限时间爆破结果。核心挑战在于克服 $Z_{11}$ 的退化性，而这通过引入谱不确定性原理得以实现。

主要定理或主要结果

定理 1.1（自相似奇异解的存在性）
设 $A$ 是 $\mathbb{R}^2$ 中的一个有界可测集，则方程

$$ \omega_t = Z_{11} \omega \cdot \omega $$

存在如下形式的自相似奇异解：

$$ \omega(x,t) = \frac{1}{T - t} Q(x), \quad T > 0, $$

其中剖面函数 $Q$ 满足：

  • $Q \in L^2(\mathbb{R}^2)$，
  • $Q|_{\mathbb{R}^2 \setminus A} = 0$（即 $Q$ 支撑在 $A$ 上），
  • $Z_{11} Q|_A = 1$（即在 $A$ 上 $Z_{11} Q \equiv 1$）。

定理 1.2（有限时间爆破结果）
设初始数据 $\omega_0$ 满足：

  • $\omega_0 \in C^\alpha(\mathbb{R}^2)$ 对某个 $\alpha > 0$，
  • $\omega_0$ 在 $\mathbb{R}^2$ 上紧支撑，
  • $\int_{\mathbb{R}^2} \omega_0 \, dx > 0$，

则对应于初值 $\omega(x,0) = \omega_0$ 的解会在有限时间内爆破。

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

a) 核心分析学工具

本文证明中的核心技术工具是 Logvinenko-Sereda 定理 及其推论，这是一种深刻的调和分析工具，属于谱不确定性原理的范畴。该定理定量描述了函数与其傅里叶变换之间在支撑集上的互斥性质。

定理 2.1（Logvinenko-Sereda）
设 $f \in L^2(\mathbb{R}^n)$，$\hat{f}$ 为其傅里叶变换。若 $\mathrm{supp}(\hat{f}) \subset A \Subset \mathbb{R}^n$，且 $B \subset \mathbb{R}^n$ 满足 $|B| > 0$，则存在仅依赖于 $A, B, n$ 的常数 $C > 0$，使得

$$ \|f\|_{L^2(B)} \geq C^{-1} \|f\|_{L^2(\mathbb{R}^n)}. $$

这一定理说明：如果一个函数的傅里叶变换被限制在一个紧集中（即频率受限），那么它在任意正测度集合上的 $L^2$ 范数不能太小。换句话说，频域局域化迫使函数在空间域上“扩散”。

由此可推出以下关键推论：

推论 2.1
设 $f \in L^2(\mathbb{R}^n)$，$\mathrm{supp}(f) \subset A \Subset \mathbb{R}^n$，$B \subset \mathbb{R}^n$ 且 $|B| > 0$，则存在常数 $C = C(A,B,n) > 0$，使得

$$ \|\hat{f}\|_{L^2(B)} \geq C^{-1} \|\hat{f}\|_{L^2(\mathbb{R}^n)}. $$

该推论指出：若函数本身具有紧支撑，则其傅里叶变换在任意正测度频率集上的能量占比是可控的，不会完全集中在零附近或其他孤立点上。

在本文中，这一工具被用于处理 $Z_{11}$ 算子的退化性。由于 $Z_{11}$ 的符号 $\frac{\lambda_1^2}{|\lambda|^2}$ 在 $\lambda_1 = 0$ 时为零，若 $\hat{f}$ 集中在 $\lambda_1=0$ 附近，则 $Z_{11} f$ 会非常小。但推论 2.1 保证：只要 $f$ 有紧支撑，$\hat{f}$ 就不可能完全集中在 $\lambda_1=0$ 平面上，从而 $Z_{11} f$ 仍保留足够的“非退化”成分。

基于此，作者构造了一个从 $L^2(A)$ 到自身的算子：

$$ L: L^2(A) \to L^2(A), \quad L\phi = (Z_{11} \tilde{\phi})|_A, $$

其中 $\tilde{\phi}$ 是 $\phi$ 在 $\mathbb{R}^2$ 上的零延拓。为了证明 $L$ 可逆，需验证其强制性（coercivity），即存在 $\delta > 0$ 使得

$$ \langle L\phi, \phi \rangle_{L^2(A)} \geq \delta \|\phi\|_{L^2(A)}^2. $$

计算内积：

$$ \langle L\phi, \phi \rangle = \int_{\mathbb{R}^2} Z_{11} \tilde{\phi} \cdot \tilde{\phi} \, dx = \int_{\mathbb{R}^2} \frac{\lambda_1^2}{|\lambda|^2} |\hat{\tilde{\phi}}(\lambda)|^2 \, d\lambda. $$

令 $\Omega_k = \left\{ \lambda \in \mathbb{R}^2 : \frac{\lambda_1^2}{|\lambda|^2} > \frac{1}{k^2 + 1} \right\}$，则在 $\Omega_k$ 外该符号小于等于 $\frac{1}{k^2+1}$。于是：

$$ \langle L\phi, \phi \rangle \geq \int_{\Omega_k} \frac{\lambda_1^2}{|\lambda|^2} |\hat{\tilde{\phi}}|^2 d\lambda \geq \frac{1}{k^2 + 1} \int_{\Omega_k} |\hat{\tilde{\phi}}|^2 d\lambda. $$

应用推论 2.1 到 $\tilde{\phi}$（其支撑在 $A$），存在 $C = C(A,\Omega_k)$ 使得：

$$ \int_{\Omega_k} |\hat{\tilde{\phi}}|^2 d\lambda \geq C^{-1} \int_{\mathbb{R}^2} |\hat{\tilde{\phi}}|^2 d\lambda = C^{-1} \|\tilde{\phi}\|_{L^2}^2 = C^{-1} \|\phi\|_{L^2(A)}^2. $$

因此

$$ \langle L\phi, \phi \rangle \geq \frac{1}{(k^2 + 1)C} \|\phi\|_{L^2(A)}^2, $$

即 $L$ 是强制的。由此可得 $L$ 是单射、逆算子有界，并最终通过闭值域定理证明 $L$ 是满射，从而存在 $Q = L^{-1} 1$，完成构造。

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

本文最大的技术难点是 $Z_{11}$ 的退化性：其符号在 $\lambda_1=0$ 上为零，导致传统的椭圆估计或强制性论证失效。例如，若允许 $\hat{\phi}$ 集中在 $\lambda_1=0$ 附近，则即使 $\|\phi\|_{L^2} = 1$，也可能使 $\langle Z_{11} \phi, \phi \rangle$ 接近零。

作者的突破性技巧是利用紧支撑性与谱不确定性原理相结合来克服这一退化。具体而言：

  • 假设剖面函数 $Q$ 支撑在有界集 $A$ 上，则其傅里叶变换 $\hat{Q}$ 不可能在频率空间中过度集中在退化区域 $\{\lambda_1 = 0\}$。
  • 通过 Logvinenko-Sereda 型不等式，量化这种“扩散”程度，确保 $\int \frac{\lambda_1^2}{|\lambda|^2} |\hat{Q}|^2 d\lambda$ 有一个正的下界。
  • 将此思想嵌入到算子 $L$ 的变分结构中，证明其强制性，进而获得可逆性。

这一技巧本质上是将空间局域化（紧支撑）转化为频率非集中性，从而激活原本退化的算子。这是本文最核心的创新点。

与之前工作的比较

a) 突破与改进

本文相较于已有工作做出了若干重要推进，主要体现在以下几个方面：

首先，在自相似解的存在性方面，本文是首个在二维非局部模型中建立一般性自相似奇异解存在结果的工作。如文中 Remark 1.1 所强调：“This is the first general result on the existence of self-similar singular solutions for the 2D nonlocal models.” 此前虽然已有针对特定几何（如椭圆域）或特定对称性的显式奇异解构造（见文献 [6]），但缺乏适用于任意有界可测集 $A$ 的统一存在性理论。本文通过泛函分析框架和谱不确定性原理，实现了这一推广，具有更强的普适性。

其次，在爆破准则方面，本文提出的条件比现有结果更为简洁且易于验证。例如，Bianchini 与 Elgindi [1] 的爆破判据需要构造一个辅助函数 $W \in L^2$ 满足多个复杂的积分不等式（包括 $Z_{11} W^2 \omega_0 \geq 0$ 等），这些条件在实际应用中较难检验。而本文的定理 1.2 仅要求初值 $\omega_0$ 属于某个 Hölder 空间 $C^\alpha$、紧支撑且总积分为正——这三个条件直观且自然。正如作者在 Remark 1.2 中指出：“The blow-up criterion in our theorem is also simpler.” 这种简化不仅提升了结果的实用性，也揭示了爆破机制的本质驱动力可能是初值的整体正性而非复杂的局部结构。

第三，在技术路径上，本文首次将谱不确定性原理系统地应用于流体模型的奇点分析。以往研究多依赖于最大值原理、能量估计或精细的渐近分析。而本文通过调和分析中的深奥工具，绕开了对非线性项的逐点控制，转而从频域角度整体把握解的行为。这种方法为未来研究其他具有退化算子的 PDE 提供了新范式。

最后，在假设条件方面，本文的结果在一定程度上放宽了对初值的要求。例如，定理 1.2 中仅需 $\omega_0 \in C^\alpha$（$\alpha>0$）和紧支撑，而不要求特定符号或对称性。相比之下，许多已知的全局存在性结果（如 Constantin 和 Sun [3]）依赖于初值为负，而爆破结果往往需要精心设计的初值。本文表明，只要初值整体为正且局部化，就足以触发爆破，这扩展了已知爆破初值的范围。

综上所述，本文不仅在结果上取得了突破——建立了首个一般的自相似解存在性定理并提出了更简洁的爆破准则——更在方法论上引入了全新的分析工具，为理解非局部流体模型中的奇点形成开辟了新路径。

7. Moisture dynamics with phase changes coupled to heat-conducting, compressible fluids

相变耦合热传导可压缩流体的水分动力学

作者: Felix Brandt, Matthias Hieber, Lin Ma, Tarek Zöchling

类别: math.AP

发布日期: 2025-12-16

arXiv ID: 2512.14931v1

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

详细分析

中文标题: 相变耦合热传导可压缩流体的水分动力学

研究对象和背景

本文研究的是一个将相变过程与热传导、可压缩流体动力学耦合的数学模型，具体描述了大气中水汽、云水和雨水在气流作用下的演化行为。该系统结合了经典的 Navier-Stokes 方程组（用于描述可压缩、非等温流体）与一个微物理水分模型（micro-physics moisture model），后者刻画了水在不同相态之间的转换机制，如蒸发、凝结、自转化和收集等过程。

所研究的控制方程定义在一个三维区域 $\Omega = G \times (0,1)$ 上，其中 $G = (0,1)^2$ 表示水平方向的单位正方形，并假设具有水平周期性边界条件。垂直方向则为有界区间 $(0,1)$，上下边界分别记为 $\Gamma_u = G \times \{1\}$ 和 $\Gamma_l = G \times \{0\}$。时间变量 $t \in (0,\tau)$，$\tau > 0$ 可以是任意大的有限值或无穷。

设：

  • $\rho_d$：干空气密度；
  • $\rho_v, \rho_c, \rho_r$：分别为水蒸气、云水和雨水的密度；
  • 混合比定义为 $q_j := \rho_j / \rho_d$，其中 $j \in \{v,c,r\}$；
  • 总湿空气密度为 $\rho = \rho_d(1 + q_v + q_c + q_r)$；
  • $u = (v,w)$：大气速度场，$v$ 为水平速度，$w$ 为垂直速度；
  • $T$：热力学温度；
  • $p$：压强；
  • $\theta$：位温；
  • $V_r$：雨滴沉降速度，假设为光滑函数。

此外引入若干热力学参数：比热容 $c_{pd}, c_{pv}, c_1$，气体常数 $R_d, R_v$，潜热 $L(T)$，以及饱和混合比 $q_{vs}(p,T)$，它是关于压强和温度的非负、一致有界且Lipschitz连续函数。

系统的完整形式如下（经简化后）：

$$ \begin{cases} \partial_t \rho_d + \mathrm{div}(\rho_d u) = 0, \\ \rho_d Q_m [\partial_t u + (u \cdot \nabla) u] - \mu \Delta u - (\mu + \lambda)\nabla \mathrm{div}\, u + \nabla p = \rho_d Q_m g e_3, \\ \rho_d Q_m \left[\partial_t T + (u \cdot \nabla) T\right] - \Delta T = -\frac{Q_1}{Q_m} T \cdot \frac{1+q_v}{1+q_v+q_c+q_r}(q_{vs} - q_v)^+ q_r \\ \quad + (T + 1)\left[(q_v - q_{vs}) q_c + (q_v - q_{vs})^+\right] - (1 + T)(S_{ac} + S_{cr}), \\ \partial_t q_v + (u \cdot \nabla) q_v - \Delta q_v = T \cdot \frac{1+q_v}{Q_m}(q_{vs} - q_v)^+ q_r - (q_v - q_{vs}) q_c - (q_v - q_{vs})^+, \\ \partial_t q_c + (u \cdot \nabla) q_c - \Delta q_c = (q_v - q_{vs}) q_c + (q_v - q_{vs})^+ - (q_c - 1)^+ - q_c q_r, \\ \partial_t q_r + (u \cdot \nabla) q_r - \Delta q_r = (q_c - 1)^+ + q_c q_r - T \cdot \frac{1+q_v}{Q_m}(q_{vs} - q_v)^+ q_r + V_r \sum_{l=1}^3 \partial_{z_l} q_r Z_{l,3} + \cdots, \end{cases} $$

其中源项 $S_{ev}, S_{cd}, S_{ac}, S_{cr}$ 分别表示蒸发、凝结、自动转化和收集速率，其具体结构包含开关型非线性（switch terms），例如 $(\cdot)^+$ 或 $(\cdot)^-$，这导致这些项在某些点不可微，从而破坏了Fréchet可微性，给数学分析带来显著困难。

边界条件设定如下：

  • 水平方向：所有变量满足周期性边界条件；
  • 垂直边界 $\Gamma_u \cup \Gamma_l$：速度满足无滑移边界条件：
    $$ u = 0 \quad \text{on } \Gamma_u \cup \Gamma_l; $$
  • 温度和混合比满足齐次Neumann边界条件：
    $$ \partial_z T = 0, \quad \partial_z q_j = 0, \quad j \in \{v,c,r\}, \quad \text{on } \Gamma_u \cup \Gamma_l. $$

初始条件为：

$$ \rho_d(0) = \rho_{d,0},\quad u(0) = u_0,\quad T(0) = T_0,\quad q_j(0) = q_{j,0},\quad j \in \{v,c,r\}. $$

值得注意的是，该模型未采用常见的 Boussinesq 近似，而是保留了完整的可压缩、非等温 Navier-Stokes 结构，因此更加贴近真实大气流动，但也带来了更强的非线性和耦合性。

另一个关键特征是相变释放大量潜热，直接反馈到能量方程中，形成强耦合效应。这种“相变触发热源”的机制使得系统的热力学部分高度敏感于水分状态的变化，进一步加剧了解析难度。

从函数空间角度看，作者采用了 $L^p$-$L^q$ 最大正则性框架，而非传统的Hilbert空间方法。具体而言：

  • 初始数据属于空间：
    $$ (\rho_{d,0}, u_0, T_0, q_{j,0}) \in H^{1,q}(\Omega) \times B^{2(1-1/p)}_{q,p}(\Omega)^3 \times B^{2(1-1/p)}_{q,p}(\Omega), $$
    其中 $B^s_{q,p}$ 是Besov空间，$H^{1,q}$ 是Sobolev空间。
  • 解的空间为：
    $$ \rho_d \in H^{1,p}(0,a; H^{1,q}(\Omega)),\quad u \in H^{1,p}(0,a; L^q(\Omega)^3) \cap L^p(0,a; H^{2,q}(\Omega)^3), $$
    对温度和混合比类似。

这一选择避免了对高阶正则性的要求以及复杂的兼容性条件，适用于更一般的初值情形。

特别地，由于连续性方程中含有对流项 $\mathrm{div}(\rho_d u)$，呈现出双曲性，难以直接处理。为此，作者采用拉格朗日坐标变换（Lagrangian coordinates），通过特征线法重写系统，将材料导数变为单纯的时间导数，从而解耦部分非线性影响。

令 $X(t,y)$ 为由速度场 $u$ 生成的特征流：

$$ \begin{cases} \partial_t X(t,y) = u(t,X(t,y)), \\ X(0,y) = y, \end{cases} $$

并设 $Y(t,\cdot) = X(t,\cdot)^{-1}$ 为其逆映射。定义变换后的变量：

$$ \tilde{\rho}_d(t,y) = \rho_d(t,X(t,y)),\quad \tilde{u}(t,y) = u(t,X(t,y)),\quad \text{etc.} $$

则原系统被转化为 Lagrangian 坐标下的新系统（见第7页），其中线性主部保持不变，但出现了一系列依赖于梯度变换 $Z = [\nabla X]^{-1}$ 的扰动项，如 $(\mathcal{L}_1 - \Delta)$、$(\mathcal{L}_2 - \nabla \mathrm{div})$ 等，这些构成了非线性估计中的主要挑战。

综上所述，本研究的核心在于建立这样一个强耦合、含非光滑源项、完全可压缩且考虑热传导的大气流体系统的局部与全局适定性理论，尤其关注如何克服由相变引起的非微分性和双曲-抛物混合结构带来的技术障碍。

主要定理或主要结果

定理 3.1（局部强适定性）

设 $2/p + 3/q < 1$，且初始数据

$$ (\rho_{d,0}, u_0, T_0, q_{j,0}) \in Y_\gamma := H^{1,q}_{\mathrm{per}}(\Omega) \times B^{2(1-1/p)}_{q,p,\mathrm{per}}(\Omega)^3 \times B^{2(1-1/p)}_{q,p,\mathrm{per}}(\Omega),\quad j \in \{v,c,r\}, $$

满足边界条件 $u_0 = 0$，$\partial_z T_0 = \partial_z q_{j,0} = 0$ on $\Gamma_u \cup \Gamma_l$，并且存在常数 $M_1,M_2 > 0$ 使得

$$ M_1 \leq \rho_{d,0}, q_{j,0} \leq M_2 \quad \text{in } \Omega. $$

则存在 $a = a(\rho_{d,0}, u_0, T_0, q_{j,0}) > 0$，使得系统 (2.1) 在 $(0,a)$ 上存在唯一强解 $(\rho_d, u, T, q_v, q_c, q_r)$，满足：

$$ \rho_d \in H^{1,p}(0,a; H^{1,q}(\Omega)),\quad u \in H^{1,p}(0,a; L^q(\Omega)^3) \cap L^p(0,a; H^{2,q}(\Omega)^3), $$

$$ T, q_j \in H^{1,p}(0,a; L^q(\Omega)) \cap L^p(0,a; H^{2,q}(\Omega)),\quad j \in \{v,c,r\}, $$

且存在 $M^*_1, M^*_2 > 0$ 使得 $\rho_d(t,x) \geq M^*_1$, $q_j(t,x) \leq M^*_2$ 对所有 $(t,x) \in (0,a)\times\Omega$ 成立。

定理 3.2（小初值在任意 $[0,\tau]$ 上的强适定性）

设 $\tau > 0$ 任意，$2/p + 3/q < 1$，并假设平衡态 $q_{vs} \equiv 0$，此时 $(\rho_d, 0, T, 0, 1, 0)$ 是系统的一个稳态解。若初值满足

$$ \|(\rho_{d,0} - \rho_d, u_0, T_0 - T, q_{v,0}, q_{c,0} - 1, q_{r,0})\|_{H^{1,q} \times B^{2(1-1/p)}_{q,p}{}^5} \leq \delta $$

对某个小 $\delta > 0$，则系统在 $[0,\tau]$ 上存在唯一强解，且满足一致先验估计：

$$ \|\rho_d - \rho_d\|_{L^\infty(H^{1,q})} + \|\nabla \rho_d\|_{H^{1,p}(L^q)} + \|\partial_t \rho_d\|_{L^p(H^{1,q})} + \|u\|_{L^p(H^{2,q})} + \cdots \leq C\delta, $$

并保证 $\rho_d(t,x) \geq \rho_d/2 > 0$ 对所有 $(t,x) \in [0,\tau]\times\Omega$。

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

a) 核心分析工具

  1. Lagrangian 坐标变换
    针对连续性方程的双曲性，作者采用特征线法引入 Lagrangian 坐标，将原系统的对流导数 $\partial_t + u\cdot\nabla$ 转化为单纯的 $\partial_t$，从而消除一部分非线性。此技巧最早由 Danchin 等人在可压缩流中使用，本文将其推广至多组分、含复杂源项的情形。

  2. 最大 $L^p$-$L^q$ 正则性理论
    作者基于 Da Prato–Grisvard、Weis 等发展的向量值调和分析理论，在 Banach 空间框架下建立最优线性估计。对于线性化系统（如 Lamé 系统、热方程），利用 operator sum method 和 $\mathcal{R}$-有界性 得到解的先验估计：

    $$ \|u\|_{H^{1,p}(L^q) \cap L^p(H^{2,q})} \lesssim \|f_u\|_{L^p(L^q)} + \|u_0\|_{B^{2(1-1/p)}_{q,p}}, $$

    类似估计对密度、温度和混合比也成立。

  3. 不动点原理（Banach 不动点定理）
    将求解非线性系统转化为在适当函数空间中寻找映射 $\Phi: (\tilde{\rho}_d, \tilde{u}, \tilde{T}, \tilde{q}_j) \mapsto (\hat{\rho}_d, \hat{u}, \hat{T}, \hat{q}_j)$ 的不动点。通过精细估计非线性项 $G = (G_d, G_u, G_T, G_j)$ 的大小，证明 $\Phi$ 是自映射且压缩，从而得到局部解的存在唯一性。

  4. 能量估计与 Besov 空间嵌入
    利用 Sobolev 嵌入 $B^{2(1-1/p)}_{q,p} \hookrightarrow H^{1,q}$（当 $2/p + 3/q < 1$ 时成立）来控制初值正则性；同时结合非线性复合函数估计（如 Moser-type 不等式）处理 $Q_m, Q_{th}, \sigma$ 等依赖于未知量的系数。

b) 技术难点与应对策略

  1. 非光滑源项问题
    源项如 $(q_v - q_{vs})^+$、$(q_c - 1)^+$ 缺乏 Fréchet 可微性，无法直接应用隐函数定理。解决办法是将其视为外力项（forcing term），在非线性估计中作为已知项处理。尽管不光滑，但在 $L^p$ 框架下仍可通过 Lipschitz 控制（如 $(a^+ - b^+) \leq |a-b|$）进行估计。

  2. 变换后项的粗糙性
    在 Lagrangian 坐标下，拉普拉斯算子变为 $\mathcal{L}_1 f = \sum \partial_{y_k y_l} f Z_{k,j} Z_{l,j} + \cdots$，其与标准 Laplacian 的差 $(\mathcal{L}_1 - \Delta)f$ 含有高阶导数乘以低阶项，需用乘积估计（product estimates in Besov spaces）加以控制。

  3. 全局小解的稳定性分析
    对于任意 $\tau > 0$ 的全局适定性，必须在接近平衡态时进行精细谱分析。作者通过对线性化算子的渐近稳定性分析，结合非线性项的二次小量性质（quadratic smallness），得出解不会偏离太远，从而延拓到任意时间。

  4. 边界条件的变换处理
    Neumann 条件在坐标变换下变为斜导数边界条件（如 $\partial_z T = B_T[\nabla T, Z]$），需借助 trace 理论和边界算子的有界性来封闭估计。

与之前工作的比较

a) 突破与改进

  1. 首次实现任意时间区间的全局适定性
    文献中已有工作（如 Doppler et al. [13]）仅在 $H^2$ 框架下证明了局部适定性。本文首次在 $L^p$-$L^q$ 设置下获得对任意 $\tau > 0$ 的强解存在性（对小初值），这是目前最接近“全局解”的结果。

  2. 放宽初始数据正则性要求
    以往多数结果依赖高阶 Sobolev 空间（如 $H^2$）或 Hilbert 空间设置，而本文允许初值在更低正则性的 Besov 空间 $B^{2(1-1/p)}_{q,p}$ 中，只要满足 $2/p + 3/q < 1$ 即可。这意味着可以处理更广泛的初始扰动类型。

  3. 无需 Boussinesq 或弱解框架
    多数大气模型采用 Boussinesq 近似简化密度变化，或在弱解框架下研究（如 Lions, Feireisl）。本文坚持强解路径，保留全可压缩结构，更具物理真实性。

  4. 处理非微分源项的新方法
    相变模型中普遍存在的开关函数长期被视为障碍。本文通过将非光滑项视为强制项，并结合现代最大正则性工具，绕开了传统需要 Fréchet 可微性的限制，为同类问题提供了新范式。

综上，本文在模型完整性、数学严格性、适用范围三方面均实现了显著推进，标志着对“相变-流体”强耦合系统的分析进入新阶段。

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

1. Effective Equations for a Compressible Liquid-Vapor Flow Model with Highly Oscillating Initial Density

具有高度振荡初始密度的可压缩液-气流模型的有效方程

作者: Christian Rohde, Florian Wendt

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15535v1

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

摘要翻译

中文标题: 具有高度振荡初始密度的可压缩液-气流模型的有效方程

摘要翻译:
本文推导并论证了一种针对喷雾尺度下可压缩黏性液-气流动的新有效模型，即适用于存在大量相界面情形的模型。在细观尺度上，我们以Navier-Stokes-Korteweg系统的抛物型松弛模型为出发点。考虑一组初始数据序列，其中初始密度序列呈高度振荡状态，以模拟初始时刻存在大量相边界的物理情景。随后，我们考察对应于该初始数据序列的一族有限能量弱解。预期有效方程将在初始相变数量趋于无穷的极限情况下出现，因此我们将密度解释为Young测度，并证明了解序列收敛至该有效模型。该有效模型由两部分组成：一部分描述流体水动力学量的确定性演化，另一部分为描述极限Young测度的输运方程，用以刻画混合动力学行为。通过将Young测度用其对应的累积分布函数进行表征，我们将关于Young测度的动能方程重写为关于累积分布函数的动能方程，从而使所得方程组可适用于标准的逼近方法处理。

2. A plethora of fully localised solitary waves for the full-dispersion Kadomtsev-Petviashvili equation

全色散Kadomtsev-Petviashvili方程的大量完全局域化孤立波

作者: Mats Ehrnström, Mark D. Groves

类别: math.AP, math-ph

发布日期: 2025-12-17

arXiv ID: 2512.15487v1

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

摘要翻译

中文标题: 全色散Kadomtsev-Petviashvili方程的大量完全局域化孤立波

摘要翻译: KP-I方程是Bond数 $β>1/3$（即强表面张力）情形下重力-毛细表面波的一个弱非线性模型方程。最近研究表明，该方程具有一族非退化的、对称的“完全局域化”孤立波，亦称为“lump”孤立波，其在所有空间方向上均趋于零。全色散KP-I方程（FDKP-I）通过在建模过程中保留水波问题中的精确色散关系而得到。本文证明了FDKP-I方程同样具有一族对称的完全局域化孤立波；这些解可通过将该方程视为KP-I方程的扰动，并应用适当形式的隐函数定理而获得。

3. Linear systems, determinants and solutions of the Kadomtsev-Petviashvili equation

线性系统、行列式与Kadomtsev-Petviashvili方程的解

作者: Gordon Blower, Simon J. Malham

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15245v1

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

摘要翻译

中文标题: 线性系统、行列式与Kadomtsev-Petviashvili方程的解

摘要翻译: 设 $(-A,B,C)$ 是一个定义在连续时间 $t>0$ 上的线性系统，其输入与输出空间为 ${\mathbb C}$，状态空间为 $H$。散射（或脉冲响应）函数 $φ_{(x)}(t)=Ce^{-(t+2x)A}B$ 确定了一个Hankel积分算子 $Γ_{φ_{(x)}}$；若 $Γ_{φ_{(x)}}$ 为迹类算子，则Fredholm行列式 $τ(x)=\det (I+Γ_{φ_{(x)}})$ 给出了系统 $(-A,B,C)$ 的tau函数。本文研究了包含算子 $R_x = \int_x^\infty e^{-tA}BCe^{-tA}\,dt$ 在 $H$ 上所生成的代数结构，并由此得到了Kadomtsev-Petviashvili偏微分方程的解。Pöppe的半可加算子被识别为积分核上某种平移作用的轨道，而Pöppe的括号运算则可用Fedosov积来表示。文章进一步指出，Fredholm行列式 $\det (I+R_x)$ 为Kadomtsev-Petviashvili方程解的数值计算提供了一种有效方法。

4. Modeling of a micropolar thin film flow with rapidly varying thickness and non-standard boundary conditions

具有快速变化厚度和非标准边界条件的微极性薄膜流动建模

作者: María Anguiano, Francisco J. Suárez-Grau

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15202v1

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

摘要翻译

中文标题: 具有快速变化厚度和非标准边界条件的微极性薄膜流动建模

摘要翻译: 本文研究微极性流体在薄区域中流动的渐近行为，假设上边界为零Dirichlet边界条件且具有快速振荡的形貌，下边界为平坦表面并施加非标准边界条件。在“雷诺粗糙度 regime”下，即区域厚度远小于粗糙度波长（表现为极轻微的粗糙度）的情形，我们严格推导出一个广义的雷诺压力方程，清晰地揭示了由粗糙度诱导的效应。此外，我们还给出了平均速度和微转动的表达式。

不相关论文（基本信息）

1. What does it mean for a 3D star-shaped scatterer to be small in the time domain?

一个三维星形散射体在时域中“小”意味着什么？

作者: Maryna Kachanovska, Adrian Savchuk

类别: math.AP, math-ph

发布日期: 2025-12-17

arXiv ID: 2512.15629v1

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

摘要: In the frequency domain wave scattering problems, obstacles can be effectively replaced by point scatterers as soon as the wavelength of the incident wave exceeds significantly their diameter. The situation is less clear in the time domain, where recent works suggest the presence of an additional temporal scale that quantifies the smallness of the obstacle. In this paper we argue that this is not necessarily the case, and that it is possible to construct asymptotic models with an error that does not deteriorate in time, at least in the case of a sound-soft scattering problem by a star-shaped obstacle in 3D.

2. A cheap way to closed operator sums

一种构造闭算子和的简便方法

作者: Bernhard H. Haak, Peer Christian Kunstmann

类别: math.FA, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15594v1

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

摘要: Let $A$ and $B$ be sectorial operators in a Banach space $X$ of angles $ω_A$ and $ω_B$, respectively, where $ω_A+ω_B<π$. We present a simple and common approach to results on closedness of the operator sum $A+B$, based on Littlewood-Paley type norms and tools from several interpolation theories. This allows us to give short proofs for the well-known results due to Da~Prato-Grisvard and Kalton-Weis. We prove a new result in $\ell^q$-interpolation spaces and illustrate it with a maximal regularity result for abstract parabolic equations. Our approach also yields a new proof for the Dore-Venni result.

3. A modified Bakry-Émery $Γ_2$ criterion inequality and the monotonicity of the Tsallis entropy

一个修正的Bakry-Émery $Γ_2$ 判别不等式与Tsallis熵的单调性

作者: Xiaohan Cai, Xiaodong Wang

类别: math.DG, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15525v1

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

摘要: The Bakry-Émery $Γ_2$ criterion inequality provides a method for establishing the logarithmic Sobolev inequality. We prove a one-parameter family of weighted Bakry-Émery $Γ_2$ criterion inequalities which in the limit case yields the improved constant due to Ji \cite{Ji24}. Furthermore, we establish a modified weighted $Γ_2$ criterion inequality which could be interpreted as a monotonicity of the Tsallis entropy under the heat flow and yields a family of sharp Sobolev inequalities.

4. A Nonlinear elliptic PDE with curve singularity on the boundary

具有边界曲线奇性的非线性椭圆型偏微分方程

作者: Mamadou Ciss, Abdourahmane Diatta, El Hadji Abdoulaye Thiam

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15475v1

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

摘要: Let $Ω$ be a bounded domain of $\mathbb{R}^{N+1}$ ($N \geq 3$) with smooth boundary $\partial Ω$ and $Σ$ be a closed submanifold contained on $\partial Ω$ and containing $0$. We are interesting in the existence of positive $H^1(Ω)$-solution of the following Hardy-Sobolev trace type equation \begin{equation*} \begin{cases} -Δu+u=0 \qquad & \textrm{ in $Ω$}\\ \displaystyle\frac{\partial u}{\partial ν}= ρ_Σ^{-s} u^{q_s-1} \qquad & \textrm{ on $\partial Ω$}, \end{cases} \end{equation*} where $ν$ is the unit outer normal of $\partial Ω$, $ρ_Σ: \partial Ω\to \mathbb{R}$ is the distance function in $\partial Ω$ to the curve $Σ$:

$$ ρ_Σ(x):= \inf_{y \in Σ} d_{\tilde{g}}(x, y) $$

and for $0\leq s <1$, $q_s:=\frac{2(N-s)}{N-1}$ is the critical Hardy-Sobolev exponent. The existence of solution may depend on the local geometry of the boundary $\partial Ω$ and $Σ$ at $0$ or in the shapes of the domain $Ω$ and its boundary $\partial Ω$.

5. Capillary $L_p$-Christoffel-Minkowski problem

毛细 $L_p$-Christoffel-Minkowski 问题

作者: Yingxiang Hu, Mohammad N. Ivaki

类别: math.AP, math.DG

发布日期: 2025-12-17

arXiv ID: 2512.15464v1

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

摘要: We solve the capillary $L_p$-Christoffel–Minkowski problem in the half-space for $1<p<k+1$ in the class of even hypersurfaces. A crucial ingredient is a non-collapsing estimate that yields lower bounds for both the height and the capillary support function. Our result extends the capillary Christoffel–Minkowski existence result of \cite{HIS25}.

6. On a relaxed Cahn-Hilliard tumour growth model with single-well potential and degenerate mobility

关于具单势阱和退化迁移率的松弛Cahn-Hilliard肿瘤生长模型

作者: Cecilia Cavaterra, Matteo Fornoni, Maurizio Grasselli, Benoît Perthame

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15426v1

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

摘要: We consider a phase-field system modelling solid tumour growth. This system consists of a Cahn-Hilliard equation coupled with a nutrient equation. The former is characterised by a degenerate mobility and a singular potential. Both equations are subject to suitable reaction terms which model proliferation and nutrient consumption. Chemotactic effects are also taken into account. Adding an elliptic regularisation, depending on a relaxation parameter $δ>0$, in the equation for the chemical potential, we prove the existence of a weak solution to an initial and boundary value problem for the relaxed system. Then, we let $δ$ go to zero, and we recover the existence of a weak solution to the original system.

7. Trace theory for parabolic boundary value problems with rough boundary conditions

粗边界条件下抛物型边值问题的迹理论

作者: Robert Denk, Floris B. Roodenburg

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15382v1

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

摘要: We characterise the trace spaces arising from intersections of weighted, vector-valued Sobolev spaces, where the weights are powers of the distance to the boundary. These weighted function spaces are particularly suitable for treating boundary value problems where derivatives of the solution blow up at the boundary. As an application of our trace theory, we prove well-posedness for the heat equation with rough inhomogeneous boundary data in Sobolev spaces of higher regularity in domains of fixed regularity $C^{1,κ}$, with $κ\in [0,1)$.

8. Spatial Approximation for Evolutionary Equations

演化方程的空间逼近

作者: Andreas Buchinger, Christian Seifert, Sascha Trostorff, Marcus Waurick

类别: math.FA, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15354v1

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

摘要: We consider evolutionary equations as introduced by R.\ Picard in 2009 and develop a general theory for approximation which can be seen as a theoretical foundation for numerical analysis for evolutionary equations. To demonstrate the approximation result, we apply it to a spatial discretisation of the heat equation using spectral methods.

9. Fully non-linear elliptic equations on noncompact complex manifolds

非紧复流形上的完全非线性椭圆方程

作者: Hanzhang Yin

类别: math.DG, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15350v1

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

摘要: In this paper, we obtain the a priori estimates and the existence results for solutions of a general class of fully non-linear equations on noncompact Kähler and Hermitian manifolds. As geometric applications, we can construct Kähler metrics with some prescribed volume forms on strictly pseudoconvex domains; and we can construct complete Chern-Einstein metrics on Hermitian manifolds with bounded geometry.

10. Weak curvature conditions on metric graphs

度量图上的弱曲率条件

作者: Juliane Krautz

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15329v1

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

摘要: Starting from pointwise gradient estimates for the heat semigroup, we study three characterizations of weak lower curvature bounds on metric graphs. More precisely, we prove the equivalence between a weak notion of the Bakry-Émery curvature condition, a weak Evolutionary Variational Inequality and a weak form of geodesic convexity. The proof is based on a careful regularization of absolutely continuous curves together with an explicit representation of the Cheeger energy. We conclude with a brief discussion on possible applications to the Schrödinger bridge problem on metric graphs.

11. A Nonhomogeneous Boundary-Value Problem For The Nonlinear KdV Equation on Star Graphs

星形图上非线性KdV方程的一类非齐次边值问题

作者: Roberto de A. Capistrano Filho, Hugo Parada, Jandeilson Santos da Silva

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15307v1

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

摘要: This paper investigates a boundary-value problem for the Korteweg-de Vries (KdV) equation on a star-graph structure. We develop a unified framework introducing the notion of $s$-compatibility, which generalizes classical compatibility conditions to star-shaped and more complex graph configurations, inspired by the works of Bona, Sun, and Zhang [14]. By combining analytical techniques with a fixed-point argument, we establish sharp global well-posedness for both the linear and nonlinear problems at the $H^s$ level. In this setting, our results extend the classical analysis for a single KdV equation [14] to star-shaped graphs composed of $N$ equations. These results provide the first comprehensive well-posedness theory for KdV equations with coupled boundary conditions on graphs. Although control issues are not treated in this article, the analytic results obtained here address several open problems, which will be addressed in a forthcoming

12. Weak Error on the densities for the Euler scheme of stable additive SDEs with Besov drift

具有Besov漂移的稳定加性SDEs在Euler格式下的密度弱误差

作者: Mathis Fitoussi, Elena Issoglio, Stéphane Menozzi

类别: math.AP, math.PR

发布日期: 2025-12-17

arXiv ID: 2512.15299v1

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

摘要: We are interested in the Euler-Maruyama dicretization of the formal SDE, $dX_t=b(t,X_t)dt+dZ_t$, where $Z$ is a symmetric isotropic d dimensional stable process of index $α\in (1,2)$, and $b$ is distributional. It belongs to a mix Lebesgue-Besov space. The associated parameters satisfy some constraints which guarantee weak-well posedness. Defining an appropriate Euler scheme, we obtain a convergence rate for the weak error on the densities. The rate depends on the parameters.

13. Carleson-type removability for $p$-parabolic equations

关于 $p$-抛物方程的Carleson型可去性

作者: Michał Borowski, Theo Elenius, Leah Schätzler, David Stolnicki

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15277v1

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

摘要: We characterize removable sets for Hölder continuous solutions to degenerate parabolic equations of $p$-growth. A sufficient and necessary condition for a set to be removable is given in terms of an intrinsic parabolic Hausdorff measure, which depends on the considered Hölder exponent. We present a new method to prove the sufficient condition, which relies only on fundamental properties of the obstacle problem and supersolutions, and applies to a general class of operators. For the necessity of the condition, we establish the Hölder continuity of solutions with measure data, provided the measure satisfies a suitable decay property. The techniques developed in this article provide a new point of view even in the case $p=2$.

14. Very weak solutions to degenerate parabolic double-phase systems

退化抛物型双相系统的极弱解

作者: Wontae Kim, Lauri Särkiö

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15255v1

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

摘要: We prove a local self-improving property for the gradient of very weak solutions to degenerate parabolic double-phase systems. The result is based on a reverse Hölder inequality with constants that are independent of the solution. Delicate methods are required to avoid a self-referential argument. In particular, we develop a new phase analysis method.

15. Strichartz estimates in Wiener amalgam spaces for Schrödinger equations with at most quadratic potentials

关于具有至多二次势的薛定谔方程在维纳 amalgam 空间中的Strichartz估计

作者: Shun Takizawa

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15218v1

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

摘要: For Schrödinger equations with potentials which grow at most quadratically at spatial infinity, we prove Strichartz estimates in Wiener amalgam spaces. These estimates provide a stronger recovery of local-in-space regularity than the classical Strichartz estimates in Lebesgue spaces. Our result is a generalization of the results on Strichartz estimates in Wiener amalgam spaces by Cordero and Nicola, which are stated for the potentials $V(x) = 0,|x|^2/2, -|x|^2/2$.

16. A multiscale framework integrating within-host infection kinetics with airborne transmission dynamics

整合宿主内感染动力学与空气传播动态的多尺度框架

作者: Andrew Omame, Sarafa Iyaniwura

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

发布日期: 2025-12-17

arXiv ID: 2512.15209v1

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

摘要: Coupling within-host infection dynamics with population-level transmission remains a major challenge in infectious disease modeling, especially for airborne pathogens with potential to spread indoor. The frequent emergence of such diseases highlight the need for integrated frameworks that capture both individual-level infection kinetics and between-host transmission. While analytical models for each scale exist, tractable approaches that link them remain limited. In this study, we present a novel multiscale mathematical framework that integrates within-host infection kinetics with airborne transmission dynamics. The model represents each host as a patch and couples a system of ordinary differential equations (ODEs) describing in-host infection kinetics with a diffusion-based partial differential equation (PDE) for airborne pathogen movement in enclosed spaces. These scales are linked through boundary conditions on each patch boundary, representing viral shedding and inhalation. Using matched asymptotic analysis in the regime of intermediate diffusivity, we derived a nonlinear ODE model from the coupled ODE-PDE system that retains spatial heterogeneity through Neumann Green’s functions. We established the existence, uniqueness, and boundedness of solutions to the reduced model and analyzed within-host infection kinetics as functions of the airborne pathogen diffusion rate and host spatial configuration. In the well-mixed limit, the model recovers the classical target cell limited viral dynamics framework. Overall, the proposed multiscale modeling approach enables the simultaneous study of transient within-host infection dynamics and population-level disease spread, providing a tractable yet biologically grounded framework for investigating airborne disease transmission in indoor environments.

17. Well-Posedness of Pseudo-Parabolic Gradient Systems with State-Dependent Dynamics

具有状态依赖动力学的拟抛物梯度系统的适定性

作者: Harbir Antil, Daiki Mizuno, Ken Shirakawa, Naotaka Ukai

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15164v1

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

摘要: This paper develops a general mathematical framework for pseudo-parabolic gradient systems with state-dependent dynamics. The state dependence is induced by variable coefficient fields in the governing energy functional. Such coefficients arise naturally in scientific and technological models, including state-dependent mobilities in KWC-type grain boundary motion and variable orientation-adaptation operators in anisotropic image denoising. We establish two main results: the existence of energy-dissipating solutions, and the uniqueness and continuous dependence on initial data. The proposed framework yields a general well-posedness theory for a broad class of nonlinear evolutionary systems driven by state-dependent operators. As illustrative applications, we present an anisotropic image-denoising model and a new pseudo-parabolic KWC-type model for anisotropic grain boundary motion, and prove that both fit naturally within the abstract structure of $(\mathrm{S})_ν$.

18. A New Fast Finite Difference Scheme for Tempered Time Fractional Advection-Dispersion Equation with a Weak Singularity at Initial Time

一种求解初始时刻具弱奇性的截断时间分数阶对流-扩散方程的新快速有限差分格式

作者: Liangcai Huang, Shujuan Lü

类别: math.NA, math.AP

发布日期: 2025-12-17

arXiv ID: 2512.15141v1

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

摘要: In this paper, we propose a new second-order fast finite difference scheme in time for solving the Tempered Time Fractional Advection-Dispersion Equation. Under the assumption that the solution is nonsmooth at the initial time, we investigate the uniqueness, stability, and convergence of the scheme. Furthermore, we prove that the scheme achieves second-order convergence in both time and space. Finally, corresponding numerical examples are provided.

19. Exact solution structures on some nonlocal overdetermined problems

一些非局部超定问题的精确解结构

作者: Kazuki Sato, Futoshi Takahashi

类别: math.AP

发布日期: 2025-12-17

arXiv ID: 2512.14987v1

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

摘要: In this paper, we study the solution structures of Serrin-type overdetermined problems with Kirchhoff-type nonlocal terms. We prove that the exact number of solutions is the same as those of some transcendental equations defined by the nonlocal terms. We also obtain the explicit form of solutions by using the unique solutions of the overdetermined problems without the nonlocal terms.

20. Unbounded Branches of Non-Radial Solutions to Semilinear Elliptic Systems on a Unit Ball in $\mathbb R^3$ and Their Patterns

单位球上半线性椭圆系统非径向解的无界分支及其模式研究

作者: Casey Crane, Ziad Ghanem

类别: math.AP

发布日期: 2025-12-16

arXiv ID: 2512.14957v1

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

摘要: We investigate symmetry-breaking phenomena in semilinear elliptic systems on the unit ball in $\mathbb{R}^3$, focusing on the emergence of non-radial solution branches with prescribed spatial and internal symmetries. Extending previous scalar results, we develop a framework for systems equivariant under $G := O(3) \times Γ\times \mathbb{Z}_2$, where $Γ$ is a finite group encoding coupling symmetries. Using the $G$-equivariant Leray–Schauder degree and Burnside ring techniques, we derive computable criteria for the existence of unbounded branches of non-radial solutions and classify their isotropy types. Our approach accommodates non-simple eigenvalue multiplicities and provides explicit bifurcation conditions in terms of spectral resonance between coupling eigenvalues and spherical Laplacian modes. Applications to coupled spherical oscillators illustrate how Platonic symmetries and internal permutations interact to produce complex patterns. These results establish a general method for detecting and characterizing symmetry-breaking in high-dimensional elliptic systems.

21. Asymptotic formulas for $L^2$ bifurcation curves of nonlocal logistic equation of population dynamics

种群动力学中非局部逻辑斯蒂方程的$L^2$分支曲线的渐近公式

作者: Tetsutaro Shibata

类别: math.AP

发布日期: 2025-12-16

arXiv ID: 2512.14955v1

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

摘要: The one-dimensional nonlocal Kirchhoff type bifurcation problems which are derived from logistic equation of population dynamics are studied. We obtain the precise asymptotic shapes of $L^2$ bifurcation curves $λ= λ(α)$ as $α\to \infty$, where $α:= \Vert u_λ\Vert_2$.

22. Nonlinear oscillators at resonance with periodic forcing

处于共振状态的非线性振子在周期激励下的响应

作者: Philip Korman, Yi Li

类别: math.DS, math-ph, math.AP, math.CA

发布日期: 2025-12-16

arXiv ID: 2512.14940v1

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

摘要: In this note we unify the results of A.C. Lazer and P.O. Frederickson [3], A.C. Lazer [6], A.C. Lazer and D.E. Leach [7], J.M. Alonso and R. Ortega [1], and P. Korman and Y. Li [4] on periodic oscillations and unbounded solutions of nonlinear equations with linear part at resonance and periodic forcing. We give conditions for the existence and non-existence of periodic solutions, and obtain a rather detailed description of the dynamics for nonlinear oscillations at resonance, in case periodic solutions do not exist.

23. Relations between principal eigenvalue and torsional rigidity with Robin boundary conditions

关于Robin边界条件下主特征值与扭转刚度之间关系的研究

作者: Giuseppe Buttazzo, Simone Cito, Francesco Solombrino

类别: math.AP

发布日期: 2025-12-16

arXiv ID: 2512.14927v1

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

摘要: We consider the torsional rigidity and the principal eigenvalue related to the Laplace operator with Dirichlet and Robin boundary conditions. The goal is to find upper and lower bounds to products of suitable powers of the quantities above in the class of Lipschitz domains. The threshold exponent for the Robin case is explicitly recovered and shown to be strictly smaller than in the Dirichlet one.

24. Chemotaxis models with signal-dependent sensitivity and a logistic-type source, I: Boundedness and global existence

具有信号依赖性敏感性和逻辑斯蒂型源的趋化模型，I：有界性与整体存在性

作者: Le Chen, Ian Ruau, Wenxian Shen

类别: math.AP

发布日期: 2025-12-16

arXiv ID: 2512.14858v1

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

摘要: We study, in Part I of this series, boundedness and global existence of positive classical solutions to a parabolic-elliptic chemotaxis system with signal-dependent sensitivity and a logistic-type source on a bounded smooth domain $Ω\subset\mathbb{R}^N$: \begin{equation*} \begin{cases} \displaystyle u_t=Δu-χ_0\nabla\cdot\left(\frac{u^m}{(1+v)^β}\nabla v\right)+au-bu^{1+α}, & x\inΩ, \cr \displaystyle 0=Δv-μv+νu^γ, & x\inΩ, \cr \displaystyle \frac{\partial u}{\partial n}=\frac{\partial v}{\partial n}=0, & x\in\partialΩ. \end{cases} \end{equation*} Here, $u$ denotes the population density and $v$ the chemical concentration. The parameters $α,γ,m,μ,ν$ are positive, $χ_0$ is real, and $a,b,β$ are nonnegative. We analyze boundedness from three viewpoints: negative chemotaxis ($χ_0<0$), the strength of the nonlinear cross diffusion rate $\frac{u^m}{(1+v)^β}$, and the strength of the logistic-type damping $u(a-bu^α)$. Under explicit conditions reflecting these mechanisms, all positive classical solutions remain bounded. Moreover, when $m\ge 1$, boundedness implies global existence. Although the decay of $χ(v) = \dfrac{χ_0}{(1+v)^β}$ for large $v$ has a damping effect, it also introduces new analytical difficulties; our techniques yield, for example, global existence for $m=1$ provided that \begin{equation*} β>\max\left{1,\frac12+\frac{χ_0}{4}\max{2,γN}\right}. \end{equation*} Several known results for special cases are recovered. Part II is devoted to the asymptotic behavior of globally defined solutions, including uniform persistence as well as stability and bifurcation of positive constant equilibria.