Initial activation comparisons for ITER DT flux spectra

iter/README.md

@@ -0,0 +1,19 @@
+# README for ITER Activation Comparisons
+
+This directory runs activation calculations for Be and W, two proposed ITER first wall materials, with ITER flux spectra provided by FISPACT [[1]](#1).
+
+These calculations are broken up into a series of 2, 4, 8, 32, and 64 pulses, with duty cycles of 100% (steady-state), 90%, 50%, and 25%. These duty cycles
+have the effect of increasing the dwell time between successive pulses, while the pulse length is only a function of the number of pulses and the total burn time.
+For these calculations, the total active burn time is 4 years. The table below summarizes the dwell times as a function of duty cycle and total number of pulses.
+
+| # Pulses | 90% | 50% | 25% |
+| ---------|-----|-----|-----|
+|  2 | 0.222 | 2 | 6 |
+|  4 | 0.111 | 1 | 3 |
+|  8 | 0.056 | 0.5 | 1.5 |
+| 32 | 0.01389 | 0.125 | 0.375 |
+| 64 | 0.006944 | 0.0625 | 0.1875 |
+
+## References
+<a id="1">[1]</a>
+https://fispact.ukaea.uk/wiki/Reference_input_spectra (accessed Nov. 24, 2025).

iter/duty_cycle_100/iter_dt_ss_2_pulses

@@ -0,0 +1,38 @@
+# This problem runs activation on two proposed ITER FW materials (Be and W)
+# Steady-state operation with no dwell time is used for this model, which will be used as the basis for comparison
+geometry rectangular
+volume
+     1    Be
+     1    W
+end
+mat_loading
+    Be    mix_1
+    W    mix_2
+end
+mixture mix_1
+    element be 1 1.0
+end
+mixture mix_2
+    element w 1 1.0
+end
+#Use default ALARA material and element libraries
+material_lib /root/data/matlib.sample
+element_lib /root/data/elelib.std
+data_library alaralib /opt/groupspace/shared/n/nukecode/ALARA/data/fendl2bin
+output zone
+ 	units Bq cm3
+        number_density
+ 	specific_activity
+end
+#     flux name    fluxin file   norm   shift   unused
+flux  my_flux       ../../ref_flux_files/iter_dt_flux_2.0986E14 1     0      default
+# Run for 4 FPY with 2 pulses
+schedule    my_schedule
+    2 y my_flux ph_1 0  s
+end
+pulsehistory  ph_1
+    2    0.0    s
+end 
+#other parameters
+truncation 1e-5
+dump_file dump_iter_dt_ss_2_pulses.file

iter/duty_cycle_100/iter_dt_ss_32_pulses

@@ -0,0 +1,38 @@
+# This problem runs activation on two proposed ITER FW materials (Be and W)
+# Steady-state operation with no dwell time is used for this model, which will be used as the basis for comparison
+geometry rectangular
+volume
+     1    Be
+     1    W
+end
+mat_loading
+    Be    mix_1
+    W    mix_2
+end
+mixture mix_1
+    element be 1 1.0
+end
+mixture mix_2
+    element w 1 1.0
+end
+#Use default ALARA material and element libraries
+material_lib /root/data/matlib.sample
+element_lib /root/data/elelib.std
+data_library alaralib /opt/groupspace/shared/n/nukecode/ALARA/data/fendl2bin
+output zone
+ 	units Bq cm3
+        number_density
+ 	specific_activity
+end
+#     flux name    fluxin file   norm   shift   unused
+flux  my_flux       ../../ref_flux_files/iter_dt_flux_2.0986E14 1     0      default
+# Run for 4 FPY with 32 pulses
+schedule    my_schedule
+    0.125 y my_flux ph_1 0  s
+end
+pulsehistory  ph_1
+    32    0.0    s
+end 
+#other parameters
+truncation 1e-5
+dump_file dump_iter_dt_ss_32_pulses.file

iter/duty_cycle_100/iter_dt_ss_4_pulses

@@ -0,0 +1,38 @@
+# This problem runs activation on two proposed ITER FW materials (Be and W)
+# Steady-state operation with no dwell time is used for this model, which will be used as the basis for comparison
+geometry rectangular
+volume
+     1    Be
+     1    W
+end
+mat_loading
+    Be    mix_1
+    W    mix_2
+end
+mixture mix_1
+    element be 1 1.0
+end
+mixture mix_2
+    element w 1 1.0
+end
+#Use default ALARA material and element libraries
+material_lib /root/data/matlib.sample
+element_lib /root/data/elelib.std
+data_library alaralib /opt/groupspace/shared/n/nukecode/ALARA/data/fendl2bin
+output zone
+ 	units Bq cm3
+        number_density
+ 	specific_activity
+end
+#     flux name    fluxin file   norm   shift   unused
+flux  my_flux       ../../ref_flux_files/iter_dt_flux_2.0986E14 1     0      default
+# Run for 4 FPY with 4 pulses
+schedule    my_schedule
+    1 y my_flux ph_1 0  s
+end
+pulsehistory  ph_1
+    4    0.0    s
+end 
+#other parameters
+truncation 1e-5
+dump_file dump_iter_dt_ss_4_pulses.file

iter/duty_cycle_100/iter_dt_ss_64_pulses

@@ -0,0 +1,38 @@
+# This problem runs activation on two proposed ITER FW materials (Be and W)
+# Steady-state operation with no dwell time is used for this model, which will be used as the basis for comparison
+geometry rectangular
+volume
+     1    Be
+     1    W
+end
+mat_loading
+    Be    mix_1
+    W    mix_2
+end
+mixture mix_1
+    element be 1 1.0
+end
+mixture mix_2
+    element w 1 1.0
+end
+#Use default ALARA material and element libraries
+material_lib /root/data/matlib.sample
+element_lib /root/data/elelib.std
+data_library alaralib /opt/groupspace/shared/n/nukecode/ALARA/data/fendl2bin
+output zone
+ 	units Bq cm3
+        number_density
+ 	specific_activity
+end
+#     flux name    fluxin file   norm   shift   unused
+flux  my_flux       ../../ref_flux_files/iter_dt_flux_2.0986E14 1     0      default
+# Run for 4 FPY with 64 pulses
+schedule    my_schedule
+    0.0625 y my_flux ph_1 0  s
+end
+pulsehistory  ph_1
+    64    0.0    s
+end 
+#other parameters
+truncation 1e-5
+dump_file dump_iter_dt_ss_64_pulses.file

iter/duty_cycle_100/iter_dt_ss_8_pulses

@@ -0,0 +1,38 @@
+# This problem runs activation on two proposed ITER FW materials (Be and W)
+# Steady-state operation with no dwell time is used for this model, which will be used as the basis for comparison
+geometry rectangular
+volume
+     1    Be
+     1    W
+end
+mat_loading
+    Be    mix_1
+    W    mix_2
+end
+mixture mix_1
+    element be 1 1.0
+end
+mixture mix_2
+    element w 1 1.0
+end
+#Use default ALARA material and element libraries
+material_lib /root/data/matlib.sample
+element_lib /root/data/elelib.std
+data_library alaralib /opt/groupspace/shared/n/nukecode/ALARA/data/fendl2bin
+output zone
+ 	units Bq cm3
+        number_density
+ 	specific_activity
+end
+#     flux name    fluxin file   norm   shift   unused
+flux  my_flux       ../../ref_flux_files/iter_dt_flux_2.0986E14 1     0      default
+# Run for 4 FPY with 8 pulses
+schedule    my_schedule
+    0.5 y my_flux ph_1 0  s
+end
+pulsehistory  ph_1
+    8    0.0    s
+end 
+#other parameters
+truncation 1e-5
+dump_file dump_iter_dt_ss_8_pulses.file