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

Make sure correct direction is applied when fission neutrons are anisotropic

Browse source
This commit is contained in:
Paul Romano 2022-11-30 14:21:21 -06:00
parent 990b5449f5
commit 956e3f4341
2 changed files with 10 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/physics.cpp
View file

@ -212,7 +212,7 @@ void create_fission_sites(Particle& p, int i_nuclide, const Reaction& rx)
site.surf_id = 0;
// Sample delayed group and angle/energy for fission reaction
sample_fission_neutron(i_nuclide, rx, p.E(), &site, p.current_seed());
sample_fission_neutron(i_nuclide, rx, &site, p);
// Store fission site in bank
if (use_fission_bank) {
@ -1031,9 +1031,13 @@ Direction sample_cxs_target_velocity(
return vt * rotate_angle(u, mu, nullptr, seed);
}
void sample_fission_neutron(int i_nuclide, const Reaction& rx, double E_in,
SourceSite* site, uint64_t* seed)
void sample_fission_neutron(
int i_nuclide, const Reaction& rx, SourceSite* site, Particle& p)
{
// Get attributes of particle
double E_in = p.E();
uint64_t* seed = p.current_seed();
// Determine total nu, delayed nu, and delayed neutron fraction
const auto& nuc {data::nuclides[i_nuclide]};
double nu_t = nuc->nu(E_in, Nuclide::EmissionMode::total);
@ -1096,9 +1100,7 @@ void sample_fission_neutron(int i_nuclide, const Reaction& rx, double E_in,
}
// Sample azimuthal angle uniformly in [0, 2*pi) and assign angle
// TODO: account for dependence on incident neutron?
Direction ref(1., 0., 0.);
site->u = rotate_angle(ref, mu, nullptr, seed);
site->u = rotate_angle(p.u(), mu, nullptr, seed);
}
void inelastic_scatter(const Nuclide& nuc, const Reaction& rx, Particle& p)