DOI QR코드

DOI QR Code

Incorporation of anisotropic scattering into the method of characteristics

  • Rahman, Anisur (Department of Nuclear Engineering, Ulsan National Institute of Science and Technology) ;
  • Lee, Deokjung (Department of Nuclear Engineering, Ulsan National Institute of Science and Technology)
  • Received : 2021.08.30
  • Accepted : 2022.03.29
  • Published : 2022.09.25

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.

Keywords

Acknowledgement

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT). (No. NRF-2019M2D2A1A03058371).

References

  1. S. Choi, D. Lee, "Three-dimensional method of characteristics/diamonddifference transport analysis method in STREAM for whole-core neutron transport calculation". Comput. Phys. Commun., 260, 107332. https://doi.org/10.1016/j.cpc.2020.107332.
  2. D. Knott, A. Yamamoto, Lattice physics computations, Handbook of Nuclear Engineering (2010) 913-1239.
  3. X-5 Monte Carlo Team, "MCNP - A General N-Particle Transport Code, Version 5" Volume I: Overview and Theory, U. R.-03-1987, 2003. updated 2005.
  4. H. Lee, W. Kim, P. Zhang, M. Lemaire, A. Khassenov, J. Yu, Y. Jo, J. Park, D. Lee, "MCS - a Monte Carlo particle transport code for large-scale power reactor analysis," Ann. Nucl. Energy, 139, 107276. https://doi.org/10.1016/j.anucene.2019.107276.
  5. S. Choi, K. Smith, H.C. Lee, D. Lee, Impact of inflow transport approximation on light water reactor analysis, Comp. Physiol. 299 (2015) 352-373, https://doi.org/10.1016/j.jcp.2015.07.005.
  6. M.J. Halsall, CACTUS, A Characteristic Solution to the Neutron Transport Equation in Complicated Geometries, AEEW-R, UK Atomic Energy Authority Establishment, Winfrith, 1980.
  7. H.G. Joo, J.Y. Cho, K.-S. Kim, C.C. Lee, S.Q. Zee, Methods and performance of a three-dimensional whole-core transport code DeCART, in: Proceedings of PHYSOR 2004, 2004. Chicago, IL.
  8. Y.S. Jung, C.B. Shim, C.H. Lim, H.G. Joo, Practical numerical reactor employing direct whole core neutron transport and subchannel thermal/hydraulic solvers, Ann. Nucl. Energy 62 (2015) 357-374, https://doi.org/10.1016/j.anucene.2013.06.031.
  9. B. Collins, S. Stimpson, B. Kochunas, T. Downar, Assessment of the 2D/1D implementation in MPACT, in: Proceedings of PHYSOR 2014, Sep. 29, 2014. Kyoto, Japan.
  10. A. Marin-Lafleche, M.A. Smith, C.H. Lee, PROTEUS-MOC: a 3D deterministic solver incorporating 2D method of characteristics, in: Proceedings of M&C 2013, 2013. Sun Valley.
  11. W. Boyd, S. Shaner, L. Li, B. Forget, K. Smith, The OpenMOC method of characteristics neutral particle transport code, Ann. Nucl. Energy 68 (2014) 43-52, https://doi.org/10.1016/j.anucene.2013.12.012.
  12. M. Jarrett, B. Kochunas, E. Larsen, T. Downar, Progress in characterizing 2D/1D accuracy in MPACT, in: Proceedings of M&C 2017, 2017. Jeju, South Korea.
  13. A.M. Graham, B.S. Collins, T. Downar, Improvement of the 2D/1D method in MPACT using the subplane scheme, in: Proceedings of M&C 2017, 2017. Jeju, South Korea.
  14. Z. Liu, H. Wu, L. Cao, Q. Chen, Y. Li, A new three-dimensional method of characteristics for the neutron transport calculation, Ann. Nucl. Energy 38 (2011) 447-454, https://doi.org/10.1016/j.anucene.2010.09.021.
  15. M.A. Shafii, Solution methods of neutron transport equation in nuclear reactors, J. I. Dasar 14 (2013) 59-65, https://doi.org/10.19184/jid.v14i2.320.
  16. Liangzhi Cao, Hongchun Wu, Spherical harmonics method for neutron transport equation based on unstructured-meshes, Nucl. Sci. Tech. 15 (2004).
  17. K. Kobayashi, H. Oigawa, H. Yamagata, The spherical harmonics method for the multigroup transport equation in x-y geometry, Ann. Nucl. Energy 13 (1986) 663-678, https://doi.org/10.1016/0306-4549(86)90045-9.
  18. A. Rahman, A. Cherezov, D. Lee, Higher Order PN Solver in Method of Characteristics Code STREAM, Transactions of the Korean Nuclear Society Autumn Meeting Changwon, Korea, 2020.
  19. G. Palmiotti, E. E. Lewis, C. B. Carrico, VARIANT: VARIational Anisotropic Nodal Transport for Multidimensional Cartesian and Hexagonal Geometry Calculation, ANL-95/40, Argonne National Laboratory, USA.
  20. A.T. Godfrey, Vera Core Physics Benchmark Progression Problem Specifications, Fourth Revision, CASL, 2014. CASL-U-2012-0131-004.
  21. M.N. Baldwin, G.S. Hoovler, R.L. Eng, F.G. Welfare, Critical Experiments Supporting Close Proximity Water Storage of Power Reactor Fuel, The Babcock & Wilcox Company, 1979.