[发明专利]一种Pin-by-Pin分析压水堆堆芯瞬态的方法

说明书

一种Pin‑by‑Pin分析压水堆瞬态的方法，步骤如下：1、计算堆芯初始状态，作为瞬态分析的初始点；2、对下一个时间步，考虑外部扰动信息，同时利用上一个时间步的中子通量密度更新截面；得到SP3稳态方程的系数和源项；3、利用指数函数展开节块方法空间上Pin‑by‑Pin数值求解时间离散的第t_n+1步的稳态SP3方程，迭代求解收敛后得到新时间步下栅元的零阶和二阶中子通量密度；4、更新缓发中子先驱核浓度；5、判断当前时间步是否为最后时间步，若不是，则重复2‑4；若是，则停止计算，输出瞬态参数。本发明相对于组件均匀化方法，不需要考虑功率重构引入误差；相对于扩散近似只能得到零阶中子通量密度，保留了更多的堆芯信息。同时求解时间并没有显著增加。

技术领域

本发明涉及核反应堆设计和反应堆物理计算领域，具体涉及一种Pin-by-Pin分析压水堆瞬态的方法。

背景技术

在反应堆运行过程中，需要实时掌控堆内的真实情况，以便于操作员可以及时进行反应性调节，避免事故的发生，或在事故发生后控制事故的范围和等级，保障工作人员与周围群众的生命财产安全。这需要研究一种实时且准确分析压水堆堆芯瞬态过程的方法，即中子动态学计算。堆芯与事件相关的现象有很多种，根据其过程时间常数的不同大致可以分成三类：短时间瞬态、中等时间瞬态和长时间瞬态。短时间瞬态持续时间通常为毫秒到秒量级，主要是由于人为地或者事故地引入反应性而导致的中子通量密度迅速改变，比如瞬态实验，启动，停堆过程等。中等时间瞬态通常为持续数小时或者数天的过程，主要是热中子堆中裂变产物(Xe135，Sm149等)的积累、消耗以及β衰变所导致的瞬态过程。长时间瞬态的时间跨度达到数个月甚至数年，主要是由于易裂变产物的积累和燃耗，和可裂变核素的积累、燃耗以及衰变。

用确定论方法直接求解中子输运方程是很耗时的，无法满足实时性要求，因此传统方法都是求解角度P1近似的中子扩散方程，空间上通常采用组件均匀化思想，先求得组件均匀化的中子通量密度，再通过功率重构得到组件内部的功率分布。但这样求得的中子通量密度和功率分布往往与实际的真实结果误差较大，特别是对于局部的非均匀效应。SP3近似方程是相对于中子扩散方程更加精确的中子输运方程形式，通过空间上Pin-by-Pin(以栅元为典型网格大小)求解SP3近似方程可以获得每个栅元内部的中子通量密度和功率分布。西安交通大学NECP实验室已经开发堆芯程序EFEN(指数函数展开节块方法)可以求解SP3近似方程。

因此，关键在于如何Pin-by-Pin求解瞬态的SP3近似方程，得到每个栅元内部的中子学信息，实时反映和分析反应堆运行瞬态过程中堆芯内部的真实情况。

发明内容

为了克服上述现有技术存在的问题，本发明的目的是为了给出一种Pin-by-Pin分析压水堆瞬态的方法，采用差分方法和积分方法将瞬态中子输运方程离散成稳态中子输运方程，再将它转化成稳态SP3近似方程，再在每个时间步下Pin-by-Pin迭代求解稳态方程，得到各个时间步的栅元内部的中子通量密度和功率分布，同时计算时间依旧满足实际运用需求。

为了达到上述目的，本发明采取了以下技术方案予以实施：

一种Pin-by-Pin分析压水堆瞬态的方法，步骤如下：

步骤1：计算堆芯初始状态，作为瞬态分析的初始点；得到每个栅元内部的零阶和二阶中子通量密度和修正后的多群微观总截面，微观散射截面和微观裂变截面；具体包括如下步骤：

1)读取输入卡片中堆芯几何参数进行几何建模，读取核素的原始宏观总截面，宏观散射截面，缓发中子份额，衰变常数，能群信息，中子能谱，单次裂变产生中子数，读取瞬态时间步长信息和扰动信息；

2)基于输入参数和堆芯几何计算初始的缓发中子先驱核浓度、SP3稳态方程的系数和源项；计算公式如下所示：

式中:

g'——碰撞前中子能群；