哈工程 核反应堆的核物理第9章 核反应堆动力学

第9章核反应堆动力学9.1缓发中子的作用尽管缓发中子所占份额很小，但这些缓发中子缓发的时间很长，它对核反应堆动态特性有重要的影响。假设所有裂变中子都是瞬发的。设t时刻堆内中子平均密度为n(t)，有效增殖系数为keff，瞬发中子平均寿命为l，得堆内中子密度的变化率：积分得：1()()effkdntntdtl01()exp()effkntntl反应堆周期：反应堆内中子密度变化e倍所需要的时间，也称为反应堆时间常数。如果只考虑瞬发中子，中子的平均寿命等于瞬发中子的慢化时间和热中子扩散时间之和，其数量级为10-4s。反应堆周期非常小。考虑缓发中子后，反应堆周期值将大大增加。考虑缓发中子后，裂变中子的平均寿命：1efflTk6611(1)()iiiiiilltllt9.2点堆中子动力学方程核反应堆动力学方程单群中子扩散方程中考虑缓发中子后的中子源项：将其代入（4-2），得：61(,)(1)(,)(,)aiiiSrtkrtCrt2611(,)(,)(,)(1)(,)(,)aaiiirtDrtrttkrtCrt先驱核通过衰变而消失，假定先驱核在衰变之前没有空间运动，得到：上两式即为核反应堆的时空动力学方程。点堆动力学方程假定当反应堆偏离临界状态不远时，空间形状函数满足：(,)(,)(,)iiaiiCrtkrtCrtt(,)()()(,)()()iiirtntrCrtCtgr22()()0rBr由以上几式得：假设与空间位置无关，则将（9-9）和（9-10）代入（9-8）得：6221()[(1)1]()()()iaiiigrdnkLBntCttr()/()igrr6eff1(1)1()()iiikdnntCttleff()()iiiidCkntCttl定义中子每代时间得点堆模型动态方程：这两个方程是相互耦合的一阶微分方程组。点堆模型的适用范围点堆动力学只适用于反应堆偏离临界状态不远和扰动不太大的问题。eff/lk61()()()iiidntntCtt()()iiiidCntCtt9.3阶跃扰动时点对模型动态方程的解不考虑反馈效应，ρ＝常数的情况下，点堆模型动态方程的解为：求解得：()tntAe()tiiCtCe()iiiCA6611111iiiiiilll方程组的解为：第i组先驱核的浓度写成确定Aj和Cij：7127012701()()jttttjjntnAeAeAenAe71()(0)jtiiijjCtCe0()(0)jiijjiiAnCC16021()iijjijjA在临界反应堆中引入阶跃反应性扰动后，中子密度开始时突然迅速变化，经过一小段时间后，开始以稳定的速度较缓慢地变化。9.4反应堆周期中子密度在一段时间的急剧变化后，其时间特性最终表现为：周期定义：中子密度按指数规律变化e倍所需的时间。上式定义的周期通常称为反应堆的稳定周期或渐进周期。用中子密度的相对变化率来直接定义反应堆周期：1/()()ttTntente或11T()/ntTdndt倍周期（倍增周期，Td）：堆内中子通量密度增长一倍所需的时间。反应堆周期根据定义可有反应性方程确定，将代入得：或dln20.693TTT1/Td0d()/exp(/)2nTnTT1iiiTTeff1iiilkTT不同反应性引入时反应堆的响应特性•当引入的反应性很小时得：可见，当引入很小的反应性时，反应堆周期与瞬发中子寿命无关，而与引入的反应性成正比，且取决于缓发中子寿命。1112l601111/iiil6110011[/]iiilTl•当引入的反应性很小时反应堆的响应特性主要决定于瞬发中子的每代时间。•当时，仅靠瞬发中子即可使反应堆保持临界，成为瞬发临界；当时，反应堆要达到临界需要缓发中子的贡献，反应堆的时间特性很大程度上由先驱核衰变的时间决定；当时，称为缓发临界，当时，称为瞬发超临界。011001T000瞬发临界条件：eff(1)1k或•当为很大的负反应性时，稳定周期T将接近于1/λ1，即约等于80s。9.5点堆动力学方程的近似解法单组缓发中子近似单组缓发中子的份额，先驱核密度和先驱核衰变常数：反应性方程简化为：666111;;/iiiiiiiCC2()0方程的解：在一定条件下，两个根简化为：点堆动态方程的解为：21,21[()()4]212,120120()[][exp()exp()]ttntnAeAentt常数缓发中子源近似研究控制棒插入时间t0以后很短的时间内中子密度对反应性的时间响应特性，由（9-29）得：反应堆动态方程为：当ρ等于常数ρ0的阶跃变化时，方程的解为：6010(0)iiiCn0()()(0)dntntndt0000()[exp()]nnt瞬跳近似认为中子的每代时间等于零，即瞬跳近似。点堆中子动力学方程在此近似下简化为：61()()()iiiNtct()(()())iiidctNtctdt9.6点堆动力学方程的数值解法点堆动力学方程的矩阵表示点堆动力学方程：初始条件：方程用矩阵表示：1()()()IiiidntntCtt()()iiiidCntCtt000()|,()|(0)titintnCtC()()()dytFtytt00()|tyty其中：11()/00()/00IiIItFt1()[(),(),,()]TIytntCtCt001[,(0),,(0)]TIynCC