查學(xué)軍, 鐘德俊, 王福瓊, 陳一平, 盧洪偉, 胡立群（.東華大學(xué)理學(xué)院應(yīng)用物理系,上海 060；.中國(guó)科學(xué)院等離子體物理研究所,合肥 3003）

EAST托卡馬克雜質(zhì)輸運(yùn)蒙特卡洛模擬

查學(xué)軍1, 鐘德俊1, 王福瓊2, 陳一平2, 盧洪偉1, 胡立群2
（1.東華大學(xué)理學(xué)院應(yīng)用物理系,上海 201620；2.中國(guó)科學(xué)院等離子體物理研究所,合肥 230031）

介紹DIVIMP程序的基本思想和采用的物理模型,研究EAST托卡馬克歐姆放電條件下碳雜質(zhì)的產(chǎn)生與輸運(yùn).模擬結(jié)果與實(shí)驗(yàn)測(cè)量具有較好的一致性,說(shuō)明將該程序應(yīng)用于EAST裝置中雜質(zhì)輸運(yùn)模擬的合理性.并利用DIVIMP程序預(yù)測(cè)EAST在8MW加熱功率條件下鎢雜質(zhì)的輸運(yùn)特性.

托卡馬克；雜質(zhì)輸運(yùn)；模擬；蒙特卡洛方法；DIVIMP

0 引言

托卡馬克[1]是一種利用磁約束來(lái)實(shí)現(xiàn)受控?zé)岷司圩兊沫h(huán)形裝置.托卡馬克邊界等離子體與第一壁材料發(fā)生相互作用,所產(chǎn)生的雜質(zhì)能通過(guò)邊界等離子體的輸運(yùn)進(jìn)入主約束區(qū)等離子體,并對(duì)其性能造成極大影響.一方面,輸運(yùn)進(jìn)入聚變反應(yīng)堆芯部的雜質(zhì)通過(guò)稀釋燃料濃度來(lái)降低聚變功率,同時(shí)通過(guò)雜質(zhì)輻射降低等離子體溫度,縮小等離子體電流通道,從而引發(fā)不穩(wěn)定性,影響等離子體約束；另一方面,存在于托卡馬克裝置偏濾器區(qū)域的雜質(zhì),其輻射過(guò)程有利于降低到達(dá)偏濾器靶板的能流,進(jìn)而又會(huì)減小偏濾器靶板材料的腐蝕及雜質(zhì)的產(chǎn)生.因此,雜質(zhì)輸運(yùn)一直是磁約束等離子體最重要的研究課題之一[2].托卡馬克裝置中雜質(zhì)產(chǎn)生及輸運(yùn)特性的研究以及雜質(zhì)控制是磁約束受控?zé)岷司圩兡茉瓷逃没M(jìn)程中必不可少的環(huán)節(jié).大科學(xué)工程裝置EAST是由中國(guó)科學(xué)院等離子體物理研究所建立的世界上第一個(gè)具有非圓截面全超導(dǎo)托卡馬克實(shí)驗(yàn)裝置.EAST邊界等離子體雜質(zhì)輸運(yùn)研究的目的是為了合理設(shè)計(jì)并優(yōu)化第一壁和偏濾器靶板部件,更好地控制等離子體中的雜質(zhì).

實(shí)驗(yàn)診斷和程序模擬是研究托卡馬克裝置中雜質(zhì)特性的兩種主要方法[3].雜質(zhì)輸運(yùn)與邊界等離子體輸運(yùn)強(qiáng)烈的相互耦合：等離子體施加給雜質(zhì)粒子的力控制了雜質(zhì)的輸運(yùn)過(guò)程,而雜質(zhì)輻射又會(huì)影響等離子體參數(shù)的分布.雜質(zhì)與等離子體之間的這種相互作用最終決定了雜質(zhì)密度剖面.由于,大多數(shù)雜質(zhì)相關(guān)實(shí)驗(yàn)診斷,如光譜儀和熱輻射儀等,僅測(cè)量了沿著觀測(cè)線的積分信號(hào),測(cè)量過(guò)程掩蓋了許多雜質(zhì)分布信息.刮削層中的等離子體溫度及流速場(chǎng)分布是決定雜質(zhì)輸運(yùn)的重要因素,但是目前大多數(shù)托卡馬克裝置中幾乎沒(méi)有有關(guān)這些物理量在整個(gè)刮削層中完整分布的實(shí)驗(yàn)測(cè)量.雜質(zhì)分布及等離子體參數(shù)剖面信息的缺失,必須通過(guò)模擬手段予以彌補(bǔ).

當(dāng)然,近年來(lái)雜質(zhì)相關(guān)實(shí)驗(yàn)診斷也得到了較好的發(fā)展.測(cè)量等離子體流速的馬赫探針[4－5]以及測(cè)量偏濾器區(qū)域等離子體溫度、密度的湯姆遜散射系統(tǒng)[6]相繼在各托卡馬克裝置中成功開(kāi)展.在許多托卡馬克中,已經(jīng)通過(guò)多通道熱輻射測(cè)量?jī)x（Bolometer）信號(hào)[7]成功建立了總輻射的二維（如圖1所示的二維極平面）剖面,并通過(guò)光譜診斷[8]建立了雜質(zhì)輻射的二維剖面.此外,在DIII-D等托卡馬克裝置中,已經(jīng)能夠通過(guò)大量的偏濾器測(cè)量推斷雜質(zhì)密度在極向平面的二維分布[9].上述各二維剖面信息均將為我們提供理解雜質(zhì)輸運(yùn)的有效信息.但是,一些診斷技術(shù)的使用受到物理或工程限制,如Bolometer所采用的斷層分析算法具有很大的不確定性,且僅能針對(duì)一些特定的電荷態(tài).此外,目前邊界等離子體溫度測(cè)量主要限于測(cè)量電子溫度,對(duì)離子溫度的直接測(cè)量結(jié)果很少,而離子溫度梯度力通常又主導(dǎo)了雜質(zhì)離子向芯部的輸運(yùn)過(guò)程.因此,為了解釋實(shí)驗(yàn)信息,詳細(xì)理解托卡馬克邊界物理活動(dòng),程序模擬包括對(duì)雜質(zhì)輸運(yùn)特性的模擬仍然必不可少.

綜上所述,程序模擬不僅可以克服雜質(zhì)相關(guān)實(shí)驗(yàn)診斷信息和診斷區(qū)域有限的缺點(diǎn),還是解釋實(shí)驗(yàn)診斷雜質(zhì)信息的有效工具.另外,實(shí)驗(yàn)與模擬結(jié)果的對(duì)比分析,可以幫助證實(shí)模擬中所采用模型的合理性.雜質(zhì)輸運(yùn)程序主要有流體程序和蒙特卡洛程序兩大類[10].蒙特卡洛雜質(zhì)輸運(yùn)程序可以方便地模擬發(fā)生在托卡馬克邊界的復(fù)雜分子、原子過(guò)程,不需要對(duì)托卡馬克位形作出除環(huán)向?qū)ΨQ以外的任何假定和近似.同時(shí),由于雜質(zhì)粒子尤其是低價(jià)態(tài)雜質(zhì)離子的特征電離時(shí)間往往短于粒子之間的平均碰撞自由時(shí)間,流體近似條件[11]可能得不到很好的滿足.此外,對(duì)于一些核電荷數(shù)較高的雜質(zhì)粒子（如74鎢）的輸運(yùn)過(guò)程,如采用流體程序模擬,需要考慮的離子種類及需要求解的方程總數(shù)較多,計(jì)算時(shí)間往往較長(zhǎng),而抽樣跟蹤單粒子運(yùn)動(dòng)軌道及電荷態(tài)的蒙特卡洛雜質(zhì)輸運(yùn)程序可克服該缺點(diǎn).因此,蒙特卡洛程序在模擬雜質(zhì)輸運(yùn)特性方面,有其優(yōu)越性.

實(shí)驗(yàn)先進(jìn)超導(dǎo)托卡馬克（EAST）[12]于2006首次投入運(yùn)行,并于2010年成功實(shí)現(xiàn)高約束運(yùn)行模式.目前,該裝置尚處于運(yùn)行初期階段,一些診斷手段,特別是邊界過(guò)程相關(guān)診斷,還不完善,如邊界湯姆遜系統(tǒng)、測(cè)量等離子體密度剖面的束發(fā)射光譜診斷等尚未成功在EAST裝置中予以開(kāi)展.同時(shí),盡管已成功開(kāi)展一些光譜診斷,如Da信號(hào)及CIII信號(hào)的診斷等,但雜質(zhì)相關(guān)診斷信息仍非常有限,并且沒(méi)有邊界等離子體參數(shù)二維分布信息.目前,CIII信號(hào)觀測(cè)通道僅涵蓋了外偏濾器區(qū)域,且沒(méi)有除CIII以外的其它C雜質(zhì)線輻射信號(hào)（如CII輻射等）的測(cè)量.此外,現(xiàn)有的實(shí)驗(yàn)測(cè)量也難免會(huì)有一些誤差.因此,針對(duì)EAST裝置開(kāi)展有關(guān)雜質(zhì)輸運(yùn)特性的模擬研究,不僅可在實(shí)驗(yàn)條件尚不完善的條件下,理解該裝置雜質(zhì)輸運(yùn)基本特性,還可通過(guò)模擬結(jié)果與現(xiàn)有實(shí)驗(yàn)數(shù)據(jù)之間的比較,檢驗(yàn)?zāi)M所采用的模型,為模擬工作的進(jìn)一步開(kāi)展奠定基礎(chǔ).

本文運(yùn)用DIVIMP程序[13]模擬研究EAST托卡馬克在歐姆加熱情況下碳雜質(zhì)的輸運(yùn)特性,并得到其相應(yīng)的碳雜質(zhì)二維空間分布.基于碳雜質(zhì)密度分布信息和背景等離子的溫度密度信息,計(jì)算碳雜質(zhì)輻射分布.實(shí)驗(yàn)測(cè)量的CIII信號(hào)與相關(guān)模擬結(jié)果之間有很好的一致性,說(shuō)明將DIVIMP用于EAST裝置上對(duì)雜質(zhì)輸運(yùn)特性的解釋和預(yù)測(cè)的合理性.由于鎢材料具有熔點(diǎn)高、耐腐蝕性強(qiáng)等優(yōu)點(diǎn),國(guó)際熱核聚變實(shí)驗(yàn)堆（ITER）、聚變示范堆（DEMO）等運(yùn)行參數(shù)更高的下一代聚變裝置的第一壁材料均考慮采用鎢材料.EAST托卡馬克偏濾器材料也將逐步升級(jí)改造為鎢材料.于是,采用DIVIMP程序預(yù)測(cè)了EAST高功率加熱條件下鎢雜質(zhì)的輸運(yùn)特性.模擬結(jié)果可為EAST第一壁和偏濾器的優(yōu)化升級(jí)和實(shí)驗(yàn)的優(yōu)化設(shè)計(jì)提供參考[14].

1 模型與方法

DIVIMP是通過(guò)抽樣跟蹤偏濾器位形下托卡馬克邊界等離子體中雜質(zhì)離子狀態(tài)和軌道的蒙特卡洛程序.與傳統(tǒng)的流體程序相比,DIVIMP考慮的是單個(gè)粒子的輸運(yùn),模擬它們?cè)谖镔|(zhì)中運(yùn)動(dòng)的歷史,記錄其在運(yùn)動(dòng)中相關(guān)物理模擬量的貢獻(xiàn),在對(duì)單個(gè)粒子運(yùn)動(dòng)歷史進(jìn)行大量的重復(fù)模擬之后,再對(duì)物理模擬量進(jìn)行統(tǒng)計(jì)平均,得到所需要的物理結(jié)果.它與傳統(tǒng)的解法相比,程序結(jié)構(gòu)簡(jiǎn)單,易于實(shí)現(xiàn),誤差較小.以下是對(duì)DIVIMP中雜質(zhì)產(chǎn)生和輸運(yùn)計(jì)算模型的詳細(xì)描述.

圖1 單零偏濾器位形Fig.1 Skematic of a single null divertor configuration

1.1 雜質(zhì)的產(chǎn)生

圖1是單零偏濾器位形下,托卡馬克磁約束等離子體的二維剖面.附加環(huán)向（圖1中垂直紙面方向）導(dǎo)體中有與等離子體電流（IP）同向的電流（ID）,由ID和IP所產(chǎn)生的極向磁場(chǎng)零點(diǎn)稱為X點(diǎn),經(jīng)過(guò)X點(diǎn)的磁面為分界面.在分界面外,磁力線開(kāi)放的區(qū)域?yàn)楣蜗鲗樱⊿crape-off Layer,SOL）.刮削層中的開(kāi)放磁力線與偏濾器靶板直接相交,如圖2所示[15].當(dāng)背景離子沿著磁力線運(yùn)動(dòng)到達(dá)偏濾器靶板表面時(shí),靶板材料會(huì)因?yàn)闉R射等過(guò)程向等離子體中發(fā)射雜質(zhì).同樣,返回偏濾器表面的雜質(zhì)離子也可能引發(fā)自濺射.

濺射過(guò)程可用濺射產(chǎn)額Y來(lái)描述,其定義為平均每入射一個(gè)粒子從靶板表面濺射出來(lái)的原子數(shù),即

圖2 偏濾器區(qū)域雜質(zhì)產(chǎn)生和輸運(yùn)過(guò)程Fig.2 Production and transport processes of impurities in the divertor region

模擬中,濺射產(chǎn)額Y對(duì)入射粒子能量的依賴采用2007年由Eckstein等人所推薦的關(guān)系式[16].假定背景離子的入射能量Eimp為

其中,ZB為等離子體離子的電荷數(shù),Te0和Ti0分別為偏濾器靶板處的電子和離子溫度.由背景等離子體中離子的轟擊所引起雜質(zhì)粒子的總濺射通量為

1.2 雜質(zhì)粒子的輸運(yùn)

1.2.1 中性雜質(zhì)粒子的輸運(yùn)

等離子體與偏濾器靶板相互作用產(chǎn)生的雜質(zhì)粒子均以中性粒子的形式進(jìn)入等離子體中.與雜質(zhì)離子不同,中性雜質(zhì)粒子的運(yùn)動(dòng)不受電磁場(chǎng)約束,在DIVIMP模擬中分別由不同的程序模塊對(duì)二者的運(yùn)動(dòng)狀態(tài)進(jìn)行單獨(dú)跟蹤.在t＝0時(shí)刻,中性雜質(zhì)粒子以上述計(jì)算所得概率從偏濾器部件發(fā)射,并遵從湯普森能量分布和余弦角分布.在中性粒子被電離、或打擊到固體表面之前,他們保持直線運(yùn)動(dòng).模擬中假定,當(dāng)中性粒子打擊到偏濾器靶板表面,其被反射,而當(dāng)中性粒子打擊到除靶板以外的第一壁表面則被吸收.在時(shí)間步長(zhǎng)Δtn時(shí)間內(nèi),中性粒子發(fā)生電離的概率

1.2.2 雜質(zhì)離子的輸運(yùn)

雜質(zhì)離子產(chǎn)生后,可能經(jīng)歷電離、復(fù)合、電荷交換等原子、分子過(guò)程.定義：

雜質(zhì)離子沿著磁力線方向的運(yùn)動(dòng)采用經(jīng)典理論描述.模擬中,假定離子沿磁力線方向受力情況為

其中,等式右邊的五項(xiàng)分別為雜質(zhì)壓力梯度力、與背景等離子體之間的摩擦力、電場(chǎng)力、電子溫度梯度力、離子溫度梯度力等,s為沿磁力線的距離,p為雜質(zhì)離子的壓強(qiáng),VB為背景等離子體的流速,v為雜質(zhì)離子運(yùn)動(dòng)速度,Z為雜質(zhì)離子的電荷數(shù),αe和βi分別為電子和離子溫度梯度系數(shù).從方程（6）可見(jiàn),雜質(zhì)離子的輸運(yùn)與背景等離子體流速場(chǎng)和電場(chǎng)、電子溫度梯度、離子溫度梯度等物理量有關(guān).模擬中,這些背景等離子體信息由SOLPS（B2－EIRENE）[17]計(jì)算得到.假定雜質(zhì)離子在垂直磁力線方向上的運(yùn)動(dòng)為反常擴(kuò)散.反常擴(kuò)散系數(shù)一般不能由第一性原理導(dǎo)出.根據(jù)相關(guān)參考文獻(xiàn)和國(guó)際上其它托卡馬克中雜質(zhì)輸運(yùn)特性模擬所采用的相關(guān)假定[18],模擬中規(guī)定雜質(zhì)離子垂直磁力線反常擴(kuò)散系數(shù)D＝1.0 m2·s－1,并且不考慮雜質(zhì)離子的徑向內(nèi)箍縮效應(yīng).在經(jīng)歷每個(gè)時(shí)間步長(zhǎng)Δt后,DIVIMP分別記錄每個(gè)雜質(zhì)離子的電荷態(tài)、位置、速度等信息.達(dá)到穩(wěn)態(tài)以后,計(jì)算得到雜質(zhì)密度和雜質(zhì)輻射等信息.

1.3 雜質(zhì)密度

圖3 DIVIMP的計(jì)算網(wǎng)格Fig.3 Grid for DIVIMPmodeling

DIVIMP計(jì)算將邊界等離子體區(qū)域劃分為計(jì)算網(wǎng)格,如圖3所示,對(duì)雜質(zhì)、粒子運(yùn)動(dòng)狀態(tài)的跟蹤和雜質(zhì)相關(guān)信息的計(jì)算基于網(wǎng)格進(jìn)行.網(wǎng)格基于實(shí)際托卡馬克放電時(shí)的平衡位形.可以看出,DIVIMP計(jì)算區(qū)域包括整個(gè)刮削層區(qū)域和分界面以內(nèi)的芯部外圍區(qū)域.

每一個(gè)網(wǎng)格由沿著磁力線的x坐標(biāo)和垂直磁力線的y坐標(biāo)來(lái)描述.t＝0時(shí)刻發(fā)射的N個(gè)雜質(zhì)離子,經(jīng)過(guò)ts后位于網(wǎng)格（x,y）內(nèi)并具有ZI電荷態(tài)雜質(zhì)離子歸一化密度為

當(dāng)發(fā)射的第i個(gè)離子進(jìn)入（x,y）網(wǎng)格并處于ZI電荷態(tài)時(shí),Ci（x,y, ZI,t）＝1（對(duì)于自濺射Ci（x,y,ZI,t）＜1）；否則,Ci（x,y,ZI,t）＝0.達(dá)到穩(wěn)態(tài)時(shí),（x,y）網(wǎng)格內(nèi)的電荷態(tài)ZI雜質(zhì)離子的歸一化密度為

因此,（x,y）網(wǎng)格內(nèi)ZI電荷態(tài)雜質(zhì)離子的絕對(duì)密度為（x,y,Z1）R,考慮到部分中性粒子直接沉積到壁表面而對(duì)離子數(shù)沒(méi)有貢獻(xiàn),R為電離產(chǎn)生離子的中性粒子份額.

1.4 雜質(zhì)輻射功率

DIVIMP中一項(xiàng)重要的輸出信息是,單位時(shí)間內(nèi)由某一電荷態(tài)雜質(zhì)離子在等離子體中給定位置附近單位體積內(nèi)所發(fā)射的具有特定波長(zhǎng)（頻率）光子數(shù)ΦPLRP.根據(jù)ΦPLRP信息和ADAS數(shù)據(jù)庫(kù)中有關(guān)一個(gè)特定電荷態(tài)雜質(zhì)離子在其被電離到其它電荷態(tài)之前所能釋放的光子數(shù)PE（光子發(fā)射效率）,可以推斷該電荷態(tài)雜質(zhì)離子的產(chǎn)生速率

1）電子碰撞引起的雜質(zhì)離子從基態(tài)的躍遷.這是最基本的過(guò)程；

2）電子－離子復(fù)合,該過(guò)程也可能導(dǎo)致激發(fā)態(tài)的產(chǎn)生,從而會(huì)釋放光子；

3）從上一電荷態(tài)的電離,形成激發(fā)態(tài)并伴隨光子發(fā)射；

4）電荷交換復(fù)合,此時(shí)需要輸入二維中性氫粒子分布相關(guān)信息（如來(lái)源于ERIENE的計(jì)算等）.

在實(shí)際計(jì)算中,可根據(jù)求解問(wèn)題的情況僅僅考慮上述過(guò)程中的某些過(guò)程.如只考慮過(guò)程1）,則

其中η為發(fā)射一個(gè)光子所發(fā)生的平均電離事件次數(shù),Iiz（x,y,z1）為單位時(shí)間發(fā)生電離的雜質(zhì)離子數(shù).

其中τiz為雜質(zhì)離子的電離特征時(shí)間.

本文首先采用上述介紹的計(jì)算模型,模擬了EAST托卡馬克第12303炮放電,歐姆加熱條件下C雜質(zhì)的輸運(yùn)特性.與模擬相關(guān)的主要參數(shù)如表1所示.表1中,安全因子定義為磁面上所躺磁力線在極向轉(zhuǎn)一圈后,沿著環(huán)向所繞行的圈數(shù),q95為95％磁通面上的安全因子,環(huán)向比壓βT＝〈P〉/[/（2μ0],極向比壓βP＝〈P〉/[/（2μ0）],BT和BP分別為環(huán)向磁場(chǎng)和極向磁場(chǎng),〈P〉為等離子體體積平均熱壓力,線平均密度為等離子體密度在沿R方向過(guò)磁軸的直線上的平均值.其它等離子體參數(shù)的定義如圖4所示.然后,本文采用DIVIMP程序預(yù)測(cè)了EAST裝置在8MW高功率加熱條件下鎢雜質(zhì)的輸運(yùn)特性.

表1 與模擬研究相關(guān)的主要參數(shù)Table 1 M ain parameters used in simulation

圖4 托卡馬克磁面位形極向剖面Fig.4 Poloidal cross-section ofmagnetic surfaces in a tokamak

2 結(jié)果與討論

圖5是模擬得到的碳雜質(zhì)密度在邊界區(qū)域的二維分布.由圖5可見(jiàn),雜質(zhì)離子的二維分布沿著磁力線方向和垂直磁力線方向均存在著不均勻性.雜質(zhì)密度垂直磁力線方向（徑向）不均勻性可由Engelhard模型定性解釋.而雜質(zhì)密度沿著磁力線的不對(duì)稱性主要由存在于邊界區(qū)域中的許多復(fù)雜雜質(zhì)過(guò)程,如電離、復(fù)合等所導(dǎo)致.另外,由圖5可見(jiàn),歐姆加熱條件下由于加熱功率非常低,即使對(duì)于較容易被腐蝕的碳雜質(zhì)在芯部等離子體中的最高密度約為1015量級(jí).

如1.4節(jié)所述.基于DIVIMP計(jì)算得到的雜質(zhì)密度信息和SOLPS所提供的等離子體參數(shù)信息,DIVIMP可輸出由C＋2離子輻射產(chǎn)生的CIII信號(hào)（波長(zhǎng)465.0nm）的二維分布.圖6是CIII輻射信號(hào)沿著經(jīng)過(guò)偏濾器區(qū)域的18個(gè)觀測(cè)通道的積分信號(hào),橫軸1－18分別對(duì)應(yīng)于不同的觀測(cè)通道.值得注意的是,圖6中的實(shí)驗(yàn)測(cè)量值為未絕對(duì)標(biāo)定的相對(duì)值,而計(jì)算值為實(shí)際值.由圖6可見(jiàn),模擬結(jié)果與實(shí)驗(yàn)結(jié)果有較好的符合,說(shuō)明了將DIVIMP程序用于EAST裝置雜質(zhì)輸運(yùn)特性模擬的合理性.

圖5 DIVIMP計(jì)算所得碳雜質(zhì)密度在邊界的二維分布Fig.5 Two-dimensional distribution of C impurities in the edge of EAST from DIVIMP calculation

圖6 CIII信號(hào)的計(jì)算值和測(cè)量值對(duì)比Fig.6 Comparison against calculated CIIIsignals and measured ones

利用DIVIMP程序預(yù)測(cè)了未來(lái)EAST裝置在8MW加熱功率條件下鎢雜質(zhì)的輸運(yùn)特性,結(jié)果如圖7所示.在8MW的高功率加熱條件下,芯部等離子體中鎢雜質(zhì)的最高密度約為1015量級(jí).盡管隨著加熱功率的升高,偏離靶板處的等離子體溫度及離子轟擊靶板時(shí)的入射能量也將升高,但具有最高束縛能的鎢原子（8.68eV）,其濺射產(chǎn)額沒(méi)有明顯的增加.此外,由于鎢雜質(zhì)粒子的質(zhì)量較大,其從偏濾器靶板濺射的出射速度較低,而其特征電離時(shí)間又非常短,也就是說(shuō),濺射產(chǎn)生的鎢雜質(zhì)原子,大多數(shù)將在偏濾器靶板附近電離,存在于偏濾器區(qū)域的高強(qiáng)粒子流的沖刷,以及靶板表面鞘層電壓對(duì)鎢離子的加速作用,使得這些鎢離子重新返回偏濾器靶板[19].上述各原因使得,即使當(dāng)加熱功率大幅提高（～50倍）時(shí),如將偏濾器靶板材料由碳材料改變?yōu)殒u材料,等離子體中雜質(zhì)水平不會(huì)有太明顯的提高.

圖7 DIVIMP計(jì)算所得鎢雜質(zhì)在邊界層的密度分布Fig.7 Two-dimensional distribution ofW impurities in the edge of EAST from DIVIMP calculation

3 結(jié)論

用DIVIMP程序模擬了EAST歐姆加熱條件下的邊界層碳雜質(zhì)的輸運(yùn),得到了相對(duì)應(yīng)的雜質(zhì)分布,并將計(jì)算得到的雜質(zhì)輻射信息與實(shí)驗(yàn)測(cè)量信息進(jìn)行比較.結(jié)果表明：（1）歐姆加熱條件下,EAST裝置中碳雜質(zhì)的最高水平約為1015量級(jí).（2）模擬結(jié)果與實(shí)驗(yàn)結(jié)果之間符合較好,證明將蒙特卡洛程序DIVIMP用于研究EAST裝置雜質(zhì)輸運(yùn)特性模擬的合理性.采用DIVIMP程序模擬預(yù)測(cè)了未來(lái)EAST在8 MW高功率加熱條件下鎢雜質(zhì)的輸運(yùn)特性.結(jié)果發(fā)現(xiàn),高加熱功率條件下,如將偏濾器靶板由碳材料改為鎢材料,雜質(zhì)含量不會(huì)有明顯的提高.

[1] 石秉仁.磁約束聚變?cè)砼c實(shí)踐[M].北京：原子能出版社,1999.

[2] 徐偉,萬(wàn)寶年,謝紀(jì)康.HT-6M托卡馬克裝置雜質(zhì)輸運(yùn)[J].物理學(xué)報(bào),2003,52（8）：1970－1909.

[3] Shimizu K,Takizuka T,Sakasai A.A review on impurity transport in divertors[J].Journal of Nuclear Materials,1997,241－243：167－181.

[4] Stansfield B L,Boucher C,Gunn J P,et al.Measured plasma parameters in TdeV's closed poloidal divertors[J].JNucl Mater,1995,220－222：1121－1125.

[5] LaBombard B,Jablonski D,Lipschultz B,McCracken G,Goetz J.Scaling of plasma parameters in the SOL and divertor for Alcator C-Mod[J].JNucl Mater,1995,220－222：976－981.

[6] Allen SL,Brooks N H,Campbell R B.Developmentof a radiative divertor for DIII-D[J].JNuclMater,1995,220－222：336 －341.

[7] Reichle R,Clement S,Gottardi N.Controlled fusion and plasma physics[C]∥Proc 22nd EPS Conf,The European Physical Society,Bournemouth,UK,July 3-July 7,1995,19C（Ill）：85.

[8] Allen S L,Bozek AS,Brooks NH.Recent DIII-D divertor research[J].Plasma Phys Control Fusion,1995,A191.

[9] Fenstermacher M E,Porter GD,Rensink M E,etal.UEDGE and DEGASmodeling of the DIII-D scrape-off layer plasma[J].JNucl Mater,1995,220－222：330－335.

[10] Stangeby PC.Dynamics of impurities in the scrape-off layer[J].Contrib Plasma Phys,1988 28：507.

[11] Wesson J.Tokamaks[M].International Series of Monographs on Physics,Oxford University Press,2011.

[12] Wang FQ,Chen Y P,Hu LQ,et al.Onion skinmodel（OSM）analysis of EAST SOL plasmas[J].Nucl Fusion,2014,54：093002..

[13] Stangeby PC,Elder JD.Calculation of observable quantities using a divertor impurity interpretive code,DIVIMP[J].JNucl Mater,1992,258：196－198.

[14] 郭后揚(yáng).JET中偏濾器幾何結(jié)構(gòu)和化學(xué)濺射對(duì)雜質(zhì)行為和等離子體性能的影響[J].國(guó)外核聚變,2001,4：48－63.

[15] Stangeby PC.The plasma boundary ofmagnetic fusion devices[M].Bristol：Institute of Physics Publishing,2000.

[16] Behrisch R,Eckstein W.Sputtering by particle bombardment：Experiments and computer calculations from threshold to MeV energies[M].1st ed.Topics in Applied Physics.Springer,Softcover Reprint of Hardcover,2007.

[17] Schneider R,Bonnin X,Borrass K.Plasma edge physicswith B2-eirene[J].Contrib Plasma Phys,2006,46（1/2）：3－191.

[18] Schmid K,Krieger K,Kukushkin K,Loarte A.DIVIMPmodeling of tungsten impurity transport in ITER[J].JNucl Mater, 2007,363－365：674－679.

[19] J?rvinen A,Giroud C,Groth M,Krieger K,Moulton D,Wiesen S,Brezinsek S,JET-EFDA Contributors.DIVIMP simulation ofW transport in the SOL of JET H-mode plasmas[J].Phys Scr,2011,T145：014013.

M onte-Carlo M odelling of Im purity Transport in EAST Tokamak

ZHA Xuejun1, ZHONG Dejun1, WANG Fuqiong2, CHEN Yiping2, LU Hongwei1, HU Liqun2
（1.Department ofApplied Physics,Donghua University,Shanghai 201620,China；
2.Institute of Plasma Physics,Chinese Academy ofSciences,Hefei 230031,China）

Basic idea and physicalmodels in DIVIMP are introduced and DIVIMPmodelling of C impurity production and transport properties in edge plasma of EAST machine with ohmic heating are provided.Simulation results show favorable consistency with experimentalmeasurements,demonstrating that themodels coupled in DIVIMP are suitable for interpretive and predictivemodelling of impurity transport in EASTmachine.Transport properties ofW impurity in EASTwith heating power as high as8MW are also predicted by DIVIMPmodelling.

tokamak；impurity transport；simulation；Monte Carlomethod,DIVIMP

1001-246X（2015）06-0715-07

TL61＋2

A

2014－11－09；

2015－04－09

國(guó)家自然科學(xué)基金（11175045）和科技部ITER計(jì)劃專項(xiàng)課題（2013GB109001）資助項(xiàng)目

查學(xué)軍（1971－）,男,博士,副研究員,主要從事等離子體物理理論和數(shù)值模擬研究,E-mail：xjzha＠dhu.edu.cn