crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 06:05:58 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Added restriction energy constraint on sampling Maxwell fission spectrum.

Browse source
This commit is contained in:
Paul Romano 2011-11-29 15:16:29 -05:00
parent 627870fec8
commit 68d229dbf8
1 changed files with 23 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

27
src/physics.f90
View file

@ -1711,13 +1711,32 @@ contains
! =======================================================================
! MAXWELL FISSION SPECTRUM
! read number of interpolation regions and incoming energies
NR = int(edist % data(1))
NE = int(edist % data(2 + 2*NR))
! determine nuclear temperature from tabulated function
T = interpolate_tab1(edist % data, E_in)
! sample maxwell fission spectrum
E_out = maxwell_spectrum(T)
! TODO: Add restriction energy constraint??
! determine restriction energy
lc = 2 + 2*NR + 2*NE
U = edist % data(lc + 1)
n_sample = 0
do
! sample maxwell fission spectrum
E_out = maxwell_spectrum(T)
! accept energy based on restriction energy
if (E_out <= E_in - U) exit
! check for large number of rejections
n_sample = n_sample + 1
if (n_sample == MAX_SAMPLE) then
message = "Too many rejections on evaporation spectrum."
call fatal_error()
end if
end do
case (9)
! =======================================================================