Nghiên cứu lý thuyết và thực nghiệm ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm thùng giao thực phẩm

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA CƠ KHÍ BỘ MƠN CƠNG NGHỆ NHIỆT LẠNH -o0o - NGUYỄN TIẾN CẢNH NGHIÊN CỨU LÝ THUYẾT VÀ THỰC NGHIỆM ỐNG NHIỆT TẬN DỤNG NHIỆT THẢI CỦA XE MÁY ĐỂ GIỮ ẤM THÙNG GIAO THỰC PHẨM Chuyên ngành: KỸ THUẬT NHIỆT Mã số: 60.52.80 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, THÁNG 07 NĂM 2014 CƠNG TRÌNH ÐƯỢC HỒN THÀNH TẠI TRƯỜNG ÐẠI HỌC BÁCH KHOA –ÐHQG – HCM Cán hướng dẫn khoa học: GS.TS LÊ CHÍ HIỆP Cán chấm nhận xét 1: TS BÙI NGỌC HÙNG Cán chấm nhận xét 2: TS BÙI TRUNG THÀNH Luận văn thạc sĩ bảo vệ Trường Ðại học Bách Khoa, ÐHQG Tp HCM ngày 15 tháng nam 2014 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: Chủ tịch: TS NGUYỄN VĂN TUYÊN Thư ký: TS HÀ ANH TÙNG Ủy viên – Phản biện 1: TS BÙI NGỌC HÙNG Ủy viên – Phản biện 2: TS BÙI TRUNG THÀNH Ủy viên: GS.TS LÊ CHÍ HIỆP Xác nhận Chủ tịch Hội đồng đánh giá luận văn Trưởng khoa quản lý chuyên ngành sau luận văn sửa chữa (nếu có) CHỦ TỊCH HỘI ÐỒNG TRƯỞNG KHOA CƠ KHÍ i ĐẠI HỌC QUỐC GIA TP.HCM TRƢỜNG ĐẠI HỌC BÁCH KHOA _ CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập – Tự Do – Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên: Năm sinh: Chuyên ngành: NGUYỄN TIẾN CẢNH 1974 CÔNG NGHỆ NHIỆT Phái: Nơi sinh: MSHV: Nam BRVT 11060417 Tên đề tài: NGHIÊN CỨU LÝ THUYẾT VÀ THỰC NGHIỆM ỐNG NHIỆT TẬN DỤNG NHIỆT THẢI CỦA XE MÁY ĐỂ GIỮ ẤM THÙNG GIAO THỰC PHẨM Nhiệm vụ luận văn - Khảo sát loại ống nhiệt nói chung đặc biệt loại ống nhiệt sử dụng nước để gia nhiệt - Nghiên cứu ống nhiệt tận dụng nhiệt thải - Tính tốn, thiết kế, chế tạo ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm cho thùng giao thực phẩm quy mô nhỏ - Thực nghiệm mơ hình để xây dựng thơng số ảnh hưởng đến hiệu suất ống nhiệt - Rút kết luận kiến nghị Ngày giao nhiệm vụ : 21/01/2014 Ngày hoàn thành nhiệm vụ: 21/06/2014 Cán hƣớng dẫn : GS TS Lê Chí Hiệp Chủ nhiệm mơn Tp.HCM, ngày … tháng …… năm…… Cán hƣớng dẫn GS.TS Lê Chí Hiệp GS.TS Lê Chí Hiệp Trƣởng khoa khí ii LỜI CẢM ƠN Học viên xin bày tỏ lòng biết ơn sâu sắc:  Ðến Thầy hướng dẫn – GS.TS Lê Chí Hiệp – dành nhiều thời gian định hướng, hướng dẫn, sửa chữa giúp đỡ học viên để học viên hoàn thành nội dung nghiên cứu luận văn  Ðến Thầy – TS Bùi Trung Thành – động viên, giúp đỡ, tạo nhiều điều kiện thuận lợi để học viên hoàn thành luận văn  Ðến quý Thầy Cô Bộ môn Công Nghệ Nhiệt Lạnh – Trường ÐH Bách Khoa – ÐHQG TP.HCM nhiệt tình giảng dạy, truyền đạt kiến thức, kinh nghiệm suốt thời gian học viên học cao học  Ðến quý Thầy, Cô đồng nghiệp Khoa Công Nghệ Nhiệt Lạnh – Trường ÐH Công Nghiệp TP.HCM giúp đỡ, động viên, chia sẻ công việc thời gian học viên học cao học học viên thực luận văn  Ðến thành viên gia đình người bạn thân đồng hành, giúp đỡ, động viên học viên suốt thời gian qua Học Viên NGUYỄN TIẾN CẢNH iii TÓM TẮT NỘI DUNG LUẬN VĂN Luận văn Nghiên cứu lý thuyết thực nghiệm ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm cho thùng giao thực phẩm bao gồm nội dung sau: - Nghiên cứu sở lý thuyết liên quan đến tính tốn, thiết kế ống nhiệt từ cơng trình nghiên cứu ngồi nước - Xây dựng sở lý thuyết viết chương trình thiết kế ống nhiệt tận dụng nhiệt thải xe máy chương trình khác tính tỷ lệ nạp theo tải, xác định phụ tải nhiệt theo tỉ lệ nạp ống nhiệt chương trình khác phục vụ cho việc đánh giá hoạt động ống nhiệt nhằm lựa chọn thông số tối ưu ống nhiệt - Chế tạo ống nhiệt tận dụng nhiệt thải xe máy tiến hành nghiên cứu thực nghiệm sau: - Nghiên cứu, thực nghiệm tỷ lệ nạp tối ưu môi chất (nước) làm việc bên ống nhiệt - Nghiên cứu thực nghiệm đánh giá hiệu giữ ấm thực phẩm ống nhiệt tận dụng nhiệt thải học viên chế tạo - Rút kết luận kiến nghị ABSTRACT The thesis “Theoretical & Experimental Study on Heat Pipe used vehicle wast heat to keep warm for dilivery foods” consists of the following main items: - Synthesizing the research results from the related references - Writing the theoretical to design the working parameters of heat pipe then writing the program to design the heat pipe by Matlab program In addition, writing some sub - programs such as calculate the filling ratio by load, calculate the heat load by filling ratio, as well as other sub programs to find out the optimine parameters of heat pipe - Manufacturing the heat pipe use vehicle wast heat And searching some iterms as follow: + Studying the optimine filling ratio of working fluid (pure water) inside the heat pipe + Studying the efficiency of manufactured heat pipe thought keep warm efficiency foods in the delivery tank - Conclutions and suggestions iv LÝ LỊCH TRÍCH NGANG Họ tên: NGUYỄN TIẾN CẢNH Năm sinh: 1974 Nơi sinh: Phƣớc Tỉnh, Long Điền, Bà Rịa – Vũng Tàu Địa liên lạc: Bộ môn Kỹ Thuật Lạnh – Khoa Công nghệ Nhiệt lạnh – Trƣờng Đại học Công nghiệp TP.HCM QUÁ TRÌNH ĐÀO TẠO: Từ năm 1994 – 1999: Sinh viên trường Đại học Thủy Sản Nha Trang - Chuyên ngành: Công nghệ chế biến Thủy Sản Từ năm 1995 – 2000: Sinh viên trường Đại học Tổng hợp TP.HCM, trường Đại Học Khoa Học Xã Hội Nhân Văn - Chuyên ngành: Ngoại ngữ Tiếng Anh Từ năm 2008 – 2010: Học viên lớp Bồi dưỡng Sau Đại Học, trường Đại học Công Nghiệp TP.HCM - Chuyên ngành: Quản trị kinh doanh Từ năm 2011 đến năm 2013: Học viên Cao học trường Ðại Học Học Bách Khoa – ÐHQG TP HCM - Chuyên ngành: Công Nghệ Nhiệt Q TRÌNH CƠNG TÁC: Từ năm 2000 đến 2003: Trưởng Bộ môn Kỹ Thuật Lạnh, khoa Điện – Điện tử - Điện lạnh, Trường TH KT NV Nguyễn Hữu Cảnh, Quận 7, TP.HCM Từ năm 2004 đến nay: Giảng viên Bộ môn Kỹ Thuật Lạnh, khoa Công Nghệ Nhiệt Lạnh, trường Đại học Công nghiệp TP.HCM Tp.HCM, ngày 15 tháng 07 năm 2014 Người khai NGUYỄN TIẾN CẢNH v LỜI CAM ĐOAN Học viên xin cam đoan nội dung trình bày luận văn học viên tìm hiểu nghiên cứu Các nội dung trích dẫn từ nguồn tài liệu khác có ghi rõ nguồn gốc Tác giả NGUYỄN TIẾN CẢNH vi MỤC LỤC Chƣơng Tựa đề Nhiệm vụ luận văn Lời cảm ơn Tóm tắt luận văn Lý lịch trích ngang Lời cam đoan Mục lục Danh mục ký hiệu chữ viết tắt Danh mục bảng Danh mục hình vẽ đồ thị Mở đầu Tổng quan 1.1 Nghiên cứu ống nhiệt 1.1.1 Tổng quan ống nhiệt 1.1.2 Ống nhiệt vòng 1.2 Các nghiên cứu ống nhiệt giới 1.3 Các nghiên cứu ống nhiệt nước Cơ sở lý thuyết 2.1 Khảo sát phụ tải hệ thống 2.1.1 Nhiệt lượng cần để trì nhiệt độ thực phẩm, Q1 2.1.2 Nhiệt lượng cần để trì nhiệt độ thùng, Q2 2.1.3 Tổn thất nhiệt môi trường 2.1.4 Một số vấn đề khác liên quan 2.2 Cơ sở truyền nhiệt 2.2.1 Truyền nhiệt sôi 2.2.2 Truyền nhiệt ngưng 2.3 Giới hạn làm việc ống nhiệt 2.4 Các vấn đề khác liên quan đến thiết kế ống nhiệt 2.4.1 Môi chất làm việc bên ống nhiệt 2.4.2 Vật liệu chế tạo ống nhiệt 2.4.3 Tỷ lệ nạp ống nhiệt 2.5 Kết luận chương Tính tốn, thiết kế ống nhiệt tận dụng nhiệt thải xe máy 3.1 Nguyên lý hoạt động hệ thống 3.2 Số liệu thiết kế 3.3 Tính tốn, thiết kế 3.3.1 Trường hợp 3.3.2 Trường hợp Trang i ii iii iv v vi viii x xi xii 1 22 24 27 27 27 27 28 29 30 30 32 34 36 36 38 39 41 42 42 43 43 43 48 vii 3.3.3 Trường hợp 48 3.4 Thiết lập phần mềm tính tốn thiết kế 50 3.4.1 Chương trình thiết kế ống nhiệt 50 3.4.2 Chương trình tính tỷ lệ nạp theo tải 51 3.5 Phân tích số ảnh hưởng hệ thống theo tính tốn lý thuyết 52 3.5.1 Tỷ lệ nạp thay đổi theo tải 52 3.5.2 Khối lượng thực phẩm thay đổi theo tỷ lệ nạp 54 3.6 Kết nghiên cứu lý thuyết 55 Nghiên cứu thực nghiệm ống nhiệt 56 4.1 Gia công chế tạo 56 4.1.1 Đặc điểm 57 4.1.2 Các bước gia công ống nhiệt 57 4.1.3 Hút chân không ống nhiệt 57 4.1.4 Nạp môi chất vào ống nhiệt 59 4.2 Phương pháp thực nghiệm 59 4.2.1 Phương pháp nghiên cứu 59 4.2.2 Các thiết bị đo phương pháp đo 60 4.2.3 Phương pháp tiến hành thí nghiệm mơ hình 63 4.3 Kết thực nghiệm 64 4.3.1 Thực nghiệm khảo sát hoạt động hệ thống theo tỷ lệ nạp 58% 64 4.3.2 Thực nghiệm khảo sát hoạt động hệ thống theo tỷ lệ nạp 68% 65 4.3.3 Thực nghiệm khảo sát hoạt động hệ thống theo tỷ lệ nạp 78% 67 4.3.4 Thực nghiệm khảo sát hoạt động hệ thống theo tỷ lệ nạp 88% 68 4.3.5 Thực nghiệm khảo sát hoạt động hệ thống theo tỷ lệ nạp 98% 69 4.3.6 Đánh giá hoạt động ống nhiệt thông qua nhiệt trở 71 4.3.7 Đánh giá ảnh hưởng tỷ lệ nạp đến hoạt động hệ thống 77 4.4 Kết nghiên cứu thực nghiệm 81 Kết luận kiến nghị 82 5.1 Kết luận chung 82 5.1.1 Các kết đạt 82 5.1.2 Một số hạn chế 82 5.2 Kiến nghị hướng nghiên cứu 83 Phụ lục 84 Tài liệu tham khảo 114 viii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT Ký hiệu a C D d F G H kv kl L M’ n Nq p Pd Ptđ Pr Q Q1 Q2 Q3 Q31 Q32 q r R Re t Te Tc Tth Tk V W  Ý nghĩa Thứ nguyên Hệ số khuếch tán nhiệt m2/s Nhiệt dung riêng J/kgK Đường kính ngồi mm Đường kính mm Diện tích m2 Lưu lượng khối lượng kg/s Chiều cao mm Hệ số truyền nhiệt qua vách phẳng W/m2.K Hệ số truyền nhiệt đường W/m.K Chiều dài mm Hệ số Merit kg/(K3/4.s5/2) Số ống dàn ngưng ống Tiêu hao điện riêng Wh/kg ẩm Áp suất bão hịa mơi chất sơi bar Áp suất dư đo bình bay cmHg Áp suất tuyệt đối tính từ áp suất dư cmHg Tiêu chuẩn Prandtl Tổng phụ tải nhiệt hệ thống kJ Nhiệt lượng cần để trì nhiệt độ thực phẩm kJ Nhiệt lượng cần để trì nhiệt độ thùng kJ Tổn thất nhiệt môi trường kJ Tổn thất nhiệt từ thùng môi trường kJ Tổn thất nhiệt đường ống dẫn kJ Mật độ dịng nhiệt kJ/kg Nhiệt ẩn hóa kJ/kg Nhiệt trở C/W Tiêu chuẩn Reynolds Nhiệt độ C Nhiệt độ trung bình đo bề mặt thiết bị bay C Nhiệt độ trung bình đo bề mặt thiết bị ngưng tụ C Nhiệt độ trung bình đo thùng chứa thực phẩm C Nhiệt độ trung bình đo đầu ống khói C Thể tích m3 Chiều rộng mm Hệ số tỏa nhiệt đối lưu W/m2.K CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH % Hints: get(hObject,'String') returns contents of edit11 as text % str2double(get(hObject,'String')) returns contents of edit11 as a double % - Executes during object creation, after setting all properties function edit11_CreateFcn(hObject, eventdata, handles) % hObject handle to edit11 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit12_Callback(hObject, eventdata, handles) % hObject handle to edit12 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit12 as text % str2double(get(hObject,'String')) returns contents of edit12 as a double % - Executes during object creation, after setting all properties function edit12_CreateFcn(hObject, eventdata, handles) % hObject handle to edit12 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit13_Callback(hObject, eventdata, handles) % hObject handle to edit13 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit13 as text % str2double(get(hObject,'String')) returns contents of edit13 as a double % - Executes during object creation, after setting all properties function edit13_CreateFcn(hObject, eventdata, handles) % hObject handle to edit13 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end Phụ lục Trang 101 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH function edit14_Callback(hObject, eventdata, handles) % hObject handle to edit14 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit14 as text % str2double(get(hObject,'String')) returns contents of edit14 as a double % - Executes during object creation, after setting all properties function edit14_CreateFcn(hObject, eventdata, handles) % hObject handle to edit14 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit15_Callback(hObject, eventdata, handles) % hObject handle to edit15 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit15 as text % str2double(get(hObject,'String')) returns contents of edit15 as a double % - Executes during object creation, after setting all properties function edit15_CreateFcn(hObject, eventdata, handles) % hObject handle to edit15 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit16_Callback(hObject, eventdata, handles) % hObject handle to edit16 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit16 as text % str2double(get(hObject,'String')) returns contents of edit16 as a double % - Executes during object creation, after setting all properties function edit16_CreateFcn(hObject, eventdata, handles) % hObject handle to edit16 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called Phụ lục Trang 102 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes on button press in pushbutton1 function pushbutton1_Callback(hObject, eventdata, handles) % hObject handle to pushbutton1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) format short %Khoi luong thuc pham, kg: mtp=str2double(get(handles.edit16,'String')); %du lieu mac dinh tinh toan: %Nhiet dung rieng cua thuc pham, kJ/kgK: ctp=4.34; %Nhiet ban dau cua thuc pham, degC: ttp=50; %Nhiet bên thung, degC: tf1=58; %Nhiet moi truong, degC: tf2=30; %Nhiet khoi, degC: tk1=120; %Nhiet khoi dau ong khoi, degC: tk2=70; %Van toc khoi, m/s: vk=5; %Chieu dai ben thung, m: Ltr=0.6; %Chieu rong ben thung, m: Wtr=0.4; %Chieu cao ben thung, m: Htr=0.6; %Chieu day lop cach nhiet thung, m: Tcn1=0.05; %Chieu dày ong dan, m: Tod=0.0007; %Chieu dày ong cach nhiet ong dan, m: Tcn2=0.013; %Chieu dài duong ong dan moi chat, m: Lod=3; %He so dan nhiet cua lop cach nhiet, W/mK: lamda1=0.041; %He so dan nhiet cua ong cach nhiet, W/mK: lamda2=0.071; %He so dan nhiet cua ong dan moi chat, W/mK: lamda3=41.8; %Duong kinh ong lon cua binh bay hoi, m: Dbh=0.067; %Duong kinh ong lon cua binh bay hoi, m: dbh=0.028; %Duong kinh ong dan moi chat, m: dod=0.008; %Duong kinh ong khoi xe, m: dok=0.022; %Thoi gian can giu am, s: time=1800; %Tính tốn: Phụ lục Trang 103 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH % Truong hop 1: %1.Nhiet luong can de tri nhiet thuc pham, Q1: Q1=mtp*ctp*(tf1-ttp-5); %kJ %Nhiet khong trung binh tinh toan: tkk=0.5*(tf1+tf2); %Thông so cua khong o nhiet trung binh: A=[40 1.128 1.005 2.76e-02 2.43e-05 1.91e-05 16.96e-06 0.699 50 1.093 1.005 2.83e-02 2.57e-05 1.96e-05 17.95e-06 0.698 60 1.06 1.005 2.9e-02 2.72e-05 2.01e-05 18.97e-6 0.696]; rho_kk=interp1(A(:,1),A(:,2),tkk); %kg/m3 cp_kk=interp1(A(:,1),A(:,3),tkk); %kJ/kgK lamda_kk=interp1(A(:,1),A(:,4),tkk); %W/mK a_kk=interp1(A(:,1),A(:,5),tkk); %m2/s m_kk=interp1(A(:,1),A(:,6),tkk); %Pa.s u_kk=interp1(A(:,1),A(:,7),tkk); %m2/s Pr_kk=interp1(A(:,1),A(:,8),tkk); %The tich thung,Vth: Vth=Ltr*Wtr*Htr; %m3 %Khoi luong khong thung: mkk=rho_kk*Vth; %kg %2.Nhiet luong can de tri nhiet thung, Q2: Q2=mkk*cp_kk*(tf1-tf2); %kJ %3 Ton that nhiet moi truong: %3.1 Ton that nhiet tu thung moi truong, Q31: alpha1=10; %W/m2K alpha2=15; %W/m2K %He so truyen nhiet qua vach thung: kv=1/(1/alpha1+Tcn1/lamda1+1/alpha2); %W/m2K %Tong dien tich be mat ngoai cua thung: Fn=2*((Ltr+2*Tcn1)*(Wtr+2* Tcn1)+ (Ltr+2*Tcn1)*(Htr+2* Tcn1)+ (Htr+2*Tcn1)*(Wtr+2* Tcn1)); %m2 %3.1 Ton that nhiet tu thung moi truong: Q31=kv/1000*Fn*(tf1-tf2)*time; %Thong so cua hoi nuoc bao hoa: th=tk1-10; %degC %Cac thong so cua nuoc bao hoa: B=[100 1.013 958.4 0.598 2256.8 4.22 2.135 0.683 0.02372 2.83E-04 1.18E-05 1.75 1.08 110 1.43 950.600 0.8263 2230 4.233 2.071 6.82E-01 2.62E-02 2.55E-04 1.26E-05 1.6 1.09 120 1.98 943.400 1.121 2203 4.25 2.12 6.83E-01 2.75E-02 2.32E-04 1.30E-05 1.47 1.09]; p_h=interp1(B(:,1),B(:,2),th); % bar rho_l=interp1(B(:,1),B(:,3),th); % kg/m3 rho_g=interp1(B(:,1),B(:,4),th); % kg/m3 r_g=interp1(B(:,1),B(:,5),th); % kJ/kg cp_l=interp1(B(:,1),B(:,6),th); % kJ/kgK cp_g=interp1(B(:,1),B(:,7),th); % kJ/kgK lamda_l=interp1(B(:,1),B(:,8),th); % W/mK lamda_g=interp1(B(:,1),B(:,9),th); % W/mK m_l=interp1(B(:,1),B(:,10),th); % Pa.s m_g=interp1(B(:,1),B(:,11),th); % Pa.s Pr_l=interp1(B(:,1),B(:,12),th); Pr_g=interp1(B(:,1),B(:,13),th); %Nhiet khoi trung binh tinh toan: tk=0.5*(tk1+tk2); %Thong so cua khoi: C=[0 1.295 1.042 2.28E-02 1.69E-05 1.58E-05 1.22E-05 0.72 100 0.950 1.068 3.13E-02 3.08E-05 2.04E-05 2.15E-05 0.69 200 0.748 1.097 0.0401 4.89E-05 2.45E-05 3.28E-05 0.67]; rho_k=interp1(C(:,1),C(:,2),tk); % kg/m3 cp_k=interp1(C(:,1),C(:,3),tk); % kJ/kgK Phụ lục Trang 104 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH lamda_k=interp1(C(:,1),C(:,4),tk); % W/mK a_k=interp1(C(:,1),C(:,5),tk); % m2/s m_k=interp1(C(:,1),C(:,6),tk); % Pa.s u_k=interp1(C(:,1),C(:,7),tk); % m2/s Pr_k=interp1(C(:,1),C(:,8),tk); %Luu luong the tich va khoi luong khoi: Vk=pi*dok^2*vk/4; %m3/s Gk=Vk*rho_k; % kg/s %Nang suat nhiet toa tu khoi: Qk=Gk*cp_k*(tk1-tk2)*1000; % W %Mat dong nhiet truyen tu khoi den moi chat: Lbh=0.15; % m qk=Qk/(pi*dbh*Lbh); %He so toa nhiet cua moi chat soi: alpha_e=3.45*(p_h)^0.5*qk^0.7; % W/m2K %He so truyen nhiet duong: kl=1/(1/(alpha_e*dod)+1/(2*lamda3)*log((dod+2*Tod)/dod)+1/(2*lamda2)*log((d od+2*Tod+2*Tcn2)/(dod+2*Tod))+1/(alpha2*(dod+2*Tod+2*Tcn2))); % W/m2K %Tong dien tich be mat ngoai cua ong dan moi chat: Fod=pi*(dod+2*Tcn2)*Lod; % m2 %Ton that nhiet tren duong ong dan moi chat: Q32=pi*Fod*kl*(th-tf2)*time/1000; %Tong ton that moi truong: Q3=Q31+Q32; %kJ %Phu tai nhiet can thiet: Q=Q1+Q2+Q3; %kJ %Luong moi chat nap can thiet: mn=Q/cp_l/(th-tf2)*1000; %kg %The tich nuoc can nap: Vn=mn*1; % ml %He so toa nhiet ngung: tw=th-10; alpha_ng=0.72*(lamda_l^3*r_g*1000*9.81*rho_l*(rho_l-rho_g)/(m_l*dod*(thtw)))^0.25; %W/m2K %He so truyen nhiet voi ong tron: kng1=1/(1/alpha1+Tod/lamda3+1/alpha_ng); % W/m2K nu=0.7; Qng=Q*1000/(nu*time); % W/m2 %Tong dien tich be mat truyen nhiet cua dan ngung: Fng1=Qng/(kng1*(th-tf1)); %m2 %Tong chieu dai duong ong tron: Lng1=Fng1/(pi*dod); % m %He so truyen nhiet voi ong co canh: hc_ng=0.005; % m sc_ng=0.003; % m kng2=kng1*(1+0.25*(hc_ng/sc_ng)^0.5); % W/m2K %Tong chieu dai duong ong dan ngung co canh: Lng2=Lng1*kng1/kng2; % m %Chon loai dan ngung lop co chieu dai 0,4m va khoang cach giua ong la 40mm Lng=400; %mm Wng=(0.04+hc_ng*2+2*dod)*1000; %mm %So ong cua dan: n_ng=round(Lng2/0.4+0.5); % ?ng %Chieu cao cua dan: Hng=round((0.5*((n_ng-1)*0.04+n_ng*dod)*1000)); % m % The tich ong nhiet: Vhp=1000000*(pi*(Dbh^2-dbh^2)*Lbh/4+pi/4*dod^2*(0.5*Lod+1/3*Lng2)); % ml %Ty le nap: TLN=Vn/Vhp*100; % % % Truong hop 2: %Phu tai nhiet can thiet: Q=Q1+Q2+Q3; %kJ Phụ lục Trang 105 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH %Luong moi chat nap can thiet: mn=Q/cp_l/(th-tf2)*1000; %kg %The tich nuoc can nap: Vn=mn*1; % ml %He so toa nhiet ngung: tw=th-10; alpha_ng=0.72*(lamda_l^3*r_g*1000*9.81*rho_l*(rho_l-rho_g)/(m_l*dod*(thtw)))^0.25; %W/m2K %He so truyen nhiet voi ong tron: kng1=1/(1/alpha1+Tod/lamda3+1/alpha_ng); % W/m2K nu=0.7; Qng=Q*1000/(nu*time); % W/m2 %Tong dien tich be mat truyen nhiet cua dan ngung: Fng1=Qng/(kng1*(th-tf1)); %m2 %Tong chieu dai duong ong tron: Lng1=Fng1/(pi*dod); % m %He so truyen nhiet voi ong co canh: hc_ng=0.005; % m sc_ng=0.003; % m kng2=kng1*(1+0.25*(hc_ng/sc_ng)^0.5); % W/m2K %Tong chieu dai duong ong dan ngung co canh: Lng2=Lng1*kng1/kng2; % m %Chon loai dan ngung lop co chieu dai 0,4m va khoang cach giua ong la 40mm Lng=400; %mm Wng=(0.04+hc_ng*2+2*dod)*1000; %mm %So ong cua dan: n_ng=round(Lng2/0.4+0.5); % ?ng %Chieu cao cua dan: Hng=round((0.5*((n_ng-1)*0.04+n_ng*dod)*1000)); % m % The tich ong nhiet: Vhp=1000000*(pi*(Dbh^2-dbh^2)*Lbh/4+pi/4*dod^2*(0.5*Lod+1/3*Lng2)); % ml %Ty le nap: TLN=Vn/Vhp*100; % % % Truong hop 2: %Phu tai nhiet can thiet: Q_2=Q3; %kJ %Luong moi chat nap can thiet: mn_2=Q_2/cp_l/(th-tf2)*1000; %kg %The tich nuoc can nap: Vn_2=mn_2*1; % ml Qng_2=Q_2*1000/(nu*time); % W/m2 %Tong dien tich be mat truyen nhiet cua dan ngung: Fng1_2=Qng_2/(kng1*(th-tf1)); %m2 %Tong chieu dai duong ong tron: Lng1_2=Fng1_2/(pi*dod); % m %Tong chieu dai duong ong dan ngung co canh: Lng2_2=Lng1_2*kng1/kng2; % m %So ong cua dan: n_ng2=round(Lng2_2/0.4+0.5); % ?ng %Chieu cao cua dan: Hng_2=round((0.5*((n_ng2-1)*0.04+n_ng2*dod)*1000)); % m % The tich ong nhiet: Vhp_2=1000000*(pi*(Dbh^2-dbh^2)*Lbh/4+pi/4*dod^2*(0.5*Lod+1/3*Lng2_2)); % ml %Ty le nap: TLN_2=Vn_2/Vhp_2*100; % % % Truong hop 3: %Phu tai nhiet can thiet: Q_3=1.1*(Q1+Q2+Q3); %kJ %Luong moi chat nap can thiet: Phụ lục Trang 106 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH mn_3=Q_3/cp_l/(th-tf2)*1000; %kg %The tich nuoc can nap: Vn_3=mn_3*1; % ml Qng_3=Q_3*1000/(nu*time); % W/m2 %Tong dien tich be mat truyen nhiet cua dan ngung: Fng1_3=Qng_3/(kng1*(th-tf1)); %m2 %Tong chieu dai duong ong tron: Lng1_3=Fng1_3/(pi*dod); % m %Tong chieu dai duong ong dan ngung co canh: Lng2_3=Lng1_3*kng1/kng2; % m %So ong cua dan: n_ng3=round(Lng2_3/0.4+0.5); % ?ng %Chieu cao cua dan: Hng_3=round((0.5*((n_ng3-1)*0.04+n_ng3*dod)*1000)); % m % The tich ong nhiet: Vhp_3=1000000*(pi*(Dbh^2-dbh^2)*Lbh/4+pi/4*dod^2*(0.5*Lod+1/3*Lng2_3)); % ml %Ty le nap: TLN_3=Vn_3/Vhp_3*100; % % % Hien thi: set(handles.edit1,'String',Q); set(handles.edit2,'String',Q_2); set(handles.edit3,'String',Q_3); set(handles.edit4,'String',Vn); set(handles.edit5,'String',Vn_2); set(handles.edit6,'String',Vn_3); set(handles.edit7,'String',TLN); set(handles.edit8,'String',TLN_2); set(handles.edit9,'String',TLN_3); set(handles.edit10,'String',Lng2); set(handles.edit11,'String',Lng2_2); set(handles.edit12,'String',Lng2_3); set(handles.edit13,'String',Hng); set(handles.edit14,'String',Hng_2); set(handles.edit15,'String',Hng_3); % Ve thi: X=[Q_2 Q Q_3]; Y=[TLN_2 TLN TLN_3]; plot(handles.axes2,X,Y,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('SO SANH TY LE NAP THEO CAC TRUONG HOP','Fontsize',14); xlabel('Phu tai nhiet, kJ','Fontsize',14); ylabel('Ty le nap, %','Fontsize',14); grid on Phụ lục 3: Chương trình vẽ đồ thị biểu diễn kết thí nghiệm PL3.1 Với tỷ lệ nạp 58%: P=[0 26 66 Tth=[36 44 Tbh=[36 84 87 Tng=[36 64.7 78 81 83 86]; 50 55 59 61 63]; 91 97.5 99 101]; 69.3 73.3 80 81.3 82.3]; plot(P,Tth,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet he thong theo ap suat bay hoi','Fontsize',14); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on Phụ lục Trang 107 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH plot(P,Tbh,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize', 10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tng,'.bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize' ,10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); legend('Nhiet thung','Nhiet Te', 'Nhiet Tc'); grid on hold on PL3.2 Với tỷ lệ nạp 68%: P=[0 35.5 54.5 63.3 67.1 74.7 75.2]; Tbh=[30.5 80.3 87 90.3 87.8 91.1 92]; Tng=[30.5 76.9 80.9 86.7 84.5 87.3 91.1]; Tth=[30.5 47.3 57.5 62.3 65 68 69.5]; plot(P,Tth,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet he thong theo ap suat bay hoi','Fontsize',14); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tbh,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize', 10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tng,'.bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize' ,10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); legend('Nhiet thung','Nhiet Te', 'Nhiet Tc'); grid on hold on PL3.3 Với tỷ lệ nạp 78%: P=[0 36 Tbh=[30.4 Tng=[30.4 Tth=[30.4 56 66 68 73 89]; 74 78.2 83.8 85.5 87.6 89.8]; 68.4 72.8 78.2 82.8 85.3 87.6]; 46.2 55.3 64 66.6 68.9 71.1]; plot(P,Tth,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet he thong theo ap suat bay hoi','Fontsize',14); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on Phụ lục Trang 108 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH plot(P,Tbh,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize', 10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tng,'.bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize' ,10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); legend('Nhiet thung','Nhiet Te', 'Nhiet Tc'); grid on hold on PL3.4 Với tỷ lệ nạp 88%: P=[0 39 Tbh=[29.7 Tng=[29.7 Tth=[29.7 76 82.4 90.4 93.6 97.2]; 54.9 83.7 93.5 92.6 96.6 102]; 41.7 63.2 75.7 83.1 87 93.1]; 37.8 43.9 51.1 53.7 55.1 57.4]; plot(P,Tth,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet he thong theo ap suat bay hoi','Fontsize',14); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tbh,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize', 10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tng,'.bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize' ,10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); legend('Nhiet thung','Nhiet Te', 'Nhiet Tc'); grid on hold on PL3.5 Với tỷ lệ nạp 98%: P=[0 25.3 88.4 98.8 112.7 119.1 125.4]; Tbh=[28.9 48.9 93.4 101.4 101.3 107.3 114.6]; Tng=[28.9 29.3 53.6 70.1 83.4 88.7 99.5]; Tth=[28.9 29.4 32.4 38.2 40.8 41.2 43.7]; plot(P,Tth,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet he thong theo ap suat bay hoi','Fontsize',14); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on Phụ lục Trang 109 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH plot(P,Tbh,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize', 10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); grid on hold on plot(P,Tng,'.bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize' ,10); xlabel('Ap suat tuyet doi, cmHg','Fontsize',14); ylabel('Nhiet do, ^oC','Fontsize',14); legend('Nhiet thung','Nhiet Te', 'Nhiet Tc'); grid on hold on Phụ lục 4: Chương trình vẽ đồ thị biểu diễn thay đổi nhiệt trở theo thời gian PL4.1 Với tỷ lệ nạp 58%: Time=[10 R=[0.2056 20 30 0.1885 40 50 0.1885 60]; 0.1864 0.1885 0.1992]; plot(Time,R,'ks','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 58%'); grid on hold on PL4.2 Với tỷ lệ nạp 68%: Time=[10 R=[0.0362 20 30 0.0650 40 50 0.0383 60]; 0.0351 0.0405 0.0096]; plot(Time,R,'bs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 68%'); grid on hold on PL4.3 Với tỷ lệ nạp 78%: Time=[10 R=[0.0596 20 30 0.0575 40 50 0.0596 60]; 0.0288 0.0245 0.0234]; plot(Time,R,'rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 78%'); grid on hold on Phụ lục Trang 110 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH PL4.4 Với tỷ lệ nạp 88%: Time=[10 R=[0.1406 20 30 0.2183 40 50 0.1896 60]; 0.1012 0.1023 0.0948]; plot(Time,R,'gs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 88%'); grid on hold on PL4.5 Với tỷ lệ nạp 98%: Time=[10 R=[0.2088 20 30 0.4239 40 50 0.3334 60]; 0.1907 0.1981 0.1608]; plot(Time,R,'ms','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 98%'); grid on hold on PL4.6 Đồ thị tổng hợp nhiệt trở Time=[10 20 30 R58=[0.2056 0.1885 40 50 0.1885 60]; 0.1864 0.1885 0.1992]; plot(Time,R58,'ks','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','c','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); grid on hold on R68=[0.0362 0.0650 0.0383 0.0351 0.0405 0.0096]; plot(Time,R68,'bo','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','c','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); grid on hold on R78=[0.0724 0.0128 0.0469 0.0107 0.0064 0.0288]; plot(Time,R78,'rd','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','c','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); Phụ lục Trang 111 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); grid on hold on R88=[0.1406 0.2183 0.1896 0.1012 0.1023 0.0948]; plot(Time,R88,'g^','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','c','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); grid on hold on R98=[0.2088 0.4239 0.3334 0.1907 0.1981 0.1608]; plot(Time,R98,'mp','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','c','MarkerSize', 10); title('Su thay doi nhiet tro theo thoi gian','Fontsize',14); xlabel('Thoi gian, Phut','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Ty le nap 58%','Ty le nap 68%','Ty le nap 78%','Ty le nap 88%','Ty le nap 98%'); grid on hold on Phụ lục 5: Chương trình vẽ đồ thị tổng hợp kết thí nghiệm PL5.1 Sự thay đổi áp suất TLN=[58 68 78 88 98]; P=[60 52.9 55.4 68.4 81.4]; Tbh=[85.1 79.9 72.9 79 85.1]; Tng=[69.6 76.8 70.5 67.6 64.8]; Tth=[52.6 57.2 57.5 47 36.4]; R=[0.1652 0.0321 0.0254 0.1210 0.2165]; plot(TLN,P,'rs','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi ap suat trung binh theo ty le nap','Fontsize',14); xlabel('Ty le nap, %','Fontsize',14); ylabel('Ap suat, cmHg','Fontsize',14); legend('Ap suat trung binh'); grid on hold on PL5.2 Sự thay đổi nhiệt độ Te TLN=[58 68 78 88 98]; P=[60 52.9 55.4 68.4 81.4]; Tbh=[85.1 79.9 75.6 79 85.1]; Tng=[69.6 76.8 72.2 67.6 64.8]; Tth=[52.6 57.2 57.5 47 36.4]; R=[0.1652 0.0321 0.0362 0.1210 0.2165]; plot(TLN,Tbh,'rs','LineWidth',2,'MarkerEdgeColor','b','MarkerFaceColor','g','MarkerSize', 10); Phụ lục Trang 112 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH title('Su thay doi nhiet Te trung binh theo ty le nap','Fontsize',14); xlabel('Ty le nap, %','Fontsize',14); ylabel('Nhiet Te, ^oC','Fontsize',14); legend('Nhiet Te trung binh'); grid on hold on PL5.3 Sự thay đổi nhiệt độ Tc TLN=[58 68 78 88 98]; P=[60 52.9 55.4 68.4 81.4]; Tbh=[85.1 79.9 72.9 79 85.1]; Tng=[69.6 76.8 70.5 67.6 64.8]; Tth=[52.6 57.2 57.5 47 36.4]; R=[0.1652 0.0321 0.0254 0.1210 0.2165]; plot(TLN,Tng,'bs','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet Tc trung binh theo ty le nap','Fontsize',14); xlabel('Ty le nap, %','Fontsize',14); ylabel('Nhiet Tc, ^oC','Fontsize',14); legend('Nhiet Tc trung binh'); grid on hold on PL5.4 Sự thay đổi nhiệt độ thùng TLN=[58 68 78 88 98]; P=[60 52.9 55.4 68.4 81.4]; Tbh=[85.1 79.9 72.9 79 85.1]; Tng=[69.6 76.8 70.5 67.6 64.8]; Tth=[52.6 57.2 57.5 47 36.4]; R=[0.1652 0.0321 0.0254 0.1210 0.2165]; plot(TLN,Tth,'ms','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet thung trung binh theo ty le nap','Fontsize',14); xlabel('Ty le nap, %','Fontsize',14); ylabel('Nhiet thung, ^oC','Fontsize',14); legend('Nhiet thung trung binh'); grid on hold on PL5.5 Sự thay đổi nhiệt trở TLN=[58 68 78 88 98]; P=[60 52.9 55.4 68.4 81.4]; Tbh=[85.1 79.9 75.6 79 85.1]; Tng=[69.6 76.8 72.2 67.6 64.8]; Tth=[52.6 57.2 57.5 47 36.4]; R=[0.1652 0.0321 0.0362 0.1210 0.2165]; plot(TLN,R,'ks','LineWidth',2,'MarkerEdgeColor','r','MarkerFaceColor','g','MarkerSize', 10); title('Su thay doi nhiet tro trung binh theo ty le nap','Fontsize',14); xlabel('Ty le nap, %','Fontsize',14); ylabel('Nhiet tro, ^oC/W','Fontsize',14); legend('Nhiet tro trung binh'); grid on hold on Phụ lục Trang 113 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH TÀI LIỆU THAM KHẢO Tiếng Việt: [1] Bùi Hải – Dương Đức Hồng – Hà Mạnh Thư, Thiết bị trao đổi nhiệt, NXB Bách Khoa Hà Nội, 2000 [2] Bùi Hải – Trần Văn Vang, Ống nhiệt ứng dụng ống nhiệt, NXB Bách Khoa Hà Nội, 2008 [3] Hà Mạnh Thư, Nghiên cứu đặc tính nhiệt tới hạn ống nhiệt trọng truờng có xẻ rãnh bên trong, Luận án tiến sĩ kỹ thuật, truờng Đại học Bách Khoa Hà Nội, 1993 [4] Hồng Đình Tín, Cơ sở Truyền Nhiệt, NXB ĐHQG TpHCM, 2003 [5] Hồng Đình Tín – Lê Chí Hiệp, Nhiệt động lực học kỹ thuật, NXB ĐHQG TpHCM [6] Hồng Đình Tín, Truyền nhiệt tính tốn thiết bị trao đổi nhiệt, NXB Khoa học kỹ thuật, 2001 [7] Hoàng An Quốc – Lê Chí Hiệp – Hồng Dương Hùng, "Nghiên cứu ảnh hưởng lượng nạp đến chế độ khởi động ống nhiệt," Khoa Học Công Nghệ Nhiệt , vol 74, pp - 11, 2007 [8] Hoàng An Quốc – Lê Chí Hiệp – Hồng Dương Hùng, Nghiên cứu thử nghiệm ống nhiệt mặt trời loại chân không điều kiện Việt Nam Khoa Học Công Nghệ Nhiệt , vol 74, pp 8, 2007 [9] Trần Văn Vang, Nghiên cứu tính chất nhiệt chưa tới hạn tới hạn ống nhiệt trọng trường có tách dịng kiểu giao pha, Luận án tiến sĩ kỹ thuật, trường Đại học Bách Khoa Ðà Nẵng, 2004 Tiếng Anh: [10] ASHRAE Handbook, Refrigeration, 2006 [11] Bui Hai, A Study of Characteristics of Gravity-assisted of the Two-phrase Closed Thermosyphon, Ph.D Dissertation, Svuss, Czechoslovakia, 1984 [12] David Reay – Peter Kew, Heat pipes, Theory, Design and Application, Butterworth - Heinemann Press, 2006 [13] Dongxing Gai, Jinguo Yang, Zhichun Liu, Wei Liu (2009), Operational characteristics of miniature loop heat pipe with flat evaporator, Springer-Verlag 2009 [14] G.F.Hewitt - G.L.Shires - T.R.Bott, Process of Heat Transfer, CRC Press, 1993 [15] G.P.Peterson, An Introduction to Heat Pipes - Modeling, Testing, and Application, A Wiley - Interscience Publication [16] H.Imura, K.Sasaguchi and H.Kozai, "Critical heat flux in a closed two-phase thermosyphon," Int.J.Heat Mass Transfer, vol 26, no 8, pp 1181-1983, 1983 [17] H.Nguyen – Chi and M.Groll, "Entrainment of flooding limit in a closed two phase thermosyphon," Heat Recovery System, vol 1, pp 275 - 286, 1981 [18] P.D Dunt, D.A Reay (1994), Heat pipes, Pergamon Press, Oxford [19] S Mostafa Ghiaasiaan, Two-Phase Flow, Boiling and Condensation, Cambridge University Press, 2008 Tài liệu tham khảo Trang 114 CBHD: GS.TS LÊ CHÍ HIỆP HVTH: NGUYỄN TIẾN CẢNH [20] Satish G Kandlikar – Masahiro Shoji – Vijay K Dhir, Handbook of phase change: boiling and condensation, Taylor & Francis publisher, 1999 [21] Yunus A Çengel, Heat transfer 2nd edition, McGraw-Hill, 2003 [22] Yunus A Çengel – Michael A Boles, Thermodynamics An Engineering Approach 5th edition, McGraw-Hill, 2004 Tài liệu tham khảo Trang 115 ... chọn đề tài: ? ?Nghiên cứu lý thuyết thực nghiệm ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm cho thùng giao thực phẩm? ?? Có thể, ngày khơng xa nữa, thấy xe máy giao hàng thực phẩm nóng, ấm gà rán,... loại ống nhiệt nói chung đặc biệt loại ống nhiệt sử dụng nước để gia nhiệt - Nghiên cứu ống nhiệt tận dụng nhiệt thải - Tính tốn, thiết kế, chế tạo ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm. .. Nghiên cứu lý thuyết thực nghiệm ống nhiệt tận dụng nhiệt thải xe máy để giữ ấm cho thùng giao thực phẩm bao gồm nội dung sau: - Nghiên cứu sở lý thuyết liên quan đến tính tốn, thiết kế ống nhiệt