安全评价软件Ecolego解法器选取研究

辐射防护 ›› 2019, Vol. 39 ›› Issue (2): 95-104.

安全评价软件Ecolego解法器选取研究

刘兴伟, 王旭宏, 吕涛, 李星宇

中国核电工程有限公司,北京 100840

收稿日期:2018-08-30 出版日期:2019-03-20 发布日期:2019-11-08
作者简介:刘兴伟(1990—),男,2018年博士毕业于中国地质大学(北京)水文地质学专业。E-mail:liuxingwei0912@163.com

Study on selection of safety assessment software Ecolego solver

LIU Xingwei, WANG Xuhong, LV Tao, LI Xingyu

China Nuclear Power Engineering Co. Ltd., Beijing 100840

Received:2018-08-30 Online:2019-03-20 Published:2019-11-08

摘要/Abstract

摘要： Ecolego软件已广泛应用于放射性废物处置场的安全评价中,其提供了十余种解法器供用户使用,但不同解法器的性能相差较大,在安全评价中如何选取适宜的解法器给用户造成了困扰。本研究在对各解法器特点进行论述的基础上,利用Ecolego 6.5提供的解法器对瑞典高放废物处置安全评价模型进行测试,结果表明NDF解法器在地质圈及总模型计算耗时最少,BDF及RADAU5解法器分别在近场模型及生物圈模型耗时最少;Rosenbrock与TR-BDF2解法器计算结果稳定。综合计算耗时与求解精度考虑,在类似瑞典高放废物处置场地的安全评价模拟中,推荐使用TR-BDF2解法器。

关键词: Ecolego, 解法器选取, 放射性废物处置, 安全评价, 耗时分析, 求解精度

Abstract: Ecolego software has been widely used in safety assessment of radioactive waste disposal sites. The software provides lots of solvers, while the performance of different solvers varies greatly. How to select a suitable solver has caused troubles to users. On the basis of discussing the characteristics of each solver, the classical Swedish safety assessment model was tested with different solvers of Ecolego version 6.5. The results show that the NDF solver requires least calculating time in geosphere and whole models. BDF and RADAU5 solvers require the least time in near-field and biosphere models respectively. The Rosenbrock and TR-BDF2 solvers are accurate and stable. Considering the time-consumption and precision of calculation, the TR-BDF2 solver is recommended for safety assessment of high-level radioactive waste disposal sites which are similar to Sweden’s site.

Key words: Ecolego, solver selection, radioactive waste disposal, safety assessment, time-consumption analysis, resolution

中图分类号:

TL942+.1

刘兴伟, 王旭宏, 吕涛, 李星宇. 安全评价软件Ecolego解法器选取研究[J]. 辐射防护, 2019, 39(2): 95-104.

LIU Xingwei, WANG Xuhong, LV Tao, LI Xingyu. Study on selection of safety assessment software Ecolego solver[J]. RADIATION PROTECTION, 2019, 39(2): 95-104.

参考文献

[1] AVILA R, BROE R, PEREIRA A. Ecolego—A toolbox for radioecological risk assessment [C]//Proceedings of the International Conference on the Protection from the Effects of Ionizing Radiation, IAEA-CN-109/80. Stockholm: International Atomic Energy Agency. 2003: 229-232.
[2] 李星宇, 刘翔宇, 王旭宏, 等. 我国岩洞型低中放废物处置核素迁移计算库室模型研究 [J]. 工业建筑, 2018, 48(04): 18-22.
LI Xingyu, LIU Xiangyu, WANG Xuhong, et al. Research on nuclide transport compartment model for low and intermediate level waste disposal in rock cavern type [J]. Industrial Construction, 2018, 48(04):18-22.
[3] 赵杨军, 李洋, 顾志杰. 应用Ecolego软件计算低、中放固体废物处置场核素在地下水中的迁移 [J]. 辐射防护, 2013, 33(4): 230-234.
ZHAO Yangjun, LI Yang, GU Zhijie. Calculation of nuclides transfer at solid ILRW disposal facility using Ecolego software [J]. Radiation Protection, 2013, 33(4):230-234.
[4] Krisanangkura P, Itthipoonthanakorn T, Udomsomporn S. Environmental dose assessment using Ecolego: case study of soil from Japan [J]. Journal of Radioanalytical & Nuclear Chemistry, 2013, 297(3): 443-450.
[5] Lindgren M, Widen H. Discretization in COMP23 for SR97 [R]. Sweden: Swedish Nuclear Fuel and Waste Management Co., SR98-03, 1998.
[6] Maul P, Robinson P, Avila R, et al. AMBER and Ecolego intercomparisons using calculations from SR 97 [R]. Sweden: Swedish Nuclear Power Inspectorate, Ski Report 2003:28.
[7] Broed R, Ekstr M P A. Sensitivity analysis methods and a biosphere test case implemented in EIKOS [R]. Finland: Facilia AB, 2006.
[8] 东方. 浅析SIMULINK中的仿真算法 [J]. 轻工科技, 2012, 28(06): 68.
DONG Fang. Analysis of simulation algorithm in Simulink [J]. Light Industry Science and Technology, 2012, 28(06): 68.
[9] 戴嘉尊,邱建贤,编著.微分方程数值解法[M] 南京:东南大学出版社,2012:1-40.
DAI Jiazun, QUI Jianxian. Numerical solution of differential equations[M] Nanjing: Southeast University Press,2012:1-40.
[10] Bogacki P, Shampine L F. A 3(2) pair of Runge - Kutta formulas [J]. Applied Mathematics Letters, 1989, 2(4): 321-325.
[11] Shampine L F, Reichelt M W. The MATLAB ode suite[J]. SIAM Journal on Scientific Computing, 1997, 18(1): 1-22.
[12] Hairer E, Wanner G. Solving ordinary differential equations II: stiff and differential-algebraic problems [M]. Switzerland:Springer, 2008.
[13] Podhaisky H, Weiner R, Schmitt B A. Rosenbrock-type ‘Peer’two-step methods[J]. Applied Numerical Mathematics, 2005, 53(2-4): 409-420.
[14] Shampine L F, Reichelt M W, Kierzenka J A, et al. Solving Index-1 DAEs in MATLAB and Simulink [J]. Siam Review, 1999, 41(3): 538-552.
[15] Maria L, Kemakta K, Fredrik L, et al. SR 97 Radionuclide transport calculations [R]. Sweden: Swedish Nuclear Fuel and Waste Management Co, TR99-23, 1999.