Download - Forschungszentrum Jülich
DE06F2007
Forschungszentrum Jülichin der Helmholtz-Gemeinschaft
Institut für Sicherheitsforschung und Reaktortechnik
V.S.O.P. (99/05) Computer Code System
HJ. Rütten, K.A. Haas, H. Brockmann, W. Scherer
Berichte des Forschungszentrums Jülich 4189
V.S.O.P. (99/05) Computer Code System
H. J. Rütten, K. A Haas, H. Brockmann, W. Scherer
Berichte des Forschungszentrums Jülich ; 4189ISSN 0944-2952Institut für Sicherheitsforschung und Reaktortechnik Jül-4189
Zu beziehen durch: Forschungszentrum Jülich GmbH • ZentralbibliothekD-52425 Jülich • Bundesrepublik DeutschlandQ 02461 61-5220 • Telefax: 02461 61-6103 • e-mail: [email protected]
ForschungszentrumJülich GmbH Jül - 4189 October 2005ISR
V. S. O. P. (99/05)
Computer Code System for Reactor Physicsand Fuel Cycle Simulation
by
HJ. Rütten, K.A. Haas,
H. Brockmann, W. Scherer
MINI I*DEO2184 9179*
Abstract
V.S.O.P. is a computer code system for the comprehensive numerical simulation of the physicsof thermal reactors. It implies the setup of the reactor and of the fuel element, processing ofcross sections, neutron spectrum evaluation, neutron diffusion calculation in two or threedimensions, fuel burnup, fuel shuffling, reactor control, thermal hydraulics and fuel cycle costs.The thermal hydraulics part (steady state and time-dependent) is restricted to HTRs and to twospatial dimensions. The code can simulate the reactor operation from the initial core towardsthe equilibrium core.
V.S.O.P.(99 / 05) represents the further development of V.S.O.P. (99). Compared to itsprecursor, the code system has been improved in many details. Major improvements andextensions have been included concerning the neutron spectrum calculation, the 3-d neutrondiffusion options, and the thermal hydraulic section with respect to 'multi-pass'-fuelled pebble-bed cores.
This latest code version was developed and tested under the WINDOWS-XP - operatingsystem. The storage requirement for the executables and the basic libraries associated with thecode amounts to about 15 MB. Another 5 MB are required - if desired - for storage of thesource code (~ 65000 Fortran statements).
11
Preface
The V.S.O.P.- system in its present form is the result of code development, which has
continuously been performed since about three decades. The first edition of the V.S.O.P.-code
was published in 1980 III. Multiple and continuous application for reactor development and
safety research as well as aspects of the fuel cycle of thermal reactors (waste assessment,
proliferation resistance) lead to the demand for various extensions of its capability. In
particular, High Temperature Reactor programs for electricity generation, for nuclear process
heat production and for the incineration of actinides under a non-proliferation aspect gave
important impulses for improvements of the code-system, specially in view of the research on
safety features, reactor operation and licensing requirements. Thus, a continuous upgrading of
the code lead to the code versions V.S.O.P.(94) 111, V.S.O.P.(97) 131 and V.S.O.P.(99) 141
Recent safety research for advanced Modular High Temperature Reactors, primarily in
cooperation with the PBMR / ESKOM-company, gave rise to some major extension of the
code concerning the thermal hydraulics as well as the modeling of neutronics for some special
operating conditions. In the frame of validation and verification a review of the neutron
spectrum calculation methods of the code system was performed, also leading to further
improvements of the respective code sections. The result of these efforts is the code version at
issue, i.e. V.S.O.P.(99/05).
The authors express their thanks to their former colleague Mrs. U. Ohlig, who participated in
updating the code libraries and in code testing and benchmarking, and to the director of the
"Institut für Sicherheitsforschung und Reaktortechnik", Professor K. Kugeler, who promoted
their development work by many suggestions and ideas.
They also express their gratitude for many contributions and suggestions of Mr. F. Reitsma
and his colleagues of PBMR-company concerning the new features of the code system and the
way of their implementation.
in
IV
CONTENTS
Page1. Introduction 1
2. New Features of the Code compared to V.S.O.P.(99) 2
2.1 Neutron Spectrum calculation 22.1.1 Modifications in GAM-I 22.1.2 Modifications in THERMOS 32.1.3 Recommendations for code application 42.1.4 Various additional improvements 4
2.2 Diffusion calculation 4
2.3 Thermal hydraulics 52.3.1 Preparation of heat sources for the thermal hydraulics calculation 52.3.2 Separate calculation of temperatures inside fuel elements for mixtures of
different element types or different numbers of fuel core passages 52.3.3 Heat conductivity of reflector graphite 6
3. Program Organization and Data Base 7
3.1 Structure of the code and program tasks 7
3.2 Nuclear data 123.2.1 Libraries 123.2.2 Identification of the nuclides 133.2.3 Fission products 133.2.4 Preparing a THERMOS-library by means of TTTT 14
3.3 Properties of reactor materials (particularly Graphite) and of the Pebble Bed 21
3.4 List of data sets 24
4. Input Manual V.S.O.P. - MS 26
4.1 Steering the execution mode. SI - S3 26
4.2 Geometric reactor design4.2.1 2-dimensional (r-z - geometry)4.2.2 3-dimensional
4.3 Fuel element design4.3.1 Specifications4.3.2 Design of fuel element-types and -variants
4.3.2.1 Coated particles4.3.2.2 Spherical fuel elements4.3.2.3 Prismatic fuel elements4.3.2.4 Additional nuclides
BI1-BI9TR1-TR5
Dl - D17D1-D4D5-D17D7-D11D12,D13D14-D16D17
282832
35353637394142
4.4 Reactor and fuel cycle4.4.1 Set up dimensions4.4.2 Definition of materials4.4.3 Design and operations
4.4.3.1 Case identification4.4.3.2 Definition of reactor batches4.4.3.3 Data for the burnup calculation4.4.3.4 Control poison search4.4.3.5 Print-out options and steering4.4.3.6 Steering the performance for spectrum and
diffusion calculation4.4.4 Fast and epithermal neutron spectrum4.4.5 Thermal cell spectrum4.4.6 Diffusion calculation
4.4.6.1 Title card4.4.6.2 General control4.4.6.3 Description of neutron flux problem4.4.6.4 Simulation of control devices and void areas4.4.6.5 Fixed source, specified by zones
4.4.7 Fuel cycle costs calculation4.4.8 Fuel management
4.4.8.1 General definitions4.4.8.2 Data for individual fuel types4.4.8.3 Aging boxes for discharged fuel4.4.8.4 Instructions for the burnup cycles4.4.8.5 Criticality search for the reloads4.4.8.6 Redefinition of CITATION edit options4.4.8.7 Extracted nuclides for printout4.4.8.8 'Status of core' - data set for TINTE
4.4.9 Fuel power histogram for decay power evaluation4.4.10 Fuel irradiation histogram for entire isotope generation4.4.11 Preparing THERMOS-library4.4.12 2d-Thermal hydraulics
5. Input Manual V.S.O.P. - ZUT
5.1 Steering the execution mode
5.2 Fuel element design
5.3 Resonance integral calculation5.3.1 Short input5.3.2 Resonance parameters5.3.3 Explicit fuel element design5.3.4 Opening of a new resonance integral data set ('resint')
V1-TX26VIV2-V5V6-V17V6V7-V9V10,VllV12-V14V15
V16, V17G1-G12Tl -T13C1-C21ClC2-C6C7-C10C11-C17C18-C21K1-K12R1-R34R1-R2R3-R4R5R6 - R27R28-R31R32R33R34LF1 - LF3P111 11 -111 15TX1-TX26
ZS
DZ1 - DZ9
Zl -Z17Zl - Z 6Z7 - Z 9Z10-Z16Z17
434344464647495051
51525763636365687072797980828298
100100100101102103105
119
119
119
123123127129132
VI
A. Appendix (Comments) 133
A.I Neutron spectrum calculation 133A. 1.1 Resonance integrals 133A. 1.2 Coated particle grain structure 133A. 1.3 Selfshielding factors in the epithermal energy range 135A. 1.4 Leakage feedback 137
A.2 Design specifications 139A.2.1 Reactor layout (input cards BI and TR) 139A.2.2 Out-of-pile fuel positions 141
A.3 Simulation of reactor operation 143A.4 Restart 143A.5 Fuel cycle costs 144A.6 Thermal hydraulics 146
A.6.1 Structure of the THERMIX code 146A.6.2 Decay power for transient THERMIX calculations 148
A.6.2.1 Power histogram of the fuel batches 148A.6.2.2 Decay power of the fission products 149A.6.2.3 Decay power of 239U and 239Np 150A.6.2.4 Contribution of neutron capture in fission products and
in actinides 151
References 153
Vll
Vlll
List of Figures
Page
Fig. 1 Installation scheme of code executables and libraries 7
Fig. 2: Calculation tasks, V.S.O.P.(99/05) 8
Fig. 3: The basic programs of the two code sections 9
Fig. 4: Built-in fission product chain 20
Fig. 5: Thermal conductivity 21
Fig. 6: Break down of the neutron escape probability 134
Fig. 7: Overlay of reactor-regions and CITATION-material compositions 140
Fig. 8: Out-of-pile fuel positions 142
Fig. 9: Lead and lag times of payments 145
Fig. 10: Coupling between neutronics and thermal hydraulics 146
Fig. 11: Flow scheme of the THERMIX 147
Fig. 12: Neutron capture in fission products and in actinides 152
Tab.Tab.Tab.
Tab.
Tab.
Tab.
Tab.
Tab.
Tab.
I:E:
m:IV:
V:
VI:
VII:
vm:IX:
List of Tables
Page
Available THERMOS libraries 14
GAM-I-library 15
THERMOS-library 16
Sequence of nuclides 17
Fission product yields 18
Available formulae of heat capacity 22
Available formulae of thermal conductivity 23
Alternative specifications of fuel rods 42
Alternative specifications of spherical fuel elements 122
IX
1. Introduction
V.S.O.P. is a computer code system for the comprehensive numerical simulation of the
physics of thermal reactors. The code has widely been used for research in the course of the
development of the High Temperature Reactor with spherical fuel elements. Thus, many tools
have been included to cover very specific features of this reactor type. The application of the
code implies processing of cross sections, the set-up of the reactor and of the fuel element,
neutron spectrum evaluation, neutron diffusion calculation, fuel burnup, fuel shuffling, reactor
control, and thermal hydraulics of steady states and transients. The neutronic calculations can
be performed in up to three dimensions. Thermal hydraulics is restricted to HTRs in two
spatial dimensions.
The V.S.O.P.-code enables the user to follow the reactor life from the initial core towards the
equilibrium core. Repeated calculation of the different physics features ensures consistency in
their feedback during the proceeding burnup, the simulation of the fuel shuffling, and
variations of the core power rating. Temperature transients can be followed doing a quasi-
static nuclear evaluation in parallel. A detailed power history of the fuel elements is used for
the calculation of their individual decay power . Evaluation of fuel cycle costs over the reactor
life time is made using the present worth method. Reprocessing and closure of the fuel cycle
can be simulated under consistent control of the mass flow of the fuel, including the isotopic
decay during periods of intermediate storage. The burnup section of the code also allows the
assessment of the impact of the most important "minor actinides" on the physics of the
reactors and on the characteristics of spent fuel.
Chapters 2 and 3 of this report give a survey of new features of the code, of its structure and
organization. Chapters 4 and 5 are the manuals of the required data input for the code sections
VSOP-MS and VSOP-ZUT, respectively. The APPENDIX contains more detailed
information and comments on special aspects and features of the code.
2. New Features of the Code compared to V.S.O.P.(99)
2.1 Neutron Spectrum calculation
2.1.1 Modifications in GAM-I
In VSOP-99 the GAM-I code /13/ is used in the PI-mode, i.e. for all spectrum zones the Pl-
approximation to the transport equation is used. A leakage term (-B2) may be individually
applied for each (fine) energy group. Once the fluxes (Po) and currents (Pi) have been
calculated they are used to collapse the (microscopic) cross sections of the individual nuclides
to the broad group structure to be used in the subsequent reactor diffusion calculation. All
reaction cross sections and the scattering matrices are collapsed using the fluxes, the transport
cross section is calculated using the currents in the inscatter approximation. The leakage terms
(-B2) are derived from the broad group leakages using the broad group diffusion constants
from the diffusion calculation. Constant leakage terms are used for each fine group within a
broad group.
For strongly negative B2, which may occur in the resonance group of core regions and in the
fast groups of reflector regions, it is observed that sometimes negative fluxes result from the
PI-equations. In these cases previously an adjustment of the flux to a small positive value was
made. As this seems to be quite arbitrary the procedure has been changed in the following
way: If negative or very large (>1020) positive fluxes occur the leakage recycling procedure is
assumed to fail (at least for the present fine group under consideration), the leakage term is
reset to zero and the flux and current are recalculated. Thus maximum leakage information is
retained without generating trouble by non-physical results.
As already mentioned the calculation of the (microscopic) broad group transport cross
sections involves the currents. It is sometimes observed, that this procedure leads to negative
broad-group collapsed transport cross sections for certain nuclides. There are several reasons
for that, like non- adequate PN expansion for the group-to group transfer cross section (e.g. for
hydrogen, where Pi is not adequate) or somehow exotic current spectra caused by use of non-
adequate constant (broad group) leakage terms. In extreme situations this may lead also to a
negative macroscopic transport cross section (and thus diffusion constant) which then will
lead to a failure of the subsequent diffusion calculation.
The modification introduced now checks for negative broad group microscopic transport cross
sections and if necessary replaces them by the (current weighted) total cross section. Thus the
results stay positive although a certain error with respect to the resulting diffusion constant
may remain.
Experience from other codes indicates that a leakage iteration on the basis of DB2 -terms
instead of B2-terms will improve the convergence of that process. As GAM-I uses for B2-
terms a switch-over to DB2-terms needs to define fine group transport cross sections prior to
the flux and current calculation. For this the diagonal-transport- correction has been applied
to obtain the fine group transport cross sections. The main difference of the two approaches is
now, that the use of DB2-terms is equivalent to a variable, fine-group modulated B2-term.
Test calculations have shown the improved stability of this method. However in situations
involving high resonance integrals, e.g. close to the infinitely diluted situation, the
convergence of the leakage iteration process could still be better. It was found, that by
applying a relatively high upper boundary of the resonance broad group (about 130 eV) in e.g.
a 4 group scheme, i.e. by collecting the major resonances of U238 and Th232 in this group, a
complete iteration stability was reached.
Caused by the problems described above, it was recommended for the old V.S.O.P.(99)-
version to use a single, energy-independent homogenised B2-term for all broad epithermal
groups. The modification now allows a leakage iteration process in the GAM-I module for
each individual broad energy group without stability problems.
2.1.2 Modifications in THERMOS
A major computing effort in the THERMOS code /14/ is spent for the calculation of the
transport kernel. In the original THERMOS code as well as in the VSOP(99)- THERMOS
module the integration over space and neutron flight angle is performed by a double
discretisation of space and angle. The either cylindrical or spherical cell is subdivided into
several material zones and each of these again into several spatial meshes. In calculating the
transport kernel by use of first flight escape probabilities, a neutron is assumed to start in the
middle of each mesh in one of about 20 different direction angles. The escape probabilities for
each mesh in outward direction along the chosen flight path are calculated and summed over
the different direction angles. If the outer cell boundary is reached, the remaining items of the
escape probability matrix are calculated using a leakage figure in terms of an albedo derived
from the few-group reactor calculation and making use of a general reciprocal theorem.
It has been experienced, that thereby a relative fine discretisation is necessary especially if
strong heterogeneous structures are considered. Trying e.g. to model a single coated particle
with its part of matrix graphite usually fails, because too few angles will cross the CP-kernel.
The same holds true for boron sticks in graphite.
There is an alternative integration method available which avoids these difficulties without a
larger detriment for the computing time.
In this method a variable transform with respect to the solid angle is made and the integration
is transformed from the angular integration to a 'slice' integration. Making use of a five-point
Newton integration scheme within each 'slice' the accuracy of this method is considerably
better than for the angular integration scheme and the number of spatial mesh points may
significantly be reduced.
This method already used to be applied in a stand-alone version of the THERMOS code /15/.
It is now implemented in the V.S.O.P.(99/05)- THERMOS module.
2.1.3 Recommendations for code application
The review of the neutron spectrum calculation methods in VSOP-99 has shown, that
accuracy and performance improvements are gained by the modifications as described above.
The impact of these on the thermal spectrum is generally small and significant only in very
heterogeneous cell models. In contrast to this, the application of DB2-like leakage terms for
the individual broad energy groups in the fast and epithermal energy range turns out to be a
significant improvement with respect to the leakage iteration stability and to the consistency
of results in various benchmark test. The same holds true for an optimal choice of the energy
range for the third broad energy group in a 4-group scheme for the overall reactor calculation.
The experience during the review process has shown that optimal use of the leakage iteration
method in the nuclear part of VSOP(99/05) is made, if the following recommendations are
followed:
• the leakage iteration process in GAM-I should be applied individually for each broad
energy group.
• the upper energy boundary for the third (resonance) broad group in the standard 4-group
calculation scheme should be around 130 eV.
2.1.4 Various additional improvements
Some further modification have been introduced for spectrum calculation:
• Generally each reactor region also is a spectrum zone.
• A different cell geometry may be defined for different local core positions
(regions), and the THERMOS-cells may be re-defined during the calculation of
proceeding reactor operation.
• Additional scattering matrices of carbon (up to 2200 °C) have been added to the
THERMOS-library (see Table III).
• The ZUT-section of the code has been enabled to treat the resonance integral
calculation for 240Pu.
• The average scattering angle of neutrons for a graphite moderator is now treated as
a function of the moderator temperature.
• The GAM-I-, THERMALIZATION- and THERMOS-libraries are now included
and used in ASCII-Format for easy modification, if required.
2.2 Diffusion calculation
The TRIGIT-section of the code (see chapter 3) has been extended to set up the pattern of
meshes and regions for 3-D neutron diffusion also in cp-r-z - geometry (in addition to x-y-z
geometry). Further, the treatment of voids by means of adapted diffusion constants (see
chapter 4.4.6.4) is now possible also in 3-D-geometry.
2.3 Thermal hydraulics
2.3.1 Preparation of heat sources for the thermal hydraulics calculation
In the last and in older V.S.O.P.-versions point neutron fluxes evaluated in the CITATION-
section first were averaged for VSOP-regions. Based on these flux values, an average power
density was calculated according to the average fission cross section of the region. This
average region value was then transferred - by means of a transfer matrix - to the grid of the
thermal hydraulics section for all those meshes covering the respective region, and the
temperature calculation was performed on the basis of these power density values.
In the present V.S.O.P.-version, the evaluated point power values are transferred directly from
the neutron diffusion section to the thermal hydraulics section of the code in case of a steady
state calculation, using the same grid in both sections as far as the entire reactor model is
covered also by the neutron diffusion calculation. In case of a transient temperature
calculation with zero fission power and the decay power being the only heat source, the
transfer is still done by means of the transfer matrix, because the decay heat is calculated in
code section NAKURE for core batches and regions, respectively (see also sections A.2.1 and
A.6.2).
2.3.2 Separate calculation of temperatures inside fuel elements for mixtures ofdifferentelement types or different numbers of fuel core passages.
In the earlier versions of V.S.O.P. the fuel element temperature, as calculated in the
THERMDC-section of the code, was determined for an "average" fuel element representative
for the local position of the regarded THERMIX-mesh. In case of a "multi-pass" fuelling
concept of a pebble-bed reactor this e.g. meant, that from the local mixture of elements -
having different states of age and burnup and being explicitly treated separately during the
neutronics code section with respect to power and burnup- average values were used for the
thermal hydraulics calculation in terms of power and neutron dose. The temperature
distribution inside the fuel element was then determined on the basis of the average fuel
element power, using heat conductivity values according to the average neutron dose and the
evaluated "average" temperature distribution of the mixture of elements.
In the fuel handling and burnup section of the code partial volumes of different kinds of fuel
elements -forming a "region"- are named "batches". The new approach in V.S.O.P.(99/05)
now is to calculate the temperature of each batch - participating in the power generation of a
THERMIX-mesh - separately. For this purpose the following items are transferred (by means
of a transfer matrix) from the "regiorT-image generally used in V.S.O.P. to the THERMDC-
meshes:
1) the partial power of the batch within the region, which also becomes the partial power of
the respective fuel element type in the THERMDC-mesh;
2) the partial volumes of the batches of the regions, which the THERMEX-mesh belongs to;
3) The dose of fast neutrons for each batch appearing in the mesh.
The temperature in each type of fuel element of each mesh is then calculated according to the
individual fuel element power and the individual neutron dose, i.e. the so called
"heterogeneous" THERMDC-calculation, which used to be applied for the "average" fuel
element of a mesh, is now performed as many times for each iteration, as is the number of fuel
element types mixed within the pebble-bed. The different surface temperatures of the
respective fuel types in a mesh are averaged to form the "solid temperature" to be used for the
next calculation of the gas temperature during the iteration process.
The final values of the fuel temperature and of the moderator temperature at different
positions within the elements of different types are also averaged for all the meshes belonging
to a VSOP-region to be used for the next calculation of the neutron spectrum.
2.3.3 Heat conductivity of reflector graphite
The heat conductivity of the reflector graphite not only is a function of the actual temperature
and of the accumulated neutron dose, but also of the temperature which occurred during the
irradiation of the graphite. With respect to this, the reflector dose is calculated from the
evaluated fast neutron flux and a user-defined reflector lifetime. The temperature field may be
retrieved from dataset 'tempstat', being the result of a previous calculation. Thus, the reflector
temperature and the graphite conductivity may be iterated.
2.4 Some additional, general extensions.
• A variable number of fuel batches per region in different radial flow channels is now
allowed (see chapters 4.2.1 and A.2.1).
• A comfortable data input procedure for the definition of a water ingress into the core of
an HTR has been provided.
• A redefinition of the cost data input in a restart is now possible (for parametric studies)
• Data sets "rstlib" and "rstnew" have been combined, (dataset "rstlib" has been dropped)
• Generally, the data input procedure had to be modified in several places, either to
account for code extensions or in order to make the input structure even more clearer.
So, input data sets, which used to be applied to use VSOP(99) should be thoroughly
revised for use with VSOP(99/05).
3. Program Organization and Data Base
3.1 Structure of the code and program tasks
The Fortran source of the code system is structured in two parts to be used for the construction
of two executable codes by the compilation and linking procedure. The first executable is
named VSOP-ZUT. This code section is used to set up a matrix of resonance cross sections to
be used in following applications of the MAIN-section of V.S.O.P. The latter code section is
available as the executable code VSOP-MS. The two executables should be positioned within
the same directory of the permanent memory device. We name this directory the "Working
Directory" (see Fig. 1). The working directory must have a subdirectory, which must be
named "Libraries", containing the various library data files. Necessary read- and write- data
transfer from and to these data sets - illustrated in Fig.l- is done automatically during the
execution of the code via this path. The name (and possibly the path) of the card image input
and of the printout data sets (this sequence) have to be given as arguments when starting the
program. All other data sets -unless they are automatically scratched after the end of code
execution- are written onto the working directory and they are available for further use ( e.g.
Restart-libraries, data to be graphically displayed etc.). A list of the data sets possibly
produced or used is given in chapter 3.4. Storage requirement for the executables and the
basic libraries associated with the code amounts to about 15 MB. The source codes consist of
about 65000 Fortran statements. They have been compiled, linked and used under the
WINDOWS-XP operating system.
Working Directory(arbitrary name)
VSOP(99/05)-ZUT.EXEVSOP(99/05)-MS.EXE
Subdirectory ^^~^*"Libraries"
Subdirectory:Libraries
(obligatory name)
adageresdatthresdatu8resdapuOresdapu2
resintgam_99_4
thermaltherml515therml516therml517
(therml518)(therml519)
Fig. 1: Installation scheme of code executables and libraries
User-data input processing
Resonanceintegrallibrary
Fuel element design
Resonance integrals
User-data input processing
Core design
Fuel element design
Neutron spectrum
Neutron diffusion(2- or 3-d flux distribution, reactor control)
Burnup
Thermal hydraulics(static or transient)
Fuel managementfuel cycle costs
N
o
*
O
Power histogram
: MS = Main Section
Fig. 2: Calculation tasks, V.S.O.P.(99/05)
Resonanceintegrallibrary
ZDATA-2
ZUT
BIRGIT / TRIGIT
DATA-2
Z^LGAM - 1 , THERMOS
CITATION
A.ADAGE, BURNUP
THERMIX
FUMAN, KUGELN, KOSTPW
HDSJsiO
siO
Externalcodes
Fig. 3: The basic programs of the two code sections
The structure of the V.S.O.P.(99/05) - code system in terms of calculation tasks on the one
hand and in terms of the programs and basic libraries on the other hand is illustrated in Fig. 2
and in Fig. 3, respectively.
Subroutines DATA-2 and ZDATA-2 use input data specifying the fuel element design to
prepare detailed data to be used in different subroutines of V.S.O.P.
The BIRGIT code prepares the mesh pattern for a 2-d diffusion calculation and for the thermal
hydraulics calculation. If experimental data for the flow of the fuel elements in a Pebble-bed
reactor are available to the user, the code can generate a sophisticated flow pattern of finite
batches of elements which move down to the discharge tubes in finite steps. It then provides
the transformation between this pattern and the one used in the diffusion calculation for the
macroscopic cross sections and for the neutron flux. Similar transformations may be
performed for decay heat values used in the thermal hydraulics calculations.
For a 3-d diffusion and burnup calculation the mesh pattern is prepared by use of subroutine
TRIGIT.
The neutron spectrum is evaluated by means of GAM-I /13/ and THERMOS /14,15,16/ for all
the different reactor material regions, providing broad group cross sections for neutron
absorption, neutron-induced fission and n,2n-reactions. The epithermal cross section library
has a 68 group structure, corresponding to the GAM-I code. The thermal library fits in the 30
group structure of THERMOS. These libraries have been generated using ENDF/B-PV, -V and
JEF-I data 111. For some nuclides the self-shielded cross sections within the resonance energy
range are calculated prior to a VSOP-MS calculation by means of the VSOP-ZUT code
section, which is based on the ZUT-DGL code /8,9/. This is done for 232Th , 238U ' 240Pu and242Pu using resonance data also from ENDF/B-IV and -V. The shielded resonance cross
sections are stored on a permanent data set (see also chapter 3.4).
Two different sets of temperature dependent cross sections may be generated and used for
each resonance absorber, both representing different absorber concentrations. Normally, these
two concentrations should represent the highest and the lowest concentration, respectively,
occurring in the fuel elements. The code then first constructs new resonance cross sections for
these two states of burnup by interpolation between the values for different temperatures
according to the calculated fuel temperature in each region. It then gains the final cross
sections by interpolation between these two sets for the highest concentration on the one hand
and the lowest concentration on the other hand according to the true absorber concentration.
Graphite scattering matrices are based on the Young -Koppel phonon spectrum in Graphite
/10.11.12/.
Neutron diffusion is calculated by CITATION /III in its 2- or 3- dimensional version.
10
The burnup of -at present- 28 heavy metal isotopes in the numerous fuel batches (up to 9999)
is evaluated by means of subroutine ADAGE. It includes the subroutines TERM, DECAY,
EQUIL and parts of subroutines NUDATA and FLUXO, as they have been used in ORIGEN-
JUEL-II151, which in its turn has been based on the ORIGEN code 161. ADAGE uses the
matrix exponential method in order to treat the decay and transition scheme. For this purpose
it constructs a transition matrix for the group of heavy metals, whose size is determined by the
number of isotopes defined in the ADAGE-library. This library contains an identification
number for each nuclide to be treated. It also contains the decay constants and information
about the possible nuclide transitions by a, ß\ ß+ - decay, by disintegration from an excited
nuclear state to the ground state and by spontaneous fission. At present, the ADAGE-library
defines the heavy metal chain from 232Th through 244Cm according to the cross section data
available in the GAM-I-library (chapter 3.2). As burnup equations now are no longer fixed
within the code, but the transition matrix is constructed according to the information contained
in the ADAGE-library, the number of the considered materials in principle can easily be
extended.
The burnup of the fission products and the simulation of the fuel shuffling are covered by
subroutines BURNUP, FUMAN and KUGELN. These subroutines are a further development
of the FEVER code /18/. A burnup chain for 44 fission products is included, which can
optionally be supplemented. Collapsing the cross sections into broader energy groups can be
made for 171 - at present- nuclides of the source libraries, as desired for a more detailed
calculation of the isotopic composition of the fuel by means of the code ORIGEN-JUEL-II151.
The THERMDC code /19,22,23/ is included as program part for static and for time dependent
thermal hydraulics. The resulting temperature values of the fuel and of the moderator regions
are fed back to the spectrum calculation in each spectrum zone for subsequent core neutronics
calculations.
The economics code KPD I20I performs the fuel cycle cost evaluation on the basis of the
"present worth" method and it is included as subroutine KOSTPW.
The status of the reactor at the end of each calculation can be preserved for a restart of the
code. Further preservation of calculated reactor data may be provided for the purpose of joint
evaluations beyond the capability of VSOP. The full power history of the fuel batches is
preserved for the calculation of the local and of the integral decay power of the fuel by means
of the NAKURE code /21/, which is included in VSOP as subroutine NACHW and which is
needed e.g. for the explicit following of reactor heat-up under loss-of-coolant conditions.
Internal restart facilities have been included for THERMOS and CITATION: After
convergence of the first calculation the neutron flux fields are preserved as initial values for
subsequent calculations, which results in a efficient reduction of the computing time.
11
The simulation of the fuel bumup is an alternation between the simulation of the fuel
shuffling, the evaluation of the space and energy dependent neutron flux and the resulting
local depletion and generation of the nuclides. In the nomenclature of the code the phase
between two shuffling steps is named "Burnup Cycle". It is subdivided into "large burnup time
steps". At its beginning the spectrum calculations and the diffusion calculation can be
repeated. All flux values then are re-normalized to the demanded power level of the reactor.
The large burnup time steps are subdivided again into "small burnup time steps". According to
the proceeding burnup of the fuel, the height of the neutron flux is re-adjusted to the core
power for each of these steps, whereas its spatial distribution remains unchanged till the start
of the next large burnup time step.
For each small burnup time step, the subroutines ADAGE and BURNUP are called in order to
solve the burnup equations for the fuel elements, using the re-normalized neutron flux at its
local position.
3.2 Nuclear data
3.2.1 Libraries
Neutron spectrum calculations are performed using the GAM-I IX3/ and the THERMOS
/14,15/ code. The two respective libraries have been derived from the evaluated nuclear data
files ENDF/B-V and JEF-1. The GAM-I-library data is given in an 68-energy group structure
ranging from 10 MeV through 0.414 eV. It contains cross section data for 171 materials. In
addition, there are 6 more heavy metal isotopes included (Id.-no. 185 through 190), which are
part of the actinide build-up chain. These are short-lived isotopes, and only decay constants
but no cross section data have been available to the authors for the time being. Their half-
lives, however, are short enough (some minutes through 2 days) to dominate neutron capture
rates by far. Thus, neutron capture may well be neglected. A list of the isotopes and of the
origin of the nuclear data is given in Table II.
The THERMOS-library data are given in 30 energy groups ranging from 10"5 through 2.05
eV. The library (Tab. Ill) is subdivided into 2 parts: (1) The absorbers with identification
numbers identical to those in the GAM-I-library. (2) The scatterers with identification
numbers made of 4 digits. For the scattering nuclides scattering kernels have formerly been
constructed applying different scattering laws and different temperatures. Tab. Ill gives the
respective information. The basic thermal library provides data in the form of the 96-energy
group structure of the THERMALIZATION spectrum code. It has to be condensed to a
THERMOS library using a neutron spectrum which is representative for the considered
reactor. This may be done by use of the program part TTTT (see section 3.2.4).
12
Some data concerning disintegration of the heavy metal isotopes by alpha emission, beta
decay transitions, positron emission, spontaneous fission and isomeric transitions as well as
their half-lives are provided by the ADAGE-library, which is used in the burnup section of
the code.
3.2.2 Identification of the nuclides
The number of materials used in the VSOP-calculation is presently limited to 200. This is
caused by the dimensioning of some data sets and may easily be extended, if required. The
identification of each nuclide is defined by input on cards D2 (Section 4.3.1). Each nuclide of
a VSOP-calculation is defined by the number IMAT(I), which is its Id.-no. in the GAM-I-
library, as listed in Table Ü. A certain sequence must be observed for the designation of the
nuclides, which is outlined in Table IV: The heavy metal isotopes (28 at present) are firmly
assigned. They are followed by the fission products, with 135Xe and the cumulative fission
product in the positions 29 and 30, respectively. The order of the explicitly treated fission
products must correspond with the chain definition (Fig. 4). Subsequently, control poisons
represent absorbers of variable concentration. They are followed by absorbers of fixed
concentrations. The scattering nuclides must be given at the end. Here, they are also identified
by their
GAM-I-Id.- numbers. Note that the code accepts several scattering matrices from the
THERMOS-library for one and the same scattering nuclide, each for a different temperature of
the scattering material. This information is given in the THERMOS data input.
3.2.3 Fission products
Yields are included in the code for 97 fission products (Tab. V). They have been taken partly
from ENDF/B-IV, partly from ENDF/B-V. The built-in fission product chain and the
sequence of the identification numbers are given in Fig. 4. The second number of this chain is
a "non-saturating" fission product. It stands for the sum of many low absorbing fission
products which are not included in the chain. The yields of the non saturating fission product
were adjusted by comparison with results generated by the ORIGEN-JÜL-Ü code 151 which
comprises more than 800 fission products explicitly. In a typical HTR with a fuel burnup
equal 80 MWd/kg HM the explicit fission products of the VSOP-chain were found to cover
98.02 % of the total neutron absorptions by fission products obtained by the ORIGEN-
calculation. The yields of the accumulative fission products have been adapted to cover the
remaining 1.98 %.
It is possible to extend the fission product chain by defining new isotopes, new yields and new
chain information on cards V4, V5. Similarly, the chain can be shortened, modified or even be
fully replaced by the user of the code.
13
3.2.4 Preparing a THERMOS-library by means of TTTT
The basic thermal library of VSOP(99/05) is given in 96 thermal energy groups ranging from
10"5 eV through 2.05 eV. It is made in the structure of the zero-dimensional thermal spectrum
code THERMALIZATION. For normal use within VSOP this code has been replaced by the
THERMOS code performing thermal cell calculations in one dimension and in 30 energy
groups.
THERMOS requires a specific 30 group library. This can be generated by condensing the 96
groups THERMALIZATION-library with the neutron flux belonging to the considered
problem. For this purpose the VSOP input must be prepared for the respective reactor design
case with one representative spectrum calculation. On input card Tl a blank space for variable
CIDTHER tells the code to run THERMALIZATION instead of THERMOS, and the thermal
neutron spectrum is preserved for the condensation of the cross sections (see input description
of sections 4.4.5 and 4.4.11).
Table I: Available THERMOS libraries
Data set name
therml515
therml516
therml517
therml518
therm!519
Neutron spectrum typical for:
HTR
LWR
HWR
*
*
* reserved for additional libraries to be generated by the user, if desired
14
Table II: GAM-I-library
Id.
1
2
3
4
5
e7
B9
10
11
12
13
14
IS16
17
22
23
24
25
26
27
2829
30
31
32
33
34
3536
37
38
39
4041
42
43
44
45
46
47
4849
50
51
52
53
54
55
5657
58
59
60
61
62
63
-no.
H-l
H-2Be-9B(nat)CTh-232Pa-233U-233U-234U-235U-236U-238Np-239Pu-239Pu-240Pu-241Pu-242N-14O-16MgAl-27SiCrMn-55Fe(nat)Co-59NiCuSe-82Br-81Kr-83Kr-84Kr-85Kr-86Rb-85Rb-87Sr-88Sr-90Y-89ZrZr-90Zr-91Zr-92Zr-93Zr-94Zr-96MoMo-95Mo-96Mo-97Mo-98Mo-100Tc-99Ru-100Ru-101Ru-102Ru-104Rh-103Pd-104
MatMatMat
MatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMat
MatMatMatMatMatMatMatMatMatMatMatMatMatMat
1301
4012
1289
1306
4902
1391
4923
4924
4925
4926
4928
4939
1264
4 94 0
4 941
1342
4074
4086
4120
4137
4140
4240
4255
4260
4279
4280
4290
4342
4351
4363
4364
4365
4366
4375
4377
4388
4380
4399
4409
4400
4401
4402
44034404
4406
4420
4425
9283
4427
9285
9287
1308
4440
4441
4442
4444
1310
4464
ENDFB-VJEF-1ENDFB-4JEF-1ENDFB-VJEF-1ENDFB-VJEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-4JEF-1JEF-1ENDFB-VJEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-VJEF-1ENDFB-VENDFB-VENDFB-VJEF-1JEF-1JEF-1JEF-1ENDFB-VJEF-1
Id.
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
8687
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
-no.
Pd-105Pd-106Pd-107Pd-108Pd-110Ag-109In-115CdCd-110Cd-111Cd-112Cd-113Cd-114Te-126Te-128Te-1301-1271-129Xe-128Xe-130Xe-131Xe-132Xe-134Xe-135Xe-136Ce-133Cs-135Cs-137Ba-134Ba-136Ba-137Ba-138La-139Ce-140Ce-142Pr-141Nd-142Nd-143Nd-144Nd-145Nd-146Nd-148Nd-150Pm-147Sm-147Sm-148Sm-149Sm-150Sm-151Sm-152Sm-154Eu-151Eu-153Eu-154Eu-155Gd-154Gd-155Gd-156Gd-157
MatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMatMat
446544664467938644601373449544804483448444B5448644874525452745299606960845424544454545464548454945514553455596694564456645674568970745804582974297634603460497664606976946009783980698071319980946219811981346314633463498324644464546464647
JEF-1JEF-1JEF-1ENDFB-VJEF-1ENDFB-VJEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-VENDFB-VJEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-VJEF-1JEF-1JEF-1JEF-1ENDFB-VJEF-1JEF-1ENDFB-VENDFB-VJEF-1JEF-1ENDFB-VJEF-1ENDFB-VJEF-1ENDFB-VENDFB-VENDFB-VENDFB-VENDFB-VJEF-1ENDFB-VENDFB-VJEF-1JEF-1JEF-1ENDFB-VJEF-1JEF-1JEF-1JEF-1
Id.
123
124
125
126127
128
129
130
132
133
134
135
136
137
139
140
141
142
143
144
145
146
147
148
149
151
152
153
154
155
156
160
164
165
166
167
168
169
170
171
172
173
174
175176177
178
179
180
181
182
183
184
185
186187
188
189
190
-no.
Gd-158Tb-159Au-197PbBi-209Li-6Li-7B-10U-237Np-237
MatMatMatMatMatMatMatMatMatMat
4648465947974820483940364037405049274937
JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1
V mit Cr-Streumatr. GERV
Nb mitTiAg-107Nb-93W(nat)Ru-105Rh-105Cs-134Ce-144Pr-142Pm-148Pm-148mZr-95Ru-103Xe-133Ce-141Pr-143Pm-1491-131
Mat 4230 JEF-1Zr-Streumatr. GERMatMatMat
MatMatMatMatMatMatMatMatMatMatMatMatMatMat
422044774413
44454455455445844592461246134405444345479725459346144536
JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-VJEF-1JEF-1JEF-1
Accum.fiss.prod.chain 44B-llHf-174Hf-176Hf-177Hf-178Hf-179Hf-180W-182W-183W-184W-186Pm-151U-232Pu-238Am-241Am-242Am-242mAm-243Cm-242Cm-243Cm-244Th-233U -239Np-238Np-240Pu-243Am-244
MatMatMatMatMatMatMatMatMat
MatMatMatMat
MatMat
MatMatMatMatMatMat
405147244726472747284729472047424743474447464615823213381361854213691363864213431344
JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1JEF-1ENDFB-VENDFB-VENDFB-VENDFB-VENDFB-VENDFB-VENDFB-VENDFB-VENDFB-V0-SIGMA0-SIGMA0-SIGMA0-SIGMA0-SIGMA0-SIGMA
15
Table III: THERMOS-library
Id.-no.
4
c.D
22
24
190
Id.-no.10121013101410221023
11011102110311041105
11111112111311141115
11211122112311241125112611271128
1600160116021603160416051606160716081609161016111612.1613161416151616
AbsorberH-l 1H-2 I see Scatterer
Be-9 JB(nat) JEF-1C see Scatterer
1I
1 see GAM-
J[-Library
0-16 see Scatterer
11 see GAM-'
JScattererBerylliumBerylliumBerylliumBeryllium
980
13661422
=>K
3Kin BeO
Berylliumoxyd
HydrogenHydrogenHydrogenHydrogenHydrogen
DeuteriumDeuteriumDeuteriumDeuteriumDeuterium
OxygenOxygenOxygenOxygenOxygenOxygenOxygenOxygen
CarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbonCarbon
293
323
373
473
573
293
323
373
473
573
293
323
373
473
573
900
12001350
300
400
500
600
700'800'900'
1000'1100'1200'1300'1350'1500'
[-library
P0 & PI Gas kernel 1/V Sigma0=10mbP0 & PI Gas kernel 2000F 1/V Sigma0=10mbP0 & PI Gas kernel 2100F 1/V Sigma0=10mb900oK Summit-BeO-Matr. minus O-Matr. Rest von 1012
900°K Summit-BeO-Matr. XA=0.01/V XS=9.69 S=S(O) & S(Be)
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
.6
>K
>K'K
'KK
K
KK
KK
K
KKKKK
1650°K1800°K2000°k2200°K
°K Nelkin-Kern (Gaker-Kira) IX.68 Darvas3K Nelkin-Kern (Gaker-Kira) IX.68 Darvas=K Nelkin-Kern (Gaker-Kira) IX.68 Darvas3K Nelkin-Kern (Gaker-Kira) IX.68 Darvas=K Nelkin-Kern (Gaker-Kira) IX.68 Darvas
3K Nelkin-Kern (Gaker-Kira) IX.68 Darvas°K Nelkin-Kern (Gaker-Kira) IX.68 Darvas>K Nelkin-Kern (Gaker-Kira) IX.68 Darvas'K Nelkin-Kern (Gaker-Kira) IX.68 Darvas>K Nelkin-Kern (Gaker-Kira) IX.68 Darvas
SK Brown-St-John-Freigas IX.68 Teuchert5K Brown-St-John-Freigas IX.68 Teuchert5K Brown-St-John-Freigas IX.68 Teuchert5K Brown-St-John-Freigas IX.68 Teuchert'K Brown-St-John-Freigas IX.68 Teuchert
Brown-St-John-Freigas XII.70 TeuchertBrown-St-John-Freigas XII.70 TeuchertBrown-St-John-Freigas XII.70 Teuchert
Young Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eVYoung Phon.-Spektr. Colli Punktwerte 2eV
SchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroederSchroeder
Brown-St-John-Freigas Haas / Rütten 10Brown-St-John-Freigas Haas / Rütten 10Brown-St-John-Freigas Haas / Rütten 10Brown-St-John-Freigas Haas / Rütten 10
4 .4 .4 .4 .4 .4.
4 .4 .4 .4 .4.
4 .4 .. 1.1
. 1
. 1
7
7
7
7
7
7
7
7
7
7
7
7
7
70
70
70
70
70
70
70
70
70
70
70
70
70
2002200220022002
16
Table IV: Sequence of nuclides
1.
2.
3.
4 .
5.
VSOP-Id.no. GAM-I-Id.no.
Heavy metal isotopes are firmly assigned:
123456789
10111213141516171819202122232425262728
Th-232Th-233Pa-233U -233U -234U -235U -236U -237U -238U -239Np-237Np-238Np-239Np-240Pu-238Pu-239Pu-240Pu-241Pu-242Pu-243Am-241Am-242mAm-242Am-243Am-244Cm-242Cm-243Cm-244
6185789
101113212
18613318713
18817714151617
189178180179181190182183184
Fission products of chain definition:
293031 ...
Xe-135Accumulative fission productFurther isotopes of the chainNO < 49 fission products are allowed
Control poison:
subsequent NC = 0-2 different nuclides are pos-sible with concentrationsadjustable to achieve given
Non burning absorbers:
subsequent Absorbers for which concentrationsremain unchanged during burnup, e.g.structural materials
Scatterers:
subsequent NKER = 1-5 scatterers must be givenat the end.
17
Table V: Fission product yields (values given in percentages)
Isotope
Se- 82Br- 81Kr- 83Kr- 84Kr- 8 5Kr- 8 6Rb- 85Rb- 87Sr- 88Sr- 90Y - 89Zr(nat)Zr- 90Zr- 91Zr- 92Zr- 93Zr- 94Zr- 95Zr- 96Mo- 95Mo- 96Mo- 97Mo- 98Mo-100Tc- 99Ru -101Ru-102Ru-103Ru-104Ru-105Rh-103Rh-105Pd-105Pd-106Pd-107Pd-108Pd-110Ag-10 9Cd-111Cd-112Cd -113Cd-114In-115Te-126Te-128Te-130I -127I -129I -131I -135Xe-131Xe-132Xe-133Xe-134Xe-13 5Xe-135Xe-136Cs-133Cs-134Cs-135Cs-137Ba-138
Type
ttttccittct
it
tccctcitttcttct
ici
ctttttct
tcctcccitccicciicct
U -233
0.56262
0.31171
1.0178
1.7034
2.1946
2.8581
6.5296-5
4.0088
5.4953
6.7952
6.2568
6.4467
0.05
6.5194
6.5949
7.011
6.8076
6.2478
5.6694
9.5909-4
6.5-3
5.4533
5.1587
4.4094
4.9573
3.2258
2.4492
1.7066
1.0276
0.48
1.4219-9
0.47126
3 .4998-11
0.24063
0.11417
0.061481
0.025376
0.043363
0 .020268
0 .014602
0.013152
0.012268
0 . 020052
0.24081
0.94592
2.3671
0.67853
1 .616
3 .70894.8597
8.4795-5
4.80386.0307
5.7588
1.3374
6.1971
6.7934
3.6998-5
1.1969-3
6.1
6.7889
5.8863
U -235
0.33405
0.21005
0.53076
0.98786
1 .314
1.9528
8.23-5
2 .551
3 .6228
5.9137
4.8469
5.803
0.047
5.926
5.966
6.3703
6.4228
6.4678
6.2506
1.641-4
5.85-4
5.965.7787
6.3096
6.1284
5.0501
4.2032
3.1411
1.8239
0.9
1.858-9
1.0199
9.83-11
0.37759
0.16317
0.071032
0.022338
0.029903
0.019714
0.012802
0.012425
0.011256
9.9367-3
0.057818
0.35046
1.44660 .13037
0.65911
2.8325
6.3482
1.54-6
4.2498
6.7859
7.6825
0.2541
6.6023
6.2701
5.08-5
3.57-5
6.45
6.269
6.8272
0
0
0
0
0
0
5
01
2
1
2
0
2
3
3
4
4
5
1
7
5
5
6
6
5
6
6
5
5
1
5
2
4
3
2
0
1
0
0
0
0
0
0
0
2
0
1
3
6
1
5
67
1
7
6
1
4
7
6
5
Pu-239
.21092
.1768
.29608
.48029
.56834
.75863
.85-5
.94936
.3703
.1134
.7075
.6405
. 0164
.4941
.018
.9031
.4431
. 9212
.0958
.492-3
.7-4
.608
.8542
.977
.1405
. 9135
.0201
.9845
.9539
.47
.358-7
.4261
.03-8
.6234
.2361
.2319
.62204
.4115
.27428
.10707
.078216
.046789
.040467
.19996
.85079
.4971
.49173
.5039
.738
. 3007
.652-5
.2688
.9758
.389
. 1517
.4524
.6153
.61-5
.61-4
.22
.6834
.7173
0
0
0
0
0
0
5
0
01
1
2
0
1
2
2
34
4
1
7
4
5
6
66
6
6
6
5
5
61
4
5
4
1
2
0
0
0
0
0
0
01
0
0
3
6
1
4
6
8
07
7
4
3
7
6
6
continued
Pu-241
.11602
.06469
.20498
. 35393
. 39618
.61392
.3024-7
.75709
.97473
.5363
.2146
.645
.0164
.8315
.2781
.9643
.4018
.0456
.4232
.2927-5
. 7-4
.8208
.2217
.2311
.2085
.0948
.4843
.2611
.9764
.47
.6028-5
.2183
.6908-6
.6314
.3339
.0191
.2091
.2836
.57261
.23001
. 15494
.075514
.040537
.077127
.35555
.6617
.23046
.77864
.1411
.95
.3066-6
.6411
.741
.1081
.22923
. 1792
.2871
.302-7
.5416-5
.8
.698
.4446
18
La-139Ce-14 0Ce-141Ce-142Ce-144Pr-141Pr-143Nd-142Nd-14 3Nd-144Nd-145Nd-146Nd-148Nd-150Pm-147Pm-148mPtn-148gPm-149Pm-151Sm-147Sm-148Sm-149Sm-150Sm-151Sm-152Sm-154Eu-153Eu-154Eu-155Gd-154Gd-155Gd-156Gd-157Gd-158Tb-159
Continuation of Fission Product yields
tt
t
tc
ittttccii
ii
c
tttic
itttt
ict
5.8856.43346.246.63044.51176.62245.85130.2.4799-84.64953.42482.59731.28670.498461.77532.7899-59.4395-70.769530.322932.0099-101.7999-80.2.5782-30.0.207840.045580.106863.7198-50.0212520.3.6198-70.0117376.7747-32.2298-39.2311-4
6.49336.32295.735.92475.9625.89295.9710.0099.5-115.45233.93392.99121.690.645932.27017.49-75.73-61.08880.420440.6.95-110.5.413-40.0.270570.0746890.162641.63-60.0330250.4.41-90.0135176.4651-33.2163-31.0394-3
Independent fission= Cumulative fission
Total chain yield
Accumulative Fission
U-233
112.3
Product:
U-235
94.76
5.64565.57516.115.01734.45145.36344.56130.00094.9-103.8343.08332.53331.69820.994512.07692.09-62.09-61.26170.77722.43-122.8-100.1.7009-30.0.596180.276820.372243.54-50.170820.2.83-70.119890.0762970.0409550.021205
yieldyield
Pu-239
115.6
6 .22835.8946 .114.8154.86444.85344.50170.00091.2106-104.15643.20462.74011.92571.1962.26015.4125-75.4125-71.46350.902380.4.272-110.3.9438-40.0.717040.379790.528155.5626-60.231810.1.9109-80.169550.131530.0867070.046741
110.4
19
Fission Product Chain
Xe-135
I
Xe-136
FP-44
Kr-83
Zr-95
Uo-97
"c-99
Ru-tOI
Ru-103-
Rh-105-
Pd-108s
Cd-113
1-131 -
Xt-133-
Pr-14'
Pr-H3
(Iron 1-135)
• Mo-95
• Rh-103
• Pd-105
• Ag-109
' Xe-lSl
• Cs-133I
Cs-134
N<J-143
I
Nd-145I
Nd-146
••Pm-147 • Sm-147
I• Pm-I*8m N
I \• Pm-I48g , Sm-148
1 I• Pm-149 • Sm-149
ISm-150
IPm-151 • Sm-151
1Srr-152
• Eu-153
Fig.4: Built-in fission product chain
135 Xe
FP-44136 Xe83 Kr95 Zr95 Mo97 Mo99 Tc
101 Ru103 Ru103 Rh105 Rh105 p d
108 p d
109 Ag113 Cd13, ,
131 Xe133 Xe133 Cs134 Cs141 Pr.43 p r
143 Nd144 Nd145 Nd146 Nd147 Pm1411 Pm-m148 Pm-g147 Sm148 Sm149 Pm149 Sm150 Sm1M Pm151 Sm152 Sm153 Eu154 Eu155 Eu
' " Gd156 Gd157 Gd
VSOP-No
29
3031
323334
35
3637
3839
4041
42
4344
45
4647
484950
5152
5354
55
5657
585960
61
62
6364
65
666768
69
7071
72
GAM-No
87
16088
35149
5254
57
59
151
62
14364
6769
75
15684
15289
144
99154
101
102
103104107
148147
108109
155
110111
175112
113116117
118
120121
122
20
3.3 Properties of reactor materials (particularly Graphite) and of thePebble Bed
The thermal conductivity A. of graphite is a function of four parameters:
1. The type of graphite material due to its fabrication techniques
2. Its temperature T
3. The exposure to fast neutrons D
4. The temperature at which this exposure occurred
Broad experimental research on this matter has been performed by L. Binkele /24/.
As an example Fig. 5 (left part) gives X,(T,D) for the NUKEM/A3-3 graphite having been
exposed to fast neutrons at the temperature 950 °C. Measurements were made at ten different
temperature values (100 1000 °C) and for 5 different status of fast neutron exposure
(0 6.09* 1021 cm 2 , E>0.1 MeV). For T>1000 °C experimental values have been missing.
Therefore, these data have been constructed by extrapolation. The function X (T,D) is included
in the code as default, but it can be replaced by any other function to be defined as data input.
For many other materials the functions X (T) are also included in the code. A survey of these
functions is given in Tab. VII. For more detailed information we refer to subroutine SLAMT.
Similarly, the included functions of the heat capacity (Tab. VI) are given in the subroutine
WKPT.
Lambda (NUKEM / A3-3-Graphite)
0.0
OS -i
0.4 •
°0.2
0.1 -
00 -
Lambda-efT, Pebble-bed
—Retold, min. n-doKZdiner-SdilAnder, nan. n-dis-
— Robdd. m a-doe- - — - Zdiner-Sdiliinder, max. n-do«
1 • i • i i
500 1000 1500
T(°C)
2000 2500 500 1000 1500
T(°C)
2000 2500
Fig. 5: Thermal conductivity
The heat transport through the bed of pebbles takes place partly by thermal conduction
through the pebbles and partly by thermal radiation from one pebble to the other. Theoretical
21
models have been developed for the description of the mechanism of heat transport. At the
end an effective thermal conductivity A.eff (T,D) is defined.
The model of Zehner-Schlünder accounts for the heat transport from one pebble to the next
one. His finding is represented in Ref. I25I. It has been verified in experiments, and it is
recommended for the pebble-bed at low and medium temperature levels.
The model of Robold 1211 takes special care of the heat transport through the openings
between the pebbles of the bed by radiation. Again it derives an A.eff (T,D) which is preferred
at higher temperatures, i.e. for T > 1400 °C.
The effective thermal conductivity of the pebble-bed for both the models is illustrated in
Fig. 5 (right part). It is drawn for the highest and for the lowest neutron dose corresponding to
the function X (T,D) of the graphite. In current calculations the code evaluates both models of
Zehner-Schliinder and of Robold, respectively, and it applies the respective maximum of the
two models.
Tab. VI: Available formulae of heat capacity
ld.no.
1
2
3
4
6
7
8
11
12
13
14
15
16
Material function of temperature dependent
Reactor graphite, SGL, Grade A, NBG10
Reactor graphite, SGL, Grade A, NBG10
Core barrel, SA-240 grade 316
Pressure vessel, SA-508
Steatite
Reactor graphite (HRB)
Carbon bricks (like Reakt. graph.)
V2A - Steel (Hoesch)
Thermal shield (HRB)
Reactor graphite (HRB)
Reactor graphite (HRB)
Reactor graphite (HRB)
AI2O3
heat capacity
density 1.75 gr/ cm
density 1.80 gr / cm3
density 1.75 gr. / cm3
density 1.55 gr. /cm3
DIN 4541
density 1.70 gr. / cm3
density 1.60 gr. / cm3
density 1.80gr. /cm3
no f(T)
22
Tab. VII: Available formulae of thermal conductivity
Id.no.
1
2
3
4
6
7
8
9
10
11
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
Sodium (liquid)
Graphite (Matrix)
Graphite (Reflector)
EPS > 0.: HeliumEPS = 0.: Static Helium
Zehner- Schlünder for Steatite
Reactor graphite
Carbon bricks
Static air
Thermal shield (HRB)
V2A - Steel (Thyssen
Steatite
Graphite balls
Armed concrete
Prismatic core
Carbon-felt
Carbon-felt
Anti-friction bearing steel
Static Nitrogen
Kaowool-mat
H A W-glass
Pebble bed
Ball graphite, Binkele/A3 graphite
LAMDA-eff., pebble bed
Lambda eff., pebble bed
Lambda eff., pebble bed
A12O3
Gilsonit coke (AGL-IE 1-24)
Core barrel, SA-240, grade 316
Pressure vessel, SA 508, grade 3
Reactor graphite, SGL grade A
Material function
T = irradiation temperature
Interpolation from tables (see Subroutine GFIT)
LAM0 = pressure (bar)lbar
Experiment (K. Verfondern)
LAMDA(O) = LAM (Comp.)
"Lukascewicz"
1 bar
DIN 4541
for heterogeneous calculation
Experiment (Robolt)
axial
in vacuum
in Ar- or N2 atmosphere
100CR6
in air (Jül 992RB)
Schürenkrämer (II.84),(Combination of No. 25 and No. 26 for 4xlO21 EDN)
Function of temperature and neutron dose (explicit)
Function of temperature and neutron dose (Robolt)
Function of temperature and dose (Zehner- Schlünder)
Function of temp, and dose (Robolt+Zehner-Schlünder)
linear (Salmang/Scholz "Keramik")
Irradiated at 760 °C (Binkele)
23
3.4 List of data sets
In the following all those units for intermediate or long-term storage of data are listed, whichare not automatically scrapped at the end of a VSOP calculation:
a) VSOP-MS
Log.UnitNo.
5
6
14
Data set name
not fixed
not fixed
'retold'
ASCII (A)/Binary (B)
(A)
(A)
(B)
15
19
2830
'rstnew'
'thermix'
'nucdens''resint'
(B)
161718
'gam_99_4''thermal''thermxxxx'
(A)(A)(A)
(B)
(B)(B)
3537
4258
59
61
6364
'rstcit''geom'
'macsig''therlist'
'tempstat'
'nakure'
'adage''tinte'
(B)(A)
(A)(A)
(A)
(B)
(A)(B)
Description
Card image data input. DS-name to be defined bythe user at the start of program execution.Print output. DS-name to be defined by the user atthe start of program execution.Read restart data (status of reactor operation),which were produced by a previous VSOP-calculation as data set 'rstnew' (see unit 15).Restart data written at the end of the lastcalculated burnup cycle, to be read from unit 14 asdata set 'rstold' in a following restart.GAM-I-library data.THERMALIZATION-library data.with 'xxxx' according to table I:THERMOS-library data.Results of THERMIX are stored onto this data setin a steady-state calculation. They are retrieved fromthis data set in a subsequent time-dependentcalculation.Optional storage of atom densities of each batch.Effective resonance integrals (GAM-group structure)produced and stored in case of a VSOP-ZUT-calcu-lation. To be read and used in a VSOP-MS-calcu-lation.Restart data for neutron diffusion section.Geometry data for neutron diffusion andthermodynamics section.Macroscopic cross sections are stored on this DS.Some results of the transient THERMIX calculationare stored for further use e.g. in external plotroutines.Point temperatures (solid material) of steady-stateTHERMEX-calculation.Power histogram of the fuel batches for decay powerevaluation in subroutine NACHW and in externalcodes NAKURE and TINTE.Contains the ADAGE-library.Preserves data for external use (TINTE code).
24
67686976778098
'ongen'phiform''origmod''birgexcl''trigplot''powform''tempinst'
99 'keff
(A) Preserves data for external use (ORIGEN-JUEL-Ü).(A) Point neutron fluxes.(A) ORIGEN-JÜL-II library (MEDUL-reactors).(A) Core geometry, 2-d.(A) Core geometry, 3-d.(A) Point power densities.(A) Point temperatures (solid material) during
transient THERMDC calculation.(A) Important data of operation-history.
b) VSOP-ZUT
242526
2728
'resdapuO''resdapu2''resdatu5'
'resdatth''resdatu8'
5, 6, 30: see a)
(A)(A)(A)
(A)(A)
VSOP-MS
240,Resonance parameters of 2W PuResonance parameters of 242 PuResonance parameters of 235 U
(not part of the present code package)Resonance parameters of 232 ThResonance parameters of 238 U
25
4. Input Manual V.S.O.P.-MS (log. unit 5)
4.1 Steering the execution mode. SI - S3
Card SI
1
2
3
4
MODE
JSER
I3D
OTT
Format (A4.3I4)
= vsop: Complete V.S.O.P.-calculation. Data input may consist of thefollowing parts: card types S, BI or TR (see variable I3D), D,V, G, T, C, K, R, LF, P, TTTT and TX.
= geom: Provide geometric reactor design only. Data input consists ofcards S, BI (2-d -) or TR (3-d - calculation), exclusively.
= fuel: Provide fuel elements design only. Data input consists ofcards S and D, exclusively.
The values of the following 3 items have no meaning in case of a coderestart, i.e. if JTPE7 > 0 on card S3!
= 0: Diffusion calculation, control poison adjustment, burnup.= 1: Diffusion calculation, burnup.= 2: Diffusion calculation.= 3: Same as 0, but control poison adjustment also in the reflector.= 4: Cell burnup calculation (no diffusion calculation, Kinf and group
fluxes by means of subroutine KINF).= 5: Same as 4, with control poison adjustment.
= 0: 2-dimensional (r - z) - geometry. Cards BI required.= 1: 3-dimensional (x - y- z) - geometry. Cards TR required.= 2: 3-dimensional (O - r - z) - geometry. Cards TR required.
= 0: No effect.> 0: Calculate TTTT (section 4.4.11) in order to prepare a new
30-groups THERMOS-library out of the 96-groupsTHERMALIZATION-library.
26
Card S2 Format (A72)
Literal description of case.
Card S3 only if MODE = 'vsop' on card SI.
Card S3
1
2
3
4
JTPE7
JTPE9
IRR9
IPKEFF
Format (414)
= 0: Normal start.> 0: Restart. JTPE7 is the identification number of restart data to be
retrieved from data set 'rstold'. Data input continues withcards R6 (optionally, see IRR9 below) or R7.
= 0: No effect.> 0: Prepare restart data with id. no. JTPE9 from this calculation,
write them onto data set 'rstnew' to be used for a followingrestart.
= 0: No effect.> 0: For restart only: Card R6 will be given to re-define the options
for the first cycle of the restart.
= 0: No effect.> 0: List of operation-history is displayed at the end of a calculation
and written onto formatted data set 'keff.
27
4.2 Geometrie reactor design
(only if MODE = 'vsop' or 'geom' on card SI)
4.2.1 2-dimensional (r-z - geometry). BI1 - BI9
Cards BI only if I3D = 0 on card SI.
Card BIl
1
2
3
4
NCASE
IPUT
IPLOT
KANAL
Format (416)
= 0: Parallel flow of spherical fuel elements in vertical flow channels,or no movement of the fuel during reactor operation. Read cardsBU - BI5.
= 1: Flow of fuel pebbles with different speed in various radialchannels and/or along non-vertical trajectories. Read cardsBI1 -BI9.
= 0: Normal output.= 1: Test output in addition.
= 0: No effect.> 0: Store data for plots on formatted data set 'birgexcl'.
Pebble bed: Number of flow channels inside the core. (< 15)(see NCASE)
Non-moving fuel: Number of radial coarse meshes (equal IMAXresulting from card BI3).
For each of the KANAL channels one card BI2.
Card BI2
1
2
3
KANTYP
KAN
IBATCH
Format (316)
= 1: For the outermost core channel.= 0: For all other core channels.
> 0: Number of axial regions in this channel (< 100). For eachregion a set of macroscopic cross sections will be generated. Foruse in the diffusion calculation, these data are transferred to thegrid which is defined on cards BI4. The total number of regions(core + reflector area) is restricted to 1500.
= 0: One region only.
> 0: Number of batches per region (< 15 !).= 0: One batch only.
28
Coarse meshes define the CITATION- and THERMIX- material compositions, the finemeshes define the grid for the neutron flux- and temperature calculation.
One card BI3 for each radial coarse mesh I.
CardBI3 Format (I6,E 12.5,16)
1
2
3
IOP(I),
DR(I),
MR(I),1-1,...
- 0: Coarse mesh is situated within the core area.- 1: Coarse mesh is situated outside the core, still covered by the
neutron diffusion section.= 2: Outside the neutron diffusion area, used for THERMIX only.= -1: End of input for radial coarse meshes.
Thickness of the I-th radial coarse mesh (cm). Within the core thewidth of the respective channel is a good choice in most cases.
Number of fine radial meshes in this coarse mesh. Must be =1.if the coarse mesh represents a void area in anv axial position.
The total number of coarse meshes must be < 100. Number of mesheswithin the core area must be < 50
The number of given radial coarse meshes inside the diffusion areadefines IMAX (see card BI5).
One card BI4 for each axial coarse mesh N.
Card BI4
1
2
3
NOP(N)
DZ(N)
MZ(N),N- l
Format (I6,E12.5,I6)
= 0: Coarse mesh is situated within the core area.- 1: Coarse mesh is situated outside the core, still covered by the
neutron diffusion section.= 2: Outside the neutron diffusion area, used for THERMIX only.= -1: End of input for axial coarse meshes.
Thickness of the N-th axial coarse mesh, (cm)
> 0: Number of fine axial meshes in this coarse mesh. Must be =1,if the coarse mesh represents a void area in anv radial position.
= -1: This coarse mesh represents the void above a pebble-bed core.
Total number of coarse meshes must be < 200 (< 100 inside the core).
The number of given axial coarse meshes inside the diffusion areadefines NMAX (see card BI5).
29
Each of the NMAX axial coarse meshes N requires one card BI5.
Card BI5 Format (2413)
IMAX
LAYVCd,N),
1=1, IM AX
> 0: Id.no. of CITATION material composition, which coarse mesh Iis to be assigned to, starting with no. " 1 " for the first composi-tion outside the core. (The numbers are preliminary and will berenumbered successively after the id. numbers of the core areahave been internally defined).
= 0: Core area. Id. numbers are defined by the code.
Cards BI6 - BI9 only if NCASE =1 on card BI1.
Card BI6
1
2
3
4
5
6
KONUS
IZFEIN
JRFEIN
EPSY
RKONUS
ZKONUS
Format (3I6,3E12.5)
= 0: No effect.= 1: An outer cone at the bottom of the core structural material is
present.= 2: A central reflector column with another cone towards the core
bottom is present.
Number of axial meshes of a superposed fine grid for cross sectionand flux transfer matrix (< 25000, see also section A.2.1).
Number of radial meshes of the superposed fine grid (< 10000).
Convergence criterion for the iteration on the radial position of themesh points defining the flow channel curves (about 1 .E-5).
Only if KONUS > 0:
Radial thickness of the cone(s). (cm)
Height of the cone(s). (cm)
30
For each of the KANAL channels one set of cards BI7 - BI9.
CardBI7 Format (112,E 12.5)
1
2
UR
VEKA
> 0: Number of axial mesh points for the definition of the outerlimiting curve of this channel (< 15). Only for the inner corechannels (KANTYP = 0). UR = 1 defines a straight vertical line.
= -1: The value of UR is taken from the preceding channel. Dropcard BI8.
= 0: Last core channel. Drop cards BI8, BI9. Limiting curve isinternally defined by the information of card BI3.
> 0.: Only when KANTYP = 0. Ratio of the volume of this channelper volume of the core. Radial mesh points of the limiting curvewill be adapted to meet this volume of the channel.
= 0.: No adaptation of the limiting curve.
Card BI8 only if IJR > 0 on card BI7.
CardBI8 Format (6E12.5)
1
UR
XWE(J),J=1,UR
Axial position of the coarse mesh points for the outer limiting curveof this channel (cm), starting from the top of the core (XWE = 0.)down to the bottom (positive values).
Card BI9 only if UR * 0 on card BI7.
Card BI9 Format (6E12.5)
1
UR
YWE(J),J=1,UR
Radial position of the coarse mesh points for the outer limiting curveof this channel (cm). In case of an annular core, YWE(J) must begiven as the distance from the inner limiting curve of the first corechannel.
31
4.2.2 3-dimensional. TR1 - TR5(only if MODE = 'vsop' or 'geom' on card SI)
This part of code defines the 3-dim. pattern of regions in the reactor. VSOP-regionsand CITATION-material compositions are identical. They are assigned with thesame id. numbers. (See also chapter A.2.1).
Cards TR only if I3D > 0 on card S1.
CardTRl
1
2
3
MYX
IPL
ICORE
Format (315)
= 0: 1 batch per region.> 0: Number of batches per region.
Note: The total number of regions as well as of batches is limited to9999 !!
= 0: No effect.> 0: Plot data for plane no. IPL is written onto data set 'trigplot'.
(For special purposes only).
= 0: Normal.> 0: Just for plot data and only for x-y-z-geometry:= 2: Data of Vi core-plane are transmuted to 1/1 core-plane.= 4: Data of lA core-plane are transmuted to 1/1 core-plane.
Meshes in X-direction (I3D = 1) or in <J> - direction (I3D=2).
Card TR2
1
2
MX(I>.
DX(I),1=1, ...
Format (6(I3,F9.3))
Number of fine meshes in the I-th coarse mesh in X- or in O-direction.
> 0.: Thickness of the I-th coarse X-mesh (cm) or O-mesh (degrees)= 0.: End of the input of coarse X- or <I> - meshes.
The number of given coarse X / <£-meshes defines IMX.
32
Meshes in Y-direction (I3D = 1) or in R - direction (I3D = 2).
Card TR3
1
2
MY(J),
DY(J),J=l,.. .
Format (6(I3,F9.3))
Number of fine meshes in the J-th coarse mesh in Y / R-direction.
> 0.: Thickness of the J-th coarse Y / R-mesh. (cm)= 0.: End of the input of coarse Y / R-meshes.
The number of given coarse Y / R-meshes defines JMY.
Meshes in Z-direction.
CardTR4
1
2
MZ(K),
DZ(K),K=l,...
Format (6(I3,F9.3))
Number of fine meshes in the K-th coarse mesh in Z-direction.
* 0.: |DZ(K)| gives the thickness of the K-th coarse Z-mesh. (cm)< 0.: Core regions.> 0.: Non-core regions (e.g. reflectors).
= 0.: End of the input of coarse Z-meshes.
The number of given coarse K-meshes defines KMZ.
Definition of the pattern of regions:
For each of the planes (Z) K = 1 ,KMZ one set of cards TR5.For each of the rows (R or Y) J = 1JMY one card TR5.
Card TR5 Format (1515)
1
IMX
LAY3(I,J,K)1=1,IMX
Only for the core:> 0: Region id. number of the I-th coarse mesh (O or X) in this row
and plane. In the upper plane the code evaluates the maximumnumber NLP of core compositions.
= 0: Id. no. of this region is internally defined by adding NLPto the LAY3(I,J,K-1) of the foregoing plane.
33
Continuation of card TR5
Only for the non-core-compositions (reflectors etc.):< 0: Id. no. of this composition is internally defined by adding the
maximum number of core compositions to the absolute of|LAY3(I,J,K)|. Reflector id. numbers must be given in an un-broken sequence starting with " -1 " .
34
4.3 Fuel Element Design. Dl - D17(only if MODE = 'vsop' or 'fuel' on card SI)
4.3.1 Specifications. Dl - D4
CardDl
1
2
3
KMAT
NHOM
NDANC
Format (314)
Sum of nuclides - except the heavy metal isotopes (see Table IV) - tobe included in the VSOP calculation (i.e. fission products, controlpoisons, non burning absorbers, scatterers). (< 172)
= 0: Normal.= 1: Drop heterogeneous evaluation of fuel elements. Read homo-
genized atom densities on cards D17.The input sequence is : Cards Dl, D2, D5, D6, D17, D5.
Only when NHOM = 0:= 0: No effect.> 0: Definition of lattice type according to card Z10 (variable K) of
VSOP-ZUT- input (section 5.3.3). Some values are calculated andprinted, which should be used as input data for the Dancoff factorcalculation in VSOP-ZUT.
CardD2 Format (1814)
KMAT
IMAT(I),1=1,KM AT
GAM-I-library identification numbers for the KMAT (card Dl) nu-clides in the sequence of Table IV, starting with the fission products(drop the heavy metal isotopes!).
Note:Nuclide id. numbers beyond the library can be used (i.e. IMAT(I)> 190). These nuclides must be identified on the cards V3.
35
Cards D3, D4 only if NHOM = 0 on card Dl.
Card D3
1
2
3
4
9
10
12
13
CURCY
NC
NF
FC(I),1=1,6
FF(I),1=1,3
DK
Format (A4,2I4,10F6.0)
Literal abbreviation for the monetary unit MU.
Number of different coated particle fabrication cost data. (< 6)
Number of different fuel element fabrication cost data. (< 3)
Fabrication costs of the I-th coated particle variant. (MU/ fuel ele-ment)
Fabrication costs of the I-th fuel element variant. (MU/ fuel element)
Fabrication costs of "dummy" elements. (MU/ element)
Card D4 Format (3E 12.5)
1
2
3
HK
AK
EK
Costs of head end and transportation. (MU/kg c )
Costs of reprocessing. (MU/kg HM )
Costs of waste treatment and disposal per 10% Fima. (MU/kg HM )
4.3.2 Design of fuel element-types and -variants. D5 - D17One set for each variant of each desired fuel type (limited to 27 different sets).Calculation is terminated by one last card D5.
Card D5 Format (18A4)
1
18
TITLE(I),1=1,18
Literal description of fuel element-types and -variants.
TTTLE(l) = 'stop1: This terminates the sequence of cards D5 - D17.
36
Card D6
1
2
3
4
5
6
NTYP
NFUTP
NFCP
NFBZ
NZUS
FF3
Format (5I4.E12.5)
= 0: Spherical fuel elements.= 1: Prismatic block fuel element.
Fuel elements: NFUTP is given as a positive number,Fuel-free ("dummv"-) elements: NFUTP is given as a negative
number:Identification of the elements in 4 digits IJKL:LI: Type (of design and cost data), increasing numbers, (< 10)KL: Variant ( e.g. for different enrichments), increasing numbers.
starting from 01 for each type LI.
Input option for coated particle definition:= 0: Data from preceding design, or if NFUTP < 0 (dummy elements),
orifNHOM = l (card Dl).= 1: Read card D7 only.= 2: Read cards D7-D11.
Option for pebble and block type element specific data:= 1: Read cards D12 - D13 or D14 - D16, respectively.= 0: Respective data like preceding design, or if NHOM = 1 .
= 0: No effect.> 0: Number of nuclides for which atom densities will be specified on
cards D17. (<30)Previously calculated values are replaced.
Variable is used for spherical elements only (NTYP=0).> 0: Volumetric filling fraction of spherical elements in the core.= 0: Default value (0.61) is used.
4.3.2.1 Coated particles. D7-D11
Card D7 only if NFCP > 0 on card D6.
Card D7
1
2
ANR
U53
Format (6E 12.5)
Fissile enrichment of the fuel (fissile/heavy metal).> 0.: Atom fraction.< 0.: Weight fraction = |ANR|.= 0.: If INDBS (card D8) = 7 .
= 0.: Fissile uranium is 235U.= 1.: Fissile uranium is 233U.
37
Continuation of card D7
FIMA
FRC(I),1=1,NC
Envisaged heavy metal burnup for reprocessing cost calculation.(FIMA)
> 0.: Fraction of coated particle variant I in this fuel element.= 0.: Coated particle variant I (card D3) is not present.
Cards D8 - Dl 1 only if NFCP = 2 on card D6.
Card D8
1
2
3
4
INDBS
NCT
NSIC1
NSIC2
Format (414)
Fuel identification:= 1:UO2 =2:UC= 3: UC2 = 4: UO2 - ThO2
= 5: UC - ThC = 6: UC2 - ThC2
= 7:PuO2 =8:PuO2-ThO2
= 9: PuO2 - UO2
Total number of coating layers (< 5), to be numbered with increasingradius.
Number of the 1. SiC coating layer, if present.
Number of the 2. SiC coating layer, if present.
Card D9
1
2
3
4
RK
ROBR1
ROBR2
BETA
Format (4E12.5)
Radius of the coated particle kernels, (cm)
Density of the kernels, (g/cm3)
Density of 2. type of kernels, if present, (g/cm3)Only if INDBS = 4, 5, 6, 8, 9 on card D8.
Enrichment of the uranium NU5 / Nu if INDBS = 4, 5, 6, 9 on cardD8.For U53 = 1. (card D7) program uses 233U instead of 235U.
38
Card D10 only if INDBS = 7, 8 or 9 on card D8.
CardDIO Format (4E 12.5)
1
41=1,4
Atom fractions of the isotopic composition in plutonium:239Pu,24OPu,241Pu,242Pu.If required, an additional 238Pu concentration may be defined on cardD17.
CardDll Format (6E12.5)
1,3,5
2,4,6
DCT(I),
ROCT(I),1=1, NCT
Thickness of the I-th coating layer, (cm)
Density of the I-th coating layer. (g/cm3)(Numbered with increasing radius, NCT on card D8).
4.3.2.2 Spherical fuel elements. 1)12,1)13
Cards D12 - D13 only if NTYP = 0 and NFBZ = 1 on card D6.
Card D12
1
2
3
4
5
Rl
R2
FF1
VMOD
ROSM
Format (5E12.5)
= 0: for "Dummy"-Elements> 0: Outer radius of fuel zone, (cm)
(Fuel zone consists of coated particles and graphite matrix).
Outer radius of the sphere, (cm)
Only one of the following three variables must be specified in case offuel elements. In case of "dummy" elements, they all must have a zerovalue.
Volume fraction: coat.part. / (coat.part. + matrix).
Moderation ratio Nc / NHM-
Density of the heavy metal, homogenized in the fuel zone (g / cm )
39
Card Dl3
1
2
3
4
ROMTX
ROSCH
SRO
FRF(I),1=1, NF
Format (6E12.5)
Density of graphite in the matrix, (g/cm3). (= 0 for "dummy"elements)
Density of graphite in the outer shell, (g/cm3)
Inner radius of the matrix, (cm) (normally = 0.)> 0 for "shell ball" design.
= 0.: Fabrication costs of the I-th fuel element variant (card D3) aredropped (always true for "dummy" elements).
> 0.: Fabrication costs of the I-th fuel element variant are used andmultiplied by FRF(I). (Usually =1.)
40
4.3.2.3 Prismatic fuel elements. D14-D16
Cards D14 - D16 only if NTYP = 1 and NFBZ = 1 on card D6.
CardDH
1
2
3
4
5
6
R(l)
R(2)
R(3)
R(4)
R(5)
R(6)
Format (6E12.5)
Radius of central graphite zone, (cm)
Outer radius of inner cooling channel, (cm)
Outer radius of inner graphite tube, (cm)
Outer radius of the fuel zone, (cm)
Outer radius of the outer graphite tube, (cm)
Outer radius of the outer cooling channel, (cm)
If FFUEL > 0. (card D15) insert "thickness" instead of "radius".
Only a selected set of the following parameters of the cards D15 and D16 is required. Possiblecombinations are given in Table VIII.
Card Dl5
1
2
3
4
5
6
FF1
VMOD
BETA
GKAN
FFUEL
ACTIV
Format (6E12.5)
Volume fraction: coat.part. / (coat.part. + matrix).
Moderation ratio N c / NHM-
Volume fraction of gaps (other than cooling channels) in the corerelative to the bulk graphite volume.
Number of fuel elements per square meter. (1/m2)
= 0.: No effect.> 0.: Cross section of the fuel zone (cm2). Use R(4) = 0. on card D14
and insert "thickness" instead of "radius".
Active length of the fuel rods in the core, (cm)
41
Card Dl6 Format (4E 12.5)
1
2
3
4
ROSM
ROMTX
ROSTR
ROHR
Density of heavy metal, homogenized in fuel zone, (g/cm3)
Density of graphite in the matrix, (g/cm3)
Density of graphite in the cooling channel, (g/cm3)
Density of graphite in the tubes, (g/cm3)
Table VIII: Alternative specifications of fuel rods
No.
FF1
VMOD
GKAN
ROSM
1 2
X X
X
X
3
X
X
4
X
X
5
negative guess
X
X
4.3.2.4 Additional nuclides. D17
Card(s) D17 only if NZUS > 0 on card D6.
CardD17
1
2
NRGAM
DENG
Format (14,4X,E 12.5)
GAM-I-lib. identification no. of nuclide with additionally given atomdensity.
Atom density (atoms / (barn cm), homogenized).
Note: Use 1 card for each of the NZUS (< 30) additional nuclides.
42
4.4 Reactor and fuel cycle. VI - TX26(only if MODE = 'vsop' on card SI)
4.4.1 Set up dimensions. VI
Card VI
1
2
3
4
5
6
7
8
N26
MMAF
MBATCH
MSTOB
JTYP
MREP
JABOX
KMAZ
Format (814)
Number of energy groups in the diffusion calculation. (< 33)
Maximum number of bumup cycles.
Maximum number of batches to be filled into storage boxes.(See card R21).
Maximum number of storage boxes to be filled.(See card R21).
Number of different fuel element types in the system. (< 10)(See card R3).
Number of reprocessing mixtures, if present. (^10)(See card R3).
Total number of aging boxes and jumble boxes as explicitly specifiedon card R5 (only if MREP > 0).
See also chapter A.2.2 .
Only if "I3D = (Toncard SI:> 0: Maximum number of THERMIX- (= KONVEK-) compositions.= 0: Default value = 50 .
43
4.4.2 Definition of materials. V2-V5
Card V2
1
2
3
4
NO
KETT
NLT
NC
Format (414)
Number of fission products (< 49):= 0: Default value = 44, code uses the built-in fission product chain
(see Fig. 4).0<NO<44: The code drops the last surplus ones of the built-in chain.> 44: See KETT and card V4.(See also cards Dl (KMAT) and D2).
= 0: No effect.> 0: Chain information of the last KETT fission products will be
defined on cards V4. This option can be used to extend the chainstructure or to define a new one.
= 0: No effect.> 0: Number of fission products, for which new yields and decay
constants will be defined on card V5.
Number of control poison nuclides. (< 2)
Card(s) V3 only if some nuclides of the library shall be duplicated and used with new id. numbersIMAT(I) > 190 for special purposes. One card V3 for every new id. number.
Card V3 Format (214)
JNEU
LMAT
GAM-I-Id. number to be assigned to the new nuclide.
GAM-I-Id. number of the original library nuclide of which the crosssections are to be duplicated.
44
Card(s) V4 only if KETT > 0 on card V2.A total of KETT cards required, starting with the card for the fission product nuclideN = NO-KETT+1 .
CardV4 Format (4E12.5)
2
3
4
DIRAC(N,1)
DIRAC(N,2)
DIRAC(N,3)
DIRAC(N,4)
Fractional production of nuclide N from N-l.> 0.: By capture.< 0.: By decay.
Fractional production of nuclide N from N-2.
Fractional production of nuclide N from N-3.
Fractional production of nuclide N from N-4.
Card(s) V5 only if NLT > 0 on card V2.A total of NLT cards is required, one for each fission product for which the yields are defined oraltered.
CardV5
1
2
3
4
5
6
N
YIELD 1(N)
YIELD2(N)
Y1ELD3(N)
YIELD4(N)
XLAM(N)
Format (I6,6X,5E12.5)
VSOP identification no. of a selected fission product nuclide.
233U fission yield of nuclide N.
235U fission yield of nuclide N.
239Pu fission yield of nuclide N.
241Pu fission yield of nuclide N.
Decay constant of nuclide N. (1/sec)
45
4.4.3 Design and operations. V6 - V17
4.4.3.1 Case identification. V6
Card V6
1
2
3
4
5
6
7
8
9
NRSTRT
NKOST
IBUCK
MUHU(3)
LOBNEW
IBASCH
IPRIN2
SERCON
ERR
Format (I8,6I4,2E 12.5)
= 0: No effect.= 1: Fuel shuffling.= 2: Fuel shuffling and iteration of the enrichment.= 3: Fuel shuffling and reprocessing.= 4: Fuel shuffling, reprocessing and iteration.
* (Cards R28-R31).
= 0: No effect.> 0: Fuel cycle cost calculations (cards Kl - K12).
= 0: No feedback of leakage from diffusion to spectrum calculation.= 1: Feedback of the broad group leakage to GAM-I, and thermal
leakage to THERMOS.= 2: Feedback of an average epithermal leakage to GAM-I, and
thermal leakage to THERMOS.
= 0: Drop streaming correction in pebble bed.= 2: Streaming correction LIEBEROTH /28/ in power generating
batches (only for pebble bed).
= 0: Normal.= 2: Life history is preserved for ORIGEN-JÜL-II (all NON-MEDUL-
reactors, see 151), starting from the first burnup cycle (data set'origen').
= 0: No effect.> 0: For 3 - d - geometry: Number of batches in the upper plane of
the core. Will be used only in connection with variable MULTon card V7.
> 0: Print layout of batches at startup.= -1: No output.
Convergence criterion for Keff when adjusting control poison or otheratom concentrations. (= 0.0001)
> 0.: Truncation error limit to be used for the burnup and spectrumcalculation.
= 0.: Default value = l.E-25 .
46
4.4.3.2 Definition of reactor batches. V 7 - V 9
Each batch (see output of "Geometric reactor design" - section) requires one set of cards V7-V9.The sequence has to be as follows:1) In-core batches, numbered from 1 through ....2) Cone regions (= batches), numbered from (in-core batches + 1) through ....3) Other non-power generating regions, numbered from 1. through ....
(The number of in-core + cone batches is automatically added up).
CardV7
1
2
3
4
5
NREAD
NCH1
NCH5
NFTST
WPART
Format (4I6,6X,E12.5,18X,2I6)
< 100: Number of atom densities to be specified on subsequent cardsV8.
> 100: Atom densities of fuel type IJ, variant KL are used.(NFUTP on card D6).
NCH1=O, if NREAD * 0 .NCH1*O, if NREAD = 0:
> 0: Number of a previously specified batch with the sameatom densities.
< 0: Read new atom densities for this batch from dataset 'nucdens', which have been stored in an earlier VSOPcalculation (compare variable LIB < 0 on card R7).|NCH1| is the batch no. of data set 'nucdens' to be applied.
> 0: Number of a previously specified batch with the same controlpoison data.
= 0: Use data of batch no. 1.< 0: Read card V9.
Only if NREAD < 100:Definition of fuel type id. no. of this batch. Only if fuel is defined bycards V8, otherwise the id. no. is taken from batch no. NCH1.In reflector regions the id. no. is 0 .
Fraction of the volume of this batch per region:= 0.: In the batches of the first region the code makes WPART =
1. / (number of batches per region). In the other batches of thecore the code copies WPART of the corresponding batch ofthe preceding region. In the regions of the reflectors the codemakes WPART = 1.
> 0.: Redefinition of the volume fraction of this batch. If redefinitionis specified, it must be given for all batches of this region, and itholds for all subsequent regions until redefined.
47
Continuation of card V7
6
7
MULT
KD18
= 0: No effect.> 0: The id. no. of this batch is defined by KD18 = KD18 + MULT *
IBASCH (Card V6). (Useful in 3-d- geometry for the batchesin the lower planes).
> 0: Id. no. of the batch for which the information of this card is to beapplied. It will also be applied for all subsequent batches untilredefined.Note:The sequence of the reflector batches (= regions) must correspondto the numbering defined on input cards BI5 or TR5, respectively.
< 0: Last card V7, holding for the batch |KD18|.
Atom densities: Card V8 only if 0 < NREAD < 100. A total of NREAD cards is required.
Card V8 Format (I4.4X.E12.5)
1
2
L
DEN
VSOP-identification number of the nuclide with atom density > 0.
Atom density (atoms per barn cm). All densities must be givenhomogenized.
Control poison: Card V9 only if NC > 0 on card V2, and if NCH5 < 0 on card V7. One card V9for each control poison nuclide.
CardV9 Format (2E 12.5)
1
2
POISM
POISL
The control poison nuclide(s) in all batches of one region have thesame limitations.
Minimum atom density of control poison in this region, (e.g. = 0.)
Maximum atom density of control poison in this region.
48
4.43.3 Data for the burnup calculation. V10, V l l
Cards V10-V11 only if JSER * 2 on card SI.
Card V10
1
2
3
4
DELDAY
POWER
FIWATT
ZKFIND
Format (4E12.5)
Length of large burnup time steps (time between possible diffusioncalculations), (days)
Thermal core power, (watts)
Initial value of Fissions/Ws:= 0.: Starting value = 3.087E+10 (235U)> 0.: Optional starting value.Note:In the course of the proceeding burnup, an actual value of FIWATT iscalculated by the code according to DIN 25485 /26/. This value thendepends on the fraction of the fission rates of the different fissileisotopes.
Minimum allowed value of Keff. The present burnup cycle is termina-ted, when Keff equals ZKFIND. Fuel shuffling is then performed, ifspecified. In case of control poison adjustment, ZKFIND is the targetKeff.
Card Vl l Format (214)
1
2
JNSTOP
JNUM
Last large burnup time step in one burnup cycle. (< 95)
Number of small time steps in one large step. Renormalization of theneutron flux to the specified reactor power is done for each small timestep.
49
4.43.4 Control poison search. V12 - V14
Cards VI2-V14 only if JSER = 0, 3, 5 on card V6.
Card V12
1
2
3
4
5
JSMAX
JSSMAX
LSIM
KSS
NPOIS(I),1=1,KSS
Format (1814)
Maximum number of control poison iterations for any region at onetime step. All batches of the region are treated simultaneously. (=50)
Maximum number of control poison iterations for the total core at onetime step. (= 200)
Number of regions, for which the control poison is adjusted simul-taneously. LSIM regions form a poison area for simultaneous poisonadjustment.
Length of the list of regions for control poison adjustments. The ratioKSS / LSIM gives the number of poison adjustment areas.
This list gives the sequence of regions in which the adjustments areperformed.
Card V13 Format (6E12.5)
KSS
PINMIN(I),I=1,KSS
Minimum fraction of control poison insertion in the I-th region to beadjusted, (e.g. = 0.)
Card V14 Format (6E12.5)
KSS
PINMAX(I),1=1,KSS
Maximum fraction of control poison insertion in the I-th region to beadjusted, (e.g. = 1.)
50
4.43.5 Print-out options and steering. V15
CardV15
1
2
3
4
5
IPRIN(l)
IPRIN(2)
IPRIN(3)
IPRIN(4)
IPRINO
Format (514)
Spectrum calculation:= -1 : Minimal output.= 0: Thermal selfshielding factors, only.= 1: Same as 0, plus averaged thermal cross sections.= 2: Same as 1, plus fine group neutron fluxes.= 3: Same as 2, plus broad groups averaged cross sections for
materials with concentration > 0.= 4: Same as 3, for all materials.= 5: Maximum output including details of neutron transport.
= 0: No output.= 1: Print layout of batches before shuffling.= 2: Same as 1, plus atom densities (only in combination with
IPRIN(3) > 0).
Burnup calculation:= -1: Global neutron balance.= 0: Detailed neutron balance.= 1: Same as 0, plus characteristic data for all fuel batches.
= 0: Perform spectrum calculation only at start of first burnuptime step.Instructions on card V16 are neglected.
= 1: Repeat spectrum calculation as defined on card VI6.
Burnup calculation (ADAGE):= 0: No output.= 1: Short output (cross sections + total flux).= 2: Detailed output.
4.4.3.6 Steering the performance for spectrum and diffusion calculation. V16, V17
Card V16 Format (1814)
ISPEKT(l) > 0: No. of the first large burnup time step in which the spectrumcalculation is to be repeated prior to the diffusion calculation.
51
Continuation of card V16
18
ISPEKT(I),1=2,18
> 0: No. of further time steps for spectrum calculation.= 0: If all ISPEKT = 0, spectrum calculation is performed in every
time step.
Card V17 only if JSER < 4 on card V6.
Card V17 Format (1814)
18
IDIFF(I),1=1,18
If all IDIFF(I) = 0:Diffusion calculation is performed at every time step.
If at least one IDIFF(I) * 0:The IDIFF(I) give the time steps at which diffusion calculation is to beperformed.
4.4.4 Fast and epithermal neutron spectrum. Gl - G12
CardGl Format (18X, 516)
IDESIN Number of different fuel element designs (< 10). Only for differentresonance integral data on cards G3-G5. The differentiation of fuelelement designs for the resonance calculation is mostly the same asfor the thermal cell calculation, i.e. IDESIN = NBER on card T6.
52
Continuation of card Gl
2
3
4
5
MSTU
MGHUS
NSSS
IPRSEL
Fission source spectrum:= 3:233U.= 5:235U.= ll:239Pu.= 13:24IPu.= 0: Unit fission source.
Only if MSTU = 0:GAM-I group no. in which the unit fission source is located.
= 0: No selfshielding factors applied.> 0: Number of sets of selfshielding factors (cards G7-G12).= -1: One single set of selfshielding factors to be applied in all regions
(cards G8-G12).
Output option of the selfshielding factors:= 0: Broad energy group definition.= 1: Selfshielding factors for the different nuclides.
CardG2 Format (6E 12.5)
NDR
TEMZUT(I),1=1,NDR
Temperature of the resonance absorbers in "NDR" different spec-trum calculations. (°C).= 0., if no fuel in the regarded region, e.g. for reflectors.
NDR is the total number of "regions", which is depicted in the outputof code section "Geometric reactor design". (Table: Region -
batchesin the region).
Card(s) G3 only if IDESIN > 1 on card Gl.
CardG3 Format (1216)
1
NDR
NDES(I),1=1,NDR
Fuel element design number used for the spectrum calculation inregion I.
53
For each design (IDESIN on card Gl) 4 sets of cards G4-G5, the first set for 232Th, the second setfor 238U, the third set for 24OPu, the forth set for 242Pu.
Card G4
1
2
3
SMI
SM2
NZ
2 cards G5, the first one
CardG5
1
NZ
IZUT(K),K=1,NZ
Format (2E 12.5,16)
Two values of homogenized atom densities of the resonance absorbernuclide, for which sets of resonance integrals are available on data set'resint'. These values should represent the highest and the lowestdensities, occurring within the reactor, respectively.
> 0.: Highest density of the absorber nuclide. [barn' cm']= -1.: Density is taken from data set 'resint'.
> 0.: Lowest density of the absorber nuclide. [barn1 cm"']= -1.: Density is taken from data set 'resint'.
Number of sets of resonance integrals for each SMI and SM2 den-sities (< 10 !). These sets represent different temperatures of thisabsorber nuclide.
for SMI, the second one for SM2.
Format (1216)
Id. numbers of the resonance integral sets to be read from data set'resint'.
Definition of broad energy groups.
Card G6 Format (6E 12.5)
CEG(I),I=1,N26-1
Desired lower energy limit of the fast energy group(s). (eV)
54
Individual epithermal selfshielding factors. G7 - G12
Card G7 only if NSSS > 0 on card Gl (see also section A. 1.3).
CardG7 Format (1216)
NDR
NSET(I),1=1, NDR
Id. no. of the set of selfshielding factors to be applied in the spectrumcalculation for region I.
Cards G8 - Gl2 only when NSSS * 0 on card G1:For each of the |NSSS| sets of selfshielding factors one set of cards G8-G12.
Card GS
1
2
3
Card(s)
Card G9
1
MOBG
MOBG
LSUB
NK
Format (316)
> 0: Number of broad epithermal energy groups for input of self-shielding factors (card G9). (Up to 67 groups can be defined).
= 0: Broad energy groups same as defined on card G6 (the code setsMOBG = N26 - 1).
< 0: Same broad energy groups as defined before.
> 0: Number of sets of cross section-selfshielding factors SC(cards G10). (<9)
= 0: No input of SC.
> 0: Number of sets of neutron flux-selfshielding factors SF(cards Gil) . (< 6)
= 0: No input of SF.
G9 only if MOBG > 0 and MOBG < 67.
MGBN(J),J=1,MOBG
Format (1216)
Id. number of the GAM-I group with the highest energy in the broaderenergy group J.
55
Card(s)G10only if LSUB > 0 on card G8.A set of (J=1,MOBG) cards G10 must be given for the MOBG broad energy groups.
CardGlO Format (6E 12.5)
1
LSUB
SC(L,J),L=1,LSUB
Broad energy group J:Cross section-selfshielding factor of set L.
Card(s) G11 only if NK > 0 on card G8.A set of (J=1,MOBG) cards Gl 1 must be given for the MOBG broad energy groups.
Card Gil Format (6E 12.5)
1
NK
SF(KJ),K=1,NK
Broad energy group J:Neutron flux-selfshielding factor of set K.
A card Gl 2 is required for each nuclide (simplification of input can be defined by variable JT).
CardG12 Format (I6,2I2,6E 10.4)
1
2
3
4
IDG
JT
LSC
ANT(K),K=1,NK
> 0: Id. no. of nuclide in the GAM-I library in rising sequence.Nuclides standing before the first given id. no. are assigned withselfshielding factors equal 1.0 .
< 0: This is the last card G12.
= 0: Information of this card applies also for all following nuclides,unless revised.
= 1: Information of this card applies only for this nuclide.
= 0: No cross section-selfshielding factors applied.> 0: Id. no. of cross section-selfshielding factors (SC(LSCJ), J = 1,
MOBG) to be applied for this nuclide.
Fraction of the homogenized atom density to be assigned to neutronflux-selfshielding factor set K (only if NK > 0).
56
4.4.5 Thermal cell spectrum - Tl - T13
CardTl
1
2
3
4
5
6
7
CIDTHER
NKER
NKERAB
NUTTE
NUCT
ITY
MUP
Format (A9,3X,6I6)
Data set name of THERMOS-library to be used (one of the existingaccording to chapter 3.2, Table I, or a new one, which must havebeen generated before by means of n i l ) .Blank value: Calculate TTTT with THERMALIZATION (ITTT > 0
on card SI). It shall be made for one spectrum zone, i.e.for one region, only.
The following variables must have a zero value, if no data set name isassigned to variable CIDTHER.
Number of scattering nuclides (< 5). (See also Tables III and IV)
Number of absorber materials for which a scattering matrix is calcu-lated internally (Brown St. Johnes). (< 10)
= 1: THERMOS calculation, no subsequent calculation of self-shielding factors.
= 2: Calculation and print-out of self shielding factors.
Maximum number of scattering matrices - according to differenttemperatures - per one scattering nuclide to be used for interpolation(cards T3). (< 20)
Identification of cell definition of which the geometry data are usedfor
streaming correction.= 0: Use the first cell definition (for which the first set of cards
T7-T12 is given).> 0: Use the ITY-th cell definition.= -1: Define the geometry data to be used for the streaming
correction on card T13. This is necessary if THERMOS celldefinition is different from the real fuel element size (e.g. in caseof "dummy" elements admixed to the fuel elements).
= 0: Normal.> 0: Number of broader thermal groups (to be read on card(s) T2) for
given individual selfshielding. (< 30)
57
Card(s) T2 only if MUP > 0 on card Tl.
CardT2 Format (6E 12.5)
MUP
EMU(I),1=1,MUP
Upper limit (eV) of the I-th thermal broader group for the givenindividual selfshielding, starting with the lowest thermal group.
For each of the NKER scattering nuclides one set of cards T3-T4, in sequence of VSOP nuclideson card D2.
CardT3
1
NUCT
Card T4
1
NDR
IKER(J),J=1,NUCT
TCELS(J),J=1,NDR
Format (1216)
Id. no. of the J-th scattering matrix to be used for interpolationaccording to the actual temperature of this scattering nuclide.
Format (6E12.5)
Temperature of this scattering nuclide for spectrum calculation ofregion J. (°C)
Card T5 only if NKER AB > 0 on card T1.
Card T5 Format (El2.5,1016)
TKG
IDTA(I),I=1,NKERAB
Relative temperature in the calculation of scattering matrices forabsorber nuclides. (°K/293.6)
GAM-I-identification no. of nuclide I for which a scattering matrix iscalculated internally.
58
CardT6
1
2
NBER
NTYSP(I),1=1,NDR
Format (1216)
Number of different cell definitions for the thermal spectrum calcu-lation. Mostly same as IDESIN on card Gl.
Identification no. of the cell definition used for spectrum calculationof region I.
For each of the NBER cell definitions one set of cards T7 - T12 required.
CardT7
1
2
3
4
5
6
7
NGEOM
TKG
FUELL
FUTYP
STRTO
PNORM
TLEAK
Format (I6.6F6.0)
= 0: Cylindrical fuel element.= 1: Spherical fuel element.
Temperature for the initial guess of Maxwell neutron spectrum,relative to TO. (See STRTO on this card). = 1.5
Ratio: Volume of the cell / homogenized volume.
> 0.: Cell definition according to fuel element type identification U asspecified by NFUTP on card D6.
= 0.: All cell specifications must be given on cards T8-T9.Necessary in case of "dummy" elements admixed to the fuelelements.
> 0.: Identification of the most important scattering nuclide. Itstemperature will be used as base temperature TO, as required forTKG.STRTO = 1., 2 identifies the first, second .... scatterer insequence of VSOP nuclide list (cards D2, T3).
= 0.: TO = 293.6 °K.
= 0.: Average cross sections are based on the average cell flux.> 0.: Average cross sections are based on the flux at the mesh point
PNORM.< 0.: Average cross sections are based on the flux at the outer edge of
the cell.
= -1.: Isotropie boundary condition. Read card T12.= 1.: White boundary condition.
59
Card T8
1
20
21
22
23
24
25
26
MTBL(J),J=l,20
IBRENN
ICOAT
COA(l)
COA(2)
COA(3)
COA(4)
Card T9 only if FUTYP
Card T9
1 RED(I+1),1=1, NCZ
Format (2011,2I2.4E12.5)
Cell zone no. in which mesh point J is located.E.g. 11122223330000000000. The highest digit defines the number ofcell zones NCZ (< 9). Each cell zone must contain a scatteringnuclide.
Skip if FUTYP >0.Cell zone no. in which the fuel is located.
Skip if FUTYP >0.> 0: Cell zone no. in which the coated particles are located.= 0: No calculation of coated particles heterogeneity.
Skip if FUTYP > 0. or ICOAT = 0 .Radius of the coated particle kernel, (cm)
Skip if FUTYP > 0. or ICOAT = 0 .Outer radius of the coating, (cm)
Skip if FUTYP > 0. or ICOAT = 0 .Volume fraction: coat. part. / (coat. part. + matrix).
Skip if FUTYP > 0. or ICOAT = 0 .Ratio: Nuclide density of matrix / total nuclide density of coating.
= 0. on card T7.
Format (6E12.5)
Outer radius of cell zone no. I. (cm)Inner radius of cell zone no. 1 is set to 0. (See variable MTBL on cardT8 for NCZ).
A card T10 is required for each nuclide. For simplified input see variable JT.
Card T10
1 IDISO
Format (15,14,11,10F6.3)
> 0: Id. no. of nuclide in the THERMOS-library, starting with ab-sorber nuclides in sequence of increasing THERMOS-librarynumbers. Followed by the scatterers with modified numbers1000+J. Here, J = 1, 2 .... identifies the first, second .... scattererin sequence of VSOP nuclide list, see cards D2 ,T3.
< 0: -IDISO terminates the input of cards T10.
60
Continuation of card T10
2
3
4
5
Card(s)
CardTl
1
MUP
Card(s)
MUPN
JT
VB(1)
VB(L),L=2,NCZ
VB(NCZ+1)
= 0: Normal.> 0: Individual thermal selfshieldings of this isotope are given on
cards T i l .
= 0: The fractional densities VB specified on this card are also validfor all subsequent nuclides, unless revised. This holds also for theselfshielding SFMU (card Tl 1), if defined.
= 1: The VB are valid only for this nuclide.
> 0.: Fraction of the homogenized atom density to be assigned to the1. cell zone. (Always required in case of "dummy" elementsadmixed to the fuel elements!)
< 0.: Fractions are derived from data input of cards D (only if NHOM= 0oncardDl) .= -1. : Nuclide distributed like fuel.= -2.: Nuclide distributed like moderator.
OnlyifVB(l)>0.:Fraction of the homogenized atom density to be assigned to the L-thcell zone. L = 2 ... NCZ.
The fraction of the nuclide assigned to the coated particle fuel zoneICOAT on card T8 must be further subdivided between kernel.coating and matrix. VB(NCZ+1) gives the fraction in the kernels.
Tl 1 only if MUPN > 0 on card T10.
1
SFMU(K),K=1,MUP
Format (6E12.5)
Individual selfshielding of this isotope in the MUP broader thermalenergy groups as defined on card(s) T2.
T12 only if TLEAK = -1. on card T7.
Card T12
1 ALBEDO(l)
Format (6E12.5)
Albedo at the outer edge of the cell for the lowest energy group no. 1.
61
Continuation of card T12
2
3
30
ALBEDO(2)
ALBEDO(J),J=3,30
Albedo for the group no. 2 .= 0.: Use ALBEDO(1) for all energy groups.* 0.: Read Albedos for all groups.
Skipped if ALBEDO(2) = 0. Otherwise the group dependent Albedosmust be given.
Card T13 only if ITY = -1 on card Tl.
CardT13 Format (4E12.5)
1
2
3
4
FF(1)
FF(2)
FF(3)
FF(4)
Volumetric filling fraction of fuel elements in the core.
Inner radius of the fuel zone of the elements, (cm) (normal = 0.)
Outer radius of the fuel zone, (cm)
Outer radius of the element, (cm)
62
4.4.6 Diffusion calculation. Cl - C21
Cards Cl - C21 only if item JSER < 3 on input card V6 !
4.4.6.1 Title card
CardCl Format (18A4)
1
181=1,18
Literal description of case.
4.4.6.2 General control. C 2 - C 6
CardC2 Format (13)
IOPT 001
Control options.
CardC3
1
2
3
NGC10
NGC15
NGC24
Format (313)
Type of eigenvalue problem.= 0: Effective multiplication factor calculation.= -5: Fixed source (read cards Cl 8 - C21).
Termination option (applied only to the flux iteration calculation).= 0: Terminate calculation and proceed as if converged if machine
time or iteration count is exceeded (see also card C5).= 1: If limits are exceeded, terminate calculation and proceed as if
converged only if the iterative process is converging.= 2: If limits are exceeded, terminate calculations.
= 0: No effect.= -1: Define - possibly unisotropic - diffusion constants on cards
C11-C17.
63
Edit options.
Card C4
1
2
3
4
5
6
7
IEDG3
IEDG4
EDG5
IEDG6
IEDG9
EDG10
IEDG12
Format (713)
= 0: No effect.> 0: Print macroscopic group-to-group transfer cross sections.
= 0: No effect.> 0: Print macroscopic reaction rate cross sections.
= 0: No effect.> 0: Print gross neutron balance over system by group.
= 0: No effect.> 0: Print gross neutron balance by zone by group.
= 0: No effect.> 0: Print zone average flux values by group (IEDG6 = 0).
= 0: No effect.= 2: Print point flux- and point power density values, write them onto
formatted data sets 'phiform' and 'powform', respectively.
= 0: No effect.> 0: Print zone average power densities.
General iteration count and machine time limit.Problems are terminated when the iteration count reaches the limit and the calculation proceeds asper NGC15 (see card C3).
Card C5
1 LTMX1
Format (313)
> 0: Maximum number of initial eigenvalue problem iterations.(< 999)
= 0: Default value = 200
The following items are machine time limits (min). Generally, calcu-lations continue if time is exceeded as if convergence criteria had beensatisfied.
64
Continuation of card C5
ITMX19
ITMX20
> 0: Limit for the initial eigenvalue problem.= 0: Default value = 60
> 0: Limit for all other eigenvalue problems.= 0: Default value = 30
General restraints.
CardC6
1
2
GLIM1
GLIM2
Format (2E 12.5)
Any calculation will be terminated if the following restraints are notmet.
> 0.: Maximum multiplication factor.= 0.: Default value = 1.5
> 0.: Minimum multiplication factor.= 0.: Default value = 0.1
4.4.63 Description of neutron flux problem. C7 - C10
CardC7
IOPT
Format (13)
003
65
General description.
Card C8
1
2
3
4
5
6
7
8
9
NU AC 11
NU AC 12
NUAC13
NUAC14
NUAC15
NUAC16
NUAC17
NUAC18
NUAC19
Format (1116)
Note:2-dim.: r = left —> right
z = top —> bottom
3-dim.: <D-r-z or x-y-z<f> or x = left —> right
r or y = top —» bottomz = front —> back
Left boundary condition (required for 2-d, 3-d).= -1: Periodic.= 0: Extrapolated (vacuum).= 1: Reflected.
Top boundary condition (required for 2-d, 3-d).= 0: Extrapolated.= 1: Reflected.
Right boundary condition (required for 2-d, 3-d).Set to-1 ifNUACll =-1= 0: Extrapolated.= 1: Reflected.
Bottom boundary condition (required for 2-d, 3-d).= 0: Extrapolated.= 1: Reflected.
Front boundary condition (required for 3-d).= 0: Extrapolated.= 1: Reflected.
Back boundary condition (required for 3-d).= 0: Extrapolated.= 1: Reflected.
Number of zone to be an internal black absorber and to have the non-return boundary condition applied at its edges (see XMIS2 on cardC10; this zone will be black to all groups unless additional data aresupplied).
= 0: Only positive neutron flux allowed.> 0: Option to allow negative neutron flux.
= 0: No effect.> 0: Override use of Chebychev polynomials in adjusting the accel-
eration parameters.
66
Continuation of card C8
10
11
NUAC20
NUAC23
= -1: Force alternating direction line relaxation on rows and columns,and also fore and after for 3-d.
= -2: Use only on rows and columns.> 0: Line relaxation only on rows.= 0: The code selects line relaxation on rows only with one inner
iteration for all problems involving upscattering, otherwise threeinner iterations for 3-d problems without I/O and five with dataI/O during iteration, and alternating direction line relaxation forall 2-d problems.
Number of inner iterations. Normally not specified (see NUAC20above).
Iteration convergence criteria
CardC9
1
2
EPI1
EPI2
Format (2E12.5)
> 0.: Maximum relative flux change for the last iteration of eachinitialization eigenvalue problem.
= 0.: Default value = 0.0001
> 0.: Maximum relative change in the eigenvalue for the last iterationof eigenvalue problems. This applies to the multiplication factorcalculation.
= 0.: Default value = 0.00001
Miscellaneous data
Card C10 Format (3E12.5)
XMIS1 External extrapolated boundary constant.= 0.: The code will use the built-in value for all extrapolated
boundaries. (0.4692)> 0.: Specifies the constant for all extrapolated boundaries for all
groups (see NUAC11 - 16 on card C8).
67
Continuation of card CIO
2
3
XMIS2
XMIS6
Internal black absorber boundary constant for the zone NUAC17.= 0.: In connection with NUAC17 > 0 on card C8 the code will use
the built-in value for all groups and the absorber will be blackover all energy. (0.4692)
> 0.: The constant for all groups applying to zone NUAC17.
Initial overtaxation factor. Normally calculated by the code and notspecified here.
£
4.4.6.4 Simulation of control devices and void areas. C11-C17
Cards C11 - C17 only if NGC24 = -1 on card C3
CardCll Format (216)
1
2
JH
KH
Number of regarded areas (card(s) C12). (< 200)
Number of different cross section sets to be applied (cards C13C17). (<20)
CardC12 Format (1814)
1
2
3
IZONE(J),
M1(J,2),J=1,JH
Id. no. of cross section set to be applied to the J. area.
First CITATION zone located in the J. area.
Last CITATION zone located in the J. area.
68
One set of cards C13 - C17 for each of the K = 1,KH cross section sets.
CardC13 Format (6E 12.5)
1
N26
RDK(K,I),1=1,N26
Diffusion constants of the energy groups I.
Card C14 Format (6E12.5)
1
N26
SGA(KJ),1=1, N26
Macroscopic absorption cross sections of the energy groups I.
Card(s) C15 for each energy group I = 1,N26.
Card C15 Format (6E12.5)
1
N26
SGTR(K,I,J),J=1,N26
Macroscopic transfer cross sections from energy group I to energygroups J.
CardC16 Format (3E12.5J3)
1
2
3
4
V2(K,1)
V2(K,2)
V2(K,3)
IKEN
V2(K,L), L = 1, 3 : Factors to be multiplied to the diffusion constants
2-d (r - z) : in r- direction.3-d (4> - r - z): in O - direction3-d (x - y - z) : in x - direction
2-d (r - z) : in z - direction3-d (O - r - z) : in r - direction3-d ( x - y - z ) : iny- direction
2-d (r - z) : =0.3-d (<I> - r - z) : in z - direction3-d ( x - y - z ) : in z - direction
= 0: No effect.> 0: Group dependent factors will be defined on card C17.
69
Card C17 only if IKEN > 0 on card C16.
CardC17 Format (6E 12.5)
1
N26
FKEN(KJ),I=],N26
Energy group dependent factors to be multiplied to the V2 of cardC16.
4.4.6.5 Fixed source, specified by zones. C18 - C21
Cards C18 - C21 only if NGC10 = -5 on card C3
CardC18 Format (13)
IOPT 026
Card C19 Format (13)
NFX2 = 0: Short output.> 0: Source (n/sec) will be edited by mesh points.
Card C20 Format (6E 12.5)
1
N261=1, N26
Fractions of the fixed neutron source distributed into each groupstarting with the highest energy group. These should sum to unity butare normalized to unity by the code.
70
One card C21 for each zone having a fixed neutron source. Fixed source input is terminatedby a'blank'card C21.
Card C21
1
2•
N2F(I),
V2F(I),1=1
Format (6(I3,E9.3))
> 0: Zone number.= 0: End of the 'fixed source' input
Fixed source, (n/sec-cm3)
71
4.4.7 Fuel cycle costs calculation. Kl - K12
Cards K1 -K12 only if NKOST > 0 on card V6, or - in case of a Restart - if NEWCOST = 1 oncard R6, then cards K1-K12 following card R6.
CardKl
1
2
3
4
5
6
7
8
9
10
NMAF
MXTYP
ND2O
NPUFD
IPRINT
1Q
NEWCO
$
$X$
$s
Format (8X,7I4,A4,I4,A4)
For the cost calculation the length of approach to equilibrium phase isassumed to be equal to NMAF burnup cycles.NMAF < MMAF (card VI).
Number of different fuel types in the system (< 10). For each type aset of cards K7-K10 is required.
= 0: Normal.> 0: Heavy water moderated reactor. D2O expenditures included in
cost calculation. Card K12 required. (This option may be used tosimulate capital costs of power plant).
= 0: Normal.> 0: Pu feed cycle, for each period a Pu equivalence value is calcu-
lated according to specified FCC for uranium feed cycle on cardKl 1.
= 0: Print-out without materials balance for each batch.= 1: Print-out includes materials balance.
= 0: Normal.> 0: Calculate average equilibrium FCC over the last IQ periods
to obtain representative FCC for the equilibrium cycle in caseit consists of more than 1 period (Ref. /20/, Section 9).
= 0: Neglect financing cost of fresh out-of-pile batches. (Normalwhen appropriate lead-times are used).
= 1: Calculate financing cost of fresh out-of-pile batches for the timeTOUT on card K4. This option only for cases with out-of-pilebatches, i.e. KUGL > 0 on card R2.
Monetary unit in which input is supplied. The user specifies the 4character alphanumeric designation to be used in print-out, e.g. EURor US$.
The energy cost is calculated in units of $$ (see below) and $X$ is theconversion from $ to $$. E.g. 100 means $ = 100 $$.
Monetary unit in which energy costs are calculated and printed inoutput, e.g. Cent.
72
CardK2
1
2
3
4
F
ETA
GLD
GMZ
Format (4E12.5)
Annual load factor. Same as AAAA on card R2.
Net efficiency of power plant.
> 0.: Total lifetime of the power plant (a). Average FCC are calcu-lated for GLD years assuming an approach to equilibrium phaseof NMAF periods. For the rest of the lifetime the last (or the IQlast) calculated periods are defined to be the equilibrium periodand repeated till the end of plant operation. For D2O cost calcu-lation GLD is taken as amortization time for heavy water invest-ments.
< 0.: Drop average FCC calculation.
Number of installments of electricity revenues within a period. Nor-mally = 1., but for longer operation periods monthly or quarterlyintervals of payment should be assumed.
CardK3
1
2
3
4
5
6
Zl
Z2
Z4
Z3
Z5
ZL
Format (6E12.5)
Pre-irradiation interest rate (I/a) on fuel expenditures.
Pre-irradiation interest rate (I/a) on fuel fabrication and D2O re-placement costs.
Interest rate (I/a) on all capital, incl. electricity revenues duringirradiation.
Post-irradiation interest rate (I/a) on capital to finance fuel credit (ineffect discount rate).
Post-irradiation interest rate (I/a) on reprocessing and shipping costs.
Discount rate (I/a) for present worth leveling of all expenditures andrevenues over reactor lifetime. In most cases all interest rates will bechosen the same with the possible exception of the present worthdiscount rate. The code offers the flexibility to model most of theeconomic situations arising for those special cases where this might beneeded.
73
Card K4
1
2
3
4
5
SS
RES
VERL
YPA
TOUT
Format (5E 12.5)
Tax rate (I/a) on fissile investments.
Reserve factor (RES = 1. + reserve) to account for additional fabri-cation costs for reserve elements in the initial core. In later cyclesprogram sets RES = 1. The capital charges arising from a reservestore are contained in the appropriate defined lead-times TIN andTFAB on card K9. Blank = 1.0 .
Recovery factor for reprocessing. (0.97 - 0.99)
Fraction of discharged 233Pa decaying into 233U during out-of-pilestorage (normally = 1.0). If storage and reprocessing time are definedon card R2, the amount of 233U reaching the reprocessing plant hasalready been explicitly accounted for, then YPA = 0. (239Np isassumed to decay completely into 239Pu).
Storage time (days) before reuse of out-of-pile batches. Financingcost with interest rate Z4 is calculated during time TOUT. Also forfresh fuel if not NEWCO on card Kl is specified = 0.TOUT = -1. will cause the code to specify TOUT = cycle length foreach cycle.
74
Card K5 is always required. The data for uranium ore and enrichment will be used to calculate theprice of 235U for different fuel types and the changing value during depletion. If this option is to beby-passed, CU8(K) for all fuel types > 0. and CU5(K) specified accordingly, see cards K7.
CardK5
1
2
3
4
5
6
CU3O8
CO8F6
CTRENN
TAIL
XLOSS1
XLOSS2
Format (6E12.5)
Cost of uranium ore as U3O8 ($/lb U3O8). The price is given as per lbcorresponding to common use in literature.
Cost of conversion of U3O8 to UF6 ($/kg U). The enriched end productis in the form UFÖ, the costs of converting the hexafluoride into UO2
or any other compound, are included in the fabrication costs.
Separation cost. ($/SWU)
Tail enrichment, i.e., 235U content in discarded uranium from enrich-ment plant.
Fraction of losses in conversion of U3O8 to UF6 (typically 0.005 -0.01).
Fraction of losses in conversion of enriched UF6 to UO2 or UC and infabrication (typically 0.005 - 0.01).
Card K6 is always required. The costs of fresh 232Th, 233U and fissile plutonium are assumed tobe the same for all types of fuel. The discharge value may, however, vary according to composi-tion and subsequent utilization and for each fuel type depreciation factors are specified on cardsK8.
CardK6
1
2
3
CTH232
CU233
CPUFIS
Format (3E12.5)
Costof232Th. ($/kg)
> 0.: Cost of fissile 233U. ($/kg)< 0.: Cost of 233U is calculated relative to cost of 93% enriched 235U
with |CU233| as parity value.Cost(233U) = |CU233| * Cost(93% 235U).
> 0.: Cost of fissile 239Pu and 24IPu. ($/kg)< 0.: Cost of Pufiss relative to cost of 93% 235U with parity value
|CPUFIS|.
75
One set of cards K7 - K10 for each fuel type.
Card K7
1
2
3
4
5
Card K8
1
2
ANSM
CU5
CU8
CFAB
CAUF
CHITH
CHIU3
Format (5E 12.5)
Type of heavy metal:l.:Th(met.), 2.: ThO2, 3.: ThC, 4.: ThC2
5.: U(met.), 6.: UO2, 1:. UC, 8.: UC2
IfCU8 (next variable) <0.:Initial reference enrichment of 235U in uranium. All cost calculationsare performed with the actual enrichment of a batch regardless of thereference enrichment for the type. Cost data on card K5 are used.IfCU8>0.:Cost of 235U ($/kg). Price kept constant during calculation.
< 0.: Cost of 238U = 0., and cost of 235U calculated from batch enrich-ment and card K5.
> 0.: Cost of 238U ($/kg). Supply 235U cost as specified above. Thisoption operates only if CU8 > 0. for ajl types!
> 0.: Fuel fabrication cost ($/kg HM) excluding cost of heavy metal.Monetary unit is variable $ as specified on card Kl and given perkg initial HM in fuel element.
< 0.: Data for this fuel element type is calculated by the code itselfusing basic cost data as specified on card D3.
> 0.: Total costs of reprocessing, shipping and storage ($/kg HM)payable at time TAUF (card K9) after discharge. Cost per kgdischarged HM.
< 0.: Data for this fuel element type is calculated by the code itselfusing basic cost data as specified on card D4.
Interests on HM during fabrication and reprocessing are calculatedseparately by the program for lead and lag times TIN and TEX oncard K9.
Format (4E12.5)
= 0.: Irradiated and discharged fertile material (Th, U) has no value.> 0.: Cost of discharged fertile material depreciated by the factor
CHITH.
= 0.: Discharged 233U has no value.> 0.: Discharged 233U cost depreciated by the factor CHIU3.
76
Continuation of card K8
3
4
CardKS
1
2
3
4
5
CHIU
CHIPU
TORE
TIN
TFAB
TEX
TAUF
CardKlO
1
2
3
TORE
TIN
TFAB
= 0.: Discharged 235U has no value.> 0.: Discharged 235U cost - for actual enrichment in depleted fuel -
depreciated by the factor CHIU.
= 0.: Discharged fissile 239Pu and 241Pu has no value.> 0.: Discharged Pu^s cost depreciated by the factor CHIPU.
Format (5E12.5)
Lead-time (d) for payment of uranium ore for replacement fuel. Lead-time is counted prior to the time of loading the fuel into the reactorand start of irradiation.
Lead-time (d) for payment of enrichment service and conversion costsfor fuel replacement relative to fuel loading. Also lead-time for pur-chase of 233U and fissile plutonium.
Lead-time (d) for payment of fabrication costs for replacement fuelrelative to fuel loading. The lead-time for D2O replacement is takenthe same as TFAB for fuel type 1.
Lag time (d) for credit for discharged fuel relative to time of dis-charge at end of irradiation. No difference between lag times for re-placement and initial fuel.
Lag time (d) for payment of reprocessing and shipping costs for dis-charged fuel relative to end of irradiation. Same for replacement andinitial fuel.
Format (3E12.5)
Lead-time (d) for payment of uranium ore for initial core.
Lead-time (d) for payment of enrichment service and conversion costsfor initial core.
Lead-time (d) for payment of fabrication costs for initial core.
In general the lead-times for purchase of initial core will be longerthan for replacement fuel as the amount of ore and the number ofelements to be manufactured are larger.
77
Card Kl 1 only if NPUFD > 0 on card Kl.
Card Kl1 Format (6E12.5)
NMAF
UUKOST(1)
UUKOST(NMAF)
Fuel cycle cost ($$/kWh) for the first period for the correspondinguranium feed cycle. A Pu price is evaluated for this period to yieldFCC equal to UUKOST(l).
The FCC for the U-cycle must be specified for all periods NMAF.This calculation allows for NMAF < 30 only.
Card K12 only if ND2O > 0 on card K1.
Card K12
1
2
3
4
5
6
DOUTIN
CDNEU
CDALT
ZD
SD
VD
Format (6E12.5)
Weight (kg) of total D2O in-pile and out-of-pile inventory. (D2O costcalculations may be used to simulate amortization of capital invest-ment for the plant. The capital is then paid in installments at the be-ginning of each cycle or accounting period, in such a way that the in-stallments leveled over the lifetime GLD give the total present worthvalue at the time of start-up. In this case DOUTIN could be interpretedas kWe of the power plant.)
Cost ($/kg) of new D2O. (Capital cost including interest during con-struction at time of start-up in $/kWe.)
Cost ($/kg) of old D2O (0., i.e. no value of station end-of-life).
Interest rate (I/a) on D2O investments (or capital costs).
Tax rate (I/a) on D2O investments (or capital costs).
D2O losses per year (normally 0.01). It is assumed that the D2O re-placement expenditures have the same lead-time and interest rate asfabrication costs. (0.0 in case of plant cost).
78
4.4.8 Fuel management. Rl - R34
Cards Rl - R34 only if NRSTRT > 0 on card V6.
4.4.8.1 General definitions. Rl - R2
Card Rl only if NRSTRT = 3 or 4 on card V6, i.e. for fuel management with reprocessing.
CardRl Format (6E12.5)
KTOT
XREPRO(I),1=1,KTOT
Reprocessing factor for material no. I. The reprocessing plant is simu-lated by reprocessing factors multiplied to the different nuclide quan-tities. The decay of the heavy metal isotopes is calculated for periodTREPRO (card R2). The reprocessing factors are defined as:
1.00 = No losses.0.00 = Complete removal.0.95 = 5% loss during reprocessing, etc.
Data must be specified for all nuclides [KTOT = 28 (i.e. no. of heavymetal isotopes) + KMAT (card Dl) non-heavy metal nuclides]. Thesame reprocessing factors are applied to all batches.
CardR2
1
2
3
KUGL
TDOWN
TSTORE
Format (7E10.3)
= 0.: Reactor without out-of-pile cycle.= 1.: Pebble bed reactor with out-of-pile fuel management:
Fresh fuel stores for each type, reusable discharged fuelbatches, and handling of discarded scrap fuel.
= 2.: Reactor with out-of-pile cycle:As above, fuel may be replaced in any core position.
Length of downtime during reload (days). The isotopic decay of theheavy metals is calculated for all in-core batches.
Length of out-of-pile storage time (days) before reuse of the fuel.Decay is calculated for all out-of-pile batches, except for the freshand for the scrap fuel. Financing costs are paid during timeTSTORE.
79
Continuation of card R2
4
5
6
7
TREPRO
AAAA
BRUCH
AGEBOX
Length of cooling, shipping and reprocessing time (days) of dis-charged fuel. Decay is calculated during this period for scrap fuelbatches and reprocessed fuel batches.
Load factor of power plant. When fuel cycle costs are evaluated, sameas F on card K2 in cost input.
Failure rate of discharged fuel (pebble bed reactor only). In each dis-charged batch a fraction BRUCH is assumed non-reusable and isadded to the scrap batches.
= 0.: No effect.= 1.: Aging boxes for reprocessing mixtures will be specified. Card
R5 is required.
4.4.8.2 Data for individual fuel types. R3 - R4
Cards R3-R4 only if KUGL > 0. on card R2. The cards are supplied as a set for each fuel type,i.e. JTYP (card V1) sets.
Definition of fresh fuel store.
Card R3 Format (2(I6,E 12.5))
1
2
NTP1
PARVL1
NISO
Id. no. of the fuel type.
Volume (cm3) of fuel store. The choice of fuel volume is arbitrary aslong as no financing costs of fresh fuel store are calculated. In manycases it is advantageous to define the volume equal to the volume ofone fuel element and when reloading specify the fraction of the storeas no. of elements in the batch. If more fuel from the store thanactually present is required, the store is regarded as unlimited. Fuelremoved from the store does not change the remaining volume,neither does the isotopic composition change during the reactor life.
> 0: Number of isotopes on the following card(s) R4. The compositionof fresh fuel is specified by the number NISO and by the card(s)R4. The rest of the isotope concentrations equals zero.
80
Continuation of card R3
4
Card(s)
CardR4
1
2
3
4
XMARX
R4onlyifNIS0
L
DAV(L)
L
DAV(L)
= 0: Card(s) R4 are skipped for this fuel type, and PARVL1 = 0.< -100: Isotope composition of the fresh fuel store has been defined
by variable NFUTP on card D6.|NISO| = NFUTP.
No. of the reprocessing mixture to which the scrapped fuel of this typeis transmitted. Each mixture may consist of one or more fuel types.After each reload the discharged fuel is volume averaged to form amixture, which at the next reload may be reprocessed and be used forrefueling. XMARX < 10.
If AGEBOX > 0. (card R2) this reprocessing mixture is transferred toits corresponding aging box(es). If for all types XMARX = 0., no re-processing mixtures are prepared.
> 0 on card R3.
Format (4(I6,E12.5))
Id. no. according to VSOP list of first nuclide with atom density * 0.
Atom density of nuclide L in fresh fuel store for type NTP1.
Id. no. of second nuclide.
Atom density of second nuclide.
Data for all NISO isotopes in fresh fuel store.
81
4.4.8.3 Aging boxes for discharged fuel. R5
Card R5 only if AGEBOX > 0. on card R2.
Card R5 Format (1014)
MREP
NAJB(J),J=1,MREP
Number of aging boxes including one jumble box to be defined for theJ-th reprocessing mixture. MREP is the total number of reprocessingmixtures defined by the XMARX sequence on card R3.NAJB < 10.
Note:
1. Scrap fuel discarded from the reactor is loaded into the repro-cessing mixture box J.
2. It is transferred to the corresponding first aging box.3. Aging boxes are stepwise transferred to the next higher ones.4. Those with an age > TREPRO (card R2) are transferred to the
corresponding jumble box J.5. If NAJB(J) = 1, the reprocessing mixture box J is immediately
given to the jumble box.6. A fraction FOJB(J) (card RIO) of the jumble box is loaded into the
reprocessing mixture box J. It is ready for use after that reload,which will be performed after the following burnup cycle.
7. That fraction which is not used, is returned to the jumble box J.
4.4.8.4 Instructions for the burnup cycles. R6 - R27
These cards will be read at the end of each burnup cycle. They define the fuel management priorto the subsequent cycle and give some new options for the next cycle.Card R6 only if IRR9 > 0 on card S3, and only at the beginning of a restart. The preceding runended after a fuel management performance. The restart starts at the beginning of the new burnupcycle. This card allows to change some options for this first cycle, which were given at the lastcard R7 of the preceding run.
Card R6 Format (613,11,I2.5I3.3E 12.5)
1
2
3
4
IPRIN(l)
IPRIN(2)
IPRIN(3)
IPRIN(4)
Same as on card R7.
Same as on card R7.
Same as on card R7. (-2: Doesn't work).
Same as on card R7.
82
Continuation of card R6
5
6
7
8
9
10
11
12
13
141516
NNSTOP
NNUM
NEWCOST
NIAVC
IBUC
MUHU3
NOCPA
rvspn
INGC24
XDAYXPOWXKAY
= 0: No effect.> 0: Number of large time steps per burnup cycle, i.e. redefinition of
JNSTOP (cards VI1 and R9).
= 0: No effect.> 0: Number of small time steps per large time step, i.e. redefinition of
JNUM (cards VI1 and R9).
= 0: No effect= 1: Redefinition of cost data (only if NKOST > 0 on card V6):
Read cards Kl - K12, following card R6.
Has no meaning, if preceding item NEWCOST = 1 != 0: No change of the option of average fuel cycle cost calculation.= 1: Drop average FCC calculation.= -1: Calculate average FCC.
= 0: Leakage feed back option unchanged.= 1: Feedback of broad group bucklings to GAM-I and thermal
leakageto THERMOS.
= 2: Feedback of average epithermal buckling to GAM-I and thermalleakage to THERMOS.
= 3: No feed back at all.
= 0: Streaming correction option unchanged.= 2: Streaming correction LIEBEROTH /28/ in power generating
batches (only for pebble bed).= 3: No streaming correction at all.
= 0: Option of control poison adjustment unchanged.< 0: If in the preceding run control poison adjustment was calculated,
it can be stopped here.
= 0: No change of diffusion calculation option.< 0: Drop diffusion calculation.> 0: Repeat diffusion calculation as defined by the IDIFF on cards
V17/R12.
= 0: No effect.> 0: Value of option NGC24 (if previously defined equal -1 on
card C3) is reset to zero value, starting with the next diffusioncalculation.
Same as DELDAY on cards V10 and R9.Same as POWER on cards V10 and R9.Same as ZKFIND on cards V10 and R9.
83
Card R7
1
2
3
4
5
6
7
IVSP(1)
IPRIN(l)
IPRIN(2)
IPRIN(3)
NPRINT
IPRINT
IN2
Format (6I3,3X,3I3,6I2,2X,9I2,110)
< 1: The information of this card holds only for this shuffling and forthe following burnup cycle.
> 1: The information of this card holds for IVSP( 1) shufflings andburnup cycles. The items 2 ... 10 are kept, the others are set to 0.
Output options:
Spectrum calculation (same as on card VI5):= -1: Minimal output.= 0: Thermal selfshielding factors, only.= 1: Same as 0, plus averaged thermal cross sections.= 2: Same as 1, plus fine group neutron fluxes.= 3: Same as 2, plus broad groups averaged cross sections for
materials with concentration > 0.= 4: Same as 3, for all materials.= 5: Maximum output including details of neutron transport.
= 0: No output.= 1: Printout of the irradiation time of the batches.= 2: Same as 1, plus atom densities (only in combination with
IPRIN(3) > 0).
Burnup calculation (same as on card VI5):= -2: All output dropped except Keff.= -1: Global neutron balance.= 0: Detailed neutron balance.= 1: Same as 0, plus characteristics of all batches.
Fuel management operations at this reload:= -1: No fuel management output.= 0: Short summary.= 1: List of all operations.= 2: Detailed printout including atom densities in all batches before
and after reload (very much!).
Fuel cycle costs calculation (same as on card Kl):= 0: Printout without materials balance for each batch.= 1: Printout includes materials balance.
= 0: No effect.> 0: No. of the burnup time step for which a list of performance data
will be printed.< 0: Print neutron balance averaged over the cycle.
84
Continuation of card R7
10
11
12
IPRIN(4)
IVSP(ll)
NCYC
NKEEP
IK
Steering the calculation performance:
= 0: Skip spectrum calculation.= 1: Repeat spectrum calculation as defined on cards V16 / R11.
= 0: Drop diffusion calculation.> 0: Repeat diffusion calculation as defined by the IDIFF on cards
V17/R12.= 2: Read card R32 for this fuel shuffling in order to redefine the
CITATION edit options. (Does not work for the first diffusioncalculation of a restart and for burnup cycles with THERMIX-calculation included).
= 0: No effect.= 1: After this fuel management some fuel of the jumble boxes can be
loaded into reprocessing mixtures. Here it is available for the fuelmanagement after the next burnup cycle. Read card RIO.
= 0: Use the previously defined fuel management scheme.= 1: Read a new fuel management scheme on the cards R21.= 2: Generate a new fuel management scheme with all batches staying
in their position (no cards R21 are required).= 3: Generate a new fuel management scheme. Batches with indivi-
dual fuel management instructions will be identified on cardsR21. Non identified power generating batches will be shuffled tothe next region. Non-power generating batches (e.g. reflectors)stay in their position.
= 4: Same as 3. Also the non-identified power generating batches stayin their position.
= 5: Generate a new fuel management scheme. For power generatingbatches: Fuel management instructions given for a batch on cardR21 are also valid for all subsequent batches until redefined by anew card R21. Non-power generating batches (e.g. reflectors)stay in their position.
= 6: Same as 3, but shuffling of each non identified power generatingbatch to the following batch within the same region.
= 0: Normal.= 1: Fuel management is preserved for ORIGEN after large time
step 1 (equilibrium cycle treatment for MEDUL-reactors - seeinput description of ORIGEN-JUEL-II151) on data set 'origmod'.Note: The value of variable JNUM (see card R9) is changed to _LIt must be redefined for the following burnup (shuffling-) cycles,if desired.
85
Continuation of card R7
13
14
15
16
17
18
LK
NTIK
NJ
IVOID
IREDE
IVSP(16)
= 2: Power histogram is preserved for decay heat calculation inthermal hydraulics part (THERMDC) or NAKURE-code /21/starting from next cycle (data set 'nakure'). IK must be defined onthe very first card R7. if on card LF2 variable KT5 shall be '=0'!
= 3: Stop the preservation of the data for decay heat calculation. Lastpreserved data are from preceding cycle. Read cards LF after thevery last card R (last burnup cycle, IVSP(24) > 0 on this card R7).
= 0: Normal.= 2: Life history is preserved for ORIGEN-JUEL-II (all NON-
MEDUL-reactors, see /5/), starting from next cycle (data set'origen').
= 3: Stop the preservation of the data for ORIGEN-JUEL-II. Lastpreserved data are from preceding cycle. Read card P after thevery last card R (last burnup cycle, IVSP(24) > 0 on this card R7)and cards LF, if present.
= 0: No effect.> 0: Perform temperature calculation:= 1: Read new time steps ITEMP on card R13, and read new input for
THERMIX on cards TX.= 2: Read new ITEMP on card R13, use previous THERMIX input.= 3: Use previous ITEMP, and give new THERMIX input.= 4: Use previous ITEMP, previous THERMIX input.
OnlyifNTIK>0:= 1: One single THERMIX calculation at each time step given by the
HEMP (steady state).> 1: Time dependent THERMIX parallel to the VSOP (cards TX, see
also Fig. 10):= 2: Power for THERMIX is only the decay heat.= 3: Power for THERMIX is the decay heat plus the fission power of
the individual time step.
Redefinitions:
= 0: No effect.= 1: Redefinition of void areas for CITATION on card R8.
= 0: No effect.= 1: Redefinition of time steps, power, criticality constraints etc. on
card R9.
= 0: No effect.= 1: Redefinition of time steps for spectrum calculation. Give
ISPEKT on card Rl 1 (same as card V16).
86
Continuation of card R7
19
20
21
22
23
24
25
IVSP(17)
IRETEM
LIB
MUHU(l)
IWATER
IREDEF
IVSP(24)
= 0: No effect.= 1: Redefinition of time steps for diffusion calculation. Give IDlhh
on card R12 (same as card V17).
= 0: No effect.= -1: Redefinition of the temperature of the resonance absorbers or of
a scattering nuclide in one or more regions - changed bythe same value in each zone.Read card R16, drop cards R14 + R15.
= 1: Redefinition of the temperatures of all the regions. Read cardsR14-R15 (same as card G2 and cards T4), drop card R16.
= 2: Read new resonance integral definition on cards R17-R18 (sameas cards G4-G5).
= 3: Includes both options 1 and 2.
= 0: No effect.> 0: Write data set 'tinte' ("status of core") for TINTE /29/ at time
step LIB. Read card R34.< 0: For each batch write atom densities on data set 'nucdens'.
= 0: No effect.> 0: Extracted number of nuclides (< 20) for the output of atom
densities at the end of the burnup cycle. Read id. numbersNUPRI(I) on card(s) R33.
= 0: No effect> 0: Calculation of water ingress to core and/or reflectors. Read
cards R20A and R20B.
= 0: No effect.> 0: Read new identification numbers of THERMOS-cell definitions
for the regions on cards R7a (like card T6).
= 0: No effect.> 0: Terminate the run after this fuel shuffling. The following burnup
cycle will be the first cycle in a restart case.
Card(s) R7a like card(s) T6 !
87
Card R8 only if IVOID = 1 on card R7 (only for the CITATION diffusion calculation).
CardR8
1
2
JH
IZONE(I),
M1(I,1),
Ml (1,2),1=1,JH
Format (1814)
Number of void areas. For each void area the following set of items:
Id. no. of a void cross section set to be inserted in the I-th void area(the cross section sets are defined on the cards C13 - C17 ofCITATION. - see section 4.4.6.4).
First CITATION region located in the I-th void area.
Last CITATION region located in the I-th void area.
Card R9 only if IREDEF = 1 on card R7.
CardR9 Format (3I4,3E12.5,4E6.0)
1
2
3
4
5
6
7
JNSTOP
JNUM
IVSP(27)
DELDAY
POWER
ZKFTND
HNUC
= 0: No effect.> 0: Redefinition of the number of large burnup time steps per burnup
cycle (see card VI1).
= 0: No effect.> 0: Redefinition of small burnup time steps in one large step (see
card VI I).
= 0: Streaming correction as defined before.= 2: Streaming correction by LIEBEROTH /28/.
= 0.: No effect.> 0.: Redefinition of the length of one large time step (days) (card
V10).
= 0.>0 .
= 0.>0.
= 0.>0.
No effect.Redefinition of thermal core power, Watts (card V10).
No effect.Redefinition of end of cycle - Keff (card V10).
No effect.Number of new atom densities to be read on card R19A. (< 12)
Continuation of Card R9
88
8
9
10
HPOS
XTDOWN
CONPOI
OnlyifHNUC>0.:> 0.: Number of batches to be loaded with the new atom densities
(read card R19B).= 0.: Load additional materials into all power generating batches.
= 0.: No effect.> 0.: Redefinition of length of down time during reload (TDOWN on
card R2).< 0.: Set TDOWN = 0.
= 0.: No effect.> 0.: Read control poison adjustment data on cards R24 - R27.
CONPOI gives the number of regions with control poison data.< 0.: Stop the control poison adjustments.
Card RIO only if NCYC = 1 on card R7.
Card RIO Format (6E12.5)
MREP
FOJB(l),1=1,MREP
> 0.: Volumetric fraction of the I-th jumble box which shall be loadedinto the reprocessing mixture I for the use at the reload after thefollowing burnup cycle.
< 0.: Volumetric fraction is calculated by the code. |FOJB(I)| gives theratio: Volume to be loaded into the reprocessing mixture I /volume of the scrap fuel, which is discharged to the first agingbox.
Card Rl 1 only if IVSP(16) = 1 on card R7.
CardRll
1
2
18
ISPEKT(l)
ISPEKT(I),1=2,18
Format (1814)
> 0: No. of the first large burnup time step in which the spectrum cal-culation is to be repeated prior to the diffusion calculation.
> 0: No. of further time steps for spectrum calculation.= 0: If all ISPEKT = 0, spectrum calculation is performed in every
time step.
Card R12 only if IVSP(17) = 1 on card R7.
89
Card R12 Format (1814)
1
18
IDIFF(I),1=1,18
If all IDIFF(I) = 0:Diffusion calculation is performed at every time step.
If at least one IDIFF(I) * 0:The IDIFF(I) give the time steps at which diffusion calculation is to beperformed.
Card R13 only if NTIK = 1 or 2 on card R7.
Card R13 Format (1814)
1
18
rrEMP(I),1=1,18
IfalHTEMP(I) = O:THERMDC temperature calculation at every time step.
If at least one ITEMP(I) * 0:The ITEMP(r) give the time steps at which temperature calculation isto be performed.
A set of cards R14 - R15 only if IRETEM = 1 or 3 on card R7.
CardR14 Format (6E 12.5)
1
NDR
SEMZUT(I),1=1,NDR
> 0.: Temperature of the resonance absorbers in the NDR differentregions. (°C).
= 0.: Temperature for region I stays unchanged.
For each of the NKER scattering nuclides (see card Tl) one set of cards R15.
CardR15 Format (6E12.5)
1
NDR
SCELS(I),1=1,NDR
> 0.: Temperature of this scattering nuclide for region I (°C).
= 0.: Temperature stays unchanged.
90
Card R16 only if IRETEM = -1 on card R7.
Card R16
1
2
3
4
5
NVAR
NXSA
NXSE
TV AR
TMIN
Format (3I4.2E 12.5)
- 1: Temperatures of the resonance absorbers are changedby AT = TV AR (see 4th item on this card).
= 2: Temperatures of scattering nuclide 1 are changed by AT =TV AR.
= 3: Like NVAR = 2, but for scattering nuclide 2= 4: Like NVAR = 2, but for scattering nuclide 3= 5: Like NVAR = 2, but for scattering nuclide 4= 6: Like NVAR = 2, but for scattering nuclide 5
ID-number of the first region to get a modified temperature.
ID-number of the last region to get a modified temperature.
Value of temperature change.
Lowest occurring temperature is limited to TMIN.
Cards R17 - R18 only if IRETEM = 2 or 3 on card R7.For each design (IDESIN on card Gl) 4 sets of cards R17-R18, the first set for 232Th, the secondset for 238U, the third set for 24OPu, the fourth set for 242Pu.
Card R17
1
2
3
SMI
SM2
NZ
Format (2E 12.5,16)
Two values of homogenized atom densities of the resonance absorbernuclide, for which sets of resonance integrals are available on data set'resint'. These values should represent the highest and the lowest den-sities, occurring within the reactor, respectively.
> 0.: Highest density of the absorber nuclide. [barn1 cm"1]= -1.: Density is taken from data set 'resint'.
> 0.: Lowest density of the absorber nuclide. [barn1 cm"1]= -1.: Density is taken from data set 'resint'.
Number of sets of resonance integrals for each SMI and SM2 den-sities. These sets represent different temperatures of this absorbernuclide.
91
2 cards R18, the first one corresponds to SMI, the second one to SM2.
CardR18 Format (1216)
1
NZ
IZUT(K),K=1,NZ
Id. numbers of the resonance integral sets to be read from data set'resint'.
Card R19A only if HNUC > 0. on card R9.
Card
1
2
R19A
INEW(I),
DNEW(I),1=1,HNUC
Format (4(I6,E 12
VSOP id. no. for
5))
the
Atom density for the
I-th
I-th
new
new
material.
material.
Card R19B only if HPOS > 0. on card R9.
Card R19B Format (1216)
HPOS
IBAE(I),1=1,HPOS
New materials only in the individual batches IBAE(I). (< 999)
92
Cards R20A and R20B only if IWATER > 0 on card R7.
Card R20A
1
2
3
4
5
PPC
EPSI
NRVO
NRVH
NZ
Format (2E12.5.3I6)
= 0.: No water ingress into reactor core (to reflectors only).> 0.: Water ingress into core;
PPC = Increase of partial pressure of steam in the core (bars).
Void fraction of the pebble bed.
VSOP-material no. of oxygen.
VSOP-material no. of hydrogen.
= 0: No water ingress into reflector batches.> 0: Number of reflector batches to undergo an ingress of water.
Read card(s) R20B.
Card(s) R20B only if NZ > 0 on card R20A.1 card for each of the NZ reflector batches.
Card R20B
1
2
3
4
NR
PPR
EPSI
T
Format (I12.3E12.5)
Id.-no. of the respective reflector batch.
Increase of partial pressure of steam in the respective batch (bars).
Void fraction of the batch.
Steam temperature in this batch.
93
When NKEEP = 1 on card R7 one full set of cards R21 - R23 for each of the different reloadbatches is required. This defines the FM-scheme. When NKEEP = 3, 4, 5 or 6, these cards areonly required for batches with important instructions. Batches which are only shuffled to the nextregion (NKEEP = 3) or to the next batch within the same region (NKEEP = 6) or stay in theirposition (NKEEP = 2 or 4) do not need the card R21.
FM means "Fuel Management".TBP means "This Batch Position".OPB means "Out of Pile Box".RPM means "Reprocessing Mixture".
Card R21
1
2
3
4
5
6
1X1
NRX
NSB
NREP
NSPALT
MAKEUP
Format (2I5,I2,6I4,3E12.5)
= 0: If a card R21 is included for each batch.> 0: When NKEEP = 3,4, 5 or 6: No. of batch position to which this
card R21 (and R22, R23) refers.< 0: Last card R21 holding for the batch |K1 | .
= 0: Refueling of this batch position TBP with fresh fuel specified by
1 / .
> 0: Id. no. of a batch which is shuffled into TBP.> 10000: Load storage box no. (NRX - 10000) into TBP. (This sto-
rage box must have been filled up at a previous reload!).< 0: New atom densities are loaded into TBP. A set of |NRX| den-
sities are defined on cards R22. The fuel type identification is un-changed (not a recommended option).
= 0: Normal.> 0: No. of storage box into which this batch is to be filled. Data can
be retrieved in the following reload. (See chapter A.2.2)
= 0: No reprocessing.= 1: Reprocessing before loading into TBP. This option only when
NRSTRT = 3 or 4 on card V6 and after having supplied cards Rl.
Number of materials for enrichment or re-enrichment. A card R23must follow.
If NREP > 0 and/or NSPALT > 0:VSOP id. no. of the isotope used as make up material in reprocessedand/or re-enriched fuel. The heavy metal density of the new batch isadjusted to the initial value of the loaded fuel type.If MAKEUP = 0:No material is added. Thus a new heavy metal loading is defined.
94
Continuation of card R21
7 MANAGE = 0: Load fuel without any change into TBP.= 1: Treated is the content of a reprocessing mixture which is the
considered OPB. It can optionally be loaded from a jumble box(comp. card R5 and NCYC on card R7). The content has beenformed from the disloaded fuel of one or more fuel types(XMARX on card R3) and has been summed up over previouscycles.
The following part of card R21 depends on the option MANAGE.
MANAGE = 0:
8
9
10
11
12
17
18
Rl
R2
R3
= 0: Fuel type no. same as batch NRX.> 0: Fuel type to be loaded into TBP (only with out of pile FM,
KUGL > 0 on card R2).
Dummy.
= 0.: Fissile enrichment stays unchanged.> 0.: Rl defines a new enrichment for the loaded fuel. Only relevant
if fresh fuel is used.
When NRX > 0 and use of in-core-batches:= 0.: New volume of TBP is defined by the batch NRX, which is
loaded.> 0.: R2 is the fraction of the in-core batch to be loaded into TBP.
When NRX > 0 and use of storage boxes:= 0.: The total storage box volume is used for loading into TBP. By
this way a new volume is assigned to TBP. A maximum is givenby filling up the region's volume.
> 0.: R2 is the fraction of the storage box to be loaded into TBP.When NRX = 0:> 0.: Fraction of the total volume of the out of pile fuel of type 17 to be
loaded into TBP. If KUGL = 1 (on card R2) the volume of allbatch positions in the upper region is automatically limited to theregion's volume. The fresh fuel store is unlimited.
= 0.
95
MANAGE = 1:
8
9
10
11
12
Card(s)
17
18
Rl
R2
R3
New fuel elements are performed in the following way:The identification number and the total heavy metal content are takenfrom the fresh fuel type numbered by 17. The isotopic composition ofthe new fuel elements is taken over from the reprocessing mixturenumbered by 18.
Id. no. of the used reprocessing mixture.
Enrichment Nflss/NHM for the new formed elements.
> 0.: Fraction of the total reprocessing mixture (RPM) volume to betreated and loaded into the volume of TBP.
< 0.: The RPM volume fraction |R2| is related to that part of the totalvolume, which has been left over from preceding loading proce-dures during the present fuel management step. If depletionwould be necessary, the program reduces the RPM volumefraction R2 instead.
= 0.: New volume of TBP is that one which has been made availablefrom the RPM.
> 0.: New volume of TBP is R3 * volume of the upper region whichthe presently considered batch belongs to.
= 1.: The upper region is filled up. The new definition of the TBPvolume immediately causes a corresponding change in the usedRPM volume fraction R2.
< 0.: The specified fraction of the RPM volume is reprocessed. Nofissile material is added or removed. Only make up material isadded or removed in order to achieve the required enrichmentRl for the defined fuel type 17. The presently considered TBPvolume is modified. For the presently considered TBP a volumeWERA = |R3| * volume of the upper region is made available. Ifthe prepared new volume of fuel elements is larger than WERA,the fraction R2 will be reduced. If it is smaller, the WERA willbe reduced correspondingly.
R22 only if NRX < 0 on card R21.
Card R22 Format (4(I6,E 12.5))
1,3,5
2,4,6
NPX(J),
CPX(J),J=1,|NRX|
VSOP id. no. of the J-th nuclide with the atom density *• 0.
Atom density of the J-th nuclide. If = 0., the nuclide needs not to bespecified.
96
Card(s) R23 only if NSPALT > 0 on card R21.
Card R23 Format (4(I6,E 12.5))
1,3,5
2,4,6
IDFISS(J),
FICOMP(J),J=l, NSPALT
KTYPE
HMETAV
VSOP id. no. of the J-th nuclide used for re-enrichment.
Relative fraction of the J-th nuclide in the enrichment composition.The sum of all FICOMP in the composition must be equal 1. Only thefissile isotopes in FICOMP are used to calculate enrichments, so thecomposition may also contain fertile materials, for instance 0.93 Uand 0.07 238U. The original fissile / HM ratio in the batch before re-enriching be YSPALT, the code distinguishes two different cases:
Case 1:
Rl > YSPALT, new material with the relative composition specifiedin FICOMP is added to make up the difference (Rl - YSPALT).
Case 2:
Rl < YSPALT, the original HM composition in the batch is un-changed, but the densities of all HM are reduced to obtain thefissile/HM ratio Rl. The out of pile volume fraction R2 (card R21) isreduced correspondingly. To maintain the correct HM density thedesignated make up material (normally a fertile isotope) is added.Here, the FICOMP data are obsolete.
= 0: Normal.> 0: The fuel in this batch position is given a new fuel type no. after
reprocessing and/or re-enrichment. Redefinition of types may benecessary in order to use pertinent cost data for recycled fuel.
= 0.: Normal.> 0.: New heavy metal density for use in this batch position, only of
significance in connection with the MAKEUP option (on cardR21). For some types of reactors the heavy metal loading in aparticular batch position may have to be varied during the life-time of the reactor, for instance during the running-in phase.
97
Cards R24 - R27 only if CONPOI > 0. on card R9.One card R24 for each of the CONPOI regions with control poison data.
Card R24
1
2
3
KR
POISL(l)
POISL(2)
Format (I12,2E12.5)
Id. number of the considered region.
Maximum atom density of the first control poison nuclide.
Maximum atom density of the second control poison nuclide (if de-fined).
Cards R25 - R27 like cards V12 - V14 !
4.4.8.5 Criticality search for the reloads. R28-R31
Cards R28 - R31 only if NRSTRT = 2 or 4 on card V6.
Card R28
1
2
JARIT
NCOL(I),1=1,JARTT
ITVAR
IR16
Format (1814)
= 0: No iteration for this reload, skip cards R29 - R31.> 0: Total number of batches to be iterated, length of following list of
batch id. no's. The atom densities of the materials specified oncard R30 are iterated to give reactivity K-search.
Id. no. of the I. batch.
= 0: No effect.> 0: Use different sets of materials to increase resp. decrease enrich-
ment to obtain correct Kefr. Two sets of cards R29 - R30:First set in case K-search > Kefr core, second set in caseK-search < Ke» core. Code reads both sets of cards and selectsthe required one in each case.
= 0: Normal.> 0: Read card R31 with new heavy metal densities for iteration
batches. IR16 is the number of batches in which the HM densityis redefined and used to determine the amount of make upmaterial to be added. The option may be necessary for reactorswhere the moderation ratio and HM loading in the fuel types varyduring reactor life, for instance during the running-in phase.
If J ARIT > 0 on card R28: At least one set of cards R29 - R30.
98
If ITVAR > 0 on card R28: Two sets of cards R29 - R30.
Card R29
1
2
3
4
ITMAT
XKEFF
MAKEUP
XTYPE
Card R30
1
2
3
4
IDIT(l)
COMPIT(l)
IDIT(2)
COMPIT(2)
Format (2(I6,E12.5))
Total number of materials iterated, i.e. length of materials list on cardR30. (<28)
> 0.: k(o), reactivity specification for iteration search in cycle i.= 0.: Use same k(o) as beginning of last cycle, k'(o) = k'"'(o).< 0.: Determine a k(o) value for beginning of next cycle so that the
end of cycle reactivity k(min) is reached after the specifiednumber of time steps JNSTOP. The extrapolation is made fromkj(o) = (k'-'(o) - k'-'CJNSTOP)) + k(min) * |XKEFF| and k(min)= ZKFTND (card V10). The value of XKEFF may be used toadjust for uncertainties in k(min). This is black magic.
VSOP id. no. of nuclide to adjust heavy metal density in batches toeither initial value in batch or as specified on card R31.
= 0.: Fuel type no. of batches is not altered.> 0.: New fuel type no. for batches. Same no. is given to all batches
for which iteration is performed.
Format (4(I6,E12.5))
VSOP id. no. of first nuclide used in iteration.
Relative fraction of first nuclide.
VSOP id. no. of second nuclide.
Relative fraction of second nuclide.
The relative fractions of all ITMAT (card R29) nuclides must be equal1. If all COMPIT = 0., the existing relative fractions in the batchesremain unaltered during iteration.
Card R31 only if IR16 > 0 on card R28.
CardR31
1 IR
Continuation
Format (I6.E12.5)
Batch no. for which the followingfication on this card only for thosediffers from the initial one.
ofcardR31
HM densitybatches for
is specified. The speci-which the HM density
99
HMETAV(IR)
New heavy metal density in batch no. IR, to be used when adjustingthe make up material.
One card for each specified batch: I=1,IR16.
4.4.8.6 Redefinition of CITATION edit options. R32
Card R32 only if IVSP( 11) = 2 on card R7.
Card R32 Format (713)
Redefinition of edit options. Same as card C4 (see input section4.4.6.2).
4.4.8.7 Extracted nuclides for printout. R33
Card R33 only if MUHU(l) > 0 on card R7.
Card R33 Format (1216)
NUPRI(I),1=1,
MUHU(l)
VSOP id. no. of nuclide for which printout is desired.Up to 20 nuclides can be specified.
4.4.8.8 "Status of core"- data set for TINTE. R34
Card R34 only if LIB > 0 on card R7.
Card R34 Format (18A4)
18
TrrEL(i),1=1,18
Literal description of the case to be transferred to the TINTE code1291 via data set 'tinte'.
100
4.4.9 Fuel power histogram for decay power evaluation. - LF1 - LF3
Only if IK = 3 on any card R7.
Card LF1 sets up the dimensions.
Card LF1
1
2
3
KM AX
LMAX
MTMAX
Format (316)
Number of VSOP burnup cycles required for the set-up of the fullirradiation history of the individual batches (compare KT5 > 0 oncard LF2).
Number of all VSOP time steps of the KMAX burnup cycles.
Number of graded time steps to be generated (< 49).
Card LF2
1
2
3
4
5
IOUT
ICOMPA
ICOMPE
LTO
KT5
Format (516)
Output option:= 0: Short output.= 1: Recommendable.= 2: Additional test output.
= 0: The generated library of all batches in the graded time steps isprinted out from batch no. 1.
> 0: Print out starts from batch no. ICOMPA.
= 0: Print out ends at the last batch.> 0: Print out ends at batch no. ICOMPE.
= 0:Default. Evaluation starts from the last time step of the givencycle.
> 0: Time step of the given cycle, from which the precursory historyevaluation begins.
= 0: No effect.= 1: Preserve only the last VSOP cycle and prepare KMAX identical
cycles out of it.> 1: Preserve only VSOP cycle with id. no. KT5 and prepare KMAX
identical cycles out of it.
101
Card LF3
1
2
3
4
TN
DT
TOOP
TEPS
Format (4E12.5)
Time span (days) to be covered by the coarse new intervals (> maxi-mum fuel element residence time + out of pile times).
First coarse time interval (days). Normally the same as the last VSOPtime step.
MEDUL-fueling only:Time span (days) between fuel discharge from the core and its reloadonto the core. Zero power is assumed during this shuffling procedure.
> 0.: Convergence limit for iterative calculation of the incrementalparameter of the coarse time steps.
= 0.: Default value = 0.1 .
4.4.10 Fuel irradiation histogram for entire isotope generation. - P
Only if LK = 3 on any card R7.
CardP
1
2
3
LXS
LTO
1OUT
Format (316)
Number of time steps with spectrum calculation (for dimensioningonly).
= 0: Normal. Evaluation starts from the last time step of the givencycle.
> 0: Time step of the given cycle, from which the precursory historyevaluation begins.
Output option:= 0: Recommendable.= 1: Additional test output.
102
4.4.11 Preparing THERMOS-library. TTTT1 - TTTT5
Only if ITTT > 0 on card SI.
Sequence of input cards:
SI - G6: VSOP input of case identification.
TTTT1, TTTT2: THERMALIZATION input without selfshielding factors.
T l : Blank card.
TTTT3 - TTTT5: Preparing (condensing) new THERMOS library.
CardTI
1
2
3
4
TT1
JNTAPE
NKER
IDKER(1,I),
IDKER(2,I),1=1,NKER
Format (1216)
Id. no. of the thermal 96 groups THERMALIZATION-library on dataset 'thermal'. Only choice at present: JNTAPE = 115 .
Number of different scattering nuclides for the present spectrum run.
VSOP-id.no. of the first scatterer.
Thermal library-id.no. of scattering matrix to be applied.
NKER different pairs of id. numbers and scattering matrices.
CardTI
1
2
3
TT2
TOM
EPSI
WAT
Format (3E12.5)
Temperature in calculation of Maxwellian neutron energy distributionfor starting iterations of thermal spectrum. (°K)
Criterion of convergence of flux iteration. (= 0.0001)
Acceleration of convergence: 0. = no, 1. = yes.
103
Card 1T1T3
1
2
3
4
5
6
7
DSIDTP
l l l l l
ffiBE
ITUTEU
KERNE
ITOT
T
Format (A9,I3,4I6,E12.5)
Data set name of the THERMOS library to be replaced or to be gene-rated in addition to existing libraries. For possible choices see Table I.
> 0: Reduced output.
Nuclide no. with full output.
> 0: Print out of scattering matrices.
Number of scattering matrices to be condensed for the THERMOSlibrary.= 1000: Condensing of all scattering matrices.
> 0: All absorbers.
Temperature (°K) for Maxwellian flux for eventually condensing ab-sorbers.
Card TTTT4 only if KERNE * 1000 on card TTTT3.
Card I I I 14 Format (1216)
1
2
IDKER(IJ),
IDKER(2,J),J=l, KERNE
GAM-I-id.no. of the J-th scattering matrix.
VSOP-id.no. of the J-th scattering matrix.
Card 11 115 Format (1015)
301=1,30
Number of the lowest THERMALIZATION group within the I-thTHERMOS group to be formed.
104
4.4.12 2d-Thermal hydraulics. TX1 - TX26
Only if NTIK = 1 or 3 on card R7.
CardTXl Format (18A4)
18
TITLE(I),1=1,18
Literal description.
CardTX2
1
2
3
4
IFKON
IPRDSTT
IPUN
IFRSTA
Format (1814)
Steering the calculation:
= 0: No calculation of gas temperature and gas streaming.= - 1 : Coupling between the temperatures of gas and solid material by
heat transfer coefficient a. Recommended for steady state calcu-lations, not valid in transient runs.
= 1: Coupling via the source/sink distribution.
= 0: Minimum output.(Recommended when temperature calculation isfrequently repeated during the VSOP-run).
= 1: Standard output. See also item "IPASS" (this input card).= 2: More detailed output.= 3: Maximum output (very much!).
= 0: No effect.Steady state calculation only:
= 3: Preserve temperature fields on data set 'thermix' as start-upvalues for a transient calculation.
= 0: No effect.= 3: Read temperature fields from data set 'thermix' as start-up
values for a transient calculation.
105
Continuation of card TX2
5
6
7
8
9
10
11
12
13
17
18
INTVAL
IFRED
MITMAX
IKORM
IFREL
ITLAM
NLOOP
EXPR
ICODEF(I),1=1,5
IPASS
= 0: Steady state run.= 1: Transient calculation: The time steps are given by the burnup
scheme (JNSTOP, DELDAY on card R9).
= 0: For steady state calculation.> 0: Calculation of the decay power according to DIN 25485, using
the fuel life history (see cards LF).
Iterations:
> 0: Maximum number of iterations of temperature calculation.= 0: Default value = 2000
> 0: Maximum number of changes of the relaxation factor.= 0: Default value = 100
= 0: Inner iteration in radial direction (I).= 1: Inner iteration in axial direction (N).
> 0: Repeat calculation of temperature dependent material data forevery ITLAM-th time step (only for steady state THERMEX-KONVEK iteration).
= 0: Default value = 10
> 0: Maximum number of THERMIX-KONVEK ("Loop") iterations( steady state).
= 0: Default value = 100
= 0: No effect.= 1: Write steady state temperature field onto data set 'tempstat'.= 2: Write temperature fields during transient calculation onto
data set 'tempinst'. Read card TX3.
= 0: Must be = 0 in steadv state calculations.> 0: Id. numbers of one or more (< 5) THERMIX-compositions,
whose properties shall be redefined for the transient calculation.Read a new set of cards TX8 and TX9 for each of thesecompositions, following card TX6. (Not allowed for thecomposition representing the heat (gas-) sink!).
Has effect only if item "IPRINT" (this input card) equals 1 and only forthe steady-state calculation.
= 0: Printout of fuel element temperatures for each batch, i.e. for eachfuel passage of a multi-pass fuelling strategy.
> 0: Printout for batch (passage) no. "IPASS" only.
106
Card TX3 only if EXPR = 2 on card TX2.
CardTX3 Format (2E12.5)
PHIA
PHE
Start of temperature table at axial position PHIA (lower value accord-ing to THERMDC mesh point positions), (cm)
End of temperature table at axial position PHE. (cm)
Card TX4
1
2
3
4
5
6
7
8
9
QNORM
ETHA
OVREL
ORMIN
TDIFF
EFAK
DTVOR
ZEITMI
EPSST
Format (10F6.1,E12.5)
> 0.: Total power (MW). Input power field is normalized to QNORM.In a transient run the QNORM must be the reactor power, forwhich the life histogram was calculated.
= 0.: Drop normalization.
> 0.: Convergence criterion for local THERMK temperature field.
= 0.: Default value = 0.01
> 0.: Maximum relaxation factor.= 0.: Default value = 1.7
> 0.: Minimum relaxation factor.= 0.: Default value = 0.6
> 0.: Relative convergence criterion of the time independent THER-MIX-KONVEK iteration.
= 0.: Default value = 0.0005
> 0.: Multiplication factor for maximum allowable error level, whichstops the run.
= 0.: Default value = 1.
> 0.: Maximum allowed relative temperature change AT / T in a timeinterval At of a transient run. The time intervals At are corre-spondingly adapted.
= 0.: Default value = 0.05
> 0.: Minimum length of the time intervals At in a transient run. (sec)= 0.: Default value = 60.
= 0.: Emission coefficient of graphite spheres from internal functionGREPS(T).
> 0.: New emission coefficient (e.g. for coated spheres).
107
Continuation of card
10
11
ZEITNW
FDOSE
TX4
= 0.: (Special option).
Only in case of steady-state calculation (INTVALFactor to convert the dose of fast neutrons with anE > CEG(l) (see card G6) into EDN-values.
= 0 on cardenergy
TX2):
CardTX5
1
2
3
4
5
6
7
8
9
10
11
EPSI1
EPSI2
OVM1
EPSI4
CP
PRAN
DRUCK
IFZDR
ITM1
ITM2
ITM3
Format (7E8.0,4I4)
> 0.: Relative criterion of convergence for gas temperature.= 0.: Default value = l.E-5
> 0.: Criterion of convergence for mass flow.= 0.: Default value = 0.01
> 0.: Extrapolation factor for iterations on mass flow (every 10iterations an extrapolation is provided with 1 + OVM1).
= 0.: Default value = 0.5
> 0.: Relative criterion of convergence of the avg. gas temperature inthe outer iterations between gas temperature and mass flow.
= 0.: Default value = 0.02
> 0.: Specific heat capacity of the gas. (J/kg/°K)= 0.: Default value = 5195. (He).
> 0.: Prandtl-constant of the gas.= 0.: Default value = 0.66
Pressure of the gas. (bar)
= 0: Pressure of the system is constant.= 2: Pressure changes according to temperature. Gas inventory is
constant.
> 0: Maximum number of iterations for gas temperature.= 0: Default value = 100
> 0: Maximum number of iterations for mass flow.= 0: Default value = 500
> 0: Maximum number of outer iterations between gas temperatureand mass flow.
= 0: Default value = 5
108
Card TX6 only if INTVAL = 1 (transient calculation) on card TX2.
Card TX6
1
2
3
4
5
6
DZEIT1
NPRIN1
NKONV1
ZEU
DZEIT2
PSPALT
Format (F6.1,2I2,3E 10.3))
Length of the first small time interval, (sec)
> 0: Print the fields of temperature and streaming for all NPRIN1small time steps.
= 0: Default value = 50 .
> 0: Run the calculation of gas temperature and gas streaming everyNKONV1 small time steps (only if IFKON * 0 on card TX2).
= 0: Default value = 1
End of the first large time interval, (hours)
The following large time intervals are defined by variable DELDAY(card RIO). The small intervals are defined by DZEIT2.
= 0.: Free choice of the small intervals.< 0.: Also free choice, but maximum = |DZEIT2|. (sec)
= 0.: No effect> 0.: Redefinition of the pressure in the gap compositions of
THERMIX
Cards TX7 - TX19 only if INTVAL = 0 (steady state calculation) on card TX2.
Card TX7 Format (414)
IFRFI
IFRFA
IFRFL
IFRFR
= 0: Default.= 1: Adiabatic boundary condition in the first radial mesh.
= 0: Default.= 1: Adiabatic boundary condition in the last radial mesh.
= 0: Default.= 1: Adiabatic boundary condition in the first axial mesh.
= 0: Default.= 1: Adiabatic boundary condition in the last axial mesh.
109
One card TX8 (optionally followed by TX9 - TX12) for each THERMIX composition num-bered continuously increasing.
Card TX8 Format (A3,6I3,1OE5.O,I1)
2
3
BEM
Kl
IFTV
IFWKT
IFLT
Literal description of this composition, 3 digits. Free choice, exceptfor:= HET: Temperatures are calculated in the inner of the fuel elements.
Possible only for a pebble-bed core!Analysis of temperature/volume.
= DBH: In a transient calculation the average and the maximum tem-perature of this composition is displayed as a function of time.
= END: Termination of cards TX8, drop all other items of this card.Code sets variable KMAX = number of THERMIX compo-sitions equal to highest value, which was assigned to variableKl.
Id. no. of this composition.
= -1: "Solid material zone". Temperature calculation comprises theheat exchange with the coolant by source/sink heat transferwithin the meshes of this zone.
= 0: "Solid material zone". No heat exchange with the coolant is in-volved.
= 1: "Fluid zone". No temperature calculation is performed for thiszone (Temperature fixed to initial value TVOR). Coupled withthe neighbors by the heat transfer coefficient ALP on this card.Note: For instance these zones are used as a heat sink at the outerboundary of the system. At least two meshes are required.
= 0: Heat capacity given by C on this card.> 0: Identification no. of the material function for temperature depend-
ent heat capacity (see Tab. VI).
= 0: Thermal conductivity A. given by LAM on this card.> 0: Identification no. of the material function for temperature and
dose dependent A, (see Tab. VII).= 7: The temperature dependent function of id. no. = 7 uses LAM0 of
this card as A.(T = 0°C).= 4: In case of EPS1 and EPS2 > 0. (see below) the function uses
LAM0 of this card as pressure (bar) of the gas in the gap(convection). In case of EPS 1 and EPS2 = 0. the function useshelium at the pressure 1 bar. No heat radiation.
110
Continuation of card TX8
6
7
8
9
10
11
12
13
14
15
16
IDIR
NTVAR
RHO
C
LAM
LAMO
EPS1
EPS2
R1R2
TVOR
WPR
Only if EPS1 (and EPS2) > 0.:= 0: Radiation in radial direction.= 1: Radiation in axial direction.
= 0: No effect.> 0: In case of fluid zone (IFTV = 1) provide time dependent tem-
peratures on card TX12.
Volumetric fraction of solid material in this composition. RHO is usedfor calculation of the heat capacity.
= 0.: When IFWKT > 0.> 0.: Heat capacity of the solid material. (J/cm3/°K)
= 0.: When IFLT > 0.> 0.: Thermal conductivity in solid material zones (only if IFTV = 0 or
-1). (W/cm/°K)
= 0.: Default.> 0.: If IFLT = 7, LAMO is MJ = 0°C).
If IFLT = 32, the conductivity according to function 32 ismodified by the factor LAMO.IF IFLT = 4 and EPS 1 > 0. and EPS2 > 0., LAMO is thepressure of the gas in this composition.
= 0.: No heat radiation.> 0.: Coefficient of emission for heat radiation: inner radial/ upper
axial wall of the gap.(Maximum number of compositions with heat radiation = 19).
Like EPS1. for outer / lower wall.
Horizontal gap: = 0.Radial gap: Ratio: inner radius / outer radius.
= 0.: Start-up temperature field results from the input temperaturefield of the cards TX15 - TX17.
> 0.: Start-up temperature of this composition (°C) superior to thestart-up temperatures of the cards TX15 - TX17.
= -1. : Field of fission power density results from CITATION.It will be normalized to QNORM (card TX4).
> 0.: Power density of this composition. (W/cm )
I I
Continuation of card TX8
17
18
ALP
ISTR
= 0.: No effect.> 0.: Heat transfer coefficient in fluid zones (only if IFTV = 1).
(W/cm2/°K)Note: ALP = 0.5: Temperature at the boundary close to that ofthis fluid zone. ALP = 0.01: Temperature at the boundary closeto that of the adjacent zone.
= 0: No gas streaming through this composition.= 1: Gas streaming; read card TX9 following this card TX8.
Card TX9 only if ISTR = 1 on card TX8.
Card TX9 Format (2I6,6E6.0,I6)
1
2
3
4
5
6
IFBQ
IFBR
PVOR
XKON
ALPHA
DHYD
= 0: When IFTV = -1 on card TX8. Convective heat source is com-puted in the meshes of this composition.
= -1: No convective heat source evaluation (e.g. in voids).
Type of composition:= 1: Gas streaming in the pebble bed.= 2: Gas streaming in vertical pipes. (Total number of such
compositions < 30)= 5: Gas streaming in a horizontal void (no more than one mesh over
its thickness).(10 horizontal voids at maximum!)
> 0.: Pressure at beginning of iterations, (bar)= -1.: Pressure = pressure of the gas (see DRUCK on card TX5).
Additional pressure drop relative to computed pressure drop over thelength of the channel (only if IFBR = 2). (I/cm)
= 0.: Internal calculation of the coefficient of heat transition. In voids(IFBR = 5) use ALPHA = 0.
> 0., * 1: No internal calculation of ALPHA. Use the given value.(W/cm2 K)
= 1.: For the pebble bed.
Hydraulic diameter (cm). Only if IFBR > 2.
The following 2 items are relevant for a steadv-state calculation only:
Continuation of card TX9
112
7
8
9
STZUK
TFLVOR
IFZST1
Cards TX10,TX 11 only
Card TX10
1
2
HKUG
NHZON
1= 1,NHZON cards TX
CardTXll
1
2
3
4
DI(I)
NHMATl(I)
NHMAT2(I)
XFWQZ(l)
>0.: Source of mass flow, (kg/s)< 0.: Sink of heat (mass flow). Mass flow according to
conservation law
Temperature of inlet gas. (°C)
In case of time dependent calculation only:
= 0: No forced cooling= 1: Time dependent forced mass flow. Read input cards TX25. TX26.
if BEM = "HET" on card TX8.
Format (E1O.3J5)
Diameter of the spherical fuel element, (cm)
Number of radial mesh intervals in the sphere. (< 5)
11.
Format (E10.3,2I5,E10.3)
Inner diameter of the I-th radial mesh interval, (cm)Caution: 1 = 1 . . . counts from the outer shell towards the inner.
Id. no. of temperature dependent thermal conductivity (see IFLT oncard TX8).
Id. no. of temperature dependent heat capacity (see IFWKT on cardTX8).
Shielding factor of the power density in the I-th shell (normally = 1. inthe fuel shells, = 0. in the outer graphite zone).
Card TX12 only if NTVAR > 0 on card TX8.
Card TX12
1
2
TKV(I),
ZEIV(I),1=1,NTVAR
Format (14F5.2)
Temperature. (°C)
Time, (h)
One card TX13 is required for each axial coarse mesh "N" (described on cards BI4).
Card TX 13 Format (2413)
113
1 KOC(1,N)
KOC(I,N)
Id. no. of THERMIX composition in the 1. radial coarse mesh.
Id. no. of THERMIX composition in the I-th radial coarse mesh.
(Radial coarse meshes are described on cards BI3)
Card TX14 Format (2413)
IYEAR(I),1=1, KM AX
Operating time (years) assumed for the calculation of the fast neutrondose values of the materials in the KMAX different THERMDC-compositions. These are required and applied only, if material
functionno. 32 according to Tab. VII is used for the calculation of the thermalconductivity of the respective material.
Card TX 15
1
2
3
IPOLI
IE
RE(I),1=1,IE
Format (2I5.6E10.3 / 7E10.3)
= 0: Linear interpolation of temperature input of cards TX17.= 1: Logarithmic interpolation (radial).
= 0: Drop cards TX16, TX17. Start-up temperature field read fromdataset 'tempstat'.
> 0: Number of radial mesh points for start-up temperature input.
Only if IE > 0:Radial mesh points for start-up temperature input on cards TX17.Continuation cards according to given FORMAT.
114
Cards TX16 - TX17 only if IE > 0 on card TX15.
CardTX16
1
2
3
IPOLN
NE
PHE(I),1=1,NE
Format (2I5,6E10.3 / 7E10.3)
= 0: Linear interpolation of temperature input of cards TX17.= 1: Logarithmic interpolation (axial).
Number of axial mesh points for start-up temperature input.
Axial mesh points for start-up temperature input on cards TX17.
Continuation cards according to given FORMAT.
One card TX17 for each of the N = 1,NE axial mesh points.
Card TX 17 Format (7E10.3)
1
IE1=1,IE
Start-up temperature at mesh point I, N.
Cards TX18 and TX19 only in case of a transient calculation (INTVAL = 1 on card TX2).
CardTX18
1
2
3
4
MC2
SIG
A0
SM
Format (I6,3E12.5)
= 0: No effect.> 0: Read card TX19 with definition of THERMIX compositions for
time dependent output of the "heat storage". (Also stored on dataset 'therlist').
= 0.: Default value = 1.> 0.: Factor to be multiplied with the explicitly evaluated decay
power function.
Initial fissile enrichment of the fuel. (w%)
Avg. heavy metal content per fuel element (incl. pure graphitespheres), (g/sphere)
115
Card TX19 only if MC2 > 0 on card TX18.
Card TX19 Format (7211)
KM AX
IKO(I),1=1, KM AX
"Heat storage" fraction id. no. to which the heat of THERMIX com-position I shall be added up. Possible "heat storage" fraction numbers:1 - 5 (core data is always stored). IKO(I) = 0 means no storage.
Card TX20
1
2
3
4
5
DELTAT
TU
TO
WRIT
RO
Format (5E12.5)
Desired temperature interval (AT) for the numerical integration insidethe fuel elements (°K). Up to 200 intervals are possible between TUand TO.
Lowest surface temperature of fuel elements.
Highest temperature at center of the fuel elements.
= 0.: Program uses standard data of the thermal conductivity as afunction of fast neutron dose and temperature.
> 0.: Like = 0, but various test output of temperature integration insidethe fuel elements in addition.
< 0.: Thermal conductivity as a function of fast neutron dose and tem-perature will be given on the cards TX21 - TX24.
= 0.: No effect.> 0.: Inner radius of the fuel matrix (if shell ball is considered).
Cards TX21 - TX24 only if WRIT < 0. on card TX20.
CardTX21
1
2
NSCH
KTEM(N),N=1,NSCH
LFAD(N),N=1,NSCH
Format (516)
Number of different functions of the thermal conductivity. (< 2)
Number of temperature mesh points for which the thermal conducti-vity will be given. (< 10)
Number of fast neutron dose mesh points for which the thermal con-ductivity will be given. (< 10)
116
For each of the NSCH thermal conductivity functions one set of cards TX22 - TX24.
Card TX22
1
KTEM
TSTUE(K),K=1,KTEM
Card TX23
1
LFAD
DSTUE(L),L=1,LFAD
Format (6E12.5)
Temperature mesh points.
Format (6E12.5)
Mesh points of fast neutron dose.
For each of the LFAD mesh points of the fast neutron dose one card TX24.
Card TX24
1
KTEM
WLSTUE(K),K=1,KTEM
Format (6E12.5)
Thermal conductivity at the temperature mesh points. (W/cm/°C)
117
Cards TX25, TX26 only if at least one of the IFZST1 = 1 on cards TX9. Up to 100 time steps canbe defined by cards TX26. Linear interpolation is provided between the time steps.
Card
1
2
TX25
IZK1
IZKOM(I),I=1,IZK1
1 card TX26 for each
Card
1
2
3
4
TX26
ZVOR
ZDR
ZST(I)
ZTF(I),I=1,IZK1
Format (4110)
Number of compositions with time dependent source of mass flowand / or inlet temperature. (< 3)
Id. no. of the I-th composition.
time step:
Format (8F9.3)
>0.:Time. (min)= 0.: End of the input of cards TX26.
Pressure, (bar)
Source of mass flow of the composition no. IZKOM(I). (kg/s)
Temperature of inlet gas of the composition no. IZKOM(I). (°C)
118
5. Input Manual V.S.O.P.-ZUT (log. unit 5)
5.1 Steering the execution mode. ZS
CardZS
1 MODE
Format (A8)
= zutalone: General input procedure for various kinds of fuel elementgeometry. Drop cards DZ. read complete set of cards Z.
= data2zut: More comfortable input procedure (making use of internal"auxiliary" subroutine ZDATA2):Read cards DZ and Z1-Z6.
This procedure is possible only for assemblies consistingexclusively of spherical fuel elements having coated
particlefuel, and graphite being the only moderator element
outsidethe coated particles.
5.2 Fuel element design. DZ1 - DZ9(only if MODE = 'data2zut' on card ZS)
One set for each variant of each desired fuel type (limited to 27 different sets).Calculation is terminated by one last card DZ1.
Card DZ1 Format (18A4)
1
18
TITLE(I),1=1,18
Literal description of fuel element-types and -variants.
TTTLE(l) = 'stop': This terminates the sequence of cards DZ1-DZ9.
Card DZ2
1
2
NFUTP
NFCP
Format (3I4,E12.5)
Identification of the fuel elements in 4 digits IJKL:U : Type (characterizes design), increasing numbers. (< 10)KL: Variant ( e.g. for different enrichments), increasing numbers,
starting from 01 for each type IJ.
Input option for the coated particle:= 2: Read cards DZ3 - DZ7.= 1: Read card DZ3 only.= 0: Data from preceding design.
119
Continuation of card DZ2
3
4
NFBZ
FF3
Input option for the fuel element:= 1: Read cards DZ8-DZ9.= 0: Data from preceding design.
> 0.: Volumetric filling fraction of spherical fuel elements in the core.= 0.: Default value = 0.61
Card DZ3 only if NFCP > 0 on card DZ2.
Card DZ3 Format (El2.5)
ANR Fissile enrichment of the fuel (fissile/heavy metal), atom fraction.= 0.: If INDBS (card DZ4) = 7 .
Cards DZ4 - DZ7 only if NFCP = 2 on card DZ2.
Card
1
2
3
4
5
DZ4
INDBS
NCT
NSIC1
NSIC2
IU8TH
Format (514)
Fuel identification:= 1: UO2
= 3: UC2
= 5: UC - ThC= 7: PuO2
= 9: PuO2 - UO2
= 2:UC= 4: UO2 - ThO2
= 6: UC2 - ThC2
= 8: PuO2 - ThO2
Total number of coating layers (< 5), to be numbered with increasingradius.
Number of the 1. SiC coating
Number of the 2. SiC coating
Preparation of ZUT-data for:= l:232Th= 2:238U= 3:242Pu= 4:235U= 5:24OPu
layer, if present.
layer, if present.
120
CardDZ5
1
2
3
4
RK
ROBR1
ROBR2
BETA
Format (4E12.5)
Radius of the coated particle kernels, (cm)
Density of the kernels, (g/cm3)
Density of 2. type of kernels, if present, (g/cm3)Only if INDBS = 4, 5, 6, 8, 9 on card DZ4.
Enrichment of uranium NU5/ Nv if INDBS = 4, 5, 6, 9 on card DZ4.
Card DZ6 only if INDBS = 7, 8 or 9 on card DZ4.
CardDZ6 Format (4E12.5)
1=1,4Atom fractions of the isotopic composition in plutonium:239Pu,24OPu,241Pu,242Pu.
Card DZ7 Format (6E12.5)
1,3,5
2,4,6
DCT(I),
ROCT(I),1=1,NCT
Thickness of the I-th coating layer, (cm)
Density of the I-th coating layer. (g/cm3)(Numbered with increasing radius, NCT on card DZ4).
Cards DZ8 - DZ9 only if NFBZ = 1 on card DZ2.Only a selected set of the following parameters of the cards DZ8 and DZ9 is required. Possiblecombinations are given in Table IX.
Card
1
2
3
DZ8
Rl
R2
FF1
Format (6E12.5)
Outer radius of fuel zone, (cm)(Fuel zone consists of coated particles
Outer radius of the sphere, (cm)
Volume fraction: coat.part. / (coat.part
and graphite matrix).
. + matrix)
121
Continuation of card DZ8
4
5
6
VMOD
INDBK
BK
Moderation ratio Nc/ NHM-
= 0: No "dummy" elements.= 1: "Dummy" elements existing.
Volume fraction: "dummy" elements / (fuel + "dummy") elements.
Table IX: Alternative specifications of spherical fuel elements
No.
INDBK
Rl
R2
FF1
VMOD
ROSM
BK
1
0
X
X
X
2
0
X
X
X
3
0
X
X
X
4
0
X
X
X
5
0
X
X
X
6
0
X
X
X
7
0
X
X
X
8
1
X
X
X
X
9
1
X
X
X
X
10
1
X
X
X
X
11
1
X
X
X
X
Card DZ9
1
2
3
4
5
ROSM
ROMTX
ROSCH
ROBK
SRO
Format (5E12.5)
Density of heavy metal, homogenized in the fuel zone, (g/cm3)
Density of graphite in the matrix, (g/cm3)
Density of graphite in the outer shell, (g/cm3)
Only relevant for INDBK = 1:> 0.: Density of graphite in the "dummy" elements, (g/cm3)= 0.: Density of graphite in the "dummy" elements equals ROSCH.
Inner radius of the matrix, (cm) (normally = 0.)> 0. for "shell ball" design.
5.3 Resonance integral calculation. Zl - Z17
122
Note:If MODE = ldata2zut' on card ZS and resonance parameters provided from a library, the input isreduced to the cards Zl - Z6.
In order to open a new direct access data set 'resint' for storage of the resonance integral valuesread cards Z17, Z6.
5.3.1 Short input Zl - Z6
Cards Zl - Z5 can be repeated for N different cases. The input stream is terminated by the cardZ6.
CardZl
1
2
3
4
IEND
JI
KI
RNRESO
Format (I1,I5,I1,E9.4)
9
00000
1 (Number of items on this card).
= 1.: Read resonance parameters for 235U from data set 'resdatu5'.= 2.: Read resonance parameters for 232Th from data set 'resdatth1.= 3.: Read resonance parameters for 238U from data set 'resdatu8\= 4.: Read resonance parameters for 242Pu from data set 'resdapu2'.= 5.: Read resonance parameters for 24OPu from data set 'resdapuO'.= 6.: Read resonance parameters from input cards Z7 - Z9.
From the value of RNRESO result different sequences of thefollowing input cards:
For RNRESO = I. - 5. the sequence is:a) MODE = 'data2zut': Zl - Z6b)MODE= 'zutalone': Zl, Z2, Z4, Z5, ZIO-Z16
For RNRESO = 6. the sequence is:Zl, Z2, Z4, Z7-Z8, Z5, Z9-Z16MODE must have the value 'zutalone'.
123
CardZ2
1
2
3
4
ETEMP
ENERGU
ENERGO
ESOLVE
Format (7E10.4)
Temperature of the resonance absorber. (°K)
Lower energy boundary for the set of resonance data to be respected.(eV)
Upper energy boundary for the set of resonance data to be respected.(eV)
Five digits IJKLM.O as specification of the calculation method:
I : Specifies the geometry of the absorber regions:= 0: Infinite size or numerical computation of the geometric escape
probability P(E). (Coated particle grain structure, see sectionA. 1.2)
= 1: Cylindrical geometry= 2: Slab.= 3: Spherical, analytic formula of P(E).
J : Specifies the method for scattering by the absorber:= 1: Down scattering based on the computed neutron flux.= 2: NR-approximation.= 3: IM-approximation.
K : In case of "double heterogeneous calculation" plus calculation ofunresolved resonances (IEND=6 on card Z8), the coating materialof the particles must be specified as moderator 1.= 0: Moderator no. 1 is not present.= 1: Down scattering in moderator no. 1 uses the computed
neutron flux.= 2: Down scattering assumes 1/E flux.
L = 0: Moderator no. 2 is not present (cp. item K).= 1: Down scattering in moderator no. 2 uses the computed neutron
flux.= 2: Down scattering assumes 1/E flux.
M = 0: Normal output option.> 0: More output.
124
Continuation of card Z2
5
6
7
TESTA
EW(13)
VAZVG
These three digits DK.O for output options. (Use 001.0)I > 1: Control data of the energy fine structure for each resonance.J > 1: Partial probabilities at the mesh points of the P(E) - calcu-
lation.K > 1: Mesh points of P(E), Dancoff factors and data of homo-
genization of the matrix zone.
Number of mesh points for P(E)-calculation. (= 20)
Volume of absorber / total cell volume.
Card Z3 only if MODE = ldata2zut' on card ZS.
CardZ3
1
2
3
4
5
6
FUTYP
EAMOD1
ESIGM1
EAMOD2
ESIGM2
ECDANC
Format (6E 12.5)
Fuel type and -variant specification to be applied. See variableNFUTP on card DZ2).
Atomic weight of moderator material 1, which is admixed with theabsorber, (only if K > 0 in ESOLVE on card Z2).
as of the moderator 1 (only if K > 0 in ESOLVE).
Atomic weight of moderator material 2, which is admixed with theabsorber (only if L > 0 in ESOLVE).
Gs of the moderator 2 (only if L > 0 in ESOLVE).
Dancoff factor.= 0.: For infinite absorber size or numerical computation of P(E) for
the spherical fuel elements.> 0.: Dancoff-Ginsburg factor. Required for the cylindrical fuel
element. This ECDANC has higher priority than the internallycalculated one. which can optionally be ordered bv the card Zl 1.
125
Card TA
1
2
3
11
IDSATZ
IDNUCL
LOESCH(J),J=l,9
Format (1116)
Identification number of this set of resonance integrals to be stored ondata set 'resint' for further use in GAM-I-calculations (VSOP-MS).
GAM-I-library identification no. of the absorber nuclide.
Id.-numbers of existing resonance integrals on data set 'resint' to bedeleted prior to creating the new set.
CardZ5 Format (4A3,24X,12A3)
1
4
5
16
HEAD(J),J=l,4
HEAD(J),J=5,16
Literal heading, e.g. date.
Literal heading, e.g. case identification.
Termination of the input sequence by card Z6.
CardZ6 Format (11,171)
1
2
END
Blank: Terminates the sequence of resonance integral calculations.
126
5.3.2 Resonance parameters. Z7-Z9
Cards Z7 - Z9 only if RNRESO = 6. on card Zl.
Note:The input of resonance parameters on cards Z7 - Z8 must be terminated by a blank-card.
One card Z7 for each resolved resonance. In rising sequence of the energy EZERO.
Card Z7
1
2
3
4
5
6
7
8
ffiND
JI
KI
EZERO
GAMN
GAGM
R
S
Format (I1,I5,I1,5E9.4)
2
00005
5
Energy at the center of the resonance. (eV)
rn (eV)
ry (eV)
= 0.: Mesh spacing under the resonance decided by the code.> 0.: Mesh spacing under the resonance.
Ratio: Range of integration / effective width. Give S = 5.
Card Z8 for the unresolved resonances.
CardZ8
1
2
3
4
5
END
JI
KI
EC
GAMNO
Format (11,15,11,6E9.4)
6
00000
6
Lower energy of the range of unresolved resonance evaluation. (eV)
Avg. [ rn° ] (eV)
127
Continuation of card Z8
6
7
8
9
GMGM
G
D
EO
Avg. [TY] (eV)
Statistical weight.
Avg. spacing between the resonances. (eV)
Upper energy of the range of unresolved resonance evaluation. (eV)
CardZ9
1
2
3
4
5
6
END
JI
KI
AZERO
G
SIGPZ
Format (11,15,11.3E9.4)
1
00001
3
Atomic weight of the absorber.
Statistical weight factor.
o s . Potential scattering cross section of the absorber.
128
5.3.3 Explicit fuel element design. Z10 - Z16
Cards Z1O - Z16 only if MODE = 'zutalone' on card ZS.
Card Z10 only if Dancoff factor is to be calculated.
Card Z10
1
2
3
4
5
6
7
8
9
10
END
K
KI
A
B
DHUE
DSP
SMOD
SHUE
SSP
Format (I1,I5,I1,7E9.4)
7
Type of lattice:00010: Square lattice up to the 4. neighbor.00014: 4x4bundle. }00015: 5 x 5 bundle. } Finite lattice00016: 6x6bundle. }00020: Triangular lattice up to the 4. neighbor.00022: 2 - rods bundle.00023: 3 - rods bundle.00027: 7 - rods bundle.00029: 19-rods bundle.
7
Radius of the rod. (cm)
Pitch. (Distance between the center lines of the rods), (cm)
Cladding thickness, (cm)
Thickness of the gaps between the bundles (for K = 14, 15, 16) (cm)
Xtoi of the moderator, (cm4)
Lot of the cladding, (cm1)
I,o, in the gaps, (cm1)
129
Card ZI 1
1
2
3
4
5
6
END
JI
KI
SOLVE
ABAR
C
Format (I1,I5,I1,3E9.4)
1
00010
3
Same value as ESOLVE on card Z2 must be given here.
Radius (cm), for cylindrical or spherical fuel. Half thickness (cm), forslab.= 0.: Infinite absorber size, or explicit calculation of P(aa).
Dancoff factor.= 0.: For infinite absorber size or numerical computation of P(E) for
the spherical fuel elements.> 0.: Dancoff-Ginsburg factor. Required for the cylindrical fuel
elements. This C has higher priority than the internally calculatedone, which can optionally be ordered by the card Z10.
Card Z12
1
2
3
4
5
6
7
8
9
10
END
JI
KI
EDZERO
EAMOD1
ESIGM1
EDIQU1
EAM0D2
ESIGM2
EDIQU2
Format (11,15,11.7E9.4)
1
00013
7
Nabs atom density of absorber (atoms / (barn*cm)).
Atomic weight of moderator no. 1.
<7S of moderator 1. (barn)
Nmod i / N a b s ratio.
Atomic weight of moderator no. 2.
c s of moderator 2. (barn)
Nmod 2 / N a b s ratio.
130
Cards ZI 3, ZI4 only if numerical evaluation of P(E) is required (I = 0 in ESOLVE on card Z2).
CardZ13
1
2
3
4
5
6
7
8
9
IEND
JI
KI
Rl
R2
R4
R5
F
H
Format (I1,I5,I1,6E9.4)
5
00018
6
Radius of the kernel of a coated particle, (cm)
> 0.: Outer radius of a coated particle, (cm)= 0.: Coated particles and matrix are treated homogenized.
Outer radius of the matrix, (cm)
> 0.: Outer radius of a spherical fuel element, (cm)= 0.: Cylindrical fuel element.
Volumetric filling fraction coat.part. / (coat.part. + matrix);(matrix including possible inner coolant / graphite zones)
= 0.: Default> 0.: Fraction of dummy graphite spheres in the pebble bed.
Card Z14
1
2
3
4
5
6
7
END
JI
KI
SI2
SI4
SI5
ALPH
Format (11,15,11,4E9.4)
5
00025
4
Avg. L,OI of the coatings, (cm1)
Xtot in outer shell of a spherical element
Xtot in dummy graphite spheres, (cm1)
Ratio of: (E,ot of the matrix / average Z,o
(cm1)
of the coatings).
131
Card Z15 Format (II)
ffiND 8: Termination of input of this resonance integral calculation.
Card Z16 Format (A3)
HEAD(l) REP: Another input case will follow (possible only for the sameresonance absorber).
END: Termination of the calculation.
5.3.4 Opening of a new resonance integral data set ('resint'). Z17
CardZ17 Format (1X.I5)
JI = - 30: Opening of direct access data set 'resint' for the resonanceintegrals in GAM-I-group-structure.
132
A. Appendix (Comments)
A.I Neutron spectrum calculation
A.1.1 Resonance integrals
For the isotopes 232Th, 238U, 240Pu and 242Pu resonance integrals are evaluated by the VSOP-
ZUT-section. The corresponding cross sections can be transferred to the GAM-I code. They
are added to the background absorption cross sections of the GAM-I-library. This background
is independent of lumping effects and temperature. (The GAM-I cross sections in the region of
the unresolved resonances are treated as infinitely diluted and are included in the background
data. They have to be corrected by means of selfshielding factors in case of high lumping of
the resonance absorber).
Prior to do a VSOP-MS task, the resonance absorption cross sections must be prepared for the
considered fuel assemblies for different absorber concentrations at different temperatures.
They are stored on the permanent data set 'resint'. VSOP-MS applies these sets for the
different spectrum calculations during the following of the reactor operation. The dependence
of the cross sections on the temperature and on the fuel burnup is achieved by linear
interpolation between the respective cross section sets.
A.1.2 Coated particle grain structure
At energy ranges with Xa (E) smaller than the mean cord length 1 of a coated particle, the grain
structure is of importance in spectrum calculations. This is true in the resonances and at the
lower end of the thermal spectrum. Therefore, the capability of grain structure effect has been
included in the ZUT and THERMOS codes.
The resonance integral calculation of the standard ZUT code is made for a homogeneous
distribution of the resonance absorber in the finite volume of a lump. The transport equation is
solved in very fine groups over the energy range of each resonance. The calculation also
includes the neutrons which are born in the lump, leave it, and are absorbed or scattered down
in any other lump of the same configuration. Excluded are those neutrons which leave the
lump and undergo scattering reactions outside. Nordheim /8/ excludes these neutrons by the
geometric escape probability
133
in which C is the Dancoff factor and Po (E) is the probability for a neutron to escape from thelump of its birth. This method has proven to be a good approximation for lumps of all degreesof grayness.
In case of coated particles inside a fuel element the absorber lump is the inner kernel of a
representative coated particle. Because of the smallness of a particle, the escape probability
Po(E) is close to 1, even for neutrons with energy close to the peak of strong resonances. The
neutron can travel through many coated particles without any collision, whether through the
coatings or stripping through the kernels. It can meet the boundary of the fuel matrix, pass
through the outer shell of the fuel element, enter another matrix, and undergo collision in any
of its coated particles or somewhere between them. Fig. 6 gives the different possibilities to
escape from a coated particle.
For such a "double heterogeneous" composition of the absorber lumps a Dancoff factor is hard
to define. Therefore in ZUT-DGL 191 the escape probability P(E) is directly evaluated by a
numerical method.
The possible path of a neutron is subdivided into parts for which the probability of traversing
can rigorously be evaluated by a numerical treatment. This requires the evaluation of 8
different probabilities Wi - W8 as indicated in Fig. 6. For instance, W| is the probability for a
neutron to undergo a collision in the coating of the coated particle in which it was borne.
Finally, the geometric escape probability is
P{ E) = W} + W2 (W3 + WA) + W2l-W.
This formula can replace the P(E) of Nordheim in the version ZUT-DGL. It has been derived
for spherical and for cylindrical elements. The outline of the numerical treatment is given in
Ref. 191.
Fig. 6: Break down of the neutron escape probability.
134
In the thermal energy range - i.e. in the THERMOS code- the grain structure is treated
analogously: in the fuel matrix the mean path of a neutron from one coated particle to the next
one be L. Its magnitude results from the diameter and from the volumetric filling of the
particles in the matrix. For a neutron at the energy E the probability W(E) of traversing one
coated particle and the corresponding amount of matrix material is calculated by a direct
numerical integration as outlined in Ref. /16/. Thereupon, W(E) is used to define an effective
macroscopic total cross section £e(E) by the equation
Ze(E) is used to replace the homogenized 1(E) of the mixture of matrix and coated particles in
the THERMOS calculation. The individual reaction rates of the different materials are defined
correspondingly.
A.1.3 Selfshielding factors in the epithermal energy range
The GAM-I code performs epithermal spectrum calculations for the nuclide compositions
resulting from a homogenization of the respective regions. Selfshielding factors can be
applied for the many nuclides. They can be given in coarse energy groups or in the 68 group
structure of the GAM-I library. Two types of selfshielding factors can be submitted:
SC: Cross section selfshielding factors allow modification of the microscopic cross sections
of the library. This might be desired to account for resonance shielding effects, for changes of
the neutron energy spectrum within a fine group g, or for improved measurements of cross
sections, respectively. The code allows input of LSUB different subsets of cross section-
selfshielding factors SCg . They can be applied for any nuclide in any region. The SCg are
multiplied to the respective cross sections ag of the energy group g.
SF: Neutron flux-selfshielding factors allow to bring the cell structure local neutron flux in
relation to the average cell flux.
Assuming a space and energy dependent cell calculation, the reaction rate per volume Vo of
the cell is given by
135
g=\ k=i
being:
g, NG energy groups of the GAM-I library
k, NK cell zones
Gg, SCg cross section and corresponding selfshielding factor
V* volume of cell zone k
fag neutron flux
(Note: R represents the reaction rate of any nuclide in any region. Respective subscripts have
been dropped.)
In view of the homogeneous treatment of the cell in the GAM-I code, the neutron reaction rate
per cell volume can be formed in terms of the average cell flux <(>og:
v o v I
with the following cell zone averaging applied
With the definition of the neutron selfshielding factors of the different cell zones
the reaction rate can be written in the form
where the cross section is modified by the selfshielding factors SC and SF
136
fiN V kg
k=l lyo VO
These modified cross sections a g are applied in the GAM-I spectrum calculation. They are
further applied to form the broad group cross sections which are turned to the diffusion and
burnup calculations.
The selfshielding factors SCg and SFtg are defined by the data input according to cards
G7 - G12. They may be given in few broad energy groups J. The code books them into the
respective fine groups g.
A.1.4 Leakage feedback
The 2 / 3-dimensional diffusion calculation by the CITATION module provides leakage terms
Lsi for the different regions S and coarse energy groups / (option IBUCK >0 on card V6).
They are available for the subsequent spectrum calculation. For the first one at the beginning
of the reactor operation all leakage terms have a zero value. Startup leakage terms can be
generated by running a dummy initial cycle of a very short time interval with practically zero-
burnup.
For the thermal cell code THERMOS the thermal leakage is transferred into the albedo at the
surface of the cell:
h
This is the ratio of the current density/ of neutrons entering the cell divided by j+ leaving thecell. The partial currents are given by
, _0o__A ; - 0 0 , ÄJ~ ~ 4 2 h " 4 2 '
j 0 = j + - j . being the net current leaving the cell per cm2. The net current of the cell is equal to
the ratio of the leakage Lc of the cell per surface Sc
-IOn the average, the leakage of one cell is equal to the leakage of the whole region divided by
the number of cells (nc) in region:
=^- and
137
being L :̂ Leakage of the cell
Ls: Leakage of the region
Vc: Volume of the cell
Vs: Volume of the region
e: Void fraction between the cells, if present
Approximating the neutron flux at the surface of each cell by the average neutron flux Os of
the region then leads to the albedo
<PS Vs(l-e) Sc
2 L V1 + —• ' '(f>5 V,-(l-e) Sc
138
A.2 Design specifications
A.2.1 Reactor layout (input cards BI and TR)
The geometric design of the reactor is provided in the subroutines BIRGIT (2-d) or TRIGIT
(3-d). In VSOP the basic unit of reactor material compositions is named a "batch". For the
first core loading the reactor design must be subdivided into batches. They all must be loaded
with fuel material or - outside the active core- with the materials of the reflectors etc. The
calculation is performed individually for each batch: For the in-core batches it is the follow of
the burnup, of the fuel shuffling, cost evaluation, and of the decay power production in
transient temperature calculations.
In many cases different types of fuel elements, or elements of a different irradiation age form a
nearly homogeneous mixture within the reactor core. They are exposed to the same local
neutron flux. For this purpose a mix of the respective batches can be put together forming a
"region". These regions represent volume parts V(I) of the reactor, which provide the
distribution of materials (and of their cross sections) for the calculation of the spatial neutron
flux distribution. Spectrum calculations are based on the averaged atom densities of the
regions. They provide the broad microscopic group cross sections.
In the pebble bed reactor the fuel elements move downward from the top towards the bottom
of the core. Here, the shape of the regions and their shuffling must fit into the flow pattern. On
the other hand the calculation of the neutron flux is performed by means of the code member
CITATION, and this is confined to a pattern of "CITATION-material compositions" W(J)
with perpendicular boundaries in e.g. r-z coordinates. Similarly the thermal hydraulics code
member THERMDC is subject to a mesh lattice of perpendicular r-z coordinates. Transfer of
the relevant data of the "regions" to CITATION material compositions and back is provided
by a volume matrix which is generated in subroutine BIRGIT.
The subroutine generates both reactor "regions" V(I) and CITATION material compositions
W(J). It then synthesizes a matrix of volumes VW(I,J), which is the overlapping set of the V(I)
and W(I) as shown in Fig. 7.
The transfer of the nuclear data proceeds as follows:
• Macroscopic cross sections Z are made for the "batches" and thereafter for the "regions":
The 1(1) are converted into macroscopic cross sections I(J) of the CITATION material
compositions by
139
Code section CITATION provides the criticality and neutron flux calculation.
Neutron fluxes O(J) of the CITATION material compositions are transformed to fluxes
of the "regions" by
and this is applied for the further burnup calculations of the "batches".
VSOP CITATION Overlay
t:::i
v-
t -
-i-.
• •
• - - •
1-I-4--
;
:
. . .
. .
. . .
-4
::
4
—
.. .
. . . .
'"1
i '•i
. . . .
+
• rL-
i-i-
F
_
Reqiom Compoaiiiona
W(J) VW(I.J)
Fig. 7: Overlay of reactor-regions and CITATION-material compositions
Analogously, a transfermatrix VW(I,K) between the core-"regions" V(I) and the corresponding
fine mesh volumes W(K) of the code member THERMIX is provided. Here, the fast neutron
dose, the local decay power and region identification numbers are transferred to the
THERMIX. The temperatures of the fuel and of the moderator of the different regions are
turned back to be used for the evaluation of the neutron spectrum.
140
As seen from Fig. 7, the two different mesh grids overlap in the area of the core. But they are
congruent in the reflectors.
In case that parallel movement of the fuel elements from the top to the bottom of the core may
be assumed as a sufficient approximation the data input is restricted to cards BI-1 through
BI-5, and thus is rather short and simple. A more sophisticated simulation of the movement of
the fuel elements can be defined by input according to cards BI-6 through BI-9 in addition.
Here, a more complex geometry of the core bottom of a pebble-bed core and very special flow
patterns of the fuel - resulting from experimental data or from theoretical research - can be
considered. The limiting curves of "flow channels" are defined by few coarse points, and the
curves are gained by interpolation. The radial position of the coarse points can internally be
modified in order to adjust the channel volume to a predefined value. Each channel is
subdivided into regions V(I), which are numbered by the code from top to bottom starting with
the first inner channel. Subsequent upon the highest region number within the core, the
reflector regions are numbered in the order as given on the input card BI-5.
The matrix of volumes VW(I,J) is derived in the following way: A very fine mesh grid of an
elementary volume of few cm3 is superposed over the given grids. For every small mesh the
code identifies the respective V(I) and W(J) in which it is located. It adds the elementary
volume to the corresponding element of the VW(I,J) volume matrix.
Similarly, subroutine TRIGIT prepares the geometric design in 3 dimensions. Here, "regions"
and CITATION material compositions must be identical.
A.2.2 Out-of-pile fuel positions
During the bumup cycle the data for the many batches of the reactor (core, reflectors etc.) are
stored in an array in the common. The atom densities within the core area are subject to
burnup.
Beyond the batches of the reactor itself, there are numerous further positions reserved, which
represent out-of-pile fuel positions. They serve as storage for fresh fuel, and for the disloaded
fuel which can be reloaded to the reactor, or stored, or reprocessed, or removed and sold,
respectively. Four different types of out-of-pile positions can be defined:
141
Intermediote boxes:
Storoge boxes:
Reprocessing mixture Aging boxes Jumble box
Fig. 8: Out-of-pile fuel positions
1. Fuel types. In the course of the shuffling the reactor batches can be loaded with fuel fromthese fuel positions. 1 through 10 different fuel types can be defined.
2. "Storage boxes". Shuffling specifications from cards R21 can also direct fuel into storage
boxes. Temporarily it is stored in intermediate boxes, because the content of the storage
boxes stems from the preceding shuffling and is available for reload at the present
shuffling procedure.
After the shuffling has been finished the remaining fuel of the storage boxes is turned to
the scrap fuel of the respective fuel types. Thereupon, the content of the intermediate boxes
is directed into the storage boxes, then being ready for use at the subsequent shuffling.
3. "Reprocessing mixtures". The scrap fuel of one or more fuel types can be directed into
reprocessing mixtures. Here, the amount of fuel adds up from cycle to cycle, and parts of it
can be used for recycling. Re-use of the fuel can be defined with reprocessing and with
re-fabrication of new fuel elements.
4. "Aging boxes". This is an extension of the reprocessing mixture option, which simulates an
intermediate storage for isotopic decay. This option allows real bookkeeping of the out-of-
pile inventories with given decay periods prior to reprocessing. It has been developed for
the simulation of closed fuel cycles.
The scrap fuel is first turned to a sequence of "aging boxes" in which the decay proceeds
over a given time period. During this period the fuel is bound to the out of pile storage and
is not available for re-use. The last box (named jumble box) is reserved for accumulation of
142
the fuel which has gone through the aging period. At every shuffling step its content is
turned to the respective mixture for the purpose of re-use. After the end of the shuffling
procedure the amount of unused fuel is turned back to the jumble box.
A.3 Simulation of reactor operation
As seen from the input manual, the design of reactor life time follow requires information
about many different physical events which are mutually coupled with each other. The user
should be familiar with the respective physical laws and with their representation by the
calculation model. Beyond the input description the following comments might give a further
help in setting up a reactor simulation.
For the repeated calculation of the thermal neutron spectrum an efficient acceleration is
included in the code: It preserves the field of the neutron flux and provides it as initial guess
for the subsequent calculation. The same is done for repeated diffusion calculation. If during
the burnup periods in between the isotopic concentrations change only slightly, the number of
iterations is reduced this way by about 90% or more.
A.4 Restart
When the option to prepare restart data is defined, the code prepares several data sets to be
used at the restart of the calculation (see chapter 3.4 for details). Using the restart option,
parametric research can easily be made at different time steps during the simulation of reactor
operation, i.e. at selected time steps of the running-in period or at different time steps of an
annual fueling cycle, e.g. for power transients, shut-down procedures, accident analysis or cost
studies.
143
A.5 Fuel cycle costs
Evaluation of the fuel cycle costs is based on the present worth method. It is performed by
means of the KPD code 1201, which is included as module KOSTPW. Before and after each
step of fuel shuffling, the code accepts the inventory of the heavy metal isotopes of all in-core
batches and of the out-of-pile fuel. All fuel which enters the system is regarded as bought, and
the fuel leaving the system is regarded as sold. Basing on this knowledge the code evaluates
expenditures and revenues individually for each of the examined cycles. From this result the
fuel cycle costs of each cycle dated to the start of the cycle.
Life time fuel cycle costs are derived from the individual cycle costs, which are re-dated to the
start up of the reactor operation. For this evaluation the last calculated cycle is considered as
an equilibrium cycle. It is considered to be identically be repeated until the end of reactor
operation.
An important role is played by the delay times of payments and revenues, especially for the
costs of fabrication and of spent fuel, respectively. The many choices of lead and lag times are
illustrated in Fig. 9 (left hand side). They are grouped together as shown on the right hand side
of the figure.
Break down of the fuel cycle costs is evaluated in four different terms:
1. Fuel costs include the expenditures for the fuel and working capital costs.
2. Revenue for the spent fuel assumes reprocessing. In case of final storage it can be
nullified by depreciation factors (input card K8).
3. Fabrication costs, which also include the re-fabrication in closed cycles
4. Reprocessing costs, which more accurately should be named "spent fuel handling costs",
because they could also be the costs of final disposal if an adequate value is given for
the input variable CAUF on card K7.
For the purpose of parametric studies, the cost input data may be varied in a series of code-
restarts.
144
Actual Cosh Flow during Simplified Cosh Flow in KPDLifetime of a Fuel Bolch
3.5'
Purchase
Shipping
U33g-U% Conversion
Shipping
Enrichment
Shipping
UF5-UO2 Conversion
Fabrication
Shipping
9-. Ao o6' «
26'
Cooling and Storage
Shipping
Reprocessing
Credit
Purchase
ConversionEnrichment
Fabrication
° ?9-. a2. 9,51 S
Reprocessing+ Shipping
Credit
Fig. 9: Lead and lag times of payments
145
A.6 Thermal hydraulicsA.6.1 Structure of the THERMIX code
The THERMEX-KONVEK code has been developed for the thermal hydraulics evaluation of
the pebble-bed HTR in two dimensions (r-z geometry) /19,22,23/. The code may simulate
steady state and transient conditions. It is linked into VSOP (see Fig.2 and Fig. 3). It receives
the power distribution of the reactor core from the nuclear code modules at given points in
time. It then returns the corresponding temperatures of the fuel and of the moderator averaged
over the volumes of the reactor regions to be used for further neutronics evaluation
(Fig. 10).
VSOP THERMIX
Neutron spectrum, -diffusionH K-eff (steady state)
5 Power density distribution
Burnup. 1-135. Xe-135,normalization of (lux
Spectrum, diffusion, K-eff
Control of reactor power
Burnup etc.
New time intervals
gD-Power density.
Temperature in fuelelements and Helium(transient)
2D-Decay power
Fig. 10: Coupling between neutronics and thermal hydraulics
146
Given distribution of temperatureand power over reactor and fuel
elements
Iter
atio
n
Iter
atio
n
"rat
ion
Iter
atio
n
f
Mass flow of the cooling gasPressure field
Temperature of the coolimg gasConvective sources and sinks
c r >»
Iter
atio
« • -
r
Heat transport in solid materialover reactor and fuel elements
when steady state
f
Temperatures of fuel and moderatoraveraged over VSOP-spectrum
zones
r
MAIN
STtUfcK
KONVEK
STROEM
L GASTEM
L TFELD
Fig. 11: Flow scheme of the THERMIX
147
A.6.2 Decay power for transient THERMIX calculations
In case of a shut down of the reactor under normal conditions or by a loss of cooling event the
power of the core is reduced to the decay power of the unstable isotopes. As the decay power
of each unit of fuel depends on its power histogram, it is not only a function of the decay time
but also of space and must be evaluated individually for each fuel batch within the reactor
core. For this purpose the code NAKURE /21/ has been included as subroutine NACHW into
the THERMIX code.
A.6.2.1 Power histogram of the fuel batches
The core is subdivided into regions, which again are composed of a mixture of one or more
batches (Section A.2.1). These batches go through burnup and fuel shuffling. For the analysis
of e.g. a loss-of-coolant incident the knowledge about the decay power of every batch is
required as a function of its preceding irradiation history.
In a VSOP calculation subroutine LIFE prepares a power histogram of the fuel batches, which
is the basis for the evaluation of the local decay power by means of the NAKURE code and by
subroutine NACHW of THERMIX, respectively. This histogram contains the relevant data of
the prehistory of every batch in coarser time intervals, i.e. the power, the burnup, the fractions
of fissions of 233U, 235U, 239Pu, 241Pu, and the rate of neutron capture by 232Th and 238U.
Most complex is the compiling of the irradiation history for the multiple passes of the
elements through the reactor. It is made by analogy to the fuel shuffling. In the calculations
the length of a shuffling cycle can vary, shutdown periods may be included, and the residence
time in the out-of-pile boxes normally is variable, too. In order to form averages of different
histories a standardized set of time intervals must be defined. The code transforms the history
data from the given pattern of time steps into the standardized intervals. This transformation is
made for the irradiation history of every batch.
The present version of the code allows to subdivide the whole burnup period into 49 intervals
of graded length. The first interval DT should be selected as short as the last cycle was prior to
reactor shut down at time to. The time steps are made from a geometrical progression for the
intervals I:
TG(I)=TG(I-J) + DT*(J+E1)'''
El is an incremental parameter derived iteratively by the condition
148
TG{MTMAX)-TN2 <TEPS
being
MTMAX: Given number of graded time steps (< 49)
77V: Space of time to be covered by the graded time steps
(> maximum fuel residence time)
TEPS: Criterion of convergence (e.g. = 0.1)
By this way intervals of high importance for the evaluated decay power are small and those of
minor importance become longer.
The code prepares the histogram according to any fuel management scheme applied in VSOP.
It can evaluate the running-in period of a reactor and any load follow prior to the shut-down
time to.
By option it is also possible to evaluate an artificial equilibrium cycle. In this case it prepares
an equilibrium operation by copying one designated cycle many times. Then it proceeds as
described above.
On the basis of this power histogram the values of the decay heat production of the fuel for
use in THERMIX are then calculated by means of the subroutine NACHW. It employs the
methods of the German Standard DIN 25485 /26/, which has been established for the
evaluation of the decay heat production of High-Temperature-Reactor fuel.
A.6.2.2 Decay power of the fission products
The contribution Ps of the fission products to the decay power - calculated as the sum of the
separate contributions of the considered fissile isotopes i - is given as:
T being a time period of reactor operation and t being the cooling down time.
In 1261 the decay power PS{ of the fission products resulting from fissions of each fissile
isotope i is approximated by the sum of 24 exponential functions. The energy release rating fj
(x) of the fission products of one single fission of isotope / is calculated as:
149
(1)
where x is the time elapsed since the fission occurred. The coefficients ay and Xij for fissions
of 235U, 238U, 239Pu, and 24lPu are listed in /26/. In case of fuel containing thorium, the
coefficients of the 235U are also applied for the bred 233U in subroutine NACHW.
Subdividing the power histogram of the fuel into a series of small time periods K with the
length Tk , the decay heat rating of the fission products related to one fission per second is
calculated as:
(2)
with tk being the cooling down time, i.e. the elapsed time since the end of the time interval
Integrating equation (2) using equation (1) for fj results in
and the contribution of the fission products to the decay power at time / after reactor shut
down amounts to
with
Pik Thermal power of the fissile isotope / during the time interval 7*
Qi Total thermal energy release of one fission of isotope i
T Total time of fuel operation in the reactor from the initial loading until the last
considered shut down.
A.6.2.3 Decay power of ""U and 239Np
The decay power Pb of these isotopes as precursors of the bred fissile material 239Pu,
respectively, is shown to be /26/:
150
being
~Q
and I = 239U and 239Np, respectively.
The values of F ; are to be calculated as follows:
The meanings of the variables are:
Eu mean decay energy of 239U (0.474 MeV)
EN p mean decay energy of 239Np (0.419 MeV)
Xv decay constant of 239U (4.91 * lCTVs)
XNp decay constant of 239Np (3.41 * 10"6/s)
.k i s the neutron capture rate of 238U, divided by the total fission rate during time interval K.
A.6.2.4 Contribution of neutron capture in fission products and in actinides
An example for neutron capture in the fission products and in the actinides is given in Fig. 12.
The contribution of the decay power due to neutron capture in fission products (Pe) is
calculated as a function of the decay time:
Pe(t,T) = PsU,T)*H(t)
using table 4 of/26/ for H(t).
Not included in Pe is the contribution of the long-lived fission product 134Cs. This is
calculated separately as Pcs according to equations 32 - 38 of Ref. I26I.
151
The contribution to the decay power by neutron capture in the heavy metal nuclides - except239U and 239Np - Pa(t,T), is evaluated on the basis of calculated nuclide concentration Ns (t,T):
28
pa (t, T) = £ Nt * A,. * £,. for i * 10 and 13
(i.e. 239U +239Np excluded, see A6.2.3)
being
Xi decay constant of heavy metal isotope i
E; mean decay energy of isotope i.
Cs- 1 33(stabil)
Cs- 1 34 —-—» Bo-1 34(stabil)
n,"y n,y n,7Cm-242—-L^- Cm-243—-L-* Cm-244—J-* Cm-245-
/
Am-241 n ' 7 > Am-242-
y
Am-244-
Pu-243-
Fig. 12: Neutron capture in fission products and in actinides
The total decay power is finally summed up:
0
152
References
IM E. Teuchert, U.Hansen, K.A. Haas:"VSOP-Computer Code System for Reactor Physics and Fuel Cycle Simulation"Kernforschungsanlage Julien, Jül-1649, (March 1980)
111 E. Teuchert, K.A. Haas, H.J. Rütten, H. Brockmann, H.Gerwin, U. Ohlig, W. Scherer"V.S.O.P. (94)Computer Code System for Reactor Physics and Fuel Cycle Simulation"Forschungszentrum Jülich, Jül-2897, (April 1994)
131 H. J. Rütten, K. A. Haas, H. Brockmann, U. Ohlig, W. Scherer:"V.S.O.P. (97)Computer Code System for Reactor Physics and Fuel Cycle Simulation"Forschungszentrum Jülich, Jül-3522, (April 1998)
141 H. J. Rütten, K. A. Haas, H. Brockmann, U. Ohlig, W. Scherer:"V.S.O.P. (99) for WINDOWS and UNIX"Computer Code System for Reactor Physics and Fuel Cycle Simulation"Forschungszentrum Jülich, Jül-3820, (October 2000)
151 H.J. Rütten:"The Depletion Computer Code ORIGEN-JUEL-Ü"Forschungszentrum Jülich, Jül-2739 (March 1993)
161 M.J. Bell:"ORIGEN - The ORNL isotope generation and depletion code"ORNL-4628, (May 1973)
HI H. Brockmann, U. Ohlig:"Erstellung einer 68-Gruppen-Wirkungsquerschnittsbibliothek im GAM-I Formatund einer 96-Gruppen-Wirkungsquerschnittsbibliothek im THERMALIZATION-Format"Kernforschungsanlage Jülich, Internal report, KFA-IRE-IB-3/89 (Juni 1989)
/8/ L.W. Nordheim, G.F. Kuncir:"A Program of Research and Calculations of Resonance Absorptions"General Atomic, San Diego, GA-2527 (1961)
191 R. Breitbarth, E. Teuchert:"Resonanzintegralberechnung für mehrfach heterogene Anordnungen"Kernforschungsanlage Jülich, Jül-551-RG (September 1968)
/10/ J.A. Young and J.U. Koppel:"Phonon Spectrum in Graphite"J. Chem. Physics, 42, 357 (1965)
153
/ l l / F. Carvalho:"Ineleastic scattering of thermal neutrons in graphite"Nucl. Sei. Eng. 34, 224 (1968)
/12/ G.Schröder:"COLLI: Ein Programmzyklus zur Berechnung der Streuung thermischer Neutronen"Kernforschungsanlage Jülich, JÜL-668-RG (1970)
/13/ CD. Joanou, J.S. Dudek:"GAM-I - A Consistant PI Multigroup Code for the calculation of Fast NeutronSpectra and Multigroup Constants"General Atomic GA-1850 (1961)
I\AI H.C. Honek:"THERMOS - A Thermalization Transport Theory Code for Reactor LatticeCalculation"Brookhaven National Laboratory, BNL-5826 (1961)
/15/ HJ. Neef, W. Scheren"THERMOS-JÜLModifikationen und Korrekturen am thermischen Zellprogramm Thermos"KFA Jülich, Internal Report IRE-IB-26/72
/16/ U. Hansen, E. Teuchert:"Influence of Coated-Particle Structure in Thermal Spectrum Energy Range"Nucl. Sei. Eng. 44, 12 (1971)
/17/ T.B. Fowler, D.R. Vondy, G.W. Cunningham:"Nuclear Reactor Core Analysis Code. CITATION"Oak Ridge National Laboratory, ORNL-TM-2496, Rev 2 (July 1971)
/18/ F.Todt:"FEVER - A One-Dimensional Few Group Depletion Program for Reactor Analysis"General Dynamics - General Atomic, GA-2749 (1962)
/19/ K. Petersen:"Zur Sicherheitskonzeption des Hochtemperaturreaktors mit natürlicher
Wärmeableitung aus dem Kern im Störfall"Kernforschungsanlage Jülich, Jül-1872 (Oktober 1983)
/20/ U.Hansen:"The VSOP System Present Worth Fuel Cycle Calculation; Methods and Codes;KPD"Atomic Energy Establishment Winfrith, England, Dragon project report 915 (1975)
UM H.J. Rütten, K.A. Haas:"NAKURE - Programm zur Berechnung der Nachzerfallsleistung derKernbrennstoffe von Hochtemperaturreaktoren mit kugelförmigen Brennelementen"Internal report KFA-ISR-IB-1/91, Januar 1991
154
1221 M. Schürenkrämer:"Theoretische und experimentelle Untersuchungen der Naturkonvektion im Kern desKugelhaufen-Hochtemperaturreaktors"Kernforschungsanlage Jülich GmbH, Jül-1912 (April 1984)
/23/ J.C. Cleveland, S.R. Greene:"Application of THERMEX-KONVEK Code to Accident Analyses of ModularPebble Bed High Temperature Reactors"Oak Ridge National Laboratory, ORNL/TM-9905 (August 1986)
/24/ L. Binkele, O. Hartleib, W. Biergans:"Wärmeleitfähigkeit und spez. el. Widerstand neutronenbestrahlter Graphite inAbhängigkeit von der Temperatur"Internal report, KFA-IRW-IB-13/75 (June 1975)
1251 VDI-Wärmeatlas:Berechnungsblätter für den WärmeübergangVerein Deutscher Ingenieure, Düsseldorf. 4. Auflage 1984
1261 "Berechnung der Nachzerfallsleistung der Kernbrennstoffe von Hochtemperatur-reaktoren mit kugelförmigen Brennelementen"DIN 25485 (1990), Deutsches Institut für Normung eV, Postfach 1107,D-1000, Berlin 30
/27/ K.Robold:"Wärmetransport im Inneren und in der Randzone von Kugelschüttungen"Kernforschungsanlage Jülich GmbH, Jül-1796 (August 1982)
/28/ J. Lieberoth, A. Stojadinovic:"Neutron Streaming in Pebble Beds"Nuc.Sci.Eng. 76, 336-344 (1980)
/29/ H. Gerwin:"Das zweidimensionale Reaktordynamikprogramm TINTE"Kernforschungsanlage JülichTeil 1: "Grundlagen und Lösungsverfahren", Jül-2167 (November 1987)Teil 2: "Anwendungsbeispiele", Jül-2266 (Februar 1989)
155
Forschungszentrum Jülichin der Helmholtz-Gemeinschaft
Jül-4189November 2005ISSN 0944-2952