高温气冷堆技术及堆芯中子动力学仿真计算解决方案.ppt

* * * * * * * 在发生停堆事故后，反应堆仍有大量的衰变热产生，能否将反应堆的衰变热载走是保证核电站的关键，模块式高温气冷堆采用了热辐射、自然对流的自然机制将余热排出，能可靠的防止发生堆芯熔化，保证了核电站的安全性。 * * * * * 这张图表示了HTR-10的主要部分，左面是核反应堆部分，其中心是由27000个燃料元件组成的球床堆芯，裂变产生的核能由氦冷却剂载出到蒸汽发生器，加热蒸气发生器中的水，产生的蒸汽推动透平发电。 本程序对于这个动态过程的仿真计算结果是收敛的，内迭代收敛精度与外迭代次数对动态仿真结果的综合影响在1%以内。 也就是说，用本程序求解形状函数时，做1次外迭代即可符合模拟机动态参量测试的精度要求。 计算时间比较 实时仿真算例 控制棒误提升 单根控制棒弹出事故 控制棒误提升 假设初始时反应堆处于正常运行状态，即反应性ρ = 0。这时发生了控制棒整体误提升的情况，假定这一提升过程动作时间为10s。 所谓过功率保护是指，假设堆功率上升到120%~130%时，触发过功率保护信号，使控制棒下插，引入负反应性，假定动作时间为30s。同时，氦气风机立刻关闭，假定关闭时间为50s。反应堆在一段时间内功率将为零，安全停闭。 仿真程序与集总参数法的热工反馈模型耦合，对这一动态过程进行计算，相对功率计算时考虑了裂变产物的剩余发热。 考虑过功率保护 提棒过程模拟： 插棒过程模拟： 仿真结果曲线 ATWS事故模拟 考虑上述的控制棒误提升的ATWS事故情况，假设触发过功率保护信号之后，由于事故原因，控制棒没有下插动作，仅靠氦气风机关闭使堆内温度升高，引入负反应性温度效应，使得反应堆自动停闭。 使用耦合了热工反馈模型的模块式高温气冷堆中子时空动力学仿真程序对这一过程进行模拟计算。 ATWS仿真结果 单根控制棒弹出事故 “弹棒”事故是压水堆安全分析中要考虑的严重事故，对于模块式高温气冷堆来说，这一事故情况也是存在的。假设在0.5s内，某根控制棒由于压力作用或其它原因从侧反射层弹出，通过改变θ方向1个网格( 0 ≤θ π/12 )控制棒区的第4群宏观吸收截面来模拟这一过程。 考虑过功率保护，控制棒全部下插引入负反应性，同时氦气风机关闭。通过改变整个控制棒区第4群宏观吸收截面模拟控制棒下插过程，假定其变化规律如下： 考虑过功率保护仿真结果 ATWS仿真结果 创新意义 针对模块式高温气冷堆模型，研究开发了三维圆柱几何多群中子时空动力学实时计算程序系统，为HTR-PM高温堆示范电站的实时模拟机的自主开发提供了核心模拟计算工具。 建立了基于集总参数法的热工反馈模型，初步实现了模块式高温气冷堆三维中子时空动力学的物理热工耦合模拟。该模型的特点在于简洁高效，可用来检验物理计算模型的合理性。 将求解大型稀疏线性方程组的不完全LU分解预处理广义共轭残量（GCR）迭代方法应用于核反应堆中子时空动力学计算，结合仿真计算过程的连续计算处理，实现了模块式高温气冷堆中子时空动力学的三维实时仿真。 * 谢谢大家！ * * * * * * * * * * * * * * * * * * * * * 各能群平均注量率随时间变化 单根控制棒移动 反应性 相对功率 物理热工耦合 堆内材料物性参数确定 热工反馈的集总参数法模型 计算结果与分析 主要设计参数 （HTR-PM） 参数 设计值 热功率 250MW 堆芯平均功率密度 3MW·m-3 堆芯最大功率密度 7.5 MW·m-3 堆芯直径 3m 堆芯等效高度 11m 堆芯燃料球数 400000 单球平均功率 0.56kW 单球最大功率 1.4kW 一回路压力 7Mpa 堆芯压力降 68kPa 氦气质量流量 96kg·s-3 一回路氦气入出口温度 250/750℃ 最高燃料温度 950℃ 堆芯出口最高氦气温度 850℃ 燃料元件最高温度限值 1620℃ 热工反馈的集总参数法模型 将反应堆内材料的热阻、热容量分别按时间和空间取平均值, 并假定每个量的平均值集中在几何形状的中心, 略去燃料元件的轴向导热。这就是反应堆动态分析中所使用的集总参数模型。 应用集总参数法，将模块式高温气冷堆堆芯简化为UO2燃料颗粒、等效慢化剂和反射层3部分，如右图所示。根据彼此之间的能量交换关系建立各区温度随时间变化的微分方程组。 能量守恒方程 参数说明 初始状态的温度场分布确定 计算结果与分析 反应性扰动 冷却剂丧失 反应性扰动 各物理量动态模拟曲线 各区温度动态变化模拟曲线 冷却剂丧失 堆内温度动态变化模拟曲线 堆芯相对功率P随时间变化 （r, z）网格上第1群中子注量率随时间变化 （r, z）网格上第2群中子注量率随时间变化 （r, z）网格上第3群中子注量率随时间变化 （r, z）网格上第4群中子注量率