Kết luận
Sau hơn 1 năm tiến hành nghiên cứu và khảo nghiệm trên mô hình thực tế của đề tài “Nghiên cứu khả năng sản xuất nước sạch từ nước phèn tận dụng nhiệt thải lò hơi”, luận văn đã đạt được những kết quả như sau:
- Luận văn dựa trên các nghiên cứu tổng quan, đề xuất nghiên cứu sản xuất nước sạch từ nước phèn bằng hệ thống chưng cất nước một cấp có bỉnh phân ly, tận dụng tối da nguồn nhiệt thải của lò hơi và nguồn nhiệt nước ngưng.
- Luận văn đã trình bày đầy đủ các công thức cần thiết để tính toán thiết kế hệ thống, từ đó xây dựng chương trình tính bằng ngôn ngữ Matlab.
- Nếu công suất lò hơi đốt trấu là 18 tấn/giờ, thì kết quả sản lượng nước chưng cất thu được là 0,9955m3
- Quá trình khảo sát một số ảnh hưởng của thông số vận hành như nhiệt độ nước cấp, nhiệt độ của khói và sản lượng hơi, đã chỉ ra rằng các yếu tố trên ảnh hưởng khá lớn đến hệ thống chưng cất nước phèn…và đã thực hiện tính toán kinh tế, kết quả cho thấy phương án này có tính khả thi – thời gian thu hồi vốn là 2,9 năm.
/h.
Kiến nghị
- Xây dựng một mô hình chưng cất nước hoàn chỉnh để thực nghiệm và đánh giá hiệu quả sản xuất nước giống thực tế.
- Mở rộng phạm vi ứng dụng đề tài cho các nguồn nhiệt khác.
HV: NGUYỄN VŨ BẰNG Page 65
TÀI LIỆU THAM KHẢO
[1]. Sở Tài Nguyên Môi Trường Đồng Tháp “Báo cáo thuyết minh Quy hoạch quản lý, khai thác và bảo vệ tài nguyên nước dưới đất trên khu vực Đồng Tháp đến năm 2015, định hướng đến năm 2020”.
[2]. Báo cáo hội thảo 30 năm khai thác Khu vực Đồng Tháp Mười.
[3]. Tạp chí Tổng cục thống kê VN
[4]. Cục quản lý tài nguyên nư ớc - Bộ Tài Nguyên và Môi trường - Báo cáo nhu cầu nước.
[5]. PGS TS. Đào Ngọc Chân, PGS TS Hoàng Ngọc Đồng - Lò hơi và thiết bị chất đốt –– NXB Khoa học Kỹ thuật.
[6]. Cục Thống kê Việt Nam -Tỷ lê tăng dân số Việt Nam 2014 –.
[7] QD/1208-TTG “Phê duyệt Quy hoạch phát triển điện lực quốc gia giai đoạn 2011 - 2020 có xét đến năm 2030”
[8] Nguyễn Thế Bảo “ Bài giảng bảo toàn và quản lý năng lượng”-2014
[9] Hoàng Đình Tín, Lê Chí Hi ệp – Nhiệt động lực học kỹ thuật – NXB ĐHQG TpHCM – 2011 (Tái bản lần thứ nhất) “trang 13, 369, 387”.
[10] Thông báo của Bộ Xây Dựng về dự án ‘tăng cường năng lực quản lý nhà nư ớc về tiết kiệm năng lượng cho các cơ sở công nghiệp’ do Cơ quan Quản lý Môi trư ờng và Năng lượng Pháp tài trợ thử nghiệm đánh giá tiềm năng thực hiện việc giảm lượng giảm phát thải CO2.
[11] TS Nguyễn Văn Tuyên “Bài giảng trung tâm nhà máy nhiệt điện”-2012.
[12] Ths. Lâm Vĩnh Sơn – “Bài giảng xử lý nước thải” – 2010.
[13] Viện môi trường và tài nguyên – Trường ĐH Quốc Gia Tp.HCM “ báo cáo hiện trạng môi trường 05 năm tỉnh Đồng Tháp” 2010.
[14] Ths. Đỗ Xuân Hoàn - “Đồ án Tính toán thiết kế hệ thống xử lý nước nhiễm phèn khu vực đồng bằng sông Cửu Long công suất 1000m3_ngày đêm” – 2008.
[15] GS.TS Phạm Xuân Vượng “ Giáo trình kỹ thuật Lò hơi” Trường Đại Học Nông Nghiệp Hà Nội – 2007.
[16] TS. Nguyễn Thanh Quáng “ Nghiên cứu sản xuất viên nhiên liệu từ Biomass” báo cáo hội nghị sinh viên khoa học lần thứ 6 – Trường Đại học Bách Khoa Đà Nẵng.
[17] PGS.TS Đào Ngọc Chân, PGS TS Hoàng Ngọc Đồng – Lò hơi và thiết bị đốt – Nhà xuất bản khoa học Kỹ thuật – 2008.
[18] PSG.TS Phạm Lê Dzần, TS Nguyễn Công Hân – Công nghệ lò hơi và mạng nhiệt - Nhà xuất bản khoa học Kỹ thuật – 2008.
[19] TS Nguyễn Thanh Hào – Nhà máy nhiệt điện – Nhà xuất bản Đại học Công nghiệp Tp.HCM, 2011.
HV: NGUYỄN VŨ BẰNG Page 66
[20] TS Nguyễn Sĩ Mão “Lý thuyết và thiết bị cháy” - Nhà xuất bản khoa học Kỹ thuật – 2002.
HV: NGUYỄN VŨ BẰNG Page 67
PHỤ LỤC Phụ lục 1
Code phần mền tính toán “Nghiên cứu khả năng sản xuất nước sạch từ nước phèn tận dụng nhiệt khói thải lò hơi”
C_lv=str2num(get(handles.C_lv,'string')); % thanh phan cacbon lam viec (%) H_lv=str2num(get(handles.H_lv,'string')); % thanh phan hidro lam viec (%) O_lv=str2num(get(handles.O_lv,'string')); % thanh phan oxi lam viec
S_lv=str2num(get(handles.S_lv,'string')); % thanh phan luu huynh lam viec N_lv=str2num(get(handles.N_lv,'string')); % thanh phan nito lam viec A_lv=str2num(get(handles.A_lv,'string')); % thanh phan tro
W_lv=str2num(get(handles.W_lv,'string')); % thanh phan am lam viec
% cac gia tri dau vao cua he thong lo hoi
p_lh=str2num(get(handles.p_lh,'string')); % ap suat hoi lam viec (bar) p_ilh=[5 6 7 8 9 10 11 12 13 14 15]; %gan bien ap suat lo hoi (bar)
i_ilh=[2749 2757 2763 2769 2774 2778 2782 2785 2788 2790 2792]; %gan bien the entanpi theo ap suat soi cua lo hpo (kj/kg)
p= polyfit(p_ilh,i_ilh,9); % tim phuong trinh bac 2 di qua qua p_ilh,i_ilh i_soi=polyval(p,p_lh) % entanpi cua hoi tuong ung voi ap suat soi (kj/kg) t_nclh=30; % nhiet do nuoc cap lo hoi khi gia dinh chua co bo ham nuoc cap (C)
t_inclh=[28 30 32 34 36 38 40 42 44 46 48 50]; % gan bien nhiet do nuoc cap lo hoi(C)
i_inclh=[117.4 125.8 134.1 142.5 150.8 159.2 167.5 175.9 184.2 192.6 201 209.3]; %gan bien the entanpi theo ap suat soi cua lo hpo (kj/kg)
i= polyfit(t_inclh,i_inclh,10); % tim phuong trinh bac 2 di qua qua t_inclh,i_inclh
i_nclh=polyval(i,t_nclh) % entanpi cua nuoc cap lo hoi (kj/kg)
% cac gia tri dau ra cua he thong lo hoi
D_h=str2num(get(handles.D_h,'string')); %san luong hoi cua lo hoi (kg/h)
% cac gia tri ton that cua he thong lo hoi
% ta hoan toan co the tinh toan duoc cac ton that nhiet cua toan lo hoi nh?
Q_2, Q_3, Q_4, Q_5 và Q_6
% nhung de tinh toan nhanh, ta co the dua vao cac tai lieu ky thuat va kinh nghiem de chon hieu suat nhiet toan phan cua lo hoi
hieusuat_lh=str2num(get(handles.hieusuat_lh,'string')); % hieu suat nhiet toan phan cua lo hoi
alpha_kk=str2num(get(handles.alpha_kk,'string')); % he so khong khi thua chon theo buong dot va nhien lieu
% tinh toan cac thong so dau vao va ra cua lo hoi
% Nhiet tri thap lam viec cua nhien lieu
Qt_lv=(339*37.13)+(1030*4.12)-(109*(31.6-0.05))-(25*9) % nhiet tri thap lam viec cua nhien lieu (kj/kg)
% l_ong??ng tieu hao nhien lieu cho lo hoi
B=D_h*(i_soi-i_nclh)/(Qt_lv*hieusuat_lh) % luong tieu hao nhien lieu cho lo hoi (kg/h)
% the tich san pham chay (khoi)
V0_oxi=22.4/100*((C_lv/12)+(H_lv/4)+(S_lv/32)-(O_lv/32))% the tich oxi ly thuyet (m3tc/kg)
V0_kk=100/21*V0_oxi % luong khong khi ly thuyet (m3tc/kg)
V0_nito=((N_lv/100)*(1/1.251))+(0.79*V0_kk) % the tich ly thuyet khi nito (m3tc/kg)
V0_RO2=(1.866*(C_lv/100))+(0.7*(S_lv/100)) % tong the tich ly thuyet cua CO2 va SO2 (m3tc/kg)
V0_H2O=(0.112*H_lv)+(0.0124*W_lv)+(0.0161*V0_kk) % the tich ly thuyet cua hoi nuoc (m3tc/kg)
V0_spc=V0_RO2+V0_nito+V0_H2O % the tich san pham chay ly thuyet (m3tc/kg)
HV: NGUYỄN VŨ BẰNG Page 68 V_spc=V0_spc+((alpha_kk-1)*V0_kk) % the tich khoi thuc te tinh cho 1kg nhien lieu (m3tc/kg)
V1_spc=(V_spc*B)/3600 % the tich khoi cho toan bo luong nhien lieu tieu hao (m3tc/s)
% vi khoi luong rieng cua khoi thay doi theo nhiet do, ma trong suot qua trinh trao doi nhiet thi nhiet do khoi co thay doi
% vay de don gian bai toan hon ta gia su rang khoi luong rien cua khoi ko thay doi nhieu va bang hang so (tra bang 23 sach
% nhiet dong luc hoc ky thuat va truyen nhiet cua hoang dinh tin va bui hai o(p=760 mm Hg; PCO2 =0,13;PH2O = 0,11;PN2 =0,76)
% ung voi nhiet do trung binh cua khoi vao va ra khoi bo eco Rho_k=0.83 % khoi luong rieng cua khoi (kg/m3tc)
G_k=V1_spc*Rho_k % luu luong khoi luong cua khoi (kg/s)
%cac gia tri dau vao cua cum chung cat nuoc phen
t_nc=str2num(get(handles.t_nc,'string')); %nhiet do nuoc cap (0C)
cp_nc=4.174; % nhiet dung rieng cua nuoc khi biet nhiet do cua nuoc cap (Kj/kg.do)
t_inc=[28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45]% gan bien nhiet do nuoc phen cap (C)
i_inc=t_inc*cp_nc % entanpi cua nuoc cap khi biet nhiet do va ap suat cua nuoc cap (Kj/kg)
inc= polyfit(t_inc,i_inc,15); % tim phuong trinh bac 2 di qua qua t_inc,i_inc
i_nc=polyval(inc,t_nc) % entanpi cua nuoc cap tuong ung voi nhiet do (kj/kg)
t_sp=str2num(get(handles.t_sp,'string')); %nhiet do san pham (0C)
cp_sp=4.174; % nhiet dung rieng cua san pham khi biet nhiet do cua san pham (Kj/kg.do)
t_isp=[35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50]% gan bien nhiet do nuoc san pham (C)
i_isp=t_isp*cp_sp % entanpi cua san pham khi biet nhiet do cua san pham (Kj/kg)
isp= polyfit(t_isp,i_isp,14); % tim phuong trinh bac 2 di qua qua t_isp,i_isp
i_sp=polyval(isp,t_sp) % entanpi cua san pham tuong ung voi nhiet do san pham (kj/kg)
p_pl=str2num(get(handles.p_pl,'string')); %ap suat binh phan ly (bar) p_ipl=[0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2]; %gan bien ap suat binh phan ly (bar)
t_ipl=[81.33 85.94 89.95 93.5 96.71 99.62 102.3 104.8 107.1 109.3 11.4 113.3 115.2 116.9 118.6 120.2]; %gan bien nhiet do theo ap suat binh phan ly (C)
tpl= polyfit(p_ipl,t_ipl,15); % tim phuong trinh bac 2 di qua qua p_ipl,t_ipl
t_pl=polyval(tpl,p_pl) % nhiet do binh phan ly (C)
i_ipl1=[340.5 359.8 376.7 391.6 405.1 417.4 428.8 439.3 449.1 458.4 467.1 475.3 483.2 490.7 497.8 504.7];% gan bien entanpi cua nuoc chua soi khi biet ap suat binh phan ly (Kj/kg)
ipl1= polyfit(p_ipl,i_ipl1,14); % tim phuong trinh bac 2 di qua qua p_ipl,i_ipl1
i_pl1=polyval(ipl1,p_pl) %entanpi cua nuoc chua soi khi biet ap suat binh phan ly (Kj/kg)
i_ipl2=[2646 2653 2660 2666 2671 2675 2680 2684 2687 2690 2694 2696 2699 2702 2704 2707];% gan bien entanpi cua hoi bao hoa khi biet ap suat binh phan ly (Kj/kg)
ipl2= polyfit(p_ipl,i_ipl2,14); % tim phuong trinh bac 2 di qua qua p_ipl,i_ipl2
i_pl2=polyval(ipl2,p_pl) %entanpi cua hoi bao hoa khi biet ap suat binh phan ly (Kj/kg)
i_ngung=i_pl1; %entanpi cua nuoc ngung khi biet ap suat ngung tu (Kj/kg)
HV: NGUYỄN VŨ BẰNG Page 69 i_xa=i_pl1; % entanpi nuoc xa day
p_r=p_pl;
p_ri=[0.5 1 1.43 2];
r_ipl2=[2286.8 2256.8 2230 2292.8]; % gan bien nhiet an hoa hoi ripl2= polyfit(p_ri,r_ipl2,3); % tim phuong trinh bac 2 di qua qua p_ipl,r_pl2
r_pl2=polyval(ripl2,p_r) %nhiet an hoa hoi cua hoi bao hoa khi biet ap suat binh phan ly (Kj/kg)
lambda_ipl2=[2.25e-02 2.372e-02 2.489e-02 2.593e-02];%gan bien he so trao doi nhiet cua hoi (W/m.K)
lambdaipl2=polyfit(p_ri,lambda_ipl2,3);% tim phuong trinh bac 2 di qua qua lambda_pl2=polyval(lambdaipl2,p_r) %he so trao doi nhiet cua h?i (W/m.K) nu_ipl2=[25e-06 20.02e-06 15.07e-06 11.46e-06];%gan bien do nhot dong hoc cua hoi (m2/s)
nuipl2=polyfit(p_ri,nu_ipl2,3);% tim phuong trinh bac 2 di qua qua nu_pl2=polyval(nuipl2,p_r) %do nhoc dong hoc trung binh cua hoi (m2/s)
rho_ipl2=[0.45 0.598 0.886 1.121];%gan bien khoi luong rieng trung binh cua h?i (kg/m3)
rhoipl2=polyfit(p_ri,rho_ipl2,3);% tim phuong trinh bac 2 di qua qua rho_pl2=polyval(rhoipl2,p_r) %khoi luong rieng cua hoi (kg/m3)
p_eco=str2num(get(handles.p_eco,'string')); % ap suat vao bo eco (bar) t_n41=t_pl; % nhiet do nuoc vao bo eco (C)
cp_tbeco=4.23 % nhiet dung rieng trung binh cua nuoc trong dan eco (kj/kg.C)
i_n41=t_n41*cp_tbeco %entanpi cua nuoc vao bo eco khi biet nhiet do va ap suat cua nuoc (Kj/kg)
t_n42=135; %nhiet do nuoc ra khoi bo eco (C)
t_in42=[t_n41 t_n41+10 t_n41+20 t_n41+30 t_n41+40 t_n41+50 ]% gan bien ap suat nuoc vao bo eco (C)
i_in421=t_in42*cp_tbeco;%entanpi cua nuoc vao bo eco khi nhiet do va ap suat cua nuoc (Kj/kg)
in42=polyfit(t_in42,i_in421,4);% tim phuong trinh bac 2 di qua qua t_in42,i_in421,4
i_n42=polyval(in42,t_n42) %entanpi cua nuoc ra khoi bo eco khi nhiet do va ap suat cua nuoc (Kj/kg)
nuy_tbnt=0.9;
t1_k=str2num(get(handles.t1_k,'string')); % nhiet do khoi vao (0C)
t1_ik=[200 220 240 260 280 300]; %gan bien nhiet do cua khoi vao dan eco (C)
Ct_ikk=[216.24 251.97 287.7 323.43 359.16 394.89]; % tich so nhiet dung rieng va nhiet do cua kk trong khoi (kJ/m3tc)
Ctik= polyfit(t1_ik,Ct_ikk,4); % tim phuong trinh bac 2 di qua qua t1_ik,Ct_ikk
Ct_kk=polyval(Ctik,t1_k) % tich so nhiet dung rieng va nhiet do cua kk (kj/m3tc)
Ct_ico2=[357.4 397.68 437.96 478.246 518.52 558.81]; % tich so nhiet dung rieng va nhiet do cua CO2 trong khoi (kJ/m3tc)
Ctico2= polyfit(t1_ik,Ct_ico2,4); % tim phuong trinh bac 2 di qua qua t1_ik,Ct_ico2
Ct_co2=polyval(Ctico2,t1_k) % tich so nhiet dung rieng va nhiet do cua CO2 (kj/m3tc)
Ct_in2=[259.02 267.618 276.216 284.814 293.412 302.01]; % tich so nhiet dung rieng va nhiet do cua N2 trong khoi (kJ/m3tc)
Ctin2= polyfit(t1_ik,Ct_in2,4); % tim phuong trinh bac 2 di qua qua t1_ik,Ct_in2
Ct_n2=polyval(Ctin2,t1_k) % tich so nhiet dung rieng va nhiet do cua N2 (kj/m3tc)
Ct_ih2o=[304.46 336.112 367.764 399.416 431.068 462.72]; % tich so nhiet dung rieng va nhiet do cua H2O trong khoi (kJ/m3tc)
Ctih2o= polyfit(t1_ik,Ct_ih2o,4); % tim phuong trinh bac 2 di qua qua t1_ik,Ct_ih2o
HV: NGUYỄN VŨ BẰNG Page 70 Ct_h2o=polyval(Ctih2o,t1_k) % tich so nhiet dung rieng va nhiet do cua H2O (kj/m3tc)
i1_0k=V0_RO2*Ct_co2+(V0_nito*Ct_n2)+(V0_H2O*Ct_h2o)% entanpi cua san pham chay li thuyet (kj/kg)
i1_k=i1_0k+((alpha_kk-1)*(V0_kk*Ct_kk)) % entanpi thuc te cua khoi thai (kj/kg)
t2_k=str2num(get(handles.t2_k,'string')); % nhiet do khoi ra khoi dan eco(0C)
t2_ik=[100 120 140 160]; %gan bien nhiet do cua khoi vao dan eco (C) Ct_i2kk=[129.95 147.208 164.466 181.724]; % tich so nhiet dung rieng va nhiet do cua kk trong khoi (kJ/m3tc)
Cti2k= polyfit(t2_ik,Ct_i2kk,3); % tim phuong trinh bac 2 di qua qua t2_ik,Ct_i2kk
Ct_2kk=polyval(Cti2k,t2_k) % tich so nhiet dung rieng va nhiet do cua kk (kj/m3tc)
Ct_i2co2=[170.03 207.504 244.978 282.452]; % tich so nhiet dung rieng va nhiet do cua CO2 trong khoi (kJ/m3tc)
Cti2co2= polyfit(t2_ik,Ct_i2co2,3); % tim phuong trinh bac 2 di qua qua t2_ik,Ct_i2co2
Ct_2co2=polyval(Cti2co2,t2_k) % tich so nhiet dung rieng va nhiet do cua CO2 (kj/m3tc)
Ct_i2n2=[129.58 155.468 181.356 297.244]; % tich so nhiet dung rieng va nhiet do cua N2 trong khoi (kJ/m3tc)
Cti2n2= polyfit(t2_ik,Ct_i2n2,3); % tim phuong trinh bac 2 di qua qua t2_ik,Ct_i2n2
Ct_2n2=polyval(Cti2n2,t2_k) % tich so nhiet dung rieng va nhiet do cua N2 (kj/m3tc)
Ct_i2h2o=[151.02 181.708 212.396 243.084]; % tich so nhiet dung rieng va nhiet do cua H2O trong khoi (kJ/m3tc)
Cti2h2o= polyfit(t2_ik,Ct_i2h2o,3); % tim phuong trinh bac 2 di qua qua t2_ik,Ct_i2h2o
Ct_2h2o=polyval(Cti2h2o,t2_k) % tich so nhiet dung rieng va nhiet do cua H2O (kj/m3tc)
i2_0k=V0_RO2*Ct_2co2+(V0_nito*Ct_2n2)+(V0_H2O*Ct_2h2o)% entanpi cua san pham chay li thuyet (kj/kg)
i2_k=i2_0k+((alpha_kk-1)*(V0_kk*Ct_2kk)) % entanpi thuc te cua khoi thai ra khoi bo eco (kj/kg)
cp_k=1.08975; % nhiet dung rieng cua khoi (Kj/kg.do)
% vi tri thiet bi lam lanh (so 1)
i_1nc=0.6*(i_ngung-i_sp)+i_nc % phuong trinh (2) entanpi cua nuoc cap khi di qua thiet bi lam lanh
t_1nc=i_1nc/4.174 % phuong trinh (3) nhiet do cua nuoc cap khi di qua thiet bi lam lanh
% vi tri thiet bi ngung tu (so 2)
% chon do gia tang nhiet do khong toi muc cua nuoc sau khi qua thiet bi ngung tu la 10 (C)
t_21nc=t_pl-10 % nhiet do nuoc cap sau khi ra khoi thiet bi ngung tu
% ung voi nhiet do dau ra cua nuoc cap va ap suat ta co the xac dinh duoc i_2nc=0.6*0.065*(i_pl2-i_ngung)+i_1nc % phuong trinh entanpi cua nuoc cap khi di qua thiet bi ngung tu
t_2nc=i_2nc/4.174 % phuong trinh (3) nhiet do cua nuoc cap khi di qua thiet bi ngung tu
% vi tri thiet bi Eco (so 4)
G_n4=(G_k*cp_k*(t1_k-t2_k))/(i_n42-i_n41)% luu luong nuoc vao dan eco
% vi tri thiet bi phan ly (so 3)
G_nc=G_n4*(i_n42-i_n41)/((0.6*i_pl2)+(0.4*i_pl1)-i_2nc) % luu luong nuoc vao he thong
% ma ta lai co
G_h=0.6*G_nc % luu luong hoi di ra khoi binh phan ly G_xa=0.4*G_nc % luu luong nuoc xa ra khoi binh phan ly
% tinh toan chon thiet bi cho he thong
HV: NGUYỄN VŨ BẰNG Page 71
% tinh toan thiet ke bo ECO
d_2=str2num(get(handles.d_2,'string')); %(m)%Chon duong kinh ngoai cua ong bo eco (m)
d_1=str2num(get(handles.d_1,'string')); %(m)%Chon duong kinh trong cua ong bo eco (m)
omega_n=str2num(get(handles.omega_n,'string')); %(m/s) chon van toc nuoc trong ong
omega_k=str2num(get(handles.omega_k,'string')); % (m/s)chon van toc khoi di qua chum ong
lambda_ong=22; %(w/m.K) he so dan nhiet cua ong thep
% nhiet do trung binh cua nuoc khi vao va ra khoi ong ECO t_ntb=(t_n41+t_n42)/2
% tra bang thong so vat ly cua nuoc tren duong bao hoa theo nhiet do ta co t_intb=[60 80 100 120 140 160 180 200]% gan bien nhiet do trung binh nuoc trong bo eco (C)
cp_intb=[4.179 4.195 4.22 4.25 4.287 4.346 4.417 4.5];%gan bien nhiet dung rieng trung binh nuoc (Kj/kg.C)
cpintb=polyfit(t_intb,cp_intb,6);% tim phuong trinh bac 2 di qua qua t_intb,cp_intb
cp_ntb=polyval(cpintb,t_ntb) %entanpi cua nuoc ra khoi bo eco khi nhiet do va ap suat cua nuoc (Kj/kg.C)
lambda_intb=[65.9e-02 67.4e-02 68.3e-02 68.6e-02 68.5e-02 68.3e-02 67.4e-02 66.3e-02];%gan bien he so trao doi nhiet trung binh nuoc (W/m.K)
lambdaintb=polyfit(t_intb,lambda_intb,6);% tim phuong trinh bac 2 di qua qua t_intb,lambda_intb
lambda_ntb=polyval(lambdaintb,t_ntb) %he so trao doi nhiet cua nuoc (W/kg.K)
nu_intb=[0.478e-06 0.365e-06 0.295e-06 0.252e-06 0.217e-06 0.191e-06 0.173e-06 0.158e-06];%gan bien do nhot dong hoc trung binh nuoc (m2/s) nuintb=polyfit(t_intb,nu_intb,6);% tim phuong trinh bac 2 di qua qua t_intb,nu_intb
nu_ntb=polyval(nuintb,t_ntb) %do nhoc dong hoc trung binh cua nuoc (m2/s) pr_intb=[2.98 2.21 1.75 1.47 1.26 1.1 1 0.93];%gan bien he so Pr trung binh nuoc
printb=polyfit(t_intb,pr_intb,6);% tim phuong trinh bac 2 di qua qua t_intb,pr_intb
pr_ntb=polyval(printb,t_ntb) %he so Pr trung binh cua nuoc
rho_intb=[983.2 971.8 958.4 943.1 926.1 907 886.9 863];%gan bien khoi luong rieng trung binh nuoc (kg/m3)
rhointb=polyfit(t_intb,rho_intb,6);% tim phuong trinh bac 2 di qua qua t_intb,rho_intb
rho_ntb=polyval(rhointb,t_ntb) %khoi luong rieng cua nuoc ra khoi bo eco (kg/m3)
% he so toa nhiet ve phia nuoc
Re=omega_n*d_1/nu_ntb % he so reynold cua nuoc chay trong ong
% de tinh Nu_f ta can biet nhiet do be mat trong cua vach ong tw tuy nhien
% do he so toa nhiet ve phia nuoc lon hon rat nhieu so voi he so toa nhiet
% cua khoi nen nhiet tro toa nhiet ve phia nuoc ko dang ke, nen tw khong
% khac t_ntb vi vay(P_r/P_w)=1, va l_ong/d>50, epsilon_l=1, epsilon_r=1 if Re>10^3
Nu_f=0.021*(Re^0.8)*(pr_ntb^0.43) end
% he so toa nhiet cua nuoc trong ong
alpha_n=Nu_f*lambda_ntb/d_1 % he so toa nhiet cua nuoc (W/m2.K)
% nhiet do trung binh cua khoi vao va ra khoi bo ECO t_ktb=(t1_k+t2_k)/2
% tra bang thong so vat ly cua nuoc tren duong bao hoa theo nhiet do ta co t_iktb=[100 125 150 175 200 225 250 275 300]% gan bien nhiet do trung binh nuoc trong bo eco (C)