Implementation and Comparison of Different Algorithms for on the fly Doppler Broadening in RMC Code
- Department of Engineering Physics, Tsinghua University, Beijing, 100084, (China)
Accurate cross sections are of great importance to Monte Carlo particle transport calculation codes which use pointwise cross sections. It is well known that cross sections have much dependence on temperature on account of Doppler effect which is particularly vital for resonance absorption cross sections. In general, pre-prepared continuous-energy libraries are used in Monte Carlo codes, which is feasible in many applications. However, in multi-physics calculations where neutronics and thermal hydraulics are coupled, cross sections at a fine interval of temperatures are desired, which brings a great challenge to computer memory. Consequently, different methods about on the fly Doppler broadening for cross sections are proposed, which can generate cross sections at the specific temperature points during neutron simulations. These methods can be roughly classified into two categories: one is to calculate the required cross sections by different algorithms, such as interpolation using a high order function expansion, multipole representation, and different integration algorithms based on Doppler broadening equation. The other, for instance, target motion sampling (TMS), does not calculate the Doppler broadening cross sections but takes into account Doppler effect by rejection technique. In this paper, several on the fly Doppler broadening algorithms based on Doppler broadening equation are implemented in the Reactor Monte Carlo code RMC, and then the accuracy and efficiency are discussed and compared. In this paper, several algorithms for on the fly Doppler broadening based on Doppler broadening equation have been implemented in RMC and the focus falls on the comparison of the Gauss Legendre quadrature and Gauss Hermite quadrature. Gauss Legendre quadrature is more accurate than Gauss Hermite quadrature but also more time consuming. Base temperature has a significant impact on Gauss Hermite quadrature while increasing base temperature to room temperature (300 K) can give satisfactory results. In addition, the proportion of overhead of Doppler broadening treatment in total computational time tends to decrease with an increase number of nuclides in fuel. For depletion calculation, both algorithms are in good agreement with reference and the relative errors of K{sub inf} are within 5 times of standard deviation. As for total computational time, the ratio between Gauss Legendre quadrature or Gauss Hermite quadrature and reference decreases and reaches 9.6 or 3.4 at last, respectively. Four important isotopes are chosen to compare the atom density and all relative errors are within 0.25%. After comparisons of different algorithms and different base temperatures, it can be found that Gauss Hermite quadrature broadening from 300 K can perform on the fly cross section Doppler broadening with least total computational time and acceptable precision. (authors)
- OSTI ID:
- 22992004
- Journal Information:
- Transactions of the American Nuclear Society, Journal Name: Transactions of the American Nuclear Society Journal Issue: 1 Vol. 114; ISSN 0003-018X
- Country of Publication:
- United States
- Language:
- English
Similar Records
On-the-fly Doppler Broadening with Probability Table Interpolation for Unresolved Resonance Region in RMC
Neutronics/Thermal-hydraulics coupling with RMC and CTF for BEAVRS benchmark calculation