Tema: problema con scilab!!
Ver Mensaje Individual
  #1 (permalink)  
Antiguo 22/10/2011, 04:43
phytona90
 
Fecha de Ingreso: diciembre-2010
Mensajes: 6
Antigüedad: 13 años, 5 meses
Puntos: 0
problema con scilab!!
Tengo que hace 1 programa con SCILAB que calcule la variación de Tº de la llama maxima durante la combustión de gas naturral en un proceso adiabatico en funcion de la concentración de metano,sabiendo que la concentración de este varía entorno al 95-63% del gas y etano 3-33%.Y HACER UNA GRAFICA que los confronte.
He plantEado lo siguiente pero me da error en x=fsolve y no se resolverlo.espero q alguien me eche 1 cable y veis algun error de planteamiento y demas comentarmelo por favor.gracias
Código scilab: 
Ver original
  1. Qcomb=(-212.798);
  2. Qcomb2=(-372.820);
  3. CH4=95; a1=3.381; b1=18.044*10^(-3);c1=-4.3*10^(-6);// Cantidad de entrada y constantes de metano
  4. C2H6=3 ;a2=2.247; b2=38.201*10^(-3);c2=-11.049^10^(-6);// Cantidad de entrada y constantes de etano
  5. N2=160;a3=6.524;b3=10.396*10^(-3);c3=-0.001*10^(-6);//Cantidad de entrada/salida de N2
  6. O2e=42;a4=6.148;b4=3.102*10^(-3);c4=-0.923*10^(-6);// Cantidad de entrada y constantes de OXIGENO
  7. CO2=98;a5=3.214;b5=10.396*10^(-3);c5=-3.545*10^(-6);// Cantidad de salida y constantes de CO2
  8. H2O=199;a6=7.256;b6=2.298*10^(-3);c6=0.283*10^(-6);// Cantidad de salida y constantes de H2O
  9.  
  10.  
  11. t0=298;
  12. i=0
  13.  
  14. scilab.org/product/dic-mat-sci/M2SCI_doc.htm">while i<=32//Condicion impuesta porque CH4 VARÍA DEL 95% AL 63%
  15.     CH4=CH4-i;
  16.     C2H6=C2H6+i
  17.     CO2=(CH4)+(2*C2H6)
  18.     H20=(2*CH4)+(3*C2H6)
  19.     O2s=42-((CH4*2)+(C2H6*(7/2)))
  20.     CONC=[]//Inicializamos vectores de % de ch4 y tº vacios para recoger una lista de valores y luego ilustrarlos
  21.     T_LLAMA=[]
  22.     //Calculo de la entalpia de entrada 
  23.     Cp1=(a1+(b1*t0)+(c1*(t0^2)))*10^(-3); Q1=CH4*Cp1// CP del metano y calor 
  24.     Cp2=(a2+(b2*t0)+(c2*(t0^2)))*10^(-3);Q2=C2H6*Cp2// CP del etano y calor
  25.     Cp3=(a3+(b3*t0)+(c3*(t0^2)))*10^(-3);Q3=Cp3*O2e// CP del O2 y calor
  26.     Cp4=(a4+(b4*t0)+(c4*(t0^2)))*10^(-3);Q4=Cp4*N2// CP del N2 y calor
  27.     He=Cp1+Cp2+Cp3+Cp4
  28.     //Calculo de la entalpia de salida
  29.     Cp5=(a3+(b3*t)+(c3*(t^2)))*10^(-3); Q5=N2*Cp5// CP del N2 y calor
  30.     Cp6=(a4+(b4*t)+(c4*(t^2)))*10^(-3); Q6=O2s*Cp6// CP del O2 y calor
  31.     Cp7=(a5+(b5*t)+(c5*(t^2)))*10^(-3); Q1=CO2*Cp7// CP del CO2 y calor
  32.     Cp8=(a6+(b6*t0)+(c6*(t^2)))*10^(-3); Q1=H2O*Cp8// CP del H2O y calor
  33.     //calculo de la entalpia de reacción
  34.     AHr=((CH4*Qcomb)+(C2H6*Qcomb2))
  35.     deff("fc=llama(t)",'fc=(Hs-He+AHr)'))
  36.     x=fsolve(298,llama)
  37.     printf(("La respuesta es :%2.3f"),x)
  38.     disp("Para este % de CH4:",CH4)
  39.     T_LLAMA=T_LLAMA+x
  40.     CONC=CONC+CH4
  41.     i=i+1
  42. scilab.org/product/dic-mat-sci/M2SCI_doc.htm">end
  43. plot(T_LLAMA,CONC)
  44. xtitle("Variación De La Concentración Frente a la Tº","Tº de Llama","% De CH4")