1. 引言
中国聚变工程实验堆CFETR (China Fusion Engineering Test Reactor)是我国正在进行工程设计的托卡马克装置,任务是开发设计出先进的稳态超导托卡马克,为未来建成可运行的聚变发电反应堆奠定坚实的科学基础 [1]。其第一阶段计划的科学目标是达到200 MW的聚变功率,在未来的第二阶段其目标聚变功率更是达到1 GW,并进一步开展DEMO的验证工作。为了实现这一聚变领域的目标,中国CFETR聚变装置汲取包括ITER在内的磁约束核聚变物理及工程技术经验,努力解决一些关键性的难题。例如,氚自持难题(TBR > 1.0)、包层以及偏滤器等一些未解决的问题 [2]。
偏滤器是CFETR装置的关键部件,部件的主要功能是排出热量与排灰 [3]。偏滤器排出热量来自最外封闭磁面以外开放磁面驱动的等离子体,其携带高额能量打到偏滤器上。能量由等离子体转移到偏滤器上,等离子体中性化,由偏滤器排出。根据CFETR装置的聚变反应功率目标,偏滤器靶板稳态热流高达数10 MW/m2,这要求偏滤器具有极强的排除热量能力 [4]。因此在CFETR物理和工程设计方案中,偏滤器的设计是一项重要研究课题 [5] [6] [7]。
关于CFETR偏滤器相关设计方面,目前已有诸多进展。卯鑫等人已经开展偏滤器结构设计,包括类ITER偏滤器、雪花偏滤器和类ITER雪花偏滤器三种结构,并验证合理性 [8];刘秀等人基于以前研究,提出CFETR偏滤器靶板的概念设计,进行热力水工、力学分析,验证了模型的可靠性 [9];黄文玉等人基于中国聚变工程试验堆(CFETR)偏滤器位形,初步设计了指状、平板和T形三种氦冷偏滤器模型,进一步分析之后提出了混合型CFETR氦冷偏滤器初步结构 [10];周一夫等人进行了在CFETR不同位置充气下的辐射偏滤器模拟研究,得到了相对更好的杂质屏蔽效果 [11];张传家等人使用SOLPS对CFETR长腿偏滤器的峰值热流进行了模拟研究,内、外靶板热流峰值分别为15.5 MW/m2、12.5 MW/m2 [12]。
本文使用TSC程序 [13] 对等离子体放电的爬升、平顶、破裂阶段进行了数值模拟。结合CFETR物理参数,在TSC程序中添加偏滤器靶板的几何结构,主要包括内靶板、外靶板、Dome板。同时模拟准雪花偏滤器位形、下单零偏滤器位形放电,使最外封闭磁面以外的等离子体打击到偏滤器靶板上。模拟得到了两种偏滤器位形放电下的内、外靶板及Dome板热流,重点分析了平顶阶段偏滤器靶板热流分布,对下单零偏滤器位形和准雪花偏滤器位形的热流分布进行了比较。TSC程序处理粒子与能量输运模型是一维横越磁面输运,并用粘滞系数替代对流微商项,输运模型较为简单,因此,热流模拟结果可以起到一定参考作用。
2. TSC数值模型
TSC程序具有非刚性、可变形等离子体模型的特点,因此采用TSC程序对CFETR装置进行放电模拟。1986年初版TSC在美国普林斯顿大学等离子体物理实验室(PPPL)编写开发完成,随着面临的问题不断变化,TSC不断发展。时至今日,TSC程序已经用于世界各地多个大型托卡马克装置的放电模拟。TSC是自由边界轴对称的托卡马克数值模拟程序,模拟自由边界等离子体随时间的演化,求解压力和温度的面平均输运方程,同时程序包含了辅助加热、电流驱动等物理模型,可用于控制系统、外部加热等的模拟 [14] [15]。
一个轴对称的环形几何系统中,极向通量和环向场函数可以用来表示磁场,表达式如下:
(1)
其中
为环向角,
为角向磁通,g为环向场函数。
为了便于对流场的可压缩部分和不可压缩部分分别处理,因此把等离子体的动量密度
表示为:
(2)
其中A为流函数,
为环向分量,
为一个势函数。
TSC程序计算区域分为三块:等离子体区、真空区、导体区,且是自由边界程序,等离子体区与真空区分界面随时间变化,穿越边界时方程形式不变且解连续。在MHD-Maxwell方程组中引入质量放大因子及粘滞放大因子等参数。
等离子体静态的平衡条件是:
(3)
TSC修正惯性项的方法是通过人为的增加等离子体质量,略去对流微商项,取粘度力为:
(4)
其中
、
为可压缩粘滞系数和不可压缩体积粘滞系数。修正后的等离子体力平衡方程如下所示:
(5)
TSC输运模型采用磁面平均输运方程组,由等离子体微分粒子密度
、等离子体熵密度
和电子熵密度
的一维演化方程组成,输运计算区域包括等离子体和刮削层区域。
(6)
(7)
(8)
式中
为环向磁通量,
、
、
分别是粒子、电子和离子的外源能量,
、
是离子和电子热流,
、K是环向电压和环电流,
是辐射引起的能量损失。
TSC程序中,完全电离原子的轫致辐射是最主要的辐射过程,面平均辐射功率密度可由下式计算:
(9)
式中
为有效电荷数,可以根据不同的放电情况调节参数。
3. CFETR下单零和准雪花偏滤器位形放电模拟
3.1. CFETR偏滤器结构
CFETR II期设计目标是稳态燃烧等离子体能产生1 GW聚变功率,如此高的聚变功率对第一壁和偏滤器提出了极高的要求,而偏滤器又承担着排热与排灰的功能。因此必须要探索研发具有更好性能的材料,同时开发更好的偏滤器位形,降低偏滤器受热峰值。偏滤器结构图如图1所示,图中内外靶板的板面与最外封闭磁面形成较小的夹角,使得靶板结构呈弧形。相较于直板式偏滤器靶板其具有更大的热沉面积,沉积的峰值热流相对减小,在此基础上进行偏滤器位形设计会有更好的效果。
Figure 1. Diagram of the practical teaching system of automation major
图1. CFETR偏滤器三维结构图
偏滤器在CFETR装置中位置如图2所示,在装置中偏滤器由特制的盒体支撑。图2标注了8个中心
Figure 2. Three-dimensional drawing of CFETR device
图2. CFETR装置三维图
螺线管线圈(CS线圈)、6个极向场线圈(PF线圈)以及偏滤器线圈(DC1线圈),这些线圈共同作用约束等离子体,同时可控制形成更好的偏滤器位形,线圈的参数如表1所示。图2所示的CFETTR装置拉长比为2,大半径7.2米,小半径2.2米,等离子体电流14 MA,CFETR的主要参数列于表2。
Table 1. Geometric parameters of CFETR coil
表1. CFETR线圈几何参数
Table 2. Physical parameters of CFETR device
表2. CFETR装置物理参数
3.2. 下单零偏滤器位形放电模拟
通过对极向场线圈电流的不断修改得到了下单零偏滤器位形,模拟放电总过程共432秒,等离子体放电主要参数如图3所示。从图3可以看出,等离子体电流在70秒左右达到14 MA,维持放电到400秒等离子体电流开始下降,在432秒等离子体停止放电,模拟结束。
Figure 3. Simulation parameters of lower single null divertor configuration discharge; (a) Plasma current; (b) Large radius; (c) Small radius; (d) Plasma density
图3. 下单零位形放电模拟参数;(a) 等离子体电流;(b) 大半径;(c) 小半径;(d) 等离子体密度
模拟整体过程可分为三个阶段:第一阶段,等离子体电流从0.5 MA开始呈线性增加,一直增加到14 MA并维持在14 MA左右。同时等离子体的大半径和小半径也在这一阶段达到稳定,在此期间大、小半径有多次波动,且大、小半径波动趋势相反。大半径增加则小半径减小,小半径增大则大半径减小,由图3可见大小半径在第一阶段后期基本稳定在7.2米、2.2米。图3(d)是等离子体密度,在这一阶段等离子体密度也大致呈线性增加由0.05 × 1020 m−3逐步增加到0.59 × 1020 m−3。第二阶段,等离子体电流变化不大,最大可达到14.62 MA,保持比较稳定的等离子体放电状态。在核聚变过程当中这一阶段是最重要的阶段,等离子体大部分时间都运行在平顶阶段。区别于上一阶段的多次波动,建立平衡后的大半径和小半径在这一阶段保持在7.2米、2.2米,只有很缓慢、微小的变化。等离子体密度和之前的增长趋势相比减缓许多,并且等离子体密度在200秒增长到1.35 × 1020 m−3,之后一直维持不变。第三阶段,等离子体垂直位移导致破裂,从400秒开始等离子体电流有较大幅度的减小,大半径和小半径也有相应的减小。在432秒等离子体电流降到12.11 MA、大半径减小到6.91米、小半径减小到1.715米,等离子体停止放电,模拟结束。
调节线圈电流是模拟等离子体放电以及模拟偏滤器位形的主要手段,图4给出了各个极向场线圈电流的变化图。图4(a)为中心螺线管线圈(CS Coils)电流图,图4(b)为极向场线圈(PF Coils)和线圈DC1电流图。在第一阶段CS线圈主要作用是提供欧姆加热场加热等离子体,PF线圈也提供一定的加热作用,但PF线圈更重要的作用是控制等离子体位形,DC1是偏滤器线圈起到调节等离子体偏滤器位形的作用。图5给出100秒时CS、PF线圈共同作用下的下单零磁面位形,此时等离子体电流为14.62 MA,大半径和小半径为7.23 m、2.286 m,中性束注入40 MW功率进行辅助加热。
3.3. 准雪花偏滤器位形放电模拟
准雪花偏滤器位形模拟共持续365秒,其放电主要参数如图6所示。从图6可以看出,准雪花偏滤器位形的等离子体电流同样在70秒左右达到14 MA,维持放电到340秒等离子体电流开始下降,在365秒等离子体停止放电,模拟结束。
Figure 4. Lower Single Null coil current evolution; (a) Cs coil current; (b) PF and DC1 coil current
图4. 下单零位形线圈电流演化;(a) CS线圈电流;(b) PF和DC1线圈电流
Figure 5. Lower single null divertor configuration
图5. 下单零偏滤器位形
Figure 6. Simulation parameters of Quasi-Snowflake divertor configuration discharge; (a) Plasma current; (b) Large radius; (c) Small radius; (d) Plasma density
图6. 准雪花位形放电模拟参数;(a) 等离子体电流;(b) 大半径;(c) 小半径;(d) 等离子体密度
图7给出了准雪花偏滤器位形的CS、PF、DC1线圈电流随时间变化图,图7(a)可以看出CS线圈在70秒之前呈线性变化,之后基本不变。CS线圈通过电流变化为等离子体提供欧姆加热场,CS线圈从0秒时50 kA/匝开始变化,在70 s维持不变。CS1U和CS1L从正向50kA/匝变为负向24.7 kA/匝、34.4 kA/匝,CS1U、CS1L位于装置中间位置还起到控制等离子体横向位形的作用。偏滤器线圈电流最大为63.5 kA,超出其电流限值15%。
图8为100秒时的准雪花偏滤器位形,此时等离子体电流为14.7 MA,大半径和小半径为7.28 m、2.292m,中性束注入40MW功率进行辅助加热。图中下半部分有两个零点,两个零点之间的距离为D,等离子体小半径为a,衡量准雪花与精确雪花接近度的一个简单方法是通过两者的比值,如下式:
(11)
本次模拟D约为2米,a为2.2米,可以通过简单的计算来判断接近程度。
Figure 7. Quasi-Snowflake coil current evolution; (a) Cs coil current; (b) PF and DC1 coil current
图7. 准雪花位形线圈电流演化;(a) CS线圈电流;(b) PF和DC1线圈电流
Figure 8. Quasi-Snowflake divertor configuration
图8. 准雪花偏滤器位形
精确的雪花偏滤器位形在拓扑结构上是不稳定的,我们可以把两个零点距离在一定范围内的偏滤器位形按照精确雪花进行判断,判断标准如下:
(12)
式中
为热流宽度,是热流研究的重要参数,还无法可靠地计算出,因此使用经验估算出CFETR的热流宽度约为1毫米 [16]。式12计算结果约为0.17米,当D小于0.17米时可视为精确雪花位形,图8所示的准雪花偏滤器位形为雪花加结构(Snow Flake-Plus, SF+) [17]。
CS、PF线圈所提供的伏秒数是一定的,在CS、PF线圈减少甚至不提供伏秒数时需要辅助加热、电流驱动维持等离子体持续放电。在模拟准雪花偏滤器位形时加入如中性束和离子回旋波进行辅助加热,CS线圈在基本70秒结束变化,PF线圈为了控制形成偏滤器位形在70秒后会有电流变化,集中在70秒到100秒附近。综合考虑在85秒时中性束开始注入20 MW功率,100秒时功率加大到40 WM,200秒时加到57MW维持不变。在170秒加入离子回旋波,功率为10 MW,200秒时上升为20 MW维持不变。
4. 偏滤器靶板热流模拟
4.1. 偏滤器靶板模拟
在CFETR偏滤器的设计当中,偏滤器靶板上的热流峰值以及热流分布是比较重要的参考因素。上一节数值模拟了两种偏滤器位形下的等离子体放电全过程,等离子体电流、大半径以及小半径等主要物理参数基本符合CFETR设计要求。此时在TSC程序中添加对应的偏滤器靶板,如图1所示是CFETR偏滤器第一壁的三维模型,包括内、外靶板,内、外反流板及Dome板,根据偏滤器的实际设计参数并结合TSC程序加入了偏滤器靶板。
调节CS、PF线圈电流获得准雪花偏滤器、下单零偏滤器位形,使最外封闭磁面以外的粒子打击到偏滤器靶板上,模拟计算可以获得对应的靶板热流。图5是下单零偏滤器位形,可以看到X点在Dome板上方。下单零偏滤器位形相较于限制器位形的优点:1) 真空室与等离子体主区域距离较远,在磁约束位形下杂质通过刮削层(SOL scrape-off layer)远离等离子体主区域,减轻等离子体污染;2) CFETR运行时产生的杂质进入刮削层沿着磁力线达到偏滤器靶板,下单零偏滤器位形下磁力线较长有利于杂质冷却,降低偏滤器靶板热流峰值。
随着核聚变功率、等离子体电流不断增大,等离子体区域向外辐射的热量不断加大,偏滤器靶板上的热流也在不断加大,下单零偏滤器位形已经无法满足更大聚变功率聚变堆的设计目标。此时提出了一种新的偏滤器位形即雪花偏滤器位形,在精确雪花偏滤器位形中X点附近分为六个区域,形状似雪花,顾称之为雪花偏滤器位形。图8是准雪花偏滤器位形,等离子体位置变化打在偏滤器靶板上的热流也会随之变化。
4.2. 偏滤器靶板热流结果与分析
对靶板上的热流结果进行分析,比较下单零与准雪花偏滤器位形下偏滤器靶板热流分布与峰值大小。图9和图10是准雪花、下单零偏滤器位形下的靶板热流,图9(a)、图10(a)是内靶板热流,图9(b)、图10(b)是Dome板热流,图9(c)、图10(c)是外靶板热流。0~70秒是爬升阶段,此时等离子体位形建立,未形成偏滤器位形,热流来自线圈电流提供的加热场。据图9、图10可知无论是准雪花偏滤器位形还是下单零偏滤器位形在这一阶段的热流值都比较小,且两种偏滤器位形三块靶板上的热流分布比较相似。爬升期下单零内靶板的热流峰值为1.372 MW/m2,而准雪花偏滤器位形内靶板热流峰值为0.692 MW/m2,另外两个靶板热流值要小于内靶板。
Figure 9. Heat flux of Quasi-Snowflake divertor configuration; (a) Inner target plate; (b) Dome plate; (c) Outer target plate
图9. 准雪花偏滤器位形热流分布;(a) 内靶板;(b) Dome板;(c) 外靶板
Figure 10. Heat flux of Lower Single Null divertor configuration; (a) Inner target plate; (b) Dome plate; (c) Outer target plate
图10. 下单零偏滤器位形热流分布;(a) 内靶板;(b) Dome板;(c) 外靶板
准雪花偏滤器位形放电平顶阶段为70~290秒,下单零偏滤器位形放电平顶阶段为70~400秒,由图9、图10对应的时间段可以观察内、外靶板及Dome板的热流。准雪花偏滤器位形下热流主要打击在内靶板和外靶板上,图9(b)等离子体稳态放电阶段Dome靶板上几乎没有热流,只是在93秒到100秒之间出现了一个峰值为24 MW/m2的一个波动热流,与准雪花偏滤器位形比较吻合。图9(a)是内靶板热流图,热流在靶板上并不是一直打击在固定的位置上,打击位置随着时间而变化。观察打击位置变化,可发现变化规律与等离子体大半径、小半径变化有关,随着等离子体大、小半径变化开放磁面驱动的热流也扫在靶板不同位置上。内靶板稳态热流峰值为23.34 MW/m2,其它位置处基本在20 MW/m2位左右,大部分时间热流值小于20 MW/m2。图9(c)准雪花偏滤器位形外靶板热流打击位置变化与内靶板基本相同,也是跟随着等离子体大、小半径变化开放磁面驱动的热流扫在靶板不同位置上。外靶板热流峰值为23.1 MW/m2,其220秒之前靶板热流小于20 MW/m2,220秒之后部分位置热流值要稍大一些。偏滤器靶板传热性能基本能够满足热流为15 MW/m2,稳态热流为20 MW/m2是目前CFETR偏滤器靶板的设计目标。
图10对应时间段里可以看到下单零偏滤器位形平顶阶段热流分布图,与准雪花偏滤器位形不同的是下单零位形热流主要打击靶板是内靶板与Dome靶板。图10(c)是外靶板热流分布图,可以看出下单零偏滤器位形下外靶板基本没有热流打击,仅仅在外靶板最外侧有小于1 MW/m2的热流扫了十几秒。图10(a)是内靶板热流分布图,内靶板热流峰值达到39 MW/m2,部分位置热流从130秒就已经超过20 MW/m2,174秒时热流值超过30 MW/m2,之后热流一直维持在大于30 W/m2。两种偏滤器位形在等离子体放电14 MA的情况下,准雪花偏滤器内靶板热流峰值比下单零偏滤器内靶板热流峰值小15.6 MW/m2,且大部分时间准雪花偏滤器内靶板热流在20 MW/m2左右比小单零小了十几MW/m2。可见在平顶阶段,准雪花偏滤器位形下的内靶板热流比下单零偏滤器内靶板热流小的多。图10(b)是Dome靶板热流变化图,其峰值为21.7 MW/m2,大部分热流维持在10 MW/m2左右。
5. 结论
本文采用TSC程序模拟了CFETR装置的两种偏滤器位形放电,即下单零偏滤器位形和准雪花偏滤器位形,模拟得到等离子体电流、大半径和小半径等主要物理参数很好的符合设计要求。在TSC程序中添加靶板数据,进行等离子体放电模拟,得到两种偏滤器位形下热流在靶板上的分布以及热流峰值。通过对比两种偏滤器位形下靶板热流分布和峰值,可知热流分布与等离子体大小半径变化有关,等离子体放电平顶阶段时准雪花偏滤器位形靶板热流值为20 MW/m2左右,准雪花偏滤器位形下内靶板热流比下单零内靶板热流值小十几MW/m2,优于下单零偏滤器位形。计算结果为下一步热结构多物理场耦合分析提供了热源分布,虽然TSC的输运模型较为简单,但本文的模拟对CFETR的工程设计具有一定的参考价值。
基金项目
国家自然科学基金(11975022);国家重点研发计划(2018YFE0309101, 2017YFE0300603, 2018YFE0310400);安徽省自然科学基金(1908085MA26, 2008085QA37)。