Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TENDL deuteron ACE formatting #184

Closed
psauvan opened this issue Nov 24, 2020 · 5 comments
Closed

TENDL deuteron ACE formatting #184

psauvan opened this issue Nov 24, 2020 · 5 comments

Comments

@psauvan
Copy link

psauvan commented Nov 24, 2020

There is a problem when processing new deuteron TENDL 2017 and 2019 libraries into ACE format.
These ace files produced by NJOY including photon emision do not work properly with MCNP6 deuteron transport.

This error is due to a missing reaction (MT number) in the acefc.f90 subroutine in the "if" sentence at line 5356.
In this instruction all the reactions producing the incident particle in the out channel should be listed. In this case, the reaction MT11 is missing.

To fix the problem, the correct instruction should be :

     if (mt.eq.2.or.mt.eq.11.or.mt.eq.32.or.mt.eq.35.or.&
       mt.eq.104.or.mt.eq.114.or.mt.eq.115.or.mt.eq.117.or.&
       mt.eq.157.or.mt.eq.158.or.mt.eq.169.or.&
       mt.eq.170.or.mt.eq.171.or.mt.eq.182.or.&
       mt.eq.183.or.mt.eq.185.or.mt.eq.187.or.&
       mt.eq.192.or.&
       (mt.ge.650.and.mt.le.699)) iskip=0

Patrick

@whaeck
Copy link
Member

whaeck commented May 13, 2021

@psauvan Can you point me to a TENDL file that I can use as a test case for this issue? Thanks.

@psauvan
Copy link
Author

psauvan commented May 13, 2021 via email

@jchsublet
Copy link

@whaeck @psauvan

There is some subtlety that you may need to account for regarding Acer charge particle processing. This is nicely explained from slides 25 of https://www-nds.iaea.org/index-meeting-crp/TM%20FENDL/docs/IAEA.FENDL_meeting2019.SKwon_rev3.pdf.

Regarding TENDL, TENDL-2019 supersedes TENDL-2017 (that had another more problematic format bug in the neutron production of mt=5) but may be more importantly for the matter in hand exist in more than one format, not at all equivalent regarding their respective processing for transport. For charge particles the s0 form (all implicit) is recommended as it mature nicely during the Acer processing steps that should all include a verifying 7 option, see https://tendl.web.psi.ch/tendl_2019/tar.html

Regarding the s30 forms (explicit up to 30 MeV, implicit afterward) Acer processing I gather that for charge particles this is very much in the making. My experience on some 283 stable targets of TENDL-2019 is that that it does run for the first Acer, but partially fails on the verifying ones and above all lead to weird MCNP messages when plugged into a realistic MCNP-6.2 simulation run.

@psauvan
Copy link
Author

psauvan commented May 14, 2021 via email

@whaeck
Copy link
Member

whaeck commented Nov 4, 2021

Fixed.

@whaeck whaeck closed this as completed Nov 4, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants