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

correct diffusion coefficient is now used

Browse source
This commit is contained in:
Bryan Herman 2012-04-06 06:39:54 -07:00
parent bf9ae7045a
commit 652c040e8e
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/cmfd_data.F90
View file

@ -154,10 +154,10 @@ contains
cmfd % diffusion(h,i,j,k) = t % scores(score_index,4) % sum / flux
! calculate diffusion coefficient
cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*cmfd%totalxs(h,i,j,k))
! cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*cmfd%totalxs(h,i,j,k))
! cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*(cmfd % totalxs(h,i,j,k) -&
! & cmfd % p1scattxs(h,i,j,k)))
! cmfd % diffcof(h,i,j,k) = cmfd % diffusion(h,i,j,k)
cmfd % diffcof(h,i,j,k) = cmfd % diffusion(h,i,j,k)
else if (ital == 2) then
@ -380,7 +380,7 @@ contains
! compute dtilde
dtilde = (2*cell_dc*(1-ref_albedo))/(4*cell_dc*(1+ &
& ref_albedo)+(1-ref_albedo)*cell_hxyz(xyz_idx))
dtilde = 0.0_8
! dtilde = 0.0_8
else ! not next to a reflector or no core map
! compute dtilde