CCC-0650 DOORS3.2.
DOORS, 1-,2-,3-dimensional discrete-ordinates system for deep-penetration neutron and photon transport
NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES
[ top ]
[ top ]
To submit a request, click below on the link of the version you wish to order.
Only liaison officers are authorised to submit online requests. Rules for requesters are
available here.
| Program name | Package id | Status | Status date |
|---|---|---|---|
| DOORS3.2A | CCC-0650/04 | Arrived | 24-JUL-2007 |
Machines used:
| Package ID | Orig. computer | Test computer |
|---|---|---|
| CCC-0650/04 | Linux-based PC,PC Windows,DEC ALPHA W.S.,IBM RISC6000 WS,SUN W.S. |
[ top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION
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 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.
[ top ]
4. 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 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.
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.
[ top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM
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.
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.
[ top ]
6. 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 1o minutes on an IBM RS/6000 Model 590 to several hours on slower machines.
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.
[ top ]
8. RELATED AND AUXILIARY PROGRAMS
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
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
[ top ]
10. REFERENCES
- 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)
- 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)
CCC-0650/04, included references:
- 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)
[ top ]
[ top ]
| Package ID | Computer language |
|---|---|
| CCC-0650/04 | C-LANGUAGE, FORTRAN |
[ top ]
13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED
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 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.
WinZIP 8.0 or newer is required to expand the distribution file under Windows.
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 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.
WinZIP 8.0 or newer is required to expand the distribution file under Windows.
[ top ]
[ top ]
CCC-0650/04
Doc - install.doc and other information files and reportsInstall - installation scripts
Outputs - developer's output for Cray and IBM
bin - executables for Windows and Linux
dort - dort, dos and falstf, etc. source and scripts
isoplot - isoplot source and test case
rscors - rscors source
sltutil - cmp and other Sandia National Lab. utilities
tort - tort source and scripts
[ top ]
- C. Static Design Studies
- J. Gamma Heating and Shield Design
- K. Reactor Systems Analysis
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.