Untitled
unknown
plain_text
6 months ago
3.1 kB
0
Indexable
Never
{Define your known properties} P_high_turbine = 15 {MPa} T_high_turbine = 600 {C} P_low_turbine = 1 {MPa} T_low_turbine = 500 {C} P_condenser = 5 {kPa} P_extractor1 = 0,6 {MPa} P_extractor2 = 0,2 {MPa} Eta_turbine = 1 {Calculate enthalpy and entropy values for key points in the cycle} h[1] = Enthalpy(Steam;T=T_high_turbine;P=P_high_turbine) s[1] = Entropy(Steam;T=T_high_turbine;P=P_high_turbine) T[1] = T_high_turbine h_s[2] = Enthalpy(Steam;S=s[1];P=P_extractor1) s_s[2] = s[1] T_s[2] = Temperature(Steam;S=s[1];P=P_extractor1) h[2] = h[1] - (h[1]-h_s[2]) * Eta_turbine s[2] = Entropy(Steam;P=P_extractor1;h=h[2]) T[2] = Temperature(Steam;P=P_extractor1;h=h[2]) h_s[3] = Enthalpy(Steam;S=s[2];P=P_low_turbine) s_s[3] = s[2] T_s[3] = Temperature(Steam;S=s[2];P=P_low_turbine) h[3] = h[2] - (h[2]-h_s[3]) * Eta_turbine s[3] = Entropy(Steam;P=P_low_turbine;h=h[3]) T[3] = Temperature(Steam;P=P_low_turbine;h=h[3]) h_s[4] = Enthalpy(Steam;S=s[3];P=P_condenser) s_s[4] = s[3] T_s[4] = Temperature(Steam;S=s[3];P=P_condenser) h[4] = h[3] - (h[3]-h_s[4] ) * Eta_turbine s[4] = Entropy(Steam;P=P_condenser;h=h[4]) T[4] = Temperature(Steam;P=P_condenser;h=h[4]) h[5] = Enthalpy(Steam;P=P_extractor1;X=1) {saturated liquid from the closed feedwater heater} s[5] = Entropy(Steam;P=P_extractor1;X=1) T[5] = Temperature(Steam;P=P_extractor1;X=1) h[6] = h[5] {Assuming adiabatic throttling to the open feedwater heater} s[6] = Entropy(Steam;P=P_extractor2;h=h[6]) T[6] = Temperature(Steam;P=P_extractor2;h=h[6]) {Define the states at the exit of the feedwater heaters and the boiler} h[7] = Enthalpy(Steam;P=P_extractor2;X=0) {saturated liquid at the exit of the open feedwater heater} s[7] = Entropy(Steam;P=P_extractor2;X=0) T[7] = Temperature(Steam;P=P_extractor2;X=0) h[8] = h[1] + y*(h[6]-h[7]) {Energy balance equation for the open feedwater heater} s[8] = Entropy(Steam;P=P_pump2;h=h[8]) T[8] = Temperature(Steam;P=P_pump2;h=h[8]) h[9] = Enthalpy(Steam;P=P_pump1;X=0) {Saturated liquid at the exit of the boiler} s[9] = Entropy(Steam;P=P_pump1;X=0) T[9] = Temperature(Steam;P=P_pump1;X=0) h[10] = h[9] + (h[8]-h[9])/(1-y) {Energy balance equation for the boiler} s[10] = Entropy(Steam;P=P_high_turbine;h=h[10]) T[10] = Temperature(Steam;P=P_high_turbine;h=h[10]) {Define the mass and energy balance equations} "Mass balance equations for the open feedwater heater and the closed feedwater heater" 1 = (1-y) + x + y 1 = (1-y-x) + x "Energy balance equations for the open feedwater heater and the closed feedwater heater" h[1] = h[7] + y*(h[6]-h[7]) h[8] = h[9] + (1-y-x)*(h[10]-h[9]) {Define the efficiency of the cycle and the net power output} "W_net = m_dot * ((1-y) * (h[2]-h[1] + h[3]-h[4]) + y * (h[2]-h[5]))" Eta_cycle = W_net / (m_dot * (1-y-x) * (h[10]-h[9])) {Continue with mass and energy balance equations and then solve}{Note: You will need to define your turbine and pump efficiencies, mass flow rate, and set up energy balance equations for the heaters and turbine}