Engineering topics
Zhu, Ang
Publications and source records attributed to Zhu, Ang.
MPACT Standard Input User s Manual, Version 2.2.0
The MPACT (Michigan PArallel Charactistics based Transport) code is designed to perform high-fidelity light water reactor (LWR) analysis using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. The code consists of several libraries which provide the functionality necessary to solve steady-state eigenvalue problems. Several transport capabilities are available within MPACT including both 2-D and 3-D Method of Characteristics (MOC). A three-dimensional whole core solution based on the 2D-1D solution method provides the capability for full core depletion calculations.
MPACT Performance Improvements in VERA-CS
The report summarizes the initial profile of MPACT representative of the version that was distributed with VERA 3.4. It then describes several activities recently performed to increase the overall performance of the code. For the various improvements described in this report, several have been well known to the development team and reactor physics methods community for some time; while others proved novel in their combination of other existing ideas. Still, other improvements were the result of unanticipated breakthroughs in numerical methods research. These performance improvement activities focused on the following components of MPACT. Method of Characteristics Transport Kernel Iterative Methods for Coarse Mesh Finite Difference Linear System Macroscopic Cross Section and Effective Fission Spectrum Calculation Resonance self-shielding (Subgroup) calculation A summary of the overall improvement from this work as tested on a few representative problems is provided in the table below. The incremental improvements for various changes are shown in the figure below.
MPACT Theory Manual, Version 2.2.0
This theory manual describes the three-dimensional (3-D) whole-core, pin-resolved transport calculation methodology employed in the MPACT code. To provide sub-pin level power distributions with sufficient accuracy, MPACT employs the method of characteristics (MOC) solutions in the framework of a 3-D coarse mesh finite difference (CMFD) formulation. MPACT provides a 3D MOC solution, but also a 2D/1D solution in which the 2D planar solution is provided by MOC and the axial coupling is resolved by one-dimensional (1-D) lower order (diffusion or P3) solutions. In Chapter 2 of the manual, the MOC methodology is described for calculating the regional angular and scalar fluxes from the Boltzmann transport equation. In Chapter 3, the 2D/1D methodology is described, together with the description of the CMFD iteration process involving dynamic homogenization and solution of the multigroup CMFD linear system. A description of the MPACT depletion algorithm is given in Chapter 4, followed by a discussion of the subgroup and ESSM resonance processing methods in Chapter 5. The final Chapter 6 describes a simplified thermal hydraulics model in MPACT.
MPACT VERA Input User s Manual, Version 2.2.0
The MPACT (Michigan PArallel Charactistics based Transport) code is designed to perform high-fidelity light water reactor (LWR) analysis using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. The code consists of several libraries which provide the functionality necessary to solve steady-state eigenvalue problems. Several transport capabilities are available within MPACT including both 2-D and 3-D Method of Characteristics (MOC). A three-dimensional whole core solution based on the 2D-1D solution method provides the capability for full core depletion calculations.
A Hybrid Nodal P3/SP3 Axial Transport Solver for the MPACT 2D/1D Scheme
Explore the source record for details and available documents.
MPACT User's Manual (Version 2.0.0)
The MPACT (Michigan PArallel Charactistics based Transport) code is designed to perform high-fidelity light water reactor (LWR) analysis using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. The code consists of several libraries which provide the functionality necessary to solve steady-state eigenvalue problems. Several transport capabilities are available within MPACT including both 2-D and 3-D Method of Characteristics (MOC). A three-dimensional whole core solution based on the 2D-1D solution method provides the capability for full core depletion calculations.
Demonstration of Full Core Reactor Depletion with MPACT
MPACT is a three-dimensional (3-D) whole core transport code capable of generating sub-pin level power distributions. This feature is accomplished by obtaining the integral transport solutions to the heterogeneous reactor problem in which the actual detailed geometrical configuration of fuel components such as the pellet and cladding is modelled explicitly.. The cross section data needed for the neutron transport calculation are obtained directly from a multigroup microscopic cross section library similar to those used in lattice physics codes. Hence MPACT involves neither a priori homogenization nor group condensation for the core spatial solution. The integral transport solution is obtained by means of the method of characteristics (MOC) that employs discrete ray tracing. Since the direct application of the 3-D MOC capability in MPACT to 3-D core configuration requires considerable amounts of memory and computing time for practical reactor applications, an alternative approximate 3-D solution method was implemented in MPACT based on a 2D-1D approach which employs planar MOC solutions in the framework of the 3-D coarse mesh finite difference (CMFD) formulation. The axial coupling is resolved by one-dimensional (1-D) diffusion solutions and the planar and axial problems are coupled through the transverse leakage. The use of a lower order 1-D solution in the axial direction is justified by the fact that most heterogeneity in the core occurs in the radial direction rather than the axial. However, it is possible to use one of the higher order methods in MPACT such as the 1-D SN kernel if improved solution accuracy is desired. The CMFD formulation which originally served as an efficient nodal formulation is used in MPACT to accelerate the whole core transport calculation. The basic mesh in the CMFD formulation is a pin cell which is much coarser than the flat source regions defined for MOC calculations. The concept of dynamic homogenization of group constants for the pin cells provides the basis of the CMFD formulation in whole core transport calculation. The intra-cell flux distribution determined from the MOC calculation is used to generate the homogenized cell constants while the MOC cell surface averaged currents are used to determine the radial nodal coupling coefficients. This dynamic implementation of the equivalence theory formulation ensures the same transport solution is obtained with CMFD as the one obtained using unaccelerated MOC calculations alone. The ability of MPACT to predict the Hot, Zero Power (HZP) condition of a reactor was reported in CASL-U-2014-0045-000. The purpose of this report is to demonstrate the capability of MPACT to simulate the operation of a reactor at Hot, Full Power (HFP) conditions during normal operation throughout a full operating cycle. The time scale for a reactor burnup cycle is typically one to two years and therefore short time phenomena such as the effect of delayed neutrons can be neglected. Instead the time dependence can be treated quasi-statically using the depletion algorithm presented in CASL-U-2013-0276-001. The full power operation requires thermal hydraulic feedback to be present during the entire simulation which will utilize the work in CASL-U-2014-0051-000. However, due to the significant computational burden of the coupled physics simulation, a simplified thermal hydraulic model is used for the analysis reported here which is principally intended to demonstrate the full core neutronics depletion. In addition to the thermal hydraulic feedback, the effect of Boron 10 depletion should be modeled in cases where major boron additions are not performed for extended periods of time. Each of these capabilities will be described in the following sections.