形式陈述
自由空间的点源核知道源在哪里,却不知道墙在哪里。怎样让同一个点源响应同时服从指定区域的零边界?
设 Ω ⊂ R d 是有界光滑区域,d ≥ 2 。固定 y ∈ Ω ,Dirichlet Green 函数 G Ω ( x , y ) 满足
于 − Δ x G Ω ( x , y ) = δ y 于 Ω , G Ω ( x , y ) = 0 ( x ∈ ∂ Ω ) . 方程在分布意义下理解。源点附近保留与基本解 公理库 Laplace 算子的基本解 Fundamental solution of the Laplacian · Newtonian kernel 在全空间表示单位点源响应的奇异核;通过小球通量确定常数,并在分布意义下满足负 Laplacian 等于 Dirac 质量。 Φ ( x − y ) 相同的奇性,而差
h y ( x ) = Φ ( x − y ) − G Ω ( x , y ) 在区域内部是光滑调和函数 公理库 调和函数与平均值性质 Harmonic function · Mean value property Laplace 方程的无源解等于每个内含球或球面的平均;这一局部平均规律也能反过来识别调和性。 ,在边界取值 Φ ( x − y ) 。因此构造可以分两步:保留点源,再减去专门消掉边界值的调和校正。
在上述光滑域中,这个核存在且唯一;一般存在性需要 Dirichlet 问题的解理论。本页在三维球上直接构造,并证明表示、正性与对称性。它不是只给一个形式逆符号,也不是把自由空间卷积截到区域上就算完成边界条件。
无界区域还要选择远场行为
半空间等无界开集也可使用同样的局部点源条件:固定 y ,要求 − Δ x G ( ⋅ , y ) = δ y 、边界值为零,且 Φ ( ⋅ − y ) − G ( ⋅ , y ) 在内部光滑调和。但这几项本身不再保证唯一,必须另选远离源点时的衰减或增长类别。例如在上半空间,若 G 满足这些条件,则对任意实数 c ,
G ~ ( x , y ) = G ( x , y ) + c x d y d 仍然满足它们:新增项对 x 调和、在 x d = 0 时为零,也不改变源点奇性;它还关于 x , y 对称。可见即使再要求对称性,仍不能省略远场选择。半空间镜像构造 公理库 Laplace 方程的镜像法 Method of images for Laplace equation · Image charge method 借平面反射把虚拟源放到区域外,以相反或相同符号分别满足零 Dirichlet 或零 Neumann 边界。 选择一个显式核,并分别核验点源、边界和远场;前面的有界域存在唯一性不被直接搬到这个情形。
直觉
可以把 G Ω ( x , y ) 看成在 y 放一单位源时,位置 x 的响应。自由空间核只包含源的局部规律;校正项则包含整面边界的影响。换一个区域,局部奇性不变,校正项却会改变。
因此 G 有两个空间变量:一个指观测点,一个指源点。即使两点距离相同,它们靠墙的方式不同,响应也可能不同。只有全空间这种平移不变环境,核才自然只依赖差 x − y 。
例子与边界
三维单位球的反演源
取 Ω = B ( 0 , 1 ) ⊂ R 3 ,0 < | y | < 1 ,定义球外点 y ∗ = y / | y | 2 。令
G ( x , y ) = 1 4 π | x − y | − 1 4 π | y | | x − y ∗ | . 第二个奇点在球外,所以它作为 x 的函数在整个球内调和。若 | x | = 1 ,直接计算距离得到
| x − y ∗ | = | x − y | | y | , 两项恰好抵消。第一项携带单位源,第二项没有内部源,因此这个显式式子同时通过方程与边界两道检查。
在 y = 0 的极限,核简化为
G ( x , 0 ) = 1 4 π ( 1 | x | − 1 ) . 减去的常数非常重要:只有 1 / ( 4 π | x | ) 时,球面上的值是 1 / ( 4 π ) ,并不满足零边界。
一个体源,直接算出中心值
设单位球内 − Δ u = 1 、边界 u = 0 。由核的对称性和球坐标,
u ( 0 ) = ∫ B ( 0 , 1 ) G ( 0 , y ) d y = ∫ 0 1 ( 1 r − 1 ) r 2 d r = 1 6 . 整个解为 u ( x ) = ( 1 − | x | 2 ) / 6 ,因为 − Δ u = 1 且边界为零。点源叠加和直接微分给出相同中心值。
这与全空间单位球体源不同:后者在球外仍有远场势,球内中心值是 1 / 2 。边界条件确实改变了解,而不是给同一个公式换个名称。
推论与应用
Green 表示公式的符号
在前述有界光滑域内,对足够光滑的 u ,若 − Δ u = f 于 Ω ,边界 u = g ,则
u ( x ) = ∫ Ω G ( x , y ) f ( y ) d y − ∫ ∂ Ω g ( y ) ∂ n y G ( x , y ) d S ( y ) . 这里 n y 是区域的外法向。公式来自两次分部积分 公理库 散度定理 Divergence theorem · Gauss theorem · Gauss–Ostrogradsky theorem 区域内部散度的体积分等于其边界上向外法向通量的积分。 :把 u 与 G ( x , ⋅ ) 放进 Green 第二恒等式,绕开源点的小球,再用单位源的通量极限恢复 u ( x ) 。由于 G 的边界值为零,含 G ∂ n u 的那项消失,留下所写的负号。
零边界时只剩源项积分;零源时,− ∂ n y G 就是边界数据的权重。球中的权重回到 Poisson 核,半空间中的权重则由镜像法 公理库 Laplace 方程的镜像法 Method of images for Laplace equation · Image charge method 借平面反射把虚拟源放到区域外,以相反或相同符号分别满足零 Dirichlet 或零 Neumann 边界。 算出。
为什么点源核非负
仍在上述有界域内,固定源点 y ,从区域里挖去很小的 B ( y , ε ) 。核在剩余区域调和,外边界为零;在小球面上,基本解奇性趋于正无穷,而调和校正有界,所以足够小的内球面上核为正。最小值原理 公理库 椭圆最大值原理 Elliptic maximum principle · Maximum principle for harmonic functions 有界区域内的次调和函数由边界上界控制;若在连通内部达到全局最大值,平均不等式迫使它恒定。 给剩余区域 G ≥ 0 ,再让 ε ↓ 0 。连通域内,强原理进一步给源点以外严格为正。
因此 f ≥ 0 、零边界时,Green 表示立即给 u ≥ 0 。这与直接比较原理得到的符号结论一致。
交换源点和观测点
对上述有界域核,固定不同的 x , y ,对 G ( ⋅ , x ) 与 G ( ⋅ , y ) 在同时挖去两个小球的区域使用 Green 第二恒等式。外边界上两函数均为零;每个小球上的通量极限分别读出另一个函数在该源点的值。最终得到
G ( x , y ) = G ( y , x ) . 对称性来自算子 − Δ 和相同的零 Dirichlet 边界条件,不是任意二维核都具备的形式性质。若加入非自伴漂移项,交换两变量一般会联系到伴随算子的 Green 核,而不是原核自己。
本页核描述的是一个特定边值问题。纯 Neumann 问题保留常数零模态,其兼容性和归一化要另行处理,不能直接把本页的零边界条件换成零法向导数而保留所有存在唯一性声明。
参考资料
Stanford Math 220B,Green’s Functions ,§4.1–§4.2:调和校正、Green 表示和球/半空间构造;Lemma 13:核的对称性。本文固定 − Δ 与外法向约定,并独立核算三维单位球例子。