UA PHYS515 电磁理论II 静电场问题5 用Green函数法求解interior Dirichlet问题的例子

例2
均匀金属空心外壳厚度可忽略的接地球球心位于原点，半径为$a$，用球坐标$(r,\theta ,\varphi )$描述，球面上边界条件为$\Phi =\Phi (a,\theta ,\varphi )$，计算球内的电场。
解
V = { r ⃗ : ∣ r ⃗ ∣ < R } V=\{\vec r:|\vec r| < R\} V={r :∣r ∣<R}，边界为 S = { r ⃗ : ∣ r ⃗ ∣ = R } S=\{\vec r:|\vec r|=R\} S={r :∣r ∣=R}，Dirichlet条件为$\Phi (\vec{r})=0,\forall \vec{r}\in S$；写出Green函数：
$$G(\vec{r},{\vec{r}}^{\prime })=\frac{1}{|\vec{r}-{\vec{r}}^{\prime }|}+F(\vec{r},{\vec{r}}^{\prime })$$

为简化起见，我们假设沿${\vec{r}}^{\prime }$的方向存在一个image charge，我们假设它的电荷量为$q^{\prime }$，用${\hat{n}}^{\prime }$表示与${\vec{r}}^{\prime }$平行的$S$的外法向，则image charge的位置可以表示为${\vec{r}}^{\prime \prime }=r^{\prime \prime }{\hat{n}}^{\prime }$；假设测试电荷的位置为$\vec{r}$，与它平行的$S$的外法向为$\hat{n}$，假设$\hat{n}$与${\hat{n}}^{\prime }$的夹角为$\gamma$；于是
$$G(\vec{r},{\vec{r}}^{\prime })=\frac{1}{|\vec{r}-{\vec{r}}^{\prime }|}+\frac{q^{\prime }}{|\vec{r}-{\vec{r}}^{\prime \prime }|}=\frac{1}{|r\hat{n}-r^{\prime }{\hat{n}}^{\prime }|}+\frac{q^{\prime }}{|r\hat{n}-r^{\prime \prime }\hat{n}|}$$

其中$q^{\prime }$与$r^{\prime \prime }$是未知参数，我们需要用image charge的性质，解出这两个参数。根据$G(\vec{r},{\vec{r}}^{\prime }){|}_{\vec{r}\in S}=0$，我们可以得到：
$$\begin{gathered} \frac{1}{|r\hat{n}-r^{\prime }{\hat{n}}^{\prime }|}+\frac{q^{\prime }}{|r\hat{n}-r^{\prime \prime }\hat{n}|}=\frac{1}{r|\hat{n}-\frac{r^{\prime }}{r}{\hat{n}}^{\prime }|}+\frac{q^{\prime }}{r^{\prime \prime }|\frac{r}{r^{\prime \prime }}\hat{n}-{\hat{n}}^{\prime }|}=0 \\ r=a \end{gathered}$$

一种可行的解是
$$\{\begin{array}{l}\frac{q^{\prime }}{r^{\prime \prime }}=-\frac{1}{r}=-\frac{1}{a}\\ \frac{r}{r^{\prime \prime }}=\frac{a}{r^{\prime \prime }}=\frac{r^{\prime }}{r}=\frac{r^{\prime }}{a}\end{array}\Rightarrow \{\begin{array}{l}{q}^{\prime }=-\frac{a}{r^{\prime 2}}\\ r^{\prime \prime }=\frac{a^2}{r^{\prime }}\end{array}$$

需要注意的是${\vec{r}}^{\prime }$的方向就是边界的一个外法线方向${\hat{n}}^{\prime }$，测试电荷的位置$\vec{r}$的方向也是边界的一个外法线方向$\hat{n}$，记这两个外法线方向夹角为$\gamma$，则在球坐标中：
$$\begin{gathered} \vec{r}=r(\sin \theta \cos \varphi ,\sin \theta ,\sin \varphi ,\cos \theta ) \\ \cos \gamma =\cos \theta \cos {\theta }^{\prime }+\sin \theta \sin {\theta }^{\prime }\cos (\varphi -{\varphi }^{\prime }) \end{gathered}$$

现在我们就可以把Green函数在球坐标系下的表达式写出来了，
$$G=\frac{1}{\sqrt{r^2+r^{\prime 2}-2r{r}^{\prime }\cos \gamma }}-\frac{1}{\sqrt{\frac{r^2{r}^{\prime 2}}{a^2}+a^2-2r{r}^{\prime }\cos \gamma }}$$

Dirichlet问题的积分解：
$$\Phi (\vec{r})={\int }_V\rho ({\vec{r}}^{\prime })G(\vec{r},{\vec{r}}^{\prime })d{x}^{\prime }d{y}^{\prime }d{z}^{\prime }-\frac{1}{4\pi }{\oint }_{S(V)}\Phi ({\vec{r}}^{\prime })\frac{\partial G}{\partial n}dS$$

因为这个题目球的内部没有source，所以第一项为零，于是我们只需计算
$$\begin{gathered} \frac{\partial G}{\partial n}=\frac{\partial G}{\partial r^{\prime }}=-\frac{r^{\prime }-r\cos \gamma }{(r^2+r^{\prime 2}-2r{r}^{\prime }\cos \gamma {)}^{3/2}} \\ +\frac{\frac{r^2{r}^{\prime }}{a^2}-r\cos \gamma }{(\frac{r^2{r}^{\prime 2}}{a^2}+a^2-2r{r}^{\prime }\cos \gamma {)}^{3/2}} \\ \frac{\partial G}{\partial n}{|}_{r^{\prime }=a}=\frac{r^2-a^2}{a(r^2+a^2-2ar\cos \gamma {)}^{3/2}} \end{gathered}$$

因此结果为
$$\Phi (r,\theta ,\varphi )=-\frac{1}{4\pi }{\oint }_S\Phi (a^{\prime },{\theta }^{\prime }{\varphi }^{\prime })\frac{r^2-a^2}{a(r^2+a^2-2ar\cos \gamma {)}^{3/2}}dS$$

其中$dS=a^2\sin \theta d\theta d\varphi$, 因此
$$\Phi (r,\theta ,\varphi )=-\frac{1}{4\pi }\iint \Phi (a^{\prime },{\theta }^{\prime }{\varphi }^{\prime })\frac{a(r^2-a^2)\sin \theta }{(r^2+a^2-2ar\cos \gamma {)}^{3/2}}d\theta d\varphi$$