生成时间: 2026年03月19日 04:39:27

重点关注论文数量: 8

了解领域论文数量: 7

不相关论文数量: 8

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

1. 流体力学的数学理论中能量行为的若干初步探讨

Some Preliminary Considerations on Energy Behavior in Fluid Dynamics

作者: Thomas Ruf

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17932v1

摘要翻译: 本研究针对数学流体动力学框架下的能量行为进行了试探性探讨。文中虽就特定条件下能量行为的某些模式提出观察见解，但相关发现的广泛影响尚不明确，需审慎解读。研究结果具有初步性质，其与分析性解的相关性现阶段亟待厘清。本文旨在提出思考以激发后续研究，而非确立任何确定性结论。读者应将 presented material 视为探索性内容，其中存在大量未解开的开放问题。

Comment: 无

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

1. 研究对象和背景

本文研究数学流体动力学中能量行为的某些方面，特别关注Navier-Stokes方程弱解的能量恒等式。核心是引入一个依赖于初始条件 $u_0 \in H^1(\mathbb{R}^3)$ 的局部有界“收紧”权函数 $w$，并考虑一族修正的Navier-Stokes型系统。通过构造一个变形 $r(u) = \rho(|u|)u$（其中 $\rho$ 满足 $\rho'(0)=0$ 及一定的线性增长控制）来替代非线性中的 $u$，得到一个满足特定能量恒等式的近似系统。研究的目标是当修正参数 $\varepsilon \to 0$ 时，将近似系统的解和能量性质传递到原Navier-Stokes系统，并探讨更广泛函数形式（如凸 radially symmetric 函数 $h$）下的能量恒等式。作者强调这些观察是探索性的，旨在激发进一步研究，而非得出最终结论。

2. 主要定理或主要结果

论文的主要结果包括修正系统解的存在唯一性、能量恒等式，以及当 $\varepsilon \to 0$ 时对原系统的逼近，最后给出一个适用于更广类函数的推广能量恒等式。

  • Lemma 2.1：对任意 $\varepsilon > 0$ 及满足条件的 $u_0$ 和 $w$，修正系统 (4) 存在满足 $v_\varepsilon \in H^1_{\text{loc}}(\text{cl}[0,T); H) \cap L^\infty_{\text{loc}}(\text{cl}[0,T); V_\varepsilon)$ 的解。
  • Lemma 2.2：在相同条件下，满足 (9) 的解是唯一的。
  • Lemma 2.3：上述唯一解满足两个能量恒等式 (10) 和 (11)，分别对应动能和某种“耗散”形式。
  • Theorem 2.1：对于任意 $u_0 \in H^1_\sigma(\mathbb{R}^3)$，存在唯一解 $v_\rho$ 满足 (20)，并满足能量恒等式 (21) 和 (22)。该解是通过对修正系统解序列取极限 $\varepsilon \to 0$ 得到的。
  • Lemma 2.6：设 $h$ 为凸的、径向对称且满足 $h(0)=0, \nabla h(0)=0, \|\nabla^2 h\|_\infty < \infty$ 的函数。若初始能量 $\int_{\mathbb{R}^3} h(u_0) < \infty$，则解 $v_\rho$ 满足更一般的能量恒等式 (29)。特别地，当 $h(u) = \frac{1}{r}|u|^r$ ($r \in [2,5]$) 时，(29) 成立。

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

论文的核心方法建立在无限维动力系统的变分框架上，关键技术包括：

  1. Gelfand 三元组理论：将问题置于 $(V, H, V^*)$ 的框架下（$H = L^2(\mathbb{R}^3), V = H^1(\mathbb{R}^3) \cap L^2(\varepsilon w dx)$），利用其推导弱解的存在性。Operator $A = A_1 + A_2$（线性部分 $-\Delta + \varepsilon w$ 和非线性部分 $\operatorname{div}(R(\cdot))$）被证明是伪单调的，这是应用莫纳赫（Montgomery-Smith）型存在定理 [5, Thm. 8.16] 的关键。
  2. Sobolev 链式法则与紧性：变形 $r(u)$ 及其张量积 $R(u) = r(u) \otimes r(u)$ 的连续性（从 $H^1$ 到 $W^{1,1}$ 等）依赖于 Sobolev 链式法则 [1, Thm. 4.4]。紧性由 $V_\varepsilon \hookrightarrow H$ 的紧嵌入保证，这结合了权函数 $w$ 的“收紧”性质（Lemma A.1）与 Rellich-Kondrachov 定理，以及后续证明中使用的 Aubin-Lions 引理。
  3. 能量恒等式推导：利用 Gelfand 三元组框架下测试函数性质，将方程乘以 $v_\varepsilon(t)$ 并积分得到基本能量恒等式 (10)。为得到更精细的 (11) 和 (29)，需要处理非线性项的贡献，技巧性地利用了散度自由条件 $\operatorname{div} v_\varepsilon = 0$ 和向量场的“交替性质”（alternation property，见 [4, Lem. 3.2]），即将 $\langle \operatorname{div}(R(v)), v' \rangle$ 转化为 $\langle (r(v) \cdot \nabla) r(v), v' \rangle$ 并进一步通过分部积分和恒等式 (28)（$\int f(|u|)\langle (v \cdot \nabla) u, u \rangle dx = 0$）来处理。
  4. 压力重建：为得到包含压力项 $p$ 的恒等式 (29)，从方程中消去压力需要求解一个关于 $p$ 的椭圆方程 $-\Delta p = \sum_{i,j} \partial_i \partial_j (r(v_{\rho,i} r(v_{\rho,j}))$。这依赖于 Calderón-Zygmund 奇异积分算子理论，引用 [4, Lem. 5.1] 给出 $p$ 和 $\nabla p$ 的估计，从而保证在所需函数空间中积分是合法的。
  5. 逼近与极限：通过修正参数 $\varepsilon \to 0$，利用能量估计 (23) 和 Diagonal 论证，从修正系统解序列 $v_{\varepsilon_n}$ 中提取弱* 收敛子列，极限函数 $v_\rho$ 满足原系统 (19)。非线性项 $A_2(v_{\varepsilon_n})$ 的收敛性依赖于其线性增长控制 (3) 与 Rellich-Kondrachov 紧嵌入（或 Aubin-Lions 引理）。

4. 与之前工作的比较

论文在引言和正文中并未直接与特定先前工作详细对比，但其方法明显建立在经典文献之上。与标准 Navier-Stokes 存在性理论（如 [4, 5]）相比，本文的创新在于引入了依赖于解初值的局部收紧权函数 $w$ 来构造修正系统，并系统性地推导了在变形 $r(u)$ 下的能量恒等式，特别是推广到一般凸径向函数 $h$。作者坦言“结果的初步性质”及“analytic properties 的关联尚需澄清”，暗示与经典结果（如奇性集、正则性准则）的比较和影响是未完成的开放问题。整体上，本文可视为对标准能量估计技术的一次拓展性应用，旨在探索在非线性项被 $r(u)$ 替换后能量结构的稳健性，为通向原系统提供一条可能路径，但距离与已知深刻结果（如部分正则性理论）建立清晰联系还有很大差距。

2. 定常欧拉-泊松系统亚声速流的整体唯一性

Global Uniqueness of Subsonic Flows for the Steady Euler-Poisson System

作者: Myoungjean Bae, Ben Duan, Chunjing Xie

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17739v1

摘要翻译: 本文证明了有界区域内定常欧拉-泊松系统亚声速解的整体唯一性。既往研究通过将解构造为一维背景态的小扰动，并在一维小扰动框架下应用压缩映射原理，建立了多维亚声速流的存在性与局部唯一性。与之相对，本文摒弃了小性假设，并在满足相同边界数据的一类亚声速解中证明了整体唯一性。

Comment: 无

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

1. 研究对象和背景

本文研究带自洽电场（Case 1）或引力场（Case 2）的等熵稳态Euler-Poisson系统在喷嘴区域 $\mathcal{N} := \{x = (x_1, x') \in \mathbb{R}^n : x_1 \in (0, L), x' \in \Lambda\}$ 中光滑亚声速解的全局唯一性问题，其中 $n=2$ 或 $3$，$\Lambda$ 是有界开连通集。系统在无旋假设（$\nabla \times u = 0$）下可简化为一维常微分方程组（1.3），其解 $(\bar{\rho}, \bar{u}_1 e_1, \bar{\Phi})$ 作为一维背景态。之前的工作（如文献[1-3]）通过将此背景态作为参考，在小边界扰动假设下利用压缩映射原理证明了局部存在性与唯一性。本文旨在此设定下移除小扰动假设，证明在给定相同边界数据的 admissible 函数类中，亚声速解若存在则必唯一，无论其是否靠近背景态。核心挑战在于处理系数 $w(x_1)$ 非常数时带来的拟线性椭圆系统耦合复杂性。

2. 主要定理或主要结果

本文两大核心定理分别针对两种场情形。定理1.2（电场情形）：设 $\mu_0 := \min_{[0,L]} w(x_1) > 0$，且压力律 $p(\rho)$ 满足凸性条件 $\frac{d}{d\rho}\left(\frac{\rho p''(\rho)}{p'(\rho)}\right) \le 0$（如 $p(\rho)=\rho^\gamma, \gamma \ge 1$ 满足）。则存在仅依赖于背景态、$L, \mu_0$ 的小常数 $\sigma_0 > 0$，当边界扰动 $g_0, h_0, v_L$ 的加权Hölder范数之和不超过 $\sigma_0$，且 $h_0$ 满足在端点的导数为零等条件时，问题1.1存在亚声速解 $(\rho, u, \Phi)$，且此解是唯一的。定理1.4（引力情形）：设 $n=2$，$\mu_1 := \max_{[0,L]} w(x_1) < 0$，压力取 $p(\rho)=\rho^\gamma$。则存在小常数 $\sigma_1$，在更强的边界条件（$g_0, h_0, v_L$ 在端点处一阶、二阶导数均为零）及 $\|w'\|_{C^0} \le \sigma_1$ 下，问题1.1存在亚声速解；进一步若 $\gamma \ge 3$，则解唯一。可见，唯一性在电场情形下只需压力凸性，而在引力情形下需更严格的压力指数条件 $\gamma \ge 3$ 和边界光滑性要求。

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

研究方法的核心是将无旋条件 $\nabla \times u = 0$ 通过引入速度势 $\varphi$（$u = \nabla \varphi$）和利用Bernoulli型恒等式，将原始系统（1.5）转化为关于 $(\varphi, \Phi)$ 的拟线性椭圆方程组（2.3）。密度由 $\rho = i^{-1}\left(K_0 + \Phi - \frac{1}{2}|\nabla \varphi|^2\right)$ 显式确定，其中 $i(\rho) = \int_1^\rho \frac{p'(\vartheta)}{\vartheta} d\vartheta$，$K_0$ 为伪Bernoulli常数。关键技术是定义 $\delta$-亚声速集合 $P_\delta := \left\{ (\varphi, \Phi) \in C^2 : p'(\rho) - |\nabla \varphi|^2 \ge \delta \right\}$，并证明其在适当函数空间中是凸集（§2.2）。利用此凸性，可证明若两个解 $(\varphi_1, \Phi_1), (\varphi_2, \Phi_2)$ 均属于 $P_\delta$，则它们的差 $(\psi, \Psi)$ 满足某一线性椭圆问题，通过能量积分与极小化原理可证 $(\psi, \Psi) \equiv 0$，从而得到全局唯一性，此论证不依赖于解是否接近背景态。核心工具包括：加权Hölder空间 $C^{k,\alpha}_{(k,\Gamma)}$ 用于处理边界行为；迭代方法证明存在性；线性化椭圆问题的先验估计（Lemma 2.3），其中关键步骤是推导 $H^1$ 估计时，通过(2.10)分析系数矩阵 $\partial_{q_j} A_i$ 的特征值，利用背景解的亚声速性（(2.5)）得到一致椭圆性条件（2.11），并处理变系数 $w'(x_1)$ 的影响，要求 $\|w'\|_{C^0}$ 足够小以保持估计的可行性；最后通过bootstrap技巧提升正则性。

4. 与之前工作的比较

此前工作（如文献[1-3]）主要依赖构造解作为一维背景态的小扰动，通过压缩映射原理在小的加权Hölder范数邻域内证明局部唯一性。其方法本质是局部化、依赖扰动态。本文的主要突破在于移除了小扰动假设，证明了全局唯一性。方法上，不通过迭代压缩，而是利用 $\delta$-亚声速集的凸性结构，结合能量方法，证明任意两个满足相同边界条件的亚声速解必相同。此方法适用于任意远离背景态的解，只要它们保持亚声速且满足系统。此外，本文考虑了源项系数 $w(x_1)$ 非常数的一般情形，比先前常数系数假设更广泛，但这也技术在估计中带来了额外困难（需 $\|w'\|_{C^0}$ 小），这是论文在§2.1证明存在性、在Lemma 2.5推导线性问题 $H^1$ 估计时详细处理的要点。在唯一性条件上，电场情形仅需压力凸性（1.10），而引力情形（二维）需更强的 $\gamma \ge 3$ 和边界各阶导数为零（(1.12)），这反映了两种情形下系统数学结构（如椭圆型/退化了椭圆型）的差异。

3. 均匀两相流模型的松弛极限：等温情形

The relaxation limit of a homogeneous two-phase flow model: isothermal case

作者: Huimin Yu

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17585v1

摘要翻译: 本文研究了针对均匀两相流设计的双曲松弛系统在松弛时间趋于零时的渐近行为。控制方程组包括混合物质量与动量的守恒律，并辅以一个包含刚性松弛源项的气相质量输运方程，该源项驱使系统趋向局部热力学平衡。在液相密度恒定且气相为理想等温气体的假设下，我们证明了当松弛参数 $ε \rightarrow 0$ 时，解序列 $(p^ε, u^ε)$ 的某个子列在 $L_{\mathrm{loc}}^{1}$ 中强收敛至平衡欧拉系统的一个熵解。证明综合运用了多种分析技术：构造了适当的熵对并推导相关能量估计，采用输运方程方法表征误差项，进行了交换子估计，并应用了补偿紧性理论。该工作为均匀两相流模型的松弛极限提供了严格的数学论证。

Comment: 无

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

1. 研究对象和背景

本文研究均质双相流模型在松弛参数趋于零时的极限行为。均质模型将气体-液体混合物视为均匀混合的“伪流体”，假设各相速度、压力和温度局部相等，适用于充分混合的气泡流或雾流等高压力工况。考虑一维、无粘、无热传导的等温情形，液体密度 $\rho_l$ 为常数，气体相满足理想气体状态方程 $p = \rho_g R T_0$，其中 $p$ 为压力，$R$ 和 $T_0$ 为给定正常数。混合物体积分数 $\alpha$（空占比）满足 $0 \le \alpha \le 1$。控制方程组包括混合质量守恒 $\partial_t \rho_m + \partial_x (\rho_m u) = 0$，混合动量守恒 $\partial_t (\rho_m u) + \partial_x (\rho_m u^2 + p) = 0$，以及带 stiff relaxation 源项的气相质量输运方程 $\partial_t (\rho_g \alpha) + \partial_x (\rho_g \alpha u) = \frac{1}{\varepsilon} (\alpha_{\text{eq}}(p) - \alpha)$，其中 $\rho_m = \alpha \rho_g + (1-\alpha) \rho_l$，$u$ 为共同速度，$\varepsilon > 0$ 为松弛时间，$\alpha_{\text{eq}}(p)$ 是给定压力下的平衡空占比（如由 Clausius-Clapeyron 方程确定）。当 $\varepsilon \to 0$ 时，源项强制 $\alpha$ 趋近 $\alpha_{\text{eq}}(p)$，系统形式上约化为平衡 Euler 方程组：$\partial_t \rho_{\text{eq}}(p) + \partial_x (\rho_{\text{eq}}(p) u) = 0$ 与 $\partial_t (\rho_{\text{eq}}(p) u) + \partial_x (\rho_{\text{eq}}(p) u^2 + p) = 0$，其中 $\rho_{\text{eq}}(p) = \alpha_{\text{eq}}(p) \rho_g(p) + (1-\alpha_{\text{eq}}(p)) \rho_l$。尽管物理上 widely accepted，该松弛极限的严格数学证明此前仍未解决。

2. 主要定理或主要结果

论文主要结果包含三个定理。定理1.1证明了：在合理初始条件下（保守变量 $U = (\rho_m, m, \Gamma)^T$，其中 $m = \rho_m u$，$\Gamma = \rho_g \alpha$，初始数据属于 $L^\infty \cap H^1$，且 $\rho_m$ 有正下界，$0 < \alpha \le 1$），对任意固定 $\varepsilon > 0$，Cauchy 问题存在全局弱解 $U^\varepsilon$，且满足如下 $\varepsilon$ 无关的一致估计：

  • 零阶能量：$\sup_{t} \|U^\varepsilon(t)\|_{L^2} \le C$，误差项积分 $\int_0^T \int_\mathbb{R} (\alpha^\varepsilon - \alpha_{\text{eq}}(p^\varepsilon))^2 \, dx dt \le C \varepsilon$；
  • 一阶空间导数：$\sup_t (\|\partial_x p^\varepsilon(t)\|_{L^2} + \|\partial_x u^\varepsilon(t)\|_{L^2}) \le C$；
  • 时间导数：$\|\partial_t p^\varepsilon\|_{L^2(0,T;L^2)} + \|\partial_t u^\varepsilon\|_{L^2(0,T;L^2)} \le C$；
  • 误差项 $R^\varepsilon$（由 $\alpha^\varepsilon = \alpha_{\text{eq}}(p^\varepsilon) + \varepsilon R^\varepsilon$ 定义）满足 $\|R^\varepsilon\|_{L^\infty(0,T;L^2)} + \sqrt{\varepsilon} \|\partial_x R^\varepsilon\|_{L^2(0,T;L^2)} \le C$。

定理1.2证明：取满足定理1.1估计的解序列，存在子列（仍记作 $\varepsilon$）使得当 $\varepsilon \to 0$ 时，$(p^\varepsilon, u^\varepsilon)$ 在 $L^1_{\text{loc}}((0,T) \times \mathbb{R})$ 中强收敛到一熵解 $(p^0, u^0)$，该解满足平衡 Euler 系统。

定理1.3给出收敛速度：$\|p^\varepsilon - p^0\|_{L^2(0,T;L^2)} + \|u^\varepsilon - u^0\|_{L^2(0,T;L^2)} \le C \varepsilon^{1/2}$。

这些结果共同严格证实了均质双相流模型的松弛极限。

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

证明的核心框架围绕熵方法和紧性论证展开。首先，论文基于亥姆霍兹自由能密度显式构造了一对熵 $(\eta, q)$，其中熵密度为

$$ \eta = -\left[ \frac{\Gamma^2}{2\rho_m} R \ln\frac{\rho_g}{\rho_l} + \frac{\Gamma R}{\rho_m} \ln\frac{\rho_m - \rho_l + \Gamma}{\rho_m} + (\rho_m - \Gamma) \frac{A_l}{\rho_m} \right], $$

熵流 $q = \eta u / \rho_m$。该熵的严格凸性（由 Hessian 矩阵正定保证）是能量估计的基石，且直接关联到 subcharacteristic 条件 $a_f^2 > a_e^2$（冻结声速大于平衡声速），后者是 relaxation 极限可行的关键稳定性条件。其次，为获得独立于 $\varepsilon$ 的一致估计，利用熵耗散不等式 $\partial_t \eta + \partial_x q \le -\frac{C_0}{\varepsilon} (\alpha - \alpha_{\text{eq}})^2$，结合熵的凸性及初始能量有界性，得到零阶能量估计及误差项积分 $\int (\alpha - \alpha_{\text{eq}})^2 dx dt = O(\varepsilon)$。进而，通过构造一阶熵密度 $\eta^1 = \frac{1}{2} (\partial_x U)^T H(U) (\partial_x U)$（$H = D^2 \eta$），并对空间导数系统做能量估计，结合 source 项的特殊结构（ commutator 估计），得到一阶空间导数 $\partial_x p, \partial_x u$ 的 $L^2$ 一致界。时间导数的估计稍复杂：由动量方程可表出 $\partial_t u$，利用 $\|u\|_{L^\infty}$ 和 $\|\partial_x p\|_{L^2}$ 的有界性得 $\|\partial_t u\|_{L^2}$ 有界；但对 $\partial_t p$，直接从连续方程推导会遇到 $\frac{1}{\varepsilon}(\alpha - \alpha_{\text{eq}})$ 项导致可能 $O(1/\sqrt{\varepsilon})$ 的增长，因此需更精细处理。论文发展了一种基于特征线积分的误差表示技术：设 $\alpha^\varepsilon = \alpha_{\text{eq}}(p^\varepsilon) + \varepsilon R^\varepsilon$，将误差项用 $R^\varepsilon$ 表示，并推导 $R^\varepsilon$ 满足的方程，该方程含有导数项，但通过沿特征线积分并利用交换子估计，可以控制其高阶范数，从而避免直接估计 $\partial_t p$ 时的坏项。论文将满足上述一致估计的解序列称为“可接受解族”。第三，证明收敛性（定理1.2）时，先对固定 $\varepsilon > 0$ 引入人工粘性 $\nu > 0$ 正则化系统，利用正则化系统的抛物性质得到光滑解，并通过上述熵方法（适当调整）获得独立于 $\nu$ 的一致估计，然后当 $\nu \to 0$ 时通过 Aubin-Lions 紧性引理得到原松弛系统的弱解，且该解 inherits 一致估计。接着，利用一致估计（特别是 $\partial_x U$ 在 $L^2$ 中有界）和 compensated compactness 理论：由于熵对凸，对均值序列应用 div-curl 引理，可证极限函数满足平衡 Euler 方程的弱形式，且为熵解。最后，定理1.3的收敛率证明采用相对熵方法。定义相对熵泛函并计算其时间导数，利用熵耗散和已知估计，直接证明 $\|(p^\varepsilon, u^\varepsilon) - (p^0, u^0)\|_{L^2}^2 = O(\varepsilon)$，故 $L^2$ 误差为 $O(\varepsilon^{1/2})$。该证明不依赖于极限解的光滑性。

4. 与之前工作的比较

本文工作在前人基础上对一般双相流模型的松弛极限问题做出了重要推进。首先，研究对象的复杂度高：考虑的是变密度（variable-density）的 $3 \times 3$ 双相流系统，而早期经典结果如 Chen-Levermore-Liu (1994) 仅针对 $2 \times 2$ 单相流 p-系统。Jin-Xin 模型虽是 relaxation 典范，但结构不同。本文系统更贴近实际双相流物理，且需处理混合密度与空占比的耦合。其次，在方法上，本文整合了多种分析工具：通过显式熵对构造将 stability condition 与热力学第二定律明确联系；发展了一种基于特征线积分的误差表示技术，巧妙地将 stiff 源项转化为可控制的积分扰动，结合交换子估计避免了直接处理高阶导数（如 $\partial_t p$）的困难；将人工粘性、能量估计与 compensated compactness 紧密结合，在仅 $H^1$ 正则性的弱解框架下完成了收敛证明；并运用相对熵方法获得了最优收敛速率 $O(\varepsilon^{1/2})$。特别地，误差分解技巧是克服 stiff 项导致潜在 $O(1/\sqrt{\varepsilon})$ 增长的关键创新。此外，与一些依赖极限解额外光滑性（如 $\partial_x u^0$ 有界）才能得到收敛速率的工作不同，本文的速率证明完全基于已知均匀估计，不要求极限解高于 $L^\infty$ 的光滑性。文中引用的其他工作如 Berthelin & Bouchut（BGK 模型）、Crin-Barat 等（Baer-Nunziato 到 Kapila）、Solem 等（线性色散分析）等，分别从不同角度考察双相流或 relaxation 极限，但本文是首次对等温均质双相流模型给出从 hyperbolic relaxation 系统到平衡 Euler 系统的完整极限理论，包括存在性、收敛性和收敛速率。

4. 欧拉方程中涡旋爆破率的下界

Lower bounds on the blowup rate of vorticity in the Euler equations

作者: Benjamin Ingimarson, Igor Kukavica

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17431v1

摘要翻译: 在假设三维不可压缩欧拉方程的一个解在时刻 $T_\ast$ 发生爆破，且 $T_\ast$ 为其首次爆破时刻的条件下，我们建立了涡量 $\|\omega\|_{L^\infty}$ 范数爆破率的下界。特别地，当区域为 $\mathbb{R}^3$ 或 $\mathbb{T}^3$ 时，我们给出了当 $t$ 充分接近 $T_\ast$ 时，$\|\omega\|_{L^\infty}$ 在区间 $[0,t]$ 上的累积以及其在 $[0,t]$ 上的上确界的下界。值得注意的是，这定量描述了 BKM 爆破准则。进一步，我们还提供了 $\|D^k \omega\|_{L^\infty}$ 在 $[0,t]$ 上确界的下界。当区域为 $\mathbb{T}^3$ 时，我们建立了当 $t$ 充分接近 $T_\ast$ 时 $\|D^k \omega\|_{L^\infty}$ 的逐时下界。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩Euler方程在整体空间$\mathbb{R}^3$或三维环面$\mathbb{T}^3$上光滑解在有限时间爆破时，涡度（vorticity）$\omega = \nabla \times u$的最大范数$\|\omega\|_{L^\infty}$的爆破速率。研究建立在经典的Beale-Kato-Majda (BKM) 爆破准则基础上，该准则指出解在$T_*$时刻离开光滑类当且仅当$\int_0^{T_*} \|\omega(s)\|_{L^\infty} \, ds = \infty$。然而，BKM准则仅给出了爆破的定性条件，并未提供涡度$\|\omega\|_{L^\infty}$或其积分在爆破时刻$T_*$附近的具体增长速率。在此背景下，先前工作如Chen-Pavlović研究了速度场$u$在$H^r$范数下的爆破下界，而Constantin-Ignatova-Vicol (CIV) 对$\|D\omega\|_{L^\infty}$给出了形如$\sup_{t<T_*} (T_*-t)^{1+\gamma} \|D\omega\|_{L^\infty} = \infty$（对所有$\gamma < 2/5$成立）的下界。本文旨在为$\|\omega\|_{L^\infty}$及其高阶导数$\|D^k\omega\|_{L^\infty}$提供更精确的定量下界，并揭示在$\mathbb{T}^3$上$\|D\omega\|_{L^\infty}$的爆破表现为真正的点态趋于无穷。

2. 主要定理或主要结果

本文主要结果包括四个定理：

  • 定理2.1：在$\mathbb{R}^3$或$\mathbb{T}^3$上，若解在$T_*$首次爆破，则涡度的积累$\int_0^t \|\omega(s)\|_{L^\infty} \, ds$和上确界$\sup_{s \le t} \|\omega(s)\|_{L^\infty}$在$t$充分接近$T_*$时满足显式下界。具体地，存在常数$C>0$使得当$T_*-t < 1/C$时，有$\int_0^t \|\omega\|_{L^\infty} \, ds \ge \frac{1}{C} \log\log \frac{1}{T_*-t}$以及$\limsup_{t\to T_*-} (T_*-t) \log \|\omega(t)\|_{L^\infty} \ge 1/C$。该结果首次为BKM准则中的涡度积分提供了增长的下界估计。
  • 定理2.2：在$\mathbb{R}^3$上，对任意正整数$m$及$k=1,\dots,m$，若初始正则性足够高，则$\sup_{s \le t} \|D^k\omega(s)\|_{L^\infty}$在$t$接近$T_*$时满足$\ge \frac{1}{C (T_*-t)^{\frac{2k}{5}+1} \log^{\frac{2k}{5}+1} \left(\frac{1}{T_*-t}\right)}$，其中临界指数为$\gamma=2/5$（带对数修正）。
  • 定理2.3：在$\mathbb{T}^3$上，若初始正则性足够高，则对接近$T_*$的所有$t$，有$\|D\omega(t)\|_{L^\infty} \ge \frac{1}{C (T_*-t)^{\frac{7}{5}} \log^{\frac{7}{5}} \left(\frac{1}{T_*-t}\right)}$。这表明$\|D\omega\|_{L^\infty}$的爆破是点态的（pointwise-in-time），而非仅上确界形式。
  • 定理2.4：在$\mathbb{T}^3$上，对应力张量$\sigma = \nabla u + (\nabla u)^T$的导数$\|D\sigma\|_{L^\infty}$也有类似定理2.3的点态下界。

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

证明的核心技术在于将含指数项的微分不等式转化为更易处理的形式。对于定理2.1，作者从涡度方程$\partial_t \omega + u\cdot\nabla\omega = \omega\cdot\nabla u$出发，沿拉格朗日轨迹得到$\frac{d}{dt} \|\omega\|_{L^\infty} \lesssim \|\nabla u\|_{L^\infty} \|\omega\|_{L^\infty}$。利用Kozono-Taniuchi估计$\|\nabla u\|_{L^\infty} \lesssim (1+\|\omega\|_{L^\infty})\left(1+\log^+ \|u\|_{H^3}\right)$及能量守恒，推导出形如$\dot{x} \lesssim x + x^2 \exp\left(\int_0^t x(s)\,ds\right)$的微分不等式（其中$x(t)=\|\omega(t)\_{L^\infty}$）。直接处理此不等式因含积分项而困难。关键步骤是时间重新参数化：定义累积涡度$\tau(t) = \int_0^t x(s)\,ds$，由于$x>0$，$\tau$是$[0,T_*)$到$[0,\infty)$的$C^1$同胚。令$X(\tau)=x(t(\tau))$，原不等式化为$\frac{dX}{d\tau} \lesssim 1 + X e^\tau$，由此通过Gronwall不等式得到$X(\tau)$的上界，再反推$T_*-t$与$\tau$的关系，最终导出$\tau(t) \ge \frac{1}{C} \log\log \frac{1}{T_*-t}$，即定理2.1中积分下界。对点态上确界下界（(2.2)式），通过反证法结合积分下界完成。 对于高阶导数（定理2.2），关键在于插值不等式（引理4.1）：$\|\omega\|_{L^\infty} \lesssim \|u\|_{L^2}^{1-\alpha} \|D^k\omega\|_{L^\infty}^{\alpha}$，其中$\alpha=5/(2k+5)$。结合能量守恒$\|u\|_{L^2}$有界及定理2.1中$\|\omega\|_{L^\infty}$的下界，即得$\|D^k\omega\|_{L^\infty}$的下界。 对于$\mathbb{T}^3$上的点态下界（定理2.3），核心是利用环面紧性进行$L^p$估计。由于$\mathbb{T}^3$紧致，$\|D\omega\|_{L^\infty} \gtrsim \|D\omega\|_{L^p}$。通过对$D\omega$的$L^p$范数建立微分不等式$\frac{d}{dt} \|D\omega\|_{L^p} \lesssim p (1+\|D\omega\|_{L^p})^{1+\alpha}$（其中$\alpha=5/(7-6/p)$），积分后得$\|D\omega\|_{L^p}$在$T_*$附近的下界。再通过优化参数$p$（取$p \sim \log \frac{1}{T_*-t}$）并利用$\|D\omega\|_{L^\infty} \gtrsim \|D\omega\|_{L^p}$，最终得到形如$\|D\omega(t)\|_{L^\infty} \gtrsim \frac{1}{(T_*-t)^{7/5} \log^{7/5} \left(\frac{1}{T_*-t}\right)}$的点态下界。证明中使用了环面上的Biot-Savart律估计以及$L^p$意义下导数与涡度关系的线性常数（独立于$p$）性质。

4. 与之前工作的比较

本文工作与之前文献相比有以下几点推进：

  • 相比经典的BKM准则仅知$\int_0^{T_*} \|\omega\|_{L^\infty} = \infty$，本文首次为涡度积累$\int_0^t \|\omega\|_{L^\infty} \, ds$给出了显式的、对数的下界$\frac{1}{C} \log\log \frac{1}{T_*-t}$（定理2.1），从而定量刻画了BKM准则中积分的发散速率。
  • 相比Constantin-Ignatova-Vicol (CIV) 对$\|D\omega\|_{L^\infty}$仅证明对所有$\gamma < 2/5$有$\sup (T_*-t)^{1+\gamma} \|D\omega\|_{L^\infty} = \infty$，本文在$\mathbb{R}^3$上（定理2.2）和$\mathbb{T}^3$上（定理2.3）均达到了临界指数$\gamma=2/5$（定理2.2中为$\frac{2k}{5}$，定理2.3中为$\frac{7}{5}$），尽管带有对数修正。这填补了CIV结果中的指数gap。
  • 本文首次在$\mathbb{T}^3$上证明了$\|D\omega(t)\|_{L^\infty}$的点态下界（定理2.3），表明涡度梯度的爆破是真正在每一点$t$（接近$T_*$）发生，而不仅仅是上确界发散。这比仅对上确界$\sup_{s\le t} \|D\omega\|_{L^\infty}$给出下界的定理2.2更强。
  • 方法上，时间重新参数化技巧（通过累积涡度$\tau$）是处理含指数积分项微分不等式的关键创新，避免了传统屏障方法的困难。对于环面上的点态估计，通过优化$L^p$范数中的$p$，结合紧性得到的$L^p$-$L^\infty$控制，是获得精确下界的核心。
  • 此外，结果推广到了任意高阶导数$\|D^k\omega\|_{L^\infty}$（定理2.2及备注2.4）和应力张量$\sigma$（定理2.4），显示了方法的普适性。

5. 具有摩擦阻尼的随机可压缩Euler方程：$L^\infty$ 鞅解的存在性及渐近类多孔介质行为

Stochastic Compressible Euler Equations with Frictional Damping: Existence of $L^\infty$ Martingale Solutions and Asymptotic Porous Medium-Like Behavior

作者: Rongyi Dai, Jeffrey Kuan, Krutika Tawri, Sunčica Čanić, Konstantina Trivisa

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17348v1

摘要翻译: 本文研究带线性（摩擦）阻尼的一维等熵可压缩Euler方程，受乘性时间白噪声外力驱动。该问题在有限区间上建立，具有$L^\infty$初值和施加于动量的Dirichlet边界条件。我们证明了满足适当熵不等式的全局时间存在性$L^\infty$鞅解。随后分析这些解的长时间行为，并在噪声的合适假设下，证明它们几乎必然且指数快收敛至系统的常数稳态。极限密度可被确定性多孔介质方程的渐近解很好地逼近，而动量则展现出达西定律所预测的渐近行为。

Comment: 无

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

1. 研究对象和背景

本文研究一维等熵可压缩欧拉方程组，该方程组刻画了可压缩流体在受到线性（摩擦）阻尼作用下的运动，同时受到乘性、时间白噪声的随机扰动。系统定义在有界区间 $[0,1]$ 上，状态变量为密度 $\rho$ 和动量 $m = \rho u$（其中 $u$ 为流速），压力由多方气体定律给出 $p(\rho) = \kappa \rho^\gamma$，$\gamma > 1$ 为绝热指数，阻尼系数 $\alpha > 0$。噪声项通过系数 $\sigma(x, \rho, m)$ 乘性作用于动量方程，且 $\sigma$ 假设关于状态变量 Lipschitz 连续，并在 $m=0$ 或 $\rho=0$ 时为零。系统配备 $L^\infty$ 初始数据，并在边界 $x=0,1$ 处施加动量的 Dirichlet 边界条件 $m(t,0)=m(t,1)=0$。核心目标是证明该随机双曲系统的全局时间 $L^\infty$ 鞅弱熵解的存在性，并深入分析其长时间渐近行为，具体而言是证明解几乎必然指数收敛至一个常数稳态，且该稳态的密度渐近由确定性多孔介质方程的解描述，动量则满足 Darcy 定律的渐近形式。该问题的随机版本分析极具挑战性，因为时间上的白噪声扰动可能显著改变系统的长时间统计特性，需要发展新的技术来控制随机效应并证明其衰减。

2. 主要定理或主要结果

本文的主要成果分为两部分：全局解的存在性和长时间渐近行为。

存在性定理（定理 1.1）：假设初始数据 $U_0 = (\rho_0, m_0)$ 满足 $0 \leq \rho_0(x) \leq M_1$，$|m_0(x)| \leq M_2 \rho_0(x)$（$M_1, M_2 > 0$），且噪声系数 $\sigma(x,U)$ 满足：(i) 紧支集包含于 $[0, M_1] \times [-M_2 M_1, M_2 M_1]$；(ii) $\sigma(x, \rho, 0)=0$ 对所有 $\rho>0$，$\sigma(x, 0, m)=0$ 对所有 $m \in \mathbb{R}$；(iii) 关于 $(\rho, m)$ Lipschitz 连续，常数为 $A_0 > 0$；(iv) 边界处 $\sigma(0,\rho,m)=\sigma(1,\rho,m)=0$。则存在全局时间 $[0,\infty)$ 上的鞅 $L^\infty$ 弱熵解（在定义 1.1 意义下）。

长时间行为定理：

  • 定理 1.2：在定理 1.1 的相同假设下，若噪声 Lipschitz 常数 $A_0 < \tilde{C} \min(\alpha,1)$（$\tilde{C}$ 为某正常数），则存在概率 1 的子集 $\bar{\Omega}_3 \subset \bar{\Omega}$，使得对任意 $\omega \in \bar{\Omega}_3$，存在常数 $C_0(\omega, \alpha, \gamma, A_0, \|U_0\|_{L^\infty}) > 0$，使得
    $$ \int_0^1 (\rho(t) - \rho_*)^2 + m^2(t) \, dx \leq C_0 e^{-r t}, \quad \forall t \geq 0. $$
    其中 $\rho_* = \int_0^1 \rho_0(x) \, dx$ 为初始密度平均值，$r = r(\alpha, \gamma, A_0, M_1) > 0$。这表明解几乎必然指数收敛至常数稳态 $(\rho_*, 0)$。
  • 定理 1.3：在相同假设下，设 $(\rho_{IE}, m_{IE})$ 为随机系统的鞅解，$(\rho_{PME}, m_{PME})$ 为具有相同初始密度平均的确定性多孔介质问题
    $$ \begin{cases} \partial_t \rho = \partial_x^2 p(\rho), & p(\rho) = c \rho^\gamma, \\ m = -\frac{1}{\alpha} p(\rho)_x, \end{cases} $$
    配备 $\rho(0,x) = \tilde{\rho}_0(x)$ 满足 $\int_0^1 \tilde{\rho}_0 dx = \int_0^1 \rho_0 dx$，及 $m(t,0)=m(t,1)=0$。则存在几乎必然的指数衰减估计：
    $$ \int_0^1 (\rho_{IE}(t) - \rho_{PME}(t))^2 + (m_{IE}(t) - m_{PME}(t))^2 \, dx \leq C_0 e^{-r t}. $$
    这证明随机解的渐近态被确定性多孔介质方程的解良好近似。

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

研究核心方法为双层近似方案结合熵估计与随机分析工具。

存在性证明：

  1. ε-正则化：对固定大时间区间 $[0,T]$，向原双曲系统添加人工粘性项 $\epsilon \partial_x^2 U$，得到抛物逼近系统，其解具有更高正则性。
  2. τ-时间分裂：对每个 $\epsilon > 0$，将 $[0,T]$ 等分为步长 $\tau$ 的子区间。采用 Lie-Trotter 型分裂格式：在每个子区间 $(t_n, t_{n+1})$ 上，构造连续随机过程，它插值于由 $t_n$ 处解作为初值求解的确定性子问题与由 $t_{n+1}$ 处解作为初值求解的随机子问题之间。这是一种连续变体的 Lie-Trotter 格式，不同于 Berthelin-Vovelle [7] 的分裂。
  3. 紧性与极限 $\tau \to 0$：证明近似解序列的测度律的紧性（tightness），通过 Arzelà-Ascoli 型论证，得到 $\epsilon$-正则化系统在 $[0,T]$ 上的解。
  4. 全局延拓：利用 Gyöngy-Krylov 定理证明 $\epsilon$-正则化系统的路径唯一性，结合标准粘贴论证，将 $[0,T]$ 上的解延拓至全局 $[0,\infty)$。
  5. 极限 $\epsilon \to 0$：这是关键且精细的一步。使用随机补偿紧性理论（基于熵-熵通量对）、Young 测度以及确定性可压缩欧拉方程的经典紧性方法（如 DiPerna-Lions），证明当 $\epsilon \to 0$ 时，$\epsilon$-解收敛至一个鞅 $L^\infty$ 弱熵解。噪声的紧支集假设在此步确保极限解保持 $L^\infty$ 界。

长时间行为证明：

  1. 熵耗散方法：将确定性理论中（如 Huang-Pan [42,43]）用于 damped Euler 方程的熵耗散技术扩展到 $\epsilon$-正则化抛物系统，并适配随机项。
  2. 高阶矩估计：在 $\epsilon$-近似水平，推导熵（或能量）的 sharp 矩估计。具体地，证明 $L^2$ 范数下密度相对其空间平均的偏差和动量的衰减：即 $\mathbb{E} \|\rho_\epsilon(t) - \int_0^1 \rho_\epsilon^0 dx\|_{L^2}^2 \to 0$，$\mathbb{E} \|m_\epsilon(t)\|_{L^2}^2 \to 0$ 当 $t \to \infty$。由于 $\epsilon$-系统满足熵等式（而非极限系统的不等式），可应用 Itô 公式并递推得到高阶矩衰减，最终导出这些量的指数衰减期望估计。
  3. 几乎必然收敛：期望衰减结合 Markov 不等式与 Borel-Cantelli 引理，可推出 $\|\rho_\epsilon(t_n) - \int \rho_\epsilon^0\|_{L^2} \to 0$ 和 $\|m_\epsilon(t_n)\|_{L^2} \to 0$ 几乎必然沿某子列 $\{t_n\}$。但需控制所有路径而非仅子列。期望指数衰减为此提供了基础。
  4. 随机积分处理：利用 Yuskovych [58] 关于随机积分在满足特定矩条件时的几乎必然收敛结果，证明当 $\epsilon \to 0$ 时，$\epsilon$-近似系统中随机项的贡献消失。这允许将 $\epsilon$-水平的衰减结论传递到极限鞅解上，得到定理 1.2 和 1.3 的几乎必然指数衰减。同时，多孔介质方程的近似（定理 1.3）通过比较熵耗散率或直接估计解距获得。

关键工具：凸熵-熵通量对（由系统结构自然导出）、随机补偿紧性（Feng-Nualart [32] 框架）、人工粘性、时间分裂 schemes（Lie-Trotter 连续变体）、Gyöngy-Krylov 路径唯一性定理、Itô 公式、高阶矩估计技术、Yuskovych 随机积分收敛引理。确定性熵耗散方法（如 [42, 52]）是长时间分析的基础，但随机情形需额外精细处理噪声的长期影响。

4. 与之前工作的比较

存在性方面：

  • 确定性理论：经典结果包括 DiPerna (1983) 在 $\mathbb{R}$ 上对 $L^\infty$ 初始数据（特定 $\gamma$）的全局有界熵解；Lions-Perthame-Souganidis (1996) 对所有 $\gamma>1$ 使用动能公式；Ding-Chen-Luo 对一般异质系统。这些工作依赖 vanishing viscosity 或动能公式加补偿紧性。
  • 随机理论：最接近的是 Berthelin-Vovelle [7]，但他们考虑环面 $\mathbb{T}$ 上的 Cauchy 问题，初始数据 $L^2$，噪声满足 $|\sigma|^2 \leq A \rho(1+u^2+\rho^{2\theta})^{1/2}$，解属于 $C(0,T; H^{-2})$。本文不同之处：(1) 区域为有界区间 $[0,1]$ 并带 Dirichlet 边界条件（而非周期边界）；(2) 目标解空间为 $L^\infty$，这得益于噪声的紧支集假设和初始 $L^\infty$ 界；(3) 噪声结构更特殊（紧支集且在边界为零），从物理上合理（如血流模型）；(4) 采用连续 Lie-Trotter 分裂而非 [7] 的分裂。
  • 其他随机工作：Feng-Nualart [32] 考虑空间时间白噪声；Debussche-Vovelle [26] 用动能公式；Bauzet 等 [4] 处理 Dirichlet 边界下的测度值熵解；Breit-Mensah [16] 得局部强路径解；Breit-Feireisl-Hofmanová [13] 显示非唯一性；Hofmanová 等 [37] 引入耗散测度值鞅解；Chen-Huang-Wang [22] 在 $\mathbb{R}$ 上处理有限相对能初值。本文在随机可压缩 Euler 方程中首次获得 $L^\infty$ 鞅解。

长时间行为方面：

  • 确定性理论：Hsiao-Liu [39,40] 对小光滑初值证明收敛至非线性扩散波；Huang-Pan [42,43] 允许真空，证明熵解强 $L^1$ 收敛至多孔介质 Barenblatt 自相似解，动量符合 Darcy 定律；Pan-Zhao [52] 在有界区域 Dirichlet 边界下，建立全局熵解存在唯一并指数收敛至稳态，关联多孔介质。这些工作核心是熵耗散方法和能量估计。
  • 随机理论：此前研究主要集中在 Navier-Stokes 方程上，关注 stationary 解或 invariant 测度的存在性，而非路径收敛。本文是首个对随机 damped Euler 方程严格证明几乎必然路径收敛至稳态的结果。技术上前所未有：需要处理随机积分长时间行为，通过 $\epsilon$-近似的高阶矩估计获得期望衰减，再用 Yuskovych 结果传递几乎必然收敛；同时证明密度和动量同时指数接近确定性多孔介质解，从而在随机框架下验证了 Darcy 定律的渐近有效性。这标志从确定性到随机系统的实质性跨越，且方法具有推广潜力。

6. 基于有限观测的三维NSE正则性判据的再表述与诠释

Reformulation and Interpretation of the Regularity Criterion for 3D NSE Based on Finitely Many Observations

作者: Abhishek Balakrishna, Animikh Biswas

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17322v1

摘要翻译: 我们重新审视并细化了一个近期提出的、关于周期立方体上三维纳维-斯托克斯方程的可观测正则性判据，该判据仅需在给定时间区间内对流动进行有限次测量。本文处理两种数据模型：(i) 模态观测（有限的低傅里叶模集合）；(ii) 节点观测，即速度场在均匀网格上有限个点处的采样值。关键的改进在于，基于每个网格立方体上固定的四面体细分构建了分段线性插值，这一方法消除了先前使用的磨光步骤，并能仅依据测量数据对插值算子的导数实现显式控制。同时，证明了该判据对于正则性既是必要条件也是充分条件。

Comment: 无

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

1. 研究对象和背景

本文聚焦于三维不可压缩Navier-Stokes方程（3D NSE）在周期性域上的正则性问题。所谓正则性，是指弱解是否在时间区间$[0,T]$上具有$H^1$有界性，从而成为强解并满足点态方程。经典正则性判据（如Prodi-Serrin-Ladyzhenskaya准则、Beale-Kato-Majda准则）均要求解在整体空间范数（如$L^s(0,T;L^q(\Omega))$或$L^1(0,T;L^\infty(\Omega))$）下的控制，这在实际实验或观测中往往难以实现，因为这些范数需要全局空间信息。本文的作者先前工作[3,4]提出了基于有限观测（finely many observations）的可观测正则性判据，即仅依赖于有限个傅里叶模态（modal observations）或有限个空间网格点（nodal observations）上的流速测量值。该工作借助数据同化（data assimilation）中的nudging技术，构造一个基于观测数据的插值函数$Iu$，并利用其$H^1$范数作为判据。然而，先前的构造使用分段常数插值后施加光滑化（mollification），导致$H^1$范数与观测数据之间的显式关系不清晰。本文旨在重构并 sharpen 这一判据，通过直接采用分段线性插值消除光滑化步骤，从而建立观测数据与正则性判据之间的透明联系，并证明该判据既必要又充分。

2. 主要定理或主要结果

本文的主要贡献体现在两个定理中。首先，定义与观测尺度$h$（对应于模态数$N$或网格尺寸）及时间区间$[t_0,T]$相关的量$M_{h,[t_0,T]}$，其显式表达式由公式(5.1)给出：在模态观测下，$M_{h,[t_0,T]}^2$本质上是$\|\nabla P_N u\|_{L^2(\Omega)}^2$在$[t_0,T]$上的ess sup，其中$P_N$是前$N$个Stokes特征函数空间上的正交投影，$h \sim 1/\lambda_1^{1/2}$；在节点观测下，$M_{h,[t_0,T]}^2$是$\sum_{\alpha,i} \|D_{\alpha,i}(u,t)\|_{L^2(\Omega)}^2$的ess sup，其中$D_{\alpha,i}(u)$是由四面体顶点值差定义的粗尺度方向导数矩阵（详见公式(4.4)）。定理5.1给出正则性的充分条件：设$u$是满足$u(0)\in V$的Leray-Hopf弱解，若存在$h>0$使得

$$ \max\left\{ \nu\lambda_1, \frac{c\|\nabla u(0)\|_{L^2(\Omega)}^4}{\nu^3}, \frac{c M_{h,[0,T]}^2}{\nu^2\lambda_1} \right\} \leq \frac{\nu}{c h^2} $$

其中$W_h^2 = (\|f\|_{L^2(\Omega)}^2 + c M_{h,[0,T]}^2)/(\nu^2\lambda_1)$，则$u$在$[0,T]$上正则。定理5.2进一步证明此条件既必要又充分：设$u(0)\in V$且$t_0\in(0,T^*)$（$T^*$由[11, Theorem 7.2]保证$u(t)\in V$对$t\in[0,T^*]$），则$u$在$[0,T]$上正则当且仅当存在$h>0$使得

$$ \max\left\{ \nu\lambda_1, \frac{c\|u(t_0)\|_{H^1(\Omega)}^4}{\nu^3}, \frac{c M_{h,[t_0,T]}^2}{\nu^2\lambda_1} \right\} \leq \frac{\nu}{c h^2}, \quad \tilde{W}_h^2 = \frac{\|f\|_{L^2(\Omega)}^2 + c M_{h,[t_0,T]}^2}{\nu^2\lambda_1}. $$

这等价性表明，正则性可通过在足够精细的观测尺度$h$下观测数据是否满足上述不等式来判定。

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

本文的核心方法是数据同化框架，通过构造辅助解$w$追踪参考解$u$，其中关键设计是插值算子$I$，它将有限观测数据嵌入到$w$的方程中。对于模态观测，$I_h u = P_N u$，即投影到前$N$个Stokes特征函数张成的空间。对于节点观测，作者采用了一种精心设计的几何剖分：将周期立方域$\Omega$划分为边长为$h$的小立方体$Q_\alpha$，每个立方体再固定地划分为五个四面体$T_{\alpha}^1,\dots,T_{\alpha}^5$（如图1所示）。参考解$u$在四面体的顶点（即立方体的顶点）上被采样，然后在整个域上构造分段线性函数$\tilde{I}u$，使得在每个四面体上$\tilde{I}u$是线性的，且由顶点值通过重心坐标表示（公式(3.2)和(3.3)）。最后，$\tilde{I}u$被校正为零平均值，得到$Iu$。这一线性插值的关键优势在于可直接计算其梯度，无需光滑化。通过分析四面体的几何结构，引入边向量$w_j = v_{\alpha}^i - v_{\alpha}^4$及由这些边向量生成的常系数矩阵$M$（独立于$h$），作者在公式(4.3)中推导出对于$x\in T_{\alpha}^i$有$\nabla Iu(x) = D_{\alpha,i}(u) M$，其中$D_{\alpha,i}(u)$是公式(4.4)定义的矩阵，其元素为顶点值差除以边长，可视为沿四面体边的方向差商。由此建立了关键等价关系（公式(4.5)）：

$$ \|M\|^{-2} \sum_{\alpha,i} \|D_{\alpha,i}(u,t)\|_{L^2(\Omega)}^2 \leq \|\nabla Iu(t)\|_{L^2(\Omega)}^2 \leq \|M\|^2 \sum_{\alpha,i} \|D_{\alpha,i}(u,t)\|_{L^2(\Omega)}^2. $$

这表明$\|\nabla Iu\|_{L^2}$的范数可被观测数据（顶点值）通过粗尺度方向导数的$L^2$范数之和直接控制。模态情形下，利用傅里叶系数的Parseval恒等式，同样有$\|\nabla P_N u\|_{L^2}^2 \sim \sum_{\lambda_k \leq N^2} \lambda_k |\hat{u}(k)|^2$。这些显式表达式是本文的核心技术突破，使得正则性判据完全由有限观测数据决定。

4. 与之前工作的比较

本文工作直接延续并显著改进了作者之前的结果[4]。在[4]中，插值函数$Iu$首先被构造为分段常数函数（在每个立方体上取平均值），然后通过光滑化（mollification）得到一个光滑函数，后者用于估计$H^1$范数。这种光滑化过程引入了额外的参数，且使得$H^1$范数与原始观测数据之间的关系变得间接和不透明。本文则彻底摒弃了光滑化步骤，采用固定五面体剖分的分段线性插值，使得$Iu$本身已具有$H^1$正则性，其梯度可直接从顶点观测值计算。这带来三重改进：一是建立了$\nabla Iu$与沿四面体边的方向导数$D_{\alpha,i}(u)$之间的直接等价性（公式(4.5)），提供了清晰的物理解释——插值算子的梯度本质上是流速在粗尺度上的方向变化率，这比光滑化后的抽象量更易于理解；二是显式地给出了$M_{h,[t_0,T]}$的表达式，使得正则性判据完全以观测数据（傅里叶系数或节点值）的形式呈现，无需任何隐藏参数；三是证明了该判据的充要性。在[4]中，仅证明了充分性（若数据满足某条件则正则），而本文定理5.2通过反向论证（若解正则，则$M_{h,[t_0,T]}$当$h\to0$时有界，从而存在足够小的$h$满足不等式）完成了必要性证明，从而确立了观测判据与正则性的完全等价。此外，作者还利用[11, Theorem 7.2]说明只需在$(t_0,T]$上验证条件，进一步减少了对初始全区间观测的依赖。总体而言，本研究将可观测正则性判据从理论可能性转化为一个可操作、可解释且充要的准则，为通过有限实验数据验证Navier-Stokes方程解的正则性提供了坚实基础。

7. 科氏力驱动柱状涡旋的渐近线性稳定性

Asymptotic linear stability of columnar vortices driven by Coriolis force

作者: Shuang Miao, Siqi Ren, Zhifei Zhang

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17268v1

摘要翻译: 本文建立了三维轴对称欧拉方程中一类科氏力驱动柱状涡旋的渐近线性稳定性。该结果是证明此类涡旋非线性渐近稳定性的关键一步。所采用的核心且适用范围广泛的策略是构造一个扭曲傅里叶基，这通过求解与系统线性化算子相关的、依赖于双参数 $(c, ξ)$ 的薛定谔方程来实现。为精确刻画解的渐近行为，我们将 $c-ξ$ 平面划分为不同区域，其划分依据是薛定谔方程在不同参数 regime 下的领先阶次剖面。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压缩欧拉方程在柱坐标下的一类特殊定常解——柱状涡旋的渐近线性稳定性问题，重点考虑了受科里奥利力驱动的非均匀旋转流。柱状涡旋是流体中一种基本且普遍存在的 organized rotational motion 结构，其速度场在轴对称假设下具有形式 $v = U(r) e_\theta$，其中 $e_\theta$ 是方位角方向的单位向量，$U(r)$ 是仅依赖径向坐标 $r$ 的角速度函数，对应的涡量 $\omega = \nabla \times v = W(r) e_z$ 沿轴向 $e_z$ 对齐。这类结构在自然和工程环境中广泛存在，如飞机翼尖涡、大气龙卷风及海洋中尺度涡旋。理解其稳定性机制对于揭示涡旋结构如何失稳、破碎并最终过渡到湍流至关重要。在数学理论方面，三维欧拉方程的全局稳定性研究极具挑战性，与二维情形（全局适定性良好）形成鲜明对比。已有的稳定性工作相对 scarce，而稳定性问题不仅是流体动力学稳定性的核心课题，对理解湍流形成机制也具有根本意义。本文特别关注由科里奥利力 $F = -e_z \times v$ 驱动的涡旋，其物理本质体现在旋转参考系中的地转平衡条件，在柱坐标下表现为径向动量方程中的 $- U(r)^2/r + \partial_r P(r) = U(r)$，这反映了旋转带来的离心力与压力梯度的平衡。研究的动机源于近年来 Guo、Pausader 和 Widmayer 对均匀旋转（$U(r) = r$）在轴对称扰动下渐近稳定性的突破性证明，该工作展示了旋转如何阻止奇性形成。本文旨在将此结果推广到更一般的非均匀旋转情形（$U(r)$ 为一般函数），这引入了显著的分析复杂性。

2. 主要定理或主要结果

本文的核心成果是建立了一类满足特定假设（A1)-(A3)）的非均匀旋转柱状涡旋对于三维轴对称扰动的渐近线性稳定性。更精确地，论文证明了线性化后的 Euler-Coriolis 系统（公式 (1.7)）的解具有如下的均匀时间衰减估计：对于 $t \gtrsim 1$，有

$$ \|\partial_z^n \partial_r^m (u_r, u_\theta, u_z)(t)\|_{L^\infty(\mathbb{R}_+ \times \mathbb{R})} \lesssim t^{-1} \| \omega_{\theta,0} \|_{W_\omega^{2+m+n,1}(\mathbb{R}_+ \times \mathbb{R})} + \| u_{\theta,0} \|_{W_\omega^{3+m+n,1}(\mathbb{R}_+ \times \mathbb{R})}. $$

这里 $W_\omega^{k,1}$ 是加权 Sobolev 空间，其范数定义为 $\|F\|_{W_\omega^{k,1}} = \sum_{0 \le i+j \le k} \| \omega(r) \partial_r^i \partial_z^j F \|_{L^1_{r,z}(rdrdz)}$，其中权函数 $\omega(r)$ 具体形式未在摘要中给出，通常与背景流的衰减行为相关。该衰减率 $t^{-1}$ 被认为是优的（optimal），因为在均匀旋转情形（$U \equiv 1$）下，第二作者通过直接傅里叶变换计算已获得相同衰减率。定理 1.1 的证明假设了背景流通过势函数 $V(r) = 2U(r)\Omega(r)$（其中 $\Omega(r) = 2U(r) + rU'(r)$）满足三个条件：(A1) $V'(r) > 0$（ radially strictly increasing）；(A2) $\mathrm{Ran} V = [V(0), 1) \subset (0,1)$（值域有界）；(A3) $V(r) = 1 - a_0 r^{-3} + O(r^{-4})$ 当 $r \to +\infty$，$a_0 > 0$（在无穷远处趋近于 1，具有特定代数衰减）。一个满足这些假设的例子是 $u(r) = \frac{1}{2}(1 - \frac{1}{3}\langle r\rangle^{-3})$，其中 $\langle r \rangle = r+1$。该结果是证明非线性渐近稳定性的关键一步，其证明策略具有更广泛的适用性。

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

论文采用的核心策略是构造一个与线性化算子相关的扭曲傅里叶基（distorted Fourier basis），其构建过程归结为求解一个关于双参数 $(c, \xi)$ 的薛定谔方程（公式 (1.13)）。具体而言，通过对线性化轴对称 Euler-Coriolis 系统（公式 (1.11)）作关于 $z$ 的傅里叶变换（波数 $k \neq 0$），系统约化为关于径向变量 $r$ 的常微分方程组，其中出现一个自伴算子 $A = -V \Delta_{1,k}^{-1}$（在加权空间 $L^2(\mathbb{R}_+; V(r)^{-1}r dr)$ 中）。该算子具有含于 $[0,1]$ 的连续谱，且谱充满整个区间 $[k^{-2}, 1]$（根据上下文，可能为 $[0, k^{-2}]$ 或类似，需核对原文）。系统的解可 formally 表示为（公式 (1.12)）涉及 $A$ 的三角函数作用在初始数据上。因此，点态衰减估计依赖于 $A$ 的谱性质，尤其是其扭曲傅里叶基 $\phi(r, c, k)$，该基满足薛定谔方程 $A \phi_k = c \phi_k$，等价于

$$ \phi''(r) + r^{-1} \phi'(r) - r^{-2} \phi(r) - k^2 \phi(r) + \frac{k^2 V(r)}{c} \phi(r) = 0, $$

其中 $c \in (0,1)$ 是谱参数，$k$ 是傅里叶频率。为分析 $\phi$ 的全局行为并捕捉解的精确渐近行为，论文将双参数 $(c, k)$ 平面（通过变量替换 $\xi = k\sqrt{c^{-1}-1}$，转化为 $(c, \xi)$ 平面）精细地分解为七个不同的区域 $\mathcal{D}_1 - \mathcal{D}_7$，划分依据是薛定谔方程在不同参数 regime 下的主导阶 Profile。这种分区分析是技术核心，因为 $\phi$ 的行为在参数空间的各个区域由不同的特殊函数主导：在低频率区域（$\mathcal{D}_1$）主导函数是贝塞尔函数（Bessel functions）；在特定高频率和谱参数远离 1 的区域（$\mathcal{D}_2, \mathcal{D}_3$）主导行为涉及线性函数 $\xi r$ 和阿贝尔型积分；在谱参数接近 1 的区域（$\mathcal{D}_4, \mathcal{D}_5$）则出现艾里函数（Airy functions）的振荡行为。论文针对每个区域，分别从 $r=0$ 向外和从 $r=+\infty$ 向内求解薛定谔方程，确定 $\phi$ 在 $r$ 小和 $r$ 大时的渐近形式（如公式 (1.19)、 (1.20)），然后通过匹配这些局部解在它们共同有效的重叠区间 $[r_0(\xi,c), r_\infty(\xi,c)]$ 内（该区间长度可能有限），利用 Wronskian 关系 $W(\xi,c) = W(\phi, f_+)$ 将两侧解粘合起来，从而得到全局定义的扭曲傅里叶基 $\hat{\phi}(r,\xi,c)$。最终，解的表示公式（如公式 (1.14) 和 (1.15)）涉及对参数 $c, k$ 的积分和扭曲傅里叶基，衰减估计的证明依赖于对这些积分核进行精细的 $\xi$-$c$ 一致估计，特别是对振荡积分（如公式 (1.16) 和 (1.17)）的估计，这需要利用每个区域中 $\phi$ 的精确渐近展开，并结合分部积分或 stationary phase 方法。整个分析过程需要跨区域进行一致估计以保证最终衰减的普适性，技术细节布满论文第 3 至 8 节。

4. 与之前工作的比较

本文的工作与已有文献既有联系又有显著区别。首先，在柱状涡旋的线性稳定性方面，Gallay 和 Smets 在其系列工作 [27, 28] 中系统研究了三维一般扰动下柱状涡旋的线性稳定性（linear stability），提供了详尽的历史回顾。他们的结果通常关注谱稳定的完整谱分析，而非本文聚焦的时间渐近衰减率。其次，在渐近稳定性（asymptotic stability）方面，Guo, Pausader 和 Widmayer 近期取得了突破，证明了均匀旋转（即 $U(r) \equiv r$）在轴对称扰动下的渐近稳定性 [36]。该工作采用了不同于扭曲傅里叶变换的直接傅里叶方法，并得到了与本文定理 1.1 在均匀旋转特例下相同的 $t^{-1}$ 衰减率（如 Remark 1.2 所述）。本文的核心目标是将均匀旋转情形推广到非均匀旋转，但这一推广带来了众多复杂的分析难题，正如引言和 1.3 节所指出的。具体而言，在均匀旋转情形下，扭曲傅里叶基可以显式地由贝塞尔函数写出，其分析相对直接；而非均匀旋转情况下，由于势函数 $V(r)$ 不再是常数，薛定谔方程 (1.13) 的系数依赖于 $r$，扭曲傅里叶基 $phi$ 无法显式构造，必须通过分区渐近分析才能获得其全局行为。这迫使论文发展出一套精细的、覆盖 $(c,\xi)$ 平面所有区域的渐近分析技术，利用贝塞尔函数、艾里函数以及 WKB 型近似在不同区域主导。从方法论上看，本文采用的扭曲傅里叶变换方法，其根源可追溯到空间时间共振（space-time resonance）方法，并被广泛用于研究色散不够强时的行波稳定性和散射问题，其发展始于 [45]，综述见 [41]。论文在构建扭曲傅里叶基和证明谱定理时，借鉴了 Gesztesy-Zinchenko [34] 对半直线上自伴薛定谔算子的处理，特别是对于强奇异势（如本文出现的 inverse square 势）需要直接验证 resolvent 在谱上的极限，而非依赖 Herglotz 性质。此外，论文也涉及对 $L^1$ 类势和矩阵算子的扭曲傅里叶理论发展，参考了 [7, 17, 18, 43, 44, 46] 等。因此，本文的主要技术创新在于：首次将扭曲傅里叶方法系统地应用于非均匀旋转柱状涡旋的渐近稳定性证明，并发展出了针对该特定薛定谔算子的、参数空间分区高度精细的渐近分析框架，以处理由于 $V(r)$ 非平凡 $r$-依赖所带来的复杂性，从而将均匀旋转情形下的显式分析拓展到更一般的物理相关模型。

8. 对称梯度p-Laplace系统中应力张量的极大Sobolev正则性

Maximal Sobolev regularity of the stress tensor for the symmetric gradient p-Laplace system

作者: Linus Behn, Andrea Cianchi, Lars Diening, Fa Peng

类别: math.AP

发布日期: 2026-03-17

arXiv ID: 2603.17036v1

摘要翻译: 对称$p$-Laplace算子出现于数学物理的各类模型中，例如具有幂型硬化不可压缩材料与非牛顿流体。本研究建立了对称$p$-Laplace系统解的二阶可微性性质，并将其表述为：对局部平方可积的右端项，非线性应力张量具有极大Sobolev正则性。

Comment: 无

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

1. 研究对象和背景

本论文研究如下形式的对称梯度p-Laplace系统：

$$ -\operatorname{div}(A_p(Eu)) = f \quad \text{在}\ \Omega \subset \mathbb{R}^n, $$

其中$u: \Omega \to \mathbb{R}^n$是未知位移场，$Eu = \frac{1}{2}(\nabla u + (\nabla u)^T)$是其对称梯度，$A_p(\xi) = |\xi|^{p-2}\xi$为应力张量，$1 < p < \infty$，$f: \Omega \to \mathbb{R}^n$为外力项。该系统源于塑性理论中的幂硬化不可压缩材料和非牛顿流体动力学模型，其中$Eu$代表无穷小应变张量，$A_p(Eu)$为偏应力张量。

与经典的完整梯度p-Laplace系统$-\operatorname{div}(A_p(\nabla u)) = f$相比，对称梯度系统的正则性理论远不完善。一个核心难点在于，$|Eu|$不像$|\nabla u|$那样是某个标量椭圆方程的子解，因此无法应用De Giorgi迭代方法直接获得$Eu$的正则性。事实上，对于$C^{1,\alpha}$正则性，仅在二维情形$(n=2)$下被证明，高维情形仍是经典公开问题。本文聚焦于二阶可微性，具体体现在应力张量$A_p(Eu)$属于局部Sobolev空间$W^{1,2}_{\text{loc}}(\Omega)$的最大正则性上。值得注意的是，即使右端项$f$是局部平方可积的，该性质也未必对所有$p>1$成立；例如当$f=0$时，存在反例表明$A_p(Eu) \notin W^{1,2}_{\text{loc}}(\mathbb{R}^n)$对于小的$p$值（如$p \leq \frac{3}{2}$）成立。因此，需要确定$p$的可行范围。

2. 主要定理或主要结果

论文的主要结果以定理1.2呈现，建立了应力张量$A_p(Eu)$属于$W^{1,2}_{\text{loc}}(\Omega)$当且仅当外力项$f$属于$L^2_{\text{loc}}(\Omega)$的等价刻画，并给出双向估计。考虑$p$满足$p_-(n) < p < p_+(n)$，其中$p_-(n)$和$p_+(n)$由下表和公式(1.7)定义：

$n$2345678
$p_+(n)$$\infty$$\infty$$\infty$$\infty$$\infty$$\infty$2.5
$p_-(n)$$\approx 1.612$$\approx 1.667$$\approx 1.691$$\approx 1.710$$\approx 1.726$$\approx 1.739$1.75

对于$n \leq 7$，$p_+(n) = \infty$；对于$n \geq 8$，$p_+(n) = 2 + \sqrt{\frac{1}{n+1-1}}$。$p_-(n)$对于$n \geq 3$为$2 - \sqrt{\frac{1}{n+1+1}}$，对于$n=2$略有不同。

设$f \in L^2_{\text{loc}}(\Omega)$，$u$为局部近似解（approximable solution）。则对任意嵌入紧的球$B_{2R} \subset \Omega$，有

$$ A_p(Eu) \in W^{1,2}(B_R) $$

且存在常数$c = c(p,n) > 0$，使得

$$ R^{-1} \|A_p(Eu)\|_{L^2(B_R)} + \|\nabla(A_p(Eu))\|_{L^2(B_R)} \leq c \left( \|f\|_{L^2(B_{2R})} + R^{-n/2-1} \|A_p(Eu)\|_{L^1(B_{2R})} \right). $$

右边的不等式是平凡的，左边是非平凡的核心估计。

由于在低维$2 \leq n \leq 7$时，$p_-(n) > \frac{3}{2}$（由反例1.1），而$p_+(n)=\infty$，因此对于这些维度，仅剩一个狭小区间$(\frac{3}{2}, p_-(n))$的正则性未解决。对于$n \geq 21$，定理要求的$p \geq \frac{2n}{n+2}$可能低于弱解定义所需的下界，此时需采用近似解概念。近似解通过正则化逼近过程定义，仅具有几何测度意义下的近似可微性，但仍可定义近似对称梯度$Eu$。

该结果首次从应力张量$A_p(Eu)$而非传统的能量应力张量$V_p(Eu)=|\ Eu|^{p-2}Eu$角度刻画正则性。同时，结论可推广到更一般的Carreau型系统$-\operatorname{div}(A^\nu_p(Eu)) = f$，其中$A^\nu_p(\xi) = (\nu + |\xi|^2)^{(p-2)/2} \xi$，$\nu \geq 0$；当$\nu > 0$时，算子在$Eu=0$处非奇异。

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

论文的核心技术围绕点态微分不等式展开，并辅以Orlicz空间框架下的正则化逼近。关键技术包括：

Orlicz空间与Young函数理论：为处理更一般的非线性增长，引入Young函数$B$及其共轭$\widehat{B}$，基于函数$a(t)$（对应$A$的标量部分）。通过Simonellen指数$i_a = \inf \theta_a(t)$和$s_a = \sup \theta_a(t)$（其中$\theta_a(t) = a'(t)t/a(t)$）来控制增长和递减。对$a_p(t)=t^{p-2}$，有$i_a = s_a = p-2$。为近似$a_p$，引入正则化函数$a_{p,\varepsilon}(t) = \frac{(\varepsilon + t^2)^{(p-2)/2} + \varepsilon}{1+\varepsilon(\varepsilon + t^2)^{(p-2)/2}}$，满足$i_{a_{p,\varepsilon}} \geq m_p-2$，$s_{a_{p,\varepsilon}} \leq M_p-2$，其中$m_p=\min\{p,2\}$，$M_p=\max\{p,2\}$。相应的$b_{p,\varepsilon}$和$B_{p,\varepsilon}$满足$B_{p,\varepsilon}(t) \approx t^{p}$（在$t$量级意义上），且当$\varepsilon \to 0$时一致收敛到原函数。

点态微分不等式：这是论文最核心的贡献。引理3.1表明，对于满足一定条件的函数$a$，有

$$ |\nabla(A(Eu))| \approx a(|Eu|) |\nabla Eu| \approx a(|Eu|) |\nabla^2 u| $$

在$Eu \neq 0$处成立，若$ a \in C^1 $且附加条件(2.17)，则全局几乎处处成立。证明利用矩阵分析：将$\nabla(A(w))$表示成涉及$w$和$\nabla w$的项，并通过分析相关矩阵的特征值界（依赖于Simonellen指数）获得等价。该估计在向量值情形需要处理张量结构，比标量情况更复杂。

对于一般Orlicz增长，论文推导出关键的微分不等式（见第3节末尾），涉及$\Delta_{\text{sym}} u = \operatorname{div}(A(Eu))$和$\nabla(A(Eu))$的关系。这些不等式是证明W^{1,2}正则性的基石。

正则化与逼近：由于直接解可能正则性不足，论文通过引入正则化算子$A_{p,\varepsilon}$来构造近似系统

$$ -\operatorname{div}(A_{p,\varepsilon}(Eu^\varepsilon)) = f. $$

先对正则化系统证明$A_{p,\varepsilon}(Eu^\varepsilon) \in W^{1,2}$，再令$\varepsilon \to 0$获得原系统解的正则性。这一步需要仔细处理函数空间的收敛和指数$i_B, s_B$的控制。

近似解框架：在第四节定义Orlicz-Sobolev空间$W^{1,B}(\Omega)$及其局部化，第五节引入“近似解”概念：$u$是近似解，如果存在近似解序列$u_k \in W^{2,2}_{\text{loc}}$满足$-\operatorname{div}(A_{p,\varepsilon}(Eu_k)) \to f$在某种意义下。弱解自动是近似解（定理5.1）。这个框架允许处理$p$较大时（$p \geq \frac{2n}{n+2}$）的低正则性解。

正则性提升的闭环：第六节证明对于具有二次增长（对应$p=2$）的非线性系统，若外力属于$L^2$，则应力张量自动属于$W^{1,2}$。这作为逼近过程中的一个关键步骤，确保在取极限时正则性能传递。

4. 与之前工作的比较

本文工作与两个主要文献脉络相关：

标准p-Laplace系统：对于$-\operatorname{div}(A_p(\nabla u)) = f$，正则性理论成熟。标量情形（$N=1$）由Caffarelli等人发展，建立了$A_p(\nabla u) \in W^{1,2}_{\text{loc}}$当且仅当$f \in L^2_{\text{loc}}$（CM18）。向量值情形（$N \geq 2$）在$p > 2(2 - \sqrt{2}) \approx 1.1715$时成立（BCDM22）。这些结果依赖于$|\nabla u|$作为子解的性质，而对称梯度系统缺乏此性质，因此方法不适用。本文首次针对对称梯度系统建立类似的最大正则性，但$p$的范围受限于$p_-(n) > \frac{3}{2}$，比向量值标准系统的下界更高，体现了耦合的更强困难。

对称梯度系统的先前工作：此前关于对称梯度p-Laplacian的正则性主要关注解$u$本身的Hölder连续性（$C^{1,\alpha}$）仅在$n=2$时已知（DKS14, KMS99）；高维情形开放。另一方向是研究“能量应力张量”$V_p(Eu) = |Eu|^{p-2}Eu$的正则性（如Serkov, Seregin等），但该量不能反映$f$的正则性，无法给出类似$A_p(Eu)$的最大正则性刻画。本文选择$A_p(Eu)$而非$V_p(Eu)$作为正则性载体，与近期标准系统中从$V_p(\nabla u)$到$A_p(\nabla u)$的范式转移一致，能更精确地匹配数据$f$的正则性。

点态不等式方法：在标准梯度情形，点态估计$\nabla(A_p(\nabla u)) \approx a_p(|\nabla u|) \nabla^2 u$可通过仅涉及矩阵的技巧证明。本文中，由于对称梯度涉及张量，证明必须处理更复杂的表达式（见(3.6)），并通过分析$w \otimes w$项与单位矩阵的组合的特征值界来克服，这是技术上的关键创新。

Orlicz方法：论文广泛使用Orlicz空间和Young函数，这允许处理Carreau型变系数和非多项式增长。虽然之前工作也使用Orlicz空间，但本文对Simonellen指数的精细控制（确保$B_{p,\varepsilon}(t) \approx t^p$）以及正则化函数$a_{p,\varepsilon}$的设计具有特色。

近似解概念：类似定义出现在De Giorgi-stype方法中（如DHM97），但本文将其应用于对称梯度系统，并特别针对低正则性情形当$p$较大时可能需要更弱解概念的情况。

总之，本文解决了对称梯度p-Laplace系统应力张量最大Sobolev正则性的一个基本问题，通过发展新的点态微分不等式和Orlicz框架下的逼近技术，在宽泛的$p$范围内建立了精确的正则性等价刻画，推进了该领域的理解。

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

1. 位于类空无穷远$i_+$附近的克尔柯西视界的非线性不稳定性

Non-linear instability of the Kerr Cauchy horizon near $i_+$

作者: Sebastian Gurriaran

类别: gr-qc, math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17911v1

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

摘要翻译

中文标题: 位于类空无穷远$i_+$附近的克尔柯西视界的非线性不稳定性

摘要翻译:
我们考虑爱因斯坦真空方程的解，这些解源于位于动态黑洞内部稍近处的光滑初始超曲面，并渐近趋于亚极端克尔黑洞，同时满足精确的非线性Price定律型估计（我们预期该估计在一般情况下成立）。我们证明，相应的最大全局双曲发展存在一片非平凡的未来类空边界——即柯西视界——自类空无穷远$i_+$发出，该边界呈现某种曲率发散现象，且时空度规在其上具有Lipschitz不可延拓性。因此，在所述Price定律型假设下，我们的结果意味着对克尔时空在类空无穷远附近的强宇宙监督假设成立一种Lipschitz版本。

分析依赖于以下工作：Dafermos与Luk关于克尔柯西视界$C^0$稳定性的证明；Giorgi-Klainerman-Szeftel的非可积形式体系及Klainerman与Szeftel在慢旋转克尔黑洞外部稳定性证明中使用的主时间规范；作者关于亚极端克尔黑洞内部Teukolsky方程的线性化分析；以及Sbierski关于Lipschitz不可延拓性的判据。更精确地说，我们通过将黑洞内部分解至带有适当规范的若干区域，得以证明稳定性估计并深入分析Teukolsky方程的非线性类比形式，进而推导出上述不稳定性结果。

2. 关于麦克斯韦-布洛赫方程单频渐近行为的纯态情形研究

On single-frequency asymptotics for the Maxwell-Bloch equations: pure states

作者: A. I. Komech, E. A. Kopylova

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

发布日期: 2026-03-18

arXiv ID: 2603.17888v1

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

摘要翻译

中文标题: 关于麦克斯韦-布洛赫方程单频渐近行为的纯态情形研究

摘要翻译: 本文研究单模麦克斯韦场与二能级分子耦合的阻尼驱动麦克斯韦-布洛赫方程组，该方程组被用于激光作用的半经典描述。我们的主要成果是在准周期泵浦情形下构造出麦克斯韦场具有单频渐近行为的解。此类渐近性适用于满足相互作用绘景下平均约化方程平衡态的正弦初值解。我们计算了所有正弦平衡态并分析其稳定性。本研究基于U(1)规范对称群下的Hopf约化方法，并运用博戈柳博夫-埃克豪斯-桑切斯-帕伦西亚平均理论推导出渐近结果。

3. 异质种群非局部平流-扩散系统

A heterogeneous nonlocal advection–diffusion system

作者: Joseph McCusker, John Christopher Meyer, Mabel Lizzy Rajendran

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17749v1

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

摘要翻译

中文标题: 异质种群非局部平流-扩散系统

摘要翻译: 对定义在 $\mathbb{R}^d$（$d \in \mathbb{N}$）上的异质种群非局部平流-扩散方程组，我们开展了关于局部与整体适定性的独立性研究。每个卷积核 $K_{ij}$ 描述物种 $i$ 依据物种 $j$ 的分布所进行的非局部平流，并假设其具有各自的正则性 $\nabla K_{ij} \in L^{q_{ij}}(\mathbb{R}^d)$（$1 < q_{ij} < \infty$）。通过半群理论与压缩映射论证，建立了温和解及其正则性的局部适定性。对于归类为正则的核族，利用纳什型不等式建立了整体有界性；而对于适合的不规则核族，则通过初始数据的小性条件获得整体有界性。最后给出一个一维数值算例，用以说明核正则性对解的影响。

4. 狄拉克型线性带交叉周期薛定谔方程的有效模型严格推导

Rigorous derivation of an effective model for periodic Schrödinger equations with linear band crossing of Dirac type

作者: Elena Danesi

类别: math.AP, math-ph

发布日期: 2026-03-18

arXiv ID: 2603.17747v1

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

摘要翻译

中文标题: 狄拉克型线性带交叉周期薛定谔方程的有效模型严格推导

摘要翻译: 本文考虑了一类含时一维三次非线性薛定谔方程（NLS），其势函数具有周期性。通过引入半经典缩放并运用多尺度分析，我们推导出一个有效的非线性狄拉克方程，该方程描述了NLS解在谱上局域于狄拉克点附近的动力学行为。

5. 反应扩散方程中非线性波速选择的最优控制方法

An optimal control approach to nonlinear wave speed selection in reaction-diffusion equations

作者: Rebecca M. Crossley, Carles Falco, Ruth E. Baker

类别: math.AP, math.OC, q-bio.TO

发布日期: 2026-03-18

arXiv ID: 2603.17601v1

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

摘要翻译

中文标题: 反应扩散方程中非线性波速选择的最优控制方法

摘要翻译: 反应扩散方程的行波解被广泛用于模拟生物学和物理学中种群的空间传播及其他现象。本文通过最优控制建模重新诠释了经典的变分原理方法，以获取非线性偏微分方程系统中行波解侵入速度的下界。我们首先分析单物种模型，其中密度演化由含密度依赖扩散项和非线性反应项的标量方程控制。证明对任意容许测试函数，针对感兴趣参数进行最大化可得到波速上界。将该框架应用于多孔Fisher方程等示例，并考察非线性选择机制何时超越经典线性 marginal 稳定性判据。进一步推广至多物种反应扩散方程组，在弱耦合条件下，以庞特里亚金型最优性系统重构后，通过变分框架推导出类似的波速上界。最后，采用数值模拟验证了所预测波速在多个示例中的准确性。

6. 动力 Sobolev 空间

Kinetic Sobolev Spaces

作者: Pascal Auscher, Lukas Niebel

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17491v1

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

摘要翻译

中文标题: 动力 Sobolev 空间

摘要翻译: 我们定义并研究了适用于科尔莫戈罗夫方程的齐次动力 Sobolev 空间。我们同时考虑了局部与非局部扩散情形。这些空间基于勒贝格空间 $ L^p $（对所有可积指数 $ p \in (1, \infty) $）构建，并在输运与扩散方向上依据科尔莫戈罗夫方程的标度施加正则性条件。该正则性标度同时适用于弱解与强解。我们证明了所提出的空间满足以下方面最佳嵌入估计：以 Bouchut-Hörmander 风格表征的正则性传递、 Lions 风格的时间连续性，以及 Sobolev 与 Hardy-Littlewood-Sobolev 类型的可积增强。核心工具是科尔莫戈罗夫算子（由基本解给出）在各向异性齐次 Sobolev 空间之间的映射性质。为此，我们证明了相关奇异积分算子的 $ L^p $ 有界性，并通过李氏-佩利分解与几何分析推导出该算子新的核估计。此外，我们提供了一个新的唯一性准则，用以证明柯西问题的适定性。

7. 局部/非局部系统的周期均质化

Periodic Homogenization of Local/Nonlocal Systems

作者: Marcone C. Pereira, Luiza C. Rosa da Silva, Julio D. Rossi

类别: math.AP

发布日期: 2026-03-17

arXiv ID: 2603.17144v1

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

摘要翻译

中文标题: 局部/非局部系统的周期均质化

摘要翻译: 本文研究一类椭圆方程的均质化问题，该方程由两部分耦合而成：局部部分由带Neumann边界条件的拉普拉斯算子给出，非局部部分则由光滑核的积分算子定义。这两个分量通过另一个同样由光滑核给出的非局部算子相互耦合。我们考虑将一个固定空间域进行一系列划分，划分为局部区域与非局部区域，二者在空间中以周期方式分布（其中一种区域由周期排列的小孔洞构成）。根据局部与非局部区域的相对位置，我们得到定性不同的极限行为。当方程的局部部分局限于这些小周期孔洞时，解序列收敛至一个极限系统的唯一解，其中局部分量消失，而非局部部分保留并分裂为两个独立的分量。另一方面，当问题的局部部分位于孔洞之外时，极限系统呈现一个均质化的局部扩散算子与非局部方程的耦合。最后，我们通过考虑由平行细条带（而非孔洞）构成的配置，分析了中间情形：此时仅部分局部扩散在极限中得以存留。

不相关论文（基本信息）

1. 一类拟线性随机演化方程的强解存在唯一性

Strong existence and uniqueness for a class of quasilinear stochastic evolution equations

作者: Sebastian Bechtel, Esmée Theewis

类别: math.PR, math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17857v1

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

摘要: We establish existence of probabilistically strong solutions and pathwise uniqueness for a class of quasilinear stochastic evolution equations on bounded domains. Our results combine recent weak existence results for quasilinear stochastic evolution equations in an $L^p$-setting (with $p > 2$) with Yamada–Watanabe theory. To establish pathwise uniqueness, we rely on an $L^1$-contraction argument.

2. 闵可夫斯基问题的定量稳定性

Quantitative Stability for Minkowski’s problem

作者: Károly Böröczky, João Miguel Machado, João P. G. Ramos

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

发布日期: 2026-03-18

arXiv ID: 2603.17726v1

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

摘要: We derive quantitative stability results for Minkowski bodies, as well as their counterparts, the $L_p$-Minkowski bodies in the range $1 \le p \neq n$. We prove that, for every pair of probability measures $μ,ν$ satisfying a quantitative form of the classical dispersion assumptions yielding existence of such bodies, we have a control of the form

$$ α(E_μ, E_ν)^2 \le C \mathrm{d}_{\mathrm{C}}(μ,ν)^{1 + \frac{1}{n}}, $$

where $α$ denotes the Fraenkel asymmetry and $\mathrm{d}_{\mathrm{C}}$ is the dual-convex distance of probability measures on the sphere. Our arguments are based on a variational problem whose optimizers are Minkowski bodies, for which we can obtain strong-concavity properties with the quantitative Brunn-Minkowski and isoperimetric inequalities.

3. 当支撑区域测度趋于零时一类混合边值问题主特征值的极限行为

Limiting behavior of principal eigenvalues for a class of mixed boundary value problems as the measure of the support domain goes to zero

作者: J. Lopez-Gomez, A. Sahuquillo

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17721v1

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

摘要: In this paper we characterize the limiting behavior of the principal eigenvalue, $\s_1[-\D,\b,Ø]$, of the boundary value problem \eqref{1.1} as the Lebesgue measure of the underlying domain, $Ø$, tends to zero. Naturally, the domains $Ø$ are assumed to be included on a fixed open set $D$ such that $\b\in\mc{C}(D)$, and they satisfy $\barØ\subset D$. Our main result establishes that, in the classical case when $\inf_{D}\b >0$,

$$ \lim_{|Ø|\da 0}\s_1[-\D,\b,Ø] =+\infty, $$

whereas

$$ \lim_{|Ø|\da 0}\s_1[-\D,\b,Ø] =-\infty\;\;\hbox{if}\;\; \sup_{D}\b <0, $$

which is a surprising result at the light of the classical existing theorems for Dirichlet boundary conditions. Furthermore, in the special case when $\b>0$ is a constant, we can prove that

$$ \lim_{R\da 0}\left( R \s_1[-\D,\b,B_R]\right)=\b \frac{\mathrm{Area}(\p B_1)}{|B_1|}, $$

where, we are denoting $B_\varrho:=\{x\in\R^N\;:\;|x|<\varrho\}$ for all $\varrho>0$.

4. 论阻尼驱动Jaynes-Cummings方程的全局动力学

On global dynamics for damped driven Jaynes-Cummings equations

作者: A. I. Komech, E. A. Kopylova

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

发布日期: 2026-03-18

arXiv ID: 2603.17562v1

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

摘要: The article concerns damped driven Jaynes-Cummings equation which describes quantised one-mode Maxwell field coupled to a two-level molecule. We consider a broad class of damping and pumping which are polynomial in the creation and annihilation operators, and their structures correspond to the theory of completely positive and trace preserving generators (CPTP) of Lindblad and Kossakowski & al. Our main result is the construction of global generalised solutions with values in the Hilbert space of nonnegative Hermitian Hilbert-Schmidt operators in the case of time-dependent pumping. The proofs rely on finite-dimensional approximations of the annihilation and creation operators.

5. 偏微分方程的层次化稳定性概念与李雅普诺夫函数

Hierarchical Stability Notions and Lyapunov Functions for PDEs

作者: Matthew M. Peet

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

发布日期: 2026-03-18

arXiv ID: 2603.17536v1

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

摘要: Unlike linear ordinary differential equations (ODEs), linear partial differential equations (PDEs) admit a multitude of non-equivalent notions of stability. This variety makes interpretation of Lyapunov stability results challenging. To simplify this interpretation, we propose a framework for hierarchical classification of notions of stability and Lyapunov conditions. To do this, for every well-posed PDE and set of boundary conditions, we define a fundamental state on $L_2$ corresponding to the minimal information needed to uniquely forward propagate the solution. Stability notions and Lyapunov functions are then defined in terms of this fundamental state. This gives rise to a hierarchy of stability notions, the weakest being fundamental state to PDE state stability. Other stability notions and Lyapunov conditions may then be interpreted relative to this weakest notion. Hierarchies are established for: Lyapunov, exponential and finite-energy stability. Sufficient Lyapunov conditions are defined in terms of operator inequalities. Illustrative examples and computational tools are provided.

6. 分片光滑介质中具有阻抗边界条件的麦克斯韦方程组的波数显式hp有限元法分析

Wavenumber-explicit $hp$-FEM analysis of Maxwell’s equations with impedance boundary conditions in piecewise smooth media

作者: Jens Markus Melenk, David Wörgötter

类别: math.NA, math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17467v1

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

摘要: We consider the time-harmonic Maxwell equations with impedance boundary conditions on a bounded Lipschitz domain $Ω$ with analytic boundary $Γ$. We suppose that $Ω$ consists of multiple subdomains, and that the permeability and permittivity tensors are analytic on every subdomain, but may jump across subdomain interfaces. Under these conditions we show that for any wavenumber $k\in\mathbb{C}$ with $|k|\geq 1$ for which Maxwell’s equations are polynomially well-posed, a Galerkin discretization based on Nédélec elements of order $p$ on a mesh with mesh width $h$ is quasi-optimal, provided that there holds the wavenumber-explicit scale resolution condition a) that $|k|h/p$ is sufficiently small and b) that $p/\log |k|$ is bounded from below.

7. 奇异扰动梯度流与无穷维空间临界点的演化

Singularly Perturbed Gradient Flows and Evolution of Critical Points in Infinite Dimensions

作者: Virginia Agostiniani, Riccarda Rossi, Giuseppe Savaré

类别: math.AP

发布日期: 2026-03-18

arXiv ID: 2603.17458v1

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

摘要: We consider singularly perturbed gradient flows in Hilbert spaces, driven by a time-dependent, nonconvex, and nonsmooth energy, and address the convergence of their solutions to curves of critical points of the driving energy functional. The degenerating nature of the estimates along the gradient-flow curves calls for novel compactness arguments, which we carefully develop by combining tools from the variational approach to Hilbert and metric gradient flows \cite{RossiSavare06,AGS08}, with fine requirements on the set of critical points of the energy. This leads us to prove that subsequential limits of singularly perturbed gradient flows are Dissipative Viscosity solutions of the limiting problem, i.e., a curves of critical points satisfying a suitable balance between the energy and a defect measure, encoding dissipation. This energy-dissipation balance encompasses information on the dynamics of the process at jump times, recording, in particular, the re-emergence of viscous behavior. Under a suitable rectifiability condition on the critical set, we show that Dissipative Viscosity solutions improve to Balanced Viscosity solutions, which have the key property that the dissipation measure is purely atomic. In the second part of the paper we show that, for smooth energies whose second differential is a Fredholm operator, the condition that the kernel of the Hessian has dimension at most one at every critical point already implies our measure-theoretic assumptions. We further relate them to the transversality conditions from bifurcation theory and show that they have a generic character.

8. 布朗运动在渐近双曲与类木星气体巨星度量曲面上的平均首次逃逸时间

Mean first escape times of Brownian motion on asymptotically hyperbolic and gas giant metric surfaces

作者: Jesse Gell-Redman, Emanuel József Godfried, Justin Tzou, Leo Tzou

类别: math.AP, math.PR

发布日期: 2026-03-18

arXiv ID: 2603.17313v1

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

摘要: This paper deals with the mean first escape time of Brownian motion on asymptotically hyperbolic and gas giant surfaces. We show that for a boundary defining function $ρ$, the mean first escape time $u_ε(x)$ from the truncated Riemannian surface with an asymptotically hyperbolic metric $(M_ε,\bar{g}/ρ^2) = (\{x\in M:ρ(x)\geq ε\},\bar{g}/ρ^2) \subset (M,\bar{g}/ρ^2)$ satisfies the asymptotic expansion $u_ε(x) = -\log ε+ \mathcal{O}(1)$ as $ε\to 0 $. Furthermore, we show that in the case of a gas giant metric $g = \bar{g}/ρ^α$, where $α\in (0,2)$, the mean first escape time from the surface $(M_ε,\bar{g}/ρ^α)$ satisfies $u_ε(x) = \mathcal{O}(1)$ as $ε\to 0 $. Using techniques from the theory of polyhomogeneous conormal functions we explain this difference between in the mean first escape time on gas giant metric surfaces and asymptotically hyperbolic surfaces on the unit disc. Finally, we confirm these results using Monte Carlo simulations and finite difference methods on the disc.