非齐次燃耗方程数值解法付元光邓力李刚Preliminarystudy-物理学报.pdf

非齐次燃耗方程数值解法 付元光 邓力 李刚 Preliminarystudyonnumericalsolverofinhomogeneous burnup equations FuYuan-Guang DengLi LiGang 引用信息Citation: ActaPhysicaSinica, ,172802(2018) DOI: 10.7498/aps.67 在线阅读Viewonline: /10.7498/aps.67 当期内容Viewtableofcontents: /CN/Y2018/V67/I17 您可能感兴趣的其他文章 Articlesyoumaybeinterestedin 基于节块展开法的Jacobian-FreeNewtonKrylov联立求解物理-热工耦合问题 Jacobian-FreeNewton-Krylovbasedonnodalexpansionmethodforneutronic-thermalhydrauliccoupling problem 物理学报.2016,65(9): 092801 /10.7498/aps.65.092801 BEAVRS基准模型热零功率状态的JMCT分析 JMCTMonteCarloanalysisofBEAVRSbenchmark: hot zeropower results 物理学报.2016,65(5): 052801 /10.7498/aps.65.052801 反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进 Modifieduniform-fission-sitealgorithminMonteCarlosimulation of reactor criticality problem 物理学报.2015,64(5): 052801 /10.7498/aps.64.052801 金属铁中空位团簇演化行为的相场研究 Phase-fieldmodelingofvacancyclusterevolutioninFe 物理学报.2013,62(18): 182801 /10.7498/aps.62.182801 Nd: YAG激光烧蚀裂解加工技术模拟分析与实验研究 FiniteelementanalysisandexperimentalstudiesonfracturesplittingprocessingbyNd: YAGlaserablation 物理学报.2012,61(9): 092801 /10.7498/aps.61.092801 物理学报 Acta Phys. Sin. Vol. 67, No. 17 (2018) 172802 非齐次燃耗方程数值解法 付元光 邓力 李刚 1) (北京应用物理与计算数学研究所, 北京 100088) 2) (中物院高性能数值模拟软件中心, 北京 100088) ( 2017 年12 月14 日收到; 2018 年4 月13 日收到修改稿) 非齐次燃耗方程常用于描述具有显著核素迁移效应的核能系统中核素含量随时间的变化规律. 国内外许 多燃耗计算程序无法求解方程非齐次项含时的情况. 本文在方程非齐次项能够被有限阶关于时间的泰勒展开 逼近这一前提下, 研究了非齐次项含时情况下方程的两种解法. 首先通过Laplace 变换推导出了方程基于线 性子链方法的解析解形式, 然后使用Carathéodory-Fejér 方法计算出了方程矩阵级数解的近最佳Chebyshev 有理逼近式. 将两种方法在燃耗计算程序JBURN 中实现, 并进行了数值计算, 绝大部分计算结果符合很好, 部分结果在较长有效数字内仍能保持一致, 验证了方法的正确性和精度. 同时为求解具有其他非齐次项形式 的燃耗方程提供了一种思路. 关键词: 非齐次燃耗方程, Laplace 变换, Carathéodory-Fejér 方法, 近最佳Chebyshev 有理逼近 PACS: 28.52.Av, 28.41.Vx, 28.41.Ak DOI: 10.7498/aps.67