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

tally photon heating with energy lost for analog estimator

Browse source
This commit is contained in:
liangjg 2019-03-27 15:26:42 -04:00
parent cfb11882a7
commit deda8ce076
1 changed files with 2 additions and 14 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/tallies/tally_scoring.cpp
View file

@ -1228,19 +1228,7 @@ score_general_ce(Particle* p, int i_tally, int start_index,
}
} else if (p->type_ == Particle::Type::photon) {
if (tally.estimator_ == ESTIMATOR_ANALOG) {
// Calculate material heating cross section
double macro_heating = 0.;
const Material& material {*model::materials[p->material_]};
for (auto i = 0; i < material.nuclide_.size(); ++i) {
auto i_element = material.element_[i];
auto atom_density = material.atom_density_(i);
auto& heating {data::elements[i_element].heating_};
auto i_grid = simulation::micro_photon_xs[i_element].index_grid;
auto f = simulation::micro_photon_xs[i_element].interp_factor;
macro_heating += std::exp(heating(i_grid) + f * (heating(i_grid+1) -
heating(i_grid))) * atom_density;
}
// All events score to an heating rate bin
// Use photon energy lost for heating deposition
if (settings::survival_biasing) {
// We need to account for the fact that some weight was already
// absorbed
@ -1248,7 +1236,7 @@ score_general_ce(Particle* p, int i_tally, int start_index,
} else {
score = p->wgt_last_;
}
score *= macro_heating * flux / p->macro_xs_.total;
score *= (E - p->E_) * flux;
} else {
// Calculate photon heating cross section on-the-fly
score = 0.;