Nghiên cứu cơ chế khí hóa sinh khối và ứng dụng thiết kế, chế tạo lò đốt khí hóa công suất nhỏ

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - Nguyễn Hữu Linh Nghiên cứu chế khí hóa nhiên liệu sinh khối ứng dụng thiết kế, chế tạo lò đốt khí hóa cơng suất nhỏ CHUN NGÀNH: Khoa học kỹ thuật tính tốn Người hướng dẫn: TS Lê Đức Dũng Hà Nội – Năm 2017 MỤC LỤC MỤC LỤC i LỜI MỞ ĐẦU iv LỜI CẢM ƠN vi LỜI CAM ĐOAN vii DANH MỤC BẢNG BIỂU viii DANH MỤC HÌNH VẼ .ix CHƯƠNG I GIỚI THIỆU 1.1 Khí hóa sinh khối 1.1.1 Khái niệm khí hóa sinh khối 1.1.2 Nhân tố khí hóa 1.1.3 Lò khí hóa sinh khối 1.2 Nhiên liệu sinh khối 1.2.1 Sinh khối gì? 1.2.2 Phân loại nhiên liệu sinh khối 1.2.3 Đặc điểm nhiên liệu sinh khối 1.2.4 Tiềm năng lượng sinh khối Việt Nam 1.3 Mơ hình hóa CFD q trình khí hóa sinh khối 1.4 Mục tiêu luận văn phương pháp nghiên cứu 12 1.4.1 Mục tiêu luận văn 12 1.4.2 Phương pháp nghiên cứu 12 CHƯƠNG II BẾP KHÍ HÓA SINH KHỐI 13 2.1 Lý thuyết khí hóa sinh khối 13 2.1.1 Lý thuyết chung q trình khí hóa sinh khối 13 2.1.2 Các phương trình phản ứng khí hóa 15 2.1.3 Động học q trình khí hóa 17 2.1.3.1 Giới thiệu .17 2.1.3.2 Động học phản ứng khí hóa char 19 i 2.1.4 Ảnh hưởng điều kiện vận hành đến thành phần sản phẩm khí 21 2.1.4.1 Ảnh hưởng nhiệt độ 21 2.1.4.2 Ảnh hưởng tác nhân khí hóa 22 2.2 Q trình cháy lò khí hóa 23 2.3 Bếp khí hóa sinh khối 24 2.3.1 Mơ tả bếp khí hóa sinh khối 24 2.3.2 Bếp khí hóa sinh khối điển hình Việt Nam 25 CHƯƠNG III CƠ SỞ MÔ HÌNH HĨA ỨNG DỤNG TRONG Q TRÌNH KHÍ HĨA SINH KHỐI 27 3.1 Các phương trình chủ đạo dòng chất lưu 27 3.2 Mơ hình rối 29 3.3 Mơ hình xạ DO 30 3.4 Mơ hình dòng phản ứng 31 3.5 Mô hình mơi trường xốp 33 CHƯƠNG IV MƠ HÌNH HĨA BẾP KHÍ HĨA SINH KHƠI 34 4.1 Mơ hình cân hóa học 34 4.2 Mơ hình CFD q trình khí hóa lớp nhiên liệu bếp khí hóa 37 4.2.1 Mơ tả chung 37 4.2.2 Xây dựng mô hình CFD 40 4.2.2.1 Nguyên lý làm việc mô hình CFD 40 4.2.2.2 Cơ chế truyền nhiệt bên lớp nhiên liệu 41 4.2.2.3 Mơ hình phản ứng lớp nhiên liệu 42 4.2.2.3 Điều kiện biên .43 4.2.2.4 Mô tả hàm macro UDF 44 CHƯƠNG V THIẾT KẾ BẾP ĐUN KHÍ HĨA TRẤU 45 5.1 Ảnh hưởng việc bố trí miệng gió thứ cấp đến q trình cháy khí 46 5.1.1 Mơ tả chung 46 5.1.2 Kết mô 49 5.1.2.1 Sự phân bố trường nhiệt độ 49 5.1.2.2 Sự phân bố khí sản phẩm cháy 52 ii 5.1.3 Nhận xét 54 5.2 Ảnh hưởng đường kính buồng đốt đến trình cháy vùng cháy thứ cấp 55 5.2.1 Mô tả chung 55 5.2.2 Kết mô 55 5.2.2.1 Sự phân bố nhiệt độ lửa 55 5.2.2.2 Sự phân bố nồng độ khí 56 5.2.3 Nhận xét 58 5.3 Thiết kế mẫu bếp đun khí hóa trấu 59 KẾT LUẬN 61 TÀI LIỆU THAM KHẢO 62 PHỤ LỤC 64 iii LỜI MỞ ĐẦU Vào năm cuối kỷ 20 đầu kỷ 21, đôi với phát triển mạnh mẽ kinh tế xã hội, nhu cầu sử dụng lượng tăng nhanh Hiện nay, giới nguồn cung cấp lượng chủ yếu từ nhiên liệu hóa thạch Tuy nhiên, trữ lượng nguồn nhiên liệu có hạn ngày giảm dần Hơn nữa, việc khai thác sử dụng dầu mỏ than đá thải khí CO2, SO2 NOx gây hiệu ứng nhà kính, làm nhiễm mơi trường môi sinh, làm thay đổi nghiêm trọng đến khí hậu tồn cầu Từ đặt cho nhân loại nhiệm vụ cấp bách phải tìm nguồn lượng có khả thay lượng từ nhiên liệu gốc khống, tái tạo thân thiện với môi trường như: lượng sinh khối, lượng hydro, địa nhiệt, sức gió lượng mặt trời, thủy triều Sinh khối chủ yếu tạo chuyển hóa lượng ánh sáng mặt trời vào trồng (thực vật) đường quang hợp Nhiên liệu sinh khối bao gồm: Gỗ, chất thải gỗ (mạt cưa, phôi bào), phân động vật, nông sản phế thải từ nông nghiệp rơm rạ, trấu, thân lõi ngô, vỏ lạc Ưu điểm nhiên liệu sinh học so với nguồn nhiên liệu dầu mỏ truyền thống có hàm lượng lưu huỳnh nitơ thấp, khơng gây hiệu ứng nhà kính có cân CO2 Việt Nam có tiềm lượng sinh khối lớn với nông nghiệp lâu đời ¾ diện tích lãnh thổ đồi núi Theo ước tính năm có khoảng 60 triệu sinh khối tạo từ phụ phẩm nơng nghiệp Tuy nhiên, có 30-40% sinh khối sử dụng cho mục đích lượng chủ yếu nhiên liệu để đun nấu hộ gia đình số lượng nhỏ sử dụng cho phát điện, chủ yếu nhà máy đường Còn lại loại sinh khối vỏ trấu dư thừa rơm xử lý cách đốt trực tiếp thành đống, kết lượng kèm theo phát thải chất ô nhiễm môi trường Bên cạnh đó, sử dụng sinh khối bếp khí hóa để nấu ăn vùng nơng thơn Việt Nam cho thấy tiềm lớn có nhiều lợi so với việc đốt cháy trực tiếp hiệu suất nhiệt cao, khói đen lượng khí thải, bụi Vì vậy, việc sử dụng sinh khối iv hiệu có ý nghĩa đặc biệt quan trọng đặc biệt nhu cầu lượng gia tăng Trên sở này, hướng dẫn TS Lê Đức Dũng, em lựa chọn đề tài: “Nghiên cứu chế khí hóa sinh khối ứng dụng thiết kế, chế tạo lò đốt khí hóa cơng suất nhỏ” Nội dung luận văn bao gồm chương sau: Chương I Trình bày khái niệm cơng nghệ khí hóa, sinh khối phương pháp nghiên cứu Chương II Trình bày lý thuyết khí hóa sinh khối q trình cháy bếp khí hóa Chương III Trình bày sở lý thuyết CFD ứng dụng khí hóa lớp nhiên liệu Chương IV Trình bày mơ hình khí hóa mơ hình cân hóa học mơ hình CFD Chương V Trình bày q trình cháy khí phần mềm mơ ANSYS Fluent thiết kế bếp khí hóa trấu v LỜI CẢM ƠN Với lòng biết ơn sâu sắc em xin gửi lời cảm ơn chân thành tới thầy giáo TS Lê Đức Dũng, người tận tình hướng dẫn em suốt trình thực luận văn Bản luận văn hồn thành, phần lớn nhờ vào ý kiến, định hướng, gợi ý dẫn dắt bảo thầy Em xin gửi lời cảm ơn đến thầy cô Viện nghiên cứu quốc tế Khoa học Kỹ thuật tính tốn ln giúp đỡ tạo điều kiện thuận lợi cho em hoàn thành luận văn Em xin gửi lời cảm ơn đến Ths Nguyễn Danh Nam cung cấp số liệu thí nghiệm nhận xét, góp ý anh cho luận văn Hà Nội, ngày 25 tháng năm 2017 Học viên Nguyễn Hữu Linh vi LỜI CAM ĐOAN Tôi xin cam đoan luận văn tơi tự tính tốn thiết kế nghiên cứu hướng dẫn thầy giáo TS Lê Đức Dũng Để hoàn thành luận văn sử dụng tài liệu ghi mục tài liệu tham khảo, ngồi khơng sử dụng tài liệu khác mà không ghi Nếu sai tơi xin chịu hình thức kỷ luật theo quy định Hà Nội, ngày 25 tháng năm 2017 Học viên Nguyễn Hữu Linh vii DANH MỤC BẢNG BIỂU Bảng 1.1 So sánh cơng nghệ khí hóa [4] Bảng 1.2 Tính chất số loại sinh khối [2] Bảng 1.3 Một số tính chất quan trọng ảnh hưởng sinh khối [2] Bảng 1.4 Tiềm sinh khối Việt Nam [8] Bảng 1.5 Một số ứng dụng CFD mơ hình hóa q trình khí hóa cháy sinh khối [10] 11 Bảng 2.1 Cơ chế phản ứng phản ứng khí hóa cháy [4] 15 Bảng 4.1 Thành phần công nghệ thành phần hóa học trấu [15] 35 Bảng 4.2 Giá trị phương trình khí hóa (1) 37 Bảng 4.3 Tỷ lệ thành phần khí hỗn hợp theo thể tích 37 Bảng 4.4 Thơng số động học trình 1, 2, [13] 43 Bảng 5.1 Điều kiện biên 47 viii DANH MỤC HÌNH VẼ Hình 1.1 Q trình khí hóa sinh khối Hình 1.2 Sơ đồ hình thành khí tổng hợp biogas với ứng dụng [5] Hình 2.1 Quá trình cháy lò khí hóa 24 Hình 2.2 Mơ hình bếp khí hóa 25 Hình 2.3 Cơ chế khí hóa bếp khí hóa TLUD 25 Hình 4.1 Q trình khí hóa lớp nhiên liệu [14] 39 Hình 4.2 Mơ hình hình học lớp nhiên liệu biên 41 Hình 4.3 Quy trình thực hàm UDF [13] 45 Hình 5.1 Mẫu bếp số 46 Hình 5.2 Mẫu bếp số 46 Hình 5.3 Miền tính tốn biên mơ hình 47 Hình 5.4 Sự phân bố nhiệt độ mẫu bếp theo mặt cắt dọc theo chiều cao bếp 49 Hình 5.5 Sự phân bố nhiệt độ mẫu bếp mặt cắt dọc theo chiều cao bếp 50 Hình 5.6 Sự phân bố nhiệt độ mặt cắt qua cụm miệng gió cấp hai 50 Hình 5.7 Sự phân bố nhiệt độ mặt cắt ngang vị trí sát đáy nồi 51 Hình 5.8 Sự thay đổi nhiệt độ đường streamline khơng khí cấp hai 51 Hình 5.9 Sự thay đổi nhiệt độ đường streamline nhiên liệu 52 Hình 5.10 Sự phân bố khí CH4 theo tỷ lệ mol 53 Hình 5.11 Sự phân bố khí CO theo tỷ lệ mol 53 Hình 5.12 Sự phân bố khí H2 theo tỷ lệ mol 53 Hình 5.13 Sự phân bố khí CO2 theo tỷ lệ mol 54 Hình 5.14 Sự phân bố khí H2O theo tỷ lệ mol 54 Hình 5.15 Kiểm tra độ hội tụ trường hợp 55 Hình 5.16 Sự phân bố nhiệt độ lửa mặt cắt qua tâm miệng gió thứ cấp 55 Hình 5.17 Sự phân bố nhiệt độ theo mặt cắt dọc thân bếp 56 ix PHỤ LỤC Phụ lục: Hàm UDF #include "udf.h" int i; /* Solid properties*/ /* Particle density [kg/m3] */ real solid_rho=400.0; /* Initial porosity */ real porosity_int=0.2; /* Rice husk composition */ real frac_moist=0.09; real frac_vol=0.562; real frac_c=0.126; real frac_ash=0.222; /* Heat of reaction C to CO [J/kg] */ real h_C_CO = 9210830.0; /* Heat of evaporation of water [J/kg] */ real h_evap = 2256900.0; /* Heat of reaction fro devolatilization [J/kg] */ real h_devol = 411100.0; /* Species enthanies */ real T_ref = 298.15; real cp_co = 1043.0; real cp_o2 = 919.31; real cp_h2o = 2014.0; real cp_vol = 1522.0; /* Volatile properties */ 64 real frac_vol_h2 = 0.0272; real frac_vol_co = 0.2123; real frac_vol_ch4 = 0.1974; real frac_vol_co2 = 0.5632; /* Devolatilization rate parameters [J/kmol] */ real A1 = 312000.0; real E1 = 7.4e+07; /* Universal gas constant */ real R = 8314.0; /* Char combustion rate parameters [J/kmol] */ real A2 = 1.45e+11; real E2 = 1.83e+08; /* Evaporation Arrhenius constant */ real A3 = 5.13e+10; real E3 = 8.8e+07; /*Defination of User Defined Memories: (Utilized for description of solid properties of cells) UDM0: Initial total mass UDM1: Total mass UDM2: C mass UDM3: Volatile mass UDM4: Moisture mass UDM5: Porosity UDM6: Total mass source UDM7: Moisture mass source UDM8: Volatile mass source UDM9: Ash mass 65 UDM10: CO mass source, from C combustion UDM11: CO mass source, from valatiles UDM12: C mass fixed char reacted */ enum { /* UDM0 */ m_init_tot, /* UDM1 */ m_tot, /* UDM2 */ m_C, /* UDM3 */ m_vol, /* UDM4 */ m_moist, /* UDM5 */ m_por, /* UDM6 */ source_total, /* UDM7 */ source_moist, /* UDM8 */ source_vol, /* UDM9 */ m_ash, /* UDM10 */ source_COC, /* UDM11 */ source_COvol, /* UDM12 */ source_C, /* UDM13 */ moist_flux, /* UDM14 */ vol_flux, /* UDM15 */ Hydrogen, }; /*******************************/ /* Initialization of start parameters */ /*********************************/ DEFINE_ON_DEMAND(reset_solid) 66 { /* declare domain pointer since it is not passed as an argument to the DEFINE marcro */ Domain *d; Thread *t; cell_t c; /* Get the domain using Fluent utility */ d = Get_Domain(1); /* Loop over all cell threads in the domain */ thread_loop_c (t,d) { if (FLUID_THREAD_P(t) && THREAD_VAR(t).fluid.porous) { /* loop over all cells */ begin_c_loop_all (c,t) { /* Reset of all UDFs */ C_UDMI(c,t,0) = 0.0; C_UDMI(c,t,1) = 0.0; C_UDMI(c,t,2) = 0.0; C_UDMI(c,t,3) = 0.0; C_UDMI(c,t,4) = 0.0; C_UDMI(c,t,5) = 0.0; C_UDMI(c,t,6) = 0.0; C_UDMI(c,t,7) = 0.0; C_UDMI(c,t,8) = 0.0; C_UDMI(c,t,9) = 0.0; 67 C_UDMI(c,t,10) = 0.0; C_UDMI(c,t,11) = 0.0; C_UDMI(c,t,12) = 0.0; C_UDMI(c,t,13) = 0.0; C_UDMI(c,t,14) = 0.0; C_UDMI(c,t,15) = 0.0; C_T(c,t) = 300; C_UDMI(c,t,20) = 0.0; C_UDMI(c,t,23) = 0.0; C_UDMI(c,t,24) = 0.0; /* Initialization of porosity */ C_UDMI(c,t,m_por) = porosity_int; /* Initial total solid mass in cell */ C_UDMI(c,t,m_init_tot) = (1.0 C_UDMI(c,t,m_por))*solid_rho*C_VOLUME(c,t); /* Time dependent total solid mass in cell */ C_UDMI(c,t,m_tot) = (1.0 C_UDMI(c,t,m_por))*solid_rho*C_VOLUME(c,t); /* Time dependent C mass in cell */ C_UDMI(c,t,m_C) = frac_c * C_UDMI(c,t,m_tot); /* Time dependent valatile mass in cell */ C_UDMI(c,t,m_vol) = frac_vol * C_UDMI(c,t,m_tot); /* Time dependent moisture mass in cell */ C_UDMI(c,t,m_moist) = frac_moist * C_UDMI(c,t,m_tot); /* Ash mass in cell */ 68 C_UDMI(c,t,m_moist) = frac_ash * C_UDMI(c,t,m_tot); } end_c_loop_all (c,t) } } } /****************************************/ /* defination of moisture source term */ /*****************************************/ DEFINE_SOURCE(source_h2o, c, t, dS, eqn) { real s_h2o; real dm_h2o; real k_moist; /* Evaporation rate of water content [s^-1] */ k_moist = A3 * exp(-E3/(R*C_T(c,t))); /* max evaporation during timestep [kg] */ dm_h2o = C_UDMI(c,t,m_moist)*k_moist*CURRENT_TIMESTEP; /* if amount of moisture for evaporation exceeds cell content of moisture */ if (dm_h2o >= C_UDMI(c,t,m_moist)) { /* Evaporation of all cell moisture content */ C_UDMI(c,t,source_moist) = C_UDMI(c,t,m_moist)/(CURRENT_TIMESTEP*C_VOLUME(c,t)); 69 } else{ /* Arrhenius determined moisture release */ C_UDMI(c,t,source_moist) = dm_h2o/(CURRENT_TIMESTEP*C_VOLUME(c,t)); } s_h2o = C_UDMI(c,t,source_moist); return s_h2o; } /****************************************************************/ /* Main of determination of valatiles realease and H2 source term */ /****************************************************************/ DEFINE_SOURCE(source_h2, c, t, dS, eqn) { real s_h2; real k_vol; real dm_vol; k_vol = A1*exp(-E1/(R*C_T(c,t))); dm_vol = C_UDMI(c,t,m_vol)*k_vol*CURRENT_TIMESTEP; /* if release of volatiles exceeds valatiles content */ if (dm_vol >= C_UDMI(c,t,m_vol)) { /* Realease based on present valatile content */ C_UDMI(c,t,source_vol) = C_UDMI(c,t,m_vol)*k_vol / (CURRENT_TIMESTEP*C_VOLUME(c,t)); 70 } else { /* Realease based on devolatilization rate */ C_UDMI(c,t,source_vol) = C_UDMI(c,t,m_vol)*k_vol / C_VOLUME(c,t); } s_h2 = C_UDMI(c,t,source_vol)*frac_vol_h2; C_UDMI(c,t,Hydrogen)=s_h2; return s_h2; } /* Defination of CH4 source term */ DEFINE_SOURCE(source_ch4, c, t, dS, eqn) { real s_ch4; s_ch4 = C_UDMI(c,t,source_vol)*frac_vol_ch4; return s_ch4; } /* Defination of CO2 source term */ DEFINE_SOURCE(source_co2, c, t, dS, eqn) { real s_co2; s_co2 = C_UDMI(c,t,source_vol)*frac_vol_co2; return s_co2; 71 } /* Main caculation of C combustion and defination of oxy source term */ DEFINE_SOURCE(source_o2, c, t, dS, eqn) { real s_o2; real m_o2; real dm_o2; real dm_c; real k_char; real f_o2 = 0.0; int n; Thread *tf; face_t f; Thread *t1; cell_t c1; /* Define index for O2 */ i=1; /* Mass rate of O2 based on cell content */ m_o2 = C_VOLUME(c,t)*C_UDMI(c,t,m_por)*C_R(c,t)*C_YI(c,t,i); /* Combustion rate of char [s^-1] */ k_char = A2*exp(-E2/(R*C_T(c,t))); 72 /*kinetic determined max char combustion during timeslep [kg] */ dm_c = C_UDMI(c,t,m_C)*k_char*CURRENT_TIMESTEP; /* Net to O2 mass flux into cell */ c_face_loop(c,t,n) { f = C_FACE(c,t,n); tf = C_FACE_THREAD(c,t,n); /* At faces with out flow */ if(F_FLUX(f,tf) > 0.0) { f_o2 += F_FLUX(f,tf) * C_YI(c,t,i); } else { /* At faces with flow */ c1 = F_C1(f,tf); t1 = F_C1_THREAD(f,tf); f_o2 += -F_FLUX(f,tf) * C_YI(c1,t1,i); } } /* Available amount of O2 */ dm_o2 = m_o2 + f_o2 * CURRENT_TIMESTEP; /* if kinetic char combustion is higher than stoichiometric amount of available O2 */ /* Limiting factor for combustion is O2 diffusion/convection */ if (dm_c*1.0/12.0 >= dm_o2 * 1.0/16.0) 73 { /* Source of CO from combustion */ C_UDMI(c,t,source_COC) = (dm_o2 * 28.0/16.0)/(C_VOLUME(c,t)*CURRENT_TIMESTEP); /* reacted amount of fixed C */ C_UDMI(c,t,source_C) = (dm_o2 * 12.0/16.0)/(C_VOLUME(c,t)*CURRENT_TIMESTEP); /* Source of O2 from combustion */ s_o2 = -dm_o2/(C_VOLUME(c,t)*CURRENT_TIMESTEP); } else /* Else kinetic of char combustion is lower than stoichiometric rate of available O2 */ /* Limiting factor for char combustion is char kinetic rate */ { /* Source of CO from combustion */ C_UDMI(c,t,source_COC) = C_UDMI(c,t,m_C)*28.0/12.0*k_char/C_VOLUME(c,t); /* Reacted amount of fixed C */ C_UDMI(c,t,source_C) = C_UDMI(c,t,m_C)*k_char/C_VOLUME(c,t); /* Source of O2 from combustion */ s_o2 = - C_UDMI(c,t,source_COC) * 16.0/28.0; } return s_o2; } /* Definition of CO source term */ 74 DEFINE_SOURCE(source_co, c, t, dS, eqn) { real s_co; C_UDMI(c,t,source_COvol) = C_UDMI(c,t,source_vol)*frac_vol_co; s_co = C_UDMI(c,t,source_COC) + C_UDMI(c,t,source_COvol); return s_co; } /* Defination of resulting mass source term */ DEFINE_SOURCE(source_mass, c, t, dS, eqn) { real s_m; /* total mass source */ C_UDMI(c,t,source_total) = C_UDMI(c,t,source_C) + C_UDMI(c,t,source_vol) + C_UDMI(c,t,source_moist); s_m = C_UDMI(c,t,source_total); return s_m; } /* Defination of energy source term */ DEFINE_SOURCE(source_energy, c, t, dS, eqn) { real s_energy; real h_react; real h_species; real rad; 75 real energy; real area = 0.01; real radtemp; h_react = C_UDMI(c,t,source_C)*h_C_CO C_UDMI(c,t,source_moist)*h_evap + C_UDMI(c,t,source_vol)*h_devol; h_species = C_UDMI(c,t,source_COC)*(cp_co - (16.0/28.0)*cp_o2) + C_UDMI(c,t,source_vol)*cp_vol + C_UDMI(c,t,source_moist)*cp_h2o*(C_T(c,t)T_ref); s_energy = h_react + h_species; /* radtemp=C_UDMI(c,t,24); rad=(epsilon * sigma * area * pow(radtemp,4.0)/0.000245; energy = -s_energy-rad; */ return s_energy; } /* Defination of porosity */ DEFINE_PROFILE(porosity, t, i) { cell_t c; /* loops over cells in a cell thread */ begin_c_loop(c,t) { F_PROFILE(c,t,i) = C_UDMI(c,t,m_por); } end_c_loop(c,t) } 76 /* Update of user defined memories*/ DEFINE_EXECUTE_AT_END(execute_at_end) { Domain *d; Thread *t; cell_t c; int count = 0; real ricehusk_radtemp; d = Get_Domain(1); thread_loop_c (t,d) { if (FLUID_THREAD_P(t) && THREAD_VAR(t).fluid.porous) { begin_c_loop(c,t) { /* Update of C cell mass */ C_UDMI(c,t,m_C) = C_UDMI(c,t,m_C) C_UDMI(c,t,source_C) * (CURRENT_TIMESTEP*C_VOLUME(c,t)); if (C_UDMI(c,t,m_C) < 0.0) { C_UDMI(c,t,m_C) = 0.0; } /* Update of volatile cell mass */ C_UDMI(c,t,m_vol) = C_UDMI(c,t,m_vol) C_UDMI(c,t,source_vol) * (CURRENT_TIMESTEP*C_VOLUME(c,t)); if (C_UDMI(c,t,m_vol) < 0.0) { C_UDMI(c,t,m_vol) = 0.0; 77 } C_UDMI(c,t,23) = C_UDMI(c,t,23) + 1; C_UDMI(c,t,24) = (C_UDMI(c,t,23)+C_T(c,t))/C_UDMI(c,t,23); ricehusk_radtemp = C_UDMI(c,t,24); /* Update of moisture cell mass */ C_UDMI(c,t,m_moist) = C_UDMI(c,t,m_moist) C_UDMI(c,t,source_moist) * (CURRENT_TIMESTEP*C_VOLUME(c,t)); if (C_UDMI(c,t,m_moist) < 0.0) { C_UDMI(c,t,m_moist) = 0.0; } /* Update of total solid cell mass */ C_UDMI(c,t,m_tot) = C_UDMI(c,t,m_C) + C_UDMI(c,t,m_vol) + C_UDMI(c,t,m_moist) + C_UDMI(c,t,m_ash); } /* Update of porosity */ C_UDMI(c,t,m_por) = 1.0 C_UDMI(c,t,m_tot)/(solid_rho*C_VOLUME(c,t)); end_c_loop (c,t) } } } 78 ... văn: Nghiên cứu chế khí hóa sinh khối thiết kế, chế tạo lò đốt khí hóa cơng suất nhỏ sau:  Nghiên cứu lý thuyết cơng nghệ khí hóa sinh khối  Ứng dụng phương pháp số CFD để mô tả trình khí hóa. .. tài: Nghiên cứu chế khí hóa sinh khối ứng dụng thiết kế, chế tạo lò đốt khí hóa cơng suất nhỏ Nội dung luận văn bao gồm chương sau: Chương I Trình bày khái niệm cơng nghệ khí hóa, sinh khối. .. trình khí hóa sinh khối Hình 1.2 Sơ đồ hình thành khí tổng hợp biogas với ứng dụng [5] Hình 2.1 Quá trình cháy lò khí hóa 24 Hình 2.2 Mơ hình bếp khí hóa 25 Hình 2.3 Cơ chế khí hóa bếp khí hóa