A least-squares finite-element S-n method for solving first-order neutron transport equation

作者：Ju Hai Tao^*; Wu Hong Chun; Zhou Yong Qiang; Cao Liang Zhi; Yao Dong; Xian Chun Yu

来源：Nuclear Engineering and Design, 2007, 237(8): 823-829.

DOI：10.1016/j.nucengdes.2006.10.024

A discrete ordinates finite-element method for solving the two-dimensional first-order neutron transport equation is derived using the least-squares variation. It avoids the singularity in void regions of the method derived from the second-order equation which contains the inversion of the cross-section. Different from using the standard Galerkin variation to the first-order equation, the least-squares variation results in a symmetric matrix, which can be solved easily and effectively. To eliminate the discontinuity of the angular flux on the vacuum boundary in the spherical harmonics method, the angle variable is discretized by the discrete ordinates method. A two-dimensional transport simulation code is developed and applied to some benchmark problems with unstructured geometry. The numerical results verified the validity of this method.

• 出版日期2007-4
• 单位中国核动力研究设计院; 西安交通大学