生成时间: 2026年06月23日 06:37:48

运行状态: 完成

计划处理论文数量: 41

已完成论文数量: 41

未完成论文数量: 0

重点关注论文数量: 8

了解领域论文数量: 16

不相关论文数量: 17

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

1. 关于二维广义不可压模型

On a generalized incompressible model in two dimensions

作者: Hantaek Bae, Kyungkeun Kang, James P Kelliher, Woojae Lee

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22891v1

摘要翻译: 我们分析 Ohkitani [19] 提出的一个广义不可压模型。该模型基于如下观察：二维 Burgers 方程可通过将速度梯度旋转 90 度而与不可压 Navier-Stokes 方程相关联。我们给出若干初始数据属于 $H^{3}$ 的结果。首先，我们考察无粘性模型，并证明存在唯一一个局部时间解；若初始涡量含有负部，则该解在有限时间内爆破。在有粘性情形下，我们证明存在唯一全局时间解，且无需对初始涡量施加符号条件；建立两个解之差的长时间行为；并在初始涡量非正时导出速度场的时间衰减率。

Comment: 16 pages

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

1. 研究对象和背景

本文研究 Ohkitani 在 [19] 中提出的二维广义不可压流体模型，其核心方程为

v_t + P(R_{-\alpha}(v) \cdot \nabla v) - \nu \Delta v = 0, \quad \mathrm{div}\, v = 0, \quad 0 < \alpha < \frac{\pi}{2}, \tag{1.5}

其中 $$ P $$ 是 Leray 投影算子， $R_{-\alpha}$ 是旋转矩阵作用于速度场：

R_{-\alpha}(v) = \begin{pmatrix} \cos\alpha & \sin\alpha \\ -\sin\alpha & \cos\alpha \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix}, \quad D_\alpha = R_\alpha \nabla = \begin{pmatrix} (\cos\alpha)\partial_1 - (\sin\alpha)\partial_2 \\ (\sin\alpha)\partial_1 + (\cos\alpha)\partial_2 \end{pmatrix}.

该模型通过连续参数 $\alpha$ 插值两个经典极限：当 $\alpha = 0$ 时退化为标准二维不可压 Navier–Stokes 方程（NS）

u_t + u \cdot \nabla u + \nabla p - \nu \Delta u = 0, \quad \mathrm{div}\, u = 0; \tag{1.1}

当 $\alpha = \frac{\pi}{2}$ 时对应由二维 Burgers 方程经正交梯度变换导出的不可压系统

e_u = \nabla^\perp \psi, \quad (e_u)_t + e_u \cdot \nabla^\perp e_u - \nu \Delta e_u = 0, \quad \mathrm{div}\, e_u = 0. \tag{1.4}

关键物理量为涡量 $\omega = \partial_1 v_2 - \partial_2 v_1$ ，其满足非标准输运-扩散方程

\omega_t + R_{-\alpha}(v) \cdot \nabla \omega - \nu \Delta \omega = - (\sin\alpha) |\nabla v|^2. \tag{2.7}

与 NS 不同，该方程右侧含负定源项 $- (\sin\alpha)|\nabla v|^2$ ；与 Burgers 不同，其对流速度 $R_{-\alpha}(v)$ 不满足 $\mathrm{div}\, R_{-\alpha}(v) = 0$ ，而仅满足 $\mathrm{div}\, R_{-\alpha}(v) = (\sin\alpha) \omega$ —— 这一结构破坏了最大值原理和经典涡量估计。问题的重要性在于：它提供了一个可解析操控的“桥梁”系统，用于探测从标量守恒律（Burgers）到向量不可压流（NS）之间正则性机制的根本转变。已有理论表明：(i) NS 在 $\nu > 0$ 、二维、 $$ H^3 $$ 初值下全局适定；(ii) 无粘 NS（Euler）在二维亦全局适定，依赖于涡量 $L^\infty$ 控制（Beale–Kato–Majda 准则）；(iii) 无粘 Burgers 在二维存在有限时间爆破（因 $\omega$ 不满足输运方程）。然而，对插值模型 (1.5)，两个基本问题悬而未决：(A) 当 $\nu = 0$ 时，是否存在初值导致有限时间爆破？(B) 当 $\nu > 0$ 时，是否无需符号条件（如 $\omega_0 \geq 0$ 或 $\omega_0 \leq 0$ ）即可获得全局解？本文旨在填补这两处关键空白：首次严格证明 $\nu = 0$ 下负涡量初值必然导致爆破，并首次建立 $\nu > 0$ 下对任意 $$ H^3 $$ 初值的全局适定性——这突破了以往类似模型（如 $\alpha$ -Euler 或 SQG 型方程）对初值符号或小初值的依赖。

2. 主要定理或主要结果

本文在初始数据 $v_0 \in H^3(\mathbb{R}^2)$ 、 $\mathrm{div}\, v_0 = 0$ 下建立以下主要结果：

定理 3.1（无粘局部适定性）：存在唯一解 $v \in C([0,T); H^3)$ 满足

\|v(t)\|_{H^3} \leq \frac{\|v_0\|_{H^3}}{1 - C t \|v_0\|_{H^3}}, \quad 0 \leq t < T = \frac{1}{C \|v_0\|_{H^3}},

且最大存在时间 $T^* < \infty$ 当且仅当

\limsup_{t \nearrow T^*} \int_0^t \|\omega(\tau)\|_{L^\infty} \, d\tau = \infty.

定理 3.2（无粘有限时间爆破）：若 $\min_{x \in \mathbb{R}^2} \omega_0(x) < 0$ ，则解不能全局存在，且爆破时间满足上界

T^* \leq -\frac{2}{(\sin\alpha) \min_{x \in \mathbb{R}^2} \omega_0(x)}.

定理 4.1（粘性全局适定性）：对任意 $\nu > 0$ （文中取 $\nu = 1$ ），存在唯一全局解 $v \in C([0,\infty); H^3)$ 满足

\|v(t)\|_{H^3}^2 + \int_0^t \|\nabla v(\tau)\|_{H^3}^2 \, d\tau \leq C \|\omega_0\|_{H^2}^2 \exp\left[ C \|v_0\|_{L^2}^2 \left(1 + \frac{\sin\alpha}{2} \|\omega_0^+\|_{L^\infty} t \right)^2 \right] + \|v_0\|_{L^2}^2 \left(1 + \frac{\sin\alpha}{2} \|\omega_0^+\|_{L^\infty} t \right)^2.

定理 4.2（参数差的长时间行为）：设 $$ v_i $$ 对应角度 $\alpha_i$ 的解（ $$ i = 1,2 $$ ），定义 $\tilde{v} = v_1 - v_2$ ，则

\|\tilde{v}(t)\|_{H^3}^2 + \int_0^t \|\nabla \tilde{v}(\tau)\|_{H^3}^2 \, d\tau \leq \left[ \|\tilde{v}_0\|_{H^3}^2 + C_\alpha t J^{(1)}(t) \right] \exp\left[ C \left( J^{(1)}(t) + J^{(2)}(t) \right) \right],

其中 $C_\alpha = C \max\{ |\cos\alpha_1 - \cos\alpha_2|, |\sin\alpha_1 - \sin\alpha_2| \}$ ，且

J^{(i)}(t) = \|\omega^{(i)}\|_{H^2}^2 \exp\left[ C \|v^{(i)}\|_{L^2}^2 \left(1 + \frac{\sin\alpha_i}{2} \|(\omega^{(i)})^+\|_{L^\infty} t \right)^2 \right] + \|v^{(i)}\|_{L^2}^2 \exp\left[ C \left(1 + \frac{\sin\alpha_i}{2} \|(\omega^{(i)})^+\|_{L^\infty} t \right)^2 \right].

定理 5.1 与 5.2（衰减率）：若 $\omega_0 \leq 0$ 且 $v_0 \in H^3 \cap L^1$ ，则 $\|v(t)\|_{L^2} \leq C (1+t)^{-1/2}$ ；若 $v_0 \in H^3 \cap \dot{H}^{-\sigma}$ （ $\sigma > 0$ ），则 $\|v(t)\|_{L^2} \leq C (1+t)^{-\sigma/2}$ 。

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

证明的核心障碍在于模型 (1.5) 缺乏两类经典结构：(i) 无粘情形下， $R_{-\alpha}(v)$ 的散度非零（ $\mathrm{div}\, R_{-\alpha}(v) = (\sin\alpha)\omega$ ），导致能量估计中出现危险的二次项 $\frac{\sin\alpha}{2} \int \omega |D^k v|^2$ （见 (3.3) 式），无法直接控制；(ii) 粘性情形下，涡量方程 (2.7) 右侧含负定源项 $-(\sin\alpha)|\nabla v|^2$ ，但该源项无法直接用于控制 $\omega^-$ ，且 $\omega^+$ 的演化受 $R_{-\alpha}(v)$ 输运影响，传统最大值原理失效。

作者的第一个创新点是构造并利用一个精巧的比较 ODE 来控制 $\omega^+$ 。具体而言，由 (2.7) 推出不等式 $\omega_t - \Delta \omega + R_{-\alpha}(v) \cdot \nabla \omega + \frac{\sin\alpha}{2} |\omega|^2 \leq 0$ （引理 4.3），再与 ODE $y' + \frac{\sin\alpha}{2} y^2 = 0$ 比较，通过 Stampacchia 型截断（引理 4.1 和 4.2）证明 $(\omega - y(t))^+ \equiv 0$ ，从而得到显式衰减估计 $\|\omega^+(t)\|_{L^\infty} \leq \frac{2\|\omega_0\|_{L^\infty}}{2 + (\sin\alpha)\|\omega_0\|_{L^\infty} t}$ 。这一估计是全局适定性的基石，它将 $\omega^+$ 的增长限制为对数型，而非指数型。

第二个创新点是发展了一套针对 $\alpha$ -依赖非线性项的高阶能量估计框架。在定理 4.1 的证明中，作者并未对 $$ v $$ 直接做 $$ H^3 $$ 估计，而是分层估计涡量 $\omega$ 的 Sobolev 范数：首先用 $\omega^+$ 的衰减控制 $$ L^2 $$ 范数（(4.9)），再利用 $\|\omega\|_{L^4} \lesssim \|\omega\|_{L^2}^{1/2} \|\nabla \omega\|_{L^2}^{1/2}$ （Ladyzhenskaya 不等式）控制 $\|\nabla \omega\|_{L^2}$ （(4.10)），最后对 $\Delta \omega$ 做估计时，关键处理了 $-\int \Delta(R_{-\alpha}(v) \cdot \nabla \omega) \Delta \omega$ 项，通过展开三阶导数并反复应用插值不等式 (2.2) 和乘积法则 (2.3)-(2.4)，将其控制为 $C \|\omega\|_{L^2}^2 \|\Delta \omega\|_{L^2}^2 + \frac{1}{2} \|\nabla \Delta \omega\|_{L^2}^2$ ，从而闭合估计。整个高阶估计的权重函数 $I(t) = \|v_0\|_{L^2}^2 (1 + \frac{\sin\alpha}{2} \|\omega_0^+\|_{L^\infty} t)^2$ 直接源于 $\omega^+$ 的衰减率，体现了方法的内在一致性。

第三个创新点是分析两个不同 $\alpha$ 解之差 $\tilde{v} = v_1 - v_2$ 的稳定性。作者将差方程 (4.13) 中的非线性项精确分解为三部分：(I) 参数差引起的强迫项 $P[(R_{-\alpha_1} - R_{-\alpha_2})(v_1) \cdot \nabla v_1]$ ，其 $$ H^3 $$ 范数被 $C_\alpha \|v_1\|_{H^3}$ 控制；(II) $\tilde{v}$ 与 $$ v_1 $$ 的耦合项，用乘积法则 (2.4) 控制；(III) $$ v_2 $$ 与 $\tilde{v}$ 的耦合项，利用 $\mathrm{div}\, R_{-\alpha_2}(v_2) = (\sin\alpha_2) \omega_2$ 将其转化为含 $\omega_2$ 的二次型，并最终吸收进耗散项。这一分解使得差估计的常数 $C_\alpha$ 显式依赖于 $|\alpha_1 - \alpha_2|$ ，为参数连续性的定量刻画提供了基础。

4. 与之前工作的比较

本文结果与既有工作形成多维度对比：

与 Ohkitani [19] 的对比：Ohkitani 提出了模型 (1.5) 并进行了数值模拟，观察到其行为随 $\alpha$ 连续变化，但未给出任何解析适定性或爆破结果。本文首次为其提供了严格的数学基础，将数值现象上升为定理：定理 3.2 证实了其数值观察到的“负涡量导致奇异性”的猜想；定理 4.2 则从解析层面验证了其“参数连续性”主张，给出了差解的显式 $\alpha$ -依赖估计。

与经典 NS 和 Burgers 的对比：在无粘情形，NS（Euler）在二维全局适定，而 Burgers 在二维爆破；本文揭示了一个精细的相变：当 $\alpha \in (0, \frac{\pi}{2})$ ，模型既非完全 NS 型（因缺乏 $\mathrm{div}\, R_{-\alpha}(v) = 0$ ），亦非完全 Burgers 型（因仍具不可压约束），其爆破判据 $\min \omega_0 < 0$ 介于两者之间——NS 无此条件，Burgers 的爆破条件更弱（仅需初值斜率负）。在粘性情形，NS 的全局适定性依赖于二维涡量守恒 $\|\omega(t)\|_{L^\infty} = \|\omega_0\|_{L^\infty}$ ，而本文的 $\|\omega^+(t)\|_{L^\infty}$ 衰减（引理 4.3）是全新的机制，它不依赖于 $\omega$ 的整体守恒，而仅利用源项 $-(\sin\alpha)|\nabla v|^2$ 的负定性，这是对经典理论的重要补充。

与类似广义模型的对比：与 $\alpha$ -Euler 方程（ $v_t + P((\nabla^\perp)^{1-\alpha} \omega) \cdot \nabla v = 0$ ）相比，后者需小初值或符号条件保证全局解；与 SQG 方程（ $\theta_t + u \cdot \nabla \theta + (-\Delta)^{1/2} \theta = 0$ ）相比，后者存在有限时间爆破的开放问题。本文模型的优势在于其线性结构更清晰（ $D_\alpha$ 是一阶常系数微分算子），劣势在于其非局部分析难度更高（旋转破坏各向同性）。本文的代价是：全局解的 $$ H^3 $$ 范数估计（定理 4.1）含有双重指数因子 $\exp[C I(t)]$ ，比 NS 的多项式增长更弱；且衰减率（定理 5.1/5.2）仅对非正初涡量成立，而 NS 的衰减对任意初值成立。此外，作者未讨论 $\alpha \to 0$ 或 $\alpha \to \frac{\pi}{2}$ 的奇摄动极限，也未涉及 $$ L^p $$ 初值或临界空间（如 $B^{2}_{p,1}$ ）下的适定性，这些是未来工作的自然延伸。

2. 关于物理真空中粘性气体恒星演化方程的局部存在性定理

On a Local Existence Theorem for the Evolution Equation of Viscous Gaseous Stars in a Physical Vacuum

作者: Demin Wang, Jiawen Zhang, Shengguo Zhu

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22822v1

摘要翻译: 本文研究三维可压 Navier-Stokes-Poisson 方程组在自引力粘性气体恒星情形下的自由边界问题，其中粘性系数具有退化性。对球对称等熵运动，我们建立了经典解的局部适定性。此处所得解在直至移动边界的整个区域内光滑，并对所有绝热指数 $γ>\frac{4}{3}$ 捕获 Lane-Emden 恒星构型的物理真空边界行为。

Comment: 61 pages

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

1. 研究对象和背景

本文研究三维自引力粘性气体恒星在物理真空边界下的自由边界问题，其数学模型为带退化粘性系数的可压Navier–Stokes–Poisson方程组（CNSP）：

\begin{cases} \rho_t + \operatorname{div}(\rho u) = 0 & \text{in } \Omega(t), \\ (\rho u)_t + \operatorname{div}(\rho u \otimes u) + \nabla P + \rho \nabla \Psi = \operatorname{div} T & \text{in } \Omega(t), \\ \rho > 0 & \text{in } \Omega(t), \\ \rho = 0 & \text{on } \partial \Omega(t), \\ V(\partial \Omega(t)) = u \cdot n(t) & \text{on } \partial \Omega(t), \\ (\rho, u)|_{t=0} = (\rho_0, u_0) & \text{in } \Omega := \Omega(0) = B_1. \end{cases}

其中密度 $\rho$ 、速度 $$u$$ 、压力 $P = A \rho^\gamma$ （ $$A>0$$ , $\gamma > 1$ ）、引力势 $\Psi$ 满足 $\Delta \Psi = 4\pi G \rho$ 且 $\lim_{|x|\to\infty} \Psi = 0$ ；粘性应力张量 $T = 2\mu(\rho) D(u) + \lambda(\rho) \operatorname{div} u \, I_3$ ，其系数为退化幂律形式： $\mu(\rho) = a_1 \rho^\delta$ , $\lambda(\rho) = a_2 \rho^\delta$ ，参数满足 $$a_1 > 0$$ , $2a_1 + 3a_2 \ge 0$ , $\delta \in (0,1)$ 。物理真空边界由条件 $-\infty < \partial_n c_0^2 < 0$ 刻画，其中 $c_0 = \sqrt{P'(\rho_0)} = \sqrt{A \gamma \rho_0^{\gamma-1}}$ 为初始声速，该条件等价于 $\rho_0^{\gamma-1} \sim \operatorname{dist}(x, \partial \Omega)$ 在边界附近成立——即对球对称情形， $\rho_0^{\gamma-1}(r) \sim 1 - r$ （ $$r = |x|$$ ）。此退化结构导致方程在边界处呈现强奇异性： $\rho^{\gamma-1}$ 的一阶导数在边界非零，而 $\rho$ 自身消失，致使经典能量方法失效，且标准Sobolev嵌入在边界附近完全失效。关键未知量为 $(\rho, u)$ 在移动边界 $\partial \Omega(t)$ 上的正则性：要求解在闭包 $\overline{\Omega(t)}$ 上处处光滑（classical solution），而非仅属 $$H^1$$ 或强解（strong solution）；尤其需验证 $u \in C^2(\overline{\Omega(t)})$ ，且满足精确的物理真空边界行为（如 $D_\eta u + \frac{2a_2}{2a_1 + a_2} \frac{u}{\eta}|_{r=1} = 0$ ）。已有理论中，Coutand–Shkoller与Jang–Masmoudi解决了无重力、无粘性（Euler）情形的局部适定性，但其框架依赖于非退化椭圆结构；Luo–Xin–Zeng（2016）及Li–Wang–Xin（2014）虽处理了粘性情形，但仅得到强解（ $u \in L^\infty_t H^1_x \cap L^2_t H^2_x$ ），不保证 $$u$$ 在边界连续可微；Xin–Zhang–Zhu（2021）与Chen–Zhang–Zhu（2023）在1D/2D Saint-Venant系统中实现了经典解，但其加权空间构造严重依赖低维几何与特殊算子分解，无法直接推广至3D CNSP。本文试图填补的核心空白是：在三维、球对称、含自引力、退化粘性、且 $\gamma > 4/3$ 的全参数范围内，建立首个满足物理真空边界条件的经典解局部适定性，即证明存在 $$T^* > 0$$ 使得 $(\rho, u) \in C^2([0,T^*] \times \overline{\Omega(t)})$ ，且所有导数在 $\partial \Omega(t)$ 处有界并满足物理边界律。该问题的重要性在于：它构成了恒星演化数学理论的基础环节——物理真空对应真实恒星大气的急剧稀薄化，其动力学稳定性、激波形成与坍缩机制均依赖于解在真空界面的精细正则性；而 $\gamma > 4/3$ 恰好是恒星流体静力学平衡存在的临界指数（Lane–Emden解存在当且仅当 $\gamma > 4/3$ ），因此该范围具有明确的天体物理意义。

2. 主要定理或主要结果

在参数范围 $\delta \in (0,1)$ , $\gamma \in (4/3, \infty)$ 下，对球对称初值 $(\rho_0, u_0)$ 满足 $\rho_0^{\gamma-1}(r) \in H^3([0,1])$ 且存在常数 $$K_2 > K_1 > 0$$ 使得 $K_1(1-r) \le \rho_0^{\gamma-1}(r) \le K_2(1-r)$ （即物理真空条件），以及 $$u_0$$ 满足加权能量有限性 $E(0,U) < \infty$ （其中 $E(t,U) = E_{\text{in}}(t,U) + E_{\text{ex}}(t,U)$ 定义为： $E_{\text{in}}(t,U) = \|\zeta r(U, D_\eta U, D_\eta U_t, D_\eta^2 U, D_\eta^3 U)(t)\|_{L^2([0,1])}^2$ ， $E_{\text{ex}}(t,U) = \|\rho_0^{1/2}(U, U_t)(t)\|_{L^2([1/2,1])}^2 + \|\rho_0^{\delta/2}(D_\eta U, D_\eta U_t)(t)\|_{L^2([1/2,1])}^2 + \|\rho_0^{(3/2 - \varepsilon_0)(\gamma-1)}(D_\eta^2 U, D_\eta^3 U)(t)\|_{L^2([1/2,1])}^2$ ， $\varepsilon_0 \in (0, \min\{ (3 - 1/(\gamma-1))/2, (1-\delta)/(\gamma-1), 1/100 \})$ ），本文证明：存在 $$T^* > 0$$ （仅依赖于 $a_1, a_2, A, \delta, \gamma, \varepsilon_0, \rho_0, u_0, K_1, K_2, G$ ），使得问题 (1.22) 存在唯一经典解 $(U, \eta)(t,r)$ 满足

E(t,U) + t D(t,U) \in L^\infty(0,T^*), \quad D(t,U) \in L^1(0,T^*),

其中耗散泛函 $D(t,U) = D_{\text{in}}(t,U) + D_{\text{ex}}(t,U)$ 定义为 $D_{\text{in}}(t,U) = \|\zeta r(U_{tt}, D_\eta^2 U_t, D_\eta^4 U)(t)\|_{L^2([0,1])}^2$ ， $D_{\text{ex}}(t,U) = \|\rho_0^{1/2} U_{tt}(t)\|_{L^2([1/2,1])}^2 + \|\rho_0^{(3/2 - \varepsilon_0)(\gamma-1)}(D_\eta^2 U_t, D_\eta^4 U)(t)\|_{L^2([1/2,1])}^2$ ；且满足几何约束 $\eta_r, \eta/r \in [1/2, 3/2]$ 于 $[0,T^*] \times [0,1]$ ，以及精确边界条件 $U|_{r=0} = 0$ , $\left(D_\eta U + \frac{2a_2}{2a_1 + a_2} \frac{U}{\eta}\right)|_{r=1} = 0$ 。进一步，若初始映射 $\eta_0(r)$ 满足 $((\eta_0)_r, \eta_0/r) \in [\delta_*, \delta^*]$ 且加权范数 $\|\zeta r D_r^k \eta_0\|_{L^2} + \|\rho_0^{1/2} \eta_0\|_{L^2} + \|\rho_0^{\delta/2} (\eta_0)_r\|_{L^2} + \|\chi^\sharp \rho_0^{(\gamma-1)(3/2 - \varepsilon_0)} \partial_r^l \eta_0\|_{L^2} < \infty$ （ $k=0,\dots,4$ , $l=2,\dots,4$ ），则结论仍成立，仅 $$T^*$$ 依赖于 $\delta_*, \delta^*$ 。

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

证明的核心障碍源于三重耦合奇异性：(1) 坐标奇异性（ $$r=0$$ 处的球坐标拉普拉斯项 $$2u/r$$ ）；(2) 边界退化（ $\rho \to 0$ 导致椭圆算子 $\operatorname{div}(\rho^\delta \nabla)$ 系数消失）；(3) 物理真空导致的非均匀权重（ $\rho^{\gamma-1} \sim 1-r$ 使标准 $$H^k$$ 范数无法控制边界行为）。作者的真正创新点在于构建了一套分区域、多尺度、自适应权重的能量框架，并发展出针对退化高阶椭圆算子的新型估计技术。

第一，内区（ $r \in [0,1/2]$ ）采用基于三维拉格朗日坐标的 $$H^k$$ 能量：定义 $E_{\text{in}}$ 中的算子 $D_\eta f = f_r / \eta_r$ 及其迭代 $D_\eta^k f$ ，其本质是将三维梯度 $\nabla$ 在球对称下投影为径向微分算子（引理B.1证实 $|D_\eta f|^2 + |f/\eta|^2 \sim |\nabla f|^2$ ）。由此 $E_{\text{in}}$ 控制的是 $\|\zeta U\|_{H^3} + \|\zeta U_t\|_{H^1}$ 类范数，其时间权重 $t D_{\text{in}}$ 用于吸收 $U_{tt}$ 项的低正则性。

第二，外区（ $r \in [1/2,1)$ ）引入精细加权 Sobolev 空间：关键观察是，为控制最高阶导数 $D_\eta^4 U$ ，需将其嵌入 $$C^2$$ ，故利用 Hardy 不等式 $\|f\|_{W^{3,1}} \lesssim \|\rho_0^{\alpha} \partial_r^4 f\|_{L^2}$ 要求 $\alpha \ge 1/2$ ；而由动量方程乘以 $\rho^{\gamma-\delta-1}$ 后的椭圆重构（公式(1.34)）：

(2a_1 + a_2) \rho^{\gamma-1} D_\eta^2 U = \rho^{\gamma-\delta} U_t + \text{低阶项},

推得 $\rho^{\alpha} D_\eta^4 U$ 的 $$L^2$$ 可控性等价于 $\rho^{\alpha + 2 - 2\delta} U_{tt} \in L^2$ ，结合 $\rho_0^{\gamma-1} \sim 1-r$ ，得出 $\alpha = (3/2 - \varepsilon_0)(\gamma-1)$ 是最优选择（ $\varepsilon_0$ 用于规避 $\gamma = 4/3$ 的临界情形）。由此 $E_{\text{ex}}$ 和 $D_{\text{ex}}$ 中的权重 $\rho_0^{1/2}, \rho_0^{\delta/2}, \rho_0^{(3/2 - \varepsilon_0)(\gamma-1)}$ 分别匹配 $U, D_\eta U, D_\eta^2 U$ 的退化阶。

第三，克服线性化困难的关键工具是“交叉导数嵌入”（Appendix D）：证明 $\|D_\eta^2 f\|_{L^2_w} \lesssim \|D_\eta f\|_{L^2_{w'}} + \|f\|_{L^2_{w''}}$ ，其中权重 $$w,w',w''$$ 满足 $w \sim w' \rho_0^{\beta}, w'' \sim w' \rho_0^{\gamma}$ ，这允许将高阶 $D_\eta$ 导数用低阶加权范数控制，避免了传统方法中因 $\rho$ 退化导致的 Gronwall 爆炸。

第四，Picard 迭代的收敛性依赖于线性化问题的全局适定性（§2.3）与一致估计（§3）：作者设计了修正 Galerkin 方案，在截断区间 $$[0,R_n]$$ 上构造近似解，并通过先验估计 $\sup_{t \in [0,T]} E(t) + \int_0^T D(t) dt \le C(E(0))$ 确保迭代序列在加权空间中强收敛，其中常数 $$C$$ 独立于截断参数 $$n$$ 。

4. 与之前工作的比较

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

（1）解的正则性维度：此前所有三维结果均为强解（如 Luo–Xin–Zeng [28] 的 $u \in L^\infty_t H^1_x \cap L^2_t H^2_x$ ），其 $$u$$ 在边界仅 $$H^1$$ 连续，无法定义 $\nabla u \cdot n$ ；而本文首次实现 $u \in C^2(\overline{\Omega(t)})$ ，即 $$u$$ 及其一、二阶导数在闭包上连续，从而严格验证物理真空边界律 $D_\eta u + \frac{2a_2}{2a_1 + a_2} \frac{u}{\eta}|_{r=1} = 0$ 。这一提升代价是必须采用分区域加权能量 $E_{\text{in}} + E_{\text{ex}}$ ，而不能使用单一 Sobolev 范数。

（2）参数适用范围：Jang [17] 在常粘性（ $\delta = 0$ ）下建立了强解局部适定性，但要求 $\gamma > 3/2$ ；本文将 $\gamma$ 下限降至 $$4/3$$ （Lane–Emden 解存在阈值），且允许 $\delta \in (0,1)$ 的全退化范围（此前 Li–Wang–Xin [22] 仅处理 $\delta = 1$ 的1D Saint-Venant）。代价是 $\delta < 1$ 时需 $\varepsilon_0$ 满足 $\varepsilon_0 < (1-\delta)/(\gamma-1)$ ，即退化越强（ $\delta$ 越大），允许的 $\gamma$ 范围越窄。

（3）时空尺度与初值要求：Coutand–Shkoller [9] 的 Euler 解要求初值 $c_0 \in H^3$ （即 $\rho_0^{\gamma-1} \in H^3$ ），但需 $\rho_0^{\gamma-1} \sim (1-r)^k$ 且 $k \ge 4$ 才能保证 $c_0 \in H^3$ ，此时气体无法加速（违反物理真空定义）；本文仅需 $$k=1$$ （即 $\rho_0^{\gamma-1} \sim 1-r$ ），完美匹配物理真空，但为此必须放弃全局 $$H^3$$ 控制，转而采用外区加权 $H^4_{\rho_0^{(3/2-\varepsilon_0)(\gamma-1)}}$ 估计。时间尺度 $$T^*$$ 依赖于初值加权能量 $$E(0,U)$$ ，而非 $$H^k$$ 范数，这反映了退化系统的本质：小 $$L^2$$ 初值未必导致长时存在，而加权小初值才保证局部适定。

（4）方法论差异：Xin–Zhang–Zhu [40] 在1D中利用“交叉导数”技巧控制 $\partial_x^2 u$ ，但其权重为 $\rho^{\theta}$ （ $\theta > 0$ ）；本文将该思想升维并精细化，发展出 $D_\eta$ 算子族与 $\rho_0^{\alpha}$ 权重的协同机制（见Appendix D），这是处理3D球对称退化椭圆性的首创工具。从提供内容中无法确认本文结果能否推广至非球对称情形，因全文分析严格依赖 $D_\eta$ 算子的径向结构及引力势 $\Psi$ 的显式积分表达式（公式(1.15)）。

3. 静置液滴问题的整体适定性：二维Navier-Stokes流

Global Well-Posedness of sessile drop problem: 2D Navier-Stokes Flow

作者: Xiaoding Yang

类别: math.AP

发布日期: 2026-06-21

arXiv ID: 2606.22384v1

摘要翻译: 我们证明了在表面张力、Navier滑移边界条件及动态接触点律作用下，二维不可压Navier–Stokes方程的静置液滴平衡态的小扰动具有整体适定性。主要困难在于平衡流形存在水平平移退化性，为此我们采用由正交性条件确定的移动极坐标系以消除该退化性。随后，我们通过无压弱解的Galerkin构造、压力的恢复、高阶估计以及压缩映射论证，建立了局部适定性。将该局部理论与我们在先前工作中获得的整体能量-耗散估计相结合，得到唯一整体解及相应的指数衰减估计。

Comment: 无

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

1. 研究对象和背景

本文研究二维不可压纳维-斯托克斯方程在具有表面张力、Navier滑移边界条件与动态接触点律的静置液滴自由边界问题中的全局适定性。研究对象是位于一维平坦固体基底上的二维粘性不可压流体液滴，其占据随时间演化的区域

\Omega(t) := \{z \in \mathbb{R}^2 : 0 < r < \rho(t, \theta)\},

其中 $z = (\theta, r)$ 是以基底中点为极点的极坐标， $\rho(t, \cdot) : [0, \pi] \to \mathbb{R}_+$ 为未知自由表面径向函数；自由表面定义为

\Sigma(t) := \{(\theta, r) : \theta \in [0, \pi],\, r = \rho(t, \theta)\},

固体-流体交界（即接触点）位于 $\Sigma_s := \{(\theta, r) : \theta = 0 \text{ 或 } \theta = \pi\}$ 。流体状态由速度场 $u(t, \cdot)$ 与压力 $P(t, \cdot)$ 描述，满足黏性应力张量 $S(P, u) = P I - \mu D u$ （其中 $\mu > 0$ ， $Du = \nabla u + (\nabla u)^T$ ），并服从如下完整系统：

\begin{cases} \partial_t u + u \cdot \nabla u + \operatorname{div} S(P, u) = 0 & \text{in } \Omega(t), \\ \operatorname{div} u = 0 & \text{in } \Omega(t), \\ S(P, u) \nu = (g \rho \sin \theta - \sigma H(\rho)) \nu & \text{on } \Sigma(t), \\ \partial_t \rho = \dfrac{1}{\rho} u \cdot N & \text{on } \Sigma(t), \\ (S(P, u) \cdot \nu - \beta u) \cdot \tau = 0 & \text{on } \Sigma_s, \\ u \cdot \nu = 0 & \text{on } \Sigma_s, \\ W(\partial_t \rho(\pi, t)) = \sigma \dfrac{\rho'}{(\rho^2 + (\rho')^2)^{1/2}}(\pi, t) + [[\gamma]], \\ W(\partial_t \rho(0, t)) = -\sigma \dfrac{\rho'}{(\rho^2 + (\rho')^2)^{1/2}}(0, t) + [[\gamma]], \end{cases}

其中 $\nu$ , $\tau$ 分别为外法向与单位切向， $\sigma > 0$ 为表面张力系数， $H(\rho) = \dfrac{2(\rho')^2 - \rho \rho'' + \rho^2}{(\rho^2 + (\rho')^2)^{3/2}}$ 是自由表面平均曲率的两倍， $\beta > 0$ 为固壁滑移摩擦系数， $N = -\rho'(\theta) \hat{e}_\phi + \rho(\theta) \hat{e}_r$ ， $W \in C^2(\mathbb{R})$ 是接触点速度响应函数的反函数（严格递增微分同胚且 $$W(0)=0$$ ）， $[[\gamma]] = \gamma_{sv} - \gamma_{sf}$ 为固-气与固-液界面能差，且假设 Young 条件 $|[[\gamma]]| < \sigma$ 成立以保证平衡态存在。

该问题的核心困难源于静置液滴平衡流形的水平平移退化性：对任意常数 $$c$$ ，若 $(u=0, P=P_0, \rho=\rho_0)$ 是平衡解，则 $(u=0, P=P_0, \rho=\rho_0(\cdot - c))$ 亦为平衡解，导致重力-毛细能量的二阶变分在全扰动空间上非强制。这一退化使得标准线性化方法失效，且无法直接应用基于谱间隙或正则性提升的标准能量方法。已有工作如 Guo–Tice [10] 建立了先验能量耗散估计，Schweizer [18] 处理了固定接触角 $\pi/2$ 下的 Navier-Stokes 接触点问题，Tice–Zheng [24] 研究了 Stokes 流中移动接触点的局部适定性，Solonnikov [20,21] 在加权 Hölder 空间中处理了静态与动态接触点的 Navier-Stokes 解存在性。但所有这些工作均未解决“同时包含动态接触点、动态接触角、平移退化性及非线性 Navier-Stokes 效应”的二维静置液滴全局适定性问题——这正是本文试图填补的根本空白。

2. 主要定理或主要结果

在平衡接触角 $\omega_{\text{eq}} \in (0, \pi)$ 且 $\varepsilon_{\max} := \min\{1, -1 + \pi/\omega_{\text{eq}}\} \in (0,1)$ 的假设下，设初始数据 $(u_0, p_0, \xi_0, \partial_t u(0), \partial_t p(0), \partial_t \xi(0), \partial_t^2 u(0), \partial_t^2 \xi(0))$ 满足如下兼容性条件：对 $$j = 0,1,2$$ ，

\begin{cases} \operatorname{div}_{A_0} D_t^j u(0) = 0 & \text{in } \Omega, \\ D_t^j u(0) \cdot \nu = 0 & \text{on } \Sigma_s, \\ D_t^j u(0) \cdot N(0) = \partial_{t}^{j+1} (\rho(0) \xi(0) - n'(0) \xi_3) & \text{on } \Sigma, \end{cases}

以及零均值条件：对 $$k = 0,1,2$$ ，

\int_0^\pi \partial_t^k (\rho(0) \xi(0)) \, d\theta = 0.

定义能量泛函与耗散泛函：

E(t) := \sum_{k=0}^2 \left( \|\partial_t^k u\|_{L^2}^2 + \|\partial_t^k \xi\|_{H^1}^2 \right) + \|u\|_{W^{2,q_+}}^2 + \|\partial_t u\|_{H^{1+\varepsilon_-/2}}^2 + \|p\|_{W^{1,q_+}}^2 + \|\partial_t p\|_{L^2}^2 + |n'(t)|^2 + |n''(t)|^2 + \|\xi\|_{W^{3-1/q_+,q_+}}^2 + \|\partial_t \xi\|_{H^{3/2 + (\varepsilon_- - \alpha)/2}}^2 + \|\partial_t^2 \xi\|_{H^1}^2,

D(t) := \|u\|_{W^{2,q_+}}^2 + \|\partial_t u\|_{W^{2,q_-}}^2 + \sum_{k=0}^2 \left( [\partial_t^{k+1} \xi]_\theta^2 + [\partial_t^k \partial_\theta \xi]_\theta^2 \right) + \|p\|_{W^{1,q_+}}^2 + \|\partial_t p\|_{W^{1,q_-}}^2 + \sum_{k=0}^2 |n^{(k+1)}|^2 + \sum_{k=0}^2 \left( \|\partial_t^k \xi\|_{H^{3/2 - \alpha}}^2 + \|\xi\|_{W^{3-1/q_+,q_+}}^2 + \|\partial_t \xi\|_{W^{3-1/q_-,q_-}}^2 + \|\partial_t^3 \xi\|_{H^{1/2 - \alpha}}^2 \right) + \sum_{k=0}^2 \left( \|\partial_t^k u\|_{H^1}^2 + \|\partial_t^k u\|_{L^2}^2 + [\partial_t^k u \cdot N]_\theta^2 \right),

其中 $0 < \alpha < \varepsilon_- < \varepsilon_+ < \varepsilon_{\max}$ ， $q_\pm = 2/(2 - \varepsilon_\pm) \in (1,2)$ 。则存在充分小常数 $\delta_0 > 0$ ，使得若初始能量满足 $\sqrt{E(0)} + \|\partial_t^2 \xi(0)\|_{W^{2-1/q_+, q_+}} \le \delta_0$ ，则系统 (1.10) 存在唯一全局解 $(u,p,\xi)$ 定义于 $[0,\infty)$ ，满足初值且满足指数衰减估计：

\sup_{t \ge 0} e^{\mu t} E(t) + \int_0^\infty D(t) \, dt \le C E(0),

其中 $$C > 0$$ 、 $\mu > 0$ 为普适常数。特别地，解满足 $\|u(t)\|_{H^1} + \|\xi(t)\|_{H^1} \lesssim e^{-\mu t/2}$ 。

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

本文的核心技术障碍有三：一是平移退化性导致线性化算子在全空间上缺乏谱间隙，致使标准能量估计无法闭合；二是移动极点坐标系引入额外未知函数 $$n(t)$$ ，其演化耦合于速度与自由表面，破坏了经典自由边界问题的解耦结构；三是初始相容性条件因 $$n(t)$$ 和辅助参数 $$m(t)$$ 的介入而高度非平凡，传统构造失效。

作者的真正创新点在于构建了一个兼具几何适配性与分析可控性的伪线性框架。首先，通过正交性条件 $\int_0^\pi \xi \rho_0 \, d\theta = 0$ 与 $\int_0^\pi \xi \xi_s \, d\theta = 0$ （其中 $\xi_s$ 为平移模态）确定移动极点位置 $$n(t)$$ ，从而获得强制性双线性型 $(\xi, \xi)_{1,\Sigma_0} \gtrsim \|\xi\|_{H^1}^2$ ，其显式形式为：

(\rho_1, \rho_2)_{1,\Sigma_0} = g \int_0^\pi \rho_0 \rho_1 \rho_2 \sin \theta \, d\theta + \sigma \int_0^\pi \frac{\rho_0 \rho_1' \rho_2'}{(\rho_0^2 + (\rho_0')^2)^{3/2}} \, d\theta - \sigma \int_0^\pi \frac{\rho_0' \rho_1' \rho_2}{(\rho_0^2 + (\rho_0')^2)^{3/2}} \, d\theta - \sigma \int_0^\pi \frac{\rho_0' \rho_1 \rho_2'}{(\rho_0^2 + (\rho_0')^2)^{3/2}} \, d\theta + \sigma \int_0^\pi \frac{\rho_0'' \rho_0 - (\rho_0')^2 - \rho_0^2}{(\rho_0^2 + (\rho_0')^2)^{3/2}} \rho_1 \rho_2 \, d\theta.

其次，在线性化步骤中，作者摒弃直接分离 $$n(t)$$ 的做法，转而引入辅助参数 $$m(t)$$ 满足 $\int_0^\pi \xi_s \partial_t^2 \xi \, d\theta = 0$ ，从而在伪线性系统 (1.22) 中显式嵌入正交约束，确保线性算子保持强制性。该 $$m(t)$$ 具有显式表达式：

m(t) = \frac{1}{\int_0^\pi \xi_s^2 \, d\theta} \int_0^\pi \xi_s \Big( \frac{1}{\rho} v_d \cdot N + \cdots \Big) \, d\theta,

其存在性与正则性由高阶 Sobolev 嵌入与 $\xi_s \in H^1(0,\pi)$ 保证。第三，针对初始数据构造，作者引入辅助函数 $\xi_5$ ，证明若 $(u_0, \xi_0, p_0)$ 满足 $$t=0$$ 时的方程，则 $(u_0, \xi_0 + C \xi_5, p_0)$ 仍为解，从而将零均值约束 $\int_0^\pi \rho(0) \xi(0) \, d\theta = 0$ 转化为关于 $$C$$ 的线性方程求解问题。

证明框架采用双层 Galerkin 极限：先对光滑逼近序列 $(\eta_k, u_k)$ 和 $(\eta_n, u_n, p_n)$ 构造压力无关弱解 $(v_d, \xi_d)$ 满足弱形式 (1.36)，再通过椭圆估计恢复压力 $$q_d$$ ，继而通过 ODE 积分重构 $(v,q,\xi)$ ，最后执行 $k \to \infty$ 后 $n \to \infty$ 的双重极限，并利用收缩映射原理建立局部解，结合前作 [29] 的全局先验估计得到整体解。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度：

几何框架与退化处理：Guo–Tice [10] 仅建立先验估计，未解决适定性；Schweizer [18] 假设接触角恒为 $\pi/2$ ，从而规避了接触角动态演化与平移退化耦合；本文首次在完全动态接触角 $\omega_{\text{eq}} \in (0,\pi)$ 下，通过移动极坐标系与正交性条件彻底消除平移退化，代价是引入额外未知 $$n(t)$$ 及其导数项 $n'(t)\partial_x u$ 、 $n''(t)\xi_s$ 等非标准低阶项，显著增加能量估计复杂度。
线性化策略：传统方法（如 [30] 的容器问题）在线性化时固定几何，而本文的伪线性系统 (1.22) 显式包含 $$m(t)$$ 及其积分项（如 $\int_0^t R_{1a}(\rho_0, \partial_\theta \eta, \eta) \partial_\theta \partial_t^2 \xi \, ds$ ），这些项在能量估计中需借助 Grönwall 不等式与 $$L^2_t$$ 时间积分控制，而非简单吸收。作者为此发展了新的积分核估计引理（Theorem 4.5），量化了 $$m(t)$$ 对 $\|\partial_t^2 \xi\|_{H^{1/2-\alpha}}$ 的依赖。
函数空间选择：本文采用非均匀 Sobolev 指标 $q_\pm = 2/(2-\varepsilon_\pm)$ ，允许 $\varepsilon_\pm$ 在 $(0,\varepsilon_{\max})$ 内精细调节，以平衡 $W^{2,q_+}$ 速度正则性与 $W^{2,q_-}$ 时间导数正则性，克服 Navier-slip 边界导致的 $W^{2,2}$ 不可及性；而 Solonnikov [20,21] 使用加权 Hölder 空间，无法处理 $$L^2$$ -base 的能量耗散结构。
初始数据构造：[30] 中初始扰动可加任意常数，故零均值易实现；本文因 $\rho = \rho_0 + \xi$ 出现在曲率 $H(\rho)$ 的分母中，破坏平移不变性，迫使作者构造特殊辅助函数 $\xi_5$ 满足 $\int_0^\pi \rho_0 \xi_5 \, d\theta \neq 0$ ，从而将约束转化为可解线性方程。此构造未见于任何先前文献。

从适用范围看，本文结果要求 $\omega_{\text{eq}} \in (0,\pi)$ 且 $\varepsilon_{\max} > 0$ ，即排除 $\omega_{\text{eq}} = 0,\pi$ （完全润湿/去润湿情形），而 Tice–Wu [23] 的 Stokes 稳定性分析允许更广角度但仅限线性流。代价是解的存在性依赖于小初值 $\delta_0$ ，且指数衰减率 $\mu$ 隐含依赖于 $\omega_{\text{eq}}$ 与 $\sigma/\mu$ 比值，作者未给出显式下界。

4. 一种经典的椭圆正则性方法用于几乎调和映射及相关系统

A Classical Elliptic Regularity Approach to Almost Harmonic Maps and Related Systems

作者: Giovanni Di Fratta

类别: math.AP, cond-mat.mtrl-sci, math-ph, math.DG, math.FA

发布日期: 2026-06-21

arXiv ID: 2606.22468v1

摘要翻译: 我们建立了一大类二维非线性椭圆系统的内正则性结果。我们的方法在一个统一的抽象框架内分离出核心可积性机制，该框架围绕一个 Campanato 型离散迭代方案与一个 Caccioppoli 型估计耦合构建。具体而言，我们证明：在任意一类关于尺度变换稳定且满足离散振荡衰减公理的容许对 $(\boldsymbol{u}, \boldsymbol{f})$ 中，映射 $\boldsymbol{u}$ 自动局部 Hölder 连续。此外，所得 Hölder 指数是显式的，并最优地达到由源项 $\boldsymbol{f}$ 的 Lebesgue 可积性所决定的经典 Morrey–Campanato 阈值。这一纯分析框架系统性地规避了现有正则性理论所依赖的 $\mathcal{H}^1$ – $\mathrm{BMO}$ 对偶、Wente 不等式、活动标架法以及共形一致化技术。

Comment: 59 pages

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

1. 研究对象和背景

本文研究二维非线性椭圆系统的内正则性，核心对象是定义在单位开圆盘 $B_1 \subset \mathbb{R}^2$ 上、取值于单位球面 $\mathbb{S}^n \subset \mathbb{R}^{n+1}$ 的弱解 $\boldsymbol{u} \in H^1(B_1, \mathbb{S}^n)$ ，满足如下几乎调和映射方程：

-\Delta \boldsymbol{u} = |\nabla \boldsymbol{u}|^2 \boldsymbol{u} + \boldsymbol{f} \quad \text{in } \mathcal{D}'(B_1),

其中扰动项 $\boldsymbol{f} \in L^q(B_1, \mathbb{R}^{n+1})$ ， $$ q > 1 $$ ，且满足几何相容性条件 $\boldsymbol{f} \cdot \boldsymbol{u} = 0$ a.e.（该条件由 $|\boldsymbol{u}| = 1$ 自动导出：对原方程与 $\boldsymbol{u}$ 取内积得 $-\boldsymbol{u}\cdot\Delta\boldsymbol{u} = |\nabla\boldsymbol{u}|^2$ ，而 $-\boldsymbol{u}\cdot\Delta\boldsymbol{u} = |\nabla\boldsymbol{u}|^2 + \boldsymbol{f}\cdot\boldsymbol{u}$ ，故 $\boldsymbol{f}\cdot\boldsymbol{u}=0$ ）。该方程是Dirichlet能量 $E(\boldsymbol{u}) = \frac{1}{2}\int_{B_1} |\nabla \boldsymbol{u}|^2\,dx$ 在约束 $|\boldsymbol{u}|=1$ 下的Euler–Lagrange方程之扰动形式，广泛出现在液晶Oseen–Frank理论、微磁学（micromagnetics）中——例如Brown静态微磁能量泛函

G(\boldsymbol{u}) = \int_U \left[ \frac{1}{2}|\nabla \boldsymbol{u}|^2 + \kappa \, \boldsymbol{u} \cdot \operatorname{curl} \boldsymbol{u} + \varphi(\boldsymbol{u}) \right] dx - \frac{1}{2} \int_U h[\boldsymbol{u}] \cdot \boldsymbol{u} \, dx,

其临界点即满足上述几乎调和映射结构（见Theorem 1.5）。

该问题的关键分析障碍源于非线性项 $|\nabla \boldsymbol{u}|^2 \boldsymbol{u}$ 的临界二次增长：在标准伸缩变换 $\boldsymbol{u}_r(x) := \boldsymbol{u}(rx)$ 下， $\|\nabla \boldsymbol{u}_r\|_{L^2(B_1)} = r^{-1}\|\nabla \boldsymbol{u}\|_{L^2(B_r)}$ ，而 $\||\nabla \boldsymbol{u}_r|^2 \boldsymbol{u}_r\|_{L^1(B_1)} \sim r^{-2} \|\nabla \boldsymbol{u}\|_{L^2(B_r)}^2$ ，与 $\|-\Delta \boldsymbol{u}_r\|_{L^1(B_1)}$ 同尺度，使方程处于Calderón–Zygmund理论的临界边界。在二维情形下，已知 $$ H^1 $$ 弱解连续，但经典证明严重依赖于刚性代数结构：Hélein通过将非线性重写为 $\Omega \cdot \nabla \boldsymbol{u}$ 形式（其中 $\Omega = \boldsymbol{u} \wedge \nabla \boldsymbol{u}$ 是反对称矩阵），再利用Uhlenbeck规范分解构造散度为零的联络 $\operatorname{div} \Omega = 0$ ，最终诉诸Hardy空间 $\mathcal{H}^1$ –BMO对偶性及Wente不等式。Rivière将此推广至一般反对称 $\Omega \in L^2(B_1; \mathfrak{so}(n+1)\otimes \mathbb{R}^2)$ 的系统 $-\Delta \boldsymbol{u} = \Omega \cdot \nabla \boldsymbol{u}$ ；Müller–Schikorra将其扩展至含源项情形。然而，所有这些方法均以 $\Omega$ 的几何反对称性为前提，其本质作用是生成 $\operatorname{div} \Omega = 0$ 的规范，从而激活补偿紧性机制。

本文试图填补的根本空白在于：是否存在一种不依赖任何几何规范构造、不调用 $\mathcal{H}^1$ –BMO对偶、不使用Wente不等式或共形化的纯分析正则性理论？具体而言，四个未解问题驱动全文：(Q1) 对 $$ L^q $$ 型张力场 $\boldsymbol{f}$ ， $\mathcal{H}^1$ –BMO对偶是否必要？能否仅用经典线性椭圆理论获得显式Hölder指数？(Q2) 反对称性是否本质？能否被纯分析条件 $\operatorname{div} \Omega \in L^q$ 替代？(Q3) 当主部变为各向异性算子 $-\operatorname{div}(A\nabla \boldsymbol{u})$ 且 $\det A \not\equiv \text{const}$ （即不可共形化）时，正则性是否仍成立？(Q4) 从连续性到无穷可微性的提升过程能否避免繁琐中间步骤，形成一个维度无关、无需平稳性假设的直接引导机制？这些问题共同指向一个核心诉求：将正则性机制从“几何代数结构”还原为“底层积分可积性”，并建立一个统一、鲁棒、可扩展的抽象框架。

2. 主要定理或主要结果

本文建立了四个层次递进的主要定理，其结论均在二维单位圆盘 $B_1 \subset \mathbb{R}^2$ 上给出，且所有Hölder指数 $\eta$ 均显式依赖于源项可积性参数 $$ q $$ ：

Theorem 1.1（几乎调和映射的Hölder连续性）：设 $\boldsymbol{u} \in H^1(B_1, \mathbb{S}^n)$ 是方程 $-\Delta \boldsymbol{u} = |\nabla \boldsymbol{u}|^2 \boldsymbol{u} + \boldsymbol{f}$ 的弱解，其中 $\boldsymbol{f} \in L^q(B_1, \mathbb{R}^{n+1})$ ， $$ q > 1 $$ 。则 $\boldsymbol{u}$ 局部Hölder连续，即存在 $\eta = \eta(q) \in (0,1)$ ，使得 $\boldsymbol{u} \in C^{0,\eta}_{\mathrm{loc}}(B_1, \mathbb{R}^{n+1})$ 。更精确地， $\eta$ 达到经典的Morrey–Campanato阈值：若 $\boldsymbol{f} \in L^q$ ，则 $\boldsymbol{u} \in C^{0,\alpha}$ 对所有 $\alpha < 1 - 2/q$ 成立（当 $$ q > 2 $$ ），或 $\alpha < 2(1 - 1/q)$ （由Campanato空间嵌入 $\mathcal{L}^{2,p} \hookrightarrow C^{0,1-2/p}$ ，此处 $$ p = q^* = 2q/(2-q) $$ ）。

Theorem 1.2（广义Rivière型系统的正则性）：设 $\boldsymbol{u} \in H^1 \cap L^\infty(B_1, \mathbb{R}^{n+1})$ 是 $-\Delta \boldsymbol{u} = \Omega \cdot \nabla \boldsymbol{u} + \boldsymbol{f}$ 的弱解，其中 $\Omega \in L^2(B_1; \mathbb{R}^{(n+1)\times(n+1)} \otimes \mathbb{R}^2)$ ， $\boldsymbol{f} \in L^q(B_1, \mathbb{R}^{n+1})$ ，且 $\operatorname{div} \Omega \in L^q(B_1; \mathbb{R}^{(n+1)\times(n+1)})$ ， $$ q > 1 $$ 。则 $\boldsymbol{u} \in C^{0,\eta}_{\mathrm{loc}}(B_1, \mathbb{R}^{n+1})$ ， $\eta$ 同Theorem 1.1。

Theorem 1.3（各向异性系统的正则性）：设 $\boldsymbol{u} \in H^1 \cap L^\infty(B_1, \mathbb{R}^{n+1})$ 满足 $-\operatorname{div}(A\nabla \boldsymbol{u}) = \Omega \cdot \nabla \boldsymbol{u} + \boldsymbol{f}$ ，其中系数矩阵 $A: B_1 \to \mathbb{R}^{2\times 2}$ 满足：(A1) 对称、一致椭圆（ $\lambda|\xi|^2 \leq A(x)\xi\cdot\xi \leq \Lambda|\xi|^2$ ）且Hölder连续 $A \in C^{0,\gamma_0}$ ， $\gamma_0 \in (0,1]$ ；(A2)–(A4) 同Theorem 1.2。则 $\boldsymbol{u} \in C^{0,\eta}_{\mathrm{loc}}(B_1, \mathbb{R}^{n+1})$ ， $\eta$ 仍由 $$ q $$ 显式决定。此结果突破了 $\det A \equiv \text{const}$ 的共形化限制。

Theorem 1.4（维度无关的正则性引导）：设 $U \subset \mathbb{R}^m$ 为开集， $\boldsymbol{u} \in W^{1,2}_{\mathrm{loc}}(U, \mathbb{S}^n) \cap C^0(U, \mathbb{S}^n)$ 是 $-\Delta \boldsymbol{u} = |\nabla \boldsymbol{u}|^2 \boldsymbol{u} + \boldsymbol{f}$ 的弱解， $\boldsymbol{f} \in L^q_{\mathrm{loc}}(U, \mathbb{R}^{n+1})$ 。则：(1) 若 $2 \leq q \leq m$ ，则 $\boldsymbol{u} \in W^{1,q^*}_{\mathrm{loc}}$ ，其中 $$ q^* = mq/(m-q) $$ （ $$ m=q $$ 时 $$ q^* $$ 任意大）；(2) 若 $$ q > m $$ ，则 $\boldsymbol{u} \in C^{1,\eta}_{\mathrm{loc}}$ ，其中 $\eta = 1 - m/q$ ；(3) 若 $\boldsymbol{f} \in C^{k-1,\beta}_{\mathrm{loc}}$ ， $k \geq 1$ ， $\beta \in (0,1)$ ，则 $\boldsymbol{u} \in C^{k+1,\beta}_{\mathrm{loc}}$ 。特别地， $\boldsymbol{f} \in C^\infty \implies \boldsymbol{u} \in C^\infty$ 。该引导链完全独立于平稳性假设，且在 $$ m=2 $$ 时由Theorem 1.1提供初始连续性，故为无条件成立。

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

本文的核心创新在于构建了一个去几何化的、基于Campanato迭代的抽象正则性原理，其技术支柱有三，彻底绕开了传统方法的代数瓶颈：

1. 抽象振荡衰减公理与Campanato迭代闭合机制：作者定义了变分振荡泛函 $\Psi_p(\boldsymbol{u}, B_r) := \inf_{c \in \mathbb{R}^{n+1}} \int_{B_r} |\boldsymbol{u}(y) - c|^p \, dy$ ，并确立了一个抽象类 $\mathcal{S}$ 的公理化框架：若 $(\boldsymbol{u}, \boldsymbol{f}) \in \mathcal{S}$ 在伸缩下封闭，且满足离散振荡衰减公理

\Psi_p(\boldsymbol{u}, B_\theta) \leq \gamma \, \Psi_p(\boldsymbol{u}, B_1) + \kappa \, \|\boldsymbol{f}\|_{L^q(B_1)}^p, \quad \theta \in (0,1), \, \gamma \in (0,1),

则 $\boldsymbol{u}$ 自动局部Hölder连续，且指数 $\eta$ 由 $\gamma$ 和 $$ p $$ 显式确定（Theorem 3.3）。关键突破在于，该公理的验证不依赖于任何规范构造，而是通过一个精巧的耦合Caccioppoli估计实现。

2. 耦合Caccioppoli型估计：这是全文最核心技术。对任意半径 $$ r $$ ，令 $$ h_r $$ 为 $\boldsymbol{u}$ 在 $$ B_r $$ 上的调和延拓（即 $\Delta h_r = 0$ in $$ B_r $$ ， $h_r = \boldsymbol{u}$ on $\partial B_r$ ）。作者证明了如下关键不等式（式(1.11)）：

\|\nabla(\boldsymbol{u} - h_r)\|_{L^q(B_r)} \leq C \, \|\nabla \boldsymbol{u}\|_{L^2(B_r)} \, \left( r^m \, \Psi_p(\boldsymbol{u}, B_r) \right)^{1/p} + C r \, \|\boldsymbol{f}\|_{L^q(B_r)},

其中指数满足 $$ 1/q = 1/2 + 1/p $$ 。在二维情形下，这强制 $$ p = q^* = 2q/(2-q) $$ ，恰好是Sobolev共轭指数。这一“巧合”是迭代闭合的根源：它使得 $\Psi_p$ 的尺度行为与 $\|\nabla \boldsymbol{u}\|_{L^2}$ 的能量控制完美匹配，从而无需Wente不等式即可完成Campanato迭代。该估计的证明依赖于对方程 $-\Delta \boldsymbol{u} = \Omega \cdot \nabla \boldsymbol{u} + \boldsymbol{f}$ 的代数重写（式(1.10)）：利用球面约束 $|\boldsymbol{u}|=1$ ，将 $|\nabla \boldsymbol{u}|^2 \boldsymbol{u}$ 表达为 $\Omega \cdot \nabla \boldsymbol{u}$ ，再对任意常向量 $c \in \mathbb{R}^{n+1}$ 应用Leibniz法则，得到

-\Delta \boldsymbol{u} = \operatorname{div}(\Omega(\boldsymbol{u} - c)) + (\boldsymbol{u} \wedge \boldsymbol{f})(\boldsymbol{u} - c) + \boldsymbol{f}.

此式将非线性项转化为散度项与低阶项，使测试函数技巧得以应用。

3. 散度条件 $\operatorname{div} \Omega \in L^q$ 的解析替代机制：在Theorem 1.2的证明中，作者并未尝试构造散度为零的规范，而是直接利用 $\operatorname{div} \Omega \in L^q$ 这一纯分析条件。关键引理（Fundamental Lemma 3.7）表明，该条件足以控制 $\Omega \cdot \nabla \boldsymbol{u}$ 在 $$ L^q $$ 意义下的“伪散度”效应，从而将 $\Omega \cdot \nabla \boldsymbol{u}$ 视为一个 $$ L^q $$ 型扰动源，融入前述耦合Caccioppoli框架。这揭示了反对称性的真正作用——仅为生成 $\operatorname{div} \Omega = 0$ 的代数工具——而正则性机制本身仅需 $\operatorname{div} \Omega$ 具有足够可积性。

4. 各向异性情形的局部系数冻结法：对于 $-\operatorname{div}(A\nabla \boldsymbol{u})$ ，作者采用局部冻结系数策略：在小球 $$ B_r $$ 上，用常系数矩阵 $$ A(x_0) $$ 近似 $$ A(x) $$ ，误差由 $[A]_{C^{0,\gamma_0}} r^{\gamma_0}$ 控制。这导致一个带有小参数 $r^{\gamma_0}$ 的扰动项，其贡献被吸收进迭代的收缩因子 $\gamma$ 中，从而维持公理成立。此方法严格要求维度 $$ m < 4 $$ （文中Section 5注明），因高维下Sobolev嵌入失效。

5. 维度无关引导的Calderón–Zygmund/Schauder级联：Theorem 1.4的证明完全脱离临界非线性，仅依赖标准线性理论：一旦 $\boldsymbol{u}$ 连续，则 $|\nabla \boldsymbol{u}|^2 \boldsymbol{u} \in L^{q/2}$ ，故右端 $\in L^{q/2}$ ；由Calderón–Zygmund， $\boldsymbol{u} \in W^{2,q/2}$ ；再由Sobolev嵌入与Schauder理论逐级提升，形成 $L^q \to W^{2,q/2} \to W^{1,q^*} \to C^{1,\eta} \to C^{k+1,\beta}$ 的清晰链条。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开，其推进本质是将正则性理论从‘几何代数依赖’转向‘纯分析可积性依赖’：

1. 与Hélein/Rivière/Müller–Schikorra的对比（绕过H¹–BMO与Wente）：Hélein [Hél90] 和 Rivière [Riv06] 的方法以 $\Omega = \boldsymbol{u} \wedge \nabla \boldsymbol{u}$ 的反对称性为基石，通过Uhlenbeck规范分解获得 $\tilde{\Omega}$ 满足 $\operatorname{div} \tilde{\Omega} = 0$ ，再利用 $\mathcal{H}^1$ –BMO对偶处理 $\tilde{\Omega} \cdot \nabla \tilde{\boldsymbol{u}}$ 。本文则证明，只要 $\operatorname{div} \Omega \in L^q$ ，即可直接启动Campanato迭代，完全规避规范构造与 $\mathcal{H}^1$ 理论。代价是：Hélein方法可处理 $\boldsymbol{f} \in \mathcal{H}^1$ （比 $$ L^1 $$ 更优），而本文要求 $\boldsymbol{f} \in L^q$ ， $$ q>1 $$ ；但本文给出了显式Hölder指数 $\eta(q)$ ，而传统方法对此无定量刻画。

2. 与Bethuel/Bet92的对比（散度条件的强化与普适化）：Bethuel [Bet92] 在平均曲率方程中提出 $\operatorname{div} \Omega \in L^{2,1}$ （Lorentz空间），本文将其推广至标准Lebesgue空间 $$ L^q $$ ， $$ q>1 $$ ，并证明该条件在广义Rivière系统中已充分。更重要的是，本文将此条件置于一个统一的抽象框架中，使其成为可移植的正则性判据，而非针对特定方程的孤立技巧。

3. 与共形化方法的对比（突破各向异性壁垒）：当主部为 $-\operatorname{div}(A\nabla \boldsymbol{u})$ 且 $\det A \not\equiv \text{const}$ 时，共形变换 $y = \Phi(x)$ 无法将方程化为 $-\Delta_y$ 形式，因 $\det D\Phi$ 无法抵消 $\det A$ 的变化。现有规范理论（如Rivière）在此失效。本文通过局部系数冻结法，在 $$ m=2 $$ 下成功处理此情形（Theorem 1.3），适用范围显著扩大。代价是：该方法依赖于 $$ A $$ 的Hölder连续性 $C^{0,\gamma_0}$ ，而共形化方法对 $$ A $$ 的正则性要求更低（仅需 $L^\infty$ ）；且本文证明中隐含维度上限 $$ m<4 $$ ，而共形化无此限制。

4. 与Evans/Bethuel的高维部分正则性对比（二维全局正则性 vs 高维奇异集估计）：Evans [Eva91] 和 Bethuel [Bet93] 针对 $m \geq 3$ 的平稳解证明奇异集 $\mathcal{S}$ 满足 $\dim_{\mathcal{H}} \mathcal{S} \leq m-2$ ，但无法排除奇异点。本文聚焦 $$ m=2 $$ ，目标是全局Hölder连续性，其结论更强（无奇异点），但适用维度受限。然而，Theorem 1.4的引导机制是维度无关的，这意味着本文框架可与高维部分正则性结果结合：先用Evans方法获得 $\boldsymbol{u} \in C^0_{\mathrm{loc}}(U \setminus \mathcal{S})$ ，再用本文引导升级光滑性，从而在正则区域获得最优正则性。文中未显式展开此比较，但从逻辑依赖看，这是自然延伸。

综上，本文的推进是范式性的：它将正则性归因于源项 $\boldsymbol{f}$ 和联络散度 $\operatorname{div} \Omega$ 的可积性，而非几何结构的刚性。其代价是牺牲了某些最优性（如 $\mathcal{H}^1$ 源项），但换来了前所未有的灵活性、显式性与物理适用性（如微磁学Theorem 1.5）。

5. 奇异定常Navier-Stokes流：例子与稳定性

Singular stationary Navier-Stokes flows: examples and stability

作者: Zachary Bradshaw, Dakota Palmer

类别: math.AP

发布日期: 2026-06-21

arXiv ID: 2606.22291v1

摘要翻译: Landau解若沿竖直轴取向，则构成一族单参数的精确、自相似、轴对称、无旋（swirl-free）解，满足受原点处向上指向的动量点源所驱动的三维定常Navier-Stokes方程；其在原点处具有孤立奇点。其他奇异定常解亦可由类似框架导出。其中某些变体已被提议作为物理场景的模型。例如，Squire首先发现的一个解被提议为油类径向排出表面层所夹带流体的模型；另一解——Serrin的旋转涡旋（swirling vortex）——则展现出与某些龙卷风共有的定性特征，如包含中心下沉气流与外围上升气流的双单元结构，以及旋涡运动。本文第一目标是对这些及其他例子进行细致分析，尤其当存在边界时。在此方向上，我们发现了若干新的、具有物理动机的解类，并识别出物理学文献与数学文献之间若干新的联系。这些例子中许多定义在半空间 $\mathbb{R}^3_+$ 上，但关于奇异定常解的大量数学文献则针对全空间 $\mathbb{R}^3$ 。本文第二目标是在含边界的区域上建立这些解的渐近稳定性。我们最一般的结论依赖于一种基于最终正则性（eventual regularity）的全新渐近稳定性方法。作为特例，我们证明了一类受Serrin旋转涡旋启发的定常解在轴对称扰动下是稳定的。这一结果要求一个新颖的观察：这些解过于奇异——文献中称之为“Type III”型——以致无法直接适用于现有方法。

Comment: 52 pages, 6 figures

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

1. 研究对象和背景

本文研究三维定常不可压纳维-斯托克斯方程（SNS）在具有边界（特别是半空间 $\mathbb{R}^3_+$ ）和无界区域（如 $\mathbb{R}^3$ ）上的奇异稳态解及其渐近稳定性。核心方程为

\begin{cases} -\Delta V + (V \cdot \nabla) V + \nabla q = f & \text{in } \Omega, \\ \nabla \cdot V = 0 & \text{in } \Omega, \\ V|_{\partial \Omega} = V_*(x) & \text{on } \partial \Omega, \end{cases}

其中 $\Omega \subseteq \mathbb{R}^3$ 是开连通区域（典型取 $\mathbb{R}^3$ 、 $\mathbb{R}^3_+ = \{x_3 > 0\}$ 或具 Lipschitz 边界的有界/外部域）， $V : \Omega \to \mathbb{R}^3$ 为速度场， $q : \Omega \to \mathbb{R}$ 为压力， $$f$$ 为外力（常取为狄拉克型源项）， $$V_*$$ 为给定边界值。研究对象是满足 $$(-1)$$ -齐次性（即 $V_\lambda(x) = \lambda V(\lambda x)$ ）、轴对称性（ $\partial_\phi V = 0$ ）且可能含孤立奇点的弱解，其奇异性集中于点（如原点）、半轴（如 $\{x_h = 0, x_3 \ge 0\}$ ）或整条轴（如 $$x_h = 0$$ ）。关键未知量是速度场 $$V$$ 在奇点附近的渐近结构（如 $|V(x)| \sim |x|^{-1}$ 或 $|V_h(x)| \sim |x_h|^{-1}, |V_3(x)| \sim |x_h|^{-2/3+}$ ）及其对扰动的鲁棒性。

此类奇异解在流体力学中具有明确物理动机：Landau 解建模淹没射流；Squire 表面排放流描述油层表面径向扩散所卷吸的流体；Serrin 涡旋以双胞结构（中心下沉流+外围上升流+涡量）模拟龙卷风。然而，数学上这些解长期受限于全空间 $\mathbb{R}^3$ 的框架，而物理场景天然涉及边界（如自由表面、固壁）。现有理论存在三重空白：（1）解的存在性与分类缺失：既有工作（如 Li-Li-Yan [28–30]）聚焦于 $\mathbb{R}^3$ 上沿轴或半轴的 Type II 奇异解（ $|V| \lesssim \rho^{-1+\varepsilon}$ ， $\rho = |x_h|$ ），但未系统处理半空间 $\mathbb{R}^3_+$ 中满足物理边界条件（如无穿透 $V \cdot e_3 = 0$ 或非零径向边界流）的对应解；（2）奇异性尺度鸿沟：Serrin 涡旋等 Type III 解（ $|V_h| \sim \rho^{-1}$ ，严格 $O(\rho^{-1})$ ）因不属 $L^{3,\infty}(\mathbb{R}^3_+)$ 而被排除在所有已知稳定性理论之外（如 [34, 18, 19] 要求 $\|V\|_{L^{3,\infty}} < \epsilon^*$ ）；（3）稳定性理论域限制：现有渐近稳定性结果（如 Karch-Pilarczyk [18] 的半群法、Kukavica-Ziane [19] 的 Fourier 分裂法）或仅适用于 $\mathbb{R}^3$ ，或要求解足够正则（ $L^{3,\infty}$ 小），无法覆盖 $\mathbb{R}^3_+$ 上 Type III 奇异解的轴对称扰动。本文旨在填补上述空白：一方面通过统一框架（Li-Li-Yan 框架）重构并推广 Squire、Serrin 等物理模型至带边界情形，识别新解类（如 Type II.5）；另一方面发展新稳定性方法，突破 $L^{3,\infty}$ 小性假设与全空间限制。

2. 主要定理或主要结果

本文证明了三类核心结果：

（一）定常解的全局 $$L^2$$ 渐近稳定性（Theorem 1.3）：设 $\Omega \subseteq \mathbb{R}^3$ 为 $\mathbb{R}^3$ 、 $\mathbb{R}^3_+$ 或具 Lipschitz 边界的有界/外部域， $V \in L^{3,\infty}(\Omega)$ 为散度自由定常解。则存在 $\epsilon^* > 0$ ，使得当 $\|V\|_{L^{3,\infty}(\Omega)} < \epsilon^*$ 时，对任意初值 $u_0 \in L^2_\sigma(\Omega)$ ，扰动方程（PNS）的弱解 $$u$$ 满足

\lim_{t \to \infty} \|u(t)\|_{L^2(\Omega)} = 0.

该结果将 [18, 19] 的 $\mathbb{R}^3$ 结论推广至 $\mathbb{R}^3_+$ 等带边界域，首次确立 Squire 表面排放流及其带涡扩展的稳定性。

（二）Type II.5 解的轴对称 $$L^2$$ 渐近稳定性（Theorem 1.5）：设 $$V$$ 为定义在 $\mathbb{R}^3$ 或 $\mathbb{R}^3_+$ 上的 Type II.5 解（Definition 1.4），即满足轴对称、自相似，且在 $$x_h = 0$$ 附近有

|V_{x_3}(x_h, \pm 1)| = O(|x_h|^{-2/3+}), \quad |V_{x_i}(x_h, \pm 1)| = O(|x_h|^{-1}) \ (i=1,2),

并存在常数 $c_{i,\pm}$ 使得 $V_i(x_1, 0, \pm 1) = c_{i,\pm}/x_1 + O(|x_1|^{-2/3+})$ 。则对任意满足广义能量不等式（5.3）的轴对称弱解 $$u$$ ，有

\lim_{t \to \infty} \|u(t)\|_{L^2(\Omega)} = 0.

此为首个针对奇点沿直线且水平分量严格 $O(\rho^{-1})$ 的解的稳定性结果，无需 $$V$$ 本身小，仅依赖扰动的全局能量不等式。

（三）非定常解的渐近稳定性（Theorem 1.6）：设 $\Omega$ 为“可容许域”（admissible domain，即 $L^{3,\infty}(\Omega)$ 中存在小数据整体适定性，如 $\mathbb{R}^3_+$ ）， $V \in L^\infty(0,\infty; L^{3,\infty}(\Omega))$ 为散度自由非定常解。则存在 $\epsilon^* > 0$ ，使得当 $\|V\|_{L^\infty_t L^{3,\infty}_x(\Omega)} < \epsilon^*$ 时，扰动解 $$u$$ 满足

\lim_{t \to \infty} \|u(t)\|_{L^2(\Omega)} = 0.

该结果通过“最终正则性”（eventual regularity）实现：弱解 $$u$$ 在有限时间后满足 $\|u(t)\|_{L^{3,\infty}(\Omega)} < \epsilon^*$ ，从而触发小数据适定性，迫使所有项衰减。

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

证明的核心障碍在于：传统稳定性方法（半群法、Fourier 分裂法）均无法处理 Type III 奇异解——前者要求 $$V$$ 定常且 $L^{3,\infty}$ 小，后者依赖 $\mathbb{R}^3$ 上的 Fourier 变换且需 $$V$$ 小。作者提出三项创新性技术：

（一）Type II.5 解的分解与导数控制（Section 5）：针对 $$V = a + b$$ 的分解，其中 $b \in L^{3,\infty}(\Omega)$ 且 $$a$$ 的水平分量主导奇异性，但其 $$x_3$$ -导数受控：

|\partial_{x_3} a|(x) \le \frac{M}{|x|^2}, \quad \forall x \notin \{x_h = 0\}.

此性质源于 Type II.5 解的自相似与轴对称结构—— $$V_h$$ 的奇异性主部与 $$x_3$$ 无关，故 $\partial_{x_3} V_h$ 具有更高阶衰减。该控制使线性项 $(V \cdot \nabla) u$ 在能量估计中可被 $\|\nabla u\|_{L^2}^2$ 吸收，规避了 $V_h \sim \rho^{-1}$ 导致的标准乘子估计失效。

（二）轴对称下的端点二维 Hardy 不等式（Section 5）：在轴对称设定下，利用柱坐标 $(\rho, \phi, z)$ 将 $\mathbb{R}^3$ 视为 $\mathbb{R}^2_+ \times \mathbb{S}^1$ ，对 $$u$$ 的轴对称分量应用 Hardy 不等式

\int_{\mathbb{R}^2_+} \frac{|u_h(\rho,z)|^2}{\rho^2} \, d\rho dz \le C \int_{\mathbb{R}^2_+} |\nabla_h u_h|^2 \, d\rho dz,

其中 $\nabla_h$ 为水平梯度。此不等式在一般三维情形失效，但轴对称性消除了 $\phi$ -导数项，使 $\rho^{-1}$ 奇异性可被梯度控制，成为处理 $V_h \sim \rho^{-1}$ 关键。

（三）基于最终正则性的非定常稳定性框架（Section 6）：放弃直接估计 $$V$$ ，转而证明扰动 $$u$$ 自身在有限时间 $$T_0$$ 后进入小数据区域：

\exists T_0 > 0, \quad \|u(t)\|_{L^{3,\infty}(\Omega)} < \epsilon^*, \quad \forall t \ge T_0.

此由 Kato 型论证实现：利用 Lorentz-Sobolev 不等式 $\|u\|_{L^{6,2}} \le C_S \|\nabla u\|_{L^2}$ 及能量不等式，证明 $\int_t^{t+1} \|\nabla u(\tau)\|_{L^2}^2 \, d\tau \to 0$ ，结合插值得 $L^{3,\infty}$ 范数衰减。此后， $$u$$ 满足小数据适定性（Theorem 6.5），其 mild 解表示

u u(t) = e^{t\Delta} u(T_0) - \int_{T_0}^t e^{(t-s)\Delta} \mathbb{P} \big[ (u\cdot\nabla)u + (u\cdot\nabla)V + (V\cdot\nabla)u \big](s) \, ds

中所有非线性项与线性项同阶小，最终全项趋于零。

4. 与之前工作的比较

本文与既有工作的差异体现在四个维度：

（一）奇异性尺度与解类拓展：[34] 仅处理 Type II 解（ $|V| \lesssim \rho^{-1+\varepsilon}$ ），其 $L^{3,\infty}$ 可积性保证半群法适用；本文定义 Type II.5 类（ $|V_h| \sim \rho^{-1}$ ， $|V_3| \sim \rho^{-2/3+}$ ），涵盖 Serrin 涡旋及 [28] 中部分 $\mathbb{R}^3$ 解，填补了 $O(\rho^{-1})$ 奇异性的稳定性空白。相较 [17] 对 $\mathbb{R}^3_+$ 上无滑移解不存在性的证明，本文通过允许内部奇点（如 Serrin 解）绕过该障碍。

（二）稳定性假设的放宽：[18, 19] 要求 $\|V\|_{L^{3,\infty}} < \epsilon^*$ ，而 Theorem 1.5 完全取消此假设，仅需扰动满足全局能量不等式，代价是限定于轴对称扰动；Theorem 1.6 将小性要求从 $$V$$ 转移至 $$u$$ 的最终状态，适用范围更广但需额外验证 $\Omega$ 的可容许性（如 $\mathbb{R}^3_+$ 需引用 [14, 23] 的小数据适定性）。

（三）域适应性突破：[19] 的 Fourier 分裂法严格限于 $\mathbb{R}^3$ ；本文 Theorem 1.3 将 [18] 的半群法推广至 $\mathbb{R}^3_+$ 等边界域，首次给出半空间中 Squire 流的稳定性；Theorem 1.6 的最终正则性法亦适用于 $\mathbb{R}^3_+$ ，克服了 Fourier 方法的域限制。

（四）物理模型与数学框架的对接：既有数学工作（如 [28–30]）未系统联系 Squire、Serrin 的物理模型；本文 Section 3.3–3.5 证明 Squire 流属于 Li-Li-Yan 框架，并指出其奇点位于半轴（Type III），进而构造带涡扩展（Section 3.4）；同时将 Serrin 涡旋纳入 Type II.5 类，建立龙卷风模型与现代分析理论的严格接口。代价是：Theorem 1.5 的轴对称限制排除了非对称扰动；Theorem 1.6 的可容许域条件尚未在一般 Lipschitz 域上验证；广义能量不等式（5.3）的普适性仍待确认（文中注明‘未来工作’）。

6. Navier-Stokes 方程的有限链 CKN-坏尺度计数：标准偏微分方程闭合与典范探测器实现

Finite-Chain CKN-Bad Scale Counting for Navier-Stokes: Standard PDE Closure and Canonical Detector Realization

作者: Runlong Yu

类别: math.AP

发布日期: 2026-06-19

arXiv ID: 2606.21783v1

摘要翻译: 本文证明了三维不可压 Navier-Stokes 方程的适定弱解的 Caffarelli–Kohn–Nirenberg 坏尺度的有限链计数定理。主要的标准偏微分方程结果给出了有限个 CKN-坏尺度的加权尺寸上界，其上界由非负通道代价组成，包括垂直单分量集中度、环形泄漏、压力尾项以及压力–通量–能量残差。所证明的闭合机制是在完全局部临界界条件下的定性单分量紧性：垂直分量充分小时，可推出更小半径处的 CKN 小性。随后，本文给出了同一有限窗口计数思想的典范探测器实现。原始抽象探测器并未与标准偏微分方程通道相认同；相反，我们定义了一种新的修正典范探测器，其坐标由能量、通量、压力尾、保留的低频压力模态以及有限维残差构成。对此修正探测器，我们证明了上实现性、下审计性、CKN 提取性以及有限链坏尺度计数性。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压Navier-Stokes方程

\partial_t u + (u \cdot \nabla) u - \Delta u + \nabla p = 0, \quad \nabla \cdot u = 0

的适合弱解（suitable weak solution）在局部尺度上的奇点结构。核心对象是Caffarelli–Kohn–Nirenberg（CKN）意义下的‘坏尺度’（bad scale）：对给定中心点 $$z_0 = (x_0, t_0)$$ 和尺度序列 $r_k = \lambda^k r_0$ （其中 $0 < \lambda < 1$ , $k = 0,1,\dots,K$ ），定义归一化场 $$u_k(y,s) = r_k u(x_0 + r_k y, t_0 + r_k^2 s)$$ , $$p_k(y,s) = r_k^2 p(x_0 + r_k y, t_0 + r_k^2 s)$$ ，并考察CKN量

\Phi_k(\rho) := \rho^{-2} \int_{Q_\rho} |u_k|^3 \, dy ds + \rho^{-2} \int_{-\rho^2}^{0} \int_{B_\rho} |p_k - (p_k)_{B_\rho}(s)|^{3/2} \, dy ds,

其中 $Q_\rho = B_\rho \times (-\rho^2, 0)$ , $(p_k)_{B_\rho}(s) = |B_\rho|^{-1} \int_{B_\rho} p_k(y,s) \, dy$ 。CKN经典epsilon正则性定理断言：存在绝对常数 $\varepsilon_{\text{CKN}} > 0$ ，使得若 $\Phi_k(\rho) < \varepsilon_{\text{CKN}}$ ，则解在更小的物理柱体 $Q_{c \rho r_k}(z_0)$ 内正则。因此，满足 $\Phi_k(\rho_M) \ge \varepsilon_{\text{CKN}}$ 的尺度 $$r_k$$ 被称为CKN坏尺度，其集合记为 $B = \{ 0 \le k \le K : \Phi_k(\rho_M) \ge \varepsilon_{\text{CKN}} \}$ 。问题本质是逆向书账（converse bookkeeping）：若一个有限几何链 $\{r_k\}_{k=0}^K$ 中存在多个坏尺度，哪些标准PDE通道（standard PDE channels）必然发生集中？这直接关联到奇点形成的动力学机制——是否可由单一分量主导？是否压力尾部或通量缺陷可独立触发奇异性？现有理论中，CKN准则本身是充分条件而非必要条件；而Leray-Hopf解的全局正则性仍是未解难题。关键空白在于：缺乏一个非循环（non-tautological）、可量化、且与标准PDE能量学直接对应的坏尺度计数框架。已有工作如[3,4,5,7]确立了CKN准则，[13,14,16–20]发展了单分量正则性判据，但均未将坏尺度数量与可测量的物理通道成本（如垂直速度分量、环形泄漏、压力尾部、压力-通量-能量残差）建立显式、加权、有限链尺度的不等式联系。尤其，以往估计常隐含地将右端项设为CKN量本身（如 $\Phi_k(\rho_M)$ ），导致计数定理空洞化。本文旨在填补这一空白：构造一组彼此独立、非负、且不包含CKN核心范数的‘诚实通道成本’（honest channel costs），并证明其加权和能控制坏尺度总数。研究设定严格限定于局部归一化柱体 $Q_1 = B_1 \times (-1,0)$ 上满足全局部临界界（full local critical bound）的适合弱解：

\Psi_k(1) := A_k(1) + E_k(1) + C_k(1) + D_k(1) \le M,

其中 $A_k(1) = \operatorname{ess\,sup}_{-1<s<0} \int_{B_1} |u_k|^2$ , $E_k(1) = \int_{Q_1} |\nabla u_k|^2$ , $C_k(1) = \int_{Q_1} |u_k|^3$ , $D_k(1) = \int_{-1}^{0} \int_{B_1} |p_k - (p_k)_{B_1}(s)|^{3/2}$ 。此假设强于仅控制 $\Phi_k(1)$ ，因它同时约束能量、耗散、三阶速度及压力振荡，为紧性论证提供足够支承。本文试图推进的空白正是：在该标准紧性框架下，给出首个非循环、可计算、且分离各物理通道贡献的有限链坏尺度上界。

2. 主要定理或主要结果

本文在经典输入（Classical Inputs 4.1 和 4.2）成立的前提下，证明了以下主要结果：

定性单分量epsilon定理（Theorem 4.3）：对任意 $M < \infty$ ，存在 $\varepsilon_3(M) > 0$ 与 $0 < \rho_M < 3/4$ ，使得若归一化解 $$(u,p)$$ 满足 $\Psi(1) \le M$ 且 $C_3 = \int_{Q_1} |u_3|^3 \le \varepsilon_3(M)$ ，则必有 $\Phi(\rho_M) < \varepsilon_{\text{CKN}}$ 。
非循环局部闭合引理（Theorem 5.1）：令 $\varepsilon_{\text{close}}(M) := \varepsilon_3(M)$ 。若 $$(u,p)$$ 满足 $\Psi(1) \le M$ 且标准成本 $C_{\text{std}} = C_3 + L^{\text{ann}} + P^{\text{tail}} + R^{\text{PFE}} \le \varepsilon_{\text{close}}(M)$ ，则 $\Phi(\rho_M) < \varepsilon_{\text{CKN}}$ ；其逆否形式为：
$\Phi(\rho_M) \ge \varepsilon_{\text{CKN}} \implies C_{\text{std}} \ge \varepsilon_{\text{close}}(M).$
加权有限链坏尺度计数定理（Theorem 6.3）：对尺度链 $r_k = \lambda^k r_0$ ， $k=0,\dots,K$ ，设 $\Psi_k(1) \le M$ 对所有 $$k$$ 成立， $B = \{k : \Phi_k(\rho_M) \ge \varepsilon_{\text{CKN}}\}$ ，则对任意非负权重 $w_k \ge 0$ ，有
$\sum_{k \in B} w_k \le \frac{1}{\varepsilon_{\text{close}}(M)} \sum_{k=0}^{K} w_k \left( C_{3,k} + L^{\text{ann}}_k + P^{\text{tail}}_k + R^{\text{PFE}}_k \right).$
无权重计数推论（Corollary 6.4）：取 $$w_k = 1$$ ，得坏尺度个数上界
$\# B \le \frac{1}{\varepsilon_{\text{close}}(M)} \sum_{k=0}^{K} \left( C_{3,k} + L^{\text{ann}}_k + P^{\text{tail}}_k + R^{\text{PFE}}_k \right).$
有限链CKN小尺度准则（Corollary 7.1）：若加权总成本满足
$\sum_{k=0}^{K} w_k \left( C_{3,k} + L^{\text{ann}}_k + P^{\text{tail}}_k + R^{\text{PFE}}_k \right) < \varepsilon_{\text{close}}(M) \sum_{k=0}^{K} w_k,$
则必存在 $k^* \in \{0,\dots,K\}$ 使得 $w_{k^*} > 0$ 且 $\Phi_{k^*}(\rho_M) < \varepsilon_{\text{CKN}}$ ，从而由CKN准则知解在 $Q_{c \rho_M r_{k^*}}(z_0)$ 内正则。

所有结果中， $\rho_M$ 仅依赖于 $$M$$ ，且 $\varepsilon_{\text{close}}(M) = \varepsilon_3(M)$ 的具体形式未显式给出，但其存在性由紧性论证保证；通道成本定义明确： $C_{3,k} = \int_{Q_1} |(u_k)_3|^3$ ；环形泄漏 $L^{\text{ann}}_k = \int_{-\lambda^2}^{0} \int_{\lambda A_\chi} \left( |u_k|^2 + |\nabla u_k|^2 + |u_k|^3 + |p_k||u_k| \right)$ ；压力尾部成本 $P^{\text{tail}}_k = \| [\eta, R_i R_j](u_{k,i} u_{k,j}) \|_{L^{3/2}(Q_{\rho_M})}^{3/2} + P^{\text{harm}}_k + P^{\text{trunc}}_k$ ；PFE残差 $R^{\text{PFE}}_k$ 是非负的、不包含CKN核心项的剩余量。

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

证明的核心障碍在于避免循环论证（tautology）：任何将CKN量 $\Phi_k(\rho_M)$ 直接作为右端项的估计都丧失物理意义。作者的真正创新点在于构建并验证一组‘诚实’（honest）、‘非循环’（non-tautological）的通道成本，并证明它们能协同驱动CKN小性。整个证明框架由三层次递进构成：

第一层：经典紧性输入的精确调用与分离。作者并未重新证明紧性，而是将两个经典结论显式列为‘Classical Input’：(i) 在全临界界 $\Psi_k(1) \le M$ 下，适合弱解序列在 $L^3_{\text{loc}}(Q_1)$ 中强收敛，且其压力分解的均值自由部分在 $L^{3/2}(Q_\rho)$ 中强收敛；(ii) 所有满足 $\Psi_v(1) \le M$ 且 $$v_3 = 0$$ 的紧性极限 $$(v,q)$$ ，必在某个仅依赖 $$M$$ 的半径 $\rho_M$ 处满足 $\Phi_{v,q}(\rho_M) < \varepsilon_{\text{CKN}}/2$ 。此分离至关重要——它将技术难点（压力紧性、零分量极限的正则性）外包给已知理论，使本文焦点纯粹落在‘如何从通道成本推出紧性’上。

第二层：单分量闭合机制的定量实现。定理4.3的证明采用反证法：假设 $$C_3$$ 小但 $\Phi(\rho_M)$ 不小，则存在序列 $\int_{Q_1} |u^{(j)}_3|^3 \to 0$ ，由Classical Input 4.1得 $u^{(j)} \to v$ 强于 $L^3_{\text{loc}}$ ，故 $$v_3 = 0$$ ；再由Classical Input 4.2得 $\Phi_{v,q}(\rho_M) < \varepsilon_{\text{CKN}}/2$ ；最后利用速度与压力的强收敛性得 $\Phi_{u^{(j)},p^{(j)}}(\rho_M) \to \Phi_{v,q}(\rho_M)$ ，矛盾。此处的关键估计是压力收敛的传递性：由于 $u^{(j)} \to v$ 强于 $L^3_{\text{loc}}$ ，Calderón-Zygmund算子 $$R_i R_j$$ 在 $L^{3/2}$ 上有界，故主动压力 $R_i R_j(\eta u^{(j)}_i u^{(j)}_j)$ 强收敛；而调和余项 $h^{(j)}$ 在 $L^{3/2}(B_\rho)$ 中紧致（由调和函数内估计），故整体压力振荡收敛。这确保了 $\Phi$ 的连续性，是闭合得以成立的基石。

第三层：通道成本的构造性定义与非负性保障。作者对每个成本项施加严格禁令以杜绝循环：(i) 环形泄漏 $L^{\text{ann}}_k$ 仅积分在支撑集 $\lambda A_\chi$ 上，该区域与核心球 $B_{\rho_M}$ 分离（因 $\rho_M < 3/4$ 且 $\eta = 1$ 在 $B_{3/4}$ 上），故不接触CKN核心；(ii) 压力尾部 $P^{\text{tail}}_k$ 的交换子项 $\| [\eta, R_i R_j](u_i u_j) \|_{L^{3/2}(Q_{\rho_M})}^{3/2}$ 由命题3.5控制为 $C_{\eta,\rho_M} \int_{(-1,0) \times A_\eta} |u|^3$ ，其支撑亦在 $A_\eta$ ，与 $B_{\rho_M}$ 无交；调和尾 $P^{\text{harm}}_k$ 和截断误差 $P^{\text{trunc}}_k$ 被明确定义为非负且不得等于完整的尺度不变压力振荡；(iii) PFE残差 $R^{\text{PFE}}_k$ 被定义为‘压力分解失配、通量平衡残差、局域能量不等式缺陷、有限秩PFE截断误差或简化探测器比较残差’，并强制其不包含 $\rho_M^{-2} \int_{Q_{\rho_M}} |u|^3$ 或 $\rho_M^{-2} \int_{-\rho_M^2}^{0} \int_{B_{\rho_M}} |p - (p)_{B_{\rho_M}}(s)|^{3/2}$ 。这些构造性定义确保了 $C_{\text{std},k}$ 是一个真正的、可物理诠释的成本 ledger，其每一项都对应一个明确的PDE现象（垂直动量、环形能量逸出、压力非局部效应、系统建模误差），而非CKN量的伪装。

4. 与之前工作的比较

本文与既有工作的核心区别在于‘非循环性’（non-tautology）与‘通道分离性’（channel separation）的严格实现，而非单纯的技术改进。与CKN原始工作[3,4,5]相比，后者建立了正则性的充分条件 $\Phi < \varepsilon_{\text{CKN}}$ ，但未提供坏尺度的定量计数；本文则首次给出了坏尺度数量的显式上界，且右端项完全剥离了CKN量本身。与Lin的简化证明[7]及后续综述[15]相比，这些工作聚焦于重述和推广CKN准则，而本文转向其逆问题——坏尺度的成因分析。在单分量正则性方向，[13,14,16–20]证明了如 $u_3 \in L^{2,\infty}$ 或 $\partial_3 u \in L^{2,\infty}$ 可保证正则，但这些是全局或时间-空间混合范数的判据，未涉及局部尺度链上的坏尺度分布；本文则在局部归一化框架下，将 $C_3 = \int_{Q_1} |u_3|^3$ 这一纯空间-时间三阶范数作为闭合通道，并量化了其与坏尺度的关系。与压力正则性工作[6,11]相比，本文的压力尾部成本 $P^{\text{tail}}_k$ 并非直接使用压力本身，而是通过交换子 $[\eta, R_i R_j]$ 捕捉压力生成的非局部性，且明确禁止其退化为CKN压力项，这比单纯要求 $p \in L^{3/2}$ 更精细。与Onsager/反常耗散文献[8,9,10]中的‘通量缺陷’概念相比，本文的 $R^{\text{PFE}}_k$ 将此类缺陷形式化为一个非负的、可添加的ledger项，而非仅作启发式讨论。与作者前期工作[22–28]相比，本文刻意‘窄化’（narrower）：前期探讨了更广泛的‘有限尺度审计框架’（finite-scale audit framework），可能包含更强的尺度一致正则性主张；而本文明确声明‘不声称尺度一致正则性定理’，仅记录一个保守的、基于标准PDE通道的计数定理及其典范探测器实现。代价与局限在于：(i) 闭合机制目前仅由垂直分量 $$C_3$$ 驱动；定理5.1及6.3虽将 $L^{\text{ann}}_k, P^{\text{tail}}_k, R^{\text{PFE}}_k$ 列入ledger，但未证明它们可独立闭合（Remark 5.2, 5.3明确指出‘更强的PFE闭合引理尚未被证明’）；(ii) 全临界界 $\Psi_k(1) \le M$ 的假设强于CKN所需的 $\Phi_k(1) \le M$ ，这虽为紧性所必需，但也限制了适用范围；(iii) $\varepsilon_{\text{close}}(M)$ 的具体衰减率未知，仅知其存在，故定量预测能力受限。从提供内容中无法确认本文结果能否推广至二维或更高维，亦无法确认 $\rho_M$ 关于 $$M$$ 的具体依赖关系。

7. 欧拉方程势流自相似解中伪声速曲线的几何结构

Geometric Structures of Pseudo-Sonic Curves in Self-Similar Solutions of the Euler Equations for Potential Flow

作者: Gui-Qiang G. Chen, Mikhail Feldman, Xin Gao, Wei Xiang

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

发布日期: 2026-06-19

arXiv ID: 2606.21793v1

摘要翻译: 我们研究二维欧拉方程势流自相似解中伪声速曲线的几何结构，允许非均匀超声速状态存在。数学上，所满足的二阶势流方程为混合型双曲-椭圆型方程，其退化性沿伪声速曲线发生。本文发展了严格的分析方法以研究此类自相似解中伪声速曲线的几何结构。首先，我们证明：若在每一点处伪速度均垂直于该曲线，则伪声速曲线必为圆。接着，我们分析伪声速点处伪速度不垂直于曲线的一般情形，并研究伪声速曲线邻域内流线的几何性质。随后，我们在解的局部行为满足自然假设的前提下，建立两个定理，给出伪声速点处伪速度垂直于曲线的充分条件。这些结果给出了伪声速曲线几何结构的精确刻画。最后，我们将所发展的理论应用于具有非均匀入流的激波反射-绕射问题。我们证明：若解在伪超声速区域或伪亚声速区域内是均匀入流解的 $$C^2$$ -小扰动，则伪声速曲线必为一段圆弧。特别地，密度与速度分别为该伪声速圆弧的半径与圆心。此外，我们证明解在伪亚声速区域内直至声速圆弧（除点 $$P_1$$ 外）具有 $C^{2,α}$ 正则性。本文所发展的技巧与思想有望适用于其他涉及类似混合型退化性的非线性问题。

Comment: 39 pages, 8 figures

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

1. 研究对象和背景

本文研究二维自相似势流中伪音速曲线 $\Gamma_{\text{sonic}}$ 的几何结构，其数学模型为二阶非线性偏微分方程：

\operatorname{div}\big(\rho(|D\varphi|^2, \varphi) D\varphi\big) + 2\rho(|D\varphi|^2, \varphi) = 0, \quad (\xi,\eta) \in \mathbb{R}^2,

其中密度函数由绝热状态方程给出：

\rho(|D\varphi|^2, \varphi) = \left(K - \varphi - \frac{1}{2}|D\varphi|^2\right)^{\frac{1}{\gamma-1}}, \quad \gamma > 1,

$$K$$ 为伯努利常数， $D = (\partial_\xi, \partial_\eta)$ ， $\varphi$ 为伪速度势。相应地，伪速度场定义为 $(U,V) = D\varphi$ ，物理速度为 $(u,v) = D\varphi + (\xi,\eta)$ ，音速平方为

c^2(|D\varphi|^2, \varphi) = (\gamma - 1)\rho^{\gamma-1} = (\gamma - 1)\left(K - \varphi - \frac{1}{2}|D\varphi|^2\right),

故伪音速点满足退化条件：

|D\varphi| = c(|D\varphi|^2, \varphi).

该方程在 $|D\varphi| < c$ 区域为椭圆型（伪亚音速区 $\Omega_+$ ），在 $|D\varphi| > c$ 区域为双曲型（伪超音速区 $\Omega_-$ ），而在 $\Gamma_{\text{sonic}}$ 上发生混合型退化——系数矩阵的符号从正定突变为不定，且退化阶为一阶（即特征值以线性速率趋于零）。此退化导致经典椭圆/双曲理论失效：解在 $\Gamma_{\text{sonic}}$ 处无法期望 $$C^2$$ 正则性；事实上，[1,9] 中构造的均匀入流解已证明第二阶导数存在跳跃间断。因此， $\Gamma_{\text{sonic}}$ 的几何形状（如是否为圆弧、曲率是否恒定、凸性方向）直接决定了解的局部适定性、正则性延拓能力及跨音速奇性传播机制。

现有理论几乎全部建立在“均匀超音速入流”假设之上：此时 $\Gamma_{\text{sonic}}$ 被先验设为圆弧，且中心与半径由常数音速 $$c_0$$ 和常数速度 $$(u_2,v_2)$$ 确定（见 [1,9,12]）。例如，在正则激波反射-衍射问题中，若入流状态为常数 $(\rho_0,u_0,v_0)$ ，则 $\Gamma_{\text{sonic}}$ 必为圆弧，且 $c \equiv c_0$ 、 $(U,V) \equiv (u_0 - \xi, v_0 - \eta)$ 在其上成立，从而可借助坐标平移与缩放将问题化归为标准圆对称情形，并利用分离变量、Fredholm 理论或 Schauder 估计获得 $C^{2,\alpha}$ 解（除端点 $$P_1$$ 外）。然而，当入流非均匀时（如含小扰动 $\delta \rho_0(\xi,\eta), \delta u_0(\xi,\eta), \delta v_0(\xi,\eta)$ ），该先验圆假设不再成立； $\Gamma_{\text{sonic}}$ 可能弯曲、扭曲甚至分裂，导致退化边界不再是等距曲面，传统正则性分析框架彻底崩溃。更严峻的是，非均匀性破坏了方程的自相似不变性—— $\Gamma_{\text{sonic}}$ 不再是相似不变集，其几何必须作为未知自由边界被同时求解。本文旨在填补这一根本性空白：不预设 $\Gamma_{\text{sonic}}$ 的任何几何形式（仅假设 $C^{1,1}$ 正则性），而通过解本身的局部行为反推其必为圆弧，并精确刻画其曲率、中心与半径如何由局部物理量 $$(c(Q), U(Q), V(Q))$$ 决定。

2. 主要定理或主要结果

本文核心结论分为五组定理，均建立在解满足基本适定性假设之上： $\varphi \in C^2(\Omega_- \cup \Omega_+ \cup \Gamma^0_{\text{sonic}})$ ，且在 $\Gamma_{\text{sonic}}$ 上连续匹配 $\varphi$ , $D\varphi$ （即 (2.14)）。

定理 2.1（法向速度蕴含圆形结构）：若对每个 $Q \in \Gamma_{\text{sonic}}$ ，伪速度 $$(U,V)(Q)$$ 均为 $\Gamma_{\text{sonic}}$ 在 $$Q$$ 处的单位法向量，则 $\Gamma_{\text{sonic}}$ 必为圆弧；进一步，若 $$c > 0$$ 且 $\Gamma_{\text{sonic}}$ 连通，则 $c \equiv c_2$ 为常数，且 $\Gamma_{\text{sonic}}$ 位于半径为 $$c_2$$ 的圆上，其圆心为 $$(u_2, v_2)$$ ，其中 $(U,V) = (u_2 - \xi, v_2 - \eta)$ 在 $\Gamma_{\text{sonic}}$ 上成立。

定理 2.2 与 2.3（切向速度情形的刚性刻画）：若 $$(U,V)$$ 处处切于 $\Gamma_{\text{sonic}}$ ，则 $\partial_n c = 0$ 当且仅当 $\Gamma_{\text{sonic}}$ 为直线段；若 $\partial_n c \neq 0$ ，则 $\Gamma_{\text{sonic}}$ 在伪超音速侧 $\Omega_-$ 视角下严格凸（即曲率 $\kappa > 0$ ）。

定理 2.4（流线几何）：设 $$Q$$ 为伪音速点且 $U(Q) \neq 0$ ，局部流线 $\eta = F(\xi) \in C^2$ 。若 $\partial_\perp c(Q) = 0$ ，则 $$F''(Q) = 0$$ （流线在 $$Q$$ 处无曲率）；若 $\partial_\perp c(Q) \neq 0$ 且 $\Omega_- \subset \{\eta > F(\xi)\}$ ，则 $$F''(Q) > 0$$ （流线在 $$Q$$ 处向上凸）。

定理 2.5（局部二次接触判据）：对任意 $Q \in \Gamma_{\text{sonic}}$ ，定义辅助函数 $\psi_Q = \varphi - \varphi^Q_2$ ，其中 $\varphi^Q_2$ 是以 $$Q$$ 为中心、具有常数伪速度 $$(U,V)(Q)$$ 的二次近似势（见 (2.17)）。若 $\partial_\perp c(Q) = 0$ ，且存在 $a_1 \in (0, \frac{1}{\gamma+1})$ 与 $\varepsilon > 0$ ，使得沿单位法向 $$n_Q$$ 或单位切向 $\tau_Q$ 满足

\psi_Q(Q + t n_Q) < a_1 \left( \frac{n_Q \cdot (U,V)(Q)}{c(Q)} \right)^2 t^2 \quad \text{或} \quad \psi_Q(Q + s \tau_Q) < a_1 \left( \frac{\tau_Q \cdot (U,V)(Q)}{c(Q)} \right)^2 s^2,

对 $t \in (0,\varepsilon)$ 或 $s \in (0,\varepsilon)$ 成立，则 $$Q$$ 必为例外点（即 $$(U,V)(Q)$$ 法向于 $\Gamma_{\text{sonic}}$ ）。

定理 2.6（小角度判据）：设 $\theta_1$ 为 $\Gamma_{\text{sonic}}$ 在 $$Q$$ 处内法向 $$n_Q$$ 与伪速度方向 $$e_v$$ 的夹角。若 $\psi_Q$ 在伪亚音速区 $\Omega_+$ 满足二次下界

\psi_Q > a_2 \big((\xi - \xi_*) \cos \sigma_* + (\eta - \eta_*) \sin \sigma_*\big)^2 - \epsilon \big((\eta - \eta_*) \cos \sigma_* - (\xi - \xi_*) \sin \sigma_*\big)^2,

其中 $$a_2 > 0$$ , $\epsilon \geq 0$ , 且 $\epsilon \leq \epsilon_0$ , $|\theta_1| \leq \varepsilon_1$ （ $\epsilon_0, \varepsilon_1$ 依赖于 $$r$$ , $$a_2$$ , $\alpha$ , $\|f\|_{C^{1,1}}$ , $\|\psi_Q\|_{C^{1,\alpha}(B_r(Q) \cap \Omega_+)}$ ），则 $$Q$$ 必为例外点。

应用结果（§8）：对非均匀入流的正则激波反射-衍射问题，若解在伪超音速区或伪亚音速区为 [9] 中均匀解的 $$C^2$$ 小扰动，则 $\Gamma_{\text{sonic}}$ 必为圆弧，且其半径 $$R$$ 与圆心 $$(x_c, y_c)$$ 满足 $$R = c(Q)$$ , $(x_c, y_c) = (U(Q) + \xi_*, V(Q) + \eta_*)$ 对所有 $Q \in \Gamma_{\text{sonic}}$ 成立（即 $$c$$ 与 $$(U,V)$$ 在 $\Gamma_{\text{sonic}}$ 上为常数）；进而，解在伪亚音速区 $\Omega_+$ 中 $C^{2,\alpha}$ 正则，直至 $\Gamma_{\text{sonic}}$ （除端点 $$P_1$$ 外）。

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

本文的核心障碍在于：伪音速曲线 $\Gamma_{\text{sonic}}$ 是一个未知自由边界，其几何与解的正则性相互耦合，且方程在该曲线上发生一阶退化，导致标准椭圆/双曲估计完全失效。作者的真正创新点在于构建了一套基于“局部二次接触分析”的几何反演框架，其关键工具与技术如下：

1. 例外点判据与法向性转化：作者并未直接处理 $\Gamma_{\text{sonic}}$ 的曲率方程，而是引入“例外点”概念（Definition 2.1），将几何问题转化为速度场与曲线的正交性判定。这一步骤的关键在于发现：若 $$(U,V)$$ 法向于 $\Gamma_{\text{sonic}}$ ，则由伯努利律 (3.2) 沿切向求导可得 $c \partial_\tau c + \frac{\gamma-1}{\gamma+1}(U + f' V) = 0$ ，而法向性即 $$U + f' V = 0$$ ，从而强制 $\partial_\tau c = 0$ ，即 $$c$$ 为常数。此代数关系绕开了对高阶导数的直接估计，成为后续所有定理的基石。

2. 辅助函数 $\psi_Q$ 的构造与精细二阶展开：为在无先验法向性假设下证明 $$Q$$ 必为例外点，作者定义 $\psi_Q = \varphi - \varphi^Q_2$ ，其中 $\varphi^Q_2$ 是以 $$Q$$ 为中心、匹配 $\varphi(Q)$ 与 $D\varphi(Q)$ 的二次多项式（(2.17)）。此构造的精妙之处在于： $\psi_Q$ 在 $$Q$$ 处自动满足 $\psi_Q(Q) = 0$ , $D\psi_Q(Q) = 0$ ，其 Hessian 矩阵 $D^2 \psi_Q(Q)$ 编码了 $\varphi$ 在 $$Q$$ 处的真实二阶信息。定理 2.5 的证明依赖于对 $D^2 \psi_Q(Q)$ 的两侧极限（来自 $\Omega_-$ 与 $\Omega_+$ ）进行比较，利用原方程 (2.4) 的结构，导出 $\partial_\tau c(Q)$ 与 $\partial_n c(Q)$ 的显式表达式，最终将 $\psi_Q$ 的二次增长速率与 $\frac{n_Q \cdot (U,V)}{c}$ 的平方联系起来。

3. 障碍函数法与角度控制：定理 2.6 的证明面临更大挑战：需在 $\theta_1$ 很小时证明法向性。作者构造了一个精心设计的障碍函数（(7.12)），其形式为 $w = a_2 ((\xi - \xi_*) \cos \sigma_* + (\eta - \eta_*) \sin \sigma_*)^2 - \epsilon ((\eta - \eta_*) \cos \sigma_* - (\xi - \xi_*) \sin \sigma_*)^2$ ，该函数在 $\Omega_+$ 内满足 $\mathcal{L} w \geq 0$ （其中 $\mathcal{L}$ 为线性化算子），并与 $\psi_Q$ 进行比较。通过矛盾法，若 $\theta_1$ 过大，则 $$w$$ 无法在 $\Gamma_{\text{sonic}}$ 附近控制 $\psi_Q$ ，违反最大值原理，从而迫使 $\theta_1$ 必须足够小，即 $$(U,V)$$ 必须接近法向。

4. 局部坐标系下的曲率 ODE 求解：在定理 2.1 的证明中，作者将几何约束 $$U + f' V = 0$$ 与运动学约束 $$U^2 + V^2 = c_0^2$$ 结合，代入方程 (3.6)，最终导出关于 $f'(\xi)$ 的常微分方程 (3.11)： $\frac{F'(\xi)}{(1+F^2(\xi))^{3/2}} = \pm \frac{1}{c_0}$ 。此方程的积分直接给出 $f(\xi)$ 满足 $(\xi \mp a c_0)^2 + (f(\xi) - b)^2 = c_0^2$ ，从而严格证明局部圆性，并将圆心坐标 $(\pm a c_0, b)$ 显式表示为初始速度分量 $$U(0,f(0)), V(0,f(0))$$ 的函数。这一解析求解避免了泛泛的“曲率恒定”论证，提供了精确的几何参数映射。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度：

1. 几何假设的彻底解放：此前所有工作（[1,2,5–7,9,11,12,16,17]）均将 $\Gamma_{\text{sonic}}$ 预设为圆弧（或直线），其半径与中心由全局常数 $$c_0, (u_2,v_2)$$ 确定。本文首次在无任何先验几何假设（仅 $C^{1,1}$ 正则性）下，证明了在小扰动框架内 $\Gamma_{\text{sonic}}$ 必为圆弧，且其几何参数由局部解值 $$(c(Q), U(Q), V(Q))$$ 决定。这是从“验证已知几何”到“反演未知几何”的范式转变。

2. 退化分析的尺度深化：以往工作（如 [1,9]）聚焦于在已知圆 $\Gamma_{\text{sonic}}$ 上建立 $C^{2,\alpha}$ 正则性，其核心是克服退化带来的系数消失，采用加权 Schauder 估计或 Morrey 空间方法。本文则将分析尺度推进至退化边界本身：通过 $\psi_Q$ 的二次接触分析，揭示了退化不仅影响解的光滑性，更严格约束了边界自身的曲率——任何偏离圆弧的扰动都会在 $\psi_Q$ 的二阶项中产生不可忽略的误差，从而被定理 2.5/2.6 排除。这表明退化不仅是正则性障碍，更是几何选择器。

3. 适用范围的实质性扩展：[1,9] 的结果仅适用于严格均匀入流；本文将结论推广至 $$C^2$$ 小扰动情形，其扰动幅度 $\delta$ 受限于 $\epsilon_0, \varepsilon_1$ （依赖于 $\|f\|_{C^{1,1}}$ , $\|\psi_Q\|_{C^{1,\alpha}}$ 等），但未要求扰动具有特定频率或对称性。这意味着只要入流扰动足够小，无论其空间分布如何（如局部峰、振荡）， $\Gamma_{\text{sonic}}$ 的全局圆结构仍被保持。代价是：结论仅对“小扰动”成立，且正则性结论 $C^{2,\alpha}$ 仅在伪亚音速区成立（双曲区 $\Omega_-$ 的正则性未被讨论）。

4. 技术工具的原创性跃迁：既有文献主要依赖特征线方法（[21,22,24]）、能量估计（[15]）或变分方法（[25–27]）。本文则开创性地引入 $\psi_Q$ 函数与例外点判据，将几何问题转化为局部二次形式的比较问题。这种“几何-代数-分析”三重耦合方法，与 [1] 中基于圆对称性的分离变量法、[9] 中的迭代构造法有根本不同。它不依赖全局对称性，因而原则上可推广至三维或多维情形（作者在摘要中明确指出其技术“expected to be applicable to other nonlinear problems involving similar mixed-type degeneracies”）。然而，本文未处理 $\Gamma_{\text{sonic}}$ 的拓扑变化（如分裂、闭合）或大扰动情形，这些仍是开放问题。

8. 旋转磁流体动力学与旋转流体系统3D-2D渐近性的改进收敛速率及更大的初值数据

Improved convergence rates and larger initial data for the 3D-2D asymptotics of the rotating magnetohydrodynamicand rotating fluids systems

作者: Frédéric Charve

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22944v1

摘要翻译: 我们最近证明：若其初始速度场与磁场均可分解为经典三维部分与某个二维部分（即仅依赖于水平空间变量）之和，则三维旋转磁流体动力学（MHD）系统的解（当Rossby数趋于零时）收敛于一个六分量二维MHD系统与一个由二维极限速度所输送的附加三维磁场的解。我们还针对大尺度非适定初值数据给出了显式的收敛速率。本文中，借助新的估计，我们改进了收敛速率，并允许更大的初值数据。最后，这些改进被推广至旋转流体系统。

Comment: 无

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

1. 研究对象和背景

本文研究三维旋转磁流体（rotating magnetohydrodynamic, MHD）系统在Rossby数 $\varepsilon \to 0$ 极限下的强解渐近行为，其控制方程为：

\begin{cases} \partial_t u^\varepsilon - \nu \Delta u^\varepsilon + \dfrac{1}{\varepsilon} u^\varepsilon \wedge e_3 + u^\varepsilon \cdot \nabla u^\varepsilon - b^\varepsilon \cdot \nabla b^\varepsilon = -\nabla p^\varepsilon, \\ \partial_t b^\varepsilon - \nu' \Delta b^\varepsilon + u^\varepsilon \cdot \nabla b^\varepsilon - b^\varepsilon \cdot \nabla u^\varepsilon = 0, \\ \operatorname{div} u^\varepsilon = 0, \\ (u^\varepsilon, b^\varepsilon)|_{t=0} = (u_{0,\varepsilon}, b_{0,\varepsilon}). \end{cases} \tag{MHD_\varepsilon}

其中 $u^\varepsilon = (u_1^\varepsilon, u_2^\varepsilon, u_3^\varepsilon)$ 为速度场， $b^\varepsilon = (b_1^\varepsilon, b_2^\varepsilon, b_3^\varepsilon)$ 为磁场， $$e_3 = (0,0,1)$$ ， $\nu > 0$ 、 $\nu' > 0$ 分别为运动粘性系数与磁扩散率。该模型描述地球物理与天体物理中受强旋转影响的导电流体（如液态外核、太阳对流层），其数学核心挑战在于：当 $\varepsilon \to 0$ 时，科里奥利项 $\frac{1}{\varepsilon} u^\varepsilon \wedge e_3$ 产生高频振荡，导致解在标准 Sobolev 空间中不紧致，传统能量方法失效；同时，三维初值与二维极限结构之间的尺度失配引发“ill-prepared”初值问题——即初值不能被完全投影到二维子空间，其三维分量 $v_{0,\varepsilon}, c_{0,\varepsilon}$ 需以可控方式衰减或增长。

作者此前工作（[11]）已确立如下关键结构：若初值形如

(u_{0,\varepsilon}, b_{0,\varepsilon}) = (\tilde{u}_0(x_h) + v_{0,\varepsilon}(x), \tilde{b}_0(x_h) + c_{0,\varepsilon}(x)), \quad x_h = (x_1,x_2), \tag{1.1}

其中 $\tilde{u}_0, \tilde{b}_0 \in H^\delta(\mathbb{R}^2)$ 为二维主部， $(v_{0,\varepsilon}, c_{0,\varepsilon})$ 为三维扰动，则解 $U^\varepsilon = (u^\varepsilon, b^\varepsilon)$ 弱收敛至由两部分组成的极限：(i) 六分量二维 MHD 系统 $(\tilde{u}, \tilde{b})$ ，仅依赖水平变量 $$x_h$$ ，满足

\begin{cases} \partial_t \tilde{u} - \nu \Delta_h \tilde{u} + \tilde{u}_h \cdot \nabla_h \tilde{u} - \tilde{b}_h \cdot \nabla_h \tilde{b} = -\nabla_h q_0, \\ \partial_t \tilde{b} - \nu' \Delta_h \tilde{b} + \tilde{u}_h \cdot \nabla_h \tilde{b} - \tilde{b}_h \cdot \nabla_h \tilde{u} = 0, \\ \operatorname{div}_h \tilde{u} = 0, \\ (\tilde{u}, \tilde{b})|_{t=0} = (\tilde{u}_0, \tilde{b}_0), \end{cases} \tag{2D-MHD_3}

(ii) 一个额外的三维磁场 $c^\varepsilon$ ，由极限速度 $\tilde{u}$ 输运并满足线性方程

\begin{cases} \partial_t c^\varepsilon - \nu' \Delta c^\varepsilon + \tilde{u} \cdot \nabla c^\varepsilon - c^\varepsilon \cdot \nabla \tilde{u} = 0, \\ c^\varepsilon|_{t=0} = c_{0,\varepsilon}. \end{cases} \tag{3D-M_\varepsilon}

此分解揭示了旋转 MHD 的深层几何结构：二维主部主导大尺度动力学，而三维扰动 $c^\varepsilon$ 在 $\tilde{u}$ 的拉伸下演化，形成“被动标量型”磁场分量。然而，[11] 中的收敛速率仅为 $\varepsilon^{1/18(\delta/2 - \gamma)}$ ，且对初值规模施加严苛限制： $\|v_{0,\varepsilon}\|_{\dot{H}^{1/2 + \delta}} \lesssim \varepsilon^{-5\delta/12}$ 、 $\|c_{0,\varepsilon}\|_{H^{1/2 + \delta}} \lesssim |\ln \varepsilon|^{1/4}$ ，且 $\delta < 1/6$ 。这些限制源于关键非线性项 $c^\varepsilon \cdot \nabla W^\varepsilon$ （记为 $G_{14}$ ）及 $c^\varepsilon \cdot \nabla d^\varepsilon$ （ $G_{12}$ ）、 $b^\varepsilon \cdot \nabla d^\varepsilon$ （ $F_{13}$ ）的能量估计瓶颈——它们无法被经典 Strichartz 或 Sobolev 嵌入充分控制，尤其当 $c^\varepsilon$ 具有对数增长初值时。因此，本文旨在突破这一分析障碍：一方面将收敛速率提升至理论最优阶 $\varepsilon^{k(\delta/2 - \gamma)}$ （ $k \to 1^-$ ），另一方面将初值容许范围扩展至 $\|c_{0,\varepsilon}\|_{H^{1/2 + \delta}} \lesssim |\ln \varepsilon|^{1/2}$ 与 $\|v_{0,\varepsilon}\|_{\dot{H}^{1/2 + \delta}} \lesssim \varepsilon^{-\delta/2}$ ，从而逼近物理上更自然的“大初值”情形。

2. 主要定理或主要结果

本文在两类初值设定下建立了严格收敛结果：

第一类结果（Theorem 2.5）：对任意 $C_0 \geq 1$ 、 $\delta \in ]0, 1/4[$ 、 $k \in [0,1[$ （可任意接近1）、 $\gamma \in [0, \delta/2[$ ，存在 $\varepsilon_0, K_0 \in ]0,1]$ ，使得对任意 $\varepsilon \in ]0, \varepsilon_0]$ 及初值满足 (1.1) 式，其中 $\tilde{u}_0, \tilde{b}_0 \in H^\delta(\mathbb{R}^2)$ 、 $(v_{0,\varepsilon}, c_{0,\varepsilon}) \in (\dot{H}^{1/2 - \delta} \cap \dot{H}^{1/2 + \delta})(\mathbb{R}^3) \times H^{1/2 + \delta}(\mathbb{R}^3)$ 且

\|\tilde{u}_0\|_{H^\delta}, \|\tilde{b}_0\|_{H^\delta} \leq C_0, \quad \|v_{0,\varepsilon}\|_{\dot{H}^{1/2 + \delta}} \leq C_0 \varepsilon^{-\gamma}, \quad \|c_{0,\varepsilon}\|_{H^{1/2 + \delta}} \leq K_0 |\ln \varepsilon|^{1/2},

则存在唯一全局强解 $U^\varepsilon = (u^\varepsilon, b^\varepsilon)$ ，且对任意 $s \in [1/2, 1/2 + l_0(\delta/2 - \gamma)]$ （ $$l_0 > 0$$ 依赖于参数），有

\|U^\varepsilon - (\tilde{u} + W^\varepsilon, \tilde{b} + c^\varepsilon)\|_{\dot{F}^s} = \|(\delta^\varepsilon, d^\varepsilon)\|_{\dot{F}^s} \leq F_0 \varepsilon^{k(\delta/2 - \gamma)}.

进一步，对任意 $p \in ]2, 2/\delta[$ ，存在常数 $$F'_0 > 0$$ 使得

\|U^\varepsilon - (\tilde{u}, \tilde{b} + c^\varepsilon)\|_{L^p L^{3/(1-2/p)}} \leq F'_0 \varepsilon^{k(\delta/2 - \gamma)}.

若 $v_{0,\varepsilon} \in \dot{H}^{1/2 - \delta} \cap \dot{H}^{1/2 + \delta}$ ，则上述估计对 $s \in [1/2 - l_0(\delta/2 - \gamma), 1/2 + l_0(\delta/2 - \gamma)]$ 成立，并可取 $$p = 2$$ 得

\|U^\varepsilon - (\tilde{u}, \tilde{b} + c^\varepsilon)\|_{L^2 L^\infty} \leq F_0 \varepsilon^{k(\delta/2 - \gamma)}.

第二类结果（Theorem 2.6）：对任意 $C_0 \geq 1$ 、 $\delta \in ]0, 1/6]$ 、 $\mu \in [0,1[$ ，存在 $m_0, \varepsilon_0, F_0 > 0$ ，使得对任意 $\varepsilon \in ]0, \varepsilon_0]$ 及初值满足 (1.1) 式，其中 $\tilde{u}_0, \tilde{b}_0 \in H^\delta(\mathbb{R}^2)$ 、 $(v_{0,\varepsilon}, c_{0,\varepsilon}) \in (\dot{H}^{1/2} \cap \dot{H}^{1/2 + \delta})(\mathbb{R}^3) \times H^{1/2 + \delta}(\mathbb{R}^3)$ 且

\|\tilde{u}_0\|_{H^\delta}, \|\tilde{b}_0\|_{H^\delta} \leq C_0, \quad \|v_{0,\varepsilon}\|_{\dot{H}^{1/2 + \delta}} \leq m_0 \varepsilon^{-\delta/2}, \quad \|c_{0,\varepsilon}\|_{H^{1/2 + \delta}} \leq C_0,

则存在唯一全局强解，且

\|U^\varepsilon - (\tilde{u} + W^\varepsilon, \tilde{b} + c^\varepsilon)\|_{\dot{F}^{1/2}} \leq F_0 \max(m_0, \varepsilon^{\delta/2}).

相应地，对任意 $p \in ]2, 2/\delta[$ ，有

\|U^\varepsilon - (\tilde{u}, \tilde{b} + c^\varepsilon)\|_{L^p L^{3/(1-2/p)}} \leq F_0 \max(m_0, \varepsilon^{\delta/2}).

该结果表明：当 $$m_0$$ 被替换为趋于零的函数 $m(\varepsilon)$ 时，仍保证收敛性，从而允许初值 $v_{0,\varepsilon}$ 的 $\dot{H}^{1/2 + \delta}$ 范数达 $\varepsilon^{-\delta/2}$ 量级，远超 [11] 中的 $\varepsilon^{-5\delta/12}$ 。

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

证明的核心障碍在于处理系统 (2.11) 中的非线性项 $G_{14} = c^\varepsilon \cdot \nabla W^\varepsilon$ 、 $F_{13} = c^\varepsilon \cdot \nabla d^\varepsilon$ 与 $G_{12} = c^\varepsilon \cdot \nabla \delta^\varepsilon$ 。这些项在 $\dot{H}^s$ 内积估计中产生形如 $|(c^\varepsilon \cdot \nabla W^\varepsilon | d^\varepsilon)_{\dot{H}^s}|$ 的表达式，经典方法（如 $$L^2$$ -Hölder 或 $$L^3$$ - $$L^6$$ 乘积估计）会引入 $\|c^\varepsilon\|_{\dot{H}^{s+1/2}} \|W^\varepsilon\|_{L^6}$ ，但 $\|W^\varepsilon\|_{L^6}$ 的 Strichartz 衰减阶不足以抵消 $c^\varepsilon$ 的对数增长或 $\varepsilon^{-\delta/2}$ 增长。本文的真正创新点是系统性引入非局部算子 $$|D|^s$$ （ $s \in ]0,1[$ ）及其双线性分解（Proposition 2.8）：

$$ |D|^s(fg) = (|D|^s f) g + f (|D|^s g) + M_s(f,g), $$

其中余项 $M_s(f,g)(x) = \int_{\mathbb{R}^3} \frac{(f(x)-f(x-y))(g(x)-g(x-y))}{|y|^{3+s}} dy$ 满足精确定量的 Besov 型估计 $\|M_s(f,g)\|_{L^r} \lesssim \|f\|_{\dot{B}^{s_1}_{r_1,q_1}} \|g\|_{\dot{B}^{s_2}_{r_2,q_2}}$ （ $$s_1+s_2=s$$ , $$1/r = 1/r_1 + 1/r_2$$ ）。对 $G_{14}$ ，作者将内积重写为

|(c^\varepsilon \cdot \nabla W^\varepsilon | d^\varepsilon)_{\dot{H}^s}| = |(\operatorname{div}(c^\varepsilon \otimes W^\varepsilon) | d^\varepsilon)_{\dot{H}^s}| = |(c^\varepsilon \otimes W^\varepsilon | \nabla d^\varepsilon)_{\dot{H}^s}| = |(|D|^s(c^\varepsilon \otimes W^\varepsilon) | |D|^s \nabla d^\varepsilon)_{L^2}|,

再应用分解，得到三项贡献： $(|D|^s c^\varepsilon) W^\varepsilon$ 、 $c^\varepsilon (|D|^s W^\varepsilon)$ 和 $M_s(c^\varepsilon, W^\varepsilon)$ 。关键在于，通过精细选择 $\beta, \beta'$ （如 $\beta = \beta' = 1/4$ ），利用 Sobolev 嵌入 $\dot{H}^{s+1/2} \hookrightarrow L^{3/(1-s)}$ 与 Bernstein 不等式，将每项控制为

\lesssim \|c^\varepsilon\|_{\dot{H}^{s+1/2}} \left( \|W^\varepsilon\|_{L^6} + \||D|^s W^\varepsilon\|_{L^{3/(s+1/2)}} + \|W^\varepsilon\|_{\dot{B}^{1/4}_{4,2}} \right) \|\nabla d^\varepsilon\|_{\dot{H}^s}.

随后，通过 Young 不等式吸收 $\|\nabla d^\varepsilon\|_{\dot{H}^s}^2$ 项，并将剩余因子代入 Strichartz 估计（Proposition 2.11）： $\|W^\varepsilon\|_{L^4 L^6} \lesssim \varepsilon^{\eta_0 \delta (1-l)}$ 、 $\||D|^s W^\varepsilon\|_{L^4 L^{3/(s+1/2)}} \lesssim \varepsilon^{\eta_0 \delta (1-l)}$ 、 $\|W^\varepsilon\|_{\dot{B}^{1/4}_{4,2}} \lesssim \varepsilon^{\eta_0 \delta (1-l)}$ ，结合 $\|c^\varepsilon\|_{\dot{H}^{s+1/2}} \lesssim |\ln \varepsilon|^{1/2}$ ，最终获得 $G_{14}$ 的控制阶为 $\varepsilon^{2\eta_0 \delta (1-l)} |\ln \varepsilon|$ ，这足以支撑 $k \to 1^-$ 的收敛速率。类似地，对 $F_{13}$ 和 $G_{12}$ ，作者采用相同框架，但利用 $c^\varepsilon$ 的更高正则性（ $\dot{H}^{3/2}$ 或 $\dot{B}^{3/2}_{2,1}$ ）与 $\delta^\varepsilon$ 、 $d^\varepsilon$ 的 $\dot{H}^s$ 范数配对，避免出现 $\varepsilon^{-\delta/2}$ 因子，从而将 $c_{0,\varepsilon}$ 的容许规模从 $|\ln \varepsilon|^{1/4}$ 提升至 $|\ln \varepsilon|^{1/2}$ 。此外，为处理 $v_{0,\varepsilon}$ 达 $\varepsilon^{-\delta/2}$ 量级的情形，作者在 Theorem 2.6 的证明中，对 $F_6 = -P(W^\varepsilon \cdot \nabla \tilde{u})$ 、 $F_8 = -P(\tilde{u} \cdot \nabla W^\varepsilon)$ 等项实施更精细的时空积分估计，利用 Proposition 2.3 中 $\tilde{u}$ 的 $L^p \dot{H}^\sigma$ 可积性（ $p \in [2, 2/(\sigma - \sigma_+)]$ ）与 $W^\varepsilon$ 的 Strichartz 范数相乘，确保所有项的总和仍可被 $\max(m_0, \varepsilon^{\delta/2})$ 控制。

4. 与之前工作的比较

本文相较前作 [11] 实现了三重实质性推进：

收敛速率维度：[11] 中定理 1.3 给出的收敛阶为 $\varepsilon^{1/18(\delta/2 - \gamma)}$ ，受限于 $G_{14}$ 项的粗略估计（使用 $$L^2$$ - $L^\infty$ 乘积），其 $\gamma$ 范围被压缩至 $[0, 5\delta/12]$ 。本文通过非局部算子分解，将 $G_{14}$ 的控制精度提升至 $\varepsilon^{2\eta_0 \delta (1-l)} |\ln \varepsilon|$ ，从而将收敛阶优化为 $\varepsilon^{k(\delta/2 - \gamma)}$ （ $k \to 1^-$ ），逼近理论预期的最优速率。这一改进直接源于对振荡解 $W^\varepsilon$ 的 Strichartz 范数（如 $\|W^\varepsilon\|_{L^4 L^6}$ ）与 $c^\varepsilon$ 的对数增长之间耦合关系的精确量化。

初值规模维度：[11] 要求 $\|c_{0,\varepsilon}\|_{H^{1/2 + \delta}} \lesssim |\ln \varepsilon|^{1/4}$ 且 $\delta < 1/6$ ，而本文将前者放宽至 $|\ln \varepsilon|^{1/2}$ 并将 $\delta$ 上界延拓至 $$1/4$$ 。这一突破的关键在于 $F_{13}$ 与 $G_{12}$ 的新估计：[11] 中这两项的控制依赖于 $\|c^\varepsilon\|_{\dot{H}^{3/2}}$ ，而 $\|c^\varepsilon\|_{\dot{H}^{3/2}}$ 的增长由 (2.5) 式中的指数因子 $\exp(C \|c_{0,\varepsilon}\|_{\dot{H}^{1/2}}^2)$ 主导，导致 $|\ln \varepsilon|^{1/4}$ 成为临界阈值；本文则通过 $$M_s$$ 项的 Besov 估计，将控制降阶至 $\|c^\varepsilon\|_{\dot{B}^{3/2}_{2,1}}$ ，其增长被 $\|c_{0,\varepsilon}\|_{H^{1/2 + \delta}}$ 线性主导，故 $|\ln \varepsilon|^{1/2}$ 成为新阈值。同样， $v_{0,\varepsilon}$ 的规模从 [11] 的 $\varepsilon^{-5\delta/12}$ 提升至 $\varepsilon^{-\delta/2}$ ，这得益于对 $F_6, F_8, G_9, G_{13}$ 等含 $\tilde{u}$ 的项采用 $\tilde{u} \in L^p \dot{H}^\sigma$ （Proposition 2.3）与 $W^\varepsilon$ 的 Strichartz 范数的匹配估计，而非依赖 $\tilde{u}$ 的 $L^\infty$ 范数（后者无时间衰减）。

中文标题: 一类一维拟线性波动方程族在变分情形之外的奇点形成

摘要翻译: 我们考虑如下参数化非线性波动方程的经典解的有限时间奇点形成问题：

中文标题: 具有无限黏弹性记忆的非对称分数阶耦合可磁化压电梁：多项式衰减及其衰减速率的尖锐性

摘要翻译: 本文研究具有无限黏弹性记忆的可磁化压电梁耦合双曲系统在长时间尺度下的动力学行为。记忆耗散由 $$A^α$$ 刻画，其中 $$A$$ 是一个正自伴算子， $α\in[0,1)$ 。采用非对称分数阶磁电耦合：力学–磁学耦合使用整数阶算子 $$A$$ ，而磁反馈至力学则由分数阶算子 $$A^β$$ （ $β\in[0,1)$ ）主导。该模型弥合了已被深入研究的整数阶耦合情形（ $$β=1$$ ）与尚未解决的全分数阶对称耦合问题之间的鸿沟，为相关耦合系统提供了一个普适框架。在记忆核与系统参数满足适当温和假设的前提下，我们借助半群理论证明系统的适定性，并对光滑初值导出显式的多项式衰减估计：

\|X(t)\|_{\mathcal H} \le C t^{-\frac{1}{4-2β-2α}}\|X_0\|_{D(\mathcal A)},\quad \forall\,t\ge 1,

其中衰减指数显式依赖于 $$α$$ 和 $$β$$ 。当记忆核呈指数衰减时，若刚度系数满足 $α_1\neα_2$ ，则该衰减速率是尖锐的；当 $$α_1=α_2$$ 时，仅能得到衰减指数 $$δ$$ 的上界估计 $δ\le \frac{1}{3-β-2α}$ ，最优速率仍属开放问题。与整数阶反馈耦合（ $$β=1$$ ）的对比表明，分数阶反馈（ $$β<1$$ ）会延缓能量衰减，揭示了 $$β$$ 削弱间接阻尼效应并劣化结构稳定性的机制。结果揭示了强耦合耗散系统中分数阶记忆耗散与分数阶耦合之间内在的相互作用。

10. Hagstrom-Warburton 无反射边界条件对一阶时变 Maxwell 方程组的 Kreiss 稳定性分析

Kreiss stability analysis of Hagstrom-Warburton nonreflecting boundary conditions for the first-order time-dependent Maxwell equations

作者: Kaifang Liu, Matthias Schlottbom, J. J. W. van der Vegt, Yan Xu

类别: math.AP, math.NA

发布日期: 2026-06-22

arXiv ID: 2606.22754v1

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

摘要翻译

14. 抛物型Sobolev空间中零延拓的失效

Failure of zero extension in parabolic Sobolev spaces

作者: Jongkeun Choi, Doyoon Kim, Kwan Woo

类别: math.AP

发布日期: 2026-06-20

arXiv ID: 2606.22059v1

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

摘要翻译

中文标题: 抛物型Sobolev空间中零延拓的失效

摘要翻译: 我们证明：在抛物型Sobolev空间 $\mathring{\mathcal{H}}^1_p((0,T) \times Ω)$ 中，沿空间边界的零延拓可能失效；该空间亦可表征为

L_p(0,T;\mathring{W}^1_p(Ω))\cap W^1_p(0,T; W^{-1}_{p}(Ω)).

更精确地，对任意 $p\in [1, \infty)$ ，我们构造一个函数 $u\in \mathring{\mathcal{H}}^1_p((0,T)\times \mathbb{R}^d_+)$ ，其零延拓不属于 $\mathcal{H}^1_p((0,T)\times \mathbb{R}^d)$ 。这一障碍即使在平坦边界情形下亦发生，其根源在于集中于初始-边界角点处的自相似边界层，该边界层在零延拓后产生一个支撑于边界的法向通量缺陷。我们还讨论了各类Sobolev型空间作为散度型抛物方程解空间的适用性。

15. 严格凸障碍物下耗散算子的本征值不存在性

Absence of eigenvalues of dissipative operator for strictly convex obstacles

作者: Vesselin Petkov

类别: math.AP, math-ph

发布日期: 2026-06-20

arXiv ID: 2606.21958v1

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

摘要翻译

中文标题: 严格凸障碍物下耗散算子的本征值不存在性

摘要翻译: 我们研究 $d \geq 3$ （奇数）维欧氏空间 ${\mathbb R}^d$ 中严格凸有界区域 $K \subset {\mathbb R}^d$ 外部的波动方程，其边界 $$Γ$$ 上满足耗散边界条件 $\partial_ν u - γ(x) \partial_t u = 0$ ，且 $0 < γ(x) <1, \:\forall x \in Γ.$ 解由压缩半群 $V(t) = e^{tG}, \: t \geq 0$ 描述。在文献 [10] 中，我们证明了当 $γ\equiv const$ 且 $K = \{x \in {\mathbb R}^3: \:|x| \leq 1\}$ 时，算子 $$G$$ 无本征值，并猜想该结论对任意严格凸障碍物均成立。本文证明了这一猜想。

16. Wasserstein 梯度流的有限维约化与尖锐衰减率

Finite-dimensional reduction of a Wasserstein gradient flow and sharp decay rates

作者: Dohyun Kim, Hansol Park, Woojoo Shim

类别: math.AP

发布日期: 2026-06-20

arXiv ID: 2606.21918v1

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

摘要翻译

中文标题: Wasserstein 梯度流的有限维约化与尖锐衰减率

摘要翻译: 我们研究由一族扩展广义方差泛函所生成的 Wasserstein 梯度流，该泛函定义为单纯形的期望平方 $$n$$ -维体积，其中包含经典方差型相互作用与广义方差作为特例。关键结构观察在于：该泛函仅依赖于协方差矩阵。因此，Wasserstein 梯度流可约化为关于协方差矩阵特征值的有限维系统，且全测度值解可通过一个显式的线性前推表示式恢复。利用该表示式，我们在不假设紧支集条件下，对任意初始数据 $\mu_0 \in \mathcal P_2(\mathbb R^d)$ 建立了整体适定性。我们还研究了该流的长时间行为：对任意初始数据 $\mu_0 \in \mathcal P_2(\mathbb R^d)$ ，解收敛于一个极限平衡测度，其协方差矩阵的秩严格小于 $$n$$ ；此外，我们在所有谱区域中获得了尖锐的收敛速率：在非退化情形下为指数速率，在退化情形下为具有最优指数的代数速率。

不相关论文（基本信息）

1. 切换机制的Volterra算子：模态稳定性与淬火放大

Regime-Switching Volterra Operators: Modal Stability and Quenched Amplification

作者: Mauricio Herrera-Marín

类别: stat.AP, math.AP, math.DS, math.PR

发布日期: 2026-06-22

arXiv ID: 2606.23558v1

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

摘要: We develop an operator-theoretic framework for finite-dimensional, regime-dependent Volterra equations with completely monotone memory kernels, dissipative network coupling, and Hawkes-type self-excitation. For each fixed regime we construct the associated Volterra resolvent family and prove global well-posedness, continuity across regime switches, and explicit a priori bounds. The main stability result is sharp in the commuting case: after simultaneous diagonalization of the network Laplacian and the excitation operator, each mode obeys a scalar characteristic equation, and global asymptotic stability holds exactly when every modal branching ratio lies below the intensity damping threshold. We also give a norm-based sufficient condition for noncommuting operators and a Perron–Frobenius spectral criterion for nonnegative intensity blocks, showing when norm estimates are conservative. Beyond mean stability, we prove a pathwise finite-range power law for burst amplitudes generated by residence in a Hurwitz but nonnormal regime: under a cone-alignment event, the survival exponent is the ratio of the regime exit rate to a cone-corrected finite-time growth rate bounded above by the logarithmic norm of a fixed Markovian realization in the chosen Euclidean metric. A complementary idealized-feedback result shows how a logarithmic-norm contraction caps the amplification band. Finally, we derive the deterministic intensity block as a mean-field limit of a relaxing long-memory Hawkes system with regimes. Numerical experiments on modal equations, a small-world network, and a switched nonnormal ODE validate the sharp threshold and the finite-range amplification mechanism without using the closed-form tail formula as input.

2. 黎曼流形正质量定理在所有维数下低余维奇点存在时的成立

Riemannian Positive Mass Theorem in All Dimensions in the Presence of Low-Codimension Singularities

作者: Marcus Khuri, Jian Wang, Jinmin Wang

类别: math.DG, gr-qc, math.AP

发布日期: 2026-06-22

arXiv ID: 2606.23529v1

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

摘要: We prove the Riemannian positive mass theorem in all dimensions for asymptotically flat $L^\infty$ -metrics with subcritical singular sets. More precisely, we consider complete asymptotically flat manifolds whose metrics are smooth away from a compact singular set of Minkowski dimension less than $n-3+\frac{2}{n}$ , and whose scalar curvature is nonnegative on the regular set. We show that the ADM mass of each asymptotically flat end is nonnegative, and that the mass vanishes in some end only in the Euclidean case. For the rigidity statement, we require additionally that the Minkowski dimension of the singular set is not larger than $n-3+\frac{1}{n-1}$ . This gives an asymptotically flat analogue of Schoen’s codimension-three conjecture for positive scalar curvature. The proof combines a density theorem for singular asymptotically flat metrics, capacity estimates across the singular set, conformal blow-up inspired by Bi-Hao-He-Shi-Zhu [3], and a $$μ$$ -bubble dimension-descent argument adapted from Brendle-Wang [6].

3. 分布意义下数量曲率有下界的流形：最优正则性与正质量

Spaces with distributional scalar curvature bounded from below: Optimal regularity and positive mass

作者: Man-Chun Lee, Florian Litzinger, Miles Simon

类别: math.DG, math.AP

发布日期: 2026-06-22

arXiv ID: 2606.23272v1

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

摘要: In this work, we study the positive mass theorem under critical low regularity assumptions using Ricci flow smoothing. We show that asymptotically flat manifolds $$(M^n,g)$$ of regularity $L^\infty\cap W^{1,n}$ with non-negative distributional scalar curvature have non-negative ADM mass. Furthermore, when the ADM mass vanishes, the manifold is globally isometric to Euclidean space with respect to an integral distance introduced by De~Cecco-Palmieri. This extends the recent work of Hafemann to the critical regularity case. Our approach is based on showing that Riemannian metrics of regularity $L^\infty\cap W^{1,n}$ , whose scalar curvature is bounded from below in the distributional sense, admit a Ricci flow smoothing whose scalar curvature is bounded from below by the same initial lower bound in the classical sense. In contrast, Cecchini-Frenck-Zeidler constructed examples of metrics which are in $L^\infty\cap W^{1,p}$ for all $$2<p<n$$ , and whose distributional scalar curvature is bounded from below, that cannot be approximated by smooth metrics with the same scalar curvature lower bound. In this sense, our result is optimal.

4. 光滑度量测度空间的质量型不变量及其与分数阶Yamabe问题的关系

A mass-type invariant for smooth metric measure spaces and its relation with the fractional Yamabe problem

作者: Sergio Almaraz, Levi Lopes de Lima, Shaodong Wang

类别: math.DG, math.AP

发布日期: 2026-06-22

arXiv ID: 2606.23248v1

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

摘要: On an asymptotically Schwarzschild smooth metric measure space modelled on the Euclidean half-space, we define a mass-type quantity similar to the classical ADM mass in General Relativity and show its geometric invariance properties. Such quantity has a close relation with the fractional Yamabe problem and the relevant Green’s function.

5. 分形神经算子：通过素谐波Weierstrass编码克服混沌吸引子中的谱偏差

The Fractal Neural Operator: Overcoming Spectral Bias in Chaotic Attractors via Prime-Harmonic Weierstrass Encodings

作者: Kanishk Awadhiya

类别: cs.LG, math.AP

发布日期: 2026-06-22

arXiv ID: 2606.23123v1

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

摘要: Deep learning models, particularly Transformers and Neural Operators, exhibit a well-documented “spectral bias,” effectively acting as low-pass filters that smooth out high-frequency information. While benign in fluid dynamics, this bias is catastrophic for Chaotic Dynamical Systems, where the underlying strange attractor is characterized by fractal geometry and infinite spectral density. We introduce the Fractal Neural Operator (FNO), a novel architecture that utilizes a non-resonant prime number basis to approximate continuous dynamical systems. Unlike geometric encodings ( $$2^k$$ ), which suffer from spectral gaps and resonance, our Harmonic Weierstrass Encoder injects infinite spectral resolution into the latent space. We demonstrate that FNO extends the valid prediction horizon of the Lorenz-63 system to 347 Lyapunov times, exceeding state-of-the-art Reservoir Computing baselines by a factor of 2.3x. These results suggest that “chaos” is not inherently unpredictable to neural networks, but rather requires non-differentiable, fractal embedding manifolds.

6. 概率测度空间上的Mortensen观测器

The Mortensen observer on the space of probability measures

作者: Martin Morange

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

发布日期: 2026-06-22

arXiv ID: 2606.23034v1

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

摘要: We study a deterministic filtering problem formulated directly on the Wasserstein space of probability measures with finite second moment. Motivated by the Mortensen minimum-energy observer, we consider the reconstruction of an evolving probability density from partial observations by minimizing an action functional combining a kinetic transport cost and a time-dependent observation mismatch. The resulting value function is defined on the infinite-dimensional manifold $$(P_2(R^d), W_2)$$ and satisfies a Hamilton-Jacobi-Bellman equation involving the Wasserstein gradient. Under suitable regularity and growth assumptions on the observation functional, we establish dynamic programming principles, continuity of the value function, existence of minimizing trajectories, and viscosity solution properties of the associated Hamilton-Jacobi equation. We provide two complementary notions of viscosity solutions: a geometric formulation based on subdifferentials in Wasserstein space, and a Hilbertian formulation inspired by Lions’ lifting approach. This allows us to prove a comparison principle and uniqueness of solutions. Extensions to transport equations with drift are also discussed. Finally, we introduce a semi-Lagrangian scheme in order to approximate the value function, and show $$Γ$$ -convergence of the scheme.

7. 一类黎曼流形上 $$p$$ -拉普拉斯方程的诺伊曼问题的比较结果

Comparison Results for a class of Neumann Problems of the $$p$$ -Laplace Equation on Riemannian Manifolds

作者: Wentao Liu, Anqiang Zhu

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22854v1

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

摘要: We consider Neumann boundary value problems for the $$p$$ -Laplace equation on Riemannian manifolds with nonnegative Ricci curvature. Using spherical symmetrization under appropriate constraints, we derive Talenti-type comparison results in Lorentz spaces. We further show that, in contrast to the Robin case, the Neumann setting admits weaker constraints, which yields stronger comparison principles.

8. 凸区域中 Hessian 特征值的 Brunn–Minkowski 不等式

A Brunn–Minkowski inequality for the Hessian eigenvalue in convex domain

作者: Jiahuan Li, Xi-Nan Ma, Paolo Salani

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22847v1

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

摘要: We use the deformation methods to obtain the strictly log concavity of solution of a class Hessian equation in bounded convex domain in $\mathbb{R}^{n}$ , as an application we get the Brunn–Minkowski inequality for the Hessian eigenvalue and characterize the equality case in bounded strictly convex domain in $\mathbb{R}^{n}$ .

9. 关于有界区域上的 Yamabe 流

On the Yamabe flow in a bounded domain

作者: Fei Fang, Zhong Tan

类别: math.AP

发布日期: 2026-06-22

arXiv ID: 2606.22815v1

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

摘要: This paper investigates the dynamical behaviors of solutions to the Yamabe flow via the modified potential well method. We first establish the local existence and regularity of weak solutions for the flow. Several new results concerning global existence and blowup are obtained by classifying initial data into stable and unstable sets. Specifically, solutions with initial data in the stable set exist globally and extinguish in finite time, whereas those originating from unstable initial data blow up in infinite time. For certain high-energy initial data, we show that the solution decays to zero as time tends to infinity and undergoes finite-time blowup. In addition, we analyze Palais-Smale sequences to reveal the intrinsic relationship between the long-time asymptotic behavior of solutions and steady states. Finally, we derive the Pohozaev identity for the equation and prove the corresponding nonexistence theorem.

10. 偏置电场下薄平面盒中弯核液晶SmAP $_\text{F}$ 相的极向指向矢结构

Polar director structure of SmAP $_\text{F}$ phase of bent-core liquid crystals in thin planar cells with bias electric field

作者: Alec D. Wendland, Xiaodong Yan

类别: math.AP, cond-mat.soft

发布日期: 2026-06-21

arXiv ID: 2606.22668v1

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

摘要: We study the polar director structure in thin planar cells filled with bent-core liquid crystals in the ferroelectric smectic-A phase (SmAP $_\text{F}$ ). We analyze a continuum phenomenological model proposed in the physics literature and present rigorous proofs of the existence and uniqueness of the equilibrium solutions. We further investigate the qualitative properties of nontrivial solutions and examine the effects of a bias electric field, surface anchoring, and cell thickness on the polar director configuration. Our results are consistent with previous experimental and numerical simulations reported in the physics literature. In addition, our analysis reveals new parameter-dependent behaviors supported by our numerical simulations and extends results reported from previous literature.

11. 反应性Maxwell–Stefan膜输运的熵几何与增强迁移率（有限占据）

Entropy Geometry and Augmented Mobility for Reactive Maxwell–Stefan Membrane Transport with Finite Occupancy

作者: E. Yu. Shchetinin, A. A. Shevchuk, S. I. Salpagarov

类别: math.AP

发布日期: 2026-06-21

arXiv ID: 2606.22624v1

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

摘要: We study a reactive Maxwell–Stefan-type membrane transport system under a finite-occupancy constraint with explicit vacancies. The admissible state space is $\mathcal{D} = \{u \in (0,1)^n : \varrho(u) := \sum_{i=1}^n u_i < 1\}$ , where the vacancy fraction $1-\varrho(u)$ represents the local free volume. This bounded-occupancy geometry induces a Boltzmann–Fermi entropy and a global parametrization by entropy variables. The associated mobility is assumed to split into a composition channel, corresponding to redistribution at fixed total occupancy, and a mass channel, corresponding to variation of the filling fraction. The main structural difficulty is that the unaugmented mobility may lose coercivity in the mass channel. We show that a single rank-one augmentation of the form $γ_0 \mathbf{1} \otimes \mathbf{1}$ restores full coercivity while leaving the composition block unchanged. On this basis, we prove four results: a quantitative channel-wise coercivity estimate for the augmented mobility; global existence of entropy weak solutions via an implicit Rothe scheme in entropy variables; a weak–strong stability estimate in relative entropy with uniqueness in the strong class; and convergence of a fully implicit finite-volume approximation that preserves the bounded-occupancy structure and satisfies a discrete entropy inequality.

12. 关于前 $$N$$ 个非零 Steklov 特征值的尖锐定容积乘积不等式

A sharp fixed-volume product inequality for the first $$N$$ nonzero Steklov eigenvalues

作者: Haiqi Zhang, Quanyu Tang, Yanyang Li

类别: math.AP, math.SP

发布日期: 2026-06-21

arXiv ID: 2606.22571v1

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

摘要: We prove a sharp fixed-volume product inequality for the first $$N$$ nonzero Steklov eigenvalues of bounded Lipschitz domains in $\mathbb R^N$ . More precisely, if $N\ge2$ and $Ω\subset\mathbb R^N$ is a bounded Lipschitz domain, then

\prod_{j=1}^N σ_j(Ω)\le \frac{ω_N}{|Ω|},

where $0=σ_0(Ω)<σ_1(Ω)\leσ_2(Ω)\le\cdots$ are the Steklov eigenvalues of $$Ω$$ , and $$ω_N$$ denotes the volume of the unit ball in $\mathbb R^N$ . This extends the convex-domain theorem of Henrot, Philippin, and Safoui to arbitrary bounded Lipschitz domains, and in particular settles the remaining higher-dimensional case of a problem posed by Henrot.

13. Brakke 流的抛物型可 rectifiability

Parabolic rectifiability of the Brakke flow

作者: Yu Tong Liu, Myles Workman

类别: math.DG, math.AP

发布日期: 2026-06-21

arXiv ID: 2606.22441v1

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

摘要: We prove that the support of the canonical space-time measure for a Brakke flow is a parabolic $$(k+2)$$ -rectifiable set. As a consequence, we obtain that at almost all points along the flow, with respect to this canonical space-time measure, there exists a unique, static, planar tangent flow, and that various notions of density for the flow agree at these points. Moreover, following on from our previous work The space-time-Grassmann measure of the Brakke flow', we continue to develop the approach to the Brakke flow as a space-time-Grassmann measure. We prove that the standard notion of convergence for Brakke flows, coming from the compactness theorem of Ilmanen (7.1 of Elliptic regularization and partial regularity for motion by mean curvature’), is equivalent to the convergence of these space-time-Grassmann Radon measures. This gives an alternate notion of varifold convergence to the one exhibited in 7.1(ii) of `Elliptic regularization and partial regularity for motion by mean curvature'.

14. 具有非局部非线性的线性椭圆问题的解的多重性与构造

Multiplicity and solution construction in linear elliptic problems with nonlocal nonlinearities

作者: Chiun-Chang Lee

类别: math.AP

发布日期: 2026-06-21

arXiv ID: 2606.22295v1

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

摘要: This work investigates linear elliptic equations with multiple nonlocal nonlinearities in a bounded domain, which takes the form \begin{equation*} -\nabla\cdot(\od\nabla u) +λfu= η\sum_{j=1}^kh_j\boldsymbol{\mathsf{N}}_j[u]+h_0. \end{equation*} Here $$λ$$ and $$η$$ are positive parameters, and $\boldsymbol{\mathsf{N}}_j[u]$ represents a nonlocal term dependent on the unknown solution~ $$u$$ . All coefficients are defined in the context of a wide range of applications. As such equations generally lack a variational structure, a new approach is developed that combines fixed-point arguments with asymptotic techniques. This method establishes the existence and multiplicity of solutions under specific conditions. Of particular interest is the role of the nonlocal nonlinearities, which lead to diverse structures of solutions. To the best of our knowledge, this property is a new finding not typically observed in related nonlocal elliptic problems.

15. 如何最优地煮溏心蛋：一个双区域热方程的拉普拉斯变换解

How to Cook a Soft-Boiled Egg Optimally: A Laplace-Transform Solution of a Two-Domain Heat Equation

作者: Matthew Lorig

类别: math.AP, math.OC, physics.class-ph

发布日期: 2026-06-20

arXiv ID: 2606.22156v1

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

摘要: We study the problem of cooking the yolk and albumen of a hen’s egg to their respective optimal temperatures of $T_Y^* = 65^\circ$ C and $T_W^* = 85^\circ$ C, subject to the requirement that neither temperature ever exceed its target at any time during cooking, since temporary overshoot still overcooks the egg even if the final reading is correct. We model the egg as a two-domain sphere with distinct thermal diffusivities, and take the Laplace transform of the heat equation in each domain, reducing the problem to a $3 \times 3$ linear system in the transform variable $$s$$ with hyperbolic-trigonometric solutions. The resulting transform is inverted numerically via Talbot’s method and validated against a finite-difference solver. A single boiling phase cannot satisfy the no-overshoot requirement: the thin outer albumen heats far faster than the insulated yolk and necessarily overshoots $$T_W^*$$ before the yolk approaches $$T_Y^*$$ . We show that a three-phase protocol resolves this: a sous-vide pre-soak at exactly $65^\circ$ C (which cannot overshoot since the bath temperature equals the target), a short boil to bring the albumen toward $$T_W^*$$ , and an ice-water bath that arrests the albumen’s residual overshoot while residual heat continues raising the yolk to its target. Optimizing the phase durations gives $$17.26$$ minutes of sous-vide, $$66$$ seconds of boiling, and an ice bath, achieving both targets at $T^* \approx 20.67$ minutes with neither constraint violated at any time. This compares favorably with the periodic protocol of Di Lorenzo et al. (2025), which requires 32 minutes and misses both targets substantially.

16. Wasserstein 空间中凸 Hamilton-Jacobi 方程的定量均匀化

Quantitative homogenization of convex Hamilton-Jacobi equations in the Wasserstein space

作者: Zhiyan Ding, Ibrahim Ekren, Yuxi Han, Antonios Zitridis

类别: math.AP

发布日期: 2026-06-20

arXiv ID: 2606.22103v1

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

摘要: We study a homogenization problem for first-order Hamilton-Jacobi equations in the Wasserstein space with a convex Hamiltonian. We show that the solution $U^\varepsilon$ , which is the value function of a mean field control problem, converges uniformly as $\varepsilon \to 0$ to the solution of a limiting Hamilton-Jacobi equation whose Hamiltonian is obtained through a suitable cell problem. Furthermore, we establish quantitative rates of convergence. Under general assumptions with multiscale dependence, we prove that the rate of convergence is $O(\sqrt{\varepsilon})$ . When the Hamiltonian depends only on the fast variable and the momentum, we establish the sharp convergence rate $O(\varepsilon)$ . To the best of our knowledge, this is the first quantitative convergence result extending the optimal rate for first-order Hamilton-Jacobi equations in finite dimensions to the Wasserstein space. Finally, we show that our analysis extends to dynamic optimal transport problems, where the terminal condition imposes a constraint on the final distribution.

17. 拟线性保守系统中能量给定的连接轨道

Prescribed–Energy Connecting Orbits for Quasilinear Conservative Systems

作者: Renan J. S. Isneri, Piero Montecchiari

类别: math.AP, math.CA

发布日期: 2026-06-20

arXiv ID: 2606.22031v1

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

摘要: We consider quasilinear conservative systems

(φ(|\dot q|)\dot q)'=\nabla V(q), \qquad q\in\R^{N},

with $$Φ$$ -growth kinetic term and potential $V\in C^{1}(\R^{N};\R)$ . Assuming that for some $c\in\R$ the sublevel set $\{V\le c\}$ splits into two disjoint closed subsets $\mathcal V_c^{-}$ and $\mathcal V_c^{+}$ , we prove the existence of trajectories $$q_c$$ with prescribed energy $$-c$$ connecting $\mathcal V_c^{-}$ and $\mathcal V_c^{+}$ , obtained through an energy-constrained variational method. Although the construction yields weak solutions in an Orlicz-Sobolev setting, minimal $$c$$ -connections are shown to be classical $$C^2$$ trajectories satisfying the strong energy identity $E_{q_c}\equiv -c$ . The resulting entire trajectories include heteroclinic, homoclinic, and brake-type orbits. Applications to double-well, Duffing-type, and multiple pendulum systems are discussed.

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

1. 关于二维广义不可压模型#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

2. 关于物理真空中粘性气体恒星演化方程的局部存在性定理#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

3. 静置液滴问题的整体适定性：二维Navier-Stokes流#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

4. 一种经典的椭圆正则性方法用于几乎调和映射及相关系统#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

5. 奇异定常Navier-Stokes流：例子与稳定性#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

6. Navier-Stokes 方程的有限链 CKN-坏尺度计数：标准偏微分方程闭合与典范探测器实现#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

7. 欧拉方程势流自相似解中伪声速曲线的几何结构#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

8. 旋转磁流体动力学与旋转流体系统3D-2D渐近性的改进收敛速率及更大的初值数据#

1. 研究对象和背景#

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

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

4. 与之前工作的比较#

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

1. 2\times 2 非线性双曲系统解的强迹#

摘要翻译#

2. 图上的分数阶多孔介质方程#

摘要翻译#

3. 一类一维拟线性波动方程族在变分情形之外的奇点形成#

摘要翻译#

4. 各向异性二维分形不确定性原理与量子开放面包师映射#

摘要翻译#

5. 对数Sobolev不等式中尖锐稳定性常数的最优值存在性#

摘要翻译#

6. 与球Banach函数空间相关的Sobolev空间的紧性及其应用#

摘要翻译#

7. 具有 Mercer 核的空间凝结系统：复制子动力学与凝胶化#

摘要翻译#

8. Schrödinger–Intermediate Long Wave 系统的适定性#

摘要翻译#

9. 具有无限黏弹性记忆的非对称分数阶耦合可磁化压电梁：多项式衰减及其衰减速率的尖锐性#

摘要翻译#

10. Hagstrom-Warburton 无反射边界条件对一阶时变 Maxwell 方程组的 Kreiss 稳定性分析#

摘要翻译#

11. 非线性Dirac方程在非齐次非线性项下的归一化解#

摘要翻译#

12. 球面上不可数多族非旋转对称的II型古代Yamabe流#

摘要翻译#

13. 球面上二次导数非线性波动方程的稳定爆破#

摘要翻译#

14. 抛物型Sobolev空间中零延拓的失效#

摘要翻译#

15. 严格凸障碍物下耗散算子的本征值不存在性#

摘要翻译#

16. Wasserstein 梯度流的有限维约化与尖锐衰减率#

摘要翻译#

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

1. 切换机制的Volterra算子：模态稳定性与淬火放大#

2. 黎曼流形正质量定理在所有维数下低余维奇点存在时的成立#

3. 分布意义下数量曲率有下界的流形：最优正则性与正质量#

4. 光滑度量测度空间的质量型不变量及其与分数阶Yamabe问题的关系#

5. 分形神经算子：通过素谐波Weierstrass编码克服混沌吸引子中的谱偏差#

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

1. 关于二维广义不可压模型

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 关于物理真空中粘性气体恒星演化方程的局部存在性定理

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 静置液滴问题的整体适定性：二维Navier-Stokes流

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 一种经典的椭圆正则性方法用于几乎调和映射及相关系统

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 奇异定常Navier-Stokes流：例子与稳定性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. Navier-Stokes 方程的有限链 CKN-坏尺度计数：标准偏微分方程闭合与典范探测器实现

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 欧拉方程势流自相似解中伪声速曲线的几何结构

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 旋转磁流体动力学与旋转流体系统3D-2D渐近性的改进收敛速率及更大的初值数据

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. 2\times 2 非线性双曲系统解的强迹

摘要翻译

2. 图上的分数阶多孔介质方程

摘要翻译

3. 一类一维拟线性波动方程族在变分情形之外的奇点形成

摘要翻译

4. 各向异性二维分形不确定性原理与量子开放面包师映射

摘要翻译

5. 对数Sobolev不等式中尖锐稳定性常数的最优值存在性

摘要翻译

6. 与球Banach函数空间相关的Sobolev空间的紧性及其应用

摘要翻译

7. 具有 Mercer 核的空间凝结系统：复制子动力学与凝胶化

摘要翻译

8. Schrödinger–Intermediate Long Wave 系统的适定性

摘要翻译

9. 具有无限黏弹性记忆的非对称分数阶耦合可磁化压电梁：多项式衰减及其衰减速率的尖锐性

摘要翻译

10. Hagstrom-Warburton 无反射边界条件对一阶时变 Maxwell 方程组的 Kreiss 稳定性分析

摘要翻译

11. 非线性Dirac方程在非齐次非线性项下的归一化解

摘要翻译

12. 球面上不可数多族非旋转对称的II型古代Yamabe流

摘要翻译

13. 球面上二次导数非线性波动方程的稳定爆破

摘要翻译

14. 抛物型Sobolev空间中零延拓的失效

摘要翻译

15. 严格凸障碍物下耗散算子的本征值不存在性

摘要翻译

16. Wasserstein 梯度流的有限维约化与尖锐衰减率

摘要翻译

不相关论文（基本信息）

1. 切换机制的Volterra算子：模态稳定性与淬火放大

2. 黎曼流形正质量定理在所有维数下低余维奇点存在时的成立

3. 分布意义下数量曲率有下界的流形：最优正则性与正质量

4. 光滑度量测度空间的质量型不变量及其与分数阶Yamabe问题的关系

5. 分形神经算子：通过素谐波Weierstrass编码克服混沌吸引子中的谱偏差