Xây dựng mô đun phần mềm tính toán các thông số của không khí ẩm và các ứng dụng

Bài viết trình bày phương pháp xây dựng mô-đun phần mềm xác định các thông số của không khí ẩm, ứng dụng để xây dựng đồ thị I-d và giải một số bài toán chuyên ngành, cụ thể liên quan đến bài toán dùng nước ngầm làm mát chuồng trại.

Thái Ngọc Sơn 64 XÂY DỰNG MƠ-ĐUN PHẦN MỀM TÍNH TỐN CÁC THƠNG SỐ CỦA KHƠNG KHÍ ẨM VÀ CÁC ỨNG DỤNG BUILDING SOFTWARE MODULE TO DETERMINE PARAMETERS OF MOIST AIR AND ITS APPLICATIONS Thái Ngọc Sơn Trường Đại học Bách khoa – Đại học Đà Nẵng; thaingocson@gmail.com Tóm tắt - Khơng khí ẩm mơi chất phổ biến ngành Kỹ thuật Nhiệt Để tính tốn thơng số khơng khí ẩm, ta sử dụng cơng thức giải tích sử dụng đồ thị I-d Hiện tại, việc sử dụng máy tính để giải tốn phổ biến; việc xây dựng mơ-đun phần mềm tính tốn thơng số vật lý, thông số trạng thái nhiệt động môi chất khác cần thiết, đặc biệt toán nhiệt thường xuyên phải sử dụng phương pháp tính lặp Bài báo trình bày phương pháp xây dựng mô-đun phần mềm xác định thông số khơng khí ẩm, ứng dụng để xây dựng đồ thị I-d giải số toán chuyên ngành, cụ thể liên quan đến toán dùng nước ngầm làm mát chuồng trại Abstract - Moist air is a widely used working fluid in thermal engineering To calculate the parameters of moist air, we can use analytical formulas or I-d graph Currently, the use of computers to solve technical problems is very common; building a software module that can calculate the physical parameters and thermodynamic state parameters of different refrigerants is very necessary, especially in thermal problems which normally require the use of iterative methods This article presents the method of building a software module to determine the parameters of moist air, its application to construct the Id graph, and solves some specific problems related to the problem of using ground water for cooling animal sheds Từ khóa - khơng khí ẩm; phương pháp lặp; nước ngầm; làm mát chuồng trại; mô-đun phần mềm; đồ thị I-d Key words - moist air; iterative methods; ground water; cooling animal sheds; software module; I-d graph Đặt vấn đề Ngày việc sử dụng máy tính để giải tốn kỹ thuật phổ biến Việc tính tốn ngành Kỹ thuật Nhiệt thường gắn liền với việc tra cứu tính chất nhiệt vật lý, thơng số trạng thái mơi chất Ngồi ra, kỹ thuật tính tốn thường hay sử dụng phương pháp lặp, khiến cho công việc thêm nặng nề Khơng khí ẩm mơi chất phổ biến ngành Kỹ thuật Nhiệt Việc dạy học, thực tế sản xuất đòi hỏi phải xây dựng phần mềm để tính tốn thơng số khơng khí ẩm Có nhiều phần mềm để tính tốn thơng số khơng khí ẩm nước ngoài, phần lớn phần mềm có quyền; việc trích xuất liệu để đưa vào tính tốn chu trình theo ý người sử dụng khó khăn Việc xây dựng đồ thị theo ý người dùng lại khó khăn Giảng viên sinh viên Khoa Công nghệ Nhiệt – Điện lạnh, Trường Đại học Bách khoa – Đại học Đà Nẵng tiến hành xây dựng hệ thống tính tốn thông số nhiệt vật lý môi chất phổ biến phục vụ giảng dạy, học tập sản xuất, cụ thể xây dựng thành công mô-đun phần mềm Water_IF97 để xác định thông số nhiệt động vật lý nước nước [1, 2]; mô-đun phần mềm Gaspro dùng để xác định tính chất nhiệt vật lý khói, khơng khí khơ, khí đơn chất Việc xây dựng mơđun phần mềm xác định thơng số khơng khí ẩm nằm hệ thống hoạt động nói Bài báo trình bày phương pháp xây dựng mô-đun phần mềm xác định thơng số khơng khí ẩm; ứng dụng để xây dựng đồ thị I-d cụ thể giải toán phun nước ngầm để làm mát chuồng trại 2.1.1 Cơ sở lý thuyết tính tốn thơng số khơng khí ẩm Tính tốn thơng số khơng khí ẩm liên quan đến áp suất bão hòa nước nhiệt độ khơng khí ẩm Giá trị xác định theo cơng thức Antoine [3]: 4026,42 𝑃𝑏 = 𝑒𝑥𝑝 (12 − ) (1) 235,5 + 𝑡 Cơng thức Antoine có sai số cực đại nhỏ 1% so với bảng thông số nhiệt động nước nước [2]; tính tốn kỹ thuật, giảng dạy học tập chấp nhận Các thơng số khơng khí ẩm bao gồm nhiệt độ t, độ ẩm tương đối φ, độ chứa (ẩm dung) d, enthalpy I Trạng thái không khí ẩm xác định biết thơng số nói Mối quan hệ chúng thể qua công thức đây: Độ chứa d [kg hơi/kg KKK] xác định theo công thức: 𝜑𝑝𝑏 𝑑 = 0,622 (2) 𝐵 − 𝜑𝑝𝑏 đó: Pb: Áp suất bão hòa nước nhiệt độ khơng khí ẩm Pa, xác định theo (1) B: Áp suất khí quyển, Pa Giá trị áp suất khí khác vùng giới, mơ-đun lấy B=105 Pa Enthalpy khơng khí ẩm I [kJ/kg KKK] xác định theo công thức: Kết nghiên cứu khảo sát 2.1 Xây dựng mô-đun phần mềm xác định thông số không khí ẩm 𝐼 = 𝑡 + 𝑑(2500 + 1,83𝑡) (3) Đối với nhiều tốn kỹ thuật, thơng thường cho trước nhiệt độ t độ ẩm tương đối φ khơng khí ẩm Từ (1) nhận thấy áp suất bão hòa nước pb hàm số theo nhiệt độ t, nên độ chứa d theo (2) hàm số ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 3(124).2018 𝑑 = 𝑓1 (𝑡, 𝜑); Enthalpy khơng khí ẩm theo (3) 𝐼 = 𝑓2 (𝑡, 𝑑) = 𝑓2 [𝑡, 𝑓1 (𝑡, 𝜑)] nên xác định biết t, φ Khi cho trước cặp thơng số khác, có trường hợp ta biến đổi từ phương trình (1) - (3) cho sau: 𝐼−𝑡 𝑑 = 𝑓3 (𝑡, 𝐼) = (4) 2500 + 1,83𝑡 𝐼 − 2500𝑑 (5) 𝑡 = 𝑓4 (𝐼, 𝑑) = + 1,83𝑑 𝑑𝐵 (6) 𝜑 = 𝑓5 (𝑡, 𝑑) = 𝑝𝑏 (0,622 + 𝑑) 4026,42 𝑡 = 𝑓6 (𝑑, 𝜑) = − 235,5 (7) 𝑑𝐵 12 − 𝑙𝑛 [ ] 𝜑(0,622 + 𝑑) Vậy cho trước cặp thơng số khơng khí ẩm, ta tính thơng số lại Bảng Bảng Hàm số xác định thông số khơng khí ẩm t, φ thường sử dụng để giải tốn Nó thuận tiện việc tính tốn, kể tính tốn theo vòng lặp chọn lựa điều kiện; dễ dàng trình bày cơng thức, bảng biểu Ngồi dùng Visual basic for application (VBA) để xử lý nhiều vấn đề phức tạp mặt kỹ thuật lập trình Do vậy, mơ-đun phần mềm tác giả viết ngôn ngữ VBA; file thu lưu giữ dạng add-in Excel dạng file.bas Visual Basic, tích hợp vào Excel phần mềm người dùng khác Mô-đun gồm hàm xác định thơng số khơng khí ẩm, xây dựng theo công thức nêu Mục 2.1.1 Để thống cách sử dụng với mô-đun xác định thông số môi chất phổ biến giảng viên Khoa Công nghệ Nhiệt – Điện lạnh nêu phần đặt vấn đề, tác giả quy ước đặt tên hàm mô-đun với cú pháp KKA_ThơngSốCầnTìm_TìmTheoĐốiSố (ĐốiSố1, ĐốiSố2) Bảng Danh mục hàm mô-đun KKA Cú pháp / Nhiệm vụ hàm Cặp thông số cho trước Thơng số cần tìm t, d t - - φ - 𝑓5 (𝑡, 𝑑) d 𝑓1 (𝑡, 𝜑) - I 𝑓2 (𝑡, 𝑑) = 𝑓2 [𝑡, 𝑓1 (𝑡, 𝜑)] 𝑓2 (𝑡, 𝑑) Cặp thông số cho trước Thông số cần tìm t, I I, d t - 𝑓4 (𝐼, 𝑑) φ 𝑓5 (𝑡, 𝑑) = 𝑓5 [𝑡, 𝑓3 (𝑡, 𝐼)] 𝑓5 (𝑡, 𝑑) = 𝑓5 [𝑓4 (𝐼, 𝑑), 𝑑] d 𝑓3 (𝑡, 𝐼) - I - Cặp thông số cho trước Thơng số cần tìm d, φ I, φ t 𝑓6 (𝑑, 𝜑) 𝑓6 (𝑑, 𝜑) = 𝑓6 [𝑓3 (𝑡, 𝐼), 𝜑] φ - - d - 𝑓3 (𝑡, 𝐼) I 𝑓2 (𝑡, 𝑑) = 𝑓2 [𝑓6 (𝑑, 𝜑), 𝑑] - Trong trường hợp đầu tiên, cho trước cặp thông số (t, φ), (t, d), (t, I), (I, d), (d, φ), cần thay vào hàm số đơn tính thơng số thứ Sau đó, thay thơng số tính vào hàm hợp lại cột, ta có thơng số thứ Trường hợp cuối cùng, cho trước cặp thơng số (I, φ), cần giải hệ phương trình liên quan đến ẩn số t, d Vấn đề giải cách lập trình theo phương pháp lặp với sai số tương đối 10-5 Phần mềm xây dựng phạm vi nhiệt độ 10≤t≤400°C; độ ẩm tương đối 0≤φ≤1, tất thông số đầu vào kiểm tra tính tương thích với phạm vi tính tốn phần mềm 2.1.2 Xây dựng mô-đun phần mềm Excel công cụ tin học mạnh phổ biến không lĩnh vực văn phòng, mà kỹ thuật 65 Đối số X t, [°C] φ, [-] KKA_d_XY(X,Y), [kg hơi/kg KKK] t, [°C] Xác định độ chứa I, [kJ/kg KKK] d theo X, Y I, [kJ/kg KKK] φ, [-] t, [°C] I, [kJ/kg KKK] t, [°C] φ, [-] φ, [-] I, [kJ/kg KKK] t, [°C] φ, [-] φ, [-] KKA_I_XY(X,Y), t, [°C] [kJ/kg KKK] Xác định enthalpy I d, [kg hơi/kg KKK] theo X, Y d, [kg hơi/kg KKK] t, [°C] d, [kg hơi/kg KKK] t, [°C] φ, [-] φ, [-] d, [kg hơi/kg KKK] I, [kJ/kg KKK] d, [kg hơi/kg KKK] d, [kg hơi/kg KKK] KKA_t_XY(X,Y), d, [kg hơi/kg KKK] [oC] Xác định nhiệt độ t φ, [-] theo X, Y I, [kJ/kg KKK] KKA_Fi_XY(X,Y), [_] Xác định độ ẩm tương đối φ theo X, Y Đối số Y I, [kJ/kg KKK] φ, [-] d, [kg hơi/kg KKK] φ, [-] φ, [-] I, [kJ/kg KKK] t, [°C] d, [kg hơi/kg KKK] d, [kg hơi/kg KKK] t, [°C] t, [°C] I, [kJ/kg KKK] I, [kJ/kg KKK] t, [°C] I, [kJ/kg KKK] d, [kg hơi/kg KKK] d, [kg hơi/kg KKK] I, [kJ/kg KKK] Ví dụ, để xác định enthalpy I (Thơng Số Cần Tìm) theo nhiệt độ t (Đối Số 1) độ ẩm tương đối φ (Đối Số 2), dùng hàm KKA_I_tFi(t, Fi) Ngoài ra, để tiện việc sử dụng, tác giả bổ sung thêm hàm theo thứ tự nhập đối Thái Ngọc Sơn 66 số ngược lại, ví dụ để xác định enthalpy I theo nhiệt độ t độ ẩm tương đối φ dùng hàm KKA_I_Fit(Fi, t) Danh mục hàm tác giả viết liệt kê Bảng Mô-đun đưa vào mơi trường Excel để kiểm tra tính xác Số liệu để so sánh toán khơng khí ẩm sách giáo trình, tập [4, 5, 6, 9, 10], đó, thơng số trạng thái phần lớn xác định theo đồ thị Kết cho thấy sai số lớn mô-đun tài liệu nhỏ 5% Tác giả so sánh kết sử dụng mô-đun với phần mềm hoạt động trực tuyến [7] Công ty Remak chuyên hệ thống điều hòa khơng khí Cộng hòa Czech Kết kiểm tra cho thấy sai số trung bình 1,3%; sai số lớn 2,86% Các sai số số cơng thức (1) - (7) có sai lệch định Tuy nhiên, sai số hồn tồn chấp nhận toán kỹ thuật nhiệt 2.2 Ứng dụng mơ-đun phần mềm để giải tốn chun ngành 2.2.1 Xây dựng đồ thị I-d khơng khí ẩm Đồ thị I-d tài liệu tham khảo, giáo trình cơng cụ đắc lực cho sinh viên, kỹ sư tính tốn khơng khí ẩm Tuy nhiên, đồ thị thường in với kích cỡ nhỏ, tra cứu gặp phải nhiều khó khăn, sai số lớn Mơ-đun khơng khí ẩm viết cài đặt vào Excel, xuất số liệu trang tính Tuy nhiên, trục tọa độ đồ thị I-d hệ tọa độ vng góc, nghĩa tung độ để vẽ đồ thị công cụ Excel giá trị I tính tốn mà cần biến đổi thành tung độ phù hợp Ta chuyển đổi giá trị enthlpy I tính tốn thành tung độ để vẽ đồ thị sau: Từ (3), ta có enthalpy điểm nằm đường đẳng nhiệt t1 xác định theo công thức: 𝐼1 = 2500𝑑1 + 𝑡1 + 1,83𝑡1 𝑑1 (8) Giá trị enthalpy chia làm phần (Hình 1) Vậy xem đường t=0 trùng với trục hồnh tọa độ dùng để vẽ đồ thị điểm 𝑥1 = 𝑑1 { (9) 𝑦1 = 𝐼1 − 2500𝑑1 = 𝑡1 + 1,83𝑡1 𝑑1 đường φ=const, t=const, I=const; biến đổi qua tọa độ để vẽ đồ thị theo (9); dùng công cụ Charts Excel, ta xây dựng đồ thị Hình Hình Đồ thị I-d Tùy theo nhu cầu người sử dụng, đặt kích thước hiển thị đồ thị, độ lớn trục, vị trí nhãn đường đẳng trị Ngồi ra, người sử dụng cho tọa độ điểm liên tiếp q trình khơng khí ẩm, mơ-đun hiển thị q trình đồ thị Giao diện sử dụng thể Hình 2.2.2 Ứng dụng mơ-đun để giải tốn chun ngành Hình Giao diện mơ-đun Hình Phân tích thành phần enthalpy khơng khí ẩm Dùng mơ-đun xây dựng, tính điểm a Bài tốn: Cho khơng khí nước trao đổi nhiệt – khối thiết bị trao đổi nhiệt hỗn hợp ngược chiều Biết khơng ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 3(124).2018 khí có lưu lượng vào thiết bị 35,1 kgKKK/s; nhiệt độ 34,5°C, φ = 55% Nước ngầm vào thiết bị có nhiệt độ 24°C Xác định: lưu lượng nước để không khí ẩm khỏi thiết bị có độ ẩm 85%; nhiệt độ khơng khí ẩm khỏi thiết bị b Giải: Đây toán cần giải nghiên cứu triển khai sử dụng nước ngầm để làm mát chuồng trại Độ ẩm khơng khí khỏi thiết bị bị khống chế theo điều kiện chăn nuôi Lưu lượng khơng khí vào thiết bị lưu lượng khơng khí tươi cần thiết để cung cấp cho chuồng trại Nhiệt độ đầu vào nước ngầm nhiệt độ khảo sát địa phương Giả thiết nhiệt độ nước khỏi thiết bị t" n Nhiệt độ t"n bị giới hạn nhiệt độ nhiệt kế ướt, phụ thuộc vào thơng số khơng khí ẩm Giả sử ta chia buồng phun thành n cấp, biến thiên enthalpy khơng khí ẩm cấp (Hình 4) Phương trình cân nhiệt cho toàn thiết bị: ′ 𝐺𝑘 𝐼 + 𝐺𝑛 𝐶𝑛 𝑡 ′ 𝑛 = 𝐺𝑘 𝐼" + (𝐺𝑛 + ∆𝑊)𝐶𝑛 𝑡"𝑛 (10) (𝐺𝑛 + ∆𝑊)𝐶𝑛 𝑡"𝑛 − 𝐺𝑛 𝐶𝑛 𝑡 ′ 𝑛 (11) 𝐼 ′ − 𝐼" = 𝐺𝑘 đó, Cn: nhiệt dung riêng nước, kJ/kgK; ΔW: lượng nước bay vào khơng khí ngưng tụ từ khơng khí, kgh/s Trong tốn thực tế thường xem ΔW=0 [10] I', I": enthalpy khơng khí ẩm vào khỏi thiết bị, kJ/kg KKK Biến thiên enthalpy cấp: 𝐼 ′ − 𝐼" ∆𝐼 = (12) 𝑛 đó, n số cấp tự chọn để tính tốn; n lớn kết xác Tuy nhiên, với mức độ yêu cầu xác tốn kỹ thuật nhiệt, ta cần chọn số cấp quãng từ đến 20 độ chênh nhiệt độ đầu vào khơng khí thiết bị từ đến 100°C 67 𝐼𝐴n = 𝐼𝐴n−1 − ∆𝐼 = 𝐼" Nhiệt độ nước tiết diện 0: 𝑡𝐾0 = 𝑡𝑛′′ (14) Nhiệt độ nước tK1 tiết diện xác định dựa phương trình cân nhiệt cho đoạn 0-1 thiết bị trao đổi nhiệt sau: 𝐺𝑘 𝐼𝐴0 + 𝐺𝑛 𝐶𝑛 𝑡𝐾1 = 𝐺𝑘 𝐼𝐴1 + 𝐺𝑛 𝐶𝑛 𝑡𝐾0 (15) từ 𝑡𝐾1 = 𝑡𝐾0 − 𝐺𝑘 𝐺𝑘 (𝐼𝐴0 − 𝐼𝐴1 ) = 𝑡𝐾0 − ∆𝐼 𝐺𝑛 𝐶𝑛 𝐺𝑛 𝐶𝑛 (16) Đặt 𝐺𝑘 ∆𝐼 𝐺𝑛 𝐶𝑛 = 𝑡𝐾0 − ∆𝑡𝑛 ∆𝑡𝑛 ≜ ⇒ 𝑡𝐾1 (17) (18) Tương tự 𝑡𝐾2 = 𝑡𝐾1 − ∆𝑡𝑛 𝑡𝐾𝑛 = 𝑡𝐾𝑛−1 − ∆𝑡𝑛 = (19) 𝑡𝑛′ (20) Khơng khí ẩm vào đoạn 0-1 (Điểm A0) hòa trộn với khơng khí ẩm trạng thái bão hòa bề mặt sát mặt nước nhiệt độ tK0 (Điểm K0) Điểm A1 giao đường A0K0 đường IA1 Đây điểm xác định trạng thái khơng khí ẩm khỏi đoạn 0-1 vào đoạn 1-2 Tương tự ta xác định tất điểm lại Bảng Tính tốn biến thiên thơng số khơng khí ẩm sơ đồ ngược chiều Diễn giải Tiết diện Tiết diện Điểm A0 A1 A2 A3 A4 A5 Nhiệt độ KKA tk' 34,5 29,6 29 28,7 28,5 28,3 Enthalpy KKA Ik' 83,7 83,3 83,0 82,6 82,2 81,9 Ẩm dung KKA Ik' 19,1 21,0 21,1 21,1 21,0 20,9 Độ ẩm tương đối KKA φk' 0,55 0,79 0,82 0,83 0,84 0,85 Điểm K0 K1 K2 K3 K4 K5 Nhiệt độ nước tn" 26,5 26 25,5 25 24,5 24,0 Enthalpy KKA bão hòa In" 83,1 80,9 78,7 76,6 74,5 Ẩm dung KKA bão hòa dn" 22,2 21,5 20,9 20,2 19,6 Tính: Điểm A1 A2 A3 A4 A5 Ẩm dung điểm hòa trộn dk" 21,0 21,1 21,1 21,0 20,9 Hình Tính biến thiên nhiệt độ sơ đồ ngược chiều Enthalpy KKA k" 83,3 83,0 82,6 82,2 81,9 Enthalpy khơng khí ẩm tiết diện 0, 1, , n: 𝐼 ′ = 𝐼𝐴0 𝐼𝐴1 = 𝐼𝐴0 − ∆𝐼 (13) ……… Nhiệt độ KKA k" 29,6 29 28,7 28,5 28,3 Độ ẩm tương đối KKA φk" 0,79 0,82 0,83 0,84 0,85 Việc tính tốn thơng số khơng khí ẩm điểm nút thực theo phương pháp lặp: Chọn lưu Thái Ngọc Sơn 68 lượng nước, tính tốn thơng số nước khơng khí ẩm tiết diện dựa quan hệ (13) - (20) môđun phần mềm; Kiểm tra độ ẩm khơng khí ẩm đầu thiết bị đạt 85% Khi lưu lượng nước Gn=6kg/s Khơng khí ẩm khỏi thiết bị với thơng số điểm A5 Khi nhiệt độ khơng khí ẩm 28,3°C Số liệu tính tốn biến thiên thơng số khơng khí ẩm thiết bị trao đổi nhiệt hỗn hợp với sơ đồ ngược chiều biểu thị Bảng 5, đơn vị nhiệt độ t [°C]; enthalpy I [kJ/kgKKK]; ẩm dung d [g hơi/kgKKK]; Q trình biến đổi trạng thái khơng khí ẩm thể đồ thị I-d Hình d với tỷ lệ xích mong muốn, làm cho việc phân tích q trình từ đồ thị trở nên dễ dàng Kết luận Tác giả xây dựng dược mô-đun phần mềm xác định thông số không khí ẩm với độ xác cao, dễ sử dụng, tích hợp vào chương trình người dùng khác vào Excel, cơng cụ hữu ích phục vụ giảng viên, sinh viên giảng dạy học tập kỹ sư cần thiết giải tốn liên quan đến khơng khí ẩm Tác giả xây dựng đồ thị I-d dựa mơ-đun phần mềm xác định thơng số khơng khí ẩm Người dùng thay đổi linh hoạt tỷ lệ xích đồ thị, hiển thị q trình tính tốn, tăng tính xác, giảm cơng sức lao động Tác giả áp dụng giải toán phức tạp, liên quan đến việc sử dụng nước ngầm để làm mát khơng khí Đây sở để tính tốn thiết kế xây dụng hệ thống làm mát chuồng trại nguồn nước ngầm ngày hè nóng nực miền Trung TÀI LIỆU THAM KHẢO Hình Biểu diễn q trình làm mát khơng khí nước ngầm sơ đồ ngược chiều đồ thị I-d Đường trình đồ thị xây dựng cách nối điểm với nhau; điều chỉnh hệ tọa độ hiển thị cho phù hợp với yêu cầu tốn Đây điểm mạnh mơ-đun phần mềm Nó giúp cho người dùng dễ dàng xây dụng trình cần thiết đồ thị I- [1] Thái Ngọc Sơn, Xây dựng chương trình tính tốn thơng số kỹ thuật nước nước, Đề tài nghiên cứu khoa học cấp sở, Trường Đại học Bách khoa – Đại học Đà Nẵng, 2010 [2] Hoàng Ngọc Đồng, Thái Ngọc Sơn, Nguyễn Quốc Huy, “Xây dựng module xác định thông số nhiệt động nước nước”, Năng lượng nhiệt, Số 103-1/2012, Hà Nội, 2012, trang 14-16 [3] Trần Văn Phú, Kỹ thuật sấy, NXB Giáo dục, Hà Nội, 2009 [4] Bùi Hải, Dương Đức Hồng, Hà Mạnh Thư, Thiết bị trao đổi nhiệt, NXB Khoa học Kỹ thuật, 1999 [5] Богословский В.Н и др, Кондиционирование воздуха и холодоснбжение, Москва, Стройиздат, 1985 [6] Hoàng Ngọc Đồng, Thái Ngọc Sơn, Kỹ thuật nhiệt, NXB Xây dựng, 2015 [7] Remak, Calculation of moist air properties, http://www.remak.eu/en/calculation-moist-air-properties [Truy cập 17/01/2018] [8] Phan Tự Hướng, Lập trình VBA (Visual Basic For Application) Excel, NXB Thống kê, Hà Nội, 2010 [9] Ерофеев В.Л., Семенов П.Д., Пряхин А.С., Теплотехника, Москва, Академкнига, 2008 [10] Лебедев П Д., Теплообменные, сушильные и холодильные установки, Москва, Энергия, 1972 (BBT nhận bài: 24/01/2018, hoàn tất thủ tục phản biện: 22/02/2018) ... Vậy cho trước cặp thông số khơng khí ẩm, ta tính thơng số lại Bảng Bảng Hàm số xác định thông số không khí ẩm t, φ thường sử dụng để giải tốn Nó thuận tiện việc tính tốn, kể tính tốn theo vòng... dựng phạm vi nhiệt độ 10≤t≤400°C; độ ẩm tương đối 0≤φ≤1, tất thông số đầu vào kiểm tra tính tương thích với phạm vi tính tốn phần mềm 2.1.2 Xây dựng mơ -đun phần mềm Excel công cụ tin học mạnh phổ... 2.2 Ứng dụng mơ -đun phần mềm để giải toán chuyên ngành 2.2.1 Xây dựng đồ thị I-d khơng khí ẩm Đồ thị I-d tài liệu tham khảo, giáo trình cơng cụ đắc lực cho sinh viên, kỹ sư tính tốn khơng khí ẩm