From e7b7cd66c201865d89a128d6038c12d03d7b8cb4 Mon Sep 17 00:00:00 2001 From: Yunhuang Zhang Date: Mon, 20 Apr 2020 00:40:02 -0500 Subject: [PATCH] Raised THERMR Emax upper Limit from 10eV to 100eV We found that the internal limit of Emax<=10eV prevented us from obtaining accurate neutron slowing-down calculations (for instance, in Borated HDPE). The initial findings with numerical experiments were documented in our M&C2017 conference paper "NJOY Maximum Energy Limit For Thermal Neutron". A more in-depth proof based on two-body kinematic analyses is currently under review (2rd round) with Progress in Nuclear Energy. The added grid points were proposed by Dr. A. C. (Skip) Kahler. Raising emax to 100eV resulted in satisfying numerical simulation results compared to experiments, and we think it should be adequate for most applications. This limit can be further lifted by extending the grid but should be used with caution. --- src/thermr.f90 | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/src/thermr.f90 b/src/thermr.f90 index 6c40fe20..5e78a886 100644 --- a/src/thermr.f90 +++ b/src/thermr.f90 @@ -1508,7 +1508,7 @@ subroutine calcem(temp,itemp,iold,inew,ne,nex) real(kr)::b,diff,enow,ep,sabmin,tev,ylast real(kr)::u,xl,yl,sum real(kr)::tone,elo - integer,parameter::ngrid=118 + integer,parameter::ngrid=134 integer,parameter::nlmax=65 integer,parameter::nemax=5000 integer,parameter::mumax=300 @@ -1540,7 +1540,11 @@ subroutine calcem(temp,itemp,iold,inew,ne,nex) 1.42e0_kr,1.55e0_kr,1.70e0_kr,1.855e0_kr,2.02e0_kr,2.18e0_kr,& 2.36e0_kr,2.59e0_kr,2.855e0_kr,3.12e0_kr,3.42e0_kr,3.75e0_kr,& 4.07e0_kr,4.46e0_kr,4.90e0_kr,5.35e0_kr,5.85e0_kr,6.40e0_kr,& - 7.00e0_kr,7.65e0_kr,8.40e0_kr,9.15e0_kr,9.85e0_kr,10.00e0_kr/) + 7.00e0_kr,7.65e0_kr,8.40e0_kr,9.15e0_kr,10.00e0_kr,& + 11.00e0_kr,12.50e0_kr,14.00e0_kr,16.50e0_kr,19.50e0_kr,& + 23.00e0_kr,27.00e0_kr,31.50e0_kr,36.50e0_kr,42.00e0_kr,& + 48.00e0_kr,54.50e0_kr,61.50e0_kr,69.00e0_kr,77.00e0_kr,& + 88.50e0_kr,100.0e0_kr/) real(kr),parameter::unity=1.0e0_kr real(kr),parameter::sabflg=-225.e0_kr real(kr),parameter::eps=1.e-4_kr