生成时间: 2025年10月16日 04:47:57

重点关注论文数量: 6

了解领域论文数量: 15

不相关论文数量: 16

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

1. Strong solution for polymeric fluid-structure interaction with small initial acceleration

小初始加速度下聚合物流体-结构相互作用的强解

作者: Prince Romeo Mensah

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13753v1

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

详细分析

中文标题: 小初始加速度下聚合物流体-结构相互作用的强解

1. 研究对象和背景

本文研究的是一个三维-三维-二维（3D-3D-2D）相互耦合的多相系统，即聚合物流体（solute）—溶剂（solvent）—柔性结构（structure）三态耦合系统。该系统描述了一种经典Oldroyd-B型聚合物流体（无中心质量扩散）与一个二维粘弹性壳结构之间的动态相互作用。该问题本质上是流体-结构相互作用（Fluid-Structure Interaction, FSI）问题的高维推广，其特殊性在于引入了聚合物分子的微观应力张量作为独立演化变量，从而形成一个“溶质-溶剂-结构”三重耦合的非线性偏微分方程组。

系统的核心方程组由以下四个方程构成，定义在时间区间 $ I = (0,T) $ 与空间域 $ \Omega_\eta(t) \subset \mathbb{R}^3 $ 上（其中 $ \Omega_\eta(t) $ 是由结构位移 $ \eta(t,y) $ 所定义的时变域），而结构方程则定义在二维参考边界 $ \omega \subset \mathbb{R}^2 $ 上：

$$ \operatorname{div}_x \mathbf{u} = 0, \quad \text{在 } I \times \Omega_\eta \text{ 上} $$
$$ \partial_t \mathbf{u} + (\mathbf{u} \cdot \nabla_x) \mathbf{u} = \Delta_x \mathbf{u} - \nabla_x p + \mathbf{f} + \operatorname{div}_x \mathbf{T}, \quad \text{在 } I \times \Omega_\eta \text{ 上} $$
$$ \partial_t^2 \eta - \partial_t \Delta_y \eta + \Delta_y^2 \eta = g - (\mathbf{S} \mathbf{n}) \circ \varphi_\eta \cdot \mathbf{n} \, |\det(\nabla_y \varphi_\eta)|, \quad \text{在 } I \times \omega \text{ 上} $$
$$ \partial_t \mathbf{T} + (\mathbf{u} \cdot \nabla_x) \mathbf{T} = (\nabla_x \mathbf{u}) \mathbf{T} + \mathbf{T} (\nabla_x \mathbf{u})^\top + 2 \mathbf{D}(\nabla_x \mathbf{u}) - 2 \mathbf{T}, \quad \text{在 } I \times \Omega_\eta \text{ 上} $$

其中：

  • $ \mathbf{u}: I \times \Omega_\eta \to \mathbb{R}^3 $ 是溶剂（牛顿流体）的速度场；
  • $ p: I \times \Omega_\eta \to \mathbb{R} $ 是压力；
  • $ \mathbf{T}: I \times \Omega_\eta \to \mathbb{R}^{3 \times 3} $ 是聚合物溶质的额外应力张量（Oldroyd-B模型）；
  • $ \eta: I \times \omega \to \mathbb{R} $ 是结构的垂直位移函数，描述二维壳的变形；
  • $ \mathbf{f} $ 和 $ g $ 是外部体积力与表面力；
  • $ \mathbf{D}(\nabla_x \mathbf{u}) = \frac{1}{2} (\nabla_x \mathbf{u} + (\nabla_x \mathbf{u})^\top) $ 是应变率张量；
  • $ \mathbf{S} = 2\mathbf{D}(\nabla_x \mathbf{u}) - p \mathbf{I} + \mathbf{T} $ 是总应力张量；
  • $ \varphi_\eta(y) = \varphi(y) + \mathbf{n}(y) \eta(t,y) $ 是参考边界 $ \partial \Omega \subset \mathbb{R}^3 $ 的参数化变形，其中 $ \varphi $ 是参考构型的光滑嵌入，$ \mathbf{n}(y) $ 是单位法向量；
  • $ \det(\nabla_y \varphi_\eta) = \partial_{y_1} \varphi_\eta \times \partial_{y_2} \varphi_\eta $ 是曲面的面积元。

系统边界条件为：

$$ \mathbf{u} \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} \quad \text{在 } I \times \omega \text{ 上} $$

初始条件为：

$$ \eta(0,\cdot) = \eta_0(\cdot), \quad \partial_t \eta(0,\cdot) = \eta_\star(\cdot) \quad \text{在 } \omega \text{ 上} $$

$$ \mathbf{u}(0,\cdot) = \mathbf{u}_0(\cdot) \quad \text{在 } \Omega_{\eta_0} \text{ 上} $$

$$ \mathbf{T}(0,\cdot) = \mathbf{T}_0(\cdot) \quad \text{在 } \Omega_{\eta_0} \text{ 上} $$

该系统在数学上的挑战性源于以下几个特殊性：

  1. 移动边界：流体域 $ \Omega_\eta(t) $ 随结构位移 $ \eta $ 动态变化，导致偏微分方程定义在非固定的时空域上，必须使用Hanzawa变换等工具将其映射到固定参考域进行分析。

  2. 无中心质量扩散：与以往研究（如 [19] 中的 2D-2D-1D 模型）不同，本文中溶质方程（方程 (1.7)）不含扩散项（即无 $ \Delta \mathbf{T} $ 项），这使得系统在溶质子系统中呈现纯双曲型结构，缺乏耗散机制，极大增加了正则性控制与解的存在性证明的难度。

  3. 高维耦合：从 2D-2D-1D（二维流体、二维溶质、一维结构）推广至 3D-3D-2D，空间维度提升导致插值不等式（如 Ladyzhenskaya 不等式）的估计失效，且 Sobolev 空间嵌入阈值升高，对解的正则性提出更高要求。

  4. 非线性与非对易项：溶质方程中的 $ (\nabla_x \mathbf{u}) \mathbf{T} + \mathbf{T} (\nabla_x \mathbf{u})^\top $ 项由于矩阵乘法的非交换性，导致能量估计中出现难以控制的交叉项，必须通过坐标变换（如将张量方程转化为向量形式）来恢复可处理性。

  5. 非平凡边界条件：流体速度在移动边界上满足 $ \mathbf{u} \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} $，而非零位移条件，这要求对Stokes方程在移动域上的最大正则性理论进行全新构建，而该理论此前仅在零边界条件下成立（如 [1]）。

因此，本研究对象是一个具有高度非线性、高维耦合、无耗散溶质、移动边界与非平凡边界条件的复杂多物理场系统，其建模背景源于聚合物材料在柔性结构（如血管壁、人工瓣膜）中的流动行为，具有重要的生物力学与工程应用价值。

2. 主要定理或主要结果

本文的核心成果是建立了该3D-3D-2D耦合系统在小初始加速度条件下的唯一局部强解的存在性与正则性。其数学表述由 定理 1.4 给出。

定理 1.4（主定理）：设 $ T > 0 $ 给定。令 $ (\eta, \mathbf{u}, p) $ 是溶剂-结构子系统（即方程 (1.1)-(1.2)）在定义 1.1 意义下的唯一局部强解，其数据 $ (f, g, \eta_0, \eta_\star, \mathbf{u}_0, \mathbf{T}) $ 满足正则性条件 (1.3)，且满足以下小性条件：

$$ \int_I \left( \|f\|_{L^2(\Omega_\eta)}^2 + \|\mathbf{T}\|_{W^{1,2}(\Omega_\eta)}^2 + \|g\|_{L^2(\omega)}^2 \right) dt + c \left( \|\eta_\star\|_{\hat{W}^{1,2}(\omega)}^2 + \|\eta_0\|_{W^{3,2}(\omega)}^2 + \|\mathbf{u}_0\|_{W^{1,2}(\Omega_{\eta_0})}^2 \right) \leq \epsilon $$

其中 $ \epsilon > 0 $ 是一个足够小的常数（依赖于参考几何与初始结构），$ c > 0 $ 是常数。则存在一个时间 $ T^* \in (0, T] $，使得整个耦合系统 (1.4)–(1.11) 在定义 1.3 意义下存在唯一的强解 $ (\eta, \mathbf{u}, p, \mathbf{T}) $，且其数据 $ (f, g, \eta_0, \eta_\star, \mathbf{u}_0, \mathbf{T}_0) $ 满足正则性条件 (1.12)。

定义 1.3 对强解的正则性要求极高，具体为：

  • 结构位移 $ \eta \in W^{1,\infty}(I; W^{3,2}(\omega)) \cap L^\infty(I; W^{5,2}(\omega)) \cap W^{3,2}(I; L^2(\omega)) \cap W^{2,2}(I; W^{2,2}(\omega)) \cap L^2(I; W^{6,2}(\omega)) $，且 $ \|\eta\|_{L^\infty(I \times \omega)} < L $，确保变形不自交；
  • 速度场 $ \mathbf{u} \in W^{2,2}(I; L^2_{\operatorname{div}}(\Omega_\eta)) \cap W^{1,2}(I; W^{2,2}(\Omega_\eta)) \cap L^2(I; W^{4,2}(\Omega_\eta)) $，并满足 $ \mathbf{u} \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} $；
  • 压力 $ p \in W^{1,2}(I; W^{1,2}(\Omega_\eta)) \cap L^2(I; W^{3,2}(\Omega_\eta)) $；
  • 应力张量 $ \mathbf{T} \in W^{1,\infty}(I; W^{2,2}(\Omega_\eta)) \cap L^\infty(I; W^{3,2}(\Omega_\eta)) $。

该定理的结论是局部在时间上的，但其正则性远超以往任何工作。特别地，$ \mathbf{T} \in L^\infty(I; W^{3,2}(\Omega_\eta)) $ 表明应力张量在空间上具有三阶弱导数，时间上具有一阶强导数，这种高阶正则性是实现点态满足方程（a.e.）的关键。

关键假设：

  • 初始结构位移 $ \eta_0 \in W^{5,2}(\omega) $，初始速度 $ \mathbf{u}_0 \in W^{3,2}_{\operatorname{div}}(\Omega_{\eta_0}) $，初始应力 $ \mathbf{T}_0 \in W^{3,2}(\Omega_{\eta_0}) $，均要求极高正则性；
  • 初始加速度 $ \eta_\star \in W^{3,2}(\omega) $，且其范数与 $ \mathbf{u}_0 $、$ \mathbf{T}_0 $ 的范数共同被控制在某个小常数 $ \epsilon $ 内，这是保证解存在性的关键小性条件；
  • 外力 $ f, g $ 与初始应力 $ \mathbf{T} $ 的时间积分需足够小，以抵消非线性项带来的增长。

该定理的突破性在于：它首次在无扩散、三维流体与三维溶质、二维结构的框架下，建立了具有如此高正则性的强解。此前工作（如 [19]）虽在 2D-2D-1D 情形下获得强解，但依赖于溶质中的扩散项提供耗散，而本文在完全无扩散情形下仍能获得类似结果，其难度显著提升。

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

本文的证明策略是解耦-迭代-不动点方法，其核心思想是将三态系统分解为两个子系统：溶剂-结构子系统（S-S）与溶质子系统（Solute），分别求解后再通过Banach不动点定理耦合。

a) 核心分析学工具

  1. Banach不动点定理：这是全文的骨架。作者定义一个映射 $ \mathcal{F} $，其输入为一个近似应力张量 $ \tilde{\mathbf{T}} $，输出为由 $ \tilde{\mathbf{T}} $ 驱动的S-S子系统的解 $ (\eta, \mathbf{u}, p) $，再将此解代入溶质方程得到新的 $ \mathbf{T} $。若该映射在某个Banach空间中是压缩的，则存在唯一不动点，即为原系统的强解。

  2. 最大正则性理论（Maximal Regularity）：这是本文的技术基石。作者在第3节独立证明了移动域上Stokes方程在非零边界条件下的最大正则性估计（定理 3.1），其结论为：

$$ \int_I \left( \|\partial_t \mathbf{u}\|_{L^2(\Omega_\eta)}^2 + \|\mathbf{u}\|_{W^{2,2}(\Omega_\eta)}^2 + \|p\|_{W^{1,2}(\Omega_\eta)}^2 \right) dt \lesssim \|\mathbf{u}_0\|_{W^{1,2}(\Omega_{\eta_0})}^2 + \|\eta_\star\|_{W^{1,2}(\omega)}^2 + \int_I \left( \|f\|_{L^2(\Omega_\eta)}^2 + \|\partial_t^2 \eta\|_{L^2(\omega)}^2 + \cdots \right) dt $$

其中右侧包含高阶结构项（如 $ \|\partial_t \eta \nabla_y^2 \eta\| $），这些项源于非零边界条件 $ \mathbf{u} \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} $ 的非齐次性。该结果是对 [1] 中零边界条件工作的关键推广，为S-S子系统的高阶估计提供了理论保障。

  1. Hanzawa变换：用于将移动域 $ \Omega_\eta(t) $ 映射到固定参考域 $ \Omega $，从而将偏微分方程转化为固定域上的带时变系数的方程。变换 $ \Psi_\eta $ 满足良好的Sobolev空间有界性（公式 (2.3)-(2.5)），使得正则性估计可在参考域中进行。

  2. 特征线方法（Method of Characteristics）：用于求解溶质方程（1.7）。由于无扩散项，该方程是双曲型的，作者通过张量向量化变换（将 $ \mathbf{T} $ 视为 $ \mathbb{R}^9 $ 中的向量），将非对易项 $ (\nabla \mathbf{u}) \mathbf{T} + \mathbf{T} (\nabla \mathbf{u})^\top $ 转化为线性算子作用于向量，从而可沿特征线积分求解，且保持正则性。

  3. 椭圆估计与插值不等式：用于控制压力与高阶导数。例如，利用 $ \operatorname{div} \mathbf{u} = 0 $ 与Stokes方程，通过椭圆正则性得到 $ p \in W^{1,2}(W^{1,2}) $；利用Sobolev嵌入 $ W^{2,2}(\Omega_\eta) \hookrightarrow L^\infty(\Omega_\eta) $（在3D中成立）控制非线性项。

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

  • 非对易项 $ (\nabla \mathbf{u}) \mathbf{T} + \mathbf{T} (\nabla \mathbf{u})^\top $：作者未直接估计该张量乘积，而是将其写为 $ \mathbf{Q}(\mathbf{T}, \nabla \mathbf{u}) = \mathbf{T} \mathbf{D}(\nabla \mathbf{u}) + \mathbf{D}(\nabla \mathbf{u}) \mathbf{T} + \mathbf{W}(\nabla \mathbf{u}) \mathbf{T} - \mathbf{T} \mathbf{W}(\nabla \mathbf{u}) $，其中 $ \mathbf{W} $ 是反对称部分。然后通过向量化，将方程转化为 $ \partial_t \vec{\mathbf{T}} + (\mathbf{u} \cdot \nabla) \vec{\mathbf{T}} = \mathbf{A}(\nabla \mathbf{u}) \vec{\mathbf{T}} + \vec{\mathbf{D}}(\nabla \mathbf{u}) $，其中 $ \mathbf{A} $ 是一个线性矩阵算子，从而可用特征线法求解。

  • 无扩散导致的正则性缺失：在无扩散情形下，无法使用Galerkin方法（因缺乏紧性），作者改用逐点构造与高阶能量估计，依赖于初始数据的极高正则性与小性条件来控制非线性增长。

  • 移动边界与高阶导数耦合：在Stokes最大正则性证明中，作者使用投影法：令 $ \Phi = F_{\text{div}}(\partial_t \eta \, \mathbf{n}) $ 为满足 $ \Phi \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} $ 的溶剂化延拓，令 $ \mathbf{w} = \mathbf{u} - \Phi $，则 $ \mathbf{w} $ 满足零边界条件，可应用 [1] 的结果，再通过三角不等式反推 $ \mathbf{u} $ 的估计。

4. 与之前工作的比较

本文在多个维度上实现了对先前工作的实质性突破。

a) 维度推广：从 2D-2D-1D 到 3D-3D-2D

现有文献（如 [19]）仅在二维流体、二维溶质、一维结构的情形下建立了强解。但在该情形下，Ladyszhenskaya不等式允许 $ W^{1,2} \hookrightarrow L^4 $，从而控制非线性项 $ (\nabla \mathbf{u}) \mathbf{T} $。在三维情形下，$ W^{1,2} \not\hookrightarrow L^\infty $，且 $ \nabla \mathbf{u} \in L^2 $、$ \mathbf{T} \in L^2 $ 时，乘积项 $ \nabla \mathbf{u} \cdot \mathbf{T} \in L^1 $，无法闭合能量估计。本文通过引入三阶空间正则性（$ \mathbf{T} \in W^{3,2} $）和高阶时间导数，成功克服了这一障碍。

b) 无扩散情形下的首次强解构造

[19] 及其后续工作 [17,18] 均依赖于溶质方程中的扩散项（$ \Delta \mathbf{T} $）提供耗散，从而获得弱解、强解与奇异极限。本文首次在完全无扩散情形下（即 $ \Delta \mathbf{T} = 0 $）构造了强解，这在数学上是革命性的，因为此时系统为纯双曲型，解可能产生冲击或奇异性，而本文证明了在小初始加速度条件下，解仍保持光滑。

c) 最大正则性理论的突破

此前所有移动域Stokes问题的工作（如 [1], [2]）均假设边界条件为零（$ \mathbf{u} = 0 $）或为刚性位移。本文首次处理了非零、时变、依赖于结构加速度的边界条件 $ \mathbf{u} \circ \varphi_\eta = (\partial_t \eta) \mathbf{n} $，并建立了包含 $ \partial_t^2 \eta $、$ \partial_t \eta \nabla_y^2 \eta $ 等高阶项的完整估计（定理 3.1）。该结果本身即为独立贡献，为后续研究（如可变形血管中的血液流动）提供了基础工具。

d) 正则性强度的提升

在 [4] 中，作者对S-S子系统（无溶质）已建立了强解，但其正则性为 $ \mathbf{u} \in W^{1,2}(W^{2,2}) $，而本文要求 $ \mathbf{u} \in W^{2,2}(L^2) \cap L^2(W^{4,2}) $，即对速度的时空二阶导数均有控制。这种高阶正则性是保证溶质方程（1.7）中 $ \nabla \mathbf{u} \in L^\infty(L^\infty) $ 的关键，从而使得特征线方法可应用。

综上，本文不仅将已有结果从低维推广到高维，更在无耗散这一极端情形下，通过构建全新的分析框架（最大正则性+特征线+不动点），首次实现了强解的存在性与唯一性，标志着聚合物流体-结构耦合系统研究的一个里程碑。

2. Partial Regularity for the Three-dimensional Stochastic Ericksen–Leslie equations

三维随机Ericksen–Leslie方程的部分正则性

作者: Hengrong Du, Chuntian Wang

类别: math.AP, math.PR, 35Q35, 60H15

发布日期: 2024-01-08

arXiv ID: 2401.03662v2

链接: http://arxiv.org/abs/2401.03662v2

详细分析

中文标题: 三维随机Ericksen–Leslie方程的部分正则性

1. 研究对象和背景

本文研究的是三维环面 $\mathbb{T}^3$ 上带加性白噪声的随机Ericksen–Leslie方程（Stochastic Ericksen–Leslie Equations, SEL），该方程组用于描述向列型液晶流体在热涨落和环境扰动下的动力学行为。液晶是一种介于各向同性液体与固体之间的各向异性相态，其分子呈棒状并具有取向有序性（即分子倾向于指向同一方向），但无位置有序性。Ericksen–Leslie方程是描述此类液晶流体宏观动力学的经典数学模型，它耦合了流体的速度场 $u$、压力 $P$ 以及表征分子取向的“指向矢”场 $d$。

在确定性情形下，简化版本的Ericksen–Leslie方程组（见论文第2页）可写为：

$$ \begin{aligned} & \partial_t u - \Delta u + (u \cdot \nabla) u + \nabla P = -\operatorname{div}(\nabla d \otimes \nabla d), \\ & \operatorname{div} u = 0, \\ & \partial_t d - \Delta d + (u \cdot \nabla) d = -f(d), \end{aligned} $$

其中非线性项 $f(d) = \nabla F(d)$，势能函数为Ginzburg–Landau型：

$$ F(d) = \frac{1}{4}(|d|^2 - 1)^2. $$

该系统由两部分构成：左侧为受Ericksen应力张量 $\nabla d \otimes \nabla d$ 驱动的Navier–Stokes方程，右侧为受对流项 $(u \cdot \nabla)d$ 和势能梯度 $f(d)$ 影响的热流调和映射方程。特别地，当 $|d| \equiv 1$ 时，该系统退化为调和映射热流与Navier–Stokes方程的耦合，这在三维空间中是高度非线性且具有潜在奇点的系统。

本文研究的是其随机版本，即在动量方程中引入加性高斯白噪声 $dW$，得到：

$$ \begin{aligned} & du - \Delta u\,dt + (u \cdot \nabla) u\,dt + \nabla P\,dt = -\operatorname{div}(\nabla d \otimes \nabla d)\,dt + dW, \\ & \operatorname{div} u = 0, \\ & dd - \Delta d\,dt + (u \cdot \nabla) d\,dt = -f(d)\,dt, \end{aligned} \tag{4} $$

其中 $W(t)$ 是定义在概率空间 $(\Omega, \mathcal{F}, \mathbb{P})$ 上、取值于 $D(A^\delta)$（$A$ 为Leray投影下的Stokes算子）的 cylindrical Brownian motion，其协方差算子 $L$ 为非负迹类算子。空间区域为三维环面 $\mathbb{T}^3 = \mathbb{R}^3 / \mathbb{Z}^3$，即周期性边界条件，这避免了边界效应，使分析更聚焦于内在的非线性与随机性相互作用。初始条件满足 $u_0 \in H$，$d_0 \in H^1(\mathbb{T}^3; \mathbb{R}^3)$，且 $|d_0(x)| = 1$ 几乎处处成立，这保证了指向矢始终为单位向量，符合物理约束。

该模型的物理动机源于液晶实验中观察到的动态散射（dynamic scattering）与电液动力不稳定性（electrohydrodynamic instabilities），如Williams域不稳定性，这些现象表明流体速度的随机波动是导致分子取向缺陷与湍流的主要驱动力。因此，将噪声引入速度方程而非取向方程，符合物理机制。此外，三维空间中确定性Ericksen–Leslie方程的全局适定性仍是开放问题，而随机情形下更缺乏系统性研究，尤其在正则性方面几乎空白。

本文研究的核心挑战在于：如何在存在白噪声的干扰下，证明解的存在性并刻画其奇点集合的“小性”。由于噪声的无界性（白噪声在时间上不可积，空间上仅为分布），传统的能量估计与紧性论证面临巨大困难。因此，作者引入“鞅合适弱解”（martingale suitable weak solution）这一概念，它既保留了弱解的全局存在性（允许奇点），又通过局部能量不等式排除了非物理解，从而为正则性分析提供了合理框架。

2. 主要定理或主要结果

本文的核心成果是定理2（Theorem 2），其内容如下：

定理 2（主定理）：假设初始数据 $(u_0, d_0) \in H \times H^1(\mathbb{T}^3; \mathbb{R}^3)$ 满足 $|d_0(x)| = 1$ 对所有 $x \in \mathbb{T}^3$ 成立（即条件 (9)），且白噪声 $W(t)$ 取值于 $D(A^{3/4 + \delta})$，其中 $\delta > 0$。则存在一个全局在时间上的鞅合适弱解 $(\mathcal{B}, u_r, P_r, d_r)$，使得在全概率集上，其奇点集合 $\Sigma$ 的一维抛物Hausdorff测度为零，即：

$$ \mathcal{P}^1(\Sigma) = 0, $$

其中 $\mathcal{P}^k$ 表示在 $\mathbb{T}^3 \times \mathbb{R}_+$ 上关于抛物距离 $\delta((x,t), (y,s)) := \max\{|x - y|, |t - s|\}$ 定义的 $k$-维Hausdorff测度。

该定理包含两个核心断言：一是全局存在性，即在任意时间区间 $[0, T]$ 上，存在满足定义3的鞅合适弱解；二是部分正则性，即解的奇点集合在抛物度量下“几乎不占空间”，其维数不超过1。

为实现该结论，作者首先引入“合适弱解”的概念（Definition 2）：设 $(u, P, d)$ 为弱解，令 $v = u - z$，$\pi = P - Q$，其中 $(z, Q)$ 是线性随机Stokes方程（10）的解：

$$ \begin{aligned} & dz - \Delta z\,dt + \nabla Q\,dt = dW, \\ & \operatorname{div} z = 0, \\ & z(0) = 0. \end{aligned} $$

则 $(v, d)$ 满足非线性系统（11），并满足如下两个关键性质：

  1. 弱连续性与正则性：$v \in L^\infty(0,T; H) \cap L^2(0,T; V)$，$d \in L^\infty(0,T; H^1) \cap L^2(0,T; H^2)$，$\pi \in L^{5/3}(0,T; L^{5/3})$；
  2. 局部能量不等式：对任意 $t \leq T$ 及几乎所有 $s < t$，有：
$$ \begin{aligned} & \int_{\mathbb{T}^3} \left( \frac{1}{2}|v|^2 + \frac{1}{2}|\nabla d|^2 + F(d) \right)(t) + \int_s^t \int_{\mathbb{T}^3} \left( |\nabla v|^2 + |\Delta d - f(d)|^2 \right) \\ & \leq \int_{\mathbb{T}^3} \left( \frac{1}{2}|v|^2 + \frac{1}{2}|\nabla d|^2 + F(d) \right)(s) + \int_s^t \int_{\mathbb{T}^3} \left[ z \cdot ((z + v) \cdot \nabla)v + (z \cdot \nabla)d \cdot (\Delta d - f(d)) \right]. \end{aligned} \tag{14} $$
  1. 局部能量不等式的测试形式（式 (15)）：对任意非负光滑紧支集测试函数 $\varphi \in C_0^\infty(\mathbb{T}^3 \times (0,T))$，有：
$$ \begin{aligned} & \int_{\mathbb{T}^3 \times (0,T)} \left( \frac{|v|^2}{2} + \frac{|\nabla d|^2}{2} + F(d) \right) \partial_t \varphi + \int_{\mathbb{T}^3 \times (0,T)} \left( |\nabla v|^2 + |\Delta d|^2 + |f(d)|^2 \right) \varphi \\ & \leq \int_{\mathbb{T}^3 \times (0,T)} \left( \frac{|v|^2}{2} + \frac{|\nabla d|^2}{2} + F(d) \right) \left[ \partial_t \varphi + (z + v) \cdot \nabla \varphi \right] + \int_{\mathbb{T}^3 \times (0,T)} \pi v \cdot \nabla \varphi \\ & + \int_{\mathbb{T}^3 \times (0,T)} \left( \nabla d \otimes \nabla d - \frac{1}{2} |\nabla d|^2 I_3 \right) : \nabla^2 \varphi + \int_{\mathbb{T}^3 \times (0,T)} (z \cdot \nabla)d \cdot (\Delta d - f(d)) \varphi \\ & + \int_{\mathbb{T}^3 \times (0,T)} \left[ (z + v) \cdot \nabla d \right] \cdot (\nabla \varphi \cdot \nabla)d - \int_{\mathbb{T}^3 \times (0,T)} f(d) \cdot (\nabla \varphi \cdot \nabla)d - 2 \int_{\mathbb{T}^3 \times (0,T)} \nabla f(d) : \nabla d \, \varphi. \end{aligned} \tag{15} $$

上述不等式是“合适弱解”的标志性特征，它控制了能量耗散与非线性项的交互，是证明部分正则性的基石。而“鞅合适弱解”（Definition 3）则是指：存在一个概率空间 $(\Omega_r, \mathcal{F}_r, \{\mathcal{F}_t^r\}, \mathbb{P}_r, \widetilde{W})$，使得其路径几乎必然（a.s.）满足上述确定性合适弱解的定义。这一定义将随机性与确定性正则性分析成功分离，是本文方法论的核心创新。

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

a) 核心分析工具

本文的证明主要依赖于以下分析工具：

  1. 局部能量不等式（Local Energy Inequality）：这是Scheffer、Caffarelli–Kohn–Nirenberg等人在Navier–Stokes方程部分正则性理论中开创的工具，本文将其推广至随机耦合系统。它允许在局部区域控制能量耗散，从而避开全局估计的困难。

  2. 抛物Hausdorff测度与尺度分析（Parabolic Hausdorff Measure）：由于方程具有时间与空间的非对称缩放（$u \sim r^{-1}, d \sim r^0, t \sim r^2$），奇点集合的“大小”需用抛物距离 $\delta((x,t),(y,s)) = \max\{|x-y|, \sqrt{|t-s|}\}$ 来定义。本文证明奇点集合的 $\mathcal{P}^1$ 测度为零，意味着奇点在时空上“稀疏”。

  3. 鞅理论与概率紧性：为处理随机噪声，作者构造了“鞅解”（martingale solution），即解的过程在某个概率空间上满足弱收敛与适应性。这避免了对噪声路径的强收敛要求，是处理随机PDE的标准方法，但本文将其与“合适弱解”结合，属创新性整合。

  4. Sobolev空间与分数阶算子：噪声被假设取值于 $D(A^{3/4+\delta})$，其中 $A = -P\Delta$ 为Stokes算子。这利用了 $A^\alpha$ 的谱分解与插值理论，确保 $z(t)$ 的路径几乎必然有界（Lemma 1），这是整个分析的起点。

  5. Gagliardo–Nirenberg不等式与嵌入定理：用于控制非线性项如 $(v \cdot \nabla)d$、$(z \cdot \nabla)d$ 在 $L^p$ 空间中的估计，尤其在处理 $d$ 的 $H^2$ 估计时至关重要。

b) 关键技术与技巧

  1. 解的分解（Split-up Scheme）：作者将总解 $u = z + v$ 分解为线性随机部分 $z$（由噪声驱动）与非线性确定性部分 $v$。这一技巧源自 [26,79] 对随机Navier–Stokes方程的研究，但本文将其首次应用于耦合的Ericksen–Leslie系统。由于 $z$ 满足 $\|z\|_{L^\infty_{\text{loc}}(T^3 \times [0,\infty))} < \infty$ a.s.（Lemma 1），其“噪声扰动”被控制为有界项，从而将随机问题转化为“有界强迫项”的确定性问题。

  2. 局部能量保持逼近（Energy-inequality-preserving approximation）：为构造全局解，作者采用正则化逼近：

$$ \Phi_\sigma[f] = \psi_\sigma * f, \quad \psi_\sigma(x) = \sigma^{-3} \psi(x/\sigma), $$

并构造近似系统（17）：

$$ \begin{aligned} & \partial_t v_\sigma - \Delta v_\sigma + (z + \Phi_\sigma[v_\sigma]) \cdot \nabla (z + v_\sigma) + \nabla \pi_\sigma = -\nabla \Phi_\sigma[d_\sigma] \cdot (\Delta d_\sigma - f(d_\sigma)), \\ & \partial_t d_\sigma - (z + v_\sigma) \cdot \nabla \Phi_\sigma[d_\sigma] = \Delta d_\sigma - f(d_\sigma). \end{aligned} $$

此逼近的关键在于：保留了局部能量不等式。这与传统Galerkin逼近不同，后者在极限过程中可能丢失能量不等式。作者通过 mollifier 的自守性（$\|\Phi_\sigma[f]\|_{L^p} \leq \|f\|_{L^p}$）与对偶估计，确保逼近解满足与极限解一致的能量不等式（14）–（15），从而保证极限解为“合适弱解”。

  1. 爆破论证（Blow-up Argument）：为证明部分正则性，作者采用反证法：若某点 $z_0 = (x_0,t_0)$ 为奇点，则在缩放尺度 $r \to 0$ 下构造序列 $(v_r, d_r)$，通过能量估计与紧性（如Aubin–Lions引理）证明其收敛至一个“非平凡”全局解，该解在原点处具有“非零能量集中”，但由能量不等式（15）可导出矛盾。此方法借鉴了Caffarelli–Kohn–Nirenberg对Navier–Stokes方程的分析，但本文需处理耦合项与随机噪声的“有界背景”项，技术上更复杂。

4. 与之前工作的比较

本文的工作在多个维度上显著推进了现有研究：

  1. 首次在三维随机Ericksen–Leslie系统中建立部分正则性：尽管确定性三维Ericksen–Leslie方程的部分正则性已在 [65] 中被证明（奇点集 $\mathcal{P}^1$ 测度为零），但随机情形下，此前仅在二维有适定性结果（如 [17,90]）。本文首次将该正则性结论推广至三维随机情形，填补了该领域关键空白。

  2. 放宽噪声正则性要求：文献 [26,79] 在随机Navier–Stokes方程中证明了相同结论，但其噪声需满足更强的条件（如 $W(t) \in D(A^{1/2})$）。本文仅要求 $W(t) \in D(A^{3/4 + \delta})$，这是由耦合项 $(z \cdot \nabla)d$ 的非线性性决定的。由于 $d$ 的梯度需控制在 $L^2$，而 $z$ 的 $L^\infty$ 有界性（依赖于 $\delta > 0$）是关键，故 $3/4$ 是临界指数。本文证明该条件足以支撑整个分析，是理论上的优化。

  3. 构造性逼近方法的创新：文献中处理随机PDE部分正则性多依赖于Galerkin逼近或Feynman–Kac公式，但这些方法难以保持局部能量不等式。本文采用的“局部能量保持正则化”（energy-inequality-preserving regularization）是直接继承自 [60,11] 对确定性Navier–Stokes方程的技巧，但首次将其成功应用于随机耦合系统，且证明了其在极限下仍成立，这为其他复杂流体模型（如磁流体、聚合物流体）的随机分析提供了新范式。

  4. 突破“唯一性”依赖的限制：以往对随机PDE的正则性研究常依赖于强解或唯一性假设。本文完全避免此假设，仅通过“合适弱解”的框架，证明了所有路径（几乎必然）都具有相同的正则性结构，这是对“随机正则性”本质的深刻揭示。

综上，本文不仅首次建立了三维随机Ericksen–Leslie方程的部分正则性，更重要的是，它提供了一套可推广至其他随机复杂流体系统的通用方法论，其核心思想——“噪声背景有界化 + 局部能量保持逼近 + 爆破论证”——将成为该领域未来研究的标准工具。

3. Quasi-Gaussianity of the 2D stochastic Navier-Stokes equations

二维随机Navier-Stokes方程的准高斯性

作者: James Coe, Martin Hairer, Leonardo Tolomeo

类别: math.PR, math.AP, 60H15, 35Q30

发布日期: 2025-10-15

arXiv ID: 2510.13460v1

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

详细分析

中文标题: 二维随机Navier-Stokes方程的准高斯性

1. 研究对象和背景

本文研究的是二维环面 $\mathbb{T}^2 = (\mathbb{R}/2\pi\mathbb{Z})^2$ 上带加性噪声的随机Navier-Stokes方程（SNS）的长期行为，特别关注其不变测度的结构特性。该系统由如下随机偏微分方程（SPDE）描述：

$$ \partial_t u + (u \cdot \nabla) u = \Delta u - \nabla p + \sqrt{2} \xi, \quad \nabla \cdot u = 0, $$

其中 $u: \mathbb{R} \times \mathbb{T}^2 \to \mathbb{R}^2$ 是速度场，$p$ 是压强，$\xi$ 是时间白噪声、空间有色噪声的高斯随机场，形式为 $\xi = |\nabla|^{-\alpha} \xi_0$，$\xi_0$ 是空间-时间白噪声，$\alpha > 0$ 控制噪声的正则性。该方程在物理上被广泛用于模拟二维湍流，其数学结构兼具非线性对流项 $(u \cdot \nabla)u$ 与线性耗散项 $\Delta u$，并受高斯噪声驱动，因此是研究随机动力系统不变测度性质的经典模型。

本文进一步推广至“低粘性”（hypoviscous）情形，即用分数阶拉普拉斯算子 $-|\nabla|^{2\gamma}$ 替代标准拉普拉斯算子 $\Delta$，得到更一般的方程：

$$ \partial_t u + (u \cdot \nabla) u = -|\nabla|^{2\gamma} u - \nabla p + \sqrt{2} \xi, \quad \nabla \cdot u = 0, $$

其中 $\gamma \in (2/3, 1]$，噪声正则性参数 $\alpha$ 与 $\gamma$ 满足 $\alpha > 2 - \gamma$，以确保解在空间 $C^{\alpha - \kappa}$（$\kappa > 0$）中具有良好的局部正则性。此处，$|\nabla|^\beta$ 表示由傅里叶乘子 $|n|^\beta$ 定义的分数阶导数算子，作用于傅里叶系数 $\hat{f}(n)$ 上。系统所有变量均假设具有零均值，以保证 $|\nabla|^{-\alpha}$ 等算子定义良好，且热半群呈指数衰减。

研究区域为二维环面 $\mathbb{T}^2$，这在数学上避免了边界效应，使傅里叶分析成为主要工具。函数空间方面，作者使用Hölder-Besov空间 $C^{\beta}(\mathbb{T}^2; \mathbb{R}^2)$，其范数定义为：

$$ \|f\|_{C^\beta} = \sup_{N \text{ dyadic}} N^\beta \|P_N f\|_{L^\infty}, $$

其中 $P_N$ 是Littlewood-Paley投影算子，用于分解函数在不同频率尺度上的成分。该空间能精确刻画解的局部正则性，尤其适用于处理非线性项的正则性损失问题。此外，作者考虑的是发散自由（divergence-free）且均值为零的向量场空间，这保证了Leray投影 $P$ 可以有效消除压强项，将方程转化为无压强形式。

该研究的核心背景是：尽管随机Navier-Stokes方程的不变测度的存在性与唯一性在较早工作中（如 [FM95, DPD02, HM06]）已得到充分研究，但其结构特性——特别是其是否“近似高斯”——长期悬而未决。一个自然的问题是：该系统的不变测度 $\rho$ 是否与对应的线性 Ornstein-Uhlenbeck（OU）过程的不变测度 $\mu_\alpha$ 等价？OU过程由下式给出：

$$ \partial_t u = \Delta u + \sqrt{2} \xi. $$

当 $\alpha = 0$ 时（即空间-时间白噪声），[DPD02] 已证明 $\rho = \mu_0$，两者完全相同。然而，对于 $\alpha > 0$ 的有色噪声情形，传统方法（如 [MS05]）表明，仅当耗散项为超粘性（hyperviscous）形式 $-(-\Delta)^{1+\delta}$ 且 $\alpha$ 足够大时，$\rho \sim \mu_\alpha$ 才成立。近期研究 [HKN24] 更指出，对一般非线性SPDE，若噪声正则性低于某一临界阈值，不变测度通常与高斯测度奇异（singular）。但本文发现，Navier-Stokes非线性项具有特殊的对称性与守恒结构（如能量守恒、涡量输运），使得即使在 $\alpha > 0$ 且为标准粘性（$\gamma=1$）或低粘性（$\gamma < 1$）情形下，该等价性仍能成立——这打破了先前认为的“临界阈值不可逾越”的认知。

2. 主要定理或主要结果

本文的核心成果由两个主要定理构成，分别处理标准粘性与低粘性情形下的不变测度等价性问题。

定理 1.1（标准粘性情形）：
对任意 $\alpha > 0$，方程 (1.1) 存在唯一的不变测度 $\rho$，且该测度与对应的 Ornstein-Uhlenbeck 过程的不变测度 $\mu_\alpha$ 等价，即 $\rho \sim \mu_\alpha$。

定理 1.2（低粘性情形）：
设 $\gamma \in (2/3, 1]$ 且 $\alpha > 2 - \gamma$，则方程 (1.3) 存在全局解，并在空间 $C^{\alpha - \kappa}$（$\kappa \in (0,\alpha)$）中存在唯一的不变测度 $\rho$，且 $\rho \sim \mu_\alpha$。

这两个定理的结论均建立在一系列先验条件之上：

  • 噪声结构：$\xi = |\nabla|^{-\alpha} \xi_0$，$\xi_0$ 为空间-时间白噪声，且为散度自由的向量场。
  • 初始条件：$u_0$ 与 $\xi$ 均为零均值，确保解在零均值空间中演化。
  • 正则性要求：$\alpha > 0$ 保证非线性项在 $C^{\alpha - \kappa}$ 中局部适定；$\gamma > 2/3$ 保证线性化算子满足Bismut-Elworthy-Li公式所需的正则性；$\alpha > 2 - \gamma$ 保证在低粘性情形下能利用涡量守恒（enstrophy）获得全局先验界。
  • 函数空间：所有解与测度均定义在发散自由、零均值的Hölder-Besov空间 $C^{\alpha - \kappa}$ 中。

定理 1.1 的突破在于：它首次证明了在标准粘性（$\gamma=1$）且任意 $\alpha > 0$ 下，SNS的不变测度仍与高斯测度等价。而此前的 [MS05] 仅在超粘性（$\gamma > 1$）与大 $\alpha$ 下成立，[HKN24] 则暗示在 $\alpha \leq 1$ 时应为奇异。定理 1.2 更进一步，将等价性推广至能量超临界情形（$\gamma < 1$），此时 $L^2$ 能量不再守恒，必须依赖涡量 $w = \nabla \wedge u$ 的 $L^2$ 范数（即 enstrophy）来控制解的全局行为。作者证明，即使在 $\gamma \in (2/3, 1)$ 的低粘性区间，只要 $\alpha > 2 - \gamma$，仍能建立全局存在性与唯一不变测度，并维持与高斯测度的等价性。

此外，作者在 Remark 1.5 中提出一个更强的局部结果：当 $\gamma \in (2/3, 1)$ 且 $\alpha > 2 - 2\gamma$ 时，即使全局解尚未被证明，对任意 $t > 0$，解在爆破前的分布 $ \text{Law}(u(t); T_* > t) $ 仍绝对连续于 $\mu_\alpha$。这表明等价性可能在更低正则性下也成立，为后续研究留下空间。

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

a) 核心分析学工具

本文的证明建立在一系列现代SPDE分析工具之上：

  1. Girsanov变换：这是本文最核心的工具。Girsanov定理允许通过改变概率测度，将一个带漂移的SDE/SPDE的分布与一个无漂移过程（如OU过程）的分布相比较。作者将其推广为“时间平移Girsanov方法”（time-shifted Girsanov method），该方法最早由 [MS05] 提出，用于比较非线性方程与线性方程的分布。本文的关键创新在于将其推广至两个非线性方程之间的直接比较。

  2. 时间平移技巧：作者引入一个关键的构造：考虑一个辅助过程 $Z$，使得方程

    $$ \partial_t w = -|\nabla|^{2\gamma} w - (Kw \cdot \nabla) w + Z + \sqrt{2} \xi $$

    在固定时间 $T$ 与另一非线性方程（“扭曲”系统）的解重合。通过将 $Z(t)$ 替换为 $2 P^{(\gamma)}_{T-t} Z(2t - T) \mathbf{1}_{t > T/2}$，利用热半群的半群性质，使该漂移项在时间上“对称化”，从而使其属于Cameron-Martin空间 $L^2([0,T]; H^{\alpha - \gamma})$，满足Girsanov定理的应用条件。

  3. 线性化与积分路径法：为处理非线性项带来的复杂性，作者不直接对原方程应用Girsanov，而是对线性化方程（即解的变分方程）进行分析。令 $J_{s,t}$ 为解的Jacobi算子，其满足：

    $$ \partial_t J_{s,t} = -|\nabla|^{2\gamma} J_{s,t} - (Kw \cdot \nabla) J_{s,t} - (K J_{s,t} \cdot \nabla) w, \quad J_{s,s} = \text{id}. $$

    作者通过对 $J_{s,t}$ 的逐点估计，构造“无穷小Girsanov位移”，并积分得到整体的Radon-Nikodym导数，从而绕过直接处理非线性漂移项的困难。

  4. 先验估计与紧性：为保证不变测度存在，作者建立了一系列先验界，特别是对涡量 $w = \nabla \wedge u$ 的估计：

    $$ \partial_t \|w\|_{L^2}^2 \leq -\||\nabla|^\gamma w\|_{L^2}^2 + C \|w\|_{C^{1-\gamma}}^2 \|\eta\|_{C^{1-\gamma}}^2 + C(1 + \|\eta\|_{C^{1-\gamma}}^4), $$

    其中 $\eta$ 是随机卷积项。结合Fernique定理与Kolmogorov连续性准则，作者证明解在 $C^{\alpha - \kappa}$ 中一致紧，从而通过Krylov-Bogolyubov方法获得不变测度。

  5. 分数阶Leibniz法则：用于估计非线性项的正则性损失，如：

    $$ \||\nabla|^{1-\gamma}(u \eta)\| \lesssim \|u\|_{C^{1-\gamma}} \|\eta\|_{L^\infty} + \|u\|_{L^4} \||\nabla|^{1-\gamma} \eta\|_{L^4}. $$

b) 关键技术与难题突破

本文的核心技术难题在于：非线性项 $(Kw \cdot \nabla)w$ 的正则性不足，无法直接满足Girsanov所需的Cameron-Martin空间条件。传统方法要求漂移项 $Z = (Kw \cdot \nabla)w - \text{线性化项}$ 属于 $L^2 H^{\alpha - \gamma}$，但直接估计表明其正则性仅为 $C^{(2\alpha - 2) \wedge (\alpha - 1)}$，在 $\gamma < 1$ 时远低于要求。

作者的突破性技巧是：

  • 构造“扭曲”系统：引入一个“伪线性”方程：

    $$ \partial_t v = -|\nabla|^{2\gamma} v - |\nabla|^{-\alpha}(Kv \cdot \nabla)|\nabla|^\alpha v + \sqrt{2} \xi, $$

    此系统虽非物理真实，但其非线性项恰好使白噪声测度成为其不变测度（因 $\langle \hat{v}, (K\hat{v} \cdot \nabla) \hat{v} \rangle = 0$），类似于 $\alpha=0$ 时的DPD02情形。

  • 比较两个非线性系统：不比较非线性系统与线性系统，而是比较两个非线性系统——原系统与“扭曲”系统。通过构造一个适配于噪声滤波的漂移 $Z$，使二者在时间 $T$ 的分布一致，再应用Girsanov定理，从而绕过对“线性近似”的依赖。

  • 正则性损失的精细控制：由于线性化算子 $J_{s,t}$ 的正则性受 $w$ 限制（最多提升 $2\gamma - 1$ 阶导数），作者证明：只有当 $(2\gamma - 1) + (\alpha - 1) > \alpha - \gamma$ 时，即 $\gamma > 2/3$，才能保证Girsanov密度的平方可积性。该条件成为定理1.2中 $\gamma > 2/3$ 的关键来源。

4. 与之前工作的比较

本文在多个维度上显著推进了前人工作，其突破性体现在对正则性阈值的彻底放宽与方法论的根本性拓展。

首先，与 [MS05] 的比较：[MS05] 首次证明在超粘性情形（$-\Delta \to -(-\Delta)^{1+\delta}, \delta > 0$）下，SNS的不变测度与高斯测度等价，但要求 $\alpha$ 足够大（$\alpha > 1$ 以上）。本文则在标准粘性（$\gamma=1$）下，对任意 $\alpha > 0$ 均成立。这意味着，即使噪声正则性极低（如 $\alpha = 0.1$），只要非线性项为Navier-Stokes型，等价性仍成立。这直接反驳了 [HKN24] 所暗示的“临界阈值不可逾越”的普遍性，揭示了Navier-Stokes非线性项的特殊结构（如散度自由性、涡量输运守恒）具有超越一般SPDE的“抗奇异”能力。

其次，与 [MRS22, HKN24] 的比较：[HKN24] 通过构造一个在高斯测度下收敛、在非线性测度下发散的“观测量”（observable），证明了在一般非线性SPDE中，等价性仅在满足特定正则性条件时成立。但本文指出，Navier-Stokes非线性项具有天然的“抵消结构”：在 $L^2$ 内积中，$\langle (u \cdot \nabla)u, u \rangle = 0$，且其涡量形式 $\langle (Kw \cdot \nabla)w, w \rangle = 0$，导致 [HKN24] 中的发散量在本系统中恒为零。这表明，非线性项的结构而非仅其正则性，才是决定测度等价性的关键。

第三，对低粘性情形的首次突破：在 $\gamma < 1$ 的能量超临界情形下，$L^2$ 能量无界，传统方法失效。[Fla94] 仅证明 $\alpha > 1/2$ 时全局解存在，但未涉及测度等价性。本文首次在 $\gamma \in (2/3, 1]$ 且 $\alpha > 2 - \gamma$ 下，建立全局存在性、唯一不变测度与高斯等价性三重结果，这是对随机流体力学理论的重大拓展。

最后，方法论革新：作者将“时间平移Girsanov方法”从“非线性 vs 线性”推广至“非线性 vs 非线性”，并结合线性化路径积分，构建了一套全新的“结构敏感型”比较框架。该框架不仅适用于Navier-Stokes，也为研究其他具有守恒律的SPDE（如Euler方程、MHD系统）提供了通用工具。

综上，本文不仅在技术上解决了长期悬而未决的“准高斯性”问题，更在理论上揭示了非线性结构对随机系统不变测度的深刻调控作用，标志着从“正则性导向”向“结构导向”的研究范式转变。

4. Asymptotic stability of the Kolmogorov flow at high Reynolds numbers

高雷诺数下Kolmogorov流的渐近稳定性

作者: Qi Chen, Hao Jia, Dongyi Wei, Zhifei Zhang

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13181v1

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

详细分析

中文标题: 高雷诺数下Kolmogorov流的渐近稳定性

1. 研究对象和背景

本文研究的是二维不可压Navier-Stokes方程（2D incompressible Navier-Stokes equations）在非正方形环面 $ \mathbb{T}_p \times \mathbb{T}_{2\pi} $ 上的Kolmogorov流的渐近稳定性，其中 $ p = 2\pi\kappa $，$ \kappa < 1 $。该方程组在物理上描述了粘性流体在周期性边界条件下的运动，其数学形式为：

$$ \begin{cases} \partial_t \mathbf{U} - \nu \Delta \mathbf{U} + \mathbf{U} \cdot \nabla \mathbf{U} + \nabla P = 0, \\ \nabla \cdot \mathbf{U} = 0, \\ \mathbf{U}|_{t=0} = \mathbf{U}_0(x,y), \end{cases} $$

其中 $ \mathbf{U} = (U_x, U_y) $ 为速度场，$ P $ 为压力，$ \nu = 1/\mathrm{Re} \in (0,1) $ 为粘性系数，$ \mathrm{Re} $ 为雷诺数。当 $ \nu \ll 1 $ 时，系统处于高雷诺数 regime，此时粘性效应微弱，流体行为主要由惯性主导，表现出丰富的非线性动力学特性。

通过引入涡量 $ \Omega = \partial_x U_y - \partial_y U_x $，上述方程可等价转化为涡量方程：

$$ \partial_t \Omega - \nu \Delta \Omega + \mathbf{U} \cdot \nabla \Omega = 0, \quad \Omega|_{t=0} = \Omega_0(x,y). $$

Kolmogorov流是该系统的一个重要稳态解，其形式为 $ \mathbf{U}_{\text{bar}} = (e^{-\nu t} \cos y, 0) $，对应的涡量为 $ \Omega_{\text{bar}} = e^{-\nu t} \sin y $。该流动是二维湍流中典型的“条带态”（bar state），因其在空间上仅沿 $ y $ 方向变化、频率最低，故在粘性耗散下衰减最慢，成为系统在中长期演化中的主要“亚稳态”（metastable state）。在非正方形环面（$ \kappa \ne 1 $）上，Kolmogorov流比偶极子态（dipole state）更易被观测到，是理解层流向湍流转捩（laminar-to-turbulent transition）的关键模型。

本文研究的区域 $ \mathbb{T}_p \times \mathbb{T}_{2\pi} $ 具有特殊几何结构：x方向周期为 $ p = 2\pi\kappa < 2\pi $，y方向周期为 $ 2\pi $，即长宽比 $ \kappa < 1 $。这一非对称性至关重要，因为它避免了偶极子态的主导作用，使Kolmogorov流成为唯一稳定的条带结构。边界条件为周期性，无物理边界，这排除了边界层效应的干扰，使研究聚焦于内部剪切流的稳定性机制。

初始扰动 $ \Omega_0 $ 被限制在 Sobolev 空间 $ H^3(\mathbb{T}_p \times \mathbb{T}_{2\pi}) $ 中，且要求其与Kolmogorov涡量 $ \sin y $ 的距离满足 $ \|\Omega_0 - \sin y\|_{H^3} \ll \nu^{1/3} $。这一函数空间的选择具有深刻意义：$ H^3 $ 是足够低的正则性空间，使得结果具有实际可验证性；同时，它又足以控制非线性项和高阶导数，确保分析的可行性。该空间的低正则性要求是本文的重要突破之一。

本文研究的核心现象是：即使在高雷诺数（低粘性）条件下，若初始扰动足够小（以 $ \nu^{1/3} $ 为尺度），系统将经历一个三阶段的渐近演化过程：（1）快速收敛到一个接近Kolmogorov流的剪切流；（2）缓慢趋近于Kolmogorov流本身；（3）最终因粘性耗散而衰减至零。这一过程跨越多个时间尺度，尤其在 $ 0 < t \leq \nu^{-1} $ 的过渡期内，展现出由“无粘阻尼”（inviscid damping）、“增强耗散”（enhanced dissipation）和“涡量耗尽”（vorticity depletion）共同主导的复杂动力学行为。这些机制在Kolmogorov流中尤为复杂，因其速度剖面 $ \cos y $ 存在临界点（即 $ \cos y = 0 $ 的位置，对应 $ y = \pi/2, 3\pi/2 $），导致涡量在这些点附近被“耗尽”或“抑制”，从而显著削弱了传统增强耗散的效果。

2. 主要定理或主要结果

本文的核心结果由定理1.1给出，其表述如下：

定理 1.1（渐近稳定性）
存在常数 $ \epsilon_1 > 0 $，使得对所有 $ 0 < \nu \leq \epsilon_1 $，若初始涡量 $ \Omega_0 $ 满足

$$ \|\Omega_0 - \sin y\|_{H^3(\mathbb{T}_p \times \mathbb{T}_{2\pi})} \leq \epsilon_1 \nu^{1/3}, $$


则方程 (1.2) 的解 $ \Omega(t) $ 全局存在，且满足如下渐近估计：

对于所有 $ t \geq 0 $，有：

$$ \|\Omega_{\ne}(t)\|_{L^2(\mathbb{T}_p \times \mathbb{T}_{2\pi})} \leq C\left( \nu^{1/3} + \langle \nu t^3 \rangle^{-1} \right) e^{-\epsilon \nu^{1/2} \min(t, \nu^{-1}) - \nu t} \|\Omega_0 - \sin y\|_{H^3}, $$
$$ \|U_{\ne}(t)\|_{L^2(\mathbb{T}_p \times \mathbb{T}_{2\pi})} \leq C\left( \nu^{1/3} + \langle t \rangle^{-1} \right) e^{-\epsilon \nu^{1/2} \min(t, \nu^{-1}) - \nu t} \|\Omega_0 - \sin y\|_{H^3}, $$
$$ \|\Omega(t) - e^{-\nu t} \sin y\|_{L^2(\mathbb{T}_p \times \mathbb{T}_{2\pi})} \leq C e^{-\nu t / 2} \|\Omega_0 - \sin y\|_{H^3}. $$

其中，$ \langle a \rangle := \sqrt{1 + |a|^2} $ 为日本括号，$ h_{\ne} = h - \frac{1}{p} \int_{\mathbb{T}_p} h(x', y) dx' $ 表示对 $ x $ 的零模投影（即去除平均值），$ C, \epsilon > 0 $ 为与 $ \nu $ 无关的常数。

该定理的结论包含三个层次的渐近行为：
（1）快速衰减阶段（$ t \ll \nu^{-1/3} $）：涡量 $ \Omega_{\ne} $ 以 $ \langle \nu t^3 \rangle^{-1} $ 的速率衰减，远快于传统增强耗散 $ e^{-c\nu^{1/2}t} $ 的预期；
（2）过渡阶段（$ t \sim \nu^{-1/3} $ 至 $ t \sim \nu^{-1} $）：非零模涡量和速度场在 $ \nu^{1/3} $ 尺度下被控制，且衰减率由 $ \langle \nu t^3 \rangle^{-1} $ 和 $ \langle t \rangle^{-1} $ 共同决定；
（3）长期衰减阶段（$ t \gg \nu^{-1} $）：x-依赖部分已可忽略，系统由粘性主导，整体涡量以 $ e^{-\nu t} $ 指数衰减。

关键突破在于阈值 $ \nu^{1/3} $ 的获得。此前，对于非单调剪切流（如Kolmogorov流），由于存在临界点，传统增强耗散速率仅为 $ \nu^{1/2} $，因此人们自然预期最优阈值应为 $ \nu^{1/2} $，远小于 $ \nu^{1/3} $。但本文证明，即使在Kolmogorov流存在临界点的情况下，通过引入“涡量耗尽”机制，仍可将稳定性阈值提升至与Couette流相同的 $ \nu^{1/3} $。这表明，尽管Kolmogorov流的线性增强耗散较弱，但其非线性稳定机制（特别是涡量在临界点处的抑制）足以补偿这一劣势，从而实现与单调剪切流相当的稳定性。

此外，定理中 $ H^3 $ 正则性要求极低，远低于以往对非线性稳定性所需的Gevrey正则性（如 $ \mathcal{G}^s $，$ s < 1 $），这使得结果更具物理可实现性。同时，衰减因子 $ \langle \nu t^3 \rangle^{-1} $ 是首次在非单调剪切流中被严格证明，揭示了涡量耗尽与增强耗散之间的深层耦合效应。

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

本文的证明体系构建于一个多层次、多机制协同的分析框架，融合了线性与非线性分析、空间-时间估计、拟线性近似与多时间尺度分析。

a) 核心分析工具：

  1. 拟线性近似方法（Quasilinear Approximation）：作者将解分解为线性部分 $ \omega_L $ 和非线性部分 $ \omega_e $，即 $ \omega = \omega_L + \omega_e $，其中 $ \omega_L $ 是线性化Navier-Stokes方程的解，$ \omega_e $ 满足一个高阶小量（$ O(\nu^{2/3}) $）的非线性方程。该方法首次在非单调流中成功应用，将原本复杂的非线性稳定性问题转化为对高阶小扰动的控制问题，极大降低了非线性项的耦合强度。

  2. 空间-时间估计（Space-time Estimates）：作者建立了一组新的空间-时间范数 $ X_I $，定义为：

    $$ \|f\|_{X_I} = \|f\|_{L^\infty(I;L^2)} + \nu^{1/2} \|\nabla f\|_{L^2(I;L^2)} + \||D_x|^{1/2} \nabla \Delta^{-1} f\|_{L^2(I;L^2)} + \||\sin y|^{1/2} \partial_x \nabla \Delta^{-1} f\|_{L^2(I;L^2)}. $$

    该范数结合了能量范数、粘性耗散项和关键的“临界点加权项”，能精确捕捉涡量在 $ \sin y = 0 $ 附近的衰减行为。其有效性依赖于对线性化Navier-Stokes算子的预解式估计（resolvent estimate），特别是对形如 $ (i\tau - \nu \Delta + \cos y \partial_x (1 + \Delta^{-1})) $ 的算子的谱分析。

  3. 多时间尺度分析（Multiple-timescale analysis）：作者区分了三个时间尺度：$ t \sim \nu^{-1/3} $（涡量耗尽主导）、$ t \sim \nu^{-1/2} $（增强耗散开始生效）、$ t \sim \nu^{-1} $（粘性主导）。通过构造依赖时间的权重函数，将不同机制在不同时间尺度的贡献分离并耦合，最终统一到统一的衰减估计中。

b) 关键技术与难题突破：

  1. 涡量耗尽估计（Vorticity Depletion Estimates）：这是本文最核心的创新。作者首次在非单调剪切流中建立了点态估计：

    $$ |\omega(t, x, y_c)| \lesssim \langle t \rangle^{-1}, \quad \text{当 } \sin y_c = 0, $$

    即在临界点 $ y_c = \pi/2, 3\pi/2 $ 处，涡量随时间衰减。该现象源于涡量在临界流线上的“相混合”（phase mixing）被抑制，导致其在这些点处被“耗尽”。作者通过构造一个几乎守恒量：

    $$ \omega_{1k} := \left( \Delta_k + ikt[\Delta_k, B] - k^2 t^2 (1 - B^2) \right) \omega_k, $$

    其中 $ B = \cos y (1 + \Delta_k^{-1}) $，并证明其满足一个强强制不等式：

    $$ \Lambda_1 (-\Delta_k)^s \Lambda_1 + 2(1 - B^2)(-\Delta_k)^{1+s} + 2(-\Delta_k)^{1+s}(1 - B^2) \geq 4k^2 A (-\Delta_k)^s A + \delta(s)(-\Delta_k)^s, $$

    其中 $ \delta(s) \in [0.52, 4] $，$ A = \sin y (1 + \Delta_k^{-1}) $。该不等式将涡量的高阶导数与临界点附近的权重 $ (1 - B^2) = \sin^2 y $ 耦合，从而在傅里叶空间中实现对涡量的精细控制。

  2. 改进的近似解构造：作者不直接求解线性化NS方程，而是通过修改线性化Euler方程的解构造 $ \omega_L $：

    $$ (\omega_L)_k(t,y) = e^{-\nu k^2 (\sin y)^2 t^3 / 3} (\omega_E)_k(t,y), $$

    其中 $ \omega_E $ 是线性化Euler方程的解。指数项 $ e^{-\nu k^2 (\sin y)^2 t^3 / 3} $ 模拟了临界点附近涡量的“局部耗散”，其根源正是涡量耗尽现象。该构造直接将 $ \langle \nu t^3 \rangle^{-1} $ 衰减嵌入近似解，为后续非线性估计提供基础。

  3. 傅里叶分析与实分析工具的结合：为证明强制不等式 (1.15)，作者将问题转化为对两个正序列的不等式证明（引理5.11），通过精细的傅里叶系数估计和组合分析完成，这是纯分析技巧的典范。

4. 与之前工作的比较

本文在多个维度上实现了对已有工作的根本性突破。

a) 阈值的显著提升：在Kolmogorov流的研究中，此前文献（如[72]）仅能证明稳定性阈值为 $ \nu^{2/3} $，而本文首次将阈值提升至 $ \nu^{1/3} $，与Couette流的最优阈值完全一致。这一结果完全颠覆了传统认知：即非单调剪切流因存在临界点，其稳定性应比单调流（如Couette）更弱。本文证明，尽管增强耗散速率从 $ \nu^{1/3} $ 降为 $ \nu^{1/2} $，但涡量耗尽机制的引入足以补偿这一损失，使整体稳定性阈值不变。

b) 正则性要求的降低：以往对非线性稳定性（如Couette流）的证明通常要求初始扰动在Gevrey空间（如 $ \mathcal{G}^s $，$ s < 1 $）中，以控制非线性项的指数增长。而本文仅要求 $ H^3 $ 正则性，远低于Gevrey空间，这使得结果更贴近实际物理实验和数值模拟的初始数据条件。

c) 方法论的普适性突破：本文首次将“涡量耗尽”估计作为核心工具，用于非线性稳定性分析。此前，涡量耗尽仅被视为一种观测现象（如Bouchet-Morita的数值发现），从未被用于控制非线性误差。本文通过构造几乎守恒量 $ \omega_{1k} $，将其转化为可估计的数学对象，开创了“从现象到工具”的新范式。该方法不仅适用于Kolmogorov流，还可推广至其他含临界点的剪切流（如平面Poiseuille流），甚至可能用于三维流动中类似机制的分析。

d) 与Couette流的对比：虽然Couette流的稳定性阈值 $ \nu^{1/3} $ 已由Masmoudi-Zhao[55]和Bedrossian[10]等建立，但其证明依赖于单调性（无临界点）和傅里叶乘子法。本文则在完全非单调、含临界点的情形下，采用完全不同的空间-时间估计+拟线性逼近方法，实现了相同阈值，说明 $ \nu^{1/3} $ 是更普适的稳定性阈值，而非单调性的特例。

综上，本文不仅是对Kolmogorov流稳定性的首次完整证明，更是一次方法论的革命——它证明了即使在最困难的非单调情形下，通过深入挖掘物理机制并将其数学化，仍可达到与最理想情形（Couette流）相同的稳定性水平。这一成果为理解二维湍流中的亚稳态动力学提供了坚实的理论基石。

5. Enhanced dissipation and Taylor dispersion by a parallel shear flow in an infinite cylinder with unbounded cross section

无限横截面圆柱中平行剪切流增强耗散与泰勒色散研究

作者: Te Li, Le Zhang

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13097v1

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

详细分析

中文标题: 无限横截面圆柱中平行剪切流增强耗散与泰勒色散研究

1. 研究对象和背景

本文研究的是在无限长圆柱域 $\Sigma = \mathbb{R} \times \Omega \subset \mathbb{R}^{d+1}$ 上，一个被动标量 $f(x,y,t)$ 在平行剪切流 $\vec{u}(y) = (v(y), 0, \dots, 0)^\top$ 作用下的输运与扩散行为，其中 $\Omega \subset \mathbb{R}^d$ 是光滑的横截面区域，且允许为无界域（如 $\mathbb{R}$、$(-\infty, L)$ 或 $(L, \infty)$ 的乘积）。该过程由如下非定常对流-扩散方程描述：

$$ \partial_t f(x,y,t) + \vec{u}(y) \cdot \nabla f(x,y,t) = \nu \Delta f(x,y,t), \quad (x,y) \in \Sigma, \; t > 0, $$

其中 $\nu > 0$ 为分子扩散系数，$\Delta = \partial_x^2 + \Delta_y$ 为全空间拉普拉斯算子，边界条件为齐次诺依曼条件：$\partial_n f = 0$ 在 $\partial\Sigma = \mathbb{R} \times \partial\Omega$ 上。该方程描述了在剪切流驱动下，标量浓度如何因扩散与平流竞争而随时间演化，是流体动力学中经典问题——“增强耗散”（enhanced dissipation）与“泰勒色散”（Taylor dispersion）的典型模型。

由于系统在 $x$ 方向具有平移不变性，作者对 $x$ 进行傅里叶变换，定义：

$$ \hat{f}(k,y,t) = \int_{\mathbb{R}} f(x,y,t) e^{-ikx} \, dx, \quad k \in \mathbb{R}, \; y \in \Omega. $$

变换后，每个傅里叶模式满足如下方程：

$$ \partial_t \hat{f}(k,y,t) + ik v(y) \hat{f}(k,y,t) = \nu (-k^2 + \Delta_y) \hat{f}(k,y,t). $$

为分离水平扩散项 $-\nu k^2$ 的影响，作者引入时间相位因子：

$$ \hat{f}(k,y,t) = e^{-\nu k^2 t} g(k,y,t), $$

从而得到“拟椭圆”演化方程：

$$ \partial_t g(k,y,t) + i k v(y) g(k,y,t) = \nu \Delta_y g(k,y,t), \quad y \in \Omega, \; t > 0, $$

并附加齐次诺依曼边界条件 $N \cdot \nabla_y g = 0$ 在 $\partial\Omega$ 上。此处 $k \in \mathbb{R}$ 作为参数，决定了频率尺度，而 $\nu \ll 1$ 是研究的核心小参数。该系统的本质是：剪切流 $v(y)$ 通过其空间变化将能量从低频模式（$|k| \ll \nu$）转移到高频模式（$|k| \gg \nu$），从而加速耗散。

与经典有限横截面情形（如 $\Omega$ 为有界域）不同，本研究中 $\Omega$ 可为无界区域，例如 $\Omega = \mathbb{R}$、$\mathbb{R}^d$ 或其半空间组合。这意味着 $k$ 的取值范围为整个实轴 $\mathbb{R}$，而非周期情形下的离散整数集 $\mathbb{Z}$。在周期情形下（$x \in \mathbb{T}$），所有非零频率 $k \neq 0$ 均可触发增强耗散；而在本研究中，由于存在 $k \to 0$ 的低频模式，系统行为呈现显著差异：当 $|k| \leq \nu$ 时，主导机制为泰勒色散（即标量沿流线扩散，有效扩散系数为 $k^2/\nu$）；而当 $|k| \geq \nu$ 时，剪切流的非均匀性引发增强耗散，其衰减速率远超经典扩散的 $\nu k^2$。

特别地，作者指出，若横截面 $\Omega$ 无界，且剪切流 $v(y)$ 在无穷远处退化（如 $v'(y) \to 0$），则即使 $|k| \gg \nu$，系统也可能不呈现指数衰减。例如，当 $v(y) = 1/y^2$ 且 $\Omega = [1, \infty)$ 时（即泰勒-库埃特流），已知线性系统无法指数衰减。因此，本工作的关键突破在于引入“无穷远处非退化条件”：要求 $\liminf_{|y|\to\infty} |v'(y)| \geq c > 0$，以确保即使在无界域中，剪切流仍能有效“搅动”标量，从而驱动增强耗散。

该研究不仅拓展了经典理论至无界域，还首次系统揭示了在无限横截面情形下，增强耗散与泰勒色散的严格分界机制及其对流场渐近行为的依赖性，为理解复杂流体中混合效率提供了新的数学框架。

2. 主要定理或主要结果

本文的核心结果由两个主要定理构成，分别处理一维与高维情形。

定理 1.1（一维情形）：设 $d=1$，$\Omega$ 为 $\mathbb{R}$、$(-\infty, L_1)$ 或 $(L_2, \infty)$，$v \in C^m(\Omega)$，$m \in \mathbb{N}^*$，满足两个条件：

  1. $m$ 阶非退化条件：对任意 $y \in \Omega$，
    $$ |v'(y)| + |v''(y)| + \cdots + |v^{(m)}(y)| > 0; $$
  2. 无穷远处非退化条件：存在常数 $c_0 > 0$，使得
    $$ \liminf_{|y|\to\infty} |v'(y)| \geq c_0 > 0. $$

则存在正常数 $C, c > 0$，使得对任意 $\nu > 0$、$k \neq 0$ 和初值 $g_0 \in L^2(\Omega)$，方程 (1.5) 的解满足：

$$ \|g(k,t)\|_{L^2(\Omega)} \leq C e^{-c \lambda_{\nu,k} t} \|g_0\|_{L^2(\Omega)}, \quad \forall t \geq 0, $$

其中衰减率参数 $\lambda_{\nu,k}$ 定义为：

$$ \lambda_{\nu,k} = \begin{cases} \nu^{\frac{m}{m+2}} |k|^{\frac{2}{m+2}}, & \text{若 } 0 < \nu \leq |k|, \\ \frac{k^2}{\nu}, & \text{若 } 0 < |k| \leq \nu. \end{cases} $$

该结果表明：当频率 $|k|$ 超过扩散尺度 $\nu$ 时，系统进入“增强耗散”区域，衰减速率由 $ \nu^{\frac{m}{m+2}} |k|^{\frac{2}{m+2}} $ 主导，该速率依赖于剪切流的光滑度 $m$；而当 $|k| \leq \nu$ 时，系统处于“泰勒色散”区域，衰减速率为经典扩散形式 $k^2/\nu$。

定理 1.2（高维情形）：设 $d \geq 2$，横截面 $\Omega = \prod_{j=1}^d \Omega_j$，其中每个 $\Omega_j$ 为上述一维区间之一，且剪切流为可分离形式：

$$ v(y) = v_1(y_1) + v_2(y_2) + \cdots + v_d(y_d). $$

对每个 $j$，假设：

  • （i）$v_j \in C^{m_j}(\Omega_j)$ 且满足 $m_j$ 阶非退化条件；
  • （ii）若 $\Omega_j$ 无界，则 $\liminf_{|y_j|\to\infty} |v_j'(y_j)| \geq c_j > 0$。

则存在常数 $C, c > 0$，使得对任意 $t \geq 0$，

$$ \|g(k,t)\|_{L^2(\Omega)} \leq C \exp\left( -c \sum_{j=1}^d \lambda_{\nu,k,j} t \right) \|g_0\|_{L^2(\Omega)}, $$

其中每一项 $\lambda_{\nu,k,j}$ 与一维情形形式相同：

$$ \lambda_{\nu,k,j} = \begin{cases} \nu^{\frac{m_j}{m_j+2}} |k|^{\frac{2}{m_j+2}}, & \nu \leq |k|, \\ \frac{k^2}{\nu}, & |k| \leq \nu. \end{cases} $$

注意：此处 $\lambda_{\nu,k,j}$ 不依赖于 $j$ 的具体形式，而是统一以 $|k|$ 为参数，表明在可分离结构下，各方向的增强耗散效应可叠加，总衰减率为各方向贡献之和。

两个定理均依赖于一个关键前提：无穷远处非退化条件。若该条件不成立（如 $v'(y) \to 0$），则即使 $|k| \gg \nu$，系统也可能不呈指数衰减（如反例 $v(y)=1/y^2$）。这与经典有界域情形（如 [8]）形成鲜明对比：在有界域中，仅需 $m$ 阶非退化条件（即无临界点）即可保证增强耗散，无需无穷远行为控制。

此外，作者通过伪谱谱上界 $\Psi(H_{\nu,k})$ 的估计，将衰减率转化为算子 $H_{\nu,k} = -\nu \Delta_y + i k v(y)$ 的逆范数下界：

$$ \Psi(H_{\nu,k}) := \sup_{z \in i\mathbb{R}} \| (H_{\nu,k} - z)^{-1} \|_{L^2 \to L^2}, $$

并利用 [17] 中的引理，将指数衰减与 $\Psi(H_{\nu,k})$ 的正下界等价关联，从而将问题转化为对算子谱的精细估计。

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

本文的核心分析工具是伪谱理论与分区域估计法，结合截断函数构造与层集测度控制，以克服无界域带来的非紧性挑战。

a) 核心分析工具

  1. 伪谱上界估计：作者采用 Gallagher–Gallay–Nier [11] 引入的伪谱谱上界 $\Psi(H_{\nu,k})$，并利用 Wei [17] 的结果，将半群衰减 $\|e^{-tH_{\nu,k}}\|$ 与 $\Psi(H_{\nu,k})$ 关联：
$$ \|e^{-tH_{\nu,k}}\| \leq e^{-t \Psi(H_{\nu,k}) + \pi/2}. $$

因此，证明指数衰减等价于证明 $\Psi(H_{\nu,k}) \geq C \lambda_{\nu,k}$。

  1. 分区域估计：为估计 $\|(H_{\nu,k} - ik\lambda)^{-1}\|$，作者将空间 $\Omega$ 分为两部分：

    • 外区域 $\Omega \setminus E_{\lambda,\delta}^m$：其中 $|v(y) - \lambda| \geq \delta^m$，此处虚部 $\Im \langle H_{\nu,k,\lambda} g, g \rangle = k \int (v-\lambda)|g|^2$ 提供强 coercive 项；
    • 内区域 $E_{\lambda,\delta}^m = \{ y \in \Omega : |v(y) - \lambda| < \delta^m \}$：此处 $v(y) \approx \lambda$，剪切效应弱，仅能依赖扩散项 $-\nu \Delta_y$ 控制。

  2. 截断函数构造：作者采用 [8] 中的光滑截断函数：

$$ \chi(y) = \phi\left( \mathrm{sign}(v(y)-\lambda) \cdot \mathrm{dist}(y, E_{\lambda,\delta}^m) / \delta \right), $$

其中 $\phi(t) = \min\{t,1\}$（$t \geq 0$）的奇延拓。该函数满足：

  • $\|\chi\|_{L^\infty} \leq 1$，$\|\nabla_y \chi\|_{L^\infty} \leq 1/\delta$；
  • $\chi(y)(v(y)-\lambda) \geq 0$；
  • 在 $\Omega \setminus E_{\lambda,\delta}^m$ 上，$\chi(y) = \mathrm{sign}(v(y)-\lambda)$。

利用该函数，作者对 $\Im \langle Hg, \chi g \rangle$ 进行估计，导出：

$$ |k| \int_{\Omega \setminus E} |g|^2 \, dy \lesssim \|Hg\| \|g\| + \frac{\nu}{\delta} \|\nabla_y g\| \|g\|. $$

结合能量估计 $\nu \|\nabla_y g\|^2 \leq \|Hg\| \|g\|$，最终得到：

$$ \|g\|^2 \lesssim \max\left\{ |k|^{-1} \delta^{-m}, \nu |k|^{-2} \delta^{-2m-2}, \nu^{-1} m(E_{\lambda,\delta}^m)^2 \right\} \|Hg\| \|g\|. $$

b) 技术性难题的解决技巧

  • 无界域中层集测度控制：这是本工作的核心难点。在有界域中，$m(E_{\lambda,\delta}^m) \leq C \delta$ 自动成立（由 $v$ 的非退化性与紧性）；但在无界域中，若 $v'(y) \to 0$，则 $E_{\lambda,\delta}^m$ 可能无限延伸，导致 $m(E) \to \infty$。作者通过引入无穷远处非退化条件 $\liminf |v'| \geq c_0 > 0$，结合分段单调性引理（Lemma 2.1），证明 $v$ 在无穷远处行为类似线性函数，故 $E_{\lambda,\delta}^m$ 的测度被 $\delta$ 控制，即 $m(E_{\lambda,\delta}^m) \leq C \delta$。

  • 参数选择策略：为使三项平衡，作者选择 $\delta^{m+2} \sim \nu / |k|$，从而统一 $\lambda_{\nu,k} \sim \nu^{m/(m+2)} |k|^{2/(m+2)}$。该选择在 $|k| \geq \nu$ 时有效；而当 $|k| \leq \nu$ 时，直接取 $\delta$ 为常数，利用泰勒色散机制。

  • 高维可分离结构处理：通过将 $v(y)$ 写为 $v_j(y_j)$ 之和，作者将高维问题分解为 $d$ 个独立的一维问题，利用张量积结构，将总衰减率叠加，避免了高维层集的复杂几何分析。

4. 与之前工作的比较

本文最显著的突破在于将增强耗散理论从有界横截面推广至无界横截面，并首次揭示了无穷远处非退化条件的必要性。

在 Coti Zelati 与 Gallay [8] 的开创性工作中，他们研究了横截面 $\Omega$ 为有界域的情形，证明了若 $v(y)$ 为光滑的 Morse 函数（即临界点有限、非退化），则对所有 $k \neq 0$，系统均呈现增强耗散，速率 $\nu^{m/(m+2)} |k|^{2/(m+2)}$。其核心依赖于：(1) 临界点有限；(2) 层集 $E_\lambda$ 为紧集，测度自动可控。

然而，在无界域中，上述假设完全失效：$v(y)$ 可能有无限多个临界点，或趋于常数（如 $v(y) \to 1$ 当 $y \to \infty$），导致层集 $E_{\lambda,\delta}^m$ 无界，测度发散。作者在 [15] 中已证明：当 $v(y) = 1/y^2$ 时，即使 $m$ 任意高，系统仍不指数衰减。这表明：在无界域中，仅靠 $m$ 阶非退化条件不足以保证增强耗散。

本文的突破在于：

  1. 引入必要条件：首次提出并证明“无穷远处非退化”条件 $\liminf_{|y|\to\infty} |v'(y)| \geq c_0 > 0$ 是保证指数衰减的必要且充分条件。该条件确保 $v(y)$ 在无穷远处保持“足够陡峭”，从而阻止层集无限扩张。

  2. 放宽对临界点的要求：在有界域中，[8] 要求 $v$ 为 Morse 函数（临界点有限）；而本文允许 $v$ 在无界域中具有无限临界点，只要其导数在无穷远处不趋于零。例如，$v(y) = y + \sin y$ 在 $\mathbb{R}$ 上有无限临界点，但因 $|v'(y)| \geq 1 - 1 = 0$ 不满足条件；但若 $v(y) = y + \sin(\sqrt{y})$，则 $v'(y) \to 1$，满足条件，系统仍可指数衰减。

  3. 统一处理泰勒色散与增强耗散：本文首次在统一框架下刻画了低频（泰勒色散）与高频（增强耗散）两种机制的切换，且其切换阈值 $\nu \sim |k|$ 与有界域一致，但其在无界域中的适用性由新条件保障。

  4. 维度推广：将结果推广至任意维度 $d \geq 1$，并允许横截面为半无限区间乘积（如 $\mathbb{R} \times [1,\infty)$），极大扩展了模型的物理适用性（如管道流、大气边界层等）。

综上，本文不仅解决了 [8] 的遗留问题，更建立了无界域中增强耗散的完整理论框架，其结论具有深刻的物理意义：即使在无限延伸的流体中，只要剪切流在无穷远处保持“非退化”强度，混合仍可高效进行。这一发现为工程与地球物理中的扩散混合设计提供了坚实的数学基础。

6. The local regularity theory for the Stokes and Navier–Stokes equations near the curved boundary

曲边界附近Stokes与Navier-Stokes方程的局部正则性理论

作者: Hui Chen, Su Liang, Tai-Peng Tsai

类别: math.AP

发布日期: 2025-10-14

arXiv ID: 2510.13028v1

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

详细分析

中文标题: 曲边界附近Stokes与Navier-Stokes方程的局部正则性理论

1. 研究对象和背景

本文研究的是在曲边界附近，满足无滑移（no-slip）或Navier边界条件的非稳态Stokes方程及其对应的Navier-Stokes方程的局部正则性问题。研究的核心对象是不可压缩粘性流体的运动，其动力学由以下Stokes方程组描述：

$$ \begin{cases} \partial_t u - \Delta u + \nabla p = f + \operatorname{div} F, \\ \operatorname{div} u = 0, \end{cases} $$

其中 $ u = (u_1, \dots, u_d) $ 是速度场，$ p $ 是压力，$ f $ 是体积力，$ F = (F_{ij})_{1 \leq i,j \leq d} $ 是表面力张量。当将 $ F $ 替换为 $ -u \otimes u $ 并令 $ f = 0 $ 时，该系统退化为经典的Navier-Stokes方程：

$$ \begin{cases} \partial_t u + u \cdot \nabla u - \Delta u + \nabla p = 0, \\ \operatorname{div} u = 0. \end{cases} $$

这些方程描述了不可压缩粘性流体的运动，是流体力学的基础模型。本文关注的区域 $ \Omega \subset \mathbb{R}^d $（$ d \geq 2 $）具有 $ C^{1,1} $ 或更强的光滑边界 $ \partial \Omega $，且原点位于边界上。为研究局部行为，作者定义了半球形区域 $ B_R = B(0,R) $，$ B_R^+ = B_R \cap \Omega $，以及时间-空间柱体 $ Q_R^+ = B_R^+ \times (-R^2, 0) $，用于分析边界附近解的正则性。

关键的边界条件有两种：其一是无滑移边界条件（no-slip BC）：

$$ u = 0 \quad \text{在} \quad \Sigma = \Gamma \times (-1,0), \quad \text{其中} \quad \Gamma = B_1 \cap \partial \Omega, $$

这表示流体在固体边界上速度完全为零，是工程与物理中最常见的边界条件。其二是Navier边界条件（Navier BC）：

$$ [(2Du + F)n]_{\text{tan}} + \alpha u_{\text{tan}} = 0, \quad u \cdot n = 0 \quad \text{在} \quad \Sigma, $$

其中 $ Du = \frac{1}{2} (\nabla u + (\nabla u)^T) $ 是变形张量，$ n $ 是单位外法向量，$ v_{\text{tan}} = v - (v \cdot n)n $ 表示向量在边界上的切向分量，$ \alpha \geq 0 $ 是摩擦系数。特别地，当 $ \alpha = 2\kappa $（$ \kappa $ 为边界主曲率张量）时，Navier条件简化为涡量边界条件 $ \omega n = 0 $，其中 $ \omega $ 是梯度 $ \nabla u $ 的反对称部分（即涡量矩阵），在三维情形下等价于 $ (\nabla \times u) \times n = 0 $。该条件允许流体在边界上具有滑移，更适用于微流体、多孔介质或非牛顿流体等场景。

本文研究的区域具有曲边界（curved boundary），这是与以往大量研究（如Seregin-Shilkin、Chang-Kang等）的显著区别。以往工作大多集中于平面边界（flat boundary），即 $ \Gamma $ 为平面区域的情形。在平面边界下，通过反射、延拓等对称性技巧，可将边界问题转化为全空间问题，从而利用傅里叶分析或热核估计获得正则性。然而，当边界为曲面时，几何曲率引入非平凡的高阶扰动项，使得标准的延拓方法失效。为处理此问题，作者采用边界拉直（boundary straightening）技术：通过坐标变换 $ T(x) = (x', x_d - \gamma(x')) $ 将曲边界 $ x_d = \gamma(x') $ 映射为平面 $ x_d = 0 $，其中 $ \gamma \in C^{1,1} $ 满足 $ \gamma(0) = 0 $, $ \nabla \gamma(0) = 0 $。此变换后，原方程在拉直后的空间中产生额外的高阶偏微分算子项，例如：

$$ T^{-1} \Delta = \Delta - 2 \nabla \gamma \cdot \nabla \partial_d + |\nabla \gamma|^2 \partial_d^2 - \Delta \gamma \cdot \partial_d, $$

这些项成为分析中的主要障碍。此外，本文在极低正则性假设下开展研究：仅假设速度 $ u $ 和压力 $ p $ 属于 $ L^{q,r}(Q_1^+) $，而无需预先假设其梯度或二阶导数存在。这与以往研究（如Seregin-Shilkin [29]）要求 $ \nabla u \in L^{q,r} $ 作为前提形成鲜明对比。函数空间方面，作者使用混合范数空间 $ L^{q,r}(Q_1^+) = L^r((-1,0); L^q(B_1^+)) $，并定义了临界指数 $ q^* $：

$$ q^* = \begin{cases} \max\{1, \frac{dq}{q+d}\}, & q \neq d, \\ 1+, & q = d, \end{cases} $$

以适应不同维度下的嵌入关系。特别地，当 $ q = d $ 时，$ q^* $ 略大于1，以保证非线性项 $ u \cdot \nabla u $ 的可积性。这些设定使本文的框架能够处理比以往更弱的初始数据，并为Navier-Stokes方程的边界正则性理论提供新基础。

2. 主要定理或主要结果

本文的核心贡献体现在三个关键定理中，分别对应无滑移与Navier边界条件下Stokes方程的局部正则性估计。

定理1.1（无滑移边界条件下的导数估计）：假设 $ 1 < q, r < \infty $，边界部分 $ \Gamma \in C^{1,1} $，且 $ u, p, F \in L^{q,r}(Q_1^+) $，$ f \in L^{q^*,r}(Q_1^+) $，若 $ (u,p) $ 是Stokes方程在 $ Q_1^+ $ 上满足无滑移边界条件 $ u=0 $ 的弱解对（见定义2.3），则其为强解对（见定义2.4），且满足：

$$ \|\nabla u\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q^*,r}(Q_1^+)} + \|F\|_{L^{q,r}(Q_1^+)}, $$

若进一步有 $ f + \operatorname{div} F \in L^{q,r}(Q_1^+) $，则时间导数与二阶空间导数也存在：

$$ \|\partial_t u, \nabla^2 u, \nabla p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f + \operatorname{div} F\|_{L^{q,r}(Q_1^+)}. $$

该定理的突破性在于：它首次在曲边界上证明，即使仅假设 $ u,p \in L^{q,r} $（无任何梯度先验信息），也能推出 $ \nabla u \in L^{q,r} $，甚至 $ \nabla^2 u, \partial_t u, \nabla p \in L^{q,r} $。此前，Seregin-Shilkin [29] 在平面边界下已证明类似结论，但需假设 $ \nabla u \in L^{q,r} $ 作为前提；而本文无需此假设，且适用于任意曲率边界。

定理1.2（Navier边界条件下的梯度估计）：在相同 $ q,r $ 条件下，设 $ \Gamma \in C^{1,1} $，$ u,F \in L^{q,r}(Q_1^+) $，$ f \in L^{q^*,r}(Q_1^+) $，且 $ \nabla u \in L^1(Q_1^+) $，若 $ u $ 是满足Navier边界条件（1.2）的弱解（定义2.5），则：

$$ \|\nabla u\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q^*,r}(Q_1^+)} + \|F\|_{L^{q,r}(Q_1^+)}. $$

该结果将作者先前在平面边界下的结论 [5, Theorem 1.1] 推广至任意维度 $ d \geq 2 $ 的曲边界。值得注意的是，尽管在平面情形下该估计对“非常弱解”成立，但在曲边界情形下，由于技术困难，作者被迫假设 $ \nabla u \in L^1 $ 作为起始正则性，但右侧不包含 $ \|\nabla u\|_{L^1} $，这是极强的“自启动”估计。

定理1.3（Navier边界条件下的高阶导数估计）：若 $ \Gamma \in C^{2,1} $，$ F=0 $，$ u,p,f \in L^{q,r}(Q_1^+) $，$ \nabla u \in L^1(Q_1^+) $，且 $ (u,p) $ 为Navier边界条件下的弱解对（定义2.6），则：

$$ \|\partial_t u, \nabla^2 u, \nabla p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^{q,r}(Q_1^+)}. $$

若进一步 $ \Gamma \in C^{3,1} $ 且 $ \nabla f \in L^{q,r}(Q_1^+) $，则可得到三阶导数估计：

$$ \|\nabla \partial_t u, \nabla^3 u, \nabla^2 p\|_{L^{q,r}(Q_{1/2}^+)} \lesssim \|u\|_{L^{q,r}(Q_1^+)} + \|p\|_{L^{q,r}(Q_1^+)} + \|f\|_{L^r(-1,0; W^{1,q}(B_1^+))}. $$

该结果首次在曲边界上实现了对Navier边界条件下三阶空间导数和二阶压力导数的全局估计，且无需预先假设高阶导数的存在。此前Dong-Kwon [9] 在类似设定下仅在假设 $ \nabla^2 u $ 有界的前提下获得二阶估计，而本文完全摆脱了此类先验假设。

这三个定理共同构成了曲边界附近Stokes方程局部正则性理论的完整框架，其核心突破在于：在极低初始正则性（仅 $ L^{q,r} $）下，通过边界拉直后复杂的扰动项控制，实现了从“弱解”到“强解”的跃迁，并首次在非平坦几何中建立了完整的高阶导数估计。

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

本文的证明体系高度依赖于局部化分析与正则化逼近的结合，其技术核心可归纳为三大创新工具。

a) 核心分析工具：

  1. 边界拉直与坐标变换（Boundary Straightening）：如公式 (2.1)-(2.4) 所示，通过变换 $ T(x) = (x', x_d - \gamma(x')) $ 将曲边界 $ \partial \Omega $ 映射为平面 $ x_d = 0 $。该变换将原方程中的拉普拉斯算子 $ \Delta $ 转化为带曲率扰动的算子：

    $$ T^{-1} \Delta = \Delta - 2 \nabla \gamma \cdot \nabla \partial_d + |\nabla \gamma|^2 \partial_d^2 - \Delta \gamma \cdot \partial_d, $$

    这使得原问题转化为在全空间 $ \mathbb{R}^d_+ \times (-1,0) $ 上的非齐次热方程，但右端出现高阶非线性扰动项。

  2. 偶-偶-奇延拓与磨光（Even-Even-Odd Extension and Mollification）：为处理边界条件，作者定义了速度场的延拓形式：

    $$ \tilde{u}_k = \begin{cases} u_k, & k = 1,\dots,d-1, \\ -\sum_{i=1}^{d-1} \gamma'_{i} u_i + u_d, & k = d, \end{cases} $$

    使得 $ \operatorname{div} T^{-1} \tilde{u} = 0 $ 且在 $ x_d = 0 $ 上满足无滑移条件。随后，作者使用水平-垂直-时间三重磨光（见公式2.10–2.12）：

    $$ E_{\varepsilon_1,\varepsilon_2,\varepsilon_3}^{\pm}(g) = E_{\varepsilon_1}^h \left( E_{\varepsilon_2}^{\pm} \left( E_{\varepsilon_3}(g) \right) \right), $$

    其中 $ E^h $ 为水平方向磨光，$ E^{\pm} $ 为垂直方向偶-奇延拓磨光，$ E $ 为时间磨光。这种构造保证了磨光后函数在边界上仍满足原始边界条件，同时光滑化高阶扰动项。

  3. 能量估计与嵌入理论：作者在磨光后的方程上进行 $ L^{q,r} $ 能量估计，利用Sobolev嵌入（如 $ W^{1,q} \hookrightarrow L^{q^*} $）和插值不等式，建立逐阶导数控制。

b) 关键技术突破：

  1. “正则形式”（Normal Form）的引入：作者定义了磨光后速度场 $ \tilde{v} = E^{\pm}(T^{-1}u \cdot \zeta) $ 的“正则形式”，其比直接拉直后的函数更光滑，作为中间桥梁，使梯度估计可在更弱的假设下完成。

  2. 水平导数控制垂直导数：这是本文最精妙的技巧之一。在平面边界情形下，由 $ \Delta u = 0 $ 可直接得到 $ \partial_d^2 u = -\sum_{i=1}^{d-1} \partial_i^2 u $。在曲边界下，作者通过方程结构（如 $ \partial_t \tilde{v} - \Delta \tilde{v} = \text{perturb} $）证明：水平导数 $ \nabla' \tilde{v} $ 的 $ L^{q,r} $ 控制可推出垂直导数 $ \partial_d \tilde{v} $ 的估计，例如：

    $$ \|\partial_d^2 \tilde{v}\|_{L^{q,r}} \lesssim \|\nabla' \nabla \tilde{v}\|_{L^{q,r}} + \|\text{perturb}\|_{L^{q,r}}, $$

    从而避免了直接估计 $ \partial_d^2 $ 的困难。

  3. 时间导数向空间导数转移：为处理 $ \partial_t u $ 项，作者利用方程结构将时间导数转化为拉普拉斯项与非线性项：

    $$ \partial_t u = \Delta u - \nabla p + f + \operatorname{div} F, $$

    从而将 $ \|\partial_t u\|_{L^{q,r}} $ 的估计转化为对 $ \|\nabla^2 u\|_{L^{q,r}} $ 的控制，再通过迭代闭合估计。这一技巧在公式 (3.32)、(5.34) 中被反复使用，是克服时间导数“非局部性”的关键。

  4. 扰动项吸收：磨光后方程的右端包含形如 $ \gamma' \nabla \partial_d u $、$ |\nabla \gamma|^2 \partial_d^2 u $ 的项，作者通过小参数 $ \|\nabla \gamma\|_{L^\infty} \ll 1 $（由 $ \rho \ll 1 $ 保证）将其吸收至左侧，利用Young不等式与磨光算子的有界性（引理2.2）完成闭合。

4. 与之前工作的比较

本文在多个维度上显著超越了以往研究，其突破性体现在假设的弱化、几何的推广与结果的强化三个方面。

首先，在边界几何上，以往工作如Seregin-Shilkin [29]、Chang-Kang [2]、Dong-Kwon [9] 均集中于平面边界（flat boundary），其核心依赖于反射对称性与傅里叶分析。而本文首次在曲边界（$ C^{1,1} $ 甚至 $ C^{3,1} $）下建立完整正则性理论，克服了曲率带来的非对称性与高阶扰动项，填补了该领域的空白。

其次，在正则性假设上，本文大幅放宽了前提。例如，Seregin-Shilkin [29] 在平面边界下证明 $ \nabla^2 u \in L^{q,r} $ 需假设 $ \nabla u \in L^{q,r} $，而本文定理1.1仅假设 $ u,p \in L^{q,r} $，即从“非常弱解”（very weak solution）出发，直接推出“强解”（strong solution）。这一提升在数学上是本质性的：它意味着解的正则性并非依赖于预先假设，而是方程结构本身的内在性质。

第三，在结果强度上，本文首次在曲边界上获得Navier边界条件下的三阶导数估计（定理1.3），而此前文献（如[5]）仅在平面情形下证明 $ \nabla^3 u \in L^{q,r} $。此外，定理1.2中梯度估计的右侧不含 $ \|\nabla u\|_{L^1} $，这在平面情形下虽已知，但在曲边界下从未实现，体现了“自正则化”（self-improving）的强性质。

特别地，作者指出，压力项 $ \|p\|_{L^{q,r}} $ 的出现是必要的（见评论(iii)），这排除了Kang [16]、Seregin-Sverák [26] 等人构造的反例——这些反例表明，若无压力控制，即使 $ u \in L^\infty $，其梯度也可能无界。本文通过引入压力项，首次在曲边界上建立了“压力-速度”耦合的稳定估计。

最后，本文的应用层面（定理1.5–1.6）具有深远意义：它首次为Navier-Stokes方程的边界正则点（boundary regular point）给出了一个可验证的定义（定义1.4）：若 $ u \in L^\infty $ 且 $ p \in L^{q_0,r_0} $，则该点为正则点，可推出更高阶正则性。这为Caffarelli-Kohn-Nirenberg [1] 的内点部分正则性理论提供了边界版本，是偏微分方程领域的重要进展。

综上，本文不仅将经典理论从平面推广至曲面，更在假设的“最弱化”与结论的“最强化”之间建立了前所未有的桥梁，为不可压缩流体边界层理论奠定了新的分析基础。

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

1. How smooth are restrictions of Besov functions?

Besov 函数的限制有多光滑？

作者: Julien Brasseur

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

发布日期: 2025-09-09

arXiv ID: 2509.07420v3

链接: http://arxiv.org/abs/2509.07420v3

摘要翻译

中文标题: Besov 函数的限制有多光滑？

摘要翻译: 在先前的研究中，我们证明了 Besov 空间仅当 $q \leq p$ 时才具有限制性质。具体而言，我们表明：若 $p < q$，则总能构造出一个函数 $f \in B_{p,q}^s(\mathbb{R}^N)$，使得对几乎所有的 $y \in \mathbb{R}^{N-d}$，其截面函数 $f(\cdot, y)$ 均不属于 $B_{p,q}^s(\mathbb{R}^d)$；而当 $q \leq p$ 时，这种“病态”现象则不会出现。我们进一步证明，当函数 $\Psi$ 满足一个仅依赖于 $p$ 和 $q$ 的简单可和性条件时，这些截面函数实际上属于具有广义光滑性的 Besov 空间 $B_{p,q}^{(s,\Psi)}(\mathbb{R}^d)$。本简短笔记通过证明该刻画的最优性，完成了这一理论图景。

2. Quantitative estimates for flows of regular Lagrangian flows for Hörmander singular kernels and LD vector fields

关于Hörmander奇性核与LD向量场的正则拉格朗日流的定量估计

作者: Henrique Borrin

类别: math.AP, math.DS

发布日期: 2025-10-15

arXiv ID: 2510.13690v1

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

摘要翻译

中文标题: 关于Hörmander奇性核与LD向量场的正则拉格朗日流的定量估计

摘要翻译: 本文针对一类向量场，其导数可表示为满足“Hörmander”条件的基本奇性核与时空可积函数的卷积，获得了与之相关的正则拉格朗日流的定量估计。

3. The Cauchy problem for the integrable RZQ equation

可积RZQ方程的Cauchy问题

作者: John Holmes, Katie Massey, Ryan C. Thompson

类别: math.AP, 35Q35

发布日期: 2025-07-17

arXiv ID: 2507.13535v2

链接: http://arxiv.org/abs/2507.13535v2

摘要翻译

中文标题: 可积RZQ方程的Cauchy问题

摘要翻译: 本文研究了由Reyes、Zhu和Qiao导出的一种新的可积五阶Camassa-Holm（CH）型方程，我们称之为RZQ方程。该方程的m形式与其对应的CH方程的m形式具有显著的相似性。然而，与CH方程不同的是，RZQ方程的非局部形式无法被解释为Burgers方程的非局部扰动。我们证明了RZQ方程对应的初值问题在Sobolev空间$H^s$（$s>7/2$）中满足Hadamard意义下的适定性。进一步地，我们表明，尽管数据到解的映射在较弱拓扑下是Hölder连续的，但在$H^s$拓扑下并不一致连续。此外，RZQ方程对应的初值问题在$H^s$（$s<7/2$）中是不适定的。

4. Rate of convergence for numerical $α$-dissipative solutions of the Hunter-Saxton equation

Hunter-Saxton方程数值α-耗散解的收敛速率

作者: Thomas Christiansen, Katrin Grunert

类别: math.NA, cs.NA, math.AP, 65M15, 65M25 (Primary) 35Q35 (Secondary)

发布日期: 2024-11-12

arXiv ID: 2411.07712v4

链接: http://arxiv.org/abs/2411.07712v4

摘要翻译

中文标题: Hunter-Saxton方程数值α-耗散解的收敛速率

摘要翻译: 我们证明，对于Cauchy问题的α-耗散解，其中α ∈ W^{1, ∞}(ℝ, [0, 1))，在满足如下条件时，其在L^∞(ℝ)范数下的数值逼近收敛阶可达O(Δx^{1/8} + Δx^{β/4})：存在常数C > 0及β ∈ (0, 1]，使得初始空间导数ū_x满足对所有h ∈ (0, 2]有‖ū_x(· + h) − ū_x(·)‖_2 ≤ Ch^β。所导出的收敛速率通过若干数值实验予以验证。

5. Analytic Fourier Integral Operators and a Problem from Seismic Inversion

解析傅里叶积分算子与地震反演中的一个难题

作者: Leonard Busch

类别: math.AP

发布日期: 2025-05-13

arXiv ID: 2505.08472v2

链接: http://arxiv.org/abs/2505.08472v2

摘要翻译

中文标题: 解析傅里叶积分算子与地震反演中的一个难题

摘要翻译: 我们建立了一个一般性结果，表明在满足包括Bolker条件在内的若干条件下，一个经典椭圆型解析傅里叶积分算子（FIO）作用后所得变换的解析波前集，可用于重构原分布的解析波前集。此外，我们为一般的解析二阶双曲型微分算子构造了一个显式的解析参数基本解，其形式为一个经典椭圆型解析FIO。最后，我们将上述结果结合微局部解析延拓与层剥离论证，应用于地震反演中的一个具体问题，证明了在解析框架下线性化算子的单射性。我们之所以能够精确刻画所考虑的FIO如何作用于解析波前集，关键在于采用了波包技术；而显式的解析参数基本解则源于对Hadamard与Hörmander参数基本解构造方法的重新表述。

6. The Neumann problem for the fractional Laplacian: optimal regularity via the Mellin transform

分数阶拉普拉斯算子的诺伊曼问题：通过梅林变换获得最优正则性

作者: Serena Dipierro, Xavier Ros-Oton, Enrico Valdinoci, Marvin Weidner

类别: math.AP, 47G20, 35B65

发布日期: 2025-10-15

arXiv ID: 2510.13340v1

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

摘要翻译

中文标题: 分数阶拉普拉斯算子的诺伊曼问题：通过梅林变换获得最优正则性

摘要翻译: 我们建立了分数阶拉普拉斯算子诺伊曼问题 $(-\Delta)^s u = h$ 在区域 $\Omega$ 内、外部条件 $\mathcal{N}^s u = 0$ 在 $\Omega^c$ 上的解的最优正则性。为此，关键在于建立一个关于具有增长性函数的一维刘维尔定理，该定理我们通过复分析与梅林变换加以证明。
更精确地，我们证明了一个“元定理”，它将 $(0,\infty)$ 上形如 $Lu = 0$ 的一般线性齐次方程的一维解的分类，与一个依赖于算子 $L$ 的显式亚纯函数 $f(z)$ 的（复）零点联系起来。在分数阶拉普拉斯算子配以诺伊曼条件的情形下，我们证明：当 $s \leq 1/2$ 时，所有解属于 $C^{2s+\alpha}$ 类；当 $s \geq 1/2$ 时，所有解属于 $C^{s+\frac{1}{2}+\alpha}$ 类。
更令人惊讶的是，我们证明即使在一维情形下，也存在高度振荡的解，其形式为 $u(x) = x^a \cos(b \log x)$（$x > 0$），其中 $a > 0$、$b > 0$ 依赖于参数 $s$，且当 $s \sim 1$ 时满足 $a < 2s$。

7. An Allen-Cahn tumor growth model with temperature

一种考虑温度效应的Allen-Cahn肿瘤生长模型

作者: Stefania Gatti, Erica Ipocoana, Alain Miranville

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13283v1

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

摘要翻译

中文标题: 一种考虑温度效应的Allen-Cahn肿瘤生长模型

摘要翻译: 本文提出了一种新的非等温Allen-Cahn（Ginzburg-Landau）模型，用于描述肿瘤生长过程。该模型通过微观力方法推导得出，并对其适定性进行了系统研究。特别地，我们证明了该偏微分方程组在局部及全局时间范围内解的存在性与唯一性。

8. Universal Potential Estimates for Mixed Local and Nonlocal Nonlinear Measure Data Problems

混合局部与非局部非线性测度数据问题的通用势估计

作者: Lingwei Ma, Qi Xiong, Zhenqiu Zhang

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13269v1

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

摘要翻译

中文标题: 混合局部与非局部非线性测度数据问题的通用势估计

摘要翻译: 本文建立了具有系数与测度数据的混合局部与非局部 $p$-拉普拉斯型方程的非线性势理论，涵盖超二次与次二次两种情形。我们通过黎兹势与沃尔夫势，证明了关于解及其梯度的一类通用逐点估计。这些估计的实现依赖于对局部项系数施加多种低正则性条件，而对非局部项的核系数仅假设其为可测函数。证明的核心在于引入一种新型的分数阶极大函数，该函数能够同时捕捉局部与非局部特征，并在此基础上建立解及其梯度的此类极大算子的逐点估计。值得注意的是，我们的通用势估计不仅精确刻画了解的振荡行为，还明确了界定其大小的临界情形，从而显著改进了先前工作中已有的逐点势估计。

9. The Boltzmann equation on smooth and cylindrical domains with Maxwell boundary conditions

在光滑与圆柱形域上满足麦克斯韦边界条件的玻尔兹曼方程

作者: Richard Medina Rodriguez

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13260v1

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

摘要翻译

中文标题: 在光滑与圆柱形域上满足麦克斯韦边界条件的玻尔兹曼方程

摘要翻译: 本文研究了玻尔兹曼方程在其流体动力学极限附近于有界域上的适定性问题。我们考虑两类域：一类是具有麦克斯韦边界条件的 $C^2$ 域，其中 Accommodation 系数为连续的空间依赖函数 $\iota \in [\iota_0,1]$，其中 $\iota_0 \in (0,1]$；另一类是圆柱形域，其底面为扩散反射，其余边界为镜面反射。此外，我们采用多项式权重、拉伸指数权重以及逆高斯权重，构建了接近平衡态的柯西理论。需指出的是，所有方法均为定量性的，因此所有常数均为可构造且可追踪的。

10. Large deviation principles for the Gross Pitaevskii Gibbs measure at low temperature

低温下Gross-Pitaevskii Gibbs测度的大偏差原理

作者: Liam Packer, Kihoon Seong, Philippe Sosoe

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

发布日期: 2025-10-15

arXiv ID: 2510.13206v1

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

摘要翻译

中文标题: 低温下Gross-Pitaevskii Gibbs测度的大偏差原理

摘要翻译: 我们证明了在低温区域下，与聚焦型Gross-Pitaevskii方程相关的条件Gibbs测度满足大偏差原理。该条件测度具有混合型结构：在能量上为正则系综，在粒子数上为微正则系综。特别地，我们的结果将Ellis、Jordan、Otto与Turkington所研究的混合系综的大偏差原理推广至更为奇异的情形，其中相互作用势为无界，且条件事件涉及发散的重整化常数。作为大偏差原理的推论，在低温极限下，Gibbs测度集中于孤子流形之上。

11. The local well-posedness, global existence and ill-posedness for the fifth order Camassa-Holm model

五阶Camassa-Holm模型的局部适定性、全局存在性与不适定性

作者: Xiaoxin Chen, Zhaoyang Yin

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13165v1

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

摘要翻译

中文标题: 五阶Camassa-Holm模型的局部适定性、全局存在性与不适定性

摘要翻译: 本文研究五阶Camassa-Holm模型。首先，我们改进了文献\cite{TangLiu2015,FOCH2021}中关于局部适定性的结果；其次，我们给出了爆破准则及全局存在性的充分条件；最后，当模型参数$b=\frac{5}{3}$时，我们证明了该模型在Besov空间$B^1_{\infty,1}$以及$B^{\frac{3}{2}}_{2,q}$（其中$q\in(1,+\infty]$）中以范数膨胀意义下是不适定的。

12. The local well-posedness, blow-up phenomena and ill-posedness of a new fifth-order Camassa-Holm type equation

一类新的五阶Camassa-Holm型方程的局部适定性、爆破现象与不适定性

作者: Yiyao Lian, Zhaoyang Yin

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13164v1

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

摘要翻译

中文标题: 一类新的五阶Camassa-Holm型方程的局部适定性、爆破现象与不适定性

摘要翻译: 本文研究了由李志（Li \cite{Li.Z}）推导出的一类新的五阶Camassa-Holm型方程。首先，我们在Besov空间中建立了该方程Cauchy问题在Hadamard意义下的局部适定性。其次，我们获得了爆破准则；在此基础上，通过利用守恒律并建立局部有界性，我们推导出一个精确确定爆破时间的爆破结果。最后，通过范数膨胀论证，我们证明了该方程在临界Sobolev空间 $H^{\frac{1}{2}}$ 中的不适定性。

13. Weak solutions to the parabolic $p$-Laplace equation in a moving domain under a Neumann type boundary condition

移动域上满足Neumann型边界条件的抛物型$p$-Laplace方程的弱解

作者: Tatsu-Hiko Miura

类别: math.AP, 35K20, 35K92, 35R37

发布日期: 2025-05-19

arXiv ID: 2505.12598v2

链接: http://arxiv.org/abs/2505.12598v2

摘要翻译

中文标题: 移动域上满足Neumann型边界条件的抛物型$p$-Laplace方程的弱解

摘要翻译: 本文研究了在移动域上、满足与总质量守恒相对应的Neumann型边界条件的抛物型$p$-Laplace方程（$p>2$）。我们通过在演化Bochner空间中运用Galerkin方法与单调性论证，建立了弱解的存在性与唯一性。主要困难在于刻画非线性梯度项的弱极限，其中需处理一个源自边界条件的项，该术语无法被纳入单调算子框架。为克服这一困难，我们在含时间依赖基函数的移动域上证明了统一时间的Friedrichs型不等式，并利用其推导出近似解的强收敛性。此外，当初始数据具有更高正则性时，我们进一步证明了时间导数在$L^2$意义下存在。

14. A Convergent Finite Element Scheme for the Q-Tensor Model of Liquid Crystals Subjected to an Electric Field

电场作用下液晶Q张量模型的收敛有限元格式

作者: Max Hirsch, Franziska Weber

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

发布日期: 2023-07-20

arXiv ID: 2307.11229v3

链接: http://arxiv.org/abs/2307.11229v3

摘要翻译

中文标题: 电场作用下液晶Q张量模型的收敛有限元格式

摘要翻译: 本文研究了在电场作用下的Landau-de Gennes Q张量液晶模型，并为其构建了一种全离散数值格式。该格式基于体势能的凸分解方法，并引入了一种针对Q张量的截断算子，以确保问题的适定性。我们证明了该数值格式的稳定性与适定性。进一步地，在对格式允许参数施加适当限制的条件下，我们证明了当时间步长与网格尺度趋于零时，全离散格式的解在某个子序列意义下收敛于Q张量模型的弱解。最后，我们给出了由该数值格式计算得到的数值结果，其中表明该格式能够有效模拟Freedericksz相变过程。

15. Dynamics of periodic fractional discrete NLS in the continuum limit

连续极限下周期性分数阶离散非线性薛定谔方程的动力学

作者: Brian Choi

类别: math.AP, cs.NA, math.CA, math.NA, 35Q55, 81T27, 65M06, 34A34, 37K60

发布日期: 2024-01-24

arXiv ID: 2401.13152v2

链接: http://arxiv.org/abs/2401.13152v2

摘要翻译

中文标题: 连续极限下周期性分数阶离散非线性薛定谔方程的动力学

摘要翻译: 本文从分析与动力学的角度，通过调节网格尺寸 $h>0$ 与非局部Lévy指数 $\alpha \in (0,2]$，研究了周期性格点上的分数阶离散非线性薛定谔方程（fDNLS）。我们证明，当 $h \rightarrow 0$ 时，该离散系统在能量空间以下的范数意义下收敛于分数阶非线性薛定谔方程（fNLS），其收敛性通过直接估计离散解与连续解在环面 $L^2(\mathbb{T})$ 上的差值，并借助离散周期性Strichartz估计得以建立。通过有限差分法（FDM）所获得的收敛速率在能量空间中被严格确定为 $O(h^{\frac{\alpha}{2+\alpha}})$。为进一步揭示fDNLS的收敛行为，我们系统考察了连续波（CW）解在调制不稳定性背景下的各类动力学特性，特别强调了由非局部格点耦合所表征的线性色散（或晶格衍射）与非线性效应之间的相互作用。尤其地，我们通过解析与数值方法共同证实：当 $h \rightarrow 0$ 时，CW解的最大增益 $\Omega_m$ 对其振幅 $A$ 的依赖关系由线性过渡为二次依赖。

不相关论文（基本信息）

1. Sharp decay estimates and numerical analysis for weakly coupled systems of two subdiffusion equations

两类亚扩散方程弱耦合系统的精确衰减估计与数值分析

作者: Zhiyuan Li, Yikan Liu, Kazuma Wada

类别: math.AP, cs.NA, math.NA, 35R11, 35K51, 35B40, 65M06

发布日期: 2025-04-25

arXiv ID: 2504.18295v2

链接: http://arxiv.org/abs/2504.18295v2

摘要: This paper investigates the initial-boundary value problem for weakly coupled systems of time-fractional subdiffusion equations with spatially and temporally varying coupling coefficients. By combining the energy method with the coercivity of fractional derivatives, we convert the original partial differential equations into a coupled ordinary differential system. Through Laplace transform and maximum principle arguments, we reveal a dichotomy in decay behavior: When the highest fractional order is less than one, solutions exhibit sublinear decay, whereas systems with the highest order equal to one demonstrate a distinct superlinear decay pattern. This phenomenon fundamentally distinguishes coupled systems from single fractional diffusion equations, where such accelerated superlinear decay never occurs. Numerical experiments employing finite difference methods and implicit discretization schemes validate the theoretical findings.

2. Liouville properties for differential inequalities with $(p,q)$ Laplacian operator

关于具有 $(p,q)$-Laplace 算子的微分不等式的刘维尔性质

作者: Mousomi Bhakta, Anup Biswas, Roberta Filippucci

类别: math.AP, 35J60, 35J92, 35B08, 35J70, 35A01

发布日期: 2025-10-15

arXiv ID: 2510.13576v1

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

摘要: In this paper, we establish several Liouville-type theorems for a class of nonhomogenenous quasilinear inequalities. In the first part, we prove various Liouville results associated with nonnegative solutions to \begin{equation*}\tag{$P_s$} -\Delta_p u-\Delta_q u\geq u^{s-1} , \text{ in }, \Omega, \end{equation*} where $1<q<p$, $s>1$ and $\Omega$ is any exterior domain of $\mathbb{R}^N$. In particular, we prove that for $q<N$, inequality $(P_s)$ does not admit any positive solution when $s<q_*$ and $(P_s)$ admits a positive solution if $s>q_*$, where $q_*=\frac{q(N-1)}{N-q}$ is the Serrin exponent for the $q$-Laplacian. Further, we show that when $s=q_*$ and $p<s$ the only nonnegative solution to $(P_s)$ is the trivial solution. On the other hand, for $q\geq N$ we prove that $u\equiv 0$ is the only nonnegative solution for $(P_s)$ for any $s>1$. In the second part, we consider the inequality \begin{equation*}\tag{$P_{sm}$} -\Delta_p u-\Delta_q u \geq u^s |\nabla u|^m \quad \text{ in }\mathbb{R}^N, \end{equation*} where $1<q<p$, $N>q$ and $s, \, m\geq 0$. We prove that, for $\{0\leq m\leq q-1\}\cup\{m>p-1\}$, the only positive solution to $(P_{sm})$ is constant, provided $s(N-q)+m(N-1)<N(q-1)$. This, in particular, proves that if $\Omega=\mathbb{R}^N$ then any nonnegative solution to $(P_s)$ with $1<q<N$ and $1<s<q_*$ is the trivial solution. To prove Liouville in the range $0\leq m<q-1$, we first prove an almost optimal lower estimate of any nonnegative supersolution of $(P_{sm})$ and then leveraging this estimate we prove Liouville result. To the best of our knowledge, this technique is completely new and provides an alternative approach to the capacity method of Mitidieri-Pohozaev provided higher regularity is available.

3. Normalized solutions to polyharmonic equations with Hardy-type potentials and exponential critical nonlinearities

具有Hardy型势与指数临界非线性的多调和方程的归一化解

作者: Bartosz Bieganowski, Olímpio Hiroshi Miyagaki, Jacopo Schino

类别: math.AP

发布日期: 2024-10-08

arXiv ID: 2410.05885v2

链接: http://arxiv.org/abs/2410.05885v2

摘要: Via a constrained minimization, we find a solution $(\lambda,u)$ to the problem \begin{equation*} \begin{cases} (-\Delta)^m u+\frac{\mu}{|x|^{2m}}u + \lambda u = \eta u^3 + g(u)\ \int_{\mathbb{R}^{2m}} u^2 , dx = \rho \end{cases} \end{equation*} with $1 \le m \in \mathbb{N}$, $\mu,\eta \ge 0$, $\rho > 0$, and $g$ having exponential critical growth at infinity and mass supercritical growth at zero.

4. Convexity and rigidity of hypersurfaces in Cartan-Hadamard manifolds

卡坦-哈达玛流形中超曲面的凸性与刚性

作者: Mohammad Ghomi

类别: math.DG, math.AP, math.MG, Primary: 53C20, 58J05, Secondary: 53C44, 52A15

发布日期: 2023-08-29

arXiv ID: 2308.15454v3

链接: http://arxiv.org/abs/2308.15454v3

摘要: We show that in Cartan-Hadamard manifolds $M^n$, $n\geq 3$, closed infinitesimally convex hypersurfaces $\Gamma$ bound convex flat regions, if curvature of $M^n$ vanishes on tangent planes of $\Gamma$. This encompasses Chern-Lashof-Sacksteder characterization of compact convex hypersurfaces in Euclidean space, and some results of Greene-Wu-Gromov on rigidity of Cartan-Hadamard manifolds. It follows that closed simply connected surfaces in $M^3$ with minimal total absolute curvature bound Euclidean convex bodies, as stated by Gromov in 1985. The proofs employ the Gauss-Codazzi equations, a generalization of Schur comparison theorem to CAT($k$) spaces, and other techniques from Alexandrov geometry outlined by Petrunin.

5. What Uniqueness for the Holst-Nagy-Tsogtgerel–Maxwell Solutions to the Einstein Conformal Constraint Equations?

霍尔斯特-纳吉-措格特格勒-麦克斯韦解在爱因斯坦共形约束方程中的唯一性问题

作者: Romain Gicquaud

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

发布日期: 2024-01-14

arXiv ID: 2401.07225v2

链接: http://arxiv.org/abs/2401.07225v2

摘要: This paper addresses the issue of uniqueness of solutions in the conformal method for solving the constraint equations in general relativity with arbitrary mean curvature as developed initially by Holst, Nagy, Tsogtegerel and Maxwell. We show that the solution they construct is unique amongst those having volume below a certain threshold.

6. Recovery of a matrix valued potential for the wave equation on stationary spacetimes

平稳时空中波动方程的矩阵值势的恢复

作者: Spyridon Filippas, Lauri Oksanen, Miika Sarkkinen

类别: math.AP, math.DG, 35R30, 53C50, 53C65

发布日期: 2025-10-15

arXiv ID: 2510.13410v1

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

摘要: We study the problem of recovering a time dependent matrix valued potential on a globally hyperbolic manifold from the knowledge of the source to solution map of a wave equation including a connection 1-form term. We exhibit sufficient conditions for solving this inverse problem under the assumption that the the manifold is stationary and that the connection term is time independent. The proof is based on two ingredients. The first is reduction of the problem to the study of a non-Abelian light ray transform and holds assuming global hyperbolicity only. The second is the study of this transform and establishing a link with a Riemannian analogue.

7. A rigorous formulation of Density Functional Theory for spinless fermions in one dimension

一维无自旋费米子密度泛函理论的严格表述

作者: Thiago Carvalho Corso

类别: math-ph, math.AP, math.MP, math.SP, quant-ph, Primary: 35J10 Secondary:81Q05, 81V74, 46N50

发布日期: 2025-04-07

arXiv ID: 2504.05501v2

链接: http://arxiv.org/abs/2504.05501v2

摘要: In this paper, we present a completely rigorous formulation of Kohn-Sham density functional theory for spinless fermions living in one dimensional space. More precisely, we consider Schr"odinger operators of the form $H_N(v,w) = -\Delta + \sum_{i\neq j}^N w(x_i,x_j) + \sum_{j=1}^N v(x_i)$ acting on $\wedge^N \mathrm{L}^2([0,1])$, where the external and interaction potentials $v$ and $w$ belong to a suitable class of distributions. In this setting, we obtain a complete characterization of the set of pure-state $v$-representable densities on the interval. Then, we prove a Hohenberg-Kohn theorem that applies to the class of distributional potentials studied here. Lastly, we establish the differentiability of the exchange-correlation functional and therefore the existence of a unique exchange-correlation potential. We then combine these results to provide a rigorous formulation of the Kohn-Sham scheme. In particular, these results show that the Kohn-Sham scheme is rigorously exact in this setting.

8. Metaplectic operators with quasi-diagonal kernels

具有准对角核的元辛算子

作者: Gianluca Giacchi, Luigi Rodino

类别: math.AP, math.FA, math.SG, quant-ph, 22E46, 47G30, 47F05, 81S30

发布日期: 2025-02-07

arXiv ID: 2502.05089v3

链接: http://arxiv.org/abs/2502.05089v3

摘要: Metaplectic operators form a relevant class of operators appearing in different applications, in the present work we study their Schwartz kernels. Namely, diagonality of a kernel is defined by imposing rapid off-diagonal decay conditions, and quasi-diagonality by imposing the same conditions on the smoothing of the kernel through convolution with the Gaussian. Kernels of metaplectic operators are not diagonal. Nevertheless, as we shall prove, they are quasi-diagonal under suitable conditions. Motivation for our study comes from problems in time-frequency analysis, that we discuss in the last section.

9. Bifurcation and multiplicity results for critical Grushin-Choquard problems

关于临界Grushin-Choquard问题的分支与多重性结果

作者: Suman Kanungo, Pawan Kumar Mishra, Giovanni Molica Bisci

类别: math.AP, 2020 Mathematics Subject Classification. 35J70, 35H20, 35A15

发布日期: 2025-10-15

arXiv ID: 2510.13299v1

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

摘要: We consider the following nonlocal Br'ezis-Nirenberg type critical Choquard problem involving the Grushin operator \begin{equation*} \left{ \begin{aligned} -\Delta_\gamma & u =\lambda u + \left(\displaystyle\int_\Omega \frac{|u(w)|^{2^_{\gamma, \mu}}}{d(z-w)^\mu}dw\right) |u|^{2^_{\gamma, \mu}-2}u \quad &&\text{in} \ \Omega, u &= 0 \quad &&\text{on} , \partial \Omega, \end{aligned} \right. \end{equation*} where $\Omega$ is an open bounded domain in $\mathbb{R}^N$, with $N \geq 3$, and $\lambda >0$ is a parameter. Here, $\Delta_\gamma$ represents the Grushin operator, defined as

$$ \Delta_\gamma u(z) = \Delta_x u(z) +(1+\gamma)^2 |x|^{2\gamma} \Delta_y u(z), \quad \gamma \geq 0, $$

where $z=(x,y)\in \Omega \subset \mathbb{R}^m\times \mathbb{R}^n$, $m+n=N \geq 3$ and $2^*_{\gamma,\mu}= \frac{2N_\gamma-\mu}{N_\gamma-2}$ is the Sobolev critical exponent in the Hardy-Littlewood context with $N_\gamma= m+(1+\gamma)n$ is the homogeneous dimension associated to the Grushin operator and $0<\mu<N_\gamma$. The homogeneous norm related to the Grushin operator is denoted by $d(\cdot)$. In this article, we prove the existence of bifurcation from any eigenvalue $\lambda^*$ of $-\Delta_\gamma$ under Dirichlet boundary conditions. Furthermore, we show that in a suitable left neighborhood of $\lambda^*$, the number of nontrivial solutions to the problem is at least twice the multiplicity of $\lambda^*$.

10. Curvature penalization of strongly anisotropic interfaces models and their phase-field approximation

强各向异性界面模型的曲率惩罚及其相场逼近

作者: Jean-François Babadjian, Blanche Buet, Michael Goldman

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13275v1

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

摘要: This paper studies the effect of anisotropy on sharp or diffuse interfaces models. When the surface tension is a convex function of the normal to the interface, the anisotropy is said to be weak. This usually ensures the lower semicontinuity of the associated energy. If, however, the surface tension depends on the normal in a nonconvex way, this so-called strong anisotropy may lead to instabilities related to the lack of lower semicontinuity of the functional. We investigate the regularizing effects of adding a higher order term of Willmore type to the energy. We consider two types of problems. The first one is an anisotropic nonconvex generalization of the perimeter, and the second one is an anisotropic nonconvex Mumford-Shah functional. In both cases, lower semicontinuity properties of the energies with respect to a natural mode of convergence are established, as well as $\Gamma$-convergence type results by means of a phase field approximation. In comparison with related results for curvature dependent energies, one of the original aspects of our work is that, in the context of free discontinuity problems, we are able to consider singular structures such as crack-tips or multiple junctions.

11. Non variational type critical growth nonlocal system

非变分型临界增长非局部系统

作者: Ashutosh Dixit, Hichem Hajaiej, Tuhina Mukherjee

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13242v1

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

摘要: This study investigates the existence, uniqueness, and multiplicity of positive solutions for a system of fractional differential equations given by: \begin{equation*} (-\Delta)^{s_i} u_{i}+\lambda_{i} u_{i}=\sum_{j=1}^{n} \alpha_{i j}\left|u_{j}\right|^{q_{i j}}\left|u_{i}\right|^{p_{i j}-2} u_{i} , u_i\in {\mathscr{D}^{s_i,2}\left(\mathbb{R}^{N}\right)}, i=1,2,\cdots,n, \end{equation*} where $N>2s=\max\{2s_i\}$, $s_i\in(0,1)$, $n\geq 2$, $\lambda_{i} \geq 0$, $\alpha _{ij}>0$, $p_{ij}<2^{*}_{s}$, and $p_{ij}+q_{ij}=2^{*}_{s}=\min\{{\frac{2N}{N-2s_i}\}}$ for $i\neq j \in \{1,2,...,n\}$. $2^{*}_s$ called the fractional critical sobolev exponent and $2^{*}_s=2 N /(N-2s)$ for $N > 2s$ and $2^{*}_s=+\infty$ for $N=2s$ or $N<2s$. Our work establishes novel uniqueness and multiplicity results for positive solutions, applicable whether the system possesses a variational structure or not. We provide a comprehensive characterization of the exact number of positive solutions under specific parameter configurations. Our analysis shows that the positive solution set behaves differently across three distinct regimes: $p_{ij}<2$, $p_{ij}=2$, and $2<p_{ij}<2^{*}_{s}$.

12. Infinitely many solutions for the prescribed scalar curvature problem with volcano-like curvature

具火山状曲率的预定标量曲率问题的无穷多解

作者: Tuoxin Li, Juncheng Wei, Haidong Yang

类别: math.AP

发布日期: 2025-10-15

arXiv ID: 2510.13239v1

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

摘要: In this paper, we consider the following prescribed scalar curvature problem: \begin{equation*} -\Delta u = K(x) u^{\frac{n+2}{n-2}}, \quad u>0\quad\hbox{in}\quad \mathbb{R}^n, \quad u \in D^{1,2}(\mathbb{R}^n), \end{equation*} where $K(x)$ is a volcano-like positive function such that

$$ K(x)= K(r_0)- c_0 | |x|- r_0|^m + O( | |x|- r_0|^{m+\theta}),\quad r_0- \delta <|x| <r_0+\delta $$

with $K(r_0), c_0, \delta>0, \theta >2, \min \{\frac{n-2}{2}, 2\} < m< n-2$. We first prove the existence of infinitely many positive solutions. A consequence of our proof yields that the infinitely many solutions constructed in \cite{WY} are non-degenerate in the whole $D^{1, 2}(\mathbb{R}^{n})$ space. To our knowledge, it seems to be the first result of infinitely many solutions of prescribed scalar curvature problem when the potential function $K(x)$ is not radial. Our non-degeneracy results are also more complete and improve the result in \cite{GMPS}.

13. The Combinatorial Loewner Property and super-multiplicativity inequalities for symmetric self-similar metric spaces

组合洛wner性质与对称自相似度量空间的超乘性不等式

作者: Riku Anttila, Sylvester Eriksson-Bique

类别: math.MG, math.AP, 30L10, 20F65, 51F99, 53C23, 28A78

发布日期: 2024-08-28

arXiv ID: 2408.15692v2

链接: http://arxiv.org/abs/2408.15692v2

摘要: This paper introduces a general construction of self-similar metric spaces as limits of discrete graphs. Our framework produces many classical examples, such as the Sierpi'nski carpet and the higher dimensional Menger sponges, but also a rich class of new examples. The main result of the work roughly speaking states: If the construction is sufficiently symmetric then the limiting object supports useful moduli estimates, namely the Combinatorial Loewner property of Bourdon–Kleiner and the super-multiplicativity inequalities. The latter are established on Menger sponges for which it had not been previously known. The main new technique the work offers is a general framework of flows and resistance estimates.

14. Super Log-concavity of the First Eigenfunctions for Horo-convex Domains in Hyperbolic Space

双曲空间中水平凸域的第一特征函数的超对数凹性

作者: Guofang Wei, Ling Xiao

类别: math.AP, math.DG, 58, 35

发布日期: 2025-10-15

arXiv ID: 2510.13072v1

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

摘要: In this paper, we prove that the first eigenfunction of the Laplacian for a horo-convex domain $\Omega\subset\mathbb H^n$ is super log-concave when $\text{diam}(\Omega)$ is not large. Our result is optimal in the sense that there are counterexamples %are constructed for the cases when $\Omega$ is not horo-convex or when $\text{diam}(\Omega)$ is large respectively

15. Strong multiplicity one theorems and homological min-max theory

强单重性定理与同调极小-极大理论

作者: Adrian Chun-Pong Chu, Yangyang Li

类别: math.DG, math.AP

发布日期: 2023-09-14

arXiv ID: 2309.07741v3

链接: http://arxiv.org/abs/2309.07741v3

摘要: It was asked by Marques-Neves which min-max $p$-widths of the unit $3$-sphere lie strictly between $2\pi^2$ and $8\pi$. We show that the 10th to the 13th widths do. More generally, we prove stronger versions of X. Zhou’s multiplicity one theorem.

16. On operator fields in the upper triangular Toeplitz form

关于上三角Toeplitz形式的算子场

作者: M. M. Chernin, A. Yu. Konyaev

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

发布日期: 2025-10-14

arXiv ID: 2510.12923v1

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

摘要: In this work, we solve the fundamental problem of describing the coordinate transformations that preserve the upper triangular Toeplitz form of the given operator field. Surprisingly, this problem is closely related to the description of all Nijenhuis operators in the same form. This description, as well as the formulas for the aforementioned coordinate transformations, are given by the implicit formulas involving matrix-valued functions.