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

Skip atomic relaxation if binding energy is larger than photon energy (#3391)

Browse source 
Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Amanda Lund 2025-05-03 00:46:14 -05:00 • committed by GitHub
parent a921280fa9
commit 512df2f4ff
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
2 changed files with 32 additions and 31 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/photon.cpp
View file

@ -770,8 +770,9 @@ void PhotonInteraction::pair_production(double alpha, double* E_electron,
void PhotonInteraction::atomic_relaxation(int i_shell, Particle& p) const
{
// Return if no atomic relaxation data is present
if (!has_atomic_relaxation_)
// Return if no atomic relaxation data is present or if the binding energy is
// larger than the incident particle energy
if (!has_atomic_relaxation_ || shells_[i_shell].binding_energy > p.E())
return;
// Stack for unprocessed holes left by transitioning electrons

58
tests/regression_tests/photon_production_fission/results_true.dat
View file

@ -1,12 +1,12 @@
k-combined:
2.272421E+00 4.418825E-02
2.270911E+00 4.568134E-02
tally 1:
2.665301E+00
2.369839E+00
2.663476E+00
2.366726E+00
0.000000E+00
0.000000E+00
2.665301E+00
2.369839E+00
2.663476E+00
2.366726E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -18,52 +18,52 @@ tally 1:
0.000000E+00
0.000000E+00
tally 2:
2.637532E+00
2.320051E+00
4.231377E+08
5.971498E+16
2.636709E+00
2.318646E+00
4.230052E+08
5.967879E+16
0.000000E+00
0.000000E+00
2.637532E+00
2.320051E+00
4.231377E+08
5.971498E+16
2.636709E+00
2.318646E+00
4.230052E+08
5.967879E+16
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.448745E+06
2.005297E+12
2.455009E+06
2.015314E+12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.448745E+06
2.005297E+12
2.455009E+06
2.015314E+12
0.000000E+00
0.000000E+00
tally 3:
2.660000E+00
2.358558E+00
4.231377E+08
5.971498E+16
2.660004E+00
2.358564E+00
4.230052E+08
5.967879E+16
0.000000E+00
0.000000E+00
2.660000E+00
2.358558E+00
4.231377E+08
5.971498E+16
2.660004E+00
2.358564E+00
4.230052E+08
5.967879E+16
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.448745E+06
2.005297E+12
2.455009E+06
2.015314E+12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.448745E+06
2.005297E+12
2.455009E+06
2.015314E+12
0.000000E+00
0.000000E+00