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

Update physical constants to CODATA 2014

Browse source
This commit is contained in:
Paul Romano 2016-11-30 07:28:05 -06:00
parent 08852a410b
commit ab72510d20
65 changed files with 927 additions and 925 deletions
Download patch file Download diff file Expand all files Collapse all files

15
openmc/data/data.py
View file

@ -208,16 +208,17 @@ def atomic_weight(element):
return None if weight == 0. else weight
# The value of the Boltzman constant in units of eV / K
# Values here are from the Committee on Data for Science and Technology
# (CODATA) 2010 recommendation (doi:10.1103/RevModPhys.84.1527).
K_BOLTZMANN = 8.6173324e-5
# (CODATA) 2014 recommendation (doi:10.1103/RevModPhys.88.035009).
# The value of the Boltzman constant in units of eV / K
K_BOLTZMANN = 8.6173303e-5
# Used for converting units in ACE data
EV_PER_MEV = 1.0e6
# Avogadro's constant from CODATA 2010
AVOGADRO = 6.02214129E23
# Avogadro's constant
AVOGADRO = 6.022140857e23
# Neutron mass from CODATA 2010 in units of amu
NEUTRON_MASS = 1.008664916
# Neutron mass in units of amu
NEUTRON_MASS = 1.00866491588

7
openmc/data/resonance.py
View file

@ -5,6 +5,7 @@ import numpy as np
from numpy.polynomial import Polynomial
import pandas as pd
from .data import NEUTRON_MASS
from .endf import get_head_record, get_cont_record, get_tab1_record, get_list_record
try:
from .reconstruct import wave_number, penetration_shift, reconstruct_mlbw, \
@ -196,7 +197,7 @@ class ResonanceRange(object):
ap = Polynomial((items[1],))
# Calculate channel radius from ENDF-102 equation D.14
a = Polynomial((0.123 * (1.008665*ev.target['mass'])**(1./3.) + 0.08,))
a = Polynomial((0.123 * (NEUTRON_MASS*ev.target['mass'])**(1./3.) + 0.08,))
return cls(target_spin, energy_min, energy_max, {0: a}, {0: ap})
@ -346,7 +347,7 @@ class MultiLevelBreitWigner(ResonanceRange):
competitive = items[3]
# Calculate channel radius from ENDF-102 equation D.14
a = Polynomial((0.123 * (1.008665*awri)**(1./3.) + 0.08,))
a = Polynomial((0.123 * (NEUTRON_MASS*awri)**(1./3.) + 0.08,))
# Construct scattering and channel radius
if nro == 0:
@ -606,7 +607,7 @@ class ReichMoore(ResonanceRange):
awri = items[0]
# Calculate channel radius from ENDF-102 equation D.14
a = Polynomial((0.123 * (1.008665*awri)**(1./3.) + 0.08,))
a = Polynomial((0.123 * (NEUTRON_MASS*awri)**(1./3.) + 0.08,))
# Construct scattering and channel radius
if nro == 0: