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

被引：9
作者：
Latimer, C. [1 ]
Kophazi, J. [1 ]
Eaton, M. D. [1 ]
McClarren, R. G. [2 ]
机构：
[1] Imperial Coll London, Dept Mech Engn, Nucl Engn Grp, City & Guilds Bldg,Exhibit Rd, London SW7 2AZ, England
[2] Univ Notre Dame, Dept Aerosp & Mech Engn, Fitzpatrick Hall, Notre Dame, IN 46556 USA
来源：
ANNALS OF NUCLEAR ENERGY | 2020年 / 136卷
基金：
英国工程与自然科学研究理事会;
关键词：
Isogeometric analysis; SAAF; Discrete ordinates; Variational principle; SPECTRAL ELEMENT METHOD; PETROV-GALERKIN METHODS; DIFFUSION EQUATION; FINITE-ELEMENTS; 1ST-ORDER FORM; NURBS; REFINEMENT;
D O I：
10.1016/j.anucene.2019.107049
中图分类号：
TL [原子能技术]; O571 [原子核物理学];
学科分类号：
0827 ; 082701 ;
摘要：
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. (C) 2019 Elsevier Ltd. All rights reserved.
引用
收藏
页数：16
相关论文
共 63 条
  • [11] A spherical harmonics - Finite element discretization of the self-adjoint angular flux neutron transport equation
    Cao, Liangzhi
    Wu, Hongchun
    [J]. NUCLEAR ENGINEERING AND DESIGN, 2007, 237 (23) : 2232 - 2239
  • [12] The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers
    Collier, Nathan
    Pardo, David
    Dalcin, Lisandro
    Paszynski, Maciej
    Calo, V. M.
    [J]. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2012, 213 : 353 - 361
  • [13] Isogeometric analysis of structural vibrations
    Cottrell, J. A.
    Reali, A.
    Bazilevs, Y.
    Hughes, T. J. R.
    [J]. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2006, 195 (41-43) : 5257 - 5296
  • [14] BDDC PRECONDITIONERS FOR ISOGEOMETRIC ANALYSIS
    da Veiga, L. Beirao
    Cho, D.
    Pavarino, L. F.
    Scacchi, S.
    [J]. MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES, 2013, 23 (06) : 1099 - 1142
  • [15] GeoPDEs: A research tool for Isogeometric Analysis of PDEs
    de Falco, C.
    Reali, A.
    Vazquez, R.
    [J]. ADVANCES IN ENGINEERING SOFTWARE, 2011, 42 (12) : 1020 - 1034
  • [16] Drumm C, 2015, ANS M C 2015
  • [17] Fan W.C., 2002, SAND20021880 SAND NA
  • [18] Multigrid methods for isogeometric discretization
    Gahalaut, K. P. S.
    Kraus, J. K.
    Tomar, S. K.
    [J]. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2013, 253 : 413 - 425
  • [19] MOOSE: A parallel computational framework for coupled systems of nonlinear equations
    Gaston, Derek
    Newman, Chris
    Hansen, Glen
    Lebrun-Grandie, Damien
    [J]. NUCLEAR ENGINEERING AND DESIGN, 2009, 239 (10) : 1768 - 1778
  • [20] The application of isogeometric analysis to the neutron diffusion equation for a pincell problem with an analytic benchmark
    Hall, S. K.
    Eaton, M. D.
    Williams, M. M. R.
    [J]. ANNALS OF NUCLEAR ENERGY, 2012, 49 : 160 - 169
1234567