(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211129779.8 (22)申请日 2022.09.16 (71)申请人 西北核技 术研究所 地址 710024 陕西省西安市灞桥区平峪路 28号 (72)发明人 商鹏　王学栋　朱金辉　刘利　 牛进林　 (74)专利代理 机构 西安智邦专利商标代理有限 公司 6121 1 专利代理师 徐秦中 (51)Int.Cl. G06F 30/25(2020.01) G06F 111/10(2020.01) G06F 119/02(2020.01) (54)发明名称 冲击波影响下的伽马辐射大气输运自适应 建模方法 (57)摘要 本发明为解决现有在冲击波影响条件下的 伽马辐射大气输运模拟中， 存在精细高效建模困 难的技术问题， 而提供了一种冲击波影 响下的伽 马辐射大气输运自适应 建模方法。 具体为： 首先， 根据大气密度分布函数建立模型并以给定步长 进行均匀粗网格剖分； 然后， 据大气密度的相对 变化程度进行网格细化插值； 其次， 在细化网格 中查找大气密度变化不大的节点进行合并； 最 后， 得到伽马辐射大气输运的精细高效自适应网 格模型， 再利用蒙特卡罗程序开展计算模拟。 本 发明适用于大气密度分布的自适应网格划分， 还 适用于对时间、 能量等进行自适应网格划分， 是 冲击波影响下的伽马辐射大气输运蒙特卡罗模 拟模型构建的新方法。 权利要求书2页 说明书7页 附图3页 CN 115422816 A 2022.12.02 CN 115422816 A 1.一种冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征在于， 包括以下步 骤： 步骤1)、 以爆点为圆心， 根据大气密度径向分布函数和时刻， 按照设定的给定步长S0对 整个区域进行均匀网格剖分， 得到多层均匀介质组合模型， 并计算各个网格节点处的大气 密度； 步骤2)、 跳转到第一个网格节点； 步骤3)、 计算第一个网格节点处与第二个网格节点处的大气密度相对误差； 步骤4)、 判断第一个网格节点与第 二个网格节点之间的网格间步长是否小于自适应网 格的最小分辨 率SN， 若是， 则跳转至下一个网格节点； 若否， 则执 行步骤5)； 步骤5)、 判断相对误差是否大于设定的最大相对误差A， 若是， 则在两个节点的中间位 置插入一个新的网格节点， 计算新插入网格节点处的密度值， 将该新插入网格节点作为第 二个网格节 点， 将原剖分的第二个网格节点及其之后的网格节点依次顺延， 并返回步骤3)， 直至第一个网格节点处与第二个网格节点处的大气密度相对误差小于最大相对误差A， 或 者第一个网格节点与第二个网格节点之间的网格间步长小于自适应网格的最小分辨率SN； 若否， 则跳转至下一个网格节 点； 步骤6)、 依次切换网格节 点， 采用步骤3) ‑步骤5)相同的方 法， 直至所有网格节点处与下一个相 邻的网格节 点处的大气密度相对误差小于最大相对误 差A， 或者与下一个相邻的网格节点之间的网格间步长小于自适应网格的最小分辨 率SN； 步骤7)、 按照所有网格节点的最终排序编号， 并再次跳转到第一个网格节点； 步骤8)、 计算第一个网格节点处与第二个网格节点处的大气密度相对误差， 判断相对 误差是否小于设定的最小相对误差B， 若是， 则将两个网格节点合并为一个网格节点， 删除 第二个网格节点， 将之后的网格节点依次按顺序重新编号， 执行步骤9)； 若否， 执行步骤 10)； 步骤9)、 返回步骤8)， 重新计算第一个网格节点处与下一个网格节点处的大气密度相 对误差， 直至第一个网格节点处与其相邻的网格节点处的大气密度相对误差大于最小相对 误差B； 步骤10)、 依次切换网格节点， 采用步骤8) ‑步骤9)相同的方法， 直至所有网格节点处与 其相邻的下一个网格节点处的大气密度相对误差大于设定的最小相对误差B； 步骤11)、 根据最终得到的所有自适应网格节点的位置， 建立伽马辐射大气输运的精细 高效自适应几何模型， 沿径向设置探测器， 进而开展伽 马辐射大气 输运的蒙特卡罗模拟， 计 算伽马辐射组织剂量。 2.根据权利要求1所述的冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征 在于： 步骤1)中， 所述选取的给定步长S0能够初步反应大气密度分布函数的变化趋势。 3.根据权利要求2所述的冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征 在于： 步骤4)中， 所述自适应网格最小分辨 率SN为： SN＝S0/2N 其中， N为自适应网格的最大分裂次数。 4.根据权利要求3所述的冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征权　利　要　求　书 1/2 页 2 CN 115422816 A 2在于， 步骤3)、 步骤8)中， 所述大气密度相对误差的计算过程 为： 第i个网格节点处与第i+1个网格节点处的大气密度相对误差 Erel为： 其中， i为网格节点的顺序编号， ρi为第i个网格节点处的大气密度。 5.根据权利要求1 ‑4任一所述的冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征在于： 步骤5)中， 所述 最大相对误差A＝5％； 步骤9)中， 最小相对误差B＝0.1％。 6.根据权利要求5所述的冲击波影响下的伽马辐射大气输运自适应建模方法， 其特征 在于： 步骤1)中， 所述大气密度径向分布函数为连续 函数或分段函数。权　利　要　求　书 2/2 页 3 CN 115422816 A 3