徑向電場對離子溫度梯度模穩定性的影響*

陳凝飛 魏廣宇 仇志勇2)?

1) (浙江大學物理學院,聚變理論與模擬中心,杭州 310027)

2) (Center for Nonlinear Plasma Science and ENEA C.R.Frascati,Frascati,Italy)

為了理解托卡馬克裝置中給定徑向電場對離子溫度梯度模(ITG)穩定性的影響,基于非線性回旋動理學理論和氣球模表象推導了環位形下包含徑向電場引起的極向流和密度擾動影響的ITG 的本征模方程,并分別在長/短波長極限下研究了高能量粒子誘發測地聲模(EGAM)所伴隨的徑向電場對ITG 的本征頻率、增長率和平行模結構的影響.不僅對該本征模方程進行了理論研究,還使用本征矩陣法對其進行數值求解,以便對理論結果進行驗證.研究發現EGAM 伴隨的電場引起的極向轉動會大幅降低ITG 的增長率,而極向模數m=1 的密度擾動對ITG 的線性穩定性影響很小.這一結果與一般認為的帶狀流通過極向流剪切抑制湍流的結果是一致的.除此之外,使用本文發展的一般性方法也可以研究高能量粒子激發的阿爾芬不穩定性與漂移波湍流通過阿爾芬不穩定性激發帶狀結構發生的間接非線性相互作用,即帶狀結構所伴隨的徑向電場通過極向轉動和密度擾動影響ITG 的穩定性.該間接非線性通道可以作為對主導背景等離子體輸運的微觀湍流和主導高能量粒子輸運的阿爾芬不穩定性之間的直接相互作用通道的補充.

1 引言

由磁約束聚變裝置中普遍存在的壓強梯度驅動的漂移波湍流(drift wave,DW)被認為是引起等離子體反常輸運的可能原因之一[1],從而導致等離子體約束性能降低.在眾多漂移波模式中,由離子溫度梯度驅動的離子溫度梯度模(ion-temperature gradient driven mode,ITG)是被研究得最多的漂移波之一,這是由于它會造成離子的熱輸運,而離子的約束對于未來磁約束裝置極為重要.因此,需要對ITG 的線性穩定和非線性飽和機制進行研究.其中,ITG 非線性飽和機制的研究對提升磁約束聚變裝置的熱離子約束性能和加深對低-高約束模轉換的理解具有重要意義.

徑向電場引起的E×B極向剪切流可以破壞ITG 的徑向相干結構,從而降低其幅度及相應的離子熱輸運[2].Hahm等[3]通過兩點非線性理論分別研究了大縱橫比或有限縱橫比的托卡馬克位形中的E×B極向剪切流對湍流的抑制作用,并發現當流剪切率大于背景湍流的去相關率時(正比于增長率),湍流就會被大幅抑制.隨后,他們還發現雖然極向流的振蕩分量(如本文中重點討論的測地聲模)對剪切率的貢獻很大,但是它對背景湍流的抑制效果不如零頻帶狀流.這可能是由于有限頻率的剪切流在湍流渦旋還未被扭曲時就改變方向,從而使其平均效果沒有零頻分量強.但是,這些徑向電場與DW 相互作用的模型都沒有包含可能會產生定性影響的環耦合效應.因此,需要研究環耦合效應存在時,徑向電場對ITG 穩定性和平行模結構的影響.

這些徑向電場有多種來源,包括湍流自發激發的帶狀流(zonal flows,ZFs)和外源產生的徑向電場.帶狀流是環向模數n=0、極向模數m ≈0 的介觀尺度徑向電場.它包含零頻帶狀流(zero-frequency ZF,ZFZF)和有限頻率的測地聲模(geodesic acoustic mode,GAM).ITG 可以通過非線性雷諾脅強自發激發帶狀流,這是ITG 非線性飽和的重要通道.在該過程中ITG 一方面會將一部分能量非線性地傳遞給帶狀流,另一方面,也會被帶狀流散射到線性穩定的短波長區間[4,5].因此,ITG 具有自調制的過程,即ITG 自發激發的帶狀流會對ITG自身性質產生影響.另外,未來及現有的眾多托卡馬克中通過中性束加熱或聚變反應產生的高能量粒子(energetic particle,EPs)對自持燃燒具有重要意義.由于GAM 具有有限頻率,可被速度空間各向異性的高能量粒子共振激發,即高能量粒子誘發測地聲模(energetic-particle-induced GAM,EGAM)[6?9].由于帶狀流可以通過非線性雷諾脅強抑制湍流,Qiu等[10]提出用中性束注入激發EGAM,進而對湍流進行主動控制.但是,Zarzoso等[11]的回旋動理學模擬表明,EGAM 被EPs 注入激發后,反而將ITG 湍流從Dimits 臨界穩定區間中激發起來[12],從而破壞等離子體約束,這與人們普遍認為的帶狀流抑制湍流的圖像矛盾.一個可能的機制是,對于臨界穩定的背景湍流,當EGAM 被EPs 強烈激發時,耦合的非線性系統將會把能量傳遞給背景湍流,這一機制可類比化學中的勒夏特列原理[5,13,14].

除此之外,EPs 的特征拉莫半徑尺度與ITG湍流尺度之比?1 (其中TE為EP 溫度,Ti為離子溫度),因此之前都認為EPs 和微觀湍流之間的相互作用較弱.同時,未來燃燒等離子體中產生的α 粒子的速度拋射角各向同性,因此,EGAM與湍流之間的相互作用和主要依靠α 粒子加熱的燃燒等離子體不完全相關.但是大量實驗[15,16]和模擬結果發現EPs 的注入可能會對湍流幅度進行調制,使得背景等離子體約束改善[17?19].針對這一現象,人們提出了各種假說,其中高能量粒子直接或間接產生的徑向電場可能會在其中發揮重要作用.一方面,EPs 可以通過軌道損失機制直接產生邊緣徑向電場[20],以此抑制湍流及其輸運,并可能導致低-高約束模轉換.另一方面,未來的燃燒等離子體中,EPs 的速度接近阿爾芬速度vA=其中ρ=nimi為等離子體質量密度(mi是離子質量,ni是離子密度),因此會激發各種阿爾芬本征模式(Alfvén eigenmodes,AEs)[21]或高能量粒子模(energetic particle mode,EPM)[22].這些阿爾芬本征模式可能通過參量衰變或調制不穩定性激發帶狀結構(zonal structures,ZSs)[23,24],并通過帶狀結構與ITG 發生間接相互作用.這一間接非線性相互作用通道與二者的直接相互作用通道一起構成了主導熱等離子體輸運的微觀湍流和主導高能量粒子輸運的阿爾芬湍流的非線性相互作用圖像[25].因此,在環位形中研究ITG 與徑向電場,尤其是與高能量粒子激發的徑向電場之間的非線性相互作用對理解ITG 的非線性飽和,以及ITG 與AEs 之間的間接相互作用都具有重要意義.

綜上所述,本文基于非線性回旋動理學理論在環位形中研究給定的徑向電場對ITG“線性”穩定性的影響.首先推導了描述一般徑向電場影響下的ITG“線性”局域本征模方程,并將其變換到氣球模空間以便對環耦合效應進行處理.這里的局域指本文只考慮沿磁力線的平行模結構,并且包含了環效應和平行可壓縮性,但忽略了徑向波包效應.研究發現,在本模型中,徑向電場通過密度擾動和電勢擾動引起的極向轉動對ITG 的線性穩定性產生影響.為了契合本文的主題,我們以EGAM 的徑向電場為例進行了具體分析.研究發現無論在短波長還是長波長極限下,EGAM 的電勢擾動引起的極向轉動都會使得ITG 的線性增長率和實頻的絕對值大幅降低.與此同時,EGAM 的密度擾動對ITG 的線性穩定性影響較小.本文的結構安排如下: 第2 節引入一般理論模型;第3 節在短波長和長波長極限下求解并分析ITG 的本征模方程;第4 節給出簡單的總結和討論.本文發展的一般模型,適用于頻率遠低于ITG 頻率的徑向電場對其穩定性的影響,因此,不僅適用于研究線性激發的帶狀流對ITG 穩定性的影響,也可以應用于理解平衡的平均流與湍流的相互作用,還可以應用于高能量粒子激發的阿爾芬不穩定性與漂移波湍流通過帶狀流結構的“間接”相互作用.

2 基本模型

為了突出徑向電場對ITG 穩定性的影響,本文采用具有同心圓磁面的托卡馬克位形,并使用右手系(r,θ,ξ),其中r,θ 和ξ 分別是小半徑、極向和環向角.在此位形中,平衡磁場由B=B0[(1-εcosθ)eξ+εeθ/q] 給出,其中ε ≡r/R是逆環徑比,R是大半徑,q ≡rBξ/(RBθ) 是安全因子,它代表磁力線在極向和環向環繞的圈數之比.由于ITG一般具有較高的環向模數n,平衡參數剖面(密度、溫度等)的特征長度遠大于有理面之間的距離(Δr=1/(nq')),因此可以將擾動電勢寫為

其中,s ≡(r-r0)/Δr=nq-m0表示用相鄰有理面間距離Δr=1/(nq') 歸一化的徑向坐標,q'=?q/?r表示安全因子的徑向梯度,r0為參考磁面的徑向位置,它滿足nq(r0)-m0=0,且滿足m0為遠大于j的整數;表示ITG 的徑向精細結構.

ITG 和EGAM 都是頻率遠低于回旋頻率的低頻擾動,因此二者之間的非線性相互作用可用非線性回旋動理學理論進行研究.為了更加關注離子溫度梯度模的物理,這里假設托卡馬克密度剖面是平的,即ηi≡Lni/LTi→∞,其中Lni≡-ni/(?ni/?r),LTi≡-Ti/(?Ti/?r)分別為密度和溫度不均勻性的特征長度.同時,假設在高約束模等離子體的臺基區內,徑向電場的頻率一般遠低于ITG 的頻率.描述ITG 離子響應的非線性回旋動理學方程由下式給出:

其中,左右兩邊的第一項分別為電子和離子的絕熱響應,即對應無慣性的玻爾茲曼分布,〈·〉 代表速度空間 積分,N0表示粒子數密度,Ts表示粒子s的溫度.典型的ITG 模通常滿足k//v//,e?ω～ω?Ti?ωd,k//v//,i,即忽略電子的慣性,其擾動分布函數的非絕熱分量 δHI,e=0.在該量級下,離子的擾動分布函數的非絕熱分量為

其中,Λ ≡ickrkθ/B0為非線性耦合系數.將ITG的離子非絕熱響應(4)代入準中性條件(3),即可得到ITG 的色散關系.同時,垂直波數可以分解為徑向和極向分量,即,則第j個極向模式在實空間內的本征模方程可寫為

在方程(5)中,s?≡r(?q/?r)/q是磁剪切,τ ≡Te/Ti是電子和離子溫度之比,bθ ≡,z ≡s-j=nq-m.方程(5)右邊第一項為磁漂移引起的相鄰極向模式之間的環耦合效應,這是托卡馬克這樣的環形磁約束裝置上特有的效應.另外,離子對徑向電場的非絕熱響應 〈δHE,i〉 可能也有極向不均勻性,因此可能造成額外的環耦合.例如,對于本文研究的EGAM 來說,其密度擾動在環向截面為上下反對稱,正比于 sinθ,而ZFZF 的密度擾動正比于 cosθ[26].方程(5)是一個涉及徑向和極向的二維微分-差分方程,很難進行理論分析.為了更好地研究環耦合效應帶來的影響,采用氣球模表象將方程(5)分解為兩個耦合的一維本征值方程,即描述ITG 平行模結構的方程和緩變的徑向包絡的方程.氣球模表象基于平移不變性,即有理面之間的距離 Δr遠小于等離子體平衡參數變化的尺度或,可以認為不同有理面上的極向分量的模結構形狀一致.此時,擾動電勢可以寫為傅里葉積分的形式:

其中,η 為沿磁力線的極向坐標,又被稱為擴展極向角,其取值范圍為η ∈(-∞,+∞) ;Φ(η) 為ITG的平行模結構;A(r)為極向分量的緩變徑向包絡.本文將基于局域模型研究ITG 的徑向模結構,并不包含其徑向波包效應[27],因此也不會涉及EGAM的全局模結構[8].基于(6)式給出的傅里葉積分形式,ITG 在氣球模空間的本征模方程可寫為

3 徑向電場對ITG 線性穩定性的影響

基于方程(8)給定的EGAM 的徑向電場,方程(7)可進一步寫為如下類薛定諤方程的形式:

3.1 短波長極限

從文獻 [28]可知,ITG 的平行模結構在實空間的寬度正比于b1/2,這說明短波長極限下(b～O(1)),ITG 的模結構在徑向較寬,相鄰極向模式之間的耦合較強,即在氣球模空間中局域在η=0 附近,形成氣球模結構.因此,短波長極限也被稱為強耦合/氣球近似(strong coupling/ballooning approximation).在這一近似下,周期性勢阱可在η=0 附近泰勒展開.值得注意的是,EGAM 正比于 sinη的貢獻會改變原方程的偶對稱性,從而使ITG 的本征模結構不再關于η=0 對稱.當偏離過大時,上述在原點附近進行的展開將可能不再適用.這里先驗地假設該偏離 Δη ?1,從而使得上述展開依然有效.此時,ITG 的本征模方程可寫為

其中,λ=(2l+1) 為該本征值方程的本征值,l為非負整數.將方程(11)代入方程(12)即可得出能級為l的本征態的色散關系:

此時,韋伯方程中的諧振子勢阱存在一系列本征值和本征函數,如圖1 所示.從圖1(a)和圖1(b)可以看出,無論是否存在徑向電場,最不穩定的本征態(本征值虛部最大)都是l=0 的基態,其本征函數由高斯函數給出為Φ(η)=exp[-σ(η+Δη)2],其中

圖1 短波長極限下,eδφE/Ti=0 (a)和 0.1 (b)時,ITG 本征值的實部(Ωr)和虛部(Ωi)的分布.在兩種情況下,最不穩定的都是l=0 的基態Fig.1.Distribution of the real (Ωr) and imaginary (Ωi)parts of eigenvalues of ITG when eδφE/Ti=0 (a) and 0.1(b) in the short-wavelength limit.In both cases,the ground state with l=0 is the most unstable eigenstate.

因此,在不存在徑向電場時(Δη=0),ITG的基態本征模在氣球模空間的平行模結構是局域在原點附近的高斯函數,其寬度正比于1/b,即實空間模寬度正比于b;相反,在系統中存在有限的徑向電場時(Δη/=0),ITG 的平行模結構相對原點有有限的偏移,如圖2 所示.此外,從圖2 發現由EGAM 引起的ITG 平行模結構的偏移 Δη ?1,這說明了前文中將勢阱在η=0 附近展開的自洽性.同時,基態對應的色散關系(l=0)為

圖2 短波長極限下,eδφE/Ti=0.1 和0 時,ITG 最不穩定的基態的平行模結構.其中,藍色實線表示徑向電場為零時ITG 的模結構,紅色虛線表示有限徑向電場時ITG 的模結構Fig.2.Mode structure of the most unstable mode of ITG.The blue solid and red dashed lines represent the cases with eδφE/Ti=0 and eδφE/Ti=0.1,respectively.

色散關系(14)除了大括號里的最后兩項表示由EGAM 的徑向電場引起的修正,其他項與文獻 [28]中的ITG 線性色散關系一致.通過求解基態色散關系(14),可以得到基態本征值的解析解,它的實部和虛部分別為ITG 的實頻和增長率;同時,本文也使用本征矩陣法直接求解方程(9)從而得到本征值的數值解,并將二者所得的結果進行比較.圖3(a)和圖3(b)分別是ITG 的增長率和實頻對EGAM幅度的依賴關系,可以看到紅線代表的理論解和藍線代表的數值解符合得很好.從圖3(a)可以看出,當徑向電場從0 增強到0.1 時,ITG 的增長率大幅下降約50%,同時,實頻的絕對值也下降約50%.從方程(2)可知,徑向電場對ITG 穩定性的影響有兩個通道: 一是電勢擾動引起的極向旋轉;二是密度擾動.為了研究兩個效應各自對ITG 增長率降低的貢獻,可在研究一個效應時將另一個效應關閉.首先,當開啟密度擾動,同時關閉電勢擾動時,ITG 的增長率基本不隨EGAM 幅度變化而變化,如圖4 中的綠色實線所示.考慮相反的情形,即只包含極向旋轉的影響而不包含密度擾動時,ITG的增長率與將二者都包含的情況十分接近,分別如圖4 中的紅色虛線和藍色實線所示.綜上,我們發現徑向電場的電勢擾動(即極向旋轉)是引起ITG 增長率下降的主要原因.

圖3 短波長極限下,ITG 的增長率 (a)和實頻(b)與EGAM幅度 eδφE/Ti 的關系.藍色圓點表示直接數值求解本征值方程(9)得到的數值解,紅色叉號表示求解基態色散關系(14)得到的理論解.實頻和增長率都是用 Cs/ 進行歸一化,其中,=2Te/mi 表示聲速;這兩幅圖所用的參數為=0.2,b=1Fig.3.Dependence of the growth rate (a) and real frequency (b) of ITG on the amplitude of EGAM eδφE/Ti in short-wavelength limit.The blue circles represent the numerical value obtained by directly solving the eigenmode equation (9);while the red crosses represent the theoretical value obtained by solving the dispersion relation (14).The real frequency and growth rate of ITG are normalized to Cs/,with =2Te/mi representing sound velocity.The parameters used here are =0.2 and b=1.

3.2 長波長極限

3.1 節討論的強耦合近似要求b ?1,這是一個很強的約束條件.當b逐漸減小時,ITG 極向分量的徑向寬度也逐漸減小,因此在b?1 的長波長極限下,強耦合近似不再成立.由于在長波長近似下ITG 極向分量的實空間寬度很小,相鄰極向分量之間的耦合相比強耦合近似弱,因此也被稱為中等/弱耦合近似.長波長近似下存在平板ITG 和環形ITG 兩個分支,而本文更關注激發閾值更低的環形分支.環形ITG 的平行模結構是環效應引起的快變勢阱(η～O(1))的響應和慢變的平衡剖面的響應的疊加.自洽性分析給出環分支是由平行可壓縮性(正比于η 二階導的項)與絕熱電子響應(正比于Ω3的項)平衡得到,由此可以得出Ω=O(b-1/3).此時,可以定義環分支的模結構為

圖5 長波長極限下,eδφE/Ti=0 (a)和 0.1 (b)時,ITG 的本征值分布.在兩種情況下,最不穩定的都是l=0 的基態Fig.5.Distribution of eigenvalues of ITG wheneδφE/Ti=0(a) and 0.1 (b) in the long-wavelength limit.In both cases,the ground state with l=0 is the most unstable eigenstate.

方程(19)左邊第二項由徑向電場電勢m=0分量貢獻,而密度擾動對長波長極限下的ITG 色散關系沒有貢獻,這可能是由于在長波長極限下ITG的平行模結構很寬,如圖6 所示,從而將EGAM的m=1 的密度擾動平均了.圖7 給出了ITG 的增長率和頻率對徑向電場幅度的依賴關系,可以發現理論解和數值解符合得較好.具體來說,當徑向電場強度eδφE/Ti從0 增強到0.1 時,ITG 的增長率從0.028 下降到0.020,下降約 28%.與此相比,短波長情形下,ITG 的增長率大幅下降約50%.這說明,在同樣幅度下,徑向電場對短波長極限的ITG 的致穩效果比對長波長極限的ITG 強,但是我們并不清楚造成該現象的原因.方程(19)中不包含密度調制的貢獻,這說明長波長極限下密度調制可能不會對ITG 的增長率造成影響.同樣,解析的色散關系(19)說明密度擾動不對ITG 的線性色散關系造成影響,為了驗證這一點,采用與短波長極限時相同的做法,即在研究一個效應時將另一個關閉,結果如圖8 所示.結果表明,只存在極向旋轉和兩種非線性效應都存在時,ITG 增長率對EGAM 幅度的依賴關系基本一致,同時,只存在密度擾動時,ITG 增長率不隨EGAM 幅度變化.因此,在長波長極限下,EGAM 的電勢擾動會導致ITG 增長率下降,而密度擾動對ITG 的增長率沒有影響.長波長極限下,ITG 的平行模結構的偶對稱性也會被徑向電場破壞,使其從壞曲率區產生微小的偏移,如圖6 所示.值得一提的是,這些結果是在b=0.01 時得到的,這會導致ITG的頻率小于EGAM 的頻率,從而破壞理論的自洽性,但是為了將快慢尺度充分分離,依然采用了這一非典型參數.

圖6 長波長極限下,eδφE/Ti=0.1 和0 時,ITG 最不穩定的基態的平行模結構.其中,藍色實線表示徑向電場為零時ITG 的模結構,紅色虛線表示有限徑向電場時ITG 的模結構Fig.6.Mode structure of the most unstable mode of ITG.The blue solid and red dashed lines represent the cases with eδφE/Ti=0 and eδφE/Ti=0.1,respectively.

圖7 長波長極限下,ITG 的增長率 (a)和實頻(b)與EGAM 幅度 eδφE/Ti 的關系.藍色圓點表示直接數值求解本征值方程(9) 得到的數值解,紅色叉號表示求解基態色散關系(19)得到的理論解.這兩幅圖所用的參數為εTi=0.2,b=0.01Fig.7.Dependence of the growth rate (a) and real frequency (b) of ITG on the amplitude of EGAM eδφE/Ti in long-wavelength limit.The blue circles represent the numerical value obtained by directly solving the eigenmode equation (9);while the red crosses represent the theoretical value obtained by solving the dispersion relation (19).The parameters used here are εTi=0.2 and b=0.01.

圖8 長波長極限下,極向旋轉和密度調制共同作用和分別作用時 ITG 增長率對 EGAM 幅度的依賴關系.圖中綠色實線表示只有密度擾動的情形,紅色虛線表示只有電勢擾動的情形,藍色實線表示兩種擾動都存在的情形Fig.8.Dependence of the growth rate of ITG in presence of poloidal rotation and/or density modulation in longwavelength limit.The green solid and red dashed lines represent the cases with only density modulation and poloidal rotation,respectively;while the blue line represents the case with both effects.

4 結論

為了研究徑向電場對ITG“線性”穩定性的影響,本文基于非線性回旋動理學理論和氣球模表象,在環位形下推導了可以描述在一般徑向電場影響下的ITG 本征模方程.該模型適用于頻率遠低于ITG 的徑向電場.另外,為了理解高能量粒子激發的EGAM 的電場對ITG 的影響,將EGAM 作為示例進行研究,并在短/長波長極限下分別對其進行解析求解.不僅對該本征模方程進行理論研究,還基于本征矩陣法對其進行數值求解以便對解析結果進行驗證.

從一般模型可知徑向電場主要通過密度擾動和電勢擾動引起的極向轉動對ITG 造成影響.無論是在短波長還是長波長極限,EGAM 的電勢擾動引起的極向轉動都會導致ITG 的線性增長率和實頻的絕對值大幅降低.同時,EGAM 的密度擾動對ITG 的線性穩定性幾乎沒有影響,尤其在長波長極限下,快變的密度擾動在慢變平行模結構尺度上被平均,從而不會進入ITG 的色散關系.另外,EGAM 的密度擾動正比于 sinη,因此會破壞原有勢阱的偶對稱性,并最終導致ITG 的平行模結構偏離η=0.雖然本文得到的結果與一般的極向剪切流抑制湍流的模型所得到的的結果一致,但是本文中極向流可以直接降低ITG 的線性增長率,且不依賴流的剪切.

本研究不僅給出了徑向電場為EGAM 時ITG的增長率、實頻和平行模結構等的具體結果,也提供了研究徑向電場對ITG 不穩定性的影響的一般理論模型,這些徑向電場不僅可以是由外部注入的高能量粒子直接激發的電場,也可以是由阿爾芬本征模等不穩定性自發激發的介觀尺度帶狀結構.帶狀結構同樣可通過極向轉動和密度擾動對ITG 的穩定性產生影響,因此,ITG 與高能量粒子激發的阿爾芬不穩定性可以通過介觀尺度帶狀結構為中介發生間接相互作用.該間接相互作用通道將與兩者的直接非線性相互作用通道一起構成主導背景等離子體輸運的微觀湍流和主導高能量粒子輸運的阿爾芬不穩定性之間的相對完整的非線性作用圖像,有助于人們理解等離子體中多尺度的非線性相互作用.同時,本工作也存在一些不足之處.首先,所做的分析都是假設了ITG 的頻率遠大于徑向電場的頻率,這一點對靠近芯部的EGAM 來說并不一直成立.因此,未來可以將此模型推廣到有限頻率的徑向電場的情形.另外,只考慮了徑向的局域理論,沒有研究ITG 的徑向包絡效應的影響.但是,托卡馬克中的徑向電場具有徑向不均勻性,如本文中研究的EGAM 就具有全局模結構.因此,本文給出的局域模型無法包含徑向電場的梯度及其產生的極向流的剪切對ITG 造成的影響.因此,需要繼續在本工作的基礎上包含ITG 的徑向包絡并進而研究徑向電場的梯度和流剪切對ITG 穩定性的影響.

感謝浙江大學陳騮教授分享的原始想法.