A time-dependent neutron transport method of characteristics formulation with time derivative propagation.

*(English)*Zbl 1351.82077Summary: A new time-dependent Method of Characteristics (MOC) formulation for nuclear reactor kinetics was developed utilizing angular flux time-derivative propagation. This method avoids the requirement of storing the angular flux at previous points in time to represent a discretized time derivative; instead, an equation for the angular flux time derivative along 1D spatial characteristics is derived and solved concurrently with the 1D transport characteristic equation. This approach allows the angular flux time derivative to be recast principally in terms of the neutron source time derivatives, which are approximated to high-order accuracy using the backward differentiation formula (BDF). This approach, called Source Derivative Propagation (SDP), drastically reduces the memory requirements of time-dependent MOC relative to methods that require storing the angular flux. An SDP method was developed for 2D and 3D applications and implemented in the computer code DeCART in 2D. DeCART was used to model two reactor transient benchmarks: a modified TWIGL problem and a C5G7 transient. The SDP method accurately and efficiently replicated the solution of the conventional time-dependent MOC method using two orders of magnitude less memory.

##### MSC:

82C80 | Numerical methods of time-dependent statistical mechanics (MSC2010) |

65M25 | Numerical aspects of the method of characteristics for initial value and initial-boundary value problems involving PDEs |

82D75 | Nuclear reactor theory; neutron transport |

PDF
BibTeX
XML
Cite

\textit{A. J. Hoffman} and \textit{J. C. Lee}, J. Comput. Phys. 307, 696--714 (2016; Zbl 1351.82077)

Full Text:
DOI

##### References:

[1] | Sutton, T. M.; Aviles, B. N., Diffusion theory methods for spatial kinetics calculations, Prog. Nucl. Energy, 30, 119-182, (1996) |

[2] | Goluoglu, S.; Dodds, H. L., A time-dependent, three-dimensional neutron transport methodology, Nucl. Sci. Eng., 139, 248-261, (2001) |

[3] | Pautz, A.; Birkhofer, A., DORT-TD: a transient neutron transport code with fully implicit time integration, Nucl. Sci. Eng., 145, 299-319, (2003) |

[4] | Cho, J. Y., Transient capability for a MOC-based whole core transport code decart, Trans. Am. Nucl. Soc., 92, 721, (2005) |

[5] | J.B. Taylor, A.J. Baratta, A time-dependent method of characteristics for 3D nuclear reactor kinetics applications, in: Proc. Intl. Conf. on Math. and Comput. Methods Applied to Nucl. Sci. Eng. (M&C2009), Saratoga Springs, USA, May 3-7, 2009. |

[6] | Seubert, A., The time-dependent 3-D transport code TORT-TD: recent advances and applications, Trans. Am. Nucl. Soc., 104, 879-882, (2011) |

[7] | Talamo, A., Numerical solution of the time dependent neutron transport equation by method of characteristics, J. Comput. Phys., 240, 248-267, (2013) |

[8] | K. Tsujita, et al., Higher order treatment on temporal derivative of angular flux for time-dependent MOC, in: Proc. Intl. Conf. on Math. and Comput. Methods Applied to Nucl. Sci. Eng. (M&C2013), Sun Valley, USA, May 5-9, 2013. |

[9] | Askew, J. R., A characteristics formulation of the neutron transport equation in complicated geometries, (1972), United Kingdom Atomic Energy Authority, AAEW-M 1108 |

[10] | Halsall, M. J., Cactus, a characteristic solution to the neutron transport equations in complicated geometries, (1980), United Kingdom Atomic Energy Authority, AEEW-R-1291 |

[11] | H.G. Joo, et al., Methods and performance of a three-dimensional whole-core transport code DeCART, in: Proc. PHYSOR2004, Chicago, USA, 2004. |

[12] | M. Hursin, et al., DeCART v2.05 theory manual, 2008. |

[13] | T.M. Pandya, M.L. Adams, Method of long characteristics applied in space and time, in: Proc. Intl. Conf. on Math. and Comput. Methods Applied to Nucl. Sci. Eng. (M&C2009), Saratoga Springs, USA, May 3-7, 2009. |

[14] | Tsujita, K.; Endo, T.; Yamamoto, A., Kinetic calculation method in space-time frame using characteristic line, Trans. Am. Nucl. Soc., 106, 743-746, (2012) |

[15] | A.J. Hoffman, J.C. Lee, A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics, in: Proc. Intl. Conf. on Math. and Comput. Methods Applied to Nucl. Sci. Eng. (M&C2013), Sun Valley, USA, May 5-9, 2013. |

[16] | Hoffman, A. J., A time-dependent method of characteristics formulation with time derivative propagation, (2013), University of Michigan Department of Nuclear Engineering and Radiological Sciences, Dissertation |

[17] | Yasinsky, J. B.; Henry, A. F., Some numerical experiments concerning space-time reactor kinetics behavior, Nucl. Sci. Eng., 22, 171, (1965) |

[18] | Lewis, E., Benchmark specification for deterministic 2D/3D MOX fuel assembly transport calculations without spatial homogenization, Nuclear Energy Agency, (2001) |

[19] | Ginestar, D., High order backward discretization of the neutron diffusion equation, Ann. Nucl. Energy, 25, 47-64, (1998) |

[20] | Garcia, V. M., Parallel resolution of the two-group time dependent neutron diffusion equation with public domain ODE codes, (Lecture Notes in Computer Science, vol. 3402, (2005)), 457-497 |

[21] | Shim, C. B., Application of backward differentiation formula to spatial reactor kinetics calculation with adaptive time step control, Nucl. Eng. Tech., 43, 531-546, (2011) |

[22] | Hoffman, A. J.; Lee, J. C., Low-order approximations to the angular flux time derivative for transport-based reactor kinetics, Trans. Am. Nucl. Soc., 108, 777-780, (2013) |

[23] | Hamilton, S.; De Oliveira, C.; Clarno, K. T., Error control in a time-dependent slice balance method, Trans. Am. Nucl. Soc., 97, 533-535, (2007) |

[24] | Fornberg, B., Generation of finite difference formulas on arbitrarily spaced grids, Math. Comput., 51, 184, 699-706, (1988) · Zbl 0701.65014 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.