fix some variable references in CO2_capture script

gmlc-dispatches · Dec 22, 2023 · 0fe948f · 0fe948f
1 parent b029d2e
commit 0fe948f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/dispatches/case_studies/fossil_case/ultra_supercritical_plant/storage/fe_ccs/CO2_capture.py b/dispatches/case_studies/fossil_case/ultra_supercritical_plant/storage/fe_ccs/CO2_capture.py
@@ -323,7 +323,7 @@ def SRD_eqn(b, t):
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.PZ_97_capture:
                     return b.SRD[t] == - 912.90813190561812007218 * CO2leanloading + 963.87626280230563224904 * exp(CO2leanloading) - 840.95461087585590576055 * (molality/2.)**2 - 614.38277546147764951456 * CO2leanloading**2 + 274.40532024830611135258 * (molality/2.)**3 - 0.61984641382653959951199 * (molality*CO2leanloading/2.)**3 + 0.26430328623580695568407 * (molality/CO2leanloading/2.) - 11.869291574687059309667 * (CO2leanloading/molality/0.5) - 0.23630953358941128236714E-002 * (molality/CO2leanloading/2.)**2
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.PZ_90to97_capture:
-                    return b.SRD[t] ==  - 0.10526756950236890175709 * Capture + 0.52057178809591142520929 * lean_loading + 0.98788205977359255793857 * molality + 2.0990364304677582296677 * Capture**0.5 + 3.5938248110847506033849 * lean_loading**0.5 - 6.3111812164050720141972 * molality**0.5 + 2.4029805693688151002618 * lean_loading**3 - 0.10937466124959506139080E-001 * Capture/molality - 41.496993880686794398116 * lean_loading/Capture + 0.36946342949766633467767E-001 * molality/lean_loading
+                    return b.SRD[t] ==  - 0.10526756950236890175709 * Capture + 0.52057178809591142520929 * CO2leanloading + 0.98788205977359255793857 * molality + 2.0990364304677582296677 * Capture**0.5 + 3.5938248110847506033849 * CO2leanloading**0.5 - 6.3111812164050720141972 * molality**0.5 + 2.4029805693688151002618 * CO2leanloading**3 - 0.10937466124959506139080E-001 * Capture/molality - 41.496993880686794398116 * CO2leanloading/Capture + 0.36946342949766633467767E-001 * molality/CO2leanloading
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.rel_PZAS_5mol_capture:
                     return (- 22.312295301354438947783 * CO2leanloading + 117.51384451486856619340 * CO2leanloading**2 - 269.05987870514132964672 * CO2leanloading**3 + 234.98153128275288281657 * CO2leanloading**4 + 3.9380879255782788028739)
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.PZ_ref1_90_capture:
@@ -365,7 +365,7 @@ def LL_flowrate(b, t):
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.PZ_97_capture:
                     return - 11073017307.372232437134 * CO2leanloading + 868882.20868137793149799 * log(CO2leanloading) - 73653243483.647430419922 * exp(molality/2.) + 11052450134.832271575928 * exp(CO2leanloading) + 141797181109.30328369141 * (molality/2.)**2 - 5395322021.3101339340210 * CO2leanloading**2 - 920354.31067807460203767 * (molality*CO2leanloading/2.) - 117144381.90464735031128 * (molality*CO2leanloading/2.)**3 - 6153045658.4519367218018 * (CO2leanloading/molality/0.5)**3
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.PZ_90to97_capture:
-                    return 21493.893918343124823878 * Capture - 18837259.090658042579889 * lean_loading - 303798.35064064717153087 * molality - 547220.78199295536614954 * Capture**0.5 + 7997911.0828296802937984 * lean_loading**0.5 + 1594229.8231671806424856 * molality**0.5 + 16593654.797202151268721 * lean_loading**2 + 58165.759382280644786078 * Capture*lean_loading + 60595.207259644026635215 * Capture/molality + 60446137.647265411913395 * (lean_loading/molality)**2
+                    return 21493.893918343124823878 * Capture - 18837259.090658042579889 * CO2leanloading - 303798.35064064717153087 * molality - 547220.78199295536614954 * Capture**0.5 + 7997911.0828296802937984 * CO2leanloading**0.5 + 1594229.8231671806424856 * molality**0.5 + 16593654.797202151268721 * CO2leanloading**2 + 58165.759382280644786078 * Capture*CO2leanloading + 60595.207259644026635215 * Capture/molality + 60446137.647265411913395 * (CO2leanloading/molality)**2
                 # No data for PZAS, using same as PZ 95 capture
                 elif b.config.CO2_capture_rate_surrogates == Surrogates.rel_PZAS_5mol_capture:
                     return - 8128346023.5960435867310 * CO2leanloading + 8097036304.4302654266357 * exp(CO2leanloading) + 143943.79958957055350766 * (molality/4.)**2 - 3906104299.9280681610107 * CO2leanloading**2 - 1694548695.5476403236389 * CO2leanloading**3 - 326827.12907762406393886 * (molality*CO2leanloading/4.)**2 + 14089575.472042093053460 * (CO2leanloading/molality/0.25) - 4257.3482804741406653193 * (molality/CO2leanloading/4.)**2 - 22772347.982399865984917 * (CO2leanloading/molality/0.25)**2 + 131.22409336241202026940 * (molality/CO2leanloading/4.)**3 + 22174561.087608262896538 * (CO2leanloading/molality/0.25)**3 - 8094525387.7033281326294