>	restart:#"m05_p31"

Determinar la potencia máxima obtenible de una corriente de 2 kg/s a 0,6 MPa, según sea de agua o de aire y esté a a 150 ºC o a 200 ºC (cuatro casos).

Datos:

>	read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):

>	su1:="Aire":su2:="H2O":dat:=[m=2*kg_/s_,p1=0.6e6*Pa_,T1=(150+273)*K_,T2=(200+273)*K_];

Esquema:

>

>

Eqs. const.:

>	Adat:=get_gas_data(su1):dat1:=op(dat),Const,Adat,SI2,SI1:gdat:=get_gas_data(su2):ldat:=get_liq_data(su2):dat2:=op(dat),Const,gdat,ldat,SI2,SI1:get_pv_data(su2):

Casos de H2O a T1 y T2

>

Phi:=m*(h1-h0-T0*(s1-s0));'hl(T)'=hl(T);'hv(T)'=hv(T);'sl(T)'=sl(T);'sv(T,p)'=sv(T,p);h0_:=subs(dat2,T=T0,dat2,hl(T)):'h0'=evalf(%/(1000*J_/kJ_));h1_:=subs(dat2,T=T1,dat2,hl(T)):'h1'=evalf(%/(1000*J_/kJ_));h2_:=subs(dat2,T=T2,dat2,hv(T)):'h2'=evalf(%/(1000*J_/kJ_));s0_:=evalf(subs(dat2,T=T0,dat2,sl(T))):'s0'=evalf(subs(m_^2=J_*s_^2/kg_,%));s1_:=evalf(subs(dat2,T=T1,dat2,sl(T))):'s1'=evalf(subs(m_^2=J_*s_^2/kg_,%));s2_:=evalf(subs(dat2,T=T1,p=p1,dat2,sv(T,p))):'s0'=evalf(subs(m_^2=J_*s_^2/kg_,%));Phi1_:=subs(h1=h1_,h0=h0_,s1=s1_,s0=s0_,dat2,Phi):'Phi1'=evalf(%/(1e3*W_/kW_));Phi2_:=subs(h0=h0_,h1=h2_,s0=s0_,s1=s2_,dat2,Phi):'Phi2'=evalf(%/(1e3*W_/kW_));

Casos de Aire a T1 y T2

>

h:=c[p]*(T-T0);s:=c[p]*ln(T/T0)-R*ln(p/p0);h0_:=subs(dat1,T=T0,dat1,h);h1_:=subs(dat1,T=T1,dat1,h):'h1'=evalf(subs(m_^2=J_*s_^2/kg_,%)/(1e3*J_/kJ_));h2_:=subs(dat1,T=T2,dat1,h):'h2'=evalf(subs(m_^2=J_*s_^2/kg_,%)/(1e3*J_/kJ_));s0_:=evalf(subs(dat1,T=T0,p=p0,dat1,s));s1_:=evalf(subs(dat1,T=T1,p=p1,dat1,s)):'s1'=evalf(subs(m_^2=J_*s_^2/kg_,%));s2_:=evalf(subs(dat1,T=T2,p=p1,dat1,s)):'s2'=evalf(subs(m_^2=J_*s_^2/kg_,%));Phi1_:=subs(h1=h1_,h0=h0_,s1=s1_,s0=s0_,dat2,Phi):'Phi1'=evalf(subs(m_^2=J_*s_^2/kg_,%)/(1e3*W_/kW_));Phi2_:=subs(h0=h0_,h1=h2_,s0=s0_,s1=s2_,dat2,Phi):'Phi2'=evalf(subs(m_^2=J_*s_^2/kg_,%)/(1e3*W_/kW_));

>