生成时间: 2026年07月16日 05:30:37

运行状态: 完成

计划处理论文数量: 34

已完成论文数量: 34

未完成论文数量: 0

重点关注论文数量: 9

了解领域论文数量: 16

不相关论文数量: 9

重点关注论文(完整分析)

1. 关于不可压缩磁流体动力学方程的粘性消失与磁扩散消失极限

On the Vanishing Viscosity and Magnetic Diffusion Limit for Incompressible Magnetohydrodynamics

作者: Sarka Necasova, Tong Tang, Yuanjing Tang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13813v1

摘要翻译: 本文研究有界区域中不可压缩磁流体动力学(MHD)方程组的无粘且无电阻极限。在内部一致正则性和边界处一致等连续性的假设下,我们推导出 Leray-Hopf 解序列在 $μ,ν\rightarrow0$ 时收敛于理想 MHD 系统弱解的充分条件。

Comment: 无

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

1. 研究对象和背景

本文研究三维有界光滑区域 $\Omega \subset \mathbb{R}^3$ 上不可压磁流体动力学(MHD)方程组在粘性系数 $\mu > 0$ 与磁扩散系数 $\nu > 0$ 同时趋于零时的奇性极限问题。核心对象是如下带耗散的 MHD 系统:

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u - \mu \Delta u + \nabla P = (B \cdot \nabla) B, \\ \partial_t B + (u \cdot \nabla) B - \nu \Delta B = (B \cdot \nabla) u, \\ \nabla \cdot u = \nabla \cdot B = 0, \\ (u, B)|_{t=0} = (u_0, B_0), \end{cases} \tag{1.1} $$

配以边界条件:$u|_{\partial \Omega} = 0$, $B \cdot n|_{\partial \Omega} = 0$, $(\nabla \times B) \times n|_{\partial \Omega} = 0$,其中 $n$ 为外法向单位向量。该系统描述导电流体(如等离子体、液态金属)在磁场作用下的运动,其耗散项 $\mu \Delta u$ 和 $\nu \Delta B$ 分别刻画动量粘性耗散与磁场电阻扩散。当 $\mu, \nu \to 0$ 时,形式极限为理想 MHD 系统:

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u + \nabla P = (B \cdot \nabla) B, \\ \partial_t B + (u \cdot \nabla) B = (B \cdot \nabla) u, \\ \nabla \cdot u = \nabla \cdot B = 0, \\ (u, B)|_{t=0} = (u_0, B_0), \end{cases} \tag{1.2} $$

配以无穿透边界条件:$u \cdot n|_{\partial \Omega} = 0$, $B \cdot n|_{\partial \Omega} = 0$。该极限问题本质是验证 Leray–Hopf 弱解序列 $\{(u^{\varepsilon}, B^{\varepsilon})\}$(对应 $\mu = \varepsilon^p$, $\nu = \varepsilon^q$)是否强收敛于 (1.2) 的弱解。其困难远超 Navier–Stokes 的零粘性极限:首先,速度场满足无滑移边界条件 $u = 0$,而磁场满足理想导体边界条件 $B \cdot n = 0$, $(\nabla \times B) \times n = 0$,二者诱导出尺度不同、强耦合的边界层;其次,当 $p \ne q$ 时,$\mu$ 与 $\nu$ 衰减速率不同,导致边界层厚度分别为 $O(\varepsilon^{p/2})$ 与 $O(\varepsilon^{q/2})$,需统一处理;再者,非线性项 $(B \cdot \nabla) B$ 与 $(u \cdot \nabla) B$ 通过 Lorentz 力与磁感应机制深度耦合,破坏能量结构的对称性。物理上,此极限直接关联 MHD 湍流中能量级联与耗散机制——Kolmogorov 理论预言高雷诺数下能量耗散率与粘性无关,Onsager 猜想则指出弱解能量守恒依赖于 $C^{\alpha}$ 正则性阈值。但对 MHD,Faraco–Lindberg [12] 已证在 $L^2_t L^2_x$ 框架下,即使 $\mu, \nu \to 0$,弱极限 $(u,B)$ 未必满足 (1.2),因边界层能量未被有效控制。因此,本文试图填补的关键空白在于:建立适用于 MHD 的 Kato 型判据,即给出边界层内总耗散 $\int_0^T \int_{\Gamma_{4\delta}} (\varepsilon^p |\nabla u^{\varepsilon}|^2 + \varepsilon^q |\nabla B^{\varepsilon}|^2) \, dx dt \to 0$ 足以保证全局耗散消失及解收敛;并首次严格处理磁普朗特数 $Pm = \mu / \nu = \varepsilon^{p-q}$ 在 $\varepsilon \to 0$ 下趋于 $0$ 或 $\infty$ 的双重极限,此前数学文献中对此无任何严格结果。

2. 主要定理或主要结果

在 $\Omega \subset \mathbb{R}^3$ 为 $C^2$ 有界区域、$T \in (0,\infty)$ 的设定下,设 $\mu = \varepsilon^p$, $\nu = \varepsilon^q$ ($p,q > 0$),记 $\delta(\varepsilon) = \min\{\varepsilon^p, \varepsilon^q\}$,$\Gamma_h = \Omega \setminus \Omega_h$ 为边界层区域。若初始数据 $(u^{\varepsilon}_0, B^{\varepsilon}_0) \to (u_0, B_0)$ 于 $L^2(\Omega)$,且 Leray–Hopf 解 $(u^{\varepsilon}, B^{\varepsilon})$ 满足:(i) 内部正则性:$(u^{\varepsilon}, B^{\varepsilon})$ 在 $L^3(0,T; B^{\alpha_1,\infty}_3(\Omega_\delta) \times B^{\alpha_2,\infty}_3(\Omega_\delta))$ 中一致有界,其中 $\alpha_1,\alpha_2$ 满足:若 $p \ge q$,则 $\alpha_i \in \left(\frac{5}{6} - \frac{q}{2p}, \frac{5}{6}\right)$;若 $p < q$,则 $\alpha_i \in \left(\frac{5}{6} - \frac{p}{2q}, \frac{5}{6}\right)$;(ii) 边界层控制:$u^{\varepsilon}, B^{\varepsilon}$ 在 $L^4(0,T; L^\infty(\Gamma_{4\delta}))$ 中一致有界,压力 $P^{\varepsilon}$ 在 $L^2(0,T; L^\infty(\Gamma_{4\delta}))$ 中一致有界;(iii) 边界层耗散消失:

$$ \lim_{\varepsilon \to 0^+} \int_0^T \int_{\Gamma_{4\delta}} \left( \varepsilon^p |\nabla u^{\varepsilon}|^2 + \varepsilon^q |\nabla B^{\varepsilon}|^2 \right) dx dt = 0, \tag{1.8} $$

则必有全局耗散消失:

$$ \lim_{\varepsilon \to 0^+} \int_0^T \int_{\Omega} \left( \varepsilon^p |\nabla u^{\varepsilon}|^2 + \varepsilon^q |\nabla B^{\varepsilon}|^2 \right) dx dt = 0, \tag{1.9} $$

且存在子列使得 $(u^{\varepsilon}, B^{\varepsilon}) \to (u,B)$ 强收敛于 $L^3(0,T; L^3(\Omega))$,极限 $(u,B)$ 是理想 MHD 系统 (1.2) 的弱解。特别地,当 $p > q$ 时 $Pm \to \infty$,当 $p < q$ 时 $Pm \to 0$,这是首个从纯数学角度确立的磁普朗特数极限结果。

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

证明的核心障碍在于:(1) 速度与磁场边界层厚度不匹配($\varepsilon^{p/2}$ vs $\varepsilon^{q/2}$)导致标准 Kato 边界层分析失效;(2) 非线性耦合项 $(u \cdot \nabla) B$ 与 $(B \cdot \nabla) u$ 在边界附近无法分离,破坏经典能量估计的对称性;(3) 压力梯度项 $\nabla P^{\varepsilon}$ 在边界层内缺乏先验控制。作者的真正创新点有三:

第一,构造“统一尺度边界层叶状分解”:定义 $\delta(\varepsilon) = \min\{\varepsilon^p, \varepsilon^q\}$,并依递推公式 (2.10) 构造严格递增序列 $\{\beta_n\}_{n=1}^N$ 满足 $0 = \beta_0 < \beta_1 < \cdots < \beta_N = 1$,进而将边界层 $\Gamma_{\delta}$ 分解为重叠带 $V_n = \Omega_{2\delta^{\beta_n}} \setminus \Omega_{2\delta^{\beta_{n-1}} + 2\delta^{\beta_n}}$。此分解确保在每一 $V_n$ 上可对 $u^{\varepsilon}, B^{\varepsilon}$ 同时进行尺度为 $\delta^{\beta_n}$ 的局部磨光(见 (2.24)),且关键性质 $\nabla(\xi_n + \xi_{n+1})^2 = 0$ 在交集上成立,避免磨光后函数越界。

第二,引入 Elsässer 变量并重构非线性通量:定义 $z^+ = u + B$, $z^- = u - B$,将原系统重写为 $\partial_t z^\pm + (z^\mp \cdot \nabla) z^\pm = \pm \varepsilon^p \Delta u \mp \varepsilon^q \Delta B + \nabla P$。此变换解耦了对流与磁伸展项,使非线性通量 $z^\mp \otimes z^\pm$ 成为分析焦点。作者据此构造对称非线性函数 (3.58)–(3.64),并在每一叶 $V_n$ 上实施空间局部化与平均,从而将总误差分解为 $\sum_{i=1}^5 R_i$(见 (3.3))。

第三,发展分层磨光与梯度补偿估计:对 $R_2 = \varepsilon^q \int_0^T \int_\Omega \theta (\sum_n \xi_n B^{\varepsilon}_n) \cdot (\sum_n \xi_n \Delta B^{\varepsilon}_n) dx dt$,利用局部正交性 $\xi_k \xi_m = 0$(当 $|k-m| \ge 2$)及分部积分,将其拆为 $R_{21} + R_{22} + R_{23}$。其中 $R_{21}$ 主导项被控为 $\lesssim \varepsilon^q \sum_{|k-m|\le 1} \int_0^T \int_{V_k \cap V_m} |\nabla B^{\varepsilon}_k| |\nabla B^{\varepsilon}_m| dx dt$,再结合 Besov 空间磨光估计 $\|\nabla f_\tau\|_{L^r} \le C \tau^{\alpha-1} \|f\|_{B^{\alpha,\infty}_r}$(见 (2.4))与 Hölder 指数阈值 (1.6),最终由假设 (1.8) 控制其趋于零。类似技术亦用于处理 $R_1$(速度耗散)与 $R_3$–$R_5$(压力与非线性误差)。

4. 与之前工作的比较

与已有工作相比,本文在三个维度实现实质性推进:

(1) 相较 Kato [16] 及 Bardos–Titi–Wiedemann [1] 的 Navier–Stokes 零粘性极限:Kato 假设边界层内 $\int_0^T \int_{\Gamma_\delta} \mu |\nabla u^{\mu}|^2 dx dt \to 0$ 即可推出 Euler 解收敛,但其框架仅适用于单一耗散系数与单一边界层;[1] 要求 $u^{\mu} \in L^3_t B^{\alpha,\infty}_3(\Omega_\delta)$ 且 $\alpha > 5/6$,而本文将阈值降至 $\alpha_i > 5/6 - \min\{q/(2p), p/(2q)\}$,且必须同时控制 $u$ 与 $B$ 的正则性。代价是引入更复杂的叶状分解与双尺度磨光,而非单一尺度 $\delta = \mu^{1/2}$。

(2) 相较 Faraco–Lindberg [12] 的 MHD 同时极限:[12] 仅设 $\mu = \nu \to 0$,证明磁螺旋度守恒但弱极限不满足理想 MHD;本文通过设定 $\mu = \varepsilon^p$, $\nu = \varepsilon^q$,首次严格获得 $Pm \to 0$ 或 $\infty$ 的极限,揭示了不同耗散速率下湍流尺度分离的新机制。代价是需更强的边界层耗散消失条件 (1.8),而 [12] 仅依赖 $L^2_t L^2_x$ 积分有界性。

(3) 相较 Chen–Liang–Wang [5] 的 Navier–Stokes 边界层叶状分解:[5] 为单一速度场设计 $\beta_n$ 序列以处理 $\mu^{1/2}$ 尺度,本文将其推广至双参数情形,递推式 (2.10) 显含 $p/q$ 或 $q/p$ 比值,且叶数 $N$ 依赖于 $p,q$ 关系(见 (2.12)–(2.15))。此外,[5] 的非线性分析基于 $u \otimes u$,而本文必须处理耦合张量 $u \otimes B$ 与 $B \otimes u$,导致误差项 $R_4,R_5$ 结构更复杂,需借助 Elsässer 变量重构才能分离。从提供内容中无法确认本文结果是否适用于无界域或周期边界,亦未显式展开与 Constantin–E–Titi [8] 关于 Onsager 正则性阈值的比较。


2. 通过热流与定量微分获得输运方程的正则性估计

Regularity estimates in transport equations via heat flow and quantitative differentiation

作者: Elia Brué, Maria Colombo, Guido De Philippis, Carl Johan Peter Johansson

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13751v1

摘要翻译: 本文的目的有二:第一,我们证明了 Ambrosio-Trevisan 交换子的定量估计,该估计蕴含对数 Sobolev 正则性的传播;第二,我们借助定量微分方法,证明了 DiPerna-Lions 交换子的类似估计。

Comment: 17 pages

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

1. 研究对象和背景

本文研究欧几里得空间 $\mathbb{R}^d$ 上的线性输运方程

$$ \partial_t u_t + b_t \cdot \nabla u_t = 0 \quad \text{in } [0,T] \times \mathbb{R}^d \tag{Tr} $$

的解 $u_t$ 的正则性传播问题,其中速度场 $b_t: \mathbb{R}^d \to \mathbb{R}^d$ 是给定的向量场。该方程在流体力学、动力系统、偏微分方程弱解理论及最优传输中具有根本地位:其解描述了由 $b_t$ 生成的流形上标量量 $u$ 的守恒演化;当 $\operatorname{div} b_t = 0$(不可压)时,经典光滑解满足 $L^q$ 范数守恒 $\|u_t\|_{L^q} = \|u_0\|_{L^q}$,且对任意 Lipschitz 函数 $\beta$ 满足 $\int_{\mathbb{R}^d} \beta(u_t)\,dx = \int_{\mathbb{R}^d} \beta(u_0)\,dx$,此恒等式蕴含解的唯一性与稳定性。然而,在 DiPerna–Lions 理论框架下,$b_t$ 仅被假设为 $b \in L^1([0,T]; L^\infty \cap W^{1,p}(\mathbb{R}^d; \mathbb{R}^d))$,$u \in L^\infty([0,T]; L^q(\mathbb{R}^d))$,且 $1/p + 1/q = 1$,此时 $u_t$ 仅为分布解,上述恒等式不再自动成立,必须通过正则化与误差控制来恢复。核心障碍在于:对标准磨光子 $\rho_\varepsilon$ 或热半群 $P_\varepsilon$,定义的正则化解 $u^\varepsilon_t$ 不再精确满足 (Tr),而引入余项——即 DiPerna–Lions 与 Ambrosio–Trevisan 两类交换子:

$$ R_\varepsilon(b,u)(x) := \int_{\mathbb{R}^d} u(y)\big(b(x)-b(y)\big)\cdot \nabla \varphi_\varepsilon(x-y)\,dy, $$

$$ C_\varepsilon(b,u)(x) := b_t \cdot \nabla(P_\varepsilon u_t) - P_\varepsilon(b_t \cdot \nabla u_t), $$

其中 $\varphi_\varepsilon(z) = \varepsilon^{-d}\varphi(z/\varepsilon)$,$P_s u = u * G_s$ 为热半群,$G_s(z) = (4\pi s)^{-d/2} e^{-|z|^2/(4s)}$。已知若 $\|R_\varepsilon(b,u)\|_{L^1} \to 0$ 或 $\|C_\varepsilon(b,u)\|_{L^1} \to 0$,则可证明 (Tr) 解的唯一性与 $\beta$-守恒律成立。但传统估计仅给出 $\|R_\varepsilon\|_{L^1} \le C \varepsilon^{1-1/p} \|\nabla b\|_{L^p} \|u\|_{L^q}$,无法刻画多尺度累积误差;更关键的是,现有工作(如 [DL89, Amb04, AT14])未提供对 $\int_0^1 \|R_\varepsilon\|_{L^1} d\varepsilon/\varepsilon$ 或 $\int_0^1 \|C_\varepsilon\|_{L^1} d\varepsilon/\varepsilon$ 的定量控制,而这正是刻画“对数型正则性”(log-Sobolev regularity)传播所必需的尺度不变积分。本文旨在填补这一空白:一方面,利用热流 $P_s$ 的精细结构,建立 Ambrosio–Trevisan 交换子 $C_\varepsilon$ 的 $L^p$-$L^q$-$L^r$ 多重范数积分估计;另一方面,首次将 Dorronsoro 定量微分理论(quantitative differentiation)系统应用于 DiPerna–Lions 交换子 $R_\varepsilon$,以克服其非卷积结构带来的技术瓶颈,从而获得 $L^2(d\varepsilon/\varepsilon)$ 型平方函数控制。

2. 主要定理或主要结果

本文建立了两类交换子的定量估计,并由此导出输运方程解的对数 Sobolev 正则性传播定理。首先,对 Ambrosio–Trevisan 交换子 $C_\varepsilon(b,u)$,在 $\operatorname{div} b = 0$ 下证明:对任意 $1 < p \le \infty$, $1 < q < \infty$, $1 < r < \infty$ 满足 $1/p + 1/q + 1/r = 1$,有

$$ \int_0^1 \left| \int_{\mathbb{R}^d} C_\varepsilon(b,u)\, \phi \,dx \right| \frac{d\varepsilon}{\varepsilon} \le C(d,p,q,r) \|\nabla b\|_{L^p} \|u\|_{L^q} \|\phi\|_{L^r} \tag{1.3} $$

;进而推出对解 $(u_t,b_t)$ 满足 $b \in L^1_t(L^\infty \cap W^{1,p}_x)$, $u \in L^\infty_t(L^q \cap L^r_x)$ 的对数正则性传播:

$$ [u_t]^2_{\log} - [u_0]^2_{\log} \le C(d,p,q,r) \|\nabla b\|_{L^1_t L^p_x} \|u\|_{L^\infty_t L^q_x} \|u\|_{L^\infty_t L^r_x}, $$

其中对数半范数定义为

$$ [u]^2_{\log} := \int_0^1 \|P_\varepsilon u - u\|^2_{L^2} \frac{d\varepsilon}{\varepsilon}. $$

该估计亦推广至 $\nabla b \in L^1_t H^1_x$(Hardy 空间)情形(见附录定理 A.1)。其次,对 DiPerna–Lions 交换子 $R_\varepsilon(b,u)$,在相同散度为零假设下证明:对 $1<p,q<\infty$ 满足 $1/p+1/q=1$,有

$$ \left( \int_0^\infty \left| \int_{\mathbb{R}^d} R_\varepsilon(b,u)\, \phi \,dx \right|^2 \frac{d\varepsilon}{\varepsilon} \right)^{1/2} \le C(d,p,q) \|b\|_{W^{1,p}} \big( \|u\|_{L^q} \|\phi\|_{L^\infty} + \|u\|_{L^\infty} \|\phi\|_{L^q} \big) \tag{1.4} $$

;此平方函数估计蕴含正则性传播(命题 4.10),其本质是将 $R_\varepsilon$ 分解为三项并分别用定量微分控制:第一、三项由 $\Omega^1_b(x,t)$(Dorronsoro $\beta$-数)与 Hardy–Littlewood 极大算子 $M$ 控制,第二项则利用 $\nabla^{\mathrm{app}}_{B_{2\varepsilon}(\bar x)} b$ 的迹为零性质及引理 4.7(一类广义 Littlewood–Paley 平方函数有界性)完成估计。

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

本文的核心创新在于将两种深刻但此前未用于输运方程交换子估计的分析工具——热流结构与定量微分——创造性地结合。首先,针对 Ambrosio–Trevisan 交换子 $C_\varepsilon$,作者的关键突破是挖掘其内在的双参数热流表示:由命题 2.2,$\int C_\varepsilon(b,u)\phi = 2 \int_0^\varepsilon \int_{\mathbb{R}^d} \nabla b \cdot \nabla P_s u \cdot \nabla P_{\varepsilon-s} \phi \,dx ds$。此恒等式非平凡,它源于热半群 $P_s$ 与拉普拉斯算子 $\Delta$ 的交换性及 $\operatorname{div} b = 0$ 下的分部积分技巧,将 $C_\varepsilon$ 转化为一个关于 $\nabla b$、$\nabla P_s u$ 和 $\nabla P_{\varepsilon-s} \phi$ 的双线性形式。随后,作者构建了一个精巧的三重积分估计框架:将 $\int_0^1 \cdots d\varepsilon/\varepsilon$ 与内层 $\int_0^\varepsilon \cdots ds$ 交换次序,再通过变量替换与 Marcinkiewicz 插值,将问题归结为对平方函数 $Sf = (\int_0^\infty |\nabla P_s f|^2 ds)^{1/2}$ 的 $L^p$ 有界性(引理 3.1)的应用。这避免了传统方法中对 $\varepsilon$ 的粗略幂次估计,实现了对 $d\varepsilon/\varepsilon$ 尺度不变测度的精确积分控制。其次,针对 DiPerna–Lions 交换子 $R_\varepsilon$,其非卷积结构使热流方法失效,作者转而采用定量微分。其真正创新点在于提出一种“中心漂移分解”:在表达式 $\int R_\varepsilon(b,u)\phi$ 中引入辅助变量 $\bar x$,并将 $b(x)-b(y)$ 分解为 $b(x)-P^{1}_{B_{2\varepsilon}(\bar x)}b(x) + \nabla^{\mathrm{app}}_{B_{2\varepsilon}(\bar x)}b \cdot (x-y) + P^{1}_{B_{2\varepsilon}(\bar x)}b(y)-b(y)$。此分解巧妙规避了直接估计 $P^{1}_{B_{\varepsilon}(y)}b(x)-b(y)$ 的困难;其中中间项 $\nabla^{\mathrm{app}}_{B_{2\varepsilon}(\bar x)}b$ 的迹为零性(引理 4.6)至关重要——它允许将 $u(y)$ 替换为 $u(y)-f(x,\varepsilon)$ 并选择 $f$ 使该差分满足引理 4.7 的条件,从而启用广义平方函数估计;首尾两项则直接由 Dorronsoro $\beta$-数 $\Omega^1_b(\bar x, \varepsilon)$ 与极大函数控制。整个证明依赖于对 $\nabla^{\mathrm{app}}_B b$ 的精细估计(引理 4.4)、其迹零性(引理 4.6)以及一类具有零均值核的 Calderón–Zygmund 平方函数算子的 $L^p$ 有界性(引理 4.7),构成了一套全新的、面向交换子结构的定量微分分析框架。

4. 与之前工作的比较

本文与既有工作存在三重实质性推进。第一,相较于 DiPerna–Lions [DL89] 和 Ambrosio [Amb04] 的定性交换子收敛结果 $\|R_\varepsilon\|_{L^1} \to 0$、$\|C_\varepsilon\|_{L^1} \to 0$,本文提供了首个 $\int_0^1 \|C_\varepsilon\|_{L^1} d\varepsilon/\varepsilon$ 的显式上界(定理 1.1),该尺度不变积分是刻画 log-Sobolev 半范数 $[\cdot]_{\log}$ 变化的必要且充分工具,而此前 [Lég18, BN21, MS24] 的相关工作仅依赖于 Littlewood–Paley 或加权范数,未触及交换子本身的多尺度积分结构。第二,与 Huysmans–Said [HS24] 利用调和分析(Calderón–Zygmund 理论)对 $R_\varepsilon$ 所做的定量估计相比,本文的定理 1.3 在技术路径与结论强度上均有差异:[HS24] 给出的是 $\|R_\varepsilon\|_{L^1} \le C \varepsilon^{\alpha} \|\nabla b\|_{L^p} \|u\|_{L^q}$ 形式的单尺度估计,而本文的 $L^2(d\varepsilon/\varepsilon)$ 型平方函数估计(1.4)天然兼容多尺度叠加,且不依赖于 $b$ 的更高阶正则性(如 $W^{2,p}$),仅需 $W^{1,p}$,适用范围更广;代价是估计常数依赖于 $\|b\|_{W^{1,p}}$ 而非 $\|\nabla b\|_{L^p}$,且对测试函数 $\phi$ 需 $L^q \cap L^\infty$ 条件。第三,本文首次将定量微分(Dorronsoro [Dor85])系统引入输运方程领域,此前该工具主要用于几何测度论(如 Reifenberg 型定理、奇异集分析)与椭圆方程正则性(如 [AT15]),其在双线性、非卷积型算子(如 $R_\varepsilon$)上的应用属首创;相较 Meyer–Seis [MS24] 的 Littlewood–Paley 方法,定量微分更直接刻画函数在各尺度下的仿射近似质量,无需傅里叶变换,因而天然适用于非齐次或非平移不变设置。从假设看,本文严格继承 DiPerna–Lions 框架($b \in L^1_t W^{1,p}_x$, $\operatorname{div} b = 0$),未放宽可积性指数,但显著加强了结论的定量性;局限在于目前结果均要求散度为零,而 Remark 2.3 指出 $\operatorname{div} b \in L^p$ 的推广虽可行但未详述,且对 $p=1$ 或临界情形(如 $b \in \operatorname{BMO}$)尚未覆盖。


3. 大普朗特数下三维无扩散Boussinesq系统中温度块的全局正则性

Global regularity of temperature patches for the 3D non-diffusive Boussinesq system with large Prandtl number

作者: Qianyun Miao, Jiakun Yang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13803v1

摘要翻译: 迄今为止,三维无扩散Boussinesq系统在大初值情形下的强解全局适定性仍是一个突出的公开问题。本文在大普朗特数(Prandtl number)参数区域中解决了该问题。更准确地说,我们证明:当普朗特数足够大(其阈值仅依赖于初值 $(u_0,θ_0)$ 的某个尺度不变范数)时,该三维Boussinesq系统对初值 $(u_0,θ_0)\in H^{\frac{1}{2}}(\mathbb{R}^3) \times (L^1\cap L^s(\mathbb{R}^3))$(其中 $s>3$)存在唯一全局强解;此外,对非常数温度块初值,我们建立了演化温度块边界在 $C^{1,γ}$、$W^{2,\infty}$ 及 $C^{2,γ}$($0<γ<1$)空间中的全局正则性保持性,并给出相应估计在大普朗特数区域中的一致性。进一步地,我们严格论证了当普朗特数趋于无穷时的极限过程:三维Boussinesq系统的块解收敛于三维Stokes-transport系统的唯一块解,且块边界在 $C^{1,γ}$、$W^{2,\infty}$ 和 $C^{2,γ}$ 中的正则性在全局时间内得以保持。特别地,我们关于三维Stokes-transport系统的结果可视为Grayer II [ARMA 2023] 关于二维Stokes-transport系统主结果的三维类比。

Comment: 31 pages

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

1. 研究对象和背景

本文研究三维无扩散Boussinesq系统(B)的Cauchy问题,其方程组为:

$$ \begin{cases} \partial_t \theta + u \cdot \nabla \theta = 0, \\ \frac{1}{\mathrm{Pr}}\left(\partial_t u + u \cdot \nabla u\right) - \nu \Delta u + \nabla p = \theta e_3, \\ \nabla \cdot u = 0, \\ (\theta, u)|_{t=0}(x) = (\theta_0, u_0)(x), \end{cases} \tag{B} $$

其中 $(t,x) \in \mathbb{R}_+ \times \mathbb{R}^3$,$e_3 = (0,0,1)^\top$,$\nu > 0$ 为运动粘性系数,$\mathrm{Pr} > 0$ 为无量纲普朗特数(定义为 $\nu / \kappa$,$\kappa$ 为热扩散率),未知量为标量温度 $\theta$、速度场 $u = (u_1,u_2,u_3)^\top$ 和压力 $p$。该系统在海洋与大气对流建模、Rayleigh-Bénard对流理论中具有核心地位。数学上,它介于Navier-Stokes方程与Stokes-transport系统之间:当 $\theta \equiv 0$ 时退化为不可压Navier-Stokes方程;当 $\mathrm{Pr} \to \infty$ 时形式上趋于Stokes-transport系统

$$ \begin{cases} \partial_t \theta + u \cdot \nabla \theta = 0, \\ -\nu \Delta u + \nabla p = \theta e_3, \\ \nabla \cdot u = 0, \\ \theta|_{t=0} = \theta_0. \end{cases} \tag{ST} $$

关键难点在于:对 $\mathrm{Pr}=1$ 的三维粘性Boussinesq系统,全局强解的存在唯一性仍是未解难题——其困难程度堪比千禧年难题中的三维Navier-Stokes方程。Danchin–Paicu [21] 在初始数据 $(u_0,\theta_0) \in L^2 \times L^2$ 下构造了有限能量Leray-Hopf型弱解,并在尺度不变小初值假设 $\|u_0\|_{L^{3,\infty}} + \|\theta_0\|_{L^1} \le c\nu$ 下得到全局强解;但该小性条件无法移除。而其极限系统(ST)却在无小性假设下已知全局适定:Mecherbet [50] 对 $\theta_0 \in L^1 \cap L^\infty$(具有限一阶矩)证明了全局存在唯一性;Mechebet–Sueur [51] 将初值空间延拓至 $L^1 \cap L^p$($p \ge 3$);Grayer II [30] 在二维情形下建立了Yudovich型初值 $L^1 \cap L^\infty$ 下的全局适定性及 $C^{1,\gamma}$、$W^{2,\infty}$、$C^{2,\gamma}$ 边界正则性的全局保持。然而,三维情形下(ST)的边界正则性传播此前仅知球形斑块的行波解 [50],一般非恒定温度斑块(即 $\theta_0 = \theta_0(x) \mathbf{1}_{D_0}(x)$,其中 $D_0 \subset \mathbb{R}^3$ 为有界单连通区域)的 $C^{1,\gamma}$ 等高阶边界正则性是否全局保持,仍属空白。本文旨在填补这一根本性空白:一方面,利用大普朗特数 $\mathrm{Pr} \gg 1$ 提供的额外耗散抑制机制,突破传统小初值限制,建立三维Boussinesq系统在 $H^{1/2} \times (L^1 \cap L^s)$($s>3$)初值下的全局强解存在唯一性;另一方面,针对温度斑块解,严格证明其边界 $C^{1,\gamma}$、$W^{2,\infty}$、$C^{2,\gamma}$ 正则性在时间演化中全局保持,且所有估计关于 $\mathrm{Pr} \in [\mathrm{Pr}_*, \infty)$ 一致成立;最终,通过紧性论证严格验证 $\mathrm{Pr} \to \infty$ 极限,将上述斑块正则性结果传递至三维Stokes-transport系统(ST),从而确立其作为二维Grayer II [30] 结果之真正三维类比的地位。

2. 主要定理或主要结果

定理1.1(大Prandtl数下Boussinesq系统的全局适定性与斑块正则性保持):设 $u_0 \in H^{1/2}(\mathbb{R}^3)$ 为无散度向量场,$\theta_0 \in L^1 \cap L^s(\mathbb{R}^3)$,$s > 3$。存在普适小常数 $c_* > 0$(由式(3.4)与(3.6)确定),使得若

$$ \mathrm{Pr} \ge \mathrm{Pr}_* \triangleq \frac{\|u_0\|_{L^{3,\infty}(\mathbb{R}^3)} + \|\theta_0\|_{L^1(\mathbb{R}^3)}}{c_*}, $$

则系统(B)存在唯一全局解 $(u, \theta)$,满足对任意 $T > 0$ 及 $q \in (3,6)$,

$$ u \in C([0,T]; H^{1/2}) \cap L^2([0,T]; H^{3/2} \cap L^\infty) \cap \widetilde{L}^1([0,T]; \dot{B}^{1+3/q}_{q,\infty}), \quad \theta \in L^\infty([0,T]; L^1 \cap L^s). $$

进一步,对非恒定温度斑块初值 $\theta_0(x) = \theta_0(x) \mathbf{1}_{D_0}(x)$($D_0 \subset \mathbb{R}^3$ 为有界单连通区域),在条件(1.3)下: (1) 若 $\partial D_0 \in C^{1,\gamma}$($0<\gamma<1$),$\theta_0 \in L^\infty(D_0)$,且

$$ u_0 \in \begin{cases} H^1(\mathbb{R}^3), & \gamma \in (0,1/2), \\ H^1 \cap W^{1,p}(\mathbb{R}^3), & \text{对某 } p \in \left(\frac{3}{2-\gamma}, \infty\right), \; \gamma \in [1/2,1), \end{cases} $$

则 $\theta(x,t) = \theta_0(X_t^{-1}(x)) \mathbf{1}_{D(t)}(x)$,其中 $D(t) = X_t(D_0)$,且 $\partial D(t) \in L^\infty([0,T]; C^{1,\gamma})$ 关于 $\mathrm{Pr} \in [\mathrm{Pr}_*, \infty)$ 一致成立。 (2) 若额外 $\partial D_0 \in W^{2,\infty}$,$\theta_0 \in C^\mu(D_0)$($\mu \in (0,1)$),$u_0 \in H^1 \cap W^{1,p}(\mathbb{R}^3)$(对某 $p \in (3,\infty)$),则 $\partial D(t) \in L^\infty([0,T]; W^{2,\infty})$ 一致成立。 (3) 若额外 $\partial D_0 \in C^{2,\gamma}$($0<\gamma<1$),$\theta_0 \in C^\gamma(D_0)$,$u_0 \in H^1 \cap W^{2,p}(\mathbb{R}^3)$(对某 $p \in \left(\frac{3}{2-\gamma}, \infty\right)$),则 $\partial D(t) \in L^\infty([0,T]; C^{2,\gamma})$ 一致成立。 定理1.4(无穷Prandtl数极限与Stokes-transport系统斑块正则性):对任意 $T>0$ 及 $\mathrm{Pr} \in [\mathrm{Pr}_*, \infty)$,设 $(u^{\mathrm{Pr}}, \theta^{\mathrm{Pr}})$ 为定理1.1所构全局解,对应初值 $\theta_0 = \theta_0(x) \mathbf{1}_{D_0}(x)$ 满足 $\partial D_0 \in C^{1,\gamma}$($0<\gamma<1$)等条件,且 $u_0^{\mathrm{Pr}} \to u_0 \in H^1(\mathbb{R}^3)$ 强收敛并满足相容性条件 $-\Delta u_0 + \nabla p_0 = \theta_0 e_3$。则当 $\mathrm{Pr} \to \infty$ 时,$(u^{\mathrm{Pr}}, \theta^{\mathrm{Pr}})$ 弱收敛于 $(u, \theta)$,后者为系统(ST)的唯一全局解,且 $\theta(x,t) = \theta_0(X_t^{-1}(x)) \mathbf{1}_{D(t)}(x)$,其中 $D(t) = X_t(D_0)$,并满足: (1) 若 $\partial D_0 \in W^{2,\infty}$ 等附加条件成立,则 $\partial D(t) \in L^\infty([0,T]; W^{2,\infty})$; (2) 若 $\partial D_0 \in C^{2,\gamma}$ 等附加条件成立,则 $\partial D(t) \in L^\infty([0,T]; C^{2,\gamma})$。

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

本文的核心技术障碍在于:三维Boussinesq系统在 $\mathrm{Pr}=1$ 时缺乏控制 $\nabla u$ 的全局先验估计,导致无法处理大初值;而在大 $\mathrm{Pr}$ 下,动量方程中 $\frac{1}{\mathrm{Pr}}\partial_t u$ 项虽提供额外耗散,但其系数衰减亦削弱了对 $u$ 时间正则性的直接控制。作者的真正创新点在于构建了一个普适的、关于 $\mathrm{Pr}$ 一致的全局能量框架,并引入一个精巧的“好未知量”(good unknown)$\Gamma$ 来解耦速度与温度的耦合效应。具体而言: 第一,全局适定性证明的关键是建立 $u \in L^\infty_T(H^{1/2}) \cap L^2_T(H^{3/2})$ 的一致估计。作者通过经典能量法结合 $L^\infty_T(L^{3,\infty})$ 型全局估计(类比Danchin–Paicu [21] 的小初值情形)实现:将动量方程乘以 $(-\Delta)^{1/2} u$ 并积分,利用 $\frac{1}{\mathrm{Pr}}\|\partial_t u\|_{H^{1/2}}^2$ 项提供主导耗散,同时通过插值与Gagliardo-Nirenberg不等式控制非线性项 $\int (u \cdot \nabla u) \cdot (-\Delta)^{1/2} u$,最终导出关于 $\|u\|_{H^{1/2}}^2 + \int_0^t \|u\|_{H^{3/2}}^2 ds$ 的闭合微分不等式,其右端仅含 $\|\theta_0\|_{L^1}$ 和 $\|u_0\|_{L^{3,\infty}}$,从而确定 $\mathrm{Pr}_*$ 的阈值。 第二,斑块边界正则性保持的核心工具是“好未知量” $\Gamma := \Omega - (R_{-1,2}\theta, -R_{-1,1}\theta, 0)^\top$,其中 $\Omega = \nabla \wedge u$ 为涡量,$R_{-1,j} = \partial_j (-\Delta)^{-1}$。该量满足一个关键性质:在极限系统(ST)中 $\Gamma \equiv 0$,且其演化方程(4.5)不含 $\theta$ 的一阶导数项,从而避免了 $\theta$ 的低正则性对 $\nabla u$ 估计的污染。作者据此将 $\nabla u$ 分解为

$$ \nabla u = (-\Delta)^{-1} \nabla \nabla \wedge \Omega = \Lambda^{-2} \nabla \nabla \wedge \Gamma + \nabla^2 \partial_3 \Lambda^{-4} \theta + \Lambda^{-2} \nabla \theta \otimes e_3, $$

其中 $\Lambda = (-\Delta)^{1/2}$。此分解将对 $\nabla u$ 的估计转化为对 $\Gamma$、$\theta$ 及其组合的估计。对 $C^{1,\gamma}$ 保持,作者证明 $\nabla u \in L^1_T(C^\gamma)$,这依赖于:(i) 利用 $\Gamma$ 方程与交换子估计获得 $\Gamma \in L^\infty_T(L^2) \cap L^2_T(H^1)$ 的一致界;(ii) 利用 $\theta \in L^\infty_T(L^1 \cap L^\infty)$ 及分数阶积分算子的Sobolev嵌入,得到 $\nabla^2 \partial_3 \Lambda^{-4} \theta$ 与 $\Lambda^{-2} \nabla \theta$ 均属于 $L^1_T(C^\gamma)$。 第三,对 $W^{2,\infty}$ 与 $C^{2,\gamma}$ 保持,作者采用切向向量场系统 $W = \{W^i\}_{i=1}^5$ 的方法(引理2.3)。关键创新在于对 $\nabla W$ 的估计:对 $W^{2,\infty}$,作者对方程 $\partial_t W^i + u \cdot \nabla W^i = \partial_{W^i} u$ 应用 $\nabla$ 后使用最大值原理,其难点在于控制 $\|\partial_W \nabla u\|_{L^1_T(L^\infty)}$。作者将此项按前述分解展开,并重点处理交换子 $[\nabla^3 \Lambda^{-4}, W \cdot \nabla]\theta$,为此建立关键交换子估计(2.11);对 $C^{2,\gamma}$,作者需控制 $\|\partial_W \nabla u\|_{L^1_T(C^\gamma)}$,其创新在于将 $\partial_W \Gamma$ 分解为零初值与零强迫两部分,分别应用平滑估计(2.18)-(2.19),并结合最优的乘积估计(2.9)与交换子估计(2.15),最终在更弱的初值假设 $p > \frac{3}{2-\gamma}$ 下完成估计,较前人工作[42]显著放宽了对 $u_0$ 的正则性要求。

4. 与之前工作的比较

本文与既有工作的核心推进体现在三个维度: (1)全局适定性假设的突破:此前所有三维Boussinesq系统的全局强解结果均依赖初值小性。Danchin–Paicu [21] 要求 $\|u_0\|_{L^{3,\infty}} + \|\theta_0\|_{L^1} \le c\nu$;Dong–Sun [23] 仅在带边界的薄层 $\mathbb{R}^2 \times (0,1)$ 中证明渐近稳定性;Hmidi et al. [2,34] 仅对无旋(no swirl)轴对称初值成立。本文首次在全空间 $\mathbb{R}^3$ 中,对任意大的 $H^{1/2} \times (L^1 \cap L^s)$($s>3$)初值,通过大普朗特数 $\mathrm{Pr} \ge \mathrm{Pr}_*$ 实现全局适定,将“小初值”范式替换为“大参数”范式,且 $\mathrm{Pr}_*$ 仅依赖于尺度不变范数 $\|u_0\|_{L^{3,\infty}} + \|\theta_0\|_{L^1}$,具有最优性。 (2)斑块正则性保持的维度与精度提升:在二维情形,Danchin–Zhang [22] 与Gancedo–García-Juárez [25] 证明了 $C^{1,\gamma}$ 与 $W^{2,\infty}$ 的全局保持;Chae–Miao–Xue [8] 将其推广至任意 $C^{k,\gamma}$($k \ge 1$)。但在三维,此前仅[42]在小初值下证明了 $C^{1,\gamma}$ 保持,且对 $u_0$ 要求 $W^{1,p}$($p>3$)。本文不仅将 $C^{1,\gamma}$、$W^{2,\infty}$、$C^{2,\gamma}$ 的保持推广至大初值与大 $\mathrm{Pr}$,更关键的是:对 $C^{1,\gamma}$ 保持,当 $\gamma \in [1/2,1)$ 时,将 $u_0$ 的正则性要求从 $p>3$ 放宽至 $p > \frac{3}{2-\gamma}$(例如 $\gamma = 0.6$ 时仅需 $p > 3.75$),这是本质性改进;对 $C^{2,\gamma}$ 保持,同样实现了 $p > \frac{3}{2-\gamma}$ 的最优假设。此外,本文首次在三维建立了 $W^{2,\infty}$ 的全局保持,此前文献中未见。 (3)极限过程的严格性与几何结构保持:Grayer II [30] 在二维证明了 (ST) 的 $C^{1,\gamma}/W^{2,\infty}/C^{2,\gamma}$ 保持,但其极限论证基于 $\mathbb{T}^2$ 上的正则解。本文在 $\mathbb{R}^3$ 全空间中,首次对斑块解严格证明 $\mathrm{Pr} \to \infty$ 极限,并完整保持斑块的几何结构(即 $\theta = \theta_0(X_t^{-1}(x)) \mathbf{1}_{D(t)}(x)$)与全部三类边界正则性。技术上,这依赖于 $\dot{H}^1(\mathbb{R}^3) \hookrightarrow L^6(\mathbb{R}^3)$ 的嵌入(Remark 1.5明确指出此嵌入是 $\mathbb{R}^3$ 相对于 $\mathbb{R}^2$ 的关键优势,因 $\dot{H}^1(\mathbb{R}^2)$ 不嵌入任何 $L^p$),从而保证 $u^{\mathrm{Pr}}$ 在 $L^{3/2}_T(W^{1,6})$ 中的一致有界性,为Aubin-Lions引理应用及 $\theta^{\mathrm{Pr}}$ 的强收敛提供基础。代价在于:所有结果均要求 $\mathrm{Pr} \ge \mathrm{Pr}_*$,且 $\mathrm{Pr}_*$ 依赖于初值;当 $\mathrm{Pr}$ 接近阈值时,解的某些常数可能退化。从提供内容中无法确认本文结果是否可推广至 $\mathrm{Pr} < \mathrm{Pr}_*$ 的中间参数范围,或是否适用于其他物理边界条件(如无滑移边界)。


4. 密度与速度均具非局部性的非局部守恒律的奇异极限

Singular limit for nonlocal conservation laws with non-locality in density and velocity

作者: Immanuel Ben Porat

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.14002v1

摘要翻译: 我们研究一类同时在密度和速度中引入非局部性的非局部守恒律的奇异极限。我们证明了底层 Burgers 方程 $u_t + u u_x = 0$ 的经典解可作为该奇异极限得到。该推导在激波形成前的最大时间之前均成立。主要结果回答了文献 [10] 中提出的问题。

Comment: 8 pages

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

1. 研究对象和背景

本文研究如下标量非局部守恒律的奇异极限问题:

$$ \begin{cases} \partial_t u^\varepsilon(t, x) + \partial_x\big(V_1\big(V_2(u^\varepsilon) * \eta^\varepsilon\big)\, u^\varepsilon(t, x)\big) = 0, & (t,x) \in (0,\infty) \times \mathbb{R}, \\ u^\varepsilon(0,x) = u_0(x), & x \in \mathbb{R}, \end{cases} \tag{1.1} $$

其中 $V_1,V_2 \in C^3_{\mathrm{loc}}(\mathbb{R})$ 是给定的速度函数,$\eta \in L^1(\mathbb{R}) \cap L^\infty(\mathbb{R})$ 满足 $\int_\mathbb{R} \eta(x)\,dx = 1$,且定义 $\eta^\varepsilon(x) := \frac{1}{\varepsilon}\eta\big(\frac{x}{\varepsilon}\big)$ 为标准磨光核族。该方程是交通流建模中提出的更一般模型(见文献[10]),其非局部性同时体现在密度 $u^\varepsilon$ 的复合映射 $V_2(u^\varepsilon)$ 和速度场 $V_1(V_2(u^\varepsilon)*\eta^\varepsilon)$ 中——这区别于经典非局部守恒律 $\partial_t u^\varepsilon + \partial_x(V_1(u^\varepsilon * \eta^\varepsilon) u^\varepsilon) = 0$(即 $V_2 = \mathrm{Id}$ 的退化情形,见(1.4)),也区别于仅在速度中含非局部性的模型 $\partial_t u^\varepsilon + \partial_x(V_1(u^\varepsilon) * \eta^\varepsilon) = 0$。其局部对应方程为

$$ \begin{cases} \partial_t u(t,x) + \partial_x\big(V_1(V_2(u(t,x)))\, u(t,x)\big) = 0, & (t,x) \in (0,\infty) \times \mathbb{R}, \\ u(0,x) = u_0(x), & x \in \mathbb{R}, \end{cases} \tag{1.3} $$

由于 $\eta^\varepsilon \to \delta_0$ 在分布意义下成立,(1.3) 形式上是 (1.1) 在 $\varepsilon \to 0$ 下的极限;但严格证明该奇异极限——即验证 $u^\varepsilon \to u$ 在适当范数下一致收敛——构成核心分析挑战。问题的重要性在于:一方面,(1.1) 作为交通流模型能更好刻画驾驶员对前方车流状态(密度与速度)的综合响应;另一方面,其数学结构揭示了双重非局部耦合对解的正则性传播、激波形成机制及紧性失效路径的深刻影响。已有理论中,当 $V_2 = \mathrm{Id}$ 时,文献[1,3,4,8,9,13,14] 建立了完备的适定性理论(包括多维、低正则核情形),而 Kružkov 熵解理论[15] 保证了 (1.3) 的全局弱解存在唯一性。然而,对非线性项 $V_1(u^\varepsilon * \eta^\varepsilon) u^\varepsilon$ 取极限需克服两项根本障碍:(i) $L^\infty$ 弱-* 收敛不足以控制乘积项的非线性振荡,尤其当 $u^\varepsilon$ 缺乏一致有界变差(TV)时;(ii) 即便对经典解,$\varepsilon \to 0$ 过程中 $u^\varepsilon * \eta^\varepsilon$ 与 $u^\varepsilon$ 的尺度分离会诱发误差累积,需精细估计 $\|u^\varepsilon * \eta^\varepsilon - u\|$。文献[10] 观察到 (1.1) 具单调性保持性(若 $u_0$ 单调,则 $u^\varepsilon(t,\cdot)$ 保持单调),从而导出一致 TV 界并借助 Helly 定理获得紧性;但此法严重依赖初值单调性,无法处理一般光滑初值。文献[11] 在 $V_1 = \mathrm{Id}$ 且 $\eta(x) = e^x \mathbf{1}_{(-\infty,0]}(x)$(指数核)的特殊设定下,对熵解证明了奇异极限,其关键在于该核满足 $\eta' = \eta$ 导致的微分恒等式,使能量估计可闭合;但该方法无法推广至一般各向异性核。本文旨在填补这一空白:在不假设初值单调、不限制核的具体形式(仅要求各向异性条件 (I))、且允许 $V_1,V_2$ 为一般非线性函数的前提下,对经典解证明 $\varepsilon \to 0$ 奇异极限成立,并给出最优误差阶 $O(\varepsilon)$。

2. 主要定理或主要结果

定理 1.1(主结果):设初值 $u_0 \in W^{2,\infty}(\mathbb{R})$ 且 $u_0 \geq 0$;核 $\eta$ 满足各向异性条件 (I)(即 $\eta \in L^1 \cap L^\infty$, $\operatorname{supp}(\eta) \subset \mathbb{R}_-$, $\int_\mathbb{R} \eta = 1$, $\eta \geq 0$ 且在 $\mathbb{R}_-$ 上非减);速度函数 $V_1,V_2 \in W^{3,\infty}_{\mathrm{loc}}(\mathbb{R})$ 满足 $V_1'(\xi) \geq 0$, $V_2'(\xi) \leq 0$ 对所有 $\xi \in [\lambda_{\min}, \lambda_{\max}]$ 成立,其中 $\lambda_{\min} := \min\{ \inf u_0, \min_{r \in [\inf u_0,\sup u_0]} V_2(r) \}$, $\lambda_{\max} := \max\{ \sup u_0, \max_{r \in [\inf u_0,\sup u_0]} V_2(r) \}$;且 $V_2$ 在 $[\inf u_0, \sup u_0]$ 上双 Lipschitz,即存在 $L > 0$ 使得 $\frac{1}{L}|\xi - \xi'| \leq |V_2(\xi) - V_2(\xi')| \leq L|\xi - \xi'|$ 对所有 $\xi,\xi' \in [\inf u_0, \sup u_0]$ 成立。令 $u^\varepsilon$ 为 (1.1) 的唯一经典解(由文献[10]定理 2.1 保证),$u$ 为 (1.3) 的唯一经典解,$T^* > 0$ 为其最大存在时间(即 $u \in C([0,T^*); W^{2,\infty}(\mathbb{R}))$ 且 $\limsup_{t \to T^*} \|u(t,\cdot)\|_{W^{2,\infty}} = \infty$)。则对任意 $T < T^*$,存在常数 $C = C(T, u_0, V_1, V_2, \eta) > 0$,使得

$$ \sup_{t \in [0,T]} \|u^\varepsilon(t,\cdot) - u(t,\cdot)\|_{L^\infty(\mathbb{R})} \leq C \varepsilon. $$

该误差阶 $O(\varepsilon)$ 是最优的,源于 $\|U * \eta^\varepsilon - U\|_{L^\infty} = O(\varepsilon)$(由 $U \in W^{1,\infty}$ 及核一阶矩有限性导出,见脚注1),且证明中所有估计均未损失 $\varepsilon$ 阶。此外,定理要求 $V_2$ 的单调递减性($V_2' \leq 0$)与 $V_1$ 的单调递增性($V_1' \geq 0$)共同确保复合函数 $V_1 \circ V_2$ 的单调性,这是控制特征线方向与避免反向耦合失稳的关键;而双 Lipschitz 条件则将 $u^\varepsilon$ 与 $u$ 的 $L^\infty$ 收敛性等价于 $V_2(u^\varepsilon)$ 与 $V_2(u)$ 的收敛性,从而绕过直接处理 $u^\varepsilon$ 的非线性耦合困难。

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

证明的核心障碍在于:双重非局部性导致标准能量估计失效——直接对 $u^\varepsilon - u$ 做 $L^2$ 或 $L^\infty$ 估计时,非线性项 $\partial_x(V_1(V_2(u^\varepsilon)*\eta^\varepsilon)u^\varepsilon) - \partial_x(V_1(V_2(u))u)$ 展开后出现 $V_2(u^\varepsilon)*\eta^\varepsilon - V_2(u)$ 与 $u^\varepsilon - u$ 的强耦合,而 $V_2(u^\varepsilon)*\eta^\varepsilon$ 的正则性低于 $u^\varepsilon$,无法用 Sobolev 嵌入控制。作者的根本创新点是引入变量代换 $U^\varepsilon := V_2(u^\varepsilon)$,将原方程转化为关于 $U^\varepsilon$ 的新方程(3.1),其非局部性仅存在于 $V_1(U^\varepsilon * \eta^\varepsilon)$ 一项,从而将问题降维为单非局部性估计。具体而言,步骤1构造 $U^\varepsilon$ 的演化方程:利用链式法则 $\partial_t U^\varepsilon = V_2'(u^\varepsilon) \partial_t u^\varepsilon$ 并代入 (1.1),得到

$$ \partial_t U^\varepsilon + V_2'(u^\varepsilon) u^\varepsilon \partial_x(V_1(U^\varepsilon * \eta^\varepsilon)) + V_1(U^\varepsilon * \eta^\varepsilon) V_2'(u^\varepsilon) \partial_x u^\varepsilon = 0. $$

类似地,$U = V_2(u)$ 满足 (3.2)。步骤2–3 对 $|U^\varepsilon - U|$ 的 $L^\infty$ 范数进行微分不等式估计:在 $|U^\varepsilon - U|(t,\cdot)$ 的最大值点 $x_{\max}$ 处,利用 $\partial_x |U^\varepsilon - U| = 0$ 消去部分导数项,将 $\frac{d}{dt} \|U^\varepsilon - U\|_{L^\infty}$ 拆分为 $I + J$(见(3.4))。关键突破在于对 $I$ 的精细分解(3.5):$I = I_1 + I_2$,其中 $I_1$ 控制 $V_2'(u)u - V_2'(u^\varepsilon)u^\varepsilon$ 的差异,利用 $V_2$ 的双 Lipschitz 性将其控制为 $C \|U^\varepsilon - U\|_{L^\infty}$;$I_2$ 进一步拆为 $I_2^1 + I_2^2$(3.9),其中 $I_2^1$ 处理 $V_1(U) - V_1(U * \eta^\varepsilon)$ 的局部近似误差,通过 $\|U - U * \eta^\varepsilon\|_{L^\infty} = O(\varepsilon)$(因 $U \in W^{1,\infty}$ 且 $\int |y| \eta(y) dy < \infty$)导出 $O(\varepsilon)$ 项;而 $I_2^2$ 利用各向异性核的支撑性质($\operatorname{supp}(\eta) \subset \mathbb{R}_-$)和 $V_1', V_2'$ 的符号条件($V_1' \geq 0, V_2' \leq 0$)证明其第二项非正(见(3.12)末行推导),从而仅需估计第一项,最终得到 $I_2 \leq O(\varepsilon) + C \|U^\varepsilon - U\|_{L^\infty}$。步骤3对 $J$ 的估计(3.17)同样依赖 $\|U - U^\varepsilon * \eta^\varepsilon\|_{L^\infty} \leq \|U - U * \eta^\varepsilon\|_{L^\infty} + \|(U - U^\varepsilon) * \eta^\varepsilon\|_{L^\infty} = O(\varepsilon) + \|U^\varepsilon - U\|_{L^\infty}$。综上,得到 Grönwall 型微分不等式 $\frac{d}{dt} \|U^\varepsilon - U\|_{L^\infty} \leq O(\varepsilon) + C \|U^\varepsilon - U\|_{L^\infty}$,积分后得 $\|U^\varepsilon - U\|_{L^\infty} = O(\varepsilon)$。最后,由 $V_2$ 的双 Lipschitz 性(1.5)立即推出 $\|u^\varepsilon - u\|_{L^\infty} \leq L \|U^\varepsilon - U\|_{L^\infty} = O(\varepsilon)$。为处理一般 $\eta$(未必 Lipschitz),作者设计逼近引理 2.3:取 Lipschitz 核列 $\{\eta_n\}$ 满足 $\|\eta_n - \eta\|_{L^1} \to 0$,证明对应的解 $u^{\varepsilon,n} \to u^\varepsilon$ 一致收敛,从而将一般核情形归约为 Lipschitz 核情形。整个框架的枢纽在于变量代换 $U = V_2(u)$ 与各向异性核的定向支撑特性($\mathbb{R}_-$)的协同作用,后者在估计 $I_2^2$ 时消除了一个潜在的坏项,这是此前工作(如[11]依赖指数核微分性质)所不具备的普适技术。

4. 与之前工作的比较

与既有工作相比,本文在多个维度实现实质性推进。首先,相较于文献[10](提出模型并建立适定性),本文首次严格回答了其明确提出的奇异极限问题(摘要首句及引言末段),且适用范围远超其单调初值限制。其次,对比文献[11](在 $V_1 = \mathrm{Id}$、指数核 $\eta(x) = e^x \mathbf{1}_{(-\infty,0]}$ 下对熵解证明极限),本文的推进体现在三方面:(i) 假设放宽:[11] 严重依赖 $\eta' = \eta$ 导致的 $\partial_x(\eta^\varepsilon * f) = \eta^\varepsilon * \partial_x f + \frac{1}{\varepsilon} \eta(0^-) f$ 恒等式,而本文仅需各向异性条件 (I),涵盖所有非减、左支撑、单位质量的 $L^1 \cap L^\infty$ 核(如截断线性核、B-spline 核等),显著扩大适用类;(ii) 模型推广:[11] 处理 $V_1 = \mathrm{Id}$ 的特例,本文允许 $V_1,V_2$ 为一般非线性函数,且要求 $V_1' \geq 0$, $V_2' \leq 0$,覆盖更广的交通流速度-密度关系(如 Greenshields 模型 $V_2(u) = 1-u$, $V_1(v) = v$);(iii) 方法论差异:[11] 基于熵条件与 Kružkov 技巧,目标是全局弱解;本文采用经典解框架,在 $T < T^*$(激波形成前)获得 $O(\varepsilon)$ 误差阶,精度更高。再与文献[7,5](研究 $V_2 = \mathrm{Id}$ 的单非局部情形)比较:[7] 在 BV 核下证明 $L^1$ 弱收敛,但未给出速率;[5] 发现某些核会导致 TV 爆炸,表明紧性不总成立;本文通过 $V_2$ 双 Lipschitz 性与各向异性核结合,规避 TV 问题,直接获得 $L^\infty$ 一致收敛及定量误差。代价方面:本文结果局限于经典解($T < T^*$),未解决熵解的全局极限(作者明确声明“留待未来工作”);且要求 $V_2$ 严格单调($V_2' \leq 0$)与双 Lipschitz,排除了 $V_2$ 为常数或平坦区间的可能性,而[10]的单调性方法对此无要求。此外,本文未处理多维情形(文献[8]已拓展单非局部情形至多维),亦未讨论 $\varepsilon \to 0$ 后激波的渐近结构(文献[2]研究 Oleinik 型估计)。综上,本文以精巧的变量代换与各向异性核分析,为双重非局部守恒律建立了首个普适、定量、经典解层面的奇异极限理论,填补了从模型提出([10])到极限验证之间的关键空白。


5. 晨辉云模型强解的整体适定性

Global well-posedness of strong solutions to a model for the morning glory cloud

作者: Jinkai Li, Yuan Ma, Dong Wang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13706v1

摘要翻译: 本文研究 Constantin-Johnson \cite{CaJo} 最近导出的一个描述对流层中非线性波传播(特别是晨辉云现象)的模型的强解整体适定性。假设初始速度 $v_0\in H^1$ 且热力学强迫项 $K\in L^2(0,T;L^2)$,我们证明该系统初边值问题存在唯一的整体强解。在现有文献中,类似结果仅对充分小的初始数据成立。

Comment: 无

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

1. 研究对象和背景

本文研究由Constantin-Johnson于2023年导出的晨光云(morning glory cloud)大气动力学简化模型,其核心方程组为:

$$ \begin{cases} \partial_t v + v\partial_x v + w\partial_z v - \mu \Delta v + \alpha v + \beta w + K = 0, \\ \partial_x v + \partial_z w = 0, \end{cases} \tag{1.1} $$

其中定义域为 $\Omega = \mathbb{T} \times (0,1)$(水平方向周期、垂直方向单位区间),未知量为水平速度 $v = v(x,z,t)$ 与垂直速度 $w = w(x,z,t)$;$\mu > 0$ 为粘性系数,$\alpha, \beta \in \mathbb{R}$ 为线性阻尼与耦合常数,$K = K(x,z,t)$ 为给定的热力学强迫项(如辐射加热或潜热释放)。由不可压缩约束 $\partial_x v + \partial_z w = 0$ 及边界条件 $v|_{z=1} = \partial_z v|_{z=0} = w|_{z=0} = 0$(见(1.2)),可显式消去 $w$:

$$ w(x,z,t) = -\int_0^z \partial_x v(x,\xi,t)\,d\xi. \tag{1.5} $$

代入后,(1.1) 等价于关于 $v$ 的单一方程:

$$ \partial_t v - \mu \Delta v + v\partial_x v - \left(\int_0^z \partial_x v\,d\xi\right) \partial_z v + \alpha v - \beta \int_0^z \partial_x v\,d\xi + K = 0, $$

具有非局部非线性项 $-\left(\int_0^z \partial_x v\,d\xi\right) \partial_z v$ 与积分型线性耦合项。该模型物理上描述对流层中大尺度非线性重力波的传播,尤其刻画澳大利亚晨光云所呈现的横跨地平线的管状滚轴云结构——其动力学本质是二维化的(因垂直尺度远小于水平尺度),但又不等价于标准二维Navier-Stokes或原始方程(PEs):关键区别在于缺失压力项 $\partial_x p$,导致无法施加双侧垂直边界条件 $w|_{z=0}=w|_{z=1}=0$(否则系统过定);而只能施加单侧 $w|_{z=0}=0$,另一侧 $w|_{z=1}$ 自由。这一结构差异使能量估计中出现致命的边界项 $\int_{\mathbb{T}} w|_{z=1} |\partial_z v|^2 dx$(见第3页分析),在PEs中因 $w|_{z=1}=0$ 而自动消失,但在本模型中必须被控制。已有工作包括:Matioc-Roberti (2021) 建立了 $v_0 \in L^2$ 下的整体弱解与局部强解;Alonso-Orán & Granero-Belinchón (2022) 在 $\|v_0\|_{L^\infty}$ 充分小时证明了整体强解存在;但对任意大的 $H^1$ 初值,整体强解的存在性仍是空白。本文旨在填补这一关键空白:在仅要求 $v_0 \in H^1(\Omega)$(周期且满足 $v_0|_{z=1}=0$)及 $K \in L^2_{\mathrm{loc}}([0,\infty); L^2(\Omega))$ 的条件下,证明强解的全局存在唯一性。该问题的重要性在于:它代表了一类具有物理真实性但数学结构脆弱的“准原始方程”模型的适定性突破——其非局部非线性、不对称边界、缺失压力项共同构成三重障碍,是检验现代PDE分析技术边界的理想测试平台。

2. 主要定理或主要结果

定理1.1(全局强解存在唯一性):设 $\Omega = \mathbb{T} \times (0,1)$,$K \in L^2_{\mathrm{loc}}([0,\infty); L^2(\Omega))$,初值 $v_0 \in H^1(\Omega)$ 满足周期性及 $v_0|_{z=1} = 0$。则系统(1.1)–(1.4)存在唯一的全局强解 $v$,即对任意 $T \in (0,\infty)$,有

$$ v \in C([0,T]; H^1(\Omega)) \cap L^2(0,T; H^2(\Omega)), \quad \partial_t v \in L^2(0,T; L^2(\Omega)). $$

该解满足:(i) $v$ 在 $\Omega \times (0,T)$ 上几乎处处满足(1.1),其中 $w$ 由(1.5)定义;(ii) 满足初始条件(1.4)及边界条件(1.2)–(1.3)。进一步,该结论在两种推广情形下仍成立(见Remark 1.1):(i) 将 $\partial_z v|_{z=0} = 0$ 替换为 $v|_{z=0} = 0$;(ii) 定义域改为 $\Omega = \mathbb{R} \times (0,1)$。证明依赖于一系列关键先验估计,其常数 $C$ 均仅依赖于 $\|v_0\|_{H^1}, \|w_0\|_{L^2}, \alpha, \beta, \mu, \|K\|_{L^2(0,T;L^2)}$ 和 $T$,且不依赖于 $T$ 的大小(从而支持全局性):

  • $L^\infty_t L^2_xL^2_z$ 能量估计(Prop. 4.1):
    $$ \sup_{0\le t\le T} \|v\|_2^2 + \int_0^T \|\nabla v\|_2^2 \,dt \le C. $$
  • $L^\infty_t L^4$ 与 $L^2_t W^{1,2}(v^2)$ 估计(Prop. 4.2):
    $$ \sup_{0\le t\le T} \|v\|_4^4 + \int_0^T \||v|\nabla v\|_2^2 \,dt \le C. $$
  • $w$ 的 $L^\infty_t L^2$ 与 $L^2_t H^1$ 估计(Prop. 4.3):
    $$ \sup_{0\le t\le T} \|w\|_2^2 + \int_0^T \|\nabla w\|_2^2 \,dt \le C. $$
  • $H^1$ 模与 $L^2_t H^2$ 模估计(Prop. 4.4):
    $$ \sup_{0\le t\le T} \|\nabla v\|_2^2 + \int_0^T \|\nabla^2 v\|_2^2 \,dt \le C. $$

其中关键尺度关系体现在Prop. 4.4的证明中:通过对方程(1.1)乘以 $-\Delta v$ 并分部积分,得到

$$ \frac{1}{2} \frac{d}{dt} \|\nabla v\|_2^2 + \mu \|\Delta v\|_2^2 \le C \|w\|_2^2 \|\nabla w\|_2^2 \|\nabla v\|_2^2 + \text{低阶项}, $$

再应用Gronwall不等式,其指数因子为 $\exp\left( \int_0^t \|w\|_2^2 \|\nabla w\|_2^2 \,ds \right)$,而由Prop. 4.3知 $\|w\|_2^2 \|\nabla w\|_2^2 \in L^1(0,T)$,故该因子有界,从而闭环估计。

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

证明的核心障碍是获取 $L^\infty(0,T; H^1)$ 先验估计。标准原始方程(PEs)技术(如对 $\partial_z v$ 能量估计)在此失效,因其在分部积分 $\int_\Omega w \partial_z v \partial_{zz} v \,dx dz$ 时产生边界项 $\int_{\mathbb{T}} w|_{z=1} |\partial_z v|^2 dx$,而 $w|_{z=1}$ 无界控(因仅知 $w|_{z=0}=0$)。作者的真正创新点在于提出一个两阶段非标准估计框架,绕过直接控制 $\partial_z v$ 的困难:

第一阶段:对 $v$ 本身实施高阶 $L^4$ 估计(Prop. 4.2)。将方程(1.1)乘以 $|v|^2 v$,利用引理2.1中精细的Gagliardo-Nirenberg插值:

$$ \|v^3\|_{L^2_x L^1_z} \le C \left( \|v \partial_x v\|_2^{1/4} + \|v\|_4^{1/2} \right) \|v\|_4^{3/2} \left( \|\partial_x v\|_2^{1/4} + \|v\|_2^{1/4} \right) \|v\|_2^{3/4}, $$

以及 $\|v^3\|_2 \le C \left( \|v \nabla v\|_2^{1/2} + \|v\|_4 \right) \|v\|_4^2$,将非线性项 $\int w \partial_z v |v|^2 v$ 控制为 $\|\partial_x v\|_2 \|v^3\|_{L^2_x L^1_z}$,最终导出 $\|v\|_4^4$ 的微分不等式,结合Prop. 4.1的能量估计,得到 $L^\infty_t L^4 \cap L^2_t W^{1,2}(v^2)$ 控制。

第二阶段:推导并分析 $w$ 的演化方程(4.1)。由(1.5)及(1.1)直接计算得:

$$ \partial_t w - \mu \Delta w + 2 \int_0^z \partial_x(v \partial_x v)\,d\xi + \partial_x(w v) + \alpha w + \beta \int_0^z \partial_x w\,d\xi + \int_0^z \partial_x K\,d\xi = 0, $$

配以边界条件 $w|_{z=0}=0$, $\partial_z w|_{z=1}=0$。此方程的关键在于其非线性项 $2 \int_0^z \partial_x(v \partial_x v)\,d\xi$ 和 $\partial_x(w v)$ 均可通过 $v$ 的 $L^4$ 与 $H^1$ 范数控制。将(4.1)乘以 $w$,利用 $\|w\|_{L^4_x L^\infty_z} \le C \|\partial_x w\|_2$(由Sobolev嵌入)及前述 $v$ 的估计,得到 $\|w\|_2^2$ 的微分不等式,从而获得 $L^\infty_t L^2 \cap L^2_t H^1$ 对 $w$ 的控制(Prop. 4.3)。

第三阶段:闭环 $H^1$ 估计。将(1.1)乘以 $-\Delta v$,关键处理项 $\int_\Omega w \partial_z v \Delta v \,dx dz$。利用 $w$ 的 $L^2$ 和 $H^1$ 估计及 $\partial_z v$ 的 $L^2$ 估计,通过Hölder与插值:

$$ \left| \int_\Omega w \partial_z v \Delta v \,dx dz \right| \le \|w\|_4 \|\partial_z v\|_4 \|\Delta v\|_2 \le C \|w\|_2^{1/2} \|\nabla w\|_2^{1/2} \|\partial_z v\|_2^{1/2} \|\Delta v\|_2^{3/2}, $$

再用Young不等式吸收至耗散项 $\mu \|\Delta v\|_2^2$,余项由 $\|w\|_2^2 \|\nabla w\|_2^2 \|\nabla v\|_2^2$ 主导,而该乘积属于 $L^1(0,T)$(由Prop. 4.3),故Gronwall引理给出 $\|\nabla v\|_2^2$ 的一致上界。整个框架的枢纽是 $w$ 的独立演化方程(4.1)及其能量估计,这本质上是将原系统的非局部耦合转化为一个可独立控制的辅助变量,是区别于PEs分析的标志性技术。

4. 与之前工作的比较

本文与既有工作的本质推进体现在三个维度:

1. 初值大小假设的彻底解除:此前唯一已知的整体强解结果是Alonso-Orán & Granero-Belinchón (2022) 的小初值定理,要求 $\|v_0\|_{L^\infty} < \varepsilon_0$($\varepsilon_0$ 依赖于 $\mu,\alpha,\beta$)。本文完全取消此限制,仅需 $v_0 \in H^1$(含任意大的 $L^2$ 或 $H^1$ 范数),实现了从“小数据理论”到“大数据理论”的跃迁。代价是:必须采用更精细的 $L^4$ 估计与 $w$ 的辅助方程框架,而非直接的小参数展开。

2. 边界条件结构的严格匹配:Matioc-Roberti (2021) 在相同方程(1.1)下建立了弱解,但其设定为无限通道 $\mathbb{R} \times I$ 且边界条件为 $v|_{z=0}=w|_{z=0}=v|_{z=1}=0$,这隐含了 $w|_{z=1}$ 的某种衰减性;而本文在紧致化水平方向 $\mathbb{T} \times (0,1)$ 下,采用更物理的混合边界:$v|_{z=1}=\partial_z v|_{z=0}=w|_{z=0}=0$,此设定下 $w|_{z=1}$ 完全自由,导致前述边界项障碍。本文的估计框架(特别是 $w$ 的演化方程(4.1)及对应边界 $\partial_z w|_{z=1}=0$)正是为此非对称边界定制,而Matioc-Roberti的弱解框架无法提供强解所需的 $H^2$ 正则性。

3. 与原始方程(PEs)的对比揭示结构性鸿沟:Remark 1.2明确指出,若在(1.1)中添加压力项 $\partial_x p$ 并补全双侧 $w$ 边界条件(即Holst-Rademacher (2020) 的修改模型),则问题退化为经典2D PEs,其全局强解早已解决。本文恰恰反向工作:在移除压力项这一关键物理简化下,证明其仍保持全局适定性。这表明,晨光云模型的全局性并非源于PEs的压力梯度平衡机制,而是由其特定的非局部耦合 $\beta w$ 与阻尼 $\alpha v$ 协同实现。局限性在于:本文结果尚未推广至三维或带温度方程的完整晨光云模型;且 $K \in L^2_{\mathrm{loc}}$ 的要求比某些PEs结果(允许 $K \in L^p$)略强;此外,解的长时间行为(如衰减率、吸引子)未被探讨。从提供内容中无法确认是否改进了弱解的正则性提升阈值,或与Constantin-Johnson原始行波解的稳定性分析存在联系。


6. 超临界空间中无粘与完全耗散Boussinesq系统的范数膨胀

Norm Inflation for Inviscid and Fully Dissipative Boussinesq Systems in Supercritical Spaces

作者: Qionglei Chen, Yaowei Xie

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13694v1

摘要翻译: 我们证明了二维Boussinesq系统在超临界Besov空间中的范数膨胀现象,即强不适定性。对于无粘系统,范数膨胀发生在 $\dot B^β_{p,q}(\mathbb R^2)\times \dot B^β_{p,r}(\mathbb R^2)$ 中,其中 $β\neq0$,$1<p\leq\infty$,$1\leq q,r\leq\infty$,且 $-2<β-\frac{2}{p}<1$。对于完全耗散系统,相同结论成立的范围为 $-2<β-\frac{2}{p}<-1$。在两种情形下,结果覆盖了几乎所有满足局部可积性条件的超临界Besov空间。范数膨胀出现在密度分量 $ρ$ 中,而速度分量 $u$ 保持有界。在完全耗散情形下,膨胀所发生的函数空间对 $u$ 是超临界的,但对 $ρ$(相对于其自身标度)却是次临界的。这并非矛盾:密度由处于超临界区域的速度场输运,而正是这一输运机制导致了 $ρ$ 中的范数膨胀。

Comment: 35 pages

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

1. 研究对象和背景

本文研究二维Boussinesq系统在超临界齐次Besov空间中的强不适定性,核心对象是如下初值问题:

$$ \begin{cases} \partial_t u - \mu \Delta u + u \cdot \nabla u + \nabla p = \rho e_2, \\ \partial_t \rho - \nu \Delta \rho + u \cdot \nabla \rho = 0, \\ \mathrm{div}\, u = 0, \\ u(x,0) = u_0(x),\quad \rho(x,0) = \rho_0(x), \end{cases} \tag{1.1} $$

其中 $u:\mathbb{R}^2\times\mathbb{R}_+\to\mathbb{R}^2$ 为速度场,$\rho:\mathbb{R}^2\times\mathbb{R}_+\to\mathbb{R}$ 为密度(或温度)标量,$p$ 为压力,$e_2=(0,1)$,$\mu,\nu\ge0$ 分别为粘性与扩散系数。本文聚焦两个极端情形:(i) 无粘情形 $\mu=\nu=0$;(ii) 完全耗散情形 $\mu=\nu=1$。该系统是流体力学中描述大气锋面、海洋环流、Rayleigh–Bénard对流与层化流动的经典耦合模型,其数学本质在于速度与标量通过输运项 $u\cdot\nabla\rho$ 和浮力项 $\rho e_2$ 相互作用,形成非线性反馈回路。

适定性理论的分水岭由尺度不变性决定。对无粘系统,其自然尺度为欧拉型:$(u_\lambda,\rho_\lambda,p_\lambda)(t,x):=(\lambda u(\lambda t,\lambda x),\, \lambda \rho(\lambda t,\lambda x),\, \lambda^2 p(\lambda t,\lambda x))$,对应临界正则性为 $u_0\in \dot B^{1+2/p}_{p,1},\, \rho_0\in \dot B^{1+2/p}_{p,1}$($1<p<\infty$),此即文献[8]所证局部适定性空间;而超临界区域定义为 $\beta - 2/p < 1$,即 $\dot B^\beta_{p,q}\times \dot B^\beta_{p,r}$ 中 $\beta$ 低于临界指数 $1+2/p$ 的部分。对完全耗散系统,其自然尺度为抛物型:$u_\lambda(t,x)=\lambda u(\lambda^2 t,\lambda x),\, \rho_\lambda(t,x)=\lambda^3 \rho(\lambda^2 t,\lambda x),\, p_\lambda(t,x)=\lambda^2 p(\lambda^2 t,\lambda x)$,导出速度临界空间为 $\dot B^{-1+2/p}_{p,q}$,密度临界空间为 $\dot B^{-3+2/p}_{p,r}$(见(1.2)式)。因此,同一 $\dot B^\beta_{p,q}\times \dot B^\beta_{p,r}$ 空间对 $u$ 超临界当且仅当 $\beta - 2/p < -1$,而对 $\rho$ 则可能处于亚临界(因 $\rho$ 临界指数为 $-3+2/p$,故 $\beta > -3+2/p$ 即亚临界)。本文揭示的关键现象正是:即使 $\rho$ 在自身尺度下亚临界,只要其被一个超临界速度场输运,仍可触发范数膨胀。

现有不适定性理论存在明显断层。在临界空间,Bourgain–Li[3]对二维Euler方程在 $H^2$ 中建立范数膨胀;Elgindi–Masmoudi[18]及Bianchini–Hientzsch–Iandoli[2]将此类结果推广至无粘Boussinesq系统的 $L^\infty$(即 $\omega,\nabla\rho\in L^\infty$)临界框架,但机制局限于 $\partial_x\rho$ 或 $\omega$ 的增长。在超临界空间,Luo[24]对三维Euler方程在 $H^s$($0<s<5/2$)中实现范数膨胀;Córdoba–Martínez-Zoroa–Ożański[14]对二维Euler在 $H^s$($1<s<2$)中构造瞬时正则性损失解。然而,对Boussinesq系统,超临界不适定性长期空白:既无针对无粘情形的系统性结果,更无对完全耗散情形的处理。尤其困难在于,Euler型膨胀依赖于速度自身的涡量拉伸,而Boussinesq中膨胀必须由耦合机制驱动——即 $u$ 的低正则性导致 $\rho$ 输运失稳。本文旨在填补这一空白,精确刻画超临界Besov空间中范数膨胀发生的完整范围,并阐明其物理根源:非退化剪切(nondegenerate shear)诱导的相位振荡放大。

2. 主要定理或主要结果

本文在两类情形下严格证明了范数膨胀(norm inflation),即强不适定性:存在任意小初值,其解在极短时间内范数爆炸。所有结论均针对 $\mathbb{R}^2$ 上的初值问题(1.1),并明确给出参数范围、初值正则性、爆炸时间尺度与范数下界。

定理1.1(无粘情形):设 $1<p\le\infty$,$\beta\neq0$ 满足 $-2<\beta-2/p<1$。则对任意 $\varepsilon>0$,存在光滑紧支初值 $(u_0,\rho_0)\in C_c^\infty(\mathbb{R}^2)$ 满足

$$ \|u_0\|_{\dot B^\beta_{p,1}} + \|\rho_0\|_{\dot B^\beta_{p,1}} < \varepsilon, $$

且存在时间 $0<t^*<\varepsilon$ 使得

$$ \|\rho(t^*)\|_{\dot B^\beta_{p,\infty}} > \frac{1}{\varepsilon}. $$

此结论对所有 $1\le q,r\le\infty$ 成立,因 $\dot B^\beta_{p,1}\hookrightarrow \dot B^\beta_{p,q}$ 且 $\|f\|_{\dot B^\beta_{p,r}}\ge\|f\|_{\dot B^\beta_{p,\infty}}$。

定理1.2(完全耗散情形):设 $1<p\le\infty$,$\beta\neq0$ 满足 $-2<\beta-2/p<-1$。则对任意 $\varepsilon>0$,存在光滑紧支初值 $(u_0,\rho_0)\in C_c^\infty(\mathbb{R}^2)$ 满足

$$ \|u_0\|_{\dot B^\beta_{p,1}} + \|\rho_0\|_{\dot B^\beta_{p,1}} < \varepsilon, $$

且存在时间 $0<t^*<\varepsilon$ 使得

$$ \|\rho(t^*)\|_{\dot B^\beta_{p,\infty}} > \frac{1}{\varepsilon}. $$

同样对所有 $1\le q,r\le\infty$ 成立。

关键技术细节:(i) 爆炸时间尺度为 $t^*=\lambda^{-2/p-1+\beta}\log\log\lambda$,其中大参数 $\lambda\gg1$ 控制初值大小;(ii) 初值构造中含对数修正因子 $(\log\log\lambda)^{-|\beta|/(2+|\beta|)}$,这是为平衡不同正则性阶下的尺度效应所必需;(iii) 当 $\beta<0$ 时,初值 $\rho_0$ 需经 $k=-\lfloor\beta\rfloor+2$ 阶径向导数预混合(premixed),以利用对偶性证明其小性,而爆炸下界则由 $L^p$ 下界与 $\dot B^1_{p,\infty}$ 上界插值得到;(iv) 速度分量 $u$ 在整个过程中保持有界:$\|u(t)\|_{\dot B^\beta_{p,\infty}}\lesssim \varepsilon$,膨胀仅发生在 $\rho$ 分量;(v) 参数范围 $-2<\beta-2/p$ 保证空间局部可积(locally integrable),使初值为真函数而非广义函数。

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

证明的核心障碍在于:如何在超临界正则性下,从耦合非线性系统中分离并放大输运机制的不稳定性,同时控制所有误差项。作者的真正创新点在于将剪切机制(shear mechanism)与预混合构造(premixed construction)相结合,并建立一种新型的“输运主导”稳定性分析框架,而非依赖传统能量估计。

第一,非退化剪切近似解的构造。作者放弃直接处理原系统,转而构造一个精确满足输运方程的近似解 $(\bar u, \bar \rho)$。关键步骤是:(i) 取径向涡量 $\bar w(r)=\lambda^{2/p+1-\beta}(\log\log\lambda)^{-|\beta|/(2+|\beta|)} f(\lambda r)$,其中 $f\in C_c^\infty((0,\infty))$ 满足 $\int_0^\infty f(r)\,dr=\int_0^\infty r f(r)\,dr=0$(Lemma 2.6),这保证其生成的Biot-Savart速度 $\bar u=u[\bar w]$ 是纯角向的 $\bar u=\bar u_\theta(r)e_\theta$,且角速度 $\Omega_\lambda(r):=\bar u_\theta(r)/r$ 具有非退化径向导数 $|\Omega_\lambda'(r)|\ge c_0>0$ 在 $g$ 的支集上成立;(ii) 初始密度 $\bar\rho_0$ 根据 $\beta$ 符号分两种构造:当 $\beta>0$ 时,$\bar\rho_0(r,\theta)=\lambda^{2/p-\beta}(\log\log\lambda)^{-\beta/(2+\beta)} g(\lambda r)\cos\theta$,其输运解 $\bar\rho(t,r,\theta)=\bar\rho_0(r,\theta-t\Omega_\lambda(r))$ 在 $t^*$ 时刻因相位 $t^*\Omega_\lambda(r)$ 的快速变化而产生高频振荡;当 $\beta<0$ 时,采用预混合:$\bar\rho_0$ 包含 $k$ 阶径向导数作用于振荡相位 $\cos(\theta+t^*\Omega_\lambda(r))$,这使得 $\bar\rho_0$ 在负Besov范数下极小(利用Proposition 2.3对偶性),而 $\bar\rho(t^*)$ 则恢复为低频振荡,从而获得 $\dot B^\beta_{p,\infty}$ 下界。

第二,扰动系统的精细能量估计。令 $\zeta=\bar u-u$,$\xi=\bar\rho-\rho$,则 $(\zeta,\xi)$ 满足以 $\bar u\cdot\nabla \bar w - \partial_1 \bar\rho$(来自涡量方程(2.1))为源项的线性系统。作者的关键突破在于:不采用标准 $L^2$ 能量法,而是利用齐次Besov空间的插值与乘子估计(Lemma 2.2, Lemma 2.5),证明 $\|\zeta(t)\|_{\dot B^\beta_{p,1}}+\|\xi(t)\|_{\dot B^\beta_{p,1}}\lesssim \lambda^{-\delta}$ 对某 $\delta>0$ 成立,直至 $t^*$。这依赖于源项的特殊结构:$\bar u\cdot\nabla \bar w$ 为零(因 $\bar w$ 径向静止),故主源项仅为 $-\partial_1 \bar\rho$,其 $\dot B^{\beta-1}_{p,1}$ 范数被精确控制为 $\lesssim \lambda^{-1}$,远小于初值尺度 $\varepsilon\sim \lambda^{-2/p-1+\beta}$。

第三,耗散项的摄动处理。在完全耗散情形,近似解需添加 $-\Delta \bar u$ 和 $-\Delta \bar \rho$ 项。作者证明这些项作为摄动误差,其 $\dot B^\beta_{p,1}$ 范数为 $\lesssim \lambda^{\beta-2/p-2}(\log\log\lambda)^{-|\beta|/(2+|\beta|)}$,比无粘误差高两阶导数。这直接导致允许的 $\beta$ 范围收缩:由无粘的 $\beta-2/p<1$ 变为耗散的 $\beta-2/p<-1$,完美匹配抛物尺度下速度临界性的提升。

整个框架的统一思想是:将不适定性归结为线性输运算子 $\partial_t + \bar u\cdot\nabla$ 在超临界速度场下的失稳,而所有非线性与耗散效应均被降阶为可控扰动。

4. 与之前工作的比较

本文与既有工作存在本质差异,主要体现在机制、范围、技术与代价四个维度。

机制维度:此前所有Boussinesq不适定性结果均基于线性化围绕稳态(如Elgindi–Masmoudi[18]围绕 $\rho_{\mathrm{eq}}=-y$)或临界正则性失稳(如Bianchini–Hientzsch–Iandoli[2]在 $L^\infty$ 中),其增长源于 $\omega$ 或 $\partial_x\rho$ 的线性不稳定性。本文首次建立纯粹由输运主导的超临界范数膨胀,其驱动力是速度场的非退化剪切 $|\Omega_\lambda'(r)|\ge c_0$,而非任何线性算子的谱性质。这种机制与Euler/SQG中的剪切构造[11–14,19]同源,但首次成功嫁接到耦合系统,并证明 $\rho$ 的膨胀可独立于 $u$ 的爆炸发生。

范围维度:在无粘情形,本文覆盖 $-2<\beta-2/p<1$,而此前最强结果[2]仅限于 $\beta-2/p=0$(即 $L^\infty$ 对应的临界点)。特别地,本文包含负正则性 $\beta<0$,这在以往工作中完全缺失;在完全耗散情形,本文范围 $-2<\beta-2/p<-1$ 与Navier-Stokes的超临界结果[25]一致,但[25]仅处理 $u$ 分量,而本文首次证明 $\rho$ 分量可在其自身亚临界尺度下仍被超临界 $u$ 触发膨胀,这揭示了耦合系统中尺度分离的新现象。

技术维度:此前超临界结果(如Luo[24], Córdoba–Martínez-Zoroa–Ożański[14])依赖于精细的傅里叶局部化与相位校正,但均未处理耦合项。本文引入预混合构造(premixed construction)处理负正则性,这是对[25]中反向能量级联思想的重大发展:不是让解自发产生负正则性,而是主动设计初值使其在特定时刻“解缠”(unwind)预存的振荡,从而获得下界。此外,本文对Biot-Savart算子在Besov空间的估计(Lemma 2.5)被用于精确量化剪切强度,而非仅作有界性断言。

代价维度:本文结果的代价是爆炸仅发生在 $\rho$ 分量,而 $u$ 保持有界,这与Euler方程中 $u$ 自身爆炸不同。同时,初值必须满足支撑分离条件($\mathrm{supp}\,f\cap\mathrm{supp}\,g=\emptyset$)及矩条件($\int f=\int rf=0$),这是为确保 $\bar u$ 的纯角向性与非退化剪切所必需,表明该机制具有几何特异性。从提供内容中无法确认是否可推广至其他边界条件(如周期域 $\mathbb{T}^2$)或分数耗散情形(Remark 1.3提及但未展开证明),亦未讨论三维推广的可能性。


7. 全局分歧与悬垂周期定常水波的构造性存在性

Global Bifurcation and the Constructive Existence of Overhanging Periodic Steady Water Waves

作者: Matthieu Cadiot, Susanna Haziot

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13567v1

摘要翻译: 我们为具有常数涡度的悬垂周期重力水波提供了构造性的存在性证明。通过引入共形映射公式,我们将流体区域参数化以容纳多值的自由表面高度,从而绕过了传统区域压平框架中出现的坐标奇点。我们在固定的宏观物理参数下进行分析,并采用全局分歧框架,辅以基于牛顿–康托罗维奇定理的计算机辅助证明,从而构造性地证明了精确解分支的存在性。

Comment: 无

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

1. 研究对象和背景

本文研究二维稳态、周期性、重力驱动、常涡度 $\gamma \in \mathbb{R}$ 的自由表面水波问题。设波以恒定相速度 $c$ 相对于静止床 $Y = 0$ 传播,在随波移动坐标系中,流体域 $D = \{ (X,Y) \in \mathbb{R}^2 : 0 < Y < \eta(X) \}$ 与自由面 $S = \{ (X, \eta(X)) : X \in \mathbb{R} \}$ 静止。物理变量为相对速度场 $U = (U,V)$、压力 $P$ 及重力常数 $g > 0$。其动力学由定常不可压 Euler 方程控制:

$$ \begin{cases} U_X + V_Y = 0, \\ U U_X + V U_Y = -P_X, \\ U V_X + V V_Y = -P_Y - g, \end{cases} \quad \text{in } D. $$

自由面 $S$ 上满足大气压边界条件 $P = P_{\mathrm{atm}}$ 与运动学条件 $U \cdot N = 0$(即 $S$ 为流线),且全域满足常涡度约束 $V_X - U_Y = \gamma$。由不可压性,引入相对流函数 $\Psi$ 满足 $U = \Psi_Y$, $V = -\Psi_X$;再定义修正流函数 $\zeta(X,Y) \coloneqq \Psi(X,Y) + \frac{\gamma}{2} Y^2$,则 $\Delta \zeta = 0$ 在 $D$ 内成立,且 $U = \zeta_Y - \gamma Y$, $V = -\zeta_X$。因 $S$ 是流线,有 $\zeta - \frac{\gamma}{2} \eta^2 = m$($m$ 为相对质量通量);结合 Bernoulli 积分得动态边界条件:

$$ U^2 + V^2 = Q - 2g \eta \quad \text{on } S, $$

其中 $Q \coloneqq 2(E - P_{\mathrm{atm}} - \gamma m)$ 为修正 Bernoulli 常数。为处理悬垂(overhanging)情形——即自由面非单值函数 $Y = \eta(X)$,而需允许 $X$ 关于弧长参数非单调——作者采用共形映射 $X + iY = \xi(x,y) + i\eta(x,y)$ 将一个物理周期的流体域 $D$ 共形地映到固定矩形域 $\Omega = \{ (x,y) \in \mathbb{R}^2 : -\pi < x < \pi,\, -h < y < 0 \}$,其中 $y = -h$ 对应床 $Y = 0$,$y = 0$ 对应 $S$。$\xi, \eta$ 满足 Cauchy–Riemann 方程 $\xi_x = \eta_y$, $\xi_y = -\eta_x$,故均为调和函数。由此导出等价椭圆边值系统:

$$ \begin{aligned} &\Delta \zeta = 0, \quad \Delta \eta = 0 && \text{in } \Omega, \\ &\zeta - \tfrac{\gamma}{2} \eta^2 - m = 0 && \text{on } y = 0, \\ &(\zeta_y - \gamma \eta \eta_y)^2 - (Q - 2g \eta)(\eta_x^2 + \eta_y^2) = 0 && \text{on } y = 0, \\ &\eta = 0,\; \zeta = 0 && \text{on } y = -h, \end{aligned} $$

并辅以 $2\pi$-周期性、对称性 $\eta(\pm x + 2\pi, y) = \eta(x,y)$、平均水深约束 $\frac{1}{2\pi} \int_{-\pi}^{\pi} \eta(x,0)\,dx = h$ 与零漂移规范 $\frac{1}{2\pi} \int_{-\pi}^{\pi} (\xi(x,0)-x)\,dx = 0$。关键几何适定性条件包括:半周期内严格垂直单调性 $\eta_x(x,0) < 0$ for $x \in (0,\pi)$;全局单射性(防止自交)要求 $0 < \xi(x,0) < \pi$ for $x \in (0,\pi)$;以及非退化性 $\eta_x^2 + \eta_y^2 \ge \varepsilon > 0$ in $\Omega$ 与 $Q - 2g \eta \ge \delta > 0$ on $y = 0$。该问题的核心难点在于:传统拉平域(domain-flattening)方法在悬垂处产生坐标奇点,无法参数化多值自由面;而奇异摄动或小重力展开(如 Hur–Wheeler [34])仅适用于 $g \ll 1$ 的渐近区域,无法在 $O(1)$ 宏观参数下追踪有限深度全局分歧支。Constantin–Strauss–Varvaruca [17] 提出的全局分歧框架虽能容纳翻转几何,但无法确定分歧支的拓扑终结方式——是形成最大高度波、角点奇点,抑或物理自交?本文旨在填补这一空白:在固定宏观参数下,严格证明存在一条从平凡状态 $\eta \equiv 0$ 出发的全局分歧支,其几何演化路径为:图形式波 $\to$ 悬垂波 $\to$ 槽线自交,从而确证 [17] 中关于“槽线自交终结”的猜想。

2. 主要定理或主要结果

定理 1.1(孤立 Ω 形波的存在性):取无量纲重力常数 $g = 1$、常涡度 $\gamma = -5$、平均水深 $h = 2$,对相对质量通量 $m = -0.85$,存在唯一精确解 $(\tilde{Q}, \tilde{\eta}) \in \mathbb{R} \times C^\infty[-\pi,\pi]$ 满足方程组 (1.4),且严格满足非退化条件 (1.7);该解对应一个物理上有效的、严格悬垂的、单射的波形,其自由面在槽线处凸起呈 Ω 形状。定理 1.2(全局分歧与拓扑跃迁):在相同参数下,存在一条连续、无穷可微的精确 $2\pi$-周期水波解分支,参数化为相对质量通量 $m \in [-1, \bar{m}]$,其中局部分歧点 $\bar{m} \approx -0.061$。该分支的物理几何经历三阶段严格验证的拓扑跃迁:(1) 局部分歧与图形式波:当 $m$ 接近 $\bar{m}$ 时,波形为严格单调图形式;(2) 向悬垂波过渡:存在严格有界临界通量 $\tilde{m}_1 \in [-0.698, -0.6707]$,使得对所有 $m \in [-1, \tilde{m}_1)$,波形严格悬垂(即存在 $x^* \in (0,\pi)$ 满足 $\xi_x(x^*,0) < 0$);(3) 非物理解极限:存在第二临界通量 $\tilde{m}_2 \in [-0.8884, -0.8781]$,当 $m < \tilde{m}_2$ 时,物理单射性条件 (1.6) 失效,自由面在槽线 $x = \pi$ 处发生物理自交(即 $\xi(x_1,0) = \xi(x_2,0)$ 且 $\eta(x_1,0) = \eta(x_2,0)$ 对某 $x_1 \ne x_2$ 成立),此时数学分支终结。所有结论均在 Banach 空间 $H^1 = \mathbb{R} \times X_e$ 中以构造性方式确立,其中 $X_e = \{ a = (a_n)_{n\in\mathbb{Z}} \in \ell^1_e : \|a\|_X = \sum_{n\in\mathbb{Z}} (1+|n|)|a_n| < \infty \}$,且误差由 Newton–Kantorovich 定理给出显式 $H^1$-范数界 $\|U - \tilde{U}\|_{H^1} \le \theta$,其中 $\theta$ 为计算机辅助验证所得严格正数。

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

核心障碍在于:在 $O(1)$ 宏观参数下,既要避免拉平域导致的坐标奇点以参数化悬垂,又要排除其他可能的退化终结(如角点、停滞点、多重停滞)以唯一确定槽线自交。作者的真正创新点在于将共形映射参数化、全局分歧理论与计算机辅助证明(CAP)三者深度融合。第一,共形框架消除了坐标奇异性:通过将自由面表示为 $\eta(x,0)$ 与 $\xi(x,0)$,并利用 Cauchy–Riemann 关系 $\xi_x = \eta_y$,将悬垂判定转化为 $\xi_x(x^*,0) < 0$ 这一纯符号条件,完全绕开 $Y = \eta(X)$ 的多值性困境。第二,Fourier-解析公式化将无限维 PDE 降维为序列空间上的零点问题:定义偶 Fourier 序列空间 $X_e$ 与 $\ell^1_e$,将解 $U = (Q,a)$ 映射至非线性算子 $F : H^1 = \mathbb{R} \times X_e \to H^0 = \mathbb{R} \times \ell^1_e$,其中 $F(U) = 0$ 精确编码了 (1.4) 的所有边界条件;关键在于 $X_e$ 是 Banach 代数(Lemma 3.1),保证非线性项 $\eta^2$、$\eta_x^2$ 等在 Fourier 域中可闭合估计。第三,Newton–Kantorovich 定理的构造性应用:对近似解 $\hat{U}$,显式构造线性化算子 $DF(\hat{U})$ 的近似逆 $A$,并严格验证三项收缩条件:$\|A F(\hat{U})\|_{H^1} \le \alpha$,$\|I - A DF(U)\|_{B(H^1)} \le \beta < 1$(对所有 $U$ 在 $\hat{U}$ 的 $r$-邻域内),及 $\alpha \beta \le \frac{1}{2}$;此处 $\alpha, \beta$ 的计算被分解为有限维截断(用区间算术严格计算)与无穷维余项(用解析衰减估计控制),最终归结为有限个多项式不等式验证。第四,几何性质的解析锁定:定理 2.1 利用调和函数的极值原理与 Hopf 边界引理,证明只要初始满足 $\eta_x \le 0$ 且非停滞,则必有 $\eta_x < 0$ 在 $(0,\pi)$ 内严格成立——此“单向陷阱”将单调性作为拓扑不变量,使数值验证只需聚焦于标量条件 $\xi_x > 0$ 与 $Q - 2g\eta > 0$。第五,全局延续框架:基于均匀收缩映射定理(Theorem 6.1),将分支表示为 Chebyshev 插值,利用其谱精度高效追踪长分支,并在每个步长内重复 Newton–Kantorovich 验证,确保整个分支的 $C^\infty$ 正则性与几何完整性。

4. 与之前工作的比较

与 Constantin–Strauss–Varvaruca [17] 相比:[17] 建立了常涡度周期波的全局分歧框架,但其结论为“分歧支必以某种退化方式终结”,未确定具体终结类型;本文则在相同框架下,通过 CAP 精确锁定终结为槽线自交,解决了其公开猜想。与 Hur–Wheeler [34] 相比:[34] 在 $g \ll 1$ 下对零重力 Crapper 解做摄动,得到悬垂波,但其解位于小参数邻域,无法刻画从平凡态出发的有限深度全局支;本文固定 $g = 1$,首次实现 $O(1)$ 参数下的构造性全局追踪。与 Davila–del Pino–Musso–Wheeler [20] 相比:[20] 针对孤立波使用精细拼接法,在小重力下构造悬垂解,其几何依赖于薄颈匹配,适用范围受限于渐近尺度;本文处理周期情形,无需小参数假设,且给出明确的通量阈值 $\tilde{m}_1, \tilde{m}_2$。与 Chen–Walsh–Wheeler [16] 相比:[16] 证明内波会翻转,但使用拉平域方法,无法处理自由面自交;本文共形框架天然兼容自交,且终结点具有明确物理意义(槽线接触)。代价方面:本文方法高度依赖 CAP,计算复杂度远高于纯解析工作,且结果局限于特定参数三元组 $(g,\gamma,h) = (1,-5,2)$;而 [17] 的结论具有一般性,但缺乏几何分辨率。适用范围上,本文证实在此参数下自交是主导终结模式,但文中明确指出“其他物理参数可能导致最大高度波等不同行为”,故其结论非普适,而是对 [17] 猜想在典型参数下的确认。技术层面,本文将几何约束(单调性、单射性)转化为可验证的解析先验,大幅降低数值负担,这是区别于多数 CAP 工作的关键设计。


8. 多维半线性随机偏微分方程在局部Lipschitz系数与粗糙初值下的存在性、唯一性及路径正则性

Existence, Uniqueness, and Pathwise Regularity for Multidimensional Semilinear SPDEs with Locally Lipschitz Coefficients and Rough Initial Data

作者: Jitendra Nath Naik, Lok Pati Tripathi

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

发布日期: 2026-07-15

arXiv ID: 2607.13341v1

摘要翻译: 我们研究由乘性噪声驱动且具有粗糙初值的多维半线性随机演化方程。漂移项与扩散项系数被假定取值于负阶分数 Sobolev 空间,且可能在初始时刻表现出时间奇异性。当这些系数全局 Lipschitz 时,我们的主要结果建立了解的整体适定性及路径依赖的时空正则性;当系数仅为局部 Lipschitz 时,我们证明了极大局部解的存在性、唯一性及路径正则性。该局部 Lipschitz 条件是相对于某一特定的时间加权范数所定义的。这使得我们的理论结果可应用于线性随机偏微分方程,以及具有非全局 Lipschitz 非线性项的模型,例如随机 Burgers 方程、Allen–Cahn 方程、Fisher–KPP 方程、Burgers–Fisher 方程和 Ginzburg–Landau 系统。此外,我们的框架可容纳奇异初值,例如在维度 $d=1$ 下的 Dirac 测度,以及在维度 $d \in \{2,3\}$ 下的非光滑初值。

Comment: 46 pages

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

1. 研究对象和背景

本文研究如下形式的多维半线性随机发展方程(stochastic evolution equation, SEE):

$$ dX(t) + AX(t)\, dt = F(t, X(t))\, dt + G(t, X(t))\, dW(t), \quad t \in (0, T], \\ X(0) = X_0, $$

其中 $A: D(A) \subset H \to H$ 是定义在可分Hilbert空间 $H$ 上的稠定、自伴、正定线性算子,且具有紧逆;$\{S(t) = e^{-tA}\}_{t \ge 0}$ 是其生成的解析 $C_0$-半群;$W(t)$ 是作用于另一可分Hilbert空间 $U$ 上的 $Q$-Wiener过程(可能为柱形),$U_0 = Q^{1/2}(U)$,$HS(U_0, H)$ 表示从 $U_0$ 到 $H$ 的Hilbert–Schmidt算子空间。关键特征在于:初值 $X_0$ 允许取值于负阶分数Sobolev空间 $\dot{H}^\mu := D(A^{\mu/2})$,即 $\mu < 0$,例如一维情形下 $X_0$ 可为Dirac测度 $\delta_0 \in \dot{H}^{-1/2-\varepsilon}(\mathbb{R})$;非线性项 $F: [0,T] \times \Omega \times \dot{H}^\mu \to \dot{H}^{-\alpha}$ 与 $G: [0,T] \times \Omega \times \dot{H}^\mu \to HS(U_0, \dot{H}^{-\beta})$ 在时间零点附近允许出现代数奇性(temporal singularity),其Lipschitz常数被赋予时间权重 $t^{-\hat{\alpha}}$ 与 $t^{-\hat{\beta}}$;且 $F,G$ 仅需满足局部Lipschitz条件(相对于时间加权范数 $\|u\|_{\lambda,t} = \|u\|_{\dot{H}^\mu} + t^\lambda \|u\|_{\dot{H}^\nu}$,其中 $\lambda = (\nu - \mu)/2 \ge 0$),而非全局Lipschitz——这使得模型能容纳物理中根本性的非全局Lipschitz非线性,如Burgers方程中的对流项 $u \partial_x u$、Allen–Cahn方程中的三次项 $u - u^3$、Fisher–KPP方程中的反应项 $u(1-u)$ 等。该问题的重要性源于:一方面,经典理论(如Da Prato–Debussche、Jentzen–Röckner、van Neerven等)要求初值至少属于 $H = \dot{H}^0$ 或更高正则空间,无法处理 $\mu < 0$ 的奇异初值;另一方面,现有局部Lipschitz框架(如Seidler [41]、Brzeźniak [6]、Agresti–Veraar [1,2])虽能处理多项式非线性,但其解空间依赖于初值正则性 $\mu$ 与非线性映射空间指标 $\nu$ 的耦合约束 $\mu \ge \nu$,从而排除了 $\mu < 0$ 且 $\nu > 0$ 的关键情形(例如 $X_0 = \delta_0 \in \dot{H}^{-1/2}$,而 $F(u) = u^2$ 需映入 $\dot{H}^{-1}$,故 $\nu = 1$, $\mu = -1/2$)。因此,本文旨在填补这一结构性空白:构建一个统一框架,同时兼容粗糙初值($\mu < 0$)与局部Lipschitz非线性(尤其是含导数的非线性),并建立其路径正则性(pathwise regularity)——即对几乎每个样本轨道 $\omega$,解 $X(\cdot, \omega)$ 在时空上具有Hölder连续性,而非仅在 $L^p(\Omega)$ 意义下。这一目标直接关联到后续数值分析(如Wang–Wang [48])的可行性,以及对白噪声驱动下奇异初值演化(如Gyöngy [24] 在1D Burgers方程中的开创性工作)的严格数学刻画。

2. 主要定理或主要结果

本文核心结果分为全局与局部两类,均建立在参数约束(Assumption 2.2)之上:$0 \le \nu - \mu < 1$,$\alpha < 2 - \nu$,$\beta < 1 - \nu$,$0 \le \hat{\alpha} < 1 - \frac{1}{2} \max\{\nu - \mu, \alpha + \nu\}$,$0 \le \hat{\beta} < \frac{1}{2} - \frac{1}{2} \max\{\nu - \mu, \beta + \nu\}$。

全局结果(Theorem 2.1):设 $X_0 \in L^p(\Omega; \dot{H}^\mu)$,$p \ge 2$,且 $F,G$ 满足全局Lipschitz条件(Assumption 2.4),即存在常数 $L_j$ 使得对所有 $u,v \in \dot{H}^\nu$,

$$ \|F(t,u) - F(t,v)\|_{\dot{H}^{-\alpha}} \le L_1 t^{-\hat{\alpha}} \|u-v\|_{\dot{H}^\nu}, \quad \|G(t,u) - G(t,v)\|_{HS(U_0, \dot{H}^{-\beta})} \le L_3 t^{-\hat{\beta}} \|u-v\|_{\dot{H}^\nu}. $$

则存在唯一 mild 解 $X \in V_p$,其中解空间 $V_p = \{ \phi \in C([0,T]; L^p(\Omega; \dot{H}^\mu)) : t^\lambda \phi \in C_b((0,T]; L^p(\Omega; \dot{H}^\nu)) \}$,范数为 $\|\phi\|_{V_p} = \sup_{t \in (0,T]} (\|\phi(t)\|_{L^p(\Omega; \dot{H}^\mu)} + t^\lambda \|\phi(t)\|_{L^p(\Omega; \dot{H}^\nu)})$,$\lambda = (\nu - \mu)/2$。且满足先验估计 $\|X\|_{V_p} \le C(1 + \|X_0\|_{L^p(\Omega; \dot{H}^\mu)})$;对初值连续依赖:$\|X-Y\|_{V_p} \le C \|X_0 - Y_0\|_{L^p(\Omega; \dot{H}^\mu)}$;时间正则性:对任意 $a \in (0,T)$,有 $X \in C^{0,\delta}([a,T]; L^p(\Omega; \dot{H}^\nu))$,其中 Hölder 指数 $\delta < \min\{1 - \frac{1}{2}(\nu + \alpha)_+, \frac{1}{2}(1 - (\nu + \beta)_+)\}$;若 $X_0 \in L^p(\Omega; \dot{H}^{\mu+\varepsilon})$,$0 < \varepsilon \le 1$,则 $X \in C^{0,\delta_\varepsilon}([0,T]; L^p(\Omega; \dot{H}^\mu))$,且 $\delta_\varepsilon < \min\{\varepsilon/2, 1 - \hat{\alpha} - \lambda - \frac{1}{2}(\mu + \alpha)_+, \frac{1}{2} - \hat{\beta} - \lambda - \frac{1}{2}(\mu + \beta)_+\}$。

路径正则性(Theorem 2.2):若 $p > 2$,则解具有更强的路径性质:$X \in L^p(\Omega; C([0,T]; \dot{H}^\mu))$,当 $(\mu + \alpha)_+ < 2 - 2/p$ 且 $(\mu + \beta)_+ < 1 - 2/p$;对任意 $a > 0$,$X \in L^p(\Omega; \mathring{C}^{0,\delta}([a,T]; \dot{H}^\nu))$,其中 $\delta < \min\{1 - 1/p - \frac{1}{2}(\nu + \alpha)_+, \frac{1}{2} - 1/p - \frac{1}{2}(\nu + \beta)_+\}$;若 $X_0 \in L^p(\Omega; \dot{H}^{\mu+\epsilon})$,则 $X \in L^p(\Omega; \mathring{C}^{0,\delta}([0,T]; \dot{H}^\mu))$,其中 $\delta < \min\{\epsilon/2, 1 - \hat{\alpha} - \lambda - \frac{1}{2}(\mu + \alpha)_+, 1 - 1/p - \frac{1}{2}(\mu + \alpha)_+, \frac{1}{2} - \hat{\beta} - \lambda - \frac{1}{2}(\mu + \beta)_+, \frac{1}{2} - 1/p - \frac{1}{2}(\mu + \beta)_+\}$。

局部结果(Theorem 2.4):设 $X_0 \in L^0(\Omega; \dot{H}^\mu)$,且 $F,G$ 满足局部Lipschitz条件(Assumption 2.5),即对每个 $R > 0$,存在 $L_j(R)$ 使得对所有 $u,v$ 满足 $\|u\|_{\lambda,t}, \|v\|_{\lambda,t} \le R$,有

$$ \|F(t,u) - F(t,v)\|_{\dot{H}^{-\alpha}} \le L_1(R) t^{-\hat{\alpha}} \|u-v\|_{\lambda,t}, \quad \|G(t,u) - G(t,v)\|_{HS(U_0, \dot{H}^{-\beta})} \le L_3(R) t^{-\hat{\beta}} \|u-v\|_{\lambda,t}. $$

则存在唯一极大局部 mild 解 $(X, \tau_{\max})$,其中 $\tau_{\max}: \Omega \to [0,T]$ 是停时,且存在局部化序列 $\tau_n \uparrow \tau_{\max}$ a.s.,使得对每个 $n$,截断解 $X(\cdot \wedge \tau_n)$ 满足:$X(\cdot \wedge \tau_n) \in C([0,T]; \dot{H}^\mu) \cap \mathring{C}^{0,\delta}([a,T]; \dot{H}^\nu)$,$0 < \delta < \min\{1 - \frac{1}{2}(\nu + \alpha)_+, \frac{1}{2} - \frac{1}{2}(\nu + \beta)_+\}$;且 $\lim_{t \searrow 0} \|X(t)\|_{\lambda,t} = \|X_0\|_{\dot{H}^\mu}$ a.s. on $\{\tau_n > 0\}$;爆炸准则为 $\limsup_{t \uparrow \tau_{\max}} \|X(t)\|_{\lambda,t} = +\infty$ a.s. on $\{\tau_{\max} < T\}$。

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

证明的核心障碍在于:如何在 $\mu < 0$(即初值极度粗糙)与局部Lipschitz非线性(导致解可能在有限时间内爆破)双重约束下,构造一个兼具存在性、唯一性与路径正则性的解框架。作者的真正创新点体现在三个相互支撑的技术层面:

第一,解空间与时间加权范数的精巧设计。标准方法(如Seidler [41])将解空间设为 $C([0,T]; \dot{H}^\nu)$ 或 $C^{0,\delta}([0,T]; \dot{H}^\nu)$,这隐含要求初值 $X_0 \in \dot{H}^\nu$,从而强制 $\mu \ge \nu$。本文彻底解耦初值正则性 $\mu$ 与非线性作用空间指标 $\nu$,引入时间加权范数 $\|u\|_{\lambda,t} = \|u\|_{\dot{H}^\mu} + t^\lambda \|u\|_{\dot{H}^\nu}$($\lambda = (\nu - \mu)/2$)作为局部Lipschitz条件的度量基准,并将解空间 $V_p$ 定义为 $\{\phi: t^\lambda \phi \in C_b((0,T]; L^p(\Omega; \dot{H}^\nu))\}$。此设计的关键在于:当 $t \to 0^+$ 时,$t^\lambda \to 0$,故 $\|X(t)\|_{\lambda,t} \to \|X_0\|_{\dot{H}^\mu}$,自然容纳 $\mu < 0$;而当 $t$ 远离零时,$t^\lambda$ 提供正则性增益,使 $X(t)$ 获得 $\dot{H}^\nu$ 正则性。这使得局部Lipschitz常数 $L_j(R)$ 仅需控制在球 $\{u: \|u\|_{\lambda,t} \le R\}$ 内,而该球在 $t=0$ 附近退化为 $\dot{H}^\mu$ 中的球,在 $t>0$ 时则包含 $\dot{H}^\nu$ 信息,从而规避了 $\mu \ge \nu$ 的硬性限制。

第二,分数卷积估计的实质性改进(Lemma 3.3)。证明解的正则性高度依赖于半群卷积算子 $\int_0^t S(t-\sigma) f(\sigma) \, d\sigma$ 和 $\int_0^t S(t-\sigma) g(\sigma) \, dW(\sigma)$ 的估计。经典估计(如Hong–Liu [26, Prop. 4.1])在边界 $t=0$ 处损失Hölder正则性。本文引理3.3给出了更精细的分数阶卷积估计:对任意 $\eta, \rho \in \mathbb{R}$,存在常数 $C$ 使得

$$ \left\| \int_0^t A^\eta S(t-\sigma) \phi(\sigma) \, d\sigma \right\|_{\dot{H}^\gamma} \le C \sup_{\sigma \in (0,t]} \sigma^{\rho} \|\phi(\sigma)\|_{\dot{H}^\theta}, \quad \text{只要 } \gamma - \theta < 1 + \eta - \rho. $$

此估计消除了边界损失,并适用于全实轴 $\eta, \rho \in \mathbb{R}$,是论证 $X \in \mathring{C}^{0,\delta}([a,T]; \dot{H}^\nu)$ 的基石。它直接支撑了定理2.1(iii)与2.2(ii)中Hölder指数的最优性。

第三,局部化与停时构造的严格化。为避免Seidler [41]等工作中因错误使用 $S(t \wedge \tau_n - \sigma)$ 导致的适应性(adaptedness)失效(见Remark 2.3(ii)),本文采用基于随机区间 $J0, \tau_nK$ 的指示函数 $\mathbf{1}_{J0, \tau_nK}(\sigma)$ 构造截断方程:

$$ \mathbf{1}_{J0, \tau_nK}(t) X(t) = \mathbf{1}_{J0, \tau_nK}(t) \left( S(t) X_0 + \int_0^t \mathbf{1}_{J0, \tau_nK}(\sigma) S(t-\sigma) F(\sigma, X(\sigma)) \, d\sigma + \int_0^t \mathbf{1}_{J0, \tau_nK}(\sigma) S(t-\sigma) G(\sigma, X(\sigma)) \, dW(\sigma) \right). $$

此形式确保所有被积函数均为 $\mathcal{F}_\sigma$-可料(predictable),从而保证Itô积分良定义。这一技术细节是证明极大局部解存在性与唯一性的逻辑支点。

4. 与之前工作的比较

本文与既有文献的比较可从四个维度展开:

1. 与全局Lipschitz框架(Andersson et al. [4])的比较:Andersson等人建立了 $X_0 \in L^p(\Omega; \dot{H}^\mu)$,$\mu < 0$,且 $F,G$ 全局Lipschitz时的解理论,但其框架要求 $F,G$ 映入 $\dot{H}^0$(即 $\alpha = \beta = 0$),且时间权重 $\hat{\alpha}, \hat{\beta}$ 的范围较窄。本文通过引入更一般的负阶空间 $\dot{H}^{-\alpha}, \dot{H}^{-\beta}$ 和更宽泛的参数约束(Assumption 2.2),显著拓展了适用范围。例如,对一维随机Burgers方程 $du + (-\partial_x^2)u \, dt = u \partial_x u \, dt + \sigma(u) \, dW$,其非线性 $F(u) = u \partial_x u$ 从 $\dot{H}^\nu$ 映入 $\dot{H}^{-1}$,故需 $\alpha = 1$;而Andersson框架无法容纳此情形。本文Remark 2.2(i)明确指出,[4]的Theorem 2.9是本文Theorem 2.1在 $\nu = 0, \mu = -2\max\{\delta,0\}, \alpha = 2\tilde{\alpha}, \beta = 2\tilde{\beta}$ 下的特例,凸显了本文参数自由度的提升。

2. 与局部Lipschitz框架(Seidler [41], Brzeźniak [6])的比较:Seidler在Hilbert空间中证明了局部Lipschitz系数下极大局部解的存在性,但其解空间为 $C([0,T]; H^\delta)$,要求初值 $X_0 \in H^\zeta$ 且 $\zeta \ge \delta$(即 $\mu \ge \nu$)。本文Remark 2.4明确指出,这是其框架的根本局限,而本文通过时间加权范数 $\|\cdot\|_{\lambda,t}$ 实现了解耦,使得 $\mu < \nu$ 成为可能。例如,对一维随机Allen–Cahn方程 $du + (-\partial_x^2)u \, dt = (u - u^3) \, dt + \sigma(u) \, dW$,初值 $X_0 = \delta_0 \in \dot{H}^{-1/2}$($\mu = -1/2$),非线性 $F(u) = u - u^3$ 从 $\dot{H}^1$ 映入 $\dot{H}^{-1}$($\nu = 1, \alpha = 1$),此时 $\mu = -1/2 < \nu = 1$,Seidler框架完全失效,而本文可处理。

3. 与高正则性理论(Jentzen–Röckner [29])的比较:Jentzen与Röckner利用半群的平滑效应,证明了当 $X_0 \in H$ 且 $F,G$ 全局Lipschitz时,解具有 $C^{0,\delta}([0,T]; H^\gamma)$ 正则性($\gamma > 0$)。本文Theorem 2.1(iii)与Theorem 2.2(iii)在 $\mu = 0$ 时可还原为此结果(Remark 2.2(ii)),但本文的贡献在于将此高正则性推广至 $\mu < 0$ 的粗糙初值情形,并给出了精确的Hölder指数范围,揭示了正则性增益与初值粗糙度 $\mu$、非线性强度 $\alpha,\beta$ 及时间权重 $\hat{\alpha},\hat{\beta}$ 的定量关系。

4. 与1D奇异初值理论(Gyöngy [24])的比较:Gyöngy针对1D随机Burgers方程,证明了Dirac初值下解的存在性与Hölder正则性,但其方法高度依赖于一维特定结构(如Burgers方程与热方程的联系)和白噪声的特殊性质。本文提供了一个抽象、多维、适用于广泛SPDEs(包括Burgers、Allen–Cahn、Ginzburg–Landau等)的统一框架,并明确指出其结果在 $d=1$ 时“complementing foundational 1D results such as those by Gyöngy [24]”,即本文并非取代而是补充和推广了此类具体模型的结果。代价在于,本文对噪声的要求更一般(trace-class noise for $d=2,3$),且其正则性结论是路径意义下的($\mathring{C}^{0,\delta}$),而非Gyöngy中可能更精细的逐点估计。


9. 含超线性源项的时间非局部发展方程的解的存在性与非存在性

The existence and nonexistence for time nonlocal evolution equations with superlinear sources

作者: Xi Huang, Li Peng

类别: math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13177v1

摘要翻译: 我们研究具有 $ (\mathcal{PC}) $ 型核的时间非局部发展方程初值问题的可解性行为。我们的目标是分析若干充分条件,以保证当非线性项呈现快速增长时,其温和解的局部存在性与可积性;此外,我们也建立了不可解结果的一个充分条件。结果表明,可解性行为与初始值的指标密切相关,并出现临界维数现象。证明依赖于从属原理(subordination)与单调迭代法(monotone iterative method)。最后,我们给出若干例子以说明该结果的广泛适用性。

Comment: 无

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

1. 研究对象和背景

本文研究如下时间非局部扩散型演化方程的初值问题:

$$ \begin{cases} \partial_t \big[ k * (u - u_0) \big](t,x) - \Delta u(t,x) = f(u(t,x)), & t > 0,\ x \in \mathbb{R}^N, \\ u(0,x) = u_0(x), & x \in \mathbb{R}^N, \end{cases} \tag{1.2} $$

其中卷积定义为 $(k * v)(t) = \int_0^t k(t-\theta) v(\theta)\, d\theta$,核对 $(k,l)$ 满足完全正性条件(PC):$k,l \in L^1_{\mathrm{loc}}(\mathbb{R}_+)$ 非负、非增,且满足 $k * l \equiv 1$ 在 $(0,\infty)$ 上几乎处处成立。该类核涵盖大量重要模型:当 $k(t) = t^{\alpha-1}/\Gamma(\alpha)$($\alpha \in (0,1)$)时,(1.2) 退化为时间分数阶扩散方程 $\partial_t^{\alpha} u - \Delta u = f(u)$;当 $k(t) = \int_0^1 t^{\alpha-1}/\Gamma(\alpha)\, d\mu(\alpha)$ 时,对应分布式阶分数阶方程。物理上,该方程建模强记忆效应下的反常输运过程,如高度无序多孔介质中的溶质迁移、粘弹性材料中具有长程时间关联的热传导等。

非线性项 $f$ 满足超线性快速增长条件(Af):$f \in C^1[0,\infty)$,$f(u)>0$,$f'(u)>0$,且其逆函数的积分收敛:

$$ G(u) := \int_u^\infty \frac{1}{f(\tau)}\, d\tau < \infty, \quad \forall u > 0, $$

并进一步要求渐近临界指标:

$$ \lim_{u \to \infty} f'(u) G(u) = 1. $$

该条件精确刻画了超多项式增长行为,典型例子包括 $f(\tau) = e^{\tau^p}$($p \ge 1$)与 $f(\tau) = e^{e^\tau}$。值得注意的是,$G(u)$ 是 $f$ 的“爆破时间函数”——若考虑常微分方程 $u' = f(u)$,其解在有限时间内爆破当且仅当 $G(u_0) < \infty$;而此处 $G(u_0)^{-\rho} \in L^1_{\mathrm{ul}}(\mathbb{R}^N)$ 成为关键可解性判据。

已有工作主要聚焦于线性情形($f \equiv 0$ 或 $f$ 线性):Vergara–Zacher 利用能量方法与新型积分微分不等式得到有界域上解的尖锐衰减估计;Kemppainen 等借助 Fourier 乘子法与松弛函数理论获得 $L^r$ 估计,并首次指出解的大时间行为存在“临界维数现象”;Pozo–Vergara 则利用从属原理构造基本解,结合调和分析揭示解的渐近结构。然而,对于超线性(尤其是超多项式)非线性源项,现有文献极度匮乏:文献 [20] 仅针对时间分数阶情形(即特定 $k$)使用单调迭代法建立局部正解存在性,其证明严重依赖 Wright 函数 $\Phi_\alpha$ 的精细渐近与高斯热核的逐点控制;而对一般 (PC) 核,既无系统性存在/不存在判据,亦无适用于超指数增长的泛函框架。本文旨在填补这一空白:建立适用于任意 (PC) 核、覆盖 $e^{u^p}$ 及 $e^{e^u}$ 等极端增长的统一存在性与不存在性理论,揭示初始数据 $u_0$ 的“临界指数” $\rho$ 如何通过 $G(u_0)^{-\rho}$ 的可积性决定解的局部可解性,并证实该现象本质上是 $N$ 维空间中记忆核 $(1*l)(t)$ 的尺度行为与非线性增长速率之间精细平衡的结果。

2. 主要定理或主要结果

本文核心结果分为存在性与不存在性两部分,均在 $L^1_{\mathrm{ul}}(\mathbb{R}^N)$ 框架下建立。

存在性(定理 4.1):设 $N \ge 1$,$(k,l) \in (\mathrm{PC})$,$l$ 非增,$u_0 \ge 0$,且 $f$ 满足 (Af)。若存在 $\rho > N/2$ 使得 $G(u_0)^{-\rho} \in L^1_{\mathrm{ul}}(\mathbb{R}^N)$,则初值问题 (1.2) 在某个 $T > 0$ 上存在非负整体温和解 $u \in L^1_{\mathrm{ul}}(\mathbb{R}^N)$,且满足:

  • $u(t,\cdot) \in L^1_{\mathrm{ul}}(\mathbb{R}^N)$ 对所有 $t \in (0,T)$;
  • $f(u(t,\cdot)) \in L^1_{\mathrm{ul}}(\mathbb{R}^N)$ 对几乎每个 $t \in (0,T)$;
  • 温和解公式成立:
    $$ u(t,x) = Z(t,\cdot) *_x u_0 + \int_0^t Y(t-\theta,\cdot) *_x f(u(\theta,\cdot))\, d\theta, $$
    其中传播函数 $Z(t,x) = -\int_0^\infty H(\tau,x)\, w(t,d\tau)$,$Y(t,x) = \int_0^\infty H(\tau,x)\, \eta(t,d\tau)$,$H$ 为 $N$ 维高斯热核,$w,\eta$ 由核 $l$ 的拉普拉斯变换 $\hat{l}(\lambda)$ 通过 $\hat{w}(\lambda,d\tau) = (\lambda \hat{l}(\lambda))^{-1} e^{-\tau \hat{l}(\lambda)^{-1}} d\tau$ 等确定。

不存在性(定理 4.2,文中未完整给出但摘要与引言明确声明):在相同 (PC) 假设下,若 $G(u_0)^{-\rho} \notin L^1_{\mathrm{ul}}(\mathbb{R}^N)$ 对所有 $\rho \le N/2$ 成立,则问题 (1.2) 不存在任何非负温和解。该临界阈值 $\rho_c = N/2$ 构成“临界维数现象”:当 $N=1$ 时 $\rho_c = 1/2$,$N=2$ 时 $\rho_c = 1$,依此类推。特别地,对 $f(u)=e^{u^p}$,$G(u) \sim u^{1-p} e^{-u^p}/p$,故 $G(u_0)^{-\rho} \in L^1_{\mathrm{ul}}$ 等价于 $u_0(x)^{\rho(p-1)} e^{\rho u_0(x)^p} \in L^1_{\mathrm{ul}}$,凸显初始数据必须以特定指数衰减才能避免瞬时爆破。

所有结论均依赖于记忆核的尺度函数 $(1*l)(t)$:其小时间行为 $\sim t^\beta$($\beta>0$)直接决定 $L^p_{\mathrm{ul}}$-$L^q_{\mathrm{ul}}$ 估计中的衰减速率 $[(1*l)(t)]^{-N/2(1/p-1/q)}$,进而将 $\rho$ 与空间维数 $N$ 耦合为 $\rho > N/2$。

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

证明的核心障碍在于:(1) 一般 (PC) 核 $k$ 不具备显式表达式(如分数阶核的幂律),导致经典格林函数或 Mittag-Leffler 函数技术失效;(2) 超线性源项 $f(u)$ 的增长远超多项式,使得标准 $L^p$ 或 Sobolev 估计无法控制非线性项;(3) 温和解公式中卷积核 $Z,Y$ 本身由抽象测度 $w,\eta$ 定义,其逐点估计极为困难。

作者的真正创新点在于构建了一套基于“广义逆函数 $G$ 的加权凸性”与“松弛函数尺度分析”的全新框架:

第一,$G$-加权凸变换:摒弃传统 $u^p$ 或 $e^u$ 作为比较函数,引入凸函数 $K(u) = G(u)^{-\rho/c}$($c > \max\{\rho,1\}$)。关键观察是:由 $\lim_{u\to\infty} f'(u)G(u)=1$ 可推出 $K$ 凸、$f/K$ 单调增、且 $K^{-1}$ 凹。这使得 Jensen 不等式可应用于 $Z$ 和 $Y$ 的加权平均(引理 3.6),从而将非线性项 $f(u)$ 控制为 $K(u)$ 的幂次:$f(K^{-1}(v)) \lesssim v^{c(\kappa+1)/\rho}$。

第二,松弛函数尺度嵌入:利用引理 3.4 与 3.5,将 $Z(t,\cdot)*_x v$ 的 $L^q_{\mathrm{ul}}$ 范数控制为 $[(1*l)(t)]^{-N/2(1/p-1/q)} \|v\|_{L^p_{\mathrm{ul}}}$ 加一个趋于零的因子 $C^*(t)$。此估计不依赖 $k$ 的具体形式,仅需 $(k,l)\in (\mathrm{PC})$ 及 $l$ 非增。特别地,在小时间 $t\in(0,T_0)$,$(1*l)(t) \gtrsim t^{\beta}$($\beta>0$),故 $[(1*l)(t)]^{-N/2(1/p-1/q)}$ 提供了必需的奇异性抑制。

第三,单调迭代的 $G$-标度构造:定义上解 $\phi(t) = K^{-1}\big((\delta+1) Z(t,\cdot)*_x K(u_1)\big)$,其中 $u_1 = \max\{u_0,C\}$。利用 $K^{-1}$ 的凹性及中值定理(式 (4.2)),将 $\phi(t) - Z(t,\cdot)*_x u_0$ 下控为 $c\delta/(\rho(\delta+1)) \cdot f(\phi(t)) G(\phi(t))$;再结合 $f/K$ 单调性与 $Y$ 的积分性质(式 (4.5)-(4.6)),将非线性积分项上控为 $\int_0^t l(t-\theta) f(\phi(\theta))\, d\theta$,最终利用 $l$ 非增性及 $(1*l)$ 的积分恒等式证明该积分 $\lesssim [(1*l)(t)]^{1-N(\kappa+1)/(2\rho)} \to 0$(因 $N(\kappa+1)/(2\rho) < 1$)。这迫使 $\phi$ 成为有效上解,从而启动单调迭代(引理 3.1)。

整个框架将非线性增长、记忆核尺度、空间维数三者通过 $G$ 函数与 $\rho$ 指标精密耦合,彻底摆脱了对特定核(如分数阶)或特定非线性(如幂函数)的依赖。

4. 与之前工作的比较

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

1. 核类推广性:此前工作 [20] 仅处理 $k(t)=t^{\alpha-1}/\Gamma(\alpha)$ 这一特殊分数阶核,其证明依赖 Wright 函数 $\Phi_\alpha$ 的尖锐渐近 $\Phi_\alpha(-t) \sim t^{-(1+\alpha)/(2+2\alpha)} e^{-c t^{\alpha/(\alpha+1)}}$ 及其与高斯核的卷积恒等式。本文将核类推广至全部 (PC) 类,涵盖分布式阶核 $k(t)=\int_0^1 t^{\alpha-1}/\Gamma(\alpha)\, d\mu(\alpha)$、指数衰减核 $k(t)=e^{-t}$ 等,其证明不涉及任何特殊函数,仅依赖松弛函数 $s_\mu,r_\mu$ 的抽象性质及 $w,\eta$ 测度的 Laplace 变换结构(引理 2.1)。

2. 非线性强度:文献 [8,11,19,22] 主要处理 $f(u)=u^{1+\sigma}$($\sigma>0$)或 $f(u)=e^u$,其存在性依赖 $\sigma < 2/N$ 或 $\|u_0\|_{L^\infty}$ 小性。本文处理 $f(u)=e^{u^p}$($p\ge1$)乃至 $f(u)=e^{e^u}$,其增长远超指数,传统 $L^\infty$ 或 $L^p$ 方法完全失效。本文通过 $G(u)$ 的临界指标 $\rho$ 将初始数据可积性 $G(u_0)^{-\rho} \in L^1_{\mathrm{ul}}$ 作为判据,实现了对任意超多项式增长的统一刻画。

3. 临界现象本质:Kemppainen 等 [14] 发现线性方程解的衰减存在临界维数,但未与非线性关联。本文首次证明:非线性问题的可解性临界点 $\rho_c = N/2$ 直接源于记忆核尺度 $(1*l)(t)$ 与 $N$ 维空间中高斯核的扩散尺度 $t^{-N/2}$ 的共振。具体而言,引理 3.4 中的 $[(1*l)(t)]^{-N/2(1/p-1/q)}$ 项是 $\rho > N/2$ 的根源;若 $N$ 增大,相同 $u_0$ 需更强衰减(更大 $\rho$)才能满足条件,这与物理直觉一致:高维空间中扩散更快,需更强初始衰减以抗衡超线性增长。

4. 代价与局限:本方法的代价是牺牲了对解的精细正则性(如 Hölder 连续性)或爆破速率的定量刻画;它仅保证 $L^1_{\mathrm{ul}}$ 框架下的存在性,未涉及 $C^\alpha$ 或 $W^{1,p}$ 正则性提升。此外,不存在性结论目前仅对非负解建立,未覆盖变号解情形;且对 $l$ 的非增性假设(用于引理 3.3 的 $b_\delta(t)$ 估计)是否可削弱,文中未讨论。从提供内容中无法确认作者是否对 $\rho = N/2$ 的临界情形进行了分析,亦未显式展开与 [21,17] 中能量方法或 Fourier 乘子法的定量比较。


了解领域论文(摘要翻译)

1. 从大初值散射行为中恢复薛定谔方程的非线性项

Recovery of Schrödinger nonlinearity from the large-data scattering behavior

作者: Yi Sun

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.14080v1

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

摘要翻译

中文标题: 从大初值散射行为中恢复薛定谔方程的非线性项

摘要翻译: 本文研究如下三维非线性薛定谔方程的散射与反散射问题:

$$ \mathrm{i}\partial_t u + Δu = a(x)|u|^2u。 $$

当 $a(x)$ 属于某一适当函数类时,该方程整体适定且具有大初值散射理论。本文证明:大初值散射行为唯一确定了非均匀系数 $a(x)$。


2. 一千零一阶高阶几何网络流

One thousand and one higher-order geometric flows of networks

作者: Luca Benatti, Alessandra Pluda

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13894v1

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

摘要翻译

中文标题: 一千零一阶高阶几何网络流

摘要翻译: 我们证明了作为涉及曲率导数的高阶泛函的 $L^2$-梯度流所导出的一大类曲线与网络的高阶几何流在重参数化意义下的短时存在性与唯一性。此外,我们借助能量方法证明了一个长时存在性结果。


3. Fokker–Planck 半群的短时参数构造及其应用

Short-time parametrix for the Fokker–Planck semigroup and applications

作者: Paul Alphonse, Jean-Marc Bouclet, Matthieu Léautaud

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13814v1

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

摘要翻译

中文标题: Fokker–Planck 半群的短时参数构造及其应用

摘要翻译: 我们在欧氏空间中构造了 Fokker–Planck 半群的短时参数(parametrix)。在可能的应用中,我们得到了该半群的光滑性与局部化性质、半群的近似短时极分解(polar decomposition)的推导、预解式(resolvent)参数的构造、拟谱估计(pseudospectral estimates),以及 Fokker–Planck 算子特征值个数的渐近估计。作为证明中的一步,我们引入了一类算子,称之为亚扇形算子(subsectorial operators),Fokker–Planck 算子属于此类,并描述了其若干泛函分析与谱性质。


4. Boltzmann 方程在中等软势下的定量混沌传播

Quantitative propagation of chaos for the Boltzmann equation with moderately soft potentials

作者: Chenguang Liu, Liping Xu, An Zhang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13825v1

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

摘要翻译

中文标题: Boltzmann 方程在中等软势下的定量混沌传播

摘要翻译: 我们研究与非截断碰撞核的齐次 Boltzmann 方程相关联的 Kac 粒子系统,其中势函数属于中等软势情形 $-1<γ<0$。我们证明了该粒子系统的解律唯一性,并特别地,在初始数据 $f_0$ 具有有限 Fisher 信息和多项式矩的假设下,建立了 Kac 粒子系统的经验测度到 Boltzmann 方程弱解在 Wasserstein-2 距离下的显式收敛速率。据我们所知,这是软势情形下 Kac 粒子系统的首个定量收敛速率结果。


5. 一维线性波动方程带局部分布阻尼的严格Lyapunov函数

Strict Lyapunov function for the one-dimensional linear wave equation with a locally distributed damping

作者: Jean-Michel Coron

类别: math.AP, math.OC

发布日期: 2026-07-15

arXiv ID: 2607.13811v1

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

摘要翻译

中文标题: 一维线性波动方程带局部分布阻尼的严格Lyapunov函数

摘要翻译: 我们为带局部分布阻尼和Dirichlet边界条件的一维线性波动方程构造了一个显式的严格Lyapunov函数。


6. 从部分数据重构对流-扩散方程的系数

Recovery of coefficients for a convection-diffusion equation from partial data

作者: Liam Buisson

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13778v1

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

摘要翻译

中文标题: 从部分数据重构对流-扩散方程的系数

摘要翻译: 本文致力于研究反问题:从狄利克雷激励所生成的通量的部分边界测量中,确定一个抛物型方程中依赖于时间和空间变量的零阶与一阶系数。更准确地说,我们证明了从相应抛物型方程的、仅在边界上被位于区域闭包外部某点所照亮的部分上获取的偏狄利克雷-诺伊曼映射(partial Dirichlet-to-Neumann map),唯一确定一个依赖于时间的对流项和势能项。我们的主要目标是将迄今仅限于椭圆情形的若干观察推广至抛物型方程,并利用这些性质来重构一类一般性的、同时依赖于时间与空间变量的系数。为此,我们构造了一类适用于该抛物型方程的特殊解,引入了一族具有非线性权函数的新卡尔曼估计(Carleman estimates),并将这些工具与X射线变换的部分数据结果相结合。本工作的主要困难在于:与现有抛物型方程部分数据反问题文献不同,此处特殊解的相位函数及相应卡尔曼估计的权函数均为非线性。作为分析的副产品,我们建立了一类具有非线性权函数的抛物型方程卡尔曼估计,该估计可视为椭圆情形下极限卡尔曼权(limiting Carleman weights)概念的抛物类比。


7. 面向非线性Dirac系统的旋量适配几何方法及其在Minkowski时空邻域上一类张量波-Dirac系统中的应用

A spinor-adapted geometric approach for nonlinear Dirac systems and its application to a tensorial wave-Dirac system near Minkowski spacetime

作者: Seokchang Hong

类别: math.AP, math-ph

发布日期: 2026-07-15

arXiv ID: 2607.13775v1

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

摘要翻译

中文标题: 面向非线性Dirac系统的旋量适配几何方法及其在Minkowski时空邻域上一类张量波-Dirac系统中的应用

摘要翻译: 我们在 $(1+3)$ 维渐近平坦时空上研究一类非线性张量波-Dirac系统,该系统作为半线性模型,其动机来自Maxwell-Dirac系统与Einstein-Dirac系统。本模型旨在分离反称张量场的零几何结构与Dirac方程固有的首阶几何结构之间的相互作用,同时规避Einstein方程的导数损失机制与Maxwell方程的规范结构。 我们的分析在整个非线性论证中保持Dirac方程的首阶性质。Dirac流提供了基本的能量恒等式,而定量的时空估计则由Dirac算子平方所诱导的波动方程获得。将这些要素与积分局部能量衰减估计及 $r^p$-加权能量层级相结合,我们为张量分量与旋量分量建立了耦合能量方法。一个关键观察是:Clifford代数与反称张量场的零分解相容,并排除了最奇异的非线性相互作用。 由此,我们建立了小初值解的全局存在性,并给出了定量的加权能量估计与衰减估计。值得注意的是,将零结构与二进制论证相结合,弱衰减如 $t^{-\frac{1}{2}-δ}$ 即足以获得该系统的非线性稳定性,这体现了 \cite{DHRT} 的精神。 我们预期本文所发展的几何思想可为更一般非线性Dirac系统(包括Maxwell-Dirac方程与Einstein-Dirac方程)的研究提供一个有益的起点。


8. 伪相对论薛定谔方程带对数非线性项的归一化行波解的存在性与极限剖面

Existence and Limiting Profiles of Normalized Travelling Wave Solutions for the Pseudo-Relativistic Schrödinger Equation with Logarithmic Nonlinearity

作者: Pietro d’Avenia, Qihan He, Alessio Pomponio, Lianfeng Yang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13628v1

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

摘要翻译

中文标题: 伪相对论薛定谔方程带对数非线性项的归一化行波解的存在性与极限剖面

摘要翻译: 我们研究如下带对数非线性项的伪相对论薛定谔方程的归一化解的存在性与渐近行为:

$$ (\sqrt{-Δ+m^2 }-m )u+i(v\cdot \nabla )u+λu = u\log|u|^2+|u|^{p-2}u, \qquad \text{in }\mathbb{R}^N, $$

并满足质量约束

$$ \|u\|_2^2=a, $$

其中 $m,a>0$, $2<p\le\frac{2N}{N-1}$ 且 $N\ge 2$, $v\in \mathbb{R}^N$ 为行波速度,满足 $|v|<1$, 而 $λ\in\mathbb{R}$ 作为拉格朗日乘子出现,对应于约束上相应能量泛函的极小值。通过变分方法,我们首先给出了此类极小值存在性与不存在性的完整分类。特别地,在质量临界情形 $p=2+\frac{2}{N}$ 下,我们证明存在一个常数 $a^\ast_v$,它构成存在性的阈值。基于此,我们分析了当 $a$ 从下方趋近 $a^\ast_v$ 时此类极小化子的爆破行为。最后,我们分别在质量次临界情形 $2<p<2+\frac{2}{N}$ 下考虑 $\lim\limits_{n\to\infty}a_n=a_0\in(0,+\infty)$ 且 $\{a_n\}\subset(0,+\infty)$,以及在质量临界情形 $p=2+\frac{2}{N}$ 下考虑 $\lim\limits_{n\to\infty}a_n= a_0\in(0,a^\ast_v)$ 且 $\{a_n\}\subset(0,a^\ast_v)$ 时,该问题极小化子的极限剖面。


9. 奇异相互作用核下动理学McKean–Vlasov随机微分方程的定量混沌传播与涨落

Quantitative Propagation of Chaos and Fluctuations for Kinetic McKean–Vlasov SDEs with Singular Interaction Kernels

作者: Zimo Hao, Xicheng Zhang, Xianliang Zhao

类别: math.PR, math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13379v1

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

摘要翻译

中文标题: 奇异相互作用核下动理学McKean–Vlasov随机微分方程的定量混沌传播与涨落

摘要翻译: 我们证明了一类具有外部漂移项和Kato类奇异相互作用核的退化动理学McKean–Vlasov随机微分方程(SDEs)所对应的粒子系统的定量混沌传播估计与中心极限定理(CLT)。特别地,该相互作用核可属于混合 $L^q_tL^{p_v}_vL^{p_x}_x$-空间,其中满足 $\frac{2}{q}+\frac{3d}{p_x}+\frac{d}{p_v}<1$。对于相应的 $N$-粒子系统,在仅假设初始数据具有熵混沌性(entropic chaoticity)的前提下,我们得到了前 $k$ 个粒子在路径空间上的相对熵上界为 $k/N$ 阶。关键技术工具包括动理学Krylov–Khasminskii估计以及关于经验相互作用场的条件Hilbert空间次高斯估计。对于中心极限定理,我们还证明了有限维投影的Berry–Esseen型误差上界。


10. 径向Sturm-Liouville算子的谱分析与Liouville-Green渐近

Spectral Analysis and Liouville-Green Asymptotics for a Radial Sturm-Liouville Operator with Variable Coefficients

作者: Victor S. Gerasimchuk, Bohdan O. Yevdokymenko, Igor V. Gerasimchuk

类别: math-ph, math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13349v1

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

摘要翻译

中文标题: 径向Sturm-Liouville算子的谱分析与Liouville-Green渐近

摘要翻译: 我们研究一个具有径向对称性的初边值问题,其变量分离导出一个正则自伴Sturm-Liouville问题,该问题的系数显式确定且权函数为正。由此得到的Sturm-Liouville算子由一种特殊几何诱导的系数组$ (p(r), w(r)) $定义,从而产生一种非标准的加权谱结构。相应的径向算子在对应的加权Hilbert空间中具有离散实谱及一组完备的标准正交特征函数系,由此给出演化问题的精确谱表示。我们推导出一个精确的超越谱方程以刻画特征值。为分析高频区域,我们将Liouville-Green变换直接应用于径向Sturm-Liouville方程$ -\frac{d}{dr}\left(p(r)\frac{du}{dr}\right) + q(r)u = \lambda w(r)u $,从而获得特征值与特征函数的显式渐近公式,以及相应的Liouville-Green拟模态(quasimodes)。我们在原算子的精确谱框架内建立了这些拟模态的加权正交性及其渐近完备性。所得渐近谱数据被用于构造原边值问题的近似解,并给出谱重构的误差估计。对精确Sturm-Liouville谱的数值计算结果与Liouville-Green预测高度吻合,从而验证了该渐近量子化公式,并确认了所提谱逼近方法的精度。


11. 状态约束的凸Nash均衡问题与线性双曲偏微分方程的耦合

State constrained convex Nash equilibrium problems coupled with linear hyperbolic PDEs

作者: Marcelo Bongarti, Michael Hintermüller

类别: math.OC, math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13299v1

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

摘要翻译

中文标题: 状态约束的凸Nash均衡问题与线性双曲偏微分方程的耦合

摘要翻译: 我们研究状态约束的凸广义Nash均衡问题(GNEPs)与双曲偏微分方程(PDEs)耦合情形下均衡解的存在性。类似的问题已分别针对状态约束的GNEPs与椭圆型和抛物型PDEs的耦合情形得到处理,但在双曲情形下,集值约束映射的正则性问题成为建立存在性理论的主要障碍。这主要源于策略到状态映射的紧性问题。除存在性外,我们还为一类相当一般的线性双曲PDEs提供了首阶最优性条件,并展示了其在诸如波动方程、广告动力学以及网络上线性化等温Euler系统等应用中的相关性。


12. Fokas-Olver-Rosenau-Qiao 方程在临界 Sobolev 空间中的不适定性

Ill-posedness in the critical Sobolev space for the Fokas-Olver-Rosenau-Qiao equation

作者: Dan-Andrei Geba

类别: math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13259v1

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

摘要翻译

中文标题: Fokas-Olver-Rosenau-Qiao 方程在临界 Sobolev 空间中的不适定性

摘要翻译: 本文证明了 Fokas-Olver-Rosenau-Qiao 方程(一类具有三次非线性的修正 Camassa-Holm 型方程)在临界 Sobolev 空间 $H^{5/2}(\bR)$ 中的范数膨胀现象。该结果补全了该方程的适定性理论:此前已知该方程在 $H^{s}(\bR)$ 中对 $s>5/2$ 局部适定。证明依赖于构造满足 Fokas-Olver-Rosenau-Qiao 方程已有爆破判据的显式初值,这一构造步骤本身似乎具有独立意义。


13. 热半群在度量测度空间上的加权可观测性的最优几何障碍

Optimal geometric barriers for weighted observability of heat semigroups on metric measure spaces

作者: Vincent Boulard, Amaury Hayat, Emmanuel Trélat

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

发布日期: 2026-07-14

arXiv ID: 2607.13279v1

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

摘要翻译

中文标题: 热半群在度量测度空间上的加权可观测性的最优几何障碍

摘要翻译: 热方程的加权积分可观测性不等式通常包含一个形如 $e^{-γ/t}$ 的小时间因子。我们证明该尺度并非 Carleman 方法或谱方法的人为产物:它由观测集 $ω$ 的几何结构所强制决定。

设 $A$ 是定义在加倍度量测度空间上有限秩欧氏向量丛截面上的非负自伴算子,满足超收缩性、Davies-Gaffney 估计(即等价于波动方程具有有限传播速度)以及逐点局部 Weyl 律。若在可测集 $ω$ 上对固定时间视界 $T\in(0,+\infty]$ 和容许权函数 $h$ 成立加权积分可观测性不等式,则对任意 $0<κ<\frac{1}{2}$,有

$$ h(t)\leq A_{T,κ}\exp\left(-κ\frac{\mathcal{L}(ω)^2}{t}\right),\qquad 0<t<T, $$

其中 $\mathcal{L}(ω)$ 表示到 $ω$ 的本质极大距离;当 $\mathcal{L}(ω)=+\infty$ 时,该量被替换为任一有限半径。因此,对 $h(t)=e^{-γ/t}$,必有 $γ\geq\mathcal{L}(ω)^2/2$。这以最优阈值解决了早先工作中遗留的关于无穷时间常数的极大距离下界问题。在控制范数约定下,快速控制速率至少为 $\mathcal{L}(ω)^2/4$,从而复现了 Miller 的界。

证明基于谱包 $(\cosh(r\sqrt A)-1)e^{-tA}$。逐点 Weyl 律给出其精确的下增长阶,而有限传播速度与弱核 Kannai 变换公式则使其在 $ω$ 上呈指数小。在不假设核连续性或预解式紧致的前提下,我们发展了逐点谱测度与弱波核理论。该框架涵盖紧黎曼流形上的拉普拉斯型算子、耦合热系统、$\mathbb{R}^d$ 上的薛定谔算子、等正则次拉普拉斯算子与 Grushin 模型,以及度量图上的 $δ'$-耦合拉普拉斯算子。


14. 多孔介质中通量约束流动的对偶框架:分析与数值方法

Duality Framework for Flux Constrained Flow in Porous Media: Analysis and Numerics

作者: Harbir Antil, Keegan L. A. Kirk, Felipe Pérez

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

发布日期: 2026-07-14

arXiv ID: 2607.13256v1

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

摘要翻译

中文标题: 多孔介质中通量约束流动的对偶框架:分析与数值方法

摘要翻译: 我们引入并分析了饱和多孔介质中的达西流,其在边界一部分上受到法向通量的双边约束。该问题被表述为在 $H(\mathrm{div};Ω)$ 的一个凸子集上最大化一个基于速度的对偶凹能量;Fenchel 对偶性识别出一个基于压力的预对偶形式,给出强对偶性,并在受约束边界上提供具有互补结构的凸最优性条件。原始–对偶间隙满足一个后验误差恒等式,该恒等式不含任何通用常数,且对任意容许近似均成立。该对偶结构被 Raviart–Thomas/Crouzeix–Raviart 离散化所继承;由此我们导出了离散误差恒等式,并在解及通量约束具有分数阶正则性假设下,得到先验误差衰减速率。数值实验(包括由局部化原始–对偶间隙指示器驱动的自适应加密)验证了理论结果。


15. 双曲抛物面 $L^p$ Fourier 延拓极大化子的存在性

Existence of maximizers for $L^p$ Fourier extension from the hyperbolic paraboloid

作者: Benjamin Bruce, Diogo Oliveira e Silva, Betsy Stovall

类别: math.CA, math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13181v1

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

摘要翻译

中文标题: 双曲抛物面 $L^p$ Fourier 延拓极大化子的存在性

摘要翻译: 我们证明了极大化子的存在性,并且证明:对于三维环境空间中双曲抛物面相关的 $L^p \to L^q$ Fourier 延拓不等式,其极大化序列存在子列,该子列在模去对称性意义下收敛于极大化子。


16. 具有 Lévy 扩散的非局部演化方程:适定性与极限行为

Non-local evolution equations with Lévy diffusion: Well-posedness and limiting behavior

作者: Xi Huang, Li Peng, Yong Zhou

类别: math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13163v1

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

摘要翻译

中文标题: 具有 Lévy 扩散的非局部演化方程:适定性与极限行为

摘要翻译: 在本注记中,我们关注一类具有 Lévy 扩散的非局部(关于时间)演化方程;这类方程作为单向粘弹性流体流动及具有记忆效应的物理现象的模型而出现。我们首先在所涉记忆核满足完全正性条件的情形下,考虑一个非局部线性演化问题的经典解的存在性。接着,我们研究该模型在 Lévy 扩散指数集中于 $2$ 附近时向广义 Rayleigh-Stokes 方程的极限行为;我们证明:具有 Lévy 扩散的非局部(关于时间)问题的解一致收敛于广义 Rayleigh-Stokes 方程的解,并揭示其收敛速率。最后,我们建立了一个含非线性的非局部演化问题的温和解的存在性及其极限行为。证明基于从属原理(subordination principle)与松弛函数理论(relaxation function theory)。


不相关论文(基本信息)

1. 边界驱动非谐链的对数Sobolev不等式

Log-Sobolev inequalities for boundary-driven anharmonic chains

作者: Jianfeng Lu

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

发布日期: 2026-07-15

arXiv ID: 2607.13953v1

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

摘要: We study the non-equilibrium steady state of a weakly anharmonic chain of $N$ oscillators driven at its boundary by Langevin thermostats at unequal temperatures. Under a perturbative weak-anharmonicity condition, we prove a full-gradient logarithmic Sobolev inequality whose constant is independent of the chain length $N$. For homogeneous pinned chains, an additional quantitative regularity assumption yields a boundary space-time logarithmic Sobolev inequality and relative-entropy decay on the same $O(N^3)$ relaxation time scale as the harmonic chain. The proof extracts a finite-dimensional Gaussian component from the boundary noise and compares conditional terminal-state laws by a change of variables. The estimates are uniform over bounded positive temperatures and require no near-equilibrium assumption on their difference.


2. 具有退化肿瘤扩散的酸介导侵袭模型的传播动力学

Propagation dynamics of an acid-mediated invasion model with degenerate tumor diffusion

作者: Xinyue Cao, Quentin Griette, Xiong Li, Yang Wang

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13944v1

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

摘要: We investigate traveling wave fronts for an acid-mediated tumor invasion model with density-dependent degenerate diffusion. The model is a partially diffusive PDE–ODE system of Gatenby–Gawlinski type, in which the tumor diffusion coefficient $D(U)$ is allowed to be a general decreasing function satisfying $D(1)=0$. This degeneracy causes the traveling wave equation for the tumor component to lose uniform ellipticity near the healthy state, and hence standard arguments for nondegenerate reaction diffusion systems are not directly applicable. To overcome this difficulty, we introduce a nonlinear change of variables which removes the degeneracy from the highest-order term of the tumor equation. For each fixed admissible tumor profile, the acid profile is represented by a Green kernel, while the healthy-tissue profile is obtained from an explicit integral formula. The transformed tumor profile is then constructed as the stationary limit of a uniformly parabolic auxiliary problem. By combining comparison principles, local Schauder estimates, carefully chosen super- and sub-solutions, and Schauder’s fixed point theorem, we prove the existence of traveling wave fronts for every wave speed $θ\ge 2\sqrt{rD(0)}$. The resulting wave connects the tumor-dominant state $(0,1,1)$ at $z=-\infty$ to the healthy state $(1,0,0)$ at $z=+\infty$. We further establish strict pointwise bounds, monotonicity of all wave components, and one-sided exponential asymptotic estimates in both the transformed variable and the original traveling-wave variable.


3. Heal Kernel 与幂零覆盖上的闭测地线渐近性

Heat Kernel and Closed Geodesic Asymptotics for Nilpotent Coverings

作者: Atsushi Katsuda

类别: math.SP, math.AP, math.DG, math.DS, math.NT

发布日期: 2026-07-15

arXiv ID: 2607.13890v1

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

摘要: We establish all order long-time asymptotic expansions for heat kernels on nilpotent coverings and for prime closed geodesics in fixed central classes of nilpotent quotients of compact hyperbolic surfaces. The exact lattice-side input is the finite-dimensional rational Floquet-Bloch theory of the companion paper: rational Kirillov restrictions give exact finite-dimensional fibers, and a generalized Pytlik functional gives exact Fourier-inversion and normalized-trace identities. At a rational parameter $p/q$ the decomposition is exact, and the fluctuation of the fiber integrand is controlled only by $q$. Hence the large-denominator comparison with the smooth Kirillov or Schrödinger normal form is uniform on the rational support of the Pytlik functional; irrational parameters do not enter the rigorous trace argument. For general nilpotent models, coefficient-weighted spectral sums are justified to every fixed order by positive Rockland estimates, the Plancherel-Mellin formula, and a trace-level order-balance argument. In contrast with approaches which usually give leading terms or integrated Edgeworth-type asymptotics, the method gives genuinely local, pointwise higher-order heat-kernel expansions. The same representation-theoretic quantity governs the leading term in the closed-geodesic asymptotics, producing a nilpotent Chebotarev-type phenomenon. The Heisenberg model is computed to the first correction term, and the Engel model is represented through the resolvent and heat-kernel calculus of the quartic oscillato


4. Poincaré–Einstein 流形共形填充的尖锐相对比较不等式

A sharp relative comparison inequality for conformal fillings of Poincaré–Einstein manifolds

作者: Nan Wu

类别: math.DG, math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13742v1

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

摘要: Let $(X^{n+1},g_+)$ be a Poincaré–Einstein manifold with conformal infinity $(M^n,[h])$ of positive Yamabe type. We prove the sharp relative comparison inequality

$$ \frac{Y_1(X,M,[\bar g])}{Y_1(\mathbb{S}^{n+1}_+,\mathbb{S}^n,[g_{\mathbb{S}_+^{n+1}}])} \geq \left(\frac{Y(M,[h])}{Y(\mathbb{S}^n,[g_{\mathbb{S}^n}])}\right)^{\frac{n}{n+1}} $$

for the type-I Escobar–Yamabe compactification, and establish the rigidity. This confirms a conjecture proposed by Sun-Yung A. Chang.


5. 基于残差神经网络的含参问题解的逼近

Approximation of solutions of parameter-dependent problems by residual neural networks

作者: Ana Carpio

类别: math.NA, cs.LG, math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13574v1

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

摘要: We develop a convergent scheme to train neural networks involving analytic activation functions based on gradient flows. Convergence properties are guaranteed by Lojasiewicz theory. The main advantage of this approach is its simplicity of implementation. The coefficients of the network are approximated by solving a system of ordinary differential equations. We test the method by constructing residual neural network approximations of solutions of parametric problems. The dependence of the solutions of simple ordinary differential equations on a few parameters is correctly reproduced. The solutions of inverse problems involving wave constraints which depend on a few parameters can be reasonably approximated, even in regions in which the problem is severely ill posed.


6. 菱形环面上的几何分析:具有两个奇点的Green函数

Geometric analysis on rhombus torus: Green function with two singularities

作者: Zhijie Chen, Erjuan Fu, Chang-Shou Lin, Zhen Song

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13392v1

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

摘要: Let $G(z)$ be the Green function on the flat torus $E_τ=\mathbb{C}/(\mathbb{Z}+\mathbb{Z}τ)$ with the singularity at $0$. Lin and Wang (Ann. Math. 2010) proved that $G(z)$ has at most one pair of nontrivial critical points. This is the third of a series of papers to study the sum of two Green functions which can be reduced to $G_p(z):=\frac12(G(z+p)+G(z-p))$. We study how the geometry of the torus and the location of singularities $\pm p$ affect the structure of critical points of $G_p(z)$. In Part I \cite{CFL}, we proved that $G_p(z)$ has at most three pairs of nontrivial critical points for all tori. In Part II \cite{CFL-II} (Proc. Lond. Math. Soc. 2026), we studied the important case that $E_τ$ is a rectangular torus. In this paper, we study the other important but more challenging case that $E_τ$ is a rhombus torus, by developing different approaches from \cite{CFL, CFL-II}. As applications, we show that the curvature equation $Δu+e^{u}=4π(δ_p+δ_{-p})$ on $E_τ$ has exactly either $0$, $1$ or $2$ even axisymmetric solutions and each number really occurs.


7. 极小面积流形奇点分析,第三部分:平面频率 $\neq 2$ 的分支点、高阶渐近性及局部拓扑

Analysis of singularities of area-minimizing currents, Part III: branch points of planar frequency $\neq$ 2, higher order asymptotics, and the local topology

作者: Brian Krummel, Neshan Wickramasekera

类别: math.DG, math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13356v1

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

摘要: This is the third part in a series of papers developing a new framework to study the local structure of $n$-dimensional area-minimizing rectifiable currents $T$ of codimension $\geq 2$. Parts I and II introduced an intrinsic frequency function for $T$ – planar frequency – and used its monotonicity properties, among other things, to establish that ${\mathcal H}^{n-2}$-a.e. branch point is a rapid-decay branch point where the planar frequency is at least $1 + α$. This paper analyses branch points of planar frequency $\neq 2$. It establishes: (1) higher order asymptotics: at ${\mathcal H}^{n-2}$-a.e. such point, the current admits an expansion of finite order $>1$, with precise decay estimates for the remainder term; (2) branch set decomposition: the set of such branch points locally decomposes into finitely many pairwise disjoint, locally $n-2$ rectifiable sets (of locally finite measure); (3) topological control: near any branch point satisfying a specific planar-frequency criterion, the support of $T$ is homeomorphic to an $n$-dimensional disk and admits a $C^{1, μ}$ parametrization. (Classical complex algebraic examples show that when this frequency criterion fails, the current need not be locally homeomorphic to an $n$-disk). The work here (as well as in parts I & II) avoids the use of center manifolds – a technically demanding foundational component of the classical Almgren framework – and uses instead intrinsic geometric arguments based on the monotonicity formula for planar frequency. In part IV, a center manifold is utilised to analyse planar frequency 2 points, where the center manifold becomes necessary and geometrically canonical, satisfying additional simplifying properties. Reduced reliance on center manifolds in our framekwork is necessitated by the structural results it establishes for $T$.


8. 具有梯度超线性增长的完全非线性椭圆方程的Liouville型定理

Liouville type theorems for fully nonlinear elliptic equations with superlinear growth in gradient

作者: Amit Kumar Acharya, Ram Baran Verma

类别: math.AP

发布日期: 2026-07-15

arXiv ID: 2607.13355v1

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

摘要: This article investigates positive supersolutions of the fully nonlinear elliptic system

$$ \begin{cases} -\mathcal{M}_{λ,Λ}^{+}(D^{2}u)+|\nabla u|^{q} \geqλ_{1}f_{1}(v)~~\text{in}~~\mathbb{R}^{n}\setminus B_{R_0},\\ -\mathcal{M}_{λ,Λ}^{+}(D^{2}v)+|\nabla v|^{q} \geqλ_{2}f_{2}(u)~~\text{in}~~\mathbb{R}^{n}\setminus B_{R_0}, \end{cases} $$

where $q>1,λ_{1},λ_{2}>0,$ and the nonlinearities exhibit power-type behaviour either near the origin or at infinity. Introducing the effective dimension $\widetilde n_{+}=\fracλΛ(n-1)+1$ associated to extremal Pucci operator, we identify the critical exponent $q_{c}=\frac{\widetilde n_{+}}{\widetilde n_{+}-1},$ which governs the qualitative behaviour of positive supersolutions. Using this framework, we establish sharp Liouville-type nonexistence theorems in exterior domains and determine optimal nonexistence regions through the interaction between the gradient exponent $q,$ the effective dimension $\widetilde n_{+}$ and the nonlinear couplings. In the prototype case $f_{1}(t)=t^{p_{1}},$ $f_{2}(t)=t^{p_{2}}$ the obtained conditions are shown to be optimal. The analysis is carried out under the natural regularity assumption $u,v\in W^{2,p}_{\mathrm{loc}}(\mathbb{R}^{n}\setminus B_{R_0}),$ for $p>n,$ which is the regularity available for fully nonlinear uniformly elliptic equations, rather than within a classical $C^{2}$ framework. Our results provide the fully nonlinear Pucci analogue of the Liouville theory for semilinear elliptic systems involving nonlinear gradient terms and reveal the fundamental role of the effective dimension in determining the critical nonexistence thresholds.


9. 非局部梯度、非局部拉普拉斯算子与极值原理

Nonlocal gradient, the nonlocal Laplacian and maximum principles

作者: José Carlos Bellido, Guillermo García-Sáez, José Camilo Rueda-Niño

类别: math.AP

发布日期: 2026-07-14

arXiv ID: 2607.13161v1

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

摘要: We study the nonlocal $ρ$-Laplacian, defined as the composition of the nonlocal divergence and gradient operators associated with a general radial kernel $ρ$: $Δ_ρu=\mbox{div}_ρ\left(D_ρu\right)$. Our first main contribution is to establish a precise connection between this operator and the class of integro-differential elliptic operators studied by Fernández-Real and Ros-Oton (\cite{FernandezRos}), identifying explicit conditions on the kernel $ρ$ that guarantee membership in this class. Our second main contribution concerns maximum and comparison principles for the $ρ$-Laplacian. We establish both a strong and a weak maximum principle under conditions on $ρ$ that are strictly weaker than those required for membership in the integro-differential class, thereby covering a genuinely broader family of operators. The results require only minimal assumptions on the kernel, and in particular do not rely on any fractional-type comparability condition.