OpenMC StatePoint无Tally问题排查及验证代码请求

问题背景 #

使用OpenMC Python绑定进行中子输运仿真，已通过Python定义Tally并生成geometry.xml、settings.xml、materials.xml、tallies.xml，但运行后生成的statepoint.h5中无Tally对象（sp.tallies_present返回False），以下针对疑问逐一解答：

1. 正确定义Tally的要求是什么？ #

• 必须指定得分类型（scores）：如flux（通量）、absorption（吸收）、fission（裂变）等，无得分类型的Tally不会被执行
• 必须配置过滤条件（filters）：明确Tally的计算范围，比如通过CellFilter指定单元格、SurfaceFilter指定表面、EnergyFilter指定能区间，无过滤条件的Tally无法确定计算对象
• 需将Tally添加到模型并导出：Python代码中需通过model.tallies.append(tally)将Tally加入模型，再执行model.export_to_xml()生成合法的tallies.xml；手动编写XML时，需保证<tally>标签下包含<scores>、<filters>等必要子标签，结构符合规范
• 可选但建议设置Tally名称：便于后续结果识别，避免默认ID带来的混淆

2. 存在tallies.xml但StatePoint无Tally的情况有哪些？ #

• tallies.xml格式错误：标签拼写错误（如<score>代替<scores>）、标签嵌套错误，导致OpenMC解析时忽略无效Tally
• Tally过滤条件无效：指定的单元格ID/表面ID不存在，或过滤参数错误（如EnergyFilter的能区间不合法），OpenMC无法定位计算对象，直接跳过该Tally
• 仿真配置关闭Tally输出：settings.xml中设置output = {'tallies': False}，或statepoint的tallies参数设为False，导致Tally结果不写入StatePoint
• 仿真未正常完成：StatePoint文件为中途中断生成，未执行到Tally结果写入步骤
• 得分类型不兼容：如临界计算中指定current得分但未设置合适的边界条件，或得分类型拼写错误（如flx代替flux），导致Tally无法执行

3. 排查Tally配置错误的方法 #

• 查看OpenMC运行日志：日志会显示Tally加载情况，如Loaded X tallies from tallies.xml，若X为0说明XML解析失败；若有警告信息（如Invalid filter for tally 1），可针对性排查对应Tally
• 用Python加载模型验证：执行model = openmc.Model.from_xml()后，打印model.tallies，若为空则说明tallies.xml解析失败；若有Tally对象，检查其scores、filters属性是否正确
• 手动校验tallies.xml结构：对照规范确保每个<tally>标签包含必要子元素，标签名、属性拼写无误
• 极简测试法：先编写仅计算全局通量的极简Tally（如示例代码），验证能否正常写入StatePoint，再逐步添加复杂配置排查问题

4. 通用测试用Tally代码片段与XML #

Python代码片段（可嵌入任意模型）

import openmc

# 创建全局通量Tally
tally = openmc.Tally(name="global_flux")
# 指定得分类型为通量
tally.scores = ["flux"]
# 应用于所有单元格
tally.filters = [openmc.CellFilter(openmc.Cell.all)]

# 假设已定义model的materials、geometry、settings
model = openmc.Model()
# 将Tally加入模型并导出XML
model.tallies.append(tally)
model.export_to_xml()

对应的tallies.xml（手动编写版）

<?xml version='1.0' encoding='UTF-8'?>
<tallies>
  <tally id="1" name="global_flux">
    <scores>
      <score>flux</score>
    </scores>
    <filters>
      <filter type="cell" bins="all" />
    </filters>
  </tally>
</tallies>

内容的提问来源于stack exchange，提问作者PackSciences