Self-similar solutions of the three-dimensional Muskat problem with surface tension

技术摘要：带有表面张力的三维 Muskat 问题的自相似解

问题表述
本文研究了具有表面张力的三维单相 Muskat 问题，该问题描述了多孔介质中不可压缩流体在自由上边界下的行为。界面 $\Sigma_t$ 被建模为定义在 $\mathbb{R}^2$ 上的函数 $\eta(t, x)$ 的图象。控制方程结合了 Darcy 定律、不可压缩性以及一个动态边界条件，其中压力跳变由表面张力 $\sigma \kappa(\eta)$ 平衡。该系统简化为一个关于界面高度的三阶拟线性抛物型方程：
$$\partial_t \eta = -\frac{\sigma}{\mu} G(\eta)\kappa(\eta),$$
其中 $G(\eta)$ 是 Dirichlet–Neumann 算子，$\kappa(\eta)$ 是平均曲率。

作者寻求从无界圆锥形初始数据 $\eta(0, x) = \varepsilon|x|$ 出发的自相似解。这类解具有形式 $\eta_\varepsilon(t, x) = t^{1/3}U_\varepsilon(t^{-1/3}x)$，其中剖面 $U_\varepsilon$ 满足一个包含非平移不变扩张项 $x \cdot \nabla$ 的定常方程。核心挑战在于，当初始数据是无界的且在原点处存在奇异性时，如何在小 $\varepsilon$ 情况下构造这些解；这种情况超出了标准有界或临界 Sobolev 空间中的小数据 Cauchy 理论的范围。

方法论与策略
证明过程采用了摄动法，将剖面 $U_\varepsilon$ 分解为一个领先线性项和一个修正项：
$$U_\varepsilon = \varepsilon u_* + v_\varepsilon,$$
其中 $u_*$ 是“正则化线性圆锥”，定义为 $u_* = e^{-\frac{\sigma}{\mu}|D|^3}|x|$。修正项 $v_\varepsilon$ 在特定的函数空间 $X_r = \dot{H}^r(\mathbb{R}^2) \cap \dot{H}^{r+3}(\mathbb{R}^2)$ 中求解，其中 $3/2 < r < 2$。

该方法依赖于三个主要的数学技术支柱：

1. 相似算子的线性分析：
   通过沿射线的傅里叶变换对线性算子 $L = \frac{3\sigma}{\mu}|D|^3 - x \cdot \nabla + 1$ 进行分析。作者利用加权 Hardy 估计构造了一个射线逆 $L^{-1}$。该逆算子在高频处获得三个导数增益，同时保留了圆锥形的低频结构，从而实现了在模去常数意义下对剖面方程的求逆。

2. 有利的二次抵消：
   处理具有无界数据的非线性项 $G(\eta)\kappa(\eta)$ 时，一个关键难点在于二次项的行为。作者识别了双线性形式 $B(f, g) = \frac{1}{2}(D_\eta G(0)[f]\Delta g + D_\eta G(0)[g]\Delta f)$ 符号中的一种特定代数抵消。在“高-高到低”频率相互作用区域（即输入频率很大但输出频率很小的区域），其符号满足估计 $|b(\alpha, \beta)| \lesssim |\alpha+\beta||\beta|^3$。这种额外的输出频率增益确保了低频求和的收敛性，否则对于圆锥形数据，该求和将会发散。

3. 有限阶圆锥适配的 Dirichlet–Neumann 估计：
   由于高 Sobolev 范数无法控制低频能量，标准的 Dirichlet–Neumann 估计对于无界图象失效。作者引入了混合半空间范数（$X_r$ 和 $Y_r$），这些范数仅在正水平频率上施加高 Sobolev 权重。他们建立了 Dirichlet–Neumann 算子及其形状导数在这些空间上的 tame 估计，证明了该算子在圆锥剖面的邻域内是 $C^3$ 的。这使得控制非线性展开中的三次余项成为可能。

主要结果
论文确立了以下主要定理（定理 1.3）：

• 存在性： 对于足够小的 $\varepsilon$，存在一个带有表面张力的 Muskat 问题的自相似解 $\eta_\varepsilon(t, x) = t^{1/3}U_\varepsilon(t^{-1/3}x)$。
• 剖面结构： 剖面 $U_\varepsilon$ 是正则化圆锥的摄动：$U_\varepsilon = \varepsilon u_* + v_\varepsilon$，其中修正项满足 $\|v_\varepsilon\|_{X_r} \leq C \varepsilon^2$。
• 二阶展开： 作者确定了领先的二次修正项 $v^{(2)} = \frac{3\sigma}{\mu} L^{-1} B(u_*, u_*)$。当 $\varepsilon \to 0$ 时，解具有展开式 $U_\varepsilon = \varepsilon u_* + \varepsilon^2 v^{(2)} + O_{X_r}(|\varepsilon|^3)$。
• 初始迹与正则化： 解在局部一致意义下达到初始迹 $\eta_\varepsilon(0, x) = \varepsilon|x|$。尽管在 $t=0$ 时存在圆锥奇异性，但对于任何 $t > 0$，解都是瞬时光滑的；剖面 $U_\varepsilon$ 是 $C^1$ 的，且奇异性在 $O(t^{1/3})$ 的空间尺度上被“圆滑化”。

意义与主张
本文声称提供了第一个关于带有表面张力的三维 Muskat 问题、从精确圆锥初始数据出发的自相似解的构造。其意义在于克服了标准 Cauchy 理论的局限性，该理论通常要求有界或临界 Sobolev 初始数据。通过分离无界圆锥成分并利用专门设计的适配于问题几何结构的函数空间，作者证明了圆锥奇异性会被表面张力瞬时正则化。

这项工作将先前关于重力驱动 Muskat 问题（其主算子是一阶）的研究扩展到了表面张力情况（三阶算子）。作者指出，虽然所使用的代数抵消是维度无关的，但特定 Sobolev 范数的选择和参数 $r$ 的范围取决于空间维度。论文得出结论，只要初始圆锥足够小，圆锥奇异性就不会成为自相似情形下适定性的障碍。