核工程仿真软件：OpenMC二次开发_（6）.源定义与输运模拟.docx

PAGE1

PAGE1

源定义与输运模拟

在核工程仿真软件OpenMC中，源定义与输运模拟是进行核反应堆分析和中子输运计算的关键步骤。源定义决定了仿真中的中子或光子的初始位置、方向和能量分布，而输运模拟则负责追踪这些粒子在几何结构中的运动轨迹，直至它们被吸收或逃逸出系统。本节将详细介绍如何在OpenMC中定义源，并进行输运模拟。

源定义

1.初始源分布

初始源分布是指在仿真开始时，中子或光子的初始位置、方向和能量分布。OpenMC提供了多种方式来定义初始源分布，包括点源、体源、表面源和文件源。

点源

点源是最简单的源定义方式，用于模拟中子或光子从一个特定点发射的情况。定义点源时，需要指定中子的初始位置、方向分布和能量分布。

#导入OpenMC库

importopenmc

#定义点源的位置

point=openmc.stats.Point((0,0,0))

#定义中子的方向分布

direction=openmc.stats.Isotropic()

#定义中子的能量分布（例如，2MeV的单能分布）

energy=openmc.stats.Discrete([2.0e6],[1.0])

#创建点源对象

source=openmc.Source(space=point,angle=direction,energy=energy)

#将源添加到设置文件中

settings=openmc.Settings()

settings.source=source

settings.run_mode=eigenvalue#设置运行模式为特征值计算

settings.batches=100#设置批次数

settings.inactive=10#设置非活动批次数

settings.particles=1000#每批粒子数

#保存设置文件

settings.export_to_xml()

2.体源

体源用于模拟中子或光子在特定几何体内的均匀分布。定义体源时，需要指定几何体的范围和粒子的初始位置、方向和能量分布。

#导入OpenMC库

importopenmc

#定义体源的范围

box=openmc.stats.Box((-1,-1,-1),(1,1,1),only_fissionable=True)

#定义中子的方向分布（各向同性）

direction=openmc.stats.Isotropic()

#定义中子的能量分布（例如，热中子分布）

energy=openmc.stats.Maxwellpectrum(temperature=900.0)

#创建体源对象

source=openmc.Source(space=box,angle=direction,energy=energy)

#将源添加到设置文件中

settings=openmc.Settings()

settings.source=source

settings.run_mode=fixedsource#设置运行模式为固定源计算

settings.batches=100#设置批次数

settings.particles=1000#每批粒子数

#保存设置文件

settings.export_to_xml()

3.表面源

表面源用于模拟中子或光子从特定表面发射的情况。定义表面源时，需要指定表面的类型、位置和粒子的初始方向和能量分布。

#导入OpenMC库

importopenmc

#定义表面源的平面

plane=openmc.Plane(a=1,b=0,c=0,d=0)

#定义中子的方向分布（垂直于平面）

direction=openmc.stats.Discrete([1,0,0],[1.0])

#定义中子的能量分布（例如，250keV的单能分布）

energy=openmc.stats.Discrete([250.0e3],[1.0])

#创建表面源对象

surface_source=openmc.Source(space=openmc.stats.SurfaceSource([plane]),angle=direction,energy=energy)

#将源添加到设置文件中

settings=openmc.Settings()

settings.source=surf