UA PHYS515 电磁理论II 静电场问题5 用Green函数法求解interior Dirichlet问题的例子
UA PHYS515 电磁理论II 静电场问题5 用Green函数法求解interior Dirichlet问题的例子
例2
均匀金属空心外壳厚度可忽略的接地球球心位于原点,半径为aaa,用球坐标(r,θ,ϕ)(r,\theta,\phi)(r,θ,ϕ)描述,球面上边界条件为Φ=Φ(a,θ,ϕ)\Phi=\Phi(a,\theta,\phi)Φ=Φ(a,θ,ϕ),计算球内的电场。
解
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条件为Φ(r⃗)=0,∀r⃗∈S\Phi(\vec r)=0,\forall \vec r \in SΦ(r)=0,∀r∈S;写出Green函数:
G(r⃗,r⃗′)=1∣r⃗−r⃗′∣+F(r⃗,r⃗′)G(\vec r,\vec r')=\frac{1}{|\vec r - \vec r'|}+F(\vec r, \vec r')G(r,r′)=∣r−r′∣1+F(r,r′)
为简化起见,我们假设沿r⃗′\vec r'r′的方向存在一个image charge,我们假设它的电荷量为q′q'q′,用n^′\hat n'n^′表示与r⃗′\vec r'r′平行的SSS的外法向,则image charge的位置可以表示为r⃗′′=r′′n^′\vec r'' = r''\hat n'r′′=r′′n^′;假设测试电荷的位置为r⃗\vec rr,与它平行的SSS的外法向为n^\hat nn^,假设n^\hat nn^与n^′\hat n'n^′的夹角为γ\gammaγ;于是
G(r⃗,r⃗′)=1∣r⃗−r⃗′∣+q′∣r⃗−r⃗′′∣=1∣rn^−r′n^′∣+q′∣rn^−r′′n^∣G(\vec r, \vec r')=\frac{1}{|\vec r - \vec r'|}+\frac{q'}{|\vec r - \vec r''|} = \frac{1}{|r \hat n-r'\hat n'|}+\frac{q'}{|r \hat n-r''\hat n|}G(r,r′)=∣r−r′∣1+∣r−r′′∣q′=∣rn^−r′n^′∣1+∣rn^−r′′n^∣q′
其中q′q'q′与r′′r''r′′是未知参数,我们需要用image charge的性质,解出这两个参数。根据G(r⃗,r⃗′)∣r⃗∈S=0G(\vec r,\vec r')|_{\vec r \in S}=0G(r,r′)∣r∈S=0,我们可以得到:
1∣rn^−r′n^′∣+q′∣rn^−r′′n^∣=1r∣n^−r′rn^′∣+q′r′′∣rr′′n^−n^′∣=0r=a\frac{1}{|r \hat n-r'\hat n'|}+\frac{q'}{|r \hat n-r''\hat n|} = \frac{1}{r|\hat n - \frac{r'}{r}\hat n'|}+\frac{q'}{r''| \frac{r}{r''}\hat n-\hat n'|}=0 \\ r=a∣rn^−r′n^′∣1+∣rn^−r′′n^∣q′=r∣n^−rr′n^′∣1+r′′∣r′′rn^−n^′∣q′=0r=a
一种可行的解是
{q′r′′=−1r=−1arr′′=ar′′=r′r=r′a⇒{q′=−ar′2r′′=a2r′\begin{cases} \frac{q'}{r''}=-\frac{1}{r} = -\frac{1}{a} \\ \frac{r}{r''} = \frac{a}{r''}=\frac{r'}{r}=\frac{r'}{a} \end{cases} \Rightarrow \begin{cases} q' = -\frac{a}{r'^2} \\ r'' = \frac{a^2}{r'} \end{cases}{r′′q′=−r1=−a1r′′r=r′′a=rr′=ar′⇒{q′=−r′2ar′′=r′a2
需要注意的是r⃗′\vec r'r′的方向就是边界的一个外法线方向n^′\hat n'n^′,测试电荷的位置r⃗\vec rr的方向也是边界的一个外法线方向n^\hat nn^,记这两个外法线方向夹角为γ\gammaγ,则在球坐标中:
r⃗=r(sinθcosϕ,sinθ,sinϕ,cosθ)cosγ=cosθcosθ′+sinθsinθ′cos(ϕ−ϕ′)\vec r = r(\sin \theta\cos \phi,\sin \theta,\sin \phi,\cos \theta) \\ \cos \gamma= \cos \theta \cos \theta'+\sin \theta \sin \theta' \cos (\phi-\phi')r=r(sinθcosϕ,sinθ,sinϕ,cosθ)cosγ=cosθcosθ′+sinθsinθ′cos(ϕ−ϕ′)
现在我们就可以把Green函数在球坐标系下的表达式写出来了,
G=1r2+r′2−2rr′cosγ−1r2r′2a2+a2−2rr′cosγG =\frac{1}{\sqrt{r^2+r'^2-2rr'\cos \gamma}}-\frac{1}{\sqrt{\frac{r^2r'^2}{a^2}+a^2-2rr'\cos \gamma}}G=r2+r′2−2rr′cosγ1−a2r2r′2+a2−2rr′cosγ1
Dirichlet问题的积分解:
Φ(r⃗)=∫Vρ(r⃗′)G(r⃗,r⃗′)dx′dy′dz′−14π∮S(V)Φ(r⃗′)∂G∂ndS\Phi(\vec r) = \int_V \rho(\vec r')G(\vec r,\vec r')dx'dy'dz'-\frac{1}{4\pi}\oint_{S(V)} \Phi(\vec r')\frac{\partial G}{\partial n}dSΦ(r)=∫Vρ(r′)G(r,r′)dx′dy′dz′−4π1∮S(V)Φ(r′)∂n∂GdS
因为这个题目球的内部没有source,所以第一项为零,于是我们只需计算
∂G∂n=∂G∂r′=−r′−rcosγ(r2+r′2−2rr′cosγ)3/2+r2r′a2−rcosγ(r2r′2a2+a2−2rr′cosγ)3/2∂G∂n∣r′=a=r2−a2a(r2+a2−2arcosγ)3/2\frac{\partial G}{\partial n}=\frac{\partial G}{\partial r'}=-\frac{r'-r \cos \gamma}{(r^2+r'^2-2rr' \cos \gamma)^{3/2}}\\ + \frac{\frac{r^2r'}{a^2}-r\cos \gamma}{(\frac{r^2r'^2}{a^2}+a^2-2rr'\cos \gamma)^{3/2}} \\ \frac{\partial G}{\partial n}|_{r'=a} = \frac{r^2-a^2}{a(r^2+a^2-2ar\cos \gamma)^{3/2}}∂n∂G=∂r′∂G=−(r2+r′2−2rr′cosγ)3/2r′−rcosγ+(a2r2r′2+a2−2rr′cosγ)3/2a2r2r′−rcosγ∂n∂G∣r′=a=a(r2+a2−2arcosγ)3/2r2−a2
因此结果为
Φ(r,θ,ϕ)=−14π∮SΦ(a′,θ′ϕ′)r2−a2a(r2+a2−2arcosγ)3/2dS\Phi(r,\theta,\phi)=-\frac{1}{4 \pi}\oint_S \Phi(a',\theta'\phi')\frac{r^2-a^2}{a(r^2+a^2-2ar\cos \gamma)^{3/2}}dSΦ(r,θ,ϕ)=−4π1∮SΦ(a′,θ′ϕ′)a(r2+a2−2arcosγ)3/2r2−a2dS
其中dS=a2sinθdθdϕdS = a^2\sin \theta d \theta d\phidS=a2sinθdθdϕ, 因此
Φ(r,θ,ϕ)=−14π∬Φ(a′,θ′ϕ′)a(r2−a2)sinθ(r2+a2−2arcosγ)3/2dθdϕ\Phi(r,\theta,\phi)=-\frac{1}{4 \pi}\iint \Phi(a',\theta'\phi')\frac{a(r^2-a^2)\sin \theta}{(r^2+a^2-2ar\cos \gamma)^{3/2}}d\theta d\phiΦ(r,θ,ϕ)=−4π1∬Φ(a′,θ′ϕ′)(r2+a2−2arcosγ)3/2a(r2−a2)sinθdθdϕ
UA PHYS515 电磁理论II 静电场问题5 用Green函数法求解interior Dirichlet问题的例子相关推荐
- UA PHYS515 电磁理论II 静电场问题4 用Green函数法求解Dirichlet问题
UA PHYS515 电磁理论II 静电场问题4 用Green函数法求解Dirichlet问题 上一讲我们讨论过Dirichlet问题的积分解: Φ(r⃗)=∫Vρ(r⃗′)G(r⃗,r⃗′)dx′d ...
- UA PHYS515 电磁理论II 静电场问题6 正交函数系简介
UA PHYS515 电磁理论II 静电场问题6 正交函数系简介 完备标准正交函数系 常用的正交系 正交系与Laplace方程 完备标准正交函数系 Poisson方程的解可以用正交函数系表示,在不同的 ...
- UA PHYS515 电磁理论II 静电场问题3 边值问题及其解的唯一性
UA PHYS515 电磁理论II 静电场问题3 边值问题及其解的唯一性 Dirichlet问题与Neumann问题 边值问题解的唯一性 Green函数解释Image charge method 上一 ...
- UA PHYS515 电磁理论II 静电场问题2 电荷与静电场的几何: Green函数法的物理背景
UA PHYS515 电磁理论II 静电场问题2 电荷与静电场的几何: Green函数法的物理背景 单个电荷形成的静电场 Green函数的一些数学结果 Green恒等式与Green定理 Green定理 ...
- UA PHYS515 电磁理论II 静电场问题1 对称性与Image Charge Method
UA PHYS515 电磁理论II 静电场问题1 对称性与Image Charge Method 对称法与Image Charge Method的理论基础 例1:接地导电板外的自由电荷 例2:均匀电场 ...
- UA PHYS515 电磁理论I 麦克斯韦方程组基础1 库仑定律与毕奥-萨伐尔定律
UA PHYS515 电磁理论I 麦克斯韦方程组基础1 库仑定律与毕奥-萨伐尔定律 Coulomb定律与电场 Biot-Savart定律与磁场 写在前面 电磁理论是物理学中最优美的理论之一,因为它完备 ...
- UA PHYS515A 电磁理论II 静电学问题的一个例子
UA PHYS515A 电磁理论II 静电学问题的一个例子 例 假设有一个中空球形导体,中空部分也是一个球形,半径为aaa,球心与导体相同,导体半径为bbb:球心处有一个+q+q+q的点电荷,距离圆心 ...
- UA PHYS515 电磁理论I 麦克斯韦方程组基础5 电动力学的四类问题与对应的麦克斯韦方程
UA PHYS515 电磁理论I 麦克斯韦方程组基础5 电动力学的四类问题与对应的麦克斯韦方程 Electrostatics Magnetostatics Electromagnetic wave 一 ...
- UA PHYS515 电磁理论I 麦克斯韦方程组基础4 介质中的麦克斯韦方程
UA PHYS515 电磁理论I 麦克斯韦方程组基础4 介质中的麦克斯韦方程 推导介质中的麦克斯韦方程 电位移矢量与辅助磁场强度 推导介质中的麦克斯韦方程 前三讲我们介绍了真空中的麦克斯韦方程的建立, ...
最新文章
- 使用openssl给web站点颁发证书
- 使用JQuery筛选sharepoint日历里的超链接
- OpenStack Placement Project(Rocky to Train)
- VC:隐藏CTabCtrl标签按钮
- hasLayout和filter滤镜的细节
- 常见特征检测算法介绍
- wps怎么画网络图_尼玛,WPS竟然悄悄搞了个免费的公文写作AI...【石说】
- 关于北大陈秉乾教授的一些趣事
- android p获取通话记录_Android 底层的进程间同步机制
- 懒人神器——自制自动写字机
- java sjis_java乱码分析
- 使用OTDR进行双向光纤测试
- 学习三部曲:WHAT、HOW、WHY
- iview中Page分页组件添加首页尾页按钮
- python离线环境迁移_Python离线项目迁移部署
- excel 如何去除两列重复项
- PyTorch搭建卷积神经网络(ResNet-50网络)进行图像分类实战(附源码和数据集)
- FF新推荐的本地文件
- 关于数据存储的三道面试题,你会吗?
- python:selenium库进行网易云歌曲匹配播放