Cinética espacial na teoria da difusão de nêutrons: uma solução via abordagens noclais e exponenciais matriciais
DOI:
https://doi.org/10.5540/03.2021.008.01.0382Palabras clave:
Cinética Espacial, Teoria da Difusão de Nêutrons, Exponencial Matricial.Resumen
Neste trabalho, uma solução para o problema da cinética espacial na teoria da difusão de nêutrons multigrupo em geometria cartesiana unidimensional é derivada. A partir da divisão do domínio espacial em nodos homogêneos, uma integração nodal é aplicada nas equações da cinética espacial, obtendo variáveis médias em cada nodo: fluxos, concentrações de precursores e densidades de correntes. As densidades de correntes para cada interface são aproximadas pela média ponderada dos coeficientes de difusão e pelos fluxos médios. Uma solução analítica na variável temporal na forma de uma exponencial matricial é proposta e, essa exponencial matricial, é avaliada pela aproximação de Padé e pelo método de Schur-Parlett. Os resultados numéricos obtidos são comparáveis aos resultados existentes na literatura e o método de Schur-Parlett mostra uma maior eficiência à aproximação de Padé.
Descargas
Citas
Aboanber, A. E. and Nahla, A.A. Solution of two-dimensional space-time multigroup reactor kinetics equations by generalized Pad´e and cut-product approximations, Annals of Nuclear Energy, 33: 209–222, 2006. DOI:10.1016/j.anucene.2005.11.003.
Al-Mohy, A. and Higham, N. A New Scaling and Squaring Algorithm for the Matrix Exponential, SIAM Journal on Matrix Analysis and Applications, 31: 970–989, 2009. DOI:10.1137/09074721X.
Anderson, E., Bai, Z., Bischof, C., Blackford, S., Demmel, J., Dongarra, J., Du Croz, J., Greenbaum, A., Hammarling, S., McKenney, A. and Sorensen, D. LAPACK Users’ guide. SIAM, Philadelphia, 1999.
Ceolin, C., Schramm M., Vilhena, M. T. and Bodmann, B. E. J. On the Neutron multi-group kinetic diffusion equation in a heterogeneous slab: An exact solution on a finite set of discrete points, Annals of Nuclear Energy, 76:271–282, 2015. DOI:10.1016/j.anucene.2014.09.038.
Davies, P. and Higham, N. A Schur-Parlett Algorithm for Computing Matrix Functions, SIAM Journal on Matrix Analysis and Applications, 25:464–485, 2003. DOI:10.1137/S0895479802410815.
Duderstadt, J. J. and Hamilton, L. J. Nuclear Reactor Analysis. John Wiley, New York, 1976. [7] Evard, J. C. On matrix functions which commute with their derivative, Linear Algebra and its Applications, 68:145–178, 1985. DOI:10.1016/0024-3795(85)90212-5.
Higham, N. The Scaling and Squaring Method for the Matrix Exponential Revisited, SIAM Journal on Matrix Analysis and Applications, 26:1179–1193, 2005. DOI:10.1137/04061101X.
Li, M., Chen, W., and Hao, J. A quasi-static orthogonal method for solving spatial-time neutron kinetics equations in non-steady processes, Progress in Nuclear Energy, 120:103229, 2020. DOI:10.1016/j.pnucene.2019.103229.
Moler, C. and Van Loan, C. Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later, SIAM Review, 45:3–49, 2003. DOI:10.1137/S0036144502418.
Parlett, B. N. A recurrence among the elements of functions of triangular matrices, Linear Algebra and its Applications, 14:117–121, 1976. DOI:10.1016/0024-3795(76)90018-5.
Pollard, J. P. AUS Diffusion Module POW Checkout - 1- and 2- Dimensional Kinetics Cal- culations. Australian Atomic Energy Commission AAEC/387, 1977.
Zanette, R., Barichello, L. B. e Petersen, C. Z. C´alculo de criticalidade pela teoria de difusão de nêutrons: uma análise comparativa de aproximação da densidade de corrente, REMAT: Revista Eletrônica da Matemática, 6:e4006, 2020. DOI: 10.35819/remat2020v6i2id4248.
