氙动态分析的几种数值解法研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 氙动态分析的几种数值解法 李富 韩松周旭华 清华大学核能与新能源技术研究院，北京，100084 摘要：氤动态过程是影响反应堆运行特性的一个重要动态环节，可用于反应堆运行监测、负荷跟随、优化 控制等方面．氤动态过程可以通过多种方法进行数值模拟、分析，如时间差分法、准稳态龙格库塔法等．本 文实现了两种数值解法。并针对一个压水堆模型，比较了不同的方法的优缺点，模拟了氙动态行为。分析了 氤振荡发生的条件，为进一步研究氙动态建立了坚实的基础． 关键词：氤动态；微分方程；扩散方程；本征值方程 1引言 氤动态是反应堆动态的一个重要方面。氙动态会影响反应堆的稳定性、安全裕厦和反应 堆的控制性能。本课题组试图对氙动态进行建模，预测氙动态，通过芯外探测器读数监视氤 动态，并发展控制氙动态的算法，进而用于反应堆运行监测、负荷跟随、优化控制等方面。 但这些工作的基础和前提是通过数值化方法把氙动态进行分析。 氙动态本身可以通过下面微分方程进行描述，这里先忽略温度反馈和人为调节等因素：1 专等=B俨吼一z嵋。g+薹∑g··暑①譬·+磊薹舛一r，S·①S--听培廊暑c·， 詈吖，丢∑“①g’一乃， ’(2) 詈=h∑∑几一譬·+乃，一以x—x∑％心一暑· (3) vI t-i|-1 式中g代表中子能群，札，D暑，z¨，zg，z，暑，听’暑和y分别代表中子速度、扩散系数、 移出截面、裂变谱、裂变截面、Xe-135的微观吸收截面、每裂变释放的中子数；G，71，h， ①g，，，X代表中子通量、1-135和Xe-135的原子密度；Z譬·．譬代表从g。到g群的散射截面． 上式有一定特殊性·首先式中存在昙。，昙，，旦atO／ Ot x，是一个动态方程；由于存在扩散算 符D。V2①。，又是一个空间扩散的Boltzmann方程：中子存在着时间变化项与空间扩散项之 间的竞争问题，且式1中存在非线性乘积项xo，，因此是一个非线性、关于时间和空间都 动态变化的动态方程． 如果再考虑在反应堆实际运行过程中，反应堆实际上维持在一种临界或稍微次临界的状 态，如果出现氙振荡情况，反应堆操纵员还会通过移动控制棒或调节硼浓度来使反应堆维持 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 在一个总体稳定的状态，同时氙动态还与温度反馈等效应互相作用，反应堆实际动态过程还 与上式的理论解不完全相同，因此氙动态的求解过程是一个很复杂的过程。 本文主要讨论上式的几种数值解法。基于这些数值解法，一则可以迸一步研究在温度反 馈和操纵员手动干预下的氙动态行为；--贝,J可以研究和预测反应堆的各种氙动态特性，进而 研究如何用芯外探测器探测氙动态和实时重构反应堆内的功率分布；三则可以研究监测和控 制氤振荡的算法。这些都是本课题组的长期研究方向之一。 2数值解法 式(1)彤)是关于时间和空间的方程。关于空间扩散的Boltzmann方程的求解有成熟的方 法，关于时间的一阶微分方程的求解也有成熟的方法，但如何同时求解并没有定法。本课题 组试验的解法的关键是分离时间与空间。经过分析，氙动态过程本身是一个缓慢的过程(时间 常数在小时量级)，即时间项是一个慢过程；而中子扩散过程是一个较快的过程．根据边界条 件和物质组成能很快达到空间平衡，即扩散项是一个主要的调节中子通量分布的因素。因此， 求解的关键是把时间项处理成空间分布的源项，再通过特定边界条件下扩散方程求解，求得 不同时间点的空间分布。不同的组合方式下具体解法也略为不同． 在通常的反应堆分析程序中，采用本征值方程求解中子方程。再依据中子通量分布根据 式(2)和式(3)求下一时间步的氤、碘分布，在已知氙、碘浓度分布的情况下又用本征值方程求 新时刻点的中子通量分布。这相当于一阶显式差分，要求时间步长不能太长，处理相对比较 粗糙。 本文试验了如下三种解法。一种是直接把时间微分项变成差分项，把整个方程化成扩散 方程的固定源问题：第二种解法是直接求解关于时间的微分方程，把扩散项按上一时刻的分 布来求解各区域的泄漏；第三种方法是只求解Xe和I的时间微分方程，把关于通量的时间和 空间的方程转化