DOI QR코드

DOI QR Code

Delta-form-based method of solving high order spatial discretization schemes for neutron transport

  • Zhou, Xiafeng (Department of Nuclear Engineering and Technology, School of Energy and Power Engineering, Huazhong University of Science and Technology) ;
  • Zhong, Changming (Department of Nuclear Engineering and Technology, School of Energy and Power Engineering, Huazhong University of Science and Technology) ;
  • Li, Fu (Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education, Institute of Nuclear and New Energy Technology, Tsinghua University)
  • Received : 2020.06.08
  • Accepted : 2020.12.29
  • Published : 2021.07.25

Abstract

Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor SN transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional SN codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the SN transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.

Keywords

Acknowledgement

This research is supported by the National Natural Science Foundation of China (12005073), the National Key R&D Program of China (2018YFE0180900), Science and Technology on Reactor System Design Technology Laboratory Foundation Project (HT-KFKT-10-2018005), Key Laboratory Foundation Project of Advanced Reactor Engineering and Safety of Ministry of Education (ARES-2018-07).

References

  1. Y. Azmy, E. Sartori, Nuclear Computational Science: a Century in Review [M], Springer Science & Business Media, 2010.
  2. K.D. Lathrop, F.W. Brinkley, TWOTRAN-II: an Interfaced, Exportable Version of the TWOTRAN Code for Two-Dimensional Transport [R], Los Alamos Scientific Laboratory, New. Mexico, USA, 1973.
  3. M.L. Adams, E.W. Larsen, Fast iterative methods for discrete-ordinates particle transport calculations, Prog. Nucl. Energy 40 (1) (2002) 3-159. https://doi.org/10.1016/S0149-1970(01)00023-3
  4. B. Collins, T. Downar, J. Gehin, et al., MPACT Theory Manual [R]. Technical Report CASL-U-2015-0078-000, Consortium for Advanced Simulation of LWRs, 2015.
  5. A. Prosperetti, G. Tryggvason, Computational Methods for Multiphase Flow [M], Cambridge university press, 2009.
  6. C.W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws [M], in: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations, Springer, Berlin, Heidelberg, 1998, pp. 325-432.
  7. R. Wang, H. Feng, R.J. Spiteri, Observations on the fifth-order WENO method with non-uniform meshes [J], Appl. Math. Comput. 196 (1) (2008) 433-447. https://doi.org/10.1016/j.amc.2007.06.024
  8. W.F. Godoy, P.E. DesJardin, On the use of flux limiters in the discrete ordinates method for 3D radiation calculations in absorbing and scattering media, J. Comput. Phys. 229 (2010) 3189-3213. https://doi.org/10.1016/j.jcp.2009.12.037
  9. C.C. Pain, M.D. Eaton, J. Bowsher, et al., Finite element based Riemann solvers for time-dependent and steady-state radiation transport, Transport Theor. Stat. Phys. 32 (5-7) (2003) 699-712. https://doi.org/10.1081/TT-120025073
  10. T.A. Brunner, J.P. Hollowa, Two-dimensional time dependent Riemann solvers for neutron transport [J], J. Comput. Phys. 210 (1) (2005) 386-399. https://doi.org/10.1016/j.jcp.2005.04.011
  11. D. Wang, T. Byambaakhuu, High-order Lax-Friedrichs WENO fast sweeping methods for the SN Neutron transport equation, Nucl. Sci. Eng. 193 (9) (2019) 982-990. https://doi.org/10.1080/00295639.2019.1582316
  12. Y.J. Zhu, C.W. Zhong, K. Xu, An implicit unified gas-kinetic scheme for unsteady flow in all Knudsen regimes, J. Comput. Phys. 386 (2019) 190-217. https://doi.org/10.1016/j.jcp.2019.01.033
  13. D.X. Pan, C.W. Zhong, C.S. Zhuo, An implicit discrete unified gas-kinetic scheme for simulations of steady flow in all flow regimes, Commun. Comput. Phys. 25 (5) (2019) 1469-1495.
  14. C. Zhang, Z.L. Guo, S.Z. Chen, Unified implicit kinetic scheme for steady multiscale heat transfer based on the phonon Boltzmann transport equation, Phys. Rev. E 96 (2017), 063311. https://doi.org/10.1103/PhysRevE.96.063311
  15. J. Smit, M. van Sint Annaland, J.A.M. Kuipers, Grid adaptation with WENO schemes for non-uniform grids to solve convection-dominated partial differential equations, Chem. Eng. Sci. 60 (10) (2005) 2609-2619. https://doi.org/10.1016/j.ces.2004.12.017
  16. J.B. Taylor, The Development of a Three-Dimensional Nuclear Reactor Kinetics Methodology Based on the Method of Characteristics, the Pennsylvania State University, 2007. PhD Thesis.
  17. X.F. Zhou, C.M. Zhong, J. Hou, et al., A unified space-time fully-transversed general nodal expansion method for multidimensional and multigroup steady and transient nuclear reactor neutronics problems [J], Ann. Nucl. Energy 132 (2019) 427-441. https://doi.org/10.1016/j.anucene.2019.04.058
  18. K. Kobayashi, N. Sugimura, 3D radiation transport benchmark problems and results for simple geometries with void region [J], Prog. Nucl. Energy 39 (2) (2001) 119-144. https://doi.org/10.1016/S0149-1970(01)00007-5