các yếu tố ảnh hưởng thời gian cấp đông

KHOA HỌC CÔNG NGHỆ NGHIÊN CỨU ĐÁNH GIÁ CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN THỜI GIAN CẤP ĐÔNG TRONG QUÁ TRÌNH CẤP ĐƠNG CÁ TRA FILLET Đỗ Hữu Hồng 1, Huỳnh Bảo Long 2, Hoàng Thị Nam Hương 1,2 Trường Đại học Công nghiệp Thực PhẩmTP.HCM, 3Trường Đại học Bách khoa Tp HCM Ngày gửi bài: 12/5/2016 Ngày chấp nhận đăng: 06/6/2016 TĨM TẮT: Bài báo trình bày kết nghiên cứu đánh giá yếu tố ảnh hưởng đến thời gian cấp đông bao gồm: vận tốc không khí, nhiệt độ khơng khí mơi trường cấp đơng chiều dày sản phẩm, từ xây dựng phương trình hồi quy thiết lập mối quan hệ thời gian cấp đông với yếu tố ảnh hưởng làm nền tảng cho việc chọn chế độ cấp đông hợp lý q trình cấp đơng cá tra fillet Từ khóa: thời gian cấp đơng, vận tớc khơng khí, nhiệt động khơng khí, chiều dày sản phẩm, phương trình hồi qui Phương pháp phần tử hữu hạn ABSTRACT This paper presents the results of studies assessing the factors affecting the freezing time includes: air velocity, air temperature and product thickness, thereby building the regression equation establish relationships between freezing time with the factors affecting the foundation for selecting appropriate freezing regime of the frozen pangasius fillet Keywords: freezing time, air velocity, air temperature, product thickness, the regression equation The element method ĐẶT VẤN ĐỀ Cá tra mặt hàng xuất khẩu chủ lực nước ta Tuy nhiên, năm gần sản lượng khai thác liên tục tăng tổng kim ngạch xuất khẩu không tăng tương ứng số sản phẩm phải bán hạ giá khơng đạt u cầu chất lượng Một yếu tố chủ yếu định chất lượng sản phẩm thủy, hải sản sau đánh bắt công nghệ kết đông bảo quản lạnh Về nguyên tắc tốc độ cấp đông nhanh, thời gian cấp đơng ngắn chất lượng sản phẩm tớt thời gian bảo quản được kéo dài Tuy nhiên điều dẫn tới nhiệt độ môi trường cấp đông phải thấp, công suất hệ thống lạnh tăng lên, hiệu suất lượng giảm đi, giá thành sản phẩm tăng cao [6,15,16] Do đó, nghiên cứu đánh giá yếu tố ảnh hưởng đến thời gian cấp đông cá tra fillet có ý nghĩa quan trọng việc đánh giá chất lượng sản phẩm, công nghệ cấp đơng tiêu hao lượng q trình làm lạnh cấp đông NỘI DUNG NGHIÊN CỨU Quá trình cấp đơng q trình phức tạp, số tượng xảy đồng thời trình truyền nhiệt, truyền chất, trình biến đổi pha, sư thay đổi tính…Đặc biệt điểm kết đơng, tính chất nhiệt vật lý cá (nhiệt dung riêng C, hệ số dẫn nhiệt , khối lượng riêng) biến đổi đột ngột, phương trình vi phân mơ tả q trình phi tuyến trở nên vơ khó giải Cho đến nhiều tác giả sử dụng phương pháp gần sai phân hữu hạn, phần tử hữu hạn với cách tính lặp để giải toán biến đổi pha theo sơ đồ biên giới phân pha cố định hay di động, xác định thời gian cấp đông theo phương pháp giải tích, nhiên, độ xác khơng cao hầu hết chỉ dừng ở toán chiều [1-3, 10-13] TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 41 KHOA HỌC CÔNG NGHỆ Trong báo này, tác giả sử dụng phương pháp phần tử hữu hạn để giải tốn dẫn nhiệt có biến đổi pha để xác định thời gian cấp đông ở chế độ Trên sở thời gian cấp đông được xác định, tác giả tiến hành đánh giá chi tiết yếu tố ảnh hưởng đến thời gian cấp đơng xây dựng phương trình hồi quy thiết lập mối quan hệ thời gian cấp đông với yếu tố ảnh hưởng 2.1 Đối tượng phạm vi nghiên cứu Đối tượng khảo sát mẫu sản phẩm cá tra dạng hình hộp có kích thước bảng [15] Bảng 1: Thơng sớ hình học cá tra fillet STT Khối lượng (g) Dày Rộng (mm) W(mm) 85,5124,5 10-14 60-80 124,5199,5 14-16 80-90 199,5256,5 16-18 Dài L (mm) 150200 200- 250 90-110 250270 1102704 256,5-313,5 18-21 120 290 1202805 320 24 130 300 Cá tra thực tế được cấp đông băng chuyền IQF dạng thẳng với miền giới hạn thông số làm việc sau: Bảng 2: Miền giới hạn thơng sớ làm việc q trình cấp đơng cá tra fillet băng chuyền IQF [15] STT Thơng sớ Giới hạn Nhiệt độ khơng khí te (-451 (0C) 35) Vận tớc khơng khí  (515) (m/s) Áp dụng kết tính chất nhiệt vật lý [7] được sử để mơ tốn dẫn nhiệt khơng ổn định nghiên cứu Hệ số tỏa nhiệt bề mặt được tính theo [10]    25  0,6 W m2 K  (2) 2.2 Phương pháp nghiên cứu 2.2.1 Xây dựng mơ hình tốn cho tốn cấp đơng Do mẫu sản phẩm cá có kích thước dài lớn nhiều bề rộng dày, truyền nhiệt dọc theo chiều dài sản phẩm nhỏ có thể bỏ qua, nên coi nhiệt độ chỉ thay đổi theo hai hướng bề rộng (hướng x) bề dày (hướng y): T = f(x,y,) ( hình 1) TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 42 KHOA HỌC CƠNG NGHỆ Tỏa nhiệt Khơng khí y Khơng khí Tỏa nhiệt x Tỏa nhiệt Tỏa nhiệt Khơng khí Khơng khí Hình Mơ tả sơ đồ vị trí xếp đặt sản phẩm Bài tốn trường hợp được mơ tả bởi phương trình vi phân chủ đạo điều kiện sau:    2T  2T  T (T)   q V  (T)C P (T) y    x   T   T         x  x 0  y  y 0     T  T  Tf x  x       T   x  x   x   T      T  Tf y  y  T    y  y    y     T  Tx, y,0  const   1 2.2.2 Phương pháp giải Khi áp dụng phương pháp PTHH, phần tử được chọn hai chiều chữ nhật bậc nút, lưới được chia theo cấu trúc đều Áp dụng phương pháp Galerkin [14] Phương trình ma trận đặc trưng PTHH cho phần tử khơng có nguồn xạ có dạng tổng quát sau: C T   KT  f  (3)    Trong đó: C    .c P N T N dV ma trận nhiệt dung V K   B DBdV   hN NdS ma trận nhiệt dẫn T V  f    hT K N T dS T S véc tơ phụ tải nhiệt S 2.2.3 Kết Bằng phần mềm Ansys: Q trình tính lặp thể nghiệm hội tụ tuyệt đới, hình TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 43 KHOA HỌC CƠNG NGHỆ Hình Q trình lặp giải toán truyền nhiệt bằng phầm mềm Ansys Biến thiên nhiệt độ cá tra fillet theo thời gian (nhiệt độ tam sản phẩm) được thể hình Phân bố nhiệt độ tâm sản phẩm theo tiết diện ngang Nhiệt độ sản phẩm t (0C) 10m/s, t=-400C Thời gian cấp đông (s) g sản u rộn Chiề (m phẩm m) Hình biến thiên nhiệt độ tâm sản phẩm Quá trình biến đổi pha mẫu sản phẩm được thể hình (a) (b) (c) Hình Quá trình hình thành phát triển pha rắn mẫu sản phẩm cá a) Bắt đầu hình thành pha rắn mặt ngồi (màu xanh lá), bên chưa đơng đặc (màu đỏ), vùng biến đổi pha (vàng) (b) Pha rắn phát triển, biên giới phân pha di chuyển dần vào sâu vật (c) Pha rắn chiếm toàn vật Cá đơng lạnh hồn tồn Theo [8], với sớ liệu chọn ở trên, sai số thời gian cấp đông lý thuyết thực nghiệm chưa đến 3% sai số trường nhiệt độ chưa đến 1,9(K) Do phương pháp hồn tồn có thể chấp nhận được để giải tốn dẫn nhiệt khơng ổn định cho nhiều chế độ cấp đông khác Kết xác định thời gian cấp đơng được trình bày bảng 10 2.3 Đánh giá các yếu tố ảnh hưởng đến thời gian cấp đơng Đới với q trình cấp đông cá tra fillet, thông số ảnh hưởng trực tiếp đến thời gian cấp đông công suất điện hệ thớng cấp đơng, là: nhiệt độ khơng khí mơi trường cấp đơng, vận tớc khơng khí mơi trường cấp đơng, chiều dày sản phẩm tính chất nhiệt vật lý cá tra Chính vậy phải xây dựng phương trình hồi qui thiết lập mới quan hệ yếu tố ảnh hưởng đến thời gian cấp đơng dựa kết trình bày ở bảng 10 Thực tế khảo sát nhà máy kích thước cá fillet chủ yếu tập trung ở khối lượng M =(57)oz tương ứng với kích thước (15x85)mm Cho nên phần nghiên cứu ảnh hưởng vận tốc nhiệt độ khơng khí đến thời gian cấp đơng suất tiêu hao lượng tác giả chọn kích thước (15x85)mm để đánh giá 2.3.1 Ảnh hưởng vận tốc khơng khí đến thời gian cấp đơng TẠP CHÍ KHOA HỌC CÔNG NGHỆ & THỰC PHẨM SỐ 09/2016 44 KHOA HỌC CÔNG NGHỆ Để đánh giá cách chi tiết việc thay đổi tốc độ ảnh hưởng đến thời gian cấp đông, tác giả thực hồi quy mối quan hệ vận tớc khơng khí với thời gian cấp đông khoảng  =(515)m/s ở chế độ nhiệt độ khác khoảng te=(-45-35)0C, kết hồi quy được trình bày bảng hình Bảng 3: Phương trình hồi quy thể mới quan hệ vận tớc khơng khí đến thời gian cấp đơng Nhiệt độ Phương trình hồi quy = 2269-169,949. + 5.29143. =2100-158,423.+ 4,93714.2 te= -35 C te= -37,0C = 1951,6-148,023.+ 4.61714. =1825,6-139,337. + te= -400C te= -42,5 C 4.34286. =1713,2-131,646. + 4,11429.2 te= -450C Quan hệ vận tớc khơng khí với thời gian cấp đơng ở chế độ nhiệt độ khác te=-350C te=-37,50C te=-400C te=-42,50C te=-450C T hờ i gian cấp đông, (s) 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 10 11 12 Vận tớc khơng khí , (m/s) 13 14 15 Hình 5: Quan hệ thời gian cấp đông với vận tốc khơng khí Ta thấy tăng tớc độ khơng khí làm thời gian đông kết giảm đáng kể Độ giảm thời gian ứng với độ tăng m/s ( =1m/s) được ghi bảng Bảng 4: Độ giảm thời gian cấp đông (s) ứng với  =1m/s ở chế độ khác Nhiệt độ te= -350C te = 00 CC te 37,5 = -40 te = 00 CC te 42,5 = -45 Vận tớc khơng khí, (m/s) 7,5 10 12, 15 -117 -91 -64 -38 -11 -109 -84 -60 -35 -10 -102 -79 -56 -33 -10 -96 -74 -52 -31 -9 -91 -70 -49 -29 -8 Nhận xét:  Với tốc độ khoảng (5 15) m/s, tăng tớc độ khơng khí lạnh làm thời gian cấp đông giảm đáng kể Tuy nhiên độ giảm thời gian cấp đông không đồng đều chế độ nhiệt độ với chế độ vận tốc khác nhau, cụ thể sau:  Trong khoảng vận tốc =(57,5)m/s, tăng vận tốc, độ giảm thời gian trung bình 7,5% vận tớc khơng khí tăng =1m/s TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 45 KHOA HỌC CÔNG NGHỆ  Trong khoảng vận tốc =(7,510) m/s, tăng vận tớc, độ giảm thời gian trung bình 5,5% vận tớc khơng khí tăng =1m/s  Trong khoảng vận tớc =(1012,5) m/s, tăng vận tốc, độ giảm thời gian trung bình 3,9% vận tớc khơng khí tăng =1m/s  Trong khoảng vận tốc =(12,515) m/s, tăng vận tớc, độ giảm thời gian trung bình 2% vận tớc khơng khí tăng =1m/s 2.3.2 Ảnh hưởng nhiệt độ khơng khí đến thời gian cấp đơng Tương tự phần đánh giá ảnh hưởng vận tớc khơng khí đến thời gian cấp đông cá tra fillet, tác giả sử dụng kết xác định thời gian cấp đông ở bảng 10 để xây dựng phương trình hồi qui thiết lập mối quan hệ thời gian cấp đông với nhiệt độ, kết thể bảng hình Bảng 6: Quan hệ thời gian cấp đông với nhiệt độ môi trường cấp đông Phương trình hồi quy =4688,26 + 128,166.te + 1,10857.t = 3862,23 e + 106,023.te + 0,914286te2+ 94,8686.te + =3414,09 0,822857.t = 3098,43e2+ 86,7429.te + 0,754286.t = 2852,31e2+ 80,1314.te + 0.697143.t Quan hệ nhiệteđộ khơng khí với thời gian cấp Vận tốc =5m/s =7,5m/s =10m/s =12,5m/s =15m/s đông các chế độ vận tốc khác 1600  =5m/s  =7,5m/s  =10m/s 1500  =12,5m/s Thời gian cấp đông, (s) 1400  =15m/s 1300 1200 1100 1000 900 800 700 600 -45 -44 -43 -42 -41 -40 -39 -38 -37 -36 -35 Nhiệt độ khơng khí, (0C) Hình 6: Quan hệ thời gian cấp đông với nhiệt độ khơng khí Dựa vào mới quan hệ trình bày bảng 5, ta thấy rằng giảm nhiệt độ khơng khí làm thời gian đơng kết giảm Độ giảm thời gian ứng với độ (te =10C) được ghi bảng Bảng Độ giảm thời gian cấp đông (s) ứng với te =-10C chế độ vận tốc khác Vận tốc -45 -42,5 -40 28 34 39 37,5 45 =5m/s 28 33 37 =7,5m/ 24 21 25 29 33 s=10m/ 19 23 26 30 s =12,5 17 21 24 28 m/s =15m/ s tác giả có sớ nhận xét sau: Qua kết tính tốn, 35 51 42 37 34 31 TẠP CHÍ KHOA HỌC CÔNG NGHỆ & THỰC PHẨM SỐ 09/2016 46 KHOA HỌC CƠNG NGHỆ - Với nhiệt độ khơng khí mơi trường cấp đơng khoảng te= (-45 -35)0C, giảm nhiệt độ khơng khí lạnh làm thời gian cấp đông giảm, nhiên mức độ giảm không đồng đều ở nhiệt độ khác nhau, cụ thể:  Trong khoảng nhiệt độ te=(-37,5-35) 0C, nhiệt độ môi trường cấp đông giảm te =10C, độ giảm thời gian trung bình 3,2%  Trong khoảng nhiệt độ te=(-40-37,5) 0C, nhiệt độ môi trường cấp đông giảm te=10C, độ giảm thời gian trung bình 3%  Trong khoảng nhiệt độ te=(-42,5-40) 0C, nhiệt độ môi trường cấp đơng giảm te=10C, độ giảm thời gian trung bình 2,85%  Trong khoảng nhiệt độ te=(-42,5-45) 0C, nhiệt độ môi trường cấp đông giảm te=10C, độ giảm thời gian trung bình 2,5% 2.3.3 Ảnh hưởng chiều dày sản phẩm đến thời gian cấp đông Trên sở thời gian cấp đông cá tra fillet được xác định ở chế độ khác cho kích thước khác Tương tự trên, tác giả tiến hành xây dựng phương trình hồi quy thời gian cấp đông với chiều dày sản phẩm ở chế độ khác nhau, kết thể bảng hình Bảng 7: Quan hệ thời gian cấp đơng với kích thước sản phẩm Thời gian cấp đơng, (s) CĐ (a) (b) (c) (d) (e) 3000 2500 Phương trình hồi quy =-90,0571 + 101,271. + = -71,6857 + 75,7905.+ 0,595238. =-49,6571 + 59,3714. + 0,47619. = -6,57143 + 44,5143. + 0,428571. =14,9429 + 34,7714. + 0,507937. 0.539683. (a) =5m/s, te=-350C (b) =7,5m/s, te=-37,50C (c) =10m/s, te=-400C (d) =12,5m/s, te=-42,50C (e) =15m/s, te=-450C (e) (d) (c) 2000 (a) (b) 150 100 50 012 14 16 18 20 22 24 Chiều dày cá tra fillet, (mm) Hình Quan hệ thời gian cấp đơng với kích thước sản phẩm ở chế độ cấp đơng khác Từ phương trình hồi quy, ta tìm được phụ thuộc độ giảm thời gian  vào độ tăng tốc độ  theo hàm trình bày bảng Bảng Độ tăng Thời gian cấp đông ,(S) ứng với =1mm chế độ khác CĐ Kích thước sản phẩm, (mm) TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 47 KHOA HỌC CÔNG NGHỆ (a) (b) (c) (d) (e) 12x70 15x85 18x100 21x115 24x125 116 119 123 126 130 87 90 93 96 99 70 72 75 77 80 57 60 63 66 69 48 51 54 57 61 Nhận xét  Theo kết trình bày hình bảng cho thấy rằng, đới với kích thước sản phẩm tăng ở chế độ nhiệt độ cao, vận tốc thấp thời gian cấp đông tăng đáng kể  Xét chế độ nhiệt độ vận tớc, kích thước sản phẩm tăng, thời gian cấp đông tăng với mức độ tăng gần theo quy ḷt tuyến tính (hình 8)  Ở nhiệt độ cao, vận tốc thấp, chiều dày sản phẩm tăng, thời gian cấp đông tăng nhanh ở nhiệt độ thấp, vận tớc cao (hình bảng 8)  Ở chế độ (=5m/s,te=-450C) thời gian cấp đông cho sản phẩm có kích thước lớn khơng phù hợp Thời gian cấp đông lớn, cụ thể đối với mẫu 24x125mm có thời gian cấp đơng =1975 s 2.3.4 Xây dựng phương trình hồi quy thiết lập mối quan hệ thời gian cấp đông với yếu tố ảnh hưởng Như phần đánh giá ở chỉ đơn đánh tác động riêng lẻ đến thời gian cấp đông Để đánh giá cách tổng quát, cần thiết phải xây dựng phương trình hồi quy đa yếu tớ thiết lập mối quan hệ thời gian cấp đông đến yếu tố ảnh hưởng Phần mềm Statghraphics Centurion XV.II được sử dụng để phân tích phương sai phân tích hồi quy kết quả, theo quy hoạch trực giao, tác giả xây dựng bảng số liệu được mã hóa yếu tớ ảnh hưởng đến thời gian cấp đông với giá trị tương ứng 1 (bảng 9) Bảng 9: Phạm vi biến đổi nhân tớ độc lập q trình cấp đơng Nhân tớ Vận tớc Nhiệt độ Chiều dày Nhân tớ mã hố X1 X2 X3 Đơn vị m/s C mm Mức nghiên cứu -1 -45 12 10 -40 18 15 -35 24 Mơ hình tốn biểu diễn phụ thuộc thời gian cấp đông vào nhân tớ được mã hóa phương trình đa thức bậc hai sau: n k 1 k k y  0   i Xi    ijXi X j  ii Xi2 i1 i 1 ji 1 (2) i1 Từ kết tính tốn theo phương án thiết kế, bằng chương trình Statghraphics Centurion XV.II phương trình hồi quy mơ tả q trình nghiên cứu dạng biến mã hóa sau: Phương trình hồi quy mơ tả q trình nghiên cứu dạng biến mã hóa: TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 48 KHOA HỌC CÔNG NGHỆ  = 2586,33 - 185,385.X1 + 98,7099.X + 206,536.X - 1,84448.X1 X (3) - 3,94533.X1 X + 2,6368X X + 5,72709X1 + 1,11451X + 0,509841X 2 Phương trình hồi quy mơ tả q trình nghiên cứu dạng biến thực:  = 2586,33 - 185,385. + 98,7099.t e + 206,536. - 1,84448..t e - 3,94533.. + 2,6368.t e  (4) + 5,727092 + 1,11451t e + 0,509841  2 Độ tương quan thời gian cấp đơng thu được từ phương trình hồi quy sớ liệu tính tốn R2=99,7375% Như vậy có thể kết ḷn rằng: mơ hình hồi qui mơ tả với kết tính tốn Bảng 10: Kết quả xác định thời gian cấp đông cá tra fillet Kích thước (dàyrộng)(mm) 12x70 15x85 18x100 21x115 21x115 Nhiệt độ te (0C) -35 -37,5 -40 -42,5 -45 -35 -37,5 -40 -42,5 -45 -35 -37,5 -40 -42,5 -45 -35 -37,5 -40 -42,5 -45 -35 -37,5 -40 -42,5 -45 1212 1119 1039 969 909 1561 1440 1335 1245 1165 1925 1773 1643 1530 1430 2302 2118 1960 1823 1703 2682 2465 2279 2117 1975 Vận tớc khơng khí  (m/s) 7,5 10 12,5 Thời gian cấp đông (s) 985 851 761 907 783 699 842 725 647 784 675 601 734 631 561 1272 1102 987 1171 1013 905 1084 936 836 1009 869 775 942 811 722 1573 1366 1226 1447 1253 1122 1336 1159 1035 1242 1072 959 1158 999 892 1886 1641 1477 1731 1504 1352 1598 1386 1245 1483 1285 1153 1382 1197 1073 2204 1924 1737 2021 1761 1589 1864 1622 1463 1728 1504 1354 1610 1399 1260 15 695 638 589 547 510 902 827 763 706 658 1123 1028 946 876 815 1357 1241 1143 1058 984 1601 1464 1347 1247 1160 Kết luận TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 49 KHOA HỌC CÔNG NGHỆ  Ảnh hưởng vận tớc khơng khí Với tớc độ khoảng (5 15) m/s, tăng tốc độ khơng khí lạnh làm thời gian cấp đơng giảm đáng kể Tuy nhiên độ giảm thời gian cấp đông không đồng đều chế độ nhiệt độ với chế độ vận tốc khác nhau, ở vận tốc thấp tăng vận tốc thời gian cấp đông giảm lớn ở vận tốc cao, =5m/s tăng =1m/s độ giảm thời gian cấp đông 7,5% (   103s ), =15m/s tăng =1m/s độ giảm thời gian cấp đông 2%(   10s )  Ảnh hưởng nhiệt độ khơng khí Với nhiệt độ khơng khí mơi trường cấp đơng khoảng te= (-45 -35)0C, giảm nhiệt độ khơng khí lạnh làm thời gian cấp đông giảm, nhiên mức độ giảm không đồng đều ở nhiệt độ khác nhau, ở nhiệt độ cao, giảm nhiệt độ thời gian cấp đông giảm lớn ở nhiệt độ thấp, te=-350C giảm te=10C thời gian cấp đông giảm trung bình 3,2%(   42s ), te=-450C giảm te=10C thời gian cấp đơng giảm trung bình 2,5%(   22s )  Ảnh hưởng kích thước sản phẩm Khi kích thước cá tra fillet tăng, thời gian cấp đơng tăng theo quy ḷt tuyến tính với độ tăng trung bình khoảng   82s chiều dày cá tra fillet tăng =1mm đối với loại cá tra fillet có chiều dày >24mm thời gian cấp đơng lớn, không phù hợp cho băng chuyền IQF Trên sở thời gian cấp đông ở chế độ khác cho kích thước khác nhau, tác giả xây dựng được phương trình hồi qui thiết lập mối quan hệ thời gian cấp đông với yếu tớ ảnh hưởng bao gồm: vận tớc khơng khí, nhiệt độ khơng khí chiều dày sản phẩm TÀI LIỆU THAM KHẢO [1] AC Cleland, RL Earle (1977), A comparison of analytical and numerical methods of predicting the Freezing times of foods Journal of Food Science (42), p1390–1395 [2] AC Cleland, R L Earle (1984), Freezing time prediction for different final product temperatures Journal of Food Science (49), p1230–1232 [3] AC Cleland (1990), Food Refrigeration Processes – Analysis, Design and Simulation New York: Elsevier Applied Science, p 79–152 [4] Bùi Minh Trí (2003), Mơ hình tốn kinh tế, NXB Khoa học Kỹ thuật, HN [5] Đặng Quốc Phú-Trần Thế Sơn-Trần Văn Phú (2004), Truyền nhiệt, Nhà xuất Giáo Dục [6] Đỗ Hữu Hồng, Đặng Q́c Phú (02/2011) Công nghệ chế biến cá da trơn khu vực đồng bằng sơng Cửu Long, tạp chí KH&CN Trường Đại học Cơng nghiệp Tp Hồ Chí Minh, tr.54  tr.60 [7] Hoang Huu DO, Phu Quoc DANG, Duy Khanh HOANG (10/2013), Determining the thermal properties of the catfish in freezing and chilling using a dual needle heat-pulse sensor, the 3rd International Conference on Sustainable Energy – “RISE towards a green future”, Ho Chi Minh University of Technology, Vietnam, p.153  p.159 [8] Hoang Huu DO, Phu Quoc DANG, Dung Viet NGUYEN (10/2013), Determining the time to freeze catfish by applying the finite element analysis, the 3rd International TẠP CHÍ KHOA HỌC CÔNG NGHỆ & THỰC PHẨM SỐ 09/2016 50 KHOA HỌC CÔNG NGHỆ Conference on Sustainable Energy – “RISE towards a green future”, Ho Chi Minh University of Technology, Vietnam, p.160  p.165 [9] Hồng Đình Tín (1996), Trùn nhiệt thiết bị trao đổi nhiệt, Nhà xuất Khoa học Kỹ thuật [10] Industrial refrigeration handbook, (2004) the Mcgraw-hill Companies p (567-599) [11] Nguyễn Bốn, Trần Quốc Minh (2004), Tính Tốn Thời Gian Đơng Lạnh Thực Phẩm, tạp chí KH&CN Nhiệt, tr.11÷tr.13 [12] N.Scheerlinck, P.Verboven, K.A Fikiin, J.DeBaerdemaeker, B.M Nicolai, (2001), Finite element computation of unsteady phase change heat transfer during freezing or thawing of food using a combined enthalpy and Kirchhop transform Method, Trasactions of the ASAE, American of Agricultural Engineers 44(2), p429-438 [13] QT Pham, (1984) Extension to Plank’s equation for predicting freezing times of foodstuffs of simple shapes International Journal of Refrigeration (7), p.377–383 [14] Roland W Lewis Fundamentals of the Finite Element Method for Heat and Fluid Flow, pp 154-170 [15] Số liệu thu tập từ nhà máy chế biến thủy hải sản (2010-2012) (Cửu long, Nam việt, Vĩnh Hồn, Afiex, Cơng ty Cataco, Công ty Cafatex, Afiex….) [16] Cá tra Việt Nam – Câu chuyện giành 99,9% thị phần (07 /8 /2009) http:// tintuc.xalo.Vn /041800223749 /ca_tra _viet_nam_cau_chuyen_gianh.html TẠP CHÍ KHOA HỌC CƠNG NGHỆ & THỰC PHẨM SỐ 09/2016 51 ... giải tốn dẫn nhiệt có biến đổi pha để xác định thời gian cấp đông ở chế độ Trên sở thời gian cấp đông được xác định, tác giả tiến hành đánh giá chi tiết yếu tố ảnh hưởng đến thời gian cấp đông. .. chế độ cấp đông khác Kết xác định thời gian cấp đơng được trình bày bảng 10 2.3 Đánh giá các yếu tố ảnh hưởng đến thời gian cấp đơng Đới với q trình cấp đơng cá tra fillet, thông số ảnh hưởng... giảm thời gian trung bình 2% vận tớc khơng khí tăng =1m/s 2.3.2 Ảnh hưởng nhiệt độ khơng khí đến thời gian cấp đông Tương tự phần đánh giá ảnh hưởng vận tớc khơng khí đến thời gian cấp đông