自然循环向强迫循环过渡时的反应性变化 transitional reactivity variation from natural to forced circulation.pdfVIP

第29卷第3期 核动力工程 Vbl．29．No．3 0 2O08 年6月 NuclearPower Jun．208 Engineering 文章编号：0258—0926(2008)03—0018-06 自然循环向强迫循环过渡时的反应性变化 姚青生1，王章立1,2，韩伟实1，贺克羽1 (1．哈尔滨工程大学核科学与技术学院，哈尔滨，150001；2．西安交通大学核科学与技术学院。西安，710049) 摘要：在Matlab平台中对反应堆及主冷却剂系统建立数学模型，研究自然循环向强迫循环过渡过程中 反应性的变化规律。计算结果表明：过渡过程中流龟按双曲正切规律变化，时间系数k和流量增加倍数肌 ／wo决定了流量的稳定时问和梯度；在反应堆没计中选择合适的k值，可有效改善反应堆周期；尽量降低过 渡过程的流量变化速度，可以减小反应性峰值，有助于提高反应堆的安全性。 关键词：自然循环；强迫循环；过渡过程；反应性变化；Matlab 中图分类号：TL329+．2文献标识码：A 1 前言 中子的点堆模型【1-3】： 核动力装置自然循环向强迫循环过渡时，反 应堆堆芯物理与一回路热工水力的耦合效应非常 o(1) 明显。由于主冷却剂流量突然增加以及蒸汽发生 簧2Tp(0-p相+酚6一 器热交换状态的改变，主冷却剂温度和燃料棒温 掣=鲁刀¨qo) 度也会随之变化，同时反应堆的负温度反应性效 式中，n(O为堆内平均中子密度，m一；Ci(力为第 应又将导致堆芯反应性的波动，甚至可能在瞬间 f组缓发中子先驱核平均浓度，m～；p为反应性； 产生很短的反应堆性周期，对反应堆的安全运行 允，为第f组缓发中子先驱核的衰变常数；A为瞬 十分不利。因此，研究自然循环向强迫循环过渡 发中子每代时间，s；卢为有效缓发中子份额；∥ 时反应性的变化具有重要意义。 ，为第i组有效缓发中子份额。 本文采用三阶Hermite插值多项式方法求解 由于问题本身的物理性质参数决定了式(1) 中子点堆方程，应用四阶龙格．库塔法求解主冷却 在求解过程中需要很小的步长才能够稳定，使得 剂传热方程，采用一维流动模型研究主冷却剂流 应用一些常规的数值方法求解该方程组会遇到严 量的变化规律，并在Matlab平台上编制了计算程 重的困难【4曲】。本文采用Hermite插值多项式法求 序。通过研究自然循环向强迫循环过渡过程中反 解点堆方程，能够很好地解决点堆方程的收敛问 应性和反应堆周期的变化规律，讨论影响反应性 题，而且中子动力学计算的时间步长与热工水力 和反应堆周期的边界条件，从而在反应堆设计中 计算相适应，可以节约计算时间【7J。 减小反应性峰值，改善反应堆周期，提高反应堆 将时间离散化，令：h=t2一tl，r=(t-t1)／h，方 的安全性。 程r1)可改写为： 2基本模型 2．1 中子动力学模型 (2) 言ldn=4(f如G)+莩以c，G) 反应堆中子动力学主要是讨论反应