|Program name||Package id||Status||Status date|
|Package ID||Orig. computer||Test computer|
|CCC-0650/04||Linux-based PC,PC Windows,DEC ALPHA W.S.,IBM RISC6000 WS,SUN W.S.|
The DOORS3.2a discrete ordinates transport code system includes the most recent versions of CCC-0543/TORT-DORT, CCC-0254/ANISN-ORNL, CCC-0628/GBANISN and CCC-0351/FALSTF. It also includes the ISOPLOT code from the PSR-0155/DOGS package and various utility programs listed below 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 of 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 BOT3P2.0 and 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.
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.
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 relationship 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.
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 groups, problem size, etc., are applicable. External forces and nonlinear physical effects cannot be treated.
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 1o minutes on an IBM RS/6000 Model 590 to several hours on slower machines.
The code system contains the following modules:
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 Uncolided 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
|Package ID||Status date||Status|
CMP Code Maintenance Package, SAND85-0825 Sandia National Laboratory Working Document (April 8, 1991)
The RSCORS Graphics System, SAND99-XXXX Sandia National Laboratory Working Document (October 19, 1991)
A Definition of the Basic Graphics Package (BGP) Intermediate File Format, Based on unpublished letter (March 1978)
|Package ID||Computer language|
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 at RSICC 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 188.8.131.52 and Cray J90 Unicos 9.2.0.03, f90 184.108.40.206
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.
change line with X1
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.
WinZIP 8.0 or newer is required to expand the distribution file under Windows.
Keywords: complex geometry, cylindrical geometry, discrete ordinate method, gamma-ray, multigroup, neutron transport equation, one-dimensional, photon transport, shielding, slabs, spherical geometry, three-dimensional, two-dimensional.