蒙特卡罗自动建模中的空腔简化算法研究
|
摘要
考虑到粒子跟踪效率,在基于面描述的蒙特卡罗粒子输运模拟程序的计算模型中,除描述系统各部件的几何实体外,还需对非实体空间即空腔进行几何描述。由于空腔几何需严格填满实体间空隙,其建模复杂度通常较实体几何高,所以复杂系统常因空腔复杂度过高而严重制约蒙特卡罗建模和计算效率。本文研究了一种基于凸实体的自动空腔简化算法,该算法首先将栅元实体分解为凸实体集合,利用与子空间相交的凸实体的补集描述该子空间区域的空腔。该算法已在多物理耦合分析自动建模软件SuperMC/MCAM中实现,并使用国际热核聚变实验堆ITER中子学基准模型进行了测试。测试结果表明,相比基准算法,该算法可简化空腔表达式到1/60,提升转换效率25倍及生成模型的计算效率1.4倍。
Considering the transport efficiency,all the spaces including solids and void space need to be described for Monte Carlo transport codes.Since the void space needs to be described accurately,the void space modeling is more complicated than the solids modeling.For system with complex geometry,the void filling may generate too complicated Monte Carlo description for efficient modeling and calculation.Based on the multi-physics coupling analysis modeling program SuperMC/MCAM,an improved void filling algorithm based on convex solids was developed.This algorithm would test all the convex solids generated by decomposing the CAD solids with sub-spaces and generate the void description by complementary describing convex volumes intersected with the sub-spaces.This algorithm was tested with ITER reference neutronics model.The result shows that this algorithm simplifies the description of the void to 1/60,accelerates the void space modeling by 25 times and speeds up calculation by 1.4times compa-ring to benchmark algorithm.
Considering the transport efficiency,all the spaces including solids and void space need to be described for Monte Carlo transport codes.Since the void space needs to be described accurately,the void space modeling is more complicated than the solids modeling.For system with complex geometry,the void filling may generate too complicated Monte Carlo description for efficient modeling and calculation.Based on the multi-physics coupling analysis modeling program SuperMC/MCAM,an improved void filling algorithm based on convex solids was developed.This algorithm would test all the convex solids generated by decomposing the CAD solids with sub-spaces and generate the void description by complementary describing convex volumes intersected with the sub-spaces.This algorithm was tested with ITER reference neutronics model.The result shows that this algorithm simplifies the description of the void to 1/60,accelerates the void space modeling by 25 times and speeds up calculation by 1.4times compa-ring to benchmark algorithm.
引文
[1]TSIGE-TAMIRAT H,FISCHER U,SERIKOV A,et al.Automatic generation and validation of an ITER neutronics model from CAD data[J].Fusion Engineering and Design,2007,82:1 956-1 959.
[2]SATO S,IIDA H,NISHITANI T.Development of CAD/MCNP interface program prototype for fusion reactor nuclear analysis[J].Fusion Engineering and Design,2006,81:2 767-2 772.
[3]WU Y.CAD-based interface programs for fusion neutron transport simulation[J].Fusion Engineering and Design,2009,84(7-11):1 987-1 992.
[4]吴宜灿,李莹,卢磊,等.蒙特卡罗粒子输运计算自动建模程序系统的研究与发展[J].核科学与工程,2006,26(1):20-27.WU Yican,LI Ying,LU Lei,et al.Research and development of the automatic modeling system for Monte Carlo particle transport simulation[J].Chinese Journal of Nuclear Science and Engineering,2006,26(1):20-27(in Chinese).
[5]WU Y C,SONG J,ZHENG H Q,et al.CADbased Monte Carlo program for integrated simulation of nuclear system SuperMC[J].Annals of Nuclear Energy,doi:10.1016/j.anucene.2014.08.058.
[6]SONG Jing,SUN Guangyao,CHEN Zhenping,et al.Benchmarking of CAD-based SuperMC with ITER benchmark model[J].Fusion Engineering and Design,2014,89(11):2 499-2 503.
[7]吴宜灿,胡丽琴,龙鹏程,等.先进核能软件发展与核信息学实践[M]∥中国科研信息化蓝皮书.北京:科学出版社,2013:232-244.
[8]HU H,WU Y,CHEN M,et al.Benchmarking of SNAM with the ITER 3D model[J].Fusion Engineering and Design,2007,82:2 867-2 871.
[9]BOTH J P,MAZZOLO A,PENELIAU Y,et al.User manual for version 4.3of the TRIPOLI-4 Monte-Carlo method particle transport computer code[M].France:CEA,2003.
[10]LI Y,LU L,DING A,et al.Benchmarking of MCAM 4.0with the ITER 3Dmodel[J].Fusion Engineering and Design,2007,82(15-24):2 861-2 866.
[11]吴宜灿,李静惊,李莹,等.大型集成多功能中子学计算与分析系统VisualBUS的研究与发展[J].核科学与工程,2007,27(4):365-373.WU Yican,LI Jingjing,LI Ying,et al.An integrated multi-functional neutronics calculation and analysis code system:VisualBUS[J].Chinese Journal of Nuclear Science and Engineering,2007,27(4):365-373(in Chinese).
[12]WU Y,FDS Team.Conceptual design activities of FDS series fusion power plants in China[J].Fusion Engineering and Design,2006,81(23-24):2 713-2 718.
[13]WU Y,QIAN J,YU J.The fusion-driven hybrid system and its material selection[J].Journal of Nuclear Materials,2002,307-311:1 629-1 636.
[14]WU Y,FDS Team.Design analysis of the China Dual-Functional Lithium Lead(DFLL)test blanket module in ITER[J].Fusion Engineering and Design,2007,82:1 893-1 903.
[15]CHEN Y,WU Y.Conceptual study on high performance blanket in a spherical Tokamak fusiondriven transmitter[J].Fusion Engineering and Design,2000,49-50:507-512.
[16]吴宜灿,柏云清,宋勇,等.中国铅基研究反应堆概念设计研究[J].核科学与工程,2014,34(2):56-63.WU Yican,BAI Yunqing,SONG Yong,et al.Conceptual design of China Lead-based Research Reactor CLEAR-Ⅰ[J].Nuclear Science and Engineering,2014,34(2):56-63(in Chinese).
[17]王国忠,党同强,熊健,等.MCAM4.8在ITER建筑大厅中子学建模中的应用[J].核科学与工程,2011,31(4):351-355.WANG Guozhong,DANG Tongqiang,XIONG Jian,et al.Application of MCAM4.8in creating neutronics model for ITER building[J].Chinese Journal of Nuclear Science and Engineering,2011,31(4):351-355(in Chinese).
[18]BRIESMEISTER J F.MCNP:A general Monte Carlo N-particle transport code,version 4C,LA-13709-M[M].US:LANL,2000.
[19]李莹.蒙特卡罗粒子输运计算自动建模系统预处理关键技术研究[D].合肥:中国科学院合肥物质科学研究院,2009.
[20]LU L,LEE Y K,ZHANG J,et al.Development of Monte Carlo automatic modeling functions of MCAM for TRIPOLI-ITER application[J].Nuclear Instruments and Methods in Physics Research Section A,2009,605(3):384-387.
[2]SATO S,IIDA H,NISHITANI T.Development of CAD/MCNP interface program prototype for fusion reactor nuclear analysis[J].Fusion Engineering and Design,2006,81:2 767-2 772.
[3]WU Y.CAD-based interface programs for fusion neutron transport simulation[J].Fusion Engineering and Design,2009,84(7-11):1 987-1 992.
[4]吴宜灿,李莹,卢磊,等.蒙特卡罗粒子输运计算自动建模程序系统的研究与发展[J].核科学与工程,2006,26(1):20-27.WU Yican,LI Ying,LU Lei,et al.Research and development of the automatic modeling system for Monte Carlo particle transport simulation[J].Chinese Journal of Nuclear Science and Engineering,2006,26(1):20-27(in Chinese).
[5]WU Y C,SONG J,ZHENG H Q,et al.CADbased Monte Carlo program for integrated simulation of nuclear system SuperMC[J].Annals of Nuclear Energy,doi:10.1016/j.anucene.2014.08.058.
[6]SONG Jing,SUN Guangyao,CHEN Zhenping,et al.Benchmarking of CAD-based SuperMC with ITER benchmark model[J].Fusion Engineering and Design,2014,89(11):2 499-2 503.
[7]吴宜灿,胡丽琴,龙鹏程,等.先进核能软件发展与核信息学实践[M]∥中国科研信息化蓝皮书.北京:科学出版社,2013:232-244.
[8]HU H,WU Y,CHEN M,et al.Benchmarking of SNAM with the ITER 3D model[J].Fusion Engineering and Design,2007,82:2 867-2 871.
[9]BOTH J P,MAZZOLO A,PENELIAU Y,et al.User manual for version 4.3of the TRIPOLI-4 Monte-Carlo method particle transport computer code[M].France:CEA,2003.
[10]LI Y,LU L,DING A,et al.Benchmarking of MCAM 4.0with the ITER 3Dmodel[J].Fusion Engineering and Design,2007,82(15-24):2 861-2 866.
[11]吴宜灿,李静惊,李莹,等.大型集成多功能中子学计算与分析系统VisualBUS的研究与发展[J].核科学与工程,2007,27(4):365-373.WU Yican,LI Jingjing,LI Ying,et al.An integrated multi-functional neutronics calculation and analysis code system:VisualBUS[J].Chinese Journal of Nuclear Science and Engineering,2007,27(4):365-373(in Chinese).
[12]WU Y,FDS Team.Conceptual design activities of FDS series fusion power plants in China[J].Fusion Engineering and Design,2006,81(23-24):2 713-2 718.
[13]WU Y,QIAN J,YU J.The fusion-driven hybrid system and its material selection[J].Journal of Nuclear Materials,2002,307-311:1 629-1 636.
[14]WU Y,FDS Team.Design analysis of the China Dual-Functional Lithium Lead(DFLL)test blanket module in ITER[J].Fusion Engineering and Design,2007,82:1 893-1 903.
[15]CHEN Y,WU Y.Conceptual study on high performance blanket in a spherical Tokamak fusiondriven transmitter[J].Fusion Engineering and Design,2000,49-50:507-512.
[16]吴宜灿,柏云清,宋勇,等.中国铅基研究反应堆概念设计研究[J].核科学与工程,2014,34(2):56-63.WU Yican,BAI Yunqing,SONG Yong,et al.Conceptual design of China Lead-based Research Reactor CLEAR-Ⅰ[J].Nuclear Science and Engineering,2014,34(2):56-63(in Chinese).
[17]王国忠,党同强,熊健,等.MCAM4.8在ITER建筑大厅中子学建模中的应用[J].核科学与工程,2011,31(4):351-355.WANG Guozhong,DANG Tongqiang,XIONG Jian,et al.Application of MCAM4.8in creating neutronics model for ITER building[J].Chinese Journal of Nuclear Science and Engineering,2011,31(4):351-355(in Chinese).
[18]BRIESMEISTER J F.MCNP:A general Monte Carlo N-particle transport code,version 4C,LA-13709-M[M].US:LANL,2000.
[19]李莹.蒙特卡罗粒子输运计算自动建模系统预处理关键技术研究[D].合肥:中国科学院合肥物质科学研究院,2009.
[20]LU L,LEE Y K,ZHANG J,et al.Development of Monte Carlo automatic modeling functions of MCAM for TRIPOLI-ITER application[J].Nuclear Instruments and Methods in Physics Research Section A,2009,605(3):384-387.