一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统的制作方法

一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统的制作方法
【专利摘要】本发明是一种用在反应堆设计中，用以进行堆芯部件设计的建模系统，本系统通过分析裂变堆芯几何特殊性，包含大量重复构造的几何，通过构建参数自动创建裂变堆堆芯的CAD工程模型和蒙特卡罗计算模型，而蒙特卡罗计算模型是辐射输运蒙特卡罗计算程序的输入，得到此模型后就能计算出反应堆堆芯的关键物理量。通过不断对计算结果中物理参数的评价，此发明可以用来对裂变堆芯模型的设计方案进行反复迭代修改，直至形成用户满意堆芯方案。此发明意义在于避免人手工编写蒙特卡罗计算程序输入文件的繁琐性，使整个建模过程变得极其直观，极大降低出错的概率，加快堆芯的设计时间。
【专利说明】一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统

【技术领域】
[0001] 本发明涉及一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，应用在基于蒙特 卡罗计算的裂变堆堆芯设计中，目标是快速自动化构建裂变堆堆芯CAD工程模型和蒙特卡 罗计算模型。

【背景技术】
[0002] 随着反应堆设计和建造技术的进步，各种新型反应堆的设计频繁面世，蒙特卡罗 计算程序由于能适应任意复杂几何和计算结果精确等特点，开始在堆芯概念设计和工程 设计上得到广泛应用，而完成复杂几何的建模是编写计算程序输入文件中的最核心部分。 目前主要采用以下几种方法：第一，手工建模，耗时耗力，容易出错；第二，脚本程序自动建 模，无法生成工程CAD模型，不能进行迭代设计操作，脚本程序需要经常修改；第三，专业的 建模软件，例如中国科学院核能安全技术研究所开发的多物理耦合建模分析自动建模软件 MCAM。
[0003] 目前大多数专业建模软件中处理复杂模型的方法是对每个栅元进行拆分解析，如 果直接将此种方法应用在裂变堆芯模型上，在时间性能和计算机资源消耗上根本无法接 受，而且转换得到的模型也丢失了原有层级信息，在处理大规模裂变堆芯模型时显得力不 从心。
[0004] 由于裂变堆芯存在大量简单几何重复并互相填充的特点，蒙特卡罗计算程序一般 都会提供一种层级几何的描述方式来构建堆芯的几何，这种几何描述方式结构十分清晰， 可以描述几何之间互相嵌套填充的关系，大大缩减几何文件长度，加快几何导入速度，提高 计算效率。于是，人们开始针对裂变堆芯这种特点，针对性的开发新的建模软件。早先 MCAM 中开发的裂变堆芯快速建模模块，无法解析蒙特卡罗计算模型中的层级结构，也无法完成 堆芯几何迭代修改，无法应用在实际堆芯的设计中。随后北京应用物理与计算数学研究所， 美国阿贡国家实验室都开发了类似的专业建模系统，不过均没有证明这些系统能够支持堆 芯的反复迭代设计，也无法解析现有蒙特卡罗计算模型，并且由于堆芯中对模型可视化和 计算的精度是灵活要求的，现有的技术也无法提供可供用户选择的功能。


【发明内容】

[0005] 本发明技术解决问题：克服现有技术的不足，提供一种基于蒙特卡罗计算的反应 堆堆芯迭代设计系统，自动快速构建裂变反应堆堆芯CAD工程模型和蒙特卡罗计算模型， 并支持在整个设计过程中迭代设计，其中反应堆CAD工程模型可以用来满足后续工程设计 样本参照，而蒙特卡罗计算模型则可以输入辐射输运蒙特卡罗计算程序中进行计算，得到 反应堆堆芯功率、K有效，燃耗等十分重要物理量，这些物理量也是评价一个堆芯是否满足 设计要求和安全性的必要物理量。本发明避免人手工编写蒙特卡罗计算程序输入文件的繁 琐性，使整个建模过程变得极其直观，极大降低出错的概率，加快堆芯的设计时间。
[0006] 本发明技术解决方案：一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，此系 统实现了基于反应堆堆芯迭代设计的建模过程，其能够生成堆芯的CAD工程设计模型和堆 芯的蒙特卡罗计算模型，其中CAD工程设计模型可以用于肉眼观察堆芯结构和工程设计， 蒙特卡罗计算模型能够导入蒙特卡罗计算程序，计算出反应堆堆芯关键物理参数，达到堆 芯的设计检验作用；
[0007]本发明包括核心转换模块、参数控制模块、计算模型解析模块和可视化交互模 块；
[000S]核心转换模块，主要包含一套所有模块共用的满足反应堆堆芯设计参数体系，适 用于反应堆堆芯部件存在多层几何嵌套填充的特点，逻辑上为一种树形结构，称为参数树； 基于此参数树此模块主要完成三维CAD模型和蒙特卡罗计算程序输入模型的生成和相互 转换；所述参数控制模块与计算模型解析模块所提供参数全部输入到核心转换模块中，由 核心转换模块构建参数树，可视化交互模块从核心转换模块中提取参数树，基于参数树，可 视化交互模块对模型进行三维渲染，而参数控制模块能够继续对界面中渲染出的三维模型 进行修改操作，修改操作产生的新的参数会传递给核心转换模块，核心转换模块负责完成 对参数树的修改，在整个过程中核心转换模块随时能够将参数树中的信息转换成蒙特卡罗 计算程序输入文件，蒙特卡罗计算程序输入文件用于基于蒙特卡罗方法的输运计算，得到 堆芯的关键物理量，评价设计是否满足各方面要求。
[0009]参数控制模块，实现通过系统界面提供修改操作和通过系统界面提供对话框输入 新的参数来构建、修改和控制堆芯的模型，并将所有参数输入到核心转换模块；
[0010]计算模型解析模块，自动识别、检测和拆分用户已有的反应堆堆芯蒙特卡罗计算 模型，最终得到堆芯各部件具体参数，将参数输入到核心转换模块之中；
[0011] 可视化交互模块，渲染出通过核心转换模块构建好的计算机三维CAD模型，让用 户实时对模型的设计状态可见；
[0012] 所述核心转换模块具体实现过程为：
[0013] (1)获取参数控制模块和计算模型解析模块参数，如果获取参数来源于参数控制 模块判定参数树是否为空，如果非空，说明已经存在上一代参数树，则通过界面上存在限定 约束关系直接找到新的参数对应节点，在此节点中替换对应旧参数，完成对上一代参数树 的更新；如果为空则根据界面上存在的约束关系直接到步骤（5);如果获取参数来源于是 计算模型解析模块，需经历步骤（2)，（3)，（4);
[0014] ⑵构建参数子树，参数子树中包含蒙特卡罗计算模型中限定的一组栅元集合，选 此集合的标号为根节点，集合内部的所有栅元为叶子节点，集合内部的单个栅元成为叶子 栅元，遍历所有叶子栅元，如果解析得到叶子栅元仍然被另一个栅元集合给填充，则将填充 此叶子栅元所使用的栅元集合的标号作为此叶子栅元的子节点，从而完成单棵参数子树构 建；
[0015] (3)合并参数子树，通过各参数子树中根节点和最下层叶子节点信息，就能讲所有 子树链接起来，构建成一颗完整参数树；
[0016] (4)归并参数树，由于蒙特卡罗计算模型中栅元中几何有多个半空间面组成，几何 信息分散零碎，通过从底向上归并参数树节点；
[0017] (5)完成最终参数树构建此参数树，根据部件之间的包含关系划分出对应的层级 属性；
[0018] (6)遍历参数树中的信息，通过几何信息构造三维几何实体，输入到可视化模块进 行渲染；
[0019] (7)遍历参数树中的信息，通过参数树中划分的层级属性，按层级的描述堆芯的每 个栅元，包括几何描述和其他物理描述，最终输出保留了层级关系的蒙特卡罗计算模型，用 以输入蒙特卡罗计算程序进行计算。
[0020]所述参数控制模块具体实现过程为：
[0021] (1)通过扑捉用户在系统界面操作，获取用户输入参数和操作信息具体包含构建 堆芯各组件几何，进一步细化堆芯组件内部燃料棒阵列，改变组件和组件内部元件的阵列 规则，改变控制棒位置，改变控制棒插入状态，燃料棒轴向生长，燃料棒径向膨胀和改变各 部分几何对应物理属性；
[0022] (2)从核心转换模块中获取已有参数树，如果参数树为空，则按照用户输入参数和 操作信息自动构建参数树；
[0023] (3)如果参数树不为空，则将用户新输入的控制参数和操作产生的结果反馈进参 数树之中，并继续保存在核心转换模块。
[0024] 所述计算模型解析模块具体实现过程为：
[0025] (1)用户导入已有的蒙特卡罗计算模型；
[0026] (2)按照蒙特卡罗计算模型语法规则，进行解析，得到蒙特卡罗计算模型中每一个 栅元的参数，栅元中包含几何信息与物理信息；
[0027] (3)遍历获取的栅元,解析栅元所属集合和栅元是否被其它栅元集合所填充；
[0028] (4)将所有解析参数输入到核心转换模块中继续处理。
[0029] 所述可视化模块具体实现过程为：
[0030] (1)从转换核心模块参数树中获取不同部件填充关系；
[0031] (2)根据此种填充关系构建可视化分段渲染树结构；
[0032] ⑶实时渲染处于树结构中不同层级的分段，关闭不需要显示的分段。
[0033]本发明与现有技术相比的优点在于：
[0034] (1)本发明能够直观的得到堆芯的蒙特卡罗计算模型，避免过去手写蒙特卡罗计 算模型的繁琐性，提高了建模效率，极大减少了建模所消耗时间和发送错误的概率，将以往 数天才能完成并需要反复验证的建模任务，缩短至一个小时以内；
[0035] (2)本发明支持反复迭代设计，不仅仅是建模，然后计算，并且开发出了很多对模 型进行修改的功能，使整个反应堆设计过程一站化实施，不需要反复建模，增加了模型的重 用效率，消除了重复性的建模工作，使整个堆芯设计过程清晰明了，减轻设计者负担；
[0036] (3)本发明支持蒙特卡罗计算模型的解析和可视化，原有技术和方法，在得到蒙特 卡罗计算模型后，因为其大多数都由固定格式文本形式进行描述，很难想象其具体模型结 构，也很难在此基础上进行重用，本系统计算模型解析模块和可视化模块能够完美解决上 述问题，不仅能重用原有的蒙特卡罗计算模型，并且能够将其完整可视化，生成三维的面片 模型显示在用户面前，最大限度避免了出错的可能，因为肉眼能清晰观察到堆芯模型的结 构和任何细节；
[0037] (4)本发明可视化模块实现了对模型的分段渲染，这样避免一次性处理过多数据 导致可视化效率过低，又能在需要时候渲染出最底层细节，又避免了可视化图形由于面片 数量太多反而不能突出用户想要关注的重点等。

【专利附图】

【附图说明】
[0038] 图1为本发明总体流程图；
[0039]图2为本发明的组成框图；
[0040]图3为图1中参数控制模块的实现流程图；
[0041]图4为图1中计算模型解析模块的实现流程图；
[0042]图5为图1中核心转换模块的实现流程图；
[0043]图6为图1中可视化模块的实现流程图；
[0044]图7是堆芯组件设计参数示例；
[0045] 图8是参数树结构图；
[0046]图9是参数树合并过程；
[0047]图1〇是通过实施例1中表格参数，构建出来的堆芯CAD工程模型线框图；
[0048]图11是实施例1中生成的部分MCNP程序计算模型； '
[0049]图12是实施例2导入的部分MCNP程序计算模型；
[0050] 图13是通过导入的计算模型文件，生成的CAD工程模型线框图。

【具体实施方式】
[0051] 本发明能够实现在基于蒙特卡罗计算方法下对反应堆堆芯迭代设计的建模过程， 其能够生成堆芯的CAD工程设计模型和堆芯的蒙特卡罗计算模型，其中CAD工程设计模型 可以用于肉眼观察堆芯结构和工程设计，蒙特卡罗计算模型能够导入蒙特卡罗计算程序， 计算出反应堆堆芯关键物理参数，用以评判堆芯设计方案是否满足预期，如果没有满足，可 以在系统界面中进行可视化修改，直到满足设计预期为止。
[0052] 如图1所示，本发明总体流程为，首先可以解析己有蒙特卡罗计算模型，或者获取 用户界面输入得到参数信息，在通过参数构建完成参数树后，得到了栅兀建的层级信息，通 过层级信息，可以直接构建对应栅元部件的三维CAD几何，同时也可以直接输出蒙特卡罗 计算文件，在生成每个部件的三维CAD几何后，可以对所有的几何按照构建的分段信息进 行分段渲染，在界面中按照用户操作实时可视化不同的分段，用户可以在界面中使用系统 提供的相应功能对这些三维渲染模型进行修改操作，包括几何属性和物理属性，在修改完 成后可以生成新的蒙特卡罗计算模型，蒙特卡罗计算模型用于蒙特卡罗输运程序计算，最 终得到想要的结果。
[0053] 如图2所示，本发明系统包括参数控制模块、计算模型解析模块、核心转换模块和 可视化交互模块。参数控制模块负责获取用户的操作，并转换为模型的控制参数，计算模型 解析模块负责解析已有的蒙特卡罗计算模型，获取模型的参数，核心转换模块负责接收前 面两个模块得到的参数，并构建新的参数树或者修改参数树的部分节点，并负责将参数树 转换为CAD工程实体几何模型和新的蒙特卡罗计算模型，可视化模块负责接收到核心转换 模块得到的CAD工程实体几何，并更具参数树结构对对应的实体几何进行分段渲染，对整 个堆芯模型进行可视化。
[0054] 如图3所示，本发明中参数控制模块具体实现如下：
[0055] (1)用户通过界面输入堆芯内部包含的每一种组件的几何参数信息，例如形状、高 度、半径、密度、中子重要性、材料号，材料号对应具体的材料卡内部包含各种同位素和其不 同的核子密度。
[0056] (2)输入堆芯各个不同组件的装配信息，在径向上每一圈进行设置，例如前四圈是 燃料组件，后四圈是反射组件，最外层两圈是屏蔽组件。
[0057] (3)对燃料组件内部进一步进行细化，设计每一根燃料棒的结构，装备进燃料棒的 阵列。
[0058] (4)完成堆芯的所有参数设计，参数自动会被保存在核心转换模块数据结构中，核 心转换模块生成堆芯CAD模型，并在界面进行可视化，选中可视化图形，例如选中一根燃料 组件，选择组件替换功能，替换此组件为安全棒组件，在替换后核心转换模块中参数也同步 发生改变。
[0059]如图4所示，本发明中计算模型解析模块联合核心转换模块构建参数树具体实现 如下：
[0060] (1)导入一个已有的蒙特卡罗计算模型，例如压水堆模型的MCNP计算程序的输入 文件（即计算模型）。
[0061] (2)自动按照MCNP输入文件语法格式对整个输入文件进行解析，并按照模型之间 相互嵌套填充关系构建参数树，参数树的构建步骤如图8,例如根据MCNP计算模型语法规 则，先构找到一个U卡为父节点，将标记为此U卡的所有栅元放入此父节点的叶子节点，再 找到每个栅元里描述的填充信息，例如有一个叶子栅元的Fill卡为2,则U卡为2的父节点 为此叶子节点的叶子节点。
[0062] ⑶在完成所有栅元的遍历后，将得到很多棵这样的子树，最后根据子树父节点和 最下层叶子节点对应的其它子树的父节点连接起来，构成完整的参数结构。
[0063] (4)参数树自动保存在核心转换模块中。
[0064]如图5所示，本发明中核心转换模块具体实现如下：
[0065] (1)通过接收参数控制模块和计算模型解析模块的参数信息，构建好参数树，
[0066] (2)按照参数树从最底层开始向上层遍历，在得到一个叶子节点，按照蒙特卡 罗计算模型的语法规则输出此节点信息，例如按照一种新的国产超级蒙特卡罗计算程序 SuperMC的语法规则，类似于XML格式输出此节点。
[0067] ⑶通过参数树开始向上归并每个叶子节点里的几何体，例如发现2个平行面，则 手工补足缺失的另外4个面，4个面的面参数都默认给定一个足够大值，从而归并出一个立 方体，如果这个节点的上层节点还能找到另外方向的平行面，这自动将最新得到的平行面 与最开始的2个平行面组合起来，剩余的两个面仍然是由默认给定的足够大的值构建的面 组成，然后继续往上归并。
[0068] ⑷通过不断的归并，最终会得到每个部件的几何参数（如果描述的是一个正确 的反应堆，一定会达到这个结果），并能找到部件与部件的填充关系。
[0069] 在参数树结构中同时保留了不同部件的阵列信息，例如燃料棒在燃料组件内部阵 列，燃料组件在堆芯围桶内部阵列，通过这些信息核心转换模块开始构建实体几何，并根据 完成归并的参数树构建分段，将对应的实体几何划分进对应的分段之中。
[0070] 如图6所示，本发明中可视化模块具体实现如下：
[0071] (1)获取核心转换模块中构建的实体和分段信息；
[0072] (2)通过分段节点构建实体模型渲染单元；
[0073] (3)按照用户的需求显示渲染单元，呈现出不同的可视化面片模型给用户。
[0074] 如图7所示，为堆芯一个组件的几何参数实例，所有几何参数可以通过参数控制 模块提供的界面进行输入，在输入后会根据不同部件的结构自动完成参数树的构建。
[0075] 如图8所示，为一颗完成的参数树示意图，其实父节点表示一系列栅元的组合，在 此图中UX节点表示堆芯总体几何，堆芯总体几何由Cl、C2和C3三个栅元合并组成，UX节 点中记录了三个栅元的组成关系，其中C1栅元为堆芯内部燃料区域，由不同的燃料组件通 过阵列排布组成，第一种燃料组件标记为 υγ，第二种标记为UZ，其中UY和UZ又分别由其下 属叶子栅元节点组成，一直到参数树最底层，这样就形成了一个包括三层几何结构的完整 参数树。
[0076] 如图9所示，为参数树构建的方法，首先解析UX，UY和UZ三种栅元集合，得到三颗 参数子树，并得到UX下属栅元C1被UY和UZ填充，从而根据此种关系，可以合并三颗子树。
[0077] 实施例1 :
[0078] 目标是利用本发明，完成一个快中子反应堆堆芯的构建，其上下栅格参数如表L 燃料棒参数如表2,其它参数见表3,在完成建模，进行部分设计优化。
[0079] 表1 :上下栅格参数
[0080]

【权利要求】
1. 一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，其特征在于：包括核心转换模 块、参数控制模块、计算模型解析模块和可视化交互模块； 核心转换模块，构建一套满足反应堆堆芯设计参数体系，适用于反应堆堆芯部件存在 多层几何嵌套填充的特点，逻辑上为一种树形结构，本发明中统称为参数树，基于此参数树 核心转换模块完成三维CAD模型和蒙特卡罗计算程序输入模型的生成和相互转换；所述 参数控制模块与计算模型解析模块所提供参数全部输入到核心转换模块中，由核心转换模 块构建参数树，可视化交互模块从核心转换模块中提取参数树，基于参数树，可视化交互模 块对模型进行三维渲染，而参数控制模块能够继续对界面中渲染出的三维模型进行修改操 作，修改操作产生的新的参数会传递给核心转换模块，核心转换模块负责完成对参数树的 修改，在整个过程中核心转换模块随时能够将参数树中的信息转换成蒙特卡罗计算程序输 入文件，蒙特卡罗计算程序输入文件用于基于蒙特卡罗方法的输运计算，得到堆芯的关键 物理量，评价设计是否满足要求； 参数控制模块，实现通过系统提供修改操作和通过系统界面提供对话框输入新的参数 来构建、修改和控制堆芯的模型，并将所有参数输入到核心转换模块； 计算模型解析模块，自动识别、检测和拆分用户已有的反应堆堆芯蒙特卡罗计算模型， 最终得到堆芯各部件具体参数，将参数输入到核心转换模块之中； 可视化交互模块，渲染出通过核心转换模块构建好的计算机三维CAD模型，让用户实 时对模型的设计状态可见，该模块从转换核心模块参数树中获取不同部件填充关系，根据 此种填充关系构建可视化分段渲染树结构，实时渲染处于树结构中不同层级的分段，关闭 不需要显示的分段。
2. 根据权利要求1所述一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，其特征在 于：所述核心转换模块具体实现过程为： (1) 获取参数控制模块和计算模型解析模块参数，如果获取参数来源于参数控制模块 判定参数树是否为空，如果不为空，说明已经存在上一代参数树，则通过界面上存在限定约 束关系直接找到新的参数对应节点，在此节点中替换对应旧参数，完成对上一代参数树的 更新；如果为空，则根据界面上存在的约束关系直接到步骤（5);如果获取参数来源于是计 算模型解析模块，转至步骤（2)，（3)，（4); (2) 构建参数子树，参数子树中包含蒙特卡罗计算模型中限定的一组栅元集合，选此集 合的标号为根节点，集合内部的所有栅元为叶子节点，集合内部的单个栅元成为叶子栅元， 遍历所有叶子栅元，如果解析得到叶子栅元仍然被另一个栅元集合填充，则将填充此叶子 栅元所使用的栅元集合的标号作为此叶子栅元的再下一层子节点，从而完成单棵参数子树 构建； (3) 合并参数子树，通过各参数子树中根节点和最下层叶子节点信息，就能讲所有子树 链接起来，构建成一颗完整参数树； (4) 归并参数树，由于蒙特卡罗计算模型中栅元中几何有多个半空间面组成，几何信息 分散零碎，遍历参数树，归并参数树节点，合并零碎的半空间面，组合成更加完整几何； (5) 完成最终参数树构建，基于此参数树，根据部件之间的包含关系划分出对应的层级 属性； (6) 遍历参数树中的信息，通过几何信息构造三维几何实体，输入到可视化模块进行渲 染； (7)遍历参数树中的信息，通过参数树中划分的层级属性，按层级的描述堆芯的每个栅 元，包括几何描述和其他物理描述，最终输出保留了层级关系的蒙特卡罗计算模型，用以导 入蒙特卡罗计算程序进行计算。
3. 根据权利要求1所述一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，其特征在 于：所述参数控制模块具体实现过程为： (1) 通过扑捉用户在系统界面操作，获取用户输入参数和操作信息，构建堆芯各组件几 何，进一步细化堆芯组件内部燃料棒阵列，改变组件和组件内部元件的阵列规则，改变控制 棒位置，改变控制棒插入状态，燃料棒轴向生长，燃料棒径向膨胀和改变各部分几何对应物 理属性； (2) 从核心转换模块中获取已有参数树，如果参数树为空，则按照用户输入参数和操作 信息自动构建参数树； (3) 如果参数树不为空，则将用户新输入的控制参数和操作产生的结果反馈进参数树 之中，并继续保存在核心转换模块。
4. 根据权利要求1所述一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，其特征在 于：所述计算模型解析模块具体实现过程为： (1) 用户导入已有的蒙特卡罗计算模型； (2) 按照蒙特卡罗算模型语法规则，进行解析，得到蒙特卡罗计算模型中每一个栅元的 参数，栅元中包含几何信息与物理信息； (3) 遍历获取的栅元，解析栅元所属集合和栅元是否被其它栅元所填充； (4) 将所有解析参数输入到核心转换模块中继续处理。
5. 根据权利要求1所述一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统，其特征在 于：所述可视化模块具体实现过程为： (1) 从转换核心模块参数树中获取不同部件填充关系； (2) 根据此种填充关系构建可视化分段渲染树结构； (3) 实时渲染处于树结构中不同层级的分段，关闭不需要显示的分段。
【文档编号】G06T15/00GK104298836SQ201410625030
【公开日】2015年1月21日 申请日期:2014年11月6日 优先权日:2014年11月6日
【发明者】甘佺, 俞盛朋, 程梦云, 汪冬, 吴宜灿 申请人:中国科学院合肥物质科学研究院