CN104573279A - 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 - Google Patents
基于深度剥离的核辐射屏蔽计算的网格快速生成方法 Download PDFInfo
- Publication number
- CN104573279A CN104573279A CN201510044611.0A CN201510044611A CN104573279A CN 104573279 A CN104573279 A CN 104573279A CN 201510044611 A CN201510044611 A CN 201510044611A CN 104573279 A CN104573279 A CN 104573279A
- Authority
- CN
- China
- Prior art keywords
- solid
- grid
- geo
- radiation shield
- nuclear radiation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000005855 radiation Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 title abstract description 20
- 238000005516 engineering process Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 239000007787 solid Substances 0.000 claims description 44
- 239000002245 particle Substances 0.000 claims description 13
- 230000032258 transport Effects 0.000 description 29
- 238000011960 computer-aided design Methods 0.000 description 20
- 238000011160 research Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 230000018109 developmental process Effects 0.000 description 6
- 238000012800 visualization Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 241001581492 Attila Species 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000005025 nuclear technology Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- PJOHVEQSYPOERL-SHEAVXILSA-N (e)-n-[(4r,4as,7ar,12br)-3-(cyclopropylmethyl)-9-hydroxy-7-oxo-2,4,5,6,7a,13-hexahydro-1h-4,12-methanobenzofuro[3,2-e]isoquinoline-4a-yl]-3-(4-methylphenyl)prop-2-enamide Chemical compound C1=CC(C)=CC=C1\C=C\C(=O)N[C@]1(CCC(=O)[C@@H]2O3)[C@H]4CC5=CC=C(O)C3=C5[C@]12CCN4CC1CC1 PJOHVEQSYPOERL-SHEAVXILSA-N 0.000 description 1
- 102100023126 Cell surface glycoprotein MUC18 Human genes 0.000 description 1
- 101000623903 Homo sapiens Cell surface glycoprotein MUC18 Proteins 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 239000004614 Process Aid Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- RFKMCNOHBTXSMU-UHFFFAOYSA-N methoxyflurane Chemical compound COC(F)(F)C(Cl)Cl RFKMCNOHBTXSMU-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000010458 rotten stone Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T5/00—Recording of movements or tracks of particles; Processing or analysis of such tracks
- G01T5/02—Processing of tracks; Analysis of tracks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Image Analysis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Architecture (AREA)
- Software Systems (AREA)
Abstract
本发明公开了一种基于深度剥离的核辐射屏蔽计算网格快速生成方法,针对复杂几何体难以快速生成核辐射屏蔽计算网格的问题提出解决方法,该方法特征在于,采用深度剥离技术生成几何体的轮廓像素矩阵,并在图像空间中变换得到几何体轮廓网格,进而基于轮廓网格根据扫描线方法得到几何体内部网格,最终实现核辐射屏蔽计算网格的快速生成。
Description
技术领域
本发明涉及核科学与技术领域中的核辐射屏蔽数值模拟计算方向,涉及计算机辅助设计(CAD)生成的几何模型到核辐射屏蔽计算网格的转换生成方法。
背景技术
中子物理学(又称中子学)是近代物理学的一个重要分支,它主要研究中子与各种物质间的相互作用。在核反应堆物理研究、工程设计、安全评估、燃料管理优化、环境监测与辐射屏蔽防护等各个领域,中子学分析都是其必不可少的重要分析方法。中子学分析主要包括中子输运、活化、燃耗、剂量分析等内容。其中,中子输运分析是重要基础和核心内容。
中子输运(Neutron Transport)是指中子在介质内运动的过程。中子输运理论是研究中子在介质内运动的过程和规律的基础理论。随着核能利用的蓬勃发展,它在核科学技术领域中已成为一个独立的基础理论科学,并在核反应堆物理、屏蔽和核技术的工程应用和军事技术等领域中获得广泛的应用。中子输运分析能够获得中子密度、中子通量密度、中子流密度等物理量结果。
研究中子输运问题的方法一般分为二类:一类称为确定论方法;另一类称为蒙特卡罗方法。在确定论方法中,根据问题的物理性质所建立起的数学模型可以用一个或一组确定的数学物理方程来表示,然后,对这些方程可以采用数学方法求出它的精确或近似的解。对于一般的实际问题,都采用数值方法来求解。
离散纵标法就是求解中子输运方程最常用的确定论方法之一。近年来,随着计算机技术的迅速发展以及研究工作的不断深入,离散纵标法已经成为研究中子输运问题的有效数值方法之一,特别是在屏蔽计算中获得了广泛的应用。
近年来,特别是近40年来,各国研究单位已先后发展了一系列辅助领域专家进行中子学分析工作的离散纵标法程序系统,在计算精度和速度上都达到了相当完善的程度,并发展成为国际通用程序,在工程上获得了广泛的应用。如美国橡树岭国家实验室(Oak Ridge National Laboratory,ORNL)研发的DOORS系统(包含一维ANISN、二维DORT和三维TORT等独立输运程序系统),美国洛斯阿拉莫斯国家实验室(Los AlamosNational Laboratory,LANL)研发的DANTSYS系统(包含一维ONEDANT、二维TWODANT、三维THREEDANT等独立输运程序系统)和PARTISN系统,日本大阪大学(Osaka University)研制的NSHEX系统,美国Transpire公司的旗舰产品Attila,美国HSW Technologies公司商业软件PENTRAN系统等。
上述软件系统中,除Attila、PENTRAN是由商业公司研发运作外,其余系统均为研究机构或大学研制。当前,广泛流行于各国研究机构与研究所的软件系统还主要是研究机构研制的系统,比如DOORS系统,PARTISN系统。这些系统研发历史悠久,且采用的算法经典、可靠,是其他系统有效性验证的参考对象。它们也是科研人员基于中子输运进行活化、剂量计算分析的重要基础工具。
然而,这些经典的离散纵标法中子输运分析系统,由于更多关注计算方法研究,在程序易用性等方面考虑不足,且开发起始年代较早,未能很好吸收后续计算机图形学、科学计算可视化等学科发展的成果,以至它们目前对大型复杂例题的计算建模和后处理的支持方面表现得较为不足,使得计算前后处理环节突出为输运分析工作的瓶颈。这使得在计算机计算能力显著提高的今天,不能很好地满足用户对计算例题复杂度、精度不断提高的需求,已严重限制了上述经典离散纵标法中子输运分析系统的扩展应用。
上述问题主要表现在于,对于当前常见的工程实际计算例题的规模、复杂度与精度要求,一方面,常规的手工文本方式描述已经变得不可实现,另一方面,基于日益复杂的计算几何产生的结果数据的空间分布日益复杂,常规的基于区域平均值、轴向分布的分析方式将会丢失大量数据细节,已不再满足要求。
随着计算机学科的不断发展,计算机辅助设计技术(Computer-Aided Design)、科学计算可视化技术日益得到国内外广大科研人员的深入研究,并已逐步发展成熟,成为计算机学科中的重要研究方向。CAD技术关注的是几何的表示、生成、操作与显示;而科学计算可视化技术则是为了帮助分析人员直观高效地分析日益庞大复杂的计算结果,并直观展现数据分布规律以及细节信息的一门技术。
因此,如何针对离散纵标法中子输运分析过程的具体需求,结合CAD、科学计算可视化等相关技术,发展易学易用、可视直观、自动化程度较高的高效可视化系统,简化前后处理流程,降低用户要求与工作量,整体改善分析环境、提升工作效率,就成为离散纵标法中子输运计算分析工具研究领域的重要课题。
辐射输运(或称粒子输运)是指辐射粒子(如中子、光子、电子、质子等)在介质内运动的过程。辐射输运理论是研究粒子在介质内运动的过程和规律的基础理论。辐射输运模拟计算能够获得粒子积分流量、通量密度、核热沉积等物理量结果。随着核能利用的蓬勃发展,它在核科学与技术领域中已成为一个独立的基础理论科学,并在核反应堆物理、辐射屏蔽与防护、以及核技术的工程应用等领域中获得广泛的应用。
辐射输运计算程序(如MCNP、TRIPOLI等)已广泛应用于反应堆物理、辐射屏蔽与防护、核探测、发射剂量学等领域。常规的手工文本建立蒙特卡罗计算模型的方法耗时耗力且易出错,而且难以实现对复杂系统几何的建模。
在核领域的中子学分析和计算中,采用离散纵标法的粒子输运计算程序得到了越来越广泛的应用。该方法主要优点是计算效率高,可以准确计算具有各向异性、中子注量率分布不均匀的反应堆,并可获得分布量的计算结果。
然而,复杂几何的数以百万的三维计算网格难以有效生成。发展的算法虽然实现了网格生成的自动化,却仍然存在耗时严重的问题。
针对复杂几何体难以快速生成离散纵标法粒子输运计算网格的问题,本方法采用深度剥离技术生成几何体的轮廓像素矩阵,并在图像空间中变换得到几何体轮廓网格,进而基于轮廓网格根据扫描线方法得到几何体内部网格,最终实现离散纵标法粒子输运计算网格的快速生成。
离散纵标法粒子输运计算网格快速生成方法与工具主要有欧洲核能机构核数据中心发展的BOT3P以及中科院FDS团队发展的SNAM。
BOT3P(BOlonga Transport analysis Pre-Post-Processors)是欧洲核能机构核数据中心针对DOORS、DANTSYS、PARTISN系统研发的计算前后处理辅助工具。其几何定义仍然需要用户以手动文本方式描述,其主要功能定位于简化网格划分工作,针对一份原始几何,可以方便地对其生成多种不同的网格划分。
中国科学院等离子体物理研究所胡海敏在充分调研和分析SN方法粒子输运计算程序自动建模方法的基础上,对建模过程中的模型文件格式识别、属性编辑、空腔处理及自动划分离散网格等关键技术问题进行了研究,并提出了合理可行的解决方法,发展了SNAM程序。
中国科学院等离子体物理研究所龙鹏程,提出了结合模型规则化与网格化过程的直接网格化建模思想,避免了规则化精度、网格化精度相互干扰的混乱局面。在建模加速算法方面,引入基于包围盒的分块策略和中子学合法几何识别机制,并结合轴向面自动切割、大包围盒体积非法实体自动分解等机制,发展了基于包围盒和辅助面切割的逐几何判断算法,将来自于工程结构设计的复杂CAD几何的建模时间从数十小时降低到一小时量级,解决了自动建模耗时严重的问题。
SNAM(SN Auto-Modeling System)是FDS团队针对离散纵标法粒子输运计算程序研发的具有自主知识产权的自动建模系统。它采用同MCAM类似的基于CAD技术的CAD几何到中子输运计算模型的自动转换建模思想,支持现有CAD几何的直接利用,将几何建模工作转移到功能强大、丰富、灵活的常规CAD环境中,从而改变传统建模方式,改善建模环境,为保证建模精度,提升建模效率提供保障。目前,SNAM设计思想的原理性验证工作已经完成。本论文的快速自动化可视建模工作就是在此基础上展开的,期望通过显著缩短建模时间来实现该思路的实用化。
仅SNAM实现了基于复杂几何体CAD模型自动生成计算网格,但网格生成耗时仍然较大。
深度剥离(Depth Peeling)是一种对深度值进行排序的技术。它的原理也比较直观,标准的深度检测使场景中Z值最小的点输出到屏幕上,也就是离我们最近的顶点。但还有离我们第二近的顶点,第三近的顶点存在。要想显示它们,可以用多遍渲染的方法。depthpeeling基于z-buffer的多层渲染,每一层渲染是基于上一层渲染的深度值基础上进行的。
深度剥离技术(Depth Peeling)的基本思想是对场景进行遍渲染,每遍渲染都比上一遍更深入场景一层,从而可以依次获取每个像素的最近、次近、第三近的、……、第N近的片段,最后使用混合技术合成所有片段获取像素的最后颜色,从而得到最终图像。深度剥离就是为了处理物体的透明渲染而发明的技术,该技术不需要三角面片排序就可以实现正确的透明渲染效果。该技术在Nvidia Geforce3上实现后,基于深度剥离技术的透明渲染已越来越受重视。
深度剥离通过多次绘制几何体以积累最终结果来进行工作。绘制的每次迭代剥离掉通过每个像素可见的单个面深度层。该核心算法相对于画家算法的优势在于:其无需对几何基元预先排序。
随着个人计算机系统的能力和复杂性的提高,越来越多地使用称为图形处理单元(GPU)的专用图形学绘制设备来执行图形学操作,诸如Z缓冲、深度剥离、光栅化等所需的那些图形学操作。GPU通常在经由标准总线(例如,AGP或者PCI Express)耦合至主板的可移动图形卡中使用。尽管GPU的能力有所不同,其通常可以用来执行诸如以下任务:二维(2D)图形数据的绘制、三维(3D)图形数据的绘制、图形用户界面(GUI)显示元素的加速绘制。GPU可以实现一个或多个应用程序接口(API),该API允许程序员调用GPU的功能。
GPU可以包括用于将数字图像绘制到成像设备的各种内置及可配置结构。
现有技术具有网格生成效率低的问题。现有技术主要通过逐网格判断法来实现基于CAD几何模型的计算网格生成,单个网格的生成依赖点与CAD几何的相互关系判断,开销大;网格生成总体时间开销依赖网格数量,当网格数越多时,时间开销特别大。在核反应堆全空间屏蔽计算时,一般网格数需要细分到数百上千的网格,计算网格的生成需要数个小时的时间。
发明内容
本发明技术克服了现有技术时间开销大的不足,提出一种基于深度剥离的核辐射屏蔽计算网格快速生成系统与方法,该方法将三维网格生成问题降低维度为二维问题,将耗时的三维图形空间中的几何关系判断问题改为二维图像空间中的像素处理问题,具有网格生成速度快的优点,可显著降低网格生成时间开销。
本发明技术解决方案:一种基于深度剥离的核辐射屏蔽计算网格快速生成系统与方法,它包含如下实现步骤:
(1)读入计算机辅助结构设计软件生成的CAD(计算机辅助设计)几何模型,解析出模型中的各几何体,统计出几何体个数N,并针对第一个几何体Geo-i(i=1)进行下一步操作;
(2)计算出几何体Geo-i的轴向包围盒Box-i;
(3)采用深度剥离技术从轴向包围盒Box-i轴向正方向往下生成几何体Geo-i的轮廓像素矩阵GeoFragPos-i;
(4)采用深度剥离技术从轴向包围盒Box-i轴向负方向往下生成几何体Geo-i的轮廓像素矩阵GeoFragNeg-i;
(5)在图像空间中,根据用户输入的网格化划分Mesh3D,结合轮廓像素矩阵GeoFragPos-i确定几何体Geo-i轴向上方向的轮廓边界网格BoundaryMeshPos-i;结合轮廓像素矩阵GeoFragNeg-i确定几何体Geo-i轴向下方向的轮廓边界网格BoundaryMeshNeg-i;
(6)基于轮廓边界网格BoundaryMeshPos-i与BoundaryMeshNeg-i,根据扫描线方法得到几何体内部网格,最终实现几何体Geo-i的离散纵标法粒子输运计算网格的生成;
(7)转到下一个几何体,即i值递增1;
(8)如i>N,则结束;否则,转到第(2)步。
本发明与现有技术相比的优点在于:现有技术具有网格生成效率低的问题。现有技术主要通过逐网格判断法来实现基于CAD几何模型的计算网格生成,单个网格的生成依赖点与CAD几何的相互关系判断,开销大;网格生成总体时间开销依赖网格数量,当网格数越多时,时间开销特别大。在核反应堆全空间屏蔽计算时,一般网格数需要细分到数百上千的网格,计算网格的生成需要数个小时的时间。而本发明的优点在于:
(1)本发明将耗时的三维图形空间中的几何关系判断问题改为二维图像空间中的像素处理问题,从而处理难度与时间代价显著下降。
(2)本发明采用图形计算处理器GPU提取三维几何模型的边界面,并转换为二维像素图像,充分利用了GPU硬件技术与高效处理性能,使得三维几何模型的边界面信息能够得以快速获得。
(3)本发明三维几何体内部空间对应的网格,无需直接判断,可通过几何体边界对应的网格扫描生成,避免了大量的耗时的几何体关系判断,较少了计算量,使得几何体内部空间与计算网格的对应关系能够快速生成。
总之,本发明最终实现核辐射屏蔽计算网格的快速生成。
附图说明
图1为本发明流程图。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示,本发明中的方法,实现过程如下:
(1)读入计算机辅助结构设计软件生成的CAD(计算机辅助设计)几何模型,解析出模型中的各几何体,统计出几何体个数N,并针对第一个几何体Geo-i(i=1)进行下一步操作;
(2)计算出几何体Geo-i的轴向包围盒Box-i(xmin,xmax;ymin,ymax;zmin,zmax);
(3)采用深度剥离技术从轴向包围盒Box-i轴向正方向往下生成几何体Geo-i的二维轮廓像素矩阵GeoFragPos-i(x,y),矩阵记录的是几何体轮廓像素在三维世界坐标系中的Z轴向坐标值zPos3D;
(4)采用深度剥离技术从轴向包围盒Box-i轴向负方向往下生成几何体Geo-i的二维轮廓像素矩阵GeoFragNeg-i(x,y),矩阵记录的是几何体轮廓像素在三维世界坐标系中的Z轴向坐标值zNeg3D;
(5)在图像空间中,根据用户输入的网格化划分Mesh3D,结合轮廓像素矩阵GeoFragPos-i(x,y)的Z轴向坐标值zPos3D,确定几何体Geo-i轴向上方向的轮廓边界网格BoundaryMeshPos-i(x,y)的Z轴向网格索引值ZPos;结合轮廓像素矩阵GeoFragNeg-i(x,y)的Z轴向坐标值zNeg3D,确定几何体Geo-i轴向下方向的轮廓边界网格BoundaryMeshNeg-i(x,y)的Z轴向网格索引值ZNeg;
(6)基于轮廓边界网格BoundaryMeshPos-i与BoundaryMeshNeg-i,根据扫描线方法得到几何体内部网格,最终实现几何体Geo-i的离散纵标法粒子输运计算网格的生成;
(a)对于几何体Geo-i轴向包围盒Box-i对应的计算网格标记三维矩阵MeshFlag-i(x,y,z)(x=1,……,X;y=1,……,Y,z=1,……,Z),初始化标记全部为0;
(b)赋值x=1;
(c)赋值y=1;
(d)查找BoundaryMeshPos-i(x,y)与BoundaryMeshNeg-i(x,y)网格标记的z值ZPos、ZNeg;如果ZPos<ZNeg,则转到下一步;否则,将MeshFlag-i(x,y,z)(z=ZNeg,……,ZPos)的标记值修改为1;
(e)赋值y=y+1;如果y>Y,则转到下一步;否则,转到(d);
(f)赋值x=x+1;如果x>X,则转到下一步;否则,转到(c);
(g)网格标记集合MeshFlag-i(x,y,z)中标记为1的网格即为几何体Geo-i的计算网格。
(7)转到下一个几何体,即i值递增1;
(8)如i>N,结束;否则,转到第(2)步。
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (2)
1.基于深度剥离的核辐射屏蔽计算的网格快速生成方法,特征在于实现步骤如下:
(1)读入计算机辅助结构设计软件生成的CAD几何模型,解析出几何模型中的各几何体,统计出几何体个数N,并针对第一个几何体Geo-i,i=1进行下一步操作;
(2)计算出几何体Geo-i的轴向包围盒Box-i;
(3)采用深度剥离技术从轴向包围盒Box-i轴向正方向往下生成几何体Geo-i的轮廓像素矩阵GeoFragPos-i;
(4)采用深度剥离技术从轴向包围盒Box-i轴向负方向往下生成几何体Geo-i的轮廓像素矩阵GeoFragNeg-i;
(5)在图像空间中,根据用户输入的网格化划分Mesh3D,结合轮廓像素矩阵GeoFragPos-i确定几何体Geo-i轴向上方向的边界网格BoundaryMeshPos-i;结合轮廓像素矩阵GeoFragNeg-i确定几何体Geo-i轴向下方向的边界网格BoundaryMeshNeg-i;
(6)基于轮廓边界网格BoundaryMeshPos-i与BoundaryMeshNeg-i,根据扫描线方法得到几何体内部网格,最终实现几何体Geo-i的离散纵标法粒子输运计算网格的生成;
(7)转到下一个几何体,即i值递增1;
(8)如i>N,则结束;否则,转到第(2)步。
2.根据权利要求1所述的基于深度剥离的核辐射屏蔽计算的网格快速生成方法,特征在于:所述方法适用于笛卡尔坐标系下的凸几何体。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510044611.0A CN104573279B (zh) | 2015-01-28 | 2015-01-28 | 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 |
PCT/CN2015/098019 WO2016119537A1 (zh) | 2015-01-28 | 2015-12-21 | 基于深度剥离的核辐射屏蔽计算网格生成方法和系统 |
US15/543,154 US20180018407A1 (en) | 2015-01-28 | 2015-12-21 | Depth peeling based nuclear radiation shield calculation grid generation method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510044611.0A CN104573279B (zh) | 2015-01-28 | 2015-01-28 | 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104573279A true CN104573279A (zh) | 2015-04-29 |
CN104573279B CN104573279B (zh) | 2016-05-18 |
Family
ID=53089332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510044611.0A Expired - Fee Related CN104573279B (zh) | 2015-01-28 | 2015-01-28 | 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20180018407A1 (zh) |
CN (1) | CN104573279B (zh) |
WO (1) | WO2016119537A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016119537A1 (zh) * | 2015-01-28 | 2016-08-04 | 中国科学院合肥物质科学研究院 | 基于深度剥离的核辐射屏蔽计算网格生成方法和系统 |
CN106126929A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 基于离散纵标法处理大规模内真空粒子输运问题的方法 |
CN106951683A (zh) * | 2017-02-28 | 2017-07-14 | 西安交通大学 | 一种用于核电厂厂房屏蔽计算的高效并行扫描方法 |
CN107437270A (zh) * | 2017-08-09 | 2017-12-05 | 安徽中科超安科技有限公司 | 一种基于包围盒截面图像的三维几何截面图像生成方法 |
CN109190144A (zh) * | 2018-07-12 | 2019-01-11 | 哈尔滨工程大学 | 一种任意形状放射源辐射屏蔽计算仿真方法 |
CN113886961A (zh) * | 2021-09-30 | 2022-01-04 | 中国科学院国家空间科学中心 | 基于航天器三维屏蔽的辐射效应计算方法、装置及设备 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109101704B (zh) * | 2018-07-24 | 2019-09-24 | 南华大学 | 贫化铀用作船用核动力舱室防护材料的可行性评价方法 |
CN113536617B (zh) * | 2021-03-31 | 2024-04-02 | 长安大学 | 一种复杂结构的精细有限元模型快速生成方法 |
CN116502373B (zh) * | 2023-06-30 | 2023-10-20 | 陕西空天信息技术有限公司 | 流道计算网格生成方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106301A (zh) * | 2013-01-22 | 2013-05-15 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗方法与特征线方法耦合的计算辐射屏蔽的方法 |
CN104268945A (zh) * | 2014-10-09 | 2015-01-07 | 中国科学院合肥物质科学研究院 | 一种基于等位面的聚变反应堆中子输运计算中的样条面处理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573279B (zh) * | 2015-01-28 | 2016-05-18 | 中国科学院合肥物质科学研究院 | 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 |
-
2015
- 2015-01-28 CN CN201510044611.0A patent/CN104573279B/zh not_active Expired - Fee Related
- 2015-12-21 WO PCT/CN2015/098019 patent/WO2016119537A1/zh active Application Filing
- 2015-12-21 US US15/543,154 patent/US20180018407A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106301A (zh) * | 2013-01-22 | 2013-05-15 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗方法与特征线方法耦合的计算辐射屏蔽的方法 |
CN104268945A (zh) * | 2014-10-09 | 2015-01-07 | 中国科学院合肥物质科学研究院 | 一种基于等位面的聚变反应堆中子输运计算中的样条面处理方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016119537A1 (zh) * | 2015-01-28 | 2016-08-04 | 中国科学院合肥物质科学研究院 | 基于深度剥离的核辐射屏蔽计算网格生成方法和系统 |
CN106126929A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 基于离散纵标法处理大规模内真空粒子输运问题的方法 |
CN106126929B (zh) * | 2016-06-24 | 2018-10-19 | 西安交通大学 | 基于离散纵标法处理大规模内真空粒子输运问题的方法 |
CN106951683A (zh) * | 2017-02-28 | 2017-07-14 | 西安交通大学 | 一种用于核电厂厂房屏蔽计算的高效并行扫描方法 |
CN106951683B (zh) * | 2017-02-28 | 2020-04-10 | 西安交通大学 | 一种用于核电厂厂房屏蔽计算的高效并行扫描方法 |
CN107437270A (zh) * | 2017-08-09 | 2017-12-05 | 安徽中科超安科技有限公司 | 一种基于包围盒截面图像的三维几何截面图像生成方法 |
CN107437270B (zh) * | 2017-08-09 | 2018-09-25 | 安徽中科超安科技有限公司 | 一种基于包围盒截面图像的三维几何截面图像生成方法 |
CN109190144A (zh) * | 2018-07-12 | 2019-01-11 | 哈尔滨工程大学 | 一种任意形状放射源辐射屏蔽计算仿真方法 |
CN109190144B (zh) * | 2018-07-12 | 2022-12-13 | 哈尔滨工程大学 | 一种任意形状放射源辐射屏蔽计算仿真方法 |
CN113886961A (zh) * | 2021-09-30 | 2022-01-04 | 中国科学院国家空间科学中心 | 基于航天器三维屏蔽的辐射效应计算方法、装置及设备 |
CN113886961B (zh) * | 2021-09-30 | 2022-08-16 | 中国科学院国家空间科学中心 | 基于航天器三维屏蔽的辐射效应计算方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
WO2016119537A1 (zh) | 2016-08-04 |
CN104573279B (zh) | 2016-05-18 |
US20180018407A1 (en) | 2018-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104573279B (zh) | 基于深度剥离的核辐射屏蔽计算的网格快速生成方法 | |
Wu et al. | CAD-based Monte Carlo program for integrated simulation of nuclear system SuperMC | |
Bremer et al. | Interactive exploration and analysis of large-scale simulations using topology-based data segmentation | |
Wilson et al. | Acceleration techniques for the direct use of CAD-based geometry in fusion neutronics analysis | |
Leppänen | Response matrix method–based importance solver and variance reduction scheme in the serpent 2 Monte Carlo code | |
Yang et al. | MOCUM: A two-dimensional method of characteristics code based on constructive solid geometry and unstructured meshing for general geometries | |
Gan et al. | Accelerating solvers for global atmospheric equations through mixed-precision data flow engine | |
CN103065056B (zh) | 一种基于数据场分割的移动人体剂量蒙特卡罗模拟方法 | |
CN113886961B (zh) | 基于航天器三维屏蔽的辐射效应计算方法、装置及设备 | |
García et al. | Development of the automatic void generation module in GEOUNED conversion tool | |
Zhang et al. | Fast solution of neutron transport SP3 equation by reduced basis finite element method | |
CN106202865B (zh) | 一种计算中子输运离散节块法中任意阶系数的方法 | |
Wang et al. | A novel and efficient method for calculating beam shadows on exterior surfaces of buildings in dense urban contexts | |
James | A Hardware Track-Trigger for CMS: At the High Luminosity LHC | |
Äkäslompolo et al. | Serpent neutronics model of Wendelstein 7-X for 14.1 MeV neutrons | |
Cottin | Reconstructing particle masses in events with displaced vertices | |
CN115114825B (zh) | 一种核反应堆退役方案的数值仿真方法和系统 | |
CN103065358A (zh) | 一种基于影像体元运算的器官几何重建方法 | |
Ladshaw et al. | Algorithms and algebraic solutions of decay chain differential equations for stable and unstable nuclide fractionation | |
Li et al. | MIGSHIELD: A new model-based interactive point kernel gamma ray shielding package for virtual environment | |
Yang et al. | A fast gamma-ray dose rate assessment method for complex geometries based on stylized model reconstruction | |
Fischer et al. | Use of CAD generated geometry data in Monte Carlo transport calculations for ITER | |
Zavorka et al. | Hybrid unstructured mesh geometry model of the next-generation spallation neutron source at LANSCE | |
Wu et al. | Comprehensive Neutronics Simulations | |
Chen et al. | Volume data numerical integration and differentiation using CUDA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160518 Termination date: 20210128 |
|
CF01 | Termination of patent right due to non-payment of annual fee |