优化时间步长的数值方法解核反应堆点动态学方程 solving point reactor kinetic equations by time step-size adaptable numerical methods.pdfVIP

第28卷第2期 核动力工程 Vbl．28．No．2 O 2OO7 NuclearPower A O 7 年4月 Engineering pr．2 文章编号：0258．0926(2007)02-0008一05 优化时间步长的数值方法解核反应堆 点动态学方程 廖茶清 (德克萨斯A＆M大学核工程系，77843-3133) 摘要：通过分析步长与数值误差的关系，论证了调整和优化步长的必要性。从步长与误差的关系出发， 介绍了数值求解初值微分方程问题的两种优化步长的方法：两步计算法和嵌入式龙格一库塔法。在用隐式欧拉 法对反应堆点动态学方程进行数值求解时，采用两步计算法对步长进行优化，计算的结果表明控制误差对步 长及数值解的精确性有决定性的影响；通过求解反应堆点动态学方程，对M枷，AB自带的嵌入式龙格．库塔 法的微分方程求解函数ode23及ode45的精确性及使用情况进行了探讨。 关键词：核反应堆点动态学方程；优化时间步长；误差；欧拉法；嵌入式龙格-库塔法 中图分类号：TL333文献标识码：A 1 引 言 I．I 旦． f + C 心 乃 ∽ 基于不同的条件和假设，核反应堆点动态学 Kmetic 方程P砌①s(PonReactor Equations)可以 ll 汁 ，∑一以 力 得到简化，相应的建立起多种分析模型，进而获 学掣 等邯儿 0一～ ．卜．所百 得解析解；也可以通过数值的方法来求解。求解 方程组中，尸(O为功率(幅值)；p∽为反应性； PIⅨEs，优化时间步长是非常有必要的。步长对肛萨o．0065；么为平均中子代时间，s；^为用冼驱 数值解的精确性有很大的影响，且对数值计算速 核素衰变常数，s～；“D为，群先驱核素浓度(幅 度起决定性的作用，步长取得小，收敛所需的计 值)；历为，群缓发中子份额；卢为缓发中子有效份 算时间就长；步长取得大，收敛所需要的计算时 额。 间短。 方程组涉及到多群缓发中子，每群缓发中子 从误差与步长关系出发，本文详细推导和介 的行为特性又各不相同。必须做一些假设，相应 绍了两步计算法及龙格．库塔．芬尔格法，前者被 地建立起多种分析模型，并对这些模型和简化后 用来优化和调整隐式欧拉法求解PRKEs的步长。的方程进行研究，以获得不同的近似解【1～。本 文采用了欧拉方法和龙格．库塔法对其数值求解。 MATLAB本身就使用龙格．库塔．芬尔格法，即嵌 入式龙格．库塔法来解ODE问题的。通过求解 方程组(1)有如下的矩阵形式：y=4】，。若用隐 P砌@s，本文对不同阶数的龙格．库塔．芬尔格法的 式欧拉法进行数值求解，其迭代的矩阵型式为 精确性及使用条件进行了探讨。 E+1=(，一4办)_1E (2) 式中，J是单位矩阵；办是时间步长；变量矩阵】， 2核反应堆点动态学方程