JEFFDOC-1160
Download
1 / 26

essing with NJOY-99 - PowerPoint PPT Presentation


  • 598 Views
  • Uploaded on

JEFFDOC-1160. Transport and Activation files processing with NJOY-99 J-Ch. Sublet & C Jouanne 1 CEA Cadarache, DEN/DER/SPRC , 13108 St Paul Lez Durance, France  1 CEA Saclay, DEN/DM2S/SERMA, 91191 Gif-sur-Yvette, France. A Tale of the unexpected. The main libraries; neutron and gamma.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'essing with NJOY-99' - ivanbritt


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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Slide1 l.jpg

JEFFDOC-1160

Transport and Activation files processing with NJOY-99J-Ch. Sublet& C Jouanne1 CEA Cadarache, DEN/DER/SPRC,13108 St Paul Lez Durance, France 1 CEA Saclay, DEN/DM2S/SERMA, 91191 Gif-sur-Yvette, France

A Tale of the

unexpected


The main libraries neutron and gamma l.jpg
The main libraries; neutron and gamma

  • JENDL-3.3/GP 337 isotopes (2 elements)

  • JEFF-3.1/GP 381 isotopes (7 elements)

  • JEFF-3.1/Th 9 compounds

  • JEFF-3.1/A 774 isotopes

  • ENDF/B-VII/GP 393 isotopes (3 elements)

  • ENDF/B-VII/Th 20 compounds

  • ENDF/B-VII/Ph 163 isotopes

  • + CENDL, BROND, …

  • + special purpose files, proton, deuteron, ..

  • + “preliminary” evaluations


The codes monte carlo and deterministic l.jpg
The codes: Monte Carlo and deterministic

  • The word “processing” encompasses a large set of processes

  • not all equivalent nor identical to one another

  • Not all Monte Carlo codes read in and/or need the same set of data or files, the processing level may differ

  • COG, MERCURY, MONK, MCBEND, MCNP5, MCNPX, TART, TRIPOLI, VIM, etc…

  • Not one cell or lattice codes use the same input libraries

  • APOLLO, CASMO, DRAGON, ERANOS, PANTHER, SCALE, TORT, DORT, WIMS, etc…

    Different modules, sequences, parameters are used every time


Monte carlo processing l.jpg
Monte Carlo processing

  • TRIPOLI-4.4.1: Monte Carlo, pointwise or groupwise

  • BROADR PENDF or GROUPR GENDF

  • CALENDF probability tables in the URR

  • 293.6, 300, 473.6, 600, 673.6, 873.6, 900, 1073.6, 273.6, 1473.6, 1673.6, 1873.6K

  • Modified THERMR

  • “TRIPOLI-4.4 JEFF-3.1 Based Libraries” CEA-R-6125

  • MCNP4c3&5: Monte Carlo, pointwise + Photonuclear (ENDF/B-VII)

  • ACER PENDF

  • PURR probability tables in the URR

  • 293.6, 300, 473.6, 600, 673.6, 873.6, 900, 1073.6, 273.6, 1473.6, 1673.6, 1873.6K

  • THERMR (up113, continuous secondary-energy distribution not yet plugged in)

  • “MCNP4c JEFF-3.1 Based Libraries” CEA-R-6110


Deterministic processing l.jpg
Deterministic processing

  • ERANOS2: 3D deterministic code, fast

  • NJOY GENDF* (+ MF-50) ( MERGE & GECCO)

  • ecco 1968, Xmas 172, VitJ 175, ecco 33 groups structures

  • 293.6, 573.6, 973.6, 1473.6, 2973.6K

  • CALENDF probability tables (MF-50 in all energy ranges)

  • P0-P1 to P5

  • “ECCOLIB-JEFF-3.1 Libraries” CEA-R-6100

  • APOLLO2: 2D deterministic cell code

  • NJOY GENDF* (+ as/afend,) up to 50 dilutions (LIBNJ90 & N2A2)

  • Universal 11276, Xmas 172, SHEM 281 groups structures

  • 273.6, 293.6, 573.6, 773.6, 973.6, 1473.6, 2273.6, 3473.6K

  • CALENDF probability tables (self shielding & mutual shielding)

  • P0-P1 to P5

  • DARWIN: activation/transmutation code

  • NJOY GENDF

  • 300 & 923K

  • 3, ecco33, sailor47, GAM-II 100, Xmas 172, 315 groups structures


Tripoli 4 4 1 specific processing l.jpg
Tripoli-4.4.1, specific processing

  • Treatment of thermal neutron scattering data, scattering kernels

  • in THERMR

    • Modification: results in terms of equi-probable cosines instead of equi-probable angles, this is not an update

    • With 32 cosines, done in up124

    • MT’s 221 and 222 hard coded

  • CALENDF-2005 probability table in, and only, the URR

    • addition of probability tables in the URR of each isotopes that contained one

      • 140 in JEFF-3.1

      • 252 in ENDF/B-VII

      • 209 in JENDL-3.3

  • PURR and CALENDF PT’s differs (order, steps, production methods, scope, …


  • Pointwise cross section comparison total l.jpg
    Pointwise cross section comparison: total

    A Cubic interpolation

    requires less points

    than a linear one

    But many more points

    exists in the CALENDF

    pointwise file in the URR,

    # smooth NJOY data

    CALENDF 115156 pts

    NJOY 72194 pts

    PT’s 11276 Gprs


    Pt s impact on the icsbep benchmarks l.jpg
    PT’s impact on the ICSBEP benchmarks

    Excellent way to

    test the influence

    of the URR

    PT’s 11276 Gprs

    2-3000 groups in

    The URR


    Slide9 l.jpg

    Tripoli-4.4.1 Libraries, internal processing

    TRIPOLI-4.4.1 generates XDR portable binary and angular

    distribution files the first time it reads in a PENDF and

    and an ENDF file

    Seven File Types are used :

    • Evaluation (ascii)

    • PENDF (binary)

    • Dictionary of cross sections (ascii)*

    • Binary XDR pointwise file*

    • Anisotropy file (ascii)*

    • Probability tables (ascii)

    • Thermal file (ascii)

    *: generated by TRIPOLI

    the first time it access the

    data file.

    140 URR JEFF-3.1

    252 URR ENDF/B-VII

    Tests for evaluations :

    ENDF-102 format

    Normalization (angular and energetic distribution)


    Tripoli 4 4 1 angular distribution l.jpg
    TRIPOLI-4.4.1, angular distribution

    • The angular distribution is computed by TRIPOLI-4.4.1 the first time a pendf file is accessed and the code require to be linked to the original endf file

    • The computed angular distribution is variable, file dependant and can contain up to 256 equally probable cosine bins

    • For the JEFF-3.1 U238 the angular distribution have been produced by the ECIS optical model code and tabulated in 91 bins from 0 to 180 degree. It is used as such by TRIPOLI-4.4.1

    • Elastic and inelastic channels structures are preserved

    • It is a much better representation than in Legendre coefficients, but not strictly identical to a 64 bins uniform processing (MCNP)



    Njoy 99 112 modules processing scheme l.jpg
    NJOY-99.112+ modules processing scheme

    • Modules Mainparameters Outputs

    • moder

    • reconr err = 0.001 (0.1%)

    • broadr errthn = 0.001 Thermal quantities

    • unresr

    • moder pendf files

    • thermr+ nbin = 32 tol = 0.001 emax = 4.95

    • moder

    Evaluations

    NJJOY-99.112+

    One single C-shell scripts

    Tripoli-4.4.1 library


    Njoy 99 161 modules processing scheme l.jpg
    NJOY-99.161 modules processing scheme

    • Modules Mainparameters Outputs

    • reconr err = 0.001 (0.1%)

    • broadr errthn = 0.001 Thermal quantities

    • unresr

    • moder pendf files

    • thermr nbin = 16 tol = 0.001 emax = 4.0

    • heatr mtk = 302 303 304 318 402 443 444

    • purr nbin = 20 nladr = 64

    • gaspr mt = 203-207

    • viewr QA graphs

    • moder

    • acer newfor = 1 iopp = 1 Ace files

    • acer iopt = 7

    • viewr QA graphs

    NJJOY-99.161

    MCNP4c3&5 library


    Qa graphs and thermal quantities l.jpg
    QA Graphs and thermal quantities

    NJOY-99

    Inter 7.0 or NJOY-99 Broadr




    Qa benchmarking independent l.jpg
    QA Benchmarking (independent)

    Regardless of the evaluations quality the same information is extracted from the data file with two different Monte Carlo code on very different benchmarks


    Reaction types and mt s in eaf 2005 a l.jpg
    Reaction types and MT’s in EAF-2005/A

    ENDF-102 MT: 1 …117

    new MT: 152…200

    declared unassigned


    Mt values l.jpg
    MT Values

    • Grid of reactions including all 36 MT numbers defined in ENDF

    • and some defined in EAF-2005


    Njoy 99 activation file processing l.jpg
    NJOY-99 Activation file processing

    • In NJOY-99.125 reconr has been modified to handle the new MT’s allowing pendf files generation

    • NJOY-99.68 has an automatic loop (10/) that will process all the nuclide production sections found in File 8

    • Results can be extracted from the output listing or read from the GENDF output

    • MATXSR knows how to add the multiple production sections together and generate a single production cross section for each product

    • Specialized output routines could be written for EAF

    • + all fissile isotopes format changed to be handled by /10 in JEFF-3.1/A and EAF-2005

    U230 pendf file

    9.223000+4 2.280580+2 0 0 1 19213 8 18 1

    0.000000+0 0.000000+0 3 0 0 259213 8 18 2

    Changed to

    9.223000+4 2.280580+2 0 0 1 19213 8 18 1

    1.000000+0 0.000000+0 3 0 0 259213 8 18 2


    Njoy 99 processing l.jpg
    NJOY-99 processing

    • Note that there may be multiple sections generating the same ZAm, and they have to be added together for total nuclide production

    ZAm of product Fe54 to Mn52m (n,t)

    

    2.605400+4 2.505210+5 1 1 0 12625 3105 37

    2.931600+2 0.000000+0 2 1 2 12625 3105 38

    4.39594+16 1.172175-7 2625 3105 39

    

    cross section

    • An update for groupr, matxsr, etc.. is needed for those modules to be able to handle properly the new MT’s


    Activation libraries the way forward l.jpg
    Activation libraries: the way forward

    • Elastic scattering and inelastic channels to be added, it allows to reconstruct and compare the total with experimental data, a must ..

    • Positive Q, (n,p) and (n,a) branching ratio in MF-9

    • Upper energy limit; 20 Mev , 60 Mev, … 150 MeV

    • MT-5/MF-6 yields

      • above 60 Mev

      • or for all, but the recognized ENDF-6 MT’s channels

         pure ENDF-6 format

    • MT-5 and activation yields in a separate file, for each isotopes

    • Deuteron, Proton, Gamma activation files


    Activation libraries the way forward23 l.jpg
    Activation libraries: the way forward

    • Uncertainties

      • In a single for EAF-2005

      • MF-33 like

        • + one comment line

        • + isomeric MAT

      • Provide a measure of the accuracies without correlations

      • simple MF-33, groups variances without correlations between the cross sections and/or the adjacent groups

      • Variable group structure

        It provide a complete set of “exploitable” data for activation,

        transmutation calculation uncertainties related to and uniquely

        the cross sections, not the reaction rates nor the neutron flux


    Processing pits l.jpg
    Processing pits

    • If the processing fails this may be due to:

      • the inputted data file (evaluation format, numeric,…

      • the process code input parameter set (P9, 128 bins,…

      • the processing code itself

    • If it succeed it may well still fail to be interpreted properly or adequately (i.e. negative PT’s, wrong sum, non positive coefficient, …) but can this be spotted, checked, and finally stamped by external verification

    • But if it succeed and lead to different results (processed file) on different platforms (computers) when using a same input deck and evaluated file who will see it but a sturdy, exemplified, thorough QA and issue tracker as it has been developed and openly shared at:

    • http://t2.lanl.gov/codes/njoy99/index.html


    Conclusions l.jpg
    Conclusions

    • NJOY is a complex code that may be dragged (less than before) into the “grey” areas of the different OS/compilers handling of number

    • Its installation requires care and knowledge with a high level of QA, Verification and follows up

    • Its input parameter sets may influence the end results, but this is a requirement

    • There are other processing codes that can be used to ascertain and verify some aspect of the processing

    • It evolved significantly over the past years, but still requires care and attention from its “knowledgeable” users


    A spell expressed during bob s symposium l.jpg
    A spell, expressed during Bob’s Symposium

    • It has been, and still is, a rare pleasure to follow the footsteps of such exceptional Wizard, as like to call him the man from the Dorset shore, he used in fact the word Magician but we differ on those delicate language manners. I would prefer the word Enchante(u)r that suits more my mind, although who am I to spell it with a u? any how it encompasses in many better ways all I know he have given us. Other did say Virtuoso, Maestro …

    • The Old Continent still needs you Bob, and may be your last spell of the 99 series, a 1XX should do fine. A rune that will make happy the many few able to fell its power, and use it wisely in accordance with your custodian masterminding.

    • Take pleasure in knowing that your legacy, did and will for many years to come entice the four corners of the world, at least the ones I know …


    ad