hcsb-mae示范堆包层子模块结构优化.docxVIP

hcsb-mae示范堆包层子模块结构优化

1中国建议的聚变示范堆

随着国际热核实验验证的进入，iter计划进入实施阶段，也被列入iter下一个目标聚变模型（梦想）的项目研究议程。欧洲、俄罗斯、日本和美国等早已开始进行聚变DEMO堆的概念设计研究,中国目前也正在积极推进具有先进托卡马克位形、稳态操作、先进偏滤器及超导线圈等特色的DEMO研究。研究内容主要包括氚增殖特性、排热和发电、集成和远程操作、安全、环境和燃料管理及经济性等。采用氦冷固态增殖包层(HCSB)的聚变示范堆是中国DEMO的技术路线之一,其主要核设计参数见文献。

在HCSB-DEMO聚变示范堆中,增殖包层位于中心等离子体与真空室壁之间,是非常重要的部件之一。其基本功能主要有两个:(1)通过聚变中子和氚增殖剂的核反应,在包层中增殖聚变反应需要的氚,由相应配套的提取系统进行氚提取;(2)沉积由中子和堆芯热辐射等方式进入包层的聚变反应热量,由冷却系统带出利用。这就需要对包层中的氚增殖剂、中子倍增剂及结构材料等进行合理的布置。前面曾针对HCSB包层进行过相关的研究,得到了初步的结果。为了在满足氚自持条件下使包层的温度分布更为合理,各材料都有足够的安全裕度,本文结合热工水力及中子学两方面对HCSB-DEMO示范堆的增殖包层进行初步的计算和优化分析。

2分析和计算流程

2.1共同冷却方案

由于面对等离子体的第一壁有较大的热负载,包层子模块的优化设计需要结合第一壁进行考虑。其冷却方式有两种:一种是采用第一壁、包层子模块单独冷却方式;另一种是采用第一壁、包层子模块共同冷却方式。本文首先对单独冷却情况进行分析,再根据优化得到的结果对共同冷却的情况进行讨论。

参考前面类似的设计,材料方面采用8MPa的氦气作为冷却剂,低活化铁素体钢(RAFM)作为结构材料,单体硅酸锂球床(38%的孔积率)作为氚增殖材料,双体铍球床(20%的孔积率)作为中子倍增材料。这些材料的热力学性能列于表1中。在结构方面,采用硅酸锂球床和铍球床沿径向间隔排列的方式,其间布置有冷却管板。

2.2热力场建模与网格划分

考虑到包层子模块结构和部件的安全性及反应堆的经济性,包层子模块的设计需要结合热工水力和中子学进行迭代优化。

本文中的热工水力和中子学分析分别采用FLUENT软件和MCNP软件。迭代计算模型如图1所示。图1a为中子学分析所用到的全堆模型示意图。由于包层子模块整体结构具有一定的上下周期性,热工水力分析是基于具有代表意义的一层冷却板来进行模拟。图1b为单独冷却情况下第一壁及包层子模块的模型示意图。鉴于包层子模块主要受面向等离子体部分第一壁的影响,虚线所示为简化后的分析模型。其各区径向分布与中子学分析模型采用相同的尺寸参数。

迭代计算的主要流程如下:根据模型的初步设计,利用MCNP软件进行中子学分析优化,使其得到尽可能高的中子学性能参数,例如合理的TBR;由中子学分析得到的平均功率密度分布如图3所示,利用FLUENT软件进行热工水力分析,得到各区的温度分布;结合中子学和热工水力的计算结果,在保证满足氚自持要求的情况下,对模型结构布局进行调整,重新进行分析计算,经过多次反复迭代优化,确保低活化铁素体钢、硅酸锂球床及铍球床的最高温度均低于许用温度限值,同时确保包层子模块的冷却剂出口温度达到500℃以上,冷却剂压降也要符合工程上的要求,即低于工作压力的4.9%。

热工水力分析模型的网格划分是基于ICEMCFD软件进行的。其网格划分方式主要有两种:四面体网格划分和六面体网格划分。与高纵横比的六面体单元相比,四面体单元不允许存在大的纵横比,否则会影响单元的偏斜,可能妨碍计算的精确与收敛,而降低其纵横比又会大大增加网格单元数量,计算经济性较差,尤其是对于本文中采用的细长管道;此外,对于长管道的流动,依靠六面体网格单元还可以尽可能地降低数值耗散。因此本文中的网格划分均采用六面体单元,图2为计算模型整体的六面体网格及其质量示意图。

由于反应堆为D型环结构,外侧第一壁的中子壁负载比内侧第一壁中子壁负载要大,致使外侧包层的功率密度比内侧包层要高,冷却相对比较困难,是反应堆内的设计难点之一。本文是基于HCSB-DEMO外侧包层的尺寸参数和功率密度分布进行热工水力和中子学的迭代优化。

3不同增殖区的优化结果分析

包层子模块中平均功率密度径向分布示意如图3所示。可以看到,热沉积主要集中在硅酸锂区,且随着距第一壁越来越远,硅酸锂区和铍区的平均功率密度均呈明显的下降趋势。这是由于随着中子在材料中的慢化吸收,各区的中子通量逐渐减少,引起反应率降低所致。鉴于后侧各区的功率密度较低,设计的模型尺寸中硅酸锂区和铍区的厚度可适当随第一壁距离而增加。由表1可见,硅酸锂区的导热系数非常低,不宜过厚,否则很容易导致材料最高温度超过限值,因此除最后