UA PHYS515A 电磁理论IV 时变电磁场理论6 用含时Green函数求解时变电磁场问题的例子
UA PHYS515A 电磁理论IV 时变电磁场理论6 用含时Green函数求解时变电磁场问题的例子
在一个nuclear中有一些photon,photon受激产生e−,e+e^-,e^+e−,e+两个electron,其中e−e^-e−沿k^\hat kk^方向运动,e+e^+e+沿−k^-\hat k−k^方向运动,速度大小记为vvv,现在我们试图计算这个过程激发的电磁场。
回顾一下相关方程为:
(∇2−1c2∂2∂t2)Φ=−4πρ(∇2−1c2∂2∂t2)A⃗=−4πcJ⃗\left( \nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2} \right)\Phi = -4 \pi \rho \\ \left( \nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2} \right)\vec A = -\frac{4 \pi }{c}\vec J (∇2−c21∂t2∂2)Φ=−4πρ(∇2−c21∂t2∂2)A=−c4πJ
以及
B⃗=∇×A⃗E⃗=−∇Φ−1c∂A⃗∂t\vec B = \nabla \times \vec A \\ \vec E =-\nabla \Phi - \frac{1}{c}\frac{\partial \vec A}{\partial t}B=∇×AE=−∇Φ−c1∂t∂A
这个问题没有边界条件,并且是过去的事件造成未来某个时刻某个位置可以观测到电磁场,所以需要retarded solution
Φ(−)(r⃗,t)=∬G+(r⃗,t,r⃗′,t′)ρ(r⃗′,t′)d3r⃗′dt′A⃗(−)(r⃗,t)=1c∬G+(r⃗,t,r⃗′,t′)J⃗(r⃗′,t′)d3r⃗′dt′\Phi^{(-)}(\vec r ,t)=\iint G^{+}(\vec r, t , \vec r' ,t') \rho(\vec r ' ,t ')d^3 \vec r' dt' \\ \vec A^{(-)}(\vec r ,t)=\frac{1}{c}\iint G^{+}(\vec r, t , \vec r' ,t') \vec J(\vec r ' ,t ')d^3 \vec r' dt'Φ(−)(r,t)=∬G+(r,t,r′,t′)ρ(r′,t′)d3r′dt′A(−)(r,t)=c1∬G+(r,t,r′,t′)J(r′,t′)d3r′dt′
Green函数为
G±(r,t−t′)=δ(t−t′−rc)rG^{\pm}(r,t-t')=\frac{\delta(t-t' - \frac{r}{c})}{r}G±(r,t−t′)=rδ(t−t′−cr)
下面计算电磁场的source:
ρ(r⃗′,t′)=−eδ(3)(r⃗′−(0,0,vt′))+eδ(3)(r⃗′−(0,0,−vt′))J⃗(r⃗′,t′)=−evδ(3)(r⃗′−(0,0,vt′))k^+evδ(3)δ(3)(r⃗′−(0,0,−vt′))(−k^)\rho(\vec r',t')=-e\delta^{(3)}(\vec r'-(0,0,vt'))+e\delta^{(3)}(\vec r'-(0,0,-vt')) \\ \vec J(\vec r',t')=-ev\delta^{(3)}(\vec r' - (0,0,vt'))\hat k+ \\ ev\delta^{(3)}\delta^{(3)}(\vec r'-(0,0,-vt'))(-\hat k)ρ(r′,t′)=−eδ(3)(r′−(0,0,vt′))+eδ(3)(r′−(0,0,−vt′))J(r′,t′)=−evδ(3)(r′−(0,0,vt′))k^+evδ(3)δ(3)(r′−(0,0,−vt′))(−k^)
然后计算积分
∬G+(r⃗,t,r⃗′,t′)ρ(r⃗′,t′)d3r⃗′dt′=∬δ(t−t′−rc)r[−eδ(3)(r⃗′−(0,0,vt′))+eδ(3)(r⃗′−(0,0,−vt′))]dt′d3r⃗′=∫d3r⃗′r[−eδ(3)(r⃗′−(0,0,vt−vr/c))+eδ(3)(r⃗′−(0,0,−vt+vr/c))]=−ex2+y2+(z−z0)2+ex2+y2+(z−z02)z0=vt−vcx2+y2+(z−z0)2\iint G^{+}(\vec r, t , \vec r' ,t') \rho(\vec r ' ,t ')d^3 \vec r' dt' \\ = \iint\frac{\delta(t-t' - \frac{r}{c})}{r}[-e\delta^{(3)}(\vec r'-(0,0,vt'))+e\delta^{(3)}(\vec r'-(0,0,-vt'))]dt'd^3 \vec r' \\ = \int \frac{d^3 \vec r'}{r}[-e\delta^{(3)}(\vec r'-(0,0,vt-vr/c))+e\delta^{(3)}(\vec r'-(0,0,-vt+vr/c))] \\ = -\frac{e}{\sqrt{x^2+y^2+(z-z_0)^2}}+\frac{e}{\sqrt{x^2+y^2+(z-z_0^2)}} \\ z_0 = vt - \frac{v}{c}\sqrt{x^2+y^2+(z-z_0)^2}∬G+(r,t,r′,t′)ρ(r′,t′)d3r′dt′=∬rδ(t−t′−cr)[−eδ(3)(r′−(0,0,vt′))+eδ(3)(r′−(0,0,−vt′))]dt′d3r′=∫rd3r′[−eδ(3)(r′−(0,0,vt−vr/c))+eδ(3)(r′−(0,0,−vt+vr/c))]=−x2+y2+(z−z0)2e+x2+y2+(z−z02)ez0=vt−cvx2+y2+(z−z0)2
基于这个结果可以确定电磁场的标量势,类似的操作可以计算出向量势,
A⃗(r⃗,t)=−1cevk^x2+y2+(z−z0)2−1cevk^x2+y2+(z+z0)2\vec A(\vec r,t)=-\frac{1}{c}\frac{ev\hat k}{\sqrt{x^2+y^2+(z-z_0)^2}}-\frac{1}{c}\frac{ev\hat k}{\sqrt{x^2+y^2+(z+z_0)^2}}A(r,t)=−c1x2+y2+(z−z0)2evk^−c1x2+y2+(z+z0)2evk^
然后计算磁感应强度和电场强度太难算了就懒得写了。
UA PHYS515A 电磁理论IV 时变电磁场理论6 用含时Green函数求解时变电磁场问题的例子相关推荐
- UA PHYS515A 电磁理论V 电磁波与辐射7 运动点电荷的辐射
UA PHYS515A 电磁理论V 电磁波与辐射7 运动点电荷的辐射 实际问题中辐射的source都是比较复杂的charge density与currency density,但作为比较简单直观易于理 ...
- UA PHYS515A 电磁理论IV 时变电磁场理论2 Helmholtz方程与含时的Green函数
UA PHYS515A 电磁理论IV 时变电磁场理论2 Helmholtz方程与含时的Green函数 上一讲的末尾我们介绍了Lorentz Gauge下的含时麦克斯韦方程: (∇2−1c2∂2∂t2) ...
- UA PHYS515A 电磁理论IV 时变电磁场理论3 电磁场的能量守恒
UA PHYS515A 电磁理论IV 时变电磁场理论3 电磁场的能量守恒 时变电磁场的传播成为电磁波,它可以携带能量,这一讲我们讨论电磁场的能量.回忆一下电荷守恒,麦克斯韦在推导他的方程组的时候,用电 ...
- UA PHYS515A 电磁理论IV 时变电磁场理论5 电磁场的角动量
UA PHYS515A 电磁理论IV 时变电磁场理论5 电磁场的角动量 回顾一下angle momentum的定义 L⃗=r⃗×p⃗\vec L= \vec r \times \vec pL=r×p ...
- UA PHYS515A 电磁理论V 电磁波与辐射9 简单辐射系统
UA PHYS515A 电磁理论V 电磁波与辐射9 简单辐射系统 前文讨论了单个直线运动的带电粒子的辐射,但在实践中只有在实验室才能观察到这种现象,应用中遇到的情况会更加复杂,比如多个粒子做直线运动, ...
- UA PHYS515A 电磁理论V 电磁波与辐射8 单个粒子的辐射 匀速运动与匀加速运动的情况
UA PHYS515A 电磁理论V 电磁波与辐射8 单个粒子的辐射 匀速运动与匀加速运动的情况 单个粒子的辐射场满足: E=q((n^−β⃗)(1−β⃗2)(1−n^⋅β⃗)3R2+n^×[n^−β⃗ ...
- UA PHYS515A 电磁理论V 电磁波与辐射4 反射与折射
UA PHYS515A 电磁理论V 电磁波与辐射4 反射与折射 假设平面z=0z=0z=0是两种不导电介质的交界面,z<0z<0z<0的空间中介质的介电常数与磁导率为ϵ1,μ1\ep ...
- UA PHYS515A 电磁理论V 电磁波与辐射3 偏振
UA PHYS515A 电磁理论V 电磁波与辐射3 偏振 之前谈到过电动力学讨论的就是作为源头的电荷与电流如何产生电磁场,以及由源头激发出的电磁场又如何反作用于这些电荷与电流以及存在于介质中的电荷.因 ...
- UA PHYS515A 电磁理论V 电磁波与辐射5 电磁波在介质中的传播
UA PHYS515A 电磁理论V 电磁波与辐射5 电磁波在介质中的传播 在介绍麦克斯韦方程的时候,我们提到过 D⃗=E⃗+4πP⃗\vec D = \vec E + 4 \pi \vec PD=E+ ...
最新文章
- 【python教程入门学习】Django应用开发框架的优点和缺点!!
- php测试网络通不通函数,PHP检查库或函数是否可用的方法
- Nagios 安装配置
- 关于Mybatis,我总结了 10 种通用的写法
- python 网络编程 问题记录
- php ignore special characters,PHP htmlspecialchars() 函數--防注入字符轉義函數
- QT的QDialogButtonBox类的使用
- ugui unity 取消选择_关于Unity中的UGUI优化,你可能遇到这些问题
- js生成图片_网页骨架屏自动生成方案(dps)
- 以太坊2.0合约质押新增4.23万ETH
- Kotlin入门(23)适配器的进阶表达
- eclipse哪个版本最稳定_新版鹿鼎记开播啦,你最喜欢哪个版本?最还原原著的原来是他...
- 【证明】—— 矩阵秩的相关证明
- 毕业论文页码及目录设置方法
- [ISITDTU 2019]EasyPHP
- 手写代码详解Android Hook入门demo,android应用案例开发大全第四版源码
- 超好用的图片压缩网站
- 远程办公:在外随时远程控制公司电脑
- 【Arduino实验11 LED点阵屏显示】
- opencv-图像金字塔