UNCLASSIFIED AEEW-R 2133 WINFRITH United Kingdom Atomic Energy Authority The ‘1986’ WIMS Nuclear Data Library MJ Halsall CJ Taubman | NEA DATA BANK POCEIVED ¢ = *OMAGE IDENTIFICATH + NEA 1803/02 fos Reactor Physics Division AEE Winfrith September 1986‘© UNITED KINGDOM ATOMIC ENERGY AUTHORITY Enquiries about copyright and reproduction should be addressed to the ‘Winfrith Secretariat, Atomic Energy Establishment, Winfrith, DORCHESTER, Dorset, England.AEEW-R 2133 This work has been partially funded by the CEGB under the Thermal Reactor Agreement THE '1986' WIMS NUCLEAR DATA LIBRARY M 3 HALSALL C J TAUBMAN SUMMARY The WIMS nuclear data library recently used for most WIMS calculations at Winfrith is known as the '1981' library. An updated version to be known as the '1986' library has been derived from the '1981' library, JEF-1 data files, and a number of integral adjustments. It contains data in 69 energy groups for 129 nuclides including fuel, moderator and structural materials, fission products and miscellaneous detectors. This document gives a brief summary of the data, concentrating on 2200 m/sec cross sections, Maxwellian averages and resonance integrals, but also including other important data and references where possible. Reactor Physics Division AEE Winfrith September 1986 sc205. CONTENTS Summary INTRODUCTION THE 69-GROUP 11986’ WIMS LIBRARY Energy Group Structure Sources of Basic Nuclear Data Resonance Integrals Thermal Scattering Cross Sections at 2200 m/sec Transport Cross Sections Fission Spectrum Burn-up Data Fission Product Treatment 0 Goldstein-Cohen \-Values and Potential Seattering Cross-Section o, NUCLIDES FROM THE JEF-1 LIBRARY FILES 3.1 New Nuclides 3.2 Updated Nuclides AMENDMENTS TO EXISTING NUCLIDES Copper (3063.0) Fission Products and Pseudo Fission Product Uranium-238 (8238.0) Plutonium-239 (6239.0) Uranium-234° (234) DELETIONS AND REPLACEMENT NUCLIDES ACKNOWLEDGEMENT REFERENCES Tables Energy Boundaries for WIMS 69-Group Libraries 69-Group Weighting Spectra for Nuclides Derived from UKAEA Nuclear Data Files List of Nuclides with Details of Data Tabulation Resonance Tabulations 2200 m/sec and Maxwellian Averaged Cross-Sections and Resonance Integrals (69-Group) AEEW-R 2133 (ii) Oo annnwwnn Swowwm @ er 10 a 14 15 17 25 ar6. Fission Spectrum 7. Burn-up Chains 8. Energy Release per Fission Including Contribution from Capture 9. Fission Product Data in WIMS 69-Group Library 10. Goldstein-Cohen \-Values and Average Potential Scattering Cross-Sections 11, Group Dependent A-Values for Oxygen, Iron and Zirconium for 69-Group Library AEEW-R 2133 (uit) Page 32 33 35 36 37 38UNCLASSIFIED 1, INTRODUCTION The '1986' library has been derived from the '1981' library and JEF-1"data files. In addition five existing nuclides have been amended, and seven have been deleted, of which all but one were redundant versions of multi-version nuclides. This report describes the '1986' WIMS library and includes sections on the new, updated and deleted nuclides. 2. THE 69-GROUP_'1986' WIMS LIBRARY 2.1 ENERGY GROUP STRUCTURE The basic WIMS library contains 69 energy groups (defined in Table 1) of which the first 14 ‘fast’ groups have equal lethargy intervals of 0.5 in the range 10 MeV - 9.118 KeV. By using a suitable choice of average condensation spectrum for deriving the cross sections in these groups, it has been possible to demonstrate that the flexibility of the 69-group model provides adequate accuracy in the computation of both effective fast fission cross sections and diffusion coefficients. for thermal reactor applications it is not necessary to include a special resonance self-shielding treatment in this fast region, other than that which is inherently contained in the average condensation spectrum. Below 9.118 KeV there are 13 'resonance' energy groups down to 4 eV for which the isotopes with significant resonance behaviour have their effective resonance integrals tabulated in each group as a function of effective potential scattering cross section and temperature. For admixed scattering materials in the fuel, and also for cladding materiais, the libvary includes Ao, where A is the Goldstein-Cohen factor appropriate to an iktermediate resonance approximation. The group boundaries in the lower part of the resonance region are chosen so that the important resonances in the fissile and fertile isotopes are well centred in their energy group. Below 4 eV in the 'thermal' region there are 42 groups whose energy boundaries are chosen to ensure adequate definition of neutron captures in the Pu-239 resonance at 0.29 eV, and the Pu-240 resonance at 1.0 eV. The scattering and absorption cross sections in the thermal range are computed from an assumed Maxwellian weighting spectrum within each group. In the case of Pu-240 the tabulation is temperature dependent, and the within-group weighting spectrum was based on a self-shielded flux distribution calculated from a mixture of Pu-240 and a typical amount of scattering material. AEEW-R 2133 12.2 SOURCES OF BASIC NUCLEAR DATA Subject to exceptions aescribea in %.3 and 2.4 below ana in Section 3 (which aescribes the new nuclides obtained from the JEF-1 library), the cross sections are obtained primarily from the UKAEA Nuclear Data Library 11] which contains multipoint Gross sections, or their equivalent functional form with appropriate parametric values. Recommended DFN's have been used in conjection with the processing code GALAXY L2) to produce group average cross sections. The various spectra used for the in-group weighting in GALAXY are detailed in Table 2. Table 3 gives a list of nucliaes contained in the 69-group library and their origins. 2.3 RESONANCE INTEGRALS In the case of the principal resonant absorbers U-238, U-235 and Pu-239, the resonance integrals in the library have been caleulatea using point energy cross sections determined from the multi-level, multi-channel formulation [3] contained in the GENEX code L¥]. Both s and p-wave neutron contributions are ineludea with a statistical treatment for defining resonances in the unresolved region. These cross sections have been used in the SDR slowing down code 14) to obtain a shielded neutron enersy spectrum defined by some 120,000 energy points. Less onerous calculations involving the use of the ERIC intermediate resonance method 15) or the EPITHET Loy slowing down coe have usually providea the basis for the other resonant absorbers. Several sets of resonance integrals can be stored in a WIMSD library for each nuclice identifier. These are identified by appending a cecimal digit to the stanaard nuclide identifier. However, in the '1986' library there is only one set for each nuclide, and with only one exception (2232.1), the iaentifying digit is zero. Historically, the work of Kemshell L7 and 8) and Chawla 19] supported the use of a uniform reduction of the U-238 microscopic cross section of 0.1 barn in the resonance region. This was consistent with the results of an American stuay [10] recommenaing modification of ENDF/B U-238 resonance «ata. Consequently, an adjustea set of U-258 resonance integrals was providea in the library (3256.5). work by Halsall Li1j led to further adjustments and a library dataset referred to as 7238.0. Subsequent re-evaluation of the R1/100H series of experiments in DIMPLE showed that the changes’in the resonance integrals to produce the 7236.0 dataset were not justified. Consequently a new dataset has now been produced (8238.0) which retains the 7238 data combined with the Kemshell/Chawla resonance integrals. Table 3 shows which of the nucliaes have associated resonance integral tabulations and these are detailed in Table 4, For these nucliaes, the resonance absorption data used in a WIMS calculation are cerived from the resonance integral tabulations rather than from the cross section tabulations (unless the code wora NORESONANCE appears in the wIMS input).Table 5 includes absorption and fission resonance integrals for all nuclides in the energy range 0.55 eV - 2 MeV. The information contained in this Table has been derived directly from the group constants in the library and because the energy range for the integral does not coincide with group boundaries, there is some approximation in computing contributions from the ends of the range. 2.4 THERMAL SCATTERING In earlier WIMS libraries, data have been supplied for several models of each moderating material. In the '1986' library special thermal scattering cross sections are supplied only for hydrogen (Nelkin model, references 12 and 13), deuterium (Honeck, reference 14), and graphite (Butland model, reference 15). The process for deriving the 42 x 42 matrix of scattering cross sections has been to supply the scattering law [16] S(a,8) distribution as input to the PIXSE code [17] which then derives the required scattering matrices at specified temperatures. Linear interpolation in the library is used to evaluate moderator properties at any required temperature in a WIMS calculation. All other materials have had their scattering laws derived from a free gas model [18]. Table 3 includes the thermal scattering models and temperatures and also the free atom scattering cross section for all nuclides in the library. Table 5 includes Maxwellian averaged thermal transport (section 2.6) and absorption cross sections obtained directly from the library cross sections. 2.5 GROSS SECTIONS AT 2200 M/SEC Table 5 includes the 2200 m/sec absorption and fission cross sections for each library nuclide. These values have been obtained directly from the library’ tabulations. Where these data tabulations were derived from data files of the Winfritn Nuclear Data Library, the values have been checked against values calculated for the data files and very close agreement found. 2.6 TRANSPORT CROSS SECTIONS The transport corrections applied to WIMS group constants are derived from an examination of the spherical harmonics formulation of the steady-state Boltzmann transport equation. A scattering cross section with P, and P, terms is assumed and it is shown that with certain idealised assumptions about the current spectrum the equation can be reduced to the isotropic scattering form with ‘transport corrected’ total scatter and self-scatter cross sections. During the development of WIMS, the assumptions about current spectrum have varied, leading to different forms of the transport correction ABEW-R 2133 3(oF (414i) later a weighted column sum correction was applied to groups above 4 eV gig 1, 3 tr Zo 3 gr ge where J(E) is assumed to have a y form from 4 eV to 10 Mev. (itt) ad current is not very appropriate at high energies and the present practice is to weight the column sum correction for the top 14 groups (above 9.1 KeV) by current spectra derived from age calculations in light water, heavy water or graphite. At the pre-library stage P, and P, scattering data are available for every isotope but for most isotopes the P, data are only stored implicitly in the library in the form of transport corrections. It is important to note that early on in a WIMS cell-calculation the code sets the self-scatter cross section equal to (transport - absorption - outscatter). Thus for elements without explicit P, data the self-scatter cross section on the library may or may not be ‘corrected. For the principal moderator explicit P, matrices are attached to the end of the library in the form: 1 488" 3 %s, starting with all scatters from group 1, all scatters from group 2 ete. Only four P, matrices are permitted. These are one each for H, D, 0, C, and are for a single’ temperature (293°K) only. The B, leakage calculations in WIMSD that use these data should therefore be made only for systems at 293°K, ie, most high leakage experimental assemblies. For these elements the P, self-scatter must be uncorrected on the library to allow calculation of total cross section. Finally, group 3 for oxygen is a special case in which manual adjustments were made to account for the "oxygen window" [19]. AEEW-R 2133 42.7 FISSION SPECTRUM The fission spectrum, new to wIMS libraries, but used extensively for fast reactor work’ [20], is equivalent to a Maxwellian with an average temperature of about 1.37 MeV. The spectrum is aetailea in Table 6, and its adoption’is justifiea in reference 11. 2.8 BURN-UP DATA A section in the WIMS library defines the links in the burn-up chains. The isotopes B-10 (identifier 10), Gd-154, Ga-155, Ga-156; Gd-157, ana Ga-t56, and the hafnium isotopes ‘are burnable inaddition to the heavy fuel nuclides and fission proauets. within the limitations of a single capture and a single aecay proauct, complete flexibility in linkage is allowed, ana is controlled entirely by the library data. It is assumed for example that a capture in am-241 produces 0.12 atoms of Am-Z4zm and that the Am-242g branch of the chain can be neglected. The burn-up chains are listed in Table 7. (See section 3.1 for the novel linkage of U-237.) The energy release per fission for each of the fissile isotopes is given in Table 8. The value for Pu-240 is probably about 7% low (214 pout consequent errors in power proauction are completely negligible. 2.9 FISSION PRODUCT TREATMENT In the '1986' WIMS burn-up library, 35 fission products are represented explicitly whilst the remainaer are combinea in a single pseudo-fission proauct (nuclide 4902). It shoula be noted that I-135 (ID 6235) ana Pm-149 (ID 5149) have been aaded to the earlier representation to facilitate more accurate calculations of Xe-135 ana Sm-149 poisoning. Table 9 lists these nucliaes giving their yields per fission for each of the Principal fissile nuclides. The decay constants were those recommenced by walker (22) when the durn-up library was set up. The cross sections were largely compiled from resonance Parameters with an appropriate amount of 1/V absorber adaed to obtain an acceptable absorption cross sevtion at 2200 m/sec. Where resonance parameters were not available tne total resonance integral was divided in proportion to the lethargy widths of the resonance groups. Other elements have been treated as pure 1/V avsorbers.’ Amenaments to the fission Product data are described in section 4.2. The fission proauct chains in the library are given in Table 7. In arriving at these chains, isotopes with short half-lives were neglectea and portions of chains associated with neutron capture, for which the capture cross section is small, were also neglectea. AEEW-R 2133 5Also contained in the library is the pseudo-fission product representation from earlier libraries. This is of no practical significance to the average SD user because the bdurn-up tables point to the explicit fission products. However there is an option in LWRWIMS [23] to switch to the pseudo-fission product treatment if desired. (See the warning in section 4.2.) 2.10 GOLDSTEIN-COHEN A-VALUES AND POTENTIAL SCATTERING CROSS SECTION 9 Each nuclide in the WIMS library has a tabulation for each of the resonance groups which includes a A-value and a potential scattering cross section, o,. These are used in calculating the contribution , \o,, per atom of that nuclide to the effective background scattering, which in turn is needed in the calculation of the resonance cross sections. Most nuclides are assumed to have energy independent values which vary smoothly with atomic weight [24], but detailed calculations with the SRD code have been made to deduce group dependent values for oxygen, iron and zirconium [24]. Table 10 gives the A-values for each nuclide and Table 11° contains’ the group dependent A-values. Also included in Table 10 are average 0, values for each nuclide, obtained by weighting the group vBlues in the WIMS library with the group lethargy interval. In many cases the o value {s constant or almost constant over the groups. Thos® cases where there is a significant variation from the average value are clearly marked in the table. 3. NUCLIDES FROM THE JEF-1 LIBRARY FILES Cross-section data in the WIMS 69-group format were produced from the JEF-1 library files for fourteen nuclides, eight new and six to update existing nuclides. NJOY was used to process JEF-1 library files (in the ENDF/B-V format) to broadened point files at 900°K. These files were then processed through GROUPR to produce 69 group data at infinite dilution. The weighting spectrum used was up to 0.2 eV Maxwellian (T = 0.0253eV) 0.2 eV to 0.82085 MeV 1/E 0:82085 MeV'to 10 MeV Fission Spectrum (T = 1.3539 MeV) To translate the output from GROUPR into the WIMS library format, the Canadian codes WIMSR and WILMA were acquired from Chalk River National Laboratory and implemented on the CRAY computer at Harwell (26]. AEEW-R 2133The resulting data were then transferred to the ICL 2988 at Winfrith and added to the 69-group library. 3.1 NEW NUCLIDES Eight new nuclides have been added to the library. These are Titanium (natural) nuclide identifier 48 Molybdenum (natural) nuclide identifier 96 Gadolinium-154 nuclide identifier 2154 Gadolinium-156 nuclide identifier 2156 Gadolinium-158 nuclide identifier 2158 Uranium-237 nuclide identifier 927 Neptunium-237 nuclide identifier 937 Plutonium-238 nuclide identifier gus The burn-up chain was extended to include U*?’, Np??? and Pu2?* as follows: Pure? (6239) 8 (wy) Po ues uae ue? uate (2235) (236) (927) (8238) (n,¥) (n,¥) (n,2n) << The half lifes used were U2?” 6.75 days Put** 87.74 years (decay product, u?** neglected) To model (n,2n) production from U??*, it has been assumed that the ratio of (n,2n) to fission in U?** is problem independent. DEN 160 + fission spectrum average o, 5, = 0.01725 op’ = 0.2865 Effective fission yield of U2" tron ut « 294725. 0,06 ABEW-R 2133 7WIMSD does not permit fission products to be fissile and therefore the fission cross sections of U**? derived from the JEF-1 files were removed. 3.2 UPDATED NUCLIDES Six nuclides have had their data updated from JEF-1 data files. These are: Neptunium-239 new nuclide identifier 1939 (old identifier 939) Plutonium-241 new nuclide identifier 1241 (old identifier 241) Plutonium-242 new nuclide identifier 1242 (old identifier 242) Americium-241 new nuclide identifier 951 (old identifier 941) Americium-242 new nuclide identifier 952 (old identifier 942) Americium-243 new nuclide identifier 953 (old identifier 943) In addition the branching ratio of Am?*' to Am?** has been altered from 25% to 12%. Harwell evaluations [27] of differential nuclear ‘data for Am?*! show recommended percentages between 10.0#1.0 (thermal) to 15.8+2.0 (fast). 4, AMENDMENTS TO EXISTING NUCLIDES Forty nuclides have had their data amended. These are: Copper new identifier 3063 36 fission products new identifiers +4000 Uranium-238 new identifier 8238 Plutonium-239 new identifier 6239 Uranium-234 . 4.1 COPPER 3063.0 To take account of resonance shielding in copper, a resonance data file has been created. This was formed from the existing copper data (nuclide 2063) in all groups, ie consistent for all op, except group 20 which was formed from » | o-f 10 0.246 30 02385 100 0.635 300 0.991 1000 12438 - 10" 1.889 sal 1:968 The resonance integrals in this table were derived from EPITHET calculations using explicit values of resonance parameters. The basic nuclide data for 2063 have been retained except for group 20 which has the infinitely dilute value for Capsorption AEEW-R 2133 8(1,968) taken from the resonance data file, and a new nuclide 3063.0 has been created to replace 2063. 4.2 FISSION PRODUCTS AND PSEUDO FISSION PRODUCT Evidence has been accumulated to indicate that WIMS overpredicted fission product absorption by at least 10%: A reduction of 15% in all fission product capture cross sections, except Xe-135, was required in order to achieve the required reduction in absorptions. Also, all scatter cross-sections have been changed to 5 barns and a number of spurious small negative values have been changed to zero (including Xe-135). The amended nuclides have identification numbers derived from the old identification number plus 4000, eg Caesium-134 old 134, new 4134 Samarium-147 old 2147, new 6147 Also retained in the '1986' library are the seven old Pseudo-fission products. ‘These are 251, 252, 253, 254, 500, 105, 135. These data have had none of the above changes made to them and are therefore no longer an acceptable approximation to the Present explicit scheme. They have been retained to facilitate comparisons with earlier calculations using unadjusted fission product data, and also because there is at present a large uncertainty associated with the 15% adjustment. 4.3 URANIUM-238 8238.0 The adjusted set of U-238 data produced for the '1981' library, 7238.0, has subsequently been shown to have been based on & flawed evaluation of experiments in the Winfrith DIMPLE reactor. Consequently, a new nuclide, 8238.0, has been Produced which retains the nuclide data from 7238 and combines this with the resonance data from 3238.5. (Note that this is a departure from past policy of numbering’ resonance data, which Would have identified the new set as 7238.5.) 4.4 PLUTONTUM=239 6239.0 The '1981' recommended Pu-239 nuclide data, 5239.0 [11], has been shown to have inconsistencies (oe pansport in resonance groups less than 9 ). A new nuclide 6239.0 ‘absorption ~ “scatter has been produced from 5239.0 with the resonance transport cross sections produced by adding 0, ..,49, to the infinitely dilute absorption cross section in the resonance tabulation. AEEW-R 2133 94,5 URANIUM-234 234.0 Two errors were found in the U-234 cross-section data, in groups 20 and 27. These resonance groups should have their values taken from the infinitely dilute value in the resonance data file. The errors would only be significant when the code word NORESONANCE was used in the WIMS input. The old identifier has therefore been retained, despite the amendments. 5. DELETIONS AND REPLACEMENT NUCLIDES Sixteen nuclides have been deleted from the '1981' set of nuclides to produce the '1986' set, nine of which have a replacement nuclide. There is now only one version of each nuclide available on the library, this being the recommended one. The relevant nuclides are: Deleted '1986' Equivalent 16 6016 10056 2056 2063 3063.0 112 2113, 235} 1235, 2235.0 3238 7238} 8238.0 939 1939 3239} 5239} 6239.0 241 1241 242 1242 gut 951 942 952 9u3 953 6. ACKNOWLEDGEMENT The authors are grateful to Mr K Deadman and Mrs C Eaton for providing the WIMS data that were derived from the JEF files. 10REFERENCES ‘lo 10. mW, 12. STORY, J S et al Evaluation, storage and processing of nuclear data for reactor calculations. 3rd Geneva Conference on Peaceful Uses of Atomic Energy, P 168 (1964) BELL, V J et al A user's guide to GALAXY 3. AEEW-R 379 (1964) JAMES, M F Recommended formulae and formats for a resonance parameter library AEEW-R 621 (1968) BRISSENDEN, R J and DURSTON, C A user's guide to GENEX, SRD and related computer codes AEEW-R 622 (1968) SUMNER, MRS HM ERIC 2, A FORTRAN program to calculate resonance integrals and from them effective capture and fission cross-sections AEEW-R 323 (1964) DRIGGERS, F E A method of calculating neutron absorptions and flux spectra at epithermal energies AECL-1996 (1964) KEMSHELL, P B An analysis of heavy water moderated cluster lattices using the WIMS code AEEW-R 549 (1969) KEMSHELL, P B Some integral properties of nuclear data deduced from WIMS analysis of well thermalised uranium lattices AEEW-R 786 CHAWLA, R. An assessment of methods and data for Predicting integral properties for thermal reactor physics experiments AEEW-R 797 HARDING, R S, GAVIN, P H and HELLENS, R L Systematic study of thermal reactor benchmark experiments ANS Trans, 1969. Winter Meeting, vol 12, No 2, p 744 HALSALL, M J Recent adjustments to the WIMS nuclear data library AEEW-R 1492 (1982) NELKIN, MS Scattering of slow neutrons by water. Physics Rev., 119, 741 (1960) AEEW-R 2133 WW13. Ts 17. 18. 19. 20. eile 22. 23. 24. BUTLAND, A T D and OLIVER, S M The preparation of WIMS light and heavy water thermal scattering data. AEEW-R 950 (1970) HONECK, H An incoherent thermal scattering model for heavy water Trans. Am. Nucl. Soc. 5, 47 (1962) BUTLAND, A TD Thermal neutron scatter by reactor moderators Ph D Thesis, University of Batn (1972) EGELSTAFF, P A and POOLE, M J Experimental neutron thermalisation. International Series of Monographs in Nuclear Energy, Vol 106 (1969) MACDOUGALL, J D PIXSE AEEW-M 318 (1963) BECKURTS, K H and WIRTZ, K Neutron Physics Springer-Verlag (9164) VEALE, M Monte Carlo calculations of neutron age in some moderating media AEEW-M 833 (1968) ROWLANDS, J L et al Production and performance of the adjusted cross section set FGL5. Proc International Symposium on Physics of Fast Reactors Tokyo, Dec 16-19, 1973 JAMES, M. Energy release in fission. Journal of Nuclear Energy 23 517 (1969) and 25 513 (1971) WALKER, W Yields and effective cross-sections of fission products and pseudo-fission products AECL 1054 (1960) HALSALL, M J LWRWIMS, a computer code for light water reactor lattice calculations AEEW-R 1498 (1982) ALDOUS, AC Numerical studies of the hydrogen equivalent of some structural materials in their effect on U-238 resonance capture AEEW-M 860 (1969) AEEW-R 2133 1225. HUGHES, D J et al Neutron cross sections BNL 325 26. SMITH, R and DEAN C J JEF-1 To be published 27. LYNN, J E et al Evaluation of differential nuclear data for americium isotopes. Part 1: ?*? Am AERE-R 8528 (1979) AEEW-R 2133 13BusRoy souwoas! FoR WIMS §9-GROUP LIBRAnIES oro eneray ceeray | teemarer [oe [tera feces we # 6.0855 aigsas | o.nsp97 | 28 o.1se3T 3.873 aije6s | ostoees | Ola3eet 2231 itis] o'sc019 | 30 2138 33 oes | orsoors | o:336Hr get Sse | oressse | o.tigio 8500, on3at sass | is pare) > 83005 | otters 50253 | de ocaeste = ort65 elias | ocsceso | 3 raete = ouitte orore | olasss | 36 0102389 2 blosnge | lotsee | ovsssrs | ar sroeess 2 Growoas | taaseg | olusses | 8 areata = oleaire | ocoreor | ovsser | 8 or02361 = Gietsos | olovg7s | ovagss9 | to o102e39 12] bioteds = afecsris | osoassiz | 9.49960 | wt perry \ 2 olosee 3 928821 Po a oc oas94 5 c22ise 1s Janec = $530.9 3508.0 0.50008 | 86 olazta 1 | S300 > i981 olasigs | a7 et it | isn > Zaggis oasis | 8 0.15353 ta | daslis > tas oasis | 89 808381 19 | twas > ‘306.98 onser | 50 cosase ne ee o'50395 | 51 2.08899 | 3eroe > tres trgoass | 5 0.11333 a | falas = ‘tigate ccerist | 33 onze B [ipo = tose auasier | 5 olaoost HR | fc! aniteo a oczsi3t & | aio > (15.968 sissees | 56 o:336tr | sg > Serr otéoss | 57 oiaagte a | ‘Serr = oa asm | ours 5 oc1has % orsaese & alias & oribese 5 olistis & 318232 & oczaais ra 0.28758 7 ocaope7 a |} Stes5is & > AEEW-R 2133 4TABLE 2 69-GROUP WEIGHTING SI TRA FOR NUCLIDES DERIVED FROM UKAEA NUCLEAR DATA FILES Dry wet Graphite Flux | Graphite Current Energy Spectrum Spectrum Spectrum Spectrum 11.3315. Mev | 000077 0003 S e 882497 3000288 10012 0.0004 0.0026 6.87289 :001082 0043, 0.00179 0.0126 5.35261 1003415 20125 0:00467 0.0225 4316862 2007253 0250 0.00897 0.0263 3524652 1015299 20477 0.0137 0.0323 2.52839 £032235 £0821 0:0332 0.0945 7.96912 040766 20925 0.0584 021610 1.53355 068804 £1270 0:0786 0.1919 1219433 092740 1446 0:1017 0.2156 0.930145 5151457 21831 0.1326 0.2534 0.724397 +293853 32321 0.1646 0: 2803 0.564161 +359617 12432 0.2029 0.3128 0.439639 :273385 12287 0.2496 03535 0.342181 3579860 +3326 0.3067 0.4060 0:266491 +797033 £3665 0.3805 0.4739 0.207543 ¥6020714 +4096 0.4728 0.5617 0.161635 15108708 54687 0.5907 0.6757 03125881 13392871 25421 0.7416 0.8237 980635 Kev | 1.538518 16346 0:9324 1.0097 76.3509 21025529 £1520 151815, 1.2587 59:4621 2.256529 £8839 154964 1.5689 46.3092 2.717019 1.0847 1.8971 1.9608 36 :0656 2:929716 152975 2.4161 224744 28:0879 2.913247 1.6002 350825 3.1337 21.8749 5: 728346 2:0012 3:9372 3:9800 17.0362 5.345382 2.4685 5.0363 5.0700 13:2678 6.613002 3.1004 4502 6.4716 10.3330 8.260331 358690 +2603 8.2603 8.04733 8.314061 4.8358 - - The DRY ‘spectrum was the original weighting spectrum used with GALAXY to produce group averaged data, for the 14 fast groups for all nuclides derived from the UKAEA data files except the following: 3001 16 for which the WET spectrum was used, and AEEW-R 2133 52212 3063 «181 «93 8238 for which the two GRAPHITE spectra were used. The current-spectrum is used for averaging the P? cross-sections In groups 15 to 56 inclusive the weighting is proportional to 1 ge "E/KT kT where kT = 2.522 x 107" MeV (NOTE: For exceptional treatment of Pu-240 refer to section 2.1). AEEW-R 2133 16TABLE 3 LIST OF NUCLIDES WITH DETAILS OF DATA TABULATION Nuclide Thermal Scattering Identification Nuclide Basic Resonance Number DEN Model Temperatures °k | Tabulation 3001 Hydrogen 10/D | Nelkin 293, 333, 373, 423, - 473, 560, 620 3002 Deuterium 905 | Honeck 293, 333, 373, 423, - 473, 550, 620 3 Helium-3 220 | cas 300 - 4 Helium-4 221 Gas 300 - 6 Lithium-6 214 | cas 300 - 7 Lithium-7 215 | Gas 300 - 9 Beryllium-9 50 | Gas 300 - 10 Boron~10 13/D | Gas 300 - 1010 Boron-10 (unburnable) | 13/D | Gas 300 - 1011 Natural Boron 15/E | Gas 300 - (purely 1/v < 100 KeV) 2012(1) [Carbon 68/D | Butiand 14 temp (1) - 2212 Graphite Damage - - - 5 17 Mitrogen 147A | Gas 300 - 6016(2) | Oxygen aT Butland ET | 9 temp (2) - 19 Fluorine 23a | Gas 300 - 23 Sodium 182/D | Gas 300 - 24(3) |Magnesium-24 gut Gas 300, 600, 900 - at Aluminium 35/0 | Gas 293 - 29 Silicon 25/a | Gas 300 - AEEW-R 2133TABLE 3 (Continued) Nuclide Thermal Scattering { Identification Nuclide Basic Resonance Number DEN Model. Temperatures °k | Tabulation 48 Titanium ger | - 900 a 51 Vanadium ge | - = - 52 Chromium 1a | Gas 293 = 1054 **Fe(n,p)°*Mn - I I 2 55 Manganese - - - - 2056 Iron out Gas 300 - 9056 (4) |Stainless Steel = - - 5 58 Nickel 46 cas 293 a 1058 S°Ni(n,p)**Co - - - - 3058 SeFe(n,y) Fe - a - - 8058(5) Inconel 750 e - - - 9058(6) | Inconel 718 - S - - 1059 *9Co(n,¥)*9Co C - - e 3063 Copper 3 Gas 300 3063.0 eu Keypton-84 - - - - a Zirconium ga | Gas 300, 600, 1200 - 1091 Zircaloy Damage = - . - = 93 Niobium 19 Gas 300 S 96 Molybdenum ver | - 900 a 2103 teeRh(nynt) !9mRh = - - I 3109 Silver 5 a = 3109.0 ABEW-R 2133 18TABLE 3 (Continued) Nuclide Thermal Scattering Identification Nuoliae Basic Resonance Number DEN Model Temperatures °k | Tabulation 2113(7) [Natural Cadmium - - - - 2115(8) | Indium - - 5 5 121 Antimony=121 - - - 121.0 123 JAntimony~123 . - 7 123.0 2154 Gadolinium-154 ger | - 900 a 2155(9) | Gado.inium-155 . a a = 2156 Gadolinium-156 ver | - 900 - 2157(10) | Gadolinium-157 7 7 a a 2158 Gadolinium-158 ver | - 909 Fe 164 Dysprosium-164 = - - 167 Erbium - - - - 176 Lutecium-176 - - - - 2176 Hafnium-176 1202 | - - 2176.0 auqT Hafnium-177 1203 | - a 2177.0 2178 Hafnium-178 roa | - - 2178.0 2179 Hafnium-179 1205 | - a 2179.0 2180 Hafnium-280 1206 | - a - 178 Natural Hafnium e - - - 181 Tantalum 328 | Gas 300 - 183 Tungsten 213. | Gas 300 - 197 Gold - - - - AEEW-R 2733 19LE 3 (Continued) Nuelide Thermal Scattering Identification Nuclide Basio Resonance Number DEN Model Temperatures °« | Tabulation ast}aat) - E e 7 Pseudo asa} can) - - e - } Fission 253) (11) - - e - Products asu}crt) - - e - 105(11) |Rnodium-105 - - 7 E 135(11) | Xenon-135 - - . - 500(11) |Samarium pseudo-fp - - . . 4083 Krypton-83 - - - . 4095 Molybdenum-95 - - . - 4099 Technet iun-99 - - . e 4101 Ruthenium-101 - - = - 5103 Rutheniun-103 - - - e 4103 Rhodium-103 - - - - 4105 Rhodium-105 - - S : 5105 Palladium-105 - - e - 4108 Palladium-108 - - : - 4109 Silver-109 - - a - ann cadmiua-113 - - eI - 4nn5 Indium-115 - - . : 427 Todine-127 - - . E 4131 Xenon=131 - - . : 4133 Caesium-133 = - - - SK G ye? ABEW-R 2133 20LE 3 (Continued) Wuclide Thermal Scattering Identification Nuclide Basic Resonance Nunber DFW Model Temperatures *k {Tabulation 6135 Todine-135 - . = - 4135 Xenon-135 - . - a 5135 Caesium-135 - - 2 - 4143, Neodymium-143 - - : - unas, Neodymium-145 - - - e aver Promethium-187 - - - _ StaT Promethium-147 - - - 5 eur Samar iun-147 - - . = 4148 Promethium-148a - - = = 5148 Promethium-148 - - S - 5149 Promethiun-149 - . - = ara Samar ium-149 - - . - 4150 Samar iun-150 - - - : ans Samar {un-151 - - - - ais2 Samar iun-152 - - - - 4153 Europiun-153 - - = e asa Europiun-154 - S 5 - 4155 Europium-155 - - . - 4197 Gadolinium-157 - - . - 4902 Pseudo - - S : 207 Lead 26 cas 300 - 2232 Thorium-232 2e/a | Gas 300 2232.1 AEEW-R 2133 21TABLE 3 (Continued) Nuclide Thermal Scattering Identification Nuclide Basic Resonance Number DEN Model Temperatures *k | Tabulation 9233 Uranium-233 202/0 | Gas 300 - 1233 Protactinium-233 > | Gas 300 - 234 Uranium-234 1TH | Gas 300 234.0 2235 Uranium-235 48 cas 293 2235.0 236 Uraniun-236 173, | Gas 300 236.0 ger Uranium-237 JEF = 900 - 8238 Uraniun-238 401/50] Gas 293 8238.0 937 Neptunium-237 JEF = 900 - 1939 Neptunium-239 JEF 900 - 948 Plutonium-238 JEF = 900 - - 6239 Plutonium-239 329 | Gas 300 6239.0 120 Plutonium-240 201 Gas 300, 600, 1200 - 1241 Plutonium-2ut BE - 900 - 1242 Plutonium-242 EF = 900 - 951 Amer icium-241 JE a 900 - 952 Amer icium-242m Jer - 900 - 983 Amer icium-243 EF - 300 - 1000 + 1/¥ absorber - 5 - - 1003 235 9 (2235) = = - t 199 Resonance part of 1/v - - - - 2000 = 1/¥ absorber - - - - 3000 Pure absorber = - - - 4000 Inverse Lethargy - a a Intervals ABEW-R 2133 22Notes to Table 3 Carbon, 2012, which combines epithermal data from DFN 68/D with a Butland improved thermal scattering model, has thermal data tabulated at 14 temperatures 293, 313, 400, 500, 600, 655, 700, 900, 1100, 1300, 1600, 2000, 2500, 3273%k Oxygen, 6016 which combines epithermal data from DEN 37 with a Butland extended effective temperature thermal scattering model, has thermal data tabulated at 9 temperatures 293.6, 450, 600, 750, 900, 1050, 1200, 1350, 1500°K Magnesium, 24, combines epithermal data from DEN 941 with thermal scattering data at 300, 600, 900°K obtained from PIXSE calculations on a free gas scattering model and a free atom scattering cross-section of 3.42 barns. Stainless steel 304, 9056, was formed from the following weight composition: 56(49.1%), 1056(21.0%), 58(10.0%), 52(18.5%), 55(0.8%), 29(0.6%) where 56 and 1056 are old versions of iron with different resonance integrals and 2056 on the present library is approximately equivalent to 56 and 1056 mixed in this proportion of 7:3. Inconel 750, 8058, was formed from the following weight compositions: 58(75.993%), 56(4.841%), 1056(2.075%), 27(1.1854) 52(14.819%), 29(0.395%)), 55(0.692%) Inconel 718, 9058, was formed from the following weight compositions: 58(59.4%), 52( 20.3%), 56(14.2%), 1056(6.1%) Data for 2113 are Cd-113 fission product data x 0.1226 plus "reasonable" fast data [25]. 2113 is not labelled a fission product. Data for 2115 are adjusted In-115 fission product data. In-115 is 96% of natural. 2115 is not labelled as a fission product. 2155 is basically Gd-155 fission product data with "reasonable" values included for og and fast groups [25]. 2155 is not labelled as a fission product and will burn. AEEW-R 2133 2310. WW. 2157 is basically Gd-157 fission product data with "reasonable" values included for o, and fast groups [25]. 2157 is not labelled as a fission product and will burn. These fission products are the original data, as on the "1981" library, without the 15% reduction applied to all the fission product capture cross-sections on this '1986' library. They will not appear in the normal burnup chains. AEEW-R 2133 24€€12 u-Mgay sz NB: All tabulations are in the form NIN 3063.0 3109.0 121.0 123.0 2176.0 2177.0 2178.0 2179.0 2232.1 Temp *k 300 300, 1200 173 173 300 300 300 300 300, 750, 1300 TABLE 4 RESONANCE TABULATIONS Barns 1.0 x 10! 3.0 x 10° 1.0 iOpen 10 infinite dilution 3.2 x 10! 6.2 x 108 1.02 6:02 x 107 2:002 x 10° 1:0002 1.0 x 108 1.0 x 108 1.0 1.0 x 108 1.0 x 107 1.0 5.0 x 10° 5.0 x 10? 5.0 infinite dilution Boe 5.0 x 10? 5.0 infinite dilution 5.0 x 108 5.0 x 107 5.0 infinite dilution 5.0 x 10° 5.0 x 107 5.0 infinite dilution pecan xel0x 1,02 x 107 1.62 15662 x 107 3208 x 10° 1:05 xx where t is the group lethargy width. 10% 102 10" 10? 10? 10? 10? 10° 107 10% 3.0 1.62 1:0 x 10? x 108 x 10* 02 x 10% infinite dilution infinite dilution infinite dilution 5.0 5.0 x 10" x 10° x 10° x 10° x 10? x 108 5.0 x 10° 5.0 x 10° 5.0 x 104 8.42 x 107€€12 u-magy 92 8238.0 6239.0 Temp _*k 300 300 300 300, 600, 900 300, 900 TABLE 4 (continued) p_Barns 1,0 x 107 5.0 x 10° 1.0 x infinite dilution 4.969 x 107 1.007 x 10% 1,701 x 8.356 x 10° infinite dilution 5.0 x 10? le Ope exe) 0c 5.0 x infinite dilution 1.553 x 10! 3.149 x 108 5.34 x 2.613 x 107 if Clete x10. 3:6 x 1.0 x 10% 3.5 x 10? 1.0 x 5.305 x 10? 6:48 x 10° infinite dilution 10% 10% 10" 108 10% 102 2,089 1.0 6.534 1:0 1,563 10* 10% 108 10! 10* 10? 5.0 x 10° 4.678 x 10? 1.0 x 10° 1.462 x 107 infinite dilution 3.136 x 107€€i2 u-Ngav lz wean Averaged Trea SECIS tte 20 we Cra sectlon ecnanoe 1tegrat enti ticatien Bel 2m (0S 446) 2.0 et = 0 (CS 28-69) Naber ecrptien | Tarepot | seater | risstn | toarptin | taeprt | scatter | Fission | moaxption isslan 32 c ot = on Ey : = ‘5 an | “si : = - ssa ox ‘| = : - 00 wm} cm | = : - 353 cy wr | = : : 08 nao | nas : - - ‘080 oa | ow - - - 0 aa 38 : : = 33 38 - e : 0 1 io - : - 090 sa | Ror : = 3 - oo oo Ei : : 2 : 1.88 16 30 = = 8 - 09 ap | - - 73 : ‘ovo yo | ie : : 3 = 36 sa | aos - : sk : ‘oon, we | oe : : 6 = 2 st 138 : : 159 - neo zp | an : : 233 : ew we De : : oie : 5.05 5.05 = : : 7 - 3310 ror 3 : : a : oo oo - = : oo = 88 38 - - = Re 2 -€€l2 u-maay 2 LE 5 (Cotirass) ‘meta Fexnarce Tegra seliian Averaged Tersl SECTS ertittation Bat Bho (OS 46) 10 ef = 0 (GS 29) inter Fister Wrarpere | seater | Fission | merption ‘seater | Fission 26 > zr 13 = zat Wo ~ 8 - % 13 - 251 103 - 2 = a 2 : ve m3 = a = one |= : 00 2 : ye - 1% = - 105 s : ae : 1% 13 : x6 mt : 8 : x m 6 - zn 132 : cy % > 3 39 - : Bo < - xe x 10 96 12 1 : 32 re |= ot x oom |e aa we 105 : an to = 31 ta | ee | = or wa 10 : ne ow |= et oo - = sea 1 - rr) - 3 om | su | = ny 1 105, : he = % sa | sar |= 2 7 - 23 = a oo ” : é = oa : So Bw: : pa : a5 = an zo : Waa 2 - amp : 21 28 : 0 02 : 183 = ct 9.23 : 22 as - 3. - i Tot : ™ as : 5 25k Be = 6 oa : Z 255 coz : tt ie : : 26 1.8 - 104 2a e : 2s zn : Eo i : - 28 6.01 : a 182 : : 16 zn : 2 00 : : wi 1% = = " e a .HELE 5 (Content) €€l2 y-magy wate a0 wa tr Stn cee age Tian neve tr aD statin Boren cae he) Na ero tas 0) oo me | am a5 alee @ iss ea res mela ae a ao f 2] 2 & ia 2]: tes 2 |: amr x B Df | a ae Df of | gs |: ane en ea | 2 | 2 | be Ee iE Se oa A ba er | 2 ft Eo ro ele as Bie EB a ao | 2 |: 2 me ff dt is |e "2 we ‘s yo | owe we lee 23 Sie m ‘De er | sa |: m8 a | : 3 fon | = ‘3 a ay | tr] 2 | ge & | ie : a so} o 23 ee)“ ] ol | ge m [| 2 m3 ae 2 a an aie a e Df ol] oe Bree a 35 BPs ie a 3 2 fe 7 ft s a i ee x oe : 2] t 3 ‘3 3 Do] of] gs 2s : > |: | wo wa Se | oe ae : eis 8 ‘| aero | aaro Dy | ‘en : E so | tar ne) Se 5 ae relive | es ca a we |e ry me : so |: melas |e "e9 mee eeea |e Bs ie : eae zs] s 83 ei | te] t . a Bale bm |: si] ote 383 to | ve | 2 oD |e |: 3 E298] ots 3 gh] Me |: a3 w fas | t em | my cer] ae is 3 te | t e m= |e |: sm | wre gee| te | owe tip we] f | ane ae | bs |: so | > by ad 3105 82 30 496 - 93 1 BS a 5.00, 7 Be 28h 08 9.28 we 4.96 - vr 6 BS : 5.00 7 sen | oS 83 es |e | 2 |e wm |b |” be |€€12 y-MasV walle "200 w/ace Crasrsxtien " ein Averaged Taal XE entitictien so 2 na (03 6) 1.) et 0 (G3 a9) abe Tmerrpticn Jane [aoa | marption | trampoe [sant [risen | tomrpin | trwewrt_| sete | Hsin pax [ans est «|: a BS = | we = pao] a me we |e 216 BS : % : ye) Mr 5.95 ke | > 10 BS : 5.2 = gto] mm B3 we | > 15 BS : 63 = Bae | MB aA % ) > © BS = aq = re a a we | > 1a Bs : a : os soe age |e 5 | BS : ‘ : roe | up | aren age | > our BS - - a 5135 7.80 awe |= 3 Bs = 656 : oa ay at ase |e 13 BS : 26 = ae at aa 4 = 25 BS : DA = det | cme 1 ae | > 29 Bs - 6 : 7) tow 1 a | = 8 BS : 1% : cud “eur a8 age |= BS : 6S = gyri} me za ke | = ae BS 2 | me > ‘Bs06] ae 12 4% 5085 185 - 13 : isbe| 190 ae |= x8 1 Bs - 105 : syn) 8 aie tse |S 2s ae BS 2] ame > 534) Ne 5 ke | > 21 Ea BS : 16.1 = nt 56 ne |= wit 3 BS = saw - ue "15 a 25 10 BS - wt - 53 a awe |: 136 vat BS - = : st 125 ae | = oi 3 BS : 13 - 3 ng ae |= 5108 sek BS = | tose 5. - ust ae | > a Ma BS =| ses x 5 we ws |: 31.9 wr BS : 15k zx : a 09 : 8 8 1s : a2 mi {mo fee€l2 u-mgay aaa Dame abn EE 5 (continu! elie enarce Integral MenelLan Averaged Toren (SECT eottFLatien 5 of = 2m (2S 6) A ero (a3 209) waptin | trwmpat | Sowte | Fastin | mxption | tramrt | sotue | Flasien ze Bs a 20 oml 6.53 wo | 1s 93 an wom 1% | sa sa tot 233 sat oH in 1a | “no wo | no att oe eI a vm] ina at | a ase me oT i | ze 2 a 30 aie 33 a 18 sit} 5a 15, tot . sar mn ‘ou 193 : x 5 a |= ese Fe a5 vo | 1 2M we | sis | 00 om % ca 18 2 1 1 We “ors 9B 43, on 14 zo | zo gs owe | 98 1 Eo ay | oa | an Ps ws | we esr 318 ie ir | oe 5 ca ws | ae 20 15 13 m aa | 2 zo 13 ome| 1a 168 2 3 1283 133 90 | 953 ee na st 26 | 16 ao | ax |” ors Bt 185 wor | 5 ser we | 2s oe vga vsi5 816 au or era ra aio 2a | Ba Ba tz | oms| 100 ) + ‘Ber on | ot = 1003 an : : Ey a - > 9 as {= : 00. 00. : : 200 sw] : = wr | = : 50 Ba - - 10 vo : ‘exo eh : - xa ve |= 5TABLE 6 FISSION SPECTRUM 69-Group library Group Spectrum 1 0.027100 2 0:116000 3 0:214800 4 0229300 5 03175800 6 0.110700 7 0:062100 8 0.032500 9 0.016300 10 0008000 W 0.003900 12 0:001900 13 0:000900 4 0.000400 15 0.000200 16 0000100 17 0:0 18 0.0 19 0:0 20 0:0 21 0:0 22 0.0 23 0.0 2 0:0 25 0:0 Total . |. 1.0 The fission spectrum is the one used generally in fast reactor work. It lies between the old Bonner spectrum, 1.3 MeV, and the later 1.43 MeV spectrum. AEEW-R 2133 32TABLE 7 BURN-UP_CHAINS Initial Nuclide | Capture Product | Half Life Decay Product 10 2154 CiEI 2155 2156 2156 2157 2157 2158 2158 2176 2177 2177 2178 2178 2179 2179 2180 2180 251 252 252 253 254 105 135 500 4083 4095 4099 4101 5103 39.9 days 4103 4103 4105 35.0 hours 5105 5105 4108 4109 4109 4113 4115 4127 4331 4133 4134 4134 A135 6135 4135 5135 yh 4145 4a47 4148 S147 5148 e147 4448 5149 42.0 days 5148 5149 5:3 days +1 years hours 4135 1 hours 5135 65 years 614a7 years 6147 AEEW-R 2133 33TABLE 7 (Continued) Initial Nuclide | Capture Product | Half Life | Decay Product 5149 2.3 days 4149 drug 4150 4150 4151 4151 4152 4152 4153, 4153 4154 4154 4155 4155 4157 4902 2232 1233 9233 234 1233 234 27.4 days 9233 234 2235 2235 236 236 927 927 8238 6.75 days 937 8238 1939 937 9he 1939 1240 2.35 days 6239 948 6239 87:74 years 6239 12k 12h0 12m 1241 1242 14.6 years 951 y2u2 953 951 952% 952 953 953 * Note: only 0.12 atoms of 952 are produced per 951 capture. AEEW-R 2133 34TABLE 8 ENERGY RELEASE PER FISSION INCLUDING CONTRIBUTION FROM CAPTURE Energy Release per Fission Isotope WIMS Nuclide Mev Joules/mole x 107" Th-232 2232 189.6 0.183 x 10718 u-233 9233 197.9 0.191 x 10788 u-234 234 189.6 83. x 10788 u-235 2235 201.7 0.1946 x 10788 u-238 8238 205.0 0.1978 x 10788 Pu-239 6239 209.9 0.2026 x 107" Pu-240 1240 196.9 0.19 x 1078 Pu-241 1241 212.8 0.2054 x 10788 Am-242 952 207.3 0.2 x 10718 Note: The energy releases from nuclides with low fission cross section in thermal reactors are known to be significantly in error (21). AEEW-R 2133 35TABLE 9 FISSION PRODUCT DATA IN WIMS 69-GROUP LIBRARY WIMS Fission Yield Isotope | ID‘ Number ve33 u235 u238 Pu2s9 Pu24o | Puewt Kr83 4083 0.0114 ]0.0055 {0.0009 |0.0029 |0.0025 |0.0021 Mo95 4095 0:0602 |0:0655 |0.0560 |0:0500 |0:0%90 }0:0460 Teg 4099 0:0480 |0.0025 0.0030 |o-:0610 |0:0595 |0-0580 Ru101 4101 0:0264 }0:0500 }0:0030 |0:0590 |0:0600 j0-:0610 Ru103 5103 020166 |0.0290 }0-:0620 |0:0560 |0:0595 |0:0630 Rh1035 4103 Rn105 4105 0.0048 }0.0065 |0.0570 }0.0550 |0.0585 0.0620 Pal05 5105 Palos 4108 0.0007 |0.0009 |0.006% Jo.oz60 |0.0295 |0.0330 Ag109 4109 0.0004 0.0004 ]0:0030 }0.0155 |0:0177 |0.0200 carts 4113, 0:0002 |0.00012 |0:00095 |0.0007 |o.0014 |0:0021 Ini5 4115 0.0002 |0.00011 |0;0004 |0:00035 |o:00042 |0:0005 112? 4127 0.0061 |0.0025° |o-:0022 }0:0037 |0:0028 0.0020 Xe131 4151 0:0340 |0:0290 |0:0340 }0:0380 |0:0340 |0:0300 C8133 4135 050583 0.0660 |0:0600 |0:0690 |0:0695 |0:0700 csi 34 4134 1135 6135 0.04594 |0.06015 [0.00671 |0.05878 |0.06207 [0.07173 xel35 4135 0:01388 }0:00558 |0:00033 |0:01252 |0:00562 |0:00283 cs135 5155 wards 4143 0.0589 |0.0575 |0.0485 |0.0400 |0.0455 | 0.0450 Nat4 4145 020332 |0:0380 |0:0400 }0:0310 |0:0310 [0.0310 pmi47* 4147 0:0094 Jo.o102 0:0137 }0:0095 |o.0102 }0:0109 Pmt 47 5147 0.0094 }0:0113 ]0-0153 Jo.0105 }O.0113 }0-0121 smi 47 e147 Pmi 48m wie pm 48 5148 Pmiag 51ug 0.00769 |0.0108 0.0167 ]o.0128 _]0.0151 |0.0146 smi 49 aiag sm150 4150 sm151 4151 0.0043 |c.0040 |o.0090 |0.0080 |0.0085 0.0090 sm15z 4152 0:0017 |0.0026 |0.:0057 |0:0000 |0:0065 {0.0070 Bu153 4153 0:0009 |0:0015 |0:0038 J0.0034 |o.0044 0.0055 Bud) 4154 Bu155 4155 0.00017 }0.00029 }0.0014 |0.0016 |0.0021 |0.0027 Gat57 4157 0:00006 |0:00007 |0.:00017 |0-0007 |0:0005 |0:0004 Pseudo 4907 0:3170 |0:3040 |0:3070 }0:3290 }|o:3190 0.3090 uesT 927 0:0600 * Two iaentical nuclides, Pmi47, are included to model the branching ratio to Pmi46m and Pmi48. ABEW-R 2133 36TABLE 10 GOLDSTEIN-COHEN A-VALUES AND AVERAGE POTENTIAL SCATTERING CROSS-SECTIONS FOR THE PRINCIPAL NUCLIDES WIMS op Identification A barns Number 3001 1,00 | 20.3 3002 1:00 3:35 10/1010/1011 0:99 430 2012 0:98 aT 2212 0:0 020 14 0:98 8:9" 6016 GD 3.7 27 0.88 V4 29 0:85 2.3 52 0:56 6:2" 55 0.0 0:0 2056 GD 10.24% 58 0.52 | 16.9% ca cD 6:9 4083-4902 (F.P.) | 0.0 020 2154/2156 /2158 0:0 0.2 2155/2157 0.28 7.0 207 o:22 | att 2232 0:20 | 11.2* 9233 0.20 | 11:9% 1233 0.20 | 12.2% 234 0:20 | 16:0 2235 0.20 | 12.8 236 0.20 | 14:7 8238 0.20 | 10:6 6239 0:20 | 12:78 1240 0:20 | 11.8* 1241 0.20 | 11:6* 1242 0:0 0.0 951/952/953 0.0 0-0 1000 0.0 0:0 GD values are group dependent. See Table 11. * Group values vary significantly from values quoted.9S TABLE 11 GROUP DEPENDENT A-VALUES FOR OXYGEN IRON AND ZIRCONIUM FOR 6 ;ROUP LIBRARY A-Values WIMS Group Oxygen Iron Zirconium 6016 2756 9 15 +9822 «9857 +9853 16 +9930 29914 19879 17 +9938 +9851 +9830 18 1.0040 +9898 +9803 19 +9908 +9862 +9671 20 1.0007 29875 29737 21 +9831 +9309 :8796 22 +9805 29038 +8158 23 +9978 -8391 +6833 24 1.0247 «5480 «3386 25 1.0448 2S4Tk «3309 26 1.4852 (0, =1,max_ group number 13:24 ‘Downscattering cross section. ‘The above described output on units nt and 19 appears for ICAR=0 and ICAR=1. For ICAR=1 additional output is written to unit 14. However, this is meaningful only under two conditions: - number of collapsed groups equals 2, - the system calculated includes: bboron-10 in water, ‘Xenon-135 and Samarium-149 in fuel. The library numbers of those isotopes can be specified by the user on the ELEM input card if the default values (see Section 4) are not appropriate. The cross sections written on unit 14 are described in Table 3.3, The ANSWERS Software Service WwimMs AEA TechnologyWIMSD NEW OPTIONS Table 3.3, Cross sections written on unit 14. Card number | Columns Description 7 72 ‘Bumup in Maite 13-72 | Alphanumeric tite defining the materials of the current cross section set z 12 3 absorption 13-24 6 fission 25-36 transport 37-48 v 3 1-12 1120.0 always 13-24 12 ‘The unit 14 cross section set has been developed for calculations of PWR operation. The cross. section set from Table 3.1 is therefore repeated in each burnup step and each region for the following materials: 1. Macroscopic cross sections of the fuel assembly (the whole system treated in the WIMS run) 2. Microscopic cross sections of Boron-10 3. Microscopic cross sections of Xenon-135 4, Microscopic cross sections of Samarium-149 5. Macroscopic cross sections of light water In set number 1 the contributions of the Boron-10, Xenon-135, and Samarium-149 are subtracted from the macroscopic cross sections, ‘The procedure followed by WIMS to generate the cross section library consists of the following steps: ())_ In CHAIN 2 the microscopic cross sections for the nuclides H, B-10, O, Xe-135, Sm-149 are stored in 69 groups on a scratch unit 10. (il) In CHAIN 5 these microscopic cross sections are condensed to the number of groups of the main transport routine, using the fuel spectra for Xe-135 and Sm-149, and the coolant spectrum for H, B- 10, O. The condensed cross sections are stored on scratch unit 13. All the cross sections are condensed by a direct weighting, but the transport cross sections is obtained following the equation: 1 1(6,,-6,) where o, and ,, —a, are condensed by a direct weighting. (ii) In CHAIN 14, these condensed cross sections are condensed again to the two groups of energy specified by the THERMAL card and for each region. A direct weighting is used for all cross sections except transport where the Bondarenko equations are applied, to preserve the reaction rates and leakages. It should be noted that for the option ICAR additional tape units are used - unit 10 -scratch, unformatted, - unit 13-scratch, unformatted, The ANSWERS Software Service wims AEA TechnologyNEW OPTIONS WIMSD = unit 14 - formatted, for the library if ICAR=1, - unit 19 - formatted, output in CITATION format, - unit nt - formatted, output in LAMS format. ‘The user should decide how to define the above units in his ‘datsets’. It should be stressed that the option does not work together with the CALEB option as both these options use the same variables IGP and IGPCUT specifying the group condensation. 3.2. Modifications taken from the WIMSD-4.1 version (33). A set of additional options has been introduced: (An approximate approach has been introduced into the burnup calculation which assumes that I-136 decays instantly to Xe-136. In this case the Xe-135 number density is independent of the number density of I-135. A further assumption is that Xe-136 is formed only by fission and I-135 decay, while 1-135 is formed only by fission. With these assumptions. the equilibrium number density of Xe-135 is, given by: Yq.¢FR Ny, = Age * hye 61) where: Yr ~ cumulative Xe-136 fission yield, including that of 1-135, FR - total fission rate, Aye total Xe-135 absorption rate divided by Xe-135 number density, 2. ye~ decay constant of Xe-136. Equation 3.1 is used under the option EQUILIBRIUM for the first intermediate burnup calculation between lattice calculations. (ii) Cluster geometries have hitherto been treated as isolated in Dancoff calculations. However, in the case of LWRs it is more appropriate to treat the cluster as a part of an infinite system. This can be achieved through the LATTICE option. (ii) An alternative resonance integral interpolation method has been developed by Segev [34] and included as an option SEGEV. (iv) During a burnup calculation all isotopes with a decay constant greater than 10 have been treated as if they were in equilibrium. A new keyword NONE (NONEQUILIBRIUM) has been defined that causes this test to be bypassed. It should be noted that the keyword SATURATE, which forces a specified isotope to be in equilibrium, remains in effect and overrides the new keyword for the specific isotope. It should also be noted that for long burnup calculations, attempts to follow xenon explicitly can lead to oscillations caused by loss of numerical accuracy. The option is therefore recommended only for short xenon transients, 3.3. Specific features of the WIMSD-2B version (i) The BINPUT subroutine (together with the associated ALPINT, RECOG, CALPIN) have been replaced by a READIN subroutine. This is part of the ‘READDATA’ package which consists of: - READIN subroutine called by the WIMSD-5A program, - auxiliary subroutines called by READIN: ERROR LNFLOW READLN The ANSWERS Software Service wims AEA TechnologyWIMSD NEW OPTIONS READNO SKIPSP TYPE READJN -alternative to READIN, not used in WIMSD currently. Itis based on reading 8 instead of 4 characters of the keywords. ‘Subroutine OPENCH is an addition to WIMSD and is used to assign unit numbers and names of the files such as: library, input, output. They are specified in the fixed name file ‘datsets', where they can be modified according to the user's needs. (ii) The addresses of arrays are passed at the level of the subroutine using those arrays and not all of them are passed at the level of ICHNm subroutines. This causes a different structure of COMMON blocks and arguments as compared to WIMSD-4 (ii) In some subroutines double precision has been introduced (iv) Numerous miscellaneous modifications and improvements have been introduced, eg, more ‘meaningful codewords to replace SEQUENCE and CELL, The user will observe a reduced output and a set of new error messages. The use of the new READDATA input package affects the input user image, as follows: (i) No spaces are allowed in the keywords, eg., the card for transport group specification has to be written: FEWGROUPS or FEWG The form FEW GROUPS will cause an error. (i) All the cards specifying PRIZE input, given previously without a keyword, must now have the keyword PRIZE at the beginning of each line. i) Continuation symbols $ are not required. (iv) Repeated data items are permitted, thus: SUPPRESS 13 @ 101 would give output from chain 14 only, or CODEWORD 4 (1 2) would give CODEWORD 12121212 Note that now the characters ( ) must not be used as separators, ‘The ARIADNE DIFFUSION COEFFICIENTS calculation follows the WIMSD-2B ie, the original WIMS algorithm. This may cause differences in those coefficients calculated with previously distributed NEA versions of WIMSD-6. 3.4 New features ‘The possibilty of microscopic library creation on unit 14 has been extended by allowing input IDs for 5 elements (keyword ELEM) This permits different libraries to be used The ANSWERS Software Service wims AEA TechnologyNEW OPTIONS WIMSD Provision has been made for the icar=1 option to give macroscopic fission cross sections in user chosen regions without creating the unit 14 library. Informative and error detection messages connected to this option have been added. In particular, the option of additional group and region condensation specified by keywords MOMOD and VECTOR in the the edit section of input data is not allowed with the CALEB option active (ICALB=1) as in this case the program variables IGP and IGPCUT (needed for additional condensation) are used by the CALEB option. ‘An efor in cluster plate calculations has been removed. The ANSWERS Software Service wiMs AEA TechnologyWIMSD INPUT 4.1. General Structure of WIMSD-5A Input 0} (i) ‘iy (w) “ The WIMS input is in a general free format. Each item consists of a keyword and, usually, of consecutive digital entries divided by one or more spaces. The input routines of the code read 4 alphabetic characters and, if these are recognized as a keyword, the digital entries following the keyword are read and interpreted. The exception is a card with a star's’on it: the text after the star is interpreted as a comment extending up to the end of that card. The keyword must not contain spaces. There are occasional exceptions to the above rules where a keyword has to be followed by a formatted set of data. The exceptions relate to more advanced options of the code and they are noted specifically in the following description: ‘Anew keyword usually starts on a new line but the data prescribed to it may be continued on the next card(s). (Keywords may continue on the same line but this is usually harder for the user to read.) Repeated data items are permitted in two forms: CODEWORD n @kIm ives n times the value of k followed by |, metc., €.9 SUPPRESS 13@ 101 gives the output of chain 14 only, CODEWORD 4 (1 2) is equivalent to CODEWORD 12121212 ‘The WIMSD-SA input is divided into 3 sections read consecutively in 3 different parts of the code. ‘The sections are divided by keywords. The first section, known as the Prelude Data, serves for fixing the problem dimensions. It is read in subroutine PRELUDE and it ends with a keyword PREOUT or ENDP ‘The second section contains the main data describing the case. It is read in subroutine DATAG of CHAIN 4 and, in the first lattice calculation, starts with a keyword: INITIATE and ends with a keyword: BEGINC ‘The third section contains the Edit Data. It is read in subroutine DATAR of CHAIN 13 and it too lends with the keyword BEGINC. Multiple cases may be executed in the same run only if the Prelude Data remains unchanged. All that is required is a set of those Main Data items which are to be changed, followed by BEGINC, and a set of Edit Data items which are to be changed, also followed by BEGINC, Ifo items in any section are changed then BEGINC is sufficient. Thus a sequence BEGINC BEGINC placed after the last BEGINC of the previous calculation causes a repetition of the previous case. ‘This is widely used for burnup calculations where a double BEGINC is placed to execute the next burnup step with the same calculational parameters, but starting with the depleted isotope number densities, The ANSWERS Software Service WiMS AEA TechnologyINPUT WIMSD 42 List of Input Data In this section, a full input description for WIMSD-5A is given, based on the input description of WIMSD-4 distributed previously by NEA Data Bank. Section 1 Prelude Data HOMOGENOUS, PINCELL CLUSTER ‘An alternative to the above three keywords is: CELL i NPLATE SEQUENCE i DSN PERSEUS Pi PRIZE NGROUP ijig nticar Homogeneous medium calculation Pincell with or without condensation or single plate calculation Cluster calculation (of pins or plates) Homogeneous Pincell, no energy condensation Pincell, with energy condensation Cluster iti positive: a slab geometry to be divided into i unit cells for resonance shielding and multigroup calculations, if i negative: a spherical geometry; in that case, only DSN is avaliable and all geometry is input as for cylindrical geometry except that the annulus radii must be r= /4eR? /3 so that x «72s the required spherical volume. choice of the main transport routine: i=1-DSN, i=2-PERSEUS, fori = 1,2- in "CELL 7" a negative sign causes the unsmearing process to be omitted. (Do not use in Burnup calculations). i=4 -PlJ - PERSEUS, 5 - PRIZE - PERSEUS. To change sequence option during a WIMS tun, see SWITCH data, The SEQUENCE card can be replaced by one of the 4 cards specified below except when a negative iis required DSN chosen. PERSEUS chosen, PJ chosen. PRIZE chosen i - the number of groups in the main transport solution, j- the number of reaction edit groups, ig - number of groups for additional condensation, equal to the number of entries on the VECTOR card, igs10, AEA Technology The ANSWERS Software Service WiMsWIMSD INPUT NMESH ij NREGION ikigr NMATERIAL ijk NRODS ijkimnan nit - logical unit for condensed output in LAMS format, int=0 - the condensed cross sections are not determined, int=8 - besides the LAMS format, the CITATION format cross sections are written on file 19, intz8 - only the LAMS format cross sections are written on fle nt. icar - if iear = 1 the macroscopic fission cross sections will be calculated in the second Condensation; in the presence of isotopes with ID = 3001, 1010, 6016, 4135, 4149 (or other for isotopes specified on the ELEM card) and second condensation to 2 groups. a library on unit 14 will also be created Remark: if icar is not equal 1, it should be zero or blank. Also i must be < the number of library groups for PRIZE and icar = 1 options. Items ig, nt and icar may be omitted i- the number of mesh points applied in the ‘main transport routine: for PRIZE calculations i ‘must be equal to the number of PRIZE regions; for multicell cases, iis the sum over all cells, j- the number of edit mesh points (default = i + 3 times number of annuli containing rods). If the SWITCH option is used to change to PlJ, the number of edit meshes will be increased automatically to the number of Pld meshes. - the number of slabs or annul, j-the number of annuli containing rods (default ), (Serious errors in PLPERSEUS if this is not correct), k-the number of edit regions (default = 34), For multicell cases, k is the total number of edit regions and the default is likely to be too small If the PRIZE option is chosen k must be present on the card and calculated as described in the PRIZE option description. igr - number of regions for the second condensation (default zero) For multicell cases, iis the maximum in any one cell type and j is the sum over all cell types. k is also the total for all cells, so the default is likely to be too small = the total number of materials, j-the number of materials undergoing burnup, k- the number of materials used in PERSEUS. If k= 0 or blank the code supplies an adequate default value, - an explicit value for k could be useful for saving space in multicell jobs). Data required only for PIU or for DSN/PERSEUS calculations with more than 5 rod-subs in any type: The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD NTUBE NREACT NCELLS POISON NTHERMAL, NISOTOPES TUT ntube i - number of rods in the cluster (i21), j-cluster symmetry factor > 0 for j fold rotational symmetry, J <0 for +2 fold rotational and j fold reflectional symmetry, with =0 being a reflection plane, k - number of lines in Pid mesh (50 or 100 mostly used depending on the case, consult printout of CHAIN 8 with exact and numerical volumes in particular cases), I number of angles in PlJ mesh (best chosen to be co-prime with problem symmetry) m - number of rod types, 1n- maximum number of rod-subs in any type, nn - maximum number of sectors in any rod- sub (The last 3 default to 5,5,2 if none is input). ‘An approximate resonance treatment for tubular fuel elements uses the plate fuel algorithms. tube = number of concentric tubes in a single fuel pin (see main RTUBE data). i number of nuclide reaction rate edits. i - number of cell types in mutticell calculations. = number of regions in a subsidiary burnable poison pin cell calculation. (Default 7 or 2 +n of the NRODS data). i - the maximum number of temperatures at which thermal library data is tabulated for nuclides in the WIMSE-type file output (see SAVE option). The default value of 7 is, adequate unless the following identifiers have been used: 1012(8), 2012(14), 4016, 5016, 6016, 4235, 6235, 4238, 6238, 6238 (all 9). 1 the number of nuclides used (default = ‘number of nuciides in the library). This allocates sufficient data storage space to calculate heat flux tit across cluster fuel pins. No main data required. AEA Technology The ANSWERS Software Service WiMSWIMSD PRIME iddimmlyy SAVE ijkima INPUT Card required with i = 0 if the job is writing WIMSE-type output to a_new tape. A non-zero value of i may be specified if the user wishes to save only the output from his first i WIMS runs. An expiry date in the form ddimmiyy is optional but should be used as a safeguard against overwriting if the card is not_ removed for subsequent jobs. i- the number of previous sets of WIMS output to be skipped before the present run writes its output. For the first run with a new tape, i=0 is required. For subsequent runs the code automatically skips to the end of the previous ‘output if this card is not supplied. If used this must be the first PRELUDE (Section 1) data card Option to write WIMSE-type files [26] i= -1 - no action taken except execution of the PROLOG subroutine i = 0 - WIMSE-type files 1-5 will be computed and printed if required i= 1 - Files 0-6 wil be output to tape with or without the CALEB edit card = 0/1 - the full WIMS input data for the job will novwill be printed before execution. ‘commences, k-- the length of a File 6 record (default = 100). 1 any value (this parameter is no longer used) m - if non-zero, extensive printout will be ‘obtained from the datablock input routine. n= 1 - a call to subroutine PROLOG is suppressed n= 0 (or blank) - the call to PROLOG is executed (PROLOG js a subroutine which automatically supplies most of the PRELUDE data including the CELL option and all options from Keyword NGROUP to NTHERM inclusive. ‘The user is still required to supply those cards which define a route of calculation - SEQUENCE, TILT, STORE, LIBRARY, PREOUT - and also. NPLATE and NISOTOPES if required. The Keyword LINES in the Main Data (Section 2) should be used to complete the NRODS specification. All user- supplied PRELUDE data will overwrite PROLOG data so that existing data decks require no modification) The SAVE card must be followed by a card containing a 20 character run identifier (cf. ref 25). Execution will be terminated if File 6 output is requested without a coolant card in Section 2 and Xe-135, U-235, Pu-239 reaction edits requested in Section 3 The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD STORE ijk LIBRARY PREOUT or ENDP Section 2 - Main Data INITIATE Geometry data ANNULUS irm SLAB irm ‘SQUARE irm If k > 64 and a reaction edit is requested for ID 1000, fast and thermal cell averaged tv reaction rates in (secicm) will be stored in positions 63 and 64 of File 6 When all data items are greater than 100, the ‘pseudo-tape’ option of the code (ex-WIMS3) will be used and all execution will be in-core Typical values for a 28 group library are 50000, 10000, 10000, 1000. The required value of i for the Winfrith 1986" 69 group library is 184000. i- the storage allocated for library data, j- the storage for dataset 3, k-- the storage for dataset 6, Datasets 3 and 6 are scratch discfles if this option is not used (ex WIMSD2). I the storage allocated for data transfers. The code prints the required library storage and indicates the required values of j and k by comments of the form: “3>mmmm" and *é>mmmm" in the output Used only with ‘pseudo-tape’ option (invoked by the STORE card)) ‘A map of library data storage is produced - indispensable if library modifications are to be specified at run time. The last card ending Section 1 ‘The first card of a complete set of Main Data. j-the annulus number, counting from the centre, r the outer radius in cm, m - the material number for all parts of the annulus not occupied by rods. Note that if intermediate annuli as are usually required for PlJ calculations are not defined in the input, the code automatically supplies appropriate data by linear interpoiation in r. For ‘multicell cases, annul starting at 1 are specified for each cell type. j- the slab number, counting from the centre, tr- the distance from the centre to the siab outer edge in cm, m- the material number. Note that symmetry is assumed about the centre line. j- the outermost annulus of the problem, AEA Technology The ANSWERS Software Service wiMsWIMSD INPUT CRAIG POLYGON inmr ARRAY Nimnre) RODSUB Nkr(me) r - radius of the inscribed circle of a square boundary which is orientated so that its corners lie on the x and y coordinate axes, ‘m - material number for all parts of the outer region not occupied by rods, If annuli are specified with radii greater than + and less than r2, segments will be correctly masked off by the square boundary. Rods of spectrum type 4 materials may be located at comers of the square and will not count as, rods in the averaging procedures at the level of multigroup calculations, The problem will be suitably homogenized and cylindricalized for DSN or PERSEUS but will be calculated explicitly by PIJ ifthe number of PJ annul is i For cluster geometries with one annulus per ring of fuel pins, the radii of the annuli containing pins will be adjusted to preserve the fuel to coolant ratio. throughout the cluster, and to minimise the differences between the fuel positions and the midpoints of the annul j-the "annulus" number, m - the number of sides of the polygon, 1n- the material number, += the “radius” to the midpoint of a side. This option merely generates appropriate ANNULUS data and has no influence on the calculation flow. Al rods of type N have their positions given by one ARRAY card. Up to 12 sets of the four numbers (m nr 0) are allowed, with m being an option on the coordinates used: m = 1: a ring with radius r of n rods, one with angular coordinate @ in radians, m= 2; (r,8) specify (xy) coordinates of one rod, m = 3: a ring of n rods with separation +, one with angular coordinate 0, m = 4: a hollow square of n rods, with separation r, 9 is the angular coordinate of the first corner of the square, m = 6: a hollow hexagon of n rods, with separation r, 0 is the angular coordinate of the first corner of the hexagon. If the cell boundary is square, non-fuel rods may be placed at the corners (cf, SQUARE option). They will be smeared for DSN. or PERSEUS. N-- the rod type (see ARRAY card), k - the rod subdivision number counting from the rod centre, The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD PLOT MATERIAL TEMPERATURE rs ‘Composition data matniist m-Atniist -mdtnlist mM list r- the outer radius, (m 0) - pairs of numbers giving material number and angle at the end of each sector. For DSN and PERSEUS only one m needs to be specified After the Main Data has been read, the geometry is "plotted" on the line-printer. This option may be used to redefine the maximum outer radius of the plot, r, and the scale, 6, to which the pins are plotted. Default values are: r = outer annulus radius, s = 0.7 (for improved readability). m - the material number, = the density in g/cm? t-the temperature in °K, fn - the spectrum type: 1 - fuel, 2 - can, 3 - coolant, 4 - moderator, ‘A negative n causes the material to be excluded from the resonance treatment and the SPECTROX spectrum calculation. This means that the representative cell will not include materials with negative n, although these materials will get the multigroup spectrum calculated for those materials with |n|, list - pairs of nuclide identifier and weight per cent (the code renormalises these to total 100%). mtn- defined as above, list consists of pairs of nuctide identifier and number density in atomsi(b.cm). If the material has been specified in a previous lattice calculation, any nuclide not specified jin list will retain its existing number density The material number may be specified negative, d is redundant, t and in are defined as above, list consists of pairs of nuclide identifier and a constant by which the number density of a ‘nuclide in an existing material will be multiplied If only two numbers are specified, material m will be given the same temperature, spectrum type and composition as material M (which must be specified previously). Note that subsequent redefinition of M does not affect m (used in burnup calculations with several burnable zones) - string of temperature increments, one for each material as defined by MATERIAL data AEA Technology The ANSWERS Software Service WiIMSWIMSD INPUT DENSITY ENRICHMENT RESTART COOLANT DTDR ELEMENTS. BELL DANCOFF The ANSWERS Software Service wims list mdes nd list (ghn), list - string of constants by which existing material densities are multiplied. (The code will fil out a truncated list with 1s), m- material number, d - density [g/cm ], (temperature and spectrum type must be defined by a MATERIAL card), e- enrichment (wt % of total heavy atoms, of U- 236 if uranium fuel, or of total plutonium if plutonium enriched fuel. In the latter case the input ratio of U-235/U-238 will be unaltered), s-the stoichiometric ratio (eg., 2.0 for UO>). ‘An enrichment specification will over-ride the isotopic specification of a MATERIAL card but ‘an appropriate default_value of s will be derived from the MATERIAL data if zero is, input. Several sets of (mde s) values can be input on one ENRICHMENT card. Use of this keyword requests that all material number densities be read from Fortran dataset 3, which has been saved from a previous job. ‘All MATERIAL data will be ignored. n - the material number used for evaluating coolant density and temperature for File 6. d_- the historic coolant density for File 6 (defaulted to the intial density of material n if no value is supplied). This card is mandatory if File 6 output has been requested (SAVE data of section 1) ‘The temperature of all materials with spectrum type 1 (fuel) will be overwritten by T + derq, where T is the temperature of the "coolant" ‘material and rq is the rating taken from the second item on the POWERC card. The subsequent changes in rating will have no effect on fuel temperature unless another DTDR card is supplied list - up to 5 identification numbers from the library defining nuclides used for unit 14 generation under iear = 1 option, denoting H, B 10, O, Xe-135, and Sm-149. Defauit numbers are: 3001, 1010, 6016, 4135, 4149. The list may be truncated at any point. (all items optional), ‘a - the Bell factor (default = 1.16) list of triads (g hn), with the following meaning: 9 - the infinite lattice Dancoff factor, AEA TechnologyINPUT WIMSD REGULAR an PIJDANCOFF 9 LAMBDA ist RESXSECS rm list NORESONANCE h - the average Dancoff factor including the outermost ring of rods (h>g) n- the library resonance group number @ and h apply to all library groups (n.4#1) tom with ng=14 in a 69 group library. If n is absent single values of g and h apply to all library resonance groups. a- usually 1. By specifying a + 1 a correction is introduced to take into account the irregularities of the lattice, 1 - can be biank or 6, the latter denotes a hexagonal geometry. This option is recommended for accurate calculations of regular lattices, especially for tight pitches, but should not be used with DANCOFF. ‘Automatic computation of cluster Dancoff and Bell factors. The job must be set up as a Pl run and g defines the main transport energy ‘group for which factors are evaluated. The recommended numbers of lines and angles in the PlJ mesh are 100 and 11, respectively. No edit data (Section 3) is required for this option and the computed factors are printed and saved for the next lattice calculations. 1 - can 4 value for the PIJDANCOFF option (default 10.3923 - appropriate for Zircaloy cans), list - can contain up to 6 numbers defining the lJ mesh regions containing cans. Defaulted values are correct unless a gas gap between fuel and can is specified by RODSUB data. To input resonance cross sections: nuclide resonance identification, eg: 2235, m - type of cross sections: 1 is absorption, 2 is fission, list - inner cross sections for each resonance group for each type in tum, followed by a similar set of outer cross sections. ‘The input cross sections replace those from the output of CHAIN 3, appearing in the first and third column under the tile "MICROSCOPIC. CROSS SECTIONS" The infinite dilution cross sections will be used the resonance shielding calculations being bypassed. AEA Technology The ANSWERS Software Service WwiMSWIMSD INPUT LATTICE SEGEV RTUBE FEWGROUPS, MESH DBSQUARED n n= -1 - cancels a previous request for that option ‘n= 0 (or blank) - the Dancoft factor calculated for inner fuel pins is used for both inner and outer pins of the cluster. a n= -1 - cancels a previous request for that option n= 0 (or blank) - the Segev interpolation is Used in resonance integral 6, interpolation 11,12,03,f For an approximate representation of resonance shielding in tubular fuel elements, t1 = fuel thickness clad thickness (inner and outer) 3 = coolant thickness (between consecutive tubes) f= volume fraction of fuel which takes the ‘outer’ resonance integrals, Gross sections and flux solution list list - the division of library groups into main transport routine groups. The length of the list is the first number on the NGROUP card of Section 1. The last number is the number of library groups. This may be omitted if no condensation is required list list - the number of mesh intervals of equal volume in each of the annul The length of the list is the first number of the NREGION card of Section 1 The sum of entries on this card should be equal to the first entry on the NMESH card. In PiJ and PRIZE calculations, the individual values are irrelevant; only the total must be correct. Separate MESH data is required for each cell type in multicell calculations. Leakage term DB? added to 5, in few groups: i= -1-- list of g values of B?, where g is the umber of main transport groups. The term B2(32,,) is added to £4 for each material. = 1 - list of g values of B? and g values of D. DB? is added to ©, for each material 2-- list of g values of B2, g values of D,and a trigger j for every material: J =-1, B2/(33,) is added to Z, , 0-5, is unaltered, J=1, DB? isadded to x, Note that material 0 (void) is never altered. The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD DIFFERENTIAL, CARDS HPS FREE TOLERANCE OMEGA Note also that if |i] is increased by | 10], ie., to 11, 11 0F 12, the DB? terms will be removed before the edit data is executed If i is omitted thermal hyperfine disadvantage factors will be calculated for each fuel ring in a cluster, i= -1-a request for the option will be cancelled, Cross sections in DSN format are output to unit 1 by CHAIN 5. This unit may be defined in DATSETS to ensure that itis not a ‘scratch’ file. This option has 2 independent effects: 1 - omitted - burnup calculations to be carried ‘out in a k-infinity spectrum, regardless of input buckiings 11> 0-ifn > 1000, m=n-1000, macroscopic cross sections in DSN format, ie as produced by the CARDS option, are read for the first _m materials and over-ride the respective material data from the input Data to be punched after the Main Data BEGINC card includes the cross sections in DSN format for materials 1 to n-1000, and the material code for each mesh. The material code should follow the DSN 5's when a condensed group structure is used (because they are read in Chain 5), and precede them otherwise (because they are fead in Chain 6 when there is no group condensation). If PlJ andlor PERSEUS are used 5's are read in the same way but only if a condensed group structure is used. No material code is required. A free (black) boundary condition (-> 0 at some extrapolation length). Note that group dependent albedoes can be input to a collision probability calculation by use of the PCELL ‘option, even in single cell cases. ‘a - convergence tolerance of the main transport, routine, default value = 0.0001 a - acceleration parameter for the collision probability solution routine, default value = 1.25, which is appropriate for water systems; for graphite systems higher values (eg 1.6) are better. AEA Technology The ANSWERS Software Service WIMsSWIMSD INPUT s A NPIJAN iti LINES ijk Burnup POWERC. {irq rtau indb rmaxt indg list The ANSWERS Software Service wims n - order of DSN calculation (SEQUENCE 1), default value = 4, but any even value > 4 may be input for higher accuracy j - number of annuli (including a square outer region) to be treated by Pld, all outer annul to be treated by PERSEUS. The first j numbers on the MESH card are arbitrary except that they must add up to the total number of zones in the PlJ calculation (consult printout of CHAIN 13° in more sophisticated cases) t- maximum number of mean free paths for which neutrons will be tracked, default value = 10, recommended value = 7 to save computing time i - if nonzero, track renormaiisation will be suppressed. The renormalisation is used to correct for inaccurate volume integration. i - number of lines in PIJ integration mesh, |= number of angles in PJ integration mesh, (These values of i and j will override values on the Prelude NRODS card) k = maximum number of boundary reflections which are tracked in a square PIV calculation, (default = 2). i the units of rating for rq, 4 - rq is MWite initial heavy elements, 2- rq is fissions/om? /sec cell averaged, i = 3 - rq is fissionsicm? /sec averaged over burnup materials, 4-1q is total fux, i= 5 - rq is flux in group indg, If rq = -1 the value reached by the previous burnup step is assumed. rtau - the timestep between criticality calculations and flux renormalisation to the input power level (days) indb - the number of such timesteps between full lattice calculations, rmaxt - if positive is the maximum fractional ‘nuclide density change per integration step (default= 0.05), if negative is a constant integration step length in days. indg - the group used to normalise the flux level in the i= 5 option (default = 1) In options i=4 and i=5 list contains a trigger for each material which is 1 if the material is to be included in the flux normalisation and 0 otherwise. If no values are specified in the list, all are assumed to be 1 AEA TechnologyINPUT WIMSD opps ij BUCKLING list WRITE FUEL ablist POISON ist 1+ 0 - 2-group microscopic effective absorption and fission cross sections are printed for all unable nuclides. Group 1 contains the first li] main transport groups, group 2 the remainder, If |i) equals the number of main transport groups, 1-group data will be obtained j +0 - multigroup fluxes will be read for the fuel, coolant, moderator and can from the input and will substitute the spectra calculated for a representative cell, The reading is carried out in CHAIN 5 and, thus, the fluxes must be given after the frst BEGINC in the format 612.5 (see ‘SPU option). lit consists of B2, , B2, , 82, , B2, where 82,,, 82, enable a if Spectrum to be ised for buthup. If 82, , BZ, are speciied critical spectrum will be used with bucklings B2, + ie B2, and B2,+16B2, with i chosen to give ett = 1 If the list contains more than 4 numbers it is assumed that values have been specified in sets of 4 for each group. Without this card Ariadne diffusion coefficients are used in the burnup calculation if available (SEQUENCE 2), otherwise transport averaged D's are used. This option results in the use of D, and D, as computed by the leakage edit, the order of preference being Ariadne, Benoist and Sy, ‘As bumnup is normally specified in MWdite of initial heavy atoms it may be desirable, on restarting a_burnup run, to specify the original fuel mass. The appropriate value of a is printed in the burnup output of the code under the heading "INITIAL INVENTORY". If a=-1 burup will be evaluated as MWdlte of current heavy atoms. The fuel mass is recalculated at every burnup step, b - an irradiation increment for File 6 output and is required only if a file output job is restarted from a non-zero irradiation list - remaining percentage weights of each burnup material and the cluster average at the restart irradiation, Appropriate values can be obtained from the WIMS bumup_ printout Default values are 100.0 The burnup of highly absorbing poison pins in a cluster can be refined by means of a pin cell colision probability calculation at each criticality calculation in the burnup module. AEA Technology The ANSWERS Software Service wimsWIMSD INPUT NONEQUILIBRIUM EQUILIBRIUM a i - number of regions in the pin cell calculation (including poison, can, coolant, and paste regions) (i < Section 1 POISON data value), jj - number of regions in the poison pin itself list defines the geometry of the poison cell, and contains i-1 sets of region number region radius = material number = volume fraction (= 1.0) followed by data for the paste region containing region number region radius = up to 3 pairs of material number and volume fraction (totalling 1.0) n option, 11 = 0 (or blank) - all isotopes are treated as not in equllibrium. The option requires a sufficiently short time step (s 1 day), but is not recommended except for short xenon transient calculations. = cancels a previous request for that fn + cancels a previous request for that option. n= 0 (or blank) - a modified equation is solved for isotopes in equilibrium (Xe-135) [32] (See also ALPHA and SATURATE in the Edit Data Section) Fuel cycle calculations CYCLE i ky, X 8X rtau indnb w j To perform automatic fuel cycle calculations i=1 oF 2 (see definition of X) kn - the irradiation averaged value of Key at which burnup will terminate, X #8X - the required initial key if i = target irradiation if =2. rtau indnb - replace the standard quantities on the POWERC card after the first burnup step, and determine the length of all subsequent steps. w = the irradiation at which plutonium will be extracted for recycle purposes (end of cycle if zero or not specified) 1= 0, ky, defined by averaging Key = 1, kp defined as the ratio of integrated neutron "production to integrated neutron losses. ‘An entire fuel cycle is treated as a single WIMS lattice calculation for data preparation, i. e., only one set of main and edit data is supplied or the The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD RECYCLE PROCESS ENRICHMENT DIFFERENTIAL CELL AEA Technology ir tit list list list It must be noted that several assumptions are made in the CYCLE options concerning the nuclide identifiers. These must be checked if the Winfrith 1986 library is not used. To continue a CYCLE run by recycling fuel J =1 - to the basic fuel (as specified by MATERIAL cards) will be added the plutonium from the previous cycle + specified enrichment to meet the objective specified by the CYCLE card. j=2 - to the basic fuel will be added sufficient plutonium of the composition extracted from the previous cycle to meet the objective, J =3 - to the basic fuel will be added only that plutonium extracted from the previous cycle, the objective is ignored and the cycle terminated when the average kai iS Km r - an estimate of the required enrichment relative to that required for the previous cycle (defauit = 1.0) ‘The composition of the recycled plutonium should be corrected for Pu-241 decay and processing losses. ty the time in days between removal of the fuel from the reactor and processing. ty - the time in days between processing and refuelling. f - the fractional efficiency of plutonium extraction (The assumed Pu-241 half life is 14.4 years). Note that in all cases other than RECYCLE options this card is interpreted as under the MATERIAL DATA Section In a RECYCLE calculation, list specifies the relative amounts of U-235, Pu-239, Pu-241, Pu- 242, Am-241 in the extemal enrichment for RECYCLE option 1 list contains up to 6 numbers specifying the required enrichment ratio in the fuel materials. For example, if materials 7, 8, 9 were fuel the first three items in the list would determine the proportional enrichment in those materials (default is uniform enrichment) In a multicell calculation all ANNULUS, MESH and ARRAY data refer to the immediately preceding CELL card, but MATERIAL and RODSUB data apply globally i-celltype, The ANSWERS Software Service wimMsWIMSD INPUT n- number of cells of this type, default =1 To minimise round-off error it may be advisable to number cell types in order of increasing n values. SPECTRUM m (or blank) If this card is supplied with any value m # i on the CELL card, the condensation spectra and microscopic resonance cross sections for cell | will be set equal to those of cell type i PCELL Pj for allj i-celltype, Pij - probabilities that a neutron leaving cell type i will immediately enter a cell of type j, without necessarily colliding in that cell. The list Of Pij values can be repeated for each group in tum. When the list is terminated the last set will be used for all remaining groups. ‘SWITCH ijk This option may be used for changing the ‘SEQUENCE option in a burnup run: i - new SEQUENCE data (eg. 1= DSN, PERSEUS) J - new number of annuli (replaces first itern on. the NREGION card), k- new number of meshes (replaces first item on the NMESH card), Jand k may be zero or blank ifthe new values equal to the old ones. If j changes the number of annul, the SWITCH card must be followed by a complete new set of ANNULUS and MESH cards, REPLACE i, ig list The following library modification options are available only when the in-core option is used (see STORE data of Section 1), The required information on storage locations is obtained with the LIBRARY option. ‘The library tape itself snot affected. i, and iz are the addresses of the first and last library locations to be modified list is ig -iy +1 new values for the library data. The list will be expanded with values equal to the last one specified if it is too short or truncated if too long. MULTIPLY, iy iglist The list will be expanded or truncated to i, - iy + 1 values as above and will then be used to multiply the library data stored in addresses i,to ig inclusive The ANSWERS Software Service wims AEA TechnologyINPUT WIMSD INCREASE iy i list Printout control TESTPRINTS SUPPRESS list End of data BEGINC n PRIZE calculations PRIZE or nnr nge REGIONS. PRIZE or zazbratbm GEOMETRY The list will be expanded or truncated as above and will then be used to increment (by negative amounts if required) the existing data in addresses iy t Note that only REAL data can be (easily) changed. More extensive printout is produced, particularly in the resonance and transport calculations (see Section 6.3 for details). list contains up to 16 integers, iy i, = 0 - normal printout in CHAIN n, (see Section 6) i, = 1 - suppresses all printout in CHAIN n except for input data listings, some (not all!) error messages, but including FORTRAN error messages. Default values are 0 for n = 1-15 and 1 forn = 16. It may be useful in a failed job to know that the output of the last suppressed CHAIN to be executed can be found in dataset 4 The final card of section 2. If n is omitted a normal lattice calculation proceeds. The only recommended alternative is n = 13 which causes execution to jump straight to Chain 13 and is used to generate more than one edit of a single lattice calculation (often used with the REGION option of Section 3). If SEQUENCE § (PRIZE rz collision probabilities) is chosen the data for PRIZE calculations must be read immediately after the first BEGINC. Each line must start with the PRIZE or alternative keyword and the data ‘must be given in the order specified below: innr - number of regions in PRIZE calculations, Inge - definition of geometry’ ge <0- symmetry in the plane ge = 0 - infinite cylinder, Inge = + 1 - finite cylinder, ge = + 2-slab 0; A set of nnr cards specifying dimensions of the regions and materials filing them. For each region (after the keyword PRIZE) za - bottom height (cm) AEA Technology The ANSWERS Software Service wIMsWIMSD INPUT PRIZE or TRACKS BUCKLING NOBUCKLING THERMAL, BEEONE DNB fanrr nth nnzz rt r2 Section 3 - Edit Data Leakage calculations B?, B2, n D1 D2D3 D4 The ANSWERS Software Service WIMS zb - top height (om) (for infinite cylinder za - 0, zb = 1.0). ra - inner radius (cm) rb - outer radius (cm) (for slabs, ra and rb are omitted) m- material 1 and r2 are inner and outer radii in om for the mesh of the “line set’ or tracks used for integration of colision probabilities; if 2 = 0 the maximum radius of the system is taken. nner - number of equally spaced tracks between r, and t2 nth - number of polar tracking angles between Oand "fp nnzz - number of tracks distributed axially between za and zb. For a slab system only one line set card is read, and the values of nner and nnzz on it are redundant. For a cylindrical system the code computes the contribution of a line set to the integrals and looks for another line set card unless the last value of r2_ was the maximum outer radius of the system The number of mesh intervals on the NMESH card must be equal to the number of PRIZE regions defined by the first PRIZE card, B2, is the radial buckling, [om] (or xy buckling in slabs) B2, - is the axial buckling. [om2] No search for critical buckling will be made. n= the number of thermal groups in the fewgroup scheme for the two group leakage edit (default = (number of groups + 1)/2). i= 1 ~ B, flux solution in the leakage edit, or By if all D's on the DNB cards are zero, i= 0 = transport corrected diffusion theory solution in the leakage edit, i=-1- both the above. Default = zero. Buckling values must be set. Zero values will be the default for any material not defined if BEEONE is requested 11 the material number, D1 to D4 - number densities of those nuclides whose P, matrices are on the tape (usually H, D, 0, C or U238 at 293K). AEA TechnologyINPUT WIMSD DIFFUSION TESTTUBES NOCORRELATION ALPHA ikl Note that non-zero D's should be specified only i all the principal nuclides in the material have , matrices. This particularly applies to oxygen in UO, if there is no U238 P, data available. Leakage flux calculations will be made by the following methods: j= 1 -Benoist. j=2- Transport, j= 3 - Ariadne (only for SEQUENCE 2), 4 -Benoist+Transport+Ariadne, 5 - Benoist+Transport, 6 - BenoisttAriadne, j= 7 - Transport+Ariadne. Default j = 2 For the Benoist treatment: If is an integer itis taken to be the number of edit regions forming the outer region of a 3- region cell. fl is not integral, itis taken to be an inner radius of the outer region. This mode of input eliminates the need to know how edit regions are defined ‘The middle Benoist region is defined by k: ~ alr gap with a tube on both sides, - air gap with a tube on the inside, ~ air gap with a tube on the outside, - air gap only, k= 5 - tube only, divided between fuel and moderator, k= 6 - tube only with infinitesimal air gap on the outside, k= 7 - tube only with infinitesimal air gap on the inside, k = 8 - two air gaps with a sleeve between them, k them. - three air gaps with two sleeves between If this card is included the tubes specified by the DIFFUSION card will be smeared with adjacent fuel or moderator in those groups where the tube thickness is greater than the tube material diffusion coefficient. ‘An earlier Benoist prescription without correlation terms will be used. Pu-240 selfshielding can be varied continuously through burnup. i - the first (highest energy) main transport group occupied by the Pu-240 resonance at 1,06eV, J -the lowest such group (default: j=i) AEA Technology The ANSWERS Software Service wimsTESTCASES WIMSD “Test Case 12 continued PRIZE 18.636 19.736 0 1.555 PRIZE 18.44 18.473 0 0.635 4 PRIZE 18.473 18.506 0 0.635 4 PRIZE 19.736 21.49 0 0.635 6 PRIZE 21.49 34 0 0.635 6 PRIZE 34 48.6 0 0.6356 PRIZE 17.68 18.18 0 0.635 7 PRIZE 18.31 18.44 0 0.635 7 PRIZE 18.606 18,636 0 0.635 7 PRIZE 17.68 18.636 0.635 1.9057 PRIZE 18.636 19.736 1.55 1.9057 PRIZE 19.736 21.49 0.635 1.905 7 PRIZE 12.6 15.14 6.35 103 PRIZE 12.6 1.14 10 153 PRIZE 15.14 21.49 6.35 103 PRIZE 15.14 21.49 10 153 PRIZE 21.49 34 0.635 103, PRIZE 21.49 34 10 153 PRIZE 34 43.6 0.635 103 PRIZE 34 43.6 10 153 PRIZE 43.6 48.6 0.635 108 PRIZE 43.6 48.6 10 158 PRIZE 5 6 50 0 0.635 PRIZE 15 6 50 0.635 6.35 PRIZE 5 6 506.35 15 THERMAL 1 NOBUCK PART 27 69 “PART 8 28 BEGIN AEA Technology * NUMBERS OF TRACKS FOR R=0 TO 635 * NUMBERS OF TRACKS FOR R=0.635 TO 6.35 * NUMBERS OF TRACKS FOR R=6.35 TO 15.0 * START OF EDIT DATA The ANSWERS Software Service WIMsWIMSD TESTCASES “Test Case 13, CLUSTER PU NMESH 45 NGROUP 82281 * TESTS ADDITIONAL OUTPUTS (LAMSICITATION/ICAR) NMATERIAL 21 17 NREGION 13 10 40 11 * 11 REGIONS IN SECOND CONDENSATION NRODS 331-1207 1171 POISON 9 PREOUT INITIATE “IAEA-CRP BENCHMARK: HEXAG ASSEMBLY, RADIUS GAD/8 POISON NO BUCK * “TEST FOR ICARME=1 OPTION WITH LIBRARIES ON UNITS 19, 8, 14* FEWGROUP 5 15 27 33 40 45 5 69 MESH 31333311333333 DANCOFF 0.605 0.605 27 BELL 1.1388 MATERIAL 1 -1 11001 2236 9.00E-04 8238.0 2.03E-02 6016 4.24-02 6239.0 0.0 MATERIAL 2-1 6102 91 4.23E-02 MATERIAL 3-1 5753 3001 4.82E-02 6016 2.41E-02 1010 1.06-09 MATERIAL 4 -1 900-1 2235 4.00E-04 8238.0 2.08E-02 6239.0 0.0 6016 4.24E-02 4155 1.92E-04 4157 2.03E-04 MATERIAL 5 -1 600-2 91 4236-02 MATERIAL 6 -1 575 -3 3001 4.82E-02 6016 2.41E-02 MATERIAL 7=4 MATERIAL 8=4 MATERIAL 9=4 MATERIAL 10=4 MATERIAL 1 MATERIAL 1: MATERIAL 1 MATERIAL 1 MATERIAL 1! MATERIAL 16=1 MATERIAL 1 MATERIAL 11 MATERIAL 1! MATERIAL 2 MATERIAL 21= ANNULUS 1 0,386 21 ANNULUS 2 0.455 5 ANNULUS 3 0.6694 6 ANNULUS 41.7711 3 ANNULUS 5 2.9178 3 ANNULUS 6 4.0718 3 ANNULUS 7 5.2282 3 ANNULUS 8 6.3857 3 ANNULUS 9 7.5438 3 ANNULUS 10 8.7022 3 The ANSWERS Software Service wims AEA TechnologyTESTCASES WIMSD “Test Case 13 continued ANNULUS 11 9.8609 3 ANNULUS 12 11.0197 3 ANNULUS 13 12.1787 3 NPIJAN 13 RODSUB 1 1 0.386 1 RODSUB 1 2 0.455 2 RODSUB 2 1 0.386 12 RODSUB 2 2 0.455 2 RODSUB 3 1 0.386 13 RODSUB 3 2 0.455 2 RODSUB 4 1 0.386 14 RODSUB 4 2 0.455 2 RODSUB 5 1 0.386 15 RODSUB 5 2 0.455 2 RODSUB 6 1 0.386 16 RODSUB 6 2 0.455 2 RODSUB 7 1 0.386 17 RODSUB 7 2 0.455 2 RODSUB 8 1 0.386 18 RODSUB 8 20.4552 RODSUB 9 1 0.386 19 RODSUB 9 20.4552 RODSUB 10 1 0.386 20 RODSUB 10 2 0.4552 RODSUB 11 1 0.1576 4 RODSUB 11 2 0.22297 RODSUB 11 3 0.27298 RODSUB 11 4 0.31529 RODSUB 11 5 0.3524 10 RODSUB 11 6 0.3860 11 RODSUB 11 7 0.45505 ARRAY 1 16 1.2203 0.0 ARRAY 2 1 122.3445 0.0 ARRAY 3 1 18 3.4948 0.0 ARRAY 4 16 4.6500 0.2618 16 4.6500 0.5236 16 4.6500 0.7854 ARRAY 5 1 306.8070 0.0 ARRAY 6 1 36 6.9648 0.0 ARRAY 7 1 42 8.1230 0.0 ARRAY 8 1 48 9.2816 0.0 ARRAY 9 1 54 10.4403 0.0 ARRAY 10 1 60 11.5992 0.0 ARRAY 11 1 6 4.6500 0.0 POISON 96 10.1876 41.0 20.2229 71.0 30.2729 81.0 40.3152 91.0 50,3524 10 1.0 60.386 111.0 70.455 51.0 80.6694 31.0 92.9178 120.3325 201295 30.5380 AEA Technology The ANSWERS Software Service wimsTESTCASES WIMSD “Test Case 14, ADDITIONAL OUTPUT ON UNIT 19 AND A USER CHOSEN UNIT (HERE 8) CLUSTER PERSEUS NMESH 27 NGROUP 82280 * ADDITIONAL OUTPUT OPTIONS, ICAR=0 NMATERIAL 6 NREGION 8 4 205 * 5 MOMOD REGIONS IN ADDITIONAL EDIT NREACT 2 PREOUT INITIATE * BENCHMARK IAEA - CRP:ZR-6, SUPERCELL:1+6+12+18+12, HEXAGONAL * “TEST FOR INFINITE LATTICE AND ICARME=0 OPTIONS * FEWGROUP 5 15 27 33 40 45 55 69 MESH6 1134444 SUPPRESS 11011 1111111100 REGULAR 16 MATERIAL 1 -1 2941 2235 8.255E-04 8238.0 2.2105E-02 6016 4.5861E-02 MATERIAL 2-1 2942 91 4.273E-02 1000 4.968E-4 178 6.6E-6 MATERIAL 3-1 294 3 3001 6.6714E-02 6016 3.3357E-02 MATERIAL 4 ~1 294 -2 6016 6.556E-02 27 4,3706E-02 MATERIAL 5 ~1 294 -2 91 4.273E-02 1000 4,968E-04 178 6.6E-06 MATERIAL 6 -1 294 -3 6016 3.3357E-02 3001 6.6714E-02 ANNULUS 1 0.38 4 ANNULUS 2 0.3875 0 ANNULUS 3 0.4525 5 ANNULUS 4 0.6668 6 ANNULUS 5 1.7642 3 ANNULUS 6 2.9065 3 ANNULUS 7 4.0560 3 ANNULUS 8 4.6676 3 RODSUB 1 10.381 RODSUB 1 20.3875 0 RODSUB 1 3 0.4525 2 ARRAY 116 1.210.0 11222700 1183400 1124.210.0 LATTICE * SET OUTER’ DANCOFF = ‘INNER’ BEGINC MOMOD 4 8 12 16 20 * ADDITIONAL SPATIAL SMEAR VECTOR 68 * GROUPS FOR ADDITIONAL EDIT REACTION 2235 294 8238 294 PARTITION 45 69 THERMAL 2 LEAKAGE 5 BUCKLING 0.4989E-02 0.4795E-02 BEGINC The ANSWERS Software Service wims AEA TechnologyWIMSD TESTCASES “Test Case 13 continued SUPPRESS 1111111111 101011 POWERC 1 45.24210.01 BEGINC. THERMAL 2 LEAKAGE 5 BUCKLING 1.0E-07 1.0€-07 VECTOR 6 8 MOMOD 1 3 12 13 15 21 22 37 38.39 40 * ADDITIONAL SPATIAL SMEAR ALPHA 6 BEGINC SUPPRESS 111111414111114 POWERC 1 45.2421 2.21033 10 BEGINC BEGINC *4000MWD/T * BEGINC BEGINC *2000MWD/T* POWERC 1 45.2421 8.84132 10 BEGINC BEGINC SUPPRESS 1111111111 10104 BEGINC BEGINC POWERC 1 45.2421 0.01 BEGINC BEGINC The ANSWERS Software Service WiMsS: AEA TechnologyWIMSD TESTCASES “Test Cas CELL SEQUENCE 4 NCELL 3 * MULTICELL TEST CASE NGROUP 27 NMESH 45, NREGION 9345 NRODS 3150103121 NMATERIAL 6 NTUBE 3 *3 TUBE TYPE REGIONS PREOUT INITIATE * CONTAINMENT WITH DIFFERENT TYPES OF FUEL NESTS BY MULTICELL * 3 TUBULAR FUEL ELEMENTS IN EACH NEST BY PIJ* * FUEL: UO-2 + AL* MATERIAL 1 3.504 293 12236 0.1166 8238 0.2047 27 0.6358 6016 0.04295 15 * CANNING * MATERIAL 2 -1 293 2 27 0.06025 “WATER * MATERIAL 3 0.998 293 3 3001 0.1119 6016 0.8881 wAL" MATERIAL 4 -1 293 -2 27 0.06025 “WATER + AL* MATERIAL 5-1 293 4 3001 6642-02 6016 3.3215E-02 27 4.1€-04 *WODA + AL* MATERIAL 6-1 293 3 3001 6.2666E-02 6016 3.1337E-02 27 3.6786E-03 RODSUB 1 10.3 3 RODSUB2 10.3 3 RODSUB3 10.3 3 RODSUB 1 2 0.3762 RODSUB 2 2 0.3762 RODSUB 3 2 0.3762 RODSUB 1 3 0.4741 RODSUB 2 3 0.4741 RODSUB 3 3 0.4741 RODSUB 1 4 0.55 2 RODSUB 2 4 0.55 2 RODSUB 3 4 0.55 2 RODSUB 1 5 0.853, RODSUB 2 5 0.853 RODSUB 3 5 0.853 RODSUB 1 6 0.9262 RODSUB 2 6 0.9262 RODSUB 3 6 0.9262 RODSUB 1 7 1.0241 RODSUB 2 7 1.0241 RODSUB 3 7 1.0241 RODSUB 1 8 1.1.2 RODSUB2 8 1.1.2 RODSUB 3 8 1.1.2 RODSUB 1 9 1.423 RODSUB2 9 1.423 RODSUB 3 9 1.423 The ANSWERS Software Service wims AEA TechnologyTESTCASES WIMSD “Test Case 15 continued RODSUB 1 10 1.502 RODSUB 2 10 1.502 RODSUB 3 10 1.502 RODSUB 111 1.601 RODSUB2 11 1.601 RODSUB 311 1.601 RODSUB 1 12 1.682 RODSUB 2 12 1.682 RODSUB 3 12 1.682 FEWGROUP 23456789 10 12 14 16 20 22 25 27 40 45 50 53 55 58 60 64 66 68 69 SUPPRESS 1101111111 110011 RTUBE 0.098 0.076 0.3 0.0884 * FUELICAN/COOL THICKNESS + F CELL 15 * 8 CELLS OF TYPE 1 NPIJAN 3 10 ANNULUS 1 4.3 6 ANNULUS 2 4.5 4 ANNULUS 3 6.2785 5 PCELL 10.50.50.0 MESH 555 *3 ANNULI + 12 RODSUBS -> 15 MESHES ARRAY 1 132.0200 CELL 210 NPIJAN 3 10 ANNULUS 1 4.3 3 ANNULUS 2 4.5 4 ANNULUS 3 6.2785 5 PCELL 20.25 0.50.25 MESH 555 ARRAY 2 132.020.0 CELL 3.25 NPIJAN 3 10 CSPECTRUM 2 * USE SPECTROX SOLUTION FROM PREVIOUS CELL ANNULUS 1 4.33 ANNULUS 2 4.5 4 ANNULUS 3 6.2785 5 PCELL 30.00.109 MESH 555 ARRAY 3 132.020.0 BEGINC LEAKAGE 5 BUCKLING 0.1E-07 1.75E-03 NOBUCKLING. THERMAL 9 BEGINC The ANSWERS Software Service wims AEA TechnologyWIMSD TESTCASES “Test Case 16 PRIME 0 SAVE 1164 FUEL-PIN-3. 1WI0-600- SINGLE PWR FUEL PIN ZERO LENGTH BURNUP TO GENERATE OFF-8ASE VALUES DSN 9 GROUP OPTIONS : EQUILIBRIUM BY JOUBERT METHOD : SEGEV INTERPOLATION “TEST 16 PINCELL + SAME CODEWORD FOR CELL 5 OR 6 DSN ENDPRELUDE * ALTERNATIVE TO PREOUT INITIATE ANNULUS1 0.4095 1 ANNULUS2 0.4178 0 ANNULUSS 0.4748 2 ANNULUS4 0.7109 3 MATERIAL 1 -1 978 1 2235 7.2904E-04 8238.0 2.2501E-02 6016 4.6462E-02 6239.0 1.00E-15 MATERIAL 2 -1 607 2 91 4.3153E-02 MATERIAL 3 -1 579 3 3001 4.7742E-02 6016 2.3871E-02 1011 2.3545-05 FEWGROUP 5 11 14 21 27 45 54 61 69 SUPPRESS 1011111111 1010114 * CHAIN 2 PRINTS TEMPERATURES MESH 3 2 2 3 POWERC 1 37.2 0.0 1 * ZERO TIME STEP -> EQUILIBRIUM XE135 SEGEV * ALTERNATIVE RESONANCE INTEGRAL INTERPOLATION EQUILIBRIUM * ALT FORMULA FOR EQUILIBRIUM XE135 DTDR 10.7 COOLANT 3 * CALC FUEL TEMP WITH VARYING RATING BEGINC LEAKAGE 5 THERMAL 3 BUCKLING 1.0E-20 7.018E-05 PARTITION 45 69 REACT 2235 978 8238 978 4135 978 6239 978 BEGINC POWERC 1 46.845 0.01 DTDR 10.7 * NEED TO RESPECIFY DTDR FOR NEW RATING BEGINC BEGINC POWERC 1 58.49 0.01 DTDR 10.7 * NEED TO RESPECIFY DTDR FOR NEW RATING BEGINC BEGINC POWERC 1 66.135 0.01 DTDR 10.7 * NEED TO RESPECIFY DTDR FOR NEW RATING BEGINC BEGINC The ANSWERS Software Service wims AEA TechnologyWIMSD OUTPUT 6.1. General remarks. ‘A full output description has been given for an earlier code version [21]. This section contains updates arising from the many modifications introduced into more recent WIMSD versions. It is worth noticing that the output volume depends strongly on the input options: the output of each chain can be omitted through the SUPPRESS input card; it can also be extended through the input option TEST under which additional information is printed concerning mainly - formulas leading to resonance integrals for a heterogeneous medium, - collision probabilities and matrices of equation coefficients, - history of convergence for iterative solution of equations. The default and TEST option outputs are described in separate sections to facilitate the usage of the present description ("TEST* output is in some chains rather long and confusing and intended, as its name suggests, primarily for debugging purposes). Both types of output are described by Chain: the introductory output produced by the "main" and “prelude” routines is followed by that obtained in consecutive “chains”. The output obtained under ‘specific options (different from TEST) is described together with the default one. ‘The actual printed text is denoted by " “in the description below ‘The first item printed in each chain is the statement: “entry into chain ..” with the given chain umber and the time of the entry. 6.2. Standard output. All equations cross-referenced in this Section are to be found in Appendix 4. INTRODUCTORY OUTPUT List of files opened in "datsets”. Under SAVE option an informative self-explanatory text. Library parameters: number of nuclides, number of energy groups, number of fission spectrum groups, number of fast groups, number of resonance groups, number of fissile nuclides, number of fission product nuclides. ‘A copy of the first section of input data (Prelude Data). ‘Available data storage = length of blank COMMON IQ. Unused storage for each chain Non-essential data storage = the size of excess blank COMMON IQ. Under STORE option, informative control text: “libin entered", “final base address", “libin completed" CHAIN 1 ‘A copy of the second section of input data (Main Data). Under CRAIG option "new annuli" for each cell The ANSWERS Software Service wims AEA TechnologyOUTPUT Under STORE option informative text on library reading: "Libin entere« “Iibin completed”. Under LIBRARY option a map of library data storage. If the CYCLE option has been defined at the next entrance to this chain the array "RECYC(1- 100)". The contents of that array is described as a comment in the subroutine RESTAR where the printout takes place. After the RECYCLE option, in addition to "RECYC( 1-100)" , "recycled plutonium composition”, ie., “number density" and “new fuel composition” Under PLOT option, a geometry illustration of the calculated system (not printed for slab geometry), CHAIN 2 Under STORE option required library storage for dataset 6 in the form: "6 > ....". "Isotope number densities for each material”. First column - a complete list of library nuclide identifiers in the order of library tape, then as many columns as there are materials with material temperature in the first ine and nuclide number densities in respective lines, CHAIN 3 ‘The output bypassed if there is a NORESONANCE card in the input. The output printed for all types of cells (unless the cell is excluded from multigroup calculations by a SPECTRUM input card) in multicell cases. "Volumes" (divided by 2) of fuel, void, can, coolant, air and moderator (the 2 last ones. may be different from zero only in the case of a cluster). In the case of a cylindrical pin cell there are simply r(i+1)?-1()2, for a slab r(it1}-r(), in the case of a cluster the volumes are calculated as_sums of volumes with respective spectrum types divided by the number of rods deciared in the input (for slabs divided by the number of rods +1). The actual number of rods used in that calculation is printed below. in cell radii for ... rods” - radii of successive regions of a representative pin cell calculated according to spectrum types prescribed to materials excluding those with negative spectrum indices. ‘The following output gives results of resonance calculations and is repeated in all resonance groups: Group number. “Total cross-sections for pin cell regions (except fuel)" - omitted in a homogeneous problem. Dancoff factors: inner followed by outer and Bell factor - omitted in a homogeneous problem. Microscopic cross sections for resonance isotopes: for each resonance nuclide cross sections are given for an “array” (infinite lattice) and for a “cluster” (in cluster cases only). The latter are given as inner and “outer” (for the outer ring of fuel rods or plates, cf. Eqs. 1.4 to 1.7 in Appendix 4), - resonance escape probability “up to the volume fraction ... of the coolant" - a fraction of the coolant volume associated with the fuel to increase the slowing down power in correspondence with the increase of the potential scattering in the resonance integral "Flux depression factor” - flux depression due to the resonance absorption, "F(p)" for “fuel, can, coolant, moderator” - the correction to the removal cross section for the flux depression due to resonances, The ANSWERS Software Service wims AEA Technology