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

revert the function calulcating fissionable mass

Browse source
This commit is contained in:
liangjg 2019-01-02 10:39:23 -05:00
parent 0a4a1e7e81
commit 3846168e70
1 changed files with 11 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

18
openmc/material.py
View file

@ -249,11 +249,15 @@ class Material(IDManagerMixin):
@property
def fissionable_mass(self):
mass = 0.0
for nuc, _, _ in self._nuclides:
if openmc.data.zam(nuc)[0] >= 90:
mass += self.get_mass(nuc)
return mass
if self.volume is None:
raise ValueError("Volume must be set in order to determine mass.")
density = 0.0
for nuc, atoms_per_cc in self.get_nuclide_atom_densities().values():
Z = openmc.data.zam(nuc)[0]
if Z >= 90:
density += 1e24 * atoms_per_cc * openmc.data.atomic_mass(nuc) \
/ openmc.data.AVOGADRO
return density*self.volume
@classmethod
def from_hdf5(cls, group):
@ -715,9 +719,9 @@ class Material(IDManagerMixin):
"""
mass_density = 0.0
for nuc, atoms_per_cc in self.get_nuclide_atom_densities().values():
density_i = 1e24 * atoms_per_cc * openmc.data.atomic_mass(nuc) \
/ openmc.data.AVOGADRO
if nuclide is None or nuclide == nuc:
density_i = 1e24 * atoms_per_cc * openmc.data.atomic_mass(nuc) \
/ openmc.data.AVOGADRO
mass_density += density_i
return mass_density