• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    運(yùn)動(dòng)介質(zhì)洛倫茲協(xié)變電磁理論

    2018-02-05 04:50:35王雯宇
    物理與工程 2018年2期
    關(guān)鍵詞:麥克斯韋電磁學(xué)本構(gòu)

    王雯宇 許 洋

    (北京工業(yè)大學(xué)應(yīng)用數(shù)理學(xué)院,北京 100124)

    1 運(yùn)動(dòng)介質(zhì)電磁學(xué)理論的協(xié)變性問題

    物理學(xué)專業(yè)學(xué)習(xí)的宏觀經(jīng)典電磁理論通常以電磁學(xué)和電動(dòng)力學(xué)兩門課程的形式來完成。電磁學(xué)課程從基本的電磁現(xiàn)象開始,逐步講解電磁場(chǎng)的各種場(chǎng)屬性,所有的屬性總結(jié)為麥克斯韋方程組作為結(jié)束;而電動(dòng)力學(xué)則開始于麥克斯韋方程組,利用各種數(shù)學(xué)理論研究電磁場(chǎng)相關(guān)應(yīng)用,包括電磁波,電磁輻射,狹義相對(duì)論等。

    相比于電磁學(xué),從某種意義上講,電動(dòng)力學(xué)唯一的“新”物理在于狹義相對(duì)論。因?yàn)殡姶挪ǖ认嚓P(guān)理論實(shí)際上只是麥克斯方程組的演繹,而真正概念上的革新則在于相對(duì)性原理、光速不變?cè)淼囊牒拖鄬?duì)論時(shí)空觀以及力學(xué)理論的修正等。相對(duì)性原理指的是所有物理規(guī)律在不同的參考系下的形式是一樣的,不存在一個(gè)特殊的參考系。因此,平直時(shí)空中,原則上所有理論都可以寫成洛倫茲協(xié)變即狹義相對(duì)論性理論。真空中的麥克斯韋方程組自然可以寫為洛倫茲協(xié)變的四維形式,但是對(duì)于介質(zhì)存在時(shí)麥克斯韋方程組則不盡然。雖然在部分教科書中有所提及,比如著名電動(dòng)力學(xué)教材杰克遜的《經(jīng)典電動(dòng)力學(xué)》[1]中就簡單說明了介質(zhì)麥克斯韋方程組洛倫茲變換的處理辦法;但是書中并沒有給出方程組的協(xié)變形式。具體地講就是,為了描述電磁場(chǎng)對(duì)介質(zhì)中微觀粒子的極化與磁化,電磁學(xué)中引入了與介質(zhì)性質(zhì)相關(guān)的介電常數(shù)與磁導(dǎo)率兩個(gè)物理量。宏觀電磁場(chǎng)論又引入了電位移矢量D、磁場(chǎng)強(qiáng)度H、極化強(qiáng)度P與磁化強(qiáng)度M等4個(gè)場(chǎng)量。它們和電場(chǎng)強(qiáng)度E與磁感應(yīng)強(qiáng)度B由本構(gòu)關(guān)系相聯(lián)系。但是通常的教科書本構(gòu)關(guān)系一般表示為分量形式,并不是明顯洛倫茲協(xié)變的。雖然實(shí)際應(yīng)用中介質(zhì)電磁相關(guān)理論應(yīng)用已經(jīng)足夠廣泛,但是從理論的角度看,寫下協(xié)變的介質(zhì)麥克斯韋方程組是理解狹義相對(duì)論理論所必須的。另外,介質(zhì)中的光速也不再是真空極限速度c0,同時(shí)介質(zhì)中電磁波波動(dòng)方程中的速度也因此不是明顯洛倫茲協(xié)變的。那如何將此時(shí)的介質(zhì)麥克斯韋方程組以及波動(dòng)方程寫為協(xié)變形式,還有介質(zhì)運(yùn)動(dòng)時(shí)的電磁波多普勒效應(yīng)具體形式等,都是非常有趣且有意義的問題。

    綜上,我們給出運(yùn)動(dòng)介質(zhì)中協(xié)變形式的電磁學(xué)方程,包括麥克斯韋方程組和波動(dòng)方程,就是本文的目的。本文分以下幾部分論述:第1部分,給出運(yùn)動(dòng)介質(zhì)中電磁學(xué)現(xiàn)象的相關(guān)論述,并列舉歷史上的一些運(yùn)動(dòng)介質(zhì)中的電磁學(xué)模型;第2部分,嘗試給出運(yùn)動(dòng)介質(zhì)電磁學(xué)理論的洛倫茲協(xié)變形式,并研究介質(zhì)中電磁波波速疊加、多普勒頻移等狹義相對(duì)論現(xiàn)象。最后一節(jié)給出結(jié)論。

    2 運(yùn)動(dòng)介質(zhì)電磁學(xué)

    從現(xiàn)代物理學(xué)角度來看宏觀電磁現(xiàn)象,我們的認(rèn)知其實(shí)并沒有超出麥克斯韋多少。雖然麥克斯韋當(dāng)時(shí)不可能了解狹義相對(duì)論,但是他寫下的方程組卻自然滿足(形式上并不是明顯的)洛倫茲協(xié)變性,這歸功于麥克斯韋對(duì)于電磁場(chǎng)理論深刻的認(rèn)識(shí)。到19世紀(jì)中葉,物理學(xué)家已經(jīng)發(fā)現(xiàn)并總結(jié)電磁現(xiàn)象的3個(gè)最基本的實(shí)驗(yàn)定律: 庫侖定律(1785年),畢奧-薩伐爾定律(1820年),法拉第定律(1831—1845年)[2]。法拉第的“電力線”和“磁力線”等概念已發(fā)展成“電場(chǎng)”和“磁場(chǎng)”概念。1855年至1865年期間,麥克斯韋在全面地審視以上電磁學(xué)3個(gè)定律的基礎(chǔ)上,把數(shù)學(xué)分析方法引入電磁學(xué)理論中,應(yīng)用嚴(yán)謹(jǐn)?shù)臄?shù)學(xué)形式總結(jié)了前人的工作。他提出了位移電流假說和感生電場(chǎng)概念,將電磁場(chǎng)基本定律歸結(jié)為4個(gè)偏微分方程。至此完整地描述宏觀電磁現(xiàn)象的麥克斯韋方程組就誕生了。真空麥克斯韋方程組為[3]

    相關(guān)符號(hào)是大多數(shù)教科書上通用的,所以就不一一說明了。這是有電荷和電流、沒有磁荷的麥克斯韋方程組。電荷密度與電流密度應(yīng)當(dāng)滿足電荷守恒定律,即

    (5)

    需要補(bǔ)充說明的是,宏觀電磁學(xué)還應(yīng)該有另外一個(gè)基本公式,即廣義洛倫茲力密度

    f=ρE+J×B

    (6)

    洛倫茲力實(shí)際上是給出了電場(chǎng)和磁場(chǎng)的定義,而麥克斯韋方程組描述的則是電場(chǎng)和磁場(chǎng)的動(dòng)力學(xué)規(guī)律。

    麥克斯韋方程組深刻地揭示了源與場(chǎng)之間的關(guān)系:電荷在真空中激發(fā)了電場(chǎng),電流則激發(fā)了磁場(chǎng),變化的電場(chǎng)與磁場(chǎng)又互相激發(fā),形成電磁波向外傳播。因此麥克斯韋電磁理論與牛頓力學(xué)之間的關(guān)鍵差別在于局域性概念的引入,即它是一個(gè)場(chǎng)的理論。隨后愛因斯坦建立的狹義相對(duì)論與之緊密相關(guān)。狹義相對(duì)論要求場(chǎng)具有相應(yīng)的四維時(shí)空對(duì)稱性,其中包括時(shí)空平移對(duì)稱性與洛倫茲變換對(duì)稱性,以及空間反射對(duì)稱性(P宇稱)與時(shí)間反演對(duì)稱性(T宇稱)等。電場(chǎng)和磁場(chǎng)滿足的對(duì)稱性對(duì)理論有嚴(yán)格的限制。例如,電場(chǎng)強(qiáng)度E是一個(gè)矢量場(chǎng),空間反射變號(hào)。根據(jù)公式(1)可知磁感應(yīng)強(qiáng)度B是一個(gè)軸矢量場(chǎng)。類似地,電流密度J是一個(gè)矢量場(chǎng),ρ是一個(gè)標(biāo)量場(chǎng)等。

    能明顯體現(xiàn)理論所包含的對(duì)稱性的方法是使用電磁場(chǎng)的分析力學(xué)拉格朗日密度。經(jīng)典電磁理論的拉格朗日密度

    (7)

    具有相應(yīng)的龐加萊時(shí)空對(duì)稱性。其中,Aμ是電磁勢(shì)矢量,F(xiàn)μ ν=?μAν-?νAμ是電磁場(chǎng)張量,具體形式為

    (8)

    Jμ是電流四矢量。注意這里我們采用的是自然單位制取真空中的光速為c0=1,(下同)這里使用的度規(guī)是

    gμ ν=diag(1,-1,-1,-1)

    在張量的表述下,麥克斯韋方程組變?yōu)?/p>

    將Aμ作為基本變量,上式第一個(gè)方程是歐拉-拉格朗日運(yùn)動(dòng)方程的直接結(jié)果; 而第二個(gè)方程是一個(gè)恒等式,方程右邊等于零意味著不存在磁荷。這兩個(gè)方程的分量形式就是麥克斯韋方程組,此即宏觀電磁場(chǎng)分析力學(xué)形式。既然真空中的麥克斯韋方程組可以寫為滿足狹義相對(duì)論要求的四維形式,我們不禁要問,介質(zhì)存在時(shí)的麥克斯韋方程組是否也是一個(gè)狹義相對(duì)論理論?此時(shí),由于電磁場(chǎng)作用下介質(zhì)微觀結(jié)構(gòu)改變,電磁理論變得復(fù)雜多了,下面我們就來專門進(jìn)行討論。

    外加恒定電磁場(chǎng)的靜止介質(zhì)中,微觀層次上發(fā)生著微妙的電磁現(xiàn)象。一方面,電場(chǎng)會(huì)使得介質(zhì)中原本對(duì)外呈電中性的正負(fù)電荷發(fā)生微小位移,定向排列形成電偶極矩,在介質(zhì)的內(nèi)部或者表面形成束縛電荷;另一方面,磁場(chǎng)會(huì)使得介質(zhì)中分子或原子內(nèi)由于電子或正電荷運(yùn)動(dòng)形成的不規(guī)則磁矩發(fā)生定向排列,微觀上形成一個(gè)個(gè)定向的電流環(huán),宏觀上會(huì)造成磁化電流。為了描述介質(zhì)微觀的極化與磁化,電磁學(xué)中引入極化強(qiáng)度P與磁化強(qiáng)度M兩個(gè)物理量。一般來說,它們與電場(chǎng)強(qiáng)度E和磁感應(yīng)強(qiáng)度B之間的關(guān)系為

    此時(shí),電磁學(xué)方程中的兩個(gè)變量E、B之外又增加了4個(gè),即D、H、P、M。若介質(zhì)處于變化的電磁場(chǎng)中,那么由于變化的電場(chǎng)與磁場(chǎng)會(huì)造成變化的束縛電荷與磁化電流,而變化的束縛電荷與磁化電流又將激發(fā)變化的電磁場(chǎng),與之前的電磁場(chǎng)疊加,這就使得情況變得非常復(fù)雜。如何寫下一組運(yùn)動(dòng)介質(zhì)中麥克斯韋方程組一直存在爭(zhēng)議。歷史上存在著一系列運(yùn)動(dòng)介質(zhì)中電磁學(xué)模型,各個(gè)模型在實(shí)驗(yàn)上的觀測(cè)效應(yīng)差別很小。它們之間的主要區(qū)別在于對(duì)待束縛電荷與極化電流的處理方式不同,即外加電磁場(chǎng)導(dǎo)致的電荷與電流在不同模型中使用的電磁學(xué)變量是不同的。

    現(xiàn)在使用最廣泛的介質(zhì)中的麥克斯韋方程組是閔可夫斯基(Minkowski)公式,它也是最早出現(xiàn)的介質(zhì)中的電磁學(xué)方程。這組方程形式優(yōu)美,對(duì)微觀現(xiàn)象的描述清晰明確,加上與實(shí)驗(yàn)符合得也很好,因此不管是在教科書中還是工程應(yīng)用上都被廣為采用。在閔可夫斯基公式中,束縛電荷ρp被認(rèn)為是由介質(zhì)的極化造成,而額外的電流則包含磁化電流JM,極化電流密度JP,它們與電磁學(xué)變量之間的關(guān)系為

    這種介質(zhì)極化與磁化描述在教科書中有詳細(xì)的介紹,在此就不再解釋了。因此由此導(dǎo)致的麥克斯韋方程組形式為

    ρf與Jf分別為自由電荷與自由電流。閔可夫斯基公式選擇了兩組電磁學(xué)變量,即(D,H),(E,B),分別描述了介質(zhì)的微觀規(guī)律與電磁場(chǎng)的性質(zhì)。各向同性介質(zhì)中,它們之間由本構(gòu)關(guān)系

    相聯(lián)系。注意在下文中,我們討論的都是各向同性的介質(zhì)。杰克遜的《經(jīng)典電動(dòng)力學(xué)》書中將閔可夫斯基公式中的麥克斯韋方程組寫為了四維協(xié)變形式

    其中

    (24)

    是Fμ ν的對(duì)偶形式。(εμ ν α β是四階反對(duì)稱張量。)方程(23)說明介質(zhì)存在時(shí)恒等式依然成立。二階全反對(duì)稱張量Gμ ν的定義為將電磁張量Fμ ν中的E替換為D,B替換為H[1],即

    (25)

    此時(shí)介質(zhì)中的麥克斯韋方程組雖然寫成了協(xié)變形式,但是仍然不能認(rèn)為這是介質(zhì)電磁學(xué)的協(xié)變形式理論。原因有二:第一,這一公式是通過與真空中的麥克斯韋方程組類比而得到的,不能由拉格朗日密度的變分得到。在真空電磁場(chǎng)的拉格朗日密度中,基本變量是Aμ。而介質(zhì)存在時(shí),基本變量是Gμ ν,因此很難構(gòu)造相應(yīng)的拉格朗日量。第二,更為重要的是,本構(gòu)關(guān)系(20)、(21)仍然是分量形式,只是在介質(zhì)靜止的參考系中成立,不是洛倫茲協(xié)變的。因此即使運(yùn)動(dòng)方程協(xié)變,由于不清楚本構(gòu)關(guān)系的協(xié)變性質(zhì),理論也不是明顯協(xié)變的。第3節(jié)將對(duì)以上兩點(diǎn)再次討論并把方程修改為協(xié)變形式。在此之前,先來了解一下,其他的介質(zhì)電磁學(xué)模型。

    另外一個(gè)影響較大的介質(zhì)電磁學(xué)模型,是在20世紀(jì)50年代末,朱蘭成(Lan-Jen Chu*朱蘭成(1913—1973),江蘇淮陰人,美籍,國際著名電磁波專家,麻省理工學(xué)院教授,美國科學(xué)院院士。)提出的一組運(yùn)動(dòng)介質(zhì)中的介質(zhì)電磁學(xué)理論公式[4]。這組公式與閔可夫斯基公式不完全相同。在基本變量上,閔可夫斯基公式選擇了E,B,D,H,而朱公式選擇了E,H,M與P。在對(duì)待磁化的觀點(diǎn)上,朱公式認(rèn)為磁化的原因是磁偶極矩所造成的,即μ0M=q*n*d*,q*是磁荷量,n*是磁偶極矩密度,d*是由負(fù)磁荷指向正磁荷的長度矢量。若定義磁標(biāo)量φm,滿足

    H=-φm

    有關(guān)系

    ·H=-2φm=·M

    (26)

    因此

    2φm=·M

    (27)

    為了與靜電勢(shì)滿足的泊松方程

    形式保持一致,朱引入了磁荷密度

    ρm=-μ0·M

    (28)

    因此

    (29)

    朱重新分析了介質(zhì)中微觀極化與磁化的規(guī)律,認(rèn)為束縛電荷密度形式為

    極化電流密度形式為

    (30)

    磁荷密度為

    磁化磁流密度形式為

    (31)

    其中,Jf是自由電流,ρf是自由電荷。相比閔可夫斯基公式,朱公式稍顯復(fù)雜。利用本構(gòu)關(guān)系(20)、(21),上式最后兩式其實(shí)與閔可夫斯基公式相應(yīng)公式完全相同,而前兩式卻明顯不同,不能由場(chǎng)的重新定義而被吸收掉。前兩式的方程右邊表示與閔可夫斯基公式描述不同的極化與磁化方式。取=0,朱公式將簡化為閔可夫斯基公式。同時(shí)應(yīng)該注意到,朱電磁理論是一個(gè)非相對(duì)論理論。朱公式中出現(xiàn)的極化電流項(xiàng)×(P×)和磁化磁流項(xiàng)×(μ0M×)是這個(gè)模型中獨(dú)特的兩項(xiàng)。與此類似,許多其他的電磁學(xué)模型(例如Boffi公式,Amperian公式[4]等,本文不再敘述,有興趣讀者可以查閱相關(guān)文獻(xiàn)。)也引入含有速度的項(xiàng)。這是非常有意思的。下面將利用空間平均方法來說明一下怎樣理解不同的介質(zhì)極化和磁化模型。

    朱公式是在微觀極化與磁化的基礎(chǔ)上考慮介質(zhì)的運(yùn)動(dòng)效應(yīng)得到的,因此包含著介質(zhì)的運(yùn)動(dòng)速度。通常,由于偶極子對(duì)外呈電中性,其定向運(yùn)動(dòng)似乎不會(huì)帶來電流或磁流。其實(shí)不是這樣的,下面說明到底怎么回事。由于迄今為止實(shí)驗(yàn)上沒有觀測(cè)到真正的磁單極,下文將不再關(guān)注朱模型的磁荷部分,而僅討論介質(zhì)運(yùn)動(dòng)產(chǎn)生的極化電流,這已經(jīng)足以說明問題。宏觀電動(dòng)力學(xué)不是考察個(gè)別粒子產(chǎn)生的微觀電磁場(chǎng),而是關(guān)注于它們的宏觀平均值。微觀的極化與磁化因模型而異,為了研究介質(zhì)在外加電磁場(chǎng)下可能出現(xiàn)哪些額外的電荷項(xiàng)與電流項(xiàng),一種相對(duì)嚴(yán)格的處理方法是對(duì)微觀電荷與電流取空間平均(spatial average)。對(duì)函數(shù)F(x,t)取空間平均定義為

    (36)

    在這里,F(xiàn)(x,t)代表著微觀電磁物理量??臻g平均操作是對(duì)x′進(jìn)行的,因此它與時(shí)間導(dǎo)數(shù)算符對(duì)易,即

    (37)

    其中權(quán)函數(shù)f(x)被假設(shè)為連續(xù)光滑,量綱為L-3的非負(fù)實(shí)函數(shù),在全空間歸一化為1。在后面計(jì)算中,只需要用到它的泰勒展開式,即

    (38)

    其中i,j=1,2,3,重復(fù)的上下指標(biāo)代表求和??紤]介質(zhì)中由電子、離子、分子組成的微觀系統(tǒng),視這些帶電粒子為點(diǎn)電荷。這些點(diǎn)電荷以某種空間構(gòu)型組成宏觀介質(zhì),因?yàn)槲⒂^上就沒有介質(zhì)存在,利用真空的麥克斯韋方程組即可。因此在這個(gè)層面上的電磁學(xué)規(guī)律由下列方程表述

    其中,小寫的e,b代表微觀的電場(chǎng)和磁場(chǎng),η,j代表微觀的電荷密度與電流密度。相應(yīng)微觀量與宏觀物理量E,B,ρ,J的關(guān)系定義為

    注意,這里的處理辦法仍然是經(jīng)典物理的方法,即我們認(rèn)為物體可以分割成非常小的微元,微元的物理是經(jīng)典的,它們的宏觀效應(yīng)可以通過大量微元疊加或者積分的方式獲得?,F(xiàn)代物理,在微觀的尺度上,實(shí)際上應(yīng)當(dāng)用量子力學(xué)描述電子與核子的狀態(tài)。利用量子力學(xué)中算符的期望值代替經(jīng)典量時(shí),二者近乎相等[1]??臻g平均方法的詳細(xì)內(nèi)容,讀者可參閱Russakoff[5];而Ian Murdoch更多的討論了權(quán)函數(shù)的選擇[6];對(duì)于極化與磁化的相對(duì)論修正,以及統(tǒng)計(jì)力學(xué)方法的引入,de Groot與Suttorp進(jìn)行了全面的討論[7]; 利用量子力學(xué)的方法處理微觀電磁現(xiàn)象的內(nèi)容,可參閱Mahan[8],Adler[9]。

    圖1 介質(zhì)微觀電磁結(jié)構(gòu)的經(jīng)典模型

    進(jìn)行空間平均計(jì)算前,最好通過定義坐標(biāo)系將物理圖像闡釋清楚。如圖1令xn為坐標(biāo)原點(diǎn)O指向第n個(gè)分子質(zhì)心的位置矢量,xl為O點(diǎn)指向第l個(gè)帶電粒子ql的位置矢量。這些帶電粒子包括傳導(dǎo)電流的自由電子與分子中形成多極矩的電子和核子。xln為第n個(gè)分子中心指向分子中第l個(gè)電子或核子的位置矢量

    xln=xl-xn

    對(duì)時(shí)間求導(dǎo)有

    (45)

    它是權(quán)函數(shù)的零階近似,代表著將分子看成點(diǎn)電荷的情況。分子偶極矩定義為

    (46)

    這與通常教科書中定義的偶極矩p=qd是一致的,d指的是兩個(gè)帶等量異種電荷q的點(diǎn)電荷分開的一小段微小距離。如果分子帶電重心xn相互重合,分子偶極矩將為零。四極矩定義為

    (47)

    作為權(quán)函數(shù)的二階近似,電四極矩代表著分子內(nèi)電荷偏離球?qū)ΨQ分布的程度。微觀電荷分為傳導(dǎo)電流的自由電荷與束縛在分子內(nèi)部的電子與核子,可以寫為

    (48)

    我們對(duì)自由電荷與束縛電荷進(jìn)行空間平均計(jì)算,利用權(quán)函數(shù)的泰勒展開公式,對(duì)所有分子n取和即可得到宏觀電荷表達(dá)式。宏觀電荷密度最終的表達(dá)式為

    (49)

    其中,ρ自由(x,t)是宏觀自由電荷密度

    (50)

    P(x,t)是宏觀極化矢量

    (51)

    以及宏觀電四極矩張量

    (52)

    注意,這里應(yīng)該區(qū)分牛頓力學(xué)與場(chǎng)論的不同之處。位置矢量xl,xln、速度矢量l,ln不是場(chǎng),不隨坐標(biāo)的改變而發(fā)生改變。而宏觀極化矢量P是一個(gè)場(chǎng),它在不同的空間點(diǎn)取不同的值。原因在于引入的權(quán)函數(shù)f(x-xn)是坐標(biāo)的函數(shù),它攜帶著場(chǎng)的信息。這一點(diǎn)是牛頓力學(xué)和場(chǎng)論區(qū)別的微妙之處,在以下的推導(dǎo)過程中非常重要,即所有對(duì)空間的微分運(yùn)算都是對(duì)于權(quán)函數(shù)來進(jìn)行,而此時(shí)點(diǎn)電荷的位置、速度矢量相當(dāng)于一個(gè)常數(shù)。

    通過空間平均的方法,可以視極化現(xiàn)象為分子內(nèi)部由于電荷分布不均勻而產(chǎn)生的多極矩效應(yīng),并且最后的結(jié)果并不依賴權(quán)函數(shù)的選擇。xln是分子半徑的尺度,大約為10-10m,因此電四極矩含有一個(gè)10-20m-2的壓低因子,讀者或許因此認(rèn)為忽略電四極矩的貢獻(xiàn)可以的。但是下面的計(jì)算會(huì)顯示,微觀電四極矩貢獻(xiàn)了一部分的分子磁矩。因此在這里我們選擇將權(quán)函數(shù)泰勒展開到二階。將式(49)代入式(41)得到

    ·(ε0E+P-Q+…)=ρ自由

    (53)

    其中,矢量Q分量為

    上式意味著可以定義宏觀電位移矢量

    D=ε0E+P-Q

    (54)

    即它的散度等于介質(zhì)中自由電荷密度。注意,到目前為止的電位移矢量中沒有介質(zhì)的速度,根據(jù)以上公式看不出來介質(zhì)的運(yùn)動(dòng)會(huì)不會(huì)造成額外的電荷項(xiàng)。

    電荷移動(dòng)時(shí)便會(huì)形成電流。在電動(dòng)力學(xué)教科書中已經(jīng)知道,介質(zhì)中的電流項(xiàng)來自于自由電流密度、極化電流密度與磁化電流密度,而朱公式中出現(xiàn)了額外的電流項(xiàng):×(P×),顯然這一項(xiàng)是由于介質(zhì)的運(yùn)動(dòng)速度造成的。利用空間平均的方法求出這一電流項(xiàng)就是接下來的任務(wù)。為此,假設(shè)微觀電流密度分為自由電流密度與分子電流密度

    (55)

    (56)

    很難看出上式方程右邊是對(duì)應(yīng)怎樣的物理量。下面逐項(xiàng)處理來得到朱公式(33)中的所有項(xiàng)以及其他更多的項(xiàng)。

    (57)

    推導(dǎo)第二步用到了關(guān)系

    (58)

    推導(dǎo)第二步用到了關(guān)系

    為了得到極化電流密度的表達(dá)式,對(duì)泰勒展開式(56)加上下面的式(59),之后再減去式(59)

    (59)

    上式第一項(xiàng)與式(56)第三項(xiàng)一并貢獻(xiàn)了磁化電流密度JM,即

    (60)

    (61)

    公式(59)最后一項(xiàng)與泰勒展開式(56)最后一項(xiàng)合并對(duì)n求和后的領(lǐng)頭階為

    -×(Q×n)

    (62)

    過程與上面推導(dǎo)過程類似,這里就省略了。

    再減去的項(xiàng)(59)與公式(57)第一項(xiàng)?P/?t組成了

    (63)

    宏觀自由電流密度定義為

    (64)

    M為宏觀磁化強(qiáng)度矢量

    (65)

    (66)

    依據(jù)上式,可以重新定義

    (67)

    這說明電偶極矩與電四極矩在介質(zhì)運(yùn)動(dòng)時(shí)貢獻(xiàn)了一個(gè)等效磁矩,它們將以此貢獻(xiàn)一個(gè)分子電流環(huán)。注意這一項(xiàng)并不違反電流守恒定律。極化電荷密度ρp=-·P,根據(jù)電流守恒定律式(5),它對(duì)時(shí)間的導(dǎo)數(shù)將導(dǎo)致極化電流密度,即

    (68)

    因?yàn)槿我馐噶繄?chǎng)的旋度是無源的

    任意電流守恒方程都存在一個(gè)額外的自由度,電流可以重新定義為

    (69)

    其中A為任意矢量。使用空間平均的方法

    A=(P-Q)×

    同時(shí),由于(P-Q)×是一個(gè)軸矢量,×[(P-Q)×]就是一個(gè)矢量,這也不違反空間反射對(duì)稱性的要求。如果極化電流Jp在xy平面內(nèi)向x軸正方向流動(dòng),那么電流×[(P-Q)×]代表著在xz平面內(nèi)的閉合電流環(huán),不影響流守恒。

    以上我們依照空間平均的方法分析極化與磁化的一般規(guī)律,發(fā)現(xiàn)運(yùn)動(dòng)介質(zhì)的極化與磁化可以認(rèn)為是由電荷與電流的多極矩造成的。介質(zhì)運(yùn)動(dòng)效應(yīng)對(duì)麥克斯韋方程組貢獻(xiàn)了額外的電流項(xiàng)×[(P-Q)×]。綜合分析,近似至二階,運(yùn)動(dòng)介質(zhì)中的麥克斯韋方程組為

    其中本構(gòu)關(guān)系為

    應(yīng)該注意到以上理論相比與朱公式(32)~(35)多了四極矩的貢獻(xiàn)。其實(shí)可以類比上文的處理辦法給出更高階的修正,其形式一般都可以化為某個(gè)多出來的矢量叉乘介質(zhì)速度的形式。這是因?yàn)楹暧^電荷密度的高階修正正是權(quán)函數(shù)的各階泰勒展開式之和,即

    (76)

    其中a1,a2,…,ak=1,2,…,k,Ra1a2…ak就是通常定義的宏觀電荷的2k極矩張量:

    (77)

    式(76)包含著權(quán)函數(shù)零頭階(分子電荷)、一階(偶極矩)、二階(四極矩)…第k階(2k極矩)以及更高階的貢獻(xiàn)。類似地,對(duì)于每一階我們都可以利用電流守恒,公式(58)以及上一階計(jì)算留下來的項(xiàng)湊出一個(gè)矢量Rk貢獻(xiàn)的電流密度項(xiàng):

    (-1)k+1×(Rk×)

    (78)

    代入方程(42),便可重新定義H。最后可以得到更高階修正的本構(gòu)關(guān)系為

    與閔可夫斯基公式中的D,H不同,這里的D′,H′包含著介質(zhì)運(yùn)動(dòng)的貢獻(xiàn)。這樣運(yùn)動(dòng)介質(zhì)中麥克斯韋方程組的形式與閔可夫斯基公式就相同了。而實(shí)際兩種模型D,H與D′,H′不相等造成的觀測(cè)效應(yīng)非常小(偶極矩、四極矩效應(yīng)等)。這也意味著,介質(zhì)運(yùn)動(dòng)產(chǎn)生的額外效果可以統(tǒng)一吸收到本構(gòu)關(guān)系中,麥克斯韋方程組形式保持不變。而將D′,H′組成協(xié)變張量Gμ ν就可以在做不同參考系變換時(shí)處理介質(zhì)效應(yīng)。如果要求理論是相對(duì)論理論的話,那本構(gòu)關(guān)系的參考系變換其實(shí)就成了問題的關(guān)鍵。在下一節(jié)中,將據(jù)此給出運(yùn)動(dòng)介質(zhì)中協(xié)變的本構(gòu)關(guān)系。

    3 介質(zhì)電磁理論的協(xié)變形式

    上節(jié)主要介紹了兩類模型,一類是閔可夫斯基公式,一類是朱公式,二者的不同點(diǎn)在于對(duì)待束縛電荷與磁化電流的方式有所不同。正如引文所說,狹義相對(duì)論對(duì)于電動(dòng)力學(xué)來說至關(guān)重要,那么介質(zhì)中的麥克斯韋方程組的相對(duì)論形式是什么樣的?這是個(gè)非常有趣的問題。從理論的觀點(diǎn)看,上節(jié)運(yùn)動(dòng)介質(zhì)中的電磁學(xué)方程存在的問題有:

    (1) 協(xié)變的麥克斯韋方程組。雖然閔可夫斯基公式可以總結(jié)為Gμ ν,Fμ ν四維協(xié)變的張量方程,但是其拉格朗日密度等很難處理。Gμ ν,Fμ ν只是分量形式清晰,Gμ ν與基本變量Aμ之間關(guān)系如何是不清楚的。

    (2) 本構(gòu)關(guān)系的變換性質(zhì)不明確。式(20)、(21)是在靜止系中得出的,它的協(xié)變性質(zhì)不明確。介質(zhì)非相對(duì)論運(yùn)動(dòng)造成的極化與磁化規(guī)律導(dǎo)致本構(gòu)關(guān)系的修正公式(75),如上節(jié)所講,它的協(xié)變性質(zhì)是問題的關(guān)鍵。

    (3) 如何檢驗(yàn)運(yùn)動(dòng)介質(zhì)中的電磁學(xué)理論的正確性?一個(gè)明顯的檢驗(yàn)方式即是運(yùn)動(dòng)介質(zhì)中的光速。介質(zhì)中的光速需要通過計(jì)算波動(dòng)方程來得到,而運(yùn)動(dòng)介質(zhì)中的光速必然是介質(zhì)中的光速與介質(zhì)運(yùn)動(dòng)速度的相對(duì)論相加,波動(dòng)方程也應(yīng)該是協(xié)變的,但是現(xiàn)在的介質(zhì)中的波動(dòng)方程是不能寫成明顯協(xié)變的。

    下面討論介質(zhì)中的波動(dòng)方程:利用矢量公式

    即可得到無源波動(dòng)方程

    (81)

    利用洛倫茲不變的微分算符

    可以將波動(dòng)方程進(jìn)一步寫為

    (82)

    εμ≠1,上式左邊括號(hào)第二項(xiàng)時(shí)間的二次偏導(dǎo)的系數(shù)不為零。因?yàn)橛新懵对谕獾穆鍌惼澴儞Q改變算符,波動(dòng)方程(82)僅僅在靜止參考系成立,它的協(xié)變性質(zhì)也不明確。問題的關(guān)鍵就是本構(gòu)關(guān)系一個(gè)分量方程,它僅在靜止系成立,因而波動(dòng)方程僅在靜止系成立。若希望得到一個(gè)協(xié)變性質(zhì)明顯的介質(zhì)中的波動(dòng)方程,必須將閔可夫斯基公式中的本構(gòu)關(guān)系改寫為四維協(xié)變的形式。

    其實(shí),對(duì)于介質(zhì)電磁理論的不同參考系的處理,已經(jīng)有一個(gè)辦法,那就是把所有介質(zhì)運(yùn)動(dòng)參考系的物理量變換到介質(zhì)靜止參考系,然后利用本構(gòu)關(guān)系式(20),式(21)處理即可。這其實(shí)就是杰克遜《經(jīng)典電動(dòng)力學(xué)》中的處理辦法。雖然這么做已經(jīng)很實(shí)用,但是從理論的角度講,我們還是應(yīng)該探索介質(zhì)電磁學(xué)理論的協(xié)變性質(zhì)。因?yàn)榻橘|(zhì)效應(yīng)其實(shí)就是大量微觀粒子真空麥克斯韋方程的宏觀效應(yīng),如果宏觀反而失去了相對(duì)論的形式,這是不自洽的。而當(dāng)我們明白了相對(duì)論的基本理論之后,這樣的拓展工作其實(shí)并不是很難,具體計(jì)算如下。

    首先說明得到協(xié)變形式的方法。本構(gòu)關(guān)系式(20)、(21)的四維形式,實(shí)際上是兩個(gè)張量Gμ ν與Fμ ν的之間關(guān)系。其實(shí)此時(shí)還有另外兩個(gè)協(xié)變張量,即介質(zhì)中的光速eμ與介質(zhì)的運(yùn)動(dòng)速度uμ

    eμ=γe(1,e),uμ=γ(1,β)

    (83)

    其中

    由于我們采用的是自然單位制,介質(zhì)中的光速為

    (84)

    介質(zhì)靜止時(shí),

    uμ=(1,0)

    (85)

    我們要做的工作就是利用這兩個(gè)張量以及一些協(xié)變常張量gμ ν,εμ ν α β湊出來本構(gòu)關(guān)系式(20)、(21)而已。為了從二階張量Fμ ν中得到E,即Fμ ν的第一列Fμ0,利用靜止介質(zhì)的四速度uμ=(1,0,0,0)就可以做到這一點(diǎn),即

    Di=εuαFiα

    (86)

    這就是本構(gòu)關(guān)系式(20)。Fμ ν的第一行也可以用相同的方法得到。改寫本構(gòu)關(guān)系的第二式需要一點(diǎn)技巧。首先先取Fμ ν的對(duì)偶形式Fμ ν,這樣B會(huì)出現(xiàn)在Fμ ν的第一行和第一列,然后用與公式(86)類似的方式得到

    (87)

    當(dāng)然,這仍然不是協(xié)變形式,為了得到完全的Gμ ν和Fμ ν的關(guān)系,需要用εμ ν α β再把Hi變換回到原來位置上。最后利用uμ就可以湊出協(xié)變本構(gòu)關(guān)系為

    (88)

    再加上介質(zhì)四維協(xié)變關(guān)系公式(為了清楚起見,這里重新寫一次式(22)、(23))

    就得到了非常漂亮的介質(zhì)電磁學(xué)洛倫茲協(xié)變理論。也就是說,介質(zhì)存在時(shí),需要處理兩個(gè)協(xié)變電磁學(xué)場(chǎng)量Gμ ν與Fμ ν。當(dāng)介質(zhì)運(yùn)動(dòng)時(shí),利用介質(zhì)運(yùn)動(dòng)的協(xié)變速度uμ和協(xié)變本構(gòu)關(guān)系(88)即可解決相關(guān)問題,而不需要專門把兩個(gè)場(chǎng)量換算到介質(zhì)靜止參考系去處理。

    到此,應(yīng)該說協(xié)變理論已經(jīng)基本完成,但是,如上文所說,大部分教材在處理介質(zhì)變換時(shí)采用的都是分量形式,下面用以上協(xié)變公式來重復(fù)一些杰克遜書中的相關(guān)公式,以驗(yàn)證該協(xié)變理論。為了簡單起見,假定介質(zhì)沿x軸正方向運(yùn)動(dòng),速度為β。這時(shí),如果處于介質(zhì)相對(duì)靜止參考系,本構(gòu)關(guān)系就是式(20)、(21)。此時(shí)洛倫茲平動(dòng)(boost)方向?yàn)閤軸負(fù)方向。利用Fμ ν的協(xié)變性質(zhì),從式(88)可以得到一般參考系中(D,H)與(E,B)的關(guān)系。寫成三維分量的形式

    這就是介質(zhì)運(yùn)動(dòng)參考系中處理本構(gòu)關(guān)系的分量公式,它可以代替之前不協(xié)變的本構(gòu)關(guān)系,與麥克斯韋方程組聯(lián)立起來構(gòu)成描述運(yùn)動(dòng)介質(zhì)中的電磁學(xué)理論。這樣,(D,H)與(E,B)的變換關(guān)系更加明確,有了協(xié)變的本構(gòu)關(guān)系,可以說Gμ ν與Fμ ν的地位是平等的?,F(xiàn)在修改過的閔可夫斯基公式應(yīng)當(dāng)寫為

    這一組方程就與杰克遜電動(dòng)力學(xué)中的形式一致了。這也是大多數(shù)教科書處理介質(zhì)電磁場(chǎng)變換的公式,需要強(qiáng)調(diào)一點(diǎn)的是上文協(xié)變關(guān)系式(88)~(90)形式更為簡捷明了。

    但是,正如上文所述,介質(zhì)電磁學(xué)理論只是真空電磁學(xué)的應(yīng)用。本構(gòu)關(guān)系也只是大量微觀帶電粒子的宏觀效應(yīng)。實(shí)際上,真正的電磁學(xué)變量只有電場(chǎng)強(qiáng)度E和磁感應(yīng)強(qiáng)度B,即Fμ ν。因此,原則上,我們應(yīng)該可以把本構(gòu)關(guān)系約掉,進(jìn)而去掉D和H,將方程組減少到真空麥克斯韋方程組一樣的4個(gè)方程,兩個(gè)描述電場(chǎng),兩個(gè)描述磁場(chǎng)。也就是說可以把介質(zhì)的運(yùn)動(dòng)效應(yīng)從本構(gòu)關(guān)系上轉(zhuǎn)移到麥克斯韋方程組上,而且也應(yīng)該是形式協(xié)變的。事實(shí)上,的確可以做到這一點(diǎn),下面就來實(shí)現(xiàn)這一點(diǎn)。首先將本構(gòu)關(guān)系式(91)、(92)代入麥克斯韋方程組中,應(yīng)用矢量公式

    方程(94)、(95)的前兩個(gè)方程可以被化簡為

    (99)

    (100)

    利用協(xié)變的本構(gòu)關(guān)系式(91)、(92)消除了麥克斯韋方程組中的兩個(gè)變量D,H,可得到的方程還很復(fù)雜,需要化簡。因?yàn)槁鍌惼澠絼?dòng)方向?yàn)閤負(fù)方向,關(guān)于B的矢量式中,第三項(xiàng)的中括號(hào)內(nèi),第一、二、三項(xiàng)只有x方向上有分量,第四、五項(xiàng)則沒有x方向上的分量,因此在式(99)中,

    將上式代入式(99)中得

    (101)

    上式揭示了運(yùn)動(dòng)介質(zhì)中束縛電荷的形式。通過式(99)與式(100)的互相代入,得到了運(yùn)動(dòng)介質(zhì)中負(fù)責(zé)電場(chǎng)的源,它不僅僅與自由電荷相關(guān),還與介質(zhì)的運(yùn)動(dòng)效應(yīng)相關(guān)。磁場(chǎng)的源也可同樣得到,為此,在式(99)兩邊乘以速度矢量β

    (102)

    應(yīng)用矢量公式

    (103)

    因?yàn)槭?102)僅在x方向上有分量,yz方向不需要考慮,將式(102)、(103)代入式(100)中,經(jīng)整理得到

    (104)

    上式揭示了運(yùn)動(dòng)介質(zhì)中的磁化電流形式相當(dāng)復(fù)雜,方程右邊的各項(xiàng)中,含有B的項(xiàng)都沒有x方向上的分量,平動(dòng)速度β的項(xiàng)只有x方向上的分量,因此不容易整理得到一個(gè)完整的公式。但是若把每一方向上的分量都列出,就可總結(jié)出它的一般表達(dá)式。

    x方向上,式(104)大括號(hào)內(nèi)只有E的貢獻(xiàn),然后寫出

    (105)

    整理一下,式(100)x方向上的分量方程可以寫為

    (106)

    而式(100)yz方向上的表達(dá)式比較簡單,可直接推出,結(jié)果為

    (107)

    以上兩式可以合并為

    (108)

    經(jīng)過復(fù)雜的推導(dǎo),將式(104)化簡與式(101)形式相似,得到了運(yùn)動(dòng)介質(zhì)中B與它的源之間的關(guān)系。經(jīng)過以上計(jì)算即可得到運(yùn)動(dòng)介質(zhì)中完全關(guān)于(E,B)的協(xié)變的電磁學(xué)方程

    這樣就將靜止介質(zhì)中的本構(gòu)關(guān)系式(20)、(21)寫成了協(xié)變形式(91)、(92),代入麥克斯韋方程組中,便消除了原來公式中的兩個(gè)電磁學(xué)變量D、H,同時(shí)原來的6個(gè)方程,現(xiàn)在只剩下了4個(gè)。在麥克斯韋方程組(109)~(112)中,只有兩個(gè)電磁學(xué)變量E、B,因此不需要額外的本構(gòu)關(guān)系??梢郧宄乜闯?,介質(zhì)的運(yùn)動(dòng)額外造成的束縛電荷與磁化電流表述在式(109)、(110)中。它們是關(guān)于β以及γ的項(xiàng),當(dāng)β=0時(shí),便回到了閔可夫斯基公式。

    方程組(109)~(112)不是明顯協(xié)變形式的,我們還可以將其寫為四維形式。使用方法和我們前面推導(dǎo)協(xié)變本構(gòu)關(guān)系的方法相同,用協(xié)變張量、矢量湊出介質(zhì)中的麥克斯韋方程組。放棄Gμ ν只使用協(xié)變張量Fμ ν、Jμ、介質(zhì)中的光速eμ與介質(zhì)的運(yùn)動(dòng)速度uμ等湊出協(xié)變形式。注意到介質(zhì)靜止時(shí),光速只能用到γe與e,而與e的方向無關(guān)。γe可以由

    γe=e·u

    (113)

    得到,而e則可以用γe算出

    (114)

    其他過程與協(xié)變本構(gòu)關(guān)系的湊法類似,這里就省略了。最終得到的運(yùn)動(dòng)介質(zhì)中麥克斯韋方程組的四維協(xié)變形式為

    這個(gè)僅關(guān)于Fμ ν的方程正是我們期望的,讀者可以驗(yàn)證該協(xié)變方程組與式(109)~(112)完全一致,且形式簡潔。在規(guī)范理論中,F(xiàn)μ ν就是規(guī)范場(chǎng)場(chǎng)強(qiáng)。這個(gè)運(yùn)動(dòng)方程的分析力學(xué)推導(dǎo)比較復(fù)雜,因?yàn)榇藭r(shí)有一個(gè)拉格朗日密度,又有一個(gè)統(tǒng)一的介質(zhì)運(yùn)動(dòng)速度uμ(構(gòu)成有限個(gè)自由度拉格朗日量),所以我們就不討論其分析力學(xué)形式了??梢钥闯觯?115)后邊除了第一項(xiàng)外,其他項(xiàng)一部分源自于自由電流,一部分源來自于電磁場(chǎng)本身,且二者都與介質(zhì)的運(yùn)動(dòng)相關(guān)。第2節(jié)中所討論的本構(gòu)關(guān)系(20)則沒有速度項(xiàng),這說明從微觀的角度分析介質(zhì)效應(yīng)也存在一定的局限性。原因在于第2節(jié)微觀方法是非相對(duì)論的,因此得到的宏觀本構(gòu)關(guān)系也是非相對(duì)論的,其協(xié)變性質(zhì)不清楚,用本節(jié)宏觀辦法,則可以彌補(bǔ)這個(gè)缺陷。

    處理完了麥克斯韋方程組的問題,接下來看協(xié)變波動(dòng)方程的處理辦法。由式(109)、(110)式再做叉乘運(yùn)算×(×E)或×(×B),即可得到運(yùn)動(dòng)介質(zhì)中無源的波動(dòng)方程為

    [?2+(εμ-1)(uμ?μ)2]Fμ ν=0

    (117)

    這時(shí),波動(dòng)方程寫成了協(xié)變的形式。當(dāng)β=0時(shí),式(117)回到式(82)。為了驗(yàn)證此波動(dòng)方程的正確性,計(jì)算運(yùn)動(dòng)介質(zhì)中的波速公式。波動(dòng)方程式(117)是一個(gè)二階偏微分方程,若洛倫茲平動(dòng)方向?yàn)閤方向,則方程化簡之后形式為

    (118)

    這是一個(gè)雙曲型方程。令

    則它的特征方程為

    adt2-cdtdx+bdx2=0

    兩族積分曲線為

    x+At=c1,x+Bt=c2

    其中,c1,c2為兩個(gè)常數(shù),A和B是方程

    (119)

    的兩個(gè)解??梢则?yàn)證,作坐標(biāo)變換

    (120)

    方程(118)即化簡為標(biāo)準(zhǔn)型

    (121)

    方程(118)化簡為一個(gè)非常簡單的波動(dòng)方程。我們可以將電磁場(chǎng)的表達(dá)式寫為

    Fμ ν(u,v)=C0exp[i(ωv-ku)]

    (122)

    其中ω=k;C0為任意常數(shù)。將u,v的表達(dá)式代入,即可得到(t,x)坐標(biāo)下的電磁場(chǎng)的表達(dá)式

    (123)

    可見運(yùn)動(dòng)介質(zhì)中的電磁波波動(dòng)的相位比真空中的情況復(fù)雜很多,當(dāng)ε=ε0,μ=μ0時(shí),電磁波的表達(dá)式就回到真空中較為簡單的形式;而當(dāng)β=0時(shí),電磁波的表達(dá)式回到了靜止介質(zhì)中形式。由式(123)得到的波速為

    (124)

    其中,e′是運(yùn)動(dòng)介質(zhì)中的光速。此式正是運(yùn)動(dòng)介質(zhì)中波速的相對(duì)論疊加形式。至此,我們寫出了協(xié)變的本構(gòu)關(guān)系,得到了運(yùn)動(dòng)介質(zhì)中關(guān)于(E,B)協(xié)變的方程組,進(jìn)而得到了協(xié)變的波動(dòng)方程,并檢驗(yàn)了運(yùn)動(dòng)介質(zhì)中的光速公式,是符合狹義相對(duì)論的。

    圖2 運(yùn)動(dòng)介質(zhì)多普勒效應(yīng)

    最后需要補(bǔ)充說明一下,介質(zhì)運(yùn)動(dòng)對(duì)電磁波多普勒效應(yīng)的影響,因?yàn)檫@是狹義相對(duì)論創(chuàng)立之前很多物理學(xué)家關(guān)注的問題。為了簡化,只討論一維運(yùn)動(dòng)。圖2是多普勒效應(yīng)的示意圖,由于波源(速度為vS,發(fā)射頻率為fS)和接收者(速度為vR,接收頻率為fR)相對(duì)介質(zhì)運(yùn)動(dòng)而造成了接收者接收到的波頻率改變的現(xiàn)象。在牛頓力學(xué)中的頻移規(guī)律為

    (125)

    其中u是介質(zhì)(靜止)中波速。注意在牛頓力學(xué)中,不必關(guān)注介質(zhì)運(yùn)動(dòng)效應(yīng),因?yàn)榻^對(duì)時(shí)空觀,存在絕對(duì)時(shí)間,把相關(guān)速度都轉(zhuǎn)換到介質(zhì)靜止參考系即可。但是在狹義相對(duì)論中,情況就不同了。首先必須選擇接收者靜止βR=0(自然單位制)參考系,因?yàn)閾Q參考系計(jì)算的頻率是沒有實(shí)際意義的。真空中波源運(yùn)動(dòng)時(shí),由于相對(duì)論時(shí)間膨脹效應(yīng),頻移公式應(yīng)該修改為

    (126)

    (127)

    注意,由于時(shí)間膨脹造成的波源發(fā)射頻率改變?nèi)匀皇?1/γS)fS,所以上式就不能化成公式(126)的形式了。如果介質(zhì)運(yùn)動(dòng),速度為β,頻移公式則為

    (128)

    可以看到,這里就必須區(qū)分介質(zhì)是否在運(yùn)動(dòng),這正是波速不是極限速度時(shí)的相對(duì)論結(jié)果。同時(shí)也可以看到,如果|βS|≥e,電磁波就出現(xiàn)了激波的現(xiàn)象。此時(shí)情況比較復(fù)雜,這里就不做進(jìn)一步討論了。

    4 結(jié)語

    本文主要討論了運(yùn)動(dòng)介質(zhì)電磁學(xué)的協(xié)變形式。雖然麥克斯韋方程組可以寫成協(xié)變形式,但由于本構(gòu)關(guān)系協(xié)變性質(zhì)不明顯,所以大多數(shù)教科書中的介質(zhì)電磁學(xué)理論不是明顯洛倫茲協(xié)變的。由于處理極化和磁化的方式不同,因而物理學(xué)家給出了不同的介質(zhì)麥克斯韋方程組。在第2節(jié)中我們列舉了其中兩個(gè)最廣泛使用的形式:閔可夫斯基公式和朱公式。然后利用空間平均的方法分析了運(yùn)動(dòng)介質(zhì)微觀極化與磁化的相關(guān)規(guī)律,并指出了不同模型的來源。由于介質(zhì)的運(yùn)動(dòng)效應(yīng)都可以被吸收到D,H的重新定義中,而麥克斯韋方程組形式保持不變,因此本構(gòu)關(guān)系的洛倫茲協(xié)變形式對(duì)于運(yùn)動(dòng)介質(zhì)中的電磁理論來說極其重要。因此在第3部分,將本構(gòu)關(guān)系寫成了洛倫茲協(xié)變形式,借此得到了關(guān)于E,B的介質(zhì)中協(xié)變的麥克斯韋方程組,重復(fù)出一些教科書中的介質(zhì)電磁學(xué)理論。我們采用的宏觀處理辦法可以彌補(bǔ)第2節(jié)微觀非相對(duì)論的方法的缺陷。然后將介質(zhì)中不明顯協(xié)變的波動(dòng)方程推廣為洛倫茲協(xié)變的波動(dòng)方程,并得到了相對(duì)論速度相加的運(yùn)動(dòng)介質(zhì)中的光速公式以及運(yùn)動(dòng)介質(zhì)中多普勒頻移公式??梢?,麥克斯韋方程組(109)~(112)清晰地描述了運(yùn)動(dòng)介質(zhì)中的電磁學(xué)現(xiàn)象。

    本文的意義在于明確了運(yùn)動(dòng)介質(zhì)中的麥克斯韋方程組以及本構(gòu)關(guān)系的協(xié)變形式,這對(duì)理解狹義相對(duì)論的電磁理論有一定的幫助。

    [1] JACKSON J. Classical electrodynamics[M].3rd.北京: 高等教育出版社,2004.

    [2] 郭奕玲,沈慧君.物理學(xué)史[M].北京:清華大學(xué)出版社,2005.

    [3] 郭碩鴻.電動(dòng)力學(xué)[M].3版.北京:高等教育出版社,2008.

    [4] PENFIELD P, HERMANN A. Electrodynamics of moving media[M]. Massachusetts: THE M.I.T.PRESS, 1967.

    [5] RUSSAKOFF G. A derivation of the macroscopic Maxwell equations[J]. American Journal of Physics, 1970, 38(10): 1188-1195.

    [6] MURDOCH AI. On spatially-averaged electrokinetics of point charges and Maxwell’s equation[J]. Journal of Elasticity, 2017: 1-35.

    [7] DEGROOT SR, LINDSAY RB. The Maxwell equations:non-relativistic and relativistic derivations from electron theory[J]. Physics Today, 1971, 24(7): 47-47.

    [8] MAHAN GD. Local-field corrections to Coulomb interactions[J]. Physical Review, 1967, 153(3): 983-988.

    [9] ADLER SL. Quantum theory of the dielectric constant in real solids[J]. Physical Review, 1962, 126(2): 413-420.

    猜你喜歡
    麥克斯韋電磁學(xué)本構(gòu)
    麥克斯韋: 掌控電磁之力的大魔法師
    Maxwell Loses a Tooth 麥克斯韋掉牙了
    雙麥克斯韋分布下極區(qū)中層塵埃粒子帶電研究
    基于設(shè)計(jì)的研究在探究式學(xué)習(xí)單開發(fā)中的應(yīng)用——以中國科技館電磁學(xué)學(xué)習(xí)單《奇妙的原電池》為例
    中國博物館(2019年2期)2019-12-07 05:41:24
    離心SC柱混凝土本構(gòu)模型比較研究
    抓住麥克斯韋妖的尾巴——重新定義能源
    能源(2018年10期)2018-12-08 08:02:48
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    高中電磁學(xué)學(xué)習(xí)常見問題及解決策略之我見
    一種新型超固結(jié)土三維本構(gòu)模型
    基于校園網(wǎng)的電磁學(xué)試題庫建設(shè)
    蜜桃亚洲精品一区二区三区| 怎么达到女性高潮| 脱女人内裤的视频| 国产极品精品免费视频能看的| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 久久久久九九精品影院| 尤物成人国产欧美一区二区三区| av天堂中文字幕网| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利免费观看在线| 色播亚洲综合网| 午夜激情欧美在线| 日韩欧美在线乱码| 99久久久亚洲精品蜜臀av| 美女黄网站色视频| 久久亚洲真实| 亚洲av不卡在线观看| 色5月婷婷丁香| 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品久久久久久毛片777| 欧美国产日韩亚洲一区| 少妇熟女aⅴ在线视频| ponron亚洲| 男女床上黄色一级片免费看| 久久午夜福利片| 女人被狂操c到高潮| 成人精品一区二区免费| 欧美+亚洲+日韩+国产| 少妇的逼好多水| 午夜免费激情av| 又爽又黄无遮挡网站| 18+在线观看网站| 亚洲第一区二区三区不卡| 人人妻,人人澡人人爽秒播| 一区二区三区激情视频| 免费一级毛片在线播放高清视频| 国产av一区在线观看免费| 国产精品野战在线观看| 亚洲,欧美,日韩| 一区二区三区激情视频| 91av网一区二区| 久久久成人免费电影| 国产三级黄色录像| 一个人免费在线观看电影| 在现免费观看毛片| 午夜精品久久久久久毛片777| 欧美另类亚洲清纯唯美| 国产av麻豆久久久久久久| h日本视频在线播放| 人人妻人人澡欧美一区二区| 欧美bdsm另类| 日韩欧美一区二区三区在线观看| 波多野结衣高清作品| 男女之事视频高清在线观看| 国产真实乱freesex| 国产单亲对白刺激| 窝窝影院91人妻| 国产老妇女一区| 久久国产乱子免费精品| 99在线人妻在线中文字幕| 91在线精品国自产拍蜜月| 成人性生交大片免费视频hd| 天堂影院成人在线观看| 少妇的逼水好多| 波野结衣二区三区在线| 国产午夜精品久久久久久一区二区三区 | 亚洲人成电影免费在线| 老司机午夜福利在线观看视频| 波多野结衣巨乳人妻| 身体一侧抽搐| 亚洲欧美日韩无卡精品| av国产免费在线观看| 久久人人爽人人爽人人片va | 国产在视频线在精品| 757午夜福利合集在线观看| 三级国产精品欧美在线观看| 少妇被粗大猛烈的视频| 久久精品久久久久久噜噜老黄 | 在线观看舔阴道视频| 久久香蕉精品热| 欧美激情久久久久久爽电影| 成人国产一区最新在线观看| 悠悠久久av| 日韩中字成人| 亚洲av.av天堂| 国产精品免费一区二区三区在线| 国产一区二区亚洲精品在线观看| 日日干狠狠操夜夜爽| 在现免费观看毛片| 别揉我奶头 嗯啊视频| 欧美日本视频| 三级男女做爰猛烈吃奶摸视频| 又粗又爽又猛毛片免费看| 婷婷六月久久综合丁香| 免费av毛片视频| 国产精品日韩av在线免费观看| 成人午夜高清在线视频| 免费人成视频x8x8入口观看| 国产美女午夜福利| 老女人水多毛片| 国产精品不卡视频一区二区 | 最近中文字幕高清免费大全6 | 少妇丰满av| 91av网一区二区| 日本 av在线| 最好的美女福利视频网| 久久精品国产自在天天线| 国产成人福利小说| 国产精品自产拍在线观看55亚洲| 最后的刺客免费高清国语| 男女床上黄色一级片免费看| 99在线人妻在线中文字幕| 99国产综合亚洲精品| 一区二区三区四区激情视频 | 给我免费播放毛片高清在线观看| 天堂av国产一区二区熟女人妻| 99久久久亚洲精品蜜臀av| 一个人免费在线观看的高清视频| 蜜桃久久精品国产亚洲av| 欧美区成人在线视频| www.999成人在线观看| 国产精品美女特级片免费视频播放器| 黄色配什么色好看| 高清日韩中文字幕在线| 一级黄片播放器| 欧美绝顶高潮抽搐喷水| 亚州av有码| 国产人妻一区二区三区在| 欧美高清成人免费视频www| 国产真实伦视频高清在线观看 | 亚洲欧美精品综合久久99| 熟妇人妻久久中文字幕3abv| 色吧在线观看| 夜夜夜夜夜久久久久| 亚洲一区二区三区不卡视频| avwww免费| 亚洲黑人精品在线| 91麻豆av在线| 在线观看免费视频日本深夜| 国产精品三级大全| 性欧美人与动物交配| 亚洲最大成人中文| 国产高清激情床上av| 久久天躁狠狠躁夜夜2o2o| av黄色大香蕉| 日韩免费av在线播放| or卡值多少钱| 日韩欧美在线二视频| 国产精品1区2区在线观看.| 免费看a级黄色片| 制服丝袜大香蕉在线| 91麻豆av在线| 午夜福利在线观看免费完整高清在 | 亚洲激情在线av| 午夜免费激情av| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| 性色av乱码一区二区三区2| 老司机福利观看| 精品福利观看| 久久精品影院6| 亚洲成人久久爱视频| netflix在线观看网站| 国产真实伦视频高清在线观看 | 国产亚洲av嫩草精品影院| 国内精品久久久久精免费| 国产精品综合久久久久久久免费| 一夜夜www| 欧美另类亚洲清纯唯美| 好男人电影高清在线观看| 亚洲av第一区精品v没综合| 亚洲不卡免费看| 精品一区二区三区视频在线观看免费| 九九热线精品视视频播放| 国产69精品久久久久777片| 国产精品乱码一区二三区的特点| 国产主播在线观看一区二区| www.www免费av| 淫妇啪啪啪对白视频| 午夜福利高清视频| 日韩欧美免费精品| av在线老鸭窝| 婷婷精品国产亚洲av在线| 免费大片18禁| 又黄又爽又刺激的免费视频.| 麻豆成人av在线观看| 天堂动漫精品| 成熟少妇高潮喷水视频| 日日摸夜夜添夜夜添av毛片 | 91九色精品人成在线观看| 色在线成人网| 最后的刺客免费高清国语| 国产探花在线观看一区二区| 午夜福利18| 欧美黑人巨大hd| 午夜久久久久精精品| 色哟哟哟哟哟哟| www.999成人在线观看| 非洲黑人性xxxx精品又粗又长| 熟妇人妻久久中文字幕3abv| 国产高清视频在线观看网站| 成人鲁丝片一二三区免费| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇中文字幕五十中出| 欧美极品一区二区三区四区| 国产免费一级a男人的天堂| 欧美日韩中文字幕国产精品一区二区三区| 美女 人体艺术 gogo| 高潮久久久久久久久久久不卡| 女同久久另类99精品国产91| 美女免费视频网站| 亚洲av第一区精品v没综合| 国产精品伦人一区二区| 国内精品一区二区在线观看| 国产乱人视频| 国产一区二区激情短视频| 精品人妻熟女av久视频| 亚洲在线观看片| 午夜福利在线观看免费完整高清在 | 亚洲中文字幕一区二区三区有码在线看| 免费看a级黄色片| 听说在线观看完整版免费高清| 嫁个100分男人电影在线观看| aaaaa片日本免费| 日本在线视频免费播放| 免费观看精品视频网站| 精品一区二区三区视频在线观看免费| 精品人妻视频免费看| 深爱激情五月婷婷| 亚洲成人久久爱视频| 国产精品人妻久久久久久| 我的老师免费观看完整版| 90打野战视频偷拍视频| 国产日本99.免费观看| 欧美激情久久久久久爽电影| 最近视频中文字幕2019在线8| 别揉我奶头 嗯啊视频| 欧美又色又爽又黄视频| 最近视频中文字幕2019在线8| 又爽又黄a免费视频| 91午夜精品亚洲一区二区三区 | 校园春色视频在线观看| 欧美又色又爽又黄视频| 精品熟女少妇八av免费久了| 国产精品1区2区在线观看.| 国产一区二区三区在线臀色熟女| 看免费av毛片| 91麻豆精品激情在线观看国产| 麻豆av噜噜一区二区三区| 午夜福利在线观看免费完整高清在 | 九色国产91popny在线| 日本与韩国留学比较| 桃红色精品国产亚洲av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲无线在线观看| 欧美一级a爱片免费观看看| 日本五十路高清| 国产亚洲精品久久久com| 99国产精品一区二区三区| 长腿黑丝高跟| 午夜久久久久精精品| 少妇人妻一区二区三区视频| 在线看三级毛片| 全区人妻精品视频| 夜夜爽天天搞| 黄色配什么色好看| 性色av乱码一区二区三区2| 国产av麻豆久久久久久久| 他把我摸到了高潮在线观看| 99久国产av精品| 日韩人妻高清精品专区| 99在线视频只有这里精品首页| 老熟妇乱子伦视频在线观看| 91在线精品国自产拍蜜月| 真人做人爱边吃奶动态| 久久久精品欧美日韩精品| 在线播放国产精品三级| 亚洲片人在线观看| 成年女人毛片免费观看观看9| 999久久久精品免费观看国产| 99久久无色码亚洲精品果冻| 日本熟妇午夜| 免费看美女性在线毛片视频| 亚洲中文日韩欧美视频| 欧美精品国产亚洲| 级片在线观看| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 男人舔奶头视频| 欧美乱妇无乱码| 亚洲av五月六月丁香网| 我要搜黄色片| 免费av观看视频| 亚洲av成人av| 亚洲天堂国产精品一区在线| avwww免费| 国产精品亚洲av一区麻豆| 在线观看av片永久免费下载| 久久九九热精品免费| 网址你懂的国产日韩在线| 99国产综合亚洲精品| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 亚洲欧美日韩无卡精品| 成人鲁丝片一二三区免费| 成人国产一区最新在线观看| 国产美女午夜福利| 日本免费一区二区三区高清不卡| 国产人妻一区二区三区在| 久久精品91蜜桃| 亚洲精品影视一区二区三区av| 内射极品少妇av片p| 90打野战视频偷拍视频| 亚洲中文日韩欧美视频| 亚洲无线在线观看| 午夜免费男女啪啪视频观看 | 欧美在线黄色| 特大巨黑吊av在线直播| 深夜a级毛片| 免费在线观看亚洲国产| 日韩欧美国产一区二区入口| 少妇人妻精品综合一区二区 | 亚洲av第一区精品v没综合| 99久久99久久久精品蜜桃| 精品久久久久久久久av| 看免费av毛片| 亚洲精品粉嫩美女一区| 啦啦啦观看免费观看视频高清| 一进一出好大好爽视频| 国产精品永久免费网站| 搡女人真爽免费视频火全软件 | 一夜夜www| 国产熟女xx| 国产亚洲欧美在线一区二区| a级一级毛片免费在线观看| 欧美精品国产亚洲| 可以在线观看毛片的网站| 国产高清视频在线观看网站| 日韩大尺度精品在线看网址| 国产av不卡久久| 精品一区二区三区av网在线观看| 深夜精品福利| 免费在线观看日本一区| 色哟哟·www| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 日韩欧美国产在线观看| 在线免费观看的www视频| av黄色大香蕉| www.色视频.com| 亚洲av不卡在线观看| 成人三级黄色视频| 中出人妻视频一区二区| 97人妻精品一区二区三区麻豆| 99久久精品国产亚洲精品| 波多野结衣高清无吗| 好男人在线观看高清免费视频| 亚洲 欧美 日韩 在线 免费| 亚洲最大成人手机在线| 亚洲中文字幕一区二区三区有码在线看| 男人的好看免费观看在线视频| 久久国产精品影院| 午夜免费激情av| 亚洲成人免费电影在线观看| 亚洲成av人片免费观看| 人妻夜夜爽99麻豆av| 最近在线观看免费完整版| 999久久久精品免费观看国产| 午夜福利视频1000在线观看| 欧美日韩黄片免| 日日干狠狠操夜夜爽| 午夜日韩欧美国产| 免费电影在线观看免费观看| 男人狂女人下面高潮的视频| 麻豆久久精品国产亚洲av| 深夜a级毛片| 午夜免费激情av| 国内精品一区二区在线观看| 免费电影在线观看免费观看| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 国产一区二区三区在线臀色熟女| 成年人黄色毛片网站| 天堂av国产一区二区熟女人妻| 日本免费一区二区三区高清不卡| 色在线成人网| 日本a在线网址| 天堂网av新在线| 1000部很黄的大片| 色哟哟·www| 亚洲成a人片在线一区二区| 国产高清三级在线| 三级国产精品欧美在线观看| 精品久久国产蜜桃| 亚洲第一欧美日韩一区二区三区| 在线观看66精品国产| 免费av毛片视频| 亚洲人与动物交配视频| 熟女人妻精品中文字幕| 99国产精品一区二区三区| 国产在视频线在精品| 亚洲欧美日韩无卡精品| 欧美乱妇无乱码| .国产精品久久| 性色av乱码一区二区三区2| 国产探花极品一区二区| 天天躁日日操中文字幕| 亚洲专区国产一区二区| 性色avwww在线观看| 全区人妻精品视频| 亚洲人成电影免费在线| 欧美日本亚洲视频在线播放| 老司机福利观看| 久久婷婷人人爽人人干人人爱| 中文字幕人妻熟人妻熟丝袜美| 欧美午夜高清在线| 久久久精品大字幕| 欧美日韩中文字幕国产精品一区二区三区| 99久久99久久久精品蜜桃| 日韩欧美免费精品| 亚洲经典国产精华液单 | 久久久久久九九精品二区国产| 午夜激情福利司机影院| 欧美日本亚洲视频在线播放| 99热这里只有精品一区| 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 亚洲狠狠婷婷综合久久图片| 日本熟妇午夜| 午夜老司机福利剧场| 999久久久精品免费观看国产| 国产人妻一区二区三区在| 午夜福利在线在线| 国产精品久久久久久久久免 | 国内精品久久久久久久电影| 一级a爱片免费观看的视频| 亚洲国产精品久久男人天堂| 国产高清视频在线观看网站| 免费av观看视频| 天美传媒精品一区二区| 日本成人三级电影网站| 一区二区三区四区激情视频 | 伦理电影大哥的女人| 成年版毛片免费区| 欧美一区二区亚洲| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 少妇人妻精品综合一区二区 | 在线国产一区二区在线| 美女黄网站色视频| 欧美另类亚洲清纯唯美| 久久国产乱子免费精品| 少妇人妻精品综合一区二区 | 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 国产精品亚洲美女久久久| 十八禁网站免费在线| 国产精品久久电影中文字幕| 村上凉子中文字幕在线| 能在线免费观看的黄片| 999久久久精品免费观看国产| 成人性生交大片免费视频hd| 日韩欧美精品v在线| 99热这里只有精品一区| 亚洲av成人av| 欧美乱色亚洲激情| 中亚洲国语对白在线视频| 国产一区二区三区视频了| 国产麻豆成人av免费视频| 久久国产乱子免费精品| 国产色爽女视频免费观看| 亚洲国产欧美人成| 少妇熟女aⅴ在线视频| 十八禁人妻一区二区| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 丰满的人妻完整版| 免费人成在线观看视频色| 欧美一区二区亚洲| 在线观看免费视频日本深夜| 51国产日韩欧美| 精华霜和精华液先用哪个| 级片在线观看| 露出奶头的视频| 激情在线观看视频在线高清| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 级片在线观看| 露出奶头的视频| 欧美不卡视频在线免费观看| 一本一本综合久久| 午夜福利18| 级片在线观看| 淫秽高清视频在线观看| 国产毛片a区久久久久| 久久九九热精品免费| 99国产极品粉嫩在线观看| 日韩欧美在线二视频| 99在线人妻在线中文字幕| 特级一级黄色大片| 午夜福利视频1000在线观看| 免费看a级黄色片| 丰满人妻熟妇乱又伦精品不卡| 精品久久国产蜜桃| 成人国产综合亚洲| 麻豆国产av国片精品| 亚洲欧美日韩无卡精品| 久久国产精品人妻蜜桃| 极品教师在线免费播放| 欧美区成人在线视频| 国产中年淑女户外野战色| 国内揄拍国产精品人妻在线| 丁香六月欧美| 91在线精品国自产拍蜜月| 国产极品精品免费视频能看的| 国产亚洲欧美98| 亚洲av.av天堂| 9191精品国产免费久久| 中亚洲国语对白在线视频| 亚洲专区中文字幕在线| 亚洲久久久久久中文字幕| 变态另类丝袜制服| 搞女人的毛片| a在线观看视频网站| 嫩草影视91久久| 蜜桃久久精品国产亚洲av| 精品久久久久久久久av| 伦理电影大哥的女人| 别揉我奶头~嗯~啊~动态视频| 婷婷色综合大香蕉| 啦啦啦韩国在线观看视频| 国产高清有码在线观看视频| 两个人的视频大全免费| 91麻豆av在线| 九色成人免费人妻av| 亚洲va日本ⅴa欧美va伊人久久| 日本一二三区视频观看| 免费看a级黄色片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 免费看光身美女| 真人一进一出gif抽搐免费| 亚洲成人中文字幕在线播放| 99热精品在线国产| 啪啪无遮挡十八禁网站| 日韩人妻高清精品专区| 丰满人妻熟妇乱又伦精品不卡| 综合色av麻豆| 国产视频一区二区在线看| 国产日本99.免费观看| 免费电影在线观看免费观看| 欧美日韩国产亚洲二区| 欧美成狂野欧美在线观看| 一级a爱片免费观看的视频| 国产色婷婷99| 丝袜美腿在线中文| 欧美+亚洲+日韩+国产| 色尼玛亚洲综合影院| 欧美性猛交黑人性爽| 小说图片视频综合网站| 国产69精品久久久久777片| 亚洲自偷自拍三级| 亚洲国产精品合色在线| 国产精品美女特级片免费视频播放器| 男女做爰动态图高潮gif福利片| 赤兔流量卡办理| 久久午夜福利片| 2021天堂中文幕一二区在线观| 最近最新中文字幕大全电影3| 欧美日韩乱码在线| 两个人的视频大全免费| 国产成人啪精品午夜网站| 90打野战视频偷拍视频| 中文亚洲av片在线观看爽| 久久精品国产99精品国产亚洲性色| 88av欧美| 国产精品一区二区三区四区免费观看 | 一进一出好大好爽视频| 国产成+人综合+亚洲专区| 欧美色视频一区免费| 日本 av在线| 91在线观看av| 一本一本综合久久| 一区福利在线观看| 赤兔流量卡办理| 一本久久中文字幕| 国产不卡一卡二| 亚洲三级黄色毛片| 日本 欧美在线| 久久人人精品亚洲av| 日韩大尺度精品在线看网址| 啪啪无遮挡十八禁网站| 亚洲真实伦在线观看| 亚洲欧美日韩高清在线视频| 丝袜美腿在线中文| av中文乱码字幕在线| 午夜精品一区二区三区免费看| 深爱激情五月婷婷| 每晚都被弄得嗷嗷叫到高潮| 国产精品永久免费网站| 欧美日韩乱码在线| 日韩免费av在线播放| 中文在线观看免费www的网站| 日本a在线网址| 成人国产一区最新在线观看| 国产高清三级在线| 丝袜美腿在线中文| 亚洲真实伦在线观看| 啦啦啦韩国在线观看视频| 国模一区二区三区四区视频| 亚洲国产精品sss在线观看|