反应堆物理数值计算课件

1、第五章核科学与技术学院1内容5.5 解析节块法5.5.1 AM法求一维中子扩散方程5.5.2 系数矩阵的推导5.5.3 AM法求多维中子扩散方程5.6 节块内精细中子通量密度分布的计算核科学与技术学院2先进节块法的共同特点1、横向积分 三维转化成一维求解，其他方向的影响用“横向泄漏”来描述，且将它作为附加源项来处理。2、表面流技术 先求系统中相邻节块表面上中子流所满足的空间耦合方程组，得到粗网格表面上的中子流；然后利用表面中子流计算节块内平均中子通量密度。3、变分原理或剩余权重法核科学与技术学院35.5 解析节块法核科学与技术学院45.5.1 AM法求一维中子扩散方程其中：核科学与技术学院5一

2、、处理方程求解析解一维稳态中子扩散方程核科学与技术学院6定义则（5-132）可表示成（5-136）对于均匀核参数区域 ，方程（5-136）有以下解析解： 核科学与技术学院7由 到 积分上式，并除以 得到（5-137）（5-138）（5-140）（5-139）此处同理得核科学与技术学院8（5-144）利用三角恒等式处理并化简方程（5-138）、（5-140）得因而（5-145）核科学与技术学院9（5-146）式中扩散方程积分形式将（5-145）代入（5-146）得5.5.2 系数矩阵的推导 对最普遍的双群进行推导 双群方程（5-132）可表示成核科学与技术学院10（5-148）方程普遍解（5-1

3、49）式中（5-150）（5-151）r, s 为双群中子通量密度耦合系数，核科学与技术学院11（5-152）（5-153）通过（5-149）式，可导出净中子流表达式把（5-149）和（5-154）归并成矩阵形式可得（5-155）式核科学与技术学院12从而得（5-158）（5-159）核科学与技术学院13则（5-159）式可表示为（5-161）（5-162）核科学与技术学院14比较上述两式，可知（5-163）上式结合双曲函数定义及（5-135）式通过更复杂的运算，由（5-138）及（5-140）可得（5-145）式，其系数矩阵为核科学与技术学院15对单群情况结果进一步化简为5.5.3 AM法求

4、多维中子扩散方程一、横向积分处理方程核科学与技术学院16多维稳态中子扩散方程对于直角坐标系中，均匀长方体节块（l,m,n），三维中子扩散方程可表示为：核科学与技术学院17对上式应用横向积分，便得设以 表示v方向的中子泄漏：（5-174）（5-175）这是v方向上，节块表面上的平均中子泄漏。核科学与技术学院18令则同理可得v、w方向上类似的横向积分中子通量密度方程。（5-177）式是非齐次方程，对泄漏项源，采用二次函数逼近。二、求方程解析解并推导节块平均通量密度和表面平均净泄漏之间的关系式核科学与技术学院19上式表明，u方向上净中子泄漏和u方向上相邻三个节块平均中子通量密度以及五个相邻节块的横向

5、泄漏相耦合。核科学与技术学院20在节块（l,m,n）内，对中子扩散方程积分得将（5-180）式代入上式并变换得核科学与技术学院21基本方程（5-184）中平均通量密度项间耦合很弱而主要是和表面平均中子泄漏项发生耦合。这与实际大型反应堆不符，若对方程（5-184）作变换得以下形式核科学与技术学院22将大大加强节块平均中子通量密度项间耦合，从而有利于数值方法求解。5.6 节块内精细中子通量密度分布的计算 先进节块法不能提供节块内功率的精细分布，以计算功率不均匀系数，确定热点及燃料棒的最大辐照点。 计算精细中子通量密度分布的两种方法：1、调制法 2、嵌入法核科学与技术学院23高阶内插法求节块内光滑中子通量密度近似分布核科学与技术学院24核科学与技术学院25角点值的确定 每个角点四个相邻节块，每个节块的中子通量密度分布均用二阶多项式近似。（5-192）四个节块共32个系数，需32个方程。核科学与技术学院26四个相邻节块平均中子通量密度，可建立四个方程。节块间平均偏通通量密，可建立四个方程。交界面上平均偏中子通量密度对应的净中子流，确定四个方程核科学与技术学院27节块界面偏中子通量连续，确定四个方程交界面上偏净中子流连续，确定四