| ²é¿´: 478 | »Ø¸´: 0 | ||
happytanгæ (СÓÐÃûÆø)
|
[ÇóÖú]
ÇëÄÇλ°ïÎÒ¿´¿´Õâ¸öMATLAB³ÌÐòÄÇÀïÓÐÎÊÌ⣡ÏÈлл
|
|
function RH_Cat_Kinetics %ÓÃÊýÖµ»ý·Ö·¨½øÐз´Ó¦ËÙÂÊ·ÖÎöµÃµ½ËÙÂʳ£Êý£¬²¢½«½á¹ûÓëʵÑéÖµ¶Ô±È·ÖÎö clear all;clc global C kk0 t=[0:300:2400]; % t:ʱ¼ä£¬s C=[0.00 0.11 0.19 0.25 0.31 0.44 0.54 0.62 0.65; 0.00 0.10 0.18 0.24 0.28 0.39 0.48 0.52 0.53; 0.10 3.20 1.70 1.30 2.50 3.90 5.20 8.20 9.7]; C=C';%ζÈ150¡æµÄ¶¯Á¦Ñ§ÊµÑéÊý¾ÝC(1) &C(2): kmol/m^3, C(3):% %·ÇÏßÐÔÄâºÏ kk0=[5.1095e-004; 2.8939;9.5453e+003]; %ÀûÓÃsureshÄ£ÐÍËã³öͬÑùζÈÌõ¼þϵÄËÙÂʳ£Êý×÷ΪÆä³õʼֵ tspan=[0:300:2400]; %ʱ¼ä½×ÌÝ C0=[0;0;0.1]; % 0.1:ÑõÆøÔÚ·´Ó¦ÌåÏµÆøÏàÖеijõʼ·ÖÂÊ£¬¿¼ÂÇÁË»·¼ºÍéµÄ±¥ºÍÕô %Æûѹ£¬ÔÚʵÑéÌõ¼þÏÂΪ50%×óÓÒ£¬³ýÈ¥¸ÃÕôÆûѹÒÔ¼°¶èÐÔÆøÌåËù %Õ¼µÄ±ÈÀýºóµÃ³öµÄÉÏÊöÊýÖµ kk=lsqnonlin(@myfunc,kk0,0,inf,[ ],tspan,C0,C); %·ÇÏßÐÔ×îС¶þ³Ë·¨ÄâºÏ¡£µ÷Óú¯Êýmyfunc ci=nlparci(beta,resid,jacobian); %»Ø¹éϵÊý£¬²Ð²î£¬ÑŸó±È¾ØÕó %ÄâºÏЧ¹ûͼ(ʵÑéÖµÓëÄâºÏÖµµÄ±È½Ï) M=[ 1 0 0 0; %½â΢·Ö´úÊý·½³Ì×éʱµÄϵÊý¾ØÕó 0 1 0 0; 0 0 0 0; 0 0 0 0]; options=odeset; %²úÉú/¸Ä±ä²ÎÊý½á¹¹ options.Mass=M; options.TolRol=le-7;%¶¨Ò徫¶È tspan=[5:40]; C0=[0 0 0.1 ];%³õÖµC P=1.2;%MPa G=1.5;%L/min kla=0.14;%1/s; apxl=input('Pls.input gas holdup you got according reaction conditions using hydro_mass_ c:','s'); [t_sim,C_sim]=ode45(@lj_ssm,tspan,CO,options,kk,P,G,kla,apxl); plot(t,C_exp(:,1),'-ks',t_sim,C_sim(:,1),':ko',t,C_exp(:,2),'-kd',t_sim,C_sim(:,2),':k+'); legend('Exp.CP','Model.CP','Exp.CI','Mode1.CI'),xlabel('ʱ¼ät/s'),ylabel('Ũ¶ÈC/(kmol/m^3)'); %¶¨Òå·ÇÏßÐÔÓÅ»¯µÄÄ¿±êº¯Êý function f=myfunc(kk,tspan,C0,C) P=12; G=1.5; kla=0.14; apxl=0.2; [t,CC]=ode45(@lj_ssm,tspan,C0,[],kk,P,G,kla,apxl); f=CC-C; %¶¨Òå´ýÇóµÄ¶¯Á¦Ñ§·½³Ì function dCdt=lj_ssm(t,C,kk,P,G,kla,apxl) G=G/(60*22.4);%¸Ä±äÆøÌåÁ÷Á¿µÄµ¥Î»£¬±ãÓÚºóÃæµÄµÈʽ³ÉÁ¢£¬,mol/s yin=0.21; GI=G*(1-yin);%°´ÕÕ½øÆøÖÐÑõÆøº¬Á¿21%¼ÆËã¶èÐÔÆøÌåµÄÌå»ýÁ÷Á¿ Vm=22.4;%ÆøÌåµÄĦ¶ûÌå»ý£¬,L Vl=0.3; % VL:ÒºÌåÌå»ý£¬°´300mL¼ÆËã VG1=Vl/(1-apxl); %·´Ó¦ÌåϵµÄ×ÜÌå»ý£¬L P=P*10; % PÔÀ´µÄµ¥Î»ÊÇMPa£¬ÎªÁËÄܹ»ÀûÓúóÃæµÄ¹«Ê½£¬½«µ¥Î»×ª±äΪbar H=1.09e-2; % H:ÑõÆøµÄHenryϵÊý kk=kk0.*exp(-Ea/(R*T)); % T:K k0:ָǰÒò×ÓEa:»î»¯ÄÜ dCdt=zeros(1,4);%ÏÈÔ¤·ÖÅä¿Õ¼ä¸øÅ¨¶ÈC£¬ÕâÑù¿ÉÒÔ¼Ó¿ìµ÷ÓÃËÙ¶È dCPdt=kk(1)*kk(3)*C(1)*kk(4)/(kk(1)+kk(2)*C(1)+kk(3)*kk(4)); dC1dt=kk(3)*C(1)*kk(4)*(kk(1)-kk(2)*C(3))/(kk(1)+kk(2)*C(1)+kk(3)*kk(4)); gas_mass=Vm*(1-apxl)/(Vl*apxl)*(GI*(yin/(1-yin)-C(4)/(1-C(4)))-kla*(H*P*C(4)-kk(4))*VG1); liquid_mass=kla*(H*P*C(4)-C(2))-(kk(3)*C(1)*C(2)*(kk(1)+kk(2)*C(1))/(kk(1)+kk(2)*C(1)+kk(3)*C(2))); % gas mass & liquidse mass·Ö±ð´ïµ½ÎÈ̬ʱÑõÆøÔÚÆøÒºÏàÖеÄÎïÁÏÆ½ºâ·½³Ì×ó¶Ë dCdt=[dCPdt;dC1dt;gas_mass;liquid_mass]; |
» ²ÂÄãϲ»¶
Ò»Ö¾Ô¸°²´óÉúÎïѧ07³õÊÔ322¡¢±¾¿Æ¶þ±¾¡¢µ÷¼ÁÇóÖú
ÒѾÓÐ7È˻ظ´
0856ר˶Çóµ÷¼Á Ï£ÍûÊÇaÇøÔºÐ£
ÒѾÓÐ25È˻ظ´
085600²ÄÁÏÓ뻯¹¤£¬Çóµ÷¼Á
ÒѾÓÐ12È˻ظ´
Ò»Ö¾Ô¸2110£¬»¯Ñ§Ñ§Ë¶310·Ö£¬±¾¿ÆÖصãË«·ÇÇóµ÷¼Á
ÒѾÓÐ19È˻ظ´
²ÄÁϹ¤³Ì085601£¬270Çóµ÷¼Á
ÒѾÓÐ33È˻ظ´
¸´ÊÔµ÷¼Á
ÒѾÓÐ12È˻ظ´
ÐÂÒ»´úµç×ÓÐÅÏ¢294Çóµ÷¼Á ²»ÌôѧУ
ÒѾÓÐ8È˻ظ´
290µ÷¼ÁÉúÎï0860
ÒѾÓÐ29È˻ظ´
279ѧ˶ʳƷרҵÇóµ÷¼ÁԺУ
ÒѾÓÐ15È˻ظ´
310Çóµ÷¼Á
ÒѾÓÐ14È˻ظ´
ÕÒµ½Ò»Ð©Ïà¹ØµÄ¾«»ªÌû×Ó£¬Ï£ÍûÓÐÓÃŶ~
°ïæ¿´¿´Õâ¸öÆ½ÃæÍ¼ÓÃmatlabÔõô»
ÒѾÓÐ19È˻ظ´
ʦ¸çʦ½ã°ïÎÒ¿´¿´ÎÒÕâ¸öÑûÇ뺯·ûºÏÒªÇóÂð~~ÏÈллÁË
ÒѾÓÐ12È˻ظ´
Çë´ó¼Ò°ïÎÒ¿´¿´Õâ¸öºË´ÅÆ×ͼ°É£¬ÊµÔÚÊÇûÓа취ÁË
ÒѾÓÐ30È˻ظ´
¡¾ÇóÖú¡¿Çë½ÌÒ»¸ömatlab³ÌÐò£¬Ð»Ð»°ï棡£¡
ÒѾÓÐ15È˻ظ´
¡¾Çë½Ì¡¿´ó¼Ò°ïÎÒ¿´¿´ÕâÖÖ±¡Ä¤ÑùÆ·µÄ½ØÃæSEMÕÕÆ¬ÊÇÔõô»ñµÃµÄ£¬Ð»Ð»ÁË
ÒѾÓÐ18È˻ظ´
¡¾ÇóÖú¡¿¼±¼±£¡£¡£¡ÄÄλ¶®xps£¬°ïÎÒ¿´Ò»ÏÂÎÒµÄÊý¾Ý°É¡£Ð»Ð»
ÒѾÓÐ5È˻ظ´
¡¾ÇóÖú¡¿´ó¼Ò°ïÎÒ¿´¿´Õâ¸öMATLAB³ÌÐò£¿
ÒѾÓÐ6È˻ظ´
¿ÆÑдÓСľ³æ¿ªÊ¼£¬ÈËÈËΪÎÒ£¬ÎÒΪÈËÈË













»Ø¸´´ËÂ¥
µã»÷ÕâÀïËÑË÷¸ü¶àÏà¹Ø×ÊÔ´
5