UA PHYS515A 电磁理论IV 时变电磁场理论2 Helmholtz方程与含时的Green函数
UA PHYS515A 電磁理論IV 時(shí)變電磁場理論2 Helmholtz方程與含時(shí)的Green函數(shù)
上一講的末尾我們介紹了Lorentz Gauge下的含時(shí)麥克斯韋方程:
(?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
與靜電學(xué)與靜磁學(xué)問題相比,含時(shí)的電磁場問題僅僅多了potential關(guān)于時(shí)間的二階導(dǎo)項(xiàng),因此類似的方法比如正交函數(shù)法、Green函數(shù)法在含時(shí)問題中依然適用。下面討論
第一種常用的工具是Fourier變換:
Φ(r?,t)=12π∫?∞+∞Φ(r?,w)e?iwtdw\Phi(\vec r,t)=\frac{1}{2 \pi} \int_{-\infty}^{+\infty} \Phi(\vec r,w)e^{-iwt}dwΦ(r,t)=2π1?∫?∞+∞?Φ(r,w)e?iwtdw
這樣做的好處是把時(shí)間與空間分離了成了兩個(gè)互不干擾的因式,對(duì)所有物理量都可以做Fourier變換,這樣關(guān)于Φ\PhiΦ的方程就變成了
(?2+w2c2)Φ(r?,w)=?4πρ(r?,w)\left( \nabla^2 +\frac{w^2}{c^2}\right)\Phi(\vec r,w)=-4\pi \rho(\vec r,w)(?2+c2w2?)Φ(r,w)=?4πρ(r,w)
其中www的物理意義是頻率。對(duì)A?\vec AA的所有分量我們也可以做類似的變化,并最終得到
(?2+w2c2)A?(r?,w)=?4πcJ?(r?,w)\left( \nabla^2 +\frac{w^2}{c^2}\right)\vec A(\vec r,w)=-\frac{4\pi}{c} \vec J(\vec r,w)(?2+c2w2?)A(r,w)=?c4π?J(r,w)
我們稱這種形式的方程為Helmo(UA PHYS515A 電磁理論IV 時(shí)變電磁場理論2 Fourier變換與含時(shí)的Green函數(shù))
上一講的末尾我們介紹了Lorentz Gauge下的含時(shí)麥克斯韋方程:(?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與靜電學(xué)與靜磁學(xué)問題相比,含時(shí)的電磁場問題僅僅多了potential關(guān)于時(shí)間的二階導(dǎo)項(xiàng),因此類似的方法比如正交函數(shù)法、Green函數(shù)法在含時(shí)問題中依然適用。
第一種常用的工具是Fourier變換:Φ(r?,t)=12π∫?∞+∞Φ(r?,w)e?iwtdw\Phi(\vec r,t)=\frac{1}{2 \pi} \int_{-\infty}^{+\infty} \Phi(\vec r,w)e^{-iwt}dwΦ(r,t)=2π1?∫?∞+∞?Φ(r,w)e?iwtdw這樣做的好處是把時(shí)間與空間分離了成了兩個(gè)互不干擾的因式,對(duì)所有物理量都可以做Fourier變換,這樣關(guān)于Φ\PhiΦ的方程就變成了(?2+w2c2)Φ(r?,w)=?4πρ(r?,w)\left( \nabla^2 +\frac{w^2}{c^2}\right)\Phi(\vec r,w)=-4\pi \rho(\vec r,w)(?2+c2w2?)Φ(r,w)=?4πρ(r,w)其中www的物理意義是頻率。對(duì)A?\vec AA的所有分量我們也可以做類似的變化,并最終得到(?2+w2c2)A?(r?,w)=?4πcJ?(r?,w)\left( \nabla^2 +\frac{w^2}{c^2}\right)\vec A(\vec r,w)=-\frac{4\pi}{c} \vec J(\vec r,w)(?2+c2w2?)A(r,w)=?c4π?J(r,w)我們稱這種形式的方程為Helmholtz方程。Helmholtz方程的求解方法可以完全照搬Poisson方程的求解。
第二種常用工具是含時(shí)的Green函數(shù)(time-dependent Green’s function)。用GwG_wGw?表示Green函數(shù)(關(guān)于頻率的與位置的),我們希望它滿足Helmholtz方程:
(?2+w2c2)Gw(r?,r?′)=?4πδ(3)(r??r?′)\left( \nabla^2 +\frac{w^2}{c^2}\right)G_w(\vec r,\vec r')=-4 \pi \delta^{(3)}(\vec r - \vec r ')(?2+c2w2?)Gw?(r,r′)=?4πδ(3)(r?r′)
定義r=∣r??r?′∣r=|\vec r - \vec r'|r=∣r?r′∣,則
1rd2dr2rGw+w2c2Gw=?4πδ(r)\frac{1}{r}\frac{d^2}{dr^2} rG_w +\frac{w^2}{c^2}G_w = -4 \pi \delta(r)r1?dr2d2?rGw?+c2w2?Gw?=?4πδ(r)
之所以做這個(gè)替換是因?yàn)樵陟o電學(xué)問題中,我們已經(jīng)知道Green函數(shù)雖然對(duì)電磁場的幾何的描述,但方向的作用并不大,真正影響強(qiáng)度的是與source的距離。這個(gè)方程的解是
Gw=e±ikrr,k=wcG_w=\frac{e^{\pm ikr}}{r},k=\frac{w}{c}Gw?=re±ikr?,k=cw?
分母與靜電學(xué)問題中的Green函數(shù)一樣,它表示庫侖potential,也就是potential與距離成反比;在靜電學(xué)問題中,分子為1,因?yàn)殪o態(tài)問題意味著空間中各處都遵循庫侖potential的規(guī)則,但在電磁學(xué)問題中,分子為e±ikre^{\pm ikr}e±ikr,它表示的是在空間中電磁場以波的形式進(jìn)行傳播的特征。
現(xiàn)在我們從Fourier形式的Green函數(shù)推導(dǎo)含時(shí)的Green函數(shù),
G(r?,t,r?′,t′)=12π∫?∞+∞e±ikrre?iwtdw=12π∫?∞+∞e±ik∣r??r?′∣∣r??r?′∣e?iw(t?t′)dw=12π∣r??r?′∣∫?∞+∞e?iw(t?t′)e±iwc∣r??r?′∣dw=12π∣r??r?′∣∫?∞+∞e?iw[(t?t′)?∣r??r?′∣c]dw?G±(r,t?t′)=δ(t?t′?rc)rG(\vec r,t,\vec r',t') = \frac{1}{2 \pi }\int_{-\infty}^{+\infty}\frac{e^{\pm ikr}}{r}e^{-iwt}dw \\ = \frac{1}{2 \pi }\int_{-\infty}^{+\infty}\frac{e^{\pm ik|\vec r - \vec r'|}}{|\vec r - \vec r'|}e^{-iw(t-t')}dw \\ = \frac{1}{2 \pi |\vec r - \vec r'|} \int_{-\infty}^{+\infty}e^{-iw(t-t')}e^{\pm i \frac{w}{c}|\vec r - \vec r'|}dw \\ = \frac{1}{2 \pi |\vec r - \vec r'|} \int_{-\infty}^{+\infty}e^{- iw [(t-t')\mp\frac{|\vec r - \vec r'|}{c}]}dw \\ \Rightarrow G^{\pm}(r,t-t')=\frac{\delta(t-t' \mp \frac{r}{c})}{r}G(r,t,r′,t′)=2π1?∫?∞+∞?re±ikr?e?iwtdw=2π1?∫?∞+∞?∣r?r′∣e±ik∣r?r′∣?e?iw(t?t′)dw=2π∣r?r′∣1?∫?∞+∞?e?iw(t?t′)e±icw?∣r?r′∣dw=2π∣r?r′∣1?∫?∞+∞?e?iw[(t?t′)?c∣r?r′∣?]dw?G±(r,t?t′)=rδ(t?t′?cr?)?
這個(gè)表達(dá)式說明兩點(diǎn):
另外需要說明的是時(shí)變Green函數(shù)中的?\mp?,如果是?-?,那么Green函數(shù)中的Dirac函數(shù)項(xiàng)為δ(t?t′?rc)\delta(t-t'-\frac{r}{c})δ(t?t′?cr?),這表示的是在觀察者觀測的時(shí)刻ttt之前的某個(gè)時(shí)刻t′t't′電磁場開始傳播,傳播到觀察者的位置需要的時(shí)間為r/cr/cr/c,這個(gè)解被稱為retarded solution;如果是+++,那么Dirac函數(shù)為δ(t?t′+rc)\delta(t-t'+\frac{r}{c})δ(t?t′+cr?),它可以理解為在觀察者觀測時(shí)刻ttt之后的某個(gè)時(shí)刻t′t't′開始傳播的電磁場回溯到觀察者的位置時(shí)間為r/cr/cr/c,也可以理解為觀察者從ttt時(shí)刻開始光速運(yùn)動(dòng)r/cr/cr/c時(shí)間后正好可以觀測到t′t't′時(shí)傳播到r?′\vec r'r′的電磁場,這個(gè)解叫advanced solution。雖然advanced solution的兩種解釋看上去都很奇怪,但它并沒有違背因果律,后續(xù)討論邊界條件的時(shí)候我們會(huì)討論advanced solution的物理含義。
對(duì)于retarded solution,根據(jù)Green函數(shù)的作用,無邊界條件時(shí),電磁場的potential可以表示為
Φ(?)(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′
總結(jié)
以上是生活随笔為你收集整理的UA PHYS515A 电磁理论IV 时变电磁场理论2 Helmholtz方程与含时的Green函数的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: Paper Review: Bayesi
- 下一篇: UA PHYS515A 电磁理论IV 时