Computing Applications Division Nuclear Engineering Applications Section Oak Ridge National Laboratory
Number 8 , January 1994
Click here for previous issues of the SCALE Newsletter
In This Issue ...
SCALE-4.2 Available from RSIC
The first version of SCALE designed for an IBM RS-6000 workstation has been tested and packaged by the Radiation Shielding Information Center (RSIC) at Oak Ridge National Laboratory (ORNL) and is now available for distribution.The modifications and enhancements included in Version 4.2 should substantially improve the portability of SCALE to all types of UNIX workstations. One master FORTRAN source that can be easily ported to Cray UNICOS, IBM MVS, or UNIX workstations is now maintained at ORNL under our quality assurance and configuration control guidelines. System-dependent routines for UNIX and MVS systems are similarly maintained.
Significant improvements in SCALE-4.2 are listed below. Most of these appeared in previous issues of the SCALE Newsletter, but are included here for completeness.
SAS4 Added an option to input axial source profiles for both radial and axial dose calculations. Also added an option to calculate point-detector estimates from collisions in both halves of the geometry for axial dose calculations.
ORIGEN-S, SAS2, COUPLE Modified programs to accept the new updated and expanded decay data and fission-product yield libraries (January 1993 Newsletter) and to use combined binary libraries that include multicycle cross sections.
ORIGEN-S LIBRARIES Replaced the six standard ORIGEN-S card image libraries with two new libraries, END6DEC and XSECTPHO. END6DEC contains the updated and expanded decay data library based on ENDF/B-VI data. XSECTPHO contains the basic cross-section and photon-spectra data and updated fission-product yield data based on ENDF/B-V data. In addition, the BASICLWR binary library has been replaced with a binary library that has updated cross sections from the ENDF/B-IV and -V 27BURNUPLIB for all nuclides in that library.
HEATING7, HTAS1, OCULAR Version 7.2 of HEATING, which operates on both UNIX and IBM/MVS operating systems, was added. HTAS1 and OCULAR were updated to run with this new version.
A program to convert the source files between IBM mainframe and workstations is included to facilitate code conversion. Language flags allow some of the source codes to compile on DEC/Ultrix and Cray/UNICOS systems. Suggestions are provided for Vax/VMS, Sun, HP, and ALPHA users. The files come configured for an IBM RS-6000 workstation.
SCALE-4.2 Minor Corrections
The following minor corrections were made to SCALE-4.1 to produce SCALE-4.2. Most of these have appeared in previous issues of the SCALE Newsletter, but are included here for completeness.
KENO V.a: Module was fixed to correctly calculate matrix keff by array number for a problem using a bare global array (i.e., there is no global unit).
KENO V.a: Corrections made for applying differential albedo boundary conditions to supergrouped problems where the global unit contains only an array specification. Note that this combination of options is rare. The effect of this error on keff is very large for this type of problem and will give physically unrealistic results.
CSAS and KENO V.a: Error checking during input processing was added so that these modules terminate with an error message if input errors are encountered.
MIPLIB: Corrected calculation of Dancoff correction factor for cylindrical and spherical lattice cells. ("See Words to the Wise..." below for more information.) Also corrected calculation of Dancoff factor for multiregion slab cell with vacuum boundary conditions to appropriately set the value to zero.
MIPLIB: Made a correction to accept an input Dancoff factor of zero and use it. Added a check to ensure that all Dancoff factors input are between zero and one.
MIPLIB: Corrected the dimensional consistency test for annular cells and asymmetric slab cells. This change allows correct problems that previously failed to run properly.
NITAWL-II: Parameter added to prevent exponent underflows for very dilute resonance calculations on workstation.
NITAWL-II: Corrected discrepancies in maximum fractional energy loss of neutron in admixed moderator calculation and Simpson rule calculation of collision density as documented in "Improved Calculation of Flux Shapes with the Resonance Shielding Code NITAWL," by J. Oppe, ECN-I--93-003. Slight effect on all calculations; impact will vary, but is insignificant for hydrogen-moderated systems.
MORSE and SAS4: (a) Fixed the option that allows the geometry media number to differ from the cross-section media number. Prior to this correction, MORSE did not select the proper cross- section media number when this option was selected. Previously run problems with geometry and cross- section media numbers that differ should be rerun with the new version. (b) Error recovery procedure was added so that when a fatal geometry error occurs, the results will be normalized and printed for n-1 batches, where n is the batch number in which the error occurred.
MORSE-SGC: Modifications made for compatibility with the new options in SAS4 (i.e., the axial source profile input option and the option for estimation to point detectors from collisions in both top and bottom halves of geometry).
Standard Composition Library: Updated resonance data flags for nuclides that have resonance data in ENDF/B-V.
XSDRNPM: Corrected calculation of number of direct access data blocks needed to weight the cross sections to prevent occasional failure. Improved calculation of balance tables.
SAS2: Corrected so that 'parm=check' would work on workstation.
27GROUP, 27BURNUP, and 218GROUP: These libraries have been updated to correct an error found in the chlorine cross sections.
BONAMI: Corrected so that a case with a number density of zero for a nuclide that has Bondarenko data will run without failing.
AIM: ft47ft001 is no longer required for AIM to execute on workstations.
All modules: The Program Verification Table, printed by Subroutine Library member LISTQA, was modified to print the correct version numbers for modules with version numbers of 10 or greater.
No KENO-V.a Training Class
Due to a reduction in funding from the Department of Energy, a KENO-V.a training course will not be offered this year. To have your name placed on a waiting list for a future training course for KENO- V.a or SCALE, please send a FAX to Steve Bowman at 615-576-3513 indicating which training course you would be interested in attending.
Words to the Wise...
Don't Let Your Holes Touch in KENO-V.a Geometry
If a KENO-V.a geometry model uses the HOLE option (Sect. F11.4.4 of the SCALE documentation, EXTENDED GEOMETRY DESCRIPTIONS), the user must ensure that the outer boundary of a unit used as a hole does not share a boundary (including points of tangency) with (1) another hole or (2) a region of the unit containing the hole. Either condition can produce erroneous results, either conservative or nonconservative, without triggering warning or error messages. In a recent analysis, tangent rows of fuel assemblies placed as holes produced a calculated keff value that was 0.7% lower than when a slight gap was placed between the rows of assemblies.This restriction is stated in Sect. F11.4.4 of the SCALE documentation near Figures F11.4.2 and F11.4.3. It applies to all versions of KENO-V.a and has been in the documentation since the SCALE-4.0 release. The SCALE staff at ORNL are currently investigating a programming change to solve this problem. In the meantime, users should make sure there is a slight gap in the KENO-V.a geometry model between each hole and any adjacent region. The gap can be as small as 1 in the fifth significant digit.
Error in Dancoff Factor Corrected
An programming error was recently discovered in MIPLIB in the calculation of the distance to the first untreated row (i.e., not treated explicitly in the analytic formula) of cylinders for use in the estimation of correction to the Dancoff factor for all the untreated rows. The code variable used as the pitch of the first row was one that had been changed to the pitch of the second row. This error caused previous calculations of small cylindrical cells (fuel O.D. < 0.4 cm) to provide slightly nonconservative k values, although the SCALE documentation and the American Nuclear Society Transactions, Vol. 56, p.324, indicated that slightly conservative results would be obtained. More details on the correction to the Dancoff factor and its potential impact on different lattice cell configurations is provided in the new SCALE-4.2 documentation in Sect. M7.2.This error has been corrected in SCALE-4.2 and SCALE-PC, as indicated on the previous page. All CSAS, SAS1, SAS2, SAS3, and SAS4 cases with cylindrical lattice cells will give slightly different results due to this correction. For consistency, the spherical case was modified to use the same variable as the cylindrical case, although the original variable had not been modified in this case. A second change that involved reducing the distance to the first untreated row by the diameter of the pin or sphere to get a surface-to-surface distance was not implemented after testing showed it to be overly conservative.
Trying to Stack SAS1 Input Doesn't Work
An error has been discovered when executing SAS1 within an analysis stream involving other SCALE modules. When a user stacks input to other SCALE modules following a SAS1 input file, the subsequent modules are not executed by the SCALE driver. SAS1 fails to set its termination flag properly, and the driver can not call the subsequent programs. This problem can be circumvented on a UNIX workstation by immediately calling the driver a second time with the input for the remaining modules. This problem has been corrected in the local version of SCALE-4.2 on the ORNL computers.
Error in SAS2H Using Multiple Libraries per Cycle
An error was discovered in SAS2H for cases using the multiple libraries per cycle option (i.e., specifying NLIB/CYC > 1). The document states that the boron content of the moderator is decreased during the cycle using a simple linear interpolation. Instead, the code actually uses the average boron content in the density data for producing all of the libraries in a cycle. This feature works correctly in SCALE-3. It is incorrect in versions 4.0 through 4.2 of SCALE. The problem may be circumvented by specifying NLIB/CYC=1 and increasing the number of cycles to the total number of libraries. Then, change BURN and BFRAC for each cycle to make the uptime and boron content what the code should have generated for NLIB/CYC > 1. For example, consider the following SAS2H input:NLIB/CYC=2 NCYCLES=2 BURN=400 BFRAC=1.0 BURN=300 BFRAC=0.8 The adjusted input to circumvent the error in SAS2H would look like this:
NLIB/CYC=1 NCYCLES=4 BURN=200 BFRAC=1.45 BURN=200 BFRAC=0.55 BURN=150 BFRAC=1.16 BURN=150 BFRAC=0.44
SCALE-PC NEWS
Update of SCALE-PC Available Through RSIC
An updated version of SCALE-PC with some minor corrections has been released to RSIC for distribution. This version is an update created on 11/15/93. It includes the following modifications:CSAS - Corrected the Dancoff correction factor calculation for cylindrical cells. The previous version gave nonconservative k values for small cylindrical cells (O.D. < 0.4 cm). (See Words to the Wise... for more information.) Also corrected the type of function E3 from single to double precision. This error caused some problems to fail with underflows or overflows. Corrected the Dancoff factor for a multiregion slab cell with vacuum boundaries to be set to zero. These errors were also in SCALE-4.1, but were corrected in SCALE-4.2.
KENO-V.a - Corrected the type of three variables from single to double precision that caused a problem with a very small number density (1E-20) to fail with an arithmetic overflow. Changed the default memory for a KENO-V.a stand-alone problem to 400,000 words. The default is 1,000,000 words in the mainframe/workstation version. Modified the program to write the keff by generation number to the PC screen as it is calculated so that the user can monitor the progress of the calculation.
All Modules - Relinked with new SUBLIB subroutine ALOCAT that dynamically allocates the requested memory for a module when it is called by the SCALE driver. Previous versions of SCALE-PC allocated 100,000 words of memory for all modules, regardless of the amount requested.
Notice to SCALE-PC Users with MS-DOS6
MS-DOS 6 made a change to in the way it accesses extended memory with EMM386. This change resulted in increased load times. We have observed that each functional module takes about 40 seconds longer and CSAS takes about 2 minutes longer to run on a 486/33MHz machine with DOS 6. Lahey has informed us that they have made a fix in Version 5.2 of their F77/EM32 Fortran compiler to handle this problem. We plan to compile the next release of SCALE-PC with this new compiler version.
For additional information contact: S. M. Bowman
Click here to return to the NEAS Home Page
Click here to return to the ORNL CPE Home Page
Click here to return to the ORNL Home Page
Last Modified Thu Feb 13 15:31:34 EST 1997