This program is based on the Method of Characteristic(MOC), and the idea is from an experience when learning the numerical calculation and undertaking a relative SRT in Tsinghua University. While the program only fulfill the basic functions now and maybe in progress updating someday, I provided a simply easy-taking way to try building the program systematically and hand-by-hand, which benefits the understanding of the mainstream way of neutron transport equation calculating and the precautions when programming and directions for future efforts,especially for the freshmen in this field.
the report is provided in Chinese.It provided the mathmatic background, the algorithm and module functions in this report.
在终端运行
python3 main.py
会自动调用src中的模块,生成名为results的文件夹,并在终端输出目前进程和$k_{eff}$的变化,推荐记录和保存终端输出,因为results文件夹中并未直接保存$k_{eff}$。
(1)3*3正方形生成0.5大小的Delunay三角形网格,可修改。
(2)2D平面中的中子输运计算,但可以拓展到三维情况下在z向在几何和材料上均是各项同性的情况。故强制只使用了$\varphi$为0,30,60,90,120,150,180,210,240,270,300,330(度)共12个均匀分布的方位角,不可修改。
(3)该程序只设置了真空边界条件去进行计算。
(4)该程序方便起见,默认选择只精确到$k_{eff}$小数点后三位。
(1)推荐如果验证多群计算,使用GEO中默认参数,但此时在6群计算中会出现psi_avg_mesh<0的情况,这是由于网格密度太低引起的,但仍会正确收敛。
(2)如果验证网格密度较大的情况时,如mesh_size = 0.2的情况,不推荐使用params_jsons文件夹中的params4.json。因为6群的迭代收敛速率会很慢,尽管此时并不会出现(1)中所述情况。mesh_size = 0.2的情况的计算结果在examples文件夹中保存。
(3)运行main程序时,程序会自动检查results文件夹并进行删除文件夹重建。
(4) INTERSECTION_2.pdf和INTERSECTION.pdf给出两个不同角度插值的示意图。