Nghiên cứu mô phỏng trong chế độ không ổn định hệ thống cung cấp nước nóng sử dụng bơm nhiệt kết hợp với bộ thu năng lượng mặt trời áp dụng cho các hộ gia đình tại việt nam

LỜI CẢM ƠN Trước hết, xin bày tỏ lòng biết ơn lời cảm ơn sâu sắc tới thầy hướng dẫn, TS Nguyễn Nguyên An, người thầy tận tình hướng dẫn thực hoàn thành luận văn tốt nghiệp Xin gửi lời cảm ơn chân thành tới thầy cô Viện Khoa học Công nghệ Nhiệt – Lạnh, đặc biệt thầy cô môn Kỹ thuật Lạnh Điều hoà Không khí, môn công tác tạo điều kiện giúp đỡ trình thực luận văn Cuối cùng, xin cảm ơn bạn bè gia đình bên tôi, cổ vũ động viên hoàn thành tốt luận văn tốt nghiệp Tôi xin chân thành cảm ơn! Hà Nội, ngày tháng năm 2014 Học viên Trịnh Viết Thiệu Học viên: Trịnh Viết Thiệu i LỜI CAM ĐOAN Bản luận văn nghiên cứu thực hướng dẫn thầy giáo TS Nguyễn Nguyên An Để hoàn thành luận văn này, tài liệu tham khảo liệt kê, cam đoan không chép công trình nghiên cứu, sử dụng tài liệu khác mà chưa liệt kê Nếu sai phạm, xin chịu hình thức kỷ luật theo qui định Hà Nội, ngày tháng năm 2014 Học viên Trịnh Viết Thiệu Học viên: Trịnh Viết Thiệu ii DANH MỤC CÁC KÍ HIỆU Các ký hiệu sử dụng luận văn A – Diện tích, m2 Ta - Nhiệt độ môi trường xung quanh bình COP - Hệ số bơm nhiệt chứa, oC C- Nhiệt dung riêng đẳng áp, J/kg.K F – Hệ số hấp thụ nhiệt, diện tích trao đổi nhiệt, m T - Nhiệt độ nước, oC fsol - Hệ số đáp ứng Vdisp - Thể tích quét pittong, m3/phút HB - Tổng trực xạ, W/m2  -Thời gian, s Ho - Bức xạ mặt trời đến mặt phẳng nằm ngang khí quyển, W/m2  - Gia số thời gian, s Hs - Tán xạ lên mặt phẳng, W/m2 () - Tích số truyền dẫn hấp thụ Hr - Phản xạ, W/m U – Hệ số tổn thất, W/m2.K p1 - Hiệu áp suất thuỷ tĩnh thu i - Số thứ tự phân tố, Enthalpy, kJ/kg lượng mặt trời, Pa IT – Cường độ xạ mặt phẳng nghiêng, W/m2 p2 - Hiệu áp suất thuỷ tĩnh đường ống k - Hệ số truyền nhiệt, W/m.K kết nối, Pa M - Khối lượng phân tố, kg pf - Tổn thất áp suất ma sát đường m - Lưu lượng khối lượng, kg/s ống thu lượng mặt trời Nr - Công suất nén thực, W  Hiệu suất thu phẳng N el - Công suất điện tiêu thụ, W 1 - Thể tích riêng môi chất đầu vào máy Ns - Nén đoạn nhiệt, W lt - Chiều cao phần trụ đứng vỏ bình chứa, m Q - Nhiệt lượng, W nén, m3/kg me - Hiệu suất el - Hiệu suất điện máy nén QLOAD Lượng nhiệt cần cung cấp cho nhu cầu sử dụng, kJ/ngày Học viên: Trịnh Viết Thiệu iii Các ký hiệu chân a– Không khí r – Ra c – Dàn ngưng, thu tth – toàn phần e – Dàn bay v – Vào hi - hữu ích w – Nước l - Tổn thất i – bên n - Nhận o – bên Học viên: Trịnh Viết Thiệu iv DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1 Sơ đồ hệ thống với môi chất nhận nhiệt trực tiếp từ thu NLMT Hình 1.2 Sơ đồ hệ thống với môi chất nhận nhiệt gián tiếp từ thu NLMT Hình 1.3 Sơ đồ BN kết hợp với NLMT với bình chứa nước nóng Hình 1.4 Sơ đồ BN kết hợp với NLMT với hai bình chứa Hình 1.5 Sơ đồ hệ thống nước nóng lượng mặt trời tuần hoàn tự nhiên .9 Hình 1.6 Phân tầng nhiệt độ bình chứa nước nóng 10 Hình 2.1 Đồ thị lgp-i chu trình bơm nhiệt .20 Hình 2.2 Chia lớp bình ngưng ống xoắn ruột gà theo chiều cao bình chứa 26 Hình 2.3 Cân nhiệt cho phân tố thứ i thời điểm  27 Hình 2.4 Sơ đồ hệ thống lượng mặt trời tuần hoàn tự nhiên 32 Hình 2.5 Quan hệ góc hình học tia xạ mặt trời mặt phẳng nghiêng 36 Hình 2.6 Sơ đồ phân bố thành phần xạ khuếch tán 38 Hình 2.7 Thành phần xạ mặt phẳng nghiêng 39 Hình 2.8 Kết nối TRNSYS với EES .48 Hình 3.1 Mô hình bơm nhiệt sản xuất nước nóng 49 Hn Lưu đồ thuật toán mô bơm nhiệt sản xuất nước nóng 50 Hình 3.3 Sơ đồ hệ thống bơm nhiệt với bình tích nhiệt 51 Hình 3.4 Mô hình bơm nhiệt với bình tích nhiệt 51 Hn Lưu đồ thuật toán mô bơm nhiệt chế độ không ổn định 52 Hình 3.6 Mô hình lượng mặt trời tuần hoàn tự nhiên 54 Hình 3.7 Lưu đồ thuật toán mô tuần hoàn tự nhiên thu lượng mặt trời .55 Học viên: Trịnh Viết Thiệu v Hình 3.8 Giao diện phần mềm mô bơm nhiệt sản xuất nước nóng 58 Hình 3.9 Giao diện phần mềm mô bơm nhiệt với bình chứa tích nhiệt .59 Hình 3.10 Sơ đồ kết nối phận mô hình thu lượng mặt trời 60 Hình 3.11 Sơ đồ kết nối mô-đun mô thu lượng mặt trời TRNSYS 61 Hình 3.12 Sơ đồ kết nối phận mô hình BN kết hợp với NLMT 62 Hình 3.13 Sơ đồ kết nối mô-đun mô xác định lượng nhiệt cần bù cho hệ thống TRNSYS 63 Hình 4.1 Phân tầng nhiệt độ ống tăng cường theo thời gian .68 Hình 4.2 Phân tầng nhiệt độ điểm ống tăng cường theo thời gian 69 Hình 4.3 Sự thay đổi áp suất ngưng tụ, bay theo thời gian .70 Hình 4.4 Sự thay đổi suất gia nhiệt theo thời gian 70 Hình 4.5 Sự thay đổi COP theo thời gian 71 Hình 4.6 Sai số nhiệt độ phân tố cao ống tăng cường T2 .72 Hình 4.7 Sai số nhiệt độ phân tố cao ống tăng cường T10 .72 Hình 4.8 Kết thực nghiệm phân tầng nhiệt độ theo thời gian 73 Hình 4.9 Kết mô phân tầng nhiệt độ theo thời gian 73 Hình 4.10 Kết đo đạc mô T2 (vị trí cao nhất) 74 Hình 4.11 Kết đo đạc mô T5 (Vị trí thấp nhất) .74 Hình 4.12 Hiệu suất thu phụ thuộc vào hệ số (T5-Ta)/I .75 Hình 4.13 Năng lượng hữu ích thu theo thời gian .76 Hình 4.14 Nhu cầu sử dụng theo .79 Hình 4.15 Cường độ xạ ngày 79 Hình 4.16 Kết mô hệ thống BN kết hợp với NLMT .79 Học viên: Trịnh Viết Thiệu vi DANH MỤC CÁC BẢNG Bảng 4.1 Bảng thông số cấu tạo thu phẳng 65 Bảng 4.2 Thống kê đặc tính kỹ thuật thiết bị thí nghiệm 67 Bảng 4.3 Thông số mô tả mô hình bơm nhiệt với bình tích nhiệt .67 Bảng PL1.1 Bảng kết mô phân tầng nhiệt độ BN với bình tích nhiệt 86 Bảng PL1.2 Bảng kết mô BN với bình tích nhiệt .87 Bảng PL1.3 ảng kết thực nghiệm phân tầng nhiệt độ 89 Bảng PL2.1 Bảng kết thực nghiệm NLMT tuần hoàn tự nhiên 90 Bảng PL2.2 Bảng kết mô NLMT tuần hoàn tự nhiên 92 Học viên: Trịnh Viết Thiệu vii MỤC LỤC LỜI CẢM ƠN i LỜI CAM ĐOAN ii DANH MỤC CÁC KÍ HIỆU iii DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ v DANH MỤC CÁC BẢNG vii MỤC LỤC viii MỞ ĐẦU CHƢƠNG - TỔNG QUAN VỀ CÁC PHƢƠNG PHÁP KẾT HỢP BƠM NHIỆT VỚI BỘ THU NĂNG LƢỢNG MẶT TRỜI SẢN XUẤT NƢỚC NÓNG VÀ MÔ PHỎNG HỆ THỐNG 1.1 Các phương pháp kết hợp bơm nhiệt với thu lượng mặt trời sản xuất nước nóng 1.1.1 Sử dụng thu lượng mặt trời để tăng nhiệt độ nguồn lạnh 1.1.2 Sử dụng thu lượng mặt trời sản xuất nước nóng với bơm nhiệt nguồn nhiệt dự phòng 1.2 Các yếu tố ảnh hưởng vai trò mô hệ thống bơm nhiệt kết hợp với thu lượng mặt trời sản xuất nước nóng 1.2.1 Các yếu tố ảnh hưởng đến hiệu lượng hệ thống nước nóng 1.2.2 Vai trò mô hệ thống bơm nhiệt kết hợp với thu lượng mặt trời sản xuất nước nóng 12 1.3 Tình hình nghiên cứu nước liên quan tới đề tài .13 1.3.1 Nghiên cứu nước 13 1.3.2 Nghiên cứu nước .13 1.4 Mục tiêu luận văn 18 Học viên: Trịnh Viết Thiệu viii CHƢƠNG - CƠ SỞ LÝ THUYẾT MÔ PHỎNG HỆ THỐNG BƠM NHIỆT KẾT HỢP VỚI NĂNG LƢỢNG MẶT TRỜI SẢN XUẤT NƢỚC NÓNG 20 2.1 Xây dựng mô hình toán học bơm nhiệt sản xuất nước nóng 20 2.1.1 Mô hình toán học máy nén pittong .21 2.1.2 Mô hình toán học dàn ngưng dàn bay 23 2.1.3 Mô hình toán học thiết bị tiết lưu 25 2.1.4 Hệ số bơm nhiệt 25 2.2 Xây dựng mô hình toán học phân tầng nhiệt độ nước bình ngưng chế độ tính nhiệt 25 2.3 Xây dựng mô hình toán học tuần hoàn tự nhiên thu lượng mặt trời 31 2.4 Các mô hình khác TRNSYS sử dụng để mô 34 2.4.1 Mô hình tính toán xạ mặt trời mặt phẳng nghiêng (Type 16a) .34 2.4.2 Mô hình xác định lượng hữu ích thu phẳng (Type 1a) 40 2.4.3 Mô hình khác .47 CHƢƠNG - MÔ PHỎNG HỆ THỐNG SẢN XUẤT NƢỚC NÓNG 49 3.1 Xây dựng bổ sung mô-đun mô hệ thống sản xuất nước nóng .49 3.1.1 Mô-đun mô bơm nhiệt sản xuất nước nóng .49 3.1.2 Mô-đun mô bơm nhiệt với bình tích nhiệt .51 3.1.3 Mô-đun mô thu lượng mặt trời tuần hoàn tự nhiên 54 3.2 Lựa chọn ngôn ngữ lập trình 55 3.2.1 Giới thiệu EES 55 3.2.2 Giới thiệu TRNSYS .56 3.3 Xây dựng chương trình phần mềm mô hệ thống sản xuất nước nóng 57 Học viên: Trịnh Viết Thiệu ix 3.3.1 Xây dựng phần mềm mô bơm nhiệt sản xuất nước nóng 57 3.3.2 Xây dựng phần mềm mô bơm nhiệt bình chứa tích nhiệt 58 3.3.3 Xây dựng phần mềm mô thu lượng mặt trời tuần hoàn tự nhiên 59 3.3.4 Xây dựng phần mềm mô bơm nhiệt kết hợp với thu lượng mặt trời 61 CHƢƠNG – KẾT QUẢ VÀ THẢO LUẬN .64 4.1 Mô tả thí nghiệm kiểm chứng 64 4.1.1 Hệ thống lượng mặt trời .64 4.1.2 Hệ thống bơm nhiệt sản xuất nước nóng 65 4.2 Đánh giá kết thu .67 4.2.1 So sánh với kết thực nghiệm hệ thống bơm nhiệt sản xuất nước nóng .67 4.2.2 So sánh với kết thực nghiệm lượng mặt trời sản xuất nước nóng .72 4.3 Ứng dụng phần mềm mô hệ thống sản xuất nước nóng .76 CHƢƠNG - KẾT LUẬN VÀ ĐỀ XUẤT 80 TÀI LIỆU THAM KHẢO 82 PHỤ LỤC 86 PL1 Kết bơm nhiệt với bình tích nhiệt 86 PL2 Kết NLMT tuần hoàn tự nhiên 90 PL3 Bảng kết BN kết hợp NLMT sản xuất nước nóng 94 PL4 Code chương trình 101 Học viên: Trịnh Viết Thiệu x Tiêu thụ nƣớc nóng t eo ngƣời ngày PL4 Code c ƣơng tr n Bơm n iệt sản xuất nƣớc nóng Procedure HP(m_a,m_w,t_a, t_wv:t_c,t_e,Q_e,Q_c,COP) {Mo phong bom nhiet} k_c=939 {w/m2.K} F_c=0.216{[m2]} C_p=4180 t_c=t_wv+5 k_e=439 F_e=0.205 C_Pair=1005 V_swept=10.29{[cm^3/rev]} m=2900 {[rev/min]} R$='R134a' t_e=t_a-10 repeat repeat p_c=PRESSURE(R$,T=t_c,x=0) i_3=ENTHALPY(R$,P=p_c,x=0) i_4=i_3 eta_m=0.8{[-], mechanical efficiency} eta_e=0.85{[-], electrical efficiency} eta_v=0.9011 Học viên: Trịnh Viết Thiệu 101 eta_s=0.84636 p_e=Pressure(R$,T=t_e,x=1) i_1=Enthalpy(R$,x=1,P=p_e) s_1=Entropy(R$,T=t_e,x=1) v_1=Volume(R$,T=t_e,x=1) s_2s=s_1 {V_disp=32.5 {[m^3/h]} V_disp_dot=V_disp/3600 {[m^3/s]} G_r=V_disp_dot*eta_v/v_1 {[kg/s]}} G_r=eta_v*V_swept*m/(60*v_1*10^6){[kg/s]} i_2s=Enthalpy(R$,s=s_1,P=p_c) t_2s=Temperature(R$,p=p_c,s=s_2s) L_s_dot=G_r*(i_2s-i_1) {[kW]} L_r_dot=L_s_dot/eta_s{[kW]} L_e_dot=L_r_dot/eta_m/eta_e {[kW]} i_2r=i_1+L_r_dot/G_r {[kJ/kg]} t_2r=Temperature(R$,p=p_c,h=i_2r){[oC]} Q_e=G_r*(i_1-i_3) T_outair=T_a-Q_e*1000/(m_a*C_Pair) Hs4=k_e*F_e/(m_a*C_Pair) Hs5=exp(Hs4) t_e1=(Hs5*t_outair-t_a)/(Hs5-1) T_e2=T_e T_e=T_e1 Until (abs(T_e2-T_e)n+1) and (i10) m_dot_1=m_dot m_dot=m_dot_cal a=m_dot*Delta_Tau/M_out b=Delta_Tau*k_l*l_tank*T_air/(M_out*C_w*(n-1)) c=m_dot*Delta_Tau/M_in d=Delta_Tau*Q_dot/(M_in*C_w*(n-1)) "2"until (abs( m_dot_1-m_dot )< 0.01) k_c=0.643 F_c= 0.215{[m2]} t_e=t_air-12 eta_m=0.8{[-], mechanical efficiency} eta_e=0.85{[-], electrical efficiency} V_swept=10.29{[cm^3/rev]} m=2900 {[rev/min]} m_air=0.2 C_Pair=1.005 k_e=0.1 F_e=0.4 Repeat Hs1=m_dot*C_w*(t[10]-t[n])/(k_c*F_c) Hs2=(t[10]-t[n])/Hs1 Hs3=exp(Hs2) t_c=(Hs3*t[10]-t[n])/(Hs3-1) R$='R134a' p_c=Pressure(R$,T=t_c,x=1) p_e=Pressure(R$,T=t_e,x=1) i_1=Enthalpy(R$,x=1,P=p_e) s_1=Entropy(R$,T=t_1,P=p_1) v_1=Volume(R$,T=t_1,P=p_1) s_2s=s_1 eta_v=0.4 {0.9011+0.4303/p_e-1.1129/p_e^2-0.10761*ln(p_c){tim ham}} eta_s=0.4 {0.64636+0.016857*p_e-0.00503*p_e^2-0.1279/p_c{tim ham}} G_r=eta_v*V_swept*m/(60*v_1*10^6){[kg/s]} i_2s=Enthalpy(R$,s=s_1,P=p_c) t_2s=Temperature(R$,p=p_c,s=s_2s) L_s_dot=G_r*(i_2s-i_1) {[W]} L_r_dot=L_s_dot/eta_s{[W]} L_e_dot=L_r_dot/eta_m/eta_e {[W]} i_2r=i_1+L_r_dot/G_r {[J/kg]} t_2r=Temperature(R$,p=p_c,h=i_2r){[oC]} t_3= t_c i_3=Enthalpy(R$,x=0,P=p_c) Q_dot_cal=G_r*(i_2r-i_3) Q_e=G_r*(i_1-i_3) T_outair=T_air-Q_e/(m_air*C_Pair) Hs4=m_air*C_Pair/(k_e*F_e) Hs5=exp(Hs4) t_e1=(Hs5*t_outair-t_air)/(Hs5-1) T_e2=T_e T_e=T_e1 Until (abs(T_e2-T_e)Tau) End {$import 'CLIPBOARD' Tau, Delta_Tau,n,D_outside,D_out_insul,D_inside,l_tank,rho_w,C_w,P_air,T_air ,alpha_w,alpha_air,lambda_insul ,g, K,f_inside,f_outside} {Tau=50} Delta_Tau=10 n=5 D_outside = 0.318 "m" D_out_insul = 0.348 "m" D_inside = 0.09 "m" l_tank = 0.6 "m" rho_w = 1000 "kg/m^3" C_w= 4.180 "J/kg.k" P_air = "bar" T_air = 24 "^0 C" alpha_w = 25/1000 "W/m^2 K" alpha_air = 10/1000 "W/m^2 K" lambda_insul = 0.05/1000 "W/m.K" g=9.81 K=1.5 f_inside=1.8 f_outside=1.9 Call Stratification(lambda_insul,D_out_insul,alpha_air,alpha_w,rho_w,C_w,D_inside,D_outside ,K,f_inside,f_outside,g,l_tank,P_air,Delta_Tau,n,T_air,Tau:t [1 10],m_dot,Q_dot_cal,t_c,p_c,t_e,p_e,COP) {$Export 'CLIPBOARD'M_in, M_out , k_l} Mô-đun NLMT tuần hoàn tự nhiên {I=400 T_air=25 T_in=27 t2=25 t3=25 t4=25 t5=25} $Import 'CLIPBOARD' T_air,T_in,t2,t3,t4,t5,I q=5 D_outside = 0.348 "m" D_inside = 0.318 "m" l_tank = 0.58 "m" rho_w = 1000 "kg/m^3" alpha_w = 50 "W/m^2 K" alpha_air = 15 "W/m^2 K" lambda_insul = 0.05 "W/m.K" M_in= rho_w*pi*D_inside^2*l_tank/(4*(q-1)) k_l= 1/(1/alpha_w + 1/alpha_air +(D_outside -D_inside)/(2*lambda_insul)) A_c=0.842 {[m^2]} d_c=0.014 {[m]} N=10 T_o=t5{[oC]} L_c=1.5{[m]} Phi=20 alpha=0.912 C_p=4180{J/kg-K]} Học viên: Trịnh Viết Thiệu 105 g=9.81 {u_wind=3 {[m/s]} h=(2.8+3*u_wind) {[W/m^2]}} DELTAtau=300 Gamma=(1/(0.5155+0.0192*T_o)-0.12)*10^(-4) Beta=(0.3+0.116*T_o-0.0004*T_o^2)*10^(-4) rho_o=1001-0.08832*T_o-0.003417*T_o^2{DENSITY(R718,T=T_o,P=1) {Note: Khong dung H2O thay cho R718}} a1=256*gamma*L_c/(pi*N*d_c^4) b1=rho_o*g*Beta*(L_c*sin(PHI*pi/180)/2+l_tank+0.15) m_dot=b1*abs(T_in-T_o)/a1 a=DELTAtau*m_dot/M_in b=DELTAtau*k_l*l_tank*T_air/(M_in*C_p*(q-1)) t1=t_in {Phan tang nhiet do} t1new=t_in t2new=(a*T_in+(1-a)*t2+b)/(1+b/T_air) t3new=(a*t2new+(1-a)*t3+b)/(1+b/T_air) t4new=(a*t3new+(1-a)*t4+b)/(1+b/T_air) t5new=(a*t4new+(1-a)*t5+b)/(1+b/T_air) {Hieu suat bo thu} eta=0.7-4.17*(T_o-T_air)/I {Ton that ap suat} {d_rp=14/1000{[mm/m]} d_sp=14/1000{[mm/m]} L_rp=1.4{[m]} L_sp=1.2{[m]} T_m=(t1+T5)/2 rho_m=DENSITY(R718,T=T_m,P=1) {Note: Khong dung H2O thay cho R718} rho1=1001-0.08832*t1-0.003417*t1^2{DENSITY(R718,T=T_o,P=1) {Note: Khong dung H2O thay cho R718}} rho5=1001-0.08832*t5-0.003417*t5^2{DENSITY(R718,T=T_o,P=1) {Note: Khong dung H2O thay cho R718}} Gamma1=(1/(0.5155+0.0192*t1)-0.12)*10^(-6) u=(4*m_dot)/(N*rho_m*pi*d_c^2) f=64*gamma/(d_c*u) DELTAP_c=(1/2)*rho_m*f*u^2*L_c/(2*d_c){[N/m^2]} {Fiction pressure drop in the collector tubes} u_sp=(4*m_dot)/(rho5*pi*d_sp^2) f_sp=64*gamma/(d_sp*u_sp) DELTAP_sp=(1/2)*rho1*f_sp*u_sp^2*L_sp/(2*d_sp){[N/m^2]} {pressure drop in the supply pipe} u_rp=(4*m_dot)/(rho1*pi*d_rp^2) f_rp=64*gamma1/(d_rp*u_rp) DELTAP_rp=(1/2)*rho5*f_rp*u_rp^2*L_rp/(2*d_rp){[N/m^2]} {pressure drop in the return pipe}} out[1]=t2new out[2]=t3new out[3]=t4new out[4]=t5new out[5]=m_dot out[6]=eta $export 'CLIPBOARD' out[1],out[2],out[3],out[4],out[5],out[6] NLMT tuần hoàn tự nhiên VERSION 15 ******************************************************************************* *** TRNSYS input file (deck) generated by IISiBat *** on Th? Hai, Tháng Ba 03, 2014 at 13:30 *** from IISiBat project: Học viên: Trịnh Viết Thiệu 106 *** *** If you edit this file, use the File/Import TRNSYS Input File function in *** IISiBat to update the project *** *** If you have problems, questions or suggestions please contact your local *** TRNSYS distributor or mailto:iisibat@cstb.fr *** ******************************************************************************* ASSIGN "Project1.LST" ******************************************************************************* *** Control cards ******************************************************************************* * START, STOP and STEP CONSTANTS START=1 STOP=24 STEP=1 *SIMULATION Start time End time Time step SIMULATION START STOP STEP * User defined CONSTANTS * TOLERANCES 0.001 0.001 * LIMITS 30 30 30 * DFQ * WIDTH 80 * LIST * MAP * SOLVER Integration Convergence Max iterations Max warnings Trace limit TRNSYS numerical integration solver method TRNSYS output file width, number of characters NOLIST statement MAP statement Solver statement ******************************************************************************* *** Units ******************************************************************************* * Model "TYPE1a" (Type 1) * UNIT TYPE TYPE1a *$UNIT_NAME TYPE1a *$MODEL \Solar Thermal Collectors\Quadratic Efficiency Collector\No Incidence Angle Modification\TYPE1a.tmf *$POSITION 288 178 *$LAYER Main PARAMETERS * Number in series * Collector area Học viên: Trịnh Viết Thiệu 107 2.0 * Fluid specific heat 4.19 * Efficiency mode * Tested flow rate 50.0000000000001 * Intercept efficiency 0.7 * Efficiency slope 15 * Efficiency curvature * Optical Mode 1 INPUTS * [unconnected] Inlet temperature 0,0 * [unconnected] Inlet flowrate 0,0 * TYPE9a:Output ->Ambient temperature 2,2 * TYPE9a:Output ->Incident radiation 2,1 *** INITIAL INPUT VALUES 30 100 10.0 * -* Model "TYPE9a" (Type 9) * UNIT TYPE TYPE9a *$UNIT_NAME TYPE9a *$MODEL \Utility\Weather and Other Data Readers\All Other Data Files\1st Line = Simulation Start\Free Format Read\TYPE9a.tmf *$POSITION 149 242 *$LAYER Weather / Data Files PARAMETERS 12 * Mode * Header Lines to Skip * No of values to read * Time interval of data 1.0 * Interpolate or not?-1 * Multiplication factor-1 1.0 * Addition factor-1 * Interpolate or not?-2 * Multiplication factor-2 1.0 * 10 Addition factor-2 * 11 Logical unit Học viên: Trịnh Viết Thiệu 108 14 * 12 Not used -1 *** External files ASSIGN "C:\Users\Administrator\Desktop\EX1DAT.DAT" 14 *|? Which file contains the data to be read by this component? |1000 * -* Model "TYPE25c" (Type 25) * UNIT TYPE 25 TYPE25c *$UNIT_NAME TYPE25c *$MODEL \Output\Printer\Don't Print Units to File\TYPE25c.tmf *$POSITION 472 189 *$LAYER Outputs PARAMETERS * Printing interval * Start time * Stop time 8760 * Logical unit 17 INPUTS * TYPE1a:Outlet temperature ->Input to be printed-1 1,1 * TYPE9a:Output ->Input to be printed-2 2,1 * TYPE9a:Output ->Input to be printed-3 2,2 *** INITIAL INPUT VALUES Outlet Output Output *** External files ASSIGN "EX1" "" 17 *|? Which file should contain the printed results? |1000 * -END Bơm n iệt kết hợp với t u lƣợng mặt trời sản xuất nƣớc nóng VERSION 15 ******************************************************************************* *** TRNSYS input file (deck) generated by IISiBat *** on Th? Hai, Tháng Tý 28, 2014 at 22:59 *** from IISiBat project: C:\Users\Administrator\Desktop\BNNNMT\BNNLMT.TPF *** *** If you edit this file, use the File/Import TRNSYS Input File function in *** IISiBat to update the project *** *** If you have problems, questions or suggestions please contact your local *** TRNSYS distributor or mailto:iisibat@cstb.fr *** ******************************************************************************* ASSIGN "Project1.LST" Học viên: Trịnh Viết Thiệu 109 ******************************************************************************* *** Control cards ******************************************************************************* * START, STOP and STEP CONSTANTS START=0 STOP=24 STEP=0.0833323333533329 *SIMULATION Start time End time Time step SIMULATION START STOP STEP * User defined CONSTANTS * Integration Convergence TOLERANCES 0.001 0.001 * Max iterations Max warnings Trace limit LIMITS 30 30 30 * TRNSYS numerical integration solver method DFQ * TRNSYS output file width, number of characters WIDTH 80 * NOLIST statement LIST * MAP statement MAP * Solver statement SOLVER ******************************************************************************* *** Units ******************************************************************************* * Model "TYPE9a" (Type 9) * UNIT TYPE TYPE9a *$UNIT_NAME TYPE9a *$MODEL \Utility\Weather and Other Data Readers\All Other Data Files\1st Line = Simulation Start\Free Format Read\TYPE9a.tmf *$POSITION 288 189 *$LAYER Weather / Data Files PARAMETERS 12 * Mode * Header Lines to Skip * No of values to read * Time interval of data 0.0833323333533329 * Interpolate or not?-1 * Multiplication factor-1 1.0 * Addition factor-1 * Interpolate or not?-2 * Multiplication factor-2 1.0 * 10 Addition factor-2 * 11 Logical unit 14 Học viên: Trịnh Viết Thiệu 110 * 12 Not used -1 *** External files ASSIGN "C:\Users\Administrator\Desktop\BNNNMT\EX11DAT.DAT" 14 *|? Which file contains the data to be read by this component? |1000 * -* Model "TYPE1a" (Type 1) * UNIT TYPE TYPE1a *$UNIT_NAME TYPE1a *$MODEL \Solar Thermal Collectors\Quadratic Efficiency Collector\No Incidence Angle Modification\TYPE1a.tmf *$POSITION 363 413 *$LAYER Main PARAMETERS * Number in series * Collector area * Fluid specific heat 4.19 * Efficiency mode * Tested flow rate 50.0000000000001 * Intercept efficiency 0.7 * Efficiency slope 15 * Efficiency curvature * Optical Mode 1 INPUTS * [unconnected] Inlet temperature 0,0 * NDF:Output-1 ->Inlet flowrate 7,1 * TYPE9a:Output ->Ambient temperature 1,2 * TYPE16a:Total radiation on surface ->Incident radiation 5,7 *** INITIAL INPUT VALUES 21 108 21 * -* Model "TYPE65" (Type 65) * UNIT TYPE 65 TYPE65 *$UNIT_NAME TYPE65 *$MODEL \Output\Online Plotter\TYPE65.tmf *$POSITION 832 338 *$LAYER Outputs PARAMETERS 10 * Nb of left-axis variables * Nb of right-axis variables * Left axis minimum 20 Học viên: Trịnh Viết Thiệu 111 * Left axis maximum 70 * Right axis minimum * Right axis maximum 25000 * Number of plots per simulation * X-axis gridpoints * Shut off Online w/o removing * 10 Logical Unit for ouput file -1 INPUTS * TYPE1a:Outlet temperature ->Left axis variable 2,1 * TYPE24:Result of integration-3 ->Right axis variable-1 8,3 * TYPE24:Result of integration-2 ->Right axis variable-2 8,2 *** INITIAL INPUT VALUES Tout Qhi Q_load LABELS C kW Temperature Power Result * -* Model "TYPE25c" (Type 25) * UNIT TYPE 25 TYPE25c *$UNIT_NAME TYPE25c *$MODEL \Output\Printer\Don't Print Units to File\TYPE25c.tmf *$POSITION 835 104 *$LAYER Outputs PARAMETERS * Printing interval 0.0833323333533329 * Start time * Stop time 19 * Logical unit 17 INPUTS * load_T14h:Average value of function ->Input to be printed-1 6,1 * NDF:Output-1 ->Input to be printed-2 7,1 * TYPE24:Result of integration-2 ->Input to be printed-3 8,2 * TYPE24:Result of integration-3 ->Input to be printed-4 8,3 * TYPE1a:Outlet temperature ->Input to be printed-5 2,1 *** INITIAL INPUT VALUES Average m_dot Q_load Q_hi Outlet *** External files Học viên: Trịnh Viết Thiệu 112 ASSIGN "EX1" "" 17 *|? Which file should contain the printed results? |1000 * -* Model "TYPE16a" (Type 16) * UNIT TYPE 16 TYPE16a *$UNIT_NAME TYPE16a *$MODEL \Physical Phenomena\Radiation Processors\Total Horiz Only Known (Mode=1)\No Radiation Smoothing\TYPE16a.tmf *$POSITION 366 104 *$LAYER Weather / Data Files *$# *$# *$# PARAMETERS * Horiz radiation mode * Tracking mode * Tilted surface mode * Starting day 200 * Latitude 21.5 * Solar constant 4871.00000000002 * Shift in solar time 0.0 * Not used * Solar time? -1 INPUTS * TYPE9a:Output ->Radiation on horizontal 1,1 * TYPE9a:Time of last read ->Time of last data read 1,99 * TYPE9a:Time of next read ->Time of next data read 1,100 * [unconnected] Ground reflectance 0,0 * [unconnected] Slope of surface 0,0 * [unconnected] Azimuth of surface 0,0 *** INITIAL INPUT VALUES 0 0.0833323333533329 0.2 20 0.0 * -* Model "load_T14h" (Type 14) * UNIT TYPE 14 load_T14h *$UNIT_NAME load_T14h *$MODEL \Utility\Forcing Functions\General\TYPE14h.tmf *$POSITION 476 456 *$LAYER Main PARAMETERS 28 * Initial value of time Học viên: Trịnh Viết Thiệu 113 * Initial value of function 0.0 * Time at point-1 * Value at point -1 * Time at point-2 * Value at point -2 0.005 * Time at point-3 * Value at point -3 0.005 * Time at point-4 * 10 Value at point -4 * 11 Time at point-5 11 * 12 Value at point -5 * 13 Time at point-6 11 * 14 Value at point -6 0.005 * 15 Time at point-7 12 * 16 Value at point -7 0.005 * 17 Time at point-8 12 * 18 Value at point -8 * 19 Time at point-9 17.5 * 20 Value at point -9 * 21 Time at point-10 17.5 * 22 Value at point -10 0.025 * 23 Time at point-11 19 * 24 Value at point -11 0.025 * 25 Time at point-12 19 * 26 Value at point -12 * 27 Time at point-13 24 * 28 Value at point -13 * -* Model "NDF" (Type 66) UNIT TYPE 66 NDF *$UNIT_NAME NDF Học viên: Trịnh Viết Thiệu 114 *$MODEL \Utility\Calling External Routines\EES Routines\Type66.tmf *$POSITION 620 292 *$LAYER Main INPUTS * load_T14h:Average value of function ->Input-1 6,1 * TYPE9a:Output ->Input-2 1,2 * TYPE1a:Outlet temperature ->Input-3 2,1 * NDF:Output-3 ->Input-4 7,3 * NDF:Output-4 ->Input-5 7,4 * NDF:Output-5 ->Input-6 7,5 * NDF:Output-6 ->Input-7 7,6 * [unconnected] Input-8 0,0 *** INITIAL INPUT VALUES 0.03 21 21 21 21 21 21 60 LABELS 1 c:\ees32\ees.exe C:\Users\Administrator\Desktop\BNNNMT\NLMT1.ees * -* Model "TYPE24" (Type 24) UNIT TYPE 24 TYPE24 *$UNIT_NAME TYPE24 *$MODEL \Utility\Integrators\Quantity Integrator\TYPE24.tmf *$POSITION 555 82 *$LAYER Main PARAMETERS * Reset time 8760.0 INPUTS * TYPE16a:Total radiation on surface ->Input to be integrated-1 5,7 * NDF:Output-7 ->Input to be integrated-2 7,7 * TYPE1a:Useful energy gain ->Input to be integrated-3 2,3 *** INITIAL INPUT VALUES 0.0 0.0 0.0 * -END Học viên: Trịnh Viết Thiệu 115 ... n-ớc nóng (bộ thu l-ợng mặt trời) Van chiều N-ớc lạnh bổ sung Bức xạ mặt trời Bơm tuần hoàn Bình bảo ôn n-ớc nóng bơm nhiệt Bộ thu dạng phẳng N-ớc nóng sử dụng Hệ thống bơm nhiệt sản xuất n-ớc nóng. .. cho nc ton b h thng núng lờn, quỏ trỡnh lu thụng tip tc ti chng no nhit ca nc b thu cũn ln hn nhit nc bỡnh cha Tới tải Bình bảo ôn chứa n-ớc nóng (bộ thu l-ợng mặt trời) N-ớc cấp Bức xạ mặt. .. nhit ca b thu mt tri tng hiu qu thu nhit ca b thu NLMT thỡ bm nhit c s dng cho quỏ trỡnh trao i nhit phớa trờn ca bỡnh, b thu nng lng mt tri tớch nhit phớa di bỡnh Ch hnh thit b BN v b thu NLMT