余江妹 余 俊 龔學余 謝寶藝 張 能

1（南華大學 數理學院 衡陽 421001）

2（南華大學 核科學技術學院 衡陽 421001）

托卡馬克等離子體極向轉動對帶狀流的影響

余江妹1余 俊1龔學余2謝寶藝1張 能1

1（南華大學 數理學院 衡陽 421001）

2（南華大學 核科學技術學院 衡陽 421001）

采用理想磁流體力學模型，利用絕熱近似方法，推導出托卡馬克等離子體極向轉動體系中帶狀流的色散方程，分析了極向轉動對帶狀流產生的影響。結果表明，等離子體極向轉動體系下，測地聲模頻率和聲波頻率都會隨極向轉動馬赫數而增大。帶狀流的頻率與采用等溫近似時有一定的數值差別；測地聲模的密度擾動量形式也發(fā)生變化，由靜態(tài)體系中純駐波的形式，變?yōu)橛梢粋€駐波和一個小振幅行波疊加的形式。

托卡馬克，帶狀流，頻率，等離子體極向轉動

在托卡馬克等離子體中，帶狀流是一種環(huán)形和極向都對稱的靜電勢漲落結構[1–2]。帶狀流有兩個分支：低頻帶狀流和測地聲模。低頻帶狀流的頻率接近0，因此又稱零頻帶狀流。測地聲模是帶狀流的高頻分支，其頻率ω～cs/R，其中cs為聲速。廣泛的理論與實驗研究表明，帶狀流是由漂移波湍流的非線性作用激發(fā)而產生，可以從湍流中獲得能量，并通過流剪切來抑制湍流所導致的反常輸運[3–5]。因此帶狀流的研究對托卡馬克聚變裝置的穩(wěn)態(tài)運行有著重要意義。

在帶狀流抑制湍流的過程中，經常會伴隨著等離子體的轉動，可以說，轉動在托卡馬克等離子體中是一種普遍存在的本質現象。對于等離子體環(huán)向轉動，王少杰等[6]理論研究了等離子體環(huán)向轉動對帶狀流的模特征的影響，結果表明，測地聲模的頻率和低頻帶狀流的頻率都隨著環(huán)向轉動馬赫數而增大，且密度擾動量的形式也與靜態(tài)等離子體中有所不同，由靜態(tài)等離子體中的駐波變?yōu)檗D動體系中行波的形式。對于等離子體極向轉動，其產生原因主要是強徑向電場、stringer自旋加快效應以及外部動量注入等[2,7]，目前國內外的研究熱點主要集中在等離子體轉動對各種等離子體模的影響上[7–10]。周登[7]采用等溫近似，理論研究了托卡馬克等離子體極向轉動體系中，帶狀流的頻率隨極向轉動馬赫數的變化情況。Pamela等[10]研究了等離子體極向轉動體系下的邊界局域模，得出了等離子體極向轉動對邊界局域模穩(wěn)定性的影響。

本文采用理想磁流體力學模型，利用絕熱近似方法，推導出帶狀流在極向轉動體系下的色散方程。通過對色散方程的求解，得出低頻帶狀流和測地聲模頻率的解析式，以及帶狀流頻率隨極向轉動馬赫數的變化情況，并分析得出托卡馬克等離子體極向轉動對帶狀流的影響。

1 理論模型

采用理想磁流體模型，基本方程組為：

式中，ρ、P和u分別是等離子體質量密度、壓力和速度；Φ、B和j分別是電勢、磁場強度和電流密度；γ是絕熱系數。

在磁流體力學中，托卡馬克等離子體速度的一般形式為[7]：

式中，K(ψ)和Φ'都是磁面ψ的函數。將平衡速度分成極向速度和環(huán)向速度，即，其中?θ和 ?φ分別為極向方向單位矢量和環(huán)向方向單位矢量。則環(huán)向速度uT=ρ–1KBT+RΦ'，極向速度up=ρ–1KBp，這說明當等離子體極向速度存在時，環(huán)向速度必然也存在，即單純的等離子體極向轉動是不存在的，環(huán)向轉動必然會伴隨著極向轉動而存在。本文采用文獻[7]所提出的轉動體系，即在環(huán)向速度中的ρ–1KBT項被另外一項RΦ'幾乎抵消掉的情況下，考慮極向轉動對帶狀流所產生的效應。

假設ρ0、p0、u0、Φ0、j0為平衡量，當體系處于靜態(tài)時，式(2)可寫作：

將式(2)分別乘以磁場B0和▽φ，可以得到：

將密度ρ0寫作的形式，由式(7)和(8)可得：

式中，Mp=BTup/Bpcs為極向轉動馬赫數，cs為聲速；，q為安全因子。

引入等離子體質量密度擾動量ρ1、壓力擾動量p1、溫度擾動量T1、速度擾動量v、靜電勢擾動量Φ1、電流密度擾動量j1，假設壓力擾動量，等離子體比壓β～ε2。

其中，式(16)為動量方程式(12)平行磁場方向的分量方程，式(17)由動量方程作磁面平均而得到。

式(15)–(18)可簡化為：

2 結果分析

對式(15)進行求解，可以得出色散方程：

方程的解為：

式(28)，10ω=，為低頻帶狀流。式(29)是由低頻帶狀流衍生出的模，其頻率隨極向轉動馬赫數而增加，且與極向轉動馬赫數成簡單的正比關系。式(30)中的兩個解，分別表示極向轉動情況下的測地聲模(+)和聲波(–)。

圖1為絕熱系數γ=5/3時，不同的安全因子下，測地聲模的歸一化頻率隨極向轉動馬赫數的變化情況。由圖1，測地聲模的頻率隨極向轉動馬赫數而逐漸增大。圖2為絕熱系數γ=5/3時，聲波的歸一化頻率隨極向轉動馬赫數的變化情況，與測地聲模相類似，聲波的頻率也隨著極向轉動馬赫數而逐漸增大。

由圖1和圖2可看出，在等離子體極向轉動體系中，測地聲模的頻率和聲波的頻率與靜態(tài)等離子體中有所不同，都會隨著極向轉動馬赫數而增大。本文采用絕熱近似方法，所得結果與文獻[7]采用等溫近似所得結論相類似，說明不論采用絕熱近似，還是等溫近似，在極向轉動體系中測地聲模和聲波的頻率都會隨著極向轉動馬赫數而增大。

圖1 GAM頻率隨極向轉動馬赫數的變化Fig.1 Frequencies of GAM plotted with respect to the poloidal Mach number.

圖2 SW頻率隨極向轉動馬赫數的變化Fig.2 Frequencies of the sound wave plotted with respect to the poloidal Mach number.

但采用不同的近似方法，帶狀流的頻率會有不同的數值結果。圖3為安全因子q=5時，絕熱系數γ=1和γ=5/3的情況下，測地聲模歸一化頻率隨極向轉動馬赫數的變化情況。圖4為安全因子為5時，絕熱系數取不同數值情況下，聲波的歸一化頻率隨極向轉動馬赫數的變化情況。由圖3和圖4可以看出，當絕熱系數取不同數值時，即采用不同的近似方法時，測地聲模頻率和聲波的頻率都會有一定的數值差別。

圖3 GAM頻率變化比對Fig.3 Frequency comparison of GAM.

圖4 SW頻率變化比對Fig.4 Frequency comparison of the sound wave.

測地聲模密度擾動量形式也發(fā)生一定的變化，與靜態(tài)等離子體中有所不同。將測地聲模密度擾動量寫作cos部分和sin部分，這兩部分之間的比值與呈一定的關系：

圖5表示測地聲模密度擾動量cos部分與sin部分的比值，隨著極向轉動馬赫數的變化情況。當時，即體系處于靜態(tài)時，式(31)的值為0，這表明測地聲模的密度擾動量為純駐波的形式。在等離子體極向轉動體系中，由圖5可看出，式(31)為一個極小的值，這與在靜態(tài)等離子體體系中值為0有所不同，表明測地聲模的密度擾動量為一個駐波和一個小振幅行波相疊加的形式。

圖5 GAM密度擾動量cos部分和sin部分的比值Fig.5 Ratio between the cos and sin components of density perturbation for GAM.

3 結語

本文采用磁流體力學模型，絕熱近似方法，理論研究了托卡馬克等離子體極向轉動對帶狀流的影響。結果表明在等離子體極向轉動體系下，測地聲模頻率和聲波頻率都隨極向轉動馬赫數而增大。采用絕熱近似方法時，壓力擾動量的形式與采用等溫近似時的并不完全相同，因此得到的帶狀流頻率會有數值上的差別。此外，測地聲模的密度擾動量形式也與靜態(tài)等離子體時有所不同，在靜態(tài)等離子體中，測地聲模的密度擾動量為駐波的形式，而在極向轉動等離子體體系中，測地聲模的密度擾動量是一個駐波和一個小振幅行波相疊加的形式。

1 Diamond P H, Itoh S I, Itoh K, et al. Zonal flows in plasma- a r eview[J]. P lasma Physics and Controlled Fusion, 2005, 47(5): R35–R50

2 Hassam A B, Drake J F. Spontaneous poloidal spin-up of tokamak plasmas: reduced equations, physical mechanism, and sonic regimes[J]. Ph ysics Plasmas, 1993, 5(11): 4022–4029

3 Yu J, Gong X Y. Geodesic acoustic mode excitation by drift waves in tokamak plasmas with toroidal rotation[J]. Nuclear Fusion, 2013, 53(12): 123027

4 Yu J, Gong X Y, Chen Y. Lo w frequency zonal flow excitation by drift waves in tokamak plasmas with toroidal rotation[J]. Nuclear Fusion, 2015, 55(2): 023015

5 Guo W F, Wang S J, Li J G. Damping effects of finite parallel thermal conductivity o n zonal f lows[J]. Plasma Physics and Controlled Fusion, 2008, 50(9): 095006

6 Wang S J. Zonal flows in tokamak plasmas with toroidal rotation[J]. Physical Review Letters, 2006, 97(8): 085002

7 Zhou D. Zon al flow modes in a tokamak plasma with dominantly poloidal mean f lows[J]. Phy sics Pla smas, 2010, 17(10): 102505

8 Elfimov A G, Galvao R M O, Sgalla R J F. Rotation effect on geodesic and zonal f low modes in tokamak plasmas with isothermal magnetic surfaces[J]. Plasma Physics and Controlled Fusion, 2011, 53(10): 105003

9 Elfimov A G. Alfvén continuum deformation by kinetic geodesic ef fect in rotating tokamak plasmas[J]. Ph ysics Plasmas, 2010, 17(2): 022102

10 Pamela S, Huysmans G, Benkadda S. Influence of poloidal equ ilibrium rotat ion i n MHD simulations of edge-localized modes[J]. Plasma Physics and Controlled Fusion, 2010, 52(7): 075006

CLC TL61+2.4

Effects of poloidal rotation on zonal flows in tokamak plasmas

YU Jiangmei1YU Jun1GONG Xueyu2XIE Baoyi1ZHANG Neng1
1(School of Mathematics and Physics, University of South China, Hengyang 421001, China)
2(School of Nuclear Science and Technology, University of South China, Hengyang 421001, China)

Background: Zonal flows are the electrostatic modes with structures of toroidal symmetry and poloidal symmetry, which play an important role in restraining turbulences. In the process of turbulence suppression, the poloidal rotation widely exists in tokamak plasmas. Purpose: In this study, the effects of the poloidal rotation on zonal flows in tokamak plasmas are concluded. Methods: With the dispersion relation derived by using ideal MHD model with the adiabatic approximation, the frequencies of zonal flows are deduced. Results: Firstly, the frequency of the geodesic acoustic mode (GAM) increases with respect to the poloidal Mach number when adiabatic coefficient is 5/3. Secondly, the frequency of the sound wave (SW) increases with respect to the poloidal Mach number when adiabatic coefficient is 5/3. Third ly, the frequency of GAM with γ=1 i s dif ferent from th at when γ=5/3. Th e fourth, the frequency of SW with γ=1 is different from that when γ=5/3. The last, the ratio between the cos and sin components of density perturbation for GAM is a very small value in the poloidally rotating plasma. Conclusion: It is shown that the frequencies of the zonal flows with the adiabatic approximation, increase with respect to the poloidal Mach number, but differ in numerical value of the frequencies of the zonal flows with the isothermal approximation. In contrast to the pure standing wave form in static plasmas, the density perturbations consist of a standing wave superimposed with a small amplitude traveling wave in the poloidally rotating plasma.

Tokamak, Zonal flows, Frequency, The poloidal rotation of plasmas

TL61+2.4

10.11889/j.0253-3219.2015.hjs.38.100602

No.11105071、No.11375085、No.11475083、No.41104094、No.11405082)、南華大學核聚變與等離子體物理創(chuàng)新團隊(No.NHXTD03)資助

余江妹，女，1989年出生，2012年畢業(yè)于洛陽師范學院，現為南華大學碩士研究生，研究領域為核聚變與等離子體物理

2015-06-19，

2015-08-23