Considering the inversion of the source term for neutron transport in such as fission reaction
the numerical solution of the effective multiplication factor related to k-eigenvalue problem is investigated. Based on the spherical harmonic functions expansion and finite difference for discretization of transport equation
an inverse power scheme is presented for calculationof the k-eigenvalue problem. This scheme can obviously speed up convergence rate under the proper initial value condition. The numerical results show that the error of the inverse power method gets only 0.045 45% at 3rd iterative step while the error of the standard power method is still reaching to 0.109% after 20 iterative steps
therefore the inverse power algorithm is endowed with higher computing rate and accuracy than the standard power method if an available estimate of the effective multiplication factor is pre-provided.
关键词
Keywords
references
Modak R S, Gupta A. A simple scheme for the direct evaluation of time-eigenvalues of neutron transport equation [J]. Annals of Nuclear Energy, 2003,30(1):211-222.
Quarteroni A, Sacco R, Saleri F. Numerical mathematics [M]. Berlin: Springer Verlag,2002:196-202.
杜书华. 输运问题的计算机模拟 [M]. 长沙: 湖南科学出版社,1989:36-73.
Lewis E E, Miller W F. Computational methods of neutron transport [M]. New York: John Wileyand Sons, 1984.
Mei Liquan. A spectral streamline diffusion finite element coupled method of unsteady transportequation in the field of neutron logging [J].J Appl Math and Mech, 1999, 20(7):739-749.
Allen J, Berry R M. The inverse method for calculation of multiplication factors [J]. Annals ofNuclear Energy, 2002,29(8):929-935.