生成时间: 2026年01月23日 03:59:21

重点关注论文数量: 5

了解领域论文数量: 3

不相关论文数量: 12

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

1. 二维各向异性Navier-Stokes方程在具有水平分数耗散下的稳定性与衰减

Stability and Decay for the 2D Anisotropic Navier-Stokes Equations with Fractional Horizontal Dissipation on $\mathbb{R}^2$

作者: Zhibin Wang, Jiahong Wu, Ning Zhu

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.16110v1

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

详细分析

中文标题: 二维各向异性Navier-Stokes方程在具有水平分数耗散下的稳定性与衰减

研究对象和背景

这篇论文研究的是二维不可压缩Navier-Stokes方程的一个特殊变体，该方程描述了在$ \mathbb{R}^2 $（整个二维欧几里得平面）上流动的牛顿流体。与标准的各向同性Navier-Stokes方程不同，其耗散项仅作用在一个方向上，并且是分数阶的。具体来说，研究的方程组如下：

$$ \begin{cases} \partial_t u + \nu \Lambda_1^{2s} u + u \cdot \nabla u = -\nabla p, \quad (x_1, x_2) \in \mathbb{R}^2, \ t > 0, \\ \nabla \cdot u = 0, \\ u|_{t=0} = u_0. \end{cases} $$

其中，$ u = (u_1, u_2) $ 是速度场，$ p $ 是压力，$ \nu > 0 $ 是粘性系数。关键的算子是 $ \Lambda_1^{2s} $，它定义为方向 $ x_1 $ 上的分数阶导数（水平方向）。通过傅里叶变换，其作用方式为 $ \mathcal{F}(\Lambda_1^{2s} f)(\xi) = |\xi_1|^{2s} \hat{f}(\xi) $，其中 $ \xi = (\xi_1, \xi_2) $ 是对偶变量，$ 0 \leq s < 1 $。当 $ s = 1 $ 时，算子退化为标准的拉普拉斯算子 $ \partial_{x_1}^2 $。

这篇论文研究的核心问题是：在 $ s \in [0, 1) $ 的范围内，该系统是否具有稳定性？解随时间的衰减行为如何？

背景与物理意义：

  1. 中间状态：该系统处于完全耗散（$ s=1 $，标准Navier-Stokes）和无耗散（$ \nu=0 $，欧拉方程）之间的中间状态。标准的Navier-Stokes方程在二维情况下，解随时间呈现代数衰减，系统是渐近稳定的。相反，二维欧拉方程虽然能保证整体存在性和唯一性（根据Yudovich理论），但解的高阶Sobolev范数（如涡度梯度）可能会随时间快速增长（例如双指数增长），表现出不稳定性。这篇论文探讨在单方向弱耗散下，系统的长期行为更偏向于Navier-Stokes（稳定、衰减）还是欧拉方程（不稳定、增长）。
  2. 各向异性与实际应用：这类模型源于地球物理流体动力学、薄层流体和边界层流。在这些情境下，水平方向和垂直方向的尺度差异巨大，耗散机制往往只在一个方向上占主导地位。此外，单方向的分数阶耗散可以作为非局部混合或异常扩散的有效描述。
  3. 数学难点：
    • 非线性项的控制：由于耗散只在一个方向起作用，传统能量估计无法直接控制所有非线性项的高阶导数。例如，在估计涡度梯度的$L^2$范数时，会出现无法积分控制的项。
    • 全空间的困难：之前的研究多集中在部分周期区域（如 $ \mathbb{T} \times \mathbb{R} $），可以利用Poincaré不等式分离零模。但在全空间 $ \mathbb{R}^2 $ 上，缺乏这种结构，使得稳定性分析更具挑战性。

区域与函数空间：

  • 区域：全空间 $ \mathbb{R}^2 $。
  • 边界条件：无（由于是全空间，通常要求解在无穷远处有适当的衰减）。
  • 函数空间：主要使用Sobolev空间 $ H^k(\mathbb{R}^2) $，以及各向异性空间（结合水平和垂直方向的不同导数性质）。为了处理特定的衰减率，还引入了带有截断多项式权重的加权 $ L^2 $ 空间。

主要定理或主要结果

论文根据分数阶耗散指数 $ s $ 的不同范围，给出了三个主要的定理。

1. 定理 1.1 (针对 $ 0 \leq s \leq \frac{3}{4} $)：

  • 假设条件：$ s $ 在此范围内，初始数据 $ u_0 \in H^k(\mathbb{R}^2) $ 且 $ k \geq 3 $，散度为零（$ \nabla \cdot u_0 = 0 $）。初始数据足够小（$ \|u_0\|_{H^k} \leq \epsilon $）。
  • 结论：系统存在唯一的全局解 $ u \in C([0, \infty); H^k(\mathbb{R}^2)) $，并且满足均匀估计：
    $$ \|u(t)\|_{H^k}^2 + \nu \int_0^t \| \Lambda_1^{s+1} u(\tau) \|_{H^k}^2 d\tau \leq C\epsilon, \quad \forall t \geq 0. $$
    这表明在此范围内，即使耗散较弱（$ s $ 小），非线性项也能通过能量方法得到控制。

2. 定理 1.2 (针对 $ \frac{3}{4} < s < \frac{11}{12} $)：

  • 假设条件：$ \frac{3}{4} < s < \frac{11}{12} $，且引入参数 $ \sigma $ 满足 $ \frac{3}{4} < s < \frac{5}{12} + \sigma $ 和 $ \frac{1}{3} < \sigma < \frac{1}{2} $。初始数据 $ u_0 \in H^k(\mathbb{R}^2) $，散度为零，且存在流函数 $ \psi_0 $ 使得 $ u_0 = \nabla^\perp \psi_0 $。
  • 关键假设：初始数据具有负阶水平Sobolev正则性，即 $ \Lambda_1^{-\sigma} u_0, \Lambda_1^{-\sigma} \psi_0, \Lambda_1^{-\sigma} \partial_{x_1} u_0 \in L^2 $，且范数足够小（见公式 1.4）。
  • 结论：存在唯一的全局解，且满足：
    $$ \|u(t)\|_{L^2} \leq C \epsilon (1+t)^{-\sigma/2s}. $$
    特别地，第二分量 $ u_2 $ 比第一分量 $ u_1 $ 衰减得更快。

3. 定理 1.5 (针对 $ \frac{3}{4} < s < 1 $，特别是 $ \frac{11}{12} \leq s < 1 $)：

  • 假设条件：$ s \in (\frac{3}{4}, 1) $。引入截断幂权重函数 $ r_{x_2}^\gamma $（当 $ |x_2| \leq 1 $ 时为1，当 $ |x_2| > 1 $ 时为 $ |x_2|^\gamma $）。假设初始数据 $ u_0 $ 及其导数在 $ x_2 $ 方向上具有适当的加权衰减（见公式 1.6）。
  • 结论：在加权空间中，系统存在唯一的全局解。并且解具有加权衰减估计，例如：
    $$ \| r_{x_2}^{\frac{3\gamma+4}{2s}} u(t) \|_{L^2} \leq C_0 \epsilon (1+t)^{-\sigma/2s}. $$
    同时，对于 $ u_2 $ 分量，有更强的加权衰减：
    $$ \| r_{x_2}^{s\gamma} u_2(t) \|_{L^2} \leq C_0 \epsilon (1+t)^{-(\sigma+1)/2s}. $$
    此结果成功地将稳定性分析扩展到了 $ s \in [0, 1) $ 的全范围。

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

这篇论文的证明体现了现代偏微分方程分析中处理各向异性非线性问题的精妙技巧。作者针对不同的 $ s $ 范围采用了截然不同的策略，特别是引入了加权能量估计来攻克 $ s $ 接近 1 时的难点。

a) 核心分析学工具：

  1. 各向异性能量估计 (Anisotropic Energy Estimates)： 这是处理该方程的基础。由于耗散仅在 $ x_1 $ 方向（记为水平方向），能量估计通常沿着 $ \partial_{x_1} $ 进行。 在定理 1.1 ($ 0 \leq s \leq \frac{3}{4} $) 中，作者直接利用高频部分的时间积分可控性。通过在 $ H^k $ 框架下做能量估计，利用 $ s $ 较小使得需要控制的导数阶数较低（由 Leibniz 法则），从而使得某些非线性项的时间积分在 $ s \leq 3/4 $ 范围内可被吸收。

  2. Sobolev 嵌入与各向异性乘积法则 (Anisotropic Product Estimates)： 为了控制非线性项 $ u \cdot \nabla u $，论文引用了关键的引理（如 Lemma 2.1 和 2.2）：

    • $ \|fg\|_{L^1_{x_1} L^2_{x_2}} \leq C \|f\|_{L^2}^{1/2} \|\partial_{x_1} f\|_{L^2}^{1/2} \|g\|_{L^2} $
    • $ \|fg\|_{L^2} \leq C \|f\|_{L^2}^{1/2} \|\partial_{x_1} f\|_{L^2}^{1/2} \|g\|_{L^2}^{1/2} \|\partial_{x_1} g\|_{L^2}^{1/2} $ 这些估计利用了 Gagliardo-Nirenberg 不等式在不同方向上的混合范数性质，是处理非线性项的核心。

  3. 傅里叶分裂法 (Fourier Splitting Method) 与 衰减估计： 在定理 1.2 中，为了处理 $ s > 3/4 $ 的情况，作者结合了衰减分析。这种方法通过将频率空间划分为高频和低频部分，利用分数阶热核的衰减性质来获得解的 $ L^2 $ 衰减率。这需要初始数据具有负阶正则性（$ \Lambda_1^{-\sigma} \in L^2 $），以保证低频部分的衰减足够快。

  4. 加权 $ L^2 $ 空间与 $ A_2 $ 权重理论 (Weighted $ L^2 $ Spaces and $ A_2 $ Weight Theory)： 这是解决 $ s \geq \frac{11}{12} $ 困难区域的“杀手锏”。

    • Riesz 变换的有界性：在压力项的估计中，压力由速度场通过 Biot-Savart 定律给出，涉及到 Riesz 变换。在标准 $ L^2 $ 空间中 Riesz 变换是有界的，但在各向异性耗散下，直接的能量估计会丢失信息。引入权重 $ r_{x_2}^\gamma $ 后，必须证明 Riesz 变换在加权 $ L^2(r_{x_2}^\gamma dx) $ 空间上的有界性。这直接依赖于 Muckenhoupt $ A_2 $ 权重理论。
    • 截断多项式权重：定义 $ r_{x_2}^\gamma(x) = 1 $ 当 $ |x_2| \leq 1 $，否则 $ |x_2|^\gamma $。这种权重在单位带内为常数，在外部增长，既保证了局部估计的平滑性，又利用了外部的大权重来吸收非线性项。

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

  1. 针对 $ 0 \leq s \leq \frac{3}{4} $ 的导数重分布： 在这个范围，虽然 $ s $ 小意味着耗散弱，但也意味着在能量估计中，为了控制时间积分，需要移动到非线性项上的导数阶数较低。作者通过精细的积分估计，发现当 $ s \leq 3/4 $ 时，非线性产生的项可以被耗散项 $ \|\Lambda_1^{s+1} u\|^2 $ 的时间积分所控制，从而不需要借助衰减估计即可关闭能量不等式。

  2. 针对 $ \frac{3}{4} < s < \frac{11}{12} $ 的混合稳定性-衰减分析： 当 $ s $ 超过 $ 3/4 $，纯能量估计失效。此时，必须利用解的衰减性来辅助控制非线性项。

    • 流函数表示：利用 $ u = \nabla^\perp \psi $，将速度场转化为流函数，有助于利用涡度方程的结构。
    • 负阶正则性传播：假设 $ \Lambda_1^{-\sigma} u_0 \in L^2 $，并证明该性质在演化中保持。这是获得代数衰减率的关键。
    • 各向异性衰减：证明中发现 $ u_2 $ 分量比 $ u_1 $ 衰减得更快。这是因为压力梯度项主要平衡 $ u_1 $ 方向的对流，而 $ u_2 $ 分量受此影响较小，从而获得了更好的衰减率。这种不对称衰减对于控制高阶非线性项至关重要。

  3. 针对 $ \frac{11}{12} \leq s < 1 $ 的加权能量估计： 这是最难的部分，因为现有的衰减估计（定理 1.2）无法提供足够快的衰减来覆盖 $ s $ 接近 1 的情形。

    • 引入权重克服 Riesz 变换难题：在标准 $ L^2 $ 空间中，压力项 $ p $ 的估计涉及 Riesz 变换 $ R_j $。在全空间且仅水平耗散的情况下，$ R_j $ 与 $ \Lambda_1^{s} $ 的交换子估计非常困难，导致能量估计无法闭合。通过引入垂直方向 $ x_2 $ 的多项式权重 $ r_{x_2}^\gamma $，作者利用 $ A_2 $ 权重理论证明了 Riesz 变换在加权空间中的有界性。这使得压力项可以被有效控制。
    • 权重的选择：权重的指数（如 $ \frac{3\gamma+4}{2s} $）是经过精心计算的，旨在抵消非线性项在演化过程中产生的“损失”。
    • 截断技巧：权重在 $ |x_2| \leq 1 $ 处截断为 1，这保证了在局部区域（物理上最重要的区域）的正则性不受权重奇异性的影响，同时在远处利用权重增长来控制大范围的非线性相互作用。

  4. 不动点方法 (Fixed Point Argument)： 在所有定理的证明中，最终都归结为在一个适当的函数空间中构造映射，并证明其是压缩映射或满足先验估计。这需要精心选取的函数空间（如带有时间权重或空间权重的空间），并建立相应的先验估计。

与之前工作的比较

这篇论文在多个方面对现有文献做出了显著的改进，特别是在处理全空间 $ \mathbb{R}^2 $ 上的各向异性耗散问题上取得了突破。

a) 突破与改进：

  1. 从部分周期域到全空间 $ \mathbb{R}^2 $ 的推广：

    • 先前工作：Dong, Wu, Xu, 和 Zhu (2019) 在 $ \mathbb{T} \times \mathbb{R} $（水平周期，垂直全空间）上证明了仅水平耗散 Navier-Stokes 方程的稳定性。
    • 局限性：该方法依赖于分离零水平傅里叶模式（$ \xi_1 = 0 $）并利用 Poincaré 不等式控制垂直方向的导数。在全空间 $ \mathbb{R}^2 $ 上，傅里叶变换是连续的，无法简单分离零模，且缺乏 Poincaré 不等式，导致直接平移该方法失败。
    • 本文改进：本文在 $ \mathbb{R}^2 $ 上解决了这一问题。对于 $ s \leq 3/4 $，通过精细的各向异性能量估计直接处理；对于更大的 $ s $，通过引入加权空间和衰减估计，巧妙地绕过了对零模分离的依赖。

  2. 耗散范围的全面覆盖 ($ s \in [0, 1) $)：

    • 先前结果：对于 $ s $ 接近 1 的情况（强各向异性，接近标准水平拉普拉斯），标准的能量方法和衰减方法往往失效，因为非线性项的增长速度超过了耗散能提供的衰减速度。
    • 本文改进：论文创新性地引入了 $ A_2 $ 加权空间理论。这在Navier-Stokes方程的研究中并不常见，通常用于处理椭圆或抛物型方程的边界正则性。本文将其用于解决全空间非线性项的控制问题，特别是解决了压力项中 Riesz 变换在加权空间的有界性问题，从而将存在性和稳定性结果扩展到了 $ 0 \leq s < 1 $ 的整个区间。

  3. 解的衰减率与各向异性结构的揭示：

    • 先前认识：对于各向异性耗散，解的衰减行为通常难以精确刻画，且容易受到各方向衰减不匹配的影响。
    • 本文改进：论文不仅给出了衰减估计，还揭示了显著的各向异性衰减现象。在定理 1.2 和 1.5 中，明确指出 $ u_2 $ 分量比 $ u_1 $ 分量衰减得更快（例如，$ \|u_2\| \lesssim (1+t)^{-(\sigma+1)/2s} $，而 $ \|u_1\| \lesssim (1+t)^{-\sigma/2s} $）。这种精细的结构分析不仅有助于关闭非线性估计，也加深了对物理系统动力学的理解。

  4. 对初始数据正则性的灵活要求：

    • 在高 $ s $ 区域（定理 1.2 和 1.5），论文允许初始数据在水平方向具有负阶 Sobolev 正则性（即 $ \Lambda_1^{-\sigma} u_0 \in L^2 $）。这放宽了对初始数据光滑性的要求，使得结果适用于更广泛的初值类别（如某些 $ L^2 $ 函数），只要它们在适当的尺度下具有一定的可积性。

总结： 该论文通过分阶段、多策略的方法，成功攻克了二维各向异性 Navier-Stokes 方程在全空间上的稳定性难题。特别是利用加权 $ L^2 $ 空间和 $ A_2 $ 权重理论处理 $ s $ 接近 1 的困难情形，是方法论上的一个重要创新。这不仅解决了该模型的长期存在性问题，还详细刻画了其代数衰减的性质，为理解弱耗散流体系统的稳定性提供了重要的理论工具。

2. 球体上的无粘流体动力学与气动声学实时计算

Real-Time Inviscid Fluid Dynamics and Aero-acoustics on a Sphere

作者: Madhusraba Sinha, Jan Stratmann

类别: math.AP, physics.flu-dyn

发布日期: 2026-01-22

arXiv ID: 2601.15982v1

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

详细分析

中文标题: 球体上的无粘流体动力学与气动声学实时计算

1. 研究对象和背景

本研究的核心对象是球面上带有嵌入式障碍物的无粘性（inviscid）不可压缩流体动力学及其伴随的气动声学（aero-acoustics） 生成过程。目标是为复杂曲面（以球体为例）上的流体与声学现象提供一套实时（real-time） 的计算框架。这项工作在科学可视化、虚拟现实（VR）和教育工具等领域具有潜在的交互式应用价值。

传统的基于网格的求解器在处理表面奇异性（如球坐标中的极点）时往往面临数值不稳定性问题，而基于网格的方法在求解具有稳定高阶精度的偏微分方程（PDEs）时缺乏直接的路径。本文提出了一种统一的框架，结合了最近点方法（Closest Point Method, CPM）、基于投影的 Navier-Stokes 求解器以及 Ffowcs Williams-Hawkings (FWH) 类比，旨在克服这些挑战。

1.1 数学模型与方程系统

文章研究的流体动力学主要基于欧拉方程（Euler equations），描述了无粘性不可压缩流体的运动。尽管欧拉方程是双曲型方程组，但在低马赫数近似下，其守恒形式常用于描述流体的惯性与压力驱动的运动。PDF内容第6页给出了具体的控制方程：

1.1.1 流体动力学方程

对于不可压缩、无粘性流体，其动量方程和连续性方程如下：

动量方程：

$$ \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} = -\frac{1}{\rho_0} \nabla p + \mathbf{f} $$

连续性方程（不可压缩条件）：

$$ \nabla \cdot \mathbf{u} = 0 $$

其中：

  • $\mathbf{u}$ 是流体速度场。
  • $p$ 是压力场。
  • $\rho_0$ 是流体密度（假设为常数）。
  • $\mathbf{f}$ 是局部外力项，用于在球体顶部附近启动并维持流体运动。

在该模型中，研究者采用了无粘性近似（$\nu = 0$），重点关注由惯性和压力驱动的流体动力学，因为这在低马赫数下是声生成的主要机制。

1.1.2 PDE 定义域与特殊性

本研究并未直接在欧几里得空间中求解上述方程，而是在**曲面（球体）**上求解。这意味着上述 PDE 的定义域 $\Omega$ 不是通常的三维体区域，而是嵌入在 $\mathbb{R}^3$ 中的二维流形 $\mathcal{S}$（球面）。由于球面具有非平凡的几何结构，直接在球坐标系下离散化会遇到极点奇异性（singularity）和网格间距不均匀的问题。

因此，本文采用了嵌入式方法（embedding method）。PDE 的求解区域实际上是一个包围球体的窄带（narrow band） 区域 $\Omega$，位于 $\mathbb{R}^3$ 中。该方法的核心在于将曲面上的解延伸到体积域中，利用标准的笛卡尔网格进行计算，从而避免了曲面坐标映射的复杂性。

1.1.3 边值条件与障碍物

系统中包含嵌入的表面障碍物。这些障碍物通过符号距离函数（Signed Distance Functions, SDFs） 进行隐式建模。

  • 无滑移条件（No-slip）：在障碍物内部，速度被强制设为零。
  • 无穿透条件（No-penetration）：在障碍物表面，速度被投影到切平面上。
  • 伯努利压力调整：为了保持与真实世界边界相互作用的一致性，模型结合了基于伯努利原理的压力调整。

1.1.4 气动声学模型

气动声学部分基于 Ffowcs Williams-Hawkings (FW-H) 类比。该类比将流体中的非定常波动转化为声辐射方程。在低马赫数和紧凑源（compact source）假设下，远场压力波动 $p'(x,t)$ 由下式给出：

$$ p'(x,t) \approx \frac{1}{4\pi c_0} \frac{i}{r} \cdot \frac{d\mathbf{F}(t)}{dt} $$

其中：

  • $c_0$ 是环境声速。
  • $r$ 是观测者距离。
  • $\mathbf{i}$ 是辐射方向单位向量。
  • $\mathbf{F}(t) = \int_S p(y,t)\mathbf{n}(y)dS$ 是作用在表面 $S$ 上的非定常气动力。
  • 声源主要来自表面压力梯度的时变导数。

1.2 现象描述

该模型旨在模拟和验证以下物理现象：

  1. 球面上的无粘性流动：观察流体在球体表面绕过嵌入障碍物的流动模式。
  2. 压力与声的耦合：验证流体表面的压力波动（通过 FW-H 类比）如何实时转换为声学信号。
  3. 数值稳定性与几何一致性：验证在球面几何结构下，使用 CPM 方法结合投影求解器是否能保持流体运动的数值稳定性和几何约束的满足。

2. 主要定理或主要结果

本文作为一篇计算数学与应用物理的论文，其主要“结果”更侧重于数值方法的性能验证、收敛性分析以及计算效率的展示，而非传统分析 PDE 中的存在性或正则性定理。基于论文内容，主要结果可归纳如下：

2.1 数值收敛性与稳定性 (Manufactured Solutions)

为了验证求解器的正确性，论文使用了制造解方法（Method of Manufactured Solutions, MMS）。这是验证数值算法收敛阶的主要手段。

  • 核心结果：通过构造满足欧拉方程和不可压缩条件的解析解（速度场 $\mathbf{u}$ 和压力场 $p$），验证数值解的误差随网格细化的收敛速率。
  • 预期结论：在理想情况下，对于二阶精度的有限差分格式，数值误差 $E$ 应与网格步长 $h$ 满足 $E = O(h^2)$。
  • 验证指标：论文（第10页及之后，虽未完全展示，但根据 MMS 章节标题推断）展示了 L2 范数误差的收敛曲线，证明了 CPM 投影求解器在球面几何下的正确性。

2.2 实时计算能力与效率

  • 结果：所提出的框架能够实现实时（Real-Time） 或近实时的流体与声学模拟。
  • 技术支撑：通过使用稀疏网格（Sparse Grids） 技术和针对小规模系统的 Cholesky 分解预计算，以及针对大规模系统的共轭梯度法（Conjugate Gradient），在保持精度的同时显著降低了计算成本，使其适用于交互式应用（如 VR）。

2.3 音频生成的物理一致性

  • 结果：通过频谱分析提取的气动力主导频率 $\omega_k$ 与合成的音频信号频率一致，且振幅 $A(t)$ 通过非线性映射 $A(t) = \tanh(\alpha|p'(t)|)$ 受控于瞬时声压。
  • 结论：该方法不仅在视觉上模拟了流体，还在听觉上提供了物理上一致的声学反馈，证明了流体-声学耦合模型的有效性。

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

本研究融合了数值分析、流体力学计算和信号处理技术。以下是核心方法的详细描述。

3.1 核心分析学与数值工具

3.1.1 最近点方法 (Closest Point Method, CPM)

这是处理曲面 PDE 的核心数学工具。

  • 原理：CPM 将定义在流形 $\mathcal{S}$ 上的 PDE 嵌入到包围该流形的欧几里得空间 $\mathbb{R}^3$ 的一个窄带 $\Omega$ 中。
  • 映射算子：对于带内的任意点 $x \in \Omega$，其上的场值定义为曲面上最近点 $CP(x)$ 处的值：
    $$ u(x) = u(CP(x)), \quad CP(x) = R \frac{x}{\|x\|} $$
    其中 $R$ 是球半径。
  • 优势：这允许直接在嵌入空间应用标准的笛卡尔有限差分算子（如拉普拉斯算子），而无需在曲面上生成复杂的参数化网格，从而避免了极点奇异性。

3.1.2 投影方法 (Projection Method / Chorin’s Method)

用于求解不可压缩欧拉方程，将速度场更新分解为两个步骤，这是一种经典的分裂格式：

  1. 预测步 (Advection)：忽略压力项，计算中间速度 $\mathbf{u}^*$。
    $$ \frac{\mathbf{u}^* - \mathbf{u}^n}{\Delta t} = -(\mathbf{u}^n \cdot \nabla) \mathbf{u}^n + \mathbf{f} $$
    这里采用了 半拉格朗日格式 (Semi-Lagrangian scheme) 进行对流项离散。半拉格朗日法通过沿特征线回溯来寻找当前速度场，具有无条件稳定性，非常适合大时间步长的实时模拟。
  2. 校正步 (Projection)：求解压力泊松方程以强制不可压缩条件 $\nabla \cdot \mathbf{u}^{n+1} = 0$。
    $$ \nabla^2 p^{n+1} = \frac{\rho_0}{\Delta t} \nabla \cdot \mathbf{u}^* $$
    得到压力后，更新速度：
    $$ \mathbf{u}^{n+1} = \mathbf{u}^* - \frac{\Delta t}{\rho_0} \nabla p^{n+1} $$

3.1.3 线性代数求解器

压力泊松方程的求解是计算瓶颈。文中采用了混合策略：

  • Cholesky 分解：对于小网格系统，预先计算矩阵的 Cholesky 分解。这是一种直接法，当系统规模较小时速度极快。
  • 共轭梯度法 (Conjugate Gradient, CG)：对于较大的系统，使用迭代法 CG 求解，以平衡内存占用和计算时间。

3.1.4 频谱分析与合成 (Fourier Analysis & Synthesis)

在声学部分，利用离散傅里叶变换 (DFT) 分析气动力的时间导数 $d\mathbf{F}(dt)$。

  • 工具：滑动时间窗 + FFT。
  • 应用：提取主导频率 $\omega_k$ 和幅度 $|\hat{F}(\omega_k)|$，用于驱动合成器的正弦波振荡器，实现物理基础的声音增强。

3.2 解决技术性难题的技巧

3.2.1 处理曲面奇异性与参数化

  • 难题：传统的球坐标网格在两极（Poles）存在度量奇异性，导致网格点聚集，计算成本高且易产生数值误差。
  • 技巧：采用 CPM + 笛卡尔窄带。不再直接在球面上定义网格，而是在 $\mathbb{R}^3$ 的均匀立方网格上定义计算域。通过限制计算在表面附近的窄带内（通常几个网格宽度），既保留了笛卡尔网格的规则性（无奇异性），又大幅减少了计算量。
  • 投影：将嵌入空间的算子（如拉普拉斯算子 $\Delta_{CPM}$）投影回曲面的切空间，确保几何一致性。

3.2.2 隐式障碍物处理 (Signed Distance Functions)

  • 难题：在传统网格方法中，处理复杂形状的障碍物需要贴体网格或繁琐的切割网格（Cut-cell）技术，这在实时应用中难以实现。
  • 技巧：使用 SDF 进行几何建模。
    • 定义函数 $\phi(\mathbf{x})$ 表示点 $\mathbf{x}$ 到障碍物表面的有符号距离。
    • 内部约束：若 $\phi(\mathbf{x}) < 0$，强制 $\mathbf{u} = 0$（或在切平面上的投影）。
    • 这种隐式表示使得添加或移动障碍物无需重新生成网格，只需更新 SDF 场即可，极大地提升了灵活性和实时性。

3.2.3 压力泊松方程的边界条件处理

  • 难题：在 CPM 的窄带方法中，如何正确施加 Neumann 边界条件（如压力梯度在边界上的条件）是一个挑战。
  • 技巧：通过构建带有边界约束的稀疏拉普拉斯矩阵来实现。在窄带边界处，算子的构建被修正以反映 Neumann 条件（即 $\nabla p \cdot \mathbf{n} = 0$ 或根据障碍物调整）。通过矩阵的预处理（如 Cholesky），保证了边界条件的精确满足。

3.2.4 音频生成中的伪影抑制 (Artifact Suppression)

  • 难题：直接从流体模拟数据生成音频可能会产生断续的噪声或突变。
  • 技巧：
    • 滞后平滑 (Hysteresis Smoothing)：在频域提取中引入时间窗和滞后机制，防止频率跳变。
    • 非线性映射：使用双曲正切函数 $A(t) = \tanh(\alpha|p'(t)|)$ 将声压映射到振幅。这不仅防止了数字信号的削波（Clipping），还模拟了人耳对响度的非线性感知（响度饱和）。

4. 与之前工作的比较

本文旨在弥合高保真流体模拟与实时交互应用之间的差距。以下是与现有技术的详细对比：

4.1 突破与改进概述

论文声称的主要突破在于构建了一个统一的、非参数化的框架，该框架能够在带有几何障碍的球面上实现实时的、稳定的无粘性流体动力学和气动声学模拟。这一目标通过结合 CPM、半拉格朗日对流和隐式 SDF 障碍物表示来实现。

4.2 针对特定方法的详细比较

4.2.1 相比于传统球坐标网格方法 (Latitude-Longitude Grids)

  • 先前工作：如 Yang 等人（参考文献 [2]）在球面上使用经纬度网格进行流体模拟。
  • 局限性：
    1. 极点奇异性：在两极网格点无限聚集，导致数值不稳定和计算资源浪费。
    2. 非均匀性：网格间距随纬度变化，导致各向异性，影响空间导数的精度，特别是在流动以平流为主时。
  • 本文的改进：
    • 消除奇异性：使用 CPM 和笛卡尔嵌入网格，完全消除了坐标奇异性，保证了计算域内的均匀分辨率。
    • 统一算子：所有微分算子（梯度、拉普拉斯）在笛卡尔网格上使用标准的有限差分公式，无需引入度量张量修正，简化了实现并提高了稳定性。

4.2.2 相比于传统贴体网格/非结构网格方法 (Body-Fitted/Unstructured Meshes)

  • 先前工作：传统的 CFD 求解器通常使用贴体网格来精确捕捉复杂几何形状。
  • 局限性：
    1. 网格生成成本：生成高质量的贴体网格在计算上是昂贵的，且难以实时更新。
    2. 障碍物移动/添加困难：如果流场中需要添加动态障碍物，通常需要重新生成网格（Remeshing），这对于实时交互应用是不可接受的。
  • 本文的改进：
    • 隐式几何：引入 SDF（符号距离函数） 来建模障碍物。障碍物的形状通过标量场定义，而非显式的网格拓扑。
    • 无需重网格化：当障碍物位置或形状改变时，只需更新 SDF 值，而不需要改变底层的流体计算网格（窄带笛卡尔网格）。这极大地提升了交互性和灵活性。
    • 边界条件处理：通过在网格点上强制速度为零或投影到切平面来处理无滑移/无穿透条件，这种方法比在非结构网格上处理复杂边界条件更简单、更鲁棒。

4.2.3 相比于 Catmull-Clark 细分曲面

  • 先前工作：Stam 等人利用细分曲面进行图形学渲染。
  • 局限性：细分曲面主要用于几何建模和渲染，虽然提供了光滑的几何表示，但并不直接支持 PDE 求解器。在细分曲面上求解流体方程通常需要将场投影回参数域或处理复杂的度量张量。
  • 本文的改进：虽然论文提到了结合细分曲面进行几何表示，但在求解层面，论文指出 CPM 提供了更直接的 PDE 求解路径。CPM 允许在不进行参数化的情况下直接在光滑几何上定义算子，相比细分曲面方法，它在数值离散化方面更具优势，且更容易与标准的有限差分格式结合。

4.2.4 相比于传统气动声学方法 (Computational Aeroacoustics, CAA)

  • 先前工作：传统的 CAA 方法（如直接数值模拟 DNS 或线性化欧拉方程）通常需要高分辨率的体积网格来捕捉声波的传播，计算成本极高，难以实时运行。
  • 局限性：对于实时应用，计算远场声压通常过于昂贵。
  • 本文的改进：
    • FWH 类比与紧凑源假设：论文采用 Ffowcs Williams-Hawkings 类比，将声源简化为表面力的导数。这种“混合方法”将流体求解（近场）与声传播解耦。
    • 实时音频合成：并非直接求解声波方程，而是将计算出的气动力谱转化为驱动正弦波振荡器的参数。这种方法牺牲了严格的声传播物理（如衍射、反射），但换来了极高的计算效率，使得实时生成听感上合理的声学信号成为可能。这与需要完整求解波动方程的高保真 CAA 方法形成了鲜明对比。

4.2.5 相比于 Stam 的 Stable Fluids

  • 先前工作：Stam 提出的 Stable Fluids 方法是实时流体模拟的里程碑，利用半拉格朗日方法实现了大时间步长的稳定性。
  • 局限性：Stam 的方法主要针对平面欧几里得域。
  • 本文的改进：论文将 Stam 的稳定对流技术与 CPM 相结合，并将投影方法扩展到曲面。这不仅保留了 Stam 方法的数值稳定性，还将其推广到了复杂的球面几何结构中，解决了在曲面上维持质量守恒和处理几何约束的问题。

4.3 总结

本文相对于先前工作，核心优势在于**“实时性”与“复杂几何”**的结合。通过 CPM 解决了几何表示的数学难题，通过 SDF 解决了障碍物处理的灵活性难题，通过 半拉格朗日投影法 保证了数值稳定性，最后通过 参数化音频合成 解决了声学计算的效率难题。这一系列的技术集成，使得在标准硬件上实现球面上的流体-声学交互模拟成为可能。

3. Riesz算子在$L^\infty$中无界性的新证明及其在Euler型方程$W^{1,\infty}$中的适度不适定性中的应用

A new proof of unboundedness of Riesz operator in $L^\infty$ and applications to mild ill-posedness in $W^{1,\infty}$ of the Euler type equations

作者: Jinlu Li, Yanghai Yu

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15922v1

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

详细分析

中文标题: Riesz算子在$L^\infty$中无界性的新证明及其在Euler型方程$W^{1,\infty}$中的适度不适定性中的应用

1. 研究对象和背景

本论文主要研究流体力学中描述不可压流体运动的欧拉型方程的适度不适定性，其核心理论工具是Riesz变换在无界函数空间中的性质。文章首先研究了Riesz变换在$L^\infty(\mathbb{R}^d)$空间上的无界性，这一性质与经典的Calderon-Zygmund理论（要求$p>1$）有着本质区别。基于这一性质，作者将其应用于两个具体的物理模型：三维旋转欧拉方程和二维带有部分阻尼的欧拉方程，证明了它们在Lipschitz函数空间$W^{1,\infty}$中是适度不适定的。

1.1 方程与系统描述

论文在引言部分详细介绍了研究对象，具体数学公式如下：

  1. 广义的欧拉型方程组： 考虑如下柯西问题：

    $$ \begin{cases} \partial_t u + u \cdot \nabla u + \nabla p = A u, & (t,x) \in \mathbb{R}^+ \times \mathbb{R}^d, \\ \operatorname{div} u = 0, & (t,x) \in \mathbb{R}^+ \times \mathbb{R}^d, \\ u(0,x) = u_0(x), & x \in \mathbb{R}^d, \end{cases} $$

    其中 $A$ 是一个常数矩阵。

  2. 三维旋转欧拉方程： 当 $A u = -\Omega e_3 \times u$ 且 $d=3$ 时，上述方程组退化为带有科里奥利力的旋转欧拉方程：

    $$ \begin{cases} \partial_t u + \Omega e_3 \times u + u \cdot \nabla u + \nabla P = 0, & (t,x) \in \mathbb{R}^+ \times \mathbb{R}^3, \\ \operatorname{div} u = 0, \\ u(0,x) = u_0(x), \end{cases} $$

    其中 $\Omega \in \mathbb{R}$ 是旋转速度（科里奥利参数），$e_3 = (0,0,1)$ 是垂直单位向量。

  3. 二维带有部分阻尼的欧拉方程： 当 $d=2$ 且 $A$ 为特定矩阵时，考虑仅在速度方程的第一分量上有阻尼的情形：

    $$ \begin{cases} \partial_t u + u \cdot \nabla u + \nabla p = (-u_1, 0), & (t,x) \in \mathbb{R}^+ \times \mathbb{R}^2, \\ \operatorname{div} u = 0, \\ u(0,x) = u_0(x). \end{cases} $$

1.2 函数空间与物理背景

文章研究的核心区域是Lipschitz函数空间，具体定义为 $W^{1,\infty}(\mathbb{R}^d)$，其范数为 $\|f\|_{W^{1,\infty}} = \|f\|_{L^\infty} + \|\nabla f\|_{L^\infty}$。对于不可压流体，速度场 $u$ 属于该空间意味着流体具有有界的Lipschitz正则性。这与经典的柯尔莫哥洛夫湍流理论中的对数正态分布以及局部适定性研究中常用的 Sobolev 空间 $H^s (s > d/2 + 1)$ 截然不同。$W^{1,\infty}$ 位于 $L^\infty$ 的框架下，而 Riesz 变换在 $L^\infty$ 上的有界性失效，这导致了 PDE 理论中的根本困难。

  • 物理背景：
    • 旋转欧拉方程：模拟了地球物理流体动力学中的大尺度流动，如大气和海洋中的环流。科里奥利力的引入带来了色散效应，这通常被认为有助于解的平滑化或长时间存在性。然而，作者指出这种色散效应具有各向异性和退化性，并未保证在 $L^\infty$ 框架下的适定性。
    • 带阻尼的欧拉方程：阻尼项 $(-u_1, 0)$ 模拟了流体与底部边界或耗散介质的相互作用，导致系统能量衰减。然而，研究表明这种局部阻尼破坏了传统的全局适定性证明框架（如Yudovich理论），在临界空间中仍可能导致不适定性。

2. 主要定理或主要结果

本节列举文章证明的核心定理和命题。

2.1 线性部分：Riesz变换的无界性

命题 1.2 (Proposition 1.2)：对于 $d \ge 1$，令 $1 \le i \le d$，设 $f_n(x)$ 由特定构造给出。则有：

$$ \|f_n\|_{L^\infty(\mathbb{R}^d)} \le C \Gamma_n, \quad \|R_i f_n\|_{L^\infty(\mathbb{R}^d)} \approx \Gamma_n \ln n, $$

其中 $R_i$ 是第 $i$ 个经典的 Riesz 变换。这直接表明 Riesz 变换 $R_i$ 在 $L^\infty$ 空间上是无界算子。

推论 1.1 (Corollary 1.1, 线性不适定性)：设 $d \ge 1$。线性方程

$$ \begin{cases} \partial_t f = R_i f, \\ f(0,x) = f_0(x), \end{cases} $$

在 $L^\infty(\mathbb{R}^d)$ 空间上是强不适定的。这意味着存在初始数据 $f_0$ 在 $L^\infty$ 中任意小，但其对应的解在任意短的时间 $t > 0$ 内 $L^\infty$ 范数可能大幅增长。

2.2 主要非线性结果：欧拉方程的适度不适定性

定理 1.1 (Theorem 1.1, 旋转欧拉方程)：设 $\Omega \neq 0$。存在一个普适常数 $\delta > 0$ 以及一列属于 Schwartz 空间 $S(\mathbb{R}^3)$ 的初始数据序列 $\{u_{0,n}\}$，满足：

$$ \|u_{0,n}\|_{W^{1,\infty}} \to 0 \quad \text{当 } n \to \infty, $$

但这些初始数据生成的对应光滑解序列 $\{u_n\}$ 满足：

$$ \|u_n(t_n)\|_{W^{1,\infty}} \ge \delta, \quad \text{其中 } t_n \to 0. $$

推论 1.2 (Corollary 1.2)：三维旋转欧拉方程在 $W^{1,\infty}(\mathbb{R}^3)$ 和 $C^1(\mathbb{R}^3)$ 空间中是适度不适定的。

定理 1.2 (Theorem 1.2, 带阻尼欧拉方程)：存在一个普适常数 $\delta > 0$ 以及一列属于 $S(\mathbb{R}^2)$ 的初始数据序列 $\{u_{0,n}\}$，满足：

$$ \|u_{0,n}\|_{W^{1,\infty}} \to 0, $$

但生成的对应光滑解序列 $\{u_n\}$ 满足：

$$ \|u_n(t_n)\|_{W^{1,\infty}} \ge \delta, \quad \text{其中 } t_n \to 0. $$

2.3 进一步应用

定理 6.1 (Theorem 6.1, 表面准地转方程)：考虑二维表面准地转方程，存在 $C > 0$ 和一列初始数据序列，使得 $C^1$ 范数趋于0，但对应解的 $C^1$ 范数在短时间 $t_n \to 0$ 内有下界 $C$，即方程在 $C^1$ 空间中不适定。

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

本论文的分析方法融合了傅里叶调制分析、Littlewood-Paley 理论以及非线性 PDE 的精细估计。以下是核心工具和技术的详细解析。

3.1 核心分析学工具

  1. Littlewood-Paley (L-P) 分解与 Besov 空间：

    • 文章大量使用 L-P 分解来处理函数的频率局部化。通过径向光滑函数 $\varphi, \phi$ 将函数分解为不同频率的块 $\Delta_j u$ 和 $\dot{\Delta}_j u$。
    • 定义了非齐次 Besov 空间 $B^s_{p,r}$ 和齐次 Besov 空间 $\dot{B}^s_{p,r}$。特别地，对于 $s > d/p$，嵌入 $B^s_{p,r} \hookrightarrow L^\infty$ 成立（Corollary 2.3），这为处理 $L^\infty$ 估计提供了正则性基准。
    • Bernstein 不等式 (Lemma 2.2)：用于控制局部傅里叶变换的导数与范数之间的关系，是处理频率截断误差的基础。

  2. Bony 仿积分解 (Bony’s Para-product)：

    • 由于 $L^\infty$ 中的乘积没有连续性，Bony 分解（Lemma 2.1）是处理非线性项 $u \cdot \nabla v$ 的关键工具。它将乘积分解为：
      $$ uv = T_u v + T_v u + R(u,v), $$
      其中 $T$ 是正则化项，$R$ 是剩余项。这种分解允许在不同频率之间建立精确的交换子估计。

3.2 解决技术性难题的技巧

A. Riesz 变换在 $L^\infty$ 无界性的新证明 (Proposition 1.2)

这是全文的理论基石。作者构造了一个特殊的函数序列 $f_n$，其傅里叶支撑在特定环域 $\{ \xi : 4/3 \le |\xi| \le 3/2 \}$ 上：

$$ f_n(x) = \sum_{k=1}^n \Gamma_n^{-1} \chi(2^k x), \quad \Gamma_n = \frac{1}{k \ln \ln n}. $$
  • 构造思想：利用对数对数的慢增长函数来抵消 Riesz 变换的平滑效应。
  • 估计策略：
    1. 利用对偶性或直接估计证明 $\|f_n\|_{L^\infty} \lesssim \Gamma_n$。
    2. 关键步骤在于证明 $\|R_i f_n\|_{L^\infty}$ 的下界。由于 $R_i$ 是奇异积分算子，在 $x=0$ 附近呈现奇异性。通过计算 $R_i \chi(2^k x)$ 在 $x=0$ 处的渐近行为，发现其贡献约为 $\ln 2^k \approx k$。
    3. 求和 $\sum_{k=1}^n k \cdot \frac{1}{k \ln \ln n} \approx \frac{n}{n \ln \ln n}$ 被对数项控制。作者通过精细的相位分析和 $L^\infty$ 振荡积分估计，证明了 $R_i f_n$ 的和产生了 $\ln n$ 的增长因子，从而得到 $\|R_i f_n\|_{L^\infty} \approx \Gamma_n \ln n$。这表明算子范数 $\|R_i\|_{L^\infty \to L^\infty}$ 发散。

B. 欧拉方程适度不适定性的证明 (Theorem 1.1 & 1.2)

这是将线性算子性质推广到非线性系统的核心。作者采用构造性反例的方法。

  1. 初始数据的构造： 为了将 Riesz 变换的无界性“嫁接”到欧拉方程的解算子上，作者选择初始数据与上述命题中的 $f_n$ 或 $R_i f_n$ 密切相关。 对于旋转欧拉方程，作者考虑一列高频振荡的初始数据 $u_{0,n}$，其形式为 $u_{0,n} = \text{div}(v_n)$ 或直接与 $R_i f_n$ 相关的向量场。由于 Riesz 变换对应于方程中的压力项投影 $P = I - (-\Delta)^{-1}\nabla \operatorname{div}$ 的组成部分，初始数据的微小 $W^{1,\infty}$ 范数可以通过 Riesz 变换转化为非线性项 $u \cdot \nabla u$ 中的放大效应。

  2. 能量估计的困难与突破：

    • 在 $W^{1,\infty}$ 框架下，经典的能量估计 $\|\nabla u(t)\|_{L^\infty}$ 通常依赖于估计 $\|\nabla u \cdot \nabla u\|_{L^\infty}$。在 $L^\infty$ 空间中，乘积 $\|\nabla u \cdot \nabla u\|_{L^\infty}$ 没有良好的下界控制，且 Riesz 变换在 $L^\infty$ 上的无界性导致压力项 $\nabla p$（包含 $R$ 变换）无法被控制。
    • 关键技巧：作者利用了Magnus 效应或色散效应的逆过程。对于旋转欧拉方程，色散项 $\Omega e_3 \times u$ 通常用于衰减解的范数，但在临界情况下，其各向异性可能导致解在特定方向上的梯度集中。
    • 证明中，作者没有直接求解方程，而是通过共轭变换或格林函数法（如果可能）来追踪解的动力学。在本论文中，作者主要依赖于对线性部分（包含 Riesz 变换）的显式估计。通过将非线性项视为微扰，或者在特定的时间尺度 $t_n \to 0$ 上，利用解的 Duhamel 公式：
      $$ u(t) = e^{t(\text{投影项})} u_0 + \int_0^t e^{(t-s)(\text{投影项})} Q(u \cdot \nabla u)(s) ds. $$
    • 核心难点攻克：在 $W^{1,\infty}$ 中，$Q(u \cdot \nabla u)$ 包含 $(-\Delta)^{-1}\nabla \operatorname{div}(u \cdot \nabla u)$。由于 $u \cdot \nabla u$ 包含两个梯度的乘积，若 $u_0$ 构造得当，该项的散度 $\operatorname{div}(u \cdot \nabla u)$ 会包含高频率成分，其 Riesz 变换的 $L^\infty$ 估计会引入 $\ln(1/t)$ 或类似的对数发散项。作者正是通过精细的估计捕捉到了这种在短时间 $t_n \to 0$ 内的放大现象，从而证明 $\|u(t_n)\|_{W^{1,\infty}}$ 无法趋于 0。

C. Littlewood-Paley 交换子估计 (Lemma 2.3)

论文推导了交换子估计：

$$ \|[R, u \cdot \nabla] f\|_{B^s_{p,r}} \le C \|\nabla u\|_{L^\infty \cap B^{s-1}_{p,\infty}} \|f\|_{B^s_{p,r}}. $$

这项估计是技术核心之一。它表明当 $u$ 具有 $B^{s-1}_{p,\infty}$ 正则性时，Riesz 变换与输运项的交换子是可控的。然而，在 $L^\infty$ 框架下（对应 $p=\infty$），正则性 $s$ 需要非常严格。作者通过此引理说明，在临界正则性下，非线性项与 Riesz 变换的耦合会导致控制失效，从而支持了不适定性的结论。

4. 与之前工作的比较

本论文在流体力学 PDE 的适定性理论，特别是 $L^\infty$ 框架下的理论，做出了显著的推进。

4.1 理论框架的突破：从涡度方程到速度方程

  • 先前工作 (Elgindi & Masmoudi, 2016 等)： Elgindi 和 Masmoudi 在其开创性工作中，利用 $L^\infty$ 中 Riesz 变换的不连续性（通过对数奇异性），建立了 2D Euler 方程涡度 $\omega = \operatorname{curl} u$ 方程的适度不适定性。然而，他们的结果主要局限在涡度空间 $L^\infty(\omega)$。 对于物理上更直观的速度空间 $W^{1,\infty}$（即 Lipschitz 连续速度场），由于 $u = P(\text{div}(u \otimes u))$ 涉及 Riesz 变换作用于非线性项，问题变得极具挑战性。先前的结果往往依赖于涡度的 Yudovich 理论（要求 $\omega \in L^\infty$ 且 $u$ 在 Besov 空间中具有额外正则性）。

  • 本文的改进： 本论文首次在不依赖涡度公式的情况下，直接针对 3D 旋转欧拉方程在 $W^{1,\infty}$ 中的不适定性进行了研究。

    1. 维度的推广：将 $L^\infty$ 适定性分析从 2D 推广到了 3D。
    2. 空间的直接性：直接处理速度场 $u \in W^{1,\infty}$，避免了通过涡度 $\omega$ 间接控制带来的限制。这要求对压力项（即 Leray 投影 $P$）在 $L^\infty$ 中的作用有更深刻的理解。
    3. 物理机制的揭示：证明了即使在 3D 旋转欧拉方程中存在色散效应，这种色散效应也不足以克服由 Riesz 变换无界性导致的非线性增长，从而使得在 Lipschitz 空间中的解不连续。

4.2 模型的推广：旋转与阻尼效应

  • 与经典欧拉方程的比较：

    • Kato-Ponce (1988) 等经典结果在 Sobolev 空间 $W^{s,p} (1<p<\infty)$ 中建立了局部适定性，但这依赖于 Riesz 变换在 $L^p$ 上的有界性。
    • Pak-Park (2009) 证明了在 Besov 空间 $B^{1}_{\infty,1}$ 中的局部存在性，但这仍是高正则性的。
    • 本文的比较：本文关注的是临界情形 $p=\infty$。作者指出，先前关于旋转欧拉方程的工作（如 Koh-Lee-Takada, Jia-Wan 等）均基于 $L^p (1<p<\infty)$ 框架。本文首次证明了，尽管科里奥利力提供了 $O(t^{-1})$ 的 $L^\infty$ 衰减率（在没有粘性的情况下），但在 Lipschitz 范数下，解的连续依赖性依然失效。这意味着旋转不能完全“修复” $L^\infty$ 框架下的适定性问题。

  • 与带阻尼欧拉方程的比较：

    • Elgindi & Masmoudi (2016) 研究了 2D 带阻尼欧拉方程的涡度方程，证明了不适定性。
    • 本文的改进：本文将这一结论推广到了速度场 $W^{1,\infty}$。阻尼项 $(-u_1, 0)$ 虽然耗散能量，但在临界 Sobolev 指数下，其正则化效果不足以抵消 $u \cdot \nabla u$ 在 $L^\infty$ 框架下的非线性奇异性。这与 Elgindi-Masmoudi 的工作形成互补，完整刻画了该模型在临界空间的适定性面貌。

4.3 技术工具的创新：从模态分析到非线性交换子

  • 先前工具：传统的证明常使用能量方法（Benergy estimates）或 Bochner 积分。对于 $L^\infty$ 问题，标准的能量方法通常失效，因为 $L^\infty$ 范数不是 Hilbert 空间结构，且缺乏自反性。
  • 本文工具：
    1. Littlewood-Paley 精细分解：论文没有简单地引用 Riesz 变换的无界性，而是重新给出了一个基于 L-P 分解的初等证明（Proposition 1.2）。这种构造性证明不仅展示了 $\ln n$ 的增长来源，而且其频率支撑结构（在 $4/3$ 到 $3/2$ 的环域）可以直接用于非线性方程的解算子构造中。
    2. 交换子的导数损失：在证明过程中，作者处理了诸如 $[R, u \cdot \nabla] \nabla$ 这类高阶交换子。与标准的 $L^p$ 理论相比，这里的困难在于 $L^\infty$ 中的正则性损失。本文通过引入 Besov 空间的插值和 $(p, r)$ 指数的精细选择，控制了在临界正则性下的误差项，这是对现有 PDE 估计技术的重要补充。

4.4 总结

本论文的核心贡献在于：

  1. 澄清了 $L^\infty$ 框架下的机制：通过简明的 Riesz 变换无界性证明，揭示了欧拉方程在 Lipschitz 空间中不适定性的本质来源。
  2. 扩展了物理模型的范围：将 Elgindi-Masmoudi 的 2D 涡度理论推广到了 3D 旋转系统和 2D 阻尼系统，并直接处理了速度空间 $W^{1,\infty}$，这在物理上比涡度空间更自然，但在数学上更困难。
  3. 提供了普适的反例构造方法：利用频率截断和对数增长构造初始数据，为研究其他临界流体动力学方程（如表面准地转方程）在 $W^{1,\infty}$ 中的稳定性提供了通用框架。

相较于前人的工作，本文在理论深度上填补了 $L^\infty$ 适定性理论在 3D 和旋转/阻尼系统中的空白，在技术上展示了如何利用 Littlewood-Paley 理论和仿积分解来克服 $L^\infty$ 空间中乘积估计失效的难题。

4. 具有摩擦型滑移边界条件的随机强迫可压缩Navier-Stokes方程

Stochastically forced compressible Navier-Stokes equations with slip boundary conditions of friction type

作者: Reo Tsuboya

类别: math.PR, math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15768v1

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

详细分析

中文标题: 具有摩擦型滑移边界条件的随机强迫可压缩Navier-Stokes方程

1. 研究对象和背景

1.1 控制方程与物理模型

本文研究由随机外力驱动的可压缩粘性流体在有界光滑域上的运动，该流体服从具有摩擦型滑移边界条件的模型。控制方程是随机可压缩Navier-Stokes方程。设 $D \subset \mathbb{R}^3$ 为一个有界开集，具有局部Lipschitz边界 $\partial D$，$n$ 为边界的单位外法向量。时间区间为 $(0,T)$。未知量为流体的密度 $\rho: (0,T) \times D \to \mathbb{R}$ 和速度场 $u: (0,T) \times D \to \mathbb{R}^3$。

方程组由连续性方程、动量方程、初值条件及边界条件组成：

  1. 连续性方程 (质量守恒)：

    $$ \mathrm{d}\rho + \mathrm{div}(\rho u) \mathrm{d}t = 0 \quad \text{in } (0,T) \times D $$

  2. 动量方程：

    $$ \mathrm{d}(\rho u) + \mathrm{div}(\rho u \otimes u) \mathrm{d}t = \mathrm{div} \mathbb{T} \mathrm{d}t + G(\rho, \rho u) \mathrm{d}W \quad \text{in } (0,T) \times D $$

    其中 $\mathbb{T}$ 为牛顿粘性流体的Cauchy应力张量，对于等熵流体，其构型为：

    $$ \mathbb{T} = \mathbb{S}(\nabla u) - p(\rho)I $$

    $$ \mathbb{S}(\nabla u) = \mu \left( \nabla u + (\nabla u)^T - \frac{2}{3} (\mathrm{div} u) I \right) + \lambda (\mathrm{div} u) I, \quad p(\rho) = a \rho^\gamma $$

    这里 $\mu > 0$ 和 $\lambda \ge 0$ 是粘性系数，$\gamma > 1$ 是绝热指数，$a > 0$ 是常数，$I$ 是单位矩阵。

  3. 初值条件：

    $$ (\rho(0), (\rho u)(0)) = (\rho_0, q_0) \quad \text{in } D $$

  4. 边界条件 (摩擦型滑移条件)：

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

    $$ |(\mathbb{T} n)_\tau| \le g, \quad (\mathbb{T} n)_\tau \cdot u_\tau + g |u_\tau| = 0 \quad \text{on } (0,T) \times \partial D $$

    这里 $u_n = u \cdot n$，$u_\tau = u - (u \cdot n)n$ 为切向速度，$(\mathbb{T} n)_\tau$ 为切向应力，$g = g(x, t) \ge 0$ 是摩擦系数（也称为摩擦模量）。

    驱动过程 $W$ 是定义在概率空间 $(\Omega, \mathcal{F}, (\mathcal{F}_t)_{t \ge 0}, \mathbb{P})$ 上的柱状Wiener过程（Cylindrical Wiener Process），取值于可分Hilbert空间 $U$。系数 $G$ 是非线性算子，满足特定的增长条件，例如：

    $$ |G_k(x, \rho, q)| \le g_k (\rho + |q|) $$

    $$ |\nabla_{\rho, q} G_k(x, \rho, q)| \le g_k $$

1.2 研究背景与动机

Navier-Stokes方程在各种边界条件下已被广泛研究。在确定性设定下，弱解的存在性理论已相当成熟（参见Lions的经典著作）。对于可压缩流体，传统的边界条件通常是无滑移条件（No-slip）或周期性边界条件。然而，无滑移条件在某些物理情境下会导致悖论（Hess, Hil），因此非标准边界条件（如滑移边界条件）近年来受到更多关注。

Navier滑移边界条件 允许切向滑移，且切向滑移与切向应力成正比。在确定性设定下，关于弱解的存在性已有很多结果（如[FKN+13]）。另一种允许大应力下滑移的边界条件是摩擦型滑移边界条件（Slip boundary conditions of friction type），它引入了摩擦不等式。对于不可压缩流体，已有关于强解局部/全局存在性的研究；对于可压缩流体，确定性情形下的弱解存在性最近由[Nečasová, Ogorzaly, Scherz]（简称[NOS23]）建立。

随机设定下的挑战：对于随机可压缩Navier-Stokes方程，已有工作主要集中在无滑移或周期性边界条件（如[BH16], [Smi17]）。然而，关于滑移型边界条件在概率框架下的弱解理论存在空白。本文旨在填补这一空白，在有界光滑域 $C^{2+\nu}$ ($\nu > 0$) 上建立具有摩擦型滑移边界条件的随机可压缩Navier-Stokes方程的弱解存在性。

1.3 空间设定与正则性要求

  • 定义域：$D \subset \mathbb{R}^3$ 是有界域，边界 $\partial D$ 具有 $C^{2+\nu}$ 正则性。这种正则性是证明迹定理和边界近似算子（如Bogovskii算子）性质所必需的。
  • 函数空间：
    • 密度 $\rho$ 的主要能量空间：$L^\infty(0,T; L^\gamma(D))$，其中 $\gamma > 3/2$（在主定理中要求 $\gamma > 3/2$）。
    • 速度 $u$ 的主要能量空间：$L^2(0,T; W^{1,2}(D))$，且满足 $u \cdot n|_{\partial D} = 0$。
    • 动量 $\rho u$ 的空间：$L^\infty(0,T; L^{\frac{2\gamma}{\gamma+1}}(D))$。
    • 能量的下界控制：为了处理随机积分和强解的等价性，常假设密度有正下界 $\rho \ge \underline{\rho} > 0$。

2. 主要定理或主要结果

论文的主要结果是Theorem 1.3，即随机强迫可压缩Navier-Stokes方程在摩擦型滑移边界条件下**耗散鞅解（Dissipative Martingale Solution）**的存在性。

定理 1.3 (存在性定理)

设 $T > 0$，$D \subset \mathbb{R}^3$ 是 $C^{2+\nu}$ 类的有界域（$\nu > 0$）。令 $\gamma > 3/2$，$\Lambda$ 是定义在 $L^1(D) \times L^1(D)$ 上的Borel概率测度，满足：

  1. $\Lambda\{\rho \ge 0\} = 1$。
  2. 存在常数 $\underline{\rho}, \bar{\rho} > 0$ 使得 $\Lambda\{ \int_D \rho \mathrm{d}x \le \bar{\rho} \} = 1$。
  3. 初始能量有限：$\int_{L^1 \times L^1} \left( \int_D \left( \frac{1}{2} \frac{|q|^2}{\rho} + P(\rho) \right) \mathrm{d}x \right)^r \mathrm{d}\Lambda(\rho, q) < \infty$，对某个 $r \ge 4$。其中 $P(\rho) = \rho \int_1^\rho \frac{p(z)}{z^2} \mathrm{d}z$ 是压力势。

假设扩散系数 $G = G(x, \rho, q)$ 连续可微，满足增长条件 (1.6)-(1.7)，摩擦系数 $g \in L^2((0,T) \times \partial D)$ 满足 $g \ge 0$ 且 $\int_0^T \int_{\partial D} g \mathrm{d}\Gamma \mathrm{d}t > 0$。

结论：存在一个耗散鞅解 $((\Omega, \mathcal{F}, (\mathcal{F}_t)_{t \ge 0}, \mathbb{P}), \rho, u, W)$ 满足定义 1.1，且解满足如下能量估计：

$$ \mathbb{E} \left[ \sup_{t \in [0,T]} \|\rho(t)\|_{L^\gamma}^r \right] + \mathbb{E} \left[ \|u\|_{L^2(0,T; W^{1,2})}^r \right] < \infty $$

$$ \mathbb{E} \left[ \sup_{t \in [0,T]} \|\rho |u|^2(t)\|_{L^1}^r \right] + \mathbb{E} \left[ \sup_{t \in [0,T]} \|\rho u(t)\|_{L^{\frac{2\gamma}{\gamma+1}}}^r \right] < \infty $$

附注：强解的等价性 (Proposition 1.5)

为了建立弱解理论的适定性，论文还证明了在密度有正下界且解具有足够正则性的情况下，弱解（Definition 1.1）与解析强解（Definition 1.4）是等价的。强解定义为 $(\rho, u) \in C([0,T]; W^{s,2})$ 且 $u \in L^2(0,T; W^{s+1,2})$ ($s > 3/2$)，满足方程的强形式及边界条件 $(\mathbb{T}n)_\tau + g u_\tau = 0$（这里应该是力的平衡，原文有详细推导）。

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

本文的证明策略融合了确定性流体力学中的近似方法与随机偏微分方程中的紧性技术。核心难点在于：1) 边界条件的非线性（摩擦不等式）；2) 随机项带来的低正则性；3) 可压缩性导致的质量守恒与动量方程的耦合。

3.1 整体框架：四层近似方案 (Four-Layer Approximation Scheme)

论文采用并扩展了 [BFH18b] 中提出的四层近似方案，该方案最初是为周期性边界条件设计的。对于滑移边界条件，作者对其进行了重要的修改和扩展。主要的近似参数包括：

  1. 时间离散化：将时间区间分割，构造时间离散逼近解。
  2. 截断与正则化：
    • 密度截断：引入参数 $\alpha$，处理方程在 $\rho=0$ 附近的退化。
    • 人工压力：引入参数 $\beta$，添加人工压力项 $\beta \rho^\beta$ 以保证密度的正性（在变分法框架下，通常需处理 $\rho^\gamma$ 的非凸性，这里使用 $p(\rho) = a\rho^\gamma + \beta \rho^\beta$，其中 $\beta > \gamma$）。
    • 人工粘性：引入参数 $\varepsilon$，添加 $\varepsilon$-Laplacian 以获得更强的紧性。
  3. 噪声近似：
    • 有限维逼近：由于 $W$ 是柱状过程，需通过正交投影 $P_m$ 逼近噪声项。
    • 系数光滑化：对扩散系数 $G$ 进行正则化，使其满足Lipschitz条件（在逼近层次上）。

证明路线图： 作者在 Section 3.1 中概述了证明路径：首先求解基本逼近问题（参数 $\varepsilon, \alpha, \beta$ 和 $R$ 截断），证明其适定性（Section 3.2, 3.3）。然后依次取极限：

  • $R \to \infty$（截断消除，Section 3.4）。
  • $\alpha \to 0$（消除人工阻尼，Section 4）。
  • Galerkin 近似参数极限（Section 5）。
  • $\varepsilon \to 0$（人工粘性消失，Section 6）。
  • $\beta \to 0$（人工压力消失，Section 7）。

3.2 核心分析工具与技术细节

3.2.1 能量估计与压缩性控制

在可压缩流体中，控制密度 $\rho$ 是难点。作者利用了压力势 $P(\rho)$：

$$ P(\rho) = \rho \int_1^\rho \frac{a z^{\gamma-1}}{z^2} \mathrm{d}z $$

通过乘以测试函数并积分，结合随机积分的Itô公式（见 Proposition 3.16），得到形如：

$$ \mathbb{E} \left[ \sup_t \int_D \left( \frac{1}{2}\rho|u|^2 + P(\rho) \right) \mathrm{d}x + \int_0^T \int_D |\nabla u|^2 \mathrm{d}x \mathrm{d}t \right] \le C $$

的关键先验估计。

3.2.2 处理摩擦型滑移边界条件

这是本文技术上的主要创新点之一。边界条件 $|(\mathbb{T}n)_\tau| \le g$ 和 $(\mathbb{T}n)_\tau \cdot u_\tau + g|u_\tau| = 0$ 包含障碍（Barrier）结构。

  • 确定性技巧：在[NOS23]中，利用凸逼近方法处理绝对值函数 $|u|$。在本论文中，作者需要在随机背景下实现这一点。
  • 变分不等式的导出： 作者在定义 1.1 中结合了动量方程和能量不等式，构造了一个动量和能量不等式 (Momentum and Energy Inequality)，即方程 (1.10)。该不等式在形式上包含了边界项：
    $$ \int_0^\tau \int_{\partial D} \left( g|\phi \varphi| - g|u| \right) \mathrm{d}\Gamma \mathrm{d}t + \dots \ge 0 $$
    其中 $\phi$ 是截断函数，$\varphi$ 是测试函数。这本质上是一个变分不等式，反映了边界上的摩擦阻力。
  • 凸逼近方法： 在 Section 4 中，为了处理 $\alpha \to 0$（对应动量方程中的阻尼项消失），作者引入了凸逼近。困难在于，当速度 $u$ 的确定性界丢失（仅能控制能量），且边界条件非线性时，如何保持紧性。作者利用了 [NOS23] 中处理绝对值函数的凸逼近技术，结合随机紧性理论，证明了逼近解的收敛性。

3.2.3 噪声系数的逼近与有限维极限

在 Section 3.2 和 Section 5 中，作者面临处理柱状噪声 $W$ 的问题。

  • 标准困难：通常使用正交投影 $P_m$ 将方程投影到有限维空间。然而，在滑移边界条件下，$P_m$（通常是 $L^2$ 投影）在 Sobolev 空间 $W^{1,p}$ 中的连续性依赖于维度，这在 $m \to \infty$ 时难以控制。
  • 本文技术：作者结合了 [Smi17] 的思想，对扩散系数 $G(\rho, \rho u)$ 进行了双重逼近：
    1. 正则化：$G_\varepsilon$，使其在空间上 Lipschitz 连续。
    2. 有限维投影：$[G_\varepsilon(\rho, \rho u)]^*_m$（投影到有限维子空间）。 关键在于利用 $C_c^\infty(D)$ 在 $W^{1,p}(D)$ 中的稠密性，避免了直接依赖投影算子在 Sobolev 范数下的范数连续性估计，从而绕过了维数依赖常数的问题。

3.2.4 密度的强收敛与有效粘性通量 (Effective Viscous Flux)

为了获得密度 $\rho$ 的强收敛（从而处理非线性压力项 $p(\rho)$），作者使用了经典的有效粘性通量技术，由 DiPerna-Lions 引入，并在随机设定下由 [FN09] 等发展。

  • 工具：对于可压缩流，有效粘性通量定义为 $V = (\mathrm{div} \mathbb{S}(\nabla u) - p(\rho))I - \mathbb{S}(\nabla u)$ 的某种形式（或其分布意义下的极限）。
  • 适用性验证：在 Section 6 和 7 中，作者需要验证在滑移边界条件下，强解定义的 Reynolds 输运定理 和 有效粘性通量恒等式 依然成立。这要求对边界项进行细致的处理，确保边界摩擦力不破坏内部的正则性机制。
  • 收敛策略：通过 Aubin-Lions-Simon 引理和有效粘性通量的弱连续性，证明 $\rho_n \to \rho$ 在 $L^1_{\mathrm{loc}}$ 意义下强收敛，进而通过能量法控制弱收敛误差。

3.2.5 Korn-Poincaré 不等式的加权推广

在滑移边界条件下，速度的控制比无滑移条件更困难，因为切向速度在边界上非零。

  • 挑战：标准的 Poincaré 不等式不适用，且通常的能量估计不能直接控制 $\nabla u$ 的全范数。
  • 解决方案：作者引用并推广了带有密度或摩擦系数作为权重的广义 Korn-Poincaré 不等式（Corollary 3.22, Lemma 3.21）。这使得即使在 $\rho$ 随机变化的情况下，也能结合粘性应力项 $\mathbb{S}(\nabla u)$ 和边界摩擦项来控制速度场的正则性。这要求对摩擦系数 $g$ 施加一定的条件（如 $g \in L^2$ 且 $\int g > 0$）。

4. 与之前工作的比较

本文在随机可压缩流体力学和滑移边界条件的交叉领域做出了显著的贡献，主要体现在以下几个方面：

4.1 边界条件的推广：从无滑移到摩擦型滑移

  • 先前工作：
    • 确定性流体：Navier滑移、完全滑移已有研究（[FKN+13], [NOS23]）。
    • 随机流体：已知结果主要集中在无滑移边界（[Smi17]）或周期性边界（[BH16]）。
  • 本文突破：这是第一篇在随机设定下处理摩擦型滑移边界条件（Slip boundary conditions of friction type）的论文。这种边界条件引入了障碍问题结构（$| \cdot | \le g$），比线性的Navier滑移条件更难处理，因为动量方程的解与边界摩擦力之间存在非线性耦合。
    • 技术意义：作者成功地将确定性情形下的凸逼近方法（[NOS23]）移植到随机背景下，并解决了随机积分带来的额外低正则性问题。

4.2 域的几何约束：从环面到有界光滑域

  • 先前工作：在随机流体领域，许多基础工作（如 [BFH18b]）基于三维环面 $\mathbb{T}^3$。虽然环面在分析上（傅里叶分析、无边界）很方便，但它缺乏物理真实性。
  • 本文突破：论文在有界光滑域 $D \subset \mathbb{R}^3$（$C^{2+\nu}$ 正则性）上建立了结果。
    • 困难：有界域上的 Sobolev 嵌入和紧性依赖于域的几何形状，且边界项的处理（迹定理、Bogovskii 算子）比环面情形复杂得多。作者在 Section 2 和 Appendix A 中详细处理了这些工具，使其适应滑移边界条件的弱解框架。

4.3 概率框架下的适定性与正则性

  • 先前工作：随机可压缩Navier-Stokes方程的弱解理论在无滑移条件下已较成熟，但在滑移条件下几乎是空白。
  • 本文突破：
    1. 解的概念：定义了兼容于摩擦边界条件的耗散鞅解。这不仅要求解满足方程的分布形式，还要求满足一个特定的变分不等式（包含边界摩擦项），这在随机背景下是新的尝试。
    2. 强-弱等价性：证明了在密度正则且有下界的情况下（$\rho \ge \underline{\rho} > 0$），所定义的弱解与解析强解（Definition 1.4）等价。这填补了随机滑移流体中强解理论的空白，并为验证数值方法提供了理论基准。
    3. 不加额外的高阶正则性假设：与一些仅在高正则性空间（如 $W^{s,p}, s$ 很大）中证明局部解的工作不同，本文在 $\gamma > 3/2$ 的自然能量空间中证明了全局存在性，这是可压缩流体力学中具有挑战性的结果。

4.4 技术工具的创新组合

  • 逼近策略：本文展示了如何将基于投影的噪声逼近（Smi17）与基于凸性的边界条件逼近（Nečasová et al.）结合。这种结合在以前的文献中没有出现过，因为它需要同时处理随机项的弱收敛和边界障碍项的非线性收敛。
  • 通用性：论文最后（Section 1.5）指出，其框架可以扩展到Navier滑移边界条件 $(\mathbb{S}(\nabla u)n)_\tau + g u_\tau = 0$。这表明作者发展的一套先验估计（特别是加权 Korn 不等式和边界项处理）具有很好的通用性，不仅限于摩擦型滑移，也适用于更一般的滑移模型。

综上所述，Reo Tsuboya 的这篇论文在随机偏微分方程领域是一项扎实且具有创新性的工作，它将确定性流体力学中成熟的滑移边界处理技术与现代随机分析相结合，解决了有界光滑域上随机可压缩流体在复杂边界条件下的弱解存在性这一开放性问题。

5. 纳维-斯托克斯方程的整体正则性及其在欧拉方程整体可解性中的应用

Global regularity for the Navier-Stokes equations with application to global solvability for the Euler equations

作者: Myong-Hwan Ri

类别: math.AP, math-ph

发布日期: 2026-01-22

arXiv ID: 2601.15685v1

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

详细分析

中文标题: 纳维-斯托克斯方程的整体正则性及其在欧拉方程整体可解性中的应用

研究对象和背景

本文研究的是高维（$d \geq 3$）不可压缩流体的运动规律，具体关注两个紧密相关的方程组：纳维-斯托克斯方程（Navier-Stokes equations）与欧拉方程（Euler equations）。

1. 纳维-斯托克斯方程 (Navier-Stokes Equations)

论文的核心研究对象是 $d$ 维欧氏空间 $\mathbb{R}^d$（$d \geq 3$）上的纳维-斯托克斯方程的柯西问题：

$$ \begin{cases} \partial_t u - \nu \Delta u + (u \cdot \nabla) u + \nabla p = 0, & \text{in } (0, \infty) \times \mathbb{R}^d, \\ \text{div } u = 0, & \text{in } (0, \infty) \times \mathbb{R}^d, \\ u(0, x) = u_0(x), & \text{in } \mathbb{R}^d. \end{cases} $$

这里，$u = u(t, x)$ 代表流体的速度场，$p = p(t, x)$ 代表压力场，$\nu > 0$ 是流体的运动粘性系数（kinematic viscosity）。该方程组由质量守恒（不可压缩条件 $\text{div } u = 0$）和动量守恒（动量方程）导出。

背景与历史沿革： 纳维-斯托克斯方程是流体力学中最基本的偏微分方程模型。1933年，Jean Leray 证明了对于三维空间中的任意 $L^2$ 初始速度场，存在全局弱解（Leray-Hopf weak solution）。然而，关于这些弱解是否光滑（正则性）以及是否在有限时间内出现奇点（blow-up），即著名的纳维-斯托克斯方程正则性问题，至今仍是克莱数学研究所的千禧年大奖难题之一。对于 $d=3$ 的情形，已知如果弱解满足 $u \in L^\infty(0, T; L^3(\mathbb{R}^3))$（临界空间），则解是光滑的。

函数空间与尺度不变性： 文章特别关注方程的尺度不变性（scaling invariance）：

$$ u_\lambda(t, x) = \lambda u(\lambda^2 t, \lambda x). $$

在该变换下，范数保持不变的空间被称为临界空间（critical spaces）。对于 $d$ 维方程，临界齐次索伯列夫空间（homogeneous Sobolev space）是 $\dot{H}^{-1 + d/2}(\mathbb{R}^d)$。如果初始数据 $u_0$ 的正则性高于该临界指标（即 $s > -1 + d/2$），问题被称为亚临界（subcritical）；如果正则性低于该指标，则为超临界（supercritical）。

2. 欧拉方程 (Euler Equations)

作为纳维-斯托克斯方程的应用，文章也研究了不可压缩理想流体的欧拉方程：

$$ \begin{cases} \partial_t u + (u \cdot \nabla) u + \nabla p = 0, & \text{in } (0, \infty) \times \mathbb{R}^d, \\ \text{div } u = 0, & \text{in } (0, \infty) \times \mathbb{R}^d, \\ u(0, x) = u_0(x), & \text{in } \mathbb{R}^d. \end{cases} $$

这可以看作是纳维-斯托克斯方程中粘性系数 $\nu = 0$ 的情形。

3. 研究区域与边值条件

文章研究的是全空间 $\mathbb{R}^d$ 上的柯西问题（Cauchy problem），不涉及具体的边界，属于无界域上的局部适定性与全局适定性问题。解空间主要涉及索伯列夫空间 $H^s(\mathbb{R}^d)$ 及其齐次版本 $\dot{H}^s(\mathbb{R}^d)$。

主要定理或主要结果

本文的主要定理建立了高维纳维-斯托克斯方程 Leray-Hopf 弱解在更宽的正则性指数 $s$ 下的整体正则性，并由此推导出欧拉方程的整体可解性。

定理 1.1 (纳维-斯托克斯方程的整体正则性) 设 $u$ 是 $d$ 维 ($d \geq 3$) 纳维-斯托克斯方程 (1.1) 的 Leray-Hopf 弱解，初始速度 $u_0 \in H^s(\mathbb{R}^d)$，其中 $s \geq -1 + \frac{d}{2}$，且 $\text{div } u_0 = 0$。则 $u$ 具有以下性质：

$$ u \in L^\infty(0, \infty; H^s(\mathbb{R}^d)) \cap L^2(0, \infty; H^{s+1}(\mathbb{R}^d)). $$

此外，存在仅依赖于 $s$ 和维数 $d$，且与粘性系数 $\nu$ 无关的常数 $C(s) > 0$，使得对于所有 $t \in (0, \infty)$，以下估计成立：

$$ \|u(t)\|_{H^s(\mathbb{R}^d)}^2 + \nu \int_0^t \|u(\tau)\|_{H^{s+1}(\mathbb{R}^d)}^2 d\tau \leq C(s) \|u_0\|_{H^s(\mathbb{R}^d)}^2. $$

特别地，该弱解是全局正则的（即解是光滑的）且在给定的函数空间中是唯一的。

推论 (欧拉方程的整体可解性) 设 $u_0 \in H^s(\mathbb{R}^d)$，$s \geq -1 + \frac{d}{2}$，且 $\text{div } u_0 = 0$。则不可压缩欧拉方程存在全局解 $u \in C([0, \infty), H^s(\mathbb{R}^d))$。当 $s > 1 + \frac{d}{2}$ 时，该解是唯一的。

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

本文的核心贡献在于构造了一种新的加权空间并结合了重缩放论证（re-scaling argument）。证明过程主要分为三个步骤：构造超临界空间、建立高频部分的能量估计、利用重缩放和能量叠加获得独立于粘性的估计。

a) 核心分析学工具

  1. 频率分解与Littlewood-Paley理论： 文章利用傅里叶变换将解在频率空间进行分解。记 $\Delta_j$ 为频率在 $[2^{j-1}, 2^j)$ 附近的投影算子。解 $u$ 被分解为低频部分 $u_{<k}$ 和高频部分 $u_k$（在文中常指 $u_{k, k+1}$ 频带）。这使得研究者可以分别处理方程在不同尺度下的行为。

  2. 能量估计 (Energy Estimates)： 这是处理耗散型方程的标准工具。通过在动量方程中乘以解 $u$ 或其频率分量 $u_k$ 并进行积分，利用内积 $((u \cdot \nabla)u, u) = 0$（在 $\text{div } u = 0$ 条件下）带来的对流项抵消性质，导出关于 $L^2$ 能量的微分不等式。

  3. 尺度变换 (Re-scaling Argument)： 这是本文将局部正则性转化为全局正则性的关键工具。作者引入了一族依赖于时间 $t$ 的伸缩变换 $u_\lambda(t, x) = \lambda u(\lambda^2 t, \lambda x)$。通过选择合适的缩放参数 $\lambda$，可以控制解在特定时间区间上的范数（特别是构造的超临界范数），使其保持在一个足够小的常数范围内。

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

  1. 构造超临界空间与对数加权： 传统的临界空间是 $\dot{H}^{-1+d/2}$。为了处理更广泛的初始数据（即 $s$ 接近 $-1+d/2$），作者构造了一个新的空间 $X_1$。该空间的范数定义为：

    $$ \|v\|_{X_1} = \left( \sum_{j \in \mathbb{Z}} \left( 2^{j(-1+d/2)} a(j)^{-1} \|\Delta_j v\|_{L^2} \right)^2 \right)^{1/2}, $$

    其中系数 $a(j)$ 是一个具有稀疏对数权重的序列：

    $$ a(j) = \begin{cases} \log j & \text{若 } j = i + 2^{2^k} \text{ 对某些 } k \in \mathbb{N}, i \in \mathbb{Z}, -k \leq i \leq k, \\ 1 & \text{其他情况}. \end{cases} $$

    与 $\dot{H}^{-1+d/2}$ 范数相比，$X_1$ 范数在大多数频率上是相同的，仅在一系列极其稀疏的频率尺度（$j = 2^{2^k}$ 附近）上引入了 $\frac{1}{\sqrt{\log j}}$ 的衰减因子。这使得 $X_1$ 是一个超临界空间（比临界空间稍弱），且满足特定的“稀疏性”条件。

  2. 高频部分的能量估计与系数衰减： 证明的核心在于对高频部分 $u_k$（或 $u_{1}$ 表示高频率累积）建立能量不等式。通过测试高频部分，非线性项 $(u \cdot \nabla)u_k$ 被重新估计。 文中指出，通过精细观察，可以得到如下形式的估计：

    $$ \frac{d}{dt} \|u_1\|_{H^{1/2}}^2 + \nu \|\nabla u_1\|_{H^{1/2}}^2 \lesssim c(k) \|u\|_{X_1} \|\nabla u_{1/2}\|_{H^{1/2}}^2. $$

    这里的关键是处理系数 $c(k)$。如果 $c(k)$ 增长过快，直接求和会导致发散。作者利用 $a(j)$ 的构造和引理 3.1 (定义了 $b(j)$)，证明了对于频率分量的系数求和满足：

    $$ \sum_{k=1}^n c(k) \lesssim n. $$

    这种“平均化”性质（Averaging condition）使得高频能量的累积不会导致控制失效。

  3. 通过重缩放控制超临界范数的全局小性： 在能量估计中，会出现项 $\|u\|_{X_1} \int_0^t \|\nabla u\|^2 d\tau$。如果 $\|u\|_{X_1}$ 无界，无法直接吸收误差。作者利用 Lemma 2.1 的性质：对于任意小的 $\epsilon > 0$，存在足够大的伸缩因子 $\lambda = 2^{2^{2^l}}$（$l$ 依赖于 $t$），使得 $\|\lambda u(t, \lambda \cdot)\|_{X_1} \leq \epsilon$。 这意味着，虽然在原始尺度下解的 $X_1$ 范数可能很大，但在经过适当重缩放的尺度下，它是小量。通过在每个时间点 $t$ 选择合适的缩放，作者确保了非线性项的影响被控制住，从而使得能量估计中的关键项可以被耗散项吸收，最终得到独立于 $\nu$ 的先验估计。

与之前工作的比较

本文在纳维-斯托克斯方程的正则性理论和欧拉方程的全局存在性方面做出了显著的改进，主要体现在初始数据正则性要求的放宽以及对超临界情形的处理上。

a) 正则性指数的改进

  1. 相对于经典 Leray 理论： Leray 最初证明了 $u_0 \in L^2(\mathbb{R}^3)$ 的全局弱解存在性，但未涉及 $L^2$ 以外的高阶正则性。 本文证明了对于 $d \geq 3$，只要初始数据 $u_0 \in H^s(\mathbb{R}^d)$ 且 $s \geq -1 + d/2$，解就是全局正则的。这实际上是亚临界到临界边界的正则性结果。

  2. 相对于临界空间理论 (Critical Spaces)： 过去二十年的许多工作（如 Escauriaza-Seregin-Šverák 对 $L^3$，Gallagher-Koch-Planchon 对 Besov 空间 $\dot{B}^{-1+3/p}_{p,q}$）集中在临界空间或稍高于临界的空间（Log-Lipschitz 空间）。 本文的结果不仅涵盖了这些临界情形，还推广到了更弱的超临界情形。虽然初始数据 $H^s$ ($s < -1+d/2$) 通常对应于超临界问题，但本文通过引入具有对数权重的空间 $X_1$，实际上在超临界框架下获得了类似于亚临界的控制。这比 Seregin (2020) 和 Pan (2019) 关于对数超临界准则（针对轴对称解）的工作更具一般性（适用于任意 $d \geq 3$ 且非轴对称）。

  3. 相对于近期超临界结果： 相比于 Barker 和 Prange (2020, 2021) 的工作，他们利用 Tao 的对数超临界估计给出了 3D 弱解的正则性准则（涉及 $\log \log \log$ 的积分条件）。本文提供了一个更统一的框架，不仅给出正则性准则，还直接给出了能量估计 (1.3)。特别重要的是，估计常数 $C(s)$ 独立于粘性系数 $\nu$。这是许多传统方法难以做到的，因为通常的能量估计在 $\nu \to 0$ 时会失效。

b) 欧拉方程全局可解性的推进

  1. 从 Navier-Stokes 到 Euler 的过渡： 欧拉方程是纳维-斯托克斯方程的极限情况（$\nu \to 0$）。传统上，获得欧拉方程的全局解非常困难，因为缺乏耗散机制，非线性项可能导致有限时间爆破。 本文利用纳维-斯托克斯方程估计的一致有界性（uniform bound with respect to $\nu$）。由于定理 1.1 中的先验估计 (1.3) 不依赖于 $\nu$，作者可以通过范数收敛论证（vanishing viscosity argument），从纳维-斯托克斯方程的强解序列过渡到欧拉方程的解。

  2. 初始数据范围的推广： 对于欧拉方程，经典结果通常要求 $s > 1 + d/2$（保证速度场 Lipschitz 连续，从而满足经典存在唯一性）。 本文证明了只要 $s \geq -1 + d/2$，欧拉方程就存在全局解（尽管唯一性仅在 $s > 1 + d/2$ 时保证）。这在 $s < 1 + d/2$ 的情形下是一个重要的全局存在性结果，覆盖了临界空间 $\dot{H}^{-1+d/2}$ 及其附近的弱正则性数据。

c) 方法论上的突破

之前的工作（如 Gallagher-Koch-Planchon）多依赖于 Profile Decomposition（剖面分解）和 Bona 方法，或者依赖于已知的 $L^3$ 估计。本文采用了完全不同的策略：

  • 内部能量测试：不再仅仅将非线性项 $(u \cdot \nabla)u$ 视为普通的二次非线性项，而是利用高频部分的分解和对流项的抵消性质 $((u \cdot \nabla)u_k, u_k) = 0$。
  • 稀疏加权空间：构造 $X_1$ 空间的技巧非常新颖。通过对数权重 $\log j$，仅在稀疏的频率带上降低范数要求，从而在整体上获得了一个足够“弱”以至于可以被控制的超临界范数，同时避免了求和系数爆炸的问题。

总结来说，这篇论文通过引入一种新型的稀疏对数加权空间，结合精细的频率分解和重缩放技术，在 $d \geq 3$ 的高维情形下，显著放宽了纳维-斯托克斯方程整体正则性和欧拉方程整体存在性对初始数据正则性的要求，解决了超临界情形下的一致估计难题。

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

1. 具有奇异快扩散和边界效应的Burgers方程向激波轮廓的收敛

Convergence to shock profiles for Burgers equation with singular fast-diffusion and boundary effect

作者: Xiaowen Li, Ming Mei

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15900v1

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

摘要翻译

中文标题: 具有奇异快扩散和边界效应的Burgers方程向激波轮廓的收敛

摘要翻译: 本文研究半空间 $(0,+\infty)$ 上具有边界条件 $u|_{x=0}=u_->0$ 和 $u|_{x=+\infty}=0$ 的 Burgers 方程 $u_t +f(u)_x = (\frac{u_{x}}{u^{1-m}})_x$ 的粘性激波轮廓的渐近稳定性。其中，参数 $\frac{1}{2}<m<1$ 度量快扩散的强度。由于在 $u=0$ 处扩散率 $\left(\frac{u_x}{u^{1-m}} \right)_x$ 的显著奇异性以及边界层的存在，我们面临一个关键挑战。我们证明了 $u$ 的长时间行为会收敛到一个平移后的激波轮廓 $U(x-st-d(t))$，其中 $d(t)$ 受 $x=0$ 处边界层动力学的支配，并由初始数据 $u(x,0)$ 驱动。为了克服由快扩散引起的奇异性以及边界层对波稳定性带来的不良影响，我们巧妙地引入了一些新的加权能量估计技术。

2. 时间分数阶非线性扩散方程 $L^1$ 解的存在唯一性

Existence and uniqueness of $L^1$-solutions to time-fractional nonlinear diffusion equations

作者: Mikiya Kametaka, Tatsuki Kawakami

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15618v1

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

摘要翻译

中文标题: 时间分数阶非线性扩散方程 $L^1$ 解的存在唯一性

摘要翻译: 我们建立了时间分数阶多孔介质型非线性扩散方程 Cauchy 问题 $L^1$ 解的全局存在唯一性。此外，我们给出了时间分数阶快速扩散方程 $L^1$ 解的质量守恒律，并证明了对于任意非负 $L^1$ 解，有限时间绝灭现象不会发生。

3. 线性密度相关 BGK 方程的非流体动力学解

Non-Hydrodynamic Solutions to the linear Density-dependent BGK equation

作者: Florian Kogelbauer

类别: math.AP, math-ph

发布日期: 2026-01-21

arXiv ID: 2601.15259v1

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

摘要翻译

中文标题: 线性密度相关 BGK 方程的非流体动力学解

摘要翻译: 本文证明了 $d$ 维线性密度相关 BGK 方程非流体动力学解的存在性。具体而言，我们证明了对于任意克努森数 $τ$，存在一种初始条件，使得宏观质量密度的耗散率发散且满足 $\sim 1/τ$ 的渐近行为。我们的结果依赖于对线性 BGK 算子的详尽谱分析，结合傅里叶级数与拉普拉斯变换得到时间相关问题的显式解公式，并后续利用复分析中的围道积分方法进行论证。

不相关论文（基本信息）

1. 关于拟线性图斑形成反应-扩散-对流系统的Ginzburg-Landau逼近研究

On the Ginzburg-Landau approximation for quasilinear pattern forming reaction-diffusion-advection systems

作者: Théo Belin, Guido Schneider

类别: math.AP, math.DS

发布日期: 2026-01-22

arXiv ID: 2601.16145v1

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

摘要: We prove that the Ginzburg-Landau equation correctly predicts the dynamics of quasilinear pattern-forming reaction-diffusion-advection systems, close to the first instability. We present a simple theorem which is easily applicable for such systems and relies on key maximal regularity results. The theorem is applied to the Gray-Scott-Klausmeier vegetation-water interaction model and its application to general reaction-diffusion-advection systems is discussed.

2. 点云上超图 $p$-拉普拉斯方程的连续极限

Continuum limit of hypergraph $p$-Laplacian equations on point clouds

作者: Kehan Shi

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.16063v1

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

摘要: This paper studies a class of $p$-Laplacian equations on point clouds that arise from hypergraph learning in a semi-supervised setting. Under the assumption that the point clouds consist of independent random samples drawn from a bounded domain $Ω\subset\mathbb{R}^d$, we investigate the asymptotic behavior of the solutions as the number of data points tends to infinity, with the number of labeled points remains fixed. We show, for any $p>d$ in the viscosity solution framework, that the continuum limit is a weighted $p$-Laplacian equation subject to mixed Dirichlet and Neumann boundary conditions. The result provides a new discretization of the $p$-Laplacian on point clouds.

3. 量子欧几里得空间上的藤田指数

Fujita exponents on quantum Euclidean spaces

作者: Edward McDonald, Michael Ruzhansky, Serikbol Shaimardan, Kanat Tulenov

类别: math.AP, math.OA

发布日期: 2026-01-22

arXiv ID: 2601.16053v1

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

摘要: We study the well-posedness of a non-linear heat equation with power nonlinearity with positive initial data on quantum Euclidean spaces. We prove a noncommutative analogue of the classical Fujita theorem by identifying the critical exponent separating finite-time blow-up from global existence for small initial data. Moreover, we establish a fundamental inequality in general semifinite von Neumann algebras that is of independent interest and plays a crucial role in the study of global existence and local well-posedness of solutions of nonlinear equations in noncommutative setting.

4. 从散射数据中恢复满足零条件的半线性位势

The Recovery of Semilinear Potentials Satisfying Null Conditions From Scattering Data

作者: Joel Nathe, Antônio Sá Barreto

类别: math.AP, math-ph

发布日期: 2026-01-22

arXiv ID: 2601.15997v1

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

摘要: We construct oscillatory solutions of fully semilinear wave equations in Minkowski space satisfying a null condition of the form

$$ \square u:=(-\partial_{x_0}^2 +\sum_{j=1}^n \partial_{x_j}^2 )u= q(x,u)((\partial_{x_0}u)^2-|\nabla_{x'}u|^2), $$

$$ x=(x_0,x'), \;\ x'=(x_1,\ldots, x_n) \text{ and } x_0=t \text{ is the time variable,} $$

on an interval $x_0\in [-T,T]$, $T<\infty$ arbitrary, which consist of the superposition of a non-oscillatory background solution and a single phase train of highly oscillatory waves of wave length $h\ll1$ and amplitudes given by powers of $h$; the waves interact with the nonlinearity and we measure the response $u(x_0,x')|_{x_0=T'}$ at a fixed time $x_0=T'<T$. We show that the coefficient of amplitude $h$ of the oscillatory part of the nonlinear geometric optics expansion of the solution determines the light-ray transform of a vector field associated with $q(x,u)$, which determines $q(x,u)$ uniquely in the maximal region determined by the data. Our methods also work for systems of semilinear wave equations satisfying null conditions, but in this paper we focus on the scalar case.

5. 超临界维数中趋化方程组解的整体存在性、有限时间爆破的尖锐判据与完全分类

A sharp criterion and complete classification of global-in-time solutions and finite time blow-up of solutions to a chemotaxis system in supercritical dimensions

作者: Yuri Soga

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15990v1

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

摘要: We consider the chemotaxis system with indirect signal production in the whole space, \begin{equation}\label{abst:p}\tag{$\star$} \begin{cases} u_t = Δu - \nabla \cdot (u\nabla v),\ 0 = Δv + w,\ w_t = Δw + u \end{cases} \end{equation} with emphasis on supercritical dimensions. In contrast to the classical parabolic-elliptic Keller–Segel system, where the analysis can be reduced to a single equation, the above system is essentially parabolic-parabolic and does not admit such a reduction. In this paper, we establish a sharp threshold phenomenon separating global-in-time existence from finite time blow-up in terms of scaling-critical Morrey norms of the initial data. In particular, we prove the existence of singular stationary solutions and show that their Morrey norm values serve as the critical thresholds determining the long-time behavior of solutions. Consequently, we identify new critical exponents at which the long-time behavior of solutions changes. This yields a complete classification of the long-time behavior of solutions, providing the first such results for the essentially parabolic-parabolic chemotaxis system \eqref{abst:p} in supercritical dimensions.

6. 非线性材料参数在拟线性Lamé系统中的恢复方法研究

Recovery of nonlinear material parameters in a quasilinear Lamé system

作者: David Johansson, Yavar Kian

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15881v1

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

摘要: We investigate the inverse problem of determining nonlinear elastic material parameters from boundary stress measurements corresponding to prescribed boundary displacements. The material law is described by a nonlinear, space-independent elastic tensor depending on both the displacement and the strain, and gives rise to a general class of quasilinear Lamé systems. We prove the unique and stable recovery of a wide class of space-independent nonlinear elastic tensors, including the identification of two nonlinear isotropic Lamé moduli as well as certain anisotropic tensors. The boundary measurements are assumed to be available at a finite number of boundary points and, in the isotropic case, at a single point. Moreover, the measurements are generated by boundary displacements belonging to an explicit class of affine functions. The analysis is based on structural properties of nonlinear Lamé systems, including asymptotic expansions of the boundary stress and tensorial calculus.

7. 具无界电势的磁 Schrödinger 算子的定量 Borg-Levinson 定理

Quantitative Borg-Levinson theorem for the magnetic Schödinger operator with unbounded electrical potential

作者: Mourad Choulli, Hiroshi Takase

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15847v1

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

摘要: The first author established in [8] a quantitative Borg-Levinson theorem for the Schrödinger operator with unbounded potential. In the present work, we extend the results in [8] to the magnetic Schrödinger operator. We discuss both the isotropic and anisotropic cases. We establish Hölder stability inequalities of determining the electrical potential or magnetic field from the corresponding boundary spectral data.

8. 具有阶跃状振荡背景的非局部修正Korteweg-de Vries方程的Riemann-Hilbert方法

Riemann-Hilbert approach for the nonlocal modified Korteweg-de Vries equation with a step-like oscillating background

作者: Yan Rybalko

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15841v1

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

摘要: This work focuses on the Cauchy problem for the nonlocal modified Korteweg-de Vries equation

$$ u_t(x,t)+6u(x,t)u(-x,-t)u_x(x,t)+u_{xxx}(x,t)=0, $$

with the oscillating step-like boundary conditions: $u(x,t)\to 0$ as $x\to-\infty$ and $u(x,t)\backsimeq A\cos(2Bx+8B^3t)$ as $x\to\infty$, where $A,B>0$ are arbitrary constants. The main goal is to develop the Riemann-Hilbert formalism for this problem, paying a particular attention to the case of the ``pure oscillating step’’ initial data, that is $u(x,0)=0$ for $x<0$ and $u(x,0)=A\cos(2Bx)$ for $x\geq0$. Also, we derive three new families of two-soliton solutions, which correspond to the values of $A$ and $B$ satisfying $B<\frac{A}{4}$, $B>\frac{A}{4}$, and $B=\frac{A}{4}$.

9. 具有自由边界的各向异性极小图的半空间Liouville定理

A half-space Liouville theorem for anisotropic minimal graph with free boundary

作者: Guofang Wang, Wei Wei, Chao Xia, Xuwen Zhang

类别: math.DG, math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15788v1

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

摘要: In this paper we prove the following Liouville-type theorem: any anisotropic minimal graph with free boundary in the half-space must be flat, provided that the graph function has at most one-sided linear growth. This extends the classical results of Bombieri-De Giorgi-Miranda and Simon to an appropriate free boundary setting.

10. 双线性傅里叶积分算子的局部光滑性估计

Local smoothing estimates for bilinear Fourier integral operators

作者: Duván Cardona

类别: math.AP, math.FA

发布日期: 2026-01-22

arXiv ID: 2601.15667v1

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

摘要: We formulate a local smoothing conjecture for bilinear Fourier integral operators in every dimension $d \ge 2,$ derived from the celebrated linear case due to Sogge, which we refer to as the \emph{bilinear smoothing conjecture}. We show that the linear local smoothing conjecture implies this bilinear version. As a consequence of our approach and due to the recent progress on the subject, we establish local smoothing estimates for Fourier integral operators in dimension $d=2,$ that is, on $\mathbb{R}^2_x \times \mathbb{R}_t$. Also, a partial progress is presented for the high-dimensional case $d\geq 3.$ In particular, our method allows us to deduce that the bilinear local smoothing conjecture holds for all odd dimensions $d$.

11. 通过局部分岔结构研究肿瘤生长自由边界问题中首次分岔的线性稳定性

Linear stability of the first bifurcation in a tumor growth free boundary problem via local bifurcation structure

作者: Junying Chen, Ruixiang Xing

类别: math.AP

发布日期: 2026-01-22

arXiv ID: 2601.15647v1

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

摘要: In this paper, we consider a 3-dimensional free boundary problem modeling tumor growth with the Robin boundary condition. The system involves a positive parameter $μ$ which reflects the intensity of tumor aggressiveness. Huang, Zhang and Hu [Nonlinear Anal. Real World Appl. 2017(35), 483-502] have shown that for each $μ_n$ ($n$ even) in a strictly increasing sequence $\{ μ_n \}(n\geq 2)$, there exists a stationary bifurcation solution $(σ_n(\varepsilon),p_n(\varepsilon),r_n(\varepsilon))$ with $μ= μ_n(\varepsilon)$ bifurcating from $μ_n$. We first derive that the bifurcation curve $(r_2(\varepsilon),μ_2(\varepsilon))$ exhibits a transcritical bifurcation with $μ_2'(0)<0$. Moreover, we show that the stationary bifurcation solution $(σ_2(\varepsilon),p_2(\varepsilon),r_2(\varepsilon))$ is linearly unstable for small $|\varepsilon|$ under non-radially symmetric perturbations. In contrast to the linear stability of the radially symmetric stationary solution, the lack of explicit expressions for bifurcation solutions adds great difficulty in analyzing their linear stability. The novelty of this paper lies in the use of the bifurcation curve’s structure to overcome the above difficulties. Moreover, this linear stability result is not established using the standard method, due to an eight-dimensional generalized kernel at eigenvalue 0 for the linearized operator.

12. 椭圆、抛物和动力方程的De Giorgi正则性理论

De Giorgi’s regularity theory for elliptic, parabolic and kinetic equations

作者: Cyril Imbert

类别: math.AP

发布日期: 2026-01-21

arXiv ID: 2601.15238v1

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

摘要: This book presents a comprehensive regularity theory for solutions of elliptic, parabolic, and kinetic equations. The foundation of this theory was laid by E. De Giorgi’s groundbreaking resolution of Hilbert’s nineteenth problem in 1956. The innovative tools he developed to tackle this problem proved to be remarkably versatile. In 1957, just one year later, J. Nash independently developed analogous techniques for parabolic equations, concurrently with De Giorgi’s research. By the year 2000, these techniques had been extended to address elliptic and parabolic equations featuring integral diffusion, such as the fractional Laplacian. More recently, the theory has evolved to encompass kinetic equations, accommodating both local and integral diffusions. This book aims to present these results in a unified and coherent manner, beginning with the classical elliptic framework and progressing through to the most recent advancements in kinetic equations.