| 查看: 277 | 回复: 0 | |||
xb490284223金虫 (正式写手)
蘑菇头
|
[求助]
密封舱浓度足迹法 求解方程问题
|
|
各位大神，对于下图的足迹法，我编写了如下的程序，但计算结果不满意，求各位帮忙。程序是 二分法求根的。 Cequ=1.5; slope=-0.21/3600; Intercept=-1.75; hm=0.0025; L=0.0038; B=0.33333; a=pi/2+0.00000001; b=pi*3/2-0.000000001; m=1; while(abs(m)>10^-10) q1=(a+b)/2; Dm=-slope*L*L/(q1*q1); %%求Dm Bim=hm*L/Dm; %%求Bim K=q1/((q1*q1/Bim-B)*tan(q1)); %%求K A1=(K*B-q1*q1*K/Bim+1)*cos(q1)-(1+2*K/Bim)*q1*sin(q1); m=log(-2*(K*B+1)*sin(q1)/q1/A1)-Intercept; Cm0=Cequ*K+1/B; %%求Cm0 if(m>0) a=q1; else b=q1; end if(abs(m)<10^-10) q1end=q1; Kend=K; Dmend=Dm; Cm0end=Cm0; m=abs(m); break; end end disp('q1是') q1end disp('初始散发浓度C0是') Cm0end disp('分离系数K是') Kend disp('扩散系数是') Dmend disp('m是') m  密封舱浓度足迹法流程图.png |
回复此楼» 猜你喜欢
有没有学校材料专业收跨调(一志愿085410)
已经有7人回复
-
311求调剂
已经有12人回复
-
268分085602化学工程调剂
已经有9人回复
-
初试分332，一志愿报考西北工业大学，
已经有5人回复
-
085801 总分275 本科新能源 求调剂
已经有9人回复
-
一志愿矿大，材料工程专硕314分，0856可调都可以
已经有4人回复
-
288求调剂
已经有6人回复
-
材料专硕初试分332一志愿西北工业大学，
已经有4人回复
-
材料专硕(0856) 339分求调剂
已经有8人回复
-
化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文
已经有11人回复
