CARR全堆芯扩散计算研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 CARR全堆芯扩散计算 吕征孙志勇 沈峰 肖诗刚 中国原子能科学研究院，北京275信箱33分箱，102413 摘要t本文采用少群扩散方法对中国先进研究堆(cA非)进行计算。计算采用基于物理启动、同位素生产 以及空间燃耗不同考虑的全堆芯模型．通过对选定的换料方案进行堆芯第1～3炉的临界计算，得到了各炉 初末期的反应性、中子注量率以及功率分布．并且采用隐式燃耗方法模拟了cA职堆芯燃耗过程，得到了堆 芯燃耗分布，分析了堆运行过程中产生的裂变产物份额并给出了中毒反应性效应。确定了CARR的基本物理 特性参数，完成了CARR物理详细设计工作为CARR最终安全分析报告的编写以及将来的运行提供必要的数据． 关键词：中国先进研究堆(CARR)，全堆芯，少群扩散，反应性，中子注量率 1引言 中国先进研究堆(CARR)将在2l世纪初建成，作为我国新一代研究堆的代表，为满足 21世纪我国国民经济建设及核科学技术发展的需要，它应该是一座多用途、高性能的反应堆。 2001年12月，CARR初步设计和初步安全分析完成后，随着工程的进展，在堆本体结构和 堆内材料成分上有了一些变化，引起了堆芯物理特性参数相应的变化，因此必须对CARR堆 芯物理特性参数重新进行计算。与初步设计相比，本计算考虑了如下因素： CARR堆在物理调试试验和实际运行中调节棒和补偿棒的高度不完全一致，应该采 用新的计算模型对调节棒和补偿棒加以区分； 作为CARR重要用途的同位素生产孔道布置不对称，为了满足生产应用的需要， 新的计算模型应该能对不同位置的垂直孔道加以描述： 由于堆芯布置的不对称以及控制棒的不同高度，导致中子注量率空间分布的不对 称，相应的燃耗空间分布也是不对称的，而且在物理启动中逐渐增加燃料元件时堆 芯也是不对称的，因此必须对堆芯每盒组件加以区分。 而这些问题都是初步设计采用的1／4对称堆芯模型无法解决的，因此有必要对CARR进 行全堆芯计算。由于几何描述功能强大的MCNP程序目前基本上用于校核计算，很少作为设 计工具，而且MCNP本身无法考虑燃耗，所以本文使用扩散程序在新的计算条件下采用全新 的计算模型对CARR堆芯进行系统的物理分析，给出一整套稳态物理特性参数，为最终安全 分析报告提供参数，在此基础之上，探索研究适合物理启动和同位素生产应用计算的模型， 为CARR工程的下一步工作做好准备。 2计算条件和计算模型 2-，f1计算条件 CARR是一座稍加压的轻水冷却、重水反射的反中子阱型紧凑结构池式反应堆。图l为 481 第十一届反成堆数值*算#植f输t学术食没Ⅱ2006年反R堆蜘目会议 CAP．R堆芯布置圈。CARP,采用2”U的富集度为19 75wPA，含铀密度为4．39／era’的u，s12．AI 弥敢体作燃料芯体、铝台金作包壳的平板型燃料组件． 燃料活性区高度gSem。 本计算沿用初步设计采用的少群扩散程序系统 WIMSD4+cITATION进行．其中W1MSD-4用来作橱元 计算以获得少群参数．CITATION则利用W1MSD-4计 算所得的少群参数作整堆扩散计算，该计算系统的可靠 性在初步设计中通过对参考堆的计算和萤特卡罗程序 MCNP的枝核计算得到了很好的验证。 相对于初步设计．本计算的计算条件有以下改变： 图1 CARR堆芯布置圉 (I)初步设计中所有堆内材料构成分取理论值．本 计算采用样品实铡值．其中燃料组件包壳材料6061A1和燃料芯体材料OjSirAI杂质含量的 变化比较大．具体鼓值见表1。 (2)堆本体一些结构发生了变化，例如堆芯容器壁厚由12mm变为10ram，重永箱外壁 厚由15ram变为24ram，本计算所有结构尺寸全部采用由本体提供的堆新数据。 袭1堆内材料成分变化对照表(表中数值为核素重量百分比，单位为％) 22计算模型 为了更加真实准确的描述CARP．堆甚结构，使计算结果更加真实可信．在进行全堆芯计 算过