生成时间: 2026年07月28日 05:37:41

运行状态: 完成

计划处理论文数量: 51

已完成论文数量: 51

未完成论文数量: 0

重点关注论文数量: 12

了解领域论文数量: 19

不相关论文数量: 20

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

1. 三维拟耗散 Navier–Stokes 方程的正则性范数激增

Positive-regularity norm inflation for the 3D hypodissipative Navier–Stokes equations

作者: Guirong Tang, Shiyang Xiong

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24635v1

摘要翻译: 我们证明了具有耗散算子 $(-Δ)^α$(其中 $0<α<1$)的三维拟耗散 Navier–Stokes 方程的同空间范数激增现象。设 $2\le p<\infty$ 且

Comment: 无

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

1. 研究对象和背景

本文研究三维弱耗散 Navier-Stokes 方程(hypodissipative Navier–Stokes equations)

$$ \begin{cases} \partial_t u + u \cdot \nabla u + \nabla \pi + (-\Delta)^\alpha u = 0, \\ \operatorname{div} u = 0, \\ u|_{t=0} = u_0, \end{cases} \quad (t,x) \in [0,T] \times \mathbb{R}^3, $$

其中耗散阶数 $0 < \alpha < 1$,非局部算子 $(-\Delta)^\alpha$ 定义为傅里叶乘子 $|\xi|^{2\alpha}$。该模型介于经典 Navier-Stokes($\alpha = 1$)与 Euler($\alpha = 0$)之间,是研究耗散临界性、尺度不变性与解的适定性边界的重要桥梁。其自然尺度变换为

$$ u(t,x) \mapsto u_\lambda(t,x) := \lambda^{2\alpha - 1} u(\lambda^{2\alpha} t, \lambda x), $$

诱导齐次 Sobolev 范数的缩放律:

$$ \|u_{0,\lambda}\|_{\dot{W}^{s,p}} = \lambda^{s + 2\alpha - 1 - 3/p} \|u_0\|_{\dot{W}^{s,p}}. $$

由此导出尺度临界指标

$$ s_c(p,\alpha) := 1 - 2\alpha + \frac{3}{p}. $$

当 $s < s_c(p,\alpha)$ 时,空间 $W^{s,p}$(或 $B^s_{p,q}$)处于严格尺度超临界(strictly scaling-supercritical)区域;而 $s > 0$ 意味着仍属“正则性”范畴(positive-regularity),即初值具有可微性而非仅分布意义。这一区域既非能量类空间($s=1$),亦非端点临界空间(如 $\dot{B}^{1/3}_{3,\infty}$),其适定性长期悬而未决。已有工作表明:在负正则性空间(如 $\dot{B}^{-1}_{\infty,\infty}$)中,通过 Fourier 高频向低频级联可证范数膨胀(norm inflation),代表工作包括 Bourgain-Pavlović(2008)、Germain(2009)、Cheskidov-Dai(2014)等;而在正则性区域,Luo(2020, 2022)首次对经典 Navier-Stokes($\alpha = 1$)构造了基于各向异性涡环(anisotropic vortex ring)的物理空间混合机制:取细长环状支撑集,利用二维稳态 meridional 流驱动 swirl 在环面横截方向快速振荡,使 $\|u(t)\|_{W^{s,p}}$ 在短时内爆炸。该机制依赖几何误差(curvature error)主导、耗散项可忽略,但其稳定性在 $\alpha < 1$ 下完全未知——因 $(-\Delta)^\alpha$ 是非局部算子,其作用于紧支近似解 $\bar{u}$ 时产生全域尾部贡献,且当 $\alpha \to 0^+$ 时耗散减弱,传统能量估计失效。本文旨在填补这一关键空白:是否对全部 $0 < \alpha < 1$,正则性超临界区域 $0 < s < s_c(p,\alpha)$ 均存在范数膨胀?换言之,弱耗散是否足以抑制物理空间混合引发的强不适定性?该问题直接关联弱耗散方程的强解存在性阈值、临界正则性猜想的鲁棒性,以及非局部耗散对流体湍流生成机制的影响本质。

2. 主要定理或主要结果

定理 1.1(范数膨胀):设 $0 < \alpha < 1$,$2 \le p < \infty$,且

$$ 0 < s < s_c(p,\alpha) = 1 - 2\alpha + \frac{3}{p}. $$

令 $X = W^{s,p}$,或 $X = B^s_{p,q}$(其中 $1 \le q \le \infty$ 固定)。则对任意 $\varepsilon > 0$,存在无散度光滑紧支初值 $u_0 \in C_c^\infty(\mathbb{R}^3)$ 及时间 $0 < t_\varepsilon \le \varepsilon$,使得方程 (1.1) 的唯一光滑局部解 $u$ 在 $[0,t_\varepsilon]$ 上存在,并满足

$$ \|u_0\|_X \le \varepsilon, \quad \|u(t_\varepsilon)\|_X \ge \varepsilon^{-1}. $$

该结论为 same-space 膨胀:初值与终值范数在同一函数空间 $X$ 中测量,且 Besov 空间中求和指标 $q$ 保持固定。证明中关键参数由超临界间隙 $\delta := s_c(p,\alpha) - s > 0$ 控制:取 $b = \delta / 100$,$\nu = \mu^{1-b}$,$A = \zeta^2 \mu^{-s + 2/p} \nu^{1/p}$,$t_* = \zeta^{-N} A^{-1} \mu^{-1}$,其中 $\zeta \in (0,1]$ 固定,$\mu \gg 1$ 大参数。构造的近似解 $\bar{u}(t)$ 满足: (i)初值估计:$\|u_0\|_X \lesssim \zeta^2$; (ii)近似解增长:$\|\bar{u}(t_*)\|_{W^{s,p}} \gtrsim A K_*^s V_p = \zeta^{2-Ns}$,其中 $K_* = t_* A \mu^2 = \zeta^{-N} \mu$ 为相位梯度主导频率,$V_p = \mu^{-2/p} \nu^{-1/p}$ 为体积因子; (iii)残差控制:几何误差 $E_{\text{geom}}$ 满足 $\|\nabla^k E_{\text{geom}}(t)\|_{L^j} \lesssim_{\zeta,k,j} S_{k,j} A \nu$,耗散误差 $E_{\text{diss}} = (-\Delta)^\alpha \bar{u}$ 满足 $\|\nabla^k E_{\text{diss}}(t)\|_{L^j} \lesssim_{\zeta,k,j} S_{k,j} \mu^{2\alpha}$,且二者均被 $S_{k,j} A \nu$ 主导(因 $\mu^{2\alpha}/(A\nu) \to 0$); (iv)时间尺度:$t_* \to 0$,$t_* \mu^{2\alpha} \to 0$,确保耗散在膨胀时间内可忽略。特别地,当 $\delta = 0$(即 $s = s_c$)时,所有关键渐近关系失效,故结论不覆盖临界线。

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

本文核心障碍在于:非局部耗散 $(-\Delta)^\alpha$ 作用于紧支涡环近似解 $\bar{u}$ 时,其全域尾部 $\operatorname{supp}((-\Delta)^\alpha \bar{u}) = \mathbb{R}^3$ 不再紧支,导致标准局部化估计失效。作者的真正创新点在于建立了一种双重超临界间隙利用机制:严格超临界间隙 $\delta = s_c - s > 0$ 同时服务于两个目的——分离时间尺度与压制非局部误差。

第一,时间尺度分离。定义膨胀时间 $t_* = \zeta^{-N} A^{-1} \mu^{-1}$,其中 $A \sim \mu^{-s + 2/p} \nu^{1/p} = \mu^{-s + 2/p + (1-b)/p} = \mu^{-s + 3/p - b/p}$。由 $s < s_c = 1 - 2\alpha + 3/p$ 得 $-s + 3/p > 2\alpha - 1 + b/p$,故 $t_* \mu^{2\alpha} \sim \mu^{-\delta + b/p} \to 0$(引理 3.6)。这保证耗散项 $\int_0^{t_*} \|(-\Delta)^\alpha u\|_{X} dt$ 在膨胀时间内积分可忽略,从而 $\bar{u}$ 的动力学由 Euler 输运主导。

第二,非局部耗散误差的全局控制。关键引理 3.9 证明:对任意 $k,j$,

$$ \|\nabla^k (-\Delta)^\alpha \bar{u}(t)\|_{L^j} \lesssim S_{k,j} \mu^{2\alpha}, $$

其中 $S_{k,j} = A \mu^k V_j$ 为近似解的 $W^{k,j}$ 尺度。而几何误差满足 $\|\nabla^k E_{\text{geom}}\|_{L^j} \lesssim S_{k,j} A \nu$(引理 3.8)。比较两者得比值

$$ \frac{\mu^{2\alpha}}{A \nu} \sim \mu^{-\delta + b(1+1/p)} \to 0, $$

故耗散误差被几何误差所控,且二者均属高阶小量。此估计依赖于 $\Lambda^{2\alpha} = (-\Delta)^\alpha$ 作为傅里叶乘子的全局有界性(Riesz 变换在 $L^j$ 上有界),并通过插值 $\|\Lambda^{k+2\alpha} \bar{u}\|_{L^j} \lesssim \|\bar{u}\|_{L^j}^{1-(k+2\alpha)/m} \|\bar{u}\|_{W^{m,j}}^{(k+2\alpha)/m}$(取 $m > k+2\alpha$)实现,避免了局部化截断。

第三,相位振荡的精确下界。为证 $\|\bar{u}(t_*)\|_{W^{s,p}}$ 膨胀,需在环面区域 $R_* = \{\mu \rho \in I_g, \varphi \in J_\varphi\}$ 上控制 $\partial_z \Phi_* = \cos \varphi / \rho + t_* A / \rho^2 \sin \varphi$。当 $K_* = \zeta^{-N} \mu \gg 1$ 时,$|\partial_z \Phi_*| \sim K_*$(式 (4.4))。引理 4.2(加权非抵消引理)给出:对 $a \ge 0$、$\phi'$ 单调且 $\inf |\phi'| \ge \lambda$,有

$$ \int_I a(y) |\sin(\phi(y)+\beta)|^p dy \ge m_p \int_I a - C_{p,C_0} \lambda^{-1} (\|a'\|_{L^1} + |I|^{-1} \|a\|_{L^1}). $$

应用于 $a = g(\mu \rho)^p$、$\phi = \Phi_*$,利用 $\lambda \sim K_*$ 及 $a$ 的紧支性质,得到 $\|\bar{u}_\theta(t_*)\|_{L^p}^p \gtrsim A^p K_*^p \mu^{-2} \nu^{-1} = A^p K_*^p V_p^p$,再经整数阶导数增长与插值得到 $W^{s,p}$ 下界 $\gtrsim A K_*^s V_p$。

第四,稳定性分析。令 $w = u - \bar{u}$,则 $w$ 满足含残差 $E$ 的方程。作者采用 Lipschitz 引导的 $H^k$ 归纳法(而非标准能量法):由 $\|\nabla u\|_{L^\infty} \lesssim L = A \mu$(引理 3.7),结合 $\|E\|_{W^{k,j}} \lesssim S_{k,j} A \nu$,证明 $\|w(t)\|_{W^{k,j}} \lesssim t \cdot S_{k,j} A \nu \cdot L$,故在 $t_*$ 内 $\|w(t_*)\|_X \ll \|\bar{u}(t_*)\|_X$,确保真实解继承近似解的范数增长。

4. 与之前工作的比较

本文与既有工作存在三重本质区别:

(1)与 Fourier 级联型范数膨胀工作的比较:Bourgain-Pavlović(2008)、Germain(2009)及 Cheskidov-Dai(2014)均在负正则性空间(如 $\dot{B}^{-1}_{\infty,\infty}$ 或 $\dot{B}^{s}_{p,q}$ with $s<0$)中构造高频向低频的能量转移,其机制依赖于 Fourier 支集分离与共振条件,适用于广义耗散 $(-\Delta)^\alpha$ with $\alpha \ge 1$。但此类构造无法进入 $s>0$ 区域,因其要求初始数据在低频处几乎为零,而正则性空间要求高频衰减。本文的物理空间涡环机制则天然适配 $s>0$,且首次将该机制推广至 $0<\alpha<1$ 全范围,突破了 Fourier 方法的正则性天花板。

(2)与 Luo(2020, 2022)经典情形的比较:Luo 对 $\alpha=1$ 构造了 $W^{s,p}$ 和 $B^s_{p,q}$ 的范数膨胀,其关键在于经典 Laplacian 的局部性使几何误差 $E_{\text{geom}}$ 成为主导,而 $\Delta \bar{u}$ 的紧支性简化了残差估计。本文则必须处理 $(-\Delta)^\alpha \bar{u}$ 的非局部尾部,为此引入全新技术:利用超临界间隙 $\delta$ 同时控制时间尺度 $t_* \mu^{2\alpha} \to 0$ 和非局部误差比 $\mu^{2\alpha}/(A\nu) \to 0$。当 $\alpha \to 1^-$ 时,$\delta$ 的容许范围收缩,但构造仍连续过渡至 Luo 的结果(文中 Remark 1.2 明确指出);而 $\alpha \to 0^+$ 时,$s_c \to 1 + 3/p$,膨胀区域扩大,体现弱耗散对不适定性的‘放大’效应。

(3)与分数阶 Navier-Stokes 临界理论的比较:对于 $\alpha < 1$,已知 $\dot{H}^{\frac{5}{2} - 2\alpha}$ 是临界能量空间(Dong-Du, 2009),且 $\dot{B}^{\frac{5}{2} - 2\alpha}_{2,1}$ 中存在弱解非唯一性(Dai-Liu, 2021)。但这些工作聚焦于临界或亚临界空间,且解类为弱解或 mild 解。本文则严格限定于 $C_c^\infty$ 初值与唯一光滑解,测量同一正则性超临界空间中的范数爆炸,揭示了即使在最‘好’的解类下,方程仍强不适定。代价在于:结论仅适用于 $s < s_c$ 的开区间,无法触及 $s = s_c$(临界线),且 $p \ge 2$ 的限制源于插值中 Sobolev 嵌入 $\|f\|_{\dot{W}^{s,p}} \lesssim \|f\|_{\dot{H}^{s + d_p}}$($d_p = 3(1/2 - 1/p)$)的要求;对 $1 < p < 2$,文中未提供相应估计,从提供内容中无法确认是否可推广。


2. 爆破与集中但无量子化:尖锐的 Harnack 型不等式

Blow up and Concentration without Quantization: sharp Harnack type inequalities

作者: Daniele Bartolucci, Paolo Cosentino, Lina Wu

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24346v1

摘要翻译: 受 Onsager 统计力学对具有点奇性的湍流 Euler 流的描述所启发,我们改进了一类扰动奇异 Liouville 方程解序列的爆破分析,这类方程共享“爆破与集中但无量子化”现象。该问题极为精细,因为我们正处理恰好位于众所周知的“集中但无量子化”现象发生之临界阈值处的情形,这一方向最近在 [C.S. Lin, G. Tarantello, C. R. Math. Acad. Sci. Paris (2016)] 与 [Y. Lee, C.S. Lin, G. Tarantello, W. Yang, Comm. PDE. (2017)] 中被推进。首先,我们需要为此类特别丰富的奇异极限建立一种新的尖锐 Harnack 型不等式。然而这尚不足够,因为由奇性继承而来的共形因子的增长阻碍了经典量子化论证的应用。我们亦通过针对“快爆破”与“慢爆破”的不同策略解决此困难:精心适配基于 Pohozaev 恒等式、椭圆估计以及“Sup+CInf”不等式的论证,其精神与 [C.C. Chen, C.S. Lin, Comm. An. Geom. (1998)] 一致。

Comment: 39 pages. Comments are welcome!

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

1. 研究对象和背景

本文研究的是带奇异权的扰动Liouville方程序列的爆破分析,其原型为

$$ -\Delta v_n = H_n e^{v_n} \quad \text{in } B_2, $$

其中权函数 $H_n(x) = (\epsilon_n^2 + |x|^2)^{\alpha_{n,\sigma}} K_n(x)$,参数 $\alpha_{n,\sigma} = \frac{\lambda_n}{4\pi \sigma}$,且满足

$$ \lambda_n = \int_{B_2} (\epsilon_n^2 + |x|^2)^{\alpha_{n,\sigma}} K_n e^{v_n} \to \lambda_\infty > 0, \quad \epsilon_n \to 0^+, \quad K_n \in C^1(B_2), \; K_n \to K \text{ in } C^1_{\mathrm{loc}}(B_2), \; K \ge a > 0. $$

关键未知量是解序列 $\{v_n\}$ 的渐近行为,特别是当其在原点附近发生爆破(即 $\sup_{B_2} v_n = v_n(x_n) \to +\infty$, $x_n \to 0$)且满足边界振荡控制 $\max_{\partial B_1} v_n - \min_{\partial B_1} v_n \le C$ 时的质量集中性质。该模型源于Onsager统计力学对含点涡奇点的湍流Euler流的描述:$v_n$ 可解释为流函数,而 $\rho_n = \frac{H_n e^{v_n}}{\lambda_n}$ 是归一化涡度分布;原点处的奇异性对应一个固定的宏观反向旋转涡心,爆破则表征正负涡旋在该点的“碰撞”。经典理论中,当 $\alpha_{n,\sigma} \equiv 0$(正则Liouville)或 $\epsilon_n \equiv 0$(纯锥形奇点)时,爆破必然导致质量量子化——即极限测度 $\eta = \lim_{n\to\infty} H_n e^{v_n} \,dx$ 在爆破点 $0$ 处的原子质量 $\beta = \eta(\{0\})$ 必为 $8\pi \mathbb{N}$ 的整数倍(如 $8\pi, 16\pi, \dots$)。然而,在 $0 < \sigma < \frac{1}{2}$ 且 $\lambda_\infty \in (8\pi, \frac{4\pi}{\sigma})$ 时,[3] 已发现“无量子化的爆破与集中”现象:$\beta$ 不再被锁定为离散值,而可在非退化区间 $[8\pi, 8\pi(1+\alpha_{\infty,\sigma})]$ 内连续变化,即 $\beta \in [8\pi, \lambda_\infty]$,且 $\beta = \lambda_\infty$(完全集中)当且仅当边界振荡受控(Theorem 2.4)。本文聚焦于该现象的临界阈值——即 $\alpha_{\infty,\sigma} = 1$,亦即 $\lambda_\infty = \frac{4\pi}{\sigma}$,此时 $\sigma \in (0, \frac{1}{2}]$。此情形下,权函数渐近为 $|x|^{2\alpha_{\infty,\sigma}} = |x|^2$,恰为文献[18]所指出的“无集中爆破”(blow up without concentration)的分界线:当 $\alpha_{\infty,\sigma} > 1$ 时,爆破不再蕴含局部一致趋于 $-\infty$(即集中失效);而 $\alpha_{\infty,\sigma} = 1$ 则处于“集中但无量子化”的最精细边缘。现有工具在此阈值处全部失效:经典的“Sup+CInf”不等式(如[2])要求 $\alpha_{n,\sigma} < 1$ 以保证 $C_1 \inf_{B_1} v_n \le -v_n(x_n) + C_2$,但当 $\alpha_{n,\sigma} \equiv 1$ 时该不等式退化为平凡的 $\inf_{B_1} v_n \le C$,无法控制爆破强度;同时,Pohozaev恒等式中由奇异性贡献的余项 $\int_{\partial B_r} (x\cdot \nu) \nabla v_n \cdot \nabla v_n$ 在 $\alpha_{\infty,\sigma} = 1$ 下产生对数发散积分,使传统量化论证崩溃。因此,本文旨在填补这一关键空白:建立适用于 $\alpha_{\infty,\sigma} = 1$ 的尖锐Harnack型不等式,并据此完成“快/慢爆破”两类情形下的质量量化,从而给出 $\lambda_\infty = \frac{4\pi}{\sigma}$ 时爆破解的完整分类。

2. 主要定理或主要结果

在假设 $K_n \in C^1(B_2)$, $K_n \to K$ in $C^1_{\mathrm{loc}}(B_2)$, $\min_{B_1} K \ge a > 0$, $\|\nabla K_n\|_{L^\infty(B_1)} \le A$, 且解序列 $\{v_n\}$ 满足 $-\Delta v_n = H_n e^{v_n}$ in $B_2$, $\lambda_n = \int_{B_2} H_n e^{v_n} \to \lambda_\infty > 0$, $\sup_{B_2} v_n = v_n(x_n) \to +\infty$, $x_n \to 0$, $\max_{\partial B_1} v_n - \min_{\partial B_1} v_n \le C$, 以及 $\alpha_{n,\sigma} = \frac{\lambda_n}{4\pi \sigma} \le 1$, $\alpha_{\infty,\sigma} = \frac{\lambda_\infty}{4\pi \sigma} = 1$ 的条件下,本文证明了以下核心结论:(i) 尖锐Harnack不等式(Theorem 3.1):对三种爆破尺度,存在仅依赖于 $(a,A,\|K\|_{L^\infty})$ 的常数 $C_I,C_{II},C_{III} > 0$,使得

$$ v_n(x_n) + \inf_{B_1} v_n \le C_I - 4\log \epsilon_n, \quad \text{若 } \frac{\epsilon_n}{t_n} \to +\infty \text{(快爆破I)}; $$

$$ v_n(x_n) + \inf_{B_1} v_n \le C_{II} - 4\log |x_n|, \quad \text{若 } \frac{\epsilon_n}{t_n} \le C, \frac{|x_n|}{\delta_n} \to +\infty \text{(快爆破II)}; $$

$$ v_n(x_n) + \inf_{B_1} v_n \le C_{III} - 4\log \delta_n, \quad \text{若 } \frac{\epsilon_n}{t_n} \le C, \frac{|x_n|}{\delta_n} \le C \text{(慢爆破)}, $$

其中 $t_n = \max\{\delta_n, |x_n|\}$, $\delta_n^{2(1+\alpha_{n,\sigma})} = e^{-v_n(x_n)}$;(ii) 完整分类定理(Theorem 1.3):在上述假设下,必有 $\sigma \in (0,\frac{1}{2}]$ 且 $\lambda_\infty = \frac{4\pi}{\sigma}$,且存在子列使得仅以下三种情形之一成立:

$$ \sigma = \frac{1}{4}, \; \lambda_\infty = 16\pi; \quad \text{或} \quad \sigma = \frac{1}{2}, \; \lambda_\infty = 8\pi, \; \frac{\epsilon_n}{\delta_n} \to +\infty, \; \frac{|x_n|}{\delta_n} \to +\infty; \quad \text{或} \quad \sigma \in \left(\frac{1}{4}, \frac{1}{2}\right), \; \lambda_\infty \in (8\pi, 16\pi), \; \frac{\epsilon_n}{\delta_n} \to \epsilon_0 > 0, \; \frac{|x_n|}{\delta_n} \to 0. $$

此外,在第三种情形(慢爆破)下,解具有渐近剖面 $v_n(x) = v_n(x_n) + \widetilde{U}_n(\delta_n^{-1}x) + O(1)$,其中 $\widetilde{U}_n(y) \to \widetilde{w}(y)$ 且 $\widetilde{w}$ 是 $-\Delta \widetilde{w} = c_0(\epsilon_0^2 + |y|^2) \widetilde{w}$ 在 $\mathbb{R}^2$ 上的解,满足 $\int_{\mathbb{R}^2} c_0(\epsilon_0^2 + |y|^2) e^{\widetilde{w}} = \lambda_\infty$;(iii) 量化结果:在快爆破I情形($\sigma = \frac{1}{2}, \lambda_\infty = 8\pi$)下,证明了 $\int_{B_r} H_n e^{v_n} = 8\pi + o(1)$ 且 $\int_{B_2\setminus B_r} H_n e^{v_n} = o(1)$;在快爆破II情形($\sigma = \frac{1}{4}, \lambda_\infty = 16\pi$)下,通过椭圆估计引理(Lemma 4.2)排除了质量区间中的非极值点,最终确认 $\beta = \lambda_\infty = 16\pi$。

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

本文的核心技术障碍在于 $\alpha_{\infty,\sigma} = 1$ 导致的经典工具全面失灵:一方面,“Sup+CInf”不等式失效,因其依赖 $\alpha_{n,\sigma} < 1$ 以获得对 $\inf v_n$ 的线性控制;另一方面,Pohozaev恒等式中由奇异性权重 $|x|^{2\alpha_{n,\sigma}}$ 贡献的边界项在 $\alpha_{n,\sigma} \to 1$ 时产生 $\int_{\partial B_r} |\nabla v_n|^2$ 级别的对数发散,使质量量化无法进行。作者的真正创新点在于针对不同爆破尺度设计三套互补策略:(1) 对“快爆破”(Cases I & II),构造基于移动平面法的尖锐Harnack不等式(Theorem 3.1)。具体而言,在极坐标下定义变换 $u_n(t,\phi) = v_n(e^t \cos\phi, e^t \sin\phi) + 2t - \frac{A}{a} e^t$,利用 $K_n$ 的梯度有界性($|\partial_t \widehat{K}_n| \le A e^t$)证明 $\partial_t(W_n e^\xi + \frac{A}{a} e^t) > 0$,从而应用Alexandrov-Bakelman-Pucci型比较原理,导出 $u_n$ 在半无限条带 $Q = \{t \le 0\}$ 上的对称性控制,最终得到 $v_n(x_n) + \inf_{B_1} v_n$ 的精确上界,其主导项 $-4\log \epsilon_n$ 或 $-4\log |x_n|$ 恰好抵消权重增长;(2) 对快爆破I($\frac{\epsilon_n}{t_n} \to +\infty$),将Pohozaev恒等式应用于缩放解 $w_n(z) = v_n(\epsilon_n z) + 2(1+\alpha_{n,\sigma})\log \epsilon_n$,其极限方程为 $-\Delta w = c_0(1+|z|^2) e^w$,此时Pohozaev余项 $\int_{\partial B_R} (z\cdot \nu)(\partial_\nu w)^2$ 在 $R\to\infty$ 下收敛,结合 $\int_{\mathbb{R}^2} c_0(1+|z|^2) e^w = \lambda_\infty$ 的约束,唯一可能解为 $\lambda_\infty = 8\pi$ 且 $\sigma = \frac{1}{2}$(Theorem 3.5);(3) 对快爆破II($\frac{\epsilon_n}{t_n} \le C, \frac{|x_n|}{\delta_n} \to +\infty$),Pohozaev分析给出质量 $\beta$ 仍可取区间值,此时作者引入精细椭圆估计:利用Green函数表示 $v_n$ 并分离主部与余部,证明余部在 $B_2\setminus B_{|x_n|/2}$ 上的一致有界性(Lemma 4.2),从而将不确定形式 $\int_{B_2\setminus B_r} H_n e^{v_n}$ 控制为 $o(1)$,最终由 $\beta \ge 8\pi$ 及 $\lambda_\infty = 16\pi$ 推出 $\beta = 16\pi$(Proposition 4.1);(4) 对“慢爆破”(Case III),Harnack不等式(3.5)退化为 $\inf v_n \le C$,故转而利用该情形下极限质量 $\lambda_\infty > 8\pi$ 的事实,适配[11]的论证:在缩放域 $D_n = B_{r/\delta_n}(0)$ 上,对 $\widetilde{w}_n(y) = v_n(\delta_n y) + 2(1+\alpha_{n,\sigma})\log \delta_n$ 应用最大值原理与梯度估计,证明存在 $C>0$ 使得 $\sup_{D_n} \widetilde{w}_n - C \inf_{D_n} \widetilde{w}_n \le C$(Lemma 5.1),进而沿用[20]的标准论证得出外部质量消失(Lemma 5.3)。

4. 与之前工作的比较

本文与既有工作的关键推进体现在三个维度:(i) 阈值突破:此前[3]处理的是 $\alpha_{\infty,\sigma} < 1$(即 $\lambda_\infty < \frac{4\pi}{\sigma}$)的严格亚临界情形,其“Sup+CInf”不等式([2])和Pohozaev量化均依赖 $\alpha_{\infty,\sigma} < 1$ 以保证衰减;而本文首次攻克 $\alpha_{\infty,\sigma} = 1$ 这一精确临界值,揭示了该阈值处爆破行为的相变本质——当 $\alpha_{\infty,\sigma} < 1$ 时集中必然发生,而 $\alpha_{\infty,\sigma} = 1$ 时集中虽存但量化结构剧变。(ii) 方法革新:对比[20](经典质量量化)依赖 $\alpha_{n,\sigma} \equiv 0$ 下的指数衰减,本文在 $\alpha_{\infty,\sigma} = 1$ 下发展出基于移动平面的尖锐Harnack不等式,其衰减率 $-4\log \epsilon_n$ 精确匹配权重 $|x|^2$ 的增长,这是此前所有工作(如[1],[5],[19])均未触及的尺度;对比[18]对 $\alpha_{\infty,\sigma} > 1$ 的“无集中”分析,本文证明 $\alpha_{\infty,\sigma} = 1$ 仍保持集中,但需全新工具规避对数发散。(iii) 分类完备性:[3]的Theorem 2.5仅覆盖 $\lambda_\infty < \frac{4\pi}{\sigma}$,本文Theorem 1.3将其延拓至等号情形,并给出 $\sigma$ 与 $\lambda_\infty$ 的精确耦合关系:$\sigma = \frac{1}{4} \iff \lambda_\infty = 16\pi$,$\sigma = \frac{1}{2} \iff \lambda_\infty = 8\pi$,这与[16]的径向解存在性条件完全吻合。代价在于:对 $\sigma = \frac{1}{4}$ 和 $\sigma = \frac{1}{2}$ 的两种快爆破情形,本文未能给出完整剖面(Remark after Theorem 1.3),因缩放后权重 $|y|^2$ 导致积分 $\int_{\mathbb{R}^2\setminus B_R} |y|^2 e^{w}$ 对数发散,暗示可能存在非标准多重爆破(non-standard multiple bubbling),此问题留待后续研究。此外,本文假设 $K_n \to K$ in $C^1_{\mathrm{loc}}$ 且 $K>0$,而[7]等处理更一般的 $K$ 变号情形,故本文适用范围更窄但精度更高。


3. 液态Lane-Emden恒星的转折点原理

A turning point principle for liquid Lane-Emden stars

作者: King Ming Lam

类别: math-ph, math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24530v1

摘要翻译: 在固定绝热指数的情形下,由Euler-Poisson系统与‘刚化气体’状态方程 $p=ρ^γ-1$ 所支配的球对称液态Lane-Emden恒星构成一个单参数族,自然地以中心密度 $κ=\barρ_κ(0)\in(1,\infty)$ 为参数。与气态情形不同,液态自由边界破坏了该族的自相似性,并使质量-半径曲线发生弯曲;当 $γ<2(d-1)/d$ 时,该弯曲导致质量出现极值点。我们仿照Zel’dovich、Wheeler与Thorne的思想,以及Hadžić、Lin与Rein关于相对论情形的严格对应工作,提出一个转折点原理:径向增长模的个数等于线性化算子的负Morse指标,该数沿参数族局部为常数,且仅能在质量-半径曲线的转折点处发生变化;该处弯曲的方向决定了是获得还是失去一个增长模。我们刻画了大中心密度极限情形,在此极限下,气态尾部(由一个平面动力系统读出)决定了质量-半径曲线是否呈螺旋状,从而决定了增长模个数是有界还是趋于无穷。作为推论,我们重新获得并加强了已知的液态Lane-Emden恒星的径向(不)稳定性结果;并结合既有的非径向分析,得到了一个非径向稳定性的转折点判据。

Comment: 32 pages

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

1. 研究对象和背景

本文研究的是三维及更高维空间中、由Euler-Poisson系统描述的自引力可压缩流体在液态物态方程下的球对称稳态解——液态Lane-Emden恒星。其控制方程为自由边界Euler-Poisson系统:

$$ \partial_t \rho + \nabla \cdot (\rho u) = 0 \quad \text{in } \Omega(t), $$

$$ \rho(\partial_t + u \cdot \nabla)u + \nabla p + \rho \nabla \phi = 0 \quad \text{in } \Omega(t), $$

$$ \Delta \phi = 4\pi \rho \quad \text{in } \mathbb{R}^d, \quad \lim_{|x|\to\infty} \phi(t,x) = 0, $$

其中 $d \geq 3$,$\Omega(t) := \{ x : \rho(t,x) > 0 \}$ 为星体支撑集,压力由液态(stiffened gas)物态方程封闭:

$$ p = \rho^\gamma - 1, \quad \gamma \in [1,2). $$

该方程在物理上对应于近不可压缩流体模型;当 $\gamma = 1$ 时退化为 $p = \rho - 1$,即Christodoulou两相模型中硬相的物态方程,亦与Zel’dovich在极高密度下提出的刚性物态方程一致。真空自由边界 $\partial \Omega(t)$ 上施加 $p = 0$ 及法向速度匹配条件 $u \cdot n = V_{\partial \Omega}$。特别地,常数项 $-1$ 的归一化使得密度在液面处恰好跳变为 $\bar{\rho}_\kappa(R_\kappa) = 1$,从而形成非平凡的密度间断——这是液态模型区别于气态模型的根本特征。

时间无关、球对称稳态解 $\bar{\rho}_\kappa(r)$ 满足静力学平衡方程(径向Ode):

$$ 0 = \begin{cases} \displaystyle \frac{\gamma}{\gamma - 1} \Delta \bar{\rho}_\kappa^{\gamma - 1} + 4\pi \bar{\rho}_\kappa, & \gamma > 1, \\ \Delta \ln \bar{\rho}_\kappa + 4\pi \bar{\rho}_\kappa, & \gamma = 1, \end{cases} \quad \nabla \bar{\rho}_\kappa(0) = 0, $$

其中拉普拉斯算子为径向形式 $\Delta = r^{1-d} \partial_r (r^{d-1} \partial_r)$。对每个中心密度 $\kappa := \bar{\rho}_\kappa(0) \in (1,\infty)$,该ODE存在唯一递减解 $f_{0,\kappa}$,液态恒星定义为截断解:

$$ \bar{\rho}_\kappa(r) := f_{0,\kappa}(r) \mathbf{1}_{[0,R_\kappa]}(r), \quad \text{where } f_{0,\kappa}(R_\kappa) = 1. $$

质量与半径分别为:

$$ M_\kappa := 4\pi \int_0^{R_\kappa} y^{d-1} \bar{\rho}_\kappa(y) \, dy, \quad R_\kappa > 0 \text{ 满足 } \bar{\rho}_\kappa(R_\kappa) = 1. $$

由此构成单参数族 $\kappa \mapsto (\bar{\rho}_\kappa, R_\kappa)$,其在 $(R,M)$ 平面上描绘出质量-半径曲线 $C = \{ (R_\kappa, M_\kappa) : \kappa \in (1,\infty) \}$。与气态Lane-Emden星不同,液态模型因密度跳变破坏了尺度不变性:气态族满足精确自相似性 $\bar{\rho}_\kappa(r) = \kappa \bar{\rho}_*(\kappa^{(2-\gamma)/2} r)$,导致 $M_\kappa \propto R_\kappa^{d\gamma - 2(d-1)}$ 为单一幂律,无极值;而液态截断引入边界效应,使 $C$ 弯曲并产生质量极值,当且仅当 $\gamma < \gamma_* := 2(d-1)/d$。该阈值正是经典气态稳定性临界值(Lin–Jang–Makino),但在液态情形下成为转折点出现的必要条件。

线性稳定性分析聚焦于径向扰动模态 $\zeta(y,t) = e^{\lambda t} \chi(y)$,在共动(Lagrangian)坐标系下导出奇异Sturm-Liouville算子:

$$ L_\kappa \chi = -\lambda^2 \chi \quad \text{on } [0,R_\kappa], $$

$$ d\chi(R_\kappa) + R_\kappa \partial_y \chi(R_\kappa) = 0, $$

其中

$$ L_\kappa \chi := -\frac{\gamma}{y^{d+1} \bar{\rho}_\kappa} \partial_y \left( \bar{\rho}_\kappa^\gamma y^{d+1} \partial_y \chi \right) + d(\gamma_* - \gamma) \frac{\partial_y \bar{\rho}_\kappa^\gamma}{y \bar{\rho}_\kappa} \chi. $$

该算子定义在加权Sobolev空间 $H_\kappa = H^1([0,R_\kappa], y^{d+1} \bar{\rho}_\kappa)$ 上,其能量型为:

$$ \langle L_\kappa \chi_1, \chi_2 \rangle_{L_\kappa} = \int_0^{R_\kappa} \left[ \gamma \bar{\rho}_\kappa^\gamma y^{d+1} (\partial_y \chi_1)(\partial_y \chi_2) + d(\gamma_* - \gamma) y^d \chi_1 \chi_2 \partial_y \bar{\rho}_\kappa^\gamma \right] dy + d\gamma R_\kappa^d \chi_1(R_\kappa) \chi_2(R_\kappa). $$

最后一项 $d\gamma R_\kappa^d \chi_1(R_\kappa) \chi_2(R_\kappa)$ 是液态边界贡献,不存在于气态情形,正是它导致质量极值与谱变化。增长模个数 $n_u(\kappa)$ 定义为负特征值个数(计重数),等价于该能量型的负Morse指标 $n_-(L_\kappa)$。问题核心在于:如何将动力学不稳定性(即 $L_\kappa$ 负谱的存在性与个数)完全编码于可观测的几何量——质量-半径曲线 $C$ 的弯曲行为中?已有工作如Lin–Zeng(气态)、Hadžić–Lin–Rein(相对论)已建立类似原理,但均依赖约束变分框架与绕数指标修正;而液态模型因共动公式自动守恒总质量,有望实现更简洁的指标等价。本文旨在填补这一Newtonian自由边界领域的严格空白,并揭示其与相对论螺旋现象的深刻类比。

2. 主要定理或主要结果

定理2.7(Morse指标计数与边缘模):对任意 $d \geq 3$、$\gamma \in [1,2)$ 及 $\kappa \in (1,\infty)$,线性化径向算子 $L_\kappa$ 在加权空间 $H_\kappa$ 上具有实谱,特征值有下界且仅在 $+\infty$ 处聚点,故增长模个数有限且精确等于负Morse指标:

$$ n_u(\kappa) = n_-(L_\kappa) < \infty. $$

边缘模 $\nu_\kappa$(即沿稳态族的形变方向)属于 $H_\kappa$,并满足分布恒等式:

$$ \langle L_\kappa \nu_\kappa, \chi \rangle = c_\kappa (\partial_\kappa M_\kappa) \chi(R_\kappa), \quad c_\kappa = -\frac{M_\kappa}{4\pi R_\kappa^{d-2}} \neq 0, $$

等价于 $L_\kappa \nu_\kappa = c_\kappa \partial_\kappa M_\kappa R_\kappa^{-d-1} \delta_{R_\kappa}$(Dirac测度)。进而有

$$ \dim \ker L_\kappa = \begin{cases} 1, & \text{若 } \partial_\kappa M_\kappa = 0, \\ 0, & \text{若 } \partial_\kappa M_\kappa \neq 0. \end{cases} $$

定理2.8(转折点原理):

(i) (静态判据)$n_u(\kappa)$ 在 $M_\kappa$ 的任意非临界区间上为常数;在临界点 $\kappa_0$(即 $\partial_\kappa M_\kappa|_{\kappa_0} = 0$ 且 $\partial_\kappa^2 M_\kappa|_{\kappa_0} \neq 0$)处,当 $\kappa$ 从小于 $\kappa_0$ 增至大于 $\kappa_0$ 时,$n_u$ 的跳跃量为:

$$ \Delta n_u = \begin{cases} +1, & \text{若 } (\partial_\kappa M_\kappa)(\partial_\kappa R_\kappa) \text{ 由负变正}, \\ -1, & \text{若 } (\partial_\kappa M_\kappa)(\partial_\kappa R_\kappa) \text{ 由正变负}. \end{cases} $$

等价于 $\Delta n_u = \Delta i_\kappa$,其中 $i_\kappa$ 为定义2.6的转折指数。

(ii) (阈值二分与失稳起始):

  • 若 $\gamma \geq \gamma_* = 2(d-1)/d$,则 $\kappa \mapsto M_\kappa$ 严格单调,$C$ 无转折点,且 $n_u(\kappa) \equiv 0$(全局线性稳定);
  • 若 $\gamma < \gamma_*$ 且 $\gamma \geq \gamma_\sharp = 2d/(d+2)$ 或判别式 $D(\gamma,d) < 0$(对所有 $d < 10$ 自动成立),则存在首个严格局部质量极大值点 $\kappa_1$,满足 $\partial_\kappa R_\kappa|_{\kappa_1} < 0$,且
    $$ n_u(\kappa) = \begin{cases} 0, & 1 < \kappa \leq \kappa_1, \\ 1, & \kappa_1 < \kappa < \kappa_1', \end{cases} $$
    其中 $\kappa_1'$ 为下一个质量临界点;此外对所有充分大的 $\kappa$,$n_u(\kappa) \geq 1$。

(iii) (大中心密度渐近与质量-半径螺旋):设 $\gamma < \gamma_*$,

  • (a) 若 $\gamma < \gamma_\sharp$ 且 $D(\gamma,d) < 0$($d < 10$ 时恒成立),则气态尾部动力系统平衡点 $v_*$ 为稳定焦点,$C$ 当 $\kappa \to \infty$ 时螺旋收敛至奇异液态星点 $(R_\infty, M_\infty)$,具有无穷多个转折点,且
    $$ \lim_{\kappa \to \infty} n_u(\kappa) = \infty; $$
  • (b) 若 $\gamma < \gamma_\sharp$ 且 $D(\gamma,d) \geq 0$(仅当 $d \geq 10$ 可能),则 $v_*$ 为稳定结点,$C$ 以渐近切向收敛,仅有有限个转折点,$n_u(\kappa)$ 收敛至某有限值(该值未确定);
  • (c) 若 $\gamma_\sharp \leq \gamma < \gamma_*$,则 $C$ 为单弧线,从原点出发并返回原点,具有限个转折点,$n_u(\kappa)$ 收敛至有限奇数值;特别地,在 $\gamma = \gamma_\sharp$ 时族显式可解,质量恰有一个临界点(严格最大值)位于 $\kappa_1 = \left( \frac{2(d-1)}{d-2} \right)^{(d+2)/2}$,且稳定化值为 $1$。判别式为
    $$ D(\gamma,d) := (d-2)^2 (2-\gamma)^2 - 8(2(d-1) - d\gamma). $$

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

证明的核心障碍在于:液态自由边界导致线性化算子 $L_\kappa$ 具有非标准Robin边界条件与显著的表面项,使得传统气态或相对论框架中的约束变分方法(需强制总质量守恒并引入绕数指标)不再适用;同时,$L_\kappa$ 的系数在 $y=0$ 处退化($y^{d+1}$ 权重),在 $y=R_\kappa$ 处存在跳跃,要求精细的加权Sobolev理论与弱形式处理。作者的真正创新点在于构建并运用一个无约束的共动变分框架,其关键工具与技术如下:

  1. 共动(Lagrangian)质量守恒嵌入:通过选择标记映射 $\eta(y,t) = y(1+\zeta(y,t))$ 并设定参考构型满足 $f J = \bar{\rho}_\kappa$,使得连续性方程自动蕴含总质量守恒 $\int_0^{R_\kappa} y^{d-1} \bar{\rho}_\kappa(y) \, dy = \text{const}$。这直接消除了气态/相对论分析中必需的“零均值约束”,从而避免了将 $L_\kappa$ 约化到子空间 $\{ \chi \in H_\kappa : \int y^{d-1} \bar{\rho}_\kappa \chi \, dy = 0 \}$ 的步骤,也无需引入绕数指标来修正Morse指标——$n_u(\kappa)$ 直接等于完整空间上的 $n_-(L_\kappa)$。

  2. 边缘模的精确构造与核刻画:利用稳态族对参数 $\kappa$ 的光滑依赖性,定义边缘模 $\nu_\kappa$ 为满足 $\partial_\kappa \bar{\rho}_\kappa + y^{1-d} \partial_y (y^d \bar{\rho}_\kappa \nu_\kappa) = 0$ 的解。通过直接计算其与 $L_\kappa$ 的配对,证明 $L_\kappa \nu_\kappa$ 是集中在 $y=R_\kappa$ 的Dirac测度,比例系数为 $c_\kappa \partial_\kappa M_\kappa$。此恒等式是转折点原理的枢纽:它表明 $\ker L_\kappa \neq \{0\}$ 当且仅当 $\partial_\kappa M_\kappa = 0$,且此时 $\nu_\kappa$ 即为零模;进而通过摄动理论证明 $\partial_\kappa M_\kappa$ 的符号变化率决定 $L_\kappa$ 负谱的增减。

  3. 气态尾部的平面动力系统提取与判别式分析:针对大 $\kappa$ 极限,作者将液态族的渐近行为映射至气态Lane-Emden方程的“尾部”——即当 $\kappa \to \infty$ 时,$\bar{\rho}_\kappa$ 在远离边界的区域趋近于缩放后的气态解 $\bar{\rho}_*$。通过变量替换 $s = \kappa^{(2-\gamma)/2} y$ 和渐近展开,导出支配尾部的二维自治系统(Proposition 3.7): $$ \frac{d}{ds} \begin{pmatrix} v_1 \ v_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \ -\frac{d}{d-2}(2-\gamma) & -\frac{d-2}{2} \end{pmatrix} \begin{pmatrix} v_1 \ v_2 \end{pmatrix}

  • \text{高阶项}, $$ 其线性化矩阵的特征方程判别式恰为 $D(\gamma,d)$。当 $D < 0$ 时,特征根为共轭复数,平衡点 $v_*$ 为稳定焦点,导致 $C$ 螺旋;当 $D \geq 0$ 时,为稳定结点,$C$ 单调收敛。此分析将几何曲线形态 $C$ 的全局行为(螺旋/节点)与谱渐近 $n_u(\kappa) \to \infty$ 或收敛完全耦合。
  1. 加权Sobolev嵌入与紧性论证:为确立 $L_\kappa$ 的离散谱,作者利用引理2.2给出的剖面估计 $1 \leq \bar{\rho}_\kappa(y) \leq \kappa$ 及 $|\partial_y \bar{\rho}_\kappa^\gamma| / y \sim \bar{\rho}_\kappa^2$,证明权重 $y^{d+1} \bar{\rho}_\kappa$ 满足Muckenhoupt $A_2$ 条件,从而 $H_\kappa \hookrightarrow L_\kappa$ 紧嵌入;结合 $L_\kappa$ 的强制性(coercivity)估计(由 $\gamma_* - \gamma$ 符号及边界项控制),最终确认其为下半有界自伴算子,谱离散。

4. 与之前工作的比较

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

  1. 模型范畴的实质性拓展:此前严格的转折点原理仅适用于两类模型:(i) Lin–Zeng (2017) 对气态Euler-Poisson系统的证明,其依赖自相似性与无自由边界,$M_\kappa$ 为幂律,无真实转折点;(ii) Hadžić–Lin–Rein (2021) 对Einstein-Euler系统的证明,其框架基于分离哈密顿结构与约束变分,需引入绕数指标 $w(\kappa)$ 并令 $n_u(\kappa) = n_-(\mathcal{L}_\kappa) - w(\kappa)$。本文首次为Newtonian自由边界液态模型建立严格原理,其物理设定(密度跳变、刚性物态)更贴近白矮星外层或强耦合等离子体,且因共动公式自动守恒质量,彻底规避了约束与绕数——这是方法论上的本质简化,而非技术微调。

  2. 假设条件的显著放宽与新现象揭示:气态理论要求 $\gamma > 1$ 且 $\gamma \neq \gamma_*$,而本文涵盖 $\gamma = 1$(线性物态 $p = \rho - 1$),并首次系统刻画 $\gamma < \gamma_\sharp = 2d/(d+2)$(无限支撑气态剖面)与 $\gamma \geq \gamma_\sharp$(紧支撑)的全域行为。尤其重要的是,本文发现维度 $d = 10$ 是判别式 $D(\gamma,d)$ 符号的临界阈值:当 $d < 10$,$D < 0$ 对所有 $\gamma \in [1,\gamma_\sharp)$ 成立,故物理维度($d = 3$)下液态星必螺旋且 $n_u \to \infty$;而 $d \geq 10$ 时存在“节点情形”,$n_u$ 收敛但极限值未知——此现象在气态或相对论工作中均未被识别,揭示了高维非线性椭圆理论中 $d = 10$ 的超临界角色(类比Lane-Emden方程径向解振荡性的经典结果)。

  3. 适用范围与代价的精确量化:本文结论覆盖整个 $\gamma \in [1,2)$,但存在明确局限:(i) 仅处理径向模态,非径向稳定性需依赖[37]的独立分析(Corollary 2.10将其与 $C$ 结合);(ii) “节点情形”($d \geq 10$, $\gamma < \gamma_\sharp$, $D \geq 0$)下 $n_u$ 的稳定化值未确定(Remark 2.8.3),此为当前技术瓶颈;(iii) 非线性稳定性未涉及,仅提供线性输入(如Hao–Miao (2023) 所用)。相较而言,相对论工作[21]虽也得出 $n_u \to \infty$,但其螺旋源于TOV方程的相对论效应,而本文螺旋纯属Newtonian液态边界与气态尾部动力学的耦合产物,且首次给出判别式 $D(\gamma,d)$ 的显式表达与维度依赖。此外,本文明确指出:对于Vlasov-Poisson系统(如球状星系),即使 $C$ 螺旋,也可能无动力学不稳定性(Ramming–Rein, 2015),故转折点原理绝非几何必然,而是模型特有性质——本文为此Newtonian自由边界族提供了首个完整证明,补全了该原理的模型谱系。


4. Lorentz 数据下无粘 SQG 方程的弱解

Weak Solutions for Inviscid SQG with Lorentz Data

作者: Peter Constantin, Mihaela Ignatova, Quoc-Hung Nguyen

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24500v1

摘要翻译: 我们为二维欧氏空间 $\mathbb R^2$ 及光滑有界区域上的无粘表面拟地转(SQG)方程构造了全局弱解,初始数据取自临界 Lorentz 空间 $L^{4/3,2}$ 中的任意函数。这些解在所有时刻守恒哈密顿量 $\|Λ^{-1/2}θ(t)\|_{L^2}^2$。第二个 Lorentz 指标由算子 $Λ^{-1/2}:L^{4/3,2}\to L^2$ 的尖锐有界性所确定;当 $q>2$ 时,该估计对 $L^{4/3,q}$ 不再成立。我们采用一种专为 Lorentz 空间设计的逼近方案:该方案对输运速度场与初始数据进行光滑化,并保持分布函数的序结构。通过高振幅截断解的 Lorentz 范数的一致估计,我们得以移除该逼近。

Comment: 无

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

1. 研究对象和背景

本文研究二维欧几里得空间 $\mathbb{R}^2$ 或光滑有界区域 $\Omega \subset \mathbb{R}^2$ 上的无粘表面准地转(Surface Quasi-Geostrophic, SQG)方程:

$$ \partial_t \theta + u \cdot \nabla \theta = 0, \quad u = \nabla^\perp \Lambda^{-1} \theta, \quad \Lambda = (-\Delta)^{1/2}. $$

该方程是地球物理流体中锋面形成的重要模型,亦被广泛视为三维不可压欧拉方程在二维的类比:其速度场 $u$ 由标量 $\theta$ 通过零阶奇异积分(即 $\nabla^\perp \Lambda^{-1}$)生成,而 $\nabla^\perp \theta$ 在结构上类比于三维欧拉方程的涡度 $\omega = \nabla \times u$。关键数学特征在于:(i) 传输项 $u \cdot \nabla \theta$ 是二次非线性且非局部;(ii) 速度算子 $\nabla^\perp \Lambda^{-1}$ 的核为 $K_{\text{SQG}}(x-y) = c (x-y)^\perp / |x-y|^3$,在对角线附近具有 $|x-y|^{-2}$ 奇异性;(iii) 方程保持哈密顿能量 $H(\theta) = \frac{1}{2} \| \Lambda^{-1/2} \theta \|_{L^2}^2 = \frac{1}{2} \int_\Omega \theta \, \Lambda^{-1} \theta \, dx$,该能量等价于负阶索博列夫范数 $\|\theta\|_{\dot{H}^{-1/2}}^2$。

从函数空间角度看,该能量自然要求 $\theta \in \dot{H}^{-1/2}$,而由二维 Riesz 位势理论,$\Lambda^{-1/2}: L^{p,q} \to L^2$ 的有界性临界情形发生在 $p = 4/3$。事实上,核 $|x|^{-3/2} \in L^{4/3,\infty}(\mathbb{R}^2)$,结合 Lorentz 卷积不等式 $\|f * g\|_{L^{r,c}} \lesssim \|f\|_{L^{p,a}} \|g\|_{L^{q,b}}$(其中 $r^{-1} = p^{-1} + q^{-1} - 1$, $c^{-1} = a^{-1} + b^{-1}$),可得 $\Lambda^{-1/2}: L^{4/3,2} \to L^2$ 连续,且该映射在次指数 $q = 2$ 处是尖锐的:对任意 $q > 2$,估计 $\|\Lambda^{-1/2} f\|_{L^2} \lesssim \|f\|_{L^{4/3,q}}$ 一般不成立(见附录 A)。因此,$L^{4/3,2}$ 是在主指数 $4/3$ 固定下、能控制哈密顿能量的最大 Lorentz 空间,即所谓的“临界 Lorentz 端点”。

已有工作集中于更规则的初值空间。例如:在 $L^2$ 中构造全局弱解([20]);推广至 $L^p$ 框架([18]);在有界域 Dirichlet 边界条件下建立 $L^2$ 弱解及其哈密顿守恒([11]);以及利用凸积分构造非唯一弱解([3,16]),甚至证明在 $C([0,1]; L^{\bar{p}})$ 中($\bar{p} = 4/3 + 10^{-5}$)存在非唯一性([2])。另一方面,[13,14] 利用 vanishing viscosity 方法在环面 $\mathbb{T}^2$ 上对 $L^{4/3}$ 和 Sobolev 初值构造了哈密顿守恒解。然而,$L^{4/3} = L^{4/3,4/3} \subsetneq L^{4/3,2}$,故 $L^{4/3,2}$ 允许更粗糙的初值(如具有更强奇异性或更慢衰减的函数)。但这一提升带来根本性困难:(1) $L^{4/3,2}$ 中的有界序列可能在点处集中(concentrate),弱收敛不足以保证哈密顿能量的下半连续性;(2) 非线性项 $u \cdot \nabla \theta$ 在 $\theta \in L^{4/3,2}$ 下无法直接定义,因 $u \sim \nabla^\perp \Lambda^{-1} \theta$ 仅属于 $L^{4,2}$(由 Lorentz Hölder 不等式),乘积 $u \theta$ 缺乏自然意义;(3) 标准紧性方法(如 Aubin–Lions)在端点 Lorentz 空间中失效。本文旨在填补这一空白:在最粗的、能控制哈密顿能量的函数空间 $L^{4/3,2}$ 中,对任意初值 $\theta_0 \in L^{4/3,2}(\Omega)$,构造一个全局弱解,并严格证明其哈密顿能量守恒。

2. 主要定理或主要结果

定理 1.1(临界 Lorentz 正则性下的弱解):设 $\theta_0 \in L^{4/3,2}(\Omega)$,其中 $\Omega = \mathbb{R}^2$ 或 $\Omega \subset \mathbb{R}^2$ 为光滑有界区域。则存在一个全局弱解 $\theta$ 满足:

(i) 正则性:$\theta \in L^\infty_{\text{loc}}([0,\infty); L^{4/3,2}(\Omega)) \cap C([0,\infty); \mathcal{D}'(\Omega))$,且在弱拓扑下连续,即 $\theta \in C([0,\infty); X_w)$,其中 $X = L^{4/3,2}(\Omega)$;

(ii) 高振幅截断模控制:对任意 $M \geq 0$ 及 $t \geq 0$,有

$$ \tau_M(\theta(t)) \leq \tau_M(\theta_0), \quad \text{其中 } \tau_M(f) := \|(|f| - M)_+\|_{L^{4/3,2}}; $$

(iii) 加权分布函数估计:存在仅依赖于 $\theta_0$ 的非负非减函数 $\Phi = \Phi_{\theta_0}$,满足 $\Phi(\lambda) \to \infty$(当 $\lambda \to \infty$),且 $1 + \Phi(2\lambda) \leq 2(1 + \Phi(\lambda))$(即 $1+\Phi$ 是 doubling 权),使得对任意 $t \geq 0$,

$$ \int_0^\infty \lambda (1 + \Phi(\lambda)) \, d\theta(t)(\lambda)^{3/2} \, d\lambda \lesssim \int_0^\infty \lambda \, d\theta_0(\lambda)^{3/2} \, d\lambda \simeq \|\theta_0\|_{L^{4/3,2}}^2; $$

(iv) 哈密顿能量守恒:对任意 $t \geq 0$,

$$ \|\theta(t)\|_{\dot{H}^{-1/2}(\Omega)} = \|\theta_0\|_{\dot{H}^{-1/2}(\Omega)}, \quad \text{即 } \int_\Omega \theta(t) \, \Lambda^{-1} \theta(t) \, dx = \int_\Omega \theta_0 \, \Lambda^{-1} \theta_0 \, dx. $$

此外,Proposition 2.1 给出关键的端点映射估计:对 $\Omega = \mathbb{R}^2$ 或有界光滑域,有

$$ \|\Lambda^{-1/2} f\|_{L^2(\Omega)} \lesssim \|f\|_{L^{4/3,2}(\Omega)} \simeq \left( \int_0^\infty \lambda \, d_f(\lambda)^{3/2} \, d\lambda \right)^{1/2}, $$

且该估计在次指数 $q=2$ 处最优:对任意 $q>2$,$\|\Lambda^{-1/2} f\|_{L^2} \lesssim \|f\|_{L^{4/3,q}}$ 不成立。

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

本文的核心技术障碍在于:如何在 $L^{4/3,2}$ 这一端点 Lorentz 空间中克服浓度效应(concentration)与对角线奇异性(diagonal singularity),从而实现近似解到弱解的紧性收敛并保持哈密顿守恒。作者的真正创新点并非引入新工具,而是将多个经典分析对象以一种前所未有的方式组合、精炼并赋予新的功能:

  1. 高振幅截断模 $\tau_M(f)$ 的系统性运用:定义 $\tau_M(f) = \|(|f| - M)_+\|_{L^{4/3,2}}$,它不仅是范数的截断,更是一个兼具凸性、弱下半连续性、正算子收缩性(Lemma 2.1(ii))和测度保持不变性(Lemma 2.1(iii))的泛函。在近似方案中,由于近似解 $\theta_{\epsilon_1,\epsilon_2}(t)$ 的分布函数完全继承自热正则化初值 $e^{\epsilon_2 \Delta} \theta_0$,利用 Lemma 2.1 对热半群和保测度流的两次应用,立即得到一致估计 $\sup_{\epsilon_1,\epsilon_2,t} \tau_M(\theta_{\epsilon_1,\epsilon_2}(t)) \leq \tau_M(\theta_0)$。这为控制大值部分提供了统一、稳健的框架,是避免浓度的关键。

  2. 加权分布函数估计与 $\Phi$-函数构造(Proposition 2.2):仅靠 $\tau_M$ 不足以处理小值部分的弱紧性。作者构造了一个显式的 doubling 权 $\Phi_{\theta_0}(\lambda)$,其定义依赖于 $\theta_0$ 的分布函数尾部衰减速率:选取递增序列 $R_k \to \infty$ 使得 $\sup_{h \in A_f} \int_{R_k}^\infty \lambda \, d_h(\lambda)^{3/2} \, d\lambda \leq 2^{-k} \int_0^\infty \lambda \, d_{\theta_0}(\lambda)^{3/2} \, d\lambda$,再令 $\Phi_{\theta_0}(\lambda) = \sum_{k=1}^\infty \mathbf{1}_{[R_k,\infty)}(\lambda)$。该构造确保了加权积分 $M_{\Phi}(h) = \int_0^\infty \lambda (1+\Phi(\lambda)) \, d_h(\lambda)^{3/2} \, d\lambda$ 在集合 $A_{\theta_0} = \{ h \in L^{4/3,2} : \tau_M(h) \leq \tau_M(\theta_0), \forall M \geq 0 \}$ 上一致有界。这个加权估计是获得 $L^{4/3,2}$ 中相对紧性的核心,它比标准的 $L^p$ 或 Orlicz 空间紧性条件更精细,专门针对端点 Lorentz 空间的结构缺陷而设计。

  3. 为 Lorentz 空间量身定制的近似方案:不同于标准的粘性逼近(vanishing viscosity),本文采用双参数正则化:(i) 截断热表示 $\Lambda^{-1}_{\epsilon_1} f = c_0 \int_{\epsilon_1}^\infty e^{\rho \Delta} f \, \rho^{-1/2} \, d\rho$,使速度 $u_{\epsilon_1,\epsilon_2} = \nabla^\perp \Lambda^{-1}_{\epsilon_1} \theta_{\epsilon_1,\epsilon_2}$ 光滑且有界;(ii) 独立热正则化初值 $\theta_{\epsilon_1,\epsilon_2}(0) = e^{\epsilon_2 \Delta} \theta_0$。此方案的关键优势在于:$\Lambda^{-1}_{\epsilon_1}$ 是自伴、非负算子,从而对近似解有精确的哈密顿恒等式 $\frac{d}{dt} H_{\epsilon_1}(\theta_{\epsilon_1,\epsilon_2}) = 0$,且 $H_{\epsilon_1}(\theta_{\epsilon_1,\epsilon_2}(t)) = \langle e^{\epsilon_2 \Delta} \theta_0, \Lambda^{-1}_{\epsilon_1} e^{\epsilon_2 \Delta} \theta_0 \rangle$。当 $\epsilon_1, \epsilon_2 \to 0$ 时,右侧收敛到 $\|\theta_0\|_{\dot{H}^{-1/2}}^2$,这是最终守恒律的基石。

  4. 对称化弱公式与对角线奇异性控制:为定义非线性项,作者采用对称化技巧:对测试函数 $\phi$,定义

    $$ N_\Omega(\theta, \phi) = \frac{1}{2} \iint_{\Omega \times \Omega} \left[ K_\Omega(x,y) \cdot \nabla \phi(x) + K_\Omega(y,x) \cdot \nabla \phi(y) \right] \theta(x) \theta(y) \, dx dy, $$

    其中 $K_\Omega$ 是 $\nabla^\perp \Lambda^{-1}$ 的核。在 $\mathbb{R}^2$ 中,$K_{\mathbb{R}^2}(x-y) = c (x-y)^\perp / |x-y|^3$,对称化后产生因子 $\nabla \phi(x) - \nabla \phi(y)$,将奇异性从 $|x-y|^{-2}$ 降为 $|x-y|^{-1}$。利用 Riesz 势估计 $\iint \frac{|\theta(x)| |\theta(y)|}{|x-y|} \, dx dy \lesssim \|\Lambda^{-1/2} |\theta|\|_{L^2}^2 \lesssim \|\theta\|_{L^{4/3,2}}^2$,从而 $|N_\Omega(\theta, \phi)| \lesssim \|\nabla^2 \phi\|_{L^\infty} \|\theta\|_{L^{4/3,2}}^2$。该估计的成立,正是依赖于前述 $\tau_M$ 控制所提供的对大值部分的精细裁剪,确保了对角线附近的积分收敛。

4. 与之前工作的比较

本文与既有工作的比较需从多个维度展开:

  1. 与 $L^p$ 类结果([11,18,20])的比较:这些工作均在 $L^p$ ($p>1$) 或 $L^2$ 空间中构造弱解,其核心依赖于 $L^p$ 空间中标准的紧性引理(如 Banach–Alaoglu 定理)和 $L^p$ 范数守恒。然而,$L^p$ 空间在临界尺度上无法控制哈密顿能量:$\Lambda^{-1/2}: L^p \to L^2$ 仅在 $p > 4/3$ 时有界,而在 $p = 4/3$ 时失败。本文将初值空间推进至 $L^{4/3,2}$,这是 $L^{4/3}$ 的真扩张($L^{4/3} = L^{4/3,4/3} \subsetneq L^{4/3,2}$),从而容纳了更多样化的粗糙初值,如 $|x|^{-3/2} \log^{-\alpha}(|x|)$(当 $\alpha > 1/2$ 时属于 $L^{4/3,2}$ 但不属于 $L^{4/3}$)。代价是放弃了 $L^p$ 范数本身的守恒,转而依赖更精细的 $\tau_M$ 和加权估计来维持紧性。

  2. 与 vanishing viscosity 方法([13,14])的比较:[13,14] 在环面 $\mathbb{T}^2$ 上对 $L^{4/3}$ 初值使用粘性逼近 $\partial_t \theta_\nu + u_\nu \cdot \nabla \theta_\nu + \nu \Lambda^\alpha \theta_\nu = 0$($\alpha > 0$)。其成功依赖于粘性项提供的额外正则性和能量耗散,从而获得 $\theta_\nu$ 在更高阶空间中的有界性。本文则在完全无粘($\nu = 0$)的情形下工作,通过截断 $\Lambda^{-1}$ 而非添加 $\Lambda^\alpha$ 项来正则化速度,这避免了任何能量耗散,从而得以严格证明哈密顿守恒。这是一个本质区别:[13,14] 的解是粘性解的极限,其能量可能随 $\nu \to 0$ 而损失;而本文的解是纯传输解的极限,能量严格不变。

  3. 与凸积分非唯一性结果([2,3,16])的比较:[2] 在略高于临界的空间 $C([0,1]; L^{\bar{p}})$($\bar{p} = 4/3 + 10^{-5}$)中构造了非唯一弱解,这表明在稍高正则性下,解的刚性已不足。本文的结果与之形成鲜明互补:它在真正的临界端点 $L^{4/3,2}$ 上构造了一个特定的、哈密顿守恒的解。这暗示了在端点处,解的“刚性”与“柔性”可能共存——守恒解存在,但可能不唯一。本文并未解决唯一性问题,但从方法论上,其 $\tau_M$ 控制和加权估计为未来研究该临界点的适定性边界提供了关键工具。

  4. 与有界域处理([11])的比较:[11] 在有界域 Dirichlet 边界条件下建立了 $L^2$ 弱解。本文将其推广至 $L^{4/3,2}$,并发展了一套适用于有界域的对称化弱公式(Section 4.2),其核心是利用 Dirichlet 热核 $H_D(t,x,y)$ 的对称性 $H_D(t,x,y) = H_D(t,y,x)$ 来保证 $G_{1,\epsilon}(x,y) = G_{1,\epsilon}(y,x)$,从而使得对称化核 $K_\epsilon(x,y) \cdot \nabla \phi(x) + K_\epsilon(y,x) \cdot \nabla \phi(y)$ 具有相同的 $|x-y|^{-1}$ 奇异性。这表明本文的方法具有良好的几何适应性,不局限于全空间。

综上,本文的推进在于:将初值正则性降至哈密顿能量所允许的最粗端点 $L^{4/3,2}$,在无粘前提下严格证明能量守恒,并为此发展了一套以 $\tau_M$ 和加权分布函数为核心的、专为端点 Lorentz 空间设计的紧性分析框架。


5. 由输运噪声驱动的随机相分离

Stochastic phase separation driven by transport noise

作者: Andrea Di Primio, Andrea Papini, Luca Scarpa

类别: math.AP, math.PR

发布日期: 2026-07-27

arXiv ID: 2607.24326v1

摘要翻译: 复杂流体中发生的相分离现象已知对流体力学效应敏感:当二元流体混合物经历旋节分解时,所形成的畴结构不再仅由界面动力学决定,而是被对流作用重塑。在若干感兴趣的应用中,对流效应甚至可能引发湍流,从而在宽广的空间尺度范围内进一步影响分离过程。受此启发,我们提出了一种由输运噪声驱动的 Cahn-Hilliard 模型,以捕捉湍流固有的随机性。该模型在其热力学相关框架下进行分析,即采用奇异 Flory-Huggins 势和可能退化的迁移率,并且噪声同时考虑 Itô 形式与 Stratonovich 形式。在本工作中,作为初步研究,我们在二维与三维空间中建立了该系统的适定性及其热力学一致性。

Comment: 49 pp

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

1. 研究对象和背景

本文研究如下两类随机Cahn–Hilliard方程:在二维或三维有界Lipschitz区域 $\Omega \subset \mathbb{R}^d$($d=2,3$)上,以$\varphi: [0,T] \times \Omega \to [-1,1]$为未知相场变量,其演化由以下带输运型噪声的四阶退化抛物系统描述:

$$ \begin{cases} d\varphi - \operatorname{div}\big(m(\varphi)\nabla\mu\big)\,dt + \kappa L_{\text{Strat}}\varphi\,dt = \nabla\varphi \cdot \Sigma\, dW & \text{in } (0,T) \times \Omega, \\ \mu = -\Delta\varphi + F'(\varphi) & \text{in } (0,T) \times \Omega, \\ \partial_n \varphi = m(\varphi)\nabla\mu \cdot n = 0 & \text{on } (0,T) \times \partial\Omega, \\ \varphi(0) = \varphi_0 & \text{in } \Omega, \end{cases} \tag{1.8} $$

其中$\kappa \in \{0,1\}$编码噪声类型:$\kappa=0$对应Itô型输运噪声(即$d\varphi = \operatorname{div}(m(\varphi)\nabla\mu)\,dt + \nabla\varphi \cdot \Sigma\, dW$),$\kappa=1$对应Stratonovich型输运噪声(即$d\varphi = \operatorname{div}(m(\varphi)\nabla\mu)\,dt + \nabla\varphi \cdot \Sigma \circ dW$);$L_{\text{Strat}}$为经典Stratonovich校正算子,显式为$L_{\text{Strat}}\varphi = \frac{1}{2}\sum_{k\in\mathbb{N}} \operatorname{div}\big( (\sigma_k \cdot \nabla)\nabla\varphi + \nabla\varphi \cdot D\sigma_k + \sigma_k (\Delta\varphi) \big)$,其结构依赖于$\Sigma = \sum_k \sigma_k u_k$的Hilbert–Schmidt表示;$W$为定义在可分Hilbert空间$U$上的柱状Wiener过程;$\Sigma: \Omega \times [0,T] \to \mathcal{L}_2(U; W^r)$为线性Hilbert–Schmidt算子,其值域位于满足$\operatorname{div}\sigma_k = 0$且$\sigma_k \cdot n|_{\partial\Omega}=0$的Sobolev空间$W^r := \{ u \in H^r(\Omega)^d : \operatorname{div} u = 0,\, u \cdot n = 0 \text{ on } \partial\Omega \}$中,$r > r^* = d/2$;$m: [-1,1] \to [0,+\infty)$为迁移率函数,允许退化情形(如$m_{\text{deg}}(x) = 1-x^2$);$F: (-1,1) \to [0,+\infty)$为奇异Flory–Huggins势,形如$F_{\log}(x) = \frac{\theta}{2}[(1+x)\log(1+x)+(1-x)\log(1-x)] - \frac{\theta_0}{2}x^2$,满足$0 < \theta < \theta_0$,并延拓为真凸下半连续函数$F: \mathbb{R} \to [0,+\infty]$,使得$F(r)=+\infty$当$r \notin [-1,1]$;初始数据$\varphi_0 \in H^1(\Omega)$满足$\|\varphi_0\|_{L^\infty} \leq 1$(物理约束)。

该模型的物理动机源于复杂流体中相分离对流体动力学效应的高度敏感性:当二元流体混合物经历旋节分解(spinodal decomposition)时,畴结构不再仅由界面动力学决定,而被对流显著重塑;在强湍流情形下,小尺度速度场的混沌性导致相场受到高频、白噪声型输运扰动。传统确定性模型(如带对流项的Cahn–Hilliard方程$\partial_t \varphi = \operatorname{div}(m(\varphi)\nabla\mu) + \nabla\varphi \cdot u$)将$u$视为给定或耦合求解的确定性速度场,但无法刻画湍流固有的随机性与多尺度特性。本文采用“湍流输运的随机建模”范式:通过Wong–Zakai逼近论证,当小尺度速度场$u^\varepsilon$的时间相关性在$\varepsilon \to 0$极限下退化时,其对相场的作用收敛于Stratonovich输运噪声$\nabla\varphi \cdot u \circ \dot{W}$。这一建模不仅物理自洽(直接源于湍流多尺度能量级联),且保持原系统的两大核心结构:质量守恒(因$\sigma_k$无散且切向边界条件保证$\int_\Omega \nabla\varphi \cdot \sigma_k \,dx = 0$)与能量耗散结构(虽含噪声,仍可建立随机能量不等式)。然而,现有文献完全缺失对此类输运噪声模型的严格分析:已知结果集中于加性噪声(Cook模型)或非输运型乘性噪声(如保守型$\varphi\,dW$),而对$\nabla\varphi \cdot \Sigma\,dW$形式,尤其在$F$奇异、$m$退化的热力学相关框架下,尚无任何适定性理论。本文旨在填补这一根本性空白——建立该模型在$\mathbb{R}^2$与$\mathbb{R}^3$中的概率弱解存在性,并验证其热力学一致性(即质量守恒与随机能量耗散律),从而为随机相分离的严格数学理论奠基。

2. 主要定理或主要结果

本文在$\Omega \subset \mathbb{R}^d$($d=2,3$)上建立了模型(1.8)的概率弱解适定性理论,主要结果分为非退化与退化迁移率两类:

(i)非退化迁移率情形(假设(A2)ₙd成立):对任意$\kappa \in \{0,1\}$及满足(A3)ₙd的噪声系数$\Sigma$(即存在$r > d/2$使得$\sup_{t\in[0,T]} \sum_k \|\sigma_k(t)\|_H^2 \in L^{q/2}(\Omega)$且$\exp\big(\int_0^T (\sum_k \|\sigma_k(t)\|_{W^r}^2)^2 dt\big) \in L^\ell(\Omega),\, \forall \ell \geq 1$),对任意初值$\varphi_0 \in V_1 = H^1(\Omega)$满足$\|\varphi_0\|_{L^\infty} \leq 1$,存在一个概率空间$(\widetilde{\Omega}, \widetilde{\mathcal{F}}, \widetilde{\mathbb{P}})$、其上定义的滤链$(\widetilde{\mathcal{F}}_t)_t$、柱状Wiener过程$\widetilde{W}$及过程$(\widetilde{\varphi}, \widetilde{\mu})$,使得$(\widetilde{\Omega}, \widetilde{\mathcal{F}}, \widetilde{\mathbb{P}}, (\widetilde{\mathcal{F}}_t)_t, \widetilde{W}, \widetilde{\varphi}, \widetilde{\mu})$构成模型(1.8)的概率弱解,且满足正则性

$$ \widetilde{\varphi} \in L^2\big(\widetilde{\Omega}; C^0([0,T]; H) \cap L^2(0,T; V_2)\big),\quad \widetilde{\mu} \in L^2\big(\widetilde{\Omega}; L^2(0,T; V_1)\big), $$

并满足几乎必然的物理约束$\widetilde{\varphi}(t,x) \in [-1,1]$对所有$(t,x)$成立。此外,该解满足质量守恒律$\int_\Omega \widetilde{\varphi}(t) = \int_\Omega \varphi_0$及随机能量不等式:存在常数$C>0$使得 $$ \mathbb{E}\left[ \frac{1}{2}|\nabla\widetilde{\varphi}(t)|H^2 + \int\Omega F(\widetilde{\varphi}(t)) \right]

  • \mathbb{E}\int_0^t |\nabla\widetilde{\mu}(s)|_H^2 ,ds \leq C\left( |\nabla\varphi_0|H^2 + \int\Omega F(\varphi_0) + t \right). $$

(ii)退化迁移率情形(假设(A2)ₑ与(A3)ₑ成立):当迁移率$m$在$\pm1$处退化(如$m_{\text{deg}}(x)=1-x^2$)且噪声满足小ness条件$\sum_k \|\sigma_k(t)\|_{W^r}^2 \leq C_\Sigma^2$ a.e.,其中$C_\Sigma < 2k_m / C_\infty^2$($k_m$来自(A2)ₑ中$m\Psi'' \geq k_m$,$C_\infty$为$H^r \hookrightarrow L^\infty$嵌入常数),则对任意初值$\varphi_0 \in V_1$满足$\|\varphi_0\|_{L^\infty} \leq 1$,存在概率弱解$(\widetilde{\varphi}, \widetilde{\mu})$,其正则性为

$$ \widetilde{\varphi} \in L^2\big(\widetilde{\Omega}; C^0_w([0,T]; V_1) \cap L^2(0,T; V_2)\big),\quad \widetilde{\mu} \in L^2\big(\widetilde{\Omega}; L^2(0,T; V_1)\big), $$

且满足$\widetilde{\varphi}(t,x) \in [-1,1]$ a.s.,以及改进的能量估计:存在$\alpha>0$使得 $$ \mathbb{E}\left[ \frac{1}{2}|\nabla\widetilde{\varphi}(t)|H^2 + \int\Omega F(\widetilde{\varphi}(t)) \right]

  • \mathbb{E}\int_0^t \int_\Omega m(\widetilde{\varphi}) |\nabla\widetilde{\mu}|^2 ,dx,ds \leq C\left( |\nabla\varphi_0|H^2 + \int\Omega F(\varphi_0) \right) e^{\alpha t}. $$

(iii)路径唯一性:当迁移率为常数(即$m \equiv m_0 > 0$)时,在非退化情形下,对任意满足(A3)ₙd的$\Sigma$,概率弱解具有路径唯一性(pathwise uniqueness),从而由Yamada–Watanabe定理推出强解存在性,即模型(1.8)在常迁移率情形下是完全适定的(well-posed)。该结论不依赖于$\kappa$取值,故同时覆盖Itô与Stratonovich两种噪声形式。

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

证明的核心障碍在于三重耦合:(a) 奇异势$F$强制解保持在$[-1,1]$内,导致$F'$在端点发散,传统Galerkin逼近需处理非线性项的紧性缺失;(b) 输运噪声$\nabla\varphi \cdot \Sigma\,dW$是梯度型乘性噪声,其Itô–Stratonovich转换引入高阶校正项$L_{\text{Strat}}\varphi$,该算子含二阶导数与$D\sigma_k$,破坏标准能量估计的封闭性;(c) 退化迁移率$m(\varphi)$在$\pm1$处 vanish,使扩散系数退化,导致解的正则性下降,且$\mu = -\Delta\varphi + F'(\varphi)$中$F'$的奇异性与$m$的零点相互作用,形成双重退化结构。

作者的真正创新点在于构建一套分层逼近与补偿紧性框架,具体包含以下关键技术:

第一,双正则化逼近方案:为同时处理$F$的奇异性与$m$的退化性,作者引入双重正则化参数$(\varepsilon, \delta)$。首先用光滑截断势$F_\varepsilon(r) = F(\rho_\varepsilon(r))$($\rho_\varepsilon$为$[-1,1]$上的光滑截断)替代$F$,再用非退化迁移率$m_\delta(r) = m(r) + \delta$替代$m$,得到正则化系统

$$ d\varphi_{\varepsilon,\delta} = \operatorname{div}\big(m_\delta(\varphi_{\varepsilon,\delta})\nabla\mu_{\varepsilon,\delta}\big)\,dt + \kappa L_{\text{Strat}}^{\delta}\varphi_{\varepsilon,\delta}\,dt + \nabla\varphi_{\varepsilon,\delta} \cdot \Sigma_\delta\,dW, \quad \mu_{\varepsilon,\delta} = -\Delta\varphi_{\varepsilon,\delta} + F_\varepsilon'(\varphi_{\varepsilon,\delta}). $$

此系统在Galerkin框架下可证强解存在唯一,关键在于利用$F_\varepsilon'$的有界性与$m_\delta$的一致正定性获得先验估计。

第二,统一的随机能量-熵联合估计:针对输运噪声,作者发展了一种新型估计技术——将经典能量估计与对数型熵估计耦合。对正则化解$\varphi_{\varepsilon,\delta}$,除标准能量不等式外,还建立熵估计 $$ \mathbb{E}\int_\Omega \Psi_\varepsilon(\varphi_{\varepsilon,\delta}(t))

  • \mathbb{E}\int_0^t \int_\Omega m_\delta(\varphi_{\varepsilon,\delta}) |\nabla F_\varepsilon’(\varphi_{\varepsilon,\delta})|^2 ,dx,ds \leq C\left( \int_\Omega \Psi_\varepsilon(\varphi_0) + \mathbb{E}\int_0^t |\nabla\varphi_{\varepsilon,\delta}|_H^2 ,ds \right), $$ 其中$\Psi_\varepsilon$为$\Psi$的光滑逼近。该估计的关键在于利用$\nabla\varphi \cdot \Sigma\,dW$的散度结构,通过Itô公式展开$\int_\Omega \Psi_\varepsilon(\varphi)$时,噪声项贡献为$\sum_k \int_0^t \int_\Omega \Psi_\varepsilon''(\varphi) |\nabla\varphi \cdot \sigma_k|^2 \,dx\,ds$,其符号由$\Psi_\varepsilon'' > 0$保证非负,从而可被耗散项控制。此技术绕开了$L_{\text{Strat}}$的复杂结构,直接利用噪声的几何特性。

第三,基于广义Korn不等式的补偿紧性:在退化情形下,$m_\delta \to m$导致$\nabla\mu_{\varepsilon,\delta}$在$L^2$中弱收敛失效。作者引入广义Korn不等式:对任意$\psi \in V_2$,有$\|\nabla\psi\|_{L^2}^2 \leq C \int_\Omega m(\psi) |\nabla(-\Delta\psi + F'(\psi))|^2 \,dx$,该不等式由假设(A2)ₑ中$m\Psi'' \geq k_m$保证。结合熵估计,可证得$\{\varphi_{\varepsilon,\delta}\}$在$C^0_w([0,T]; V_1)$中弱-*紧,且$\{m_{\delta}(\varphi_{\varepsilon,\delta})\nabla\mu_{\varepsilon,\delta}\}$在$L^2(0,T; H)$中弱紧。最终通过单调算子理论与Minty技巧识别极限,证明$\mu = -\Delta\varphi + F'(\varphi)$在分布意义下成立。

第四,Stratonovich校正项的精细处理:对$L_{\text{Strat}}\varphi$,作者将其分解为三部分:$L_{\text{Strat}}^1\varphi = \frac{1}{2}\sum_k \operatorname{div}((\sigma_k \cdot \nabla)\nabla\varphi)$,$L_{\text{Strat}}^2\varphi = \frac{1}{2}\sum_k \operatorname{div}(\nabla\varphi \cdot D\sigma_k)$,$L_{\text{Strat}}^3\varphi = \frac{1}{2}\sum_k \sigma_k (\Delta\varphi)$。利用$\sigma_k \in W^r$($r>d/2$)的Sobolev嵌入性质($W^r \hookrightarrow L^\infty$),证明$L_{\text{Strat}}^1$与$L_{\text{Strat}}^2$属于$V_1^*$,而$L_{\text{Strat}}^3$属于$H$,从而整体属于$V_1^*$,确保校正项在弱解框架下可积。

4. 与之前工作的比较

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

(1)噪声结构的根本性差异:现有文献(如[27,29,31,35,67,90,92])均处理加性噪声($d\varphi = \cdots + \Sigma\,dW$)或非输运型乘性噪声(如$\varphi\,dW$或$\operatorname{div}(\varphi\,\Sigma\,dW)$),其数学性质与物理含义迥异。加性噪声破坏质量守恒(因$\int_\Omega \Sigma\,dW \neq 0$),非输运乘性噪声不具散度自由结构,无法反映湍流对流的本质。本文首次系统分析输运噪声$\nabla\varphi \cdot \Sigma\,dW$,其核心优势在于:(i) 自动保持质量守恒(因$\int_\Omega \nabla\varphi \cdot \sigma_k \,dx = 0$);(ii) 在Stratonovich意义下与确定性对流项$\nabla\varphi \cdot u$连续,为湍流建模提供严格基础。此前唯一相关工作[42]仅作数值探索,且限于多项式势与常迁移率,无理论分析。

(2)热力学框架的严格实现:既有随机Cahn–Hilliard研究(如[29,31,67,90])在奇异势下多采用反射测度或变分不等式方法,但未处理输运噪声;而[35,90]在退化迁移率下考虑保守噪声,其噪声形式为$\operatorname{div}(\Sigma\,dW)$,与本文的$\nabla\varphi \cdot \Sigma\,dW$结构不同,前者不依赖于相场梯度,无法体现对流效应。本文首次在$F$奇异、$m$退化的完整热力学框架下,证明输运噪声模型的适定性,并验证其能量耗散律与质量守恒律,填补了该交叉领域的理论空白。

(3)时间尺度与适用范围的推进:在确定性退化Cahn–Hilliard中([38,39]),解的存在性已建立,但唯一性仍是开放问题(因四阶退化抛物性导致解可能不唯一)。本文在常迁移率情形下证明路径唯一性,将适定性结果从存在性提升至强解层面;而在退化情形下,通过小ness条件$C_\Sigma < 2k_m / C_\infty^2$,将噪声强度控制在能量耗散可主导的范围内,从而突破确定性理论的局限。该条件物理上对应湍流强度低于临界阈值,此时噪声不破坏相分离的热力学驱动力,而仅调制其动力学。

(4)代价与局限:本文的代价在于对噪声正则性的较高要求:非退化情形需$\sigma_k \in L^4(0,T; W^r)$($r>d/2$),退化情形需$\sigma_k \in L^\infty(0,T; W^r)$且满足小ness条件。这比标准Navier–Stokes解的正则性更强(如在$d=3$时需$u \in L^\infty(0,T; H^{3/2+})$,强于经典强解的$u \in L^\infty(0,T; H^1) \cap L^2(0,T; H^2)$)。此外,本文未解决退化迁移率下的路径唯一性,亦未处理主动相分离(即$\Sigma$依赖于解本身),这些是未来工作的核心挑战。作者明确指出,主动情形将导致解依赖型输运噪声,属未探索领域,其分析需全新工具。


6. 时间精化Triebel–Lizorkin估计及其在Keller–Segel型方程中的应用

Time-refined Triebel–Lizorkin Estimates and Applications to Keller–Segel Type Equations

作者: Wenhai Shan, Xiao-song Yang

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24089v1

摘要翻译: 我们发展了临界齐次Triebel–Lizorkin空间中的热流估计,并将其应用于二维Keller–Segel型方程。对 $q>2$,我们证明热流从 $\dot F^0_{1,q}(\mathbb R^2)$ 到 $L^2(0,T;\dot F^1_{1,q}(\mathbb R^2))$ 是无界的。这促使我们引入一类新空间——时间精化Triebel–Lizorkin空间,其中时间范数取在二进制求和之前。在此框架下,我们建立了一族热光滑估计及相应的极大正则性估计。我们进一步证明了Banach值Peetre估计、Banach值Jawerth型估计、齐次Poisson估计,以及Keller–Segel漂移项的端点双线性估计。作为应用,我们得到了二维抛物–椭圆Keller–Segel方程在临界空间 $\dot F^0_{1,q}(\mathbb R^2)$($2\le q<\infty$)中对任意初值的局部适定性,以及对充分小初值的全局适定性。同样的解析框架亦可推出一类相关两组分漂移–扩散系统的类似适定性结果。

Comment: 无

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

1. 研究对象和背景

本文研究二维抛物-椭圆Keller–Segel方程

$$ \partial_t u - \Delta u + \nabla \cdot \big( u \nabla (-\Delta)^{-1} u \big) = 0, \quad u|_{t=0} = u_0 \tag{1.1} $$

在空间端点 $p = 1$ 的临界正则性框架下的适定性问题。该方程建模细胞趋化聚集过程,其数学核心在于扩散项 $-\Delta u$ 与非线性漂移项 $\nabla \cdot (u \nabla (-\Delta)^{-1} u)$ 之间的竞争:前者抑制集中,后者通过自诱导势场驱动质量坍缩。方程具有抛物标度不变性 $u_\lambda(t,x) = \lambda^2 u(\lambda^2 t, \lambda x)$,由此导出临界函数空间需满足 $\|u_0\|_{X} \sim \|u_{0,\lambda}\|_{X}$。在 $\mathbb{R}^2$ 中,对应于 $L^1$ 空间尺度的齐次Triebel–Lizorkin临界空间为 $\dot F^0_{1,q}(\mathbb{R}^2)$,其中 $2 \le q < \infty$;当 $q=2$ 时,$\dot F^0_{1,2}(\mathbb{R}^2) \cong H^1(\mathbb{R}^2)$(Hardy空间),已有Ogawa–Shimizu在 $H^1$ 中建立大初值局部适定性与小初值整体适定性。然而,对 $q > 2$,标准Bochner空间框架 $L^2(0,T; \dot F^1_{1,q}(\mathbb{R}^2))$ 下的热流估计失效——Proposition 1.1严格证明:对任意 $q > 2$、$T > 0$,不存在常数 $C > 0$ 使得

$$ \|e^{t\Delta} u_0\|_{L^2(0,T; \dot F^1_{1,q})} \le C \|u_0\|_{\dot F^0_{1,q}} \quad \forall u_0 \in \dot F^0_{1,q}(\mathbb{R}^2). $$

其根源在于标准Bochner范数中时间 $L^2$ 积分与频率 $\ell^q$ 求和的不可交换性:对由 $N$ 个频带分离的初值块 $f_{k_m}$ 构成的 $u_0^N = \sum_{m=1}^N f_{k_m}$,其 $\dot F^0_{1,q}$ 范数为 $\|u_0^N\|_{\dot F^0_{1,q}} \sim N^{1/q} \|\phi\|_{L^1}$,而热流在各频带上的贡献在时间上近似正交,导致 $\|e^{t\Delta} u_0^N\|_{L^2(0,T; \dot F^1_{1,q})} \gtrsim N^{1/2} \|\phi\|_{L^1}$,从而要求 $N^{1/2} \lesssim N^{1/q}$,这在 $q > 2$ 时不可能成立。因此,现有基于Besov空间(如Iwabuchi、Zhao在 $\dot B^0_{1,\rho}$ 中的工作)或Hardy空间($q=2$)的临界理论无法覆盖 $q > 2$ 的Triebel–Lizorkin端点尺度。本文旨在填补这一空白:构造一种新型函数空间框架,使热流估计在 $\dot F^0_{1,q}$ 初值下对所有 $2 \le q < \infty$ 均成立,并由此建立统一的适定性理论。

2. 主要定理或主要结果

本文核心结果是建立时间精炼Triebel–Lizorkin空间中的热流估计与Keller–Segel方程的临界适定性。首先,定义时间精炼范数:对 $1 \le q < \infty$、$1 \le \rho \le \infty$,

$$ \|u\|_{\dot F^s_{1,q}(L^\rho_T)} := \left\| \left( 2^{js} \|\dot \Delta_j u(\cdot, x)\|_{L^\rho(0,T)} \right)_{j \in \mathbb{Z}} \right\|_{L^1_x \ell^q_j}. $$

该范数将时间 $L^\rho$ 范数置于每个频带 $\dot \Delta_j$ 内部,再进行空间 $L^1_x$ 与频率 $\ell^q_j$ 的合成,从而规避Proposition 1.1的障碍。在此框架下,Theorem 1.2给出广义热估计:设 $1 \le q < \infty$,$1 \le a \le r \le \infty$,$s \in \mathbb{R}$,令 $\sigma = 2\left(1 + \frac{1}{r} - \frac{1}{a}\right)$,若 $u$ 满足 $\partial_t u - \Delta u = F$,$u|_{t=0} = u_0$,且 $u_0 \in \dot F^{s+2-2/a}_{1,q}$,$F \in \dot F^s_{1,q}(L^a_T)$,则

$$ \|u\|_{\dot F^{s+\sigma}_{1,q}(L^r_T)} \lesssim \|u_0\|_{\dot F^{s+2-2/a}_{1,q}} + \|F\|_{\dot F^s_{1,q}(L^a_T)}, $$

$$ \|\Delta u\|_{\dot F^{s+\sigma-2}_{1,q}(L^r_T)} \lesssim \|u_0\|_{\dot F^{s+2-2/a}_{1,q}} + \|F\|_{\dot F^s_{1,q}(L^a_T)}, $$

$$ \partial_t u \in \dot F^{s+\sigma-2}_{1,q}(L^r_T) + \dot F^s_{1,q}(L^a_T), $$

其中隐含常数独立于 $T$。特别地,取 $(a,r,s) = (1,\infty,0)$ 与 $(a,r,s) = (1,2,0)$,得关键离对角估计:

$$ \left\| \int_0^t e^{(t-\tau)\Delta} F(\tau) \, d\tau \right\|_{\dot F^0_{1,q}(L^\infty_T)} \lesssim \|F\|_{\dot F^0_{1,q}(L^1_T)}, $$

$$ \left\| \int_0^t e^{(t-\tau)\Delta} F(\tau) \, d\tau \right\|_{\dot F^1_{1,q}(L^2_T)} \lesssim \|F\|_{\dot F^0_{1,q}(L^1_T)}. $$

这两个估计直接支撑非线性分析。其次,Proposition 1.3给出Keller–Segel漂移的端点双线性估计:对 $2 \le q < \infty$,

$$ \left\| \nabla \cdot \big( f \nabla (-\Delta)^{-1} g \big) \right\|_{\dot F^0_{1,q}(L^1_T)} \lesssim \|f\|_{\dot F^1_{1,q}(L^2_T)} \|g\|_{\dot F^1_{1,q}(L^2_T)}. $$

结合Theorem 1.2,可控制Duhamel项:

$$ \left\| \int_0^t e^{(t-\tau)\Delta} \nabla \cdot \big( f \nabla (-\Delta)^{-1} g \big)(\tau) \, d\tau \right\|_{X^q_T} \lesssim \|f\|_{\dot F^1_{1,q}(L^2_T)} \|g\|_{\dot F^1_{1,q}(L^2_T)}, $$

其中解空间 $X^q_T := \dot F^0_{1,q}(L^\infty_T) \cap \dot F^1_{1,q}(L^2_T)$。最终,Theorem 1.4断言:对任意 $2 \le q < \infty$ 及 $u_0 \in \dot F^0_{1,q}(\mathbb{R}^2)$,存在 $T > 0$,方程(1.1)在 $X^q_T$ 中存在唯一温和解 $u$;且 $u \in C([0,T]; \dot F^0_{1,q}(\mathbb{R}^2))$,数据到解映射连续;进一步,存在 $\varepsilon_q > 0$,若 $\|u_0\|_{\dot F^0_{1,q}} \le \varepsilon_q$,则解整体存在且属于 $X^q_\infty$。

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

本文的核心技术障碍有二:其一,标准Bochner范数下 $\dot F^0_{1,q}$ 初值的热流估计在 $q > 2$ 时失效,源于时间 $L^2$ 积分与频率 $\ell^q$ 求和顺序冲突;其二,在空间端点 $L^1_x$ 下,经典向量值极大函数不适用,导致频率局域化热核估计无法提升至全空间范数。作者的真正创新点在于构造并系统开发“时间精炼Triebel–Lizorkin空间” $\dot F^s_{1,q}(L^\rho_T)$,并建立一套适配该空间的端点调和分析工具链。具体而言:(1) 针对障碍一,作者摒弃标准Bochner框架,转而定义 $\|u\|_{\dot F^s_{1,q}(L^\rho_T)}$ 将时间范数嵌入每个频带内部,使时间积分与频带自然热时间尺度 $t \sim 2^{-2j}$ 同步,从而恢复热流的逐频带衰减结构。(2) 针对障碍二,作者发展Banach值Peetre极大估计(Lemma 2.1):对E值函数 $u_j$ 满足 $\operatorname{supp} \hat u_j \subset \{ |\xi| \le C 2^j \}$,有

$$ u^*_j(x) := \sup_{y \in \mathbb{R}^2} \frac{ \|u_j(x-y)\|_E }{ (1 + 2^j |y|)^M } \lesssim \big[ M(\|u_j\|^\theta_E)(x) \big]^{1/\theta}, $$

其中 $0 < \theta < \min\{1,q\}$,$M > 2/\theta$,$M$ 为Hardy–Littlewood极大算子。此估计通过标量化(scalarization)而非直接应用E值极大不等式,绕过 $L^1_x$ 上极大算子无界的困境。(3) 利用该Peetre估计与Fefferman–Stein向量值极大不等式(作用于指数 $>1$ 的 $\ell^{q/\theta}_j L^{1/\theta}_x$ 空间),导出序列估计(Corollary 2.2):

$$ \|(u^*_j)_j\|_{L^1_x \ell^q_j} \lesssim \big\| \big( \|u_j(x)\|_E \big)_j \big\|_{L^1_x \ell^q_j}, $$

并结合频带有限重叠性(Lemma 2.3)实现频带和的重建:

$$ \left\| \sum_k u_k \right\|_{\dot F^0_{1,q}(\mathbb{R}^2; E)} \lesssim \|(u^*_k)_k\|_{L^1_x \ell^q_k}. $$

(4) 对热核 $K_j(t,x)$(限制于环形区域 $|\xi| \sim 2^j$),利用其逐点衰减 $|K_j(t,x)| \lesssim_N 2^{2j} e^{-c t 2^{2j}} (1 + 2^j |x|)^{-N}$,结合上述Peetre估计与时间Young不等式,得到关键离对角估计:对 $1 \le a \le r \le \infty$,

$$ \|\dot \Delta_j e^{t\Delta} u_0\|_{L^r_t} \lesssim 2^{-2j (1/a - 1/r)} \|\dot \Delta_j u_0\|_{L^1_x}, $$

此即Theorem 1.2中 $\sigma$ 参数的来源。最后,对非线性项 $\nabla \cdot (f \nabla (-\Delta)^{-1} g)$,作者采用Bony分解,并分别处理:低-高项 $T_f b$($b = \nabla (-\Delta)^{-1} g$)用Banach值Peetre估计控制;高-低与同频项 $T_b f + R(f,b)$ 则借助Banach值Jawerth型估计(Section 4)将其嵌入 $\dot B^0_{1,1}(L^1_T)$,再利用 $\dot B^0_{1,1} \hookrightarrow \dot F^0_{1,q}$($q \ge 1$)完成闭环。

4. 与之前工作的比较

本文与既有工作的比较需从空间框架、估计强度、适用范围及代价四维度展开。首先,与Besov空间工作(如Iwabuchi [19]、Zhao [33])相比:Iwabuchi在 $\dot B^0_{1,\rho}(\mathbb{R}^2)$($1 < \rho \le \infty$)中建立小初值整体适定性,但 $\dot B^0_{1,\rho}$ 与 $\dot F^0_{1,q}$ 在端点 $p=1$ 本质不同——前者范数为 $\|f\|_{\dot B^0_{1,\rho}} = \|(\dot \Delta_j f)_j\|_{\ell^\rho_j L^1_x}$,频带求和在空间积分之后,故无Proposition 1.1的障碍;Zhao推广至更广Besov指标,但仍限于Besov尺度。本文首次在Triebel–Lizorkin端点 $\dot F^0_{1,q}$($q > 2$)实现突破,且 $\dot F^0_{1,q}$ 严格包含 $\dot B^0_{1,1}$(因 $\ell^1_j \hookrightarrow \ell^q_j$),但被 $\dot B^0_{1,\infty}$ 包含,故本文结果与Besov工作互补而非覆盖。其次,与Hardy空间工作(Ogawa–Shimizu [26])相比:其在 $H^1 = \dot F^0_{1,2}$ 中的结果是本文 $q=2$ 的特例,但本文将 $q$ 推广至 $[2,\infty)$ 全范围,揭示了 $q=2$ 的特殊性源于 $\ell^2$ 与 $L^2_t$ 的Hilbert空间兼容性,而本文通过时间精炼范数消除了这一限制。第三,与UMD型极大正则性理论(如Weis [31])相比:标准理论要求底空间具有UMD性质且算子族R-有界,而 $L^1$ 不具UMD,故无法适用;本文绕过UMD框架,以Peetre估计与标量化机制在端点 $L^1_x$ 实现类似效果,代价是仅适用于Littlewood–Paley可分解空间且估计常数依赖于 $q$。第四,与Guo–Li对Euler方程的工作 [15] 相比:其亦用Peetre估计处理 $\dot F^0_{1,q}$ 端点,但针对稳态问题且未引入时间精炼结构;本文将该机制拓展至Banach值函数空间 $L^\rho_T$,并首创时间精炼范数以解决热流估计的根本障碍,技术深度显著提升。最后,本文方法的局限在于:目前仅适用于二维Keller–Segel及同类漂移-扩散系统;对更高维或更复杂非线性(如抛物-抛物系统)的推广尚待验证;且时间精炼空间的插值理论与乘子刻画尚未建立,制约其在其他PDE中的应用广度。


7. 随机抛物-抛物Keller-Segel系统的全局存在性与路径唯一性

Global Existence and Pathwise Uniqueness for a Stochastic Parabolic-Parabolic Keller-Segel System

作者: Jinhuan Wang, Qian Li, Hui Huang

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24128v1

摘要翻译: 本文研究定义在二维有界区域上的、由非局部非线性乘性噪声驱动的随机抛物-抛物Keller-Segel系统。在适当假设下,我们对任意初值建立了强解的全局存在性与路径唯一性,且未施加任何小性条件。这与确定性二维Keller-Segel系统形成鲜明对比——后者通常要求初值质量满足小性条件。

Comment: 无

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

1. 研究对象和背景

本文研究二维有界光滑区域 $\Omega \subset \mathbb{R}^2$ 上如下随机抛物-抛物Keller-Segel系统:

$$ \begin{cases} du = \big(\Delta u - \nabla \cdot (u \nabla v)\big)\,dt + \sigma(u)\,dW(t), & x \in \Omega,\ t > 0, \\ dv = (\Delta v - v + u)\,dt, & x \in \Omega,\ t > 0, \\ \frac{\partial u}{\partial \nu} = \frac{\partial v}{\partial \nu} = 0, & x \in \partial\Omega,\ t > 0, \\ u(x,0) = u_0(x),\quad v(x,0) = v_0(x), & x \in \Omega, \end{cases} \tag{1.3} $$

其中 $u = u(x,t)$ 表示细菌密度,$v = v(x,t)$ 表示趋化物质浓度;$\nu$ 为外单位法向量;$W(t)$ 是定义在可分Hilbert空间 $U$ 上的柱状Wiener过程;噪声项 $\sigma(u)\,dW(t)$ 具有非局部、非线性乘性结构:

$$ (\sigma(u)e_k)(x) = b_k \|u\|_{L^p(\Omega)}^r u(x),\quad k \in \mathbb{N}, $$

且系数序列 $(b_k)_{k\ge1} \in \ell^2(\mathbb{R})$,参数满足

$$ 2 < p < \infty,\quad 1 < r < \frac{p}{p-2}. $$

该模型源于趋化现象的随机建模:确定性Keller-Segel系统(1.1)刻画了细胞在自分泌化学信号驱动下的定向运动,其核心机制是扩散-趋化竞争;当 $d=2$ 时,经典结果表明解的全局存在性由初始质量 $M_0 = \int_\Omega u_0(x)\,dx$ 决定——若 $M_0 < 4\pi$ 则全局存在且有界,若 $M_0 > 4\pi$(且 $M_0 \notin \{4k\pi\mid k\in\mathbb{N}\}$),则解在有限或无限时间内爆破。这一临界质量现象根植于自由能泛函 $\mathcal{F}[u,v] = \int_\Omega u\log u\,dx + \frac{1}{2}\int_\Omega |\nabla v|^2\,dx - \int_\Omega uv\,dx$ 的凹性与梯度流结构,其耗散率 $-\frac{d}{dt}\mathcal{F} = \int_\Omega u|\nabla(\log u - v)|^2\,dx$ 在 $\int u = 4\pi$ 处失去控制能力。

相比之下,随机扰动引入了新的物理与数学维度:环境噪声(如温度、光照波动)作为公共噪声源,在粒子系统平均场极限下不可忽略,导致宏观方程含非退化随机项。已有工作集中于抛物-椭圆情形(即 $-\Delta v + v = u$ 或 $-\Delta v = u$),例如 Chen et al. [7] 与 Tang–Wang [33] 证明:当 $\sigma(u) = b_k \|u\|_{L^p}^r u$ 且 $r < p/(p-2)$ 时,非局部噪声可通过增强耗散抑制爆破,实现任意初值下的全局解;但该结论依赖椭圆变量 $v$ 的瞬时响应特性($v = (-\Delta + 1)^{-1}u$),从而将系统降阶为单一方程并利用 $(-\Delta)^{1/2}$ 型分数阶估计。而本文研究的抛物-抛物耦合(1.3)中,$v$ 满足完整热方程,其时间演化与 $u$ 同步且存在记忆效应,导致:(i)无法消去 $v$ 得到闭合的 $u$-方程;(ii)$\nabla v$ 的正则性需通过联合能量估计获得,而标准 $L^\infty$ 或 $W^{1,q}$ 估计(见引理2.2)在二维中仅给出 $\int_0^t \|\nabla v(s)\|_{L^\infty}^2 ds \lesssim t + t^2 \int_0^t \|u(s)\|_{H^1}^2 ds$,缺乏对 $\|\nabla v\|_{L^\infty}$ 的一致控制;(iii)趋化项 $-\nabla \cdot (u \nabla v)$ 缺乏强制性(coercivity):$\int u \nabla v \cdot \nabla \phi$ 在 $\phi = u$ 时不产生 $\|\nabla u\|_{L^2}^2$ 型耗散,反而可能放大 $\|u\|_{L^2}$;(iv)噪声系数 $\sigma(u)$ 既不满足全局Lipschitz条件(因 $\|\sigma(u)-\sigma(v)\|_{\mathcal{L}_2} \sim \|u\|_{L^p}^r \|u-v\|_{L^2} + \cdots$,参见引理2.1(2.2)),也不满足标准单调性或Hemicontinuity假设,致使经典SPDE存在性理论(如[22, Theorem 4.2.4])失效。因此,本文旨在填补关键空白:在完全抛物耦合、无小初值假设、无强制性与全局Lipschitz性的三重困难下,建立非局部非线性噪声驱动的全局强解与路径唯一性,揭示随机扰动如何实质性地改变二维Keller-Segel系统的临界质量范式。

2. 主要定理或主要结果

在假设 $2 < p < \infty$、$1 < r < p/(p-2)$ 及初值 $(u_0,v_0) \in H^1(\Omega) \times W^{1,\infty}(\Omega)$ 下,本文证明:系统 (1.3) 存在唯一的全局强解 $(u,v)$,满足

$$ u \in L^2_{\mathrm{loc}}([0,\infty); H^1(\Omega)) \cap L^\infty_{\mathrm{loc}}([0,\infty); L^2(\Omega)),\quad v \in L^\infty_{\mathrm{loc}}([0,\infty); W^{1,\infty}(\Omega)), $$

且对任意 $T > 0$,有

$$ \sup_{t\in[0,T]} \|u(t)\|_{L^2}^2 + \int_0^T \|\nabla u(s)\|_{L^2}^2\,ds \le C(T, \|u_0\|_{H^1}, \|v_0\|_{W^{1,\infty}}) < \infty. $$

更精确地,存在常数 $\alpha \in \left(\max\left\{0, \frac{2-r}{2}\right\}, \frac{1}{2}\right)$,使得Lyapunov泛函 $\Phi(x) = (1+x)^\alpha$ 满足:对截断解 $(u_N,v_N)$,有

$$ \mathbb{E}\left[\Phi\left(\|u_N(t)\|_{L^2}^2\right)\right] + \mathbb{E}\int_0^t \left[ c_1 \|\nabla u_N(s)\|_{L^2}^2 + c_2 \|u_N(s)\|_{L^2}^2 \right] \Phi'\left(\|u_N(s)\|_{L^2}^2\right) ds \le \mathbb{E}\left[\Phi\left(\|u_0\|_{L^2}^2\right)\right] + C_\alpha B_0^2 \mathbb{E}\int_0^t \|u_N(s)\|_{L^p}^{2r} \|u_N(s)\|_{L^2}^2 \Phi''\left(\|u_N(s)\|_{L^2}^2\right) ds, $$

其中 $c_1,c_2 > 0$,$B_0^2 = \sum_{k=1}^\infty b_k^2$。利用 $\Phi'' < 0$ 及Gagliardo–Nirenberg插值 $\|u\|_{L^p}^{2r} \le \eta \|\nabla u\|_{L^2}^2 + C_\eta \|u\|_{L^2}^{\beta}$($\beta = 4r/(p-r(p-2))$),结合 $r < p/(p-2)$ 保证 $\beta < 2$,可导出 $\mathbb{E}\left[\Phi(\|u_N(t)\|_{L^2}^2)\right]$ 的一致有界性。进而由Markov不等式,截断停止时间 $\tau_N = \inf\{t > 0 : \|u_N\|_{X_t}^2 \ge N\}$ 满足 $\mathbb{P}(\tau_N \le T) \le \frac{C_T}{N}$,故 $\lim_{N\to\infty} \mathbb{P}(\tau_N > T) = 1$,从而极大局部解可延拓至全局。路径唯一性由差值估计保证:对任意两个局部强解 $(u^{(1)},v^{(1)})$, $(u^{(2)},v^{(2)})$,有

$$ \mathbb{E}\left[\sup_{s\in[0,t\wedge\tau]} \|u^{(1)}(s)-u^{(2)}(s)\|_{L^2}^2\right] \le C \mathbb{E}\int_0^{t\wedge\tau} \left[ \|u^{(1)}-u^{(2)}\|_{L^2}^2 + \|v^{(1)}-v^{(2)}\|_{H^1}^2 \right] ds, $$

结合Gronwall引理及 $v$-方程的线性性,得 $u^{(1)} = u^{(2)}$ a.s.,从而确立路径唯一性。

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

本文的核心技术障碍在于:(1)完全抛物耦合导致 $u$ 与 $v$ 的正则性相互制约,无法分离估计;(2)趋化漂移项 $-\nabla\cdot(u\nabla v)$ 和非局部噪声 $\sigma(u)$ 均缺失强制性与全局Lipschitz性,标准SPDE框架失效;(3)二维临界情形下,$L^p$ 范数对 $p>2$ 的控制不足以直接控制 $\|\nabla v\|_{L^\infty}$,而后者是控制趋化项的关键。

作者的真正创新点体现在三层次分析框架:

第一,构造基于联合能量范数的截断方案。定义 $\|\rho\|_{X_t}^2 := \sup_{s\in[0,t]} \|\rho(s)\|_{L^2}^2 + \int_0^t \|\nabla \rho(s)\|_{L^2}^2 ds$,并设 $\theta_N(r) = \theta(r/N)$($\theta \in C^2$,支集在 $[0,2]$,$\theta\equiv1$ on $[0,1]$)。截断系统为

$$ \begin{cases} du_N = \Delta u_N dt - \theta_N(\|u_N\|_{X_t}^2) \nabla\cdot(u_N \nabla v_N) dt + \theta_N(\|u_N\|_{X_t}^2) \sigma(u_N) dW(t), \\ dv_N = (\Delta v_N - v_N + u_N) dt, \end{cases} $$

此设计的关键在于:当 $\theta_N \neq 0$ 时,自动有 $\|u_N\|_{X_t}^2 < 2N$,从而将非线性项的大小控制在截断水平内,为后续收缩映射提供先验界。值得注意的是,该截断并非作用于单个函数(如 $\|u\|_{L^2}$),而是作用于耗散型范数 $\|\cdot\|_{X_t}$,精准匹配抛物系统的能量结构。

第二,发展“冻结系数+收缩映射”的局部存在性策略。由于截断后系统仍不满足Lipschitz条件,作者在Banach不动点框架中冻结 $v_N$ 和 $\sigma(u_N)$:对给定 $\bar{v} \in C([0,T]; H^1)$ 和 $\bar{u} \in L^2(\Omega_T)$,求解线性化方程 $du = \Delta u dt - \theta_N(\|\bar{u}\|_{X_t}^2) \nabla\cdot(\bar{u} \nabla \bar{v}) dt + \theta_N(\|\bar{u}\|_{X_t}^2) \sigma(\bar{u}) dW(t)$,并证明映射 $\Gamma:(\bar{u},\bar{v}) \mapsto (u,v)$ 在合适加权空间中为严格收缩。此步骤依赖于Burkholder–Davis–Gundy不等式(公式2.1)和引理2.1的局部Lipschitz型估计(2.3),后者通过Gagliardo–Nirenberg插值将 $\|\sigma(u)-\sigma(v)\|_{\mathcal{L}_2}^2$ 控制为 $\varepsilon \|\nabla(u-v)\|_{L^2}^2 + C_{\varepsilon,\eta}(u,v) \|u-v\|_{L^2}^2$,其中 $\varepsilon$ 可任意小,$C_{\varepsilon,\eta}$ 依赖于 $u,v$ 的 $H^1$ 范数,从而允许在小时间区间上吸收耗散项。

第三,设计适配系统结构的凹性Lyapunov泛函。区别于文献[3,7,33]中常用的 $\Phi(x)=x^\alpha$ 或 $\log(1+x)$,本文选取 $\Phi(x)=(1+x)^\alpha$($\alpha \in ((2-r)/2, 1/2)$),其二阶导数 $\Phi''(x) = \alpha(\alpha-1)(1+x)^{\alpha-2} < 0$。应用Itô公式于 $\Phi(\|u_N\|_{L^2}^2)$,关键项为

$$ \mathbb{E}\int_0^t \Phi''(\|u_N\|_{L^2}^2) \cdot 2\langle u_N, \sigma(u_N) dW\rangle^2 \approx -C_\alpha B_0^2 \mathbb{E}\int_0^t \|u_N\|_{L^p}^{2r} \|u_N\|_{L^2}^2 (1+\|u_N\|_{L^2}^2)^{\alpha-2} ds, $$

该负项与趋化项产生的正贡献 $\mathbb{E}\int_0^t \langle u_N, \nabla\cdot(u_N \nabla v_N)\rangle \Phi'(\|u_N\|_{L^2}^2) ds$ 竞争。利用 $\|\nabla v_N\|_{L^\infty} \le C(1+\|u_N\|_{L^p})$(引理2.2(i))及Young不等式,可证此负项主导正项,从而获得 $\mathbb{E}[\Phi(\|u_N\|_{L^2}^2)]$ 的一致上界。此构造的精妙之处在于:$\alpha < 1/2$ 确保 $\Phi''$ 衰减足够快以平衡 $\|u_N\|_{L^p}^{2r}$ 的增长,而 $\alpha > (2-r)/2$ 则保证 $\Phi'$ 的衰减速率足以控制噪声的非线性强度,二者共同由参数约束 $r < p/(p-2)$ 所保障。

4. 与之前工作的比较

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

(1)模型设定与适用范围:此前关于随机Keller-Segel的研究几乎全部集中于抛物-椭圆情形(如[7,16,23,33]),即 $v$ 满足稳态方程 $-\Delta v + v = u$。此类模型可消去 $v$ 得到单一方程 $du = (\Delta u - \nabla\cdot(u\nabla (-\Delta+1)^{-1}u))dt + \sigma(u)dW$,其非线性项具有弱于 $\nabla u$ 的正则性(属 $H^{-1}$ 类),便于利用分数阶算子估计。而本文首次系统处理完全抛物-抛物耦合(1.3),保留了 $v$ 的动力学记忆,使问题升维至强耦合系统,技术难度显著提升。文献[11]虽研究一维抛物-抛物情形,但采用Stratonovich噪声且仅得martingale解(弱解),而本文在二维获得强解与路径唯一性。

(2)初值条件与临界性突破:确定性抛物-抛物系统(1.1)在 $\Omega \subset \mathbb{R}^2$ 上要求 $\int_\Omega u_0 < 4\pi$ 才能保证全局存在(Nagai et al. [26])。本文彻底摆脱此限制,对任意 $u_0 \in H^1$、$v_0 \in W^{1,\infty}$ 均成立,揭示非局部噪声可完全消除确定性临界质量阈值。对比Chen et al. [7](抛物-椭圆情形),其全局性亦无需小初值,但依赖椭圆结构提供的即时正则化;本文则证明该效应在更物理的抛物-抛物框架下依然成立,且噪声参数范围 $r < p/(p-2)$ 与[7]一致,表明该约束是维持耗散优势的本质条件。

(3)技术代价与局限:本文的全局存在性依赖于非局部噪声结构 $\sigma(u) = b_k \|u\|_{L^p}^r u$,其强度随全局 $L^p$ 范数增长。若采用局部噪声(如 $\sigma(u) = u \dot{W}$),则引理2.1的估计失效,$\|\sigma(u)-\sigma(v)\|_{\mathcal{L}_2}^2 \sim \|u-v\|_{L^2}^2$ 不足以控制趋化项,故无法推广。此外,证明要求 $v_0 \in W^{1,\infty}$(而非仅 $H^1$),源于引理2.2(i)中 $\|\nabla v\|_{L^\infty}$ 估计对初值正则性的依赖;若放宽至 $v_0 \in H^1$,则 $\|\nabla v\|_{L^\infty}$ 仅具 $L^2_t L^{2+}_x$ 可积性,不足以控制 $\int u \nabla v \cdot \nabla u$。此为抛物-抛物耦合固有的正则性代价。

(4)方法论推进:相较于Huang–Qiu [16](使用Galerkin逼近与紧性)或Misiats et al. [23](基于熵方法),本文发展出针对强耦合SPDE的新型截断-Lyapunov框架。特别是,联合能量范数 $\|\cdot\|_{X_t}$ 的截断设计、冻结系数的收缩映射构造、以及凹性泛函 $\Phi(x)=(1+x)^\alpha$ 的参数精调,均为处理非强制性非线性SPDE提供了可迁移的范式。从提供内容中无法确认本文结果能否推广至三维或无界区域 $\mathbb{R}^2$,因引理2.2的估计在无界域中失效,且 $r < p/(p-2)$ 在 $d\ge3$ 时变为 $r < p/(p-d)$,可能导致 $\beta \ge 2$,破坏Lyapunov论证。


8. 具有小初值的流体-薄板模型解的指数衰减

Exponential Decay of Solutions to a Fluid-Plate Model with Small Initial Data

作者: Mario Bukal, Igor Kukavica, Linfeng Li, Boris Muha

类别: math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23756v1

摘要翻译: 我们考虑一个三维流体-结构相互作用问题,该问题将不可压 Navier-Stokes 方程耦合于一个随时间变化的区域中,并与定义在流体运动上边界的平方根阻尼薄板方程相耦合。我们在先验意义下证明:当初始数据在某个适当的 Sobolev 空间中充分小时,强解具有指数衰减性。证明结合了高阶能量估计、Stokes 型正则性估计以及一种在小初值假设下闭合的非线性 bootstrap 方法。

Comment: 无

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

1. 研究对象和背景

本文研究三维不可压 Navier-Stokes 流体与移动上边界上四阶平方根阻尼板方程耦合的流体-结构相互作用(FSI)系统。设 $ \mathbb{T}^2 = (0,1)^2 \subset \mathbb{R}^2 $ 为单位环面,流体域随时间演化为 $ \Omega_{\eta}(t) = \{ (x,z) \in \mathbb{R}^3 : x = (x_1,x_2) \in \mathbb{T}^2,\, z \in (0, \eta(x,t)) \} $,其中 $ \eta : \mathbb{T}^2 \times [0,+\infty) \to \mathbb{R} $ 表征上边界高度。流体满足不可压 Navier-Stokes 方程:

$$ \partial_t u - \Delta u + u \cdot \nabla u + \nabla p = 0 \quad \text{in } \Omega_{\eta}(t) \times (0,T), \tag{2.1} $$

$$ \mathrm{div}\, u = 0 \quad \text{in } \Omega_{\eta}(t) \times (0,T). \tag{2.2} $$

上边界 $ \Gamma(t) = \{ (x,\eta(x,t)) : x \in \mathbb{T}^2 \} $ 上的结构运动由四阶平方根阻尼板方程描述:

$$ \eta_{tt} - \Delta_x \eta_t + \Delta_x^2 \eta = g \quad \text{in } \mathbb{T}^2 \times (0,T), \tag{2.3} $$

其中 $ \Delta_x = \partial_{11} + \partial_{22} $ 为水平拉普拉斯,外力密度 $ g $ 由流体应力在法向 $ n_\eta = (1+|\nabla \eta|^2)^{-1/2}(-\partial_1 \eta, -\partial_2 \eta, 1) $ 上的垂直分量给出:

$$ g = - S_\eta (-p I + \nabla u)(x, \eta(x,t), t) \, n_\eta \cdot e_3, \quad S_\eta = \sqrt{1 + |\nabla \eta|^2}. \tag{2.4} $$

边界条件为底部无滑移 $ u(x,0,t) = 0 $,顶部速度匹配 $ u(x,\eta(x,t),t) = (0,0,\eta_t) $。该模型物理上刻画了粘性流体驱动弹性薄板振动的典型场景,如血流-血管壁、微流控芯片中液滴驱动柔性膜等,在生物力学与微机电系统中具有基本建模意义。

从数学角度看,该问题的核心困难在于几何非线性与强耦合:流体域 $ \Omega_{\eta}(t) $ 随 $ \eta $ 变形,导致方程系数含 $ \nabla \eta $ 的高阶非线性;界面力 $ g $ 同时依赖于压力 $ p $ 和速度梯度 $ \nabla u $,而 $ p $ 又由不可压约束隐式决定,缺乏显式控制;此外,板方程 $ \eta_{tt} - \Delta_x \eta_t + \Delta_x^2 \eta $ 兼具双曲($ \eta_{tt} $)、抛物($ -\Delta_x \eta_t $)与四阶椭圆($ \Delta_x^2 \eta $)特性,其耗散机制弱于标准阻尼波方程(如 $ \eta_{tt} + \eta_t + \Delta_x^2 \eta $),仅通过 $ \Delta_x \eta_t $ 提供一阶耗散,对高频模态抑制能力有限。已有工作多聚焦于弱解存在性(如 [CDEG, Gr1, LR, MC1, MS, TW]),或强解局部存在性([Be, Le1, Le2, DT])。在全局强解方面,[IKLT2] 和 [KO1] 对带结构阻尼的流体-弹性体系统证明了小初值下的指数衰减,但要求初始位移 $ \eta_0 \equiv 1 $(即初始域为平坦参考域 $ \Omega = \mathbb{T}^2 \times (0,1) $),从而规避了 Lagrangian 坐标下几何系数 $ J_0 = \det(\nabla \psi_0) = 1 + h_0(x) $ 和变形矩阵 $ a = (\nabla \psi)^{-1} $ 的非平凡初值问题;而本文允许 $ \eta_0(x) = 1 + h_0(x) $ 非恒定,即初始板位移 $ h_0 \not\equiv 0 $,导致 $ J_0 $ 和 $ a $ 在 $ t=0 $ 时已含非零低阶项,破坏了标准能量泛函的正定性与耗散结构,尤其在推导 equipartition estimate(能量均分估计)时,需精细处理 $ J_0 $、$ a $ 及其时间导数对高阶能量演化的扰动。因此,本文试图填补的关键空白是:在初始几何非平坦($ h_0 \neq 0 $)且无额外稳定化项(如 [IKLT2] 中的边界摩擦或 [QGY] 中的边界耗散)的前提下,建立三维流体-四阶平方根阻尼板耦合系统的全局强解及其指数衰减性。

2. 主要定理或主要结果

定理 2.1 给出核心结论:设 $ (v,q,a,\Psi,h) $ 是 Lagrangian 坐标下系统 (2.11)–(2.15) 在时间区间 $ [0,T) $($ 0 < T \leq +\infty $)上的光滑解,定义高阶能量范数

$$ Y(t) := \|v\|_{H^3} + \|v_t\|_{H^1} + \|v_{tt}\|_{L^2} + \|h\|_{H^4} + \|h_t\|_{H^3} + \|h_{tt}\|_{H^2} + \|h_{ttt}\|_{L^2}. \tag{2.16} $$

则存在常数 $ C \geq 1 $ 与 $ \epsilon > 0 $,使得若初值满足 $ Y(0) \leq \epsilon $,则解在整个存在区间内满足指数衰减估计

$$ Y(t) \leq C Y(0) e^{-t/C}, \quad \forall t \in [0,T). \tag{2.17} $$

该结果为先验估计,不涉及解的存在性构造(见 Remark 2.3)。关键技术细节包括:(i) 能量空间为 Sobolev 空间 $ H^k $,而非更弱的 $ L^2 $ 或 Besov 空间;(ii) 衰减率 $ 1/C $ 为绝对常数,不依赖于所预设的最终时间 $ T $(区别于 [Le1, Le2] 中依赖 $ T $ 的小初值条件);(iii) 初值 $ Y(0) $ 包含位移 $ h $、速度 $ h_t $、加速度 $ h_{tt} $ 及其高阶导数,且 $ h $ 的 $ H^4 $ 范数直接控制几何变形程度;(iv) 小初值阈值 $ \epsilon $ 依赖于物理参数(文中归一化为 1,Remark 2.4 指出其一般依赖于流体粘度、板密度、弯曲刚度等);(v) 时间导数 $ v_t, v_{tt}, h_t, h_{tt}, h_{ttt} $ 在 $ t=0 $ 处由方程 (2.11)–(2.13) 迭代求导并代入初值得到,隐含高阶相容性条件(Remark 2.3)。

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

证明的核心障碍有三:(1) 几何非线性主导的系数扰动:由于允许非零初位移 $ h_0 \neq 0 $,Lagrangian 坐标下的 Jacobian $ J_0 = 1 + h_0 $ 和逆变形矩阵 $ a = (\nabla \psi)^{-1} $ 在 $ t=0 $ 时已含 $ \nabla h_0 $ 项,其时间导数 $ a_t $、$ (J_0 a)_t $ 等在能量估计中产生非耗散型高阶耦合项;(2) 压力 $ q $ 的隐式性与低正则性:$ q $ 由不可压约束和界面耦合确定,基本能量 $ \int |v|^2 + |h_t|^2 + |\Delta_x h|^2 $ 无法直接控制 $ \nabla q $,而 $ \nabla q $ 又出现在动量方程的源项中,形成闭合循环;(3) 平方根阻尼的弱耗散:板方程 $ h_{tt} - \Delta_x h_t + \Delta_x^2 h $ 的耗散项 $ -\Delta_x h_t $ 仅提供 $ \|\nabla_x h_t\|_{L^2}^2 $ 级耗散,不足以直接控制 $ \|h_{tt}\|_{H^2} $ 或 $ \|\Delta_x^2 h\|_{L^2} $,需借助流体耗散与界面耦合实现能量转移。

作者的真正创新点在于构建了一个兼具几何鲁棒性与耗散封闭性的非线性 bootstrap 框架:首先,引入修正的 tangential 微分算子集 $ \mathcal{F} = \{ \partial_t, \bar{\partial}, \partial_{tt}, \bar{\partial} \partial_t, \bar{\partial}^2 \} $(其中 $ \bar{\partial} = \partial_1 $ 或 $ \partial_2 $),对系统进行高阶切向微分,避免法向导数带来的边界奇异性;其次,设计两类互补的能量恒等式:(i) 对 $ S \in \mathcal{F} $,将 (2.11) 与 $ Sv_i $ 做内积、(2.13) 与 $ Sh_t $ 做内积,得到恒等式 (4.5),其左边含 $ \frac{d}{dt} \left( \int_\Omega J_0 |Sv|^2 + \int_{\mathbb{T}^2} |Sh_t|^2 + \int_{\mathbb{T}^2} |\Delta_x Sh|^2 \right) $ 与耗散项 $ \int_\Omega |S\nabla v|^2 + \int_{\mathbb{T}^2} |\nabla_x Sh_t|^2 $,右边则包含所有非线性耦合误差;(ii) 对同一 $ S $,将 (2.13) 与 $ Sh $ 做内积得 (4.6),将 (2.11) 与精心构造的测试函数 $ \Phi_i = S\Psi_i - S\Psi_i(x,z,\tau) + z \delta_{i3} Sh(x,\tau) $ 做内积得 (4.10),二者组合消去关键的坏项 $ \int_{\mathbb{T}^2} Sq|_{z=1} Sh $,揭示界面处的 cancellation structure(抵消结构)。第三,发展 Stokes 型正则性估计:利用 Piola 恒等式 $ \partial_i(J_0 a_{ij}) = 0 $ 将 (2.11)–(2.12) 改写为标准 Stokes 形式 (3.19),再应用经典 Stokes 正则性引理 3.2,得到关键估计 (3.18):

$$ \|v_t\|_{H^2} + \|v\|_{H^3} + \|\nabla q\|_{H^1} + \|\nabla q_t\|_{L^2} \lesssim \|v_{tt}\|_{L^2} + \|v_t\|_{L^2} + \|h_t\|_{H^3} + \|h_{tt}\|_{H^2}, $$

该估计将压力高阶导数完全由速度与板变量控制,打破了压力的隐式困境。最后,通过非线性 bootstrap:假设在 $ [0,t] $ 上 $ Y(s) \leq 2\epsilon $,利用 (3.9)–(3.10) 控制几何系数 $ a, J_0 $ 及其导数($ \|a_t\|_{L^\infty H^2} \lesssim \epsilon $),代入 (4.5)、(4.6)、(4.10) 的右边误差项,结合 Poincaré 不等式获得 coercivity($ \int_\Omega |Sv|^2 \lesssim \int_\Omega |S\nabla v|^2 $),最终在 $ \epsilon $ 足够小时吸收所有误差,闭合 bootstrap 循环,得到 (2.17)。

4. 与之前工作的比较

本文与既有工作的主要推进体现在四个维度:

  1. 初始几何假设的实质性放宽:[KO1] 要求 $ \eta_0 \equiv 1 $(即 $ h_0 \equiv 0 $),使 $ J_0 \equiv 1 $、$ a(0) = I $,几何系数扰动可忽略;而本文允许 $ \|h_0\|_{H^4} \leq \gamma $(Lemma 3.1),首次在三维 FSI 中处理非平坦初值下的指数衰减,代价是需建立更精细的几何系数估计 (3.9)–(3.10),并引入 shifted variables $ h = \eta -1 $、$ \Psi = \psi - (x,z) $ 来分离几何扰动。

  2. 耗散机制的纯粹性:[IKLT2] 依赖于边界摩擦项 $ \alpha \partial_t u \cdot \tau $($ \tau $ 为切向量)提供额外耗散;[QGY] 依赖边界耗散;[DT] 要求 Navier 边界条件。本文仅依靠流体粘性 $ -\Delta u $ 和板的平方根阻尼 $ -\Delta_x \eta_t $,证明二者耦合足以实现指数衰减,拓展了纯流体-结构耗散机制的适用范围。

  3. 时间尺度与小初值性质:[Le1, Le2] 的全局强解要求初值 $ \|u_0\|_{H^{5/2}} + \|\eta_0\|_{H^3} + \|\eta_1\|_{H^{5/2}} \leq \epsilon(T) $,其中 $ \epsilon(T) \to 0 $ 当 $ T \to \infty $,即小初值阈值依赖于目标时间;本文的 $ \epsilon $ 为绝对常数,与 $ T $ 无关,属 uniform-in-time smallness,这是更强的结论,代价是需更高阶能量 $ Y(t) $(含 $ v_{tt}, h_{ttt} $ 等)来控制长期行为。

  4. 模型复杂度与物理保真度:[BM] 研究二维简化薄膜模型;[BKMT] 处理无粘流体(无 $ -\Delta u $ 耗散);[Co] 考虑无阻尼线性波方程(无 $ -\Delta_x \eta_t $)。本文采用三维完整 Navier-Stokes 与四阶平方根阻尼板,更贴近实际物理(如血管壁的 Kelvin-Voigt 类型阻尼),但数学难度显著增加。文中明确指出,与 [IKLT1, IKLT2, KO1] 相比,非零初位移导致 equipartition estimate 的证明成为 major difficulty,此点未被先前工作覆盖。从提供内容中无法确认本文结果是否最优(如衰减速率 $ 1/C $ 是否可达 $ 1/2 $),也未显式展开与 [TW](广义 Kirchhoff/von Kármán 板)或 [MS](非线性 Koiter 壳)在正则性提升方面的定量比较。


9. 基于部分切向边界观测的二维Navier-Stokes方程连续数据同化

Continuous Data Assimilation for the 2D Navier-Stokes Equations from Partial Tangential Boundary Observations

作者: Gianmarco Del Sarto, Buddhika Priyasad

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

发布日期: 2026-07-26

arXiv ID: 2607.23707v1

摘要翻译: 我们在光滑、有界、连通区域 $Ω$ 上研究二维Navier–Stokes方程在Navier-滑移边界条件下的连续数据同化问题,且不使用任何内部观测。可用数据仅包含非空相对开子集 $Γ\subset\partialΩ$ 上切向速度的有限维测量。我们证明:若边界反馈足够强、观测足够精细,则同化误差具有强制性谱隙;该谱隙的极限被识别为在 $Γ$ 上施加齐次Dirichlet条件所得混合边界问题的谱隙,并被证明具有量级 $ν$。将此由反馈诱导的强制性与参考解长时间平均对称梯度能量所控制的非线性误差生成估计相结合,我们得到了指数同步化的充分判据。我们在无外力情形下验证了该判据;当未受控Navier-滑移形式具有 $\mathrm{L}^2$ 谱隙时,对足够小的外力成立;在不存在切向刚体运动的区域上,对足够大的粘性系数成立;在存在正边界摩擦的情形下,对足够大的粘性系数亦成立。我们还处理了允许切向刚体运动的区域上的完美滑移情形,此时同步化在强迫项的非刚体不可压分量满足小性条件时成立。

Comment: 无

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

1. 研究对象和背景

本文研究二维有界光滑连通区域 $\Omega \subset \mathbb{R}^2$ 上带Navier-slip边界条件的不可压Navier-Stokes方程的连续数据同化问题,其核心挑战在于:所有观测信息严格局限于边界子集 $\Gamma \subset \partial\Omega$ 上的有限维切向速度测量,且无任何内部观测。该设定源于海洋学、微流控等实际场景——内部传感器部署不可行,而沿岸或壁面切向流速(如ADCP、壁面剪切应力仪)可被高分辨率获取。研究对象是参考解 $u$ 与同化解 $v$ 的误差 $w = v - u$,二者均满足不可压约束 $\nabla \cdot u = 0$、$\nabla \cdot v = 0$ 及法向无穿透条件 $u \cdot n = 0$、$v \cdot n = 0$,但边界动力学存在本质差异。

参考解 $u$ 满足标准Navier-slip系统:

$$ \begin{cases} \partial_t u - \nu \Delta u + (u \cdot \nabla) u + \nabla p = g, & \text{in } \Omega, \\ \nabla \cdot u = 0, & \text{in } \Omega, \\ u D(u)n \cdot \tau + \alpha u_\tau = 0, & \text{on } \partial\Omega, \end{cases} $$

其中 $D(u) = \frac{1}{2}(\nabla u + \nabla u^T)$ 为对称梯度,$\tau = n^\perp$ 为单位切向量,$u_\tau = u \cdot \tau$,$\alpha \geq 0$ 为边界摩擦系数;当 $\alpha = 0$ 时退化为完美滑移(perfect slip),此时 $\nu D(u)n \cdot \tau = 0$,允许刚性切向运动(tangential rigid motions)存在,即非零解 $r \in V$ 满足 $D(r) = 0$,例如圆盘上的旋转场 $r(x) = (-x_2, x_1)$。该自由度在经典粘性耗散中不可阻尼,构成同步障碍。

同化解 $v$ 由边界反馈系统驱动:

$$ \begin{cases} \partial_t v - \nu \Delta v + (v \cdot \nabla) v + \nabla q = g, & \text{in } \Omega, \\ \nabla \cdot v = 0, & \text{in } \Omega, \\ \nu D(v)n \cdot \tau + \alpha v_\tau + \mu P_\delta^* (P_\delta T_\Gamma v - y_\delta) = 0, & \text{on } \Gamma, \\ \nu D(v)n \cdot \tau + \alpha v_\tau = 0, & \text{on } \partial\Omega \setminus \Gamma, \end{cases} $$

其中 $y_\delta(t) = P_\delta T_\Gamma u(t)$ 是 $u_\tau|_\Gamma$ 的有限维观测,$P_\delta : L^2(\Gamma) \to Y_\delta$ 为有界线性算子(如截断傅里叶系数或边界单元平均),满足假设1.1(稠密观测):$\|P_\delta f\|_{Y_\delta} \leq C_P \|f\|_{L^2(\Gamma)}$ 且 $\|P_\delta f\|_{Y_\delta} \to \|f\|_{L^2(\Gamma)}$ 当 $\delta \to 0$;$\mu > 0$ 为反馈强度;$P_\delta^*$ 将观测残差 $P_\delta T_\Gamma v - y_\delta \in Y_\delta$ 提升为作用于 $\Gamma$ 的边界力。关键在于,该反馈不修改内部方程,仅通过修正Navier-slip律中的切向应力项实现闭环,与传统体积分布式nudging(如AOT算法)有本质区别。

误差 $w = v - u$ 满足齐次方程:

$$ \partial_t w - \nu \Delta w + (v \cdot \nabla) w + (w \cdot \nabla) u + \nabla \pi = 0, $$

并继承混合边界条件:在 $\Gamma$ 上 $\nu D(w)n \cdot \tau + \alpha w_\tau + \mu P_\delta^* P_\delta T_\Gamma w = 0$,在 $\partial\Omega \setminus \Gamma$ 上 $\nu D(w)n \cdot \tau + \alpha w_\tau = 0$。其能量平衡(Lemma 2.5)为

$$ \frac{1}{2} \frac{d}{dt} \|w(t)\|_{L^2(\Omega)}^2 + a_{\mu,\delta}(w(t), w(t)) = - \int_\Omega (w(t) \cdot \nabla) u(t) \cdot w(t) \, dx, $$

其中双线性型 $a_{\mu,\delta}(z,z) = 2\nu \|D(z)\|_{L^2(\Omega)}^2 + \alpha \|z_\tau\|_{L^2(\partial\Omega)}^2 + \mu \|P_\delta T_\Gamma z\|_{Y_\delta}^2$ 集成三重耗散:粘性、边界摩擦、反馈。问题的核心空白在于:如何从纯边界切向观测中提取足够强的$L^2$-尺度 coercivity,以压制非线性项 $-(w \cdot \nabla)u \cdot w$ 的增长? 现有理论(如Foias-Prodi确定模态、AOT算法)依赖内部观测或全局边界信息,而边界反馈稳定化工作(如[3,4,20])针对预设平衡态设计控制器,不处理未知参考轨迹的实时同化。本文首次严格建立:仅凭 $\Gamma$ 上有限维 $u_\tau$ 观测,能否实现指数同步?其根本障碍是Navier-slip算子在 $L^2$ 中缺乏固有谱隙(尤其当 $\alpha = 0$ 且存在刚性运动时),而边界反馈必须克服此缺陷,生成一个“人工谱隙” $\kappa > 0$ 满足 $a_{\mu,\delta}(z,z) \geq \kappa \|z\|_{L^2(\Omega)}^2$。该谱隙的大小、可达性及与物理参数($\nu, \alpha, \Gamma$)的定量关系,构成理论基石。

2. 主要定理或主要结果

本文核心结论是建立了在纯边界切向观测下二维Navier-Stokes方程连续数据同化的指数同步充分条件,其技术细节高度精确:

定理1.2(反馈诱导的coercivity):设 $\Gamma \subset \partial\Omega$ 非空且相对开,Assumption 1.1成立。则存在仅依赖于 $\Omega$ 和 $\Gamma$ 的正常数 $c_\Gamma, C_\Gamma > 0$,使得极限混合边值问题的谱隙满足

$$ c_\Gamma \nu \leq \kappa_{\infty,\Gamma} := \inf_{\substack{z \in V \setminus \{0\} \\ T_\Gamma z = 0}} \frac{a_\alpha(z,z)}{\|z\|_{L^2(\Omega)}^2} \leq C_\Gamma \nu, \quad \forall \nu > 0, \, \alpha \geq 0. $$

进一步,对任意 $0 < \kappa < \kappa_{\infty,\Gamma}$,存在 $\mu_0 > 0$ 和 $\delta_0 > 0$,使得当 $\mu \geq \mu_0$ 且 $0 < \delta \leq \delta_0$ 时,

$$ a_{\mu,\delta}(z,z) \geq \kappa \|z\|_{L^2(\Omega)}^2, \quad \forall z \in V. $$

该定理表明:反馈强度 $\mu$ 与观测分辨率 $\delta$ 的联合调控可逼近极限谱隙 $\kappa_{\infty,\Gamma}$,且其量级为 $O(\nu)$,与粘性系数线性相关。

定理1.3(基于平均对称梯度控制的指数同步):设 $u$ 为(1.1)的解,$M^D_u := \limsup_{t \to \infty} \frac{1}{t} \int_0^t \|D(u(s))\|_{L^2(\Omega)}^2 \, ds < \infty$,$\Theta_{\kappa,\nu} := C_4^2 C_{\text{Korn}} (\kappa^{-1} + (2\nu)^{-1})$($C_4$ 为Ladyzhenskaya常数,$C_{\text{Korn}}$ 为Korn不等式常数)。若存在 $\kappa$ 满足

$$ 0 < \kappa < \kappa_{\infty,\Gamma} \quad \text{且} \quad \Theta_{\kappa,\nu} M^D_u < \kappa, $$

则存在 $\mu_0 > 0$, $\delta_0 > 0$,使得对任意 $\mu \geq \mu_0$, $0 < \delta \leq \delta_0$, 以及任意 $0 < \gamma < \frac{1}{2}(\kappa - \Theta_{\kappa,\nu} M^D_u)$,存在常数 $C > 0$(独立于 $t, v_0, \mu, \delta$)满足

$$ \|v(t) - u(t)\|_{L^2(\Omega)} \leq C e^{-\gamma t} \|v_0 - u_0\|_{L^2(\Omega)}, \quad \forall t \geq 0. $$

该结论将同步性归结为谱隙 $\kappa$ 与参考流平均形变能 $M^D_u$ 的竞争。

具体可验证准则(Propositions 4.3–4.9): (i) 无外力情形(Corollary 1.4):若 $g = 0$,则 $M^D_u = 0$,故对任意 $\nu > 0$ 和 $\kappa \in (0, \kappa_{\infty,\Gamma})$,存在 $\mu_0, \delta_0$ 使同步成立。 (ii) 大粘性无刚性运动:若唯一满足 $\nabla \cdot r = 0$ 且 $r \cdot n = 0$ 的刚性场为零,且 $\nu^4 > \frac{C_4^2 C_\Omega^3}{8} \|g\|_{L^2(\Omega)}^2$($C_\Omega$ 为Korn-Poincaré常数),则同步成立。 (iii) 大粘性正边界摩擦:若 $\alpha > 0$,且 $\nu > \max\left\{ \frac{C_4^2 C_{\alpha,\nu}}{\lambda_{\alpha,\nu}^2} \|g\|_{L^2(\Omega)}^2, \, \nu \right\}$($\lambda_{\alpha,\nu}, C_{\alpha,\nu}$ 来自引理A.5),则同步成立;此条件适用于圆盘、同心环等含刚性运动的域。 (iv) 小外力与未受控谱隙:若存在 $\lambda_\alpha > 0$ 使得 $a_\alpha(z,z) \geq \lambda_\alpha \|z\|_{L^2(\Omega)}^2$(即未加反馈时已有$L^2$谱隙),且 $\|g\|_{L^2(\Omega)}^2 < \frac{2\nu \lambda_\alpha \kappa}{\Theta_{\kappa,\nu}}$,则同步成立。 (v) 完美滑移含刚性运动:若 $\alpha = 0$ 且 $R_\Omega = \{r \in V : D(r) = 0\} \neq \{0\}$,令 $g^\perp = (I - \Pi_R) P_H g$($\Pi_R$ 为到 $R_\Omega$ 的正交投影),若 $\frac{C_R \Theta_{\kappa,\nu}}{4\nu^2} \|g^\perp\|_{L^2(\Omega)}^2 < \kappa$($C_R$ 为模刚性运动的Korn-Poincaré常数),则同步成立。

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

证明的核心障碍在于:Navier-slip算子 $a_\alpha(z,z)$ 在 $L^2(\Omega)$ 中通常不coercive,尤其当 $\alpha = 0$ 且存在非零刚性场 $r \in R_\Omega$ 时,$a_\alpha(r,r) = 0$,导致谱隙 $\kappa_{\infty,\Gamma}$ 可能退化。本文的真正创新点在于构建并量化了边界反馈生成的“人工谱隙”机制,并将其与参考流的形变能精确耦合,而非依赖传统内插型观测的全局控制。

1. 三重谱隙分析框架(Section 3):作者引入三个递进的Rayleigh商:实际反馈谱隙 $\kappa_{\mu,\delta} = \inf_{z \in V \setminus \{0\}} \frac{a_{\mu,\delta}(z,z)}{\|z\|_{L^2}^2}$、理想全迹谱隙 $\kappa_{\mu,\Gamma} = \inf_{z \in V \setminus \{0\}} \frac{a_{\mu,\Gamma}(z,z)}{\|z\|_{L^2}^2}$($a_{\mu,\Gamma}(z,z) = a_\alpha(z,z) + \mu \|T_\Gamma z\|_{L^2(\Gamma)}^2$)、极限混合边值谱隙 $\kappa_{\infty,\Gamma} = \inf_{\substack{z \in V \setminus \{0\} \\ T_\Gamma z = 0}} \frac{a_\alpha(z,z)}{\|z\|_{L^2}^2}$。关键引理3.1证明 $\kappa_{\mu,\Gamma} \nearrow \kappa_{\infty,\Gamma}$ 当 $\mu \to \infty$,即强反馈趋近于强制 $z_\tau|_\Gamma = 0$ 的混合问题;Proposition 3.4则利用Assumption 1.1(稠密观测)证明 $\kappa_{\mu,\delta} \to \kappa_{\mu,\Gamma}$ 当 $\delta \to 0$,从而确立 $\kappa_{\mu,\delta} \geq \kappa$ 对足够大 $\mu$ 和足够小 $\delta$ 成立。此框架将离散观测的逼近性与连续谱隙的渐近行为严格链接。

2. 刚性运动的定向阻尼技术(Section 4.4):当 $\alpha = 0$ 且 $R_\Omega \neq \{0\}$ 时,经典粘性无法耗散刚性分量。作者的关键洞察是:边界反馈 $\mu P_\delta^* P_\delta T_\Gamma w$ 虽不直接作用于 $R_\Omega$ 内部,但通过分解 $w = w^\perp + w^R$($w^R \in R_\Omega$),利用 $T_\Gamma w^R$ 的非零性(因刚性场在 $\Gamma$ 上一般非零)和 $P_\delta$ 的稠密性,证明 $\mu \|P_\delta T_\Gamma w\|_{Y_\delta}^2$ 对 $w^R$ 产生有效阻尼。引理A.3(模刚性运动的Korn-Poincaré不等式)提供 $\|w^\perp\|_{L^2} \leq C_R \|D(w^\perp)\|_{L^2}$,而 $\|w^R\|_{L^2}^2 \lesssim \|T_\Gamma w^R\|_{L^2(\Gamma)}^2 \lesssim \|P_\delta T_\Gamma w\|_{Y_\delta}^2 + \|P_\delta T_\Gamma w^\perp\|_{Y_\delta}^2$,最终导出 $a_{\mu,\delta}(w,w) \geq c_\nu \|w^\perp\|_{L^2}^2 + c_\mu \|w^R\|_{L^2}^2$,彻底解决刚性运动不可观问题。

3. 非线性项的最优控制:误差方程中的非线性项 $-\int_\Omega (w \cdot \nabla) u \cdot w \, dx$ 被精细估计为 $\leq C \|D(u)\|_{L^2} \|w\|_{L^2} \|w\|_{H^1}$。结合Korn不等式 $\|w\|_{H^1} \leq C_K \|D(w)\|_{L^2}$ 和Young不等式,得到

$$ -\int_\Omega (w \cdot \nabla) u \cdot w \, dx \leq \frac{1}{2} a_{\mu,\delta}(w,w) + \Theta_{\kappa,\nu} \|D(u)\|_{L^2}^2 \|w\|_{L^2}^2, $$

其中 $\Theta_{\kappa,\nu}$ 显式依赖于 $\kappa$ 和 $\nu$。此估计避免了对 $u$ 的高正则性假设,仅需 $M^D_u < \infty$,并通过长时平均Gronwall引理(Appendix A)完成最终同步论证。

4. 与之前工作的比较

本文与既有工作的根本区别在于观测模式、反馈机制与数学结构的三重革新,而非对现有框架的参数优化。

1. 与AOT类内插同化算法的对比:Azouani-Olson-Titi(AOT)算法[1]及其拓展[6,9,18]依赖内部观测(如Fourier模、节点值、局部平均),其插值算子 $I_h$ 作用于 $\Omega$ 内部,构造的反馈项 $\mu (I_h u - I_h v)$ 是体积分佈力。本文则完全摒弃内部观测,仅用 $\Gamma \subset \partial\Omega$ 上的 $P_\delta T_\Gamma u$,反馈项 $\mu P_\delta^* (P_\delta T_\Gamma v - y_\delta)$ 是边界应力修正,直接嵌入Navier-slip律。这导致分析对象从椭圆/抛物型内插误差转变为边界混合问题的谱隙,技术工具从Sobolev嵌入转向边界迹理论与混合边值问题谱分析。文中明确指出:“no previous rigorous result establishes exponential synchronisation […] under this precise combination of boundary-only observations, finite-dimensional tangential measurements, localisation to a boundary subset, and feedback acting directly through the boundary condition”(第3页)。

2. 与边界反馈稳定化工作的对比:文献[3,4,20,21]研究的是状态反馈镇定(state-feedback stabilization),即针对预设平衡态 $\bar{u}$ 设计控制器,目标是 $v \to \bar{u}$,其可观测性不等式是控制设计工具,不对应真实观测流。本文则是轨迹同化(trajectory assimilation),$u$ 是未知的、由外力 $g$ 驱动的演化解,$y_\delta$ 是其真实切向迹,目标是 $v \to u$。这要求反馈必须适应 $u$ 的时变特性,且非线性项 $-(w \cdot \nabla)u \cdot w$ 的处理远比稳定化中围绕平衡态的线性化复杂。此外,[30]虽考虑边界滑移,但其nudging仍为体积分佈,与本文边界律修正有本质不同。

3. 与边界观测反问题的对比:文献[2]基于速度与Cauchy应力的边界迹做稳定性估计,属适定性反问题;文献[31]的Navier-Stokes观测器使用内部观测。本文不重构初值或参数,而是实时同步轨迹,且观测仅为标量切向速度 $u_\tau$,不含应力信息,难度更高。

4. 推进与代价:本文将同步适用范围从内部推广至边界,时间尺度上实现全局指数同步($t \geq 0$),而不仅限于大时间行为。代价在于:(i) 需要足够强的反馈 $\mu$ 和足够细的观测 $\delta$($\mu_0, \delta_0$ 依赖于 $\Omega, \Gamma, \nu$);(ii) 同步速率 $\gamma$ 受限于 $\kappa - \Theta_{\kappa,\nu} M^D_u$,当 $M^D_u$ 很大时需更大 $\kappa$,即更强反馈;(iii) 完美滑移情形需额外小性条件 $\|g^\perp\|$,因刚性运动分量 $g^R$ 无法被反馈抑制。文中未显式展开与[27](海洋学边界观测应用)的数值比较,亦未讨论噪声观测(如[5,8]),这些是未来方向。


10. 局部有限图上带幂-对数非线性的热方程正解的寿命

The lifespan of positive solutions of heat equation with power-logarithmic nonlinearity on locally finite graph

作者: Pengxiu Yu, Yiping Zhang

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23251v1

摘要翻译: 在局部有限连通图 $G=(V,E)$ 上,利用 Kaplan \cite{MR160044} 引入的首特征值方法以及 Hu-Wang \cite{cvhuyuanyang} 发展的离散 Phragmén-Lindelöf 原理,我们首先建立了带幂-对数非线性项 $u^p|\log u|^q$ 的半线性热方程正解寿命的渐近行为,前提是初始数据有正的下界。这些结果将 Hu-Wang \cite{cvhuyuanyang} 的结论推广至具有幂-对数源项的方程。此外,通过一种更为直接的论证,我们证明:当非负初始数据 $u(x,0)$ 在某个顶点 $x_i\in V$ 处满足 $u(x_i,0)$ 充分大时,相应的寿命估计依然成立。

Comment: 无

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

1. 研究对象和背景

本文研究定义在无限局部有限连通加权图 $ G = (V, E, w, \mu) $ 上的半线性热方程

$$ u_t - \Delta u = u^p |\log u|^q, \quad (x,t) \in V \times (0,T), \tag{1.3} $$

其中图拉普拉斯算子定义为

$$ \Delta u(x) = \frac{1}{\mu(x)} \sum_{y \sim x} w_{xy} \big( u(y) - u(x) \big), $$

$ \mu: V \to \mathbb{R}_+ $ 为顶点测度,$ w_{xy} = w_{yx} > 0 $ 为边权,且满足 $ D_\mu := \sup_{x \in V} \frac{1}{\mu(x)} \sum_{y \sim x} w_{xy} < \infty $(即图具有有界度量几何)。初始数据取为 $ u(x,0) = \lambda \psi(x) $,其中 $ \lambda > 0 $ 为尺度参数,$ \psi: V \to [0,\infty) $ 为有界非负函数,不恒为零。研究对象是该初值问题的正古典解 $ u \in C^1_V((0,T)) \cap C^0_V([0,T)) $,其生命期(lifespan)定义为最大存在时间

$$ T_\lambda := \sup \big\{ T > 0 : \text{方程 } (1.3) \text{ 在 } [0,T) \text{ 上存在唯一非负古典解} \big\}. $$

若 $ T_\lambda < \infty $,则称解在有限时间内爆破(blow-up),且由比较原理可知爆破必为 $ L^\infty(V) $-范数意义下的全局爆破。本文聚焦于两类临界非线性结构:(i) $ p = 1, q > 1 $,对应对数超线性源项 $ u |\log u|^q $;(ii) $ p > 1, q \ge 1 $,对应幂-对数耦合源项 $ u^p |\log u|^q $。这两类非线性处于经典 Fujita 型阈值 $ p = 1 + \frac{2}{N} $ 的离散类比之外——因图无固有维度,传统欧氏尺度分析失效;而 $ p=1 $ 时纯线性热方程 $ u_t = \Delta u $ 解全局存在,故引入对数扰动旨在刻画线性与幂型非线性之间的精细过渡行为。物理上,此类非线性出现在网络扩散模型中受对数增长抑制或增强的种群动力学、信息传播饱和机制等场景;数学上,其凸性结构 $ f(u) = u^p |\log u|^q $ 在 $ u \ge e^{1-q} $(当 $ p=1 $)或 $ u \ge 1 $(当 $ p>1 $)时成立,直接决定 Jensen 不等式与能量泛函单调性的适用范围。已有工作如 Hu-Wang (2024) 在相同图框架下处理了纯幂型 $ u^p $ 情形,证明了

$$ \lim_{\lambda \to \infty} \lambda^{p-1} T_\lambda = \frac{1}{(p-1) \|\psi\|_{L^\infty(V)}^{p-1}}, $$

其证明依赖于第一特征值方法与离散 Phragmén–Lindelöf 原理。然而,当源项含对数因子时,标准特征值方法失效:因 $ f'(u) = u^{p-1}|\log u|^q + q u^{p-1} |\log u|^{q-1} \operatorname{sgn}(\log u) $ 在 $ u \to 1 $ 处不可控,且 $ f $ 的增长速率不再由单一幂次主导,导致测试函数法、Duhamel 表示与热核估计难以直接移植。本文试图填补的关键空白在于:在无欧氏尺度、无体积增长假设的抽象图结构上,建立幂-对数非线性热方程生命期的精确渐近公式,并阐明对数指数 $ q $ 如何通过 $ (\log \lambda)^{q-1} $ 修正经典幂律标度,同时将结果从全局正初值推广至仅在单点充分大的非负初值——这触及离散框架下“局域爆破触发全局爆破”的机制本质。

2. 主要定理或主要结果

本文在 $ D_\mu < \infty $ 的局部有限连通图上,对初值 $ u(x,0) = \lambda \psi(x) $,建立了四组主定理:

Theorem 1.1($ p = 1, q > 1 $,全局正初值):若 $ \eta(0) := \int_\Omega \lambda \psi \phi \, d\mu \ge 1 $(其中 $ \Omega \subset V $ 有限连通,$ \phi > 0 $ 为 $ -\Delta_\Omega $ 对应第一特征值 $ \lambda_1(\Omega) $ 的归一化特征函数),且 $ \Lambda \inf_{x\in V} \psi(x) \ge \max\{ C(q,V), e^{1-q} \} $,则存在 $ \Lambda > 0 $ 使得当 $ \lambda > \Lambda $ 时 $ T_\lambda < \infty $,且

$$ \lim_{\lambda \to +\infty} \big( \log (\lambda \|\psi\|_{L^\infty(V)}) \big)^{q-1} T_\lambda = \frac{1}{q-1}. $$

Theorem 1.2($ p = 1, q > 1 $,局域大初值):若 $ \deg(x_i) = m \ge 1 $ 且 $ \lambda \psi(x_i) > C(V) \exp\{ m^{1/q} \} $,则解 $ u(x_i,t) $ 在有限时间内爆破,且

$$ \limsup_{\lambda \to +\infty} \big( \log(\lambda \psi(x_i)) \big)^{q-1} T_\lambda \le \frac{1}{q-1} \le \liminf_{\lambda \to +\infty} \big( \log(\lambda \|\psi\|_{L^\infty(V)}) \big)^{q-1} T_\lambda. $$

Theorem 1.4($ p > 1, q \ge 1 $,局域大初值):若 $ \deg(x_i) = m \ge 1 $ 且 $ (\lambda \psi(x_i))^{p-1} (\log (\lambda \psi(x_i)))^q > C(V) m $ 以及 $ \lambda \psi(x_i) \ge 1 $,则 $ T_\lambda < \infty $,且

$$ \limsup_{\lambda \to +\infty} (\lambda \psi(x_i))^{p-1} (\log(\lambda \psi(x_i)))^q T_\lambda \le \frac{1}{p-1} \le \liminf_{\lambda \to +\infty} (\lambda \|\psi\|_{L^\infty(V)})^{p-1} (\log(\lambda \|\psi\|_{L^\infty(V)}))^q T_\lambda. $$

Theorem 1.5($ p > 1, q \ge 1 $,全局正初值):若 $ \eta(0) \ge 1 $ 且 $ \Lambda \inf_{x\in V} \psi(x) \ge \max\{ C(p,q,V), 1 \} $,则当 $ \lambda > \Lambda $ 时 $ T_\lambda < \infty $,且

$$ \lim_{\lambda \to +\infty} (\lambda \|\psi\|_{L^\infty(V)})^{p-1} (\log(\lambda \|\psi\|_{L^\infty(V)}))^q T_\lambda = \frac{1}{p-1}. $$

所有常数 $ C(\cdot) $ 仅依赖于所列参数及图的几何量 $ D_\mu $,不显含维度或体积增长阶。

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

证明的核心障碍在于:传统基于热核或尺度变换的方法在离散无尺度图上失效,而幂-对数非线性破坏了纯幂型问题中 $ f(u)/u^p $ 的齐次性,导致第一特征值方法无法直接应用。作者的真正创新点在于构建了一套双层比较框架,将连续 ODE 技术与离散谱理论深度耦合。

第一层:加权限制能量泛函的构造与凸性驱动的微分不等式。针对 $ p=1 $ 情形,定义加权平均 $ \eta(t) = \int_\Omega u(\cdot,t) \phi \, d\mu $,其中 $ \phi $ 是有限子图 $ \Omega $ 上拉普拉斯的第一特征函数。关键突破在于发现:当 $ u \ge e^{1-q} $ 时,函数 $ g(u) = u |\log u|^q $ 在 $ (e^{1-q},\infty) $ 上严格凸,从而可对 $ \eta(t) $ 应用 Jensen 不等式得到

$$ \eta'(t) \ge -\lambda_1 \eta(t) + \eta(t) |\log \eta(t)|^q. $$

此微分不等式不依赖图的全局结构,仅需 $ D_\mu < \infty $ 保证 $ \Delta $ 的有界性(引理 2.4),且其右端项精确捕获了对数非线性的主导增长。求解该 ODE 得到爆破时间上界 $ T_\lambda \le \frac{1}{q-1} (\log \eta(0))^{1-q} $,进而导出定理 1.1 的上极限估计。

第二层:显式上下解的构造与离散 Phragmén–Lindelöf 原理的精巧运用。为获得下极限,作者构造了仅依赖 $ \lambda \|\psi\|_{L^\infty} $ 的显式上解

$$ \bar{u}(t) = \exp\left\{ \big[ (\log(\lambda \|\psi\|_{L^\infty}))^{1-q} - (q-1)t \big]^{-1/(q-1)} \right\}, $$

验证其满足 $ \bar{u}_t - \Delta \bar{u} \ge \bar{u} |\log \bar{u}|^q $(利用 $ D_\mu < \infty $ 控制离散拉普拉斯作用于常值函数的误差),再通过引理 2.3(离散 Phragmén–Lindelöf)确保 $ u \le \bar{u} $,从而得到 $ T_\lambda \ge \frac{1}{q-1} (\log(\lambda \|\psi\|_{L^\infty}))^{1-q} $。此构造避免了连续情形所需的热核逐点估计,纯粹基于图的度量有界性。

第三层:局域爆破的“单点触发”机制分析。对于定理 1.2 与 1.4,作者放弃全局加权平均,转而利用顶点度 $ \deg(x_i) = m $ 构造以 $ x_i $ 为中心的星形子图 $ \Omega_i = \{x_i\} \cup \{ y \sim x_i \} $,并证明:若 $ u(x_i,0) $ 足够大,则 $ u(x_i,t) $ 满足 ODE 不等式 $ u_t(x_i,t) \ge \frac{1}{\mu(x_i)} \sum_{y\sim x_i} w_{x_i y} (u(y,t) - u(x_i,t)) + u^p |\log u|^q $。通过忽略邻点贡献(仅保留源项主导)并利用 $ \deg(x_i) $ 提供的离散梯度上界,导出 $ u_t(x_i,t) \ge C m^{-1} u^p |\log u|^q $,再积分得爆破时间上界。此技术揭示了图结构中度数作为局部几何约束如何量化“大初值”的阈值,是纯连续分析中不存在的新现象。

4. 与之前工作的比较

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

1. 非线性结构的实质性推进:Hu-Wang (2024) 处理 $ u^p $ 时,其生命期标度 $ \lambda^{1-p} $ 由幂次唯一决定;本文首次确立 $ u |\log u|^q $ 的标度为 $ (\log \lambda)^{1-q} $,证明对数指数 $ q $ 主导渐近行为——这超越了 Chabi-Souplet (2019) 在欧氏空间对 $ u^p L(u) $($ L $ 缓变)的定性描述,给出了离散框架下首个精确对数标度公式。特别地,当 $ q \to 1^+ $,$ (\log \lambda)^{q-1} \to 1 $,恢复 Hu-Wang 的幂律,验证了连续性。

2. 几何假设的极大放宽:Euclidean 空间中 Lee-Ni (1980) 与 Gui-Wang (2005) 的精确常数结果要求 $ \psi(x) \to \psi_\infty > 0 $(无穷远处收敛),而本文仅需 $ \psi $ 有界且不恒为零,且完全不依赖图的体积增长、曲率或热核衰减——仅靠 $ D_\mu < \infty $ 这一最弱测度有界性即可,适用于任意快速/慢速生长图(如树、Cayley 图、随机图)。

3. 初值条件的降维突破:以往所有图上结果(如 [10,21,22])均要求 $ \psi \ge c > 0 $ 全局成立;本文定理 1.2 与 1.4 首次证明:只要存在单个顶点 $ x_i $ 满足 $ \lambda \psi(x_i) $ 超过由度数 $ \deg(x_i) $ 决定的阈值 $ \exp\{ \deg(x_i)^{1/q} \} $,即可触发全局爆破。这揭示了离散爆破的“薄弱环节”本质,是连续框架中不可能的现象(因欧氏空间无孤立高权重点)。

4. 方法论代价与局限:精确常数 $ \frac{1}{q-1} $ 的获得以牺牲初值正则性为代价——要求 $ \psi $ 有界但无需连续或衰减;然而,本文无法处理 $ \psi $ 在无穷远处趋于零的情形(如 $ \psi(x) \sim \operatorname{dist}(x,x_0)^{-\alpha} $),因 Jensen 不等式失效且局域方法无法控制远距离耦合。此外,所有结果均未涉及爆破解的形态(如是否集中爆破),亦未讨论 $ q \le 1 $ 或 $ p < 1 $ 的亚线性情形,这些是未来工作方向。


11. 关于非齐次 Dirichlet 边界数据下不规则区域上 Navier-Stokes 方程的极弱解适定性理论

On well-posedness theory of very weak solutions to Navier-Stokes equations on irregular domains with nonhomogeneous Dirichlet boundary data

作者: Xiaojin Bai, Siran Li, Xiangxiang Su

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23234v1

摘要翻译: 极弱解的适定性理论是数学流体力学中的核心课题,尤其在 Navier-Stokes 方程的正则性理论中占据关键地位。该理论已在 $\mathbb{R}^3$ 中具有 $C^{2,1}$-正则性的有界区域上对不可压流体流动得到充分发展。本文基于 [D. Breit 和 A. Gaudin, ArXiv 预印本: 2511.19091 (2025)] 与 [V.G. Maz’ya 和 T.O. Shaposhnikova, Vol.337, Grundlehren der mathematischen Wissenschaften (2009)] 中的解析理论,建立了 Navier-Stokes 方程在有界 Lipschitz 区域上的极弱解适定性理论;这些区域的边界具有局部图表示函数,且其 Sobolev 乘子范数足够小,特别地,包含 Lipschitz 常数足够小的有界 Lipschitz 区域作为特例。

Comment: 22 pages

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

1. 研究对象和背景

本文研究三维有界Lipschitz区域 $\Omega \subset \mathbb{R}^3$ 上非定常不可压缩Navier-Stokes方程(NSE)的极弱解(very weak solutions)的适定性问题,其具体形式为:

$$ \begin{cases} \partial_t u - \Delta u + \operatorname{div}(u \otimes u) + \nabla \pi = \operatorname{div} F & \text{in } (0, T) \times \Omega, \\ \operatorname{div} u = 0 & \text{in } (0, T) \times \Omega, \\ u = g & \text{on } (0, T) \times \partial \Omega, \\ u(0) = u_0 & \text{in } \Omega. \end{cases} \tag{1.1} $$

其中外力项 $F \in L^s(0,T; L^r(\Omega))$,边界数据 $g \in L^s(0,T; W^{1-1/q,q}(\partial \Omega))$,初值 $u_0 \in L^q_\sigma(\Omega)$,且满足相容性条件

$$ \int_{\partial \Omega} g(t) \cdot N \, dS = 0 \quad \text{a.e. } t \in (0,T), \tag{1.3} $$

其中 $N$ 是 $\partial \Omega$ 上的单位外法向量场。参数范围被严格限定为:$3 < q < \infty$, $2 < s < \infty$, $1 < r < q$,并满足尺度关系

$$ \frac{1}{3} + \frac{1}{q} = \frac{1}{r}, \quad \frac{2}{s} + \frac{3}{q} = 1. \tag{1.2'} $$

该尺度关系正是经典Ladyzhenskaya–Prodi–Serrin正则性准则中 $L^s_t L^q_x$ 类解的临界指数对,亦是极弱解理论中自然出现的时空可积性阈值。

极弱解的核心特征在于:解 $u$ 仅要求属于 $L^s(0,T; L^q(\Omega))$,而不要求任何空间导数的存在性;相应地,测试函数空间需大幅收缩以保证所有配对有意义。在经典 $C^{2,1}$ 区域上,Farwig–Galdi–Sohr 等人已建立完备理论:利用 $W^{2,q}$-正则性,可将边界项 $\langle g, N \cdot \nabla w \rangle_{\partial \Omega}$ 显式定义,并通过分部积分将方程转化为对 $w \in C^1_c([0,T); C^2_{0,\sigma}(\Omega))$ 成立的弱形式。然而,在Lipschitz区域上,$\partial \Omega$ 仅为局部Lipschitz连续,导致边界迹算子 $\operatorname{Tr}: W^{1,q}(\Omega) \to W^{1-1/q,q}(\partial \Omega)$ 虽仍有定义,但 $\nabla w$ 在边界上的点态意义丧失,传统分部积分公式失效;更严重的是,Stokes算子 $A_q$ 的定义域 $D(A_q)$ 不再等于 $W^{2,q}(\Omega) \cap W^{1,q}_0(\Omega) \cap L^q_\sigma(\Omega)$,标准Calderón–Zygmund理论完全崩溃。Shen (2007) 的里程碑工作表明:在 $d$ 维Lipschitz区域上,Stokes半群在 $L^q_\sigma(\Omega)$ 上解析当且仅当 $|1/q - 1/2| < 1/(2d) + \varepsilon$,即对 $d=3$,仅允许 $q \in (3/2 - \varepsilon, 3 + \varepsilon)$,而极弱解所需的 $q > 3$ 正好落在该区间的右侧补集内——二者几乎不交,构成根本性障碍。因此,如何在 $q > 3$ 的范围内重建Stokes算子的全谱性质与半群估计,是发展极弱解理论的关键瓶颈。本文旨在填补这一空白:在比经典Lipschitz区域更精细的几何分类下,建立 $q > 3$ 时极弱解的局部存在唯一性与小数据全局存在唯一性。该分类并非基于经典Hölder或Sobolev范数,而是基于边界图函数的Sobolev乘子范数,其物理含义是:局部扰动对Sobolev空间乘法作用的‘微小性’,从而在不提升整体光滑度的前提下,保障了关键线性估计的成立。

2. 主要定理或主要结果

定理1.1(主定理):存在普适常数 $\varepsilon_0 > 0$,使得以下结论成立:设 $\Omega \subset \mathbb{R}^3$ 为有界Lipschitz区域,其边界 $\partial \Omega$ 属于Sobolev乘子类 $M^{1+\alpha,\rho}_W(\varepsilon)$(定义见2.2),其中 $\alpha \in (0,1), \rho \in [1,\infty], \varepsilon < \varepsilon_0$。假设数据 $F, g, u_0$ 满足正则性假设 (1.2) 及相容性条件 (1.3),且 $q,s,r$ 满足 (1.2’)。则存在依赖于数据的 $T^* = T^*(F,g,u_0) \in (0,T]$,使得方程 (1.1) 存在唯一的极弱解 $u \in L^s(0,T^*; L^q(\Omega))$,并满足先验估计

$$ \|u\|_{L^s(0,T; L^q(\Omega))} \leq C \left( \|F\|_{L^s(0,T; L^r(\Omega))} + \|g\|_{L^s(0,T; W^{1-1/q,q}(\partial \Omega))} + \|u_0\|_{L^q(\Omega)} \right), \tag{1.4} $$

其中常数 $C = C(\Omega, \alpha, \rho, q)$ 与 $T$ 无关。进一步,存在仅依赖于 $\Omega, \alpha, \rho, q$ 的正数 $\mu > 0$,使得若数据满足小数据条件

$$ \|F\|_{L^s(0,T; L^r(\Omega))} + \|g\|_{L^s(0,T; W^{1-1/q,q}(\partial \Omega))} + \|u_0\|_{L^q(\Omega)} \leq \mu, \quad \forall T > 0, \tag{1.5} $$

则存在唯一的全局极弱解 $u \in L^s(0,\infty; L^q(\Omega))$,且满足 (1.4)。注1.2指出,该结果可推广至非齐次Besov空间框架:将 $L^q(\Omega)$ 替换为 $B^s_{p,q}(\Omega)$,其中 $1<p<\infty, 1\leq q \leq \infty, s \in (-1+1/p, 1/p)$,并相应调整数据空间,结论依然成立。核心技术支撑来自命题2.3:对任意 $p \in \mathcal{E}_q := \{q, q', r, r'\}$,在 $\varepsilon < \varepsilon_0$ 的乘子类假设下,Stokes算子 $A_p$ 满足:(i) Helmholtz投影 $P_p: L^p(\Omega) \to L^p_\sigma(\Omega)$ 有界;(ii) $A_p^{1/2}$ 生成等价范数 $\|A_p^{1/2} v\|_{L^p} \simeq \|\nabla v\|_{L^p}$,且 $D(A_p^{1/2}) = W^{1,p}_{0,\sigma}(\Omega)$;(iii) 线性初值问题 $\partial_t z + A_p z = h, z(0)=0$ 在 $W^{1,\tau}(0,T; L^p_\sigma(\Omega)) \cap L^\tau(0,T; D(A_p))$ 中有唯一解,并满足最大正则性估计 $\|\partial_t z\|_{L^\tau L^p_\sigma} + \|A_p z\|_{L^\tau L^p_\sigma} \leq C \|h\|_{L^\tau L^p_\sigma}$;(iv) Stokes半群 $e^{-tA_q}$ 满足衰减估计 $\|\nabla e^{-tA_{q'}} \psi\|_{L^{(q/2)'}(\Omega)} \leq C t^{-1/2 - 3/(2q)} e^{-\omega t} \|\psi\|_{L^{q'}(\Omega)}$,$t>0$。

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

证明的核心障碍在于:在 $q > 3$ 的极弱解框架下,无法使用传统的 $W^{2,q}$-Stokes正则性来处理边界数据 $g$ 和散度型外力 $\operatorname{div} F$,因为它们会直接进入解的二阶导数项。作者的根本创新在于彻底放弃对Stokes系统本身进行高正则性求解的思路,转而构建两个关键的‘静态辅助场’(stationary auxiliary fields),将非齐次项从原方程中剥离,从而将问题还原为一个齐次初值的、纯线性化的极弱问题。第一个辅助场 $v_F \in W^{1,r}_{0,\sigma}(\Omega)$ 由命题3.1定义,它满足变分恒等式

$$ \int_\Omega \nabla v_F : \nabla \phi \, dx = - \int_\Omega F : \nabla \phi \, dx, \quad \forall \phi \in W^{1,r'}_{0,\sigma}(\Omega), \tag{3.1} $$

其本质是将 $\operatorname{div} F$ 项通过分部积分‘吸收’进 $v_F$ 的梯度中,避免了在原始方程中出现 $\operatorname{div} F$ 的显式作用。由于 $r < q$ 且满足 Sobolev嵌入 $W^{1,r}(\Omega) \hookrightarrow L^q(\Omega)$,$v_F$ 自然属于目标空间 $L^q_\sigma(\Omega)$。第二个辅助场 $U_g \in W^{1,q}(\Omega)$ 由命题3.2构造,它是一个满足 $\operatorname{div} U_g = 0$、$\operatorname{Tr} U_g = g$ 且在 $W^{1,q'}_{0,\sigma}(\Omega)$ 上具有零能量的调和延伸,即

$$ \int_\Omega \nabla U_g : \nabla \phi \, dx = 0, \quad \forall \phi \in W^{1,q'}_{0,\sigma}(\Omega). \tag{3.6} $$

这一性质至关重要,因为它确保了在后续的极弱公式中,$U_g$ 与测试函数的耦合项能精确抵消边界贡献,从而绕过对 $\nabla \phi$ 边界值的依赖。在此基础上,作者重新定义了极弱解(定义3.3):对任意 $\Phi \in C^1_c([0,T); D(A_{q'}))$,要求

$$ -\int_0^T \int_\Omega u \cdot \partial_t \Phi \, dxdt + \int_0^T \int_\Omega (P_q u - P_q U_g) \cdot A_{q'} \Phi \, dxdt - \int_0^T \int_\Omega (u \otimes u) : \nabla \Phi \, dxdt \\ = \int_\Omega u_0 \cdot \Phi(0) \, dx + \int_0^T \int_\Omega v_F \cdot A_{q'} \Phi \, dxdt. \tag{3.12} $$

此定义的精妙之处在于:所有涉及 $A_{q'} \Phi$ 的项都通过Stokes算子的对偶性得以定义,而 $\nabla \Phi$ 的出现则由估计 (3.14) 控制:$\|\nabla \Phi\|_{L^{(q/2)'}(\Omega)} \leq C \|A_{q'} \Phi\|_{L^{q'}(\Omega)}$,这依赖于命题2.3(iv)给出的半群核估计。最终,通过将 $u = w + P_q U_g + v_F$ 代入,原问题被等价转化为关于 $w$ 的积分方程 $w(t) = e^{-tA_q} (u_0 - P_q U_g(0) - v_F) + \int_0^t e^{-(t-\tau)A_q} \mathcal{P}_q \operatorname{div}( (w + P_q U_g + v_F) \otimes (w + P_q U_g + v_F) ) \, d\tau$,其中 $\mathcal{P}_q$ 是Helmholtz投影。该方程在Banach空间 $X_T := \{ w \in L^s(0,T; L^q(\Omega)) : \|w\|_{X_T} < \infty \}$ 上通过压缩映射原理求解,其不动点存在性完全依赖于命题2.3(iii)(iv)提供的线性估计与非线性项的双线性结构。

4. 与之前工作的比较

本文与既有工作的主要区别体现在三个维度:几何假设、线性理论基础与解的概念框架。首先,在几何假设上,本文显著弱于 $C^{2,1}$ 区域(Farwig–Galdi–Sohr [7,10]),但强于一般Lipschitz区域。经典Lipschitz区域仅要求边界图函数 $\phi_j$ 属于 $\operatorname{Lip}(\mathbb{R}^2)$,而本文要求 $\phi_j \in M^{1+\alpha,\rho}_W(\mathbb{R}^2, \varepsilon)$,即其梯度 $\nabla \phi_j$ 作为乘子作用于 $W^{\alpha,\rho}(\mathbb{R}^2)$ 时的算子范数小于 $\varepsilon$。这一条件蕴含了 $\phi_j \in C^{0,\alpha}(\mathbb{R}^2)$(由Maz’ya–Shaposhnikova [22]),故该类区域严格介于Lipschitz与 $C^{1,\alpha}$ 之间,且当 $\varepsilon \to 0$ 时退化为 $C^{1,\alpha}$ 区域。其次,在线性理论层面,本文依赖Breit–Gaudin [5] 和 Maz’ya–Shaposhnikova [22] 的突破性工作,而非Shen [23] 的 $L^q$-谱理论。Shen的结果限定了 $q$ 的狭窄区间 $(3/2-\varepsilon, 3+\varepsilon)$,而本文通过乘子类假设,首次在 $q > 3$ 的整个Ladyzhenskaya–Prodi–Serrin临界线上建立了Stokes算子的完整最大正则性(命题2.3(iii))与半群核估计(命题2.3(iv))。这是对Tolksdorf [24] 工作的重要推进,后者虽在Lipschitz区域上建立了 $L^q$-NSE理论,但其 $q$ 范围受限于Shen的谱条件。最后,在解的概念上,本文彻底重构了极弱解的定义。以往工作(如[8])依赖于 $\langle g, N \cdot \nabla w \rangle_{\partial \Omega}$ 这一边界配对,这在Lipschitz区域上无定义;而本文通过辅助场 $U_g$ 和 $v_F$,将边界与体积分力完全内化到解的结构中,使定义 (3.12) 仅依赖于 $L^q$-空间内的对偶配对 $\langle \cdot, A_{q'} \Phi \rangle_\Omega$,从而实现了概念上的自洽性。代价在于:该理论目前仅适用于小乘子范数 $\varepsilon < \varepsilon_0$ 的区域,且全局存在性依赖于小数据条件 (1.5),而 $C^{2,1}$ 区域上的理论允许大初值(只要 $L^s_t L^q_x$ 范数有限)。此外,从提供内容中无法确认本文是否改进了时间存在性 $T^*$ 的下界估计(如是否达到 $T^* \gtrsim \|u_0\|_{L^q}^{-s/(s-2)}$),也未显式展开与Coscia [6] 关于稳态NSE在Lipschitz区域上工作的比较。


12. Dacorogna-Moser 型传输映射在 $\mathbb{R}^d$ 上的构造及其在耦合空间测地线中的应用

A Dacorogna-Moser construction of transport maps on $\mathbb{R}^d$ with application to geodesics on the space of couplings

作者: Louis-Pierre Chaintron, Matteo Picco

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23241v1

摘要翻译: Dacorogna 与 Moser 的开创性工作提出了一种从有界区域上一个概率分布到另一个概率分布构造正则传输映射的方法。本文将该构造推广至整个 $\mathbb{R}^d$,适用于严格渐近对数凹测度——这是一类广泛的分布,包含对数凹测度的 Lipschitz 扰动。随后我们利用该构造研究定义在乘积集上的概率测度空间中具有给定边缘律(即耦合)的测地线,并推导其最优性条件,从而回答了 Conforti、Lacker 与 Pal 最近一项工作中提出的公开问题。受 Brenier 关于不可压流体的变分模型及其正则化的启发,我们进一步引入测地线问题的熵正则化形式,该形式可视为耦合空间上的 Schrödinger 桥问题,并同样推导其最优性条件。最后,我们研究当正则化参数趋于零时极小元的收敛性,并证明相关 Lagrange 乘子的收敛性。我们的方法涉及在 $\mathbb{R}^d$ 上对椭圆型与抛物型方程建立关于时间一致的全局正则性估计,其关键在于利用概率论中的反射耦合(reflection coupling)概念来挖掘渐近对数凹测度的结构。

Comment: 无

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

1. 研究对象和背景

本文研究的核心对象是定义在全空间 $\mathbb{R}^d$ 上的概率测度之间的正则传输映射,以及由这些映射诱导的、受边际约束的耦合空间 $\Pi(\mu, \nu) := \{ \rho \in \mathcal{P}_2(\mathbb{R}^d \times \mathbb{R}^d) : (\pi_1)_\# \rho = \mu,\, (\pi_2)_\# \rho = \nu \}$ 上的测地线问题。研究动机源于两个相互交织的方向:一方面,Dacorogna-Moser 构造是构建正则传输映射的经典 PDE 方法,其原始形式 [DM90] 要求定义域 $\Omega \subset \mathbb{R}^d$ 有界且具有光滑边界,并通过求解带 Neumann 边界条件的 Poisson 方程

$$ abla \cdot v = \varphi, \quad \frac{\partial u}{\partial \nu} = 0 \text{ on } \partial \Omega $$

来获得速度场 $v = -\nabla u$,进而由流 $\partial_s \rho_s + \nabla \cdot (\rho_s v) = 0$ 生成传输映射 $T$ 满足 $T_\# \mu = \nu$。该方法的关键在于其对密度扰动 $\nu = (1+\varphi)\mu$ 的显式、可控制的响应性,即 $\|T - \mathrm{Id}\|_{W^{n,\infty}} \lesssim \|\varphi\|_{W^{n,\infty}}$。然而,[DM90] 中所有估计均显式依赖于 $|\Omega|$,导致其无法直接延拓至无界域 $\mathbb{R}^d$。另一方面,Conforti-Lacker-Pal [CLP25] 提出的耦合空间测地线问题是一个新兴的、具有深刻几何与概率意义的变分问题:给定固定边际 $\mu, \nu \in \mathcal{P}_2(\mathbb{R}^d)$,定义子流形 $\Pi(\mu, \nu)$ 上的形式 Riemann 度量为

$$ d_{\Pi(\mu,\nu)}^2(\rho_0, \rho_1) = \inf \left\{ \int_0^1 \frac{1}{2} \|v_t\|_{L^2(\rho_t)}^2 \, dt : \partial_t \rho_t + \nabla \cdot (\rho_t v_t) = 0,\, \rho_t \in \Pi(\mu, \nu),\, (\rho_0, \rho_1) \text{ fixed} \right\}. $$

这一问题旨在刻画 $\Pi(\mu, \nu)$ 的内在几何结构,其困难根源在于该子流形既非凸亦非测地凸 [CLP25, Prop. 2.2],致使标准梯度流理论 [AGS08] 失效;同时,[CLP25, Open Problem 6] 明确提出:如何为该测地线问题建立严格的、基于偏微分方程的最优性条件?此前 [CL25] 仅在强对数凹(strongly log-concave)且 $\log \mu, \log \nu$ 具有有界 Hessian 的强假设下证明了测地线的存在性,但未给出其 PDE 描述。因此,本文试图填补两大空白:其一,将 Dacorogna-Moser 构造从有界域推广至 $\mathbb{R}^d$,为此需克服无界域上椭圆/抛物方程全局正则性估计的根本障碍;其二,在此新构造基础上,为 $\Pi(\mu, \nu)$ 上的测地线问题提供首个严格、全局、基于 Lagrange 乘子的 PDE 刻画,从而回答 [CLP25] 的核心开放问题。关键的数学挑战在于,$\mathbb{R}^d$ 上的全局正则性要求必须由测度 $\mu$ 的尾部衰减性质来保证。作者引入“严格渐近对数凹性”(strictly asymptotically log-concave)作为核心假设:设 $\mu(dx) = e^{-U(x)} dx$,则 $-U$ 需满足 $\liminf_{r \to +\infty} \kappa_{-U}(r) > 0$,其中弱半凸性剖面 $\kappa_f(r) := r^{-2} \inf_{|x-y|=r} \langle \nabla f(x) - \nabla f(y), x-y \rangle$。这一假设远弱于经典对数凹性(即 $\nabla^2 U \geq \alpha I$),它允许 $U$ 在无穷远处渐近平坦,例如 $U(x) = |x|^p$ ($p>1$) 或 $U(x) = |x|^2 + \sin(|x|)$,从而涵盖了 Lipschitz 扰动的对数凹测度。在此框架下,传输映射的构造被归约为求解加权 Laplace 方程

$$ \Delta_\mu \psi := \nabla \cdot (\nabla \psi) + \nabla U \cdot \nabla \psi = \varphi, \quad \int_{\mathbb{R}^d} \psi \, d\mu = 0, $$

其中 $\Delta_\mu$ 是关于 $\mu$ 的 Ornstein-Uhlenbeck 型算子。而测地线问题的变分分析则要求对由该方程解 $\psi$ 所生成的流进行精细的、一致时间的正则性控制,这正是本文技术的核心。

2. 主要定理或主要结果

本文的主要结果分为两个层次,均建立在严格渐近对数凹性假设 (A.1) 之上:首先,定理 1.2 给出了 $\mathbb{R}^d$ 上 Dacorogna-Moser 构造的严格推广。假设 $\mu$ 满足 (A.1),即 $\mu$ 是严格渐近对数凹的,且 $\nabla^j \log \mu$ 对所有 $1 \leq j \leq n+1$ 是 Lipschitz 连续的 ($n \geq 1$);再设扰动 $\varphi \in W^{n,\infty}(\mathbb{R}^d)$ 满足 $\int_{\mathbb{R}^d} \varphi \, d\mu = 0$、$\|\varphi\|_{W^{n,\infty}} \leq M$ 且 $\varphi \geq -1 + \delta$ ($M, \delta > 0$),则存在一个可逆映射 $T: \mathbb{R}^d \to \mathbb{R}^d$,使得 $(1+\varphi) \cdot \mu = T_\# \mu$,且 $T - \mathrm{Id}, T^{-1} - \mathrm{Id} \in W^{n,\infty}(\mathbb{R}^d; \mathbb{R}^d)$,并满足全局范数估计

$$ \|T - \mathrm{Id}\|_{W^{n,\infty}} + \|T^{-1} - \mathrm{Id}\|_{W^{n,\infty}} \leq C \|\varphi\|_{W^{n,\infty}}, $$

其中常数 $C > 0$ 仅依赖于 $\log \mu$ 的界、弱半凸性剖面 $\kappa_{-\log \mu}$、阶数 $n$ 以及 $M, \delta$。该结论的证明基石是命题 2.1,它给出了加权算子 $\frac{1}{2} \Delta_\mu$ 对应半群 $S_t$ 的尖锐估计:对任意 $t > 0$,有 $\|\nabla S_t \varphi\|_{W^{n,\infty}} \leq c e^{-\alpha t} \|\varphi\|_{W^{n,\infty}}$(当 $t \geq 1$)和 $\|\nabla S_t \varphi\|_{W^{n,\infty}} \leq c t^{-1/2} \|\varphi\|_{W^{n,\infty}}$(当 $t \in (0,1]$),其中 $c, \alpha > 0$ 仅依赖于 $\mu$ 和 $n$。由此可得椭圆方程 $\frac{1}{2} \Delta_\mu \psi = \varphi$ 的唯一解 $\psi = -\int_0^{+\infty} S_t \varphi \, dt$ 满足 $\|\nabla \psi\|_{W^{n,\infty}} \leq C \|\varphi\|_{W^{n,\infty}}$。其次,定理 2.6 和定理 2.9 解决了耦合空间测地线问题。在 $\mu, \nu$ 均满足 (A.1) 的前提下,测地线问题(即 $\varepsilon = 0$ 的 (15) 式)存在至少一个最小化器;进一步,若初始/终态 $\rho_0, \rho_1$ 关于 $\mu \otimes \nu$ 的相对熵有限,则对任意 $\varepsilon \in (0,1]$,正则化问题 (15) 是有限的,且存在唯一最小化器,并满足 $\Gamma$-收敛性:

$$ \Gamma\text{-}\lim_{\varepsilon \to 0} A_\varepsilon + i_{01} = A_0 + i_{01}, $$

其中 $i_{01}$ 是指示函数,收敛在 $C([0,1], \Pi(\mu, \nu))$ 的一致收敛拓扑及点态弱收敛拓扑下成立。最关键的是,定理 2.9 为测地线建立了严格的最优性条件:存在唯一的 Lagrange 乘子 $\alpha \in (W^{n,\infty}(\mathbb{R}^d))'$ 和 $\beta \in (W^{n,\infty}(\mathbb{R}^d))'$(定义 2.8),使得对任意测试函数 $\zeta \in C_c^\infty((0,1) \times \mathbb{R}^d \times \mathbb{R}^d)$,有

$$ \int_0^1 \int_{\mathbb{R}^d \times \mathbb{R}^d} \left( \partial_t \zeta + \nabla_x \zeta \cdot \nabla_x \phi + \nabla_y \zeta \cdot \nabla_y \phi \right) d\rho_t dt = 0, $$

且 $\alpha, \beta$ 与速度势 $\phi$ 满足关系

$$ \alpha(x) + \beta(y) = \partial_t \phi(t,x,y) + \frac{1}{2} \left( |\nabla_x \phi|^2 + |\nabla_y \phi|^2 \right). $$

定理 2.11 进一步在更高正则性假设下,赋予了动量守恒方程 $\rho \nabla (\alpha \oplus \beta) = \partial_t (\rho \nabla \phi) + \nabla \cdot (\rho (\nabla \phi \otimes \nabla \phi))$ 以严格的分布意义。

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

本文的证明框架围绕两个核心障碍展开:其一,如何在无界域 $\mathbb{R}^d$ 上对加权椭圆算子 $\Delta_\mu = \Delta + \nabla U \cdot \nabla$ 的解建立全局、一致的高阶 Sobolev 估计;其二,如何在缺乏经典凸性的情况下,为一个高度非线性的、受双重边际约束的变分问题建立严格的 Lagrange 乘子理论。作者的真正创新点在于将概率论中的“反射耦合”(reflection coupling)这一强大工具,系统性地嵌入到 PDE 正则性理论中,从而打通了从测度的几何结构到方程解的定量控制的路径。具体而言,第一个障碍的突破依赖于命题 2.1 中的半群估计。传统方法(如 Schauder 理论或 Calderón-Zygmund 理论)在 $\mathbb{R}^d$ 上失效,因为它们无法处理 $\nabla U$ 的无界增长。作者转而利用 $\Delta_\mu$ 生成的扩散过程 $dX_t = -\nabla U(X_t) dt + \sqrt{2} dB_t$ 的 Feynman-Kac 表示,并构造一对由反射耦合驱动的进程 $(X_t, Y_t)$。关键洞察在于,严格渐近对数凹性 $\liminf_{r \to \infty} \kappa_{-U}(r) > 0$ 保证了该耦合具有指数收缩率,即存在 $\alpha > 0$ 使得 $\mathbb{E}|X_t - Y_t| \leq e^{-\alpha t} \mathbb{E}|X_0 - Y_0|$。这一概率性质被转化为对半群 $S_t$ 的导数的尖锐控制:通过 [PW06, EZ19] 的经典结果,$\|\nabla S_t \varphi\|_{L^\infty}$ 的衰减速率直接由耦合的收缩率决定,从而得到 $t \geq 1$ 时的指数衰减和 $t \leq 1$ 时的 $t^{-1/2}$ 奇异性。这种将 PDE 估计还原为随机过程行为的方法,是本文最根本的技术创新。第二个障碍的解决则依赖于定理 1.2 提供的、可精确控制的传输映射。在 Brenier 的变分框架中,为证明最优性,需要构造一族满足边际约束的扰动曲线 $\rho_t^s$。在有界域上,[DM90] 的构造能保证 $\rho_t^s$ 的正则性;而在 $\mathbb{R}^d$ 上,任何粗糙的扰动都会破坏 $\Pi(\mu, \nu)$ 的约束。作者的精妙之处在于,他们利用定理 1.2 构造的映射 $T_s$ 来生成扰动:设 $\rho_t$ 是候选测地线,则定义 $\rho_t^s = (\mathrm{Id}, T_s)_\# \rho_t$,其中 $T_s$ 是从 $\mu$ 到 $(1+s\varphi)\mu$ 的映射。由于 $T_s$ 是全局 $W^{n,\infty}$ 的且其范数线性依赖于 $s$,该扰动自动保持在 $\Pi(\mu, \nu)$ 内,并且其变分导数可以被精确计算。这使得作者能够绕过传统方法中对测地凸性的依赖,直接对 Lagrange 泛函进行变分,从而导出定理 2.9 中的乘子方程。此外,对于正则化问题,作者摒弃了文献中常用的热流正则化(因其不保持 $\Pi(\mu, \nu)$),而是创新性地采用相对熵 $H(\cdot | \mu \otimes \nu)$ 在 Wasserstein 度量下的梯度流作为正则化方案,这保证了恢复序列的可行性,并借助 [MTV23] 发展的广义 Schrödinger 桥理论获得了所需的收敛性。

4. 与之前工作的比较

本文与既有工作在多个维度上实现了实质性推进。首先,与 Dacorogna-Moser 原始工作 [DM90] 及其后续改进 [RY96] 相比,本文首次将该构造成功延拓至全空间 $\mathbb{R}^d$。[DM90] 的估计明确依赖于区域测度 $|\Omega|$,其核心是 Neumann 边界条件带来的紧性;而本文则完全放弃了边界,转而利用测度 $\mu$ 自身的渐近几何(即严格渐近对数凹性)来获得全局正则性。代价是正则性损失:[DM90] 给出 $C^{k+1,\alpha}$ 映射,而本文仅得到 $W^{n,\infty}$ 映射,但这是全局性所必需的交换。其次,与基于 Langevin 动力学的构造 [KM12, FMS24] 相比,本文方法提供了更清晰的扰动响应性。[KM12] 构造的是从 $\mu$ 到 $\nu = e^{-V}\mu$ 的 Lipschitz 映射,但其 Lipschitz 常数依赖于 $V$ 的凸性参数,而非一个直接的 $\|V\|_{\mathrm{Lip}}$ 控制;而本文的定理 1.2 显式给出了 $\|T - \mathrm{Id}\|_{W^{n,\infty}} \lesssim \|\varphi\|_{W^{n,\infty}}$,这使其在测地线变分中成为无可替代的工具。再次,与 Gaussian 设置下的工作 [MS23, CCE25, Sté26] 相比,本文的设定更为一般:这些工作均假设参考测度 $\mu$ 是标准 Gaussian,而本文允许 $\mu$ 是任意严格渐近对数凹测度,极大地拓宽了适用范围。最后,与测地线问题的先驱工作 [CLP25, CL25] 相比,本文的贡献是奠基性的。[CL25] 仅在强对数凹且 Hessian 有界的强假设下证明了存在性,而本文将其推广至更弱的严格渐近对数凹性;更重要的是,[CLP25] 仅提出了问题,而本文首次给出了完整的、严格的 PDE 最优性条件(定理 2.9),这在概念上将一个纯概率/几何问题转化为了一个可分析的 PDE 问题。然而,本文也存在明确的局限:它未解决 Lagrange 乘子 $\alpha, \beta$ 的最优正则性问题(例如,它们是否属于某个 Sobolev 空间?),而 Brenier 的不可压缩流工作 [FA07] 已对此进行了深入探讨;此外,本文的正则性假设 (A.1) 要求 $\nabla^{n+1} \log \mu$ Lipschitz,这比 [FMS24] 中仅需 $\nabla^3 U$ 有界的条件更强,反映了 Dacorogna-Moser 方法对参考测度光滑性的固有要求。


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

1. 最优传输中的全局 $W^{2,p}$ 正则性

Global $W^{2,p}$ Regularity in Optimal Transport

作者: Shibing Chen, Yuanyuan Li, Xianduo Wang

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24666v1

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

摘要翻译

中文标题: 最优传输中的全局 $W^{2,p}$ 正则性

摘要翻译: 本文建立了有界凸区域之间二次最优传输的凸势函数的全局 $W^{2,p}$ 估计。源密度与目标密度仅被假设为连续且在零与无穷大之间有界;对区域边界的正则性不作任何假设。


2. Sharp kinetic trace theory

Sharp kinetic trace theory

作者: Lukas Niebel, Lisa Valentini

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24708v1

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

摘要翻译

中文标题: Sharp kinetic trace theory

摘要翻译: 我们建立了多个速度模型下的尖锐动能迹估计及其反例。对于半空间位置区域,在速度支撑无任何公共界的情形下,我们证明了关于Lebesgue测度与标准Gaussian测度的自然迹估计。密度性导出相应动能能量空间上的自然迹算子与Green公式。对于 $d\ge2$ 中的有界空间区域,在有界支撑Euclidean速度模型与球面速度模型中,我们确定了迹权 $\min\{|v \cdot n|,|v \cdot n|^p\}$(其中 $1\le p<\infty$)的尖锐边界正则性阈值。记 $α_p=1/(p+1)$,该估计在每个有界 $\mathrm{C}^{1,α}$ 区域上成立当且仅当 $α\geα_p$;而对每个 $0<α<α_p$,存在一个正则性恰为 $\mathrm{C}^{1,α}$ 的严格凸有界区域,使得该估计失效。特别地,自然迹($p=1$)的正则性阈值为 $\mathrm{C}^{1,1/2}$。在 $d\ge2$ 中的有界 $\mathrm{C}^{1,1/2}$ 区域上,密度性在有界支撑Euclidean速度模型(配以Lebesgue测度或标准Gaussian测度)及球面速度模型中均导出自然迹算子与Green公式。保持范数的速度平移排除了在任意有界 $\mathrm{C}^1$ 区域上无限制的Lebesgue迹估计。在无限制Gaussian模型中,对每个 $1\le p<2$,我们在任意 $d\ge2$ 维有界 $\mathrm{C}^{1,1}$ 区域上构造反例,从而否定地回答了Albritton、Armstrong、Mourrat与Novack(2024)中问题1.8。对 $2\le p<\infty$,Gaussian $ω_2$ 估计与密度性共同导出 $ω_p$-迹算子。


3. 非线性算子第一Robin特征值的渐近行为

Asymptotic behavior of the first Robin eigenvalue of nonlinear operators

作者: Rosa Barbato, Francesco Della Pietra, Alba Lia Masiello

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

发布日期: 2026-07-27

arXiv ID: 2607.24488v1

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

摘要翻译

中文标题: 非线性算子第一Robin特征值的渐近行为

摘要翻译: 设 $Ω$ 是 $\mathbb R^N$ 中的一个有界Lipschitz区域,其中 $N\geq 2$。本文研究 $p$-Laplace 算子的第一 Robin 特征值在 $β\to 0$ 和 $β\to +\infty$ 时的渐近行为,并推导该特征值的尖锐渐近展开式;在 Dirichlet 极限 $β\to+\infty$ 下的展开式是在附加假设 $\partialΩ$ 属于 $C^{1,1}$ 类的条件下得到的。


4. Alexandrov-Bakelman-Pucci 估计的尖锐性与稳定性:一种凸几何方法

Sharpness and Stability of the Alexandrov-Bakelman-Pucci Estimate: a Convex Geometric Approach

作者: Antonio Porricelli

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24388v1

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

摘要翻译

中文标题: Alexandrov-Bakelman-Pucci 估计的尖锐性与稳定性:一种凸几何方法

摘要翻译: 本文基于 G. Talenti 于 1981 年的工作,该工作建立了 Monge-Ampère 方程光滑解的比较原理以及相关的二维先验估计;进一步探讨了其与 Blaschke-Santaló 不等式及其 K. Mahler(1939)提出的反向形式之间的联系。更具体地,本文证明了 Talenti 所得先验估计之一的一个定量版本。我们通过凸几何,挖掘 Talenti 的分析框架与 Mahler/Blaschke-Santaló 不等式之间的相互作用,并展示这一联系如何被用于获得椭圆型偏微分方程 Alexandrov-Bakelman-Pucci (ABP) 极大值原理的尖锐刻画。


5. Vlasov–Poisson 方程组在 $L^{d+}$ 初值密度与速度变量分数阶正则性下的局部适定性

Local Well-Posedness for Vlasov–Poisson with $L^{d+}$ Initial Density and Fractional Velocity Regularity

作者: Quoc-Hung Nguyen

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24400v1

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

摘要翻译

中文标题: Vlasov–Poisson 方程组在 $L^{d+}$ 初值密度与速度变量分数阶正则性下的局部适定性

摘要翻译: 我们在 $\R^d_x\times\R^d_v$ 上、对 $d\geq2$ 的 Vlasov–Poisson 系统,证明了一个局部适定性判据,该判据基于初值分布的各向异性假设。初值具有有限质量,其加权速度上确界属于某个 $p>d$ 的 $L^p_x$ 空间,且在速度变量上具有一致于速度的、任意小的正 Hölder 正则性,并保持相同的关于空间变量的 $L^p$ 控制。主要估计是一个非线性混合型估计

$$ [ρ(t)]_{C^α_x}\lesssim t^{-d/p-ε}C(f_0), \qquad ε>0~~\text{small}. $$

因此密度关于时间可积,取值于一个正的空间 Hölder 类,相应电场属于 $L^1_tC^{1,α}_x$。我们通过在密度上应用 Schauder 不动点定理构造解,并利用 Loeper 型稳定性估计证明唯一性。


6. 一个带有阶数接近零的非局部算子的逻辑斯蒂方程

A logistic equation with non-local operators of order near zero

作者: Luigi Appolloni, Serena Dipierro, Enrico Valdinoci

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24357v1

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

摘要翻译

中文标题: 一个带有阶数接近零的非局部算子的逻辑斯蒂方程

摘要翻译: 本文研究一个由阶数“接近零”的非局部算子驱动、并具有混合 Dirichlet-Neumann 边界条件的方程,该方程建模了一个在封闭环境中自我竞争资源的种群,并受到极端重尾扩散过程的影响。我们建立了若干关于种群可存活状态存在性的结果,这些结果依赖于非线性反应项的不同假设,以及 Dirichlet 条件或 Neumann 条件所作用的集合的各种配置。


7. 通过Bruhat-Tits树扩张得到的广义Vladimirov导数的能量关系

Energy relations for the generalised Vladimirov derivative via Bruhat-Tits tree extensions

作者: An Huang, Yaojia Sun

类别: math.AP, math.FA

发布日期: 2026-07-27

arXiv ID: 2607.24316v1

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

摘要翻译

中文标题: 通过Bruhat-Tits树扩张得到的广义Vladimirov导数的能量关系

摘要翻译: 我们在将边界 $\partial T_p$ 与 $\mathbb{P}^1(\mathbb{Q}_p)$ 等同后,以能量关系为框架,在Bruhat-Tits树 $T_p$ 上构造了分数阶Laplacian的非阿基米德类比——Caffarelli-Silvestre扩张,并证明了该扩张所诱导的边界函数经尺度变换后的广义Vladimirov导数。


8. 超解构造与分数阶拉普拉斯算子的最优Hardy不等式

Supersolution Construction and Optimal Hardy Inequality for Fractional Laplacians

作者: Philipp Hake, Matthias Keller, Felix Pogorzelski

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

发布日期: 2026-07-27

arXiv ID: 2607.24169v1

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

摘要翻译

中文标题: 超解构造与分数阶拉普拉斯算子的最优Hardy不等式

摘要翻译: 我们给出一个新准则,通过超解构造来证明图上一般算子的Hardy权的最优性。对于图上无杀死项的拉普拉斯算子,该准则总能借助Green函数给出一个最优Hardy权,且无需任何额外假设。此外,与早期结果不同,我们的结果不限于局部有限图。这使得我们尤其能够获得一般图上分数阶拉普拉斯算子的最优Hardy权。对于欧氏格点 $\mathbb{Z}^d$ 上的分数阶拉普拉斯算子,我们进而得到一个最优Hardy权,其常数及渐近行为与连续情形所预期的一致。


9. Airy 方程在直线上的 Hausdorff 型时间迹可观测性与环面上的点可观测性

Hausdorff type Time-Trace Observability for Airy Equations on the Line and Point Observability on the Torus

作者: Ze Li

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24076v1

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

摘要翻译

中文标题: Airy 方程在直线上的 Hausdorff 型时间迹可观测性与环面上的点可观测性

摘要翻译: 本文的主要结果有三方面。第一,我们证明了 Airy 方程 $\partial_t u + \partial_x^3 u = 0$ 在实直线 $\mathbb{R}$ 上对任意观测时间 $T>0$,从 Hausdorff-thick 集合出发、在时间迹意义下的可观测性不等式。观测泛函定义为在 Hausdorff-thick 集上取 $L^2(0,T)$ 时间迹的块上确界。第二,我们证明了 Airy 方程在实直线上对某些周期集合的可观测性不等式,特别地,这蕴含了一类空间点列上的点可观测性。第三,我们给出了带实值有界势 $V \in L^\infty(\mathbb{T})$ 的 Airy 方程在环面 $\mathbb{T}$ 上有限点可观测性的充要条件:对有限观测集 $F$,存在一个有限维不变子空间上的 Kalman 秩条件。作为推论,我们得到了 Airy 方程与线性 KdV 方程在环面上的尖锐点可观测性结果。此外,对满足有限阶正则性假设的势 $V$,我们证明:任何具有聚点的观测集(特别是任何具有正 Hausdorff 维数的集合)均对每个观测时间 $T>0$ 给出可观测性不等式。 由于 Airy 方程既不具备高频指数衰减性,也不具备逐点光滑效应——而这两者是近期关于热方程的 Hausdorff 型可观测性研究中的关键性质——我们引入了若干适配 Airy 情形的新思想。


10. 关于Boltzmann方程的一个自由边界问题

On a free boundary problem of the Boltzmann equation

作者: Shengchuang Chang, Renjun Duan, Shuangqian Liu

类别: math.AP, math-ph

发布日期: 2026-07-27

arXiv ID: 2607.24053v1

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

摘要翻译

中文标题: 关于Boltzmann方程的一个自由边界问题

摘要翻译: 本文研究Boltzmann方程的一个自由边界问题,该问题建模稀薄气体与运动壁面的相互作用——即动理学理论中经典的活塞问题。活塞运动由牛顿定律支配,其所受阻力来自气体施加的拖曳力,且可附加或不附加Hookean恢复力。此类问题的核心难点在于动理学方程与自由运动边界的强耦合;由于解的低正则性,经典的Lagrangian表述不可用。我们的方法依赖于两个新要素:引入一个共形变换,将运动区域约化为固定区域;并揭示了一种耦合的能量结构,该结构联结了动理分布函数与自由边界变量,从而在界面处展现出内禀的耗散与相消机制。这些结构性观察导出了全耦合系统的全局非线性理论。结果,我们在全局Maxwellian附近建立了解的全局存在性、唯一性、非线性稳定性以及指数收敛至平衡态。这为活塞型全耦合Boltzmann自由边界问题提供了首个严格的全局适定性理论。


11. 欧几里得球上的局部化框架

Localized frames on Euclidean balls

作者: Kevin Hughes, Arie Israel, Azita Mayeli

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

发布日期: 2026-07-27

arXiv ID: 2607.23953v1

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

摘要翻译

中文标题: 欧几里得球上的局部化框架

摘要翻译: 我们构造了适配于欧几里得球的显式波包框架,并利用该框架获得时空–谱截断算子的定量特征值估计。设 $d\geq 2$,令 $B_d(R)\subset \mathbb{R}^d$ 为半径为 $R$ 的欧几里得球,令 $S\subset \mathbb{R}^d$ 为一个可测集,使得 $\partial S$ 具有有限的 $(d-\eta)$-上Minkowski容度,其中 $0 < \eta\leq 1$。我们构造了 $L^2(B_d(R))$ 上的一个单位范数框架,其框架界仅依赖于维数 $d$,且其元素适配于该球的径向与角向几何结构。我们证明了该框架的定量Fourier局域化估计:相对于 $S$,该框架分解为三类元素——集中于 $S$ 的波包、集中于 $\mathbb{R}^d\setminus S$ 的波包,以及一个例外族,其基数被显式地控制在 $R$ 与 $\partial S$ 的Minkowski容度的函数之下。作为应用,我们导出了与集合 $B_d(R)$ 和 $S$ 相关联的时空–谱截断算子之陡降区(plunge region)的一个上界。


12. 相空间延拓的动理学Cucker-Smale模型的弱稳定性与随机平均场极限

Weak stability and random mean-field limit of the phase-spatially extended kinetic Cucker-Smale model

作者: Seung-Yeal Ha, Xinyu Wang

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.23906v1

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

摘要翻译

中文标题: 相空间延拓的动理学Cucker-Smale模型的弱稳定性与随机平均场极限

摘要翻译: 我们在相空间延拓的框架下研究动理学Cucker-Smale(简记为KCS)模型的弱稳定性,该模型可形式地由无穷粒子Cucker-Smale模型在平均场极限下导出。对有界Lipschitz连续的通讯权重函数,我们导出了具有指数速度尾和有限空间二阶矩的测度值解的有限时间Osgood型弱稳定性。与相空间受限情形不同,KCS模型的解算子关于初始数据并非Lipschitz连续;这是由于速度尾无界,从而导致对齐力缺乏一致的Lipschitz界。作为弱稳定性的应用,我们得到一个i.i.d.采样推论:由独立初始样本生成的经验测度,在任意有限时间区间上,以期望意义收敛于相应动理学模型的测度值解。


13. 一类具有部分反平方势的色散方程的衰减估计

Decay estimates for a class of dispersive equations with partial inverse-square potentials

作者: Jiabin Qian, Manli Song

类别: math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23578v1

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

摘要翻译

中文标题: 一类具有部分反平方势的色散方程的衰减估计

摘要翻译: 令 $\mathcal{L}_a=-Δ_x-Δ_y+\frac{a}{2}|x|^{-2}$ 表示定义在 $L^2(\mathbb{R}^2_x\times \mathbb{R}^n_y)$ 上的薛定谔算子,其中包含一个奇异的部分反平方势。本文目的有二:第一,借助 Zhang-Zhang [J. Geom. Anal. \textbf{35}(3), Paper No. 71, 27pp (2025)] 所建立的关于算子 $\mathcal{L}_a$ 的谱测度的显式表达式,我们研究形如 $e^{itφ(\sqrt{\mathcal{L}_a})}$ 的一类色散半群的衰减估计,其中 $φ: \mathbb{R}^+ \to \mathbb{R}$ 是光滑函数;为处理由相位函数 $φ$ 的非齐次性所引发的技术困难,我们采用频率局域化与平稳相位法。本文第二部分首先推导分数阶薛定谔算子 $e^{it\mathcal{L}_a^ν}$($0<ν\neq\frac{1}{2}$)的边界 Strichartz 估计;作为所建立衰减估计的应用,我们进一步获得若干与算子 $\mathcal{L}_a$ 相关联的具体波动方程的 Strichartz 估计,对应于 $φ(r)=r, r^2, r^2+r^4, \sqrt{1+r^2}, \sqrt{1+r^4}$ 以及 $r^μ,0<μ\neq 1$。尤为值得注意的是,我们的结果统一并简化了现有针对算子 $\mathcal{L}_a$ 的色散估计,同时将相关理论推广至更一般的情境。


14. 环面 $\mathbb{T}$ 上带阻尼的 BBM 方程的周期解及其稳定性

Periodic solution and its stability of a damped BBM equation posed on $\mathbb{T}$

作者: Chun Ho Lau, Taige Wang

类别: math.AP, math.DS

发布日期: 2026-07-26

arXiv ID: 2607.23431v1

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

摘要翻译

中文标题: 环面 $\mathbb{T}$ 上带阻尼的 BBM 方程的周期解及其稳定性

摘要翻译: 本文研究环面 $\mathbb{T}$ 上一类带阻尼的 Benjamin-Bona-Mahony (BBM) 方程的时间周期解的存在性与稳定性,其内部施加了具有时间周期 $θ$ 的周期外力 $f(x, t)$。此类解在 Sobolev 空间 $H^{\ell}, \ell\ge 0$ 中被建立;值得注意的是,当追求低正则性情形 $\ell\in[0, 1)$ 时,采用了 I-能量方法。


15. 高维多孔介质方程带燃烧反应的尖锐阈值

Sharp Thresholds for the Porous Medium Equation with a Combustion Reaction in Higher Dimensions

作者: Maolin Zhou

类别: math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23429v1

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

摘要翻译

中文标题: 高维多孔介质方程带燃烧反应的尖锐阈值

摘要翻译: 我们研究带燃烧型反应的多孔介质方程

$$ u_t=Δu^m+f(u),\qquad x\in\mathbb R^N,\ t>0, $$

在径向、非负、紧支集初值条件下的解。我们建立了有界解长时间行为的完整分类。在维数 $N=2$ 时,每个这样的解局部一致收敛于常数 $0$、$θ$ 或 $1$ 之一;对有序初值族,存在唯一阈值参数将熄灭(vanishing)与传播(spreading)分离开,且临界解收敛于点火温度 $θ$。在维数 $N\ge3$ 时,在基态(ground state)中心值集合满足自然的全不连通性(total-disconnectedness)条件下,每个有界解局部一致收敛于 $0$、$θ$、$1$ 或某个径向基态 $U\in\mathcal S$ 之一。此外,借助一种新颖的归一化环形扰动论证(normalized annular perturbation argument),我们在 $N\ge3$ 时排除了点火态 $θ$ 作为过渡极限(transition limit)的可能性。对于过渡解(transition solutions),我们给出了自由边界传播速度的估计:在二维情形下,

$$ b(t)\asymp \frac{\sqrt t}{(\log t)^{\frac{m-1}{2m}}}, $$

而在维数 $N\ge3$ 且极限为基态时,

$$ b(t)\asymp t^{\frac{m}{N(m-1)+2}}. $$

这些结果揭示了退化燃烧动力学中尖锐的维数二分现象(sharp dimensional dichotomy)。


16. Oseen–Frank 模型中液晶与电场相互作用的变分原理

Variational principles for the interaction of liquid crystals and electric fields in the Oseen–Frank model

作者: Giovanni Di Fratta, Valeriy Slastikov, Arghir Zarnescu

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23315v1

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

摘要翻译

中文标题: Oseen–Frank 模型中液晶与电场相互作用的变分原理

摘要翻译: 我们为单轴向向列相液晶在单常数 Oseen–Frank 近似下与外加电场相互作用建立了严格的变分框架。平衡构型由一个具有挑战性的非局部、非线性能量所支配,该能量具有 min-max 鞍点结构。我们的第一个主要结果将该问题重述为一个纯粹的双重最小化问题:利用凸对偶性和 Hodge 分解,我们将标量静电势替换为矢量势,从而得到一个具有唯一最小化子的闭式对偶泛函。这一直接能量最小化原理对理论分析和数值模拟均具优势。我们的第二个结果严格量化了介电各向异性趋于零时静电反作用的解耦行为:我们在精确的非局部能量与其标准局部近似之间建立了统一的二次能量估计,形式上证明了物理界广泛采用的忽略感应去极化场这一启发式做法的合理性。最后,在严格的强制性条件下,我们结合 $Γ$-收敛、一致 Sobolev 正则性以及微扰意义下的强制性传递,证明物理最小化子以最优的、定量的线性速率收敛于极限调和映射。


17. 图上的广义多孔介质方程:适定性、解的熄灭与质量守恒

The generalized porous medium equation on graphs: well-posedness, extinction, and mass conservation

作者: Davide Bianchi, Bobo Hua, Alberto G. Setti, Radosław K. Wojciechowski

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23091v1

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

摘要翻译

中文标题: 图上的广义多孔介质方程:适定性、解的熄灭与质量守恒

摘要翻译: 我们研究无限加权图上广义多孔介质方程的柯西问题。对一般的非线性项,我们在有限子图上建立狄利克雷比较原理与弱极大值原理,并通过穷竭论证构造出任意 $\ell^\infty$ 初值数据的极小与极大逐点解,这些解被显式的、可能依赖于时间的上下界所控制。 对多孔型非线性项 $φ(s)=s|s|^{m-1}$,在假设存在 $ν$-Sobolev 不等式且 $ν>2$ 的条件下,我们为 $\ell^1$-温和解导出定量的能量估计。这些估计给出快扩散范围 $0<m<2/ν$ 内的有限时间熄灭现象,以及在 $m>2/ν$ 范围内的 $\ell^1$-$\ell^q$ 光滑效应。有趣的是,我们在若干模型图上恢复了欧氏情形的临界指数。最后,我们在无穷远处随机完备的图上,对非负 $\ell^1$-温和解证明了一个精确的广义质量平衡公式,该公式允许存在任意的 killing 项;对适当的经典解与有界逐点解,同样的质量平衡亦成立。当不存在 killing 项时,该平衡退化为质量守恒。


18. 格点上的渐近尖锐 Hardy-Rellich 不等式

Asymptotically sharp Hardy-Rellich inequalities on lattices

作者: Xia Huang, Dong Ye

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23057v1

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

摘要翻译

中文标题: 格点上的渐近尖锐 Hardy-Rellich 不等式

摘要翻译: 我们确定了格点 $\mathbb{Z}^d$ 上离散 Hardy-Rellich 不等式中最优常数的渐近尖锐行为。对每个固定的整数 $m\ge 1$,令 ${\mathcal C}_{m,d}$ 为 $m$ 阶 Hardy-Rellich 不等式的最佳常数。我们证明

$$ \lim_{d\to\infty}\frac{\mathcal C_{m,d}}{d^m}=2^m. $$

我们的方法结合了 Fourier 约化至带权不等式(定义在平环面 $\mathbb{T}^d$ 上)与一阶、二阶广义带权 Hardy-Rellich 恒等式。一个关键的新颖之处在于运用概率集中估计——特别是 Hoeffding 不等式与熵方法——以维度一致的方式处理涉及各向异性权 $ω^γ$ $(γ\geq 1)$ 的估计,其中

$$ ω(x)=\sum_{j=1}^{d} \Big(\sin\frac{x_j}{2}\Big)^2. $$

这些工具导出了环面上的渐近尖锐带权估计,再经迭代即得格点上的不等式。


19. 具有域壁的狄拉克方程的色散估计

Dispersive estimates of Dirac equations with a domain wall

作者: M. Burak Erdogan, Joseph Kraisler, Amir Sagiv

类别: math.AP, math-ph

发布日期: 2026-07-24

arXiv ID: 2607.22879v1

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

摘要翻译

中文标题: 具有域壁的狄拉克方程的色散估计

摘要翻译: 我们建立了二维具有有界与无界域壁的狄拉克方程的色散衰减估计,并同时给出了相应一维问题的估计。这些是拓扑绝缘体中体-边对应关系的典型模型,但其色散动力学此前尚未被研究。我们证明了拓扑等价的模型可能展现出定性上不同的动力学行为。我们的新方法使得我们能够分析超出通常常系数狄拉克算子局部摄动框架的模型。


不相关论文(基本信息)

1. 通过Morse理论定义的曲面上的Yang-Mills测度

The Yang-Mills measure on surfaces via Morse theory

作者: Reda Chhaibi, Nguyen Viet Dang, Yannick Guedes Bonthonneau, Gabriel Rivière, Tat Dat Tô

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

发布日期: 2026-07-27

arXiv ID: 2607.24640v1

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

摘要: We introduce a Morse theoretical approach to the construction of the Yang–Mills measure on the space of connections of a compact Riemannian surface. This provides a direct continuous version of this measure which was previously obtained through lattice approximations by Chevyrev in the case of the flat torus and by one of the authors and Nohra for general compact Riemannian surfaces. The starting point is the new notion of a Morse gauge together with the resolution of random cohomological equations associated to Morse–Smale vector fields. This is achieved by improving exponential convergence to equilibrium results for Morse–Smale gradient flows that were obtained by two of the authors in the context of the study of Ruelle spectra and by Jia, Stewart and Sverak in the context of simplified models from fluid mechanics. Combining these random solutions with the data given by the Morse complex, we introduce a free Yang-Mills measure on space of connections and, using classical tools from stochastic differential equations, we show how to make sense of holonomies for random connections along a large class of curves. Finally, by setting a proper conditioning of this free measure through these random holonomies, we define the Yang–Mills measure and we compute its partition function together with the law of random holonomies with respect to this measure, recovering the formulas from the works of Migdal, Witten and Lévy.


2. Derangetropy 算子

Derangetropy Operators

作者: Masoud Ataei, Sepideh Forouzi

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

发布日期: 2026-07-27

arXiv ID: 2607.24705v1

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

摘要: A derangetropy operator reweighs a probability density by a fixed profile of its own cumulative distribution function, acting through ranks alone. We prove that these operators are precisely the transformations of absolutely continuous laws equivariant under monotone changes of variable, and that they compose through interval maps, making their dynamics exactly solvable: iteration condenses every law onto its median with a universal Koenigs limit law, the continuous flow is solvable in closed form with the Cauchy family as invariant hyperbolic manifold, and balanced against diffusion the distribution function obeys an overdamped sine-Gordon equation whose unique steady law, the hyperbolic secant, is a globally stable kink. A variational principle selects the canonical kernel, the squared ground state of the Dirichlet Laplacian on the unit interval, whose update is a Bayesian posterior costing exactly one bit for every law, and a unitary lift equips each law with an isospectral Sturm-Liouville system and a quantum carpet: a strong law of large numbers for the critical Sobolev mass, with sharp constant given by Wiener’s jump statistic, proves the fractality of the Schrödinger density for arbitrary real bounded-variation data with a jump, previously known only for rational step data. A lift to the dual of the Virasoro algebra, at central charge one, a normalization, identifies laws with disconjugate Hill potentials, tails with conformal weights, and the Cauchy family with the exceptional orbit. Randomized phases yield an exact sampler and a multiplicative chaos, log-correlated in the weakly tempered regime; in several dimensions, dependence stratifies into a conserved interaction, a torsion driven by maximal correlation, and a flat Sinkhorn transport.


3. Tikhonov 正则化的物理信息神经网络用于抛物型偏微分方程终态分布式最优控制

Tikhonov-Regularized Physics-Informed Neural Networks for Terminal-State Distributed Optimal Control of Parabolic Partial Differential Equations

作者: Nguyen Thanh Quang, Ta Thi Thanh Mai

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24572v1

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

摘要: In this study, we propose a Physics-Informed Neural Networks (PINNs) framework that incorporates Tikhonov regularization to solve terminal-state tracking optimal control constrained by parabolic partial differential equations (PDEs). This problem is inherently ill-posed, as infinitely many distributed controls may drive the system to the same desired state, so the regularization guides the optimizer toward the minimum-energy control, restoring numerical stability and yielding a smooth, physically meaningful solution. On the theoretical side, we establish a consistency result showing that PINNs minimizers nearly attain the continuous regularized objective under residual and quadrature approximation assumptions, and a novel error estimate that bounds the deviation of the learned control from the minimum-energy solution in terms of the PINNs training residuals and the regularization parameter. Numerical experiments on the linear heat equation and the nonlinear Burgers’equation demonstrate that the regularized PINNs framework accurately achieves the target terminal state while producing controls with significantly lower energy and smoother profiles compared to unregularized baselines.


4. Gompertz 增长模型下肿瘤总种群数量的优化

Optimization of the total tumor population under Gompertz growth

作者: Iulia Martina Bulai, Francesca Gladiali, Benedetta Pellacci

类别: math.AP, math.NA

发布日期: 2026-07-27

arXiv ID: 2607.24468v1

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

摘要: We study optimal control problems for a stationary reaction–diffusion model describing the spatial distribution of a tumor cell population with Gompertz growth. The control $m(x)$ represents a treatment term acting as a density-dependent removal rate and it is subject to $L^{1}-L^{\infty}$ constraints. When the intrinsic growth rate is constant, the uniform distribution of the treatment is shown to be the unique minimizer. For the maximization problem, we prove that every optimal control is of bang-bang type. In addition, we show that in the one dimensional case and for sufficiently large diffusion rates, the positivity set of optimal controls is an interval sticking to one of the extrema of the domain. Finally, numerical simulations complement the theoretical analysis and explore regimes that are not fully covered by the results proved in the paper. The computations confirm the bang-bang structure of maximizers, and illustrate how the shape of optimal controls and the associated states are affected by spatial heterogeneity in the growth rate, localized admissible treatment regions, and the diffusion coefficient. Moreover, they reveal a monotone dependence of the optimized total population on the diffusion coefficient: this is a new phenomenon with respect to the logistic setting.


5. 一维椭圆方程(含奇异一阶散度项)解的新定义

New definitions of a solution\break of a one-dimensional elliptic equation\break with a singular first order divergence term

作者: Daniela Giachetti, Pedro Jesús Martínez-Aparicio, François Murat, Francesco Petitta

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24308v1

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

摘要: In this paper, which is a continuation of our recent paper [1], we give two new definitions of a weak solution of the one-dimensional, elliptic, nonlinear singular problem which is formally written as

$$ \begin{cases} \displaystyle -\frac{d}{dx}\left(a(x) \frac{d u}{dx}\right) = - \frac{d φ(u) }{dx}- \frac{d g(x) }{dx}& \text{in}\;(0,L),\\ \newline u(0)=u(L)=0\,, & \end{cases} $$

where $φ:\mathbb{R}\mapsto\mathbb{R}\cup\{+\infty\}$ is a singular function whose model is given by $\displaystyleφ(s)=φ_γ(s)=\frac{1}{|s|^γ}$ with $γ>0$. In these two definitions the solutions belong respectively to the space $u\in W_0^{1,1}(0,L)$ with $φ(u)\in L^1(0,L)$, and to the space $u\in W_0^{1,m}(0,L)$ with $φ(u)\in L^m(0,L)$ for $m>1$. In these frameworks we state and prove results of non-existence, existence, and non-isolation of the solutions.


6. 紧密排列的高折射率纳米粒子之间的非奇异热点

Non-singular hotspots between closely spaced high-index nanoparticles

作者: Konstantinos Alexopoulos, Bryn Davies, Pierre Millien

类别: math.AP, math-ph, math.SP, physics.optics

发布日期: 2026-07-27

arXiv ID: 2607.24204v1

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

摘要: We study the concentration of the field between two nearly touching high-index dielectric resonators in three dimensions. The model is a scalar Helmholtz transmission problem in the resonant regime, wherein the wavelength inside the resonators of the same order as their typical diameter. The material contrast enters only a lower-order term of the equation and not its principal part. As a consequence, the gradient of the field stays bounded independently of the distance separating the particles, and does not blow up. Nevertheless, simulations and experiments show that the gradient still concentrates in the gap: as the particles approach, it grows like the inverse of their separation over an intermediate range of distances and then saturates once they are very close. We explain this pre-asymptotic effect through a weak-coupling regime, in which the resonant modes of the pair are perturbations of the modes of each isolated particle. When such a mode keeps a nonzero contrast between the two facing boundaries, a mean value argument across the gap accounts for the growth of the gradient; this requires strengthening the standard spectral perturbation estimates from an average to a pointwise control. The growth stops once the interaction between the particles is no longer weak and the perturbation theory fails. Numerical experiments confirm the transition from amplification to saturation.


7. 多重与奇异界面附近的边缘态传播

Edge State Propagation Near Multiple and Singular Interfaces

作者: Nicolas Frantz, {É}ric Vacelet

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24069v1

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

摘要: We study the semiclassical propagation of edge-states for a two-dimensional Dirac operator with a spatially varying mass. The zero set of the mass models the interface between topological phases and is allowed either to consist of two disjoint smooth curves or to possess an isolated singular point. For disjoint interfaces, we prove that a wave packet initially localized on one component remains confined near this interface over long semiclassical time scales and we quantify the influence of the distance between the two components on the accuracy of the approximation. For singular interfaces, we determine the regime in which the wave-packet approximation remains valid as the packet approaches the singularity and identify the scaling at which our estimates lose their uniformity.


8. 拟线性 Hénon 型 $N$-拉普拉斯李维-方程的非径向解

Non-radial solutions for the quasi-linear Hénon type $N$-Laplacian Liouville equation

作者: Wei Dai, Lixiu Duan, Changfeng Gui, Yuan Li

类别: math.AP

发布日期: 2026-07-27

arXiv ID: 2607.24012v1

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

摘要: In this paper, we investigate the following quasi-linear weighted $N$-Laplacian Liouville equation \begin{equation*}\label{0} -Δ_N u=|x|^{Nα}e^{u}, \qquad x\in \R^N, \end{equation*} where $N \geq 2$. For $\al>0$, by carefully studying the linearized problem and applying the approximation method and bifurcation theory, we prove that, when the parameter $α$ equals to the critical values $α(k):=\frac{\sqrt{k(N-1)(k+N-2)}}{N-1}-1$ for $k \geq 2$, there exist non-radial solutions $u$ (bifurcating from $U_{α(k)}$) to the above quasi-linear Hénon type Liouville equation such that $u\sim \ln|x|$, $|\nabla u|= O(|x|^{-1})$ at $\infty$ and $\int_{\R^N}|x|^{Nα}e^{u}\md x=N\left(\frac{N^2}{N-1}\right)^{N-1}(α+1)^{N-1}ω_N$. One should note that, $α(k)=k-1$ for $k\geq2$ when $N=2$. Our results successfully extend the existence result of J. Prajapat and G. Tarantello in \cite{PT} concerning the $2$-dimension and Laplacian case (i.e., $N=2$) to the more general $N$-dimension and $N$-Laplacian cases ($N\geq 2$), and extend the results of F. Gladiali, M. Grossi, and S. L. N. Neves in \cite{GGN} and the authors in \cite{DDGL} from $1<p<N$ to the much more complicated limiting case $p=N$. We introduced some new ideas and overcame a series of crucial difficulties, including the nonlinearity nature of the $N$-Laplacian $Δ_N$, the lack of Green integral representation formula and critical weighted Sobolev embedding inequality, the absence of Kelvin type transforms for linearized/difference equations, the invariance of the total mass under scalings of $u$, and the signs-changing and divergence (to $-\infty$) at $\infty$ of the solutions, which makes the suitable choices of the approximate problems, the (normalized) approximate function sequences and the working space to be quite difficult.


9. 具有给定切流的闭平均曲率流

Closed mean curvature flows with prescribed tangent flows

作者: Jingwen Chen, Tang-Kai Lee, Ao Sun, Xinrui Zhao

类别: math.DG, math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23879v1

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

摘要: Given an embedded shrinker $Σ$ in $\mathbb{R}^{n+1}$ that is either closed, asymptotically conical, or a Cartesian product of such a shrinker with $\mathbb{R}^k$, we construct a closed embedded mean curvature flow whose tangent flow at the first singularity is modeled on $Σ$. We also prescribe the first-order asymptotics of the tangent flow. This result is a consequence of a more general theorem that allows us to construct mean curvature flows with an additional force whose tangent flow and first-order asymptotics at the first singularity are prescribed.


10. 紧致 Hermitian 流形上复 Hessian 方程的一致 $L^ fty$ 估计

Uniform $L^\infty$ estimates for complex hessian equations on compact Hermitian manifolds

作者: Truong Dinh Dat

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

发布日期: 2026-07-26

arXiv ID: 2607.23724v1

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

摘要: We develop a pluripotential approach to complex Hessian equations on compact Hermitian manifolds. In this setting, the lack of closedness of the background metric introduces torsion terms that prevent a direct extension of the Kähler theory. Our main result is a uniform $L^\infty$ estimate for bounded $ω$-$m$-subharmonic solutions of the equation

$$ (ω+ dd^c u)^m \wedge ω^{n-m} = cf\,ω^n, $$

under the assumption that $f \in L^p$, $f \ge 0$ for some $p>1$. The proof combines a weak comparison principle with torsion error, a capacity theory adapted to the Hermitian setting, and a nonlinear iteration scheme controlling the decay of sublevel sets. As applications, we obtain existence, stability and compactness results for weak solutions with $L^p$ densities. These results extend several aspects of the pluripotential theory. of complex Hessian equations beyond the Kähler framework.


11. 通过覆叠空间求解Plateau问题

Plateau’s Problem via covering spaces

作者: James Tissot

类别: math.DG, math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23703v1

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

摘要: In 1995, Brakke proposed a formulation of the Plateau problem for a given boundary $Γ$ that allows for triple junctions and tetrahedral singularities. Let $Γ$ be a smooth closed curve and let $G=π_1(\mathbb{R}^3\setminus Γ)$. For each proper finite index subgroup $N$ of $G$, Brakke constructs a $(\mathrm{\mathbf{M}}, 0, \infty)$-minimal surface $Σ_N$, which is obtained as the projection of the boundary of a perimeter-minimising fundamental domain in the covering space associated to $N$. We advance the theory in two ways. Firstly, we extend Brakke’s construction to include all normal subgroups $N\triangleleft G$, i.e. possibly with infinite index. Secondly, we prove a compactness result which implies that there exists a proper normal subgroup $N_0\triangleleft G$ such that

$$ \mathrm{Area}(Σ_{N_0})=\inf_{\{N\triangleleft G, N \neq G\}}\ \mathrm{Area}(Σ_N). $$

A similar result holds when $Γ$ has many connected components. Furthermore, we study the spanning and minimising properties of the $Σ_N$s and $Σ_{N_0}$.


12. 非凸区域上扭转函数的 Hessian 退化性与非星形超水平集

Hessian degeneracy and non-star-shaped superlevel sets of the torsion function on non-convex domains

作者: Xiuda Liang, Peng Luo, Wenjie Wang

类别: math.AP

发布日期: 2026-07-26

arXiv ID: 2607.23655v1

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

摘要: We construct a sequence of bounded $C^\infty$-smooth, simply connected and non-convex planar domains $Ω_{a,ε}$ for which the torsion function

$$ -Δu_{a,ε}=1\quad\text{in }Ω_{a,ε},\qquad u_{a,ε}=0\quad\text{on }\partialΩ_{a,ε} $$

has a strict global maximum point at origin $0$ satisfying

$$ λ_{max}\bigl(D^2u_{a,ε}(0)\bigr)\longrightarrow 0^-. $$

Moreover, the ratios $\text{diam}(Ω_{a,ε}) / \text{inrad}(Ω_{a,ε})$ remain uniformly bounded. Hence the Hessian estimate proved by Steinerberger (J. Funct. Anal. 274, 1611–1630, 2018) for convex planar domains cannot be extended to the class of smooth simply connected non-convex domains. \vskip0.2cm Moreover, the domains we construct are doubly symmetric, only one-directionally convex, but star-shaped. We show that some superlevel sets of the torsion function fail to be star-shaped when the slits are sufficiently long, hence we can also give some comments on the question raised by Gladiali and Grossi (Amer. J. Math. 144, 1221–1240, 2022) from the perspective of strengthened geometric conditions.


13. 关于特殊拉格朗日曲率方程的一些反例

Some counterexamples for the special lagrangian curvature equation

作者: Guohuan Qiu, Guanyu Tao

类别: math.AP, math.DG

发布日期: 2026-07-26

arXiv ID: 2607.23592v1

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

摘要: We construct three counterexamples for the special Lagrangian curvature equation (SLCE). First, in dimension two, we use a post-focal branch of a parallel surface with constant positive Gauss curvature to construct an explicit Lipschitz viscosity solution which is not $C^1$. Second, still in dimension two, we construct a sequence of smooth admissible solutions on a fixed rectangle with uniformly bounded $C^1$-norm but unbounded curvature at one point; furthermore, we show that any uniform $C^{1,β}$ estimate fails for $β> 1/3$. Third, in dimension three and in the subcritical phase, we construct a Mooney-Savin type Lipschitz viscosity solution whose gradient has a jump discontinuity across an analytic surface. These examples demonstrate the sharpness of the recent a priori estimates by Qiu and Zhou, revealing that their structural assumptions-convexity and the critical phase-are strictly necessary.


14. 一个 $p$-进乘性微积分的伪微分算子

A Pseudo-Differential Operator of the $p$-Adic Multiplicative Calculus

作者: Alexandra V. Antoniouk, Anatoly N. Kochubei, Mariia V. Serdiuk

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23256v1

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

摘要: We introduce and study a pseudo-differential operator $W^α$, defined via the $p$-adic Mellin transformation, acting on real- or complex-valued functions on the multiplicative group of the field of $p$-adic numbers. The operator $W^α$ is maximally dissipative. An estimate of its heat kernel is given.


15. Hölder 不等式的方差表示公式

A Variance Representation Formula for Hölder’s Inequality

作者: Li-Chang Hung

类别: math.PR, math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23255v1

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

摘要: We establish an exact representation formula for the deficit in Hölder’s inequality. Rather than estimating the deficit by auxiliary quantities, we show that it can be represented as an accumulated variance along the natural exponential interpolation connecting the two endpoint densities. More precisely, the logarithmic Hölder deficit is expressed as the integral of the variance of the logarithmic density ratio weighted by the Green kernel associated with the one-dimensional interpolation parameter. This identity reveals that Hölder’s inequality is a consequence of the convexity of a logarithmic partition function and provides an intrinsic interpretation of the deficit as an interpolation energy. The representation suggests a broader framework for studying functional inequalities through variance identities.


16. 几何结构的流 II

Flows of geometric structures II

作者: Daniel Fadel, Udhav Fowdar, Eric Loubeau, Andrés J. Moreno, Henrique N. Sá Earp

类别: math.DG, math.AP

发布日期: 2026-07-25

arXiv ID: 2607.23231v1

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

摘要: We advance the general theory of flows of tensorial $\mathrm{H}$-structures, focusing on non-isometric flows and on the case $\mathrm{H}=\mathrm{SU}(m)\subset\mathrm{SO}(2m)$. After developing the relevant $\mathrm{SU}(m)$ algebra, we compare two natural evolutions: the unrestricted negative gradient flow of the intrinsic-torsion energy and a Ricci-harmonic flow. We prove short-time existence and uniqueness for the Ricci-harmonic $\mathrm{H}$-flow, with arbitrary lower-order torsion-quadratic terms, for every closed subgroup $\mathrm{H}\subset\mathrm{SO}(n)$. For groups for which the projection to $\mathfrak{h}^\perp$ defines a $4$-form, including $\{1\}$, $\mathrm{SU}(2)$, $\mathrm{G}_2$, and $\mathrm{Spin}(7)$, we express the negative gradient flow in Ricci-harmonic form up to explicit lower-order torsion terms and prove short-time existence and uniqueness by a modified DeTurck argument. We treat the genuinely different $\mathrm{SU}(m)$ case by a separate principal-symbol computation, proving short-time existence and uniqueness for the unrestricted negative gradient flow of $\mathrm{SU}(m)$-structures. The same computation identifies the natural negative gradient flow of $\mathrm{U}(m)$-structures as a borderline case, which cannot be made strictly parabolic by first-order diffeomorphism gauges. For the modified Ricci-harmonic flow, we derive heat-type evolution equations for the intrinsic torsion, a doubling-time estimate and Shi-type derivative estimates for $(|\mathrm{Rm}|^2+|\nabla T|^2+|T|^4)^{1/2}$, and a finite-time continuation criterion. In dimension six, we translate the formalism into the standard torsion forms of an $\mathrm{SU}(3)$-structure and describe, to highest order, the corresponding family of second-order quasilinear $\mathrm{SU}(3)$-flows.


17. $\mathcal{P}(Φ)_2$ 理论源自多体量子Gibbs态

$\mathcal{P}(Φ)_2$ Theory from many-body quantum Gibbs states

作者: Phan Thành Nam, Zhilin Yang, Xiangchan Zhu

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

发布日期: 2026-07-25

arXiv ID: 2607.23084v1

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

摘要: We derive the $\mathcal{P}(Φ)_2$ measure on the two-dimensional unit torus as the rigorous limit of many-body quantum Gibbs states. In particular, the quantum model corresponding to the $Φ^{2p}_2$ measure is formulated in the grand-canonical ensemble with a general symmetric $p$-body interaction potential which is not required to have a factorized form. Unlike in the $Φ^4_2$ case investigated by Fröhlich–Knowles–Schlein–Sohinger in \cite{FKSS25}, in the treatment of higher-order interactions, Wick renormalization generates a full hierarchy of lower-order interactions which we organize systematically using a graphical formalism. One key ingredient of our analysis is a logarithmic stability estimate for both the nonlocal Hartree functional and the many-body Hamiltonian that is uniform in the interaction range, allowing us to use the positivity of the leading $p$-body interaction to control all lower-order terms generated by the Wick counterterms.


18. 来自具有三体相互作用的玻色子Gibbs态的非局部三次密度Gibbs测度

Nonlocal Cubic Density Gibbs Measures from Bosonic Gibbs States with Three-Body Interactions

作者: Hao Liang

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

发布日期: 2026-07-25

arXiv ID: 2607.23041v1

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

摘要: We study the high-temperature mean-field limit of grand-canonical bosonic Gibbs states on the torus with renormalized nonlocal three-body interactions. In dimensions two and three, we construct the limiting nonlinear classical Gibbs measure and prove convergence of the relative free energy and of the reduced density matrices of every fixed order; in three dimensions, a smallness condition on the interaction is imposed. The proof combines the density-channel representation of the interaction with a coherent-state variational method based on the upper-symbol representation of the free Gibbs state. The same framework also contains, as a special case, the homogeneous positive-type model studied by Lewin, Nam, and Rougerie (2021).


19. 抛物-椭圆约化抑制具有级联趋化性的觅食者-捕食者系统的 Hopf 分岔

Parabolic-elliptic reduction suppresses Hopf bifurcation in a forager-exploiter system with cascade taxis

作者: Huaizhi Cao, Jiawei Chu, Hai-Yang Jin

类别: math.AP

发布日期: 2026-07-25

arXiv ID: 2607.22986v1

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

摘要: We study a parabolic-elliptic forager-exploiter model describing the cascade interaction between two species through a shared environmental resource with a constant renewal rate. We first establish the global existence and uniform boundedness of classical solutions in arbitrary dimensions {\color{black}for large initial population data and large taxis sensitivity coefficients. We then investigate the long-time behavior of solutions.} When the resource renewal rate is small, all solutions will converge exponentially to the unique constant steady state. Furthermore, for large resource renewal rates, the parabolic-elliptic system does not undergo Hopf bifurcation, which is in sharp contrast to the corresponding fully parabolic case, where temporal oscillations arise via Hopf bifurcation.


20. 关于一般线性退化椭圆型偏微分方程组

On General Linear Degenerate Elliptic PDE Systems

作者: Nikos Katzourakis, Frederick Temple

类别: math.AP

发布日期: 2026-07-24

arXiv ID: 2607.22888v1

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

摘要: Let $Ω\Subset \mathbb{R}^{n}$ be a strictly convex bounded domain. Suppose $\mathbf{A} : \mathbb{R}^{Nn} \longrightarrow \mathbb{R}^{Nn}$, $\mathbf{B}: \mathbb{R}^{Nn} \longrightarrow \mathbb{R}^{N}$, $\mathbf{C}: \mathbb{R}^{N} \longrightarrow \mathbb{R}^{N}$ are linear maps, where $\mathbf{A}$ is symmetric and non-negative definite. Given $f \in L^2(Ω, \mathbb{R}^N)$, we consider the problem of existence of solutions $u: Ω\longrightarrow \mathbb{R}^N$ to the PDE system

$$ \left\{ \begin{array}{rl} \displaystyle\sum_{β= 1}^{N}\sum_{i, j = 1}^{n} \mathbf{A}_{αi βj}\mathrm{D}_{ij}^{2}u_β + \sum_{β= 1}^{N} \sum_{i=1}^{n} \mathbf{B}_{αβi}\mathrm{D}_{i}u_β + \sum_{β= 1}^{N} \mathbf{C}_{αβ}u_β = f_α, &\text{ in $Ω$}, \\ u = 0,\ \,& \text{ on $\partial Ω$}. \end{array} \right. $$

This is a linear \textit{degenerate elliptic} system, and it has not been considered before without the assumption of strict rank-one convexity. In general, it may not possess not even distributional solutions. By introducing some natural structural assumptions, we prove the existence of an appropriately defined unique generalised solution $u\in L^2(Ω, \mathbb{R}^N)$, satisfying additional partial regularity properties. This paper extends earlier work of the first appearing author [\textit{N. Katzourakis, On linear degenerate elliptic PDE systems with constant coefficients}, Adv.\ in Calc.Var.\ 9:3, 283-291 (2016)] to include lower-order terms.