2026-09-18T19:23:02Zhttps://dspace.mit.edu/server/oai/requestoai:dspace.mit.edu:1721.1/874902026-06-16T18:14:42Zcom_1721.1_7582com_1721.1_7581col_1721.1_131022
Benoit Forget.
Roberts, Jeremy Alyn
Massachusetts Institute of Technology. Department of Nuclear Science and Engineering.
Massachusetts Institute of Technology. Department of Nuclear Science and Engineering
2014-05-23T19:37:24Z
2014-05-23T19:37:24Z
2014
2014
http://hdl.handle.net/1721.1/87490
879667137
Thesis: Ph. D., Massachusetts Institute of Technology, Department of Nuclear Science and Engineering, 2014.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 195-200).
Modeling full reactor cores with high fidelity transport methods is a difficult task, requiring the largest computers available today. This thesis presents work on an alternative approach using the eigenvalue response matrix method (ERMM). The basic idea of ERMM is to decompose a reactor spatially into local "nodes." Each node represents an independent fixed source transport problem, and the nodes are linked via approximate boundary conditions to reconstruct the global problem using potentially many fewer explicit unknowns than a direct fine mesh solution. This thesis addresses several outstanding issues related to the ERMM based on deterministic transport. In particular, advanced transport solvers were studied for application to the relatively small and frequently repeated problems characteristic of response function generation. This includes development of preconditioners based on diffusion for use in multigroup Krylov linear solvers. These new solver combinations are up to an order of magnitude faster than competing algorithms. Additionally, orthogonal bases for space, angle, and energy variables were investigated. For the spatial variable, a new basis set that incorporates a shape function characteristic of pin assemblies was found to reduce significantly the error in representing boundary currents. For the angular variable, it was shown that bases that conserve the partial current at a boundary perform very well, particularly for low orders. For the deterministic transport used in this work, such bases require use of specialized angular quadratures. In the energy variable, it was found that an orthogonal basis constructed using a representative energy spectrum provides an accurate alternative to few group calculations. Finally, a parallel ERMM code Serment was developed, incorporating the transport and basis development along with several new algorithms for solving the response matrix equations, including variants of Picard iteration, Steffensen's method, and Newton's method. Based on results from several benchmark models, it was found that an accelerated Picard iteration provides the best performance, but Newton's method may be more robust. Furthermore, initial scoping studies demonstrated good scaling on an [omicron](100) processor machine.
by Jeremy Alyn Roberts.
Ph.D.
200 pages
eng
Massachusetts Institute of Technology
M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.
http://dspace.mit.edu/handle/1721.1/7582
Nuclear Science and Engineering.
Advanced response matrix methods for full core analysis
Thesis
application/pdf
false
Publication
unknown
unknown
<Publication xmlns="https://www.openaire.eu/cerif-profile/1.1/" id="f5bab53d-3044-423a-92b6-7dc49bcc0f0b">
<Type xmlns="https://www.openaire.eu/cerif-profile/vocab/COAR_Publication_Types">http://purl.org/coar/resource_type/c_1843</Type>
<Language>eng</Language>
<Title>Advanced response matrix methods for full core analysis</Title>
<PublishedIn>
<Publication>
</Publication>
</PublishedIn>
<PublicationDate>2014</PublicationDate>
<Authors>
<Author>
<DisplayName>Roberts, Jeremy Alyn</DisplayName>
<Affiliation>
<OrgUnit>
</OrgUnit>
</Affiliation>
</Author>
</Authors>
<Editors>
</Editors>
<Publishers>
<Publisher>
<DisplayName>Massachusetts Institute of Technology</DisplayName>
<OrgUnit />
</Publisher>
</Publishers>
<License>http://dspace.mit.edu/handle/1721.1/7582</License>
<Keyword>Nuclear Science and Engineering.</Keyword>
<Abstract>Modeling full reactor cores with high fidelity transport methods is a difficult task, requiring the largest computers available today. This thesis presents work on an alternative approach using the eigenvalue response matrix method (ERMM). The basic idea of ERMM is to decompose a reactor spatially into local "nodes." Each node represents an independent fixed source transport problem, and the nodes are linked via approximate boundary conditions to reconstruct the global problem using potentially many fewer explicit unknowns than a direct fine mesh solution. This thesis addresses several outstanding issues related to the ERMM based on deterministic transport. In particular, advanced transport solvers were studied for application to the relatively small and frequently repeated problems characteristic of response function generation. This includes development of preconditioners based on diffusion for use in multigroup Krylov linear solvers. These new solver combinations are up to an order of magnitude faster than competing algorithms. Additionally, orthogonal bases for space, angle, and energy variables were investigated. For the spatial variable, a new basis set that incorporates a shape function characteristic of pin assemblies was found to reduce significantly the error in representing boundary currents. For the angular variable, it was shown that bases that conserve the partial current at a boundary perform very well, particularly for low orders. For the deterministic transport used in this work, such bases require use of specialized angular quadratures. In the energy variable, it was found that an orthogonal basis constructed using a representative energy spectrum provides an accurate alternative to few group calculations. Finally, a parallel ERMM code Serment was developed, incorporating the transport and basis development along with several new algorithms for solving the response matrix equations, including variants of Picard iteration, Steffensen's method, and Newton's method. Based on results from several benchmark models, it was found that an accelerated Picard iteration provides the best performance, but Newton's method may be more robust. Furthermore, initial scoping studies demonstrated good scaling on an [omicron](100) processor machine.</Abstract>
<Access xmlns="http://purl.org/coar/access_right"
>
</Access>
</Publication>