基于非结构网格粒子输运的蒙特卡罗模拟计算研究  

Study on the Particle Transport Calculation with Monte Carlo Method Based on Unstructured Grid

作　　者：王立鹏[1] 陈森 张信一[1] 姜夺玉 胡田亮 曹璐 李达[1] 陈立新[1] WANG Lipeng;CHEN Sen;ZHANG Xinyi;JIANG Duoyu;HU Tianliang;CAO Lu;LI Da;CHEN Lixin(Northwest Institute of Nuclear Technology,State Key Laboratory of Intense Pulsed Radiation Simulation and Effect,Xi an of Shaanxi Prov.710024,China)

机构地区：[1]西北核技术研究所强脉冲辐射环境模拟与效应国家重点实验室,陕西西安710024

出　　处：《核科学与工程》2024年第2期286-294,共9页Nuclear Science and Engineering

基　　金：自然科学基金的资助(12275219,12205237)

摘　　要：传统的构造实体几何的建模方式不能满足复杂几何输运、多物理耦合计算和大尺度核辐射效应模拟等问题,目前MCNP程序允许在非结构化的网格下跟踪粒子,它是作为一个网格空间嵌入在传统的构造实体几何之中,非结构网格几何采用有限元程序ABAQUS创建,使用MCNP的径迹长度估计器计算非结构网格单元内的输运结果,并在非结构网格上执行单元对单元的跟踪,结果可以输出到ABAQUS或者Paraview进行可视化和其他物理分析。本文分析了两种典型的基准实验:Ueki固定源和Godiva临界源问题,模拟结果表明,UM几何的运行时间和与CSG几何的计算偏差与输入网格单元的数量、单元的阶数成正比,同样阶数的四面体比六面体结果误差更小,本文的非结构网格计算结果与传统的构造实体几何结构结果具有很好的一致性。The traditional modeling method of constructing solid geometry cannot meet the problems of complex geometric transport,multi-physical coupling calculation,and large-scale nuclear radiation effect simulation.As a mesh space embedded in the traditional structural entity geometry,the unstructured grid geometry is created by using the finite element program ABAQUS,and the transport results in the unstructured mesh grid element are calculated by using the track length estimator of MCNP,and the element tracking is performed on the unstructured mesh grid.The results can be exported to ABAQUS or Paraview for visualization and other physical analysis.In this paper,we analyze two typical benchmark experiments:the Ueki fixed source and the Godiva critical source.The results show that the calculation deviation and the running time are proportional to the number of the input mesh grid cells and the order of the elements,and the results of the same order tetrahedron are smaller than that of hexahedron.The results of unstructural mesh grid calculation in this paper are in good agreement with those of traditional constructing solid geometry.

关 键 词：非结构网格 粒子输运 蒙特卡罗 有限元 

分 类 号：TL32[核科学技术—核技术及应用]