生成时间: 2026年08月27日 13:08:19

运行状态: 完成

计划处理论文数量: 29

已完成论文数量: 29

未完成论文数量: 0

重点关注论文数量: 8

了解领域论文数量: 5

不相关论文数量: 16

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

1. 涡丝动力学与涡环运动的再审视

Vortex filament dynamics and vortex ring motion revisited

作者: Andrew D Gilbert

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

发布日期: 2026-08-26

arXiv ID: 2608.25827v1

摘要翻译: 理想不可压缩流体中细长涡的运动是水力学中的一个经典问题。十九世纪，Kelvin、Helmholtz、Hicks 和 Dyson 的工作已经给出了涡环速度的公式；近年来，人们又建立了若干用于模拟一般形状细长涡度管运动的模型。对于局部化在曲线 $$C(t)$$ 附近的涡管或涡丝而言，衡量其细长程度的相关小参数是管半径与曲线 $$C$$ 的曲率半径之比。本文通过引入一套与细长涡中涡面几何密切相关的坐标系，重新审视了这些经典问题。在该坐标系下，流体质点的运动具有作用-角（action–angle）形式，从而使涡度方程得以极大简化。然而与此同时，涡的形状或演化的大部分特征都被转移到了坐标系的描述中，特别是相应的度量与体积形式。

Comment: 无

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

1. 研究对象和背景

理想不可压缩流体中细长涡的运动是流体力学中的经典问题。十九世纪，Kelvin, Helmholtz, Hicks 和 Dyson 等人给出了涡环速度的公式；近年来，针对一般形状细长涡管（vortex tube）或涡丝（vortex filament）的动力学模型也不断涌现。对于一个局域在曲线 $$C(t)$$ 附近的涡管或涡丝，衡量其细长性的小参数为涡管半径 $$a$$ 除以 $$C$$ 的曲率半径，即 $\epsilon = a\kappa \ll 1$ 。

本文旨在重新审视这一系列经典问题，通过引入一个与细长涡中涡面几何紧密相关的坐标系，来简化并重构涡动力学的数学框架。在此坐标系下，流体元素的运动具有作用-角（action-angle）形式，从而使得涡量方程发生根本性的简化。然而，这种简化的代价是涡的形状及其演化的大部分信息被转移到了坐标系的描述中，特别是相应的度量（metric）和体积形式（volume form）上。

由于该坐标系是非正交且随时间变化的，微分几何工具成为描述涡量结构和坐标系结构的最自然选择。作者建立了一个通用的数学框架，并在仅假设涡保持细长性的前提下，尽可能推导出允许涡核和中心曲线 $$C(t)$$ 进行任意运动和变形的方程组。随后，作者将该模型应用于计算具有任意轴向流的细长涡环的运动，直接求解受扰涡面的形状，并与早期研究结果进行了对比。该通用框架为未来开发涡运动和相互作用的简化方程奠定了基础。

2. 主要定理或主要结果

论文的主要成果在于构建了一套基于微分几何的、适用于任意细长涡演化的通用数学框架，并将其成功应用于细长涡环运动的精确计算。

通用数学框架的建立：
- 作者引入了三个关键的坐标系：环境坐标系（Ambient coordinate system）、涡坐标系（Vortex coordinate system）和拉格朗日坐标系（Lagrangian coordinate system）。其中，拉格朗日坐标系 $(r, \theta, z, t)$ 是与涡线螺旋所在的表面最紧密适配的坐标系，涡面由常数坐标 $$r$$ 的表面给出。
- 利用微分几何语言（外微分、内积、李导数等），作者将理想不可压缩流体的涡量方程写成纯几何形式： $\partial_t \xi + d(v \lrcorner \xi) = 0, \quad \xi = d\nu, \quad \nu = u^\flat, \quad d(u \lrcorner \mu) = 0, \quad u = U + v.$ 这里 $\xi$ 是涡量二形式， $\nu$ 是动量一形式， $$u$$ 是全速度场， $$U$$ 是背景流速度， $$v$$ 是相对速度场， $\mu$ 是体积形式。
- 在细长极限（slender vortex limit）下，作者推导出了度量和背景流性质的简化形式，从而获得了一组仅依赖于细长性假设的通用方程系统（详见附录 C）。这些方程允许对涡核内部复杂的涡量分布和形状变形进行描述，而无需预先假设涡量剖面。
细长涡环运动的计算：
- 作为框架的应用示例，作者计算了具有任意轴向流的细长涡环的运动。
- 通过直接求解控制涡量和动量的两个微分方程，得到了扰动涡面的几何形状。
- 通过与标准的 Biot-Savart 定律计算相匹配，确定了涡环的传播速度。
- 结果与 Saffman (1970), Widnall, Bliss & Zalay (1971) 等早期研究相符，但新之处在于能够直接解析地求解涡面的形状，而不仅仅是得到速度表达式。

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

本文的核心方法是将微分几何工具系统地应用于涡动力学，以处理非正交、时变坐标系下的复杂几何问题。

微分几何框架：
- 作者采用微分形式的语言来表述流体动力学方程，这使得方程在不同坐标系之间的转换变得自然且紧凑。特别是，利用 Cartan 公式将涡量输运项表示为李导数 $L_v \xi$ ，清晰地揭示了其几何意义。
- 通过引入 Frenet-Serret 标架 $${n, b, t}$$ 来描述中心曲线 $$C(t)$$ 的局部几何，并建立了曲线动力学方程，包括曲率 $\kappa$ 、扭转 $\tau$ 和雅可比 $$J$$ 的时间演化。
混合欧拉-拉格朗日坐标系：
- 作者设计了一种“混合”坐标系，其中背景流 $$U$$ 吸收了涡丝中心线及其邻域的运动，而相对流 $$v$$ 则捕捉流体元素绕中心线的旋转。这种分解允许更灵活地处理涡丝的平移、拉伸和旋转。
- 在拉格朗日坐标系中，涡面被定义为常数 $$r$$ 的表面，这极大地简化了涡量方程的结构，因为涡线在这些表面上螺旋缠绕，使得涡量的输运主要体现为角度坐标 $\theta$ 的变化。
细长极限渐近展开：
- 在细长极限 $\epsilon = a\kappa \ll 1$ 下，作者对度量和背景流性质进行了系统展开，分离出主导阶和高阶修正项。这种方法避免了传统匹配渐近展开中对涡核内部结构的强假设，而是直接从几何出发推导内部涡量分布。
Biot-Savart 匹配：
- 对于涡环运动，作者将内部流动的解与外部无旋流动通过 Biot-Savart 定律进行匹配，从而确定涡环的整体运动速度。这一过程确保了内外解的一致性，并提供了验证模型正确性的手段。

4. 与之前工作的比较

与现有文献相比，本文的主要推进体现在以下几个方面：

方法论的创新：
- 传统的涡丝模型（如 Moore & Saffman, 1972; Marshall, 1991-1993）通常假设涡管横截面为圆形，并预设涡量剖面（如均匀轴向涡量），从而抑制了核心动力学的复杂性。相比之下，本文的方法不假设核心结构，而是通过系统展开来获得内部涡量剖面和曲线运动的方程。这与 Ting 及其合作者（Tung & Ting, 1967; Ting, 1971; Callegari & Ting, 1978; Klein & Ting, 1992）的方法更为接近，但本文采用了更强大的微分几何工具来处理非正交和时变坐标系，使得模型能够容纳更一般的涡核变形和运动。
对强相互作用和复杂变形的适应性：
- 现有模型在处理涡丝强烈相互作用（如发夹形偶极子形成、涡量剥离、有限时间奇异性等）时往往遇到困难。本文的框架仅要求涡始终保持细长，因此理论上可以应用于这些极端情况，尽管具体的应用将是未来工作的重点。
涡环运动的直接求解：
- 在涡环运动方面，早期工作（Kelvin, Helmholtz, Dyson, Hicks, Saffman, Widnall 等）主要关注传播速度的计算，尤其是包含对数项 $\log(R/a)$ 的结果。Saffman (1970) 和 Widnall, Bliss & Zalay (1971) 研究了轴向流的影响，发现轴向流会降低涡环速度甚至导致方向反转。Fraenkel (1970, 1972) 进行了高阶计算，Kopev & Chernyshev (1991, 1995) 研究了涡环对等涡扰动的稳定性。本文的新颖之处在于，它不仅重现了已知的速度结果，还直接求解了受扰涡面的形状，提供了更完整的涡结构图像。
与面积坐标的联系：
- 本文的坐标系与 Childress (1987) 和 Childress & Gilbert (2025) 中使用的面积坐标密切相关，后者在对轴向速度和涡量进行轮廓平均后得到。这表明本文的工作是该系列研究的一部分，旨在揭示涡动力学的几何基础。
局限性与未来方向：
- 尽管本文建立的框架非常通用，但要将其应用于具体的新型涡配置或动力学过程，仍需针对特定几何和过程进行进一步的建模和近似。这是未来研究的重点。此外，目前的应用仅限于无粘情形，粘性效应的纳入将是后续工作需要考虑的问题。

2. 半空间中定常 Navier–Stokes 方程的强迫自相似解

Forced self-similar solutions to the stationary Navier–Stokes equations in a half-space

作者: Yun Wang, Chunjing Xie, Shaoheng Zhang

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25432v1

摘要翻译: 我们研究在半空间中、带无滑移边界条件的轴对称定常 Navier–Stokes 方程的自相似解，其驱动力为一个轴对称的 $$(-3)$$ -齐次外力。若该外力在单位球面上的切向旋度足够小，则证明存在唯一的小解；当外力无旋时，解自动无旋且唯一。对于无旋外力 $\boldsymbol{F}$ ，我们引入一个缩放参数 $$λ$$ ，并考虑带有外力 $λ\boldsymbol{F}$ 的系统；我们证明解恰好存在于包含零的一个开区间内的 $$λ$$ 值上。同样的方法可推广到带无滑移边界条件的实心锥体情形，其中更小的张角允许在更大的外力下存在解。

Comment: 23 pages; comments are welcome

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

1. 研究对象和背景

本文研究三维上半空间 $\mathbb{R}^3_+ = \{x=(x_1,x_2,x_3) \in \mathbb{R}^3: x_3 > 0\}$ 中带有无滑移边界条件的稳态不可压缩 Navier-Stokes 方程:

\begin{cases} -\Delta u + (u \cdot \nabla)u + \nabla p = F & \text{在 } \mathbb{R}^3_+, \\ \operatorname{div} u = 0 & \text{在 } \mathbb{R}^3_+, \\ u = 0 & \text{在 } \partial \mathbb{R}^3_+ \setminus \{0\}. \end{cases}

其中 $$u$$ 和 $$p$$ 分别为速度场和压力， $$F$$ 为外部力。该方程组在标度变换 $u_\lambda(x) = \lambda u(\lambda x)$ , $p_\lambda(x) = \lambda^2 p(\lambda x)$ , $F_\lambda(x) = \lambda^3 F(\lambda x)$ 下保持不变。自相似解是指满足 $u(\lambda x) = \lambda^{-1}u(x)$ 且 $p(\lambda x) = \lambda^{-2}p(x)$ 的解，相应地，保持标度不变的外部力必须满足齐次性条件 $F(\lambda x) = \lambda^{-3}F(x)$ 。

自相似解在理解流体流动的局部与整体结构中起着基础性作用，它们自然地出现在射流、尾迹以及奇点附近或远距离处的涡旋运动的领头阶渐近描述中。此外，自相似解为各种区域中解的正则性理论、Liouville 型性质以及渐近行为提供了关键见解。在无边界（全空间）情形下，经典的 Landau 解构成了一个单参数的轴对称、无旋涡（swirl-free）、 $$(-1)$$ -齐次解族。然而，当存在边界时，问题发生了实质性变化。Kang, Miura, and Tsai 证明，在半空间中，任何满足无滑移或 Navier 滑移边界条件的轴对称 $$(-1)$$ -齐次解恒为零。这意味着，与全空间的 Landau 解族不同，半空间中非平凡的轴对称齐次流动必须由非势型（non-potential）外部力来维持。

因此，一个自然的问题是：给定一个 $$(-3)$$ -齐次的外部力 $$F$$ ，半空间是否存在相应的自相似解？本文在假设 $$F$$ 是轴对称的条件下，研究轴对称自相似解的存在性与唯一性。主要关注点在于外部力的切向旋度（tangential curl）对解的影响，特别是在无旋涡（swirl-free）情形下解对力幅值的依赖关系。

2. 主要定理或主要结果

本文的主要结果包括以下几个方面：

小切向旋度下的存在性与唯一性 (Theorem 1.1): 设 $$F$$ 是 $\mathbb{R}^3_+$ 上的 $$C^1$$ 类轴对称向量场，且关于尺度为 $$-3$$ 齐次。令 $(\operatorname{curl} F)_{tan}$ 表示 $\operatorname{curl} F$ 在单位球面 $$S^2$$ 上的切向分量。若
$\|(\operatorname{curl} F)_{tan}\|_{C^0(S^2 \cap \mathbb{R}^3_+)} \le \varepsilon,$
其中 $\varepsilon > 0$ 足够小，则系统 (1) 在有界角剖面类中存在唯一的轴对称自相似小解。具体地，对于任意 $\alpha \in (0,1)$ ，有
$u \in C^{2,\alpha}_{loc}(\mathbb{R}^3_+ \setminus \{0\}), \quad p \in C^{1,\alpha}_{loc}(\mathbb{R}^3_+ \setminus \{0\}).$
此外，如果 $$F$$ 是无旋涡的（即没有绕轴旋转的分量），则任何轴对称自相似解必然是无旋涡的且唯一的。
无旋涡外力下解对幅值参数 $\lambda$ 的依赖关系 (Theorem 1.2): 设 $$F$$ 是 $\mathbb{R}^3_+$ 上的 $$C^1$$ 类轴对称、无旋涡外部力，且关于尺度为 $$-3$$ 齐次。引入参数 $\lambda$ ，考虑外部力为 $\lambda F$ 的 Navier-Stokes 方程。存在常数 $\underline{\lambda} < 0 < \bar{\lambda}$ （可能为 $\pm \infty$ ），使得对于每个 $\lambda \in (\underline{\lambda}, \bar{\lambda})$ ，方程存在唯一的轴对称自相似解；而对于 $\lambda \notin (\underline{\lambda}, \bar{\lambda})$ ，不存在这样的有界剖面解。在球坐标下写 $F = \rho^{-3}(X(\varphi)e_\rho + Y(\varphi)e_\varphi)$ 。端点 $\underline{\lambda}$ 和 $\bar{\lambda}$ 的有限性由判别式
$$ N(\varphi) = \int_0^\varphi (G(\tau) + 2C) \sin \tau (\cos \varphi - \cos \tau) d\tau, \quad \varphi \in [0, \pi/2] $$
的符号决定，其中 $G(\varphi) = X(\varphi) + 2\int_0^\varphi Y(\tau)d\tau$ 且 $C = -\int_0^{\pi/2} G(\tau) \sin \tau \cos \tau d\tau$ 。具体而言：
- 若 $N \ge 0$ 且 $N \not\equiv 0$ ，则 $\underline{\lambda} \in (-\infty, 0)$ 且 $\bar{\lambda} = +\infty$ ；
- 若 $N \le 0$ 且 $N \not\equiv 0$ ，则 $\underline{\lambda} = -\infty$ 且 $\bar{\lambda} \in (0, +\infty)$ ；
- 若 $$N$$ 变号，则 $\underline{\lambda} \in (-\infty, 0)$ 且 $\bar{\lambda} \in (0, +\infty)$ ；
- 若 $N \equiv 0$ ，则 $\underline{\lambda} = -\infty$ 且 $\bar{\lambda} = +\infty$ 。
锥体情形的推广: 同样的方法可以推广到带有无滑移边界条件的实心锥体，其中较窄的开口角允许在较大的外部力下存在解。

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

本文的核心技术障碍在于处理带边界条件下的非线性耦合，特别是如何将偏微分方程系统简化为可分析的一维常微分方程，并精确刻画解的存在区间。

无旋涡情形的化简: 对于无旋涡（swirl-free）情况，作者首先证明了速度场的旋涡分量 $h \equiv 0$ 。随后，通过消去压力并利用散度自由条件，将系统转化为关于经线速度分量 $g(\varphi)$ 的二阶标量方程。进一步通过变量代换 $t = \cos \varphi$ 和 $L(t) = g(\varphi)\sin \varphi$ ，该方程被转化为 Riccati 初值问题：
$\frac{dM}{dt} + \frac{1}{2}M(t)^2 = \lambda N(t), \quad M(0) = 0.$
这里 $$M(t) = L(t)/(1-t^2)$$ 。解的存在性等价于该 Riccati 方程解在 $$t=1$$ 处的延拓性。利用比较原理，作者证明了可容许的幅值集合是一个开区间，而判别式的符号决定了端点的有限性。
Riccati 方程的可解性引理 (Lemma 3.1): 这是证明 Theorem 1.1 的关键工具。作者建立了一个关于一般区间 $$[t_0, 1]$$ 上 Riccati 方程初始值问题的可解性引理。若 $\|N\|_{C^0} < 1/(2\ell^2)$ （其中 $\ell = 1-t_0$ ），则存在唯一解 $M \in C^1[t_0, 1]$ ，并给出了具体的范数估计：
$\|M\|_{C^0} \le 2\ell \|N\|_{C^0}, \quad \|M\|_{C^1} \le 2(\ell+1)\|N\|_{C^0}.$
此外，还证明了两个不同右端项对应的解之间的 Lipschitz 连续性。这一引理不仅用于半空间情形，也为后续锥体情形的推广奠定了基础。
含旋涡情形的压缩映射论证: 当存在旋涡分量时，方位角方程先被求解以用 $$g$$ 和力的旋涡分量表示 $$h$$ 。剩余的方程定义了关于 $$g$$ 的非线性映射。作者证明，当切向旋度足够小时，该映射在有界角剖面的球上是压缩映射，从而保证了唯一小解的存在性。
正则性提升: 一旦获得了解的基本结构，作者利用标准的局部 Stokes 正则性理论，证明了角剖面的 $$C^1$$ 正则性足以保证整体解在 $\mathbb{R}^3_+ \setminus \{0\}$ 上具有 $C^{2,\alpha}$ 的速度和 $C^{1,\alpha}$ 的压力正则性。

4. 与之前工作的比较

本文与既有工作的主要区别和推进体现在以下几个方面：

边界条件的引入与非平凡解的存在性: 与全空间中的 Landau 解不同，半空间中的无外力轴对称齐次解恒为零（由 Kang, Miura, and Tsai 证明）。因此，半空间中的非平凡自相似解必须由外部力驱动。本文首次系统地研究了半空间中受迫自相似解的存在性问题，填补了这一空白。
外部力兼容性条件的差异: 在全空间中，Shi 的研究表明，轴对称 $$(-3)$$ -齐次外力的存在性需要满足特定的兼容性条件（如 Remark 1.1 中所述）。然而，在半空间或实心锥体中，由于边界条件的约束，不需要这种全局的兼容性条件。本文的结果表明，只要切向旋度足够小，解就存在，这放宽了对外力形式的全局限制。
切向旋度与全旋度的区分: 作者明确区分了切向旋度 $(\operatorname{curl} F)_{tan}$ 和全旋度 $\operatorname{curl} F$ 。通过 Example 1.1，作者构造了一族光滑轴对称力，其切向旋度的 $$C^0$$ 范数趋于零，但全旋度的 $$C^0$$ 范数发散。这表明 Theorem 1.1 适用于全旋度可能很大的力，而不仅仅是小全旋度的力。这是一个重要的技术细节，因为它意味着即使力的“强度”（由全旋度衡量）很大，只要其切向部分很小，解仍然存在。
无旋涡情形下解对幅值的精细刻画: 之前的工作通常只给出小力下的存在性。本文通过 Theorem 1.2 进一步刻画了无旋涡情形下解对力幅值参数 $\lambda$ 的依赖关系。作者证明了存在一个开区间 $(\underline{\lambda}, \bar{\lambda})$ ，使得解恰好在此区间内存在，并且端点的有限性由一个显式的判别式 $N(\varphi)$ 的符号决定。这种精细的刻画在之前的文献中尚未见到。
推广到锥体: 本文将方法推广到实心锥体，并分析了开口角对解存在性的影响。结果显示，较窄的开口角允许在更大的外部力下存在解。这为理解几何形状对流体解的影响提供了新的视角。

3. 二维欧拉方程无穷通道中定常解的灵活性与刚性

Flexibility and rigidity of steady states of the two-dimensional Euler equations in an infinite channel

作者: Yupei Huang, Chunjing Xie, Chilin Zhang

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25859v1

摘要翻译: 我们研究二维不可压欧拉方程在无穷通道中的定常解，其远场极限为均匀非滞止剪切流。在光滑范畴中，对于一大类给定的远场剪切剖面，通过极小极大方法构造流函数的半线性椭圆方程的二维解，存在非剪切的定常解。在解析范畴中，我们建立了针对解析定常解的比较原理，并证明对于稠密的一族解析均匀非滞止剪切剖面，具有给定远场的每个解析定常解本身必须是一个剪切流。特别地，存在一些远场剪切剖面在光滑范畴中表现出灵活性，而在解析范畴中表现出刚性。此外，这种稠密刚性是锐利的，因为存在一个解析剪切剖面，它在解析范畴中允许灵活性。

Comment: 无

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

1. 研究对象和背景

本文研究二维不可压缩欧拉方程在无限通道 $\Omega = \mathbb{R} \times (-1, 1)$ 中的稳态解。具体而言，我们关注那些在远场极限下收敛于均匀非停滞剪切流（uniformly non-stagnant shear flows）的稳态。二维不可压缩欧拉方程由以下系统给出：

\begin{cases} \partial_t v + v \cdot \nabla v + \nabla P = 0, \\ \nabla \cdot v = 0. \end{cases}

其中 $$v$$ 为速度场， $$P$$ 为压力。引入流函数 $\psi$ 使得 $v = \nabla^\perp \psi$ ，则稳态方程可以改写为

\nabla^\perp \psi \cdot \nabla \Delta \psi = 0.

问题设定要求流函数 $\psi$ 满足边界条件 $\psi(x, \pm 1) = \varphi(\pm 1)$ ，且在 $|x| \to \infty$ 时， $\psi(x, y)$ 以 $$C^2$$ 范数收敛于给定的剪切流剖面 $\varphi(y)$ 。

该问题的核心在于回答一个自然的问题：如果一个稳态在远场收敛于均匀非停滞剪切流，它是否必须本身就是一个剪切流？即，是否存在非剪切（non-shear）的稳态结构？此前的研究表明，在扰动区域内，背景剪切流会将非剪切分量输运至空间无穷远，导致解弱收敛于背景剪切流。然而，在非扰动区域，这种松弛机制是否依然成立尚不明确。此外，已有工作如 Hamel 和 Nadirashvili 证明了如果速度模长远离零，则稳态必须是剪切流；而后续工作考虑了存在停滞点的情况，但在层流假设下仍得到刚性结论。

然而，这些结果主要依赖于特定的几何或动力学约束。对于一般的远场剪切剖面，特别是当剖面具有更复杂的曲率性质时，稳态的分类变得复杂。本文试图填补这一空白，通过区分光滑类别（smooth category）和分析类别（analytic category），探讨远场剪切剖面对稳态结构的决定性影响。具体地，我们将证明在光滑类别下，对于广泛的一类远场剪切剖面，存在非剪切稳态；而在分析类别下，对于稠密的一类分析剪切剖面，所有分析稳态必须本身就是剪切流。这种“柔性”与“刚性”的对比揭示了正则性在流体动力学稳态分类中的关键作用。

2. 主要定理或主要结果

本文的主要结果概括如下：

定理 1.1 (柔性与刚性) 设 $\varphi \in C^\infty[-1, 1]$ 是一个单调递增函数，且满足以下条件：

\varphi'(y) > 0, \quad \varphi'''(y) > 0, \quad \text{以及} \quad \varphi'(y)\varphi''''(y) < \varphi''(y)\varphi'''(y), \quad \forall y \in [-1, 1].

(i) 光滑柔性: 在上述条件下，存在一个非剪切的光滑稳态解 $\psi$ 满足方程 (1.2)，即在通道内 $\nabla^\perp \psi \cdot \nabla \Delta \psi = 0$ ，边界条件 $\psi(x, \pm 1) = \varphi(\pm 1)$ ，且远场行为 $\lim_{|x| \to \infty} \|\psi(x, \cdot) - \varphi(\cdot)\|_{C^2([-1,1])} = 0$ 。

(ii) 分析刚性: 存在一族分析函数 $\varphi$ ，它们在拓扑 $C^\omega([-1, 1]) \cap \{\varphi | \varphi' > 0\}$ 中是稠密的，使得对于这族函数中的任意一个，任何满足上述远场条件的分析稳态解 $\psi$ 必然是一个剪切流（即 $\psi(x, y) = \varphi(y)$ ）。

推论/备注:

锐利性: 上述分析刚性结果是锐利的。存在一个没有停滞点的分析剪切剖面，其对应的分析稳态中存在非剪切解（见命题 2.2）。这意味着并非所有分析剪切剖面都导致刚性，刚性仅在稠密集合上成立。
具体例子: 考虑 $\Phi_0(t) = \frac{t^2}{2} + t + e^{-t} - 1$ ，并令 $\varphi$ 满足 $\varphi'(y) = \sqrt{2\Phi_0(\varphi(y))}$ 且 $\varphi(-1)=1$ 。可以验证该 $\varphi$ 满足定理 1.1 的条件 (1.3)。在此情况下，光滑类别下存在非剪切稳态，但分析类别下唯一解为剪切流 $\varphi$ 。

命题 2.2 (分析柔性) 存在一个分析函数 $\varphi(y)$ 满足 $\varphi' > 0$ 在 $$[-1, 1]$$ 上，使得问题 (2.1) 存在一个非切片的分析解 $\psi$ 。具体构造取 $\varphi$ 为方程 $\varphi'(y) = \sqrt{2(\varphi^2 - \varphi^3)}$ 的解，此时非线性项 $$F(t) = 2t - 3t^2$$ 是解析的，从而允许在分析框架下应用变分方法构造非平凡解。

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

本文的证明策略基于将流体力学问题转化为半线性椭圆方程的变分问题，并利用不同正则性类别下的结构差异来推导结果。

1. 半线性椭圆方程的转化 首先，利用远场剖面 $\varphi$ 定义涡度-流函数关系。由于 $\varphi' > 0$ ，可以在区间 $\mathcal{R}_\varphi = \varphi([-1, 1])$ 上定义函数 $F_\varphi$ ，使得 $\varphi''(y) = F_\varphi(\varphi(y))$ 。对于任意稳态 $\psi$ ，若其在某个区域内满足相同的涡度-流函数关系 $\omega = F(\psi)$ ，则流函数满足半线性椭圆方程：

\Delta \psi = F(\psi).

为了寻找非剪切解，我们设 $\psi = \varphi + u$ ，其中 $$u$$ 是小扰动。代入方程后，得到关于 $$u$$ 的方程：

\Delta u = F(\varphi + u) - F(\varphi).

定义势能函数 $G(u, y) = \Phi(\varphi + u) - \Phi(\varphi) - \Phi'(\varphi)u$ ，其中 $\Phi$ 是 $$F$$ 的原函数。则方程变为：

\Delta u = G_u(u, y).

相应的能量泛函定义为：

I(u) = \int_\Omega \left( \frac{|\nabla u|^2}{2} + G(u, y) \right) dx dy, \quad u \in H_0^1(\Omega).

2. 光滑类别下的 Min-Max 方法 在光滑类别下，关键在于 $$F$$ 在 $\mathcal{R}_\varphi$ 之外的取值是自由的。我们可以选择 $$F$$ 的一个光滑延拓到整个实轴，使得能量泛函 $$I$$ 具有合适的变分结构。

局部极小值: 根据假设 (1.3)， $F'(\varphi) > 0$ 且 $F''(\varphi) < 0$ 。这保证了剪切流 $\varphi$ 对应于能量泛函 $$I$$ 的局部极小值点（ $$u=0$$ ）。
竞争者: 我们选择 $$F$$ 的延拓使得存在另一个函数 $e \in C$ （其中 $$C$$ 是适当的函数空间），其能量低于鞍点水平，并且被能量势垒与 $$u=0$$ 分离。
Min-Max 构造: 考虑连接 $$0$$ 和 $$e$$ 的所有连续路径 $\Upsilon = \{\gamma \in C([0,1], C) \mid \gamma(0)=0, \gamma(1)=e\}$ ，定义 min-max 水平： $\sigma = \inf_{\gamma \in \Upsilon} \max_{s \in [0,1]} I(\gamma(s)).$ 在适当的紧性条件下，min-max 论证产生能量水平为 $\sigma$ 的临界点。该临界点对应于非零解 $$u$$ ，从而给出非剪切稳态 $\psi = \varphi + u$ 。

技术难点：无界域上的紧性 由于通道在水平方向上是无限的，能量泛函对水平平移不变。因此，min-max 序列可能漂移至空间无穷远。解决这一问题的关键是通过对序列进行水平重新中心（recentering）并证明重排后的序列收敛到一个非平凡的临界点。这需要建立 Sobolev 不等式（附录 A）以及热方程解的性质（附录 B）来辅助紧性论证。

3. 分析类别下的比较原理 与分析类别相比，分析函数的性质更强。如果 $\psi$ 是分析稳态，且远场剖面 $\varphi$ 也是分析的，那么 $\psi$ 在整个通道内必须满足由 $\varphi$ 唯一确定的半线性椭圆方程。也就是说，分析正则性消除了我们在光滑类别中“自由选择 $$F$$ 延拓”的自由度。

对于稠密的一类分析剪切剖面，作者建立了分析稳态的比较原理。具体来说，如果两个分析稳态具有相同的远场行为，并且在某处相等，则它们在整个域内相等。进一步地，通过最大原理或其他椭圆估计，可以证明对于这类剖面，唯一的分析解就是剪切流本身。这导致了分析类别下的刚性结果。

4. 锐利性的证明 为了证明刚性结果是锐利的，作者构造了一个具体的分析剖面 $\varphi$ ，使得对应的非线性项 $$F(t)$$ 在整个实轴上都是解析的（例如 $$F(t) = 2t - 3t^2$$ ）。在这种情况下，光滑类别中的 min-max 论证同样适用于分析类别，因为所需的变分结构在解析延拓下依然保持。因此，存在非切片的分析解，表明并非所有分析剖面都导致刚性。

4. 与之前工作的比较

本文的工作与已有的关于二维欧拉方程稳态分类的研究有显著的区别和推进，主要体现在以下几个方面：

1. 从扰动区域到一般区域的推进 早期的工作（如 [12]）主要在背景剪切流的扰动区域内研究长时间行为，证明了解弱收敛于背景剪切流。这些结果依赖于背景剪切流的均匀非停滞性以及扰动的小性。本文则直接研究稳态的存在性和分类，不依赖于时间演化或小扰动假设。通过考察稳态，本文为理解非扰动区域内的松弛机制提供了静态基础。如果稳态只能是剪切流，那么任何偏离剪切流的扰动最终必须消散；反之，如果存在稳定的非剪切稳态，则长期行为可能更加复杂。

2. 正则性类别的关键区分 以往关于稳态分类的工作（如 Hamel-Nadirashvili [14], [4], [18], [11]）通常假设速度模长远离零或流动是层流的，并在这些额外假设下证明刚性（即稳态必为剪切流）。然而，这些结果并未明确区分光滑和分析正则性。本文的核心贡献在于揭示了正则性在稳态分类中的根本性作用：

光滑类别: 存在广泛的远场剪切剖面，允许非剪切稳态的存在。这与某些物理直觉或早期刚性结果形成对比，表明在缺乏分析正则性时，流体动力学允许更丰富的稳态结构。
分析类别: 对于稠密的一类分析剖面，稳态必须是剪切流。这表明分析正则性施加了极强的约束，迫使涡度-流函数关系在整个域内保持一致。

这种“光滑柔性 vs. 分析刚性”的现象在之前的文献中未被充分探讨。虽然 [5] 指出在单连通域中，非径向分析稳态的流函数满足半线性椭圆方程，而 [7] 指出存在光滑 Morse 稳态不满足此性质，但本文将这些观察扩展到无限通道且具有特定远场行为的场景，并给出了明确的稠密刚性结果。

3. 与相混合和无粘阻尼的联系 在有限域或周期性域中，相混合（phase mixing）导致无粘阻尼，使速度场收敛于剪切流。在无限通道中，输运机制取代了相混合。本文的结果暗示，如果远场剖面属于刚性类（分析类别下的稠密集），则不存在稳定的非剪切结构，这可能支持输运机制在所有尺度上的有效性。相反，如果剖面属于柔性类，则可能存在阻碍完全松弛的稳定结构。

4. 技术方法的创新

变分方法与无界域紧性: 作者发展了在无界通道上应用 min-max 方法的技术，特别是处理水平平移不变性导致的紧性损失问题。通过结合 Sobolev 不等式和热流方法，他们成功构造了非平凡临界点。
比较原理的应用: 在分析类别中，作者利用了分析函数的唯一延拓性质和比较原理，这是光滑类别中不可用的工具。这一方法的选择直接导致了刚性结果的得出。

5. 局限性与未来方向

稠密性而非普遍性: 分析刚性仅在稠密集合上成立，而非对所有分析剖面成立。这意味着存在例外情况（如命题 2.2 所示），在这些情况下分析稳态也可以是柔性的。未来的工作可能需要刻画这些例外剖面的具体特征。
高维推广: 本文仅限于二维情况。三维欧拉方程的稳态分类更为复杂，且缺乏类似的流函数简化。将此框架推广到高维是一个开放问题。
动态稳定性: 本文仅研究稳态的存在性，未涉及这些非剪切稳态的动态稳定性。如果非剪切稳态是不稳定的，则它们可能在长时间演化中消失，从而恢复刚性图像。这是一个值得进一步研究的动力学问题。

综上所述，本文不仅在技术上解决了无限通道中欧拉方程稳态的构造与分类问题，更在概念上深化了对正则性、几何结构和流体动力学行为之间关系的理解。

4. 球对称下具有退化黏度的可压缩 Navier-Stokes 方程的膨胀解：全局存在性与无黏极限

Expanding solutions to the compressible Navier-Stokes equations with degenerate viscosities in spherical symmetry: global existence and inviscid limit

作者: Shuying Hu, Zhouping Xin, Yuan Yuan

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25394v1

摘要翻译: 本文致力于研究球对称条件下密度依赖黏度的等熵可压缩 Navier-Stokes 方程真空自由边界问题的强解，该问题描述了被真空包围的等熵可压缩黏性流动的运动。当初始数据是绝热指数 $$γ>1$$ 且黏度系数与 $$ρ^γ$$ 成正比的膨胀仿射解的小扰动时（其中 $$ρ$$ 为流体密度），我们构造了一类膨胀全局解。此外，当黏度系数趋于零时，证明了受扰动的解以黏度系数的显式收敛率收敛到可压缩 Euler 方程的解。

Comment: 无

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

1. 研究对象和背景

本文研究的是等熵可压缩 Navier-Stokes 方程 (CNS) 在真空自由边界问题 (VFBP) 下的强解, 特别是在球对称假设下, 当流体被真空包围时的运动. 具体而言, 考虑定义在随时间变化的开单连通区域 $\Omega(t) \subset \mathbb{R}^3$ 中的等熵黏性气体流动, 其动力学由以下方程组描述:

\begin{cases} \rho_t + \text{div}(\rho u) = 0 & \text{in } \Omega(t), \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \text{div} S = 0 & \text{in } \Omega(t), \\ \rho > 0 & \text{in } \Omega(t), \\ \rho = 0 \quad \text{and} \quad S \cdot n = 0 & \text{on } \Gamma(t) := \partial \Omega(t), \\ V(\Gamma(t)) = u \cdot n, \\ (\rho, u) = (\rho_0, u_0) & \text{on } \Omega := \Omega(0). \end{cases}

其中 $\rho$ 为密度, $$u$$ 为速度, $$S$$ 为应力张量, $\Gamma(t)$ 为流体与真空的界面, $$n$$ 为外法向量. 压力遵循多项式状态方程 $p = A \rho^\gamma$ , 其中 $$A>0$$ 且绝热指数 $\gamma > 1$ . 本文特别关注黏度系数依赖于密度的情形, 即 $\mu = \bar{\mu}\rho^\delta$ 和 $\lambda = \bar{\lambda}\rho^\delta$ , 其中 $\delta > 0$ . 这种设置源于气体动力学的动理论, 特别是 Chapman-Enskog 展开, 其中黏度系数是温度的幂函数, 而在等熵流中转化为密度的依赖关系.

由于黏度在真空边界处退化 ( $\rho=0$ 时 $\mu, \lambda \to 0$ ), 该问题具有双重退化特性: 既有时域上的演化退化, 又有空间耗散项的退化. 这使得标准的椭圆估计或 div-curl 论证难以直接应用. 此外, 在球对称设定下, 原点处的坐标奇异性带来了额外的技术挑战. 传统的局部适定性理论通常依赖于物理真空条件 $-\infty < \nabla_n(c_0^2) < 0$ , 其中 $c_0 = \sqrt{P'(\rho_0)}$ 为声速, 这意味着密度在真空界面附近仅具有 $C^{1/2}$ Hölder 连续性.

已有工作中, Hadžić 和 Jang [20] 以及 Shkoller 和 Sideris [59] 建立了三维不可压缩 Euler 方程在仿射解附近的全球膨胀解的存在性. 然而, 对于具有退化黏性的黏性流动, 相关理论尚不完善. 特别是, 围绕仿射解的全局存在性和无黏极限问题仍未解决. 仿射解的形式为 $r_\alpha(x,t) = \alpha(t)x$ , 其存在性要求 $\delta = \gamma$ . 本文旨在填补这一空白, 研究初始数据为仿射解小扰动时, 具有退化黏性的 CNS 方程的全局强解的存在性, 并证明当黏度系数趋于零时, 这些解收敛于对应的 Euler 方程解, 且给出显式的收敛速率.

2. 主要定理或主要结果

本文的主要结果包括两个核心定理: 全局存在性和无黏极限.

定理 2.1 (全局存在性): 假设绝热指数 $\gamma > 1$ , 黏度指数 $\delta = \gamma$ , 且满足条件 (2.9): 对于 $2\bar{\mu} + 3\bar{\lambda} \in [0, \varpi]$ , 有 $\alpha_1 > 0$ , $\alpha_0 - \frac{\varpi}{A}\alpha_1 > 0$ , 以及 $\alpha_0^{(3\gamma-1)/2} - \frac{\varpi}{A}c_1^{1/2} > 0$ . 若初始扰动 $(\eta_0, \eta_1)$ 足够小, 则存在唯一的全局强解 $(\rho, u)$ , 使得扰动变量 $\eta(x,t) := \frac{r(x,t) - r_\alpha(x,t)}{r_\alpha(x,t)}$ 满足 Lyapunov 稳定性. 具体地, 存在常数 $$C$$ 使得:

\|\eta(\cdot, t)\|_{H^s} + \|\eta_\tau(\cdot, t)\|_{H^{s-1}} \le C \left( \|\eta_0\|_{H^s} + \|\eta_1\|_{H^{s-1}} \right), \quad \forall t \ge 0,

其中 $$s$$ 为适当的 Sobolev 阶数. 这表明仿射解在扰动下是结构稳定的, 且保持近线性膨胀行为.

定理 2.2 (无黏极限): 当黏度系数 $\bar{\mu}, \bar{\lambda} \to 0$ 时, 上述黏性解 $(\rho_\epsilon, u_\epsilon)$ 收敛于对应不可压缩 Euler 方程的解 $(\rho_0, u_0)$ . 收敛速率显式依赖于黏度系数, 具体形式为:

\|\rho_\epsilon - \rho_0\|_{L^\infty([0,T]; L^q)} + \|u_\epsilon - u_0\|_{L^\infty([0,T]; L^q)} \le C \epsilon^\beta,

其中 $\epsilon$ 代表黏度尺度, $\beta > 0$ 为某个正指数, $$C$$ 独立于 $\epsilon$ . 这一结果不仅提供了黏性到无黏过渡的严格数学基础, 还作为不可压缩 Euler 方程相应结果的替代证明.

关键假设包括: 初始数据是仿射解的小扰动; 黏度系数满足 $\delta = \gamma$ ; 且参数满足特定的大小关系以确保膨胀机制的稳定作用.

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

证明的核心障碍在于处理黏度的强退化和球对称原点处的奇异性. 作者采用了以下创新方法:

非平凡微分算子的引入: 由于标准法向导数估计会产生不可控项 (如 $\partial_N^{N+1}\eta$ 和 $x\partial_N^N\eta$ ), 作者引入了一个特殊的微分算子 $\mathcal{D}_N = x\partial_x^N + (N+1)\partial_x^{N-1}$ . 该算子的设计避免了边界项的出现, 并确保对特定线性组合的估计能够单独控制高阶导数. 这一技巧是克服原点奇异性的关键.
加权能量估计与 Hardy 不等式: 尽管黏度退化导致椭圆结构减弱, 作者仍通过精细的法向导数估计结合 Hardy 不等式来恢复正则性. 与局部理论不同, 全局存在性要求估计在时间上一致或有界, 因此需要仔细平衡时间增长项和耗散项.
仿射解的稳定化机制: 利用仿射解的快速膨胀性质, 作者证明了扰动方程中存在新的稳定化机制. 这通过对流映射的 Jacobian 进行强控制来实现, 从而防止内爆奇性和自由界面的自相交.
时间变换: 引入新时间变量 $\tau = \int_0^t \frac{ds}{\alpha(s)}$ , 将原问题转化为固定边界问题, 简化了分析. 在此坐标系下, 仿射解的形式更为简单, 便于进行扰动分析.
先验估计的统一性: 建立独立于黏度系数的统一先验估计, 这是证明无黏极限的基础. 通过细致分析方程结构, 特别是退化黏性项在法向微分下的行为, 作者获得了所需的一致界.

4. 与之前工作的比较

与既有工作相比, 本文在以下几个方面取得了显著推进:

从不可压缩到黏性流动: Hadžić-Jang [20] 和 Shkoller-Sideris [59] 的结果仅限于不可压缩 Euler 方程. 本文将这些结果扩展到具有退化黏性的三维球对称黏性流动, 填补了黏性情形下的理论空白.
处理双重退化: 现有文献中, 大多数关于退化 CNS 的工作集中在弱解的全局存在性上, 而强解的理论较少. 本文首次系统性地处理了时间演化和空间耗散的双重退化问题, 并通过引入特殊微分算子克服了由此带来的技术困难.
全局存在性与无黏极限的结合: 以往工作通常分别研究全局存在性或无黏极限, 而本文同时建立了两者之间的联系, 并给出了显式的收敛速率. 这不仅深化了对黏性-无黏过渡的理解, 还为不可压缩 Euler 方程的全局解提供了新的证明路径.
球对称设定的特殊性: 虽然球对称假设限制了适用范围, 但它允许更精细的分析, 特别是在处理原点奇异性方面. 与 Luo-Xin-Zeng [46] 等工作相比, 本文的方法更侧重于法向导数估计而非切向导数, 从而更好地适应了全局存在性的需求.

局限方面, 本文的结果依赖于初始数据是仿射解的小扰动, 且黏度指数必须满足 $\delta = \gamma$ . 对于更大范围的初始数据或不同的黏度依赖关系, 结论是否成立仍需进一步研究.

5. 带体-表面相互作用的不压缩两相流扩散界面模型的 Voigt 正则化的全局准强解

Global Quasi-strong Solutions to the Voigt Regularization of a Diffuse Interface Model for Incompressible Two-phase Flows with Bulk-surface Interaction

作者: Patrik Knopf, Maoyin Lv, Hao Wu

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25578v1

摘要翻译: 我们分析了在光滑有界区域 $Ω\subset\mathbb{R}^d$ $$(d=2,3)$$ 中，密度不匹配的不压缩两相粘性流动的热力学一致扩散界面模型。该模型由速度场的 Navier-Stokes-Voigt 方程和相场变量的带有奇异势的对流 Cahn-Hilliard 方程组成。所得到的流体动力学系统对速度场施加广义 Navier 滑移边界条件，并对相场变量和化学势施加 Cahn-Hilliard 型动态边界条件。借助 Voigt 正则化，我们建立了满足能量等式的全局准强解的存在性。这是通过巧妙的近似方案与半 Galerkin 方法的结合实现的。

Comment: 无

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

1. 研究对象和背景

本文研究的是在光滑有界区域 $\Omega \subset \mathbb{R}^d$ ( $$d=2,3$$ ) 中，具有不相等密度的不可压缩粘性两相流的扩散界面模型。该模型旨在克服传统无滑移边界条件和不均匀Neumann边界条件在描述接触线动力学、动态接触角以及体-表面物质交换方面的局限性。

研究对象是一个由Navier-Stokes-Voigt方程和带有奇异势的平流Cahn-Hilliard方程组成的耦合系统。设 $v: \Omega \times (0,T) \to \mathbb{R}^d$ 为体积平均速度场， $\phi: \Omega \times (0,T) \to [-1,1]$ 为序参量（即两种流体体积分数之差）， $$p$$ 为压力。密度 $\rho(\phi)$ 和粘度 $\nu(\phi)$ 依赖于相场变量，具体形式为：

\rho(\phi) = \rho_1 \frac{1+\phi}{2} + \rho_2 \frac{1-\phi}{2}, \quad \nu(\phi) = \nu_1 \frac{1+\phi}{2} + \nu_2 \frac{1-\phi}{2}

其中 $\rho_1 > \rho_2 > 0$ 且 $\nu_i > 0$ 。化学势 $\mu$ 由 $\mu = -\Delta \phi + F'(\phi)$ 给出，其中 $$F$$ 是具有双井结构的自由能密度函数，物理上最相关的选择是Flory-Huggins势：

F_{log}(r) = \frac{\Theta}{2} [(1+r)\ln(1+r) + (1-r)\ln(1-r)] - \frac{\Theta_0}{2} r^2, \quad 0 < \Theta < \Theta_0.

为了处理不相等密度的情况并改善时间正则性，作者引入了Voigt正则化项。修正后的应力张量为 $T(v) = 2\nu(\phi)Dv + 2\alpha D\partial_t v$ ，其中 $\alpha > 0$ 是正则化参数。系统的主体方程如下：

\begin{cases} \partial_t(\rho(\phi)v) + \text{div}\left(v \otimes (\rho(\phi)v + J)\right) - \text{div}\left(2\nu(\phi)Dv + 2\alpha D\partial_t v\right) + \nabla p = \mu \nabla \phi & \text{in } Q_T, \\ \text{div } v = 0 & \text{in } Q_T, \\ \partial_t \phi + v \cdot \nabla \phi = \Delta \mu & \text{in } Q_T, \\ \mu = -\Delta \phi + F'(\phi) & \text{in } Q_T, \end{cases}

其中通量项 $J = -\rho' \nabla \mu$ ，且 $\rho' = \frac{\rho_1 - \rho_2}{2}$ 。

与以往工作不同，本文考虑了一类新的边界条件，允许体与表面之间的物质传输以及动态接触角的存在。这些边界条件基于局部质量守恒定律、能量耗散原理和拉格朗日乘子法推导而来。在边界 $\Gamma = \partial \Omega$ 上，引入表面相场变量 $\psi$ 和表面化学势 $$L$$ 。边界条件包括：

速度的广义Navier滑移条件，涉及一个作为拉格朗日乘子的表面压力 $$q$$ ，以强制满足 $\text{div}_\Gamma v_\tau = 0$ 。
相场变量 $\phi$ 的Robin型或Neumann型边界条件：若 $K \in [0, +\infty)$ ，则 $K\partial_n \phi = \psi - \phi$ ；若 $K = +\infty$ ，则 $\partial_n \phi = 0$ 。
化学势 $\mu$ 的Robin型或Neumann型边界条件：若 $L \in [0, +\infty)$ ，则 $L\partial_n \mu = L - \mu$ ；若 $L = +\infty$ ，则 $\partial_n \mu = 0$ 。
表面相场 $\psi$ 的动态边界条件： $\partial_t \psi + v_\tau \cdot \nabla_\tau \psi = \Delta_\tau L - \partial_n \mu$ ，其中 $L = -\Delta_\tau \psi + \partial_n \phi + G'(\psi)$ 。

系统的总能量 $E_{tot}$ 包含动能、体自由能和表面自由能，并满足能量恒等式。对于足够正则的解，该系统还满足体-表面总质量守恒性质。本文的核心目标是证明该Voigt正则化系统在二维和三维空间中全局准强解的存在性，并验证其满足能量恒等式。这一结果填补了在不相等密度条件下，结合复杂动态边界条件时，高于弱解正则性的理论空白。

2. 主要定理或主要结果

本文的主要结果是定理2.1，证明了上述Voigt正则化的Navier-Stokes-Cahn-Hilliard系统在全局时间尺度上的准强解存在性。

主要结论： 假设初始数据 $(v_0, \phi_0, \psi_0)$ 满足适当的相容性条件，且势函数 $$F$$ 和 $$G$$ 满足特定的凸性和增长条件（例如Flory-Huggins型奇异势）。对于任意 $\alpha > 0$ 和参数 $K \in (0, +\infty]$ , $L \in [0, +\infty]$ ，存在一个定义在 $$[0, T]$$ 上的全局准强解 $(v, \phi, \psi, \mu, L, q, p)$ ，使得：

正则性：
- 速度场 $$v$$ 具有较高的时间正则性，特别是 $\partial_t v \in L^2(0, T; H^1(\Omega))$ 。
- 相场变量 $\phi$ 和表面相场 $\psi$ 满足严格的分离性质，即它们的值始终位于区间 $(-1+\sigma, 1-\sigma)$ 内，从而避免了奇异势在 $\pm 1$ 处的奇异性影响估计。
- 化学势 $\mu$ 和表面化学势 $$L$$ 的时间导数满足 $(\partial_t \mu, \partial_t L) \in L^2(0, T; H^1)$ 。
能量恒等式： 解严格满足以下能量耗散关系：
$\frac{d}{dt} E_{tot}(v(t), \phi(t)) + \int_\Omega 2\nu(\phi)|Dv|^2 dx + \int_\Gamma \beta(\psi)|v_\tau|^2 dS + \int_\Omega |\nabla \mu|^2 dx + \int_\Gamma |\nabla_\tau L|^2 dS + \chi(L) \int_\Gamma |\mu - L|^2 dS = 0.$
其中 $\chi(L)$ 是一个指示函数，当 $L \in (0, +\infty)$ 时为 $$1/L$$ ，否则为0。
质量守恒： 解满足体-表面耦合的质量守恒律。具体地，若 $L \in [0, +\infty)$ ，则 $\int_\Omega \phi(t) dx + \int_\Gamma \psi(t) dS$ 为常数；若 $L = +\infty$ ，则体质量和表面质量分别守恒。

技术细节：

证明首先针对情形 $(K, L) \in (0, +\infty] \times (0, +\infty)$ 进行，利用半Galerkin逼近方案构造近似解序列。
通过建立关于逼近参数的一致先验估计，特别是利用Voigt项提供的 $\partial_t v$ 的 $$H^1$$ 正则性，控制边界项中的高阶非线性项。
最后，通过紧性论证取极限得到原问题的解。对于 $$L=0$$ 和 $L=+\infty$ 的情形，通过对相应渐近极限的分析得出存在性结果。

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

证明的核心障碍在于处理不相等密度带来的低正则性以及复杂动态边界条件引起的非线性耦合项。作者采用了以下关键技术和创新点：

Voigt正则化的利用： Voigt项 $\text{div}(2\alpha D\partial_t v)$ 提供了 $\partial_t v$ 在 $H^1(\Omega)$ 中的估计。这是获得“准强”正则性的关键。在推导先验估计时，作者选取 $w = \partial_t v$ 作为速度方程的测试函数。这导致边界上出现涉及 $\partial_t v$ 的积分项，如 $\int_\Gamma \beta(\psi)v_\tau \cdot \partial_t v_\tau dS$ 。由于 $\alpha > 0$ ，可以利用Sobolev嵌入定理将 $\|\partial_t v\|_{L^4(\Gamma)}$ 控制在内部能量项 $\int_\Omega \rho(\phi)|\partial_t v|^2 dx + 2\alpha \int_\Omega |D\partial_t v|^2 dx$ 之内。这种技巧使得能够闭合能量估计，从而提升时间正则性。
精细的逼近方案与半Galerkin方法： 直接求解原问题困难重重，因此作者构建了一个多层逼近过程。
- 初值正则化： 对初始数据进行平滑处理，以确保后续估计的有效性。
- 粘性正则化： 对体-表面平流Cahn-Hilliard子系统引入额外的粘性正则化项。这一步至关重要，因为它允许作者严谨地计算高阶Sobolev估计，特别是获得 $(\partial_t \mu, \partial_t L) \in L^2(0, T; H^1)$ 的正则性。
- 半Galerkin逼近： 在有限维子空间上求解逼近问题，确保近似解的存在性。
严格分离性质（Strict Separation Property）的证明： 为了处理奇异势 $$F$$ 和 $$G$$ （其在 $\pm 1$ 处发散），作者引入了一个辅助的奇异势 $F_\sigma$ ，其奇点位于 $\pm(1-\sigma)$ 。通过比较原理，证明了相场变量 $\phi$ 和 $\psi$ 的值被限制在 $(-1+\sigma, 1-\sigma)$ 内。这使得在该区间内，势函数的导数可以被视为正则函数，从而避免了奇异性对微分方程分析的干扰。这一技巧借鉴了文献中的方法，但在此处结合动态边界条件进行了适配。
边界项的控制与拉格朗日乘子： 边界条件中引入了表面压力 $$q$$ 作为拉格朗日乘子，以强制满足 $\text{div}_\Gamma v_\tau = 0$ 。这在分析中简化了边界项的处理，因为不需要显式求解 $$q$$ ，而是将其视为保证边界兼容性存在的未知量。此外，Robin型边界条件 $K\partial_n \phi = \psi - \phi$ 在 $$K > 0$$ 时提供了 $\partial_n \phi$ 的高阶正则性，这对于推导 $\partial_t \mu$ 的 $$H^1$$ 正则性是不可或缺的。
紧性论证与极限过程： 在获得一致估计后，作者使用Aubin-Lions引理等紧性工具，从逼近解序列中提取收敛子列。特别地，需要处理非线性项的强收敛性，这通常依赖于相场的严格分离性质和速度场的时间正则性。最后，通过让逼近参数趋于零或无穷大，得到原问题的解。

4. 与之前工作的比较

本文将现有理论与以下方面进行了对比和推进：

与经典AGG模型及Model H的比较： 传统的AGG模型（Abels-Garcke-Grün）通常假设无滑移边界条件和齐次Neumann边界条件。这些条件限制了模型对移动接触线和动态接触角的描述能力。本文采用的边界条件（1.7）允许表面扩散、动态接触角以及体-表面间的物质传输，因此在物理适用性上显著优于经典模型。此外，经典模型在不相等密度情况下通常只能证明弱解的存在性，而本文通过Voigt正则化获得了更高正则性的准强解。
与先前动态边界条件工作的比较： 文献[31]和[24]研究了带有Cahn-Hilliard型动态边界条件的模型。然而，这些工作主要关注弱解的存在性，且在处理不相等密度时面临巨大挑战。本文的创新在于：
- 引入Voigt正则化： 这是获得时间方向更高正则性的关键手段。先前的工作大多没有利用这一正则化来获得准强解。
- 新的边界条件结构： 本文考虑的边界条件（1.7）比[31]中的（1.5）多了一个表面压力项 $$q$$ ，并强制 $\text{div}_\Gamma v_\tau = 0$ 。这种结构在数学分析上更为复杂，但也更符合某些物理情境下的动量平衡。
- 适用范围： 本文的结果适用于 $K \in (0, +\infty]$ 和 $L \in [0, +\infty]$ 的广泛参数范围，并通过渐近分析覆盖了极端情形。
与Voigt正则化其他应用的比较： Voigt正则化常用于粘弹性流体或湍流模型的数学分析。在不可压缩两相流领域，其应用相对较少。本文展示了如何利用Voigt项提供的额外耗散来克服非线性耦合项的低正则性问题，特别是在边界上的交叉项。这与单纯依靠能量估计的方法不同，提供了一种新的正则化路径。
局限性与代价：
- 正则化依赖： 结果依赖于Voigt参数 $\alpha > 0$ 。当 $\alpha \to 0$ 时，能否恢复原始Navier-Stokes系统的解是未来的研究方向。目前，本文仅证明了正则化系统的全局存在性，并未直接解决原始系统的正则性问题。
- 奇异势的限制： 虽然使用了严格分离技巧，但仍要求势函数满足特定的凸性和增长条件。对于更一般的非凸势或非标准增长条件，结果可能不成立。
- 唯一性： 本文未讨论解的唯一性。在高维（ $$d=3$$ ）非线性PDE中，唯一性通常是开放问题。

6. 三维 $rac{1}{2}|u|^2+p$ 的几何细化

Geometric refinements of Liouville-type theorems for the stationary Navier–Stokes equations in $\mathbb{R}^3$

作者: Juhyeong Lee

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25342v1

摘要翻译: 我们证明了在 $rac{1}{2}|u|^2+p$ 条件下，满足有限 Dirichlet 能量条件和一致衰减条件的光滑解 $$(u,p)$$ 的 Liouville 型定理。我们的结果对两个近期的 Osgood 型判据进行了几何细化，即以头部压力 $Q= rac{1}{2}|u|^2+p$ 表述的相对衰减判据和加权可积性判据。证明利用了头部压力的若干性质以及 $(u imes rac{1}{2}|u|^2+p)

Comment: 10 pages

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

1. 研究对象和背景

本文研究三维欧氏空间 $\mathbb{R}^3$ 中稳态 Navier-Stokes 方程的 Liouville 型问题。具体而言，我们考虑如下方程组：

-\Delta u + (u \cdot \nabla)u + \nabla p = 0, \quad \nabla \cdot u = 0,

其中 $u: \mathbb{R}^3 \to \mathbb{R}^3$ 是速度场， $p: \mathbb{R}^3 \to \mathbb{R}$ 是压力。Liouville 问题的核心在于：如果一个光滑解 $$(u, p)$$ 满足无穷远处的均匀衰减条件

u(x) \to 0 \quad \text{当 } |x| \to +\infty,

以及有限的 Dirichlet 能量条件

\int_{\mathbb{R}^3} |\nabla u|^2 dx < +\infty,

那么该解是否必然为平凡解（即 $u \equiv 0$ ）？

这一问题在数学流体力学中具有重要意义。Galdi 在其专著中指出，对于 $n \geq 4$ 维情况，结合 Sobolev 嵌入和衰减条件可以证明解的平凡性；二维情况已由 Gilbarg 和 Weinberger 利用涡量方程的最大值原理解决。然而，三维情形至今仍未完全解决，是该领域的一个著名开放问题。

近年来，基于“头压”（head pressure） $Q := \frac{1}{2}|u|^2 + p$ 的方法成为研究此问题的有力工具。头压满足椭圆方程

\Delta Q - u \cdot \nabla Q = |\omega|^2,

其中 $\omega = \nabla \times u$ 是涡量。Chae 等人发展了基于水平集方法和 coarea 公式的 Osgood 型判据。例如，Chae [7] 证明了如果存在满足 Osgood 条件 $\int_0^1 \frac{ds}{\xi(s)} = +\infty$ 的函数 $\xi$ ，使得

\frac{|u|^2}{\xi(|Q|)^{\frac{3-2q}{q}}} \in L^{\frac{q}{q-1}}(\mathbb{R}^3), \quad \frac{|\nabla u|}{\xi(|Q|)^{\frac{3q-3}{2q}}} \in L^{2q}(\mathbb{R}^3),

则解必为平凡。此外，还导出了相对衰减判据，如

\sup_{x \in \mathbb{R}^3} \frac{|u(x)|^2}{\xi(|Q(x)|)} < +\infty \quad \text{或} \quad \sup_{x \in \mathbb{R}^3} \frac{|p(x)|}{\xi(|Q(x)|)} < +\infty.

尽管这些结果推进了对问题的理解，但它们通常涉及速度场 $$u$$ 或其梯度 $\nabla u$ 的整体范数，未能充分利用方程中的几何结构。特别是，标量三重积 $(u \times \omega) \cdot \nabla Q$ 在 Lamb 形式 $\nabla Q = \Delta u + u \times \omega$ 中扮演关键角色，但其几何意义未被充分挖掘。本文旨在通过引入沿 $$Q$$ 的水平面的切向分解和角度因子，对上述 Osgood 型判据进行几何细化，从而放宽对解的增长率或可积性的要求，特别是在允许亚临界增长方面取得进展。

2. 主要定理或主要结果

本文的主要结果是两个关于 Liouville 型定理的几何细化版本，分别对应于相对衰减设定和加权可积性设定。

定理 1.1 (Osgood 尺度小性判据): 设 $$(u, p)$$ 是满足均匀衰减条件和有限 Dirichlet 能量条件的稳态 Navier-Stokes 方程的光滑解，且 $\xi$ 满足 Osgood 条件。假设以下任一条件成立：

速度条件: $\liminf_{\lambda \to 0} \left( \int_\lambda^1 \frac{ds}{\xi(s)} \right)^{-1} \left\| \frac{|u_T|^2 \Theta^2}{\xi(|Q|)} \right\|_{L^\infty(\{|Q|>\lambda\})} < 1,$
涡量条件: 存在足够小的常数 $$c > 0$$ ，使得 $\liminf_{\lambda \to 0} \left( \int_\lambda^1 \frac{ds}{\xi(s)} \right)^{-1} \left\| \frac{|\omega_T|^2 \Theta^2}{\xi(|Q|)} \right\|_{L^{3/2}(\{|Q|>\lambda\})} < c.$ 则 $u \equiv 0$ 。

这里， $$u_T$$ 和 $\omega_T$ 分别是 $$u$$ 和 $\omega$ 在 $$Q$$ 的水平面上的切向分量， $\Theta$ 是定义在 $\{|u_T||\omega_T||\nabla Q| > 0\}$ 上的角度因子：

\Theta := \frac{(u_T \times \omega_T) \cdot n_Q}{|u_T||\omega_T|},

其中 $n_Q = \nabla Q / |\nabla Q|$ 。注意 $-1 \leq \Theta \leq 1$ ，且 $(u \times \omega) \cdot \nabla Q = |u_T||\omega_T||\nabla Q|\Theta$ 。

推论 (压力替代条件): 由于在非平凡情况下 $$Q < 0$$ ，有 $|u_T|^2 \leq |u|^2 \leq 2|p|$ 。因此，上述速度条件可被以下压力条件替代：

\liminf_{\lambda \to 0} \left( \int_\lambda^1 \frac{ds}{\xi(s)} \right)^{-1} \left\| \frac{|p|\Theta^2}{\xi(|Q|)} \right\|_{L^\infty(\{|Q|>\lambda\})} < \frac{1}{2}.

定理 1.2 (加权可积性判据): 设 $$(u, p)$$ 是满足相同基本假设的光滑解，且 $\xi$ 满足 Osgood 条件。假设

\frac{|u_T||\Theta|}{\xi(|Q|)^a} \in L^q(\mathbb{R}^3), \quad \frac{|\omega_T|}{\xi(|Q|)^b} \in L^r(\mathbb{R}^3),

其中 $\frac{1}{q} + \frac{1}{r} = \frac{1}{2}$ ， $2 \leq q, r \leq +\infty$ ，且 $a + b = \frac{1}{2}$ ， $a, b \in \mathbb{R}$ 。则 $u \equiv 0$ 。

与先前结果的比较: 定理 1.1 的条件比 Chae [7] 中的相对衰减条件更弱。后者要求 $\| \frac{|u|^2}{\xi(|Q|)} \|_{L^\infty(\{|Q|>\lambda\})} = O(1)$ 当 $\lambda \to 0$ ，而本文允许在 Osgood 尺度下的亚临界增长，即存在序列 $\lambda_j \to 0$ 使得

\left\| \frac{|u_T|^2 \Theta^2}{\xi(|Q|)} \right\|_{L^\infty(\{|Q|>\lambda_j\})} \leq c \int_{\lambda_j}^1 \frac{ds}{\xi(s)}, \quad c < 1.

类似地，定理 1.2 细化了 Chae [7] 的加权可积性条件，仅涉及切向分量和角度因子，而非全速度或全梯度。

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

本文的证明核心在于利用头压 $$Q$$ 的几何结构和标量三重积 $(u \times \omega) \cdot \nabla Q$ 的分解。主要技术步骤如下：

头压的性质与水平集分析: 首先，利用最大值原理确定非平凡解的头压 $$Q$$ 在 $\mathbb{R}^3$ 上严格为负，且趋于 0。这使得超水平集 $\Omega_\lambda = \{x \in \mathbb{R}^3 : |Q(x)| > \lambda\}$ 有界，且由 Sard 定理，对于几乎处处的 $\lambda$ ，其边界是光滑的，从而可以应用散度定理。
Osgood 型表示: 引用 Chae [7] 的引理 2.1，它建立了总涡量能量与水平集积分之间的关系：
$\lim_{\lambda \to 0} \left( \int_\lambda^1 \frac{ds}{\xi(s)} \right)^{-1} \int_{\{|Q|>\lambda\}} \frac{|\nabla Q|^2}{\xi(|Q|)} dx = \int_{\mathbb{R}^3} |\omega|^2 dx.$
这一恒等式是将局部估计转化为全局矛盾的关键桥梁。
切向分解与角度因子: 引入法向量 $n_Q = \nabla Q / |\nabla Q|$ ，并将速度和涡量分解为法向和切向部分：
$u_T = u - (u \cdot n_Q)n_Q, \quad \omega_T = \omega - (\omega \cdot n_Q)n_Q.$
定义角度因子 $\Theta$ ，使得 $(u \times \omega) \cdot \nabla Q = |u_T||\omega_T||\nabla Q|\Theta$ 。这一分解揭示了非线性项的几何本质：只有垂直于 $$Q$$ 水平面的相互作用才贡献于能量平衡。
Lamb 形式的测试与 H"older 不等式: 将 Lamb 形式 $\nabla Q = \Delta u + u \times \omega$ 乘以 $\nabla Q / \xi(|Q|)$ 并在 $\Omega_\lambda$ 上积分。利用 $\Delta u$ 的无散度性质，消去扩散项，得到：
$\int_{\Omega_\lambda} \frac{|\nabla Q|^2}{\xi(|Q|)} dx = \int_{\Omega_\lambda} \frac{(u \times \omega) \cdot \nabla Q}{\xi(|Q|)} dx = \int_{\Omega_\lambda} \frac{|u_T||\omega_T||\nabla Q|\Theta}{\xi(|Q|)} dx.$
随后，应用 H"older 不等式，将右端拆分为三个因子的乘积，并选择适当的指数 $$q, r$$ 和权重 $$a, b$$ 以满足缩放关系。这导致了引理 2.2 中的关键估计：
$\int_{\mathbb{R}^3} |\omega|^2 dx \leq \liminf_{\lambda \to 0} \left( \int_\lambda^1 \frac{ds}{\xi(s)} \right)^{-1} \left\| \frac{|u_T||\Theta|}{\xi(|Q|)^a} \right\|_{L^q}^2 \left\| \frac{|\omega_T|}{\xi(|Q|)^b} \right\|_{L^r}^2.$
矛盾论证: 在定理 1.1 的证明中，假设解非平凡，则 $\int |\omega|^2 dx > 0$ 。利用上述估计和给定的小性条件，导出形如 $E \leq \alpha E$ 的不等式，其中 $E = \int |\omega|^2 dx$ 且 $\alpha < 1$ ，从而产生矛盾。在定理 1.2 中，利用全局可积性和 Osgood 积分的发散性，直接导出 $E \leq 0$ ，同样导致矛盾。

4. 与之前工作的比较

本文与既有工作的主要区别体现在以下几个方面：

从整体范数到几何切向量的细化: Chae [7] 的 Osgood 型判据涉及全速度 $$|u|$$ 和全梯度 $|\nabla u|$ （或涡量 $|\omega|$ ）。例如，条件 $\frac{|u|^2}{\xi(|Q|)^{\dots}} \in L^{\dots}$ 要求在整个空间中对速度的控制。相比之下，本文的结果仅依赖于速度 $$u$$ 和涡量 $\omega$ 在 $$Q$$ 的水平面上的切向分量 $$u_T$$ 和 $\omega_T$ 。由于 $|u_T| \leq |u|$ 且 $|\omega_T| \leq |\omega|$ ，这些条件是严格更弱的。这意味着即使速度场在某些方向上有较大的法向分量，只要其在水平面上的切向相互作用足够小，解仍必须是平凡的。
引入角度因子 $\Theta$ : 本文首次明确引入了角度因子 $\Theta$ ，它度量了切向速度 $$u_T$$ 和切向涡量 $\omega_T$ 之间的相对取向。标量三重积 $(u \times \omega) \cdot \nabla Q$ 的大小不仅取决于 $$|u_T|$$ 和 $|\omega_T|$ 的模长，还取决于它们是否正交于法向量 $$n_Q$$ 以及彼此之间的夹角。如果 $$u_T$$ 和 $\omega_T$ 平行（即 $\Theta = 0$ ），则非线性项对能量平衡没有贡献。这一几何洞察使得判据更加精细，能够捕捉到流动结构中与能量耗散直接相关的部分。
允许亚临界增长: 在相对衰减设定下，Chae [7] 的条件要求 $\frac{|u|^2}{\xi(|Q|)}$ 在 $\{|Q|>\lambda\}$ 上的一致范数为 $$O(1)$$ 当 $\lambda \to 0$ 。本文的定理 1.1 允许该量以 Osgood 尺度 $\int_\lambda^1 \frac{ds}{\xi(s)}$ 的速度增长，只要其系数小于 1（或压力情况下小于 1/2）。这是一种显著的放宽，因为它允许解在接近 $$Q=0$$ 的区域有更剧烈的行为，只要这种增长被 Osgood 积分的发散性所抑制。
加权可积性的灵活性: 定理 1.2 提供了更一般的加权可积性框架，参数 $$a, b, q, r$$ 满足 $$a+b=1/2$$ 和 $$1/q+1/r=1/2$$ 。这不仅涵盖了 Chae [7] 中的特定情形（如取特定的 $$s$$ 值），还允许在不同的函数空间中进行权衡。例如，可以将更多的权重分配给 $$u_T$$ 或 $\omega_T$ ，或者选择不同的可积性指数，从而适应不同正则性假设下的解。
局限性: 尽管本文的判据在几何上更为精细，但它们仍然依赖于头压 $$Q$$ 的结构。如果 $$Q$$ 的水平面非常复杂或不规则，计算或验证这些条件可能具有挑战性。此外，这些结果仍然是 Liouville 型的，即证明解的平凡性，而不是提供解的存在性或稳定性。对于更一般的非平稳或高维情况，这些方法可能需要进一步的推广。

7. 可压缩 Navier-Stokes 方程 ALE 格式中熵稳定的移动壁面边界条件

Entropy-stable moving-wall boundary conditions for the ALE formulation of the compressible Navier-Stokes equations

作者: Luca Galimberti, Roberto Nuca, Lisandro Dalcin, Alberto Guardone, Matteo Parsani

类别: math.NA, math.AP

发布日期: 2026-08-25

arXiv ID: 2608.25146v1

摘要翻译: 我们提出了一个用于运动域上可压缩 Euler 方程和 Navier-Stokes 方程的高阶熵稳定框架。描述域运动的时空映射被重新表述为任意拉格朗日-欧拉 (ALE) 形式，其中物理无粘通量和由网格运动引起的贡献以统一的方式处理。在连续层面，我们证明了所提出的移动壁面边界条件对 Euler 方程是熵守恒的，对 Navier-Stokes 方程是熵稳定的。无滑移条件以相对于运动壁面的速度来表述，从而产生有界的无粘项对熵平衡的贡献，而粘性项仅贡献熵耗散。利用对角范数分部求和 (SBP) 算子及适当的数值通量，这些性质被扩展到半离散格式，从而得到 $$L^2$$ 意义下的非线性稳定性。通过一系列广泛的数值实验展示了所提方法的精度、鲁棒性和可扩展性，范围从经典二维验证算例到涉及运动边界和流固相互作用的大规模湍流和高超声速模拟。结果证实了高阶熵稳定格式适用于广泛流动机制和多物理场应用中的复杂运动域问题。由于分析依赖于 SBP 性质而非特定离散化，该框架自然推广到基于对角范数 SBP 算子的广泛方法类，包括有限体积法、有限元法和通量重构方案。

Comment: 57 pages, 15 figures

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

1. 研究对象和背景

在计算流体力学中，任意拉格朗日-欧拉（ALE）框架为模拟运动和变形域上的流动提供了自然的数学基础。这类问题广泛存在于包含运动或旋转部件、流固耦合以及预设体运动的工程配置中。从计算角度看，ALE 方法通过保持贴体离散化并允许自由度的重新分配（即 r-适应），能够显著改善感兴趣区域（特别是固体壁面和边界层附近）的分辨率分布。

本文研究的核心对象是定义在运动曲率网格上的可压缩 Navier-Stokes 方程。设物理域由时空坐标 $[t, \mathbf{x}]$ 描述，参考域由 $[\tau, \boldsymbol{\xi}]$ 描述。映射 $\chi(\tau, \boldsymbol{\xi})$ 将参考域映射到物理域，其逆映射为 $\boldsymbol{\xi}(t, \mathbf{x})$ 。网格速度定义为 $\mathbf{V} = \frac{\partial \chi}{\partial \tau}$ 。基于该映射，原始的笛卡尔坐标系下的守恒变量方程被转化为 ALE 形式：

\frac{\partial (J \mathbf{Q})}{\partial \tau} + \sum_{l,m=1}^{3} \frac{\partial}{\partial \xi_l} \left( J \frac{\partial \xi_l}{\partial x_m} (\mathbf{F}_{xm}^{(ALE)} - \mathbf{F}_{xm}^{(V)}) \right) = 0

其中 $$J$$ 是映射的雅可比行列式， $\mathbf{Q}$ 是守恒变量向量， $\mathbf{F}_{xm}^{(I)}$ 是无粘通量， $\mathbf{F}_{xm}^{(V)}$ 是粘性通量。ALE 通量定义为 $\mathbf{F}_{xm}^{(ALE)} = \mathbf{F}_{xm}^{(I)} - V_m \mathbf{Q}$ ，它将无粘通量和网格运动贡献统一处理。

对于高阶数值格式，熵稳定性提供了一种与热力学第二定律一致的强非线性稳定性机制。然而，在 ALE 设置下，熵平衡耦合了物理通量、网格速度和坐标变换引入的度量项。现有的熵稳定壁面处理方法主要针对静态网格开发，无法直接推广到一般的 ALE 边界。尽管文献中已提出了一些针对可压缩 Navier-Stokes 方程的熵保守和熵稳定 ALE 公式 [10, 11]，但这些分析要么局限于周期性域，要么未对移动壁面边界条件的处理提供完整的分析和证明。这一空白使得在处理复杂的运动边界问题时，缺乏具有严格稳定性保证的高阶数值方案。本文旨在填补这一空白，推导连续水平下移动壁面边界条件的熵保守性质，并在半离散 DG-SBP-SAT 框架内证明其熵稳定性和熵保守性。

2. 主要定理或主要结果

本文的主要结果包括连续水平和半离散水平的熵稳定性理论证明以及广泛的数值验证。

连续熵平衡推导：作者推导了运动壁面边界的连续熵平衡方程。通过收缩 ALE 系统与熵变量 $\mathbf{W} = \frac{\partial S}{\partial \mathbf{Q}}$ ，得到了标量熵方程的微分形式：
$\frac{\partial (J S)}{\partial \tau} = -\sum_{l,m=1}^{3} \frac{\partial}{\partial \xi_l} \left[ J \frac{\partial \xi_l}{\partial x_m} (F_{xm} - V_m S) \right] + \sum_{l,n=1}^{3} \left[ \frac{\partial}{\partial \xi_l} \left( \mathbf{W}^T \hat{C}_{ln} \frac{\partial \mathbf{W}}{\partial \xi_n} \right) - \frac{\partial \mathbf{W}^T}{\partial \xi_l} \hat{C}_{ln} \frac{\partial \mathbf{W}}{\partial \xi_n} \right]$
其中后两项中的度量恒等式项在满足几何守恒律时消失。积分后得到全局熵守恒陈述：
$\frac{d}{d\tau} \int_{\Omega} J S d\Omega \leq \sum_{l,m=1}^{3} \int_{\Gamma} \left( \mathbf{W}^T \mathbf{F}_{xm}^{(V)} - J \frac{\partial \xi_l}{\partial x_m} (F_{xm} - V_m S) \right) n_{\xi_l} d\Gamma - D_T$
其中 $D_T \geq 0$ 是粘性耗散项。这表明数学熵的增加只能来源于边界贡献，而粘性项仅贡献熵耗散。
移动壁面边界条件的熵性质：对于无滑移绝热移动壁面，边界条件以相对于移动壁面的速度形式表述。作者证明了这种边界条件对于 Euler 方程是熵保守的，而对于 Navier-Stokes 方程是熵稳定的。具体而言，无粘部分对熵平衡的贡献是有界的，而粘性部分仅产生熵耗散。
半离散熵稳定性：使用对角范数求和分部（SBP）算子和适当的数值通量，上述性质被扩展到半离散公式中，导致 $$L^2$$ 意义下的非线性稳定性。这一构造不依赖于特定的离散化方法，而是基于 SBP 属性，因此自然适用于有限体积、有限元和通量重构等基于一类对角范数 SBP 算子的方法。
数值验证：通过一系列数值实验展示了方法的精度、鲁棒性和可扩展性，包括规范二维验证案例、大规模湍流和超音速模拟，涉及运动边界和流固相互作用。结果证实了高阶熵稳定方案在复杂运动域问题和多种流态及多物理场应用中的适用性。

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

本文的核心技术障碍在于如何在 ALE 框架下统一处理物理通量和网格运动诱导的贡献，同时保持熵稳定性。作者的真正创新点在于推导了适用于一般运动和变形 ALE 边界的熵保守边界条件，并在半离散水平上证明了其稳定性。

ALE 公式的统一处理：作者将空间-时间映射重新表述为 ALE 形式，其中物理无粘通量和网格运动贡献通过 ALE 通量 $\mathbf{F}_{xm}^{(ALE)} = \mathbf{F}_{xm}^{(I)} - V_m \mathbf{Q}$ 统一处理。这种统一处理简化了熵平衡的推导，使得可以清晰地分离出无粘、网格运动和粘性部分的贡献。
熵变量和势函数：利用 Godunov 定理，定义了熵势 $\Phi(\mathbf{W}) = \mathbf{W}^T \mathbf{Q}(\mathbf{W}) - S(\mathbf{Q}(\mathbf{W}))$ 和熵通量势 $\Psi_{xm}(\mathbf{W}) = \mathbf{W}^T \mathbf{F}_{xm}^{(I)}(\mathbf{Q}(\mathbf{W})) - F_{xm}(\mathbf{Q}(\mathbf{W}))$ 。这些势函数在推导熵平衡时起到了关键作用，使得可以将熵方程写成守恒形式。
SBP 算子和 SAT 技术：在半离散水平上，作者使用对角范数 SBP 算子来离散导数项，并使用同时逼近项（SAT）技术来处理边界条件。SBP 算子的性质确保了离散算子与连续算子在熵意义上的兼容性。通过选择适当的数值通量和 SAT 参数，作者证明了半离散格式的熵保守性和熵稳定性。
几何守恒律：作者在推导过程中利用了时间度量恒等式 $\frac{\partial J}{\partial \tau} - \sum_{l,m=1}^{3} \frac{\partial}{\partial \xi_l} \left( J \frac{\partial \xi_l}{\partial x_m} V_m \right) = 0$ 和空间度量恒等式 $\sum_{l=1}^{3} \frac{\partial}{\partial \xi_l} \left( J \frac{\partial \xi_l}{\partial x_m} \right) = 0$ 。这些恒等式保证了在运动网格上守恒律的正确性，并且在熵平衡推导中消除了额外的源项。
正定性保持：虽然本文主要关注熵稳定性，但作者指出熵估计的解释依赖于密度和温度的正定性。在实际应用中，通常需要额外的正定性保持机制来确保解保持在物理可行集内。这一点在讨论中有所提及，但未作为主要证明内容展开。

4. 与之前工作的比较

本文与既有工作的主要区别在于其对移动壁面边界条件的完整熵稳定性分析，以及对 ALE 框架下统一处理的深入探讨。

与静态网格方法的比较：现有的熵稳定壁面处理方法主要针对静态网格开发，例如 [15, 16, 17, 18, 19, 20, 21, 22, 23, 24]。这些方法在处理运动边界时无法直接应用，因为熵平衡中引入了网格速度和度量项。本文通过推导适用于一般 ALE 边界的熵保守边界条件，解决了这一问题。与静态网格方法相比，本文的方法能够处理更广泛的运动边界问题，包括非刚体运动和变形域。
与现有 ALE 方法的比较：文献中已有一些针对可压缩 Navier-Stokes 方程的熵保守和熵稳定 ALE 公式 [10, 11]。然而，这些工作要么局限于周期性域，要么未对移动壁面边界条件的处理提供完整的分析和证明。本文填补了这一空白，提供了针对一般运动壁面的完整熵稳定性证明。此外，本文的方法不仅限于特定的离散化方案，而是基于 SBP 属性，因此可以推广到有限体积、有限元和通量重构等多种方法。
与低阶方法的比较：高阶方法在给定精度下通常比低阶方法需要更少的自由度，并具有更好的色散和耗散特性。然而，高阶方法的鲁棒性依赖于可证明稳定的公式设计，特别是在欠分辨流动、间断、尖锐梯度和复杂边界存在的情况下。本文通过提供熵稳定的高阶方案，解决了这一挑战，使其适用于大规模湍流和超音速模拟等复杂流动问题。
代价与局限：尽管本文的方法具有广泛的适用性，但其实现复杂度较高，需要精确处理几何守恒律和 SBP 算子的构造。此外，熵稳定性分析主要适用于空间半离散化，完全离散的熵估计还依赖于时间积分器的选择。作者指出，隐式格式可能贡献数值熵耗散，而显式时间积分方法可能会产生熵，即使与熵保守的空间离散化耦合。因此，完全离散的熵估计需要熵兼容的时间积分器，如适当的空间-时间离散化或松弛 Runge-Kutta 方法 [39]。这些限制在实际应用中可能需要额外的考虑和优化。

8. 线性化周期 Whitham–Boussinesq 系统的最佳多项式稳定

Optimal Polynomial Stabilization of the Linearized Periodic Whitham–Boussinesq System

作者: Roberto de A. Capistrano Filho, William Artiles Roqueta

类别: math.AP

发布日期: 2026-08-25

arXiv ID: 2608.25226v1

摘要翻译: 我们研究一维环面上线性化周期 Whitham–Boussinesq 系统的稳定性。我们在自然能量空间中建立了保守动力学与阻尼动力学的适定性，并描述了保守生成元的谱结构，其频率在高频处呈现 $|k|^{1/2}$ 阶次线性增长。随后，我们证明了阻尼半群的强稳定性，并获得了关于谱参数呈线性增长的高频预解式估计。在真正局部化的阻尼条件下，一族高频准模提供了相匹配的下界，表明该预解式增长是最佳的。依据 Borichev–Tomilov 定理 [5]，我们推导出半群在生成元定义域上的 $t^{-1}$ 衰减率，这对应于能量的 $t^{-2}$ 衰减率。在相同的局部化假设下，这些多项式衰减率是最佳的。

Comment: 38 pages. Comments are welcome

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

1. 研究对象和背景

本文研究一维环面 $\mathbb{T} = \mathbb{R}/(2\pi \mathbb{Z})$ 上线性化周期 Whitham–Boussinesq 系统的长期动力学行为, 特别是其在空间依赖阻尼作用下的稳定性与能量衰减率. 该模型源于有限深度不可压缩无旋流体表面波的双向传播问题. 在无量纲化后, 考虑如下非线性 Whitham–Boussinesq 系统:

\begin{cases} u_t + h_x + \alpha u u_x + \alpha \beta^2 (T_\beta u)(T_\beta u)_x = 0, \\ h_t - T_\beta u + \alpha \left[ \beta^2 T_\beta(h T_\beta u) + uh \right]_x = 0, \end{cases}

其中 $\alpha = a/h_0$ 衡量非线性强度, $\beta = h_0/\lambda$ 为水深与特征波长之比. 算子 $T_\beta$ 是一个非局部傅里叶乘子, 定义为

\widehat{T_\beta v}(k) = -i m_\beta(k) \hat{v}(k), \quad m_\beta(k) = \frac{\tanh(\beta k)}{\beta},

它保留了完整水波问题的精确线性色散关系. 本文关注的是该系统在零平衡点附近的线性化动力学, 即忽略非线性项后得到的保守系统:

\begin{cases} h_t - T_\beta u = 0, \\ u_t + h_x = 0. \end{cases}

利用恒等式 $\partial_x P_\beta = -T_\beta$ , 其中 $P_\beta$ 是符号为 $p_\beta(k) = \tanh(\beta k)/(\beta k)$ ( $k \neq 0$ ) 且 $p_\beta(0)=1$ 的正自伴傅里叶乘子, 上述系统可等价写为:

\begin{cases} h_t + P_\beta u_x = 0, \\ u_t + h_x = 0. \end{cases}

自然能量空间取为 $\mathcal{H} = L^2(\mathbb{T}) \times H^{-1/2}(\mathbb{T})$ , 配备内积

\langle (h,u), (\phi, \psi) \rangle_{\mathcal{H}} = \langle h, \phi \rangle_{L^2} + \langle P_\beta u, \psi \rangle_{H^{1/2}, H^{-1/2}},

对应能量泛函为 $E(t) = \frac{1}{2}\|h(t)\|_{L^2}^2 + \frac{1}{2}\langle P_\beta u(t), u(t)\rangle_{H^{1/2}, H^{-1/2}}$ . 保守生成元 $$A$$ 定义为

A \binom{h}{u} = \binom{T_\beta u}{-h_x}, \quad D(A) = H^{1/2}(\mathbb{T}) \times L^2(\mathbb{T}),

它在 $\mathcal{H}$ 中是斜自伴的, 因此生成一个酉 $$C_0$$ -群.

为了研究稳定性, 引入空间依赖的非负阻尼函数 $a(x), b(x) \in C^\infty(\mathbb{T}; \mathbb{R})$ , 满足 $a \not\equiv 0, b \not\equiv 0$ . 采用反馈控制律 $f = -ah, g = -bP_\beta u$ , 得到闭环阻尼系统:

\begin{cases} h_t - T_\beta u + a^2 h = 0, \\ u_t + h_x + b^2 P_\beta u = 0. \end{cases}

其抽象形式为 $$Z_t = A_d Z$$ , 其中 $$A_d = A - BB^*$$ , $B(f,g)^\top = (af, bg)^\top$ . 阻尼算子 $$D = BB^*$$ 使得能量耗散:

\frac{d}{dt} E(t) = -\|ah(t)\|_{L^2}^2 - \|bP_\beta u(t)\|_{L^2}^2 \le 0.

核心问题是确定能量 $$E(t)$$ 是否趋于零, 以及最优衰减率是什么. 由于色散关系的次线性增长特性, 经典的指数稳定性通常不成立, 而多项式稳定性成为关键研究对象. 高频色散频率满足 $\omega_k \sim \beta^{-1/2}|k|^{1/2}$ ( $|k| \to \infty$ ), 这种次线性增长导致经典 Ingham 型方法失效, 必须通过高频留数估计来确定衰减率.

2. 主要定理或主要结果

本文的主要结果由定理 1.1 (最优多项式稳定) 给出. 设 $$S_d(t)$$ 是由 $$A_d$$ 在 $\mathcal{H}$ 上生成的收缩半群, 假设局域化函数 $$a, b$$ 满足 $a, b \ge 0$ 且均不恒为零. 则:

谱性质: $i\mathbb{R} \subset \rho(A_d)$ , 即虚轴包含在 $$A_d$$ 的预解式定义域中.
半群衰减率: 存在常数 $$C > 0$$ 使得对于所有 $t \ge 1$ ,
$\|S_d(t) A_d^{-1}\|_{\mathcal{L}(\mathcal{H})} \le \frac{C}{t}.$
相应地, 对于任意初始数据 $Z_0 \in D(A_d)$ , 解 $$Z(t) = S_d(t)Z_0$$ 满足
$\|Z(t)\|_{\mathcal{H}} \le \frac{C}{t} \|A_d Z_0\|_{\mathcal{H}}, \quad t \ge 1.$
特别地,
$\|Z(t)\|_{\mathcal{H}} \le \frac{C}{t} \|Z_0\|_{D(A_d)}, \quad t \ge 1.$
从而能量衰减率为
$E(t) \le \frac{C}{t^2} \|Z_0\|_{D(A_d)}^2, \quad t \ge 1.$
最优性: 如果作用于第一分量 $$h$$ 的阻尼函数 $$a$$ 是真正局域化的, 即存在非空开集 $\mathcal{O} \Subset \mathbb{T} \setminus \text{supp}(a)$ , 则上述 $t^{-1}$ 半群衰减率在多项式尺度下是最优的. 相应地, $t^{-2}$ 的能量衰减率也是最优的, 即不存在更快的多项式衰减率.
留数估计: 高频留数满足线性增长上界:
$\|(i\tau I - A_d)^{-1}\|_{\mathcal{L}(\mathcal{H})} \lesssim |\tau|, \quad |\tau| \to \infty.$
在真正局域化假设下, 存在高频拟模序列 $\omega_N \to \infty$ 使得
$\|(i\omega_N I - A_d)^{-1}\|_{\mathcal{L}(\mathcal{H})} \gtrsim \omega_N,$
从而证明留数增长的线性阶数是尖锐的.
强稳定性: 由于 $i\mathbb{R} \subset \rho(A_d)$ , 根据 Arendt-Batty-Lyubich-Vũ 定理, 阻尼半群 $$S_d(t)$$ 是强稳定的, 即对于任意 $Z_0 \in \mathcal{H}$ , $\|S_d(t)Z_0\|_{\mathcal{H}} \to 0$ 当 $t \to \infty$ .

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

证明的核心框架分为四个主要步骤, 每一步都依赖于特定的分析工具和技术创新.

第一步: 功能分析与适定性建立. 首先建立保守算子 $$A$$ 的斜自伴性和阻尼算子 $$A_d$$ 的最大耗散性. 通过 Parseval 恒等式和傅里叶系数计算, 验证了 $$A$$ 在能量内积下的斜对称性. 关键观察是 $P_\beta$ 作为阶数为 $$-1$$ 的傅里叶乘子, 定义了 $H^s(\mathbb{T}) \to H^{s+1}(\mathbb{T})$ 的同构, 这使得能量范数与标准乘积范数等价. 阻尼算子 $$BB^*$$ 是有界的耗散扰动, 因此 $$A_d$$ 生成收缩半群. 这一步为后续的谱分析和留数估计奠定了严格的功能分析基础.

第二步: 保守算子的谱结构分析. 通过傅里叶分解, 每个非零空间模态满足振荡方程 $\partial_t^2 \hat{h}(k,t) + \omega_k^2 \hat{h}(k,t) = 0$ , 其中频率为

\omega_k = \sqrt{\frac{k \tanh(\beta k)}{\beta}}.

高频渐近行为为 $\omega_k \sim \beta^{-1/2}|k|^{1/2}$ . 这一次线性增长是本文技术难点的根源: 频率序列不满足均匀间隙条件, 因此经典的 Ingham 不等式或非调和傅里叶级数方法无法直接应用. 作者明确指出, 这迫使定量稳定性的研究转向高频留数估计方法. 此外, 证明了保守算子的留数是紧的, 谱由离散的本征值组成.

第三步: 高频留数分析 (核心技术创新). 这是论文最技术密集的部分, 分为上界和下界两部分.

上界证明: 关键技巧是通过变量替换 $v = P_\beta^{1/2} u$ 将二阶耦合系统约化为标量算子. 令 $L_\beta = T_\beta \partial_x$ , 则标量算子满足
$\beta L_\beta = |D| + R_\beta,$
其中 $$|D|$$ 是分数阶导数算子 (符号为 $$|k|$$ ), 而 $R_\beta$ 是一个光滑傅里叶乘子. 这一表示允许作者借用 Alazard, Marzuola 和 Wang [1] 针对分数阶波动方程 $|D| - \tau^2 - i\tau \chi$ 的高频驻波留数估计. 具体而言, 他们处理了形如 $|D| - \tau^2 - i\tau a(x)$ 的算子, 并证明了留数的线性增长上界. 在本文中, $R_\beta$ 作为一个紧致微扰出现, 通过紧性论证可以吸收其影响, 从而将 Alazard-Marzuola-Wang 的结果推广到当前的非局部 Whitham–Boussinesq 设置. 最终得到 $\|(i\tau I - A_d)^{-1}\| \lesssim |\tau|$ .
下界构造 (拟模): 为了证明上界的尖锐性, 作者在真正局域化假设下构造了一族高频拟模. 这些拟模集中在未受阻尼的区域 $\mathcal{O}$ 中. 通过精心选择频率 $\omega_N$ 和振幅分布, 使得残余项足够小, 从而得到
$\|(i\omega_N I - A_d)^{-1}\| \gtrsim \omega_N.$
这证明了留数不能以比线性更慢的速度增长, 从而确立了线性增长的最优性.

第四步: Borichev-Tomilov 定理的应用. 最后, 利用 Borichev 和 Tomilov [5] 的定理, 将留数的多项式增长转化为半群的衰减率. 该定理指出, 如果留数满足 $\|(i\tau I - A)^{-1}\| = O(|\tau|^m)$ , 则半群在生成元定义域上的衰减率为 $O(t^{-1/m})$ . 在这里 $$m=1$$ , 因此得到 $\|S_d(t) A_d^{-1}\| = O(t^{-1})$ , 进而能量衰减为 $O(t^{-2})$ . 结合拟模的下界, 证明了该多项式速率在多项式尺度下是最优的, 即无法获得 $o(t^{-1})$ 或 $o(t^{-2})$ 的衰减.

4. 与之前工作的比较

本文将线性化 Whitham–Boussinesq 系统的稳定性问题置于更广泛的偏微分方程稳定性理论背景下, 并与几类相关工作进行了比较.

1. 与经典阻尼波动方程稳定性的比较: 传统上, 阻尼波动方程的多项式稳定性研究主要集中在几何控制条件失败的情形. Anantharaman 和 Léautaud [2] 研究了环面上缺乏均匀几何控制时的尖锐多项式衰减率, 发现衰减由高频留数估计、捕获轨迹以及阻尼系数在零集附近的行为决定. 尽管机制不同, 但本文采用了类似的基于留数的策略: 高频拟模提供留数下界以识别指数稳定性的障碍, 而多项式留数估计确定相应的衰减率. 然而, 本文的障碍并非来自经典测地线流或几何控制条件的失败, 而是源于非局部系统的次线性色散关系. 这意味着即使阻尼在整个环面上非平凡, 由于色散的低频增强效应, 指数稳定性依然不可能, 这与经典双曲方程的行为有本质区别.

2. 与分数阶波动方程及非局部方程的比较: 本文的技术核心依赖于 Alazard, Marzuola 和 Wang [1] 关于阻尼分数阶波动方程的高频驻波估计. 他们的结果针对算子 $|D| - \tau^2 - i\tau \chi$ 提供了留数的线性增长上界. 本文的创新在于展示了如何将这一结果适配到 Whitham–Boussinesq 系统. 通过变量替换, 标量算子被表示为 $$|D|$$ 加上一个光滑乘子 $R_\beta$ . 这表明, 尽管 Whitham–Boussinesq 系统是非局部的, 但其高频行为本质上类似于分数阶波动方程, 差异仅体现在一个紧致微扰中. 这一观察使得作者能够直接借用现有的分数阶工具, 避免了从头推导复杂的非局部留数估计.

3. 与 Whitham–Boussinesq 系统其他分析工作的比较: 近年来, Whitham–Boussinesq 类型系统的许多分析方面已被研究. Dinvay [9], Dinvay-Selberg-Tesfahun [11], Emerald [13], Paulsen [19] 以及 Pilod et al. [21] 分别建立了局部适定性、低正则性全局存在性、长时间适定性以及解的生命期估计. 这些工作主要关注非线性动力学和适定性, 而非线性化系统的稳定性. 此外, Hur-Pandey [15] 和 Creedon-Deconinck-Trichtchenko [8] 研究了小振幅周期行波的谱稳定性, 发现了纯虚本征值的碰撞和高频谱不稳定性. 虽然这些工作与本文都涉及高频谱结构, 但它们的对象是非平凡周期行波周围的线性化算子, 而本文研究的是零平衡点周围的线性化动力学. 本文的贡献在于首次系统地确定了零平衡点处阻尼线性化系统的最优多项式衰减率, 揭示了次线性色散对长期动力学的根本限制.

arXiv ID: 2608.25458v1

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

摘要: Let $$H_0$$ be the standard discrete Laplacian on $\mathbb Z^d$ , $d\geq4$ , let $R_0(z)=(H_0-z)^{-1}$ , and let $$p'$$ denote the Hölder conjugate of $$p$$ , with $p'=\infty$ when $$p=1$$ . We establish uniform diagonal resolvent estimates from $\ell^p(\mathbb Z^d)$ to $\ell^{p'}(\mathbb Z^d)$ by proving local Fourier-decay bounds for surface measure on the Fermi surfaces of the lattice dispersion relation. At a regular point, the decay is governed by the number of coordinate directions in which the quadratic term vanishes: after a local change of variables, these directions produce cubic one-dimensional phases, while the remaining directions are quadratic. For $d\geq5$ , we prove the sharp global range

1\leq p\leq\frac{2(d+2)}{d+5}.

Away from the threshold energies, this range remains sharp in odd dimensions, whereas for even $d\geq6$ the sharp range improves to

1\leq p\leq\frac{2(2d+5)}{2d+11}.

In dimension four, we obtain uniform estimates for $1\leq p<4/3$ and an $\ell^{4/3}$ -to- $\ell^4$ endpoint estimate with a square-root logarithmic loss. Matching anisotropic Knapp-type constructions give the corresponding necessary conditions. We also derive Birman–Schwinger and Kato smoothing bounds, complex-eigenvalue estimates, and improved thin spectral projection estimates.

8. Minkowski 双曲面剖分上的三维分部求和（SBP）格式

A 3D Summation-by-Parts scheme on a Hyperboloidal Foliation of Minkowski

作者: Shalabh Gautam

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

发布日期: 2026-08-26

arXiv ID: 2608.25363v1

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

摘要: This paper summarises our previous work on a fully $$3$$ D Summation-by-Parts scheme, derived for a class of linear wave equations on hyperboloidal slices on a fixed Minkowski background. The scheme is derived in spherical polar coordinates, and allows having grid points at the origin and on the $$z$$ -axis, despite coordinate singularities, and at infinity, by introducing compactification followed by rescaling, and is proved to be stable. Reducing it to the standard Cauchy problem, or to finite spacelike slices with an outer boundary, will follow a similarly. Second-order accurate finite-difference methods are used to implement this scheme numerically, but higher-order finite-difference or spectral methods could also be used. Kreiss-Oliger dissipation operators are generalized to curvilinear coordinates and are defined everywhere in the domain, including at the boundary points, such that they satisfy the dissipative property in the energy norms. We also propose new norm convergence tests that include all the grid points at all resolutions and produce more accurate results. Promising results are obtained, giving hope for application to fully nonlinear systems, like the Einstein Field Equations, and extracting the resulting gravitational waves free of systematic errors or gauge ambiguities.

9. 具有Soler型非线性的非线性Dirac方程的周期解

Periodic Solutions for a Nonlinear Dirac Equation with Soler-Type Nonlinearity

作者: Fuping Zhang

类别: math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25320v1

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

摘要: We study nonlinear Dirac equations on the three-dimensional torus with subcritical Soler-type nonlinearities. The variational problem is strongly indefinite, while the non-coercive Soler potential prevents direct compactness. We recover the Palais–Smale condition by introducing a small coercive perturbation. Refined spectral estimates and a linking argument yield perturbed solutions with uniform energy bounds. Passing to the limit as the perturbation vanishes gives a nontrivial continuously differentiable periodic solution. This existence result holds for every spectral parameter in an explicit interval below the mass.

10. 杨-米尔斯-希格斯-薛定谔流

Yang-Mills-Higgs-Schrödinger flow

作者: Bo Chen, Chong Song

类别: math.DG, math.AP

发布日期: 2026-08-26

arXiv ID: 2608.25297v1

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

摘要: In this paper, we initiate the study of the Yang-Mills-Higgs-Schrödinger(YMHS) flow, i.e. the Hamiltonian flow of the Yang-Mills-Higgs functional defined on a symplectic fiber bundle. The YMHS flow provides a natural gauge-theoretic extension of the classical Schrödinger flow within the framework of symplectic reduction theory, and generalizes the Chern-Simons-Schrödinger equations to a non-Abelian gauge and non-linear fiber setting. We study its geometric structures and establish the local well-posedness of the corresponding Cauchy problem on compact Riemann surfaces.

11. $\Z^3$ 上的精确离散 Hardy 不等式

The sharp discrete Hardy inequality on $\Z^3$

作者: Natanael Alpay

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

发布日期: 2026-08-26

arXiv ID: 2608.25262v1

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

摘要: We determine the sharp constant in the nearest-neighbor Hardy inequality on $\Z^3$ with the Euclidean inverse-square weight. For every finitely supported function $u:\Z^3\to\C$ , we prove

\sum_{x\in\Z^3}\sum_{j=1}^3 |u(x+e_j)-u(x)|^2 \geq \frac14\sum_{x\in\Z^3\setminus\{0\}} \frac{|u(x)|^2}{|x|^2}.

The coefficient $$1/4$$ is sharp, and equality is not attained by a nonzero finitely supported function. The proof uses an explicit reciprocal edge field and an edgewise completion of squares, together with a concavity argument for the associated vertex weight.

12. 关于分数阶 $$1$$ -拉普拉斯演化方程

On fractional $$1$$ -Laplacian evolution equation

作者: Rakesh Arora, Nitin Kumar Maurya

类别: math.AP

发布日期: 2026-08-25

arXiv ID: 2608.25105v1

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

摘要: We study a nonlocal evolution problem driven by the fractional $$1$$ -Laplacian with a Carathéodory nonlinearity satisfying subcritical growth conditions. We develop a potential well framework to investigate the existence and qualitative behavior of solutions at different initial energy levels. By combining a modified potential well method with Galerkin approximations, we establish the global existence of weak and strong solutions under appropriate conditions on the initial energy and the sign of the associated Nehari functional. To address the singular nature of the fractional $$1$$ -Laplacian, we approximate the problem by a family of fractional $$p$$ -Laplacian equations, derive estimates uniform with respect to $$p>1$$ , and pass to the limit as $p\to 1^{+}$ . Furthermore, in the low-dimensional regime $$N<2s$$ , we employ a subdifferential approach to establish the local existence of strong solutions and investigate their qualitative behavior.

13. 振荡子半群的 Schrödinger 刻画

A Schrödinger characterization of the oscillator semigroup

作者: Gianluca Giacchi

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

发布日期: 2026-08-25

arXiv ID: 2608.25082v1

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

摘要: The oscillator (or metaplectic) semigroup is classically defined as the two-fold cover of the semigroup of positive complex symplectic matrices. Although geometrically precise, this definition does not provide an intrinsic operator-theoretic characterization of the bounded operators corresponding to positive symplectic matrices. This is in contrast with the real metaplectic group, whose relation with the symplectic group can be expressed directly in terms of the Schrödinger representation of the Heisenberg group through its intertwining property. In the complex setting, such a relation appears in the existing literature mainly in infinitesimal form or on suitable classes of Gaussian functions. In this work we prove that the (complexified) Schrödinger intertwining relation holds for every function in $L^2(\mathbb{R}^d)$ . The main point is the converse statement: if an arbitrary complex symplectic matrix $$S$$ admits a nonzero bounded operator on $L^2(\mathbb{R}^d)$ satisfying this intertwining relation, then $$S$$ is necessarily positive and the operator coincides, up to a nonzero scalar, with the corresponding element of the oscillator semigroup. Consequently, positive complex symplectic matrices are exactly those admitting a nonzero bounded Schrödinger intertwiner, and the associated intertwining space is one-dimensional. This provides a Schrödinger-representation characterization of the oscillator semigroup, parallel to the classical one for the real metaplectic group.

14. 超二次Hamilton-Jacobi方程的最大正则性：Lions猜想中的端点情形

Maximal Regularity of Superquadratic Hamilton-Jacobi Equations: The Endpoint Case in Lions’ Conjecture

作者: Fanze Kong

类别: math.AP

发布日期: 2026-08-25

arXiv ID: 2608.25066v1

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

摘要: A celebrated conjecture of P.-L. Lions concerns maximal regularity for viscous Hamilton-Jacobi equations. In this paper, we study the endpoint case. We consider normalized strong solutions of $$-Δu+|Du|^γ=f$$ on $\mathbb T^d$ , where $d\geq 2$ , $$γ>2$$ , and $f\in L^{q_c}(\mathbb T^d)$ with $$q_c=d(γ-1)/γ$$ . At this critical exponent, the main difficulty is possible concentration under the critical scaling. Assuming that the source terms form a uniformly equi-integrable subset of $L^{q_c}$ , we rule out this concentration and prove maximal $L^{q_c}$ regularity for strong solutions. The proof combines a two-stage blow-up argument with a Liouville rigidity theorem. A critical local energy estimate excludes concentration of the gradient energy in $L^{γq_c}$ at the first blow-up scale, while small-drift regularity upgrades weak compactness of the blow-up sequence to strong local compactness.

15. Transformer的均场理论：耦合数据-参数动力学的适定性与训练的全局收敛

A Mean-Field Theory of Transformers: Well-Posedness of the Coupled Data–Parameter Dynamics and Global Convergence of Training

作者: Michael Herty, Hailiang Liu

类别: math.AP

发布日期: 2026-08-25

arXiv ID: 2608.25055v1

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

摘要: We develop a rigorous mean-field theory for transformer networks that captures two large-scale limits inherent in the architecture: the number of tokens $N\to\infty$ in the input sequence and the number of attention heads $H\to\infty$ in each layer. It further considers an infinite number of layers leading to a time-continuous formulation. The resulting framework couples two interacting mean-field objects: a token distribution $μ_t\in\mathcal P(\mathbb R^d)$ , which evolves through network depth $t\in[0,T]$ according to a McKean–Vlasov transport equation, and the attention-parameter distribution $ρ_s\in\mathcal P(Θ)$ , that evolves through training time $s\ge0$ according to a Wasserstein gradient flow of the empirical risk, with optional entropic or Tikhonov regularization. We establish a comprehensive analytical foundation for the system coupling the transformer and the training dynamics; in particular, we establish global well-posedness of the resulting nonlinear Fokker–Planck system describing the coupled mean-field and training dynamics. Beyond well-posedness, we connect the mean-field formulation to optimization. For shallow, single-layer attention models, we prove exponential convergence to the entropy-regularized global optimum under a log-Sobolev condition. For genuinely deep, compositional transformers, we establish local linear convergence under a Neural Tangent Kernel non-degeneracy condition.

16. 无 Hessian 高分辨率蒙特卡洛采样的改进分析

Improved Analysis for Hessian-free High-resolution Monte Carlo Sampling

作者: Wujun Lv, Xiaoyu Wang, Yingli Wang, Lingjiong Zhu

类别: stat.ML, cs.LG, math.AP, math.PR

发布日期: 2026-08-25

arXiv ID: 2608.25052v1

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

摘要: Hessian-free high-resolution (HFHR) dynamics augments underdamped Langevin dynamics (ULD) with reversible position diffusion for sampling problems that arise in machine learning. We establish an explicit quantitative contraction rate for HFHR dynamics under a position Poincaré inequality, weighted Hessian and Laplacian bounds, and a compact Sobolev embedding, where the potential function is not necessarily convex. An adapted time-augmented Poincaré inequality yields an explicit rate that improves upon the contraction rate of the underdamped Langevin dynamics. We also give a weak-solution construction and a self-contained spectral proof of the divergence lemma underlying the argument. For HFHR Monte Carlo (HFHRMC) algorithm, which is based on a discretization scheme of HFHR dynamics, we use a path-space Girsanov argument to obtain a non-asymptotic convergence bound and an explicit iteration complexity in total variation distance. The bounds hold for every $α\geq0$ and $$γ>0$$ and remain regular at the ULD endpoint. Optimizing the iteration complexity bound yields a positive, accuracy-dependent position-diffusion parameter at finite accuracy, while its leading high-accuracy order coincides with that of the optimized ULD endpoint. Our iteration complexity bound improves upon the existing work on HFHR algorithms. Numerical experiments including Bayesian learning problems on real data are provided to illustrate the effect of positive $$α$$ and its benefit.

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

1. 涡丝动力学与涡环运动的再审视#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

2. 半空间中定常 Navier–Stokes 方程的强迫自相似解#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

3. 二维欧拉方程无穷通道中定常解的灵活性与刚性#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

4. 球对称下具有退化黏度的可压缩 Navier-Stokes 方程的膨胀解：全局存在性与无黏极限#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

5. 带体-表面相互作用的不压缩两相流扩散界面模型的 Voigt 正则化的全局准强解#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

6. 三维 $ rac{1}{2}|u|^2+p$ 的几何细化#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

7. 可压缩 Navier-Stokes 方程 ALE 格式中熵稳定的移动壁面边界条件#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

8. 线性化周期 Whitham–Boussinesq 系统的最佳多项式稳定#

1. 研究对象和背景#

2. 主要定理或主要结果#

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

4. 与之前工作的比较#

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

1. Heisenberg 群上临界拟线性 Schrödinger 方程：存在性与非存在性#

摘要翻译#

2. $λ$-耗散解的存在性、唯一性与长时间行为：二分量 Hunter-Saxton 系统#

摘要翻译#

3. 二维应力扩散型 Oldroyd-B 模型在蠕动流区间的初边值问题强解的全局适定性#

摘要翻译#

4. 去聚焦 mKdV 方程的几乎周期解#

摘要翻译#

5. Kirchhoff–Pohozaev 波动方程解的存在性与等连续性#

摘要翻译#

不相关论文（基本信息）#

1. 关于散度方程及相关函数不等式的统一方法#

2. 亚极端 Kerr-de Sitter 时空上标量波动方程的模式稳定性#

3. Kerr时空中的双模与Whiting变换：模稳定性再探#

4. 完全渐近锥形流形上的 $G_2$-拉普拉斯流#

5. 具有自由边界和径向对称性的高维单稳态反应扩散方程#

6. $σ_k$ 曲率方程在 $k rac{}{}geq n/2$ 时的整体曲率估计#

7. 高维离散 Schrödinger 算子的均匀预解式估计#

8. Minkowski 双曲面剖分上的三维分部求和（SBP）格式#

9. 具有Soler型非线性的非线性Dirac方程的周期解#

10. 杨-米尔斯-希格斯-薛定谔流#

11. $\Z^3$ 上的精确离散 Hardy 不等式#

12. 关于分数阶 $1$-拉普拉斯演化方程#

13. 振荡子半群的 Schrödinger 刻画#

14. 超二次Hamilton-Jacobi方程的最大正则性：Lions猜想中的端点情形#

15. Transformer的均场理论：耦合数据-参数动力学的适定性与训练的全局收敛#

16. 无 Hessian 高分辨率蒙特卡洛采样的改进分析#

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

1. 涡丝动力学与涡环运动的再审视

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

2. 半空间中定常 Navier–Stokes 方程的强迫自相似解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

3. 二维欧拉方程无穷通道中定常解的灵活性与刚性

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

4. 球对称下具有退化黏度的可压缩 Navier-Stokes 方程的膨胀解：全局存在性与无黏极限

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

5. 带体-表面相互作用的不压缩两相流扩散界面模型的 Voigt 正则化的全局准强解

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

6. 三维 $rac{1}{2}|u|^2+p$ 的几何细化

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

7. 可压缩 Navier-Stokes 方程 ALE 格式中熵稳定的移动壁面边界条件

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

8. 线性化周期 Whitham–Boussinesq 系统的最佳多项式稳定

1. 研究对象和背景

2. 主要定理或主要结果

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

4. 与之前工作的比较

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

1. Heisenberg 群上临界拟线性 Schrödinger 方程：存在性与非存在性

摘要翻译

2. $$λ$$ -耗散解的存在性、唯一性与长时间行为：二分量 Hunter-Saxton 系统

摘要翻译

3. 二维应力扩散型 Oldroyd-B 模型在蠕动流区间的初边值问题强解的全局适定性

摘要翻译

4. 去聚焦 mKdV 方程的几乎周期解

摘要翻译

5. Kirchhoff–Pohozaev 波动方程解的存在性与等连续性

摘要翻译

不相关论文（基本信息）

1. 关于散度方程及相关函数不等式的统一方法

2. 亚极端 Kerr-de Sitter 时空上标量波动方程的模式稳定性

3. Kerr时空中的双模与Whiting变换：模稳定性再探

4. 完全渐近锥形流形上的 $$G_2$$ -拉普拉斯流

5. 具有自由边界和径向对称性的高维单稳态反应扩散方程

6. $$σ_k$$ 曲率方程在 $$k rac{}{}geq n/2$$ 时的整体曲率估计

7. 高维离散 Schrödinger 算子的均匀预解式估计

8. Minkowski 双曲面剖分上的三维分部求和（SBP）格式

9. 具有Soler型非线性的非线性Dirac方程的周期解

10. 杨-米尔斯-希格斯-薛定谔流

11. $\Z^3$ 上的精确离散 Hardy 不等式

12. 关于分数阶 $$1$$ -拉普拉斯演化方程

13. 振荡子半群的 Schrödinger 刻画

14. 超二次Hamilton-Jacobi方程的最大正则性：Lions猜想中的端点情形

15. Transformer的均场理论：耦合数据-参数动力学的适定性与训练的全局收敛

16. 无 Hessian 高分辨率蒙特卡洛采样的改进分析