Nghiên cứu thực nghiệm collector mặt trời cấp nhiệt cho máy lạnh hấp thụ NH3 h2o loại gián đoạn để sản xuất nước đá

1 LỜI CÁM ƠN Trong suốt thời gian thực luận văn tác giả nhận nhiều giúp đỡ mặt vật chất lẫn tinh thần thầy cô, anh chị em môn Công Nghệ Nhiệt Lạnh, bạn lớp gia đình Trước tiên tác giả xin gửi lời cảm ơn đến PGS.TS Lê Chí Hiệp, người hướng dẫn kỹ thuật hỗ trợ phần kinh phí cho luận văn này, tác giả xin gửi lời cảm ơn đến TS Bùi Ngọc Hùng, người hỗ trợ kinh phí thiết bị cho mô hình thực nghiệm Cám ơn thầy Nguyễn Tấn Phương, bạn Nguyễn Văn Hạp, anh em công ty nhiệt lạnh Bách Khoa hỗ trợ việc chế tạo hoàn tất mô hình Tác giả xin cảm ơn thầy cô anh chị đồng nghiệp môn Công Nghệ Nhiệt Lạnh đóng góp ý kiến quý báo Và sau xin cảm ơn gia đình động viên tinh thần góp phần cho thành công luận văn TÓM TẮT NỘI DUNG LUẬN VĂN Trong luận văn tác giả nghiên cứu, trình bày phương pháp tính toán, thiết kế chế tạo mô hình máy lạnh hấp thụ NH3-H2O cấp nhiệt lượng mặt trời loại gián đoạn để sản xuất nước đá dùng collector loại phẳng Từ suất số liệu lựa chọn ban đầu, tác giả trình bày phương pháp tính toán thông số chu trình như: suất lạnh, khối lượng tác nhân lạnh, khối lượng dung dịch, Trên sở số liệu vừa tính toán tác giả trình bày phương pháp tính toán kích thước bề mặt trao đổi nhiệt mô hình như: dàn ngưng tụ, collector,… Sau tính toán nhiệt chu trình kiểm tra lại diện tích trao đổi nhiệt collector Từ đưa mô hình thực nghiệm Trong luận văn tác giả nghiên cứu thực nghiệm mô hình máy lạnh hấp thụ NH3-H2O loại gián đoạn cấp nhiệt collector dạng phẳng Dựa số liệu thực nghiệm đo từ mô hình tác giả so sánh số liệu thực nghiệm với lý thuyết lựa chọn ban đầu để làm sở cho việc tính toán đánh giá Trong luận văn tác giả trình bày phương pháp tính toán nhiệt collector từ xây dựng chương trình mô máy lạnh hấp thụ NH3-H2O cấp nhiệt lượng mặt trời loại gián đoạn Từ số liệu chương trình mô tác giả so sánh với số liệu đo đạc từ mô hình thực nghiệm đưa nhận xét kết luận Ngoài tác giả trình bày quy trình vận hành thông qua kinh nghiệm thực nghiệm đánh giá khả ứng dụng máy lạnh hấp thụ NH3-H2O loại gián đoạn điều kiện Việt Nam Từ mô hình thực tác giả đề xuất số cải tiến nhằm nâng cao hiệu làm việc hệ thống MỤC LỤC Trang Chương 1: TỔNG QUAN VỀ MÁY LẠNH HẤP THỤ SỬ DỤNG NĂNG LƯNG MẶT TRỜI 1.1 Tổng quan xạ mặt trời öùng duïng 1.2 Tiềm lượng mặt trời Việt Nam 13 1.3 Tình hình sử dụng lượng mặt trời giới Việt Nam 19 1.4 Các hệ thống máy lạnh hấp thụ NH3-H2O sản xuất nước đá vận hành lượng mặt trời thực 21 1.5 Mục tiêu phương pháp nghiên cứu 24 Chương 2: CÁC DẠNG COLLECTOR MẶT TRỜI VÀ CÁC ƯU NHƯC ĐIỂM 26 2.1 Collector phẳng 26 2.2 Collector daïng ống hút chân không 28 2.3 Collector dạng tập trung 29 2.4 Phaân tích lựa chọn phương án 30 Chương 3: MÁY LẠNH HẤP THỤ NH3-H2O 33 3.1 Tổng quát máy hấp thụ NH3-H2O 33 3.1.1 Sự khác biệt máy lạnh hấp thụ máy lạnh có máy nén .33 3.1.2 Dung dịch làm việc, nồng độ dung dịch .33 3.1.3 Nguyên lý làm việc máy lạnh hấp thụ NH3-H2O 34 3.1.4 Đồ thị i-c dung dịch 35 3.1.5 Caùc công thức tính toán thông số nhiệt động 36 3.2 Máy lạnh hấp thụ NH3-H2O cấp nhiệt lượng mặt trời 39 3.2.1 Nguyên lý hoạt động 39 3.2.2 Máy lạnh hấp thụ NH3-H2O loại liên tục cấp nhiệt lượng mặt trời 40 3.2.3 Máy lạnh hấp thụ NH3-H2O loại gián đoạn cấp nhiệt lượng mặt trời 40 3.2.4 So sánh chu trình liên tục chu trình gián đoạn 42 Chương 4: THIẾT KẾ, LẮP ĐẶT MÔ HÌNH VÀ ĐO ĐẠC THỰC NGHIỆM .43 4.1 Thiết kế, lắp đặt mô hình 43 4.1.1 Collector: 43 4.1.2 Dàn ngưng tụ: .43 4.1.3 Bình chứa: 44 4.1.4 Van tiết lưu: 44 4.1.5 Dàn bay hơi: 45 4.2 Tính toán thiết keá 45 4.2.1 Chọn thông số ban đầu 45 4.2.2 Xác định suất lạnh .45 4.2.3 Xác định áp suất ngưng tuï 46 4.2.4 Xác định áp suất bay 46 4.2.5 Xác định phụ tải lạnh 46 4.2.6 Xác định thể tích bình chứa NH3 sau ngưng tụ .47 4.2.7 Xác định phụ tải nhiệt ngưng tụ 47 4.2.8 Tính toán kích thước bình nước giải nhiệt: 48 4.2.9 Tính toán diện tích trao đổi nhiệt dàn ngưng tụ 48 4.2.10 Tính kích thước bình bay 52 4.2.11 Tính suất lạnh dàn lạnh .52 4.2.12 Chọn nồng độ dung dịch đậm đặc nồng độ dung dịch loãng 52 4.2.13 Tính khối lượng dung dịch chứa collector 53 4.2.14 Tính toán kích thước collector 53 4.3 Tổng nhiệt lượng mà collector nhận từ xạ mặt trời: 56 4.3.1 Tính toán thông số trạng thái dung dịch thời điểm: 56 4.3.2 Mô tả trình làm việc máy lạnh hấp thụ NH3-H2O loại gián đoạn đồ thị i-c 64 4.4 Tính nhiệt lượng nhận vào nhả dung dịch chứa collector 69 4.5 Cân nhiệt toàn hệ thoáng 71 4.6 Tính toán kiểm tra diện tích trao đổi nhiệt collector: 74 4.6.1 Tính kiểm tra diện tích trao đổi nhiệt collector vào ban ngày .74 4.6.2 Tính kiểm tra khả giải nhiệt vào ban đêm collector: 76 4.6.3 Cách bố trí dụng cụ đo số liệu đo từ mô hình thực nghiệm 77 4.6.4 Kết thực hiện: 85 4.7 So sánh kết đo đạc thực nghiệm lý thuyết 86 Chương : CHƯƠNG TRÌNH MÔ PHỎNG 89 5.1 Giới thiệu chương trình 89 5.2 Cơ sở lý thuyết bước tính toán 89 5.2.1 Tính toán nhiệt collector 89 5.2.2 Các bước thực .100 5.3 Chương trình mô 104 5.4 Caùc kết từ chương trình mô nhận xét 105 Chương 6: PHÂN TÍCH ĐÁNH GIÁ KHẢ NĂNG ỨNG DỤNG MÁY LẠNH HẤP THỤ NH3-H2O LOẠI GIÁN ĐOẠN TRONG ĐIỀU KIỆN VIỆT NAM 109 KẾT LUẬN 114 QUY TRÌNH VẬN HÀNH 117 TÀI LIỆU THAM KHẢO 121 10 PHUÏ LUÏC .122 Các hệ thống máy lạnh hấp thụ NH3-H2O loại gián đoạn thực giới 122 Caùc kết đo từ thực nghiệm 125 Code chương trình mô phoûng 128 Code chương trình tính toán entanpi 144 Chương trình tính thông số vật lý dung dịch 146 MỞ ĐẦU Ngày nay, lượng mối quan tâm hàng đầu quốc gia giới Cùng với phát triển kinh tế nhu cầu tiêu thụ lượng tăng theo Trong nguồn nhiên liệu truyền thống như: dầu mỏ, than đá, khí thiên nhiên có giới hạn cạn kiệt dần Theo ước tính với tình hình sử dụng lượng khoảng 50 đến 60 năm nguồn nhiên liệu cạn kiệt Bên cạnh vấn đề ô nhiễm môi trường việc sử dụng nhiên liệu hoá thạch đặt toán khó cho người phát triển bền vững Để giải vấn đề này, người tìm kiếm nguồn lượng thay như: lượng hạt nhân, lượng mặt trời, lượng gió, lượng địa nhiệt, lượng thuỷ triều,…Trong lượng mặt trời ứng dụng rộng rãi Năng lượng mặt trời nguồn lượng vô tận có mặt khắp nơi trái đất, nguồn lượng hoàn toàn miễn phí Trong bối cảnh sử dụng lượng việc nghiên cứu ứng dụng lượng mặt trời vào thực tế việc làm cần thiết có ý nghóa Đối với nước ta nước có khí hậu nhiệt đới, gần đường xích đạo nên lượng xạ mặt trời năm tương đối cao Thời gian chiếu sáng trung bình năm từ 2000 đến 2600 giờ, với mật độ nằm khoảng từ 420 đến 730 kJ/cm2/năm Đặc biệt tỉnh phía Nam lượng xạ trung bình tương đối cao không thay đổi nhiều thích hợp cho việc ứng dụng lượng mặt trời vào nhiều mục đích khác Hiện việc ứng dụng lượng mặt trời vào thực tế nước ta phần lớn tập trung hệ thống cung cấp nước nóng lượng mặt trời quy mô nhỏ Tuy nhiên việc ứng dụng lượng mặt trời vào mục đích cung cấp nước nóng thực chưa mang lại hiệu cao hiệu suất thiết bị thấp nhiệt độ nước nóng đạt không cao, thích hợp cho sinh hoạt hàng ngày Ngoài nhược điểm khác thiết bị phụ thuộc vào thời tiết, việc ứng dụng thiết bị cung cấp nước nóng lượng mặt trời chưa quan tâm đáng kể Do việc nghiên cứu tìm giải pháp, công nghệ ứng dụng lượng mặt trời vào mục đích khác như: sử dụng lượng mặt trời cấp nhiệt cho máy lạnh hấp thụ, chưng cất nước, bơm nước,… hướng nghiên cứu nhà khoa học quan tâm Hệ thống máy lạnh hấp thụ sử dụng lượng mặt trời nghiên cứu chế tạo thành công số nước giới Các hệ thống máy lạnh hấp thụ sử dụng lượng mặt trời chủ yếu dùng để điều hoà không khí với thu collector phẳng Ngoài có số hệ thống dùng để sản xuất nước đá với thu collector tập trung Tuy nhiên vấn đề xa lạ mẻ nước ta Để tránh tụt hậu so với nước khu vực giới, cần phải đẩy mạnh việc nghiên cứu ứng dụng lượng mặt trời vào thực tiễn Bên cạnh tảng hệ thống cung cấp nước nóng lượng mặt trời, cần tập trung vào ứng dụng khác máy lạnh hấp thụ, pin mặt trời,… Trong luận văn tác giả nghiên cứu thực nghiệm collector mặt trời cấp nhiệt cho máy lạnh hấp thụ loại gián đoạn để sản xuất nước đá Chương 1: TỔNG QUAN VỀ MÁY LẠNH HẤP THỤ SỬ DỤNG NĂNG LƯNG MẶT TRỜI 1.1 Tổng quan xạ mặt trời ứng dụng Mặt trời cầu khổng lồ có đường kính 1392.106km cách trái đất khoảng 149,6.106km, nhiệt độ bề mặt khoảng 5800K Tại tâm mặt trời nhiệt độ lên đến 15 triệu độ áp suất khoảng hàng trăm triệu atm Ở nhiệt độ áp suất hạt nhân nguyên tố nhẹ Hydro thu động lớn đủ để thắng lực đẩy tónh điện kết hợp chặt với tạo thành hạt nhân Heli (được gọi phản ứng nhiệt hạt nhân hay nổ hạt nhân) Chính phản ứng giải phóng lượng vô lớn trì hoạt động mặt trời Với nhiệt độ bề mặt mặt trời xạ lượng từ bề mặt 62MW/m2 Điều kiểm chứng lại định luật Stefan-Boltzmann xem mặt trời vật đen tuyệt đối ta có: E=σ.T4=5,67.10-8.57604=62,412MW/m2 Nếu đường kính mặt trời 1,392.109m tổng xạ mà mặt trời toả không trung là: ΣE=F.E=πd2.E=3,14159.(1,392.109)2.62,412=3,8.1020MW Hình 1.1: Vị trí tương đối mặt trời trái đất DS=1,392.109m RE=6,367.106m d=1,495.1011m Hình khoảng cách mặt trời trái đất 1,495.1011m Như xạ mặt trời đến bên bầu khí tính sau: D s2 = 6,241.10 =E 4.π d 4.π E out ⎛ 1,392.10 ⎞ ⎜⎜ ⎟ = 1353W/m 11 ⎟ ⎝ 1,495.10 ⎠ Khoảng cách trung bình từ trái đất đến mặt trời 150 triệu km Tuy nhiên quỹ đạo quay trái đất quanh mặt trời hình elip khoảng cách không cố định mà thay đổi theo mùa năm Vào lúc trái đất cách xa mặt trời khoảng cách mặt trời trái đất 152 triệu km Vào lúc gần 147 triệu km Do giá trị Eout thay đổi theo, nhiên thay đổi không lớn nên giá trị Eout xem số gọi số mặt trời Tổng lượng mà trái đất nhận xạ mặt trời tính toán sau: ( P = E out π R 2E = 1353.3,14159 6,367.10 ) = 1,7232.1011 MW Trong ngày tổng lượng xạ mặt trời đến bên bầu khí trái đất là: E day = E out 24.3600 = 1,489.1016 MJ Cường độ xạ mặt trời mặt đất thay đổi theo vó độ, vị trí địa lý, mùa năm, tỷ lệ mây bao phủ, ô nhiễm bầu khí quyển, độ cao so với mặt nước biển vó độ mặt trời Góc nghiêng 23,5o trục trái đất ảnh hưởng đến góc tới của tia xạ mặt trời đến bề mặt trái đất, thời tiết mùa năm độ dài ngày năm Năng lượng mặt trời truyền đến trái đất dạng xạ nhiệt Bức xạ mặt trời dạng lượng truyền qua môi trường chân không 10 Bức xạ mặt trời truyền dạng sóng ánh sáng với vận tốc 300.000km/s Bức xạ mặt trời xem hạt nhỏ gọi photon chuyển động với vận tốc ánh sáng Photon đơn vị nhỏ xạ mặt trời Tập họp nhiều photon tạo thành sóng ánh sáng, tập họp sóng tạo thành chùm tia xạ Sự khác biệt loại sóng bước sóng chúng Bước sóng xạ mặt trời tỷ lệ nghịch với tần số xạ mặt trời Ta có quan hệ bước sóng tần số sau: λ= Trong đó: c f (1.1) λ bước sóng, m c vận tốc ánh sáng, m/s f tần số, s-1 Năng lượng đơn vị xạ mặt trời (photon) phụ thuộc vào bước sóng tần số sau: E = hf = hc λ (1.2) Với h số Planck, h=6,63.10-34J.s Năng lượng xạ mặt trời tỷ lệ nghịch với với bước sóng Tia có bước sóng ngắn có mức lượng cao Tuy nhiên xạ mặt trời phát phần lớn nằm vùng có bước sóng dài Bức xạ điện từ phát từ mặt trời có dãy bước bước sóng rộng Nó chia làm vùng chủ yếu vùng xạ ion hoá (tia X, tia gamma) vùng xạ không ion hoá (bức xạ cực tím, xạ nhìn thấy xạ hồng ngoại) Những tia xạ ion hoá không thâm nhập qua tầng khí trái đất Bức xạ mặt trời thường chia làm nhiều vùng dãy bước soùng 135 Dim v2 As Double v1 = Text30 + (Text28 - Text31) * (Text36 - Text30) / (Text37 - Text31) v2 = Text33 + (Text28 - Text34) * (Text39 - Text33) / (Text40 - Text34) Dim v As Double v = v1 + (Text29 - Text32) * (v2 - v1) / (Text35 - Text32) 'The tich bay hoi Dim TT As Double TT = 0.016 Dim GNH3 As Double GNH3 = TT / v 'Nong cua dung dich sau bay hoi Dim cx As Double cx = Round((c * 15.5 - GNH3) / 15.5, 3) End If If (PAA - Pk) > Then P = Pk Dim acc As Double Dim bcc As Double Dim CAA As Double For CAA = 0.01 To Step 0.00001 acc = -2103.5 + 4669.96 * CAA - 20228.3 * (CAA ^ 2) + 56507 * (CAA ^ 3) 80989.9 * (CAA ^ 4) + 55286.5 * (CAA ^ 5) - 14361.4 * (CAA ^ 6) bcc = 5.65208 - 7.0317 * CAA + 37.9018 * (CAA ^ 2) - 102.912 * (CAA ^ 3) + 135.789 * (CAA ^ 4) - 82.7106 * (CAA ^ 5) + 18.4113 * (CAA ^ 6) Dim CCCC As Double CCCC = acc / ((Log(P) / Log(10)) + 0.00847711 - bcc) If (CCCC - T1) > And (CCCC - T1) < 0.01 Then Exit For End If Next CAA cx = Round(CAA, 4) GNH3 = 15.5 * c - 15.5 * cx End If 'TINH SO C1 Dim c1 As Double c1 = -6571.06 + 39.9544 * T1 - 0.00243781 * (T1 ^ 2) + 0.00000039792 * (T1 ^ 3) - 22.7722 * P - 0.00000004979 * (P ^ 2) + 8286.02 * (P / T1) + 0.0172363 * P * T1 + 0.00000077344 * P * (T1 ^ 2) - 3.62962 * T1 * (Log(T1)) - T1 * (Log(P)) 136 Dim d As Double d = (-1467 * (cx ^ 2) * (1.25 * cx ^ - 5.58 * cx ^ + 5.96 * cx - 0.42)) + ((17.3 * (cx ^ 3) - 13.8 * (cx ^ 4)) * T) 'TINH ENTANPI CUA NH3 Dim inh3 As Double inh3 = Round(-1058.43 + 3.02148 * T1 + 0.00312072 * T1 ^ + 0.15536 * P 0.0000123079 * P ^ - 0.0000012946 * P * T1 ^ 2, 2) 'TINH ENTANPI CUA NUOC Dim ih2o As Double ih2o = Round(-1039.5 + 3.56393 * T1 + 0.000919594 * T1 ^ + 0.124376 * P 0.00000222979 * P ^ - 0.00000035696 * P * T1 ^ 2, 2) 'TINH qt Dim qtt As Double qtt = Round(-839.87 * (1 - 0.071 * cx) * ((1 - cx) * cx + (1 - cx ^ 2) * cx ^ 2), 2) 'ENTANPI CUA DUNG DICH LONG Dim i1 As Double i1 = Round((cx * inh3 + (1 - cx) * ih2o + qtt), 3) 'TINH NONG DO CUA HOI Dim cv As Double cv = - (1 - cx) * Exp((c1 + d) / T1) 'TINH IVNH3 Dim ivnh3 As Double ivnh3 = Round(770.761 + 1.86947 * T1 + 0.0000578293 * T1 ^ + 0.000000731509 * T1 ^ + 8.98074 * P - 4580.15 * P / T1, 2) 'TINH IVH2O Dim ivh2o As Double ivh2o = Round(1993.19 + 1.88878 * T1 - 0.000205512 * T1 ^ + 0.000000367295 * T1 ^ + 10.6342 * P - 7648.34 * P / T1, 2) 'ENTANPI CUA HOI Dim i2 As Double i2 = Round((cv * ivnh3 + (1 - cv) * ivh2o), 3) Dim Q1 As Double Q1 = (15.5 - GNH3) * i1 + GNH3 * i2 'Entanpi cua dung dich o trang thai tdd = + 273 - 0.5 Dim a2 As Double Dim b2 As Double 137 a2 = -2103.5 + 4669.96 * c - 20228.3 * (c ^ 2) + 56507 * (c ^ 3) - 80989.9 * (c ^ 4) + 55286.5 * (c ^ 5) - 14361.4 * (c ^ 6) b2 = 5.65208 - 7.0317 * c + 37.9018 * (c ^ 2) - 102.912 * (c ^ 3) + 135.789 * (c ^ 4) - 82.7106 * (c ^ 5) + 18.4113 * (c ^ 6) Dim P2 As Double Dim PA As Double PA = Round(10 ^ ((a2 / tdd) + b2 - 0.00847711), 3) If (PA - Pk) Then P2 = Pk Dim CA As Double For CA = 0.01 To Step 0.00001 Dim ac As Double Dim bc As Double Dim CCC As Double ac = -2103.5 + 4669.96 * CA - 20228.3 * (CA ^ 2) + 56507 * (CA ^ 3) - 80989.9 * (CA ^ 4) + 55286.5 * (CA ^ 5) - 14361.4 * (CA ^ 6) bc = 5.65208 - 7.0317 * CA + 37.9018 * (CA ^ 2) - 102.912 * (CA ^ 3) + 135.789 * (CA ^ 4) - 82.7106 * (CA ^ 5) + 18.4113 * (CA ^ 6) CCC = ac / ((Log(P2) / Log(10)) + 0.00847711 - bc) If (CCC - tdd) > And (CCC - tdd) < 0.01 Then Exit For End If Next CA cxx = Round(CA, 4) GNH3x = 15.5 * c - 15.5 * cxx End If 'TINH SO C1 Dim c2 As Double c2 = -6571.06 + 39.9544 * tdd - 0.00243781 * (tdd ^ 2) + 0.00000039792 * (tdd ^ 3) - 22.7722 * P2 - 0.00000004979 * (P2 ^ 2) + 8286.02 * (P2 / tdd) + 0.0172363 * P2 * tdd + 0.00000077344 * P2 * (tdd ^ 2) - 3.62962 * tdd * (Log(tdd)) - tdd * (Log(P2)) Dim d2 As Double d2 = (-1467 * (cxx ^ 2) * (1.25 * cxx ^ - 5.58 * cxx ^ + 5.96 * cxx - 0.42)) + ((17.3 * (cxx ^ 3) - 13.8 * (cxx ^ 4)) * tdd) 'TINH ENTANPI CUA NH3 Dim inh3x As Double inh3x = Round(-1058.43 + 3.02148 * tdd + 0.00312072 * tdd ^ + 0.15536 * P2 0.0000123079 * P2 ^ - 0.0000012946 * P2 * tdd ^ 2, 2) 'TINH ENTANPI CUA NUOC Dim ih2ox As Double ih2ox = Round(-1039.5 + 3.56393 * tdd + 0.000919594 * tdd ^ + 0.124376 * P2 - 0.00000222979 * P2 ^ - 0.00000035696 * P2 * tdd ^ 2, 2) 'TINH qt Dim qttx As Double 139 qttx = Round(-839.87 * (1 - 0.071 * cxx) * ((1 - cxx) * cxx + (1 - cxx ^ 2) * cxx ^ 2), 2) 'ENTANPI CUA DUNG DICH LONG Dim i3 As Double i3 = Round((cxx * inh3x + (1 - cxx) * ih2ox + qttx), 3) 'TINH NONG DO CUA HOI Dim cvx As Double cvx = - (1 - cxx) * Exp((c2 + d2) / tdd) 'TINH IVNH3 Dim ivnh3x As Double ivnh3x = Round(770.761 + 1.86947 * tdd + 0.0000578293 * tdd ^ + 0.000000731509 * tdd ^ + 8.98074 * P2 - 4580.15 * P2 / tdd, 2) 'TINH IVH2O Dim ivh2ox As Double ivh2ox = Round(1993.19 + 1.88878 * tdd - 0.000205512 * tdd ^ + 0.000000367295 * tdd ^ + 10.6342 * P2 - 7648.34 * P2 / tdd, 2) 'ENTANPI CUA HOI Dim i4 As Double i4 = Round((cvx * ivnh3x + (1 - cvx) * ivh2ox), 3) Dim Q2 As Double Q2 = (15.5 - GNH3x) * i3 + GNH3x * i4 Dim Qu As Double Qu = (Q2 - Q1) + TT * (P - P2) * 100 If (Qu + Ql + Qs > Q) Then Exit For End If Next For tdd = + 273 - 1.5 To + 273 - 0.5 Step 0.01 a2 = -2103.5 + 4669.96 * c - 20228.3 * (c ^ 2) + 56507 * (c ^ 3) - 80989.9 * (c ^ 4) + 55286.5 * (c ^ 5) - 14361.4 * (c ^ 6) b2 = 5.65208 - 7.0317 * c + 37.9018 * (c ^ 2) - 102.912 * (c ^ 3) + 135.789 * (c ^ 4) - 82.7106 * (c ^ 5) + 18.4113 * (c ^ 6) PA = Round(10 ^ ((a2 / tdd) + b2 - 0.00847711), 3) If (PA - Pk) Then P2 = Pk For CA = 0.01 To Step 0.00001 ac = -2103.5 + 4669.96 * CA - 20228.3 * (CA ^ 2) + 56507 * (CA ^ 3) - 80989.9 * (CA ^ 4) + 55286.5 * (CA ^ 5) - 14361.4 * (CA ^ 6) bc = 5.65208 - 7.0317 * CA + 37.9018 * (CA ^ 2) - 102.912 * (CA ^ 3) + 135.789 * (CA ^ 4) - 82.7106 * (CA ^ 5) + 18.4113 * (CA ^ 6) CCC = ac / ((Log(P2) / Log(10)) + 0.00847711 - bc) If (CCC - tdd) > And (CCC - tdd) < 0.01 Then Exit For End If Next CA cxx = Round(CA, 4) GNH3x = 15.5 * c - 15.5 * cxx End If c2 = -6571.06 + 39.9544 * tdd - 0.00243781 * (tdd ^ 2) + 0.00000039792 * (tdd ^ 3) - 22.7722 * P2 - 0.00000004979 * (P2 ^ 2) + 8286.02 * (P2 / tdd) + 0.0172363 141 * P2 * tdd + 0.00000077344 * P2 * (tdd ^ 2) - 3.62962 * tdd * (Log(tdd)) - tdd * (Log(P2)) d2 = (-1467 * (cxx ^ 2) * (1.25 * cxx ^ - 5.58 * cxx ^ + 5.96 * cxx - 0.42)) + ((17.3 * (cxx ^ 3) - 13.8 * (cxx ^ 4)) * tdd) 'TINH ENTANPI CUA NH3 inh3x = Round(-1058.43 + 3.02148 * tdd + 0.00312072 * tdd ^ + 0.15536 * P2 0.0000123079 * P2 ^ - 0.0000012946 * P2 * tdd ^ 2, 2) 'TINH ENTANPI CUA NUOC ih2ox = Round(-1039.5 + 3.56393 * tdd + 0.000919594 * tdd ^ + 0.124376 * P2 - 0.00000222979 * P2 ^ - 0.00000035696 * P2 * tdd ^ 2, 2) 'TINH qt qttx = Round(-839.87 * (1 - 0.071 * cxx) * ((1 - cxx) * cxx + (1 - cxx ^ 2) * cxx ^ 2), 2) 'ENTANPI CUA DUNG DICH LONG i3 = Round((cxx * inh3x + (1 - cxx) * ih2ox + qttx), 3) 'TINH NONG DO CUA HOI cvx = - (1 - cxx) * Exp((c2 + d2) / tdd) 'TINH IVNH3 ivnh3x = Round(770.761 + 1.86947 * tdd + 0.0000578293 * tdd ^ + 0.000000731509 * tdd ^ + 8.98074 * P2 - 4580.15 * P2 / tdd, 2) 'TINH IVH2O ivh2ox = Round(1993.19 + 1.88878 * tdd - 0.000205512 * tdd ^ + 0.000000367295 * tdd ^ + 10.6342 * P2 - 7648.34 * P2 / tdd, 2) 'ENTANPI CUA HOI i4 = Round((cvx * ivnh3x + (1 - cvx) * ivh2ox), 3) Q2 = (15.5 - GNH3x) * i3 + GNH3x * i4 Qu = (Q2 - Q1) + TT * (P - P2) * 100 If (Qu + Ql + Qs - Q > 0) And (Qu + Ql + Qs - Q < 1) Then Exit For End If Next tdd Text58 = Round(qt, 2) Text59 = Round(Qb, 2) Text60 = Round(Qe, 2) Text61 = Round(Qs, 2) Text16 = Round(Qu, 2) Text14 = Round(tp, 2) Text15 = Round((tdd - 273), 2) 142 Text62 = Round(Ql, 2) Text17 = cxx * 100 Text18 = Round(cvx * 100, 2) Text19 = P2 Text20 = Round(i3, 2) Text21 = Round(i4, 2) Text22 = Round(GNH3x, 3) Text23 = Round(Q, 2) Text24 = Round(((Qu + Qs) / Q), 2) Text77 = Text16 Data15.Recordset.AddNew Data15.Refresh Text81 = Text23 Data16.Recordset.AddNew Data16.Refresh 'Text3 = 15.5 - Text22 'Text2 = Text17 Text63 = Text14 Text80 = tdd - 273 Text64 = Text64 + Data13.Recordset.AddNew Data13.Refresh Data14.RecordSource = "select*from bangd,bange where bangd.entan=bange.entr" Data14.Refresh Text1 = Text65 Text4 = Text66 Text69 = Text74 Text70 = Text75 Text78 = Q1 Text79 = Q2 Dim tam As Double tam = Text77 Text76 = Text76 + tam Dim tam1 As Double tam1 = Text81 Text82 = Text82 + tam1 Text83 = tp1 143 Text84 = tk Text85 = Text14 Text86 = Text15 Text87 = Text17 Text88 = Text19 Data17.Recordset.AddNew End Sub Private Sub Command3_Click() Unload Me End Sub Private Sub Form_Load() Text64 = Data14.RecordSource = "select*from bangd,bange where bange.entr=1" Data14.Refresh Text67 = Text72 Text71 = Text73 Text69 = Text74 Text70 = Text75 Text1 = Text65 Text2 = 48 Text3 = 15.5 Text4 = Text66 Text5 = 0.08 Text6 = 32 Text7 = 15 Text8 = 0.04 Text9 = 2.25 Text10 = Text11 = 0.95 Text12 = 0.88 Text25 = 0.1 Text26 = 1.5 Text27 = 1.5 Text43 = 51.36 Text44 = 0.46 Text63 = Text4 Text68 = 10 Text80 = Text63 144 Text76 = Text82 = End Sub Code chương trình tính toán entanpi Private Sub Command1_Click() Dim X As Variant X = InputBox("NHAP VAO NONG DO CUA DUNG DICH") Text1 = X Dim Y As Variant Y = InputBox("NHAP VAO NHIET DO CUA DUNG DICH") Text2 = Y Dim Z As Variant Z = InputBox("NHAP VAO AP SUAT CUA DUNG DICH") Text3 = Z End Sub Private Sub Command2_Click() If Text1 "" And Text2 "" And Text3 "" Then If Text1 >= And Text1 And Text3 > Then Dim c As Double c = Text1 / 100 Dim T As Double T = Text2 Dim P As Double P = Text3 P = Text3 Dim cs As Double 'TINH SO C1 c1 = -6571.06 + 39.9544 * T - 0.00243781 * (T ^ 2) + 0.00000039792 * (T ^ 3) 22.7722 * P - 0.00000004979 * (P ^ 2) + 8286.02 * (P / T) + 0.0172363 * P * T + 0.00000077344 * P * (T ^ 2) - 3.62962 * T * (Log(T)) - T * (Log(P)) Dim d As Double d = (-1467 * (c ^ 2) * (1.25 * c ^ - 5.58 * c ^ + 5.96 * c - 0.42)) + ((17.3 * (c ^ 3) - 13.8 * (c ^ 4)) * T) 'TINH ENTANPI CUA NH3 Dim inh3 As Double inh3 = Round(-1058.43 + 3.02148 * T + 0.00312072 * T ^ + 0.15536 * P 0.0000123079 * P ^ - 0.0000012946 * P * T ^ 2, 2) 'TINH ENTANPI CUA NUOC 145 Dim ih2o As Double ih2o = Round(-1039.5 + 3.56393 * T + 0.000919594 * T ^ + 0.124376 * P 0.00000222979 * P ^ - 0.00000035696 * P * T ^ 2, 2) 'TINH qt Dim qt As Double qt = Round(-839.87 * (1 - 0.071 * c) * ((1 - c) * c + (1 - c ^ 2) * c ^ 2), 2) 'ENTANPI CUA DUNG DICH LONG Text5 = Round((c * inh3 + (1 - c) * ih2o + qt), 3) 'TINH NONG DO CUA HOI Dim cv As Double cv = - (1 - c) * Exp((c1 + d) / T) Text4 = Round(cv * 100, 3) 'TINH IVNH3 Dim ivnh3 As Double ivnh3 = Round(770.761 + 1.86947 * T + 0.0000578293 * T ^ + 0.000000731509 * T ^ + 8.98074 * P - 4580.15 * P / T, 2) 'TINH IVH2O Dim ivh2o As Double ivh2o = Round(1993.19 + 1.88878 * T - 0.000205512 * T ^ + 0.000000367295 * T ^ + 10.6342 * P - 7648.34 * P / T, 2) 'ENTANPI CUA HOI Text6 = Round((cv * ivnh3 + (1 - cv) * ivh2o), 3) 'TINH NHIET DUNG RIENG CUA DUNG DICH Text7 = Round(3.56393 + 1.83918 * T / 1000 + 7.1392 * P * T / 10000000 (0.54245 - 4.40232 * T / 1000 + 1.87528 * P * T / 1000000) * c, 3) 'TINH HE SO DAN NHIET CUA DUNG DICH Dim LAMDANH3 As Double Dim LAMDAH2O As Double LAMDANH3 = 0.528 - 1.669 * (T - 273) / 1000 - 6.2 * (T - 273) ^ / 1000000 LAMDAH2O = 0.562 + 1.893 * (T - 273) / 1000 - 7.11 * (T - 273) ^ / 1000000 Text8 = Round(c * LAMDANH3 + (1 - c) * LAMDAH2O, 3) End If End If End Sub Private Sub Command3_Click() Text1 = "" Text2 = "" Text3 = "" 146 Text4 = "" Text5 = "" Text6 = "" Text7 = "" Text8 = "" Command2.Enabled = True End Sub Chương trình tính thông số vật lý dung dòch 'TINH NHIET DO SOI 'NHAP Private Sub Command1_Click() Dim X As Variant X = InputBox("NHAP VAO NONG DO CUA DUNG DICH") Text1 = X Dim Y As Variant Y = InputBox("NHAP VAO AP SUAT CUA DUNG DICH") Text2 = Y End Sub 'XUAT Private Sub Command2_Click() Dim a As Double Dim B As Double If Text1 "" And Text2 "" Then If Text1 >= And Text1 = Then 'NONG DO Dim c As Double c = Text1 / 100 'AP SUAT Dim P As Double P = Text2 'TINH SO a VA b a = -2103.5 + 4669.96 * c - 20228.3 * (c ^ 2) + 56507 * (c ^ 3) - 80989.9 * (c ^ 4) + 55286.5 * (c ^ 5) - 14361.4 * (c ^ 6) B = 5.65208 - 7.0317 * c + 37.9018 * (c ^ 2) - 102.912 * (c ^ 3) + 135.789 * (c ^ 4) - 82.7106 * (c ^ 5) + 18.4113 * (c ^ 6) 'NHIET DO SOI T Text3 = Round(a / ((Log(P) / Log(10)) + 0.00847711 - B), 3) Command1.Enabled = False 147 End If If Text1 < Or Text1 > 100 Or Text2 < Then MsgBox ("YOUR INPUT DATA-ERROR") End If End If End Sub 'XOA Private Sub Command3_Click() Text1 = "" Text2 = "" Text3 = "" Command1.Enabled = True End Sub 'TINH AP SUAT 'NHAP Private Sub Command4_Click() Dim X As Variant X1 = InputBox("NHAP VAO NONG DO CUA DUNG DICH") Text4 = X1 Dim Y As Variant Y1 = InputBox("NHAP VAO NHIET DO CUA DUNG DICH") Text5 = Y1 End Sub 'XUAT Private Sub Command5_Click() Dim A1 As Double Dim B1 As Double If Text4 "" And Text5 "" Then If Text4 >= And Text4 = Then Dim c1 As Double c1 = Text4 / 100 Dim T As Double T = Text5 'TINH THEO CONG THUC 'If Option1 = True Then A1 = -2103.5 + 4669.96 * c1 - 20228.3 * (c1 ^ 2) + 56507 * (c1 ^ 3) - 80989.9 * (c1 ^ 4) + 55286.5 * (c1 ^ 5) - 14361.4 * (c1 ^ 6) 148 B1 = 5.65208 - 7.0317 * c1 + 37.9018 * (c1 ^ 2) - 102.912 * (c1 ^ 3) + 135.789 * (c1 ^ 4) - 82.7106 * (c1 ^ 5) + 18.4113 * (c1 ^ 6) Text6 = Round(10 ^ ((A1 / T) + B1 - 0.00847711), 3) 'End If 'TINH THEO CONG THUC 'If Option2 = True Then 'A1 = 7.44 - 1.767 * c1 + 0.9823 * c1 ^ + 0.3627 * c1 ^ 'B1 = 2013.8 - 2155.7 * c1 + 1540.9 * c1 ^ - 194.7 * c1 ^ 'Text6 = Round((10 ^ (A1 - (B1 / T))) / 100, 3) 'End If Command4.Enabled = False End If If Text4 < Or Text4 > 100 Or Text5 < Then MsgBox ("YOUR INPUT DATA-ERROR") End If End If End Sub 'XOA Private Sub Command6_Click() Text4 = "" Text5 = "" Text6 = "" Command4.Enabled = True End Sub 'TINH NONG DO DUNG DICH 'NHAP Private Sub Command7_Click() Dim x2 As Variant Dim y2 As Variant x2 = InputBox("NHAP VAO AP SUAT CUA DUNG DICH") y2 = InputBox("NHAP VAO NHIET DO CUA DUNG DICH") Text7 = x2 Text8 = y2 End Sub 'XUAT Private Sub Command8_Click() Dim a2 As Double Dim b2 As Double 149 If Text7 "" And Text8 "" Then If Text7 >= And Text8 >= Then Dim T2 As Double T2 = Text8 Dim P2 As Double P2 = Text7 Dim c2 As Double Dim CC As Double CC = For c2 = 0.01 To Step 0.00001 'CC = CC + C2 a2 = -2103.5 + 4669.96 * c2 - 20228.3 * (c2 ^ 2) + 56507 * (c2 ^ 3) - 80989.9 * (c2 ^ 4) + 55286.5 * (c2 ^ 5) - 14361.4 * (c2 ^ 6) b2 = 5.65208 - 7.0317 * c2 + 37.9018 * (c2 ^ 2) - 102.912 * (c2 ^ 3) + 135.789 * (c2 ^ 4) - 82.7106 * (c2 ^ 5) + 18.4113 * (c2 ^ 6) CC = a2 / ((Log(P2) / Log(10)) + 0.00847711 - b2) If (CC - T2) > And (CC - T2) < 0.01 Then Exit For End If Next c2 Text9 = Round(c2 * 100, 2) Command7.Enabled = False End If If Text7 < Or Text8 < Then MsgBox ("YOUR INPUT DATA-ERROR") End If End If End Sub Private Sub Command9_Click() Text7 = "" Text8 = "" Text9 = "" Command7.Enabled = True End Sub ... thụ, pin mặt trời, … Trong luận văn tác giả nghiên cứu thực nghiệm collector mặt trời cấp nhiệt cho máy lạnh hấp thụ loại gián đoạn để sản xuất nước đá 8 Chương 1: TỔNG QUAN VỀ MÁY LẠNH HẤP THỤ... hấp thụ NH3- H2O loại gián cấp nhiệt lượng mặt trời để sản xuất nước đá thực hiện: - Hệ thống máy lạnh hấp thụ chạy lượng mặt trời sản xuất nước đá để trữ thực phẩm cá làng Maruata bang Michoacan,... thu loại phẳng dùng để cấp nhiệt cho máy lạnh hấp thụ • Xây dựng thực nghiệm mô hình máy lạnh hấp thụ loại gián đoạn • Xây dựng chương trình mô cho máy lạnh hấp thụ loại gián đoạn So sánh với thực