slutrapport - Kemiteknik

Skriptfil: comb_prim.mfunction [ut] = comb_prim()Fch4=4.7356*3.6;%kmol/hFco=0.1275*3.6;%kmol/hFh2=53.2822*0.2*3.6;%kmol/hbm=Fch4*22.414; %m3/hbc=Fco*22.414; %m3/hbv=Fh2*22.414; %m3/hQm=bm*35.33; %MJ/hQc=bc*12.48;%MJ/hQv=bv*10.62;%MJ/hQ=1.582502*3600; %MJ/hCpmc=1686.74e-6; %MJ/m3ºCCpmv=1480.91e-6; %MJ/m3ºCCpmm=1517.98e-6; %MJ/m3ºCCpml=1396e-6;%MJ/m3ºCg0c=2.885; %m3/m3g0v=2.885; %m3/m3g0m=10.54; %m3/m3ex=0.8;l0c=2.385; %m3/m3l0v=2.385; %m3/m3l0m=9.54; %m3/m3res=@(T)(bc*Cpmc*g0c*T+bv*Cpmv*g0v*T+bm*Cpmm*g0m*T+Q-Qm-Qc-Qv+bc*ex*Cpml*l0c*T+bv*ex*Cpml*l0v*T+bm*ex*Cpml*l0m*T);Tut=fsolve(res,800);Vn2=(g0c*bc*0.6529+g0v*bv*0.6529+g0m*bm*0.7149+bc*l0c*ex*0.79+bv*l0v*ex*0.79+bm*l0m*ex*0.79)/22.414; %kmol/hVco2=(g0c*bc*0.3471+g0m*bm*0.0950)/22.414; %kmol/hVh2o=(g0v*bv*0.3471+g0m*bm*0.1900)/22.414; %kmol/hVo2=(bc*l0c*ex*0.21+bv*l0v*ex*0.21+bm*l0m*ex*0.21)/22.414;%kmol/hut=[Tut Vn2 Vco2 Vh2o Vo2];endSkriptfil: n_eff_prim.mfunction [n] = n_eff_prim(W,F)R=8.3145;%J/molKrhokat=2835; %kg/m3Dp=0.0001; %mn=zeros(size(W,1),1);for k=1:size(W,1)Fch4=F(k,1);Fco=F(k,2);Fco2=F(k,3);Fh2=F(k,4);Fh2o=F(k,5);T=F(k,6);Ptot=F(k,7);Tk=T+273.15;Ftot=Fch4+Fco+Fco2+Fh2+Fh2o;k1=(4.23e15)*exp(-240.1e3/(R*Tk));k2=(2e6)*exp(-67.13e3/(R*Tk));k3=(1.02e15)*exp(-243.9e3/(R*Tk));kj1=exp(30.114-26830/Tk);kj2=exp(-4.036+4400/Tk);kco=(8.23e-5)*exp(70.65e3/(R*Tk));kch4=(6.65e-4)*exp(38.28e3/(R*Tk));kh2o=(1.77e5)*exp(-88.68e3/(R*Tk));kh2=(6.12e-9)*exp(82.90e3/(R*Tk));62