|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
程序如下：3 r. T% R% E7 p. E' H
function out=mult_syg7(z,x)$ _- d/ Q4 [! X$ ]8 i
Pi=3.14159;) `( m K1 A3 B" |
r0= 1.24*10^-3; % Outer diameter of the membrane, [m]
$ S3 H6 q- g, }3 j% _! iri= 9.4*10^-4;% Inner diameter of the membrane, [m]3 |8 X9 ~5 O$ j& `5 K
r=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]7 s3 L( E. e7 }8 j
L=0.05;% membrane length,[m]- p/ }8 G( ^* L9 b q7 D' A
Ac=Pi*r*L;% membrane area [m^2]7 ], |4 e7 p w; ^ S5 I; ]5 e
R=8.314;
3 v- s! F9 U) F& B" z0 D/ F; x1 zT0=1000;
V. G ~- C5 `, q. uT=T0+273;
8 G/ m# n; a; J$ h/ CP=101.325;%kpa
. e+ o/ L. S0 F7 ]: \F=2.678*10^-2*2;0 s s3 P: Y) N* A
m=3;
* W+ x$ F* C/ C2 |8 k% On=1;
/ `/ T% A9 S- ]& C8 ]8 Y5 |2 iS=(1+m+n+2*x(1))/P;
7 D; ^( Y/ U8 {2 yPCH4=(1-x(1))/S;
, w. a8 Z/ ~, y" ~7 KPH2O=(m-x(1)-x(2))/S;
# D$ ^7 V9 }0 ]9 h4 Yout= ...* O4 B& N+ y h: ~, F7 L; X, a H
[Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08# Y# u! I! T/ v$ X
(3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];) [: p& {9 x8 x/ P- b. S+ o
" k" n9 L) H5 P
调用如下:' w: y& e5 @' x# I7 m3 U, f" ]
>> x0 = [0; 0];
; }2 v) M9 i6 k0 i0 G Rtspan = [0:0.1:1];
. e" K& G$ ^$ ^: s. yM = [1 0 ; 0 0];( E* E- o" R: c# Z4 D
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
- y' d( N' {! u/ e3 o3 D: F7 {[t,x] = ode15s(@mult_syg7,tspan,x0,options)3 Y X" W% H* u! t! Z% w1 a
; X5 |( {6 j8 t报错如下：
3 {/ g. y6 K& k/ p+ }7 g错误使用 daeic12 (line 76)6 k2 M0 o- _& G$ V) a
此 DAE 的索引大于 1。
# V7 c- @/ _, C0 \- s1 q) P- z2 o9 p
出错 ode15s (line 310)0 J) k8 E, ^( b. a" ~2 ?# ?
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...& e5 b& H. `: E) y: P
|
|
/1 
发表于 2021-3-19 14:48
收藏
淘帖
支持！
反对！