基于adi的一种堆芯物理实时仿真算法.pdfVIP

第32卷第3期 核动力工程 VbI．32．NO．3 20l l 0 年6月 NuclearPower JUn．2ll Engineering 文章编号：0258．0926(201 1)03-0130·04 基于AD I的一种堆芯物理实时仿真算法 T／J,Jll1，路海晋2，赵强1 (1．哈尔滨工程大学核安全与仿真技术国防霞点学科实验室．哈尔滨，150001； 2．哈尔滨工程大学核科学与技术学院，哈尔滨，150001) 摘要：针对传统的交替方向隐式(ADI)方法求解中子扩散方程，在仿真中难以满足实时性要求这一问 题，提出了以组件为单位划分网格以保证计算速度，并通过对ADI方法进行调整，引入调整因子使扩散方 程强制守恒。为了验证该方法的可行性，以秦山二期核电厂为对象建立了仿真模型，结果表明，调整后的 ADI方法能够满足计算精度和速度要求，可以用于反应堆堆芯物理实时仿真。本文结合在线校正技术对仿真 结果进行在线校正，实现了在线仿真。 关键词：堆芯物理；实时仿真；ADI；Douglas格式；在线仿真 中图分类号：TL329．2文献标识码：A 1 引 言 上丝掣：VD：o，，)v妒：㈥一 在核电厂仿真中，堆芯物理仿真是全电厂仿 (2) 真的重要组成部分。目前，堆芯物理仿真普遍基 乞。：◇，f溉◇，f)+五：O，f砑p，f) 于中子扩散模型，对此模型的求解方法有节块展 ’ 坠型一4c，㈤+ 优 一 开法、节块格林函数法等粗网格计算方法。这些 粗网计算方法可以很好地满足计算精度要求，但 局【v。。A◇，f)+v：毛：如◇，t)l(3) 式中，魂、九为快中子、热中子注量率，cm～·s～； 耗费的时间长。由于实时仿真要求每一步的计算 时间小于仿真时间，因此，在对堆芯物理进行实 Dl、岛为快中子、热中子扩散系数，cm；乞．1、 时仿真时，节块展开法、节块格林函数法等方法 乞．2为快中子、热中子吸收截面，cm～；五2为转 较难满足实时性要求。本文以交替方向隐式 移截面，cm～；邑l、五2为快中子、热中子裂变 截面，cm-1；Cf为缓发中子先驱核密度，cm。3；以 (ADI)方法为基础，以调整后的ADI方法求解 中子扩散方程。 为缓发中子先驱核衰变常数，s～；s(，．，，)为外加 中子源，cm一·s-1；vl、v2为每次裂变释放平均快、 2堆芯物理仿真模型 热中子数；q、02为快、热中子平均速度；屈为 2．1扩散方程 各组缓发中子份额；∥为缓发中子总份额。 反应堆堆芯物理数学模型采用三维两群带 2．2改进的准静态近似 六组缓发中子的扩散模型。其表达