The spectral element method for static neutron transport in A approximation. Part II

Barbarino, Andrea; Dulla, Sandra; Mund, E. H.; Ravetto, Piero

doi:10.1016/j.anucene.2015.04.015

Cited by 2 publications

(1 citation statement)

References 10 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In fact most of the unstructured mesh neutron transport literature is devoted to the development of FE based discretisation schemes. However, more recently high-order accurate spectral element (SE) [29,30,31,32] and isogeometric analyis (IGA) methods have been developed for modelling neutron transport [33,34,35,36,37,38]. The IGA method is a spatial discretisation methodology first introduced in 2005 by T. J. R Hughes [39].…”

Section: Introductionmentioning

confidence: 99%

A geometry conforming isogeometric method for the self-adjoint angular flux (SAAF) form of the neutron transport equation with a discrete ordinate (SN) angular discretisation

Latimer

Kópházi

Eaton

et al. 2020

Annals of Nuclear Energy

View full text Add to dashboard Cite

This paper presents the application of isogeometric analysis (IGA) to the spatial discretisation of the multi-group, self-adjoint angular flux (SAAF) form of the neutron transport equation with a discrete ordinate (S N) angular discretisation. The IGA spatial discretisation is based upon non-uniform rational B-spline (NURBS) basis functions for both the test and trial functions. In addition a source iteration compatible maximum principle is used to derive the IGA spatially discretised SAAF equation. It is demonstrated that this maximum principle is mathematically equivalent to the weak form of the SAAF equation. The rate of convergence of the IGA spatial discretisation of the SAAF equation is analysed using a method of manufactured solutions (MMS) verification test case. The results of several nuclear reactor physics verification benchmark test cases are analysed. This analysis demonstrates that for higher-order basis functions, and for the same number of degrees of freedom, the FE based spatial discretisation methods are numerically less accurate than IGA methods. The difference in numerical accuracy between the IGA and FE methods is shown to be because of the higher-order continuity of NURBS basis functions within a NURBS patch as well as the preservation of both the volume and surface area throughout the solution domain within the IGA spatial discretisation. Finally, the numerical results of applying the IGA SAAF method to the OECD/NEA, seven-group, two-dimensional C5G7 quarter core nuclear reactor physics verification benchmark test case are presented. The results, from this verification benchmark test case, are shown to be in good agreement with solutions of the first-order form as well as the second-order even-parity form of the neutron transport equation for the same order of discrete ordinate (S N) angular approximation.

show abstract