/** Water heater synchronization determines the time to next synchronization state and the power drawn since last synch **/ TIMESTAMP waterheater::sync(TIMESTAMP t0, TIMESTAMP t1) { double internal_gain = 0.0; double nHours = (gl_tohours(t1) - gl_tohours(t0))/TS_SECOND; double Tamb = get_Tambient(location); // use re_override to control heat_needed state // runs after thermostat() but before "the usual" calculations if(re_override == OV_ON){ heat_needed = TRUE; } else if(re_override == OV_OFF){ heat_needed = FALSE; } if(Tw > 212.0 - thermostat_deadband){ // if it's trying boil, turn it off! heat_needed = FALSE; is_waterheater_on = 0; } TIMESTAMP t2 = residential_enduse::sync(t0,t1); // Now find our current temperatures and boundary height... // And compute the time to the next transition... //Adjusted because shapers go on sync, not presync set_time_to_transition(); // determine internal gains if (location == INSIDE){ if(this->current_model == ONENODE){ internal_gain = tank_UA * (Tw - get_Tambient(location)); } else if(this->current_model == TWONODE){ internal_gain = tank_UA * (Tw - Tamb) * h / height; internal_gain += tank_UA * (Tlower - Tamb) * (1 - h / height); } } else { internal_gain = 0; } // determine the power used if (heat_needed == TRUE){ /* power_kw */ load.total = (heat_mode == GASHEAT ? gas_fan_power : heating_element_capacity); is_waterheater_on = 1; } else { /* power_kw */ load.total = (heat_mode == GASHEAT ? gas_standby_power : 0.0); is_waterheater_on = 0; } //load.total = load.power = /* power_kw */ load.power; load.power = load.total * load.power_fraction; load.admittance = load.total * load.impedance_fraction; load.current = load.total * load.current_fraction; load.heatgain = internal_gain; waterheater_actual_power = load.power + (load.current + load.admittance * load.voltage_factor )* load.voltage_factor; actual_load = waterheater_actual_power.Re(); if (actual_load != 0.0) { prev_load = actual_load; power_state = PS_ON; } else power_state = PS_OFF; // gl_enduse_sync(&(residential_enduse::load),t1); if(re_override == OV_NORMAL){ if (time_to_transition >= (1.0/3600.0)) // 0.0167 represents one second { TIMESTAMP t_to_trans = (TIMESTAMP)(t1+time_to_transition*3600.0/TS_SECOND); return -(t_to_trans); // negative means soft transition } // less than one second means never else return TS_NEVER; } else { return TS_NEVER; // keep running until the forced state ends } }
/** oven synchronization determines the time to next synchronization state and the power drawn since last synch **/ TIMESTAMP range::sync(TIMESTAMP t0, TIMESTAMP t1) { double internal_gain = 0.0; double nHours = (gl_tohours(t1) - gl_tohours(t0))/TS_SECOND; double Tamb = get_Tambient(location); double dt = gl_toseconds(t0>0?t1-t0:0); if (oven_check == true || remainon == true) time_oven_operation +=dt; if (remainon == false) time_oven_operation=0; enduse_queue_oven += enduse_demand_oven * dt/3600/24; if (t0>TS_ZERO && t1>t0) { // compute the total energy usage in this interval load.energy += load.total * dt/3600.0; } if(re_override == OV_ON){ heat_needed = TRUE; } else if(re_override == OV_OFF){ heat_needed = FALSE; } if(Tw > 212.0 - thermostat_deadband){ // if it's trying boil, turn it off! heat_needed = FALSE; is_range_on = 0; } // determine the power used if (heat_needed == TRUE){ /* power_kw */ load.total = heating_element_capacity * (heat_mode == GASHEAT ? 0.01 : 1.0); is_range_on = 1; } else { /* power_kw */ load.total = 0.0; is_range_on = 0; } TIMESTAMP t2 = residential_enduse::sync(t0,t1); set_time_to_transition(); if (location == INSIDE){ if(this->current_model == ONENODE){ internal_gain = oven_UA * (Tw - get_Tambient(location)); } } else { internal_gain = 0; } dt = update_state(dt, t1); //load.total = load.power = /* power_kw */ load.power; load.power = load.total * load.power_fraction; load.admittance = load.total * load.impedance_fraction; load.current = load.total * load.current_fraction; load.heatgain = internal_gain; range_actual_power = load.power + (load.current + load.admittance * load.voltage_factor )* load.voltage_factor; actual_load = range_actual_power.Re(); if (heat_needed == true) total_power_oven = actual_load; else total_power_oven =0; if (actual_load != 0.0) { prev_load = actual_load; power_state = PS_ON; } else power_state = PS_OFF; // gl_enduse_sync(&(residential_enduse::load),t1); if(re_override == OV_NORMAL){ if (time_to_transition < dt) { if (time_to_transition >= (1.0/3600.0)) // 0.0167 represents one second { TIMESTAMP t_to_trans = (t1+time_to_transition*3600.0/TS_SECOND); return -(t_to_trans); // negative means soft transition } // less than one second means never else return TS_NEVER; } else return (TIMESTAMP)(t1+dt); } else { return TS_NEVER; // keep running until the forced state ends } }