slutrapport - Kemiteknik

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=(3/th)*(1/tanh(th)-1/th);r1=n*k1*(pch4*ph2o*(ph2^-2)-pco*ph2/kj1)/(z^2);r2=n*k2*(pco*ph2o*(ph2^-1)-pco2/kj2)/(z^2);r3=n*k3*(pch4*(ph2o^2)*(ph2^-3.5)-pco2*(ph2^0.5)/(kj1*kj2))/(z^2);dFdW(1)=-r1-r3;dFdW(2)=r1-r2;dFdW(3)=r2+r3;dFdW(4)=3*r1+r2+4*r3;dFdW(5)=-r1-r2-2*r3;dFdW(6)=((4*U*(Tvvx-T)/(dtub*rhokat))-(r1*dH1+r2*dH2+r3*(dH1+dH2)))/(Ftot*cpm);dFsum=dFdW(1)+dFdW(2)+dFdW(3)+dFdW(4)+dFdW(5);dFdW(7)=(Ptot0/(Ftot0*Tk0))*(Ftot*dFdW(6)+Tk*dFsum);dFdW=dFdW';endendSkriptfil: sek.mfunction [wut,Fut]=sek(Fin,Ptot0)R=8.3145;%J/molKrhokat=1701;cpm=32.4e3;%J/kmoldH1=205.9e6; %J/kmoldH2=-41e6;%J/kmolDp=0.0001; %mTk0=Fin(7);Ftot0=sum([Fin(1) Fin(2) Fin(3) Fin(4) Fin(5) Fin(6)]);init=[Fin(1) Fin(2) Fin(3) Fin(4) Fin(5) Fin(6) Fin(7) Ptot0];%init=[kmol/h kmol/h kmol/h kmol/h kmol/h Kelvin bar][wut,Fut]=ode15s(@diffen,[0 5],init);function [dFdW]=diffen(w,F)Fch4=F(1);Fco=F(2);Fco2=F(3);Fh2=F(4);Fh2o=F(5);Fn2=F(6);T=F(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);60