Nghiên cứu xác định một số thông số đặc tính của quá trình khử muối nước biển bằng phương pháp phun tách ẩm

ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA VÕ KIẾN QUỐC NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THƠNG SỐ ĐẶC TÍNH CỦA Q TRÌNH KHỬ MUỐI NƯỚC BIỂN BẰNG PHƯƠNG PHÁP PHUN – TÁCH ẨM 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 VÕ KIẾN QUỐC NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THƠNG SỐ ĐẶC TÍNH CỦA Q TRÌNH KHỬ MUỐI NƯỚC BIỂN BẰNG PHƯƠNG PHÁP PHUN – TÁCH ẨM Chuyên ngành: Công nghệ thiết bị nhiệt Mã số chuyên ngành: 62 52 80 01 Phản biện độc lập 1: PGS.TS NGUYỄN NGUYÊN AN Phản biện độc lập 2: PGS.TS TRẦN VĂN VANG Phản biện 1: PGS.TS HOÀNG AN QUỐC Phản biện 2: PGS.TS LÊ ANH ĐỨC Phản biện 3: PGS.TS HỒNG ĐÌNH TÍN NGƯỜI HƯỚNG DẪN KHOA HỌC GS.TS LÊ CHÍ HIỆP TS NGUYỄN VĂN TUYÊN 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 Võ Kiến Quốc i TÓM TẮT LUẬN ÁN Luận án trình bày kết nghiên cứu lý thuyết thực nghiệm hệ thống khử muối phương pháp phun – tách ẩm Dựa phân tích lý thuyết truyền nhiệt, truyền chất hệ thống kết hợp với ứng dụng công nghệ Pinch thu hồi nhiệt, chương trình mơ xây dựng để xác định tỷ số lưu lượng khối lượng tối ưu nước phun khơng khí theo nhiệt độ nước phun, nhiệt độ nước cấp chênh lệch nhiệt độ tối thiểu điểm Pinch (Tmin) Kết cho thấy nhiệt độ nước phun cao tỷ lệ lưu lượng khối lượng tối ưu nước phun khơng khí cao, nhiệt độ nước cấp thấp tỷ lệ lưu lượng khối lượng tối ưu nước phun khơng khí thấp Nếu chênh lệch nhiệt độ tối thiểu điểm Pinch thiết bị phun ẩm thiết bị tách ẩm tỷ tỷ lệ lưu lượng khối lượng tối ưu nước phun khơng khí phụ thuộc vào nhiệt độ nước phun nhiệt độ nước cấp Từ kết mô phương trình hồi quy xây dựng để xác định tỷ lệ lưu lượng tối ưu nước phun khơng khí theo nhiệt độ nước phun nhiệt độ nước cấp Từ kết mô phỏng, tỷ lệ thu hồi nhiệt hệ số suất GOR tối ưu mặt kỹ thuật hệ thống xác định theo chênh lệch nhiệt độ tối thiểu Tmin Từ xác định nhiệt độ nước phun hợp lý cho hệ thống Nếu Tmin cao nhiệt độ nước phun hợp lý cao ngược lại Với chênh lệch nhiệt độ tối thiểu Tmin=5oC nhiệt độ nước phun hợp lý 70-75oC Dựa việc phân tích q trình trao đổi nhiệt chất thiết bị phun ẩm, luận án trình bày phương pháp phân tích q trình trao đổi nhiệt chất thiết bị phun ẩm theo phương pháp Merkel phương pháp Poppe Từ đó, chương trình mơ trình truyền nhiệt truyền chất thiết bị phun ẩm xây dựng để đánh giá ảnh hưởng thông số đầu vào (nhiệt độ nước phun, nhiệt độ khơng khí, số Me tỷ lệ lưu lượng khối lượng nước phun không khí) đến thơng số đầu thiết bị phun ẩm (nhiệt độ khơng khí, nhiệt độ nước) Từ kết mô phỏng, giá trị hợp lý số Me thiết kế thiết bị phun ẩm xác định Me1,2 ii Sau phân tích lý thuyết mơ hệ thống, luận án trình bày kết nghiên cứu thực nghiệm trình trao đổi nhiệt trao đổi chất thiết bị phun ẩm Bằng cách thay đổi thông số đầu vào bao gồm nhiệt độ nước phun, nhiệt độ khơng khí vào, chiều cao lớp đệm tỷ lệ lưu lượng khối lượng nước phun khơng khí nhiệt độ khơng khí nước khỏi thiết bị phun ẩm xác định Kết hợp với mô lý thuyết thực trước đó, phương trình hồi quy thể mối quan hệ số Me với chiều cao lớp đệm, nhiệt độ nước phun tỷ lệ lưu lượng khối lượng nước phun khơng khí xác định Kết thực nghiệm cho thấy nhiệt độ khơng khí vào thiết bị phun ẩm khơng ảnh hưởng đến số Me khoảng nhiệt độ 34-38oC Ngoài ra, phương pháp thực nghiệm mật độ xối tưới thiết bị phun ẩm xác định 0,3kg/(m s) ≤ M ≤ 0,75kg/(m s) So sánh với kết mơ lý thuyết kết thực nghiệm có sai số khoảng 4% Điều cho thấy phù hợp mơ hình thực nghiệm mơ hình tốn học Để đánh giá suất hệ thống khử muối phun – tách ẩm cấp nhiệt lượng mặt trời, luận án trình bày kết thực nghiệm loại thu điều kiện thời tiết Tp.HCM Kết thực nghiệm cho thấy thu tập trung dạng máng trụ không phù hợp để cấp nhiệt cho hệ thống khử muối phun – tách ẩm hiệu suất loại thu thấp Đối với thu ống dầu hiệu suất tốt loại thu phù hợp để cấp nhiệt cho hệ thống khử muối phun – tách ẩm Ngoài ra, với mục đích thiết kế loại thu phù hợp với điều kiện kinh tế - kỹ thuật điều kiện Việt Nam, thu ống nhiệt trọng trường tách dòng độc lập chế tạo thực nghiệm Kết thực nghiệm cho thấy môi chất nạp phù hợp ống nhiệt nước tỷ lệ nạp hợp lý 42%-62% hiệu suất nhiệt ống nhiệt trọng trường tách dòng độc lập 45%-60% Dựa liệu cường độ xạ mặt trời Tp.HCM, suất hệ thống khử muối phun – tách ẩm cấp nhiệt lượng mặt trời đánh giá Bằng phương pháp ước lượng, suất hệ thống khử muối phun – tách ẩm sử dụng lượng mặt trời có suất trung bình từ 4-5,4kg/(m2ngày) thu ống nhiệt trọng trường tách dòng độc lập từ 4,3-5,9kg/(m2ngày) thu ống dầu iii 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 Tp.HCM Nhà B5, 268 Lý Thường Kiệt, Q.10, TP.HCM ĐT: (08) 3864-7256 Ext 5897 iv ABSTRACT The study results on theoretical and experimental problems of humidification – dehumidification desalination systems are presented in the thesis Based on the heat and mass transfer analysis and the Pinch technology, the modeling program has been done in order to determine the optimum mass flow rate between sprayed water and air in relationship with temperature of sprayed water, temperature of cooling water and minimum temperature difference at the Pinch point (Tmin) The done results show that the higher sprayed water temperature, the higher optimum mass flow rate between sprayed water and air, the lower cooling water temperature, the lower optimum mass flow rare between sprayed water and air If Tmin of the humidifier equals with Tmin of the dehumidifier, the optimum mass flow rate between sprayed water and air depends only on temperature of sprayed water and temperature of cooling water Based on the simulation results, it is concluded that the heat recovery rate and the technical optimum gain output ratio of the system, and then the suitable water temperature for the system, can be determined by the minimum temperature difference Tmin It is clear that, the higher Tmin, the higher suitable water temperature With minimum temperature difference Tmin = 5oC, the suitable water temperature is 7075oC Based on the analysis of the heat and mass transfers in the humidifier, two methods of analyzing the heat and mass transfers in the humidifier done by using Merkel method and Poppe method have been presented Hence, the simulation program for heat and mass transfers in the humidifier has been developed to evaluate the effect of input parameters such as water temperature, air temperature, Me number and mass flow rate between sprayed water and air to output parameters of the humidifier such as air temperature and water temperature The simulation results show that the suitable value of Me number used to design the humidifier should be higher than 1.2 Besides theoretical study, experimental study on heat and mass transfers in the humidifier has been also carried out by varying the input parameters such as water temperature, air temperature, packing bed height and mass flow rate between sprayed v water and air in order to determine the most appropriate values of air temperature and water temperature at the output of the humidifier Combined with the previous theoretical simulations, the regression equations showing the relationships of the Me number, the packing bed height, the sprayed water temperature and the mass flow rate between sprayed water and air have been determined Experimental results show that the air temperature at the input of the humidifier does not affect the Me number in the range of 34-38oC In addition, the experimental results show that the density of irrigation M in the humidifier should be from 0.3kg/(m2s) to 0.75kg/(m2s) which are nearly 4%-error compared with theoretical simulation results To evaluate the production of the desalination systems heated by solar energy, some solar collectors have been added and tested in weather conditions of HCMC Experimental results show that concentrating solar collector is not suitable for heating the desalination system because its total efficiency is very low For oil tube solar collector, its performance is very good so this collector is very suitable for heating the desalination system In addition, for the purpose of designing a new type of solar collector that is suitable with economic and technical conditions in Vietnam, a new model working based on heat pipe principle has been developed and tested In this prototype, experimental results show that the most suitable working fluid should be pure water and the suitable filling ratio is from 42% to 62%, the thermal efficiency of the prototype is 45%-60% Based on the measured solar radiation of HCMC, the capacity of the solar thermal humidifier-dehumidifier desalination systems has been evaluated By the estimation method, the capacities of the solar thermal humidifierdehumidifier desalination systems connected with the new model and oil pipe collector are from to 5.4 kg/(m2day) and from 4.3 to 5.9 kg/(m2day), respectively vi LỜI CÁM ƠN Trước tiên, xin trân trọng cám ơn GS.TS Lê Chí Hiệp TS Nguyễn Văn Tuyên tận tình hướng dẫn, giúp đỡ tạo điều kiện thuận lợi cho suốt thời gian thực luận án Tiếp theo, xin chân thành cảm ơn đến Thầy Nguyễn Nhựt, Thầy Lê Chí Hiệp Công ty Châu Hiệp Phát tài trợ phần kinh phí thiết bị cho mơ hình thực nghiệm dụng cụ đo Sau cùng, xin cám ơn đồng nghiệp Bộ môn Công nghệ nhiệt lạnh đóng góp ý kiến hỗ trợ việc tìm kiếm báo khoa học có liên quan đến đề tài Luận án Tiến sĩ vii MỤC LỤC DANH MỤC CÁC HÌNH ẢNH xii DANH MỤC BẢNG BIỂU xv DANH MỤC KÝ HIỆU VÀ CÁC TỪ VIẾT TẮT xvi MỞ ĐẦU xviii CHƯƠNG 1.1 TỔNG QUAN Kỹ thuật khử muối truyền thống 1.1.1 Phương pháp khử muối bay nhiều tầng 1.1.2 Phương pháp khử muối đa hiệu ứng 1.1.3 Phương pháp khử muối sử dụng máy nén 1.1.4 Phương pháp thẩm thấu ngược 1.1.5 Phương pháp điện phân 1.1.6 Nhược điểm phương pháp khử muối truyền thống 1.1.7 Phương pháp chưng cất lượng mặt trời 1.2 Khử muối phương pháp phun – tách ẩm 1.2.1 Giới thiệu phương pháp phun – tách ẩm 1.2.2 Phạm vi ứng dụng phương pháp phun – tách ẩm 1.2.3 Phân loại hệ thống khử muối phun – tách ẩm 1.3 Tình hình nghiên cứu ứng dụng khử muối nước biển Việt Nam giới 1.3.1 Tình hình Việt Nam 1.3.2 Tình hình giới 11 1.4 Tổng quan phương pháp phun – tách ẩm 12 1.4.1 Các thông số dùng để đánh giá hệ thống khử muối phun – tách ẩm 12 1.4.2 Chu trình khơng khí tuần hồn kín - nước hở, gia nhiệt nước (Close air, open water, water heated CAOW-WH) 14 1.4.3 Chu trình nước tuần hồn kín, khơng khí hở, gia nhiệt nước (Close water, open air, water heated CWOA-WH) 17 1.4.4 Chu trình khơng khí tuần kín, nước hở, gia nhiệt khơng khí (Close air, open water, air heated CAOW-AH) 20 viii = (15,875 − 6,86) 10 + 0,805 lg 1,28 15,875 15,8766 − 0,3 6.86 15,875 l=11,876 [mm] = Cánh nhơm có hệ số dẫn nhiệt = 203,5 [W/m K] = 26,82 → = 96,74% ∗ = = + Diện tích ống phần khơng gắn cánh ứng với 1m ống = 1− = 0,040557 [m /m] Diện tích cánh m ống: = − 0,785 = 0,314462 [m /m] Tổng diện tích bề mặt ứng với 1m ống: ∗ = + = 0,35522 [m /m] = 8,9 Chọn nhiệt trở bụi bẩn bám phía cánh Rc=0,0003 [m2.K/W] nhiệt trở cáu bẩn bám phía ống Rt=0,0003 [m2.K/W] Hệ số truyền nhiệt thiết bị ngưng tụ = 97,76 [W/m2.K] 201 Diện tích thiết bị truyền nhiệt = ∗ Δ = 42,5 [m] Chọn L=45 [m] Lưu lượng thể tích khơng khí ẩm qua dàn ngưng tụ: = = 85,45 [m /h] Diện tích tiết diện ngang chỗ hẹp nhất: = = 0,0505 [m ] Chọn kết cấu dàn ngưng sau: Số dãy ống theo chiều chuyển động khơng khí n1=15 Số ống dãy n2=5 Chiều dài đoạn ống: = 45/(5.15) = 0,6 Kiểm tra diện tích tiết diện ngang chỗ hẹp nhất: = ( − ) − = 0,0508 [m ] 202 PHỤ LỤC 7: CODE CHƯƠNG TRÌNH MƠ PHỎNG Private Sub Command1_Click() Text1 = 30 Text2 = 70 Text3 = Text4 = Text5 = End Sub Private Sub Command2_Click() 'Nhiet nuoc cap vao thiet bi tach am Dim t1 As Double t1 = Text1 'Nhiet nuoc phun Dim t3 As Double t3 = Text2 ' Ty le luu luong khoi luong Dim m As Double m = Text3 ' Nhiet diem Pinch Dim DT1 As Double DT1 = Text4 Dim DT2 As Double DT2 = Text5 'Nhiet nuoc cap khoi thiet bi tach am Dim t2 As Double 'Nhiet nuoc khoi thiet bi phun am Dim t4 As Double 'Nhiet khong vao thiet bi phun am Dim t5 As Double 'Nhiet khong khoi thiet bi phun am Dim t6 As Double ' Nhiet khong Dim ta As Double ' Nhiet nuoc Dim tw As Double ta = t1 tw = t1 ' Entanpi cua nuoc Dim Iw As Double ' Nhiet dung rieng Dim Cp As Double Cp = 4.186 Iw = m * Cp * tw 203 ' Entanpi cua khong Dim Ia As Double Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) ' Chenh lech entanpi cua nuoc va khong Dim DeltaI1 As Double Dim Delta As Double If Iw > Ia Then Do Delta = Iw - Ia tw = tw + 0.001 ta = tw Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) DeltaI1 = Iw - Ia Loop Until DeltaI1 - Delta < Text6 = Round(tw, 2) If t3 > tw Then 'Entanpi nuoc phun Dim I3 As Double I3 = m * Cp * t3 'Entanpi khong khoi thiet bi phun am Dim I6 As Double I6 = I3 - DeltaI1 - m * Cp * DT1 For t6 = t1 To t3 Step 0.001 I6 = Exp(2.39329 + 0.10648 * t6 - 0.00135 * t6 ^ + 0.000010058 * t6 ^ 3) If (I6 - I3 + DeltaI1 + m * Cp * DT1) > And (I6 - I3 + DeltaI1 + m * Cp * DT1) < 0.1 Then Exit For End If Next t6 'Do chua hoi khong khoi thiet bi phun am Dim w6 As Double w6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 End If If t3 DT2 Then t4 = t3 - Qreco / (m * Cp) End If If t5 - t1 < DT2 Then t5 = t1 + DT2 I5 = Exp(2.39329 + 0.10648 * t5 - 0.00135 * t5 ^ + 0.000010058 * t5 ^ 3) Qreco = I6 - I5 t2 = t1 + Qreco / (m * Cp) t4 = t3 - Qreco / (m * Cp) End If Text8 = Round(t2, 2) Text9 = Round(t5, 2) Text10 = Round(t4, 2) Text11 = Round(Qreco, 2) Dim Qadd As Double Qadd = m * Cp * (t3 - t2) Text12 = Round(Qadd, 2) Dim HRR As Double HRR = Qreco / (Qadd + Qreco) Text13 = Round(HRR * 100, 2) Dim d5 As Double d5 = (2.19 * (t5 ^ 3) / 10 ^ 6) - (1.85 * (t5 ^ 2) / 10 ^ 4) + (7.06 * t5 / 10 ^ 3) - 0.077 Dim d6 As Double d6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 Dim mw As Double mw = d6 - d5 Dim ttb As Double ttb = (t6 + t5) / Dim hfg As Double hfg = 2501.8971440379 - 2.407064037 * ttb + 1.192217 * 10 ^ (-3) * ttb ^ - 1.5863 * 10 ^ (-5) * ttb ^ Dim GOR As Double GOR = mw * hfg / Qadd 205 Text14 = Round(mw, 2) Text15 = Round(GOR, 2) Text19 = t1 Text20 = Round(t2, 2) Text21 = Round(t3, 2) Text22 = Round(t5, 2) Text23 = Round(t4, 2) Text24 = Round(t2, 2) End Sub Private Sub Command3_Click() Dim t1 As Double t1 = Text1 'Nhiet nuoc phun Dim t3 As Double t3 = Text2 ' Ty le luu luong khoi luong Dim m As Double m = Text3 - 0.1 ' Nhiet diem Pinch Dim DT1 As Double DT1 = Text4 Dim DT2 As Double DT2 = Text5 'Nhiet nuoc cap khoi thiet bi tach am Dim t2 As Double 'Nhiet nuoc khoi thiet bi phun am Dim t4 As Double 'Nhiet khong vao thiet bi phun am Dim t5 As Double 'Nhiet khong khoi thiet bi phun am Dim t6 As Double ' Nhiet khong Dim ta As Double ' Nhiet nuoc Dim tw As Double ta = t1 tw = t1 ' Entanpi cua nuoc Dim Iw As Double ' Nhiet dung rieng Dim Cp As Double Cp = 4.186 Iw = m * Cp * tw ' Entanpi cua khong Dim Ia As Double 206 Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) ' Chenh lech entanpi cua nuoc va khong Dim DeltaI1 As Double Dim Delta As Double If Iw > Ia Then Do Delta = Iw - Ia tw = tw + 0.001 ta = tw Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) DeltaI1 = Iw - Ia Loop Until DeltaI1 - Delta < If t3 > tw Then 'Entanpi nuoc phun Dim I3 As Double I3 = m * Cp * t3 'Entanpi khong khoi thiet bi phun am Dim I6 As Double I6 = I3 - DeltaI1 - m * Cp * DT1 For t6 = t1 To t3 Step 0.001 I6 = Exp(2.39329 + 0.10648 * t6 - 0.00135 * t6 ^ + 0.000010058 * t6 ^ 3) If (I6 - I3 + DeltaI1 + m * Cp * DT1) > And (I6 - I3 + DeltaI1 + m * Cp * DT1) < 0.1 Then Exit For End If Next t6 'Do chua hoi khong khoi thiet bi phun am Dim w6 As Double w6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 End If If t3 DT2 Then t4 = t3 - Qreco / (m * Cp) End If If t5 - t1 < DT2 Then t5 = t1 + DT2 I5 = Exp(2.39329 + 0.10648 * t5 - 0.00135 * t5 ^ + 0.000010058 * t5 ^ 3) Qreco = I6 - I5 t2 = t1 + Qreco / (m * Cp) t4 = t3 - Qreco / (m * Cp) End If Dim Qadd As Double Qadd = m * Cp * (t3 - t2) Text12 = Round(Qadd, 2) Dim HRR As Double HRR = Qreco * 100 / (Qadd + Qreco) Dim d5 As Double d5 = (2.19 * (t5 ^ 3) / 10 ^ 6) - (1.85 * (t5 ^ 2) / 10 ^ 4) + (7.06 * t5 / 10 ^ 3) - 0.077 Dim d6 As Double d6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 Dim mw As Double mw = d6 - d5 Dim ttb As Double ttb = (t6 + t5) / Dim hfg As Double hfg = 2501.8971440379 - 2.407064037 * ttb + 1.192217 * 10 ^ (-3) * ttb ^ - 1.5863 * 10 ^ (-5) * ttb ^ Dim GOR As Double GOR = mw * hfg / Qadd ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' If HRR > Text13 Then Dim HRRop As Double Do HRRop = HRR m = m - 0.01 DT1 = Text4 208 DT2 = Text5 ta = t1 tw = t1 Cp = 4.186 Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) If Iw > Ia Then Do Delta = Iw - Ia tw = tw + 0.001 ta = tw Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) DeltaI1 = Iw - Ia Loop Until DeltaI1 - Delta < If t3 > tw Then I3 = m * Cp * t3 I6 = I3 - DeltaI1 - m * Cp * DT1 For t6 = t1 To t3 Step 0.001 I6 = Exp(2.39329 + 0.10648 * t6 - 0.00135 * t6 ^ + 0.000010058 * t6 ^ 3) If (I6 - I3 + DeltaI1 + m * Cp * DT1) > And (I6 - I3 + DeltaI1 + m * Cp * DT1) < 0.1 Then Exit For End If Next t6 w6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 End If If t3 DT2 Then 209 t4 = t3 - Qreco / (m * Cp) End If If t5 - t1 < DT2 Then t5 = t1 + DT2 I5 = Exp(2.39329 + 0.10648 * t5 - 0.00135 * t5 ^ + 0.000010058 * t5 ^ 3) Qreco = I6 - I5 t2 = t1 + Qreco / (m * Cp) t4 = t3 - Qreco / (m * Cp) End If Qadd = m * Cp * (t3 - t2) HRR = Qreco * 100 / (Qadd + Qreco) d5 = (2.19 * (t5 ^ 3) / 10 ^ 6) - (1.85 * (t5 ^ 2) / 10 ^ 4) + (7.06 * t5 / 10 ^ 3) - 0.077 d6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 mw = d6 - d5 ttb = (t6 + t5) / hfg = 2501.8971440379 - 2.407064037 * ttb + 1.192217 * 10 ^ (-3) * ttb ^ - 1.5863 * 10 ^ (-5) * ttb ^ GOR = mw * hfg / Qadd Loop Until HRR < HRRop m = m + 0.01 End If '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' If HRR < Text13 Then Do HRRop = HRR m = m + 0.01 DT1 = Text4 DT2 = Text5 ta = t1 tw = t1 Cp = 4.186 Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) If Iw > Ia Then Do Delta = Iw - Ia tw = tw + 0.001 ta = tw Iw = m * Cp * tw Ia = Exp(2.39329 + 0.10648 * ta - 0.00135 * ta ^ + 0.000010058 * ta ^ 3) DeltaI1 = Iw - Ia Loop Until DeltaI1 - Delta < If t3 > tw Then I3 = m * Cp * t3 I6 = I3 - DeltaI1 - m * Cp * DT1 210 For t6 = t1 To t3 Step 0.001 I6 = Exp(2.39329 + 0.10648 * t6 - 0.00135 * t6 ^ + 0.000010058 * t6 ^ 3) If (I6 - I3 + DeltaI1 + m * Cp * DT1) > And (I6 - I3 + DeltaI1 + m * Cp * DT1) < 0.1 Then Exit For End If Next t6 w6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 End If If t3 DT2 Then t4 = t3 - Qreco / (m * Cp) End If If t5 - t1 < DT2 Then t5 = t1 + DT2 I5 = Exp(2.39329 + 0.10648 * t5 - 0.00135 * t5 ^ + 0.000010058 * t5 ^ 3) Qreco = I6 - I5 t2 = t1 + Qreco / (m * Cp) t4 = t3 - Qreco / (m * Cp) End If Qadd = m * Cp * (t3 - t2) HRR = Qreco * 100 / (Qadd + Qreco) d5 = (2.19 * (t5 ^ 3) / 10 ^ 6) - (1.85 * (t5 ^ 2) / 10 ^ 4) + (7.06 * t5 / 10 ^ 3) - 0.077 d6 = (2.19 * (t6 ^ 3) / 10 ^ 6) - (1.85 * (t6 ^ 2) / 10 ^ 4) + (7.06 * t6 / 10 ^ 3) - 0.077 mw = d6 - d5 ttb = (t6 + t5) / hfg = 2501.8971440379 - 2.407064037 * ttb + 1.192217 * 10 ^ (-3) * ttb ^ - 1.5863 * 10 ^ (-5) * ttb ^ GOR = mw * hfg / Qadd 211 Loop Until HRR < HRRop m = m - 0.01 End If Text16 = Round(m, 2) Text17 = Round(HRR, 2) Text18 = Round(GOR, 2) Text25 = Round(t2, 2) Text26 = Round(t4, 2) End Sub 212 PHỤ LỤC 8: HÌNH ẢNH TỔ CHỨC THỰC NGHIỆM Cách bố trí vật liệu đệm thiết bị phun ẩm Cách bố trí vịi phun thiết bị phun ẩm 213 Mơ hình thực nghiệm 214 Đo lưu lượng khơng khí 215 ... nghiên cứu trước hệ thống khử muối phun tách ẩm chưa có nghiên cứu thơng số đặc tính hệ thống Vì tác giả chọn đề tài ? ?nghiên cứu xx xác định số thơng số đặc tính q trình khử muối nước biển phương pháp. .. BÁCH KHOA VÕ KIẾN QUỐC NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THƠNG SỐ ĐẶC TÍNH CỦA Q TRÌNH KHỬ MUỐI NƯỚC BIỂN BẰNG PHƯƠNG PHÁP PHUN – TÁCH ẨM Chuyên ngành: Công nghệ thiết bị nhiệt Mã số chuyên ngành: 62 52... pháp phun – tách ẩm? ?? đề tài nghiên cứu luận án Mục tiêu nghiên cứu Mục tiêu nghiên cứu luận án xác định ảnh hưởng thông số đầu vào đến thông số đầu hệ thống khử muối phun tách ẩm Từ xác định