Improvement of Geant4 Neutron-HP package: Doppler broadening of the neutron elastic scattering kernel and cross sections

M. Zmeškal,L. Thulliez,E. Dumonteil
2023-03-14
Abstract:Whether it is for shielding applications or for safety criticality studies, numerically solving the neutron transport equation with a good accuracy requires to precisely estimate the Doppler broadened elastic scattering kernel in the thermal and epithermal energy range of neutrons travelling in a free gas. In Geant4, low energy neutrons are transported using evaluated data libraries handled by the Neutron High-Precision (Neutron-HP) package. Version 11.00.p03 of the code features in particular the Doppler broadened elastic scattering kernel, provided by the so-called 'Sampling of the Velocity of the Target' (SVT) method. However this latter fails for resonant heavy nuclei such as 238U and can severely impact the solving of the Boltzmann equation in fissile media. To overcome this shortcoming, the Doppler Broadened Rejection Correction (DBRC) method has been implemented in Geant4 and successfully validated with the reference Monte Carlo neutron transport code Tripoli4 (version 11). This development will be taken into account in the next release of the code. The cross section Doppler broadening process, which is performed on-the-fly, is also carefully investigated and ways to improve it on a simulation-by-simulation basis are presented. All the validations have been performed with an automated benchmark tool which has been designed to support the quality assurance of the Geant4 Neutron-HP package. This tool is currently available on an ad hoc Gitlab repository and will be included in Geant4.
Computational Physics,Nuclear Experiment
What problem does this paper attempt to address?
This paper attempts to solve the problem of inaccurate calculation of Doppler - broadened elastic scattering cross - sections due to the resonance effect of heavy nuclei in low - energy neutron transport simulations. Specifically, the paper focuses on how to improve the Neutron - HP package in Geant4 to more accurately handle the Doppler - broadening of neutron elastic scattering cross - sections in the thermal and epithermal energy ranges. ### Background and Problem Description In Geant4, the transport of low - energy neutrons is handled using evaluated databases, which are managed by the Neutron High - Precision (Neutron - HP) package. Geant4 version 11.00.p03 introduced Doppler - broadened elastic scattering cross - sections, which were implemented through the so - called "Target - Nucleus Velocity Sampling" (SVT) method. However, the SVT method has flaws when dealing with resonant heavy nuclei (such as \(^{238}\text{U}\)), which can affect the accuracy of solving the Boltzmann equation in fission media. To overcome this shortcoming, the Doppler - Broadening Rejection Correction (DBRC) method was introduced in the paper and successfully verified in the reference Monte Carlo neutron transport code Tripoli - 4® (version 11). ### Solutions 1. **Implementation of the DBRC method**: - The paper details the implementation process of the DBRC method. This method precisely samples the \((\mu, V_t)\) pairs by adding a rejection loop/criterion on top of the SVT algorithm, thus obtaining the exact Doppler - broadened elastic scattering nucleus. - The specific formula is as follows: \[ p_{\text{DBRC}}(V_t, \mu) = C' p_{\text{SVT}}(V_t, \mu) \left(\frac{\sigma_0(v_r)}{\sigma_{\text{max}}^0(v_\xi)}\right) \] where \(C'\) is a normalization constant, and \(\sigma_{\text{max}}^0(v_\xi)\) is the maximum 0 K cross - section within the range of \(v_n \pm 4\sqrt{\frac{2k_BT}{M}}\). 2. **Verification of the DBRC method**: - The accuracy of the DBRC method was verified by comparing the results of Geant4 and Tripoli - 4® through a thin - cylinder benchmark test. - The results show that the secondary neutron energy distribution of the DBRC method at different temperatures is significantly different from that of the SVT method, and is consistent with the results of Tripoli - 4® within the statistical error range. 3. **Real - time cross - section Doppler - broadening**: - The paper also studied the process of real - time cross - section Doppler - broadening in Geant4 and proposed an improvement method. - By adjusting the convergence criterion \(\epsilon_{xs}\), the accuracy of cross - section calculation can be improved, but the calculation time will be increased. The paper suggests choosing an appropriate value of \(\epsilon_{xs}\) in practical applications to balance accuracy and computational efficiency. ### Conclusion By introducing and verifying the DBRC method, the paper significantly improves the calculation accuracy of Doppler - broadened elastic scattering cross - sections in Geant4 when dealing with resonant heavy nuclei. In addition, the paper also explores the optimization strategy for real - time cross - section Doppler - broadening, providing an important reference for future neutron transport simulations.