Skip to content

Commit

Permalink
Add physical heat sink variables to IPH MSLF cmod.
Browse files Browse the repository at this point in the history
  • Loading branch information
@taylorbrown75
taylorbrown75 committed Nov 22, 2024
1 parent 366ae1a commit 71a0b4b
Showing 1 changed file with 23 additions and 2 deletions.
25 changes: 23 additions & 2 deletions ssc/cmod_fresnel_physical_iph.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -153,7 +153,17 @@ static var_info _cm_vtab_fresnel_physical_iph[] = {

// 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", "", "" },


// TES

Expand Down Expand Up @@ -911,7 +921,12 @@ class cm_fresnel_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 @@ -940,6 +955,12 @@ class cm_fresnel_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("fresnel_physical_iph", "hs_type != 0; other heat sink models are not currently supported");
}

// Electricity pricing schedule
Expand Down Expand Up @@ -1087,7 +1108,7 @@ class cm_fresnel_physical_iph : public compute_module
// Instantiate Solver
C_csp_solver csp_solver(weather_reader,
c_fresnel,
c_heat_sink,
*c_heat_sink_pointer,
storage,
tou,
dispatch,
Expand Down

0 comments on commit 71a0b4b

Please sign in to comment.