Add Electron Transport Capability

mcdc/constant.py

@@ -27,6 +27,7 @@
 SCORE_SPACE_MOMENT_CURRENT = 10
 SCORE_TIME_MOMENT_MU_SQ = 11
 SCORE_SPACE_MOMENT_MU_SQ = 12
+SCORE_EDEP = 13
 
 # Tally multipliers
 MULTIPLIER_ENERGY = 0
@@ -99,9 +100,17 @@
 REACTION_NEUTRON_FISSION = 4
 REACTION_NEUTRON_FISSION_PROMPT = 5
 REACTION_NEUTRON_FISSION_DELAYED = 6
+#
+REACTION_ELECTRON_ELASTIC_SCATTERING = 101
+REACTION_ELECTRON_ELASTIC_SMALL_ANGLE = 102
+REACTION_ELECTRON_ELASTIC_LARGE_ANGLE = 103
+REACTION_ELECTRON_IONIZATION = 104
+REACTION_ELECTRON_BREMSSTRAHLUNG = 105
+REACTION_ELECTRON_EXCITATION = 106
 
 # Particle types
 PARTICLE_NEUTRON = 0
+PARTICLE_ELECTRON = 1
 
 # Data
 DATA_NONE = 0
@@ -167,15 +176,34 @@
 AXIS_Z = 2
 AXIS_T = 3
 
+# ======================================================================================
 # Physics
+# ======================================================================================
+
+# Constants
 LIGHT_SPEED = 2.99792458e10  # cm/s
-NEUTRON_MASS = 939.565413e6  # eV/c^2
 BOLTZMANN_K = 8.61733326e-5  # eV/K
+
+# Adjustables
 THERMAL_THRESHOLD_FACTOR = 400
 
+# Neutron properties
+NEUTRON_MASS = 939.565413e6  # eV/c^2
+
+# Electron energy
+ELECTRON_MASS = 510.99895069e3  # eV/c^2
+ELECTRON_CUTOFF_ENERGY = 100  # eV (TODO: Is this adjustable?)
+FINE_STRUCTURE_CONSTANT = 7.2973525693e-3  # ~1/137
+MU_CUTOFF = 0.999999
+
+# ======================================================================================
+# Techniques
+# ======================================================================================
+
 # Weight Windows Methods
 WW_USER = 0
 WW_PREVIOUS = 1
+
 # Weight Windows Modifications
 WW_MIN = 0
 WW_WOLLABER = 1

mcdc/mcdc_get/__init__.py

@@ -44,6 +44,16 @@
 
 import mcdc.mcdc_get.tabulated_energy_angle_distribution as tabulated_energy_angle_distribution
 
+import mcdc.mcdc_get.element as element
+
+import mcdc.mcdc_get.electron_bremsstrahlung_reaction as electron_bremsstrahlung_reaction
+
+import mcdc.mcdc_get.electron_elastic_scattering_reaction as electron_elastic_scattering_reaction
+
+import mcdc.mcdc_get.electron_excitation_reaction as electron_excitation_reaction
+
+import mcdc.mcdc_get.electron_ionization_reaction as electron_ionization_reaction
+
 import mcdc.mcdc_get.gpu_meta as gpu_meta
 
 import mcdc.mcdc_get.native_material as native_material

mcdc/mcdc_get/electron_bremsstrahlung_reaction.py

@@ -0,0 +1,3 @@
+# The following is automatically generated by code_factory.py
+
+from numba import njit

mcdc/mcdc_get/electron_elastic_scattering_reaction.py

@@ -0,0 +1,3 @@
+# The following is automatically generated by code_factory.py
+
+from numba import njit

mcdc/mcdc_get/electron_excitation_reaction.py

@@ -0,0 +1,3 @@
+# The following is automatically generated by code_factory.py
+
+from numba import njit

mcdc/mcdc_get/electron_ionization_reaction.py

@@ -0,0 +1,32 @@
+# The following is automatically generated by code_factory.py
+
+from numba import njit
+
+
+@njit
+def binding_energy(index, electron_ionization_reaction, data):
+    offset = electron_ionization_reaction["binding_energy_offset"]
+    return data[offset + index]
+
+
+@njit
+def binding_energy_all(electron_ionization_reaction, data):
+    start = electron_ionization_reaction["binding_energy_offset"]
+    size = electron_ionization_reaction["binding_energy_length"]
+    end = start + size
+    return data[start:end]
+
+
+@njit
+def binding_energy_last(electron_ionization_reaction, data):
+    start = electron_ionization_reaction["binding_energy_offset"]
+    size = electron_ionization_reaction["binding_energy_length"]
+    end = start + size
+    return data[end - 1]
+
+
+@njit
+def binding_energy_chunk(start, length, electron_ionization_reaction, data):
+    start += electron_ionization_reaction["binding_energy_offset"]
+    end = start + length
+    return data[start:end]