冷冻靶的环境热辐射屏蔽技术研究.docxVIP

? ? 冷冻靶的环境热辐射屏蔽技术研究 ? ? 陶朝友，刘喜川，林 伟，杨 洪，代 飞，王 凯 (中国工程物理研究院 激光聚变研究中心，四川 绵阳 621900) 美国NIF设计的间接驱动冷冻靶在黑腔中心悬挂有1个靶丸，靶丸内含有低温冷冻的燃料[1-3]。高强激光束穿过激光入射孔(LEH)打到黑腔内壁，将激光转换为X射线来辐照靶丸，靶丸内爆并压缩、加热氘氘(DD)/氘氚(DT)燃料冰层到一定条件以实现热核聚变[4-6]。要获得足够多的聚变能量，需要初始燃料冰层足够均匀、光滑，以尽量缓解扰乱冰层压缩的瑞利-泰勒不稳定性[7-10]。根据历史研究经验，燃料冰层对热辐射环境极为敏感[11]。因此，在极低温下获得内表面光滑、厚度均匀的燃料冰层充满巨大挑战。大部分环境热辐射处于红外波段，且黑体辐射峰值在10 μm左右[11]。减弱外界环境辐射对冷冻靶靶丸温度均匀性的影响是提升燃料冰层品质的一种重要方法。美国NIF点火靶在LEH处采用封口膜+暴风窗双层膜设计方案，其中LEH膜镀铝层厚度约25 nm，主要用以减弱注入到黑腔内部的环境热辐射；暴风窗采用镀碳薄膜，主要用于缓解LEH窗口的结霜问题[12]。物理实验、冰层表征、DT冰层制备工艺特点等要求限制了其封口膜铝层厚度的进一步增加。然而，对于我国以DD为燃料的间接驱动柱腔冷冻靶，LEH结霜问题已通过其他方式解决，LEH膜环境辐射屏蔽性能成为冷冻靶的关键设计参数之一。 本文主要分析影响DD冰层的因素，并通过数值模拟方法研究具有不同透光率的镀铝聚酰亚胺(PI)膜时，靶丸外表面和冰层内表面的温度场分布。同时，研究在LEH膜上镀不同厚度铝层时，实验过程中靶丸内冰层的时空演变规律，以期为优化冷冻靶的外界环境辐射屏蔽，从而制备出高质量的冰层奠定基础。 1 实验 1.1 实验过程 本文所用冷冻靶和屏蔽罩的结构如图1所示，冷冻靶被低温屏和常温屏所包裹。黑腔由连接在导冷杆上的导冷硅臂冷却。外部腔体为铝套筒，铝套筒包裹黑腔，黑腔内充有1 kPa左右的氦气。腔的两端为LEH，其上有厚度约为0.5 μm的PI薄膜，膜上镀不同厚度的铝层。靶丸为碳氢聚合物，内含DD燃料冰层。靶丸通过石英充气管支撑在黑腔中心位置。实验过程中，由Gifford-McMahon(G-M)制冷机提供的低温源将靶丸中的氘气冷冻成液态或固态DD，同时采用水平方向的X射线相衬成像系统表征装置记录靶球内DD的变化过程。当黑腔的温度降至20 K时，靶丸内的氘气(26.5 kPa)逐渐开始液化，18.7 K左右完全液化，靶球内液面不再发生变化，由于重力及液态DD表面张力的作用，靶球中液态DD的分布呈底部多、顶部少的状态，形成弯月形亮线[13]。通过控制电阻加热器和温度传感器来控制传热腔上下的温度梯度。本文计算模型中，上述所有靶零件材料的物性参数(表1)均在18 K附近取值。 图1 冷冻靶及其防辐射屏蔽罩结构示意图Fig.1 Schematic diagram of cryogenic target and its radiation shielding cover 表1 主要材料的热物特性(18 K)[14]Table 1 Thermophysical property of main material (18 K)[14] 1.2 计算模型和数值方法 根据柱腔冷冻靶结构和运行工况，基于COMSOL软件建立了用于求解靶丸温度场的三维数值模型。在网格划分方面，考虑到靶丸外表面为热梯度变化的主要区域，靶丸与填充气体的交界区域采用多层结构化网格过渡，以建立流体边界层。靶丸温度场是此模型的关键计算区域。在参数化计算前，首选采用与冷冻靶常规实验相符的基准计算工况进行数值模型验证：外界环境辐射温度为120 K，LEH膜透射率为0.4，上下硅臂末端温度均为18.5 K，零件接触界面两侧温度自动耦合，不考虑接触热阻。通过对比不同网格密度的基准数值模型计算结果，获得满足网格无关性的模型离散设置参数：靶丸壳层及冰层网格尺寸约为10 μm，氦气中心部分网格尺寸约为50 μm，黑腔及套筒区域网格尺寸约为80 μm，氦气边界层网格由靶丸、黑腔部分按1.2倍尺寸增长比自适应过渡。模拟中连续方程、动量方程和能量方程[15]为： (1) (2) (3) (4) 式中：ρ为密度，kg/m3；t为时间，s；u、v为速度矢量u在x、y方向的分量；cp为比热容，J/(kg·K)；μ为黏度，Pa·s；T为参考温度，K；k为传热系数，W/(m2·K)；St为流体的内热源及由于黏性作用流体机械能转换为热能的部分，W。 由于黑腔内的氦气和靶丸中心的氘气温度变化较小，密度变化也小，在动量方程中使用了Boussinesq假设。对于理想气体，热膨胀系数β=-1/T，式(1)～(3)适用于气体区域，式