Untitled

mail@pastecode.io avatar
unknown
plain_text
2 years ago
3.1 kB
1
Indexable
{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}