Untitled

"given"

p_b=150 [bar]

T_max= 500 [C]

p_reheat = 30 [Bar]

T_reheat= 500 [c]

P_1= 5 [bar]


p_c = 0.1 [bar]

eta_turbine=0.85

eta_pump=0.95



{point 1} 

p[1]=p_c 

x[1]=0

h[1]=enthalpy(water,p=p[1],x=x[1])

s[1]=entropy(water,p=p[1],x=x[1])


{point 2} 

p[2]=P_1

{isentropic}

s_s[2]=s[1]

h_s[2]=enthalpy(water,p=p[2],s=s_s[2])

{Actual} 

eta_pump*(h[1]-h[2])=(h[1]-h_s[2])

{point 3}

p[3]=P_1

x[3]=0

h[3]=enthalpy(water,p=p[3],x=x[3])

s[3]=entropy(water,p=p[3],x=x[3])


{point 4}

p[4]=p_b

{isentropic}

s_s[4]=s[3]

h_s[4]=enthalpy(water,p=p[4],s=s_s[4])

{Actual} 

eta_pump * ( h[3]-h[4] ) = ( h[3] - h_s[4] )


{point 5}

p[5]=p_b

T[5]=T_max

h[5]=enthalpy(water,p=p[5],T=T[5])

S[5]=entropy(water,p=p[5],T=T[5])


{------------------------------------------------------}
{point 6} 

p[6] =p_reheat

{isentropic} 

s_s[6]=s[8]

h_s[6]=enthalpy(water,p=p[6],s=s_s[6])

{Actual} 

eta_turbine * ( h_s[6]-h[8] )=( h[6] - h[8] )
{------------------------------------------------------}


{point 7}

p[7]=p_reheat

T[7]=T_reheat

h[7]=enthalpy(water,p=p[7],T=T[7])

s[7]=entropy(water,p=p[7],T=T[7])


{------------------------------------------------------}
{point 8}

p[8] =P_1

{isentropic} 

s_s[8]=s[5]

h_s[8]=enthalpy(water,p=P_1,s=s_s[8])

s[8]=entropy(water,p=p[8],h=h[8])

{Actual} 

"eta_turbine*(h_s[8]-h[7])=(h[8]-h[7])"
eta_turbine = (h[5] - h[8]) / (h[5] - h_s[8])
{------------------------------------------------------}



{------------------------------------------------------}
{point 9}

p[9] =p_c

{isentropic} 

s_s[9]= s[7]

h_s[9]=enthalpy(water,p=p[9],s=s_s[9])

{Actual} 

eta_turbine*(h_s[9]-h[7])=(h[9]-h[7])
{------------------------------------------------------}


{mass balance}

(1-m_1)*(h[2])+(m_1)*(h[8]) = h[3]

{Work net}


{------------------------------------------------------}
qadd=( h[5]-h[4] ) + (1-m_1)*( h[7]-h[6] ) 
{------------------------------------------------------}

qrej=( (1-m_1) * (h[9]-h[1]) )
w_net= qadd-qrej

{w_net=( h[5]-h[4] ) + ( h[7]-h[6] ) - ( (1-m_1) * (h[9]-h[1]) )}



{efficiency}

{eta_th= (( (h[5]-h[4]) + (h[7]-h[6])) - ( (1-m_1)*(h[9]-h[1]) )) / ((h[5]-h[4]) + (h[7]-h[6]))}

eta_th=w_net / qadd