Homogenization for the $p$-Laplacian in a $d$-dimensional ball perforated along the unit sphere: the critical case $p=d$

技术摘要：沿单位球面穿孔的 $d$ 维球体中 $p$-Laplacian 的均匀化研究（临界情况 $p=d$）

问题陈述
本文研究了穿孔区域 $B(0, \rho) \setminus \Gamma_\varepsilon \subset \mathbb{R}^d$ 中（其中 $\rho > 1$）$p$-Laplacian 边界值问题的 Perron 解 $u_\varepsilon$ 的渐近行为。该区域由分布在单位球面 $S^{d-1}$ 附近的许多小紧致腔体 $\Gamma_\varepsilon$ 穿孔而成。这些腔体以“锚点” $S(\varepsilon)$ 为中心，这些锚点在球面上呈渐近等分布，分离尺度为 $\varepsilon$。每个腔体是紧致集 $K_s$ 的缩放副本 $s + \alpha(\varepsilon)\varepsilon K_s$，其中缩放因子 $\alpha(\varepsilon) \to 0$。

边界条件为：在 $\Gamma_\varepsilon$ 上 $u_\varepsilon = 1$，在 $\partial B(0, \rho)$ 上 $u_\varepsilon = 0$。本研究特别关注临界情况，即维度等于指数 $p = d > 1$。这一情况与先前文献中分析过的亚临界（$p < d$）和超临界（$p > d$）情形不同。在这种情况下，腔体与体势之间的相互作用受对数容量机制控制，导致腔体尺寸呈现指数级的临界缩放。

方法论
作者结合了针对 $d$-调和设置（$p=d$）定制的变分方法、容量估计和比较原理。

1. 缩放与临界窗口： 本文通过参数 $\tau = \lim_{\varepsilon \downarrow 0} [\varepsilon \log(1/\alpha(\varepsilon))]^{-1}$ 确定了临界缩放机制。这意味着腔体尺寸的行为满足 $\alpha(\varepsilon) \approx \exp(-1/(\tau\varepsilon))$。参数 $\tau \in [0, \infty]$ 决定了解的宏观行为。
2. 构造假设函数（Ansatz）： 构建了两种近似函数：
   • 一种显式的径向假设函数 $u_\varepsilon^A$，它结合了一个体 $d$-调和势 $U_\rho$（经常数 $A$ 缩放）和一系列局部的径向凸起。
   • 一种形状自适应假设函数 $\hat{u}_\varepsilon^A$，它用单个腔体 $K_s$ 的精确 $d$-容量势 $V^s_R$ 取代了径向凸起。
3. 能量最小化： 分析的核心依赖于最小化 $d$-能量泛函 $E(v) = \int |\nabla v|^d$。作者推导了假设函数的渐近能量，表明总能量收敛到一个取决于参数 $A$ 的泛函 $\phi_\tau(A)$。最优系数 $A^*$ 被确定为该泛函的极小值点。
4. 局部估计与分离：
   • 局部能量界限： 通过“好/坏锚点”分类法和局部振荡估计（引理 3.1–3.3），作者确立了解在单位球面附近的球面壳层上接近一个常数值。
   • 分离定理： 一个关键的技术工具（定理 6.1）是 $p$-Laplacian（$p=d$）的分离定理。它限制了位于腔体与单位球面之间区域内调和函数的振荡，证明了腔体的影响随远离球面而迅速衰减。这依赖于极大值原理和特定的屏障构造。
5. 收敛性证明： 作者证明了真实解 $u_\varepsilon$ 在 $L^\infty$ 范数和 Sobolev 空间 $H^{1,d}_0$ 中都收敛于具有最优系数 $A^*$ 的假设函数。这是通过将 $u_\varepsilon$ 的能量与假设函数的能量进行比较，并利用 Clarkson 不等式来控制梯度差来实现的。

主要结果
在满足锚点等分布（H1）、腔体一致正容量（H2）和临界缩放（H3）的假设下，本文建立了以下主要定理：

1. 均匀化极限（定理 1.1）： 当 $\varepsilon \to 0$ 时，$u_\varepsilon$ 在 $\bar{B}(0, \rho) \setminus S^{d-1}$ 的紧子集上一致收敛于一个缩放的径向势：
   $$u_\varepsilon(x) \to A^* U_\rho(x)$$
   其中 $U_\rho(x) = \min\{1, 1 - \log|x|/\log\rho\}$ 是单位球的 $d$-调和势。缩放因子 $A^*$ 由 $\tau$ 显式确定：
   $$A^* = \begin{cases} \frac{\tau \sigma^{1/(d-1)} \log \rho}{1 + \tau \sigma^{1/(d-1)} \log \rho} & 0 \le \tau < \infty \\ 1 & \tau = \infty \end{cases}$$
   这里 $\sigma$ 与锚点的密度相关。穿孔集的 $d$-容量收敛至 $\omega_{d-1} (\log \rho)^{-(d-1)} (A^*)^{d-1}$。

2. 强收敛与一致逼近（定理 1.2）： 解 $u_\varepsilon$ 在 $L^\infty$ 范数和 $H^{1,d}_0$ 范数下均收敛于形状自适应假设函数 $\hat{u}_\varepsilon^{A^*}$（从而也收敛于显式径向假设函数 $u_\varepsilon^{A^*}$）。这证实了该假设函数不仅捕捉到了宏观极限，还捕捉到了腔体附近的微观轮廓。

意义与主张
本文声称解决了在临界情况 $p=d$ 下，穿孔集中在超曲面上的 $p$-Laplacian 的均匀化问题。

• 临界缩放： 与缩放为代数形式的亚临界情况（$p<d$）不同，作者证明了对于 $p=d$，其临界缩放是指数级的。这导致了一种性质不同的极限行为，其中有效容量取决于对数缩放参数 $\tau$。
• 共形不变性： 作者指出，在临界情况下，由于 $d$-能量的共形不变性，腔体的具体几何细节（除了其一致的非退化性之外）不会影响领先项的渐近行为。这与需要要求腔体容量值必须一致的亚临界情况形成对比。
• 显式构造： 本文提供了一个能够高精度逼近解（在能量和逐点范数下）的显式假设函数，提供了从微观腔体结构到宏观均匀化势的完整描述。

这项工作扩展了作者之前的研究结果（特别是关于 $p<d$ 的亚临界情况），并填补了边界指数 $p=d$ 的空白，为理解分布在球面上的消失孔洞如何影响临界维数下非线性椭圆方程的全局解提供了严谨的框架。