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

Enable the LegendreFilter filter to be used in photon tallies for orders greater than P0. (#3245)

Browse source 
Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Adam Nelson 2025-01-14 09:51:47 -06:00 • committed by GitHub
parent d39a414011
commit 549cc0973c
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
1 changed files with 7 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/physics.cpp
View file

@ -293,8 +293,8 @@ void sample_photon_reaction(Particle& p)
// Coherent (Rayleigh) scattering
prob += micro.coherent;
if (prob > cutoff) {
double mu = element.rayleigh_scatter(alpha, p.current_seed());
p.u() = rotate_angle(p.u(), mu, nullptr, p.current_seed());
p.mu() = element.rayleigh_scatter(alpha, p.current_seed());
p.u() = rotate_angle(p.u(), p.mu(), nullptr, p.current_seed());
p.event() = TallyEvent::SCATTER;
p.event_mt() = COHERENT;
return;
@ -303,10 +303,10 @@ void sample_photon_reaction(Particle& p)
// Incoherent (Compton) scattering
prob += micro.incoherent;
if (prob > cutoff) {
double alpha_out, mu;
double alpha_out;
int i_shell;
element.compton_scatter(
alpha, true, &alpha_out, &mu, &i_shell, p.current_seed());
alpha, true, &alpha_out, &p.mu(), &i_shell, p.current_seed());
// Determine binding energy of shell. The binding energy is 0.0 if
// doppler broadening is not used.
@ -322,9 +322,9 @@ void sample_photon_reaction(Particle& p)
double E_electron = (alpha - alpha_out) * MASS_ELECTRON_EV - e_b;
int electron = static_cast<int>(ParticleType::electron);
if (E_electron >= settings::energy_cutoff[electron]) {
double mu_electron = (alpha - alpha_out * mu) /
double mu_electron = (alpha - alpha_out * p.mu()) /
std::sqrt(alpha * alpha + alpha_out * alpha_out -
2.0 * alpha * alpha_out * mu);
2.0 * alpha * alpha_out * p.mu());
Direction u = rotate_angle(p.u(), mu_electron, &phi, p.current_seed());
p.create_secondary(p.wgt(), u, E_electron, ParticleType::electron);
}
@ -338,7 +338,7 @@ void sample_photon_reaction(Particle& p)
phi += PI;
p.E() = alpha_out * MASS_ELECTRON_EV;
p.u() = rotate_angle(p.u(), mu, &phi, p.current_seed());
p.u() = rotate_angle(p.u(), p.mu(), &phi, p.current_seed());
p.event() = TallyEvent::SCATTER;
p.event_mt() = INCOHERENT;
return;