Altmetric

A geometry preserving, conservative, mesh-to-mesh isogeometric interpolation algorithm for spatial adaptivity of the multigroup, second-order even-parity form of the neutron transport equation

File Description SizeFormat 
1-s2.0-S0021999117304606-main.pdfPublished version1.16 MBAdobe PDFView/Open
Title: A geometry preserving, conservative, mesh-to-mesh isogeometric interpolation algorithm for spatial adaptivity of the multigroup, second-order even-parity form of the neutron transport equation
Authors: Welch, J
Kophazi, J
Owens, A
Eaton, M
Item Type: Journal Article
Abstract: In this paper a method is presented for the application of energy-dependent spatial meshes applied to the multigroup, second-order, even-parity form of the neutron transport equation using Isogeometric Analysis (IGA). The computation of the inter-group regenerative source terms is based on conservative interpolation by Galerkin projection. The use of Non-Uniform Rational B-splines (NURBS) from the original computer-aided design (CAD) model allows for efficient implementation and calculation of the spatial projection operations while avoiding the complications of matching different geometric approximations faced by traditional finite element methods (FEM). The rate-of-convergence was verified using the method of manufactured solutions (MMS) and found to preserve the theoretical rates when interpolating between spatial meshes of different refinements. The scheme’s numerical efficiency was then studied using a series of two-energy group pincell test cases where a significant saving in the number of degrees-of-freedom can be found if the energy group with a complex variation in the solution is refined more than an energy group with a simpler solution function. Finally, the method was applied to a heterogeneous, seven-group reactor pincell where the spatial meshes for each energy group were adaptively selected for refinement. It was observed that by refining selected energy groups a reduction in the total number of degrees-of-freedom for the same total L2 error can be obtained.
Issue Date: 27-Jun-2017
Date of Acceptance: 7-Jun-2017
URI: http://hdl.handle.net/10044/1/48627
DOI: https://dx.doi.org/10.1016/j.jcp.2017.06.015
ISSN: 0021-9991
Publisher: Elsevier
Start Page: 129
End Page: 146
Journal / Book Title: Journal of Computational Physics
Volume: 347
Copyright Statement: © 2017 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Sponsor/Funder: Engineering & Physical Science Research Council (EPSRC)
Rolls-Royce Plc
Engineering & Physical Science Research Council (E
Funder's Grant Number: EP/J002011/1
PO 5001626145
EP/K503733/1
Keywords: Applied Mathematics
01 Mathematical Sciences
02 Physical Sciences
09 Engineering
Publication Status: Published
Appears in Collections:Faculty of Engineering
Mechanical Engineering



Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.

Creative Commons