slutrapport - Kemiteknik

63Di=(4.9781*Tk^1.75)/(Ptot*10^5);Deff=Di*0.3/4;pch4=Fch4*Ptot/Ftot;pco=Fco*Ptot/Ftot;pco2=Fco2*Ptot/Ftot;ph2=Fh2*Ptot/Ftot;ph2o=Fh2o*Ptot/Ftot;z=1+kco*pco+kch4*pch4+kh2*ph2+kh2o*ph2o/ph2;Ch2o=Ptot*101.325/(R*Tk);r1=k1*(pch4*ph2o*(ph2^-2)-pco*ph2/kj1)/(z^2);r2=k2*(pco*ph2o*(ph2^-1)-pco2/kj2)/(z^2);r3=k3*(pch4*(ph2o^2)*(ph2^-3.5)-pco2*(ph2^0.5)/(kj1*kj2))/(z^2);th=Dp*sqrt((r1+r3)*rhokat/(Deff*Ch2o))/2;n(k,1)=(3/th)*(1/tanh(th)-1/th);endendSkriptfil: n_eff_sek.mfunction [n] = n_eff_sek(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);Fn2=F(k,6);T=F(k,7);Ptot=F(8);Tk=T+273.15;Ftot=Fch4+Fco+Fco2+Fh2+Fh2o+Fn2;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));Di=(4.9781*Tk^1.75)/(Ptot*10^5);Deff=Di*0.3/4;pch4=Fch4*Ptot/Ftot;pco=Fco*Ptot/Ftot;pco2=Fco2*Ptot/Ftot;ph2=Fh2*Ptot/Ftot;ph2o=Fh2o*Ptot/Ftot;z=1+kco*pco+kch4*pch4+kh2*ph2+kh2o*ph2o/ph2;Ch2o=Ptot*101.325/(R*Tk);r1=k1*(pch4*ph2o*(ph2^-2)-pco*ph2/kj1)/(z^2);r2=k2*(pco*ph2o*(ph2^-1)-pco2/kj2)/(z^2);r3=k3*(pch4*(ph2o^2)*(ph2^-3.5)-pco2*(ph2^0.5)/(kj1*kj2))/(z^2);th=Dp*sqrt((r1+r3)*rhokat/(Deff*Ch2o))/2;n(k,1)=(3/th)*(1/tanh(th)-1/th);endend