Added parametric fuel pin size and driving script for 37-element CANDU bundle

1 files changed, 25 insertions, 0 deletions

diff --git a/driver-37.sh b/driver-37.sh
new file mode 100755
index 0000000..31a81f5
--- /dev/null
+++ b/driver-37.sh
@@ -0,0 +1,25 @@
+#!/bin/bash
+
+# Driver script for parametric pin size analysis of the 
+# 37-element CANDU fuel bundle. AnnCon 2026 paper.
+# Connor Moore, 2026. <connor.moore@ontariotechu.net>
+
+green="\033[32m"
+reset="\033[0m"
+
+mkdir -p 37-results/
+touch 37-results/37_keff_table
+
+for r in $(seq 0.1 0.05 0.8); do
+    echo -e ${green}Running rf = $r cm... ${reset}
+    filename=37-results/$r.log
+    python3 candu-37.py -rf $r > $filename
+    keff=$(cat $filename | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
+    pm=$(cat $filename | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
+    flow_area=$(cat $filename | grep "Flow area" | cut -d "=" -f 2)
+    fuel_mass=$(cat $filename | grep "Fuel mass" | cut -d "=" -f 2)
+    clad_mass=$(cat $filename | grep "Cladding mass" | cut -d "=" -f 2)
+    dtu=$(cat $filename | grep "DTU ratio" | cut -d "=" -f 2)
+
+    echo $keff +/- $pm, $flow_area cm² flow, $fuel_mass g/cm fuel, $clad_mass g/cm clad, $dtu DTU ratio | cowsay
+done