11. 可视化

CRANE内建有非常丰富的可视化功能,这些功能根据内容可以分为建模的可视化和计算结果的可视化两类。

Note

CRANE内建的可视化功能都是基于两维矢量图的,如果你对堆芯三维可视化有需求,可关注CRANE的进展, 我们正在研发基于web的堆芯三维可视化工具。

为了能够方便用户使用,所有可视化命令都采用交互式命令行方式,只要执行 crane visualize 命令进入到交互式命令行环境, 按 Tab 键就可提示出可视化命令及选项。

对于建模可视化,crane visualize 后面可以是 用户输入文件,也可以是 HDF5 结果文件。但如果要进行计算结果可视化,crane visualize 后面必须是 HDF5 结果文件

比如,你可以执行以下命令

crane visualize VERA_7.yaml

或者

crane visualize VERA_7.h5

进入到 VERA_7 问题的交互式可视化环境,前者只支持建模可视化,而后者建模和计算结果可视化都支持。

除此之外,CRANE还支持对不同计算结果的比较可视化功能,如执行以下命令:

crane visualize VERA_7.h5 VERA_7_another.h5

或者:

crane visualize VERA_7.h5 VERA_7_another.xlsx

可进入到 VERA_7 问题的两个计算结果比较的可视化环境,结果比较的偏差是前者减去后者,更为详细的说明可参阅 结果比较可视化

以下我们以 VERA_7 基准题为例,对所有可视化功能进行逐一介绍。

11.1. 模型可视化

模型可视化通过 plot-model 命令实现,是最基本的建模可视化功能,该功能会对不同材料渲染不同颜色来绘制模型两维图, 支持径向、轴向两种切面。用户可以通过模型可视化来检查几何建模是否正确,材料区输入是否正确等。

执行以下交互式命令行:

plot-model --direction radial --axial-height 200

Note

--direction 选项用来选择切面类型,默认为径向切面,--axial-height 选项用来选择径向切面的高度。

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行可视化:

../_images/plot_model_radial.png

将选项 --axial-height 设置为 422

plot-model --direction radial --axial-height 422

可对距离堆芯底部(包含反射层)距离422cm处的径向切面进行可视化:

../_images/plot_model_radial_422.png

可见这个高度是活性区之外了,没有燃料栅元,但能够看到导向管、以及插入导向管中的可燃毒物棒和控制棒。

Note

这里如果没有 --axial-height 选项,则会对所有轴向层从底到顶逐一可视化,过程中可以按 Ctrl+C 终止。

执行以下交互式命令行:

plot-model --direction axial --axial-slice-coef 1 1 0

Note

--axial-slice-coef 选项用来选择轴向切面的方程系数。

可对方程为 x+y=0 的轴向切面进行可视化:

../_images/plot_model_axial.png

11.2. 网格可视化

网格可视化通过 plot-meshes 命令实现,该可视化功能会对不同的网格渲染不同的颜色,支持材料区、MOC、 CMFD、共振计算、燃耗计算、和热工水力通道共六种类型的网格。

所有网格可视化都支持径向、轴向两类切面,以下以径向切面网格可视化为例。

11.2.1. 材料区网格可视化

--mesh-type 设置为 material 时则是对材料区网格的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type material --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行材料区网格的可视化:

../_images/plot_meshes_material_radial.png

放大至组件规模:

../_images/plot_meshes_material_radial_zoom.png

可见材料网格可视化是根据材料区来渲染颜色的。

11.2.2. MOC网格可视化

--mesh-type 设置为 moc 时则是对MOC网格的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type moc --axial-height 200

Note

由于MOC网格数量较多,首次绘制时会需要些时间,后续绘制都会使用缓存,速度不会受到网格数量影响。

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行MOC网格的可视化:

../_images/plot_meshes_moc_radial.png

放大至组件规模:

../_images/plot_meshes_moc_radial_zoom.png

可见MOC网格可视化是根据MOC网格来渲染颜色的,MOC网格是整个模拟计算中最细的网格。

11.2.3. CMFD网格可视化

--mesh-type 设置为 cmfd 时则是对CMFD网格(栅元网格)的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type cmfd --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行CMFD网格的可视化:

../_images/plot_meshes_cmfd_radial.png

放大至组件规模:

../_images/plot_meshes_cmfd_radial_zoom.png

可见CMFD网格可视化是根据CMFD网格(在CRANE中就是栅元)来渲染颜色的。

11.2.4. 共振区网格可视化

--mesh-type 设置为 resonance 时则是对共振区网格的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type resonance --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行共振区网格的可视化:

../_images/plot_meshes_resonance_radial.png

放大至组件规模:

../_images/plot_meshes_resonance_radial_zoom.png

可见共振区网格可视化是根据共振区来渲染颜色的,非共振区不渲染颜色。

11.2.5. 燃耗区网格可视化

--mesh-type 设置为 depletion 时则是对燃耗区网格的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type depletion --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行燃耗区网格的可视化:

../_images/plot_meshes_depletion_radial.png

放大至组件规模:

../_images/plot_meshes_depletion_radial_zoom.png

可见燃耗区网格可视化是根据燃耗区来渲染颜色的,非燃耗区不渲染颜色。

11.2.6. 热工水力通道网格可视化

--mesh-type 设置为 channel 时则是对热工水力通道网格的可视化。

执行以下交互式命令行:

plot-meshes --mesh-type channel --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行热工水力通道网格的可视化:

../_images/plot_meshes_channel_radial.png

放大至组件规模:

../_images/plot_meshes_channel_radial_zoom.png

可见热工水力通道网格可视化是根据热工水力计算的通道(在CRANE中是栅元通道)来渲染颜色的,非通道位置不渲染颜色。

11.3. 组件布置图可视化

组件布置图可视化通过 plot-assembly-map 命令实现,该可视化功能以组件为单位绘制堆芯布置图,支持组件装载、 组件历史、组件类型、组件编号、控制棒组、分离式可燃毒物、中子测量管共七种。

这里我们用多循环计算基准题 BEAVRS_C2 作为示例来进行组件布置图的可视化功能。

11.3.1. 组件装载布置图

--map-type 设置为 loading_pattern 时则是组件装载布置图。

执行以下交互式命令行:

plot-assembly-map --map-type loading_pattern

可绘制出以下组件装载布置图:

../_images/plot_assembly_map_lp.png

这里有颜色渲染的组件是旧料,且根据燃耗过的循环数来渲染不同颜色。

11.3.2. 组件历史布置图可视化

--map-type 设置为 assembly_history 时则是组件历史布置图。

执行以下交互式命令行:

plot-assembly-map --map-type assembly_history

可绘制出以下组件历史布置图:

../_images/plot_assembly_map_his.png

可见组件历史布置图与组件装载布置图有点类似,主要的区别是旧料组件会把该组件经历的所有循环的详细信息(每个循环所在位置, 燃耗)都打印出来。

11.3.3. 组件类型布置图可视化

--map-type 设置为 assembly_type 时则是组件类型布置图。

执行以下交互式命令行:

plot-assembly-map --map-type assembly_type

可绘制出以下组件类型布置图:

../_images/plot_assembly_map_type.png

可见组件类型布置图是根据组件的类型来渲染不同颜色,而不区分该组件是新料还是旧料。

11.3.4. 组件编号布置图可视化

--map-type 设置为 assembly_index 时则是组件装载布置图。

执行以下交互式命令行:

plot-assembly-map --map-type assembly_index

可绘制出以下组件编号布置图:

../_images/plot_assembly_map_index.png

方形组件堆型的组件编号布置图是结合组件行、列编号得到的,而六角形堆型的组件则按从底到顶、 从左到右的顺序来对组件进行编号,如下图为VVER-1000的组件编号布置图

../_images/plot_assembly_map_index_vver.png

11.3.5. 控制棒组布置图可视化

--map-type 设置为 RCCA 时则是控制棒组布置图。

执行以下交互式命令行:

plot-assembly-map --map-type RCCA

可绘制出以下控制棒组布置图:

../_images/plot_assembly_map_rcca.png

可见控制棒组件布置图是以棒组名称来进行插棒组件不同颜色的渲染。

11.3.6. 分离式可燃毒物布置图可视化

--map-type 设置为 inserted_BP 时则是分离式可燃毒物布置图。

执行以下交互式命令行:

plot-assembly-map --map-type inserted_BP

可绘制出以下分离式可燃毒物布置图:

../_images/plot_assembly_map_bp.png

可见分离式可燃毒物布置图是以分离式可燃毒物在组件内的不同布置来进行不同颜色的渲染。

11.3.7. 中子测量管布置图可视化

--map-type 设置为 inserted_IT 时则是中子测量管布置图。

执行以下交互式命令行:

plot-assembly-map --map-type inserted_IT

可绘制出以下中子测量管布置图:

../_images/plot_assembly_map_it.png

11.4. 标记区域可视化

标记区域可视化通过 plot-tagged-regions 命令实现,该可视化功能只对标记区域渲染颜色, 用来检查用户所标记的区域是否正确,有关标记区域的功能请参阅 tally_tagged_regions

与模型可视化一样,标记区域可视化也支持径向、轴向两类切面。

执行以下交互式命令行:

plot-tagged-regions --axial-height 200

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行标记区域的可视化,这里我们对围板、 围桶等反射层区域进行了标记:

../_images/plot_tagged_regions_radial.png

11.5. MOC特征线可视化

MOC特征线可视化通过 plot-moc-tracks 命令实现,该可视化功能在模型径向切面基础上绘制出不同层级的特征线, 支持堆芯、组件、特征线段三个层级,其 --track-level 选项分别为 core, assembly, segments

Note

层级表示特征线贯穿的区域,如组件特征线指的是一个端点在组件一边,另一个端点在组件另一边上的特征线。

执行以下交互式命令行:

plot-moc-tracks --axial-height 200 --track-level assembly

Note

segments 层级(贯穿MOC细网)的特征线太多,在大规模问题上绘制该层级特征线会非常耗时。

可对距离堆芯底部(包含反射层)距离200cm处的径向切面进行组件层级的特征线可视化。

../_images/plot_moc_tracks_assembly.png

由于特征线太密,只看到黑色,当放大至栅元规模时才能看清特征线:

../_images/plot_moc_tracks_assembly_zoom.png

11.6. 功率分布可视化

功率分布可视化通过 plot-power-distribution 命令实现,能够对以下 8 种功率分布进行可视化:

11.6.1. 组件积分功率分布可视化

组件积分功率分布为最常用的功率分布,是对每个组件进行轴向积分后的两维功率分布,其 --power-type 选项为 assembly-integral

执行以下交互式命令行:

plot-power-distribution --power-type assembly-integral

可绘制出组件积分功率分布图:

../_images/plot_power_distribution_assembly_integral.png

11.6.2. 全堆组件轴向功率偏移分布可视化

该命令能够将全堆每个组件的轴向功率偏移进行可视化。 其 --power-type 选项为 assembly-axial-offset-full-core

执行以下交互式命令行:

plot-power-distribution --power-type assembly-axial-offset-full-core

可绘制出全堆组件轴向功率偏移分布图:

../_images/plot_axial_offset_distribution_assembly_full_core.png

11.6.3. 全堆组件积分功率分布可视化

通常我们会根据对称性做部分堆芯计算,而功率分布图很多时候需要全堆芯的,因此该功能是对组件积分功率分布根据对称性进行补全。 其 --power-type 选项为 assembly-integral-full-core

执行以下交互式命令行:

plot-power-distribution --power-type assembly-integral-full-core

可绘制出全堆组件积分功率分布图:

../_images/plot_power_distribution_assembly_integral_full_core.png

11.6.4. 以组件为单位的轴向功率分布可视化

以组件为单位的轴向功率分布,即把每个组件的轴向(层积分)功率分布绘制在一张图上,以比较不同组件轴向功率分布的特点。 其 --power-type 选项为 assembly-axial

执行以下交互式命令行:

plot-power-distribution --power-type assembly-axial

可绘制出以组件为单位的轴向功率分布图:

../_images/plot_power_distribution_assembly_axial.png

Note

该图在可视化客户端下只能对单个子图进行放大,建议增加 --save-fig 选项保存图片后放大查看。

11.6.5. 径向切面的组件功率分布可视化

径向切面的组件功率分布是在三维组件功率分布下,选取某径向切面绘制的两维功率分布。其 --power-type 选项为 assembly-radial-slice

可通过 --map-color-by 来选择是以当前层还是全堆的功率数值范围来渲染颜色。

执行以下交互式命令行:

plot-power-distribution --power-type assembly-radial-slice --axial-height 200

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的组件功率分布图:

../_images/plot_power_distribution_assembly_radial_slice.png

11.6.6. 棒积分功率分布可视化

棒积分功率分布,是对每根燃料棒的进行轴向积分后的两维功率分布,其 --power-type 选项为 pin-integral

执行以下交互式命令行:

plot-power-distribution --power-type pin-integral

可绘制出棒积分功率分布图:

../_images/plot_power_distribution_pin_integral.png

11.6.7. 径向切面的棒功率分布可视化

径向切面的棒功率分布是在三维棒功率分布下,选取某径向切面绘制的两维功率分布。其 --power-type 选项为 pin-radial-slice

可通过 --map-color-by 来选择是以当前层还是全堆的功率数值范围来渲染颜色。

执行以下交互式命令行:

plot-power-distribution --power-type pin-radial-slice --axial-height 200

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的棒功率分布图:

../_images/plot_power_distribution_pin_radial_slice.png

11.6.8. 轴向切面的棒功率分布可视化

轴向切面的棒功率分布是在三维棒功率分布下,选取某轴向切面绘制的两维功率分布。其 --power-type 选项为 pin-axial-slice

执行以下交互式命令行:

plot-power-distribution --power-type pin-axial-slice --axial-slice-coef 1 1 0

可绘制出方程为 x+y=0 轴向切面的棒功率分布图:

../_images/plot_power_distribution_pin_axial_slice.png

11.6.9. 堆芯轴向功率分布可视化

堆芯轴向功率分布,是对层积分后的整个堆芯的轴向功率分布,其 --power-type 选项为 axial

执行以下交互式命令行:

plot-power-distribution --power-type axial

可绘制出堆芯轴向功率分布图:

../_images/plot_power_distribution_axial.png

Note

默认会把堆芯轴向切面图绘制在左边,以方便观察功率分布对应的材料区域。

11.7. 通量分布可视化

通量分布可视化通过 plot-flux-distribution 命令实现,能够对以下 5 种功率分布进行可视化:

11.7.1. 径向切面的栅元通量分布可视化

径向切面的栅元通量分布是指在三维栅元通量分布下,选取某径向切面绘制的两维通量分布。其 --flux-type 选项为 cell-radial-slice

执行以下交互式命令行:

plot-flux-distribution --flux-type cell-radial-slice --axial-height 200 --energy thermal

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的栅元热通量分布图:

../_images/plot_flux_distribution_cell_radial_slice_thermal.png

通过 --energy 选项可以选择中子通量的能量范围,支持 thermal|epithermal 两个选项,默认为 thermal

Note

thermalepithermal 的能量分界线为 0.625eV

选择绘制非热群通量分布:

plot-flux-distribution --flux-type cell-radial-slice --axial-height 200 --energy epithermal

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的栅元非热通量分布图:

../_images/plot_flux_distribution_cell_radial_slice_epithermal.png

11.7.2. 轴向切面的栅元通量分布可视化

轴向切面的栅元通量分布是指在三维栅元通量分布下,选取某轴向切面绘制的两维通量分布。其 --flux-type 选项为 cell-axial-slice

执行以下交互式命令行:

plot-flux-distribution --flux-type cell-axial-slice --axial-slice-coef 1 1 0 --energy thermal

可绘制出方程为 x+y=0 轴向切面的栅元热通量分布图:

../_images/plot_flux_distribution_cell_axial_slice.png

11.7.3. 径向切面的细网通量分布可视化

径向切面的细网通量分布是指在三维细网(MOC计算网格)通量分布下,选取某径向切面绘制的两维通量分布。其 --flux-type 选项为 mesh-radial-slice

Note

绘制细网通量分布的前提是必须保存计算结果文件必须保存MOC计算网格通量数据,由于MOC计算网格通量数据量较大, 只有当输出类型 type 设置为 MAXIMUM 才会保存。

执行以下交互式命令行:

plot-flux-distribution --flux-type mesh-radial-slice --axial-height 200 --energy thermal

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的细网热群通量分布图:

../_images/plot_flux_distribution_mesh_radial_slice.png

我们选取一个4x4栅元区域并放大,可以看到细网通量分布是基于MOC计算网格的,比栅元网格更加精细:

../_images/plot_flux_distribution_mesh_radial_slice_zoom.png

11.7.4. 轴向切面的细网通量分布可视化

轴向切面的细网(MOC计算网格)通量分布是指在三维细网通量分布下,选取某轴向切面绘制的两维通量分布。其 --flux-type 选项为 mesh-axial-slice

执行以下交互式命令行:

plot-flux-distribution --flux-type mesh-axial-slice --axial-slice-coef 1 1 0 --energy thermal

可绘制出方程为 x+y=0 轴向切面的细网热通量分布图:

../_images/plot_flux_distribution_mesh_axial_slice.png

我们选择有格架且高度为205cm处的一个区域,放大后可以看到比栅元更精细的轴向切面网格热通量分布:

../_images/plot_flux_distribution_mesh_axial_slice_zoom.png

11.7.5. 堆芯轴向通量分布可视化

类似于 以组件为单位的轴向功率分布可视化,堆芯轴向通量分布是指对层积分后的整个堆芯的轴向通量分布。 其 --flux-type 选项为 axial

执行以下交互式命令行:

plot-flux-distribution --flux-type axial --energy thermal

可绘制出堆芯轴向热通量分布图:

../_images/plot_flux_distribution_axial.png

Note

默认会把堆芯轴向切面图绘制在左边,以方便观察功率分布对应的材料区域。

由于格架对热中子的吸收,可以看到在格架位置处,热通量分布曲线有明显凹陷。

11.8. 探测器活度分布可视化

探测器活度分布可视化通过 plot-detector-reaction-distribution 命令实现, 该功能对以下 3 组件探测器的活度分布进行可视化:

Note

绘制探测器活度分布两维图的前提是要有探测器活度分布的计算结果,也就是用户输入文件中要有对探测器建模的所有数据, 即 instrument_bundledetector 对象 等。

11.8.1. 组件积分探测器活度分布可视化

该功能根据组件轴向积分后的探测器活度(反应率)进行颜色渲染产生两维图。其 --reaction-type 选项为 assembly-integral

执行以下交互式命令行:

plot-detector-reaction-distribution --reaction-type assembly-integral

Note

由于探测器一般会非对称布置,因此当CRANE进行对称性堆芯计算时(如1/4旋转对称堆芯), 全堆探测器活度分布是根据对称性恢复出来的。

可绘制出全堆组件积分的探测器活度分布图:

../_images/plot_detector_reaction_distribution_assembly_integral.png

11.8.2. 以组件为单位的轴向探测器活度分布可视化

以组件为单位的轴向探测器活度分布,即把每个组件的轴向(层积分)探测器活度分布绘制在一张图上,以比较不同组件轴向探测器活度分布的特点。 没有探测器组件位置则为空,其 --reaction-type 选项为 assembly-axial

执行以下交互式命令行:

plot-detector-reaction-distribution --reaction-type assembly-axial

可绘制出以组件为单位的轴向探测器活度分布图:

../_images/plot_detector_reaction_distribution_assembly_axial.png

Note

该图在可视化客户端下只能对单个子图进行放大,建议增加 --save-fig 选项保存图片后放大查看。

11.8.3. 径向切面组件探测器活度分布可视化

径向切面的组件探测器活度分布是在三维组件探测器活度分布下,选取某一层绘制的两维探测器活度分布。其 --reaction-type 选项为 assembly-radial-slice,可通过 --plane-index 来选择画哪一层,该选项缺省情况下默认为所有层都进行可视化。

执行以下交互式命令行:

plot-detector-reaction-distribution --reaction-type assembly-radial-slice --plane-index 0

可绘制出第0层(最底层)探测器活度分布图:

../_images/plot_detector_reaction_distribution_assembly_radial_slice_plane_0.png

Note

该探测器活度分布图的颜色是以全堆三维数据来渲染的。

设置 --plane-index 为18:

plot-detector-reaction-distribution --reaction-type assembly-radial-slice --plane-index 18

可绘制出第18层探测器活度分布图:

../_images/plot_detector_reaction_distribution_assembly_radial_slice_plane_18.png

11.9. 燃耗分布可视化

燃耗分布可视化通过 plot-burnup-distribution 命令实现,能够对以下 5 种燃耗分布进行可视化:

由于 VERA_7 不带燃耗计算,这里我们以 VERA_9 基准题为例来进行燃耗分布可视化。

11.9.1. 全堆组件积分燃耗分布可视化

全堆组件积分燃耗分布为最常用的燃耗分布,是对每个组件进行轴向积分后的两维燃耗分布,其 --burnup-type 选项为 assembly-integral-full-core

执行以下交互式命令行:

plot-burnup-distribution --burnup-type assembly-integral-full-core --efpd 441

可绘制出在441个等效满功率天时的全堆组件积分燃耗分布图:

Note

可以通过 --burnup 或者 --epfd 选项来选择哪个状态点,下同。

_images/plot_burnup_distribution_assembly_integral_full_core.png

11.9.2. 组件积分燃耗分布可视化

同上,只对计算区域的组件积分燃耗进行可视化,其 --burnup-type 选项为 assembly-integral

执行以下交互式命令行:

plot-burnup-distribution --burnup-type assembly-integral --efpd 441

可绘制出在441个等效满功率天时的组件积分燃耗分布图:

Note

可以通过 --burnup 或者 --epfd 选项来选择哪个状态点,下同。

../_images/plot_burnup_distribution_assembly_integral.png

11.9.3. 棒积分燃耗分布可视化

棒积分燃耗分布,是对每根燃料棒的进行轴向积分后的两维燃耗分布,其 --burnup-type 选项为 pin-integral

执行以下交互式命令行:

plot-burnup-distribution --burnup-type pin-integral --efpd 441

可绘制出在441个等效满功率天时的棒积分燃耗分布图:

../_images/plot_burnup_distribution_pin_integral.png

11.9.4. 径向切面的棒燃耗分布可视化

径向切面的棒燃耗分布是在三维棒燃耗分布下,选取某径向切面绘制的两维燃耗分布。其 --burnup-type 选项为 pin-radial-slice

可通过 --map-color-by 来选择是以当前层还是全堆的燃耗数值范围来渲染颜色。

执行以下交互式命令行:

plot-burnup-distribution --burnup-type pin-radial-slice --axial-height 200 --efpd 441

可绘制出在441个等效满功率天时距离堆芯底部(包含反射层)距离200cm处径向切面的棒燃耗分布图:

../_images/plot_burnup_distribution_pin_radial_slice.png

11.9.5. 轴向切面的棒燃耗分布可视化

轴向切面的棒燃耗分布是在三维棒燃耗分布下,选取某轴向切面绘制的两维燃耗分布。其 --burnup-type 选项为 pin-axial-slice

执行以下交互式命令行:

plot-burnup-distribution --burnup-type pin-axial-slice --axial-slice-coef 1 1 0 --efpd 441

可绘制出在441个等效满功率天时方程为 x+y=0 轴向切面的棒燃耗分布图:

../_images/plot_burnup_distribution_pin_axial_slice.png

11.9.6. 堆芯轴向燃耗分布可视化

堆芯轴向燃耗分布,是对层积分后的整个堆芯的轴向燃耗分布,其 --burnup-type 选项为 axial

执行以下交互式命令行:

plot-burnup-distribution --burnup-type axial --efpd 441

可绘制出在441个等效满功率天时堆芯轴向燃耗分布图:

../_images/plot_burnup_distribution_axial.png

Note

默认会把堆芯轴向切面图绘制在左边,以方便观察燃耗分布对应的材料区域。

11.10. 核素核子密度分布可视化

核素核子密度可视化通过 plot-number-density-distribution 命令实现, 该功能根据指定核素每个燃耗区的核子密度进行颜色渲染产生两维图,支持径向、轴向两种切面。 用户可以通过核素核子密度可视化来查看核素成分在堆芯内的空间分布。

Note

绘制核子密度分布的前提是计算结果文件中必须保存各个燃耗区的核素核子密度数据,由于该数据量较大, 只有当输出类型 type 设置为 MAXIMUM 才会保存。

执行以下交互式命令行:

plot-number-density-distribution --nuclide-name Xe-135 --direction radial --axial-height 200

Note

目前支持的核子密度可视化的核素可查看 plot-number-density-distribution--nuclide-name 选项。

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的Xe-135核子密度分布图:

../_images/plot_number_density_distribution_radial.png

或者可以执行:

plot-number-density-distribution --nuclide-name Xe-135 --direction axial --axial-slice-coef 1 1 0

绘制出方程为 x+y=0 轴向切面的Xe-135核子密度分布图:

../_images/plot_number_density_distribution_axial.png

11.11. Dancoff因子分布可视化

Dancoff因子可视化通过 plot-dancoff-factors 命令实现,该功能以每个共振区的Dancoff因子进行颜色渲染产生两维图, 支持径向、轴向两种切面。

Note

由于Dancoff因子是能量相关的,如果保存所有共振能群的Dancoff因子会导致存储量较大, 因此这里的Dancoff因子是每个共振区不同能群的Dancoff因子的一个平均值。

用户可以通过Dancoff因子分布可视化来查看共振互屏效应,Dancoff因子越小,共振互拼效应越强。

Note

绘制Dancoff因子分布的前提是计算结果文件中必须保存各个共振区的Dancoff因子,由于该数据量较大, 只有当输出类型 type 设置为 MAXIMUM 才会保存。

执行以下交互式命令行:

plot-dancoff-factors --direction radial --axial-height 200

可绘制出距离堆芯底部(包含反射层)距离200cm处径向切面的Dancoff因子分布图:

../_images/plot_dancoff_factors_radial.png

选择最外围的一个组件区域放大:

../_images/plot_dancoff_factors_radial_zoom.png

可见最外围燃料棒的Dancoff因子最大,因为更多慢化剂减弱了燃料棒之间的共振互屏效应。

或者可以执行:

plot-dancoff-factors --direction axial --axial-slice-coef 1 1 0

绘制出方程为 x+y=0 轴向切面的Dancoff因子分布图:

../_images/plot_dancoff_factors_axial.png

选择高度为205cm左右包含有格架的一个区域放大:

../_images/plot_dancoff_factors_axial_zoom.png

同样可看到,共振区附近慢化剂越多,其Dancoff因子越大,共振互屏效应越小。

11.12. 汇总数据随燃耗变化可视化

汇总数据随燃耗可视化通过 plot-result-summary 命令实现,该功能以燃耗(或时间)为X轴,汇总数据( /statepoint<index>/summary/ 中的数据,如Keff、临界硼浓度等)为Y轴绘制折线(或散点图)。

这里我们以 VERA_9 基准题为例来进行汇总数据随燃耗可视化。

执行以下交互式命令行:

plot-result-summary --x-axis effective_full_power_days --y-axis-first boron_concentration

可绘制出临界硼浓度随等效满功率天的变化折线图:

../_images/plot_result_summary_efpd_bor.png

也可以通过 --plot-type-first 设置为 scatter 将对应的数据绘制成散点图:

../_images/plot_result_summary_efpd_bor_scatter.png

可支持同时绘制最多 3 个Y轴的数据,如执行以下交互式命令行:

plot-result-summary --y-axis-first boron_concentration --y-axis-second control_rod_withdrawn_steps --rod-bank D --y-axis-third relative_power

Note

这里 control_rod_withdrawn_steps 需要 --rod-bank 来指定是哪一组棒,默认为 banks 中的第一组。

可绘制出临界硼浓度、D棒组棒位和相对功率随燃耗的变化关系:

../_images/plot_result_summary_multi_y_crd.png

或者把D棒组棒位替换成轴向功率偏移AO:

plot-result-summary --y-axis-first boron_concentration --y-axis-second axial_offset --y-axis-third relative_power

可绘制出临界硼浓度、轴向功率偏移和相对功率随燃耗的变化关系:

../_images/plot_result_summary_multi_y_ao.png

11.13. 低功率物理实验结果可视化

11.13.1. 临界实验结果可视化

临界实验结果可视化通过 plot-criticals 命令实现,该功能将 临界实验计算keff1.0 的差值绘制成条形图。

这里我们以 VERA_5 基准题为例来进行临界实验结果可视化。

执行以下交互式命令行:

plot-criticals

可绘制出临界实验结果的条形图:

../_images/plot_criticals.png

11.13.2. 控制棒组价值可视化

控制棒组价值(积分价值)可视化通过 plot-control-rod-worths 命令实现,该功能将 控制棒组价值计算 的结果绘制成条形图。

这里我们以 VERA_5 基准题为例来进行控制棒组价值可视化。

执行以下交互式命令行:

plot-control-rod-worths

可绘制出控制棒组价值的条形图:

../_images/plot_control_rod_worths.png

11.13.3. 控制棒组价值曲线可视化

控制棒组价值曲线(微分价值)可视化通过 plot-control-rod-worths-curve 命令实现,该功能将 控制棒组价值曲线计算 的结果绘制成折线图或条形图。

这里我们以 VERA_5 基准题为例来进行控制棒组价值曲线可视化。

执行以下交互式命令行:

plot-control-rod-worths-curve --rod-bank D --worth-type IRW(pcm)

可绘制出 D 棒组的积分价值的曲线:

../_images/plot_control_rod_worths_curve_IRW.png

--worth-type 选项设置为 DRW(pcm)

plot-control-rod-worths-curve --rod-bank D --worth-type DRW(pcm)

可绘制出 D 棒组的微分价值的条状图:

../_images/plot_control_rod_worths_curve_DRW.png

--worth-type 选项设置为 DRW(pcm/cm)

plot-control-rod-worths-curve --rod-bank D --worth-type DRW(pcm/cm)

可绘制出 D 棒组的微分价值的曲线图,单位为 pcm/cm

../_images/plot_control_rod_worths_curve_DRW_1.png

11.14. 分项耗时统计可视化

分项耗时统计可视化通过 plot-timing-statistics 命令实现,该功能以各个分项(或者模块)的计算耗时来绘制条状图或者饼图。

这里我们以 VERA_7 基准题为例来进行分项耗时统计可视化。

执行以下交互式命令行:

plot-timing-statistics --plot-type bar

可绘制出分项耗时统计的条状图:

../_images/plot_timing_statistics_bar.png

也可以通过 --plot-type 设置为 pie

plot-timing-statistics --plot-type pie

可绘制成饼图,方便查看各个分项耗时的占比:

../_images/plot_timing_statistics_pie.png

可见对于 VERA_7 基准题,计算耗时占比最大的是热工水力反馈,占20%左右,第二是两维MOC模块,占比17%, 第三是几何处理,占比10%。