Skip to content

Commit

Permalink
Add physical heat sink variables to IPH trough cmod.
Browse files Browse the repository at this point in the history
  • Loading branch information
@taylorbrown75
taylorbrown75 committed Nov 22, 2024
1 parent 901ddf9 commit 8e67452
Showing 1 changed file with 22 additions and 2 deletions.
24 changes: 22 additions & 2 deletions ssc/cmod_trough_physical_iph.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -169,7 +169,16 @@ static var_info _cm_vtab_trough_physical_iph[] = {
{ SSC_INPUT, SSC_NUMBER, "p_start", "Collector startup energy, per SCA", "kWhe", "", "solar_field", "*", "", "" },

// Heat Sink
{ SSC_INPUT, SSC_NUMBER, "pb_pump_coef", "Pumping power to move 1kg of HTF through PB loop", "kW/kg", "", "Heat Sink", "*", "", "" },
{ SSC_INPUT, SSC_NUMBER, "pb_pump_coef", "Pumping power to move 1kg of HTF through PB loop", "kW/kg", "", "Heat Sink", "*", "", "" },

{ SSC_INPUT, SSC_NUMBER, "hs_type", "0: ideal model, 1: physical steam model", "", "", "Heat Sink", "?=0", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_N_sub", "Number physical heat sink HX nodes", "", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_tol", "Physical heat sink solve tolerance", "", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_f_mdot_steam_min", "Min steam mdot fraction for physical heat sink", "", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_f_mdot_steam_max", "Max steam mdot fraction for physical heat sink", "", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_T_steam_cold_des", "Steam inlet temperature for physical heat sink", "C", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_P_steam_hot_des", "Steam outlet (and inlet) pressure for physical heat sink", "bar", "", "Heat Sink", "hs_type=1", "", "" },
{ SSC_INPUT, SSC_NUMBER, "hs_phys_Q_steam_hot_des", "Steam outlet quality for physical heat sink", "", "", "Heat Sink", "hs_type=1", "", "" },

// Parallel heater parameters
{ SSC_INPUT, SSC_NUMBER, "heater_mult", "Heater multiple relative to design cycle thermal power", "-", "", "Parallel Heater", "is_parallel_htr=1", "", "" },
Expand Down Expand Up @@ -1492,7 +1501,12 @@ class cm_trough_physical_iph : public compute_module
}

// Heat Sink
int hs_type = as_integer("hs_type");
C_csp_power_cycle* c_heat_sink_pointer;
C_pc_heat_sink c_heat_sink;

// Ideal heat sink
if(hs_type == 0)
{
//size_t n_f_turbine1 = 0;
//ssc_number_t* p_f_turbine1 = as_array("f_turb_tou_periods", &n_f_turbine1); // heat sink, not turbine
Expand Down Expand Up @@ -1522,6 +1536,12 @@ class cm_trough_physical_iph : public compute_module
c_heat_sink.mc_reported_outputs.assign(C_pc_heat_sink::E_M_DOT_HTF, allocate("m_dot_htf_heat_sink", n_steps_fixed), n_steps_fixed);
c_heat_sink.mc_reported_outputs.assign(C_pc_heat_sink::E_T_HTF_IN, allocate("T_heat_sink_in", n_steps_fixed), n_steps_fixed);
c_heat_sink.mc_reported_outputs.assign(C_pc_heat_sink::E_T_HTF_OUT, allocate("T_heat_sink_out", n_steps_fixed), n_steps_fixed);

c_heat_sink_pointer = &c_heat_sink;
}
else
{
throw exec_error("trough_physical_iph", "hs_type != 0; other heat sink models are not currently supported");
}

// Electricity pricing schedule
Expand Down Expand Up @@ -1665,7 +1685,7 @@ class cm_trough_physical_iph : public compute_module
// Instantiate Solver
C_csp_solver csp_solver(weather_reader,
c_trough,
c_heat_sink,
*c_heat_sink_pointer,
*storage_pointer,
tou,
dispatch,
Expand Down

0 comments on commit 8e67452

Please sign in to comment.