基于通量展开节块法的六角形时空动力学方程数值解.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 基于通量展开节块法的六角形时．空动力学 方程数值解 夏榜样1，谢仲生1，成春雨2，姚栋2 (1西安交通大学能源与动力工程学院，710049；成都核动力研究设计院，610041) 摘要；提出了一种基于节块内瞬态中子通量展开的六角形几何时．空动力学方程数值解法。在该方法中，各 群中子通量分布用解析基函数和二阶正交多项式近似展开，而包含各组缓发中子先驱核浓度的固定源项则利 用多项式进行近似。面平均偏流及其一次矩作为节块之间的耦合条件，不但明显改善了节块耦合关系，而且 使得响应矩阵技术比较容易地应用于迭代求解过程．对二维、三维基准问题计算表明，该方法能高效：准确 的给出各时间步内的堆芯总功率和节块功率分布． 1前言 反应堆安全事故分析往往需要求解多维多群时．空中子动力学方程组，以确定堆内瞬态中 子密度的空间分布。在空间处理上，各种现代节块法得到了广泛的应用，尤其是应用于正方形 组件的方法已经相当成熟，但对于六角形组件应用的还比较少Il捌。最近，通量展开节块法【4】 成功的应用于稳态六角形几何中子扩散方程的求解，并取得较高的计算精度。与传统的节块法 相比，该方法利用解析基函数把节块内通量分布近似展开来求解稳态中子扩散方程，有效的避 免了应用横向积分技术产生的奇异项问题，提高了计算精度，同时计算结果可直接用于节块内 的精细功率恢复。为了将该方法应用于瞬态计算并保留上述优点，本文将节块内的瞬态中子通 量分布用解析基函数和在六角形几何节块内正交的二阶多项式近似展开，其中多项式部分用权 重法求解，而包含各组缓发中子的先驱核浓度的固定源项则利用二阶多项式进行近似。 2瞬态中子动力学方程的离散 在如图l所不的六角形几何节块内，多维多群时·空动力学方程组。司写成： 毒掣二VDs∽V呲沪∑“附力+妊峨(∽ (2·1) +z妇∑4c，(r，，)+(1一∥)如∑vE店．o)西g．(r，r) (2．2) v∑石．(f)①譬．(r，f)一^cj(r，f)，i=1，2，…，， 皇竺夸}盟=屈∑G 上式中所用的符号都是习惯表示方法。 354 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 图1．六角形节块及其表面编号 考虑到收敛性和稳定性的要求，对于式(2·1)中时间变量的离散，本文采用全隐向后Eulcr 格式来处理，在|L。时刻其时间导数项可近似为： IZJ —ig①g—(r,t)I，=竺墨坠!二竺墨坠二!! (2) 西 ”。 △f一 一-==。·_-________。-__-·__--____-____-·_一 实践证明，隐式差分格式对数值计算过程是稳定的15】。对于先驱核浓度方程(2．2)‘，本 文采用时间积分法嘲来处理·假设在一个时间步内，V∑唐(r)①g(r，r)是线性变化的，然后对 先驱核浓度方程式(2-2)在一个时间步内积分，可得到乙时刻先驱核浓度的近似表达式： G(r，乙)=q(r，in_l弦一^M+磁∑V∑J窖．(，川)o∥(r，，州) g’Il +E∑’，∑fg,以)①g．(r，f一) (3) 譬’-l 式中：碟和露为常量，其值取决于材料特性。 把式(2)、(3)代入式(卜1)中，可以得到时间变量离散tn时刻的时·空动力学方程： -D(t一)V2西(r，t一)+A(‘)m(r，，一)=S谚(r，，一)