RSICC CODE PACKAGE CCC-708
1. NAME and TITLE
REBUS-PC 1.4: Code System for Analysis of Research Reactor Fuel Cycles.
REBUS-PC contains DIF3D 9.0/VARIANT9.0, a code system using variational
nodal methods and finite difference methods to solve neutron diffusion
and transport theory problems. This release is an update of CCC-649 DIF3D8.0/VARIANT8.0,
which was provided by the authors of that code. REBUS-PC provides new capabilities
beyond that of CCC-653, REBUS3/VARIANT8.0. REBUS-PC was developed for analysis
of research reactor cores and fuel cycles, but it remains generally useful
for any reactor type.
RERTR Program, Technology Development Division, Argonne National Laboratory,
Argonne, Illinois. The included DIF3D/VARIANT9.0 code was provided by the
Reactor Analysis and Engineering Division, Argonne National Laboratory,
3. CODING LANGUAGE and COMPUTER
Fortran 77, PC under either WINDOWS or Linux (C00709/PC586/00).
4. NATURE of PROBLEM SOLVED
REBUS-PC Version 1.4 is a system of codes designed for the analysis of research reactor fuel cycles. It is based on an updated 9.0 version of DIF3D that is similar to RSICC Code Package CCC-653, REBUS-3/VARIANT 8.0 (which is intended for use on unix workstations). The full capabilities of the workstation version are retained and enhanced for use in a PC environment. Two basic types of analysis problems are solved: 1) the infinite-time, or equilibrium, conditions of a reactor operating under a fixed fuel management scheme, or 2) the explicit cycle-by-cycle, or non-equilibrium operation of a reactor under a specified periodic or non-periodic fuel management program. For the equilibrium type problems, the code uses specified external fuel supplies to load the reactor. Optionally, reprocessing may be included in the specification of the external fuel cycle and discharged fuel may be recycled back into the reactor. For non-equilibrium cases, the initial composition of the reactor core may be explicitly specified or the core may be loaded from external feeds and discharged fuel may be recycled back into the reactor as in equilibrium problems.
Four types of search procedures may be carried out in order to satisfy user-supplied constraints: 1) adjustment of the reactor burn cycle time to achieve a specified discharge burnup, 2) adjustment of the fresh fuel enrichment to achieve a specified multiplication constant at a specified point during the burn cycle, 3) adjustment of the control poison density to maintain a specified value of the multiplication constant throughout the reactor burn cycle, and 4) adjustment of the reactor burn cycle time to achieve a specified value of the multiplication constant at the end of the burn step.
REBUS-PC will handle both equilibrium and non-equilibrium problems using a number of different core geometries including triangular and hexagonal mesh. The neutronics solution may be obtained using finite difference or nodal diffusion-theory methods. Other features include: fully automatic restart capability, no restrictions on number of neutron energy groups, and general external cycle with no restrictions on number of external feeds, reprocessing plants, etc. Fuel management is completely general for non-equilibrium problems.
Microscopic cross sections are permitted to vary as a function of the atom density of various reference isotopes in the problem as appropriate for thermal reactor systems. The previous capability where neutron capture and fission processes were fitted to low-order polynomials as functions of burnup is retained. In addition, the user now may select cubic spline interpolation for (n,gamma), (n,fission), (n,alpha), (n,p), (n,d), and (n,t) reactions as functions of burnup. The user may specify control rod positions at each time node in the problem. Output edits have been extensively revised and better organized for use in a PC environment. A number of ASCII format datasets containing various types of summary results are available for use in tailoring reports with the aid of auxiliary PC software such as spreadsheet or word processor programs.
This is a standalone and expanded version of the modular REBUS-3 code system described in Refs. 1-7. It utilizes the 9.0 update of CCC-649/DIF3D code to obtain the neutronics solution. All but the main program "path driver" of DIF3D 9.0/VARIANT 9.0 (with finite difference and hexagonal nodal option) is included in this code package. Dataset A.ISO is extended to retain full single-precision accuracy for neutron cross sections in BCD format. REBUS-PC operation is fully compatible with the DOE Committee on Computer Code Coordination coding standards and interface data sets [Ref. 10].
Changes were made to power-related edits in the fuel cycle computational
module FCC004 such that they are now consistent with those from DIF3D.
The cross section homogenization module HMG4C was extended to provide the
necessary information. HMG4C is now called after a neutronics solution
in order to write the necessary data to file COMPXS. FCC004 is then called
to create power-related edits (using COMPXS) and to predict the next burn
step isotopic compositions.
5. METHOD of SOLUTION
The total reactor burn cycle time is divided into one or more subintervals,
the number of which is specified by the user. An explicit burnup is performed
in each region of the reactor over each of these subintervals using the
average reaction rates over the subinterval. These average reaction rates
are based on fluxes obtained from an explicit 1-, 2-, or 3-dimensional
finite-difference or nodal diffusion theory solution computed at both the
beginning and end of the subinterval, iterated to a prescribed convergence.
The transmutation equations are solved by the matrix-exponential technique.
The isotopes to be considered in the burnup equations, as well as their
transmutation reactions, are specified by the user. Burnup-dependent transmutation
reaction cross sections can be used.
6. RESTRICTIONS or LIMITATIONS
Very large problems can be solved. LF95 v6.1 is limited to 4 Gigabytes
of available memory, 250 files open concurrently, and a maximum file size
of 2 Gigabytes. Variable dimensioning is used throughout.
7. TYPICAL RUNNING TIME
Minutes to many hours depending on size and complexity of the problem,
and on CPU speed and memory utilized. An equilibrium problem with 84 burnable
regions, 7 neutron groups, an X-Y-Z mesh of 78x69x23, requiring 9 neutronics
solutions over 4 burn steps, took 34 minutes on a 2.0GHz Pentium IV PC
running under Linux. Each sample problem required over 1 GB memory and
about 2 hours execution time on a P3 650mHz with 128 MB RAM running Windows
8. COMPUTER HARDWARE REQUIREMENTS
The PC version of the code is in production use at Argonne National Laboratory on PC's using Pentium IV processors under the Red Hat Linux 7.2 operating system, using up to 1 Gigabyte of random access memory. The code can be compiled for use on any PC from Intel 80486, Pentium, Pentium MMX, Pentium Pro, Pentium II, Pentium III, Pentium IV, and Celeron processors or their generic counterparts.
At least 32 MB of RAM are needed for compilation. 128 Mbytes of RAM
are recommended for program and file buffer storage, and internal data,
when running the code. More RAM is better, as wall clock time is generally
reduced when all or most of the problem is contained in core memory (limiting
paging). External data storage must be available for approximately 40 scratch
and interface files. Fourteen of these files are random access scratch
files (grouped into 6 file groups). The remainder are sequential access
files with formatted or unformatted record types. The code is of modular
design, but standalone form. Dynamic memory allocation uses two containers
whose sizes are dictated by the needs of a given problem and by the user's
choice depending upon code options and hardware limitations.
9. COMPUTER SOFTWARE REQUIREMENTS
The code is written entirely in Fortran 77. The program can be compiled
with Lahey Fortran 95 Pro V6.1 for Linux or Lahey Fortran 95 V5.60g or
V5.56 for Windows. Executables created with these Lahey compilers are included
in the package for Linux and Windows. Other advanced Fortran compilers
could be used but would require minor changes in dynamic memory management
calls and clock timer routines, or anywhere else where the particular Fortran
dialect differs from Lahey Fortran 95. The standalone source code contains
approximately 1140 subroutines and 262,000 Fortran statements. The operating
system can be any of Windows 95/98/NT4.0/2000 or newer, or any linux variation
such as Red Hat Linux 7.2 that is compatible with the user's Fortran compiler.
An interpreter or manager or executive routine is not needed. No nonstandard
library routines are used. No special I/O table is required. Output file
FT06 is actually created on unit 4, in order to reserve unit 6 for screen
a Included in documentation:
A. P. Olson, "A User's Guide for the REBUS-PC Code, Version 1.4," Argonne
National Laboratory, (December 21, 2001).
b. Background information:
1. J. Hoover, G. K. Leaf, D. A. Meneley, and P. M. Walker, "The Fuel
Cycle Analysis System, REBUS," Nucl. Sci. Eng. 45, p. 53, 1971.
2. A. P. Olson, J. P. Regis, D. A. Meneley, and L. J. Hoover, "A User's
Manual for the Reactor Burnup System, REBUS," FRA-TM-41, Argonne National Laboratory, September 28, 1972.
3. A. P. Olson, "A User's Manual for the Reactor Burnup System,
REBUS-2," FRA-TM-62, Argonne National Laboratory, March 1, 1974.
4. R. P. Hosteny, "The ARC System Fuel Cycle Analysis Capability, REBUS-2," ANL-7721 (October 1978) distributed in file rebus-pc/documents/rebus3_docs/ANL-7221.pdf.
5. B. J. Toppel, "A User's Guide for the REBUS-3 Fuel Cycle Analysis Capability," ANL-83-2 (1983) distributed in file rebus-pc/documents/rebus3_docs/REBUS-3_document.pdf.
6. K. L. Derstine, "DIF3D, A Code to Solve One-, Two-, and Three-
Dimensional Finite-Difference Diffusion Theory Problems," ANL-82-64 (April 1984) distributed in directory rebus-pc/documents/rebus3_docs/D3D.
7. G. Palmiotti, E. E. Lewis, and C. B. Carrico, "VARIANT: VARIational
Anistropic Nodal Transport for Multidimensional Cartesian and Hexagonal Geometry Calculation," ANL-95/40 (October 1995) distributed in file rebus-pc/documents/rebus3_docs/Var.manual.pdf.
8. J. F. Breismeister, ed., "MCNP - A General Monte Carlo N-Particle
Transport Code, Version 4B," LA-12625-M. Los Alamos National Laboratory, 1997.
9. N. A. Hanan, A. P. Olson, R. B. Pond, and J. E. Matos, "A Monte Carlo
Burnup Code Linking MCNP and REBUS," presented at the 1998 International Meeting on Reduced Enrichment for Research and Test Reactors, October 18-23, Sao Paolo, Brazil.
10. R. D. O'Dell, "Standard Interface Files and Procedures for Reactor
Physics Codes, Version IV," LA-6941-MS, UC-32, Los Alamos National Laboratory, September 1977.
11. R. D. Lawrence, "The DIF3D Nodal Neutronics Option for Two-and Three-Dimensional Diffusion Theory Calculations in Hexagonal Geometry," ANL-83-1 (March 1983) distributed in directory rebus-pc/documents/rebus3_docs/D3DN.
12. C. H. Adams, et.al., "The Utility Subroutine Package Used by Applied
Physics Division Export Codes," ANL-83-3 (May 1992) distributed in file
11. CONTENTS of CODE PACKAGE
Included are the reference in 10.a and the references in 10.b which
indicate they are distributed in electronic files and a CD which contains
a self-extracting compressed Windows file and a GNU compressed tar file.
Approximately 360 MB is required for installation. Included in the distribution
files are the source code, Lahey-compiled Windows and Linux executables,
sample problem input and output, and code dependent BCD and binary card-image
12. DATE of ABSTRACT
KEYWORDS: BURNUP; FUEL MANAGEMENT; DIFFUSION THEORY; CRITICALITY CALCULATIONS; REACTOR PHYSICS; COMPLEX GEOMETRY; CCCC INTERFACE FORMAT