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

disallow target energy greater than is used in SIGMA1

Browse source
This commit is contained in:
walshjon 2018-02-17 14:09:37 -08:00
parent 4d240dc303
commit 9ec044ec2e
1 changed files with 7 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/physics.F90
View file

@ -925,11 +925,14 @@ contains
maxval(nuc % elastic_0K(i_E_low + 1 : i_E_up)), xs_up)
DBRC_REJECT_LOOP: do
! sample target velocity with the constant cross section (cxs) approx.
call sample_cxs_target_velocity(nuc, v_target, E, uvw, kT)
TARGET_ENERGY_LOOP: do
! sample target velocity with the constant cross section (cxs) approx.
call sample_cxs_target_velocity(nuc, v_target, E, uvw, kT)
E_rel = dot_product((v_neut - v_target), (v_neut - v_target))
if (E_rel < E_up) exit TARGET_ENERGY_LOOP
end do TARGET_ENERGY_LOOP
! perform Doppler broadening rejection correction (dbrc)
E_rel = dot_product((v_neut - v_target), (v_neut - v_target))
xs_0K = elastic_xs_0K(E_rel, nuc)
R = xs_0K / xs_max
if (prn() < R) exit DBRC_REJECT_LOOP