林麥麥 王明月 蔣蕾

(西北師范大學(xué)物理與電子工程學(xué)院,蘭州 730070)

采用Sagdeev 勢方法研究了多組分塵埃等離子體中的非線性塵埃聲孤波的傳播特征.推導(dǎo)得到了在含有塵埃顆粒、電子流、質(zhì)子流及Kappa 電子和離子的多組分塵埃等離子體中的非線性塵埃聲波所對應(yīng)的Sagdeev 勢函數(shù).并利用定性分析方法,確定該系統(tǒng)同時存在線性周期波解軌道、非線性周期波解軌道和孤立波解軌道.與此同時,借助數(shù)值模擬方法發(fā)現(xiàn):該多組分塵埃等離子體系統(tǒng)中僅存在振幅小于零的稀疏型孤立波,且非線性塵埃聲孤波的振幅、寬度及波形等傳播特征均與馬赫數(shù)M、多種粒子數(shù)密度、溫度、電荷量、質(zhì)量及Kappa 電子數(shù)、離子數(shù)等各種系統(tǒng)參數(shù)存在緊密關(guān)聯(lián).

1 引言

目前,關(guān)于多組分塵埃等離子體中多種非線性波動過程的研究已成為等離子體物理的重要前沿研究領(lǐng)域,其研究內(nèi)容廣泛,包括天體物理學(xué)、半導(dǎo)體制造、聚變反應(yīng)堆等[1?8].實際條件下塵埃顆粒廣泛存在于實驗室等離子體和星際空間、太陽系、地球電離層,以及暫星尾和行星環(huán)空間等離子體環(huán)境中[9?14].1990 年Rao等[15]首次從理論上預(yù)測了未磁化塵埃等離子體中的非線性塵埃聲波的存在性;1995 年,Barkan等[16]在實驗室中觀測到了非線性塵埃聲波.現(xiàn)在,非線性塵埃聲波的相關(guān)研究課題已成為了等離子體物理學(xué)的熱門課題之一,眾多科研工作者對該課題進行了深入探討,Tiwari等[17]研究了塵埃等離子體中大振幅的非線性離子聲孤波的傳播特點;Adhikary等[18]利用約化攝動法研究了有正、負離子的塵埃等離子體中小振幅的非線性塵埃聲孤波和沖擊波的傳播特征;Mamun[19]依靠贗勢方法研究了絕熱條件下熱塵埃等離子體中非線性塵埃聲孤立波的特點;Sinha 和Sahu[20]研究了非磁化四分量等離子體中的塵埃離子聲波.

研究塵埃等離子體中非線性波動過程的方法比較多樣,常見的約化攝動法可以用于研究小振幅的非線性波;如果考慮邊界效應(yīng),則可利用贗勢法進行非線性波動過程的討論;除此之外,Sagdeev勢方法可以研究等離子體中大振幅的非線性波動過程,已廣泛應(yīng)用于多種課題的理論研究中.例如,Mamun 和Shukla[21]利用Sagdeev 勢方法研究了強耦合塵埃等離子體中任意振幅的塵埃聲孤波;Hatami 和Nikknam[22]則使用Sagdeev 勢方法研究了雙溫非廣延電子對離子雙層結(jié)構(gòu)的影響;El-Hanbaly等[23]利用Sagdeev 勢方法研究了具有超熱電子和非廣延離子的塵埃等離子體中非線性塵埃聲波的傳播特征;Sebastian等[24]則通過Sagdeev勢方法研究了五分量等離子體中壓縮孤立波和稀疏孤立波的存在及孤立波的傳播特征.由此不難看出,Sagdeev 勢方法是等離子體中非線性波研究的重要方法.

本文利用Sagdeev 勢方法研究了非線性塵埃聲波在含有帶正電的塵埃顆粒、電子流、質(zhì)子流、呈Kappa 分布的電子和離子多組分塵埃等離子體中的基本傳播特征.第2 節(jié)給出原始方程并利用Sagdeev 勢方法推導(dǎo)得到二維自治動力系統(tǒng)的具體表達式.第3 節(jié)對系統(tǒng)相圖進行分析討論,確定該系統(tǒng)同時存在線性波解、非線性波解和孤立波解,且多種波動行為的存在與多種系統(tǒng)因素存在緊密關(guān)聯(lián).第4 節(jié)求得Sagdeev 勢函數(shù)和孤立波解表達式,并借助數(shù)值模擬的方法獲得Sagdeev 勢的基本演化特征及孤立波的振幅、寬度和波形隨系統(tǒng)參數(shù)的基本變化規(guī)律.從而得出結(jié)論:多組分塵埃等離子體中的多種系統(tǒng)參數(shù),諸如塵埃顆粒、電子、質(zhì)子、離子等多種粒子的質(zhì)量、數(shù)密度、溫度、電荷量等均對該系統(tǒng)中的非線性塵埃聲孤波的傳播特征存在重要影響.

2 原始方程及其求解

為了研究多組分塵埃等離子體中非線性塵埃聲波的傳播特征,假設(shè)多組分塵埃等離子體是由太陽風(fēng)中帶正電的塵埃顆粒、電子流、質(zhì)子流,以及呈現(xiàn)Kappa 分布的電子和離子共同組成.以磁流體力學(xué)中理想磁流體力學(xué)(magneto hydro dynamic,MHD)模型為理論基礎(chǔ),忽略粒子間的相互碰撞,對塵埃等離子體受到的庫侖力、重力、拖曳力等作用力暫時忽略不計,可以得到如下形式的無量綱化方程組[25]:

其中,Kappa 分布的離子和電子的數(shù)密度形式如下[25]:

在這里Kappa 分布表達式僅適用于κe(i)＞3/2 和極限κe(i)→∞的條件[25],其中κe和κi分別為Kappa分布電子數(shù)和離子數(shù).

接下來,利用Sagdeev 勢方法研究非線性塵埃聲波的傳播特征.首先,對(1)式—(6)式中的自變量進行行波變換:令ξx-Mt,M為馬赫數(shù).(1)式—(6)式可以化為

考慮邊界條件:當ξ→∞時,Nd→1,Nc→1,Nb→1,ud→0,uc→0,ub→0,φ→0,則通過化簡(10)式—(15)式可得Nd,Nc和Nb的表達式:

將其代入(7)式后可得

其中

將(16)式轉(zhuǎn)化為二維自治動力系統(tǒng)如下:

其中,φ表示靜電勢,ψ表示靜電勢φ對ξ的一階導(dǎo)數(shù).

3 系統(tǒng)相圖分析

為了研究多組分塵埃等離子體系統(tǒng)的動力學(xué)特征,利用數(shù)值模擬方法研究不同系統(tǒng)參數(shù)對非線性二維自治動力系統(tǒng)(17)相圖的具體影響.

圖1 給出了二維自治動力系統(tǒng)(17)的相圖隨著系統(tǒng)參數(shù)馬赫數(shù)M的變化規(guī)律,其他系統(tǒng)參數(shù)取值分別為[25,26]:Hb0.8,Hc0.5,μb0.7,μc0.6,δe0.08,δi0.1,σd0.0009,κiκe2.0 .圖1(a)顯示:當馬赫數(shù)M=1.2 時,多組分塵埃等離子體系統(tǒng)中同時存在線性周期波軌道、非線性周期波軌道和同宿軌道.這表明該系統(tǒng)同時存在線性周期波解、非線性周期波解和孤立波解.與此同時,圖1(b)和圖1(c)顯示:隨著馬赫數(shù)M的逐步增大,系統(tǒng)的不同軌道特征不會改變,但多種不同軌道的存在范圍會依次增大.

圖1 系統(tǒng)相圖隨馬赫數(shù)的變化(a)M=1.2;(b)M=1.3;(c)M=1.4Fig.1.Variations of system phase diagram with Mach number:(a)M=1.2;(b)M=1.3;(c)M=1.4.

二維自治動力系統(tǒng)(17)的相圖隨著系統(tǒng)參數(shù)Hb的變化規(guī)律見圖2,其中,Hc0.5,M1.2,μb0.7,μc0.6,δe0.08,δi0.1,σd0.0009,κiκe2.0.結(jié)果顯示:當Hb變化時,系統(tǒng)中始終同時存在線性周期軌道、非線性周期軌道和同宿軌道,但各種軌道存在的范圍有所變化.與此同時,考 慮到Hb,而μbmd/(Zd0mb),γb(nb0Zd0)/nd0,σbTb/Teff,這表明多組分塵埃等離子體中的各種系統(tǒng)參數(shù),諸如塵埃顆粒和質(zhì)子流的質(zhì)量、數(shù)密度、電荷量和溫度等,均對自治系統(tǒng)的相軌線存在一定的影響.

圖2 系統(tǒng)相圖隨Hb 的變化(a)Hb=0.4;(b)Hb=0.6;(c)Hb=0.8Fig.2.Variations of system phase diagram withHb:(a)Hb=0.4;(b)Hb=0.6;(c)Hb=0.8.

圖3 給出了二維自治動力系統(tǒng)(17)的相圖隨著系統(tǒng)參數(shù)Hc的變化規(guī)律,其中Hb0.8,M1.2,μb0.7,μc0.6,δe0.08,δi0.1,σd0.0009,κiκe2.0.圖3(a)顯示:當Hc0.4 時,多組分塵埃等離子體系統(tǒng)中仍然同時存在線性周期波軌道、非線性周期波軌道和同宿軌道.通過對圖3(a)—(c)進行對比,系統(tǒng)相圖顯示:隨著Hc逐步增大,系統(tǒng)的不同軌道特征不會改變,但多種不同軌道的存在范圍會依次增大.在這里Hc,而μcmd/(Zd0me),γc(nd0Zd0)/nc0,σcTc/Teff,這意味著對于含有Kappa 電子和離子的多組分塵埃等離子體而言,系統(tǒng)平衡狀態(tài)時的塵埃顆粒和電子的質(zhì)量、數(shù)密度、電荷量和溫度的不同,只會使系統(tǒng)相圖中不同類型的塵埃聲波解的存在范圍發(fā)生改變,而不會改變系統(tǒng)相圖中不同類型塵埃聲波解共同存在的基本屬性.

圖3 系統(tǒng)相圖隨Hc 的變化(a)Hc=0.4;(b)Hc=0.6;(c)Hc=0.8Fig.3.Variations of system phase diagram withHc:(a)Hc=0.4;(b)Hc=0.6;(c)Hc=0.8.

圖4 給出了系統(tǒng)相圖隨著系統(tǒng)參數(shù)μb的變化規(guī)律,其他系統(tǒng)參數(shù)取值分別為:μc0.6,M1.2,Hb0.8,Hc0.5,δe0.08,δi0.1,σd0.0009,κiκe2.0 .圖4 顯示:當系統(tǒng)參數(shù)μb取值不同時,多組分塵埃等離子體系統(tǒng)中仍然同時存在線性周期波軌道、非線性周期波軌道和同宿軌道,而且隨著μb的逐步增大,系統(tǒng)的多種軌道的存在范圍會發(fā)生改變,而μbmd/(Zd0mb),這說明當多組分塵埃等離子體中所含的塵埃顆粒質(zhì)量較大而質(zhì)子流的質(zhì)量較小、且平衡狀態(tài)下塵埃顆粒的電荷量較小時,該自治系統(tǒng)的線性周期波解、非線性周期波解和孤立波解的存在范圍將減小.

圖5 給出了系統(tǒng)相圖隨著參數(shù)μc的變化規(guī)律,其中μb0.7,M=1.2,Hb0.8,Hc0.5,δe0.08,δi0.1,σd0.0009,κiκe2.0 .可以看出,當系統(tǒng)參數(shù)μc從0.3 增加到0.6 再增加到0.9 時,系統(tǒng)相圖中不同類型的塵埃聲波解的存在范圍隨著系統(tǒng)參數(shù)μc的增大而增大,這與圖4 系統(tǒng)相圖隨著系統(tǒng)參數(shù)μb的變化有所不同,但系統(tǒng)相圖中的線性周期波解、非線性周期波解和孤立波解三種不同類型的波解結(jié)構(gòu)不會發(fā)生改變.

圖4 系統(tǒng)相圖隨μb 的變化(a)μb=0.3;(b)μb=0.6;(c)μb=0.9Fig.4.Variations of system phase diagram withμb :(a)μb=0.3;(b)μb=0.6;(c)μb=0.9.

圖5 系統(tǒng)相圖隨μc 的變化(a)μc=0.3;(b)μc=0.6;(c)μc=0.9Fig.5.Variations of system phase diagram withμc:(a)μc=0.3;(b)μc=0.6;(c)μc=0.9.

圖6 給出了自治動力系統(tǒng)(17)的相圖隨著系統(tǒng)參數(shù)σd的變化規(guī)律,其他系統(tǒng)參數(shù)取值分別為:M1.2,Hb0.8,Hc0.5,μb0.7,μc0.6,δe0.08,δi0.1,κiκe2.0 .不難看出,當其他參數(shù)取定值時,而多組分塵埃等離子體中所含的塵埃顆粒的溫度與有效溫度的比值σd取值不同時,系統(tǒng)相圖中各種軌道的存在范圍發(fā)生變化,由于σdTd/Teff,也就是說塵埃顆粒的溫度越高,系統(tǒng)相圖中線性周期波軌道、非線性周期波軌道和同宿軌道的存在范圍越小.

圖6 系統(tǒng)相圖隨σd 的變化(a)σd=0.0004;(b)σd=0.0006;(c)σd=0.0008Fig.6.Variations of system phase diagram withσd:(a)σd=0.0004;(b)σd=0.0006;(c)σd=0.0008.

圖7 給出了系統(tǒng)參數(shù)κi對自治動力系統(tǒng)相圖的影響,其中M1.2,Hb0.8,Hc0.5,μb0.7,μc0.6,δe0.08,δi0.1,σd0.0009,κe2.0 .從圖7 可以看到,當Kappa 分布離子數(shù)κi分別為2.0,2.5,3.0 時,該自治系統(tǒng)同時存在線性周期波軌道、非線性周期波軌道和同宿軌道,但多種不同軌道的存在范圍會依次減小.二維自治動力系統(tǒng)的相圖隨系統(tǒng)參數(shù)κe的變化規(guī)律由圖8 給出,其他系統(tǒng)參數(shù)分別為:M1.2,Hb0.8,Hc0.5,μb0.7,μc0.6,δe0.08,δi0.1,σd0.0009,κi2.0.結(jié)果表明:當系統(tǒng)參數(shù)κe取值不同時,該系統(tǒng)同時存在線性周期波解、非線性周期波解和孤立波解,不難看出,當Kappa 分布電子數(shù)κe從2.0 增加到2.5 時,系統(tǒng)相圖中不同類型的塵埃聲波解的存在范圍增大,但當系統(tǒng)參數(shù)κe從2.5 增加到3.0 時,系統(tǒng)相圖中塵埃聲波解的存在范圍變化不大.

圖7 系統(tǒng)相圖隨κi 的變化(a)κi=2.0;(b)κi=2.5;(c)κi=3.0Fig.7.Variations of system phase diagram withκi:(a)κi=2.0;(b)κi=2.5;(c)κi=3.0.

圖8 系統(tǒng)相圖隨κe 的變化(a)κe=2.0;(b)κe=2.5;(c)κe=3.0Fig.8.Variations of system phase diagram withκe:(a)κe=2.0;(b)κe=2.5;(c)κe=3.0.

4 Sagdeev 勢函數(shù)及孤立波解

為了研究多組分塵埃等離子體中的非線性塵埃聲孤波的傳播特征,將(16)式等式兩邊同乘以dφ/dξ,并對ξ積分一次可得Sagdeev 勢方程:

其中Sagdeev 勢函數(shù)為

接下來,通過數(shù)值模擬對Sagdeev 勢進行分析討論,從而進一步明確多組分塵埃等離子體的多種系統(tǒng)參數(shù)對非線性塵埃聲波的具體影響.圖9 給出了Sagdeev 勢函數(shù)隨不同系統(tǒng)參數(shù)的變化規(guī)律.可以看出存在φm＜0,使得0且,這意味著該系統(tǒng)僅存在振幅小于0的稀疏型孤立波.

圖9(a)給出了系統(tǒng)參數(shù)Hb對Sagdeev勢V(φ)的具體影響,結(jié)果表明孤立波的振幅φm隨著系統(tǒng)參數(shù)Hb的增加而減小.考慮到Hb,而μbmd/(Zd0mb),γb(nb0Zd0)/nd0,σbTb/Teff,這說明多組分塵埃等離子體中的各種系統(tǒng)參數(shù),諸如塵埃顆粒和質(zhì)子流的質(zhì)量、數(shù)密度、電荷量和溫度等系統(tǒng)因素,均對稀疏孤立波的振幅存在顯著影響.圖9(b)顯示隨著參數(shù)Hc的逐步增加,孤立波的振幅φm將依次增大.而考慮到Hc,其中μcmd/(Zd0me),γc(nd0Zd0)/nc0,σcTc/Teff,這表明當多組分塵埃等離子體中所含的塵埃顆粒質(zhì)量大、數(shù)密度高且電荷量低時,該系統(tǒng)中的非線性塵埃聲孤波的振幅就較大;而當電子質(zhì)量大、電子流數(shù)密度高且溫度高時,該系統(tǒng)中的非線性塵埃聲孤波則振幅較小.圖9(c)—(f)給出了系統(tǒng)參數(shù)μb,μc,M和σd取不同值時,Sagdeev 勢函數(shù)的變化規(guī)律,由圖9(c)和圖9(f)可知,孤立波的振幅φm隨著系統(tǒng)參數(shù)μb和σd的增大而減小,這與圖9(a)中的結(jié)果相同.而圖9(f)說明σdTd/Teff越大孤波振幅越小,即塵埃顆粒的溫度越高孤立波的振幅φm越小.同理從圖9(d)和圖9(e)可知,孤立波振幅φm的變化情況與圖9(b)類似,即隨著系統(tǒng)參數(shù)μc和馬赫數(shù)M的增大,孤立波的振幅φm也隨之增大.

圖9 Sagdeev勢V(φ)隨不同參數(shù)的變化規(guī)律Fig.9.Variations of Sagdeev potentialV(φ)with different parameters.

圖10 給出了Sagdeev勢V(φ)隨Kappa 分布電子數(shù)κe和離子數(shù)κi的變化規(guī)律,表明多組分塵埃等離子體系統(tǒng)中所含的Kappa 電子分布數(shù)κe的增加和Kappa 離子分布數(shù)κi的減少,將使孤立波的振幅φm增大.這說明多組分塵埃等離子體系統(tǒng)中電子和離子的Kappa 分布形式對Sagdeev 勢函數(shù)及稀疏型孤立波的傳播特征均存在重要影響.接下來對Sagdeev 勢方程(18)進行求解,可得到非線性塵埃聲孤波的形式如下:

其中θkξ-ω0τ,k為ξ方向的波數(shù),ω0為波速;φm-2A1/(3A2)和D(2/A1)-1/2分別為孤立波的振幅和寬度.

圖11 給出了多組分塵埃等離子體系統(tǒng)中非線性塵埃聲孤波的波形隨著系統(tǒng)參數(shù)的變化趨勢.結(jié)果顯示:當參數(shù)Hb,Hc,μb,μc,M和σd取值不同時,該系統(tǒng)僅存在振幅小于0 的稀疏型孤立波.該結(jié)論與Sagdeev 勢函數(shù)的數(shù)值模擬結(jié)果保持一致.圖11(a),(c),(f)表示,當系統(tǒng)參數(shù)Hb,μb和σd增加時,孤立波振幅減小而寬度增大;而當系統(tǒng)參數(shù)Hc,μc和M增加時,孤立波的振幅增大且寬度減小.這些結(jié)論與圖10 中Sagdeev 勢函數(shù)所反映的結(jié)論相一致.

圖10 Sagdeev勢V(φ)隨κe和κi 的變化規(guī)律Fig.10.Variations of Sagdeev potentialV(φ)withκeandκi .

圖11 孤立波φ 的波形隨不同參數(shù)的變化規(guī)律Fig.11.Waveform variation law of the solitary wavesφ with different parameters.

非線性塵埃聲孤波的波形φ隨電子和離子Kappa 分布數(shù)κe和κi的變化規(guī)律如圖12 所示,孤立波的振幅φm隨著κe的增大而增大、隨著κi的增大而減小.這表明當多組分塵埃等離子體系統(tǒng)中含有較多的Kappa 分布電子和較少的Kappa 分布離子時,非線性塵埃聲孤波將具有更大的振幅.

圖12 孤立波φ 的波形變化規(guī)律Fig.12.Waveform variation law of the solitary wavesφ .

5 結(jié)論

本文研究了由塵埃顆粒、電子流、質(zhì)子流、Kappa 電子和離子所組成的多組分塵埃等離子體中非線性塵埃聲孤波的傳播特征.利用Sagdeev 勢方法求解得到二維自治動力系統(tǒng)具體表示形式,并通過數(shù)值模擬方法得到相圖,結(jié)果顯示:多組分塵埃等離子體系統(tǒng)中同時存在線性波解軌道、非線性波解軌道和同宿軌道;而Sagdeev 勢函數(shù)的變化規(guī)律則顯示該系統(tǒng)僅存在振幅小于0 的稀疏型孤立波.通過進一步孤立波波形特征的討論,不難發(fā)現(xiàn)含有Kappa 電子和離子分布的多組分塵埃等離子體系統(tǒng)中的多種系統(tǒng)參數(shù)對非線性塵埃聲孤波的振幅、寬度和波形等傳播特征均存在不可忽視的重要影響.