CN104361150B - 一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化系统 - Google Patents

Abstract

本发明公开了一种用于中子输运的基于MCAM‑Geant4自动建模方法的转换简化方法，根据基本体特征，根据查找和组合的思想，对从半空间几何到基本体几何的转换进行简化的方法。该转换方法简化了从半空间描述转换为基于基本体的蒙特卡罗计算模型，从而使Geant4蒙特卡罗计算程序在转换时提供简化。该方法对于核辐射评估等必须利用蒙特卡罗程序的领域提供了使用基于基本体描述的蒙特卡罗计算程序的简便方法，具有重要的科学研究与工程应用价值。

Description

一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化系统

技术领域

本发明涉及一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化系统，属于辐射输运计算领域，是一种从半空间描述转化为基于基本体描述的计算模型中的简化方法。

背景技术

蒙特卡罗输运程序Geant4是使用蒙特卡罗方法通过计算机模拟粒子运动进而模拟粒子输运过程的技术。该技术是在核辐射防护、反应堆物理、高能物理领域有广泛的应用。这种计算程序从每个粒子运行的角度对真实场景进行模拟，在计算的过程中需要知道粒子运行路径中与哪些物体发生了碰撞，碰撞后飞行进入下一个粒子的位置等等。为了能够快速判断粒子的位置，传统的蒙特卡罗输运程序使用了半空间方法描述所模拟场景的几何。但是随着计算几何技术的发展，为了直观地表示物体的几何，Geant4使用基于基本体的表达方式来表示几何。半空间几何和基本体几何的区别如图1.所示，二者主要在表达方式上有区别，如图1所示。

但是单纯使用文本对场景进行建模的效率是十分低下的，但现在已经有了成熟的从CAD模型建立基于半空间几何的蒙卡模型的软件MCAM-Geant4。这个软件大大提升了蒙特卡罗模拟的效率，但它直接将模型转换后的半空间描述方式的每一个半空间描述为一个基本体。这种方法比较简单，直接将半空间描述的模型转换为一个个的单独的基本体，但是增加了基本体的数量，造成了基本体的冗余，降低了计算速度，增加了计算时的内存开销等。

由于蒙特卡罗方法中，在判断粒子与几何体的相对位置时需要比较粒子和组成几何体的所有表面的位置，所以几何体的组成越简单，计算速度越快，所以在这种情况下，基本体的数量对几何描述的实用程度有密切的关系。本发明提出的基本体几何精简方法对实现CAD模型到蒙卡模型的有效转换具有重要意义。

发明内容

本发明技术解决问题：克服现有技术的不足，提出一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化系统，采用查找和组合思想，尽可能的查找半空间直到形成一个完整的基本体，或基于半空间尽可能的组合为一个基本体，已达到减少基本体，提高建模质量，减少相同精度下模型模拟时间。

本发明的技术方案：一种基于MCAM-Geant4自动建模方法的转换简化系统，由输入模块，查找模块、组合模块、输出模块共四个模块组成；输入模块输入半空间描述的模型，送至查找模块进行查找；查找模块根据输入的半空间描述的模型查找能否构成为一个完整的基本体，如果不能传递数据给组合模块；组合模块根据输入的半空间描述的模型查找能否使用尽可能多的半空间构成一个基本体，并在剩余的半空间中继续描述；最后查找模块构成的基本体、组合模块组成的基本体，统一送至输出模块，进行输出；

输入模块，读取MCAM-Geant4模块中读取半空间描述的模型数据，根据几何构造器的数据结构，然后送到组合模块中；

查找模块，根据读取到的半空间模型数据进行查找；

读取平面半空间，得到其法向量，平面上一点，方向等几何数据；读取圆柱半空间，得到其轴向法向量，轴线上一点，半径，内外半空间指示等几何数据；读取圆锥空间，得到其轴向法向量，顶点，圆锥与轴线夹角，内外半空间指示等几何数据；读取球半空间，得到球心点，半径，内外半空间等几何数据；读取圆环半空间，得到轴线向量，轴心点，轴向大半径，小半径，内外半空间等几何数据；

判断一个实体下面的所有半空间是否能够组合成一个基本体，对于实体，判断是否有八个半空间，且对于这八个半空间对应的八个面中存在两个面平行并垂直与其它面，存在其余六个面两两相邻平面成120°角度；同时半空间的正负方向符合六棱柱体的方向，平行的两个面对应的两个半空间交集不为空；据此构建六棱柱体；

判断是否有且仅有六个半空间，且六个半空间对应的六个面存在两个面平行并垂直与其余四个面，剩余的四个面两两垂直或两两平行；同时半空间的正负方向符合长方体的方向，平行的两个面对应的两个半空间交集不为空；据此构建长方体；

判断是否有且仅有三个半空间，且三个半空间对应的三个面有两个平面一个圆柱面；两个平面平行并垂直圆柱面；三个半空间的交集不为空；据此构建圆柱体；

判断是否有且仅有两个半空间，对应的两个面是平面和圆锥面；平面垂直圆锥面；两个半空间的交集不为空；据此构建圆锥体；

组合模块，根据读取到的半空间模型数据，根据优先级进行组合；主要优先级原则是优先组合为含较多半空间的基本体；输入数据为得到进行查找，但得不到基本体的半空间描述的模型；输出为基本体集合；

f)对于第一优先级六棱柱，遍历所有半空间a，进行判断：

i.首先搜索几何体中与之垂直或平行或成120°的平面，如果有这样的半空间b，判断找到的半空间与a包围成的空间是否包含原几何体，如果包含，那么以b作为将代替a的六棱柱的另一个面，当前的平面集合A＝{a,b}；

ii.继续在几何体中搜索其它半空间，新的半空间需要与集合A中最多一个半空间平行且与其他的半空间垂直，并且与这些半空间一起交出的空间能够包含原几何体，找到后将新的半空间加入集合A，最后集合A的半空间的基础上补齐的六棱柱可以成为在基本体表达式中代替这些半空间的基本体，从而使最终生成的基本体表达式减少了新找到半空间数量的基本体；

g)对于第二优先级长方体，搜索几何体中与之构成长方体的半空间，补齐缺少的面形成长方体，这样能够减少找到面数量的基本体；

i.首先搜索几何体中与之垂直或平行或成90°的平面，如果有这样的半空间b，判断找到的半空间与a包围成的空间是否包含原几何体，如果包含，那么以b作为将代替a的六棱柱的另一个面，当前的平面集合A＝{a,b}；

ii.继续在几何体中搜索其它半空间，新的半空间需要与集合A中最多一个半空间平行且与其他的半空间垂直，并且与这些半空间一起交出的空间能够包含原几何体，找到后将新的半空间加入集合A，最后集合A的半空间的基础上补齐的长方体可以成为在基本体表达式中代替这些半空间的基本体，从而使最终生成的基本体表达式减少了新找到半空间数量的基本体；

h)对于第三优先级圆柱面，搜索几何体中与之构成圆柱体的面，采用类似2.a的方法，最后补齐缺少的面形成圆柱体，这样能够减少找到面数量的基本体；

i.首先找到圆柱体内半空间；找到后再查找垂直于这个圆柱体的两个平面半空间；

ii.找到后，使用圆柱体代替找到的圆柱内半空间和平面半空间；

i)对于第四优先级圆锥面，搜索几何体中与之构成圆锥体的面，采用类似2.a的方法，最后补齐缺少的面形成圆锥体，这样能够减少找到面数量的基本体；

i.首先找到圆锥体内半空间；找到后再查找垂直于这个圆锥体的一个平面半空间；找到后，使用圆锥体代替找到的圆锥体内半空间和平面半空间；

本发明与现有技术相比的优点在于：

(1)通过直接寻找基本体，将一个实体描述为一个基本体，对比其他方式中采用的一个半空间使用一个基本体描述的方式，减少了基本体的数量；例如六棱柱，按照原来的方法，需要使用八个布尔实体才能描述，现在只需要一个六棱柱实体即可，实体数量减少了七倍；对比其他方式中采用的一个半空间使用一个基本体描述的方式，本发明减少了基本体的数量，减少了相同精度下模型模拟时间；

(2)通过查找基本体的基本组成元素，最大可能的组合成基本体，对比其他方式中采用的一个半空间使用一个基本体描述的方式，尽可能的减少了基本体数量；例如对于一个六面体，主要有三个半空间两两垂直，就可以使用一个基本体来代替这三个半空间，基本体从原来的六个减少到四个，效果提高了1.5倍。

附图说明

图1为基于半空间几何和基于基本体的几何的示意图；其中左图是基于半空间几何，右图是基于基本体的几何；

图2为本发明系统的组成框；

图3为图2中查找模块的实现流程图；

图4为图2中组合模块的实现流程图。

具体实施方式

如图2所示，本发明包括输入模块、查找模块、组合模块、输出模块等。给出每个模块的功能和作用及相互的连接关系。

输入模块读取已分解后的模型信息，以半空间表达形式存在；

在得到输入模块传递的输入后，本发明中查找模块具体实现如下，并如图3所示：

得到半空间描述的实体列表，对于其中每一个实体进行遍历；例如一个实体模型，含有六个半空间面；

查找实体列表中的每一个实体是否是一个基本体；根据权利要求1中列出的不同几何体的各种不同的参数，判断是否能够构成一个基本几何体，缺乏部分条件也可以构成一个基本体；以上面的模型为例，经过判断可以构成一个长方体，则输出为一个长方体；

在得到查找模块传递的输入后，本发明中组合模块具体实现如下，并如图4所示：

对于查找之后发现不能直接构成一个基本几何体的实体，采取第二种方式：组合，尽可能将更多的半空间组成成一个栅元，以一个四面体为例，其中三个平面半空间为两两垂直相交的平面半空间，最后一个是平面半空间：

a)基于优先级，先查找有无平面；再查找有无其它平面满足平行或垂直等能构成基本体的条件的平面：若有，加入其中，再继续检索下一个面；若无，判断其它面；添加一个平面半空间进行判断；

b)根据优先级来，判断完平面，再判断圆柱，圆锥；判断可以组成一个平面半空间；

c)尽最大可能判断之后，得到的半空间列表转换为一个基本体；再去判断其它平面半空间，最终发现两两垂直的三个平面半空间可以构成为一个长方体的三个面；

d)继续对剩余的面进行判断，重复abc步骤，并判断能否增加新的半空间到集合中；判断发现已经遍历所有的半空间，但不能再增加新的半空间到半空间集合中；

e)结束判断：

■如果所有的半空间遍历完，且所有的半空间都被处理，则结束；

■如果所有的半空间遍历完，但存在半空间未被处理，则保存目前的基本体到输出列表中，删去已经处理过的半空间，重新处理其它半空间；

■如果不是所有的半空间都被遍历完成，但已经构成完整的基本体，则保存该基本体到输出列表中，删去已经处理过的半空间，重新处理其它半空间；

■如果不是所有的半空间都被遍历完成，且基本体不完整，继续处理；

■根据以上判断，上述四面体例子在判断三个半空间构成一个长方体后，将其加入到输出列表中等待输出，并继续处理剩余的一个平面半空间；

在结束查找模块和组合模块的处理后，本发明中，输出模块输出所需要的Geant4模型文件；使用Geant4模拟代码进行模拟；由于降低Geant4模型复杂度，减少了Geant4模型中的实体数量，因此Geant4模拟速度提高。

Claims (1)

1.一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化系统，其特征在于包括输入模块、查找模块、组合模块和输出模块；其中：

输入模块，输入半空间描述的模型，送至查找模块进行查找；

查找模块，根据输入的半空间描述的模型查找能否构成为一个完整的基本体，如果能构成一个完整的基本体，输出到输出模块中；如果不能构成为一个完整的基本体，则传递数据给组合模块；所述基本体是指长方体、圆柱体、圆锥体、六棱柱体、球体、圆环体；

组合模块，根据查找模块传递来的半空间描述的模型，查找能否使用尽可能多的半空间组成一个基本体，并将结果输出到输出模块中；

输出模块，将查找模块构成的基本体、组合模块组成的基本体统一输出；

所述查找模块具体实现如下：

(1)根据读取到的半空间模型数据进行查找；

读取平面半空间，得到平面半空间的法向量、平面上一点和方向的几何数据；读取圆柱半空间，得到圆柱半空间的轴向法向量、轴线上一点、半径、内外半空间指示的几何数据；读取圆锥空间，得到圆锥空间的轴向法向量、顶点、圆锥与轴线夹角、内外半空间指示的几何数据；读取球半空间，得到球心点，半径，内外半空间的几何数据；读取圆环半空间，得到轴线向量，轴心点，轴向大半径，小半径，内外半空间的几何数据；

判断一个实体下面的所有半空间是否能够组合成一个基本体，对于实体，判断是否有八个半空间，且对于这八个半空间对应的八个面中存在两个面平行并垂直于其它面，存在其余六个面两两相邻平面成120°角度；同时半空间的正负方向符合六棱柱的方向，平行的两个面对应的两个半空间交集不为空；据此构建六棱柱；

判断是否有且仅有六个半空间，且六个半空间对应的六个面存在两个面平行并垂直于其余四个面，剩余的四个面两两垂直或两两平行；同时半空间的正负方向符合长方体的方向，平行的两个面对应的两个半空间交集不为空；据此构建长方体；

判断是否有且仅有三个半空间，且三个半空间对应的三个面有两个平面一个圆柱面；两个平面平行并垂直圆柱面；三个半空间的交集不为空；据此构建圆柱体；

判断是否有且仅有两个半空间，对应的两个面是平面和圆锥面；平面垂直圆锥面；两个半空间的交集不为空；据此构建圆锥体；

(2)对于不能组成一个基本体的实体，保存数据并传递给组合模块；

所述组合模块具体实现如下：

根据读取到的半空间模型数据，根据优先级进行组合；主要优先级原则是优先组合为含较多半空间的基本体；根据输入数据进行查找，但得不到基本体的半空间描述的模型；输出为基本体集合；

a)对于第一优先级六棱柱，遍历所有半空间a，进行判断：

i.首先搜索几何体中与之垂直或平行或成120°的平面，如果有这样的半空间b，判断找到的半空间与a包围成的空间是否包含原几何体，如果包含，那么以b作为将代替a的六棱柱的另一个面，当前的平面集合A＝{a,b}；

ii.继续在几何体中搜索其它半空间，新的半空间需要与集合A中最多一个半空间平行且与其他的半空间垂直，并且与这些半空间一起交出的空间能够包含原几何体，找到后将新的半空间加入集合A，最后集合A的半空间的基础上补齐的六棱柱可以成为在基本体表达式中代替这些半空间的基本体，从而使最终生成的基本体表达式减少了新找到半空间数量的基本体；

b)对于第二优先级长方体，搜索几何体中与之构成长方体的半空间，补齐缺少的面形成长方体，这样能够减少找到面数量的基本体；

i.首先搜索几何体中与之垂直或平行或成90°的平面，如果有这样的半空间b，判断找到的半空间与a包围成的空间是否包含原几何体，如果包含，那么以b作为将代替a的六棱柱的另一个面，当前的平面集合A＝{a,b}；

ii.继续在几何体中搜索其它半空间，新的半空间需要与集合A中最多一个半空间平行且与其他的半空间垂直，并且与这些半空间一起交出的空间能够包含原几何体，找到后将新的半空间加入集合A，最后集合A的半空间的基础上补齐的长方体可以成为在基本体表达式中代替这些半空间的基本体，从而使最终生成的基本体表达式减少了新找到半空间数量的基本体；

c)对于第三优先级圆柱面，搜索几何体中与之构成圆柱体的面，采用类似a)的方法，最后补齐缺少的面形成圆柱体，这样能够减少找到面数量的基本体；

i.首先找到圆柱体内半空间；找到后再查找垂直于这个圆柱体的两个平面半空间；找到后，使用圆柱体代替找到的圆柱体内半空间和平面半空间；

d)对于第四优先级圆锥面，搜索几何体中与之构成圆锥体的面，采用类似a)的方法，最后补齐缺少的面形成圆锥体，这样能够减少找到面数量的基本体；

i.首先找到圆锥体内半空间；找到后再查找垂直于这个圆锥体的一个平面半空间；找到后，使用圆锥体代替找到的圆锥体内半空间和平面半空间；

e)对于全部实体检索完，输出给输出模块进行输出。