肖　湧

(浙江大學 聚變理論與模擬中心，浙江 杭州 310027)

?

磁約束聚變等離子體的大規(guī)?；匦齽永韺W模擬

肖湧

(浙江大學 聚變理論與模擬中心，浙江 杭州 310027)

【摘要】文章介紹了磁約束核聚變能源科學及其密切相關(guān)的聚變等離子體物理，以及研究磁約束聚變等離子體非線性物理所必需的大規(guī)模回旋動理學模擬的方法及其前沿進展.

【關(guān)鍵詞】等離子體；核聚變；磁約束；大規(guī)模數(shù)值模擬；回旋動理學

1聚變能源科學

1.1聚變能源

核聚變能源是資源取之不盡、清潔的理想能源，將徹底解決人類社會能源問題.如圖1，太陽通過氫原子聚變生成氦原子釋放能量，在太陽核心部位，每秒鐘約有6.2億噸的氫原子發(fā)生聚變反應. 核聚變能源具有以下幾個特點.

1.1.1能量巨大

核能是能量巨大的能源，一次核反應釋放的能量是MeV量級，而平時我們廣泛使用的煤、石油、天然氣等化石能源涉及的是化學反應，釋放的能量是eV量級，只有前者的千分之一.核能有重核裂變和輕核聚變兩種方式，其中核聚變比核裂變釋放出更多的能量.例如，鈾235的裂變反應，將千分之一的物質(zhì)變成了能量；而氘的聚變反應，將近千分之四的物質(zhì)變成了能量.

1.1.2資源豐富

重核裂變使用的主要原料是鈾，目前探明的儲量僅夠使用約1 000年；而輕核聚變使用的燃料是海水中的氘，1升海水能提取30毫克氘，在聚變反應中能產(chǎn)生約等于300升汽油的能量，即“1升海水約等于300升汽油”，地球上海水中就有45萬億噸氘，足夠人類使用百億年.聚變所需的另外一種燃料是氚，在地球上儲量并不多，但是可以通過鋰轉(zhuǎn)變獲得，而地球上尤其是月球上都有大量的鋰礦.一旦氘氚聚變實現(xiàn)商用，我們也可以采用其他聚變過程.總之，聚變能基本可以算是取之不盡的能源.

1.1.3成本低廉

1千克氘的價格只為1千克濃縮鈾的1/40.

1.1.4安全、無污染

核聚變不會產(chǎn)生長半衰期的放射性污染物，萬一發(fā)生事故，反應堆會自動冷卻而停止反應，不會發(fā)生爆炸.

圖1　太陽的聚變反應示意圖Figure 1　Fusion reaction schematic of the sun

1.2磁約束聚變

1.2.1約束方式

聚變原料發(fā)生聚變反應時，最佳溫度高達一千萬度.在這種高溫下，電子被完全剝離出原子，形成物質(zhì)第四態(tài)等離子體.當沒有約束時，帶電粒子的熱運動會使等離子體迅速擴散，無法發(fā)生核聚變反應，所以需要外界約束使高溫等離子體聚在一起，保持一定的密度，從而發(fā)生可觀的核聚變反應.等離子體目前有兩種約束方式，磁場約束和慣性約束.

1) 磁約束

磁場可以約束粒子，使其沿磁力線運動，如圖2.等離子體作為帶電粒子，會受到電磁場的影響，所以可以通過外加磁場約束等離子體.磁場是目前已知最好的絕緣體，其隔熱性能比航天飛機的瓦片還要強1 000倍.

圖2　磁場對粒子的約束Figure 2　Constraint of magnetic field on the particles

2) 磁約束中存在的問題

磁場只能約束粒子在垂直磁場方向的運動，但并不能約束粒子在平行磁場方向的運動，如圖3，帶電粒子在垂直磁場的方向做回旋運動，磁場不均勻性會產(chǎn)生梯度漂移.為此可以將磁場彎曲成環(huán)形，首尾相接，如圖4，這樣就可以把粒子約束在環(huán)中運動.

圖3　帶電粒子在磁場中的運動圖Figure 3　Movement of charged particles in magnetic field

圖4　環(huán)形磁場中螺旋形狀磁場線的方向Figure 4　Direction of magnetic field lines with circular spiral shape in magnetic field

但是，當把磁場彎曲成環(huán)形時，磁場中存在垂直于磁場方向的梯度，粒子會發(fā)生垂直于磁場方向的梯度漂移.同時粒子沿著環(huán)形磁力線上的運動，會受到沿磁場曲率方向的離心力作用，從而產(chǎn)生一個垂直于離心力和外加磁場為方向的曲率漂移.因此需要再增加一個沿極向的磁場，使得粒子在不同位置處的漂移作用相互抵消，從而保證約束的穩(wěn)定.

目前國際上比較主流的磁約束聚變裝置有托卡馬克和仿星器等，如圖5(a).其中托卡馬克由于其磁場設計的相對簡單和約束性能的良好成為各國進行磁約束聚變研究的首選.仿星器有在穩(wěn)態(tài)運行方面有其獨特的優(yōu)點，但是其復雜磁場的設計和建造也極有挑戰(zhàn)性.

3) 慣性約束

除了磁約束的方法之外，還可以使用慣性來約束等離子體.宇宙中的許多星體都是依靠自身巨大質(zhì)量產(chǎn)生的引力作用來約束粒子的，如圖5(b)所示，上面為太陽依靠其巨大的慣性引力約束等離子體，下面為星云依靠其巨大的慣性引力約束等離子體.在實驗室里可以通過強激光壓縮氘氚靶丸實現(xiàn)慣性約束，如美國的國家點火裝置(NIF)和我國的神光系列裝置都是采用這一原理的.

圖5　等離子體的約束Figure 5　Plasma confinement

1.2.2磁約束裝置

1) 托卡馬克(Tokamak)磁約束聚變實驗裝置

托卡馬克是前蘇聯(lián)科學家發(fā)明的一種利用磁約束來實現(xiàn)磁約束聚變的環(huán)性容器，依靠環(huán)形線圈和等離子體電流產(chǎn)生的環(huán)狀螺旋形強磁場，將極高溫等離子體狀態(tài)的聚變物質(zhì)約束在環(huán)形容器里，以此實現(xiàn)聚變反應.

2) 國際磁約束核聚變實驗裝置

自從20世紀50年代以來，人類開始研究聚變能源及其密切相關(guān)的等離子體物理，測試了很多不同的實驗裝置，發(fā)現(xiàn)到目前為止托卡馬克還是表現(xiàn)最好的.評判聚變裝置約束效果的重要參數(shù)包括聚變等離子體的溫度、密度以及約束時間，理論分析表明只有當這三者的乘積超過某一臨界值時，燃料才能燃燒發(fā)生聚變，這就是著名的勞遜判據(jù).自20世紀90年代以來，歐盟的JET、美國的TFTR、DIII-D和日本的JT-60、JT-60SA這幾個大型托卡馬克裝置在磁約束核聚變研究中獲得許多重要成果.圖6為不同托卡馬克試驗裝置大半徑對比圖.我們現(xiàn)在的實驗裝置越來越靠近勞遜判據(jù)的臨界值，如圖7和圖8.另外，一個經(jīng)驗判定是，托卡馬克尺寸越大越穩(wěn)定，約束時間越長，越接近勞遜判據(jù)的臨界值，這就是為什么現(xiàn)在托卡馬克越做越大的原因.

圖6　不同托卡馬克試驗裝置大半徑對比Figure 6　Major radius comparison among different Tokamak experiment devices

圖7　不同托卡馬克實驗裝置約束時間對比Figure 7　Confinement time comparison among different Tokamak experiment devices

圖8　不同托卡馬克實驗裝置約束效果對比Figure 8　Confinement effect comparison among different Tokamak experiment devices

隨著研究進展的深入，各國之間的學術(shù)交流愈加頻繁，國際熱核聚變實驗反應堆(International Thermonuclear Experimental Reactor，ITER)將成為世界上最大的磁約束等離子體物理實驗裝置，如圖9. 這是目前正在建設的世界上最大的實驗性托卡馬克核聚變反應堆，由七個成員實體資助和運行，包括歐盟、印度、日本、中國、俄羅斯、韓國和美國.它建立在TFTR、歐洲聯(lián)合環(huán)形加速器(JET)、JT-60和T-15等裝置的研究之上，并將顯著的超越所有前者.ITER實驗成功后將是人類歷史上第一個長時間實現(xiàn)核聚變反應的實驗裝置，對整個聚變能源科學的影響是巨大的.

國際熱核反應試驗堆計劃(ITER)是商用化核聚變能源的關(guān)鍵一步.作為ITER成員，中國在國家中長期科技規(guī)劃中把ITER計劃作為重點高技術(shù)，投入近百億人民幣進行研究.我國關(guān)于ITER采購包的研制在我國科研和工程技術(shù)人員的協(xié)同努力下進展順利，在七方成員中排在前列，體現(xiàn)了我國的科學技術(shù)實力和綜合國力.

圖9　國際熱核反應試驗堆(ITER)結(jié)構(gòu)示意圖Figure 9　Structure of ITER

ITER計劃中的等離子體的參數(shù)：

·等離子體大半徑：6.2 m

·等離子體小半徑：2.0 m

·等離子體容量：840 m3

·等離子體電流：15.0 MA

·軸向環(huán)形磁場強度：5.3T

·聚變功率：500 MW

·等離子體維持時間：>400s

·能量倍增因數(shù)(Q值)：>10

(3) 國內(nèi)磁約束核聚變實驗裝置

目前，國內(nèi)主要的托卡馬克實驗裝置有中科院等離子體所的EAST(圖10)和核工業(yè)部西南物理研究院的HL-2A(圖11).

圖10　中科院等離子體所EAST超導托卡馬克Figure 10　EAST superconducting tokamak in Institute of Plasma Physics, Chinese Academy of Sciences

圖11　核工業(yè)部西南物理研究院HL-2A托卡馬克Figure 11　HL-2A tokamak in Southwestern Institute of Physics, Department of Nuclear Industry

先進實驗超導托卡馬克實驗裝置(Experimental Advanced Superconducting Tokamak，EAST)是一個具有非圓小截面的全超導托卡馬克.它可以實現(xiàn)高參數(shù)和長脈沖的實驗運行；科研人員在以上條件下可以開發(fā)先進運行模式，并進行熱流平衡和粒子流平衡控制的實驗研究.

環(huán)流器二號A(HL-2A)是我國第一個具有先進偏濾器位形的非圓截面的托卡馬克核聚變實驗研究裝置，可以進行高參數(shù)等離子體條件下的改善約束實驗，并利用其獨特的大體積封閉偏濾器結(jié)構(gòu)，開展核聚變領(lǐng)域許多前沿物理課題以及相關(guān)工程技術(shù)的研究.

1.3磁約束聚變的目標和難點

磁約束聚變的目標是使粒子在擴散之前發(fā)生聚變，并讓聚變產(chǎn)生的alpha粒子加熱燃料，從而實現(xiàn)自持燃燒.實現(xiàn)這些目標需要對磁約束聚變發(fā)生的等離子體物理有深刻的理解，主要涉及等離子體的約束、輸運、加熱和電流控制等復雜物理過程.同時這些聚變等離子體物理需要研究波-粒子、波-波以及粒子-粒子等多時空尺度下復雜的非線性行為，需要設計復雜的計算機模擬來輔助裝置設計，物理實驗設計以及實驗解釋和預測.從20世紀90年代開始，磁約束聚變等離子體模擬就開始使用日益興起的超級計算技術(shù)進行大規(guī)模并行模擬，開發(fā)了針對聚變等離子體宏觀微觀不穩(wěn)定性、湍流輸運和加熱物理研究的大規(guī)模數(shù)值模擬程序，并興起了以回旋動理學模擬為主導的模擬科學.

2大規(guī)模并行數(shù)值模擬

2.1大規(guī)模并行數(shù)值模擬的起源

大規(guī)模并行數(shù)值模擬是指將一個復雜的計算問題分成若干部分，通過以網(wǎng)絡(光纖)相互聯(lián)接的大量的計算機核(幾百到20萬核)的超級計算機，同時計算同一個任務的不同部分.其目的是提高問題求解速度，同時擴大求解問題規(guī)模.

大規(guī)模并行數(shù)值模結(jié)合了最尖端的計算機科學、應用數(shù)學和具體科技應用領(lǐng)域的跨學科全新研究方法，適應于那些對科學研究和技術(shù)創(chuàng)新有重大革新性意義的課題，適合于團隊作業(yè).

2.2大規(guī)模并行數(shù)值模擬的發(fā)展

數(shù)值模擬已經(jīng)成為除理論、實驗外，支持科學探索和技術(shù)革新的第三個主要研究手段.大規(guī)模數(shù)值模擬能夠幫助理解和預測那些科學技術(shù)上極其重要的系統(tǒng)的行為，從而加速科技創(chuàng)新的速度.超級計算能力已成為衡量一個國家(地區(qū))科技實力和創(chuàng)新能力的重要標志.

2.3超級計算機的世界分布

大規(guī)模數(shù)值模擬從20世紀90年代中葉開始興起，最早的超級計算機計算能力約為1G flops，其計算能力以每年2～3倍的速度增加.現(xiàn)在世界上最強的超級算計機的速度在1015flop量級，正在向1018flop進軍.目前世界上最強的超級算計主要集中在美國、中國、歐洲和日本.其中美國的超級計算資源占到全世界一半左右.

1) 國際超級計算機

近幾年，世界上排名前500的超級計算機中，較著名的有美國橡樹嶺國家實驗室的“泰坦”，美國勞倫斯利福摩爾國家實驗室的“紅杉”，日本理化學研究所的“京”，美國阿貢國家實驗室的“米拉”等.

2) 國內(nèi)超級計算機

中國硬件發(fā)展上已經(jīng)成為僅次于美國的超級計算大國，深圳星云超級計算機，曾排名世界第二；天津天河超級計算機，2011排名世界第一；廣州天河-2超級計算機，2014排名世界第一.

現(xiàn)階段我國超級計算機發(fā)展中的主要問題是軟件發(fā)展和硬件發(fā)展不相稱.超級計算平臺輔助軟件和大規(guī)模數(shù)值模擬應用程序的發(fā)展嚴重滯后于超級計算機的建造.

2.4大規(guī)模數(shù)值模擬的應用

大規(guī)模數(shù)值模擬廣泛應用于新能源、新材料、物理、化工、生物、氣候環(huán)境、物流、金融等領(lǐng)域，對國民經(jīng)濟和科技創(chuàng)新產(chǎn)生廣泛而深刻的影響，如圖12.

圖12　大規(guī)模數(shù)值模擬應用舉例Figure 12　Application example of large-scale numerical simulation

3回旋動理學模擬

基于磁約束聚變等離子物理的高度復雜性和多時空尺度特性，回旋動理學模擬已成為研究聚變中不穩(wěn)定性和湍流輸運的主要工具.近幾十年來的研究表明，托卡馬克中普遍存在的電磁場微擾是造成托卡馬克反常輸運的主要原因，一般也稱之為湍流輸運.這些電磁場微擾是由各種宏觀和微觀的電磁不穩(wěn)定性激發(fā)的.聚變等離子體約束的好壞直接決定于其輸運過程，其產(chǎn)生和控制的物理機制是目前研究的前沿課題.

3.1等離子體的輸運過程

3.1.1輸運過程的種類

在等離子體中，輸運過程可以分為碰撞輸運與湍流輸運.

1) 碰撞輸運

等離子體中粒子之間會存在相互作用，一般稱為碰撞，可以用碰撞算符來表示.碰撞可以使粒子回旋運動的導心或者沿著磁力線運動的香蕉軌道中心產(chǎn)生偏移，從而導致碰撞輸運，如圖13.碰撞輸運即使在沒有電磁場擾動的情況下仍然存在，是托卡馬克裝置最低水平的輸運過程.

圖13　帶有相反電荷的粒子由于導心漂移發(fā)生碰撞Figure 13　collision of particles with opposite charge due to guide heart drift

2) 湍流輸運

湍流是托卡馬克中粒子和熱輸運的主要原因.湍流輸運可以分為兩種，靜電性和電磁性.根據(jù)其不穩(wěn)定性來源又可細分為各種不同模，比如離子梯度模(ITG)、電子梯度模(ETG)、捕獲電子模(TEM)、動理學氣球模(KBM)，等等.每個模線性不穩(wěn)定性原因不一樣，同時飽和機制和輸運機制也往往不一用，從而導致輸運系數(shù)不一樣.一種最常見湍流輸運原因就是帶電粒子的在靜電場微擾下的EXB漂移，這種垂直于磁場的漂移由于電場的退相關(guān)過程而導致隨機運動，從而導致徑向的湍流輸運，如圖14.

圖14　導心的EXB漂移導致湍流輸運Figure 14　Turbulent transport because of EXB drift of guide heart

3.2回旋動理學簡介

3.2.1Vlasov方程與Maxwell方程

等離子體有動理學和流體力學兩種描述方式，其中流體力學描述可通過對動理學方程求矩獲得，但是它忽略了朗道阻尼等重要動理學內(nèi)容，因而不能準確描述聚變等離子體中一些重要的動理學不穩(wěn)定性及其湍流輸運.最一般的動理學方程可以從非平衡態(tài)統(tǒng)計力學中的BBGKY遞推關(guān)系獲得如下的帶電磁場的玻爾茲曼方程

(1)

-▽2φ=4π(qi∫d3vfi-e∫d3vfe).

(2)

其中，f表示單粒子的相空間分布函數(shù).第一個方程在等離子體物理被稱之為含碰撞的Vlasov方程，方程右邊為代表兩粒子相關(guān)函數(shù)的碰撞算符.為簡單起見，我們用第二個方程即泊松方程代替求解電磁場的Maxwell(麥克斯韋)方程組.這組方程可以視為靜電極限下的Vlasov-Maxwell方程組.上面這組方程可以描述時間尺度從等離子體震蕩到輸運，空間尺度從德拜長度到系統(tǒng)尺寸等大范圍多尺度物理.依靠現(xiàn)在的超級計算機的計算能力是不可能模擬這么大范圍的多尺度物理.

但在托卡馬克中引起湍流輸運最危險的模是低頻漂移波湍流，比如ITG、TEM、KBM和ETG等，在燃燒等離子體中可能會激發(fā)比較高頻的Alfven波，但是仍然低于離子回旋頻率，如圖15.

圖15　托卡馬克中不穩(wěn)定性的時空頻率特點Figure 15　Space and time frequency characteristics of instability in Tokamak

(3)

在這種條件下，我們可以將式(3)的Vlasov方程進行坐標變換：

(x,v)→(R,μ,v||),

(4)

將方程從粒子坐標移到導心坐標，然后再回旋平均得到如下的回旋動理學(gyrokinetics)-麥克斯韋方程組[1-6]：

=C[f].

(5)

相應的回旋動理學泊松方程整理如下：

(6)

其中

代表從導心到粒子位置的反演[7].強磁場下的導心運動如圖16.

圖16　強磁場下的導心運動Figure 16　Movement of guide heart in strong magnetic field

3.3回旋動理學模擬

3.3.1模擬方法和算法特點

回旋動理學模擬是研究磁約束等離子體微觀不穩(wěn)定性中非線性物理的重要工具.它將六維相空間的動理學方程變成五維相空間，大大提高計算效率.相對流體模擬而言，它又保留了磁約束聚變等離子體物理中如朗道阻尼重要的動理學效應，可以比較真實地描述等離子體物理.等離子體物理模擬一般通過將模擬結(jié)果與理論分析進行對比來驗證程序真實反映了理論模型，再用模擬與實驗結(jié)果進行比較來驗證物理模型的可信度，最終達到模擬的可預測性.

回旋動理學湍流模擬的理論基礎(chǔ)起源于20世紀70年代，從90年代以來隨著超級計算機和大規(guī)模數(shù)值模擬的發(fā)展而獲得蓬勃發(fā)展.回旋動理學-麥克斯韋方程體系有兩種主要數(shù)值方法可求解其動力學行為，一種叫粒子模擬方法，即用蒙特卡洛的方法對相空間取樣；另外一種是用解多維偏微分方程的網(wǎng)格法，也叫Vlasov方法或者連續(xù)法.在算法上，為了去除模擬中的粒子數(shù)或網(wǎng)格數(shù)與物理要求相差很遠而導致的噪聲，我們一般采用所謂部分分布函數(shù)法，關(guān)注分布函數(shù)中偏離平衡分布的部分，從而提高信噪比.

3.3.2模擬程序

從20世紀90年代以來，涌現(xiàn)了很多回旋動理學程序.這些程序的發(fā)展代表了磁約束聚變模擬這二十年來最富代表性的發(fā)展歷程.代表性的程序美國有GTC、GYRO、GS2、GEM、XGC等，歐洲有GENE、ORB5等，日本有GT3D等.

其中，GTC(Gyrokinetic Toroidal Code)程序[8]是美國加州大學的林志宏教授90年代從普林斯頓讀研究生時開始發(fā)展的三維全域性粒子回旋動理學程序.歷經(jīng)二十年的發(fā)展，已成長為一個包含了微觀湍流、動理學磁流體、射頻加熱等多功能的模擬程序或者模擬平臺.近年來浙江大學、中科院物理所和北京大學也積極參與了其程序開發(fā).

GTC程序的創(chuàng)新之一體現(xiàn)在其采用磁坐標，構(gòu)造了沿著磁力線的三維電磁場網(wǎng)格，如圖17，圖18和圖19.在某一個極面上GTC采用非結(jié)構(gòu)化網(wǎng)格，保證徑向和極向計算精度大致相近；在同一個磁面上，格點沿著磁力線成為一條直線.沿著磁力線的格點設置契合磁約束聚變等離子體各項異性的特點.回旋平均和導心到粒子位置的反演都采用了所謂4點平均的辦法[9].這些計算算法的創(chuàng)新使得GTC能夠進行真正三維的實空間的全域模擬，同時大大提高了計算效率.

圖17　GTC磁面上的格點設置Figure 17　Lattice point set on magnetic surface of GTC

圖18　GTC極面上的格點設置Figure 18　Lattice point set on polar surface of GTC

圖19　GTC的3維格點設置和湍流特點Figure 19　3D point set and turbulence characteristics of GTC

3.3.3粒子推動方程

回旋動理學粒子模擬方法主要由兩部分組成，導心(或者回旋中心)運動積分器和電磁場解法器.導心的運動方程可以描述如下：

(7)

(8)

假設磁場擾動可以寫成如下形式：

δB=▽×αB,

(9)

ρc=ρ||+α,

(10)

D=ρc(gI′-Ig′)+I+qg,

(11)

磁坐標下的運動方程可以表述為[10]：

(12)

(13)

(14)

(15)

式(12)～(15)這組導心運動方程可以用來描述導心在托卡馬克中運動，一般用2階龍格庫塔方法來積分導心軌道，因為它是一種顯式算法，可以和描述碰撞的蒙特卡洛方法相匹配.

3.3.4泊松方程

(16)

式(16)的關(guān)鍵在于計算分子中的拉普拉斯算符(Laplacian).

3.3.5一般幾何位型下的拉普拉斯算符

磁坐標下的拉普拉斯算符(Laplacian)展開可得

(17)

其中，協(xié)變張量gxy可以通過B樣條方法計算得出.逆變張量gxy可以通過協(xié)變張量gxy計算得出.

3.4基于回旋動理學模擬的發(fā)現(xiàn)

回旋動理學模擬能夠幫助解釋實驗現(xiàn)象和促進理論的發(fā)展，并為實驗建模提供必需的參數(shù).

3.4.1帶狀流對湍流的抑制

回旋動理學模擬最有影響的結(jié)果之一是模擬發(fā)現(xiàn)湍流會自發(fā)激發(fā)帶狀流(zonal flow)，從而反過來抑制湍流[7]，如圖20.這一模擬上的發(fā)現(xiàn)極大促進人們對于湍流理論和實驗的理解，比如促進了漂移波參量不穩(wěn)定性理論的發(fā)展[11].我國的實驗也多次發(fā)現(xiàn)了帶狀流對湍流的調(diào)制作用[12-13].

3.4.2湍流輸運系數(shù)

關(guān)于磁約束聚變等離子體輸運過程的集成建模迫切需要在給定等離子體徑向剖面的條件下湍流輸運的強度，基于第一性原理的回旋動理學模擬能夠提供湍流輸運的系數(shù)[15]，如圖21.而對于湍流強度的徑向分布目前很難用現(xiàn)有理論預測.

圖20　GTC模擬發(fā)現(xiàn)捕獲電子模湍流中帶狀流對湍流有抑制作用[14]Figure 20　GTC simulation found that capture banded in electronic mode turbulence flow has inhibitory effect on turbulence

圖21　回旋動理學模擬提供的熱傳導系數(shù)的徑向分布Figure 21　Radial distribution of heat transport coefficient by gyrokinetic simulation

4結(jié)語

本文介紹了磁約束核聚變中的基于回旋動理學的大規(guī)模數(shù)值模擬.最近20年來這個領(lǐng)域已經(jīng)取得很大進步，但是還有很多關(guān)鍵性問題沒有得到解決，比如聚變等離子體邊界電磁模擬和長時間尺度模擬，仍需要研究人員的持續(xù)努力.同時，這個領(lǐng)域隨著超級計算技術(shù)的進步，將會進一步促進可控磁約束聚變的實現(xiàn).

(致謝：本文作者感謝袁泉準備部分初稿并仔細校對文稿，同時感謝P.J. Catto和Z. Lin博士提供部分插圖.)

【參考文獻】

[1]CATTO P J. Linear gyrokinetics[J]. Plasma Physics,1978,20:719.

[2]ANTONSEN T, LANE B. Kinetic equations for low frequency instabilites in inhomogeneous plasmas[J].Physics of Fluids,1980,23(6):061205.

[3]FRIEMAN E A, CHEN L. Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria[J]. Physics of Fluidsvol,1982,25(3):502-508.

[4]HAHM T S. Nonlinear gyrokinetic equations for tokamak microturbulence[J]. Physics of Fluids,1988,31(9):2670-2673.

[5]QIN H W, TANG M, REWOLDT G. Gyrokinetic theory for arbitrary wavelength electromagnetic modes in tokamaks[J]. Physics of Plasmas,1998,5(4):1035-1049.

[6]LEE W W. Nonlinear gyrokinetic equations[J]. Physics of Fluids,1983,26(12):3524.

[7]LIN Zhihong, HAHM T S, LEE W W, et al.Turbulent transport reduction by zonal flows: Massively parallel simulations[J]. Science,1998,281(5384):1835-1837.

[8]LEE W W. Gyrokinetic particle simulation model[J]. Journal of Computational Physics,1987,72(1):243-269.

[9]LIN Z, LEE W W. Method for solving the gyrokinetic Poisson equation in general geometry[J]. Phys Rev E,1995,52(5):5646-5652.

[10]WHITE R B, CHANCE M S. Hamiltonian guiding center drift orbit calculation for plasmas of arbitrary cross-section[J]. Physics of Fluids,1984,27(10):2455-2467.

[11]CHEN L, LIN Z, WHITE R B. Excitation of zonal flow by drift waves in toroidal plasmas[J].Physics of Plasmas,2000,7(8):3129.

[12]ZHAO K J, LAN T, DONG J Q, et al. Toroidal symmetry of the geodesic acoustic mode zonal flow in a tokamak plasma[J]. Physical Review Letters,2006,96(25):255004.

[13]XU G, WAN B N, SONG M, et al. Direct measurement of poloidal long-wavelength E× B flows in the HT-7 tokamak[J]. Physical Review Letters,2003,91(12):125001.

[14]XIAO Y, LIN Z. Turbulent transport of trapped-electron modes in collisionless plasmas[J]. Physical Review Letters, 2009,103(8):085004.

[15]XIAO Y, HOLOD I, WANG Z X, et al. Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles[J]. Physics of Plasmas,2015(22):022516.

Massively parallel gyrokinetic simulation for magnetic fusion plasmas

XIAO Yong

(Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, China)

Abstract:This article introduces the magnetic fusion energy science and its associated fusion plasma physics. Massively parallel gyrokinetic simulation has emerged as a major tool to study the nonlinear physics. Its research frontier is also discussed.

Key words:plasma; fusion; magnetic confinement; large scale data simulation; gyrokinetics

【文章編號】1004-1540(2016)01-0001-10

DOI:10.3969/j.issn.1004-1540.2016.01.001

【收稿日期】2016-01-18《中國計量學院學報》網(wǎng)址：zgjl.cbpt.cnki.net

【作者簡介】肖湧(1977- )，男，湘南省湘潭人，教授，主要研究方向為磁約束核聚變等離子體中的不穩(wěn)定性，輸運和加熱以及大規(guī)模并行回旋動理學模擬等. E-mail:yxiao@zju.edu.cn

【中圖分類號】O532+.11

【文獻標志碼】A