生成时间: 2026年04月22日 03:42:21

重点关注论文数量: 5

了解领域论文数量: 9

不相关论文数量: 12

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

1. 广义SQG方程在无粘极限下的哈密顿紧性与耗散

Hamiltonian compactness and dissipation for the generalized SQG equation in the inviscid limit

作者: Luigi De Rosa, Utku Kemal Yuzbasioglu

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19732v1

摘要翻译: 我们考虑具有任意阶拉普拉斯分数幂耗散的耗散型广义表面拟地转（SQG）方程。在无粘极限下，我们证明：解在使非线性项良定义的最低范数下的强紧性可阻止哈密顿量的反常耗散。事实上，仅特定频率上的动力学起作用。该论证具有高度鲁棒性，适用于任意临界性情形，且不依赖于外部强迫项是否存在（即使该强迫项可能非紧）。这揭示了近期关于纳维–斯托克斯方程及临界耗散型SQG方程所得若干结果背后更一般的机制。由于非唯一性问题，在我们更广泛的框架中，必须采用满足适当高阶先验估计的Leray解；我们证明了此类解的存在性，该结果本身可能具有独立兴趣。最后，我们证明：对任一具有临界可积性的初值，强紧性均被保证，从而推出无粘问题存在全局保守的弱解——尽管这些解处于Onsager超临界范畴。这是迄今已知能保证全局存在的最大初值类，且与Delort在端点情形所考虑的初值类相吻合。

Comment: 29 pages

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

1. 研究对象和背景

本文研究广义表面准地转方程（gSQG）在无粘极限下的哈密顿紧性与耗散问题。考虑参数 $\alpha \in (0, 1]$ 和 $\gamma \in (0, \infty)$ ，方程为：

\begin{cases} \partial_t \theta^\nu + u^\nu \cdot \nabla \theta^\nu + \nu (-\Delta)^\gamma \theta^\nu = f^\nu \\ u^\nu = \mathcal{R}^\perp_\alpha \theta^\nu \\ \theta^\nu(\cdot, 0) = \theta^\nu_0 \end{cases} \quad (\text{gSQG}_\nu)

其中 $\mathcal{R}^\perp_\alpha := \nabla^\perp (-\Delta)^{-\alpha}$ 为正交广义Riesz变换， $\nu \in (0, 1)$ 为粘性参数， $\theta^\nu_0$ 为零平均初始数据， $f^\nu$ 为外部力，定义在 $\mathbb{T}^2 \times [0, \infty)$ 上。通过与 $\theta^\nu$ 在 $\dot{H}^{-\alpha}(\mathbb{T}^2)$ 内积，得哈密顿量平衡：

\frac{1}{2} \frac{d}{dt} \|\theta^\nu(t)\|_{\dot{H}^{-\alpha}}^2 + \nu \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 = \langle f^\nu(t), \theta^\nu(t) \rangle_{\dot{H}^{-\alpha}}

$\dot{H}^{-\alpha}(\mathbb{T}^2)$ 是使非线性项良定义的最小正则性空间。核心问题是：当 $\nu \to 0^+$ 时，是否发生反常耗散，即 $\lim_{\nu \to 0^+} \nu \int_0^T \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 dt = 0$ ？根据 $\alpha$ 与 $\gamma$ 的关系，方程分为超临界（ $\alpha < 1-\gamma$ ）、临界（ $\alpha = 1-\gamma$ ）和次临界（ $\alpha > 1-\gamma$ ）情形，其中 $\alpha = 1$ 对应Euler涡度形式， $\alpha = 1/2$ 对应标准SQG。已有工作表明，在 $\alpha = \gamma = 1$ （Navier-Stokes）[54] 和 $\alpha = \gamma = 1/2$ （临界SQG）[26] 时，强紧性与反常耗散不相容，揭示了二维与三维湍流的本质差异。本文旨在建立适用于任意 $\alpha \in (0,1]$ 和 $\gamma \in (0,\infty)$ 的统一框架，解决非唯一性问题，填补现有理论在一般参数下强紧性与耗散关系的空白。

2. 主要定理或主要结果

定理1.1（紧性与耗散）：设 $\alpha \in (0, 1]$ 和 $\gamma \in (0, \infty)$ 任意。设 $\{\theta^\nu_0\}_\nu$ 为零平均初始数据序列， $\{f^\nu\}_\nu$ 为外部力序列（几乎处处时间有零空间平均）。假设：(i) $\{\theta^\nu_0\}_\nu$ 在 $\dot{H}^{-\alpha}(\mathbb{T}^2)$ 中强紧；(ii) 若 $\alpha \leq \gamma$ ，则 $\{f^\nu\}_\nu$ 在 $L^{\frac{\alpha+\gamma}{\gamma}}([0, T]; \dot{H}^{-\alpha}(\mathbb{T}^2))$ 中有界；若 $\alpha > \gamma$ ，则 $\{f^\nu\}_\nu$ 在 $L^\infty([0, T]; \dot{H}^{-\gamma}(\mathbb{T}^2))$ 中有界。对由命题2.8构造的 (gSQG $_\nu$ ) 的弱解序列 $\{\theta^\nu\}_\nu$ ，有

\{\theta^\nu\}_\nu \subset L^2([0, T]; \dot{H}^{-\alpha}(\mathbb{T}^2)) \text{ 强紧} \quad \Longrightarrow \quad \lim_{\nu \to 0^+} \nu \int_0^T \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 dt = 0

定理1.3（全局存在性）：设 $\alpha \in (0, 1)$ 任意，令 $p_\alpha := \frac{2}{1+\alpha}$ 。设 $\theta_0 \in L^{p_\alpha}(\mathbb{T}^2)$ 为零平均初始数据， $f \in L^1_{\text{loc}}([0, \infty); L^{p_\alpha}(\mathbb{T}^2)) \cap L^{1+\alpha}_{\text{loc}}([0, \infty); \dot{H}^{-\alpha}(\mathbb{T}^2))$ 且几乎处处时间有零空间平均。则存在全局弱解

\theta \in C^0([0, \infty); \dot{H}^{-\alpha}(\mathbb{T}^2)) \cap C^0_w([0, \infty); L^{p_\alpha}(\mathbb{T}^2))

满足

\|\theta(t)\|_{\dot{H}^{-\alpha}}^2 = \|\theta_0\|_{\dot{H}^{-\alpha}}^2 + 2 \int_0^t \langle f(\tau), \theta(\tau) \rangle_{\dot{H}^{-\alpha}} d\tau

和

\|\theta(t)\|_{L^{p_\alpha}} \leq \|\theta_0\|_{L^{p_\alpha}} + \int_0^t \|f(\tau)\|_{L^{p_\alpha}} d\tau

对所有 $t \geq 0$ 成立。推论2.13（频率依赖紧性）：若存在 $N(\nu) \to \infty$ 使得 $\{\theta^\nu_{\leq N(\nu)}\}_\nu$ 在 $L^2([0,T]; \dot{H}^{-\alpha}(\mathbb{T}^2))$ 中强紧，且 $\nu N(\nu)^{2\gamma} \to 0$ ，则 $\lim_{\nu \to 0^+} \nu \int_0^T \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 dt = 0$ 。特别地，当 $N(\nu) \sim \nu^{-1/(2\gamma)}$ 时条件满足。

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

证明包含五个关键步骤：1. 高阶估计：对任意 $\delta > 0$ ，证明 $\nu^{\frac{\alpha+\gamma}{\gamma}} \int_\delta^T \|\theta^\nu(t)\|_{\dot{H}^\gamma}^2 dt \lesssim 1$ 。由能量平衡 $\frac{d}{dt} \|\theta^\nu(t)\|_{L^2}^2 = -2\nu \|\theta^\nu(t)\|_{\dot{H}^\gamma}^2$ 出发，结合插值 $\|\theta^\nu(t)\|_{L^2} \lesssim \|\theta^\nu(t)\|_{\dot{H}^{-\alpha}}^{\frac{\gamma}{\alpha+\gamma}} \|\theta^\nu(t)\|_{\dot{H}^\gamma}^{\frac{\alpha}{\alpha+\gamma}}$ 及 $\|\theta^\nu(t)\|_{\dot{H}^{-\alpha}} \leq \|\theta^\nu_0\|_{\dot{H}^{-\alpha}}$ ，推得 $\frac{d}{dt} \|\theta^\nu(t)\|_{L^2}^2 \lesssim -\nu \|\theta^\nu(t)\|_{L^2}^{\frac{2(\alpha+\gamma)}{\alpha}}$ ，积分得 $\|\theta^\nu(t)\|_{L^2}^2 \lesssim (\nu t)^{-\alpha/\gamma}$ ，进而 $\nu \int_\delta^T \|\theta^\nu(t)\|_{\dot{H}^\gamma}^2 dt \leq \|\theta^\nu(\delta)\|_{L^2}^2 \lesssim (\nu \delta)^{-\alpha/\gamma}$ 。当 $\alpha > \gamma$ 时，需借助Grönwall论证。2. 正时间耗散分裂：对任意 $\delta > 0$ 和 $$N > 0$$ ，证明 $\nu \int_\delta^T \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 dt \lesssim \nu N^{2\gamma} + \Phi^{\frac{\alpha}{\alpha+\gamma}}(N)$ ，其中 $\Phi(N) := \sup_{\nu > 0} \int_\delta^T \|\theta^\nu_{>N}(t)\|_{\dot{H}^{-\alpha}}^2 dt$ 。通过频率分解 $\theta^\nu = \theta^\nu_{\leq N} + \theta^\nu_{>N}$ ，低频部分由 $\|\theta^\nu_{\leq N}(t)\|_{\dot{H}^{\gamma-\alpha}}^2 \leq N^{2\gamma} \|\theta^\nu(t)\|_{\dot{H}^{-\alpha}}^2 \lesssim N^{2\gamma}$ 控制，高频部分用插值 $\|\theta^\nu_{>N}(t)\|_{\dot{H}^{\gamma-\alpha}} \lesssim \|\theta^\nu_{>N}(t)\|_{\dot{H}^{-\alpha}}^{\frac{\alpha}{\alpha+\gamma}} \|\theta^\nu(t)\|_{\dot{H}^\gamma}^{\frac{\gamma}{\alpha+\gamma}}$ ，再结合Hölder不等式与高阶估计完成。3. 无瞬时耗散：证明 $\{\theta^\nu_0\}_\nu \subset \dot{H}^{-\alpha}(\mathbb{T}^2)$ 强紧 $\Rightarrow$ $\lim_{\delta \to 0^+} \sup_{\nu > 0} \nu \int_0^\delta \|\theta^\nu(t)\|_{\dot{H}^{\gamma-\alpha}}^2 dt = 0$ ，通过弱解公式中选取适当测试函数实现，不依赖非线性结构。4. 特殊Leray解的构造：为克服非唯一性[2,62]，构造满足高阶估计的Leray解（命题2.8）。对任意 $\theta^\nu_0 \in \dot{H}^{-\alpha}(\mathbb{T}^2)$ ，构造正则化近似序列，使其满足 $\nu^{\frac{\alpha+\gamma}{\gamma}} \int_\delta^T \|\theta^\nu(t)\|_{\dot{H}^\gamma}^2 dt \lesssim 1$ ，再通过Aubin-Lions引理获得 $L^2([0,\infty); \dot{H}^{-\alpha}(\mathbb{T}^2))$ 紧性。5. 临界初始数据的紧性证明：对定理1.3，需证明 $\{\theta^\nu\}_\nu$ 在 $L^2_{\text{loc}}([0,\infty); \dot{H}^{-\alpha}(\mathbb{T}^2))$ 中强紧。由于 $L^{p_\alpha}(\mathbb{T}^2) \subset \dot{H}^{-\alpha}(\mathbb{T}^2)$ 仅为连续嵌入（ $p_\alpha = 2/(1+\alpha)$ 为临界指数），作者通过引理3.1将解分裂为 $\theta^\nu = \theta^\nu_{\text{small}} + \theta^\nu_{\text{reg}}$ ，其中 $\|\theta^\nu_{\text{small}}\|_{L^{p_\alpha}} \leq \varepsilon$ ， $\theta^\nu_{\text{reg}}$ 具更高正则性，利用 $\varepsilon$ 任意性消除高频集中，实现强紧性。

4. 与之前工作的比较

普遍机制的揭示：本文证明 $\dot{H}^{-\alpha}(\mathbb{T}^2)$ 是唯一满足以下三性质的空间：(1) 使非线性项在分布意义下良定义；(2) 在无粘情况下形式守恒；(3) 其强紧性阻止反常耗散。该机制统一了[15,26,45,54]等工作的核心思想。先前结果仅限于特定参数：[54]处理 $\alpha = \gamma = 1$ （Navier-Stokes），[26]处理 $\alpha = \gamma = 1/2$ （临界SQG）。本文推广至任意 $\alpha \in (0,1]$ 和 $\gamma \in (0,\infty)$ ，并允许非紧致外力，显著扩展适用范围。2. 与凸积分构造解的对比：Onsager猜想对无粘活性标量成立[23,40,42,72]，凸积分构造的’wild’解正则性为 $L^3([0,T]; B^{2-4\alpha/3}_{3,\infty}(\mathbb{T}^2))$ ，紧嵌入于 $L^2([0,T]; \dot{H}^{-\alpha}(\mathbb{T}^2))$ ，但此类解无法通过物理合理的弱解在 $\nu \to 0^+$ 极限中实现，因其正则性不关于粘性参数一致，突显二维与三维湍流本质区别。3. 临界性阈值的区分：(gSQG $_\nu$ ) 的临界性阈值为 $\alpha = 1-\gamma$ ，而高阶估计所需阈值为 $\alpha = \gamma$ ，二者不同。当 $\alpha = \gamma$ 时，高阶估计自然导出；一般情形下需额外操作，且 $\alpha > \gamma$ 时因耗散不足需更强外力假设（ $L^\infty$ 而非 $L^{(\alpha+\gamma)/\gamma}$ ）。4. 非唯一性问题的处理：与[24,45,54]中Navier-Stokes不同，(gSQG $_\nu$ ) 的Leray解可能不唯一[2,62]。本文仅考虑满足高阶估计(2.2)的解（命题2.8），否则定理1.1可能失效（备注2.1），此为关键创新。5. 初始数据范围的扩展：定理1.3提供目前已知最大初始数据类，使无粘gSQG全局弱解存在。当 $\alpha = 1/2$ 时恢复[26, Theorem 1.1]；当 $\alpha \to 1^-$ 时， $p_\alpha \to 1^+$ ，与Euler方程Delort结果[27-30]在端点匹配。 $\alpha = 1$ 端点因涉及带符号奇异测度初始数据被排除，因方法不同。6. 外力项的鲁棒性：序列 $\{f^\nu\}_\nu$ 仅需有界，无需紧致，体现方法强鲁棒性，表明[9,10,16,17,46]中三维受迫Navier-Stokes的反常耗散构造不能扩展至二维gSQG族。7. 与磁流体动力学（MHD）的类比：Taylor猜想[35,36]指出理想MHD中磁螺旋度在零粘性/电阻率极限下守恒，超越Onsager阈值[37,48]，与本文一致。备注1.2指出，二维湍流的’逆级联’（能量向大尺度传递）可能是定理1.1的根本原因，而三维湍流的’正向级联’允许强紧性与反常耗散共存。8. 频率依赖紧性的推广：推论2.13证明仅需特定频率（ $N(\nu)$ ）上强紧且 $\nu N(\nu)^{2\gamma} \to 0$ 即可排除耗散异常，条件弱于定理1.1。当 $N(\nu) \sim \nu^{-1/(2\gamma)}$ 时满足，揭示耗散由特定频率动力学主导，为理解湍流能量传递提供新视角。

2. 二维重力-毛细水波中带涡度的小尺度生成

Small scale creation in 2D gravity-capillary water waves with vorticity

作者: Yuanpeng Tu

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19380v1

摘要翻译: 本文考虑定义在有限深度、具有自由表面与固定底部的无界区域上的二维不可压Euler方程；流体运动受重力与表面张力共同作用。我们构造了具有平坦自由表面和小初速度的初始数据，使得相应解在其存在时间区间内，涡度梯度的 $L^\infty$ 范数至少以双重指数速率增长。本工作将Zlatos关于固定边界的结论推广至自由表面情形，并将Hu–Luo–Yao关于有界区域的结果推广至无界区域情形。

Comment: 21 pages, 2 figures

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

1. 研究对象和背景

本文研究二维不可压缩Euler方程在无界域中带自由表面和固定底部（有限深度）的流体运动，受重力和表面张力共同作用。流体区域 $$D_t$$ 位于固定底部 $Γ_b = \mathbb{R} \times \{x_2 = 0\}$ 和上自由表面 $$Γ_t$$ 之间，初始流体域为 $D_0 = \mathbb{R} \times (0, 2)$ 。流体速度 $$u(t, x)$$ 和压力 $$p(t, x)$$ 满足控制方程：

\begin{cases} \partial_t u + u \cdot \nabla u + \nabla p + g e_2 = 0, & \text{in } D_t, \\ \nabla \cdot u = 0, & \text{in } D_t, \end{cases}

其中 $$x = (x_1, x_2)$$ ， $\nabla = (\partial_1, \partial_2)$ ， $$e_2 = (0, 1)$$ 。自由表面运动满足：

V = u \cdot N \quad \text{on } \u0393_t,

其中 $$N$$ 是 $\u0393_t$ 的外法向量；表面张力效应由边界条件：

p = \sigma H \quad \text{on } \u0393_t,

描述， $\sigma > 0$ 为表面张力常数， $$H$$ 为 $\u0393_t$ 的平均曲率；固定底部满足滑移边界条件：

u \cdot e_2 = 0 \quad \text{on } \u0393_b.

关键未知量为速度场 $$u(t, x)$$ 、压力场 $$p(t, x)$$ 和自由边界 $\u0393_t$ ，涡度定义为 $\omega(t, x) := \partial_1 u_2 - \partial_2 u_1$ 。该问题的重要性在于：一方面，水波系统的长时间行为是流体力学中的基本问题；另一方面，涡度对水波稳定性的影响尚未被充分理解。已有理论表明，无旋水波（ $\omega \equiv 0$ ）在小初值下具有长时间适定性（如 $O(\varepsilon^{-2})$ 尺度），但非常数涡度的长时间存在性结果极少。Kiselev–Šverák [15] 在有界域中构造了涡度梯度双指数增长的例子，Zlatoš [33] 将其推广到半平面（无界域但无自由表面），Hu–Luo–Yao [10] 则推广到有界自由表面域。然而，无界自由表面域（水平方向为 $\mathbb{R}$ ）中的小尺度生成问题仍未解决，因为：(1) 自由表面演化导致域变形难以控制；(2) 无界域使得传统基于周长的能量估计失效；(3) 自由边界设置缺乏精确Biot-Savart定律。本文旨在填补这一空白，构造初始自由表面平坦、速度任意小的光滑初值，使得涡度梯度 $\|\nabla \omega(t)\|_{L^\infty}$ 在解的生命跨度内至少具有双指数增长速率。

2. 主要定理或主要结果

在 $\sigma > 0$ ， $$g > 0$$ 情形下，主要结果为定理1.1：存在光滑速度场 $v_0 \in C^\infty(D_0)$ 和通用常数 $\varepsilon_0, c_1, c_2 > 0$ ，使得对任意 $\varepsilon \in (0, \varepsilon_0)$ ，以初值 $u_0 = \varepsilon v_0$ 和初始域 $D_0 = \mathbb{R} \times (0, 2)$ 的解满足：

\|\nabla \omega(t, \cdot)\|_{L^\infty(D_t)} \geq \varepsilon \exp\big(c_1 \exp(c_2 \varepsilon t)\big)

对所有 $t \in [0, T)$ 成立，其中 $$T$$ 为解的生命跨度， $\omega(t, x) = \partial_1 u_2 - \partial_2 u_1$ 。该估计的尺度关系为 $t \sim O(\varepsilon^{-1} \log \log \varepsilon^{-1})$ 时， $\|\nabla \omega\|_{L^\infty}$ 增长至 $$O(1)$$ 量级。推论1.2进一步指出，此时 $\|u\|_{W^{2,\infty}}$ 也增长至 $$O(1)$$ ，与无旋水波的 $O(\varepsilon^{-2})$ 长时间存在性形成尖锐对比。证明中关键依赖于对称性假设：初始速度满足 $$-u_1^0(-x_1, x_2) = u_1^0(x_1, x_2)$$ 和 $$u_2^0(-x_1, x_2) = u_2^0(x_1, x_2)$$ ，这保证了涡度 $\omega(t, x)$ 在 $$x_1$$ 上保持奇性，并通过近似Biot-Savart定律导出速度场的点态估计：

|u(x) - u(y)| \leq C \|\omega\|_{L^\infty} |x - y| \left(1 + \log^+ \frac{\text{diam}(D_t)}{|x - y|} \right),

结合能量守恒 $$E(t) = E(0)$$ 和自由表面与底部的分离估计 $\u0393_t \subset \mathbb{R} \times (3/2, 5/2)$ ，最终通过追踪涡度极值点附近的演化得到双指数增长下界。

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

核心障碍在于：(1) 非零涡度导致流体域 $$D_t$$ 随时间显著变形，自由表面 $\u0393_t$ 可能接触底部；(2) 无界域（水平方向 $\mathbb{R}$ ）使传统周长估计失效；(3) 自由边界设置缺乏精确Biot-Savart定律。作者的创新点如下：

Steiner对称化处理无界域：针对无界域问题，作者将注意力限制在 $$|x_1| < R$$ 的截断域内，对上区域 $D^+(R) = D_t(R) \backslash D_0$ 和下区域 $D^-(R) = (D_0 \backslash D_t(R)) \cap \{|x_1| < R\}$ 分别应用两次Steiner对称化。第一次 $S_{e_1}$ 将截断域内最高/最低点重定位至 $$x_1 = 0$$ ，第二次 $S_{e_2}$ 将不规则区域转换为边界为函数图形的规则区域。关键在于通过重定位保证相对高度在对称化下不变，从而控制自由边界变形。具体地，定义对称化后区域：

\tilde{D}^+(R) := \{(x_1, x_2) \in \mathbb{R}^2 : (x_1, 2 + |x_2|) \in D^+(R)\},

D^*_+(R) := S_{e_2}(S_{e_1}(\tilde{D}^+(R))) \cap (\mathbb{R} \times \mathbb{R}^+),

并证明其边界为函数图形，从而将一般情形转化为图形情形进行能量估计。

重归一化能量泛函：为克服无界域问题，作者构造了重归一化能量：

E(t) = K(t) + P(t), \quad K(t) = \frac{1}{2} \int_{D_t} |u|^2 dx,

P(t) = \sigma \lim_{R \to \infty} \left( \int_{\u0393_t \cap \{|x_1| < R\}} dS_t - 2R \right) + \lim_{R \to \infty} \int_{D_t \cap \{|x_1| < R\}} g(x_2 - 1) dx.

通过精细分析边界项 $$R_K(R, t)$$ 和 $R_\sigma(R, t)$ （见(2.3)-(2.5)），证明当 $R \to \infty$ 时：

\frac{d}{dt}E_R(t, R) = R_K(R, t) + \sigma R_\sigma(R, t) \to 0,

从而得到能量守恒 $$E(t) = E(0)$$ 。该能量泛函的关键在于表面张力部分控制自由边界超出初始长度的额外部分，重力部分控制流体区域的垂直位移。

自由表面与底部的分离：利用能量守恒，作者证明当初始动能 $K(0) < C_1 \varepsilon^2$ 足够小时，自由表面始终满足 $\u0393_t \subset \mathbb{R} \times (3/2, 5/2)$ 。证明分两步：(i) 当 $\u0393_t$ 为图形时，通过分解势能 $P(t) = P_1(\phi) + P_2(\phi)$ ，其中 $P_1(\phi) = \int \sigma (\sqrt{1 + \phi'^2} - 1) dx_1$ 控制表面张力， $P_2(\phi) = \frac{g}{2} \int \phi^2 dx_1$ 控制重力，导出 $P(t) \geq c > 0$ 与 $E(0) \leq C_1 \varepsilon^2$ 矛盾；(ii) 一般情形通过Steiner对称化转化为图形情形，利用对称化不增加势能的性质完成证明。

近似Biot-Savart定律的应用：尽管自由边界设置无精确Biot-Savart定律，但作者采用[10]中的近似形式，获得速度场的点态估计：

u(x) = \frac{1}{2\pi} \int_{D_t} \frac{(x - y)^\perp}{|x - y|^2} \omega(y) dy + E(x),

其中误差项 $$E(x)$$ 满足 $\|E\|_{L^\infty} \leq C \|\omega\|_{L^\infty}$ 。结合对称性（涡度 $\omega$ 为 $$x_1$$ 奇函数）和自由表面与底部的分离，导出关键估计：

|u_1(x)| \leq C \|\omega\|_{L^\infty} x_1 \left(1 + \log^+ \frac{1}{x_1} \right), \quad |u_2(x)| \leq C \|\omega\|_{L^\infty} x_1,

从而得到涡度梯度演化方程 $\partial_t \nabla \omega + u \cdot \nabla \nabla \omega = -\nabla u \cdot \nabla \omega$ 沿特征线的双指数增长下界。

4. 与之前工作的比较

与Zlatoš [33]的比较：该文证明了半平面（无界域）中2D Euler方程的涡度梯度双指数增长，但未考虑自由表面。本文将其推广到自由表面设置，主要推进在于处理了自由边界演化带来的域变形问题。关键差异在于：(i) Zlatoš [33] 的固定边界允许直接使用半平面Green函数，而本文需通过Steiner对称化控制自由表面；(ii) 本文证明了自由表面与底部的分离 $\u0393_t \subset \mathbb{R} \times (3/2, 5/2)$ ，这是Zlatoš [33] 不需要的；(iii) 适用范围从固定边界扩展到自由边界系统，代价是需要更复杂的能量估计和对称化技巧。时间尺度上两者均为双指数增长，但本文的常数 $$c_2$$ 依赖于表面张力 $\sigma$ 和重力 $$g$$ 。

与Hu–Luo–Yao [10]的比较：该文将Kiselev–Šverák [15] 推广到周期性自由表面域（水平方向 $\mathbb{T}$ ），但域仍有界。本文的主要推进在于处理了无界域情形（水平方向 $\mathbb{R}$ ）。关键差异在于：(i) [10] 利用周期性将势能 $$P(t)$$ 与自由边界周长成正比，而本文因无界性需重定义能量泛函并引入极限过程；(ii) [10] 直接假设自由表面为图形，本文需通过Steiner对称化将一般情形转化为图形情形；(iii) 适用范围从有界周期域扩展到无界域，代价是技术上更复杂（需两次Steiner对称化）且无法获得全局周长控制。值得注意的是，当 $R \to \infty$ 时，本文的重归一化能量退化为[10]的周期情形。

与无旋水波长时间适定性结果的比较：文献[3,6,8-9,11-14,24-25,28-29,34]证明了小初值无旋水波的 $O(\varepsilon^{-2})$ 甚至更长的生命跨度。本文的推论1.2表明，有旋情形下 $\|u\|_{W^{2,\infty}}$ 在 $O(\varepsilon^{-1} \log \log \varepsilon^{-1})$ 内增长至 $$O(1)$$ ，这比无旋情形短得多。关键差异在于：(i) 无旋水波的涡度 $\omega \equiv 0$ 消除了涡度梯度增长机制；(ii) 本文构造的初值具有特殊对称性（ $$u_1$$ 奇、 $$u_2$$ 偶），导致涡旋对向原点加速运动。这解释了为何非常数涡度的长时间存在性结果极少——涡度可能通过小尺度生成机制破坏正则性。

与Su [23]的比较：该文证明了2D无限深度水波中特定初始涡对（向下移动对称涡对）的 $O(\varepsilon^{-2})$ 生命跨度。本文的推进在于：(i) 考虑有限深度情形，更具物理现实性；(ii) 构造了导致更快正则性损失的初值（ $O(\varepsilon^{-1} \log \log \varepsilon^{-1})$ vs $O(\varepsilon^{-2})$ ）；(iii) 不依赖特定涡对结构，而是利用自由表面与底部的相互作用。局限在于Su [23] 的结果表明某些特殊初值可能有较长生命跨度，而本文仅证明存在导致快速增长的初值，两者互补说明涡度对水波稳定性的影响具有初值依赖性。

3. 球面上欧拉方程的涡量梯度增长

Growth of vorticity gradient for the Euler equation on the sphere

作者: Daomin Cao, Junhong Fan, Guolin Qin

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19358v1

摘要翻译: 我们证明：球面上欧拉方程的解，其涡量梯度在时间上至多呈双重指数增长；并通过构造具有奇对称性的显式解，在半球上实现了该双重指数增长，从而表明该上界是尖锐的。我们还将该结果推广至旋转球面情形。据我们所知，这是关于非平凡几何紧致流形上理想流体涡量梯度增长的首个结果。

Comment: 无

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

1. 研究对象和背景

本文研究单位球面 $S^2 := \{x = (x_1, x_2, x_3) \in \mathbb{R}^3 : x_1^2 + x_2^2 + x_3^2 = 1\}$ 上不可压缩理想流体的运动，其动力学由球面上的不可压缩欧拉方程描述：

\begin{cases} \partial_t v + \nabla_v v = -\nabla P, & \text{on } S^2 \times \mathbb{R}^+, \\ \operatorname{div} v = 0, & \text{on } S^2 \times \mathbb{R}^+, \end{cases}

其中 $v(x,t) \in T_x S^2$ 是切于球面的速度场， $$P$$ 为标量压力， $\nabla_v$ 表示沿 $$v$$ 的协变导数， $\nabla$ 和 $\operatorname{div}$ 分别为球面上的标准度量下的梯度和散度算子。等价地，用涡度 $\omega$ 表示的方程为：

\begin{cases} \partial_t \omega + u \cdot \nabla \omega = 0, & \text{on } S^2 \times \mathbb{R}^+, \\ u(x,t) = \nabla^\perp G\omega(x,t), & \text{on } S^2 \times \mathbb{R}^+, \end{cases}

其中 $G\omega(x,t) = \int_{S^2} G(x,y)\omega(y,t) \, d\sigma(y)$ ， $G(x,y) = \frac{1}{2\pi} \ln |x-y|$ 为格林函数， $\omega$ 需满足高斯约束 $\int_{S^2} \omega(x,t) \, d\sigma(x) = 0$ ，且 $$u$$ 需满足不可穿透条件 $\tilde{u}_\theta(\varphi, \pm\pi/2, t) = 0$ 。

在二维欧拉方程研究中，最大涡度幅度 $\|\omega\|_{L^\infty}$ 是守恒量，但其梯度 $\|\nabla \omega\|_{L^\infty}$ 可能增长，反映了小尺度结构的持续生成。量化涡度梯度的增长率是流体动力学中的经典问题。Wolibner、H"older 和 Yudovich 证明了梯度增长不超过双指数速率：

\|\nabla \omega(\cdot, t)\|_{L^\infty} \leq (1 + \|\nabla \omega_0\|_{L^\infty})^{C \exp(Ct)},

其中常数 $$C$$ 仅依赖于 $\|\omega_0\|_{L^\infty}$ 。这一上界在具有光滑边界的区域上成立，但在边界不光滑的区域（如仅有两个角点）可能发生有限时间爆破（Kiselev 和 Zlato\v{s}）。然而，该上界是否紧致（即是否存在解确实以双指数速率增长）长期是核心开放问题。

Kiselev 和 \v{S}ver’ak [19] 在圆盘上构造了具有奇对称性的光滑初始涡度，证明了双指数上界可实现，其机制依赖于边界和 Bahouri-Chemin 涡的稳定性。Denisov [7,8] 在无边界但紧致的环面 $$T^2$$ 上证明了无限时间超线性增长，并后来表明双指数增长可在任意长的有限时间间隔内实现。Zlato\v{s} [37] 通过考虑稳态 $\omega^*(x) = \operatorname{sgn}(x_1)\operatorname{sgn}(x_2)$ 的扰动，证明了 $C^{1,\alpha}$ 类初始数据（ $\alpha < 1$ ）的指数增长。最近，Zlato\v{s} [38] 在半平面上利用奇对称性和边界效应，首次在无界域上实现了双指数增长。

然而，球面 $$S^2$$ 这一自然几何结构在涡度梯度增长研究中几乎未被探索。作为行星表面流体流动的基本模型，球面欧拉方程对地球物理和气象学具有重要意义。从数学角度看，球面是具有非平凡曲率的紧致无边流形，其几何特性与平面、环面或半平面有本质区别。本文旨在填补这一空白，研究球面上涡度梯度的增长机制，探究平面域中已知的双指数增长机制能否在球面上运作，并确定最优增长率。这一研究不仅推进了二维欧拉动力学在曲面流形上的理论，也对理解行星尺度大气和海洋流动中的小尺度生成机制具有重要价值。

2. 主要定理或主要结果

在 $\tilde{\omega}(\cdot, 0) \in W^{1,\infty}([\!-\pi, \pi) \times (0, \pi/2])$ 或 $\tilde{\omega}(\cdot, 0) \in W^{1,\infty}([\!-\pi, \pi) \times [\!-\pi/2, \pi/2])$ 且 $\|\tilde{\omega}(\cdot, 0)\|_{L^\infty} = 1$ 的假设下，定理 1.1 证明了球面上欧拉方程解的涡度梯度满足：

\lim_{t \to \infty} \frac{\ln \ln \|\nabla \tilde{\omega}(\cdot, t)\|_{L^\infty([\!-\pi,\pi)\times(0,\pi/2])}}{t} \leq \frac{2}{\pi}.

该结果表明球面上涡度梯度增长至多为双指数速率，且增长率上界为 $2/\pi$ 。此上界通过精确估计球面格林函数和速度场得到，无需任何对称性假设。

定理 1.2 进一步证明了该上界是紧致的：存在 $\tilde{\omega}(\cdot, 0) \in C^\infty([\!-\pi, \pi) \times (0, \pi/2])$ 且 $\|\tilde{\omega}(\cdot, 0)\|_{L^\infty} = 1$ ，使得半球 $$S^2_+$$ 上的解满足：

\lim_{t \to \infty} \frac{\ln \ln \|\nabla \tilde{\omega}(\cdot, t)\|_{L^\infty([\!-\pi,\pi)\times(0,\pi/2])}}{t} = \frac{2}{\pi}.

这里半球上的欧拉方程为：

\begin{cases} \partial_t \omega + u \cdot \nabla \omega = 0, & \text{on } S^2_+ \times \mathbb{R}^+, \\ u(x,t) = \nabla^\perp G_+\omega(x,t), & \text{on } S^2_+ \times \mathbb{R}^+, \end{cases}

其中 $G_+\omega(x,t) = \int_{S^2} G_+(x,y)\omega(y,t) \, d\sigma(y)$ ， $G_+(x,y) = \frac{1}{2\pi} \ln \frac{|x-y|}{|x-\bar{y}|}$ ， $\bar{y} = (y_1,y_2,-y_3)$ 。证明中通过施加赤道和子午线的奇-奇对称性（即 $\omega(x) = -\omega(\tilde{x}) = -\omega(\bar{x})$ ，其中 $\tilde{x} = (x_1,-x_2,x_3)$ ），将动力学简化到四分之一球面 $$S^2_q$$ ，并构造了实现最优增长率的显式解。

注 1.3 指出，若初始涡度满足 $\|\tilde{\omega}(\cdot, 0)\|_{L^\infty} = b > 0$ ，则双指数增长率 $2/\pi$ 被替换为 $2b/\pi$ 。注 1.4 指出，构造的解在半球边界（赤道）处存在奇异性，整个球面上光滑解能否实现双指数增长仍是开放问题。此外，第 4 节将结果扩展到了旋转球面情形，这对地球物理应用具有重要意义。

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

主要技术障碍在于球面几何的复杂性：球面作为具有正曲率的紧致无边流形，其 Biot-Savart 定律和速度场表达式比平面域更为复杂，且缺乏真正的边界来压缩涡度梯度。作者的核心创新在于将平面域中的奇对称性策略适应到球面几何，通过赤道和子午线的奇-奇对称性在四分之一球面 $$S^2_q$$ 上创建有效的’角’压缩机制。

关键工具是引理 2.1，它给出了奇-奇对称假设下球坐标中速度分量的精确渐近表达式。定义 $$x_0 = (1,0,0)$$ 和区域 $Q(k_1,k_2) = \{(\cos \theta' \cos \varphi', \cos \theta' \sin \varphi', \sin \theta') : \varphi' \in (k_1,\pi), \theta' \in (k_2,\pi/2)\}$ ，则对 $(\varphi,\theta) \in [0,1]^2 \setminus \{(0,0)\}$ 有：

\tilde{u}_\varphi(\varphi, \theta) = \frac{x_2}{\sqrt{1-x_3^2}} \left( -\frac{4}{\pi} \int_{Q(3\varphi/2, 3\theta/2)} \frac{y_2 y_3}{|y-x_0|^4} \omega(y) \, d\sigma(y) + B_\varphi(\varphi, \theta) \right),

\tilde{u}_\theta(\varphi, \theta) = \frac{x_3}{\sqrt{1-x_3^2}} \left( \frac{4x_1}{\pi} \int_{Q(3\varphi/2, 3\theta/2)} \frac{y_2 y_3}{|y-x_0|^4} \omega(y) \, d\sigma(y) + B_\theta(\varphi, \theta) \right),

其中余项满足：

|B_\varphi(\varphi, \theta)| \leq C_1 \|\tilde{\omega}\|_{L^\infty} \left( 1 + \min\left\{ \ln\left(1 + \frac{\theta}{\varphi}\right), \frac{\|\nabla \tilde{\omega}\|_{L^\infty([0,3\theta/2]^2)}}{\|\tilde{\omega}\|_{L^\infty}} \theta \right\} \right),

|B_\theta(\varphi, \theta)| \leq C_1 \|\tilde{\omega}\|_{L^\infty} \left( 1 + \min\left\{ \ln\left(1 + \frac{\varphi}{\theta}\right), \frac{\|\nabla \tilde{\omega}\|_{L^\infty([0,3\varphi/2]^2)}}{\|\tilde{\omega}\|_{L^\infty}} \varphi \right\} \right).

该引理的证明涉及将积分区域划分为主导区域 $Q(3\varphi/2, 3\theta/2)$ 和余项区域 $$Q_2$$ , $$Q_3$$ ，并利用球面距离的等价关系 $|x-y| \sim \sqrt{(\varphi-\varphi')^2 + (\theta-\theta')^2}$ 进行精细估计。在主导区域，通过展开 $|x-y|^{-2}|x-\tilde{y}|^{-2} - |y-x_0|^{-4} = O(|x-x_0||y-x_0|^{-5})$ 得到主部积分；在余项区域，通过变量替换和对称性估计余项大小，例如对 $$Q_3$$ 区域有：

\int_{Q_3} \frac{2y_2(x_3-y_3)}{|x-y|^2|x-\tilde{y}|^2} \omega(y) \, d\sigma(y) \lesssim \|\tilde{\omega}\|_{L^\infty} \left( \ln\left(1 + \frac{\theta}{\varphi}\right) + 1 \right).

当 $\varphi < \theta$ 时，进一步将 $$Q_3$$ 分为 $$Q_3'$$ 和 $$Q_4$$ ，并在 $$Q_4$$ 上利用 $\omega(y) = \omega_1(y_1,y_2) + \omega_2(y)$ 的分解（其中 $\omega_2(y)$ 由梯度控制）得到更精确估计。

定理 1.1 的证明基于引理 2.1 和粒子轨迹分析：通过追踪两点间距离的演化，结合速度差的估计 $|u(x,t) - u(y,t)| \lesssim \|\omega\|_{L^\infty} \ln \frac{1}{|x-y|} + C$ ，推导出 $\|\nabla \omega(\cdot,t)\|_{L^\infty}$ 与初始梯度的双指数关系。

定理 1.2 的证明核心是构造一族收缩域 $\alpha(t)\Omega_{\varepsilon(t)}$ 。在奇-奇对称下，四分之一球面 $$S^2_q$$ 中的对称轴 $\{\varphi = 0\}$ 附近，速度场的渐近行为为 $\tilde{u}_\varphi \sim -\frac{4}{\pi} \theta \ln \frac{1}{\theta}$ ， $\tilde{u}_\theta \sim \frac{4}{\pi} \varphi \ln \frac{1}{\theta}$ 。选择 $\alpha(t) = \exp(-c_1 t)$ 和 $\varepsilon(t)$ 满足 $\frac{d\varepsilon}{dt} = -\frac{4}{\pi} \varepsilon \ln \frac{1}{\varepsilon} + O(\varepsilon)$ ，则解 $\varepsilon(t) \sim \exp(-c_2 e^{(4/\pi)t})$ 表明涡度梯度以最优速率 $2/\pi$ 呈双指数增长。通过验证 $\omega \equiv 1$ 的区域边界不进入收缩域，确保了构造的有效性。

4. 与之前工作的比较

本文与先前工作的比较可从多个维度展开：

1. 几何结构的拓展：

Kiselev 和 \v{S}ver’ak [19] 在圆盘（有界域）上证明了双指数增长的可实现性，依赖于边界和 Bahouri-Chemin 涡的稳定性。球面 $$S^2$$ 作为紧致无边流形，缺乏真正的边界，作者通过奇-奇对称性在赤道创建’有效边界’，将动力学简化到半球 $$S^2_+$$ ，并在四分之一球面 $$S^2_q$$ 上构造压缩机制。
Denisov [7,8] 在环面 $$T^2$$ （平坦紧致流形）上证明了超线性增长和有限时间双指数增长，但环面无曲率，而球面具有正曲率，导致 Biot-Savart 核心和速度估计更为复杂。例如，球面上的弦距离 $$|x-y|$$ 与角度距离的关系不同于平面欧氏距离。
Zlato\v{s} [38] 在半平面上利用奇对称性实现了双指数增长，但半平面有直线边界，而球面的’边界’（赤道）是弯曲的，引入了额外的几何因子（如 $\sqrt{1-x_3^2}$ ）。本文将半平面技术适应到弯曲边界，证明了 $$S^2$$ 上的增长率 $2/\pi$ 与半平面的 $1/\pi$ 不同，体现了几何对增长率的影响。

2. 增长率的精确确定：

本文首次在具有非平凡曲率的紧致流形上确定了涡度梯度的最优双指数增长率 $2/\pi$ ，而先前工作仅在平面域给出上界或有限时间结果。注 1.3 表明该增长率与初始涡度幅度 $$b$$ 成正比（ $2b/\pi$ ），推广了 Zlato\v{s} [38] 中半平面的 $b/\pi$ 结果。
与环面 $$T^2$$ 上 Denisov [8] 的结果相比，球面增长率 $2/\pi$ 是全局的（对所有时间成立），而环面上的双指数增长仅在任意长的有限时间间隔内实现，且未确定最优常数。
与 Kiselev 和 \v{S}ver’ak [19] 在圆盘上的结果相比，球面增长率 $2/\pi$ 小于圆盘的 $1/\pi$ （归一化后），反映了球面几何对涡度压缩的抑制效应。

3. 技术方法的推进：

引理 2.1 中的速度估计比平面域更精细：余项 $B_\varphi$ , $B_\theta$ 同时包含对数项和梯度控制项，而 [19] 仅需对数估计。这是由于球面曲率导致的额外耦合效应，例如在 $$Q_4$$ 区域需分解 $\omega(y) = \omega_1 + \omega_2$ 并利用 $|\omega_2| \leq \|\nabla \tilde{\omega}\|_{L^\infty} |\theta-\theta'|$ 。
本文无需 [19] 中 Bahouri-Chemin 涡的稳定性分析，而是直接通过收缩域 $\alpha(t)\Omega_{\varepsilon(t)}$ 的构造证明最优性，其中 $\varepsilon(t)$ 满足微分不等式 $\frac{d\varepsilon}{dt} \leq -\frac{4}{\pi} \varepsilon \ln \frac{1}{\varepsilon} + C\varepsilon$ ，解为 $\varepsilon(t) \leq \exp(-c e^{(4/\pi)t})$ 。
与 [38] 相比，球面速度场的渐近行为 $\tilde{u}_\varphi \sim -\frac{4}{\pi} \theta \ln \frac{1}{\theta}$ 包含几何因子 $x_2/\sqrt{1-x_3^2}$ ，导致增长率常数从半平面的 $1/\pi$ 变为 $2/\pi$ 。

4. 适用范围与局限：

本文结果适用于球面和旋转球面（第 4 节），对地球物理流体动力学（如大气和海洋流动）有直接应用价值，而 [19,38] 仅处理静态域。
构造的实现双指数增长的解在半球边界（赤道）处有奇异性（注 1.4），而整个球面上光滑解能否实现双指数增长仍是开放问题，这与 [19] 中圆盘上的光滑解形成对比。
与 [17] 相比，本文需要更强的奇-奇对称性（关于赤道和子午线），而 [17] 仅需小扰动和鞍点，但 [17] 未确定最优增长率。

5. 代价与推进：

本文放宽了 [19] 对边界的依赖，但引入了对称性假设；推进了 [8] 的有限时间结果到全局时间，并确定了最优常数。
与 [37] 相比，本文处理的是光滑解而非 $C^{1,\alpha}$ 解，但需对称性；而 [37] 无需对称性但仅得指数增长（非双指数）。
代价是构造的解在边界有奇异性，且结果局限于球面对称几何，而 [20] 在非光滑边界域证明了有限时间爆破，表明几何正则性对增长率上界至关重要。

综上，本文首次建立了具有非平凡几何结构的紧致流形上涡度梯度的双指数上界并证明其紧性，确定了球面特有的最优增长率 $2/\pi$ ，为二维欧拉动力学在曲面流形上的研究开辟了新方向。

4. 粘性β-平面方程的整体适定性与时间衰减估计

Global well-posedness and temporal decay estimates for the viscous $$β$$ -plane equations

作者: Tomoaki Yoshizawa

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19130v1

摘要翻译: 我们考虑粘性β-平面方程的柯西问题。首先，对初值数据相对于Rossby参数足够小的情形，建立了该系统的整体适定性。继而证明了流的正则化效应，并指出所得整体解在经典意义上满足该方程。我们还揭示：在若干附加假设下，解的衰减速率与对应线性解一致，且其渐近行为等价于线性化方程积分核的一个常数倍。特别地，该解的衰减速率快于由热方程流所预期的衰减速率。

Comment: 66 pages

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

1. 研究对象和背景

本文研究粘性β平面方程的Cauchy问题：

\begin{cases} \partial_t\omega + (u \cdot \nabla)\omega = \Delta\omega + \beta L_1\omega & \text{in } (0, \infty) \times \mathbb{R}^2, \\ u = \nabla^\perp(-\Delta)^{-1}\omega & \text{in } (0, \infty) \times \mathbb{R}^2, \\ \omega(0, \cdot) = \omega_0 & \text{in } \mathbb{R}^2, \end{cases}

其中 $\omega = \omega(t, x): (0, \infty)\times\mathbb{R}^2 \to \mathbb{R}$ 是未知标量涡度，算子 $L_1 := \partial_1(-\Delta)^{-1}$ 。该系统描述了受科里奥利力影响的二维流体动力学，其中科里奥利参数随纬度变化，Rossby参数 $\beta \in \mathbb{R}$ 表示科里奥利参数的变化率。通过Duhamel公式，该系统可重写为积分方程：

\omega(t) = T_\beta(t)\omega_0 - \int_0^t T_\beta(t-\tau) \text{div}(u(\tau)\omega(\tau)) d\tau,

其中 $T_\beta(t) = e^{t(\Delta+\beta L_1)}$ 为线性传播子。当 $\beta = 0$ 时，系统退化为二维Navier-Stokes方程的涡度形式，Giga-Miyakawa-Osada [GMO88]证明了初始涡度 $\omega_0$ 位于有限Radon测度空间时全局解存在，Gallagher-Gallay [GG05]补充了唯一性。值得注意的是，方程中同时包含耗散算子 $\Delta$ 和色散算子 $$L_1$$ ，导致线性化解可能比热方程预期衰减更快，这对非线性问题的全局适定性分析有益。然而， $\beta \neq 0$ 时的非线性结构 $(u \cdot \nabla)\omega = \text{div}(u\omega)$ 中，一个因子比另一个多至少一阶导数，而 $\dot{H}^1(\mathbb{R}^2) \not\hookrightarrow L^\infty(\mathbb{R}^2)$ 的嵌入失败导致处理非线性估计时面临关键困难。本文旨在填补以下空白：(1) 建立初始数据相对于Rossby参数足够小时的全局适定性框架；(2) 证明解的平滑效应并验证其满足经典意义下的方程；(3) 证明解衰减速度比热方程更快，且渐近行为与线性化解一致，特别是将衰减估计推广到所有 $2 \leq p \leq \infty$ 的 $$L^p$$ 范数范围，这在先前关于旋转流体的研究中尚未实现。

2. 主要定理或主要结果

定理1.1（全局适定性）：设 $\beta \neq 0$ ， $2 < p_1, p_2, r_1 < \infty$ ， $2 \leq r_2 < \infty$ ， $0 \leq \delta \leq 1/5$ 满足：

\delta < \min\left\{\frac{2}{p_1}, \frac{2}{p_2}\right\}, \quad \max\left\{1-\frac{1}{p_1}, 1-\frac{1}{p_2}\right\} \leq \frac{1}{p_1} + \frac{1}{p_2} - \frac{\delta}{2},

\max\left\{\frac{1}{2} - \frac{1}{p_1} + \frac{\delta}{2} - \frac{3}{2}\left(1-\frac{2}{p_2}\right), \frac{1}{2}\left(1-\frac{2}{p_1}\right)\right\} \leq \frac{1}{r_1} \leq \min\left\{\frac{1}{2} - \frac{1}{p_1} + \frac{\delta}{2}, \frac{5}{4}\left(1-\frac{2}{p_1}\right)\right\},

\frac{1}{2}\left(1-\frac{2}{p_1}\right) \leq \frac{1}{r_2} \leq 1 - \frac{1}{p_2} + \frac{\delta}{2},

且满足 $$r_2 > 2$$ 时 $\frac{1}{2}(1-2/p_2) \leq 1/r_2 \leq 5/4(1-2/p_2)$ ，或 $$r_2 = 2$$ 时 $\frac{1}{2}(1-2/p_2) < 1/2 < 5/4(1-2/p_2)$ 。则存在 $$C_0 > 0$$ ，当 $\omega_0 \in \dot{H}^{-1+\delta} \cap \dot{H}^\delta$ 满足

|\beta|^{-\delta/3} \|\omega_0\|_{\dot{H}^{-1+\delta}} + |\beta|^{-1+\delta/3} \|\omega_0\|_{\dot{H}^\delta} \leq C_0

时，积分方程(1.2)存在唯一解 $\omega \in C([0, \infty); \dot{H}^{-1+\delta}) \cap L^{r_1}([0, \infty); \dot{W}^{-1+\delta,p_1}) \cap L^{r_2}([0, \infty); \dot{W}^{\delta,p_2})$ 。

定理1.6（平滑效应）：在定理1.1假设下，若 $1/r_1 < 1/2 - 1/p_1 + \delta/2$ 且 $1/r_2 < 1 - 1/p_2 + \delta/2$ ，则对任意 $m \geq 0$ ，解 $\omega \in C((0, \infty); \dot{H}^{-1+\delta+m}) \cap C^1((0, \infty); \dot{H}^{\delta+m})$ 且对 $$t > 0$$ 给出经典解。若 $\omega_0 \in \dot{H}^{-1+\delta} \cap \dot{H}^{\delta+m}$ ，则 $\omega \in C([0, \infty); \dot{H}^{-1+\delta+m})$ 。

定理1.7（时间衰减与渐近行为）：设 $0 \leq \delta < 1/13$ ， $A = (\delta, p, r_1, p, r_2)$ 满足定理1.1假设及(1.6)式，存在 $$C_A > 0$$ ，当 $\omega_0 \in \dot{H}^{-1+\delta} \cap \dot{H}^{3-2/p} \cap L^1 \cap \dot{W}^{-1+\delta, p/(p-1)}$ 满足

|\beta|^{-\delta/3} \|\omega_0\|_{\dot{H}^{-1+\delta}} + |\beta|^{-1+\delta/3} \|\omega_0\|_{\dot{H}^\delta} \leq C_A,

则对任意 $s \geq 0$ ， $a \in [2, \infty]$ ，存在常数 $$C$$ （依赖于 $s, A, \|\omega_0\|_{L^1}, |\beta|^{-\delta/3} \|\omega_0\|_{\dot{H}^{-1+\delta}}$ 等），使得

\|\omega(t)\|_{\dot{W}^{s,a}} \leq C t^{-s/2 -1 + 1/a} \min\left\{1, |\beta|^{-1+2/a} t^{-3/2(1-2/a)}\right\}, \quad \forall t > 0.

若 $\frac{1}{2} - \frac{1}{p} + \frac{\delta}{2} - \frac{3}{2}(1-2/p) < \frac{1}{r_1}$ ，则渐近行为成立：

\lim_{t\to\infty} \left(t^{-s/2 -1 + 1/a} \min\left\{1, |\beta|^{-1+2/a} t^{-3/2(1-2/a)}\right\}\right)^{-1} \left\|\omega(t) - K_{\beta,t} \int_{\mathbb{R}^2} \omega_0(y) dy \right\|_{\dot{W}^{s,a}} = 0,

其中 $K_{\beta,t} = e^{t\beta L_1} G_t$ 为线性化解的积分核， $$G_t$$ 为Gauss核。

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

主要技术障碍源于非线性项 $(u \cdot \nabla)\omega = \text{div}(u\omega)$ 的特殊结构：一个因子比另一个多至少一阶导数，而 $\dot{H}^1(\mathbb{R}^2) \not\hookrightarrow L^\infty(\mathbb{R}^2)$ 导致标准嵌入失效，使得非线性估计面临临界困难。作者的创新点在于：(1) 通过直接估计积分方程而非微分方程的能量估计来证明平滑效应；(2) 引入双辅助时空范数 $L^{r_1}([0, \infty); \dot{W}^{-1+\delta,p_1})$ 与 $L^{r_2}([0, \infty); \dot{W}^{\delta,p_2})$ ，精确匹配非线性项中导数阶数的不匹配；(3) 利用Littlewood-Paley分解与色散估计获得线性传播子的精细时空估计。

关键估计包括：引理2.1证明了 $\|e^{t\beta L_1}P_k f\|_{L^\infty} \leq C 2^{3k}|\beta t|^{-1}\|P_k f\|_{L^1}$ ，揭示了色散效应导致的 $|\beta t|^{-1}$ 衰减因子。命题2.2-2.3建立了线性传播子在齐次Besov空间中的关键估计：

\|T_\beta(t)f\|_{\dot{B}^{s'}_{p,r}} \leq C t^{-\frac{s'-s}{2}-1+\frac{2}{p}} \min\left\{1, |\beta|^{-1+\frac{2}{p}} t^{-\frac{3}{2}(1-\frac{2}{p})}\right\} \|f\|_{\dot{B}^s_{p',r}},

其中 $\min$ 项体现了色散效应带来的额外衰减。命题2.5证明了Strichartz型估计：

\|T_\beta(t)f\|_{L^r([0,\infty); \dot{W}^{s,p})} \leq C |\beta|^{-\frac{1}{3}(\frac{2}{r} + \frac{2}{p} -1 -s)} \|f\|_{L^2},

这对处理非线性项至关重要。在证明定理1.1时，作者通过分解非线性项并应用Hölder不等式、Sobolev嵌入及上述线性估计，得到收缩映射所需的估计：

\left\| \int_0^t T_\beta(t-\tau) \text{div}(u\omega) d\tau \right\|_{X} \leq C (\|\omega\|_X^2 + \|\omega\|_X^3),

其中 $$X$$ 为构造的时空函数空间。对于衰减估计（定理1.7），作者通过比较线性解与非线性解，结合命题2.4的渐近展开：

\left\| T_\beta(t)\omega_0 - K_{\beta,t} \int \omega_0 \right\|_{\dot{W}^{s,a}} \lesssim t^{-\frac{s}{2}-1+\frac{1}{a}} \min\left\{1, |\beta|^{-1+\frac{2}{a}} t^{-\frac{3}{2}(1-\frac{2}{a})}\right\},

并利用非线性项的衰减性质完成证明。特别地，通过 $W^{s+4-4/a,1}$ 范数控制 $\|G_1(\cdot - t^{-1/2}y) - G_1\|_{B^{s+4-4/a}_{1,\infty}}$ ，确保了渐近行为的严格成立。

4. 与之前工作的比较

与二维Navier-Stokes方程比较：当 $\beta=0$ 时，Giga-Miyakawa-Osada [GMO88]证明了有限Radon测度初值的全局适定性且无需小性条件，而本文需施加 $|\beta|^{-\delta/3} \|\omega_0\|_{\dot{H}^{-1+\delta}} + |\beta|^{-1+\delta/3} \|\omega_0\|_{\dot{H}^\delta} \leq C_0$ 的小性条件（定理1.1），这是利用 $$L_1$$ 色散效应的必要代价。衰减方面， $\beta=0$ 时解满足 $\|\partial_t^\ell \partial_x^\alpha \omega(t)\|_{L^p} \leq C t^{-\ell - |\alpha|/2 -1 + 1/p}$ （热方程速率），而 $\beta \neq 0$ 时获得额外衰减因子 $\min\{1, |\beta|^{-1+2/a} t^{-3/2(1-2/a)}\}$ （定理1.7），显著快于热方程。

与三维旋转Navier-Stokes方程比较：Chemin等 [CDGG02, CDGG06]证明当 $|\Omega| \geq C$ 时 $L^2 + H^{1/2}$ 初值存在全局解；Iwabuchi-Takada [IT13]将范围扩展至 $\dot{H}^s$ （ $$1/2 < s < 9/10$$ ）且 $\|u_0\|_{\dot{H}^s} \leq C|\Omega|^{s/2 - 1/4}$ 。关键差异在于非线性结构：三维情况下散度自由条件允许将非线性项重写为 $\text{div}(u \otimes u)$ （无导数乘积），而β平面方程中为 $\text{div}(u\omega)$ （一阶导数差异），导致本文需引入额外正则性 $\omega_0 \in \dot{H}^\delta$ 及双辅助范数（见Remark 1.3）。Koh-Lee-Takada [KLT14]将 $$s$$ 范围扩展至 $$1/2 < s < 9/10$$ ，Ahn-Kim-Lee [AKL22]进一步推广至分数阶Laplacian情形，但均受限于三维空间特性。

与已有衰减结果比较：Ahn-Kim-Lee [AKL22]对三维旋转Navier-Stokes方程得到 $\|u(t)\|_{L^p} \leq C \|u_0\|_{L^q} t^{-3/2(1/q - 1/p)}(1 + |\Omega|t)^{-1 + 2/p}$ ，但仅适用于 $1 < q \leq p/(p-1) \leq 2 \leq p < 3$ ；Egashira-Takada [ET23]证明当 $u_0 \in L^1 \cap \dot{H}^s$ （ $$1/2 < s < 9/10$$ ）时， $\|u(t)\|_{L^p} \leq C t^{-3/2(1-1/p)}(1 + |\Omega|t)^{-1 + 2/p}$ ，但 $p \in [2, p^*]$ 且 $$p^* < 3$$ ；Kim [Kim22]对MHD方程得到 $\|u(t)\|_{L^p} \leq C t^{-3/2(1-1/p)-1+2/p} |\Omega|^{-1+2/p}$ ，同样限于 $$p < 3$$ 。本文突破性在于：(1) 衰减估计适用于所有 $2 \leq p \leq \infty$ （含 $L^\infty$ ），而先前结果均要求 $$p < 3$$ ；(2) 允许任意高阶导数（ $s \geq 0$ ），而现有文献通常限于低阶导数；(3) 严格证明渐近行为 $\omega(t) \sim K_{\beta,t} \int \omega_0$ ，而非仅衰减估计。代价方面：需小性条件（二维Navier-Stokes无需），且难以扩展至时间导数（Remark 1.8指出 $\partial_t \omega$ 涉及 $L_1\omega$ ，其衰减更慢）。此外，本文结果依赖于 $\beta \neq 0$ 的色散效应，而三维旋转流体研究中 $|\Omega|$ 大时的适定性源于科里奥利力的稳定化作用，机制不同。

5. Vlasov(-Fokker-Planck)-Navier-Stokes 方程组的定量粒子系统逼近

Quantitative approximation of the Vlasov(-Fokker-Planck)-Navier-Stokes system by stochastic particle systems

作者: Ludovic Goudenège, Christian Olivera, Gabriela Planas, Alexandre Richard

类别: math.PR, math.AP

发布日期: 2026-04-20

arXiv ID: 2604.18871v1

摘要翻译: 本文研究一类由不可压 Navier-Stokes 方程与 Vlasov(-Fokker-Planck) 方程通过拖曳力耦合而成的流体-粒子系统；此类模型自然出现于气溶胶、喷雾以及更一般的两相流研究中。在空间维数 $d\in \{2,3\}$ 下，我们建立了由 $$N$$ 个相互作用的随机粒子与流体耦合所构成的系统，在 $N\to \infty$ 时向 Vlasov(-Fokker-Planck)-Navier-Stokes 系统收敛的定量收敛速率。特别地，考虑粒子所受噪声强度随 $N\to \infty$ 而消失的情形，此时对应于 Vlasov-Navier-Stokes 系统。更精确地，我们证明：该粒子系统的经验测度收敛于耦合系统的 Vlasov(-Fokker-Planck) 分量，而流体速度场收敛于其 Navier-Stokes 分量。证明结合了随机分析与偏微分方程方法，对离散系统和连续系统分别建立能量估计与交换子估计。

Comment: 无

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

1. 研究对象和背景

本文研究不可压缩Navier-Stokes方程与Vlasov(-Fokker-Planck)方程通过拖曳力耦合的流体-粒子系统，该模型源于气溶胶、喷雾及两相流研究。系统在环面 $\mathbb{T}^d$ （ $d \in \{2,3\}$ ）上定义，包含流体速度 $$u$$ 、压力 $$p$$ 和相空间分布函数 $$F$$ ，其耦合项 $$(u_t - v)F_t$$ 基于Stokes定律。微观模型由 $$N$$ 个粒子的随机微分方程组描述，粒子位置 $X^{i,N}_t$ 与速度 $V^{i,N}_t$ 通过截断流体速度 $\chi_A(u^N_t)$ 与噪声 $\sigma_N dB^i_t$ 演化。目标是证明经验测度 $F^N_t = S^N_t \ast \vartheta^N$ （ $S^N_t = \frac{1}{N}\sum_{i=1}^N \delta_{(X^{i,N}_t, V^{i,N}_t)}$ ）与流体速度 $$u^N_t$$ 在 $N \to \infty$ 时收敛至宏观系统解 $$F$$ 与 $$u$$ 。尽管系统(1.1)的全局弱解存在性已在 $$d=2,3$$ 维由[5,9,13,34]建立，且Flandoli等[13]在 $$d=2$$ 、 $\sigma > 0$ 下证明了经验测度收敛，但尚未有工作建立定量收敛率，亦未覆盖 $\sigma \geq 0$ 全范围及强混沌传播的轨迹收敛分析。

2. 主要定理或主要结果

在假设H1和H2下，特别地 $p \in (d, +\infty)$ ， $\gamma \in (d/p, 1)$ ，整数 $k \geq 3$ ，且截断参数 $A \geq \|u\|_{L^\infty_{t,x}}$ ，若 $$d = 3$$ 则需满足 $\gamma/2 + 3/2(1/2 - 1/p) < 1$ ，有以下结果：

• 定理2.1：当 $\sigma > 0$ 时，对任意 $q \in [1, +\infty)$ ，存在 $C = C(q, d, \gamma, p, A, \sigma) > 0$ ，使得对任意 $N \in \mathbb{N}^*$ ，

\left(\mathbb{E} \sup_{t \in [0,T]} \|u^N_t - u_t\|^{2q}_{\gamma,p}\right)^{1/q} + \left(\mathbb{E} \sup_{t \in [0,T]} \|\langle v \rangle^k (F^N_t - F_t)\|^{2q}_{L^2_{x,v}}\right)^{1/q}

\leq C \left(\left(\mathbb{E} \|u^{\circ,N} - u^\circ\|^{2q}_{\gamma,p}\right)^{1/q} + \left(\mathbb{E} \|\langle v \rangle^k (F^N_0 - \vartheta^N \ast F^\circ)\|^{2q}_{L^2_{x,v}}\right)^{1/q} + N^{-2(\alpha-\beta)} + N^{-2\beta(\gamma-d/p)} + N^{-(1/2 - d\beta - (d+1)\alpha)} + \tilde{\rho}^2_N\right) + 2\rho^2_N,

其中 $\rho^2_N = \sup_{t \in [0,T]} \|u^{(N)}_t - u_t\|^2_{\gamma,p} + \sup_{t \in [0,T]} \|\langle v \rangle^k (F^{(N)}_t \ast \vartheta^N - F_t)\|^2_{L^2_{x,v}}$ ， $\tilde{\rho}^2_N = \sup_{t \in [0,T]} \|\langle v \rangle^k (F^{(N)}_t \ast \vartheta^N - F^{(N)}_t)\|^2_{L^2_{x,v}}$ 。当 $\sigma = 0$ 时，对任意 $q \in [1, +\infty)$ 和 $\varepsilon > 0$ ，存在 $C = C(q, d, \gamma, p, A, \varepsilon) > 0$ ，使得

\left(\mathbb{E} \sup_{t \in [0,T]} \|u^N_t - u_t\|^{2q}_{\gamma,p}\right)^{1/q} + \left(\mathbb{E} \sup_{t \in [0,T]} \|\langle v \rangle^k (F^N_t - F_t)\|^{2q}_{L^2_{x,v}}\right)^{1/q}

\leq C \left(\cdots + N^{-2\beta(\gamma-d/p)+\varepsilon} + N^{-(1/2 - d\beta - (d+1)\alpha)}\right) \frac{1}{1-\tilde{\gamma}} e^{CT(1+\sigma^{-2}_N + \sigma^{-1/(1-\tilde{\gamma})}_N)} + 2\rho^2_N,

其中 $\tilde{\gamma} = \frac{1}{2}(\gamma + d(1/2 - 1/p)) \vee \gamma + 1/2 \in (0, 1)$ 。

• 定理2.5：在相同假设下，若 $\sup_{t \in [0,T]} \|m_0(F_t)\|_{L^\infty_x} < \infty$ ，则有能量范数估计：

\left(\mathbb{E} \sup_{t \in [0,T]} \|u^N_t - u_t\|^{2q}_{L^2_x}\right)^{1/q} + \left(\mathbb{E} \left\|\int_0^T \|\nabla(u^N_s - u_s)\|^2_{L^2_x} ds\right\|^q\right)^{1/q}

\leq C \left(\cdots \right) + C\rho^2_N,

范数替换为 $$L^2_x$$ ，具体形式与定理2.1类似。

• 定理2.7：在相同假设下，粒子轨迹的强混沌传播满足：

\mathbb{E} \sup_{t \in [0,T]} \left(|X^{i,N}_t - X^i_t| + |V^{i,N}_t - V^i_t|\right) \leq |\sigma - \sigma_N| + \mathbb{E} \sup_{t \in [0,T]} \|u^N_t - u_t\|_{\gamma,p},

其中 $$(X^i_t, V^i_t)$$ 为极限流体中的独立粒子轨迹。

根据注释2.2，收敛率具体为：当 $\sigma > 0$ 时， $\rho_N \lesssim |\sigma_N - \sigma|^{1/2} + N^{-\frac{1}{2}\beta(\gamma \wedge 1)} + N^{-\beta(\gamma - d/p)} + N^{-(\alpha - \beta)}$ ；当 $\sigma = 0$ 时， $\rho_N \lesssim \sigma_N^{1/2}$ 。注释2.4指出， $\sigma > 0$ 时收敛率为代数阶 $N^{-\min\{2(\alpha-\beta), 2\beta(\gamma-d/p), 1/2 - d\beta - (d+1)\alpha\}}$ ，而 $\sigma = 0$ 时因 $\sigma_N \gtrsim (\log N)^{-1/2}$ （由(2.11)），收敛率为对数阶 $(\log N)^{-1/2}$ 。

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

证明的核心障碍在于流体与粒子系统的强耦合：粒子运动依赖于流体速度，而流体方程又受粒子分布反馈。作者的创新框架包含三个关键环节：

辅助PDE系统的构造与分析：引入辅助系统 $(u^{(N)}, F^{(N)})$ ，其为系统(1.1)中 $\sigma$ 替换为 $\sigma_N$ 且对 $$u$$ 施加截断 $\chi_A$ 的解。关键引理5.1证明其适定性，引理5.2建立 $F^{(N)}$ 的矩估计 $\sup_{t \in [0,T]} M_{p+d(p-1)}(F^{(N)}_t) < \infty$ 。在 $\sigma > 0$ 情形（第5.4节），利用Fokker-Planck方程的抛物正则性，引理5.8证明 $\|\langle v \rangle^k (F^{(N)}_t - F_t)\|_{L^2_{x,v}} \lesssim |\sigma_N - \sigma|^{1/2}$ ；在 $\sigma = 0$ 情形（第5.5节），通过能量估计与引理5.13对 $\|\nabla_v F\|_{L^2_{x,v}}$ 的控制，得到 $\rho_N \lesssim \sigma_N^{1/2}$ 。
离散-连续系统的误差分解：定义 $\tilde{F}^N_t = F^{(N)}_t \ast \vartheta^N$ ，将 $$F^N_t - F_t$$ 分解为 $(F^N_t - \tilde{F}^N_t) + (\tilde{F}^N_t - F^{(N)}_t) + (F^{(N)}_t - F_t)$ 。关键估计（第5.3节）通过Itô公式与交换子技术，建立：
$\mathbb{E} \sup_{t \in [0,T]} \|\langle v \rangle^k (F^N_t - \tilde{F}^N_t)\|^{2q}_{L^2_{x,v}} \lesssim N^{-2(\alpha-\beta)} + N^{-2\beta(\gamma-d/p)} + N^{-(1/2 - d\beta - (d+1)\alpha)} + \mathbb{E} \int_0^T \|u^N_s - u^{(N)}_s\|^{2q}_{\gamma,p} ds,$
其中 $N^{-(1/2 - d\beta - (d+1)\alpha)}$ 来自鞅项收敛率， $N^{-2(\alpha-\beta)}$ 源于粒子速度的光滑化误差。
Grönwall论证的精细设计：在第6节，将 $u^N - u^{(N)}$ 表示为Leray投影下的积分形式，关键创新在于处理交叉项 $\langle F^N_s - F^{(N)}_s, u^N_s - v \rangle$ ：利用Bessel空间嵌入 $H^\gamma_p \hookrightarrow L^\infty$ （因 $\gamma > d/p$ ）与引理4.3的先验界，结合交换子估计
$\left\| \int_0^t \nabla e^{(t-s)\Delta} \mathbb{P} \langle F^N_s - \tilde{F}^N_s, u^N_s - v \rangle ds \right\|_{\gamma,p} \lesssim \int_0^t (t-s)^{-\frac{1}{2}(\gamma + d(1-1/p))} \|\langle v \rangle^k (F^N_s - \tilde{F}^N_s)\|_{L^2_{x,v}} ds,$
将粒子误差转化为流体源项。最终通过Grönwall不等式 $\|u^N - u^{(N)}\|_{\gamma,p} \lesssim \eta_N + \int_0^t \|u^N - u^{(N)}\|_{\gamma,p} ds$ 闭合估计，其中 $\eta_N$ 包含初始误差与 $$N$$ 相关项。特别地，在 $\sigma = 0$ 情形，因缺乏抛物正则性，引入 $\tilde{\gamma} = \frac{1}{2}(\gamma + d(1/2 - 1/p)) \vee \gamma + 1/2$ ，并利用 $\sigma_N^{-2}$ 增长控制，通过引理6.2得到 $\|u^N - u\|_{\gamma,p} \lesssim e^{C T \sigma_N^{-2}} (\cdots)$ ，结合 $\sigma_N \gtrsim (\log N)^{-1/2}$ 导出对数收敛率。

4. 与之前工作的比较

与既有文献相比，本文在三个维度实现突破：

适用范围的扩展：

与Flandoli等[13]（仅 $$d=2$$ ， $\sigma > 0$ ）相比，本文覆盖 $d \in \{2,3\}$ 全维度与 $\sigma \geq 0$ 全范围。特别地， $\sigma = 0$ 对应Vlasov-Navier-Stokes系统，其缺乏扩散正则性，现有文献[20]仅建立弱解存在性，本文首次给出微观推导的定量分析。
与[16]（ $$d=3$$ 光滑解存在性）和[10]（唯一性）相比，本文不依赖小初值假设，仅需H1中的正则性条件（如 $u^\circ \in H^\gamma_p$ ， $\gamma > d/p$ ），这比 $$H^1$$ 正则性更弱。

收敛率的精确刻画：

与[17,18]（PDE到PDE收敛）不同，本文处理从随机粒子系统到PDE的收敛。在 $\sigma > 0$ 情形，收敛率 $N^{-\min\{2(\alpha-\beta), 2\beta(\gamma-d/p), 1/2 - d\beta - (d+1)\alpha\}}$ 由三部分主导： $N^{-2(\alpha-\beta)}$ 源于速度空间光滑化（ $\alpha > \beta$ 保证主导）， $N^{-2\beta(\gamma-d/p)}$ 来自流体Bessel范数逼近， $N^{-(1/2 - d\beta - (d+1)\alpha)}$ 为鞅项随机收敛率（要求 $d\beta + (d+1)\alpha < 1/2$ ），优于[26]中直接粒子相互作用的 $N^{-1/2}$ ，因间接耦合减弱了相关性。
在 $\sigma = 0$ 情形，收敛率退化为对数阶 $(\log N)^{-1/2}$ ，这是本质的：由假设(2.11) $\lim_{N\to\infty} \sqrt{\log N} \sigma_N = +\infty$ ，若 $\sigma_N$ 衰减过快（如 $N^{-\varepsilon}$ ），则辅助系统 $(u^{(N)}, F^{(N)})$ 无法逼近 $$(u, F)$$ ，与[20]中 $\sigma = 0$ 的弱解理论一致——缺乏正则性导致定量分析更困难。

技术框架的创新：

相比[1,8,11]的均匀化方法（处理周期性介质中的流体），本文采用随机粒子系统方法，更贴近物理模型。关键区别在于：均匀化推导Brinkman项 $-\mu u$ ，而本文得到非线性项 $-\langle F, u - v \rangle$ ，后者包含粒子分布的动态反馈。
与[24,26]的适度相互作用粒子系统不同，本文粒子仅通过流体间接耦合（无直接相互作用），导致标准传播混沌技术失效。作者通过引入辅助PDE和误差分解，将问题转化为确定性逼近与随机逼近的耦合，这一框架可推广至其他间接耦合系统。

代价与局限：

截断函数 $\chi_A$ 的引入是技术必需（见第4节注释），因无截断时系统(1.2)的适定性未解决。这限制了结果对大速度流体的适用性，但物理上 $A \geq \|u\|_{L^\infty_{t,x}}$ 合理。
$$d=3$$ 时要求 $\gamma/2 + 3/2(1/2 - 1/p) < 1$ ，这是三维Navier-Stokes方程能量估计的典型限制（类似[32]中的条件）。若 $$p > 3$$ ， $\gamma < 1$ 自动满足；若 $p \to 3^+$ ，则 $\gamma < 1/2$ ，弱于 $$H^1$$ 正则性。
$\sigma = 0$ 情形的对数收敛率可能最优：若 $\sigma_N = 0$ （即无噪声粒子系统），则强混沌传播不成立（粒子轨迹确定性依赖），必须保留 $\sigma_N > 0$ 且缓慢衰减。

中文标题: Hartree–Fock 均匀平衡态在 $\mathbb{R}^d$ 中的渐近稳定性

摘要翻译: 本文在全空间 $\mathbb{R}^d$ （ $d \ge 3$ ）中，针对含 { em 非对角交换算子} 的含时 Hartree–Fock 方程，建立了大类平移不变稳态解的非线性 Landau 阻尼与渐近稳定性。该交换算子自然出现于大费米子系统的平均场理论中。尽管交换项仅为低阶项，其引入却扰动了经典 Schrödinger 色散结构，并导致背景电子对空间密度产生复杂的线性响应：其色散关系不再如经典 Vlasov 或 Hartree 理论中那样表现为 Fourier 乘子。此外，每一基本波的群速度均耦合所有其他 Fourier 模式，从而引发精细的 { em 动量依赖回波共振}。为克服上述困难，我们发展了一种非线性迭代格式：该格式依赖于细致的预解式分析，利用 Fourier 空间中的输运型色散结构，并在加权范数 $L^\infty_{k,p}$ 中传播相混合与 Landau 阻尼。

不相关论文（基本信息）

1. 球面Baouendi-Grushin抛物方程的零控最小时间

Minimal time for null controllability of the parabolic spherical Baouendi-Grushin equation

作者: Cyprien Tamekue

类别: math.AP, math.OC

发布日期: 2026-04-21

arXiv ID: 2604.19700v1

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

摘要: We study null controllability for the parabolic equation on $\mathbb{S}^{2}$ endowed with its canonical almost-Riemannian structure. For a spherical crown $ω=\{α<x_3<β\}$ , where $0\le α<β\le1$ , we prove the sharp minimal time formula $T_{\min}(ω)=\ln(1/\sqrt{1-α^{2}})$ for null controllability in $$ω$$ . We also prove that, whenever the control region contains the equator, null controllability holds in every positive time. The proof combines two complementary tools. First, after Fourier decomposition with respect to the periodic variable, we establish observability estimates for a family of one-dimensional singular parabolic equations, with constants uniform with respect to the Fourier mode; the singularities at the poles are handled via a Hardy-Poincaré inequality. Second, for crowns away from the equator, we use the moment method to construct controls on the pole-touching crown $$α<x_3< 1$$ from sharp weighted lower bounds on associated Legendre functions, and then pass to a general crown $$α<x_3<β$$ by a cut-off argument on the full domain combined with the arbitrary-time controllability of crowns containing the equator. The result closes the large-time gap left in earlier work and gives the exact null-controllability threshold for the canonical almost-Riemannian heat equation on $\mathbb S^2$ .

2. 生命系统的多尺度动力学结构

Multiscale Kinetic Structures for Living Systems

作者: Diletta Burini, Damian A. Knopoff

类别: math.DS, math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19659v1

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

摘要: This paper develops a conceptual extension of the Kinetic Theory of Active Particles, building upon the framework introduced in [2]. Living systems cannot be adequately described within classical single-scale paradigms, even when refined. To overcome this limitation, we introduce a Multiscale Kinetic Theory of Active Particles (MS-KTAP), in which a sub-microscopic scale of interacting entities is incorporated into the description of collective dynamics. In this framework, the activity variable is interpreted as an emergent quantity arising from lower-scale regulatory mechanisms and influenced by interactions across higher scales. The proposed framework captures key features of living systems, including heterogeneity, adaptive decision-making, nonlinear and non-conservative interactions, spatial dynamics, and cross-scale feedback, within a unified mathematical structure. Competition and cooperation are thus described across multiple levels of organization. The first part of the paper derives the mathematical framework, while the second presents how specific models can be obtained. The paper concludes with perspectives on further developments, including possible integrations with scientific machine learning.

3. 极小图增长的上界估计

An upper bound on the growth of minimal graphs

作者: Allen Weitsman

类别: math.DG, math.AP, math.CV

发布日期: 2026-04-21

arXiv ID: 2604.19630v1

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

摘要: Graphs of solutions to the minimal surface equation over simply connected domains with boundary values 0 can have at most exponential growth.

4. 强锚定条件下带Flory-Huggins势的Cahn-Hilliard能量泛函的极小元

Minimizers for the Cahn-Hilliard energy functional with the Flory-Huggins potential under strong anchoring conditions

作者: Shibin Dai, Abba Ramadan, Natasha Sharma

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19590v1

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

摘要: In this paper, we theoretically and numerically study the minimizers for the Cahn-Hilliard energy with the Flory-Huggins potential under the strong anchoring condition, i.e., the Dirichlet boundary condition. We reveal bifurcation phenomena mediated by the boundary condition, the transition layer thickness, and the temperature of the system. Numerical simulations are also presented to approximate the minimizers of this energy by solving a gradient-flow equation, namely the Allen-Cahn equation constrained with strong anchoring conditions and random initial data. The effects of varying the transition layer thickness and temperature are presented to confirm the theoretical analysis.

5. 相分离模型中的奇异性：非局部Cahn-Hilliard方程的谱元法

Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation

作者: Andrés Miniguano-Trujillo, Andrea Poiatti, Maurizio Grasselli, Benjamin Goddard, John Pearson

类别: math.NA, math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19521v1

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

摘要: The nonlocal Cahn-Hilliard equation provides a natural extension of the classical model for phase separation by incorporating long-range interactions through a singular convolution kernel. While this formulation admits a rich existence and regularity theory, its numerical approximation remains challenging: discretisation of the nonlocal term leads to dense operators, and the singularity of the kernel requires special treatment in collocation-based schemes. In this work, we develop an efficient and error-controlled numerical framework for the nonlocal Cahn-Hilliard system with constant mobility, logarithmic potential, Newtonian interaction kernel, and no-flux boundary conditions. Our approach is based on a pseudospectral multishape method that accurately approximates the action of singular convolution operators. We present high-resolution numerical solutions for this nonlocal system of equations that can be achieved with limited computational resources.

6. 具空间非均匀非线性的半线性热方程的稳定爆破型貌

Stable blowup profile for a semilinear Heat Equation with spatially inhomogeneous nonlinearity

作者: Irfan Glogić, Sarah Kistner, Birgit Schörkhuber

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19389v1

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

摘要: We study the focusing semilinear heat equation with an additional defocusing Hénon-type nonlinearity, the coupling of which is measured by a constant $$c >0$$ . For $c \in (0,c^*)$ , the model admits a closed-form self-similar blowup solution in every space dimension $d \geq 1$ . Restricting ourselves to the three-dimensional case, we study the stability of this solution under small non-radial perturbations. By working in intersection Sobolev spaces with additional angular regularity, we prove finite co-dimension stability for all admissible values of $$c$$ . Furthermore, we analyze the spectrum of the underlying linearized operator and we prove stable blowup for the cubic-quintic case and $$c$$ sufficiently close to $$c^*$$ . Finally, we discuss the situation for small values of $$c$$ and use a modified version of the classical GGMT criterion to give an upper bound on the number of unstable eigenvalues.

7. 具非局部相互作用的等周问题中的分岔

Bifurcations in Isoperimetric Problems with Nonlocal Interactions

作者: Fabio De Regibus, Massimo Grossi, Monica Musso

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19170v1

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

摘要: We study isoperimetric problems modeled on the liquid drop model, with nonlocal interactions under a volume constraint. While balls are natural critical points, we show that, for an unbounded sequence of radii, non-spherical solutions bifurcate from the family of balls. These new solutions lie arbitrarily close to balls and can have arbitrarily large volume. Conversely, at radii outside this sequence, no bifurcation occurs, and nearby solutions are trivial, arising only from rigid motions.

8. 具有给定梯度像的平均曲率型超曲面

The mean curvature type hypersurfaces with prescribed gradient image

作者: Jiguang Bao, Rongli Huang, Qinfeng Jiang

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19020v1

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

摘要: In this paper, we consider the existence of mean curvature type hypersurfaces with prescribed gradient image. Let $$Ω$$ and $\tildeΩ$ be uniformly convex bounded domains in $\mathbb{R}^n$ with smooth boundary. We show that there exists unique convex solutions for the second boundary value problem of mean curvature type equations.

9. 关于拉格朗日平均曲率型方程的第二边值问题

On the second boundary value problem for Lagrangian mean curvature type equation

作者: Jiguang Bao, Qinfeng Jiang

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19006v1

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

摘要: This article is concerned with the second boundary value problem of the Lagrangian mean curvature type equation arising from special Lagrangian geometry. By the parabolic method, we consider a fully nonlinear parabolic equation with oblique derivative boundary condition, and show the long time existence and convergence of the flow. It follows that the existence and uniqueness of the smooth uniformly convex solution are obtained, which generalizes the Brendle–Warren’s theorem about minimal Lagrangian diffeomorphism in Euclidean metric space.

10. 具有对称性的四边形上第二阶诺伊曼特征函数的临界点

Critical points of the second Neumann eigenfunctions on the quadrangles with symmetry

作者: Haiyun Deng, Changfeng Gui, Xuyong Jiang, Xiaoping Yang, Ruofei Yao, Jun Zou

类别: math.AP

发布日期: 2026-04-21

arXiv ID: 2604.19003v1

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

摘要: In this paper, we focus primarily on the symmetry properties of the second Neumann eigenfunction $$u$$ with respect to the symmetry axis or symmetry center of the relevant domain $$Q$$ , such as isosceles trapezoids, parallelograms, kite domains, and we provide some affirmative answers to the Hot Spots Conjecture for these domains. Our proofs combine symmetry decomposition, comparison of eigenvalues, and the continuity method. Precisely, we have the following three aspects of results. (1) when $$Q$$ is an isosceles trapezoid, if the base angle $α\le \fracπ{3}$ , $$u$$ is antisymmetric about the symmetric axis; if the base angle $α> \fracπ{3}$ , there exists a critical height $\hat{h}$ , when height $h<\hat{h}$ , $$u$$ is antisymmetric about the symmetric axis; when height $h>\hat{h}$ , $$u$$ is symmetric about the symmetric axis; when height $h=\hat{h}$ , the multiplicity of second Neumann eigenvalue is 2. Meanwhile, we fully characterize the location of non-vertex critical points of $$u$$ on $\overline{Q}$ . (2) When $$Q$$ is a parallelogram, $$u$$ is centrally antisymmetric about the center of $$Q$$ and does not have any non-vertex critical points. In particular, when $$Q$$ is a rhombus, $$u$$ is symmetric with respect to the longer diagonal and is antisymmetric with respect to the short diagonal. (3) When $$Q$$ is a kite $$P_1P_2P_3P_4$$ , where $$P_1$$ is the origin, $$P_2=(a,-h)$$ lies in the four quadrant, $$P_3=(1,0)$$ lies on the positive $$x$$ -axis, and $$P_4=(a,h)$$ is symmetric with $$P_2$$ about $$x$$ -axis which lies in the first quadrant. If $a\ge 2$ , $$u$$ is antisymmetric about $$x$$ -axis; if $$0<a<2$$ , there exist two constants $$h_0$$ and $$h_1$$ ( $h_0\le h_1$ ), when $$h<h_0$$ , $$u$$ is symmetric about $$x$$ -axis; when $$h>h_1$$ , $$u$$ is antisymmetric about $$x$$ -axis. Meanwhile, we fully characterize the location of non-vertex critical points of $$u$$ on $\overline{Q}$ .

11. 分层Anderson-Bernoulli模型在 $\bZ^d$ 上的Anderson局域化

Anderson Localization for the hierarchical Anderson-Bernoulli model on $\mathbb{Z}^d$

作者: Shihe Liu, Yunfeng Shi, Zhifei Zhang

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

发布日期: 2026-04-21

arXiv ID: 2604.18989v1

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

摘要: In this paper, we prove Anderson localization for a hierarchical Anderson-Bernoulli model on lattice with arbitrary dimension, where the potential is characterized by a geometric hierarchical structure combined with fluctuations induced by independent and identically distributed (i.i.d.) Bernoulli random variables. Our method is also applicable to proving a probabilistic unique continuation result on $\mathbb{Z}^d$ .

12. 非局部PDE–ODE细菌持留细胞模型的持续性阈值分析

Analysis of persistence thresholds for a nonlocal PDE–ODE model of bacterial persister cells

作者: Chongming Li, Tyler Meadows, Troy Day

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

发布日期: 2026-04-20

arXiv ID: 2604.18809v1

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

摘要: Within many bacterial colonies, persister cells exist as a subpopulation that is tolerant to antibiotics and other stressors, yet not genetically distinct from the rest of the colony. A recent study has proposed epigenetic inheritance as a mechanism that leads to the presence of persister cells. We analyze a nonlocal PDE–ODE model introduced in that study to describe the epigenetic inheritance process and establish its mathematical well-posedness, including existence, uniqueness, and nonnegativity of solutions. We identify a sharp parameter threshold delineating extinction from persistence of the colony: below this threshold the washout equilibrium is globally asymptotically stable, while above it a unique positive equilibrium exists and the population is weakly persistent. Notably, this threshold is independent of the internal community structure.

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

1. 广义SQG方程在无粘极限下的哈密顿紧性与耗散#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

2. 二维重力-毛细水波中带涡度的小尺度生成#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

3. 球面上欧拉方程的涡量梯度增长#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

4. 粘性β-平面方程的整体适定性与时间衰减估计#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

5. Vlasov(-Fokker-Planck)-Navier-Stokes 方程组的定量粒子系统逼近#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

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

1. 全局各向异性微局部分析中的奇点滤子#

摘要翻译#

2. 由电荷交换驱动的粒子动力学#

摘要翻译#

3. 具有变号扩散系数的抛物型方程的激波波前#

摘要翻译#

4. Allard 正则性理论在阿贝尔 Yang–Mills–Higgs 方程中的应用#

摘要翻译#

5. 二维空间上带Rashba型自旋-轨道耦合的三次非线性薛定谔系统的小微涡解的存在性#

摘要翻译#

6. 二维薛定谔算子弱耦合现象再探#

摘要翻译#

7. Orlicz框架下的随机Σ-收敛及其应用#

摘要翻译#

8. 非椭圆Fredholm理论讲义#

摘要翻译#

9. Hartree–Fock 均匀平衡态在 $\mathbb{R}^d$ 中的渐近稳定性#

摘要翻译#

不相关论文（基本信息）#

1. 球面Baouendi-Grushin抛物方程的零控最小时间#

2. 生命系统的多尺度动力学结构#

3. 极小图增长的上界估计#

4. 强锚定条件下带Flory-Huggins势的Cahn-Hilliard能量泛函的极小元#

5. 相分离模型中的奇异性：非局部Cahn-Hilliard方程的谱元法#

6. 具空间非均匀非线性的半线性热方程的稳定爆破型貌#

7. 具非局部相互作用的等周问题中的分岔#

8. 具有给定梯度像的平均曲率型超曲面#

9. 关于拉格朗日平均曲率型方程的第二边值问题#

10. 具有对称性的四边形上第二阶诺伊曼特征函数的临界点#

11. 分层Anderson-Bernoulli模型在$\bZ^d$上的Anderson局域化#

12. 非局部PDE–ODE细菌持留细胞模型的持续性阈值分析#

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

1. 广义SQG方程在无粘极限下的哈密顿紧性与耗散

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 二维重力-毛细水波中带涡度的小尺度生成

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 球面上欧拉方程的涡量梯度增长

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 粘性β-平面方程的整体适定性与时间衰减估计

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. Vlasov(-Fokker-Planck)-Navier-Stokes 方程组的定量粒子系统逼近

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 全局各向异性微局部分析中的奇点滤子

摘要翻译

2. 由电荷交换驱动的粒子动力学

摘要翻译

3. 具有变号扩散系数的抛物型方程的激波波前

摘要翻译

4. Allard 正则性理论在阿贝尔 Yang–Mills–Higgs 方程中的应用

摘要翻译

5. 二维空间上带Rashba型自旋-轨道耦合的三次非线性薛定谔系统的小微涡解的存在性

摘要翻译

6. 二维薛定谔算子弱耦合现象再探

摘要翻译

7. Orlicz框架下的随机Σ-收敛及其应用

摘要翻译

8. 非椭圆Fredholm理论讲义

摘要翻译

9. Hartree–Fock 均匀平衡态在 $\mathbb{R}^d$ 中的渐近稳定性

摘要翻译

不相关论文（基本信息）

1. 球面Baouendi-Grushin抛物方程的零控最小时间

2. 生命系统的多尺度动力学结构

3. 极小图增长的上界估计

4. 强锚定条件下带Flory-Huggins势的Cahn-Hilliard能量泛函的极小元

5. 相分离模型中的奇异性：非局部Cahn-Hilliard方程的谱元法

6. 具空间非均匀非线性的半线性热方程的稳定爆破型貌

7. 具非局部相互作用的等周问题中的分岔

8. 具有给定梯度像的平均曲率型超曲面

9. 关于拉格朗日平均曲率型方程的第二边值问题

10. 具有对称性的四边形上第二阶诺伊曼特征函数的临界点

11. 分层Anderson-Bernoulli模型在 $\bZ^d$ 上的Anderson局域化

12. 非局部PDE–ODE细菌持留细胞模型的持续性阈值分析