>	restart:#"m16_p25"

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

Se ha llenado con metanol un pequeño matraz de cuello largo, de 1 cm de diámetro, quedando 5 cm libres a la atmósfera, que está a 25 ºC, 93 kPa y 50% de humedad relativa. Se pide:
a) Composición esperada de la fase gaseosa en las proximidades del líquido (molar y másica).
b) Velocidad media de escape en régimen casi-estacionario, suponiendo que la convección del aire exterior permite suponer que la concentración de alcohol en la boca es prácticamente nula, y aproximando la difusividad másica por el valor de la difusividad térmica.

c) Composición estequiométrica para la combustión y límites de ignición, en % de combustible en la mezcla..

d) Indicar qué ocurriría al acercar una llama a la boca del matraz. Y si se agitara antes (e.g. tapando la boca y dándole la vuelta un par de veces).

Datos:

>	fuel:=CH4O;fstr:=convert(fuel,string):su0:="Aire":dat:=[D=1e-2*m_,L=5e-2*m_,T0=(25+273)*K_,p0=93e3*Pa_];

Eqs. const.:

>	dat:=op(dat),op(subs(g=g0,[Const])),SI2,SI1:Adat:=get_gas_data(su0):Fdat:=get_gas_data(fstr),get_liq_data(fstr):get_pv_data(fstr):fuel_data:=subs(Fdat,dat,[rho_=rho,M_=M,a_=k*R*T0/(c[p]*p0),nu_=mu*R*T0/p0]);

a) Composición esperada de la fase gaseosa en las proximidades del líquido.

Aire saturado de alcohol, que se podrá aproximar como mezcla ideal de aire seco (de 79%N2 más 21%O2 en volumen) y alcohol a Tamb (despreciando la pequeña proporción de vapor de agua, del orden del 1%, y el pequeño enfriamiento evaporativo del alcohol, que se podrían calcular después).

>

x[F,0]:='pv(T0)/p0';xF0:=evalf(subs(dat,%)):'x[F,0]'=evalf(%,2);x[N,0]:='c79*(1-x[F,0])';xN0:=evalf(subs(dat,%)):'x[N,0]'=evalf(%,2);x[O,0]:='c21*(1-x[F,0])';xO0:=evalf(subs(dat,%)):'x[O,0]'=evalf(%,2);eqxy:=y[i]=x[i]*M[i]/Sum(x[i]*M[i],i=1..n);y[F,0]:='x[F,0]*M[F]/(x[F,0]*M[F]+(1-x[F,0])*M[A])';yF0:=evalf(subs(M[A]=Ma,M[F]=M,Fdat,Ma=M,Adat,dat,%)):'y[F,O]'=evalf(%,2);y[N,0]:=subs(Fdat,xN0*(0.028*kg_/mol_)/(xN0*(0.028*kg_/mol_)+xF0*M+xO0*(0.032*kg_/mol_))):'y[N,0]'=evalf(%,2);'y[O,0]'=evalf(1-y[N,0]-yF0,2);

i.e. 18% en volumen de alcohol, o 19% en masa, y el resto aire.

b) Velocidad media de escape en régimen casi-estacionario, suponiendo que la convección del aire exterior permite suponer que la concentración de alcohol en la boca es prácticamente nula, y aproximando la difusividad másica por el valor de la difusividad térmica.

Pese a llegar a un 19% en masa de alcohol, supondremos densidad constante e igual a la del aire.

Las fórmulas son:

>

eqvG:=v[gas]='Di*(y[F,0]-y[F,L])/L';eqTCG:=Di=k/(rho*c[p]);Di_:=subs(Fdat,dat,k*R*T0/(c[p]*p0)):'Di'=evalf(%,2);eqvG_:=v[gas]=subs(dat,Di_*(yF0-0)/L):evalf(%,2);eqvG_*60*s_/min_;eqvL:=v[liq]='v[gas]*rho[gas]/rho[liq]';eqvL:=subs(eqvG_,dat,Adat,Fdat,dat,v[liq]=v[gas]*(p0/(R*T0))/rho);v[liq]=evalf(rhs(%)*86400*s_/day_);

i.e. sale gas a una velocidad de 2,2 mm/min, retrocediendo el nivel de líquido a una velocidad 4,3 mm/día.

Nótese que todo lo que sale es alcohol y no hay movimiento de aire seco.

La paradoja de que sólo sale alcohol y se supone que fuera no hay, i.e. y[F,L]=0 se explica por el balance másico rho*v=rho[i]*v-Di*drho[i]/dz.

La deducción es:

>

x:='x':eqBM:=rho*V=constant;eqBMi:='diff(rhoi(z,t),t)=Di*diff(rhoi(z,t),z,z)-diff(rhoi(z,t)*v(z,t),x)';eqBMi_:=0=Di*diff(yi(z),z,z)-diff(yi(z)*V,z);eqBC1:=yi(0)=yi0;eqBC2:=yi(L)=yi1;dsolyi:=dsolve({eqBMi_,eqBC1,eqBC2},yi(z));dsolyilin:=yi(z)=collect(simplify(convert(series(rhs(dsolyi),V=0,3),polynom)),z);eqFick:=rho*V=-rho*Di*Diff(yi(z),z);eqFick_:=V=-Di*(yi0-yi1)/L;eqFick__:=subs(yi0=yF0,yi1=0,Di=Di_,dat,eqFick_);plot(subs(Di=Di_,eqFick__,yi0=yF0,yi1=0,dat,SI0,[rhs(dsolyi),rhs(dsolyilin)]),z=0..subs(dat,SI0,L),y=0..0.2);

Nótese que la aproximación lineal es muy buena, y que se ha tomado Di=a=k/(rho*cp)=9e-6 m2/s, cuando el valor experimental es Di=14e-6 m2/s, y la otra aproximación sería Di=mu/rho=12e-6 m2/s (la teoría cinética de gases predice valore iguales para los tres valores).

Aunque en la práctica hubiera un 2% o un 4% de alcohol a la salida, apenas influiría en el cálculo de la velocodad de difusión.

c) Composición estequiométrica para la combustión y límites de ignición, en % de combustible en la mezcla..

>	eq1:=CH4O+(3/2)*O2=CO2+2*H2O;evalf(Ateo(fuel));x[STQ]:=1/(1+A[0]);xSTQ_:=subs(A[0]=rhs(Ateo(fuel)),x[STQ]):'x[STQ]'=evalf(%);x[LFL]:='x[STQ]/2';'x[LFL]'=evalf(xSTQ_/2,2);

En efecto, de las tables se obtiene un límite de ignición pobre LFL=6,7% y un límite de ignición rico UFL=36% (hasta 50% según las fuentes).

d) Indicar qué ocurriría al acercar una llama a la boca del matraz. Y si se agitara antes (e.g. tapando la boca y dándole la vuelta un par de veces).

La respuesta es que depende de las condiciones de salida reales, ya que la suposición inicial de que la concentración de alcohol en la boca es nula por estar bien aireada es sólo una aproximación para estimar el gasto, pero no valdrá aquí, donde habría que estimar y[F,L] por acoplamiento entre la difusión en el tubo y la convección exterior. Lo que sí está claro es que la concentración nunca podrá superar el 18% en volumen, pero este valor está dentro del margen de flamabilidad [6,7%..36%].

En conclusión, es de esperar que en el problema difusivo estacionario, a la salida sea x[F,L]<6,7% y por tanto que no se encienda.

Sin embargo, si se agitara el contenido, la concentración sería uniforme en la fase gaseosa e igual a la de saturación (si hay suficiente líquido), luego x[F,L]=18% comprendido entre [6,7%..36%], luego se formaría una llama de premezcla que viajaría hacia el interior a unos 0,5 m/s (o más, por el efecto pared) hasta alcanzar el líquido y apagarse por falta de comburente, produciendo el típico sonido de retroceso de llama, pero sin peligro de explosión porque hay escape libre (sí podría ocurrir si el cuello fuese muy estrecho y largo, y hubiese sólo un poco de alcohol en el fondo del matraz, por imposibilidad de salida rápida de los gases de combustión).

>