生成时间: 2026年02月25日 04:36:21

重点关注论文数量: 5

了解领域论文数量: 6

不相关论文数量: 17

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

1. 等涡度流中平坦底床上周期性行波下方的压力

Pressure beneath a periodic travelling water-wave in constant-vorticity flow over a flat bed

作者: Adrian Constantin, Nicolas Gindrier, Otmar Scherzer

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

发布日期: 2026-02-24

arXiv ID: 2602.21077v1

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

详细分析

中文标题: 等涡度流中平坦底床上周期性行波下方的压力

研究对象与背景

本文研究二维周期性行波在常涡度、水平平坦底床上的流体压力特性。研究背景源于海洋工程中对水下极端压力估计的迫切需求，底部压力传感器常用于通过静水压力分布估算波高，但非静水修正可能显著，这需依赖对波浪下方精确压力分布的理解。绝大多数海洋表面波浪（如涌浪）源于风场，远离生成区后趋于周期性，传播距离可达数百公里。在无洋流情形下，这类波浪的流动近似无旋；然而，波浪常与变深厚流（如潮汐流、 Lyapunov 流）相互作用。为兼顾物理真实性与数学可处理性，本文考虑 常涡度流（constant vorticity）——即涡量 $\Omega$ 为常数，这是风生波与剪切流相互作用的常用模型，且在波长远大于水深时，非零平均涡度比其具体分布更重要。底床设定为水平平坦，适用于运河、深海平原或大陆架等场景。假设流动无粘、二维、周期性，采用 线性理论（wave amplitude $a$ 远小于水深 $d$，即 $\varepsilon = a/d \ll 1$），以聚焦小振幅波对背景流的扰动。

控制方程组为二维不可压缩 Euler 方程，重力为唯一恢复力。设波以恒定速度 $c>0$ 沿 $X$ 方向传播，水深 $d$，平坦底面位于 $Y=-d$。引入常涡度 $\Omega$，背景流速为 $(\Omega(Y+d), 0)$，故总速度场为 $(U(X,Y,T)+\Omega(Y+d), V(X,Y,T))$。在流体域 $-d \leq Y \leq \eta(X,T)$ 中（$\eta$ 为自由表面），方程为：

$$ \begin{cases} U_X + V_Y = 0, & \text{(质量守恒)} \\[6pt] U_T + [\Omega(Y+d)+U]U_X + V[\Omega + U_Y] = -\frac{1}{\rho}P_X, & \\[6pt] V_T + [\Omega(Y+d)+U]V_X + V V_Y = -\frac{1}{\rho}P_Y - g, \end{cases} $$

其中 $P$ 为总压力，$\rho$ 为密度，$g$ 为重力加速度，涡量 $\Omega = U_Y - V_X$ 为常数。边界条件：底部无穿透 $V=0$ 于 $Y=-d$；自由表面 $Y=\eta$ 处运动学条件 $V = \eta_T + [\Omega(Y+d)+U]\eta_X$ 和动力学条件 $P=P_{\text{atm}}$（大气压常数）。定义 动压力 $\hat{p}$ 使

$$ P(X,Y,T) = P_{\text{atm}} - \rho g Y + \rho g d \hat{p}(x,y,t), \quad \text{其中 } \hat{p} \text{ 无量纲化后为 } p(x,y,t)。 $$

引入无量纲变量：水平尺度 $L$（波长），垂直尺度 $d$，波幅 $a$。令 $\delta = d/L$（浅水参数），$\varepsilon = a/d$（幅值参数）。缩放：

$$ X = L x,\quad Y = d y,\quad T = \frac{L}{c_0}t,\quad c_0 = \sqrt{gd}, $$

$$ U = c_0 u,\quad V = \frac{c_0 d}{L} v,\quad \eta = a h,\quad \omega = \Omega \sqrt{\frac{d}{g}} $$

得无量纲方程组（在 $-1<y<\varepsilon h$）：

$$ \begin{aligned} u_x + v_y &= 0,\\ u_t + [u + \omega(1+y)]u_x + v[\omega + u_y] &= -p_x,\\ \delta^2 v_t + [u + \omega(1+y)]v_x + v v_y &= -\delta^2 p_y,\\ v &= 0 \quad\text{on } y=-1,\\ v &= \varepsilon h_t + [u + \omega(1+\varepsilon h)]h_x \quad\text{on } y=\varepsilon h,\\ p &= \varepsilon h \quad\text{on } y=\varepsilon h,\\ u_y &= \delta^2 v_x. \end{aligned} $$

线性化取 $\varepsilon \to 0$（波幅趋于零），同时保留 $\delta$（波长与水深比）。设 $u=\varepsilon \hat{u}$，$v=\varepsilon \hat{v}$，$p=\varepsilon \hat{p}$，得线性系统（自由表面退到 $y=0$）：

$$ \begin{cases} \hat{u}_x + \hat{v}_y = 0, & -1<y<0 \\[6pt] \hat{u}_t + \omega(1+y)\hat{u}_x + \omega \hat{v} = -\hat{p}_x, & \\[6pt] \delta^2 \hat{v}_t + \omega(1+y)\hat{v}_x = -\delta^2 \hat{p}_y, & \\[6pt] \hat{v}=0 \text{ on } y=-1, & \\[6pt] \hat{u}\big|_{y=0} = \delta^2 \hat{v}_y\big|_{y=0}, & \\[6pt] \hat{p}\big|_{y=0} = h, & \\[6pt] h_t = \hat{v}\big|_{y=0}. & \end{cases} $$

引入 行波 Ansatz：$\hat{u},\hat{v},\hat{p},h$ 依赖于 $\xi = x - \hat{c}t$，其中 $\hat{c}>0$ 为无量纲波速（物理波速 $c = \hat{c}\sqrt{gd}$）。系统化为常微分方程组（$y$ 为自变量）：

$$ \begin{cases} \hat{u}_\xi + \hat{v}_y = 0, \\ [\hat{c} - \omega(1+y)]\hat{u}_\xi = \hat{p}_\xi, \\ \delta^2 \hat{v}_\xi = \hat{p}_y, \\ \hat{v}=0 \text{ on } y=-1, \\ \hat{u}\big|_{y=0} = \delta^2 \hat{v}_\xi\big|_{y=0}, \\ \hat{p}\big|_{y=0} = h, \\ h_\xi = -\hat{v}\big|_{y=0}. \end{cases} $$

利用 $x$ 方向周期性（周期 1），进行 Fourier 模式展开。设

$$ \hat{v}(x,y) = \sum_{k\in\mathbb{Z}} \alpha_k(y) e^{2\pi i k x},\quad -1\le y\le 0. $$

由 $\hat{u}_\xi + \hat{v}_y=0$ 和周期性得 $\alpha_k''(y) - 4\pi^2 k^2 \delta^2 \alpha_k(y)=0$。边界 $\hat{v}(x,-1)=0$ 给出

$$ \alpha_k(y) = A_k e^{2\pi k \delta (y+1)} + B_k e^{-2\pi k \delta (y+1)},\quad A_k + B_k=0. $$

结合实值条件，最终得

$$ \hat{v}(x,y) = \sum_{k\ge1} \gamma_k \sinh[2\pi k \delta (y+1)] \cos(2\pi k x), $$
$$ \hat{u}(x,y) = -\delta \sum_{k\ge1} \gamma_k \cosh[2\pi k \delta (y+1)] \sin(2\pi k x), $$

其中 $\gamma_k$ 为实系数。定义 流函数 $\hat{\psi}$ 使 $\hat{u}=\hat{\psi}_y$, $\hat{v}=-\hat{\psi}_x$，则

$$ \hat{\psi}(x,y) = -\sum_{k\ge1} \frac{\gamma_k}{2\pi k} \sinh[2\pi k \delta (y+1)] \sin(2\pi k x). $$

将 $\hat{u},\hat{\psi}$ 代入第二、三方程积分得动压力：

$$ \hat{p}(x,y) = [\hat{c}-\omega(1+y)]\hat{u}(x,y) + \omega \hat{\psi}(x,y) + \beta, $$

$\beta$ 为常数。利用 $\hat{p}(x,0)=h(x)$ 及平均海平面为零的条件 $\int_0^1 h(x)\,dx=0$ 推出 $\beta=0$。因此

$$ \hat{p}(x,y) = [\hat{c}-\omega(1+y)]\hat{u}(x,y) + \omega \hat{\psi}(x,y). $$

波廓 $h(x)$ 由 $\hat{p}(x,0)=h(x)$ 和 $h_\xi = -\hat{v}(x,0)$ 联立，导出色散关系（对每个 $k\ge1$）：

$$ \gamma_k \left\{ [\hat{c}-\omega]^2 + \frac{\omega \tanh(2\pi k \delta)}{2\pi k \delta} [\hat{c}-\omega] - \frac{\tanh(2\pi k \delta)}{2\pi k \delta} \right\} = 0. $$

函数 $S_\pm(s)= -\frac{\omega s}{2} \pm \sqrt{\frac{\omega^2 s^2}{4}+s}$ 在 $s>0$ 下单调性表明，至多一个 $k$ 使 $\gamma_k \ne 0$。故考虑 主周期波（单位周期），取 $k=1$，$\gamma_k=0$ for $k\ge2$，得

$$ h(x) = A \sin(2\pi x),\quad A>0, $$

$$ \hat{p}(x,y) = -\delta \gamma_1 \left\{ [\hat{c}-\omega(1+y)] \cosh[2\pi \delta (y+1)] + \frac{\omega}{2\pi} \sinh[2\pi \delta (y+1)] \right\} \sin(2\pi x). $$

波速 $\hat{c}$ 由色散关系（取负号以允许流反转）确定：

$$ \hat{c}-\omega = -\frac{\omega \tanh(2\pi \delta)}{4\pi \delta} \pm \sqrt{\left(\frac{\omega \tanh(2\pi \delta)}{4\pi \delta}\right)^2 + \frac{\tanh(2\pi \delta)}{2\pi \delta}}. $$

论文在区域 $\{(x,y): -1\le y\le 0,\ -1/4 \le x \le 1/4\}$（半周期矩形）中研究动压力 $\hat{p}$ 和总压力 $P = P_{\text{atm}} - \rho g d y + \rho g d \hat{p}$ 的极值。

主要定理与结论

论文未明确编号定理，但核心结果集中在第 3、4 节，可归纳为：

结果 1（动压力极值分布）：对于常涡度流中的主周期行波，动压力 $\hat{p}(x,y)$ 的最大值和最小值分别出现在 $x=1/4$（波峰）和 $x=-1/4$（波谷）的某些深度 $y$ 上。当涡度 $\omega=0$（无旋）时，函数 $Q(y) = \frac{A}{\cosh(2\pi \delta)} \cosh[2\pi \delta (y+1)]$ 在 $[-1,0]$ 单调递增，故最大值在 $(x,y)=(1/4,0)$（表面波峰），最小值在 $(-1/4,0)$（表面波谷）。当 $\omega \ne 0$ 时，$Q(y)$ 的表达式变复杂，可能出现非单调或变号，导致全局极值点 $y$ 坐标改变。具体地：

  • 若 $\omega>0$（顺流）且波长较长（$\delta$ 较大），或 $\omega<0$（逆流）且波长足够长，流反转可能出现（$\hat{c}-\omega<0$），临界层位于 $y_c = (\hat{c}-\omega)/\omega \in (-1,0)$。此时 $Q(y)$ 可能在 $y=-1$（平坦底床）或 $y=y_c$ 达到极值。
  • 极值点可能“along the flat bed”（即 $y=-1$）或“along the critical level”（即 $y=y_c$），而 $x$ 仍为 $\pm 1/4$，但最大值与最小值对应的 $x$ 可能因 $Q(y)$ 变号而交换（例如 $Q(y)<0$ 时，最大值在波谷 $x=-1/4$，最小值在波峰 $x=1/4$）。

结果 2（总压力极值）：总压力 $P$ 包含静水项 $-\rho g d y$，故随深度 $y$ 减小（$y$ 更负）而增加，该趋势不受涡度影响。但在固定深度 $y=y_0$，$P$ 关于 $x$ 的极值位置仍为 $x=\pm 1/4$，然而极值的大小顺序可能随 $\omega$ 改变。更重要的是，全局极值的 $y$ 坐标可显著偏离表面 $y=0$，例如出现在底床 $y=-1$，这对底部压力传感器的安装和解读有实际意义。

结果 3（流反转与临界层）：当 $\hat{c}-\omega < 0$ 且 $\omega>0$ 时，流反转发生，水平速度 $ \hat{c}-\omega(1+y)$ 在 $y \in (-1,0)$ 内变号，临界层 $y_c = (\hat{c}-\omega)/\omega$ 是速度为零的深度。极值位置与 $y_c$ 密切相关，尤其当 $Q(y)$ 在 $y_c$ 处有奇性或极值。

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

a) 核心分析学工具

  1. 线性化与行波约化：将 Euler 方程组通过小参数 $\varepsilon$ 线性化，并施加行波 Ansatz $\xi=x-\hat{c}t$，将原初边值问题转化为常微分方程组（仅依赖 $y$）。这极大简化了分析，使问题可显式求解。
  2. Fourier 模式分离：利用 $x$ 方向周期性，将速度场展开为 Fourier 级数。由连续性方程和底部无穿透条件，导出各模态 $\alpha_k(y)$ 的二阶线性常微分方程，解为双曲函数。通过实值条件和平均流假设（$\gamma_0=0$），将解表示为余弦级数形式。
  3. 色散关系推导：由自由表面条件 $\hat{p}(x,0)=h(x)$ 和 $h_\xi=-\hat{v}(x,0)$ 关联波廓 $h$ 与速度系数 $\gamma_k$，得到对每个 $k$ 的代数方程（30）。通过函数 $S_\pm(s)$ 的单调性论证，证明至多一个 $k$ 激发，从而聚焦于主模态 $k=1$。这是关键性分析，将无限维问题约化为有限维参数确定。
  4. 显式表达式求导：得到 $\hat{p}(x,y)$ 的闭形式（34）后，通过计算偏导数 $\hat{p}_x$ 和 $\hat{p}_y$，结合矩形区域边界（$x=\pm 1/4$, $y=-1$, $y=0$）及内部临界点条件，定位极值。特别地，$\hat{p}_x=0$ 恒给出 $x=0$ 或 $1/2$，但极值实际在 $x=\pm 1/4$，需验证二阶条件；$\hat{p}_y=0$ 给出关于 $y$ 的超越方程，其解可能为 $y=-1$, $y=0$ 或 $y=y_c$。

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

  1. 常数项 $\beta$ 的确定：由 $h$ 平均为零导出 $\beta=0$，这确保了压力定义与海平面基准一致，避免冗余常数。
  2. 流函数 $\hat{\psi}$ 的构造：通过 $\hat{u}=\hat{\psi}_y$, $\hat{v}=-\hat{\psi}_x$ 定义流函数，将压力表达式中的 $\omega \hat{\psi}$ 项自然引入，并利用 $\psi$ 的显式级数简化计算。
  3. 流反转的判据：通过分析 $\hat{c}-\omega(1+y)$ 在 $y\in[-1,0]$ 的符号变化，结合色散关系（31）中负号的选取，精确给出流反转条件（37）。这使得临界层 $y_c$ 成为分析压力极值的重要候选位置。
  4. 极值位置的分类讨论：基于 $\omega$ 的符号和大小，分无旋、顺流、逆流数种情形。通过研究 $Q(y)$ 的导数符号和边界值，判断其最大值点和最小值点在 $[-1,0]$ 上的位置（端点或内部临界点）。特别地，当 $\hat{c}-\omega(1+y)$ 在 $y_c$ 处为零时，$Q(y)$ 可能具有特殊行为（如导数发散），需单独检查。
  5. 与无旋情况的对比：通过比较 $Q(y)$ 在 $\omega=0$ 时双曲余弦的单调性与 $\omega\ne0$ 时混合双曲函数的非单调性，清晰展示涡度如何扭曲压力分布。技巧上，将 $Q(y)$ 重写为 $-\delta \gamma_1 \left\{[\hat{c}-\omega] \cosh[2\pi\delta(y+1)] - \omega y \cosh(...) + \frac{\omega}{2\pi} \sinh(...)\right\}$，突出线性项 $-\omega y \cosh$ 的非平庸贡献。

与之前工作的比较

a) 突破与改进

本文的核心突破在于 系统性研究常涡度对周期性行波下方压力分布的显著影响，而此前文献主要集中于无旋流情形。具体改进之处：

  • 假设条件的放宽：经典理论（如 [Constantin(2011), ConstantinandStrauss(2010), Oliverasetal.(2023), Umeyama(2018), Zhangetal.(2023)]）均假设背景流无旋（$\Omega=0$），此时流速仅由波势导出，均匀或水平变化。本文引入常涡度 $\Omega \ne 0$，背景流为线性剪切流 $(\Omega(Y+d),0)$，这更符合实际海洋中风生浪与潮流相互作用的场景（如 Ewing(1990) 所支持）。数学上，常涡度保留了 Euler 方程的线性可解性（与均匀涡度类似），但增加了背景流的水平变化，使压力公式中出现 $\omega y$ 项，导致新现象。

  • 结果强度的提升：无旋情况下，已知动压力和总压力的极值（在固定波长下）总出现在波峰和波谷（所有深度），且总压力随深度线性增加。本文证明：

    • 涡度可改变动压力极值的 深度位置：极值可能出现在底床 $y=-1$ 或临界层 $y=y_c$，而非仅表面 $y=0$。
    • 涡度可改变极值对应的 水平位置：当 $Q(y)$ 变号时，最大值可能出现在波谷而非波峰，反之亦然。
    • 涡度不改变总压力随深度增加的趋势（静水项主导），但能显著改变固定深度处极值的 $x$ 位置顺序。 这些结论比无旋情形更精细，揭示了波-流相互作用中涡度作为关键参数的作用。

  • 对物理应用的启示：底部压力传感器常假设静水分布来反演波高，本文表明在强涡度流中，非静水修正可能导致传感器在底床测量到极大值（而非表面波峰对应的极小修正），影响波高估计精度。特别是，当极值出现在临界层时，传感器需避开该深度。

b) 与具体文献的对比

  • 与 [Constantin(2011)] 针对无旋流的研究比较：后者证明动压力在波峰最大、波谷最小，且随深度单调变化。本文显示，涡度引入后，$Q(y)$ 可能非单调，故极值深度可转移至底床或临界层。
  • 与 [Clamond(2013), ClamondandHenry(2020)] 关于非静水修正的工作比较：他们考虑无旋流，重点在波陡效应。本文在线性框架下已出现极值位置的剧变，表明涡度效应可比波陡更显著。
  • 与 [Constantinetal.(2015)] 对潮汐流中波浪的研究比较：后者也考虑常涡度，但关注行波存在性与速度剖面，未深入压力极值。本文补充了压力场的关键分析。
  • 与 [AbrashkinandConstantin(2023)] 对南极环极流涡度的讨论比较：他们提供涡度数量级估计（约 $10^{-4}\,\text{s}^{-1}$），本文则展示即使如此小的涡度，在浅水参数 $\delta$ 适中时，仍可能使 $Q(y)$ 产生可观测的非单调性。

总之，本文通过精细的线性分析，揭示了常涡度对波浪压力分布的深刻影响，放宽了无旋假设，为实际海洋中的压力测量和反演提供了理论基础。

2. 实轴上中间长波方程浅水收敛的奇异性

On the singular nature of shallow-water convergence of the intermediate long wave equation on the real line

作者: Andreia Chapouto, Benjamin Harrop-Griffiths, Guopeng Li, Tadahiro Oh

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20695v1

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

详细分析

中文标题: 实轴上中间长波方程浅水收敛的奇异性

  1. 研究对象和背景
    本文研究对象是中间长波方程（Intermediate Long Wave Equation, ILW）在实轴$ \mathbb{R} $上的解映射正则性性质，特别关注其浅水极限（shallow-water limit, $\delta \to 0$）下的行为。ILW是描述分层流体中界面内波传播的模型，其深度参数$\delta \in (0,\infty)$连续连接Benjamin-Ono方程（BO，深水极限$\delta \to \infty$）和Korteweg-de Vries方程（KdV，浅水极限$\delta = 0$）。原始ILW方程写作

    $$ \partial_t u - G_\delta \partial_x^2 u = \partial_x (u^2), \quad (1.1) $$


    其中非线性算符$ G_\delta $定义为

    $$ G_\delta f(\xi) = -i \left( \coth(\delta\xi) - \delta^{-1}\xi^{-1} \right) \hat{f}(\xi). \quad (1.2) $$


    为在浅水极限下观测有意义的收敛行为，需对振幅进行缩放$ v(t,x) = \delta^{-1} u(\delta^{-3}t, \delta^{-1}x) $，得到缩放ILW方程：

    $$ \partial_t v - \widetilde{G}_\delta \partial_x^2 v = \partial_x (v^2), \quad (1.4) $$


    其中$ \widetilde{G}_\delta $的傅里叶乘子为

    $$ p_\delta(\xi) = \xi^2 \coth(\delta\xi). \quad (2.6) $$


    当$\delta \to 0$时，$ p_\delta(\xi) \to \xi^3 $，故(1.4)收敛到KdV方程：

    $$ \partial_t v + \partial_x^3 v = \partial_x (v^2). \quad (1.5) $$


    此前工作（如[10]）已证明缩放ILW在$L^2(\mathbb{R})$中收敛到KdV，但解映射的正则性呈现奇异现象：已知缩放ILW的解映射在$H^s(\mathbb{R})$中不是$C^2$（[36]），而KdV的解映射在$s \geq -3$时是解析的（[21]）。这种正则性增益（从$C^2$到解析）的机制尚未解释。本文的核心贡献在于通过频率分解揭示这一机制，并阐明浅水极限的奇异性质——即正则性在极限过程中非连续变化，源于不同频率相互作用在$\delta \to 0$时的different singular behaviors。

  2. 主要定理或主要结果
    定理1.1（主要定理）
    (i) 低频率动力学的解析性：设$0 \leq \delta \leq 1$，$s \geq 0$。对任意$\phi \in H^s(\mathbb{R})$，存在$T = T(\|\phi\|_{H^s}) > 0$（独立于$\delta$），使得映射$\phi \mapsto v_{\text{low}}$（低频部分解，见(1.7)）在$H^s$中是解析的。
    (ii) 剩余动力学的$C^2$失效：设$0 < \delta \leq 1$，$s \in \mathbb{R}$。对任意$T > 0$，映射$\phi \mapsto v_{\text{res}}$（剩余部分解，见(1.8)）在零函数处不是$C^2$可微（关于$H^s$正则性）。
    (iii) 低频动力学对KdV的浅水收敛：当$\delta \to 0$时，低频解$v_{\text{low}}$在时间紧开拓扑下的$C(\mathbb{R}; H^s(\mathbb{R}))$中收敛到KdV解$v$。

其中，通过频率投影算子$P_{\delta^{-1}}$（局化于频率$|\xi| \leq \delta^{-1}$），系统分解为：
低频部分：

$$ \begin{cases} \partial_t v_{\text{low}} - \widetilde{G}_\delta \partial_x^2 v_{\text{low}} = P_{\delta^{-1}} \partial_x \left( (P_{\delta^{-1}} v_{\text{low}})^2 \right), \\ v_{\text{low}}|_{t=0} = P_{\delta^{-1}} \phi. \end{cases} \quad (1.7) $$


剩余部分：

$$ \begin{cases} \partial_t v_{\text{res}} - \widetilde{G}_\delta \partial_x^2 v_{\text{res}} = \partial_x (v_{\text{res}}^2) + 2 \partial_x (v_{\text{res}} v_{\text{low}}) + P_{\delta^{-1}}^\perp \partial_x \left( (P_{\delta^{-1}} v_{\text{low}})^2 \right), \\ v_{\text{res}}|_{t=0} = P_{\delta^{-1}}^\perp \phi. \end{cases} \quad (1.8) $$


这里$P_{\delta^{-1}}^\perp = \text{Id} - P_{\delta^{-1}}$。关键性质是$v_{\text{low}}$的频率支集满足$\text{supp} \widehat{v_{\text{low}}}(t,\cdot) \subset \{|\xi| \leq \delta^{-1}\}$（见(1.9)），故其$L^2$范数守恒，从而对任意$s \geq 0$局部适定且可全局延拓（Remark 3.2）。

  1. 研究方法、关键技术和核心工具
    a) 核心分析工具：
    • $X^{s,b}$范数空间：定义
      $$ \|v\|_{X_{\delta}^{s,b}} = \|\langle \xi \rangle^s \langle \tau - p_\delta(\xi) \rangle^b \hat{v}(\tau,\xi)\|_{L^2_{\tau,\xi}}, \quad (3.3) $$

      其中$p_\delta(\xi) = \xi L_\delta(\xi)$为缩放ILW的乘子（见(2.11)），$L_\delta(\xi)$由级数
      $$ L_\delta(\xi) = 6\xi^2 \sum_{k=1}^\infty \frac{1}{k^2\pi^2 + \delta^2\xi^2} \quad (2.12) $$

      给出。$X^{s,b}$空间用于控制线性演化算子的行为，并通过斯特林artz估计获取非线性项的控制。
    • 线性斯特林artz估计：对$0 < \delta \leq 1$，有估计
      $$ \|v\|_{L^4_{t,x}} \lesssim \|v\|_{X_{\delta}^{0,3/8}}, \quad (3.10) $$

      其证明通过周期性版本[30]的直接修改得到，关键在于共振函数的齐次性在低频区$|\xi| \lesssim \delta^{-1}$与KdV一致（Lemma 2.2(i)）。
    • 频率投影技术：通过$P_{\delta^{-1}}$严格分离低频区（$|\xi| \leq \delta^{-1}$）与高频区，低频动力学被简化为与KdV类似的方程（但乘子非均匀收敛），而剩余部分捕获了高频贡献。

b) 关键技术技巧：

  • 共振函数分析：缩放ILW的共振函数
    $$ \Xi_\delta(\xi,\xi_1,\xi_2) = p_\delta(\xi) - p_\delta(\xi_1) - p_\delta(\xi_2) \quad (\xi = \xi_1+\xi_2) $$

    通过展开
    $$ \Xi_\delta(\bar{\xi}) = \Xi_{\text{KdV}}(\bar{\xi}) - \left( \xi^3 h(\delta,\xi) + \xi_1^3 h(\delta,\xi_1) + \xi_2^3 h(\delta,\xi_2) \right), \quad (2.17) $$

    其中$h(\delta,\xi)$由级数(2.14)定义且$h(\delta,\xi) \leq \delta^2 \xi^2$。这显示当$|\xi_{\max}| \lesssim \delta^{-1}$时，$\Xi_\delta \sim \Xi_{\text{KdV}}$（Lemma 2.2(i)），从而低频动力学具有与KdV相同的多线性色散；而当$|\xi_{\max}| \gg \delta^{-1}$时，$|\Xi_\delta| \sim \delta^{-1} \xi_{\min} \xi_{\max}$（Lemma 2.2(ii)），导致色散较弱，这是$C^2$失效的根源。
  • 关键双线性估计（Lemma 3.1）：低频部分的核心估计为
    $$ \left\| P_{\delta^{-1}} \partial_x \left( (P_{\delta^{-1}} v_1)(P_{\delta^{-1}} v_2) \right) \right\|_{X^{s,b'-1}} \lesssim \prod_{j=1}^2 \|P_{\delta^{-1}} v_j\|_{X^{s,b}}, \quad (3.6) $$

    需在$0 \leq \delta \leq 1$上一致成立。证明采用经典[21]的柯西-施瓦茨框架，但需处理$\Xi_\delta$的复杂形式。关键步骤是变量变换
    $$ \mu = \tau - p_\delta(\xi) + \Xi_\delta(\xi,\xi_1,\xi_2) \quad (3.16) $$

    以利用共振函数的非退化性。$d\mu/d\xi_1$的计算（(3.18)）揭示当$|\xi-2\xi_1| \gtrsim |\xi|$时，主导项$3\xi(\xi-2\xi_1)$与KdV情形一致，而剩余项涉及$A_{\delta,k}(\bar{\xi})$（(3.18)末项），需利用$|A_{\delta,k}| \lesssim_k 1$及其对称性进行积分控制。难点在于确保所有估计在$\delta \in [0,1]$上一致，尤其是当$\delta \to 0$时$L_\delta(\xi)$的非一致收敛性（(2.13)参见）。
  • $C^2$失效的证明：适应[36]的方法，考察解映射$\Phi_\delta$在零函数处的二阶Gâteaux导数。取特定方向（如高斯数据）并分析剩余非线性项$2\partial_x(v_{\text{res}} v_{\text{low}})$和$P_{\delta^{-1}}^\perp \partial_x((P_{\delta^{-1}} v_{\text{low}})^2)$。关键观察是“低×高→高”相互作用（即低频与高频碰撞产生高频）由$P_{\delta^{-1}}^\perp$截断，其大小依赖于$\delta^{-1}$，且当$|\xi_{\max}| \gg \delta^{-1}$时共振函数行为如$\delta^{-1} \xi_{\min} \xi_{\max}$，导致二次型无界，从而二阶导数不存在。
  1. 与之前工作的比较
    a) 突破与改进：
    • 此前研究（如[36]）已知缩放ILW解映射整体不是$C^2$，但未解释为何在浅水极限中会获得KdV的正则性（解析性）。本文通过分解机制首次阐明：低频部分（$\text{supp} \hat{v}_{\text{low}} \subset \{|\xi| \leq \delta^{-1}\}$）的解映射在固定$\delta$下已是解析的（Theorem 1.1(i)），且该解析性在$\delta \to 0$下一致保持，从而低频动力学收敛到解析的KdV（Theorem 1.1(iii)）；而剩余部分的$C^2$失效（Theorem 1.1(ii)）源于高频区“低×高→高”的弱色散相互作用，其影响幅度随$\delta \to 0$（高频阈值$\delta^{-1} \to \infty$）而消失。这定量解释了正则性增益：整体解映射的正则性由低频部分主导，而剩余部分的不规则性在极限中被压制。
    • 假设条件的放宽：Theorem 1.1(i)要求数据$\phi \in H^s$（$s \geq 0$），但存在时间$T$仅依赖于$\|\phi\|_{H^s}$而非$\delta$，这比局部适定性理论中$T_\delta \to 0$（当$\delta \to 0$）的通常行为更强。此前类似结果（如[21]对KdV）依赖具体色散关系，本文则将分析扩展到缩放ILW乘子$p_\delta(\xi)$的非一致极限情形，通过Lemma 3.1的δ-一致估计克服了$p_\delta$在无穷频的奇异收敛（$p_\delta(\xi)/\xi^3 \to 0$当$|\xi| \gg \delta^{-1}$）。
    • 现象验证的深化：浅水极限的“奇异性质”不仅体现在解映射正则性的跳跃（从低频部分的解析到整体非$C^2$），还体现在守恒律和测度收敛的塌缩（见Remark 1.5）。本文通过动力系统分解，将现象归结于频率空间的几何效应：低频区色散强（$\Xi_\delta \sim \xi_1\xi_2\xi$），高频区色散弱（$\Xi_\delta \sim \delta^{-1}\xi_{\min}\xi_{\max}$），且两者通过线性投影耦合。
    • 与深水极限的对比：深水极限（$\delta \to \infty$）中，ILW收敛到BO，且收敛是正则性的（解映射连续依赖），因为BO的乘子$|\xi|\xi$与ILW乘子$\xi^2\coth(\delta\xi)$的差为$O(\delta^{-2})$（(2.8)），一致收敛。而浅水极限中，$p_\delta$的收敛仅在固定频率成立，无穷频表现不同，导致非一致估计和奇异行为，这是浅水问题固有难度。

b) 结果强度比较：

  • 正则性结论：KdV解映射在$s \geq -3$解析（[17]），本文Theorem 1.1(i)将解析性扩展到缩放ILW的低频部分，但仅适用于$s \geq 0$（因$X^{s,b}$嵌入需要$s \geq 0$）。对于$s < 0$，ILW全局适定性未知，故受限。
  • 适定性范围：KdV在$H^{-3/4}$以下适定（[23]），而缩放ILW在负正则性上已知 ill-posed（[6] for $s < -1/2$）。本文剩余部分的$C^2$失效（Theorem 1.1(ii)）对**所有$s \in \mathbb{R}$**成立，比[36]的$s \geq 0$更广，因为证明仅依赖线性估计和局部行为。
  • 收敛模式：此前[10]证明缩放ILW整体在$L^2$中收敛到KdV，但未涉及正则性机制。本文Theorem 1.1(iii)表明低频部分本身收敛到KdV，且收敛在$C(\mathbb{R; H^s})$中成立（紧开拓扑），这强于点态收敛。证明使用了半线性论证（无需求完整可积性），而[10]依赖ILW的可积结构（散射理论）。这提供了不同视角：本文聚焦局部动力学，[10]处理全局散射。
  • 周期性情形：Remark 1.3指出，在圆$\mathbb{T}$上，BO解映射是解析的（[34]），故缩放ILW在周期边界下无“低×高→高”不稳定性，正则性增益现象仅为实轴特有。这突出了实轴无穷频区的影响。

综上，本文通过精细频率分解和共振分析，揭示了中间长波方程浅水收敛中正则性增益的微观机理，将整体解映射的奇异极限归因于高频部分的不规则性，而低频部分正则性一致保持并收敛到KdV。这些结果深化了对色散方程极限行为的理解，突显了频率几何与乘子渐近性的关键作用。

3. 关于最小压力梯度原理的数学分析与物理启示

On the Mathematical Analysis and Physical Implications of the Principle of Minimum Pressure Gradient

作者: Haithem Taha

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

发布日期: 2026-02-24

arXiv ID: 2602.20637v1

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

详细分析

中文标题: 关于最小压力梯度原理的数学分析与物理启示

1. 研究对象和背景

本文研究的核心对象是**不可压缩Navier-Stokes方程（INSE）与最小压力梯度原理（PMPG）**之间的数学等价关系。不可压缩Navier-Stokes方程描述粘性不可压缩流体的运动，其标准形式为：

$$ \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} = -\frac{1}{\rho} \nabla p + \nu \nabla^2 \mathbf{u} + \frac{1}{\rho} \mathbf{f}, \quad \nabla \cdot \mathbf{u} = 0 $$

其中 $\mathbf{u}(\mathbf{x},t)$ 表示速度场，$p(\mathbf{x},t)$ 为压力场，$\rho > 0$ 为常数密度，$\nu > 0$ 为运动粘度，$\mathbf{f}(\mathbf{x},t)$ 为外加体积力（如重力）。定义在光滑有界区域 $\Omega \subset \mathbb{R}^3$ 上，并满足无穿透边界条件：

$$ \mathbf{u} \cdot \mathbf{n} = h(\mathbf{x}) \quad \text{on} \quad \partial\Omega $$

这里 $\mathbf{n}$ 为边界外法向，$h(\mathbf{x})$ 为 prescribed 光滑函数，通常对于 impermeable 边界取 $h=0$。不可压缩性条件 $\nabla \cdot \mathbf{u} = 0$ 要求速度场属于散度自由函数空间，结合边界条件，自然引出 Helmholtz-Leray 投影的几何框架：任何平方可积向量场可唯一分解为散度自由分量（满足无穿透条件）与梯度分量的正交和（在 $L^2$ 内积下）。压力梯度 $\nabla p$ 正是该正交补空间（梯度场）的投影，起到约束力作用，强制维持连续性约束。

最小压力梯度原理（PMPG）是 Gauss 最小约束原理在连续介质力学中的推广。Gauss 原理指出，受约束力学系统在每瞬时加速时，其运动方向应使约束力范数最小化。类比地，对于不可压缩流，压力力是维持 $\nabla \cdot \mathbf{u}=0$ 的约束力，因此 PMPG 断言：流场从一瞬态到下一瞬态的演化，在满足不可压缩性和边界条件的所有运动学可接纳加速度中，实际发生的演化使压力梯度的 $L^2$ 范数最小。数学上，定义 PMPG 成本泛函为：

$$ A\left(\frac{\partial \mathbf{u}}{\partial t}\right) = \frac{1}{2} \int_\Omega \rho \left| \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} - \frac{1}{\rho} (\nabla \cdot \boldsymbol{\tau} + \mathbf{f}) \right|^2 d\mathbf{x} $$

其中 $\boldsymbol{\tau}$ 为粘性应力张量（对牛顿流体，$\nabla \cdot \boldsymbol{\tau} = \mu \nabla^2 \mathbf{u}$，$\mu$ 为动力粘度），项 $\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} - \frac{1}{\rho} (\nabla \cdot \boldsymbol{\tau} + \mathbf{f})$ 表示“自由运动”与“实际加速度”的偏差，其平方积分最小化等价于最小化压力梯度范数。优化在** admissible 加速度空间** $\Theta$ 上进行：

$$ \Theta = \left\{ \mathbf{v} \in L^2(\Omega) \ \middle| \ \nabla \cdot \mathbf{v} = 0 \ \text{in} \ \Omega, \ \mathbf{v} \cdot \mathbf{n} = 0 \ \text{on} \ \partial\Omega \right\} $$

该空间确保任何候选加速度都保持不可压缩性和齐次法向边界条件（因边界值固定，变化需满足 $\delta \mathbf{u} \cdot \mathbf{n} = 0$）。

PMPG 的物理启示在于：它提供了 Navier-Stokes 动力学的变分因果解释。例如，流动分离、涡旋生成等复杂现象，可视为系统通过调整流场以最小化所需压力力，任何替代演化若需要更大压力梯度，则被物理过程排除。这区别于传统的平稳作用量原理（如 Hamilton 原理），PMPG 是瞬时、局部的最小化，每时刻独立优化，不涉及时间积分。

本文背景还涉及经典 Galerkin 投影：在有限维子空间（如模态展开）中，NSE 的 Galerkin 近似等价于将残差正交投影到子空间。PMPG 将此思想推广至非线性、非模态表示，适用于模型降阶和复杂基函数选择。此外，论文探讨 PMPG 与升力变分理论（如 Kutta-Joukowski 定理）的联系，并猜想其在稳定性判据和零粘度极限收敛性中的潜在应用。

2. 主要定理或主要结果

本文的核心数学结果是建立了 Navier-Stokes 方程与 PMPG 之间的双向（if and only if）等价性。具体定理陈述如下（基于 Lemma 1 和充分性论证）：

定理（双向等价定理）：设 $\Omega \subset \mathbb{R}^3$ 为有界光滑区域，$\rho > 0$ 常数。假设速度场 $\mathbf{u}(\mathbf{x},t)$ 满足以下正则性条件：对任意 $t \in [0,T]$，有 $\mathbf{u} \cdot \nabla \mathbf{u} \in L^2(\Omega)$，$\nabla \cdot \boldsymbol{\tau} \in L^2(\Omega)$，且外力 $\mathbf{f} \in L^2(\Omega)$。初始条件满足 $\nabla \cdot \mathbf{u}(\mathbf{x},0) = 0$ 于 $\Omega$ 及 $\mathbf{u}(\mathbf{x},0) \cdot \mathbf{n} = h(\mathbf{x})$ 于 $\partial\Omega$，其中 $h$ 为光滑给定函数。则 $\mathbf{u}$ 是Navier-Stokes方程

$$ \rho \left( \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} \right) = -\nabla \lambda + \nabla \cdot \boldsymbol{\tau} + \mathbf{f}, \quad \nabla \cdot \mathbf{u} = 0 $$

在 $\Omega \times [0,T]$ 上的解（$\lambda$ 为压力函数，属于 $H^1(\Omega)$），并保持边界条件 $\mathbf{u} \cdot \mathbf{n} = h$ 于 $\partial\Omega \times [0,T]$，当且仅当在每瞬时 $t \in [0,T]$，$\frac{\partial \mathbf{u}}{\partial t}$ 是 PMPG 成本泛函

$$ A\left(\frac{\partial \mathbf{u}}{\partial t}\right) = \frac{1}{2} \int_\Omega \rho \left| \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} - \frac{1}{\rho} (\nabla \cdot \boldsymbol{\tau} + \mathbf{f}) \right|^2 d\mathbf{x} $$

在 admissible 空间 $\Theta = \{ \mathbf{v} \in L^2(\Omega) \mid \nabla \cdot \mathbf{v} = 0 \text{ in } \Omega, \mathbf{v} \cdot \mathbf{n} = 0 \text{ on } \partial\Omega \}$ 上的唯一最小化点。

该定理的证明分为两部分：

  • 必要性：若 $\frac{\partial \mathbf{u}}{\partial t}$ 最小化 $A$，则通过变分法构造 Lagrangian，取一阶变分为零并利用 Green 公式与边界条件，导出 Euler-Lagrange 方程即 Navier-Stokes 方程。这已在 Lemma 1 中严格证明。
  • 充分性：若 $\mathbf{u}$ 满足 Navier-Stokes 方程，则 $\frac{\partial \mathbf{u}}{\partial t}$ 自动满足 admissible 空间 $\Theta$ 中的约束，且由成本泛函 $A$ 的强凸性（作为 $\frac{\partial \mathbf{u}}{\partial t}$ 的二次型）及欧拉-拉格朗日方程的唯一性，$\frac{\partial \mathbf{u}}{\partial t}$ 必为唯一最小化点。关键依赖：约束集 $\Theta$ 是 $L^2$ 的闭线性子空间，且成本泛函在 $\Theta$ 上严格凸。

此外，论文还表明 PMPG 精确对应于 Leray-Helmholtz 投影的变分形式：Navier-Stokes 方程可视为将“自由运动” $\frac{1}{\rho} \mathbf{f} - \mathbf{u} \cdot \nabla \mathbf{u} + \nu \nabla^2 \mathbf{u}$ 正交投影到散度自由空间，投影的梯度分量即 $-\nabla p$。因此，PMPG 提供了投影机制的最小化解释。

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

a) 核心分析学工具

  1. 变分法与 Lagrange 乘子：PMPG 的优化问题为带线性等式约束的二次最小化。论文采用经典 Lagrange 乘子法处理不可压缩性约束 $\nabla \cdot \frac{\partial \mathbf{u}}{\partial t} = 0$。构造 Lagrangian：
$$ \mathcal{L}\left(\frac{\partial \mathbf{u}}{\partial t}, \lambda\right) = A\left(\frac{\partial \mathbf{u}}{\partial t}\right) - \int_\Omega \lambda(\mathbf{x}) \, \nabla \cdot \frac{\partial \mathbf{u}}{\partial t} \, d\mathbf{x} $$

其中 $\lambda(\mathbf{x})$ 为 Lagrange 乘子（即压力函数）。一阶必要条件要求对 admissible 变分 $\delta \frac{\partial \mathbf{u}}{\partial t} \in \Theta$，有 $\delta \mathcal{L} = 0$。通过积分 by parts 将散度项转化为梯度项，并利用边界条件 $\delta \frac{\partial \mathbf{u}}{\partial t} \cdot \mathbf{n} = 0$ 消去边界贡献，最终得到：

$$ \rho \left( \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} \right) + \nabla \lambda - \nabla \cdot \boldsymbol{\tau} - \mathbf{f} = 0 \quad \text{in} \ \Omega $$

即 Navier-Stokes 方程。此推导依赖于光滑性假设以确保积分运算合法性。

  1. Helmholtz 分解与正交投影理论：对于任意 $\mathbf{v} \in L^2(\Omega)$，Helmholtz 分解给出唯一分解 $\mathbf{v} = \mathbf{w} + \nabla \phi$，其中 $\mathbf{w} \in H(\text{div}^0; \Omega)$（散度自由且满足适当边界条件），$\nabla \phi$ 属于梯度空间，且 $(\mathbf{w}, \nabla \phi)_{L^2} = 0$。在 PMPG 中，令 $\mathbf{v} = \frac{1}{\rho} \mathbf{f} - \mathbf{u} \cdot \nabla \mathbf{u} + \nu \nabla^2 \mathbf{u}$，则 $\frac{\partial \mathbf{u}}{\partial t}$ 正是 $\mathbf{v}$ 的散度自由分量。最小化 $A$ 等价于寻求 $\frac{\partial \mathbf{u}}{\partial t}$ 使余项（梯度分量）的范数最小，即正交投影。这揭示了 PMPG 本质是 Leray 投影的变分表述。

  2. 凸优化理论：成本泛函 $A$ 关于 $\frac{\partial \mathbf{u}}{\partial t}$ 是强凸二次型（积分被积函数为平方范数），优化在闭凸集 $\Theta$ 上进行。因此，问题有唯一解，且一阶必要条件（Euler-Lagrange 方程）既是必要也是充分的。这直接支持双向等价：满足 Euler-Lagrange 方程且满足约束的候选必为唯一最小化点。

  3. 函数空间方法：论文涉及 $L^2$、$H^1$ 等 Sobolev 空间。压力 $\lambda$ 自然属于 $H^1(\Omega)$（由方程 $-\Delta \lambda = \text{div} \cdots$ 诱导），速度场在 $L^2$ 中考虑无散度条件。正则性假设（$\mathbf{u} \cdot \nabla \mathbf{u} \in L^2$ 等）确保所有项平方可积，使成本泛函良好定义。

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

  1. 处理非线性项 $\mathbf{u} \cdot \nabla \mathbf{u}$：在瞬时最小化中，$t$ 固定时 $\mathbf{u}$ 视为已知，因此 $\mathbf{u} \cdot \nabla \mathbf{u}$ 作为给定场参与优化，不引入额外非线性。这避免了传统 NSE 理论中时间连续性的非线性耦合难题，使每步优化独立且线性于 $\frac{\partial \mathbf{u}}{\partial t}$。

  2. 边界条件处理： admissible 空间 $\Theta$ 强制 $\mathbf{v} \cdot \mathbf{n} = 0$ 于 $\partial\Omega$，而原始速度场边界条件 $\mathbf{u} \cdot \mathbf{n} = h$ 通过初始条件和 $\frac{\partial \mathbf{u}}{\partial t} \in \Theta$ 自动保持。变分时，$\delta \frac{\partial \mathbf{u}}{\partial t} \cdot \mathbf{n} = 0$，故在 Green 公式产生的边界项 $\int_{\partial\Omega} \lambda \, \delta \frac{\partial \mathbf{u}}{\partial t} \cdot \mathbf{n} \, dS$ 为零，无需额外边界条件处理。这简化了推导，并确保投影与边界兼容。

  3. 双向等价的完整证明：文献 [6,17] 仅证明必要性（最小化 $\Rightarrow$ NSE）。本文补充充分性（NSE $\Rightarrow$ 最小化），其技巧依赖于成本泛函的严格凸性和约束的线性性。若 $\mathbf{u}$ 满足 NSE，则 $\frac{\partial \mathbf{u}}{\partial t}$ 满足 Euler-Lagrange 方程和 $\Theta$ 约束。由于 $A$ 在 $\Theta$ 上强凸，任何满足约束且使 $\delta A = 0$ 的点必为全局唯一最小点。论文可能显式验证了 $A$ 的 Hessian 在 $\Theta$ 上的正定性（因 $A$ 是 $\frac{\partial \mathbf{u}}{\partial t}$ 的 $L^2$ 范数平方，在 $\Theta$ 上等价于范数）。

  4. 有限维推广与 Galerkin 联系：在有限维子空间 $V_h \subset \Theta$（由基函数张成），PMPG 最小化等价于求解

$$ \min_{\mathbf{v}_h \in V_h} \left\| \frac{\partial \mathbf{u}_h}{\partial t} + \mathbf{u}_h \cdot \nabla \mathbf{u}_h - \frac{1}{\rho} (\nabla \cdot \boldsymbol{\tau}_h + \mathbf{f}) - \mathbf{v}_h \right\|_{L^2} $$

这恰好导出经典 Galerkin 投影：将残差正交于 $V_h$。论文强调此框架不要求基函数线性或模态，允许非线性基（如神经网络）或数据驱动基底，这为模型降阶提供新视角。技巧在于将 Galerkin 条件 $\langle \text{residual}, \phi_i \rangle = 0$ 重新解释为 PMPG 最小化在子空间上的离散形式。

  1. 与变分升力理论的联系：论文提及 PMPG 与升力变分理论（如 Weizsäcker 的“最小阻力”原理）的类比，但未深入。潜在技巧是将定常流问题中的功能最小化（如最小化阻力满足 lift 约束）与瞬时 PMPG 对比，指出前者是后者在 asymptotic 稳态下的特例。

4. 与之前工作的比较

a) 突破或改进

本文的主要突破在于确立了 PMPG 与不可压缩 Navier-Stokes 方程之间的双向严格等价，而此前工作仅知单向蕴含关系。

  • 先前工作：Taha 等人在 [6] 首次提出 PMPG，并证明必要条件：若 $\frac{\partial \mathbf{u}}{\partial t}$ 最小化 $A$，则 $\mathbf{u}$ 满足 Navier-Stokes 方程（即 Lemma 1）。这显示 PMPG 是 NSE 的变分基础，但未证明充分性。换言之，已知最小化演化必为 NSE 解，但 NSE 解是否一定最小化 $A$ 当时未定。
  • 本文改进：通过强调成本泛函的强凸性和约束集的闭凸性，严格证明充分性：若 $\mathbf{u}$ 是 NSE 解，则 $\frac{\partial \mathbf{u}}{\partial t}$ 自动最小化 $A$。因此，等价关系成立。这一改进将 PMPG 从“NSE 的一个性质”提升为 “NSE 的等价变分表述”，意味着：
    • 因果解释：流场行为（如分离、涡脱落）可归因于系统局部选择压力梯度最小的演化路径。任何非最小化演化都不可能对应于真实物理过程。
    • 唯一性保证：对光滑解，最小化点唯一，故 NSE 解在瞬时意义下唯一（符合经典理论）。
    • 数值方法启示：直接求解 PMPG 最小化问题（如通过梯度下降）应产生与 NSE 相同的解，为替代数值格式提供理论基础。

在假设条件方面，本文采用光滑性假设：要求 $\mathbf{u} \cdot \nabla \mathbf{u} \in L^2$，$\nabla \cdot \boldsymbol{\tau} \in L^2$，这比经典 NSE 存在性理论中的弱解假设（如 $\mathbf{u} \in L^2(0,T; H^1)$）更强，但适合证明等价性。对于弱解（如 Leray 解），PMPG 可能需在平均意义或广义函数框架下解释，本文未涉及，这是未来方向。

  • 结果加强：

    1. 变分框架明确化：论文清晰将 PMPG 识别为 Leray-Helmholtz 投影的最小化表述。传统投影方法（如 Chorin 投影法）中，压力 Poisson 方程源于离散不可压缩性，而 PMPG 从连续优化角度统一理解：压力梯度是“自由运动”到散度自由空间的投影残差。
    2. 推广 Galerkin 投影：经典 Galerkin 方法要求基底函数线性且满足散度自由条件。本文显示 PMPG 自然适用于非线性基（如径向基函数、神经网络）和非模态表示（如本征正交分解的非线性扩展），因为最小化不依赖基的线性结构，仅需 admissible 空间 $\Theta$ 的闭合性。这为非线性模型降阶提供理论基础。
    3. 物理洞察深化：通过双向等价，本文论证瞬时 PMPG 最小化是 NSE 的因果机制。例如，流动分离发生是因为在分离点附近，最小压力梯度路径导致流线弯曲；若采用其他路径，需更大压力力，违反物理原理。这超越了传统稳定性分析（如线性模态Growth），提供能量最小化视角。
    4. 猜想提出：基于 PMPG，论文提出两个重要猜想：
      • 稳定性必要条件：若某流场是 NSE 解，则其瞬时演化必最小化 $A$，因此对于扰动，若另一演化需更大 $A$，则不稳定。这可能导致新的能量型稳定性判据，超越线性不稳定性（如 Tollmien-Schlichting 波）。
      • 零粘度极限收敛性：当 $\nu \to 0$，NSE 解是否收敛到 Euler 解？PMPG 视角下，最小化压力梯度在极限下应保持，可能给出收敛的充分条件（如压力梯度序列的有界性）。这连接变分问题与奇异扰动理论。

  • 与之前文献对比：

    • Gauss 原理（1828）应用于粒子力学，Newton 方程是其欧拉-拉格朗日方程。本文将此原理系统地移植到连续介质力学，完成类比如下表：
粒子力学连续介质力学（不可压缩流）
Gauss 最小约束原理PMPG（最小压力梯度原理）
约束力（如张力、法向力）压力梯度（约束力）
自由运动 $\mathbf{F}/m$自由运动 $\frac{1}{\rho}\mathbf{f} - \mathbf{u}\cdot\nabla\mathbf{u} + \nu\nabla^2\mathbf{u}$
加速度 $\mathbf{a}$局部加速度 $\frac{\partial \mathbf{u}}{\partial t}$
配置流形散度自由速度场空间
Newton 方程Navier-Stokes 方程
  • 传统变分原理：Hamilton 原理（平稳作用量）给出 NSE 的弱形式，但涉及时间积分，非瞬时。PMPG 是瞬时变分原理，与 Gauss 原理一致，更适合局部动力学解释。
  • Dirichlet 原理：在投影法中，压力满足 Poisson 方程 $\Delta p = \text{div}(\cdots)$，可视为 Dirichlet 问题。但 PMPG 不同于 Dirichlet 最小化能量；它是关于加速度的二次最小化。

总之，本文通过严格数学证明，将 PMPG 确立为 Navier-Stokes 方程的等价变分形式，在理论深度、物理解释和数值应用上均超越先前工作，为不可压缩流研究提供了新范式。

4. 重力作用下定常轴对压缩性自由面流动中退化点处的奇点

The singularity at degenerate points in steady axisymmetric compressible free surface flows with gravity

作者: Lili Du, Chunlei Yang

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20506v1

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

详细分析

中文标题: 重力作用下定常轴对压缩性自由面流动中退化点处的奇点

一、研究对象和背景

本文研究在重力作用下三维轴对称可压缩流体的定常自由面流动中，自由边界在退化点（degenerate points）处的奇异形状（singular shape）。物理模型描述了一个无旋（irrotational）、无 swirl（swirl-free）的轴对称可压缩流体（例如从喷嘴射出的可压缩射流）在重力场中的运动，其流体区域记为 $D \subset \mathbb{R}^3$，自由界面为 $\partial_a D$。流动由欧拉方程（Euler equations）支配：

$$ \begin{cases} (\mathbf{v}\cdot\nabla)\mathbf{v} = -\dfrac{1}{\rho}\nabla p - g\mathbf{e}_z & \text{in } D, \\ (\mathbf{v}\cdot\nabla)\rho + \rho\,\text{div}\,\mathbf{v} = 0 & \text{in } D, \end{cases} $$

其中 $\mathbf{v}$ 是速度场，$\rho$ 是密度，$p$ 是压力，$g>0$ 是重力加速度，$\mathbf{e}_z = (0,0,1)^\top$。对于可压缩流体，状态方程 $p = p(\rho)$ 满足正则性条件 $p \in C^{1,\alpha}([0,\infty))$ 且

$$ p'(\rho) > 0, \quad 2p'(\rho) + \rho p''(\rho) > 0 \quad \text{for } \rho > 0. $$

例如 $\gamma$-律气体 $p(\rho) = A\rho^\gamma$（$\gamma>1, A>0$）满足该条件。

由于关注轴对称流动，速度场可写为 $\mathbf{v} = v_r(x_1,x_2)\mathbf{e}_r + v_z(x_1,x_2)\mathbf{e}_z$，其中 $(x_1,x_2) \in \mathbb{R}^2$ 是轴对称坐标（$x_1 = \sqrt{X^2+Y^2}$ 为径向，$x_2=Z$ 为轴向），$\mathbf{e}_r, \mathbf{e}_z$ 是基向量。系统在 $(v_r,v_z,\rho)$ 下简化为：

$$ \begin{cases} (v_r\partial_{x_1} + v_z\partial_{x_2})\rho + \rho\left(\dfrac{1}{x_1}\partial_{x_1}(x_1 v_r) + \partial_{x_2}v_z\right) = 0 & \text{in } D, \\ (v_r\partial_{x_1} + v_z\partial_{x_2})v_r - \dfrac{v_\theta^2}{x_1} = -\dfrac{1}{\rho}\partial_{x_1}p & \text{in } D, \\ (v_r\partial_{x_1} + v_z\partial_{x_2})v_z = -\dfrac{1}{\rho}\partial_{x_2}p - g & \text{in } D, \end{cases} $$

其中 $v_\theta \equiv 0$（无 swirl）。引入声速 $c(\rho) = \sqrt{p'(\rho)}$ 和马赫数 $M = \sqrt{v_r^2+v_z^2}/c(\rho)$。流动在点 $x=(x_1,x_2)$ 处称为亚声速（subsonic）若 $M<1$（即 $\rho > \rho_{\text{cr},x_2}$），或跨声速/超声速（sonic-supersonic）若 $M \geqslant 1$（$\rho \leqslant \rho_{\text{cr},x_2}$），其中 $\rho_{\text{cr},x_2}$ 是局部临界密度，由

$$ \frac{p'(\rho_{\text{cr},x_2})}{2} + h(\rho_{\text{cr},x_2}) + g\frac{x_2}{2} = \frac{p'(\bar{\rho}_0)}{2}, $$

定义，这里 $\bar{\rho}_0 > 0$ 是大气密度（常数），$h(\rho) = \int_{\bar{\rho}_0}^\rho \frac{p'(s)}{s}\,ds$ 是焓（enthalpy），且 $h(\bar{\rho}_0)=0$。自由边界处密度为 $\bar{\rho}_0$，由 Bernoulli 定律在自由面上给出：

$$ \frac{|\nabla\psi|^2}{2x_1^2\rho_{\bar{\rho}_0}^2} + h(\bar{\rho}_0) + g\frac{x_2}{2} = \frac{p'(\bar{\rho}_0)}{2} \quad \text{on } \partial_a D, $$

或等价地

$$ \frac{|\nabla\psi|^2}{x_1^2} = p'(\bar{\rho}_0) - 2gx_2 \quad \text{on } \partial_a D, $$

其中 $\psi$ 是Stokes 流函数，满足 $\nabla \psi = (x_1\rho v_z, -x_1\rho v_r)$，且 $\psi>0$ 在 $D$ 内，$\psi=0$ 在 $\partial_a D$ 上。流体域可识别为 $\{\psi>0\}$。

将 Bernoulli 定律写成

$$ F(\rho;x_2) := \frac{|\nabla\psi|^2}{2x_1^2} - h(\rho) - g\frac{x_2}{2} = 0 \quad \text{在 } \{\psi>0\}, $$

利用 $F_\rho < 0$（由状态方程条件保证）可从 $F(\rho;x_2)=0$ 隐式解出 $\rho = \rho(|\nabla\psi|^2/x_1^2, x_2)$。代入无旋条件 $\text{div}(\mathbf{v})=0$（即 vorticity $\omega=0$）后，问题转化为关于 $\psi$ 的拟线性 Bernoulli 型自由边界问题：

$$ \begin{cases} \text{div}\left( \dfrac{\nabla\psi}{x_1} \left( \partial_1 \rho\left( \dfrac{|\nabla\psi|^2}{x_1^2}; x_2 \right) , \partial_2 \rho\left( \dfrac{|\nabla\psi|^2}{x_1^2}; x_2 \right) \right) \right) = 0 & \text{在 } \{\psi>0\}, \\ x_1^2 |\nabla\psi|^2 = p'(\bar{\rho}_0) - 2gx_2 & \text{在 } \partial\{\psi>0\}. \end{cases} $$

该问题在自由边界上退化，因为当 $x_2 \to x_2^{\text{st}} := p'(\bar{\rho}_0)/(2g)$ 时，右边趋于零，导致 $|\nabla\psi| \to 0$。退化点出现在三类位置：停滞点（stagnation points）$\{x_2 = x_2^{\text{st}}\}$ 但 $x_1 \neq 0$，非滞轴点（non-stagnation axis points）$\{x_1=0\}$ 但 $x_2 < x_2^{\text{st}}$，以及原点 $(0,0)$（对应物理坐标 $(0,x_2^{\text{st}})$）。在这些点，$\nabla\psi=0$，自由边界条件退化，且方程系数可能奇异（因 $\rho$ 依赖于 $|\nabla\psi|^2$）。

为分析方便，做坐标变换 $y_1 = x_1$, $y_2 = x_2^{\text{st}} - x_2$，并缩放：

$$ u(y_1,y_2) = \frac{\psi(y_1, x_2^{\text{st}}-y_2)}{\sqrt{2\rho_{\bar{\rho}_0}^2 g}}, \quad H(t;y_2) = \rho\left(2\rho_{\bar{\rho}_0}^2 g\, t, x_2^{\text{st}}-y_2\right). $$

则问题化为（在区域 $\Omega \subset \{x_1 \geqslant 0\}$ 相对开集）：

$$ \begin{cases} \text{div}\left( \dfrac{\nabla u}{x_1} \left( \partial_1 H\left( \dfrac{|\nabla u|^2}{x_1^2}; x_2 \right) , \partial_2 H\left( \dfrac{|\nabla u|^2}{x_1^2}; x_2 \right) \right) \right) = 0 & \text{在 } \{u>0\}, \\ x_1^2 |\nabla u|^2 = x_2 & \text{在 } \partial\{u>0\}. \end{cases} $$

解的定义采用亚声速变分解（subsonic variational solution，Definition 1.3）：$u \in W^{1,2}_{\text{w,loc}}(\Omega)$（加权 Sobolev 空间，权重 $x_1$），非负，$u=0$ 在 $x_1=0$ 上，$C^0$ 连续，$C^2$ 在正区域，满足兼容性（轴上导数为零），亚声速条件：存在 $\varepsilon_0>0$ 使得

$$ H\left( \frac{|\nabla u|^2}{x_1^2}; x_2^{\text{st}}-x_2 \right) - \rho_{\text{cr},x_2} \geqslant \varepsilon_0 \quad \text{局部在 } \Omega, $$

且能量泛函（由问题导出）的第一变分公式（1.24）为零。亚声速弱解（Definition 1.9）进一步要求自由边界在非退化点（$x_1>0, x_2>0, \nabla u \neq 0$）附近是 $C^{2,\alpha}$ 曲线。

函数空间 $W^{1,p}_{\text{w,loc}}$ 由加权 $L^p$ 空间 $L^p_{\text{w}}(U) = \{ v : \int_U x_1^{-1}|v|^p\,dx < \infty \}$ 构成，这源于能量泛函中 $x_1^{-1}|\nabla u|^2$ 的积分有界性需求。退化点处 $\nabla u \to 0$，导致自由边界条件 $x_1^2|\nabla u|^2 = x_2$ 的右端也趋于零，故缩放分析需匹配这一阶。

二、主要定理

本文核心成果是完整刻画了上述问题自由边界在所有可能退化点处的渐近形状。主要定理包括：

定理 2.1（停滞点）。设 $u$ 是问题 (1.20) 的亚声速弱解，$\Omega \subset \{x_1 \geqslant 0\}$ 连通相对开，且存在常数 $C$ 使得

$$ \frac{|\nabla u|^2}{x_1^2} \leqslant C x_2 \quad \text{局部在 } \Omega \cap \{u>0\}. $$

对任意停滞点 $x^\circ = (x_1^\circ, 0)$ 满足 $x_1^\circ \neq 0$，考虑爆破序列

$$ u_r(x) := \frac{u(x_1^\circ + r x_1, r x_2)}{r^{3/2}}. $$

则以下之一成立：

  1. 爆破极限 $u_0 := \lim_{r\to 0^+} u_r \equiv 0$，此时加权密度

    $$ D_{x_2}(0+) := \lim_{r\to 0^+} r^{-3} \int_{B_r(x^\circ)} x_2 \chi_{\{u>0\}} \,dx $$

    满足 $D_{x_2}(0+) = 0$ 或 $D_{x_2}(0+) = \dfrac{2}{3}$。

  2. $u_0 \not\equiv 0$，此时 $u_0$ 唯一地由下式给出（Stokes 角流）：

    $$ u_0(x_1,x_2) = \dfrac{3}{4} (x_1^2+x_2^2)^{3/2} \cos\left( \dfrac{3}{2} \arctan \dfrac{x_1}{x_2} \right) \chi_{\{ -\pi < \arctan(x_1/x_2) < \pi \}}, $$

    且加权密度 $D_{x_2}(0+) = \dfrac{\sqrt{3}}{3}$。

收敛均为 $W^{1,2}_{\text{loc}}(\mathbb{R}^2)$ 强收敛和局部一致收敛。

定理（非滞轴点，摘要及第 6 节）。设 $x^\circ = (0, x_2^\circ)$ 是自由边界点，$x_2^\circ > 0$ 且 $x_2^\circ \neq x_2^{\text{st}}$（即非停滞轴点）。则自由边界在 $x^\circ$ 附近的形状只能是向下尖点（downward pointing cusp）或向上尖点（upward pointing cusp），即存在 $\alpha \in (0,1)$ 使得在局部坐标下 $\partial\{u>0\}$ 由 $x_1 = c |x_2|^{\alpha}$ 描述（$c>0$ 常数，向下或向上取决于符号）。

定理（原点，摘要及第 7 节）。在原点 $(0,0)$ 处，自由边界的渐近形状只有两种非平凡可能：

  • Garabedian 尖泡（Garabedian’s pointed bubble）：一种特定的尖点形状。
  • 水平平坦面（horizontal flat surface）：自由边界在原点附近是水平的（即 $x_2 \equiv 0$ 附近）。

注：定理 2.1 中停滞点的爆破缩放选择 $r^{-3/2}$ 源于自由边界条件 $x_1^2|\nabla u|^2 = x_2$ 在 $x_2=0$ 处的退化：若 $u \sim r^\beta$，则 $|\nabla u| \sim r^{\beta-1}$，代入得 $r^{2(\beta-1)} \sim r$，故 $\beta = 3/2$。非滞轴点和原点的缩放可能不同（例如原点处若为尖点，可能 $u \sim r^\alpha$ 且 $\alpha<3/2$）。

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

a) 核心分析学工具

  1. 变分方法与第一变分公式：问题源于变分结构，能量泛函 (1.23) 定义为

    $$ E_F(v;\Omega) = \int_\Omega \left[ \frac{|\nabla v|^2}{x_1} F\left( \frac{|\nabla v|^2}{x_1^2}; x_2 \right) + \lambda(x_2) \chi_{\{v>0\}} \right] dx, $$

    其中 $F$ 和 $\lambda$ 由流体参数确定。亚声速变解要求该泛函在 $u$ 处对任意光滑耗散场 $\phi$ 的第一变分为零，推导出公式 (1.24)。该公式是研究解正则性和奇点的基石，通过选择特定 $\phi$ 可得到局部几何信息。

  2. 单调性公式（Monotonicity formulas）：这是本文首要创新。对于经典的不可压缩 Bernoulli 自由边界问题（如 Alt-Caffarelli-Friedman），Weiss 单调性公式在非退化点处成立，但本文问题因拟线性性（$\rho$ 依赖于梯度）和自由边界退化， Weiss 公式失效。作者在 Section 4 首次建立适用于此类拟线性 Bernoulli 型问题的单调性公式，形式类似：

    $$ \frac{d}{dr} \left( r^{-2} \int_{B_r^+(x^\circ) \cap \{u>0\}} \Phi\left( \frac{|\nabla u|^2}{x_1^2}; x_2 \right) x_1^{-1} dx \right) \geqslant \text{（几何项）}, $$

    其中 $\Phi$ 是与 $F$ 相关的 Enrichment 函数（由 $\Phi = 2t \partial_t F - F$ 定义），$t = |\nabla u|^2/x_1^2$。该公式在方程奇异（马赫数趋近 1，即 $\rho \to \rho_{\text{cr}}$）和自由边界退化（$x_2 \to 0$ 或 $x_1 \to 0$）时均成立，是分析奇点的核心工具，用于控制爆破序列的极限行为。

  3. 频率函数与频率公式：对于原点处的水平平坦点（$u$ 在原点附近恒为零），标准单调性公式失效。作者在 Section 8 建立非线性频率公式，类似于 Almgren 频率函数，但针对当前拟线性结构。定义“频率” $N(r)$ 为某种增长比率的极限，证明其收敛到某整数或半整数，从而区分 Garabedian 尖泡（对应特定频率）与水平平坦（频率可能不同）。这结合补偿紧性理论，确保变分解序列的强收敛。

  4. 补偿紧性理论（Compensated compactness）：在 Section 10，为证明变分解序列的强收敛（排除模糊收敛带来的额外奇点），将问题与可压缩 Euler 方程联系起来。利用 divergence-free 结构（质量守恒）和熵条件，应用 compensated compactness 技术（如 Tartar 定理）证明序列 $\rho^n = H(|\nabla u_n|^2/x_1^2; x_2)$ 在 $L^1$ 中强收敛，从而保证极限解满足原方程。这对处理可压缩流体的密度振荡至关重要。

  5. 爆破分析（Blow-up analysis）：对每类退化点，构造适当的缩放序列（如停滞点用 $r^{-3/2}$，轴点可能不同），利用单调性公式证明极限函数 $u_0$ 存在且满足某个简化方程（常微分方程或已知驻波）。通过分析极限方程的解空间，得出 $u_0$ 只能是特定形式（Stokes 角流、尖点、平坦等）。

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

  1. 处理拟线性与退化边界：问题中 $\rho$ 通过 $|\nabla u|^2/x_1^2$ 隐式依赖，导致方程系数在 $|\nabla u| \to 0$ 时奇异。作者通过引入辅助函数 $F(t;s)$ 和 $\Phi(t;s)$，将变分公式转化为可控制的单调性量。关键步骤是计算第一变分 (1.24) 并识别出 $\Phi$ 的单调性贡献，即使 $t \to 0$ 时 $F$ 和 $\partial_t F$ 可能无界，但组合 $\Phi$ 保持良好性质。

  2. 退化点分类与缩放选择：三类退化点（停滞、非滞轴、原点）的物理和数学本质不同：

    • 停滞点 ($x_2=0, x_1 \neq 0$)：自由边界条件要求 $|\nabla u|^2 \sim x_2$，故缩放 $u \sim r^{-3/2}$ 使极限方程非平凡。证明中需排除平凡极限 $u_0 \equiv 0$ 时加权密度以外的可能性，并证明非平凡极限必为 Stokes 角流（一个已知的 Euler 角流解）。
    • 非滞轴点 ($x_1=0, x_2>0$)：轴上条件 $u=0$ 和 $x_1\partial_{x_1}u \to 0$ 暗示解在轴上平缓。通过不同缩放（如 $u \sim r^\alpha$），从边界条件 $x_1^2|\nabla u|^2 = x_2$ 和方程推导出 $\alpha \in (0,1)$，且 $u_0$ 必须满足 $u_0(x_1,x_2) = c x_1^{\beta} \phi(x_2/x_1^{\alpha})$，最终化为 ODE 解为尖点。
    • 原点：最困难，因两种退化同时发生。若原点处 $u_0 \not\equiv 0$ 但 $\nabla u_0(0)=0$，则可能为尖点或平坦。通过频率公式，证明若频率为特定值则对应 Garabedian 尖泡，否则必须水平平坦。水平平坦点的证明需结合补偿紧性，密度序列强收敛到常数 $\bar{\rho}_0$，无振荡。

  3. 加权密度的几何意义：加权密度 $D_{x_2}(0+)$ 本质上是自由边界“体积”在 $x_2$ 权重下的渐近比例。对于 Stokes 角流，计算得 $D_{x_2}(0+)=\sqrt{3}/3$；平凡极限时 $0$ 或 $2/3$ 对应不同几何（如 $0$ 可能表示自由边界在球体内测度趋于零，$2/3$ 可能对应某种退化尖点）。这量化了奇点的“强度”。

  4. 排除其他可能性：通过单调性公式控制爆破极限，并结合极限方程的解的分类（如停滞点极限满足 $u_0$ 为 Stokes 角流或零），严格排除了其他渐近行为（如其他角流或更奇异的形状）。

四、与之前工作的比较

本文工作标志着可压缩流体自由边界奇点理论的重要进展，主要体现在以下方面：

  1. 问题框架的推广：

    • 从不可压缩到可压缩：经典结果（如 Vǎrvǎrucǎ-Weiss [Comm. Pure Appl. Math., (67), 2014]）针对二维不可压缩 Euler 方程（密度常数），处理非退化自由边界点的正则性。本文研究三维轴对称可压缩情形，密度变化，状态方程满足 (1.2)，引入了局部临界密度 $\rho_{\text{cr},x_2}$ 和亚声速条件。这是对可压缩流体自由边界问题的首次系统奇点分析。
    • 从非退化到退化点：Alt-Caffarelli-Friedman 系列工作 [ACF84, ACF85] 奠定了 Bernoulli 型自由边界问题的正则性理论，但假设在自由边界上 $|\nabla u|>0$（非退化），即边界 $C^{1,\alpha}$ 且无角点或尖点。本文专门针对退化点（$|\nabla u|=0$）研究，这些点物理上对应驻点或轴上点，是流动中自然出现的奇点（如 jets 的端点、水滴形状等）。此前对这类点的分析仅在特定模型（如二维水波）中有零星结果，本文给出完整分类。

  2. 假设条件的放宽：

    • 允许方程奇异：当马赫数 $M \to 1$（即 $\rho \to \rho_{\text{cr}}$）时，拟线性方程系数退化（因 $\partial_t H \to \infty$）。本文方法不要求 $M$ 远离 1，亚声速条件仅要求 $\rho \geqslant \rho_{\text{cr}}+\varepsilon$，允许接近临界密度。
    • 允许自由边界退化：在停滞点 ($x_2=0$) 和轴点 ($x_1=0$)，自由边界条件 $x_1^2|\nabla u|^2 = x_2$ 本身退化（右端为零）。此前工作通常避开这些点或假设边界光滑。本文的单调性公式和爆破分析专门处理这种退化。
    • 更一般的解概念：采用亚声速弱解（Definition 1.9），弱于经典解甚至粘性解，允许自由边界在非退化点外只有局部 $C^{2,\alpha}$，且无全局正则性假设。这覆盖了可能的奇点。

  3. 结果的强度与完整性：

    • 精确渐近形状：不仅证明奇点存在，而且给出每种退化点处自由边界的显式渐近形式：停滞点处 Stokes 角流（解决了 Stokes 在可压缩框架下的猜想），非滞轴处尖点（向上或向下），原点处 Garabedian 尖泡或水平平坦。这比仅知“存在角点”强得多。
    • 解决 Stokes 猜想的广义版本：Stokes 在 1880 年猜测二维不可压缩无旋重力波中，角点处的自由边界形状必为 Stokes 角流。本文在三维轴对称可压缩框架下证明类似结论（停滞点），且适用于可能跨声速的流动，是猜想的重大推广。
    • 加权密度量化：引入加权密度 $D_{x_2}(0+)$ 等量，提供奇点强度的数值不变量（如 Stokes 角流对应 $\sqrt{3}/3$），可进一步用于分类或稳定性分析。

  4. 方法论的创新：

    • 首次建立拟线性 Bernoulli 型问题的单调性公式：这是核心技术突破，克服了 Weiss 公式失效的障碍。公式构造依赖于函数 $F$ 的凸性（由状态方程 (1.2) 保证）和坐标变换后的结构，适用于广泛的一类拟线性椭圆自由边界问题。
    • 频率公式的新发展：在原点水平平坦点，建立了非线性频率公式，结合补偿紧性，这是处理此类奇点的新工具，可能推广到其他退化椭圆问题。
    • 系统化奇点分析框架：论文提供了一套从单调性公式 → 爆破分析 → 频率分析 → 补偿紧性的完整流程，为研究其他奇异自由边界问题（如一般拟线性方程、其他退化类型）提供了模板。

  5. 与 Vǎrvǎrucǎ-Weiss 工作的关系：Vǎrvǎrucǎ-Weiss 研究二维不可压缩自由表面流动，证明非退化点正则，并在某些退化点发现角点。本文可视为其可压缩对应：将不可压缩（密度 $\equiv1$，$F \equiv t/2$）推广到可压缩（$F$ 复杂），且从二维到三维轴对称（引入 $x_1$ 权重）。更重要的是，本文处理了 Vǎrvǎrucǎ-Weiss 未触及的停滞点和原点的完整分类，并证明了 Stokes 角流在可压缩中的唯一性。

总之，本文在模型、方法和结果上均超越先前工作：从不可压缩/非退化到可压缩/退化，从局部正则性到全局奇点分类，并发展出全新的分析工具，为可压缩流体自由边界理论树立了新标杆。

5. 单涡环的全局动力学

Global dynamics of a single vortex ring

作者: Dengjun Guo, In-Jee Jeong, Lifeng Zhao

类别: math.AP

发布日期: 2026-02-23

arXiv ID: 2602.20131v1

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

详细分析

中文标题: 单涡环的全局动力学

论文研究三维不可压缩欧拉方程下涡环的全局时间动力学。系统描述在轴假设无旋（axisymmetric flow without swirl）条件下，速度场和涡度场在柱坐标下具有特定形式：速度 $U = u_r(r,z,t)e_r + u_z(r,z,t)e_z$，涡度 $\Omega = w(r,z,t)e_\theta$。对应的涡度方程（1.3）为：

$$ \partial_t w + \nabla \cdot (u w) = 0, $$

其中 $\nabla \cdot (u w) = \partial_r(u_r w) + \partial_z(u_z w)$，定义在半平面 $H = \{ x = (r,z) : r > 0, z \in \mathbb{R} \}$ 上。速度由Biot-Savart律给出（公式(2.1)）：

$$ u(x,t) = \int_H K(x,x') w(x',t) \, dx', $$

核函数 $K = (K_r, K_z)$ 的具体形式涉及椭圆积分。研究关注涡环沿对称轴（$z$ 方向）的传播行为，其宏观不变量包括：涡量总量 $M_0(w) = \int_H w \, dx$，二阶矩 $M_2(w) = \int_H r^2 w \, dx$，以及具有对数奇性的能量 $E(w)$（公式(1.4)）。初始数据被限制为“薄涡环”构型，即满足 $w_0(x) \lesssim \epsilon^{-2}/r$，且 $M_0, M_2, E$ 与半径为 $r_0$、环量为 $\mu$ 的理想薄环（其涡度集中在 $r \approx r_0$ 附近）的宏观量相差 $o(1)$ 或 $O(|\log \epsilon|^{-1})$。核心现象是经典的 Kelvin-H Hicks 公式 预测的传播速度：

$$ V_z = \frac{\mu}{4\pi r_0} (|\log \epsilon| + O(1)). $$

涡丝猜想 预言：此类集中涡结构将保持集中，并以上述速度沿 $z$ 轴运动。验证此猜想需解决两大技术障碍：一是克服 Biot-Savart 核在 $r=0$ 附近的奇异性导致的轴向追踪误差累积；二是证明即使存在潜在 丝化现象（filamentation），涡环芯仍能保持全局时间的主导传播。

主要定理

论文证明了以下核心结果：

定理 1.1：在假设 (i)-(v) 下（见下文），解 $w_\epsilon(\cdot, t)$ 对所有 $t \in \mathbb{R}$ 保持集中在某点 $x^*(t) = (r^*(t), z^*(t)) \in H$ 附近，满足：

  • (a) $|r^*(t) - r_0| \lesssim |\log \epsilon|^{-1}$，且 $\int_{|x-x^*| \geq \epsilon^{1/2}} (1+r^2) w \, dx \lesssim |\log \epsilon|^{-1}$。
  • (b) 若进一步假设 $\int_H \langle z \rangle r^2 w_0 \, dx < +\infty$，则
    $$ \left| z^*(t) - \frac{\mu}{4\pi r_0} |\log \epsilon| \, t \right| \lesssim 1 + A_0 + |t|, $$
    其中 $A_0 = \int_H \langle z \rangle r^2 w_0 \, dx$。这给出了 Kelvin-Hicks 速度的全局时间验证。

定理 1.3 （普遍丝化）：若初始涡度在径向方向有显著分布（存在子集 $U \subset \text{supp}(w_0)$ 满足 (1.5) 的几何条件），则其像在速度场下的轴向直径满足：

$$ \text{diam}_z(X(U,t)) \gtrsim \frac{\mu}{8\pi r_0} |\log \epsilon| \, t - C, $$

即 线性时间增长。

推论 1.5 & 定理 1.6：丝化导致涡度二阶径向导数范数 $\|\partial_{rr} w_\epsilon\|_{L^\infty}$ 呈 $\sim \sqrt{t}$ 增长，进而证明 薄涡环构型在 $W^{2,\infty}$ 范数下动力学不稳定。

假设条件 (i)-(v) 的核心是宏观不变量约束：

  • (i) $0 \le w_0(x) \lesssim \epsilon^{-2}/r$ (有界集中性)；
  • (ii) $|M_0(w_0) - \mu| \lesssim |\log \epsilon|^{-1}$；
  • (iii) $|M_2(w_0) - r_0^2 \mu| \lesssim |\log \epsilon|^{-1}$；
  • (iv) $\left| E(w_0) - \frac{r_0 \mu^2}{2\pi} |\log \epsilon| \right| \lesssim 1$；
  • (v) $\int_H \langle z \rangle r^2 w_0 \, dx < +\infty$ (关于 $z$ 的弱衰减控制)。

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

证明框架分为两大支柱：全局集中性与轴向追踪 及 普遍丝化机制。

1. 全局集中性与轴向追踪的克服

A. 加权动量障碍 $A(t)$ 的引入：传统方法追踪涡度质心 $\bar{z}(t) = \frac{\int z w \, dx}{\int w \, dx}$，其在时间演化中因 Biot-Savart 核在 $r=0$ 附近的奇异性，需进行径向截断（如取 $r \ge c$），这会破坏对称性并产生形如 $\int |z| w \, dx$ 的误差项。由于方程不提供轴向二阶矩 $M_z = \int z^2 w \, dx$ 的守恒，这些项对长时不可控。关键创新 是定义加权量：

$$ A(t) := \int_H r^2 \langle z - V_\epsilon t \rangle w(x,t) \, dx, \quad \langle z \rangle := \sqrt{1+z^2}, $$

其中 $V_\epsilon = \frac{\mu}{4\pi r_0} |\log \epsilon|$ 是期望速度。$r^2$ 的权重至关重要：

  • 消除径向奇点：$r^2$ 精确抵消了 Biot-Savart 核在 $r \to 0$ 时的 $1/r$ 奇异性，使积分无需径向截断即可良好定义。
  • 实现误差抵消：计算 $A'(t)$ 时，主部与能量 $E(w(t))$ 精确成比例（见引理 2.8）。更精妙的是，$r^2$ 权重触发了一种 反对称结构，使得所有含 $\int |z| w \, dx$ 的误差项在 $\nabla \cdot (u w) = 0$ 的框架下 精确抵消。若换为其他权重（如 $r$），将遗留不可控项。

由此得到 关键估计（引理 2.8）：存在 $\epsilon_0$ 和常数 $C$ 使得对所有 $t \ge 0$ 和 $0 < \epsilon \le \epsilon_0$ 有 $|A'(t)| \le C$。这建立了 尺度分离：涡环以 $O(|\log \epsilon|)$ 速度平移，但其内部变形或“泄漏”在移动参考系中仅为 $O(1)$。

B. 全局径向与轴向的加权浓度：利用能量 $E$ 和守恒量 $M_0, M_2$ 的耦合关系，通过精细的加权估计（引理 2.10, 2.11），证明：

$$ \int_{|x-x^*| \ge R\epsilon} (1+r^2) w(x,t) \, dx \lesssim \frac{1}{\log R}, $$

其中 $x^*(t)$ 是某种加权质心。这强于 Jerrard 的原始估计（ proposition 2.5 仅控制 $\int_{|x-x^*| \ge \epsilon^{1/2}} w \, dx$）。该估计蕴含：

  • 径向严格约束：涡度几乎全集中在 $r \approx r_0$ 的 $O(|\log \epsilon|^{-1})$ 邻域内；
  • 轴向约束增强：结合 $A(t)$ 的有界性，可推导出 $z^*(t)$ 的线性追踪（定理 1.1(b)）。

证明中使用了经典的 浓度-紧性论证：若大半径处有显著质量，则能量会高于集中构型，而守恒的能量约束了这种扩散。

2. 普遍丝化机制：双速度框架与几何阈值

A. 物理直觉与几何阈值：涡环芯（$d \sim \epsilon$）的平移速度 $V_\epsilon \sim |\log \epsilon|$。距芯距离 $d$ 处的局部诱导速度 $\sim d^{-1}$。若某部分涡度初始位于 $d \gg |\log \epsilon|^{-1}$，则其局部诱导速度远小于芯速度。在芯的移动参考系中，这些外围部分将 滞后，导致涡丝被拉伸。

B. 双速度定义与丝化证明：定义沿流形 $X(\cdot,t)$ 的pull-back速度场 $v(x,t)$，其满足 $v(X(x,t),t) = u(X(x,t),t) - V_\epsilon e_z$，表示在平移芯的参考系中的速度。利用测度保持性 $\det(DX) = 1$ 和能量估计，证明（定理 1.8）：

  • 存在初始时刻的“慢部分” $w_{\text{d},0,\epsilon}$（其支持在 $r \ge c_5$ 且满足 (1.6) 的质量与尺度条件），其在流形下的像满足轴向直径线性增长。
  • 条件 (1.6) 本质上保证初始支持在径向方向延伸至 $O(|\log \epsilon|^{-1})$ 以上，从而其粒子速度 $\sim d^{-1} \ll |\log \epsilon|$。

C. 从几何拉伸到范数增长：支撑直径的线性拉伸 $\Delta z \sim t$，由不可压缩性（体积守恒）导致径向压缩 $\Delta r \sim t^{-1/2}$，从而径向梯度增大。通过构造光滑闭曲线 $\Gamma$ 包围高涡度区域（推论 1.5 的假设），将拉伸转化为 $\|\partial_{rr} w\|_{L^\infty}$ 的 $\sqrt{t}$ 增长。

D. 不稳定性结论（定理 1.6）：通过选择适当的扰动项 $w_{\star,0,p}$ 进行叠加，使扰动后初始数据满足丝化条件。由于扰动可任意小（在 $L^p$ 或 $W^{k,p}$ 意义下），这证明了薄涡环构型是 动力不稳定的。

与之前工作的比较

论文声称的突破主要体现在 全局时间验证涡丝猜想 和 发现普遍丝化机制 两方面，显著放宽并推广了前人的结果。

1. 全局时间验证涡丝猜想：从短时到全局

  • 前人工作局限：
    • BCM00 (Benedetto–Caglioti–Marchioro)：证明了对于广义集中初始数据，涡度始终集中在某个点附近（全局 $L^1$ 集中），但无法追踪该点 $z^*(t)$ 的精确轨迹。其分析隐含的时间尺度限制为 $t \lesssim |\log \epsilon|^{-1}$（见论文中公式 (2.5) 的讨论），因为其依赖的质心追踪在长时会累积误差。
    • Jer20 (Jerrard)：改进了径向浓度估计，但仍未解决轴向长时追踪问题。其方法同样受困于 $r=0$ 奇点导致的误差。
    • 有限时间行波研究：如 [Nor73, FB74, Ben76] 仅构造了特解（行波），不涉及通用初始数据。
  • 本文改进：
    • 假设条件大幅放宽：仅需 宏观不变量（$M_0, M_2, E$）与理想薄环一致，无需 初始数据接近任何特定行波解，也无需 预设涡度剖面（如 Fraenkel 的剖面假设 (1.7)）或强的 $z$ 方向衰减假设（仅需 $\int \langle z \rangle r^2 w_0 \, dx < \infty$）。
    • 时间尺度突破：首次将验证扩展到 所有 $t \in \mathbb{R}$，无需 $|\log \epsilon|^{-1}$ 限制。核心工具是 $A(t)$ 及其导数估计，从根本上解决了轴向误差累积问题。
    • 精度：给出传播速度 $z^*(t)$ 的主部精确符合 Kelvin-Hicks 公式，误差界为 $O(1 + A_0 + |t|)$，虽非最优（因 $O(1)$ 项依赖于剖面），但在仅知宏观量的框架下是Sharp的（Remark 1.2）。

2. 普遍丝化：从对称性假设到通用构型

  • 前人工作局限：
    • 二维类比：丝化/不稳定性的严格证明通常依赖 附加对称性。例如：
      • 在无界半平面，[ISG99, Ift07] 需 $z$ 方向的奇对称性（$w(r,z) = -w(r,-z)$）证明线性丝化。
      • 在周期域或有界域，[Den09, KvS14, Zla15] 等常对大涡旋的 奇/偶扰动 进行分析。
      • 对于稳定行波（如 Lamb 偶极子），[CJ22a, JYZ25] 需族中特定不稳定扰动。
    • 三维情况：[CJ23] 证明 Hill 涡旋（另一类行波）附近的 $W^{2,\infty}$ 不稳定；[Elg21, EGM21, CJ25, KJ22, CMZZ25] 等研究 反平行涡环对（本质是 $z$ 方向奇对称），其不稳定性源于环间相互作用。
    • 共同点：所有严格丝化结果要么针对 特定行波，要么假设 全局 $z$ 方向奇对称性（强制正负涡度配对），这些条件排除了单孤立涡环的普遍行为。
  • 本文改进：
    • 通用性：定理 1.3 和 1.8 仅要求初始数据满足 宏观不变量 和 一个几何条件：即存在相当质量 ($\gtrsim \mu |\log \epsilon|^{-2}$) 位于 $r \ge c_5$（与 $|\log \epsilon|^{-1}$ 同阶）。这适用于 任意 满足宏观约束的单涡环，不依赖 任何对称性或接近特定行波。
    • 机制普适：证明揭示了丝化是 速度失配（芯高速平移 vs 外围弱局部诱导）的必然几何后果，是一种 普适的非线性不稳定性，而非依赖精细谱分析或特定扰动模式。
    • 强不稳定结论：直接证明 $\|\partial_{rr} w\|_{L^\infty} \sim \sqrt{t}$ 增长，且通过任意小扰动即可触发（定理 1.6），表明几乎任何薄环构型在 $W^{2,\infty}$ 意义下都是 结构不稳定的。

3. 技术框架的独立价值

  • 加权量 $A(t)$：作为追踪平移涡环芯的新工具，其抗误差累积的特性独立于具体问题，可能适用于其他具有相干平移结构的系统（如孤子、激波）。
  • 分解思想：将初始涡度分解为“核心”（$w_{m,0,\epsilon}$）和“慢部分”（$w_{d,0,\epsilon}$）的框架（定理 1.8），为分析多尺度相互作用提供了清晰构造。此分解稳健，可迁移到其他存在 相干运动与剪切竞争 的输运方程。

总结：本文通过发展 加权动量障碍 和 双速度几何分析，将涡环动力学研究从 短时、特殊数据 推进到 全局时间、通用初始数据，同时首次严格证明了 无对称性假设下单涡环的普遍丝化。其核心假设仅涉及宏观不变量，与经典变分行波构造或对称性假设的方法形成鲜明对比，代表了该领域方法论上的重要进展。

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

1. 论热渗与热分子压差

On thermal transpiration and thermomolecular pressure difference

作者: Kai-Li Wang, I-Kun Chen

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.21079v1

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

摘要翻译

中文标题: 论热渗与热分子压差

摘要翻译: 本文论证了在有界凸区域内的热渗现象。我们采用带截断势的稳态玻尔兹曼方程。对于边界条件，我们将边界划分为漫反射区与入射区。我们在加权L^∞空间中建立了解的存在性。此外，我们又考虑了具有中间漫反射边界条件及两端入射边界条件的凸区域。 我们首先考虑两端具有相同压强但不同温度的麦克斯韦分布。我们证明总流U(x)指向热端。此外，我们得到了总流的估计： \begin{align} U(x)\geq C\left(1-\frac{1}{\sqrt{T_2}}\right). \end{align} 进一步，我们证明当两端的压强与温度满足关系 \begin{align} \frac{P_1}{P_2}=\sqrt{\frac{T_1}{T_2}}, \end{align} 时，解的总流阶为$\mathcal{O}(\frac{1}κ)$。该结果与Knudsen在1909年关于热分子压差的发现相一致。

2. 一维组合非线性薛定谔方程的爆破准则

Blow-up criteria for the 1D NLS with combined nonlinearities

作者: Alex D Rodriguez

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.21040v1

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

摘要翻译

中文标题: 一维组合非线性薛定谔方程的爆破准则

摘要翻译: 本文针对一维组合非线性薛定谔方程的有限时间爆破解，建立了两类不同的准则。第一类为针对三重组合非线性以及包含指数型非线性的（收敛）无穷级数组合非线性情形所发展的负能量准则。第二类为针对具有散焦-聚焦系数的双重非线性问题建立的正能量准则。同时，本文也给出了满足上述任一类准则的初值示例。

3. 圆周上半波映射方程：在 $H^{1/2}$ 中的全局适定性与概周期性

The half-wave maps equation on $\mathbb{T}$: Global well-posedness in $H^{1/2}$ and almost periodicity

作者: Patrick Gérard, Enno Lenzmann

类别: math.AP, math-ph, nlin.SI

发布日期: 2026-02-24

arXiv ID: 2602.20895v1

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

摘要翻译

中文标题: 圆周上半波映射方程：在 $H^{1/2}$ 中的全局适定性与概周期性

摘要翻译: 我们考虑定义在 $\mathbb{R} \times \mathbb{T}$ 上、取值于 $\mathbb{S}^2$ 的半波映射方程

$$ \partial_t \mathbf{u} = \mathbf{u} \times |D| \mathbf{u} $$

，其中 $\mathbb{T}=\mathbb{R}/2 π\mathbb{Z}$ 为一维圆周，$\mathbb{S}^2 \subset \mathbb{R}^3$ 表示单位球面。通过从有理初始数据进行延拓，我们在临界能量空间 $H^{1/2}(\mathbb{T}; \mathbb{S}^2)$ 中构造了唯一且连续的流映射。此外，我们证明了这些解在时间上具有概周期性。对于稠密的有理初始数据子集，我们建立了时间上的拟周期性，并对任意 $s >0$ 建立了 $\| \mathbf{u}(t) \|_{H^s(\mathbb{T})}$ 的先验估计。 我们的分析关键在于利用一个显式公式，该公式源于作用在单位圆盘上向量值全纯函数Hardy空间上的Lax对结构。作为核心工具，我们发展了一个针对可积偏微分方程（在Hardy空间上具有Lax对结构）显式公式的通用稳定性原理，该原理适用于Benjamin–Ono方程、Calogero–Sutherland DNLS以及定义在 $\mathbb{T}$ 上的半波映射方程。 我们的结果推广到矩阵值半波映射方程

$$ \partial_t \mathbf{U} = -\frac{i}{2} [ \mathbf{U}, |D| \mathbf{U} ] $$

，其目标流形为复格拉斯曼流形 $\mathsf{Gr}_k(\mathbb{C}^d)$，从而将特例 $\mathbb{S}^2 \cong \mathbb{CP}^1 \cong \mathsf{Gr}_1(\mathbb{C}^2)$ 进行了推广。在配套研究中，我们通过在上半平面 $\mathbb{C}_+$ 的Hardy空间上为显式公式建立稳定性原理，证明了半波映射方程在 $\mathbb{R}$ 上于标度临界能量空间 $\dot{H}^{1/2}$ 中的全局适定性。

4. 逆随机源问题的赫尔德-对数稳定性与收敛速率

Hölder-Logarithmic Stability and Convergence Rates for an Inverse Random Source Problem

作者: Philipp Mickan, Thorsten Hohage

类别: math.NA, math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20822v1

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

摘要翻译

中文标题: 逆随机源问题的赫尔德-对数稳定性与收敛速率

摘要翻译: 本文研究一个逆随机源问题，旨在根据发射出的时谐声波的相关性测量值，重建一个随机、非相关声源的强度。此类问题在气动声学和地震成像等领域具有实际应用背景。与确定性逆问题不同，已知在无噪声情形下逆随机源问题具有唯一可解性。然而，由于其内在不适定性，需采用正则化方法以实现源强度的稳定重建。
在索博列夫光滑性假设下，本文利用复几何光学解推导出条件赫尔德-对数稳定性估计。进一步，通过建立变分源条件，获得了谱正则化方法的赫尔德-对数收敛速率。在固定频率时，随着源函数索博列夫正则性的提高，对数稳定性和收敛估计中的指数将无限增长。最后，本文通过数值实验验证了理论结论。

5. 具有线性与拟线性扩散的Hardy-Hénon方程的全局自相似解

Global self-similar solutions for Hardy-Hénon equations with linear and quasilinear diffusion

作者: Razvan Gabriel Iagar, Ariel Sánchez, Erik Sarrion-Pedralva

类别: math.AP, math.DS

发布日期: 2026-02-24

arXiv ID: 2602.20699v1

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

摘要翻译

中文标题: 具有线性与拟线性扩散的Hardy-Hénon方程的全局自相似解

摘要翻译: 在指数范围 $m\geq1$，$p>m$ 且 $\sigma>\max\{-2,-N\}$ 内，对抛物型Hardy-Hénon方程

$$ u_t=Δu^m+|x|^σu^p, \quad (x,t)\in\mathbb{R}^N\times(0,\infty) $$

的全局自相似解进行了分类。该分类强烈依赖于著名的Fujita临界指数与Sobolev临界指数：

$$ p_F(σ)=m+\frac{σ+2}{N}, \quad p_S(σ)= \begin{cases} \frac{m(N+2σ+2)}{N-2}, & \mbox{若 } N\geq3, \\[1mm] \infty, & \mbox{若 } N\in\{1,2\}. \end{cases} $$

具体而言，当 $p\in(p_F(σ),p_S(σ))$ 时，两类方程均存在自相似解，其既有紧支集（若 $m>1$），也有当 $|x|\to\infty$ 时呈类高斯衰减的解（若 $m=1$），还参有一族满足渐近行为

$$ u(x,t)\sim C|x|^{-(σ+2)/(p-m)}, \quad \text{当 } |x|\to\infty $$

的单参数解族。若 $p\geq p_S(σ)$，则仅存在具有上述代数衰减尾部的自相似解；而当 $m<p\leq p_F(σ)$ 时，全局解不存在。这些结果为深入研究此类解在Hardy-Hénon方程动力学中所起的作用奠定了基础。

6. 具有一般势函数的二维Gross-Pitaevskii模型中沿一维曲线的浓度解

Solutions with one dimensional concentration for a two dimensional Gross-Pitaevskii model with general potential

作者: Lipeng Duan, Suting Wei, Jun Yang

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20554v1

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

摘要翻译

中文标题: 具有一般势函数的二维Gross-Pitaevskii模型中沿一维曲线的浓度解

摘要翻译: 本文研究在单位质量约束下，具有阱势的二维Gross-Pitaevskii方程（用于描述具有吸引相互作用的玻色-爱因斯坦凝聚体）的频率为$λ$的驻波解。首先，我们探讨了浓度现象沿封闭光滑曲线方向产生的解存在的必要条件。其次，在曲线关于涉及阱势的辅助加权长度满足平稳性和非退化条件的基础上，并补充若干技术性假设，我们选取频率$λ$的一个序列$\{λ_j\}$（满足$-λ_j\rightarrow +\infty$），构造了浓度沿曲线方向集中的解。该结果部分回答了[A. Ambrosetti, A. Malchiodi, W.-M. Ni, Comm. Math. Phys. 2003]中关于解在子流形上集中的必要条件所提出的猜想。本文构造的解所集中的曲线长度是非一致有界的，因此与[M. del Pino, M. Kowalczyk, J. Wei, Comm. Pure Appl. Math. 2007]中的情形存在显著差异。

不相关论文（基本信息）

1. 一般开集上p-拉普拉斯算子的低特征值

Low eigenvalues of the $p-$Laplacian in general open sets

作者: Lorenzo Brasco, Luca Briani, Francesca Prinari

类别: math.AP, math.SP

发布日期: 2026-02-24

arXiv ID: 2602.21118v1

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

摘要: We consider the minmax Ljusternik-Schnirelmann levels of the constrained $p-$Dirichlet integral, on a general open set of the Euclidean space. We show that, whenever one of these levels lies below the threshold given by the $L^p$ Poincaré constant ``at infinity’’, it actually defines an eigenvalue of the Dirichlet $p-$Laplacian. We also prove an exponential decay at infinity for the relevant eigenfunctions: this can be seen as a Šnol-Simon–type estimate for the nonlinear case. Finally, we exhibit some peculiar examples of unbounded open sets to which our main result applies.

2. MEMS问题中半线性椭圆方程稳定径向解的正则性

Regularity of stable radial solutions to semilinear elliptic equations in MEMS problems

作者: Fa Peng, Salvador Villegas

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.21115v1

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

摘要: This paper investigates the regularity of stable radial solutions to semilinear elliptic equations arising in MEMS problems, modeled by the Dirichlet problem $-Δu=f(u)$ in the unit ball $B_1$, where the nonlinearity $f\in C^1([0,1))$ is nonnegative and satisfies $\int^1_0f(s)\,ds=+\infty$. We focus on the case where $f$ blows up as $u\to 1^{-}$. Micro-electro-mechanical systems (MEMS) are widely used devices in engineering and technology. Our main result establishes for dimensions $2\le n\le 6$, every stable radial solution is regular, meaning $\|u\|_{L^{\infty}(B_1)}<1$. This result gives a positive answer to an open problem posed by Bruera and Cabré concerning the regularity of stable solutions for singular nonlinearities without requiring a Crandall-Rabinowitz type condition, at least in the radial case.

3. 狭长凸域中的边界驱动模式

Boundary-driven patterns in elongated convex domains

作者: Maicon Sonego

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20938v1

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

摘要: We consider the heat equation in a smooth bounded convex domain $Ω\subset \mathbb{R}^2$ with nonlinear Neumann boundary condition $\partial_νu = λ(u - u^3)$. Stable non-constant stationary solutions do not exist when $Ω$ is a ball. We show that this behavior is not a consequence of convexity alone. More precisely, if the inradius of $Ω$ is fixed and its diameter is sufficiently large, then there exists $λ>0$ for which the problem admits such a solution. The result reveals a geometric mechanism for the emergence of stable non-constant stationary solutions in elongated convex domains.

4. 高余维带自由边界的曲线缩短流

High Codimension Curve Shortening Flow with Free Boundary

作者: Huy The Nguyen, Artemis A. Vogiatzi

类别: math.DG, math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20865v1

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

摘要: We study curve shortening flow in high codimension for arcs with free boundary meeting a fixed smooth barrier orthogonally. We prove dilation-invariant curvature and higher-derivative estimates up to the boundary using a Stahl-type localised maximum principle and an adapted cut-off. Using a reflected Gaussian entropy and blow-up analysis, Type I boundary singularities yield a shrinking semicircle model after reflection. Type II blow-ups give a Grim Reaper translator, which is ruled out under a free-boundary entropy bound $<2$. Hence in the low-entropy regime the flow either converges to the orthogonal chord or has only semicircle boundary singularities.

5. 巴拿赫空间上强不定泛函的莫尔斯同调

Morse homology for strongly indefinite functionals on Banach spaces

作者: L. Asselle, S. Cingolani, M. Starostk

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20863v1

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

摘要: In this paper we lay the foundations for the Morse theoretical study of strongly indefinite functionals on Banach manifolds by developing the local theory for a specific model class that captures several key analytical features also arising in the variational formulations of geometric problems such as Dirac-harmonic maps. As a corollary, we obtain existence results of solutions to certain systems of quasilinear elliptic problems involving the $p$-area functional. Abstracting from the concrete setting, we then formulate general conditions ensuring that Morse homology is well-defined for strongly indefinite functionals on a Banach space.

6. 一类非光滑泛函的多重临界点定理及其在1-拉普拉斯算子和不连续非线性项驱动问题中的应用

Multiple critical points theorems for a class of nonsmooth functionals and applications to problems driven by 1-Laplacian and discontinuous nonlinearities

作者: Ismael Sandro da Silva, Marcos T. Oliveira Pimenta, Pedro Fellype Pontes

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20858v1

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

摘要: In this paper, we present a novel approach to investigate the existence of multiple critical points for a class of nonsmooth functionals. This method provides a robust framework to analyze the existence of solutions for problems involving the $1$-Laplacian operator with discontinuous nonlinearities. Our results contribute to advancing the study of nonsmooth variational problems, by establishing new nonsmooth multiple critical point theorems.

7. 爱因斯坦-波多尔斯基-罗森论证的数学模型

A mathematical model for the Einstein-Podolsky-Rosen argument

作者: Riccardo Adami, Luigi Barletti, Alessandro Teta

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

发布日期: 2026-02-24

arXiv ID: 2602.20827v1

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

摘要: We study a nonrelativistic system made of two quantum particles constrained to move on a line and a spin located at a fixed point of the line. Initially the two particles are in a maximally entangled state and the spin is down. The first particle interacts with the spin while the second particle is free, i.e., it does not interact neither with the first particle nor with the spin. We rigorously prove that there is a correlation between the state of the spin and the state of the second particle. More precisely, we show that, in a suitable scaling limit, if the first particle flips the spin, then the second particle possesses a definite momentum in the direction opposite to the spin.

8. 半线性抛物方程的连续数据同化：基于演化方程的一般方法

Continuous Data Assimilation for Semilinear Parabolic Equations: A General Approach by Evolution Equations

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

类别: math.AP, math-ph

发布日期: 2026-02-24

arXiv ID: 2602.20769v1

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

摘要: This article develops a general framework for continuous deterministic data assimilation for semilinear parabolic equations by means of evolution equations. Introducing a nudged model driven by partial observations, the global well-posedness of the reference and the approximating systems is established under natural assumptions. In addition, it is shown that the approximating solution converges exponentially to the solution of the reference system, provided the observational resolution and the nudging parameter are suitably chosen. The approach allows us to consider many systems, such as the Allen-Cahn, Cahn-Hilliard, Sellers-type energy balance, and bidomain systems, for the first time.

9. Toda系统的流方法

A flow approach to the Toda system

作者: Yong Luo, Linlin Sun, Guofang Wang

类别: math.DG, math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20707v1

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

摘要: In this paper we introduce a flow to study the Toda system, which we call {\it Toda flow.} More generally, we introduce a flow of the Liouville systems, formulated as a coupled parabolic system with nonlocal interactions. Finite-time singularities are characterized and both necessary and sufficient conditions for convergence are provided in this general setting, even when the prescribed functions are allowed to change sign. As an application, we prove a global existence for the Toda flow in the critical case without restricting the sign of the prescribed functions. We provide a detailed description of blow-up behavior at infinity and obtain a sharp lower bound for the functional in cases where global convergence fails. By constructing appropriate test functions, we further establish a sufficient condition for the global convergence of the flow. These results are not affected by the sign-changing nature of the prescribed functions, and extend the theorem of Jost, Lin and Wang (Comm. Pure Appl. Math. 59, 526-558, 2006) to systems of multiple equations under this more general and physically relevant condition.

10. 基于梯度敏感性分析的超弹性介质降阶计算均质化方法

Reduced-order computational homogenization for hyperelastic media using gradient based sensitivity analysis of microstructures

作者: Vladimír Lukeš, Eduard Rohan

类别: math.NA, math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20697v1

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

摘要: We propose an algorithm for the computational homogenization of locally periodic hyperelastic structures undergoing large deformations due to external quasi-static loading. The algorithm performs clustering of macroscopic deformations into subsets called “centroids”, and, as a new ingredient, approximates the homogenized coefficients using sensitivity analysis of micro-configurations with respect to the macroscopic deformation. The novel “model-order reduction” approach significantly reduces the number of microscopic problems that must be solved in nonlinear simulations, thereby accelerating the overall computational process. The degree of reduction can be controlled by a user-defined error tolerance parameter. The algorithm is implemented in the finite element framework SfePy, and its performance effectiveness is demonstrated using two-dimensional test examples, when compared with solutions obtained by the proper orthogonal decomposition method, and by the full “FE-square” simulations. Extensions beyond the present implementations and the scope of tractable problems are discussed.

11. 基于各向同性松弛微极连续体的长圆柱壳径向变形解析解

Analytical solutions for long cylindrical shells under radial deformations based on the isotropic relaxed micromorphic continuum

作者: Esmaeal Ghavanloo, Pierre Fritsch, Patrizio Neff

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20675v1

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

摘要: This study presents a closed-form analytical solution for the elastostatic response of long cylindrical shells composed of microstructured materials within the framework of the isotropic relaxed micromorphic continuum. The formulation accounts for microstructural effects by introducing an independent micro-distortion tensor field in addition to the classical displacement field. Under the assumptions of axisymmetric deformation and plane strain conditions, the governing equilibrium equations reduce to a coupled system of ordinary differential equations in the radial coordinate. By introducing suitable auxiliary variables, the system is reformulated into a non-homogeneous modified Bessel equation, which admits an exact analytical solution. Explicit expressions are derived for the radial displacement field and the non-zero components of the micro-distortion tensor. Numerical examples are presented to illustrate the influence of material parameters and the characteristic length on the displacement. The results demonstrate that the relaxed micromorphic model predicts deviations from classical elasticity where microstructural effects are more pronounced. The obtained solution provides valuable physical insight into the mechanics of cylindrical shells and serves as a benchmark for validating numerical implementations of relaxed micromorphic models.

12. 状态依赖时滞与快速增长非线性项反应扩散方程的全局存在性

Global Existence for Reaction-diffusion Equations with State-Dependent Delay and Fast-growing Nonlinearities

作者: Ruijing Wang

类别: math.AP

发布日期: 2026-02-24

arXiv ID: 2602.20473v1

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

摘要: This work aims to study the initial-boundary value problem of the reaction-diffusion equation with state-dependent delay $\pa_{t}u-Δu=f(u)+g(u,u(t-τ(t,u_t)))+h(t,x)$ in a bounded domain. We establish the global existence of the problem under suitable dissipative-type structural conditions, allowing both nonlinear terms $f$ and $g$ to have arbitrary polynomial growth rates. Another highlight in this work is that, we significantly relax the continuity assumptions imposed on the delay functions.

13. 改进的分形 Weyl 界限匹配双曲曲面与开放量子映射中的改进谱隙

Improved fractal Weyl bounds matching improved spectral gaps for hyperbolic surfaces and open quantum maps

作者: Travis Cunningham

类别: math.SP, math.AP

发布日期: 2026-02-23

arXiv ID: 2602.20392v1

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

摘要: We prove a new fractal Weyl upper bound for the high-energy distribution of resonances of convex co-compact hyperbolic surfaces which matches the improved spectral gap given by Fourier decay. This improves upon the fractal Weyl bound of Dyatlov which matches the Patterson-Sullivan spectral gap. We also give a new resolvent estimate improving the ones given by Dyatlov-Zahl and Dyatlov. Analogous results are obtained for quantum open baker’s maps, improving an estimate of Dyatlov-Jin, where we also give an improved fractal Weyl bound matching a spectral gap given by additive energy estimates. We refine known methods for proving fractal Weyl bounds which reduce the problem to an estimate of a certain determinant function; however, we use a different determinant function which allows us to make sharper estimates by applying the methods of proof of the fractal uncertainty principle in each setting.

14. 线性弹性中脆性断裂的均匀可微性

Uniform rectifiability of brittle fractures in linear elasticity

作者: Camille Labourie

类别: math.AP

发布日期: 2026-02-23

arXiv ID: 2602.20381v1

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

摘要: We prove the uniform rectifiability of brittle fractures in arbitrary dimension. The existing approach for the Mumford-Shah functional, which relies on separation-type properties of the singular set, faces serious obstacles in the Griffith setting due to the lack of coarea formula for the symmetric gradient. We present an alternative route to uniform rectifiability for free-discontinuity problems by proving that cracks have ``plenty of big projections’'.

15. 球面上的新极小曲面：源自毛细极小锥

New minimal surfaces in the sphere from capillary minimal cones

作者: Benjy Firester, Raphael Tsiamis

类别: math.DG, math.AP

发布日期: 2026-02-23

arXiv ID: 2602.20124v1

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

摘要: For every $p,q\geq 1$, we construct minimal embeddings of $\mathbb{S}^p \times \mathbb{S}^q \times \mathbb{S}^1$ in $\mathbb{S}^{p + q + 2}$ by doubling the links of free-boundary minimal cones in $\mathbb{R}^{p+q+3}$ with bi-orthogonal symmetry. This solves problems posed by Hsiang-Lawson and Hsiang-Hsiang. The equivariance reduces the minimal surface equation to an ODE, and we prove the existence of capillary minimal cones for every contact angle. We obtain free-boundary solutions as limits of capillary surfaces via a singular shooting problem with infinite initial slope. As the contact angle degenerates to $0$, rescalings of the capillary cones converge to a homogeneous solution of the one-phase Bernoulli problem, further illustrating the connection between one-phase free boundaries and minimal surfaces through the capillary functional.

16. 爱因斯坦-标量场共形约束方程的球对称解

Spherically symmetric solutions to the Einstein-scalar field conformal constraint equations

作者: Philippe Castillon, The-Cang Nguyen

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

发布日期: 2026-02-23

arXiv ID: 2602.20099v1

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

摘要: Recent works by the second author and Dilts et al. have shown that the Einstein-scalar field conformal constraint equations are highly complex and generally intractable, even in the vacuum case. In this article, to gain a clearer understanding and offer a new perspective, we study these equations under special assumptions: the manifold $(M,g)$ is harmonic and data is radial. In this setting, the system reduces to a single nonlinear equation and is completely resolved in the standard cases. In particular, on the sphere, our results reveal phenomena that contrast with the well-known achievements on compact manifolds without conformal Killing vector fields, including nonexistence of solutions in the near-CMC regime and instability when the mean curvature is non-constant. By contrast, on Euclidean or hyperbolic manifolds, the equations are always solvable, with all expected properties of solutions satisfied. These findings support the view that, although the conformal method appears to present some drawbacks on compact manifolds, it remains a promising tool for parametrizing solutions to the constraint equations on asymptotically flat and hyperbolic manifolds in arbitrary mean curvature regimes. In this article, we also investigate the sign of mass, showing that the ADM and asymptotically hyperbolic mass can take arbitrary sign when the decay rate of symmetric $(0,2)$-tensor $k$ at infinity is critical. Finally, most solution classes in our framework are explicit, providing a variety of models in general relativity and offering insights into the structure and behavior of initial data, particularly in numerical applications.

17. 三维空间中2-蒙日-安培方程的刘维尔型定理

A Liouville-type theorem for $2$-Monge-Ampère equation in dimension three

作者: Weisong Dong

类别: math.AP

发布日期: 2026-02-23

arXiv ID: 2602.20090v1

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

摘要: We prove that every entire solution with quadratic growth, lying in a suitable cone, to the 2-Monge-Ampère equation on $\mathbb{R}^3$ is a quadratic polynomial. The proof proceeds by first establishing a concavity inequality, and then deriving a Pogorelov-type interior $C^2$ estimate.