生成时间: 2025年12月05日 03:52:00

重点关注论文数量: 4

了解领域论文数量: 4

不相关论文数量: 18

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

1. Time-periodic solutions to an energy balance model coupled with an active fluid under arbitrary large forces

在任意大幅力作用下与活性流体耦合的能量平衡模型的时间周期解

作者: Gianmarco Del Sarto, Matthias Hieber, Filippo Palma, Tarek Zöchling

类别: math.AP, math-ph

发布日期: 2025-12-04

arXiv ID: 2512.04800v1

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

详细分析

中文标题: 在任意大幅力作用下与活性流体耦合的能量平衡模型的时间周期解

研究对象和背景

本文研究的是一个将二维 Sellers型能量平衡模型（Energy Balance Model, EBM）与三维 原始方程（Primitive Equations, PE）通过动态边界条件耦合而成的偏微分方程系统。该系统的物理背景源于气候建模，旨在描述地球大气与海洋之间大规模热力与流体动力学相互作用的动力学过程。

能量平衡模型是一类简化的气候模型，其核心思想是通过平衡地球接收到的太阳辐射与向外发射的长波辐射来模拟全球温度分布。尽管EBM在空间分辨率上远不如完整的一般环流模型（GCMs），但它们能够捕捉到气候系统中的关键反馈机制，例如冰反照率反馈（ice-albedo feedback），从而被广泛用于理论分析、长期气候变化预测以及临界点（tipping points）的研究。

本文所考虑的EBM属于Sellers类型，即其共反照率函数 $ \beta(\rho) $ 是Lipschitz连续的，具体形式为：

$$ \beta(\rho) = \beta_1 + (\beta_2 - \beta_1)\frac{1 + \tanh(\rho - \rho_{\text{ref}})}{2}, $$

其中 $ 0 < \beta_1 < \beta_2 $ 分别对应冰覆盖区与无冰区的反照率，$ \rho_{\text{ref}} \approx 263K $ 是相变参考温度。这种光滑过渡避免了Budyko型模型中因不连续反照率导致的非唯一性问题。

而原始方程则是从Navier-Stokes方程在静力平衡假设下推导出的大尺度地转流体运动方程，适用于描述大气和海洋中的大尺度流动。PE包含水平动量方程、垂向静力平衡、不可压缩性条件等，已被Lions、Temam和Wang等人系统研究。

本文的关键创新在于：将EBM作为表面温度演化方程，与内部流体的PE进行双向耦合，且这种耦合发生在边界上——即通过动态边界条件实现。这使得表面温度 $ \rho = T|_{\Gamma_u} $ 不再是被动标量，而是由一个独立的偏微分方程驱动，从而形成“活性流体”（active fluid）结构。

整个系统定义在柱形区域 $ \Omega = G \times (0,1) $，其中 $ G = (0,1)^2 $ 表示水平域，$ \Gamma_u = G \times \{1\} $ 为其上表面边界（海洋-大气界面）。设未知函数如下：

  • $ u = (v,w): (0,\tau) \times \Omega \to \mathbb{R}^3 $：海洋速度场，其中 $ v $ 为水平速度，$ w $ 为垂直速度；
  • $ p: (0,\tau) \times \Omega \to \mathbb{R} $：压力；
  • $ T: (0,\tau) \times \Omega \to \mathbb{R} $：三维温度场；
  • $ \rho = T|_{\Gamma_u}: (0,\tau) \times G \to \mathbb{R} $：表面温度；
  • $ \bar{v}(x,y) = \int_0^1 v(x,y,\xi)\,d\xi $：水平速度的垂直平均。

系统控制方程如下（已归一化参数并取 $ \varepsilon = 1 $）：

$$ \left\{ \begin{aligned} &\partial_t v + u \cdot \nabla v - \Delta v + \nabla_H p = f_1, && \text{in } (0,\tau) \times \Omega, \\ &\mathrm{div}_H\, \bar{v} = 0, && \text{in } (0,\tau) \times G, \\ &\partial_t T + u \cdot \nabla T - \Delta T = f_2, && \text{in } (0,\tau) \times \Omega, \\ &\partial_t \rho + \bar{v} \cdot \nabla_H \rho - \Delta_H \rho + (\partial_z T)|_{\Gamma_u} = R(x,\rho) + f_3, && \text{in } (0,\tau) \times G, \end{aligned} \right. $$

其中 $ R(x,\rho) = Q(x)\beta(\rho) - |\rho|^3\rho $ 包含太阳辐射输入 $ Q(x) > 0 $ 和Stefan-Boltzmann型辐射损失项。

边界条件包括：

  • 垂直速度与垂直导数在上下边界消失：$ (\partial_z v)|_{\Gamma_u \cup \Gamma_b} = w|_{\Gamma_u \cup \Gamma_b} = 0, (\partial_z T)|_{\Gamma_b} = 0 $；
  • 水平方向周期性边界条件，应用于 $ \partial G \times [0,1] $ 上所有变量；
  • 初始条件：$ v(0) = v_0, T(0) = T_0 $。

特别值得注意的是，表面温度 $ \rho $ 的演化由一个二维抛物型方程支配，并通过通量项 $ (\partial_z T)|_{\Gamma_u} $ 与三维温度场耦合，构成典型的动态边界条件。这类问题不能用传统方法处理，因为边界本身具有独立动力学。

为了数学处理方便，作者引入了压力分解：

$$ p(t,x,y,z) = p_s(t,x,y) - \theta(t,x,y,z), \quad \theta(t,x,y,z) = \int_0^z T(t,x,y,\xi)\,d\xi, $$

以及由不可压条件导出的垂直速度表达式：

$$ w(t,x,y,z) = -\int_0^z \mathrm{div}_H\, v(t,x,y,\xi)\,d\xi. $$

由此可将原系统重写为仅依赖于 $ v, T, \rho $ 的封闭系统。

函数空间方面，文章采用带有周期性边界的Sobolev空间：

  • $ H^{m,q}_{\text{per}}(\Omega) $：在水平方向具有周期性的 $ W^{m,q} $ 空间；
  • $ L^2_\sigma(\Omega) $：满足 $ \mathrm{div}_H v = 0 $ 的水平速度空间；
  • 解空间定义为插值空间：
    $$ X_w := (X_0, X_1)_{1/2,2}, \quad \text{其中 } X_0 = L^2_\sigma(\Omega;\mathbb{R}^2) \times L^2(\Omega) \times L^2(G), $$
    $$ X_1 = \{v \in H^2_{\text{per}} \cap L^2_\sigma : (\partial_z v)|_{\Gamma_u \cup \Gamma_b} = 0\} \times \{T \in H^2_{\text{per}} : (\partial_z T)|_{\Gamma_b} = 0\} \times H^2_{\text{per}}(G). $$

最终解空间为：

$$ E_{1,\tau} := H^1(0,\tau;X_0) \cap L^2(0,\tau;X_1), \quad E_{w,\tau} := C([0,\tau];X_w) \cap L^2(0,\tau;X_0). $$

这一设定反映了问题的正则性需求：强解需在时间和空间上均有足够光滑性，而弱解则允许更低的正则性，但仍保持时间连续性。

主要定理或主要结果

本文的核心成果是以下两个定理：

定理 2（存在性）
令 $ T^* > 0 $，外力 $ (f_1,f_2,f_3) \in E_{0,\tau} $ 对所有 $ \tau > 0 $ 成立，且为 $ T^* $-周期函数。则系统 (2.8) 在边界条件 (2.1) 下至少存在一个强 $ T^* $-周期解 $ (v,T,\rho) \in E_{1,\tau} $。

推论 3（稳态解存在性）
若外力 $ (f_1,f_2,f_3) \in L^2(X_0) $ 为常数，则系统存在至少一个强稳态解。

这些结论的关键突破在于：不要求外力满足任何小范数条件。相比之下，大多数关于Navier-Stokes或Keller-Segel等系统的周期解研究均需外力“足够小”以保证解的稳定性或使用压缩映射原理。而本文在任意大的周期外力下仍能证明强周期解的存在性，这是非常罕见且深刻的。

此外，文中还建立了以下辅助结果：

命题 4（全局强解存在唯一性）
对任意初值 $ (v_0,T_0,\rho_0) \in X_w $ 和外力 $ (f_1,f_2,f_3) \in E_{0,\tau} $，初边值问题存在唯一的全局强解 $ (v,T,\rho) \in E_{1,\tau} $。

命题 5（弱周期解存在性）
对任意大的 $ T^* $-周期外力，系统存在至少一个弱 $ T^* $-周期解。

最后通过弱-强唯一性原理，将弱周期解提升为强解，完成主定理的证明。

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

a) 核心分析工具

  1. Galerkin逼近法与Brouwer不动点定理
    文章首先构造有限维近似解序列 $ (v_m, T_m, \rho_m) $，选取 $ X_1 $ 中的一组正交基 $ \{\psi_k, \tilde{\psi}_k\} $，将解投影到前 $ m $ 个模态。所得系数满足常微分方程组，利用标准ODE理论得到局部解。接着定义映射 $ S: c_m(0) \mapsto c_m(T^*) $，并通过能量估计得出一致有界性，从而在某个大球 $ B_R \subset \mathbb{R}^m $ 上应用Brouwer不动点定理，获得周期初始值 $ c_m(0) = c_m(T^*) $，进而延拓为全局周期解。

  2. 紧性与收敛性论证
    通过对近似解建立统一的能量估计（如(4.2)），结合Ascoli-Arzelà定理和Banach-Alaoglu定理，提取子列使得：

    $$ (v_m, T_m, \rho_m) \overset{*}{\rightharpoonup} (v,T,\rho) \quad \text{in } L^\infty(0,\tau;X_0), $$

    $$ (v_m, T_m, \rho_m) \rightharpoonup (v,T,\rho) \quad \text{in } L^2(0,\tau;X_1), $$

    并借助Friedrichs引理证明其在 $ E_{0,\tau} $ 中强收敛，从而可通过极限过程验证弱解满足积分形式与能量不等式。

  3. 最大正则性理论（Maximal $ L^2 $-regularity）
    在证明全局强解时，作者依赖于线性算子的最大正则性性质。对于形如 $ u' + Au = f $ 的抽象Cauchy问题，若 $ A $ 生成解析半群且满足适当谱条件，则解满足：

    $$ \|u'\|_{L^2} + \|Au\|_{L^2} \leq C\|f\|_{L^2}. $$

    这一工具用于控制高阶范数 $ \|\Delta v\|_{L^2} $ 等，确保解不会爆破。

  4. 弱-强唯一性原理
    设 $ (v_1,T_1,\rho_1) $ 为弱周期解，$ (v_2,T_2,\rho_2) $ 为以相同初值出发的强解。通过比较两者差异并利用能量估计，结合Gronwall不等式，证明二者几乎处处相等。此技术绕过了直接构造强解的困难，而是“先弱后强”。

b) 关键技术难点与处理方式

  1. 动态边界条件带来的非局部性
    边界方程涉及 $ (\partial_z T)|_{\Gamma_u} $，它依赖于内部温度梯度，使边界与体域深度耦合。作者通过引入积分项 $ \theta = \int_0^z T\,d\xi $ 将其转化为源项出现在动量方程中，从而将系统统一为全空间问题。

  2. 非线性项的控制
    关键非线性项包括：

    • 流体对流项 $ u \cdot \nabla v $：利用Sobolev嵌入 $ H^1 \hookrightarrow L^6 $ 控制；
    • 边界对流 $ \bar{v} \cdot \nabla_H \rho $：由于 $ \bar{v} \in H^2 $，可用 $ \|\bar{v}\|_{L^\infty} \leq C\|\bar{v}\|_{H^2} $ 处理；
    • 反应项 $ R(x,\rho) = Q\beta(\rho) - |\rho|^3\rho $：五次非线性 $ |\rho|^3\rho $ 提供天然耗散，有助于能量闭合。

  3. 无小力假设下的存在性突破
    传统方法往往依赖于外力小以构造压缩映射。本文完全避开此路径，转而采用：

    • 先证弱周期解存在（无需小力）；
    • 再证强解全局存在唯一（即使力大也可控）；
    • 最后用弱-强唯一性推出弱解实为强解。

这种方法摆脱了对非线性扰动强度的限制，具有更强普适性。

与之前工作的比较

a) 突破性改进

  1. 首次在动态边界条件下建立强周期解存在性
    此前关于EBM的研究（如Diaz & Tello [13], Badii & Diaz [3]）仅考虑无流体或静态边界情形。而本文首次将其与活性流体（primitive equations）动态耦合，并成功处理由此产生的高度非线性交互。

  2. 彻底移除外力的小性假设
    经典文献如Arendt-Bu [2]、Galdi-Sohr [20]等要求 $ \|f\| $ 足够小才能应用隐函数定理或摄动方法。而本文允许 $ f \in L^2 $ 任意大小，这是目前在类似PDE系统中极为少见的结果。尤其值得注意的是，即使是Navier-Stokes方程，也尚无在任意大幅周期力下存在强周期解的通用结论。

  3. 推广了解的概念框架
    以往对动态边界问题的研究多集中于弱解或能量解（energy solutions）。本文不仅构造了弱周期解，还通过弱-强唯一性获得了强解，显著提升了正则性水平。

  4. 为气候模型提供新的数学基础
    本工作为“气候系统在外强迫（如季节辐射）下趋于周期响应”的直观现象提供了严格的PDE解释。特别是当太阳辐射呈周期变化时，系统存在稳定的周期温度与流场结构，这对理解年际振荡、季风周期等具有重要意义。

综上所述，本文通过巧妙组合Galerkin方法、最大正则性与弱-强唯一性，在不要求任何小参数的前提下，解决了强非线性、高维度、多尺度耦合系统的时间周期解问题，代表了在非线性偏微分方程周期解理论上的重要进展。

2. Generalized Navier-Stokes equations, associated with the Dolbeault complex

与Dolbeault复形相关的广义Navier-Stokes方程

作者: Shlapunov Alexander, Polkovnikov Alexander

类别: math.AP, math.FA

发布日期: 2025-12-04

arXiv ID: 2512.04777v1

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

详细分析

中文标题: 与Dolbeault复形相关的广义Navier-Stokes方程

研究对象和背景

本文研究的是一类在复空间 $\mathbb{C}^n$（其中 $n > 1$）上的带时间变量的偏微分方程系统，其结构形式上类似于经典的不可压缩流体的 Navier-Stokes方程，但其基本算子并非来自欧几里得空间中的梯度、散度和旋度，而是由多维 Cauchy-Riemann算子 $\overline{\partial}$ 及其伴随算子 $\overline{\partial}^*$ 所生成。这类系统被置于一个称为 Dolbeault复形 的框架下进行分析。

具体地，作者考虑的是定义在时空区域 $\mathbb{C}^n \times [0, T]$ 上的 (0,q)-型外微分形式 $u(z,t)$ 和压力项 $p(z,t)$ 满足如下非线性发展系统：

$$ \begin{cases} \partial_t u + \mu \Delta_q u + N_q\left(\partial_q u, u\right) - \partial_{q-1} p = f & \text{in } \mathbb{C}^n \times (0,T), \\ (\partial_{q-1})^* u = 0 & \text{in } \mathbb{C}^n \times (0,T), \\ u(\cdot, 0) = u_0 & \text{on } \mathbb{C}^n, \end{cases} \tag{1} $$

其中：

  • $\mu > 0$ 是粘性常数；
  • $f$ 是给定的外部力项，为 $(0,q)$-形式；
  • $u_0$ 是初始速度场，也为 $(0,q)$-形式；
  • $\Delta_q = (\partial_q)^* \partial_q + \partial_{q-1} (\partial_{q-1})^*$ 是 Dolbeault-Laplacian算子，即第 $q$ 阶的 广义Laplace算子，具有强椭圆性；
  • $N_q$ 是一个非线性双线性微分算子，表示对流项，其构造基于两个固定零阶常系数双线性算子 $M_q^1$ 和 $M_q^2$，定义为：
    $$ N_q u = M_q^1(\partial_q u, u) + \partial_{q-1} M_q^2(u, u). $$
  • 压力项 $p$ 是一个 $(0,q{-}1)$-形式，满足兼容条件 $(\partial_{q-2})^* p = 0$；
  • $\partial_t$ 表示关于时间的偏导数；
  • 所有形式的系数依赖于参数 $t \in [0,T]$。

该系统的“不可压缩”条件体现为 $(\partial_{q-1})^* u = 0$，这对应于经典Navier-Stokes中 $\mathrm{div}\,u = 0$ 的角色，但在Dolbeault复形中，这是通过 共轭外微分算子 来实现的。

值得注意的是，当 $q=1$ 时，存在一种自然的非线性选择：

$$ N_q u = \star(\star_1 u \wedge u) + \partial_0 |u|^2, $$

它在结构上类似于Navier-Stokes方程中的 Lamb形式非线性项 $\mathrm{rot}\,u \times u$，并且满足关键的能量守恒条件：

$$ \left(M_q^1(\partial_q w, v), v\right)_{L^2} = 0, \quad \forall v \in V_q, $$

其中 $V_q$ 是光滑 $(0,q)$-形式空间中满足 $(\partial_{q-1})^* u = 0$ 的子空间。这一条件是证明弱解存在的核心假设之一。

研究区域为全空间 $\mathbb{C}^n \cong \mathbb{R}^{2n}$，无边界，因此不涉及传统意义上的边值问题。然而，由于空间是非紧的，解的行为在无穷远处需加以控制，例如通过引入加权积分或局部有界性条件。文中特别强调了使用 Bochner-Sobolev空间 来处理时间和空间正则性的耦合。

这些函数空间的设计极具特色。对于整数 $s,k \in \mathbb{Z}_+$，定义三类空间：

  1. 速度空间 $B^{k,2s,s}_{\mathrm{vel},q}$：包含所有满足一定时空正则性的 $(0,q)$-形式 $u$，使得

    $$ u \in C\left(I, H^{k+2s}_q\right) \cap L^2\left(I, H^{k+1+2s}_q\right), $$

    并且所有满足 $|\alpha| + 2j \leq 2s$ 的混合导数 $\partial_x^\alpha \partial_t^j u$ 属于适当的Sobolev类。

  2. 外力空间 $B^{k,2s,s}_{\mathrm{for},q}$：类似定义，用于容纳右端项 $f$。

  3. 压力空间 $B^{k+1,2s,s}_{\mathrm{pre},q-1}$：更精细，要求压力 $p$ 满足平均值归零条件：

    $$ \int_{\mathbb{R}^{2n}} p_J(x) h_0(x)\,dx = 0, \quad \forall t \in [0,T],\ #J = q-1, $$

    其中 $h_0 \in C^\infty(\mathbb{C}^n)$ 是固定满足 $\int h_0 dx = 1$ 的截断函数，以消除常数自由度，确保唯一性。

此外，还引入了 Galerkin逼近方法、能量估计 和 投影算子 $P_q$（基于Dolbeault-Laplacian的基本解 $\varphi_q$ 构造）来处理非线性和压力项的分离。

综上所述，本研究将经典Navier-Stokes理论推广至复几何背景下的 Dolbeault复形，利用其内在的椭圆-抛物结构，在全空间上建立了一套完整的适定性理论框架，重点在于弱解的存在性、强解的正则性准则以及整体映射的开性（稳定性）。

主要定理或主要结果

定理 1（弱解的存在性）

设 $s \in \mathbb{N}$ 且 $s \leq n \leq s+1$。若非线性项满足：

$$ \left(M_q^1(\partial_q w, v), v\right)_{L^2} = 0, \quad \forall v \in V_q, $$

则对任意 $(f, u_0) \in L^2\left(I, (H^1_q)'\right) \times H^0_q$，存在一个 $(0,q)$-形式 $u \in L^\infty(I, H^0_q) \cap L^2(I, H^1_q)$，满足以下弱形式：

$$ \frac{d}{dt}(u, v)_{L^2} + \mu (\partial_q u, \partial_q v)_{L^2} = \langle f - N_q u, v \rangle, \quad \forall v \in H^s_q, $$

且初值条件 $u(\cdot, 0) = u_0$ 成立。此外，$\partial_t u \in L^{n+2/(1-s)}\left(I, (H^s_q)'\right)$。

定理 2（强解的正则性准则）

在定理1条件下，若解 $u$ 进一步属于 $L^s(I, L^r_q)$，其中 $r > 2n$，$s$ 满足

$$ \frac{2}{s} + \frac{2n}{r} = 1, $$

则当右端项 $f$ 光滑时，该解为光滑解。

定理 3（线性化系统的可解性）

设 $n \geq 2$, $0 \leq q < n$, $s \in \mathbb{N}$, $k \in \mathbb{Z}_+$, $2s + k > n$，且 $w \in B^{k,2s,s}_{\mathrm{vel},q}$。则线性化问题：

$$ \partial_t u + \mu \Delta_q u + B_q(w, u) - \partial_{q-1} p = f, \quad (\partial_{q-1})^* u = 0 $$

诱导出一个连续的双射线性映射：

$$ A^w_q : B^{k,2s,s}_{\mathrm{vel},q} \times B^{k+1,2(s-1),s-1}_{\mathrm{pre},q-1} \to B^{k,2(s-1),s-1}_{\mathrm{for},q} \times H^{2s+k}_q, $$

并具有连续逆。

定理 4（非线性映射的开性）

在相同条件下，原非线性问题 (1) 诱导出一个连续、单射的非线性映射：

$$ A_q : B^{k,2s,s}_{\mathrm{vel},q} \times B^{k+1,2(s-1),s-1}_{\mathrm{pre},q-1} \to B^{k,2(s-1),s-1}_{\mathrm{for},q} \times H^{2s+k}_q, $$

该映射是 开映射，即其像集在拓扑意义下是开的，从而保证了解的局部稳定性。

定理 5（满射性判据）

上述映射 $A_q$ 是满射当且仅当：对任意子集 $S = S_{\mathrm{vel}} \times S_{\mathrm{pre}} \subset B^{k,2s,s}_{\mathrm{vel},q} \times B^{k+1,2(s-1),s-1}_{\mathrm{pre},q-1}$，若其像 $A_q(S)$ 在目标空间中预紧，则 $S_{\mathrm{vel}}$ 在 $L^s(I, L^r_q)$ 中有界，其中 $r,s$ 满足 $\frac{2}{s} + \frac{2n}{r} = 1$。

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

a) 核心分析学工具

本文采用了多个现代偏微分方程中的经典与前沿工具。首先是 Galerkin逼近法，这是处理非线性演化方程的标准手段。通过选取 $H^s_q$ 空间的一组正交基，将无限维问题投影到有限维子空间中求解近似解序列 $\{u_m\}$，再利用紧性论证提取收敛子列。为此，需要建立一致先验估计，尤其是 能量估计。

能量估计的形式为：

$$ \frac{d}{dt} \|u\|_{L^2}^2 + 2\mu \|\partial_q u\|_{L^2}^2 = 2\langle f - N_q u, u \rangle. $$

关键在于利用非线性项的反对称性质（由条件 $\left(M_q^1(\partial_q w, v), v\right)=0$ 保证），使其在内积中消失或可控，从而得到：

$$ \frac{d}{dt} \|u\|_{L^2}^2 + \mu \|\partial_q u\|_{L^2}^2 \leq C\|f\|^2, $$

结合 Grönwall引理 即可获得解的一致有界性。

其次，为了控制高阶导数，作者引入了精心构造的 Bochner-Sobolev空间族 $B^{k,2s,s}$，它们综合了时间与空间的正则性，并允许不同阶导数之间的插值。这类空间的嵌入关系依赖于 Gagliardo-Nirenberg不等式，例如：

$$ \|D^j u\|_{L^p} \leq C \|D^m u\|_{L^r}^\theta \|u\|_{L^q}^{1-\theta}, $$

其中参数满足适当的关系。这些不等式用于估计非线性项 $N_q u$ 的范数，确保其落在合适的对偶空间中。

另一个核心技术是 投影算子 $P_q$ 的构造。定义为：

$$ P_q = \varphi_q (\partial_q)^* \partial_q, $$

其中 $\varphi_q$ 是 $\Delta_q$ 的基本解。此算子将外力项 $f$ 投影到 $\ker((\partial_{q-1})^*)$ 上，从而分离出有效部分，排除掉由压力梯度引起的无散分量。这一操作类似于经典Stokes问题中的 Helmholtz分解，但在Dolbeault复形中，它是通过伪微分算子实现的。

最后，为了证明映射的开性，作者使用了 Banach空间中的隐函数定理。这要求验证Fréchet导数的可逆性，而这正是定理3的内容——线性化算子 $A^w_q$ 是连续双射，因而其逆连续（由开映射定理）。由此推出原非线性映射在局部是开的。

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

首要挑战是如何处理非线性项 $N_q u$ 的奇异性与增长性。作者通过引入双线性算子 $M_q^1, M_q^2$ 并施加结构性条件（如能量守恒），使非线性项在弱解估计中“自我抵消”。此外，借助Gagliardo-Nirenberg不等式，可以将高阶范数用低阶范数插值，从而控制非线性项的增长速度。

其次，压力项 $p$ 不显式出现在主方程中，必须通过约束条件重构。这里的关键技巧是利用 引理3：若 $F \in B^{k,2(s-1),s-1}_{\mathrm{for},q}$ 且 $P_q F = 0$，则存在唯一的 $p \in B^{k+1,2(s-1),s-1}_{\mathrm{pre},q-1}$ 使得 $\partial_{q-1} p = F$。这实际上是Dolbeault复形上 Poincaré引理 的定量版本，依赖于基本解 $\varphi_q$ 的良好性质。

此外，为处理全空间上的无界性，作者采用 局部化函数 $h_0$ 来规范压力的平均值，避免常数模的出现，从而保证唯一性。同时，通过设定不同的范数条件（如 $L^2(I,C_b)$ 或 $C(I,C_b)$）来区分不同正则性层级下的压力行为。

最后，在证明满射性判据时，使用了 反证法+紧性论证：假设解序列无界，但数据收敛，则由像集预紧性导出矛盾；反之，若有界，则可通过弱紧性提取极限解，证明像集闭，进而结合开性得满射。

与之前工作的比较

本文的主要突破体现在三个方面：一是将Navier-Stokes型方程从 de Rham复形 推广至 Dolbeault复形，二是构建了一套适用于复结构的 Bochner-Sobolev空间体系，三是获得了比已有结果更强的 稳定性结论。

以往工作如 [18] 已在de Rham复形框架下研究过类似的广义Navier-Stokes系统，其算子为外微分 $d$ 与余微分 $d^*$，对应的Laplacian为Hodge-Laplacian。然而，Dolbeault复形在复流形上有更丰富的几何结构，且其Neumann问题具有 次椭圆性（subellipticity），而非完全椭圆，导致解空间更大（例如零阶上同调可能无限维），这增加了分析难度。

本文首次系统地将此类流体力学模型移植到复空间 $\mathbb{C}^n$ 上，揭示了即使在没有物理流体背景的情况下，Navier-Stokes的数学结构仍能在抽象复形中重现。这种推广不仅拓展了维度（从实 $n$ 到复 $n$，即实 $2n$ 维），更重要的是改变了基本算子的本质属性。

其次，在函数空间的选择上，本文提出的 $B^{k,2s,s}$ 类空间比传统的 $L^p(I,W^{k,q})$ 更精细，能够精确刻画时间与空间导数之间的平衡关系（通过权重 $|\alpha| + 2j \leq 2s$），更适合处理抛物型方程的正则性传播。相比之下，早期工作往往只关注低阶能量空间，难以获得高阶稳定性和满射性判据。

最后，本文不仅证明了弱解存在性（已有文献已知），还进一步建立了 非线性映射的开性 和 满射性判据，这是对经典Ladyzhenskaya理论的重要推进。这意味着解不仅存在，而且在小扰动下依然存在，具有鲁棒性。而此前多数结果仅限于局部适定性或唯一性缺失的弱解框架。

综上，本文通过深度融合复分析、微分复形与非线性PDE方法，开辟了Navier-Stokes理论的一个新方向，为未来研究更一般的椭圆-抛物系统提供了范式。

3. Optimal cost for the null controllability of the Stokes system with controls having $n-1$ components and applications

Stokes系统在具有$n-1$个分量的控制下的零可控性的最优代价及其应用

作者: Felipe W. Chaves-Silva, Diego A. Souza, Marcos G. Ferreira-Silva

类别: math.AP, math.OC

发布日期: 2025-12-04

arXiv ID: 2512.04721v1

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

详细分析

中文标题: Stokes系统在具有$n-1$个分量的控制下的零可控性的最优代价及其应用

研究对象和背景

本文研究的是 $n$ 维（$n = 2,3$）有界光滑区域 $\Omega \subset \mathbb{R}^n$ 上的 Stokes系统 的 零可控性（null controllability）问题，其中控制作用仅施加于速度场的 $n-1$ 个分量上。该系统的数学形式如下：

$$ \begin{cases} u_t - \Delta u + \nabla p = f \mathbf{1}_\omega & \text{in } Q := \Omega \times (0,T), \\ \nabla \cdot u = 0 & \text{in } Q, \\ u = 0 & \text{on } \Sigma := \partial\Omega \times (0,T), \\ u(\cdot,0) = u_0 & \text{in } \Omega, \end{cases} $$

其中：

  • $u = u(x,t)$ 是流体的速度场，
  • $p = p(x,t)$ 是压力函数，
  • $f = (f_1,\dots,f_n)$ 是定义在控制区域 $\omega \subset \Omega$ 上的控制输入（或源项），
  • $\mathbf{1}_\omega$ 表示特征函数，表明控制只作用于子域 $\omega$。

此系统是经典的 Navier-Stokes 方程在线性化情形下的简化模型，适用于低雷诺数下粘性主导、惯性可忽略的不可压缩流体流动。

为了研究此类系统的可控性，首先需要引入适当的函数空间。文中定义了以下关键空间：

令 $L^2_{\text{div}} := \left\{ u \in L^2(\Omega)^n : \nabla \cdot u \in L^2(\Omega) \right\}$，

再设

$$ H := \left\{ u \in L^2_{\text{div}} : \nabla \cdot u = 0 \text{ in } \Omega,~ u \cdot \nu = 0 \text{ on } \partial\Omega \right\}, $$

以及

$$ V := H \cap H^1_0(\Omega)^n. $$

这些空间构成了分析 Stokes 系统的基础：$H$ 是满足不可压缩性和无穿透边界条件的速度场空间，而 $V$ 进一步要求速度属于 $H^1_0$，即具有更高正则性且在边界为零。

所谓零可控性，是指对任意初始状态 $u_0 \in H$，存在一个控制 $f \in L^2(\omega \times (0,T))^n$，使得对应的解 $u \in L^2((0,T);V) \cap C([0,T];H)$ 满足终值条件 $u(\cdot,T) = 0$ 在 $\Omega$ 中成立。

这一性质等价于对偶系统的可观测性不等式。具体地，考虑如下 对偶系统（adjoint system）：

$$ \begin{cases} z_t - \Delta z + \nabla \pi = 0 & \text{in } Q, \\ \nabla \cdot z = 0 & \text{in } Q, \\ z = 0 & \text{on } \Sigma, \\ z(\cdot,0) = z_0 & \text{in } \Omega. \end{cases} $$

则原系统的零可控性等价于存在常数 $C_{\text{obs}}(T) > 0$，使得对所有 $z_0 \in H$，其对应解 $z$ 满足：

$$ \|z(\cdot,T)\|_H^2 \leq C_{\text{obs}}(T) \int_0^T \int_\omega |z(x,t)|^2 dx dt. $$

然而，在本文中，由于控制仅作用于 $n-1$ 个分量（例如 $f_n \equiv 0$），因此对应的观测项也仅涉及 $z$ 的前 $n-1$ 个分量。于是所需的可观测性不等式变为：

$$ \|z(\cdot,T)\|_H^2 \leq C_{\text{obs}}(T) \int_0^T \int_\omega \sum_{k=1}^{n-1} |z_k(x,t)|^2 dx dt. $$

这带来了显著的技术挑战：如何从部分分量的信息中恢复整个向量场的能量？尤其是当缺少一个完整分量时，传统的 Carleman 估计方法可能不再直接适用。

此外，本文特别关注“控制代价”（cost of controllability）随时间 $T \to 0^+$ 的渐近行为。已知对于全控制（即所有 $n$ 个分量均可控）的情形，Stokes 系统的控制代价满足最优估计：

$$ C_{\text{obs}}(T) \leq C_1 e^{C_2 / T}, $$

这与热方程的零可控性代价一致，已被 G. Lebeau 和 F. W. Chaves-Silva [5] 利用谱方法证明。

但此前关于 $n-1$ 分量控制的工作（如 Fernández-Cara et al. [9], Coron–Guerrero [8]）虽然证明了零可控性，所得到的代价上界却非最优，例如形如 $e^{C/T^4}$ 或 $e^{C/T^9}$，远劣于 $e^{C/T}$。因此一个自然的问题是：

当控制减少至 $n-1$ 个分量时，零可控性的代价是否仍然保持为 $O(e^{C/T})$？

本文的核心贡献正是回答这一问题——答案是肯定的。

值得注意的是，根据 Lions 与 Zuazua [17] 的结果，在三维情况下，若只允许一个非零控制分量（即两个分量为零），一般无法实现零可控性。因此，$n-1$ 是实现全局零可控性的最小可能控制自由度数量。这也凸显了 $n-1$ 控制情形的重要性与临界性。

综上所述，本研究聚焦于：在不牺牲控制效率的前提下，能否用最少数量的控制分量驱动 Stokes 系统归零？ 特别是，即使去掉一个控制分量，控制代价的增长速率是否仍维持最优？

主要定理或主要结果

本文的主要结果包括以下几个定理，它们共同构成了对上述问题的完整解答。

定理 1.3（$n-1$ 分量控制下的零可控性及最优代价）

设 $T > 0$，$\omega \subset \Omega$ 为非空开集，$u_0 \in H$。则存在控制 $f = (f_1,\dots,f_{n-1},0) \in L^2(\omega \times (0,T))^n$，使得系统 (1.1) 的解满足 $u(\cdot,T) = 0$。进一步，存在仅依赖于 $\Omega$ 和 $\omega$ 而与 $T$ 无关的常数 $C_1,C_2 > 0$，使得控制满足：

$$ \left( \sum_{k=1}^{n-1} \|f_k\|_{L^2(\omega \times (0,T))}^2 \right)^{1/2} \leq C_1 e^{C_2 / T} \|u_0\|_H. $$

该定理表明，尽管控制少了一个分量，控制代价仍保持为 $O(e^{C/T})$，与全控制情形相同，从而实现了最优代价。

定理 1.4（对应对偶系统的可观测性不等式）

对任意 $z_0 \in H$，其对应的对偶解 $z = (z_1,\dots,z_n)$ 满足：

$$ \|z(\cdot,T)\|_H^2 \leq C_1 e^{C_2 / T} \int_0^T \int_\omega \sum_{k=1}^{n-1} |z_k(x,t)|^2 dx dt. $$

这是零可控性的对偶刻画，也是证明定理 1.3 的关键步骤。

定理 1.5（低频谱不等式）

存在常数 $M,K,\Lambda_1 > 0$，使得对任意复数序列 $(a_j) \in \ell^2$ 和任意 $\Lambda \geq \Lambda_1$，有：

$$ \sum_{\mu_j \leq \Lambda} |a_j|^2 \leq M e^{K \sqrt{\Lambda}} \int_\omega \sum_{k=1}^{n-1} \left| \sum_{\mu_j \leq \Lambda} a_j e_{j,k}(x) \right|^2 dx, $$

其中 $\{e_j\}_{j\geq1}$ 是 Stokes 算子的特征函数系，$e_j = (e_{j,1},\dots,e_{j,n})$，$\mu_j$ 为其对应的特征值。

这个谱不等式是整篇论文的技术核心，它说明：即使只观测 $n-1$ 个分量，也能以指数型代价控制整体能量。

定理 1.6（插值不等式）

设 $\hat{v}_\Lambda$ 如 (1.10) 所定义，即由低频特征函数构造的函数：

$$ \hat{v}_\Lambda(s,x) := \sum_{\mu_j \leq \Lambda} a_j \frac{\sinh(s \sqrt{\mu_j})}{\sqrt{\mu_j}} \Delta_x (e_{j,1}(x), \dots, e_{j,n-1}(x), 0), $$

并设 $h = \delta \Lambda^{-1/2}$（$\delta > 0$ 小）。则存在 $\delta_0, \Lambda_0, \mu, C, \rho \in (0,1)$，使得对任意 $\delta \in (0,\delta_0]$，存在 $C_\delta > 0$，对所有 $\Lambda \geq \Lambda_0$ 成立：

$$ \|\hat{v}_\Lambda\|_{H^1(W)^{n-1}} \leq C_\delta \left( e^{-\mu / h} \|\hat{v}_\Lambda\|_{H^1(Z)^{n-1}} + e^{C / h} \|\hat{v}_\Lambda\|_{H^{-1}_\rho(Z)^{n-1}}^\rho \| \partial_s \hat{v}_\Lambda(0,\cdot) \|_{L^2(\omega)^{n-1}}^{1-\rho} \right), $$

其中 $Z = (0,1) \times \Omega$, $W = (1/4,3/4) \times \Omega$。

此插值不等式通过半经典分析工具建立，是推导谱不等式的关键中间步骤。

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

a) 核心分析学工具

本文采用的主要方法是谱方法结合半经典分析（semiclassical analysis），这是处理线性偏微分方程中高频渐近行为和控制代价估计的强大框架。

1. 谱展开与谱不等式策略

受 Lebeau–Robbiano [16] 启发，作者将可控性问题转化为对 Stokes 算子低频部分的谱不等式的证明。基本思想是利用特征函数展开，并通过局部观测控制整体模长。这种策略避免了直接构造 Carleman 权重函数，尤其适合处理多分量耦合系统。

具体而言，通过构造辅助函数 $u_\Lambda(s,x)$，它是低频特征函数的双曲组合：

$$ u_\Lambda(s,x) = \sum_{\mu_j \leq \Lambda} a_j \frac{\sinh(s \sqrt{\mu_j})}{\sqrt{\mu_j}} e_j(x), $$

并验证其满足扩展的 Stokes 系统：

$$ \begin{cases} -\partial_s^2 u_\Lambda - \Delta_x u_\Lambda + \nabla_x p_\Lambda = 0 & \text{in } Z, \\ \nabla_x \cdot u_\Lambda = 0 & \text{in } Z, \end{cases} $$

从而可以借助椭圆或抛物型估计来获得能量控制。

2. 半经典参数与伪微分算子理论

由于压力项难以直接处理，作者引入半经典参数 $h = \delta \Lambda^{-1/2}$，将高频问题“缩放”为小 $h$ 极限下的问题。这样可使用半经典伪微分算子（$h$-pseudodifferential operators）进行局部化分析。

关键技巧在于：应用拉普拉斯算子 $\Delta_x$ 到前 $n-1$ 个分量上，构造新函数 $\hat{v}_\Lambda$，使其满足：

$$ \hat{v}_\Lambda = \Delta_x (\text{前 } n-1 \text{ 个分量}), $$

从而消去压力影响，并利用 Sobolev 范数等价性：

$$ \|\Delta_x u\|_{L^2} \sim \|u\|_{H^2}, \quad \forall u \in D(A), $$

实现从部分信息重建整体能量。

3. 局部 Carleman 估计

作者分别在内部点和边界点附近使用不同的 Carleman 不等式。

  • 对于内点，使用 Hörmander 子椭圆条件下的半经典 Carleman 估计（Theorem 2.8）；
  • 对于边界点，则构造合适的权重函数 $\phi(s,x) = e^{D(r(x) - (s-s_0)^2)}$，满足 $\partial_r \phi \neq 0$，并应用边界附近的 Carleman 估计（Theorem 2.10）。

这些估计允许将局部观测能量放大到整体范数。

4. 插值与迭代论证

通过插值不等式（Theorem 1.6），结合指数衰减项 $e^{-\mu/h}$ 与观测项之间的平衡，最终通过迭代推出谱不等式。

b) 技术性难题的克服

1. 压力项的消除

原始系统含压力 $\nabla p$，阻碍直接观测。传统做法是作用 $\nabla_x \Delta_x$ 消压，但作者指出只需作用 $\Delta_x$ 即可，因为：

$$ \|\Delta_x u\|_{L^2} \approx \|u\|_{H^2} $$

在 Stokes 算子定义域上成立。这简化了计算并保留了足够的正则性。

2. 缺失一个分量的信息补偿

由于仅观测 $n-1$ 个分量，必须利用不可压缩条件 $\nabla \cdot z = 0$ 进行补全。作者通过构造 $\hat{v}_\Lambda$ 并利用其满足的 PDE 结构，在局部坐标下结合边界的几何结构（如法向导数控制切向导数），实现信息传递。

3. 局部唯一延拓失效的应对

Remark 1.3 指出系统 (1.8) 不具备局部唯一延拓性质，故不能期望全局 Carleman 估计。为此，作者转而使用局部 Carleman 估计拼接的方法，在覆盖 $[0,1]\times\overline{\Omega}$ 的有限多个邻域中分别应用估计，再通过截断函数与紧性完成全局控制。

与之前工作的比较

a) 突破与改进

本文相比前人工作 [8][9][5] 实现了三方面重大突破：

  1. 首次获得 $n-1$ 控制下的最优控制代价 $O(e^{C/T})$，此前最好结果为 $e^{C/T^4}$（Fernández-Cara et al.）或更差；
  2. 完全去除对控制区域的几何限制，Coron–Guerrero [8] 虽然去除了接触边界的假设，但代价估计非最优；
  3. 揭示了控制自由度减少不影响代价阶数，说明 $e^{C/T}$ 是 Stokes 系统的本质代价尺度，不受控制分量数目的轻微变化影响。

b) 方法论优势

相较于依赖复杂 Carleman 权重构造的传统方法，本文采用谱+半经典方法更具结构性与可推广性。特别是通过引入辅助变量 $\hat{v}_\Lambda$，巧妙绕开了压力耦合难题，为未来研究其他流体系统（如 Navier-Stokes）的部分控制问题提供了新路径。

此外，所得结果在应用层面也有重要意义，例如在工程中设计低成本控制器时，允许关闭一个执行器而不显著增加能耗。

4. Diffusive limit of the Boltzmann equation around Rayleigh profile in the half space

半空间中绕Rayleigh分布的Boltzmann方程的扩散极限

作者: Hongxu Chen, Renjun Duan

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04403v1

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

详细分析

中文标题: 半空间中绕Rayleigh分布的Boltzmann方程的扩散极限

研究对象和背景

本文研究的是在半空间 $\Omega = \mathbb{T}^2 \times \mathbb{R}^+$ 中，关于 Boltzmann方程 在小Knudsen数 $\varepsilon > 0$ 下的 扩散极限（diffusive limit）问题。具体而言，考虑如下带边界的初边值问题：

$$ \frac{1}{\varepsilon}\partial_t F + v \cdot \nabla_x F = \frac{1}{\varepsilon^2} Q(F, F), \quad (t,x,v) \in (0,T] \times \Omega \times \mathbb{R}^3, $$

其中 $F = F(t,x,v)$ 是气体粒子的速度分布函数，$x = (x_1,x_2,x_3) \in \Omega = \mathbb{T}^2 \times \mathbb{R}^+$ 表示位置变量，$\mathbb{T}^2$ 为二维环面（即切向周期性），$v = (v_1,v_2,v_3) \in \mathbb{R}^3$ 为速度变量。

碰撞项 $Q(F,F)$ 采用硬球模型（hard sphere model）定义：

$$ Q(F,G)(v) = \int_{\mathbb{R}^3} \int_{S^2} |(v - u)\cdot \omega| \left[ F(u')G(v') - F(u)G(v) \right] d\omega du, $$

其中 $(v',u')$ 是弹性碰撞后的速度，满足动量与能量守恒：

$$ v' = v - [(v - u)\cdot \omega]\omega, \quad u' = u + [(v - u)\cdot \omega]\omega, \quad \omega \in S^2. $$

该系统处于 扩散尺度（diffusive scaling），即时间与空间被拉伸为 $t/\varepsilon$, $x/\varepsilon$，而碰撞项放大为 $1/\varepsilon^2$，这正是从 稀薄气体动力学 过渡到 不可压缩流体动力学 的典型标度。

边界条件设定在 $x_3 = 0$ 处，即物理壁面，采用 漫反射边界条件（diffuse reflection），且壁面以切向速度运动：

$$ F(t,x_1,x_2,0,v)\big|_{v_3 > 0} = M_w \int_{v'_3 < 0} F(t,x_1,x_2,0,v') |v'_3| dv', $$

其中 $M_w$ 是 壁面Maxwellian，由下式给出：

$$ M_w = \frac{1}{2\pi} \exp\left( -\frac{|v - u_w|^2}{2} \right), $$

而壁面速度 $u_w = (\varepsilon u_b, 0, 0)$ 与小参数 $\varepsilon$ 成正比，表示一个弱切向移动的边界。远场行为则趋于静止的全局Maxwellian：

$$ \mu(v) = \frac{1}{(2\pi)^{3/2}} e^{-|v|^2 / 2}, $$

并满足远场条件：

$$ F(t,x,v) \to \mu(v), \quad \text{当 } x_3 \to \infty. $$

这一设定对应经典的 Rayleigh问题：初始时刻流体静止于无限平板之上，突然启动平板使其沿切向以恒定速度滑动，从而激发剪切流动。在经典流体力学中，此问题由Navier-Stokes方程描述，其解称为 Rayleigh剖面（Rayleigh profile），表现为一个随时间扩散的粘性层。

然而，本文关注的是 稀薄气体 情形，即Knudsen数较小但非零时，如何从微观的Boltzmann方程推导出宏观的流体行为。目标是验证当 $\varepsilon \to 0$ 时，Boltzmann解收敛到以Rayleigh剖面为背景的不可压缩Navier-Stokes解。

值得注意的是，区域具有特殊结构：$\mathbb{T}^2$ 提供了切向周期性，避免无穷远处扰动传播的问题；而法向方向 $x_3 \in \mathbb{R}^+$ 构成半无限域，导致必须仔细处理边界效应，尤其是由于移动边界引起的非平衡边界层。

此外，函数空间的选择也极具挑战性。作者使用混合范数空间，包括 $L^2_{t,x,v}$、加权 $L^\infty$ 范数（如 $\|wR\|_{L^\infty}$，其中 $w(v) = e^{\beta |v|^2}, \beta \leq 1/8$）、以及边界积分空间 $L^2_{t,\gamma_\pm}$，这些空间用于控制余项 $R$ 及其时间导数 $\partial_t R$ 的行为。

特别地，边界投影算子 $P_\gamma$ 定义为：

$$ P_\gamma R := c_\mu \sqrt{\mu(v)} \int_{v'_3 < 0} R(t,x',0,v') \sqrt{\mu(v')} |v'_3| dv', \quad c_\mu = \sqrt{2\pi}, $$

它体现了漫反射边界对分布函数的影响，是非局部且非线性的关键部分。

综上所述，本文的研究对象是一个在 半空间中带有切向移动边界的非线性Boltzmann方程，在扩散极限下试图构造围绕 Rayleigh流体解 的渐近展开，并严格证明其有效性。核心难点在于：边界运动引入非齐次项、初始时刻可能存在的奇异性、以及高阶Hilbert展开中难以控制的非线性耦合。

主要定理或主要结果

论文的主要结果由以下定理给出：

Theorem 1. 设 $0 < \varepsilon \ll 1$，令非奇异的Rayleigh剖面 $u = (u_1,0,0)$ 如(1.20)所定义：

$$ u_1(t,x_3) = \frac{2u_b}{\sqrt{\pi}} \int_{\sqrt{x_3 / 4\kappa(t+\delta)}}^{\infty} e^{-r^2} dr, $$

并设 $f_1, f_2$ 分别由(1.14)和(1.22)给出：

$$ f_1 = u \cdot v \sqrt{\mu}, \quad f_2 = \frac{1}{3}|u|^2 \sqrt{\mu} - L^{-1} A : \nabla_x u + \frac{A : (u \otimes u)}{2}, $$

其中 $A_{ij}(v) = v_i v_j - \delta_{ij} |v|^2 / 3$ 为Burnett函数。

假设初始数据形如：

$$ F_0(x,v) = \mu + \varepsilon \sqrt{\mu} f_{1,0} + \varepsilon^2 \sqrt{\mu} f_{2,0} + \varepsilon^{3/2} \sqrt{\mu} R_0 \geq 0, $$

其中 $f_{1,0}, f_{2,0}$ 为 $f_1, f_2$ 在 $t=0$ 的取值。

则对于任意给定时间 $T > 0$，存在常数 $0 < \delta_0 \ll 1$ 满足 $\delta_0 T \ll 1$，使得若 $u_b > 0$ 满足 $|u_b| T^{1/2} < \delta_0$，且余项初值 $R_0$ 满足：

$$ \|R_0\|_{L^2_{x,v}} + \|\varepsilon^{1/2} w R_0\|_{L^\infty_{x,v}} + \varepsilon^{-1} \|(I-P)R_0\|_{L^2_{x,\nu}} + |\varepsilon^{-1/2}(I-P_\gamma)R_0|_{L^2_{\gamma_+}} + \|\partial_t R_0\|_{L^2_{x,v}} + \|\varepsilon^{3/2} w \partial_t R_0\|_{L^\infty_{x,v}} < \delta_0, $$

那么存在唯一的非负解

$$ F = \mu + \sqrt{\mu}(\varepsilon f_1 + \varepsilon^2 f_2) + \varepsilon^{3/2} \sqrt{\mu} R \geq 0 $$

定义在 $[0,T]$ 上，满足初边值问题 (1.1)-(1.3)，并且余项 $R$ 满足一致估计：

$$ \begin{aligned} &\|R\|_{L^\infty_t L^2_{x,v}} + \|\varepsilon^{1/2} w R\|_{L^\infty_{t,x,v}} + \varepsilon^{-1/2}|(I-P_\gamma)R|_{L^2_{t,\gamma_+}} \\ &\quad + \varepsilon^{-1}\|(I-P)R\|_{L^2_{t,x,\nu}} + \|\partial_t R\|_{L^\infty_t L^2_{x,v}} + \|\varepsilon^{3/2} w \partial_t R\|_{L^\infty_{t,x,v}} + \varepsilon^{-1}\|(I-P)\partial_t R\|_{L^2_{t,x,\nu}} < C \delta_0, \end{aligned} $$

其中 $C$ 为与 $T$ 无关的常数。

进一步地，有如下关于 $\varepsilon$ 的 $L^2$ 收敛性：

$$ \sup_{0 \leq t \leq T} \left\| \frac{F(t) - \mu}{\varepsilon \sqrt{\mu}} - f_1(t) \right\|_{L^2_{x,v}} < C \delta_0 \varepsilon^{1/2}. $$

该定理表明：在适当的小性条件下，可以构造出围绕Rayleigh剖面的Boltzmann解，且余项在多种范数下保持一致有界，并实现向不可压缩Navier-Stokes流的收敛。

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

a) 核心分析学工具

本文的核心技术建立在 Hilbert展开法 与 $L^2-L^6-L^\infty$ 混合估计框架 之上，结合了Guo等人发展的现代Boltzmann方程边值问题处理方法。

首先，采用经典的 Hilbert渐近展开：

$$ F = \mu + \sqrt{\mu}(\varepsilon f_1 + \varepsilon^2 f_2) + \varepsilon^{3/2} \sqrt{\mu} R, $$

将原方程代入后，通过比较各阶 $\varepsilon$ 的幂次，得到逐阶方程：

  • 零阶：$Lf_1 = 0$，说明 $f_1$ 属于碰撞算子的核空间；
  • 一阶：$v \cdot \nabla_x f_1 + L f_2 = \Gamma(f_1,f_1)$；
  • 二阶：$\mathcal{P}(\partial_t f_1 + v \cdot \nabla_x f_2) = 0$，其中 $\mathcal{P}$ 为宏观投影。

由此导出的宏观量满足 不可压缩Navier-Stokes方程：

$$ \begin{cases} \partial_t u + u \cdot \nabla_x u + \nabla_x p = \kappa \Delta u, & \nabla_x \cdot u = 0, \\ \nabla_x p = \nabla_x \left( \rho_2 + \theta_2 - \frac{1}{2}|u|^2 \right). \end{cases} $$

特别地，选取 $u$ 为Rayleigh剖面，即仅 $x_3$ 方向依赖的剪切流：

$$ u = (u_1(t,x_3), 0, 0), \quad u_1(t,x_3) = \frac{2u_b}{\sqrt{\pi}} \int_{\sqrt{x_3 / 4\kappa(t+\delta)}}^\infty e^{-r^2} dr, $$

它满足热传导方程 $\partial_t u_1 = \kappa \partial_{x_3 x_3} u_1$，从而自动满足上述NS方程（压力梯度为零）。这种选择极大简化了分析。

接下来，对余项 $R$ 建立控制方程：

$$ \partial_t R + \frac{1}{\varepsilon} v \cdot \nabla_x R + \frac{1}{\varepsilon^2} L R = h + \frac{1}{\varepsilon} \tilde{L} R + \varepsilon^{-1/2} \Gamma(R,R), $$

其中 $h = \varepsilon^{-1/2} h_1 + \varepsilon^{1/2} h_2$ 为源项，$\tilde{L} R = 2 \Gamma(f_1 + \varepsilon f_2, R)$ 为线性化非齐次项。

为了估计 $R$，作者采用了 三重范数结构：

  1. $L^2$ 能量估计：利用内积 $\langle R, \cdot \rangle$ 对方程进行测试，结合线性算子 $L$ 的谱性质（$L \geq \nu(v)(I-P)$），获得耗散项 $\|(I-P)R\|_{L^2_{t,x,\nu}}$ 的控制。
  2. $L^6$ 估计：借助Sobolev嵌入与速度平均引理，在空间上提升正则性，用于控制非线性项 $\Gamma(R,R)$。
  3. $L^\infty$ 加权估计：使用特征线法结合迭代技巧，建立对 $\|w R\|_{L^\infty}$ 的先验估计，其中权重 $w(v) = e^{\beta |v|^2}$ 控制高速粒子贡献。

此外，还对方程的时间导数 $\partial_t R$ 建立平行估计，以保证解的正则性，这对处理边界项和非线性耦合至关重要。

边界条件通过投影分解处理：

$$ R|_{v_3>0} = P_\gamma R + \cdots + r, $$

其中 $r$ 为由边界展开产生的非齐次项，其大小为 $O(\varepsilon^{1/2})$。通过对 $P_\gamma$ 和 $(I-P_\gamma)$ 的分别估计，结合边界迹不等式，实现对进出流的精细控制。

整个框架本质上是一种 perturbative argument：将Boltzmann解视为Maxwellian附近的扰动，借助线性算子的良好性质与非线性项的小性，通过压缩映射原理或连续性方法构造唯一解。

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

  1. 初始奇异性消除：经典Rayleigh问题中，若边界在 $t=0$ 突然启动，则流速梯度 $\partial_{x_3} u_1 \sim t^{-1/2}$ 发散。本文通过引入光滑化参数 $\delta > 0$，将初始条件设为：

    $$ u(t,x)|_{t=0} = \frac{2u_b}{\sqrt{\pi}} \int_{\sqrt{x_3 / 4\kappa \delta}}^\infty e^{-r^2} dr, $$

    从而避免 $t=0$ 处的奇异性，使所有导数有界。这是能进行高阶展开与余项估计的前提。

  2. 移动边界的处理：由于 $u_w = \varepsilon u_b$，墙Maxwellian $M_w$ 与全局平衡态 $\mu$ 不同。利用 Taylor展开引理（Lemma 1）：

    $$ \left| \mu^{-1/4} \left( M_w - c_\mu \mu \right) \right| \lesssim \varepsilon, \quad \left| \mu^{-1/4} \left( M_w - c_\mu \mu - \varepsilon u_b v_1 \mu \right) \right| \lesssim \varepsilon^2, $$

    实现对边界非齐次项 $r$ 的精确估计，确保其为高阶小量。

  3. 非线性项的控制：$\Gamma(R,R)$ 具有双线性结构，直接估计困难。作者利用已知的 非线性估计：

    $$ |\langle \Gamma(g,h), f \rangle| \lesssim \|g\|_{L^2_\nu} \|h\|_{L^2_\nu} \|f\|_{L^2_\nu}, $$

    并结合 $\varepsilon^{-1/2}$ 因子，要求余项足够小（$\delta_0$ 小），从而实现非线性项被线性耗散主导。

  4. 无稳态解的动态处理：由于Rayleigh剖面是自相似扩散解，不存在非平凡稳态。因此不能像Couette流那样研究稳态附近稳定性。本文转而研究 有限时间区间上的动态解，并通过控制 $|u_b| T^{1/2} < \delta_0$ 来限制总扰动积累。

  5. 省略Knudsen层修正：通常高阶展开需引入边界层校正 $f_b$ 以匹配边界条件。但作者发现，即使不引入 $f_b$，只要初始数据“well-prepared”，仍可在当前精度下完成估计。这简化了构造过程，但也限制了收敛阶数。

与之前工作的比较

a) 突破与改进

本文相对于前人工作有多方面重要推进：

  1. 首次严格证明Boltzmann方程绕Rayleigh剖面的扩散极限：尽管Sone [19] 和 Kuo [15,16] 对此问题进行了物理讨论与渐近分析，但缺乏严格的数学存在性与收敛性证明。本文首次在 全非线性、非稳态、带边情形 下完成了这一极限的严格验证。

  2. 克服移动边界带来的非平衡效应：多数已有工作（如[9][11]）考虑的是静止边界或仅温度变化。本文明确处理了 切向移动边界 导致的非零滑移速度，并将其纳入Hilbert展开框架，展示了如何协调边界Maxwellian与远场平衡态之间的差异。

  3. 去除初始奇异性：Maekawa [17] 在Navier-Stokes层面研究了带奇性的Rayleigh层稳定性，但需Gevrey类正则性。本文通过引入平滑初值（$\delta > 0$），在 Sobolev空间 框架内解决了该问题，降低了对初始数据的要求，使结果更具普适性。

  4. 避免高阶展开复杂性：相比传统需要Knudsen层校正的方法（如[20][21]），本文仅用两阶展开即达成目的，简化了构造流程。虽然牺牲了更高阶收敛性，但在实际应用中已足够，且更利于后续推广。

  5. 统一处理时间和空间变量：不同于仅局部存在性结果（如Kuo [15]），本文实现了 任意有限时间 $T$ 上的整体存在性（前提是 $u_b T^{1/2}$ 小），增强了结果的实用性。

  6. 清晰揭示约束来源：文中明确指出两个关键小性条件的物理含义：

    • $|u_b| T^{1/2} < \delta_0$：防止边界运动引起过强扰动；
    • $\delta_0 T \ll 1$：源于 $L^2-L^\infty$ 方法的能量累积效应，未来可通过 $L^3$ 方法改进。

综上，本文不仅填补了稀薄气体中Rayleigh问题的理论空白，也为研究其他非平衡边界驱动流动提供了新范式。

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

1. Fractured Poroelastic Media in the Limit of Vanishing Aperture

裂缝孔弹性介质在裂缝张开度趋于零极限情形下的渐近分析

作者: Maximilian Hörl, Kundan Kumar, Christian Rohde

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04978v1

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

摘要翻译

中文标题: 裂缝孔弹性介质在裂缝张开度趋于零极限情形下的渐近分析

摘要翻译: 本文考虑一类含有薄层非均质结构（即裂缝）的孔弹性介质。介质整体区域及裂缝内部的流体流动与力学变形均由准静态Biot方程所支配。假设裂缝的材料参数（如水力传导系数和弹性模量）按裂缝宽长比 $\varepsilon$ 的幂次进行尺度变换。基于先验估计，我们严格推导了当 $\varepsilon \rightarrow 0$ 时的极限模型，并识别出不同的极限情形。对于水力传导系数，我们得到了五种不同的极限情形；对于弹性性质，则得到两种情形。在多数情况下，极限模型退化为离散裂缝模型，而在其他情形下，则演化为由法向流动或法向变形主导的双尺度极限问题。

2. On hyperbolic approximations for a class of dispersive and diffusive-dispersive equations

关于一类色散型及扩散-色散耦合方程的双曲近似

作者: Rahul Barthwal, Firas Dhaouadi, Christian Rohde

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04882v1

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

摘要翻译

中文标题: 关于一类色散型及扩散-色散耦合方程的双曲近似

摘要翻译: 我们提出了一类具有非线性通量的色散型与扩散-色散耦合方程的新型近似系统。对于纯色散方程，我们构造了一个一阶严格双曲型近似系统，并通过构造适用于任意光滑通量的唯一对称化矩阵，建立了光滑解的局部适定性。在对通量施加更强条件的情况下，我们为该双曲系统提供了一个严格凸熵函数，其对应于原始色散方程的能量结构。为了近似扩散-色散方程，我们采用了一种粘弹性阻尼系统，该系统与色散演化方程双曲近似的所获熵结构相容。对于由此得到的双曲-抛物耦合近似系统，我们证明了其整体适定性。借助相对熵方法 \cite{dafermos2005hyperbolic}，我们证明了近似系统的解收敛到原方程的解。新近似系统的结构使得我们可以应用双曲守恒律领域中的标准数值模拟方法。我们在一系列涵盖不同目标方程的测试案例中验证了近似方法的收敛性，且其有效性超出了理论结果所限定的适用范围。此外，我们展示了该方法在强非线性效应下的适用性，包括导致振荡行为或激波层形成的复杂动力学现象。

3. Homogenized limits of Stokes flow and advective transport in thin perforated domains

薄型穿孔区域中 Stokes 流动与对流输运问题的均匀化极限

作者: Markus Gahn, Vlad Revnic

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04782v1

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

摘要翻译

中文标题: 薄型穿孔区域中 Stokes 流动与对流输运问题的均匀化极限

摘要翻译: 本文严格研究了定义在厚度为 $ \varepsilon^\alpha $ 阶（其中 $ \alpha \in (0,1) $）的薄 $ \varepsilon $-周期穿孔层中的流动与输运问题的均匀化及降维过程，因此该薄层的厚度相对于其孔隙率而言较大。目标是当 $ \varepsilon \to 0 $ 时，即薄层厚度趋于零的情况下，推导出相应的有效模型。对于流动问题，我们考虑在层的顶部和底部具有压力边界条件的不可压缩 Stokes 方程；而输运问题则由一个反应-扩散-对流方程描述，其中对流项由 Stokes 模型提供的流体速度驱动，并针对水平方向的不同扩散系数进行多种尺度分析，以刻画低扩散与高扩散情形。在极限情况下，Stokes 流动呈现出一种达西型定律，其达西速度仅依赖于达西压力在垂直方向上的导数。输运问题的有效方程仍为扩散-对流型，包含均匀化后的系数，其对流速度由达西速度给出，且对流仅发生在垂直方向。当垂直方向扩散较慢时，有效扩散仅出现在垂直方向；而在水平方向具有强扩散的情形下，有效扩散则在所有空间方向上均存在。为了实现极限过程，我们采用了适用于当前微观几何结构的双尺度收敛方法，该方法基于一致的先验估计。在宏观模型推导中的关键难点包括对流体压力的控制——为此我们为薄穿孔区域构造了一个 Bogovskiǐ 算子，以及输运方程微观解的强双尺度收敛性，后者对于处理对流项的极限过程至关重要。

4. Mixing at the Batchelor Scale for White-In-Time Flows

关于白噪声时间流场在Batchelor尺度下的混合性研究

作者: Robin Chemnitz, Dennis Chemnitz

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

发布日期: 2025-12-03

arXiv ID: 2512.04297v1

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

摘要翻译

中文标题: 关于白噪声时间流场在Batchelor尺度下的混合性研究

摘要翻译: 本文研究了在二维环面上具有四个受迫模态的白噪声时间速度场作用下，对流-扩散方程解的混合性质。当扩散系数趋于零时，我们证明了几乎必然的指数耗散率有正的下界。结合Gess与Yaroslavtsev所建立的相应上界估计，这一结果给出了首个可验证Batchelor尺度猜想的速度场实例。此外，我们刻画了该系统在无扩散情形下的指数混合速率。我们的结果并不局限于二维情形，且我们构造了一个具有相同性质的三维白噪声时间速度场。

不相关论文（基本信息）

1. Cumulant expansions of operator groups of quantum many-particle systems

量子多粒子系统算符群的累积量展开

作者: V. I. Gerasimenko, I. V. Gapyak

类别: math-ph, cond-mat.stat-mech, math.AP

发布日期: 2025-12-04

arXiv ID: 2512.05036v1

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

摘要: The article presents a method of cluster expansions for groups of operators associated with the von Neumann equations for states and the Heisenberg equations for observables, aiming to construct generating operators for nonperturbative solutions to the Cauchy problem for hierarchies of evolution equations of many-particle quantum systems.

2. Geophysical intensity problems: the axisymmetric case

地球物理强度问题：轴对称情形

作者: Ralf Kaiser

类别: math.AP, astro-ph.EP

发布日期: 2025-12-04

arXiv ID: 2512.05010v1

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

摘要: Considering the earth or any other celestial body the main sources of the gravitational as well as of the magnetic field lie inside the body. Above the surface both fields are in good approximation harmonic vector fields determined by their values at the body’s surface or any other surface enclosing the body. The intensity problem seeks to determine harmonic vector fields vanishing at infinity and with prescribed intensity of the field at the surface. This problem constitutes a nonlinear boundary value problem, whose general solvability is not yet established. In this paper {\em axisymmetric} harmonic fields ${\bf H}$ outside the unit sphere $S^2$ are studied and, given an axisymmetric Hölder continuous intensity function $I\neq 0$ on $S^2$, the existence of infinitely many solutions of the intensity problem is proved. These solutions can more precisely be characterized as follows: fix a number $\de \in \nat\setminus \{1 \}$ and a meridional plane $M$ through the symmetry axis $S\!A$, and in $M$ a unit circle $S^1$ (symmetric with respect to $S\!A$) and, furthermore, $2\, N$, $N \in \nat_0$, points $z_n \in M$ (symmetric with respect to $S\!A$, avoiding $S\!A$, and outside $S^1$), then the existence of an (up to a sign) unique harmonic field ${\bf H}$ is established that vanishes at (the axisymmetric circles piercing $M$ at) $z_n$ and nowhere else, that has intensity $I$ at $S^2$ and (exact) decay order $\de$ at infinity. The proof is based on the solution of a nonlinear elliptic equation with discontinuous coefficients, which are, moreover, singular at the symmetry axis. Its combination with fixed boundary conditions was the basis of a recent treatment of the ``geomagnetic direction problem’’ \cite{KR22}. Here we have instead natural boundary conditions, which provide less information, and which require, therefore, in part new solution techniques and sharper estimates.

3. Existence and a priori bounds for fully nonlinear PDEs with a harmonic map-like structure

具有类调和映射结构的完全非线性偏微分方程的存在性与先验估计

作者: Gabrielle Nornberg, Ricardo Ziegele

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04961v1

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

摘要: In this paper, we study a new class of fully nonlinear uniformly elliptic equations with a so-called harmonic map-like structure, whose model case is given by \begin{equation*} \mathcal{M}^{\pm}_{λ,Λ}(D^2u) \pm b(x) |Du| \pm β(u)\langle M(x) Du,Du \rangle \pm c(x) u = f(x); \textrm{ in } Ω, \end{equation*} where $Ω\subset \mathbb{R}^n$ is a bounded $C^{1,1}$ domain, $\mathcal{M}^{\pm}$ are the Pucci extremal operators, $β(s) = s^k$ for some $k \in \mathbb{N} $ odd, $b \in L^{q}_{+}(Ω)$, $c,f \in L^p(Ω)$, and $n \leq p \leq q$, $q>n$. We obtain existence results under a smallness regime on the coefficients, along with some classical results such as the Aleksandrov–Bakelman–Pucci estimate and the comparison principle, as well as a priori bounds for the respective Dirichlet problem in the noncoercive case. We also establish multiplicity results and qualitative behavior, which seem to be new in the case of the Laplacian operator.

4. Quantitative rigidity of the Wasserstein contraction under convolution

卷积下Wasserstein收缩的定量刚性

作者: Max Fathi, Michael Goldman, Daniel Tsodyks

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

发布日期: 2025-12-04

arXiv ID: 2512.04928v1

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

摘要: The aim of this paper is to investigate the contraction properties of $p$-Wasserstein distances with respect to convolution in Euclidean spaces both qualitatively and quantitatively. We connect this question to the question of uniform convexity of the Kantorovich functional on which there was substantial recent progress (mostly for $p=2$ and partially for $p>1$). Motivated by this connection we extend these uniform convexity results to the case $p=1$, which is of independent interest.

5. Spectral Theory of Krein-Feller Type Operators and Applications in Stochastic Fractional Elliptic and Parabolic Equations

Krein-Feller型算子的谱理论及其在随机分数阶椭圆与抛物方程中的应用

作者: Kelvin J. R. Almeida-Sousa, Alexandre B. Simas

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

发布日期: 2025-12-04

arXiv ID: 2512.04826v1

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

摘要: It has been shown that the space $C^{\infty}_{W,V}(\mathbb{T})$, introduced in Simas and Sousa (Potential Analysis, 2025), is the natural regularity space for solutions of the eigenvalue problem $Δ_{W,V} u = λu$ on the torus $\mathbb{T}$, where $Δ_{W,V} = \frac{d^{+}}{dV}\frac{d^{-}}{dW}$ is the Krein Feller operator in the case where $W$ and $V$ are strictly increasing and right continuous (respectively left continuous), possibly with dense sets of discontinuities. In this work we provide conditions ensuring that every function in $C^{\infty}_{W,V}(\mathbb{T})$, which may be highly discontinuous, admits a series expansion that generalizes the classical Taylor expansion. A central feature of our approach is that all proofs are nonstandard, since classical analytical and spectral arguments cannot be adapted to this singular setting. Using these methods we characterize the eigenvectors of $Δ_{W,V}$ in terms of generalized trigonometric functions and obtain an asymptotic lower bound for the associated eigenvalues. We also derive a sharp upper bound for the convergence exponent of these eigenvalues, and as a consequence we prove that $C^{\infty}_{W,V}(\mathbb{T})$ is a nuclear space. Further consequences include results on the asymptotic behavior of eigenvalues of compact operators and improvements in traceability. As a final application we establish existence results for generalized fractional stochastic and deterministic differential equations, as well as for parabolic stochastic partial differential equations acting on nuclear spaces.

6. Hierarchical matrix approximability of inverse of convection dominated finite element matrices

对流主导型有限元矩阵逆的分层矩阵可逼近性研究

作者: Arthur Saunier, Leo Agelas, Ani Anciaux Sedrakian, Ibtihel Ben Gharbia, Xavier Claeys

类别: math.NA, math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04824v1

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

摘要: Several researchers have developed a rich toolbox of matrix compression techniques that exploit structure and redundancy in large matrices. Classical methods such as the block low-rank format and the Fast Multipole Method make it possible to manipulate very large systems by representing them in a reduced form. Among the most sophisticated tools in this area are hierarchical matrices (H-matrices), which exploit local properties of the underlying kernel or operator to approximate matrix blocks by low-rank factors, organized in a recursive hierarchy. H-matrices offer a flexible and scalable framework, yielding nearly linear complexity in both storage and computation. Hierarchical matrix techniques, originally developed for boundary integral equations, have recently been applied to matrices stemming from the discretization of advection-dominated problems. However, their effectiveness is limited by the loss of coercivity induced by convection phenomena, where traditional methods fail. Initial work by Le Borne addressed this by modifying the admissibility criterion for structured grids with constant convection, but challenges remain for more general grids and advection fields. In this work, we propose a novel partitioning strategy based on “convection tubes”, clusters aligned with the convection vector field. This method does not require a structured grid or constant convection, overcoming the limitations of previous approaches. We present both theoretical analyses and numerical experiments, that demonstrate the efficiency and robustness of our method for convection-dominated PDEs on unstructured grids. The approach builds on a Péclet-robust Caccioppoli inequality, crucial for handling convection-dominated problems.

7. The initial-to-final-state inverse problem with unbounded potentials and Strichartz estimates

具有无界势的初值-终值反问题及Strichartz估计

作者: Pedro Caro, Alberto Ruiz

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

发布日期: 2025-12-04

arXiv ID: 2512.04796v1

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

摘要: The initial-to-final-state inverse problem consists in determining a quantum Hamiltonian assuming the knowledge of the state of the system at some fixed time, for every initial state. We formulated this problem to establish a theoretical framework that would explain the viability of data-driven prediction in quantum mechanics. In a previous work, we analysed this inverse problem for Hamiltonians of the form $-Δ+ V$ with an electric potential $V = V({\rm t}, {\rm x})$, and we showed that uniqueness holds whenever the potentials are bounded and decay super-exponentially at infinity. In this paper, we extend this result for unbounded potentials. One of the key steps consists in proving a family of suitable Strichartz estimates – including the corresponding endpoint of Keel and Tao. In the context of the inverse Calderón problem this family of inequalities corresponds to the Carleman inequality proved by Kenig, Ruiz and Sogge. Haberman showed that this inequality can be also retrieved as an embedding of a suitable Bourgain space. The corresponding Bourgain space in our context do not capture the mixed-norm Lebesgue spaces of Strichartz inequalities. In this paper, we give a counterexample that justifies this fact, and shows the limitations of Bourgain spaces to address the initial-to-final-state inverse problem.

8. On a fuzzy Landau Equation: Part III. The grazing collision limit

关于模糊化的朗道方程：第三部分. 擦边碰撞极限

作者: Manh Hong Duong, Boris Golubkov, Zihui He

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04713v1

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

摘要: In this paper, we study the grazing limit from the non-cutoff fuzzy Boltzmann equations to the fuzzy Landau equation, where particles interact through delocalised collisions. We show the grazing limit through variational formulations that correspond to the GENERIC (General Equations for Non-Equilibrium Reversible-Irreversible Coupling) structure of the respective equations. We show that the variational formulation associated with a non-quadratic dual dissipation pair for the fuzzy Boltzmann equations converges to a variational formulation of the fuzzy Landau equation corresponding to a quadratic dissipation pair.

9. Infinity of solutions to initial-boundary value problems for linear constant-coefficient evolution PDEs on semi-infinite intervals

关于半无限区间上线性常系数演化偏微分方程初边值问题解的无穷性研究

作者: Andreas Chatziafratis, Spyridon Kamvissis

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

发布日期: 2025-12-04

arXiv ID: 2512.04670v1

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

摘要: In this short communication, we announce an algorithmic procedure for constructing non-uniqueness counter-examples of classical solutions to initial-boundary-value problems for a wide class of linear evolution partial differential equations, of any order and with constant coefficients, formulated in a quarter-plane. Our approach relies on analysis of regularity and asymptotic properties, near the boundary of the spatio-temporal domain, of closed-form integral-representation formulae derived via complex-analytic techniques and rigorous implementation of the modern PDE technique known as Fokas unified transform method. In order to elucidate the novel idea and demonstrate the proposed technique in a self-contained fashion, we explicitly present its application to two concrete examples, namely the heat equation and the linear KdV equation with Dirichlet data. New uniqueness theorems for these two models are also presented herein.

10. Critical concave-convex problems in Carnot group

Carnot群中的临界凹-凸问题

作者: Mattia Galeotti, Eugenio Vecchi

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04640v1

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

摘要: We consider a model Dirichlet problem with concave-convex and critical nonlinearity settled in Carnot groups. Our aim is to prove the existence of two positve solutions in the spirit of a famous result by Ambrosetti, Brezis and Cerami. To this aim we use a variational Perron method combined with proper estimates of a family of functions which are minimizers of the relevant Sobolev inequality. Due to the lack of boundary regularity, we also have to be careful while proving that the first solution found is a local minimizer in the proper topology.

11. Sharp stability on the second Robin eigenvalue with negative boundary parameters

负边界参数下第二Robin特征值的锐利稳定性

作者: Zhijie Chen, Zhen Song, Wenming Zou

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04584v1

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

摘要: In this paper, we prove a quantitative refinement of the isoperimetric type inequality for the second Robin eigenvalue with negative boundary parameters established by Freitas and Laugesen [Amer.J.Math.143 (2021), no.3, 969-994].Such new stability estimate is proved when the boundary parameter is not too far from 0.By constructing a suitable family of nearly spherical domains, we prove that the exponent for the Fraenkel asymmetry in this quantitative type inequality is sharp.

12. Parabolic problems whose Fujita critical exponent is not given by scaling

一类其Fujita临界指数不由尺度变换给出的抛物型问题

作者: Ahmad Z. Fino, Berikbol T. Torebek

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04506v1

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

摘要: This paper investigates the (fractional) heat equation with a nonlocal nonlinearity involving a Riesz potential: \begin{equation*} u_{t}+(-Δ)^{\fracβ{2}} u= I_α(|u|^{p}),\qquad x\in \mathbb{R}^n,,,,t>0, \end{equation*} where $α\in(0,n)$, $β\in(0,2]$, $n\geq1$, $p>1.$ We introduce the Fujita-type critical exponent $p_{\mathrm{Fuj}}(n,β,α)=1+(β+α)/(n-α)$, which characterizes the global behavior of solutions: global existence for small initial data when $p>p_{\mathrm{Fuj}}(n,β,α),$ and finite-time blow-up when $p\leq p_{\mathrm{Fuj}}(n,β,α)$. It is remarkable that the critical Fujita exponent is not determined by the usual scaling argument that yields $p_{sc}=1+(β+α)/n$, but instead arises in an unconventional manner, similar to the results of Cazenave et al. [Nonlinear Analysis, 68 (2008), 862-874] for the heat equation with a nonlocal nonlinearity of the form $\int_0^t(t-s)^{-γ}|u(s)|^{p-1}u(s)ds,\,0\leq γ<1.$ The result on global existence for $p>p_{\mathrm{Fuj}}(n,2,α),$ provides a positive answer to the hypothesis proposed by Mitidieri and Pohozaev in [Proc. Steklov Inst. Math., 248 (2005) 164-185]. We further establish global nonexistence results for the above heat equation, where the Riesz potential term $I_α(|u|^{p})$ is replaced by a more general convolution operator $(\mathcal{K}\ast |u|^p),\,\mathcal{K}\in L^1_{loc}$, thereby extending the Mitidieri-Pohozaev’s results established in the aforementioned work. Proofs of the blow-up results are obtained using a nonlinear capacity method specifically adapted to the structure of the problem, while global existence is established via a fixed-point argument combined with the Hardy-Littlewood-Sobolev inequality.

13. Irreversibility condition and stability of equilibria in the inverse-deformation approach to fracture

反变形方法中裂纹问题的不可逆性条件与平衡稳定性

作者: Arnav Gupta

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

发布日期: 2025-12-04

arXiv ID: 2512.04479v1

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

摘要: We derive the irreversibility condition in fracture for the inverse-deformation approach using the second law of thermodynamics. We consider the problem of brittle failure in an elastic bar previously solved in (Rosakis et al 2021). Despite the presence of a non-zero interfacial/surface energy, the third derivative of the inverse-deformation map is discontinuous at the crack faces. This is due to the presence of the inequality constraint ensuring the inverse strain is nonnegative and the orientation of matter is preserved. A change in the material location of a crack results in negative entropy production, violating the second law. Consequently, such changes are disallowed giving the irreversibility condition. The inequality constraint and the irreversibility condition limit the space of admissible variations. We prove necessary and sufficient conditions for local stability that incorporate these restrictions. Their numerical implementation shows that all broken equilibria found in (Rosakis et al 2021) are locally stable.

14. A note on lifespan estimates for higher-order parabolic equations

关于高阶抛物型方程解的寿命估计的一点注记

作者: Nurdaulet N. Tobakhanov, Berikbol T. Torebek

类别: math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04428v1

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

摘要: We investigate the lifespan of solutions to the higher-order semilinear parabolic equation

$$ u_t+(-Δ)^m u=|u|^p, \quad x \in \mathbb{R}^n, t>0 $$

with initial data. We focus on the precise asymptotic behavior of the lifespan of nontrivial solutions. By combining the test function method and semigroup estimates, we derive both upper and lower bounds for the lifespan of solutions

$$ T_{\varepsilon} \simeq \left\{\begin{array}{l}\varepsilon^{-\left(\frac{1}{p-1}-\frac{n}{2m}\right)^{-1}}, \,\, 1<p<p_{\text {Fuj}}, \\ \exp\left(\varepsilon^{-(p-1)}\right), \,\, p=p_{\text {Fuj}},\end{array}\right. $$

where $p_{Fuj}=1+\frac{2m}{n}$ is the critical exponent of Fujita. These estimates refine and extend the earlier results of Caristi-Mitidieri [J. Math. Anal. Appl., 279:2 (2003), 710-722] and Sun [Electron. J. Differential Equations, 17 (2010)], who obtained only upper bounds under slowly decaying initial data assumptions. In our setting, the above condition on the initial data is replaced by the assumption $L^1\cap L^\infty$, which sharpens the results of the aforementioned works.

15. The Dirichlet heat trace for domains with curved corners

具有曲角区域的狄利克雷热迹

作者: Shi Zhuo Looi, David Sher

类别: math.SP, math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04422v1

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

摘要: We study the short-time asymptotics of the Dirichlet heat trace on planar curvilinear polygons. For such domains we show that the coefficient of $t^{1/2}$ in the expansion splits into a boundary integral of $κ^2$ and a sum of local corner contributions, one for each vertex. Each curved corner contribution depends only on the interior angle $α$ and on the limiting curvatures $κ_{\pm}$ on the adjacent sides. Using a conformal model and a parametrix construction on the sector heat space, we express this contribution in the form $c_{1/2}(α)\,r_0(α,κ_+, κ_-)$, where $c_{1/2}(α)$ is given by a Hadamard finite part of an explicit trace over the exact sector. For right-angled corners we compute $c_{1/2}(π/2)=1/(16\sqrtπ)$ and obtain a closed formula for the $t^{1/2}$ coefficient. As an application we extend a previous result in the literature by showing that any admissible curvilinear polygon that is Dirichlet isospectral to a polygon must itself be a polygon with straight sides.

16. Homogenizationof non-divergence form operators in i.i.d. random environments

独立同分布随机环境中非散度型算子的均匀化

作者: Xiaoqin Guo, Timo Sprekeler, Hung V. Tran

类别: math.PR, math.AP

发布日期: 2025-12-04

arXiv ID: 2512.04410v1

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

摘要: We study random walks in a balanced, i.i.d. random environment in $\mathbb Z^d$ for $d\geq 3$. We establish improved convergence rates for the homogenization of the Dirichlet problem associated with the corresponding non-divergence form difference operators, surpassing the $O(R^{-1})$ rate, which is expected to be optimal for environments with a finite range of dependence. In particular, the improved rates are $O(R^{-3/2})$ when $d=3$, and $O(R^{-2}\log R)$ when $d\geq 4$.

17. Regularity for minimizers of degenerate, non-autonomous, orthotropic integral functionals

退化、非自治、正交各向异性积分泛函极小化子的正则性

作者: Antonio Giuseppe Grimaldi, Stefania Russo

类别: math.AP

发布日期: 2025-12-03

arXiv ID: 2512.04281v1

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

摘要: We prove the higher differentiability of integer order of locally bounded minimizers of integral functionals of the form \begin{equation*} \mathcal{F}(u,Ω):= ,\sum_{i=1}^{n} \dfrac{1}{p_i}\displaystyle \int_Ω, a_i(x) \lvert u_{x_i} \rvert^{p_i} dx- \int_Ωω(x)u(x) dx, \end{equation*} where the exponents $ p_i \geq 2 $ and the coefficients $ a_i(x) $ satisfy a suitable Sobolev regularity. The main novelty consists in dealing with non-autonomous, anisotropic functionals, which depend also on the solution.

18. Phase mixing and the Vlasov equation in cosmology

相位混合与宇宙学中的Vlasov方程

作者: Martin Taylor, Renato Velozo Ruiz

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

发布日期: 2025-12-03

arXiv ID: 2512.04214v1

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

摘要: We consider the Vlasov equation on slowly expanding isotropic homogeneous tori, described by the Friedmann–Lemaître–Robertson–Walker cosmological spacetimes. For expansion rate $t^q$, with $0< q<\frac{1}{2}$ (excluding certain exceptional values), we show that the spatial density decays at the rate $t^{-6q}$ and that, when the spatial average is removed, the density decays at an enhanced rate due to a phase mixing effect. This enhancement is polynomial for Sobolev initial data and super-polynomial, but sub-exponential, for real analytic initial data. We further show that, when the expansion rate is the borderline $t^{\frac{1}{2}}$ – the rate which describes a radiation filled universe – a degenerate phase mixing effect results in a logarithmic enhancement for Sobolev initial data and a super-logarithmic enhancement (in fact, a gain of $\exp(-μ(\log t)^ε)$ for some $μ,ε>0$) for analytic initial data. The proof is based on a collection of commuting vector fields, and certain combinatorial properties of an associated collection of differential operators. The vector fields are not explicit, but are shown to have good properties when $t$ is large with respect to the momentum support of the solution. A physical space dyadic localisation is employed to treat non-compactly supported (in particular, non-trivial real analytic) but suitably decaying solutions.