Nghiên cứu tính toán thiết kế máy lạnh hấp thụ nh3 h2o loại gián đọan cấp nhiệt bằng phụ phẩm nông nghiệp để sản xuất nước đá báo cáo tổng kết kết quả đề tài khcn cấp trường msđt t ck 2012 04

ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA FOG BÁO CÁO TỔNG KẾT KẾT QUẢ ĐỀ TÀI KHCN CẤP TRƯỜNG Tên đề tài: NGHIÊN CỨU TÍNH TỐN THIẾT KẾ MÁY LẠNH HẤP THỤ NH3-H2O LOẠI GIÁN ĐOẠN CẤP NHIỆT BẰNG PHỤ PHẨM NÔNG NGHIỆP ĐỂ SẢN XUẤT NƯỚC ĐÁ Mã số đề tài:T-CK-2012-04 Thời gian thực hiện: 12 tháng Chủ nhiệm đề tài: VÕ KIẾN QUỐC Cán tham gia đề tài: Thành phố Hồ Chí Minh – Tháng 9/2014 Danh sách cán tham gia thực đề tài (Ghi rõ học hàm, học vị, đơn vị công tác gồm môn, Khoa/Trung tâm) ThS Võ Kiến Quốc M c l c  1.  Đặt vấn đề 4  2.  Tổng quát máy hấp thụ NH3-H2O 5  2.1  Sự khác biệt máy lạnh hấp thụ máy lạnh có máy nén 5  2.2  Dung dịch làm việc 6  2.3  Ưu điểm MLHT 7  2.4  Nhược điểm MLHT: 9  2.5  Nguyên lý làm việc máy lạnh hấp thụ NH3-H2O 10  3.  Các kết nghiên cứu 11  1  Thống kê lượng phụ phẩm nông nghiệp khu vực đồng sông Cửu Long: 11  2  Xác định nhiệt trị trấu 14  3  Lập sơ đồ phương án phân tích lựa chọn phương án thiết kế 15  4.  Tính tốn chu trình thiết bị lạnh: 19  5.  Chương trình mơ 24  5.1  sinh Ảnh hưởng nhiệt độ dung dịch đến áp suất nồng độ dung dịch bình phát 25  5.2  Ảnh hưởng nhiệt độ dung dịch đến nồng độ NH3 bình phát sinh 26  5.3  đặc Ảnh hưởng áp suất dung dịch nhiệt độ giải nhiệt đến nồng độ dung dịch đậm 27  5.4  Ảnh hưởng nhiệt độ dung dịch nhiệt độ ngưng tụ đến nồng độ dung dịch loãng 27  5.5  sinh Ảnh hưởng lượng tiêu hao nhiên liệu đến nhiệt độ dung dịch bình phát 28  6.  Các kết bật 29  7.  Kết luận kiến nghị 30  Phụ lục 32  Phụ lục 32  Phụ lục 34  Đặt vấn đề Có thể nói máy lạnh có máy nén độc chiếm thị trường máy lạnh Việt Nam số nước khác Bên cạnh ưu điểm bật máy lạnh có máy nén tồn nhược điểm lớn tiêu tốn nhiều điện năng, làm việc với điện sử dụng môi chất lạnh không thân thiện với môi trường Ngày khoa học phát triển cho đời môi chất lạnh có mức độ ảnh hưởng đến mơi trường khơng nhiều Tuy nhiên Việt Nam việc ứng dụng mơi chất vào thực tiễn cịn hạn chế Máy lạnh hấp thụ đời từ lâu, nhiên ưu điểm máy lạnh bị lãng quên thời gian dài vấn đề mơi trường biến đổi khí hậu đề cao máy lạnh hấp thụ ý trở lại Hiện máy lạnh hấp thụ thay vai trị máy lạnh có máy nén nước phát triển như: Mỹ, Nhật Bản, Trung Quốc, Nga,… Ưu điểm máy lạnh hấp thụ không làm việc điện mà làm việc với nhiệt sử dụng môi chất lạnh thân thiện với môi trường Nguồn nhiệt sử dụng cho máy lạnh hấp thụ nguồn nhiệt thải cơng nghiệp, nhiên liệu hóa thạch, phụ phẩm nơng nghiệp lượng mặt trời Nước đá sản phẩm thiếu đời sống đại, toàn nước đá nước ta sản xuất máy lạnh có máy nén Thời gian làm đá 50kg 20 giờ/mẻ, với thời gian hệ thống máy lạnh có máy nén tiêu thụ lượng điện khơng nhỏ Do khu vực khơng có lưới điện việc sản xuất nước đá máy lạnh có máy nén khó khăn Làm lạnh máy lạnh hấp thụ phương pháp biết đến từ lâu Tuy nhiên ưu điểm máy lạnh có máy nén làm cho người quên lãng kỹ thuật này, người thấy tác hại đến môi trường máy lạnh có máy nén máy lạnh hấp thụ lựa chọn ưu tiên nước phát triển Ở nước ta máy lạnh hấp thụ cịn xa lạ Tuy nhiên năm gần máy lạnh hấp thụ quan tâm trở lại Năm 2004 Trần Thanh Kỳ [1] nghiên cứu thiết kế chế tạo thành công hệ thống máy lạnh hấp thụ NH3-H2O để sản xuất nước đá có suất 3,2 đá/ngày Nhiên liệu sử dụng loại nhiên liệu rẻ tiền như: than cám, trấu mùn cưa Kết nghiên cứu cho thấy giá thành sản xuất 1kg đá dao động khoảng từ 86,1-95,8 đồng Năm 2006 Võ Kiến Quốc Lê Chí Hiệp [2] nghiên cứu lý thuyết thực nghiệm máy lạnh hấp thụ NH3-H2O loại gián đoạn cấp nhiệt collector mặt trời dạng ống cánh để sản xuất nước đá Kết nghiên cứu cho thấy suất trung bình hệ thống khoảng từ 1,6-2,5kg đá/m2 collector/ngày Nồng độ dung dịch đậm đặc tối ưu khoảng từ 46% đến 48% nồng độ dung dịch loãng từ 36% đến 39% Trong thời đại vấn đề môi trường quan tâm hàng đầu việc nghiên cứu ứng dụng máy lạnh hấp thụ nước ta việc làm cần thiết Việc ứng dụng máy lạnh hấp thụ giải vấn đề thiếu điện vấn đề ô nhiễm môi trường mà nước ta phải đối mặt Việt Nam nước có nơng nghiệp phát triển lượng phụ phẩm nơng nghiệp lớn, điển hình trấu Tuy nhiên lượng phụ phẩm nông nghiệp lại thừa số lượng lớn vừa gây lãng phí vừa ảnh hưởng đến mơi trường Với mục đích bảo vệ mơi trường đẩy mạnh việc sử dụng phụ phẩm nông nghiệp, nghiên cứu tác giả trình bày phương pháp nghiên cứu lý thuyết hệ thống máy lạnh hấp thụ NH3-H2O loại gián đoạn sử dụng trấu để sản xuất nước đá cho tỉnh Đồng Bằng Sông Cửu Long Tổng quát máy hấp thụ NH3-H2O 2.1 Sự khác biệt máy lạnh hấp thụ máy lạnh có máy nén So với máy lạnh có máy nén máy lạnh hấp thụ có khác biệt lớn sau đây: - Năng lượng cấp vào cho máy lạnh hấp thụ nhiệt năng, lượng cấp vào cho máy lạnh có máy nén - Môi chất làm việc máy lạnh hấp thụ dung dịch chất khiết có nhiệt độ sôi cách biệt lớn điều kiện áp suất Mơi chất làm việc máy lạnh có máy nén thường chất khiết hỗn hợp đồng sơi.Mơi chất làm việc chu trình máy lạnh hấp thụ dung dịch, chất có nhiệt độ bay thấp gọi tác nhân lạnh chất có nhiệt độ bay cao gọi chất hấp thụ Trong thực tế cặp môi chất sử dụng rộng rãi NH3H2O LiBr-H2O 2.2 Dung dịch làm việc Cặp môi chất NH3-H2O sử dụng để làm lạnh đến nhiệt độ âm dương cặp môi chất thường dùng kỹ thuật đông lạnh Đối với cặp mơi chất LiBr-H2O làm lạnh nhiệt độ dương nên cặp môi chất sử dụng kỹ thuật điều hịa khơng khí Hiệu suất hệ thống MLHT phụ thuộc nhiều vào tính chất hóa học nhiệt động cặp dung dịch làm việc bên Dung dịch làm việc MLHT hỗn hợp chất khiết khác nhau, có nhiệt độ sơi cách biệt điều kiện áp suất khơng phản ứng hóa học với nhau, không độc hại không gây nổ Đối với tác nhân lạnh phải có ẩn nhiệt hóa lớn Theo định lý Carnot hiệu suất chu trình thuận nghịch khơng phụ thuộc vào chất mơi giới làm việc chu trình mà phụ thuộc vào nhiệt độ nguồn nóng nguồn lạnh Tuy nhiên chu trình thực tế khơng phải thuận nghịch hiệu chu trình MLHT phụ thuộc nhiều vào chất mơi giới chu trình Có nhiều cặp tác nhân sử dụng cho MLHT Theo kết nhà nghiên cứu trước có đến 40 loại tác nhân lạnh 200 loại chất hấp thụ sử dụng MLHT Tuy nhiên thực tế hai cặp dung dịch NH3-H2O H2O-LiBr sử dụng phổ biến Trong cặp dung dịch làm việc chất đứng trước môi chất lạnh, chất đứng sau chất hấp thụ Do tính chất nhiệt động tác nhân lạnh người ta dùng cặp dung dịch NH3-H2O cho ứng dụng làm lạnh sâu đến 77oC cịn cặp dung dịch H2O-LiBr dùng để làm lạnh kỹ thuật điều hòa khơng khí Dung dịch làm việc MLHT chứng minh hồn tồn khơng gây vấn đề ảnh hưởng đến mơi trường Điểm đặc biệt cặp dung dịch H2OLiBr áp suất làm việc thấp Do H2O đóng vai trị tác nhân lạnh nên để nước bay nhiệt độ thấp áp suất làm việc phải thấp Đối cặp dung dịch NH3-H2O áp suất làm việc áp suất khí cần phải ý đến vấn đề rò rỉ tác nhân NH3 môi trường Một vấn đề khác cần phải quan tâm cặp môi chất NH3-H2O NH3 bay thường có lẫn theo nước Nhiệt độ làm việc cao nồng độ NH3 bay giảm, cần phải có thêm phận tinh luyện NH3 bay để không làm giảm suất lạnh nước gây Và cuối tính chất ăn mịn tác nhân NH3 đồng kim loại chế tạo phải sử dụng thép Đối với cặp dung dịch H2O-LiBr LiBr khó bị theo nước nhiệt độ sôi LiBr nước xa điều kiện áp suất So với cặp dung dịch H2O-LiBr cặp dung dịch NH3-H2O có nhược điểm sau đây: - Hệ số COP cặp dung dịch NH3-H2O thấp cặp dung dịch H2O-LiBr - Nhiệt độ nguồn nhiệt cấp cho bình phát sinh cao cặp dung dịch H2O-LiBr - Áp suất ngưng tụ cao tiêu tốn cơng bơm lớn - Hơi nước dễ bị bay theo NH3 khỏi bình phát sinh - Nguy hiểm NH3 gây ngộ độc 2.3 Ưu điểm MLHT - Phương pháp nén tác nhân lạnh: Trong máy lạnh có máy nén tác nhân lạnh nén từ áp suất thấp lên áp suất cao máy nén Ngoài máy nén hệ thống lạnh có máy nén cịn có nhiệm vụ tạo tuần hoàn tác nhân lạnh hệ thống Máy nén làm việc với tác nhân lạnh trạng thái mà tích riêng lớn công tiêu hao cấp vào cho máy nén lớn Một vấn đề khác kích thước máy nén lớn công suất tăng điện tiêu hao cho máy nén tăng theo Trong MLHT không cần dùng đến máy nén mà chức máy nén thay bình phát sinh, bình hấp thụ bơm Tác nhân lạnh sau bay chất hấp thụ để trở thành trạng thái lỏng để làm giảm thể tích riêng cơng cấp vào cho bơm nhỏ không thay đổi nhiều suất tăng - Dạng lượng sử dụng: Năng lượng sử dụng MLHT chủ yếu lượng nhiệt Đây ưu điểm bật MLHT so với máy lạnh có máy nén Năng lượng cấp vào cho MLHT nhiệt thải từ quy trình cơng nghệ cơng nghiệp, nước từ lị hơi, lượng mặt trời, khí đốt, nhiên liệu hóa thạch, phụ phẩm nông nghiệp,… Năng lượng sử dụng lượng nhiệt góp phần giảm tải cho hệ thống điện sử dụng khu vực khơng có điện lưới quốc gia - Chi phí vận hành: Máy lạnh có máy nén vận hành với điện phụ thuộc nhiều vào điện Hiện giá điện đắt đỏ chi phí vận hành cho máy lạnh có máy nén lớn Trong MLHT sử dụng bơm dung dịch nhỏ phí vận hành cho bơm nhỏ - Tiếng ồn rung: Máy nén máy lạnh có máy nén hoạt động tốc độ cao cần bơi trơn gây tiếng ồn MLHT khơng có phận chuyển động hoạt động êm không gây rung lắc móng khơng địi hỏi gia cố chống rung - Bảo dưỡng: máy nén phận máy lạnh có máy nén có nhiều chi tiết chuyển động Các chi tiết phải bôi trơn kiểm tra thường xuyên để đảm bảo cho máy nén vận hành tốt Máy nén cần thay chi tiết bên hoạt động lâu ngày bị mòn làm hiệu suất nén giảm xuống Trong có số chi tiết đắt tiền suppap hút đẩy môtơ điện Máy nén cần làm mát cần phải có thêm bơm nước để làm mát cho máy nén Trong MLHT chi tiết chuyển động bơm dung dịch gặp hư hỏng cơng tác bảo trì cho MLHT đơn giản - Điều khiển suất: máy lạnh có máy nén suất hệ thống điều chỉnh theo cấp chi phí lượng tiêu hao cho máy nén tỷ lệ theo cấp điều chỉnh Trong MLHT việc điều chỉnh suất vô cấp chi phí tiêu hao lượng cho hệ thống với công suất hoạt động - Tác nhân lạnh chi phí: cặp tác nhân lạnh NH3-H2O NH3 tác nhân lạnh, cặp tác nhân lạnh H2O-LiBr H2O tác nhân lạnh Cả NH3 nước chất rẻ tiền va dễ tìm kiếm Trong máy lạnh có máy nén sử dụng tác nhân lạnh CFC HCFC đắt tiền - Rò rỉ tác nhân: MLHT khơng có xảy rị rỉ tác nhân lạnh có rị rỉ chi phí nạp lại khơng đáng kể Trong máy lạnh có máy nén dễ xảy rị rỉ chi phí nạp lại lớn - Hiệu ứng nhà kính : Hầu hết chất làm lạnh sử dụng hệ thống máy lạnh có máy nén gây hiệu ứng nhà kính Theo Nghị định thư Montreal, việc sử dụng tác nhân lạnh gây hiệu ứng nhà kính phải dừng lại hồn tồn vào năm 2020 Trong hệ thống MLHT không sử dụng tác tạo hiệu ứng nhà kính , việc sử dụng tác nhân không bị cản trở tương lai 2.4 Nhược điểm MLHT: - Kích thước trọng lượng: MLHT có kích thước trọng lượng lớn nhiều so với hệ thống máy lạnh có máy nén cơng suất - Hiệu suất thấp hơn: hệ số làm lạnh MLHT thấp nhiều so với máy lạnh có máy nén Hệ số làm lạnh MLHT double effect H2O-LiBr có hệ số làm lạnh tối đa khoảng 1,1 máy lạnh có máy nén có hệ số làm lạnh khoảng từ 4-5 so sánh điều kiện làm việc Tuy nhiên nhìn vào hệ số làm lạnh khơng thể kết luận MLHT làm việc không hiệu máy lạnh có máy nén Bởi dạng lượng sử dụng khác - Giá thành cao hơn: MLHT có nhiều thiết bị trao đổi nhiệt kích thước thiết bị lớn phí sản xuất cho thiết bị cao - Khó vận chuyển hơn: kích thước hệ thống MLHT lớn việc vận chuyển lắp đặt khó khăn nhiều khơng thể lắp đặt khu vực cao - Có thể bị kết tinh - Nhả nhiệt nhiều hơn: MLHT ngồi việc giải nhiệt cho bình ngưng cịn phải giải nhiệt cho bình hấp thụ, tháp tinh luyện,… tính tồn hệ thống hệ số nhả nhiệt MLHT khoảng 2,5 Trong máy lạnh có máy nén giải nhiệt cho bình ngưng tụ cho máy nén nên hệ số nhả nhiệt khoảng 1,2 Chính lý làm cho giá thành MLHT cao máy lạnh có máy nén tháp giải nhiệt lớn - Ăn mòn hệ thống làm việc nhiệt độ cao - Áp suất làm việc thấp: áp suất làm việc hệ thống MLHT H2O-LiBr thấp xảy tượng khơng khí bên ngồi xâm nhập vào hệ thống 2.5 Nguyên lý làm việc máy lạnh hấp thụ NH3-H2O Sơ đồ nguyên lý máy lạnh hấp thụ NH3-H2O trình bày hình đây: A B D C Hình 1: Sơ đồ nguyên lý máy lạnh hấp thụ NH3-H2O Gọi cs nồng độ dung dịch đậm đặc bình phát sinh A, cấp nhiệt cho dung dịch chứa bình phát sinh dung dịch sơi bay Hơi NH3 bay làm cho nồng độ dung dịch chứa bình phát sinh giảm xuống ta gọi nồng độ dung dịch loãng cw Hơi bay NH3 nhiên thường có lượng nhỏ nước bay theo nồng độ bay thường nhỏ 100% Để tách lượng nước khỏi NH3 người ta cho hỗn hợp qua thiết bị có tên tháp chưng cất Hơi NH3 sau khỏi tháp chưng cất có nồng độ cao, xem tinh khiết vào bình ngưng B để thực trình ngưng tụ Sau khỏi bình ngưng B, lỏng NH3 đưa qua van tiết lưu để vào bình bay C Hơi NH3 bay từ thiết bị bay C đưa bình hấp thụ D hấp thụ dung dịch lỗng có nồng độ cw đến từ bình phát sinh A Quá trình hấp thụ phát sinh nhiệt lượng NH3 nhả nhiệt ẩn ngưng tụ bình hấp thụ D phải 10 Dim viNH3max As Double 'Nhiet dung dich cuoi giai doan phat sinh Dim Tmax As Double Tmax = Text2 + 273.15 'Ap suat ngung tu Dim Pk As Double Pk = Text8 liNH3max = Round(-1058.43 + 3.02148 * T1 + 0.00312072 * T1 ^ + 0.15536 * Pk 0.0000123079 * Pk ^ - 0.0000012946 * Pk * T1 ^ 2, 2) viNH3max = Round(770.761 + 1.86947 * Tmax + 0.0000578293 * Tmax ^ + 0.000000731509 * Tmax ^ + 8.98074 * Pk - 4580.15 * Pk / Tmax, 2) Text12 = Round((Text11 * (viNH3max - liNH3max)), 1) 'Tinh nong 'dam dac Dim a3 As Double Dim b3 As Double If Text4 "" And Text9 "" Then If Text4 >= And Text9 >= Then Dim T3 As Double T3 = Text4 + 273.15 Dim P3 As Double P3 = Text9 Dim c3 As Double Dim CCC As Double CCC = For c3 = 0.01 To Step 0.00001 'CC = CC + C2 a3 = -2103.5 + 4669.96 * c3 - 20228.3 * (c3 ^ 2) + 56507 * (c3 ^ 3) - 80989.9 * (c3 ^ 4) + 55286.5 * (c3 ^ 5) - 14361.4 * (c3 ^ 6) b3 = 5.65208 - 7.0317 * c3 + 37.9018 * (c3 ^ 2) - 102.912 * (c3 ^ 3) + 135.789 * (c3 ^ 4) - 82.7106 * (c3 ^ 5) + 18.4113 * (c3 ^ 6) CCC = a3 / ((Log(P3) / Log(10)) + 0.00847711 - b3) If (CCC - T3) > And (CCC - T3) < 0.01 Then Exit For End If Next c3 Text14 = Round(c3 * 100, 2) End If End If 'loang Dim a4 As Double Dim b4 As Double If Text2 "" And Text8 "" Then If Text2 >= And Text8 >= Then Dim T4 As Double T4 = Text2 + 273.15 36 Dim P4 As Double P4 = Text8 Dim c4 As Double Dim CCCC As Double CCCC = For c4 = 0.01 To Step 0.00001 'CC = CC + C2 a4 = -2103.5 + 4669.96 * c4 - 20228.3 * (c4 ^ 2) + 56507 * (c4 ^ 3) - 80989.9 * (c4 ^ 4) + 55286.5 * (c4 ^ 5) - 14361.4 * (c4 ^ 6) b4 = 5.65208 - 7.0317 * c4 + 37.9018 * (c4 ^ 2) - 102.912 * (c4 ^ 3) + 135.789 * (c4 ^ 4) - 82.7106 * (c4 ^ 5) + 18.4113 * (c4 ^ 6) CCCC = a4 / ((Log(P4) / Log(10)) + 0.00847711 - b4) If (CCCC - T4) > And (CCCC - T4) < 0.01 Then Exit For End If Next c4 Text15 = Round(c4 * 100, 2) End If End If 'Tinh khoi luong dung dich Dim cs As Double Dim cw As Double Dim mNH3 As Double Dim mH2O As Double Dim m As Double m = Text11 Dim cs100 As Double Dim cw100 As Double cs100 = Text14 cw100 = Text15 cs = cs100 / 100 cw = cw100 / 100 mNH3 = (m - m * cw) / (1 - cw - (1 - cs) * cw / cs) mH2O = mNH3 * (1 - cs) / cs Text13 = Round((mNH3 + mH2O), 1) 'Tinh ap suat o trang thai ban dau Dim Tx As Double Dim Px As Double Dim muy As Double Dim roNH3 As Double Dim P1 As Double muy = 17.03 Tx = Text47 + 273.15 Dim aaa As Double Dim bbb As Double 37 aaa = -2103.5 + 4669.96 * cs - 20228.3 * (cs ^ 2) + 56507 * (cs ^ 3) - 80989.9 * (cs ^ 4) + 55286.5 * (cs ^ 5) - 14361.4 * (cs ^ 6) bbb = 5.65208 - 7.0317 * cs + 37.9018 * (cs ^ 2) - 102.912 * (cs ^ 3) + 135.789 * (cs ^ 4) - 82.7106 * (cs ^ 5) + 18.4113 * (cs ^ 6) P1 = Round(10 ^ ((aaa / Tx) + bbb - 0.00847711), 3) ' Tinh nong dung dich binh phat sinh Dim cx As Double Dim Ax As Double Dim Bx As Double Dim cc As Double If P1 < Pk Then Px = P1 For cx = cw To cs Step 0.000001 Ax = -2103.5 + 4669.96 * cx - 20228.3 * (cx ^ 2) + 56507 * (cx ^ 3) - 80989.9 * (cx ^ 4) + 55286.5 * (cx ^ 5) - 14361.4 * (cx ^ 6) Bx = 5.65208 - 7.0317 * cx + 37.9018 * (cx ^ 2) - 102.912 * (cx ^ 3) + 135.789 * (cx ^ 4) - 82.7106 * (cx ^ 5) + 18.4113 * (cx ^ 6) cc = Ax / ((Log(Px) / Log(10)) + 0.00847711 - Bx) If (cc - Tx) > And (cc - Tx) < 0.01 Then Exit For End If Next cx Text49 = Round(cx * 100, 2) End If If P1 > Pk Then Px = Pk For cx = cw To cs Step 0.0001 Ax = -2103.5 + 4669.96 * cx - 20228.3 * (cx ^ 2) + 56507 * (cx ^ 3) - 80989.9 * (cx ^ 4) + 55286.5 * (cx ^ 5) - 14361.4 * (cx ^ 6) Bx = 5.65208 - 7.0317 * cx + 37.9018 * (cx ^ 2) - 102.912 * (cx ^ 3) + 135.789 * (cx ^ 4) - 82.7106 * (cx ^ 5) + 18.4113 * (cx ^ 6) cc = Ax / ((Log(Px) / Log(10)) + 0.00847711 - Bx) If (cc - Tx) > And (cc - Tx) < 0.01 Then Exit For End If Next cx Text49 = Round(cx * 100, 2) Gh1 = (cx * Text13 - mNH3) / (cx - 1) End If roNH3 = muy * Px / (0.08314 * (Tx + 273)) Dim thetichbayhoi As Double thetichbayhoi = 0.05 Dim thetichrieng As Double thetichrieng = / roNH3 'Dim Gh1 As Double Gh1 = thetichbayhoi / thetichrieng 38 Dim G1 As Double G1 = Text13 - Gh1 Text16 = Px Text50 = Round(Gh1, 2) 'tinh entanpi cua dung dich o trang thai ban dau Dim cc1 As Double cc1 = -6571.06 + 39.9544 * Tx - 0.00243781 * (Tx ^ 2) + 0.00000039792 * (Tx ^ 3) 22.7722 * Px - 0.00000004979 * (Px ^ 2) + 8286.02 * (Px / Tx) + 0.0172363 * Px * Tx + 0.00000077344 * Px * (Tx ^ 2) - 3.62962 * Tx * (Log(Tx)) - Tx * (Log(Px)) Dim dd1 As Double dd1 = (-1467 * (cx ^ 2) * (1.25 * cx ^ - 5.58 * cx ^ + 5.96 * cx - 0.42)) + ((17.3 * (cx ^ 3) - 13.8 * (cx ^ 4)) * Tx) Dim linh31 As Double Dim vinh31 As Double Dim liH2O1 As Double Dim viH2O1 As Double Dim qt1 As Double linh31 = Round(-1058.43 + 3.02148 * Tx + 0.00312072 * Tx ^ + 0.15536 * Px 0.0000123079 * Px ^ - 0.0000012946 * Px * Tx ^ 2, 2) vinh31 = Round(770.761 + 1.86947 * Tx + 0.0000578293 * Tx ^ + 0.000000731509 * Tx ^ + 8.98074 * Px - 4580.15 * Px / Tx, 2) liH2O1 = Round(-1039.5 + 3.56393 * Tx + 0.000919594 * Tx ^ + 0.124376 * Px 0.00000222979 * Px ^ - 0.00000035696 * Px * Tx ^ 2, 2) qt1 = Round(-839.87 * (1 - 0.071 * cx) * ((1 - cx) * cx + (1 - cx ^ 2) * cx ^ 2), 2) Dim ldd1 As Double ldd1 = Round((cx * linh31 + (1 - cx) * liH2O1 + qt1), 3) Text17 = ldd1 viH2O1 = Round(1993.19 + 1.88878 * Tx - 0.000205512 * Tx ^ + 0.000000367295 * Tx ^ + 10.6342 * Px - 7648.34 * Px / Tx, 2) Dim cv1 As Double cv1 = - (1 - cx) * Exp((cc1 + dd1) / Tx) Dim vdd1 As Double vdd1 = Round((cv1 * vinh31 + (1 - cv1) * viH2O1), 3) Text18 = vdd1 'Tinh khoi luong rieng cua hoi NH3 o trang thai Dim Tx2 As Double Dim Px2 As Double Dim cx2 As Double Dim ax2 As Double Dim bx2 As Double Dim roNH32 As Double Px2 = Text8 For cx2 = 0.1 To Step 0.00001 ax2 = -2103.5 + 4669.96 * cx2 - 20228.3 * (cx2 ^ 2) + 56507 * (cx2 ^ 3) - 80989.9 * (cx2 ^ 4) + 55286.5 * (cx2 ^ 5) - 14361.4 * (cx2 ^ 6) 39 bx2 = 5.65208 - 7.0317 * cx2 + 37.9018 * (cx2 ^ 2) - 102.912 * (cx2 ^ 3) + 135.789 * (cx2 ^ 4) - 82.7106 * (cx2 ^ 5) + 18.4113 * (cx2 ^ 6) Tx2 = Round(ax2 / ((Log(Px2) / Log(10)) + 0.00847711 - bx2), 3) roNH32 = muy * Px2 / (0.08314 * (Tx2 + 273)) Dim v2 As Double v2 = / roNH32 Dim Gh2 As Double Gh2 = thetichbayhoi / v2 Dim G2 As Double G2 = Text13 - Gh2 Dim nongdo2 As Double nongdo2 = (mNH3 - Gh2) / (Text13 - Gh2) If cx2 - nongdo2 > Then Exit For End If Next cx2 Text19 = Tx2 Text22 = Round(cx2 * 100, 3) 'Tinh entanpi cua dung dich o trang thai Dim cc2 As Double cc2 = -6571.06 + 39.9544 * Tx2 - 0.00243781 * (Tx2 ^ 2) + 0.00000039792 * (Tx2 ^ 3) - 22.7722 * Px2 - 0.00000004979 * (Px2 ^ 2) + 8286.02 * (Px2 / Tx2) + 0.0172363 * Px2 * Tx2 + 0.00000077344 * Px2 * (Tx2 ^ 2) - 3.62962 * Tx2 * (Log(Tx2)) - Tx2 * (Log(Px2)) Dim dd2 As Double dd2 = (-1467 * (cx2 ^ 2) * (1.25 * cx2 ^ - 5.58 * cx2 ^ + 5.96 * cx2 - 0.42)) + ((17.3 * (cx2 ^ 3) - 13.8 * (cx2 ^ 4)) * Tx2) Dim linh32 As Double Dim vinh32 As Double Dim liH2O2 As Double Dim viH2O2 As Double Dim qt2 As Double linh32 = Round(-1058.43 + 3.02148 * Tx2 + 0.00312072 * Tx2 ^ + 0.15536 * Px2 0.0000123079 * Px2 ^ - 0.0000012946 * Px2 * Tx2 ^ 2, 2) vinh32 = Round(770.761 + 1.86947 * Tx2 + 0.0000578293 * Tx2 ^ + 0.000000731509 * Tx2 ^ + 8.98074 * Px2 - 4580.15 * Px2 / Tx2, 2) liH2O2 = Round(-1039.5 + 3.56393 * Tx2 + 0.000919594 * Tx2 ^ + 0.124376 * Px2 - 0.00000222979 * Px2 ^ - 0.00000035696 * Px2 * Tx2 ^ 2, 2) qt2 = Round(-839.87 * (1 - 0.071 * cx2) * ((1 - cx2) * cx2 + (1 - cx2 ^ 2) * cx2 ^ 2), 2) Dim ldd2 As Double ldd2 = Round((cx2 * linh32 + (1 - cx2) * liH2O2 + qt2), 3) viH2O2 = Round(1993.19 + 1.88878 * Tx2 - 0.000205512 * Tx2 ^ + 0.000000367295 * Tx2 ^ + 10.6342 * Px2 - 7648.34 * Px2 / Tx2, 2) Dim cv2 As Double cv2 = - (1 - cx2) * Exp((cc2 + dd2) / Tx2) 40 Dim vdd2 As Double vdd2 = Round((cv2 * vinh32 + (1 - cv2) * viH2O2), 3) Text20 = ldd2 Text21 = vdd2 'Tinh entanpi cua dung dich o trang thai Dim Tx3 As Double Tx3 = Text2 + 273.15 Dim Px3 As Double Px3 = Px2 Dim cx3 As Double cx3 = cw Dim cc3 As Double cc3 = -6571.06 + 39.9544 * Tx3 - 0.00243781 * (Tx3 ^ 2) + 0.00000039792 * (Tx3 ^ 3) - 22.7722 * Px3 - 0.00000004979 * (Px3 ^ 2) + 8286.02 * (Px3 / Tx3) + 0.0172363 * Px3 * Tx3 + 0.00000077344 * Px3 * (Tx3 ^ 2) - 3.62962 * Tx3 * (Log(Tx3)) - Tx3 * (Log(Px3)) Dim dd3 As Double dd3 = (-1467 * (cx3 ^ 2) * (1.25 * cx3 ^ - 5.58 * cx3 ^ + 5.96 * cx3 - 0.42)) + ((17.3 * (cx3 ^ 3) - 13.8 * (cx3 ^ 4)) * Tx3) Dim linh33 As Double Dim vinh33 As Double Dim liH2O3 As Double Dim viH2O3 As Double Dim qt3 As Double linh33 = Round(-1058.43 + 3.02148 * Tx3 + 0.00312072 * Tx3 ^ + 0.15536 * Px3 0.0000123079 * Px3 ^ - 0.0000012946 * Px3 * Tx3 ^ 2, 2) vinh33 = Round(770.761 + 1.86947 * Tx3 + 0.0000578293 * Tx3 ^ + 0.000000731509 * Tx3 ^ + 8.98074 * Px3 - 4580.15 * Px3 / Tx3, 2) liH2O3 = Round(-1039.5 + 3.56393 * Tx3 + 0.000919594 * Tx3 ^ + 0.124376 * Px3 - 0.00000222979 * Px3 ^ - 0.00000035696 * Px3 * Tx3 ^ 2, 2) qt3 = Round(-839.87 * (1 - 0.071 * cx3) * ((1 - cx3) * cx3 + (1 - cx3 ^ 2) * cx3 ^ 2), 2) Dim ldd3 As Double ldd3 = Round((cx3 * linh33 + (1 - cx3) * liH2O3 + qt3), 3) viH2O3 = Round(1993.19 + 1.88878 * Tx3 - 0.000205512 * Tx3 ^ + 0.000000367295 * Tx3 ^ + 10.6342 * Px3 - 7648.34 * Px3 / Tx3, 2) Dim cv3 As Double cv3 = - (1 - cx3) * Exp((cc3 + dd3) / Tx3) Dim vdd3 As Double vdd3 = Round((cv3 * vinh33 + (1 - cv3) * viH2O3), 3) Text23 = ldd3 Text24 = vdd3 'Tinh entanpi cua dung dich o trang thai Dim Tx4 As Double Tx4 = 25 + 273.15 Dim aa4 As Double 41 Dim bb4 As Double aa4 = -2103.5 + 4669.96 * cw - 20228.3 * (cw ^ 2) + 56507 * (cw ^ 3) - 80989.9 * (cw ^ 4) + 55286.5 * (cw ^ 5) - 14361.4 * (cw ^ 6) bb4 = 5.65208 - 7.0317 * cw + 37.9018 * (cw ^ 2) - 102.912 * (cw ^ 3) + 135.789 * (cw ^ 4) - 82.7106 * (cw ^ 5) + 18.4113 * (cw ^ 6) Dim Px4 As Double Px4 = Round(10 ^ ((aa4 / Tx4) + bb4 - 0.00847711), 3) Dim cx4 As Double cx4 = cx3 Dim cc4 As Double cc4 = -6571.06 + 39.9544 * Tx4 - 0.00243781 * (Tx4 ^ 2) + 0.00000039792 * (Tx4 ^ 3) - 22.7722 * Px4 - 0.00000004979 * (Px4 ^ 2) + 8286.02 * (Px4 / Tx4) + 0.0172363 * Px4 * Tx4 + 0.00000077344 * Px4 * (Tx4 ^ 2) - 3.62962 * Tx4 * (Log(Tx4)) - Tx4 * (Log(Px4)) Dim dd4 As Double dd4 = (-1467 * (cx4 ^ 2) * (1.25 * cx4 ^ - 5.58 * cx4 ^ + 5.96 * cx4 - 0.42)) + ((17.3 * (cx4 ^ 3) - 13.8 * (cx4 ^ 4)) * Tx4) Dim linh34 As Double Dim vinh34 As Double Dim liH2O4 As Double Dim viH2O4 As Double Dim qt4 As Double linh34 = Round(-1058.43 + 3.02148 * Tx4 + 0.00312072 * Tx4 ^ + 0.15536 * Px4 0.0000123079 * Px4 ^ - 0.0000012946 * Px4 * Tx4 ^ 2, 2) vinh34 = Round(770.761 + 1.86947 * Tx4 + 0.0000578293 * Tx4 ^ + 0.000000731509 * Tx4 ^ + 8.98074 * Px4 - 4580.15 * Px4 / Tx4, 2) liH2O4 = Round(-1039.5 + 3.56393 * Tx4 + 0.000919594 * Tx4 ^ + 0.124376 * Px4 - 0.00000222979 * Px4 ^ - 0.00000035696 * Px4 * Tx4 ^ 2, 2) qt4 = Round(-839.87 * (1 - 0.071 * cx4) * ((1 - cx4) * cx4 + (1 - cx4 ^ 2) * cx4 ^ 2), 2) Dim ldd4 As Double ldd4 = Round((cx4 * linh34 + (1 - cx4) * liH2O4 + qt4), 3) viH2O4 = Round(1993.19 + 1.88878 * Tx4 - 0.000205512 * Tx4 ^ + 0.000000367295 * Tx4 ^ + 10.6342 * Px4 - 7648.34 * Px4 / Tx4, 2) Dim cv4 As Double cv4 = - (1 - cx4) * Exp((cc4 + dd4) / Tx4) Dim vdd4 As Double vdd4 = Round((cv4 * vinh34 + (1 - cv4) * viH2O4), 3) Text25 = ldd4 Text26 = vdd4 'Tinh nhiet luong cap vao qua trinh phat sinh Dim QG As Double Dim QG1 As Double Dim QG2 As Double Dim vdd As Double vdd = 0.00119 42 Dim thetich As Double thetich = Text13 * vdd + Gh1 * thetichrieng Dim Gh3 As Double Gh3 = Text11 Text54 = Text52 Text57 = Text27 QG1 = (Gh2 * vdd2 + (Text13 - Gh2) * ldd2 - (Gh1 * vdd1 + (Text13 - Gh1) * ldd1) thetich * (Pk - Px) * 100) QG2 = Gh3 * vdd3 + (Text13 - Gh3) * ldd3 - Gh2 * vdd2 - (Text13 - Gh2) * ldd2 QG = QG1 + QG2 Text52 = QG Text27 = Round(Text54 - Text52 + Text57, 2) If Text27 < Then Text27 = End If 'Tinh luong tieu hao nhien lieu Dim nhiettri As Double nhiettri = 13500 Dim hieusuat As Double hieusuat = 0.8 Dim nhienlieu As Double nhienlieu = Text27 / (hieusuat * nhiettri) Text46 = Round(nhienlieu, 1) Text38 = Text47 Text39 = Text19 Text40 = Text2 Text41 = Tx4 - 273.15 Text34 = Text49 Text35 = Text22 Text36 = Text15 Text37 = Text15 Text30 = Px Text31 = Pk Text32 = Pk Text33 = Px4 Text48 = mNH3 Text47 = Text47 + Text61 = Text47 - Dim TAM1 As Double TAM1 = Text2 If Text47 > TAM1 Then Command2.Visible = False Text47 = Text47 - 10 Command3.Visible = True End If 'Tính entanpi cua hoi tinh khiet NH3 bay vao binh ngung 43 Text51 = Text2 Text53 = Text8 Dim Tx5 As Double Tx5 = Text2 + 273.15 Dim Px5 As Double Px5 = Px2 Dim cx5 As Double cx3 = 100 Dim vinh35 As Double vinh35 = Round(770.761 + 1.86947 * Tx3 + 0.0000578293 * Tx3 ^ + 0.000000731509 * Tx3 ^ + 8.98074 * Px3 - 4580.15 * Px3 / Tx3, 2) Text55 = vinh35 ' Tinh entanpi cua long NH3 sau ngung tu Dim Tx6 As Double Tx6 = Text10 + 273.15 Dim Px6 As Double Px6 = Px2 Dim cx6 As Double cx6 = 100 Dim linh36 As Double linh36 = Round(-1058.43 + 3.02148 * Tx6 + 0.00312072 * Tx6 ^ + 0.15536 * Px6 0.0000123079 * Px6 ^ - 0.0000012946 * Px6 * Tx6 ^ 2, 2) Text56 = Text10 Text58 = Text8 Text60 = linh36 'Tinh entanpi NH3 sau tiet luu Text66 = Text6 Text68 = Text9 Text70 = Text60 'Tinh entanpi NH sau khoi dan bay hoi Dim Tx8 As Double Tx8 = Text6 + 273.15 Dim Px8 As Double Px8 = Text9 Dim cx8 As Double cx8 = 100 Dim vinh38 As Double vinh38 = Round(770.761 + 1.86947 * Tx8 + 0.0000578293 * Tx8 ^ + 0.000000731509 * Tx8 ^ + 8.98074 * Px8 - 4580.15 * Px8 / Tx8, 2) Text71 = Text66 Text73 = Text68 Text75 = vinh38 End Sub Private Sub Command3_Click() Dim Pk As Double Pk = Text8 44 'Tinh khoi luong dung dich Dim cs As Double Dim cw As Double Dim mNH3 As Double Dim mH2O As Double Dim m As Double m = Text11 Dim cs100 As Double Dim cw100 As Double cs100 = Text14 cw100 = Text15 cs = cw100 / 100 cw = cw100 / 100 mNH3 = Text11 'Tinh ap suat dung dich binh phat sinh Dim Tx As Double Dim Px As Double Dim muy As Double Dim roNH3 As Double Dim P1 As Double muy = 17.03 Tx = Text47 + 273.15 Dim aaa As Double Dim bbb As Double aaa = -2103.5 + 4669.96 * cs - 20228.3 * (cs ^ 2) + 56507 * (cs ^ 3) - 80989.9 * (cs ^ 4) + 55286.5 * (cs ^ 5) - 14361.4 * (cs ^ 6) bbb = 5.65208 - 7.0317 * cs + 37.9018 * (cs ^ 2) - 102.912 * (cs ^ 3) + 135.789 * (cs ^ 4) - 82.7106 * (cs ^ 5) + 18.4113 * (cs ^ 6) P1 = Round(10 ^ ((aaa / Tx) + bbb - 0.00847711), 3) ' Tinh nong dung dich binh phat sinh Dim cx As Double Dim Ax As Double Dim Bx As Double Dim cc As Double If P1 < Pk Then Px = P1 roNH3 = muy * Px / (0.08314 * (Tx + 273)) Dim thetichbayhoi As Double thetichbayhoi = 0.05 Dim thetichrieng As Double thetichrieng = / roNH3 Dim Gh1 As Double Gh1 = thetichbayhoi / thetichrieng End If If P1 > Pk Then Px = Pk 45 End If Text16 = Px 'Tinh entanpi cua dung dich o trang thai Dim Tx3 As Double Tx3 = Text2 + 273.15 Dim Px3 As Double Px3 = Text8 Dim cx3 As Double cx3 = cw Dim cc3 As Double cc3 = -6571.06 + 39.9544 * Tx3 - 0.00243781 * (Tx3 ^ 2) + 0.00000039792 * (Tx3 ^ 3) - 22.7722 * Px3 - 0.00000004979 * (Px3 ^ 2) + 8286.02 * (Px3 / Tx3) + 0.0172363 * Px3 * Tx3 + 0.00000077344 * Px3 * (Tx3 ^ 2) - 3.62962 * Tx3 * (Log(Tx3)) - Tx3 * (Log(Px3)) Dim dd3 As Double dd3 = (-1467 * (cx3 ^ 2) * (1.25 * cx3 ^ - 5.58 * cx3 ^ + 5.96 * cx3 - 0.42)) + ((17.3 * (cx3 ^ 3) - 13.8 * (cx3 ^ 4)) * Tx3) Dim linh33 As Double Dim vinh33 As Double Dim liH2O3 As Double Dim viH2O3 As Double Dim qt3 As Double linh33 = Round(-1058.43 + 3.02148 * Tx3 + 0.00312072 * Tx3 ^ + 0.15536 * Px3 0.0000123079 * Px3 ^ - 0.0000012946 * Px3 * Tx3 ^ 2, 2) vinh33 = Round(770.761 + 1.86947 * Tx3 + 0.0000578293 * Tx3 ^ + 0.000000731509 * Tx3 ^ + 8.98074 * Px3 - 4580.15 * Px3 / Tx3, 2) liH2O3 = Round(-1039.5 + 3.56393 * Tx3 + 0.000919594 * Tx3 ^ + 0.124376 * Px3 - 0.00000222979 * Px3 ^ - 0.00000035696 * Px3 * Tx3 ^ 2, 2) qt3 = Round(-839.87 * (1 - 0.071 * cx3) * ((1 - cx3) * cx3 + (1 - cx3 ^ 2) * cx3 ^ 2), 2) 'Dim ldd3 As Double ldd3 = Round((cx3 * linh33 + (1 - cx3) * liH2O3 + qt3), 3) viH2O3 = Round(1993.19 + 1.88878 * Tx3 - 0.000205512 * Tx3 ^ + 0.000000367295 * Tx3 ^ + 10.6342 * Px3 - 7648.34 * Px3 / Tx3, 2) Dim cv3 As Double cv3 = - (1 - cx3) * Exp((cc3 + dd3) / Tx3) Dim vdd3 As Double vdd3 = Round((cv3 * vinh33 + (1 - cv3) * viH2O3), 3) Text23 = ldd3 Text24 = vdd3 'Tinh entanpi cua dung dich o trang thai Dim Tx4 As Double Tx4 = Text47 + 273.15 Dim aa4 As Double Dim bb4 As Double 46 aa4 = -2103.5 + 4669.96 * cw - 20228.3 * (cw ^ 2) + 56507 * (cw ^ 3) - 80989.9 * (cw ^ 4) + 55286.5 * (cw ^ 5) - 14361.4 * (cw ^ 6) bb4 = 5.65208 - 7.0317 * cw + 37.9018 * (cw ^ 2) - 102.912 * (cw ^ 3) + 135.789 * (cw ^ 4) - 82.7106 * (cw ^ 5) + 18.4113 * (cw ^ 6) Dim Px4 As Double Px4 = Px Dim cx4 As Double cx4 = cx3 Dim cc4 As Double cc4 = -6571.06 + 39.9544 * Tx4 - 0.00243781 * (Tx4 ^ 2) + 0.00000039792 * (Tx4 ^ 3) - 22.7722 * Px4 - 0.00000004979 * (Px4 ^ 2) + 8286.02 * (Px4 / Tx4) + 0.0172363 * Px4 * Tx4 + 0.00000077344 * Px4 * (Tx4 ^ 2) - 3.62962 * Tx4 * (Log(Tx4)) - Tx4 * (Log(Px4)) Dim dd4 As Double dd4 = (-1467 * (cx4 ^ 2) * (1.25 * cx4 ^ - 5.58 * cx4 ^ + 5.96 * cx4 - 0.42)) + ((17.3 * (cx4 ^ 3) - 13.8 * (cx4 ^ 4)) * Tx4) Dim linh34 As Double Dim vinh34 As Double Dim liH2O4 As Double Dim viH2O4 As Double Dim qt4 As Double linh34 = Round(-1058.43 + 3.02148 * Tx4 + 0.00312072 * Tx4 ^ + 0.15536 * Px4 0.0000123079 * Px4 ^ - 0.0000012946 * Px4 * Tx4 ^ 2, 2) vinh34 = Round(770.761 + 1.86947 * Tx4 + 0.0000578293 * Tx4 ^ + 0.000000731509 * Tx4 ^ + 8.98074 * Px4 - 4580.15 * Px4 / Tx4, 2) liH2O4 = Round(-1039.5 + 3.56393 * Tx4 + 0.000919594 * Tx4 ^ + 0.124376 * Px4 - 0.00000222979 * Px4 ^ - 0.00000035696 * Px4 * Tx4 ^ 2, 2) qt4 = Round(-839.87 * (1 - 0.071 * cx4) * ((1 - cx4) * cx4 + (1 - cx4 ^ 2) * cx4 ^ 2), 2) Dim ldd4 As Double ldd4 = Round((cx4 * linh34 + (1 - cx4) * liH2O4 + qt4), 3) If Tx4 > Tx3 Then ldd4 = ldd3 End If viH2O4 = Round(1993.19 + 1.88878 * Tx4 - 0.000205512 * Tx4 ^ + 0.000000367295 * Tx4 ^ + 10.6342 * Px4 - 7648.34 * Px4 / Tx4, 2) Dim cv4 As Double cv4 = - (1 - cx4) * Exp((cc4 + dd4) / Tx4) Dim vdd4 As Double vdd4 = Round((cv4 * vinh34 + (1 - cv4) * viH2O4), 3) 'Tinh nhiet luong dung dich nhan duoc qua trinh phat sinh Dim QGL As Double QGL = (Text13 - Gh3) * ldd3 - (Text13 - Gh1) * ldd1 'Text29 = QGL 'Tinh nhiet luong nha sau qua trinh phat sinh Dim vdd As Double 47 vdd = 0.00119 Dim thetich As Double thetich = Text13 * vdd + Gh1 * thetichrieng 'Dim Gh3 As Double Gh3 = Text11 'Text59 = Text28 Dim QR As Double QR = (Text13 - Gh3) * ldd4 - (Text13 - Gh3) * ldd3 - thetich * (Px - Pk) * 100 Text28 = Round(QR, 0) Text47 = Text47 - Text61 = Text47 + Dim TAM2 As Double TAM2 = Text5 If Text61