Nghiên cứu quá trình truyền nhiệt truyền chất trong bình hấp thụ của máy lạnh hấp thụ NH3 h2o loại liên tục phù hợp với điều kiện việt nam

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN HIẾU NGHĨA NGHIÊN CỨU QUÁ TRÌNH TRUYỀN NHIỆT-TRUYỀN CHẤT TRONG BÌNH HẤP THỤ CỦA MÁY LẠNH HẤP THỤ NH3-H2O LOẠI LIÊN TỤC PHÙ HỢP VỚI ĐIỀU KIỆN VIỆT NAM LUẬN ÁN TIẾN SĨ KỸ THUẬT TP HỒ CHÍ MINH NĂM 2017 ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN HIẾU NGHĨA NGHIÊN CỨU QUÁ TRÌNH TRUYỀN NHIỆT-TRUYỀN CHẤT TRONG BÌNH HẤP THỤ CỦA MÁY LẠNH HẤP THỤ NH3-H2O LOẠI LIÊN TỤC PHÙ HỢP VỚI ĐIỀU KIỆN VIỆT NAM Chuyên ngành: Mã số chuyên ngành: Kỹ thuật nhiệt 62520115 Phản biện độc lập 1: PGS.TS Nguyễn Việt Dũng Phản biện độc lập 2: PGS.TS Võ Chí Chính Phản biện 1: PGS.TS Hồng Đình Tín Phản biện 2: PGS.TS Bùi Trung Thành Phản biện 3: PGS.TS Đặng Thành Trung NGƯỜI HƯỚNG DẪN KHOA HỌC GS.TS LÊ CHÍ HIỆP PGS.TS HỒNG AN QUỐC LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu thân tác giả Các kết nghiên cứu, kết luận luận án trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu (nếu có) thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận án Nguyễn Hiếu Nghĩa i TÓM TẮT LUẬN ÁN Trong máy lạnh hấp thụ, trình hấp thụ quan trọng Vì lý này, cấu trúc hấp thụ trình truyền nhiệt-truyền chất hấp thụ thu hút nhiều nhà khoa học giới tập trung nghiên cứu vào lĩnh vực Nhằm mục đích phát triển ứng dụng máy lạnh hấp thụ NH3-H2O Việt Nam, hướng nghiên cứu mấu chốt chọn lựa cấu trúc hấp thụ với mong muốn chế tạo Việt Nam mà không cần đầu tư công nghệ Để giải vấn đề này, trình truyền nhiệt-truyền chất xuất dung dịch lỏng hấp thụ dạng màng chảy nghiên cứu lý thuyết thực nghiệm luận án Nhằm mục đích thực thí nghiệm hấp thụ đề xuất, máy lạnh hấp thụ hoàn chỉnh chế tạo có cơng suất lạnh khoảng kW thiết kế để sản xuất nước đá Mục tiêu thí nghiệm đề cập để thiết lập mối quan hệ hệ số truyền nhiệt hệ số truyền chất với biến số ảnh hưởng đến trình truyền nhiệt-truyền chất trình hấp thụ Để làm điều này, mơ hình tốn có liện quan đến q trình hấp thụ phát triển Dựa mơ hình tốn này, nghiên cứu lý thuyết thực cho hấp thụ máy lạnh hấp thụ để giới hạn vùng làm việc mà thí nghiệm sau tập trung vào Các kết cuối thực từ thí nghiệm sau dùng để kiểm tra lại với mô xác minh từ mơ hình tốn đề cập phía Bên cạnh đó, hệ số truyền nhiệt hệ số truyền chất nghiên cứu cịn so sánh với nghiên cứu có liên quan cho thấy phù hợp với Một chuỗi thí nghiệm xác định nồng độ từ 29,5% tới 32,5% để tìm nồng độ dung dịch nạp phù hợp Các kết cho thấy nồng độ nạp phù hợp 31% tương ứng với tính chất sau: hệ số hiệu COP = 0,436 hệ số hiệu làm lạnh nước muối COPu = 0,262, nhiệt độ trung bình NH3 rời khỏi bình phát sinh t5 = 116,5 °C, Nồng độ trung bình dung dịch loãng dung dịch đặc 29,14% 34,11% ii Nghiên cứu ảnh hưởng giảm nhiệt độ bay môi chất lạnh, tăng nhiệt độ ngưng tụ môi chất tăng, nhiệt độ hấp thụ dung dịch khỏi bình hấp thụ đến hiệu suất máy lạnh hấp thụ thực Một chương trình mơ hoạt động máy lạnh hấp thụ thiết lập, kết hợp tính tốn lý thuyết đo đạc thực tế, để kiểm tra đánh giá lại liệu thực nghiệm nhằm đảm bảo máy lạnh hấp thụ hoạt động theo thiết kế ban đầu Các thông số đo đạc điểm nút đưa vào chương trình mơ để xác định hiệu suất máy lạnh hấp thụ COP hiệu suất làm lạnh nước muối COPu để chế độ vận hành phù hợp cho mục đích sản xuất nước đá Khi nhiệt độ dung dịch nước muối từ -10 oC tới -19 oC, COP theo lý thuyết so với thực nghiệm có sai số trung bình 1,2% Các kết liệu mô máy lạnh hấp thụ NH3-H2O cấp phạm vi làm nước đá so sánh với mô từ nhiều báo khoa học khác giới kết mô tương đồng Mối tương quan của nhiệt độ phát sinh tối ưu theo nhiệt độ bay môi chất lạnh bay hơi, nhiệt độ ngưng tụ mơi chất lạnh bình ngưng tụ nhiệt độ hấp thụ dung dịch khỏi hấp thụ thiết lập Mơ hình tốn phát triển cho bình hấp thụ kiểu màng chảy chùm ống trịn nằm ngang xuất phát từ mơ hình tốn phát triển cho phần tử thể tích ống kiểm tra Phần tử thể tích ống đơn giản hố thành mơ hình vật lý chiều theo nhiều nghiên cứu trước Ống giải nhiệt có đường kính 9,6 mm Mật độ phân phối dung dịch loãng thấp để có chế độ nhỏ giọt vào ống Các trường vận tốc, nhiệt độ, nồng độ bề dày lớp màng dung dịch thay đổi theo điều kiện đầu vào dung dịch loãng nhiệt độ nước giải nhiệt chảy ống thể cho phần tử thể tích ống xét Chương trình mơ số viết ngơn ngữ lập trình Matlab để giải hệ phương trình vi phân phần mơ tả tượng truyền nhiệt-truyền chất trình hấp thụ mơ hình vật lý chiều Miền khảo sát màng dung dịch chảy theo nửa chu vi ống Hệ số tỏa nhiệt lớp màng αiw, hệ số truyền nhiệt từ mặt tiếp xúc lỏng-hơi vào nước giải nhiệt chảy ống k theo trục ε (x) Các hệ số tăng ¼ đầu nửa ống giảm dần ¼ sau cho thấy tốc độ hấp thụ giảm hệ số truyền nhiệt giảm Hệ số truyền chất tăng vọt đến vị trí dịng iii dung dịch lỗng vừa tiếp xúc với ống Sau đó, giảm nhanh phẳng trước khỏi ống Độ xác mơ hình số bình hấp thụ dạng màng chảy thực nghiệm so sánh nhiệt độ nước giải nhiệt vào khỏi chùm ống; cơng suất nhiệt bình hấp thụ Sai số tính tốn hệ thống Qa_compute so với kết thực nghiệm Qa_meas 4,3% Sai số kết mô (từ mối quan hệ hệ số truyền nhiệt k = f(C; Г; T) = f(0,308; 0,008; 306,3) = 0,863 kW/(m2.k) Qa_sim so với kết thực nghiệm Qa_meas 12,3% Ngoài ra, hệ số truyền nhiệt hệ số truyền chất nghiên cứu so sánh với nghiên cứu trước Hệ số truyền chất tìm từ phương trình mối quan hệ hệ số truyền chất hm = f(C; Г; T) = f(0,308; 0,008; 306,3) = 1,45*10-5 m/s Sangsoo Lee cộng tìm hệ số truyền nhiệt k = f(C; Г; P) = f(0,25; 0,008; 2,5) = 0,88 kW/(m2.K) hệ số truyền chất hm = f(C; Г; P) = f(0,25; 0,008; 2,5) = 1,65*10-5 m/s Mối quan hệ hệ số truyền nhiệt k [W/(m2.K)] hệ số truyền chất hm (m/s) trình hấp thụ với: (i) nồng độ dung dịch Cal khoảng từ 28% đến 31%, (ii) mật độ phân phối dung dịch theo chiều dài Г khoảng từ 0,001 [kg/(m.s)] đến 0,03 [kg/(m.s)] (iii) nhiệt độ nước giải nhiệt T khoảng từ 301 K đến 311 K thiết lập Mọi ý kiến đóng góp trao đổi, xin vui lịng liên hệ: Bộ mơn Cơng nghệ Nhiệt-Lạnh, khoa Cơ khí, Trường Đại học Bách Khoa Phịng 115B3, 268 Lý Thường Kiệt, Q.10, TP.HCM ĐT: (08) 3864-7256 Ext 5897 Email: llnghiaa@gmail.com iv ABSTRACT The absorption process has been confirmed as the most important process in absorption refrigeration machines in terms of improving their total efficiency For this reason, absorber structures in general and heat and mass transfers in absorber in particular have attracted the interest of many researchers in this field Aiming to develop the application of NH3-H2O absorption refrigeration machines in Vietnam, one of the key research directions is the selection of absorber structure which is expected to be fabricated in Vietnam without demand of new infrastructure investment In order to solve this problem, the heat and mass transfers occurring between liquid and vapor in falling film absorber have been studied theoretically and experimentally in this thesis For the purpose of performing the experiments on the proposed absorber, a complete NH3-H2O absorption refrigeration machine was fabricated, its refrigeration capacity is about 2kW and it was designed to produce ice The main target of the above-mentioned experiments is to establish the relationships of the heat transfer coefficient and mass transfer coefficient with the basic factors affecting the heat & mass transfer process of the absorption process In order to reach this target, the mathematical models related to absorption process were developed Based on these models, the theoretical studies were done for both absorber and absorption refrigeration machine in order to narrow the working area where the experiments later focused on The final results done by experiments were then verified with the simulation realized by the above-mentioned mathematical models and they were verified also with the ones published in the related literatures, it was proved that there is good conformity with them A series of experiments varying concentration from 29.5% to 32.5% was done to look for the suitable intake concentration The results showed that the suitable intake concentration is 31% which corresponds to the following properties: the coefficient of performance COP = 0.436, the brine refrigeration performance COPu = 0.262, the v average temperature of ammonia vapor leaving the generator is t5 = 116.5°C, the average concentrations of weak and strong solution are 29.14% and 34.11%, respectively Studies on the influence of the decrease of the evaporation temperature, the increase of the condensation temperature, the increase of the temperature of strong solution leaving the absorber to the performance of the absorption refrigeration machine were done as well The absorption refrigeration machine simulation program is established, is a combination of theoretical calculations and actual measurements The absorption refrigeration machine simulation confirmed as consistent with actual model in terms of design and operation The measured parameters at the state points are included in the simulation program to determine the coefficient of performance COP and brine refrigeration performance COPu to indicate the most suitable operating mode When the brine temperature is from -10 oC to -19 oC, the COP mean deviation of the computing program and practical result is 1.2% The computing program results of the single NH3-H2O absorption refrigeration machine to produce ice have been compared with other simulations from many scientific papers in the world are similar The relationship between the optimal generation temperature according to the saturation temperature of the refrigerant in the evaporator, the saturation temperature of the refrigerant in the condenser and the temperature of the strong solution leaving the absorber was set into an equation using polynomial regression method Mathematical model is developed for the falling film flowing on horizontal round tubes absorber derived from the mathematical model of the test volume element The test volume element is simplified into two-dimensional physical model in many previous studies Cooling tube diameter is 9.6 mm The liquid mass flow rate per unit length of tube is low to get droplet mode The fields of velocity, temperature, concentration and thickness of the falling film solution are varied by the input conditions of weak solution and cooling water temperature flowing in the tube represented by a test volume element of the tube The numerical simulation is written by the Matlab programming language to solve partial differential equations describing heat-mass transfer phenomenon of the absorption vi process of the 2-dimensional physical model Heat transfer coefficient of the film αiw, and the heat transfer coefficient from the liquid-vapor interface to the cooling water flowing in the k tube along ε (x) axis These coefficients increase for the first quarter of the tube and gradually decrease for the rest quarter showed the absorption rate decrease as the heat transfer coefficients decrease The mass transfer coefficient increases rapidly whereas the solution just contact the tube wall Then, it decreases sharply and quite smooth before leaving the tube wall The accuracy of numerical model and experiments are compared by the inlet, outlet the tube bundle of cooling water temperatures and absorber heat load The absorber heat load deviation of the computing program Qa_cumpute and experimental result Qa_meas is 4.3% The absorber heat load deviation of simulation result Qa_sim and experimental result Qa_meas is 12.3% The overall heat transfer coefficient k used for simulation result of absorber heat load was taken from the relationship of the heat transfer coefficient k = f(C; Г; T) = f(0.308; 0.008; 306.3) = 0.863 kW/(m2.K) In addition, the heat transfer and mass transfer coefficients of this research are compared with previous studies The mass transfer coefficient function hm = f(C; Г; T) = f(0.308; 0.008; 306.3) = 1.45*10-5 m/s Sangsoo Lee et al found heat transfer coefficient k = f(C; Г; P) = f(0.25; 0.008; 2.5) = 0.88 kW/(m2.K) and mass transfer coefficient hm = f(C; Г; P) = f(0.25; 0.008; 2.5) = 1.65*10-5 m/s The relationships of the heat transfer and the mass transfer coefficients of the absorption process to: (i) the solution concentration ranging from 28% to 31%, (ii) solution mass flow rate per unit tube length ranging from 0.001 kg/(m.s) to 0.03 kg/(m.s) and (iii) the cooling water temperature ranging from 301 K to 311 K were established If you have any comments in this PhD thesis, please feel free to contact: Department of Heat and Refrigeration, Faculty of Mechanic Engineering, Ho Chi Minh City University of Technology Room 115B3, 268 Ly Thuong Kiet St, Dist 10, Ho Chi Minh City Phone: (08) 3864-7256 Ext 5897 Email: llnghiaa@gmail.com vii LỜI CÁM ƠN Tôi xin bày tỏ lòng biết ơn chân thành sâu sắc tới GS.TS Lê Chí Hiệp, người trực tiếp hướng dẫn thực luận án tiến sĩ ý kiến mà Thầy cung cấp định hướng, nhận xét, góp ý vơ q báu Thầy suốt thời gian thực đề tài Tôi xin trân trọng cám ơn PGS.TS Hồng An Quốc tận tình hướng dẫn, giúp đỡ, nhận xét, góp ý tạo điều kiện thuận lợi cho suốt thời gian thực luận án Tôi xin chân thành cám ơn tới Thầy môn Công nghệ Nhiệt – Lạnh Thầy/Cô trường Đại Học Bách Khoa Thành Phố Hồ Chí Minh kiến thức lịng nhiệt tình với khoa học mà Thầy/Cơ truyền đạt suốt thời gian học tập nghiên cứu vừa qua Tôi xin chân thành cám ơn lãnh đạo khoa Công nghệ Nhiệt – Lạnh, trường đại học Cơng nghiệp thành phố Hồ Chí Minh ban Giám hiệu nhà trường tạo điều kiện thuận lợi hỗ trợ kinh phí học tập thời gian tơi làm nghiên cứu sinh Tôi xin gửi lời cám ơn đến anh em nghiên cứu sinh, bạn học viên cao học ngành Công nghệ Nhiệt, bạn Huỳnh Thế Vinh, em sinh viên Đại học Công Nghiệp TP HCM giúp đỡ, hỗ trợ thực mơ hình thí nghiệm Thành cơng luận án cịn thể lòng biết ơn sâu sắc đến người thân thương hỗ trợ động viên tác giả suốt thời gian thực luận án Tác giả Nguyễn Hiếu Nghĩa viii v22_20lim= sum(DataNew(ind20lim,24))/ (length(ind2010)+length(ind100)+length(ind0m10) +length(indm10m17p5) ) % %{ v8_20lim= sum(DataNew(ind20lim,20))/ (length(ind2010)+length(ind100)+length(ind0m10) ) % Tinh trung binh v13_20lim= sum(DataNew(ind20lim,21))/ (length(ind2010)+length(ind100)+length(ind0m10) ) v2_20lim= sum(DataNew(ind20lim,19))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t3_20lim= sum(DataNew(ind20lim,3))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t7_20lim= sum(DataNew(ind20lim,7))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t11_20lim= sum(DataNew(ind20lim,9))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t2_20lim= sum(DataNew(ind20lim,2))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t20_20lim= sum(DataNew(ind20lim,12))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t21_20lim= sum(DataNew(ind20lim,13))/ (length(ind2010)+length(ind100)+length(ind0m10) ) t22_20lim= sum(DataNew(ind20lim,14))/ (length(ind2010)+length(ind100)+length(ind0m10) ) v20_20lim= sum(DataNew(ind20lim,22))/ (length(ind2010)+length(ind100)+length(ind0m10) ) v21_20lim= sum(DataNew(ind20lim,23))/ (length(ind2010)+length(ind100)+length(ind0m10) ) v22_20lim= sum(DataNew(ind20lim,24))/ (length(ind2010)+length(ind100)+length(ind0m10) ) %} time = 1:length(VCOP); % Thoi gian chay, phut yyaxis left; plot(time, Data(1:r,11),'LineWidth', 2); % Ve truc t_Nacl ben trai yyaxis right plot(time,VCOP,'LineWidth', 2); % Ve truc COP ben phai %{ title('Plots with Temperature, COP') % Ten cua hinh xlabel('time, minutes (23/07/2016)') % Ten cua truc ngang ylabel('COP') % Ten cua truc dung yyaxis left ylabel('Temperature, \circC') % Ten cua truc dung ben trai %} %title('COP va COP_u thay doi theo cac nhiet do') % Ten cua hinh xlabel('Thoi gian, phut (14/07/2016)') % Ten cua truc ngang ylabel('COP') % Ten cua truc dung yyaxis left ylabel('Nhiet do, \circC') % Ten cua truc dung ben trai hold on yyaxis left plot(time, Data(1:r,10),'LineWidth', 2); % Ve truc t12 ben trai yyaxis right scatter(Vtime_u, VCOP_u,'LineWidth', 2) % Ve truc COP ben phai plot(Vtime_u, VCOP_u,'LineWidth', 2) % Ve truc COP ben phai yyaxis left plot(time, Data(1:r,5),'LineWidth', 2); % Ve truc t5 ben trai yyaxis left plot(time, Data(1:r,1),'LineWidth', 2); % Ve truc t1 ben trai %legend('t_N_a_C_l, C', 't_1_2, C','COP', 'COP_u'); legend('t_N_a_C_l, \circC', 't_1_2, \circC', 't_5, \circC', 't_1, \circC', 'COP', 'COP_u'); % X = [Data(:,11) VCOP'] xlswrite([date '.xlsx'],[[1:length(ind20lim)]' DataNew(ind20lim,[5, 10, 11, 20, 21, 19]) VCOP(ind20lim)' VXw(ind20lim)' VXs(ind20lim)' Vmws_v(ind20lim)' Vm1_v(ind20lim)' Vmss_v(ind20lim)' Vmws_meas(ind20lim)' Vm1_meas(ind20lim)' Vmss_meas(ind20lim)' Vmws(ind20lim)' Vm1(ind20lim)' Vmss(ind20lim)' ],1) % Luu du lieu PL28 xlswrite([date '.xlsx'],[Vtime_u(ind20lim_ireduced)' Vt5(ind20lim_ireduced)' Vt_NaCl(ind20lim_ireduced)' VCOP_u(ind20lim_ireduced)'],2) xlswrite([date '.xlsx'],[COP20' COP2010' COP100' COP0m10' COPm10m17p5' COPm17p5' COP_aver'; COP_u20' COP_u2010' COP_u100' COP_u0m10' COP_um10m17p5' COP_um17p5' COP_u_aver'; Xw20' Xw2010' Xw100' Xw0m10' Xwm10m17p5' Xwm17p5' Xw_aver'; Xs20' Xs2010' Xs100' Xs0m10' Xsm10m17p5' Xsm17p5' Xs_aver'],3) xlswrite([date '.xlsx'],[mws_v20' mws_v2010' mws_v100' mws_v0m10' mws_vm10m17p5' mws_vm17p5' mws_v_aver'; m1_v20' m1_v2010' m1_v100' m1_v0m10' m1_vm10m17p5' m1_vm17p5' m1_v_aver'; mss_v20' mss_v2010' mss_v100' mss_v0m10' mss_vm10m17p5' mss_vm17p5' mss_v_aver'; mws_vmeas20' mws_vmeas2010' mws_vmeas100' mws_vmeas0m10' mws_vmeasm10m17p5' mws_vmeasm17p5' mws_vmeas_aver'; m1_vmeas20' m1_vmeas2010' m1_vmeas100' m1_vmeas0m10' m1_vmeasm10m17p5' m1_vmeasm17p5' m1_vmeas_aver'; mss_vmeas20' mss_vmeas2010' mss_vmeas100' mss_vmeas0m10' mss_vmeasm10m17p5' mss_vmeasm17p5' mss_vmeas_aver'], 4) %lit/phut xlswrite([date '.xlsx'],[mws20' mws2010' mws100' mws0m10' mwsm10m17p5' mwsm17p5' mws_aver'; m120' m12010' m1100' m10m10' m1m10m17p5' m1m17p5' m1_aver'; mss20' mss2010' mss100' mss0m10' mssm10m17p5' mssm17p5' mss_aver'; mws_meas20' mws_meas2010' mws_meas100' mws_meas0m10' mws_measm10m17p5' mws_measm17p5' mws_meas_aver'; m1_meas20' m1_meas2010' m1_meas100' m1_meas0m10' m1_measm10m17p5' m1_measm17p5' m1_meas_aver'; mss_meas20' mss_meas2010' mss_meas100' mss_meas0m10' mss_measm10m17p5' mss_measm17p5' mss_meas_aver'], 5) % g/s xlswrite([date '.xlsx'],[v8_20lim' v13_20lim' v2_20lim' t3_20lim' t7_20lim' t11_20lim' t2_20lim' t20_20lim' t21_20lim' t22_20lim' v20_20lim' v21_20lim' v22_20lim'], 6) %find(VCOP==max(VCOP)) PL29 Phụ lục 2.2 Mô phần tử thể tích ống clear all clc di= 8*(10^-3); % Duong kinh ong giai nhiet, m F=6*3.5*pi*(di/2)^2; % Tiet dien ong, m2 V= 10* (10^-3)/60 ; %Luu luong nuoc giai nhiet, LPM, m3/s omega= V/F; % Van toc dong chay, m/s upsilon= 0.805*10^-6; % Do nhot dong hoc cua nuoc, m2/s Re= omega*di/V; Pr= 5.42; kw= 61.8*10^-2; % He so dan nhiet cua nuoc, W/mK hw1= 4.36*kw/di; % He so truyen nhiet doi luu cua nuoc giai nhiet, W/m2K hw2= 0.023*Re^0.8*Pr^(1/3)*(kw/di); % He so truyen nhiet doi luu cua nuoc giai nhiet, W/m2K % Tinh be day lop mang delta, m Veta=[]; Vepsilon=[]; Vdeltac=[]; Vdelta=[]; Vvc= []; Vv= []; Vuc =[]; Vu= []; Vu1c =[]; Vu1= []; VX=[]; VT=[]; Vmp = []; Vq=[]; VXint = []; VTint= []; Vhib= []; Vhbw= []; Vhiw= []; VU= []; Vhm= []; Vdelta_al= []; Vu1_al= []; VX_al= []; VT_al= []; % Xac dinh ap suat bay hoi tE=-12 ; %Nhiet bo bay hoi, C TE= tE + 273.15 ; PE=APSUATBH('apsuat',TE); P2= PE + 0.2; % Absorber pressure, bar Xin= 0.295; % Nong khoi luong dd vao BHT, Xin1= 0.295; % Nong khoi luong dd vao BHT, Xa= 0.9; % Hoi ammonia gan tinh khiet, .most important tin= 43; % Nhiet dd loang vao BHT, C Tin= tin + 273.15 ; % Tin1= tin + 273.15 ; % vector = [Tin1 Xin1] ; F=APSUATBH ('specvolume',vector); vint= double(F); % The tich rieng cua dd, m3/kg PL30 roint=1/vint; % KLR dd, kg/m3 vector = [Tin1 Xin1] ; F=APSUATBH ('specrosoli',vector); rointl= double(F); % Khoi luong rieng cua dd long, kg/m3 rho= rointl; %} mu= 3.958* 10^-4; % Do nhot dong luc, 4.8 *10^-3; 801.5 *10^-6 kg/m.s Gamma= 0.008 ; % Luu luong khoi luong phan bo, 0.0542; 0.005 kg/m.s (5 cho) g= 9.81; % Gia toc truong, m/s2 ro= 0.005; % Ban kinh ngoai ong giai nhiet, m anpha= 6.7*10^-8; % He so khuech tan nhiet cua nuoc, 15*10^-8 m2/s D= 4.4*10^-9; % He so khuech tan khoi cua nuoc, 28.2*10^-6; 1.07*10^-9 m2/s Twall= 303.15; % Nhiet vach ong giai nhiet, K R = 30; C = 60; % 60 spDepsilon = linspace(0.0001,1,C);% Chia theo chieu dong chay C cot Depsilon = spDepsilon(2) - spDepsilon(1);% Buoc doc theo dong chay spDeta = linspace(0.0001,1,R);% Chia theo chieu cat dong chay R hang Deta = spDeta(2) - spDeta(1);% Buoc cat dong chay X = ones(R , 1)*(Xin1); T1 = ones(R , 1)*(Tin1); dXdeta0 = 0.0029606; % Chenh nong tu dong hoi toi lop lien ke mat tiep xuc, tang Xint tang for j = 1:length(spDepsilon) % Toa KTN theo dong chay, j= 1, 2, 3, upsilon= mu/rho; % Do nhot dong hoc, m2/s Dep = spDepsilon(j); % Gia tri tai diem j Qx = zeros(R + 2,R); Qt = zeros(R + 2,R); if j~=1 X= VX(:,j -1); % Vector X cua lan j truoc giu co dinh cho lan j tiep theo T1=VT(:,j-1); end minEr = 10; % Sai so nho nhat gan ban dau cua moi lan chay j %for mp = 0.01: 0.001: 0.04 % assume absorbed mass flux, kg/(m2.s) mp= 0.0025; for i= 1:length(spDeta) % Toa KTN cat dong chay Dt = spDeta(i); % Gia tri tai diem i delta= ( (3*upsilon*Gamma/(rho^2 *g*sind(180*Dep)))^(1/3)); %be day lop mang delta, m u= ( (rho* g* delta^2* sind(180*Dep)/mu)*(Dt- 0.5* Dt^2))/10; %van toc doc theo dong chay u, m/s ddeltadepsilon= ( -(mu* Gamma* pi^3/(9* rho^2* g))^(1/3) * (1/(sind(180*Dep)^(1/3))* tand(180*Dep))); % Vi phan toan phan cua ti so delta/epsilon v= ( -(rho* g* delta^2* Dt^2/ (2* mu* ro))* ((1/pi)*ddeltadepsilon* pi*cosd(180*Dep) + delta*(1- (Dt/3)* cos(180*Dep)))); Vdeltac = [Vdeltac; delta]; Vuc = [Vuc; u]; % Luu u theo i cot Vu1c = [Vu1c; u1]; % Luu u theo i cot Vvc = [Vvc; v]; if delta > 1.5736e-05 % Phu thuoc luu chat Xint= Xin1 ; else PL31 Xint= (rho*D/(delta*mp ))* dXdeta0; % Du doan nong tai mat tiep xuc, *1.03 % end %} Rx= X/Depsilon; Rx(1,1) = Xin1; % Tai vach ong thay doi nong on dinh Rx(length(spDeta),1) = Xint; % Tai mat tiep xuc thay doi nong theo ham cua mass flux b_o= (1/Depsilon + 2*pi*ro*D/(u*delta^2*Deta^2)); c_o= (1/(delta*Deta))* ((pi*ro/u)* (v/2- D/(delta*Deta)) - (0.5* Dt* ddeltadepsilon)); a_o= (1/(delta*Deta))* (-(pi*ro/u)* (v/2+ D/(delta*Deta)) + (0.5* Dt* ddeltadepsilon)); Qx(i,i) = a_o; Qx(i + 1,i) = b_o; Qx(i + 2,i) = c_o; Qx(1,1) = 0; Qx(2,1) = 1; Qx(3,1) = 0; Qx(length(spDeta),length(spDeta)) = 0; Qx(length(spDeta)+ 1,length(spDeta)) = 1; Qx(length(spDeta)+ 2,length(spDeta)) = 0; % Nhiet tai mat tiep xuc duoc xac dinh theo DKCB ap suat hoi tu nong % mat tiep xuc va ap suat bo hap thu, C vector = [P2 Xint] ; F=APSUATBH ('sattemp',vector); Tintbh=double(F.d(1)); % Nhiet bao hoa tai trang thai vua dat ap suat P13, X13, K tintbh=Tintbh-273.15; %C Tint= Tintbh; % Nhiet tai mat tiep xuc, %C vector = [Tintbh P2 Xint]; F=APSUATBH ('enlisolic',vector); hintl= double(F); % Enthalpy NH3 bao hoa am, kJ/kg hint= hintl; % Enthalpy tai mat tiep xuc, kJ/kg vector = [Tintbh P2 Xint]; % Them hoi bao hoa F=APSUATBH ('envaammo',vector); %hintv= double(F); % Enthalpy, kJ/kg %} % Xac dinh NH3 enthalpy bao hoa kho Dqn= 40; % Do qua nhiet cua hoi NH3 sau BBH, C TE=273.15+tE + Dqn; % K, hintv = APSUATBH('enthalpy_v',TE); Rt= T1/Depsilon; Rt(1,1) = Twall; % Tai vach ong nhiet cua mang bang nhiet vach Rt(length(spDeta),1) = Tint; % Tai mat tiep xuc nhiet bang nhiet bao hoa theo ap suat va nong tai MTX b_t= (1/Depsilon + 2*pi*ro*anpha/(u*delta^2*Deta^2)); c_t= (1/(delta*Deta))* ((pi*ro/u)* (v/2- anpha/(delta*Deta)) - (0.5* Dt* ddeltadepsilon)); a_t= (1/(delta*Deta))* (-(pi*ro/u)* (v/2+ anpha/(delta*Deta)) + (0.5* Dt* ddeltadepsilon)); Qt(i,i) = a_t; Qt(i + 1,i) = b_t; Qt(i + 2,i) = c_t; Qt(1,1) = 0; Qt(2,1) = 1; Qt(3,1) = 0; PL32 Qt(length(spDeta),length(spDeta)) = 0; Qt(length(spDeta)+ 1,length(spDeta)) = 1; Qt(length(spDeta)+ 2,length(spDeta)) = 0; Qtempx = Qx(2:end-1, 1:end)'; Ax = diag(Qtempx,-1); % Lay duong cheo duoi Bx = diag(Qtempx); % Lay duong cheo giua Cx = diag(Qtempx,1); % Lay duong cheo tren X1 = thomas(Ax, Rx, Cx, Bx); %exQt = Qt(2:end-1, 1:end); Qtempt = Qt(2:end-1, 1:end)'; At = diag(Qtempt,-1); % Lay duong cheo duoi Bt = diag(Qtempt); % Lay duong cheo giua Ct = diag(Qtempt,1); % Lay duong cheo tren T = thomas(At, Rt, Ct, Bt) ; % Kiem tra dieu kien hoi tu tinh lai mass flux mpp kf= 0.384 ; % He so dan nhiet cua lop mang; W/m.K mpMin= 0.000611; Ta= 273.15 + tE + Dqn; % K dTdr0= 0.323; % Chenh nhiet tu cuc hoi toi mat tiep xuc, C dTdri= 0.32; % Chenh nhiet tu vach ong toi cuc hoi, C q= kf* dTdr0/ delta; % Heat flux, W/m2 ha= 1461.86; % Enthalpy cua NH3 tai P2=2.9bar, t2= 30C, kJ/kg hw= (83.7 + 167.5)/ ; % Enthalpy cua nuoc tai P2=3bar, t2= 30C, kJ/kg %hab= hintv - hint + Xint* (Xint*ha + (1- Xint)* hw); % Nhiet hap thu, kJ/kg hs= 465 ; %kJ/kg hv= 655; %kJ/kg k= 0.134; %W/mK %Xint=0.85; ha= 1461.86; % Enthalpy cua NH3 tai P2=2.9bar, t2= 30C, kJ/kg hw= (83.7 + 167.5)/ ; % Enthalpy cua nuoc tai P2=3bar, t2= 30C, kJ/kg %hab= (hv - hs + Xint* (Xint*ha + (1- Xint)* hw)); % Nhiet hap thu, kJ/kg mpp= (q)/(hab* 1000); % Tinh lai mass flux, kg/(m2.s) tang mpMin giam; nho tang nghieng Residue= abs((mpp-mp)/mp); % Sai so tuong doi X2=X1; T2=T; Xint2= Xint ; Tint2= Tint; Xsb= X2(length(spDeta)-R/2); Tsb= T2(length(spDeta)-R/2); if C/2< j && X2(length(spDeta)) < VX(end, end) % Hieu chinh phu hop; || (Xint2 - VX(end, end))> 0.05 Xint= VX(end, end) + 0.0005-delta; [X2, T2, Xint2, Tint2, Xsb, Tsb] = HIEUCHINH1(X,Xin1,spDeta,Depsilon,ro,D,delta,Deta,u,v,Dt,ddeltadepsilon,P2,Xint,T1,Twall,anpha,i,tE,Qx, Qt,R) end if delta > 1.5736e-05 % Phu thuoc luu chat 1*10^-3 hib= 0; PL33 hbw= 0; hm= ; hiw= 0; U= 0; else hib= kf* dTdr0/ ((Tint2-Tsb)* delta); % He so truyen nhiet (doi luu) tu mat tiep xuc toi cuc hoi, W/m2K hbw= kf* dTdri/ ((Tsb-Twall)* delta); % He so truyen nhiet (doi luu) tu cuc hoi toi vach, W/m2K delta_wall= 0.0008; % Do day ong Inox, m k_wall= 15.4; % He so dan nhiet cua Inox, W/mK hwc= 2068.756*(omega^0.8/di^0.2); % He so truyen nhiet (doi luu) cua nuoc giai nhiet, W/m2K hiw= (1/hib + 1/hbw)^-1; %He so truyen nhiet tong cua lop mang, W/m2K U= (1/hiw + delta_wall/k_wall + 1/hwc)^-1; %He so truyen nhiet tong tu mtx toi nuoc giai nhiet, W/m2K hm= mpMin/ (rho*(Xint2-Xsb)); % He so truyen khoi, m/s % end if hm > 6e-05 hm= 6e-05; end %Vq= [Vq qMin] %Vmp = [Vmp mpMin] % Sau moi j giu lai mp co Residue Vdelta = [Vdelta Vdeltac]; Vdeltac= []; Vu = [Vu Vuc]; Vuc = []; % Cot gan lai rong de khong bi xep noi tiep chay j tiep theo Vu1 = [Vu1 Vu1c] Vu1c = []; % Cot gan lai rong de khong bi xep noi tiep chay j tiep theo Vv = [Vv Vvc]; Vvc = []; VX = [VX X2] % Sau moi j giu lai X1 tuong ung Residue VT = [VT T2] VXint= [VXint Xint2] VTint= [VTint Tint2] Vhib= [Vhib hib] Vhbw= [Vhbw hbw] Vhiw= [Vhiw hiw] VU= [VU U] Vhm= [Vhm hm] delta_al= sum(Vdelta,1)/length(Vdelta(:,1)) % Xuat gia tri trung binh cuc bo theo chieu cat dong chay u1_al= sum(Vu1,1)/length(Vu1(:,1)) % Cong ma tran theo cot roi chia cho so hang cua cot mot X_al=sum(VX,1)/length(VX(:,1)) % Xuat gia tri trung binh cuc bo theo chieu cat dong chay T_al=sum(VT,1)/length(VT(:,1)) % Xuat gia tri trung binh cuc bo theo chieu cat dong chay end delta_aver= sum(delta_al(1:59))/59 % && isinf(delta_al) % Xuat gia tri trung binh cua phan bo u1_aver= sum(u1_al(1:59))/59 % && isinf(u1_al) % Xuat gia tri trung binh cua phan bo X_aver= sum(X_al(1:59))/59 % && isnan(X_al) % Xuat gia tri trung binh cua phan bo T_aver= sum(T_al(1:59))/59 % && isnan(T_al) % Xuat gia tri trung binh cua phan bo hib_aver= sum(Vhib(Vhib ~= 0))/length(Vhib(Vhib ~= 0)) % Xuat gia tri trung binh cua phan bo hbw_aver= sum(Vhbw(Vhbw ~= 0))/length(Vhbw(Vhbw ~= 0)) % Xuat gia tri trung binh cua phan bo hiw_aver= sum(Vhiw(Vhiw ~= 0))/length(Vhiw(Vhiw ~= 0)) % Xuat gia tri trung binh cua phan bo U_aver= sum(VU(VU ~= 0))/length(VU(VU ~= 0)) % Xuat gia tri trung binh cua phan bo hm_aver= sum(Vhm(Vhm ~= 0))/length(Vhm(Vhm ~= 0)) % Xuat gia tri trung binh cua phan bo % hold on [X_axis,Y_axis] = meshgrid(spDepsilon, spDeta); [r,c]=size(X_axis); PL34 %} figure rVu = reshape(Vu,r,c); mesh(X_axis, Y_axis, rVu); xlabel('Toa KTN nua chu vi ong, \epsilon') ylabel('Toa KTN be day mang long, \eta') zlabel('Van toc, u (m/s)') figure rVv = reshape(Vv,r,c); mesh(X_axis, Y_axis, rVv); xlabel('Toa KTN nua chu vi ong, \epsilon') ylabel('Toa KTN be day mang long, \eta') zlabel('Van toc, v (m/s)') figure rVdelta = reshape(Vdelta,r,c); mesh(X_axis, Y_axis, rVdelta); xlabel('Non-dimensional tube half-circumference, \epsilon') ylabel('Non-dimensional film thickness, \eta') zlabel('Film thickness, \delta (m)') figure rVX = reshape(VX,r,c); mesh(X_axis, Y_axis, rVX); xlabel('Toa KTN nua chu vi ong, \epsilon') ylabel('Toa KTN be day mang long, \eta') zlabel('Truong nong hap thu, \omega') figure rVT = reshape(VT,r,c); mesh(X_axis, Y_axis, rVT); xlabel('Toa KTN nua chu vi ong, \epsilon') ylabel('Toa KTN be day mang long, \eta') zlabel('Truong nhiet phan bo, T (K)') figure plot(spDepsilon, Vhiw, 'k*', 'MarkerSize', 6) hold on plot(spDepsilon, VU, 'ro', 'MarkerSize', 6) xlabel('Toa KTN nua chu vi ong, \epsilon') ylabel('He so truyen nhiet, [W/(m_2_._K)]') legend('\alpha_i_b', 'k') figure %plot(spDepsilon, Vhm, 'ro', 'MarkerSize', 5) plot(spDepsilon, Vhm) %xlabel('Toa KTN nua chu vi ong, \epsilon') %ylabel('He so truyen chat, h_m (m/s)') xlabel('Non-dimensional tube half-circumference, \epsilon') ylabel('Mass transfer coefficient, h_m (m/s)') PL35 Phụ lục 2.3 Mô chùm ống song song % Analysis of horizontal tube clear all clc global anpha D Gama Gama= 0.0051 ; % Luu luong khoi luong phan bo, 0.0542 kg/m.s Do= 9.6*10^-3; % Duong kinh ngoai cua ong giai nhiet, m Di= 8.8*10^-3; % Duong kinh ngoai cua ong giai nhiet, m kwall= 15.3; % He so dan nhiet cua ong gn, W/mK kw= 0.618; % He so dan nhiet cua nuoc, W/mK hw= 4.36*kw/ Di; % He so truyen nhiet doi luu cua nuoc ong, W/m2K DL= 2* 10^-2; % Chieu dai cua the tich kiem tra, m ATV= pi*Do*DL; % Dien tich kiem tra, m2 Ms= Gama* DL; % Luu luong khoi luong dung dich tren the chieu dai cua VT, kg/s Mw= 0.0257; % Luu luong khoi luong nuoc giai nhiet; 0.15, kg/s Cps= 1960; % Nhiet dung rieng cua dung dich; J/kgK Cpw= 4187; % Nhiet dung rieng cua nuoc giai nhiet; J/kgK for m=1:1:28 % So tang ong giai nhiet [Tws_ins, Cws_ins, Twis, Twos, hiw, U, hm, qf, mf] = BFarhaniehPro23Absorber(Twallo, hwc, Cpw, Tws_in, Cws_in, Twi, Two) kf= 0.134 ; % He so dan nhiet cua lop mang; W/m.K dTdr0= 3; % Chenh nhiet tu cuc hoi toi mat tiep xuc, C dTdri= 0.5; % Chenh nhiet tu vach ong toi cuc hoi, C q= kf* dTdr0; % Heat flux, W/m2 dQ= q* ATV; % Dong nhiet qua dien tich kiem tra, W ErTwn= 10; %Sai so thiet lap ban dau for Twkk= 303.2: 0.02 : 306.15 % Guess nhiet nuoc giai nhiet VT, K for Twalli= 303.4: 0.02 : 306.15 %Twk= Twkk - (2* dQ2/ (Mw*Cpw)) ; % Tinh Nhiet vao VT, K %dTw= Twkk- Twk; % Chenh nhiet va vao TV cua nuoc giai nhiet, K %dQ2= Mw* Cpw* dTw; % Cong suat nhiet nuoc giai nhiet mang khoi TV, W P= solve([num2str(Twkk) '- ( dQ3/ (' num2str(Mw) '*' num2str(Cpw) '))- Twk'],[num2str(Twkk) '- Twk - dTw' ], [ num2str(Mw) '*' num2str(Cpw) '* dTw- dQ3' ]); Twk= double(P.Twk) dTw= double(P.dTw) dQ3= double(P.dQ3) Tw= (Twkk+Twk)/2; % Nhiet trung binh nuoc giai nhiet TV, K %dQ3= hw* pi* Di* DL*(Twalli-Tw); % Cong suat nhiet truyen tu vach vao nuoc giai nhiet, W %dQ1= 2* pi* kwall* DL*(Twallo-Twalli)/ log(Do/Di) % Cong suat nhiet truyen qua vach vao nuoc giai nhiet, W Twalli_cal= dQ3/(hw*pi*Di*DL) + Tw; % Tu dQ3 end for n=1:1:10 Twn= Twkk; % Nhiet TV truoc la nhiet vao cua TV sau, K Twnnext= Twn+ 2*dQ/ (Mw*Cpw); % Nhiet cua TV dang xet Twme= (Twn+Twnnext)/ 2; % Nhiet trung binh cua nuoc giai nhiet cua TV, K dTwme= Twnnext- Twn; % Chenh nhiet va vao TV cua nuoc giai nhiet, K dQ2= Mw* Cpw* dTwme; % Cong suat nhiet nuoc giai nhiet mang khoi TV, W Twalli= dQ2/(hw*pi*Di*DL) + Twme; % Tu dQ3 dQ3= hw* pi* Di* DL*(Twalli-Twme); % Cong suat nhiet truyen tu vach vao nuoc giai nhiet, W PL36 Twn= Twnnext; % The Nhiet cua TV dang xet la nhiet vao TV ke tiep end Row= [1 12 16 19 22 25 28]; tw_out = 34.2 + 273; % Nhiet nuoc khoi bo HT, C, K tw_in = 31 + 273; % Nhiet nuoc vao bo HT, C, K del_tw = tw_out - tw_in; % Chenh nhiet nuoc vao va ra, C pass = 7; % So pass nuoc giai nhiet qua bo HT, theta_t = del_tw/pass; coef = pass/2; a = theta_t/ coef - 0.2*(theta_t/ coef); % So gia nhiet t = zeros(pass+2,1); % Cap phat mang cho vector t % formula t(end)= tw_out ; t(end-1) = (tw_out - (theta_t+coef*a)) ; for i = pass:-1:1 coef = coef - 1; t(i) = (t(i+1) - (theta_t+coef*a)) ; end if m== 1: t= t(4); elseif m == 5: t= t(8); elseif m == 9: 12 t= t(12); elseif m == 13: 16 t= t(16); elseif m == 17: 19 t= t(29); elseif m == 20: 22 t= t(22); elseif m == 23:25 t= t(25); elseif m == 26:28 t= t(28); end Twn4mea= t; % Nhiet khoi ong giai nhiet duoc 304, K ResTwn= abs(Twn-Twn4mea); % Sai so if Twn == Twn4mea; Twn_s = Twn; dQ2_s= dQ2; Tw_s= Tw; Twalli_s= Twalli; break; elseif ResTwn < ErTwn; % Sai so phai nho nhat ErTwn = ResTwn; Twn_s = Twn; dQ2_s= dQ2; Tw_s= Tw; Twalli_s= Twalli; end end del_t = t2r + 273 - tw; % Chenh nhiet nuoc vao va (tss - tw), C, K qf= U*del_t; % Heat flux, W/m2 PL37 rho_ss= 880; % Khoi luong rieng dd, kg/m3 Cint= 0.98; % Nong hoi NH3 tai bien mf= rho_ss*(Cint-X2).*hmz ; % Dong chat duoc HT, kg/(m2.s) do=0.0096; % Duong kinh ngoai ong giai nhiet, m l= 0.18; % Chieu dai doan ong giai nhiet, m n=6; % So ong cua moi tang ong, cai A1= pi*do*l*n; % Dien tich cua moi tang ong, m2 %Rows= [0 3*0.95 3*0.9 3*0.85 3*0.8 4*0.75 4*0.7 4*0.65 4*0.6]; Rows= [0 4*0.6 4*0.65 4*0.7 4*0.75 3*0.8 3*0.85 3*0.9 3*0.95]; A= cumsum(Rows'*A1) %Dien tich trao doi nhiet tich luy, m2 Qf= cumsum(qf.*Rows'*A1) *10^-3 * 0.6 %Cong suat bo hap thu, kW P1= 2; X2 = zeros(length(t),1); for i=length(Cs):-1:1 vector = [t(i) P1]; % Them F=APSUATBH ('satconcen',vector); Xs=double(F.c(1)); % Nong dung dich dac vao BPS, tinh duoc X2(i)=Xs; end end ErTwn dQ2_s PL38 Phụ lục Một số hình ảnh thí nghiệm Hình Sản xuất nước đá PL39 180 Hình Chế tạo máy sản xuất nước đá giai đoạn đầu PL40 PL41 Phụ lục Sai số thành phần đo nồng độ dung dịch Bảng kết tính độ khơng đảm bảo đo NH3 Bảng tổng hợp sai số dụng cụ Loại DCTT Đại lượng Bđm 250 250.0817 Buret 25 ml 24.9750 Pipet bầu 10 10.0028 Cân 200.6650 Sai số ống chuẩn H2SO4 0.1 Sai số qui trình 32.790 u(V) 0.1496 0.0948 0.0251 0.0856 0.0010 0.3868 Sai số Qui trình phân tích NH3 33.30 32.20 32.83 32.99 33.00 32.65 32.75 32.60 32.790 SS2.1 SS2.2 Trung bình ĐỘ KHÔNG ĐẢM BẢO ĐO NH3 0.529123909 (X ± 1,04)% 𝛼𝑤 = 2068,756 PL42 u/(V) 0.0006 0.0038 0.0025 0.0004 0.0100 0.0118 𝜔 0,8 𝑑𝑖0,2 0.38678 1.03708 ... NGHĨA NGHIÊN CỨU QUÁ TRÌNH TRUYỀN NHIỆT-TRUYỀN CHẤT TRONG BÌNH HẤP THỤ CỦA MÁY LẠNH HẤP THỤ NH3- H2O LOẠI LIÊN TỤC PHÙ HỢP VỚI ĐIỀU KIỆN VIỆT NAM Chuyên ngành: Mã số chuyên ngành: Kỹ thuật nhiệt. .. trình lạnh hấp thụ chủ yếu sử dụng máy lạnh công nghiệp lớn Đối với máy lạnh hấp thụ có cơng suất nhỏ cần thêm q trình nghiên cứu phận chu trình lạnh hấp thụ Điều giúp thiết kế máy lạnh hấp thụ. .. hình nghiên cứu máy lạnh hấp thụ NH 3H2O loại liên tục trình truyền nhiệt -truyền chất xuất dung dịch lỏng hấp thụ dạng màng chảy để đánh giá ưu điểm tồn nghiên cứu nước Một hướng nghiên cứu mấu