Abstract
In this paper, we introduce a novel approach to numerical solving for nuclear reactivity using the inverse equation of point kinetics as a model. Our method leverages the Gauss-Lobatto quadrature, incorporating multiple time steps and nodes derived from Legendre polynomials of varying orders. To enhance the accuracy of the differential component of the inverse equation, we employed a differential scheme based on the Gauss-Lobatto quadrature nodes. To address fluctuations in neutron population density, we implemented a second-order low-pass Butterworth filter with a minimal window length of M = 3 samples. Through extensive numerical simulations, we assessed the precision of our proposed method and filtering process by varying time steps and standard deviations associated with noise or uncertainty. Additionally, we benchmarked our results against the Savitzky-Golay filter, which uses a significantly larger sampling window of M = 225. Our findings reveal that the integration of the Gauss-Lobatto quadrature method with the Butterworth filter not only significantly reduces fluctuations but also demonstrates potential for effective implementation in digital reactivity meters.
Recommended Citation
Suescún-Díaz, Daniel; Cardoso-Páez, Luis E.; and Peña-Lara, Diego
(2024)
"Gauss-Lobatto Method for Nuclear Reactivity Calculation,"
Karbala International Journal of Modern Science: Vol. 11
:
Iss.
1
, Article 2.
Available at:
https://doi.org/10.33640/2405-609X.3384
Creative Commons License
This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License.
