RSIC CODE PACKAGE CCC-650
1. NAME AND TITLE
DOORS3.2a: One, Two- and Three-Dimensional Discrete Ordinates Neutron/Photon Transport Code System.
AUXILIARY ROUTINES
TORT: Three-dimensional neutron/photon transport.
DORT: Two-dimensional neutron/photon transport.
ANISN: One-dimensional neutron/photon transport.
GBANISN: ANISN with “group band” option.
JDOS: Script for performing calculations.
DRV: DOORS driver module called by JDOS.
ALC: Material cross-section library maintenance.
GIP: Group cross-section preparation.
TORSED: Couple R-Z DORT to X-Y-Z TORT calculation.
TORSET: Couple primary TORT to secondary TORT calculation.
VISA: Prepare a TORSED input file from a DORT output file.
GRTUNCL: Uncollided flux and first collision source.
FALSTF: Last collision source projection to point detector.
BNDRYS: Convert internal boundary flux to internal boundary source.
RTFLUM: Convert flux moment files between various formats.
ISOPLOT: Display 2-D flux or response contours.
XTORID: Prepare 2-D flux slices from TORT for display in ISPL3D.
CMP: Sandia National Lab. code system maintenance and update.
RSCORS: Sandia National Lab. subroutine library of graphical primitives.
2. CONTRIBUTOR
Oak Ridge National Laboratory, Oak Ridge, Tennessee.
3. CODING LANGUAGE AND COMPUTER
Fortran and C; IBM RS/6000, Sun, DEC Alpha, and on PC under Linux and Windows (C00650MFMWS04).
4. NATURE OF PROBLEM SOLVED
The DOORS3.2a discrete ordinates transport code system includes the most recent versions of CCC-543/TORT-DORT, CCC-254/ANISN-ORNL, CCC-628/GBANISN and CCC-351/FALSTF. It also includes the ISOPLOT code from the PSR-155/DOGS package and various utility programs listed above which were previously included in the TORT-DORT package. In this release each module is a separate executable file. Several modules, as needed, can be run in a single job by using “jdos” to call the “drv” module which interprets the sequence specified in the input.
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, and DORT is used in one- or two-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.
ANISN solves the one-dimensional Boltzmann transport equation for neutrons or gamma rays in slab, spherical, or cylindrical geometry. GBANISN is based on ANISN but was modified to allow randomizing of multibank fluxes within a group at all interfaces between dissimilar materials and a reduction in the number of outer iterations for problems involving neutron upscatter into higher energy groups. GBANISN requires fewer outer iterations than ANISN by performing "inner iterations" over energy groups within a "band" and converging those groups before moving to the next band. These "inner" iterations slightly resemble outer iterations in ANISN. Thus, a calculation with upscatter and no fission can be solved with one traditional outer iteration. GBANISN, like ANISN, includes a technique for handling general anisotropic scattering, pointwise convergence criteria, and alternate step function difference equations that effectively remove the oscillating flux distributions sometimes found in discrete ordinates solutions.
ISOPLOT uses the Sandia National Laboratory RSCORS graphical system. SNL’s CMP system for code maintenance is used to build the Fortran source files for the target computer. The RSCORS library (librscors.a) can be used to build executables for PSR-524/BOT3P2.0 and PSR-525/BOT3P3.0 (Bologna Transport Analysis Pre‑Post‑Processors)
DOORS 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”.
In October 2003, the package was updated to DOORS3.2a to meet increasing requests for support on Linux and Windows systems. The remainder of the package, other than relatively minor modifications to support PC architectures and current compilers, is unchanged.
5. METHOD OF SOLUTION
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. Time dependence is not treated. 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.
Two multitasking options are available for solving the Linear Nodal method equations in TORT on UNICOS Cray platforms: Octant Parallel (OP) which concurrently performs the right ® left and left ® right sweeps for each in the mesh; Direction Parallel (DP) which concurrently sweeps each row in all directions within an octant. Wall-clock speedup generally scales up with problem size but is very sensitive to machine loading during execution. On a lightly loaded Y/MP speedup factor of 5 was achieved using 8 tasks with a problem as small as ~ 104,000 cells.
6. RESTRICTIONS OR LIMITATIONS
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. Flexible dimensioning is used throughout so that no fixed limits on group, problem size, etc., are applicable. External forces and nonlinear physical effects cannot be treated.
7. TYPICAL RUNNING TIME
Central processor unit (CPU) time used by the flux sweep is roughly proportional to the number of flux calculations: spatial mesh cells * directions * energy groups * iterations/group. Times for running all test cases in the package range from about 10 minutes on an IBM RS/6000 Model 590 to several hours on slower machines.
8. COMPUTER HARDWARE REQUIREMENTS
DOORS3.2a modules run on IBM RS/6000, Sun, DEC Alpha and personal computers. The previous DOORS3.2 release ran on Cray, SGI, and Hewlett Packard. These computers are not available to test DOORS3.2a buy may still work. Reducing compiler optimization may be required.
9. COMPUTER SOFTWARE REQUIREMENTS
DOORS3.2a runs under Unix, Linux and Windows operating systems. Some numerical errors may occur when solving difficult problems on short-word configurations. Executables created with Portland Group, Inc. compilers are included both for Linux and Windows. All other systems require Fortran and C compilers. DOORS 3.2a was tested on the following systems:
· PC Windows XP with Portland Group Fortran 4.0‑2 and Portland Group C 4.0‑2 - requires Cygwin environment to build (see http://www.cygwin.com)
· AMD Athlon running RedHat Linux 7.3 with Portland Group, Inc. 4.0‑2 compiler
· AMD Athlon running RedHat Linux 7.3 with GNU Fortran 0.5.26 and GNU gcc 2.96
· PC Red Hat Linux 9.0 system with GNU Fortran 3.2.2 and GNU gcc 2.96/GNU gcc 3.2.2
· DEC Alpha ‑ Compaq Tru64 UNIX V5.0A ‑ Compaq Fortran Compiler V5.5‑1877‑48BBF
· DEC alpha ‑ Digital UNIX V40.D (Rev. 878) ‑ Compaq Fortran Compiler V5.5‑1877‑48BBF
· IBM RS/6000 ‑ AIX 5.1 system ‑ IBM XL Fortran for AIX Version 08.01.0000.0000
· IBM RS/6000 ‑ AIX 4.3.3 XLF 6.1 & XLC/C++ 4.4
· Sun SparcStation ‑ Sun OS 5.6 ‑ Sun WorkShop Compilers 5.0 98/12/15 FORTRAN 77 5.0
· Sun SparcStation ‑ Sun OS 5.7 ‑ Sun WorkShop 6 2000/04/07 FORTRAN 77 5.1
· UBUNTU Linux 7.04 with g77 [GNU Fortran (GCC) 3.4.6 (Ubuntu 3.4.6-5ubuntu1)]
·
The previous DOORS 3.2 release ran under the following systems that are no longer available for testing.
· Cray Unicos 8.0.3, cft77 6.0.4.10 and Cray J90 Unicos 9.2.0.03, f90 3.0.1.3
· SGI IP22, IRIX 5.3, F77 Ver 4.0.2
· HP 9000/715, HP‑UX 10.20, HP Fortran 77 Version B
****Added May 2007 **** DOORS 3.2a will build successfully and run test cases under ubuntu linux 7.04 using g77. Detailed instructions can be found in entry 26 of the DOORS electronic notebook on RSICC’s website. Only one small change to the mupopx11.lnx script is required to fix a problem with the X11 library.
Modify Install/mupopx11.lnx:
change line with X1
to
X1="-lX11"
The distributed Windows executables can be run in a command prompt window (of WindowsXP or Windows2000) in a manner similar to UNIX executables (uses redirection for input and output). The included jdos utililty is a c shell script intended for UNIX and UNIX‑like systems. It can be used under the Cygwin environment on Windows. Without the availability of the jdos script, each of the codes can be run separately, or in a batch file. See Doc/doors-3.2a-notes.txt for PC information. Installation instructions are in file Doc/install.doc.
10. REFERENCES
a. included in file C650.pdf:
W. A. Rhoades and D. B. Simpson, “The TORT Three-Dimensional Discrete Ordinates Neutron/Photon Transport Code,” ORNL/TM-13221 (October 1997).
M. B. Emmett, W. A. Rhoades, R. L. Childs, and J. O. Johnson, the DORT section of “A User’s Manual for MASH 1.0 - A Monte Carlo Adjoint Shielding Code System,” ORNL/TM-11778 (March 1992).
R. L. Childs, “FALSTF Last-Flight Computer Program,” ORNL/TM-12675 (Jan. 1996).
RSICC, “Input Modifications Required for ANISN and GBANISN” (July 1998).
W. W. Engle, Jr., “ANISN, A One-Dimensional Discrete Ordinates Transport Code with Anisotropic Scattering,” K-1693 (March 1967).
W. W. Engle, Jr., M. A. Boling, and B. W. Colston, “DTF-III, A One-Dimensional, Multigroup Neutron Transport Program,” NAA-SR-10951 (March 1966).
R. L. Childs and D. E. Cullen, “GROUPBAND-ANISN: A Code to Perform Multiband Calculations,” Informal Notes (June 1994).
W. A. Rhoades and M. B. Emmett, "DOS: The Discrete Ordinates System," ORNL/TM‑8362 (September 1982).
C. O. Slater, “The XTORID and ISPL3D Codes for Plotting TORT Activities,” ORNL Internal Memo (March 25, 1996).
R. L. Childs, "GRTUNCL: First Collision Source Program," ORNL Informal Notes (1982).
R. L. Childs, W. A. Rhoades, "Theoretical Basis of the Linear Nodal and Linear Characteristic Methods in the TORT Computer Code," ORNL/TM-12246 (January 1993).
W. A. Rhoades, "The TORSED Method for Construction of TORT Boundary Sources from External DORT Flux Files," ORNL/TM-12359 (August 1993).
W. A. Rhoades and D. B. Simpson, “Splicing and Bootstrapping Methods for DORT/TORT-to-TORT Coupling, (TORSET Section),” ORNL/TM-13350 (to be published).
W. A. Rhoades, “The ALC1 Program for Cross-Section Library Management,” ORNL/TM-4015 (December 1972).
b. background information (included in tar file):
S. L. Thompson, “CMP Code Maintenance Package,” SAND85-0825 Sandia National Laboratory Working Document (April 8, 1991).
S. L. Thompson, “The RSCORS Graphics System,” SAND99-XXXX Sandia National Laboratory Working Document (October 19, 1991).
S. L. Thompson, “A Definition of the Basic Graphics Package (BGP) Intermediate File Format,” Based on unpublished letter (March 1978).
11. CONTENTS OF CODE PACKAGE
Included are the referenced documents in PDF format and a GNU compressed Unix tar file, which is distributed on CD-ROM. The tar file contains source files for all programs in the auxiliary codes list, executables for Windows and Linux, implementation instructions, procedures, description of sample problem cases, test case input and output. WinZIP 8.0 or newer is required to expand the distribution file under Windows.
12. DATE OF ABSTRACT
June 1996, revised August 1996, February 1998, April 1998, July 1998, October 2003, May 2007.
KEYWORDS: DISCRETE ORDINATES; NEUTRON; GAMMA-RAY; MULTIGROUP; ADJOINT; SPHERICAL GEOMETRY; SLAB; CYLINDRICAL GEOMETRY; ONE-DIMENSION; TWO-DIMENSIONS; COMPLEX GEOMETRY; WORKSTATION.