考虑到粒子跟踪效率,在基于面描述的蒙特卡罗粒子输运模拟程序的计算模型中,除描述系统各部件的几何实体外,还需对非实体空间即空腔进行几何描述。由于空腔几何需严格填满实体间空隙,其建模复杂度通常较实体几何高,所以复杂系统常因空腔复杂度过高而严重制约蒙特卡罗建模和计算效率。本文研究了一种基于凸实体的自动空腔简化算法,该算法首先将栅元实体分解为凸实体集合,利用与子空间相交的凸实体的补集描述该子空间区域的空腔。该算法已在多物理耦合分析自动建模软件Super-MC/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.4 times comparing to benchmark algorithm.