Implementation of CMFD Acceleration Scheme in PROTEUS-MOC

Journal Article · Sat Jul 01 04:00:00 EDT 2017 · Transactions of the American Nuclear Society

OSTI ID:23050340

Hsieh, Albert; Yang, Won Sik ^[1]

---

1. School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907 (United States)

The PROTEUS-MOC code is a three-dimensional (3D) transport code, which is being developed at Argonne National Laboratory under the DOE-NE's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. It provides a fully consistent 3D discretization of the transport equation by combining the two-dimensional (2D) method of characteristics (MOC) for the planar direction with the discontinuous Galerkin method for the treatment of the axial variable. PROTEUS-MOC uses the unstructured finite element method (FEM) and thus allows the explicit representation of any extruded geometries. It also has a capability of massively parallel computation in energy group, angle, plane, and space using a top-level GMRES linear algebra solver. For parallel computation in space, it uses an unstructured domain decomposition by adapting METIS. The computational burden for direct heterogeneous core calculations based on unstructured geometries is excessive for practical reactor problems, and hence it is essential to implement an efficient acceleration scheme in PROTEUS-MOC. Recently, the coarse mesh finite difference (CMFD) method, which was originally applied to the nodal diffusion methods, has been successfully applied to the MOC transport solvers, resulting in a remarkable reduction in the number of fission source iterations. In addition, a modified CMFD method has been developed to overcome the limitation of the conventional CMFD that is applicable only to neutronics solvers with a conservative discretization. This CMFD method was derived to be consistent with the FEM, which is based on a nonconservative discretization. In this method, the neutron balance in each coarse mesh is enforced by introducing a pseudo absorption cross section (PAXS), and the well-established alternating solution process of CMFD and transport calculations is employed to accelerate source convergence. This modified CMFD acceleration scheme was adopted in PROTEUS-MOC to accelerate the fission source iteration. Pin-cell and assembly based CMFD schemes have been implemented by defining individual pin cells or assemblies as a coarse mesh. For this, an algorithm to map the finite elements to coarse meshes was developed, and the routines were developed to determine the coupling coefficients and PAXS of the CMFD equations from the transport solution. Preliminary performance tests were carried out using the C5G7 benchmark problem. This paper presents the main features of this CMFD acceleration scheme and the preliminary test results.

• OSTI ID: 23050340
• Journal Information: Transactions of the American Nuclear Society, Journal Name: Transactions of the American Nuclear Society Vol. 116; ISSN 0003-018X
• Country of Publication: United States
• Language: English

Related Subjects

22 GENERAL STUDIES OF NUCLEAR REACTORS
97 MATHEMATICS AND COMPUTING
ABSORPTION
ALGEBRA
ALGORITHMS
ANL
BENCHMARKS
COMPUTERIZED SIMULATION
CROSS SECTIONS
FINITE ELEMENT METHOD
FISSION
GEOMETRY
IMPLEMENTATION
NEUTRONS
NUCLEAR ENERGY
PERFORMANCE
THREE-DIMENSIONAL CALCULATIONS
TRANSPORT THEORY
TWO-DIMENSIONAL CALCULATIONS