适用于任意几何的特征线边界条件处理方法
Boundary Condition Processing Method for MOC Calculation in Arbitrary Geometry
-
摘要: 边界条件处理是特征线方法(MOC)向任意三维几何拓展时遇到的难点之一。本文提出一种边界条件处理方法,既保留循环特征线中首尾相连的特性,又能像插值方法一样适用于任意几何。首先推导了平源近似下的特征线方程,提出了一种将源项和边界角通量分离处理的内迭代解法。然后证明了该解法具有唯一解,并类似于循环特征线方法给出解的构造方法。最后借助数值积分和权重插值给出迭代计算流程。采用Takeda算例、单铀球水腔模型和C5G7算例进行验证计算,keff的最大计算误差分别为21、319和138.8 pcm,表明方法可靠。该方法可应用于任意几何,且不需存储边界通量和进行边界迭代。Abstract: Boundary condition processing is one of the difficulties encountered in the application of method of characteristics (MOC) to arbitrary three-dimensional geometry. In this paper, a boundary condition processing method was proposed, which not only preserved the track continuity as cyclic track method, but also could be applied to arbitrary geometry as the interpolation method. The MOC equation was derived under the flat source approximation and an internal iterative method was proposed in which the source term and the boundary angular flux were processed separately. It was proved that the equation had a unique solution which could be constructed similarly to the cyclic track method. The iterative calculation flow was given by numerical integration and weight interpolation. Takeda benchmark, single uranium sphere model with water cavity and C5G7 benchmark were calculated to test the accuracy. The maximum error of keff is 21, 319 and 138.8 pcm respectively, which shows that the method is reliable. This method can be applied to arbitrary geometry without storing boundary fluxes and performing boundary iteration.
-
-
[1] SANCHEZ R. Prospects in deterministic three-dimensional whole core transport calculation[J]. Nuclear Engineering and Technology, 2012, 44(2): 113-150. [2] KOCHUNAS B M. A hybrid parallel algorithm for the 3-D method of characteristics solution of the boltzmann transport equation on high performance compute clusters[D]. Michigan: University of Michigan, 2013. [3] SHANER S, GUNOW G, FORGET B, et al. Verification of the 3D method of characteristics solver in OpenMOC[C]∥Proceedings of Physics of Reactors Conference 2016 (Physor 2016). USA: American Nuclear Society, 2016. [4] SCIANNANDRONE D, SANTANDREA S, SANCHEZ R. Optimized tracking strategies for step MOC calculations in extruded 3D axial geometries[J]. Annals of Nuclear Energy, 2016, 87: 49-60. [5] 汤春桃. 中子输运方程特征线解法及嵌入式组件均匀化方法的研究[D]. 上海:上海交通大学机械与动力工程学院,2009. [6] 柴晓明,姚栋,王侃. 基于特征线方法的三维中子输运程序,Ⅰ:边界条件的插值处理[J]. 核动力工程,2010,31(2):11-15.CHAI Xiaoming, YAO Dong, WANG Kan. Three dimension neutron transport characteristics method code, Ⅰ: Interpolation treatment of boundary condition[J]. Nuclear Power Engineering, 2010, 31(2): 11-15(in Chinese). [7] 吴宏春,郑友琦,曹良志,等. 中子输运方程确定论数值解法[M]. 北京:中国原子能出版社,2018. [8] TAKEDA T, IKEDA H. Final report on the 3-D neutron transport benchmarks[R]. [S. l.]: Organisation for Economic Co-operation and Development Nuclear Energy Agency, 1991. [9] KELLER S E, OLIVEIRA D, CASSIANO R E. Two-dimensional C5G7 mox fuel assembly benchmark calculations using the FEM-PN code EVENT[J]. Progress in Nuclear Energy, 2004, 45(2-4): 255-263. [10] 马纪敏,刘永康,李茂生. 次临界能源堆用多群截面库的研制与校验[J]. 核动力工程,2012,33(5):16-21.MA Jimin, LIU Yongkang, LI Maosheng. Development and validation of multi-group cross-section library for subcritical energy reactor[J]. Nuclear Power Engineering, 2012, 33(5): 16-21(in Chinese).
计量
- 文章访问数: 253
- HTML全文浏览量: 7
- PDF下载量: 870