Mô hình hóa quá trình sản xuất hydrogen theo chu trình tuần hoàn oxy hóa khử của các chất mang oxy từ nguồn nguyên liệu mêtan và hơi nước

Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành LỜI CAM ĐOAN  Bản luận văn thạc sĩ Chuyên ngành Kỹ thuật Hóa học với đề tài: “Mơ hình hóa q trình sản xuất hydrogen theo chu trình tuần hồn oxy hóa khử chất mang oxy từ nguồn nguyên liệu khí mêtan nước” hoàn thành hướng dẫn TS Nguyễn Đặng Bình Thành - Bộ mơn Máy Thiết bị cơng nghiệp Hóa chất - Dầu khí, Viện Kỹ thuật Hóa học, Trường Đại học Bách Khoa Hà Nội Tôi xin cam đoan, luận văn không chép nội dung từ luận văn thạc sĩ luận án tiến sĩ khác Hà Nội, ngày 26 tháng 10 năm 2016 Người viết Hoàng Trung Kiên I Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành LỜI CẢM ƠN  Trước tiên tơi xin bày tỏ lịng biết ơn sâu sắc tới TS Nguyễn Đặng Bình Thành, Giảng viên Bộ môn Máy Thiết bị Công nghiệp Hóa Chất - Dầu khí người trực tiếp hướng dẫn chuyên môn, phương pháp nghiên cứu tạo điều kiện giúp đỡ tơi suốt q trình thực đề tài Xin trân trọng cảm Thầy, Cô giáo Viện Kỹ thuật Hóa học - Đại học Bách khoa Hà Nội, môn Máy Thiết bị Công nghiệp Hóa chất - Dầu khí đ tận tình dạy d , gi p đỡ tạo điều kiện thuận lợi cho tơi hồn thành nội dung học tập Cuối xin gửi lời cảm ơn chân thành tới gia đình, bạn bè, đồng nghiệp đ động viên, gi p đỡ thời gian học tập trình thực đề tài Hà Nội, ngày 26 tháng 10 năm 2016 II Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành MỤC LỤC CHƢƠNG TỔNG QUAN VỀ CHU TRÌNH TUẦN HỒN HĨA HỌC .9 1.1 Khái quát chung .9 1.2 Năng lượng tái tạo 11 1.3 Khái niệm chu trình tuần hồn hóa học 13 1.3.1 Tách nước để sản xuất Hyđrơ Ơxy 19 1.3.2 Sản xuất Hydro Peroxit 21 1.3.3 Sản xuất Maleic Anhydride từ butan 21 1.3.4 Một số ví dụ phản ứng tuần hồn hóa học liên quan tới q trình chuyển hóa cacbon 23 CHƢƠNG TỔNG QUAN VỀ CÁC CÔNG NGHỆ TUẦN HỒN HĨA HỌC VÀ CÁC PHƢƠNG PHÁP SẢN XUẤT HYĐRƠ .28 2.1 Các công nghệ tuần hồn hóa học 28 2.2 Các phương pháp sản xuất Hyđrô 42 2.2.1 Giới thiệu chung 42 2.2.2 Các phương pháp sản xuất Hyđrô .43 CHƢƠNG MƠ HÌNH HĨA Q TRÌNH SẢN XUẤT H2 BẰNG CHU TRÌNH TUẦN HỒN ƠXY HĨA KHỬ 49 3.1 Đặc điểm loại vật liệu tuần hoàn .49 3.2 Giới thiệu trình sản xuất H2 chu trình tuần hồn ơxy hóa khử .51 3.3 Mơ hình động học 54 3.4 Mơ hình tầng chuyển động 57 CHƢƠNG KẾT QUẢ VÀ THẢO LUẬN 65 KẾT LUẬN VÀ KIẾN NGHỊ 75 TÀI LIỆU THAM KHẢO 76 PHỤ LỤC – CÁC TÍNH CHẤT NHIỆT ĐỘNG CỦA CÁC CẤU TỬ 89 PHỤ LỤC – CHƢƠNG TRÌNH MATLAB .91 HVTH: Hoàng Trung Kiên Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT STT Kí hiệu Tên gọi Đơn vị A Diện tích mặt cắt ngang thiết bị chuyển hóa nhiên liệu thiết bị chuyển hóa nước m2 Ap Diện tích bề mặt hạt m2 bi Hệ số tỷ lượng cho phản ứng i C Nồng độ pha khí Cgf Nồng độ mêtan thiết bị chuyển hóa nhiên liệu Cgo,fuel Nồng độ mêtan dòng cấp đến thiết bị chuyển hóa nhiên liêu mol/m3 Cgo,steam Nồng độ nước dịng cấp đến thiết bị chuyển hóa nước mol/m3 Cgs Nồng độ nước thiết bị chuyển hóa nước mol/m3 CH2 Nồng độ Hyđrơ thiết bị chuyển hóa nước mol/m3 10 Cp,i Nhiệt dung riêng cấu tử i J/mol.K 11 dp Kích thước hạt mang ơxy 12 Ea Năng lượng hoạt hóa phản ứng J/mol 13 Eaf Năng lượng hoạt hóa phản ứng thiết bị chuyển hóa nhiên liệu J/mol 14 Eas Năng lượng hoạt hóa phản ứng thiết bị chuyển hóa nước J/mol 15 Fgo,fuel Lưu lượng mêtan cấp vào thiết bị chuyển hóa nhiên liệu mol/s 16 Fgo,steam Lưu lượng cấp vào thiết bị chuyển hóa nước mol/s 17 Fso,fuel Lưu lượng Fe2O3 cấp vào thiết bị chuyển hóa nhiên liệu mol/s 18 Fso,steam Lưu lượng Fe3O4 cấp vào thiết bị chuyển hóa nước mol/s 19 F ZrO2 Lưu lượng cấp ZrO2 chất mang ơxy mol/s HVTH: Hồng Trung Kiên mol/m3 mol/m3 m Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành STT Kí hiệu Tên gọi Đơn vị W/m2.K 20 h Hệ số trao đổi nhiệt 21 k Hằng số tốc độ phản ứng 22 k* Hằng số tốc độ phản ứng biểu kiến 23 KE Hằng số cân phản ứng: 3FeO + H2O = Fe3O4 + H2 24 kf Hằng số tốc độ phản ứng biểu kiến thiết bị chuyển hóa nhiên liệu 25 kg Hệ số dẫn nhiệt khí W/m.K 26 ko Hệ số va chạm hữu ích 1/s 27 k *o Hệ số va chạm hữu ích tính toán số tốc độ phản ứng biểu kiến mol-n m3n s-1 28 kof Hệ số va chạm hữu ích tính tốn số tốc độ phản ứng biểu kiến thiết bị chuyển hóa nhiên liệu mol-n m3n s-1 29 kos Hệ số va chạm hữu ích tính tốn số tốc độ phản ứng biểu kiến thiết bị chuyển hóa nước mol-n m3n s-1 30 ks Hằng số tốc độ phản ứng biểu kiến thiết bị chuyển hóa nước mol-n m3n s-1 31 ksurface Hằng số tốc độ phản ứng bề mặt mol1-n m3n-2 s-1 32 L Chiều dày lớp chất rắn phản ứng m 33 L* Chiều cao lớp đệm m 34 mp Khối lượng hạt kg 35 n Bậc phản ứng 36 Np Số lượng hạt đơn vị thể tích 37 R Hằng số khí 38 Sm Diện tích bề mặt trao đổi nhiệt đơn vị thể tích chất mang ôxy 39 t Thời gian phản ứng s 40 T Nhiệt độ tuyệt đối K 41 Tgf Nhiệt độ khí thiết bị chuyển hóa nhiên liệu K HVTH: Hồng Trung Kiên 1/s mol-n m3n s-1 mol-n m3n s-1 m-3 m3barK-1mol-1 m2 m-3 Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành STT Kí hiệu Tên gọi Đơn vị 42 Tgs Nhiệt độ khí thiết bị chuyển hóa nước K 43 Tgo,fuel Nhiệt độ mêtan ban đầu K 44 Tgo,steam Nhiệt độ ban đầu K 45 Tsf Nhiệt độ chất rắn thiết bị chuyển hóa nhiên liệu K 46 Tss Nhiệt độ chất rắn thiết bị chuyển hóa nước K 47 Vp Thể tích hạt m3 48 V Vận tốc dịng khí m/s 49 X, Xs Độ chuyển hóa chất rắn 50 Xgf 51 Xgs Độ chuyển hóa mêtan thiết bị chuyển hóa nhiên liệu Độ chuyển hóa nước thiết bị chuyển hóa nước 52 Xsf Độ chuyển hóa chất rắn thiết bị chuyển hóa nhiên liệu 53 Xss Độ chuyển hóa chất rắn thiết bị chuyển hóa nước 54 z Chiều dài lớp chất rắn thiết bị 55  Hệ số giãn nở khí 56 g Độ nhớt khí 57 particle Khối lượng riêng chất mang ơxy kg/m3 58 bulk Khối lượng riêng đổ đống chất rắn thiết bị chuyển hóa kg/m3 59 g Khối lượng riêng chất khí kg/m3 60 m 61 s Hàm lượng mol chất hoạt động chất mang ôxy Hàm lượng mol chất mang ôxy thiết bị chuyển hóa 62  Thời gian chuyển hóa hồn tồn chất mang ôxy 63 Hrxn1 Nhiệt phản ứng: 4Fe2O3 + CH4 = 8FeO + CO2 + 2H2O HVTH: Hoàng Trung Kiên m kg m-1 s-1 mol/m3 mol/m3 s J/mol Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành STT Kí hiệu Tên gọi 64 Hrxn2 Nhiệt phản ứng: 3FeO + H2O = Fe3O4 + H2 65 AR Air Reactor 66 CFB Circulating Fluidized Bed 67 CLC Chemical Looping Combustion 68 CLP Calcium Looping Process 69 FR Fuel Reactor 70 IGCC Integrated gasification combined cycle 71 SCL Syngas Chemical Looping 72 SR Steam Reactor 73 TRCL Three-Reactors Chemical Looping 74 WGS Water-Gas Shift HVTH: Hoàng Trung Kiên Đơn vị J/mol Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành DANH MỤC CÁC BẢNG Bảng 2.1 Tổng hợp phản ứng hóa học để sản xuất H2 đồng thời thu hồi khí thải CO2 sử dụng phản ứng khử ơxy hóa sắt oxit sắt [8] 47 Bảng 3.1 Các thông số động học 57 Bảng 3.2 Các tính chất vật lý hạt chất mang ơxy 60 Bảng 4.1 Các thông số vận hành .68 Bảng PL Tính chất nhiệt động cấu tử K 298,15 K [9] 89 Bảng PL Các hệ số dùng cơng thức thực nghiệm tính tốn thơng số động học nhiệt độ [9] 90 HVTH: Hoàng Trung Kiên Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành DANH MỤC CÁC HÌNH VẼ Hình 1.1 Dự kiến lượng tồn cầu [18] 10 Hình 1.2 Sơ đồ thu hồi lượng hữu ích hệ thống khí hóa cacbon/ tuabin khí 14 Hình 1.3 Sơ đồ thu hồi lượng hữu ích q trình chuyển dịch khí-hơi nước - khí hóa .16 Hình 1.4 Hệ thống tuần hồn hóa chất kiểu tầng sơi tuần hoàn để sản xuất Maleic Anhydride 23 Hình 2.1 Cơng nghệ sắt – nước sử dụng đệm cố định 29 Hình 2.2 Sơ đồ tuần hồn hóa học tầng sơi sản xuất CO2 từ nhiên liệu rắn 31 Hình 2.3 Sơ đồ cơng nghệ chu trình đơt cháy tuần hồn hóa học CLC 32 Hình 2.4 Sơ đồ chu trình tuần hồn hóa học khí tổng hợp SCL 32 Hình 2.5 Sơ đồ chu trình tuần hốn hóa học OSU’s 34 Hình 2.6 Sơ đồ cơng nghệ HyPr-Ring .35 Hình 2.7 Sơ đồ chu trình tuần hồn canxi CLP [23] 36 Hình 2.8 Sơ đồ công nghệ sử dụng linh hoạt loại nhiên liệu GE [86] 39 Hình 2.9 Sơ đồ khối cơng nghệ đốt khí hóa ALSTOM Hybird [5] 41 Hình 2.10 Sơ đồ nguyên lý chu trình tuần hồn hóa học thiết bị phản ứng sản xuất Hyđrô 42 Hình 3.1 Sơ đồ ngun lý chu trình tuần hồn hóa học thiết bị phản ứng để sản xuất Hyđrô 52 Hình 3.2 (a) Hình dạng thiết bị chuyển hóa nhiên liệu (b) Hình dạng thiết bị chuyển hóa nước 60 Hình 4.1a Sự thay đổi nhiệt độ độ chuyển hóa pha rắn thiết bị chuyển hóa nhiên liệu điều kiện Fe2O3 wt.% = 10, Tsf,in = 1173 K 66 Hình 4.1b Sự thay đổi nhiệt độ độ chuyển hóa pha khí thiết bị chuyển hóa nhiên liệu điều kiện Fe2O3 wt.% = 10, Tsf,in = 1173 K 66 HVTH: Hoàng Trung Kiên Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành Hình 4.2a Sự thay đổi nhiệt độ độ chuyển hóa pha rắn thiết bị chuyển hóa nước điều kiện Fe2O3 wt.% = 10, Tss,in = 1100 K, lưu lượng = 8,53 mol/s Xsf = 0,9 67 Hình 4.2b Sự thay đổi nhiệt độ độ chuyển hóa pha khí thiết bị chuyển hóa nước điều kiện Fe2O3 wt.% = 10, Tss,in = 1100 K, lưu lượng = 8,53 mol/s Xsf = 0,9 67 Hình 4.3 Lượng chất rắn tuần hoàn độ giảm nhiệt theo độ chuyển hóa chất rắn trọng thiết bị chuyển hóa nhiên liệu (Fe2O3 10 wt.%, nhiệt độ chất rắn vào 1173K) 70 Hình 4.4 Thay đổi lượng chất rắn tuần hoàn với nhiệt độ chất rắn vào thiết bị chuyển hóa nhiên liệu (Fe2O3 10 wt.%, độ chuyển hóa chất rắn thiết bị chyển hóa nhiên liệu Xsf = 0,9) 70 Hình 4.5 Thay đổi lượng chất rắn tuần hồn với Fe2O3 wt.% chất mang oxy nhiệt độ chất rắn vào thiết bị chuyển hóa nhiên liệu độ chuyển hóa thiết bị chuyển hóa nhiên liệu 0,9 71 Hình 4.6 Thay đổi lượng chất rắn tuần hồn hàm nhiệt độ chất rắn từ thiết bị chuyển hóa nhiên liệu điều kiện lưu lượng nước = 9,188 mol/s, độ chuyển hóa chất rắn thiết bị chuyển hóa nhiên liệu = 0,9 hàm lượng Fe2O3 = 10 wt.% 20 wt.% .73 Hình 4.7 Lượng chất rắn tuần hồn thiết bị chuyển hóa nước thay đổi với lưu lượng cấp nước điều kiện nhiệt độ chất rắn vào thiết bị chuyển hóa nước = 1100 K, độ chuyển hóa chất rắn thiết bị chuyển hóa nhiên liệu = 0,9 hàm lượng Fe2O3 = 10 wt.% 74 HVTH: Hoàng Trung Kiên Trang Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành +A_H2O_1(8)/T)*R*T; end if (T >= 1000 && T < 6000) H_H2 = (-A_H2_2(1)*T^(-2)+A_H2_2(2)*log(T)/T+A_H2_2(3) +A_H2_2(4)*T/2+A_H2_2(5)*T^2/3+A_H2_2(6)*T^3/4+A_H2_2(7)*T^4/5 +A_H2_2(8)/T)*R*T; H_O2 = (-A_O2_2(1)*T^(-2)+A_O2_2(2)*log(T)/T+A_O2_2(3) +A_O2_2(4)*T/2+A_O2_2(5)*T^2/3+A_O2_2(6)*T^3/4+A_O2_2(7)*T^4/5 +A_O2_2(8)/T)*R*T; H_CH4 = (-A_CH4_2(1)*T^(-2)+A_CH4_2(2)*log(T)/T+A_CH4_2(3) +A_CH4_2(4)*T/2+A_CH4_2(5)*T^2/3+A_CH4_2(6)*T^3/4+A_CH4_2(7)*T^4/5 +A_CH4_2(8)/T)*R*T; H_CO2 = (-A_CO2_2(1)*T^(-2)+A_CO2_2(2)*log(T)/T+A_CO2_2(3) +A_CO2_2(4)*T/2+A_CO2_2(5)*T^2/3+A_CO2_2(6)*T^3/4+A_CO2_2(7)*T^4/5 +A_CO2_2(8)/T)*R*T; H_H2O = (-A_H2O_2(1)*T^(-2)+A_H2O_2(2)*log(T)/T+A_H2O_2(3) +A_H2O_2(4)*T/2+A_H2O_2(5)*T^2/3+A_H2O_2(6)*T^3/4+A_H2O_2(7)*T^4/5 +A_H2O_2(8)/T)*R*T; end % % Heat of formation of each species: % (1) Heat of formation of Fe2O3 is calculated from the reaction: % 2Fe + (3/2)O2 = Fe2O3 H_Fe2O3_f = H_Fe2O3 - (2*H_Fe + 3/2*H_O2); % (2) Heat of formation of CH4 is calculated from the reaction: % C + 2H2 = CH4 H_CH4_f = H_CH4 - (H_C + 2*H_H2); % (3) Heat of formation of FeO is calculated from the reaction: % Fe + (1/2)O2 = FeO H_FeO_f = H_FeO - (H_Fe + 1/2*H_O2); % (4) Heat of formation of CO2 is calculated from the reaction: % C + O2 = CO2 H_CO2_f = H_CO2 - (H_C + H_O2); % (5) Heat of formation of H2O is calculated from the reaction: % H2 + (1/2)O2 = H2O H_H2O_f = H_H2O - (H_H2 + 1/2*H_O2); % - HVTH: Hoàng Trung Kiên Trang 109 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành % Test with HSC Chemistry 5.1 results: %dH_rx = H_Fe2O3_f %dH_rx = H_CH4_f %dH_rx = H_FeO_f %dH_rx = H_CO2_f %dH_rx = H_H2O_f % % Heat of the reaction (1): % 4Fe2O3 + CH4(g) = 8FeO + CO2(g) + 2H2O(g) (2) dH_rx = 8*H_FeO_f + H_CO2_f + 2*H_H2O_f - (4*H_Fe2O3_f + H_CH4_f); %dH_rx = 8*H_FeO + H_CO2 + 2*H_H2O - (4*H_Fe2O3 + H_CH4); end Chương trình dHrx_FeO_H2O.m: function dH_rx = dHrx_FeO_H2O(T) %This function is used to calculate the reaction heat of Wustite (FeO) %oxidation to Magnetite (Fe3O4) by Steam (H2O) by the following %reaction equation: % 3FeO + H2O(g) = Fe3O4 + H2(g) (1) %Data taken from: NASA Glenn Coefficients for Calculating Thermodynamic %Properties of Individual Species, NASA/TP-2002-211556 %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 9th 2016 % Parameters %Gas constant, unit: J/(mol-K) R = 8.31451; % %Enthalpy of formation of each species present in Reaction (1): %(1) Enthalpy of Fe3O4: %H_Fe3O4 = 0; if (T >= 298 && T < 800) H_Fe3O4 = (-A_Fe3O4_1(1)*T^(2)+A_Fe3O4_1(2)*log(T)/T+A_Fe3O4_1(3) +A_Fe3O4_1(4)*T/2+A_Fe3O4_1(5)*T^2/3+A_Fe3O4_1(6)*T^3/4 +A_Fe3O4_1(7)*T^4/5+A_Fe3O4_1(8)/T)*R*T; end if (T >= 800 && T < 850) H_Fe3O4 = (-A_Fe3O4_2(1)*T^(2)+A_Fe3O4_2(2)*log(T)/T+A_Fe3O4_2(3) +A_Fe3O4_2(4)*T/2+A_Fe3O4_2(5)*T^2/3+A_Fe3O4_2(6)*T^3/4 HVTH: Hoàng Trung Kiên Trang 110 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành +A_Fe3O4_2(7)*T^4/5+A_Fe3O4_2(8)/T)*R*T; end if (T >= 850)%Caution! Temperature boundary for Fe2O3 is 1800 K H_Fe3O4 = (-A_Fe3O4_3(1)*T^(2)+A_Fe3O4_3(2)*log(T)/T+A_Fe3O4_3(3) +A_Fe3O4_3(4)*T/2+A_Fe3O4_3(5)*T^2/3+A_Fe3O4_3(6)*T^3/4 +A_Fe3O4_3(7)*T^4/5+A_Fe3O4_3(8)/T)*R*T; end %(2) Enthalpy of Fe: H_Fe = 0; if (T >= 200 && T < 500) H_Fe = (-A_Fe_1(1)*T^(-2)+A_Fe_1(2)*log(T)/T+A_Fe_1(3) +A_Fe_1(4)*T/2+A_Fe_1(5)*T^2/3+A_Fe_1(6)*T^3/4+A_Fe_1(7)*T^4/5 +A_Fe_1(8)/T)*R*T; end if (T >= 500 && T < 800) H_Fe = (-A_Fe_2(1)*T^(-2)+A_Fe_2(2)*log(T)/T+A_Fe_2(3) +A_Fe_2(4)*T/2+A_Fe_2(5)*T^2/3+A_Fe_2(6)*T^3/4+A_Fe_2(7)*T^4/5 +A_Fe_2(8)/T)*R*T; end if (T >= 800 && T < 1042) H_Fe = (-A_Fe_3(1)*T^(-2)+A_Fe_3(2)*log(T)/T+A_Fe_3(3) +A_Fe_3(4)*T/2+A_Fe_3(5)*T^2/3+A_Fe_3(6)*T^3/4+A_Fe_3(7)*T^4/5 +A_Fe_3(8)/T)*R*T; end if (T >= 1042)%Caution! Temperature boundary for Fe is 1184K for Fe! H_Fe = (-A_Fe_4(1)*T^(-2)+A_Fe_4(2)*log(T)/T+A_Fe_4(3) +A_Fe_4(4)*T/2+A_Fe_4(5)*T^2/3+A_Fe_4(6)*T^3/4+A_Fe_4(7)*T^4/5 +A_Fe_4(8)/T)*R*T; end % (3) Enthalpy of FeO: H_FeO = (-A_FeO(1)*T^(-2)+A_FeO(2)*log(T)/T+A_FeO(3) +A_FeO(4)*T/2+A_FeO(5)*T^2/3+A_FeO(6)*T^3/4+A_FeO(7)*T^4/5 +A_FeO(8)/T)*R*T; % (4) Enthalpy of Graphite (C): H_C = 0; if (T >= 200 && T < 600) H_C = (-A_C_1(1)*T^(-2)+A_C_1(2)*log(T)/T+A_C_1(3) +A_C_1(4)*T/2+A_C_1(5)*T^2/3+A_C_1(6)*T^3/4+A_C_1(7)*T^4/5 +A_C_1(8)/T)*R*T; end if (T >= 600)%Caution! Temperature boundary for Graphite is 2000 K H_C = (-A_C_2(1)*T^(-2)+A_C_2(2)*log(T)/T+A_C_2(3) HVTH: Hoàng Trung Kiên Trang 111 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành +A_C_2(4)*T/2+A_C_2(5)*T^2/3+A_C_2(6)*T^3/4+A_C_2(7)*T^4/5 +A_C_2(8)/T)*R*T; end % (5) Enthalpies of H2, O2, CH4, CO2, and H2O: if (T >= 200 && T < 1000) H_H2 = (-A_H2_1(1)*T^(-2)+A_H2_1(2)*log(T)/T+A_H2_1(3) +A_H2_1(4)*T/2+A_H2_1(5)*T^2/3+A_H2_1(6)*T^3/4+A_H2_1(7)*T^4/5 +A_H2_1(8)/T)*R*T; H_O2 = (-A_O2_1(1)*T^(-2)+A_O2_1(2)*log(T)/T+A_O2_1(3) +A_O2_1(4)*T/2+A_O2_1(5)*T^2/3+A_O2_1(6)*T^3/4+A_O2_1(7)*T^4/5 +A_O2_1(8)/T)*R*T; H_CH4 = (-A_CH4_1(1)*T^(-2)+A_CH4_1(2)*log(T)/T+A_CH4_1(3) +A_CH4_1(4)*T/2+A_CH4_1(5)*T^2/3+A_CH4_1(6)*T^3/4+A_CH4_1(7)*T^4/5 +A_CH4_1(8)/T)*R*T; H_CO2 = (-A_CO2_1(1)*T^(-2)+A_CO2_1(2)*log(T)/T+A_CO2_1(3) +A_CO2_1(4)*T/2+A_CO2_1(5)*T^2/3+A_CO2_1(6)*T^3/4+A_CO2_1(7)*T^4/5 +A_CO2_1(8)/T)*R*T; H_H2O = (-A_H2O_1(1)*T^(-2)+A_H2O_1(2)*log(T)/T+A_H2O_1(3) +A_H2O_1(4)*T/2+A_H2O_1(5)*T^2/3+A_H2O_1(6)*T^3/4+A_H2O_1(7)*T^4/5 +A_H2O_1(8)/T)*R*T; end if (T >= 1000 && T < 6000) H_H2 = (-A_H2_2(1)*T^(-2)+A_H2_2(2)*log(T)/T+A_H2_2(3) +A_H2_2(4)*T/2+A_H2_2(5)*T^2/3+A_H2_2(6)*T^3/4+A_H2_2(7)*T^4/5 +A_H2_2(8)/T)*R*T; H_O2 = (-A_O2_2(1)*T^(-2)+A_O2_2(2)*log(T)/T+A_O2_2(3) +A_O2_2(4)*T/2+A_O2_2(5)*T^2/3+A_O2_2(6)*T^3/4+A_O2_2(7)*T^4/5 +A_O2_2(8)/T)*R*T; H_CH4 = (-A_CH4_2(1)*T^(-2)+A_CH4_2(2)*log(T)/T+A_CH4_2(3) +A_CH4_2(4)*T/2+A_CH4_2(5)*T^2/3+A_CH4_2(6)*T^3/4+A_CH4_2(7)*T^4/5 +A_CH4_2(8)/T)*R*T; H_CO2 = (-A_CO2_2(1)*T^(-2)+A_CO2_2(2)*log(T)/T+A_CO2_2(3) +A_CO2_2(4)*T/2+A_CO2_2(5)*T^2/3+A_CO2_2(6)*T^3/4+A_CO2_2(7)*T^4/5 +A_CO2_2(8)/T)*R*T; H_H2O = (-A_H2O_2(1)*T^(-2)+A_H2O_2(2)*log(T)/T+A_H2O_2(3) HVTH: Hoàng Trung Kiên Trang 112 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành +A_H2O_2(4)*T/2+A_H2O_2(5)*T^2/3+A_H2O_2(6)*T^3/4+A_H2O_2(7)*T^4/5 +A_H2O_2(8)/T)*R*T; end % % Heat of formation of each species: % (1) Heat of formation of Fe3O4 is calculated from the reaction: % 3Fe + 2O2 = Fe3O4 H_Fe3O4_f = H_Fe3O4 - (3*H_Fe + 2*H_O2); % (2) Heat of formation of CH4 is calculated from the reaction: % C + 2H2 = CH4 %H_CH4_f = H_CH4 - (H_C + 2*H_H2); % (3) Heat of formation of FeO is calculated from the reaction: % Fe + (1/2)O2 = FeO H_FeO_f = H_FeO - (H_Fe + 1/2*H_O2); % (4) Heat of formation of CO2 is calculated from the reaction: % C + O2 = CO2 %H_CO2_f = H_CO2 - (H_C + H_O2); % (5) Heat of formation of H2O is calculated from the reaction: % H2 + (1/2)O2 = H2O H_H2O_f = H_H2O - (H_H2 + 1/2*H_O2); % % Test with HSC Chemistry 5.1 results: %dH_rx = H_Fe2O3_f %dH_rx = H_CH4_f %dH_rx = H_FeO_f %dH_rx = H_CO2_f %dH_rx = H_H2O_f % % Heat of the reaction (1): % 3FeO + H2O = Fe3O4 + H2 dH_rx = H_Fe3O4_f + H_H2 - (3*H_FeO_f + H_H2O_f); %dH_rx = 8*H_FeO + H_CO2 + 2*H_H2O - (4*H_Fe2O3 + H_CH4); end PL 2.2.5 Các chương trình tính nhiệt dung riêng cấu tử Phần tính nhiệt dung riêng cấu tử bao gồm chương trình Nội dung chương trình sau: Chương trình Cp_CH4.m function y = Cp_CH4(T) HVTH: Hoàng Trung Kiên Trang 113 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành %This function is used to calculate the heat capacity of Methane gas %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 Parameters; %Gas constant, unit: J/(mol-K) R = 8.31451; % Calculation of heat capacity y = 0; if (T >= 200 && T < 1000) y = (A_CH4_1(1)*T^(-2)+A_CH4_1(2)/T+A_CH4_1(3) +A_CH4_1(4)*T+A_CH4_1(5)*T^2+A_CH4_1(6)*T^3+A_CH4_1(7)*T^4)*R; end if (T >= 1000 && T < 6000) y = (A_CH4_2(1)*T^(-2)+A_CH4_2(2)/T+A_CH4_2(3) +A_CH4_2(4)*T+A_CH4_2(5)*T^2+A_CH4_2(6)*T^3+A_CH4_2(7)*T^4)*R; End Chương trình Cp_CO.m function y = Cp_CO(T) %This function is used to calculate the heat capacity of CO gas %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 %Gas constant, unit: J/(mol-K) R = 8.31451; %The thermodynamic parameters of the species related to the reaction: A_CO_1 = [1.489045E+04 -2.922286E+02 5.724527E+00 -8.176235E-03 1.456903E-05 -1.087746E-08 3.027942E-12 -1.303132E+04 7.859241E+00]; A_CO_2 = [4.619197E+05 -1.944705E+03 5.916714E+00 -5.664283E-04 HVTH: Hoàng Trung Kiên Trang 114 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành 1.398815E-07 -1.787680E-11 9.620936E-16 -2.466261E+03 1.387413E+01]; %The values of A_CO2_1 are valid for the temperature range of 200 - 1000K, %these of A_CO_2 for the temperature range of 1000 - 6000K % Calculation of heat capacity y = 0; if (T >= 200 && T < 1000) y = (A_CO_1(1)*T^(-2)+A_CO_1(2)/T+A_CO_1(3) +A_CO_1(4)*T+A_CO_1(5)*T^2+A_CO_1(6)*T^3+A_CO_1(7)*T^4)*R; end if (T >= 1000 && T < 6000) y = (A_CO_2(1)*T^(-2)+A_CO_2(2)/T+A_CO_2(3) +A_CO_2(4)*T+A_CO_2(5)*T^2+A_CO_2(6)*T^3+A_CO_2(7)*T^4)*R; end Chương trình Cp_CO2.m function y = Cp_CO2(T) %This function is used to calculate the heat capacity of CO2 %Gas constant, unit: J/(mol-K) R = 8.31451; %The thermodynamic parameters of the species related to the reaction: A_CO2_1 = [4.943651E+04 -6.264116E+02 5.301725E+00 2.503814E-03 2.127309E-07 -7.689989E-10 2.849678E-13 -4.528198E+04 7.048279E+00]; A_CO2_2 = [1.176962E+05 -1.788791E+03 8.291523E+00 -9.223157E-05 4.863677E-09 -1.891053E-12 6.330037E-16 -3.908351E+04 2.652669E+01]; %The values of A_CO2_1 are valid for the temperature range of 200 - 1000K and %these of A_CO_2 for the temperature range of 1000 - 6000K % Calculation of heat capacity y = 0; if (T >= 200 && T < 1000) y = (A_CO2_1(1)*T^(2)+A_CO2_1(2)/T+A_CO2_1(3)+A_CO2_1(4)*T+A_CO2_1(5)*T^2+A_CO2_1(6)* T^3+A_CO2_1(7)*T^4)*R; end if (T >= 1000 && T < 6000) y = (A_CO2_2(1)*T^(2)+A_CO2_2(2)/T+A_CO2_2(3)+A_CO2_2(4)*T+A_CO2_2(5)*T^2+A_CO2_2(6)* T^3+A_CO2_2(7)*T^4)*R; end Chương trình Cp_Fe.m HVTH: Hồng Trung Kiên Trang 115 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành function y = Cp_Fe(T) %This function is used to calculate the heat capacity of Fe %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 Parameters; %Gas constant, unit: J/(mol-K) R = 8.31451; % Calculation of heat capacity y = 0; if (T >= 200 && T < 500) y = (A_Fe_1(1)*T^(-2)+A_Fe_1(2)/T+A_Fe_1(3) +A_Fe_1(4)*T+A_Fe_1(5)*T^2+A_Fe_1(6)*T^3+A_Fe_1(7)*T^4)*R; end if (T >= 500 && T < 800) y = (A_Fe_2(1)*T^(-2)+A_Fe_2(2)/T+A_Fe_2(3) +A_Fe_2(4)*T+A_Fe_2(5)*T^2+A_Fe_2(6)*T^3+A_Fe_2(7)*T^4)*R; end if (T >= 800 && T < 1042) y = (A_Fe_3(1)*T^(-2)+A_Fe_3(2)/T+A_Fe_3(3) +A_Fe_3(4)*T+A_Fe_3(5)*T^2+A_Fe_3(6)*T^3+A_Fe_3(7)*T^4)*R; end if (T >= 1042)%Caution! Temperature boundary for Fe is 1184K for Fe! y = (A_Fe_4(1)*T^(-2)+A_Fe_4(2)/T+A_Fe_4(3) +A_Fe_4(4)*T+A_Fe_4(5)*T^2+A_Fe_4(6)*T^3+A_Fe_4(7)*T^4)*R; End Chương trình Cp_Fe2O3.m function y = Cp_Fe2O3(T) %This function is used to calculate the heat capacity of Fe2O3 %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 %Gas constant, unit: J/(mol-K) R = 8.31451; HVTH: Hoàng Trung Kiên Trang 116 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành %The thermodynamic parameters of the species related to the reaction: A_Fe2O3_1 = [-2.713212E+06 2.323048E+04 -6.765702E+01 1.306052E-01 -6.970688E-05 0.000000E+00 0.000000E+00 -2.256089E+05 4.228073E+02]; A_Fe2O3_2 = [3.333493E+08 -1.859962E+06 3.885138E+03 -3.562599E+00 1.231056E-03 0.000000E+00 0.000000E+00 1.058059E+07 2.555306E+04]; A_Fe2O3_3 = [-2.844009E+07 1.253786E+05 -2.119080E+02 2.221482E-01 -1.206765E-04 3.480738E-08 -4.165298E-12 -8.483927E+05 1.433002E+03]; %The valid temperature range of A_Fe2O3_1 is 298 - 600K, %A_Fe2O3_2 600 - 960K, and A_Fe2O3_3 960 - 1800K y = 0; % Calculation of heat capacity if (T >= 298 && T < 600) y = (A_Fe2O3_1(1)*T.^(-2)+A_Fe2O3_1(2)/T+A_Fe2O3_1(3) +A_Fe2O3_1(4)*T+A_Fe2O3_1(5)*T.^2+A_Fe2O3_1(6)*T.^3 +A_Fe2O3_1(7)*T.^4)*R; end if (T >= 600 && T < 960) y = (A_Fe2O3_2(1)*T^(-2)+A_Fe2O3_2(2)/T+A_Fe2O3_2(3) +A_Fe2O3_2(4)*T+A_Fe2O3_2(5)*T^2+A_Fe2O3_2(6)*T^3 +A_Fe2O3_2(7)*T^4)*R; end if (T >= 960)%Caution! Temperature boundary for Fe2O3 is 1800 K y = (A_Fe2O3_3(1)*T^(-2)+A_Fe2O3_3(2)/T+A_Fe2O3_3(3) +A_Fe2O3_3(4)*T+A_Fe2O3_3(5)*T^2+A_Fe2O3_3(6)*T^3 +A_Fe2O3_3(7)*T^4)*R; end end Chương trình Cp_FeO.m function y = Cp_FeO(T) %This function is used to calculate the heat capacity of FeO %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 %Gas constant, unit: J/(mol-K) R = 8.31451; HVTH: Hoàng Trung Kiên Trang 117 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành %The thermodynamic parameters of the species related to the reaction: A_FeO = [-1.179194E+04 1.388393E+02 2.999842E+00 1.274527E-02 -1.883886E-05 1.274258E-08 -3.042206E-12 -3.417351E+04 1.284759E+01]; % Calculation of heat capacity Cp_FeO = (A_FeO(1)*T^(2)+A_FeO(2)/T+A_FeO(3)+A_FeO(4)*T+A_FeO(5)*T^2 +A_FeO(6)*T^3+A_FeO(7)*T^4)*R; y = Cp_FeO; end Chương trình Cp_H2.m function y = Cp_H2(T) %This function is used to calculate the heat capacity of Hydrogen gas %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 Parameters; %Gas constant, unit: J/(mol-K) R = 8.31451; % Calculation of heat capacity y = 0; if (T >= 200 && T < 1000) y = (A_H2_1(1)*T^(-2)+A_H2_1(2)/T+A_H2_1(3) +A_H2_1(4)*T+A_H2_1(5)*T^2+A_H2_1(6)*T^3+A_H2_1(7)*T^4)*R; end if (T >= 1000 && T < 6000) y = (A_H2_2(1)*T^(-2)+A_H2_2(2)/T+A_H2_2(3) +A_H2_2(4)*T+A_H2_2(5)*T^2+A_H2_2(6)*T^3+A_H2_2(7)*T^4)*R; end Chương trình Cp_H2O.m function y = Cp_H2O(T) %This function is used to calculate the heat capacity of Steam %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T HVTH: Hoàng Trung Kiên Trang 118 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành %Date: 2016 %Last changed: September 23rd 2016 Parameters; %Gas constant, unit: J/(mol-K) R = 8.31451; % Calculation of heat capacity y = 0; if (T >= 200 && T < 1000) y = (A_H2O_1(1)*T^(-2)+A_H2O_1(2)/T+A_H2O_1(3) +A_H2O_1(4)*T+A_H2O_1(5)*T^2+A_H2O_1(6)*T^3+A_H2O_1(7)*T^4)*R; end if (T >= 1000 && T < 6000) y = (A_H2O_2(1)*T^(-2)+A_H2O_2(2)/T+A_H2O_2(3) +A_H2O_2(4)*T+A_H2O_2(5)*T^2+A_H2O_2(6)*T^3+A_H2O_2(7)*T^4)*R; End Chương trình Cp_ZrO2.m function y = Cp_ZrO2(T) %This function is used to calculate the heat capacity of Methane gas %Ref.: NASA Glenn Coefficients for Calculating Thermodynamic Properties % of Individual Species (NASA/TP-2002-211556), Glenn Research Center, % Cleveland, Ohio %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 23rd 2016 Parameters; %Gas constant, unit: J/(mol-K) R = 8.31451; y = 0; % Calculation of heat capacity if (T >= 200 && T < 500) y = (A_ZrO2_1(1)*T^(-2)+A_ZrO2_1(2)/T+A_ZrO2_1(3) +A_ZrO2_1(4)*T+A_ZrO2_1(5)*T^2+A_ZrO2_1(6)*T^3+A_ZrO2_1(7)*T^4)*R; end if (T >= 500 && T < 1445) y = (A_ZrO2_2(1)*T^(-2)+A_ZrO2_2(2)/T+A_ZrO2_2(3) +A_ZrO2_2(4)*T+A_ZrO2_2(5)*T^2+A_ZrO2_2(6)*T^3+A_ZrO2_2(7)*T^4)*R; end if (T >= 1445)%Caution! Temperature boundary for ZrO2 is 2620 K y = (A_ZrO2_3(1)*T^(-2)+A_ZrO2_3(2)/T+A_ZrO2_3(3) HVTH: Hoàng Trung Kiên Trang 119 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành +A_ZrO2_3(4)*T+A_ZrO2_3(5)*T^2+A_ZrO2_3(6)*T^3+A_ZrO2_3(7)*T^4)*R; end end PL 2.2.6 Các chương trình tính tốn hệ số dẫn nhiệt cấu tử Phần tính hệ số dẫn nhiệt cấu tử bao gồm chương trình Nội dung chương trình sau: Chương trình k_CH4.m: function y = k_CH4(T) %This function returns the thermal conductivity of methane gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_CH4 = [-1.869E-03 8.727E-05 1.179E-07 -3.614E-11]; y = A_CH4(1)+A_CH4(2)*T+A_CH4(3)*T^2+A_CH4(4)*T^3;%Unit: W/m-K end Chương trình k_CO.m: function y = k_CO(T) %This function returns the thermal conductivity of CO gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_CO = [5.067E-04 9.125E-05 -3.524E-08 8.1991E-12]; y = A_CO(1)+A_CO(2)*T+A_CO(3)*T^2+A_CO(4)*T^3;%Unit: W/m-K end Chương trình k_CO2.m function y = k_CO2(T) %This function returns the thermal conductivity of CO2 gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: HVTH: Hoàng Trung Kiên Trang 120 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành A_CO2 = [-7.215E-03 8.015E-05 5.477E-09 -1.053E-11]; y = A_CO2(1)+A_CO2(2)*T+A_CO2(3)*T^2+A_CO2(4)*T^3;%Unit: W/m-K end Chương trình k_H2.m function y = k_H2(T) %This function returns the thermal conductivity of hydrogen gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_H2 = [8.099E-03 6.689E-04 -4.158E-07 1.562E-10]; y = A_H2(1)+A_H2(2)*T+A_H2(3)*T^2+A_H2(4)*T^3;%Unit: W/m-K end Chương trình k_H2O.m function y = k_H2O(T) %This function returns the thermal conductivity of steam %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_H2O = [7.341E-03 -1.013E-05 1.801E-07 -9.100E-11]; y = A_H2O(1)+A_H2O(2)*T+A_H2O(3)*T^2+A_H2O(4)*T^3;%unit: W/m-K end PL 2.2.7 Các chương trình tính tốn độ nhớt cấu tử Phần tính độ nhớt cấu tử bao gồm chương trình Nội dung chương trình sau: Chương trình mu_CH4.m: function y = mu_CH4(T) %This function returns the viscosity of methane gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % - HVTH: Hoàng Trung Kiên Trang 121 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành %Parameters: A_CH4 = [1.596E-03 3.439E-05 -8.140E-09]; y = (A_CH4(1)+A_CH4(2)*T+A_CH4(3)*T^2)*1E-3;%unit: kg/m-s %Note: % The original unit in KDB is centipose (cP) and 1cP = 0.001kg/m-s end Chương trình mu_CO.m: function y = mu_CO(T) %This function returns the viscosity of CO gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_CO = [3.228E-03 4.747E-05 -9.648E-09]; y = (A_CO(1)+A_CO(2)*T+A_CO(3)*T^2)*1E-3;%unit: kg/m-s %Note: % The original unit in KDB is centipose and 1cP = 0.001kg/m-s end Chương trình mu_CO2.m: function y = mu_CO2(T) %This function returns the viscosity of CO2 gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_CO2 = [2.545E-03 4.549E-05 -8.649E-09]; y = (A_CO2(1)+A_CO2(2)*T+A_CO2(3)*T^2)*1E-3;%unit: kg/m-s %Note: % The original unit in KDB is centipose and 1cP = 0.001kg/m-s end Chương trình mu_H2.m: function y = mu_H2(T) %This function returns the viscosity of hydrogen gas %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 HVTH: Hoàng Trung Kiên Trang 122 Luận văn tốt nghiệp GVHD : TS Nguyễn Đặng Bình Thành % %Parameters: A_H2 = [2.187E-03 2.220E-05 -3.751E-09]; y = (A_H2(1)+A_H2(2)*T+A_H2(3)*T^2)*1E-3;%unit: kg/m-s %Note: % The original unit in KDB is centipose and 1cP = 0.001kg/m-s end Chương trình mu_H2O.m: function y = mu_H2O(T) %This function returns the viscosity of steam %Data taken from Korea Thermodynamics Data Bank (KDB) %Developed by: Thanh D.B Nguyen and Kien H.T %Date: 2016 %Last changed: September 22nd 2016 % %Parameters: A_H2O = [-3.189E-03 4.145E-05 -8.272E-10]; y = (A_H2O(1)+A_H2O(2)*T+A_H2O(3)*T^2)*1E-3;%unit: kg/m-s %Note: % The original unit in KDB is centipose and 1cP = 0.001kg/m-s end HVTH: Hoàng Trung Kiên Trang 123 ... Độ chuyển hóa chất rắn 50 Xgf 51 Xgs Độ chuyển hóa mêtan thiết bị chuyển hóa nhiên liệu Độ chuyển hóa nước thiết bị chuyển hóa nước 52 Xsf Độ chuyển hóa chất rắn thiết bị chuyển hóa nhiên liệu. .. wt.% chất mang oxy nhiệt độ chất rắn vào thiết bị chuyển hóa nhiên liệu độ chuyển hóa thiết bị chuyển hóa nhiên liệu 0,9 71 Hình 4.6 Thay đổi lượng chất rắn tuần hoàn hàm nhiệt độ chất. .. ngưng tụ hoàn toàn lượng nước theo sản phẩm này, vật liệu tuần hoàn Fe FeO chuyển hóa thành Fe3O4 Fe3O4 hình thành thiết bị ôxy hóa tiếp tục chuyển hóa thiết bị hoàn nguyên vật liệu tuần hoàn tạo