Simulation for non-homogeneous transport equation by Nyström method
DOI:
https://doi.org/10.15392/bjrs.v8i3A.1507Keywords:
transport equation, integral formulation, Nyström methodAbstract
In this work we solve numerically the one-dimensional transport equation with semi-reflective boundary conditions and non-homogeneous domain. The proposed methodology consists of applying the Nyström method in order to discretize the integral formulation of this problem which is an equation involving weakly singular integral operators. For this purpose, analytical and computational techniques were applied to deal with the singularities. The Nyström method is an integral method which approximates the integral operator by a numerical quadrature and turns the integral equation into a finite dimensional linear system. This formulation allows us to use any function to describe both scattering cross section and total cross section. The algorithm is implemented in C language with the use of routines of GNU scientific library and computational techniques for code optimization. The scalar flux was calculated for two numerical quadrature, namely Gauss-Legendre quadrature and Boole's rule. The numerical results were determined for transport problem with homogeneous and non-homogeneous domains. In order to validate the proposed method-ology, our numerical results were compared with those from the literature and presented with several correct significant digits.Downloads
References
CHANDRASEKHAR, S. Radiative transfer, Dover Publications, Inc., Oxford University Press, London, 1950.
LEWIS, E. E.; MILLER, J. W. F. Computational Methods of Neutron Transport, Copyright, United States of America, 1984.
DALMOLIN, D.; DE AZEVEDO F. S.; SAUTER, E. Nyström method in transport equation, PROCEEDING OF INAC 2017 INTERNATIONAL NUCLEAR ATLANTIC CONFERENCE, Rio de Janeiro, Brazil, 2017.
DE AZEVEDO, F. S.; SAUTER, E.; KONZEN, P. H. A.; THOMPSON, M.; BARICHELLO, L. B. Integral formulation and numerical simulations for the neutron transport equation in X-Y geometry. Annals of Nuclear Energy, vol. 112, p. 735-747, 2018.
SAUTER, E.; DE AZEVEDO, F. S.; KONZEN, P. H. A. Nyström Method Applied to the Transport Equation in a Semi-Reflective Rectangle. Journal of Computational and Theoretical Transport, p. 2332-4325, 2019.
DELVES, L. M.; MOHAMED, J. L. Computational methods for integral equations, Cambridge, UK: Cambridge University Press, 1985.
NUNES, C. E. A.; BARROS, R. C. Aplicativo computacional para cálculos de blindagem com modelo de transporte Sn unidimensional e monoenergético, PROCEEDING OF INAC 2009 INTERNATIONAL NUCLEAR ATLANTIC CONFERENCE, Rio de Janeiro, Brazil, 2009.
GARCIA, R. D. M.; SIEWERT, C. E. Radiative transfer in finite inhomogeneous plane-parallel atmospheres. Journal of Quantitative Spectroscopy e Radiative Transfer, 2009.
Free Software Foundation, INC., 2009. GSL-GNU Scientific Library.
Published
Issue
Section
License
Copyright (c) 2021 Brazilian Journal of Radiation Sciences
This work is licensed under a Creative Commons Attribution 4.0 International License.
Licensing: The BJRS articles are licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
How to Cite
