解析节块法收敛性问题初步研讨.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 解析节块法收敛性问题初步研究 李志勇1，张少泓1，赵荣安2，2 l上海交通大学核科学与工程学院，上海，200030 2西屋电气公司堆芯技术部 摘要：中子扩散方程有限差分方法，已在理论上证明数值迭代收敛性．但节块法一直未在理论上证明其数 值收敛性．本文使用解析节块法(ACM)，构造一维单群情况下特殊计算条件，出现迭代不收敛现象，在某些 特定情况下，则收敛到高阶谐波。由于解析节块法(ACM)迭代过程中，迭代矩阵引入迭代变量七咖使其迭 代过程为非线性。迭代矩阵失去了有限差分迭代矩阵所具正定等性质，因而出现不收敛或收敛到高阶谐波的 现象． 关键词，解析节块法；收敛性；高阶谐波 1引言 粗网格差分方法(CMFD)，显式表示耦合系数，建立仅仅包括节块平均通量作为未知量的 中子平衡方程。CMFD广泛应用于先进节块方法的加速计算。近年，西屋公司赵荣安博士提 出了具有良好理论基础的解析粗网格差分方法(ACMFD)【1，卫，进一步完善了有关粗网格差分 方法的理论。 本文参考ACMFD方法，编写了基于最简单情形(一维单群)解析节块法(ACM)计算 模块．由于在一维单群情况下，没有引入任何近似，因此从理论上说，数值解的结果应和理 论解完全一致，但本文的研究工作表明，在某些特定的计算条件下(如节块宽度很宽)，解析 节块法将不收敛，或出现收敛到谐波解的情况。本文专门对该现象进行了研究。 2一维单群ACM理论模型 对于单一节块，一维单群扩散方程有如下形式： 掣+七场似：0 (1) 当七20时方程(1)的解析解为： 咖例=Asinkx+Bcoskx (2) 七：0时的解为指数函数，形式类似。利用式(2)EP可以获得节块左、右表面净中子流(．，工、，)、 表面平均中子通量(扩、o^)和节块平均中子通量(丕)的表达式，从而建立起净中子流、面通 量和节块平均通量之间如下的关系(假设节块宽度为口)： 矿t∥芸+c，_‘J， (3) 作者简介t李志勇．上海交通大学硕士研究生。目前在上海核工程研究设计院工作． 75 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 咖鼻l≈7茜+C，石 (4) 其中：Cj=百tan(ak／2)。c，-ak／sinak 应用相邻节块间的连续性条件．可得： 一12．rc．j瓦1+c。|七Col石n+l--c霸@ 或重新改写上式为： ．，一一，D．L+，；．“一D：石．．， (6) 称上式中的D正、D置分别为节块左、右表面的等效扩散系数(EDCs)： D：=2fc：去+co|七1．‘c{ m D：t2【C{去+col七∥c： 嘞 可以看出，等效扩散系数(EDCs)不仅仅与节块几何和材料参数有关，还与方程(1)中的参 数k有关．针对本文所研究的单群一维均匀反应堆问题，参数k可写为： 。 七；仁*屉∥一砂舱／L (9) 其中，L为反应堆的扩散长度。 此外，一维单群中子平衡方程为： (J：一JD／o．+zi．=Q 00) 式(6)代入式00)，就得到如下矩阵形式，