Merge pull request #1353 from NREL/sco2_tsf

Add two new sco2 configurations

Author: taylorbrown75

ssc/cmod_csp_heatsink.cpp

@@ -60,7 +60,7 @@ class cm_csp_heatsink : public compute_module
 
             NS_HX_counterflow_eqs::E_UA_target_type ua_type = NS_HX_counterflow_eqs::E_UA_target_type::E_constant_UA;
 
-            mc_phx.initialize(hot_fl, 10, ua_type);
+            mc_phx.initialize(hot_fl, 10, ua_type, 2024.0);
 
             double q_dot = 1000.0;      //[kWt]
             double T_co2_hot = 700.0;   //[C]

ssc/cmod_sco2_csp_system.cpp

@@ -38,6 +38,7 @@ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 #include <ctime>
 
 #include "sco2_pc_csp_int.h"
+#include "sco2_htrbypass_cycle.h"
 
 static var_info _cm_vtab_sco2_csp_system[] = {
 
@@ -363,7 +364,7 @@ class cm_sco2_csp_system : public compute_module
 		bool is_od_set_control = is_assigned("od_set_control");
         bool is_od_generate_udpc_assigned = is_assigned("od_generate_udpc");
 		if (is_od_cases_assigned && is_P_mc_in_od_sweep_assigned)
-		{
+		{ 
 			log("Both off design cases and main compressor inlet pressure sweep assigned. Only modeling off design cases");
 			is_P_mc_in_od_sweep_assigned = false;
 		}
@@ -452,11 +453,18 @@ class cm_sco2_csp_system : public compute_module
 			P_LP_comp_in_des = c_sco2_cycle.get_design_solved()->ms_rc_cycle_solved.m_pres[C_sco2_cycle_core::MC_IN] / 1000.0;		//[MPa] convert from kPa
 			delta_P = 10.0;
 		}
-		else
+		else if (cycle_config == 2)
 		{
 			P_LP_comp_in_des = c_sco2_cycle.get_design_solved()->ms_rc_cycle_solved.m_pres[C_sco2_cycle_core::PC_IN] / 1000.0;		//[MPa] convert from kPa
 			delta_P = 6.0;
 		}
+        else
+        {
+            log("Cycle configuration not supported for off-design");
+            std::string err_msg = util::format("Cycle configuration not supported for off-design");
+            throw exec_error("sco2_csp_system", err_msg);
+            return;
+        }
 
         // Get turbine inlet mode
         int T_t_in_mode = as_integer("od_T_t_in_mode");
@@ -1741,3 +1749,7 @@ class cm_sco2_comp_curves : public compute_module
 };
 
 DEFINE_MODULE_ENTRY(sco2_comp_curves, "Calls sCO2 auto-design cycle function", 1)
+
+
+
+

ssc/csp_common.cpp

@@ -83,8 +83,6 @@ extern var_info vtab_sco2_design[];
 
 int sco2_design_cmod_common(compute_module *cm, C_sco2_phx_air_cooler & c_sco2_cycle);
 
+#endif
 
 
-
-
-#endif

ssc/csp_common.h

@@ -52,10 +52,12 @@ set(TCS_SRC
 		numeric_solvers.cpp
         ptes_solver_design_point.cpp
 		sco2_cycle_components.cpp
+        sco2_htrbypass_cycle.cpp
 		sco2_partialcooling_cycle.cpp
 		sco2_pc_csp_int.cpp
 		sco2_power_cycle.cpp
 		sco2_recompression_cycle.cpp
+        sco2_turbinesplitflow_cycle.cpp
 		tcskernel.cpp
 		trnsys_weatherreader.cpp
 		typelib.cpp
@@ -131,10 +133,12 @@ set(TCS_SRC
 		sam_csp_util.h
 		sco2_cycle_components.h
 		sco2_cycle_templates.h
+        sco2_htrbypass_cycle.h
 		sco2_partialcooling_cycle.h
 		sco2_pc_csp_int.h
 		sco2_power_cycle.h
 		sco2_recompression_cycle.h
+        sco2_turbinesplitflow_cycle.h
 		storage_hx.h
 		tcskernel.h
 		tcstype.h

tcs/CMakeLists.txt

@@ -35,6 +35,7 @@ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 #include "sam_csp_util.h"
 #include <algorithm>
 #include "numeric_solvers.h"
+#include "sco2_cycle_components.h"
 
 double NS_HX_counterflow_eqs::calc_max_q_dot_enth(int hot_fl_code /*-*/, HTFProperties & hot_htf_class,
     int cold_fl_code /*-*/, HTFProperties & cold_htf_class,
@@ -1523,11 +1524,13 @@ void NS_HX_counterflow_eqs::solve_q_dot_for_fixed_UA_enth(int hot_fl_code /*-*/,
 			throw(C_csp_exception("Off-design heat exchanger method failed"));
 		}
 
-		//if (!(od_hx_code == C_monotonic_eq_solver::CONVERGED || od_hx_code == C_monotonic_eq_solver::SLOPE_POS_NO_POS_ERR || od_hx_code == C_monotonic_eq_solver::SLOPE_POS_BOTH_ERRS))
-		//{
-		//	throw(C_csp_exception("Off-design heat exchanger method failed"));
-		//}
+		/*if (!(od_hx_code == C_monotonic_eq_solver::CONVERGED || od_hx_code == C_monotonic_eq_solver::SLOPE_POS_NO_POS_ERR || od_hx_code == C_monotonic_eq_solver::SLOPE_POS_BOTH_ERRS))
+		{
+			throw(C_csp_exception("Off-design heat exchanger method failed"));
+		}*/
+
 	}
+
 	else if (test_code == 0 && diff_UA_max_eff <= 0.0)  // UA_max_eff <= UA_target)
 	{
 		// At maximum allowable heat transfer, the calculated UA is less than target
@@ -1785,16 +1788,37 @@ double C_HX_counterflow_CRM::calc_max_q_dot_enth(double h_h_in /*kJ/kg*/, double
 
 double /*M$*/ C_HX_counterflow_CRM::calculate_equipment_cost(double UA /*kWt/K*/,
     double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/,
-    double T_cold_in /*K*/, double P_cold_in /*kPa*/, double m_dot_cold /*kg/s*/)
+    double T_cold_in /*K*/, double P_cold_in /*kPa*/, double m_dot_cold /*kg/s*/,
+    double yr_inflation /*yr*/)
 {
     switch (m_cost_model)
     {
     case C_HX_counterflow_CRM::E_CARLSON_17_RECUP:
-        return 1.25*1.E-3*UA;		//[M$] needs UA in kWt/K
+    {
+        double yr_base_inflation = 2017;
+        double f_inflation = calculate_inflation_factor(yr_base_inflation, yr_inflation);
+        return 1.25 * 1.E-3 * UA * f_inflation;		//[M$] needs UA in kWt/K
+    }  
     case C_HX_counterflow_CRM::E_WEILAND_19_RECUP:
-        return 49.45*std::pow(UA*1.E3, 0.7544)*1.E-6;  //[M$] needs UA in Wt/K
+    {
+        double C_recup = 49.45 * std::pow(UA * 1.E3, 0.7544) * 1.E-6; //[M$] needs UA in Wt/K
+        double T_factor = 1;
+        double T_max_C = std::max(T_hot_in, T_cold_in) - 273.15;
+        if (T_max_C >= 550)
+        {
+            T_factor = 1.0 + 0.02141 * (T_max_C - 550);
+        }
+        double yr_base_inflation = 2017;
+        double f_inflation = calculate_inflation_factor(yr_base_inflation, yr_inflation);
+        return C_recup * T_factor * f_inflation;  //[M$]
+        break;
+    }
     case C_HX_counterflow_CRM::E_CARLSON_17_PHX:
-        return 3.5*1.E-3*UA;		//[M$] needs UA in kWt/K
+    {
+        double yr_base_inflation = 2017;
+        double f_inflation = calculate_inflation_factor(yr_base_inflation, yr_inflation);
+        return 3.5 * 1.E-3 * UA * f_inflation;		//[M$] needs UA in kWt/K
+    } 
     default:
         return std::numeric_limits<double>::quiet_NaN();
     }
@@ -1860,7 +1884,8 @@ void C_HX_counterflow_CRM::design_calc_UA(C_HX_counterflow_CRM::S_des_calc_UA_pa
 
 	ms_des_solved.m_cost_equipment = calculate_equipment_cost(ms_des_solved.m_UA_design,
 		ms_des_calc_UA_par.m_T_h_in, ms_des_calc_UA_par.m_P_h_in, ms_des_calc_UA_par.m_m_dot_hot_des,
-		ms_des_calc_UA_par.m_T_c_in, ms_des_calc_UA_par.m_P_c_in, ms_des_calc_UA_par.m_m_dot_cold_des);
+		ms_des_calc_UA_par.m_T_c_in, ms_des_calc_UA_par.m_P_c_in, ms_des_calc_UA_par.m_m_dot_cold_des,
+        m_yr_inflation);
 
     ms_des_solved.m_cost_bare_erected = calculate_bare_erected_cost(ms_des_solved.m_cost_equipment);
 
@@ -1981,7 +2006,7 @@ void C_HX_counterflow_CRM::design_calc_UA_TP_to_PH(C_HX_counterflow_CRM::S_des_c
 
     ms_des_solved.m_cost_equipment = calculate_equipment_cost(ms_des_solved.m_UA_design,
         ms_des_calc_UA_par.m_T_h_in, ms_des_calc_UA_par.m_P_h_in, ms_des_calc_UA_par.m_m_dot_hot_des,
-        ms_des_calc_UA_par.m_T_c_in, ms_des_calc_UA_par.m_P_c_in, ms_des_calc_UA_par.m_m_dot_cold_des);
+        ms_des_calc_UA_par.m_T_c_in, ms_des_calc_UA_par.m_P_c_in, ms_des_calc_UA_par.m_m_dot_cold_des, m_yr_inflation);
 
     ms_des_solved.m_cost_bare_erected = calculate_bare_erected_cost(ms_des_solved.m_cost_equipment);
 
@@ -2059,7 +2084,8 @@ void C_HX_counterflow_CRM::design_for_target__calc_outlet(int hx_target_code /*-
 
     ms_des_solved.m_cost_equipment = calculate_equipment_cost(ms_des_solved.m_UA_design,
         T_h_in, P_h_in, m_dot_h,
-        T_c_in, P_c_in, m_dot_c);
+        T_c_in, P_c_in, m_dot_c,
+        m_yr_inflation);
 
     ms_des_solved.m_cost_bare_erected = calculate_bare_erected_cost(ms_des_solved.m_cost_equipment);
 }
@@ -3370,7 +3396,8 @@ double NS_HX_counterflow_eqs::UA_scale_vs_m_dot(double m_dot_cold_over_des /*-*/
     return pow(m_dot_ratio, 0.8);
 }
 
-void C_HX_co2_to_co2_CRM::initialize(int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type)
+void C_HX_co2_to_co2_CRM::initialize(int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type,
+                                     double yr_inflation)
 {
     // Set design parameters member structure
     ms_init_par.m_N_sub_hx = N_sub_hx;
@@ -3379,9 +3406,11 @@ void C_HX_co2_to_co2_CRM::initialize(int N_sub_hx, NS_HX_counterflow_eqs::E_UA_t
     m_is_HX_initialized = true;
 
     m_od_UA_target_type = od_UA_target_type;
+    m_yr_inflation = yr_inflation;
 }
 
-void C_HX_co2_to_htf::initialize(int hot_fl, util::matrix_t<double> hot_fl_props, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type)
+void C_HX_co2_to_htf::initialize(int hot_fl, util::matrix_t<double> hot_fl_props, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type,
+                                 double yr_inflation)
 {
 	// Hard-code some of the design parameters
     ms_init_par.m_N_sub_hx = N_sub_hx;  //[-]
@@ -3427,13 +3456,16 @@ void C_HX_co2_to_htf::initialize(int hot_fl, util::matrix_t<double> hot_fl_props
 	// Class is initialized
 	m_is_HX_initialized = true;
 
+    m_yr_inflation = yr_inflation;
+
 }
 
-void C_HX_co2_to_htf::initialize(int hot_fl, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type)
+void C_HX_co2_to_htf::initialize(int hot_fl, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type,
+                                 double yr_inflation)
 {
 	util::matrix_t<double> null_fluid_props;
 
-	initialize(hot_fl, null_fluid_props, N_sub_hx, od_UA_target_type);
+	initialize(hot_fl, null_fluid_props, N_sub_hx, od_UA_target_type, yr_inflation);
 }
 
 void C_HX_co2_to_htf::design_and_calc_m_dot_htf(C_HX_counterflow_CRM::S_des_calc_UA_par& des_par,
@@ -3691,7 +3723,8 @@ bool C_CO2_to_air_cooler::design_hx(S_des_par_ind des_par_ind, S_des_par_cycle_d
 	solver_code = -1;
 	i_W_par = -1;
 
-	while (solver_code != 0 || std::abs(T_hot_in_calc_2 - T_hot_in_calc) / T_hot_in_calc < 0.01)
+    while (solver_code != 0 || std::abs(T_hot_in_calc_2 - T_hot_in_calc) / T_hot_in_calc < 0.01
+        || ms_des_par_cycle_dep.m_T_hot_in_des - T_hot_in_calc_2 > 100.0)
 	{
 		i_W_par++;
 
@@ -3774,7 +3807,7 @@ bool C_CO2_to_air_cooler::design_hx(S_des_par_ind des_par_ind, S_des_par_cycle_d
 	ms_hx_des_sol.m_W_dot_fan = ms_des_par_cycle_dep.m_W_dot_fan_des;	//[MWe]
 
 	ms_hx_des_sol.m_cost_equipment = calculate_equipment_cost(ms_hx_des_sol.m_UA_total*1.E-3, ms_hx_des_sol.m_V_total,
-		ms_hx_des_sol.m_T_in_co2, ms_hx_des_sol.m_P_in_co2, ms_hx_des_sol.m_m_dot_co2);		//[M$]
+		ms_hx_des_sol.m_T_in_co2, ms_hx_des_sol.m_P_in_co2, ms_hx_des_sol.m_m_dot_co2, des_par_ind.m_yr_inflation);		//[M$]
 
     ms_hx_des_sol.m_cost_bare_erected = calculate_bare_erected_cost(ms_hx_des_sol.m_cost_equipment);
 
@@ -4361,14 +4394,24 @@ int C_CO2_to_air_cooler::C_MEQ_target_T_hot__width_parallel::operator()(double W
 }
 
 double /*M$*/ C_CO2_to_air_cooler::calculate_equipment_cost(double UA /*kWt/K*/, double V_material /*m^3*/,
-	double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/)
+	double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/, double yr_inflation/**/)
 {
 	switch (m_cost_model)
 	{
 	case C_CO2_to_air_cooler::E_CARLSON_17:
-		return 2.3*1.E-3*UA;		//[M$] needs UA in kWt/K
+    {
+        double yr_base_inflation = 2017;
+        double f_inflation = calculate_inflation_factor(yr_base_inflation, yr_inflation);
+        return 2.3 * 1.E-3 * UA * f_inflation;		//[M$] needs UA in kWt/K
+    }
+		
     case C_CO2_to_air_cooler::E_WEILAND_19:
-        return 32.88*std::pow(UA*1.E3, 0.75)*1.E-6; //[M$] needs UA in Wt/K
+    {
+        double yr_base_inflation = 2017;
+        double f_inflation = calculate_inflation_factor(yr_base_inflation, yr_inflation);
+        return 32.88 * std::pow(UA * 1.E3, 0.75) * 1.E-6 * f_inflation; //[M$] needs UA in Wt/K
+    }
+        
 	default:
 		return std::numeric_limits<double>::quiet_NaN();
 	}

tcs/heat_exchangers.cpp

@@ -478,6 +478,7 @@ class C_HX_counterflow_CRM
 
     int m_cost_model;		//[-]
     int m_od_solution_type; //[-]
+    double m_yr_inflation = 0; //[yr]
 
     bool m_is_single_node_des_set;
     NS_HX_counterflow_eqs::S_hx_node_info ms_node_info_des;
@@ -678,7 +679,8 @@ class C_HX_counterflow_CRM
 
     double calculate_equipment_cost(double UA /*kWt/K*/,
         double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/,
-        double T_cold_in /*K*/, double P_cold_in /*kPa*/, double m_dot_cold /*kg/s*/);
+        double T_cold_in /*K*/, double P_cold_in /*kPa*/, double m_dot_cold /*kg/s*/,
+        double yr_inflation /*yr*/);
 
     double calculate_bare_erected_cost(double cost_equipment /*M$*/);
 
@@ -838,9 +840,10 @@ class C_HX_co2_to_htf : public C_HX_counterflow_CRM
 	void design_and_calc_m_dot_htf(C_HX_counterflow_CRM::S_des_calc_UA_par &des_par, 
 		double q_dot_design /*kWt*/, double dt_cold_approach /*C/K*/, C_HX_counterflow_CRM::S_des_solved &des_solved);
 
-	virtual void initialize(int hot_fl, util::matrix_t<double> hot_fl_props, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type);
+	virtual void initialize(int hot_fl, util::matrix_t<double> hot_fl_props, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type,
+                            double yr_inflation);
 
-	virtual void initialize(int hot_fl, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type);
+	virtual void initialize(int hot_fl, int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type, double yr_inflation);
 
 };
 
@@ -861,7 +864,8 @@ class C_HX_co2_to_co2_CRM : public C_HX_counterflow_CRM
 
     }
 
-    virtual void initialize(int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type);
+    virtual void initialize(int N_sub_hx, NS_HX_counterflow_eqs::E_UA_target_type od_UA_target_type,
+                            double yr_inflation);
 };
 
 namespace N_compact_hx
@@ -907,10 +911,12 @@ class C_CO2_to_air_cooler
 
         double m_eta_fan;       //[-] Fan isentropic efficiency
         int m_N_nodes_pass;     //[-] Number of nodes per pass
+        double m_yr_inflation;  //[yr] Inflation year target
 
 		S_des_par_ind()
 		{
-			m_T_amb_des = m_elev = std::numeric_limits<double>::quiet_NaN();
+			m_T_amb_des = m_elev = m_yr_inflation =
+                std::numeric_limits<double>::quiet_NaN();
 
             // Set realistic default values so model can solve without inputs for these values
             m_eta_fan = 0.5;
@@ -1414,7 +1420,8 @@ class C_CO2_to_air_cooler
 		double & k_air /*W/m-K*/, double & Pr_air);
 
 	double /*M$*/ calculate_equipment_cost(double UA /*kWt/K*/, double V_material /*m^3*/,
-		double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/);
+		double T_hot_in /*K*/, double P_hot_in /*kPa*/, double m_dot_hot /*kg/s*/,
+        double yr_inflation /*yr*/);
 
     double /*M$*/ calculate_bare_erected_cost(double cost_equipment /*M$*/);
 };