生成时间: 2026年06月19日 08:24:14

运行状态: 完成

计划处理论文数量: 17

已完成论文数量: 17

未完成论文数量: 0

重点关注论文数量: 4

了解领域论文数量: 5

不相关论文数量: 8

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

1. 通道中涡旋块的稳定性

Stability of Vortex Patches in Channels

作者: Zelin Dong, Chenyun Luo

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.19942v1

摘要翻译: 本文研究二维不可压 Euler 方程在满足“弱有限体积条件”的一类区域以及任意宽度条带区域中涡旋块的轨道稳定性。我们证明：对合适的参数 $(μ,λ)$，带罚项的动能泛函存在极小元；且每个此类极小元均满足椭圆方程 $ω= λ(ψ- W x_2 - γ)_+$。此外，我们证明该极小元集合在 Euler 演化下是轨道稳定的。本工作将 Abe 与 Choi 发展的变分框架推广至既不具有空间尺度不变性、也不具有水平平移不变性的区域。这些不变性的缺失为证明带来了实质性困难，因为经典的重排与尺度论证不再适用。我们通过将 Green 函数与半平面 Green 函数进行比较，并利用衰减条件构造一个集中-紧性论证，最终得到所期望的稳定性结果。

Comment: 23 pages

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

1. 研究对象和背景

本文研究二维不可压 Euler 方程

$$ \partial_t \omega + u \cdot \nabla \omega = 0, \quad \text{in } \Omega \times (0,T), $$

在两类非齐性无界区域上的涡旋块（vortex patch）轨道稳定性问题，其中 $\omega$ 是标量涡度，速度场 $u = \nabla^\perp \psi$ 由流函数 $\psi$ 通过 Biot–Savart 律确定：

$$ \psi(x) = \int_\Omega G(x,y) \, \omega(y) \, dy, $$

且 $G(x,y)$ 是区域 $\Omega$ 上 Dirichlet-Laplace 算子 $-\Delta$ 的 Green 函数，满足 $u \cdot N = 0$（滑移边界条件）。研究对象是满足 $\omega = \lambda(\psi - W x_2 - \gamma)_+$ 形式的特征解——即支撑集为水平剪切背景 $W x_2$ 下的截断流函数水平集——其物理图像为在通道中沿 $x_1$ 方向平移不变性被破坏后仍能维持形状的集中涡结构。本文聚焦的两类区域为：(i) 满足“弱有限体积条件”的区域 $D = D_u \cup D_l$，其中 $D_u = D \cap \{x_2 \ge 1\}$ 有有限测度，$D_l = D \cap \{x_2 < 1\}$ 满足 $\int_{D_l} x_2^p \, dx < \infty$ 对某个 $p \ge 2$ 成立；(ii) 无限条带 $S = \mathbb{R} \times (0,L)$，其 Green 函数显式为

$$ G(x,y) = \frac{1}{4\pi} \ln \left( 1 + \frac{\sin(\frac{\pi}{L}x_2) \sin(\frac{\pi}{L}y_2)}{\sinh^2(\frac{\pi}{2L}(x_1 - y_1)) + \sin^2(\frac{\pi}{2L}(x_2 - y_2))} \right). $$

这两类区域均缺失两个关键对称性：空间尺度不变性（即 $G(cx,cy) \ne G(x,y)$）与水平平移不变性（即 $G(x,y)$ 不能仅依赖于 $|x_1 - y_1|$）。这一缺失使得经典变分框架（如 Abe–Choi [1] 在上半平面 $\mathbb{R}^2_+$ 的工作）完全失效：其证明严重依赖于 $G$ 的平移不变性以施行重排（rearrangement），以及尺度不变性以建立能量泛函 $I_{\mu,\nu,\lambda} = \inf_{\omega \in K_{\mu,\nu}} \{ -E_{2,\lambda}[\omega] \}$ 关于 $\mu$ 的严格单调递减性（从而控制极小化序列的集中行为）。而在 $D$ 和 $S$ 中，重排不保能量，尺度变换扭曲动能与惩罚项的平衡（见原文第3页计算：$E^{(c)}_{2,\lambda}[\tilde{\omega}] = c^{-4} E[\omega] - \frac{1}{2}\lambda c^{-2} \|\omega\|_2^2$，无法令 $c^{-4}=c^{-2}$ 除非 $c=1$），且 Robin 函数、对称极小点等几何工具亦不可用。因此，已有理论在非齐性通道中存在根本性空白：既无适用于一般弱有限体积区域的变分存在性定理，也无条带中不依赖平移/尺度对称性的轨道稳定性机制。本文旨在填补此空白，构建一套摆脱重排与尺度论证、仅依赖 Green 函数比较与加权集中紧性（concentration-compactness）的新分析框架。

2. 主要定理或主要结果

在参数 $(\mu, \lambda)$ 满足

$$ (\mu, \lambda) \in \bigcup_{\substack{K \subseteq \Omega \\ |K| < \infty}} (0, \mu_K] \times (\lambda_K, \infty), \quad \text{其中 } \mu_K = \frac{1}{|K|} \int_K x_2 \, dx, \; \lambda_K = \frac{|K|}{\iint_{K \times K} G(x,y) \, dydx}, $$

的条件下，本文建立了如下核心结论：(i) 【存在性】对任意满足弱有限体积条件的区域 $D$（Class 1），变分问题 $I_{\mu,\lambda} = \inf_{\omega \in K_\mu} \{ -E_{2,\lambda}[\omega] \}$ 存在极小元，即集合 $S_{\mu,\lambda} = \{ \omega \in K_\mu : -E_{2,\lambda}[\omega] = I_{\mu,\lambda} \}$ 非空（Theorem 3.2）；(ii) 【结构刻画】每个极小元 $\omega \in S_{\mu,\lambda}$ 必满足椭圆自由边界方程

$$ \omega = \lambda(\psi - W x_2 - \gamma)_+, \quad \psi(x) = \int_D G(x,y) \, \omega(y) \, dy, $$

其中常数 $W, \gamma \in \mathbb{R}$ 由 $\omega$ 唯一确定（Proposition 2.6）；(iii) 【条带中的存在性与单调性】对无限条带 $S = \mathbb{R} \times (0,L)$（Class 2），在 $\mu$ 充分小时，$I_{\mu,\lambda}$ 关于 $\mu$ 严格单调递减（Proposition 4.12），并由此导出极小元存在性及一般收敛定理（Theorem 4.14）；(iv) 【轨道稳定性】集合 $S_{\mu,\lambda}$ 在 Euler 动力学下轨道稳定：对任意 $\varepsilon > 0$，存在 $\delta > 0$，使得若初值 $\omega_0 \in L^2(\Omega)$ 满足 $\inf_{\tilde{\omega} \in S_{\mu,\lambda}} \|\omega_0 - \tilde{\omega}\|_{L^2} < \delta$，则其对应解 $\omega(t)$ 满足 $\sup_{t \in [0,T]} \inf_{\tilde{\omega} \in S_{\mu,\lambda}} \|\omega(t) - \tilde{\omega}\|_{L^2} < \varepsilon$ 对所有有限时间 $T < \infty$ 成立（Theorem 5.1）。该稳定性是 $L^2$-范数意义下的，并未声称全局时间稳定性，但覆盖任意有限时间区间。

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

证明的核心障碍在于：缺失平移与尺度不变性导致三重技术断裂——重排失效、能量单调性消失、Green 函数缺乏显式可控性。作者的真正创新点在于构建了一个不依赖对称性的、基于 Green 函数比较与加权集中紧性的新框架。具体包含以下关键技术：(1) 【Green 函数上界比较】利用 $\Omega \subseteq \mathbb{R}^2_+$ 及 Dirichlet 边界条件，获得关键不等式 $0 \le G(x,y) \le G_H(x,y)$，其中 $G_H$ 是上半平面 Green 函数。此比较使所有动能估计得以移植：例如，由 $\int_{\mathbb{R}^2_+} G_H(x,y)^q \, dy \le C x_2^{q/2}$（$q \in (1,2)$）导出 Proposition 2.1 中的加权 $L^1$-$L^2$ 混合估计，如

$$ |E[\omega]| \le C \|x_2 \omega\|_1^{1/2} \|\omega\|_1 \|\omega\|_2^{1/2}, $$

此类估计是控制极小化序列 $L^2$ 一致有界性（Remark 2.4）与误差项收敛性的基石。(2) 【弱有限体积条件驱动的加权集中】对 Class 1 区域 $D$，利用 $\int_{D_l \cap \{|x_1| > k\}} x_2^p \, dx \to 0$（$k \to \infty$）及 Hölder 不等式，证明对任意极小化序列 $\{\omega_n\}$ 和 $\varepsilon > 0$，存在 $N$ 使得 $\int_{D_l \cap \{|x_1| > N\}} x_2 \omega_n \, dx < \varepsilon$（Proposition 3.1）。这迫使加权质量 $x_2 \omega_n$ 在 $x_1$ 方向集中于有界带 $D_0 = D_u \cup (D_l \cap \{|x_1| \le N\})$，而 $G \in L^2(D_0 \times D_0)$，从而在 $D_0$ 上可对弱 $L^2$ 收敛序列 $\omega_n \rightharpoonup \omega$ 通过紧性传递动能极限 $E[\omega_n] \to E[\omega]$。(3) 【条带中的重排重构与单调性重建】对 Class 2 条带 $S$，虽无平移不变性，但其 Green 函数结构满足标准重排不等式前提（Proposition 4.5），且通过嵌入 $L^i(S) \hookrightarrow L^i(\mathbb{R}^2_+)$ 并利用 $G \le G_H$，构造出上半平面流函数 $\psi_H(x) = \int_S G_H(x,y) \omega(y) \, dy$，满足 $\psi(x) \le \psi_H(x)$（1.10）。这允许将条带问题“提升”至上半平面，在后者中运用重排获得先验控制，再“拉回”条带；而 $I_{\mu,\lambda}$ 的严格单调性（Proposition 4.12）则放弃尺度论证，转而直接证明当 $\mu$ 充分小时，极小元具有一致 $L^1$ 衰减（Proposition 4.9），从而在集中-紧性框架中分离出紧支集极限。（4) 【集中-紧性框架的定制化实现】全文统一采用 Lions 型集中-紧性，但关键在于“集中”判据的加权化：不是传统 $\|\omega_n\|_{L^1(B_R)}$，而是 $\int_{B_R} x_2 \omega_n \, dx$，其权重 $x_2$ 恰好匹配第一矩约束 $\int x_2 \omega = \mu$，确保集中位置物理意义明确（靠近 $x_2$ 小的区域），且与条带高度 $L$ 或 $D$ 的几何衰减相容。

4. 与之前工作的比较

本文与既有工作的本质区别体现在四个维度：(i) 【几何适用性】Abe–Choi [1] 仅适用于上半平面 $\mathbb{R}^2_+$，其证明完全依赖 $G_H(x,y) = G_H(x-y)$ 的平移不变性与 $G_H(cx,cy) = G_H(x,y)$ 的尺度不变性；Choi–Jeong–Lim [7] 限于半柱体（half-cylinder），Turkington [22] 依赖重排与凸域假设；而本文首次将变分稳定性推广至两类无任何平移/尺度对称性的区域：弱有限体积域（含所有有界域、周期边界域、图域）与任意宽度条带 $\mathbb{R} \times (0,L)$。代价是放弃了显式解构造，转而依赖抽象存在性。(ii) 【技术路径】此前所有工作（[1],[5],[7],[8]）均以重排为核心工具，通过将 $\omega$ 重排为关于某轴对称的递减函数来控制能量；本文彻底摒弃重排，代之以 Green 函数上界比较与加权集中，这是方法论的根本转向。例如，[1] 中 $I_{\mu,\lambda}$ 的单调性源于尺度变换 $\omega_c(x) = c^2 \omega(cx)$ 下 $E[\omega_c] = E[\omega]$ 与 $\|\omega_c\|_2^2 = c^2 \|\omega\|_2^2$ 的平衡，而本文（Proposition 4.12）必须在固定尺度下，通过精细估计极小元 $L^1$ 范数的 $\mu$-依赖性来重建单调性。(iii) 【时间尺度与稳定性强度】Nadirashvili–Šverák [18] 等揭示了 Euler 方程的普遍不稳定性（双指数增长），本文的轨道稳定性仅保证有限时间 $T < \infty$ 内的 $L^2$ 接近性，未宣称长期或渐近稳定性，这与 Cao–Wang [6] 在非退化 Robin 函数极小点附近的非线性稳定性（要求初值高度集中）形成互补：本文稳定性不依赖初值是否集中，而依赖其是否接近某个极小元集合，适用范围更广但时间尺度受限。(iv) 【假设松弛与新增限制】本文放宽了对区域对称性的要求，但引入了新的几何约束：“弱有限体积条件”（对 $D$）与条带宽度 $L$ 的任意性（对 $S$）；同时，参数 $(\mu,\lambda)$ 的允许范围（2.7）显式依赖于区域几何（如 $\lambda_K$ 依赖 $\iint_{K\times K} G$），而 [1] 中 $\lambda$ 可任取正数。从提供内容中无法确认本文结果能否推广至三维或粘性情形，亦未显式展开与近期基于流形稳定性（如 [17]）或准地转模型工作的比较。

2. Bregman投影镜像方法用于正则化定常平均场博弈

Bregman-projected mirror methods for regularized stationary mean-field games

作者: Hussain Al Abdulaziz, Yuri Ashrafyan, Yeva Gevorgyan, Diogo Gomes

类别: math.NA, math.AP

发布日期: 2026-06-17

arXiv ID: 2606.19611v1

摘要翻译: 我们在其自然的Banach空间框架下，发展并分析了一类针对定常平均场博弈（MFG）系统的低阶正则化的Bregman投影镜像迭代。对于形如 $H(x,p,m)=H_0(x,p)-g(m)$ 的可分离Hamiltonian，其中Hamiltonian具有二次或超二次增长性、密度耦合项具有线性或超线性增长性，我们将定常MFG系统的低阶 $\barγ$-Laplacian 正则化表述为定义在 $L^{\barβ}(\mathbb T^d)\times W^{1,\barγ}(\mathbb T^d)$ 上的一个变分不等式。为逼近该正则化变分不等式的解，我们引入一种与问题所涉混合Lebesgue–Sobolev指数相匹配的Bregman几何，并分析一类带冻结算子求值的约束型两步镜像方法。对于精确约束迭代及每个固定的正则化参数 $\epsi>0$，我们导出一个单步Bregman不等式，并据此证明：在步长满足自然可和性条件的前提下，该约束迭代强收敛于正则化变分不等式的唯一解。我们在一维与二维模型上进行的数值实验，通过与精确测试解对比验证了其有效性，展示了残差随网格加密而衰减的行为，并在所测试的离散格式中表明两步实现具有更优的实际性能。

Comment: 无

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

1. 研究对象和背景

本文研究的是定常平均场博弈（Stationary Mean-Field Games, MFG）系统的数值求解问题，其核心数学对象是如下耦合非线性偏微分方程系统：

$$ \begin{cases} u + H_0(x, Du) \le g(m) + V(x), \\ u + H_0(x, Du) = g(m) + V(x) \quad \text{a.e. on } \{m > 0\}, \\ -\operatorname{div}\big(m D_p H_0(x, Du)\big) + m = 1. \end{cases} \tag{1.1} $$

其中 $x \in \mathbb{T}^d$（$d$ 维环面），未知量为价值函数 $u : \mathbb{T}^d \to \mathbb{R}$ 和群体密度 $m : \mathbb{T}^d \to [0,\infty)$；$H_0 : \mathbb{T}^d \times \mathbb{R}^d \to \mathbb{R}$ 是凸 Hamiltonian，满足 $\alpha$ 阶增长：$\frac{1}{C}|p|^\alpha - C \le H_0(x,p) \le C(1+|p|^\alpha)$；$g : [0,\infty) \to \mathbb{R}$ 是单调非减耦合项，满足 $\beta$ 阶密度依赖：$|g(m)| \le C(1 + m^{\bar\beta-1})$，其中 $\bar\beta := \beta + 1$；势能 $V \in L^{\bar\beta'}(\mathbb{T}^d)$。该系统源于 Lasry–Lions 和 Huang–Malhamé–Caines 的奠基性工作，建模大量理性智能体在稳态下的纳什均衡行为：第一式为 Hamilton–Jacobi–Bellman 方程（HJB），刻画个体最优控制策略；第二式为 Fokker–Planck 型输运方程（FP），描述群体密度演化；二者通过 $m$ 与 $Du$ 非线性耦合，构成强非线性、非自伴、非变分结构的鞍点型系统。

由于缺乏经典椭圆/抛物结构，未正则化系统 (1.1) 在泛函分析层面存在严重困难：HJB 方程对 $u$ 缺乏强制性（coercivity），无法在 $W^{1,\bar\gamma}$ 中建立能量估计；而 FP 方程中 $m D_p H_0(x,Du)$ 项在 $L^{\bar\beta} \times W^{1,\bar\gamma}$ 框架下不满足标准单调算子理论所需的连续性与强制性配对。为此，Gomes 等人在 [10] 中提出低阶 $\bar\gamma$-Laplacian 正则化框架：对 $u$ 引入 $\varepsilon \big(-\Delta_{\bar\gamma} u + |u|^{\bar\gamma-2}u\big)$ 项，即定义正则化算子 $A_\varepsilon : K_+ \to X^*$，其中 $X = L^{\bar\beta}(\mathbb{T}^d) \times W^{1,\bar\gamma}(\mathbb{T}^d)$，$K_+ = \{(m,u) \in X : m \ge 0 \text{ a.e.}\}$，且

$$ \langle A_\varepsilon(m,u), (\mu,v) \rangle = \int_{\mathbb{T}^d} \big[ -u - H_0(x,Du) + g(m) + V \big] \mu ,dx

  • \int_{\mathbb{T}^d} m D_p H_0(x,Du) \cdot Dv ,dx
  • \int_{\mathbb{T}^d} (m-1)v ,dx
  • \varepsilon \int_{\mathbb{T}^d} |Du|^{\bar\gamma-2} Du \cdot Dv ,dx
  • \varepsilon \int_{\mathbb{T}^d} |u|^{\bar\gamma-2} u v ,dx. \tag{3.1} $$

关键参数由物理与对偶性要求严格确定：$\bar\gamma := \frac{\alpha \beta + 1}{\beta} = \alpha \bar\beta'$，以保证所有积分项良定——例如 $H_0(x,Du) \sim |Du|^\alpha \in L^{\bar\beta'}$，$m D_p H_0(x,Du) \sim m |Du|^{\alpha-1} \in L^{\bar\gamma'}(\mathbb{T}^d;\mathbb{R}^d)$，且满足恒等式 $\frac{\alpha}{\bar\gamma} = \frac{1}{\bar\beta'}$，$\frac{\alpha-1}{\bar\gamma} + \frac{1}{\bar\beta} = \frac{1}{\bar\gamma'}$。该框架的优势在于自然匹配 MFG 解的内在正则性：$m$ 属于 $L^{\bar\beta}$（由 $g(m) \sim m^{\bar\beta-1}$ 反推），$u$ 属于 $W^{1,\bar\gamma}$（由 $H_0$ 的 $\alpha$-增长与 $\bar\gamma$-Laplacian 正则化共同决定）。然而，现有数值方法均无法在此自然 Banach 空间中建立收敛性理论：Hilbert 空间方法（如 [2,7–9]）依赖内积投影，但 $L^{\bar\beta} \times W^{1,\bar\gamma}$ 非 Hilbert；离散或有限维镜像下降（如 [23,25]）作用于概率单形，几何失配；而 [10] 虽建立了存在性与 $\varepsilon \to 0$ 极限理论，却未提供与之兼容的迭代算法。本文旨在填补这一根本性空白：构造并严格证明一个在 $X = L^{\bar\beta} \times W^{1,\bar\gamma}$ 中强收敛的 Bregman 投影镜像迭代法，使数值方案的收敛性范数与 PDE 理论的先验估计范数完全一致。

2. 主要定理或主要结果

本文的核心结果是定理 1.1：在假设 2.3（结构条件）、2.4（定量单调性）和 2.5（多项式差分界）成立的前提下，对任意固定正则化参数 $\varepsilon > 0$，令 $K_+ = \{(m,u) \in X : m \ge 0 \text{ a.e.}\}$，并取截断集 $K = K_+ \cap B_X(0,R)$（$R > 0$ 充分大），初始点 $x_0 = (m_0,u_0) \in K$。定义 $\kappa = \max\{\bar\beta, \bar\gamma\}$，$\kappa' = \kappa/(\kappa-1)$。若步长序列 $(\lambda_n)$ 满足

$$ \sum_{n=1}^\infty \lambda_n = \infty, \quad \sum_{n=1}^\infty \lambda_n^{\kappa'} < \infty, $$

则由公式 (6.1) 定义的约束两步镜像方法在 $X = L^{\bar\beta}(\mathbb{T}^d) \times W^{1,\bar\gamma}(\mathbb{T}^d)$ 中强收敛到正则化变分不等式 $\mathrm{VI}(A_\varepsilon, K)$ 的唯一解 $(m_\varepsilon, u_\varepsilon)$。该结论的关键技术细节包括：(i) 收敛是强收敛（而非弱收敛或弱-*收敛），即 $\|m_n - m_\varepsilon\|_{L^{\bar\beta}} + \|u_n - u_\varepsilon\|_{W^{1,\bar\gamma}} \to 0$；(ii) 唯一性由定理 3.1 的定量单调性保证：存在 $c > 0$（与 $\varepsilon$ 无关）使得对任意 $(m,u), (\theta,v) \in K_+$，

$$ \langle A_\varepsilon(m,u) - A_\varepsilon(\theta,v), (m-\theta, u-v) \rangle \ge c \|m-\theta\|_{L^{\bar\beta}}^{\bar\beta} + c\varepsilon \|u-v\|_{W^{1,\bar\gamma}}^{\bar\gamma}; $$

(iii) 步长条件是经典的“发散但平方可和”型，但指数 $\kappa'$ 由空间混合指数 $\bar\beta, \bar\gamma$ 决定，而非统一取 2；(iv) 该结果不涉及 $\varepsilon \to 0$ 极限，而是对每个固定 $\varepsilon > 0$ 的精确迭代收敛性。此外，定理 3.4 给出了正则化极限的强化结果：若 $\varepsilon_k \to 0$ 且 $(m_{\varepsilon_k}, u_{\varepsilon_k}) \rightharpoonup (m,u)$ 弱于 $X$，则 $m_{\varepsilon_k} \to m$ 强于 $L^{\bar\beta}(\mathbb{T}^d)$，且 Hamiltonian Bregman 缺陷消失：

$$ \lim_{k\to\infty} \int_{\mathbb{T}^d} m \big[ H_0(x,Du_{\varepsilon_k}) - H_0(x,Du) - D_p H_0(x,Du) \cdot (Du_{\varepsilon_k} - Du) \big] \,dx = 0. $$

该强密度收敛性依赖于假设 2.4 的定量单调性 $ (g(m)-g(\theta))(m-\theta) \ge c_g |m-\theta|^{\bar\beta} $，是 [10] 中弱极限结果的重要加强。

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

本文的技术障碍根植于 Banach 空间 $X = L^{\bar\beta} \times W^{1,\bar\gamma}$ 的几何异质性：两个分量具有本质不同的可积性与微分阶数，导致标准 Hilbert 方法（如正交投影、内积诱导的镜像映射）完全失效。作者的真正创新点在于构造了一个与混合指数 $\bar\beta, \bar\gamma$ 精确匹配的 Bregman 几何，并围绕此几何设计了冻结算子评估的两步镜像迭代。具体而言，核心工具与技术分为三层：

第一，混合指数 Bregman 散度的构造（第 4 节）。作者定义标量生成元 $\psi_{\bar\beta}(s) = \frac{1}{\bar\beta} s^{\bar\beta}$（$s \ge 0$）和向量生成元 $\phi_{\bar\gamma}(\xi) = \frac{1}{\bar\gamma} |\xi|^{\bar\gamma}$（$\xi \in \mathbb{R}^d$），进而构建 $X$ 上的复合 Bregman 散度：

$$ D_\Psi\big((m,u), (\theta,v)\big) = \frac{1}{\bar\beta} |m|{L^{\bar\beta}}^{\bar\beta} - \frac{1}{\bar\beta} |\theta|{L^{\bar\beta}}^{\bar\beta} - \int_{\mathbb{T}^d} \theta^{\bar\beta-1} (m-\theta) ,dx

  • \frac{1}{\bar\gamma} |u|{W^{1,\bar\gamma}}^{\bar\gamma} - \frac{1}{\bar\gamma} |v|{W^{1,\bar\gamma}}^{\bar\gamma} - \langle \partial \phi_{\bar\gamma}(v), u-v \rangle_{W^{-1,\bar\gamma’}, W^{1,\bar\gamma}}. $$ 该散度的关键性质由引理 2.1 保证：利用点wise不等式 $(a^{r-1}-b^{r-1})(a-b) \ge |a-b|^r$（$r \ge 2$），可得 $D_\Psi\big((m,u), (\theta,v)\big) \ge \frac{1}{\bar\beta} \|m-\theta\|_{L^{\bar\beta}}^{\bar\beta} + c_{\bar\gamma} \|u-v\|_{W^{1,\bar\gamma}}^{\bar\gamma}$，这为后续收敛性提供了强制性距离度量。

第二，冻结算子评估误差的局部 Lipschitz 控制（第 3.4 节）。由于镜像迭代需在每步计算 $A_\varepsilon(x_n)$，而 $A_\varepsilon$ 本身是非线性的，作者引入假设 2.5（多项式差分界）：$|D_p H_0(p)-D_p H_0(q)| \le C(1+|p|^{\alpha-2}+|q|^{\alpha-2})|p-q|$ 和 $|g(m)-g(\theta)| \le C(1+m^{\bar\beta-2}+\theta^{\bar\beta-2})|m-\theta|$。结合假设 2.3，作者在任意有界集 $B_X(0,R)$ 上证明了 $A_\varepsilon$ 的局部 Lipschitz 性：$\|A_\varepsilon(m,u) - A_\varepsilon(\theta,v)\|_{X^*} \le C_{R,\varepsilon} \big( \|m-\theta\|_{L^{\bar\beta}} + \|u-v\|_{W^{1,\bar\gamma}} \big)$。该估计是处理“冻结”近似 $A_\varepsilon(x_n)$ 与真实算子 $A_\varepsilon(x)$ 差异的核心，确保迭代中累积误差可控。

第三，一步 Bregman 不等式与收敛框架（第 6 节）。作者对两步镜像迭代（先 $m$ 更新，再 $u$ 更新）导出关键的一阶不等式：存在常数 $C > 0$，使得对任意 $n$，

$$ D_\Psi\big(z_{n+1}, z^*\big) \le D_\Psi\big(z_n, z^*\big) - \lambda_n \langle A_\varepsilon(z_n), z_n - z^* \rangle + C \lambda_n^{\kappa'} \big( \|z_n - z^*\|_X^{\kappa} + \|z_n - z^*\|_X^{\kappa'} \big), $$

其中 $z^* = (m_\varepsilon, u_\varepsilon)$。该不等式将收敛性归结为两项：$-\lambda_n \langle A_\varepsilon(z_n), z_n - z^* \rangle$ 提供负耗散（由定量单调性 (3.4) 控制），而 $C \lambda_n^{\kappa'} \|z_n - z^*\|_X^{\kappa}$ 项代表冻结误差。在步长满足 $\sum \lambda_n = \infty$ 且 $\sum \lambda_n^{\kappa'} < \infty$ 的条件下，前者主导，后者可求和，从而通过标准递推论证得到 $D_\Psi(z_n, z^*) \to 0$，再由 Bregman 散度的强制性推出强收敛。

4. 与之前工作的比较

本文与既有工作的比较需从四个维度展开：函数空间框架、正则化类型、算法几何与收敛性保证。

首先，在函数空间框架上，本文严格区别于 Hilbert 空间方法（如 [2,7–9]）。这些工作将 MFG 置于 $L^2 \times H^1$ 或 $H^1 \times H^1$ 中，虽便于使用 Riesz 表示定理和正交投影，但严重违背解的自然正则性：当 $H_0(p) \sim |p|^\alpha$ 且 $g(m) \sim m^\beta$ 时，$m$ 的最优可积性是 $L^{\beta+1} = L^{\bar\beta}$，而非 $L^2$；$u$ 的最优 Sobolev 正则性是 $W^{1, \alpha(\beta+1)/(\beta+1)-\beta/\alpha} = W^{1,\bar\gamma}$，而非 $H^1$。因此，[2] 中的 $L^2$ 单调方法虽具启发性，但其收敛性范数与 [10] 的存在性理论不兼容，无法保证数值解在 PDE 理论关心的范数下收敛。本文则首次将收敛性证明置于 [10] 所确立的自然 Banach 空间 $L^{\bar\beta} \times W^{1,\bar\gamma}$ 中，实现了理论与算法的范数统一。

其次，在正则化类型上，本文采用低阶 $\bar\gamma$-Laplacian 正则化，区别于 [7–9] 的高阶正则化（如 $\Delta^2 u$）。高阶正则化虽在 Hilbert 框架下简化分析，但引入非物理的高阶导数项，损害解的物理可解释性，且其 $\varepsilon \to 0$ 极限的收敛速率通常更慢。而本文的低阶正则化直接匹配 Hamiltonian 的 $\alpha$-增长，其 $\varepsilon$-项 $\varepsilon \int |Du|^{\bar\gamma-2} Du \cdot Dv$ 在能量估计中贡献 $\varepsilon \|Du\|_{L^{\bar\gamma}}^{\bar\gamma}$，与 $u$ 的本征范数尺度一致，从而在定理 3.4 中获得了 $L^{\bar\beta}$ 中的强密度收敛，这是高阶方法难以企及的。

第三，在算法几何上，本文的 Bregman 几何是本文最根本的创新。文献 [23,25] 使用 KL 散度作用于有限维概率单形，其几何对象是离散分布，与连续密度 $m \in L^{\bar\beta}$ 的 Banach 结构无关；[21,22] 的 primal-dual 梯度法基于有限差分离散化后的欧氏几何，丧失了连续空间的尺度不变性。而本文的 Bregman 散度 $D_\Psi$ 显式依赖于 $\bar\beta$ 和 $\bar\gamma$，其梯度映射 $\nabla \Psi : X \to X^*$ 将 $L^{\bar\beta}$ 中的密度变量映射至 $L^{\bar\beta'}$，将 $W^{1,\bar\gamma}$ 中的价值函数映射至 $W^{-1,\bar\gamma'}$，完美适配变分不等式 (3.1) 中各项的对偶配对。这种几何适配性是获得定理 1.1 中强收敛的必要条件。

最后，在收敛性保证上，本文结果是首个针对定常 MFG 的、在自然 Banach 空间中建立的、具有明确步长条件的强收敛定理。此前所有工作或仅给出弱收敛（[10] 的存在性理论），或仅给出离散格式的稳定性（[11] 的预处理外梯度法），或仅在有限维近似中证明收敛（[23]）。本文的代价是假设 2.4 的定量单调性 $(g(m)-g(\theta))(m-\theta) \ge c_g |m-\theta|^{\bar\beta}$，这排除了次线性耦合（如 $g(m) = \log(1+m)$），但涵盖了所有物理上重要的幂律耦合（如 $g(m) = m^\beta$）。从提供内容中无法确认本文方法是否可推广至非分离型 Hamiltonian $H(x,p,m)$ 或非幂律增长情形。

3. 通过穿孔域上的双尺度收敛方法对可压缩Navier-Stokes方程进行均匀化

Homogenization of the compressible Navier-Stokes equations via two-scale convergence in perforated domains

作者: Markus Gahn, Kuntal Bhandari, Šárka Nečasová, Maria Neuss-Radu

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.19872v1

摘要翻译: 我们研究在周期性穿孔域中可压缩等熵Navier-Stokes方程的均匀化问题，其中障碍物的尺寸与相邻障碍物之间的距离处于同一量级。利用可通过展开算子刻画的双尺度收敛方法，我们导出了由达西定律（Darcy’s law）确定的相应宏观模型；特别地，宏观密度满足多孔介质方程（porous medium equation）。主要难点在于识别极限中的压力项。我们通过建立密度的强双尺度收敛性来克服该困难，而该强双尺度收敛性是通过控制展开后密度的振荡缺陷测度（oscillation defect measure）实现的。本工作的关键贡献在于发展了一套适用于更复杂可压缩流体模型的方法论框架。此外，针对保守力情形，我们将文献中已有结果推广至绝热常数 $γ> \frac95$。

Comment: 无

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

1. 研究对象和背景

本文研究三维周期穿孔区域 $\Omega_\varepsilon \subset \Omega = (a,b) \subset \mathbb{R}^3$ 上的可压缩等熵Navier-Stokes方程组的均匀化极限行为，其中穿孔（障碍物）的直径与相邻穿孔中心间距同阶，即均属 $O(\varepsilon)$ 量级。该设定对应于Allaire分类中的临界穿孔尺度 $\alpha = 1$（在三维中，若穿孔直径为 $O(\varepsilon^\alpha)$，则 $\alpha = 3$ 对应Darcy尺度，$\alpha > 3$ 为无影响尺度，而本文 $\alpha = 1 < 3$，但因穿孔体积分数不趋于零，故仍属强影响、非平凡均匀化情形）。具体微观模型为：

$$ \begin{cases} \partial_t \varrho_\varepsilon + \mathrm{div}(\varrho_\varepsilon \mathbf{u}_\varepsilon) = 0 & \text{in } (0,T) \times \Omega_\varepsilon, \\ \varepsilon^\lambda \big( \partial_t(\varrho_\varepsilon \mathbf{u}_\varepsilon) + \mathrm{div}(\varrho_\varepsilon \mathbf{u}_\varepsilon \otimes \mathbf{u}_\varepsilon) \big) - \varepsilon^2 \, \mathrm{div}\big( S(\nabla \mathbf{u}_\varepsilon) \big) + \nabla p(\varrho_\varepsilon) = \varrho_\varepsilon \nabla F & \text{in } (0,T) \times \Omega_\varepsilon, \\ \mathbf{u}_\varepsilon = 0 & \text{on } (0,T) \times \partial \Omega_\varepsilon, \end{cases} $$

其中 $p(s) = a s^\gamma$，$a>0$，绝热指数 $\gamma > \frac{9}{5}$；粘性应力张量 $S(A) = \mu \big( A + A^\top - \frac{2}{3} \mathrm{tr}(A) I_3 \big) + \eta \, \mathrm{tr}(A) I_3$，$\mu > 0$，$\eta \ge 0$；外力由势函数 $F \in W^{1,\infty}((0,T) \times \Omega)$ 生成；参数 $\lambda > 1 + \frac{3}{\gamma}$ 是关键缩放因子，确保动量方程中惯性项、粘性项与压力项在均匀化过程中具有可比的渐近阶。初始数据满足 $\varrho_{\varepsilon,0} \ge 0$，$\int_{\Omega_\varepsilon} \varrho_{\varepsilon,0}^\gamma \, dx \le C$，且 $\int_{\Omega_\varepsilon} \varepsilon^\lambda \frac{|\mathbf{m}_{\varepsilon,0}|^2}{\varrho_{\varepsilon,0}} \, dx \le C$（当 $\varrho_{\varepsilon,0}=0$ 时要求 $\mathbf{m}_{\varepsilon,0}=0$）。穿孔结构定义为：设 $Y = (0,1)^3$，$Y_s \Subset Y$ 为单胞内刚性障碍物，$Y_f = Y \setminus Y_s$ 为流体区域，则 $\Omega_\varepsilon = \mathrm{int}\big( \bigcup_{k \in K_\varepsilon} \varepsilon(Y_f + k) \big)$，其中 $K_\varepsilon = \{ k \in \mathbb{Z}^3 : \varepsilon(Y + k) \subset \Omega \}$。

该问题的重要性在于它刻画了高孔隙率、大障碍物情形下可压缩流体的有效宏观输运规律——此时经典Darcy定律不再直接适用，而需耦合密度演化。Masmoudi（2001）首次在此尺度下证明了宏观密度 $\varrho_0$ 满足多孔介质型方程 $|Y_f| \partial_t \varrho_0 + \mathrm{div}_x \big( K \varrho_0^2 \nabla F - a \varrho_0 \nabla_x \varrho_0^\gamma \big) = 0$，其中渗透率张量 $K$ 由Stokes型单元问题确定。后续Feireisl–Lu（2016）、Feireisl–Novotný–Petržela（2020）等将结果推广至Navier-Stokes-Fourier系统，但均采用基于Tartar限制算子与能量法的传统框架，其技术瓶颈在于：如何从弱紧性 $\varrho_\varepsilon \rightharpoonup \varrho_0$ 推出压力项 $p(\varrho_\varepsilon) = a \varrho_\varepsilon^\gamma$ 的强收敛性，从而完成非线性项的极限识别。现有工作依赖于对扩展压力 $\tilde{p}(\varrho_\varepsilon)$（定义在全区域 $\Omega$ 上）的精细估计，但未建立与双尺度分析自然兼容的强紧性机制。本文试图填补的关键空白是：构建一个严格基于双尺度收敛（或等价的展开算子）的统一框架，以系统处理可压缩流体在临界穿孔尺度下的均匀化，并克服压力项极限识别这一核心障碍。特别地，作者旨在将已知结果中对 $\gamma$ 的要求从 $\gamma > \frac{3}{2}$（见Feireisl et al. [18]）或 $\gamma > 2$（见[19]）推进至 $\gamma > \frac{9}{5}$，同时为更复杂的耦合模型（如含化学反应、相变或多物理场）提供可扩展的方法论基础。

2. 主要定理或主要结果

定理2.3给出了全文的核心结论：在 $\gamma > \frac{9}{5}$ 且 $\lambda > 1 + \frac{3}{\gamma}$ 的假设下，微观解序列 $(\varrho_\varepsilon, \mathbf{u}_\varepsilon)$ 在双尺度意义下收敛于极限 $(\varrho_0, \mathbf{u}_0)$，其中 $\varrho_0 \in L^\infty((0,T), L^\gamma(\Omega \times Y_f)) \cap L^{\gamma+\nu}((0,T) \times \Omega \times Y_f)$，$\nu = \frac{2}{3\gamma - 1}$；$\mathbf{u}_0 \in L^2((0,T) \times \Omega; H^1(Y))^3$ 且在 $Y_s$ 上为零；并存在 $p_1 \in L^2((0,T), L^q(\Omega; L^2(Y_f)))$，$q = \min\{2, \gamma\}$，使得 $(\mathbf{u}_0, p_1)$ 满足如下单元问题（在集合 $\{\varrho_0 > 0\} \times Y_f$ 上）：

$$ -2\mu \, \mathrm{div}_y \big( D_y(\mathbf{u}_0) \big) + \nabla_y p_1 = \varrho_0 \nabla F - \nabla_x p_0, \quad \mathrm{div}_y \mathbf{u}_0 = 0, \quad \mathbf{u}_0 = 0 \text{ on } \{\varrho_0 > 0\} \times \Gamma, $$

其中 $p_0 = a \varrho_0^\gamma$，且 $(\mathbf{u}_0, p_1)$ 关于 $Y$ 周期。宏观密度 $\varrho_0 = \varrho_0(t,x)$ 满足多孔介质型方程：

$$ |Y_f| \, \partial_t \varrho_0 + \mathrm{div}_x \big( K \varrho_0^2 \nabla F - a \varrho_0 \nabla_x \varrho_0^\gamma \big) = 0 \quad \text{in } (0,T) \times \Omega, $$

初值条件为 $\varrho_0(0) = |Y_f| \, \varrho_{0,0}$。渗透率张量 $K \in \mathbb{R}^{3 \times 3}$ 由单元问题 $-2\mu \, \mathrm{div}_y(D_y \mathbf{w}_i) + \nabla_y \pi_i = e_i$（$e_i$ 为标准基向量）定义为 $K_{ij} = \int_{Y_f} \mathbf{w}_i \cdot e_j \, dy$。此外，Darcy平均速度 $\bar{\mathbf{u}}_0(t,x) = \int_{Y_f} \mathbf{u}_0(t,x,y) \, dy$ 显式表示为 $\bar{\mathbf{u}}_0 = K (\varrho_0 \nabla F - \nabla_x p_0)$。该定理的证明依赖于两个关键先验估计：Proposition 3.1 给出一致能量界 $\varepsilon^\lambda \|\varrho_\varepsilon |\mathbf{u}_\varepsilon|^2\|_{L^\infty(L^1)} + \|\mathbf{u}_\varepsilon\|_{L^2} + \varepsilon \|\nabla \mathbf{u}_\varepsilon\|_{L^2} + \|\varrho_\varepsilon\|_{L^\infty(L^\gamma)} \le C$；Proposition 3.3 给出改进的压力分解：存在 $q_\varepsilon \in L^2((0,T), H^1(\Omega_\varepsilon)) + L^\infty((0,T), W^{1,\gamma}(\Omega_\varepsilon)) + \varepsilon \big( L^2((0,T) \times \Omega_\varepsilon) + L^\infty((0,T), L^\gamma(\Omega_\varepsilon)) \big)$，使得对任意均值为零的测试函数 $\psi_\varepsilon$，有

$$ \left| \int_0^T \int_{\Omega_\varepsilon} (p(\varrho_\varepsilon) - q_\varepsilon) \psi_\varepsilon \, dx dt \right| \le C \Big( \varepsilon^{\lambda/2} \|B_\varepsilon(\partial_t \psi_\varepsilon)\|_{L^1(L^{2\gamma/(\gamma-1)})} + \varepsilon^{\lambda/s - 2\theta/(s-1)} \varepsilon \|\nabla B_\varepsilon(\psi_\varepsilon)\|_{L^\infty(L^{s'\gamma p/((p-2)\gamma-p)})} + \varepsilon^{\lambda/2} \|B_\varepsilon(\psi_\varepsilon)\|_{W^{1,1}(L^{2\gamma/(\gamma-1)})} \Big), $$

其中 $B_\varepsilon$ 为穿孔区域上的Bogovskii算子，$\theta = \frac{3p-6}{2p}$，$p \in \big( \frac{2\gamma}{\gamma-1}, \min\{6, \frac{4\gamma}{(\gamma-1)(2-s)}\} \big)$，$s \in (1,2]$。该估计的关键在于其右端所有项均以正幂次 $\varepsilon^\kappa$（$\kappa > 0$）衰减，从而保证 $p(\varrho_\varepsilon) - q_\varepsilon \to 0$ 在双尺度意义下强收敛，最终导出 $\varrho_\varepsilon \to \varrho_0$ 的强双尺度收敛性（引理4.5），这是识别非线性压力极限 $p(\varrho_\varepsilon) \to p_0 = a \varrho_0^\gamma$ 的基石。

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

本文的核心技术障碍在于：在临界穿孔尺度下，微观密度 $\varrho_\varepsilon$ 仅具备弱紧性（$L^\infty(L^\gamma)$ 有界），而压力 $p(\varrho_\varepsilon) = a \varrho_\varepsilon^\gamma$ 是非线性函数，其极限识别要求 $\varrho_\varepsilon$ 具备某种形式的强紧性。传统方法（如[19,31]）通过对压力作全局延拓并利用Tartar振荡引理，但该路径与双尺度框架不天然兼容，且难以推广。本文的真正创新点在于将“振荡缺陷测度”（oscillation defect measure）理论与“展开算子”（unfolding operator）深度耦合，构建了一套可迁移的强紧性分析范式。具体而言，作者首先将密度 $\varrho_\varepsilon$ 通过展开算子 $\mathcal{T}_\varepsilon$ 映射为定义在固定域 $\Omega \times Y$ 上的函数 $\mathcal{T}_\varepsilon(\varrho_\varepsilon)(t,x,y)$，其关于 $x$ 是分段常数。随后，作者定义并精确控制该展开密度的振荡缺陷测度 $\mathcal{D}_\varepsilon$，即

$$ \mathcal{D}_\varepsilon := \int_{\Omega \times Y} \big| \mathcal{T}_\varepsilon(\varrho_\varepsilon) - \overline{\mathcal{T}_\varepsilon(\varrho_\varepsilon)} \big|^2 \, dx dy, $$

其中 $\overline{\cdot}$ 表示在 $Y$ 上的平均。控制 $\mathcal{D}_\varepsilon \to 0$ 即等价于 $\mathcal{T}_\varepsilon(\varrho_\varepsilon)$ 在 $L^2$ 中强收敛，进而推出 $\varrho_\varepsilon$ 的强双尺度收敛。实现这一控制的关键工具是Proposition 3.3中建立的改进压力估计。该估计的本质是将压力 $p(\varrho_\varepsilon)$ 分解为 $q_\varepsilon + r_\varepsilon$，其中 $q_\varepsilon$ 属于一个具有良好正则性的函数空间之和（$H^1$、$W^{1,\gamma}$ 及其小量扰动），而余项 $r_\varepsilon$ 在对偶空间中被证明是 $o(1)$ 的。这一分解的构造高度依赖于穿孔区域上Bogovskii算子 $B_\varepsilon$ 的一致估计（见附录A.2），作者利用 $B_\varepsilon$ 将测试函数 $\psi_\varepsilon$ “提升”至具有零散度的向量场 $B_\varepsilon(\psi_\varepsilon)$，从而将压力项的积分转化为对流项、粘性项和外力项的积分。通过对这些项分别应用Gagliardo-Nirenberg插值不等式（引理A.1）、Poincaré不等式及Bogovskii算子的 $L^p$ 有界性（引理A.2），作者得到各项的最优衰减速率，特别是对 $I_2^\varepsilon$ 项的估计中引入了参数 $s \in (1,2]$ 进行插值，巧妙避开了文献[24]中因取 $s=1$ 而导致的时间导数项无法控制的困境。最终，该估计保证了 $r_\varepsilon \to 0$ 在 $L^1$ 意义下，从而 $q_\varepsilon$ 与 $p(\varrho_\varepsilon)$ 具有相同的双尺度极限，而 $q_\varepsilon$ 的强收敛性由其所在函数空间的紧嵌入定理（如 $H^1 \hookrightarrow L^2$）保证。这一整套技术——“展开+振荡缺陷测度+改进压力分解+Bogovskii算子插值估计”——构成了本文的方法论核心，其普适性体现在它不依赖于特定的边界条件或外力形式，为处理更复杂的可压缩流体模型（如含热传导、化学反应或随机穿孔）提供了坚实的基础。

4. 与之前工作的比较

本文与既有工作的比较可从四个维度展开。第一，与Allaire的经典工作[3,4]相比：Allaire针对不可压缩Stokes/Navier-Stokes方程，在 $\alpha < 3$ 的小穿孔尺度下建立了Darcy律，其核心是线性椭圆型单元问题的解的紧性；而本文处理的是可压缩、非线性、临界尺度 $\alpha = 1$ 的问题，宏观模型不再是线性Darcy律，而是由多孔介质方程驱动的非线性输运，且必须解决密度强收敛这一Allaire框架中不存在的根本性困难。第二，与Masmoudi[31]及Feireisl et al.[19]相比：这些开创性工作虽在同一临界尺度下获得了多孔介质方程，但其技术路线基于Tartar限制算子与能量法，其压力估计依赖于对延拓压力 $\tilde{p}(\varrho_\varepsilon)$ 的 $L^2$ 控制，且对 $\gamma$ 的要求为 $\gamma > \frac{3}{2}$（[19]）或隐含更强条件。本文将 $\gamma$ 的下界推进至 $\frac{9}{5}$，这得益于Proposition 3.3中更精细的插值估计：通过选择 $p$ 接近 $\frac{2\gamma}{\gamma-1}$ 和 $s$ 接近1，使得衰减指数 $\frac{\lambda}{s} - \frac{2\theta}{s-1}$ 为正的条件 $\lambda > 1 + \frac{3}{\gamma}$ 恰好允许 $\gamma > \frac{9}{5}$（因为当 $\gamma = \frac{9}{5}$ 时，$1 + \frac{3}{\gamma} = \frac{14}{9} \approx 1.555$，而 $\lambda$ 可取略大于此值）。第三，与近期定量均匀化工作Höfer et al.[24]相比：[24]亦使用双尺度方法，但其压力估计（Lemma 2.4）未显式分离出 $H^1$ 和 $W^{1,\gamma}$ 分量，且对 $I_2^\varepsilon$ 的估计取 $s=1$，导致时间导数项需用 $L^\infty$ 范数控制，这阻碍了对密度时间导数的强紧性分析；本文通过 $s \in (1,2]$ 插值，将时间导数项纳入 $W^{1,1}$ 范数，从而为后续应用Aubin-Lions引理证明 $\varrho_\varepsilon$ 在 $C([0,T]; L^\gamma_w)$ 中的强收敛铺平道路。第四，与一般可压缩均匀化文献（如Feireisl–Lu[15], Lu–Schwarzacher[28]）相比：这些工作均在“小穿孔”尺度（$\alpha > 3$）下进行，此时穿孔影响可忽略，宏观方程与原始方程形式相同；而本文处理的是“大穿孔”尺度，宏观模型发生质变，且证明难度剧增。本文的代价在于：为获得 $\gamma > \frac{9}{5}$ 的结果，必须将外力限制为保守力 $\nabla F$（Remark 2.1），若考虑一般外力 $\mathbf{f} \in L^\infty$，则需加强假设至 $\gamma \ge 2$（Remark 3.2），这反映了非线性与外力复杂性之间的固有张力。综上，本文并非简单复现已有结论，而是通过一套原创的、基于展开算子与振荡缺陷测度的强紧性框架，在更广的绝热指数范围、更自然的双尺度语境下，为可压缩流体均匀化理论奠定了新的方法论基石。

4. 初始速度属于 $VMO^{-1}$ 的三维非均匀不可压缩 Navier-Stokes 方程组的解

Solutions of the 3D inhomogeneous incompressible Navier-Stokes system with initial velocity in $VMO^{-1}$

作者: Ruilin Hu, Quoc-Hung Nguyen, Feng Shao, Dongyi Wei, Ping Zhang, Zhifei Zhang

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.20207v1

摘要翻译: 本文中，我们建立了三维非均匀不可压缩 Navier-Stokes 方程组在初值 $(ρ_0,u_0) \in C^1 \times (L^2 \cap VMO^{-1})$ 下强解的局部存在性，其中 $ρ_0$ 具有正的下界。进一步地，若 $ρ_0 \in C^2$ 且 $||ρ_0-1||_{L^\infty}+||u_0||_{BMO^{-1}}$ 充分小，则我们证明了解的整体存在性。为实现此目标，我们利用输运方程的估计来获得密度的正则性，并对动量方程采用一种新的冻结系数法。

Comment: 无

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

1. 研究对象和背景

本文研究三维非均匀不可压缩 Navier-Stokes 方程组的初值问题：

$$ \begin{cases} \partial_t \rho + u \cdot \nabla \rho = 0, \\ \partial_t (\rho u) + \mathrm{div}(\rho u \otimes u) - \Delta u + \nabla P = 0, \\ \nabla \cdot u = 0, \\ (\rho, u)|_{t=0} = (\rho_0, u_0), \end{cases} \quad (t,x) \in \mathbb{R}_+ \times \mathbb{R}^3, $$

其中 $\rho(t,x) > 0$ 表示密度，$u(t,x) \in \mathbb{R}^3$ 为速度场，$P(t,x)$ 为压力。该系统描述由若干不可混溶、密度可变的不可压缩流体组成的混合物，其物理基础见文献 [31]。与经典 Navier-Stokes（NS）方程不同，该系统中密度满足输运方程，导致动量方程系数 $\rho$ 依赖于时空变量，破坏了线性结构和标准能量守恒律；同时，$\rho$ 的非恒定性使 Leray 投影 $\mathbb{P} = \mathrm{Id} - \nabla \Delta^{-1} \nabla \cdot$ 不再与时间导数对易，且 $\mathbb{P}$ 在 $L^\infty$ 上无界，显著增加分析难度。

该系统的尺度不变性由变换

$$ (\rho(t,x), u(t,x), P(t,x)) \mapsto \left( \rho(\lambda^2 t, \lambda x),\; \lambda u(\lambda^2 t, \lambda x),\; \lambda^2 P(\lambda^2 t, \lambda x) \right) $$

给出。因此，临界函数空间需满足范数在该缩放下不变。在 $\mathbb{R}^3$ 中，典型临界空间包括齐次 Sobolev 空间 $\dot{H}^{1/2}$、Lebesgue 空间 $L^3$、Besov 空间 $\dot{B}^{-1+3/p}_{p,q}$，以及本文核心关注的 $\mathrm{BMO}^{-1}$ 与 $\mathrm{VMO}^{-1}$：前者定义为 $\|u\|_{\mathrm{BMO}^{-1}} = \sup_{t>0} t^{1/2} \|e^{t\Delta} u\|_{L^\infty} + \sup_{x,R>0} R^{-3/2} \left( \int_0^{R^2} \int_{B(x,R)} |e^{t\Delta} u(y)|^2 \, dy dt \right)^{1/2}$，后者为其闭子空间 $\mathrm{VMO}^{-1} = \{ f \in \mathcal{S}' : \lim_{b \to 0} \|f\|_{\mathrm{BMO}^{-1}_b} = 0 \}$，其中局部化范数 $\|f\|_{\mathrm{BMO}^{-1}_b}$ 将上式中 $t$ 和 $R$ 的上界限制为 $b$。

已有理论可分为三类：（1）高正则初值：Kazhikhov（1974）[28] 在 $\rho_0 \in L^\infty$ 且 $\inf \rho_0 > 0$、$u_0 \in H^1$ 下建立全局强解；（2）临界初值但密度近常：Danchin [16] 在 $\rho_0 - 1 \in \dot{B}^{3/2}_{2,1}$、$u_0 \in \dot{B}^{-1/2}_{2,1}$ 小初值下得全局解；Zhang [35] 将速度初值推广至 $\dot{B}^{1/2}_{2,1}$，[20] 进一步延拓至 $\dot{B}^{-1+3/p}_{p,1}$；（3）低正则密度：Hao–Li–Zhang [24,25] 在 $\rho_0 \in L^\infty$（仅需 $\inf \rho_0 > 0$）、$u_0 \in \dot{H}^{1/2}$ 小初值下获全局解。然而，所有前述结果均未触及 $u_0 \in \mathrm{VMO}^{-1} \cap L^2$ 这一更弱的临界初值设定——该空间严格嵌入 $\dot{B}^{-1}_{\infty,\infty}$ 但不包含于任何 $L^p$（$p<\infty$），且其元素可能具有分形奇异性（如 $|x|^{-1}$ 型奇点），使得经典能量方法与 Besov 插值技术失效。本文旨在填补这一关键空白：首次在 $u_0 \in \mathrm{VMO}^{-1} \cap L^2$、$\rho_0 \in C^1$ 且 $\inf \rho_0 > 0$ 的纯临界初值框架下，建立局部强解的存在性；并进一步在 $\rho_0 \in C^2$、$\|\rho_0 - 1\|_{L^\infty} + \|u_0\|_{\mathrm{BMO}^{-1}} \ll 1$ 下证明全局存在性。此推进不仅将速度初值的临界性推向极限，更揭示了非均匀性与临界正则性之间深刻的兼容机制。

2. 主要定理或主要结果

定理 1.1（局部存在性）：设 $\rho_0 \in C^1$ 满足 $0 < C_1 < \rho_0 < C_2 < \infty$，且 $u_0 \in L^2 \cap \mathrm{VMO}^{-1}$。则存在 $b \in (0,1)$ 与 $\varepsilon_0 > 0$，使得若 $\|u_0\|_{\mathrm{BMO}^{-1}_b} \leq \varepsilon_0$，则系统 (1.1) 在时间区间 $[0,T]$ 上存在强解 $(\rho,u,P)$，其中 $T \sim b^2$（即存在普适常数 $C'_1,C'_2$ 使得 $C'_1 b^2 \leq T \leq C'_2 b^2$），且解满足

$$ \|u\|_{X_T} \leq C \varepsilon_0^{5/6} \|u_0\|_{L^2}^{1/6} \|\rho_0\|_{C^1}, $$

其中混合范数 $\|u\|_{X_T} = \|u\|_{X^1_T} + \|u\|_{X^2_T}$ 定义为

$$ \|u\|_{X^1_T} = \sup_{t\leq T} t^{1/2} \|u(t)\|_{L^\infty} + \|u\|_{V^b_T} + \sup_{t\leq T} t^{1+\kappa/2} \|\nabla u(t)\|_{\dot{C}^\kappa} + \sup_{0<s<t\leq T} \frac{s^{1/2+\kappa_1} \|u(t)-u(s)\|_{L^\infty}}{(t-s)^{\kappa_1}}, $$

$$ \|u\|_{X^2_T} = \sup_{t\leq T} \left( t^{5/12} \|u(t)\|_{L^{12}} + \sup_{0<s<t\leq T} \frac{s^{5/12+\kappa_2} \|u(t)-u(s)\|_{L^{12}}}{(t-s)^{\kappa_2}} \right), $$

其中 $\kappa = 1/700$, $\kappa_1 = 18/35 > 1/2$, $\kappa_2 = 5/7$。此外，解还满足能量型先验估计

$$ \sup_{t\leq T} \left( \|u(t)\|_{L^2} + t^{1/2} \|\nabla u(t)\|_{L^2} + t \|u_t(t)\|_{L^2} \right) + \left\| \left( \nabla u, t^{1/2} u_t, t \nabla u_t \right) \right\|_{L^2_t L^2_x} \leq C(\varepsilon_0, \|u_0\|_{L^2}). $$

定理 1.2（全局存在性）：设 $M > 0$，$u_0 \in L^2$，$\rho_0 \in C^2$ 满足 $\|u_0\|_{L^2} + \|\rho_0\|_{C^2} \leq M$。则存在仅依赖于 $M$ 的常数 $\varepsilon_2 \in (0,1/2)$，使得若

$$ \|\rho_0 - 1\|_{L^\infty} + \|u_0\|_{\mathrm{BMO}^{-1}} < \varepsilon_2, $$

则系统 (1.1) 存在全局强解 $(\rho,u,P)$ 定义于 $[0,\infty) \times \mathbb{R}^3$。此处 $\rho_0 \in C^2$ 的正则性要求是本质的（见 Remark 1.3），源于局部存在时间 $T$ 对 $\varepsilon_0$ 的隐含依赖性；而小初值条件以 $L^\infty$-密度扰动与 $\mathrm{BMO}^{-1}$-速度范数之和度量，体现了对非均匀性与临界性的双重控制。

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

本文的核心障碍在于：（1）密度 $\rho$ 由输运方程 $\partial_t \rho + u \cdot \nabla \rho = 0$ 控制，其正则性完全依赖于速度场 $u$ 的 Lipschitz 性；但若 $u_0 \in \mathrm{VMO}^{-1}$，则 $\|\nabla u(t)\|_{L^\infty} \sim t^{-1}$，其时间积分发散，导致标准输运理论失效；（2）动量方程中 $\rho$ 作为时空调制系数，破坏了标准热半群的光滑性估计，尤其当 $\rho$ 仅具 $C^1$ 正则性时，无法直接应用 $L^p$-理论或标准冻结系数法。

作者的真正创新点在于构建了一套精细的“双尺度冻结系数—加权 Hölder—分层插值”分析框架。首先，对任意固定 $x_0 \in \mathbb{R}^3$，将动量方程在 $x_0$ 处冻结系数 $\rho_0(x_0)$，得到线性主部 $\rho_0(x_0) \partial_t u - \Delta u$，并将其余项（含 $\rho - \rho_0(x_0)$、$\partial_t \rho$ 等）视为余项 $u_{R,i,x_0}$。由此获得解的逐点分解

$$ u(t,x) = \sum_{\star=L,N,R_1,R_2,R_3} \tilde{u}_\star(t,x), \quad \text{其中 } \tilde{u}_\star(t,x) := u_{\star,x_0}(t,x)|_{x_0=x}. $$

关键突破在于对余项 $u_{R_2,x_0} = \int_0^t \partial_t K_{x_0}(t-\tau) * \mathbb{P}(\varrho u)(\tau) \, d\tau + \mathbb{P}(\varrho u)(t)$ 的处理：其末项 $\mathbb{P}(\varrho u)(t)$ 因 $\varrho = \rho - \rho_0$ 仅具 $C^1$ 正则性而缺乏空间光滑性，但作者引入时间加权 Hölder 范数 $\|u\|_{X^1_T}$ 中的 $s^{1/2+\kappa_1} \|u(t)-u(s)\|_{L^\infty}/(t-s)^{\kappa_1}$ 项，通过充分利用 $\varrho$ 的时间 Hölder 连续性（Proposition 2.2 给出 $\||D|^{\alpha'} \varrho(t) - |D|^{\alpha'} \varrho(s)\|_{L^\infty} \lesssim (t-s)^{(\alpha-\alpha')/2}$）补偿空间损失。其次，对 $u_{R_1,x_0}$ 中的奇异项 $\nabla \Delta^{-1}(\nabla \rho_0 \cdot u)(t)$，作者利用附录 7.1 的加权能量估计 $\sup_{t\leq T} t \|u_t(t)\|_{L^2} < \infty$，再经 Gagliardo–Nirenberg 插值 $\|\nabla \Delta^{-1}(\nabla \rho_0 \cdot u)\|_{\dot{C}^\kappa} \lesssim \|u\|_{L^2}^{1-\theta} \|u_t\|_{L^2}^{\theta}$（$\theta = \kappa/(1+\kappa)$）获得所需控制。最后，对密度 $\rho$ 的正则性，作者设计辅助范数 $\|u\|_T = \sup_{t\leq T} (t^{1/2} \|u(t)\|_{L^\infty} + t \|u(t)\|_{\dot{C}^1})$，并结合流映射 $X(t,x)$ 的 Gronwall 型估计 $|X(s,x) - X(s,y)| \leq \exp(C \|u\|_T \log(t/s)) |X(t,x) - X(t,y)|$，最终导出 $\rho(t) \in C^\gamma$ 一致有界性（Proposition 2.1）。整套工具环环相扣，将 $\mathrm{VMO}^{-1}$ 的振荡特性、输运方程的流形几何、热核的时空调制效应与加权插值的精细平衡融为一体。

4. 与之前工作的比较

本文相较既有工作实现了三重实质性推进：

（1）初值临界性维度：此前所有关于非均匀 NS 的全局存在性结果均要求 $u_0$ 属于比 $\mathrm{BMO}^{-1}$ 更强的空间。Danchin [16] 及后续 [1,18] 依赖 $\dot{B}^{-1/2}_{2,1}$ 或 $\dot{B}^{-1+3/p}_{p,1}$（$p<\infty$），其范数控制高频振荡的能力弱于 $\mathrm{BMO}^{-1}$；Zhang [35] 与 [20] 虽进入 $\dot{B}^{1/2}_{2,1}$，但该空间仍嵌入 $L^3$，而 $\mathrm{VMO}^{-1} \not\subset L^3$。本文首次将速度初值放宽至 $\mathrm{VMO}^{-1} \cap L^2$，这是目前已知最弱的临界初值空间（因 $\dot{B}^{-1}_{\infty,\infty}$ 已被证明不适定 [7]），且 $\mathrm{VMO}^{-1}$ 是 $\mathrm{BMO}^{-1}$ 中唯一能保证热半群 $e^{t\Delta} u_0$ 的 $L^\infty$-模具有 $t^{-1/2}$ 衰减的子空间，故此推进具有本质性。

（2）密度正则性维度：以往工作或要求 $\rho_0$ 接近常数（如 $\dot{B}^{3/2}_{2,1}$ [16]），或要求 $\rho_0 \in L^\infty$ 但需 $u_0$ 具更高正则性（如 $\dot{H}^{1/2}$ [24,25]）。本文定理 1.1 仅假设 $\rho_0 \in C^1$ 且 $\inf \rho_0 > 0$，不施加任何小ness或频域条件，这与 Kazhikhov [28] 的经典 $H^1$-速度框架形成互补，首次在纯临界速度下实现对一般 $C^1$ 密度的处理。

（3）方法论代价与局限：本文的推进伴随明确代价。首先，局部存在性未解决唯一性（Remark 1.2 明确指出：$\rho_0 \in C^1$、$u_0 \in L^2 \cap \mathrm{VMO}^{-1}$ 下的唯一性仍是开放问题），这源于 $\mathrm{BMO}^{-1}$ 中缺乏有效的 $f = f_1 + f_2$ 分解（$f_1$ 属于某亚临界空间），而 Koch–Tataru 方法 [29] 在非均匀情形下失效。其次，全局存在性要求 $\rho_0 \in C^2$（而非 $C^1$），这是为控制余项 $u_{R_3}$ 中的 $\rho_0(y) - \rho_0(x_0)$ 项所必需（见 Remark 1.3 与第 6 节），反映了非均匀性对高阶正则性的额外索取。最后，本文未处理真空情形（即 $\inf \rho_0 = 0$），此方向仍由 Lions [31] 的弱解理论主导，而强解框架下的真空问题仍是重大挑战。

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

1. 修正 Camassa-Holm 方程中的谱稳定性

Spectral stability in the modified Camassa-Holm equation

作者: Lili Fan, Hongjun Gao, Ji Li

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.20273v1

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

摘要翻译

中文标题: 修正 Camassa-Holm 方程中的谱稳定性

摘要翻译: 我们研究具有三次非线性的修正 Camassa-Holm 方程的小振幅、周期行波解的谱稳定性。更准确地说，我们在谱平面中零点的邻域内分析相应线性化算子的 $L^2(\mathbb{R})$-谱。受近期基于 Kato 摄动理论的一种新颖方法启发 [Berti 等, Full description of Benjamin-Feir instability of Stokes waves in deep water, \textit{Invent. Math.}, 230 (2022), 651-711.]，我们给出了该线性化算子——一个具有周期系数的积分-微分算子——在零点附近的谱的完整刻画，从而证明此类波不发生调制不稳定性。此外，谱分析揭示了一个显著的阈值现象：当波数满足 $k^2\leq 3$ 时，此类波具有谱稳定性；而当 $k^2>3$ 时，则出现不稳定性。

2. 具可积振荡系数的椭圆方程的广义Morrey-Campanato估计

Generalized Morrey-Campanato estimates for elliptic equations with coefficients of integrable oscillation

作者: Laurent Seppecher

类别: math.AP, math.FA

发布日期: 2026-06-18

arXiv ID: 2606.20237v1

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

摘要翻译

中文标题: 具可积振荡系数的椭圆方程的广义Morrey-Campanato估计

摘要翻译: 本工作研究散度型椭圆方程 $-\mathrm{div}(a\nabla u) = \mathrm{div}\, F$ 的弱解在系数 $a$ 与源项 $F$ 均仅满足低正则性假设下的正则性性质。我们引入广义Morrey空间与广义Campanato空间，其定义通过将经典定义中的逐点一致有界性要求替换为适当的可积性条件而得到推广。在此框架下，我们建立了弱解梯度在这些广义空间中的正则性估计。作为应用，我们重新得到了经典的Hölder估计与Lebesgue估计，并推导出分数阶Sobolev正则性结果。特别地，所提出的办法可在系数可能不连续、且解的梯度未必局部有界的场合下，给出分数阶Sobolev估计。

3. 关于 $(p,q)$-拉普拉斯方程的刘维尔定理

Liouville Theorem for $(p,q)$-Laplace Equations

作者: Yang Zhou, Hua Zhu

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.20033v1

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

摘要翻译

中文标题: 关于 $(p,q)$-拉普拉斯方程的刘维尔定理

摘要翻译: 我们运用向量场方法，建立了欧氏空间 $\mathbb{R}^n$ 中一类 $(p,q)$-拉普拉斯方程的刘维尔型定理。通过调整微分恒等式中的指数，我们在次临界范围 $p-1<α<q^*-1$ 内证明了解的不存在性，其中 $q^*=nq/(n-q)$。该方法依赖于构造适当的微分恒等式、利用截断函数进行精确的积分估计，并结合对截断误差的符号控制与衰减分析。

4. 离散Riesz变换的最优维数依赖的 $\ell^p$ 和 $\ell^{1,\infty}$ 估计

Optimal dimension-dependent $\ell^p$ and $\ell^{1,\infty}$ estimates of the discrete Riesz Transforms

作者: Junjie Shao, Hanli Tang, Zewei Xu

类别: math.CA, math.AP

发布日期: 2026-06-18

arXiv ID: 2606.19841v1

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

摘要翻译

中文标题: 离散Riesz变换的最优维数依赖的 $\ell^p$ 和 $\ell^{1,\infty}$ 估计

摘要翻译: 本文研究定义在 $\mathbb{Z}^d$ 上的离散Riesz变换 $R_{\text{dis}}^{(k)}$ 的最优维数依赖的 $\ell^p$ 范数，其由奇异卷积核 $K_k(m)=c_d m_k/|m|^{d+1}$ 给出，其中 $c_d=Γ(\frac{d+1}{2})/π^{(d+1)/2}$ 。我们证明：对固定的 $1<p<\infty$，当 $d\to \infty$ 时，

$$ \|R_{dis}^{\left( k \right)}\|_{\ell ^p\left( \mathbb{Z}^d \right) \rightarrow \ell ^p\left( \mathbb{Z}^d \right)}=2c_d\left( 1+\frac{\left( \sqrt{2}+o\left( 1 \right) \right) d}{2^{\frac{d}{2}}} \right) . $$

由于 Stirling 公式给出 $c_d\sim(\frac{d-1}{2eπ})^{\frac{d-1}{2}}\sqrt{\frac{d-1}π}$，故 $R_{\text{dis}}^{(k)}$ 的算子范数随 $d\to\infty$ 超指数增长，从而否定了 Bañuelos、Kim 与 Kwaśnicki 在 \cite{BKK} 中提出的猜想。此外，我们还建立了 $R_{\text{dis}}^{(k)}$ 的最优维数依赖的 $\ell^{1,\infty}$ 估计。

5. 四维 Zakharov 系统在基态下方的散射

Scattering for the 4D Zakharov system below the ground state

作者: Timothy Candy, Kenji Nakanishi

类别: math.AP

发布日期: 2026-06-17

arXiv ID: 2606.19634v1

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

摘要翻译

中文标题: 四维 Zakharov 系统在基态下方的散射

摘要翻译: 对于四维空间中的 Zakharov 系统，我们证明：所有位于基态势阱内的解在能量空间中整体存在且散射，无需附加对称性等其他限制条件。该证明已被 [3] 归结为排除存在一条轨道上预紧的极小非散射解。本文完成证明的最后一步，即通过结合两种依赖于轨道运动的不同论证，排除势阱内存在预紧解的可能性。

不相关论文（基本信息）

1. 约束在半空间上的库仑气体极小化子

Minimizers for Coulomb gases constrained to a halfspace

作者: Rupert L. Frank, Paata Ivanishvili, Clara Torres-Latorre

类别: math.AP, math-ph

发布日期: 2026-06-18

arXiv ID: 2606.20484v1

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

摘要: We consider a family of optimization problems, based on a mean-field description of particles interacting through Coulomb forces in a quadratic trap. In addition, the particles are constrained to lie in a halfspace and we are interested in the way the particle distribution changes as the halfspace varies. In particular, we can prove the existence of a phase transition, thereby settling a recent conjecture by Byun, Forrester, Majumdar and Schehr.

2. 含 $(1,q)$-拉普拉斯算子与不连续超线性非线性的椭圆问题解的存在性

Existence of solutions for elliptic problems involving the $(1,q)$-Laplacian operator and a discontinuous superlinear nonlinearity

作者: Marcos A. V. Costa, Olímpio H. Miyagaki, Marcos T. O. Pimenta

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.20217v1

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

摘要: In this paper, we study a class of quasilinear elliptic problems involving the $(1,q)-$Laplacian operator and a discontinuous superlinear nonlinearity governed by the Heaviside function. The main difficulty of the problem arises from the presence of the $1$-Laplacian operator, whose natural setting is the Space of Functions of Bounded Variation. Our approach is based on an approximation method involving $(p,q)-$Laplacian problems as $p\to1^+$. As a consequence, we prove the existence of a nontrivial and nonnegative solution belonging to $W^{1,p}_0(Ω)$, in an appropriate weak sense. Moreover, we investigate the asymptotic behavior of the solutions as $β\to0^+$, showing that the family of solutions converges to a solution of the limit problem without discontinuity.

3. 关于一个非齐次混合局部-非局部方程的弱解与粘性解

On weak and viscosity solutions to a nonhomogeneous mixed local-nonlocal equation

作者: R. Lakshmi, Sekhar Ghosh

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.20099v1

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

摘要: This paper explores the relationship between weak and viscosity solutions to a nonhomogeneous mixed local and non-local $p$-Laplace equation in a bounded Lipschitz domain in $\mathbb{R}^N$. Under certain conditions, we derive the comparison principle for weak subsolutions and weak supersolutions to the problem. For $1<p<\infty$, we establish that continuous weak supersolutions to the problem are viscosity supersolutions, using the comparison principle. Furthermore, we show that bounded viscosity supersolutions are weak supersolutions for $p \geq 2$.

4. 带状区域中过定毛细问题的分歧

Bifurcation of overdetermined capillary problems in a strip domain

作者: Yuanyuan Lian, Pieralberto Sicbaldi

类别: math.AP

发布日期: 2026-06-18

arXiv ID: 2606.19885v1

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

摘要: In this paper, we consider the classical overdetermined capillary problem: \begin{equation*} \begin{cases} \mathrm{div} \left(\frac{\nabla u}{\sqrt{1+|\nabla u|^2}}\right) - bu =0 &\mbox{in} Ω, \partial_ν u=κ&\mbox{on}\partialΩ, u=c &\mbox{on}\partialΩ, \end{cases} \end{equation*} where $b$, $c$ and $κ$ are positive constants, and $Ω\subset \mathbb{R}^2$. When $Ω$ is an infinite strip, i.e., a domain bounded by two parallel straight lines, there exists a unique one-dimensional solution (called the trivial solution) to this problem. By means of a bifurcation argument, we establish the existence of a critical period $T_*$ at which a branch of non-trivial solutions bifurcates from the trivial one. These solutions are genuinely two-dimensional and are defined in unbounded periodic domains $Ω$ that are diffeomorphic to an infinite strip, yet whose boundaries are no longer straight lines. This result offers a significant physical interpretation in the context of capillary phenomena.

5. 复Hessian流的Cauchy-Dirichlet问题：从先验估计到多重位势论

The Cauchy-Dirichlet Problem for Complex Hessian Flows: From A Priori Estimates to Pluripotential Theory

作者: Haoyuan Sun

类别: math.AP, math.CV

发布日期: 2026-06-18

arXiv ID: 2606.19713v1

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

摘要: We study the Cauchy–Dirichlet problem for parabolic complex Hessian equations on Hermitian manifolds and on bounded strictly m-pseudoconvex domains. In the smooth setting, we prove global existence and uniqueness of classical solutions under the presence of an admissible parabolic subsolution, by establishing a priori estimates up to the parabolic boundary. The estimates combine parabolic boundary techniques for complex Hessian equations with interior second order estimates and a blow-up argument. We then develop a general pluripotential framework for degenerate right-hand sides with L^p densities, p>n/m, and bounded Cauchy–Dirichlet data. Since the usual automorphism and Walsh-type arguments do not directly apply in a variable Hermitian background, we use approximation by smooth data, balayage, parabolic Perron envelopes, and a continuous obstacle approximation based on Harvey–Lawson–Plis subequation theory. The resulting solution is continuous for positive time, locally uniformly Lipschitz and semi-concave in time, and continuous up to the initial slice when the initial datum is continuous. We also prove a parabolic comparison principle via time regularization, Riemann sum approximations, and mixed Hessian inequalities.

6. 一种基于容度的加权图上带位势半线性椭圆不等式的方法

A Capacitary Approach to Semilinear Elliptic Inequalities with Potentials on Weighted Graphs

作者: Mohamed Jleli, Bessem Samet

类别: math.AP

发布日期: 2026-06-17

arXiv ID: 2606.19650v1

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

摘要: We develop a capacitary approach to semilinear elliptic inequalities on weighted graphs with a potential. More precisely, we study the nonexistence of nontrivial nonnegative solutions of

$$ Δu+w(x)u+v(x)u^σ\le0 \qquad\text{in }V, $$

where $(V,ω,μ)$ is a connected, locally finite weighted graph, $Δ$ is the associated graph Laplacian, $σ>1$, $v>0$, and $w$ is a real-valued potential. The potential term is handled by means of a positive solution $H$ of $ΔH+wH=0$, which transforms the operator $Δ+w$ into the $H$-Laplacian associated with a new weighted graph. Our main nonexistence criterion is formulated directly in terms of cut-off functions and the regions where their $H$-Laplacian is controlled. Unlike metric criteria based on pseudo-metric annuli, our formulation determines the capacitary sets from the support of the $H$-Laplacian estimates for the cut-off functions. We provide an example showing that our result applies in situations not covered by previous nonexistence criteria based on structural lower bounds or pseudo-metric annular volume estimates. We also show that the growth exponent in our capacitary condition is sharp by constructing an example for which the condition fails by an arbitrary power $R^\varepsilon$, while a positive nontrivial solution exists.

7. 带维数畸变的有号分形测度的最优传输：一个变分刻画

Optimal transport of signed fractal measures with dimensional distortion: a variational characterization

作者: Bwo’nyahre Baidi Barthelemy, Kouakep Tchaptchie Yannick, Houpa Danga Duplex Elvis

类别: math.AP

发布日期: 2026-06-17

arXiv ID: 2606.19631v1

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

摘要: We extend the optimal transport theory for signed measures supported on Ahlfors-regular fractal sets (Bwo’Nyahre et al., 2026) to allow a controlled dimensional distortion between source and target. A penalization term $\varepsilon Φ(d_s(x) - d_t(y))$ – where $Φ$ is a fixed smooth strictly convex function and $d_s, d_t$ are the local Hausdorff dimensions of the fractal supports – is added to the transport cost on inter-sign regions, with~$\varepsilon \ge 0$ controlling the tolerance for distortion. Under hypotheses H1–H7, we prove: the existence and uniqueness of an optimal transport map~$T^{\varepsilon}$ for every~$\varepsilon > 0$; coupled Monge–Ampère equations with a distortion correction term, generalizing the classical Brenier–Caffarelli equation; a double Legendre–Fenchel characterization of the optimal potentials, giving a complete variational description of the transport in each of the four sign regimes. The double Legendre–Fenchel system (Theorem~4.2) is the central contribution: it shows that the optimal potentials are the unique fixed points of a system of conjugacy equations, one per transport regime, and it provides the foundation for numerical algorithms and asymptotic analysis.

8. 实测度上符号间最优传输中位置惩罚函数的正则性

Regularity of the positional penalization function in inter-sign optimal transport on real measures

作者: Bwo’nyahre Baidi Barthelemy, Kouakep Tchaptchie Yannick, Houpa Danga Duplex Elvis

类别: math.AP

发布日期: 2026-06-17

arXiv ID: 2606.19621v1

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

摘要: We study the Monge–Kantorovich optimal transport problem between two signed measures~$μ$ and~$ν$ on convex compact subsets of~$\mathbb{R}^d$, with a positional penalization function~$λ(x, y)$ that modulates the cost of inter-sign transport. Using four independent positive measures~$(π^{++}, π^{+-}, π^{-+}, π^{--})$ as decision variables, we prove that the admissible set~$\mathcal{A}(μ, ν)$ is weakly-$*$ compact and non-empty if and only if $μ^+(X) = ν^+(Y)$ and~$μ^-(X) = ν^-(Y)$. Strong duality is established via the Kantorovich minimax theorem, yielding a new compatibility condition on~$λ$ at the intersection of inter-sign supports. The penalization~$λ$ is shown to be Lipschitz and to admit Alexandrov second derivatives almost everywhere. Modified Monge–Ampère equations governing inter-sign transport maps are derived in the Alexandrov sense, with well-posedness characterized by $σ\det(D^2_{yx}Λ) e > 0$. The classical Brenier equation is recovered in the limit~$λ\to 0$.