>	restart:#"m09_pp5

Un motor cohete de propulsante líquido usa keroseno (supóngase dodeceno) y N2O4. Suponiendo reacción estequiométrica completa, se pide:
a) Masa molar y relación de capacidades térmicas de los gases producto.
b) Energía liberada por unidad de masa de keroseno, y comparación con el PCI de la combustión de keroseno en aire.
c) Temperatura de combustión adiabática.
d) Velocidad de salida de los gases para una presión de cámara de 7 MPa y salida adaptada a 100 kPa.
Datos:

>	read`../therm_eq.m`:read`../therm_chem.m`:with(therm_chem);with(therm_proc):#Global Elem,Comp,C_,MfComp,cpComp,eqNX,eqBC,eqBH,eqBO,eqBN,eqBS,delta_

>	su1:="C12H24":su2:="N2O4":su3:="N2":su4:="CO2":su5:="O2":su6:="H2O":dat:=[pc=7e6*Pa_]:

Eqs. const.:

>	dat:=op(dat),op(subs(g=g0,[Const])),SI2,SI1:N2dat:=get_gas_data("N2"):CO2dat:=get_gas_data("CO2"):O2dat:=get_gas_data("O2"):H2Odat:=get_gas_data(su6):eqF:=Mf(su1);eqNTO:=Mf(su2);

a) Masa molar y relación de capacidades térmicas de los gases producto.

>

eq0:=eq_fit(C12H24+a*N2O4=b*CO2+d*H2O+e*N2);eq7_2;Mm:=x[CO2]*M[CO2]+x[H2O]*M[H2O]+x[N2]*M[N2];eq0_:=eqMIX(a*C12H24+b*N2O4=[3,4,6]);sol_:=evalf(solve({eqBC,eqBH,eqBO,eqBN,eqNX},{a,b,x[CO2],x[H2O],x[N2]}));Mm_:=subs(sol_,M[CO2]=M1,M[H2O]=M2,M[N2]=M3,M1=M,CO2dat,M2=M,H2Odat,M3=M,N2dat,Mm);eqcpmol:=lhs(sum(delta_[i]*x[Comp[i]]*cpComp[i],i=1..C_))=subs(sol_,rhs(sum(delta_[i]*x[Comp[i]]*cpComp[i],i=1..C_)));

	(1)

i.e. la masa molar de la mezcla es 0,030 kg/mol, y su capacidad térmica media a alta temperatura (i.e. con salida a unos 2000..2500 K) cp=46 J/(mol·K).

b) Energía liberada por unidad de masa de keroseno, y comparación con el PCI de la combustión de keroseno en aire.

>	eq15_5_;eqPCs;PCI_NTO:=PCI(eq0);PCIm_NTO:=%/rhs(eqF);eq1:=eq_fit(C12H24+a*O2=b*CO2+d*H2O+e*N2);PCI_O2:=PCI(eq1);PCIm_O2:=%/rhs(eqF);

	(2)

Se liberan 44,5 MJ/kg cuando se usa NTO, y 44,0 MJ/kg cuando se usa O2 (da igual si está disuelto con N2). La pequeña diferencia es debida a la entalpía de formación del N2O4. Pero la enorme ventaja del NTO es que es fácilmente licuable (a T ambiente, no a T criogénica), y además genera más gases.

c) Temperatura de combustión adiabática.

>	subs(sol_,eq0_);eq15_7_2;Ta_:=subs(PCI_=PCI_NTO,sol_,cpComp,dat,rhs(eqTa));'Ta_'=TKC(%);

	(3)