11.4. ICE: Module to Mix Multigroup Cross Sections
N. M. Greene,^{*} L. M. Petrie, S. K. Fraley^{*} [ice1]
ABSTRACT
ICE is a legacy SCALE utility program that reads microscopic cross sections from an AMPX working library and uses input mixture number densities to produce macroscopic cross sections, which are written to an output file in the AMPX working library format. User input is entered with the FIDO procedures.
11.4.1. Introduction
ICE (Intermixed Cross Sections Effortlessly) is a legacy SCALE utility program that reads microscopic cross sections from an AMPX working library and uses input mixture number densities to produce macroscopic cross sections, which are written to an AMPX working library output file. The code was originally developed to allow efficient cross section mixing with minimum user effort and with reduced core storage requirements. The SCALE version of ICE is the latest in a series [UTILS-Fra76] of versions of the program.
In previous versions of SCALE several sequences employed the ICE module as a component in the self-shielding procedure; however in modern sequences the functionality of ICE has been replaced by new routines in the XSProc module. The ICE module is retained in the modern version of SCALE mainly for use as a standalone executable module to compute macroscopic cross sections and to provide backward compatibility with legacy inputs.
11.4.2. Cross Section Mixing Expressions
The mixing operations in ICE use the simple expressions presented below.
11.4.2.1. Cross-section mixing for AMPX libraries
For the options that produce AMPX working libraries, the mixing of cross sections involves a very simple summing of constituent values times a number density for the constituent, that is, \(\Sigma\), a macroscopic value, is determined by
where the j are the individual nuclides in the mixture whose number density and microscopic cross sections are N_{j} and \(\sigma\)_{j}, respectively.
The only exceptions to the above rule are for fissionable mixtures where the number of neutrons per fission, \(\nu\)_{g}, or a fission spectrum, \(\chi\)_{g}, is required:
\(\chi\)_{g} is defined as the fraction of the fission neutrons produced by the mixture which fall in group g. By definition,
ICE uses the following scheme to determine \(\chi\). First, terms F_{g} are determined by
where \(\text { VO } \mathrm{fg}, \mathrm{j}\) is the average of the product of \(\nu\) times \(\sigma0\)_{f} for the nuclide, \(\chi\)_{g,j} is the nuclide fission spectrum, and \(\hat{\phi}_{\mathrm{g}}\) is an estimate for the integrated flux in group g’. Once the F_{g} are determined, \(\chi\)_{g} is determined by normalizing the sum of F_{g} to unity.
In many AMPX libraries, the integrals of the spectrum used to determine the multigroup values are carried on the library for each nuclide. ICE uses this nuclide-dependent spectrum to determine \(\chi\)_{g}. This option should be exercised with caution, however, for no attempt is made to ensure that the individual spectra are consistently normalized.
11.4.3. Input Instructions
The input to ICE uses the FIDO schemes described in the FIDO chapter. In the descriptions, the number of entries expected in an array is given in square brackets.
********************************************************************************
Card A (20A4)
Title card
Data Block 1
-1$ Direct-Access Specifications [4]
NB8 No longer used.
NL8 No longer used.
NB9 No longer used.
NL9 No longer used.
0$ Logical Unit Specifications for Various Cross-Section Libraries [5]
INTAPE Input AMPX working library unit; default 4.
IOT1 Output AMPX working library unit; default 3.
IOT2 No longer used.
IOT3 No longer used.
IOT4 No longer used.
1$ Problem Size and Major Options [7]
MIX Number of cross-section mixtures to be made.
NMIX Number of mixing operations (elements times density operations) to be performed.
IFLAG(1) Set greater than ten if AMPX working library output desired
IFLAG(2) No longer used..
IFLAG(3) No longer used..
IFLAG(4) No longer used.
KOPT No longer used.
T - Terminate Block 1
Data Block 2
2$ [NMIX]
1. (KM(I),I=1,NMIX) Mixture numbers in the mixture specification table – values range from 1 to MIX.
3$ [NMIX]
1. (KE(I),I=1,NMIX) Element identifiers for the mixture specification table.
4* [NMIX]
(RHO(I),I=1,NMIX) Atom densities for the mixture specification table.
5$ [MIX]
1. (NCOEF(I),I=1,MIX) Number of Legendre coefficients, including P_{o}, to be mixed for each mixture.
6* [NG+4]
No longer used.
12$ [NMIX]
`1.(NUCMX(I),I=1,NMIX) Element mixture identifiers for the mixture specification table.
7$ No longer used.
T - Terminate Data Block 2
Data Block 3
8$ [MIX] Required only if IFLAG(1) > 0
1. (MID(I),I=1,MIX) Mixture ID numbers for AMPX working library; default (MID(I)=I,I=1,MIX)
9$ [N] No longer used.
.
.
10 No longer used.
11 No longer used.
T - Terminate Data Block 3
11.4.4. Sample Problem
A simple case has been selected to demonstrate the use of ICE. In this case, it is desired to produce mixture cross sections for UO_{2} and H_{2}O using basic data from ENDF version 7 238 group SCALE library. Information pertinent to the basic data is given in the following table:
Nuclide |
Identifier |
Order of Scattering |
||
---|---|---|---|---|
^{235}U |
92235 |
5 |
||
^{238}U |
92238 |
5 |
||
O |
8016 |
5 |
||
H |
1001 |
5 |
The atom densities to be used are:
UO_{2}
N(^{235}U) = 0.01 atoms/(barn-cm)
N(^{238}U) = 0.04 atoms/(barn-cm)
N(O) = 0.08 atoms/(barn-cm)
Water
N(H) = 0.06 atoms/(barn-cm)
N(O) = 0.03 atoms/(barn-cm)
In the sample case, we have elected to make an AMPX working library on logical 61,
We have selected further to identify UO_{2} with a 111 on the AMPX working library.
CSAS-MG PARM=CHECK is run to set up the master library, then WORKER is run to produce a working library for ICE.
A listing of the input follows:
=csas-mg parm=(check)
cross sections for ice sample problem
v7-238
read composition
atom 1 1 4 1001 1 8016 1 92235 1 92238 1 end atom
end composition
end
=ice
sample ice problem
0$$ 4 61 62 63 64
1$$ 2 5 13 13 13 13 2 1t
2$$ 3r1 2r2
3$$ 92235 92238 8016 1001 8016
4** 0.01 0.04 0.08 0.06 0.03
5$$ 1 2
12$$ f1
2t
8$$ 111 222 9$$ 1 2 3 11$$ 100 1111 2222 3t
End