本文關(guān)鍵詞： 數(shù)值模擬 邊界局域模熱流 鎢偏濾器靶板 熔化 流動(dòng)　出處：《大連理工大學(xué)》2016年博士論文　論文類型：學(xué)位論文

【摘要】：等離子體面壁材料的選擇是聚變裝置托卡馬克的關(guān)鍵問題,但到目前為止,還沒有真正合適的材料。鎢材料因?yàn)榫哂袩釋?dǎo)率高、熔點(diǎn)高、濺射率低、氚滯留率低等優(yōu)點(diǎn),成為首選的面壁材料。但是鎢材料在高瞬時(shí)熱流的作用下會(huì)熔化、汽化,熔化層在各種力(如表面張力、磁場(chǎng)力、壓強(qiáng)梯度力等)作用下的流動(dòng),導(dǎo)致材料形貌變化,影響材料的傳熱能力,誘發(fā)電弧,甚至液體鎢會(huì)以液滴的形式飛濺出去。這些作用不但縮短器壁的壽命,而且產(chǎn)生的鎢雜質(zhì)(原子、塵埃顆粒、液滴)將嚴(yán)重威脅等離子體運(yùn)行。因此對(duì)瞬時(shí)熱流作用下鎢材料侵蝕的研究至關(guān)重要。本論文利用數(shù)值方法模擬研究了邊界局域模(edge localized modes, ELMs)熱負(fù)荷對(duì)鎢偏濾器靶板的侵蝕。在未來(lái)聚變裝置例如國(guó)際熱核聚變實(shí)驗(yàn)反應(yīng)堆(ITER)中強(qiáng)ELMs作用下,鎢偏濾器靶板會(huì)熔化,甚至汽化。在熔化和汽化過(guò)程中,相變界面隨時(shí)間不斷變化,熔化層流動(dòng)也會(huì)導(dǎo)致液體表面的位置隨時(shí)間不斷變化。本論文根據(jù)研究問題的需要,對(duì)所建立的模型采用有限體積法追蹤界面的位置,較理想地計(jì)算了鎢壁各個(gè)時(shí)刻的熔化速度及界面位置。論文在第一章緒論中介紹了工作的研究背景,第五章結(jié)論展望部分闡述了本工作的主要結(jié)論和創(chuàng)新點(diǎn),其它各章研究?jī)?nèi)容分別為：第二章,為研究ELMs熱流對(duì)鎢偏濾器靶板的侵蝕,把偏濾器靶板簡(jiǎn)化為一個(gè)平板,建立了一維對(duì)流傳熱模型,考慮了熔化、汽化和熱輻射效應(yīng)。利用自主開發(fā)的一維移動(dòng)邊界對(duì)流傳熱程序,首先,針對(duì)EAST,采用實(shí)驗(yàn)中取得的ELMs能流密度,計(jì)算了鎢靶板的表面溫度分布。計(jì)算結(jié)果顯示當(dāng)前第一類ELMs作用在鎢靶板上,在高約束模式運(yùn)行時(shí)間取32s情況下,靶板表面溫度從350 K增加到373 K,表明在當(dāng)前的參數(shù)范圍內(nèi),只要避免其它更嚴(yán)重的瞬時(shí)事件如破裂的發(fā)生,ELMs本身不會(huì)帶來(lái)嚴(yán)重的威脅。然后,將ELMs的能量增加到接近未來(lái)ITER的能量范圍,計(jì)算了ELMs能流密度峰值、間隙間能流密度、頻率取不同值時(shí),鎢偏濾器靶板的溫度分布和熔化厚度。研究發(fā)現(xiàn),對(duì)于不同幅值的ELMs,當(dāng)其它參數(shù)相同時(shí),偏濾器靶板的表面溫度(熔化厚度)與ELMs的能流密度峰值、間隙間能流密度以及頻率成正比。第三章,建立了二維對(duì)流傳熱的自洽模型,研究鎢偏濾器靶板瓦片在類似ITER運(yùn)行參數(shù)下ELMs過(guò)程中的侵蝕問題。針對(duì)不同形狀的鎢偏濾器靶板瓦片,考慮瓦片縫隙對(duì)等離子體行為、能流密度的影響,將利用二維邊緣等離子體動(dòng)力學(xué)程序(2d3v PIC-MCC)計(jì)算得到的不同形狀瓦片表面的能流密度分布,作為自主開發(fā)的二維移動(dòng)邊界對(duì)流傳熱程序的輸入?yún)?shù),計(jì)算了ITER中三種不同形狀瓦片(直角瓦片、斜邊瓦片、圓角瓦片)在發(fā)生ELMs過(guò)程中的侵蝕程度。研究結(jié)果表明,每一種形狀瓦片都是在面向入射等離子體一側(cè)棱邊區(qū)域的表面溫度最高,熔化厚度最大,汽化厚度最大：在相同的ELMs能流作用下,圓角瓦片的侵蝕程度最小,而且表面曲率半徑越大,侵蝕程度越小。論文中分析了這些結(jié)果的原因。第四章,通過(guò)將二維熱傳導(dǎo)方程與流體力學(xué)的納維斯托克斯方程結(jié)合,建立了二維流體動(dòng)力學(xué)模型,進(jìn)一步研究在ELMs強(qiáng)熱流轟擊下,鎢熔化層在表面張力、壓強(qiáng)梯度力、磁場(chǎng)力等力的作用下流動(dòng)時(shí),偏濾器靶板的侵蝕和形貌變化。在模型的數(shù)值求解中,采用交錯(cuò)網(wǎng)格的方法進(jìn)行離散,成功解決了液體表面追蹤的算法難點(diǎn),保證了鎢偏濾器靶板侵蝕程度計(jì)算的準(zhǔn)確性。研究結(jié)果表明,在ELMs過(guò)程中,鎢熔化層在表面張力、壓強(qiáng)梯度力、磁場(chǎng)力等各種力的作用下,會(huì)在熔化層兩側(cè)的邊緣區(qū)域形成類似山峰結(jié)構(gòu)的凸起,加重鎢偏濾器靶板的侵蝕程度。如果不考慮磁場(chǎng)力的作用,在空間分布為高斯形狀能流的作用下熔化層兩側(cè)的山峰結(jié)構(gòu)是對(duì)稱的。如果考慮磁場(chǎng)力的作用,熔化層在磁場(chǎng)力的作用下流動(dòng),使得熔化層兩側(cè)的山峰結(jié)構(gòu)不再對(duì)稱,迎著磁場(chǎng)力一側(cè)的山峰高度增加,另一側(cè)的山峰高度降低。
[Abstract]:The plasma wall material selection is a key problem in fusion devices Tokar Mark, but so far has not really suitable materials. Because the tungsten material has high thermal conductivity, high melting point, low sputtering rate, tritium retention rate is low, become the first choice of wall material. But the tungsten materials will melt in high instantaneous heat flux under the action of vaporization, melting layer in various forces (such as surface tension, magnetic force, pressure gradient force etc.) under the action of flow, resulting in material morphology, effect of heat transfer capacity, material induced arc, even liquid tungsten will in the form of droplets splashing out. These effects not only shorten the life of the wall, and the impurities (tungsten atoms, dust particles, droplets) will be a serious threat to operation. Therefore the plasma transient heat flow under the action of vital tungsten material erosion. Simulation Research on the boundary of local mode by using the numerical method in this thesis (edge localized modes, ELMs) erosion on the heat load of tungsten divertor plate. In future fusion devices such as the International Thermonuclear Experimental Reactor (ITER) in the presence of ELMs tungsten divertor plate will melt and vaporization. The melting and vaporization process, phase interface changing with time, melting the flow will lead to the liquid surface position changing with time. This paper according to the demand of research, tracking the interface of the model using the finite volume method, ideal to calculate the wall every time and melting rate of tungsten and the position of the interface. This paper introduces the research background of the work in the first chapter the fifth chapter is the conclusion and outlook section describes the main conclusions of this work and other innovations, the contents of every chapter are as follows: second chapter is study on erosion of ELMs heat flow on the tungsten divertor plate, the divertor plate is simplified A flat, one-dimensional convection heat transfer model, considering the melting, vaporization and heat radiation effect. By using the one-dimensional moving convective heat transfer process, the independent development of the first, for EAST, using ELMs to achieve current density, the surface temperature of tungsten target were calculated. The calculation results show that the current first class the role of ELMs in the tungsten target board in high confinement mode running time take the case of 32S, the surface temperature of the target plate was increased from 350 K to 373 K, that in the range of parameters in the current, as long as the other to avoid more severe transient events such as rupture, ELMs itself does not pose a serious threat. Then, ELMs the energy is increased to close to the energy range of ITER in the future, the ELMs peak value of flux density is calculated, the gap between the energy flow density, the frequency of different values, the temperature distribution and the thickness of the melting tungsten divertor plate. The study found that for different amplitude ELMs, when the other parameters are the same, the surface temperature of divertor target plate (melt thickness) and energy flow density of ELMs, the gap between energy density and frequency is proportional to. The third chapter, a self consistent model of two-dimensional convection heat transfer, erosion of tungsten divertor plate tiles in ELMs the process is similar to the ITER operating parameters in different shapes. The tungsten divertor plate tile, consider tile gap on the plasma behavior, influence of energy flow density, using two-dimensional edge plasma dynamics program (2d3v PIC-MCC) calculated by different shapes of tile surface the flux density distributions, as input parameters of two-dimensional moving boundary the convective heat transfer program developed by the three different shapes in the calculation of ITER tiles (rectangular tiles, tile edge, corner tiles) erosion degree in ELMs process. The results show that the shape of each tile The pieces are in the edge plasma facing incident side edge area, the maximum surface temperature, maximum thickness of maximum thickness of melting, vaporization: in the same ELMs can flow under the action of the minimum erosion degree of fillet tile, and the surface of the radius of curvature, the smaller the degree of erosion. The paper analyzes the reasons for these results. The fourth chapter and by combining the Navier Stokes equation of two-dimensional heat conduction equation and fluid mechanics, a two-dimensional hydrodynamic model, further research on ELMs heat under the bombardment of tungsten melting layer on the surface tension, pressure gradient force, magnetic field force under the action of force flow, erosion and topography in the divertor target plate in the numerical model, using the method of staggered grid discretization, successfully solved the difficulty of liquid surface tracking algorithm, which ensures the accurate degree of tungsten divertor plate erosion calculation. Results Show that in the ELMs process, the tungsten melting layer in the surface tension, pressure gradient force, magnetic force and other force, will form a convex structure similar peaks in the edge region of melted layer on both sides of the erosion severity of tungsten divertor target plate. If you do not consider the magnetic force for in space the distribution of Gauss shape can flow under the action of the melting layer on both sides of the mountain structure is symmetric. If considering the magnetic force, the melted layer under magnetic force flow, the structure of the melted layer on both sides of the mountain is no longer symmetric, side facing the magnetic force peak height is increased, the other side of the mountain height decreased.

【學(xué)位授予單位】：大連理工大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2016
【分類號(hào)】：O53;TL631.24

【相似文獻(xiàn)】

中國(guó)期刊全文數(shù)據(jù)庫(kù) 前10條

1 鄧柏權(quán),嚴(yán)建成,黃錦華;雙群模型在偏濾器材料濺射研究中的應(yīng)用[J];核聚變與等離子體物理;2002年03期

2 吳繼紅,張斧,許增裕,嚴(yán)建成;高熱負(fù)荷下偏濾器材料的溫度測(cè)量[J];核聚變與等離子體物理;2003年01期

3 洪文玉;嚴(yán)龍文;程均;趙開君;錢俊;王明旭;;HL-2A裝置送氣和加料的脫靶特性[J];核聚變與等離子體物理;2008年04期

4 鄧柏權(quán),謝中友,施漢文,H.H.Abou-gabal,G.A.Emmert;偏濾器靶板附近氦輸運(yùn)的蒙特卡洛模擬[J];核聚變與等離子體物理;1992年04期

5 周小兵,盛光昭,馮開明,朱毓坤,何開輝,龔學(xué)余;腐蝕和沉積對(duì)偏濾器靶板壽命的影響[J];核聚變與等離子體物理;1998年04期

6 鄧柏權(quán),嚴(yán)建成,黃錦華;自由表面液態(tài)鋰偏濾器靶板物理過(guò)程研究[J];核科學(xué)與工程;2000年04期

7 羅正明;等離子體平衡問題的兩種新的數(shù)值解法[J];核聚變與等離子體物理;1982年04期

8 潘宇東;用B2程序模擬HL-2A裝置附著狀態(tài)等離子體[J];核聚變與等離子體物理;2003年03期

9 歐靖;楊錦宏;;偏濾器運(yùn)行模式對(duì)托卡馬克邊緣區(qū)等離子體平行流的影響[J];物理學(xué)報(bào);2012年07期

10 朱毓坤,黃錦華,馮開明,鄧培智,李毅強(qiáng);聚變實(shí)驗(yàn)增殖堆FEB-E粒子抽除和抽氣系統(tǒng)(英文)[J];核聚變與等離子體物理;1999年04期

中國(guó)重要會(huì)議論文全文數(shù)據(jù)庫(kù) 前10條

1 歐靖;;數(shù)值分析東方超環(huán)的偏濾器運(yùn)行空間[A];第十五屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議會(huì)議摘要集[C];2011年

2 歐靖;;模擬分析在偏濾器脫靶狀態(tài)下的高馬赫數(shù)流[A];第十四屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議暨第五屆中國(guó)電推進(jìn)技術(shù)學(xué)術(shù)研討會(huì)會(huì)議摘要集[C];2009年

3 彭海鷗;陳銀華;時(shí)桂芬;鄭聚高;;Tokamak邊緣區(qū)磁場(chǎng)混沌結(jié)構(gòu)的研究[A];第十五屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議會(huì)議摘要集[C];2011年

4 桑超峰;孫繼忠;戴舒宇;徐倩;王德真;;托卡馬克偏濾器瓦片縫隙等離子體及燃料滯留模擬研究[A];第十六屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議暨第一屆全國(guó)等離子體醫(yī)學(xué)研討會(huì)會(huì)議摘要集[C];2013年

5 李敏;周宇璐;張寶玲;侯氫;;氦輻照鎢表面行為的分子動(dòng)力學(xué)模擬研究[A];第五屆反應(yīng)堆物理與核材料學(xué)術(shù)研討會(huì)、第二屆核能軟件自主化研討會(huì)會(huì)議摘要集[C];2011年

6 王德真;桑超峰;孫繼忠;;托卡馬克偏濾器等離子體波動(dòng)現(xiàn)象及其對(duì)壁材料侵蝕效應(yīng)[A];第十四屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議暨第五屆中國(guó)電推進(jìn)技術(shù)學(xué)術(shù)研討會(huì)會(huì)議摘要集[C];2009年

7 韓曉玉;余德良;段旭如;;HL-2A多道可見光譜系統(tǒng)測(cè)量的初步結(jié)果[A];第十三屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議論文集[C];2007年

8 李雪泓;李偉;盧杰;易萍;李棟;楊青巍;;HL-2A裝置偏濾器靶板熱沉積功率的測(cè)量[A];第十三屆全國(guó)等離子體科學(xué)技術(shù)會(huì)議論文集[C];2007年

9 趙偉;鐘光武;;用于ITER偏濾器等離子體診斷的朗繆爾探針熱分析[A];中國(guó)核科學(xué)技術(shù)進(jìn)展報(bào)告（第二卷）——中國(guó)核學(xué)會(huì)2011年學(xué)術(shù)年會(huì)論文集第7冊(cè)（核電子學(xué)與核探測(cè)技術(shù)分卷、脈沖功率技術(shù)及其應(yīng)用分卷、核聚變與等離子體物理分卷）[C];2011年

10 陳培銘;葉民友;彭學(xué)兵;卯鑫;宋云濤;萬(wàn)元熙;;CFETR偏濾器第一壁模塊熱工水力學(xué)初步分析[A];第八屆華東三省一市真空學(xué)術(shù)交流會(huì)論文集[C];2013年

中國(guó)博士學(xué)位論文全文數(shù)據(jù)庫(kù) 前5條

1 韓樂;偏濾器水冷鎢銅模塊傳熱與熱應(yīng)力問題研究[D];南京航空航天大學(xué);2015年

2 黃艷;邊界局域模熱負(fù)荷對(duì)鎢偏濾器靶板侵蝕的模擬研究[D];大連理工大學(xué);2016年

3 桑超峰;偏濾器等離子體鞘層及氫同位素在器壁材料中的滯留模擬研究[D];大連理工大學(xué);2012年

4 戴舒宇;偏濾器靶板的侵蝕及雜質(zhì)輸運(yùn)沉積的模擬研究[D];大連理工大學(xué);2014年

5 王東升;EAST偏濾器物理及雜質(zhì)注入行為研究[D];中國(guó)科學(xué)技術(shù)大學(xué);2012年

中國(guó)碩士學(xué)位論文全文數(shù)據(jù)庫(kù) 前7條

1 吳昊聲;基于SOLPS的CFETR雪花偏濾器的數(shù)值模擬[D];中國(guó)科學(xué)技術(shù)大學(xué);2015年

2 楊濤;托卡馬克中偏濾器材料的濺射行為研究[D];南華大學(xué);2015年

3 李永亮;偏濾器粗糙靶板對(duì)塵埃逃逸能力的影響[D];大連理工大學(xué);2012年

4 劉楠;水冷鎢銅偏濾器臨界熱流數(shù)值模擬方法研究[D];南京航空航天大學(xué);2013年

5 吳美薇;偏濾器靶板附近等離子體鞘層特性的模擬研究[D];大連理工大學(xué);2012年

6 陳培銘;CFETR偏濾器熱工水力設(shè)計(jì)與分析優(yōu)化[D];中國(guó)科學(xué)技術(shù)大學(xué);2015年

7 杜海龍;EAST高約束模式伴隨邊界局域模時(shí)能量在偏濾器靶板沉積模擬研究[D];大連理工大學(xué);2014年

，

本文編號(hào)：1496097