Skip to Main content Skip to Navigation

Acceleration and higher order schemes of a characteristic solver for the solution of the neutron transport equation in 3D axial geometries

Daniele Sciannandrone 1
1 LTSD - Laboratoire de Transport Stochastique et Déterministe
SERMA - Service des Réacteurs et de Mathématiques Appliquées : DEN/DM2S/SERMA
Abstract : The topic of our research is the application of the Method of Long Characteristics (MOC) to solve the Neutron Transport Equation in three-dimensional axial geometries. The strength of the MOC is in its precision and versatility. As a drawback, it requires a large amount of computational resources. This problem is even more severe in three-dimensional geometries, for which unknowns reach the order of tens of billions for assembly-level calculations.The first part of the research has dealt with the development of optimized tracking and reconstruction techniques which take advantage of the regularities of three-dimensional axial geometries. These methods have allowed a strong reduction of the memory requirements and a reduction of the execution time of the MOC calculation.The convergence of the iterative scheme has been accelerated with a lower-order transport operator (DPN) which is used for the initialization of the solution and for solving the synthetic problem during MOC iterations.The algorithms for the construction and solution of the MOC and DPN operators have been accelerated by using shared-memory parallel paradigms which are more suitable for standard desktop working stations. An important part of this research has been devoted to the implementation of scheduling techniques to improve the parallel efficiency.The convergence of the angular quadrature formula for three-dimensional cases is also studied. Some of these formulas take advantage of the reduced computational costs of the treatment of planar directions and the vertical direction to speed up the algorithm.The verification of the MOC solver has been done by comparing results with continuous-in-energy Monte Carlo calculations. For this purpose a coupling of the 3D MOC solver with the Subgroup method is proposed to take into account the effects of cross sections resonances. The full calculation of a FBR assembly requires about 2 hours of execution time with differences of few PCM with respect to the reference results.We also propose a higher order scheme of the MOC solver based on an axial polynomial expansion of the unknown within each mesh. This method allows the reduction of the meshes (and unknowns) by keeping the same precision.All the methods developed in this thesis have been implemented in the APOLLO3 version of the neutron transport solver TDT.
Complete list of metadatas

Cited literature [61 references]  Display  Hide  Download
Contributor : Abes Star :  Contact
Submitted on : Monday, January 11, 2016 - 8:30:06 AM
Last modification on : Tuesday, April 28, 2020 - 11:28:10 AM
Long-term archiving on: : Tuesday, April 12, 2016 - 11:04:38 AM


Version validated by the jury (STAR)


  • HAL Id : tel-01253553, version 1



Daniele Sciannandrone. Acceleration and higher order schemes of a characteristic solver for the solution of the neutron transport equation in 3D axial geometries. Computational Physics [physics.comp-ph]. Université Paris Sud - Paris XI, 2015. English. ⟨NNT : 2015PA112171⟩. ⟨tel-01253553⟩



Record views


Files downloads