Nghiên cứu thực nghiệm quá trình làm mát phụ cho hệ thống điều hòa không khí co2 bằng nước ngưng

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHĨA LUẬN TỐT NGHIỆP NGÀNH CÔNG NGHỆ KỸ THUẬT NHIỆT NGHIÊN CỨU THỰC NGHIỆM QUÁ TRÌNH LÀM MÁT PHỤ CHO HỆ THỐNG ĐIỀU HỊA KHƠNG KHÍ CO2 BẰNG NƯỚC NGƯNG GVHD: PGS TS ĐẶNG THÀNH TRUNG SVTH: NGUYỄN NHƯ Ý HỒ HOÀNG SƠN VÕ TRẦN THẢO VÂN NGUYỄN ĐẶNG PHI THUẦN SKL011035 Tp Hồ Chí Minh, tháng 7/2023 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP HCM KHOA CƠ KHÍ ĐỘNG LỰC ĐỒ ÁN TỐT NGHIỆP *** NGHIÊN CỨU THỰC NGHIỆM QUÁ TRÌNH LÀM MÁT PHỤ CHO HỆ THỐNG ĐIỀU HỊA KHƠNG KHÍ CO2 BẰNG NƯỚC NGƯNG GVHD: PGS TS ĐẶNG THÀNH TRUNG HỌ VÀ TÊN MSSV NGUYỄN NHƯ Ý HỒ HOÀNG SƠN VÕ TRẦN THẢO VÂN NGUYỄN ĐẶNG PHI THUẦN 19147274 19147231 19147269 19147245 Tp Hồ Chí Minh, tháng 07 năm 2023 Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp GVHD: PGS.TS Đặng Thành Trung GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp 5.1.3 Hệ thống lạnh sử dụng đường ống làm mát phụ mm Sau chạy thực nghiệm với ống phụ mm hệ thống lấy số liệu để tính tốn, ta có thơng số trạng thái hệ thống Hình 5.9 Bảng 5.3: Hình 5.9 Số liệu thực nghiệm từ hệ thống sử dụng ống làm mát phụ mm - Lưu lượng môi chất ���2 = 0,670 l/m = 0,007 kg/s - Lưu lượng gió dàn lạnh : ��� = 0,2 kg/s - Áp suất ngưng tụ tk: pk = 71,02 bar - Áp suất bay t0: p0 = 46,03 bar - Với áp suất ngưng tụ pk (CO2) = 71,02 bar áp suất bay p0 = 46,03 bar Tra bảng tính chất vật lý lỏng bão hòa CO2:[10] Nhiệt độ ngưng tụ tk (CO2) = 29,36oC Nhiệt độ bay hơi: t0(CO2) = 10,91oC Nhiệt độ nhiệt tqn= t0+ ∆tqn , với ∆tqn = 11oC ( Độ nhiệt diện tích trao đổi nhiệt thiết bị bay lớn, độ nhiệt từ dàn lạnh máy nén gần tổn thất áp suất nên độ nhiệt 0oC), tqn = 10,91 + 11 = 21,91oC Nhiệt độ lạnh ∆tql = 0,1oC (Độ lạnh môi chất nhả nhiệt với nước ngưng đường ống làm mát phụ trước vào van tiết lưu) 49 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp Ta có tỉ số nén chu trình : �= �� (��2 ) �� (��2 ) = 72,02 46,03 = 1,56 Hình 5.10 Đồ thị p-h chu trình lạnh sử dụng ống phụ mm Bảng 5.3 Thông số trạng thái điểm nút Điểm Nhiệt độ Áp suất Entanpy Entropy Thể tích t(°C) p(bar) h(kJ/kg) s(kJ/kg.K) v(m3/kg) 21,91 46,03 444,4 1,85 0,008 55,7 71,02 462,09 1,85 0,006 3’ 29,36 71,02 299,92 1,33 0,0016 29,26 70,86 299,07 1,32 0,0016 10,91 46,03 299,07 1,09 _ 10,91 46,03 422,34 1,78 0,0072 50 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp Năng suất lạnh riêng: q0 = h1 – h4 = 444,4 – 299,07 = 145,33 kJ/kg Công máy nén riêng: l = h2 – h1 = 462,09 – 444,4 = 17,69 kJ/kg Công suất nhiệt riêng: qk = h2 – h3 = 462,09 – 299,07 = 163,02 kJ/kg Công suất lạnh phía mơi chất: Cơng máy nén lạnh: Cơng suất nhiệt: Q0 = ���2 q0 = 0,007.145,33 = kW L = ���2 l = 0,007.17,69 = 0,12 kW Qk = ���2 qk = 0,007.163,02 = 1,14 kW Hệ số lạnh chu trình: ε= �0 � = 145,33 17,69 = 8,2 Hình 5.11 Biểu đồ hệ số làm lạnh sử dụng ống phụ mm 51 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp COP hệ thống lạnh sử dụng đường ống phụ mm làm mát nước ngưng thiết bị vận hành hệ thống nhiều tháng điều kiện áp suất ngưng tụ pk = 71,02 bar áp suất bay p0 = 46,03 bar lưu lượng môi chất ���2 = 0,007 kg/s Dựa vào biểu đồ Hình 5.11, ta thấy COP dao động trung bình 8,2 Hệ số làm lạnh hệ thống có xu hướng tăng dần từ lúc khởi động đến thứ giảm từ lúc thứ sau tăng từ thứ trì ổn định đến thứ Điều cho thấy trưa hệ số làm lạnh có xu hướng giảm lúc sáng chiều Dàn lạnh sử dụng hệ thống dàn lạnh có cơng suất 2,65 kW với tốc độ gió lớn đạt m/s thu bảng số liệu đồ thị Hình 5.12 Hình 5.12 Biểu đồ nhiệt độ gió vào dàn lạnh sử dụng đường ống phụ mm Khi vận hành lấy số liêu hệ thống, ta thu nhiệt độ gió vào nhiệt độ gió dàn lạnh biểu diển Hình 5.12 trình hoạt động hệ thống phịng có kích thước lớn, khơng kín dụng cụ đo có phần hạn chế nên nhiệt độ đo không đạt yêu cầu khơng xác hồn tồn Nhưng từ biểu đồ ta thấy nhiệt độ gió vào nhiệt độ gió dàn lạnh dao động khoảng từ 9oC 10oC Tuy chưa hiệu mặt nhiệt độ kết luận hệ thống hoạt động ổn định 52 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp 5.2 So sánh kết thực nghiệm 5.2.1 Kết tính tốn trường hợp thực nghiệm Sau tiến hành tính tốn lý thuyết chạy thực nghiệm hệ thống điều hịa khơng khí sử dụng mơi chất CO2 kết hợp đường ống làm mát phụ vào hệ thống ta thu số liệu Bảng 5.4 Bảng 5.4 Thông số lý thuyết q trình thực nghiệm Thơng số Hệ thống Lưu lượng ���2 (kg/s) Năng suất lạnh phía Qo mơi chất (kW) Năng suất lạnh phía Qo khơng khí (kW) Hệ số làm lạnh Áp suất Pk (bar) Áp suất Po (bar) Tận dụng nước ngưng Lý thuyết Chưa tận dụng nước ngưng 0,007 0,007 0,007 0,007 0,95 0,99 1 0,92 0,92 0,92 7,2 8,12 8,16 8,21 72 45 72,29 48 71,59 46,48 71,02 46,03 Ống 10mm Ống 8mm Nhận xét: Những thông tin từ Bảng 5.4 cho thấy lý thuyết tính tốn cịn có sai số chênh lệch với khơng đáng kể chấp nhận Các q trình thực nghiệm kết hợp đường ống phụ làm mát nước ngưng hệ thống cải thiện hệ số làm lạnh tăng đáng kể có tính hiệu ứng dụng cao Nhưng cần phải cải thiện thêm để đạt hệ số làm lạnh cao Quan sát Bảng 5.4 so sánh thông số thực nghiệm ta thấy được: - Lưu lượng Gco2 tính tốn lý thuyết cao 5% so với lưu lượng Gco2 hệ thống thực tế - Hệ số COP hệ thống thực tế cao tính tốn lý thuyết hệ thống dụng đường ống làm mát phụ mm - Công suất làm lạnh hệ thống thực nghiệm gần tính tốn lý thuyết nên ta dễ dàng xác định khả hoạt động hệ thống - Công suất làm lạnh phía gió trung bình hệ thống thực nghiệm thấp 3,8% so với lý thuyết 53 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp - Tổn thất áp suất (bar) lý luyết so với tổn thất áp suất (bar) thực nghiệm không đáng kể - Những thông tin cho thấy hệ thống thực nghiệm cho kết tốt, xấp xỉ sai số chấp nhận 5.2.2 Hệ số làm lạnh trường hợp thực nghiệm Khi chạy thực nghiệm hệ thống điều hịa khơng khí sử dụng môi chất CO2 chưa tận dụng nước ngưng làm mát trường hợp sử dụng đường ống làm mát phụ mm đướng ống làm mát phụ 10 mm nhằm sử dụng nhiệt từ nước ngưng dàn lạnh để giúp tăng độ lạnh với mục đích nâng cao suất lạnh hệ thống cao Sau tiến hành thực nghiệm ta thu số liệu thể Bảng 5.5 Bảng 5.5 Hệ số làm lạnh trường hợp thực nghiệm Thời gian COP COP (10mm) COP (8mm) 1h 2h 3h 4h 5h 6h 8,10 8,12 8,12 8,11 8,10 8,12 8,14 8,16 8,16 8,15 8,15 8,16 8,20 8,20 8,21 8,21 8,20 8,21 Từ thông số trình bày Bảng 5.5 ba trình chưa tận dụng nước ngưng với sử dụng ống làm mát phụ, ta thấy hệ số làm lạnh cải thiện so với lý thuyết tính tốn Các trường hợp thực nghiệm trì ổn định không chênh lệch nhiều suốt 6h Đặc biệt kết hợp đường ống làm mát phụ vào hệ thống kết COP khả quan hệ thống chưa tận dụng nhiệt từ nước ngưng Tuy hai đường ống phụ làm mát mm 10 mm cải thiện COP đường ống mm cao đạt COP = 8,21 5.2.3 Nhiệt độ môi trường nhiệt độ nước bể quạt làm mát Sau vận hành hệ thống ghi nhận số liệu vài tháng ta có nhiệt độ mơi trường nhiệt độ nước bể hệ thống chưa tận dụng làm mát phụ nước ngưng hệ thống sử dụng ống làm mát phụ trường hợp Từ số liệu thực nghiệm, nhóm tiến hành so sánh chọn lọc số liệu nhiều trường hợp với nhiệt độ môi trường nước giải nhiệt cho thiết bị trao đổi nhiệt chạy thực nghiệm Ta thu số liệu Bảng 5.6 Hình 5.13 sau: 54 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp Bảng 5.6 Nhiệt độ môi trường nhiệt độ nước bể Chưa tận dụng nước ngưng Thời gian Nhiệt độ nước bể (℃) 1h 2h 3h 4h 5h 6h 28 29 29,2 29 29 29 29,2 Nhiệt độ môi trường (℃) 31 31 31 31,5 31,4 31 31 Ống 10 mm Nhiệt độ nước bể (℃) 28 28,6 29 29 28,8 28,9 29 Nhiệt độ môi trường (℃) 31 31 32 32 32 32 32 Ống mm Nhiệt độ nước bể (℃) 28 28,4 28,5 28,3 28,5 28,6 28,5 Nhiệt độ môi trường (℃) 31 32 32 31,5 32 32 32 Trong Hình 5.13 biểu diễn hiệu nhiệt độ mơi trường nhiệt độ nước bể quạt làm mát ba trường hợp, khơng có chênh lệch nhiều trường hợp nhiệt độ dao động khoảng 2oC - 3oC Nhưng ta nhận thấy việc kết hợp đường ống phụ làm mát độ chênh lệch nhiệt độ môi trường nước giải nhiệt thấp hệ thống bình thường Đặc biệt đường ống làm mát phụ mm độ chệnh lệch cao 3,5oC Hình 5.13 Biểu đồ hiệu nhiệt độ mơi trường nhiệt độ nước bể 55 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp 5.2.4 Nhiệt độ đường ống làm mát Khi lắp đặt đường ống làm mát phụ 10 mm đường ống làm mát phụ mm, sau tiến hành thực nghiệm hệ thống ta thu số liệu Bảng 5.7 Hình 5.14 Bảng 5.7 Nhiệt độ đường ống làm mát Nhiệt độ CO2 với ống 10 mm Nhiệt độ CO2 với ống mm Nhiệt độ CO2 với ống mm cách nhiệt Thời gian Trước làm mát (℃) Sau làm mát (℃) Trước làm mát (℃) Sau làm mát (℃) Trước làm mát (℃) Sau làm mát (℃) 1h 2h 3h 4h 5h 6h 29,2 29,6 29,6 29,8 30 29,8 29,1 29,5 29,6 29,8 29,9 29,8 29,5 30 30,1 30,3 30,2 30,1 29,4 29,9 30 30 30 30 28,4 29,3 29,3 29,4 29,5 30,2 28,2 29,1 29 29 29,1 29,8 Thơng qua Hình 5.14 ta thấy nhiệt độ đường ống làm mát phụ hai đầu đầu vào có thay đổi chưa đáng kể Cụ thể, với đường ống phụ 10 mm độ chênh nhiệt độ chưa có thay đổi, cao có 0,1oC Cịn đường ống phụ mm độ chênh lệch dao động từ 0,1oC - 0,2oC có thay đổi nhiều sử dụng đường ống phụ 10 mm chưa đạt độ lạnh mong muốn nên ta tiến hành bọc cách nhiệt vào đường ống phụ mm thu kết độ chênh nhiệt độ đạt cao 0,4oC ổn định từ thứ Hình 5.14 Biểu đồ độ chênh lệch nhiệt độ đường ống làm mát phụ 56 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp Biểu đồ Hình 5.14, ta thấy khác biệt tận dụng nước ngưng giải nhiệt cho ống đồng đường kính mm do: - Khi dùng ống 10 mm, khoảng cách ống 10 mm ống đồng lớn nước ngưng lúc không giải nhiệt triệt để cho ống đồng, phần nước ngưng chủ yếu trôi thẳng không tiếp xúc nhiều với ống đồng nên hiệu giải nhiệt - Khi dùng ống 8mm, khoảng cách ống mm ống đồng vừa phải nước ngưng lúc tận dụng giải nhiệt triệt để so với ống 10 mm, lúc nước ngưng ống đồng tiếp xúc với nhiều nên hiệu giải nhiệt cải thiện rõ ống 10 mm - Khi sử dụng ống mm có bọc cách nhiệt, lúc hiệu tận dụng nước ngưng để giải nhiệt ổn so với hai phương án Lúc thất thoát nhiệt độ nước ngưng hạn chế Do sử dụng ống mm có bọc cách nhiệt tận dụng nước ngưng để giải nhiệt tối ưu so với hai phương án 5.2.5 Kết nhiệt độ thực nghiệm dàn lạnh Sử dụng máy ảnh hồng ngoại Fluke TiS10 để có ảnh thực nghiệm, từ ta quan sát rõ độ nhiệt dàn lạnh Cho thấy nhiệt độ môi chất lạnh vào dàn trung bình khoảng 12,4oC cịn nhiệt độ mơi chất lạnh dàn lạnh trung bình khoảng 22oC Ta suy độ q nhiệt trung bình thực tế dàn 11oC thể Hình 5.15 Hình 5.16 Kết phù hợp với tính tốn lý thuyết Hình 5.15 Nhiệt độ mơi chất vào môi chất dàn lạnh 57 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp Hình 5.16 Nhiệt độ bên cánh tản nhiệt dàn lạnh Cũng máy ảnh hồng ngoại Fluke TiS10, chúng em kiểm tra lại nhiệt độ gió khỏi dàn lạnh thực tế, trung bình khoảng 24,5oC, nhiệt độ phù hợp với yêu cầu mà nhóm đạt ban đầu Hình 5.17 Hình 5.17 Nhiệt độ gió dàn lạnh 58 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp CHƯƠNG 6: KẾT LUẬN VÀ KIẾN NGHỊ 6.1 Kết luận Qua thời gian tiến hành thực nghiên cứu đề tài “Nghiên cứu thực nghiệm trình làm mát phụ cho hệ thống điều hịa khơng khí CO2 nước ngưng ” Nhóm thực việc tính tốn, kiểm tra vận hành hệ thống, để tiến hành cho q trình xây dựng mơ hình thực nghiệm Từ kết thực nghiệm nhiều trường hợp với điều kiện áp suất ta thu thông số thực nghiệm sau: Đối với hệ thống chưa lắp đường ống làm mát phụ: Nhiệt độ phòng dao động khoảng 29oC-30oC, hệ số làm lạnh khoảng 8,12 với lưu lượng môi chất ổn định mức 0,007 kg/s Sau vận hành hệ thống thu kết ổn định, nhóm tiến hành thiết kế đường ống làm mát phụ nước ngưng vào hệ thống với hai loại đường kính ống mm 10 mm để thực nghiệm tính hiệu hai loại đường ống Qua nhiều lần vận hành hệ thống kết thu kết hợp đường ống phụ làm mát nước ngưng khả quan : Đối với đường ống phụ làm mát 10 mm: Nhiệt độ phòng dao động khoảng 29oC-30oC so với hệ thống lúc bình thường, hệ số làm lạnh khoảng 8,16 lưu lương môi chất ổn định mức 0,007 kg/s Tuy có cải thiện so với hệ thống lúc bình thường chưa đáng kể Đối với đường ống phụ làm mát mm: Nhiệt độ phòng dao động khoảng 28,5oC-29oC so với hệ thống lúc bình thường, hệ số làm lạnh khoảng 8,21 lưu lương môi chất ổn định mức 0,007 kg/s Ta thấy nhiệt độ phòng hệ số làm lạnh thu khả quan ổn định Từ kết thực nghiệm với hai đường ống phụ so với kết thực nghiệm trước kết đường ống phụ đường kính mm hiệu Qua kết thực nghiệm đạt thấy việc sử dụng đường ống phụ làm mát nước ngưng có hiệu Tuy nhiên, cịn số hạn chế định, đặc biệt thay đổi điều kiện mơi trường bên ngồi 6.2 Kiến nghị Đề tài hoàn thành mục tiêu đề tiến hành q trình thực nghiệm hệ thống thực tế để lấy thông số nhiệt độ, áp suất lưu lượng sau so sánh với kết kết hợp đường ống phụ làm mát nước ngưng, ngồi hệ thống cịn có sử 59 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp dụng xử lí tín hiệu MX100 giúp cho việc lấy thơng số đạt hiệu xác cao Từ so sánh đưa nhận xét kết luận Mặc dù nhóm có nhiều cố gắng trình nghiên cứu thực hạn chế thời gian, kiến thức chuyên sâu nên đề tài dừng lại phạm vi tính tốn, thiết kế lý thuyết, so sánh với số liệu thực nghiệm thu Vì vậy, đề tài cịn có khó khăn định cần khắc phục tiếp tục nghiên cứu sâu Ngoài ra, việc kết hợp đường ống phụ làm mát nước ngưng hiệu chưa thực đạt mong muốn cần áp dụng thêm trường hợp khác để cải thiện hệ số làm lạnh thêm cho hệ thống sau Do nhóm nghiên cứu cần góp ý đánh giá thầy để giúp đề tài hồn thiện Với hạn chế nêu trên, nhóm nghiên cứu xin đưa kiến nghị để tiếp tục phát triển hoàn thiện đề tài cách tốt nhất: Kế thừa liệu tính tốn nhóm để tiến hành cải tiến giúp cho hệ thống đạt thêm nhiều kết khả quan tương lai 60 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp TÀI LIỆU THAM KHẢO [1] PGS TS Đặng Thành Trung TS Đoàn Minh Hùng, Nghiên cứu thực nghiệm thông số nhiệt động trình giản nở tiết lưu hệ thống điều hồ khơng khí dùng mơi chất CO2, Trường Đại học Sư phạm Kỹ thuật Tp HCM, năm 2016 [2] Ketdoan V Chau, Tronghieu Nguyen, and Thanhtrung Dang, Numerical Simulation on Heat Transfer Phenomena in Microchannel Evaporator of A CO2 Air Conditioning System, American Journal of Engineering Research, Issue 2, 2017 [3] Đặng Thanh Phong, Huỳnh Trần Trúc Vương, Phạm Tuấn Kiệt, Nguyễn Tấn Phát, Trần Thành Phát, Nghiên cứu mô số thực nghiệm q trình làm mát mơi chất CO2 thiết bị ngưng tụ nước, năm 2023 [4] Gustav Lorentzen, Jostein Pettersen, A new efficient and environmentally benign system for car air-conditioning, 2003 [5] Man-Hoe Kim, Jostein Pettersen, Clark W Bullard, Fundamental process and system design issues in CO2 vapor compression systems, 2003 [6] Simarpreet Singh, MS Dasgupta, Evaluation of research on CO2 trans-critical work recovery expander using multi attribute decision making methods, 2016 [7] PanelIgnacio Peñarrocha a, Rodrigo Llopis b, Luis Tárrega a, Daniel Sánchez b, Ramón Cabell, A new approach to optimize the energy efficiency of CO2 transcritical refrigeration plants, 2014 [8] Laura Nebot-Andr'es, Daniel Sanchez, Daniel Calleja-Anta, Ramon'Cabello, Rodrigo Lopis, Experimental determination of the optimum intermediate and gas-cooler pressures of a commercial transcritical CO2 refrigeration plant with parallel compression, 2021 [9] Pradeep Bansal, A review e Status of CO2 as a low temperature refrigerant: Fundamentals and R&D opportunities, The University of Auckland, Private Bag, 92019 Auckland, New Zealand, 2021 [10] Fang Liu, Eckhard A, Groli, Daqing Li, Modeling study of an ejector expansion residential CO2 air conditioning system, 2021 [11] ASHRAE Handbook, Fundamentals, 2017 61 GVHD: PGS.TS Đặng Thành Trung Đồ án tốt nghiệp [12] PGS TS Nguyễn Đức Lợi, Sách Hướng Dẫn Thiết Kế Hệ Thống Lạnh, NXB Khoa Học Và Kĩ Thuật Hà Nội, Năm 2005 [13] PGS TS Hồng Đình Tín, Cơ sở Truyền nhiệt Thiết kế thiết bị trao đổi nhiệt, NXB Đại Học Quốc Gia TP.HCM, năm 2013 [14] Tronghieu Nguyen and Thanhtrung Dang, The Effect of Fin Shape on the Heat Transfer and the Solution Time of a Microchannel Evaporator in a CO2 Air Conditioning System—A Numerical Investigation, 2022 [15] Philip Kosky, Robert Balmer, William Keat, Mechanical Engineering, 2020 [16] E Edge, Convective Heat Transfer Coefficients Table Chart, 2020-2023 [17] GS Lê Xn Hịa, Giáo Trình Kỹ Thuật Lạnh, Thành Phố Hồ Chí Minh năm 2007 [18] PGS TS Hồng Đình Tín PGS TS Bùi Hải, Bài tập Nhiệt động học kỹ thuật Truyền nhiệt, NXB Đại Học Quốc Gia TP.HCM, năm 2018 62 S K L 0