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

Change outgoing energy adjustment for incoherent inelastic to match MCNP6

Browse source
This commit is contained in:
Paul Romano 2019-11-27 11:11:01 -06:00
parent d3251e62f8
commit eb12d612e9
1 changed files with 7 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

13
src/secondary_thermal.cpp
View file

@ -252,8 +252,8 @@ IncoherentInelasticAE::sample(double E_in, double& E_out, double& mu) const
double f;
get_energy_index(energy_, E_in, i, f);
// Sample between ith and [i+1]th bin
int l = f > prn() ? i + 1 : i;
// Pick closer energy based on interpolation factor
int l = f > 0.5 ? i + 1 : i;
// Determine endpoints on grid i
auto n = distribution_[i].e_out.size();
@ -301,11 +301,12 @@ IncoherentInelasticAE::sample(double E_in, double& E_out, double& mu) const
2.0*frac*(r1 - c_j))) - p_l_j) / frac;
}
// Now interpolate between incident energy bins i and i + 1
if (l == i) {
E_out = E_1 + (E_out - E_i_1) * (E_J - E_1) / (E_i_J - E_i_1);
// Adjustment of outgoing energy
double E_l = energy_[l];
if (E_out < 0.5*E_l) {
E_out *= 2.0*E_in/E_l - 1.0;
} else {
E_out = E_1 + (E_out - E_i1_1) * (E_J - E_1) / (E_i1_J - E_i1_1);
E_out += E_in - E_l;
}
// Sample outgoing cosine bin