核反应堆物理课设

程设计（综合实验）报告1课程设计(综合实验)报告名称：核反应堆物理分析题目：利用双群理论求解堆芯参数院系：1111111111111班级：111111111111111学号：111111111111学生姓名：11111111111指导教师：111111111设计周数：11111111成绩：程设计（综合实验）报告2一、课程设计(综合实验)的目的与要求1课程设计的要求课程设计是重要的实践教学环节。它是根据教学计划的要求，在教师指导下对学生进行的阶段基础或专业技术训练，该实践环节着重培养学生综合分析和解决实际问题的方法与能力，实现由知识向智能的初步转化；是对前期理论与实践教学效果的检验，也是对学生综合分析能力与独立工作能力的核反应堆物理分析课程设计的目的是对理论课上学过的理论知识进行践应用，进而加深对前期理论知识的学习，是对学生综合运用核反应堆物理分析知识和思想方法的综合检验过程。2课程设计要求核反应堆物理分析课程设计的要求有如下几点：（1）学生必须修完课程设计的先修课程，才有资格做课程设计。（2）明确课程设计的目的和重要性，认真领会课程设计的题目，读懂课程设计指导书的要求，学会设计的基本方法与步骤，积极认真地做好准备工作。（3）在课程设计中，学会如何运用前修知识与收集、归纳相关资料解决具体问题的方法。（4）学生必须在指导教师指导下独立完成设计任务，严禁抄袭、找人代做等，一经发现成绩记零分，按考试作弊处理。（5）课程设计报告学校有统一格式，学生必须按照此格式填写，说明书、计算书要求简洁、通顺、计算正确，图纸表达内容完整、清楚、规范。二、设计（实验）正文1.课程设计题目：二维XY几何模型，中间是堆芯（裂变区）：-50cmx+50cm,-55cmy+55cm,四周是15cm厚反射层，堆芯及反射层的群常数见下表：表1堆芯各种参数程设计（综合实验）报告3群常数堆芯反射层第1群第2群第1群第2群g1.00.01.00.0)(1cmvf0.00750.16510.00.0)(1cma0.01210.1210.00040.020)(121cms0.0241-0.0493-)(cmD1.260.3541.1300.166利用数值求解下面问题，计算其快群和热群中子通量密度空间分布及有效增殖因子。提示：（1）采用内外迭代方法求解；（2）给出源迭代过程；（3）给出内迭代过程并给出堆芯双群方程、反射层双群方程以及边界条件，详见课本第五章；（4）建立数学模型，由于问题具有对称性，因此，只计算1/4堆芯即可。（5）可用赛德尔迭代法求解；（详见课本第五章）；（6）给出计算流程（7）编写程序，给出结果，并用OREGIN画图把结果表现出来；（8）结果的分析讨论；（9）写出报告。2.1差分方程的建立为简化计算，进行以下假设：（1）二维XY几何模型是一个整体，即堆芯和反射层交界处的中子通量密度和中子流连续；（2）认为该问题中反射层厚度已包含外推边界；（3）认为该问题中堆芯和反射层对称分布。在带反射层的多区双群扩散方程的数值解法中，通常是通过外-内迭代过程来求解。外迭代是一个通程设计（综合实验）报告4过迭代求特征值的迭代过程，所以又成为源迭代。内迭代则是对源迭代中出现的扩散方程进行具体的数值求解。计算中若不考虑中子自低能群的向上散射，多群扩散方程可以写为：)()(rr21grrr''''11-n1-neffg11')(g'nggggtgrngnggrnggrrQQSGSDggGgfggngg）（）（）（，，）（）（，）（）（）（，，，），（）（）（（1）略取上下标，该方程可以写成以下标准形式：)r(rrrSD）（）（（2）本设计讨二维（x，y）情况下的差分方程。对所取平面首先用x=0x，1x，…，ix，…，nxy=0y，1y，…，jy，…，my直线族把平面分成许多矩形网格（如图一）。交点（ix，iy）称为节点或网点，如图一中共有（N+1）(M+1)个节点，iiiiiiyyyxxx11,称为网距。图1网点的划分0x1xixnx程设计（综合实验）报告5图2（i，j）节点示意图按照数值解法的相关要求和定义，讨论如图二所示的（i，j）节点，根据二维几何情况下扩散方程导出（i，j）节点的差分方程为　　　　　　　　式中系数为jijijijijijijijijijijiSedcba,,,1,,,1,,1,1,,（3）式中系数为1411,1211,/21/)(21aijjjijiijixyDyDbyxDxDjijiijjjiyxDDdxyDyD/x21/21c31-i4,312,jijijijijirjirjirjirjidcbaxxyxyxyxe,,,,114,3,12,111,,41)(4111411312111,jijijijijiyxSSyxSyxSS（4）对于边界上的点需要根据边界条件来确定，本设计由于反射层已经包括了外推距离，所以外围上的点的中子通量密度均设为0。程设计（综合实验）报告6上述（3）和（4）式便是求解单群中子扩散的差分方程组，其系数jia,，jib,，…，jie,等均可以事先计算求得，对每次外迭代，其右端源项jiS,也是已知的。因此，它是一个含有ji,的代数方程组，对于二维问题，其系数矩阵是一个五队角的矩阵，可以很方便的应用通常的线性方程组的求解方法进行求解。2.2源迭代过程根据连续性的假设条件，堆芯的双群方程与放射层得双群方程可以归结为下面两个方程：快群中子扩散方程：rrkrrDffeffsa221111211,11（5）热群中子扩散方程：)(12122,22rrrDsa（6）设rrrQff2211)(（7）式（5）是一个关于)(1r的齐次方程组。式（5）除了中子通量密度）（r1外，rQ和effk也是未知参数，并且这两个参数又与中子通量密度有关。因此我们假定一个初始的裂变源分布)()0(rQ，认为在整个堆芯都等于1或某个给定常数，并猜测一个初始的)0(effk值，同时把)0()0()(effkrQ作为初始迭代源项代入到方程（5）的右端中。这样，方程（1）的右端便是一个已知项，方程（1）便成为一个普通的非齐次方程组，它可以用通常的数值方法求解。假定由它求出的中子通量密度分布为r11，把r11代入方程（6），由方程（6）数值求解求出）（）（r12，将求得的r11、）（）（r12带入（7）式便可求得第二次迭代裂变中子源分布)()1(rQ，并由它可以求得有效增殖因数的新的估计值)1(effk和第二代迭代源项)1()1()(effkrQ来，……。依此类推，逐次的迭代下去，对于第n次迭代计算有)()()1()1(1)(1211,)(11rQrrDnneffnsan）（（8）式中：)()()()()()1(22)1(11)1(rrrQnfnfn（9）根据effk的物理意义（它等于新生一代的裂变总中子数与上一代源中子总数之比）可以由下式算出)1(nk的估计值程设计（综合实验）报告7VnneffVnneffdVrQkdVrQk)(1)()2()2()1()1(（10）方程（3），（4）和（5）便是源迭代法的计算公式。从物理上看，上面所确定的迭代过程相当于将反应堆内的中子数人为地区分为不同的“代”，每一“代”中子有裂变形成之时作为起始。每一次源迭代则相当于老一代源中子通过扩散、慢化、俘获、裂变，最后又形成新一代源中子的过程。式（5）中的)1(neffk显然反映了中子每经过一代的倍增。因而可以预期，经过了足够的循环代数后，任意的初始中子源)()0(rQ必将在堆内形成一个稳定的最终中子源分布)(rQ。多群扩散方程的特征值问题在数学上也已被详细研究过，并证明了上述迭代过程的收敛性。2.3内迭代过程对于一次内迭代过程，假定一个初始的裂变源分布)0rQ（）（及一个初始的)0(effk数值，得如下计算步骤：（1）芯部中子通量密度的计算芯部快群（第1群）芯部快群中子扩散方程为：rrkrrDffeffsa221111211,11内迭代所需群常数为：1DD211,sarrrrSffeff22111内迭代相应的系数为：)1,1()1,1()1,1()1,1()1,1()1,1()1,1()1,1(214,213,212,211,jijijijijijijijisarsarsarsar程设计（综合实验）报告8);1,1();1,1();1,1();1,1(4321jiDDjiDDjiDDjiDD;)1,1()1,1()()1,1()1,1()(;)1,1()1,1()()1,1()1,1()(;)1,1()1,1()()1,1()1,1()(;)1,1()1,1()()1,1()1,1()(221114221113221112221111jijijijikSjijijijikSjijijijikSjijijijikSffeffffeffffeffffeff);(21;21;21);(2143,32,41,21,DDdDDcDDbDDajijijijijijijijirrrrjidcbae,,,,4,3,2,1,,414321ji41SSSSS，由节点（i，j）处的差分方程为：jijijijijijijijijijijiSedcba,,,1,,,1,,1,1,,最后得相应点处的快中子通量密度）（r1：)(1,1,,,1,,1,1,,,,jijijijijijijijijijijiSdcbae边界部分的点的中子通量密度认为为0，不必计算。芯部热群（第2群）芯部热群中子扩散方程为：)()()r(12122a22rrDs，内迭代所需群常数为：程设计（综合实验）报告9rrSDDsar1212,2内迭代相应的系数为：);1,1();1,1();1,1();1,1(4321jiDDjiDDjiDDjiDD);1,1();1,1();1,1();1,1(4,3,2,1,jijijijiarararar1,1)1,1(1,1)1,1(1,1)1,1(1,1)1,1(121s4121s3121s2121s1jijiSjijiSjijiSjijiS);(21;21;21);(2143,32,41,21,DDdDDcDDbDDajijijijijijijijirrrrjidcbae,,,,4,3,2,1,,414321ji,41SSSSS由节点（i，j）处的差分方程为：jijijijijijijijijijijiSedcba,,,1,,,1,,1,1,,