>	restart:#"m15_p65"

>	read`../therm_chem.m`:with(therm_chem);with(therm_proc):

Un motor diésel de 250 kW y rendimiento energético 0,4 quema un combustible asimilable al n-dodecano, con un 80% de exceso de aire. Se pide:
a) Calcular los poderes caloríficos inferior y superior del combustible.
b) Consumos de combustible y aire.
c) Composición en el escape.
d) Emisión de CO2 en kg/h y en kg/kWh.
e) Temperatura de rocío de los gases de escape.

Datos:

>	su1:="Aire":su2:="H2O":fuel:=C12H26:dat:=[Wdot=250e3*W_,eta=0.4,e=0.80];

Eqs. const.:

>	Adat:=get_gas_data(su1):Wdat:=get_gas_data(su2),get_liq_data(su2):get_pv_data(su2):dat:=op(dat),op(subs(g=g0,[Const])),Adat,SI2,SI1:

a) Calcular los poderes caloríficos inferior y superior del combustible.

>	eq15_5;eq1:=C12H26+(37/2)*O2=12*CO2+13*H2O;PCS_:=PCS(eq1);PCI_:=PCI(eq1);Mf_:=rhs(Mf(fuel));PCSm_:=subs(Adat,PCS_/Mf_);PCIm_:=subs(Adat,PCI_/Mf_);

i.e. en base molar PCS=8,1 MJ/mol, y en base másica PCS=48 MJ/kg de combustible.

b) Consumos de combustible y aire.

>	eq15_2;eqA0_:=Ateo(fuel);eqe:=eq15_4;A_:=solve(subs(A[0]=rhs(eqA0_),dat,eqe),A);A0m_:=subs(Adat,rhs(eqA0_)*M/Mf_);Am_:=subs(Adat,A_*M/Mf_);eqBE:=Wdot=eta*mf*PCIm;mf_:=subs(PCIm=PCIm_,dat,solve(%,mf));ma:=mf*Am;ma_:=mf_*Am_;

i.e. el motor ha de tomar 14 g/s de combustible y 0,38 kg/s de aire, para generar los 250 kW.

c) Composición en el escape.

>	eqMIX(a*fuel+a*A_*(c21*O2+c79*N2)=[2,3,4,6]);eqNX;eqBC;eqBH;eqBO;eqBN;sol1:=subs(dat,solve({eqNX,eqBC,eqBH,eqBO,eqBN},{a,x[Comp[2]],x[Comp[3]],x[Comp[4]],x[Comp[5]]}));

i.e. 76% de N2, 9,0% de O2, 7,9% de H2O, y 7,3% de CO2.

d) Emisión de CO2 en kg/h y en kg/kWh.

>	ndotCO2:=(nf/a)*x[CO2];mdotCO2_:=subs(nf=mf_/Mf_,sol1,dat,%*0.044*kg_/mol_);%*3600*s_/h_;eqE:=mdotCO2/Wdot=subs(dat,mdotCO2_/Wdot);eqE:=lhs(%)=subs(dat,rhs(%)*3.6e6*J_/kWh_);

i.e. se emiten 44 g/s de CO2 (0,63 kg/kWh).

e) Temperatura de rocío de los gases de escape.

>	xH2O_:=subs(sol1,x[H2O]);eqTR:=x[v]=p[v](T[R])/p;eqTR_:=xH2O_=subs(dat,pv(T)/p0);TR_:=solve(%,T);'TR'=TKC(%);

>