Reconsidere o Prob.10-88. Usando o EES (ou outro programa), estude os efeitos da razão de pressão do ciclo a gás quando ela varia de 10 a 20 sobre a razão entre os fluxos de massa de gás e de vapor sobre a eficiência térmica do ciclo. Trace curvas para os resultados como funções da razão de pressões do ciclo a gás e discuta os resultados.

Sendo assim, o código comentado em EES é dado a seguir:

"Dados fornecidos no enunciado"

T[8] = 300 [K]

P[8] = 100 [kPa]

"Pr = 11"

T[10] = 1100 [K]

T[12] = 420 [K]

P[12] = P[8]

W_dot_liq=280000 [kW]

Eta_comp = 0,82

Eta_gas_turb = 0,86

Eta_bomba = 1

Eta_vapor_turb = 0,86

P[5] = 5000 [kPa]

T[5] =(350+273,15)

P[6] = 800 [kPa]

P[7] = 10 [kPa]

"Análise do ciclo a gás"

"Análise do compressor"

s[8]=ENTROPY(Air;T=T[8];P=P[8])

ss9=s[8]

P[9] = Pr*P[8]

Ts9=temperature(Air;s=ss9;P=P[9])

Eta_comp = w_gas_comp_isen/w_gas_comp

h[8] + w_gas_comp_isen =hs9

h[8]=ENTHALPY(Air;T=T[8])

hs9=ENTHALPY(Air;T=Ts9)

h[8] + w_gas_comp = h[9]

T[9]=temperature(Air;h=h[9])

s[9]=ENTROPY(Air;T=T[9];P=P[9])

"Análise da câmara de combustão"

h[9] + q_e = h[10]

h[10]=ENTHALPY(Air;T=T[10])

P[10]=P[9]

Q_dot_e=m_dot_gas*q_e

"Análise da turbina"

s[10]=ENTROPY(Air;T=T[10];P=P[10])

ss11=s[10]

P[11] = P[10] /Pr

Ts11=temperature(Air;s=ss11;P=P[11])

Eta_gas_turb = w_gas_turb /w_gas_turb_isen

h[10] = w_gas_turb_isen + hs11

hs11=ENTHALPY(Air;T=Ts11)

h[10] = w_gas_turb + h[11]

T[11]=temperature(Air;h=h[11])

s[11]=ENTROPY(Air;T=T[11];P=P[11])

"Análise do trocador de calor entre gás e vapor"

m_dot_gas*h[11] + m_dot_vapor*h[4] = m_dot_gas*h[12] + m_dot_vapor*h[5]

h[12]=ENTHALPY(Air; T=T[12])

s[12]=ENTROPY(Air;T=T[12];P=P[12])

"Análise do ciclo a vapor"

"Análise da bomba 1"

Fluid$='Steam'

P[1] = P[7]

P[2]=P[6]

h[1]=enthalpy(Fluid$;P=P[1];x=0)

v1=volume(Fluid$;P=P[1];x=0)

s[1]=entropy(Fluid$;P=P[1];x=0)

T[1]=temperature(Fluid$;P=P[1];x=0)

w_bomba1_s=v1*(P[2]-P[1])

w_bomba1=w_bomba1_s/Eta_bomba

h[1]+w_bomba1= h[2]

s[2]=entropy(Fluid$;P=P[2];h=h[2])

T[2]=temperature(Fluid$;P=P[2];h=h[2])

"Análise do aquecedor de água de alimentação"

y*h[6] + (1-y)*h[2] = 1*h[3]

P[3]=P[6]

h[3]=enthalpy(Fluid$;P=P[3];x=0)

T[3]=temperature(Fluid$;P=P[3];x=0)

s[3]=entropy(Fluid$;P=P[3];x=0)

"Análise da bomba 2"

P[4] = P[5]

v3=volume(Fluid$;P=P[3];x=0)

w_bomba2_s=v3*(P[4]-P[3])

w_bomba2=w_bomba2_s/Eta_bomba

h[3]+w_bomba2= h[4]

s[4]=entropy(Fluid$;P=P[4];h=h[4])

T[4]=temperature(Fluid$;P=P[4];h=h[4])

w_vapor_bombas = (1-y)*w_bomba1+ w_bomba2

"Análise da turbina a vpor"

h[5]=enthalpy(Fluid$;T=T[5];P=P[5])

s[5]=entropy(Fluid$;P=P[5];T=T[5])

ss6=s[5]

hs6=enthalpy(Fluid$;s=ss6;P=P[6])

Ts6=temperature(Fluid$;s=ss6;P=P[6])

h[6]=h[5]-Eta_vapor_turb*(h[5]-hs6)

T[6]=temperature(Fluid$;P=P[6];h=h[6])

s[6]=entropy(Fluid$;P=P[6];h=h[6])

ss7=s[5]

hs7=enthalpy(Fluid$;s=ss7;P=P[7])

Ts7=temperature(Fluid$;s=ss7;P=P[7])

h[7]=h[5]-Eta_vapor_turb*(h[5]-hs7)

T[7]=temperature(Fluid$;P=P[7];h=h[7])

s[7]=entropy(Fluid$;P=P[7];h=h[7])

h[5] = w_vapor_turb + y*h[6] +(1-y)*h[7]

"Análise do condensador"

(1-y)*h[7]=q_s+(1-y)*h[1]

Q_dot_s=m_dot_vapor*q_s

"Análise do ciclo"

RazaoMassica_GasParaVapor =m_dot_gas/m_dot_vapor

W_dot_liq=m_dot_gas*(w_gas_turb-w_gas_comp)+m_dot_vapor*(w_vapor_turb - w_vapor_bombas)

Eta=(W_dot_liq/Q_dot_e)

RazaoTrabalho=(m_dot_gas*w_gas_comp + m_dot_vapor*w_vapor_bombas)/(m_dot_gas*w_gas_turb + m_dot_vapor*w_vapor_turb)

W_dot_net_vapor = m_dot_vapor*(w_vapor_turb - w_vapor_bombas)

W_dot_net_gas = m_dot_gas*(w_gas_turb - w_gas_comp)

RazaoTrabLiquido_GasParaVapor = W_dot_net_gas/W_dot_net_vapor

Então, com o código acima, é possível estudar a influência da razão de pressão do ciclo a gás sobre a razão mássica entre gás e vapor. Essa influência é mostrada no gráfico a seguir:

E, a influência da razão de pressão sobre a eficiência é mostrada a seguir: