شبیه سازی راکتور CSTR در MATLAB

پیرجو

مدیر ارشد
مدیر کل سایت
مدیر ارشد
CONTINUOUS STIRRED TANK REACTORS​

کد:
% rate_T
function fn=rate_T(c)
global beta gamma flowvol
T=1+beta*(1-c)
rate=c*exp(gamma*(1-1/T))
fn=flowvol*(1-c)-rate
>> rate_T(0.5)
T=1.0750
rate=4.0547
fn=0.2953
ans=0.2953
>> fzero(@rate_T,0.5)
ans=0.73107565193641
% rate_T2
function fn=rate_T2(y)
global beta gamma flowvol
c=y(1);
T=y(2);
rate=c*exp(gamma*(1-1/T));
fn(1)=flowvol*(1-c)-rate;
fn(2)=flowvol*(1-T)+beta*rate;
>> feval(@rate_T2,[0.5 1.1])
ans=-3.29556351331856 0.27683452699778
global beta gamma flowvol
beta=0.25; gamma=30; flowvol=25;
fsolve(@rate_T2,[0.0 1.25])
 
بالا