BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI VŨ HUY KHUÊ NGHIÊN CỨU QUÁ TRÌNH TRUYỀN NHIỆT TRUYỀN CHẤT VÀ CÁC GIẢI PHÁP TIẾT KIỆM NĂNG LƯỢNG TRONG LẠNH ĐÔNG CÁ THU Chuyên ngành: : Kỹ thuật nhiệt Mã số: 62520115 TÓM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT NHIỆT NGƯỜI HƯỚNG DẪN KHOA HỌC: TS: NGUYỄN XUÂN TIÊN GS.TSKH: TRẦN VĂN PHÚ Hà Nội – 2015 Hi st xa b mt nhng nhóm thy hi sn có giá tr kinh t cao nht, bao gm mt cá thu, cá ng trích, cá cng nm cách xa b, o qun nht cho nhóm ng làm lnh c bin ti nhà máy tái co qu nâng cao chng, gim hao ht cbo qun, vic xây dng b d liu v thành phn, v tính cht nhit vt lý ca sn phm ch cù ht sc cn thit S bii tính cht ca thc phm ct q trình lý hóa rt phc tp S gim nhi ca thc phm xung di nhi âm s khin dung dch lng bên thc phm chuyn thành dng rn Quá trình chuyn pha s i hồn tồn tính cht nhit vt lý ca sn phm Các nghiên cu th gii hin ch mi c tính cht nhit vt lý ca mt s sn phi vi sn phm khác, nhà nghiên cu ch ng di giá tr có tính chng vi sai s có th t nhng lý dc tìm kim d liu vào cho q trình tính tốn, thit k vn hành h thng lnh hp lý, góp phn nâng cao chng, gim hao ht sn phm tit king trình c gii quyt v t tác gi la chn ng tài lun án thơng qua: Nghiên cu q trình truyn nhit truyn cht, tìm kim gii pháp tit king lơng mt loi hi sn ph bin có giá tr cao cá thu - Vi vic la chn phát trin nh tính cht nhit vt lý b hc, ch thơng qua phân tích thành phn ca thc phng nghiên cu cá thu B n có th nh tính cht nhit vt lý ca thy sn theo lơ sn xut nhanh chóng - Xây dng mơ hình mơ phng q trình làm lnh cn ng ca truyn cht, kt qu mô ph vi kt qu thc nghim - Xây dng, hoàn thin phn mm nh tính cht nhit vt lý, mơ phnh thi gian ng nhi c sn xut ch bin có th s dng phn mm d dàng - Nghiên cu ng ca trình truyn cht trình cy sn, thc nghim hao ht khng ccá thu - Nghiên cu ng ca mt s yu t n trình c cp lý vi cá thu sn xut thc ti pháp tit king trình c 1 1.1 1.1.1 Nâng cao chng ch bin loc, thc phm sau thu hot mt thách thi vi toàn cu, bi quyt v c V c th hin hai mt, th nht phm bo cung c v s ng, th hai phi bm v chng Trong ch bin thc phm cc bit công ngh ch bin lnh rt nhiu hn ch bt cp T l tn tht sn phm sau thu hoch cao, ví d i vi rau qu khoi vi ch bin thy sn khong 20% Chng sn ph t khu b hn ch n 2001-a thy sn vào GDP chung tồn qung t 3,72%3,1% Bình qn thy sn gii quyt công c làm cho khoy sn cung cp thc phm cho 80 trii dân Viy sng khong t 39,3-42,86% tng sng thc phm, góp phn quan trng vim bo an ninh thc ph ng quc gia Nhy sng chuyn du kinh t ngành nông nghiu sn xut nông, lâm, thu sn chuyn dt, chng, hiu qu, giá tr gn vi th v kinh t, phát trin thy s sâu sc v an ninh quc phòng 1.1.2 Vit Nam c xp hng 10 nhà cung cp ln nht, sn phm thy sc xut khu quc gia toàn th gii Theo tng cc thng kê, Vit Nam gc nhiu thành công vi tng st 5.918,6 ngàn tn thu v 6,8 t USD [2] Cho ti nghiên cu mang tính h th v c Chúng ta tay v tính cht nhit vt lý ca thy sng loi thy sn n rt l hoàn thin công ngh làm lnh cy sng thm bo c hai tiêu chí (im bo chng ca thy sn (ii) s dng hng Mc dù thi gian gt mnh m c Tháng 6/2013 B Khoa hc Công ngh khánh thành phịng thí nghim cơng ngh CAS (Cell Alive System)Công ngh bo qung[10], bn cht CAS công ngh kt hp làm ln t kt hp v bo qun sn phng công ngh ng dng hiu qu Vit nam cn phi tip tc nghiên cu ng nhà máy ch biên thc phm cịn ln Cn có gii pháp tit kim ng c bit h thng c ht khi ng trình cc t b mt sn phng v phi nghiên c m, t l n c trình c nhà máy ch bin, thi gian cc xnh bng thc nghim kinh nghim Cn có nghiên c d ng lý thuyt ng quan trp công nghip m cao nên tht ca cá thu d i, chóng b phân hy có th gây ng c nm ca cá Thu khin cho yêu cu v k thut c nên rt kht nhng yu t quan trng vic la chn cá tng nghiên cu ca lun án S dng vt chng, xn s gia dn nhit khuch tán ng quát ca trình truyn nhit truyn chc vit nh[19]: t t C (1.6) D C q C gii h ) nêu trên, Tri tích di s Jo i, s dng s ca toán dn nhit toán dn m tng quát Vi cách tip cn này, mt s toán ca trình sy du vi t chc gii quyt 1.2.1.2 Q trình cc phm có nhit Trong q trình cy hing bin ti mt ph n mà tn tng ca hi pha không liên tc, mc thc phng chim t 70% Quá trình chuyn pha t c phm din khong nhi rng bu t nhi n gn nhi không tuyi Mt khác xi nhit liên hp gm (i) dn nhit không nh lòng thc phi nhia b mt thc phm ng làm lng nhi ca thc phm không u Dn ti mi mt phân t ca thc phm có mt nhi ng vt nhit vt lý (Cp(T); (T); bit Vì vy h 6) có tính phi tuyn rõ rt Vì nhng lý nêu trên, nghiên cu ki tốn mơ phng q trình i nhit-i cht cc phm s suy bin v n nhit khơng nh có ngun nhit bên có k i cht C (T ) (T ) T (r , ) (T) T( r, ) q( r , ) ( hi mi ) (1.7) Q trình cc phm ln ln gn lin vi trình truyn chu t quan trng nh ng trc tin ht khng chng thc phm Dn nhit không phi nhit nht, c thc hin thông qua vic khuch tán vt cht Thành phn c biu ding nhit vii dng entanpi dng rút gn c (1.8) [75] h t hw.mw (1.8) hw mw entanpi dòng khch tán ca thành phc bên thc phm c b mt thc phm u t c thc phm b c nh trình khuch tán bên s giing bng li m lúc di Nc t th rn b mt thc phm s ng c Hing ch din lp rt mng b mt thc phm, gi hing mc b mt sn ph Dehydration u ơng thnh g q trình dehydration Psat Ts m 1/ Pa (1.9) RT g w D u dày thoát m D h s khuch tán qua lp thoát w khng kg/kmol cc Trong khuôn kh nghiên cu ca lung nghiên cu cá thu v cp ng ca trình truyn cht khơng l mơ phng q trình truyn cht trình cc phm bng lý thuyt ca nhà khoa h u kin Vit Nam khó có th thc hic Do khơng th xây dc mơ hình cho h s khuch tán m hiu dng (Dw) ca vt liu ng ca truyn cht phng bán thc nghi nim enthalpy trung bình s c s d ng nhit n sinh s ht khng trình c trên, thc phc cn nhit truyn cht xy ng thi Tuy nhiên, trình truyn cht xy v thc bii vi thc phc t, cá), phn ln tác gi xut b qua ng ca trình truyn cht Lúc trình c thun túy trình dn nhit Khi b qua trình khuch tán vt cht bên b mt khi thc phc c cn nhit thun nhit phi tuyn vt cht ng chng, tính cht nhit vt lý ca mi theo nhi trình vi phân dn nhit (Carslaw and Jaeger, 1959) mơ t q trình lc phm vit cho phân t có to r , thm T(r, ) (1.10) C(T) (T) (T) T(r,τ) q(r, ) C(T) Nhit dung riêng ph thuc vào nhi ca thc phm, kJ/kg.K Khng riêng ph thuc vào nhi ca thc phm, kg/m3 H s dn nhit ph thuc vào nhi ca thc phm, W/m.K q(r Lng nhit ngun nhit sinh s chuyn pha cc phân t ph thuc vào ta thi gian, W/m3 Nhi ca phân t ph thuc vào t thi gian, K ng nhi ca thc phm l vic gii b xung u kiu ki bao gm u kiu kin biên[6]: u kiu: Cho bit phân b ng nhi ti thu T(r , τ = 0) = Tin (r ) (1.11) u kin biên: Trong toán làm lnh cc phm, ch yu gu kin biên loi hou kin biên liên hu kin biên loi 4) u kin biên long hp b mt thc phm tip xúc trc tip vng làm lnh quy lut truyn nhit gia b mc T = α T(rn ,τ)-Ta (τ) (1.12a) n n Bài toán làm lnh cc phm vu kin biên long gp png thit b ABF ( Air Blast freezer) IQF (Individual quickly Freezer) ng hc bit cu kin biên loi 3, gu kin biên loi hou kin biên liên hp, xy b mt tip xúc trc tip vi vt rn khác (gp p xúc) T1 T2 -λ1 = -λ2 n 1,n n 2,n (1.12b) Khi nghiên cu trình dn nhit bên thc phm, nhà nghiên cng rt nhiu gii tốn trên, song có th -λ c da gi thuyt trình chuyng, s chuyn pha gii phóng nhit n nhi ng thi tn ti b mt phân pha gia nhit vc ly hng s i tích gii toán lni ting nht ca Plank.R (1913) Plank Lf PD Tf Ta RD2 (1.13) dng nhiu gi thit gn kt qu tính tốn khơng xát vi thc t (sai s n 50-70% d i gian c phát trin nhng mơ hình thc nghim d dốn thi gian ci bt công thc ca Q.T.Pham (1986a) [22,70,72], Rewtov [22], Hung, Thompson [22] f (SPHH) u hn, d dàng lp trình gii tốn vi phân bc ho nm cng ng dc vi vt th có hình dáng phc tp vic i không gian r n t hu ht nht v mt toán h c áp dng cho vt th i phc tc tp nht , c v mt thut tốn v mt lp trình ng dng n t hu h xác tốn hc cao nht, cho phép gic toán rt nhiu chiu, phi tuyt tốn phc tp, khó thc hi tích hu hn, vm thun, có kh i quyc tốn vi phân phi tuyng hình hc phc tc la chn gii h n nhit thc phm c c trình bày chi tit phn tip theo ca lun án da theo thut toán ca Onhishenko V.P.[31] ), Vì t g 28], Tchigeov(1979) [92], Latyshev(1992) [49 Trong , Tchigeov (1979)[92] xice xwo 1,105 0,7138 ln(t f t 1) (1.32) Mơ hình tu tiên d t dung riêng thc phnh ca Sibel (1982) [89] nhiên mơ hình ch xét ng cn nhit dung riêng Nh i Leniger& Beverloo (1975) Charm (1978) [27,51] Schwartzberg (1976,1981), xui tin phát trin mơ hình xác nh nhit dung riêng Thc t m ka thc phm, dòng nhit ta t sn phm bao gm: nhit hi thay ng nhi nhit i dc d Trong lun án s dng mơ hình cng phái Trumak I.G-Onhishenko V.Pc xây d lý thuyt coi thc ph nhin cân bng có bin pha[31,33,66] 1.3.3 Entanpi Nghiên cu ca Chang Tao (1981) Chen (1985a) Schwartzberg (1976) nhng mô hình c xây dng t nghiên cu thc nghim Khi s dng nhc phm khác s có sai s khác n án entanpi s c tính theo công thnh nhit dung riêng hiu dng Ce(T) ca TrumakI.G-Onhishenko V.P[31,33] c trình bp theo Khng riêng ca thc phm ph thuc ch yu vào thành phn thc phm nhi Thc phm c coi mt h nhiu pha Trong lun án s dng mơ hình ca Choi & Okos (1986 nh khng riêng ca thc pht thành phn khng ca chúng: Hin nay, có nhinh h s dn nhit khác nh h s dn nhit Mơ hình nh h s dn nhit xu xác, mt s tác gi khác Levy (1981) phát trin mơ hình nh h s dn nhit t mơ hình ca Maxwell Trong khn kh lun án tác gi dùng mơ hình ca Vi (2000)[67 nh h s dn nhit c trình by p theo i vi thy sn ca Vi nghiên cu ng ca trình truyn cht q trình ci mơ hình cc xây d n nhit không nh Do vy vic nghiên cu m ng ca truyn cht rt cn thit t mi có th khc có th b qua ng ca truyn cht, ch n yu t truyn nhit trình cy sn hay khơng Mơ hình mơ phng q trình cc xây d h n nhit phi tuyn không nh kt hp vu kin biên loi ba tích hu hn vi thu hi t ca Trumak I.G- OnhishenkoV.P[31,33c s dng lu gii h t trình truyn nhit không nh xi vi cá thu trình c Trong nghiên cu xây dng mơ hình mơ phng tính cht nhit vt lý ca thc phm tác gi chn cách tip cn cng phái Chumak I.G & Onhishenko V.P.[31,32,33], d coi thc pht h nhin, cân bng có chuyn pha theo nhi Cách tip cn cho kt qu c lý thuyt cht ch tác gi trình by mơ hình nh thơng s nhit vt lý ca thc phm 1.5 1.5.1 1.5.2 - La chn phát trin mơ hình hồn chnh thơng s nhit vt lý ca cá thu di nhi rng (-400C ÷ 400C) - Nghiên cu trình truyn nhin ng ca truyn chc nghi hao ht ca sn phm trình cng thc nghim cá thu - Xây dng mơ hình ng nhi , thi gian c hp lý ng phái ChumakOnhishenko, gii trc tip toán dn nhit phi tuyo hàm riêng b tích hu hn, có kh ng dng cho vt th có hình dng bt k Kt qu nghiên cu lý thuyt c kim chng bng thc nghim vi ng nghiên cu cá thu - Áp dng mơ hình mơ phng, kt hp vi kho sát thc t ti nhà máy, nghiên cu ng ca mt s yu t ti trình ct s gii pháp nhm tit king trình c 1.5.3 Cá thu vùng bin Vit Nam 1.5.4 Nghiên cu lý thuyt kt hp vi nghiên cu thc nghim Nghiên cu lý thuyt truyn nhit có xét ti ng ca truyn cht trình cnh ng ca trình truyn ch l hao ht khng trình c Thu Phát trin phát trin mơ hình hoàn chnh ng vt lý gm: h s dn nhit dung riêng Cp(t), kha cá thu di nhi -400n 400 nhng nghiên cu lý thuyt mi nht ca tác gi khác th gii Xây dng phn mnh tính cht nhit vt lý di nhi -40 ÷ 400C Xây dng mơ hình mơ phng q trình ccá thu vu kin biên loi gii toán dn nhit phi tuyn khơng ng chng b tích hu hn Thc nghim so sánh kt qu tính tốn thơng qua mơ hình vi s liu thí nghim cá thu ti phịng thí nghim t xác c ng dng khác S dng mơ hình phn mm nghiên cu ng ca thơng s , t gió ti thi gian cc phm, t t s gii pháp tit king cho h thng lnh Ni dung ca lun án gm: - La chn phát trin hồn chnh thơng s nhit vt lý ca cá thu di nhi rng (-400C ÷ 400C) - Nghiên cu trình truyn nhin ng ca truyn chc nghi hao ht ca sn phm trình c ng thc nghim cá thu - Xây dng mô hình ng nhi , thi gian c hp lý ng phái ChumakOnhishenko, gii trc tip tốn dn nhit phi tuyo hàm riêng b tích hu hn, có kh ng dng cho vt th có hình dng bt k Kt qu nghiên cu lý thuyt c kim chng bng thc nghim vi ng cá thu - Nghiên cu ng ca mt s yu t ti trình ct s gii pháp nhm tit kim ng q trình c CÁ THU Qua phân tích tác gi thy rng vn mt s hn ch mô , t n kt qu tính tốn có sai s ln tác gi s dng nhng mơ hình m tin c nh tính cht nhit vt lý ca cá thu phc v cho vic nghiên cu xây dng mơ hình tốn hc mơ phng nhi thi gian c p theo 2.1 thu 2.1.1 c chìa khố ca s sng, dung mơi quan trng nht tồn b i sng sinh vt Bn cht ca trình k s chuyn pha cc t lng sang rn quyn kh o v hay phá hu cu trúc sinh vt Trong thc phm cha nhiu long hp thc phm kn bao gm: t - icec liên kt (xb), u) Nu coi tng u xwo: x wo = xu+ xice+ xb (2.1) 2.1.2 Công thc nh t l a Tchigeov (1979)[92] c Fikiin (1996)[37,38] i nhiu loi thc phm di nhi rng Chính vy, lun án la chn công thc Tchigeov làm công thnh t l c phm c 2.1.3 ng nhit dung riêng hiu dng Ce c Chumak- Onhishenko [33,93 n nh ng ca c nhit n nhit hin trình làm lnh c n nhit tng quát ca trình lc phm 1.7, q(r, τ) c coi lng nhit ngun nhit sinh s chuyn pha cc phân t ph thuc vào ta thi gian[85] q r, lim V xwo T V L T V 0 lim V Q V lim V xwo LT T 0 xwo T L T d d q r, xwo T L T d T r, dT (2.3) 2.3 8, ta : C( T) ( T) T r, C( T) ( T) T r, xwo L(T ) d dT (T ) xwo ( T) L(T ) T r, d T r, dT (2.4) (T ) T r , Nt: Ce (T ) C(T) xwo L( T) d , 2.4 tr thành: dT T(r, ) Ce T (T) (2.5) (T) T(r,τ) e(T) nhit dung riêng hiu dng n ng ca thành phn nhit n nhit hin trình làm lnh c 2.1.3.1 Khi c c thc ph nhit n chuyn pha, cơng thc tính nhit dung riêng hiu dc la Choi & Okos (1986)[30]: Ce ( T) Cp ( T) xi Ci i xwo Cw xp Cp xfa Cfa xas Cas x fi C fi xCHO CCHO (2.6) 2.1.3.2 Nc thc phm b thc phm ph thuc vào nhing thi biu thc nhit dung riêng hiu dng có xut hin thành phn nhit n chuyn pha cc Nhit dung riêng hiu dc xác nh theo công thc: Ce ( T) C( T) xwo L(T) d dT xsi Csi ( T) xwo (1 ) Cw ( T) xwo Cice ( T) xwo L( T) i d dT (2.11) 2.1.4 Entanpi Xut phát t công thnh nhit dung riêng hiu dng theo nhi 2.6 2.11, c xác nh t h( T) (2.14) Ce t dT 40 2.1.5 Trong thc phm có nhiu thành phn có cu trúc nn Nu xem xét h s dn nhit ca mt thành phn bt k vi thành phn li thơng qua h s dn nhit hiu dng h gm hai thành phn s dn nhit hiu dng ca thành phn l [66,67]: Li j i Vj (2.17) j j - h s dn nhit ca thành phn th j ph thuc vào nhi, theo bng 2.3 lun án Vj x j - thành phn th tích th j j nh h s dn nhit Ki ca h cu thành t hai phn t: thành phn th i ( i , Vi L i ,1 Vi )[66,67] Ki L i (1 ) m (1 m)(1 m1/3)(1 m1/3 ) m Vi - thành phn th tích ca phn cịn li; H s dn nhit hiu dng ca toàn h s bng[66,67] (2.18) i Li (T) i i K i Li Li i (2.19) y h s dn nhit hiu dng ca thc phc thc hic sau [66]: c tính tốn theo thuc mơ t nh h s dn nhit i T ca thành phn t th i theo bng 2.3 c 2: nh h s dn nhit hiu dng Li T ca h phn t bao quanh phn t i theo công thc 2.17 c 3nh h s dn nhit Ki ca h cu thành t hai phn t: thành phn th i ( i , Vi ) thành phn li ( L i,1 Vi ) theo 2.18 nh h s dn nhit hiu dng ca tồn h theo cơng thc 2.19 2.1.6 (2.20) xi i i - thành phn rng ca vi vi tht, cá ) x i thành phn khng th i ca thc phm i - khng riêng thành phn th i ca thc phc ly theo bng 2.4 lun án 2.2 the Phn mc vit bng ngơn ng lp trình Matlab vi giao din, d s d S dng mơ hình nh tính cht nhit vt lý ca thc phc tng hp mc 2.1 D liu vào thành phn khc, protein, m, tro, hydrocacbon, nhi n mm x xut d liu dng excel ho th ngày màm hình 2.3 tin cy cnh tính cht nhit vt lý, tác gi s tin hành so sánh kt qu tính tốn bng phn mm kt qu thc nghii vi loi thc phm khác có thành phn khng tính cht nhit vt lý c công b th gii ca nhà khoa hc khác 2.3 ng nghiên cu ca lun án cá thu Vit Nam Lồi cá ch có mt vài s lic nghim ti mt s nhi khác v tính cht nhit vt lý nên r tin cy Tác gi la chn mt s loi cá bic tính gn vi cá thu ki xác ca mơ hình xánh tính cht nhit vt lý Các loi cá bao gm: cá ng, cá hi, cá vc, cá tuyt Ngồi ra, tác gi cịn kim chng vi tht bị, thc phm n tính cht nhit vc công b mt s báo quc t ASHARE 2010 2.3.2 Hình 2.7: S t sai lch gia kt qu c nghim kt qu tính tốn lý thuyt nh (cao nht 4% vi tht bò 12% i vi cá tuyt cá v sai lch chp nhc cng nht s ng ca thành phn khu xung quanh giá tr trung bình Hình 2.10 Các s liu thc nghic hiu nm vùng nhi âm sâu Kt qu tính tốn ca mơ hình thc nghim g Sai s ti nhc mơ hình khong 1% Ti vùng lân cm nhi bu kt dung riêng hiu dt ngi gim xung nhanh chóng vùng t l t ngt dn n nhit c thc hin hing chuyn pha rt ln Sau nhi gi a t l c m d gim nhi, dn nhit dung riêng hiu dng gim dn Khi nhi thc phm nh bu kng 10oC, t l i rt chm theo nhi, ng ca n nhit chuyn pha gim dn nhii ch gim nhi entanpi n thu S i ca entanpi theo nhi n nhi n nhit nm lân ci nhi bu kn nhi âm sâu Sai s gia kt qu tính tốn mơ hình khong 3% 10 n Vi nhi lân ci nhi bu k s dn nhit ngt s xut hin ca thành phi 0oC, h s dn nhit cn gp 3,9 ln so vc dng lng) Sai s gia thc nghim tính tính tốn ln nht 7% Kt qu tính cht nhit vt lý theo nhi ca cá thu Thông số nhiệt vật lý cá thu dải nhiệt độ (t= -40 0C đến 400C) xwo= 72,02 %, xp = 15,82%, xfa = 6,48%, x as = 1,18%, x cho =4,50% t f =-2,2oC T, [K] t, [oC] Tỷ lệ nước đóng băng, % Khối lượng Hệ số dẫn nhiệt, riêng, [kg/m3] [W/mK] Enthalpy tính tốn, [kJ/kg] Nhiệt dung riêng, [kJ/kgK] 233.15 -40 0.937 932.015 1.393 2.3822077 243.15 -30 0.919 931.576 1.322 24.105 2.4321517 253.15 -20 0.889 931.475 1.245 51.545 2.860829 263.15 -10 0.813 933.668 1.129 90.189 4.8660383 268.15 -5 0.659 940.914 0.958 135.755 14.400095 269.15 -4 0.563 945.871 0.86 159.3 24.564457 270.15 -3 0.336 958.075 0.639 211.813 61.285237 271.15 -2 976.86 0.404 309.361 3.5606318 272.15 -1 976.789 0.405 313.142 3.5582641 273.15 976.707 0.407 316.928 3.5562848 283.15 10 975.306 0.423 355.307 3.5503504 293.15 20 973.057 0.436 394.461 3.5488326 303.15 30 970.134 0.448 432.899 3.5411427 313.15 40 966.622 0.457 471.38 3.5523684 2.4 K T nhng nh th tác gi có th khnh tính cht nhit vt lý theo nhi t ca s bii cng nhi i So sánh vi kt qu thc nghim, kt qu tính tốn theo mơ hình y Qua nhng s liu tính tốn thc nghi th nhn thy tính cht nhit v i rõ rt Khong nhi t 5on 40oC khong nhi thc phm khơng có bii pha vy tính cht nhit vc thc nghii d dàng, kt qu thc nghim vi xác, sai s khong 5% Trong khong nhi t -5on 5oC, tính cht nhit vt lý ca thc phm có s bit ngt Nguyên nhân khong nhi c thc phm bu chuyn pha t lng sang rn xut hin mng nhit gii phóng gi nhit n chuyn pha Nhit ngt, h s dn nhit n gp ln cc Trong khong nhi t -5on -40o C tính cht nhit vt lý ca thc phi theo quy lut, khong nhi c thc phy nhit n chuyn pha ca 11 c có mà ch cịn lng nhit hic gi gim nhi sn phm Tác gi a chn, phát trin xây dc mơ hình mơ phnh tính cht nhit vt lý ca cá thu t nhng mô hình c nhiu nhà khoa hc th gii kim chng ch tin cy cao nh tính cht nhit vt lý mà tác gi n, s liu vào ch thành phn thc phm có th d nh t sn xut ch bin thy hi sn y, s dng mơ hình phn mm có th nh tính cht nhit vt lý mt cách nhanh chóng cho tng lô sn phm, sn phm ca tng vùng mit Bng ngơn ng Matlab tác gi ng hồn thic phn mc tính cht nhit vt lý ca cá thu mà ch cn bit thành phn ca thc phm Phn mm có th ng dng tra cu bt k tính cht nhit vt lý ca thy hi sn bt k nhi khong nhi -400C ÷ 400C B s liu tính cht nhit vt lý ca cá thu ph thuc vào nhi bng mơ hình thông qua phm mm mà tác gi m b tin c d liu dùng cho nghiên cu xây dng mơ hình tốn xác nh thi gian cng nhi ctp theo QUÁ TRÌNH CÁ THU 3.1 cá thu i cht gia b mt thc phng xung quanh nhân t n ht khng n chng sn phm c n u trình cp c pha lng s ng Sau c bn m l mt làm cn tr i cht gia thc phng Lúc trình bc mt sn phm tr thành hichuyn pha rn T s kt hp nghiên cu lý thuyt thc nghi kho sát ng ca i nhit truyn chn i nhit gia thc phng trình ct t khng bng lý thuynh t l nhit n trình c phn sau 3.2 hu 3.2.2 Khi xây dng mơ hình tốn hc hing dn nhit bên thc phc ct s gi thit gc chp nhn: i nhit bên thc phm ch thun túy trình dn nhit, b qua hing khung ca thc phm B qua hii cht gia thc phm ng c 3.2.3 k & Onhishenko Chumak I.G & Onhishenko V.P xut vic s dng h s hình hc chuyn t mơ hình tốn chiu vit h t sang mơ hình toán mt chiu vit h t tng quát (x* = x ) H c mơ t q trình dn nhii v h t chic vi Ce T T r,0 T T r, f r ,r T R1 , T R 2, r r T r Т r, r , (3.8) (3.9) R1 , R2 T R 1, r T R2, r Tc1 Tc2 T R1, F1 F2 T R 2, 0, n nhit vit t tng quát 12 (3.10) (3.11) u kiu u kin biên loi vit cho b mt rng bên (nu có) b mt bên ngồi vc làm lnh R1 , R2 nh ng vi b mt bên b mt bên 12 s ta nhit ng vi b mt b mt ngồi H s hình hc Г= RF -1 V (3.12) nh ca khi thc phm, F V lt din tích th tích ca khi thc phm Vi khi th n, s nhn giá tr bng vi vt làm lnh có hình dng: tm phng vơ hn, khi tr dài vô hn, khi cu Trong thc t, h s i nhi phía ca khi thc phm có th sai khác Tuy nhiên, thc t không th c t cc b ca tm b mt khi thc phm, h s i nhit i b mc xem nh theo t trung bình ti tt c mt xung quanh khi thc phm Do h s ta nhing hp bng tu ki hai i xng, nên toán ch c xét khônnh ca vt H mô t c vit l Ce T T T 0, r T r r Т r, r , (3.17) r T R, T R, T r, T r, Tc r f r ,r T R, 0, 0, R thu 3.3.1 cá thu ri r tích hu hn, hai trc t không gian r thc chia thành nhng kho (3.18) ri ri r vi i N , 1, r R N τ j+1 = τ j + τ v r (3.19) c gc thi gian j+1 A B L e j-1 f h D r i-1 N ri C r i+1 r Hinh 3.1c: Biểu diễn phân tố theo tọa độ thời gian 3.3.2 t 13 i,j-1 , , 3.29, 3.30, 3.31 e , Ce , i,j-1 i,j-1 , x j =0 T i,0 = T = 1÷N j = j+1 Tj=A -1.B , 3.29, 3.30, 3.31 e , Ce , Ti,jvà T i,j i,j-1vài,j-1 Sai Tj =j j=A -1.B T0,j < Max(|T j-j|)