Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 128 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
128
Dung lượng
4,56 MB
Nội dung
http://www.ebook.edu.vn ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƢỜNG ĐẠI HỌC BÁCH KHOA KHOA CƠ KHÍ BỘ MƠN CƠNG NGHỆ NHIỆT LẠNH -oOo - LUẬN VĂN TỐT NGHIỆP THIẾT KẾ HỆ THỐNG THU HỒI NHIỆT THẢI ĐỂ CUNG CẤP NƢỚC LẠNH BẰNG MÁY LẠNH HẤP THỤ TẠI CÔNG TY TAE KWANG VINA GVHD : TS NGUYỄN VĂN TUYÊN SVTH : TRƢƠNG QUỐC BẢO MSSV : 20500166 TP HỒ CHÍ MINH, tháng 01/2010 i Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi LỜI CẢM ƠN Lời xin chân thành cảm ơn TS Nguyễn Văn Tuyên tận tình hƣớng dẫn, truyền đạt kiến thức kinh nghiệm có nhiều góp ý quan trọng giúp tơi hồn thành luận văn Xin cảm ơn thầy cô môn Công Nghệ Nhiệt Lạnh, khoa Cơ Khí nhƣ thầy trƣờng ĐH Bách Khoa TP HCM tận tình giảng dạy suốt thời gian học tập trƣờng Xin cảm ơn gia đình, ngƣời thân bạn bè quan tâm giúp đỡ năm học qua Sinh viên Trƣơng Quốc Bảo ii http://www.ebook.edu.vn TÓM TẮT LUẬN VĂN Nội dung luận văn trình bày vấn đề thu hồi nhiệt thải: điều kiện, sở cần thiết để đƣa phƣơng án thu hồi nhiệt thải hiệu quả, hợp lý tối ƣu Thơng qua phân tích, đánh giá trạng nguồn nhiệt thải doanh nghiệp cụ thể (công ty Tae Kwang Vina), từ đƣa phƣơng án tối ƣu thu hồi nhiệt thải khói phân ly để cung cấp nƣớc lạnh máy lạnh hấp thụ Luận văn trình bày cách tính tốn thiết kế thiết bị thu hồi nhiệt thải hệ thống nhƣ chọn lựa thiết bị khác nhƣ máy lạnh hấp thụ, bơm… Bên cạnh đó, luận văn trình bày cách tính tốn chu trình máy lạnh hấp thụ Single Effect chƣơng trình tính tốn ngơn ngữ C# Ngồi ra, luận văn nghiên cứu tính tốn hiệu biện pháp thu hồi nhiệt phân ly ejector nhằm mục tiêu tiết kiệm chi phí nhiên liệu Qua việc đánh giá tính kinh tế hệ thống, luận văn đƣa số nhận định điều kiện thực tế Việt Nam ứng dụng phƣơng án thu hồi nhiệt thải để cấp nhiệt cho máy lạnh hấp thụ iii http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi MỤC LỤC ĐỀ MỤC TRANG TRANG BÌA i LỜI CẢM ƠN ii TÓM TẮT LUẬN VĂN iii MỤC LỤC iv DANH MỤC CÁC HÌNH VẼ vii DANH MỤC CÁC BẢNG BIỂU ix MỞ ĐẦU CHƢƠNG 1: TỔNG QUAN VỀ THU HỒI NHIỆT THẢI 1.1 Nhiệt thải 1.2 Cơ sở lý thuyết thu hồi nhiệt thải 1.2.1 Điều kiện cần để thu hồi nhiệt thải 1.2.2 Đặc điểm nguồn nhiệt thải 1.2.3 Khó khăn nhƣợc điểm hệ thống thu hồi nhiệt thải 1.3 Ví dụ hệ thống tận dụng nhiệt thải CHƢƠNG 2: TÌNH HÌNH SỬ DỤNG NĂNG LƢỢNG TẠI DOANH NGHIỆP 2.1 Giới thiệu doanh nghiệp 2.2 Quy trình cơng nghệ sản xuất 2.3 Đánh giá việc sử dụng nguồn lƣợng 2.3.1 Tổng quát 2.3.2 Về nguồn nhiệt thải 14 2.4 Đặt vấn đề 14 CHƢƠNG 3: LỰA CHỌN PHƢƠNG ÁN THU HỒI NHIỆT THẢI ĐỂ CUNG CẤP NƢỚC LẠNH 15 3.1 Tính cơng suất nhiệt thu hồi 15 3.1.1 Tính cơng suất nhiệt thu hồi phân ly 15 3.1.2 Tính cơng suất nhiệt thu hồi khói thải 16 3.2 Phân tích, lựa chọn phƣơng án 22 3.2.1 Phƣơng án 22 3.2.2 Phƣơng án 23 3.2.3 Phƣơng án 24 3.2.4 Kết luận, lựa chọn phƣơng án thiết kế 25 CHƢƠNG 4: TÍNH TỐN LỰA CHỌN MÁY LẠNH HẤP THỤ 26 4.1 Tổng quan máy lạnh hấp thụ 26 4.1.1 Sự khác biệt máy lạnh hấp thụ máy lạnh có máy nén 26 iv http://www.ebook.edu.vn 4.1.2 Nguyên lý làm việc máy lạnh hấp thụ H2O – LiBr loại Single Effect 27 4.1.3 Lựa chọn máy lạnh hấp thụ 28 4.2 Tính tốn chu trình máy lạnh hấp thụ 28 4.2.1 Các cơng thức dùng để tính tốn nhiệt động 28 4.2.2 Các phƣơng trình cân nhiệt trình tự tính tốn 36 CHƢƠNG 5: TÍNH TỐN THIẾT KẾ THIẾT BỊ THU HỒI NHIỆT THẢI 44 5.1 Tính toán sơ đồ nhiệt 44 5.1.1 Tính lƣợng trích bổ sung 44 5.1.2 Tính lƣợng nƣớc qua thiết bị thu hồi nhiệt 47 5.2 Thiết bị thu hồi nhiệt phân ly 48 5.3 Thiết bị thu hồi nhiệt khói thải 53 5.3.1 Thiết bị thu hồi nhiệt khói thải lị số 53 5.3.2 Thiết bị thu hồi nhiệt khói thải lị số 60 5.4 Thiết bị trao đổi nhiệt bổ sung bồn chứa nƣớc nóng: 65 5.5 Tính tốn trở lực 72 5.5.1 Trở lực đƣờng nƣớc nóng hệ thống 72 5.5.2 Trở lực đƣờng nƣớc nóng cung cấp cho MLHT 76 5.5.3 Trở lực đƣờng nƣớc lạnh cung cấp cho công nghệ 78 5.5.4 Trở lực đƣờng nƣớc giải nhiệt 79 5.6 Chọn bơm 81 CHƢƠNG 6: ĐÁNH GIÁ HIỆU QUẢ KINH TẾ 82 6.1 Chi phí đầu tƣ vận hành 82 6.1.1 Chi phí đầu tƣ 82 6.1.2 Chi phí vận hành 82 6.1.3 Chi phí tiết kiệm đƣợc hệ thống: 83 6.2 Tính toán lại phƣơng án 84 6.2.1 Chọn máy lạnh hấp thụ 84 6.2.2 Tính tốn thiết bị thu hồi nhiệt thải 85 6.2.3 Chọn thiết bị khác 87 6.2.4 Tính tốn lại chi phí 88 6.3 Nhận xét 89 PHỤ CHƢƠNG: HỆ THỐNG EJECTOR HƠI ĐỂ THU HỒI NHIỆT LƢỢNG HƠI PHÂN LY 91 P.1 Tổng quan ejector 91 P.2 Tính tốn ejector 92 P.2.1 Cơ sở lý thuyết 92 P.2.2 Tính toán ejector 95 v http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi P.3 Nhận xét 98 KẾT LUẬN 101 TÀI LIỆU THAM KHẢO 102 PHỤ LỤC 103 PHỤ LỤC 107 PHỤ LỤC 114 PHỤ LỤC 116 PHỤ LỤC 117 vi http://www.ebook.edu.vn DANH MỤC CÁC HÌNH VẼ Hình 2.1: Quy trình cơng nghệ sản xuất giày công ty Tae Kwang Vina Hình 2.2: Đồ thị phân bố tỉ lệ tiêu thụ điện thiết bị doanh nghiệp Hình 2.3: Sơ đồ hệ thống cung cấp nƣớc lạnh cho cơng nghệ Hình 2.4: Lò đốt dầu FO với hệ thống béc đốt quay 10 Hình 2.5: Nguyên lý hoạt động béc đốt quay 11 Hình 2.6: Sơ đồ hệ thống cung cấp nhiệt doanh nghiệp 11 Hình 2.7: Bố trí đƣờng ống khói nhà lị 12 Hình 2.8: Bình góp 12 Hình 2.9: Hơi phân ly khơng đƣợc thu hồi bồn nƣớc cấp 14 Hình 3.1: Đồ thị xác định nhiệt độ đọng sƣơng khói 17 Hình 3.2: Phƣơng án 22 Hình 3.3: Phƣơng án 23 Hình 3.4: Phƣơng án 24 Hình 4.1: Sơ đồ nguyên lý MLHT H2O – LiBr loại Single Effect 27 Hình 4.2: Sơ đồ khối MLHT H2O – LiBr loại Single Effect 36 Hình 4.3: Bình phát sinh 38 Hình 4.4: Bình ngƣng 39 Hình 4.5: Bình bốc 39 Hình 4.6: Bình hấp thụ 40 Hình 4.7: Bộ trao đổi nhiệt 40 Hình 4.8: Giao diện chƣơng trình 42 Hình 4.9: Giao diện “Thơng số điểm đặc trƣng” chu trình 42 Hình 4.10: Đồ thị liên hệ áp suất - nhiệt độ nồng độ dung dịch H2O - LiBr 43 Hình 5.1: Sơ đồ nhiệt phƣơng án 44 Hình 5.2: Thiết bị thu hồi nhiệt phân ly 48 Hình 5.3: Ngƣng tụ thiết bị thu hồi nhiệt phân ly 49 Hình 5.4: Mặt sàng thiết bị thu hồi nhiệt phân ly 49 Hình 5.5: Đồ thị trao đổi nhiệt thiết bị thu hồi nhiệt phân ly 51 vii http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi Hình 5.6: Thiết bị thu hồi nhiệt khói thải 53 Hình 5.7: Các kích thƣớc ống cánh 54 Hình 5.8: Đồ thị tra hiệu suất cánh 56 Hình 5.9: Đồ thị trao đổi nhiệt thiết bị thu hồi nhiệt khói thải lị số 58 Hình 5.10: Đồ thị trao đổi nhiệt thiết bị thu hồi nhiệt khói thải lị số 63 Hình 5.11: Thiết bị trao đổi nhiệt bồn nƣớc nóng 65 Hình 5.12: Bồn nƣớc nóng nhiệt độ nƣớc hoạt động 66 Hình 5.13: Q trình hịa trộn bồn nƣớc nóng 68 Hình 5.14: Đồ thị trao đổi nhiệt thiết bị trao đổi nhiệt bổ sung 70 Hình 5.15: Sơ đồ đƣờng nƣớc nóng qua thiết bị 72 Hình 5.16: Trở lực cục dịng chảy từ tiết diện nhỏ sang tiết diện lớn ngƣợc lại 75 Hình P.1: Ejector 91 Hình P.2: Quá trình giãn nở ejector 92 Hình P.3: Sự thay đổi áp suất vận tốc ejector 93 Hình P.4: Đồ thị i - s 97 Hình P.5: Sơ đồ hệ thống nhiệt Công ty Nikkico 99 viii http://www.ebook.edu.vn DANH MỤC CÁC BẢNG BIỂU Bảng 2.1: Thơng số kỹ thuật lị đốt dầu FO 13 Bảng 3.1: Thành phần dầu FO (%) 16 Bảng 4.1 28 Bảng 5.1: Kết tính lƣợng trích bổ sung 46 Bảng 5.2: Thông số chọn bơm 81 Bảng 5.3: Thông số bơm hệ thống 81 Bảng 6.1: Thông số chọn bơm 88 Bảng 6.2: Thông số bơm hệ thống 88 ix http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi MỞ ĐẦU Cùng với phát triển khoa học cơng nghệ nhu cầu sử dụng lƣợng ngày tăng Hiện nay, ngƣời phải dựa vào nguồn nhiên liệu hóa thạch (than, dầu, khí đốt…) để cung cấp lƣợng cho hoạt động sản xuất, sinh hoạt… Do mức độ khai thác ngày gia tăng nên tƣơng lai sản lƣợng nguồn nhiên liệu giảm dần cạn kiệt Để đối phó với tình trạng này, ngƣời khơng cịn cách khác phải có biện pháp sử dụng tiết kiệm nguồn lƣợng có tìm nguồn lƣợng để thay Bên cạnh vấn đề sản lƣợng, nguồn nhiên liệu hóa thạch cịn gây tác động nguy hại đến mơi trƣờng – vấn đề trở thành mối quan tâm chung tồn cầu Hiện tƣợng nóng dần lên toàn cầu gây tác động đến sinh thái, khí hậu gây hậu lớn mà nguyên nhân phần lớn khí nhà kính, số phải kể đến khí CO2 – sản phẩm q trình đốt cháy nguồn nhiên liệu hóa thạch Một vấn đề đáng quan tâm biến động giá nhiên liệu thị trƣờng Điều ảnh hƣởng không nhỏ đến sản xuất giá thành sản phẩm Tại Việt Nam, việc ứng dụng nguồn lƣợng cịn gặp khó khăn Do đó, nâng cao hiệu sử dụng tiết kiệm lƣợng mối quan tâm hàng đầu Nó khơng lợi ích doanh nghiệp, cá nhân mà cịn lợi ích quốc gia hay rộng cho toàn cầu http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi //Phụ tải nhiệt BPS double Qh = mr * qh; //Năng suất giải nhiệt đơn vị BN double qk = i2 - i2_1; //Năng suất giải nhiệt BN double Qk = mr * qk; //Năng suất giải nhiệt đơn vị BHT double qa = i3_2 + (a - 1) * i5 - a * i4; //Năng suất giải nhiệt BHT double Qa = mr * qa; //Hệ số COP double COP = Qo / Qh; txt_Qa.Text = Qa.ToString(); txt_Qh.Text = Qh.ToString(); txt_Qk.Text = Qk.ToString(); txt_COP.Text = COP.ToString(); } private void button1_Click(object sender, EventArgs e) { Thong_so_cac_diem_dac_trung show = new Thong_so_cac_diem_dac_trung(); show.Show(); } private void Thong_so_cac_diem_dac_trung_Load(object sender, EventArgs e) { txt_c1.Text = Form1.Cw.ToString(); txt_c4.Text = Form1.Cw.ToString(); txt_c7.Text = Form1.Cw.ToString(); txt_c5.Text = Form1.Cs.ToString(); txt_c6.Text = Form1.Cs.ToString(); txt_p1.Text = Form1.pk.ToString(); txt_p2.Text = Form1.pk.ToString(); txt_p2_1.Text = Form1.pk.ToString(); txt_p6.Text = Form1.pk.ToString(); 105 http://www.ebook.edu.vn txt_p7.Text = Form1.pk.ToString(); txt_p3.Text = Form1.po.ToString(); txt_p3_2.Text = Form1.po.ToString(); txt_p4.Text = Form1.po.ToString(); txt_p5.Text = Form1.po.ToString(); txt_t1.Text = Form1.t1.ToString(); txt_t2.Text = Form1.t2.ToString(); txt_t2_1.Text = Form1.tk.ToString(); txt_t3.Text = Form1.to.ToString(); txt_t3_2.Text = Form1.to.ToString(); txt_t4.Text = Form1.t4.ToString(); txt_t5.Text = Form1.t5.ToString(); txt_t6.Text = Form1.t6.ToString(); txt_i1.Text = Form1.i1.ToString(); txt_i2.Text = Form1.i2.ToString(); txt_i2_1.Text = Form1.i2_1.ToString(); txt_i3.Text = Form1.i3.ToString(); txt_i3_2.Text = Form1.i3_2.ToString(); txt_i4.Text = Form1.i4.ToString(); txt_i5.Text = Form1.i5.ToString(); txt_i6.Text = Form1.i6.ToString(); txt_i7.Text = Form1.i7.ToString(); } 106 http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi PHỤ LỤC CHƢƠNG TRÌNH CON TÍNH TỐN CÁC THƠNG SỐ NHIỆT ĐỘNG //Tính áp suất bão hòa tác nhân lạnh public static double tinh_ap_suat_TNL(double nhiet_do)//Đơn vị bar / độ C { double kq; double log_p; double ap_suat;//Đơn vị bar if (nhiet_do < 5) { log_p = 10.5380997 - 2663.91 / (nhiet_do + 273.15); ap_suat = Math.Pow(10, log_p) * 0.001; } else { log_p = 28.59051 - 8.2 * Math.Log10(nhiet_do + 273.15) + 0.0024804 * (nhiet_do + 273.15) - 3142.31 / (nhiet_do + 273.15); ap_suat = Math.Pow(10, log_p); } kq = ap_suat; return kq; } //Tính enthalpy dung dịch public static double tinh_enthalpy_dd(double nong_do, double nhiet_do) { double enthalpy=0; double[,] heso = new double[6, 3]{{1.134125,4.124891,5.743693*1E-4},{4.80045/10,-7.643903/100,5.870921/100000}, {-2.161438/1000,2.589577/1000,-7.375319/1000000},{2.336235/10000,9.500522/100000,3.277592/10000000},{-1.188679/100000,1.708026/1000000,6.062304/1000000000}, {2.291532*1E-7,-1.102363*1E-8,3.901897*1E-11}}; for (int j = 0; j < 3; j++) { for (int i = 0; i < 6; i++) { double number = heso[i, j]; enthalpy += number * Math.Pow(nong_do, i) * Math.Pow(nhiet_do, j); } } 107 http://www.ebook.edu.vn return enthalpy; } //Tính nhiệt độ bão hịa tác nhân lạnh public static double tinh_nhiet_do_bh_TNL(double ap_suat)//Đơn vị oC,bar { double nhiet_do = 0; ap_suat = ap_suat / 10; if (ap_suat < 12.33) nhiet_do = 42.6776 - 3892.7 / (Math.Log(ap_suat, Math.E) - 9.48654); else nhiet_do = -387.592 - 12587.5 / (Math.Log(ap_suat, Math.E) - 15.2578); nhiet_do = nhiet_do - 273.15; return nhiet_do; } //Tính enthalpy trạng thái bão hịa TNL public static void tinh_enthalpy_TNL(double nhiet_do,out double enthalpy_long,out double enthalpy_hoi) { enthalpy_hoi = enthalpy_long = 0; nhiet_do = nhiet_do + 273.15; double number_1 = 0, number_2 = 0; double[] heso_a = new double[8] { 8.839230108 * 1E-1, -2.67172935, 6.22640035, 13.1789573, -1.91322436, 68.7937653, -1.24819906 * 1E2, 72.1435404 }; double[] heso_b = new double[8] { 4.57874342 * 1E-1, 5.08441288, -1.48513244, 4.81351884, 2.69411782, -7.39064542, 10.4961689, -5.46840036 }; double heso_TR = (647.3 - nhiet_do) / 647.3; for (int i = 1; i < 8; i++) { number_1 += heso_a[i] * Math.Pow(heso_TR, i); } enthalpy_long = (number_1+heso_a[0]) * 2099.3; for (int j = 3; j < 8; j++) { number_2 += heso_b[j] * Math.Pow(heso_TR, j - 2); } enthalpy_hoi = (number_2+1 + heso_b[0] * Math.Pow(heso_TR, 1.0 / 3) + heso_b[1] * Math.Pow(heso_TR, 5.0 / 6) + heso_b[2] * Math.Pow(heso_TR, 7.0 / 8)) * 2099.3; } 108 http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi //Tính nhiệt độ sơi dung dịch public static double tinh_nhiet_do_soi_dd(double nong_do, double ap_suat) { double A,B,D, E, F, N; double nhiet_do = 0; ap_suat = ap_suat * 0.0193367136; A = -2.00755 + 0.16976 * nong_do - 3.133362 * 1E-3 * nong_do * nong_do + 1.97668 * 1E-5 * nong_do * nong_do * nong_do; B = 321.128 - 19.322 * nong_do + 0.374382 * nong_do * nong_do - 2.0637 * 1E-3 * nong_do * nong_do * nong_do; D = -2886.373; E = -337269.46; F = 6.21147; N = Math.Log10(ap_suat); double ngoac = (-2 * E / (D + Math.Pow((Math.Pow(D, 2) - * (F - N) * E), 0.5))) 459.72; nhiet_do = A * ngoac + B; nhiet_do = * (nhiet_do - 32) / 9; return nhiet_do; } //Tính nồng độ dung dịch public static double tinh_nong_do_dd(double nhiet_do_soi_dd, double nhiet_do_bh_TNL) { double nong_do = 0; double[] heso = new double[8] { 0.5362, 2.103 * 1E-4, -0.1335, 7.7844 * 1E-4, 4.7942 * 1E-3, -7.4752 * 1E-5, -4.5258 * 1E-5, 6.1135 * 1E-7 }; double delta = nhiet_do_soi_dd - nhiet_do_bh_TNL; nong_do = 38.3893 + heso[0] * delta + heso[1] * delta * delta + heso[2] * nhiet_do_bh_TNL + heso[3] * nhiet_do_bh_TNL * nhiet_do_bh_TNL + heso[4] * delta * nhiet_do_bh_TNL + heso[5] * delta * delta * nhiet_do_bh_TNL + heso[6] * delta * nhiet_do_bh_TNL * nhiet_do_bh_TNL + heso[7] * delta * delta * nhiet_do_bh_TNL * nhiet_do_bh_TNL; return nong_do; } //Tính enthalpy nhiệt public static double tinh_enthalpy_hqn(double ap_suat, double nhiet_do) { double enthalpy; double[] n0 = new double[9]; 109 http://www.ebook.edu.vn n0[0] = -0.969276865002 * 1E1; n0[1] = 0.10086655968018 * 1E2; n0[2] = -0.56087911283020 * 1E-2; n0[3] = 0.71452738081455 * 1E-1; n0[4] = -0.40710498223928; n0[5] = 0.14240819171444 * 1E1; n0[6] = -0.43839511319450 * 1E1; n0[7] = -0.28408632460772; n0[8] = 0.21268463753307 * 1E-1; double[] J0 = new double[9]; J0[0] = 0; J0[1] = 1; J0[2] = -5; J0[3] = -4; J0[4] = -3; J0[5] = -2; J0[6] = -1; J0[7] = 2; J0[8] = 3; double[] n = new double[43]; n[0] = -0.17731742473213 * 1E-2; n[1] = -0.17834862292358 * 1E-2; n[2] = -0.45996013696365 * 1E-1; n[3] = -0.57581259083432 * 1E-1; n[4] = -0.50325278727930 * 1E-1; n[5] = -0.33032641670203 * 1E-4; n[6] = -0.18948987516315 * 1E-3; n[7] = -0.39392777243355 * 1E-2; n[8] = -0.43797295650573 * 1E-1; n[9] = -0.26674547914087 * 1E-4; n[10] = 0.20481737692309 * 1E-7; n[11] = 0.43870667284435 * 1E-6; n[12] = -0.32277677238570 * 1E-4; n[13] = -0.15033924542148 * 1E-2; n[14] = -0.40668253562649 * 1E-1; n[15] = -0.78847309559367 * 1E-9; n[16] = 0.12790717852285 * 1E-7; n[17] = 0.48225372718507 * 1E-6; n[18] = 0.22922076337661 * 1E-5; n[19] = -0.16714766451061 * 1E-10; n[20] = -0.21171472321355 * 1E-2; n[21] = -0.23895741934104 * 1E2; n[22] = -0.59059564324270 * 1E-17; 110 http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi n[23] = -0.12621808899101 * 1E-5; n[24] = -0.38946842435739 * 1E-1; n[25] = 0.11256211360459 * 1E-10; n[26] = -0.82311340897998 * 1E1; n[27] = 0.19809712802088 * 1E-7; n[28] = 0.10406965210174 * 1E-18; n[29] = -0.10234747095929 * 1E-12; n[30] = -0.10018179379511 * 1E-8; n[31] = -0.80882908646985 * 1E-10; n[32] = 0.10693031879409; n[33] = -0.33662250574171; n[34] = 0.89185845355421 * 1E-24; n[35] = 0.30629316876232 * 1E-12; n[36] = -0.42002467698208 * 1E-5; n[37] = -0.59056029685639 * 1E-25; n[38] = 0.37826947613457 * 1E-5; n[39] = -0.12768608934681 * 1E-14; n[40] = 0.73087610595061 * 1E-28; n[41] = 0.55414715350778 * 1E-16; n[42] = -0.94369707241210 * 1E-6; double[] J = new double[43]; J[0] = 0; J[1] = 1; J[2] = 2; J[3] = 3; J[4] = 6; J[5] = 1; J[6] = 2; J[7] = 4; J[8] = 7; J[9] = 36; J[10] = 0; J[11] = 1; J[12] = 3; J[13] = 6; J[14] = 35; J[15] = 1; J[16] = 2; J[17] = 3; J[18] = 7; J[19] = 3; J[20] = 16; J[21] = 35; 111 http://www.ebook.edu.vn J[22] = 0; J[23] = 11; J[24] = 25; J[25] = 8; J[26] = 36; J[27] = 13; J[28] = 4; J[29] = 10; J[30] = 14; J[31] = 29; J[32] = 50; J[33] = 57; J[34] = 20; J[35] = 35; J[36] = 48; J[37] = 21; J[38] = 53; J[39] = 39; J[40] = 26; J[41] = 40; J[42] = 58; double[] I = new double[43]; I[0] = 1; I[1] = 1; I[2] = 1; I[3] = 1; I[4] = 1; I[5] = 2; I[6] = 2; I[7] = 2; I[8] = 2; I[9] = 2; I[10] = 3; I[11] = 3; I[12] = 3; I[13] = 3; I[14] = 3; I[15] = 4; I[16] = 4; I[17] = 4; I[18] = 5; I[19] = 6; I[20] = 6; 112 http://www.ebook.edu.vn Ket-noi.com Ket-noi.com kho kho tai tai lieu lieu mien mien phi phi I[21] = 6; I[22] = 7; I[23] = 7; I[24] = 7; I[25] = 8; I[26] = 8; I[27] = 9; I[28] = 10; I[29] = 10; I[30] = 10; I[31] = 16; I[32] = 16; I[33] = 18; I[34] = 20; I[35] = 20; I[36] = 20; I[37] = 21; I[38] = 22; I[39] = 23; I[40] = 24; I[41] = 24; I[42] = 24; ap_suat = ap_suat * 133.3224 / 1000000; //Đơn vị mmHg nhiet_do = nhiet_do + 273.15; double p0 = 1; double T0 = 540; double ti_so_ap_suat = ap_suat / p0; double ti_so_nhiet_do = T0 / nhiet_do; double delta0 = 0; for (int i = 0; i