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

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 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 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu kết nghiên cứu nêu luận án trung thực chƣa đƣợc khác công bố cơng trình Tác giả Vũ Huy Khuê LỜI CẢM ƠN Tôi xin trân trọng cám ơn Trƣờng Đại học Bách khoa Hà nội, Viện đào tạo sau đại học Trƣờng ĐHBK Hà nội, Viện Khoa hoc Công nghệ Nhiệt –Lạnh, Bộ môn Kỹ thuật lạnh điều hồ khơng khí, Giáo sƣ, Tiến sỹ, đồng nghiệp Viện Khoa học Công nghệ Nhiệt Lạnh tạo điều kiện thuận lợi, giúp đỡ, đóng góp nhiều ý kiến q báu giúp tơi hồn thành luận án Tơi xin bày tỏ lịng biết ơn chân thành đến TS Nguyễn Xuân Tiên GS TSKH Trần Văn Phú hƣớng dẫn giúp đỡ tơi q trình thực luận án Tơi xin chân thành cám ơn TS Nguyễn Việt Dũng, trƣởng mơn Kỹ thuật lạnh Điều hồ khơng khí, Viện trƣởng Viện KH&CN Nhiệt Lạnh, tạo điều kiện, hƣớng dẫn giúp đỡ thời gian thực luận án Bộ môn Tôi xin trân trọng cảm ơn Bộ môn Kỹ thuật Nhiệt, Khoa Điện Trƣờng Đại Học Công nghiệp Hà nội tạo điều kiện giúp đỡ tơi q trình làm thí nghiệm phịng thí nghiệm Bộ mơn Tơi xin trân trọng cảm ơn anh em, bạn bè đồng nghiệp ngồi nƣớc giúp đỡ tơi q trình học tập làm thí nghiệm, tiếp cận với hàng trăm báo quốc tế liên quan đến vấn đề nghiên cứu luận án Cuối xin cảm ơn Bố, Mẹ, Vợ, Con, ngƣời thân bạn bè ủng hộ, động viên suốt chặng đƣờng nghiên cứu thực luận án Tác giả Vũ Huy Khuê ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii PHỤ LỤC v DANH MỤC BẢNG vi DANH MỤC HÌNH VẼ, ĐỒ THỊ vii DANH MỤC CÁC KÝ HIỆU SỬ DỤNG TRONG LUẬN ÁN ix LỜI MỞ ĐẦU xii Chƣơng NGHIÊN CỨU TỔNG QUAN 1.1 Tổng quan công nghệ lạnh đông thực phẩm 1.1.1 Vị trí, vai trị cơng nghệ lạnh đơng 1.1.2 Thực trạng, tiềm công nghệ lạnh thủy sản Việt Nam 1.1.3 Tiêu hao lƣợng, hao hụt khối lƣợng thời gian q trình cấp đơng 1.1.4 Cá thu, loại hải sản có giá trị dinh dƣỡng cao 1.2 Tổng quan mơ hình truyền nhiệt truyền chất q trình cấp đơng thực phẩm 1.2.1 Mơ hình toán truyền nhiệt truyền chất liên hợp trình cấp đơng thực phẩm 1.2.1.1 Mơ hình tốn truyền nhiệt truyền chất tổng quát 1.2.1.2 Vấn đề truyền nhiệt truyền chất trình cấp đơng 1.2.1.3 Mơ hình truyền chất xảy q trình cấp đông 1.2.2 Các phƣơng pháp giải tốn truyền nhiệt khơng ổn định q trình cấp đơng 1.2.2.1 Mơ hình tốn dẫn nhiệt khơng ổn định q trình cấp đơng 1.2.2.2 Phƣơng pháp giải tích gần cơng thức bán thực nghiệm kinh điển 10 1.2.2.3 Ứng dụng phƣơng pháp số giải toán truyền nhiệt cấp đông thực phẩm 14 1.3 Tổng quan mơ hình xác định tính chất nhiệt vật lý 17 1.3.1 Thành phần băng thực phẩm kết đông 17 1.3.2 Nhiệt dung riêng 18 1.3.3 Entanpi 20 1.3.4 Khối lƣợng riêng 21 1.3.5 Hệ số dẫn nhiệt 22 1.4 Kết luận 25 1.5 Xác định đề tài luận án 26 1.5.1 Tên đề tài luận án 26 1.5.2 Mục tiêu nghiên cứu 26 1.5.3 Đối tƣợng nghiên cứu 27 1.5.4 Phƣơng pháp nghiên cứu 27 1.5.5 Nội dung nghiên cứu 27 Chƣơng 28 MƠ HÌNH XÁC ĐỊNH TÍNH CHẤT NHIỆT VẬT LÝ CỦA CÁ THU 28 2.1 Phát triển mơ hình xác định tính chất nhiệt vật lý cá thu 28 2.1.1 Thành phần băng thực phẩm kết đông 28 2.1.2 T lệ đóng băng nƣớc trình cấp đơng thực phẩm 30 2.1.3 Nhiệt dung riêng hiệu dụng 31 2.1.4 Entanpi 34 2.1.5 Hệ số dẫn nhiệt hiệu dụng 34 iii 2.1.6 Khối lƣợng riêng 37 2.2 Xây dựng hoàn thiện phần mềm xác định tính chất nhiệt vật lý cá thu theo nhiệt độ 38 2.3 Đánh giá độ tin cậy mô hình tính chất nhiệt vật lý 39 2.3.1 Đối tƣợng nghiên cứu kiểm chứng 39 2.3.2 So sánh kết tính tốn kết thực nghiệm t lệ nƣớc đóng băng 40 2.3.3 So sánh kết tính tốn kết thực nghiệm nhiệt dung riêng hiệu dụng 41 2.3.4 So sánh kết tính tốn kết thực nghiệm entanpi 43 2.3.5 So sánh kết tính tốn kết thực nghiệm hệ số dẫn nhiệt hiệu 45 Chƣơng 3: 48 NGHIÊN CỨU ẢNH HƢỞNG CỦA TRUYỀN CHẤT VÀ MÔ PHỎNG QUÁ TRÌNH CẤP ĐƠNG CÁ THU 48 3.1 Đánh giá ảnh hƣởng trình trao đổi chất đến q trình cấp đơng cá thu 48 3.2 Lựa chọn mơ hình mơ q trình cấp đơng cá thu 48 3.2.1 Thí nghiệm kiểm chứng mơ hình 48 3.2.2 Những giả thiết sử dụng lựa chọn mô hình tốn học tốn dẫn nhiệt thực phẩm cấp đông 50 3.2.3 Bài toán dẫn nhiệt thực phẩm cấp đơng theo mơ hình Chumak & Onhishenko 50 3.2.4 Lựa chọn mơ hình tốn học theo mơ hình cấp đơng Chumak & Onhishenko 51 3.3 Giải tốn dẫn nhiệt q trình cấp đơng cá thu 54 3.3.1 Phƣơng pháp thể tích hữu hạn giải hệ phƣơng trình dẫn nhiệt q trình cấp đơng cá thu 54 3.3.2 Thuật toán giải toán dẫn nhiệt bên q trình cấp đơng cá thu theo phƣơng pháp thể tích hữu hạn xác định trƣờng nhiệt độ 57 3.4 Xây dựng hồn thiện phầm mềm mơ q trình cấp đông cá thu 59 3.5 Xác định nhiệt ẩn trình trao đổi chất 61 3.5.1 Khái niệm nhiệt độ trung bình entanpi tốn xác định phụ tải lạnh cấp đông thực phẩm 61 3.5.2 Xác định nhiệt ẩn trình trao đổi chất 62 3.6 Đánh giá độ tin cậy mơ hình mơ q trình cấp đơng cá Thu 63 3.6.1 Xây dựng hệ thống thiết bị thí nghiệm 63 3.6.1.1 Mô tả tủ cấp đơng thí nghiệm 63 3.6.1.2 Thiết bị đo 64 3.6.1.3 Xác định thông số đặc trƣng mơ hình thí nghiệm 65 3.6.2 Nghiên cứu thực nghiệm xác định độ hụt khối lƣợng cấp đông 68 Kết luận 70 3.6.3 Kết nghiên cứu thực nghiệm trƣờng nhiệt độ thời gian cấp đông cá thu 71 3.7 Kết luận chƣơng 77 Chƣơng 79 ĐÁNH GIÁ ẢNH HƢỞNG CỦA CÁC YẾU TỐ ĐẾN Q TRÌNH CẤP ĐƠNG VÀ CÁC GIẢI PHÁP TIẾT KIỆM NĂNG LƢỢNG TRONG QUÁ TRÌNH CẤP ĐÔNG CÁ THU 79 4.1 Đánh giá ảnh hƣởng yếu tố mơi trƣờng tới q trình cấp đơng cá Thu 79 4.1.1 Đánh giá ảnh hƣởng vận tốc khơng khí tới q trình cấp đơng cá thu 79 4.1.2 Đánh giá ảnh hƣởng nhiệt độ mơi trƣờng làm lạnh tới q trình cấp đơng thực phẩm 82 4.1.3 Đánh giá ảnh hƣởng chiều dày sản phẩm tới trình cấp đơng thực phẩm 86 4.2 Nghiên cứu xác định vùng cấp đông hợp lý cho cá thu 89 iv 4.2.1 Xây dựng hàm hồi quy 89 4.2.2 Tìm vùng cấp đông hợp lý cho cá Thu 91 4.3 Đề xuất giải pháp tiết kiệm lƣợng q trình cấp đơng cá thu 92 4.3.1 Xác định chế độ cấp đông tối ƣu theo phƣơng pháp thực nghiệm 92 4.3.2 Giải pháp tiết kiệm lƣợng thiết bị cấp đông dạng IQF nhà máy 93 4.3.2.1 Ứng dụng kết nghiên cứu cho hệ thống cấp đông IQF 94 4.3.2.2 Khả tiết kiệm lƣợng thay đổi tốc độ quạt lƣu chuyển gió 95 4.3.2.3 Khả tiết kiệm lƣợng thay đổi nhiệt độ buồng cấp đông chênh lệch nhiệt độ sôi nhiệt độ buồng 96 4.4 Kết luận chƣơng 98 Chƣơng 100 KẾT LUẬN VÀ KIẾN NGHỊ 100 5.1 Kết luận 100 5.2 Những đóng góp luận án 101 5.3 Kiến nghị đề xuất hƣớng nghiên cứu 101 TÀI LIỆU THAM KHẢO 102 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 108 PHỤ LỤC v DANH MỤC BẢNG Bảng 1.1: GDP thủy sản kinh tế Quốc dân giai đoạn 2001-2011 Bảng 1.2:Các cơng thức hiệu chỉnh cơng thức tính thời gian cấp đông Plank 11 Bảng 1.3: Hằng số hình dáng cho tính tốn hệ số E 13 Bảng 1.4 Khối lƣợng riêng phụ thuộc nhiệt độ thành phần thực phẩm 22 Bảng 2.1: Quan hệ độ ẩm điểm kết đông số loại thực phẩm 30 Bảng 2.2: Nhiệt dung riêng phụ thuộc vào nhiệt độ thành phần thực phẩm 32 Bảng 2.3 Hệ số dẫn nhiệt phụ thuộc vào nhiệt độ thành phần thực phẩm 35 Bảng 2.4 Khối lƣợng riêng phụ thuộc nhiệt độ thành phần thực phẩm[30] 37 Bảng 2.5: Thành phần khối lƣợng nhiệt độ bắt đầu kết đông số loại cá 40 Bảng 3.1 Kết khảo sát tốc độ quạt buồng cấp đông 65 Bảng 3.2: Các chế độ thí nghiệm với cá thu tủ cấp đông KTE-5000 LT 68 Bảng 3.3: Kết đo thời gian cấp đông độ hụt khối lƣợng cấp đông cá Thu 70 Bảng 3.4: Ảnh hƣởng độ hụt khối tới tải nhiệt cấp đông thực phẩm 70 Bảng 3.5: Thời gian cấp đông thực nghiệm tính tốn cá Thu 76 Bảng 4.1: Kết xác định thời gian cấp đơng cá thu cắt lát có chiều dày 3mm 79 Bảng 4.2: Quan hệ vận tốc không khí đến thời gian cấp đơng 81 Bảng 4.3: Mối quan hệ  = f() 81 Bảng 4.4: Mức độ ảnh hƣởng chế độ vận tốc khác 82 Bảng 4.5: Quan hệ nhiệt độ không khí đến thời gian cấp đơng 84 Bảng 4.6: Quan hệ  =f(T) 84 Bảng 4.7: Mức độ ảnh hƣởng chế độ nhiệt độ khơng khí khác 85 Bảng 4.8: Phƣơng trình hồi quy thiết lập mối quan hệ thời gian cấp đông chiều dày sản phẩm chế độ nhiệt độ khác 86 Bảng 4.9: Mức độ ảnh hƣởng chiều dày sản phẩm đến thời gian cấp đông 87 Bảng 4.10: Phƣơng trình hồi quy thiết lập mối quan hệ thời gian cấp đông chiều dày sản phẩm chế độ vận tốc khác 88 Bảng 4.11: Mức độ ảnh hƣởng chiều dày sản phẩm đến thời gian cấp đông 88 Bảng 4.12: Phạm vi biến đổi nhân tố độc lập q trình cấp đơng 89 Bảng 4.13: Số liệu xây dựng phƣơng trình hồi quy (4.3) 90 Bảng 4.14:Kết tối ƣu thời gian cấp đông 91 Bảng 4.15: Khoảng giá trị t, ω, chiều dày hợp lý cấp đông cá Thu 92 Bảng 4.16: Công suất quạt thay đổi theo vận tốc gió 96 Bảng 4.17: Công suất máy nén thay đổi theo nhiệt độ sôi 97 Bảng 4.18: Mức độ tiết kiệm lƣợng thay đổi nhiệt độ sôi môi chất 97 vi DANH MỤC HÌNH VẼ, ĐỒ THỊ Hình 1.1: Kim ngạch xuất thủy sản Việt Nam tăng trƣởng qua giai đoạn 2003 – 2013 Hình 1.2: Mơ hình cấp đơng theo phƣơng pháp CAS Hình 1.3: a, b Hình ảnh cá thu thƣơng phẩm ;c Sản lƣợng đánh bắt cá thu đến 2010 Hình 1.4: Ảnh hƣởng thời gian cấp đông tới chiều dày thoát ẩm tốc độ dehydration Campanone Hình 1.5: Ảnh hƣởng hoạt hóa nƣớc đến thời gian cấp đơng Hình 1.6: Mơ tả nguồn gốc phƣơng trình Plank 11 Hình 1.7: Sơ đồ chia lƣới điểm nút trục không gian 0x giá trị hàm Fi tƣơng ứng với nút i lƣới 15 Hình 2.1: Lƣu đồ thuật toán xác định nhiệt dung riêng hiệu dụng 33 Hình 2.2: Lƣu đồ thuật tốn xác định entanpi thực phẩm 34 Hình 2.3: Lƣu đồ thuật tốn tính hệ số dẫn nhiệt hiệu dụng 37 Hình 2.4: Giao diện khởi động phần mềm xác định tính chất nhiệt vật lý thời gian cấp đông 38 Hình 2.5: Giao diện làm việc phần mềm 39 Hình 2.6: Đồ thị biến thiên tính chất nhiệt vật lý theo nhiệt độ File liệu dạng excel phần mềm 39 Hình 2.7: Sự biến đổi t lệ nƣớc đóng băng theo nhiệt độ với thịt bị 40 Hình 2.8: Sự biến đổi t lệ nƣớc đóng băng theo nhiệt độ với cá Tuyết 40 Hình 2.9: Sự biến đổi t lệ nƣớc đóng băng theo nhiệt độ với cá Vƣợc 41 Hình 2.10: Sự biến đổi nhiệt dung riêng theo nhiệt độ cá Ngừ 41 Hình 2.11: Sự biến đổi nhiệt dung riêng theo nhiệt độ cá Hồi 42 Hình 2.12: Sự biến đổi nhiệt dung riêng theo nhiệt độ cá thu 42 Hình 2.13: Sự biến đổi entanpi theo nhiệt độ cá Ngừ 43 Hình 2.14: Sự biến đổi entanpi theo nhiệt độ cá Hồi 43 Hình 2.15: Sự biến đổi entanpi theo nhiệt độ cá thu 43 Hình 2.16: Sự biến đổi entanpi theo nhiệt độ cá Tuyết 44 Hình 2.17: Sự biến đổi entanpi theo nhiệt độ cá Vƣợc 44 Hình 2.18: Sự biến đổi hệ số dẫn nhiệt theo nhiệt độ cá Ngừ 45 Hình 2.19: Sự biến đổi hệ số dẫn nhiệt theo nhiệt độ cá Hồi 45 Hình 2.20: Sự biến đổi hệ số dẫn nhiệt theo nhiệt độ cá thu 45 Hình 3.1a: Mơ hình mẫu cá thí nghiệm 49 Hình 3.1b: Mơ hình lát cá thu cấp đông băng chuyền lƣới 49 Hinh 3.1c: Biểu diễn phân tố theo tọa độ thời gian 54 Hình 3.2: Lƣu đồ thuật tốn xác định trƣờng nhiệt độ thực phẩm đƣợc cấp đơng 58 Hình 3.3: Lƣu đồ thuật tốn xác định trƣờng nhiệt độ bƣớc thời gian 59 Hình 3.4: Giao diện phần mềm mơ trình dẫn nhiệt thực phẩm đƣợc cấp đông 60 Hình 3.5: Dạng đồ thị nhiệt độ tâm bề mặt điển hình 61 Hình 3.6: Dữ liệu đƣợc xuất dƣới dạng file excel 61 Hình 3.7: Tủ thí nghiệm cấp đơng gió KTE-5000LT 64 Hình 3.8: Chu trình lạnh thiết bị cấp đơng 64 Hình 3.9: Khối tín hiệu vào NI 9211A Khối truyền thông cDAQ-9174 Dataloger nhiệt độ kênh National Instruments Corporation 64 Hình 3.10: Bộ đo nhiệt độ chuyển mạch tay đầu đo dải nhiệt độ đo -1000C đến 1000C 64 vii Hình 3.11: Dataloger đo nhiệt độ di động Datatrace MPIII, MPRF dải nhiệt độ đo -800C đến 140oC 64 Hình 3.12: Thiết bị đo tốc độ gió TSI 8347 65 Hình 3.13: Cân kỹ thuật UEC – 7500; 65 Hình 3.14 a) Kích thƣớc buồng cấp đơng b) Vị trí điểm đo xác định trƣờng tốc độ gió, c) Vị trí gắn điểm đo miếng cá Thu 66 Hình 3.15: Bố trí sản phẩm tủ cấp đơng (a), gắn cặp nhiệt lên miếng cá (b) 67 Hình 3.16: Quy trình thí nghiệm q trình cấp đơng cá Thu 67 Hình 3.17: Bố trí thí nghiệm xác định độ hao hụt khối lƣợng cá thu cấp đông 68 Hình 3.18: Mật độ ẩm mẫu cá thu cấp đông theo thời gian 69 Hình 3.19: Sự biến đổi nhiệt độ khối thực phẩm theo thời gian ( mẫu 5) 69 Hình 3.20: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 01 71 Hình 3.21: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 02 72 Hình 3.22: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 03 72 Hình 3.23: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 04 73 Hình 3.24: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 05 73 Hình 3.25: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 06 74 Hình 3.26: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 07 74 Hình 3.27: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 08 75 Hình 3.28: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 09 75 Hình 3.29: Sự sai khác trƣờng nhiệt độ tính tốn thực nghiệm với cá thu chế độ thí nghiệm 10 76 Hình 4.1: Ảnh hƣởng tốc độ gió đến thời gian cấp đông cá thu cắt lát 30mm 80 Hình 4.2: Quan hệ vận tốc khơng khí đến thời gian cấp đơng 81 Hình 4.3 : Ảnh hƣởng nhiệt độ đến thời gian cấp đông cá thu cắt lát 30mm 83 Hình 4.4: Quan hệ thời gian cấp đơng với nhiệt độ khơng khí 84 Hình 4.5: Quan hệ kích thƣớc sản phẩm thời gian cấp đông 87 Hình 4.6: Quan hệ chiều day sản phẩm thời gian cấp đông 88 Hình 4.7: Mối quan hệ tƣơng quan số liệu tính tốn số liệu từ mơ hình tốn (4.3) 90 Hình 4.8: Lƣợc đồ phƣơng pháp tối ƣu hóa thực nghiệm 93 viii DANH MỤC CÁC KÝ HIỆU SỬ DỤNG TRONG LUẬN ÁN Ký hiệu C Cp Cf Cu Ca C B Bi E F S f h hf L Lo Ms ΔC n T t Tin To Tr Tref T q R Rg rf V v r x m Ký tự Hy lạp α Đại lƣợng Nhiệt dung riêng Nhiệt dung riêng đẳng áp Nhiệt dung riêng thực phẩm đóng băng hồn tồn Nhiệt dung riêng thực phẩm điểm đóng băng Nhiệt dung riêng biểu kiến Ma trận nhiệt dung T N Tiêu chuẩn Biot (αR/λ) Hệ số hình dáng thời gian lạnh đơng Diện tích bề mặt Diện tích tiết diện Ma trận nguồn nhiệt Entanpi Enthanlpy khối thực phẩm nhiệt độ đông đặc, Nhiệt ẩn đóng băng Nhiệt ẩn đóng băng nƣớc, Lo = 333600 J/kg Phân tử lƣợng chất khô hòa tan Chênh lệch nhiệt dung riêng nƣớc băng Vectơ đơn vị theo phƣơng pháp tuyến với bề mặt Nhiệt độ Nhiệt độ Nhiệt độ ban đầu Nhiệt độ kết đông nƣớc, To=273,15 K Nhiệt độ tham chiếu (entanpi = 0), Tr = 233,15K (-40oC) Nhiệt độ tham chiếu phƣơng trình Cleland & Earle, Tref = 263,15K Vectơ nút nhiệt độ Nguồn nhiệt phân tố Kích thƣớc xác định Hằng số khí lý tƣởng, Rg = 8314 J/kgK Bề dày lớp băng Thể tích Vận tốc Tọa độ vị trí Thành phần khối lƣợng Khối lƣợng thực phẩm Đơn vị J.kg−1K−1 J.kg−1K−1 J.kg−1K−1 J.kg−1K−1 J.kg−1K−1 m2 m2 J.kg-1 J.kg-1 J.kg−1 J.kg-1 kg.kmol-1 K C K K K K J.m-3 m M m3 ms-1 M % kg kg Wm-2K-1 Hệ số toả nhiệt tổng ix function x_g_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function x_p_Callback(hObject, eventdata, handles) function x_p_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function T_cr_Callback(hObject, eventdata, handles) function T_cr_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function pushbutton1_Callback(hObject, eventdata, handles) Food_data=getfooddata(hObject, eventdata, handles); T=str2double(get(handles.T_property,'String'))+273.15; Trunggian=WT_Tchigeov(T,Food_data); x_ice=Trunggian(1); x_w=Trunggian(2); set(handles.x_ice,'string',x_ice) set(handles.x_w,'string',x_w) ro=klr(T,Food_data); set(handles.rho,'string',ro) ce=nhietdungrieng(T,Food_data); set(handles.ce,'string',ce) I=enthalpy(T,Food_data) set(handles.enthalpy,'string',I); hesodannhiet=hesodn1(T,Food_data); set(handles.hsdn,'string',hesodannhiet) function edit8_Callback(hObject, eventdata, handles) function edit8_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function T_property_Callback(hObject, eventdata, handles) function T_property_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end 154 function x_w_Callback(hObject, eventdata, handles) function x_w_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function rho_Callback(hObject, eventdata, handles) function rho_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function x_ice_Callback(hObject, eventdata, handles) function x_ice_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function hsdn_Callback(hObject, eventdata, handles) function hsdn_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function enthalpy_Callback(hObject, eventdata, handles) function enthalpy_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function ce_Callback(hObject, eventdata, handles) function ce_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function Untitled_1_Callback(hObject, eventdata, handles) function Untitled_2_Callback(hObject, eventdata, handles) function ten_tinh_chat_Callback(hObject, eventdata, handles) function ten_tinh_chat_CreateFcn(hObject, eventdata, handles) 155 if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function pushbutton2_Callback(hObject, eventdata, handles) Food_data=getfooddata(hObject, eventdata, handles); str=get(handles.ten_tinh_chat,'string'); val=get(handles.ten_tinh_chat,'value'); switch val case for i=1:81 t(i)=-41+i; T(i)=t(i)+273; %plot(handles.t,handles.L1,'k-'); %xlabel('t, [oC]'); %ylabel('Heat conductivity,landa[W/mK]'); end; for i=1:81 tg=WT_Tchigeov(T(i),Food_data); x_ice(i)=tg(1)*100; end plot(handles.tinhchat,t,x_ice,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'x_i_c_e, %'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') case for i=1:81 t(i)=-41+i; T(i)=t(i)+273.15; end; for i=1:81 tg=WT_Tchigeov(T(i),Food_data); x_w(i)=tg(2)*100; end plot(handles.tinhchat,t,x_w,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'x_w, %'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') case for i=1:81 t(i)=-41+i; T(i)=t(i)+273.15; end; for i=1:81 tg=klr(T(i),Food_data); x_w(i)=tg; end plot(handles.tinhchat,t,x_w,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'Khoi luong rieng, [kg/m3]'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') case for i=1:81 t(i)=-41+i; T(i)=t(i)+273.15; end; for i=1:81 156 tg=hesodn1(T(i),Food_data); x_w(i)=tg; end plot(handles.tinhchat,t,x_w,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'He so dan nhiet, [W/mK]'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') case for i=1:161 t(i)=-40.5+i*0.5; T(i)=t(i)+273.15; end; for i=1:161 tg=enthalpy(T(i),Food_data); x_w(i)=tg; end plot(handles.tinhchat,t,x_w,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'Enthalpy, [W/mK]'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') case for i=1:81 t(i)=-41+i; T(i)=t(i)+273.15; end; for i=1:81 tg=nhietdungrieng(T(i),Food_data); x_w(i)=tg; end plot(handles.tinhchat,t,x_w,'LineWidth',2) xlabel(handles.tinhchat,'T, [oC]'); ylabel(handles.tinhchat,'Nhiet dung rieng hieu dung, [kJ/kgK]'); set(handles.tinhchat,'YGrid','on') set(handles.tinhchat,'XGrid','on') end function kichthuoc_Callback(hObject, eventdata, handles) function kichthuoc_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function dientich_Callback(hObject, eventdata, handles) function dientich_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function nhietdobandau_Callback(hObject, eventdata, handles) function nhietdobandau_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function tocdogio_Callback(hObject, eventdata, handles) 157 function tocdogio_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function nhietdogio_Callback(hObject, eventdata, handles) function nhietdogio_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit22_Callback(hObject, eventdata, handles) function edit22_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function thetich_Callback(hObject, eventdata, handles) function thetich_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function pushbutton3_Callback(hObject, eventdata, handles) Food_data=getfooddata(hObject, eventdata, handles); for i=1:81 t(i)=-41+i; T(i)=t(i)+273.15; end; for i=1:81 tg=WT_Tchigeov(T(i),Food_data); data(i,1)=T(i); data(i,2)=t(i) data(i,3)=tg(1)*100; data(i,4)=tg(2)*100; data(i,5)=klr(T(i),Food_data); data(i,6)=hesodn1(T(i),Food_data); data(i,7)=enthalpy(T(i),Food_data); data(i,8)=nhietdungrieng(T(i),Food_data); end; tt=clock tg2=''; tg1=tt(1) tg2=num2str(tg1) tg1=tt(2) tg2=strcat(tg2,num2str(tg1)) tg1=tt(3) tg2=strcat(tg2,num2str(tg1)) tg1=tt(4) tg2=strcat(tg2,num2str(tg1)) tg1=tt(5) tg2=strcat(tg2,num2str(tg1),'Thermal_properties.xls') 158 a={'T, [K]','t, [oC]','x_ice, %', 'x_w, %','Khoi luong rieng, [kg/m3]','He so dan nhiet, [W/mK]','Enthalpy, [kJ/kg]','Nhiet dung rieng, [kJ/kgK]'}; xlswrite(tg2,a,'Thermal properties', 'A1'); xlswrite(tg2,data,'Thermal properties', 'A2'); 'completed' function pushbutton4_Callback(hObject, eventdata, handles) % hObject handle to pushbutton4 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) n=40; %so khoang chia kich thuoc xac dinh to=10; %A=truongnhietdo_crank_nicolson(hObject, eventdata, handles,n,to); A=truongnhietdogama(hObject, eventdata, handles,n,to); tg=size(A); i_max=tg(1); m=1; for i=1:i_max thoigian(i)=(i-1)*to/60; %chuyen sang don vi phut nhietdobemat(i)=A(i,n+1)-273.15; nhietdotam(i)=A(i,1)-273.15; data=[thoigian;nhietdobemat;nhietdotam]; end plot(handles.temp,thoigian,nhietdotam,thoigian,nhietdobemat,'LineWidth',1 5) xlabel(handles.temp,'Thoi gian, [phut]'); ylabel(handles.temp,'T, [oC]'); set(handles.temp,'YGrid','on') set(handles.temp,'XGrid','on') T_yeucau=str2double(get(handles.nhietdoyeucau,'String'))+273.15; for i=1:i_max if (A(i,1)