1 / 37

Update of uncertainty file in the EAF project

FENDL Vienna 2-5.12.2008. Update of uncertainty file in the EAF project. J . Kopecky 1 and R.A. Forrest 2 1 JUKO Research, the Netherlands 2 EURATOM/ UKAEA Association, Culham, UK. Scope of the presentation. History of EAF uncertainties EAF-2007 status

lilli
Download Presentation

Update of uncertainty file in the EAF project

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. FENDL Vienna 2-5.12.2008 Update of uncertainty file in the EAF project J. Kopecky1 and R.A. Forrest21JUKO Research, the Netherlands 2EURATOM/UKAEA Association, Culham, UK

  2. Scope of the presentation • History of EAF uncertainties • EAF-2007 status • Current actions and updates for EAF-2009 • Conclusions

  3. History of EAF uncertainties • EAF-2 (1991) - Use of error factors derived from 14.5 MeV systematic to assess the overall uncertainty for most of the reactions. These factors indicate the accuracy estimate only.

  4. History of EAF uncertainties (cont.) • EAF-3 (1992) – Uncertainty file (in the ENDF-6 MF = 33 format) for all reactions in a 1 group (threshold reactions [Eth–20MeV]) and 2 groups (non-threshold reactions [10-5eV–EH–20MeV]) structure.The uncertainty values adopted correspond to the error factorf = 1 + , where  is the relative error (2 is actually stored in the file) in the cross section . In this case the best estimate of the cross section uncertainty is /f <  < f.

  5. MF33 from EAF-3

  6. Error factors (EAF-3)Used for all reactions

  7. History of EAF uncertainties (cont.) • EAF-4 (1995) • The uncertainty file contained more information for (n,g) and (n,f) reactions. These were given in 3 groups (10-5eV – EL – EH–20MeV]) rather than the 2 groups used previously. s(th) and Ig are used for respective energy groups as an experimental information. • Estimates for data with no experimental information :

  8. History of EAF uncertainties (cont.) • Experimental uncertainties for threshold reactions (EAF data set of recommended single s(14.5 MeV) data points) were used to replace the generic values adopted from 14.5 MeV systematic.This 14.5 MeV experimental variance is used for the entire energy range of the excitation function, and because no information on uncertainties at other energies has been used, an uncertainty value of  = 3exp has been adopted if  < 0.25.

  9. History of EAF uncertainties (cont.) • EAF-97 (1997) - New treatment of uncertainties for split reactions has been adopted. If experimental values of partial cross sections are known their uncertainties are adopted. However, in general for uncertainties of cross sections to ground and isomeric states an uncertainty factor of the branching ration (fb) of 1.3 has been adopted when the branching is derived from BR systematic. • By definition sm = bs and sg = (1-b)s. Thus the uncertainty in sm is given by and Ds/s=f-1 and Db=fb -1. Thus .

  10. History of EAF uncertainties (cont.) • EAF-99 (1999) - EAF-2001 (2001) - EAF-2003 (2003) - The corresponding uncertainty files for reactions supported with experimental data includes visual estimates based on all EXFOR data points covering broader energy regions. Three energy groups have been applied also for non-threshold (n,p) and (n,a) reactions with LEA treatment below EH.

  11. Error factors (EAF-97 – EAF-2003) Used for all reactions with no experimental support

  12. History of EAF uncertainties (cont.) • EAF-2005 – The upper energy extended to 60 MeV by adding a new energy group (20 – 60 MeV). The new error factors have been adopted using the assumption of f < 20MeV = f > 20 MeV. • Threshold reactions: Eth Eth 20MeV 20MeV 60MeV 60MeV

  13. Uncertainties – non-threshold reactions DIg or 0.5 EH EL Dexpor 0.5 20MeV Ds(th) 60MeV EL

  14. EAF-2007 status (exp. data) • Error factors based on experimental data -The experimental variance is adopted from both EXFOR as well the EAF internal data base and is used for the entire energy range of the excitation function. This uncertainty has been estimated visually from the data scatter around the library excitation curve from cross section plots or from growing data base of experimental integral validations.

  15. EAF-2007 status (no exp.) • Error factors with no experimental data (based on s(14.5 MeV) systematic - These factors quantify how well the systematic of s(14.5 MeV) fits the experimental data. They are used in the generation of the uncertainty file if no experimental data are available and the adopted library data agree with s(14.5 MeV). They are assumed to be representative for the whole energy range above Eth and the conservative values are applied. • Error factors f without systematic - If no cross section systematic is available, the results of graphical treatment of calculations or adopted data (e.g. SACS analysis) and the educated conservative guess of the accuracy of the calculation (upper limit) are applied. • Error factors f > 20 MeV - If no experimental data are available to estimate f, the same uncertainty is adopted as the uncertainty of the energy group just below 20 MeV.

  16. Error factors EAF-2007 Used for all reactions with no experimental support

  17. References prior EAF-2009 • J Kopecky, D Nierop and RA Forrest, Uncertainties in the european activation file EAF-3.1, ECN-C—94-015 (March 1994). • RA Forrest, J Kopecky, Statistical analysis of cross sections, Fusion Eng. and Design, 82(2007) 73. • RA Forrest, J Kopecky, AJ Koning, Detailed analysis of (n,p) and (n,a) reactions in the EAF-2007 and TALYS-generated libraries, Fusion Eng. and Design, 83(2007) 614. • RA Forrest, J Kopecky, AJ Koning, Revisions and improvements of neutron capture cross sections for EAF-2009 and validation of TALYS calculations, UKAEA FUS 546 (March 2008). • RA Forrest et al., Validation of EASY-2007 using integral experiments, UKAEA FUS 547 (April 2008). • Graphical comparison of all EAF-2007 Score =1-6 reactions, corresponding TALYS calculations (if TALYS not adopted) with experimental data (J Kopecky – private information).

  18. Actions and updates for EAF-2009Standard treatment • Complete revision of D(exp) values is in progress. 1. Threshold reactions. Several new rules shall be applied and new error factors proposed based on TALYS validations, SACS analysis and studies and comparisons against experimental data. • If D(Eth-20 MeV) < 0.1 the value has been increased to 0.1. • If D(Eth-20 MeV) < 0.2 the value of D(20 - 60 MeV) is specially treated for excitation curves with dominant cross section part above 20 MeV. Further the starting energy of the last energy group (20 MeV) has been made free, in order to better cover the the shape of the excitation curve. The use of C/E (integral data) for cross section improvements is added. • The uncertainty visual estimate shall be applied to account for the energy range of experimental data in the shape of calculated excitation curves (see next slide). • Details for threshold reactions are described together with relevant plots in EAF-Doc-48 to be released in December 2008.

  19. Uncertainty ranges1,2 – weak arguments 3,4,(5) – strong arguments for D 3 4 2 1 5

  20. Uncertainty ranges - a dream situation -

  21. Uncertainty rangesCombination of sources smax SACS EXP SACS Emaxsmax D1/2x TALYS validation TALYS predictive power

  22. TALYS predictive power for major threshold reactions (references) • Since EAF-2007 includes 65,565 reactions and only about 1750 are supported by experimental information, the test of predictive power with TALYS calculated data is of great importance. For validation of TALYS calculations see: • RA Forrest, J Kopecky, Statistical analysis of cross sections, Fusion Eng. and Design, 82(2007) 73. • RA Forrest, J Kopecky, AJ Koning, Detailed analysis of (n,p) and (n,a) reactions in the EAF-2007 and TALYS-generated libraries, Fusion Eng. and Design, 83(2007) 614. • RA Forrest, J Kopecky, AJ Koning, Revisions and improvements of neutron capture cross sections for EAF-2009 and validation of TALYS calculations, UKAEA FUS 546 (March 2008). • J Kopecky, TALYS validation plots, private communication.

  23. TALYS predictive power for major threshold reactions (cont.) • Neutron emission channels:For (n,n′), (n,xn), (n,n′p), (n,n′a) reactions the shape of calculated cross sections is in good agreement with the trend of experimental data and the shape of the excitation curves is considered good. • Charge particle emission channels: • A similar conclusion was made for (n,p), (n,d), (n,a) reactions with two specific cases to look at manually. The FWHM of (n,p) reaction in the transitional range with 0.08 < s < 0.12 and smaxof (n,a) reaction to check with SACS analysis. • (n,t) data are now in a good agreement. • (n,h), (n,2p) reactions are the experimental data (especially activation measurements) probably often overestimated due to the interference of competing reactions from elemental or isotopic impurities.

  24. SACS D½ (n,p) supported by EXFOR data Transition range Ca-43 Ti-47 Co-59 Zn-66 Fe-58 Ni-62 Cu-65 Zn-67, Pd-102,Sn-112

  25. (n,a) EAF-2009 (preliminary) Ir-194m

  26. (n,t) agreement

  27. (n,He-3) problem

  28. Actions and updates for EAF-2009 (cont.) Standard treatment • Details for threshold reactions are described together with relevant plots in EAF-Doc-48 to be released in December 2008. 2. Non-threshold reactions. Complete revision of D(exp) values shall be performed again with improved rules and new error factors based on TALYS validations, SACS analysis and agreements with experimental data: • Details will be described together with relevant plots in EAF-Doc-49 to be released by April 1, 2009. 3. Uncertainties with TALYS. The uncertainties for a limited number of important targets will be obtained by applying a Monte Carlo method to the TALYS nuclear reaction code.

  29. Error factors (EAF-2009) – in progress * - good fit with experimental data& - shape problems close to the threshold# - shape problem > 20 MeV for (n,p) and s(max) for (n,a)

  30. TALYS Monte Carlo calculation (A. Koning) • The uncertainties are obtained by applying a Monte Carlo method on the TALYS nuclear reaction code. In the fast neutron range, a number of model parameters (between 20 and 30) are varied. These parameters concern: • the optical model • the level density model • pre-equilibrium parameters • gamma-ray model parameters • They are randomly varied following a normal distribution, centered at the default TALYS parameter values, with a default standard deviation. These two values (mean and standard deviation) are obtained from a global fitting procedure so that a default TALYS calculation is rather close to experimental data from F-19 to Bi-209.

  31. TALYS Monte Carlo calculationFrom the scatter of (n,g) curves after 100 runs the average, error bars and all covariance matrix elements can be derived.

  32. MF33 of EAF-2009 for special targetsTALYS calculations

  33. Selection of important reactions for TALYS uncertainties • Importance diagrams with factor > 1 (90% of activation) used to select a number of reactions for most materials includes 5105 reactions. • Daughter nuclide Importance H1(n,g)H2 |H3 5 H2(n,g)H3 |H3 5 He3(n,p)H3 |H3 5 He4(n,2n)He3 |H3 5 Li6(n,p)He6 |He6 3 Li6(n,a)H3 |H3 5 Li7(n,na)H3 |H3 5 Li7(n,g)Li8 |Li8 5 Li7(n,d)He6 |He6 5 Be10(n,na)He6 |He6 3 Be10(n,g)Be11 |Li8, Be11 5 Be10(n,t)Li8 |Li8 3 Be10(n,nta)H3 |H3 1 Be9(n,nt)Li6 |He6 2 Be9(n,g)Be10 |H3, Li8, Be10, Be11 5 …… etc

  34. Selection of important reactions for TALYS uncertainties • From the previous list of important reactions a selections has been made for reactions on important targets for materials relevant to present fusion applications. This resulted in the following targets: Li6,7; Be9,10; B10,11; C12-14; N14,15; O16-18; F19, Sc45,46,48; Ti46-50; V49-51; Cr50,52-54; Mn54,55; Fe54,56-59; Co57-60; Ni58-64; Cu63,65; Nb91-95; Mo92,94-98,100; Ta177,179,180m-182; W180-186; Re185-187, Os184,186-190,192,193; Pb204-209; Bi203-209 and Po210. • This selection amounts to 100 isotopic targets and all reaction channels shall be calculated. The calculation will be performed in 1, 2 and 5 MeV energy steps for data for energy regions <20MeV, 20-30MeV and 30-60MeV, respectively. The results will be used to overwrite the existing error factors for these reaction channels.

  35. EAF v TALYS uncertainties Good agreement!

  36. EAF v TALYS uncertainties EAF uncertainty smaller – use of experimental data near threshold

  37. Conclusions • A thorough revision and extension of the uncertainty data for EAF-2009 by a combination of an improvement of error-factor approach and TALYS Monte Carlo calculation will be performed. This will result in a new high-quality uncertainty file in MF33 format. • The derivation of energy dependent error factors, based on TALYS Monte Carlo calculations, will be considered to replace the single-value (energy independent) factors f. • Use the results of the recent TENDL-2008 library and especially the future TENDL-2009 release.

More Related