An application of a 3D Boundary Element Method (BEM), coupled with the Response Matrix (RM) technique, to solve the neutron diffusion and transport equations for multi-region domains is presented. The discussion is here limited to steady state problems, in which the neutrons have a wide energy spectrum, which leads to systems of several diffusion or transport equations. Moreover, the number of regions with different physical constants can be very large. The boundary integral equations concerning each region are solved via a polynomial moment expansion and, taking advantage of suitable recurrence formulas, the multi-fold integrals there involved are reduced to single or double integrals. The usual unknowns (the boundary particle density and its normal derivative) are here replaced by the partial currents entering or leaving each computational cell. The intuitive physical meaning of such quantities facilitates the application of the response matrix technique. Only eigenvalue (criticality) problems will be here considered. As it regards the transport equation, the use of the so called Simplified Spherical Harmonics method (SPN) allows, through suitable approximations, to cast the problem into a system of differential elliptic equations of the diffusion type, which can still be solved by BEM.

A Boundary Element - Response Matrix Method for 3D Neutron Diffusion and Transport Problem

GIUSTI, VALERIO;
2014-01-01

Abstract

An application of a 3D Boundary Element Method (BEM), coupled with the Response Matrix (RM) technique, to solve the neutron diffusion and transport equations for multi-region domains is presented. The discussion is here limited to steady state problems, in which the neutrons have a wide energy spectrum, which leads to systems of several diffusion or transport equations. Moreover, the number of regions with different physical constants can be very large. The boundary integral equations concerning each region are solved via a polynomial moment expansion and, taking advantage of suitable recurrence formulas, the multi-fold integrals there involved are reduced to single or double integrals. The usual unknowns (the boundary particle density and its normal derivative) are here replaced by the partial currents entering or leaving each computational cell. The intuitive physical meaning of such quantities facilitates the application of the response matrix technique. Only eigenvalue (criticality) problems will be here considered. As it regards the transport equation, the use of the so called Simplified Spherical Harmonics method (SPN) allows, through suitable approximations, to cast the problem into a system of differential elliptic equations of the diffusion type, which can still be solved by BEM.
2014
Giusti, Valerio; Montagnini, B.
File in questo prodotto:
File Dimensione Formato  
cmes.2014.102.229.pdf

accesso aperto

Descrizione: Articolo finale
Tipologia: Versione finale editoriale
Licenza: Creative commons
Dimensione 263.46 kB
Formato Adobe PDF
263.46 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11568/651663
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact