生成时间: 2026年02月02日 04:39:14

重点关注论文数量: 4

了解领域论文数量: 2

不相关论文数量: 18

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

1. 二维泰勒-格林涡旋的不稳定性

Instability of two-dimensional Taylor-Green Vortices

作者: Gonzalo Cao-Labora, Maria Colombo, Michele Dolce, Paolo Ventura

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.23040v1

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

详细分析

中文标题: 二维泰勒-格林涡旋的不稳定性

研究对象和背景

本文研究的是二维周期域 $ \mathbb{T}^2 = \mathbb{R}^2/(2\pi\mathbb{Z}^2) $ 上不可压理想流体的Euler 方程在vorticity（涡量）形式下的动力学行为。该系统描述无粘、无源、不可压缩二维流体的演化，其数学模型为：

$$ \begin{cases} \partial_t \omega + u \cdot \nabla \omega = 0, & (x,y) \in \mathbb{T}^2,\ t \geq 0, \\ u = \nabla^\perp \psi,\quad \Delta \psi = \omega, \end{cases} $$

其中 $ \omega = \nabla^\perp \cdot u = \partial_x u_2 - \partial_y u_1 $ 是涡量，$ \psi $ 是流函数，$ \nabla^\perp = (-\partial_y,\partial_x) $，且所有函数均满足零平均条件（即属于 $ L^2_0(\mathbb{T}^2) $），以保证流函数 $ \psi $ 的泊松方程有唯一解。

该系统的稳态解之一是著名的 Taylor–Green 涡旋（又称“蜂窝流”），其解析表达式为：

$$ \omega_E(x,y) = 2\sin x \sin y,\quad \psi_E(x,y) = -\sin x \sin y,\quad U_E = \nabla^\perp \psi_E. $$

注意到 $ \omega_E = F(\psi_E) $，其中 $ F(s) = -2s $，因此该稳态具有全局涡量-流函数关系，属于 Arnold 意义下的 Hamiltonian 结构。该流场具有丰富的对称性：关于 $ x = 0,\pi $ 和 $ y = 0,\pi $ 均为偶函数；同时在 $ (x,y) \mapsto (x+\pi,y+\pi) $ 下不变；且在原点 $ (0,0) $、$ (\pi,\pi) $ 等处存在双曲型驻点（hyperbolic stagnation points），这些结构对扰动的放大机制至关重要。

研究区域为紧致环面 $ \mathbb{T}^2 $，避免了边界效应，使 Fourier 分析成为自然工具；函数空间取为 $ L^2_0(\mathbb{T}^2) $（零均值平方可积函数），并进一步考虑 Sobolev 空间 $ H^m_0(\mathbb{T}^2) $（如 $ m = -1 $ 对应速度场能量空间）。特别地，线性化算子作用于 $ L^2_0 $，但其本质谱（essential spectrum）依赖于所选空间的正则性指标 $ m $ —— 这一点在后文稳定性分析中起关键作用。

论文聚焦于该稳态的谱稳定性（spectral stability）问题：即考察线性化 Euler 算子 $ L_E $ 在 $ L^2_0 $ 上是否存在具有正实部的本征值。若存在，则称该稳态谱不稳定（spectrally unstable）；进一步，若该不稳定模对应非光滑本征函数（如仅属 $ L^2 $ 而非 $ H^1 $），则传统非线性升级方法失效，需新工具处理。

这一问题长期悬而未决。Drivas–Elgindi 在 2023 年综述中仍将 Taylor–Green 涡旋的不稳定性列为公开问题。此前仅 Zhao–Protas–Shvydkoy（2022）通过高精度数值模拟发现其存在复本征值，但缺乏严格证明。本文首次给出完全严格的谱不稳定性刻画，填补了数学流体力学中一个基础性空白。

值得注意的是，该稳态不属于经典稳定性理论适用的“单向流”（unidirectional flows），如剪切流 $ u = (U(y),0) $ 或径向涡旋；其流线呈闭合胞状结构，几何复杂，导致标准 Rayleigh–Kelvin–Tollmien 理论无法应用。此外，其线性化算子 $ L_E $ 是非正规算子（non-normal），其谱不仅含离散本征值，还包含连续谱（essential spectrum），且不稳定本征值嵌入（embedded）于连续谱中 —— 特别是在速度空间 $ L^2 $（对应 $ m = -1 $）中，其本质谱为 $ i\mathbb{R} $，而所发现的不稳定本征值 $ \lambda_\star $ 满足 $ \operatorname{Re}\lambda_\star > 0 $，却仍满足 $ \lambda_\star \in i\mathbb{R} $ 的闭包，即 $ \lambda_\star \in \sigma_{\mathrm{ess}}(L_E) $。这种嵌入不稳定性（embedded instability）极大增加了分析难度，使传统扰动理论（如 Fredholm 理论或 Weyl 序列构造）难以直接使用。

因此，本文的研究对象本质上是一个具有 Hamiltonian 结构、非正规性、嵌入谱、低正则本征模的二维周期涡旋稳态，其稳定性问题横跨偏微分方程、泛函分析、动力系统与计算数学多个领域，具有典型性与挑战性。

主要定理或主要结果

论文的核心成果由以下定理与推论构成：

Theorem 1.1（Taylor–Green 涡旋的谱不稳定性刻画）
设 $ L_E : \operatorname{Dom}(L_E) \subset L^2_0(\mathbb{T}^2) \to L^2_0(\mathbb{T}^2) $ 为线性化 Euler 算子，定义为

$$ L_E(w) = -\big\{ \psi_E, w \big\} - \big\{ \Delta^{-1} w, F'(\psi_E) \big\} = -\big\{ \psi_E, (I + 2\Delta^{-1}) w \big\}, $$


其中 $ \{f,g\} = \nabla^\perp f \cdot \nabla g $ 为 Poisson 括号。则存在唯一 $ \lambda_\star \in [0.1,0.17] - i[0.57,0.63] $，使得 $ L_E $ 在虚轴外的谱为四个单重本征值：

$$ \sigma(L_E) \setminus i\mathbb{R} = \{ \lambda_\star,\ \overline{\lambda_\star},\ -\lambda_\star,\ -\overline{\lambda_\star} \}. $$

Corollary 1.2（非线性不稳定性）
设 $ \omega_\star \in L^2_0(\mathbb{T}^2) $ 是 $ L_E $ 关于 $ \lambda_\star $ 的本征函数。则

  • 正则性：$ \omega_\star \in (W^{1,p} \cap L^q)(\mathbb{T}^2) $，对任意 $ 1 \le p < \frac{1}{1-\operatorname{Re}\lambda_\star} $，$ 1 \le q < \infty $；
  • 非线性不稳定性：存在常数 $ C_0,\beta_0,\delta_0 > 0 $，使得对任意 $ 0 < \delta < \delta_0 $，初值为 $ \omega(0,\cdot) = \omega_E + \delta \operatorname{Re}(\omega_\star) $ 的 Euler 解 $ (\omega(t),u(t)) $ 满足
    $$ \sup_{0 < t < C_0 |\ln \delta|} \| u(t) - U_E \|_{L^2} \ge \beta_0. $$

Corollary 1.3（粘性扰动下的不稳定性延续）
对足够小的 $ \nu > 0 $，粘性算子 $ L_\nu = L_E + \nu \Delta $（对应受迫二维 Navier–Stokes 线性化）的谱中存在一对共轭本征值 $ \{ \lambda_\nu, \overline{\lambda_\nu} \} \subset \sigma(L_\nu) $，满足 $ \operatorname{Re} \lambda_\nu > 0 $ 且 $ \lambda_\nu \to \lambda_\star $ 当 $ \nu \to 0^+ $。

Theorem 1.5（一般 Hamiltonian 算子的不稳定性判据）
设 $ L = JH $ 为 Hilbert 空间 $ X $ 上的 Hamiltonian 算子（$ J $ 斜自伴、$ H $ 自伴有界）。若 $ H = H_s + H_u $ 满足：

  1. $ H_s $ 正定（$ \inf \langle H_s f,f\rangle > 0 $）；
  2. $ H_u $ 有限秩，$ \dim \operatorname{Ran}(H_u) = n < \infty $；
  3. $ [H_s,H_u] = 0 $；
    定义 $ T = H_s^{1/2} J H_s^{1/2} $，对 $ \lambda \in \mathbb{C}\setminus i\mathbb{R} $，令
    $$ M_\lambda = I + H_u H_s^{-1} (T - \lambda)^{-1} T,\quad \Phi(\lambda) = \det M_\lambda. $$

    则 $ \ker(L - \lambda) \cong \ker(M_\lambda) $，从而
    $$ \lambda \in \sigma(L) \setminus i\mathbb{R} \iff \Phi(\lambda) = 0. $$

该判据将无限维谱问题完全等价约化为有限维矩阵行列式零点问题，是全文技术基石。

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

a) 核心分析学工具

本文构建了一套融合Hamiltonian 结构分析、算子分解、复分析与计算机辅助证明的新型稳定性理论框架，其核心工具包括：

  • Hamiltonian 算子的正交分解与变分结构：利用 $ L = JH $ 的结构，将自伴部分 $ H $ 分解为“稳定” $ H_s $（正定）与“不稳定” $ H_u $（有限秩负方向）两部分。此分解并非任意，而是严格依据 $ H $ 的谱性质（Fourier 乘子 $ 1 - 2/(j^2+k^2) $ 的负特征空间 $ N = \operatorname{span}\{ \sin x,\sin y,\cos x,\cos y \} $）进行，确保 $ H_u $ 精确捕获所有可能的不稳定自由度。

  • 斜自伴算子的预解式强收敛理论：关键引理 Lemma 2.4 证明了对任意闭稠定斜自伴算子 $ T $，有

    $$ \lim_{\lambda \to 0^+} (T - \lambda)^{-1} T = I - P_T, $$


    其中 $ P_T $ 是到 $ \ker T $ 的正交投影。该结论基于谱定理与控制收敛定理，将预解式极限转化为投影算子，为后续 $ \Phi(\lambda) $ 在 $ \lambda = 0 $ 处的极限分析（Proposition 2.3）提供了严格基础。

  • 复分析中的幅角原理与区间算术（Interval Arithmetic）：为定位 $ \Phi(\lambda) = 0 $ 的根，作者将复平面区域 $ [0.1,0.17] \times [-0.63,-0.57] $ 划分为细网格，并在每个矩形子区域上用区间算术严格计算 $ \Phi $ 的值域（即包含所有可能取值的区间）。若某子区域上 $ \Phi $ 的值域不包含零，则根必不在其中；若覆盖整个区域的子区域均被排除，则矛盾。最终在某子区域确认 $ 0 \in \operatorname{Range}(\Phi) $，结合 $ \Phi $ 的解析性与幅角原理，严格证明零点存在。该过程在 SageMath 中实现，代码开源，确保数学严谨性。

  • 对称性约化与不变子空间分解：利用 Taylor–Green 涡旋的 $ \pi $-周期性与奇偶性，将 $ L^2_0 $ 分解为八个 $ L_E $-不变子空间：

    $$ L^2_0 = \bigoplus_{\varepsilon_x,\varepsilon_y \in \{\text{odd},\text{even}\}} L^2_{\varepsilon_x \varepsilon_y}, $$


    并进一步按 Fourier 模式 $ \sin(jx)\sin(ky) $ 的 $ j+k $ 奇偶性细分。在此框架下，$ H_u $ 的负方向仅出现在 $ \text{Even-Even} $ 子空间，且仅贡献两个负方向，从而将问题降至 $ n = 2 $ 的 $ \Phi(\lambda) $ 计算；再利用复共轭对称性，最终归结为单变量复函数 $ \Phi(\lambda) $ 的零点搜索。

b) 技术性难题的处理技巧

  • 嵌入本征值的处理：针对 $ \lambda_\star \in \sigma_{\mathrm{ess}}(L_E) $ 的困难，作者摒弃传统分离谱/连续谱的方法，转而采用 “有限秩扰动 + 预解式解析延拓” 思路。通过 Theorem 1.5，将 $ \lambda $ 视为参数，构造全纯函数 $ \Phi(\lambda) $，其零点即为离散谱点，无论是否嵌入。此法绕开了对本质谱精细结构的刻画，是处理嵌入不稳定性的一大突破。

  • 低正则本征函数的容错性设计：数值实验表明 $ \omega_\star $ 在双曲驻点处不光滑（非 $ H^1 $）。传统 Galerkin 方法要求本征函数足够光滑以保证截断误差可控，而本文的 $ \Phi(\lambda) $ 定义于 Fourier 基上，仅需计算 $ H_u $ 在有限维空间上的作用，完全规避了对 $ \omega_\star $ 全局正则性的依赖，体现了方法对奇异性的天然鲁棒性。

  • 计算机辅助证明的数学可信度保障：不同于黑箱数值模拟，本文采用区间算术，所有浮点运算均以区间（上下界）形式执行，并严格追踪舍入误差。例如，计算 $ a + b $ 时，程序返回区间 $ [a_L + b_L, a_R + b_R] $，确保真值必然落入其中。这使得“$ \Phi $ 在某区域变号”、“$ \Phi $ 的模小于某正数”等断言具有可验证的数学确定性，将计算从启发式提升至证明级别。

与之前工作的比较

a) 突破与改进

本文相较既有文献实现了三重根本性突破：

  • 适用范围的根本拓展：此前所有严格谱不稳定性判据（如 Lin–Zeng 2017 的 Pontryagin 空间理论、Friedlander–Strauss–Vishik 1997 的能量泛函方法）均要求稳态为单向流（shear flow）或径向对称涡旋，其流线为直线或同心圆，从而可分离变量并利用 Sturm–Liouville 理论。而 Taylor–Green 涡旋的流线为复杂闭合曲线，传统方法完全失效。本文提出的 $ \Phi(\lambda) = \det M_\lambda $ 判据不依赖任何几何假设，仅需 $ H $ 可分解为正定+有限秩，故可推广至任意具有类似谱结构的稳态（如广义 $ (m,n) $-Taylor–Green 涡旋，见 Section 5），标志着从“特例分析”到“普适框架”的范式跃迁。

  • 对嵌入不稳定性（Embedded Instability）的首次严格处理：以往工作（如 Bardos–Guo–Strauss 2002）要求不稳定本征值脱离本质谱（即 $ \operatorname{Re}\lambda > \mu $，$ \mu $ 为 Lyapunov 指数），以启用标准扰动理论。但本文 $ \lambda_\star $ 满足 $ \operatorname{Re}\lambda_\star \approx 0.13 < \mu = 1 $，严格处于本质谱内部。Corollary 1.2 直接引用 Lin（2004）的深刻结果，证明即使嵌入，只要存在 $ L^2 $ 本征函数，即可推出非线性不稳定性，打破了“脱离本质谱”这一长期桎梏。

  • 计算严格性与数学深度的统一：此前计算机辅助证明多用于验证已知现象（如 Schrödinger 方程基态唯一性），或依赖高度光滑性假设。本文首次将区间算术与 Hamiltonian 谱理论深度耦合，不仅定位了零点，更给出了零点所在区域的显式、可验证的数学描述（矩形框），且整个证明链（从 PDE 到矩阵到复分析）均为解析推导，计算仅作为最后一步的“存在性确认”。这为未来研究复杂流场稳定性树立了新的严格性标杆。

综上，本文不仅是对一个经典问题的解答，更是提供了一套普适、鲁棒、严格的非线性波与流体稳定性分析新范式，其影响将远超 Taylor–Green 涡旋本身。

2. 受限四阶Prandtl方程的局部适定性与解的爆破

Local Well-posedness and Blow-up for the Restricted Fourth-Order Prandtl Equation

作者: Ik Hyun Choi

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22940v1

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

详细分析

中文标题: 受限四阶Prandtl方程的局部适定性与解的爆破

1. 研究对象和背景

本文研究的核心对象是如下定义在半直线 $\mathbb{R}_+ = [0,\infty)$ 上、带夹紧边界条件（clamped boundary conditions）的一维受限四阶Prandtl型方程：

$$ a_t = -a_{yyyy} + a^2 - a_y \int_0^y a(t,\eta)\,d\eta, \quad (t,y) \in (0,T) \times \mathbb{R}_+, $$

并满足边界条件：

$$ a(t,0) = 0, \quad a_y(t,0) = 0, $$

以及当 $y \to \infty$ 时的适当衰减条件（如 $a(t,y) \to 0$ 或匹配外部流）。该方程记为 (5) 式，是从二维四阶Prandtl系统通过线性于 $x$ 的试探解假设（linear-in-$x$ ansatz）约化而来。

原始二维系统源于对经典 Prandtl 边界层方程的高阶耗散正则化，其形式为：

$$ \begin{cases} u_t + u u_x + v u_y = -u_{yyyy}, & (t,x,y) \in (0,T) \times \mathbb{R} \times \mathbb{R}_+, \\ u_x + v_y = 0, & (t,x,y) \in (0,T) \times \mathbb{R} \times \mathbb{R}_+, \end{cases} $$

配以夹紧边界条件：

$$ u(t,x,0) = 0, \quad u_y(t,x,0) = 0, \quad v(t,x,0) = 0. $$

该系统描述了在强剪切流中受四阶粘性效应主导的边界层演化，物理上可对应于某些非牛顿流体或高阶梯度扩散模型。值得注意的是，此处的四阶耗散项 $-u_{yyyy}$ 不同于经典的二阶 Navier–Stokes 或 Prandtl 型耗散，它导致更强的色散-耗散耦合，也使得相应的半群理论更为复杂。

作者采用形如 $u(t,x,y) = -x a(t,y)$ 的自相似结构试探解（ansatz），该假设将二维系统降维为关于 $a = a(t,y)$ 的一维非局部偏微分方程。由不可压缩性 $u_x + v_y = 0$ 可得 $v(t,x,y) = \int_0^y a(t,\eta)\,d\eta$，代入动量方程后即导出上述受限方程。因此，该一维模型虽为简化，却完全保留了原系统的奇异性本质：非线性项包含一个非局部积分耦合 $a_y \int_0^y a\,d\eta$，这使得标准能量方法失效，且无法直接应用经典抛物型理论。

研究区域为半直线 $\mathbb{R}_+$，而非全空间 $\mathbb{R}$，带来显著技术挑战：边界反射效应、核函数在 $x=0$ 处的奇性控制、以及Sobolev 嵌入在无界域上的失效风险。边界条件为二阶夹紧型（即函数及其一阶导数在边界处为零），对应于双调和算子 $\partial_y^4$ 的标准 Dirichlet 型边界条件，确保算子在 $L^2(\mathbb{R}_+)$ 上是正定自伴的。

函数空间选取亦具特殊性：文中主要工作于 $W^{m,1}(\mathbb{R}_+)$ 空间（即 $L^1$ 基础上的 Sobolev 空间），而非更常见的 $H^s$ 或 $L^2$-based 空间。这是因为非线性项含卷积型积分，其估计天然依赖于 $L^1$-型核范数；而 Duhamel 公式闭合要求半线性项在 $L^1$ 意义下一致有界，故必须建立半直线双调和热核 $K(t,x,y)$ 及其各阶导数的 uniform $L^1$ 估计：

$$ \int_0^\infty |\partial_x^m K(t,x,y)|\,dy \le C_m t^{-m/4}, \quad \forall t>0,\, x\ge 0. $$

这一估计远比通常的 $L^2$-核估计困难，因 $L^1$ 范数不具有 Plancherel 性质，且边界破坏傅里叶对称性。此外，初始数据需满足相容性条件（compatibility conditions），例如若 $a_0 \in W^{4,1}(\mathbb{R}_+)$，则须有 $a_0(0) = a_0'(0) = 0$，以保证解在 $t=0$ 处连续满足边界条件。

综上，该问题本质上是：在非紧、带边界的非均匀区域上，针对一个含非局部非线性项的高阶退化抛物型方程，建立其经典解的局部存在唯一性，并进一步揭示其有限时间爆破现象——这既是数学上对高阶边界层理论的深化，也对理解高雷诺数下湍流前兆的奇点形成机制具有潜在意义。

2. 主要定理或主要结果

论文确立了两个核心结论，分别对应局部适定性与爆破现象：

Theorem 1（热核表示定理）：
设 $f \in C(\mathbb{R}_+) \cap L^1(\mathbb{R}_+)$，考虑夹紧双调和热方程：

$$ \begin{cases} u_t + u_{xxxx} = 0, & t > 0,\, x > 0, \\ u(t,0) = 0,\; u_x(t,0) = 0, & t > 0, \\ u(0,x) = f(x), & x > 0, \end{cases} $$

并要求解满足 $\lim_{x\to\infty} \partial_x^m u(t,x) = 0$ 对所有 $m \ge 0$ 成立。定义辅助函数 $\Phi(r) = e^{-r} + \sin r - \cos r$，及核：

$$ K(t,x,y) = \frac{1}{\pi} \int_0^\infty e^{-k^4 t}\, \Phi(kx)\, \Phi(ky)\, dk. $$

则函数

$$ u(t,x) = \int_0^\infty K(t,x,y)\, f(y)\, dy $$

是上述初边值问题的古典解。

Lemma 1（一致 $L^1$ 核估计）：
对任意整数 $m \ge 0$，存在常数 $C_m > 0$，独立于 $t > 0$ 和 $x \ge 0$，使得

$$ \int_0^\infty |\partial_x^m K(t,x,y)|\, dy \le C_m\, t^{-m/4}. $$

该引理是后续所有分析的基础。

Main Theorem（局部适定性与爆破）：
设初始数据 $a_0 \in W^{4,1}(\mathbb{R}_+)$ 满足相容性条件 $a_0(0) = a_0'(0) = 0$。则存在 $T > 0$，使得受限四阶 Prandtl 方程

$$ a_t = -a_{yyyy} + a^2 - a_y \int_0^y a(t,\eta)\,d\eta, $$

在 $[0,T) \times \mathbb{R}_+$ 上存在唯一的古典解 $a \in C([0,T); W^{4,1}(\mathbb{R}_+)) \cap C^1((0,T); L^1(\mathbb{R}_+))$，且满足边界条件 $a(t,0) = a_y(t,0) = 0$。

进一步，存在初值 $a_0 \ge 0$，$a_0 \not\equiv 0$，使得对应解在有限时间 $T^* < \infty$ 内发生**$L^\infty$-爆破**，即

$$ \limsup_{t \to T^*-} \|a(t,\cdot)\|_{L^\infty(\mathbb{R}_+)} = +\infty. $$

该爆破由非线性项 $a^2$ 主导，通过构造合适的下解（sub-solution）并结合比较原理实现。

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

a) 核心分析学工具

本文采用Duhamel 原理 + 不动点论证（contraction mapping principle）构建局部解。具体而言，将方程写为积分形式：

$$ a(t) = e^{-t\partial_y^4} a_0 + \int_0^t e^{-(t-s)\partial_y^4} \Big[ a^2(s) - a_y(s) \int_0^y a(s,\eta)\,d\eta \Big]\, ds, $$

其中 $e^{-t\partial_y^4}$ 是夹紧双调和半群，由核 $K$ 实现。为在 Banach 空间 $X_T = C([0,T]; W^{4,1}(\mathbb{R}_+))$ 中应用压缩映射，需验证：

  • 半群作用保持正则性：$\|e^{-t\partial_y^4} f\|_{W^{4,1}} \le C \|f\|_{W^{4,1}}$；
  • 非线性项映射 $X_T \to X_T$ 且 Lipschitz 连续；
  • 积分算子在小时间尺度下为压缩。

关键支撑是 uniform $L^1$ 核估计（Lemma 1），它直接推出：

$$ \|\partial_y^m e^{-t\partial_y^4} f\|_{L^1} \le C_m t^{-m/4} \|f\|_{L^1}, $$

从而控制高阶导数的增长。此外，作者利用分部积分导出的核新表示（alternative representation via integration by parts），绕过直接估计振荡积分的困难，证明半群在相容初值下保持 $W^{m,1}$ 正则性。

另一重要工具是Riemann–Lebesgue 引理与支配收敛定理，用于验证解在无穷远处趋于零（即 $\lim_{y\to\infty} a(t,y) = 0$），这对保证物理合理性及能量守恒至关重要。

b) 技术性难题的处理技巧

  • 非局部非线性项的控制：项 $a_y \int_0^y a\,d\eta$ 既含导数又含积分，不能直接用 Hölder 不等式。作者将其拆分为乘积形式，并利用 $L^1$–$L^\infty$ 对偶性与核的 $L^1$ 有界性进行估计：

    $$ \left\| \int_0^t e^{-(t-s)\partial_y^4} \big(a_y \cdot (Ia)\big)(s)\,ds \right\|_{L^1} \lesssim \int_0^t (t-s)^{-1/4} \|a_y(s)\|_{L^1} \|Ia(s)\|_{L^\infty}\,ds. $$

  • 边界奇性处理：由于 $\Phi(kx)$ 在 $x=0$ 处展开为 $\Phi(0)=0$, $\Phi'(0)=0$，核 $K(t,x,y)$ 自动满足夹紧条件。作者通过精细分析 $\Phi^{(m)}(kx)$ 在 $x=0$ 附近的渐近行为，结合变量替换 $s = k t^{1/4}$，分离出 $t$-尺度与空间尺度，实现对 $x=0$ 附近一致估计。

  • 爆破构造：摒弃抽象泛函方法，转而构造显式下解 $\underline{a}(t,y) = \lambda(t) \phi(y)$，其中 $\phi \ge 0$ 满足 $\phi(0)=\phi'(0)=0$ 且 $\int_0^\infty \phi < \infty$。代入方程后导出关于 $\lambda(t)$ 的 ODE：

    $$ \lambda'(t) \ge -C \lambda(t) + \alpha \lambda(t)^2, $$

    当 $\lambda(0)$ 足够大时，该 Riccati 型 ODE 在有限时间爆破，从而 $\|a(t)\|_{L^\infty} \ge \lambda(t) \to \infty$。

4. 与之前工作的比较

a) 突破与改进

本文首次建立了四阶 Prandtl 型方程在半直线上的 $L^1$-框架下的局部适定性，此前工作多集中于：

  • 经典二阶 Prandtl 方程的 $H^s$ 理论（如 Alexandre–Wang–Xu, 2015）；
  • 全空间上四阶色散方程（如 biharmonic NLS）的 $L^2$ 理论；
  • 半直线上线性四阶方程的 Fokas 统一变换解法（如 Fokas–Pelloni, 2005），但未给出 $L^1$ 核估计。

本文关键改进在于：

  • 放宽正则性要求：以往高阶抛物问题常需 $H^s$ 初值（$s > 5/2$），而本文仅需 $W^{4,1}$，且允许初值无紧支集或指数衰减；
  • 突破非局部障碍：首次将含 $\int_0^y a$ 的非局部项纳入严格适定性框架，此前类似工作（如 Liu–Yang, 2020）仅处理局部非线性；
  • 爆破机制揭示：不同于二阶 Prandtl 中由逆向流动引发的“分离型爆破”，本文爆破纯由非线性自增强驱动，体现高阶耗散不足以抑制 $a^2$ 增长，具新物理含义。

综上，该文不仅拓展了高阶边界层理论的数学基础，也为非局部抛物模型的分析提供了可推广的技术范式。

3. 密度依赖黏性及允许真空情形下二维非均匀Navier-Stokes方程Cauchy问题强解的整体适定性

Global Well-posedness of Strong Solutions to the Cauchy Problem of 2D Nonhomogeneous Navier-Stokes Equations with Density-Dependent Viscosity and Vacuum

作者: Bing Yuan, Rong Zhang, Peng Zhou

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22877v1

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

详细分析

中文标题: 密度依赖黏性及允许真空情形下二维非均匀Navier-Stokes方程Cauchy问题强解的整体适定性

研究对象和背景

本文研究的是二维非均匀、不可压缩、密度依赖黏性的Navier-Stokes方程组的Cauchy问题（即全空间 $\mathbb{R}^2$ 上的初值问题），其核心方程为：

$$ \begin{cases} \partial_t \rho + \mathrm{div}(\rho u) = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) - \nabla^\perp\big( \mu(\rho)\omega \big) + \nabla P = 0, \\ \mathrm{div}\, u = 0, \end{cases} \quad (x,t) \in \mathbb{R}^2 \times (0,T), $$

其中：

  • $\rho = \rho(x,t) \geq 0$ 是流体密度，允许在某些区域（甚至大范围）退化为零，即允许真空；
  • $u = (u_1,u_2)(x,t)$ 是速度场，满足不可压缩性约束 $\mathrm{div}\,u = 0$；
  • $P = P(x,t)$ 是压力；
  • $\omega = \nabla^\perp \cdot u \triangleq \partial_2 u_1 - \partial_1 u_2$ 是涡量（vorticity）；
  • $\nabla^\perp f \triangleq (\partial_2 f, -\partial_1 f)$ 是正交梯度算子；
  • 黏性系数 $\mu(\rho)$ 是密度的正则函数，满足
    $$ \mu \in C^1[0,\infty), \quad \mu(\rho) > 0 \quad \text{for all } \rho \geq 0. $$

该模型是对经典非均匀Navier-Stokes方程的一种修正形式。值得注意的是，当 $\mu(\rho) \equiv \mu$ 为常数时，由于 $\mathrm{div}\,u = 0$，有恒等式

$$ \nabla^\perp\big( \mu \omega \big) = \mu \Delta u, $$


因此系统退化为标准的非均匀不可压缩Navier-Stokes方程。但当 $\mu = \mu(\rho)$ 时，扩散项不再是线性拉普拉斯型，而是与密度强耦合的非线性算子，导致强非线性耦合效应，显著增加分析难度。

论文考虑的是Cauchy问题，即定义在全平面 $\mathbb{R}^2$ 上，无边界条件，仅要求解在无穷远处具有特定渐近行为：

$$ (\rho, u)(x,t) \to (\rho_\infty, 0) \quad \text{as } |x| \to \infty, $$


其中远场密度 $\rho_\infty$ 允许取 $\rho_\infty = 0$（真空情形）或 $\rho_\infty > 0$（非真空情形） —— 这是本文的关键设定之一，区别于多数前期工作要求 $\inf \rho_0 > 0$。

初始数据满足：

  • $\rho_0 \geq 0$，可在某开集上恒为零（即允许大区域真空）；
  • 初始动量 $\rho_0 u_0 \in L^2(\mathbb{R}^2)$；
  • 初始速度 $u_0$ 属于Sobolev空间 $D^{1,2}(\mathbb{R}^2) \triangleq \{ v \in L^1_{\mathrm{loc}} : \nabla v \in L^2 \}$，且满足 $\mathrm{div}\,u_0 = 0$（即属于 $D^{1,2}_{0,\sigma}$）；
  • 在真空情形（$\rho_\infty = 0$）下，还要求加权正则性：$\sqrt{\rho_0}\, \bar{x}^a \in L^1 \cap H^1 \cap W^{1,q}$，其中权重 $\bar{x} \triangleq (e+|x|^2)^{1/2}\log^{1+\eta_0}(e+|x|^2)$，$a \in (1,2), q > 2, \eta_0 > 0$；
  • 在非真空情形（$\rho_\infty > 0$）下，要求 $\rho_0 - \rho_\infty \in L^r \cap L^\infty \cap D^{1,q}$，$r \in (1,\infty), q > 2$。

函数空间选取具有深刻物理与数学意义：

  • $D^{1,2}$ 是自然的能量空间（对应有限动能梯度），在 $\mathbb{R}^2$ 中不嵌入 $L^\infty$，故需借助加权范数控制低频衰减；
  • 权重 $\bar{x}^a$ 的引入是为了补偿 $\mathbb{R}^2$ 中 $L^1$ 型估计的缺失，并确保密度梯度 $\nabla\rho$ 的 $L^q$ 可控性（$q>2$）；
  • $W^{1,q}$ 正则性对后续应用 Beale–Kato–Majda 型不等式至关重要；
  • 所有估计均需独立于 $\inf \rho_0$，即不假设密度有正下界，这是处理真空的核心难点。

该问题的物理背景源于多相流、悬浮液、含气泡液体等密度剧烈变化的复杂流体系统，其中黏性随局部密度变化而变化，例如聚合物溶液、血液流动建模等。数学上，它属于强耦合、拟线性、退化抛物-双曲混合型偏微分方程组，其全局适定性长期悬而未决，尤其在二维真空情形下，此前尚无无需小初值假设的全局强解存在性结果。

主要定理或主要结果

本文建立两个核心定理，分别覆盖真空与非真空两类远场状态：

定理 1.1（真空情形，$\rho_\infty = 0$）

设 $\rho_\infty = 0$，初始数据 $(\rho_0, u_0)$ 满足：存在 $a \in (1,2), q > 2$，使得

$$ \sqrt{\rho_0}\,\bar{x}^a \in L^1 \cap H^1 \cap W^{1,q}, \quad u_0 \in D^{1,2}_{0,\sigma}, \quad \rho_0 u_0 \in L^2, $$


其中 $\bar{x} \triangleq (e+|x|^2)^{1/2}\log^{1+\eta_0}(e+|x|^2),\, \eta_0 > 0$。
则 Cauchy 问题存在唯一全局强解 $(\rho, u, P)$，满足对任意 $T > 0$：

  • $\rho \in C([0,T]; L^1 \cap H^1 \cap W^{1,q})$，且 $\rho \bar{x}^a \in L^\infty(0,T; L^1 \cap H^1 \cap W^{1,q})$；
  • $\sqrt{\rho} u, \nabla u, \bar{x}^{-1} u, \sqrt{t}\,\sqrt{\rho} u_t, \sqrt{t}\,\nabla P, \sqrt{t}\,\nabla^2 u \in L^\infty(0,T; L^2)$；
  • $\nabla u \in L^2(0,T; H^1) \cap L^{(q+1)/q}(0,T; W^{1,q})$；
  • $\nabla P \in L^2(0,T; L^2) \cap L^{(q+1)/q}(0,T; L^q)$；
  • $\sqrt{t}\,\nabla u \in L^2(0,T; W^{1,q}),\ \sqrt{t}\,\nabla P \in L^2(0,T; L^q)$；
  • $\sqrt{\rho} u_t, \sqrt{t}\,\nabla u_t, \sqrt{t}\,\bar{x}^{-1} u \in L^2(\mathbb{R}^2 \times (0,T))$；
  • 存在 $N_1 > 0$ 使得 $\inf_{0 \le t \le T} \int_{B_{N_1}} \rho(x,t)\,dx \ge \frac{1}{4}$（即质量不完全逃逸）；
  • 大时间衰减率：对任意 $p \in [2,\infty),\, t \ge 1$，
    $$ t^{\frac{p-1}{p}} |\nabla u(\cdot,t)|_{L^p}
    • t \left| \nabla\big( \mu(\rho)\omega \big)(\cdot,t) \right|_{L^2}
    • t |\nabla P(\cdot,t)|_{L^2} \le C, $$
      其中常数 $C$ 仅依赖于 $p, \mu, \bar{\mu}', \|\rho_0\|_{L^2 \cap L^\infty}, \|\sqrt{\rho_0} u_0\|_{L^2}, \|\nabla u_0\|_{L^2}$。

定理 1.2（非真空情形，$\rho_\infty > 0$）

设 $\rho_\infty > 0$，初始数据满足：存在 $r \in (1,\infty), q > 2$，使得

$$ \rho_0 - \rho_\infty \in L^r \cap L^\infty \cap D^{1,q}, \quad u_0 \in D^{1,2}_{0,\sigma}, \quad \rho_0 u_0 \in L^2. $$


则 Cauchy 问题存在唯一全局强解 $(\rho, u, P)$，满足类似正则性与衰减估计（详见原文 (1.12)–(1.13)），且常数 $C$ 依赖于 $p,r,\rho_\infty,\mu,\bar{\mu}'$, $\|\rho_0-\rho_\infty\|_{L^r \cap L^\infty}, \|\sqrt{\rho_0} u_0\|_{L^2}, \|\nabla u_0\|_{L^2}$。

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

a) 核心分析学工具

  1. 局部适定性引理（Lemma 2.1）：基于标准迭代与紧性论证，给出短时强解存在唯一性，为全局延拓提供起点。

  2. 先验估计框架：全部证明围绕构造与 $\inf \rho_0$ 无关的高阶先验估计展开，核心目标是控制：

    • $\sup_{t\in[0,T]} \|\nabla \rho\|_{L^q}$（$q>2$），此为突破真空障碍的关键；
    • $\sup_{t\in[0,T]} \|\nabla u\|_{L^\infty}$，以控制非线性对流项；
    • $\int_0^T \|\nabla u\|_{L^\infty} dt < \infty$，保证解不爆破。

  3. 材料导数（Material Derivative）技巧：定义 $f^\cdot \triangleq f_t + u \cdot \nabla f$，将动量方程改写为含 $u^\cdot$ 的形式，避免直接处理 $\mathrm{div}(\rho u \otimes u)$ 的强非线性，使能量估计更干净。

  4. 涡量-压力耦合结构利用：关键观察在于，由 $\mathrm{div}\,u = 0$，可将动量方程中耗散项 $\nabla^\perp(\mu(\rho)\omega)$ 与压力梯度分离，并通过椭圆理论获得：

    $$ \left\| \nabla\big( \mu(\rho)\omega \big) \right\|_{L^2} + \|\nabla P\|_{L^2} \le C \|\sqrt{\rho}\, u^\cdot\|_{L^2}, $$

    此不等式（见 (3.6)）成为连接密度梯度与速度时间导数的桥梁。

  5. Beale–Kato–Majda 型不等式（Lemma 2.7）：对 $q > 2$，

    $$ \|\nabla u\|_{L^\infty} \le C \left( \|\mathrm{div}\,u\|_{L^\infty} + \|\nabla^\perp \cdot u\|_{L^\infty} \right) \log\big(e + \|\nabla^2 u\|_{L^q}\big) + C \|\nabla u\|_{L^2} + C. $$

    结合 $\mathrm{div}\,u = 0$ 及 $\omega = \nabla^\perp \cdot u$，转化为对 $\|\mu(\rho)\omega\|_{L^2}$ 和 $\|\nabla(\mu(\rho)\omega)\|_{L^q}$ 的控制，最终归结为 $\|\sqrt{\rho}\,u^\cdot\|_{L^q}$ 的估计。

  6. 加权 Sobolev 不等式（Lemmas 2.4–2.5）：用于处理 $\mathbb{R}^2$ 中低维空间的紧性缺失，特别是建立加权密度控制：

    $$ \|\rho v\|_{L^r} \le C \left( \|\sqrt{\rho} v\|_{L^2(B_{N_1})} + \|\nabla v\|_{L^2} \right), $$

    对控制 $\nabla \rho$ 的 $L^q$ 范数至关重要。

  7. 对数 Gronwall 不等式：由 Beale–Kato–Majda 型估计导出形如

    $$ \frac{d}{dt} \|\nabla \rho\|_{L^q} \le C \|\nabla \rho\|_{L^q} \left( 1 + \log\big(e + \|\nabla \rho\|_{L^q}\big) \right) $$

    的微分不等式，解得 $\|\nabla \rho\|_{L^q} \le C(T)$，从而完成关键控制。

b) 技术性难题解决技巧

  • 真空处理：放弃传统 $L^\infty$ 下界假设，转而用加权 $L^1$ 范数 $\int \rho \bar{x}^a dx$ 控制质量分布，结合 Gagliardo–Nirenberg 插值与 Poincaré 型不等式（Lemma 2.5），实现 $\nabla \rho$ 的 $L^q$ 有界性。

  • 密度依赖黏性耦合：利用系统特殊结构 $\nabla^\perp(\mu(\rho)\omega) = \mathrm{div}(2\mu(\rho)D(u)) - 2\nabla u \cdot \nabla \mu(\rho)$，但本文仅取散度自由部分，回避 $\nabla u \cdot \nabla \mu(\rho)$ 项，大幅简化分析。

  • 压力估计：不直接估计 $P$，而是通过 $\nabla P$ 的椭圆方程 $\Delta P = \mathrm{div}\,F$（其中 $F$ 含 $u^\cdot$）及 Hardy–BMO 对偶性（Lemma 2.6），得到 $\|\nabla P\|_{L^2} \le C \|\sqrt{\rho}\,u^\cdot\|_{L^2}$。

  • 时间衰减：通过乘以权重 $t^i$（$i=1,2$）并应用 Grönwall，获得 $\sqrt{t}\,\nabla u \in L^\infty_t L^2_x$ 等加权估计，再结合插值与 Sobolev 嵌入推出代数衰减率。

与之前工作的比较

a) 突破与改进

本文实现了三项根本性突破：

  1. 首次在二维真空情形下取消初值小性假设：此前所有关于密度依赖黏性 Navier-Stokes 方程的全局强解结果（如 Abidi–Zhang [3,4]、Huang–Wang [17]、He–Li–Lü [13]）均要求 $\|u_0\|_{L^2}$ 或 $\|\nabla \mu(\rho_0)\|_{L^q}$ 充分小。本文完全去除任何小初值条件，允许 $u_0$ 和 $\rho_0$ 为任意大的 $L^2$-可积函数，且 $\rho_0$ 可大面积为零。

  2. 推广至全空间真空 Cauchy 问题：先前真空结果集中于有界区域（如 Huang–Wang [17]）或一维情形。本文是首个针对 $\mathbb{R}^2$ 上允许真空的密度依赖黏性系统的全局强解存在性定理，解决了文献中明确提出的开放问题（见 Introduction 第二段末）。

  3. 更广的黏性函数类：允许 $\mu(\rho) = 1 + \rho^\beta$（$\beta \ge 0$），而压缩情形中全局大解仅对 $\beta > 4/3$ 成立（[11,14]）。这表明不可压缩模型对黏性增长更宽容，物理适用性更强。

综上，本文不仅填补了二维真空情形下的理论空白，更将全局适定性理论推进至“任意大的、允许真空的、密度依赖黏性”这一最一般设定，标志着该方向的重大进展。

4. 随机欧拉方程的概率强解

Probabilistically Strong Solutions to Stochastic Euler Equations

作者: Benjamin Gess, Robert Lasarzik

类别: math.AP

发布日期: 2026-01-29

arXiv ID: 2601.22073v1

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

详细分析

中文标题: 随机欧拉方程的概率强解

研究对象和背景

本文研究的核心对象是三维不可压随机欧拉方程（stochastic incompressible Euler equations）及其粘性逼近——随机不可压纳维–斯托克斯方程（stochastic incompressible Navier–Stokes equations），二者均受加性噪声与输运噪声（transport noise）共同驱动。这些方程描述了在随机外力扰动下、满足质量守恒（即不可压性）的流体运动，是经典流体力学中确定性欧拉/纳维–斯托克斯系统在随机分析框架下的自然推广。

论文在第2.1节明确给出了所研究的统一形式的随机偏微分方程组：

$$ du + \big((u\cdot\nabla)u - \nu\Delta u\big)\,dt + d\nabla p = \sigma_1\,dW + (\sigma_2\cdot\nabla)u \circ dW, \quad \text{in } (0,T)\times D, $$

$$ \nabla\cdot u = 0, \quad \text{in } (0,T)\times D, $$

$$ \nu(I - n\otimes n)u = 0,\quad u\cdot n = 0, \quad \text{on } (0,T)\times \partial D, $$

$$ u(0,x) = u_0(x), \quad \text{in } D. $$

其中：

  • $u(\omega,t,x) \in \mathbb{R}^d$ 是速度场，$p(\omega,t,x)\in\mathbb{R}$ 是压力（作为不可压约束 $\nabla\cdot u = 0$ 的拉格朗日乘子）；
  • $D \subset \mathbb{R}^d$ 是一个具有 $C^2$ 边界的有界区域，维度取 $d = 2$ 或 $3$，本文重点处理三维情形（$d=3$），这是数学上最困难且物理意义最丰富的设置；
  • $\nu > 0$ 是粘性系数；当 $\nu \to 0$ 时，系统退化为随机欧拉方程；
  • $W = (W_t)_{t\geq 0}$ 是定义在可分希尔伯特空间 $U$ 上的柱状维纳过程（cylindrical Wiener process），其轨道几乎必然连续地取值于某个嵌入空间 $U_0 \supset U$；
  • $\sigma_1 \in L_2(U;L^2_\sigma(D))$ 是加性噪声的 Hilbert–Schmidt 算子；
  • $\sigma_2 \in L_2(U; H^2(D)\cap H^1_{0,\sigma}(D))$ 控制输运噪声项 $(\sigma_2\cdot\nabla)u \circ dW$，该符号 $\circ$ 表示 Stratonovich 积分，后通过 Itô 公式转化为等价的 Itô 形式：
    $$ du + P\Big[(u\cdot\nabla)u - \nu\Delta u - \frac{1}{2}(\sigma_2\cdot\nabla)P[(\sigma_2\cdot\nabla)u]\Big]\,dt = \big(\sigma_1 + P[(\sigma_2\cdot\nabla)u]\big)\,dW. $$
    这里 $P$ 是 Helmholtz 投影，将向量场正交投影到散度为零的空间 $L^2_\sigma(D)$ 上，从而消去压力项，得到无压形式的演化方程。

边界条件采用零通量条件（zero outflow）：$u\cdot n = 0$，并在 $\nu > 0$ 时进一步要求切向分量消失（即 $ (I - n\otimes n)u = 0$），这涵盖了标准的齐次 Dirichlet 边界条件 $u = 0$ 在 $\partial D$ 上的情形。

函数空间设定极为关键：初始数据 $u_0$ 属于 $L^p(\Omega; L^2_\sigma(D))$，其中 $p > 4$，即仅需 $L^2$ 初值的有限高阶矩，而不要求任何 Sobolev 正则性（如 $H^1$ 或更高）。这一低正则性假设直接对应物理中“粗糙初态”的建模需求，也是传统强解理论无法覆盖的区域。

此外，作者特别强调，他们考虑的是概率强解（probabilistically strong solution），即解过程 $u$ 是定义在给定概率空间 $(\Omega,\mathcal{F},(\mathcal{F}_t),\mathbb{P})$ 上的适应过程（adapted process），不依赖于对原始概率空间的修改。这与经典的 martingale 解（或称 probabilistically weak solution）形成鲜明对比——后者通常需借助 Skorokhod–Jakubowski 表示定理改变概率空间以获得几乎必然收敛，因而丧失了原始噪声结构的确定性关联。

因此，本文的研究背景可概括为：在最弱的初始数据假设（$L^2$）与最自然的概率设定（给定噪声、无需重构造空间）下，建立三维随机欧拉方程的存在性理论，并揭示其作为粘性极限的结构性起源。

主要定理或主要结果

本文最核心的成果由以下定理体现：

Theorem 3.7（Navier–Stokes 方程的概率强测度值解存在性）：
设 $D \subset \mathbb{R}^3$ 是 $C^2$ 有界区域，$(\Omega,\mathcal{F},(\mathcal{F}_t),\mathbb{P})$ 是完备概率空间，$\sigma_1,\sigma_2$ 满足 Assumption 3.1，且初始数据 $u_0 \in L^p(\Omega;\mathcal{F}_0;L^2_\sigma(D))$，$p > 4$。则存在一个 probabilistically strong, measure-valued, and probabilistically weak, weak solution $(\tilde{\Omega},\tilde{\mathcal{F}},(\tilde{\mathcal{F}}_t),\tilde{\mathbb{P}},\tilde{u},\tilde{W})$ 及其对应的 Young 测度 $\mu$，满足 Definition 3.6，并且

$$ u = \langle \mu, \mathrm{id} \rangle \in L^p\big(\Omega; L^\infty(0,T;L^2_\sigma(D)) \cap L^2(0,T;H^1_\sigma(D)) \cap W^{\alpha,p/2}(0,T;(H^2\cap H^1_{0,\sigma})^*)\big), $$


对任意 $\alpha < \frac{1}{2}$ 成立；同时 Young 测度 $\mu$ 属于

$$ L^s\big(\Omega; L^s(D\times(0,T); \mathcal{P}(\mathbb{R}^3))\big), \quad s = \frac{2(d+2)}{d} = \frac{10}{3}. $$

Theorem 4.5（Euler 方程的概率强解存在性及 vanishing viscosity 极限）：
在相同假设下，令 $\{u^\nu\}_{\nu>0}$ 为上述 Theorem 3.7 所构造的随机 Navier–Stokes 方程的概率强测度值解序列。则存在子列 $\nu_k \to 0$，使得 $u^{\nu_k} \to u$ 在适当拓扑下收敛，且极限 $u$ 是 probabilistically strong, dissipative measure-valued solution 到随机 Euler 方程，满足能量不等式：

$$ \frac{1}{2}\int_D \langle \mu_{t}, |\lambda|^2 \rangle\,dx + \int_s^t \int_D \mathrm{Tr}\big[ \mathcal{D}(\mu_r) \big] \,dx\,dr \leq \frac{1}{2}\int_D \langle \mu_{s}, |\lambda|^2 \rangle\,dx + \int_s^t \int_D u\cdot\sigma_1 \,dx\,dW_r + \frac{1}{2}\int_s^t \|\sigma_1\|_{L_2(U;L^2_\sigma)}^2 \,dr, $$

对几乎所有的 $0 \leq s < t \leq T$，$\mathbb{P}$-a.s. 成立，其中 $\mathcal{D}(\mu)$ 是缺陷测度（defect measure），刻画非线性项的弱极限损失。

特别地，该解 $u$ 是由粘性系统一致构造而来，且对任意 $L^2$ 初值均成立，从而首次解决了 Breit & Moyo 提出的公开问题：构造满足能量不等式的、概率强的随机 Euler 解。

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

a) 核心分析学工具

本文融合了现代偏微分方程与随机分析中的多项前沿技术：

  • Young 测度理论：这是全文最根本的工具。不同于传统紧性论证依赖 Skorokhod 表示定理强制更换概率空间，作者采用 Young 测度 $\mu_{\omega,t,x} \in \mathcal{P}(\mathbb{R}^3)$ 来参数化振荡行为（oscillations）与集中现象（concentrations）。Theorem 2.3 给出了在随机时空域 $\Omega \times (0,T) \times D$ 上 Young 测度的基本存在性与弱收敛性质：若 $\{u_n\}$ 在 $L^p$ 中有界，则存在子列使得对任意 Carathéodory 被积函数 $f$，

    $$ f(\omega,t,x,u_n(\omega,t,x)) \rightharpoonup \langle \mu_{\omega,t,x}, f(\omega,t,x,\cdot) \rangle \quad \text{weakly in } L^q. $$

    这一框架允许在固定概率空间上处理非线性项 $(u\cdot\nabla)u$ 的弱极限，避免了 martingale 解路径依赖的缺陷。

  • 能量变分解（Energy-variational solutions）：作者将 deterministic setting 中由 [9,22] 发展的能量变分框架推广至随机情形。此类解不直接满足方程的分布意义，而是通过能量耗散泛函的变分不等式来刻画：

    $$ \partial_t \mathcal{E}(u) + \mathcal{D}(u) \leq \text{noise terms}, $$

    其中 $\mathcal{E}(u) = \frac{1}{2}\|u\|_{L^2}^2$，$\mathcal{D}(u)$ 包含粘性耗散或缺陷测度。该结构天然兼容 Young 测度表示，并能绕过识别非线性极限 $\lim (u_n\cdot\nabla)u_n$ 的技术障碍。

  • 随机紧性与路径正则性引理：Lemma 2.4 建立了随机 Sobolev–Slobodeckij 空间嵌入 $L^q(\Omega;W^{\alpha,r}(0,T;Y^*)) \hookrightarrow L^q(\Omega;C([0,T];Y^*))$，前提是 $\alpha r > 1$。这保证了近似解 $u_N$ 的时间连续性（在弱拓扑下），是证明解过程属于 $C_w([0,T];L^2_\sigma)$ 的关键。

  • Itô–Ventzell 公式与随机积分估计：为处理 Stratonovich 噪声项 $(\sigma_2\cdot\nabla)u \circ dW$，作者反复使用 Itô 公式将其转化为 Itô 形式，并利用 Hilbert–Schmidt 算子范数控制噪声项的 $L^2$ 增长，例如：

    $$ \mathbb{E}\left[\sup_{t\in[0,T]}\left|\int_0^t \phi\cdot\sigma_1\,dW_s\right|^2\right] \leq C \|\phi\|_{L^2}^2 \|\sigma_1\|_{L_2(U;L^2_\sigma)}^2. $$

b) 处理非线性与奇异性技巧

  • Galerkin 截断与一致先验估计：采用 Stokes 算子特征函数构成的 Galerkin 子空间 $V_N$，构造近似解 $u_N(t) = \sum_{k=1}^N a_k^N(t)v_k$。通过测试函数 $\varphi = u_N$，获得统一的 $L^\infty_t L^2_x$ 与 $L^2_t H^1_x$ 能量估计（含噪声修正项），且该估计关于 $\nu$ 一致，为 vanishing viscosity 极限提供基础。

  • 缺陷测度建模非线性极限：对非线性对流项，定义缺陷测度 $\mathcal{D} \in L^\infty_w(\Omega; \mathcal{M}^+(D\times(0,T)))$，使得

    $$ \langle \mu_{\omega,t,x}, \lambda\otimes\lambda \rangle = u\otimes u + \mathcal{D}_{\omega,t,x}, $$

    从而将 $(u_n\cdot\nabla)u_n \rightharpoonup u\cdot\nabla u + \nabla\cdot\mathcal{D}$ 分离为规则部分与奇异部分。能量不等式中 $\mathrm{Tr}[\mathcal{D}]$ 即代表“额外耗散”。

  • Transport noise 的几何结构利用：注意到输运噪声 $(\sigma_2\cdot\nabla)u \circ dW$ 在 SALT 框架下保持环量守恒，作者借助 $H^2\cap H^1_{0,\sigma}$ 正则性假设，确保算子 $(\sigma_2\cdot\nabla)$ 作用于 $u$ 后仍处于合适 Hilbert–Schmidt 类，从而维持整个系统的闭包性。

与之前工作的比较

a) 突破性改进

本文相较既有文献实现了三重实质性突破：

  • 解的概念强度提升：此前，Breit & Moyo [4] 仅构造了 probabilistically weak, dissipative measure-valued solutions；Hofmanová–Zhu–Zhu [17] 利用凸积分法获得 probabilistically strong solutions，但其解是“wild”（不满足能量不等式）、非唯一、且不能作为 vanishing viscosity 极限实现。本文首次给出 probabilistically strong + energy-dissipative + vanishing viscosity limit 三位一体的解，填补了该领域的根本性空白。

  • 初始数据正则性大幅降低：经典强解理论（如 Glatt-Holtz & Vicol [14]）要求 $u_0 \in H^s$, $s > 5/2$；而本文仅需 $u_0 \in L^2$，且具 $p>4$ 阶矩，完全脱离 Sobolev 框架，进入更普适的测度值解范畴。

  • 噪声类型全面覆盖：以往工作多集中于加性噪声；本文将结果稳健地延拓至物理意义更深刻的 transport noise，并指出其在 SALT 框架与奇点正则化中的双重角色，拓展了随机流体模型的适用边界。

综上，本文不仅解决了长期悬而未决的公开问题，更建立了一套可推广至其他非线性 SPDE 的概率强测度值解范式，标志着随机流体力学基础理论的重要跃进。

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

1. 具有复特征的双曲型偏微分方程在傅里叶–勒贝格空间上的适定性

Hyperbolic partial differential equations with complex characteristics on Fourier Lebesgue spaces

作者: Duván Cardona, William Obeng-Denteh, Frederick Opoku

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.23138v1

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

摘要翻译

中文标题: 具有复特征的双曲型偏微分方程在傅里叶–勒贝格空间上的适定性

摘要翻译: 本文旨在建立双曲型偏微分方程（PDE）在傅里叶–勒贝格空间（Fourier Lebesgue spaces）上的适定性理论。我们考虑具有复特征的双曲型算子。由于本研究方法源于调和分析，我们系统建立了相位函数为复值的傅里叶积分算子（Fourier integral operators）在傅里叶–勒贝格空间、伯格空间（Besov spaces）以及特里贝尔–李佐尔金空间（Triebel–Lizorkin spaces）上的有界性。事实上，此类算子恰为所研究PDE问题的传播子（propagators）。在有界性方面，我们证明了若干新结果：特别地，当算子的典型关系（canonical relation）满足所谓的“空间光滑分解条件”（spatial smooth factorization condition）时，相关有界性结论成立。

2. 自相似随机不可压流场中的超扩散现象与反常正则性

Superdiffusion and anomalous regularization in self-similar random incompressible flows

作者: Scott Armstrong, Ahmed Bou-Rabee, Tuomo Kuusi

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

发布日期: 2026-01-29

arXiv ID: 2601.22142v1

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

摘要翻译

中文标题: 自相似随机不可压流场中的超扩散现象与反常正则性

摘要翻译:
我们研究一个粒子在 $\mathbb{R}^d$（$d \geq 2$）中受分子扩散与随机不可压流场对流共同作用下的长时间行为。该速度场定义为一平稳随机流函数矩阵 $\mathbf{k}$ 的散度，其对应的 Hurst 指数 $\gamma > 0$ 为正，因而所生成的随机环境具有多尺度结构且满足自相似性。在微扰区域 $\gamma \ll 1$ 下，我们证明了淬火意义下的幂律型超扩散：对典型环境实现，粒子位移方差在时间 $t$ 处以 $t^{2/(2-\gamma)}$ 的速率增长——该标度律恰与重整化群启发式分析所预言的结果一致。此外，我们还确定了主导前置因子（leading prefactor），其淬火相对误差为 $O\big(\gamma^{\frac{1}{2}} |\log \gamma|^3\big)$。证明采用基于无穷小生成元 $\nabla \cdot (\nu I_d + \mathbf{k}) \nabla$ 的 Wilson 型重整化群方案，其核心是沿尺度实施自相似归纳。我们证明：在每一粗粒化尺度 $r$ 上，粗粒化后的生成元可被一个具有效扩散系数的常系数拉普拉斯算子良好逼近，且该有效扩散系数随尺度按 $r^\gamma$ 增长。这一逼近本质上具有尺度局域性：反映环境的多重分形本质，其相对误差并不随尺度增大而衰减，而是恒定地保持在 $O\big(\gamma^{\frac{1}{2}} |\log \gamma|^2\big)$ 量级。此外，我们在淬火概率测度下证明了反常正则性：对几乎每一个漂移场的实现，相应椭圆方程的解均为 Hölder 连续，Hölder 指数为 $1 - C \gamma^{\frac{1}{2}}$，且所满足的估计关于分子扩散系数 $\nu$ 和空间尺度具有一致性。

不相关论文（基本信息）

1. 非一致椭圆变分问题在有界变差函数空间（BV）上的研究

Non-uniformly elliptic variational problems on BV

作者: Lisa Beck, Franz Gmeineder, Mathias Schäffner

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.23195v1

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

摘要: We establish $\mathrm{W}^{1,1}$-regularity and higher gradient integrability for relaxed minimizers of convex integral functionals on $\mathrm{BV}$. Unlike classical examples such as the minimal surface integrand, we only require linear growth from below but not necessarily from above. This typically comes with a non-uniformly degenerate elliptic behaviour, for which our results extend the presently available bounds from the superlinear growth case in a sharp way.

2. 具有磁势的度量图上的非线性薛定谔方程

Nonlinear Schrödinger Equation with magnetic potential on metric graphs

作者: Riccardo Adami, Nicolò Cangiotti, Ivan Gallo, David Spitzkopf

类别: math.AP, math-ph

发布日期: 2026-01-30

arXiv ID: 2601.23115v1

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

摘要: In this manuscript, we shall investigate the Nonlinear Magnetic Schrödinger Equation on noncompact metric graphs, focusing on the existence of ground states. We prove that the magnetic Hamiltonian is variationally equivalent to a non-magnetic operator with additional repulsive potentials supported on the graph’s cycles. This effective potential is strictly determined by the Aharonov-Bohm flux through the topological loops. Leveraging this reduction, we extend classical existence criteria to the magnetic setting. As a key application, we characterize the ground state structure on the tadpole graph, revealing a mass-dependent phase transition. The ground states exist for sufficiently small repulsion in an intermediate regime of masses while sufficiently strong flux prevents the formation of ground states.

3. 无穷作用范围FPUT晶格中行波解的存在性

Existence of Traveling Waves in Infinite Range FPUT Lattices

作者: Michael Herrmann, Karsten Matthies, Jan-Patrick Meyer

类别: math.AP, math.DS

发布日期: 2026-01-30

arXiv ID: 2601.23091v1

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

摘要: We prove the existence of solitary waves in a lattice where all particles interact with each other by pair-wise repulsive forces that decay with distance. The variational existence proof is based on constrained optimization and provides a one-parameter family of unimodal solutions. We also describe the asymptotic behavior of large, fast, high-energy waves.

4. 关于具有临界δ-壳相互作用的二维狄拉克算子

On two-dimensional Dirac operators with critical delta-shell interactions

作者: William Borrelli, Pietro Carimati, Davide Fermi

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

发布日期: 2026-01-30

arXiv ID: 2601.23053v1

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

摘要: We study two-dimensional Dirac operators with singular interactions of electrostatic and Lorentzscalar type, supported either on a straight line or a circle. For certain critical values of the interaction strengths, the essential spectrum of such operators comprises an isolated point lying within the mass gap. We clarify the nature of this point in both geometries. For the straight line model, this point is known to be an eigenvalue of infinite multiplicity, and we provide a detailed analysis of the corresponding eigenfunctions. By contrast, in the case of a circle, we show that the said point is not itself an eigenvalue, but rather an accumulation point of a double sequence of simple eigenvalues. In view of the high degree of symmetry of the configurations under analysis, this behavior is unexpected and our findings lead us to formulate some conjectures concerning critical singular interactions supported on generic smooth curves.

5. 四维单位球面上的预设 $T$-曲率流

Prescribed $T$-curvature flow on the four-dimensional unit ball

作者: Pak Tung Ho, Cheikh Birahim Ndiaye, Liming Sun, Heming Wang

类别: math.DG, math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22934v1

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

摘要: In this paper, we study the prescribed $T$-curvature problem on the unit ball $\mathbb{B}^4$ of $\mathbb{R} ^4$ via the $T$-curvature flow approach. By combining Ache-Chang’s inequality with the Morse-theoretic approach of Malchiodi-Struwe, we establish existence results under strong Morse-type inequalities at infinity. As a byproduct of our argument, we also prove the exponential convergence of the $T$-curvature flow on $\mathbb{B}^4$, starting from a $Q$-flat and minimal metric conformal to the standard Euclidean metric, to an extremal metric of Ache-Chang’s inequality whose explicit expression was derived by Ndiaye-Sun.

6. TV-Wasserstein梯度流解的存在性

Existence of a solution of the TV Wasserstein gradient flow

作者: Kexin Lin, Filippo Santambrogio

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22847v1

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

摘要: On the flat torus in any dimension we prove existence of a solution to the TV Wasserstein gradient flow equation, only assuming that the initial density $ρ_0$ is bounded from below and above by strictly positive constants. This solution preserves upper and lower bounds of the densities, and shows a certain decay of the BV norm (of the order of $t^{-1/3}$ for $t\to 0$ – if $ρ_0\notin BV$, otherwise the BV norm is of course bounded – and of the order of $t^{-1}$ as $t\to\infty$). This generalizes a previous result by Carlier and Poon, who only gave a full proof in one dimension of space and did not consider the case $ρ_0\notin BV$. The main tool consists in considering an approximated TV-JKO scheme which artificially imposes a lower bound on the density and allows to find a continuous-in-time solution regular enough to prove that the lower bounds of the initial datum propagates in time, and study on this approximated equation the decay of the BV norm.

7. 控制型单调平均场博弈主方程的无条件适定性

Unconditional well-posedness of the master equation for monotone mean field games of controls

作者: Joe Jackson, Alpár R. Mészáros

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

发布日期: 2026-01-30

arXiv ID: 2601.22845v1

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

摘要: We establish the first unconditional well-posedness result for the master equation associated with a general class of mean field games of controls. Our analysis covers games with displacement monotone or Lasry–Lions monotone data, as well as those with a small time horizon. By unconditional, we mean that all assumptions are imposed solely at the level of the Lagrangian and the terminal cost. In particular, we do not require any a priori regularity or structural assumptions on the additional fixed-point mappings arising from the control interactions; instead we show that these fixed-point mappings are well-behaved as a consequence of the regularity and the monotonicity of the data. Our approach is bottom-up in nature, unlike most previous results which rely on a generalized method of characteristics. In particular, we build a classical solution of the master equation by showing that the solutions of the corresponding $N$-player Nash systems are compact, in an appropriate sense, and that their subsequential limit points must be solutions to the master equation. Compactness is obtained via uniform-in-$N$ decay estimates for derivatives of the $N$-player value functions. The underlying games are driven by non-degenerate idiosyncratic Brownian noise, and our results allow for the presence of common noise with constant intensity.

8. 具混合Steklov-Robin边界条件的椭圆问题的特征值单重性

Simplicity of eigenvalues for elliptic problems with mixed Steklov-Robin boundary condition

作者: Marco Ghimenti, Anna Maria Micheletti, Angela Pistoia

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22829v1

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

摘要: This paper investigates the spectral properties of two classes of elliptic problems characterized by mixed Steklov-Robin boundary conditions. Our main objective is to prove that, for a generic domain, all the eigenvalues are simple. This result is established by employing domain perturbation techniques and analyzing the transversality of the associated operators.

9. Escobar泛函与Gagliardo-Nirenberg泛函的精确阈值：Escobar–Willmore质量、几何选择性及紧性三分律

Sharp thresholds for Escobar and Gagliardo-Nirenberg functionals: the Escobar-Willmore mass, geometric selection, and compactness trichotomy

作者: Mayukh Mukherjee, Utsab Sarkar

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

发布日期: 2026-01-30

arXiv ID: 2601.22665v1

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

摘要: We develop a unified quantitative framework for sharp threshold phenomena in boundary-critical variational problems on compact Riemannian manifolds, covering the Escobar quotient and Gagliardo-Nirenberg inequalities. Via transfer-stability-reduction, we obtain attainment-versus-bubbling alternatives, $H^1$-compactness, and finite-dimensional reductions. Geometric selection is governed by mean curvature $H_g$ and a Willmore-type anisotropy from $|\mathring{\mathrm{II}}|^2$. At hemisphere threshold $S_\ast=C^*_{\mathrm{Esc}}(\mathbb S^n_+)$ for $n\ge5$ on $H_g\equiv0$, we identify a renormalized boundary mass $\mathfrak R_g=κ_1(n)\,\mathrm{Ric}_g(ν,ν)+κ_2(n)\,\mathrm{Scal}_{g|\partial M}+κ_3(n)\,|\mathring{\mathrm{II}}|^2$, $κ_3(n)<0$, yielding one-bubble expansions and energy-only estimators. Threshold dichotomy: if the first nonvanishing coefficient among $\{ρ_n^{\mathrm{conf}}H_g,\mathfrak R_g,Θ_g\}$ is negative somewhere, then $C^*_{\mathrm{Esc}}(M,g)<S_\ast$ and sequences are precompact. At threshold, blow-up concentrates where $H_g$ is critical; on $H_g\equiv0$, stationarity forces $\mathfrak R_g(p)=\nabla_\partial\mathfrak R_g(p)=0$. If $H_g$ is Morse and $\mathfrak R_g>0$ at all critical points, no bubbling occurs. In multi-bubble regime ($n\ge5$), dynamics governed by $\mathcal W_k=\sum_{i=1}^k\mathfrak R_g(x_i)$ produce $k$-bubble critical points at levels $k^{1/(n-1)}S_\ast$. In the degenerate case we obtain conformal hemispherical rigidity. The GN track yields analogous dichotomies and resolves a question of Christianson et al.: the sharp constant with small Dirichlet windows diverges at optimal capacitary rate, relating threshold to spectral/isoperimetric invariants. Applications include entropy inequalities for fast diffusion, curvature-driven NLS ground states, and (in $n=2$) Euler characteristic recovery from GN measurements.

10. Hodge-Maxwell 系统的加权估计

Weighted estimates for Hodge-Maxwell systems

作者: Rohit Mahato, Swarnendu Sil

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22604v1

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

摘要: We establish up to the boundary regularity estimates in weighted $L^{p}$ spaces with Muckenhoupt weights $A_{p}$ for weak solutions to the Hodge systems \begin{align*} d^{\ast}\left(Adω\right) + B^{\intercal}dd^{\ast}\left(Bω\right) = λBω+ f \quad \text{ in } Ω \end{align*} with either $ν\wedge ω$ and $ν\wedge d^{\ast}\left(Bω\right)$ or $ν\lrcorner Bω$ and $ν\lrcorner Adω$ prescribed on $\partialΩ.$ As a consequence, we prove the solvability of Hodge-Maxwell systems and derive Hodge decomposition theorems in weighted Lebesgue spaces. Our proof avoids potential theory, does not rely on representation formulas and instead uses decay estimates in the spirit of `Campanato method’ to establish weighted $L^{p}$ estimates.

11. 周期性量子图上狄拉克方程的谱性质与束缚态

Spectral properties and bound states of the Dirac equation on periodic quantum graphs

作者: Zhipeng Yang, Ling Zhu

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22603v1

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

摘要: We investigate nonlinear Dirac equations on a periodic quantum graph $G$ and develop a variational approach to the existence and multiplicity of bound states. After introducing the Dirac operator on $G$ with a $\mathbb Z^{d}$-periodic potential, we describe its spectral decomposition and work in the natural energy space. Under asymptotically linear or superquadratic assumptions on the nonlinearity, we establish the required linking geometry and a Cerami-type compactness property modulo $\mathbb Z^{d}$-translations. As a consequence, we prove the existence of at least one bound state and, when the nonlinearity is even, infinitely many geometrically distinct bound states.

12. 具有截断周期势且支持缺陷态的薛定谔算子的透射与反射系数

Transmission and Reflection coefficients for Schrödinger Operators with Truncated Periodic Potentials that support defect states

作者: Joseph C. Stellman, Jeremy L. Marzuola

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22544v1

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

摘要: We consider scattering waves through truncated periodic potentials with perturbations that support localized gap eigenstates. In a small complex neighborhood around an assumed positive bound state of the model operator, we prove the existence of a distinct zero reflection state, or transmission resonance. We compare its location to a previously found scattering resonance and use the properties of solutions near these interesting points to analyze the behavior of transmission and reflection coefficients of scattering solutions near the assumed bound state. By example, we also discuss the truncated simple harmonic oscillator and compare the analysis to the crystalline case.

13. 具有一般非线性项的Hardy抛物型方程解的局部存在性与不存在性

Local existence and nonexistence of solutions to the Hardy parabolic equation with general nonlinearity

作者: Yo Tsusaka

类别: math.AP

发布日期: 2026-01-30

arXiv ID: 2601.22520v1

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

摘要: In this paper, we consider the Cauchy problem for the Hardy parabolic equation with general nonlinearity and establish the local existence and nonexistence results. Our results provide the optimal integrability conditions on initial function for the existence of a local-in-time nonnegative solution. The proof of the existence result is based on the supersolution method.

14. 各类凸图形平移子的动力学稳定性

Dynamical stability of various convex graphical translators

作者: Junyoung Park

类别: math.DG, math.AP

发布日期: 2026-01-29

arXiv ID: 2601.22368v1

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

摘要: In the first part of the paper, we prove the existence of longtime solution to mean curvature flow starting from a graph of a continuous function defined over a slab. Then, we establish dynamical stability results for various types of graphical translators to mean curvature flow, namely the grim reaper, two dimensional graphical translators, and asymptotically cylindrical translators.

15. 增强型地热系统的一般化解析传热模型：刻画裂隙相互作用并修正经典模型的乐观预测

A Generalized Analytical Heat Transfer Model for Enhanced Geothermal Systems: Capturing Fracture Interactions and Correcting Classical Optimistic Predictions

作者: Nelson Barros-Galvis or Christine Ehlig-Economides or Cristi Darley Guevara

类别: math.AP, physics.geo-ph

发布日期: 2026-01-29

arXiv ID: 2601.22316v1

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

摘要: Numerical analytical heat transfer models play a critical role in geothermal design and feasibility studies. Classical solutions, such as those proposed by Gringarten et al. 1975, rely on simplified assumptions and systematically overestimate thermal performance, which can lead to unrealistic engineering decisions. This study presents a generalized analytical model for enhanced geothermal systems that explicitly captures thermal interactions between fractures while preserving analytical tractability. The formulation is based on Greenś functions and reproduces realistic thermal behavior under conditions representative of fractured geothermal reservoirs. The resulting solution is computationally efficient and sufficiently simple to be implemented directly in standard spreadsheets, without requiring Laplace space transformations or numerical inversion algorithms. The model is validated against numerical simulations performed using CMG STARS and Volsung software, showing close agreement in temperature evolution, including the effects of interacting fractures. Compared with classical analytical approaches, the proposed model corrects optimistic bias and provides more reliable predictions of production temperature and energy recovery. These results have direct implications for geothermal feasibility studies, well design, and power forecasting, effectively bridging the gap between legacy analytical models and numerical or commercial engineering tools. Building on the analytical framework originally introduced by Gringarten et al. 1975, the proposed formulation generalizes classical heat transfer solutions to account for fracture interaction while retaining analytical simplicity and practical applicability.

16. 辛正则半群及其在时频分析与演化算子中的应用

The metaplectic semigroup and its applications to time-frequency analysis and evolution operators

作者: Gianluca Giacchi, Luigi Rodino, Davide Tramontana

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

发布日期: 2026-01-29

arXiv ID: 2601.22252v1

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

摘要: We develop a systematic analysis of the metaplectic semigroup $\mathrm{Mp}_+(d,\mathbb{C})$ associated with positive complex symplectic matrices, a notion introduced almost simultaneously and independently by Hörmander, Brunet, Kramer, and Howe, thereby extending the classical metaplectic theory beyond the unitary setting. While the existing literature has largely focused on propagators of quadratic evolution equations, for which results are typically obtained via Mehler formulas, our approach is operator-theoretic and symplectic in spirit and adapts techniques from the standard metaplectic group $\mathrm{Mp}(d,\mathbb{R})$ to a substantially broader framework that is not driven by differential problems or particular propagators. This point of view provides deeper insight into the structure of the metaplectic semigroup, and allows us to investigate its generators, polar decomposition, and intertwining relations with complex conjugation and with the Wigner distribution. We then exploit these structural results to characterize, from a metaplectic perspective, classes of time-frequency representations satisfying prescribed structural properties. Finally, we discuss further implications for parabolic equations with complex quadratic Hamiltonians, we study the boundedness of their propagators on modulation spaces, we obtain estimates in time of their operator norms. Finally, we apply our theory to the study of propagation of Wigner singularities.

17. 从几何观点看微局部极大次椭圆性：I

Microlocal maximal hypoellipticity from the geometric viewpoint: I

作者: Omar Mohsen

类别: math.AP, math.OA

发布日期: 2026-01-29

arXiv ID: 2601.22122v1

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

摘要: Given some vector fields on a smooth manifold satisfying Hörmander’s condition, we define a bi-graded pseudo-differential calculus which contains the classical pseudo-differential calculus and a pseudo-differential calculus adapted to the sub-Riemannian structure induced by the vector fields. Our approach is based on geometric constructions (resolution of singularities) together with methods from operators algebras. We develop this calculus in full generality, including Sobolev spaces, the wavefront set, and the principal symbol, etc. In particular, using this calculus, we prove that invertibility of the principal symbol implies microlocal maximal hypoellipticity. This allows us to resolve affirmatively the microlocal version of a conjecture of Helffer and Nourrigat.

18. 关于磁性双分量Hunter-Saxton系统的整体弱解

On Global Weak Solutions for the Magnetic Two-Component Hunter-Saxton System

作者: Levin Maier

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

发布日期: 2026-01-29

arXiv ID: 2601.22088v1

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

摘要: We study the magnetic two-component Hunter-Saxton system (M2HS), which was recently derived in \cite{M24} as a magnetic geodesic equation on an infinite-dimensional configuration space. While the geometric framework and the global weak flow were outlined there, the present paper provides the analytical foundations of this construction from the PDE perspective. First, we derive an explicit solution formula in Lagrangian variables via a Riccati reduction, yielding an alternative proof of the blow-up criterion together with an explicit expression for the blow-up time. Second, we rigorously construct global conservative weak solutions by developing the analytic theory of the relaxed configuration space and the associated weak magnetic geodesic flow, thereby realizing the geometric program proposed in \cite{M24}.