Abstract
A discrete ordinates method for three-dimensional neutron transport equation based on unstructured-meshes was derived from the first-order neutron transport equation. A set of differential equations about the spatial variables were obtained. It avoids the singularity in void regions from the second-order form of the equation, which implies the inversion of the total cross-section. The differential equations were solved iteratively using the least-squares finite element method. For the symmetric stiffness matrix, the fast iteration method can be used. A three-dimensional transport calculation code was programmed. The triangular prism and the tetrahedron elements were used in the calculation. The numerical results of some benchmark problems show that this method can be used in unstructured-meshes and can give high precise result. For most problems, the error is less than 0.3% for the eigenvalue and 3.0% for the angular flux respectively.
| Original language | English |
|---|---|
| Pages (from-to) | 363-366 |
| Number of pages | 4 |
| Journal | Hsi-An Chiao Tung Ta Hsueh/Journal of Xi'an Jiaotong University |
| Volume | 41 |
| Issue number | 3 |
| State | Published - Mar 2007 |
Keywords
- Discrete ordinate method
- Finite element
- Least-square
- Unstructured-mesh