Incorporation of anisotropic scattering into the method of characteristics
0
조회
0
다운로드
1
인용
0
참고문헌
키워드
초록
In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses
Abstract
In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses
참고문헌 (21)
- Sooyoung Choi, "Three-dimensional method of characteristics", diamond-difference transport analysis method in STREAM for whole-core neutron transport calculation, 2021, pp.260, 107332~[DOI]
- D. Knott, A. Yamamoto, "Lattice physics computations, Handbook of Nuclear Engineering", 913~1239 (2010)
- X-5 Monte Carlo Team, "MCNP-A General N-Particle Transport Code, Version 5", Volume I : Overview and Theory, U. R. -03-1987, ~ (2003)
- Hyunsuk Lee, "MCS – A Monte Carlo particle transport code for large-scale power reactor analysis", Annals of Nuclear Energy, 139, pp.107276~ (2020), Elsevier BV[DOI]
- S. Choi, "Impact of inflow transport approximation on light water reactor analysis", Comp. Physiol, 299, pp.352~373 (2015)[DOI]
- CACTUS, A Characteristic Solution to the Neutron Transport Equation in Complicated Geometries, "M. J. Halsall", UK Atomic Energy Authority Establishment, ~ (1980)
- H.G. Joo, "Methods and performance of a three-dimensional whole-core transport code DeCART", Proceedings of PHYSOR 2004 (2004)
- Y. S. Jung, "Practical numerical reactor employing direct whole core neutron transport and subchannel thermal", hydraulic solvers, 2015, pp.62, 357~374
- B. Collins, "Assessment of the 2D", 1D implementation in MPACT, 2014, ~
- A. Marin-Lafleche, "PROTEUS-MOC: a 3D deterministic solver incorporating 2D method of characteristics", Proceedings of M&C 2013 (2013)
- W. Boyd, "The OpenMOC method of characteristics neutral particle transport code", Ann. Nucl. Energy, 68, pp.43~52 (2014)[DOI]
- M. Jarrett, "Progress in characterizing 2D", 1D accuracy in MPACT, 2017, ~
- A.M. Graham, "Improvement of the 2D", 1D method in MPACT using the subplane scheme, 2017, ~
- Z. Liu, "A new three-dimensional method of characteristics for the neutron transport calculation", Ann. Nucl. Energy, 38, pp.447~454 (2011)[DOI]
- M. A. Shafii, "Solution methods of neutron transport equation in nuclear reactors", J. I. Dasar, 14, pp.59~65 (2013)[DOI]
- Liangzhi Cao, "Spherical harmonics method for neutron transport equation based on unstructured-meshes", Nucl. Sci. Tech, 15 (2004)
- K. Kobayashi, "The spherical harmonics method for the multigroup transport equation in x-y geometry", Ann. Nucl. Energy, 13, pp.663~678 (1986)[DOI]
- Higher Order PN Solver in Method of Characteristics Code STREAM, "A. Rahman", Transactions of the Korean Nuclear Society Autumn Meeting Changwon, ~ (2020)
- VARIANT: VARIational Anisotropic Nodal Transport for Multidimensional Cartesian and Hexagonal Geometry Calculation, ANL-95, "40", G. Palmiotti
- Vera Core Physics Benchmark Progression Problem Specifications, "A.T. Godfrey", CASL, ~, pp.CASL-U-2012-0131-004 (2014)
- M. N. Baldwin, "Critical Experiments Supporting Close Proximity Water Storage of Power Reactor Fuel" (1979), The Babcock &Wilcox Company