ML22277A789
Text
IDENTIFICATION REVISION
FFDC00561 E2 4.0 I 11 I AREVA NP A Fuel BU AREVA
TOTAL NUMBER OF PAGES: 15
Safety-Criticality Study Fresh UO2 fuel transportation casks Determining the uncertainties to apply to the criticality studies performed using CRISTAL
ADDITIONAL INFORMATION:
Criticality - Uncertainties - Shipping cask
PROJECT DISTRIBUTION TO PURPOSE OF DISTRIBUTION
HANDLING Public For application I I For action CATEGORY STR - Study Report For information
STATUS BPE
ROLES This text shall no be visible - Adjust frames over to ensure signature block will completely cover this text NAMES DATES ORGANIZATIONS SIGNATURES
AUTHOR I I REVIEWER I I I
APPROVER I I This text shall no be visible - Adjust frames over to ensure signature block will completely cover this text
Classification Export AL: 0E001 ECCN: N RELEASE DATA: Les marchandises portant la désignation "AL inégal N" sont soumises la réglementation européenne ou allemande en matire de contrle des exportations au sein ou hors de l'UE. Les marchandises portant la désignation "ECCN inégal N" sont soumises la réglementation américaine. Les marchandises portant les désignations "AL:N" ou "ECCN:N" peuvent, selon la destination ou l'utilisation finales du produit, également tre soumises autorisation.
Export classification AL: 0E001 ECCN: N Goods labeled with AL not equal to N are subject to European or German export authorization when being exported within or out of the EU. Goods labeled with ECCN not equal to N are subject to US reexport authorization. Even without a label, or with label AL: N or ECCN: N, authorization may be required due to the final whereabouts and purpose for which the goods are to be used.
CHANGE CONTROL RECORDS: France: Y Exportkennzeichnung AL: 0E001 ECCN: N Die mit "AL ungleich N" gekennzeichneten Güter unterliegen bei der Ausfuhr aus der EU bzw.
This document, when revised, must be USA: N innergemeinschaftlichen Verbringung der europischen bzw. deutschen Ausfuhrgenehmigungspflicht. Die mit reviewed or approved by following regions: Germany: N "ECCN ungleich N" gekennzeichneten Güter unterliegen der US-Reexportgenehmigungspflicht. Auch ohne Kennzeichen, bzw. bei Kennzeichen "AL: N" oder "ECCN: N", kann sich eine Genehmigungspflicht, unter anderem durch den Endverbleib und Verwendungszweck der Güter, ergeben.
Non-proprietary document
CW01 Rev. 4.5 - 08/03/17 N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 2/15 AREVA
REVISIONS
DATE INDEX OBSERVATIONS
19/03/2004 B Additional information (Translation of French version to English).
13/05/2013 C Update for reference to CRISTAL V1.1 (synchronisation with the French version of Revision C)
See first 4.0 Update of the template of the document page Numerical revision number
Addition of the new statistical uncertainty calculated with I I (§4.2.2.2,
§4.2.3 and §5) (synchronisation with the French version of Revision 4.0)
- Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 3/15 AREVA
CONTENTS
- 1. INTRODUCTION...................................................................................................................................................... 7
- 2. REGULATIONS........................................................................................................................................................ 7
- 3. THECRISTALV1PACKAGE....................................................................................................................................... 8 3.1 CIGALES..................................................................................................................................................................... 8 3.2 APOLLO2.................................................................................................................................................................... 8 3.3 MORET4..................................................................................................................................................................... 8
- 4. DETERMININGCALCULATIONUNCERTAINTIES........................................................................................................ 9 4.1 PRELIMINARIES............................................................................................................................................................... 9 4.2 APPLICATIONOFUNCERTAINTIESTOTHECALCULATIONRESULTS.............................................................................................. 9 4.2.1 Systematicuncertainty..........................................................................................................................................9 4.2.2 Statisticaluncertainty......................................................................................................................................... 10 4.2.3 Combinedoveralluncertainty............................................................................................................................. 10
- 5. CONCLUSION........................................................................................................................................................ 12
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 4/15 AREVA
LIST OF TABLES
Table 1 : APOLLO2-MORET4 results for UO 2 rod arrays in transportation configurations.................... 13
LIST OF FIGURES
Figure 1 : Flow diagram explai ning the CRISTAL package................................................................... 14 Figure 2 : Calculation/measurement comparisons of the various calculation codes used..................... 15
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 5/15 AREVA
REFERENCES
[1] TS-R-1 - IAEA Safety Standards Series -
Regulations for the Safe Transport of Radioa ctive Material - 1996 Edition (As Amended 2003)
[2] TS-G-1.1 (ST2) - IAEA Safety Standards Series Advisory Material for the IAEA Regulations fo r the Safe Transport of Radioactive Material
[3] Framatome-ANP France Critic ality Calculation Route - Qualification and Safety Margin Determination - M. Doucet and al - ICNC'2003 - Tokai Mura - Japan - October 20-24 2003
[4] NEA/NSC DOC(95)03 - September 2002 Edition International Handbook of Evaluated Critica lity Safety Benchmark Experiments (ICSBEP)
[5] DSU/SEC/T/2007 Indice A Contribution la qualification de la voie standard APOLLO2-MORET4 du formulaire de criticité CRISTAL Synthse des résultats de calcul obtenu avec la version V1.1 du formulaire CRISTAL.
[6] DSMR/96-393 - Eléments nécessaires l'exper tise des études de criticité des emballages chargés de matires fissiles.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 6/15 AREVA
SUMMARY
The purpose of this document is to determine the uncertainties to apply to the criticality studies performed using CRISTAL for casks used to transport fresh UO 2 fuel assemblies or rods. The method used is based on IAEA regulations, and is generally valid for all Safety Authorities that apply IAEA regulations for the transportation of fissile material.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 7/15 AREVA
- 1. INTRODUCTION Since 2002, a safety analysis report must be submitt ed for casks transporting fuel material, and must meet all regulations [1]. While no changes have been made regarding statutory mechanical and heat tests, justification of the results of criticality stud ies is, however, now formally required, as explained in
[2]. This justification must be based on comparisons between calculations and measurements that are representative of the configurations under study, in or der to establish a reacti vity bias, which is then included in the unrefined results.
This document describes the method used and the uncerta inties obtained on the basis of [2], illustrated in [3] and applied in the CRISTAL package.
- 2. REGULATIONS Appendix VII of [2] entitled "Criticality Safety Assessments" gives a detailed description of the procedure for determining this uncertainty. The me thod used is briefly described below. The Upper Safety Limit (USL) is defined as follows:
USL = 1.00 - km - ku /1/
Where km is the "administrative" safety margin and cannot be compressed. In typical practice this is 5% k, And where ku is the reactivity bias from the calculat ions/measurements comparison for experiments close to the configurations to be evaluated by th e criticality code package. This uncertainty depends on the calculation sequence used. One of the aims of this document is to describe the approach used for the CRISTAL package.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 8/15 AREVA
- 3. THE CRISTAL V1 PACKAGE The CRISTAL package is made up of three main codes: CIGALES, APOLLO2 and MORET4.
A flow diagram explaining CRISTAL is shown in Figure 1.
3.1 CIGALES CIGALES is a code for processing physical data (fissile material, geometric shape, density, concentration, etc.) and which prepares datasets for APOLLO2. Its main functionalities are:
calculation of atomic concentrations for fissile environments, based on theoretical or experimental dilution laws,
calculation of atomic concentrations for non-fissile environments and/or use of a predefined materials database (steel, concrete, etc.),
automated preparation of the APOLLO2 datasets, to perform:
o calculations based on collision probabiliti es (Pij), which produce neutron parameters such as k, B²m, etc., and homogenous and self-s hielded macroscopic cross-sections for 3D MORET4 code,
o reactivity calculations using the Sn method with 1D geometry, used ma inly for criticality standards calculations.
3.2 APOLLO2 APOLLO2 is a cell / fuel assembly deterministic c ode for spectral calculations. It can use a variety of flow calculation types, such as the collision probabili ty method (Pij) and the Sn method. The Pij solver approximates currents at the interfaces, making an exact 2D geometry. The Sn solver uses a diamond system. APOLLO2 uses an explicit self-shielding processing method, based on equivalence theories. It is also capable of providing self-shielded multig roup cross-sections for resonating isotopes. The version of APOLLO2 built into the CRISTAL V1.1 pa ckage uses sets of microscopic cross-sections from the JEF2.2 evaluation. The energy mesh used comprises 172 groups, making it possible to process correctly the neutron spectrum from 0 to 10 MeV. The homogenized and self-shielded macroscopic cross-sections feat uring 172 energy groups are then sent to the MORET4 3D Monte Carlo code.
3.3 MORET4 MORET4 is a 3D calculation code that uses the Monte Carlo method to calculate the reactivity (Keff) in configurations of varying degrees of complexity, and the reaction rates in the different volumes that make up the configuration. Thanks to the aforementioned APOLLO2 calculations, MORET4 is well suited to the safety/criticality studies needed to justify fuel installations and the transportation of fissile material.
MORET4 uses the three conventional estimators (t rack length, collision and absorption) to determine Keff. These three estimators are combined with each sampling result, and a more general combination is used at the end of the calculation.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 9/15 AREVA
- 4. DETERMINING CALCULATION UNCERTAINTIES
4.1 PRELIMINARIES Calculation uncertainties reflect the quality of the code sequence used to calculate the reactivity of a configuration. The codes used by AREVA NP to obtai n a transportation authorization for fissile material are part of the CRISTAL V1.1 package, described in the previous section, and were the subject of a qualification file [5]. The qualification file is prim arily based on the results of the international work group ICSBEP [4] and is supplemented by experiment s performed by France s CEA-Valduc and CEA-Cadarache. The file covers virtually the whole fuel cycle, from enric hment to reprocessing.
4.2 APPLICATION OF UNCERTAINTIES TO THE CALCULATION RESULTS The overall uncertainty value is obt ained using the following formula:
I bias2 t 222 M Keff sexp Where:
I M = uncertainty to be added to the calculated value Keffbias = systematic uncertainty derived from calculation/measurement comparisons (ignored if positive value) s = statistical uncertainty derived from the calculation, exp = statistical uncertainty deriv ed from laboratory measurements, t = statistical uncertainty deriv ed from manufacturing tolerances.
4.2.1 SYSTEMATIC UNCERTAINTY Systematic uncertainty is that which is deduc ed from Keff calculation/ measurement comparisons.
Regarding rod arrays in transportation configurations, the qualification file [5] includes calculation/measurement comparisons based on the re sults of the work group ICSBEP [4]. This is available for consultation for further information on laboratory conditions.
The experimental results used in this analysis concern LEU COMP THERM 010 - 017 - 027 - 040 configurations, which include 25 configurations representing fresh UO 2 fuel transportation, including various 235U enrichment levels, neutron-absorbing screens, reflectors and rod arrays.
The APOLLO2 - MORET4 results (from [5]) using these configurations are given in Table 1.
When we look at the results in this table, we see that, with one exception, the CRISTAL package overestimates the reactivity of the configurat ions under study. The only result that gave an underestimation was also analyzed using t he other calculation codes and methods:
KENO IV - 27 energy groups - ENDF/BVI library:
Keff = I I MCNP - continuous energy - ENBF/BV library:
Keff = I I This underestimation would therefore seem to apply in general, and higher uncertainties than the one mentioned in [5]( I I) must be taken into account when interpreting laboratory measurements.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 10/15 AREVA
4.2.2 STATISTICAL UNCERTAINTY There are three types of statistical uncertainty:
s = statistical uncertainty derived from the calculation, exp = statistical uncertainty derived from laboratory measurements, t = statistical uncertainty of the quadratic co mbination of the manufacturing tolerances.
4.2.2.1 EXPERIMENTAL UNCERTAINTY (exp)
The four experiment groups chosen for qualification each include a limited number of independent experiments. Statistically, we cannot therefore per form a quadratic combination of the measurement uncertainties. Instead, we use the arithmetic averag e of the uncertainties of all the experiments. An examination of Figure 2 confirms that the CRISTAL results overestimate the laboratory measurements, including uncertainties.
The average experimental uncer tainty is calculated on the basis of the values in Table 1, and is I I.
4.2.2.2 STATISTICAL UNCERTAINTY DERIVED FROM CALCULATIONS ( S)
All MORET4 results used in our studies are given with a maximum standard deviation of:
S = I I, for the studies performed before july 2013, S = I I, for the studies performed after july 2013.
4.2.2.3 MANUFACTURING TOLERANCE UNCERTAINTY (T)
This uncertainty is covered by the recommendations in document [6], which in criticality studies are applied to the pellet-gap and cladding dimensions.
When a neutron-absorbing poison is us ed, manufacturing tolerances are included conservatively in the studies.
In the special case where borated resin (FS69 type) is used for the FCC, the average density of the product is I I. Each resin batch is analyzed to det ermine the variability of the product in terms of meeting the isotopic specification and actual density per unit volume. This variability is I I. Worth -
noting in these calculations is the penalizing density of I I -
4.2.3 COMBINED OVERALL UNCERTAINTY In section 4.2.1, we saw how there was a systematic positive deviation. This is an average value that must be corrected by the experimental and theoretical values, to verify whether the systematic bias is always positive. The verification is made using the following equation:
bias Keffbias* 2 cal22 Keff exp Where:
Keff bias* = Average C-M, as shown in Table 1 (Em) (I I) exp = Average experimental standard deviat ion as shown in Table 1 (I I) cal = Average theoretical standard deviation as shown in Table 1 (I I)
This leads to a systematic corrected bias ( Keff bias) of I Ipcm. This positive value is therefore not taken into account in the combined overall uncertainty. -
To support this systematic zero bias, Figure 2 shows how the results of CRISTAL V1.1 rate against the measurements and results of other calculation codes used throughout the world.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 11/15 AREVA
Overall uncertainty IM (as defined in section 4.2), can therefore be defined as follows:
I 2 t* 222 M exp s
Where, in this case :
s represents the standard deviation of industrial MORET4 calculations :
o s = I I pcm, for the studies performed before july 2013, o s = I I pcm, for the studies performed after july 2013, exp = I I pcm.r I t I I pcm.
From this we can deduce an uncertainty of:
I I pcm,. for the studies performed before july 2013 with s = I I pcm, I I pcm,. for the studies performed after july 2013 with s = I I pcm, Here, IM represents kU from equation /1/, with a system atic zero bias, as shown above.
Overall combination from equation /1/ gives us:
USL = 1.00 - km - kU for the studies performed before july 2013 with s = I I pcm:
o USL = I I o USL = I I I for the studies performed after july 2013 with s = I I pcm o USL = I I o USL = I I I
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 12/15 AREVA
- 5. CONCLUSION This document demonstrates the safety/criticality approach used by AREVA NP in its application for authorization to transport fresh UO 2 fuel with a maximum 235U enrichment of 5%, as a means of complying with IAEA regulations [1] and [2].
The reactivity uncertainty to apply to the CRISTAL results is:
I I pcm, for the studies performed before july 2013 with s = I I pcm, I I pcm, for the studies performed after july 2013 with s = I I pcm, which validates the CRISTAL package for criticality studies with the purpose to transporting fresh UO 2 fuel with a maximum 235U enrichment of 5%.
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 13/15 AREVA
Table 1 : APOLLO2-MORET4 results for UO2 rod arrays in transportation configurations
Proprietary table
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 14/15 AREVA
Figure 1 : Flow diagram explaining the CRISTAL package
FORMULAIRE CRISTAL en
- W w z
CIGALES JEF 2.2 z Donnees 0 0
LOIS DE DILUTION nucleaires en w
EGM PREAPOL 0 Windows NT 1-z PROCEDURES APROC THEMIS/NJOY THEMIS/NJOY w
~
w I-1k'.
I-
CESAR4 DARWIN CEA93 Sections Sections Ponctuellell
...J llllultigroupe* Format pendf
~
~
1-TRIPOLI 4 w
0 Monte-Carlo z Energie continue 0
~
I-en K.11 (3D)
Taux de reactions
VOIE DE REFERENCE
Formulaire CRISTAL - CRISTAL package PC utilisateur - User PC Lois de dilution - Dilution laws Procédures APROC - APROC procedures Données nucléaires - Nuclear data Traitement des données - Data processing Station de travail - Workstation Méthode Pij - Pij method Sections multigroupes - Multigroup cross-sections Sections ponctuelles format pendf - Regular Pendf-format cross-sections
Énergie continue - Continuous energy Taux de reactions - Reaction rates Dimensions Keff données - Given Keff dimensions Flux-Curants - Flux-Current Voie standard - Standard route Voie de reference - Reference channel
Non-proprietary document
N° FFDC00561 E2 Rev. 4.0 STR - Study Report A Handling: Restricted AREVA Page 15/15 AREVA
Figure 2 : Calculation/measurement comparisons of the various calculation codes used C-M Comparison with Uncertainties
Proprietary drawing
Non-proprietary document