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

added more parameters, for non accelerated regions, albedo rejection and zero flux

Browse source
This commit is contained in:
Bryan Herman 2012-10-17 06:59:58 -07:00
parent 31c08efeee
commit e97fa4da41
2 changed files with 12 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/cmfd_data.F90
View file

@ -469,7 +469,7 @@ contains
subroutine compute_dtilde()
use constants, only: CMFD_NOACCEL
use constants, only: CMFD_NOACCEL, ZERO_FLUX, TINY_BIT
use global, only: cmfd, cmfd_coremap
integer :: nx ! maximum number of cells in x direction
@ -546,7 +546,7 @@ contains
& (1-albedo(l))*cell_hxyz(xyz_idx))
! check for zero flux
if (albedo(l) - 999.0_8 < 1.e-8_8) dtilde = 2*cell_dc/cell_hxyz(xyz_idx)
if (albedo(l) - ZERO_FLUX < TINY_BIT) dtilde = 2*cell_dc/cell_hxyz(xyz_idx)
else ! not a boundary
@ -562,7 +562,7 @@ contains
if (cmfd_coremap) then
if (cmfd % coremap(neig_idx(1),neig_idx(2),neig_idx(3)) == &
& 99999 .and. cmfd % coremap(i,j,k) /= 99999) then
& CMFD_NOACCEL .and. cmfd % coremap(i,j,k) /= CMFD_NOACCEL) then
! get albedo
ref_albedo = get_reflector_albedo(l,g,i,j,k)
@ -758,7 +758,8 @@ contains
function get_reflector_albedo(l, g, i, j, k)
use global, only: cmfd, cmfd_hold_weights
use constants, only: ALBEDO_REJECT
use global, only: cmfd, cmfd_hold_weights
real(8) :: get_reflector_albedo ! reflector albedo
integer :: i ! iteration counter for x
@ -780,7 +781,7 @@ contains
! calculate albedo
if ((shift_idx == 1 .and. (current(2*l )-0.0_8) < 1.0e-10_8) .or. &
(shift_idx == -1 .and. (current(2*l-1)-0.0_8) < 1.0e-10_8)) then
albedo = 999.0_8
albedo = ALBEDO_REJECT
cmfd_hold_weights = .true.
else
albedo = (current(2*l-1)/current(2*l))**(shift_idx)

6
src/constants.F90
View file

@ -371,4 +371,10 @@ module constants
! for non-accelerated regions on coarse mesh overlay
integer, parameter :: CMFD_NOACCEL = 99999
! constant to represent a zero flux "albedo"
real(8), parameter :: ZERO_FLUX = 999.0_8
! constant to represent albedo rejection
real(8), parameter :: ALBEDO_REJECT = 999.0_8
end module constants