crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Calculate density_ in Material exactly as before

Browse source 
Without this change, the cmfd_feed_2g and cmfd_feed_ng tests fail due to
slightly different answers.
This commit is contained in:
Paul Romano 2019-01-28 10:46:23 -06:00
parent 170f4a0a62
commit 72a7c23199
1 changed files with 3 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

7
src/material.cpp
View file

@ -352,10 +352,9 @@ void Material::finalize()
void Material::normalize_density()
{
bool percent_in_atom = (atom_density_(1) > 0.0);
bool percent_in_atom = (atom_density_(0) > 0.0);
bool density_in_atom = (density_ > 0.0);
double sum_percent = 0.0;
for (int i = 0; i < nuclide_.size(); ++i) {
// determine atomic weight ratio
int i_nuc = nuclide_[i];
@ -376,7 +375,7 @@ void Material::normalize_density()
// Change density in g/cm^3 to atom/b-cm. Since all values are now in
// atom percent, the sum needs to be re-evaluated as 1/sum(x*awr)
if (!density_in_atom) {
sum_percent = 0.0;
double sum_percent = 0.0;
for (int i = 0; i < nuclide_.size(); ++i) {
int i_nuc = nuclide_[i];
double awr = settings::run_CE ?
@ -384,7 +383,7 @@ void Material::normalize_density()
sum_percent += atom_density_(i)*awr;
}
sum_percent = 1.0 / sum_percent;
density_ *= -N_AVOGADRO / MASS_NEUTRON * sum_percent;
density_ = -density_ * N_AVOGADRO / MASS_NEUTRON * sum_percent;
}
// Calculate nuclide atom densities