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

Changed calculation of scattering angle in elastic scattering to use dot product.

Browse source
This commit is contained in:
Paul Romano 2012-07-24 17:42:27 -04:00
parent 6d0d9369a3
commit 42f0ac75d4
1 changed files with 3 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/physics.F90
View file

@ -562,8 +562,9 @@ contains
E = dot_product(v_n, v_n)
vel = sqrt(E)
! convert cosine of scattering angle from CM to LAB
mu = (ONE + awr*mu)/sqrt(awr**2 + TWO*awr*mu + ONE)
! compute cosine of scattering angle in LAB frame by taking dot product of
! neutron's pre- and post-collision angle
mu = dot_product(p % coord0 % uvw, v_n) / vel
! Set energy and direction of particle in LAB frame
p % E = E