Update fusion cross-section to Bosch-Hale for PRD (#101)

* Switch from Hively to Bosch-Hale fusion cross-section * Update notebooks * Use all significant figures from Bosch-Hale * Fix notebook error
cfs-energy · Sep 11, 2024 · 6825b2d · 6825b2d
1 parent ecbc626
commit 6825b2d
Show file tree

Hide file tree

Showing 8 changed files with 90 additions and 85 deletions.
diff --git a/cfspopcon/formulas/fusion_power/fusion_data.py b/cfspopcon/formulas/fusion_power/fusion_data.py
@@ -34,15 +34,15 @@ def __init__(self) -> None:
     def calc_rate_coefficient(ion_temp: float) -> float:
         r"""Calculate :math:`\\langle \\sigma v \rangle` for a given ion temperature.
 
-        Cross-section from :cite:`bosch_improved_1992`
+        Cross-section from :cite:`bosch_improved_1992`, equation 12, using coefficients from table VII, page 625 (first column is DT reaction)
 
         Args:
             ion_temp: [keV]
 
         Returns:
             :math:`\\langle \\sigma v \rangle` [cm^3/s]
         """
-        C = [0.0, 1.173e-9, 1.514e-2, 7.519e-2, 4.606e-3, 1.35e-2, -1.068e-4, 1.366e-5]
+        C = [0.0, 1.17302e-9, 1.51361e-2, 7.51886e-2, 4.60643e-3, 1.35000e-2, -1.06750e-4, 1.36600e-5]
         B_G = 34.3827
         mr_c2 = 1124656
 

diff --git a/docs/doc_sources/getting_started.ipynb b/docs/doc_sources/getting_started.ipynb
diff --git a/docs/doc_sources/separatrix_operational_space.ipynb b/docs/doc_sources/separatrix_operational_space.ipynb
diff --git a/docs/doc_sources/time_independent_inductances_and_fluxes.ipynb b/docs/doc_sources/time_independent_inductances_and_fluxes.ipynb
diff --git a/docs/doc_sources/understanding_algorithms.ipynb b/docs/doc_sources/understanding_algorithms.ipynb
diff --git a/example_cases/SPARC_PRD/input.yaml b/example_cases/SPARC_PRD/input.yaml
@@ -168,7 +168,7 @@ plasma_current: 8.7e+6
 fraction_of_external_power_coupled: 0.9
 
 # What fusion reaction are we using?
-fusion_reaction: DTFusionHively
+fusion_reaction: DTFusionBoschHale
 # Fraction of fuel ions which are the heavier species
 heavier_fuel_species_fraction: 0.5
 

diff --git a/tests/regression_results/PRD.json b/tests/regression_results/PRD.json
@@ -65,35 +65,35 @@
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 22.55867182,
+            "data": 23.15670789,
             "dims": []
         },
         "P_auxiliary_absorbed": {
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 7.702542635,
+            "data": 7.104506562,
             "dims": []
         },
         "P_auxiliary_launched": {
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 8.558380706,
+            "data": 7.893896180,
             "dims": []
         },
         "P_external": {
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 8.904164506,
+            "data": 8.306128433,
             "dims": []
         },
         "P_fusion": {
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 112.7933591,
+            "data": 115.7835394,
             "dims": []
         },
         "P_in": {
@@ -107,7 +107,7 @@
             "attrs": {
                 "units": "megawatt"
             },
-            "data": 90.23468727,
+            "data": 92.62683156,
             "dims": []
         },
         "P_ohmic": {
@@ -135,7 +135,7 @@
             "attrs": {
                 "units": "dimensionless"
             },
-            "data": 11.55669358,
+            "data": 12.72973555,
             "dims": []
         },
         "SOL_momentum_loss_function": {
@@ -320,14 +320,14 @@
             "attrs": {
                 "units": "dimensionless"
             },
-            "data": 0.08439528049,
+            "data": 0.08439527661,
             "dims": []
         },
         "edge_impurity_concentration_in_core": {
             "attrs": {
                 "units": "dimensionless"
             },
-            "data": 0.01687905610,
+            "data": 0.01687905532,
             "dims": []
         },
         "edge_impurity_enrichment": {
@@ -634,7 +634,7 @@
         },
         "fusion_reaction": {
             "attrs": {},
-            "data": "DTFusionHively",
+            "data": "DTFusionBoschHale",
             "dims": []
         },
         "fusion_triple_product": {
@@ -687,7 +687,7 @@
             "data": [
                 0.06000000000,
                 0.003100000000,
-                0.01687905610,
+                0.01687905532,
                 0.000000000,
                 1.500000000e-05
             ],
@@ -933,14 +933,14 @@
             "attrs": {
                 "units": "megawatt / meter ** 2"
             },
-            "data": 1.624788752,
+            "data": 1.667862312,
             "dims": []
         },
         "neutron_rate": {
             "attrs": {
                 "units": "1 / second"
             },
-            "data": 4.000004399e+19,
+            "data": 4.106045524e+19,
             "dims": []
         },
         "normalized_beta": {

diff --git a/tests/regression_results/SPARC_PRD_result.nc b/tests/regression_results/SPARC_PRD_result.nc