Singular semilinear elliptic equations in nondivergence form

技术摘要：非散度型奇异半线性椭圆方程

问题陈述
本文研究了奇异半线性 Dirichlet 问题解的存在性与唯一性：
$$\begin{cases} -Pu = \frac{f}{u^\gamma} & \text{在 } \Omega \text{ 内}, \\ u \equiv 0 & \text{在 } \partial\Omega \text{ 上}, \end{cases}$$
其中 $\Omega \subset \mathbb{R}^d$ ($d \geq 2$) 是有界区域，$f \in L^1(\Omega)$ 是非负且非平凡的，且 $\gamma > 0$。算子 $P$ 是一个二阶一致椭圆微分算子，具有非散度形式：
$$Pu(x) = \sum_{i,j=1}^d a_{ij}(x) \frac{\partial^2 u}{\partial x_i \partial x_j},$$
其中对称矩阵 $A(x) = [a_{ij}(x)]$ 满足一致椭圆性条件。奇异性源于项 $u^{-\gamma}$，该项在 $u=0$ 处无定义。边界条件被解释为在边界附近迹消失的意义：$\lim_{\epsilon \to 0} \frac{1}{\epsilon} \int_{\Omega_\epsilon} |u| = 0$，其中 $\Omega_\epsilon = \{x \in \Omega : \delta_\Omega(x) < \epsilon\}$。

方法论与技术框架
作者采用了近似技术、非线性不等式和格林函数估计相结合的方法。证明根据系数 $A$ 和解 $u$ 的正则性区分了两种解的概念：

1. 弱解：定义于 $A \in C^1(\Omega)$，此时方程在分布意义下成立，采用形式 $\langle -Pu, \phi \rangle = \int_\Omega (A\nabla u \cdot \nabla \phi + (\text{div} A \cdot \nabla u)\phi)$。
2. 分布解：定义于 $A \in C^2(\Omega)$，此时方程通过伴随算子 $P^*$ 成立：$\langle -Pu, \phi \rangle = -\int_\Omega u P^* \phi$。

主要贡献与结果

1. 存在性与正则性（定理 1.1）
在假设 $\Omega \in C^{1,1}$，$A \in C^1(\Omega)$，且漂移项 $\text{div} A$ 满足特定的小性条件（不等式 1.5）下，作者证明了弱解 $u \in W^{1,1}_{\text{loc}}(\Omega) \cap L^1(\Omega)$ 的存在性。

• 正则性：解表现出依赖于 $\gamma$ 的特定 Sobolev 正则性：
  • 当 $\gamma \geq 2$ 时，$u \in H^1_{\text{loc}}(\Omega)$。
  • 当 $\gamma \in (1, 2)$ 且 $1 \leq q < 2/\gamma$ 时，$u \in H^1_{\text{loc}}(\Omega) \cap W^{1,q}_0(\Omega)$。
  • 当 $\gamma = 1$ 时，$u \in H^1_0(\Omega)$。
  • 当 $\gamma \in (0, 1)$ 时，$u \in W^{1,q}_0(\Omega)$，其中 $q$ 为特定指数。
  • 关键在于，对于所有 $\gamma > 0$，均有 $u^{(\gamma+1)/2} \in H^1_0(\Omega)$。
• 正性：解在 $\Omega$ 内严格为正，并满足边界行为 $u \gtrsim \delta_\Omega$。

2. 唯一性（定理 1.2）
假设更高的正则性（$\Omega \in C^2$ 且 $A \in C^2(\Omega)$），作者建立了 $L^1(\Omega)$ 中分布解的唯一性。该结果适用于比存在性定理中考虑的更广泛的解类，因为它先验地不假设 $u \in W^{1,1}_{\text{loc}}(\Omega)$。

方法论创新
本文强调了与以往关于散度型算子（例如 Boccardo 和 Orsina [6]）的研究相比的三个主要技术差异：

• 非线性 Gagliardo–Nirenberg 不等式：存在性证明利用了最近推导出的先验估计，形式为 $\int_\Omega \|\nabla u\|_A^2 u^{\gamma-1} \leq C \int_\Omega |Pu|u^\gamma$。这些估计是从 [30] 中改编以适应非散度情形的。
• 格林函数估计：利用格林函数 $G_\Omega$ 及其伴随算子的尖锐估计来推导近似解的精确衰减率，这对于建立特定的 Sobolev 正则性类至关重要。
• Kato 型不等式：建立了一个针对非散度形式强解的新颖 Kato 型不等式（引理 3.5），以证明近似序列的单调性。对于唯一性，推导了一个针对极弱解的 Kato 型不等式（引理 4.10），从而克服了解空间中正则性缺失的问题。

意义与范围
作者将这项工作定位为 Emden–Fowler 方程在非散度型算子情形下的多维对应物。

• 新颖性：虽然散度型奇异半线性方程已得到充分研究，但本文处理的是 $P$ 为非散度型算子的情形，这类算子是随机过程理论和扩散生成元理论的核心。
• 技术成就：本文证明了在算子系数为 $C^1$ 的情况下可以证明存在性，而唯一性则需要 $C^2$ 系数以处理伴随算子的结构并扩展测试函数类。
• 局限性与未解决问题：作者指出，定理 1.1(b) 中建立的 $\gamma \in (1, 2)$ 情形的全局 Sobolev 正则性目前对于散度型算子尚不清楚。他们指出了关于一般非线性项 $\tau(u)$、系数正则性较低的算子以及非局部算子（例如分数阶拉普拉斯算子）的未解决问题。

证明依赖于通过施auder 不动点定理对正则化问题构造解，随后利用序列的单调性和推导出的先验估计取极限。唯一性证明涉及将测试函数类扩展到 $C^{1,P}_0(\Omega)$，并利用源项 $f/u^\gamma$ 关于距离函数 $\delta_\Omega$ 的加权可积性。