Merge pull request #2118 from NREL/sco2_py_samples

Update sco2 python functionality for new configurations

Author: taylorbrown75

.gitignore

@@ -391,3 +391,4 @@ deploy/mslf_iph_mimo_ssc_auto_exec.lk
 deploy/physical_trough_v2_mimo_ssc_auto_exec.lk
 
 deploy/physical_trough_IPH_v2_mimo_ssc_auto_exec.lk
+/samples/CSP/sco2_analysis_python_V2/.vscode

samples/CSP/sco2_analysis_python_V2/changelog.md

@@ -0,0 +1,31 @@
+# Changelog
+
+All notable changes to this project will be documented in this file.
+
+## [2025.10.21]
+### Git Commits
+- **SSC**: https://github.com/NREL/ssc/commit/1e7fc7d740bd1dbd56a578b0323af96da70a525f
+- **SAM**: https://github.com/NREL/SAM/commit/82466c3527505f0adc59bf9e8a140d3c08ed593a
+### Added
+- Two new cycle configurations:
+  - Recompression with HTR Bypass
+  - Turbine Split Flow
+- Inflation year parameter (`yr_inflation`)
+  - Modifies cycle component costs to target inflation year
+- Thermal efficiency cutoff
+  - Prevents finalization of component designs if efficiency is below the threshold
+
+### Changed
+- Recuperator cost model
+  - Adds temperature-based cost factor for temperatures above 550 °C
+  - Based on Weiland et al., 2019
+- Axial turbine cost model
+  - Adds temperature-based cost factor for inlet temperatures above 550 °C
+  - Based on Weiland et al., 2019
+- Optimization objective metric calculation
+  - Improved behavior when using a minimum temperature difference (`min_phx_deltaT`)
+
+### Fixed
+- Improved stability of air cooler convergence at low cycle efficiencies
+- Fixed bug solving heat exchangers with low effectiveness due to pinch points
+- Fixed heat exchanger bug with 2-phase cold fluid stream

samples/CSP/sco2_analysis_python_V2/core/sco2_cycle_ssc.py

@@ -351,6 +351,10 @@ def get_cycle_name_str(cycle_des_par_dict):
         return "simple"
     elif(cycle_config == 2):
         return "partialcooling"
+    elif(cycle_config == 3):
+        return "htrbp"
+    elif(cycle_config==4):
+        return "turbinesplitflow"
     else:
         return "unknowncycle"
 
@@ -370,6 +374,14 @@ def mod_dict_for_cycle_config(dict_in, cycle_config_str):
         dict_out["cycle_config"] = 2
         dict_out["is_recomp_ok"] = 1
 
+    elif(cycle_config_str == "htrbp"):
+        dict_out["cycle_config"] = 3
+        dict_out["is_recomp_ok"] = 1
+
+    elif(cycle_config_str == "turbinesplitflow"):
+        dict_out["cycle_config"] = 4
+        dict_out["is_recomp_ok"] = 0
+
     else:
         dict_out["cycle_config"] = -999
         dict_out["is_recomp_ok"] = -999
@@ -445,7 +457,18 @@ def get_default_sco2_dict():
     des_par["fan_power_frac"] = 0.02  # [-] Fraction of net cycle power consumed by air cooler fan. 2% here per Turchi et al.
     # Default
     des_par["deltaP_counterHX_frac"] = 0.0054321  # [-] Fraction of CO2 inlet pressure that is design point counterflow HX (recups & PHX) pressure drop
-	
+    
+    # HTR BP parameters
+    des_par['is_bypass_ok'] = -0.11      # 1 = Optimize, 0 = no bypass, < 0 = fix bp_frac to abs(input)
+    des_par['T_bypass_target'] = 480 # [C] Target HTF Outlet Temperature (if bypass frac is not fixed)
+    des_par['T_target_is_HTF'] = 1   # Temperature target is HTF outlet (rather than sco2)
+    des_par['deltaT_bypass'] = 0     # [dC] Temperature difference at Mixer 2
+    des_par['set_HTF_mdot'] = 0      # [kg/s or False] Do NOT set HTF mdot (model solves it with approach temps)
+
+    # TSF parameters
+    des_par['is_turbine_split_ok'] = 1   # 1 = Optimize, 0 = no bypass, < 0 = fix bp_frac to abs(input)
+    des_par['eta_isen_t2'] = des_par['eta_isen_t']       # Turbine 2 isentropic efficiency
+
     return des_par
 
 class C_sco2_sim: