>	restart:#"m08_p33"

Determinar las condiciones a la salida de una cámara adiabática donde se mezclan dos corrientes de aire húmedo, una de 0,01 kg/s que entra saturada a 5 ºC y la otra de 0,02 kg/s que entra a 40 ºC y 70% de humedad.
Datos:

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

>	su1:="Aire":su2:="H2O":dat:=[m1=0.01*kg_/s_,phi1=1,T1=(5+273)*K_,m2=0.02*kg_/s_,phi2=0.7,T2=(40+273)*K_];

Esquema:

>

>

Eqs. const.:

>	eqET:=eq1_12:Adat:=get_gas_data(su1):Adat:=subs(c[p]=c[pa],R=R[a],M=M[a],T[b]=nada,[Adat]):Wgdat:=get_gas_data(su2):Wgdat:=subs(c[p]=c[pv],R=R[v],M=M[v],[Wgdat]):Wldat:=get_liq_data(su2):Wdat:=op(Wgdat),Wldat:get_pv_data(su2):dat:=op(dat),Const,SI2,SI1:

Como la 1 está saturada, supongamos que la 3 sale sobresaturada (3=sat, 4= liq).

>

eqBMa:=ma1+ma2=ma3;eqBMw:=ma1*w1+ma2*w2=ma3*w3+m4;eqBE:=ma1*h1+ma2*h2=ma3*h3+m4*h4;w1_:=evalf(subs(dat,w(phi1,T1,p0))):'w1'=evalf(%,1);w2_:=evalf(subs(dat,w(phi2,T2,p0))):'w2'=evalf(%,2);eqm:=m=ma*(1+w);ma1_:=subs(dat,m1/(1+w1_)):'ma1'=evalf(%,2);ma2_:=subs(dat,m2/(1+w2_)):'ma2'=evalf(%,3);ma3_:=ma1_+ma2_:'ma3'=evalf(%,3);phi3_:=1;w3_:=subs(dat,w(phi3_,T3,p0));h1_:=subs(Adat,Wdat,dat,T=T1,dat,h(T,w1_)):'h1'=evalf(%/(1e3*J_/kJ_),3);h2_:=subs(Adat,Wdat,dat,T=T2,dat,h(T,w2_)):'h2'=evalf(%/(1e3*J_/kJ_),4);h3_:=subs(Adat,Wdat,dat,T=T3,dat,h(T,w3_));h4_:=subs(Wdat,T=T3,hl(T));

>	sol1_:=fsolve(subs(w1=w1_,w2=w2_,w3=w3_,ma1=ma1_,ma2=ma2_,ma3=ma3_,h1=h1_,h2=h2_,h3=h3_,h4=h4_,SI0,{eqBE,eqBMw}),{m4,T3},T3=200..400):'T3'=evalf(subs(sol1_,T3),3)*K_;'m4'=evalf(subs(sol1_,m4),2)*kg_;

luego no sale sobresaturada, y entonces:

>

m4:=0;w3_:=subs(w1=w1_,w2=w2_,w3=w3_,ma1=ma1_,ma2=ma2_,ma3=ma3_,solve(eqBMw,w3)):'w3'=evalf(%,2);h3_:=subs(Adat,Wdat,T=T3,h(T,w3_));T3_:=fsolve(subs(w1=w1_,w2=w2_,w3=w3_,ma1=ma1_,ma2=ma2_,ma3=ma3_,h1=h1_,h2=h2_,h3=h3_,SI0,eqBE),T3)*K_:'T3'=evalf(%,3);phi3_:=subs(dat,phi(w3_,T3_,p0)):'phi3'=evalf(%,2);

>