RSICC CODE PACKAGE CCC-709
1. NAME AND TITLE
TDTORT: Time-Dependent, Three-Dimensional, Discrete Ordinates, Neutron
Transport Code System.
AUXILIARY ROUTINES
GIP: Group cross-section preparation.
2. CONTRIBUTORS
The University of Tennessee, Nuclear Engineering Department, Knoxville, Tennessee.
Oak Ridge National Laboratory, Oak Ridge, Tennessee.
3. CODING LANGUAGE AND COMPUTER
Fortran 77 and C; SUN Ultra; DEC OSF/1, PC under Linux (C00709/MNYWS/00).
4. NATURE OF PROBLEM SOLVED
TDTORT solves the time-dependent, three-dimensional neutron transport equation with explicit representation of delayed neutrons to estimate the fission yield from fissionable material transients. This release includes a modified version of TORT from the C00650MFMWS01 DOORS3.1 code package plus the time-dependent TDTORT code. GIP is also included for cross-section preparation.
TORT calculates the flux or fluence of particles due to particles incident upon the system from extraneous sources or generated internally as a result of interaction with the system in two- or three-dimensional geometric systems. The principle application is to the deep-penetration transport of neutrons and photons. Reactor eigenvalue problems can also be solved. Numerous printed edits of the results are available, and results can be transferred to output files for subsequent analysis.
TDTORT reads ANISN-format cross-section libraries, which are not included in
the package. Users may choose from several available in RSICC's data library
collection which can be identified by the keyword "ANISN FORMAT."
5. METHOD OF SOLUTION
The time-dependent spatial flux is expressed as a product of a space-, energy-, and angle-dependent shape function, which is usually slowly varying in time and a purely time-dependent amplitude function. The shape equation is solved for the shape using TORT; and the result is used to calculate the point kinetics parameters (e.g., reactivity) by using their inner product definitions, which are then used to solve the time-dependent amplitude and precursor equations. The amplitude function is calculated by solving the kinetics equations using the LSODE solver. When a new shape calculation is needed, the flux is calculated using the newly computed amplitude function.
The Boltzmann transport equation is solved using the method of discrete ordinates to treat the directional variable and weighted finite-difference methods, in addition to Linear Nodal and Linear Characteristic methods in TORT to treat spatial variables. Energy dependence is treated using a multigroup formulation. Starting in one corner of a mesh, at the highest energy, and with starting guesses for implicit sources, boundary conditions and recursion relationships are used to sweep into the mesh for each discrete direction independently. Integral quantities such as scalar flux are obtained from weighted sums of the directional results. The calculation then proceeds to lower energy groups, one at a time.
Iterations are used to resolve implicitness caused by scattering between
directions within a single energy group, by scattering from an energy group to
another group previously calculated, by fission, and by certain types of
boundary conditions. Methods are available to accelerate convergence of both
inner and outer iterations. Anisotropic scattering is represented by a Legendre
expansion of arbitrary order, and methods are available to mitigate the effect
of negative scattering source estimates resulting from finite truncation of the
expansion. Direction sets can be biased, concentrating work into directions of
particular interest. Fixed sources can be specified at either external or
internal mesh boundaries, or distributed within mesh cells.
6. RESTRICTIONS OR LIMITATIONS
TORT's limitations for solving a three-dimensional, fixed source problem
apply (i.e., geometry, convergence, non-linear effects, etc.). External forces
and nonlinear physical effects cannot be treated. Penetration through large,
non-scattering regions may become inaccurate due to ray effects. Problems with
scattering ratios near unity or eigenvalue calculations with closely spaced
eigenvalues may be quite time-consuming. Although flexible dimensioning is used
in TORT so that no fixed limits on group, problem size, etc., are applicable,
TDTORT uses a fixed size container array, which may not be big enough for very
large problems. The user should change the value of isdim in comsd3.F and
recompile. When the size is not sufficient, the code will indicate how much it
needs (which can be used to determine the new size).
7. TYPICAL RUNNING TIME
Running time is roughly proportional to the number of fixed source (flux
update) calculations each being proportional to mesh cells, directions, groups,
and iterations. The codes compiled in only a few minutes. Test cases ran in
about 47 minutes on a Sun UltraSparc 60, in 19 minutes on a DEC, and in 18
minutes on Pentium III 450Mhz computer.
8. COMPUTER HARDWARE REQUIREMENTS
TDTORT was developed on Sun workstations. It was compiled, and test cases
were run on DEC Alpha Digital Unix workstations and personal computers under
Linux.
9. COMPUTER SOFTWARE REQUIREMENTS
The codes run under Unix or Linux. Fortran and C compilers are required for installation on Unix systems. Executables created with the Portland Group Inc. Workstation 3.3 f77 compiler are included for Linux users. TDTORT was tested at RSICC on the following systems:
Sun SparcStation 60, OS5.6, SUN Fortran 77 Ver 5
DEC 500au, Digital UNIX 4.00, Digital Fortran 77 Version 5.6-075
Pentium III 450Mhz under RedHat Linux 6+ with PGI Workstation 3.3 f77
10. REFERENCES
a. included in hardcopy form:
Sedat Goluoglu, "A Deterministic Method for Transient, Three-Dimensional
Neutron Transport," Ph.D. Dissertation, Nuclear Engineering Department, The
University of Tennessee, Knoxville (August 1997).
b. included in electronic pdf format:
W. A. Rhoades and D. B. Simpson, "The TORT Three-Dimensional Discrete
Ordinates Neutron/Photon Transport Code," ORNL/TM-13221 (October 1997).
W. A. Rhoades and M. B. Emmett, the GIP section of "DOS: The Discrete
Ordinates System," ORNL/TM-8362 (September 1982).
c. background information:
S. Goluoglu, H.L. Dodds, "A Time-Dependent, Three-Dimensional Neutron
Transport Methodology," Nucl. Sci. Eng., 139, 248-261 (2001).
11. CONTENTS OF CODE PACKAGE
Included are the referenced documents and a GNU compressed Unix tar file on
CD. The tar file contains the source files for TDTORT, TORT, and GIP, plus Linux
executables, test cases, information file, and scripts.
12. DATE OF ABSTRACT
March 2002.
KEYWORDS: DISCRETE ORDINATES; NEUTRON; GAMMA-RAY; MULTIGROUP; ADJOINT; SPHERICAL GEOMETRY; SLAB; CYLINDRICAL GEOMETRY; COMPLEX GEOMETRY; KINETICS; TIME-DEPENDENT; WORKSTATION