形式陈述
一个函数在原点以外处处无源,怎么还能在原点携带一个单位源?关键是把方程放到分布 公理库 Schwartz 分布 Schwartz distribution · Generalized function · 广义函数 以紧支撑光滑测试函数上的连续线性泛函容纳普通函数、跳跃导数和点源,并用显式 Green 核连接分布方程与零边界能量解。 意义下,并检查包围奇点的通量。
在 R d 、d ≥ 2 中,记 ω d − 1 = | S d − 1 | 。对算子 − Δ ,取基本解
Φ ( x ) = { | x | 2 − d ( d − 2 ) ω d − 1 , d ≥ 3 , − 1 2 π log | x | , d = 2. 它局部可积,在 x ≠ 0 时光滑且 Δ Φ = 0 ,而在整个空间满足
− Δ Φ = δ 0 . 具体含义是:对每个 φ ∈ C c ∞ ( R d ) ,
∫ Φ ( x ) ( − Δ φ ( x ) ) d x = φ ( 0 ) . 三维最熟悉的公式为 Φ ( x ) = 1 / ( 4 π | x | ) 。本页始终采用 − Δ 的符号;若方程写成 Δ u = f ,基本解要取相反数。这个核没有施加某个有界区域的边界条件。
直觉
离开源点,通量从一层球面传到下一层,没有新增或损失。半径为 r 的球面面积按 r d − 1 增长,因此场的径向强度必须按 r 1 − d 衰减,才能让总通量固定为一。
对这种场沿半径积分,就得到 r 2 − d ;二维积分的是 1 / r ,因此出现对数。维数改变的不只是一个常数,而是远场和奇点的函数形状。
例子与边界
三维单位源的常数不能靠记忆猜
设候选势为 Φ ( r ) = C / r 。在半径 r 的球面上,外向单位法线就是径向方向,∂ ν Φ = − C / r 2 。因此
− ∫ ∂ B ( 0 , r ) ∂ ν Φ d S = C r 2 4 π r 2 = 4 π C . 要让通量等于一,必须取 C = 1 / ( 4 π ) 。只写 1 / r 会多出 4 π 的点源质量。
把单位点源铺成一个实心球
在三维中,考虑 − Δ u = 1 B ( 0 , 1 ) ,并要求 u → 0 于无穷远。径向方程是
− ( r 2 u ′ ( r ) ) ′ = r 2 1 r < 1 . 原点处不出现额外点源,迫使内部积分常数为零;内外解及一阶导数在 r = 1 连续,得到
u ( r ) = { 1 2 − r 2 6 , 0 ≤ r ≤ 1 , 1 3 r , r ≥ 1. 源的总质量为 4 π / 3 ,远处势恰是 ( 4 π / 3 ) Φ ( r ) = 1 / ( 3 r ) 。边界处没有额外 Dirac 面质量,因为法向一阶导数连续;二阶导数的跳跃只反映源密度从一变到零。
奇点之外的方程不能决定奇点本身
Φ 在原点之外满足 − Δ Φ = 0 ,并不允许把这个等式直接推广为整个空间的零分布。二阶导数在奇点附近不绝对可积,绕开奇点时留下的球面项不会消失。另一方面,原点的任意点值不影响局部可积函数所定义的分布,真正携带源的是邻域中的奇异形状。
二维的 − log r / ( 2 π ) 在无穷远不趋零。给它加常数不改变方程;若要求指定远场归一化,必须与维数和总源质量一起检查。
推论与应用
先由径向方程寻找候选核
对径向函数 v ( r ) ,
Δ v = v ″ + d − 1 r v ′ = r 1 − d ( r d − 1 v ′ ) ′ . 在 r > 0 上令它为零,得到 v ′ = C r 1 − d 。积分并按单位通量规范化,就产生上述两个维数范围的公式。它们在原点附近都局部可积:高维体积积分只含 ∫ 0 ε r d r ,二维则含 ∫ 0 ε r | log r | d r 。
分布恒等式来自被挖掉的小球
对紧支撑测试函数 φ ,在其支撑外选一个大球,再挖去半径 ε 的原点小球。穿孔区域内 Δ Φ = 0 。两次应用散度定理 公理库 散度定理 Divergence theorem · Gauss theorem · Gauss–Ostrogradsky theorem 区域内部散度的体积分等于其边界上向外法向通量的积分。 ,得到内球面贡献
∫ | x | > ε Φ ( − Δ φ ) d x = ∫ | x | = ε ( φ ∂ n Φ − Φ ∂ n φ ) d S . 这里 n 是穿孔区域的外法向,在内球面上指向原点,所以 ∂ n Φ = − Φ ′ ( ε ) 为正。第一项是球面上的 φ 平均,趋向 φ ( 0 ) ;第二项在 d ≥ 3 时为 O ( ε ) ,在二维为 O ( ε | log ε | ) ,均趋零。这正好证明 − Δ Φ = δ 0 。
若 f ∈ C c ∞ ,叠加点源给 Newton 势 u = Φ ∗ f ,并在分布意义下满足 − Δ u = f ;卷积可把导数转移到光滑紧支撑的 f 上,因而此时解也光滑。要把它用于有界域的Poisson 边值问题 公理库 偏微分方程 Partial differential equation · PDE 含未知多元函数及其偏导数的方程。 ,还需要减去合适的调和校正;这一步形成Dirichlet Green 函数 公理库 Dirichlet Green 函数 Dirichlet Green function 在自由空间点源核上减去调和边界校正,得到同时满足单位点源方程和零 Dirichlet 边界的区域核。 。
参考资料