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

Only accumulate k-effective global tallies if eigenvalue calculation

Browse source
This commit is contained in:
Paul Romano 2015-06-12 00:37:39 +07:00
parent 64d3ccba9d
commit 2bdd579ba5
2 changed files with 17 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/physics.F90
View file

@ -253,15 +253,19 @@ contains
p % last_wgt = p % wgt
! Score implicit absorption estimate of keff
tally_absorption = tally_absorption + p % absorb_wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
if (run_mode == MODE_EIGENVALUE) then
tally_absorption = tally_absorption + p % absorb_wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
end if
else
! See if disappearance reaction happens
if (micro_xs(i_nuclide) % absorption > &
prn() * micro_xs(i_nuclide) % total) then
! Score absorption estimate of keff
tally_absorption = tally_absorption + p % wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
if (run_mode == MODE_EIGENVALUE) then
tally_absorption = tally_absorption + p % wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
end if
p % alive = .false.
p % event = EVENT_ABSORB

13
src/tracking.F90
View file

@ -1,5 +1,6 @@
module tracking
use constants, only: MODE_EIGENVALUE
use cross_section, only: calculate_xs
use error, only: fatal_error, warning
use geometry, only: find_cell, distance_to_boundary, cross_surface, &
@ -110,8 +111,10 @@ contains
call score_tracklength_tally(p, distance)
! Score track-length estimate of k-eff
tally_tracklength = tally_tracklength + p % wgt * distance * &
material_xs % nu_fission
if (run_mode == MODE_EIGENVALUE) then
tally_tracklength = tally_tracklength + p % wgt * distance * &
material_xs % nu_fission
end if
if (d_collision > d_boundary) then
! ====================================================================
@ -136,8 +139,10 @@ contains
! PARTICLE HAS COLLISION
! Score collision estimate of keff
tally_collision = tally_collision + p % wgt * &
material_xs % nu_fission / material_xs % total
if (run_mode == MODE_EIGENVALUE) then
tally_collision = tally_collision + p % wgt * &
material_xs % nu_fission / material_xs % total
end if
! score surface current tallies -- this has to be done before the collision
! since the direction of the particle will change and we need to use the