龙格库塔方法在求解瞬态中子扩散方程中的应用.pdfVIP

第47卷第1期 原子能科学技术 V01．47 ，No ．1 2013年1月 At omi c Ener gy Sc i e nce and Technol ogy J an ．20 13 龙 格 库 塔 方 法 在 求 解 瞬 态 中 子 扩 散 方 程 中 的 应 用 赵文博1’2 ，姚 栋2 ，王 侃 1，胡永 明3 ( 1．清华大学工程物理系，北京100084；2 ．中国核动力研究设计 院核反应堆系统设计技术重点实验室 ，四川成都6 1004 1； 3 ．清华大学核 能与新 能源技术研究院，北京 100084) 摘要：本工作开发了N FMN- K程序求解瞬态中子扩散方程 。空间离散采用第2类边界条件节块格林函 数方法，时间离散分别选取向后欧拉格式和 四阶精度对角线隐式龙格库塔( DI RK) 格式。对DI RK格式采 用嵌人三阶精度格式估计截断误差实现时间步长自调节 。数值验证结果表 明，两种格式 的计算结果与参 考程序结果符合很好，对于剧烈瞬变情况，DI RK格式较向后欧拉格式更为精确、高效。 关键词：时空动力学；节块格林函数法；向后欧拉；龙格库塔方法 中图分类号 ：TL325 文献标 志码：A 文章编号：1000—693 1( 20 1 3) 01—00 89 —0 8 do i ：10 ．753 8 ／yz k．20 13．4 7．01 ．0 089 Appl i cat i on of Runge - Kut t a Met hod t o Sol ve Tr ans i ent Neut r on Di f f usi on Equ at i on ZHAOWen—b0 1～，YAO Don92 ，WAN Kanl ，HU Yong—mi n 93 ( 1．Depar t ment of Engi ne er i ng Phys i cs ，Ts i nghua Uni ver s i t y ，Bei j i ng 100084，Chi na ； 2．Sci e nc e an d Techn ol ogy on React or Syst em Des i gn Technol ogy Lab or a t or y ， Nucl ear Power I n st i t ut e of Ch i na ，Che ngdu 6 1004 1，Chi na ； 3．I nst i t ut e of Nucl ear and Ne w En er gy Tech nol ogy ，Tsi nghua Uni ver s i t y，Bei j i ng 100084 ，Chi na) Abs t r ac t ：N FMN—K c ode was deve l ope d t o sol ve t r ans i ent neut r on di f f usi o n equat i ons ． Nodal r e en ’S f u nct i on met h od bas ed on t he se c ond boundar y condi t i on was ut i l i ze d f or spat i al di s cr et enes s ．Backwar d Eul er ( BE) an d a f our t h—or de r acc ur a t e di agon al l y i mpl i ci t Runge—Ku