Cinética espacial na teoria da difusão de nêutrons: uma solução via abordagens noclais e exponenciais matriciais

Rodrigo Zanette, Liliane B. Barichello, Cláudio Z. Petersen

Resumo


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é.


Palavras-chave


Cinética Espacial, Teoria da Difusão de Nêutrons, Exponencial Matricial.

Texto completo:

PDF

Referências


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.




DOI: https://doi.org/10.5540/03.2021.008.01.0382

Apontamentos

  • Não há apontamentos.


SBMAC - Sociedade de Matemática Aplicada e Computacional
Edifício Medical Center - Rua Maestro João Seppe, nº. 900, 16º. andar - Sala 163 | São Carlos/SP - CEP: 13561-120
 


Normas para publicação | Contato