K thut tớnh toỏn ong cụng ngh húa hc Mc Lc Li núi u Phn I: C s v Matlab Chng I:Bt u V4 Matlab 1.1 Nhp d liu qua dũng lnh 1.2 S dng help trc tuyn 1.3 ng dn .7 1.4 Lu v ti d liu Chng 2: Cỏc cu trỳc c bn v bin 2.1 tớnh toỏn vi Matlab 2.2 Gii thiu dng d liu Chng 3: Tớnh toỏn vúi ma trn v vec t .10 3.1 Vec t 10 3.1.2 Vec t ct v s chuyn v 12 3.1.3 Nhn, chia v m ca vec t .12 3.2 M a n 14 3.2.1 Cỏc ma frn c bit 15 3.2.2 Xõy dng ma n v cỏch trớch cỏc phn t ca ma trn 16 3.2.3 Cỏc phộp toỏn vi ma trn 19 Chng 4: 22 4.1 th n gin 22 4.2 V ng 23 4.3 V mt 26 Chng 5: Cỏc lung iu khin 30 5.1 Cỏc toỏn t logic 30 5.2 Lnh find 30 5.3 Cu lnh if .32 5.4 Cu trỳc lp 34 Chng 6: Phng phỏp s 35 6.1 ng cong hi quy 35 6.2 Phộp ni suy 35 6.3 Giỏ tr ca hm s 36 6.3.1 Hm 'inline' 37 6.4 Phộp tớnh tớch phõn v tớch phõn 38 6.5 Tớnh toỏn s hc v cu trỳc iu khin 38 Chng 7: Vit chng trỡnh Matlab .39 7.1 Kch bn v hm m-fle 39 7.2 Hm m-file 39 7.2.1 Nhng bin c bit ca hm s 40 7.2.2 Bin a phng v bin ton cc 41 7.2.3 Tớnh giỏ tr hm s cỏch giỏn tip 41 K thut tớnh toỏn ong cụng ngh húa hc Chng 8: Vn bõn .43 8.1 Chui kớ t 43 8.2 Vo - d liu 44 8.2.1 Vo d iu t bn p h ớm 44 8.2.2 a d iu mn hỡnh 44 8.2.3 X lớfile bn .46 Phn II: ng dng ca Matlab 51 cụng ngh húa hc 51 Chng 1: Dn nhit v i lu 52 1.1 Dn nhit 52 1.1.1 Túm tt lớ thuyt 52 1.1.2 V d 53 1.2 i lu 56 1.2.1 Túm tt lớ thuyt 56 .2.2 V d 59 Chng 2: Bc x nhit v truyn nhit 64 2.1 Trao i nhit bc x gia cỏc vt 64 2.1.1 Túm tt lớ thuyt 64 2.1.2 V .64 2.2 Truyn nhit v thit b trao i nhit 67 2.2.1 Túm tt thuyt 67 2.2.2 V d 68 Chng 3: K thut tỏch cht 74 3.1 Cõn bng lng h i 74 3.1.1 Túm tt l thuyt 74 3.1.2 Vớ d 75 3.1.3 Bi 79 3.2 Phng phỏp McCabe Thiele 80 3.2.1 Túm tt lý thuyt 80 3.2.2 V d 83 Chng 4: K thut phn ng 89 4.1 Nhit ng hc 89 4.1.1 Túm tt thuyt 89 4.1.2 Vi d 91 4.2 Cỏc k thut tin hnh phn ng 94 4.2.1 Túm tt lý thuyt 94 4.2.2 V d 95 4.2.3Bềè 99 Ph lc 100 Ti liu tham kho 125 K thut tớnh toỏn ong cụng ngh húa hc Li núi u Cỏc bi toỏn cụng ngh húa hc thc t l rt phc tp, chỳng ta cn cú cụng c h gii nú Thc t cú rt nhiu phn mn, nhiờn, ong cun sỏch ny, chỳng tụi mun gii thiu n cỏc sinh viờn nghnh cụng ngh húa mt cụng c rt mnh, ú l matlab Mat lab l mt cụng c toỏn hc rt mnh, ngoi nú cũn h tr ngụn ng lp trỡnh bc vúi cu trỳc n gin, gn gi, d tip cn, cú th lnh hot gii quyt cỏc bi toỏn thc t t c yờu cu mong mun Do ú, vic s dng nú lm cụng c gii quyt cỏc bi toỏn k thut l rt t t Ni dung cun sỏch bao gm: Phn I: C s v Matlab nhm gii thiu cho bn c s lc v Matlab Phn II: n g dng Matlab ong cụng ngh húa hc gm cú cỏc chng: Dn nhit v i lu; bc x v truyn nhit; k thut tỏch cht; k thut phn n g Trong tng chng cú túm tt c s lớ thuyt, kốm theo mi chng ú thỡ cú cỏc vớ d v bi liờn quan i kốm, mi mt vớ d minh c gii bng tay sau ú c gii li bng matlab mt cỏch sinh ng v d hiu, cú tớnh tng quỏt cao a matlab vo cụng ngh húa l mc ớch ca chỳng tụi Tuy nhiờn, c biờn son ln u nờn khụng th ỏnh c nhng sai sút, rt mong nhn c úng gúp ca bn c cun sỏch tr lờn gn gi hn vi cỏc sinh viờn cụng ngh húa cng nh nhng ngi quan tõm n lmh vc n y K thut tớnh toỏn ong cụng ngh húa hc Phn I: C s v Matlab Matlab l cụng c toỏn hc, trc ht l cỏc phộp tớnh khoa hc, ngoi nú cũn cho phộp v cỏc th, tớnh toỏn ma trn, lm vic vi cỏc a thc v cỏc hm tớch phõn ng thi,chỳng ta cú th lp trỡnh cỏc chng trỡnh tớnh toỏn mt cỏch linh hot theo nhu cu ca mỡnh c bit Matlab t rt manh m vic thit k v tớnh toỏn vector v ma trn Khi chy Matlab thỡ xut hin mt hoc nhiu ca s trờn mn hỡnh, ú quan trng nht l Command window, nú cho phộp nhp lnh v hin th kt qu K thut tớnh toỏn ong cụng ngh húa hc Chng I:Bt u vi Matab 1.1 Nhp d liu qua dũng lnh *) Cỏc lnh s c thc hin sau n phỳn Enter Ket qu s c hin th trờn mn hỡnh nu mun Tuy nhiờn, nú ch thc hin cỏc dũng lnh l hp l Cú th tỡm hiu thờm t menu Help Vớ d : ằ + 7.5 ằ 18/4 ằ * Chỳ ý : cỏc du cỏch khụng quan ng Matlab *) Kt qu tớnh toỏn sờ c t ng lu vo bin ans Vớ d : ằ 14/4 ans = 3.5000 ằ ansA(-6) ans 5.4399e-04 Chỳng ta cú th nh ngha bin cho riờng mỡnh Vớ d: ằ a = 14/4 a 3.5000 ằ b = aA(-6) b= 5.4399e-04 D liu s c lu vo bin a v b *) Khi kt thỳc dũng lnh kt thỳc bng du thỡ s khụng hin th kt qu mn hnh, cú th kim chng bng hai biu thc sau : ằ + 7.5 ằ + 7.5; *) Cú th thc thi nhiu lnh cựng mt lỳc, cỏc cõu lnh ngn cỏch bi du phy hoc du chm phy >> sỡnpỡ/4'), cos(pi); sin(O) K thut tớnh toỏn ong cụng ngh húa hc ans 0.707ỡ ans Chỳ ý : cos(pi) khụng c in mn hnh *) Mc nh Matlab ch hin th ch s Do ú, mun cú kt qu chớnh xỏc hn thỡ ta dựng lnh form at long thỡ cú th hin th c ti 15 ch s ằ 312/56 ans = 5.5714 ằ format long ằ 312/56 ans 5.57142857142857 Cú th tỡm hiu thờm: >> help format *) Khi dũng quỏ di cú th dựng du ni cõu xung dũng sau >> sin(l) + sin(2) - sin(3) + sin(4) - sin(5) + sin(6) - sin(8) + sin(9) - sin(10) + sin (ll) - sin(12) ans 1.0357 *) Dũng kớ t t sau du % ch cú tỏc dng lm dũng chỳ thớch >> sin(pi) % dựng tớnh giỏ tr ca sin (pi) 1.2 S dng help trc tuyn Chỳng ta cú th dựng cu trỳc >> help Vd: help ops Hoc >> lookfor Vd : lookfor inverse K thut tớnh toỏn ong cụng ngh húa hc 1.3 ng dn Trong Matlab, cỏc lnh hay cỏc chng trỡnh c cha m-file, cỏi m c biờn son t file bn nhng lu vúi uụi m rng l .m M-file chi chy c kbi nú nm th mc current working directory 1.4 Lu v ti d liu Rt d dng lu hoc gi bin no ú, ta ch can click vo File- menu sau ú chn Save Workspace as hoc Load Workspace V d : save >> s i = sin(pi/4); ằ c l = cos(pi/4); c2 = COS(pi/2); ằ str= hello world; ằ save ằ save data ằ save numdata si, cl ằ save strdata sfr % õy l chui % lu ton b bin vi nh dng matlab.mat % ltu cỏc bin vi pth dng nh phõn data.mat % lu cỏc bin s sl v c l ti numdata.mat % hru bin chui S vo strdata.mat ằ save allcos.dat c* -ascii allcos.dat % saves cl,c2 in 8-digit ascii format to Cũn load ằ load ằ load data s i cl ằ load S- data % ti tt c cỏc bin tfile matlab.mat % chi ti riờng tng bin file data.mat % ti tt c cỏc bin t file strdata.mat K thut tớnh toỏn ong cụng ngh húa hc Chng 2: Cỏc cu trỳc c bn v bin 2.1 Tớnh toỏn vi Matlab Cú loi bin s hc: s nguyờn, s thc, s phc Matlab cũn h tr thờm nhng loi khụng phi l s: Inf (giỏ tr vụ cc, vớ d : 1/0 s cho kt qu nh vy) hoc NaN (khụng phi l s, l kt qu ca cỏc phộp tớnh nh 0/0 hoc vụ cựng - vụ cựng) Chỳng ta cú th vit cỏc biu thc trờn dũng lnh ca Matlab : ằ (23*17)/7 Kt qu s l : ans 55.857! Matlab cú toỏn t c bn l \ , /(chia phi hoc ỏi) v A(hm m) Chỳ ý s khỏc ca toỏn t chia trỏi v chia phi : >> 19/3 % thc hin phộp chia: 19/3 ans 6.3333 ằ 19\3, 3/19 % thc hin phộp chia: 3/19 ans = 0.1579 ans 0.1579 Nhng hm c bn v lng giỏc, m cú th tỡm hiu >> help elfim Bi tp: Xỏc nh giỏ tr cỏc biu thc di õy bng tay v bng Matlab so sỏnh Chỳ ý ti s khỏc gia chia trỏi v phi, v dựng lnh help tỡm hiu v cỏc hm: round, floor, ce il 2/2-3 - 5\4 - (5\4) - -4\9 - 2/5 + 7A2 - K thut tớnh toỏn ong cụng ngh húa hc 10/2\5 - + - 3A2/4 3A2A3 + round (6/9 + - 2)12 + floor (6/9 + - 2)/2 + ceil (6/9 + - 2)12 X = pi/3, x = x - l , x = x + 5, x = abs(x)/x 2.2 Gii thiu dng d liu Mc dự tt c cỏc tớnh toỏn s Matlab u c tớnh toỏn vi chớnh xỏc kộp (double precision), nhng khuụn dng ca d liu a cú th nhdng li nh cỏc lnh nh dng ca Matlab Cỏc bin ngm nh cng nhcỏc bin ca ngi s dng nh ngha u cú th a nhiu nh dng khỏc inh dng c chn nh s dng lnh form at tip n l ch nh long hay short v cui cựng cú th l dng biu din ca d liu, c th l e i vi dng hm m, f i dng du phy ng(ngm nh l f) Mt khuụn dng ớt dựng hn l rat a d liu di dng phõn s (chớnh xỏc hoc gn ỳng ) Vớ d: >> 4*atan(l) ans = 3.1416 ằ format longe ằ ans ans = 3.141592653589793e+000 ằ format lo n g f ằ ans ans 3.14159265358979 ằ format short e ằ ans ans 3.1416c 000 ằ format rat ằ ans ans 355/113 K thut tớnh toỏn ong cụng ngh húa hc Chng 3: Tớnh toỏn vúi matrn v vec t Nờn tng ca Matlab chớnh l ma trn (hay mng).Trng hp c bit: ma trn 1X1: l s vụ hng; ma trn ch cú dũng hoc ct: l vec t 3.1 Vec t Vec t hng l danh sỏch cỏc s c ngn cỏch bi du cỏch hoc du phy S lng phn t ca vec t c gi l length ca vec t Cỏc phn t nhp vo phi t ong du [ ] Vd: >> V = [-1 sm(3) 7] V= -1.0000 0.1411 7.0000 ằ length(v) ans Mt s phộp toỏn cú th thc hin vi vec t nh nhõn vụ hng, cng hoc tr t cỏc vec t khỏc cú s phn t tng t, hoc gia s vi vec t Mi phộp toỏn u c thc hin trcn tng phn t Cỏc vec t cng cú th c xõy dng t cỏc phn t hu hn ằ V = [-1 7J; w = [2 4]; ằ z - V+ w % cng tng phn t 511 ằ vv = V + % cng vo tt c cỏc phn t vo vec t vv = 149 ằ t = [2*v, -wj ans -2 14 -2 -3 -4 ng thi phn t cú th c trớch dn ra, hoc chinh sa: ằ v(2) = -1 % sa phn t th ca V -1 -1 ằ w(2) % hin th phn t th ca w 10 K thut tớnh toỏn ong cụng ngh húa hc % k : cỏc hng s ca phng trỡnh antoine (nhp dng ma trn) %t : nhit ca h ( c ) V du : a = antoine([6.90565 1211.033 220.79], 65) = 2.6682 p = 10 = 168.82 Muc ớch : tớnh ỏp sut riờng hi bóo hũa ca mt cu t nguyờn cht da vo phng trỡnh antione: lg(/>O) = A - nhit ang x ộ t c +6 5.3.2 Hm McCabeThiele.m õy thc cht l mt script file ch khụng phi l mt hm Code : fonction McCabeThiele %chuc nang cua chuong trinh la xay dung gian bang phuong phap %McCabeThiele va tinh so dia ly thuyet theo phuong phap global fL X hv Cp R y alpha alpha = input('cho biet he so alpha:'); xd = input('hay nhap gia tri cua xd:'); xb = input('hay nhap gia tri cua xb:'); xf = input('hay nhap gia tri cua xf:'); R = input('Ty so hoi luu R:'); t = inputCnhiet soi hoac nhiet ngung tu (tuy dieu ken):'); tf = input('nhiet cua nguyen lieu dau vao:'); X = [x f l-x f]; hv = input('nhap ma tran nhiet hoa hoi tuong ung cua cac phan:'); Cp = mput('nhiet dung rieng cua hon hop fL = mputCphan mol long(chi nhap nguyen lieu hoan toan dang hoi):'); dispC -') ớprmtf('chon nhung chuoi:\n') ớprmtf('long qua lanh,diem soi,mot phan hoi,diem ngung tu,nhiet qua soi.\n') d is p C . - ') dk = inputCNguyen lieu dau vao ong dieu kien:'); 111 K thut tớnh toỏn ong cụng ngh húa hc %phuong trinh duong cong can bang for i= l:l 10 y=0.01*(i-l); ye(i)=0.01*(i-l); xe(i)=fzero('eqliq',0.5); end q = q_line(dk,t,tớ); %toa cua diem cat yi=(xf+xd*q/R)/(l+q/R); xi=(-(q-l )*(1 -R/(R+1 ))*xd-xớ)/((q-l )*R/(R+1 )-q); figure(l); hold on; AXIS([0 01]); %thiet lap khung thi plot(xe,ye,V); set(line([0 1],[0 1]),'Color',[0 0]); set(line([xd xi],[xd yi]),'Color',[1 1]); set(line( [xf xi], [xf yi]),'Color', [1 1]); set(line([xb xi],[xb yi]),'Color',[1 1]); %bo phan en cot cat i=l; xp(l)=xd; yp(l)=xd; y=xd; while (xp(i)>xi), xp(i+l )=fzero('eqliq',0.5); yp(i+l )=R/(R+1 )*xp(i+l )+xd/(R+l); y=yp(i+l); set(line([xp(i) xp(i+l)],[yp(i) yp(i)]),Color', [0 1]); if (xp(i+l)>xi) set(line([xp(i+l) xp(i+l)],[yp(i) yp(i+l)]),'Color',[0 1]); end i= i+ l; 112 K thut tớnh toỏn ong cụng ngh húa hc end %bo phan day cot cat ss=(yi-xb)/(xi-xb); yp(i)=ss*(xp(i)-xb)+xb; y=yp(i); set(line([xp(i) xp(i)],[yp(i-l) yp(i)]),'Color',[0 1]); while (xp(i)>xb), xp(i+l)=feero('eqliq',0.5); yp(i+l)=ss*(xp(i+l)-xb)+xb; y=yp(i+l); set(line([xp(i) xp(i+l)],[yp(i) yp(i)]),'Color',[0 1]); if (xp(i+l)>xb) set(line([xp(i+l) xp(i+l)],[yp(i) yp(i+l)]),'Color',[0 1]); end i= i+ l; end fprintf(' - \n') fprintf('So dia ly thuyet tinh theo phuong phap la: %f\n',i) hold off; Muc ớch : xõy dng gin h hai cu t theo phng phỏp McCabeTheile v tớnh toỏn s a lớ thuyt trờn ú b i t: h s a , xd , xb , x f , nhit húa hi ca cỏc thnh phn nguyờn cht cú pha lng hv, nhit dung riờng ca hn hrp Cp, fL nu nguyờn liu cú phn dng hi (thng l khụng cú nờn khụng nhp) , v iu kin u vo ca nguyờn liu (nh nờu lớ thuyt ca chng 5) Vi du : a = 3.2, xd = 0.97, xb = 0.02 , xf = 0.4 , hv = [3.08e4 1,6e2 , iu kin ca nguyờn liu u vo l : 'long qua lanh' 3.33e4] , Cp = Ta thu c kt qu : s a lớ thuyt l : 13 hin th ờn mn hỡnh Command Window v gin : 113 K thut tớnh toỏn cụng ngh húa hc 5.3.3 Hm eqlớq.m C ode: function f=eqliq(x) global y alpha f=y-alpha*x/( 1+(alpha-l )*x); C phỏp : y = eqliq(x) %x: thnh phn pha lng Muc ớch: xõy dng ng cõn bng lng hi da vo phng trỡnh 5.23 Do , kh bit thnh phn pha hi ( , ) thỡ cú th tỡm c thnh phn pha lng cõn bng X; = fzeroCeqliq'jO.S) 5.3.4 Hm Fenske.m C ode: function n = Fenske(xD,yD,xB,yB) %tinh so dia ly thuyet bang phuong tnnh Fenske n = log(xD*(l -xB)/(xB*(l -xD)))/log(sqrt(alpha(xD,yD)*alpha(xB,yB)))1; function y = alpha(x,y) y = y*(l-x)/(x*(y-l)); 114 K thut tớnh toỏn ong cụng ngh húa hc return Cỳ vhỏp : y = Fenske(xd,yd,xb,yb) %xd: Phn mol ca thnh phn A (A thng l cu t nh)trong phn ct %yd: phn hi cõn bng vi xd %xb: Phn moi ca thnh phn A phn luyn ỏy %yb: phn hi cõn bng vi xb V du : n = Fenske(0.95 , 0.98 , 0.05 , 0.11) = 5.5393 Muc ớch : dựng tớnh s a lý thuyt ca ct ct da vo phng trỡnh Fenske ( phng trỡnh 5.22) 5.3.5 Hm Heat.m Code : function [Qc,Qb] = Heat(D,td,B,tb,tf) %D:phan tren cot cat %td:nhiet cua phan en cot cat %B:phan day cot cat %tb:nhiet phan day cot cat global R Cp Xhv %R:ty so hoi luu %Cpl:nhiet dung rieng dang ap %x:ma tran phan mol cua cac phan long %hv:ma tran nhiet hoa hoi tuong ung cua cac phan %tuong ung Qc = D *(R+1)*lamda(x,hv); %Qc:cong suat bo ngung tu Qb = D*Cp*(td-tf)+B*Cp*(tb-tf)+Qc; %Qb:cong suat cua noi hoi Cỳ phỏp : [Qc,Qb] = Heat(D,td,B,tb,tf) %Qc : nhit lng ca thit b ngng t %Qb : cụng sut ca thit b cung cp nhit t núng ong quỏ trỡnh chng ct %D : lng sn phm trờn ct ct %B : lng sn phm phn luyn %td : nhit trờn thit b ngng t %tb : nhit cung cp luyn % tf: nhit nguyờn liu u vo K thut tớnh toỏn cụng ngh húa hc Vi u : D = 23.16 D = 355* , tB = 382.8T, tF = 293K kmol/h , B = 34.71kmol/h, Chỳ ý: hm ny ngoi cỏc bin a phng thỡ cũn cú bin cc b c khai bỏo sau t khúa gLobal: R Cp X hv hm ny chy c thỡ trc gi hm thỡ khai bo bin ton cc trc, v kh khai bỏo cng phi cú t khúa global Cho : R = 1.92 , Cp = 1.6x10 kJ/mol.K,X = [0.4 ],hv = [3.08xl04 3.33xl04] kJ/km ol: [Qc,Qb] = Heat(23.16,355,34.71, 382.8,293) = [2.1844xl06 2.9128xl06] kw Nh h ỡn h : Hỡnh 5.3 Muc ớch : Tớnh cụng sut ca b chng ct v cụng sut có lũ luyn cn thit thu c sn phm theo yờu c u 5.3.6 Hm lamda.m C ode: 116 K thut tớnh toỏn cụng ngh húa hc ^ t:\ M A I L A B \h i e u ld i s l i ll \l a m d a m File Edit D Text Window Hep Z ọr f u n c t i o n y ằ ln a a a (K ,liv j svua=0; n = le n c ii(x ) J f 01 i * l : n sum - S i m end Y - S im I ttu m ; Hỡnh 7A Cỳ vhv : y = lamda(x, hv) Trong ú : % : nhit húa hi ca hn hp lng % X : thnh phn ca cu t pha lng % h V : nhit húa hi cựa cu t i tng ng Vi d u : V = lamda(r0.4 0.6], [3.08xl04 3.33 xio4] ) = 3.23 x io kJ/mol Muc ớch : tớnh nhit húa hi ca hn hp lng nhiu thnh phn theo cụng thc : = txi-hvi /=1 5.3.7 Hm material.m C ode: function [D,B] = material(F,xfjxd,xb) %F phai CO don vi la:(kmol/h) a=[l xd xb]; b = [F xf*F]'; x = a\b; D = x(l); B = x(2); C phỏp : [D,B] = material(F,xf,xd,xb) %D : lng cht trờn ct ct(pin chung) %B : lng cht phn ỏy ct ct( phn luyn) K thut tớnh toỏn cụng ngh húa hc %F : lu lng dũng nguyờn liu u vo % x f : thnh phn phn mol ca cu t A (cu t nh)trong hn hp nguyờn liu vo %xd : thnh phn phn mol ca cu t A phn chng %xb : thnh phn phn mol ca cu t A phn luyn V du : [D , B] = material(57,87,0.4, 0.97,0.02) = [23.148 34.722] Nh hmh5.5 Muc ich : tớnh lng ca cỏc cu t phn chung v phn luyn r V MATIAB r* E* WwƠ Web D ó r B ^ @ m E: \H*T ~ Q n m ấ3 a - A l n i e s i đ Wmrlnwj IWp ? Crrrt Direetry E >MATLADHiô ijVỡỡsđ \ 1Camman >rtndow M 1i i ằ V đ l CD , B] = m a c e r ia l ,8 J , , 0.2) = H-Ci ự n to ln i t ằ [7 b t i l ằ n -c i tt-C i [T x n -c i [ J h t4 n n - i [] e q l q a H-Ci [l r e n s n e s H-Cl [jjH eac in 2 ằ n -E l landa.n W-Ci m auei.ial.iu H-Cl [H M cC abeThlele.E < B - H - t l q _ in e a B -f i, v p e q il.n i n -C i > currntDiredory _ Hỡnh 5.5 H m q ^ n e m C ode: nh ong hỡnh 5.6 Cỳ phỏp : [q,q_slope] = q_line(dk, t,tf) %q: l t s gia nhit chuyn mol nguyờn liu thnh hi bóo hũa vi nhit húa hi ca mol nguyờn liu Xem li chng %q_slope: h s gúc ca ng q % dk: iu kin ca nguyờn liu u vo (xem bng di phng trỡnh 5.19) % t: nhit sụi nu hn hp nguyờn liu vo l hn hp lụng nhit thp , v l nhit ngtmg t nu hn hp nguyờn liu ua vo iu kin quỏ sụi % tf: nhit ca nguyờn liu u vo K thut tớnh toỏn cụng ngh húa hc r V E:VMATLAB\hieuWisIill\q ire.m Fis Edit Text Window Help D c1 ô 10 11 12 13 14 15 16 1? 13 19 20 21 22 23 24 25 2ô 27 @đ E m õ * f fu n c tio n fq , c[_sJLope! - q _ lin e (d k ,t.,t;f ) ichuc nang; ciHh ỗ i a c r i < tro n g cac cruong hop g lo b a l L X hv Cp *ớ!.:p !w n o i long ùla n d a :n h ie c hos hoi cua mDl *Cp; co Hi*3 la Cpl tioaỗ Cpu vcpl : Sỡlỡaat un cua Chat lo n ; *L'pv;rtfiiet dng cua pha hoi sw itch di case ('ioncr qua la n h ') q - + C p*(c - tÊ ) /la f ld a ( x ,h v ) ; Sct:o3U5t Eo 1.39l8e+ 00S ằ > Current DlrectV - s ta rt'I Hỡnh 5.8 5.4.3 Hm BRk.m C ode: nh ong ni dung hỡnh 5.9 % x a : chuyn húa ca phn ng % t: thi gian m phn ỳng t c n chuyn húa xa %n : bc ca phn ng ( õy tụi mi ch thit lp cho phn ng bc v bc 2) % const: t l ca h s phn ng (A+2B -> C)thỡ const = ,Gi s cho Tựy thuc vo phn ng m chnh const cho phự hp Ky thuat tinh toin cong nghp h6a hpc %quad8: hm tớnh tớch phõn C phỏp : k = BRk(xa) E:\MATLAB\hieu\kirieticsWlkm File Edit Text ? (3 đ Window Help % m V d M (2) fu n c tio n y = BRkJxaJ g lo b a l t n co n st i f n ==1 S = quads(@ ylr 0#x a , le - J / t ; elsexE n ==2 y = wadS (@y2,0xx a , l- - ) / t ; end fu n c tio n y = y l(K ) 10 11 Y = l / ( l- x ) ; fu n c tio n y = y2 -H C u rren! it tlu iy iaart Hỡnh 5.12 124 K thut tớnh toỏn ong cụng ngh húa hc Ti liu tham kho 1) Dilip K Das , P E , Dr Rajaram K Prabhudesai, P E , EIT Chemical Review, 2s edn Engineering Press : USA, 1998 2) Bruce A Finlayson , Introduction to chemical engineering computing John Wilay & Sons , Inc., 2006 3) Ela P, ekalska, Maijolein van der Glas , Introduction to Matlab Pattern Recognition Group, Faculty o f Applied Sciences Delft University o f Technology , January 2002 4) http://www.mathworks.com/matlabcentral/link exchange/MATLAB/Academic Curri cula/Engineering/Chemical Engineering/ 5) http://www.che.ncsu.edu/current students/matlab che.pdf 6) http://www.owlnet.rice.edu/~chbe403/bdistmccabe.html 125 ... 28 Kỹ thuật tính toán công nghệ hóa học Điểm khác hai lệnh surf mesh là: surf tô màu bề mặt,còn mesh không 29 Kỹ thuật tính toán ừong công nghệ hóa học Chương 5: Các ỉuồng điều khiển 5.1 Các toán. .. quan tâm đến lmh vực n y Kỹ thuật tính toán ừong công nghệ hóa học Phần I: Cơ sở Matlab Matlab công cụ toán học, trước hết phép tính khoa học, cho phép vẽ đồ thị, tính toán ma trận, làm việc với... 100 Tài liệu tham khảo 125 Kỹ thuật tính toán ừong công nghệ hóa học Lời nói đầu Các toán công nghệ hóa học thực tế phức tạp, cần có công cụ hỗ ừợ để giải Thực tế có nhiều phần